格子エッジごとに値を持つ計算結果

格子エッジごとに値を持つ計算結果を出力する場合、 表 73 に示す関数を 使用します。

出力するプログラムの例は、リスト 105 を参照してください。

表 73 格子エッジごとに値を持つ計算結果の出力に利用する関数

関数

備考

cg_iric_write_sol_iface_integer_f

整数のI方向の格子エッジごとに値を持つ計算結果を出力する

cg_iric_write_sol_iface_real_f

倍精度実数のI方向の格子エッジごとに値を持つ計算結果を出力する

cg_iric_write_sol_jface_integer_f

整数のJ方向の格子エッジごとに値を持つ計算結果を出力する

cg_iric_write_sol_jface_real_f

倍精度実数のJ方向の格子エッジごとに値を持つ計算結果を出力する
リスト 105 サンプルプログラム (格子エッジごとに値を持つ計算結果)
 1program SampleProgram
 2  implicit none
 3  include 'cgnslib_f.h'
 4
 5  integer:: fin, ier, isize, jsize
 6  integer:: canceled
 7  integer:: locked
 8  double precision:: time
 9  double precision, dimension(:,:), allocatable::grid_x, grid_y
10  double precision, dimension(:,:), allocatable:: fluxi, fluxj
11  character(len=20):: condFile
12
13  condFile = 'test.cgn'
14
15  ! CGNS ファイルのオープン
16  call cg_open_f(condFile, CG_MODE_MODIFY, fin, ier)
17  if (ier /=0) STOP "*** Open error of CGNS file ***"
18
19  ! 内部変数の初期化
20  call cg_iric_init_f(fin, ier)
21  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
22
23  ! 格子のサイズを調べる
24  call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
25  ! 格子を読み込むためのメモリを確保
26  allocate(grid_x(isize, jsize), grid_y(isize, jsize))
27  ! 計算結果を保持するメモリも確保
28  allocate(fluxi(isize, jsize - 1), fluxj(isize - 1, jsize))
29  ! 格子を読み込む
30  call cg_iric_getgridcoord2d_f (grid_x, grid_y, ier)
31
32  ! 初期状態の情報を出力
33  time = 0
34  convergence = 0.1
35  call cg_iric_write_sol_time_f(time, ier)
36  call cg_iric_write_sol_iface_real_f('FluxI', fluxi, ier)
37  call cg_iric_write_sol_jface_real_f('FluxJ', fluxj, ier)
38  do
39    time = time + 10.0
40
41    ! (ここで計算を実行)
42
43    call iric_check_cancel_f(canceled)
44    if (canceled == 1) exit
45
46    ! 計算結果を出力
47    call iric_write_sol_start_f(condFile, ier)
48    call cg_iric_write_sol_time_f(time, ier)
49    call cg_iric_write_sol_iface_real_f('FluxI', fluxi, ier)
50    call cg_iric_write_sol_jface_real_f('FluxJ', fluxj, ier)
51    call cg_iric_flush_f(condFile, fin, ier)
52    call iric_write_sol_end_f(condFile, ier)
53
54    if (time > 1000) exit
55  end do
56
57  ! CGNS ファイルのクローズ
58  call cg_close_f(fin, ier)
59  stop
60end program SampleProgram