1 | \documentclass[../main/NEMO_manual]{subfiles} |
---|
2 | |
---|
3 | \begin{document} |
---|
4 | |
---|
5 | \chapter{Configurations} |
---|
6 | \label{chap:CFGS} |
---|
7 | |
---|
8 | \chaptertoc |
---|
9 | |
---|
10 | \paragraph{Changes record} ~\\ |
---|
11 | |
---|
12 | {\footnotesize |
---|
13 | \begin{tabularx}{\textwidth}{l||X|X} |
---|
14 | Release & Author(s) & Modifications \\ |
---|
15 | \hline |
---|
16 | {\em 4.0} & {\em ...} & {\em ...} \\ |
---|
17 | {\em 3.6} & {\em ...} & {\em ...} \\ |
---|
18 | {\em 3.4} & {\em ...} & {\em ...} \\ |
---|
19 | {\em <=3.4} & {\em ...} & {\em ...} |
---|
20 | \end{tabularx} |
---|
21 | } |
---|
22 | |
---|
23 | \clearpage |
---|
24 | |
---|
25 | %% ================================================================================================= |
---|
26 | \section{Introduction} |
---|
27 | \label{sec:CFGS_intro} |
---|
28 | |
---|
29 | The purpose of this part of the manual is to introduce the \NEMO\ reference configurations. |
---|
30 | These configurations are offered as means to explore various numerical and physical options, |
---|
31 | thus allowing the user to verify that the code is performing in a manner consistent with that we are running. |
---|
32 | This form of verification is critical as one adopts the code for his or her particular research purposes. |
---|
33 | The reference configurations also provide a sense for some of the options available in the code, |
---|
34 | though by no means are all options exercised in the reference configurations. |
---|
35 | Configuration is defined manually through the \nam{cfg}{cfg} namelist variables. |
---|
36 | |
---|
37 | \begin{listing} |
---|
38 | \nlst{namcfg} |
---|
39 | \caption{\forcode{&namcfg}} |
---|
40 | \label{lst:namcfg} |
---|
41 | \end{listing} |
---|
42 | |
---|
43 | %% ================================================================================================= |
---|
44 | \section[C1D: 1D Water column model (\texttt{\textbf{key\_c1d}})]{C1D: 1D Water column model (\protect\key{c1d})} |
---|
45 | \label{sec:CFGS_c1d} |
---|
46 | |
---|
47 | The 1D model option simulates a stand alone water column within the 3D \NEMO\ system. |
---|
48 | It can be applied to the ocean alone or to the ocean-ice system and can include passive tracers or a biogeochemical model. |
---|
49 | It is set up by defining the position of the 1D water column in the grid |
---|
50 | (see \path{./cfgs/SHARED/namelist\_ref}). |
---|
51 | The 1D model is a very useful tool |
---|
52 | \textit{(a)} to learn about the physics and numerical treatment of vertical mixing processes; |
---|
53 | \textit{(b)} to investigate suitable parameterisations of unresolved turbulence |
---|
54 | (surface wave breaking, Langmuir circulation, ...); |
---|
55 | \textit{(c)} to compare the behaviour of different vertical mixing schemes; |
---|
56 | \textit{(d)} to perform sensitivity studies on the vertical diffusion at a particular point of an ocean domain; |
---|
57 | \textit{(d)} to produce extra diagnostics, without the large memory requirement of the full 3D model. |
---|
58 | |
---|
59 | The methodology is based on the configuration of the smallest possible domain: |
---|
60 | a 3x3 domain with 75 vertical levels. |
---|
61 | |
---|
62 | The 1D model has some specifies. First, all the horizontal derivatives are assumed to be zero, |
---|
63 | and second, the two components of the velocity are moved on a $T$-point. |
---|
64 | Therefore, defining \key{c1d} changes some things in the code behaviour: |
---|
65 | \begin{enumerate} |
---|
66 | \item a simplified \rou{stp} routine is used (\rou{stp\_c1d}, see \mdl{step\_c1d} module) in which |
---|
67 | both lateral tendancy terms and lateral physics are not called; |
---|
68 | \item the vertical velocity is zero |
---|
69 | (so far, no attempt at introducing a Ekman pumping velocity has been made); |
---|
70 | \item a simplified treatment of the Coriolis term is performed as $U$- and $V$-points are the same |
---|
71 | (see \mdl{dyncor\_c1d}). |
---|
72 | \end{enumerate} |
---|
73 | All the relevant \textit{\_c1d} modules can be found in the src/OCE/C1D directory of |
---|
74 | the \NEMO\ distribution. |
---|
75 | |
---|
76 | % to be added: a test case on the yearlong Ocean Weather Station (OWS) Papa dataset of Martin (1985) |
---|
77 | |
---|
78 | %% ================================================================================================= |
---|
79 | \section{ORCA family: global ocean with tripolar grid} |
---|
80 | \label{sec:CFGS_orca} |
---|
81 | |
---|
82 | The ORCA family is a series of global ocean configurations that are run together with |
---|
83 | the SI3 model (ORCA-ICE) and possibly with PISCES biogeochemical model (ORCA-ICE-PISCES). |
---|
84 | An appropriate namelist is available in \path{./cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg} for ORCA2. |
---|
85 | The domain of ORCA2 configuration is defined in \textit{ORCA\_R2\_zps\_domcfg.nc} file, |
---|
86 | this file is available in tar file on the \NEMO\ community zenodo platform: \\ |
---|
87 | https://doi.org/10.5281/zenodo.2640723 |
---|
88 | |
---|
89 | In this namelist\_cfg the name of domain input file is set in \nam{cfg}{cfg} block of namelist. |
---|
90 | |
---|
91 | \begin{figure}[!t] |
---|
92 | \centering |
---|
93 | \includegraphics[width=0.66\textwidth]{CFGS_ORCA_NH_mesh} |
---|
94 | \caption[ORCA mesh conception]{ |
---|
95 | ORCA mesh conception. |
---|
96 | The departure from an isotropic Mercator grid start poleward of 20\deg{N}. |
---|
97 | The two "north pole" are the foci of a series of embedded ellipses (blue curves) which |
---|
98 | are determined analytically and form the i-lines of the ORCA mesh (pseudo latitudes). |
---|
99 | Then, following \citet{madec.imbard_CD96}, |
---|
100 | the normal to the series of ellipses (red curves) is computed which |
---|
101 | provides the j-lines of the mesh (pseudo longitudes).} |
---|
102 | \label{fig:CFGS_ORCA_msh} |
---|
103 | \end{figure} |
---|
104 | |
---|
105 | %% ================================================================================================= |
---|
106 | \subsection{ORCA tripolar grid} |
---|
107 | \label{subsec:CFGS_orca_grid} |
---|
108 | |
---|
109 | The ORCA grid is a tripolar grid based on the semi-analytical method of \citet{madec.imbard_CD96}. |
---|
110 | It allows to construct a global orthogonal curvilinear ocean mesh which has no singularity point inside |
---|
111 | the computational domain since two north mesh poles are introduced and placed on lands. |
---|
112 | The method involves defining an analytical set of mesh parallels in the stereographic polar plan, |
---|
113 | computing the associated set of mesh meridians, and projecting the resulting mesh onto the sphere. |
---|
114 | The set of mesh parallels used is a series of embedded ellipses which foci are the two mesh north poles |
---|
115 | (\autoref{fig:CFGS_ORCA_msh}). |
---|
116 | The resulting mesh presents no loss of continuity in either the mesh lines or the scale factors, |
---|
117 | or even the scale factor derivatives over the whole ocean domain, as the mesh is not a composite mesh. |
---|
118 | \begin{figure}[!tbp] |
---|
119 | \centering |
---|
120 | \includegraphics[width=0.66\textwidth]{CFGS_ORCA_NH_msh05_e1_e2} |
---|
121 | \includegraphics[width=0.66\textwidth]{CFGS_ORCA_aniso} |
---|
122 | \caption[Horizontal scale factors and ratio of anisotropy for ORCA 0.5\deg\ mesh]{ |
---|
123 | \textit{Top}: Horizontal scale factors ($e_1$, $e_2$) and |
---|
124 | \textit{Bottom}: ratio of anisotropy ($e_1 / e_2$) |
---|
125 | for ORCA 0.5\deg\ mesh. |
---|
126 | South of 20\deg{N} a Mercator grid is used ($e_1 = e_2$) so that the anisotropy ratio is 1. |
---|
127 | Poleward of 20\deg{N}, |
---|
128 | the two "north pole" introduce a weak anisotropy over the ocean areas ($< 1.2$) except in |
---|
129 | vicinity of Victoria Island (Canadian Arctic Archipelago).} |
---|
130 | \label{fig:CFGS_ORCA_e1e2} |
---|
131 | \end{figure} |
---|
132 | |
---|
133 | The method is applied to Mercator grid (\ie\ same zonal and meridional grid spacing) poleward of 20\deg{N}, |
---|
134 | so that the Equator is a mesh line, which provides a better numerical solution for equatorial dynamics. |
---|
135 | The choice of the series of embedded ellipses (position of the foci and variation of the ellipses) |
---|
136 | is a compromise between maintaining the ratio of mesh anisotropy ($e_1 / e_2$) close to one in the ocean |
---|
137 | (especially in area of strong eddy activities such as the Gulf Stream) and keeping the smallest scale factor in |
---|
138 | the northern hemisphere larger than the smallest one in the southern hemisphere. |
---|
139 | The resulting mesh is shown in \autoref{fig:CFGS_ORCA_msh} and \autoref{fig:CFGS_ORCA_e1e2} for |
---|
140 | a half a degree grid (ORCA\_R05). |
---|
141 | The smallest ocean scale factor is found in along Antarctica, |
---|
142 | while the ratio of anisotropy remains close to one except near the Victoria Island in the Canadian Archipelago. |
---|
143 | |
---|
144 | %% ================================================================================================= |
---|
145 | \subsection{ORCA pre-defined resolution} |
---|
146 | \label{subsec:CFGS_orca_resolution} |
---|
147 | |
---|
148 | The \NEMO\ system is provided with five built-in ORCA configurations which differ in the horizontal resolution. |
---|
149 | The value of the resolution is given by the resolution at the Equator expressed in degrees. |
---|
150 | Each of configuration is set through the \textit{domain\_cfg} domain configuration file, |
---|
151 | which sets the grid size and configuration name parameters. |
---|
152 | The \NEMO\ System Team provides only ORCA2 domain input file "\textit{ORCA\_R2\_zps\_domcfg.nc}" file |
---|
153 | (\autoref{tab:CFGS_ORCA}). |
---|
154 | |
---|
155 | \begin{table}[!t] |
---|
156 | \centering |
---|
157 | \begin{tabular}{p{4cm} c c c c} |
---|
158 | Horizontal Grid & \texttt{ORCA\_index} & \texttt{jpiglo} & \texttt{jpjglo} \\ |
---|
159 | \hline \hline |
---|
160 | % 4 \deg\ & 4 & 92 & 76 \\ |
---|
161 | 2 \deg\ & 2 & 182 & 149 \\ |
---|
162 | 1 \deg\ & 1 & 362 & 292 \\ |
---|
163 | 0.5 \deg\ & 05 & 722 & 511 \\ |
---|
164 | 0.25\deg\ & 025 & 1442 & 1021 \\ |
---|
165 | \hline \hline |
---|
166 | \end{tabular} |
---|
167 | \caption[Domain size of ORCA family configurations]{ |
---|
168 | Domain size of ORCA family configurations. |
---|
169 | The flag for configurations of ORCA family need to be set in \textit{domain\_cfg} file.} |
---|
170 | \label{tab:CFGS_ORCA} |
---|
171 | \end{table} |
---|
172 | |
---|
173 | The ORCA\_R2 configuration has the following specificity: starting from a 2\deg\ ORCA mesh, |
---|
174 | local mesh refinements were applied to the Mediterranean, Red, Black and Caspian Seas, |
---|
175 | so that the resolution is 1\deg\ there. |
---|
176 | A local transformation were also applied with in the Tropics in order to refine the meridional resolution up to |
---|
177 | 0.5\deg\ at the Equator. |
---|
178 | |
---|
179 | The ORCA\_R1 configuration has only a local tropical transformation to refine the meridional resolution up to |
---|
180 | 1/3\deg\ at the Equator. |
---|
181 | Note that the tropical mesh refinements in ORCA\_R2 and R1 strongly increases the mesh anisotropy there. |
---|
182 | |
---|
183 | The ORCA\_R05 and higher global configurations do not incorporate any regional refinements. |
---|
184 | |
---|
185 | For ORCA\_R1 and R025, setting the configuration key to 75 allows to use 75 vertical levels, otherwise 46 are used. |
---|
186 | In the other ORCA configurations, 31 levels are used |
---|
187 | (see \autoref{tab:CFGS_ORCA}). %\sfcomment{HERE I need to put new table for ORCA2 values} and \autoref{fig:DOM_zgr_e3}). |
---|
188 | |
---|
189 | Only the ORCA\_R2 is provided with all its input files in the \NEMO\ distribution. |
---|
190 | %It is very similar to that used as part of the climate model developed at IPSL for the 4th IPCC assessment of |
---|
191 | %climate change (Marti et al., 2009). |
---|
192 | %It is also the basis for the \NEMO\ contribution to the Coordinate Ocean-ice Reference Experiments (COREs) |
---|
193 | %documented in \citet{griffies.biastoch.ea_OM09}. |
---|
194 | |
---|
195 | This version of ORCA\_R2 has 31 levels in the vertical, with the highest resolution (10m) in the upper 150m |
---|
196 | (see \autoref{tab:CFGS_ORCA} and \autoref{fig:DOM_zgr_e3}). |
---|
197 | The bottom topography and the coastlines are derived from the global atlas of Smith and Sandwell (1997). |
---|
198 | The default forcing uses the boundary forcing from \citet{large.yeager_trpt04} (see \autoref{subsec:SBC_blk_ocean}), |
---|
199 | which was developed for the purpose of running global coupled ocean-ice simulations without |
---|
200 | an interactive atmosphere. |
---|
201 | This \citet{large.yeager_trpt04} dataset is available through |
---|
202 | the \href{http://nomads.gfdl.noaa.gov/nomads/forms/mom4/CORE.html}{GFDL web site}. |
---|
203 | The "normal year" of \citet{large.yeager_trpt04} has been chosen of the \NEMO\ distribution since release v3.3. |
---|
204 | |
---|
205 | ORCA\_R2 pre-defined configuration can also be run with multiply online nested zooms (\ie\ with AGRIF, \key{agrif} defined). |
---|
206 | This is available as the AGRIF\_DEMO configuration that can be found in the \path{./cfgs/AGRIF_DEMO/} directory. |
---|
207 | |
---|
208 | A regional Arctic or peri-Antarctic configuration is extracted from an ORCA\_R2 or R05 configurations using |
---|
209 | sponge layers at open boundaries. |
---|
210 | |
---|
211 | %% ================================================================================================= |
---|
212 | \section{GYRE family: double gyre basin} |
---|
213 | \label{sec:CFGS_gyre} |
---|
214 | |
---|
215 | The GYRE configuration \citep{levy.klein.ea_OM10} has been built to |
---|
216 | simulate the seasonal cycle of a double-gyre box model. |
---|
217 | It consists in an idealized domain similar to that used in the studies of \citet{drijfhout_JPO94} and |
---|
218 | \citet{hazeleger.drijfhout_JPO98, hazeleger.drijfhout_JPO99, hazeleger.drijfhout_JGR00, hazeleger.drijfhout_JPO00}, |
---|
219 | over which an analytical seasonal forcing is applied. |
---|
220 | This allows to investigate the spontaneous generation of a large number of interacting, transient mesoscale eddies |
---|
221 | and their contribution to the large scale circulation. |
---|
222 | |
---|
223 | The GYRE configuration run together with the PISCES biogeochemical model (GYRE-PISCES). |
---|
224 | The domain geometry is a closed rectangular basin on the $\beta$-plane centred at $\sim$ 30\deg{N} and |
---|
225 | rotated by 45\deg, 3180~km long, 2120~km wide and 4~km deep (\autoref{fig:MISC_strait_hand}). |
---|
226 | The domain is bounded by vertical walls and by a flat bottom. |
---|
227 | The configuration is meant to represent an idealized North Atlantic or North Pacific basin. |
---|
228 | The circulation is forced by analytical profiles of wind and buoyancy fluxes. |
---|
229 | The applied forcings vary seasonally in a sinusoidal manner between winter and summer extrema \citep{levy.klein.ea_OM10}. |
---|
230 | The wind stress is zonal and its curl changes sign at 22\deg{N} and 36\deg{N}. |
---|
231 | It forces a subpolar gyre in the north, a subtropical gyre in the wider part of the domain and |
---|
232 | a small recirculation gyre in the southern corner. |
---|
233 | The net heat flux takes the form of a restoring toward a zonal apparent air temperature profile. |
---|
234 | A portion of the net heat flux which comes from the solar radiation is allowed to penetrate within the water column. |
---|
235 | The fresh water flux is also prescribed and varies zonally. |
---|
236 | It is determined such as, at each time step, the basin-integrated flux is zero. |
---|
237 | The basin is initialised at rest with vertical profiles of temperature and salinity uniformly applied to |
---|
238 | the whole domain. |
---|
239 | |
---|
240 | The GYRE configuration is set like an analytical configuration. |
---|
241 | Through \np[=.false.]{ln_read_cfg}{ln\_read\_cfg} in \nam{cfg}{cfg} namelist defined in |
---|
242 | the reference configuration \path{./cfgs/GYRE_PISCES/EXPREF/namelist_cfg} |
---|
243 | analytical definition of grid in GYRE is done in mdl{usrdef\_hrg}, \mdl{usrdef\_zgr} routines. |
---|
244 | Its horizontal resolution (and thus the size of the domain) is determined by |
---|
245 | setting \np{nn_GYRE}{nn\_GYRE} in \nam{usr_def}{usr\_def}: |
---|
246 | |
---|
247 | \begin{align*} |
---|
248 | jpiglo = 30 \times \text{\np{nn_GYRE}{nn\_GYRE}} + 2 + 2 \times \text{\np{nn_hls}{nn\_hls}} \\ |
---|
249 | jpjglo = 20 \times \text{\np{nn_GYRE}{nn\_GYRE}} + 2 + 2 \times \text{\np{nn_hls}{nn\_hls}} |
---|
250 | \end{align*} |
---|
251 | |
---|
252 | Obviously, the namelist parameters have to be adjusted to the chosen resolution, |
---|
253 | see the Configurations pages on the \NEMO\ web site (\NEMO\ Configurations). |
---|
254 | In the vertical, GYRE uses the default 30 ocean levels (\forcode{jpk = 31}, \autoref{fig:DOM_zgr_e3}). |
---|
255 | |
---|
256 | \begin{listing} |
---|
257 | \begin{forlines} |
---|
258 | !----------------------------------------------------------------------- |
---|
259 | &namusr_def ! GYRE user defined namelist |
---|
260 | !----------------------------------------------------------------------- |
---|
261 | nn_GYRE = 1 ! GYRE resolution [1/degrees] |
---|
262 | ln_bench = .false. ! ! =T benchmark with gyre: the gridsize is kept constant |
---|
263 | jpkglo = 31 ! number of model levels |
---|
264 | / |
---|
265 | \end{forlines} |
---|
266 | \caption{\forcode{&namusr_def}} |
---|
267 | \label{lst:namusr_def} |
---|
268 | \end{listing} |
---|
269 | |
---|
270 | The GYRE configuration is also used in benchmark test as it is very simple to increase its resolution and |
---|
271 | as it does not requires any input file. |
---|
272 | For example, keeping a same model size on each processor while increasing the number of processor used is very easy, |
---|
273 | even though the physical integrity of the solution can be compromised. |
---|
274 | Benchmark is activate via \np[=.true.]{ln_bench}{ln\_bench} in \nam{usr_def}{usr\_def} in |
---|
275 | namelist \path{./cfgs/GYRE_PISCES/EXPREF/namelist_cfg}. |
---|
276 | |
---|
277 | \begin{figure}[!t] |
---|
278 | \centering |
---|
279 | \includegraphics[width=0.66\textwidth]{CFGS_GYRE} |
---|
280 | \caption[Snapshot of relative vorticity at the surface of the model domain in GYRE R9, R27 and R54]{ |
---|
281 | Snapshot of relative vorticity at the surface of the model domain in GYRE R9, R27 and R54. |
---|
282 | From \citet{levy.klein.ea_OM10}.} |
---|
283 | \label{fig:CFGS_GYRE} |
---|
284 | \end{figure} |
---|
285 | |
---|
286 | %% ================================================================================================= |
---|
287 | \section{AMM: atlantic margin configuration} |
---|
288 | \label{sec:CFGS_config_AMM} |
---|
289 | |
---|
290 | The AMM, Atlantic Margins Model, is a regional model covering the Northwest European Shelf domain on |
---|
291 | a regular lat-lon grid at approximately 12km horizontal resolution. |
---|
292 | The appropriate \textit{\&namcfg} namelist is available in \path{./cfgs/AMM12/EXPREF/namelist\_cfg}. |
---|
293 | It is used to build the correct dimensions of the AMM domain. |
---|
294 | |
---|
295 | This configuration tests several features of \NEMO\ functionality specific to the shelf seas. |
---|
296 | In particular, the AMM uses $s$-coordinates in the vertical rather than $z$-coordinates and |
---|
297 | is forced with tidal lateral boundary conditions using a Flather boundary condition from the BDY module. |
---|
298 | Also specific to the AMM configuration is the use of the GLS turbulence scheme (\np[=.true.]{ln_zdfgls}{ln\_zdfgls}). |
---|
299 | |
---|
300 | In addition to the tidal boundary condition the model may also take open boundary conditions from |
---|
301 | a North Atlantic model. |
---|
302 | Boundaries may be completely omitted by setting \np{ln_bdy}{ln\_bdy} to false. |
---|
303 | Sample surface fluxes, river forcing and a sample initial restart file are included to test a realistic model run. |
---|
304 | The Baltic boundary is included within the river input file and is specified as a river source. |
---|
305 | Unlike ordinary river points the Baltic inputs also include salinity and temperature data. |
---|
306 | |
---|
307 | \subinc{\input{../../global/epilogue}} |
---|
308 | |
---|
309 | \end{document} |
---|