open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,
In Telemac2D, spatially varying diffusivity coef. can be easily implemented in the Corvis.f subroutine. But in Telemac3D, this has to be done in the subroutine Viscos.f and there is no much explanation how to do it. I would be greateful if some one has already done this before and can share me his/her modified subroutine.

Thanking you in advance.

Yohannis

Hello

You can use the same loops as in 2D cases in order to use a horizontal space variation.
To make variation in vertical you need to know that nodes have the same xy coordinates and that vertical variation can be done in the same way than in the subroutine condim for instance.

I hope that this helps
with my best regards

Riadh ATA

Hi Yohannis,

you could also use a user defined variable (structure is already present in the code)
so you could read the horizontally varying velocity diffusivity from a map in your geometry file. This is done by the subroutine FONSTR.f in BIEF library. You could then introduce that variable in VISCOS.f and assign the values.

I have done this for salinity as a tracer in an estuary to introduce an initial salinity field. I then wrote the values to the vertical layers in the subroutine CONDIM.f like Riadh suggested. You can find a good scripting example there. I use blue kenue to interpolate measured values to my mesh and then add this map to my geometry file. I find this very easy to manipulate the data like this. I do the same with the BOTTOM FRICTION, but that is already completely inside the code.

Hope you can make it work!

Best regards,
Sven

Thanks Sven for the idea.