--- manual/s_examples/rotating_tank/tank.tex 2004/07/26 16:21:15 1.3 +++ manual/s_examples/rotating_tank/tank.tex 2004/07/26 19:13:08 1.6 @@ -1,4 +1,4 @@ -% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/rotating_tank/tank.tex,v 1.3 2004/07/26 16:21:15 afe Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/manual/s_examples/rotating_tank/tank.tex,v 1.6 2004/07/26 19:13:08 afe Exp $ % $Name: $ \bodytext{bgcolor="#FFFFFFFF"} @@ -13,22 +13,28 @@ %{\large May 2001} %\end{center} -This is the first in a series of tutorials describing -example MITgcm numerical experiments. The example experiments -include both straightforward examples of idealized geophysical -fluid simulations and more involved cases encompassing -large scale modeling and -automatic differentiation. Both hydrostatic and non-hydrostatic -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 -configure the MITgcm code and input files for each experiment. - \section{A Rotating Tank in Cylindrical Coordinates} \label{sect:eg-tank} \label{www:tutorials} +This section illustrates an example of MITgcm simulating a laboratory +experiment on much smaller scales than those common to geophysical +fluid dynamics. + +\subsection{Overview} +\label{www:tutorials} + + +This example experiment demonstrates using the MITgcm to simulate +a laboratory experiment with a rotating tank of water with an ice +bucket in the center. The simulation is configured for a laboratory +scale on a +$3^{\circ}$ $\times$ 20cm +cyclindrical grid with twenty-nine vertical +levels. +\\ + + This example experiment demonstrates using the MITgcm to simulate a Barotropic, wind-forced, ocean gyre circulation. The experiment @@ -187,18 +193,19 @@ \label{www:tutorials} \label{SEC:eg-baro-code_config} -The model configuration for this experiment resides under the -directory {\it verification/exp0/}. The experiment files +The model configuration for this experiment resides under the +directory {\it verification/rotatingi\_tank/}. The experiment files \begin{itemize} \item {\it input/data} \item {\it input/data.pkg} \item {\it input/eedata}, -\item {\it input/windx.sin\_y}, -\item {\it input/topog.box}, +\item {\it input/bathyPol.bin}, +\item {\it input/thetaPol.bin}, \item {\it code/CPP\_EEOPTIONS.h} \item {\it code/CPP\_OPTIONS.h}, -\item {\it code/SIZE.h}. +\item {\it code/SIZE.h}. \end{itemize} + contain the code customizations and parameter settings for this experiments. Below we describe the customizations to these files associated with this experiment. @@ -212,17 +219,23 @@ \begin{itemize} -\item Line 7, \begin{verbatim} viscAh=4.E2, \end{verbatim} this line sets -the Laplacian friction coefficient to $400 m^2s^{-1}$ +\item Line X, \begin{verbatim} viscAh=5.0E-6, \end{verbatim} this line sets +the Laplacian friction coefficient to $0.000006 m^2s^{-1}$, which is ususally +low because of the small scale, presumably.... qqq + +\item Line X, \begin{verbatim}f0=0.5 , \end{verbatim} this line sets the +coriolis term, and represents a tank spinning at qqq \item Line 10, \begin{verbatim} beta=1.E-11, \end{verbatim} this line sets $\beta$ (the gradient of the coriolis parameter, $f$) to $10^{-11} s^{-1}m^{-1}$ \item Lines 15 and 16 \begin{verbatim} -rigidLid=.FALSE., -implicitFreeSurface=.TRUE., +rigidLid=.TRUE., +implicitFreeSurface=.FALSE., \end{verbatim} -these lines suppress the rigid lid formulation of the surface + +these lines do the opposite of the following: +suppress the rigid lid formulation of the surface pressure inverter and activate the implicit free surface form of the pressure inverter. @@ -234,23 +247,17 @@ and implicitly suppresses searching for checkpoint files associated with restarting an numerical integration from a previously saved state. -\item Line 29, -\begin{verbatim} -endTime=12000, -\end{verbatim} -this line indicates that the experiment should start finish at $t=12000s$. -A restart file will be written at this time that will enable the -simulation to be continued from this point. - \item Line 30, \begin{verbatim} -deltaTmom=1200, +deltaT=0.1, \end{verbatim} -This line sets the momentum equation timestep to $1200s$. +This line sets the integration timestep to $0.1s$. This is an unsually +small value among the examples due to the small physical scale of the +experiment. \item Line 39, \begin{verbatim} -usingCartesianGrid=.TRUE., +usingCylindricalGrid=.TRUE., \end{verbatim} This line requests that the simulation be performed in a Cartesian coordinate system. @@ -312,7 +319,7 @@ notes. \begin{small} -\input{part3/case_studies/barotropic_gyre/input/data} +\input{part3/case_studies/rotating_tank/input/data} \end{small} \subsubsection{File {\it input/data.pkg}} @@ -327,28 +334,21 @@ This file uses standard default values and does not contain customizations for this experiment. -\subsubsection{File {\it input/windx.sin\_y}} +\subsubsection{File {\it input/thetaPol.bin}} \label{www:tutorials} -The {\it input/windx.sin\_y} file specifies a two-dimensional ($x,y$) -map of wind stress ,$\tau_{x}$, values. The units used are $Nm^{-2}$. -Although $\tau_{x}$ is only a function of $y$n in this experiment -this file must still define a complete two-dimensional map in order -to be compatible with the standard code for loading forcing fields -in MITgcm. The included matlab program {\it input/gendata.m} gives a complete -code for creating the {\it input/windx.sin\_y} file. +The {\it input/thetaPol.bin} file specifies a three-dimensional ($x,y,z$) +map of initial values of $\theta$ in degrees Celsius. -\subsubsection{File {\it input/topog.box}} +\subsubsection{File {\it input/bathyPol.bin}} \label{www:tutorials} -The {\it input/topog.box} file specifies a two-dimensional ($x,y$) +The {\it input/bathyPol.bin} file specifies a two-dimensional ($x,y$) map of depth values. For this experiment values are either -$0m$ or {\bf -delZ}m, corresponding respectively to a wall or to deep -ocean. The file contains a raw binary stream of data that is enumerated +$0m$ or {\bf -delZ}m, corresponding respectively to outside or inside of +the tank. The file contains a raw binary stream of data that is enumerated in the same way as standard MITgcm two-dimensional, horizontal arrays. -The included matlab program {\it input/gendata.m} gives a complete -code for creating the {\it input/topog.box} file. \subsubsection{File {\it code/SIZE.h}} \label{www:tutorials}