[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