Changeset 11435 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex
- Timestamp:
- 2019-08-14T14:45:08+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex
r11386 r11435 9 9 \label{chap:LDF} 10 10 11 \ minitoc11 \chaptertoc 12 12 13 13 \newpage … … 22 22 (3) the space and time variations of the eddy coefficients. 23 23 These three aspects of the lateral diffusion are set through namelist parameters 24 (see the \n gn{nam\_traldf} and \ngn{nam\_dynldf} below).24 (see the \nam{tra\_ldf} and \nam{dyn\_ldf} below). 25 25 Note that this chapter describes the standard implementation of iso-neutral tracer mixing. 26 26 Griffies's implementation, which is used if \np{ln\_traldf\_triad}\forcode{ = .true.}, … … 53 53 {Laplacian mixing (\protect\np{ln\_traldf\_lap}, \protect\np{ln\_dynldf\_lap})} 54 54 Setting \protect\np{ln\_traldf\_lap}\forcode{ = .true.} and/or \protect\np{ln\_dynldf\_lap}\forcode{ = .true.} enables 55 a second order diffusion on tracers and momentum respectively. Note that in \NEMO 4, one can not combine55 a second order diffusion on tracers and momentum respectively. Note that in \NEMO\ 4, one can not combine 56 56 Laplacian and Bilaplacian operators for the same variable. 57 57 … … 60 60 Setting \protect\np{ln\_traldf\_blp}\forcode{ = .true.} and/or \protect\np{ln\_dynldf\_blp}\forcode{ = .true.} enables 61 61 a fourth order diffusion on tracers and momentum respectively. It is implemented by calling the above Laplacian operator twice. 62 We stress again that from \NEMO 4, the simultaneous use Laplacian and Bilaplacian operators is not allowed.62 We stress again that from \NEMO\ 4, the simultaneous use Laplacian and Bilaplacian operators is not allowed. 63 63 64 64 % ================================================================ … … 90 90 \subsection{Slopes for tracer geopotential mixing in the $s$-coordinate} 91 91 92 In $s$-coordinates, geopotential mixing (\ie horizontal mixing) $r_1$ and $r_2$ are the slopes between92 In $s$-coordinates, geopotential mixing (\ie\ horizontal mixing) $r_1$ and $r_2$ are the slopes between 93 93 the geopotential and computational surfaces. 94 94 Their discrete formulation is found by locally solving \autoref{eq:tra_ldf_iso} when 95 95 the diffusive fluxes in the three directions are set to zero and $T$ is assumed to be horizontally uniform, 96 \ie a linear function of $z_T$, the depth of a $T$-point.96 \ie\ a linear function of $z_T$, the depth of a $T$-point. 97 97 %gm { Steven : My version is obviously wrong since I'm left with an arbitrary constant which is the local vertical temperature gradient} 98 98 … … 124 124 Their formulation does not depend on the vertical coordinate used. 125 125 Their discrete formulation is found using the fact that the diffusive fluxes of 126 locally referenced potential density (\ie $in situ$ density) vanish.126 locally referenced potential density (\ie\ $in situ$ density) vanish. 127 127 So, substituting $T$ by $\rho$ in \autoref{eq:tra_ldf_iso} and setting the diffusive fluxes in 128 128 the three directions to zero leads to the following definition for the neutral slopes: … … 230 230 To overcome this problem, several techniques have been proposed in which the numerical schemes of 231 231 the ocean model are modified \citep{weaver.eby_JPO97, griffies.gnanadesikan.ea_JPO98}. 232 Griffies's scheme is now available in \NEMO if \np{ln\_traldf\_triad}=\forcode{= .true.}; see \autoref{apdx:triad}.232 Griffies's scheme is now available in \NEMO\ if \np{ln\_traldf\_triad}\forcode{ = .true.}; see \autoref{apdx:triad}. 233 233 Here, another strategy is presented \citep{lazar_phd97}: 234 234 a local filtering of the iso-neutral slopes (made on 9 grid-points) prevents the development of … … 284 284 \textit{(a)} in the real ocean the slope is the iso-neutral slope in the ocean interior, 285 285 which has to be adjusted at the surface boundary 286 \ie it must tend to zero at the surface since there is no mixing across the air-sea interface:286 \ie\ it must tend to zero at the surface since there is no mixing across the air-sea interface: 287 287 wall boundary condition). 288 288 Nevertheless, the profile between the surface zero value and the interior iso-neutral one is unknown, … … 309 309 \textit{vw}- points for the $v$-component. 310 310 They are computed from the slopes used for tracer diffusion, 311 \ie \autoref{eq:ldfslp_geo} and \autoref{eq:ldfslp_iso}:311 \ie\ \autoref{eq:ldfslp_geo} and \autoref{eq:ldfslp_iso}: 312 312 313 313 \[ … … 323 323 The major issue remaining is in the specification of the boundary conditions. 324 324 The same boundary conditions are chosen as those used for lateral diffusion along model level surfaces, 325 \ie using the shear computed along the model levels and with no additional friction at the ocean bottom325 \ie\ using the shear computed along the model levels and with no additional friction at the ocean bottom 326 326 (see \autoref{sec:LBC_coast}). 327 327 … … 420 420 421 421 The 3D space variation of the mixing coefficient is simply the combination of the 1D and 2D cases above, 422 \ie a hyperbolic tangent variation with depth associated with a grid size dependence of422 \ie\ a hyperbolic tangent variation with depth associated with a grid size dependence of 423 423 the magnitude of the coefficient. 424 424 … … 527 527 the formulation of which depends on the slopes of iso-neutral surfaces. 528 528 Contrary to the case of iso-neutral mixing, the slopes used here are referenced to the geopotential surfaces, 529 \ie \autoref{eq:ldfslp_geo} is used in $z$-coordinates,529 \ie\ \autoref{eq:ldfslp_geo} is used in $z$-coordinates, 530 530 and the sum \autoref{eq:ldfslp_geo} + \autoref{eq:ldfslp_iso} in $s$-coordinates. 531 531 532 If isopycnal mixing is used in the standard way, \ie \np{ln\_traldf\_triad}\forcode{ = .false.}, the eddy induced velocity is given by:532 If isopycnal mixing is used in the standard way, \ie\ \np{ln\_traldf\_triad}\forcode{ = .false.}, the eddy induced velocity is given by: 533 533 \begin{equation} 534 534 \label{eq:ldfeiv} … … 539 539 \end{split} 540 540 \end{equation} 541 where $A^{eiv}$ is the eddy induced velocity coefficient whose value is set through \np{nn\_aei\_ijk\_t} \n gn{namtra\_eiv} namelist parameter.541 where $A^{eiv}$ is the eddy induced velocity coefficient whose value is set through \np{nn\_aei\_ijk\_t} \nam{tra\_eiv} namelist parameter. 542 542 The three components of the eddy induced velocity are computed in \rou{ldf\_eiv\_trp} and 543 543 added to the eulerian velocity in \rou{tra\_adv} where tracer advection is performed. … … 570 570 %-------------------------------------------------------------------------------------------------------------- 571 571 572 If \np{ln\_mle}\forcode{ = .true.} in \n gn{namtra\_mle} namelist, a parameterization of the mixing due to unresolved mixed layer instabilities is activated (\citet{foxkemper.ferrari_JPO08}). Additional transport is computed in \rou{ldf\_mle\_trp} and added to the eulerian transport in \rou{tra\_adv} as done for eddy induced advection.572 If \np{ln\_mle}\forcode{ = .true.} in \nam{tra\_mle} namelist, a parameterization of the mixing due to unresolved mixed layer instabilities is activated (\citet{foxkemper.ferrari_JPO08}). Additional transport is computed in \rou{ldf\_mle\_trp} and added to the eulerian transport in \rou{tra\_adv} as done for eddy induced advection. 573 573 574 574 \colorbox{yellow}{TBC}
Note: See TracChangeset
for help on using the changeset viewer.