バグが直ったので再計算しました。 名古屋港の潮位変化(T.P 0) 潮汐は考慮していない。 名古屋港での最高潮位と時刻は概ね再現されていると思う。 若干、過剰評価しているので海面抵抗計算を補正する必要がありそうです。

やっと間違いが分かった。 後日、修正したコードで再計算してみようと思う。 さて、間違いは単純だったがこれの間違いが分かるのには相当時間が掛かった。 簡単に言うと、、、 int foo( int W, int H, int i, int j, Foo& f ); という関数の呼び出しを間違え…

計算してみた。 当時の名古屋港の観測値と全くといっていいほど合っていない事が分かった。 コンタで見る限りはそれなりの結果だったので今までは合ってると思っていたが、、、 ふと、思って名古屋港の潮位変化をグラフにしてみたらたった５０ｃｍしか上昇し…

スペクトルでレンダリング出来るようにしています。 レンダラが大体出来るとスペクトルへの対応は８割は終わっているようなものだと思います。 というのも処理の流れは全く同じなはずで若干のスペクトルへの対応があるだけですが残りの２割は以外に壁が厚い…

ようやくPSSMLTがちゃんと動くようになった。 OpenMPとMPIで並列化も出来た。 実は最近まで実装後の結果で正しい結果が何か分からなくなっていた。 穴空間「http://kagamin.net/hole/」で公開されているのと比較しても結果が異なるし、、、 試行錯誤している…