[MITgcm-devel] problem with SEAICEadvSnow

[Dimitris Menemenlis]

Martin, welcome back.  Where have you been?  We have missed you!

> I cannot remember having any problems with my 2deg model

Ice and snow seem OK most everywhere.  Only a problem is one region right along
Antarctic coastline, maybe because of strong katabatic winds?

> However, I normally use advection scheme 1 (1st order upwind) and not 
> diffusion (DIFF1 = 0.)

Is advection scheme 1 the default?  If not, how would I set it?  Is there a 
particular advection scheme that would be more stable than the default that I 
can try?

Do you set DIFF1 to zero both for ice and snow?  Do you do that by commenting 
lines in the code or is there parameters to set in data.seaice?

> Why should the snow get unstable and the ice not, they usually have similar
> distributions (in my case).

I will take closer look at ice too and get back if I see something.  Maybe ice 
does not blow up because of 10-m cap?  The integration completes normally for 
"SEAICEadvSnow=.false."

On a different topic, please talk with JM and Chris and let me know how you 
would like me to respond to Elizabeth.

D.