open TELEMAC-MASCARET The mathematically superior suite of solvers

I am try a case which have temporal-spatial varying wind and air pressure data to be input into Telemac. According to the manual, I have to modified the Meteo.f which is not a standard subroutine of Telemac. I wonder if there is some other way which is part of standard of telemac? Could any one tell me the format of the wind and air pressure file which could be accepted by standard telemac?

Hello,

I am not sure you have correctly understood one of my previous answer:
www.opentelemac.org/index.php/kunena/21-...he-calculation#42538

What is not standard for TELEMAC is the IDWM_T2D subroutine which is not in the sources folder. METEO subroutine is a standard subroutine available in standard TELEMAC but as written in the user manual, if you want to have meteo variable varying in time and space, you have to modify it.

What is available with TELEMAC for wind and air pressure are the following (as written in the user manual for the first 3 ones) without touching Fortran files:
- no wind (OPTION FOR WIND = 0),
- wind constant in time and space (OPTION FOR WIND = 1) provided by keywords,
- wind constant in space and variable in time (OPTION FOR WIND = 2) + see ASCII ATMOSPHERIC DATA FILE,
- wind and/or air pressure and/or air temperature varying in space and time and stored in the BINARY ATMOSPHERIC DATA FILE (see wind_txy example with steering file t2d_wind_txy_bin.cas and the BINARY ATMOSPHERIC DATA FILE = atm_wind_txy.slf). It is a SERAFIN file with variables with name WIND VELOCITY U, WIND VELOCITY V, SURFACE PRESSURE, AIR TEMPERATURE (that can be changed in METEO_SET_VAR_NAMES e.g. in FORTRAN FILE), these variables are stored at every node of the mesh and a few times. Between 2 times, a linear interpolation is done.

Except if using the last solution, as soon as you want to have wind variable in time and space, you have to implement something and adapt your input data to this implementation (or vice-versa).

Hope this helps,

Chi-Tuan

Hello,



In the case where I want to vary wind and pressure in space and time, I can't find how to create the file (binary or ASCII) allowing this (I managed to create and read a file of temporal variation only of wind and pressure, but as you mention above there is no effect of pressure due to spatial homogeneity).
In the example you gave (wind_txy from telemac2d) I can't figure out how to do this (there's only wind in atm_wind_txy.slf). There is however an ASCII file (fo1_wind_txy, attached file) mentioning the spatial and temporal variability of wind over 5 stations. But how do you add pressure data? There are no column names... (the first column corresponds to time, the second column corresponds to wind speed station 1, the third column corresponds to wind direction station 1...).

Thanks for your help

Raphaël

Hi
This is not a standard implementation so you have to adapt some of the fortran file to have this feature in your case.
The solution you're looking at is a particular case where an interpolation between values at specific station had to be done to have a spatial and temporal variation. Usually, the spatial and temporal variation of atmospheric data comes in a grid form and there is python scripts in Telemac to deal with those values and transform them in a binary atmospheric data file

Hello,
Thank you for your prompt reply.
Ok, I'm a beginner in this field so I'm a long way from mastering Fortran, but I'll look into it.
But can you confirm that there isn't a "simpler" method for seeing the effect of pressure variations on water levels (coastal zone of Charente Maritime)? I saw an article (Modification of telemac2D for storm surge use) that used the keyword INVERTED BAROMETER EFFECT (a method I'd like to use), but telemac3d doesn't recognize this keyword...
Raphaël

Sounds to be a specific version as Telemac2d doesn't have this keyword by default...
Could you share the link?
We recently added (for the v8.5) an option to include the effect of atmospheric pressure on the tidal boundary (by added the barometric effect on the computed tidal value) but the pressure is supposed to be known on the model

here the article: henry.baw.de/server/api/core/bitstreams/...c4e709fd8482/content

Hello Raphaël,

Indeed unfortunately what is described in the paper you refer has not been pushed to the official release of TELEMAC-2D, contrary to what was planned in the paper.

But as Christophe has written, since the latest tag v8p5r0 recently released las Friday, the inverted barometer method has been available (see manuals for keyword ATMOSPHERIC PRESSURE AT TIDAL BOUNDARIES) but as he has written, you have to provide a pressure field for the keyword BINARY ATMOSPHERIC DATA FILE. It is available both for TELEMAC-2D and TELEMAC-3D.

Hope this helps,

Chi-Tuan

Hi Raphaël,

I have been working on a similar problem (incorporating both wind and pressure into the ASCII file).

Have you made any progress on this since the last post?

I have a setup for Wind speed and direction in the ASCII file, maybe we could share some ideas...

Ed

Hello

I've kind of given up on problems with pressure since I've had decent results without it. If you ever manage to include it (and see its effects) in the v8p4 version, I'd love to hear from you.

Raphaël