Changeset 3267
- Timestamp:
- 2012-01-18T11:28:11+01:00 (12 years ago)
- Location:
- branches/2011/dev_NEMO_MERGE_2011/DOC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/DOC/NEMO_book.tex
r3236 r3267 17 17 \usepackage{xcolor} 18 18 %\usepackage{graphics} % allows insertion of pictures 19 \usepackage{graphicx} % allows insertion of 20 % pictures21 \ usepackage[capbesideposition={top,center}]{floatrow}22 \ floatsetup[table]{style=plaintop}23 \usepackage [margin=10pt,font={small},labelsep=colon,labelfont={bf}]{caption}19 \usepackage{graphicx} % allows insertion of pictures 20 \usepackage[capbesideposition={top,center}]{floatrow} % allows captions 21 \floatsetup[table]{style=plaintop} % beside pictures 22 \usepackage[margin=10pt,font={small},labelsep=colon,labelfont={bf}]{caption} % Gives small font for captions 23 \usepackage{enumitem} % allows non-bold description items 24 24 %\usepackage{colortbl} % gives coloured panels behind table columns 25 25 … … 194 194 % ================================================================ 195 195 196 \usepackage{pstricks}196 %\usepackage{pstricks} 197 197 \title{ 198 \psset{unit=1.1in,linewidth=4pt} %parameters of the units for pstricks 199 %\rput(-1.4,2){ \includegraphics[width=0.2\textwidth]{./TexFiles/Figures/logo_CNRS.pdf} } 200 %\quad 201 %\rput(-.15,2){ \includegraphics[width=0.2\textwidth]{./TexFiles/Figures/logo_MERCATOR.pdf} } 202 %\quad 203 %\rput(.85,2){ \includegraphics[width=0.1\textwidth]{./TexFiles/Figures/logo_UKMO.pdf} } 204 %\quad 205 %\rput(1.85,2){ \includegraphics[width=0.2\textwidth]{./TexFiles/Figures/logo_NERC.pdf} } \\ 206 %\begin{figure}[t] 207 % \includegraphics[width=0.9\textwidth]{./TexFiles/Figures/logo_ALL.pdf} 208 %\end{figure} 209 \rput(0,2){ \includegraphics[width=1.1\textwidth]{./TexFiles/Figures/logo_ALL.pdf} } \\ 210 \vspace{0.1cm} 198 %\psset{unit=1.1in,linewidth=4pt} %parameters of the units for pstricks 199 % \rput(0,2){ \includegraphics[width=1.1\textwidth]{./TexFiles/Figures/logo_ALL.pdf} } \\ 200 % \vspace{0.1cm} 201 \vspace{-6.0cm} 202 \includegraphics[width=1.1\textwidth]{./TexFiles/Figures/logo_ALL.pdf}\\ 203 \vspace{5.1cm} 211 204 \includegraphics[width=0.9\textwidth]{./TexFiles/Figures/NEMO_logo_Black.pdf} \\ 212 205 \vspace{1.4cm} -
branches/2011/dev_NEMO_MERGE_2011/DOC/TexFiles/Chapters/Annex_ISO.tex
r3218 r3267 2 2 % Iso-neutral diffusion : 3 3 % ================================================================ 4 \chapter{Griffies's iso-neutral diffusion} 4 \chapter[Iso-neutral diffusion and eddy advection using 5 triads]{Iso-neutral diffusion and eddy advection using triads} 5 6 \label{sec:triad} 6 7 \minitoc 7 8 \section{Griffies's formulation of iso-neutral diffusion} 8 \pagebreak 9 \section{Choice of namelist parameters} 10 %-----------------------------------------nam_traldf------------------------------------------------------ 11 \namdisplay{namtra_ldf} 12 %--------------------------------------------------------------------------------------------------------- 13 If the namelist variable \np{ln\_traldf\_grif} is set true (and 14 \key{ldfslp} is set), \NEMO updates both active and passive tracers 15 using the Griffies triad representation of iso-neutral diffusion and 16 the eddy-induced advective skew (GM) fluxes. Otherwise (by default) the 17 filtered version of Cox's original scheme is employed 18 (\S\ref{LDF_slp}). In the present implementation of the Griffies 19 scheme, the advective skew fluxes are implemented even if 20 \key{traldf\_eiv} is not set. 21 22 Values of iso-neutral diffusivity and GM coefficient are set as 23 described in \S\ref{LDF_coef}. If none of the keys \key{traldf\_cNd}, 24 N=1,2,3 is set (the default), spatially constant iso-neutral $A_l$ and 25 GM diffusivity $A_e$ are directly set by \np{rn\_aeih\_0} and 26 \np{rn\_aeiv\_0}. If 2D-varying coefficients are set with 27 \key{traldf\_c2d} then $A_l$ is reduced in proportion with horizontal 28 scale factor according to \eqref{Eq_title} \footnote{Except in global 29 $0.5^{\circ}$ runs (\key{orca\_r05}) with \key{traldf\_eiv}, where 30 $A_l$ is set like $A_e$ but with a minimum vale of 31 $100\;\mathrm{m}^2\;\mathrm{s}^{-1}$}. In idealised setups with 32 \key{traldf\_c2d}, $A_e$ is reduced similarly, but if \key{traldf\_eiv} 33 is set in the global configurations \key{orca\_r2}, \key{orca\_r1} or 34 \key{orca\_r05} with \key{traldf\_c2d}, a horizontally varying $A_e$ is 35 instead set from the Held-Larichev parameterisation\footnote{In this 36 case, $A_e$ at low latitudes $|\theta|<20^{\circ}$ is further 37 reduced by a factor $|f/f_{20}|$, where $f_{20}$ is the value of $f$ 38 at $20^{\circ}$~N} (\mdl{ldfeiv}) and \np{rn\_aeiv\_0} is ignored 39 unless it is zero. 40 41 The options specific to the Griffies scheme include: 42 \begin{description}[font=\normalfont] 43 \item[\np{ln\_traldf\_gdia}] Default value is false. See \S\ref{sec:triad:sfdiag}. If this is set true, time-mean 44 eddy-advective (GM) velocities are output for diagnostic purposes, even 45 though the eddy advection is accomplished by means of the skew 46 fluxes. 47 \item[\np{ln\_traldf\_iso}] See \S\ref{sec:triad:taper}. If this is set false (the default), then 48 `iso-neutral' mixing is accomplished within the surface mixed-layer 49 along slopes linearly decreasing with depth from the value immediately below 50 the mixed-layer to zero (flat) at the surface (\S\ref{sec:triad:lintaper}). This is the same 51 treatment as used in the default implementation 52 \S\ref{LDF_slp_iso}; Fig.~\ref{Fig_eiv_slp}. Where 53 \np{ln\_traldf\_iso} is set true, the vertical skew flux is further 54 reduced to ensure no vertical buoyancy flux, giving an almost pure 55 horizontal diffusive tracer flux within the mixed layer. This is similar to 56 the tapering suggested by \citet{Gerdes1991}. See \S\ref{sec:triad:Gerdes-taper} 57 \item[\np{ln\_traldf\_botmix}] See \S\ref{sec:triad:iso_bdry}. If this 58 is set false (the default) then the lateral diffusive fluxes 59 associated with triads partly masked by topography are neglected. If 60 it is set true, however, then these lateral diffusive fluxes are 61 applied, giving smoother bottom tracer fields at the cost of 62 introducing diapycnal mixing. 63 \end{description} 64 \section{Triad formulation of iso-neutral diffusion} 9 65 \label{sec:triad:iso} 10 11 We define a scheme inspired by \citet{Griffies_al_JPO98}, but formulated within the \NEMO 66 We have implemented into \NEMO a scheme inspired by \citet{Griffies_al_JPO98}, but formulated within the \NEMO 12 67 framework, using scale factors rather than grid-sizes. 13 68 … … 613 668 or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ $u$-point is 614 669 masked. The associated lateral fluxes (grey-black dashed line) are 615 masked if \n lv{ln\_botmix\_grif=.false.}, but left unmasked,616 giving bottom mixing, if \n lv{ln\_botmix\_grif=.true.}.617 618 The default option \n lv{ln\_botmix\_grif=.false.}is suitable when the619 bbl mixing option is enabled (\key{trabbl}, with \n lv{nn\_bbl\_ldf=1}),670 masked if \np{ln\_botmix\_grif}=false, but left unmasked, 671 giving bottom mixing, if \np{ln\_botmix\_grif}=true. 672 673 The default option \np{ln\_botmix\_grif}=false is suitable when the 674 bbl mixing option is enabled (\key{trabbl}, with \np{nn\_bbl\_ldf}=1), 620 675 or for simple idealized problems. For setups with topography without 621 bbl mixing, \n lv{ln\_botmix\_grif=.true.}may be necessary.676 bbl mixing, \np{ln\_botmix\_grif}=true may be necessary. 622 677 % >>>>>>>>>>>>>>>>>>>>>>>>>>>> 623 678 \begin{figure}[h] \begin{center} … … 636 691 or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ $u$-point 637 692 is masked. The associated lateral fluxes (grey-black dashed 638 line) are masked if \ smnlv{ln\_botmix\_grif=.false.}, but left639 unmasked, giving bottom mixing, if \ smnlv{ln\_botmix\_grif=.true.}}693 line) are masked if \np{botmix\_grif}=.false., but left 694 unmasked, giving bottom mixing, if \np{botmix\_grif}=.true.} 640 695 \end{center} \end{figure} 641 696 % >>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 665 720 iso-neutral density flux that drives dianeutral mixing. In particular this iso-neutral density flux 666 721 is always downwards, and so acts to reduce gravitational potential energy. 667 \subsection{Tapering within the surface mixed layer} 722 \subsection{Tapering within the surface mixed layer}\label{sec:triad:taper} 723 668 724 Additional tapering of the iso-neutral fluxes is necessary within the 669 725 surface mixed layer. When the Griffies triads are used, we offer two … … 671 727 \subsubsection{Linear slope tapering within the surface mixed layer}\label{sec:triad:lintaper} 672 728 This is the option activated by the default choice 673 \n lv{ln\_triad\_iso=.false.}. Slopes $\tilde{r}_i$ relative to729 \np{ln\_triad\_iso}=false. Slopes $\tilde{r}_i$ relative to 674 730 geopotentials are tapered linearly from their value immediately below the mixed layer to zero at the 675 731 surface, as described in option (c) of Fig.~\ref{Fig_eiv_slp}, to values … … 798 854 % >>>>>>>>>>>>>>>>>>>>>>>>>>>> 799 855 800 \subsubsection{Additional truncation of skew iso-neutral flux components} 801 The alternative option is activated by setting \nlv{ln\_triad\_iso = 802 .true.}. This retains the same tapered slope $\rML$ described above for the 856 \subsubsection{Additional truncation of skew iso-neutral flux 857 components} 858 \label{sec:triad:Gerdes-taper} 859 The alternative option is activated by setting \np{ln\_triad\_iso} = 860 true. This retains the same tapered slope $\rML$ described above for the 803 861 calculation of the $_{33}$ term of the iso-neutral diffusion tensor (the 804 862 vertical tracer flux driven by vertical tracer gradients), but … … 839 897 % Skew flux formulation for Eddy Induced Velocity : 840 898 % ================================================================ 841 \section{Eddy induced advection and its formulation as a skewflux}899 \section{Eddy induced advection formulated as a skew flux}\label{sec:triad:skew-flux} 842 900 843 901 \subsection{The continuous skew flux formulation}\label{sec:triad:continuous-skew-flux} 844 902 845 When Gent and McWilliams's [1990] diffusion is used (\key{traldf\_eiv} defined),903 When Gent and McWilliams's [1990] diffusion is used, 846 904 an additional advection term is added. The associated velocity is the so called 847 905 eddy induced velocity, the formulation of which depends on the slopes of iso- … … 852 910 853 911 The eddy induced velocity is given by: 854 \begin{equation} \label{eq:triad:eiv_v} 912 \begin{subequations} \label{eq:triad:eiv} 913 \begin{equation}\label{eq:triad:eiv_v} 855 914 \begin{split} 856 u^* & = - \frac{1}{e_{3}}\; \partial_ k \left( A_{e} \; \tilde{r}_1 \right)\\857 v^* & = - \frac{1}{e_{3}}\; \partial_ k \left( A_{e} \; \tilde{r}_2 \right)\\858 w^* & = \frac{1}{e_{1}e_{2}}\; \left\{ \partial_i \left( e_{2} \, A_{e} \; \tilde{r}_1\right)859 + \partial_j \left( e_{1} \, A_{e} \;\tilde{r}_2 \right) \right\}915 u^* & = - \frac{1}{e_{3}}\; \partial_i\psi_1, \\ 916 v^* & = - \frac{1}{e_{3}}\; \partial_j\psi_2, \\ 917 w^* & = \frac{1}{e_{1}e_{2}}\; \left\{ \partial_i \left( e_{2} \, \psi_1\right) 918 + \partial_j \left( e_{1} \, \psi_2\right) \right\}, 860 919 \end{split} 861 920 \end{equation} 862 where $A_{e}$ is the eddy induced velocity coefficient, and $\tilde{r}_1$ and $\tilde{r}_2$ the slopes between the iso-neutral and the geopotential surfaces. 863 864 The traditional way to implement this additional advection is to add it to the Eulerian 865 velocity prior to computing the tracer advection. This allows us to take advantage of 866 all the advection schemes offered for the tracers (see \S\ref{TRA_adv}) and not just 867 a $2^{nd}$ order advection scheme. This is particularly useful for passive tracers 868 where \emph{positivity} of the advection scheme is of paramount importance. 869 870 \citet{Griffies_JPO98} introduces another way to implement the eddy induced advection, 871 the so-called skew form. It is based on a transformation of the advective fluxes 921 where the streamfunctions $\psi_i$ are given by 922 \begin{equation} \label{eq:triad:eiv_psi} 923 \begin{split} 924 \psi_1 & = A_{e} \; \tilde{r}_1, \\ 925 \psi_2 & = A_{e} \; \tilde{r}_2, 926 \end{split} 927 \end{equation} 928 \end{subequations} 929 with $A_{e}$ the eddy induced velocity coefficient, and $\tilde{r}_1$ and $\tilde{r}_2$ the slopes between the iso-neutral and the geopotential surfaces. 930 931 The traditional way to implement this additional advection is to add 932 it to the Eulerian velocity prior to computing the tracer 933 advection. This is implemented if \key{traldf\_eiv} is set in the 934 default implementation, where \np{ln\_traldf\_grif} is set 935 false. This allows us to take advantage of all the advection schemes 936 offered for the tracers (see \S\ref{TRA_adv}) and not just a $2^{nd}$ 937 order advection scheme. This is particularly useful for passive 938 tracers where \emph{positivity} of the advection scheme is of 939 paramount importance. 940 941 However, when \np{ln\_traldf\_grif} is set true, \NEMO instead 942 implements eddy induced advection according to the so-called skew form 943 \citep{Griffies_JPO98}. It is based on a transformation of the advective fluxes 872 944 using the non-divergent nature of the eddy induced velocity. 873 945 For example in the (\textbf{i},\textbf{k}) plane, the tracer advective … … 883 955 &= 884 956 \begin{pmatrix} 885 { - \partial_k \left( e_{2} \, A_{e} \; \tilde{r}_1 \right) \; T \;} \\886 {+ \partial_i \left( e_{2} \, A_{e} \; \tilde{r}_1 \right) \; T \;} \\957 { - \partial_k \left( e_{2} \,\psi_1 \right) \; T \;} \\ 958 {+ \partial_i \left( e_{2} \, \psi_1 \right) \; T \;} \\ 887 959 \end{pmatrix} \\ 888 960 &= 889 961 \begin{pmatrix} 890 { - \partial_k \left( e_{2} \, A_{e} \; \tilde{r}_1 \; T \right) \;} \\891 {+ \partial_i \left( e_{2} \, A_{e} \; \tilde{r}_1\; T \right) \;} \\962 { - \partial_k \left( e_{2} \, \psi_1 \; T \right) \;} \\ 963 {+ \partial_i \left( e_{2} \,\psi_1 \; T \right) \;} \\ 892 964 \end{pmatrix} 893 965 + 894 966 \begin{pmatrix} 895 {+ e_{2} \, A_{e} \; \tilde{r}_1 \; \partial_k T} \\896 { - e_{2} \, A_{e} \; \tilde{r}_1 \; \partial_i T} \\967 {+ e_{2} \, \psi_1 \; \partial_k T} \\ 968 { - e_{2} \, \psi_1 \; \partial_i T} \\ 897 969 \end{pmatrix} 898 970 \end{split} … … 902 974 \begin{equation} \label{eq:triad:eiv_skew_ijk} 903 975 \textbf{F}_\mathrm{eiv}^T = \begin{pmatrix} 904 {+ e_{2} \, A_{e} \; \tilde{r}_1 \; \partial_k T} \\905 { - e_{2} \, A_{e} \; \tilde{r}_1 \; \partial_i T} \\976 {+ e_{2} \, \psi_1 \; \partial_k T} \\ 977 { - e_{2} \, \psi_1 \; \partial_i T} \\ 906 978 \end{pmatrix} 907 979 \end{equation} … … 909 981 \begin{equation}\label{eq:triad:eiv_skew_physical} 910 982 \begin{split} 911 f^*_1 & = \frac{1}{e_{3}}\; A_{e} \; \tilde{r}_1 \partial_k T \\912 f^*_2 & = \frac{1}{e_{3}}\; A_{e} \; \tilde{r}_2 \partial_k T \\913 f^*_3 & = -\frac{1}{e_{1}e_{2}}\; A_{e} \left\{ e_{2} \tilde{r}_1 \partial_i T914 + e_{1} \ tilde{r}_2 \partial_j T \right\}. \\983 f^*_1 & = \frac{1}{e_{3}}\; \psi_1 \partial_k T \\ 984 f^*_2 & = \frac{1}{e_{3}}\; \psi_2 \partial_k T \\ 985 f^*_3 & = -\frac{1}{e_{1}e_{2}}\; \left\{ e_{2} \psi_1 \partial_i T 986 + e_{1} \psi_2 \partial_j T \right\}. \\ 915 987 \end{split} 916 988 \end{equation} 917 989 Note that Eq.~ \eqref{eq:triad:eiv_skew_physical} takes the same form whatever the 918 990 vertical coordinate, though of course the slopes 919 $\tilde{r}_i$ are relative to geopotentials.991 $\tilde{r}_i$ which define the $\psi_i$ in \eqref{eq:triad:eiv_psi} are relative to geopotentials. 920 992 The tendency associated with eddy induced velocity is then simply the convergence 921 993 of the fluxes (\ref{eq:triad:eiv_skew_ijk}, \ref{eq:triad:eiv_skew_physical}), so 922 994 \begin{equation} \label{eq:triad:skew_eiv_conv} 923 995 \frac{\partial T}{\partial t}= -\frac{1}{e_1 \, e_2 \, e_3 } \left[ 924 \frac{\partial}{\partial i} \left( e_2 A_{e} \; \tilde{r}_1 \partial_k T\right)925 + \frac{\partial}{\partial j} \left( e_1 A_{e}\;926 \ tilde{r}_2 \partial_k T\right)927 - \frac{\partial}{\partial k} A_{e} \left( e_{2} \tilde{r}_1 \partial_i T928 + e_{1} \ tilde{r}_2 \partial_j T \right) \right]996 \frac{\partial}{\partial i} \left( e_2 \psi_1 \partial_k T\right) 997 + \frac{\partial}{\partial j} \left( e_1 \; 998 \psi_2 \partial_k T\right) 999 - \frac{\partial}{\partial k} \left( e_{2} \psi_1 \partial_i T 1000 + e_{1} \psi_2 \partial_j T \right) \right] 929 1001 \end{equation} 930 1002 It naturally conserves the tracer content, as it is expressed in flux … … 976 1048 The discretization conserves tracer variance, $i.e.$ it does not 977 1049 include a diffusive component but is a `pure' advection term. This can 978 be seen either from Appendix \ref{Apdx_eiv_skew} or by considering the 1050 be seen 1051 %either from Appendix \ref{Apdx_eiv_skew} or 1052 by considering the 979 1053 fluxes associated with a given triad slope 980 1054 $_i^k{\mathbb{R}}_{i_p}^{k_p} (T)$. For, following … … 1064 1138 and $\triadt{i+1}{k}{R}{-1/2}{1/2}$ are masked when either of the 1065 1139 $i,k+1$ or $i+1,k+1$ tracer points is masked, i.e.\ the $i,k+1$ 1066 $u$-point is masked. The namelist parameter \n lv{ln\_botmix\_grif} has1140 $u$-point is masked. The namelist parameter \np{ln\_botmix\_grif} has 1067 1141 no effect on the eddy-induced skew-fluxes. 1068 1142 … … 1079 1153 option (c) of Fig.~\ref{Fig_eiv_slp}. This linear tapering for the 1080 1154 slopes used to calculate the eddy-induced fluxes is 1081 unaffected by the value of \n lv{ln\_triad\_iso}.1155 unaffected by the value of \np{ln\_triad\_iso}. 1082 1156 1083 1157 The justification for this linear slope tapering is that, for $A_e$ … … 1094 1168 1095 1169 \subsection{Streamfunction diagnostics}\label{sec:triad:sfdiag} 1096 Where the namelist parameter \n lv{ln\_botmix\_grif=.true.}, diagnosed1170 Where the namelist parameter \np{ln\_botmix\_grif}=true, diagnosed 1097 1171 mean eddy-induced velocities are output. Each time step, 1098 1172 streamfunctions are calculated in the $i$-$k$ and $j$-$k$ planes at … … 1104 1178 \begin{equation} 1105 1179 \label{eq:triad:sfdiagi} 1106 {\psi_{[i]}}_{i+1/2}^{k+1/2}={\quarter}\sum_{\substack{i_p,\,k_p}} 1107 {A_e}_{i+1/2-i_p}^{k+1/2-k_p}\:\triadd{i+1/2-i_p}{k+1/2-k_p}{R}{i_p}{k_p} 1108 \end{equation} 1109 1110 \newpage %force an empty line 1111 % ================================================================ 1112 % Discrete Invariants of the skew flux formulation 1113 % ================================================================ 1114 \subsection{Discrete Invariants of the skew flux formulation} 1115 \label{Apdx_eiv_skew} 1116 1117 1118 Demonstration for the conservation of the tracer variance in the (\textbf{i},\textbf{j}) plane. 1119 1120 This have to be moved in an Appendix. 1121 1122 The continuous property to be demonstrated is : 1123 \begin{align*} 1124 \int_D \nabla \cdot \textbf{F}_\mathrm{eiv}(T) \; T \;dv \equiv 0 1125 \end{align*} 1126 The discrete form of its left hand side is obtained using \eqref{eq:triad:allskewflux} 1127 \begin{align*} 1128 \sum\limits_{i,k} \sum_{\substack{i_p,\,k_p}} \Biggl\{ \;\; 1129 \delta_i &\left[ 1130 {e_{2u}}_{i+i_p+1/2}^{k} \;\ \ {A_{e}}_{i+i_p+1/2}^{k} 1131 \ \ \ { _{i+i_p+1/2}^k \mathbb{R}_{-i_p}^{k_p} } \quad \delta_{k+k_p}[T_{i+i_p+1/2}] 1132 \right] \; T_i^k \\ 1133 - \delta_k &\left[ 1134 {e_{2u}}_i^{k+k_p+1/2} \ \ {A_{e}}_i^{k+k_p+1/2} 1135 \ \ { _i^{k+k_p+1/2} \mathbb{R}_{i_p}^{-k_p} } \ \ \delta_{i+i_p}[T^{k+k_p+1/2}] 1136 \right] \; T_i^k \ \Biggr\} 1137 \end{align*} 1138 apply the adjoint of delta operator, it becomes 1139 \begin{align*} 1140 \sum\limits_{i,k} \sum_{\substack{i_p,\,k_p}} \Biggl\{ \;\; 1141 &\left( 1142 {e_{2u}}_{i+i_p+1/2}^{k} \;\ \ {A_{e}}_{i+i_p+1/2}^{k} 1143 \ \ \ { _{i+i_p+1/2}^k \mathbb{R}_{-i_p}^{k_p} } \quad \delta_{k+k_p}[T_{i+i_p+1/2}] 1144 \right) \; \delta_{i+1/2}[T^{k}] \\ 1145 - &\left( 1146 {e_{2u}}_i^{k+k_p+1/2} \ \ {A_{e}}_i^{k+k_p+1/2} 1147 \ \ { _i^{k+k_p+1/2} \mathbb{R}_{i_p}^{-k_p} } \ \ \delta_{i+i_p}[T^{k+k_p+1/2}] 1148 \right) \; \delta_{k+1/2}[T_{i}] \ \Biggr\} 1149 \end{align*} 1150 Expending the summation on $i_p$ and $k_p$, it becomes: 1151 \begin{align*} 1152 \begin{matrix} 1153 &\sum\limits_{i,k} \Bigl\{ 1154 &+{e_{2u}}_{i+1}^{k} &{A_{e}}_{i+1 }^{k} 1155 &\ {_{i+1}^k \mathbb{R}_{- 1/2}^{-1/2}} &\delta_{k-1/2}[T_{i+1}] &\delta_{i+1/2}[T^{k}] &\\ 1156 &&+{e_{2u}}_i^{k\ \ \ \:} &{A_{e}}_{i}^{k\ \ \ \:} 1157 &\ {\ \ \;_i^k \mathbb{R}_{+1/2}^{-1/2}} &\delta_{k-1/2}[T_{i\ \ \ \;}] &\delta_{i+1/2}[T^{k}] &\\ 1158 &&+{e_{2u}}_{i+1}^{k} &{A_{e}}_{i+1 }^{k} 1159 &\ {_{i+1}^k \mathbb{R}_{- 1/2}^{+1/2}} &\delta_{k+1/2}[T_{i+1}] &\delta_{i+1/2}[T^{k}] &\\ 1160 &&+{e_{2u}}_i^{k\ \ \ \:} &{A_{e}}_{i}^{k\ \ \ \:} 1161 &\ {\ \ \;_i^k \mathbb{R}_{+1/2}^{+1/2}} &\delta_{k+1/2}[T_{i\ \ \ \;}] &\delta_{i+1/2}[T^{k}] &\\ 1162 % 1163 &&-{e_{2u}}_i^{k+1} &{A_{e}}_i^{k+1} 1164 &{_i^{k+1} \mathbb{R}_{-1/2}^{- 1/2}} &\delta_{i-1/2}[T^{k+1}] &\delta_{k+1/2}[T_{i}] &\\ 1165 &&-{e_{2u}}_i^{k\ \ \ \:} &{A_{e}}_i^{k\ \ \ \:} 1166 &{\ \ \;_i^k \mathbb{R}_{-1/2}^{+1/2}} &\delta_{i-1/2}[T^{k\ \ \ \:}] &\delta_{k+1/2}[T_{i}] &\\ 1167 &&-{e_{2u}}_i^{k+1 } &{A_{e}}_i^{k+1} 1168 &{_i^{k+1} \mathbb{R}_{+1/2}^{- 1/2}} &\delta_{i+1/2}[T^{k+1}] &\delta_{k+1/2}[T_{i}] &\\ 1169 &&-{e_{2u}}_i^{k\ \ \ \:} &{A_{e}}_i^{k\ \ \ \:} 1170 &{\ \ \;_i^k \mathbb{R}_{+1/2}^{+1/2}} &\delta_{i+1/2}[T^{k\ \ \ \:}] &\delta_{k+1/2}[T_{i}] 1171 &\Bigr\} \\ 1172 \end{matrix} 1173 \end{align*} 1174 The two terms associated with the triad ${_i^k \mathbb{R}_{+1/2}^{+1/2}}$ are the 1175 same but of opposite signs, they cancel out. 1176 Exactly the same thing occurs for the triad ${_i^k \mathbb{R}_{-1/2}^{-1/2}}$. 1177 The two terms associated with the triad ${_i^k \mathbb{R}_{+1/2}^{-1/2}}$ are the 1178 same but both of opposite signs and shifted by 1 in $k$ direction. When summing over $k$ 1179 they cancel out with the neighbouring grid points. 1180 Exactly the same thing occurs for the triad ${_i^k \mathbb{R}_{-1/2}^{+1/2}}$ in the 1181 $i$ direction. Therefore the sum over the domain is zero, $i.e.$ the variance of the 1182 tracer is preserved by the discretisation of the skew fluxes. 1183 1184 %%% Local Variables: 1185 %%% TeX-master: "../../NEMO_book-luatex.tex" 1186 %%% End: 1180 {\psi_1}_{i+1/2}^{k+1/2}={\quarter}\sum_{\substack{i_p,\,k_p}} 1181 {A_e}_{i+1/2-i_p}^{k+1/2-k_p}\:\triadd{i+1/2-i_p}{k+1/2-k_p}{R}{i_p}{k_p}. 1182 \end{equation} 1183 The streamfunction $\psi_1$ is calculated similarly at $vw$ points. 1184 The eddy-induced velocities are then calculated from the 1185 straightforward discretisation of \eqref{eq:triad:eiv_v}: 1186 \begin{equation}\label{eq:triad:eiv_v} 1187 \begin{split} 1188 {u^*}_{i+1/2}^{k} & = - \frac{1}{{e_{3u}}_{i}^{k}}\left({\psi_1}_{i+1/2}^{k+1/2}-{\psi_1}_{i+1/2}^{k+1/2}\right), \\ 1189 {v^*}_{j+1/2}^{k} & = - \frac{1}{{e_{3v}}_{j}^{k}}\left({\psi_2}_{j+1/2}^{k+1/2}-{\psi_2}_{j+1/2}^{k+1/2}\right), \\ 1190 {w^*}_{i,j}^{k+1/2} & = \frac{1}{e_{1t}e_{2t}}\; \left\{ 1191 {e_{2u}}_{i+1/2}^{k+1/2} \,{\psi_1}_{i+1/2}^{k+1/2} - 1192 {e_{2u}}_{i-1/2}^{k+1/2} \,{\psi_1}_{i-1/2}^{k+1/2} \right. + \\ 1193 \phantom{=} & \qquad\qquad\left. {e_{2v}}_{j+1/2}^{k+1/2} \,{\psi_2}_{j+1/2}^{k+1/2} - {e_{2v}}_{j-1/2}^{k+1/2} \,{\psi_2}_{j-1/2}^{k+1/2} \right\}, 1194 \end{split} 1195 \end{equation} -
branches/2011/dev_NEMO_MERGE_2011/DOC/TexFiles/Chapters/Chap_DOM.tex
r3221 r3267 123 123 the following discrete forms in the curvilinear $s$-coordinate system: 124 124 \begin{equation} \label{Eq_DOM_grad} 125 \nabla q\equiv \frac{1}{e_{1u} } \delta _{i+1/2 } [q] \;\, {\rm {\bf i}}126 + \frac{1}{e_{2v} } \delta _{j+1/2 } [q] \;\, {\rm {\bf j}}127 + \frac{1}{e_{3w}} \delta _{k+1/2} [q] \;\, {\rm {\bf k}}125 \nabla q\equiv \frac{1}{e_{1u} } \delta _{i+1/2 } [q] \;\,\mathbf{i} 126 + \frac{1}{e_{2v} } \delta _{j+1/2 } [q] \;\,\mathbf{j} 127 + \frac{1}{e_{3w}} \delta _{k+1/2} [q] \;\,\mathbf{k} 128 128 \end{equation} 129 129 \begin{multline} \label{Eq_DOM_lap} 130 \Delta q\equiv \frac{1}{e_{1t}\,e_{2t}\,e_{3t} } 130 \Delta q\equiv \frac{1}{e_{1t}\,e_{2t}\,e_{3t} } 131 131 \;\left( \delta_i \left[ \frac{e_{2u}\,e_{3u}} {e_{1u}} \;\delta_{i+1/2} [q] \right] 132 132 + \delta_j \left[ \frac{e_{1v}\,e_{3v}} {e_{2v}} \;\delta_{j+1/2} [q] \right] \; \right) \\ … … 139 139 \begin{eqnarray} \label{Eq_DOM_curl} 140 140 \nabla \times {\rm {\bf A}}\equiv & 141 \frac{1}{e_{2v}\,e_{3vw} } \ \left( \delta_{j +1/2} \left[e_{3w}\,a_3 \right] -\delta_{k+1/2} \left[e_{2v} \,a_2 \right] \right) &\ \ rm{\bfi} \\142 +& \frac{1}{e_{2u}\,e_{3uw}} \ \left( \delta_{k+1/2} \left[e_{1u}\,a_1 \right] -\delta_{i +1/2} \left[e_{3w}\,a_3 \right] \right) &\ \ rm{\bf j} \\143 +& \frac{1}{e_{1f} \,e_{2f} } \ \left( \delta_{i +1/2} \left[e_{2v}\,a_2 \right] -\delta_{j +1/2} \left[e_{1u}\,a_1 \right] \right) &\ \ rm{\bfk}141 \frac{1}{e_{2v}\,e_{3vw} } \ \left( \delta_{j +1/2} \left[e_{3w}\,a_3 \right] -\delta_{k+1/2} \left[e_{2v} \,a_2 \right] \right) &\ \mathbf{i} \\ 142 +& \frac{1}{e_{2u}\,e_{3uw}} \ \left( \delta_{k+1/2} \left[e_{1u}\,a_1 \right] -\delta_{i +1/2} \left[e_{3w}\,a_3 \right] \right) &\ \mathbf{j} \\ 143 +& \frac{1}{e_{1f} \,e_{2f} } \ \left( \delta_{i +1/2} \left[e_{2v}\,a_2 \right] -\delta_{j +1/2} \left[e_{1u}\,a_1 \right] \right) &\ \mathbf{k} 144 144 \end{eqnarray} 145 145 \begin{equation} \label{Eq_DOM_div} -
branches/2011/dev_NEMO_MERGE_2011/DOC/TexFiles/Chapters/Chap_LDF.tex
r3221 r3267 21 21 and for tracers only, eddy induced advection on tracers). These three aspects 22 22 of the lateral diffusion are set through namelist parameters and CPP keys 23 (see the \textit{nam\_traldf} and \textit{nam\_dynldf} below). 23 (see the \textit{nam\_traldf} and \textit{nam\_dynldf} below). Note 24 that this chapter describes the default implementation of iso-neutral 25 tracer mixing, and Griffies's implementation, which is used if 26 \np{traldf\_grif}=true, is described in Appdx\ref{sec:triad} 24 27 25 28 %-----------------------------------nam_traldf - nam_dynldf-------------------------------------------- … … 128 131 $\ $\newline % force a new ligne 129 132 130 A space variation in the eddy coefficient appeals several remarks:133 The following points are relevant when the eddy coefficient varies spatially: 131 134 132 135 (1) the momentum diffusion operator acting along model level surfaces is 133 136 written in terms of curl and divergent components of the horizontal current 134 (see \S\ref{PE_ldf}). Although the eddy coefficient c anbe set to different values135 in these two terms, this option is not available.137 (see \S\ref{PE_ldf}). Although the eddy coefficient could be set to different values 138 in these two terms, this option is not currently available. 136 139 137 140 (2) with an horizontally varying viscosity, the quadratic integral constraints … … 275 278 276 279 \item[$s$- or hybrid $s$-$z$- coordinate : ] in the current release of \NEMO, 277 there is no specific treatment for iso-neutral mixing in the $s$-coordinate. 280 iso-neutral mixing is only employed for $s$-coordinates if the 281 Griffies scheme is used (\np{traldf\_grif}=true; see Appdx \ref{sec:triad}). 278 282 In other words, iso-neutral mixing will only be accurately represented with a 279 283 linear equation of state (\np{nn\_eos}=1 or 2). In the case of a "true" equation … … 332 336 \end{description} 333 337 334 This implementation is a rather old one. It is similar to the one proposed 335 by Cox [1987], except for the background horizontal diffusion. Indeed, 336 the Cox implementation of isopycnal diffusion in GFDL-type models requires 337 a minimum background horizontal diffusion for numerical stability reasons. 338 To overcome this problem, several techniques have been proposed in which 339 the numerical schemes of the ocean model are modified \citep{Weaver_Eby_JPO97, 340 Griffies_al_JPO98}. Here, another strategy has been chosen \citep{Lazar_PhD97}: 341 a local filtering of the iso-neutral slopes (made on 9 grid-points) prevents 342 the development of grid point noise generated by the iso-neutral diffusion 343 operator (Fig.~\ref{Fig_LDF_ZDF1}). This allows an iso-neutral diffusion scheme 344 without additional background horizontal mixing. This technique can be viewed 345 as a diffusion operator that acts along large-scale (2~$\Delta$x) 346 \gmcomment{2deltax doesnt seem very large scale} 347 iso-neutral surfaces. The diapycnal diffusion required for numerical stability is 348 thus minimized and its net effect on the flow is quite small when compared to 349 the effect of an horizontal background mixing. 338 This implementation is a rather old one. It is similar to the one 339 proposed by Cox [1987], except for the background horizontal 340 diffusion. Indeed, the Cox implementation of isopycnal diffusion in 341 GFDL-type models requires a minimum background horizontal diffusion 342 for numerical stability reasons. To overcome this problem, several 343 techniques have been proposed in which the numerical schemes of the 344 ocean model are modified \citep{Weaver_Eby_JPO97, 345 Griffies_al_JPO98}. Griffies's scheme is now available in \NEMO if 346 \np{traldf\_grif\_iso} is set true; see Appdx \ref{sec:triad}. Here, 347 another strategy is presented \citep{Lazar_PhD97}: a local 348 filtering of the iso-neutral slopes (made on 9 grid-points) prevents 349 the development of grid point noise generated by the iso-neutral 350 diffusion operator (Fig.~\ref{Fig_LDF_ZDF1}). This allows an 351 iso-neutral diffusion scheme without additional background horizontal 352 mixing. This technique can be viewed as a diffusion operator that acts 353 along large-scale (2~$\Delta$x) \gmcomment{2deltax doesnt seem very 354 large scale} iso-neutral surfaces. The diapycnal diffusion required 355 for numerical stability is thus minimized and its net effect on the 356 flow is quite small when compared to the effect of an horizontal 357 background mixing. 350 358 351 359 Nevertheless, this iso-neutral operator does not ensure that variance cannot increase, -
branches/2011/dev_NEMO_MERGE_2011/DOC/TexFiles/Chapters/Chap_TRA.tex
r3221 r3267 491 491 \label{TRA_ldf_iso} 492 492 493 The general form of the second order lateral tracer subgrid scale physics 493 If the Griffies trad scheme is not employed 494 (\np{ln\_traldf\_grif}=true; see App.\ref{sec:triad}) the general form of the second order lateral tracer subgrid scale physics 494 495 (\ref{Eq_PE_zdf}) takes the following semi-discrete space form in $z$- and 495 496 $s$-coordinates: … … 536 537 (see \S\ref{LDF}) allows the model to run safely without any additional 537 538 background horizontal diffusion \citep{Guilyardi_al_CD01}. An alternative scheme 538 developed by \cite{Griffies_al_JPO98} which preserves both tracer and its variance539 developed by \cite{Griffies_al_JPO98} which ensures tracer variance decreases 539 540 is also available in \NEMO (\np{ln\_traldf\_grif}=true). A complete description of 540 the algorithm is given in App.\ref{ Apdx_Griffies}.541 the algorithm is given in App.\ref{sec:triad}. 541 542 542 543 Note that in the partial step $z$-coordinate (\np{ln\_zps}=true), the horizontal
Note: See TracChangeset
for help on using the changeset viewer.