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_LDF.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_LDF.tex

    r11596 r11597  
    2222is described in \autoref{apdx:TRIADS} 
    2323 
    24 %-----------------------------------namtra_ldf - namdyn_ldf-------------------------------------------- 
    25  
    26 %-------------------------------------------------------------------------------------------------------------- 
    27  
     24 
     25 
     26%% ================================================================================================= 
    2827\section[Lateral mixing operators]{Lateral mixing operators} 
    2928\label{sec:LDF_op} 
    3029We remind here the different lateral mixing operators that can be used. Further details can be found in \autoref{subsec:TRA_ldf_op} and  \autoref{sec:DYN_ldf}. 
    3130 
     31%% ================================================================================================= 
    3232\subsection[No lateral mixing (\forcode{ln_traldf_OFF} \& \forcode{ln_dynldf_OFF})]{No lateral mixing (\protect\np{ln_traldf_OFF}{ln\_traldf\_OFF} \& \protect\np{ln_dynldf_OFF}{ln\_dynldf\_OFF})} 
    3333 
     
    3737see \autoref{subsec:DYN_adv_ubs}) and can be useful for testing purposes. 
    3838 
     39%% ================================================================================================= 
    3940\subsection[Laplacian mixing (\forcode{ln_traldf_lap} \& \forcode{ln_dynldf_lap})]{Laplacian mixing (\protect\np{ln_traldf_lap}{ln\_traldf\_lap} \& \protect\np{ln_dynldf_lap}{ln\_dynldf\_lap})} 
    4041Setting \protect\np[=.true.]{ln_traldf_lap}{ln\_traldf\_lap} and/or \protect\np[=.true.]{ln_dynldf_lap}{ln\_dynldf\_lap} enables 
     
    4243Laplacian and Bilaplacian operators for the same variable. 
    4344 
     45%% ================================================================================================= 
    4446\subsection[Bilaplacian mixing (\forcode{ln_traldf_blp} \& \forcode{ln_dynldf_blp})]{Bilaplacian mixing (\protect\np{ln_traldf_blp}{ln\_traldf\_blp} \& \protect\np{ln_dynldf_blp}{ln\_dynldf\_blp})} 
    4547Setting \protect\np[=.true.]{ln_traldf_blp}{ln\_traldf\_blp} and/or \protect\np[=.true.]{ln_dynldf_blp}{ln\_dynldf\_blp} enables 
     
    4749We stress again that from \NEMO\ 4, the simultaneous use Laplacian and Bilaplacian operators is not allowed. 
    4850 
     51%% ================================================================================================= 
    4952\section[Direction of lateral mixing (\textit{ldfslp.F90})]{Direction of lateral mixing (\protect\mdl{ldfslp})} 
    5053\label{sec:LDF_slp} 
     
    6972%gm% add here afigure of the slope in i-direction 
    7073 
     74%% ================================================================================================= 
    7175\subsection{Slopes for tracer geopotential mixing in the $s$-coordinate} 
    7276 
     
    99103and either \np[=.true.]{ln_traldf_hor}{ln\_traldf\_hor} or \np[=.true.]{ln_dynldf_hor}{ln\_dynldf\_hor}. 
    100104 
     105%% ================================================================================================= 
    101106\subsection{Slopes for tracer iso-neutral mixing} 
    102107\label{subsec:LDF_slp_iso} 
     
    273278\colorbox{yellow}{add here a discussion about the flattening of the slopes, vs tapering the coefficient.} 
    274279 
     280%% ================================================================================================= 
    275281\subsection{Slopes for momentum iso-neutral mixing} 
    276282 
     
    299305(see \autoref{sec:LBC_coast}). 
    300306 
     307%% ================================================================================================= 
    301308\section[Lateral mixing coefficient (\forcode{nn_aht_ijk_t} \& \forcode{nn_ahm_ijk_t})]{Lateral mixing coefficient (\protect\np{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    302309\label{sec:LDF_coef} 
     
    305312The way the mixing coefficients are set in the reference version can be described as follows: 
    306313 
     314%% ================================================================================================= 
    307315\subsection[Mixing coefficients read from file (\forcode{=-20, -30})]{ Mixing coefficients read from file (\protect\np[=-20, -30]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=-20, -30]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    308316 
     
    313321The provided fields can either be 2d (\np[=-20]{nn_aht_ijk_t}{nn\_aht\_ijk\_t}, \np[=-20]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t}) or 3d (\np[=-30]{nn_aht_ijk_t}{nn\_aht\_ijk\_t},  \np[=-30]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t}). They must be given at U, V points for tracers and T, F points for momentum (see \autoref{tab:LDF_files}). 
    314322 
    315 %-------------------------------------------------TABLE--------------------------------------------------- 
    316323\begin{table}[htb] 
    317324  \centering 
     
    327334  \label{tab:LDF_files} 
    328335\end{table} 
    329 %-------------------------------------------------------------------------------------------------------------- 
    330  
     336 
     337%% ================================================================================================= 
    331338\subsection[Constant mixing coefficients (\forcode{=0})]{ Constant mixing coefficients (\protect\np[=0]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=0]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    332339 
     
    345352 $U_{scl}$ and $L_{scl}$ are given by the namelist parameters \np{rn_Ud}{rn\_Ud}, \np{rn_Uv}{rn\_Uv}, \np{rn_Ld}{rn\_Ld} and \np{rn_Lv}{rn\_Lv}. 
    346353 
     354%% ================================================================================================= 
    347355\subsection[Vertically varying mixing coefficients (\forcode{=10})]{Vertically varying mixing coefficients (\protect\np[=10]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=10]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    348356 
     
    352360This profile is hard coded in module \mdl{ldfc1d\_c2d}, but can be easily modified by users. 
    353361 
     362%% ================================================================================================= 
    354363\subsection[Mesh size dependent mixing coefficients (\forcode{=20})]{Mesh size dependent mixing coefficients (\protect\np[=20]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=20]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    355364 
     
    377386\colorbox{yellow}{CASE \np{nn_aht_ijk_t}{nn\_aht\_ijk\_t} = 21 to be added} 
    378387 
     388%% ================================================================================================= 
    379389\subsection[Mesh size and depth dependent mixing coefficients (\forcode{=30})]{Mesh size and depth dependent mixing coefficients (\protect\np[=30]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=30]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    380390 
     
    383393the magnitude of the coefficient. 
    384394 
     395%% ================================================================================================= 
    385396\subsection[Velocity dependent mixing coefficients (\forcode{=31})]{Flow dependent mixing coefficients (\protect\np[=31]{nn_aht_ijk_t}{nn\_aht\_ijk\_t} \& \protect\np[=31]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    386397In that case, the eddy coefficient is proportional to the local velocity magnitude so that the Reynolds number $Re =  \lvert U \rvert  e / A_l$ is constant (and here hardcoded to $12$): 
     
    397408\end{equation} 
    398409 
     410%% ================================================================================================= 
    399411\subsection[Deformation rate dependent viscosities (\forcode{nn_ahm_ijk_t=32})]{Deformation rate dependent viscosities (\protect\np[=32]{nn_ahm_ijk_t}{nn\_ahm\_ijk\_t})} 
    400412 
     
    433445where $C_{min}$ and $C_{max}$ are adimensional namelist parameters given by \np{rn_minfac}{rn\_minfac} and \np{rn_maxfac}{rn\_maxfac} respectively. 
    434446 
     447%% ================================================================================================= 
    435448\subsection{About space and time varying mixing coefficients} 
    436449 
     
    446459(\autoref{sec:INVARIANTS_dynldf_properties}). 
    447460 
     461%% ================================================================================================= 
    448462\section[Eddy induced velocity (\forcode{ln_ldfeiv})]{Eddy induced velocity (\protect\np{ln_ldfeiv}{ln\_ldfeiv})} 
    449463 
    450464\label{sec:LDF_eiv} 
    451465 
    452 %--------------------------------------------namtra_eiv--------------------------------------------------- 
    453466 
    454467\begin{listing} 
     
    458471\end{listing} 
    459472 
    460 %-------------------------------------------------------------------------------------------------------------- 
    461473 
    462474%%gm  from Triad appendix  : to be incorporated.... 
     
    513525In case of setting \np[=.true.]{ln_traldf_triad}{ln\_traldf\_triad}, a skew form of the eddy induced advective fluxes is used, which is described in \autoref{apdx:TRIADS}. 
    514526 
     527%% ================================================================================================= 
    515528\section[Mixed layer eddies (\forcode{ln_mle})]{Mixed layer eddies (\protect\np{ln_mle}{ln\_mle})} 
    516529\label{sec:LDF_mle} 
    517530 
    518 %--------------------------------------------namtra_eiv--------------------------------------------------- 
    519531 
    520532\begin{listing} 
     
    524536\end{listing} 
    525537 
    526 %-------------------------------------------------------------------------------------------------------------- 
    527538 
    528539If  \np[=.true.]{ln_mle}{ln\_mle} in \nam{tra_mle}{tra\_mle} namelist, a parameterization of the mixing due to unresolved mixed layer instabilities is activated (\citet{foxkemper.ferrari_JPO08}). Additional transport is computed in \rou{ldf\_mle\_trp} and added to the eulerian transport in \rou{tra\_adv} as done for eddy induced advection. 
Note: See TracChangeset for help on using the changeset viewer.