open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone

I have a quick question regarding mesh transformation - value 3 (fixed levels at a given altitude)

I am using this option to set a few layers at specific depths for my models and it works fine using CONDIM subroutine . I am utilising 5 layers, and they are set as follows:

DO IPLAN = 1,NPLAN
TRANSF_PLANE%I(IPLAN)=3
ENDDO
ZPLANE%R(2)=-30.D0
ZPLANE%R(3)=-20.D0
ZPLANE%R(4)=-10.D0

with layer 1 is at the seabed and layer 5 corresponds to the surface

my questions:

1. I am focusing on 2 regions for my work. The first zone has a depth of 80m, while the second area with 25m

2. Using the code above, the vertical layers would work for the 1st zone, since the depth covers the 1st and 2nd (prescribed) layer.

3. The thing that I am not sure is what happens to the bottom vertical layer (layer 1 and 2) for area 2 (shallow region) with depths less than 25 m, since my 2nd layer were set at 30m depth (which is less than the overall depth of zone 2)?

will the vertical layers at zone 2 still be 5 layers, or will it be less?

Thanks

Hello,

You will always have 5 layers but they may be superimposed. If you look at subroutine calcot.f, you will find at the beginning 2 hardcoded parameters that give the minimum distance between planes, and on the bottom (I think it is DISBOT, not sure and cannot check) it can be set to 0. If you want to do this you will have to set parameters as if you had tidal flats, because some elements may have no volumes, like on tidal flats.

With best regards,

Jean-Michel Hervouet

Thank you JMH,

I think that make sense. I keep on getting error after a few timestep using the fixed elevation.

I will have a look on the calcot subroutine as suggested.

regards
Anas