格子セルごとに値を持つ計算結果

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

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

表 80 格子セルごとに値を持つ計算結果の出力に利用する関数

関数

備考

cg_iric_write_sol_cell_integer

整数の格子セルごとに値を持つ計算結果を出力する

cg_iric_write_sol_real

倍精度実数の格子セルごとに値を持つ計算結果を出力する
リスト 171 サンプルプログラム (格子セルごとに値を持つ計算結果)
 1program SampleProgram
 2  use iric
 3  implicit none
 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:: depth
11  integer, dimension(:,:), allocatable:: wetflag
12  character(len=20):: condFile
13
14  condFile = 'test.cgn'
15
16  ! CGNS ファイルのオープン
17  call cg_iric_open(condFile, IRIC_MODE_MODIFY, fin, ier)
18  if (ier /=0) STOP "*** Open error of CGNS file ***"
19
20  ! 格子のサイズを調べる
21  call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
22  ! 格子を読み込むためのメモリを確保
23  allocate(grid_x(isize, jsize), grid_y(isize, jsize))
24  ! 計算結果を保持するメモリも確保
25  allocate(depth(isize - 1, jsize - 1), wetflag(isize - 1, jsize - 1))
26  ! 格子を読み込む
27  call cg_iric_read_grid2d_coords(fin, grid_x, grid_y, ier)
28
29  ! 初期状態の情報を出力
30  time = 0
31  convergence = 0.1
32  call cg_iric_write_sol_start(fin, ier)
33  call cg_iric_write_sol_time(fin, time, ier)
34  call cg_iric_write_sol_cell_real(fin, 'Depth', depth, ier)
35  call cg_iric_write_sol_cell_integer(fin, 'Wet', wetflag, ier)
36  call cg_iric_write_sol_end(fin, ier)
37  do
38    time = time + 10.0
39
40    ! (ここで計算を実行)
41
42    call iric_check_cancel(canceled)
43    if (canceled == 1) exit
44
45    ! 計算結果を出力
46    call cg_iric_write_sol_start(fin, ier)
47    call cg_iric_write_sol_time(fin, time, ier)
48    call cg_iric_write_sol_cell_real(fin, 'Depth', depth, ier)
49    call cg_iric_write_sol_cell_integer(fin, 'Wet', wetflag, ier)
50    call cg_iric_write_sol_end(fin, ier)
51
52    if (time > 1000) exit
53  end do
54
55  ! CGNS ファイルのクローズ
56  call cg_iric_close(fin, ier)
57  stop
58end program SampleProgram