open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I have run some TOMAWAC simulations on a large scale grid (a bay), and I would like to focus on a smaller grid (near the shore).
I want to get back the calculated spectrum of the previous simulations (large scal) near the shore, and to impose it on the boundary of my new grid.
Can I do that ? What is the process, do I have to change the LIMWAC routine ?

Regards,
Rémi

Hello,

Apologies for the delay of the answer.
I guess that your finer model grid is nested into the large scale model grid (i.e. the boundary nodes of your finer model correspond to internal nodes of the coarser model).
The steps you should undertake are the following:

• In the large scale model : save into the "PUNCTUAL RESULTS FILE" the spectra associated to the nodes of the larger model that correspond to the boundary nodes of the finer mesh.
• In the finer model : modify LIMWAC in order to read the boundary conditions from this file (defined by the keyword BINARY FILE 1 in the steering file) and impose them as boundary conditions.

The difficulty consists of referring to the correct finer mesh node numbers when you read the spectra from the binary file and impose them as boundary conditions in the finer model.
Unluckily, there is not an automatic and easy to implement technique to do so: I am afraid you have to list "manually" the boundary node numbers and to assign to them the corresponding spectra of the binary file obtained from the larger model.

Hope it helps.
Kind regards,

Giovanni

Hello,

Thank you for your answer, and sorry for the delay of mine, too !


- Regarding the first step you adviced me, I did it. I saved into "spectre.spe" (PUNCTUAL RESULTS FILE) the spectra associated to the nodes of the large model.
- Unfortunatly the second step is more complicated, and I'm afraid I do not know how to do it. I took a look to the LIMWAC.f fortran file. But I have no idea how I can change it in order to read my "spectre.spe" file. Do I have to use a "read" function ? How can I allocate the new spectra parameters ?

Thank you for your time,

Regards,
Rémi

Hello,

the PUNCTUAL RESULTS FILE has a SERAFIN binary format.
Have a look at the Appendix 8 of the TOMAWAC guide to have some information about the structure of the file.
You will find also attached a file in which I copied a part of a fortran file used to read SERAFIN file.

Hope it helps,

Giovanni

Hi,

As we discussed previously, I changed the Limwac subroutine in order to input a spectra file (serafin format) at the boundary of a smaller domain.

It works quite well with my case files. So I attached to this message the limwac modifications (some of the comments are in French, sorry..).

I’ve deleted some parts of the initial code (of limwac.f), so it works only if you have a BINARY FILE 1 (ie. spectra.spe) to set the boundary conditions.

If you find any mistakes or modifications that need to be done, let me know.

Thanks’ again for your help,
Regards,
Rémi

Hi Remi

I used your hf22b4cd.f and modified it.
Thanks for that.

Now, it works for nesting tomawac varying ins pace and time.
you need to save the spectra at some poitns close to boundary from the coarser grid. and then run the finer grid and read those spectra output. in the limwac you need to specify, the totla number of times steps existing in the spectra output file and the coordinates of the points.

I hope this helps someone.

Homayoon

Hello
I've tried to use the limwac.f subroutine (attached in the previous post), and my problem consists in the fact that when the subroutine lit.f read the BINARY FILE 1 (my spectrum file) i have this error:

LIT : UNKNOWN STANDARD:

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2



pourriez vous m'aider? (je vous joins le fichier ires.spe)
Je vous remercie d'avance
Vassilis

Hi vasilis

which version of tomawac are you using? that one was for version v7p0.
for the new version use the attached.

and dont forget to modify the coordinates of the spectral output points.

Thank you for your help , unfortunately im still getting the same error

LT 0
AT 0.0000000000000000
limwac 231
reading spectra from coarse grid for first time
LIT : UNKNOWN STANDARD:

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2


i run the code on ubuntu 14.04, is that the cause of the problem?

i attached my files