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

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

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


Revision 1.26 - (show annotations) (download) (as text)
Fri Aug 29 17:08:21 2014 UTC (10 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint01, HEAD
Changes since 1.25: +141 -64 lines
File MIME type: application/x-tex
Error occurred while calculating annotation data.
update the list of experiments

1 % $Header: /u/gcmpack/manual/s_examples/text/model_examples.tex,v 1.25 2013/03/26 13:21:28 jmc Exp $
2 % $Name: $
3
4 \section[MITgcm Example Experiments]{Example experiments}
5 \label{sec:modelExamples}
6 \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 \end{itemize}
430
431 Once you have chosen the example you want to run, you are ready to
432 compile the code.
433

  ViewVC Help
Powered by ViewVC 1.1.22