\section{計算手法}

\subsection{計算手順}

地表気圧, 各種物理定数(比熱, 重力加速度, 気体定数)はパラメータとして与
える. 計算は以下の手順で行なう.

\begin{itemize}
  \item [\bf 1.] 適当な温度(温位)分布を初期条件として与える. 

  \item [\bf 2.] 状態方程式と静水圧平衡の式\Deqref{静水圧の式}を組合せ,  
	     密度と圧力の鉛直分布を計算する.

  \item [\bf 3.] \Deqref{基礎方程式}式を放射加熱/冷却の寄与だけを用いて, 
	     温度(温位)分布の仮値を求める.
	
  \item [\bf 4.] 温度(温位)分布の仮値から鉛直拡散係数を計算する.

  \item [\bf 5.] \Deqref{基礎方程式}式の拡散項の計算を行なう.

  \item [\bf 6.] 2. に戻る\Deqref{基礎方程式}式の拡散項の計算を行なう.
\end{itemize}

\subsection{設定}

\begin{itemize}
  \item 大気は CO${}_{2}$ のみで構成されるとする. 物性定数はその値にする.
  \item 計算領域は高度 50 km までとる.
        鉛直格子間隔は 500m , 時間格子間隔は 1 分にとる.
  \item 計算は前日との温度差が 0.2 K 程度になったところまで行なう.
	計算モデル時間はおよそ 150 日前後である.
\end{itemize}

