open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I want tos ask if there is a way to define initial cinditions like the initial bed composition and the non erodable areas via BlueKenue or Fudaa and to add these files to a selfain object like the bottom friction for Telemac.
I think this would be a nice and user friendly feature.

Best regards,
Clemens

Hello,

Sorry this is not yet done (to my knowledge), this would be to ask first to Fudaa or BlueKenue developers (Martin ?). It would just be a matter of agreeing on a variable name, and then in Sisyphe we can look for the geometry file if this variable is present.

With best regards,

Jean-Michel Hervouet

Hi Clemens,

You can add arbitrary variables to the selafin file with BlueKenue. As JMH said, it's a matter then of agreeing on an appropriate variable name and units.

Martin

Hi Clemens,

did you worked this out yet? I had the same idea since I use Blue Kenue for defining my areas of different bottom friction. now I would like to do the same with the D50 values of the sediment and the non erodable layers.

Jean-Michel,

can you tell my where T2D is looking at the geometry file to see if bottom friction values are present? I would like to use this part of the code as an example for the sisyphe part.

kind regards,

Sven

Hi Sven

In T2D, you could have a look into telemac2d.f
There is a vector called ALIRE. When a value of this vector is equal to 1, telemac try to read this value in the previous computational file or in the geo.
For the friction coefficient, you could search ALIRE(19) which is modified specificaly in the subroutine to could read the friction value.

Hope this helps

Hello,

The bottom and friction coefficients are read in subroutine fonstr.f (library BIEF), which is called by telemac2d.f, you can inspire from this.

With best regards,

Jean-Michel

Christophe and Jean-Michel,

as I understood, in sisyphe one can also use the function 'computation continued'. So the D50 values of the sediment are probably read from the previous computation file.
Would I be right to add a part of code there, which also checks the geo file for this information? (like T2D does)

Thanks for the very quick response!

Sven

Hello,

Yes we can imagine a subroutine called like FONSTR in Telemac-2D, that would be called by sisyphe.f, at the level where D50 is built, and with the same hierarchy of decisions.

JMH

Hello Sven,

in scalar mode i managed this via reading the initial grain size distribution from a text file. For definig the non erodable zone I used also text files in combination with the very nice inpoly function (in bief).

Hello Christophe,
last week I learnt how to address boundary nodes in parallel mode and it works fine.
But how to address generic node numbers since my procedures e.g. with the inpoly function don't work in parallel mode..

Best regards,
Clemens

Hi Clemens
You could follow the same way with the generic node.
KNOLG exist for all nodes, not only the boundary.

Good luck