wave series at boundary with subroutine

System Modules TOMAWAC
1

Hi,

I’m trying to set up a tomawac model and my wave data are in the format of time, Hs, Tp, and direction. I use a subroutine of someone who resolve that kind of case but I still have an error.

runcode::main:

:

|runCode: Fail to run

|C:\Users\asus\Desktop\STAGE_KWATER\DOC\TOMAWAC_eureka\17\1\model_3_eureka_unsteady.cas_2017-06-27-20h14min00s\out_princi.exe

|~~~~~~~~~~~~~~~~~~

|At line 116 of file C:\opentelemac-mascaret\v7p0\sources\utils\bief\leclim.f (unit = 7, file = ‘[name removed]’)

|Fortran runtime error: End of file

|~~~~~~~~~~~~~~~~~~

Is there someone who had that error before ?I’ve attached my files.

Thanks

princi.f (10.9 KB)
model_2_eureka_unsteady.cas (3.03 KB)

2

Hello

Apparently from the message you got i would say that you do not have enough points in your boundary file. leclim is the routine that reads boundary condition and [name removed] is your

BOUNDARY CONDITIONS FILE

hope it helps.

T.

3

Thanks, it was a mistake in my waves data file. But I still have some trouble. It does compute but when I look at the results it seems that I have only the first value which is propagated but the new boundary values are not taken in consideration (like see on the gif).

Hs.gif
POW.gif

4

Hello

I read your princi.f. It would be better to use your HM0L FPICL and TETAL

after you had found the time that interest you. (L245-253)

you have to set the spectre you want to FB_CTE by

      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 )
5

It 's the first time that I have to use fortran and actually I’m lost. Whem I looked at the code the algorithm seems good to me. I add the function. As results, I have a segmentation fault.

                 *************************************

                 *    END OF MEMORY ORGANIZATION:    *

                 *************************************

INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)

STRCHE (BIEF): NO MODIFICATION OF FRICTION

TOM_CORFON : NO MODIFICATION OF BOTTOM

================================================================================

   TEMPS :       0.0000 SECONDES                            0IEME  ITERATION

########### 0 1.6299999999999999 0.<CREDIT_CARD> 90.000000000000000

runcode::main:

:

|runCode: Fail to run

|C:\Users\asus\Desktop\STAGE_KWATER\DOC\TOMAWAC_eureka\17\2\new_try.cas_2017-06-29-14h23min42s\out_fortran_modif2.exe |~~~~~~~~~~~~~~~~~~

|Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

|

|Backtrace for this error:

|~~~~~~~~~~~~~~~~~~

fortran_modif2.f (11.1 KB)

6

Actually, it’s working. I called the wrong boundary condition file in my steering file.

Thanks for your help :wink:

7

Hi, i’m trying to do the same as you but I dont really understand what should I modify in the LIMWAC subroutine ? My wave data are also time, Hmo, Tp and direction.

Can you help me please ?

8

Hi,

In limwac subroutine, you have to read your data at the right time step and fill the arrays HM0L, FPICL and TETA1L.

Then, for each boundary points, you call speini with these parameters to build the spectra (stored in FBOR array) according the TYPE OF BOUNDARY DIRECTIONAL SPECTRUM you want.

I hope it helps.

[name removed]

9

Hi, thank you for this answer. Unfortunately this is still really unclear for me.

If I have my datas in a text file, how can I read them at the right time step and fill the arrays ?