| 25 |
the experimental configuration and detailed information on how to |
the experimental configuration and detailed information on how to |
| 26 |
configure the MITgcm code and input files for each experiment. |
configure the MITgcm code and input files for each experiment. |
| 27 |
|
|
| 28 |
\section{Barotropic Ocean Gyre In Cartesian Coordinates} |
\section[Barotropic Gyre MITgcm Example]{Barotropic Ocean Gyre In Cartesian Coordinates} |
| 29 |
\label{sect:eg-baro} |
\label{sect:eg-baro} |
| 30 |
|
\label{www:tutorials} |
| 31 |
|
\begin{rawhtml} |
| 32 |
|
<!-- CMIREDIR:eg-baro: --> |
| 33 |
|
\end{rawhtml} |
| 34 |
|
|
| 35 |
|
|
| 36 |
This example experiment demonstrates using the MITgcm to simulate |
This example experiment demonstrates using the MITgcm to simulate |
| 71 |
Figure \ref{FIG:eg-baro-simulation_config} |
Figure \ref{FIG:eg-baro-simulation_config} |
| 72 |
summarizes the configuration simulated. |
summarizes the configuration simulated. |
| 73 |
|
|
| 74 |
|
%% === eh3 === |
| 75 |
\begin{figure} |
\begin{figure} |
| 76 |
\begin{center} |
%% \begin{center} |
| 77 |
\resizebox{7.5in}{5.5in}{ |
%% \resizebox{7.5in}{5.5in}{ |
| 78 |
\includegraphics*[0.2in,0.7in][10.5in,10.5in] |
%% \includegraphics*[0.2in,0.7in][10.5in,10.5in] |
| 79 |
{part3/case_studies/barotropic_gyre/simulation_config.eps} } |
%% {part3/case_studies/barotropic_gyre/simulation_config.eps} } |
| 80 |
\end{center} |
%% \end{center} |
| 81 |
|
\centerline{ |
| 82 |
|
\scalefig{.95} |
| 83 |
|
\epsfbox{part3/case_studies/barotropic_gyre/simulation_config.eps} |
| 84 |
|
} |
| 85 |
\caption{Schematic of simulation domain and wind-stress forcing function |
\caption{Schematic of simulation domain and wind-stress forcing function |
| 86 |
for barotropic gyre numerical experiment. The domain is enclosed bu solid |
for barotropic gyre numerical experiment. The domain is enclosed bu solid |
| 87 |
walls at $x=$~0,1200km and at $y=$~0,1200km.} |
walls at $x=$~0,1200km and at $y=$~0,1200km.} |
| 89 |
\end{figure} |
\end{figure} |
| 90 |
|
|
| 91 |
\subsection{Equations Solved} |
\subsection{Equations Solved} |
| 92 |
|
\label{www:tutorials} |
| 93 |
The model is configured in hydrostatic form. The implicit free surface form of the |
The model is configured in hydrostatic form. The implicit free surface form of the |
| 94 |
pressure equation described in Marshall et. al \cite{marshall:97a} is |
pressure equation described in Marshall et. al \cite{marshall:97a} is |
| 95 |
employed. |
employed. |
| 124 |
|
|
| 125 |
|
|
| 126 |
\subsection{Discrete Numerical Configuration} |
\subsection{Discrete Numerical Configuration} |
| 127 |
|
\label{www:tutorials} |
| 128 |
|
|
| 129 |
The domain is discretised with |
The domain is discretised with |
| 130 |
a uniform grid spacing in the horizontal set to |
a uniform grid spacing in the horizontal set to |
| 133 |
model is configured with a single layer with depth, $\Delta z$, of $5000$~m. |
model is configured with a single layer with depth, $\Delta z$, of $5000$~m. |
| 134 |
|
|
| 135 |
\subsubsection{Numerical Stability Criteria} |
\subsubsection{Numerical Stability Criteria} |
| 136 |
|
\label{www:tutorials} |
| 137 |
|
|
| 138 |
The Laplacian dissipation coefficient, $A_{h}$, is set to $400 m s^{-1}$. |
The Laplacian dissipation coefficient, $A_{h}$, is set to $400 m s^{-1}$. |
| 139 |
This value is chosen to yield a Munk layer width \cite{adcroft:95}, |
This value is chosen to yield a Munk layer width \cite{adcroft:95}, |
| 187 |
of 0.5 and limits $\delta t$ to $1200s$. |
of 0.5 and limits $\delta t$ to $1200s$. |
| 188 |
|
|
| 189 |
\subsection{Code Configuration} |
\subsection{Code Configuration} |
| 190 |
|
\label{www:tutorials} |
| 191 |
\label{SEC:eg-baro-code_config} |
\label{SEC:eg-baro-code_config} |
| 192 |
|
|
| 193 |
The model configuration for this experiment resides under the |
The model configuration for this experiment resides under the |
| 207 |
to these files associated with this experiment. |
to these files associated with this experiment. |
| 208 |
|
|
| 209 |
\subsubsection{File {\it input/data}} |
\subsubsection{File {\it input/data}} |
| 210 |
|
\label{www:tutorials} |
| 211 |
|
|
| 212 |
This file, reproduced completely below, specifies the main parameters |
This file, reproduced completely below, specifies the main parameters |
| 213 |
for the experiment. The parameters that are significant for this configuration |
for the experiment. The parameters that are significant for this configuration |
| 319 |
\end{small} |
\end{small} |
| 320 |
|
|
| 321 |
\subsubsection{File {\it input/data.pkg}} |
\subsubsection{File {\it input/data.pkg}} |
| 322 |
|
\label{www:tutorials} |
| 323 |
|
|
| 324 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
| 325 |
customizations for this experiment. |
customizations for this experiment. |
| 326 |
|
|
| 327 |
\subsubsection{File {\it input/eedata}} |
\subsubsection{File {\it input/eedata}} |
| 328 |
|
\label{www:tutorials} |
| 329 |
|
|
| 330 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
| 331 |
customizations for this experiment. |
customizations for this experiment. |
| 332 |
|
|
| 333 |
\subsubsection{File {\it input/windx.sin\_y}} |
\subsubsection{File {\it input/windx.sin\_y}} |
| 334 |
|
\label{www:tutorials} |
| 335 |
|
|
| 336 |
The {\it input/windx.sin\_y} file specifies a two-dimensional ($x,y$) |
The {\it input/windx.sin\_y} file specifies a two-dimensional ($x,y$) |
| 337 |
map of wind stress ,$\tau_{x}$, values. The units used are $Nm^{-2}$. |
map of wind stress ,$\tau_{x}$, values. The units used are $Nm^{-2}$. |
| 342 |
code for creating the {\it input/windx.sin\_y} file. |
code for creating the {\it input/windx.sin\_y} file. |
| 343 |
|
|
| 344 |
\subsubsection{File {\it input/topog.box}} |
\subsubsection{File {\it input/topog.box}} |
| 345 |
|
\label{www:tutorials} |
| 346 |
|
|
| 347 |
|
|
| 348 |
The {\it input/topog.box} file specifies a two-dimensional ($x,y$) |
The {\it input/topog.box} file specifies a two-dimensional ($x,y$) |
| 354 |
code for creating the {\it input/topog.box} file. |
code for creating the {\it input/topog.box} file. |
| 355 |
|
|
| 356 |
\subsubsection{File {\it code/SIZE.h}} |
\subsubsection{File {\it code/SIZE.h}} |
| 357 |
|
\label{www:tutorials} |
| 358 |
|
|
| 359 |
Two lines are customized in this file for the current experiment |
Two lines are customized in this file for the current experiment |
| 360 |
|
|
| 377 |
\end{small} |
\end{small} |
| 378 |
|
|
| 379 |
\subsubsection{File {\it code/CPP\_OPTIONS.h}} |
\subsubsection{File {\it code/CPP\_OPTIONS.h}} |
| 380 |
|
\label{www:tutorials} |
| 381 |
|
|
| 382 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
| 383 |
customizations for this experiment. |
customizations for this experiment. |
| 384 |
|
|
| 385 |
|
|
| 386 |
\subsubsection{File {\it code/CPP\_EEOPTIONS.h}} |
\subsubsection{File {\it code/CPP\_EEOPTIONS.h}} |
| 387 |
|
\label{www:tutorials} |
| 388 |
|
|
| 389 |
This file uses standard default values and does not contain |
This file uses standard default values and does not contain |
| 390 |
customizations for this experiment. |
customizations for this experiment. |