[MITgcm-devel] problem with 59c

[Dimitris Menemenlis]

Martin, more on the 59c problem:

The "SEAICEuseDynamics=.false." problem is not "directly" related to the 59c 
crash.  That is, "SEAICEuseDynamics=.false." crashes in Arctic Ocean even with 
pre-59c code.

What causes trouble in 59c is the change from ZMIN=4e8 to ZMIN=0.

Figure http://ecco2.jpl.nasa.gov/data1/cube/cube66/crash.pdf
shows sea ice thickness around Antarctica right before crash.  If you blow up 
the pdf file, you will see that there are several regions with grid-scale noise 
in HEFF.  Also shown are time series of HEFF at and near region of maximum 
thickness showing that thickness reaches 2500 m before model crashes.

So what to do?  Go back to ZMIN=4e8?  Change sea ice advection scheme or 
diffusivity?  (I am using default advection scheme and diffusivity.) 
Artificially impose a limit on HEFF?  Something else?

D.