- 2.1 Equation of motion
- 2.2 Thermodynamics equation
- 2.3 Diagonostic equation of pressure function
- 2.4 Basic state equations

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
第2.3節.

(4) |

(5) | |||

(6) | |||

(7) |

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
第5節.

(20) |

(21) |

(23) | |||

(25) |

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.

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) |

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) |