[MITgcm-devel] bug in seaice_ocean_stress?

[Martin Losch]

Hallo,

I know it's my own code, but I think that there is a bug in the way  
the drag coefficient is interpolated to u/v points in  
seaice_ocean_stress (if useHB87stressCoupling = .false., which is the  
default).

           fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*
should be
           fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) )*

and
           fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*
should be
           fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) )*

if dwatn is defined at  C/Mass-points (which they appear to be). Do  
you agree?

Dimitris, could you try this in your Arctic configuration (with and  
without EVP) (before I check it in)?

Martin