--- manual/s_examples/baroclinic_gyre/fourlayer.tex 2002/05/16 15:54:37 1.13 +++ manual/s_examples/baroclinic_gyre/fourlayer.tex 2003/07/30 13:42:53 1.14 @@ -1,4 +1,4 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/baroclinic_gyre/fourlayer.tex,v 1.13 2002/05/16 15:54:37 adcroft Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/baroclinic_gyre/fourlayer.tex,v 1.14 2003/07/30 13:42:53 edhill Exp $ % $Name: $ \section{Four Layer Baroclinic Ocean Gyre In Spherical Coordinates} @@ -383,13 +383,13 @@ $x$ and $y$. The values specified here are read into the variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} tRef \begin{rawhtml} \end{rawhtml} } in the model code, by procedure {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. @@ -402,7 +402,7 @@ %% \file{ini_parms} \newcommand{\VARtref}{ {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} tRef \begin{rawhtml} \end{rawhtml} } @@ -417,7 +417,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -430,19 +430,19 @@ for this operator are specified later. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} viscAz \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and is copied into model general vertical coordinate variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} viscAr \begin{rawhtml} \end{rawhtml} }. At each time step, the viscous term contribution to the momentum equations @@ -455,7 +455,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -469,13 +469,13 @@ for this operator are specified later. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} viscAh \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and applied in routines {\it CALC\_MOM\_RHS} and {\it CALC\_GW}. @@ -486,7 +486,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -497,7 +497,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -512,13 +512,13 @@ $\frac{\partial v}{\partial x}$=0 along boundaries in $x$. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} no\_slip\_sides \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and the boundary condition is evaluated in routine @@ -531,7 +531,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -545,13 +545,13 @@ e.g. $u=v=0$ at $z=-H$, where $H$ is the local depth of the domain. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} no\_slip\_bottom \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and is applied in the routine {\it S/R CALC\_MOM\_RHS}. @@ -562,7 +562,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -577,13 +577,13 @@ all boundaries. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} diffKhT \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and used in routine {\it S/R CALC\_GT}. @@ -593,7 +593,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -607,19 +607,19 @@ operator is $\frac{\partial}{\partial z}$ = 0 on all boundaries. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} diffKzT \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. It is copied into model general vertical coordinate variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} diffKrT \begin{rawhtml} \end{rawhtml} } which is used in routine {\it S/R CALC\_DIFFUSIVITY}. @@ -629,7 +629,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -644,13 +644,13 @@ to $2 \times 10^{-4}\,{\rm degrees}^{-1}$ The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} tAlpha \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. The routine {\it S/R FIND\_RHO} makes use of {\bf tAlpha}. @@ -661,7 +661,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -673,13 +673,13 @@ This line selects the linear form of the equation of state. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} eosType \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. The values of {\bf eosType} sets which formula in routine @@ -691,7 +691,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -709,13 +709,13 @@ on spherical polar geometry. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} usingSphericalPolarGrid \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. When set to {\bf .TRUE.} the settings of {\bf delX} and {\bf delY} are @@ -728,7 +728,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -746,13 +746,13 @@ not alter the kernel equation discretisation. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} phiMin \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and is used in routine {\it INI\_SPEHRICAL\_POLAR\_GRID}. @@ -763,7 +763,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -776,13 +776,13 @@ in the discrete grid to $1^{\circ}$ in longitude. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} delX \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and is used in routine {\it INI\_SPEHRICAL\_POLAR\_GRID}. @@ -793,7 +793,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -806,13 +806,13 @@ in the discrete grid to $1^{\circ}$ in latitude. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} delY \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} } and is used in routine {\it INI\_SPEHRICAL\_POLAR\_GRID}. @@ -823,7 +823,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -837,20 +837,20 @@ is $2\,{\rm km}$. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} delZ \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. It is copied into the internal model coordinate variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} delR \begin{rawhtml} \end{rawhtml} } which is used in routine {\it INI\_VERTICAL\_GRID}. @@ -861,7 +861,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -883,13 +883,13 @@ bathymetry file. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} bathyFile \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. The bathymetry file is read in the routine {\it INI\_DEPTHS}. @@ -900,7 +900,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} } @@ -919,13 +919,13 @@ file. The variable {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} zonalWindFile \begin{rawhtml} \end{rawhtml} } is read in the routine {\it -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} INI\_PARMS \begin{rawhtml} \end{rawhtml} }. The wind-stress file is read in the routine @@ -937,7 +937,7 @@ \end{minipage} } {\bf -\begin{rawhtml} \end{rawhtml} +\begin{rawhtml} \end{rawhtml} goto code \begin{rawhtml} \end{rawhtml} }