open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, everyone
I want to give a space-varying CRITICAL EROSION SHEAR STRESS OF THE MUD PER LAYER. For example, if(Y(I).ge.3720000.D0) then
TOCE_VASE(1) = 0.4
else if (Y(I).gt.3660000.D0.and.Y(I).lt.3720000.D0) then
TOCE_VASE(1) = 0.6
else if (Y(I).gt.3560000.D0.and.Y(I).lt.3660000.D0) then
TOCE_VASE(1) = 0.6
endif

But I notice that the definition of CRITICAL EROSION SHEAR STRESS OF THE MUD PER LAYER aims at a layer instead of a node. Is it right? And which subroutine should I use?

Thank you very much.

Dear Zqzuoan,
thanks for your post. Effectively, the variable CRITICAL EROSION SHEAR STRESS is designed to be defined by layers. You can probably manipulate their values, for example in the subroutine suspension_erosion_coh.f

I hope it helps,

Pablo

Dear Pablo,
Thank you very much. I have read the subroutine suspension_erosion_coh.f, so I cannot give different critical erosion stress in different points?