%表題   物理過程 湿潤対流調節
%
%履歴   93/10/12 保坂征宏
%
\Drireki{93/10/12 保坂征宏}
%
\section{サブルーチンの説明 : 湿潤対流調節}

\subsection{CUMADJ ( adustment )}
\label{cumadj(adjustment)}
\begin{enumerate}
  \item 機能 
    \begin{quote}
       湿潤対流調節の計算. 
       温度減率が湿潤断熱減率よりも大きく, 
       飽和または過飽和にあるときに
       湿潤安定度が中立に, 
       また温度の鉛直積分は潜熱解放分だけ上昇するように, 
       飽和状態になるように, 
       降水を起こし
       温度・比湿を変化させる.
    \end{quote}

  \item 呼び出し方法
  
    \begin{verbatim}
      SUBROUTINE CUMADJ       !" 湿潤対流調節
     M         ( GDT   , GDQ   ,
     M           GQPR  , GQCNV , GSCNV ,
     I           GDP   , GDPM  , GDZ   , DELTA  )
    \end{verbatim}
    
  \item パラメタの説明
  
    \begin{verbatim}
*   [MODIFY] 
      REAL       GDT   ( IDIM*JDIM, KMAX )   !" 温度  Ｔ
      REAL       GDQ   ( IDIM*JDIM, KMAX )   !" 比湿  ｑ
      REAL       GQPR  ( IDIM*JDIM       )   !" 降水量
      REAL       GQCNV ( IDIM*JDIM, KMAX )   !" 凝結調節加熱率
      REAL       GSCNV ( IDIM*JDIM, KMAX )   !" 凝結調節比湿変化
*
*   [INPUT] 
      REAL       GDP    ( IDIM*JDIM, KMAX )  !" 気圧  Ｐ
      REAL       GDPM   ( IDIM*JDIM, KMAX+1 )!" 気圧（半整数）
      REAL       GDZ    ( IDIM*JDIM, KMAX )  !" 高度
      REAL       DELTA                       !" 2Δt
*
*ネームリスト
      NAMELIST   NMMADJ
      REAL       EPST ( ITRMAX )             !" 不安定の許容誤差
      REAL       RCRIT                       !" 臨界相対湿度
     \end{verbatim}
     
        
  \item 備考
     \begin{enumerate}
       \item 共通ワーク領域 \verb!COMWRK! を
             \verb!2*IDIM*JDIM!+\verb!6*IDIM*JDIM*KMAX! 使用. 
       \item 不安定の評価ために, 
             \begin{eqnarray}
               \DP{T}{z} - \frac{RT}{C_p p} \DP{p}{z} 
               + \frac{L}{C_p} \DP{q^*}{z} > 0
             \end{eqnarray}
             を離散化したものを用いている. 
             ここで $T,p$ の中間レベルでの値が必要になる. 
             この際に, 中間レベルは本来力学過程でのそれと
             一致している必要はないことに注意しよう. 
             さて, ここで
             $p$ については力学過程で用いたのと同じ
             内挿公式を用いたのに対し, 
             $T$ については算術平均の
             \begin{eqnarray}
                   T_{k-1/2} = \frac{T_{k-1}+   T_k    }{2}
             \end{eqnarray}
             を用いた. 
             これは乾燥対流調節での内挿公式とも食い違っており, 
             今後の検討を要する. 
       \item 現状では質量の重みをかけた温度の鉛直積分を
             保存するように調節しているが, 
             力学過程などとの関係上, 
             温位の鉛直積分を保存させるように
             すべきである. 
             熱のソースがあることになっていて, これはバグである. 
       \item 下から順々に下層を冷やし上層を温めるという調節を
             繰り返し行なうというスキームなので,  
             原理的には正確に調節させるのは不可能である. 
             実際, 鉛直各レベル間で最高 0.4 K の
             不安定の許容誤差を許しており, 
             これが鉛直成層を変えている可能性も少なくない. 
             うまい打開策は検討課題である. 
       \item 単位時間あたり降水量, 凝結率, 比湿変化率を求めるときに 
             \verb!DELTA! を用いているのはいいのか？
             モデル時間刻を用いるべきであるような気がとってもするが. 
     \end{enumerate}

\end{enumerate}
