[MITgcm-devel] seaice lsr

[Dimitris Menemenlis]

Go ahead.  I am assuming that later this week we will change a whole bunch of other things in lab_sea as well, so output will change yet again.  How does above change impact funnel experiements and arching?  D.

Dimitris Menemenlis
cell: 818-625-6498

-----Original Message-----

From:  Martin Losch <Martin.Losch at awi.de>
Subj:  [MITgcm-devel] seaice lsr
Date:  Tue Sep 25, 2007 11:14 am
Size:  1K
To:  MITgcm-devel at mitgcm.org

Hi there,

I have found that the LSR solver is a little inconsistent (which  
causes it to produce a slight assymetry with a symmetric domain and  
forcing):
the following lines:
>           URT(I,J)=FXY(I,J,bi,bj)+AA3
>      &          +UVRT1(I,J,bi,bj)*uIce(I,J-1,1,bi,bj)
>      &          +UVRT2(I,J,bi,bj)*uIce(I,J+1,1,bi,bj)
(and the corresponding ones for the vice equation) within a J-loop  
cause URT to be computed from the updated uIce(I,J-1,1,bi,bj) and the  
not yet updated uIce(I,J+1,1,bi,bj) (=uIce(I,J+1,3,bi,bj) at this  
time). One could replace these lines by:
>           URT(I,J)=FXY(I,J,bi,bj)+AA3
>      &          +UVRT1(I,J,bi,bj)*uIce(I,J-1,3,bi,bj)
>      &          +UVRT2(I,J,bi,bj)*uIce(I,J+1,3,bi,bj)
which changes the results (lab_sea fails with cg2d=10), but then the  
solution is symmetric, when it should be. Also this formulation  
removes a dependency from the j-loop which impairs vectorization. Can  
I change that?

Martin

_______________________________________________
MITgcm-devel mailing list
MITgcm-devel at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-devel