Changeset 11692 for NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles/chap_STO.tex
- Timestamp:
- 2019-10-12T16:08:18+02:00 (5 years ago)
- 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
-
Property
svn:externals
set to
-
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
-
Property
svn:externals
set to
-
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 8 2 *.ilg 9 *.ind10 *.log11 *.maf12 *.mtc*13 *.out14 *.pdf15 *.toc16 _minted-*
-
- Property svn:ignore
-
NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles/chap_STO.tex
r11435 r11692 2 2 3 3 \begin{document} 4 % ================================================================ 5 % Chapter stochastic parametrization of EOS (STO) 6 % ================================================================ 4 7 5 \chapter{Stochastic Parametrization of EOS (STO)} 8 6 \label{chap:STO} 9 7 8 \thispagestyle{plain} 9 10 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 } 11 24 12 25 % \vfill 13 26 % \begin{figure}[b] 27 %% ================================================================================================= 14 28 % \subsubsection*{Changes record} 15 29 % \begin{tabular}{l||l|m{0.65\linewidth}} … … 20 34 % \end{figure} 21 35 22 Authors: \\ 23 C. Levy release 4.0.1 update \\ 24 P.-A. Bouttier release 3.6 inital version 25 26 \newpage 36 \clearpage 27 37 28 38 As 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. … … 30 40 As detailed in \cite{brankart_OM13}, the stochastic formulation of the equation of state can be written as: 31 41 \begin{equation} 32 \label{eq: eos_sto}42 \label{eq:STO_eos_sto} 33 43 \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)] \} 34 44 \end{equation} … … 37 47 the scalar product of the respective local T/S gradients with random walks $\mathbf{\xi}$: 38 48 \begin{equation} 39 \label{eq: sto_pert}49 \label{eq:STO_sto_pert} 40 50 \Delta T_i = \mathbf{\xi}_i \cdot \nabla T \qquad \hbox{and} \qquad \Delta S_i = \mathbf{\xi}_i \cdot \nabla S 41 51 \end{equation} … … 44 54 $\mathbf{\xi}$ are uncorrelated over the horizontal and fully correlated along the vertical. 45 55 46 56 %% ================================================================================================= 47 57 \section{Stochastic processes} 48 58 \label{sec:STO_the_details} … … 59 69 60 70 \begin{equation} 61 \label{eq: autoreg}71 \label{eq:STO_autoreg} 62 72 \xi^{(i)}_{k+1} = a^{(i)} \xi^{(i)}_k + b^{(i)} w^{(i)} + c^{(i)} 63 73 \end{equation} … … 69 79 70 80 \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, 73 82 and the parameters $a^{(i)}$, $b^{(i)}$, $c^{(i)}$ are given by: 74 83 75 84 \[ 76 % \label{eq: ord1}85 % \label{eq:STO_ord1} 77 86 \left\{ 78 87 \begin{array}{l} … … 84 93 \] 85 94 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, 88 96 standard deviation equal to~$\sigma^{(i)}$; 89 97 correlation timescale equal to~$\tau^{(i)}$; … … 91 99 92 100 \begin{equation} 93 \label{eq: ord2}101 \label{eq:STO_ord2} 94 102 \left\{ 95 103 \begin{array}{l} … … 107 115 \noindent 108 116 In 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 of111 \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)}$. 118 The 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, 112 120 with the additional condition that the $n-1$ first derivatives of the autocorrelation function are equal to 113 121 zero at~$t=0$, so that the resulting processes become smoother and smoother as $n$ increases. … … 122 130 for any other configuration or resolution of the model. 123 131 132 %% ================================================================================================= 124 133 \section{Implementation details} 125 134 \label{sec:STO_thech_details} 126 127 135 128 136 The code implementing stochastic parametrisation is located in the src/OCE/STO directory. … … 135 143 136 144 \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 state138 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}. 139 147 % \end{description} 140 148 141 149 The \mdl{stopar} module includes three public routines called in the model: 142 150 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}, 144 152 applied at each model grid point (in 2D or 3D), and called at each model time step ($k$) to 145 153 update every autoregressive process ($i=1,\ldots,m$). … … 154 162 155 163 (\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 of164 (which suffix name is given by \np{cn_storst_out}{cn\_storst\_out} namelist parameter) containing the current value of 157 165 all autoregressive processes to allow creating the file needed for a restart. 158 166 This 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 by167 When \np{ln_rststo}{ln\_rststo} is set to \forcode{.true.}), 168 the restart file (which suffix name is given by \np{cn_storst_in}{cn\_storst\_in} namelist parameter) is read by 161 169 the initialization routine (\rou{sto\_par\_init}). 162 170 The simulation will continue exactly as if it was not interrupted only 163 when \np{ln \_rstseed} is set to \forcode{.true.},171 when \np{ln_rstseed}{ln\_rstseed} is set to \forcode{.true.}, 164 172 \ie\ when the state of the random number generator is read in the restart file.\\ 165 173 … … 170 178 Options and parameters \\ 171 179 172 The \np{ln \_sto\_eos} namelist variable activates stochastic parametrisation of equation of state.180 The \np{ln_sto_eos}{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. 173 181 By default it set to \forcode{.false.}) and not active. 174 The set of parameters is available in \nam{sto} namelist182 The set of parameters is available in \nam{sto}{sto} namelist 175 183 (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} 180 190 181 191 The variables of stochastic paramtetrisation itself (based on the global 2° experiments as in \cite{brankart_OM13} are: 182 192 183 193 \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) 191 203 \end{description} 192 204 193 205 The first four parameters define the stochastic part of equation of state. 194 \biblio 195 196 \pindex 206 207 \onlyinsubfile{\input{../../global/epilogue}} 197 208 198 209 \end{document}
Note: See TracChangeset
for help on using the changeset viewer.