/[MITgcm]/manual/s_examples/global_oce_latlon/climatalogical_ogcm.tex
ViewVC logotype

Diff of /manual/s_examples/global_oce_latlon/climatalogical_ogcm.tex

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by cnh, Thu Feb 28 19:32:19 2002 UTC revision 1.14 by molod, Tue Jun 27 19:08:22 2006 UTC
# Line 1  Line 1 
1  % $Header$  % $Header$
2  % $Name$  % $Name$
3    
4  \section{Global Ocean Simulation at 4$^\circ$ Resolution}  \section[Global Ocean MITgcm Exmaple]{Global Ocean Simulation at $4^\circ$ Resolution}
5    \label{www:tutorials}
6  \label{sect:eg-global}  \label{sect:eg-global}
7    \begin{rawhtml}
8    <!-- CMIREDIR:eg-global: -->
9    \end{rawhtml}
10    
11  \bodytext{bgcolor="#FFFFFFFF"}  \bodytext{bgcolor="#FFFFFFFF"}
12    
# Line 21  This example experiment demonstrates usi Line 25  This example experiment demonstrates usi
25  the planetary ocean circulation. The simulation is configured  the planetary ocean circulation. The simulation is configured
26  with realistic geography and bathymetry on a  with realistic geography and bathymetry on a
27  $4^{\circ} \times 4^{\circ}$ spherical polar grid.  $4^{\circ} \times 4^{\circ}$ spherical polar grid.
28    The files for this experiment are in the verification directory
29    under tutorial\_global\_oce\_latlon.
30  Twenty levels are used in the vertical, ranging in thickness  Twenty levels are used in the vertical, ranging in thickness
31  from $50\,{\rm m}$ at the surface to $815\,{\rm m}$ at depth,  from $50\,{\rm m}$ at the surface to $815\,{\rm m}$ at depth,
32  giving a maximum model depth of $6\,{\rm km}$.  giving a maximum model depth of $6\,{\rm km}$.
# Line 29  can be integrated forward for thousands Line 35  can be integrated forward for thousands
35  processor desktop computer.  processor desktop computer.
36  \\  \\
37  \subsection{Overview}  \subsection{Overview}
38    \label{www:tutorials}
39    
40  The model is forced with climatological wind stress data and surface  The model is forced with climatological wind stress data and surface
41  flux data from DaSilva \cite{DaSilva94}. Climatological data  flux data from DaSilva \cite{DaSilva94}. Climatological data
# Line 86  simulation are described in the experime Line 93  simulation are described in the experime
93    
94    
95  \subsection{Discrete Numerical Configuration}  \subsection{Discrete Numerical Configuration}
96    \label{www:tutorials}
97    
98    
99   The model is configured in hydrostatic form.  The domain is discretised with   The model is configured in hydrostatic form.  The domain is discretised with
# Line 203  elevation $\eta$ and the hydrostatic pre Line 211  elevation $\eta$ and the hydrostatic pre
211  \\  \\
212    
213  \subsubsection{Numerical Stability Criteria}  \subsubsection{Numerical Stability Criteria}
214    \label{www:tutorials}
215    
216  The Laplacian dissipation coefficient, $A_{h}$, is set to $5 \times 10^5 m s^{-1}$.  The Laplacian dissipation coefficient, $A_{h}$, is set to $5 \times 10^5 m s^{-1}$.
217  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},
218  \begin{eqnarray}  \begin{eqnarray}
219  \label{EQ:eg-global-munk_layer}  \label{EQ:eg-global-munk_layer}
220  M_{w} = \pi ( \frac { A_{h} }{ \beta } )^{\frac{1}{3}}  && M_{w} = \pi ( \frac { A_{h} }{ \beta } )^{\frac{1}{3}}
221  \end{eqnarray}  \end{eqnarray}
222    
223  \noindent  of $\approx 600$km. This is greater than the model  \noindent  of $\approx 600$km. This is greater than the model
# Line 222  $\delta t_{v}=40~{\rm minutes}$ for mome Line 231  $\delta t_{v}=40~{\rm minutes}$ for mome
231  parameter to the horizontal Laplacian friction \cite{adcroft:95}  parameter to the horizontal Laplacian friction \cite{adcroft:95}
232  \begin{eqnarray}  \begin{eqnarray}
233  \label{EQ:eg-global-laplacian_stability}  \label{EQ:eg-global-laplacian_stability}
234  S_{l} = 4 \frac{A_{h} \delta t_{v}}{{\Delta x}^2}  && S_{l} = 4 \frac{A_{h} \delta t_{v}}{{\Delta x}^2}
235  \end{eqnarray}  \end{eqnarray}
236    
237  \noindent evaluates to 0.16 at a latitude of $\phi=80^{\circ}$, which is below the  \noindent evaluates to 0.16 at a latitude of $\phi=80^{\circ}$, which is below the
# Line 299  S_{c} = \frac{c_{g} \delta t_{v}}{ \Delt Line 308  S_{c} = \frac{c_{g} \delta t_{v}}{ \Delt
308  stability limit of 0.5.  stability limit of 0.5.
309        
310  \subsection{Experiment Configuration}  \subsection{Experiment Configuration}
311    \label{www:tutorials}
312  \label{SEC:eg-global-clim_ocn_examp_exp_config}  \label{SEC:eg-global-clim_ocn_examp_exp_config}
313    
314  The model configuration for this experiment resides under the  The model configuration for this experiment resides under the
# Line 325  experiments. Below we describe the custo Line 335  experiments. Below we describe the custo
335  to these files associated with this experiment.  to these files associated with this experiment.
336    
337  \subsubsection{Driving Datasets}  \subsubsection{Driving Datasets}
338    \label{www:tutorials}
339    
340  Figures (\ref{FIG:sim_config_tclim}-\ref{FIG:sim_config_empmr}) show the  Figures (\ref{FIG:sim_config_tclim}-\ref{FIG:sim_config_empmr}) show the
341  relaxation temperature ($\theta^{\ast}$) and salinity ($S^{\ast}$) fields,  relaxation temperature ($\theta^{\ast}$) and salinity ($S^{\ast}$) fields,
# Line 337  domain. Line 348  domain.
348    
349    
350  \subsubsection{File {\it input/data}}  \subsubsection{File {\it input/data}}
351    \label{www:tutorials}
352    
353  This file, reproduced completely below, specifies the main parameters  This file, reproduced completely below, specifies the main parameters
354  for the experiment. The parameters that are significant for this configuration  for the experiment. The parameters that are significant for this configuration
# Line 348  are Line 360  are
360  \begin{verbatim} tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,  \end{verbatim}  \begin{verbatim} tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,  \end{verbatim}
361  $\cdots$ \\  $\cdots$ \\
362  set reference values for potential  set reference values for potential
363  temperature and salinity at each model level in units of $^{\circ}$C and  temperature and salinity at each model level in units of $^{\circ}\mathrm{C}$ and
364  ${\rm ppt}$. The entries are ordered from surface to depth.  ${\rm ppt}$. The entries are ordered from surface to depth.
365  Density is calculated from anomalies at each level evaluated  Density is calculated from anomalies at each level evaluated
366  with respect to the reference values set here.\\  with respect to the reference values set here.\\
# Line 626  notes. Line 638  notes.
638  \end{small}  \end{small}
639    
640  \subsubsection{File {\it input/data.pkg}}  \subsubsection{File {\it input/data.pkg}}
641    \label{www:tutorials}
642    
643  This file uses standard default values and does not contain  This file uses standard default values and does not contain
644  customisations for this experiment.  customisations for this experiment.
645    
646  \subsubsection{File {\it input/eedata}}  \subsubsection{File {\it input/eedata}}
647    \label{www:tutorials}
648    
649  This file uses standard default values and does not contain  This file uses standard default values and does not contain
650  customisations for this experiment.  customisations for this experiment.
651    
652  \subsubsection{File {\it input/windx.sin\_y}}  \subsubsection{File {\it input/windx.sin\_y}}
653    \label{www:tutorials}
654    
655  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$)
656  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}$.
# Line 646  in MITgcm. The included matlab program { Line 661  in MITgcm. The included matlab program {
661  code for creating the {\it input/windx.sin\_y} file.  code for creating the {\it input/windx.sin\_y} file.
662    
663  \subsubsection{File {\it input/topog.box}}  \subsubsection{File {\it input/topog.box}}
664    \label{www:tutorials}
665    
666    
667  The {\it input/topog.box} file specifies a two-dimensional ($x,y$)  The {\it input/topog.box} file specifies a two-dimensional ($x,y$)
# Line 657  The included matlab program {\it input/g Line 673  The included matlab program {\it input/g
673  code for creating the {\it input/topog.box} file.  code for creating the {\it input/topog.box} file.
674    
675  \subsubsection{File {\it code/SIZE.h}}  \subsubsection{File {\it code/SIZE.h}}
676    \label{www:tutorials}
677    
678  Two lines are customized in this file for the current experiment  Two lines are customized in this file for the current experiment
679    
# Line 683  the vertical domain extent in grid point Line 700  the vertical domain extent in grid point
700  \end{small}  \end{small}
701    
702  \subsubsection{File {\it code/CPP\_OPTIONS.h}}  \subsubsection{File {\it code/CPP\_OPTIONS.h}}
703    \label{www:tutorials}
704    
705  This file uses standard default values and does not contain  This file uses standard default values and does not contain
706  customisations for this experiment.  customisations for this experiment.
707    
708    
709  \subsubsection{File {\it code/CPP\_EEOPTIONS.h}}  \subsubsection{File {\it code/CPP\_EEOPTIONS.h}}
710    \label{www:tutorials}
711    
712  This file uses standard default values and does not contain  This file uses standard default values and does not contain
713  customisations for this experiment.  customisations for this experiment.
714    
715  \subsubsection{Other Files }  \subsubsection{Other Files }
716    \label{www:tutorials}
717    
718  Other files relevant to this experiment are  Other files relevant to this experiment are
719  \begin{itemize}  \begin{itemize}
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22