: 3 Turbulent parameterization : Two dimensional anelastic model : 1 Outline of discretizatoin

# 2 Atmospheric model

## 2.1 Equation of motion

Before making finite difference equations, equation (1)(3) show in Part I are transformed as follows.

where,

In this formulation, the time dependence of upper and lower boundary conditions are disappeared.

The advection terms are evaluated by the combination scheme of flux and advection forms. The time integration is performed by the forward scheme for the friction term , , and combination of the leap-frog and forward scheme for the other terms. The calculation method of pressure term are shown in $BBh(B2.3$B@a(B.

 (1) (2) (3)

 (4)

 (5) (6) (7)

 (8) (9) (10) (11) (12) (13)

 (14) (15) (16) (17) (18) (19)

## 2.2 Thermodynamics equation

The advection terms of equation (5) in Part I are evaluated by forth order centered scheme. In time integration, the forward scheme is adapted for the friction terms , , the radiative heating term and the dissipative heating term . The calculation method of radiative heating term is shown in $BBh(B5$B@a(B.

 (20)

 (21)

 (22)

 (23)

 (24)

 (25)

## 2.3 Diagonostic equation of pressure function

The diagnostic equation of nodimensional pressure function is solved by using the dimension reduction method. Before making the finite difference equation, equation (9) show in Part I is transformed as follows.

The finite difference form of the pressure equation can be written in matrix form as follows.
 (26)

where are matrixes whose elements are finite difference form of following terms.

(26) can be rewritten by using the eigenvalue matrix and the eigenvector matrix of .

where and . The final form of matrix equation is as follows.
 (27)

In calculating elements of matrix , the vertical derivative in

are evaluated by the second and forth order centered schemes because the space differencing in the continuity equation is evaluated by the forth order centered scheme while that in the pressure gradient term is evaluated by the second order centered scheme. Therefore is represented as a band matrix whose elements are given as follows.
 (28) (29) (30) (31) (32)

The boundary conditions are at the lower and upper boundary.

The horizontal differencing is evaluated by using the Fourier expansion.

 (33) (34) (35)

 (36)

## 2.4 Basic state equations

The basic state pressure () and density () are calculated by the hydrostatic equation and the equation of state given the basic state temperature .

 (37) (38)

, are calculated by using , as follows.
 (39) (40)

: 3 Turbulent parameterization : Two dimensional anelastic model : 1 Outline of discretizatoin
Odaka Masatsugu $BJ?@.(B19$BG/(B4$B7n(B26$BF|(B