[MITgcm-support] EVP stability

Dimitris Menemenlis menemenlis at sbcglobal.net
Sat Jan 27 09:36:55 EST 2007


Matt, stability criterion for EVP solver depends on internal ice pressure, 
strain rates, a damping parameter for the elastic waves, and other such things. 
  So it's not easy to compute explicitly.

In "Journal of Computational Physics 170, 18–38 (2001)", Elizabeth Hunke used a 
time step of 30 s for a domain with horizontal grid spacing of 16 km.  So 45 s 
for your 1/6-deg domain seems too large.  What is your smallest horizontal grid 
spacing, 4 km?  I suggest that you try a 5-s to 10-s time step.

Dimitris




More information about the MITgcm-support mailing list