open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

For working with random waves in TOMAWAC V6P2 ,I edited the LIMWAC.f file and included in the case file as "Fortran file". When I put wave height as 2m meter and a peak period of 0.1s, propagating from SW direction. But after the completion of simulation for 12 time steps with a interval of 1800 seconds shows decrease in wave height to almost zero.variation of wave also stops after some time steps.It seems that temporal variation of wave is not working.Here I am attaching the all files related this simulation. Help to solve this issue.

Thanks in advance


Regards,
Glejin.

Hello,

did you check the reading of your wave conditions in your fortran file? I don't see where the time dependence in the reading appears in your limwac.f routine.

As you didn't send your parameter file, I don't think I can help further.

Regards,

Elodie

Here I am attaching the parameter file casex1.txt.

I am looking for your reply.

Thank you

Glejin.

I will try to have a look at it but first I'm waiting for your answer concerning the reading of your boundary conditions file : "Did you check the reading of your wave conditions in your fortran file? I don't see where the time dependence in the reading appears in your limwac.f routine."

I have wave data set for half an hour interval at the Boundary. So I put wave data at time steps of 1800 seconds in the boundary condition file (tomawac_b_limit3.txt). I directly read the boundary condition at each time step from the boundary file and assign the values to HM0L,FPCIL,TETA1L and passed to SPEINI.f.

Ok, I understand, but this is not exactly my question. You need to check that your wave parameters are well updated at each time step by your limawc routine. Did you try to write them at each time step to check this? What does it give you for your 10 boundary points?
Then, when you plot your wave height, what are the results at each of the boundary points?

Elodie

Hello Elodie,

I did write the boundary wave parameters at each time step in the terminal. It shows that boundary wave condition is changing each time. I am attaching both modified part of the limwac.f and screen shot of boundary conditions changing with each time step. But the output file shows very small values of significant wave height.


MODIFIED PART OF LIMAWAC

IF (SPEULI) THEN

READ(26,*)TT(jj,1),HM0L1(jj,1),FPICL1(jj,1),THETA1(jj,1)

HM0L=HM0L1(jj,1)
FPICL=FPICL1(jj,1)
TETA1L = THETA1(jj,1)
WRITE(*,*),TT(jj,1),HM0L,FPICL,TETA1L,AT,LT

IF (NPB.EQ.NPTFR) THEN
CALL SPEINI
&( FBOR , SPEC , FRA , UV2D , VV2D , FREQ ,
& TETA , GRAVIT, FPMAXL , FETCHL, SIGMAL, SIGMBL, GAMMAL, FPICL,
& HM0L , APHILL, TETA1L , SPRE1L, TETA2L, SPRE2L, XLAMDL,
& NPTFR , NPLAN , NF , LIMSPE, E2FMIN, PROF , FRABL )
ELSE
CALL SPEINI
&( FB_CTE, SPEC , FRA , UV2D , VV2D , FREQ ,
& TETA , GRAVIT, FPMAXL , FETCHL, SIGMAL, SIGMBL, GAMMAL, FPICL,
& HM0L , APHILL, TETA1L , SPRE1L, TETA2L, SPRE2L, XLAMDL,
& NPB , NPLAN , NF , LIMSPE, E2FMIN, PROF , FRABL )
ENDIF


ENDIF

Hello,

I tried to run your case but it stopped immediately because of a segmentation fault (before that I had to correct the parameter file because parentheses around 'true' were not accepted). Do you manage to run this case without any problem?
I am sorry but I cannot spend time debugging your case, I can only try to give you some advice.

Regards,

Elodie

Hi Elodie,

I used Cygwin to run TOMAWAC in windows 7 system. I successfully run this case with out any error (with and without parantesis ' '). Here I am attaching the cas file with out paranthesis.
I think it would be much appreciable if tomawac include the temporal and spatial variation with in the source code. It will help people who are not handy with Fortran code.

Regards ,

Glejin

glejin wrote:
Hi, this may be irrelevant to your problem, but a wave with a height of 2.0m and peak period of 0.1s is VERY steep. I don't know what are the objectives of your test, but this is stretching the boundaries of wave mechanics!

Costas