Changeset 11577 for NEMO/trunk/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex
- Timestamp:
- 2019-09-19T19:01:38+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/apdx_DOMAINcfg.tex
r11571 r11577 51 51 52 52 The user has three options available in defining a horizontal grid, which involve the 53 namelist variable \np{jphgr \_mesh} of the \nam{dom} (\texttt{DOMAINcfg} variant only)53 namelist variable \np{jphgr_mesh}{jphgr\_mesh} of the \nam{dom} (\texttt{DOMAINcfg} variant only) 54 54 namelist. 55 55 56 56 \begin{description} 57 \item[\np{jphgr \_mesh}=0] The most general curvilinear orthogonal grids.57 \item[\np{jphgr_mesh}{jphgr\_mesh}=0] The most general curvilinear orthogonal grids. 58 58 The coordinates and their first derivatives with respect to $i$ and $j$ are provided 59 59 in a input file (\ifile{coordinates}), read in \rou{hgr\_read} subroutine of the domhgr module. 60 60 This is now the only option available within \NEMO\ itself from v4.0 onwards. 61 \item[\np{jphgr \_mesh}=1 to 5] A few simple analytical grids are provided (see below).61 \item[\np{jphgr_mesh}{jphgr\_mesh}=1 to 5] A few simple analytical grids are provided (see below). 62 62 For other analytical grids, the \mdl{domhgr} module (\texttt{DOMAINcfg} variant) must be 63 63 modified by the user. In most cases, modifying the \mdl{usrdef\_hgr} module of \NEMO\ is … … 67 67 68 68 There are two simple cases of geographical grids on the sphere. With 69 \np{jphgr \_mesh}=1, the grid (expressed in degrees) is regular in space,70 with grid sizes specified by parameters \np{ppe1 \_deg} and \np{ppe2\_deg},69 \np{jphgr_mesh}{jphgr\_mesh}=1, the grid (expressed in degrees) is regular in space, 70 with grid sizes specified by parameters \np{ppe1_deg}{ppe1\_deg} and \np{ppe2_deg}{ppe2\_deg}, 71 71 respectively. Such a geographical grid can be very anisotropic at high latitudes 72 72 because of the convergence of meridians (the zonal scale factors $e_1$ 73 73 become much smaller than the meridional scale factors $e_2$). The Mercator 74 grid (\np{jphgr \_mesh}=4) avoids this anisotropy by refining the meridional scale74 grid (\np{jphgr_mesh}{jphgr\_mesh}=4) avoids this anisotropy by refining the meridional scale 75 75 factors in the same way as the zonal ones. In this case, meridional scale factors 76 76 and latitudes are calculated analytically using the formulae appropriate for 77 a Mercator projection, based on \np{ppe1 \_deg} which is a reference grid spacing77 a Mercator projection, based on \np{ppe1_deg}{ppe1\_deg} which is a reference grid spacing 78 78 at the equator (this applies even when the geographical equator is situated outside 79 79 the model domain). 80 80 81 In these two cases (\np{jphgr \_mesh}=1 or 4), the grid position is defined by the81 In these two cases (\np{jphgr_mesh}{jphgr\_mesh}=1 or 4), the grid position is defined by the 82 82 longitude and latitude of the south-westernmost point (\np{ppglamt0} 83 83 and \np{ppgphi0}). Note that for the Mercator grid the user need only provide … … 87 87 88 88 Rectangular grids ignoring the spherical geometry are defined with 89 \np{jphgr \_mesh} = 2, 3, 5. The domain is either an $f$-plane (\np{jphgr\_mesh} = 2,90 Coriolis factor is constant) or a beta-plane (\np{jphgr \_mesh} = 3, the Coriolis factor89 \np{jphgr_mesh}{jphgr\_mesh} = 2, 3, 5. The domain is either an $f$-plane (\np{jphgr\_mesh} = 2, 90 Coriolis factor is constant) or a beta-plane (\np{jphgr_mesh}{jphgr\_mesh} = 3, the Coriolis factor 91 91 is linear in the $j$-direction). The grid size is uniform in meter in each direction, 92 and given by the parameters \np{ppe1 \_m} and \np{ppe2\_m} respectively.92 and given by the parameters \np{ppe1_m}{ppe1\_m} and \np{ppe2_m}{ppe2\_m} respectively. 93 93 The zonal grid coordinate (\textit{glam} arrays) is in kilometers, starting at zero 94 94 with the first $t$-point. The meridional coordinate (gphi. arrays) is in kilometers, … … 97 97 latitude for computation of the Coriolis parameter. In the case of the beta plane, 98 98 \np{ppgphi0} corresponds to the center of the domain. Finally, the special case 99 \np{jphgr \_mesh}=5 corresponds to a beta plane in a rotated domain for the99 \np{jphgr_mesh}{jphgr\_mesh}=5 corresponds to a beta plane in a rotated domain for the 100 100 GYRE configuration, representing a classical mid-latitude double gyre system. 101 101 The rotation allows us to maximize the jet length relative to the gyre areas … … 170 170 \end{gather} 171 171 172 If the ice shelf cavities are opened (\np{ln \_isfcav}\forcode{ = .true.}), the definition172 If the ice shelf cavities are opened (\np{ln_isfcav}{ln\_isfcav}\forcode{ = .true.}), the definition 173 173 of $z_0$ is the same. However, definition of $e_3^0$ at $t$- and $w$-points is 174 174 respectively changed to: … … 312 312 313 313 Three options are possible for defining the bathymetry, according to the namelist variable 314 \np{nn \_bathy} (found in \nam{dom} namelist (\texttt{DOMAINCFG} variant) ):314 \np{nn_bathy}{nn\_bathy} (found in \nam{dom} namelist (\texttt{DOMAINCFG} variant) ): 315 315 \begin{description} 316 \item[\np{nn \_bathy}\forcode{ = 0}]:316 \item[\np{nn_bathy}{nn\_bathy}\forcode{ = 0}]: 317 317 a flat-bottom domain is defined. 318 318 The total depth $z_w (jpk)$ is given by the coordinate transformation. 319 319 The domain can either be a closed basin or a periodic channel depending on the parameter \np{jperio}. 320 \item[\np{nn \_bathy}\forcode{ = -1}]:320 \item[\np{nn_bathy}{nn\_bathy}\forcode{ = -1}]: 321 321 a domain with a bump of topography one third of the domain width at the central latitude. 322 322 This is meant for the "EEL-R5" configuration, a periodic or open boundary channel with a seamount. 323 \item[\np{nn \_bathy}\forcode{ = 1}]:323 \item[\np{nn_bathy}{nn\_bathy}\forcode{ = 1}]: 324 324 read a bathymetry and ice shelf draft (if needed). 325 325 The \ifile{bathy\_meter} file (Netcdf format) provides the ocean depth (positive, in meters) at … … 332 332 The \ifile{isfdraft\_meter} file (Netcdf format) provides the ice shelf draft (positive, in meters) at 333 333 each grid point of the model grid. 334 This file is only needed if \np{ln \_isfcav}\forcode{ = .true.}.334 This file is only needed if \np{ln_isfcav}{ln\_isfcav}\forcode{ = .true.}. 335 335 Defining the ice shelf draft will also define the ice shelf edge and the grounding line position. 336 336 \end{description} … … 359 359 % z-coordinate with constant thickness 360 360 % ------------------------------------------------------------------------------------------------------------- 361 \subsubsection[$Z$-coordinate with uniform thickness levels (\forcode{ln_zco})]{$Z$-coordinate with uniform thickness levels (\protect\np{ln \_zco})}361 \subsubsection[$Z$-coordinate with uniform thickness levels (\forcode{ln_zco})]{$Z$-coordinate with uniform thickness levels (\protect\np{ln_zco}{ln\_zco})} 362 362 \label{subsec:DOMCFG_zco} 363 363 … … 371 371 % z-coordinate with partial step 372 372 % ------------------------------------------------------------------------------------------------------------- 373 \subsubsection[$Z$-coordinate with partial step (\forcode{ln_zps})]{$Z$-coordinate with partial step (\protect\np{ln \_zps})}373 \subsubsection[$Z$-coordinate with partial step (\forcode{ln_zps})]{$Z$-coordinate with partial step (\protect\np{ln_zps}{ln\_zps})} 374 374 \label{subsec:DOMCFG_zps} 375 375 … … 394 394 $250~m$). Two variables in the namdom namelist are used to define the partial step 395 395 vertical grid. The mimimum water thickness (in meters) allowed for a cell partially 396 filled with bathymetry at level jk is the minimum of \np{rn \_e3zps\_min} (thickness in397 meters, usually $20~m$) or $e_{3t}(jk)*$\np{rn \_e3zps\_rat} (a fraction, usually 10\%, of396 filled with bathymetry at level jk is the minimum of \np{rn_e3zps_min}{rn\_e3zps\_min} (thickness in 397 meters, usually $20~m$) or $e_{3t}(jk)*$\np{rn_e3zps_rat}{rn\_e3zps\_rat} (a fraction, usually 10\%, of 398 398 the default thickness $e_{3t}(jk)$). 399 399 … … 401 401 % s-coordinate 402 402 % ------------------------------------------------------------------------------------------------------------- 403 \subsubsection[$S$-coordinate (\forcode{ln_sco})]{$S$-coordinate (\protect\np{ln \_sco})}403 \subsubsection[$S$-coordinate (\forcode{ln_sco})]{$S$-coordinate (\protect\np{ln_sco}{ln\_sco})} 404 404 \label{sec:DOMCFG_sco} 405 405 %------------------------------------------nam_zgr_sco--------------------------------------------------- … … 411 411 \end{listing} 412 412 %-------------------------------------------------------------------------------------------------------------- 413 Options are defined in \nam{zgr \_sco} (\texttt{DOMAINcfg} only).414 In $s$-coordinate (\np{ln \_sco}\forcode{ = .true.}), the depth and thickness of the model levels are defined from413 Options are defined in \nam{zgr_sco}{zgr\_sco} (\texttt{DOMAINcfg} only). 414 In $s$-coordinate (\np{ln_sco}{ln\_sco}\forcode{ = .true.}), the depth and thickness of the model levels are defined from 415 415 the product of a depth field and either a stretching function or its derivative, respectively: 416 416 … … 426 426 since a mixed step-like and bottom-following representation of the topography can be used 427 427 (\autoref{fig:DOM_z_zps_s_sps}) or an envelop bathymetry can be defined (\autoref{fig:DOM_z_zps_s_sps}). 428 The namelist parameter \np{rn \_rmax} determines the slope at which428 The namelist parameter \np{rn_rmax}{rn\_rmax} determines the slope at which 429 429 the terrain-following coordinate intersects the sea bed and becomes a pseudo z-coordinate. 430 The coordinate can also be hybridised by specifying \np{rn \_sbot\_min} and \np{rn\_sbot\_max} as430 The coordinate can also be hybridised by specifying \np{rn_sbot_min}{rn\_sbot\_min} and \np{rn_sbot_max}{rn\_sbot\_max} as 431 431 the minimum and maximum depths at which the terrain-following vertical coordinate is calculated. 432 432 … … 435 435 436 436 The original default \NEMO\ s-coordinate stretching is available if neither of the other options are specified as true 437 (\np{ln \_s\_SH94}\forcode{ = .false.} and \np{ln\_s\_SF12}\forcode{ = .false.}).437 (\np{ln_s_SH94}{ln\_s\_SH94}\forcode{ = .false.} and \np{ln_s_SF12}{ln\_s\_SF12}\forcode{ = .false.}). 438 438 This uses a depth independent $\tanh$ function for the stretching \citep{madec.delecluse.ea_JPO96}: 439 439 … … 455 455 456 456 A stretching function, 457 modified from the commonly used \citet{song.haidvogel_JCP94} stretching (\np{ln \_s\_SH94}\forcode{ = .true.}),457 modified from the commonly used \citet{song.haidvogel_JCP94} stretching (\np{ln_s_SH94}{ln\_s\_SH94}\forcode{ = .true.}), 458 458 is also available and is more commonly used for shelf seas modelling: 459 459 … … 476 476 %% %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 477 477 478 where $H_c$ is the critical depth (\np{rn \_hc}) at which the coordinate transitions from pure $\sigma$ to479 the stretched coordinate, and $\theta$ (\np{rn \_theta}) and $b$ (\np{rn\_bb}) are the surface and478 where $H_c$ is the critical depth (\np{rn_hc}{rn\_hc}) at which the coordinate transitions from pure $\sigma$ to 479 the stretched coordinate, and $\theta$ (\np{rn_theta}{rn\_theta}) and $b$ (\np{rn_bb}{rn\_bb}) are the surface and 480 480 bottom control parameters such that $0 \leqslant \theta \leqslant 20$, and $0 \leqslant b \leqslant 1$. 481 481 $b$ has been designed to allow surface and/or bottom increase of the vertical resolution 482 482 (\autoref{fig:DOMCFG_sco_function}). 483 483 484 Another example has been provided at version 3.5 (\np{ln \_s\_SF12}) that allows a fixed surface resolution in484 Another example has been provided at version 3.5 (\np{ln_s_SF12}{ln\_s\_SF12}) that allows a fixed surface resolution in 485 485 an analytical terrain-following stretching \citet{siddorn.furner_OM13}. 486 486 In this case the a stretching function $\gamma$ is defined such that: … … 504 504 505 505 This gives an analytical stretching of $\sigma$ that is solvable in $A$ and $B$ as a function of 506 the user prescribed stretching parameter $\alpha$ (\np{rn \_alpha}) that stretches towards506 the user prescribed stretching parameter $\alpha$ (\np{rn_alpha}{rn\_alpha}) that stretches towards 507 507 the surface ($\alpha > 1.0$) or the bottom ($\alpha < 1.0$) and 508 user prescribed surface (\np{rn \_zs}) and bottom depths.508 user prescribed surface (\np{rn_zs}{rn\_zs}) and bottom depths. 509 509 The bottom cell depth in this example is given as a function of water depth: 510 510 … … 514 514 \] 515 515 516 where the namelist parameters \np{rn \_zb\_a} and \np{rn\_zb\_b} are $a$ and $b$ respectively.516 where the namelist parameters \np{rn_zb_a}{rn\_zb\_a} and \np{rn_zb_b}{rn\_zb\_b} are $a$ and $b$ respectively. 517 517 518 518 %% %>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 542 542 the critical depth $h_c$. 543 543 In this example two options are available in depths shallower than $h_c$, 544 with pure sigma being applied if the \np{ln \_sigcrit} is true and pure z-coordinates if it is false544 with pure sigma being applied if the \np{ln_sigcrit}{ln\_sigcrit} is true and pure z-coordinates if it is false 545 545 (the z-coordinate being equal to the depths of the stretched coordinate at $h_c$). 546 546 … … 553 553 % z*- or s*-coordinate 554 554 % ------------------------------------------------------------------------------------------------------------- 555 \subsubsection[\zstar- or \sstar-coordinate (\forcode{ln_linssh})]{\zstar- or \sstar-coordinate (\protect\np{ln \_linssh})}555 \subsubsection[\zstar- or \sstar-coordinate (\forcode{ln_linssh})]{\zstar- or \sstar-coordinate (\protect\np{ln_linssh}{ln\_linssh})} 556 556 \label{subsec:DOMCFG_zgr_star} 557 557
Note: See TracChangeset
for help on using the changeset viewer.