--- MITgcm_contrib/articles/ceaice/ceaice.tex 2008/02/25 19:30:56 1.11 +++ MITgcm_contrib/articles/ceaice/ceaice.tex 2008/02/25 22:06:17 1.12 @@ -1,4 +1,4 @@ -% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/articles/ceaice/ceaice.tex,v 1.11 2008/02/25 19:30:56 mlosch Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/articles/ceaice/ceaice.tex,v 1.12 2008/02/25 22:06:17 dimitri Exp $ % $Name: $ \documentclass[12pt]{article} @@ -54,22 +54,28 @@ \begin{abstract} As part of ongoing efforts to obtain a best possible synthesis of most -available, global-scale, ocean and sea ice data, dynamic and thermodynamic -sea-ice model components have been incorporated in the Massachusetts Institute -of Technology general circulation model (MITgcm). Sea-ice dynamics use either -a visco-plastic rheology solved with a line successive relaxation (LSR) -technique, reformulated on an Arakawa C-grid in order to match the oceanic and -atmospheric grids of the MITgcm, and modified to permit efficient and accurate -automatic differentiation of the coupled ocean and sea-ice model -configurations. +available, global-scale, ocean and sea ice data, a dynamic and thermodynamic +sea-ice model has been coupled to the Massachusetts Institute of Technology +general circulation model (MITgcm). Ice mechanics follow a viscous plastic +rheology and the ice momentum equations are solved numerically using either +line successive relaxation (LSR) or elastic-viscous-plastic (EVP) dynamic +models. Ice thermodynamics are represented using either a zero-heat-capacity +formulation or a two-layer formulation that conserves enthalpy. The model +includes prognostic variables for snow and for sea-ice salinity. The above +sea ice model components were borrowed from current-generation climate models +but they were reformulated on an Arakawa C-grid in order to match the MITgcm +oceanic grid and they were modified in many ways to permit efficient and +accurate automatic differentiation. This paper describes the MITgcm sea ice +model; it presents example Arctic and Antarctic results from a realistic, +eddy-permitting, global ocean and sea-ice configuration; it compares B-grid +and C-grid dynamic solvers in a regional Arctic configuration; and it presents +example results from coupled ocean and sea-ice adjoint-model integrations. \end{abstract} \section{Introduction} \label{sec:intro} -more blabla - \section{Model} \label{sec:model} @@ -367,13 +373,18 @@ \subsection{Arctic Domain with Open Boundaries} \label{sec:arctic} -The Arctic domain of integration is illustrated in Fig.~\ref{???}. It +The Arctic domain of integration is illustrated in Fig.~\ref{fig:arctic1}. It is carved out from, and obtains open boundary conditions from, the global cubed-sphere configuration of the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) project \citet{menemenlis05}. The domain size is 420 by 384 grid boxes horizontally with mean horizontal grid spacing of 18 km. +\begin{figure} +%\centerline{{\includegraphics*[width=0.44\linewidth]{\fpath/arctic1.eps}}} +\caption{Bathymetry of Arctic Domain.\label{fig:arctic1}} +\end{figure} + There are 50 vertical levels ranging in thickness from 10 m near the surface to approximately 450 m at a maximum model depth of 6150 m. Bathymetry is from the National Geophysical Data Center (NGDC) 2-minute gridded global relief