昔のソースを使った場合に走るかどうかの確認. 以下のようになった.
spml-0.8.0 を使うと走らなくなるように見える.
発散成分にdamping をかけた計算 DivDampPeriodValue が 64 [day] 以下だと走らない. DivDampPeriodValue = 128 [day] だと落ちなかった.
今後の計算は DivDampPeriodValue = 128 [day] を用いて行うことにする.
これで T21 の計算を行い, 図をかくこことにする.
sjpack を使った計算も試すことにする. sjpack を使った方が速いかもしれない.
FFLAGS=-DSJPACK
を設定して make する.