Changeset 11597 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.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/chap_ZDF.tex
r11596 r11597 12 12 %gm% Add here a small introduction to ZDF and naming of the different physics (similar to what have been written for TRA and DYN. 13 13 14 %% ================================================================================================= 14 15 \section{Vertical mixing} 15 16 \label{sec:ZDF} … … 38 39 %and thus of the formulation used (see \autoref{chap:TD}). 39 40 40 %--------------------------------------------namzdf--------------------------------------------------------41 41 42 42 \begin{listing} … … 45 45 \label{lst:namzdf} 46 46 \end{listing} 47 %-------------------------------------------------------------------------------------------------------------- 48 47 48 %% ================================================================================================= 49 49 \subsection[Constant (\forcode{ln_zdfcst})]{Constant (\protect\np{ln_zdfcst}{ln\_zdfcst})} 50 50 \label{subsec:ZDF_cst} … … 66 66 $\sim10^{-9}~m^2.s^{-1}$ for salinity. 67 67 68 %% ================================================================================================= 68 69 \subsection[Richardson number dependent (\forcode{ln_zdfric})]{Richardson number dependent (\protect\np{ln_zdfric}{ln\_zdfric})} 69 70 \label{subsec:ZDF_ric} 70 71 71 %--------------------------------------------namric---------------------------------------------------------72 72 73 73 \begin{listing} … … 76 76 \label{lst:namzdf_ric} 77 77 \end{listing} 78 %--------------------------------------------------------------------------------------------------------------79 78 80 79 When \np[=.true.]{ln_zdfric}{ln\_zdfric}, a local Richardson number dependent formulation for the vertical momentum and … … 124 123 the empirical values \np{rn_wtmix}{rn\_wtmix} and \np{rn_wvmix}{rn\_wvmix} \citep{lermusiaux_JMS01}. 125 124 125 %% ================================================================================================= 126 126 \subsection[TKE turbulent closure scheme (\forcode{ln_zdftke})]{TKE turbulent closure scheme (\protect\np{ln_zdftke}{ln\_zdftke})} 127 127 \label{subsec:ZDF_tke} 128 %--------------------------------------------namzdf_tke--------------------------------------------------129 128 130 129 \begin{listing} … … 133 132 \label{lst:namzdf_tke} 134 133 \end{listing} 135 %--------------------------------------------------------------------------------------------------------------136 134 137 135 The vertical eddy viscosity and diffusivity coefficients are computed from a TKE turbulent closure model based on … … 196 194 \np{rn_avt0}{rn\_avt0} (\nam{zdf}{zdf} namelist, see \autoref{subsec:ZDF_cst}). 197 195 196 %% ================================================================================================= 198 197 \subsubsection{Turbulent length scale} 199 198 … … 266 265 $\bar{e}$ reach its minimum value ($1.10^{-6}= C_k\, l_{min} \,\sqrt{\bar{e}_{min}}$ ). 267 266 267 %% ================================================================================================= 268 268 \subsubsection{Surface wave breaking parameterization} 269 %-----------------------------------------------------------------------%270 269 271 270 Following \citet{mellor.blumberg_JPO04}, the TKE turbulence closure model has been modified to … … 300 299 surface $\bar{e}$ value. 301 300 301 %% ================================================================================================= 302 302 \subsubsection{Langmuir cells} 303 %--------------------------------------%304 303 305 304 Langmuir circulations (LC) can be described as ordered large-scale vertical motions in … … 354 353 \] 355 354 355 %% ================================================================================================= 356 356 \subsubsection{Mixing just below the mixed layer} 357 %--------------------------------------------------------------%358 357 359 358 Vertical mixing parameterizations commonly used in ocean general circulation models tend to … … 402 401 % (\eg\ Mellor, 1989; Large et al., 1994; Meier, 2001; Axell, 2002; St. Laurent and Garrett, 2002). 403 402 403 %% ================================================================================================= 404 404 \subsection[GLS: Generic Length Scale (\forcode{ln_zdfgls})]{GLS: Generic Length Scale (\protect\np{ln_zdfgls}{ln\_zdfgls})} 405 405 \label{subsec:ZDF_gls} 406 406 407 %--------------------------------------------namzdf_gls---------------------------------------------------------408 407 409 408 \begin{listing} … … 412 411 \label{lst:namzdf_gls} 413 412 \end{listing} 414 %--------------------------------------------------------------------------------------------------------------415 413 416 414 The Generic Length Scale (GLS) scheme is a turbulent closure scheme based on two prognostic equations: … … 463 461 They are made available through the \np{nn_clo}{nn\_clo} namelist parameter. 464 462 465 %--------------------------------------------------TABLE--------------------------------------------------466 463 \begin{table}[htbp] 467 464 \centering … … 490 487 \label{tab:ZDF_GLS} 491 488 \end{table} 492 %--------------------------------------------------------------------------------------------------------------493 489 494 490 In the Mellor-Yamada model, the negativity of $n$ allows to use a wall function to force the convergence of … … 522 518 in \citet{reffray.guillaume.ea_GMD15} for the \NEMO\ model. 523 519 520 %% ================================================================================================= 524 521 \subsection[OSM: OSMosis boundary layer scheme (\forcode{ln_zdfosm})]{OSM: OSMosis boundary layer scheme (\protect\np{ln_zdfosm}{ln\_zdfosm})} 525 522 \label{subsec:ZDF_osm} 526 %--------------------------------------------namzdf_osm---------------------------------------------------------527 523 528 524 \begin{listing} … … 531 527 \label{lst:namzdf_osm} 532 528 \end{listing} 533 %--------------------------------------------------------------------------------------------------------------534 529 535 530 The OSMOSIS turbulent closure scheme is based on...... TBC 536 531 532 %% ================================================================================================= 537 533 \subsection[ Discrete energy conservation for TKE and GLS schemes]{Discrete energy conservation for TKE and GLS schemes} 538 534 \label{subsec:ZDF_tke_ene} … … 635 631 %For the latter, it is in fact the ratio $\sqrt{\bar{e}}/l_\epsilon$ which is stored. 636 632 633 %% ================================================================================================= 637 634 \section{Convection} 638 635 \label{sec:ZDF_conv} … … 645 642 or/and the use of a turbulent closure scheme. 646 643 644 %% ================================================================================================= 647 645 \subsection[Non-penetrative convective adjustment (\forcode{ln_tranpc})]{Non-penetrative convective adjustment (\protect\np{ln_tranpc}{ln\_tranpc})} 648 646 \label{subsec:ZDF_npc} … … 705 703 having to recompute the expansion coefficients at each mixing iteration. 706 704 705 %% ================================================================================================= 707 706 \subsection[Enhanced vertical diffusion (\forcode{ln_zdfevd})]{Enhanced vertical diffusion (\protect\np{ln_zdfevd}{ln\_zdfevd})} 708 707 \label{subsec:ZDF_evd} … … 728 727 a leapfrog environment \citep{leclair_phd10} (see \autoref{sec:TD_mLF}). 729 728 729 %% ================================================================================================= 730 730 \subsection[Handling convection with turbulent closure schemes (\forcode{ln_zdf_}\{\forcode{tke,gls,osm}\})]{Handling convection with turbulent closure schemes (\forcode{ln_zdf{tke,gls,osm}})} 731 731 \label{subsec:ZDF_tcs} … … 752 752 % gm% + one word on non local flux with KPP scheme trakpp.F90 module... 753 753 754 %% ================================================================================================= 754 755 \section[Double diffusion mixing (\forcode{ln_zdfddm})]{Double diffusion mixing (\protect\np{ln_zdfddm}{ln\_zdfddm})} 755 756 \label{subsec:ZDF_ddm} 756 757 757 %-------------------------------------------namzdf_ddm-------------------------------------------------758 758 % 759 759 %\nlst{namzdf_ddm} 760 %--------------------------------------------------------------------------------------------------------------761 760 762 761 This parameterisation has been introduced in \mdl{zdfddm} module and is controlled by the namelist parameter … … 840 839 This avoids duplication in the computation of $\alpha$ and $\beta$ (which is usually quite expensive). 841 840 841 %% ================================================================================================= 842 842 \section[Bottom and top friction (\textit{zdfdrg.F90})]{Bottom and top friction (\protect\mdl{zdfdrg})} 843 843 \label{sec:ZDF_drg} 844 844 845 %--------------------------------------------namdrg--------------------------------------------------------846 845 % 847 846 \begin{listing} … … 861 860 \end{listing} 862 861 863 %--------------------------------------------------------------------------------------------------------------864 862 865 863 Options to define the top and bottom friction are defined through the \nam{drg}{drg} namelist variables. … … 916 914 Note than from \NEMO\ 4.0, drag coefficients are only computed at cell centers (\ie\ at T-points) and refer to as $c_b^T$ in the following. These are then linearly interpolated in space to get $c_b^\textbf{U}$ at velocity points. 917 915 916 %% ================================================================================================= 918 917 \subsection[Linear top/bottom friction (\forcode{ln_lin})]{Linear top/bottom friction (\protect\np{ln_lin}{ln\_lin})} 919 918 \label{subsec:ZDF_drg_linear} … … 952 951 $mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 953 952 953 %% ================================================================================================= 954 954 \subsection[Non-linear top/bottom friction (\forcode{ln_non_lin})]{Non-linear top/bottom friction (\protect\np{ln_non_lin}{ln\_non\_lin})} 955 955 \label{subsec:ZDF_drg_nonlinear} … … 984 984 $mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 985 985 986 %% ================================================================================================= 986 987 \subsection[Log-layer top/bottom friction (\forcode{ln_loglayer})]{Log-layer top/bottom friction (\protect\np{ln_loglayer}{ln\_loglayer})} 987 988 \label{subsec:ZDF_drg_loglayer} … … 1007 1008 %In this case, the relevant namelist parameters are \np{rn_tfrz0}{rn\_tfrz0}, \np{rn_tfri2}{rn\_tfri2} and \np{rn_tfri2_max}{rn\_tfri2\_max}. 1008 1009 1010 %% ================================================================================================= 1009 1011 \subsection[Explicit top/bottom friction (\forcode{ln_drgimp=.false.})]{Explicit top/bottom friction (\protect\np[=.false.]{ln_drgimp}{ln\_drgimp})} 1010 1012 \label{subsec:ZDF_drg_stability} … … 1065 1067 The number of potential breaches of the explicit stability criterion are still reported for information purposes. 1066 1068 1069 %% ================================================================================================= 1067 1070 \subsection[Implicit top/bottom friction (\forcode{ln_drgimp=.true.})]{Implicit top/bottom friction (\protect\np[=.true.]{ln_drgimp}{ln\_drgimp})} 1068 1071 \label{subsec:ZDF_drg_imp} … … 1092 1095 Superscript $n+1$ means the velocity used in the friction formula is to be calculated, so it is implicit. 1093 1096 1097 %% ================================================================================================= 1094 1098 \subsection[Bottom friction with split-explicit free surface]{Bottom friction with split-explicit free surface} 1095 1099 \label{subsec:ZDF_drg_ts} … … 1105 1109 Note that other strategies are possible, like considering vertical diffusion step in advance, \ie\ prior barotropic integration. 1106 1110 1111 %% ================================================================================================= 1107 1112 \section[Internal wave-driven mixing (\forcode{ln_zdfiwm})]{Internal wave-driven mixing (\protect\np{ln_zdfiwm}{ln\_zdfiwm})} 1108 1113 \label{subsec:ZDF_tmx_new} 1109 1114 1110 %--------------------------------------------namzdf_iwm------------------------------------------1111 1115 % 1112 1116 \begin{listing} … … 1115 1119 \label{lst:namzdf_iwm} 1116 1120 \end{listing} 1117 %--------------------------------------------------------------------------------------------------------------1118 1121 1119 1122 The parameterization of mixing induced by breaking internal waves is a generalization of … … 1166 1169 % Jc: input files names ? 1167 1170 1171 %% ================================================================================================= 1168 1172 \section[Surface wave-induced mixing (\forcode{ln_zdfswm})]{Surface wave-induced mixing (\protect\np{ln_zdfswm}{ln\_zdfswm})} 1169 1173 \label{subsec:ZDF_swm} … … 1196 1200 (for more information on wave parameters and settings see \autoref{sec:SBC_wave}) 1197 1201 1202 %% ================================================================================================= 1198 1203 \section[Adaptive-implicit vertical advection (\forcode{ln_zad_Aimp})]{Adaptive-implicit vertical advection(\protect\np{ln_zad_Aimp}{ln\_zad\_Aimp})} 1199 1204 \label{subsec:ZDF_aimp} … … 1319 1324 \end{figure} 1320 1325 1326 %% ================================================================================================= 1321 1327 \subsection{Adaptive-implicit vertical advection in the OVERFLOW test-case} 1322 1328
Note: See TracChangeset
for help on using the changeset viewer.