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/chap_ZDF.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/chap_ZDF.tex

    r11596 r11597  
    1212%gm% Add here a small introduction to ZDF and naming of the different physics (similar to what have been written for TRA and DYN. 
    1313 
     14%% ================================================================================================= 
    1415\section{Vertical mixing} 
    1516\label{sec:ZDF} 
     
    3839%and thus of the formulation used (see \autoref{chap:TD}). 
    3940 
    40 %--------------------------------------------namzdf-------------------------------------------------------- 
    4141 
    4242\begin{listing} 
     
    4545  \label{lst:namzdf} 
    4646\end{listing} 
    47 %-------------------------------------------------------------------------------------------------------------- 
    48  
     47 
     48%% ================================================================================================= 
    4949\subsection[Constant (\forcode{ln_zdfcst})]{Constant (\protect\np{ln_zdfcst}{ln\_zdfcst})} 
    5050\label{subsec:ZDF_cst} 
     
    6666$\sim10^{-9}~m^2.s^{-1}$ for salinity. 
    6767 
     68%% ================================================================================================= 
    6869\subsection[Richardson number dependent (\forcode{ln_zdfric})]{Richardson number dependent (\protect\np{ln_zdfric}{ln\_zdfric})} 
    6970\label{subsec:ZDF_ric} 
    7071 
    71 %--------------------------------------------namric--------------------------------------------------------- 
    7272 
    7373\begin{listing} 
     
    7676  \label{lst:namzdf_ric} 
    7777\end{listing} 
    78 %-------------------------------------------------------------------------------------------------------------- 
    7978 
    8079When \np[=.true.]{ln_zdfric}{ln\_zdfric}, a local Richardson number dependent formulation for the vertical momentum and 
     
    124123the empirical values \np{rn_wtmix}{rn\_wtmix} and \np{rn_wvmix}{rn\_wvmix} \citep{lermusiaux_JMS01}. 
    125124 
     125%% ================================================================================================= 
    126126\subsection[TKE turbulent closure scheme (\forcode{ln_zdftke})]{TKE turbulent closure scheme (\protect\np{ln_zdftke}{ln\_zdftke})} 
    127127\label{subsec:ZDF_tke} 
    128 %--------------------------------------------namzdf_tke-------------------------------------------------- 
    129128 
    130129\begin{listing} 
     
    133132  \label{lst:namzdf_tke} 
    134133\end{listing} 
    135 %-------------------------------------------------------------------------------------------------------------- 
    136134 
    137135The vertical eddy viscosity and diffusivity coefficients are computed from a TKE turbulent closure model based on 
     
    196194\np{rn_avt0}{rn\_avt0} (\nam{zdf}{zdf} namelist, see \autoref{subsec:ZDF_cst}). 
    197195 
     196%% ================================================================================================= 
    198197\subsubsection{Turbulent length scale} 
    199198 
     
    266265$\bar{e}$ reach its minimum value ($1.10^{-6}= C_k\, l_{min} \,\sqrt{\bar{e}_{min}}$ ). 
    267266 
     267%% ================================================================================================= 
    268268\subsubsection{Surface wave breaking parameterization} 
    269 %-----------------------------------------------------------------------% 
    270269 
    271270Following \citet{mellor.blumberg_JPO04}, the TKE turbulence closure model has been modified to 
     
    300299surface $\bar{e}$ value. 
    301300 
     301%% ================================================================================================= 
    302302\subsubsection{Langmuir cells} 
    303 %--------------------------------------% 
    304303 
    305304Langmuir circulations (LC) can be described as ordered large-scale vertical motions in 
     
    354353\] 
    355354 
     355%% ================================================================================================= 
    356356\subsubsection{Mixing just below the mixed layer} 
    357 %--------------------------------------------------------------% 
    358357 
    359358Vertical mixing parameterizations commonly used in ocean general circulation models tend to 
     
    402401% (\eg\ Mellor, 1989; Large et al., 1994; Meier, 2001; Axell, 2002; St. Laurent and Garrett, 2002). 
    403402 
     403%% ================================================================================================= 
    404404\subsection[GLS: Generic Length Scale (\forcode{ln_zdfgls})]{GLS: Generic Length Scale (\protect\np{ln_zdfgls}{ln\_zdfgls})} 
    405405\label{subsec:ZDF_gls} 
    406406 
    407 %--------------------------------------------namzdf_gls--------------------------------------------------------- 
    408407 
    409408\begin{listing} 
     
    412411  \label{lst:namzdf_gls} 
    413412\end{listing} 
    414 %-------------------------------------------------------------------------------------------------------------- 
    415413 
    416414The Generic Length Scale (GLS) scheme is a turbulent closure scheme based on two prognostic equations: 
     
    463461They are made available through the \np{nn_clo}{nn\_clo} namelist parameter. 
    464462 
    465 %--------------------------------------------------TABLE-------------------------------------------------- 
    466463\begin{table}[htbp] 
    467464  \centering 
     
    490487  \label{tab:ZDF_GLS} 
    491488\end{table} 
    492 %-------------------------------------------------------------------------------------------------------------- 
    493489 
    494490In the Mellor-Yamada model, the negativity of $n$ allows to use a wall function to force the convergence of 
     
    522518 in \citet{reffray.guillaume.ea_GMD15} for the \NEMO\ model. 
    523519 
     520%% ================================================================================================= 
    524521\subsection[OSM: OSMosis boundary layer scheme (\forcode{ln_zdfosm})]{OSM: OSMosis boundary layer scheme (\protect\np{ln_zdfosm}{ln\_zdfosm})} 
    525522\label{subsec:ZDF_osm} 
    526 %--------------------------------------------namzdf_osm--------------------------------------------------------- 
    527523 
    528524\begin{listing} 
     
    531527  \label{lst:namzdf_osm} 
    532528\end{listing} 
    533 %-------------------------------------------------------------------------------------------------------------- 
    534529 
    535530The OSMOSIS turbulent closure scheme is based on......   TBC 
    536531 
     532%% ================================================================================================= 
    537533\subsection[ Discrete energy conservation for TKE and GLS schemes]{Discrete energy conservation for TKE and GLS schemes} 
    538534\label{subsec:ZDF_tke_ene} 
     
    635631%For the latter, it is in fact the ratio $\sqrt{\bar{e}}/l_\epsilon$ which is stored. 
    636632 
     633%% ================================================================================================= 
    637634\section{Convection} 
    638635\label{sec:ZDF_conv} 
     
    645642or/and the use of a turbulent closure scheme. 
    646643 
     644%% ================================================================================================= 
    647645\subsection[Non-penetrative convective adjustment (\forcode{ln_tranpc})]{Non-penetrative convective adjustment (\protect\np{ln_tranpc}{ln\_tranpc})} 
    648646\label{subsec:ZDF_npc} 
     
    705703having to recompute the expansion coefficients at each mixing iteration. 
    706704 
     705%% ================================================================================================= 
    707706\subsection[Enhanced vertical diffusion (\forcode{ln_zdfevd})]{Enhanced vertical diffusion (\protect\np{ln_zdfevd}{ln\_zdfevd})} 
    708707\label{subsec:ZDF_evd} 
     
    728727a leapfrog environment \citep{leclair_phd10} (see \autoref{sec:TD_mLF}). 
    729728 
     729%% ================================================================================================= 
    730730\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}})} 
    731731\label{subsec:ZDF_tcs} 
     
    752752% gm%  + one word on non local flux with KPP scheme trakpp.F90 module... 
    753753 
     754%% ================================================================================================= 
    754755\section[Double diffusion mixing (\forcode{ln_zdfddm})]{Double diffusion mixing (\protect\np{ln_zdfddm}{ln\_zdfddm})} 
    755756\label{subsec:ZDF_ddm} 
    756757 
    757 %-------------------------------------------namzdf_ddm------------------------------------------------- 
    758758% 
    759759%\nlst{namzdf_ddm} 
    760 %-------------------------------------------------------------------------------------------------------------- 
    761760 
    762761This parameterisation has been introduced in \mdl{zdfddm} module and is controlled by the namelist parameter 
     
    840839This avoids duplication in the computation of $\alpha$ and $\beta$ (which is usually quite expensive). 
    841840 
     841%% ================================================================================================= 
    842842\section[Bottom and top friction (\textit{zdfdrg.F90})]{Bottom and top friction (\protect\mdl{zdfdrg})} 
    843843\label{sec:ZDF_drg} 
    844844 
    845 %--------------------------------------------namdrg-------------------------------------------------------- 
    846845% 
    847846\begin{listing} 
     
    861860\end{listing} 
    862861 
    863 %-------------------------------------------------------------------------------------------------------------- 
    864862 
    865863Options to define the top and bottom friction are defined through the \nam{drg}{drg} namelist variables. 
     
    916914Note 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. 
    917915 
     916%% ================================================================================================= 
    918917\subsection[Linear top/bottom friction (\forcode{ln_lin})]{Linear top/bottom friction (\protect\np{ln_lin}{ln\_lin})} 
    919918\label{subsec:ZDF_drg_linear} 
     
    952951$mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 
    953952 
     953%% ================================================================================================= 
    954954\subsection[Non-linear top/bottom friction (\forcode{ln_non_lin})]{Non-linear top/bottom friction (\protect\np{ln_non_lin}{ln\_non\_lin})} 
    955955\label{subsec:ZDF_drg_nonlinear} 
     
    984984$mask\_value$ * \np{rn_boost}{rn\_boost} * \np{rn_Cd0}{rn\_Cd0}. 
    985985 
     986%% ================================================================================================= 
    986987\subsection[Log-layer top/bottom friction (\forcode{ln_loglayer})]{Log-layer top/bottom friction (\protect\np{ln_loglayer}{ln\_loglayer})} 
    987988\label{subsec:ZDF_drg_loglayer} 
     
    10071008%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}. 
    10081009 
     1010%% ================================================================================================= 
    10091011\subsection[Explicit top/bottom friction (\forcode{ln_drgimp=.false.})]{Explicit top/bottom friction (\protect\np[=.false.]{ln_drgimp}{ln\_drgimp})} 
    10101012\label{subsec:ZDF_drg_stability} 
     
    10651067The number of potential breaches of the explicit stability criterion are still reported for information purposes. 
    10661068 
     1069%% ================================================================================================= 
    10671070\subsection[Implicit top/bottom friction (\forcode{ln_drgimp=.true.})]{Implicit top/bottom friction (\protect\np[=.true.]{ln_drgimp}{ln\_drgimp})} 
    10681071\label{subsec:ZDF_drg_imp} 
     
    10921095Superscript $n+1$ means the velocity used in the friction formula is to be calculated, so it is implicit. 
    10931096 
     1097%% ================================================================================================= 
    10941098\subsection[Bottom friction with split-explicit free surface]{Bottom friction with split-explicit free surface} 
    10951099\label{subsec:ZDF_drg_ts} 
     
    11051109Note that other strategies are possible, like considering vertical diffusion step in advance, \ie\ prior barotropic integration. 
    11061110 
     1111%% ================================================================================================= 
    11071112\section[Internal wave-driven mixing (\forcode{ln_zdfiwm})]{Internal wave-driven mixing (\protect\np{ln_zdfiwm}{ln\_zdfiwm})} 
    11081113\label{subsec:ZDF_tmx_new} 
    11091114 
    1110 %--------------------------------------------namzdf_iwm------------------------------------------ 
    11111115% 
    11121116\begin{listing} 
     
    11151119  \label{lst:namzdf_iwm} 
    11161120\end{listing} 
    1117 %-------------------------------------------------------------------------------------------------------------- 
    11181121 
    11191122The parameterization of mixing induced by breaking internal waves is a generalization of 
     
    11661169% Jc: input files names ? 
    11671170 
     1171%% ================================================================================================= 
    11681172\section[Surface wave-induced mixing (\forcode{ln_zdfswm})]{Surface wave-induced mixing (\protect\np{ln_zdfswm}{ln\_zdfswm})} 
    11691173\label{subsec:ZDF_swm} 
     
    11961200(for more information on wave parameters and settings see \autoref{sec:SBC_wave}) 
    11971201 
     1202%% ================================================================================================= 
    11981203\section[Adaptive-implicit vertical advection (\forcode{ln_zad_Aimp})]{Adaptive-implicit vertical advection(\protect\np{ln_zad_Aimp}{ln\_zad\_Aimp})} 
    11991204\label{subsec:ZDF_aimp} 
     
    13191324\end{figure} 
    13201325 
     1326%% ================================================================================================= 
    13211327\subsection{Adaptive-implicit vertical advection in the OVERFLOW test-case} 
    13221328 
Note: See TracChangeset for help on using the changeset viewer.