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

熱膨張項の見積り (水平平均, 非断熱加熱の項のみ)

計算の概要

このページでは, deepconv で -1.0 K day-1 の一様冷却計算を行ったときの, deepconv 定式化マニュアル の (A.44) における熱膨張項の水平平均値を見積もる. ただしここでは, 非断熱加熱の項のみを見る.

実験設定

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

• モデル:deepconv (基本場の圧力を用いて地表面フラックスを計算したものはこちら)
  • 移流
    • 準圧縮方程式系
  • 乱流混合
    • 1.5 次クロージャ―
    • バルク法による下部境界フラックス
      • バルク係数 1.5d-3
      • 風速のかさ上げ 3 m/s
  • 温かい雨のパラメタリゼーション
  • 摩擦層
    • z > 25 km に適用
• 領域: 2 次元
  • 水平 256 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
• 積分時間
  • 30 日

計算に用いたスクリプト

• 熱膨張項を計算して netCDF に出力するスクリプト
  • makenetCDF_t_NetsuboutyouTerm.rb
    • 同じディレクトリに thermal-moist_PTempCond_IntValue30day.nc, thermal-moist_PTempDisp_IntValue30day.nc, thermal-moist_PTempRad_IntValue30day.nc, thermal-moist_PTempTurb_IntValue30day.nc, thermal-moist_restart.nc を置いて実行すること.
• 描画スクリプト
  • draw_z_density-thermalexpansionterm.rb

計算結果

xz 図はこちら

実線 鉛直領域 30 km, 破線 鉛直領域 48 km, 点線 鉛直領域 24 km, 一点鎖線 鉛直領域 36 km	実線 鉛直領域 30 km, 破線 鉛直領域 48 km, 点線 鉛直領域 24 km, 一点鎖線 鉛直領域 36 km (鉛直領域 0-12km のみ)	基本場の圧力を用いて地表面フラックスを計算, 実線 鉛直領域 30 km, 破線 鉛直領域 48 km, 点線 鉛直領域 24 km, 一点鎖線 鉛直領域 36 km	基本場の圧力を用いて地表面フラックスを計算, 実線 鉛直領域 30 km, 破線 鉛直領域 48 km, 点線 鉛直領域 24 km, 一点鎖線 鉛直領域 36 km (鉛直領域 0-12km のみ)
鉛直領域 30 km (標準実験)	鉛直領域 24 km	鉛直領域 36 km	鉛直領域 48 km
基本場の圧力を用いて地表面フラックスを計算, 鉛直領域 30 km (標準実験)	基本場の圧力を用いて地表面フラックスを計算, 鉛直領域 24 km	基本場の圧力を用いて地表面フラックスを計算, 鉛直領域 36 km	基本場の圧力を用いて地表面フラックスを計算, 鉛直領域 48 km