%表題 bsnsq_2d 基本ライブラリ (空間差分の計算)
%
%履歴 
\Drireki{95/08/14 竹広真一}
%
%
\section{サブルーチンの説明 : 空間差分の計算}
\label{ucdiff}

\subsection{ARAJAC}
\label{arajac}
\begin{enumerate}
  \item 機能 
  \begin{quote}
    ヤコビアン
    ${\displaystyle 
      J(\psi,G) = \DP{\psi}{x}\DP{G}{z} - \DP{\psi}{z}\DP{G}{x}}$
   の差分計算を, ARAKAWA のスキームで行う. 
  \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
          CALL ARAJAC
         O       ( AJAC ,
         I         PSI  , G   ,
         C         DX   , DZ  ,
         I         ISX  , IFX , ISZ , IFZ ,
         D         NX   , NZ    )
      \end{verbatim}
  \item パラメーターの説明 
  \begin{quote}
    \begin{tabular}{llp{11cm}}
      {\tt AJAC(NX,NZ)} & {\tt [R]} & ARAKAWA のヤコビアン {\tt J(PSI,G)}\\
      {\tt PSI(NX,NZ)}  & {\tt [R]} & 流線関数\\
      {\tt G(NX,NZ)}    & {\tt [R]} & 物理量の2次元場\\
      {\tt DX}   & {\tt [R]} & X 方向の格子点間隔\\
      {\tt DZ}   & {\tt [R]} & Z 方向の格子点間隔\\
      {\tt ISX}  & {\tt [I]} & X 方向にヤコビアンを計算し始める点\\
      {\tt IFX}  & {\tt [I]} & X 方向にヤコビアンを計算し終える点\\
      {\tt ISZ}  & {\tt [I]} & Z 方向にヤコビアンを計算し始める点\\
      {\tt IFZ}  & {\tt [I]} & Z 方向にヤコビアンを計算し終える点\\
      {\tt NX}   & {\tt [I]} & X 方向の配列の次元\\
      {\tt NZ}   & {\tt [I]} & Z 方向の配列の次元
    \end{tabular}
  \end{quote}
  \item 備考 
  \begin{enumerate}
    \item 2 $\leq$ {\tt (ISX,IFX)} $\leq$ {\tt NX-1}, 
          2 $\leq$ {\tt (ISZ, IFZ)}$\leq$ {\tt NZ-1} 
          の範囲で指定すること. 
          さもなくば範囲外の配列参照が生じる. 
  \end{enumerate}
\end{enumerate}

\subsection{CLAPLA}
\label{clapla}
\begin{enumerate}
  \item 機能 
  \begin{quote}
    ラプラシアン
     ${\displaystyle 
       \nabla^2 F =  \DP[2]{F}{x} + \DP[2]{F}{z} }$
    の差分計算を行う. 
  \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
          CALL CLAPLA
         O       ( ALAP,
         I         F   ,
         C         DX  , DZ  , 
         I         ISX , IFX , ISZ , IFZ ,
         D         NX  , NZ    )
      \end{verbatim}
  \item パラメーターの説明 
  \begin{quote}
    \begin{tabular}{llp{11cm}}
      {\tt ALAP(NX,NZ)}& {\tt [R]} & 物理量のラプラシアン\\
      {\tt F(NX,NZ)}   & {\tt [R]} & 物理量の2次元場\\
      {\tt DX}   & {\tt [R]} & X 方向の格子点間隔\\
      {\tt DZ}   & {\tt [R]} & Z 方向の格子点間隔\\
      {\tt ISX}  & {\tt [I]} & X 方向にラプラシアンを計算し始める点\\
      {\tt IFX}  & {\tt [I]} & X 方向にラプラシアンを計算し終える点\\
      {\tt ISZ}  & {\tt [I]} & Z 方向にラプラシアンを計算し始める点\\
      {\tt IFZ}  & {\tt [I]} & Z 方向にラプラシアンを計算し終える点\\
      {\tt NX}   & {\tt [I]} & X 方向の配列の次元\\
      {\tt NZ}   & {\tt [I]} & Z 方向の配列の次元
    \end{tabular}
  \end{quote}
  \item 備考 
  \begin{enumerate}
    \item 2 $\leq$ {\tt (ISX,IFX)} $\leq$ {\tt NX-1}, 
          2 $\leq$ {\tt (ISZ, IFZ)}$\leq$ {\tt NZ-1} 
          の範囲で指定すること. 
          さもなくば範囲外の配列参照が生じる. 
  \end{enumerate}
\end{enumerate}

%\newpage