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

Odaka et al. (1998) の再計算

arare3m で Odaka et al. (1998) の計算を行った.

計算設定

ソースプログラム arare3m-20050810
基本場の温度 高度 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
粘性 あり
時間フィルタ あり (係数は 0.1)
数値粘性 あり (5.0 x 10-4)
計算領域 水平 51.2 km x 鉛直 10 km
格子点間隔 水平 100 m x 鉛直 100 m
積分時間 12 時間
時間ステップ 長い: 1.5 秒, 短い: 0.15 秒
境界条件 水平: 周期, 鉛直: free-slip
数値解法 水平: 陽解法, 鉛直: 陰解法

実行環境と計算時間

実行環境 宇宙研 SX-6
計算時間 1 時間 30 分

計算結果


温位の鉛直分布

熱(温位)の式の各項の寄与


arare-odaka1998-dx100-20050811-PotTemp-prof.png
右から順に t = 0, 2 hour, 4 hour, 6 hour, 8 hour, 10 hour, 12 hour.


arare-odaka1998-dx100-20050811-PotTempTend2.png
t = 11--12 hour の平均.



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

アニメーション[0--2 hour][10--12 hour]

arare-odaka1998-dx100-20050811-PotTemp-xm-t0.png
t = 0


arare-odaka1998-dx100-20050811-PotTemp-xm-t7200.png
t = 2 hour


arare-odaka1998-dx100-20050811-PotTemp-xm-t14400.png
t = 4 hour


arare-odaka1998-dx100-20050811-PotTemp-xm-t21600.png
t = 6 hour


arare-odaka1998-dx100-20050811-PotTemp-xm-t28800.png
t = 8 hour


arare-odaka1998-dx100-20050811-PotTemp-xm-t36000.png
t = 10 hour


arare-odaka1998-dx100-20050811-PotTemp-xm-t43200.png
t = 12 hour



水平風速


arare-odaka1998-dx100-20050811-VelX-t0.png
t = 0


arare-odaka1998-dx100-20050811-VelX-t7200.png
t = 2 hour


arare-odaka1998-dx100-20050811-VelX-t14400.png
t = 4 hour


arare-odaka1998-dx100-20050811-VelX-t21600.png
t = 6 hour


arare-odaka1998-dx100-20050811-VelX-t28800.png
t = 8 hour


arare-odaka1998-dx100-20050811-VelX-t36000.png
t = 10 hour


arare-odaka1998-dx100-20050811-VelX-t43200.png
t = 12 hour



鉛直風速


arare-odaka1998-dx100-20050811-VelZ-t0.png
t = 0


arare-odaka1998-dx100-20050811-VelZ-t7200.png
t = 2 hour


arare-odaka1998-dx100-20050811-VelZ-t14400.png
t = 4 hour


arare-odaka1998-dx100-20050811-VelZ-t21600.png
t = 6 hour


arare-odaka1998-dx100-20050811-VelZ-t28800.png
t = 8 hour


arare-odaka1998-dx100-20050811-VelZ-t36000.png
t = 10 hour


arare-odaka1998-dx100-20050811-VelZ-t43200.png
t = 12 hour



乱流拡散係数(Kh)


arare-odaka1998-dx100-20050811-Kh-t0.png
t = 0


arare-odaka1998-dx100-20050811-Kh-t7200.png
t = 2 hour


arare-odaka1998-dx100-20050811-Kh-t14400.png
t = 4 hour


arare-odaka1998-dx100-20050811-Kh-t21600.png
t = 6 hour


arare-odaka1998-dx100-20050811-Kh-t28800.png
t = 8 hour


arare-odaka1998-dx100-20050811-Kh-t36000.png
t = 10 hour


arare-odaka1998-dx100-20050811-Kh-t43200.png
t = 12 hour



乱流拡散係数(Km)


arare-odaka1998-dx100-20050811-Km-t0.png
t = 0


arare-odaka1998-dx100-20050811-Km-t7200.png
t = 2 hour


arare-odaka1998-dx100-20050811-Km-t14400.png
t = 4 hour


arare-odaka1998-dx100-20050811-Km-t21600.png
t = 6 hour


arare-odaka1998-dx100-20050811-Km-t28800.png
t = 8 hour


arare-odaka1998-dx100-20050811-Km-t36000.png
t = 10 hour


arare-odaka1998-dx100-20050811-Km-t43200.png
t = 12 hour

   

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