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 11692 for NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles/chap_STO.tex – NEMO

Ignore:
Timestamp:
2019-10-12T16:08:18+02:00 (5 years ago)
Author:
francesca
Message:

Update branch to integrate the development starting from the current v4.01 ready trunk

Location:
NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc

    • Property svn:externals set to
      ^/utils/badges badges
      ^/utils/logos logos
  • NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex

    • Property svn:ignore deleted
  • NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO

    • Property svn:externals set to
      ^/utils/figures/NEMO figures
  • NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles

    • Property svn:ignore
      •  

        old new  
        1 *.aux 
        2 *.bbl 
        3 *.blg 
        4 *.dvi 
        5 *.fdb* 
        6 *.fls 
        7 *.idx 
         1*.ind 
        82*.ilg 
        9 *.ind 
        10 *.log 
        11 *.maf 
        12 *.mtc* 
        13 *.out 
        14 *.pdf 
        15 *.toc 
        16 _minted-* 
  • NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles/chap_STO.tex

    r11435 r11692  
    22 
    33\begin{document} 
    4 % ================================================================ 
    5 % Chapter stochastic parametrization of EOS (STO) 
    6 % ================================================================ 
     4 
    75\chapter{Stochastic Parametrization of EOS (STO)} 
    86\label{chap:STO} 
    97 
     8\thispagestyle{plain} 
     9 
    1010\chaptertoc 
     11 
     12\paragraph{Changes record} ~\\ 
     13 
     14{\footnotesize 
     15  \begin{tabularx}{\textwidth}{l||X|X} 
     16    Release & Author(s) & Modifications \\ 
     17    \hline 
     18    {\em   4.0} & {\em ...} & {\em ...} \\ 
     19    {\em   3.6} & {\em ...} & {\em ...} \\ 
     20    {\em   3.4} & {\em ...} & {\em ...} \\ 
     21    {\em <=3.4} & {\em ...} & {\em ...} 
     22  \end{tabularx} 
     23} 
    1124 
    1225% \vfill 
    1326% \begin{figure}[b] 
     27%% ================================================================================================= 
    1428% \subsubsection*{Changes record} 
    1529% \begin{tabular}{l||l|m{0.65\linewidth}} 
     
    2034% \end{figure} 
    2135 
    22 Authors: \\ 
    23 C. Levy release 4.0.1 update \\ 
    24 P.-A. Bouttier release 3.6 inital version 
    25  
    26 \newpage 
     36\clearpage 
    2737 
    2838As a result of the nonlinearity of the seawater equation of state, unresolved scales represent a major source of uncertainties in the computation of the large-scale horizontal density gradient from the large-scale temperature and salinity fields. Following  \cite{brankart_OM13}, the impact of these uncertainties can be simulated by random processes representing unresolved T/S fluctuations. The Stochastic Parametrization of EOS (STO) module implements this parametrization. 
     
    3040As detailed in \cite{brankart_OM13}, the stochastic formulation of the equation of state can be written as: 
    3141\begin{equation} 
    32   \label{eq:eos_sto} 
     42  \label{eq:STO_eos_sto} 
    3343  \rho = \frac{1}{2} \sum_{i=1}^m\{ \rho[T+\Delta T_i,S+\Delta S_i,p_o(z)] + \rho[T-\Delta T_i,S-\Delta S_i,p_o(z)] \} 
    3444\end{equation} 
     
    3747the scalar product of the respective local T/S gradients with random walks $\mathbf{\xi}$: 
    3848\begin{equation} 
    39   \label{eq:sto_pert} 
     49  \label{eq:STO_sto_pert} 
    4050  \Delta T_i = \mathbf{\xi}_i \cdot \nabla T \qquad \hbox{and} \qquad \Delta S_i = \mathbf{\xi}_i \cdot \nabla S 
    4151\end{equation} 
     
    4454$\mathbf{\xi}$ are uncorrelated over the horizontal and fully correlated along the vertical. 
    4555 
    46  
     56%% ================================================================================================= 
    4757\section{Stochastic processes} 
    4858\label{sec:STO_the_details} 
     
    5969 
    6070\begin{equation} 
    61   \label{eq:autoreg} 
     71  \label{eq:STO_autoreg} 
    6272  \xi^{(i)}_{k+1} = a^{(i)} \xi^{(i)}_k + b^{(i)} w^{(i)} + c^{(i)} 
    6373\end{equation} 
     
    6979 
    7080\begin{itemize} 
    71 \item 
    72   for order~1 processes, $w^{(i)}$ is a Gaussian white noise, with zero mean and standard deviation equal to~1, 
     81\item for order~1 processes, $w^{(i)}$ is a Gaussian white noise, with zero mean and standard deviation equal to~1, 
    7382  and the parameters $a^{(i)}$, $b^{(i)}$, $c^{(i)}$ are given by: 
    7483 
    7584  \[ 
    76     % \label{eq:ord1} 
     85    % \label{eq:STO_ord1} 
    7786    \left\{ 
    7887      \begin{array}{l} 
     
    8493  \] 
    8594 
    86 \item 
    87   for order~$n>1$ processes, $w^{(i)}$ is an order~$n-1$ autoregressive process, with zero mean, 
     95\item for order~$n>1$ processes, $w^{(i)}$ is an order~$n-1$ autoregressive process, with zero mean, 
    8896  standard deviation equal to~$\sigma^{(i)}$; 
    8997  correlation timescale equal to~$\tau^{(i)}$; 
     
    9199 
    92100  \begin{equation} 
    93     \label{eq:ord2} 
     101    \label{eq:STO_ord2} 
    94102    \left\{ 
    95103      \begin{array}{l} 
     
    107115\noindent 
    108116In this way, higher order processes can be easily generated recursively using the same piece of code implementing 
    109 \autoref{eq:autoreg}, and using successive processes from order $0$ to~$n-1$ as~$w^{(i)}$. 
    110 The parameters in \autoref{eq:ord2} are computed so that this recursive application of 
    111 \autoref{eq:autoreg} leads to processes with the required standard deviation and correlation timescale, 
     117\autoref{eq:STO_autoreg}, and using successive processes from order $0$ to~$n-1$ as~$w^{(i)}$. 
     118The parameters in \autoref{eq:STO_ord2} are computed so that this recursive application of 
     119\autoref{eq:STO_autoreg} leads to processes with the required standard deviation and correlation timescale, 
    112120with the additional condition that the $n-1$ first derivatives of the autocorrelation function are equal to 
    113121zero at~$t=0$, so that the resulting processes become smoother and smoother as $n$ increases. 
     
    122130for any other configuration or resolution of the model. 
    123131 
     132%% ================================================================================================= 
    124133\section{Implementation details} 
    125134\label{sec:STO_thech_details} 
    126  
    127135 
    128136The code implementing stochastic parametrisation is located in the src/OCE/STO directory. 
     
    135143 
    136144\mdl{stopts} : stochastic parametrisation associated with the non-linearity of the equation of 
    137  seawater, implementing \autoref{eq:sto_pert} so as specifics in the equation of state 
    138  implementing \autoref{eq:eos_sto}. 
     145 seawater, implementing \autoref{eq:STO_sto_pert} so as specifics in the equation of state 
     146 implementing \autoref{eq:STO_eos_sto}. 
    139147% \end{description} 
    140148 
    141149The \mdl{stopar} module includes three public routines called in the model: 
    142150 
    143 (\rou{sto\_par}) is a direct implementation of \autoref{eq:autoreg}, 
     151(\rou{sto\_par}) is a direct implementation of \autoref{eq:STO_autoreg}, 
    144152applied at each model grid point (in 2D or 3D), and called at each model time step ($k$) to 
    145153update every autoregressive process ($i=1,\ldots,m$). 
     
    154162 
    155163(\rou{sto\_rst\_write}) writes a restart file 
    156 (which suffix name is given by \np{cn\_storst\_out} namelist parameter) containing the current value of 
     164(which suffix name is given by \np{cn_storst_out}{cn\_storst\_out} namelist parameter) containing the current value of 
    157165all autoregressive processes to allow creating the file needed for a restart. 
    158166This restart file also contains the current state of the random number generator. 
    159 When \np{ln\_rststo} is set to \forcode{.true.}), 
    160 the restart file (which suffix name is given by \np{cn\_storst\_in} namelist parameter) is read by 
     167When \np{ln_rststo}{ln\_rststo} is set to \forcode{.true.}), 
     168the restart file (which suffix name is given by \np{cn_storst_in}{cn\_storst\_in} namelist parameter) is read by 
    161169the initialization routine (\rou{sto\_par\_init}). 
    162170The simulation will continue exactly as if it was not interrupted only 
    163 when \np{ln\_rstseed} is set to \forcode{.true.}, 
     171when \np{ln_rstseed}{ln\_rstseed} is set to \forcode{.true.}, 
    164172\ie\ when the state of the random number generator is read in the restart file.\\ 
    165173 
     
    170178Options and parameters \\ 
    171179 
    172 The \np{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. 
     180The \np{ln_sto_eos}{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. 
    173181By default it set to \forcode{.false.}) and not active. 
    174 The set of parameters is available in \nam{sto} namelist 
     182The set of parameters is available in \nam{sto}{sto} namelist 
    175183(only the subset for equation of state stochastic parametrisation is listed below): 
    176 %---------------------------------------namsto-------------------------------------------------- 
    177  
    178 \nlst{namsto} 
    179 %-------------------------------------------------------------------------------------------------------------- 
     184 
     185\begin{listing} 
     186  \nlst{namsto} 
     187  \caption{\forcode{&namsto}} 
     188  \label{lst:namsto} 
     189\end{listing} 
    180190 
    181191The variables of stochastic paramtetrisation itself (based on the global 2° experiments as in \cite{brankart_OM13} are: 
    182192 
    183193\begin{description} 
    184 \item[\np{nn\_sto\_eos}:]   number of independent random walks 
    185 \item[\np{rn\_eos\_stdxy}:] random walk horizontal standard deviation (in grid points) 
    186 \item[\np{rn\_eos\_stdz}:]  random walk vertical standard deviation (in grid points) 
    187 \item[\np{rn\_eos\_tcor}:]  random walk time correlation (in timesteps) 
    188 \item[\np{nn\_eos\_ord}:]   order of autoregressive processes 
    189 \item[\np{nn\_eos\_flt}:]   passes of Laplacian filter 
    190 \item[\np{rn\_eos\_lim}:]   limitation factor (default = 3.0) 
     194\item [{\np{nn_sto_eos}{nn\_sto\_eos}:}]     number of independent random walks 
     195\item [{\np{rn_eos_stdxy}{rn\_eos\_stdxy}:}] random walk horizontal standard deviation 
     196  (in grid points) 
     197\item [{\np{rn_eos_stdz}{rn\_eos\_stdz}:}]   random walk vertical standard deviation 
     198  (in grid points) 
     199\item [{\np{rn_eos_tcor}{rn\_eos\_tcor}:}]   random walk time correlation (in timesteps) 
     200\item [{\np{nn_eos_ord}{nn\_eos\_ord}:}]     order of autoregressive processes 
     201\item [{\np{nn_eos_flt}{nn\_eos\_flt}:}]     passes of Laplacian filter 
     202\item [{\np{rn_eos_lim}{rn\_eos\_lim}:}]     limitation factor (default = 3.0) 
    191203\end{description} 
    192204 
    193205The first four parameters define the stochastic part of equation of state. 
    194 \biblio 
    195  
    196 \pindex 
     206 
     207\onlyinsubfile{\input{../../global/epilogue}} 
    197208 
    198209\end{document} 
Note: See TracChangeset for help on using the changeset viewer.