New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11435 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex – NEMO

Ignore:
Timestamp:
2019-08-14T14:45:08+02:00 (5 years ago)
Author:
nicolasmartin
Message:

Various corrections on chapters

Cleaning the indexes by fixing/removing wrong entries (or appending a ? to unknown items) and
improve the classification with new index definitions for CPP keys and namelist blocks:

  • from \key{...} cmd, key_ prefix no longer precedes the index entry
  • namelist block declaration moves from \ngn{nam...} to \nam{...} (i.e. \ngn{namtra\_ldf} -> \nam{tra\_ldf}) The expected prefix nam is added to the printed word but not the index entry.

Now we have indexes with a better sorting instead of all CPP keys under 'K' and namelists blocks under 'N'.

Fix missing space issues with alias commands by adding a trailing backslash (\NEMO\, \ie\, \eg\, ...).
There is no perfect solution for this, and I prefer not using a particular package to solve it.

Review the initial LaTeX code snippet for the historic changes in chapters

Finally, for readability and future diff visualisations, please avoid writing paragraphs with continuous lines.
Break the lines around 80 to 100 characters long

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex

    r11386 r11435  
    99\label{chap:LDF} 
    1010 
    11 \minitoc 
     11\chaptertoc 
    1212 
    1313\newpage 
     
    2222(3) the space and time variations of the eddy coefficients. 
    2323These three aspects of the lateral diffusion are set through namelist parameters 
    24 (see the \ngn{nam\_traldf} and \ngn{nam\_dynldf} below). 
     24(see the \nam{tra\_ldf} and \nam{dyn\_ldf} below). 
    2525Note that this chapter describes the standard implementation of iso-neutral tracer mixing.  
    2626Griffies's implementation, which is used if \np{ln\_traldf\_triad}\forcode{ = .true.}, 
     
    5353{Laplacian mixing (\protect\np{ln\_traldf\_lap}, \protect\np{ln\_dynldf\_lap})} 
    5454Setting \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 combine  
     55a second order diffusion on tracers and momentum respectively. Note that in \NEMO\ 4, one can not combine  
    5656Laplacian and Bilaplacian operators for the same variable. 
    5757 
     
    6060Setting \protect\np{ln\_traldf\_blp}\forcode{ = .true.} and/or \protect\np{ln\_dynldf\_blp}\forcode{ = .true.} enables  
    6161a 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. 
     62We stress again that from \NEMO\ 4, the simultaneous use Laplacian and Bilaplacian operators is not allowed. 
    6363 
    6464% ================================================================ 
     
    9090\subsection{Slopes for tracer geopotential mixing in the $s$-coordinate} 
    9191 
    92 In $s$-coordinates, geopotential mixing (\ie horizontal mixing) $r_1$ and $r_2$ are the slopes between 
     92In $s$-coordinates, geopotential mixing (\ie\ horizontal mixing) $r_1$ and $r_2$ are the slopes between 
    9393the geopotential and computational surfaces. 
    9494Their discrete formulation is found by locally solving \autoref{eq:tra_ldf_iso} when 
    9595the 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.  
    9797%gm { Steven : My version is obviously wrong since I'm left with an arbitrary constant which is the local vertical temperature gradient} 
    9898 
     
    124124Their formulation does not depend on the vertical coordinate used. 
    125125Their discrete formulation is found using the fact that the diffusive fluxes of 
    126 locally referenced potential density (\ie $in situ$ density) vanish. 
     126locally referenced potential density (\ie\ $in situ$ density) vanish. 
    127127So, substituting $T$ by $\rho$ in \autoref{eq:tra_ldf_iso} and setting the diffusive fluxes in 
    128128the three directions to zero leads to the following definition for the neutral slopes: 
     
    230230To overcome this problem, several techniques have been proposed in which the numerical schemes of 
    231231the 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}. 
     232Griffies's scheme is now available in \NEMO\ if \np{ln\_traldf\_triad}\forcode{ = .true.}; see \autoref{apdx:triad}. 
    233233Here, another strategy is presented \citep{lazar_phd97}: 
    234234a local filtering of the iso-neutral slopes (made on 9 grid-points) prevents the development of 
     
    284284      \textit{(a)} in the real ocean the slope is the iso-neutral slope in the ocean interior, 
    285285      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: 
    287287      wall boundary condition). 
    288288      Nevertheless, the profile between the surface zero value and the interior iso-neutral one is unknown, 
     
    309309\textit{vw}- points for the $v$-component. 
    310310They 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}: 
    312312 
    313313\[ 
     
    323323The major issue remaining is in the specification of the boundary conditions. 
    324324The 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 bottom 
     325\ie\ using the shear computed along the model levels and with no additional friction at the ocean bottom 
    326326(see \autoref{sec:LBC_coast}). 
    327327 
     
    420420 
    421421The 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 of 
     422\ie\ a hyperbolic tangent variation with depth associated with a grid size dependence of 
    423423the magnitude of the coefficient.  
    424424 
     
    527527the formulation of which depends on the slopes of iso-neutral surfaces. 
    528528Contrary 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, 
    530530and the sum \autoref{eq:ldfslp_geo} + \autoref{eq:ldfslp_iso} in $s$-coordinates. 
    531531 
    532 If isopycnal mixing is used in the standard way, \ie \np{ln\_traldf\_triad}\forcode{ = .false.}, the eddy induced velocity is given by:  
     532If isopycnal mixing is used in the standard way, \ie\ \np{ln\_traldf\_triad}\forcode{ = .false.}, the eddy induced velocity is given by:  
    533533\begin{equation} 
    534534  \label{eq:ldfeiv} 
     
    539539  \end{split} 
    540540\end{equation} 
    541 where $A^{eiv}$ is the eddy induced velocity coefficient whose value is set through \np{nn\_aei\_ijk\_t} \ngn{namtra\_eiv} namelist parameter.  
     541where $A^{eiv}$ is the eddy induced velocity coefficient whose value is set through \np{nn\_aei\_ijk\_t} \nam{tra\_eiv} namelist parameter.  
    542542The three components of the eddy induced velocity are computed in \rou{ldf\_eiv\_trp} and 
    543543added to the eulerian velocity in \rou{tra\_adv} where tracer advection is performed. 
     
    570570%-------------------------------------------------------------------------------------------------------------- 
    571571 
    572 If  \np{ln\_mle}\forcode{ = .true.} in \ngn{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. 
     572If  \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. 
    573573 
    574574\colorbox{yellow}{TBC} 
Note: See TracChangeset for help on using the changeset viewer.