[MITgcm-devel] (not so) funny things happen in seaice_lsr and pickups

[Jinlun Zhang]

Martin, so a non-zero ZMIN helps to stabilize the system. What value do 
you use to make it stable?

Jinlun

Martin Losch wrote:
> The question marathon continues ...
> Jinlun, my observation is, that the metric terms cause all the trouble 
> (no problems so far in cubed sphere integration of any resolution, nor 
> on cartesian grid).  I have the suspicion that the symmetry of the 
> laplacian operator is broken, because the metric terms are distributed 
> unevenly between left and right hand side of the u and v equations. In 
> order to ensure the diagonal dominance of the coefficient matrix, 
> would it be possible to move all metric terms to the rhs? Just 
> wondering ...
>
> BTW, the system seems to be stable when I have a non-zero ZMIN.
>
> Martin
>