Changeset 11565
- Timestamp:
- 2019-09-18T16:58:58+02:00 (5 years ago)
- Location:
- NEMO/trunk/doc/latex/NEMO/subfiles
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/subfiles/chap_DYN.tex
r11561 r11565 211 211 % enstrophy conserving scheme 212 212 %------------------------------------------------------------- 213 \subsubsection[Enstrophy conserving scheme (\forcode{ln_dynvor_ens = .true.})]{Enstrophy conserving scheme (\protect\np{ln\_dynvor\_ens}\forcode{ = .true.})}213 \subsubsection[Enstrophy conserving scheme (\forcode{ln_dynvor_ens})]{Enstrophy conserving scheme (\protect\np{ln\_dynvor\_ens})} 214 214 \label{subsec:DYN_vor_ens} 215 215 … … 234 234 % energy conserving scheme 235 235 %------------------------------------------------------------- 236 \subsubsection[Energy conserving scheme (\forcode{ln_dynvor_ene = .true.})]{Energy conserving scheme (\protect\np{ln\_dynvor\_ene}\forcode{ = .true.})}236 \subsubsection[Energy conserving scheme (\forcode{ln_dynvor_ene})]{Energy conserving scheme (\protect\np{ln\_dynvor\_ene})} 237 237 \label{subsec:DYN_vor_ene} 238 238 … … 254 254 % mix energy/enstrophy conserving scheme 255 255 %------------------------------------------------------------- 256 \subsubsection[Mixed energy/enstrophy conserving scheme (\forcode{ln_dynvor_mix = .true.})]{Mixed energy/enstrophy conserving scheme (\protect\np{ln\_dynvor\_mix}\forcode{ = .true.})}256 \subsubsection[Mixed energy/enstrophy conserving scheme (\forcode{ln_dynvor_mix})]{Mixed energy/enstrophy conserving scheme (\protect\np{ln\_dynvor\_mix})} 257 257 \label{subsec:DYN_vor_mix} 258 258 … … 279 279 % energy and enstrophy conserving scheme 280 280 %------------------------------------------------------------- 281 \subsubsection[Energy and enstrophy conserving scheme (\forcode{ln_dynvor_een = .true.})]{Energy and enstrophy conserving scheme (\protect\np{ln\_dynvor\_een}\forcode{ = .true.})}281 \subsubsection[Energy and enstrophy conserving scheme (\forcode{ln_dynvor_een})]{Energy and enstrophy conserving scheme (\protect\np{ln\_dynvor\_een})} 282 282 \label{subsec:DYN_vor_een} 283 283 … … 489 489 % 2nd order centred scheme 490 490 %------------------------------------------------------------- 491 \subsubsection[CEN2: $2^{nd}$ order centred scheme (\forcode{ln_dynadv_cen2 = .true.})]{CEN2: $2^{nd}$ order centred scheme (\protect\np{ln\_dynadv\_cen2}\forcode{ = .true.})}491 \subsubsection[CEN2: $2^{nd}$ order centred scheme (\forcode{ln_dynadv_cen2})]{CEN2: $2^{nd}$ order centred scheme (\protect\np{ln\_dynadv\_cen2})} 492 492 \label{subsec:DYN_adv_cen2} 493 493 … … 512 512 % UBS scheme 513 513 %------------------------------------------------------------- 514 \subsubsection[UBS: Upstream Biased Scheme (\forcode{ln_dynadv_ubs = .true.})]{UBS: Upstream Biased Scheme (\protect\np{ln\_dynadv\_ubs}\forcode{ = .true.})}514 \subsubsection[UBS: Upstream Biased Scheme (\forcode{ln_dynadv_ubs})]{UBS: Upstream Biased Scheme (\protect\np{ln\_dynadv\_ubs})} 515 515 \label{subsec:DYN_adv_ubs} 516 516 … … 591 591 % z-coordinate with full step 592 592 %-------------------------------------------------------------------------------------------------------------- 593 \subsection[Full step $Z$-coordinate (\forcode{ln_dynhpg_zco = .true.})]{Full step $Z$-coordinate (\protect\np{ln\_dynhpg\_zco}\forcode{ = .true.})}593 \subsection[Full step $Z$-coordinate (\forcode{ln_dynhpg_zco})]{Full step $Z$-coordinate (\protect\np{ln\_dynhpg\_zco})} 594 594 \label{subsec:DYN_hpg_zco} 595 595 … … 636 636 % z-coordinate with partial step 637 637 %-------------------------------------------------------------------------------------------------------------- 638 \subsection[Partial step $Z$-coordinate (\forcode{ln_dynhpg_zps = .true.})]{Partial step $Z$-coordinate (\protect\np{ln\_dynhpg\_zps}\forcode{ = .true.})}638 \subsection[Partial step $Z$-coordinate (\forcode{ln_dynhpg_zps})]{Partial step $Z$-coordinate (\protect\np{ln\_dynhpg\_zps})} 639 639 \label{subsec:DYN_hpg_zps} 640 640 … … 722 722 % Time-scheme 723 723 %-------------------------------------------------------------------------------------------------------------- 724 \subsection[Time-scheme (\forcode{ln_dynhpg_imp = .{true,false}.})]{Time-scheme (\protect\np{ln\_dynhpg\_imp}\forcode{ = .\{true,false\}}.)}724 \subsection[Time-scheme (\forcode{ln_dynhpg_imp})]{Time-scheme (\protect\np{ln\_dynhpg\_imp})} 725 725 \label{subsec:DYN_hpg_imp} 726 726 … … 825 825 % Explicit free surface formulation 826 826 %-------------------------------------------------------------------------------------------------------------- 827 \subsection[Explicit free surface (\ texttt{ln\_dynspg\_exp}\forcode{ = .true.})]{Explicit free surface (\protect\np{ln\_dynspg\_exp}\forcode{ = .true.})}827 \subsection[Explicit free surface (\forcode{ln_dynspg_exp})]{Explicit free surface (\protect\np{ln\_dynspg\_exp})} 828 828 \label{subsec:DYN_spg_exp} 829 829 … … 851 851 % Split-explict free surface formulation 852 852 %-------------------------------------------------------------------------------------------------------------- 853 \subsection[Split-explicit free surface (\ texttt{ln\_dynspg\_ts}\forcode{ = .true.})]{Split-explicit free surface (\protect\np{ln\_dynspg\_ts}\forcode{ = .true.})}853 \subsection[Split-explicit free surface (\forcode{ln_dynspg_ts})]{Split-explicit free surface (\protect\np{ln\_dynspg\_ts})} 854 854 \label{subsec:DYN_spg_ts} 855 855 %------------------------------------------namsplit----------------------------------------------------------- … … 1096 1096 % Filtered free surface formulation 1097 1097 %-------------------------------------------------------------------------------------------------------------- 1098 \subsection[Filtered free surface (\ texttt{dynspg\_flt?})]{Filtered free surface (\protect\texttt{dynspg\_flt?})}1098 \subsection[Filtered free surface (\forcode{dynspg_flt?})]{Filtered free surface (\protect\texttt{dynspg\_flt?})} 1099 1099 \label{subsec:DYN_spg_fltp} 1100 1100 … … 1162 1162 1163 1163 % ================================================================ 1164 \subsection[Iso-level laplacian (\forcode{ln_dynldf_lap = .true.})]{Iso-level laplacian operator (\protect\np{ln\_dynldf\_lap}\forcode{ = .true.})}1164 \subsection[Iso-level laplacian (\forcode{ln_dynldf_lap})]{Iso-level laplacian operator (\protect\np{ln\_dynldf\_lap})} 1165 1165 \label{subsec:DYN_ldf_lap} 1166 1166 … … 1187 1187 % Rotated laplacian operator 1188 1188 %-------------------------------------------------------------------------------------------------------------- 1189 \subsection[Rotated laplacian (\forcode{ln_dynldf_iso = .true.})]{Rotated laplacian operator (\protect\np{ln\_dynldf\_iso}\forcode{ = .true.})}1189 \subsection[Rotated laplacian (\forcode{ln_dynldf_iso})]{Rotated laplacian operator (\protect\np{ln\_dynldf\_iso})} 1190 1190 \label{subsec:DYN_ldf_iso} 1191 1191 … … 1245 1245 % Iso-level bilaplacian operator 1246 1246 %-------------------------------------------------------------------------------------------------------------- 1247 \subsection[Iso-level bilaplacian (\forcode{ln_dynldf_bilap = .true.})]{Iso-level bilaplacian operator (\protect\np{ln\_dynldf\_bilap}\forcode{ = .true.})}1247 \subsection[Iso-level bilaplacian (\forcode{ln_dynldf_bilap})]{Iso-level bilaplacian operator (\protect\np{ln\_dynldf\_bilap})} 1248 1248 \label{subsec:DYN_ldf_bilap} 1249 1249 -
NEMO/trunk/doc/latex/NEMO/subfiles/chap_LDF.tex
r11561 r11565 34 34 % Lateral Mixing Operator 35 35 % ================================================================ 36 \section[Lateral mixing operators] 37 {Lateral mixing operators} 36 \section[Lateral mixing operators]{Lateral mixing operators} 38 37 \label{sec:LDF_op} 39 38 We 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}. 40 39 41 \subsection[No lateral mixing (\forcode{ln_traldf_OFF}, \forcode{ln_dynldf_OFF})] 42 {No lateral mixing (\protect\np{ln\_traldf\_OFF}, \protect\np{ln\_dynldf\_OFF})} 40 \subsection[No lateral mixing (\forcode{ln_traldf_OFF} \& \forcode{ln_dynldf_OFF})]{No lateral mixing (\protect\np{ln\_traldf\_OFF} \& \protect\np{ln\_dynldf\_OFF})} 43 41 44 42 It is possible to run without explicit lateral diffusion on tracers (\protect\np{ln\_traldf\_OFF}\forcode{=.true.}) and/or … … 47 45 see \autoref{subsec:DYN_adv_ubs}) and can be useful for testing purposes. 48 46 49 \subsection[Laplacian mixing (\forcode{ln_traldf_lap}, \forcode{ln_dynldf_lap})] 50 {Laplacian mixing (\protect\np{ln\_traldf\_lap}, \protect\np{ln\_dynldf\_lap})} 47 \subsection[Laplacian mixing (\forcode{ln_traldf_lap} \& \forcode{ln_dynldf_lap})]{Laplacian mixing (\protect\np{ln\_traldf\_lap} \& \protect\np{ln\_dynldf\_lap})} 51 48 Setting \protect\np{ln\_traldf\_lap}\forcode{=.true.} and/or \protect\np{ln\_dynldf\_lap}\forcode{=.true.} enables 52 49 a second order diffusion on tracers and momentum respectively. Note that in \NEMO\ 4, one can not combine 53 50 Laplacian and Bilaplacian operators for the same variable. 54 51 55 \subsection[Bilaplacian mixing (\forcode{ln_traldf_blp}, \forcode{ln_dynldf_blp})] 56 {Bilaplacian mixing (\protect\np{ln\_traldf\_blp}, \protect\np{ln\_dynldf\_blp})} 52 \subsection[Bilaplacian mixing (\forcode{ln_traldf_blp} \& \forcode{ln_dynldf_blp})]{Bilaplacian mixing (\protect\np{ln\_traldf\_blp} \& \protect\np{ln\_dynldf\_blp})} 57 53 Setting \protect\np{ln\_traldf\_blp}\forcode{=.true.} and/or \protect\np{ln\_dynldf\_blp}\forcode{=.true.} enables 58 54 a fourth order diffusion on tracers and momentum respectively. It is implemented by calling the above Laplacian operator twice. … … 62 58 % Direction of lateral Mixing 63 59 % ================================================================ 64 \section[Direction of lateral mixing (\textit{ldfslp.F90})] 65 {Direction of lateral mixing (\protect\mdl{ldfslp})} 60 \section[Direction of lateral mixing (\textit{ldfslp.F90})]{Direction of lateral mixing (\protect\mdl{ldfslp})} 66 61 \label{sec:LDF_slp} 67 62 … … 325 320 % Lateral Mixing Coefficients 326 321 % ================================================================ 327 \section[Lateral mixing coefficient (\forcode{nn_aht_ijk_t}, \forcode{nn_ahm_ijk_t})] 328 {Lateral mixing coefficient (\protect\np{nn\_aht\_ijk\_t}, \protect\np{nn\_ahm\_ijk\_t})} 322 \section[Lateral mixing coefficient (\forcode{nn_aht_ijk_t} \& \forcode{nn_ahm_ijk_t})]{Lateral mixing coefficient (\protect\np{nn\_aht\_ijk\_t} \& \protect\np{nn\_ahm\_ijk\_t})} 329 323 \label{sec:LDF_coef} 330 324 … … 332 326 The way the mixing coefficients are set in the reference version can be described as follows: 333 327 334 \subsection[Mixing coefficients read from file (\forcode{nn_aht_ijk_t=-20, -30}, \forcode{nn_ahm_ijk_t=-20,-30})] 335 { Mixing coefficients read from file (\protect\np{nn\_aht\_ijk\_t}\forcode{=-20, -30}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=-20, -30})} 328 \subsection[Mixing coefficients read from file (\forcode{=-20, -30})]{ Mixing coefficients read from file (\protect\np{nn\_aht\_ijk\_t}\forcode{=-20, -30} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=-20, -30})} 336 329 337 330 Mixing coefficients can be read from file if a particular geographical variation is needed. For example, in the ORCA2 global ocean model, … … 357 350 %-------------------------------------------------------------------------------------------------------------- 358 351 359 \subsection[Constant mixing coefficients (\forcode{nn_aht_ijk_t=0}, \forcode{nn_ahm_ijk_t=0})] 360 { Constant mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=0}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=0})} 352 \subsection[Constant mixing coefficients (\forcode{=0})]{ Constant mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=0} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=0})} 361 353 362 354 If constant, mixing coefficients are set thanks to a velocity and a length scales ($U_{scl}$, $L_{scl}$) such that: … … 374 366 $U_{scl}$ and $L_{scl}$ are given by the namelist parameters \np{rn\_Ud}, \np{rn\_Uv}, \np{rn\_Ld} and \np{rn\_Lv}. 375 367 376 \subsection[Vertically varying mixing coefficients (\forcode{nn_aht_ijk_t=10}, \forcode{nn_ahm_ijk_t=10})] 377 {Vertically varying mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=10}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=10})} 368 \subsection[Vertically varying mixing coefficients (\forcode{=10})]{Vertically varying mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=10} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=10})} 378 369 379 370 In the vertically varying case, a hyperbolic variation of the lateral mixing coefficient is introduced in which … … 382 373 This profile is hard coded in module \mdl{ldfc1d\_c2d}, but can be easily modified by users. 383 374 384 \subsection[Mesh size dependent mixing coefficients (\forcode{nn_aht_ijk_t=20}, \forcode{nn_ahm_ijk_t=20})] 385 {Mesh size dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=20}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=20})} 375 \subsection[Mesh size dependent mixing coefficients (\forcode{=20})]{Mesh size dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=20} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=20})} 386 376 387 377 In that case, the horizontal variation of the eddy coefficient depends on the local mesh size and … … 408 398 \colorbox{yellow}{CASE \np{nn\_aht\_ijk\_t} = 21 to be added} 409 399 410 \subsection[Mesh size and depth dependent mixing coefficients (\forcode{nn_aht_ijk_t=30}, \forcode{nn_ahm_ijk_t=30})] 411 {Mesh size and depth dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=30}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=30})} 400 \subsection[Mesh size and depth dependent mixing coefficients (\forcode{=30})]{Mesh size and depth dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=30} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=30})} 412 401 413 402 The 3D space variation of the mixing coefficient is simply the combination of the 1D and 2D cases above, … … 415 404 the magnitude of the coefficient. 416 405 417 \subsection[Velocity dependent mixing coefficients (\forcode{nn_aht_ijk_t=31}, \forcode{nn_ahm_ijk_t=31})] 418 {Flow dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=31}, \protect\np{nn\_ahm\_ijk\_t}\forcode{=31})} 406 \subsection[Velocity dependent mixing coefficients (\forcode{=31})]{Flow dependent mixing coefficients (\protect\np{nn\_aht\_ijk\_t}\forcode{=31} \& \protect\np{nn\_ahm\_ijk\_t}\forcode{=31})} 419 407 In 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$): 420 408 \colorbox{yellow}{JC comment: The Reynolds is effectively set to 12 in the code for both operators but shouldn't it be 2 for Laplacian ?} … … 430 418 \end{equation} 431 419 432 \subsection[Deformation rate dependent viscosities (\forcode{nn_ahm_ijk_t=32})] 433 {Deformation rate dependent viscosities (\protect\np{nn\_ahm\_ijk\_t}\forcode{=32})} 420 \subsection[Deformation rate dependent viscosities (\forcode{nn_ahm_ijk_t=32})]{Deformation rate dependent viscosities (\protect\np{nn\_ahm\_ijk\_t}\forcode{=32})} 434 421 435 422 This option refers to the \citep{smagorinsky_MW63} scheme which is here implemented for momentum only. Smagorinsky chose as a … … 483 470 % Eddy Induced Mixing 484 471 % ================================================================ 485 \section[Eddy induced velocity (\forcode{ln_ldfeiv=.true.})] 486 {Eddy induced velocity (\protect\np{ln\_ldfeiv}\forcode{=.true.})} 472 \section[Eddy induced velocity (\forcode{ln_ldfeiv})]{Eddy induced velocity (\protect\np{ln\_ldfeiv})} 487 473 488 474 \label{sec:LDF_eiv} … … 555 541 % Mixed layer eddies 556 542 % ================================================================ 557 \section[Mixed layer eddies (\forcode{ln_mle=.true.})] 558 {Mixed layer eddies (\protect\np{ln\_mle}\forcode{=.true.})} 559 543 \section[Mixed layer eddies (\forcode{ln_mle})]{Mixed layer eddies (\protect\np{ln\_mle})} 560 544 \label{sec:LDF_mle} 561 545 -
NEMO/trunk/doc/latex/NEMO/subfiles/chap_TRA.tex
r11561 r11565 60 60 % Tracer Advection 61 61 % ================================================================ 62 \section[Tracer advection (\textit{traadv.F90})] 63 {Tracer advection (\protect\mdl{traadv})} 62 \section[Tracer advection (\textit{traadv.F90})]{Tracer advection (\protect\mdl{traadv})} 64 63 \label{sec:TRA_adv} 65 64 %------------------------------------------namtra_adv----------------------------------------------------- … … 186 185 % 2nd and 4th order centred schemes 187 186 % ------------------------------------------------------------------------------------------------------------- 188 \subsection[CEN: Centred scheme (\forcode{ln_traadv_cen=.true.})] 189 {CEN: Centred scheme (\protect\np{ln\_traadv\_cen}\forcode{=.true.})} 187 \subsection[CEN: Centred scheme (\forcode{ln_traadv_cen})]{CEN: Centred scheme (\protect\np{ln\_traadv\_cen})} 190 188 \label{subsec:TRA_adv_cen} 191 189 … … 254 252 % FCT scheme 255 253 % ------------------------------------------------------------------------------------------------------------- 256 \subsection[FCT: Flux Corrected Transport scheme (\forcode{ln_traadv_fct=.true.})] 257 {FCT: Flux Corrected Transport scheme (\protect\np{ln\_traadv\_fct}\forcode{=.true.})} 254 \subsection[FCT: Flux Corrected Transport scheme (\forcode{ln_traadv_fct})]{FCT: Flux Corrected Transport scheme (\protect\np{ln\_traadv\_fct})} 258 255 \label{subsec:TRA_adv_tvd} 259 256 … … 298 295 % MUSCL scheme 299 296 % ------------------------------------------------------------------------------------------------------------- 300 \subsection[MUSCL: Monotone Upstream Scheme for Conservative Laws (\forcode{ln_traadv_mus=.true.})] 301 {MUSCL: Monotone Upstream Scheme for Conservative Laws (\protect\np{ln\_traadv\_mus}\forcode{=.true.})} 297 \subsection[MUSCL: Monotone Upstream Scheme for Conservative Laws (\forcode{ln_traadv_mus})]{MUSCL: Monotone Upstream Scheme for Conservative Laws (\protect\np{ln\_traadv\_mus})} 302 298 \label{subsec:TRA_adv_mus} 303 299 … … 335 331 % UBS scheme 336 332 % ------------------------------------------------------------------------------------------------------------- 337 \subsection[UBS a.k.a. UP3: Upstream-Biased Scheme (\forcode{ln_traadv_ubs=.true.})] 338 {UBS a.k.a. UP3: Upstream-Biased Scheme (\protect\np{ln\_traadv\_ubs}\forcode{=.true.})} 333 \subsection[UBS a.k.a. UP3: Upstream-Biased Scheme (\forcode{ln_traadv_ubs})]{UBS a.k.a. UP3: Upstream-Biased Scheme (\protect\np{ln\_traadv\_ubs})} 339 334 \label{subsec:TRA_adv_ubs} 340 335 … … 408 403 % QCK scheme 409 404 % ------------------------------------------------------------------------------------------------------------- 410 \subsection[QCK: QuiCKest scheme (\forcode{ln_traadv_qck=.true.})] 411 {QCK: QuiCKest scheme (\protect\np{ln\_traadv\_qck}\forcode{=.true.})} 405 \subsection[QCK: QuiCKest scheme (\forcode{ln_traadv_qck})]{QCK: QuiCKest scheme (\protect\np{ln\_traadv\_qck})} 412 406 \label{subsec:TRA_adv_qck} 413 407 … … 432 426 % Tracer Lateral Diffusion 433 427 % ================================================================ 434 \section[Tracer lateral diffusion (\textit{traldf.F90})] 435 {Tracer lateral diffusion (\protect\mdl{traldf})} 428 \section[Tracer lateral diffusion (\textit{traldf.F90})]{Tracer lateral diffusion (\protect\mdl{traldf})} 436 429 \label{sec:TRA_ldf} 437 430 %-----------------------------------------nam_traldf------------------------------------------------------ … … 465 458 % Type of operator 466 459 % ------------------------------------------------------------------------------------------------------------- 467 \subsection[Type of operator (\texttt{ln\_traldf}\{\texttt{\_OFF,\_lap,\_blp}\})] 468 {Type of operator (\protect\np{ln\_traldf\_OFF}, \protect\np{ln\_traldf\_lap}, or \protect\np{ln\_traldf\_blp}) } 460 \subsection[Type of operator (\texttt{ln\_traldf}\{\texttt{\_OFF,\_lap,\_blp}\})]{Type of operator (\protect\np{ln\_traldf\_OFF}, \protect\np{ln\_traldf\_lap}, or \protect\np{ln\_traldf\_blp})} 469 461 \label{subsec:TRA_ldf_op} 470 462 … … 500 492 % Direction of action 501 493 % ------------------------------------------------------------------------------------------------------------- 502 \subsection[Action direction (\texttt{ln\_traldf}\{\texttt{\_lev,\_hor,\_iso,\_triad}\})] 503 {Direction of action (\protect\np{ln\_traldf\_lev}, \protect\np{ln\_traldf\_hor}, \protect\np{ln\_traldf\_iso}, or \protect\np{ln\_traldf\_triad}) } 494 \subsection[Action direction (\texttt{ln\_traldf}\{\texttt{\_lev,\_hor,\_iso,\_triad}\})]{Direction of action (\protect\np{ln\_traldf\_lev}, \protect\np{ln\_traldf\_hor}, \protect\np{ln\_traldf\_iso}, or \protect\np{ln\_traldf\_triad})} 504 495 \label{subsec:TRA_ldf_dir} 505 496 … … 527 518 % iso-level operator 528 519 % ------------------------------------------------------------------------------------------------------------- 529 \subsection[Iso-level (bi-)laplacian operator (\texttt{ln\_traldf\_iso})] 530 {Iso-level (bi-)laplacian operator ( \protect\np{ln\_traldf\_iso})} 520 \subsection[Iso-level (bi-)laplacian operator (\texttt{ln\_traldf\_iso})]{Iso-level (bi-)laplacian operator ( \protect\np{ln\_traldf\_iso})} 531 521 \label{subsec:TRA_ldf_lev} 532 522 … … 564 554 %&& Standard rotated (bi-)laplacian operator 565 555 %&& ---------------------------------------------- 566 \subsubsection[Standard rotated (bi-)laplacian operator (\textit{traldf\_iso.F90})] 567 {Standard rotated (bi-)laplacian operator (\protect\mdl{traldf\_iso})} 556 \subsubsection[Standard rotated (bi-)laplacian operator (\textit{traldf\_iso.F90})]{Standard rotated (bi-)laplacian operator (\protect\mdl{traldf\_iso})} 568 557 \label{subsec:TRA_ldf_iso} 569 558 The general form of the second order lateral tracer subgrid scale physics (\autoref{eq:MB_zdf}) … … 608 597 %&& Triad rotated (bi-)laplacian operator 609 598 %&& ------------------------------------------- 610 \subsubsection[Triad rotated (bi-)laplacian operator (\textit{ln\_traldf\_triad})] 611 {Triad rotated (bi-)laplacian operator (\protect\np{ln\_traldf\_triad})} 599 \subsubsection[Triad rotated (bi-)laplacian operator (\textit{ln\_traldf\_triad})]{Triad rotated (bi-)laplacian operator (\protect\np{ln\_traldf\_triad})} 612 600 \label{subsec:TRA_ldf_triad} 613 601 … … 641 629 % Tracer Vertical Diffusion 642 630 % ================================================================ 643 \section[Tracer vertical diffusion (\textit{trazdf.F90})] 644 {Tracer vertical diffusion (\protect\mdl{trazdf})} 631 \section[Tracer vertical diffusion (\textit{trazdf.F90})]{Tracer vertical diffusion (\protect\mdl{trazdf})} 645 632 \label{sec:TRA_zdf} 646 633 %--------------------------------------------namzdf--------------------------------------------------------- … … 685 672 % surface boundary condition 686 673 % ------------------------------------------------------------------------------------------------------------- 687 \subsection[Surface boundary condition (\textit{trasbc.F90})] 688 {Surface boundary condition (\protect\mdl{trasbc})} 674 \subsection[Surface boundary condition (\textit{trasbc.F90})]{Surface boundary condition (\protect\mdl{trasbc})} 689 675 \label{subsec:TRA_sbc} 690 676 … … 759 745 % Solar Radiation Penetration 760 746 % ------------------------------------------------------------------------------------------------------------- 761 \subsection[Solar radiation penetration (\textit{traqsr.F90})] 762 {Solar radiation penetration (\protect\mdl{traqsr})} 747 \subsection[Solar radiation penetration (\textit{traqsr.F90})]{Solar radiation penetration (\protect\mdl{traqsr})} 763 748 \label{subsec:TRA_qsr} 764 749 %--------------------------------------------namqsr-------------------------------------------------------- … … 883 868 % Bottom Boundary Condition 884 869 % ------------------------------------------------------------------------------------------------------------- 885 \subsection[Bottom boundary condition (\textit{trabbc.F90}) - \forcode{ln_trabbc=.true.})] 886 {Bottom boundary condition (\protect\mdl{trabbc})} 870 \subsection[Bottom boundary condition (\textit{trabbc.F90}) - \forcode{ln_trabbc})]{Bottom boundary condition (\protect\mdl{trabbc} - \protect\np{ln\_trabbc})} 887 871 \label{subsec:TRA_bbc} 888 872 %--------------------------------------------nambbc-------------------------------------------------------- … … 924 908 % Bottom Boundary Layer 925 909 % ================================================================ 926 \section[Bottom boundary layer (\textit{trabbl.F90} - \forcode{ln_trabbl=.true.})] 927 {Bottom boundary layer (\protect\mdl{trabbl} - \protect\np{ln\_trabbl}\forcode{=.true.})} 910 \section[Bottom boundary layer (\textit{trabbl.F90} - \forcode{ln_trabbl})]{Bottom boundary layer (\protect\mdl{trabbl} - \protect\np{ln\_trabbl})} 928 911 \label{sec:TRA_bbl} 929 912 %--------------------------------------------nambbl--------------------------------------------------------- … … 961 944 % Diffusive BBL 962 945 % ------------------------------------------------------------------------------------------------------------- 963 \subsection[Diffusive bottom boundary layer (\forcode{nn_bbl_ldf=1})] 964 {Diffusive bottom boundary layer (\protect\np{nn\_bbl\_ldf}\forcode{=1})} 946 \subsection[Diffusive bottom boundary layer (\forcode{nn_bbl_ldf=1})]{Diffusive bottom boundary layer (\protect\np{nn\_bbl\_ldf}\forcode{=1})} 965 947 \label{subsec:TRA_bbl_diff} 966 948 … … 1001 983 % Advective BBL 1002 984 % ------------------------------------------------------------------------------------------------------------- 1003 \subsection[Advective bottom boundary layer (\forcode{nn_bbl_adv=[12]})] 1004 {Advective bottom boundary layer (\protect\np{nn\_bbl\_adv}\forcode{=[12]})} 985 \subsection[Advective bottom boundary layer (\forcode{nn_bbl_adv=[12]})]{Advective bottom boundary layer (\protect\np{nn\_bbl\_adv}\forcode{=[12]})} 1005 986 \label{subsec:TRA_bbl_adv} 1006 987 … … 1093 1074 % Tracer damping 1094 1075 % ================================================================ 1095 \section[Tracer damping (\textit{tradmp.F90})] 1096 {Tracer damping (\protect\mdl{tradmp})} 1076 \section[Tracer damping (\textit{tradmp.F90})]{Tracer damping (\protect\mdl{tradmp})} 1097 1077 \label{sec:TRA_dmp} 1098 1078 %--------------------------------------------namtra_dmp------------------------------------------------- … … 1153 1133 % Tracer time evolution 1154 1134 % ================================================================ 1155 \section[Tracer time evolution (\textit{tranxt.F90})] 1156 {Tracer time evolution (\protect\mdl{tranxt})} 1135 \section[Tracer time evolution (\textit{tranxt.F90})]{Tracer time evolution (\protect\mdl{tranxt})} 1157 1136 \label{sec:TRA_nxt} 1158 1137 %--------------------------------------------namdom----------------------------------------------------- … … 1189 1168 % Equation of State (eosbn2) 1190 1169 % ================================================================ 1191 \section[Equation of state (\textit{eosbn2.F90})] 1192 {Equation of state (\protect\mdl{eosbn2})} 1170 \section[Equation of state (\textit{eosbn2.F90})]{Equation of state (\protect\mdl{eosbn2})} 1193 1171 \label{sec:TRA_eosbn2} 1194 1172 %--------------------------------------------nameos----------------------------------------------------- … … 1204 1182 % Equation of State 1205 1183 % ------------------------------------------------------------------------------------------------------------- 1206 \subsection[Equation of seawater (\texttt{ln}\{\texttt{\_teso10,\_eos80,\_seos}\})] 1207 {Equation of seawater (\protect\np{ln\_teos10}, \protect\np{ln\_teos80}, or \protect\np{ln\_seos}) } 1184 \subsection[Equation of seawater (\texttt{ln}\{\texttt{\_teso10,\_eos80,\_seos}\})]{Equation of seawater (\protect\np{ln\_teos10}, \protect\np{ln\_teos80}, or \protect\np{ln\_seos})} 1208 1185 \label{subsec:TRA_eos} 1209 1186 … … 1330 1307 % Brunt-V\"{a}is\"{a}l\"{a} Frequency 1331 1308 % ------------------------------------------------------------------------------------------------------------- 1332 \subsection[Brunt-V\"{a}is\"{a}l\"{a} frequency] 1333 {Brunt-V\"{a}is\"{a}l\"{a} frequency} 1309 \subsection[Brunt-V\"{a}is\"{a}l\"{a} frequency]{Brunt-V\"{a}is\"{a}l\"{a} frequency} 1334 1310 \label{subsec:TRA_bn2} 1335 1311 … … 1385 1361 % Horizontal Derivative in zps-coordinate 1386 1362 % ================================================================ 1387 \section[Horizontal derivative in \textit{zps}-coordinate (\textit{zpshde.F90})] 1388 {Horizontal derivative in \textit{zps}-coordinate (\protect\mdl{zpshde})} 1363 \section[Horizontal derivative in \textit{zps}-coordinate (\textit{zpshde.F90})]{Horizontal derivative in \textit{zps}-coordinate (\protect\mdl{zpshde})} 1389 1364 \label{sec:TRA_zpshde} 1390 1365 -
NEMO/trunk/doc/latex/NEMO/subfiles/chap_ZDF.tex
r11561 r11565 58 58 % Constant 59 59 % ------------------------------------------------------------------------------------------------------------- 60 \subsection[Constant (\forcode{ln_zdfcst=.true.})] 61 {Constant (\protect\np{ln\_zdfcst}\forcode{=.true.})} 60 \subsection[Constant (\forcode{ln_zdfcst})]{Constant (\protect\np{ln\_zdfcst})} 62 61 \label{subsec:ZDF_cst} 63 62 … … 81 80 % Richardson Number Dependent 82 81 % ------------------------------------------------------------------------------------------------------------- 83 \subsection[Richardson number dependent (\forcode{ln_zdfric=.true.})] 84 {Richardson number dependent (\protect\np{ln\_zdfric}\forcode{=.true.})} 82 \subsection[Richardson number dependent (\forcode{ln_zdfric})]{Richardson number dependent (\protect\np{ln\_zdfric})} 85 83 \label{subsec:ZDF_ric} 86 84 … … 143 141 % TKE Turbulent Closure Scheme 144 142 % ------------------------------------------------------------------------------------------------------------- 145 \subsection[TKE turbulent closure scheme (\forcode{ln_zdftke=.true.})] 146 {TKE turbulent closure scheme (\protect\np{ln\_zdftke}\forcode{=.true.})} 143 \subsection[TKE turbulent closure scheme (\forcode{ln_zdftke})]{TKE turbulent closure scheme (\protect\np{ln\_zdftke})} 147 144 \label{subsec:ZDF_tke} 148 145 %--------------------------------------------namzdf_tke-------------------------------------------------- … … 427 424 % GLS Generic Length Scale Scheme 428 425 % ------------------------------------------------------------------------------------------------------------- 429 \subsection[GLS: Generic Length Scale (\forcode{ln_zdfgls=.true.})] 430 {GLS: Generic Length Scale (\protect\np{ln\_zdfgls}\forcode{=.true.})} 426 \subsection[GLS: Generic Length Scale (\forcode{ln_zdfgls})]{GLS: Generic Length Scale (\protect\np{ln\_zdfgls})} 431 427 \label{subsec:ZDF_gls} 432 428 … … 552 548 % OSM OSMOSIS BL Scheme 553 549 % ------------------------------------------------------------------------------------------------------------- 554 \subsection[OSM: OSMosis boundary layer scheme (\forcode{ln_zdfosm=.true.})] 555 {OSM: OSMosis boundary layer scheme (\protect\np{ln\_zdfosm}\forcode{=.true.})} 550 \subsection[OSM: OSMosis boundary layer scheme (\forcode{ln_zdfosm})]{OSM: OSMosis boundary layer scheme (\protect\np{ln\_zdfosm})} 556 551 \label{subsec:ZDF_osm} 557 552 %--------------------------------------------namzdf_osm--------------------------------------------------------- … … 569 564 % TKE and GLS discretization considerations 570 565 % ------------------------------------------------------------------------------------------------------------- 571 \subsection[ Discrete energy conservation for TKE and GLS schemes] 572 {Discrete energy conservation for TKE and GLS schemes} 566 \subsection[ Discrete energy conservation for TKE and GLS schemes]{Discrete energy conservation for TKE and GLS schemes} 573 567 \label{subsec:ZDF_tke_ene} 574 568 … … 688 682 % Non-Penetrative Convective Adjustment 689 683 % ------------------------------------------------------------------------------------------------------------- 690 \subsection[Non-penetrative convective adjustment (\forcode{ln_tranpc=.true.})] 691 {Non-penetrative convective adjustment (\protect\np{ln\_tranpc}\forcode{=.true.})} 684 \subsection[Non-penetrative convective adjustment (\forcode{ln_tranpc})]{Non-penetrative convective adjustment (\protect\np{ln\_tranpc})} 692 685 \label{subsec:ZDF_npc} 693 686 … … 754 747 % Enhanced Vertical Diffusion 755 748 % ------------------------------------------------------------------------------------------------------------- 756 \subsection[Enhanced vertical diffusion (\forcode{ln_zdfevd=.true.})] 757 {Enhanced vertical diffusion (\protect\np{ln\_zdfevd}\forcode{=.true.})} 749 \subsection[Enhanced vertical diffusion (\forcode{ln_zdfevd})]{Enhanced vertical diffusion (\protect\np{ln\_zdfevd})} 758 750 \label{subsec:ZDF_evd} 759 751 … … 781 773 % Turbulent Closure Scheme 782 774 % ------------------------------------------------------------------------------------------------------------- 783 \subsection {Handling convection with turbulent closure schemes (\forcode{ln_zdf{tke,gls,osm}=.true.})}775 \subsection[Handling convection with turbulent closure schemes (\forcode{ln_zdf{tke,gls,osm}})]{Handling convection with turbulent closure schemes (\forcode{ln_zdf{tke,gls,osm}})} 784 776 \label{subsec:ZDF_tcs} 785 777 … … 809 801 % Double Diffusion Mixing 810 802 % ================================================================ 811 \section[Double diffusion mixing (\forcode{ln_zdfddm=.true.})] 812 {Double diffusion mixing (\protect\np{ln\_zdfddm}\forcode{=.true.})} 803 \section[Double diffusion mixing (\forcode{ln_zdfddm})]{Double diffusion mixing (\protect\np{ln\_zdfddm})} 813 804 \label{subsec:ZDF_ddm} 814 805 … … 905 896 % Bottom Friction 906 897 % ================================================================ 907 \section[Bottom and top friction (\textit{zdfdrg.F90})] 908 {Bottom and top friction (\protect\mdl{zdfdrg})} 909 \label{sec:ZDF_drg} 898 \section[Bottom and top friction (\textit{zdfdrg.F90})] {Bottom and top friction (\protect\mdl{zdfdrg})} 899 \label{sec:ZDF_drg} 910 900 911 901 %--------------------------------------------namdrg-------------------------------------------------------- … … 986 976 % Linear Bottom Friction 987 977 % ------------------------------------------------------------------------------------------------------------- 988 \subsection[Linear top/bottom friction (\forcode{ln_lin=.true.})] 989 {Linear top/bottom friction (\protect\np{ln\_lin}\forcode{=.true.)}} 990 \label{subsec:ZDF_drg_linear} 978 \subsection[Linear top/bottom friction (\forcode{ln_lin})]{Linear top/bottom friction (\protect\np{ln\_lin})} 979 \label{subsec:ZDF_drg_linear} 991 980 992 981 The linear friction parameterisation (including the special case of a free-slip condition) assumes that … … 1026 1015 % Non-Linear Bottom Friction 1027 1016 % ------------------------------------------------------------------------------------------------------------- 1028 \subsection[Non-linear top/bottom friction (\forcode{ln_non_lin=.true.})] 1029 {Non-linear top/bottom friction (\protect\np{ln\_non\_lin}\forcode{=.true.})} 1030 \label{subsec:ZDF_drg_nonlinear} 1017 \subsection[Non-linear top/bottom friction (\forcode{ln_non_lin})]{Non-linear top/bottom friction (\protect\np{ln\_non\_lin})} 1018 \label{subsec:ZDF_drg_nonlinear} 1031 1019 1032 1020 The non-linear bottom friction parameterisation assumes that the top/bottom friction is quadratic: … … 1062 1050 % Bottom Friction Log-layer 1063 1051 % ------------------------------------------------------------------------------------------------------------- 1064 \subsection[Log-layer top/bottom friction (\forcode{ln_loglayer=.true.})] 1065 {Log-layer top/bottom friction (\protect\np{ln\_loglayer}\forcode{=.true.})} 1066 \label{subsec:ZDF_drg_loglayer} 1052 \subsection[Log-layer top/bottom friction (\forcode{ln_loglayer})]{Log-layer top/bottom friction (\protect\np{ln\_loglayer})} 1053 \label{subsec:ZDF_drg_loglayer} 1067 1054 1068 1055 In the non-linear friction case, the drag coefficient, $C_D$, can be optionally enhanced using … … 1089 1076 % Explicit bottom Friction 1090 1077 % ------------------------------------------------------------------------------------------------------------- 1091 \subsection{Explicit top/bottom friction (\forcode{ln_drgimp=.false.})}1092 1078 \subsection[Explicit top/bottom friction (\forcode{ln_drgimp=.false.})]{Explicit top/bottom friction (\protect\np{ln\_drgimp}\forcode{=.false.})} 1079 \label{subsec:ZDF_drg_stability} 1093 1080 1094 1081 Setting \np{ln\_drgimp} \forcode{= .false.} means that bottom friction is treated explicitly in time, which has the advantage of simplifying the interaction with the split-explicit free surface (see \autoref{subsec:ZDF_drg_ts}). The latter does indeed require the knowledge of bottom stresses in the course of the barotropic sub-iteration, which becomes less straightforward in the implicit case. In the explicit case, top/bottom stresses can be computed using \textit{before} velocities and inserted in the overall momentum tendency budget. This reads: … … 1150 1137 % Implicit Bottom Friction 1151 1138 % ------------------------------------------------------------------------------------------------------------- 1152 \subsection[Implicit top/bottom friction (\forcode{ln_drgimp=.true.})] 1153 {Implicit top/bottom friction (\protect\np{ln\_drgimp}\forcode{=.true.})} 1154 \label{subsec:ZDF_drg_imp} 1139 \subsection[Implicit top/bottom friction (\forcode{ln_drgimp=.true.})]{Implicit top/bottom friction (\protect\np{ln\_drgimp}\forcode{=.true.})} 1140 \label{subsec:ZDF_drg_imp} 1155 1141 1156 1142 An optional implicit form of bottom friction has been implemented to improve model stability. … … 1181 1167 % Bottom Friction with split-explicit free surface 1182 1168 % ------------------------------------------------------------------------------------------------------------- 1183 \subsection[Bottom friction with split-explicit free surface] 1184 {Bottom friction with split-explicit free surface} 1185 \label{subsec:ZDF_drg_ts} 1169 \subsection[Bottom friction with split-explicit free surface]{Bottom friction with split-explicit free surface} 1170 \label{subsec:ZDF_drg_ts} 1186 1171 1187 1172 With split-explicit free surface, the sub-stepping of barotropic equations needs the knowledge of top/bottom stresses. An obvious way to satisfy this is to take them as constant over the course of the barotropic integration and equal to the value used to update the baroclinic momentum trend. Provided \np{ln\_drgimp}\forcode{= .false.} and a centred or \textit{leap-frog} like integration of barotropic equations is used (\ie\ \forcode{ln_bt_fw=.false.}, cf \autoref{subsec:DYN_spg_ts}), this does ensure that barotropic and baroclinic dynamics feel the same stresses during one leapfrog time step. However, if \np{ln\_drgimp}\forcode{= .true.}, stresses depend on the \textit{after} value of the velocities which themselves depend on the barotropic iteration result. This cyclic dependency makes difficult obtaining consistent stresses in 2d and 3d dynamics. Part of this mismatch is then removed when setting the final barotropic component of 3d velocities to the time splitting estimate. This last step can be seen as a necessary evil but should be minimized since it interferes with the adjustment to the boundary conditions. … … 1199 1184 % Internal wave-driven mixing 1200 1185 % ================================================================ 1201 \section[Internal wave-driven mixing (\forcode{ln_zdfiwm=.true.})] 1202 {Internal wave-driven mixing (\protect\np{ln\_zdfiwm}\forcode{=.true.})} 1186 \section[Internal wave-driven mixing (\forcode{ln_zdfiwm})]{Internal wave-driven mixing (\protect\np{ln\_zdfiwm})} 1203 1187 \label{subsec:ZDF_tmx_new} 1204 1188 … … 1264 1248 % surface wave-induced mixing 1265 1249 % ================================================================ 1266 \section[Surface wave-induced mixing (\forcode{ln_zdfswm=.true.})] 1267 {Surface wave-induced mixing (\protect\np{ln\_zdfswm}\forcode{=.true.})} 1250 \section[Surface wave-induced mixing (\forcode{ln_zdfswm})]{Surface wave-induced mixing (\protect\np{ln\_zdfswm})} 1268 1251 \label{subsec:ZDF_swm} 1269 1252 … … 1298 1281 % Adaptive-implicit vertical advection 1299 1282 % ================================================================ 1300 \section[Adaptive-implicit vertical advection (\forcode{ln_zad_Aimp=.true.})] 1301 {Adaptive-implicit vertical advection(\protect\np{ln\_zad\_Aimp}\forcode{=.true.})} 1283 \section[Adaptive-implicit vertical advection (\forcode{ln_zad_Aimp})]{Adaptive-implicit vertical advection(\protect\np{ln\_zad\_Aimp})} 1302 1284 \label{subsec:ZDF_aimp} 1303 1285
Note: See TracChangeset
for help on using the changeset viewer.