[MITgcm-devel] bug in seaice_ocean_stress?

[Martin Losch]

I don't think so (although you never know), but I want to make sure  
that fixing this bug will not destroy everything we have right now.
M.

On 14 May 2007, at 16:01, Dimitris Menemenlis wrote:

>>           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.
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel