open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all,
I think I have understood how to modify the limwac subroutine to read temporally varying wave height at the boundary,however I can't find a test case where a spatially and temporally varying wave height is imposed. Is there an example of a modified limwac.f subroutine that fits this criteria? Or else could you provide me with some directions.I am thinking that the limwac.f must be modified to work like the meteo.f subroutine with regards to the boundary points.

Thank you for your assistance,
Kind regards,
Andreas

Hello

You will find examples of modified limwac in the test cases. For example, in case Manche there is some reading in limwac.

regards

Thierry

Thank you very much.
Could you possibly pinpoint where can I find the manche test case?

In the directory examples/tomawac you will find all the test cases and among them you will find Manche test

regards

T.

Hello,

I am also looking for a modified limwac.f routine to read time- and space-varying wave boundary conditions.

I have been looking into the latest versions (v7) of opentelemac Repository, but unfortunately, I did not find the mentionned "Manche test".

What am I missing?
Thanks for your help.

Regards,
Claire

Hello Claire

Sorry i didn't pay attention on the number of the official version, you will find the test Manche in the trunk version or the V7P3 version

svn.opentelemac.org/svn/opentelemac/bran...ples/tomawac/Manche/

regards

T.

Hello Thierry and community

In the Manche test there is a *.lar file that contains the spatio temporal information. Is there any heading for the information provided?
I suspect it is information about the spectrum but not sure.
For my case I need to provide wave height, period and direction and I am not sure what array I need to use?. In the specific example the values from the file are placed into the variables CL1 and CL2 which are 3 dimensional (directions,frequencies,number of nodes);

Also another thing is the definition of time. At line 3 of *.lar file the number is the same given at steering file after the keyword "DATE OF COMPUTATION BEGINNING". However the manual says that this keyword is of yyyymmddhhmm

Many thanks and regards

Hello,

Indeed, the *.lar file contains the 2D density spectra (frequency/direction) for 21 particular boundary points (third dimension of CL1 and CL2), and an interpolation is done after between couple of spectra for the other boundary points.
In this example, CL1 and CL2 arrays contain energy densities for 2 consecutive time steps. A time interpolation is also done.
The time format is transformed via the TEMP subroutine (at the end of manche.f) to be compliant with the standard format of the steering file.
You can find more details about the time format and a general description of test case here:
v7p3r0/documentation/tomawac/validation/tomawac_validation_v7p3.pdf
see section 5

For your case, you have to introduce your time varying parameters (Hs, Tp, direction) in the limwac subroutine too, before the call to speini. From this parameters, the speini subroutine will then build the spectra according to the TYPE OF BOUNDARY DIRECTIONAL SPECTRUM that you have chosen. The same spectrum will be imposed over your whole liquid boundary.

I hope it helps.

Regards,
Laurent

Hello Laurent,

Thanks for your reply. In my case, I have different values for each boundary point. Subroutine SPEINI loops over all points of the mesh. However, I want calculation on specific boundary points. How can I pass that information from limwac.f to speini.f?

Many thanks and regards
Nikos

Hi Nikos,

You have first to turn on the keyword "LIMIT SPECTRUM MODIFIED BY USER" in your steering file.
Then in limwac, you have to modify the section from "IF(SPEULI) THEN" so as to fill the array FBOR(IBoundaryPoints,IFrequency,IDirection) for each boundary points with your 2D spectra, and (if necessary) interpolate the spectra between 2 boundary points.

Regards,
Laurent