平均操作ルーチン avr.f90 の動作をテストする. 特にここでは x フラックス格子点とスカラー格子点の変換を行う関数のテストを行う.
テストの際の注意点を以下にあげる.
テスト結果として元の分布, 平均をとった後の分布、およびその差の図を描く. 図は計算領域内のみのものと計算領域外を含むものの2種類を描く.
| プログラムソース | arare (微分・平均演算テスト用) 2004-10-04 版 |
| 行った操作 | u_avr_s( s_avr_u(u) ) |
| 元の関数の分布 | u = sin(2πx/Lx) sin(2πz/Lz) |
| 境界条件 | x: 周期境界, z:固定壁 |
| 格子点数 | x: 20, z: 20 |
| 出力ファイル | avr01_u.nc |
|
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 出力ファイル | avr01_u_mgn.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 行った操作 | u_avr_s( s_avr_u(u) ) を 5 回繰り返す |
| 元の関数の分布 | u = sin(2πx/Lx) sin(2πz/Lz) |
| 境界条件 | x: 周期境界, z:固定壁 |
| 格子点数 | x: 20, z: 20 |
| 出力ファイル | avr02_u.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 出力ファイル | avr02_u_mgn.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 行った操作 | u_avr_s( s_avr_u(u) ) |
| 元の関数の分布 | u = cos(2πx/Lx) cos(2πz/Lz) |
| 境界条件 | x: 周期境界, z:固定壁 |
| 格子点数 | x: 20, z: 20 |
| 出力ファイル | avr04_u.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 出力ファイル | avr04_u_mgn.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 行った操作 | u_avr_s( s_avr_u(u) ) を 5 回繰り返す |
| 元の関数の分布 | u = cos(2πx/Lx) cos(2πz/Lz) |
| 境界条件 | x: 周期境界, z:固定壁 |
| 格子点数 | x: 20, z: 20 |
| 出力ファイル | avr05_u.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
| 出力ファイル | avr05_u_mgn.nc |
|
|
|
| 元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |