Changeset 13463 for NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc/latex/NEMO/subfiles/chap_STO.tex
- Timestamp:
- 2020-09-14T17:40:34+02:00 (4 years ago)
- 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 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@13382 sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc
-
Property
svn:externals
set to
^/utils/badges badges
^/utils/logos logos
-
Property
svn:externals
set to
-
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
-
Property
svn:externals
set to
-
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 8 2 *.ilg 9 *.ind10 *.log11 *.maf12 *.mtc*13 *.out14 *.pdf15 *.toc16 _minted-*
-
- Property svn:ignore
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/doc/latex/NEMO/subfiles/chap_STO.tex
r11344 r13463 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 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 } 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. … … 118 126 either first principles, model simulations, or real-world observations. 119 127 The 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°) NEMOconfiguration. where this parametrization produces a major effect on the average large-scale circulation, especilally in regions of intense mesoscale activity.128 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. 121 129 The set of parameters will need further investigation to find appropriate values 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 135 127 128 136 The code implementing stochastic parametrisation is located in the src/OCE/STO directory. 129 It contains three modules : 137 It contains three modules : 130 138 % \begin{description} 131 139 … … 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$). … … 150 158 the values $a^{(i)}, b^{(i)}, c^{(i)}$ for each autoregressive process, 151 159 as 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). 153 161 This routine also includes the initialization (seeding) of the random number generator. 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.}, 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. \\ 171 when \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 174 The implementation includes the basics for a few possible stochastic parametrisations including equation of state, 175 lateral diffusion, horizontal pressure gradient, ice strength, trend, tracers dynamics. 176 As for this release, only the stochastic parametrisation of equation of state is fully available and tested. \\ 167 177 168 178 Options and parameters \\ 169 179 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 %-------------------------------------------------------------------------------------------------------------- 180 The \np{ln_sto_eos}{ln\_sto\_eos} namelist variable activates stochastic parametrisation of equation of state. 181 By default it set to \forcode{.false.}) and not active. 182 The 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} 176 190 177 191 The variables of stochastic paramtetrisation itself (based on the global 2° experiments as in \cite{brankart_OM13} are: 192 178 193 \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) 186 203 \end{description} 187 204 188 205 The first four parameters define the stochastic part of equation of state. 189 \biblio 190 191 \pindex 206 207 \subinc{\input{../../global/epilogue}} 192 208 193 209 \end{document}
Note: See TracChangeset
for help on using the changeset viewer.