バグ取り

プログラムを徹底的に見直してみたところ、いくつかおかしなところを発見しました。
明らかな間違いも一カ所。
あと、ちょっと気になるところも一カ所。

あーーっもう、バグ取りはめんどくさいなあ。
これで修正したらすぐにうまく計算が走ってくれればまだいいのですが、なかなかそうはいきません。

ひとまずこれでプログラムミスは全て直したはずです。
ごちゃごちゃしていたところもすっきりと直したので、計算時間も大幅に短縮されたはずです。
あとは、ちゃんと上昇流が吹いてくれるかどうか・・・なんですよね、問題は。
なんで地面に向かって風が吹くんだろう？？

まずは非移流項の計算で u<0 の場合を場合分けするべきかどうかを考えてみることにします。
場合分けしないといけない気もするし、でも非移流項なんだからどっちでも同じ気もするし、、、まあ両方の場合で計算して確かめろって話ですよね。
めんどくさいけどやらなきゃ。

明日は１４時から東工大でセミナーです。
それまでにいくつか気になるところをいじってみるとしましょう。