Changeset 11597 for NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex
- Timestamp:
- 2019-09-25T20:20:19+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.tex
r11596 r11597 17 17 \chaptertoc 18 18 19 %% ================================================================================================= 19 20 \section[Choice of \forcode{namtra\_ldf} namelist parameters]{Choice of \protect\nam{tra_ldf}{tra\_ldf} namelist parameters} 20 %-----------------------------------------nam_traldf------------------------------------------------------ 21 22 %--------------------------------------------------------------------------------------------------------- 21 23 22 24 23 Two scheme are available to perform the iso-neutral diffusion. … … 63 62 \end{description} 64 63 64 %% ================================================================================================= 65 65 \section{Triad formulation of iso-neutral diffusion} 66 66 \label{sec:TRIADS_iso} … … 69 69 but formulated within the \NEMO\ framework, using scale factors rather than grid-sizes. 70 70 71 %% ================================================================================================= 71 72 \subsection{Iso-neutral diffusion operator} 72 73 … … 156 157 is evaluated at $w$-points but involves horizontal gradients defined at $u$-points. 157 158 159 %% ================================================================================================= 158 160 \subsection{Standard discretization} 159 161 … … 187 189 (\ie\ they enter the computation of density), but it does not work for a passive tracer. 188 190 191 %% ================================================================================================= 189 192 \subsection{Expression of the skew-flux in terms of triad slopes} 190 193 … … 282 285 and in \autoref{eq:TRIADS_i31} $a'_{1}={\:}_i^k{\mathbb{A}_w}_{1/2}^{1/2}$. 283 286 287 %% ================================================================================================= 284 288 \subsection{Full triad fluxes} 285 289 … … 365 369 \end{flalign} 366 370 371 %% ================================================================================================= 367 372 \subsection{Ensuring the scheme does not increase tracer variance} 368 373 \label{subsec:TRIADS_variance} … … 462 467 \] 463 468 469 %% ================================================================================================= 464 470 \subsection{Triad volumes in Griffes's scheme and in \NEMO} 465 471 … … 490 496 we can replace $\overline{A}_{\,i+1/2}^k$ by $A_{i+1/2}^k$ in the above. 491 497 498 %% ================================================================================================= 492 499 \subsection{Summary of the scheme} 493 500 … … 625 632 \end{description} 626 633 634 %% ================================================================================================= 627 635 \subsection{Treatment of the triads at the boundaries} 628 636 \label{sec:TRIADS_iso_bdry} … … 673 681 % >>>>>>>>>>>>>>>>>>>>>>>>>>>> 674 682 683 %% ================================================================================================= 675 684 \subsection{ Limiting of the slopes within the interior} 676 685 \label{sec:TRIADS_limit} … … 701 710 and so acts to reduce gravitational potential energy. 702 711 712 %% ================================================================================================= 703 713 \subsection{Tapering within the surface mixed layer} 704 714 \label{sec:TRIADS_taper} … … 707 717 When the Griffies triads are used, we offer two options for this. 708 718 719 %% ================================================================================================= 709 720 \subsubsection{Linear slope tapering within the surface mixed layer} 710 721 \label{sec:TRIADS_lintaper} … … 822 833 % >>>>>>>>>>>>>>>>>>>>>>>>>>>> 823 834 835 %% ================================================================================================= 824 836 \subsubsection{Additional truncation of skew iso-neutral flux components} 825 837 \label{subsec:TRIADS_Gerdes-taper} … … 864 876 % This may give strange looking results, 865 877 % particularly where the mixed-layer depth varies strongly laterally. 878 %% ================================================================================================= 866 879 \section{Eddy induced advection formulated as a skew flux} 867 880 \label{sec:TRIADS_skew-flux} 868 881 882 %% ================================================================================================= 869 883 \subsection{Continuous skew flux formulation} 870 884 \label{sec:TRIADS_continuous-skew-flux} … … 975 989 Since it has the same divergence as the advective form it also preserves the tracer variance. 976 990 991 %% ================================================================================================= 977 992 \subsection{Discrete skew flux formulation} 978 993 … … 1017 1032 It also ensures the following two key properties. 1018 1033 1034 %% ================================================================================================= 1019 1035 \subsubsection{No change in tracer variance} 1020 1036 … … 1039 1055 Hence the two fluxes associated with each triad make no net contribution to the variance budget. 1040 1056 1057 %% ================================================================================================= 1041 1058 \subsubsection{Reduction in gravitational PE} 1042 1059 … … 1093 1110 \beta_i^k\delta_{k+ k_p}[S^i]<0$, this PE change is negative. 1094 1111 1112 %% ================================================================================================= 1095 1113 \subsection{Treatment of the triads at the boundaries} 1096 1114 \label{sec:TRIADS_skew_bdry} … … 1104 1122 The namelist parameter \np{ln_botmix_triad}{ln\_botmix\_triad} has no effect on the eddy-induced skew-fluxes. 1105 1123 1124 %% ================================================================================================= 1106 1125 \subsection{Limiting of the slopes within the interior} 1107 1126 \label{sec:TRIADS_limitskew} … … 1111 1130 Each individual triad \rtriadt{R} is so limited. 1112 1131 1132 %% ================================================================================================= 1113 1133 \subsection{Tapering within the surface mixed layer} 1114 1134 \label{sec:TRIADS_taperskew} … … 1131 1151 (the horizontal flux convergence is relatively insignificant within the mixed-layer). 1132 1152 1153 %% ================================================================================================= 1133 1154 \subsection{Streamfunction diagnostics} 1134 1155 \label{sec:TRIADS_sfdiag}
Note: See TracChangeset
for help on using the changeset viewer.