open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi all,
I have been attempting to set a simple sediment distribution (grain size fractions in one layer) using various methods. I have tried using Private Variables, reading a .DAT file directly and reading from Selafin file directly.
After reading all the relevant posts on this forum, I settled on using Private Variables, the modified subroutine init_compo has successfully ran but the result file does not show the specified grain size distribution. Is there a general consensus on the best approach to such a problem? and are there any obvious errors in my FORTRAN code?

Update, the Private Variable approach works! Now I just need to define 3 Private Variables and set them equal to the concentration of 3 sediment types in layer 1.

Hi can someone please delete this topic

Hi,
If anyone else is wanting to implement a spatially variable sediment distribution in SISYPHE, you can find my modified subroutine init_compo here: github.com/eroome/T2D-SIS-SED-DIST/tree/main

My code allows for three size classes to be defined in the geometry file, each class should be named FRAC1, FRAC2, FRAC3. The values assigned to each class represents the proportion of sediment in each class at a given node (the sum of all classes should be 1).
The D50 of each class is defined as usual in the SISYPHE steering file...
SEDIMENT DIAMETERS = xxx; xxx; xxx (m)