open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everybody,

I would like to propagate a swell, like a wave height signal, which shows a growing wave height and then a decrease 'til 0 value.

May anyone know how to put this information in Tomawac?

Thanks everyone

Sorry. I saw that there is the possibility to modify LIMWAC subroutines and write a txt file in order to make wave height variable in time. Is there one example, please?

Thanks again

Hi,

You can have a look at the following topics for instance:
www.opentelemac.org/index.php/kunena/19-...routine?limitstart=0
www.opentelemac.org/index.php/kunena/19-...ith-subroutine#30595
www.opentelemac.org/index.php/kunena/19-...ndary?start=10#29194

Regards,
Laurent

Thank you very much.

The last url has been very useful for me, because it's the only that works and I have just seen the differences between bottom friction example and the "modified" one with the new limwac subroutine.

One more question, please. I'm new to fortran, and I don't understrand well all the added lines. May someone explain me what did the creator of the new limwac want to simulate, please?

I've attached the lines. Thanks a lot

The other urls also work fine!

In the first IF statement (corresponding to the first time step), 3 boundary points are defined (nodes n°2, 6 and 16) for which a significant wave height value is associated (hm_L).
Then, the 2D spectrum (frequency/direction) is built by calling speini and filled in FB_CTE variable.
F_tmp is an extra variable, not sure is useful for you, depending on the following of the subroutine.

Regards,
Laurent

Ok, thanks a lot.

I'm trying to run every simulation with the files in the posts which you referred to.

I've attached the steering file, the fortran file (both are quite the same than the two of one of the url you sent me) and the error pic.

Do you know the reason why the simulation does not work?

Thanks again for willingness

As mentioned in the error message, the end of your file (waves_lim.txt) is reached during its reading in limwac.
Try to comment out the "write" statements after the "read" ones to see if the reading is correctly done.

(By the way, replace the "read(26,*) by read (NFO1,*))

Regards,
Laurent

It doesn't work even if I modified the fortran file as you told me before.

Might you give me any other suggestions, please?

Thanks a lot

You have to adapt the limwac subroutine according to your wave data file (or the contrary). Is the reading of the header correct?
The number of records NO_TIME_POINTS seems not to be in agreement with the number of lines of your file (header excluded). You're trying to read more wave data than there are in the file.

Maybe you can begin with learning basic skills for fortran programming.

Regards,
Laurent