deepconv/arare4 : 主成分凝結を考慮した湿潤対流実験(中層一様冷却・下層一様加熱, 臨界飽和比 1.0 or 1.35, 地表面温度 165 K の場合)

下層の加熱と上層の冷却を釣り合わせて計算させている筈が, Scr = 1.0 の場合において雲が完全に蒸発した後でも, 温度が～ 1 K/day ずつ上昇するという現象が見られた. その原因を究明する為に, いくつか図を描き出してみた.

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 の出力に誤差が含まれているのか?

deepconv/arare4 : 主成分凝結を考慮した湿潤対流実験(中層一様冷却・下層一様加熱, 臨界飽和比 1.0 or 1.35, 地表面温度 165 K の場合)

計算設定

考察

結果一覧