open TELEMAC-MASCARET The mathematically superior suite of solvers

and of course, you can call speini from limwac to build your spectra at suited boundary points.

Hello Laurent

thanks for your replies. I attached my fortran file to have a look. In limwac i set an array with 3 points and 3 wave height (one for each boundary point).
The I call speini.f into a do loop for those points only.
My spectrum option is 7. When I check those points at bluekenue the wave height is not what i prescribed in the .f file

best regards,
Nikos

Hi Nikos,

Here is an example of limwac subroutine that you can use. No need to modify speini here.
I adapted it from "bottom_friction" test case with Jonswap spectrum instead of TMA (less parameters to define). You can run it with default geometry and boundary conditions files), and you'll normally see that you'll have the right values at the right boundary points.






I hope it helps !

Regards,
Laurent

Dear all, I am trying to simulate wave height which is constant over time but variable in space. I modified the limwac subroutine following the bottom friction test case, however I read the wave height,frequency and direction values from an external file since the values are many to be prescribed in the code. However i get a segmentation fault error in the first iteration. Could you help me locate the error?

Kind regards,
Andreas

Hello

What is the message of the error. Could you post the output of the code ?

regards.

Of course, here is the output of the code

You call speini with FB_CTE but FB_CTE is not allocated if speuli. See line 160 of your limwac : npb=nptfr.
Try to use fbor instead of FB_cte in your call to speini.

And The last three loops of limwac are useless. If you call speini with fbor, the 3 loops arounf line 334 is what you want i guess.

Regards.

Thank you very much Thierry for your suggestions.

I had to move the reading of the formatted file under the user modifications section (like the example of the bottom friction case posted above by Lebaleur) however i notice some abnormalities.

1) Some boundary nodes have zero wave height values although I think they are defined in my txt file (attached nodes.jpg). Are these nodes a part of a loop where the wave energy is set to zero?

2) Some wave directions are wrong (between the arrows in dir.jpg).According to the manual the direction is the one waves propagate to. But since the spectra is built for every boundary point and speini is called for every point direction is it calculated with respect to the local point orientation?

Once again your help is greatly appreciated,

Kind regards,
Andreas

I guess you made an error with the files you joined the limwac procedure does not take into account what i told you yesterday about F_CTE and FBOR. you still have your F_tmp....

Dear Mr. Fouquet,

I implemented the modifications you adviced yesterday, however i noticed that the wave height at the boundary is constant in space (only the last record of the wav1.txt file is kept and also the problems mentioned in the previous post are present (wrong direction and not good prescription on the nodes).

So i noticed that if i use the previous setup (ommiting the Limit spectrum modified by user keyword and using the previous limwac subroutine at least the wave height seems to be prescribed correctly.

Here is the limwac.f with the modifications you suggested though

Thank you very much for the assistance
Andreas