| 1 |
cnh |
1.2 |
% $Header: /u/gcmpack/mitgcmdoc/part1/illustration.tex,v 1.1.1.1 2001/08/08 16:16:16 adcroft Exp $ |
| 2 |
|
|
% $Name: $ |
| 3 |
adcroft |
1.1 |
|
| 4 |
|
|
\section{Illustrations of the model in action} |
| 5 |
|
|
|
| 6 |
|
|
The MITgcm has been designed and used to model a vast range of phenomena, |
| 7 |
|
|
from convection on the scale of meters in the ocean to the global pattern of |
| 8 |
|
|
atmospheric winds - see fig.2\ref{fig:all-scales}. To give a flavor of the |
| 9 |
|
|
kinds of problems the model has been used to study, we briefly describe some |
| 10 |
|
|
of them here. A more detailed description of the underlying formulation, |
| 11 |
|
|
numerical algorithm and implementation that lie behind these calculations is |
| 12 |
|
|
given later. Indeed it is easy to reproduce the results shown here: simply |
| 13 |
|
|
download the model (the minimum you need is a PC running linux, together |
| 14 |
|
|
with a FORTRAN\ 77 compiler) and follow the examples. |
| 15 |
|
|
|
| 16 |
|
|
\subsection{Global atmosphere: `Held-Suarez' benchmark} |
| 17 |
|
|
|
| 18 |
|
|
Fig.E1a.\ref{fig:Held-Suarez} is an instaneous plot of the 500$mb$ height |
| 19 |
|
|
field obtained using a 5-level version of the atmospheric pressure isomorph |
| 20 |
|
|
run at 300$km$ resolution. We see fully developed baroclinic eddies along |
| 21 |
|
|
the northern hemisphere storm track. There are no mountains (but you can |
| 22 |
|
|
easily put them in) or land-sea contrast. The model is driven by relaxation |
| 23 |
|
|
to a radiative-convective equilibrium profile, following the description set |
| 24 |
|
|
out in Held and Suarez; 1994 designed to test atmospheric hydrodynamical |
| 25 |
|
|
cores - there are no mountains or land-sea contrast. As decribed in Adcroft |
| 26 |
|
|
(2001), `cubed sphere' is used to descretize the globe permitting a uniform |
| 27 |
|
|
gridding and obviated the need to fourier filter. |
| 28 |
|
|
|
| 29 |
|
|
Fig.E1b shows the 5-year mean, zonally averaged potential temperature, zonal |
| 30 |
|
|
wind and meridional overturning streamfunction from the 5-level model. |
| 31 |
|
|
|
| 32 |
|
|
A regular spherical lat-lon grid can also be used. |
| 33 |
|
|
|
| 34 |
|
|
Results from this integration, together with a 20-level calculation, can be |
| 35 |
|
|
found here - link through to channel. |
| 36 |
|
|
|
| 37 |
|
|
This calculation takes 1 day of computation per year of intergation on a |
| 38 |
|
|
pentium IV PC. |
| 39 |
|
|
|
| 40 |
|
|
\subsection{Ocean gyres} |
| 41 |
|
|
|
| 42 |
|
|
\subsection{Global ocean circulation} |
| 43 |
|
|
|
| 44 |
|
|
Fig.E2a shows the pattern of ocean currents at the surface of a 4$^{\circ }$ |
| 45 |
|
|
global ocean model run with 15 vertical levels. The model is driven using |
| 46 |
|
|
monthly-mean winds with mixed boundary conditions on temperature and |
| 47 |
|
|
salinity at the surface. Fig.E2b shows the overturning (thermohaline) |
| 48 |
cnh |
1.2 |
circulation. Lopped cells are used to represent topography on a regular $ |
| 49 |
adcroft |
1.1 |
lat-lon$ grid extending from 70$^{\circ }N$ to 70$^{\circ }S$. |
| 50 |
|
|
|
| 51 |
|
|
\subsection{Flow over topography} |
| 52 |
|
|
|
| 53 |
|
|
\subsection{Ocean convection} |
| 54 |
|
|
|
| 55 |
|
|
Fig.E3 shows convection over a slope using the non-hydrostatic ocean |
| 56 |
|
|
isomorph and lopped cells to respresent topography. .....The grid resolution |
| 57 |
|
|
is |
| 58 |
|
|
|
| 59 |
|
|
\subsection{Boundary forced internal waves} |
| 60 |
|
|
|
| 61 |
|
|
\subsection{Carbon outgassing sensitivity} |
| 62 |
|
|
|
| 63 |
|
|
Fig.E4 shows.... |