Changeset 10442 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex
- Timestamp:
- 2018-12-21T15:18:38+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_DOM.tex
r10414 r10442 64 64 the barotropic stream function $\psi$ is defined at horizontal points overlying the $\zeta$ and $f$-points. 65 65 66 The ocean mesh ( $i.e.$the position of all the scalar and vector points) is defined by66 The ocean mesh (\ie the position of all the scalar and vector points) is defined by 67 67 the transformation that gives ($\lambda$ ,$\varphi$ ,$z$) as a function of $(i,j,k)$. 68 68 The grid-points are located at integer or integer and a half value of $(i,j,k)$ as indicated on \autoref{tab:cell}. … … 162 162 163 163 The vertical average over the whole water column denoted by an overbar becomes for a quantity $q$ which 164 is a masked field ( i.e.equal to zero inside solid area):164 is a masked field (\ie equal to zero inside solid area): 165 165 \begin{equation} 166 166 \label{eq:DOM_bar} … … 191 191 the differencing operators ($\delta_i$, $\delta_j$ and $\delta_k$) are skew-symmetric linear operators, 192 192 and further that the averaging operators $\overline{\,\cdot\,}^{\,i}$, $\overline{\,\cdot\,}^{\,k}$ and 193 $\overline{\,\cdot\,}^{\,k}$) are symmetric linear operators, 194 $i.e.$ 193 $\overline{\,\cdot\,}^{\,k}$) are symmetric linear operators, \ie 195 194 \begin{align} 196 195 \label{eq:DOM_di_adj} … … 219 218 \caption{ 220 219 \protect\label{fig:index_hor} 221 Horizontal integer indexing used in the \ textsc{Fortran}code.220 Horizontal integer indexing used in the \fortran code. 222 221 The dashed area indicates the cell in which variables contained in arrays have the same $i$- and $j$-indices 223 222 } … … 226 225 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 227 226 228 The array representation used in the \ textsc{Fortran}code requires an integer indexing while227 The array representation used in the \fortran code requires an integer indexing while 229 228 the analytical definition of the mesh (see \autoref{subsec:DOM_cell}) is associated with the use of 230 229 integer values for $t$-points and both integer and integer and a half values for all the other points. 231 230 Therefore a specific integer indexing must be defined for points other than $t$-points 232 ( $i.e.$velocity and vorticity grid-points).231 (\ie velocity and vorticity grid-points). 233 232 Furthermore, the direction of the vertical indexing has been changed so that the surface level is at $k=1$. 234 233 … … 252 251 253 252 In the vertical, the chosen indexing requires special attention since 254 the $k$-axis is re-orientated downward in the \ textsc{Fortran}code compared to253 the $k$-axis is re-orientated downward in the \fortran code compared to 255 254 the indexing used in the semi-discrete equations and given in \autoref{subsec:DOM_cell}. 256 255 The sea surface corresponds to the $w$-level $k=1$ which is the same index as $t$-level just below … … 263 262 have the same $i$ or $j$ index 264 263 (compare the dashed area in \autoref{fig:index_hor} and \autoref{fig:index_vert}). 265 Since the scale factors are chosen to be strictly positive, a \emph{minus sign} appears in the \ textsc{Fortran}264 Since the scale factors are chosen to be strictly positive, a \emph{minus sign} appears in the \fortran 266 265 code \emph{before all the vertical derivatives} of the discrete equations given in this documentation. 267 266 … … 272 271 \caption{ 273 272 \protect\label{fig:index_vert} 274 Vertical integer indexing used in the \ textsc{Fortran }code.273 Vertical integer indexing used in the \fortran code. 275 274 Note that the $k$-axis is orientated downward. 276 275 The dashed area indicates the cell in which variables contained in arrays have the same $k$-index. … … 300 299 \section{Needed fields} 301 300 \label{sec:DOM_fields} 302 The ocean mesh ( $i.e.$the position of all the scalar and vector points) is defined by the transformation that gives $(\lambda,\varphi,z)$ as a function of $(i,j,k)$.301 The ocean mesh (\ie the position of all the scalar and vector points) is defined by the transformation that gives $(\lambda,\varphi,z)$ as a function of $(i,j,k)$. 303 302 The grid-points are located at integer or integer and a half values of as indicated in \autoref{tab:cell}. 304 303 The associated scale factors are defined using the analytical first derivative of the transformation … … 352 351 \label{subsec:DOM_hgr_coord_e} 353 352 354 The ocean mesh ( $i.e.$the position of all the scalar and vector points) is defined by353 The ocean mesh (\ie the position of all the scalar and vector points) is defined by 355 354 the transformation that gives $(\lambda,\varphi,z)$ as a function of $(i,j,k)$. 356 355 The grid-points are located at integer or integer and a half values of as indicated in \autoref{tab:cell}. … … 391 390 392 391 Note that the definition of the scale factors 393 ( $i.e.$as the analytical first derivative of the transformation that392 (\ie as the analytical first derivative of the transformation that 394 393 gives $(\lambda,\varphi,z)$ as a function of $(i,j,k)$) 395 394 is specific to the \NEMO model \citep{Marti_al_JGR92}. … … 461 460 (2) the number of levels of the model (\jp{jpk}); 462 461 (3) the analytical transformation $z(i,j,k)$ and the vertical scale factors (derivatives of the transformation); and 463 (4) the masking system, $i.e.$the number of wet model levels at each462 (4) the masking system, \ie the number of wet model levels at each 464 463 $(i,j)$ column of points. 465 464 … … 563 562 The \ifile{bathy\_meter} file (Netcdf format) provides the ocean depth (positive, in meters) at 564 563 each grid point of the model grid. 565 The bathymetry is usually built by interpolating a standard bathymetry product ( $e.g.$ETOPO2) onto564 The bathymetry is usually built by interpolating a standard bathymetry product (\eg ETOPO2) onto 566 565 the horizontal ocean mesh. 567 566 Defining the bathymetry also defines the coastline: where the bathymetry is zero, … … 926 925 927 926 % ------------------------------------------------------------------------------------------------------------- 928 % z*- or s*-coordinate927 % \zstar- or \sstar-coordinate 929 928 % ------------------------------------------------------------------------------------------------------------- 930 929 \subsection{$Z^*$- or $S^*$-coordinate (\protect\np{ln\_linssh}\forcode{ = .false.}) } … … 945 944 follow the face of the model cells (step like topography) \citep{Madec_al_JPO96}. 946 945 The distribution of the steps in the horizontal is defined in a 2D integer array, mbathy, 947 which gives the number of ocean levels ( $i.e.$those that are not masked) at each $t$-point.946 which gives the number of ocean levels (\ie those that are not masked) at each $t$-point. 948 947 mbathy is computed from the meter bathymetry using the definiton of gdept as 949 948 the number of $t$-points which gdept $\leq$ bathy. … … 961 960 the cavities are performed in the \textit{zgr\_isf} routine. 962 961 The compatibility between ice shelf draft and bathymetry is checked. 963 All the locations where the isf cavity is thinnest than \np{rn\_isfhmin} meters are grounded ( $i.e.$masked).962 All the locations where the isf cavity is thinnest than \np{rn\_isfhmin} meters are grounded (\ie masked). 964 963 If only one cell on the water column is opened at $t$-, $u$- or $v$-points, 965 964 the bathymetry or the ice shelf draft is dug to fit this constrain. … … 1017 1016 \biblio 1018 1017 1018 \pindex 1019 1019 1020 \end{document}
Note: See TracChangeset
for help on using the changeset viewer.