Thank you for this information about strong currents.
The current here is not that large... but I can have a try.
In parallel, going back to my initial question, I would also test the following change in CORMAR.f subroutine in order to set current velocity to zero at each time step:
IF(PART.EQ.1) THEN
CALL OV('X=C ',UC, U_TEL%R, U_TEL%R, 0.D0, NPOIN2)
CALL OV('X=C ',VC, V_TEL%R, V_TEL%R, 0.D0, NPOIN2)
ENDIF
instead of:
IF(PART.EQ.1) THEN
CALL OV('X=Y ',UC, U_TEL%R, U_TEL%R, 0.D0, NPOIN2)
CALL OV('X=Y ',VC, V_TEL%R, V_TEL%R, 0.D0, NPOIN2)
ENDIF
Is this change correct and harmless for the rest of the TOMAWAC computations?
I think this would do the job while the sea level is still varying and should permit to emphasize the effect of currents on Hm0 by comparing with my previous results.
JBS.