[MITgcm-devel] problem with SEAICEadvSnow

[Dimitris Menemenlis]

Martin,

> Do you actually use the 10m cap for heff? This is what I find in the current
> code (and this part has not change for a long time):

I was looking at this line:

> C     Capping the actual ice thickness effectively enforces a
> C     minimum of heat flux through the ice and helps getting rid of 
> C     very thick ice.
>           HICE(I,J)    = MIN(HICE(I,J),9.0 _d +00)

We cap sea-ice thickness prior to calling SEAICE_BUDGET_ICE.  Not sure what 
impact of this will be on HEFF.  I will take a look at model output prior to 
blow-up when I get in the office.

> SEAICEadvScheme   = 1, ! this is not the default. The default is 2 (2nd order
> central differences with modified euler time step, the original code)

Is "SEAICEadvScheme = 1" supposed to be more stable than "2"?

D.