% 表題   らくらく DCPAM5 -- 軸対称 2 次元計算を行うには
%
% 履歴 
%\Drireki{2011/09/30 高橋芳幸}
%
\chapter{軸対称 2 次元計算を行うには}


\section{はじめに}

\Dmodel は, 軸対称 2 次元計算に用いることができる. 
この章では, \Dmodel を用いた軸対称 2 次元計算の実行方法について述べる. 


\section{\Dmodel の軸対称 2 次元化の概要}

軸対称 2 次元計算のための準備と実行方法について述べる前に, \Dmodel での
軸対称 2 次元化の概要について簡単に説明しておく.

軸対称 2 次元化は, 下のふたつの方法によって実装している.
%
\begin{itemize}
\item 移流計算におけるスペクトル変換に spml の wa\_zonal\_module モジュールを用いる,
\item 経度方向の格子点数を 1 にする.
\end{itemize}
%
移流計算において用いる spml のスペクトル変換モジュールは, 
コンパイル時にプリプロセッサオプションで指定することによって選択する. 
経度方向の格子点数は, 計算実行時の設定ファイル (namelist ファイル) で指定する. 

なお, 移流計算におけるスペクトル変換以外は 3 次元計算の時に用いていたモジュールを
そのまま用いている. 
このとき, 計算される子午面は, 経度 0$^\circ$ における子午面として扱われ, 
例えば短波放射計算のように惑星上の経度 (地方時) に依存するような計算においては, 
経度 0$^\circ$ における日変化, 季節変化が計算される
\footnote{
  もちろん, 日変化, 季節変化をなくす設定をすれば話は別.
}.


\section{コンパイル}

コンパイルの基本的な方法は逐次版と同じであり, 詳細は, 
「\Dmodel インストールガイド」
(\url{http://www.gfd-dennou.org/library/dcpam/dcpam5/dcpam5_latest/INSTALL.htm}) 
を参照すること.
ただし, 下の点に注意すること. 
%
\begin{itemize}
\item 環境変数の FFLAGS に -DAXISYMMETRY または -DAXISYMMETRY\_SJPACK を指定する.
\end{itemize}


\section{軸対称 2 次元計算のための設定}

軸対称 2 次元計算の特有の設定は, 格子点数の指定である. 
下に示すように指定する. 

\subsection{格子点数の指定}

既に述べたように, 軸対称 2 次元計算は, 経度方向の格子点数を 1 にする
ことで実現している. 
経度方向の格子点数は, gridset\_nml namelist ブロックにより, 下のように指定する.
%
\begin{verbatim}
&gridset_nml
  ...
  imax   =  1                  ! 経度格子点数.
                               ! Number of grid points in longitude
  ...
/
\end{verbatim}


\section{軸対称 2 次元計算の実行}

軸対称 2 次元計算の実行方法は, 3 次元モデルの実行方法と
全く同じである. 
詳細は, 「ごくらく \Dmodel」
(\url{http://www.gfd-dennou.org/library/dcpam/dcpam5/dcpam5_latest/doc/tutorial/gokuraku/}) 
を参照すること.