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 13463 for NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc/latex/NEMO/subfiles/chap_STO.tex – NEMO

Ignore:
Timestamp:
2020-09-14T17:40:34+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2195:update to trunk 13461

Location:
NEMO/branches/2019/dev_r11351_fldread_with_XIOS
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8 
         9# SETTE 
         10^/utils/CI/sette@13382        sette 
  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc

    • Property svn:externals set to
      ^/utils/badges badges
      ^/utils/logos logos
  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc/latex

    • Property svn:ignore deleted
  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc/latex/NEMO

    • Property svn:externals set to
      ^/utils/figures/NEMO figures
  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/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_r11351_fldread_with_XIOS/doc/latex/NEMO/subfiles/chap_STO.tex

    r11344 r13463  
    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 
    10 \minitoc 
     8\thispagestyle{plain} 
     9 
     10\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. 
     
    118126either first principles, model simulations, or real-world observations. 
    119127The parameters are set by default as described in \cite{brankart_OM13}, which has been shown in the paper 
    120 to give good results for a global low resolution (2°) NEMO configuration. where this parametrization produces a major effect on the average large-scale circulation, especilally in regions of intense mesoscale activity. 
     128to give good results for a global low resolution (2°) \NEMO\ configuration. where this parametrization produces a major effect on the average large-scale circulation, especilally in regions of intense mesoscale activity. 
    121129The set of parameters will need further investigation to find appropriate values 
    122130for any other configuration or resolution of the model. 
    123131 
     132%% ================================================================================================= 
    124133\section{Implementation details} 
    125134\label{sec:STO_thech_details} 
    126135 
    127  
    128136The code implementing stochastic parametrisation is located in the src/OCE/STO directory. 
    129 It contains three modules :  
     137It contains three modules : 
    130138% \begin{description} 
    131139 
     
    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$). 
     
    150158the values $a^{(i)}, b^{(i)}, c^{(i)}$ for each autoregressive process, 
    151159as a function of the statistical properties required by the model user 
    152 (mean, standard deviation, time correlation, order of the process,\ldots).  
     160(mean, standard deviation, time correlation, order of the process,\ldots). 
    153161This routine also includes the initialization (seeding) of the random number generator. 
    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.}, 
    164 \ie when the state of the random number generator is read in the restart file.\\ 
    165  
    166 The implementation includes the basics for a few possible stochastic parametrisations including equation of state, lateral diffusion, horizontal pressure gradient, ice strength, trend, tracers dynamics. As for this release, only the stochastic parametrisation of equation of state is fully available and tested. \\  
     171when \np{ln_rstseed}{ln\_rstseed} is set to \forcode{.true.}, 
     172\ie\ when the state of the random number generator is read in the restart file.\\ 
     173 
     174The implementation includes the basics for a few possible stochastic parametrisations including equation of state, 
     175lateral diffusion, horizontal pressure gradient, ice strength, trend, tracers dynamics. 
     176As for this release, only the stochastic parametrisation of equation of state is fully available and tested. \\ 
    167177 
    168178Options and parameters \\ 
    169179 
    170 The \np{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. By default it set to \forcode{.false.}) and not active. 
    171 The set of parameters is available in \ngn{namsto} namelist(only the subset for equation of state stochastic parametrisation is listed below): 
    172 %---------------------------------------namsto-------------------------------------------------- 
    173  
    174 \nlst{namsto} 
    175 %-------------------------------------------------------------------------------------------------------------- 
     180The \np{ln_sto_eos}{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. 
     181By default it set to \forcode{.false.}) and not active. 
     182The set of parameters is available in \nam{sto}{sto} namelist 
     183(only the subset for equation of state stochastic parametrisation is listed below): 
     184 
     185\begin{listing} 
     186  \nlst{namsto} 
     187  \caption{\forcode{&namsto}} 
     188  \label{lst:namsto} 
     189\end{listing} 
    176190 
    177191The variables of stochastic paramtetrisation itself (based on the global 2° experiments as in \cite{brankart_OM13} are: 
     192 
    178193\begin{description} 
    179 \item[\np{nn\_sto\_eos}:]   number of independent random walks 
    180 \item[\np{rn\_eos\_stdxy}:] random walk horizontal standard deviation (in grid points) 
    181 \item[\np{rn\_eos\_stdz}:]  random walk vertical standard deviation (in grid points) 
    182 \item[\np{rn\_eos\_tcor}:]  random walk time correlation (in timesteps) 
    183 \item[\np{nn\_eos\_ord}:]   order of autoregressive processes 
    184 \item[\np{nn\_eos\_flt}:]   passes of Laplacian filter 
    185 \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) 
    186203\end{description} 
    187204 
    188205The first four parameters define the stochastic part of equation of state. 
    189 \biblio 
    190  
    191 \pindex 
     206 
     207\subinc{\input{../../global/epilogue}} 
    192208 
    193209\end{document} 
Note: See TracChangeset for help on using the changeset viewer.