[地球流体電脳倶楽部] [deepconv | 計算結果] [dcmodel-thum.rb]

Odaka et al. (1998) の計算

arare3m を用いて Odaka et al. (1998) の計算を行った. ただし, 計算モデル時間は 24 時間である.

計算設定

ソースプログラム arare3m-20051130
基本場の温度 高度 5 km 以下は温位 245 K, 高度 5 km 以上では温度 220 K
基本場の圧力 静水圧平衡 (地表面で 700 Pa)
温度の初期擾乱 最下層に最大振幅 2 K のランダムな擾乱を与える
圧力の初期擾乱 なし
風速の初期値 u = 0 m/s , w = 0 m/s
放射加熱 高度 5 km 以下で -50 K /day の一様冷却
地表面フラックス バルク法(バルク係数は熱, 運動量ともに 0.01 で一定, 地表面温度 270 K で固定)
重力加速度 3.72 m/s2
粘性 Klemp and Wilhelmson (1978). Cm = 0.2
時間フィルタ あり (係数は 0.1)
数値粘性 あり (5.0 x 10-4)
計算領域 水平 51.2 km x 鉛直 10 km
格子点間隔 水平 100 m x 鉛直 100 m
積分時間 24 時間
時間ステップ 長い: 1.5 秒, 短い: 0.15 秒
境界条件 水平: 周期, 鉛直: free-slip
数値解法 水平: 陽解法, 鉛直: 陰解法

実行環境と CPU 時間

実行環境 宇宙研 SX-6
CPU 時間 約 3 時間 3 分

計算結果



温位の水平平均値からのずれ


arare-odaka1998-dx100-20051130-PotTemp-t0.png
t = 0 hr


arare-odaka1998-dx100-20051130-PotTemp-t7200.png
t = 2 hr


arare-odaka1998-dx100-20051130-PotTemp-t14400.png
t = 4 hr


arare-odaka1998-dx100-20051130-PotTemp-t21600.png
t = 6 hr


arare-odaka1998-dx100-20051130-PotTemp-t28800.png
t = 8 hr


arare-odaka1998-dx100-20051130-PotTemp-t36000.png
t = 10 hr


arare-odaka1998-dx100-20051130-PotTemp-t43200.png
t = 12 hr


arare-odaka1998-dx100-20051130-PotTemp-t50400.png
t = 14 hr


arare-odaka1998-dx100-20051130-PotTemp-t57600.png
t = 16 hr


arare-odaka1998-dx100-20051130-PotTemp-t64800.png
t = 18 hr


arare-odaka1998-dx100-20051130-PotTemp-t72000.png
t = 20 hr


arare-odaka1998-dx100-20051130-PotTemp-t79200.png
t = 22 hr


arare-odaka1998-dx100-20051130-PotTemp-t86400.png
t = 24 hr



水平速度


arare-odaka1998-dx100-20051130-VelX-t0.png
t = 0 hr


arare-odaka1998-dx100-20051130-VelX-t7200.png
t = 2 hr


arare-odaka1998-dx100-20051130-VelX-t14400.png
t = 4 hr


arare-odaka1998-dx100-20051130-VelX-t21600.png
t = 6 hr


arare-odaka1998-dx100-20051130-VelX-t28800.png
t = 8 hr


arare-odaka1998-dx100-20051130-VelX-t36000.png
t = 10 hr


arare-odaka1998-dx100-20051130-VelX-t43200.png
t = 12 hr


arare-odaka1998-dx100-20051130-VelX-t50400.png
t = 14 hr


arare-odaka1998-dx100-20051130-VelX-t57600.png
t = 16 hr


arare-odaka1998-dx100-20051130-VelX-t64800.png
t = 18 hr


arare-odaka1998-dx100-20051130-VelX-t72000.png
t = 20 hr


arare-odaka1998-dx100-20051130-VelX-t79200.png
t = 22 hr


arare-odaka1998-dx100-20051130-VelX-t86400.png
t = 24 hr



鉛直速度


arare-odaka1998-dx100-20051130-VelZ-t0.png
t = 0 hr


arare-odaka1998-dx100-20051130-VelZ-t7200.png
t = 2 hr


arare-odaka1998-dx100-20051130-VelZ-t14400.png
t = 4 hr


arare-odaka1998-dx100-20051130-VelZ-t21600.png
t = 6 hr


arare-odaka1998-dx100-20051130-VelZ-t28800.png
t = 8 hr


arare-odaka1998-dx100-20051130-VelZ-t36000.png
t = 10 hr


arare-odaka1998-dx100-20051130-VelZ-t43200.png
t = 12 hr


arare-odaka1998-dx100-20051130-VelZ-t50400.png
t = 14 hr


arare-odaka1998-dx100-20051130-VelZ-t57600.png
t = 16 hr


arare-odaka1998-dx100-20051130-VelZ-t64800.png
t = 18 hr


arare-odaka1998-dx100-20051130-VelZ-t72000.png
t = 20 hr


arare-odaka1998-dx100-20051130-VelZ-t79200.png
t = 22 hr


arare-odaka1998-dx100-20051130-VelZ-t86400.png
t = 24 hr



乱流拡散係数(熱)


arare-odaka1998-dx100-20051130-Kh-t0.png
t = 0 hr


arare-odaka1998-dx100-20051130-Kh-t7200.png
t = 2 hr


arare-odaka1998-dx100-20051130-Kh-t14400.png
t = 4 hr


arare-odaka1998-dx100-20051130-Kh-t21600.png
t = 6 hr


arare-odaka1998-dx100-20051130-Kh-t28800.png
t = 8 hr


arare-odaka1998-dx100-20051130-Kh-t36000.png
t = 10 hr


arare-odaka1998-dx100-20051130-Kh-t43200.png
t = 12 hr


arare-odaka1998-dx100-20051130-Kh-t50400.png
t = 14 hr


arare-odaka1998-dx100-20051130-Kh-t57600.png
t = 16 hr


arare-odaka1998-dx100-20051130-Kh-t64800.png
t = 18 hr


arare-odaka1998-dx100-20051130-Kh-t72000.png
t = 20 hr


arare-odaka1998-dx100-20051130-Kh-t79200.png
t = 22 hr


arare-odaka1998-dx100-20051130-Kh-t86400.png
t = 24 hr



乱流拡散係数(運動量)


arare-odaka1998-dx100-20051130-Km-t0.png
t = 0 hr


arare-odaka1998-dx100-20051130-Km-t7200.png
t = 2 hr


arare-odaka1998-dx100-20051130-Km-t14400.png
t = 4 hr


arare-odaka1998-dx100-20051130-Km-t21600.png
t = 6 hr


arare-odaka1998-dx100-20051130-Km-t28800.png
t = 8 hr


arare-odaka1998-dx100-20051130-Km-t36000.png
t = 10 hr


arare-odaka1998-dx100-20051130-Km-t43200.png
t = 12 hr


arare-odaka1998-dx100-20051130-Km-t50400.png
t = 14 hr


arare-odaka1998-dx100-20051130-Km-t57600.png
t = 16 hr


arare-odaka1998-dx100-20051130-Km-t64800.png
t = 18 hr


arare-odaka1998-dx100-20051130-Km-t72000.png
t = 20 hr


arare-odaka1998-dx100-20051130-Km-t79200.png
t = 22 hr


arare-odaka1998-dx100-20051130-Km-t86400.png
t = 24 hr



エネルギーの時間変化


arare-odaka1998-dx100-20051130-EnergyTotal.png
領域全体のエネルギーの時間変化


arare-odaka1998-dx100-20051130-EnergyTend.png



arare-odaka1998-dx100-20051130-EnergyTend2.png



質量の時間変化

温位の水平平均値の時間変化

温位変化の寄与

arare-odaka1998-dx100-20051130-MassTotal.png
領域全体の質量の時間変化


arare-odaka1998-dx100-20051130-PotTempProf.png
t = 0 hr (紫) → 24 hr (赤) まで 2 時間ごと.


arare-odaka1998-dx100-20051130-PotTempTend.png


     

This page is generated by thum/dcmodel-thum-make.rb and thum/odaka1998-20051130-2.txt (2005/12/05 20:41:16 kitamo)
dcmodel-thum.rb [Script | Doc | Sample] is provided by dcmodel project in GFD Dennou Club
Copyright © GFD Dennou Club 2005