open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all,

When I running the 2 layers sediment problems, I found it is OK.

But, I try to define a 5 layers in the sedimeent transport. I use the 'SUBROUTINE INIT_COMPO' to define it

Set the 'NCOUCHES(J)=5', and then set

ES (J,1)= 0.01.D0
...
ES (J,5)= 0.4.D0

After that, I set the sediment diameter on each layer.


MEAN DIAMETER OF THE SEDIMENT =0.72D-3;2.49D-3;3.57D-3;4.62D-3;5.59D-3

And also fraction on each layers.


INITIAL FRACTION FOR PARTICULAR SIZE CLASS = 0.567;0.119;0.175;0.067;0.072

But when I runing it, it show a problem in it.



I think may be some thing wrong with my set, or I misunderstand some parameters

Hello,

Just to be clear, you want 5 sediment classes AND 5 sediment layers, or just 5 sediment classes?

Here are some suggestions for your init_compo:
Regards,
Costas

Thanks cyamin and Pablo, I am read the materails you give to me.


I want 5 sediment classes AND 5 sediment layers, For the 5 layers, I think we should set NCOUCHES(J) = 5,

Yes, I missed it.

Dear Harry,
your screen capture shows serious (basic) fortran programming issues, e.g syntax problems 0.01.D0, undefined variables, etc.

Please take a look to the attached file for inspiration.

I hope it helps,

Pablo

Dear Pablo,

I changed the fortran file in the right way, for 5 layers depth defined.


NCOUCHES(J) = 6

ES(J,1) = 0.01D0

...

ES(J,5) = 0.4D0

I also define the fraction of each layers


AVAIL(J,1,1) = AVA0(1)
AVAIL(J,1,2) = AVA0(2)
AVAIL(J,1,3) = AVA0(3)
...

AVAIL(J,5,1) = 0.2D0
AVAIL(J,5,2) = 0.5D0
AVAIL(J,5,3) = 0.3D0

For the value of AVA0(1), AVA0(2), AVA0(3). I defined through the sis.cas files

INITIAL FRACTION FOR PARTICULAR SIZE CLASS = 0.6; 0.3; 0.1


For the 5 sediment types on 5 different layers are defined

MEAN DIAMETER OF THE SEDIMENT =0.72D-3;2.49D-3;3.57D-3;4.62D-3;5.59D-3

In here I don't know how to define the sediment grading curves on each layers, so I hide the 'SEDIMENT DIAMETERS '

There still an unknow errors.





I think I made errors in the parameters 'MEAN DIAMETER OF THE SEDIMENT ' , Maybe this is not for the mean diameter for different layer.

For me, I think the parameters 'SEDIMENT DIAMETERS ' is only the sediment compostion for 1 layer.

Hello,

You defined NCOUCHES(J) = 6. Shouldn't it be NCOUCHES(J) = 5?

Costas

I think for 5 defined layers the NCOUCHES(J) should be = 6. Because there always have a solid layer on the bottom.

I already check it

The solid bottom is defined separately at the noerod subroutine. Init_compo deals with the transportable layers above the sold bed. I strongly believe that NCOUCHES(J) should not include the solid bed, but I am not an expert.

I am also not sure about whether it is include the solid bed. but I know the NCOUCHES(J) must one more than the ES(J,5) = 0.4D0.

Thant is what I am sure.