Changeset 11596 for NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex
- Timestamp:
- 2019-09-25T19:06:37+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex
r11584 r11596 5 5 6 6 \begin{document} 7 % ================================================================8 % Chapter Vertical Ocean Physics (ZDF)9 % ================================================================10 7 \chapter{Vertical Ocean Physics (ZDF)} 11 8 \label{chap:ZDF} … … 15 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. 16 13 17 \newpage18 19 % ================================================================20 % Vertical Mixing21 % ================================================================22 14 \section{Vertical mixing} 23 15 \label{sec:ZDF} … … 55 47 %-------------------------------------------------------------------------------------------------------------- 56 48 57 % -------------------------------------------------------------------------------------------------------------58 % Constant59 % -------------------------------------------------------------------------------------------------------------60 49 \subsection[Constant (\forcode{ln_zdfcst})]{Constant (\protect\np{ln_zdfcst}{ln\_zdfcst})} 61 50 \label{subsec:ZDF_cst} … … 77 66 $\sim10^{-9}~m^2.s^{-1}$ for salinity. 78 67 79 % -------------------------------------------------------------------------------------------------------------80 % Richardson Number Dependent81 % -------------------------------------------------------------------------------------------------------------82 68 \subsection[Richardson number dependent (\forcode{ln_zdfric})]{Richardson number dependent (\protect\np{ln_zdfric}{ln\_zdfric})} 83 69 \label{subsec:ZDF_ric} … … 138 124 the empirical values \np{rn_wtmix}{rn\_wtmix} and \np{rn_wvmix}{rn\_wvmix} \citep{lermusiaux_JMS01}. 139 125 140 % -------------------------------------------------------------------------------------------------------------141 % TKE Turbulent Closure Scheme142 % -------------------------------------------------------------------------------------------------------------143 126 \subsection[TKE turbulent closure scheme (\forcode{ln_zdftke})]{TKE turbulent closure scheme (\protect\np{ln_zdftke}{ln\_zdftke})} 144 127 \label{subsec:ZDF_tke} … … 248 231 evaluate the dissipation and mixing length scales as 249 232 (and note that here we use numerical indexing): 250 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>251 233 \begin{figure}[!t] 252 234 \centering … … 255 237 \label{fig:ZDF_mixing_length} 256 238 \end{figure} 257 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>258 239 \[ 259 240 % \label{eq:ZDF_tke_mxl2} … … 421 402 % (\eg\ Mellor, 1989; Large et al., 1994; Meier, 2001; Axell, 2002; St. Laurent and Garrett, 2002). 422 403 423 % -------------------------------------------------------------------------------------------------------------424 % GLS Generic Length Scale Scheme425 % -------------------------------------------------------------------------------------------------------------426 404 \subsection[GLS: Generic Length Scale (\forcode{ln_zdfgls})]{GLS: Generic Length Scale (\protect\np{ln_zdfgls}{ln\_zdfgls})} 427 405 \label{subsec:ZDF_gls} … … 544 522 in \citet{reffray.guillaume.ea_GMD15} for the \NEMO\ model. 545 523 546 547 % -------------------------------------------------------------------------------------------------------------548 % OSM OSMOSIS BL Scheme549 % -------------------------------------------------------------------------------------------------------------550 524 \subsection[OSM: OSMosis boundary layer scheme (\forcode{ln_zdfosm})]{OSM: OSMosis boundary layer scheme (\protect\np{ln_zdfosm}{ln\_zdfosm})} 551 525 \label{subsec:ZDF_osm} … … 561 535 The OSMOSIS turbulent closure scheme is based on...... TBC 562 536 563 % -------------------------------------------------------------------------------------------------------------564 % TKE and GLS discretization considerations565 % -------------------------------------------------------------------------------------------------------------566 537 \subsection[ Discrete energy conservation for TKE and GLS schemes]{Discrete energy conservation for TKE and GLS schemes} 567 538 \label{subsec:ZDF_tke_ene} 568 539 569 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>570 540 \begin{figure}[!t] 571 541 \centering … … 576 546 \label{fig:ZDF_TKE_time_scheme} 577 547 \end{figure} 578 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>579 548 580 549 The production of turbulence by vertical shear (the first term of the right hand side of … … 666 635 %For the latter, it is in fact the ratio $\sqrt{\bar{e}}/l_\epsilon$ which is stored. 667 636 668 % ================================================================669 % Convection670 % ================================================================671 637 \section{Convection} 672 638 \label{sec:ZDF_conv} … … 679 645 or/and the use of a turbulent closure scheme. 680 646 681 % -------------------------------------------------------------------------------------------------------------682 % Non-Penetrative Convective Adjustment683 % -------------------------------------------------------------------------------------------------------------684 647 \subsection[Non-penetrative convective adjustment (\forcode{ln_tranpc})]{Non-penetrative convective adjustment (\protect\np{ln_tranpc}{ln\_tranpc})} 685 648 \label{subsec:ZDF_npc} 686 649 687 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>688 650 \begin{figure}[!htb] 689 651 \centering … … 704 666 \label{fig:ZDF_npc} 705 667 \end{figure} 706 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>707 668 708 669 Options are defined through the \nam{zdf}{zdf} namelist variables. … … 744 705 having to recompute the expansion coefficients at each mixing iteration. 745 706 746 % -------------------------------------------------------------------------------------------------------------747 % Enhanced Vertical Diffusion748 % -------------------------------------------------------------------------------------------------------------749 707 \subsection[Enhanced vertical diffusion (\forcode{ln_zdfevd})]{Enhanced vertical diffusion (\protect\np{ln_zdfevd}{ln\_zdfevd})} 750 708 \label{subsec:ZDF_evd} … … 770 728 a leapfrog environment \citep{leclair_phd10} (see \autoref{sec:TD_mLF}). 771 729 772 % -------------------------------------------------------------------------------------------------------------773 % Turbulent Closure Scheme774 % -------------------------------------------------------------------------------------------------------------775 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}})} 776 731 \label{subsec:ZDF_tcs} 777 778 732 779 733 The turbulent closure schemes presented in \autoref{subsec:ZDF_tke}, \autoref{subsec:ZDF_gls} and … … 798 752 % gm% + one word on non local flux with KPP scheme trakpp.F90 module... 799 753 800 % ================================================================801 % Double Diffusion Mixing802 % ================================================================803 754 \section[Double diffusion mixing (\forcode{ln_zdfddm})]{Double diffusion mixing (\protect\np{ln_zdfddm}{ln\_zdfddm})} 804 755 \label{subsec:ZDF_ddm} 805 806 756 807 757 %-------------------------------------------namzdf_ddm------------------------------------------------- … … 818 768 it leads to relatively minor changes in circulation but exerts significant regional influences on 819 769 temperature and salinity. 820 821 770 822 771 Diapycnal mixing of S and T are described by diapycnal diffusion coefficients … … 844 793 \end{align} 845 794 846 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>847 795 \begin{figure}[!t] 848 796 \centering … … 861 809 \label{fig:ZDF_ddm} 862 810 \end{figure} 863 %>>>>>>>>>>>>>>>>>>>>>>>>>>>>864 811 865 812 The factor 0.7 in \autoref{eq:ZDF_ddm_f_T} reflects the measured ratio $\alpha F_T /\beta F_S \approx 0.7$ of … … 893 840 This avoids duplication in the computation of $\alpha$ and $\beta$ (which is usually quite expensive). 894 841 895 % ================================================================896 % Bottom Friction897 % ================================================================898 842 \section[Bottom and top friction (\textit{zdfdrg.F90})]{Bottom and top friction (\protect\mdl{zdfdrg})} 899 843 \label{sec:ZDF_drg} … … 924 868 As the friction processes at the top and the bottom are treated in and identical way, 925 869 the description below considers mostly the bottom friction case, if not stated otherwise. 926 927 870 928 871 Both the surface momentum flux (wind stress) and the bottom momentum flux (bottom friction) enter the equations as … … 973 916 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. 974 917 975 % -------------------------------------------------------------------------------------------------------------976 % Linear Bottom Friction977 % -------------------------------------------------------------------------------------------------------------978 918 \subsection[Linear top/bottom friction (\forcode{ln_lin})]{Linear top/bottom friction (\protect\np{ln_lin}{ln\_lin})} 979 919 \label{subsec:ZDF_drg_linear} … … 1012 952 $mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 1013 953 1014 % -------------------------------------------------------------------------------------------------------------1015 % Non-Linear Bottom Friction1016 % -------------------------------------------------------------------------------------------------------------1017 954 \subsection[Non-linear top/bottom friction (\forcode{ln_non_lin})]{Non-linear top/bottom friction (\protect\np{ln_non_lin}{ln\_non\_lin})} 1018 955 \label{subsec:ZDF_drg_nonlinear} … … 1047 984 $mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 1048 985 1049 % -------------------------------------------------------------------------------------------------------------1050 % Bottom Friction Log-layer1051 % -------------------------------------------------------------------------------------------------------------1052 986 \subsection[Log-layer top/bottom friction (\forcode{ln_loglayer})]{Log-layer top/bottom friction (\protect\np{ln_loglayer}{ln\_loglayer})} 1053 987 \label{subsec:ZDF_drg_loglayer} … … 1073 1007 %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}. 1074 1008 1075 % -------------------------------------------------------------------------------------------------------------1076 % Explicit bottom Friction1077 % -------------------------------------------------------------------------------------------------------------1078 1009 \subsection[Explicit top/bottom friction (\forcode{ln_drgimp=.false.})]{Explicit top/bottom friction (\protect\np[=.false.]{ln_drgimp}{ln\_drgimp})} 1079 1010 \label{subsec:ZDF_drg_stability} … … 1134 1065 The number of potential breaches of the explicit stability criterion are still reported for information purposes. 1135 1066 1136 % -------------------------------------------------------------------------------------------------------------1137 % Implicit Bottom Friction1138 % -------------------------------------------------------------------------------------------------------------1139 1067 \subsection[Implicit top/bottom friction (\forcode{ln_drgimp=.true.})]{Implicit top/bottom friction (\protect\np[=.true.]{ln_drgimp}{ln\_drgimp})} 1140 1068 \label{subsec:ZDF_drg_imp} … … 1164 1092 Superscript $n+1$ means the velocity used in the friction formula is to be calculated, so it is implicit. 1165 1093 1166 % -------------------------------------------------------------------------------------------------------------1167 % Bottom Friction with split-explicit free surface1168 % -------------------------------------------------------------------------------------------------------------1169 1094 \subsection[Bottom friction with split-explicit free surface]{Bottom friction with split-explicit free surface} 1170 1095 \label{subsec:ZDF_drg_ts} … … 1180 1105 Note that other strategies are possible, like considering vertical diffusion step in advance, \ie\ prior barotropic integration. 1181 1106 1182 1183 % ================================================================1184 % Internal wave-driven mixing1185 % ================================================================1186 1107 \section[Internal wave-driven mixing (\forcode{ln_zdfiwm})]{Internal wave-driven mixing (\protect\np{ln_zdfiwm}{ln\_zdfiwm})} 1187 1108 \label{subsec:ZDF_tmx_new} … … 1245 1166 % Jc: input files names ? 1246 1167 1247 % ================================================================1248 % surface wave-induced mixing1249 % ================================================================1250 1168 \section[Surface wave-induced mixing (\forcode{ln_zdfswm})]{Surface wave-induced mixing (\protect\np{ln_zdfswm}{ln\_zdfswm})} 1251 1169 \label{subsec:ZDF_swm} … … 1278 1196 (for more information on wave parameters and settings see \autoref{sec:SBC_wave}) 1279 1197 1280 % ================================================================1281 % Adaptive-implicit vertical advection1282 % ================================================================1283 1198 \section[Adaptive-implicit vertical advection (\forcode{ln_zad_Aimp})]{Adaptive-implicit vertical advection(\protect\np{ln_zad_Aimp}{ln\_zad\_Aimp})} 1284 1199 \label{subsec:ZDF_aimp}
Note: See TracChangeset
for help on using the changeset viewer.