open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I am using VENUTI.f 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?

Thank you.
JBS.

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

Hello

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

CALL VENUTI(X,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 TV1 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.