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 11597 for NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex – NEMO

Ignore:
Timestamp:
2019-09-25T20:20:19+02:00 (5 years ago)
Author:
nicolasmartin
Message:

Continuation of coding rules application
Recovery of some sections deleted by the previous commit

File:
1 edited

Legend:

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

    r11596 r11597  
    1717\chaptertoc 
    1818 
     19%% ================================================================================================= 
    1920\section[Choice of \forcode{namtra\_ldf} namelist parameters]{Choice of \protect\nam{tra_ldf}{tra\_ldf} namelist parameters} 
    20 %-----------------------------------------nam_traldf------------------------------------------------------ 
    21  
    22 %--------------------------------------------------------------------------------------------------------- 
     21 
    2322 
    2423Two scheme are available to perform the iso-neutral diffusion. 
     
    6362\end{description} 
    6463 
     64%% ================================================================================================= 
    6565\section{Triad formulation of iso-neutral diffusion} 
    6666\label{sec:TRIADS_iso} 
     
    6969but formulated within the \NEMO\ framework, using scale factors rather than grid-sizes. 
    7070 
     71%% ================================================================================================= 
    7172\subsection{Iso-neutral diffusion operator} 
    7273 
     
    156157is evaluated at $w$-points but involves horizontal gradients defined at $u$-points. 
    157158 
     159%% ================================================================================================= 
    158160\subsection{Standard discretization} 
    159161 
     
    187189(\ie\ they enter the computation of density), but it does not work for a passive tracer. 
    188190 
     191%% ================================================================================================= 
    189192\subsection{Expression of the skew-flux in terms of triad slopes} 
    190193 
     
    282285and in \autoref{eq:TRIADS_i31} $a'_{1}={\:}_i^k{\mathbb{A}_w}_{1/2}^{1/2}$. 
    283286 
     287%% ================================================================================================= 
    284288\subsection{Full triad fluxes} 
    285289 
     
    365369\end{flalign} 
    366370 
     371%% ================================================================================================= 
    367372\subsection{Ensuring the scheme does not increase tracer variance} 
    368373\label{subsec:TRIADS_variance} 
     
    462467\] 
    463468 
     469%% ================================================================================================= 
    464470\subsection{Triad volumes in Griffes's scheme and in \NEMO} 
    465471 
     
    490496we can replace $\overline{A}_{\,i+1/2}^k$ by $A_{i+1/2}^k$ in the above. 
    491497 
     498%% ================================================================================================= 
    492499\subsection{Summary of the scheme} 
    493500 
     
    625632\end{description} 
    626633 
     634%% ================================================================================================= 
    627635\subsection{Treatment of the triads at the boundaries} 
    628636\label{sec:TRIADS_iso_bdry} 
     
    673681% >>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    674682 
     683%% ================================================================================================= 
    675684\subsection{ Limiting of the slopes within the interior} 
    676685\label{sec:TRIADS_limit} 
     
    701710and so acts to reduce gravitational potential energy. 
    702711 
     712%% ================================================================================================= 
    703713\subsection{Tapering within the surface mixed layer} 
    704714\label{sec:TRIADS_taper} 
     
    707717When the Griffies triads are used, we offer two options for this. 
    708718 
     719%% ================================================================================================= 
    709720\subsubsection{Linear slope tapering within the surface mixed layer} 
    710721\label{sec:TRIADS_lintaper} 
     
    822833% >>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    823834 
     835%% ================================================================================================= 
    824836\subsubsection{Additional truncation of skew iso-neutral flux components} 
    825837\label{subsec:TRIADS_Gerdes-taper} 
     
    864876% This may give strange looking results, 
    865877% particularly where the mixed-layer depth varies strongly laterally. 
     878%% ================================================================================================= 
    866879\section{Eddy induced advection formulated as a skew flux} 
    867880\label{sec:TRIADS_skew-flux} 
    868881 
     882%% ================================================================================================= 
    869883\subsection{Continuous skew flux formulation} 
    870884\label{sec:TRIADS_continuous-skew-flux} 
     
    975989Since it has the same divergence as the advective form it also preserves the tracer variance. 
    976990 
     991%% ================================================================================================= 
    977992\subsection{Discrete skew flux formulation} 
    978993 
     
    10171032It also ensures the following two key properties. 
    10181033 
     1034%% ================================================================================================= 
    10191035\subsubsection{No change in tracer variance} 
    10201036 
     
    10391055Hence the two fluxes associated with each triad make no net contribution to the variance budget. 
    10401056 
     1057%% ================================================================================================= 
    10411058\subsubsection{Reduction in gravitational PE} 
    10421059 
     
    10931110\beta_i^k\delta_{k+ k_p}[S^i]<0$, this PE change is negative. 
    10941111 
     1112%% ================================================================================================= 
    10951113\subsection{Treatment of the triads at the boundaries} 
    10961114\label{sec:TRIADS_skew_bdry} 
     
    11041122The namelist parameter \np{ln_botmix_triad}{ln\_botmix\_triad} has no effect on the eddy-induced skew-fluxes. 
    11051123 
     1124%% ================================================================================================= 
    11061125\subsection{Limiting of the slopes within the interior} 
    11071126\label{sec:TRIADS_limitskew} 
     
    11111130Each individual triad \rtriadt{R} is so limited. 
    11121131 
     1132%% ================================================================================================= 
    11131133\subsection{Tapering within the surface mixed layer} 
    11141134\label{sec:TRIADS_taperskew} 
     
    11311151(the horizontal flux convergence is relatively insignificant within the mixed-layer). 
    11321152 
     1153%% ================================================================================================= 
    11331154\subsection{Streamfunction diagnostics} 
    11341155\label{sec:TRIADS_sfdiag} 
Note: See TracChangeset for help on using the changeset viewer.