open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi everyone,

I would like to put in my t3d model with tpxo a wind

first, a constant wind N 300 and 3 m/s of velocity, I don't know how to manage U and V variables of the wind parameter in T3D...

secondly a varying wind according the meteo and according my simulation date(s). the main simulation date are august 2002 and august 2020




I use v8p2r1 for this simulation but I already installed v8p3 and v8p3r1 if you want I switch just tell me, what modification do I need to do with my files

if you need any other information or file just tell me...

thanks for helping

Josias

Hello Josias,

You have the answer in $HOMETEL/sources/utils/bief folder/meteo_telemac.f

DTR = pi/180, WINDD for the wind direction and WINDS for the wind magnitude.

You can have an example of ASCII ATMOSPHERIC DATA FILE in
$HOMETEL/examples/waqtel/heat_exchange example (meteo_data.txt file).

For example you can create such a file even if it is a constant wind with short names WINDS and WINDD (please read the WAQTEL user manual, in particular the Practical aspects chapter).

You can also check if it is right by adding variables X,Y for keyword VARIABLES FOR 2D GRAPHIC PRINTOUTS in your 2D RESULT FILE (1 time step is enough).

Hope this helps,

Chi-Tuan

hi everybody, hi chi-tuan,

thanks for your reply,

I got this error message at the very end of the computation.



after introducing this ascci file


with this modification of my steering file ("ATMOSPHERIC CONDITIONS")




thanks for your help

hi sorry for adding this,

as I am not a good handler of fortran files this is the modification I did in the meteo_telemac.f

DO IPOIN = 1,WINDX%DIM1
WINDX%R(IPOIN) = - WINDS%R(3)*SIN( WINDD%R(225)*pi/180 )
WINDY%R(IPOIN) = - WINDS%R(3)*COS( WINDD%R(225)*pi/180 )
ENDDO

is it right?
I ask also because I got this error message,

WINDX%R(IPOIN) = - WINDS%R(3)*SIN( WINDD%R(225)*pi/180 )
1
Error: Symbol ‘pi’ at (1) has no IMPLICIT type

thx
Josias

Hello Josias,

My previous post was only to indicate how to compute the 2 keywords WIND VELOCITY ALONG X and WIND VELOCITY ALONG Y when you have magnitude and direction of velocity:
WIND VELOCITY ALONG X = -magnitude * sin(direction)
WIND VELOCITY ALONG Y = -magnitude * cos(direction)
If you do not want to use these 2 keywords, you can manage wind with the
ASCII ATMOSPHERIC DATA FILE but you do not need to modify METEO_TELEMAC module (if you follow rules to write this file).

The error message tells you that variable PI is not defined, I detailed the formula to compute DTR but in METEO_TELEMAC it is already correctly computed.

If you are not familiar with Fortran, I would strongly encourage you to learn bases of this language which is not so difficult, in particular if you want to modify TELEMAC subroutines a little bit.

Hope this helps,

Chi-Tuan

Hi Chi-Tuan

Thank you for your reply , advices and precisions.

Bests

Josias

Hi everybody, hi Chi tuan,

thx for your reply,
I succesfully run the simulation with both ascii file first then the implementationn of X and Y.

the ascii was easier for me to catch but the X and Y took a while till I created this file

for the conversion.

bests
josias