open TELEMAC-MASCARET The mathematically superior suite of solvers

Bonjour à tous,

J'aurais besoin d'aide et d'idées !

Je modélise une onde de crue qui inonde une plaine. J'ai besoin de réaliser une cartographie de ma plaine inondée avec des isochrones (lignes) qui représentent le temps que met l'onde à atteindre chaque maille du modèle (une ligne pour les secteur en eau au bout de 10min, une autre pour les secteurs en eau au bout de 30min ... etc...). J'ai donc besoin de connaitre "la date" de première mise nen eau de chaque de mes mailles.
Et là je sèche complètement

Y a t'il un paramètre à activer pour connaitre la date/pas de temps de première mise en eau de mes mailles ?

NB : j'utilise un modèle qui tourne sous V7P1R1

Sinon peut-on extraire l'évolution de la hauteur d'eau au cours du temps en chaque point/nœud du modèle, au format ASCII par exemple, pour pouvoir ensuite créer un script d'extraction de la "date" de mise en eau ? (quand la valeur n'est plus =0)

Merci d'avance pour votre aide

For English speakers, at the end find a translation of what Laurie is trying to do and a proposition of solution.

French version :

Bonjour Laurie,
Je vous recommande de faire une sortie graphique (si c'est possible au niveau du stockage) toutes les 10 minutes par exemple(suivant le pas de temps que tu as besoin
pour tes tracés).
Ensuite, il y a des scripts python inclus avec la suite Telemac qui permettent d'extraire les hauteurs d'eau du fichier de résultat selafin.
Vous pouvez retrouver plus de détails dans les notebooks sur la récupération des données (par exemple pour représenter les hauteurs d'eau en 2D) et les adapter à l'utilisation.
Je vous joins une solution rapide, pour récupérer les coordonnées des nœuds et les hauteurs d'eau associées pour les sorties graphiques.
Reste à adapter le script python pour votre utilisation et trouver la date de mise en eau sur chaque nœud pour l'extraire dans un fichier txt ou directement faire le post-traitement voulu.


English version :
To sum up briefly, it has been asked, whether it was possible to know at which time in the simulation a node is getting wet for floods applications.

I would recommend to have the water depth as graphical output every 10 minutes for instance (choose the time step needed for your post-treatment) in the steering file.
Then, there are already notebooks included with the telemac software to extract water depth from selafin files.
All the details can be found in the jupyter notebooks on how to retrieve these date (for instance the plotting of wather depth in 2D) and adapt it to your use.
For now, I join you a fast naive solution to retrieve the coordinates of the nodes with their respective water depth for each graphical output.
Then, the script just need adaptation to get the times where the water depth is greater than zero and retrieve the dates in a txt file or directly plot some hraphs.

Good luck with these hints,
Jean-Paul

Hi Jean-Paul,

First of all, sorry for the French text, I completely forgot to translate it.
Secondly, thank you for your answer. I will try to make "my" solution based on your explanations.

If I find something else, I will share it here for sure.

Kind regards,
Laurie

My solution is based on the latest version of telemac. The notebooks and the scripts may not be present on your version, so you may need to download them from git or upload your version if you want to use this solution.

Hi,

You can also use the postTelemac plugin for QGIS.
Go to "outils", "extract max" and it allows to compute water arrival for a given value of water depth.

Alexis

Hi Alexis,

Thanks for your answer.
I'll try !
I meet some difficulties using PostTelemac in recent version of QGis. Do you ?


Kind regards,

On my laptop, Post telemac is running well with qgis 3.22.

Alexis

Hello Laurie,

Please have a look at the Malpasset example and the section in the TELEMAC-2D validation manual. You have a user Fortran subroutine (user_utimp_telemac2d) that may be interesting for you.

Hope this helps,

Chi-Tuan