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

deepconv/arare5 内部重力波の伝播

deepconv/arare5 線形モデルによる Skamarock & Klemp (1994) の再計算.

参考文献

• Skamarock,W. C., Klemp,J.B., 1994: Efficienty and Accuracy of the Klemp-Wilhelmson Time-Splitting Technique. Mon. Wea. Rev., 122, 2623--2630

計算設定

• 計算領域
  • 水平 300 km, 鉛直 10 km
• 基本場の設定
  • 浮力振動数一定 (N = 10^-2 1/sec)
  • 一様な水平風 (U = 20 m/sec)
• 物理過程の設定
  • なし
• 大気成分
  • N2 のみ
• 空間/時間格子間隔
  • 水平格子間隔 1000 m ,鉛直格子間隔 1000 m
  • 時間格子間隔
    • 長い時間ステップ: 12 sec
    • 短い時間ステップ: 2 sec
• 境界条件
  • 水平方向に周期境界, 上下境界で w = 0
• 初期条件
  • u = w = π = 0
  • 温位擾乱は θ(x,z) = A sin (πz/H) / [1 + (x-x_c)^2/a^2]
    • A = 1.0e-2 K, x_c = 100 km, a = 5 km
• 積分時間
  • 3000 sec

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

• ソースコード
  • 20120828 版ソースコードを線形化したもの.
  • tgz ファイル
• 設定ファイル
• 計算データ
• 描画スクリプト

結果一覧

初期条件(温位)
温位擾乱 (t = 600 sec)	温位擾乱 (t = 1200 sec)	温位擾乱 (t = 1800 sec)	温位擾乱 (t = 2400 sec)	温位擾乱 (t = 3000 sec)
圧力擾乱 (t = 600 sec)	圧力擾乱 (t = 1200 sec)	圧力擾乱 (t = 1800 sec)	圧力擾乱 (t = 2400 sec)	圧力擾乱 (t = 3000 sec)
水平流 (t = 600 sec)	水平流 (t = 1200 sec)	水平流 (t = 1800 sec)	水平流 (t = 2400 sec)	水平流 (t = 3000 sec)
鉛直流 (t = 600 sec)	鉛直流 (t = 1200 sec)	鉛直流 (t = 1800 sec)	鉛直流 (t = 240 sec)	鉛直流 (t = 300 sec)