Changeset 11596 for NEMO/trunk/doc/latex/NEMO/subfiles/apdx_triads.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/apdx_triads.tex
r11584 r11596 12 12 13 13 \begin{document} 14 % ================================================================15 % Iso-neutral diffusion :16 % ================================================================17 14 \chapter{Iso-Neutral Diffusion and Eddy Advection using Triads} 18 15 \label{apdx:TRIADS} 19 16 20 17 \chaptertoc 21 22 \newpage23 18 24 19 \section[Choice of \forcode{namtra\_ldf} namelist parameters]{Choice of \protect\nam{tra_ldf}{tra\_ldf} namelist parameters} … … 42 37 The options specific to the Griffies scheme include: 43 38 \begin{description} 44 \item [{\np{ln_triad_iso}{ln\_triad\_iso}}]39 \item [{\np{ln_triad_iso}{ln\_triad\_iso}}] 45 40 See \autoref{sec:TRIADS_taper}. 46 41 If this is set false (the default), … … 53 48 giving an almost pure horizontal diffusive tracer flux within the mixed layer. 54 49 This is similar to the tapering suggested by \citet{gerdes.koberle.ea_CD91}. See \autoref{subsec:TRIADS_Gerdes-taper} 55 \item [{\np{ln_botmix_triad}{ln\_botmix\_triad}}]50 \item [{\np{ln_botmix_triad}{ln\_botmix\_triad}}] 56 51 See \autoref{sec:TRIADS_iso_bdry}. 57 52 If this is set false (the default) then the lateral diffusive fluxes … … 59 54 If it is set true, however, then these lateral diffusive fluxes are applied, 60 55 giving smoother bottom tracer fields at the cost of introducing diapycnal mixing. 61 \item [{\np{rn_sw_triad}{rn\_sw\_triad}}]56 \item [{\np{rn_sw_triad}{rn\_sw\_triad}}] 62 57 blah blah to be added.... 63 58 \end{description} 64 59 The options shared with the Standard scheme include: 65 60 \begin{description} 66 \item [{\np{ln_traldf_msc}{ln\_traldf\_msc}}] blah blah to be added67 \item [{\np{rn_slpmax}{rn\_slpmax}}] blah blah to be added61 \item [{\np{ln_traldf_msc}{ln\_traldf\_msc}}] blah blah to be added 62 \item [{\np{rn_slpmax}{rn\_slpmax}}] blah blah to be added 68 63 \end{description} 69 64 … … 548 543 The diffusion scheme satisfies the following six properties: 549 544 \begin{description} 550 \item [$\bullet$ horizontal diffusion]545 \item [$\bullet$ horizontal diffusion] 551 546 The discretization of the diffusion operator recovers the traditional five-point Laplacian 552 547 \autoref{eq:TRIADS_lat-normal} in the limit of flat iso-neutral direction: … … 559 554 \] 560 555 561 \item [$\bullet$ implicit treatment in the vertical]556 \item [$\bullet$ implicit treatment in the vertical] 562 557 Only tracer values associated with a single water column appear in the expression \autoref{eq:TRIADS_i33} for 563 558 the $_{33}$ fluxes, vertical fluxes driven by vertical gradients. … … 575 570 (where $b_w= e_{1w}\,e_{2w}\,e_{3w}$ is the volume of $w$-cells) can be quite large. 576 571 577 \item [$\bullet$ pure iso-neutral operator]572 \item [$\bullet$ pure iso-neutral operator] 578 573 The iso-neutral flux of locally referenced potential density is zero. 579 574 See \autoref{eq:TRIADS_latflux-rho} and \autoref{eq:TRIADS_vertflux-triad2}. 580 575 581 \item [$\bullet$ conservation of tracer]576 \item [$\bullet$ conservation of tracer] 582 577 The iso-neutral diffusion conserves tracer content, \ie 583 578 \[ … … 587 582 This property is trivially satisfied since the iso-neutral diffusive operator is written in flux form. 588 583 589 \item [$\bullet$ no increase of tracer variance]584 \item [$\bullet$ no increase of tracer variance] 590 585 The iso-neutral diffusion does not increase the tracer variance, \ie 591 586 \[ … … 600 595 the field on which it is applied becomes free of grid-point noise. 601 596 602 \item [$\bullet$ self-adjoint operator]597 \item [$\bullet$ self-adjoint operator] 603 598 The iso-neutral diffusion operator is self-adjoint, \ie 604 599 \begin{equation} … … 753 748 described above by \autoref{eq:TRIADS_Rtilde}. 754 749 \begin{enumerate} 755 \item 756 Mixed-layer depth is defined so as to avoid including regions of weak vertical stratification in 750 \item Mixed-layer depth is defined so as to avoid including regions of weak vertical stratification in 757 751 the slope definition. 758 752 At each $i,j$ (simplified to $i$ in \autoref{fig:TRIADS_MLB_triad}), … … 766 760 output the diagnosed mixed-layer depth $h_{\mathrm{ML}}=|z_{W}|_{k_{\mathrm{ML}}+1/2}$, 767 761 the depth of the $w$-point above the $i,k_{\mathrm{ML}}$ tracer point. 768 \item 769 We define `basal' triad slopes ${\:}_i{\mathbb{R}_{\mathrm{base}}}_{\,i_p}^{k_p}$ as 762 \item We define `basal' triad slopes ${\:}_i{\mathbb{R}_{\mathrm{base}}}_{\,i_p}^{k_p}$ as 770 763 the slopes of those triads whose vertical `arms' go down from the $i,k_{\mathrm{ML}}$ tracer point to 771 764 the $i,k_{\mathrm{ML}}-1$ tracer point below. … … 790 783 one gridbox deeper than the diagnosed ML depth $z_{\mathrm{ML}})$ that sets the $h$ used to taper the slopes in 791 784 \autoref{eq:TRIADS_rmtilde}. 792 \item 793 Finally, we calculate the adjusted triads ${\:}_i^k{\mathbb{R}_{\mathrm{ML}}}_{\,i_p}^{k_p}$ within 785 \item Finally, we calculate the adjusted triads ${\:}_i^k{\mathbb{R}_{\mathrm{ML}}}_{\,i_p}^{k_p}$ within 794 786 the mixed layer, by multiplying the appropriate ${\:}_i{\mathbb{R}_{\mathrm{base}}}_{\,i_p}^{k_p}$ by 795 787 the ratio of the depth of the $w$-point ${z_w}_{k+k_p}$ to ${z_{\mathrm{base}}}_{\,i}$. … … 872 864 % This may give strange looking results, 873 865 % particularly where the mixed-layer depth varies strongly laterally. 874 % ================================================================875 % Skew flux formulation for Eddy Induced Velocity :876 % ================================================================877 866 \section{Eddy induced advection formulated as a skew flux} 878 867 \label{sec:TRIADS_skew-flux}
Note: See TracChangeset
for help on using the changeset viewer.