open TELEMAC-MASCARET The mathematically superior suite of solvers

My problem about multiple layers has been solved. I am inspired from the exampel give by Pablo, and also from the use manual about sisyphe. I also thanks the useful discuss with Cyamin.

In here I want to share my experience in here.


Firstly, we need to define the depth of each layers and each types of sediment composition. (Note: if we defined the active layer in sis.cas file. the depth of ES(J,1) should be eactly the same with the active layer)


For the multiple layers (6 layers, 4 types of sediments):

First, we need to define the depth of the each layers (we only need to define 5 layers): the depth are 0.1; 0.2; 0.5; 2m (also has one more layer on the bottom)

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


Firstly define the class of sediment types: for example, a multiple layers (6) have several types of sediments: 0.2; 0.8; 2; 6mm. that means the j=4. Has 4 types of sediment.

For the first layers, the percentage of the 4 types of sediment are 60%; 20; 10%; 10%

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


Which means the sediment percentage for 0.2; 0.8; 2; 6mm are 60%,20%, 20%,0%

….


For 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


Thanks again for the useful suggestion