[deepconv] [プログラムのテスト計算 (deepconv/arare4)]

arare4 テスト計算: z 方向の移流計算 (3次元版,y方向一様)

概要

arare4 ソースコード中の x 方向の移流を計算するモジュールのテストを行う. 具体的には, 温位の予報方程式を計算する部分を改変し, テストを行う.

• 初期の温移分布はガウス分布
• z 方向の風速は一定 (20 m/sec)
• y 方向は一様
• 温位の移流のみ計算
• 乱流拡散なし

ソースコード

• arare4-20070810.tgz を改変したもの
  • メインプログラムは arare_3d_adv.f90
  • 温位の移流のみ計算
  • 乱流拡散なし

計算設定

• 設定ファイル
• 設定概要
  • 音波減衰項の係数
    • 2.0e-7
  • 計算領域
    • x, y, z 方向ともに 10 km
  • 格子点数
    • x 方向 50 点, y 方向 10 点, z 方向 50 点
  • 時間格子間隔
    • 長い時間ステップ: 4.0 sec, 短い時間ステップ 2.0e-1 sec
  • 境界条件
    • 水平, 鉛直ともに周期境界
  • 初期条件
    • u = 0 m/sec, v = 0 m/sec, w = 20 m/sec
    • θ = 300 K の等温位大気に振幅 0.1 K のガウス分布型の擾乱を加える
      • 温位擾乱は y 方向に一様
  • 計算モデル時間
    • 500 sec

計算結果 (温位)

pottemp_x5000_y5000_t000.png x=5000, y=5000, t=0 の分布	pottemp_x5000_y5000_t500.png x=5000, y=5000, t=500 の分布	pottemp_xz_t000.png y=5000, t=0 の分布
pottemp_xz_t500.png y=5000, t=500 の分布	pottemp_zt.png x=5000, y=5000 での zt ダイアグラム