[MITgcm-devel] rotate_uv2en.F

[Jean-Michel Campin]

Hi Martin,

I did not see yet Gael answer.
I would be in favor of doing this weighting outside,
(otherwise, would need an other argument to switch this
off, since there are some fields for which we don't want
hFac weighting). But it's not a strong opinion.
Cheers,
Jean-Michel


On Wed, Aug 18, 2010 at 10:22:04AM +0200, Martin Losch wrote:
> 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
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel