更に, 南北風, ω, EP flux, 赤道上東西平均東西風 (sigma=1, 0.1, 0.01, 0.001) の時系列の図も作成する.
EP flux の作図には <URL:http://www.gfd-dennou.org/library/dcpam/related-program/util_epf-2012-11-20.tgz> を使う.
T21 はまだ実行時エラーが出る問題が解決しておらず, 実行することができていない.
この問題に対して想像されることは, T21 だと 400K 一様温度から計算を開始したときに ニュートン冷却で温度分布ができつつあるときに 不安定が起こって重力波が出て計算が壊れているのかもしれない, である.
そこで, 発散成分にdamping をかけることを試す. dcpam では既に実装されているので,
dynamics_hspl_vas83_nml
において,
FlagDivDamp = .true., DivDampPeriodValue = 2.0, DivDampPeriodUnit = 'day'
などを試してみる.