Changeset 11335 for NEMO/trunk/doc
- Timestamp:
- 2019-07-24T12:16:18+02:00 (5 years ago)
- Location:
- NEMO/trunk/doc/latex/NEMO
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/doc/latex/NEMO/main/bibliography.bib
r11333 r11335 2991 2991 } 2992 2992 2993 @article{ white.hoskins.ea_QJRMS05, 2994 title = "Consistent approximate models of the global atmosphere: shallow, deep, 2995 hydrostatic, quasi-hydrostatic and non-hydrostatic", 2996 pages = "2081--2107", 2997 journal = "Quarterly Journal of the Royal Meteorological Society", 2998 volume = "131", 2999 author = "A. A. White and B. J. Hoskins and I. Roulstone and A. Staniforth", 3000 year = "2005", 3001 doi = "10.1256/qj.04.49" 3002 } 3003 2993 3004 @article{ white.adcroft.ea_JCP09, 2994 3005 title = "High-order regridding-remapping schemes for continuous -
NEMO/trunk/doc/latex/NEMO/subfiles/annex_A.tex
r11151 r11335 29 29 \begin{equation} 30 30 \label{apdx:A_s_slope} 31 \sigma_1 =\frac{1}{e_1 } \;\left. {\frac{\partial z}{\partial i}} \right|_s31 \sigma_1 =\frac{1}{e_1 } \; \left. {\frac{\partial z}{\partial i}} \right|_s 32 32 \quad \text{and} \quad 33 \sigma_2 =\frac{1}{e_2 }\;\left. {\frac{\partial z}{\partial j}} \right|_s 34 \end{equation} 35 36 The chain rule to establish the model equations in the curvilinear $s-$coordinate system is: 33 \sigma_2 =\frac{1}{e_2 } \; \left. {\frac{\partial z}{\partial j}} \right|_s . 34 \end{equation} 35 36 The model fields (e.g. pressure $p$) can be viewed as functions of $(i,j,z,t)$ (e.g. $p(i,j,z,t)$) or as 37 functions of $(i,j,s,t)$ (e.g. $p(i,j,s,t)$). The symbol $\bullet$ will be used to represent any one of 38 these fields. Any ``infinitesimal'' change in $\bullet$ can be written in two forms: 39 \begin{equation} 40 \label{apdx:A_s_infin_changes} 41 \begin{aligned} 42 & \delta \bullet = \delta i \left. \frac{ \partial \bullet }{\partial i} \right|_{j,s,t} 43 + \delta j \left. \frac{ \partial \bullet }{\partial i} \right|_{i,s,t} 44 + \delta s \left. \frac{ \partial \bullet }{\partial s} \right|_{i,j,t} 45 + \delta t \left. \frac{ \partial \bullet }{\partial t} \right|_{i,j,s} , \\ 46 & \delta \bullet = \delta i \left. \frac{ \partial \bullet }{\partial i} \right|_{j,z,t} 47 + \delta j \left. \frac{ \partial \bullet }{\partial i} \right|_{i,z,t} 48 + \delta z \left. \frac{ \partial \bullet }{\partial z} \right|_{i,j,t} 49 + \delta t \left. \frac{ \partial \bullet }{\partial t} \right|_{i,j,z} . 50 \end{aligned} 51 \end{equation} 52 Using the first form and considering a change $\delta i$ with $j, z$ and $t$ held constant, shows that 53 \begin{equation} 54 \label{apdx:A_s_chain_rule} 55 \left. {\frac{\partial \bullet }{\partial i}} \right|_{j,z,t} = 56 \left. {\frac{\partial \bullet }{\partial i}} \right|_{j,s,t} 57 + \left. {\frac{\partial s }{\partial i}} \right|_{j,z,t} \; 58 \left. {\frac{\partial \bullet }{\partial s}} \right|_{i,j,t} . 59 \end{equation} 60 The term $\left. \partial s / \partial i \right|_{j,z,t}$ can be related to the slope of constant $s$ surfaces, 61 (\autoref{apdx:A_s_slope}), by applying the second of (\autoref{apdx:A_s_infin_changes}) with $\bullet$ set to 62 $s$ and $j, t$ held constant 63 \begin{equation} 64 \label{apdx:a_delta_s} 65 \delta s|_{j,t} = 66 \delta i \left. \frac{ \partial s }{\partial i} \right|_{j,z,t} 67 + \delta z \left. \frac{ \partial s }{\partial z} \right|_{i,j,t} . 68 \end{equation} 69 Choosing to look at a direction in the $(i,z)$ plane in which $\delta s = 0$ and using 70 (\autoref{apdx:A_s_slope}) we obtain 71 \begin{equation} 72 \left. \frac{ \partial s }{\partial i} \right|_{j,z,t} = 73 - \left. \frac{ \partial z }{\partial i} \right|_{j,s,t} \; 74 \left. \frac{ \partial s }{\partial z} \right|_{i,j,t} 75 = - \frac{e_1 }{e_3 }\sigma_1 . 76 \label{apdx:a_ds_di_z} 77 \end{equation} 78 Another identity, similar in form to (\autoref{apdx:a_ds_di_z}), can be derived 79 by choosing $\bullet$ to be $s$ and using the second form of (\autoref{apdx:A_s_infin_changes}) to consider 80 changes in which $i , j$ and $s$ are constant. This shows that 81 \begin{equation} 82 \label{apdx:A_w_in_s} 83 w_s = \left. \frac{ \partial z }{\partial t} \right|_{i,j,s} = 84 - \left. \frac{ \partial z }{\partial s} \right|_{i,j,t} 85 \left. \frac{ \partial s }{\partial t} \right|_{i,j,z} 86 = - e_3 \left. \frac{ \partial s }{\partial t} \right|_{i,j,z} . 87 \end{equation} 88 89 In what follows, for brevity, indication of the constancy of the $i, j$ and $t$ indices is 90 usually omitted. Using the arguments outlined above one can show that the chain rules needed to establish 91 the model equations in the curvilinear $s-$coordinate system are: 37 92 \begin{equation} 38 93 \label{apdx:A_s_chain_rule} 39 94 \begin{aligned} 40 95 &\left. {\frac{\partial \bullet }{\partial t}} \right|_z = 41 \left. {\frac{\partial \bullet }{\partial t}} \right|_s 42 -\frac{\partial \bullet }{\partial s}\;\frac{\partial s}{\partial t}\\96 \left. {\frac{\partial \bullet }{\partial t}} \right|_s 97 + \frac{\partial \bullet }{\partial s}\; \frac{\partial s}{\partial t} , \\ 43 98 &\left. {\frac{\partial \bullet }{\partial i}} \right|_z = 44 99 \left. {\frac{\partial \bullet }{\partial i}} \right|_s 45 -\frac{\partial \bullet }{\partial s}\;\frac{\partial s}{\partial i}=46 \left. {\frac{\partial \bullet }{\partial i}} \right|_s 47 -\frac{e_1 }{e_3 }\sigma_1 \frac{\partial \bullet }{\partial s} \\100 +\frac{\partial \bullet }{\partial s}\; \frac{\partial s}{\partial i}= 101 \left. {\frac{\partial \bullet }{\partial i}} \right|_s 102 -\frac{e_1 }{e_3 }\sigma_1 \frac{\partial \bullet }{\partial s} , \\ 48 103 &\left. {\frac{\partial \bullet }{\partial j}} \right|_z = 49 \left. {\frac{\partial \bullet }{\partial j}} \right|_s 50 -\frac{\partial \bullet }{\partial s}\;\frac{\partial s}{\partial j}=51 \left. {\frac{\partial \bullet }{\partial j}} \right|_s 52 - \frac{e_2 }{e_3 }\sigma_2 \frac{\partial \bullet }{\partial s} \\53 &\;\frac{\partial \bullet }{\partial z} \;\; = \frac{1}{e_3 }\frac{\partial \bullet }{\partial s} 104 \left. {\frac{\partial \bullet }{\partial j}} \right|_s 105 + \frac{\partial \bullet }{\partial s}\;\frac{\partial s}{\partial j}= 106 \left. {\frac{\partial \bullet }{\partial j}} \right|_s 107 - \frac{e_2 }{e_3 }\sigma_2 \frac{\partial \bullet }{\partial s} , \\ 108 &\;\frac{\partial \bullet }{\partial z} \;\; = \frac{1}{e_3 }\frac{\partial \bullet }{\partial s} . 54 109 \end{aligned} 55 110 \end{equation} 56 111 57 In particular applying the time derivative chain rule to $z$ provides the expression for $w_s$,58 the vertical velocity of the $s-$surfaces referenced to a fix z-coordinate:59 \begin{equation}60 \label{apdx:A_w_in_s}61 w_s = \left. \frac{\partial z }{\partial t} \right|_s62 = \frac{\partial z}{\partial s} \; \frac{\partial s}{\partial t}63 = e_3 \, \frac{\partial s}{\partial t}64 \end{equation}65 112 66 113 % ================================================================ … … 131 178 \end{subequations} 132 179 133 Here, $w$ is the vertical velocity relative to the $z-$coordinate system. 134 Introducing the dia-surface velocity component, 135 $\omega $, defined as the volume flux across the moving $s$-surfaces per unit horizontal area: 180 Here, $w$ is the vertical velocity relative to the $z-$coordinate system. 181 Using the first form of (\autoref{apdx:A_s_infin_changes}) 182 and the definitions (\autoref{apdx:A_s_slope}) and (\autoref{apdx:A_w_in_s}) for $\sigma_1$, $\sigma_2$ and $w_s$, 183 one can show that the vertical velocity, $w_p$ of a point 184 moving with the horizontal velocity of the fluid along an $s$ surface is given by 185 \begin{equation} 186 \label{apdx:A_w_p} 187 \begin{split} 188 w_p = & \left. \frac{ \partial z }{\partial t} \right|_s 189 + \frac{u}{e_1} \left. \frac{ \partial z }{\partial i} \right|_s 190 + \frac{v}{e_2} \left. \frac{ \partial z }{\partial j} \right|_s \\ 191 = & w_s + u \sigma_1 + v \sigma_2 . 192 \end{split} 193 \end{equation} 194 The vertical velocity across this surface is denoted by 136 195 \begin{equation} 137 196 \label{apdx:A_w_s} 138 \omega = w - w_s - \sigma_1 \,u - \sigma_2 \,v \\ 139 \end{equation} 140 with $w_s$ given by \autoref{apdx:A_w_in_s}, 141 we obtain the expression for the divergence of the velocity in the curvilinear $s-$coordinate system: 142 \begin{subequations} 143 \begin{align*} 144 { 145 \begin{array}{*{20}l} 146 \nabla \cdot {\mathrm {\mathbf U}} 147 &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ 197 \omega = w - w_p = w - ( w_s + \sigma_1 \,u + \sigma_2 \,v ) . 198 \end{equation} 199 Hence 200 \begin{equation} 201 \frac{1}{e_3 } \frac{\partial}{\partial s} \left[ w - u\;\sigma_1 - v\;\sigma_2 \right] = 202 \frac{1}{e_3 } \frac{\partial}{\partial s} \left[ \omega + w_s \right] = 203 \frac{1}{e_3 } \left[ \frac{\partial \omega}{\partial s} 204 + \left. \frac{ \partial }{\partial t} \right|_s \frac{\partial z}{\partial s} \right] = 205 \frac{1}{e_3 } \frac{\partial \omega}{\partial s} + \frac{1}{e_3 } \left. \frac{ \partial e_3}{\partial t} . \right|_s 206 \end{equation} 207 208 Using (\autoref{apdx:A_w_s}) in our expression for $\nabla \cdot {\mathrm {\mathbf U}}$ we obtain 209 our final expression for the divergence of the velocity in the curvilinear $s-$coordinate system: 210 \begin{equation} 211 \nabla \cdot {\mathrm {\mathbf U}} = 212 \frac{1}{e_1 \,e_2 \,e_3 } \left[ 148 213 \left. \frac{\partial (e_2 \,e_3 \,u)}{\partial i} \right|_s 149 214 +\left. \frac{\partial (e_1 \,e_3 \,v)}{\partial j} \right|_s \right] 150 215 + \frac{1}{e_3 } \frac{\partial \omega }{\partial s} 151 + \frac{1}{e_3 } \frac{\partial w_s }{\partial s} \\ \\ 152 &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ 153 \left. \frac{\partial (e_2 \,e_3 \,u)}{\partial i} \right|_s 154 +\left. \frac{\partial (e_1 \,e_3 \,v)}{\partial j} \right|_s \right] 155 + \frac{1}{e_3 } \frac{\partial \omega }{\partial s} 156 + \frac{1}{e_3 } \frac{\partial}{\partial s} \left( e_3 \; \frac{\partial s}{\partial t} \right) \\ \\ 157 &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ 158 \left. \frac{\partial (e_2 \,e_3 \,u)}{\partial i} \right|_s 159 +\left. \frac{\partial (e_1 \,e_3 \,v)}{\partial j} \right|_s \right] 160 + \frac{1}{e_3 } \frac{\partial \omega }{\partial s} 161 + \frac{\partial}{\partial s} \frac{\partial s}{\partial t} 162 + \frac{1}{e_3 } \frac{\partial s}{\partial t} \frac{\partial e_3}{\partial s} \\ \\ 163 &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ 164 \left. \frac{\partial (e_2 \,e_3 \,u)}{\partial i} \right|_s 165 +\left. \frac{\partial (e_1 \,e_3 \,v)}{\partial j} \right|_s \right] 166 + \frac{1}{e_3 } \frac{\partial \omega }{\partial s} 167 + \frac{1}{e_3 } \frac{\partial e_3}{\partial t} 168 \end{array} 169 } 170 \end{align*} 171 \end{subequations} 216 + \frac{1}{e_3 } \left. \frac{\partial e_3}{\partial t} \right|_s . 217 \end{equation} 172 218 173 219 As a result, the continuity equation \autoref{eq:PE_continuity} in the $s-$coordinates is: … … 178 224 {\left. {\frac{\partial (e_2 \,e_3 \,u)}{\partial i}} \right|_s 179 225 + \left. {\frac{\partial (e_1 \,e_3 \,v)}{\partial j}} \right|_s } \right] 180 +\frac{1}{e_3 }\frac{\partial \omega }{\partial s} = 0 181 \end{equation} 182 A additional term has appeared that takeinto account226 +\frac{1}{e_3 }\frac{\partial \omega }{\partial s} = 0 . 227 \end{equation} 228 An additional term has appeared that takes into account 183 229 the contribution of the time variation of the vertical coordinate to the volume budget. 184 230 … … 210 256 + w \;\frac{\partial u}{\partial z} \\ \\ 211 257 &= \left. {\frac{\partial u }{\partial t}} \right|_z 212 - \left. \zeta \right|_z v 213 + \frac{1}{e_1 \,e_2 }\left[ { \left.{ \frac{\partial (e_2 \,v)}{\partial i} }\right|_z 258 - \frac{1}{e_1 \,e_2 }\left[ { \left.{ \frac{\partial (e_2 \,v)}{\partial i} }\right|_z 214 259 -\left.{ \frac{\partial (e_1 \,u)}{\partial j} }\right|_z } \right] \; v 215 260 + \frac{1}{2e_1} \left.{ \frac{\partial (u^2+v^2)}{\partial i} } \right|_z … … 230 275 } \\ \\ 231 276 &= \left. {\frac{\partial u }{\partial t}} \right|_z 232 +\left. \zeta \right|_s \;v277 - \left. \zeta \right|_s \;v 233 278 + \frac{1}{2\,e_1}\left. {\frac{\partial (u^2+v^2)}{\partial i}} \right|_s \\ 234 279 &\qquad \qquad \qquad \quad 235 280 + \frac{w}{e_3 } \;\frac{\partial u}{\partial s} 236 -\left[ {\frac{\sigma_1 }{e_3 }\frac{\partial v}{\partial s}281 + \left[ {\frac{\sigma_1 }{e_3 }\frac{\partial v}{\partial s} 237 282 - \frac{\sigma_2 }{e_3 }\frac{\partial u}{\partial s}} \right]\;v 238 283 - \frac{\sigma_1 }{2e_3 }\frac{\partial (u^2+v^2)}{\partial s} \\ \\ 239 284 &= \left. {\frac{\partial u }{\partial t}} \right|_z 240 +\left. \zeta \right|_s \;v285 - \left. \zeta \right|_s \;v 241 286 + \frac{1}{2\,e_1}\left. {\frac{\partial (u^2+v^2)}{\partial i}} \right|_s \\ 242 287 &\qquad \qquad \qquad \quad … … 245 290 - \sigma_1 u\frac{\partial u}{\partial s} - \sigma_1 v\frac{\partial v}{\partial s}} \right] \\ \\ 246 291 &= \left. {\frac{\partial u }{\partial t}} \right|_z 247 +\left. \zeta \right|_s \;v292 - \left. \zeta \right|_s \;v 248 293 + \frac{1}{2\,e_1}\left. {\frac{\partial (u^2+v^2)}{\partial i}} \right|_s 249 294 + \frac{1}{e_3} \left[ w - \sigma_2 v - \sigma_1 u \right] 250 \; \frac{\partial u}{\partial s} \\295 \; \frac{\partial u}{\partial s} . \\ 251 296 % 252 \intertext{Introducing $\omega$, the dia- a-surface velocity given by (\autoref{apdx:A_w_s}) }297 \intertext{Introducing $\omega$, the dia-s-surface velocity given by (\autoref{apdx:A_w_s}) } 253 298 % 254 299 &= \left. {\frac{\partial u }{\partial t}} \right|_z 255 +\left. \zeta \right|_s \;v300 - \left. \zeta \right|_s \;v 256 301 + \frac{1}{2\,e_1}\left. {\frac{\partial (u^2+v^2)}{\partial i}} \right|_s 257 + \frac{1}{e_3 } \left( \omega -w_s \right) \frac{\partial u}{\partial s} \\302 + \frac{1}{e_3 } \left( \omega + w_s \right) \frac{\partial u}{\partial s} \\ 258 303 \end{array} 259 304 } … … 266 311 { 267 312 \begin{array}{*{20}l} 268 w_s\;\frac{\partial u}{\partial s}269 = \frac{\partial s}{\partial t}\; \frac{\partial u }{\partial s}270 = \left. {\frac{\partial u }{\partial t}} \right|_s - \left. {\frac{\partial u }{\partial t}} \right|_z \ quad ,313 \frac{w_s}{e_3} \;\frac{\partial u}{\partial s} 314 = - \left. \frac{\partial s}{\partial t} \right|_z \; \frac{\partial u }{\partial s} 315 = \left. {\frac{\partial u }{\partial t}} \right|_s - \left. {\frac{\partial u }{\partial t}} \right|_z \ . 271 316 \end{array} 272 317 } 273 318 \] 274 leads to the $s-$coordinate formulation of the total $z-$coordinate time derivative,319 This leads to the $s-$coordinate formulation of the total $z-$coordinate time derivative, 275 320 \ie the total $s-$coordinate time derivative : 276 321 \begin{align} … … 278 323 \left. \frac{D u}{D t} \right|_s 279 324 = \left. {\frac{\partial u }{\partial t}} \right|_s 280 +\left. \zeta \right|_s \;v325 - \left. \zeta \right|_s \;v 281 326 + \frac{1}{2\,e_1}\left. {\frac{\partial (u^2+v^2)}{\partial i}} \right|_s 282 + \frac{1}{e_3 } \omega \;\frac{\partial u}{\partial s} 327 + \frac{1}{e_3 } \omega \;\frac{\partial u}{\partial s} . 283 328 \end{align} 284 329 Therefore, the vector invariant form of the total time derivative has exactly the same mathematical form in … … 306 351 + \frac{1}{e_3} \frac{\partial \omega}{\partial s} \right] \\ \\ 307 352 &&- \frac{v}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 308 -u \;\frac{\partial e_1 }{\partial j} \right) \\353 -u \;\frac{\partial e_1 }{\partial j} \right) . \\ 309 354 \end{array} 310 355 } … … 328 373 - e_2 u \;\frac{\partial e_3 }{\partial i} 329 374 - e_1 v \;\frac{\partial e_3 }{\partial j} \right) 330 -\frac{1}{e_3} \frac{\partial \omega}{\partial s} \right] \\ \\375 + \frac{1}{e_3} \frac{\partial \omega}{\partial s} \right] \\ \\ 331 376 && - \frac{v}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 332 377 -u \;\frac{\partial e_1 }{\partial j} \right) \\ \\ … … 337 382 && - \,u \left[ \frac{1}{e_1 e_2 e_3} \left( \frac{\partial(e_2 e_3 \, u)}{\partial i} 338 383 + \frac{\partial(e_1 e_3 \, v)}{\partial j} \right) 339 -\frac{1}{e_3} \frac{\partial \omega}{\partial s} \right]384 + \frac{1}{e_3} \frac{\partial \omega}{\partial s} \right] 340 385 - \frac{v}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 341 -u \;\frac{\partial e_1 }{\partial j} \right) \\386 -u \;\frac{\partial e_1 }{\partial j} \right) . \\ 342 387 % 343 388 \intertext {Introducing a more compact form for the divergence of the momentum fluxes, … … 361 406 + \left. \nabla \cdot \left( {{\mathrm {\mathbf U}}\,u} \right) \right|_s 362 407 - \frac{v}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 363 -u \;\frac{\partial e_1 }{\partial j} \right) 408 -u \;\frac{\partial e_1 }{\partial j} \right). 364 409 \end{flalign} 365 410 which is the total time derivative expressed in the curvilinear $s-$coordinate system. … … 377 422 & =-\frac{1}{\rho_o e_1 }\left[ {\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{e_1 }{e_3 }\sigma_1 \frac{\partial p}{\partial s}} \right] \\ 378 423 & =-\frac{1}{\rho_o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s +\frac{\sigma_1 }{\rho_o \,e_3 }\left( {-g\;\rho \;e_3 } \right) \\ 379 &=-\frac{1}{\rho_o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{g\;\rho }{\rho_o }\sigma_1 424 &=-\frac{1}{\rho_o \,e_1 }\left. {\frac{\partial p}{\partial i}} \right|_s -\frac{g\;\rho }{\rho_o }\sigma_1 . 380 425 \end{split} 381 426 \] 382 427 Applying similar manipulation to the second component and 383 replacing $\sigma_1$ and $\sigma_2$ by their expression \autoref{apdx:A_s_slope}, it comes:428 replacing $\sigma_1$ and $\sigma_2$ by their expression \autoref{apdx:A_s_slope}, it becomes: 384 429 \begin{equation} 385 430 \label{apdx:A_grad_p_1} … … 391 436 -\frac{1}{\rho_o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z 392 437 &=-\frac{1}{\rho_o \,e_2 } \left( \left. {\frac{\partial p}{\partial j}} \right|_s 393 + g\;\rho \;\left. {\frac{\partial z}{\partial j}} \right|_s \right) \\438 + g\;\rho \;\left. {\frac{\partial z}{\partial j}} \right|_s \right) . \\ 394 439 \end{split} 395 440 \end{equation} … … 405 450 \[ 406 451 \begin{split} 407 p &= g\; \int_z^\eta \rho \; e_3 \; dk = g\; \int_z^\eta \ left( \rho_o \,d + 1 \right) \; e_3 \; dk \\408 &= g \, \rho_o \; \int_z^\eta d \; e_3 \; dk + g \, \int_z^\eta e_3 \; dk452 p &= g\; \int_z^\eta \rho \; e_3 \; dk = g\; \int_z^\eta \rho_o \left( d + 1 \right) \; e_3 \; dk \\ 453 &= g \, \rho_o \; \int_z^\eta d \; e_3 \; dk + \rho_o g \, \int_z^\eta e_3 \; dk . 409 454 \end{split} 410 455 \] … … 412 457 \begin{equation} 413 458 \label{apdx:A_pressure} 414 p = \rho_o \; p_h' + g \, ( z + \eta)459 p = \rho_o \; p_h' + \rho_o \, g \, ( \eta - z ) 415 460 \end{equation} 416 461 and the hydrostatic pressure balance expressed in terms of $p_h'$ and $d$ is: 417 462 \[ 418 \frac{\partial p_h'}{\partial k} = - d \, g \, e_3 463 \frac{\partial p_h'}{\partial k} = - d \, g \, e_3 . 419 464 \] 420 465 421 466 Substituing \autoref{apdx:A_pressure} in \autoref{apdx:A_grad_p_1} and 422 using the definition of the density anomaly it comes theexpression in two parts:467 using the definition of the density anomaly it becomes an expression in two parts: 423 468 \begin{equation} 424 469 \label{apdx:A_grad_p_2} … … 426 471 -\frac{1}{\rho_o \, e_1 } \left. {\frac{\partial p}{\partial i}} \right|_z 427 472 &=-\frac{1}{e_1 } \left( \left. {\frac{\partial p_h'}{\partial i}} \right|_s 428 + g\; d \;\left. {\frac{\partial z}{\partial i}} \right|_s \right) - \frac{g}{e_1 } \frac{\partial \eta}{\partial i} \\473 + g\; d \;\left. {\frac{\partial z}{\partial i}} \right|_s \right) - \frac{g}{e_1 } \frac{\partial \eta}{\partial i} , \\ 429 474 % 430 475 -\frac{1}{\rho_o \, e_2 }\left. {\frac{\partial p}{\partial j}} \right|_z 431 476 &=-\frac{1}{e_2 } \left( \left. {\frac{\partial p_h'}{\partial j}} \right|_s 432 + g\; d \;\left. {\frac{\partial z}{\partial j}} \right|_s \right) - \frac{g}{e_2 } \frac{\partial \eta}{\partial j} \\477 + g\; d \;\left. {\frac{\partial z}{\partial j}} \right|_s \right) - \frac{g}{e_2 } \frac{\partial \eta}{\partial j} . \\ 433 478 \end{split} 434 479 \end{equation} … … 463 508 - \frac{1}{e_1 } \left( \frac{\partial p_h'}{\partial i} + g\; d \; \frac{\partial z}{\partial i} \right) 464 509 - \frac{g}{e_1 } \frac{\partial \eta}{\partial i} 465 + D_u^{\vect{U}} + F_u^{\vect{U}} 510 + D_u^{\vect{U}} + F_u^{\vect{U}} , 466 511 \end{multline} 467 512 \begin{multline} … … 473 518 - \frac{1}{e_2 } \left( \frac{\partial p_h'}{\partial j} + g\; d \; \frac{\partial z}{\partial j} \right) 474 519 - \frac{g}{e_2 } \frac{\partial \eta}{\partial j} 475 + D_v^{\vect{U}} + F_v^{\vect{U}} 520 + D_v^{\vect{U}} + F_v^{\vect{U}} . 476 521 \end{multline} 477 522 \end{subequations} … … 483 528 \label{apdx:A_PE_dyn_flux_u} 484 529 \frac{1}{e_3} \frac{\partial \left( e_3\,u \right) }{\partial t} = 485 \nabla \cdot \left( {{\mathrm {\mathbf U}}\,u} \right)530 - \nabla \cdot \left( {{\mathrm {\mathbf U}}\,u} \right) 486 531 + \left\{ {f + \frac{1}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 487 532 -u \;\frac{\partial e_1 }{\partial j} \right)} \right\} \,v \\ 488 533 - \frac{1}{e_1 } \left( \frac{\partial p_h'}{\partial i} + g\; d \; \frac{\partial z}{\partial i} \right) 489 534 - \frac{g}{e_1 } \frac{\partial \eta}{\partial i} 490 + D_u^{\vect{U}} + F_u^{\vect{U}} 535 + D_u^{\vect{U}} + F_u^{\vect{U}} , 491 536 \end{multline} 492 537 \begin{multline} … … 494 539 \frac{1}{e_3}\frac{\partial \left( e_3\,v \right) }{\partial t}= 495 540 - \nabla \cdot \left( {{\mathrm {\mathbf U}}\,v} \right) 496 +\left\{ {f + \frac{1}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i}541 - \left\{ {f + \frac{1}{e_1 e_2 }\left( v \;\frac{\partial e_2 }{\partial i} 497 542 -u \;\frac{\partial e_1 }{\partial j} \right)} \right\} \,u \\ 498 543 - \frac{1}{e_2 } \left( \frac{\partial p_h'}{\partial j} + g\; d \; \frac{\partial z}{\partial j} \right) 499 544 - \frac{g}{e_2 } \frac{\partial \eta}{\partial j} 500 + D_v^{\vect{U}} + F_v^{\vect{U}} 545 + D_v^{\vect{U}} + F_v^{\vect{U}} . 501 546 \end{multline} 502 547 \end{subequations} … … 505 550 \begin{equation} 506 551 \label{apdx:A_dyn_zph} 507 \frac{\partial p_h'}{\partial k} = - d \, g \, e_3 552 \frac{\partial p_h'}{\partial k} = - d \, g \, e_3 . 508 553 \end{equation} 509 554 … … 531 576 \left[ \frac{\partial }{\partial i} \left( {e_2 \,e_3 \;Tu} \right) 532 577 + \frac{\partial }{\partial j} \left( {e_1 \,e_3 \;Tv} \right) \right] \\ 533 +\frac{1}{e_3} \frac{\partial }{\partial k} \left( Tw \right)578 - \frac{1}{e_3} \frac{\partial }{\partial k} \left( Tw \right) 534 579 + D^{T} +F^{T} 535 580 \end{multline} 536 581 537 The expression for the advection term is a straight consequence of ( A.4),582 The expression for the advection term is a straight consequence of (\autoref{apdx:A_sco_Continuity}), 538 583 the expression of the 3D divergence in the $s-$coordinates established above. 539 584 -
NEMO/trunk/doc/latex/NEMO/subfiles/annex_B.tex
r11331 r11335 53 53 { 54 54 \begin{array}{*{20}l} 55 D^T=& \frac{1}{e_1\,e_2\,e_3 }\;\left[ {\ \ \ \ e_2\,e_3\,A^{lT} \;\left. 56 {\frac{\partial }{\partial i}\left( {\frac{1}{e_1}\;\left. {\frac{\partial T}{\partial i}} \right|_s -\frac{\sigma_1 }{e_3 }\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right. \\ 57 &\qquad \quad \ \ \ +e_1\,e_3\,A^{lT} \;\left. {\frac{\partial }{\partial j}\left( {\frac{1}{e_2 }\;\left. {\frac{\partial T}{\partial j}} \right|_s -\frac{\sigma_2 }{e_3 }\;\frac{\partial T}{\partial s}} \right)} \right|_s \\ 58 &\qquad \quad \ \ \ + e_1\,e_2\,A^{lT} \;\frac{\partial }{\partial s}\left( {-\frac{\sigma_1 }{e_1 }\;\left. {\frac{\partial T}{\partial i}} \right|_s -\frac{\sigma_2 }{e_2 }\;\left. {\frac{\partial T}{\partial j}} \right|_s } \right. 59 \left. {\left. {+\left( {\varepsilon +\sigma_1^2+\sigma_2 ^2} \right)\;\frac{1}{e_3 }\;\frac{\partial T}{\partial s}} \right)\;\;} \right] 55 D^T= \frac{1}{e_1\,e_2\,e_3 } & \left\{ \quad \quad \frac{\partial }{\partial i} \left. \left[ e_2\,e_3 \, A^{lT} 56 \left( \ \frac{1}{e_1}\; \left. \frac{\partial T}{\partial i} \right|_s 57 -\frac{\sigma_1 }{e_3 } \; \frac{\partial T}{\partial s} \right) \right] \right|_s \right. \\ 58 & \quad \ + \ \left. \frac{\partial }{\partial j} \left. \left[ e_1\,e_3 \, A^{lT} 59 \left( \ \frac{1}{e_2 }\; \left. \frac{\partial T}{\partial j} \right|_s 60 -\frac{\sigma_2 }{e_3 } \; \frac{\partial T}{\partial s} \right) \right] \right|_s \right. \\ 61 & \quad \ + \ \left. e_1\,e_2\, \frac{\partial }{\partial s} \left[ A^{lT} \; \left( 62 -\frac{\sigma_1 }{e_1 } \; \left. \frac{\partial T}{\partial i} \right|_s 63 -\frac{\sigma_2 }{e_2 } \; \left. \frac{\partial T}{\partial j} \right|_s 64 +\left( \varepsilon +\sigma_1^2+\sigma_2 ^2 \right) \; \frac{1}{e_3 } \; \frac{\partial T}{\partial s} \right) \; \right] \; \right\} . 60 65 \end{array} 61 66 } … … 90 95 { 91 96 \begin{array}{*{20}l} 92 \intertext{Noting that $\frac{1}{e_1} \left. \frac{\partial e_3 }{\partial i} \right|_s = \frac{\partial \sigma_1 }{\partial s}$, itbecomes:}97 \intertext{Noting that $\frac{1}{e_1} \left. \frac{\partial e_3 }{\partial i} \right|_s = \frac{\partial \sigma_1 }{\partial s}$, this becomes:} 93 98 % 94 & =\frac{1}{e_1\,e_2\,e_3 }\left[ {\left. {\;\;\;\frac{\partial }{\partial i}\left( {\frac{e_2\,e_3 }{e_1}\,A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right)} \right|_s \left. -\, {e_3 \frac{\partial }{\partial i}\left( {\frac{e_2 \,\sigma_1 }{e_3 }A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right. \\99 D^T & =\frac{1}{e_1\,e_2\,e_3 }\left[ {\left. {\;\;\;\frac{\partial }{\partial i}\left( {\frac{e_2\,e_3 }{e_1}\,A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right)} \right|_s \left. -\, {e_3 \frac{\partial }{\partial i}\left( {\frac{e_2 \,\sigma_1 }{e_3 }A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right. \\ 95 100 & \qquad \qquad \quad -e_2 A^{lT}\;\frac{\partial \sigma_1 }{\partial s}\left. {\frac{\partial T}{\partial i}} \right|_s -e_1 \,\sigma_1 \frac{\partial }{\partial s}\left( {\frac{e_2 }{e_1 }A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right) \\ 96 101 & \qquad \qquad \quad\shoveright{ \left. { +e_1 \,\sigma_1 \frac{\partial }{\partial s}\left( {\frac{e_2 \,\sigma_1 }{e_3 }A^{lT}\;\frac{\partial T}{\partial s}} \right)+\frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 }{e_3 }A^{vT}\;\frac{\partial T}{\partial s}} \right)\;\;\;} \right] }\\ … … 99 104 & \qquad \qquad \quad \left. {+\frac{e_2 \,\sigma_1 }{e_3}A^{lT}\;\frac{\partial T}{\partial s} \;\frac{\partial e_3 }{\partial i}} \right|_s -e_2 A^{lT}\;\frac{\partial \sigma_1 }{\partial s}\left. {\frac{\partial T}{\partial i}} \right|_s \\ 100 105 & \qquad \qquad \quad-e_2 \,\sigma_1 \frac{\partial}{\partial s}\left( {A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right)+\frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 \,\sigma_1 ^2}{e_3 }A^{lT}\;\frac{\partial T}{\partial s}} \right) \\ 101 & \qquad \qquad \quad\shoveright{ \left. {-\frac{\partial \left( {e_1 \,e_2 \,\sigma_1 } \right)}{\partial s} \left( {\frac{\sigma_1 }{e_3}A^{lT}\;\frac{\partial T}{\partial s}} \right) + \frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 }{e_3 }A^{vT}\;\frac{\partial T}{\partial s}} \right)\;\;\;} \right]} 106 & \qquad \qquad \quad\shoveright{ \left. {-\frac{\partial \left( {e_1 \,e_2 \,\sigma_1 } \right)}{\partial s} \left( {\frac{\sigma_1 }{e_3}A^{lT}\;\frac{\partial T}{\partial s}} \right) + \frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 }{e_3 }A^{vT}\;\frac{\partial T}{\partial s}} \right)\;\;\;} \right]} . 102 107 \end{array} 103 108 } \\ … … 105 110 \begin{array}{*{20}l} 106 111 % 107 \intertext{Using the same remark as just above, itbecomes:}112 \intertext{Using the same remark as just above, $D^T$ becomes:} 108 113 % 109 &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ {\left. {\;\;\;\frac{\partial }{\partial i} \left( {\frac{e_2 \,e_3 }{e_1 }A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s -e_2 \,\sigma_1 A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right.\;\;\; \\114 D^T &= \frac{1}{e_1 \,e_2 \,e_3 } \left[ {\left. {\;\;\;\frac{\partial }{\partial i} \left( {\frac{e_2 \,e_3 }{e_1 }A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s -e_2 \,\sigma_1 A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right.\;\;\; \\ 110 115 & \qquad \qquad \quad+\frac{e_1 \,e_2 \,\sigma_1 }{e_3 }A^{lT}\;\frac{\partial T}{\partial s}\;\frac{\partial \sigma_1 }{\partial s} - \frac {\sigma_1 }{e_3} A^{lT} \;\frac{\partial \left( {e_1 \,e_2 \,\sigma_1 } \right)}{\partial s}\;\frac{\partial T}{\partial s} \\ 111 116 & \qquad \qquad \quad-e_2 \left( {A^{lT}\;\frac{\partial \sigma_1 }{\partial s}\left. {\frac{\partial T}{\partial i}} \right|_s +\frac{\partial }{\partial s}\left( {\sigma_1 A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right)-\frac{\partial \sigma_1 }{\partial s}\;A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s } \right) \\ 112 & \qquad \qquad \quad\shoveright{\left. {+\frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 \,\sigma_1 ^2}{e_3 }A^{lT}\;\frac{\partial T}{\partial s}+\frac{e_1 \,e_2}{e_3 }A^{vT}\;\frac{\partial T}{\partial s}} \right)\;\;\;} \right] }117 & \qquad \qquad \quad\shoveright{\left. {+\frac{\partial }{\partial s}\left( {\frac{e_1 \,e_2 \,\sigma_1 ^2}{e_3 }A^{lT}\;\frac{\partial T}{\partial s}+\frac{e_1 \,e_2}{e_3 }A^{vT}\;\frac{\partial T}{\partial s}} \right)\;\;\;} \right] . } 113 118 \end{array} 114 119 } \\ … … 120 125 the third line reduces to a single vertical derivative, so it becomes:} 121 126 % 122 & =\frac{1}{e_1 \,e_2 \,e_3 }\left[ {\left. {\;\;\;\frac{\partial }{\partial i}\left( {\frac{e_2 \,e_3 }{e_1 }A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s -e_2 \,\sigma_1 \,A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right. \\127 D^T & =\frac{1}{e_1 \,e_2 \,e_3 }\left[ {\left. {\;\;\;\frac{\partial }{\partial i}\left( {\frac{e_2 \,e_3 }{e_1 }A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s -e_2 \,\sigma_1 \,A^{lT}\;\frac{\partial T}{\partial s}} \right)} \right|_s } \right. \\ 123 128 & \qquad \qquad \quad \shoveright{ \left. {+\frac{\partial }{\partial s}\left( {-e_2 \,\sigma_1 \,A^{lT}\;\left. {\frac{\partial T}{\partial i}} \right|_s +A^{lT}\frac{e_1 \,e_2 }{e_3 }\;\left( {\varepsilon +\sigma_1 ^2} \right)\frac{\partial T}{\partial s}} \right)\;\;\;} \right]} \\ 124 129 % 125 130 \intertext{In other words, the horizontal/vertical Laplacian operator in the ($i$,$s$) plane takes the following form:} 126 131 \end{array} 127 } \\ 132 } \\ 128 133 % 129 134 {\frac{1}{e_1\,e_2\,e_3}} … … 228 233 The isopycnal diffusion operator \autoref{apdx:B4}, 229 234 \autoref{apdx:B_ldfiso} conserves tracer quantity and dissipates its square. 230 The demonstration of the first property is trivial as \autoref{apdx:B4} is the divergence of fluxes. 231 Let us demonstrate the second one:235 As \autoref{apdx:B4} is the divergence of a flux, the demonstration of the first property is trivial, providing that the flux normal to the boundary is zero 236 (as it is when $A_h$ is zero at the boundary). Let us demonstrate the second one: 232 237 \[ 233 238 \iiint\limits_D T\;\nabla .\left( {\textbf{A}}_{\textbf{I}} \nabla T \right)\,dv … … 248 253 j}-a_2 \frac{\partial T}{\partial k}} \right)^2} 249 254 +\varepsilon \left(\frac{\partial T}{\partial k}\right) ^2\right] \\ 250 & \geq 0 255 & \geq 0 . 251 256 \end{array} 252 257 } … … 365 370 the third component of the second vector is obviously zero and thus : 366 371 \[ 367 \Delta {\textbf{U}}_h = \nabla _h \left( \chi \right) - \nabla _h \times \left( \zeta \ right) + \frac {1}{e_3 } \frac {\partial }{\partial k} \left( {\frac {1}{e_3 } \frac{\partial {\textbf{ U}}_h }{\partial k}} \right)372 \Delta {\textbf{U}}_h = \nabla _h \left( \chi \right) - \nabla _h \times \left( \zeta \textbf{k} \right) + \frac {1}{e_3 } \frac {\partial }{\partial k} \left( {\frac {1}{e_3 } \frac{\partial {\textbf{ U}}_h }{\partial k}} \right) . 368 373 \] 369 374 … … 381 386 - \nabla _h \times \left( {A^{lm}\;\zeta \;{\textbf{k}}} \right) 382 387 + \frac{1}{e_3 }\frac{\partial }{\partial k}\left( {\frac{A^{vm}\;}{e_3 } 383 \frac{\partial {\mathrm {\mathbf U}}_h }{\partial k}} \right) \\388 \frac{\partial {\mathrm {\mathbf U}}_h }{\partial k}} \right) , \\ 384 389 \end{equation} 385 390 that is, in expanded form: … … 388 393 & = \frac{1}{e_1} \frac{\partial \left( {A^{lm}\chi } \right)}{\partial i} 389 394 -\frac{1}{e_2} \frac{\partial \left( {A^{lm}\zeta } \right)}{\partial j} 390 +\frac{1}{e_3} \frac{\partial }{\partial k} \left( \frac{A^{vm}}{e_3} \frac{\partial u}{\partial k} \right)\\395 +\frac{1}{e_3} \frac{\partial }{\partial k} \left( \frac{A^{vm}}{e_3} \frac{\partial u}{\partial k} \right) , \\ 391 396 D^{\textbf{U}}_v 392 397 & = \frac{1}{e_2 }\frac{\partial \left( {A^{lm}\chi } \right)}{\partial j} 393 398 +\frac{1}{e_1 }\frac{\partial \left( {A^{lm}\zeta } \right)}{\partial i} 394 +\frac{1}{e_3} \frac{\partial }{\partial k} \left( \frac{A^{vm}}{e_3} \frac{\partial v}{\partial k} \right)399 +\frac{1}{e_3} \frac{\partial }{\partial k} \left( \frac{A^{vm}}{e_3} \frac{\partial v}{\partial k} \right) . 395 400 \end{align*} 396 401 -
NEMO/trunk/doc/latex/NEMO/subfiles/chap_model_basics.tex
r11151 r11335 32 32 \begin{enumerate} 33 33 \item 34 \textit{spherical earth approximation}: the geopotential surfaces are assumed to be spheres so that 35 gravity (local vertical) is parallel to the earth's radius 34 \textit{spherical Earth approximation}: the geopotential surfaces are assumed to be oblate spheriods 35 that follow the Earth's bulge; these spheroids are approximated by spheres with 36 gravity locally vertical (parallel to the Earth's radius) and independent of latitude 37 \citep[][section 2]{white.hoskins.ea_QJRMS05}. 36 38 \item 37 39 \textit{thin-shell approximation}: the ocean depth is neglected compared to the earth's radius … … 63 65 \nabla \cdot \vect U = 0 64 66 \end{equation} 67 \item 68 \textit{Neglect of additional Coriolis terms}: the Coriolis terms that vary with the cosine of latitude are neglected. 69 These terms may be non-negligible where the Brunt-Vaisala frequency $N$ is small, either in the deep ocean or 70 in the sub-mesoscale motions of the mixed layer, or near the equator \citep[][section 1]{white.hoskins.ea_QJRMS05}. 71 They can be consistently included as part of the ocean dynamics \citep[][section 3(d)]{white.hoskins.ea_QJRMS05} and are 72 retained in the MIT ocean model. 65 73 \end{enumerate} 66 74 67 75 Because the gravitational force is so dominant in the equations of large-scale motions, 68 it is useful to choose an orthogonal set of unit vectors $(i,j,k)$ linked to the earth such that76 it is useful to choose an orthogonal set of unit vectors $(i,j,k)$ linked to the Earth such that 69 77 $k$ is the local upward vector and $(i,j)$ are two vectors orthogonal to $k$, 70 78 \ie tangent to the geopotential surfaces. … … 107 115 an air-sea or ice-sea interface at its top. 108 116 These boundaries can be defined by two surfaces, $z = - H(i,j)$ and $z = \eta(i,j,k,t)$, 109 where $H$ is the depth of the ocean bottom and $\eta$ is the height of the sea surface. 110 Both $H$ and $\eta$ are usually referenced to a given surface, $z = 0$, chosen as a mean sea surface 117 where $H$ is the depth of the ocean bottom and $\eta$ is the height of the sea surface 118 (discretisation can introduce additional artificial ``side-wall'' boundaries). 119 Both $H$ and $\eta$ are referenced to a surface of constant geopotential (\ie a mean sea surface height) on which $z = 0$. 111 120 (\autoref{fig:ocean_bc}). 112 121 Through these two boundaries, the ocean can exchange fluxes of heat, fresh water, salt, and momentum with … … 210 219 The flow is barotropic and the surface moves up and down with gravity as the restoring force. 211 220 The phase speed of such waves is high (some hundreds of metres per second) so that 212 the time step would have to be very short if they were present in the model.221 the time step has to be very short when they are present in the model. 213 222 The latter strategy filters out these waves since the rigid lid approximation implies $\eta = 0$, 214 223 \ie the sea surface is the surface $z = 0$. … … 217 226 The rigid-lid hypothesis is an obsolescent feature in modern OGCMs. 218 227 It has been available until the release 3.1 of \NEMO, and it has been removed in release 3.2 and followings. 219 Only the free surface formulation is now described in th e this document (see the next sub-section).228 Only the free surface formulation is now described in this document (see the next sub-section). 220 229 221 230 % ------------------------------------------------------------------------------------------------------------- … … 237 246 Allowing the air-sea interface to move introduces the external gravity waves (EGWs) as 238 247 a class of solution of the primitive equations. 239 These waves are barotropic because of hydrostatic assumption,and their phase speed is quite high.248 These waves are barotropic (\ie nearly independent of depth) and their phase speed is quite high. 240 249 Their time scale is short with respect to the other processes described by the primitive equations. 241 250 … … 266 275 the implicit scheme \citep{dukowicz.smith_JGR94} or 267 276 the addition of a filtering force in the momentum equation \citep{roullet.madec_JGR00}. 268 With the present release, \NEMO offers the choice between277 With the present release, \NEMO offers the choice between 269 278 an explicit free surface (see \autoref{subsec:DYN_spg_exp}) or 270 279 a split-explicit scheme strongly inspired the one proposed by \citet{shchepetkin.mcwilliams_OM05} … … 338 347 the vertical scale factor is a single function of $k$ as $k$ is parallel to $z$. 339 348 The scalar and vector operators that appear in the primitive equations 340 (\autoref{eq:PE_dyn} to \autoref{eq:PE_eos}) can be written in the tensorial form,349 (\autoref{eq:PE_dyn} to \autoref{eq:PE_eos}) can then be written in the tensorial form, 341 350 invariant in any orthogonal horizontal curvilinear coordinate system transformation: 342 351 \begin{subequations} … … 384 393 \end{gather} 385 394 386 Using the fact that the horizontal scale factors $e_1$ and $e_2$ are independent of $k$ and that395 Using again the fact that the horizontal scale factors $e_1$ and $e_2$ are independent of $k$ and that 387 396 $e_3$ is a function of the single variable $k$, 388 397 $NLT$ the nonlinear term of \autoref{eq:PE_dyn} can be transformed as follows: … … 456 465 & &= &\nabla \cdot (\vect U \, u) - (\nabla \cdot \vect U) \ u 457 466 + \frac{1}{e_1 e_2} \lt( -v^2 \pd[e_2]{i} + u v \, \pd[e_1]{j} \rt) \\ 458 \intertext{as $\nabla \cdot {\vect U} \; = 0$ (incompressibility) it comes:}467 \intertext{as $\nabla \cdot {\vect U} \; = 0$ (incompressibility) it becomes:} 459 468 & &= &\, \nabla \cdot (\vect U \, u) + \frac{1}{e_1 e_2} \lt( v \; \pd[e_2]{i} - u \; \pd[e_1]{j} \rt) (-v) 460 469 \end{alignat*} … … 516 525 % \label{eq:PE_dyn_flux_v} 517 526 \pd[v]{t} = - \lt[ f + \frac{1}{e_1 \; e_2} \lt( v \pd[e_2]{i} - u \pd[e_1]{j} \rt) \rt] \, u \\ 518 +\frac{1}{e_1 \; e_2} \lt( \pd[(e_2 \, u \, v)]{i} + \pd[(e_1 \, v \, v)]{j} \rt) \\527 - \frac{1}{e_1 \; e_2} \lt( \pd[(e_2 \, u \, v)]{i} + \pd[(e_1 \, v \, v)]{j} \rt) \\ 519 528 - \frac{1}{e_3} \pd[(w \, v)]{k} - \frac{1}{e_2} \pd[]{j} \lt( \frac{p_s + p_h}{\rho_o} \rt) 520 529 + D_v^{\vect U} + F_v^{\vect U} … … 526 535 p_s = \rho \,g \, \eta 527 536 \] 528 with $\eta$ issolution of \autoref{eq:PE_ssh}.537 and $\eta$ is the solution of \autoref{eq:PE_ssh}. 529 538 530 539 The vertical velocity and the hydrostatic pressure are diagnosed from the following equations: … … 536 545 \] 537 546 where the divergence of the horizontal velocity, $\chi$ is given by \autoref{eq:PE_div_Uh}. 538 \item \textit{tracer equations}: 539 \[ 540 %\label{eq:S} 541 \pd[T]{t} = - \frac{1}{e_1 e_2} \lt[ \pd[(e_2 T \, u)]{i} + \pd[(e_1 T \, v)]{j} \rt] 547 548 \item 549 \textbf{tracer equations}: 550 \begin{equation} 551 \begin{split} 552 \pd[T]{t} = & - \frac{1}{e_1 e_2} \lt[ \pd[(e_2 T \, u)]{i} + \pd[(e_1 T \, v)]{j} \rt] 542 553 - \frac{1}{e_3} \pd[(T \, w)]{k} + D^T + F^T \\ 543 %\label{eq:T} 544 \pd[S]{t} = - \frac{1}{e_1 e_2} \lt[ \pd[(e_2 S \, u)]{i} + \pd[(e_1 S \, v)]{j} \rt] 545 - \frac{1}{e_3} \pd[(S \, w)]{k} + D^S + F^S 546 %\label{eq:rho} 547 \rho = \rho \big( T,S,z(k) \big) 548 \] 554 \pd[S]{t} = & - \frac{1}{e_1 e_2} \lt[ \pd[(e_2 S \, u)]{i} + \pd[(e_1 S \, v)]{j} \rt] 555 - \frac{1}{e_3} \pd[(S \, w)]{k} + D^S + F^S \\ 556 \rho = & \rho \big( T,S,z(k) \big) 557 \end{split} 558 \end{equation} 549 559 \end{itemize} 550 560 … … 575 585 follows the isopycnal surfaces, \eg an isopycnic coordinate. 576 586 577 In order to satisfy two or more constrain s one can even be tempted to mixed these coordinate systems, as in587 In order to satisfy two or more constraints one can even be tempted to mixed these coordinate systems, as in 578 588 HYCOM (mixture of $z$-coordinate at the surface, isopycnic coordinate in the ocean interior and $\sigma$ at 579 589 the ocean bottom) \citep{chassignet.smith.ea_JPO03} or … … 594 604 This so-called \textit{generalised vertical coordinate} \citep{kasahara_MWR74} is in fact 595 605 an Arbitrary Lagrangian--Eulerian (ALE) coordinate. 596 Indeed, choosing an expression for $s$ is an arbitrary choice thatdetermines606 Indeed, one has a great deal of freedom in the choice of expression for $s$. The choice determines 597 607 which part of the vertical velocity (defined from a fixed referential) will cross the levels (Eulerian part) and 598 608 which part will be used to move them (Lagrangian part). … … 601 611 Its most often used implementation is via an ALE algorithm, 602 612 in which a pure lagrangian step is followed by regridding and remapping steps, 603 the lat er step implicitly embedding the vertical advection613 the latter step implicitly embedding the vertical advection 604 614 \citep{hirt.amsden.ea_JCP74, chassignet.smith.ea_JPO03, white.adcroft.ea_JCP09}. 605 615 Here we follow the \citep{kasahara_MWR74} strategy: 606 a regridding step (an update of the vertical coordinate) followed by an eulerian step with616 a regridding step (an update of the vertical coordinate) followed by an Eulerian step with 607 617 an explicit computation of vertical advection relative to the moving s-surfaces. 608 618 609 619 %\gmcomment{ 610 620 %A key point here is that the $s$-coordinate depends on $(i,j)$ ==> horizontal pressure gradient... 611 the generalized vertical coordinates used in ocean modelling are not orthogonal,621 The generalized vertical coordinates used in ocean modelling are not orthogonal, 612 622 which contrasts with many other applications in mathematical physics. 613 623 Hence, it is useful to keep in mind the following properties that may seem odd on initial encounter. … … 615 625 The horizontal velocity in ocean models measures motions in the horizontal plane, 616 626 perpendicular to the local gravitational field. 617 That is, horizontal velocity is mathematically the same regardless the vertical coordinate, be it geopotential,627 That is, horizontal velocity is mathematically the same regardless of the vertical coordinate, be it geopotential, 618 628 isopycnal, pressure, or terrain following. 619 629 The key motivation for maintaining the same horizontal velocity component is that … … 660 670 \[ 661 671 % \label{eq:PE_sco_w} 662 \omega = w - e_3 \, \pd[s]{t}- \sigma_1 \, u - \sigma_2 \, v672 \omega = w - \, \lt. \pd[z]{t} \rt|_s - \sigma_1 \, u - \sigma_2 \, v 663 673 \] 664 674 … … 671 681 % \label{eq:PE_sco_u_vector} 672 682 \pd[u]{t} = + (\zeta + f) \, v - \frac{1}{2 \, e_1} \pd[]{i} (u^2 + v^2) - \frac{1}{e_3} \omega \pd[u]{k} \\ 673 - \frac{1}{e_1} \pd[]{i} \lt( \frac{p_s + p_h}{\rho_o} \rt) +g \frac{\rho}{\rho_o} \sigma_1683 - \frac{1}{e_1} \pd[]{i} \lt( \frac{p_s + p_h}{\rho_o} \rt) - g \frac{\rho}{\rho_o} \sigma_1 674 684 + D_u^{\vect U} + F_u^{\vect U} 675 685 \end{multline*} … … 677 687 % \label{eq:PE_sco_v_vector} 678 688 \pd[v]{t} = - (\zeta + f) \, u - \frac{1}{2 \, e_2} \pd[]{j}(u^2 + v^2) - \frac{1}{e_3} \omega \pd[v]{k} \\ 679 - \frac{1}{e_2} \pd[]{j} \lt( \frac{p_s + p_h}{\rho_o} \rt) +g \frac{\rho}{\rho_o} \sigma_2689 - \frac{1}{e_2} \pd[]{j} \lt( \frac{p_s + p_h}{\rho_o} \rt) - g \frac{\rho}{\rho_o} \sigma_2 680 690 + D_v^{\vect U} + F_v^{\vect U} 681 691 \end{multline*} … … 687 697 - \frac{1}{e_3} \pd[(\omega \, u)]{k} 688 698 - \frac{1}{e_1} \pd[]{i} \lt( \frac{p_s + p_h}{\rho_o} \rt) 689 +g \frac{\rho}{\rho_o} \sigma_1 + D_u^{\vect U} + F_u^{\vect U}699 - g \frac{\rho}{\rho_o} \sigma_1 + D_u^{\vect U} + F_u^{\vect U} 690 700 \end{multline*} 691 701 \begin{multline*} … … 695 705 - \frac{1}{e_3} \pd[(\omega \, v)]{k} 696 706 - \frac{1}{e_2} \pd[]{j} \lt( \frac{p_s + p_h}{\rho_o} \rt) 697 +g \frac{\rho}{\rho_o}\sigma_2 + D_v^{\vect U} + F_v^{\vect U}707 - g \frac{\rho}{\rho_o}\sigma_2 + D_v^{\vect U} + F_v^{\vect U} 698 708 \end{multline*} 699 709 where the relative vorticity, $\zeta$, the surface pressure gradient, … … 750 760 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 751 761 752 In th atcase, the free surface equation is nonlinear, and the variations of volume are fully taken into account.762 In this case, the free surface equation is nonlinear, and the variations of volume are fully taken into account. 753 763 These coordinates systems is presented in a report \citep{levier.treguier.ea_rpt07} available on the \NEMO web site. 754 764 … … 766 776 The position (\zstar) and vertical discretization (\zstar) are expressed as: 767 777 \[ 768 % \label{eq: z-star}778 % \label{eq:PE_z-star} 769 779 H + \zstar = (H + z) / r \quad \text{and} \quad \delta \zstar 770 780 = \delta z / r \quad \text{with} \quad r 771 = \frac{H + \eta}{H} 781 = \frac{H + \eta}{H} . 782 \] 783 Simple re-organisation of the above expressions gives 784 \[ 785 % \label{eq:PE_zstar_2} 786 \zstar = H \lt( \frac{z - \eta}{H + \eta} \rt) . 772 787 \] 773 788 Since the vertical displacement of the free surface is incorporated in the vertical coordinate \zstar, … … 776 791 Also the divergence of the flow field is no longer zero as shown by the continuity equation: 777 792 \[ 778 \pd[r]{t} = \nabla_{\zstar} \cdot \lt( r \; \vect U_h \rt) (r \; w *) = 0793 \pd[r]{t} = \nabla_{\zstar} \cdot \lt( r \; \vect U_h \rt) + \pd[r \; w^*]{\zstar} = 0 . 779 794 \] 780 781 % from MOM4p1 documentation 782 To overcome problems with vanishing surface and/or bottom cells, we consider the zstar coordinate 783 \[ 784 % \label{eq:PE_} 785 \zstar = H \lt( \frac{z - \eta}{H + \eta} \rt) 786 \] 787 788 This coordinate is closely related to the "eta" coordinate used in many atmospheric models 795 This \zstar coordinate is closely related to the "eta" coordinate used in many atmospheric models 789 796 (see Black (1994) for a review of eta coordinate atmospheric models). 790 797 It was originally used in ocean models by Stacey et al. (1995) for studies of tides next to shelves, … … 798 805 These properties greatly reduce difficulties of computing the horizontal pressure gradient relative to 799 806 terrain following sigma models discussed in \autoref{subsec:PE_sco}. 800 Additionally, since \zstarwhen $\eta = 0$,807 Additionally, since $\zstar = z$ when $\eta = 0$, 801 808 no flow is spontaneously generated in an unforced ocean starting from rest, regardless the bottom topography. 802 809 This behaviour is in contrast to the case with "s"-models, where pressure gradient errors in the presence of … … 804 811 depending on the sophistication of the pressure gradient solver. 805 812 The quasi -horizontal nature of the coordinate surfaces also facilitates the implementation of 806 neutral physics parameterizations in \zstar models using the same techniques as in $z$-models813 neutral physics parameterizations in \zstar models using the same techniques as in $z$-models 807 814 (see Chapters 13-16 of \cite{griffies_bk04}) for a discussion of neutral physics in $z$-models, 808 815 as well as \autoref{sec:LDF_slp} in this document for treatment in \NEMO). 809 816 810 The range over which \zstar varies is time independent $-H \leq \zstar \leq 0$.811 Hence, all cells remain nonvanishing, so long as the surface height maintains $\eta > ?H$.817 The range over which \zstar varies is time independent $-H \leq \zstar \leq 0$. 818 Hence, all cells remain nonvanishing, so long as the surface height maintains $\eta > -H$. 812 819 This is a minor constraint relative to that encountered on the surface height when using $s = z$ or $s = z - \eta$. 813 820 814 Because \zstar has a time independent range, all grid cells have static increments ds,815 and the sum of the ver 821 Because \zstar has a time independent range, all grid cells have static increments ds, 822 and the sum of the vertical increments yields the time independent ocean depth. %k ds = H. 816 823 The \zstar coordinate is therefore invisible to undulations of the free surface, 817 824 since it moves along with the free surface. 818 This proper ty means that no spurious vertical transport is induced across surfaces of constant \zstarby825 This property means that no spurious vertical transport is induced across surfaces of constant \zstar by 819 826 the motion of external gravity waves. 820 Such spurious transpor 821 Quite generally, the time independent range for the \zstar coordinate is a very convenient property that822 allows for a nearly arbitrary ver 827 Such spurious transport can be a problem in z-models, especially those with tidal forcing. 828 Quite generally, the time independent range for the \zstar coordinate is a very convenient property that 829 allows for a nearly arbitrary vertical resolution even in the presence of large amplitude fluctuations of 823 830 the surface height, again so long as $\eta > -H$. 824 831 %end MOM doc %%% … … 870 877 \begin{equation} 871 878 \label{eq:PE_p_sco} 872 \nabla p |_z = \nabla p |_s - \ pd[p]{s} \nabla z |_s879 \nabla p |_z = \nabla p |_s - \frac{1}{e_3} \pd[p]{s} \nabla z |_s 873 880 \end{equation} 874 881 875 882 The second term in \autoref{eq:PE_p_sco} depends on the tilt of the coordinate surface and 876 introducesa truncation error that is not present in a $z$-model.883 leads to a truncation error that is not present in a $z$-model. 877 884 In the special case of a $\sigma$-coordinate (i.e. a depth-normalised coordinate system $\sigma = z/H$), 878 885 \citet{haney_JPO91} and \citet{beckmann.haidvogel_JPO93} have given estimates of the magnitude of this truncation error. … … 887 894 However, the definition of the model domain vertical coordinate becomes then a non-trivial thing for 888 895 a realistic bottom topography: 889 a envelope topography is defined in $s$-coordinate on which a full or896 an envelope topography is defined in $s$-coordinate on which a full or 890 897 partial step bottom topography is then applied in order to adjust the model depth to the observed one 891 898 (see \autoref{sec:DOM_zgr}. … … 922 929 923 930 The \ztilde -coordinate has been developed by \citet{leclair.madec_OM11}. 924 It is available in \NEMO since the version 3.4 .931 It is available in \NEMO since the version 3.4 and is more robust in version 4.0 than previously. 925 932 Nevertheless, it is currently not robust enough to be used in all possible configurations. 926 933 Its use is therefore not recommended. … … 934 941 \label{sec:PE_zdf_ldf} 935 942 936 The primitive equations describe the behaviour of a geophysical fluid at space and time scales larger than943 The hydrostatic primitive equations describe the behaviour of a geophysical fluid at space and time scales larger than 937 944 a few kilometres in the horizontal, a few meters in the vertical and a few minutes. 938 945 They are usually solved at larger scales: the specified grid spacing and time step of the numerical model. … … 984 991 All the vertical physics is embedded in the specification of the eddy coefficients. 985 992 They can be assumed to be either constant, or function of the local fluid properties 986 (\eg Richardson number, Brunt-Vais\"{a}l\"{a} frequency ...),993 (\eg Richardson number, Brunt-Vais\"{a}l\"{a} frequency, distance from the boundary ...), 987 994 or computed from a turbulent closure model. 988 995 The choices available in \NEMO are discussed in \autoref{chap:ZDF}). … … 1016 1023 both horizontal and isoneutral operators have no effect on mean (\ie large scale) potential energy whereas 1017 1024 potential energy is a main source of turbulence (through baroclinic instabilities). 1018 \citet{gent.mcwilliams_JPO90} haveproposed a parameterisation of mesoscale eddy-induced turbulence which1025 \citet{gent.mcwilliams_JPO90} proposed a parameterisation of mesoscale eddy-induced turbulence which 1019 1026 associates an eddy-induced velocity to the isoneutral diffusion. 1020 1027 Its mean effect is to reduce the mean potential energy of the ocean. … … 1033 1040 Another approach is becoming more and more popular: 1034 1041 instead of specifying explicitly a sub-grid scale term in the momentum and tracer time evolution equations, 1035 one uses a advective scheme which is diffusive enough to maintain the model stability.1042 one uses an advective scheme which is diffusive enough to maintain the model stability. 1036 1043 It must be emphasised that then, all the sub-grid scale physics is included in the formulation of 1037 1044 the advection scheme. 1038 1045 1039 1046 All these parameterisations of subgrid scale physics have advantages and drawbacks. 1040 The reare not all available in \NEMO. For active tracers (temperature and salinity) the main ones are:1047 They are not all available in \NEMO. For active tracers (temperature and salinity) the main ones are: 1041 1048 Laplacian and bilaplacian operators acting along geopotential or iso-neutral surfaces, 1042 1049 \citet{gent.mcwilliams_JPO90} parameterisation, and various slightly diffusive advection schemes. … … 1141 1148 - \nabla_h \times \big( A^{lm} \, \zeta \; \vect k \big) \\ 1142 1149 &= \lt( \frac{1}{e_1} \pd[ \lt( A^{lm} \chi \rt) ]{i} \rt. 1143 - \frac{1}{e_2 e_3} \pd[ \lt( A^{lm} \; e_3 \zeta \rt) ]{j} 1150 - \frac{1}{e_2 e_3} \pd[ \lt( A^{lm} \; e_3 \zeta \rt) ]{j} , 1144 1151 \frac{1}{e_2} \pd[ \lt( A^{lm} \chi \rt) ]{j} 1145 1152 \lt. + \frac{1}{e_1 e_3} \pd[ \lt( A^{lm} \; e_3 \zeta \rt) ]{i} \rt) … … 1164 1171 a geographical coordinate system \citep{lengaigne.madec.ea_JGR03}. 1165 1172 1166 \subsubsection{ lateral bilaplacian momentum diffusive operator}1173 \subsubsection{Lateral bilaplacian momentum diffusive operator} 1167 1174 1168 1175 As for tracers, the bilaplacian order momentum diffusive operator is a re-entering Laplacian operator with
Note: See TracChangeset
for help on using the changeset viewer.