[10598] | 1 | ************************ |
---|
| 2 | Reference configurations |
---|
| 3 | ************************ |
---|
[10186] | 4 | |
---|
[10605] | 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). |
---|
[10186] | 8 | |
---|
[10605] | 9 | .. attention:: |
---|
[10186] | 10 | |
---|
[10605] | 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 | |
---|
[10541] | 45 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
| 46 | Configuration Component(s) Input & Forcing File(s) |
---|
| 47 | ---------------------- ---------------------------------- ------------------------------------------------ |
---|
| 48 | Name OPA SI3 TOP PISCES AGRIF |
---|
| 49 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
[10605] | 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 |
---|
[10541] | 53 | GYRE_BFM_ X X *none* |
---|
| 54 | GYRE_PISCES_ X X X *none* |
---|
[10605] | 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 |
---|
[10541] | 60 | ====================== ===== ===== ===== ======== ======= ================================================ |
---|
[10186] | 61 | |
---|
[10541] | 62 | AGRIF_DEMO |
---|
| 63 | ---------- |
---|
[10186] | 64 | |
---|
[10605] | 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). |
---|
[10541] | 68 | |
---|
[10605] | 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. |
---|
[10541] | 73 | |
---|
[10605] | 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. |
---|
[10541] | 79 | |
---|
[10186] | 80 | AMM12 |
---|
| 81 | ----- |
---|
| 82 | |
---|
[10605] | 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`). |
---|
[10186] | 86 | |
---|
[10541] | 87 | This configuration allows to tests several features of NEMO specifically addressed to the shelf seas. |
---|
[10605] | 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``). |
---|
[10541] | 90 | Boundaries may be completely omitted by setting ``ln_bdy = .false.`` in ``nambdy``. |
---|
[10186] | 91 | |
---|
[10605] | 92 | Sample surface fluxes, river forcing and an initial restart file are included to test a realistic model run |
---|
| 93 | (``AMM12_v4.0.tar``). |
---|
[10186] | 94 | |
---|
[10605] | 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. |
---|
[10186] | 97 | |
---|
| 98 | C1D_PAPA |
---|
| 99 | -------- |
---|
| 100 | |
---|
[10605] | 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. |
---|
[10186] | 103 | |
---|
[10605] | 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. |
---|
[10186] | 106 | |
---|
[10541] | 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. |
---|
[10605] | 108 | Data provided with ``INPUTS_C1D_PAPA_v4.0.tar`` file account for: |
---|
[10186] | 109 | |
---|
[10541] | 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 |
---|
[10186] | 113 | |
---|
| 114 | GYRE_BFM |
---|
| 115 | -------- |
---|
| 116 | |
---|
[10605] | 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``. |
---|
[10186] | 118 | |
---|
| 119 | GYRE_PISCES |
---|
| 120 | ----------- |
---|
| 121 | |
---|
[10605] | 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. |
---|
[10186] | 126 | |
---|
[10605] | 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``: |
---|
[10186] | 129 | |
---|
[10605] | 130 | .. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_cfg |
---|
| 131 | :language: fortran |
---|
| 132 | :lines: 34-42 |
---|
[10186] | 133 | |
---|
[10605] | 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``. |
---|
[10186] | 137 | |
---|
[10605] | 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. |
---|
[10186] | 141 | |
---|
[10605] | 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: |
---|
[10186] | 143 | |
---|
[10541] | 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 | =========== ========= ========== ============ =================== |
---|
[10186] | 153 | |
---|
[10605] | 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. |
---|
[10186] | 157 | |
---|
[10605] | 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 |
---|
[10541] | 160 | |
---|
[10605] | 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. |
---|
[10186] | 164 | |
---|
| 165 | ORCA2_ICE_PISCES |
---|
| 166 | ---------------- |
---|
| 167 | |
---|
[10605] | 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 |
---|
[10186] | 175 | |
---|
[10605] | 176 | this configuration uses the three components |
---|
[10186] | 177 | |
---|
[10605] | 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` |
---|
[10186] | 181 | |
---|
[10541] | 182 | All components share the same grid. |
---|
[10186] | 183 | |
---|
[10605] | 184 | The model is forced with CORE-II normal year atmospheric forcing and |
---|
| 185 | it uses the NCAR bulk formulae. |
---|
[10186] | 186 | |
---|
[10605] | 187 | In this ``ORCA2_ICE_PISCES`` configuration, |
---|
| 188 | AGRIF nesting can be activated that includes a nested grid in the Agulhas region. |
---|
[10186] | 189 | |
---|
[10605] | 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 | |
---|
[10541] | 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. |
---|
[10186] | 209 | The Beckmann [19XX] simple bottom boundary layer parameterization is applied along continental slopes. |
---|
| 210 | A linear friction is applied on momentum. |
---|
[10541] | 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. |
---|
[10186] | 213 | |
---|
| 214 | ORCA2_OFF_PISCES |
---|
| 215 | ---------------- |
---|
| 216 | |
---|
[10605] | 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. |
---|
[10186] | 219 | |
---|
| 220 | |
---|
| 221 | ORCA2_OFF_TRC |
---|
| 222 | ------------- |
---|
| 223 | |
---|
[10605] | 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. |
---|
[10186] | 226 | |
---|
[10605] | 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. |
---|
[10186] | 230 | |
---|
[10605] | 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): |
---|
[10186] | 234 | |
---|
[10605] | 235 | .. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_ref |
---|
| 236 | :language: fortran |
---|
[11303] | 237 | :lines: 935-960 |
---|
[10186] | 238 | |
---|
[10605] | 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. |
---|
[10541] | 241 | |
---|
[10605] | 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. |
---|
[10541] | 246 | |
---|
[10605] | 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). |
---|
[10541] | 250 | |
---|
| 251 | In addition, the offline module (OFF) allows for the provision of further fields: |
---|
| 252 | |
---|
[10605] | 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: |
---|
[10541] | 255 | |
---|
| 256 | .. code-block:: fortran |
---|
| 257 | |
---|
[10605] | 258 | sn_rnf = 'dyna_grid_T', 120, 'sorunoff' , .true., .true., 'yearly', '', '', '' |
---|
[10541] | 259 | |
---|
[10605] | 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 |
---|
[10541] | 264 | |
---|
| 265 | .. code-block:: fortran |
---|
| 266 | |
---|
[10605] | 267 | sn_div = 'dyna_grid_T', 120, 'e3t' , .true., .true., 'yearly', '', '', '' |
---|
| 268 | sn_empb = 'dyna_grid_T', 120, 'sowaflupb', .true., .true., 'yearly', '', '', '' |
---|
[10541] | 269 | |
---|
| 270 | |
---|
[10605] | 271 | More details can be found by inspecting the offline data manager in |
---|
| 272 | the routine ``./src/OFF/dtadyn.F90``. |
---|
[10541] | 273 | |
---|
[10186] | 274 | ORCA2_SAS_ICE |
---|
| 275 | ------------- |
---|
| 276 | |
---|
[10605] | 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). |
---|
[10186] | 279 | |
---|
[10605] | 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. |
---|
[10186] | 284 | |
---|
[10541] | 285 | More informations about SAS can be found in NEMO manual. |
---|
[10186] | 286 | |
---|
| 287 | SPITZ12 |
---|
| 288 | ------- |
---|
| 289 | |
---|
[10605] | 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. |
---|
[10186] | 293 | |
---|
[10605] | 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. |
---|
[10186] | 297 | |
---|
[10605] | 298 | References |
---|
| 299 | ========== |
---|
[10186] | 300 | |
---|
[10605] | 301 | .. bibliography:: configurations.bib |
---|
| 302 | :all: |
---|
| 303 | :style: unsrt |
---|
| 304 | :labelprefix: C |
---|
[10541] | 305 | |
---|
[10605] | 306 | .. Links and substitutions |
---|
[10541] | 307 | |
---|