[地球流体電脳倶楽部] [dcmodel | dcmodel-tools] [dcmodel-thum.rb]

新旧ソース比較

計算の概要

arare5-takahashisan と arare5-20130125 の計算結果を比較する. 両方とも地表面フラックスは純正バルク法で計算する, すなわち強制は入れていない. 旧ソースは熱膨張項なし, 新ソースは熱膨張項なし / あり両方.

実験設定

基本方針は, Nakajima and Matsuno (1988) に近い条件での計算.

• モデル: deepconv / arare5
  • 旧ソース(arare5-takahashisan)
  • 新ソース(arare5-20130125)
  • 移流
    • 準圧縮方程式系
  • 乱流混合
    • 1.5 次クロージャ―
    • バルク法による下部境界フラックス
      • バルク係数 1.5d-3
      • 風速のかさ上げ 3 m/s
  • 温かい雨のパラメタリゼーション
  • 摩擦層
    • z > 25 km に適用
• 領域: 2 次元
  • 水平 256 km
  • 鉛直 48 km
• 格子間隔
  • 水平 1000 m
  • 鉛直 300 m
• 下部境界
  • 温度 302 K (固定)
• 大気組成
  • 乾燥成分 : N2 (空気ではない)
  • 凝結成分 : H2O
• 熱強制
  • -1 K/day for z<= 10 km
  • 0 K/day for z>= 15 km
  • 10 km < z < 15 km では線形でつなぐ
• 温かい雨のパラメタリゼーション
  • 雲水から雨水への変換の閾値 q_{cr} = 1e-3
• 積分時間
  • (ひとまず) 8 日
• タイムステップ
  • 旧ソース
    • 長いタイムステップ: 2 sec
    • 短いタイムステップ: 1 sec
  • 新ソース
    • 長いタイムステップ: 2 sec
    • 短いタイムステップ: 0.5 sec
• NAMELIST ファイルへのリンク
  • 旧ソース用
    • 初期値生成用
    • メイン計算用
  • 新ソース用(熱膨張項なし)
    • 初期値生成用
    • メイン計算用
  • 新ソース用(熱膨張項あり)
    • 初期値生成用
    • メイン計算用

基本場

以下のように与える.

• 乾燥断熱減率に沿った温度分布を求める
• 静水圧平衡から圧力分布を求める
• 水蒸気混合比は湿度の分布と飽和水蒸気圧から求める

計算結果

左: 旧ソース, 熱膨張項なし (arare5-takahashisan) 中: 新ソース, 熱膨張項なし (arare5-20130125) 右: 新ソース, 熱膨張項あり (arare5-20130125)

dcl_Exner_t_atSfc_1-2-2_z48km.gif	dcl_Exner_t_atSfc_1-2_z48km.gif
dcl_Exner_t_atSfc_1-4-2_z48km.gif	dcl_Exner_tz_1-2_z48km.gif
dcl_Exner_tz_1-4-2_z48km.gif	dcl_H2O-gAll_tz_1-2-2_z48km.gif
dcl_H2O-gAll_tz_1-2_z48km.gif	dcl_H2O-gAll_tz_1-4-2_z48km.gif
dcl_H2O-gAll_tz_z=0-10km_1-2-2_z48km.gif	dcl_H2O-gAll_tz_z=0-10km_1-2_z48km.gif
dcl_H2O-gAll_tz_z=0-10km_1-4-2_z48km.gif	dcl_H2O-g_Cond_tz_1-2-2_z48km.gif
dcl_H2O-g_Cond_tz_1-2_z48km.gif	dcl_H2O-g_Cond_tz_1-4-2_z48km.gif
dcl_H2O-l-CloudAll_tz_1-2-2_z48km.gif	dcl_H2O-l-CloudAll_tz_1-2_z48km.gif
dcl_H2O-l-CloudAll_tz_1-4-2_z48km.gif	dcl_H2O-l-CloudAll_tz_z=0-20km_1-2-2_z48km.gif
dcl_H2O-l-CloudAll_tz_z=0-20km_1-2_z48km.gif	dcl_H2O-l-CloudAll_tz_z=0-20km_1-4-2_z48km.gif
dcl_H2O-l-RainAll_tz_1-2-2_z48km.gif	dcl_H2O-l-RainAll_tz_1-2_z48km.gif
dcl_H2O-l-RainAll_tz_1-4-2_z48km.gif	dcl_H2O-l-RainAll_tz_z=0-20km_1-2-2_z48km.gif
dcl_H2O-l-RainAll_tz_z=0-20km_1-2_z48km.gif	dcl_H2O-l-RainAll_tz_z=0-20km_1-4-2_z48km.gif
dcl_Pressure_t_atSfc_1-2-2_z48km.gif	dcl_Pressure_t_atSfc_1-2_z48km.gif
dcl_Pressure_t_atSfc_1-4-2_z48km.gif	dcl_Temp_tz_1-2-2_z48km.gif
dcl_Temp_tz_1-2_z48km.gif	dcl_Temp_tz_1-4-2_z48km.gif
dcl_VelX_tz_1-2-2_z48km.gif	dcl_VelX_tz_1-2_z48km.gif
dcl_VelX_tz_1-4-2_z48km.gif	dcl_VelZ_tz_1-2-2_z48km.gif
dcl_VelZ_tz_1-2_z48km.gif	dcl_VelZ_tz_1-4-2_z48km.gif