時刻 (もしくはループ回数) の出力¶
CGNSファイルに、時刻もしくはループ回数を出力します。
その時刻での計算格子の出力や計算結果の出力を行うより前に、必ず実行してください。
また、時刻とループ回数を両方出力することはできません。必ずいずれかのみ出力してください。
関数 |
備考 |
---|---|
cg_iric_write_sol_start |
計算結果の出力を開始する |
cg_iric_write_sol_end |
計算結果の出力を終了する |
cg_iric_write_sol_time |
時刻を出力する |
cg_iric_write_sol_iteration |
ループ回数を出力する |
時刻を出力する処理の例を リスト 167 に示します。
1program Sample4
2 use iric
3 implicit none
4
5 integer:: fin, ier, i
6 double precision:: time
7
8 ! CGNS ファイルのオープン
9 call cg_iric_open('test.cgn', IRIC_MODE_MODIFY, fin, ier)
10 if (ier /=0) STOP "*** Open error of CGNS file ***"
11
12 ! 初期状態の情報を出力
13 time = 0
14
15 call cg_iric_write_sol_start(fin, ier)
16 call cg_iric_write_sol_time(fin, time, ier)
17 ! (ここで、初期の計算格子や計算結果を出力)
18 call cg_iric_write_sol_end(fin, ier)
19
20 do
21 time = time + 10.0
22 ! (ここで計算を実行)
23 call cg_iric_write_sol_start(fin, ier)
24 call cg_iric_write_sol_time(fin, time, ier)
25 ! (ここで、計算格子や計算結果を出力)
26 call cg_iric_write_sol_end(fin, ier)
27 If (time > 1000) exit
28 end do
29
30 ! CGNS ファイルのクローズ
31 call cg_iric_close(fin, ier)
32 stop
33end program Sample4