open TELEMAC-MASCARET The mathematically superior suite of solvers

Bonjour, débutant de Telemac2d, j'essaye de réaliser un cas simple, malheureusement, j'ai quelques soucis.
Je souhaite modéliser un bassin de fond plat relié à une rivière(canal) en pente. Mon bassin étant relié à la mer, je place une condition SL imposée en aval et un débit imposé (de la rivière) en amont. J'ai créé un fichier de frontière liquide pour imposer une SL variable qui correspond à la marée. J'aimerais savoir comment imposer comme condition initiale le niveau moyen d'eau 4.13m pour le bassin. Et si quelqu'un pouvait m'expliquer comment TELEMAC2d distingue les endroits "où il y a de l'eau", car ma rivière possède de l'eau mais elle est au dessus du niveau de la mer (niveau moyen:4.13m).
Je n'arrive pas à joindre mon fichier frontière liquide, si vous avez un fichier .lqd exemple à me fournir, je vous en remercierai grandement.
Cordialement,
Julien BAILLS

Hello,

It is not a simple case, because the initial conditions are not trivial. You will have to give them either in Blue Kenue of Fudaa, or to implement them in Fortran in subroutine condin.f.

You can start by prescribing INITIAL CONDITIONS as PRESCRIBED ELEVATION, giving the value of 4.13 m. This will do for the basin. Then you can, in condin.f taken in a FORTRAN FILE specific to your case, try to add a depth in the river, based on the coordinates, if it is possible.

Another more radical solution is to put an average constant depth in the river and let it flow until things settle.

With best regards,

Jean-Michel Hervouet

Hello,I thank you for your advice. I'll try the "radical solution" (the easiest one, I think).

Faithfully,

Julien Baills

Bonjour, j'ai un peu de mal avec les type BIEF_MESH et BIEF_OBJ. Je souhaite écrire dans le fichier utimp_telemac2d.f, j'aimerais sélectionner le point 14 par exemple du maillage: que faut il faire?
do i=1,NPOIN
if(MESH%I(i).EQ.14) then
... ...
end ...
ou
if(MESH%KNOLG%I(i).EQ.14) then


sauf que KNOLG est utilisé pour le parallélisme et que dans mon cas, je n'utilise pas le code parallèle. Est ce que la première formulation fonctionnerait?

Cordialement,
Julien BAILLS

Hello,

If you are not using parallelism, basically it is just:

PRINT*,'H(14)=',H%R(14)

for printing the depth at point 14 on the listing.

Regards,

JMH