1 | ************************ |
---|
2 | Reference configurations |
---|
3 | ************************ |
---|
4 | |
---|
5 | NEMO is distributed with a set of reference configurations allowing both |
---|
6 | the user to set up his own first applications and |
---|
7 | the developer to test/validate his NEMO developments (using SETTE package). |
---|
8 | |
---|
9 | .. attention:: |
---|
10 | |
---|
11 | Concerning the configurations, |
---|
12 | the NEMO System Team is only in charge of the so-called reference configurations described below. |
---|
13 | |
---|
14 | .. hint:: |
---|
15 | |
---|
16 | Configurations developed by external research projects or initiatives that |
---|
17 | make use of NEMO are welcome to be publicized through the website by |
---|
18 | filling up the form :website:`to add an associated project<projects/add>`. |
---|
19 | |
---|
20 | How to compile an experiment from a reference configuration |
---|
21 | =========================================================== |
---|
22 | |
---|
23 | A user who wants to compile the ORCA2_ICE_PISCES_ reference configuration using ``makenemo`` |
---|
24 | should use the following, by selecting among available architecture file or providing a user defined one: |
---|
25 | |
---|
26 | .. code-block:: console |
---|
27 | |
---|
28 | $ ./makenemo -r 'ORCA2_ICE_PISCES' -m 'my-fortran.fcm' -j '4' |
---|
29 | |
---|
30 | A new ``EXP00`` folder will be created within the selected reference configurations, |
---|
31 | namely ``./cfgs/ORCA2_ICE_PISCES/EXP00``, |
---|
32 | where it will be necessary to uncompress the Input & Forcing Files listed in the above table. |
---|
33 | |
---|
34 | Then it will be possible to launch the execution of the model through a runscript |
---|
35 | (opportunely adapted to the user system). |
---|
36 | |
---|
37 | List of Configurations |
---|
38 | ====================== |
---|
39 | |
---|
40 | All forcing files listed below in the table are available from |NEMO archives URL|_ |
---|
41 | |
---|
42 | .. |NEMO archives URL| image:: https://www.zenodo.org/badge/DOI/10.5281/zenodo.1472245.svg |
---|
43 | .. _NEMO archives URL: https://doi.org/10.5281/zenodo.1472245 |
---|
44 | |
---|
45 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
46 | Configuration Component(s) Input & Forcing File(s) |
---|
47 | ---------------------- ---------------------------------- ------------------------------------------------ |
---|
48 | Name OPA SI3 TOP PISCES AGRIF |
---|
49 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
50 | AGRIF_DEMO_ X X X AGRIF_DEMO_v4.0.tar, ORCA2_ICE_v4.0.tar |
---|
51 | AMM12_ X AMM12_v4.0.tar |
---|
52 | C1D_PAPA_ X INPUTS_C1D_PAPA_v4.0.tar |
---|
53 | GYRE_BFM_ X X *none* |
---|
54 | GYRE_PISCES_ X X X *none* |
---|
55 | ORCA2_ICE_PISCES_ X X X X ORCA2_ICE_v4.0.tar, INPUTS_PISCES_v4.0.tar |
---|
56 | ORCA2_OFF_PISCES_ X X ORCA2_OFF_v4.0.tar, INPUTS_PISCES_v4.0.tar |
---|
57 | ORCA2_OFF_TRC_ X ORCA2_OFF_v4.0.tar |
---|
58 | ORCA2_SAS_ICE_ X ORCA2_ICE_v4.0.tar, INPUTS_SAS_v4.0.tar |
---|
59 | SPITZ12_ X X SPITZ12_v4.0.tar |
---|
60 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
61 | |
---|
62 | AGRIF_DEMO |
---|
63 | ---------- |
---|
64 | |
---|
65 | ``AGRIF_DEMO`` is based on the ``ORCA2_ICE_PISCES`` global configuration at 2° of resolution with |
---|
66 | the inclusion of 3 online nested grids to demonstrate the overall capabilities of AGRIF in |
---|
67 | a realistic context (including the nesting of sea ice models). |
---|
68 | |
---|
69 | The configuration includes a 1:1 grid in the Pacific and two successively nested grids with |
---|
70 | odd and even refinement ratios over the Arctic ocean, |
---|
71 | with the finest grid spanning the whole Svalbard archipelago that is of |
---|
72 | particular interest to test sea ice coupling. |
---|
73 | |
---|
74 | The 1:1 grid can be used alone as a benchmark to check that |
---|
75 | the model solution is not corrupted by grid exchanges. |
---|
76 | Note that since grids interact only at the baroclinic time level, |
---|
77 | numerically exact results can not be achieved in the 1:1 case. |
---|
78 | Perfect reproducibility is obtained only by switching to a fully explicit setup instead of a split explicit free surface scheme. |
---|
79 | |
---|
80 | AMM12 |
---|
81 | ----- |
---|
82 | |
---|
83 | ``AMM12`` stands for *Atlantic Margin Model at 12 km* that is |
---|
84 | a regional configuration covering the Northwest European Shelf domain on |
---|
85 | a regular horizontal grid of ~12 km of resolution (see :cite:`ODEA2012`). |
---|
86 | |
---|
87 | This configuration allows to tests several features of NEMO specifically addressed to the shelf seas. |
---|
88 | In particular, ``AMM12`` accounts for vertical s-coordinates system, GLS turbulence scheme, |
---|
89 | tidal lateral boundary conditions using a flather scheme (see more in ``BDY``). |
---|
90 | Boundaries may be completely omitted by setting ``ln_bdy = .false.`` in ``nambdy``. |
---|
91 | |
---|
92 | Sample surface fluxes, river forcing and an initial restart file are included to test a realistic model run |
---|
93 | (``AMM12_v4.0.tar``). |
---|
94 | |
---|
95 | Note that, the Baltic boundary is included within the river input file and is specified as a river source, |
---|
96 | but unlike ordinary river points the Baltic inputs also include salinity and temperature data. |
---|
97 | |
---|
98 | C1D_PAPA |
---|
99 | -------- |
---|
100 | |
---|
101 | ``C1D_PAPA`` is a 1D configuration for the `PAPA station <http://www.pmel.noaa.gov/OCS/Papa/index-Papa.shtml>`_ located in the northern-eastern Pacific Ocean at 50.1°N, 144.9°W. |
---|
102 | See `Reffray et al. (2015) <http://www.geosci-model-dev.net/8/69/2015>`_ for the description of its physical and numerical turbulent-mixing behaviour. |
---|
103 | |
---|
104 | The water column setup, called NEMO1D, is activated with the inclusion of the CPP key ``key_c1d`` and |
---|
105 | has a horizontal domain of 3x3 grid points. |
---|
106 | |
---|
107 | This reference configuration uses 75 vertical levels grid (1m at the surface), GLS turbulence scheme with K-epsilon closure and the NCAR bulk formulae. |
---|
108 | Data provided with ``INPUTS_C1D_PAPA_v4.0.tar`` file account for: |
---|
109 | |
---|
110 | - ``forcing_PAPASTATION_1h_y201[0-1].nc`` : ECMWF operational analysis atmospheric forcing rescaled to 1h (with long and short waves flux correction) for years 2010 and 2011 |
---|
111 | - ``init_PAPASTATION_m06d15.nc`` : Initial Conditions from observed data and Levitus 2009 climatology |
---|
112 | - ``chlorophyll_PAPASTATION.nc`` : surface chlorophyll file from Seawifs data |
---|
113 | |
---|
114 | GYRE_BFM |
---|
115 | -------- |
---|
116 | |
---|
117 | ``GYRE_BFM`` shares the same physical setup of GYRE_PISCES_, but NEMO is coupled with the `BFM <http://www.bfm-community.eu/>`_ biogeochemical model as described in ``./cfgs/GYRE_BFM/README``. |
---|
118 | |
---|
119 | GYRE_PISCES |
---|
120 | ----------- |
---|
121 | |
---|
122 | ``GYRE_PISCES`` is an idealized configuration representing a Northern hemisphere double gyres system, |
---|
123 | in the Beta-plane approximation with a regular 1° horizontal resolution and 31 vertical levels, |
---|
124 | with PISCES BGC model :cite:`gmd-8-2465-2015`. |
---|
125 | Analytical forcing for heat, freshwater and wind-stress fields are applied. |
---|
126 | |
---|
127 | This configuration acts also as demonstrator of the **user defined setup** (``ln_read_cfg = .false.``) and |
---|
128 | grid setting are handled through the ``&namusr_def`` controls in ``namelist_cfg``: |
---|
129 | |
---|
130 | .. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_cfg |
---|
131 | :language: fortran |
---|
132 | :lines: 34-42 |
---|
133 | |
---|
134 | Note that, the default grid size is 30x20 grid points (with ``nn_GYRE = 1``) and |
---|
135 | vertical levels are set by ``jpkglo``. |
---|
136 | The specific code changes can be inspected in ``./src/OCE/USR``. |
---|
137 | |
---|
138 | **Running GYRE as a benchmark** : |
---|
139 | this simple configuration can be used as a benchmark since it is easy to increase resolution, |
---|
140 | with the drawback of getting results that have a very limited physical meaning. |
---|
141 | |
---|
142 | GYRE grid resolution can be increased at runtime by setting a different value of ``nn_GYRE`` (integer multiplier scaling factor), as described in the following table: |
---|
143 | |
---|
144 | =========== ========= ========== ============ =================== |
---|
145 | ``nn_GYRE`` *jpiglo* *jpjglo* ``jpkglo`` **Equivalent to** |
---|
146 | =========== ========= ========== ============ =================== |
---|
147 | 1 30 20 31 GYRE 1° |
---|
148 | 25 750 500 101 ORCA 1/2° |
---|
149 | 50 1500 1000 101 ORCA 1/4° |
---|
150 | 150 4500 3000 101 ORCA 1/12° |
---|
151 | 200 6000 4000 101 ORCA 1/16° |
---|
152 | =========== ========= ========== ============ =================== |
---|
153 | |
---|
154 | Note that, it is necessary to set ``ln_bench = .true.`` in ``namusr_def`` to |
---|
155 | avoid problems in the physics computation and that |
---|
156 | the model timestep should be adequately rescaled. |
---|
157 | |
---|
158 | For example if ``nn_GYRE = 150``, equivalent to an ORCA 1/12° grid, |
---|
159 | the timestep ``rn_rdt = 1200`` should be set to 1200 seconds |
---|
160 | |
---|
161 | Differently from previous versions of NEMO, |
---|
162 | the code uses by default the time-splitting scheme and |
---|
163 | internally computes the number of sub-steps. |
---|
164 | |
---|
165 | ORCA2_ICE_PISCES |
---|
166 | ---------------- |
---|
167 | |
---|
168 | ``ORCA2_ICE_PISCES`` is a reference configuration for the global ocean with |
---|
169 | a 2°x2° curvilinear horizontal mesh and 31 vertical levels, |
---|
170 | distributed using z-coordinate system and with 10 levels in the top 100m. |
---|
171 | ORCA is the generic name given to global ocean Mercator mesh, |
---|
172 | (i.e. variation of meridian scale factor as cosinus of the latitude), |
---|
173 | with two poles in the northern hemisphere so that |
---|
174 | the ratio of anisotropy is nearly one everywhere |
---|
175 | |
---|
176 | this configuration uses the three components |
---|
177 | |
---|
178 | - |OPA|, the ocean dynamical core |
---|
179 | - |SI3|, the thermodynamic-dynamic sea ice model. |
---|
180 | - |TOP|, passive tracer transport module and PISCES BGC model :cite:`gmd-8-2465-2015` |
---|
181 | |
---|
182 | All components share the same grid. |
---|
183 | |
---|
184 | The model is forced with CORE-II normal year atmospheric forcing and |
---|
185 | it uses the NCAR bulk formulae. |
---|
186 | |
---|
187 | In this ``ORCA2_ICE_PISCES`` configuration, |
---|
188 | AGRIF nesting can be activated that includes a nested grid in the Agulhas region. |
---|
189 | |
---|
190 | To set up this configuration, after extracting NEMO: |
---|
191 | |
---|
192 | Build your AGRIF configuration directory from ``ORCA2_ICE_PISCES``, |
---|
193 | with the ``key_agrif`` CPP key activated: |
---|
194 | |
---|
195 | .. code-block:: console |
---|
196 | |
---|
197 | $ ./makenemo -r 'ORCA2_ICE_PISCES' -n 'AGRIF' add_key 'key_agrif' |
---|
198 | |
---|
199 | By using the input files and namelists for ``ORCA2_ICE_PISCES``, |
---|
200 | the AGRIF test configuration is ready to run. |
---|
201 | |
---|
202 | **Ocean Physics** |
---|
203 | |
---|
204 | - *horizontal diffusion on momentum*: the eddy viscosity coefficient depends on the geographical position. It is taken as 40000 m^2/s, reduced in the equator regions (2000 m^2/s) excepted near the western boundaries. |
---|
205 | - *isopycnal diffusion on tracers*: the diffusion acts along the isopycnal surfaces (neutral surface) with an eddy diffusivity coefficient of 2000 m^2/s. |
---|
206 | - *Eddy induced velocity parametrization* with a coefficient that depends on the growth rate of baroclinic instabilities (it usually varies from 15 m^2/s to 3000 m^2/s). |
---|
207 | - *lateral boundary conditions* : zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries. |
---|
208 | - *bottom boundary condition* : zero fluxes of heat and salt are applied through the ocean bottom. |
---|
209 | The Beckmann [19XX] simple bottom boundary layer parameterization is applied along continental slopes. |
---|
210 | A linear friction is applied on momentum. |
---|
211 | - *convection*: the vertical eddy viscosity and diffusivity coefficients are increased to 1 m^2/s in case of static instability. |
---|
212 | - *time step* is 5760sec (1h36') so that there is 15 time steps in one day. |
---|
213 | |
---|
214 | ORCA2_OFF_PISCES |
---|
215 | ---------------- |
---|
216 | |
---|
217 | ``ORCA2_OFF_PISCES`` shares the same general offline configuration of ``ORCA2_ICE_TRC``, |
---|
218 | but only PISCES model is an active component of TOP. |
---|
219 | |
---|
220 | |
---|
221 | ORCA2_OFF_TRC |
---|
222 | ------------- |
---|
223 | |
---|
224 | ``ORCA2_OFF_TRC`` is based on the ORCA2 global ocean configuration |
---|
225 | (see ORCA2_ICE_PISCES_ for general description) along with the tracer passive transport module (TOP), but dynamical fields are pre-calculated and read with specific time frequency. |
---|
226 | |
---|
227 | This enables for an offline coupling of TOP components, |
---|
228 | here specifically inorganic carbon compounds (cfc11, cfc12, sf6, c14) and water age module (age). |
---|
229 | See ``namelist_top_cfg`` to inspect the selection of each component with the dedicated logical keys. |
---|
230 | |
---|
231 | Pre-calculated dynamical fields are provided to NEMO using |
---|
232 | the namelist ``&namdta_dyn`` in ``namelist_cfg``, |
---|
233 | in this case with a 5 days frequency (120 hours): |
---|
234 | |
---|
235 | .. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_ref |
---|
236 | :language: fortran |
---|
237 | :lines: 935-960 |
---|
238 | |
---|
239 | Input dynamical fields for this configuration (``ORCA2_OFF_v4.0.tar``) comes from |
---|
240 | a 2000 years long climatological simulation of ORCA2_ICE using ERA40 atmospheric forcing. |
---|
241 | |
---|
242 | Note that, this configuration default uses linear free surface (``ln_linssh = .true.``) assuming that |
---|
243 | model mesh is not varying in time and |
---|
244 | it includes the bottom boundary layer parameterization (``ln_trabbl = .true.``) that |
---|
245 | requires the provision of bbl coefficients through ``sn_ubl`` and ``sn_vbl`` fields. |
---|
246 | |
---|
247 | It is also possible to activate PISCES model (see ``ORCA2_OFF_PISCES``) or |
---|
248 | a user defined set of tracers and source-sink terms with ``ln_my_trc = .true.`` |
---|
249 | (and adaptation of ``./src/TOP/MY_TRC`` routines). |
---|
250 | |
---|
251 | In addition, the offline module (OFF) allows for the provision of further fields: |
---|
252 | |
---|
253 | 1. **River runoff** can be provided to TOP components by setting ``ln_dynrnf = .true.`` and |
---|
254 | by including an input datastream similarly to the following: |
---|
255 | |
---|
256 | .. code-block:: fortran |
---|
257 | |
---|
258 | sn_rnf = 'dyna_grid_T', 120, 'sorunoff' , .true., .true., 'yearly', '', '', '' |
---|
259 | |
---|
260 | 2. **VVL dynamical fields**, |
---|
261 | in the case input data were produced by a dyamical core using variable volume (``ln_linssh = .false.``) |
---|
262 | it necessary to provide also diverce and E-P at before timestep by |
---|
263 | including input datastreams similarly to the following |
---|
264 | |
---|
265 | .. code-block:: fortran |
---|
266 | |
---|
267 | sn_div = 'dyna_grid_T', 120, 'e3t' , .true., .true., 'yearly', '', '', '' |
---|
268 | sn_empb = 'dyna_grid_T', 120, 'sowaflupb', .true., .true., 'yearly', '', '', '' |
---|
269 | |
---|
270 | |
---|
271 | More details can be found by inspecting the offline data manager in |
---|
272 | the routine ``./src/OFF/dtadyn.F90``. |
---|
273 | |
---|
274 | ORCA2_SAS_ICE |
---|
275 | ------------- |
---|
276 | |
---|
277 | ORCA2_SAS_ICE is a demonstrator of the Stand-Alone Surface (SAS) module and |
---|
278 | it relies on ORCA2 global ocean configuration (see ORCA2_ICE_PISCES_ for general description). |
---|
279 | |
---|
280 | The standalone surface module allows surface elements such as sea-ice, iceberg drift, and |
---|
281 | surface fluxes to be run using prescribed model state fields. |
---|
282 | It can profitably be used to compare different bulk formulae or |
---|
283 | adjust the parameters of a given bulk formula. |
---|
284 | |
---|
285 | More informations about SAS can be found in NEMO manual. |
---|
286 | |
---|
287 | SPITZ12 |
---|
288 | ------- |
---|
289 | |
---|
290 | ``SPITZ12`` is a regional configuration around the Svalbard archipelago |
---|
291 | at 1/12° of horizontal resolution and 75 vertical levels. |
---|
292 | See `Rousset et al. (2015) <https://www.geosci-model-dev.net/8/2991/2015/>`_ for more details. |
---|
293 | |
---|
294 | This configuration references to year 2002, |
---|
295 | with atmospheric forcing provided every 2 hours using NCAR bulk formulae, |
---|
296 | while lateral boundary conditions for dynamical fields have 3 days time frequency. |
---|
297 | |
---|
298 | References |
---|
299 | ========== |
---|
300 | |
---|
301 | .. bibliography:: configurations.bib |
---|
302 | :all: |
---|
303 | :style: unsrt |
---|
304 | :labelprefix: C |
---|
305 | |
---|
306 | .. Links and substitutions |
---|
307 | |
---|