%表題 bsnsq_2d 力学ライブラリ (渦度から流線関数の計算)
%
%履歴 
\Drireki{95/08/29 竹広真一}
%
%
\section{サブルーチンの説明 : 
         渦度からの流線関数の計算(ポアソン方程式の解)}
\label{dztpsi}

\subsection{DZ2PSI}
\label{dz2psi}
\begin{enumerate}
  \item 機能 
  \begin{quote}
    ポアソン方程式 
    ${\displaystyle \left( \DP[2]{}{x} + \DP[2]{}{z} \right) \psi = \zeta}$
    を解いて, 渦度 $\zeta$ から流線関数 $\psi$ を求める. 
  \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
          CALL DZ2PSI
         O      ( PSI   , 
         I        ZETA  , 
         I        PSIBTM, PSITOP, 
         I        DX    , DZ     )
      \end{verbatim}
  \item パラメーターの説明 
  \begin{quote}
    \begin{tabular}{llp{11cm}}
      {\tt PSI(0:NX+1,0:NZ+1)}  & {\tt [R]} & 流線関数\\
      {\tt ZETA(0:NX+1,0:NZ+1)} & {\tt [R]} & 渦度\\
      {\tt PSIBTM}          & {\tt [R]} & 下端境界の流線関数値\\
      {\tt PSITOP}          & {\tt [R]} &  下端境界の流線関数値\\
      {\tt DX}              & {\tt [R]} & X 方向の格子点間隔\\
      {\tt DZ}              & {\tt [R]} & Z 方向の格子点間隔
    \end{tabular}
  \end{quote}
  \item 備考 
  \begin{enumerate}
    \item 流線関数 $\psi$ に関する境界条件は. 
          水平方向(X)は周期的境界条件, 
          鉛直方向(Z)は {\tt PSI=const.} である. 
  \end{enumerate}
\end{enumerate}

\subsection{DZ2PST}
\label{dztpst}
\begin{enumerate}
  \item 機能 
  \begin{quote}
    ポアソン方程式 
    ${\displaystyle \left( \DP[2]{}{x} + \DP[2]{}{z} \right) \psi = \zeta}$
    の係数行列の計算を行う.
  \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
          CALL DZTPST
         O       (  ALU  , 
         I          DX   , DZ  )
      \end{verbatim}
  \item パラメーターの説明 
  \begin{quote}
    \begin{tabular}{llp{11cm}}
      {\tt ALU(NX,NZ,-1:1)} & {\tt [R]} & LU 分解された係数行列(3重対角)\\
      {\tt DX}              & {\tt [R]} & X 方向の格子点間隔\\
      {\tt DZ}              & {\tt [R]} & Z 方向の格子点間隔\\
    \end{tabular}
  \end{quote}
  \item 備考 
  \begin{enumerate}
    \item 係数行列の設定時には境界条件が組みこまれている. 
          水平方向(X)は周期的境界条件, 
          鉛直方向(Z)は {\tt PSI=const.} である. 
  \end{enumerate}
\end{enumerate}

%\newpage