Swell

user_12498 · August 30, 2018, 7:20pm

Hello [name removed],

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

user_12498 · August 30, 2018, 7:36pm

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

user_97 · August 31, 2018, 7:40am

Hi,

You can have a look at the following topics for instance:

http://www.opentelemac.org/index.php/kunena/19-tomawac/11026-error-using-limwac-subroutine?limitstart=0

http://www.opentelemac.org/index.php/kunena/19-tomawac/10888-wave-series-at-boundary-with-subroutine#30595

http://www.opentelemac.org/index.php/kunena/19-tomawac/11074-spatially-and-temporally-varying-wave-height-at-boundary?start=10#29194

[name removed]

user_12498 · August 31, 2018, 8:32am

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

user_97 · August 31, 2018, 8:52am

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.

[name removed]

user_12498 · August 31, 2018, 9:43am

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

tom.cas (3.98 KB)
fortran_modif2.f (11.1 KB)

user_97 · August 31, 2018, 9:56am

As mentioned in the error message, the end of your file ([group removed]) 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,))

[name removed]

user_12498 · August 31, 2018, 1:58pm

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

fortran_modif2_2018-08-31.f (11.2 KB)

user_97 · August 31, 2018, 2:25pm

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.

[name removed]