Changeset 11692 for NEMO/branches/2019/dev_r11514_HPC-02_single-core-extrahalo/doc/latex/NEMO/subfiles/chap_model_basics_zstar.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_model_basics_zstar.tex
r11435 r11692 2 2 3 3 \begin{document} 4 % ================================================================ 5 % Chapter 1 Model Basics 6 % ================================================================ 7 % ================================================================ 8 % Curvilinear \zstar- \sstar-coordinate System 9 % ================================================================ 4 10 5 \chapter{ essai \zstar \sstar} 6 7 \thispagestyle{plain} 8 9 \chaptertoc 10 11 \paragraph{Changes record} ~\\ 12 13 {\footnotesize 14 \begin{tabularx}{\textwidth}{l||X|X} 15 Release & Author(s) & Modifications \\ 16 \hline 17 {\em 4.0} & {\em ...} & {\em ...} \\ 18 {\em 3.6} & {\em ...} & {\em ...} \\ 19 {\em 3.4} & {\em ...} & {\em ...} \\ 20 {\em <=3.4} & {\em ...} & {\em ...} 21 \end{tabularx} 22 } 23 24 \clearpage 25 26 %% ================================================================================================= 11 27 \section{Curvilinear \zstar- or \sstar coordinate system} 12 13 % -------------------------------------------------------------------------------------------------------------14 % ????15 % -------------------------------------------------------------------------------------------------------------16 28 17 29 \colorbox{yellow}{ to be updated } … … 26 38 To overcome problems with vanishing surface and/or bottom cells, we consider the zstar coordinate 27 39 \[ 28 % \label{eq: PE_}40 % \label{eq:MBZ_PE_} 29 41 z^\star = H \left( \frac{z-\eta}{H+\eta} \right) 30 42 \] … … 40 52 the surface height, it is clear that surfaces constant $z^\star$ are very similar to the depth surfaces. 41 53 These properties greatly reduce difficulties of computing the horizontal pressure gradient relative to 42 terrain following sigma models discussed in \autoref{subsec: PE_sco}.54 terrain following sigma models discussed in \autoref{subsec:MB_sco}. 43 55 Additionally, since $z^\star$ when $\eta = 0$, no flow is spontaneously generated in 44 56 an unforced ocean starting from rest, regardless the bottom topography. … … 66 78 the surface height, again so long as $\eta > -H$. 67 79 68 %%%69 80 % essai update time splitting... 70 %%% 71 72 % ================================================================ 73 % Surface Pressure Gradient and Sea Surface Height 74 % ================================================================ 75 \section[Surface pressure gradient and sea surface heigth (\textit{dynspg.F90})] 76 {Surface pressure gradient and sea surface heigth (\protect\mdl{dynspg})} 77 \label{sec:DYN_hpg_spg} 78 %-----------------------------------------nam_dynspg---------------------------------------------------- 81 82 %% ================================================================================================= 83 \section[Surface pressure gradient and sea surface heigth (\textit{dynspg.F90})]{Surface pressure gradient and sea surface heigth (\protect\mdl{dynspg})} 84 \label{sec:MBZ_dyn_hpg_spg} 79 85 80 86 %\nlst{nam_dynspg} 81 %------------------------------------------------------------------------------------------------------------ 82 Options are defined through the \nam{\_dynspg} namelist variables. 83 The surface pressure gradient term is related to the representation of the free surface (\autoref{sec:PE_hor_pg}). 87 Options are defined through the \nam{_dynspg}{\_dynspg} namelist variables. 88 The surface pressure gradient term is related to the representation of the free surface (\autoref{sec:MB_hor_pg}). 84 89 The main distinction is between the fixed volume case (linear free surface or rigid lid) and 85 90 the variable volume case (nonlinear free surface, \key{vvl} is active). 86 In the linear free surface case (\autoref{subsec: PE_free_surface}) and rigid lid (\autoref{PE_rigid_lid}),91 In the linear free surface case (\autoref{subsec:MB_free_surface}) and rigid lid (\autoref{PE_rigid_lid}), 87 92 the vertical scale factors $e_{3}$ are fixed in time, 88 while in the nonlinear case (\autoref{subsec: PE_free_surface}) they are time-dependent.93 while in the nonlinear case (\autoref{subsec:MB_free_surface}) they are time-dependent. 89 94 With both linear and nonlinear free surface, external gravity waves are allowed in the equations, 90 95 which imposes a very small time step when an explicit time stepping is used. 91 96 Two methods are proposed to allow a longer time step for the three-dimensional equations: 92 the filtered free surface, which is a modification of the continuous equations %(see \autoref{eq: PE_flt?}),97 the filtered free surface, which is a modification of the continuous equations %(see \autoref{eq:MB_flt?}), 93 98 and the split-explicit free surface described below. 94 99 The extra term introduced in the filtered method is calculated implicitly, 95 100 so that the update of the next velocities is done in module \mdl{dynspg\_flt} and not in \mdl{dynnxt}. 96 101 97 %-------------------------------------------------------------98 102 % Explicit 99 %------------------------------------------------------------- 100 \subsubsection[Explicit (\texttt{\textbf{key\_dynspg\_exp}})] 101 {Explicit (\protect\key{dynspg\_exp})} 102 \label{subsec:DYN_spg_exp} 103 %% ================================================================================================= 104 \subsubsection[Explicit (\texttt{\textbf{key\_dynspg\_exp}})]{Explicit (\protect\key{dynspg\_exp})} 105 \label{subsec:MBZ_dyn_spg_exp} 103 106 104 107 In the explicit free surface formulation, the model time step is chosen small enough to … … 106 109 The sea surface height is given by: 107 110 \begin{equation} 108 \label{eq: dynspg_ssh}111 \label{eq:MBZ_dynspg_ssh} 109 112 \frac{\partial \eta }{\partial t}\equiv \frac{\text{EMP}}{\rho_w }+\frac{1}{e_{1T} 110 113 e_{2T} }\sum\limits_k {\left( {\delta_i \left[ {e_{2u} e_{3u} u} … … 116 119 and $\rho_w =1,000\,Kg.m^{-3}$ is the volumic mass of pure water. 117 120 The sea-surface height is evaluated using a leapfrog scheme in combination with an Asselin time filter, 118 (\ie\ the velocity appearing in (\autoref{eq: dynspg_ssh}) is centred in time (\textit{now} velocity).121 (\ie\ the velocity appearing in (\autoref{eq:DYN_spg_ssh}) is centred in time (\textit{now} velocity). 119 122 120 123 The surface pressure gradient, also evaluated using a leap-frog scheme, is then simply given by: 121 124 \begin{equation} 122 \label{eq: dynspg_exp}125 \label{eq:MBZ_dynspg_exp} 123 126 \left\{ 124 127 \begin{aligned} … … 130 133 131 134 Consistent with the linearization, a $\left. \rho \right|_{k=1} / \rho_o$ factor is omitted in 132 (\autoref{eq:dynspg_exp}). 133 134 %------------------------------------------------------------- 135 (\autoref{eq:DYN_spg_exp}). 136 135 137 % Split-explicit time-stepping 136 %------------------------------------------------------------- 137 \subsubsection[Split-explicit time-stepping (\texttt{\textbf{key\_dynspg\_ts}})] 138 {Split-explicit time-stepping (\protect\key{dynspg\_ts})} 139 \label{subsec:DYN_spg_ts} 140 %--------------------------------------------namdom---------------------------------------------------- 141 142 \nlst{namdom} 143 %-------------------------------------------------------------------------------------------------------------- 138 %% ================================================================================================= 139 \subsubsection[Split-explicit time-stepping (\texttt{\textbf{key\_dynspg\_ts}})]{Split-explicit time-stepping (\protect\key{dynspg\_ts})} 140 \label{subsec:MBZ_dyn_spg_ts} 141 144 142 The split-explicit free surface formulation used in OPA follows the one proposed by \citet{Griffies2004?}. 145 143 The general idea is to solve the free surface equation with a small time step, 146 144 while the three dimensional prognostic variables are solved with a longer time step that 147 is a multiple of \np{rdtbt} in the \nam{dom} namelist (Figure III.3). 148 149 %> > > > > > > > > > > > > > > > > > > > > > > > > > > > 145 is a multiple of \np{rdtbt}{rdtbt} in the \nam{dom}{dom} namelist (Figure III.3). 146 150 147 \begin{figure}[!t] 151 \ begin{center}152 \includegraphics[width=\textwidth]{Fig_DYN_dynspg_ts}153 \caption{154 \protect\label{fig:DYN_dynspg_ts}155 156 157 158 159 160 161 162 163 164 165 166 167 168 A baroclinic leap-frog time step carries the surface height to $t+\Delta t$ using the convergence of169 the time averaged vertically integrated velocity taken from baroclinic time step t.170 }171 \ end{center}148 \centering 149 \includegraphics[width=0.66\textwidth]{MBZ_DYN_dynspg_ts} 150 \caption[Schematic of the split-explicit time stepping scheme for 151 the barotropic and baroclinic modes, after \citet{Griffies2004?}]{ 152 Schematic of the split-explicit time stepping scheme for the barotropic and baroclinic modes, 153 after \citet{Griffies2004?}. 154 Time increases to the right. 155 Baroclinic time steps are denoted by $t-\Delta t$, $t, t+\Delta t$, and $t+2\Delta t$. 156 The curved line represents a leap-frog time step, 157 and the smaller barotropic time steps $N \Delta t=2\Delta t$ are denoted by the zig-zag line. 158 The vertically integrated forcing \textbf{M}(t) computed at 159 baroclinic time step t represents the interaction between the barotropic and baroclinic motions. 160 While keeping the total depth, tracer, and freshwater forcing fields fixed, 161 a leap-frog integration carries the surface height and vertically integrated velocity from 162 t to $t+2 \Delta t$ using N barotropic time steps of length $\Delta t$. 163 Time averaging the barotropic fields over the N+1 time steps (endpoints included) 164 centers the vertically integrated velocity at the baroclinic timestep $t+\Delta t$. 165 A baroclinic leap-frog time step carries the surface height to $t+\Delta t$ using 166 the convergence of the time averaged vertically integrated velocity taken from 167 baroclinic time step t.} 168 \label{fig:MBZ_dyn_dynspg_ts} 172 169 \end{figure} 173 %> > > > > > > > > > > > > > > > > > > > > > > > > > > >174 170 175 171 The split-explicit formulation has a damping effect on external gravity waves, … … 186 182 We have 187 183 \[ 188 % \label{eq: DYN_spg_ts_eta}184 % \label{eq:MBZ_dyn_spg_ts_eta} 189 185 \eta^{(b)}(\tau,t_{n+1}) - \eta^{(b)}(\tau,t_{n+1}) (\tau,t_{n-1}) 190 186 = 2 \Delta t \left[-\nabla \cdot \textbf{U}^{(b)}(\tau,t_n) + \text{EMP}_w(\tau) \right] 191 187 \] 192 188 \begin{multline*} 193 % \label{eq: DYN_spg_ts_u}189 % \label{eq:MBZ_dyn_spg_ts_u} 194 190 \textbf{U}^{(b)}(\tau,t_{n+1}) - \textbf{U}^{(b)}(\tau,t_{n-1}) \\ 195 191 = 2\Delta t \left[ - f \textbf{k} \times \textbf{U}^{(b)}(\tau,t_{n}) … … 207 203 This is also the time that sets the barotropic time steps via 208 204 \[ 209 % \label{eq: DYN_spg_ts_t}205 % \label{eq:MBZ_dyn_spg_ts_t} 210 206 t_n=\tau+n\Delta t 211 207 \] … … 213 209 The density scaled surface pressure is evaluated via 214 210 \[ 215 % \label{eq: DYN_spg_ts_ps}211 % \label{eq:MBZ_dyn_spg_ts_ps} 216 212 p_s^{(b)}(\tau,t_{n}) = 217 213 \begin{cases} … … 222 218 To get started, we assume the following initial conditions 223 219 \[ 224 % \label{eq: DYN_spg_ts_eta}220 % \label{eq:MBZ_dyn_spg_ts_eta} 225 221 \begin{split} 226 222 \eta^{(b)}(\tau,t_{n=0}) &= \overline{\eta^{(b)}(\tau)} \\ … … 230 226 with 231 227 \[ 232 % \label{eq: DYN_spg_ts_etaF}228 % \label{eq:MBZ_dyn_spg_ts_etaF} 233 229 \overline{\eta^{(b)}(\tau)} = \frac{1}{N+1} \sum\limits_{n=0}^N \eta^{(b)}(\tau-\Delta t,t_{n}) 234 230 \] … … 236 232 Likewise, 237 233 \[ 238 % \label{eq: DYN_spg_ts_u}234 % \label{eq:MBZ_dyn_spg_ts_u} 239 235 \textbf{U}^{(b)}(\tau,t_{n=0}) = \overline{\textbf{U}^{(b)}(\tau)} \\ \\ 240 236 \textbf{U}(\tau,t_{n=1}) = \textbf{U}^{(b)}(\tau,t_{n=0}) + \Delta t \ \text{RHS}_{n=0} … … 242 238 with 243 239 \[ 244 % \label{eq: DYN_spg_ts_u}240 % \label{eq:MBZ_dyn_spg_ts_u} 245 241 \overline{\textbf{U}^{(b)}(\tau)} = \frac{1}{N+1} \sum\limits_{n=0}^N\textbf{U}^{(b)}(\tau-\Delta t,t_{n}) 246 242 \] … … 251 247 produce the updated vertically integrated velocity at baroclinic time $\tau + \Delta \tau$ 252 248 \[ 253 % \label{eq: DYN_spg_ts_u}249 % \label{eq:MBZ_dyn_spg_ts_u} 254 250 \textbf{U}(\tau+\Delta t) = \overline{\textbf{U}^{(b)}(\tau+\Delta t)} 255 251 = \frac{1}{N+1} \sum\limits_{n=0}^N\textbf{U}^{(b)}(\tau,t_{n}) … … 258 254 a baroclinic leap-frog using the following form 259 255 \begin{equation} 260 \label{eq: DYN_spg_ts_ssh}256 \label{eq:MBZ_dyn_spg_ts_ssh} 261 257 \eta(\tau+\Delta) - \eta^{F}(\tau-\Delta) = 2\Delta t \ \left[ - \nabla \cdot \textbf{U}(\tau) + \text{EMP}_w \right] 262 258 \end{equation} … … 267 263 268 264 In general, some form of time filter is needed to maintain integrity of the surface height field due to 269 the leap-frog splitting mode in equation \autoref{eq: DYN_spg_ts_ssh}.265 the leap-frog splitting mode in equation \autoref{eq:MBZ_dyn_spg_ts_ssh}. 270 266 We have tried various forms of such filtering, 271 267 with the following method discussed in Griffies et al. (2001) chosen due to its stability and … … 273 269 274 270 \begin{equation} 275 \label{eq: DYN_spg_ts_sshf}271 \label{eq:MBZ_dyn_spg_ts_sshf} 276 272 \eta^{F}(\tau-\Delta) = \overline{\eta^{(b)}(\tau)} 277 273 \end{equation} … … 279 275 280 276 \[ 281 % \label{eq: DYN_spg_ts_sshf2}277 % \label{eq:MBZ_dyn_spg_ts_sshf2} 282 278 \eta^{F}(\tau-\Delta) = \eta(\tau) 283 279 + (\alpha/2) \left[\overline{\eta^{(b)}}(\tau+\Delta t) … … 288 284 This isolation allows for an easy check that tracer conservation is exact when eliminating tracer and 289 285 surface height time filtering (see ?? for more complete discussion). 290 However, in the general case with a non-zero $\alpha$, the filter \autoref{eq: DYN_spg_ts_sshf} was found to286 However, in the general case with a non-zero $\alpha$, the filter \autoref{eq:MBZ_dyn_spg_ts_sshf} was found to 291 287 be more conservative, and so is recommended. 292 288 293 %-------------------------------------------------------------294 289 % Filtered formulation 295 %------------------------------------------------------------- 296 \subsubsection[Filtered formulation (\texttt{\textbf{key\_dynspg\_flt}})] 297 {Filtered formulation (\protect\key{dynspg\_flt})} 298 \label{subsec:DYN_spg_flt} 290 %% ================================================================================================= 291 \subsubsection[Filtered formulation (\texttt{\textbf{key\_dynspg\_flt}})]{Filtered formulation (\protect\key{dynspg\_flt})} 292 \label{subsec:MBZ_dyn_spg_flt} 299 293 300 294 The filtered formulation follows the \citet{Roullet2000?} implementation. 301 295 The extra term introduced in the equations (see {\S}I.2.2) is solved implicitly. 302 296 The elliptic solvers available in the code are documented in \autoref{chap:MISC}. 303 The amplitude of the extra term is given by the namelist variable \np{rnu} .297 The amplitude of the extra term is given by the namelist variable \np{rnu}{rnu}. 304 298 The default value is 1, as recommended by \citet{Roullet2000?} 305 299 306 \colorbox{red}{\np{rnu}\forcode{ = 1} to be suppressed from namelist !} 307 308 %------------------------------------------------------------- 300 \colorbox{red}{\np[=1]{rnu}{rnu} to be suppressed from namelist !} 301 309 302 % Non-linear free surface formulation 310 %------------------------------------------------------------- 311 \subsection[Non-linear free surface formulation (\texttt{\textbf{key\_vvl}})] 312 {Non-linear free surface formulation (\protect\key{vvl})} 313 \label{subsec:DYN_spg_vvl} 303 %% ================================================================================================= 304 \subsection[Non-linear free surface formulation (\texttt{\textbf{key\_vvl}})]{Non-linear free surface formulation (\protect\key{vvl})} 305 \label{subsec:MBZ_dyn_spg_vvl} 314 306 315 307 In the non-linear free surface formulation, the variations of volume are fully taken into account. 316 308 This option is presented in a report \citep{levier.treguier.ea_rpt07} available on the \NEMO\ web site. 317 309 The three time-stepping methods (explicit, split-explicit and filtered) are the same as in 318 \autoref{ DYN_spg_linear} except that the ocean depth is now time-dependent.310 \autoref{?:DYN_spg_linear?} except that the ocean depth is now time-dependent. 319 311 In particular, this means that in filtered case, the matrix to be inverted has to be recomputed at each time-step. 320 312 321 \biblio 322 323 \pindex 313 \onlyinsubfile{\input{../../global/epilogue}} 324 314 325 315 \end{document}
Note: See TracChangeset
for help on using the changeset viewer.