下層の加熱と上層の冷却を釣り合わせて計算させている筈が, Scr = 1.0
の場合において雲が完全に蒸発した後でも, 温度が 〜 1 K/day
ずつ上昇するという現象が見られた.
その原因を究明する為に, いくつか図を描き出してみた.
解析の結果, 以下の 2 つが原因であることが判明した.
rand_CS1_coolheatbalance_DensCloud_00000000.png
|
rand_CS1_coolheatbalance_DensCloud_00086400.png
|
rand_CS1_coolheatbalance_DensCloud_00172800.png
|
rand_CS1_coolheatbalance_DensCloud_00259200.png
|
rand_CS1_coolheatbalance_DensCloud_00345600.png
|
rand_CS1_coolheatbalance_DensCloud_00432000.png
|
rand_CS1_coolheatbalance_DensCloud_00864000.png
|
rand_CS1_coolheatbalance_DensCloud_01728000.png
|
rand_CS1_coolheatbalance_DensCloud_02592000.png
|
rand_CS1_coolheatbalance_DensCloud_03456000.png
|
rand_CS1_coolheatbalance_DensCloud_04320000.png
|
rand_CS1_coolheatbalance_PotTempDev_00000000.png
|
rand_CS1_coolheatbalance_PotTempDev_00086400.png
|
rand_CS1_coolheatbalance_PotTempDev_00172800.png
|
rand_CS1_coolheatbalance_PotTempDev_00259200.png
|
rand_CS1_coolheatbalance_PotTempDev_00345600.png
|
rand_CS1_coolheatbalance_PotTempDev_00432000.png
|
rand_CS1_coolheatbalance_PotTempDev_00864000.png
|
rand_CS1_coolheatbalance_PotTempDev_01728000.png
|
rand_CS1_coolheatbalance_PotTempDev_02592000.png
|
rand_CS1_coolheatbalance_PotTempDev_03456000.png
|
rand_CS1_coolheatbalance_PotTempDev_04320000.png
|
rand_CS1_coolheatbalance_VelX_00000000.png
|
rand_CS1_coolheatbalance_VelX_00086400.png
|
rand_CS1_coolheatbalance_VelX_00172800.png
|
rand_CS1_coolheatbalance_VelX_00259200.png
|
rand_CS1_coolheatbalance_VelX_00345600.png
|
rand_CS1_coolheatbalance_VelX_00432000.png
|
rand_CS1_coolheatbalance_VelX_00864000.png
|
rand_CS1_coolheatbalance_VelX_01728000.png
|
rand_CS1_coolheatbalance_VelX_02592000.png
|
rand_CS1_coolheatbalance_VelX_03456000.png
|
rand_CS1_coolheatbalance_VelX_04320000.png
|
rand_CS1_coolheatbalance_VelZ_00000000.png
|
rand_CS1_coolheatbalance_VelZ_00086400.png
|
rand_CS1_coolheatbalance_VelZ_00172800.png
|
rand_CS1_coolheatbalance_VelZ_00259200.png
|
rand_CS1_coolheatbalance_VelZ_00345600.png
|
rand_CS1_coolheatbalance_VelZ_00432000.png
|
rand_CS1_coolheatbalance_VelZ_00864000.png
|
rand_CS1_coolheatbalance_VelZ_01728000.png
|
rand_CS1_coolheatbalance_VelZ_02592000.png
|
rand_CS1_coolheatbalance_VelZ_03456000.png
|
rand_CS1_coolheatbalance_VelZ_04320000.png
|
rand_CS135_coolheatbalance_DensCloud_00000000.png
|
rand_CS135_coolheatbalance_DensCloud_00086400.png
|
rand_CS135_coolheatbalance_DensCloud_00172800.png
|
rand_CS135_coolheatbalance_DensCloud_00259200.png
|
rand_CS135_coolheatbalance_DensCloud_00345600.png
|
rand_CS135_coolheatbalance_DensCloud_00432000.png
|
rand_CS135_coolheatbalance_DensCloud_00864000.png
|
rand_CS135_coolheatbalance_DensCloud_01728000.png
|
rand_CS135_coolheatbalance_DensCloud_02592000.png
|
rand_CS135_coolheatbalance_DensCloud_03456000.png
|
rand_CS135_coolheatbalance_DensCloud_04320000.png
|
rand_CS135_coolheatbalance_PotTempDev_00000000.png
|
rand_CS135_coolheatbalance_PotTempDev_00086400.png
|
rand_CS135_coolheatbalance_PotTempDev_00172800.png
|
rand_CS135_coolheatbalance_PotTempDev_00259200.png
|
rand_CS135_coolheatbalance_PotTempDev_00345600.png
|
rand_CS135_coolheatbalance_PotTempDev_00432000.png
|
rand_CS135_coolheatbalance_PotTempDev_00864000.png
|
rand_CS135_coolheatbalance_PotTempDev_01728000.png
|
rand_CS135_coolheatbalance_PotTempDev_02592000.png
|
rand_CS135_coolheatbalance_PotTempDev_03456000.png
|
rand_CS135_coolheatbalance_PotTempDev_04320000.png
|
rand_CS135_coolheatbalance_VelX_00000000.png
|
rand_CS135_coolheatbalance_VelX_00086400.png
|
rand_CS135_coolheatbalance_VelX_00172800.png
|
rand_CS135_coolheatbalance_VelX_00259200.png
|
rand_CS135_coolheatbalance_VelX_00345600.png
|
rand_CS135_coolheatbalance_VelX_00432000.png
|
rand_CS135_coolheatbalance_VelX_00864000.png
|
rand_CS135_coolheatbalance_VelX_01728000.png
|
rand_CS135_coolheatbalance_VelX_02592000.png
|
rand_CS135_coolheatbalance_VelX_03456000.png
|
rand_CS135_coolheatbalance_VelX_04320000.png
|
rand_CS135_coolheatbalance_VelZ_00000000.png
|
rand_CS135_coolheatbalance_VelZ_00086400.png
|
rand_CS135_coolheatbalance_VelZ_00172800.png
|
rand_CS135_coolheatbalance_VelZ_00259200.png
|
rand_CS135_coolheatbalance_VelZ_00345600.png
|
rand_CS135_coolheatbalance_VelZ_00432000.png
|
rand_CS135_coolheatbalance_VelZ_00864000.png
|
rand_CS135_coolheatbalance_VelZ_01728000.png
|
rand_CS135_coolheatbalance_VelZ_02592000.png
|
rand_CS135_coolheatbalance_VelZ_03456000.png
|
rand_CS135_coolheatbalance_VelZ_04320000.png
|
CS1-Temp-overplot.png
Scr=1.0 のときの温度分布を 10 日おきに重ね書きしたもの. 約 1 K/day で温度が上昇していることが分かる.
|
CS1-thetatend_49-50days.png
Scr=1.0, 49 -- 50 日目における温位の tendency. 赤は移流項, 青は放射加熱項, 緑は凝結加熱項, 黒は tendency の総和を表す.
|
CS1-RadHeatIntegral_45-50days.png
Scr=1.0 のときの温度に関する放射加熱の寄与の全領域積分(45 -- 50 日). 正値になっていることが分かる.
|
CS135-Temp-overplot.png
Scr=1.35 のときの温度分布を 10 日おきに重ね書きしたもの. 10 日目以降は温度がほとんど上昇していないことが分かる.
|
CS135-thetatend_49-50days.png
Scr=1.35, 49 -- 50 日目における温位の tendency. 赤は移流項, 青は放射加熱項, 緑は凝結加熱項, 黒は tendency の総和を表す.
|
CS135-RadHeatIntegral_40-50days.png
Scr=1.35 のときの温度に関する放射加熱の寄与の全領域積分(40 -- 50 日). 正値になっていることが分かる.
|
CS1-RadHeatIntegral_900sec.png
改良した放射のプログラムを用いたテスト計算での全領域積分値(Scr=1.0, 900 sec 積分). 加熱・冷却は完全に釣り合ってはいないものの, 加熱と冷却の差は 1 / 1000 程度(温度変化に換算すると 〜 0.001 K/day )となっていることが分かる .
|
CS1-RadHeatIntegral_900sec_v2.png
改良版(第 2 弾)でのテスト計算での残領域積分値. 陽に加熱と冷却を合わせるようにしたはずなのだが, 何故か 100 W 程度の差が出る. tendency の出力に誤差が含まれているのか?
|
|
|
|