: 4 積雲パラメタリゼーション
 : DCAPM4 第1部 数理モデル化
 : 2 座標系
  
1
11
ここでは支配方程式を記し, 特に力学過程2について詳細を記す.
 
まず支配方程式の中で力学過程と認知される部分・項を示す. 
ついで各々の詳細を述べる. 
  
2
22
2
22
ここでは支配方程式を順に示す. 
この方程式の詳細に関しては, Haltiner and Williams (1980) もしくは
第A章 を参照せよ.
  | 
(1) | 
 
  | 
(2) | 
 
  | 
(3) | 
 
  | 
(4) | 
 
ここで,
  | 
  | 
  | 
(7) | 
  | 
  | 
  | 
(8) | 
  | 
  | 
  | 
(9) | 
  | 
  | 
  | 
(10) | 
  | 
  | 
  | 
(11) | 
  | 
  | 
  | 
(12) | 
  | 
  | 
  | 
(13) | 
  | 
  | 
  | 
(14) | 
  | 
  | 
  | 
(15) | 
  | 
  | 
  | 
(16) | 
  | 
  | 
  | 
(17) | 
  | 
  | 
  | 
(18) | 
  | 
  | 
  | 
(19) | 
  | 
  | 
  | 
(20) | 
  | 
  | 
  | 
(21) | 
  | 
  | 
  | 
  | 
|   | 
  | 
  | 
(22) | 
  | 
  | 
![$\displaystyle \frac{1}{a^{2}(1-\mu^{2})}
\frac{\partial^{2} }{\partial \lambda^...
...ial }{\partial \mu}
\left[ (1-\mu^{2})
\frac{\partial }{\partial \mu} \right] .$](img65.png)  | 
(23) | 
 
ただし, 
は水平拡散項であり, 3.3 で説明される. 
は小規模運動過程による力である. 
 は放射, 凝結, 小規模運動過程等による加熱・温度変化, 
は凝結, 小規模運動過程等による水蒸気ソース項, 
 は摩擦熱である. 
鉛直流に関する境界条件は
  | 
(24) | 
 
である. よって(3.1) から,
地表気圧の時間変化式と
系での鉛直速度
を求める診断式
  | 
(25) | 
 
  | 
(26) | 
 
が導かれる.
ただし熱的境界条件については 6 章において記述する. 
2
22
3
水平拡散項は, 次のように
の形で計算されるのが普通である. 
![\begin{displaymath}
{\cal D}(\zeta) = - K_{HD}
\left[ (-1)^{N_D/2} \nabla^{N_D}
- \left( \frac{2}{a^2} \right)^{N_D/2}
\right]
\zeta ,
\end{displaymath}](img77.png)  | 
(27) | 
 
![\begin{displaymath}
{\cal D}(D) = - K_{HD}
\left[ (-1)^{N_D/2} \nabla^{N_D}
- \left( \frac{2}{a^2} \right)^{N_D/2}
\right]
D ,
\end{displaymath}](img78.png)  | 
(28) | 
 
  | 
(29) | 
 
  | 
(30) | 
 
この水平拡散項は計算の安定化のための意味合いが強い.
小さなスケールに選択的な水平拡散を表すため,
としては, 4
16を用いる.
水平拡散を波数に依存しない一様な値にすることもできる.
詳細省略. 
- Haltiner, G.J. and  Williams, R.T., 
        1980:
        Numerical Prediction and Dynamic Meteorology (2nd ed.).
        John Wiley & Sons, 477pp.
 
- ... 特に力学過程2
 
-  力学過程という単語が適切かどうかは不明である.  実態は, 
  モデルにおいて各格子点で計算されない部分を指す.  
 
- ...P3
 
-  (2005/4/4 石渡) 力学過程という節が昔存在していたが, 必要か???
 
 
 
   
 : 4 積雲パラメタリゼーション
 : DCAPM4 第1部 数理モデル化
 : 2 座標系
Yasuhiro MORIKAWA
平成19年8月3日