%表題   物理過程 地表面フラックス
%
%履歴   93/09/01 保坂征宏
%

\section{サブルーチンの説明 : 地表面フラックス}


\subsection{SFCUVT}
\label{sfcuvt}
\begin{enumerate}
  \item 機能 
    \begin{quote}
       地表面・大気下端 ( 1/2 レベル ) の
       上向きフラックスバルク法で評価し, 
       各フラックス変数に足し込む. 
       また, 陰的解法で用いられる行列要素も計算する. 
   \end{quote} 
   
  \item 呼び出し方法
   \begin{verbatim}
      SUBROUTINE SFCUVT           !"  地表フラックス
     M         ( UFLUX , VFLUX , TFLUX , QFLUX ,
     M           UMTXS , TMTXS , QMTXS ,
     I           GDU   , GDV   , GDT   , GDTM  , GDTG  , GDQ   ,
     I           GDP   , GDPM  , GDZ   ,
     I           GRBET , GRZ0  , IDSFC                          )
   \end{verbatim}

  \item パラメータの説明
    \begin{verbatim}
*   [MODIFY] 
      REAL       UFLUX ( IDIM*JDIM, KMAX+1 ) !" Ｕのフラックス
      REAL       VFLUX ( IDIM*JDIM, KMAX+1 ) !" Ｖのフラックス
      REAL       TFLUX ( IDIM*JDIM, KMAX+1 ) !" Ｔのフラックス
      REAL       QFLUX ( IDIM*JDIM, KMAX+1 ) !" ｑのフラックス
*
      REAL       UMTXS ( IDIM*JDIM )             !" ｕ陰解行列：地表
      REAL       TMTXS ( IDIM*JDIM , 0:1, -1:1 ) !" Ｔ陰解行列：地表
      REAL       QMTXS ( IDIM*JDIM , 0:1, -1:1 ) !" ｑ陰解行列：地表
*
*   [INPUT] 
      REAL       GDU   ( IDIM*JDIM, KMAX )   !" 西風  ｕ
      REAL       GDV   ( IDIM*JDIM, KMAX )   !" 南風  ｖ
      REAL       GDT   ( IDIM*JDIM, KMAX )   !" 温度  Ｔ
      REAL       GDTM  ( IDIM*JDIM, KMAX )   !" 温度（半整数）
      REAL       GDTG  ( IDIM*JDIM       )   !" 地表面温度
      REAL       GDQ   ( IDIM*JDIM, KMAX )   !" 比湿  ｑ
      REAL       GDP   ( IDIM*JDIM, KMAX )   !" 気圧
      REAL       GDPM  ( IDIM*JDIM, KMAX )   !" 気圧（半整数）
      REAL       GDZ   ( IDIM*JDIM, KMAX )   !" 高度
*
      REAL       GRBET ( IDIM*JDIM )         !" 地表湿潤度β
      REAL       GRZ0  ( IDIM*JDIM )         !" 地表粗度長
      INTEGER    IDSFC ( IDIM*JDIM )         !" 地表状態
    \end{verbatim}

 \item 備考
   \begin{enumerate}
     \item 現在はフラックスの評価方法として
           バルク法のみが用意されている. 
     \item \verb!*MTXS! は陰的解法で用いられる行列要素である. 
     \item バルク係数の計算は \verb!BLKCOF! でなされる. 
   \end{enumerate}
\end{enumerate}
