open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all,

I like to simulate a wind which is various in time but not in space.
I try to include my windfile in my tomawac simulation using "winds file format" 4.

I have attached my windfile (*.txt)

Time(s), WindDirection (°), windAlongX (u, m/s),windAlongY (v, m/s)

my steering file(*.cas) and the venturi.f file.

Unfortunatly I am not fimilar with fortran at the moment. Can anyone help me, how I have to modify the venuti.f file?! Or is there any better solution for this issue?

Thanks for any help!

runkle

Hello Runkle,

I suppose you don't need wind direction if you have the UV vectors. Have a look at telemac's subroutine meteo.f. It does exactly what you want but with telemac's variables, so you have to adapt it for tomawac.

You can start from my venuti file, BUT pay attention that it is setup to read wind direction and speed and then convert to vectors! Your case is simpler. Follow meteo.f and you are done.

I hope it helps,
Costas

hey cyamin

thanks a lot, I will try and let you know....

regards,
runkle

Morning Cyamin,

General, I have wind data in format T, Dir, and Speed. I have calculated the components with a own matlab code. So, I prefer to use the venuti.f.

Your venuti script works but it read only the first time step and not for the entire simulation time. What have I to extend in the venuti file?

attached is my original windfile

Thanks for any help!

regards,
runkle

Hi Runkle,

venuti.f is just for reading the wind file every time it is called, so there shouldn't be any problem unless your computation starts at a time later than 0.

Maybe subroutine noudon.f (that calls venuti.f) could be the cause of your problem. I had issues my self at the past. Can you be more specific with your problem? You can turn debugger=1 and post your output?

Regards,
Costas

my simulation starts at Time = 0, see attached testwindfile and my casFile.
Further, I have attached my sortie file after set debugger=1.

simulations ends normal, but when i open my results file, the windvelo along xy is for the entire simulation periode 8.224 with the direction 0 (first row of my windfile). Unfortunatly, I have no idea at the moment, what I am making wrong!
Do you have any idea?

regards,
runkle

Hi Runkle,

Can you increase the duration of the computation so that it covers at least two entries in your wind file and check again?

Costas

I have tried follow thinks:

- timestep =800s and number of timestep = 10
- timestep = 400 and number of timestep = 15

but nothing helps

Further, my simulation covers all entries of the windfile:
T Dir w
(s) (°) (m/s)
0 0 8.224
400 0 9.766
800 40 6.168
1200 0 6.168
1600 0 8.224
2000 40 9.766
2400 60 7.196
2800 320 6.168
3200 0 9.252
3600 40 10.794
4000 80 7.196

a part of my *.cas
/
TITLE = 'mauri_Tomawac_v5'
TIME STEP = 400
NUMBER OF TIME STEP = 10
PERIOD FOR GRAPHIC PRINTOUTS = 1
PERIOD FOR LISTING PRINTOUTS = 1

runkle

I suppose noudon.f is your problem, because venuti.f works OK since it reads at least the first record.

I have made some changes some time ago when I had a similar problem with wind data not being read and interpolated correctly, especially in restarts when I had to parse the wind file at the first time step. I have found that (probably) noudon is not well coded when it comes to reading from user defined files, so I tried to correct this by replicating the approach to selafin wind files to user defined files.

A attach the changes to both noudon and venuti. Give it a try and see if it works for you. Please bear in mind that I haven't time to test it thoroughly and I am not sure whether it is a robust solution and that it does not cause trouble elsewhere.



I hope it helps,
Costas

unfortunatly, nothing changed. same problem. makes me crazy not to know, what the problem is because the first data line is read!

have u an example of your windfile?

regards,
runkle