37 |
% $Header$ |
% $Header$ |
38 |
% $Name$ |
% $Name$ |
39 |
|
|
40 |
\section{Introduction} |
This document provides the reader with the information necessary to |
|
|
|
|
This documentation provides the reader with the information necessary to |
|
41 |
carry out numerical experiments using MITgcm. It gives a comprehensive |
carry out numerical experiments using MITgcm. It gives a comprehensive |
42 |
description of the continuous equations on which the model is based, the |
description of the continuous equations on which the model is based, the |
43 |
numerical algorithms the model employs and a description of the associated |
numerical algorithms the model employs and a description of the associated |
47 |
both process and general circulation studies of the atmosphere and ocean are |
both process and general circulation studies of the atmosphere and ocean are |
48 |
also presented. |
also presented. |
49 |
|
|
50 |
|
\section{Introduction} |
51 |
|
\begin{rawhtml} |
52 |
|
<!-- CMIREDIR:innovations --> |
53 |
|
\end{rawhtml} |
54 |
|
|
55 |
|
|
56 |
MITgcm has a number of novel aspects: |
MITgcm has a number of novel aspects: |
57 |
|
|
58 |
\begin{itemize} |
\begin{itemize} |
87 |
computational platforms. |
computational platforms. |
88 |
\end{itemize} |
\end{itemize} |
89 |
|
|
90 |
Key publications reporting on and charting the development of the model are: |
Key publications reporting on and charting the development of the model are |
91 |
|
\cite{hill:95,marshall:97a,marshall:97b,adcroft:97,marshall:98,adcroft:99,hill:99,maro-eta:99}: |
92 |
|
|
93 |
\begin{verbatim} |
\begin{verbatim} |
|
|
|
94 |
Hill, C. and J. Marshall, (1995) |
Hill, C. and J. Marshall, (1995) |
95 |
Application of a Parallel Navier-Stokes Model to Ocean Circulation in |
Application of a Parallel Navier-Stokes Model to Ocean Circulation in |
96 |
Parallel Computational Fluid Dynamics |
Parallel Computational Fluid Dynamics |
99 |
Elsevier Science B.V.: New York |
Elsevier Science B.V.: New York |
100 |
|
|
101 |
Marshall, J., C. Hill, L. Perelman, and A. Adcroft, (1997) |
Marshall, J., C. Hill, L. Perelman, and A. Adcroft, (1997) |
102 |
Hydrostatic, quasi-hydrostatic, and nonhydrostatic ocean modeling, |
Hydrostatic, quasi-hydrostatic, and nonhydrostatic ocean modeling |
103 |
J. Geophysical Res., 102(C3), 5733-5752. |
J. Geophysical Res., 102(C3), 5733-5752. |
104 |
|
|
105 |
Marshall, J., A. Adcroft, C. Hill, L. Perelman, and C. Heisey, (1997) |
Marshall, J., A. Adcroft, C. Hill, L. Perelman, and C. Heisey, (1997) |
132 |
application to Atlantic heat transport variability |
application to Atlantic heat transport variability |
133 |
J. Geophysical Res., 104(C12), 29,529-29,547. |
J. Geophysical Res., 104(C12), 29,529-29,547. |
134 |
|
|
|
|
|
135 |
\end{verbatim} |
\end{verbatim} |
136 |
|
|
137 |
We begin by briefly showing some of the results of the model in action to |
We begin by briefly showing some of the results of the model in action to |
154 |
described in detail in the documentation. |
described in detail in the documentation. |
155 |
|
|
156 |
\subsection{Global atmosphere: `Held-Suarez' benchmark} |
\subsection{Global atmosphere: `Held-Suarez' benchmark} |
157 |
|
\begin{rawhtml} |
158 |
|
<!-- CMIREDIR:atmospheric_example --> |
159 |
|
\end{rawhtml} |
160 |
|
|
161 |
|
|
162 |
|
|
163 |
A novel feature of MITgcm is its ability to simulate, using one basic algorithm, |
A novel feature of MITgcm is its ability to simulate, using one basic algorithm, |
164 |
both atmospheric and oceanographic flows at both small and large scales. |
both atmospheric and oceanographic flows at both small and large scales. |
195 |
%% CNHend |
%% CNHend |
196 |
|
|
197 |
\subsection{Ocean gyres} |
\subsection{Ocean gyres} |
198 |
|
\begin{rawhtml} |
199 |
|
<!-- CMIREDIR:oceanic_example --> |
200 |
|
\end{rawhtml} |
201 |
|
\begin{rawhtml} |
202 |
|
<!-- CMIREDIR:ocean_gyres --> |
203 |
|
\end{rawhtml} |
204 |
|
|
205 |
Baroclinic instability is a ubiquitous process in the ocean, as well as the |
Baroclinic instability is a ubiquitous process in the ocean, as well as the |
206 |
atmosphere. Ocean eddies play an important role in modifying the |
atmosphere. Ocean eddies play an important role in modifying the |
227 |
|
|
228 |
|
|
229 |
\subsection{Global ocean circulation} |
\subsection{Global ocean circulation} |
230 |
|
\begin{rawhtml} |
231 |
|
<!-- CMIREDIR:global_ocean_circulation --> |
232 |
|
\end{rawhtml} |
233 |
|
|
234 |
Figure \ref{fig:large-scale-circ} (top) shows the pattern of ocean currents at |
Figure \ref{fig:large-scale-circ} (top) shows the pattern of ocean currents at |
235 |
the surface of a 4$^{\circ }$ |
the surface of a 4$^{\circ }$ |
248 |
%%CNHend |
%%CNHend |
249 |
|
|
250 |
\subsection{Convection and mixing over topography} |
\subsection{Convection and mixing over topography} |
251 |
|
\begin{rawhtml} |
252 |
|
<!-- CMIREDIR:mixing_over_topography --> |
253 |
|
\end{rawhtml} |
254 |
|
|
255 |
|
|
256 |
Dense plumes generated by localized cooling on the continental shelf of the |
Dense plumes generated by localized cooling on the continental shelf of the |
257 |
ocean may be influenced by rotation when the deformation radius is smaller |
ocean may be influenced by rotation when the deformation radius is smaller |
271 |
%%CNHend |
%%CNHend |
272 |
|
|
273 |
\subsection{Boundary forced internal waves} |
\subsection{Boundary forced internal waves} |
274 |
|
\begin{rawhtml} |
275 |
|
<!-- CMIREDIR:boundary_forced_internal_waves --> |
276 |
|
\end{rawhtml} |
277 |
|
|
278 |
The unique ability of MITgcm to treat non-hydrostatic dynamics in the |
The unique ability of MITgcm to treat non-hydrostatic dynamics in the |
279 |
presence of complex geometry makes it an ideal tool to study internal wave |
presence of complex geometry makes it an ideal tool to study internal wave |
293 |
%%CNHend |
%%CNHend |
294 |
|
|
295 |
\subsection{Parameter sensitivity using the adjoint of MITgcm} |
\subsection{Parameter sensitivity using the adjoint of MITgcm} |
296 |
|
\begin{rawhtml} |
297 |
|
<!-- CMIREDIR:parameter_sensitivity --> |
298 |
|
\end{rawhtml} |
299 |
|
|
300 |
Forward and tangent linear counterparts of MITgcm are supported using an |
Forward and tangent linear counterparts of MITgcm are supported using an |
301 |
`automatic adjoint compiler'. These can be used in parameter sensitivity and |
`automatic adjoint compiler'. These can be used in parameter sensitivity and |
316 |
%%CNHend |
%%CNHend |
317 |
|
|
318 |
\subsection{Global state estimation of the ocean} |
\subsection{Global state estimation of the ocean} |
319 |
|
\begin{rawhtml} |
320 |
|
<!-- CMIREDIR:global_state_estimation --> |
321 |
|
\end{rawhtml} |
322 |
|
|
323 |
|
|
324 |
An important application of MITgcm is in state estimation of the global |
An important application of MITgcm is in state estimation of the global |
325 |
ocean circulation. An appropriately defined `cost function', which measures |
ocean circulation. An appropriately defined `cost function', which measures |
326 |
the departure of the model from observations (both remotely sensed and |
the departure of the model from observations (both remotely sensed and |
327 |
in-situ) over an interval of time, is minimized by adjusting `control |
in-situ) over an interval of time, is minimized by adjusting `control |
328 |
parameters' such as air-sea fluxes, the wind field, the initial conditions |
parameters' such as air-sea fluxes, the wind field, the initial conditions |
329 |
etc. Figure \ref{fig:assimilated-globes} shows an estimate of the time-mean |
etc. Figure \ref{fig:assimilated-globes} shows the large scale planetary |
330 |
surface elevation of the ocean obtained by bringing the model in to |
circulation and a Hopf-Muller plot of Equatorial sea-surface height. |
331 |
|
Both are obtained from assimilation bringing the model in to |
332 |
consistency with altimetric and in-situ observations over the period |
consistency with altimetric and in-situ observations over the period |
333 |
1992-1997. {\bf CHANGE THIS TEXT - FIG FROM PATRICK/CARL/DETLEF} |
1992-1997. |
334 |
|
|
335 |
%% CNHbegin |
%% CNHbegin |
336 |
\input{part1/assim_figure} |
\input{part1/assim_figure} |
337 |
%% CNHend |
%% CNHend |
338 |
|
|
339 |
\subsection{Ocean biogeochemical cycles} |
\subsection{Ocean biogeochemical cycles} |
340 |
|
\begin{rawhtml} |
341 |
|
<!-- CMIREDIR:ocean_biogeo_cycles --> |
342 |
|
\end{rawhtml} |
343 |
|
|
344 |
MITgcm is being used to study global biogeochemical cycles in the ocean. For |
MITgcm is being used to study global biogeochemical cycles in the ocean. For |
345 |
example one can study the effects of interannual changes in meteorological |
example one can study the effects of interannual changes in meteorological |
355 |
%%CNHend |
%%CNHend |
356 |
|
|
357 |
\subsection{Simulations of laboratory experiments} |
\subsection{Simulations of laboratory experiments} |
358 |
|
\begin{rawhtml} |
359 |
|
<!-- CMIREDIR:classroom_exp --> |
360 |
|
\end{rawhtml} |
361 |
|
|
362 |
Figure \ref{fig:lab-simulation} shows MITgcm being used to simulate a |
Figure \ref{fig:lab-simulation} shows MITgcm being used to simulate a |
363 |
laboratory experiment inquiring in to the dynamics of the Antarctic Circumpolar Current (ACC). An |
laboratory experiment inquiring into the dynamics of the Antarctic Circumpolar Current (ACC). An |
364 |
initially homogeneous tank of water ($1m$ in diameter) is driven from its |
initially homogeneous tank of water ($1m$ in diameter) is driven from its |
365 |
free surface by a rotating heated disk. The combined action of mechanical |
free surface by a rotating heated disk. The combined action of mechanical |
366 |
and thermal forcing creates a lens of fluid which becomes baroclinically |
and thermal forcing creates a lens of fluid which becomes baroclinically |
376 |
% $Name$ |
% $Name$ |
377 |
|
|
378 |
\section{Continuous equations in `r' coordinates} |
\section{Continuous equations in `r' coordinates} |
379 |
|
\begin{rawhtml} |
380 |
|
<!-- CMIREDIR:z-p_isomorphism --> |
381 |
|
\end{rawhtml} |
382 |
|
|
383 |
To render atmosphere and ocean models from one dynamical core we exploit |
To render atmosphere and ocean models from one dynamical core we exploit |
384 |
`isomorphisms' between equation sets that govern the evolution of the |
`isomorphisms' between equation sets that govern the evolution of the |
387 |
and encoded. The model variables have different interpretations depending on |
and encoded. The model variables have different interpretations depending on |
388 |
whether the atmosphere or ocean is being studied. Thus, for example, the |
whether the atmosphere or ocean is being studied. Thus, for example, the |
389 |
vertical coordinate `$r$' is interpreted as pressure, $p$, if we are |
vertical coordinate `$r$' is interpreted as pressure, $p$, if we are |
390 |
modeling the atmosphere (left hand side of figure \ref{fig:isomorphic-equations}) |
modeling the atmosphere (right hand side of figure \ref{fig:isomorphic-equations}) |
391 |
and height, $z$, if we are modeling the ocean (right hand side of figure |
and height, $z$, if we are modeling the ocean (left hand side of figure |
392 |
\ref{fig:isomorphic-equations}). |
\ref{fig:isomorphic-equations}). |
393 |
|
|
394 |
%%CNHbegin |
%%CNHbegin |
512 |
at fixed and moving $r$ surfaces we set (see figure \ref{fig:zandp-vert-coord}): |
at fixed and moving $r$ surfaces we set (see figure \ref{fig:zandp-vert-coord}): |
513 |
|
|
514 |
\begin{equation} |
\begin{equation} |
515 |
\dot{r}=0atr=R_{fixed}(x,y)\text{ (ocean bottom, top of the atmosphere)} |
\dot{r}=0 \text{\ at\ } r=R_{fixed}(x,y)\text{ (ocean bottom, top of the atmosphere)} |
516 |
\label{eq:fixedbc} |
\label{eq:fixedbc} |
517 |
\end{equation} |
\end{equation} |
518 |
|
|
519 |
\begin{equation} |
\begin{equation} |
520 |
\dot{r}=\frac{Dr}{Dt}atr=R_{moving}\text{ \ |
\dot{r}=\frac{Dr}{Dt} \text{\ at\ } r=R_{moving}\text{ \ |
521 |
(ocean surface,bottom of the atmosphere)} \label{eq:movingbc} |
(ocean surface,bottom of the atmosphere)} \label{eq:movingbc} |
522 |
\end{equation} |
\end{equation} |
523 |
|
|
655 |
|
|
656 |
\subsection{Hydrostatic, Quasi-hydrostatic, Quasi-nonhydrostatic and |
\subsection{Hydrostatic, Quasi-hydrostatic, Quasi-nonhydrostatic and |
657 |
Non-hydrostatic forms} |
Non-hydrostatic forms} |
658 |
|
\begin{rawhtml} |
659 |
|
<!-- CMIREDIR:non_hydrostatic --> |
660 |
|
\end{rawhtml} |
661 |
|
|
662 |
|
|
663 |
Let us separate $\phi $ in to surface, hydrostatic and non-hydrostatic terms: |
Let us separate $\phi $ in to surface, hydrostatic and non-hydrostatic terms: |
664 |
|
|