open TELEMAC-MASCARET The mathematically superior suite of solvers

hi everybody,

Is it possible to do a spatial distribution of sediment classes. I mean setting a specific class in a certain area and an other class in an other area.

thx and bests


Josias

Hi Josias,

the sediment classes must be specified globally. So if you have one class in area 1 and another class in area 2 you need to specify 2 sediment classes. But then you set the initial mass of class 2 in area 1 and class 1 in area 2 to zero. This can be done in user_bed_init.f setting the RATIO_INIT(iclass,ilayer,inode) to zero or 1 (in case of 2 classes) depending of the area.

example of fortran code:
(see also in the examples hippodrome-t2d)
!loop over all nodes
DO IPOIN = 1,NPOIN
! area 1 has x-coordinate lower than 100
IF (mesh%X%R(ipoin).LT.100.) THEN
! the initial ratio of class 1 is set to one and of class 2 is set to zero
RATIO_INIT(1,1,IPOIN) = 1.0
RATIO_INIT(2,1,IPOIN) = 0.0
ENDIF
END DO

Best regards,
Rebekka

hi everybody, hi rebekha.

just for some precision, is the ESTRATUM 1 the layer above the ESTRATUM 2.

thx

bests

Josias

Hi Josias,
yes in gaia the layers are numbered in negative z-direction. Layer 1 has contact to the water body, below is layer 2 etc.
Beware that ESTRATUM is the initial thickness.
from the user manual: An active layer will be added inside this stratification at the beginning of the computation. Therefore, the actual number of layers will be equal to the number of layers for the initial stratification plus one.
The thicknesses of the layers during the computation are stored in the variable ES. This variable is written to the gaia results file.
Best regards,
Rebekka

hi everybody, hi Rebekha,

Rebekha, thanks for your reply, well detailed , appreciated.

I need someone to check this user_bed_init I modified, it doesn't work properly for me.
also can I use the 'AND' keyword for more precision when defining different areas (and How to use it ) , I mean:
IF(MESH%X%R(IPOIN).GE.6.D0 AND LT.7D0)

thx for your help

josias

hi,
sorry, I add, the gaia steering file,
maybe the problem is not from the fortran (or from the fortran only).

thx
josias

Hi
The conditions are not well organised in you fortran file.
Your code, only for the Y part is like this:The last condition is never reached as Y is necessarily in the 2 first cases...

On the other point, you could use AND in your condition. Just have a look on some fortran tutorial...
The syntax is like: IF ((Y.GT.10).AND.(X.LT.8)) THEN ...

Hope this helps

Hi everybody,

Thx c.coulet for your reply.

Based on this previous reply and on some fortran tutos, here the modifications I did.
Is it right?

thx

Josias

Hi
Hard to say as we don't know exactly what you want to have at the end but sounds better...