1つのタイムステップで1つ値を持つ計算結果

1つのタイムステップで1つ値を持つ計算結果を出力する場合、 表 78 に示す関数を 使用します。

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

表 78 1つのタイムステップで1つ値を持つ計算結果の出力に利用する関数

関数

備考

cg_iric_write_sol_baseiterative_integer

整数の計算結果を出力する

cg_iric_write_sol_baseiterative_real

倍精度実数の計算結果を出力する

cg_iric_write_sol_baseiterative_string

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