open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I am trying to run a simple 'multi-class' (1 gravel and 1 silt).

I get an error message during the first computation time step:

LAYER 1 CLASS 2 POINT 21551
AVAIL= -1.084895468465458E-009

[...]

This error message is set from line 341 to line 350 in the subroutine 'layer.f'.

Is there any reason why this truncation error is not smooth by the mean of:

AVAIL(J,1,I)=MAX(AVAIL(J,1,I),0.D0)
AVAIL(J,1,I)=MIN(AVAIL(J,1,I),1.D0)

??

I thank you in advance for your help with this issue.

With best regards.

Alexis

Hello,

we had similar problems, see the forum topic:
"multiple grain sizes + non erodable areas"

www.opentelemac.org/index.php?option=com...41&Itemid=62&lang=en

If you are working with v6p0 try to add the subroutine layer.f from the v6p1 source files to your princi file or switch entirely to v6p1.

Hope this helps,

Clemens

Hello,

layer was very difficult to write and we kept the tests to be sure that everything is OK. However I guess that you are in version 6.0. In the version 6.1 an error has been corrected (thanks to the checking) and the value of ZERO has been set to 1.D-7 because such cases as yours were actually not an error but just truncation error. You can either shift to version 6.1 (recommended) or change ZERO into 1.D-7 instead of 1.D-10.

With best regards,

Jean-Michel Hervouet