% 表題   Nakajima et al.,1992 test program  コード解説 -- グリッドと物理量の配置について
%
% 履歴 
%\Drireki{2012/11/09  大西 将徳}
%\Drireki{2012/11/09  大西 将徳}
%
\section{グリッドと物理量の配置について}


\Dmodel は\ref{sec:grid}節で述べたように、温度はグリッド上に、Flux ConvergenceはLayerに対して計算される。

\begin{verbatim}

k = 1    -----  大気上端
           :
k = i    -----  z_Press(i), z_Temp(i)

                FluxConv(i)

k = i+1  -----
           :
k = kmax -----  地表面
\end{verbatim} %\vspace{-2ex}

そのためFlux Convergenceの大きさに従って温度を上昇させる場合に、layerに対して与えられた温度上昇率から、
グリッド上で定義された温度を上昇させている。具体的には、温度が定義されたグリッドの上下の層の温度上昇分を
半分ずつ足すことで温度を上昇させている。

\begin{verbatim}

k = 1    -----  大気上端
           :
k = i-1  -----
              
                  z_FluxConv(i-1), delTemp(i-1)

k = i    -----  z_Temp(i)

                  z_FluxConv(i), delTemp(i)

k = i+1  -----
           :
k = kmax -----  地表面


  z_Temp(i) = z_Temp(i) + (delTemp(i-1) + delTemp(i))/2
  (delTemp(i): z_FluxConv(i)から計算された温度の上昇分)

  ただし、大気上端では以下のように与える。
  z_Temp(1) = z_Temp(1) + delTemp(1)
\end{verbatim} %\vspace{-2ex}

温度をグリッドで与え、温度上昇をlayerに対して計算するという、'ずれ'のために、このまま温度を上昇させると、
上下のグリッドで温度の大小が連続する鉛直分布が実現してしまうことが、明らかとなった。
このため、本プログラムではこのような温度分布が生じないように、プログラム上の制約を課して温度上昇をさせている。
それは、'隣り合うグリッドのFlux Convergence が正から負になる場合(地上方向から上空方向へ)
は温度を変化させない'というものである。これにより、温度の大小が連続するこをは抑えられた。
ただし、対流圏界面付近ではそのような現象が生じる場合がある。