% 表題  中島健介対流モデル解説  基本場の計算
%
% 履歴  1998/07/03  小高正嗣
%       1999/06/14  小高正嗣
%
\subsection{基本場の計算}

サブルーチン {\tt VBASIC} を呼び基本場の計算を行なう. 計算結果は {\tt
@T1\{CNUM(IRUN+1)\}.VB.DATA} に出力される.

\subsubsection{VBASIC} 

\begin{enumerate}
  \item 機能
    \begin{quote}
      基本場の計算を行なう.
    \end{quote}
  \item パラメータ
    \begin{quote}
      \begin{tabular}{lllp{10cm}}
        {\tt ZGRID} & (実数) & {\tt ( -2:NZB )} & 鉛直方向格子点座標 \\
        {\tt PSFC}  & (実数) &                  & 地表気圧 \\
        {\tt PRES0} & (実数) & {\tt ( -2:NZB )} & 基本場の圧力 \\
        {\tt DENS0} & (実数) & {\tt ( -2:NZB )} & 基本場の密度 \\
        {\tt TEMP0} & (実数) & {\tt ( -2:NZB )} & 基本場の温度 \\
        {\tt TPOT0} & (実数) & {\tt ( -2:NZB )} & 基本場の温位 \\
        {\tt PPAI0} & (実数) & {\tt ( -2:NZB )} & 基本場の無次元圧力関数 \\
        {\tt QVAP0} & (実数) & {\tt ( -2:NZB )} & 基本場の水蒸気混合比. \\
        {\tt RVAP0} & (実数) & {\tt ( -2:NZB )} &  \\
        {\tt GRAV}  & (実数) &                  & 重力加速度 \\
        {\tt RDRY}  & (実数) &                  & 乾燥空気気体定数 \\
        {\tt CPDRY} & (実数) &                  & 乾燥空気定圧比熱 \\
        {\tt PSD}   & (実数) &                  & 基準気圧 \\
        {\tt RMVAP} & (実数) &                  & 水蒸気分子量 \\
        {\tt NZ}    & (整数) &                  & 鉛直方向格子点数 \\
        {\tt NZB}   & (整数) &                  & 鉛直方向次元数 
      \end{tabular}
    \end{quote}
  \item 解説
    \begin{enumerate}
    \item 基本場の温度分布 $T_{0}$ を 1m 毎に与え, 内部変数 {\tt ZPROF}
      に用意する.
    \item 静水圧平衡の式
        \[
          d \ln p_{0} = - \frac{g}{RT_{0}}dz
        \]
        及び状態方程式 $p_{0}=\rho _{0}RT_{0}$ から基本場の圧力及び密
        度分布 $p_{0}, \rho _{0}$ を計算し, それぞれ内部変数 {\tt
        pprof, dprof} に代入する.
      \item 1m 毎に与えた {\tt pprof}, {\tt dprof} を補間
        して高度 {\tt ZGRID} 上の値 {\tt PRES0}, {\tt DENS0} を求める. 
      \item {\tt PRES0}, {\tt DENS0} から {\tt TEMP0}, {\tt PPAI0},
        {\tt TPOT0} を求める.
      \item 上下端における境界条件は壁で対称とする. 例えば,
        \begin{quote}
          {\tt DENS0 ( -1 ) = DENS0 ( 0 ) } \\
          {\tt DENS0 ( NZ ) = DENS0 ( NZ-1 ) }
        \end{quote}
        とする.
      \end{enumerate}
\end{enumerate}