open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

we are modelling reservoir flushing using Sisyphe coupled internally with Telemac3D.
The simulation is running well by using te default active layer thickness but
when we switch to CONSTANT ACTIVE LAYER THICKNESS = NO, then we get the error message "ERROR IN LAYER CASE 1" until approx. 3600 time steps (0.5s).
If we change the conditions to a constant active layer thickness of 0.1m, we get the same error message after approx. 13h of simulation.
Because of the calculation routine of the scin friction correction we would like to use a smaller activ elayer thickness and not the default value to take the changes in the scin friction dur to sediment transport in the flushing process into account.

Has anybody some hints for us?
Gabriele

Hello,

You can take subroutine layer.f (library Sisyphe) in your Fortran file and change parameter ZERO, which is currently at 1.D-7. We just check that a fraction is between 0 and 1 with some tolerance. Try ZERO=1.D-4 and see what happens. If a fraction happens to be neatly out of the range (0,1) we would need to have a look at the case.

With best regards,

Jean-Michel Hervouet