# Spatial sediment distribution in SISPHYE **URL:** \ **Category:** GAIA / SISYPHE\ **Created:** [August 22, 2023, 3:34pm UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357 "2023-08-22T15:34:24Z")\ **Posts on this page:** 5\ **Page:** 1
**Author:** ![user\_17461](https://opentelemac.org/forum/letter_avatar_proxy/v4/letter/u/ecb155/32.png) [@user\_17461](https://opentelemac.org/forum/u/user_17461)\ **Post date:** [August 22, 2023, 3:34pm UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357/1 "2023-08-22T15:34:24Z")
Hi [name removed], 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 [name removed] 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? [fortran.f](https://opentelemac.org/forum/uploads/short-url/h2NUO0fhND1fPpMjFq5TTTQ09XH.f) (2.09 KB) [NorthWales\_sis.cas](https://opentelemac.org/forum/uploads/short-url/woFi6NJ1eDDq1ioYSJthT47iRRe.cas) (3.29 KB) ---
**Author:** ![user\_17461](https://opentelemac.org/forum/letter_avatar_proxy/v4/letter/u/ecb155/32.png) [@user\_17461](https://opentelemac.org/forum/u/user_17461)\ **Post date:** [August 23, 2023, 9:27am UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357/2 "2023-08-23T09:27:16Z")
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. ---
**Author:** ![user\_17461](https://opentelemac.org/forum/letter_avatar_proxy/v4/letter/u/ecb155/32.png) [@user\_17461](https://opentelemac.org/forum/u/user_17461)\ **Post date:** [August 25, 2023, 1:09pm UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357/3 "2023-08-25T13:09:59Z")
Hi [name removed] someone please delete this topic ---
**Author:** ![user\_17461](https://opentelemac.org/forum/letter_avatar_proxy/v4/letter/u/ecb155/32.png) [@user\_17461](https://opentelemac.org/forum/u/user_17461)\ **Post date:** [November 27, 2023, 1:13pm UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357/4 "2023-11-27T13:13:17Z")
Hi, If anyone else is wanting to implement a spatially variable sediment distribution in SISYPHE, you can find my modified subroutine init\_compo here: [https://github.com/eroome/T2D-SIS-SED-DIST/tree/main](https://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, [name removed]. 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) ---
**Author:** ![user\_16725](https://opentelemac.org/forum/letter_avatar_proxy/v4/letter/u/c5a1d2/32.png) [@user\_16725](https://opentelemac.org/forum/u/user_16725)\ **Post date:** [December 23, 2024, 6:45pm UTC](https://opentelemac.org/forum/t/spatial-sediment-distribution-in-sisphye/7357/5 "2024-12-23T18:45:25Z")
Hi! thanks for sharing, I’m trying to spatialized grain size diameters but in GAIA, subroutine init\_compo it’s not in GAIA, I just want to spatialized D50 all over domain, but it seems that the correct way it’s by fractions in each node, so is that the only way to do it? by fractions?. On the other hand, in GAIA the subroutine USER\_BED\_INIT it’s the one who I’ve seen that people use, so I think that’s the correct subroutine to modify, is it correct? I appreciate your help, thanks.