/[MITgcm]/manual/s_examples/text/model_examples.tex
ViewVC logotype

Diff of /manual/s_examples/text/model_examples.tex

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

revision 1.6 by edhill, Mon Aug 1 23:53:01 2005 UTC revision 1.26 by jmc, Fri Aug 29 17:08:21 2014 UTC
# Line 1  Line 1 
1  % $Header$  % $Header$
2  % $Name$  % $Name$
3    
4  \section{Tutorials}  \section[MITgcm Example Experiments]{Example experiments}
5  \label{sect:tutorials}  \label{sec:modelExamples}
6  \label{www:tutorials}  \begin{rawhtml}
7    <!-- CMIREDIR:modelExamples: -->
8    \end{rawhtml}
9    
10    %% a set of pre-configured numerical experiments
11    
12    The full MITgcm distribution comes with a set of pre-configured
13    numerical experiments.  Some of these example experiments are tests of
14    individual parts of the model code, but many are fully fledged
15    numerical simulations. Full tutorials exist for a few of the examples,
16    and are documented in sections \ref{sec:eg-baro} -
17    \ref{sec:eg-tank}. The other examples follow the same general
18    structure as the tutorial examples. However, they only include brief
19    instructions in a text file called {\it README}.  The examples are
20    located in subdirectories under the directory \texttt{verification}.
21    Each example is briefly described below.
22    
23    \subsection{Full list of model examples}
24    
25    \begin{enumerate}
26    
27    \item \texttt{tutorial\_advection\_in\_gyre} - Test of various
28      advection schemes in a single-layer double-gyre experiment.
29      This experiment is described in detail in section
30      \ref{sec:eg-adv-gyre}.
31    
32    \item \texttt{tutorial\_baroclinic\_gyre} - Four layer, ocean double
33      gyre. This experiment is described in detail in section
34      \ref{sec:eg-fourlayer}.
35    
36    \item \texttt{tutorial\_barotropic\_gyre} - Single layer, ocean double
37      gyre (barotropic with free-surface).
38      This experiment is described in detail in section \ref{sec:eg-baro}.
39    
40    \item \texttt{tutorial\_cfc\_offline} - Offline form of the MITgcm to
41      study advection of a passive tracer and CFCs.
42      This experiment is described in detail in section \ref{sec:eg-offline-cfc}.
43    
44    \item \texttt{tutorial\_deep\_convection} - Non-uniformly forced
45      ocean convection in a doubly periodic box. This experiment is
46      described in detail in section \ref{sec:eg-bconv}.
47    
48    \item \texttt{tutorial\_dic\_adoffline} - Offline form of MITgcm
49      dynamics coupled to the dissolved inorganic carbon biogeochemistry model;
50      adjoint set-up.
51    
52    \item \texttt{tutorial\_global\_oce\_biogeo} - Ocean model coupled to
53      the dissolved inorganic carbon biogeochemistry model. This
54      experiment is described in detail in section
55      \ref{sec:eg-biogeochem_tutorial}.
56    
57    \item \texttt{tutorial\_global\_oce\_in\_p} - Global ocean simulation in
58      pressure coordinate (non-Boussinesq ocean model). Described in
59      detail in section \ref{sec:eg-globalpressure}.
60    
61    \item \texttt{tutorial\_global\_oce\_latlon} - 4x4 degree global ocean
62      simulation with steady climatological forcing. This experiment is
63      described in detail in section \ref{sec:eg-global}.
64    
65    \item \texttt{tutorial\_global\_oce\_optim} - Global ocean state
66      estimation at $4^\circ$ resolution.  This experiment is described in
67      detail in section \ref{sec:eg-global_state_estimate}.
68    
69    \item \texttt{tutorial\_held\_suarez\_cs} - 3D atmosphere dynamics
70      using Held and Suarez (1994) forcing on cubed sphere grid.  This
71      experiment is described in detail in section \ref{sec:eg-hs}.
72    
73    \item \texttt{tutorial\_offline} - Offline form of the MITgcm to study
74      advection of a passive tracer.  This experiment is described in
75      detail in section \ref{sec:eg-offline}.
76    
77    \item \texttt{tutorial\_plume\_on\_slope} - Gravity Plume on a
78      continental slope.  This experiment is described in detail in
79      section \ref{sec:eg-gravityplume}.
80    
81    \item \texttt{tutorial\_tracer\_adjsens} - Simple passive tracer
82      experiment. Includes derivative calculation. This experiment is
83      described in detail in section \ref{sec:eg-simple-tracer-adjoint}.\\
84      Also contains an additional set-up using Secon Order Moment (SOM) advection
85      scheme ({\it input\_ad.som81/}).
86    
87    \item \texttt{1D\_ocean\_ice\_column} - Oceanic column with seaice on top.
88    
89    \item \texttt{adjustment.128x64x1} - Barotropic adjustment problem on
90      latitude longitude grid with 128x64 grid points ($2.8^\circ$ resolution).
91    
92    \item \texttt{adjustment.cs-32x32x1} - Barotropic adjustment problem on
93      cube sphere grid with 32x32 points per face (roughly $2.8^\circ$
94      resolution).\\
95      Also contains a non-linear free-surface adjustment version ({\it input.nlfs/}).
96    
97    \item \texttt{advect\_cs} - Two-dimensional passive advection test on
98      cube sphere grid (32x32 grid points per face, roughly $2.8^\circ$ resolution)
99    
100    \item \texttt{advect\_xy} - Two-dimensional (horizontal plane) passive
101      advection test on Cartesian grid.\\
102      Also contains an additional set-up using Adams-Bashforth 3 ({\it input.ab3\_c4/}).
103    
104    \item \texttt{advect\_xz} - Two-dimensional (vertical plane) passive
105      advection test on Cartesian grid.\\
106      Also contains an additional set-up using non-linear free-surface
107       with divergent barotropic flow and implicit vertical advection ({\it input.nlfs/}).
108    
109    \item \texttt{aim.5l\_Equatorial\_Channel} -
110      5-levels Intermediate Atmospheric physics,
111      3D Equatorial Channel configuration.
112    
113    \item \texttt{aim.5l\_LatLon} - 5-levels Intermediate Atmospheric physics,
114      Global configuration, on latitude longitude grid with 128x64x5 grid
115      points ($2.8^\circ$ resolution).
116    
117    \item \texttt{aim.5l\_cs} - 5-levels Intermediate Atmospheric physics,
118      Global configuration on cube sphere grid
119      (32x32 grid points per face, roughly $2.8^\circ$).\\
120      Also contains an additional set-up with a slab-ocean and thermodynamic
121      sea-ice ({\it input.thSI/}).
122    
123    \item \texttt{bottom\_ctrl\_5x5} - Adjoint test using the bottom
124      topography as the control parameter.
125    
126    \item \texttt{cfc\_example} - Global ocean with online computation and
127      advection of CFC11 and CFC12.
128    
129    \item \texttt{cheapAML\_box} - Example using cheap atmospheric mixed layer
130       (cheapAML) package.
131    
132    \item \texttt{cpl\_aim+ocn} - Coupled Ocean - Atmosphere realistic
133      configuration on cubed-sphere cs32 horizontal grid,
134      using Intermediate Atmospheric physics ({\it pkg/aim\_v23})
135      thermodynamic seaice ({\it pkg/thsice}) and land packages.
136      on cubed-sphere cs32 in a realistic configuration.
137    
138    \item \texttt{cpl\_atm2d+ocn} - Coupled Ocean - Atmosphere realistic
139      configuration using 2-D Atmospheric Model ({\it pkg/atm2d}).
140    
141    \item \texttt{deep\_anelastic} - Convection simulation on a giant planet:
142      relax both the Boussinesq approximation (anelastic) and the thin atmosphere
143      approximation (deep atmosphere).
144    
145    \item \texttt{dome} - Idealized 3D test of a density-driven bottom current.
146    
147    \item \texttt{exp2} - Old version of the global ocean experiment (no GM,
148          no partial-cells).\\
149      Also contains an additional set-up with rigid-lid ({\it input.rigidLid/}).
150    
151    \item \texttt{exp4} - Flow over a Gaussian bump in open-water or
152      channel with open boundaries.\\
153      Also contains an additional set-up using non-linear free-surface ({\it input.nlfs/}).
154    
155    \item \texttt{fizhi-cs-32x32x40} - Global atmospheric simulation with
156      realistic topography, 40 vertical levels, a cubed sphere grid and
157      the full atmospheric physics package.
158    
159    \item \texttt{fizhi-cs-aqualev20} - Global atmospheric simulation on an
160      aqua planet with full atmospheric physics. Run is perpetual march
161      with an analytical SST distribution.  This is the configuration for
162      the APE (Aqua Planet Experiment) participation experiment.
163    
164    \item \texttt{fizhi-gridalt-hs} - Global atmospheric simulation
165      Held-Suarez (1994) forcing, with the physical forcing and the
166      dynamical forcing running on different vertical grids.
167    
168    \item \texttt{flt\_example} - Example of using float package.
169    
170    \item \texttt{front\_relax} - Relaxation of an ocean thermal front
171      (test for Gent/McWilliams scheme). 2D (y-z).\\
172      Also contains additional set-ups:
173      \begin{enumerate}
174       \item using the Boundary-Value Problem method
175              (Ferrari et al., 2010) ({\it input.bvp/}).
176       \item with Mixed-Layer Eddy parameterization
177              (Ferrari \& McWilliams, 2007) ({\it input.mxl/}).
178      \end{enumerate}
179    
180    \item \texttt{global\_ocean.90x40x15} - Global ocean simulation at 4x4
181      degree resolution. Similar to tutorial\_global\_oce\_latlon, but using
182      $z^*$ coordinates with quasi-non-hydrostatic and non-hydrostatic metric terms.
183      This experiment also illustrate the use of SBO package.
184      Also contains additional set-ups:
185      \begin{enumerate}
186       \item using down-slope package ({\it pkg/down\_slope}) ({\it input.dwnslp/})
187       \item an Open-AD adjoint set-up ({\it code\_oad/, input\_oad/}).
188       \item four TAF adjoint set-ups ({\it code\_ad/}):
189       \begin{enumerate}
190         \item standard experiment ({\it input\_ad/}).
191         \item with bottom drag as a control ({\it input\_ad.bottomdrag/}).
192         \item with kappa GM as a control ({\it input\_ad.kapgm/}).
193         \item with kappa Redi as a control ({\it input\_ad.kapredi/}).
194       \end{enumerate}
195      \end{enumerate}
196    
197    \item \texttt{global\_ocean.cs32x15} - Global ocean experiment on the
198      cubed sphere grid.\\
199      Also contains additional forward set-ups:
200      \begin{enumerate}
201       \item non-hydrostatic with biharmonic viscosity ({\it input.viscA4/})
202       \item using thermodynamic sea ice and bulk force ({\it input.thsice/})
203       \item using thermodynamic ({\it pkg/thsice}) dynamic ({\it pkg/seaice}) sea-ice
204             and {\it exf} package ({\it input.icedyn/})
205       \item using thermodynamic - dynamic ({\it pkg/seaice}) sea-ice
206             with {\it exf} package ({\it input.seaice/})
207      \end{enumerate}
208      and few additional adjoint set-ups ({\it code\_ad/}):
209      \begin{enumerate}
210       \item standard experiment without sea-ice ({\it input\_ad/}).
211       \item using thermodynamic - dynamic sea-ice ({\it input\_ad.seaice/})
212       \item same as above without adjoint sea-ice dynamics ({\it input\_ad.seaice\_dynmix/})
213       \item using thermodynamic sea-ice from {\it thsice} package ({\it input\_ad.thsice/})
214      \end{enumerate}
215    
216    \item \texttt{global\_ocean\_ebm} - Global ocean experiment on a lat-lon
217      grid coupled to an atmospheric energy balance model. Similar to
218      global\_ocean.90x40x15 experiment.\\
219      Also contains an adjoint set-up ({\it code\_ad/, input\_ad/}).
220    
221    \item \texttt{global\_with\_exf} - Global ocean experiment on a lat-lon
222      grid using the {\it exf} package. Similar to tutorial\_global\_oce\_latlon
223      experiment.\\
224      Also contains a secondary set-up with yearly {\it exf} fields ({\it input\_ad.yearly/}).
225    
226    \item \texttt{halfpipe\_streamice} - Example using package "streamice".\\
227      Also contains adjoint set-ups using TAF ({\it code\_ad/, input\_ad/})
228       and using Open-AD ({\it code\_oad/, input\_oad/}).
229    
230    \item \texttt{hs94.128x64x5} - 3D atmosphere dynamics on lat-lon grid,
231      using Held and Suarez '94 forcing.
232    
233    \item \texttt{hs94.1x64x5} - Zonal averaged atmosphere dynamics
234      using Held and Suarez '94 forcing.\\
235      Also contains adjoint set-ups using TAF ({\it code\_ad/, input\_ad/})
236       and using Open-AD ({\it code\_oad/, input\_oad/}).
237    
238    \item \texttt{hs94.cs-32x32x5} - 3D atmosphere dynamics using Held and
239      Suarez (1994) forcing on the cubed sphere, similar to tutorial\_held\_suarez\_cs
240      experiment but using linear free-surface and only 5 levels.\\
241      Also contains an additional set-up with Implicit Internal gravity waves
242      treatment and Adams-Bashforth 3 ({\it input.impIGW/}).
243    
244    \item \texttt{ideal\_2D\_oce} - Idealized 2D global ocean simulation on
245      an aqua planet.
246    
247    \item \texttt{internal\_wave} - Ocean internal wave forced by open
248      boundary conditions.\\
249      Also contains an additional set-up using {\it pkg/kl10} (see section
250       \ref{sec:pkg:kl10}, Klymak and Legg, 2010) ({\it input.kl10/}).
251    
252    \item \texttt{inverted\_barometer} - Simple test of ocean response to
253      atmospheric pressure loading.
254    
255    \item \texttt{isomip} - ISOMIP like set-up including ice-shelf cavities
256      ({\it pkg/shelfice}).\\
257      Also contains additional set-ups:
258      \begin{enumerate}
259       \item with "htd" ({\it input.htd/})
260              but only Martin knows what "htd" stands for.
261       \item using package {\it icefront} ({\it input.icefront})
262      \end{enumerate}
263      and also adjoint set-ups using TAF ({\it code\_ad/, input\_ad/, input\_ad.htd/})
264      or using Open-AD ({\it code\_oad/, input\_oad/}).
265    
266    \item \texttt{lab\_sea} - Regional Labrador Sea simulation on a lat-lon
267      grid using the sea ice package.\\
268      Also contains additional set-ups:
269      \begin{enumerate}
270       \item using the simple "free-drift" assumption for seaice ({\it input.fd/})
271       \item using EVP dynamics (instead of LSR solver) and Hibler \& Bryan (1987)
272             sea-ice ocean stress ({\it input.hb87/})
273       \item using package {\it salt\_plume} ({\it input.salt\_plume/})
274      \end{enumerate}
275      and also 3 adjoint set-ups ({\it code\_ad/, input\_ad/, input\_ad.noseaicedyn/,
276       input\_ad.noseaice/}).
277    
278    \item \texttt{matrix\_example} - Test of experimental method to
279      accelerated convergence towards equilibrium.
280    
281    \item \texttt{MLAdjust} - Simple tests for different viscosity formulations.\\
282      Also contains additional set-ups (see: {\it verification/MLAdjust/README}):
283      \begin{enumerate}
284       \item ({\it input.A4FlxF/})
285       \item ({\it input.AhFlxF/})
286       \item ({\it input.AhVrDv/})
287       \item ({\it input.AhStTn/})
288      \end{enumerate}
289    
290    \item \texttt{natl\_box} - Eastern subtropical North Atlantic with KPP
291      scheme; 1 month integration
292    
293    \item \texttt{obcs\_ctrl} - Adjoint test using Open-Boundary conditions
294      as control parameters.
295    
296    \item \texttt{offline\_exf\_seaice} - Seaice on top of oceanic surface layer
297      in an idealized channel. Forcing is computed by bulk-formulae ({\it pkg/exf})
298      with temperature relaxation to prescribed SST (offline ocean).\\
299      Also contains additional set-ups:
300      \begin{enumerate}
301       \item sea-ice dynamics-only using JFNK solver
302             and {\it pkg/thsice} advection ({\it input.dyn\_jfnk/})
303       \item sea-ice dynamics-only using LSR solver
304             and {\it pkg/seaice} advection ({\it input.dyn\_lsr/})
305       \item sea-ice thermodynamics-only using {\it pkg/seaice} ({\it input.thermo/})
306       \item sea-ice thermodynamics-only using {\it pkg/thsice} ({\it input.thsice/})
307      \end{enumerate}
308      and also 2 adjoint set-ups ({\it code\_ad/, input\_ad/, input\_ad.thsice/}).
309    
310    \item \texttt{OpenAD} - Simple Adjoint experiment (used also to test
311      Open-AD compiler)
312    
313    \item \texttt{rotating\_tank} - Rotating tank simulation in cylindrical
314      coordinates.  This experiment is described in detail in section
315      \ref{sec:eg-tank}.
316    
317    \item \texttt{seaice\_itd} - Seaice example using Ice Thickness Distribution (ITD).\\
318      Also contains additional set-ups:
319      \begin{enumerate}
320       \item ({\it input.thermo/})
321       \item ({\it input.lipscomb07/})
322      \end{enumerate}
323    
324    \item \texttt{seaice\_obcs} - Similar to "lab\_sea" ({\it input.salt\_plume/})
325      experiment with only a fraction of the domain and open-boundary conditions
326      derived from "lab\_sea" experiment.\\
327      Also contains additional set-ups:
328      \begin{enumerate}
329       \item ({\it input.seaiceSponge/})
330       \item ({\it input.tides/})
331      \end{enumerate}
332    
333    \item \texttt{short\_surf\_wave} - Short surface wave adjusment
334       (non-hydrostatic) in homogeneous 2-D vertical section (x-z).
335    
336    \item \texttt{so\_box\_biogeo} - Open-boundary Southern ocean box around
337     Drake passage, using same model parameters and forcing as experiment
338      "tutorial\_global\_oce\_biogeo" from which initial conditions
339     and OB conditions have been extracted.
340    
341    \item \texttt{solid-body.cs-32x32x1} - Solid body rotation test for cube
342      sphere grid.
343    
344    \item \texttt{tidal\_basin\_2d} - 2-D vertical section (x-z) with tidal forcing
345     (untested)
346    
347    \item \texttt{vermix} - Simple test in a small domain (3 columns) for
348      ocean vertical mixing schemes. The standard set-up ({\it input/}) uses
349      KPP scheme \cite[]{lar-eta:94}.\\
350      Also contains additional set-ups:
351      \begin{enumerate}
352       \item with Double Diffusion scheme from KPP ({\it input.dd/})
353       \item with \cite{gas-eta:90} ({\it pkg/ggl90}) scheme ({\it input.ggl90/})
354       \item with \cite{Mellor:Yamada1982} level 2. ({\it pkg/my82}) scheme ({\it input.my82/})
355       \item with \cite{pal-rom:97} ({\it pkg/opps}) scheme ({\it input.opps/})
356       \item with \cite{Pacanowski:Philander1981} ({\it pkg/pp81}) scheme ({\it input.pp81/})
357      \end{enumerate}
358    
359    \end{enumerate}
360    
361    \subsection{Directory structure of model examples}
362    
363    Each example directory has the following subdirectories:
364    
365    \begin{itemize}
366    \item \texttt{code}: contains the code particular to the example. At a
367      minimum, this directory includes the following files:
368    
369      \begin{itemize}
370      \item \texttt{code/packages.conf}: declares the list of packages or
371        package groups to be used.  If not included, the default version
372        is located in \texttt{pkg/pkg\_default}.  Package groups are
373        simply convenient collections of commonly used packages which are
374        defined in \texttt{pkg/pkg\_default}.  Some packages may require
375        other packages or may require their absence (that is, they are
376        incompatible) and these package dependencies are listed in
377        \texttt{pkg/pkg\_depend}.
378    
379      \item \texttt{code/CPP\_EEOPTIONS.h}: declares CPP keys relative to
380        the ``execution environment'' part of the code. The default
381        version is located in \texttt{eesupp/inc}.
382    
383      \item \texttt{code/CPP\_OPTIONS.h}: declares CPP keys relative to
384        the ``numerical model'' part of the code. The default version is
385        located in \texttt{model/inc}.
386    
387      \item \texttt{code/SIZE.h}: declares size of underlying
388        computational grid.  The default version is located in
389        \texttt{model/inc}.
390      \end{itemize}
391    
392      In addition, other include files and subroutines might be present in
393      \texttt{code} depending on the particular experiment. See Section 2
394      for more details.
395    
396    \item \texttt{input}: contains the input data files required to run
397      the example. At a minimum, the \texttt{input} directory contains the
398      following files:
399    
400      \begin{itemize}
401      \item \texttt{input/data}: this file, written as a namelist,
402        specifies the main parameters for the experiment.
403    
404      \item \texttt{input/data.pkg}: contains parameters relative to the
405        packages used in the experiment.
406    
407      \item \texttt{input/eedata}: this file contains ``execution
408        environment'' data. At present, this consists of a specification
409        of the number of threads to use in $X$ and $Y$ under multi-threaded
410        execution.
411      \end{itemize}
412    
413      In addition, you will also find in this directory the forcing and
414      topography files as well as the files describing the initial state
415      of the experiment.  This varies from experiment to experiment. See
416      the verification directories refered to in this chapter for more details.
417    
418    \item \texttt{results}: this directory contains the output file
419      \texttt{output.txt} produced by the simulation example. This file is
420      useful for comparison with your own output when you run the
421      experiment.
422    
423    \item \texttt{build}: this directory is initially empty and is used
424      to compile and load the model, and to generate the executable.
425    
426    \item \texttt{run}: this directory is initially empty and is used
427      to run the executable.
428    
429  \input{part3/case_studies/barotropic_gyre/baro.tex}  \end{itemize}
430    
431  \newpage  Once you have chosen the example you want to run, you are ready to
432  \input{part3/case_studies/fourlayer_gyre/fourlayer.tex}  compile the code.
433    
 \newpage  
 \input{part3/case_studies/climatalogical_ogcm/climatalogical_ogcm.tex}  
   
 \newpage  
 \input{part3/case_studies/ogcm_in_pressure/ogcm_in_pressure.tex}  
   
 \newpage  
 \input{part3/case_studies/hs_atmosphere/hs_atmos.tex}  
   
 \newpage  
 \input{part3/case_studies/doubly_periodic_convection/convection.tex}  
   
 \newpage  
 \input{part3/case_studies/plume_on_slope/plume_on_slope.tex}  
   
 \newpage  
 \input{part3/case_studies/carbon_outgassing_sensitivity/co2sens.tex}  
   
 \newpage  
 \input{part3/case_studies/biogeochem_tutorial/biogeochem.tex}  
   
 \begin{versionprivate}  
 \newpage  
 \input{part3/case_studies/rotating_tank/tank.tex}  
 \end{versionprivate}  

Legend:
Removed from v.1.6  
changed lines
  Added in v.1.26

  ViewVC Help
Powered by ViewVC 1.1.22