open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I'm working on a hydro-sedimentary model. I modified the subroutine init_compo.f for the sediment distribution because I have two classes (1:d50=40mm downstream ; 2:d50=1mm upstream), I used the geographic coordinate y to define my areas but I would like to be able to define them on bluekenue to make it more practical. I tried to create new variables and read them in the init_compo subroutine but an error message appears:

INITIALISATION EN CAS DE COUPLAGE : PART= 0
start reading sediment fractions
from file: geo_frac_iniBIS.slf
variable name:FRACLAY1_CL1
ERROR 5002 DURING CALL OF GET_DATA_TIME_SRF:READ
ERROR TEXTE: UNKONWN ERROR

In the geo file I have my two variables : FRACLAY1_CL1 and FRACLAY2_CL2.

Do you know how I can solve this error ?
Thanks,
Claire

Change the variable names to
FRACLAY 1 CL 1
So without the underscore. The variable names might have undrscores, but the model looks for variable names without them.(Even with the double space between the first number and the secon CL, I think).

Thanks now it works !

Hello,

I have an other problem... It works when my first variable = 1 and the second =0, but when I put my values (0.3/0.7) it doesn't work. For each modified nodes I have the error :
PROBLEM IN LAYER: J,TEST1 28041 0.30000001192092896
IN LAYER 1 SUM OF FRACTIONS NOT 1
PROBLEM IN LAYER: J,TEST2 28041 0.30000001192092896
IN LAYER 2 SUM OF FRACTIONS IS NOT 1
For this node, I put my first varaible to 0.7 and the second to 0.3

Do you know what is the problem ?
Thanks,
Claire

The values you have given have a round off error. They are defined in single precision, I suspect, but are read as a double precision creating a small error, whichmeans the sum of fractions is no longer 1.00000000000 in double precision.

In the file layer.f an error margin is defined, which by default is very low). If you add that file to your Fortran file and set zero to a larger value (above 1e-8 I think), it should disappear. For example
DATA ZERO/1.D-7/