熱対流の数値計算を行なうプログラムです。 性懲りも無くまたやってしまいました。 とりあえず作ってみたんですが、こういう分野は未だ FORTRAN が現役で、Ｃ言語で記述されたものは意外と 公開されていない(？)ので公開しました。とはいっても 解法はいちばん単純なＦＴＣＳのオイラー陽解法です。 研究などで利用される方、分かる方は使ってください。 実装はかなりプリリミナリですが、そこらへんの 数値計算の書籍に付録してあるプログラム より多少はマシかもしれません（と願いつつ...）。
アーカイブを展開するときは、フォルダを１つ作って、その中で展開してください。
上のプログラムには実行ファイルが含まれています。 トレンドマイクロ社の オンラインスキャンではウイルスは検出されませんでしたが、不安な方は アーカイブ展開後実行ファイルを削除し、ご自分でコンパイルを行なってください。 なお、開発環境はＧＣＣです。
以下、熱対流（自然対流）問題について簡単に説明します。
扱う問題はこんな感じです。非圧縮粘性流体を考えています。
以下、Ra=15000, Pr=7.1 で計算したアニメーション表示のスクリーンショットです。 縦横比が１：１になっていませんが、正方形での計算です （いつものクセでうっかり画面の width と height を変えてしまいました）。 上でダウンロードできるアニメーションプログラムは width : height = 1 : 1 に修正しました。
