% 表題   AGCM5 第3部 コード解説    力学過程
%
% 履歴
\Drireki{97/03/06  堀之内武}
%

\Dchapterhead

\section{力学過程の概要}

本章では力学過程のコードの構造と計算手順を説明する. ここでいう力学過程と
はスペクトルで計算される部分のことを指し, 移流項や水平拡散項, 圧力勾配項
などである.  あるいは物理過程以外と言っても良いかもしれない. スペクトル
とは球面調和関数展開の係数のことであり, スペクトルへの変換とは球関数展開
に他ならない. 

実際の時間積分は「力学過程」の中で行われる. 物理過程で計算された時間変化
項は力学過程に引き渡され, 「力学過程」の時間変化項に加えられ, セミインプ
リシット法による時間積分に組み込まれる. 但し, 降水過程などの調節過程は時
間積分後に適用される. 

\subsection{セミインプリシット時間積分の概要}\label{ss:3bu_semiimp_abs}

セミインプリシット法についてはドキュメント第2部の 2.4 節に解説があるので
参照のこと. 時刻 $t-\Delta t$ から $t+\Delta t$ へのセミインプリシットの
計算手順は以下のようにまとめられる. 

\begin{enumerate}
 \item 発散項 $\Dvect{D}$ のみに関する単一のセミインプリシット方程式,
   \begin{equation}
     \underline{M}\, \overline{\Dvect{D}}^{t} = \Dvect{f}
                                         \label{eq:3bu_Dsemiimp_schematic}
   \end{equation}
   を解く. ここで
   \( \displaystyle \overline{\Dvect{D}}^{t} \equiv 
      \frac{1}{2} (\Dvect{D}^{t+\Delta t} + \Dvect{D}^{t-\Delta t}) \)
   である. 

 \item $\overline{\Dvect{D}}^{t}$ と時刻 $t$ における時間変化の非重力波
       (NG)項から各変数の $t+\Delta t$ 
       での値を求める.

\end{enumerate}

以上の手順の説明は \ref{s:semiimp詳細} 節で行う.
