平均操作ルーチン avr.f90 の動作をテストする. 特にここでは z フラックス格子点とスカラー格子点の変換を行う関数のテストを行う.
テストの際の注意点を以下にあげる.
テスト結果として元の分布, 平均をとった後の分布、およびその差の図を描く. 図は計算領域内のみのものと計算領域外を含むものの2種類を描く.
プログラムソース | arare_test_2004-09-29 |
行った操作 | w_avr_s( s_avr_w(w) ) |
元の関数の分布 | w = sin(2πx/Lx) sin(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr01_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr01_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
行った操作 | w_avr_s( s_avr_w(w) ) を 5 回繰り返す |
元の関数の分布 | w = sin(2πx/Lx) sin(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr02_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr02_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
行った操作 | w_avr_s( s_avr_w(w) ) を 30 回繰り返す |
元の関数の分布 | w = sin(2πx/Lx) sin(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr03_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr03_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
行った操作 | w_avr_s( s_avr_w(w) ) |
元の関数の分布 | w = cos(2πx/Lx) cos(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr04_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr04_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
行った操作 | w_avr_s( s_avr_w(w) ) を 5 回繰り返す |
元の関数の分布 | w = cos(2πx/Lx) cos(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr05_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr05_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
行った操作 | w_avr_s( s_avr_w(w) ) を 30 回繰り返す |
元の関数の分布 | w = cos(2πx/Lx) cos(2πz/Lz) |
境界条件 | x: 周期境界, z:固定壁 |
格子点数 | x: 20, z: 20 |
出力ファイル | avr06_w.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |
出力ファイル | avr06_w_mgn.nc |
![]() |
![]() |
![]() |
元の関数の分布 | 平均操作後の分布 | 操作前と後の差 |