FORTRAN peculiarity in VENUTI.f

user_8512 · August 28, 2018, 10:03am

Hello,

I am using [name removed] in TOMAWAC in order to include (time-dependent) wind effects in the wave action equation. So I modified the file so that it reads a file with three columns:

t u10 v10

and assign these values to the i/o arrays U1,U2,V1,V2 for the two first time steps.

Strangely, after having assigned U1 and V1 the assignment of U2 and V2 seems to change the value of U1 and V1. In attachment a screenshot of my playing-around.

That is absolutely uncanny and illogical. And in the end, my U1,U2 and V1,V2 are equal. So, not what I precisely want.

Does anybody have an explanation?

[name removed]

user_63 · August 28, 2018, 10:25am

Hi

really strange but hard to analyse without the full test case to test on our side…

As a trial I would test to change the call OV by a local loop…

Regards

user_8473 · August 28, 2018, 10:50am

Hello

At the beginning (t=0) there are only one time old=actual. So when venuti is called by [name removed] you will see that U1 and U2 are pointing to the same vector (F1,F2 and F1,F2)

CALL [name removed],Y,NPOIN2,NDON,FFORMAT,NBOR,NPTFR,

 &                0.D0,0.D0,T1VOID,T2VOID,F1,F2,F1,F2)

That is not the case for next time step when noudon calls vebuti

CALL VENUTI(X,Y,NPOIN,NDON,FFORMAT,NBOR,NPTFR,AT,DDC,

 &                  TV1,TV2,F11,F21,F12,F22)

In your print i guess that if you write [name removed] and TV2 too, you will get the same value you are just printing your wind twice that the time 0.

So i would check your printing for iteration 1 or 2

hope it helps

T.