[MITgcm-devel] bug in seaice_ocean_stress?

[Dimitris Menemenlis]

>           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)?


OK. I will take a look when I get in to work and give it a try.  Could that be 
cause of striations?  D.