54 |
\begin{abstract} |
\begin{abstract} |
55 |
|
|
56 |
As part of ongoing efforts to obtain a best possible synthesis of most |
As part of ongoing efforts to obtain a best possible synthesis of most |
57 |
available, global-scale, ocean and sea ice data, dynamic and thermodynamic |
available, global-scale, ocean and sea ice data, a dynamic and thermodynamic |
58 |
sea-ice model components have been incorporated in the Massachusetts Institute |
sea-ice model has been coupled to the Massachusetts Institute of Technology |
59 |
of Technology general circulation model (MITgcm). Sea-ice dynamics use either |
general circulation model (MITgcm). Ice mechanics follow a viscous plastic |
60 |
a visco-plastic rheology solved with a line successive relaxation (LSR) |
rheology and the ice momentum equations are solved numerically using either |
61 |
technique, reformulated on an Arakawa C-grid in order to match the oceanic and |
line successive relaxation (LSR) or elastic-viscous-plastic (EVP) dynamic |
62 |
atmospheric grids of the MITgcm, and modified to permit efficient and accurate |
models. Ice thermodynamics are represented using either a zero-heat-capacity |
63 |
automatic differentiation of the coupled ocean and sea-ice model |
formulation or a two-layer formulation that conserves enthalpy. The model |
64 |
configurations. |
includes prognostic variables for snow and for sea-ice salinity. The above |
65 |
|
sea ice model components were borrowed from current-generation climate models |
66 |
|
but they were reformulated on an Arakawa C-grid in order to match the MITgcm |
67 |
|
oceanic grid and they were modified in many ways to permit efficient and |
68 |
|
accurate automatic differentiation. This paper describes the MITgcm sea ice |
69 |
|
model; it presents example Arctic and Antarctic results from a realistic, |
70 |
|
eddy-permitting, global ocean and sea-ice configuration; it compares B-grid |
71 |
|
and C-grid dynamic solvers in a regional Arctic configuration; and it presents |
72 |
|
example results from coupled ocean and sea-ice adjoint-model integrations. |
73 |
|
|
74 |
\end{abstract} |
\end{abstract} |
75 |
|
|
76 |
\section{Introduction} |
\section{Introduction} |
77 |
\label{sec:intro} |
\label{sec:intro} |
78 |
|
|
|
more blabla |
|
|
|
|
79 |
\section{Model} |
\section{Model} |
80 |
\label{sec:model} |
\label{sec:model} |
81 |
|
|
373 |
\subsection{Arctic Domain with Open Boundaries} |
\subsection{Arctic Domain with Open Boundaries} |
374 |
\label{sec:arctic} |
\label{sec:arctic} |
375 |
|
|
376 |
The Arctic domain of integration is illustrated in Fig.~\ref{???}. It |
The Arctic domain of integration is illustrated in Fig.~\ref{fig:arctic1}. It |
377 |
is carved out from, and obtains open boundary conditions from, the |
is carved out from, and obtains open boundary conditions from, the |
378 |
global cubed-sphere configuration of the Estimating the Circulation |
global cubed-sphere configuration of the Estimating the Circulation |
379 |
and Climate of the Ocean, Phase II (ECCO2) project |
and Climate of the Ocean, Phase II (ECCO2) project |
380 |
\citet{menemenlis05}. The domain size is 420 by 384 grid boxes |
\citet{menemenlis05}. The domain size is 420 by 384 grid boxes |
381 |
horizontally with mean horizontal grid spacing of 18 km. |
horizontally with mean horizontal grid spacing of 18 km. |
382 |
|
|
383 |
|
\begin{figure} |
384 |
|
%\centerline{{\includegraphics*[width=0.44\linewidth]{\fpath/arctic1.eps}}} |
385 |
|
\caption{Bathymetry of Arctic Domain.\label{fig:arctic1}} |
386 |
|
\end{figure} |
387 |
|
|
388 |
There are 50 vertical levels ranging in thickness from 10 m near the surface |
There are 50 vertical levels ranging in thickness from 10 m near the surface |
389 |
to approximately 450 m at a maximum model depth of 6150 m. Bathymetry is from |
to approximately 450 m at a maximum model depth of 6150 m. Bathymetry is from |
390 |
the National Geophysical Data Center (NGDC) 2-minute gridded global relief |
the National Geophysical Data Center (NGDC) 2-minute gridded global relief |