102 |
numerical algorithm and implementation that lie behind these calculations is |
numerical algorithm and implementation that lie behind these calculations is |
103 |
given later. Indeed many of the illustrative examples shown below can be |
given later. Indeed many of the illustrative examples shown below can be |
104 |
easily reproduced: simply download the model (the minimum you need is a PC |
easily reproduced: simply download the model (the minimum you need is a PC |
105 |
running linux, together with a FORTRAN\ 77 compiler) and follow the examples |
running Linux, together with a FORTRAN\ 77 compiler) and follow the examples |
106 |
described in detail in the documentation. |
described in detail in the documentation. |
107 |
|
|
108 |
\subsection{Global atmosphere: `Held-Suarez' benchmark} |
\subsection{Global atmosphere: `Held-Suarez' benchmark} |
126 |
%% CNHend |
%% CNHend |
127 |
|
|
128 |
As described in Adcroft (2001), a `cubed sphere' is used to discretize the |
As described in Adcroft (2001), a `cubed sphere' is used to discretize the |
129 |
globe permitting a uniform gridding and obviated the need to Fourier filter. |
globe permitting a uniform griding and obviated the need to Fourier filter. |
130 |
The `vector-invariant' form of MITgcm supports any orthogonal curvilinear |
The `vector-invariant' form of MITgcm supports any orthogonal curvilinear |
131 |
grid, of which the cubed sphere is just one of many choices. |
grid, of which the cubed sphere is just one of many choices. |
132 |
|
|
163 |
visible. |
visible. |
164 |
|
|
165 |
%% CNHbegin |
%% CNHbegin |
166 |
\input{part1/ocean_gyres_figure} |
\input{part1/atl6_figure} |
167 |
%% CNHend |
%% CNHend |
168 |
|
|
169 |
|
|
191 |
ocean may be influenced by rotation when the deformation radius is smaller |
ocean may be influenced by rotation when the deformation radius is smaller |
192 |
than the width of the cooling region. Rather than gravity plumes, the |
than the width of the cooling region. Rather than gravity plumes, the |
193 |
mechanism for moving dense fluid down the shelf is then through geostrophic |
mechanism for moving dense fluid down the shelf is then through geostrophic |
194 |
eddies. The simulation shown in the figure \ref{fig::convect-and-topo} |
eddies. The simulation shown in the figure \ref{fig:convect-and-topo} |
195 |
(blue is cold dense fluid, red is |
(blue is cold dense fluid, red is |
196 |
warmer, lighter fluid) employs the non-hydrostatic capability of MITgcm to |
warmer, lighter fluid) employs the non-hydrostatic capability of MITgcm to |
197 |
trigger convection by surface cooling. The cold, dense water falls down the |
trigger convection by surface cooling. The cold, dense water falls down the |
231 |
|
|
232 |
As one example of application of the MITgcm adjoint, Figure \ref{fig:hf-sensitivity} |
As one example of application of the MITgcm adjoint, Figure \ref{fig:hf-sensitivity} |
233 |
maps the gradient $\frac{\partial J}{\partial \mathcal{H}}$where $J$ is the magnitude |
maps the gradient $\frac{\partial J}{\partial \mathcal{H}}$where $J$ is the magnitude |
234 |
of the overturning streamfunction shown in figure \ref{fig:large-scale-circ} |
of the overturning stream-function shown in figure \ref{fig:large-scale-circ} |
235 |
at 60$^{\circ }$N and $ |
at 60$^{\circ }$N and $ |
236 |
\mathcal{H}(\lambda,\varphi)$ is the mean, local air-sea heat flux over |
\mathcal{H}(\lambda,\varphi)$ is the mean, local air-sea heat flux over |
237 |
a 100 year period. We see that $J$ is |
a 100 year period. We see that $J$ is |
248 |
An important application of MITgcm is in state estimation of the global |
An important application of MITgcm is in state estimation of the global |
249 |
ocean circulation. An appropriately defined `cost function', which measures |
ocean circulation. An appropriately defined `cost function', which measures |
250 |
the departure of the model from observations (both remotely sensed and |
the departure of the model from observations (both remotely sensed and |
251 |
insitu) over an interval of time, is minimized by adjusting `control |
in-situ) over an interval of time, is minimized by adjusting `control |
252 |
parameters' such as air-sea fluxes, the wind field, the initial conditions |
parameters' such as air-sea fluxes, the wind field, the initial conditions |
253 |
etc. Figure \ref{fig:assimilated-globes} shows an estimate of the time-mean |
etc. Figure \ref{fig:assimilated-globes} shows an estimate of the time-mean |
254 |
surface elevation of the ocean obtained by bringing the model in to |
surface elevation of the ocean obtained by bringing the model in to |
277 |
\subsection{Simulations of laboratory experiments} |
\subsection{Simulations of laboratory experiments} |
278 |
|
|
279 |
Figure \ref{fig:lab-simulation} shows MITgcm being used to simulate a |
Figure \ref{fig:lab-simulation} shows MITgcm being used to simulate a |
280 |
laboratory experiment enquiring in to the dynamics of the Antarctic Circumpolar Current (ACC). An |
laboratory experiment inquiring in to the dynamics of the Antarctic Circumpolar Current (ACC). An |
281 |
initially homogeneous tank of water ($1m$ in diameter) is driven from its |
initially homogeneous tank of water ($1m$ in diameter) is driven from its |
282 |
free surface by a rotating heated disk. The combined action of mechanical |
free surface by a rotating heated disk. The combined action of mechanical |
283 |
and thermal forcing creates a lens of fluid which becomes baroclinically |
and thermal forcing creates a lens of fluid which becomes baroclinically |
351 |
|
|
352 |
\begin{equation} |
\begin{equation} |
353 |
\frac{DS}{Dt}=\mathcal{Q}_{S}\text{ humidity/salinity} |
\frac{DS}{Dt}=\mathcal{Q}_{S}\text{ humidity/salinity} |
354 |
\label{eq:humidtity_salt} |
\label{eq:humidity_salt} |
355 |
\end{equation} |
\end{equation} |
356 |
|
|
357 |
Here: |
Here: |
432 |
|
|
433 |
\begin{equation} |
\begin{equation} |
434 |
\dot{r}=\frac{Dr}{Dt}atr=R_{moving}\text{ \ |
\dot{r}=\frac{Dr}{Dt}atr=R_{moving}\text{ \ |
435 |
(oceansurface,bottomoftheatmosphere)} \label{eq:movingbc} |
(ocean surface,bottom of the atmosphere)} \label{eq:movingbc} |
436 |
\end{equation} |
\end{equation} |
437 |
|
|
438 |
Here |
Here |
525 |
atmosphere)} \label{eq:moving-bc-atmos} |
atmosphere)} \label{eq:moving-bc-atmos} |
526 |
\end{eqnarray} |
\end{eqnarray} |
527 |
|
|
528 |
Then the (hydrostatic form of) equations (\ref{eq:horizontal_mtm}-\ref{eq:humidity_slainty}) |
Then the (hydrostatic form of) equations (\ref{eq:horizontal_mtm}-\ref{eq:humidity_salt}) |
529 |
yields a consistent set of atmospheric equations which, for convenience, are written out in $p$ |
yields a consistent set of atmospheric equations which, for convenience, are written out in $p$ |
530 |
coordinates in Appendix Atmosphere - see eqs(\ref{eq:atmos-prime}). |
coordinates in Appendix Atmosphere - see eqs(\ref{eq:atmos-prime}). |
531 |
|
|
562 |
\end{eqnarray} |
\end{eqnarray} |
563 |
where $\eta $ is the elevation of the free surface. |
where $\eta $ is the elevation of the free surface. |
564 |
|
|
565 |
Then equations (\ref{eq:horizontal_mtm}-\ref{eq:humidity_slainty}) yield a consistent set |
Then equations (\ref{eq:horizontal_mtm}-\ref{eq:humidity_salt}) yield a consistent set |
566 |
of oceanic equations |
of oceanic equations |
567 |
which, for convenience, are written out in $z$ coordinates in Appendix Ocean |
which, for convenience, are written out in $z$ coordinates in Appendix Ocean |
568 |
- see eqs(\ref{eq:ocean-mom}) to (\ref{eq:ocean-salt}). |
- see eqs(\ref{eq:ocean-mom}) to (\ref{eq:ocean-salt}). |
1178 |
\label{eq:non-boussinesq} |
\label{eq:non-boussinesq} |
1179 |
\end{eqnarray} |
\end{eqnarray} |
1180 |
These equations permit acoustics modes, inertia-gravity waves, |
These equations permit acoustics modes, inertia-gravity waves, |
1181 |
non-hydrostatic motions, a geostrophic (Rossby) mode and a thermo-haline |
non-hydrostatic motions, a geostrophic (Rossby) mode and a thermohaline |
1182 |
mode. As written, they cannot be integrated forward consistently - if we |
mode. As written, they cannot be integrated forward consistently - if we |
1183 |
step $\rho $ forward in (\ref{eq-zns-cont}), the answer will not be |
step $\rho $ forward in (\ref{eq-zns-cont}), the answer will not be |
1184 |
consistent with that obtained by stepping (\ref{eq-zns-heat}) and (\ref |
consistent with that obtained by stepping (\ref{eq-zns-heat}) and (\ref |