いかんね〜.
プログラムをちゃんと整理していなかったツケが,
この忙しい時期になってまわってきた!!
多分,問題なのは昔取り上げた緯度経度関係.
OFESデータはのデータ間隔は0.1度グリッドになっている.
一方で,流速データの単位はcm/secなので,
緯度経度を距離に直してやる必要があるのだ〜.
だから粒子追跡のプログラムは,
緯度経度,距離に着目すると,
1. 粒子投入位置を緯度経度で読み込み,原点からの距離(cm)に変換
2. 流速データから移流,拡散方程式を解いて粒子の挙動を計算
3. 書き出す際に原点からの距離(cm)を緯度経度に戻す
という要領になっている.
たぶんここの書き間違いによってうまくいっていないんだと思う.
少し混乱しているので,明日なんとかします!!
さてと,ジムに行ってきます!
0 件のコメント:
コメントを投稿