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 9392 for branches/2017/dev_merge_2017/DOC/tex_sub/annex_iso.tex – NEMO

Ignore:
Timestamp:
2018-03-09T16:57:00+01:00 (6 years ago)
Author:
nicolasmartin
Message:

Global replacement of patterns \np{id}=value by \forcode{id = value} for integer and booleans

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/DOC/tex_sub/annex_iso.tex

    r9389 r9392  
    1515 
    1616Two scheme are available to perform the iso-neutral diffusion.  
    17 If the namelist logical \np{ln\_traldf\_triad} is set true,  
     17If the namelist logical \np{ln_traldf_triad} is set true,  
    1818\NEMO updates both active and passive tracers using the Griffies triad representation  
    1919of iso-neutral diffusion and the eddy-induced advective skew (GM) fluxes.  
    20 If the namelist logical \np{ln\_traldf\_iso} is set true,  
     20If the namelist logical \np{ln_traldf_iso} is set true,  
    2121the filtered version of Cox's original scheme (the Standard scheme) is employed (\S\ref{LDF_slp}).  
    2222In the present implementation of the Griffies scheme,  
    23 the advective skew fluxes are implemented even if \np{ln\_traldf\_eiv} is false. 
     23the advective skew fluxes are implemented even if \np{ln_traldf_eiv} is false. 
    2424 
    2525Values of iso-neutral diffusivity and GM coefficient are set as 
     
    3131The options specific to the Griffies scheme include: 
    3232\begin{description}[font=\normalfont] 
    33 \item[\np{ln\_triad\_iso}] See \S\ref{sec:triad:taper}. If this is set false (the default), then 
     33\item[\np{ln_triad_iso}] See \S\ref{sec:triad:taper}. If this is set false (the default), then 
    3434  `iso-neutral' mixing is accomplished within the surface mixed-layer 
    3535  along slopes linearly decreasing with depth from the value immediately below 
    3636  the mixed-layer to zero (flat) at the surface (\S\ref{sec:triad:lintaper}).  
    3737  This is the same treatment as used in the default implementation \S\ref{LDF_slp_iso}; Fig.~\ref{Fig_eiv_slp}.   
    38   Where \np{ln\_triad\_iso} is set true, the vertical skew flux is further reduced  
     38  Where \np{ln_triad_iso} is set true, the vertical skew flux is further reduced  
    3939  to ensure no vertical buoyancy flux, giving an almost pure 
    4040  horizontal diffusive tracer flux within the mixed layer. This is similar to 
    4141  the tapering suggested by \citet{Gerdes1991}. See \S\ref{sec:triad:Gerdes-taper} 
    42 \item[\np{ln\_botmix\_triad}] See \S\ref{sec:triad:iso_bdry}.  
     42\item[\np{ln_botmix_triad}] See \S\ref{sec:triad:iso_bdry}.  
    4343  If this is set false (the default) then the lateral diffusive fluxes 
    4444  associated with triads partly masked by topography are neglected.  
    4545  If it is set true, however, then these lateral diffusive fluxes are applied,  
    4646  giving smoother bottom tracer fields at the cost of introducing diapycnal mixing. 
    47 \item[\np{rn\_sw\_triad}]  blah blah to be added.... 
     47\item[\np{rn_sw_triad}]  blah blah to be added.... 
    4848\end{description} 
    4949The options shared with the Standard scheme include: 
    5050\begin{description}[font=\normalfont] 
    51 \item[\np{ln\_traldf\_msc}]   blah blah to be added 
    52 \item[\np{rn\_slpmax}]  blah blah to be added 
     51\item[\np{ln_traldf_msc}]   blah blah to be added 
     52\item[\np{rn_slpmax}]  blah blah to be added 
    5353\end{description} 
    5454\section{Triad formulation of iso-neutral diffusion} 
     
    651651or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ $u$-point is 
    652652masked. The associated lateral fluxes (grey-black dashed line) are 
    653 masked if \np{ln\_botmix\_triad}=false, but left unmasked, 
    654 giving bottom mixing, if \np{ln\_botmix\_triad}=true. 
    655  
    656 The default option \np{ln\_botmix\_triad}=false is suitable when the 
    657 bbl mixing option is enabled (\key{trabbl}, with \np{nn\_bbl\_ldf}=1), 
     653masked if \forcode{ln_botmix_triad = .false.}, but left unmasked, 
     654giving bottom mixing, if \forcode{ln_botmix_triad = .true.}. 
     655 
     656The default option \forcode{ln_botmix_triad = .false.} is suitable when the 
     657bbl mixing option is enabled (\key{trabbl}, with \forcode{nn_bbl_ldf = 1}), 
    658658or  for simple idealized  problems. For setups with topography without 
    659 bbl mixing, \np{ln\_botmix\_triad}=true may be necessary. 
     659bbl mixing, \forcode{ln_botmix_triad = .true.} may be necessary. 
    660660% >>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    661661\begin{figure}[h] \begin{center} 
     
    674674      or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ $u$-point 
    675675      is masked. The associated lateral fluxes (grey-black dashed 
    676       line) are masked if \protect\np{botmix\_triad}=.false., but left 
    677       unmasked, giving bottom mixing, if \protect\np{botmix\_triad}=.true.} 
     676      line) are masked if \protect\np{botmix_triad}=.false., but left 
     677      unmasked, giving bottom mixing, if \protect\np{botmix_triad}=.true.} 
    678678 \end{center} \end{figure} 
    679679% >>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    710710\subsubsection{Linear slope tapering within the surface mixed layer}\label{sec:triad:lintaper} 
    711711This is the option activated by the default choice 
    712 \np{ln\_triad\_iso}=false. Slopes $\tilde{r}_i$ relative to 
     712\forcode{ln_triad_iso = .false.}. Slopes $\tilde{r}_i$ relative to 
    713713geopotentials are tapered linearly from their value immediately below the mixed layer to zero at the 
    714714surface, as described in option (c) of Fig.~\ref{Fig_eiv_slp}, to values 
     
    842842  components} 
    843843\label{sec:triad:Gerdes-taper} 
    844 The alternative option is activated by setting \np{ln\_triad\_iso} = 
     844The alternative option is activated by setting \np{ln_triad_iso} = 
    845845  true. This retains the same tapered slope $\rML$  described above for the 
    846846calculation of the $_{33}$ term of the iso-neutral diffusion tensor (the 
     
    917917it to the Eulerian velocity prior to computing the tracer 
    918918advection. This is implemented if \key{traldf\_eiv} is set in the 
    919 default implementation, where \np{ln\_traldf\_triad} is set 
     919default implementation, where \np{ln_traldf_triad} is set 
    920920false. This allows us to take advantage of all the advection schemes 
    921921offered for the tracers (see \S\ref{TRA_adv}) and not just a $2^{nd}$ 
     
    924924paramount importance. 
    925925 
    926 However, when \np{ln\_traldf\_triad} is set true, \NEMO instead 
     926However, when \np{ln_traldf_triad} is set true, \NEMO instead 
    927927implements eddy induced advection according to the so-called skew form 
    928928\citep{Griffies_JPO98}. It is based on a transformation of the advective fluxes 
     
    11231123and $\triadt{i+1}{k}{R}{-1/2}{1/2}$ are masked when either of the 
    11241124$i,k+1$ or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ 
    1125 $u$-point is masked. The namelist parameter \np{ln\_botmix\_triad} has 
     1125$u$-point is masked. The namelist parameter \np{ln_botmix_triad} has 
    11261126no effect on the eddy-induced skew-fluxes. 
    11271127 
     
    11381138option (c) of Fig.~\ref{Fig_eiv_slp}. This linear tapering for the 
    11391139slopes used to calculate the eddy-induced fluxes is 
    1140 unaffected by the value of \np{ln\_triad\_iso}. 
     1140unaffected by the value of \np{ln_triad_iso}. 
    11411141 
    11421142The justification for this linear slope tapering is that, for $A_e$ 
     
    11531153 
    11541154\subsection{Streamfunction diagnostics}\label{sec:triad:sfdiag} 
    1155 Where the namelist parameter \np{ln\_traldf\_gdia}=true, diagnosed 
     1155Where the namelist parameter \forcode{ln_traldf_gdia = .true.}, diagnosed 
    11561156mean eddy-induced velocities are output. Each time step, 
    11571157streamfunctions are calculated in the $i$-$k$ and $j$-$k$ planes at 
Note: See TracChangeset for help on using the changeset viewer.