1つのタイムステップで1つ値を持つ計算結果¶
1つのタイムステップで1つ値を持つ計算結果を出力する場合、 表 70 に示す関数を 使用します。
出力するプログラムの例は、リスト 102 を参照してください。
関数 |
備考 |
|---|---|
cg_iric_write_sol_baseiterative_integer_f |
整数の計算結果を出力する |
cg_iric_write_sol_baseiterative_real_f |
倍精度実数の計算結果を出力する |
cg_iric_write_sol_baseiterative_string_f |
文字列の計算結果を出力する |
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:: 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_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 call cg_iric_getgridcoord2d_f (grid_x, grid_y, ier)
29
30 ! 初期状態の情報を出力
31 time = 0
32 convergence = 0.1
33 call cg_iric_write_sol_time_f(time, ier)
34 call cg_iric_write_sol_baseiterative_real_f('Convergence', convergence, ier)
35 do
36 time = time + 10.0
37
38 ! (ここで計算を実行)
39
40 call iric_check_cancel_f(canceled)
41 if (canceled == 1) exit
42
43 ! 計算結果を出力
44 call iric_write_sol_start_f(condFile, ier)
45 call cg_iric_write_sol_time_f(time, ier)
46 call cg_iric_write_sol_baseiterative_real_f('Convergence', convergence, ier)
47 call cg_iric_flush_f(condFile, fin, ier)
48 call iric_write_sol_end_f(condFile, ier)
49
50 if (time > 1000) exit
51 end do
52
53 ! CGNS ファイルのクローズ
54 call cg_close_f(fin, ier)
55 stop
56end program SampleProgram