[地球流体電脳倶楽部] [dcmodel | deepconv | プログラムのテスト計算[deepconv/arare5] ]

deepconv/arare5 : Odaka et al. (1998) の再計算 (3)

deepconv/arare5 を用いて Odaka et al. (1998) の再計算を行う. ただし以下の設定について Odaka et al. (1998) と異なる ( [ ] 内は odaka et al (1998) の設定).

• 水平計算領域は 20 km [51.2 km]
• 格子間隔は 200 m [100 m]
• 高度 7 km より上層にレイリー摩擦 [レイリー摩擦とニュートン冷却]
• 地表フラックス計算時に密度擾乱の寄与を考慮

Odaka et al. (1998) の結果こちらを参照

計算設定

• 計算領域
  • 水平 20 km, 鉛直 10 km
• 基本場の設定
  • 地表気温 245 K
  • 高度 5 km までは等温位, 5 km より上層は等温度
• 物理過程の設定
  • 乱流拡散係数は, Klemp and Wilhelmson (1978) のクロージャモデル
  • 高度 5 km 以下に水平一様な固定冷却
  • 高度 7 km より上層にレイリー摩擦, 時定数は 300 sec
  • 地表からの運動量と熱のフラックスはバルク法で計算, バルク係数は 0.01
• 大気成分
  • CO2 のみ
• 空間/時間格子間隔
  • 水平格子間隔 200 m ,鉛直格子間隔 200 m
  • 時間格子間隔
    • 長い時間ステップ: 4 sec
    • 短い時間ステップ: 0.4 sec
• 境界条件
  • 水平方向に周期境界, 上下境界で w = 0
  • 地表面温度 270 K
• 初期条件
  • u = w = Km = 0
  • 最下層に最大振幅 3 K の温位擾乱を乱数的に与える,
• 積分時間
  • 10 日

ソースコード, 設定ファイル

• ソースコード
  • 20120511 版ソースコードを修整したもの
  • tgz ファイル
• 設定ファイル
• 計算データ
• 描画スクリプト

結果一覧

場の量の分布
velz_864000.png 鉛直流 (t=864000 sec)	velx_864000.png 水平流 (t=864000 sec)	sf_864000.png 質量流線関数 (t=864000 sec)	pt_864000.png 温位偏差 (t=864000 sec)	km_864000.png 乱流拡散係数 (t=864000 sec)
水平平均量
ptmean_864000.png 温位 (水平平均, 2日毎)	velz2_864000.png 鉛直速度の2乗平均	pttend_864000.png 加熱率 (t=860600-864000 sec の平均)	pi_sfc_mean_864000.png 圧力関数 (z=100 m での水平平均)	sfcheat_864000.png 熱フラックス
保存量
ke_864000.png 全運動エネルギー	pe_864000.png 全ポテンシャルエネルギー