open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to investigate the propagation of one sediment class imposed as initial condition value (INITIAL SUSPENDED SEDIMENTS CONCENTRATION VALUES). therefore I would like to avoid erosion of bed material. However, I like to allow imposed sediment class to deposit and erode. I followed suggestions here, using noerod.f to change CALL OV('X=Y+C ', X=ZR, Y=ZF, C=-100.D0, DIM1=NPOIN) to CALL OV('X=Y+C ', X=ZR, Y=ZF, C=0.D0, DIM1=NPOIN). I am using v8p1r2 and Gaia folder do not have noerod.f.
I tried this method anyway, but this didn't change anything. I still see additional material being generated in the domain. i know this from sortie file
FINAL BALANCE OF COH SEDIMENT1 (UNIT: g/l * M3)

INITIAL QUANTITY : 5870628.
FINAL QUANTITY : 0.4779305E+10
QUANTITY ENTERED THROUGH LIQ. BND.: 0.6819894E+10 ( IF <0 EXIT )
TOTAL QUANTITY LOST : 0.2046459E+10
RELATIVE ERROR : 0.3000720
Please advise on the methods to assign bed as non erodible.
Thanks in advance

Hello,

You should take a look at USER_BED_INIT.f and change the variable ESTRATUM.

Best regards,
P.N

Hello,

to complete the previous answer, to set a non erodible bed you can also use in your steering file the keyword: LAYERS INITIAL THICKNESS = 0. (without modifying user_bed_init.f).

Noerod.f can be used only in Sisyphe.

Finally I suggest you to switch to the last version of TELEMAC (v8p3).

Best regards,

Sara

Thanks for the answer. I was using LAYERS INITIAL THICKNESS = 0.0 already.
My major concern is how is such a large quantity of sediment entering from the liquid boundary? when I have no prescribed tracer value provided (QUANTITY ENTERED THROUGH LIQ. BND.: 0.6819894E+10 ( IF <0 EXIT )!!). Also to make it more clear, I extracted time series of COH sediment at 3 suspected nodes on liquid boundary (off-shore). The excess flux from the liquid boundary appears to start after 11-12 days of the model run and then continues for the rest of the 3-month simulation. I wonder why ? Does anyone have some suggestions?

Hello,

Can you please provide more information about your case like boundary conditions and your steering file?

Best regards

Sara

Hi Sara,

I think I managed to resolve it. I didn't use the PRESCRIBED TRACERS VALUES in the steering file (assuming that will imposes zero flux of sediment at boundary). Now with 0.0 value assigned for that keyword, the mass balance of sediment is making sens. Sorry for the confusion.