/[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.18 by jmc, Fri Aug 27 13:25:32 2010 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 Example]{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    \begin{center}
11    (in directory: {\it verification/tutorial\_global\_oce\_latlon/})
12    \end{center}
13    
14  \bodytext{bgcolor="#FFFFFFFF"}  \bodytext{bgcolor="#FFFFFFFF"}
15    
# Line 21  This example experiment demonstrates usi Line 28  This example experiment demonstrates usi
28  the planetary ocean circulation. The simulation is configured  the planetary ocean circulation. The simulation is configured
29  with realistic geography and bathymetry on a  with realistic geography and bathymetry on a
30  $4^{\circ} \times 4^{\circ}$ spherical polar grid.  $4^{\circ} \times 4^{\circ}$ spherical polar grid.
31    The files for this experiment are in the verification directory
32    under tutorial\_global\_oce\_latlon.
33  Twenty levels are used in the vertical, ranging in thickness  Twenty levels are used in the vertical, ranging in thickness
34  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,
35  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 38  can be integrated forward for thousands
38  processor desktop computer.  processor desktop computer.
39  \\  \\
40  \subsection{Overview}  \subsection{Overview}
41    \label{www:tutorials}
42    
43  The model is forced with climatological wind stress data and surface  The model is forced with climatological wind stress data and surface
44  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 96  simulation are described in the experime
96    
97    
98  \subsection{Discrete Numerical Configuration}  \subsection{Discrete Numerical Configuration}
99    \label{www:tutorials}
100    
101    
102   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 214  elevation $\eta$ and the hydrostatic pre
214  \\  \\
215    
216  \subsubsection{Numerical Stability Criteria}  \subsubsection{Numerical Stability Criteria}
217    \label{www:tutorials}
218    
219  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}$.
220  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},
221  \begin{eqnarray}  \begin{eqnarray}
222  \label{EQ:eg-global-munk_layer}  \label{EQ:eg-global-munk_layer}
223  M_{w} = \pi ( \frac { A_{h} }{ \beta } )^{\frac{1}{3}}  && M_{w} = \pi ( \frac { A_{h} }{ \beta } )^{\frac{1}{3}}
224  \end{eqnarray}  \end{eqnarray}
225    
226  \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 234  $\delta t_{v}=40~{\rm minutes}$ for mome
234  parameter to the horizontal Laplacian friction \cite{adcroft:95}  parameter to the horizontal Laplacian friction \cite{adcroft:95}
235  \begin{eqnarray}  \begin{eqnarray}
236  \label{EQ:eg-global-laplacian_stability}  \label{EQ:eg-global-laplacian_stability}
237  S_{l} = 4 \frac{A_{h} \delta t_{v}}{{\Delta x}^2}  && S_{l} = 4 \frac{A_{h} \delta t_{v}}{{\Delta x}^2}
238  \end{eqnarray}  \end{eqnarray}
239    
240  \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 311  S_{c} = \frac{c_{g} \delta t_{v}}{ \Delt
311  stability limit of 0.5.  stability limit of 0.5.
312        
313  \subsection{Experiment Configuration}  \subsection{Experiment Configuration}
314    \label{www:tutorials}
315  \label{SEC:eg-global-clim_ocn_examp_exp_config}  \label{SEC:eg-global-clim_ocn_examp_exp_config}
316    
317  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 338  experiments. Below we describe the custo
338  to these files associated with this experiment.  to these files associated with this experiment.
339    
340  \subsubsection{Driving Datasets}  \subsubsection{Driving Datasets}
341    \label{www:tutorials}
342    
343  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
344  relaxation temperature ($\theta^{\ast}$) and salinity ($S^{\ast}$) fields,  relaxation temperature ($\theta^{\ast}$) and salinity ($S^{\ast}$) fields,
# Line 337  domain. Line 351  domain.
351    
352    
353  \subsubsection{File {\it input/data}}  \subsubsection{File {\it input/data}}
354    \label{www:tutorials}
355    
356  This file, reproduced completely below, specifies the main parameters  This file, reproduced completely below, specifies the main parameters
357  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 363  are
363  \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}
364  $\cdots$ \\  $\cdots$ \\
365  set reference values for potential  set reference values for potential
366  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
367  ${\rm ppt}$. The entries are ordered from surface to depth.  ${\rm ppt}$. The entries are ordered from surface to depth.
368  Density is calculated from anomalies at each level evaluated  Density is calculated from anomalies at each level evaluated
369  with respect to the reference values set here.\\  with respect to the reference values set here.\\
# Line 568  See section \ref{SEC:cd_scheme}. Line 583  See section \ref{SEC:cd_scheme}.
583  \fbox{  \fbox{
584  \begin{minipage}{5.0in}  \begin{minipage}{5.0in}
585  {\it S/R INI\_PARMS}({\it ini\_parms.F})\\  {\it S/R INI\_PARMS}({\it ini\_parms.F})\\
586  {\it S/R CALC\_MOM\_RHS}({\it calc\_mom\_rhs.F})  {\it S/R MOM\_FLUXFORM}({\it mom\_fluxform.F})
587  \end{minipage}  \end{minipage}
588  }  }
589    
# Line 622  that are described in the MITgcm Getting Line 637  that are described in the MITgcm Getting
637  notes.  notes.
638    
639  \begin{small}  \begin{small}
640  \input{part3/case_studies/climatalogical_ogcm/input/data}  \input{s_examples/global_oce_latlon/input/data}
641  \end{small}  \end{small}
642    
643  \subsubsection{File {\it input/data.pkg}}  \subsubsection{File {\it input/data.pkg}}
644    \label{www:tutorials}
645    
646  This file uses standard default values and does not contain  This file uses standard default values and does not contain
647  customisations for this experiment.  customisations for this experiment.
648    
649  \subsubsection{File {\it input/eedata}}  \subsubsection{File {\it input/eedata}}
650    \label{www:tutorials}
651    
652  This file uses standard default values and does not contain  This file uses standard default values and does not contain
653  customisations for this experiment.  customisations for this experiment.
654    
655  \subsubsection{File {\it input/windx.sin\_y}}  \subsubsection{File {\it input/windx.sin\_y}}
656    \label{www:tutorials}
657    
658  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$)
659  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 664  in MITgcm. The included matlab program {
664  code for creating the {\it input/windx.sin\_y} file.  code for creating the {\it input/windx.sin\_y} file.
665    
666  \subsubsection{File {\it input/topog.box}}  \subsubsection{File {\it input/topog.box}}
667    \label{www:tutorials}
668    
669    
670  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 676  The included matlab program {\it input/g
676  code for creating the {\it input/topog.box} file.  code for creating the {\it input/topog.box} file.
677    
678  \subsubsection{File {\it code/SIZE.h}}  \subsubsection{File {\it code/SIZE.h}}
679    \label{www:tutorials}
680    
681  Two lines are customized in this file for the current experiment  Two lines are customized in this file for the current experiment
682    
# Line 679  the vertical domain extent in grid point Line 699  the vertical domain extent in grid point
699  \end{itemize}  \end{itemize}
700    
701  \begin{small}  \begin{small}
702  \input{part3/case_studies/climatalogical_ogcm/code/SIZE.h}  \input{s_examples/global_oce_latlon/code/SIZE.h}
703  \end{small}  \end{small}
704    
705  \subsubsection{File {\it code/CPP\_OPTIONS.h}}  \subsubsection{File {\it code/CPP\_OPTIONS.h}}
706    \label{www:tutorials}
707    
708  This file uses standard default values and does not contain  This file uses standard default values and does not contain
709  customisations for this experiment.  customisations for this experiment.
710    
711    
712  \subsubsection{File {\it code/CPP\_EEOPTIONS.h}}  \subsubsection{File {\it code/CPP\_EEOPTIONS.h}}
713    \label{www:tutorials}
714    
715  This file uses standard default values and does not contain  This file uses standard default values and does not contain
716  customisations for this experiment.  customisations for this experiment.
717    
718  \subsubsection{Other Files }  \subsubsection{Other Files }
719    \label{www:tutorials}
720    
721  Other files relevant to this experiment are  Other files relevant to this experiment are
722  \begin{itemize}  \begin{itemize}

Legend:
Removed from v.1.8  
changed lines
  Added in v.1.18

  ViewVC Help
Powered by ViewVC 1.1.22