[MITgcm-devel] verification/seaice_obcs

[Jean-Michel Campin]

Hi Dimitris,

On Mon, May 12, 2008 at 06:03:59PM -0700, Dimitris Menemenlis wrote:
> >Not sure where the problem is.
> >e.g., using obcs:
> >uwind(i=Nx) has 2 center neighbours XC(Nx-1) & XC(Nx) which should  
> >receive
> >valid values (and unaffected by the exchanges) from exf
> >and it should not be a problem to interpolate the wind to this  
> >location.
> 
> The funny edge winds are caused by vector rotation in exf_set_uv.F,  
> which requires valid xG and yG values at Nx+1 and Ny+1.
> 

This might be a good reason to switch to the pre-defined
angle rotation, stored in arrays (angleCosC, angleSinC)
and corresponding to the center (tracer point) of the grid cell.

Jean-Michel