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 11123 for NEMO/trunk/doc/latex/NEMO/subfiles/annex_iso.tex – NEMO

Ignore:
Timestamp:
2019-06-17T14:22:27+02:00 (5 years ago)
Author:
nicolasmartin
Message:

Modification of LaTeX subfiles accordingly to new citations keys

Location:
NEMO/trunk/doc/latex/NEMO/subfiles
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/doc/latex/NEMO/subfiles

    • Property svn:ignore deleted
  • NEMO/trunk/doc/latex/NEMO/subfiles/annex_iso.tex

    r10442 r11123  
    5252  the vertical skew flux is further reduced to ensure no vertical buoyancy flux, 
    5353  giving an almost pure horizontal diffusive tracer flux within the mixed layer. 
    54   This is similar to the tapering suggested by \citet{Gerdes1991}. See \autoref{subsec:Gerdes-taper} 
     54  This is similar to the tapering suggested by \citet{gerdes.koberle.ea_CD91}. See \autoref{subsec:Gerdes-taper} 
    5555\item[\np{ln\_botmix\_triad}] 
    5656  See \autoref{sec:iso_bdry}.  
     
    7171\label{sec:iso} 
    7272 
    73 We have implemented into \NEMO a scheme inspired by \citet{Griffies_al_JPO98}, 
     73We have implemented into \NEMO a scheme inspired by \citet{griffies.gnanadesikan.ea_JPO98}, 
    7474but formulated within the \NEMO framework, using scale factors rather than grid-sizes. 
    7575 
     
    194194\subsection{Expression of the skew-flux in terms of triad slopes} 
    195195 
    196 \citep{Griffies_al_JPO98} introduce a different discretization of the off-diagonal terms that 
     196\citep{griffies.gnanadesikan.ea_JPO98} introduce a different discretization of the off-diagonal terms that 
    197197nicely solves the problem. 
    198198% Instead of multiplying the mean slope calculated at the $u$-point by 
     
    473473 
    474474To complete the discretization we now need only specify the triad volumes $_i^k\mathbb{V}_{i_p}^{k_p}$. 
    475 \citet{Griffies_al_JPO98} identifies these $_i^k\mathbb{V}_{i_p}^{k_p}$ as the volumes of the quarter cells, 
     475\citet{griffies.gnanadesikan.ea_JPO98} identifies these $_i^k\mathbb{V}_{i_p}^{k_p}$ as the volumes of the quarter cells, 
    476476defined in terms of the distances between $T$, $u$,$f$ and $w$-points. 
    477477This is the natural discretization of \autoref{eq:cts-var}. 
     
    685685As discussed in \autoref{subsec:LDF_slp_iso}, 
    686686iso-neutral slopes relative to geopotentials must be bounded everywhere, 
    687 both for consistency with the small-slope approximation and for numerical stability \citep{Cox1987, Griffies_Bk04}. 
     687both for consistency with the small-slope approximation and for numerical stability \citep{cox_OM87, griffies_bk04}. 
    688688The bound chosen in \NEMO is applied to each component of the slope separately and 
    689689has a value of $1/100$ in the ocean interior. 
     
    859859\footnote{ 
    860860  To ensure good behaviour where horizontal density gradients are weak, 
    861   we in fact follow \citet{Gerdes1991} and 
     861  we in fact follow \citet{gerdes.koberle.ea_CD91} and 
    862862  set $\rML^*=\mathrm{sgn}(\tilde{r}_i)\min(|\rMLt^2/\tilde{r}_i|,|\tilde{r}_i|)-\sigma_i$. 
    863863} 
     
    865865This approach is similar to multiplying the iso-neutral diffusion coefficient by 
    866866$\tilde{r}_{\mathrm{max}}^{-2}\tilde{r}_i^{-2}$ for steep slopes, 
    867 as suggested by \citet{Gerdes1991} (see also \citet{Griffies_Bk04}). 
     867as suggested by \citet{gerdes.koberle.ea_CD91} (see also \citet{griffies_bk04}). 
    868868Again it is applied separately to each triad $_i^k\mathbb{R}_{i_p}^{k_p}$ 
    869869 
     
    925925 
    926926However, when \np{ln\_traldf\_triad} is set true, 
    927 \NEMO instead implements eddy induced advection according to the so-called skew form \citep{Griffies_JPO98}. 
     927\NEMO instead implements eddy induced advection according to the so-called skew form \citep{griffies_JPO98}. 
    928928It is based on a transformation of the advective fluxes using the non-divergent nature of the eddy induced velocity. 
    929929For example in the (\textbf{i},\textbf{k}) plane, 
     
    11391139it is equivalent to a horizontal eiv (eddy-induced velocity) that is uniform within the mixed layer 
    11401140\autoref{eq:eiv_v}. 
    1141 This ensures that the eiv velocities do not restratify the mixed layer \citep{Treguier1997,Danabasoglu_al_2008}. 
     1141This ensures that the eiv velocities do not restratify the mixed layer \citep{treguier.held.ea_JPO97,danabasoglu.ferrari.ea_JC08}. 
    11421142Equivantly, in terms of the skew-flux formulation we use here, 
    11431143the linear slope tapering within the mixed-layer gives a linearly varying vertical flux, 
     
    11531153$uw$ (integer +1/2 $i$, integer $j$, integer +1/2 $k$) and $vw$ (integer $i$, integer +1/2 $j$, integer +1/2 $k$) 
    11541154points (see Table \autoref{tab:cell}) respectively. 
    1155 We follow \citep{Griffies_Bk04} and calculate the streamfunction at a given $uw$-point from 
     1155We follow \citep{griffies_bk04} and calculate the streamfunction at a given $uw$-point from 
    11561156the surrounding four triads according to: 
    11571157\[ 
Note: See TracChangeset for help on using the changeset viewer.