open TELEMAC-MASCARET The mathematically superior suite of solvers

Thank you for the suggestion. Now the detailed output file shows that the problems comes from a conflict in Init_sediment.f file. But I cannot figure it out where the conflict is.
here are my new init_compo.f and output files.
regards and thanks
Najib

In init_avai you will see that the error message is printed when there si the line write(lu,1815)...
There is 2 possibilities, allways with the same check:
IF (NLAYER%I(J).GT.NOMBLAY) THEN
Try to see which condition is true and leads to this message

Hi Najib,
sorry but you have some mess in your subroutine?
You declare NPOIN as local variable with INTEGER NPOIN. And then you do a loop with NPOIN.
So first you should set NPOIN = MESH%NPOIN.

(Christophe, please correct me if I'm wrong.)

Best regards,
Clemens

I already did that,the suggestion came from Christophe and solved the Segmentation fault.
Now the problem is that the NLAYER becomes bigger than NOMBLAY.
look at line 502,531-534 and 659-640 of my fortran file. I'm searching to find why this NLAYER becomes greater than NOMBLAY.
najib

Please post your cas file together with the fortran file for every computation you make.

here are the files.
Najib

Can you make these changes in the cas file and try again?

Yes but you let the line do J=1, NPOIN
This is this point clemens noticed

NPOIN doesn't have any value in your subroutine!

Thank you, i set the NPOIN to MESH%NPOIN and it seems to run without problem.
thank you everybody for your help. I will post the files here later on.
Najib

your are confuss about the layer and the fraction of the sediment.

There are some disharmony in your sis.cas and the xx.f files.

In your sis.cas file, you set 5 types of sediment 0.0367;0.0193;0.0105;0.007;0.00475;

but in the fortran file, you set 9 types of the sediment.

The right set is like following. For example if you set the 6 layers, the depth of the layer is 0.1m, 0.2m, 0.5m, 2m, 3m. laster lay is equal to total depth minus the 5 layer.



NCOUCHES(J) = 6
ES(J,1) = 0.1D0
ES(J,2) = 0.2D0
ES(J,3) = 0.5D0
ES(J,4) = 2.D0
ES(J,5) = 3.D0

For the sediment types, in your case ,you have 5 types of sediment, you could set the first layer

AVAI(J,1,1) = 0.6 D0
AVAI(J,1,2) = 0.2 D0
AVAI(J,1,3) = 0.2 D0
AVAI(J,1,4) = 0.0 D0
AVAI(J,1,5) = 0.0 D0



***

***



The 5th layer

AVAIL(J,5,1) = 0.0D0
AVAIL(J,5,2) = 0.2D0
AVAIL(J,5,3) = 0.1D0
AVAIL(J,5,4) = 0.6D0
AVAIL(J,5,5) = 0.0D0



In your case, you need rewrite the codes in INIT_COMPO.f

Hope it is helpful for your case.