--- manual/s_examples/barotropic_gyre/baro.tex 2001/09/27 00:58:17 1.2 +++ manual/s_examples/barotropic_gyre/baro.tex 2001/10/22 11:55:47 1.3 @@ -1,7 +1,8 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/barotropic_gyre/baro.tex,v 1.2 2001/09/27 00:58:17 cnh Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/barotropic_gyre/baro.tex,v 1.3 2001/10/22 11:55:47 cnh Exp $ % $Name: $ \section{Example: Barotropic Ocean Gyre In Cartesian Coordinates} +\label{sec:eg-baro} \bodytext{bgcolor="#FFFFFFFF"} @@ -15,15 +16,13 @@ %{\large May 2001} %\end{center} -\subsection{Introduction} - -This document is the first in a series of documents describing +This is the first in a series of sections describing example MITgcm numerical experiments. The example experiments include both straightforward examples of idealised geophysical fluid simulations and more involved cases encompassing large scale modeling and automatic differentiation. Both hydrostatic and non-hydrostatic -experiements are presented, as well as experiments employing +experiments are presented, as well as experiments employing cartesian, spherical-polar and cube-sphere coordinate systems. These ``case study'' documents include information describing the experimental configuration and detailed information on how to @@ -81,33 +80,27 @@ \label{FIG:simulation_config} \end{figure} -\subsection{Discrete Numerical Configuration} - - The model is configured in hydrostatic form. The domain is discretised with -a uniform grid spacing in the horizontal set to - $\Delta x=\Delta y=20$~km, so -that there are sixty grid cells in the $x$ and $y$ directions. Vertically the -model is configured with a single layer with depth, $\Delta z$, of $5000$~m. -The implicit free surface form of the -pressure equation described in Marshall et. al \cite{Marshall97a} is -employed. +\subsection{Equations Solved} +The model is configured in hydrostatic form. The implicit free surface form of the +pressure equation described in Marshall et. al \cite{Marshall97a} is +employed. A horizontal laplacian operator $\nabla_{h}^2$ provides viscous dissipation. The wind-stress momentum input is added to the momentum equation for the ``zonal flow'', $u$. Other terms in the model are explicitly switched off for this experiement configuration (see section -\ref{SEC:code_config} ), yielding an active set of equations solved in this -configuration as follows +\ref{SEC:code_config} ), yielding an active set of equations solved in this +configuration as follows \begin{eqnarray} \label{EQ:model_equations} -\frac{Du}{Dt} - fv + - g\frac{\partial \eta}{\partial x} - - A_{h}\nabla_{h}^2u +\frac{Du}{Dt} - fv + + g\frac{\partial \eta}{\partial x} - + A_{h}\nabla_{h}^2u & = & \frac{\tau_{x}}{\rho_{0}\Delta z} \\ \frac{Dv}{Dt} + fu + g\frac{\partial \eta}{\partial y} - - A_{h}\nabla_{h}^2v + A_{h}\nabla_{h}^2v & = & 0 \\ @@ -117,9 +110,18 @@ \end{eqnarray} \noindent where $u$ and $v$ and the $x$ and $y$ components of the -flow vector $\vec{u}$. +flow vector $\vec{u}$. \\ + +\subsection{Discrete Numerical Configuration} + + The domain is discretised with +a uniform grid spacing in the horizontal set to + $\Delta x=\Delta y=20$~km, so +that there are sixty grid cells in the $x$ and $y$ directions. Vertically the +model is configured with a single layer with depth, $\Delta z$, of $5000$~m. + \subsubsection{Numerical Stability Criteria} The laplacian dissipation coefficient, $A_{h}$, is set to $400 m s^{-1}$.