Numerical Solution By Explicit-FTCS Method - 6
■ プログラム ■
前項の手計算をプログラムにより計算します。
| プログラム | 説明 |
| exftcs.f | 移流方程式を陽解式FTCS差分法により近似解を求めるプログラムソースです。FORTRANの文法書を参考にしつつ，十分に理解してください。 |
| プログラム中で使われている変数 | ||
| 変数 | 変数の型 | 説明 |
| C(J，N) | 倍精度型 二次元配列 |
節点上でのcの値。例えば、c12値は，プログラム中でC(2，1)と表現される。 |
| MXNX | 整数型 | 二次元配列C(J，N)におけるJの配列数の上限値。 |
| MXNT | 整数型 | 二次元配列C(J，N)におけるNの配列数の上限値。。 |
| CHARTRIAL | 文字列型 | 出力ファイルに番号をつけるために使う。2文字の文字列を入力する必要がある（日本語不可）。 |
| INPDATA | 文字列型 | 出力ファイルパス「./input__.dat」を示すもの。 |
| OUTDATA | 文字列型 | 出力ファイルパス「./output__.dat」を示すもの。 |
| プログラム中で使われているサブルーチン | |
| サブルーチン | 説明 |
| INPUT | 解析空間XX，解析時間TT，空間刻み幅⊿x，時間刻み幅⊿t，移流速度u及びファイル番号の入力を行い，総空間ステップ数NX，総時間ステップ数NT及び無次元パラメータσを計算します。 |
| SOLVER | 移流方程式を陽解式FTCS差分法により計算を行います。 |
| OUTPUT | 出力ファイルを2つ作成します（入力パラメータを記録したもの，解析結果を記録したもの）。 |
| Back | Next | |||
| Advection Equation | Exact Solution |
Numerical Solution By Ex-FTCS Method |
Numerical Solution By Ex-CTCS Method |