Home Contact Us Site Map  
    next up previous contents
Next: 1.5 Appendix OCEAN Up: 1.4 Appendix ATMOSPHERE Previous: 1.4 Appendix ATMOSPHERE   Contents


1.4.1 Hydrostatic Primitive Equations for the Atmosphere in pressure coordinates

The hydrostatic primitive equations (HPEs) in p-coordinates are:

$\displaystyle \frac{D\vec{\mathbf{v}}_{h}}{Dt}+f\hat{\mathbf{k}}\times \vec{\mathbf{v}}
_{h}+\mathbf{\nabla }_{p}\phi$ $\displaystyle =$ $\displaystyle \vec{\mathbf{\mathcal{F}}}$ (1.45)
$\displaystyle \frac{\partial \phi }{\partial p}+\alpha$ $\displaystyle =$ 0 (1.46)
$\displaystyle \mathbf{\nabla }_{p}\cdot \vec{\mathbf{v}}_{h}+\frac{\partial \omega }{
\partial p}$ $\displaystyle =$ 0 (1.47)
$\displaystyle p\alpha$ $\displaystyle =$ $\displaystyle RT$ (1.48)
$\displaystyle c_{v}\frac{DT}{Dt}+p\frac{D\alpha }{Dt}$ $\displaystyle =$ $\displaystyle \mathcal{Q}$ (1.49)

where $ \vec{\mathbf{v}}_{h}=(u,v,0)$ is the `horizontal' (on pressure surfaces) component of velocity, $ \frac{D}{Dt}=\vec{\mathbf{v}}_{h}\cdot
\mathbf{\nabla }_{p}+\omega \frac{\partial }{\partial p}$ is the total derivative, $ f=2\Omega \sin \varphi$ is the Coriolis parameter, $ \phi =gz$ is the geopotential, $ \alpha =1/\rho $ is the specific volume, $ \omega =\frac{Dp
}{Dt}$ is the vertical velocity in the $ p-$coordinate. Equation(1.49) is the first law of thermodynamics where internal energy $ e=c_{v}T$, $ T$ is temperature, $ Q$ is the rate of heating per unit mass and $ p\frac{D\alpha }{Dt}$ is the work done by the fluid in compressing.

It is convenient to cast the heat equation in terms of potential temperature $ \theta $ so that it looks more like a generic conservation law. Differentiating (1.48) we get:

$\displaystyle p\frac{D\alpha }{Dt}+\alpha \frac{Dp}{Dt}=R\frac{DT}{Dt}$    

which, when added to the heat equation (1.49) and using $ c_{p}=c_{v}+R$, gives:

$\displaystyle c_{p}\frac{DT}{Dt}-\alpha \frac{Dp}{Dt}=\mathcal{Q}$ (1.50)

Potential temperature is defined:

$\displaystyle \theta =T(\frac{p_{c}}{p})^{\kappa }$ (1.51)

where $ p_{c}$ is a reference pressure and $ \kappa =R/c_{p}$. For convenience we will make use of the Exner function $ \Pi (p)$ which defined by:

$\displaystyle \Pi (p)=c_{p}(\frac{p}{p_{c}})^{\kappa }$ (1.52)

The following relations will be useful and are easily expressed in terms of the Exner function:

$\displaystyle c_{p}T=\Pi \theta \;\;;\;\;\frac{\partial \Pi }{\partial p}=\frac...
...\theta \;\;;\;\;\frac{D\Pi }{Dt}=\frac{\partial \Pi }{\partial p} \frac{Dp}{Dt}$    

where $ b=\frac{\partial \ \Pi }{\partial p}\theta $ is the buoyancy.

The heat equation is obtained by noting that

$\displaystyle c_{p}\frac{DT}{Dt}=\frac{D(\Pi \theta )}{Dt}=\Pi \frac{D\theta }{Dt}+\theta \frac{D\Pi }{Dt}=\Pi \frac{D\theta }{Dt}+\alpha \frac{Dp}{Dt}$    

and on substituting into (1.50) gives:

$\displaystyle \Pi \frac{D\theta }{Dt}=\mathcal{Q}$ (1.53)

which is in conservative form.

For convenience in the model we prefer to step forward (1.53) rather than (1.49). Boundary conditions

The upper and lower boundary conditions are :

at the top:$\displaystyle \;\;p=0$   $\displaystyle \omega =\frac{Dp}{Dt}=0$ (1.54)
at the surface:$\displaystyle \;\;p=p_{s}$   $\displaystyle \phi =\phi _{topo}=g~Z_{topo}$ (1.55)

In $ p$-coordinates, the upper boundary acts like a solid boundary ($ \omega
=0 $); in $ z$-coordinates and the lower boundary is analogous to a free surface ($ \phi $ is imposed and $ \omega \neq 0$). Splitting the geo-potential

For the purposes of initialization and reducing round-off errors, the model deals with perturbations from reference (or ``standard'') profiles. For example, the hydrostatic geopotential associated with the resting atmosphere is not dynamically relevant and can therefore be subtracted from the equations. The equations written in terms of perturbations are obtained by substituting the following definitions into the previous model equations:

$\displaystyle \theta$ $\displaystyle =$ $\displaystyle \theta _{o}+\theta ^{\prime }$ (1.56)
$\displaystyle \alpha$ $\displaystyle =$ $\displaystyle \alpha _{o}+\alpha ^{\prime }$ (1.57)
$\displaystyle \phi$ $\displaystyle =$ $\displaystyle \phi _{o}+\phi ^{\prime }$ (1.58)

The reference state (indicated by subscript ``0'') corresponds to horizontally homogeneous atmosphere at rest ( $ \theta _{o},\alpha _{o},\phi
_{o}$) with surface pressure $ p_{o}(x,y)$ that satisfies $ \phi
_{o}(p_{o})=g~Z_{topo}$, defined:
$\displaystyle \theta _{o}(p)$ $\displaystyle =$ $\displaystyle f^{n}(p)$  
$\displaystyle \alpha _{o}(p)$ $\displaystyle =$ $\displaystyle \Pi _{p}\theta _{o}$  
$\displaystyle \phi _{o}(p)$ $\displaystyle =$ $\displaystyle \phi _{topo}-\int_{p_{0}}^{p}\alpha _{o}dp$  

The final form of the HPE's in p coordinates is then:

$\displaystyle \frac{D\vec{\mathbf{v}}_{h}}{Dt}+f\hat{\mathbf{k}}\times \vec{\mathbf{v}}
_{h}+\mathbf{\nabla }_{p}\phi ^{\prime }$ $\displaystyle =$ $\displaystyle \vec{\mathbf{\mathcal{F}}}$ (1.59)
$\displaystyle \frac{\partial \phi ^{\prime }}{\partial p}+\alpha ^{\prime }$ $\displaystyle =$ 0 (1.60)
$\displaystyle \mathbf{\nabla }_{p}\cdot \vec{\mathbf{v}}_{h}+\frac{\partial \omega }{
\partial p}$ $\displaystyle =$ 0 (1.61)
$\displaystyle \frac{\partial \Pi }{\partial p}\theta ^{\prime }$ $\displaystyle =$ $\displaystyle \alpha ^{\prime }$ (1.62)
$\displaystyle \frac{D\theta }{Dt}$ $\displaystyle =$ $\displaystyle \frac{\mathcal{Q}}{\Pi }$ (1.63)

next up previous contents
Next: 1.5 Appendix OCEAN Up: 1.4 Appendix ATMOSPHERE Previous: 1.4 Appendix ATMOSPHERE   Contents
Copyright 2002 Massachusetts Institute of Technology