粒子の座標ごとに値を持つ計算結果 (複数グループ可)¶
粒子の座標ごとに値を持つ計算結果を出力する場合、 表 83 に示す関数を 使用します。
ここで示す関数を使うと、複数のグループの粒子を出力することができます。 各グループの出力の最初と最後で、 cg_iric_write_sol_particlegroup_groupbegin と cg_iric_write_sol_particlegroup_groupend を呼び出してください。
出力するプログラムの例は、リスト 174 を参照してください。
注釈
粒子の座標ごとに値を持つ計算結果 で示した関数とは異なり、 ここで示す関数では、一度の関数呼び出しでは粒子一つ分のデータを出力します。
関数 |
備考 |
---|---|
cg_iric_write_sol_particlegroup_groupbegin |
粒子の計算結果の出力を開始する |
cg_iric_write_sol_particlegroup_groupend |
粒子の計算結果の出力を終了する |
cg_iric_write_sol_particlegroup_pos2d |
粒子の位置を出力する (2次元) |
cg_iric_write_sol_particlegroup_pos3d |
粒子の位置を出力する (3次元) |
cg_iric_write_sol_particlegroup_integer |
整数の粒子ごとに値を持つ計算結果を出力する |
cg_iric_write_sol_particlegroup_real |
倍精度実数の粒子ごとに値を持つ計算結果を出力する |
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 integer:: numparticles = 10
11 double precision, dimension(:), allocatable:: particle_x, particle_y,
12 double precision, dimension(:), allocatable:: velocity_x, velocity_y, temperature
13 integer:: i
14 integer:: status = 1
15 character(len=20):: condFile
16
17 condFile = 'test.cgn'
18
19 ! CGNS ファイルのオープン
20 call cg_iric_open(condFile, IRIC_MODE_MODIFY, fin, ier)
21 if (ier /=0) STOP "*** Open error of CGNS file ***"
22
23 ! 格子のサイズを調べる
24 call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
25 ! 格子を読み込むためのメモリを確保
26 allocate(grid_x(isize, jsize), grid_y(isize, jsize))
27 ! 計算結果を保持するメモリを確保。
28 allocate(particle_x(numparticles), particle_y(numparticles))
29 allocate(velocity_x(numparticles), velocity_y(numparticles), temperature(numparticles))
30
31 ! 格子を読み込む
32 call cg_iric_read_grid2d_coords(fin, grid_x, grid_y, ier)
33
34 ! 初期状態の情報を出力
35 time = 0
36 call cg_iric_write_sol_start(fin, ier)
37 call cg_iric_write_sol_time(fin, time, ier)
38
39 call cg_iric_write_sol_particlegroup_groupbegin(fin, 'driftwood', ier)
40 do i = 1, numparticles
41 ! (ここで particle_x, particle_x, velocity_x, velocity_y, temperature に適切な値を設定)
42 call cg_iric_write_sol_particlegroup_pos2d(fin, particle_x(i), particle_y(i), ier)
43 call cg_iric_write_sol_particlegroup_real(fin, 'VelocityX', velocity_x(i), ier)
44 call cg_iric_write_sol_particlegroup_real(fin, 'VelocityY', velocity_y(i), ier)
45 call cg_iric_write_sol_particlegroup_real(fin, 'Temperature', temperature(i), ier)
46 end do
47 call cg_iric_write_sol_particlegroup_groupend(fin, ier)
48 call cg_iric_write_sol_end(fin, ier)
49
50 do
51 time = time + 10.0
52
53 ! (ここで計算を実行)
54
55 call iric_check_cancel(canceled)
56 if (canceled == 1) exit
57
58 ! 計算結果を出力
59 call cg_iric_write_sol_start(fin, ier)
60 call cg_iric_write_sol_time(fin, time, ier)
61 call cg_iric_write_sol_particlegroup_groupbegin(fin, 'driftwood', ier)
62 do i = 1, numparticles
63 ! (ここで particle_x, particle_x, velocity_x, velocity_y, temperature に適切な値を設定)
64 call cg_iric_write_sol_particlegroup_pos2d(fin, particle_x(i), particle_y(i), ier)
65 call cg_iric_write_sol_particlegroup_real(fin, 'VelocityX', velocity_x(i), ier)
66 call cg_iric_write_sol_particlegroup_real(fin, 'VelocityY', velocity_y(i), ier)
67 call cg_iric_write_sol_particlegroup_real(fin, 'Temperature', temperature(i), ier)
68 end do
69 call cg_iric_write_sol_particlegroup_groupend(fin, ier)
70 call cg_iric_write_sol_end(fin, ier)
71
72 if (time > 1000) exit
73 end do
74
75 ! CGNS ファイルのクローズ
76 call cg_iric_close(fin, ier)
77 stop
78end program SampleProgram