[MITgcm-devel] pkg exf useRelativeWind and constant wind

[Jean-Michel Campin]

Hi Dimitris (and Gael),

There is a problem with useRelativeWind when wind is not reset by S/R EXF_SET_UV call.
This happens when either u/vwinfile is empty or when u/vwindperiod is zero
(constant in time). The only exception would be when exf-interp is used
for both components together (uvInterp_wind=T).
This is bad since the ocean surface velocity keeps accumulating (with minus sign)
into uwind/vwind.

For now, I am going to put a stop in exf_check.F but a clean way to do this 
useRelativeWind would be to have specific arrays (relWindX,relWindY) that contain the
relative wind component (if useRelativeWind) or just a copy of uwind,vwinf (it not).
Would also simplify useRelativeWind code in pkg/seaice (+ cleaner diagnostics).

There is a simlar problem with rotateStressOnAgrid (might be easier to fix ?)
but for now I am going to add an other stop there.

Cheers,
Jean-Michel