[MITgcm-devel] rotate_uv2en.F

[Martin Losch]

Hi Gael,

I was quite happy to learn that there is already a routine that rotates vector: rotate_uv2en.F. This is exactly what I need for my attempts to couple MITgcm to ECHAM5!

I'd like to use it to interpolate and rotate ocean velocities to be used by the atmosphere. I have the impression that it is necessary to do the averaging like this:

u(on c) = (uvel(i,j)*hfacw(i,j)+uvel(i+1,j)*hFacw(i+1,j))*.5*recip_hfacc(i,j)

Do you think that this is also useful for your purposes and can we somehow include this in the routine, or should I rather do this outside?

Martin