%表題   物理過程 メイン
%
%履歴   93/09/01 保坂征宏
%

\section{サブルーチンの説明 : 物理過程メインプログラム}

\subsection{PHYSCS}
\label{physcs}
\begin{enumerate}
  \item 機能 
    \begin{quote}
       物理過程の計算. 地表面過程を含む鉛直拡散, 放射. 
    \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
      SUBROUTINE PHYSCS      !" 物理過程メイン
     O         ( GTU   , GTV   , GTT   , GTQ   , 
     M           GDTG  , 
     I           GDU   , GDV   , GDT   , GDPS  , GDQ   ,
     I           GDVOR , GDDIV , 
     I           IT    , IDATE , DELT  , DELTM , 
     C           ALON  , ALAT  ,
     C           SIG   , SIGM  , DSIG  , DSIGM          )
      \end{verbatim}
      
  \item パラメーターの説明 
      \begin{verbatim}
*
*   [OUTPUT] 
      REAL       GTU   ( IDIM, JDIM, KMAX )  !" 東西運動量変化 ＵＡ
      REAL       GTV   ( IDIM, JDIM, KMAX )  !" 南北運動量変化 ＶＡ
      REAL       GTT   ( IDIM, JDIM, KMAX )  !" 温度時間変化項 Ｈ
      REAL       GTQ   ( IDIM, JDIM, KMAX )  !" 比湿時間変化項 Ｒ
*
*   [MODIFY]
      REAL       GDTG  ( IDIM, JDIM )        !" 地表面温度
*
*   [INPUT] 
      REAL       GDU   ( IDIM, JDIM, KMAX )  !" 西風  ｕ
      REAL       GDV   ( IDIM, JDIM, KMAX )  !" 南風  ｖ
      REAL       GDT   ( IDIM, JDIM, KMAX )  !" 温度  Ｔ
      REAL       GDPS  ( IDIM, JDIM       )  !" 地表気圧 Ｐｓ
      REAL       GDQ   ( IDIM, JDIM, KMAX )  !" 比湿  ｑ
      REAL       GDVOR ( IDIM, JDIM, KMAX )  !" 渦度  ζ
      REAL       GDDIV ( IDIM, JDIM, KMAX )  !" 発散  Ｄ
*
      INTEGER    IT                          !" 通し時間
      INTEGER    IDATE ( * )                 !" 時刻：年月日時分秒
      REAL       DELTM                       !" モデル時間刻みΔt
      REAL       DELT                        !" 物理過程時間刻みΔt
*
      REAL       ALON  ( IDIM )              !" 経度
      REAL       ALAT  ( JDIM )              !" 緯度
      REAL       SIG   ( KMAX )              !" σレベル(整数)
      REAL       SIGM  (KMAX+1)              !" σレベル(半整数)
      REAL       DSIG  ( KMAX )              !" Δσ
      REAL       DSIGM (KMAX+1)              !" Δσ(半整数)
*
*   [INTERNAL SAVE] 
      REAL       RFLUX  ( IDIM, JDIM, KMAX+1 )   !" 長波フラックス
      REAL       SFLUX  ( IDIM, JDIM, KMAX+1 )   !" 短波フラックス
      REAL       DRFLUX ( IDIM, JDIM, KMAX+1, 0:1 ) !" 長波温度変化
       \end{verbatim}

  \item 備考
  \begin{enumerate}
    \item 放射は数ステップ毎にしか計算しないことを
          前提とした構成をしている. 
          放射フラックスを内部的に SAVE しているのもそのためである.  
    \item 鉛直拡散は陰解法を用いることを前提としている. 
          ここで用いる時間刻はモデル時間刻と異なるように
          とることができる. ( デフォルトはモデル時間刻の 2 倍. )
  \end{enumerate}
\end{enumerate}

\subsection{PADJST}
\label{padjst}
\begin{enumerate}
  \item 機能 
    \begin{quote}
       調節過程の計算. 雲を含む対流過程. 
    \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
      ENTRY      PADJST       !" 物理調節
     M         ( GDU   , GDV   , GDT   , GDPS  , GDQ   ,
     I           GDVOR , GDDIV ,
     I           DELT  ,
     C           SIG   , SIGM  , DSIG  , DSIGM          )
      \end{verbatim}
      
  \item パラメーターの説明 
      \begin{verbatim}
*
*   [MODIFY]
      REAL       GDU   ( IDIM, JDIM, KMAX )  !" 西風  ｕ
      REAL       GDV   ( IDIM, JDIM, KMAX )  !" 南風  ｖ
      REAL       GDT   ( IDIM, JDIM, KMAX )  !" 温度  Ｔ
      REAL       GDPS  ( IDIM, JDIM       )  !" 地表気圧 Ｐｓ
      REAL       GDQ   ( IDIM, JDIM, KMAX )  !" 比湿  ｑ
*
*   [INPUT]
      REAL       GDVOR ( IDIM, JDIM, KMAX )  !" 渦度  ζ
      REAL       GDDIV ( IDIM, JDIM, KMAX )  !" 発散  Ｄ
*
      REAL       DELT                        !" 物理過程時間刻みΔt
*
      REAL       SIG   ( KMAX )              !" σレベル(整数)
      REAL       SIGM  (KMAX+1)              !" σレベル(半整数)
      REAL       DSIG  ( KMAX )              !" Δσ
      REAL       DSIGM (KMAX+1)              !" Δσ(半整数)
       \end{verbatim}
      
  \item 備考
    \begin{quote}
      なし. 
    \end{quote}
\end{enumerate}


\subsection{PHSTRG}
\label{phstrg}
\begin{enumerate}
  \item 機能 
    \begin{quote}
       物理過程時間平均出力の登録. 
    \end{quote}
  \item 呼び出し方法 
      \begin{verbatim}
      SUBROUTINE PHSTRG    !" 物理過程時間平均出力の登録
      \end{verbatim}
  
  \item パラメータの説明
    \begin{quote}
      なし. 
    \end{quote}

  \item 備考
    \begin{enumerate}
      \item 実際の登録は \verb!HISTRG! でを呼ぶことによって行なわれる. 
      \item 現在のところ, 以下の諸量が出力可能である. 
        \begin{verbatim}
         RAIN   precipitation  
         EVAP   latent heat flux 
         SENS   sensible heat flux 
         OLR    outgoing longwave
         SLR    surface longwave 
         OSR    outgoing shortwave 
         SSR    surface shortwave  
         TAUX   surface stress(x) 
         TAUY   surface stress(y) 
         QVDF   diffusive heating 
         SVDF   diffusive moistning 
         FVDFX  diffusive force(x) 
         FVDFY  diffusive force(y) 
         QCND   condensation heating
         SCND   condensation moistning
         QCNDC  cumulus heating 
         SCNDC  cumulus moistning
         SRMNQ  remove negative moist
         QRADL  longwave heating 
         FGRVX  gravity wave u force
         FGRVY  gravity wave v force
         TAUGX  gravity wave u flux
         TAUGY  gravity wave v flux
       \end{verbatim}
    \end{enumerate}
      
\end{enumerate} 

