open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Telemac users,

I have issue passing a spatial information to the setup of various grain sizes:
*I have mapped various habitat in my geometry file and the name of the variable is "habitat". Each habitat got a number (9,10 or 11).
*The later geometry file is called in both Telemac2d and Gaia steering files. Both steering files contain the following lines:
NUMBER OF PRIVATE ARRAYS = 1
NAMES OF PRIVATE VARIABLES = 'habitat'
*I amended the fortran subroutine init_compo.f with the following declaration at the beginning of the fortran file:
USE BIEF
USE DECLARATIONS_TELEMAC
USE DECLARATIONS_GAIA
USE DECLARATIONS_SPECIAL
IMPLICIT NONE
*within the fortran subroutine (which is called in both Telemac and GAIA steering files), I tried to associate a mesh node number with its habitat using the following:
PRIVE%ADR(1)%P%R(J) with J the node number. I also tried to output this variable using the WRITE and PRINT statements, but failed, and I do not know why.
*The fortran code compiles well, and got the following line in outputs:
"VARIABLE habitat FOUND IN THE GEOMETRY FILE"

But the different grain sizes are not implemented in the geometry and no sediment in present in the initial file.
Would you have any idea of what I may be missing?
Do I need to call the fortran file in both Telemac and Gaia steering files?

please find attached the geometry file and the fortran code.
Thanks for any idea!

Hello Claire,

USER_INIT_COMPO is a SISYPHE subroutine, not a GAIA one. You should use the equivalent subroutine from GAIA so that it may be taken into account I think.

As it is a question related to sediment transport, I moved your post to the dedicated module.

Hope this helps,

Chi-Tuan

Hello Claire,

To complete Chi-Tuân’s answer, user_bed_init is the subroutine in Gaia dedicated to advanced initialization of grain sizes, bed layers etc.

Best regards

Sara