DCPAM 作業ミーティング記録 (2015/12/15)

参加者

DCPAM 最新公開版 (dcapm5-20150804) を使おうとしてみている.

もともと陸惑星計算で Abe et al. (2005) と違う結果が出ていた問題を新しいソースで確かめようとしていた. しかし, そもそも地球設定で走らない問題につきあたった. 今は, なぜ走らないのか原因を究明しようとしているところ.

現象としては, 北大のパソコンで走らせると, 走ったり走らなかったりする. 走らない場合はあっと言う間に落ちる (0.25 日も走らない). 落ちるときには,
- ** ERROR (Code -1000) [cloud_utils] *** Evaporation loop is not convergent, 101, NaN.
というメッセージが出てしまう. どういう場合で走ってどういう場合で走らないのか, その切り分けはまだできていない.
環境研の SXACE でも新しいソースコードを試したみた. 最新公開版 dcpam5-20150804 版は走らない. 堀之内から報告があった 0-dev エラーが起きて止まるようだ. 20151213 時点でのソースのスナップショット (これは上記の 0-dev エラーを修正したもの) は走るようだ. これについては, まだちゃんと確認していない.
TODO さらに, コンパイラを変える・マシンを変える・ノードを変える, とした場合で同じ問題が発生するかどうかを確認する.

観測された地形を与えた場合の DCPAM 計算で得られるダストデビルによるダスト巻き上げダストフラックスと水平風の分布との関係を見ている.
風ベクトルの動画を作ってみた. Ls=270 の時期. 地形高度の等値線と地表風ベクトルの重ね描き.

斜面の勾配の向きに沿って風が吹いているように見える領域もあれば, そうでない領域もある (斜面風で決まっている領域とそうでない領域があるように見える).
TODO
- 地表風を時間平均値と時間変動成分などに分解してみる.
- 風のベクトルの unit vector の取り方を工夫する (ちゃんと図の縦横比とあうようにする).