| 1 | 
adcroft | 
1.1 | 
% $Header: $ | 
| 2 | 
  | 
  | 
% $Name: $ | 
| 3 | 
  | 
  | 
 | 
| 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 | 
  | 
  | 
circulation. Lopped cells are used to represent topography on a regular $% | 
| 49 | 
  | 
  | 
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.... |