% 表題  中島健介対流モデル解説  乱流拡散係数の計算 (TSTEP)
%
% 履歴  1998/07/07  小高正嗣
%       1999/06/14  小高正嗣
%
\subsection{乱流拡散係数の計算}

乱流拡散係数を計算する. ここでは拡散係数がある閾値よりも大きくならない
ようにリミッターをかける. 

\subsubsection{CDCOEF}

\begin{enumerate}
  \item 機能
    \begin{quote}
      乱流拡散係数を計算する.
    \end{quote}
  \item パラメータ
    \begin{quote}
      \begin{tabular}{lllp{10cm}}
        {\tt CDTURB}& (実数) &{\tt ( -2:NXB, -2:NZB, NROT )} & 乱流拡散係数 \\
        {\tt IC}  & (整数) & & 時間方向次元 \\
        {\tt CD1} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 乱流拡散係数(
        $(i+\frac{1}{2}), j+\frac{1}{2})$ 格子上) \\
        {\tt CD2} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 乱流拡散係数(制限値) \\
        {\tt NX  } & (整数) & & 水平方向格子点数 \\
        {\tt NZ  } & (整数) & & 鉛直方向格子点数 \\
        {\tt NXB } & (整数) & & 水平方向次元数 {\tt (=nx+2)} \\
        {\tt NZB } & (整数) & & 鉛直方向次元数 {\tt (=nz+2)} \\
        {\tt NROT} & (整数) & & 時間方向次元数 
      \end{tabular}
    \end{quote}
  \item 解説
    \begin{enumerate}
      \item 引数 {\tt IC} には {\tt ITC1} が与えられる.
      \item 内部で拡散係数の上限値({\tt CDMAX}) と 下限値({\tt CDMIN})
        を定義してから, 一ステップ前の計算により得られた拡散係数{\tt
        CDTURB}に {\tt CDMAX} と {\tt CDMIN} によるフィルターをかけた
        {\tt CD2}を求める. 次に {\tt CD2} の 4 点平均により
        ($(i+\frac{1}{2}), j+\frac{1}{2})$) 格子上の拡散係数 {\tt CD1}
        を求める.
      \end{enumerate}
\end{enumerate}


