既存の計算結果の読み込み¶
既存のCGNSファイルに格納されている計算結果を読み込みます。
関数 |
備考 |
---|---|
cg_iric_read_sol_count |
計算結果の数を取得する |
cg_iric_read_sol_time |
計算結果の時刻の値を取得する |
cg_iric_read_sol_iteration |
計算結果のループ回数の値を取得する |
cg_iric_read_sol_baseiterative_integer |
整数の計算結果の値を取得する |
cg_iric_read_sol_baseiterative_real |
倍精度実数の計算結果の値を取得する |
cg_iric_read_sol_grid2d_coords |
計算結果の2次元構造格子を取得する |
cg_iric_read_sol_grid3d_coords |
計算結果の3次元構造格子を取得する |
cg_iric_read_sol_node_integer |
整数の格子点ごとに値を持つ計算結果の値を取得する |
cg_iric_read_sol_node_real |
倍精度実数の格子点ごとに値を持つ計算結果の値を取得する |
cg_iric_read_sol_cell_integer |
整数の格子セルごとに値を持つ計算結果の値を取得する |
cg_iric_read_sol_cell_real |
倍精度実数の格子セルごとに値を持つ計算結果の値を取得する |
既存のCGNSファイルを読み込み、格納されている計算結果を標準出力に 出力する処理の例を リスト 109 に示します。
1program SampleX
2 use iric
3 implicit none
4
5 integer:: fin, ier, isize, jsize, solid, solcount, iter, i, j
6 double precision, dimension(:,:), allocatable::grid_x, grid_y, result_real
7
8 ! CGNS ファイルのオープン
9 call cg_iric_open('test.cgn', IRIC_MODE_READ, fin, ier)
10 if (ier /=0) STOP "*** Open error of CGNS file ***"
11
12 ! 格子のサイズを調べる
13 call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
14
15 ! 計算結果を読み込むためのメモリを確保
16 allocate(grid_x(isize,jsize), grid_y(isize,jsize))
17 allocate(result_real(isize, jsize))
18
19 ! 計算結果を読み込み出力
20 call cg_iric_read_sol_count(fin, solcount, ier)
21 do solid = 1, solcount
22 call cg_iric_read_sol_iteration(fin, solid, iter, ier)
23 call cg_iric_read_sol_gridcoord2d(fin, solid, grid_x, grid_y, ier)
24 call cg_iric_read_sol_real(fin, solid, 'result_real', result_real, ier)
25
26 print *, 'iteration: ', iter
27 print *, 'grid_x, grid_y, result: '
28 do i = 1, isize
29 do j = 1, jsize
30 print *, '(', i, ', ', j, ') = (', grid_x(i, j), ', ', grid_y(i, j), ', ', result_real(i, j), ')'
31 end do
32 end do
33 end do
34
35 ! CGNS ファイルのクローズ
36 call cg_iric_close(fin, ier)
37 stop
38end program SampleX
なお、計算結果読み込みの関数を用いて、既存のCGNSファイルの計算結果を 分析・加工することができます(計算結果分析ソルバーの開発手順 参照)。