open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

Interested to know if anybody has attempted to modify the source code to read non-stationary wave files, similar to how the atmospheric data files are read and interpolated. As i understand, the currently functionality only considers the last timestep of the wave file.

Being more familiar with the SWAN and it with the benefit of its computational speed (perhaps at the expense of decreased accuracy) I would like to create a wave file with radiation stresses output from SWAN to drive wave induced currents in TELEMAC - this would be simple enough to construct.

Any help would be greatly appreciated - under a tight timeframe and unfortunately not familiar enough with the source to make these edits from scratch at the moment.

Thanks :)

I believe I've managed to get this to work through relatively simple modifications to the 'PROSOU' subroutine by adapting the code used to interpolate the binary winds file with the 'READ_BIN_2D' subroutine.

The SWAN model outputs radiation force (FX,FY) in units of pressure (N/m^2) or Force/unit surface area, while TOMWAC provides it in units of acceleration (m/s^2). Through some simple algebra of equating these two forces I believe that the SWAN outputs should be divided by the water density adopted in SWAN to arrive at the units of acceleration required by TELEMAC.

Comparing direct TOMAWAC and SWAN radiation force outputs, the differences were of the order of 1000 which gives me confidence that this is the correct approach. Moreover, the model produced currents are of the order I'd expect.

Has anyone attempted this before and if so are my assumptions here correct?

Hi,
I am also interested in this. Would you kindly share the 'PROSOU' subroutine? Thank you so much in advanced.

Best,
Htun

Hi,

Not a problem at all.

Attached is the modified TRISOU subroutine. This is the T3D version of PROSOU but the required modifications between the two will be much the same. Lines commented with '!TJN' are where the edits have been made.

Hope this helps

Regards,

Thank you very much. It is very helpful.
Regards,
Htun.