I’m working on a wave model for a lake and would like to know if it is possible to run TOMAWAC as a standalone model, since I want to compare the results with a SWAN model.
If this is possible how does one setup the model, e.g. initial water level?
I could run a model without coupling and this went well. However, what I can not get hold of is how to give the wind information to TOMAWAC.
I used to give them as an ASCII file to Telemac (Timestep, Direction, Windspeed). From reading the manual I could not figure out if this is possible for TOMAWAC as well, since the manual and the reference manual only mentions the formats 3 and 4, which seem not to be what I am searching for. I used v8p2r1 manual.
I looked at the example turning_wind which uses a .slf file as wind input, which is not really what I am searching for.
Introducing your wind input in ASCII format is principally possible via selecting option-4 (user format), but it still requires programming Fortran subroutines. I recommend running a TELEMAC2D simulation with ASCII wind input (see wind_txy in telemac2d examples), then use the results file as your TOMAWAC’s winds file (option-3). Alternatively, I would write a Python/MATLAB script to write wind data directly into a selafin format, which is far way more efficient.
Thanks for the help. I kind of assumed that I need to run the Telemac model first.
Writing a Python script which writes wind data into a selafin format sound much more efficient to me, however I have no clue how the selafin file is structured and may need to run a Telemac model anyway to get this structure. If you already have such a script I would be happy if you could share it.
You don’t need to know exactly how a selafin file is structured as there is python script already dedicated to such tasks…
Try to have a look at the existing scripts dealing with meteo data could give you some ideas to create your own script (for example those which convert ECMWF data to telemac atmospheric data file)
Agreed - I’ve also found writing a simple script to interpolate wind data onto the computational mesh in SERAFIN format is the most straight forward and effective way of implementing non-stationary winds. Personally I do this with the Telemac matlab toolbox tools.
I create a dummy SERAFIN file in BK with the wind variables to be read. Then use the telemac matlab functions to write the interpolated wind data onto the grid for each timestep.