--- MITgcm_contrib/articles/ceaice/ceaice_model.tex 2008/02/28 20:45:29 1.5 +++ MITgcm_contrib/articles/ceaice/ceaice_model.tex 2008/04/29 14:03:31 1.8 @@ -12,7 +12,7 @@ C-grid variants are available; the C-grid code allows for no-slip and free-slip lateral boundary conditions; \item two different solution methods for solving the nonlinear - momentum equations have been adopted: LSOR \citep{zha97}, EVP + momentum equations have been adopted: LSOR \citep{zhang97}, EVP \citep{hunke97}; \item ice-ocean stress can be formulated as in \citet{hibler87}; \item ice variables are advected by sophisticated advection schemes; @@ -98,7 +98,7 @@ \left(\dot{\epsilon}_{11}^2+\dot{\epsilon}_{22}^2\right) (1+e^{-2}) + 4e^{-2}\dot{\epsilon}_{12}^2 + 2\dot{\epsilon}_{11}\dot{\epsilon}_{22} (1-e^{-2}) - \right]^{-\frac{1}{2}} + \right]^{\frac{1}{2}} \end{align*} The bulk viscosities are bounded above by imposing both a minimum $\Delta_{\min}=10^{-11}\text{\,s}^{-1}$ (for numerical reasons) and a @@ -280,7 +280,7 @@ algorithm following Archimedes' principle) turns snow into ice until the ice surface is back at $z=0$ \citep{leppaeranta83}. -Effective ich thickness (ice volume per unit area, +Effective ice thickness (ice volume per unit area, $c\cdot{h}$), concentration $c$ and effective snow thickness ($c\cdot{h}_{s}$) are advected by ice velocities: \begin{equation}