A 4-level baroclinic numerical model is designed by using the vorticity equation and Omega equation. Block-Cyclic-Reduction method is applied to the solution of the Helmholtz defferential equation, which is proved to be better than the Relaxation method from the composite viewpoint of accuracy, stability and economy. It was investigated whether the model explains the physical process influenced by voricity and temperature advection. It was also examined if the model atmosphere describes the general circulation. This examination is similar to Phillips(1956). The result of this numerical experiment shows that the model explains qualitatively the Quasi-Geostrophic theory for the development of Baroclinic wave, as throughly described in Holton(1972).