粒子の座標ごとに値を持つ計算結果

粒子の座標ごとに値を持つ計算結果を出力する場合、 表 74 に示す関数を 使用します。

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

注釈

ここで示す関数群は現在は非推奨です

粒子の座標ごとに値を持つ計算結果を出力する場合は、 粒子の座標ごとに値を持つ計算結果 (複数グループ可) に示す関数を使用する ことをおすすめします。そちらに示した関数を使用すれば、複数のグループの 粒子を出力でき、グループごとに色の設定や粒子のサイズを変えて可視化することができます。

表 74 粒子ごとに値を持つ計算結果の出力に利用する関数

関数

備考

cg_iric_write_sol_particle_pos2d_f

粒子の位置を出力する (2次元)

cg_iric_write_sol_particle_pos3d_f

粒子の位置を出力する (3次元)

cg_iric_write_sol_particle_integer_f

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

cg_iric_write_sol_particle_real_f

倍精度実数の粒子ごとに値を持つ計算結果を出力する
リスト 106 サンプルプログラム (粒子の座標ごとに値を持つ計算結果)
 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, 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  character(len=20):: condFile
14
15  condFile = 'test.cgn'
16
17  ! CGNS ファイルのオープン
18  call cg_open_f(condFile, CG_MODE_MODIFY, fin, ier)
19  if (ier /=0) STOP "*** Open error of CGNS file ***"
20
21  ! 内部変数の初期化
22  call cg_iric_init_f(fin, ier)
23  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
24
25  ! 格子のサイズを調べる
26  call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
27  ! 格子を読み込むためのメモリを確保
28  allocate(grid_x(isize, jsize), grid_y(isize, jsize))
29  ! 計算結果を保持するメモリを確保
30  allocate(particle_x(numparticles), particle_y(numparticles))
31  allocate(velocity_x(numparticles), velocity_y(numparticles), temperature(numparticles))
32
33  ! 格子を読み込む
34  call cg_iric_getgridcoord2d_f (grid_x, grid_y, ier)
35
36  ! 初期状態の情報を出力
37  time = 0
38  call cg_iric_write_sol_time_f(time, ier)
39  call cg_iric_write_sol_particle_pos2d_f(numparticles, particle_x, particle_y, ier)
40  call cg_iric_write_sol_particle_real_f('VelocityX', velocity_x, ier)
41  call cg_iric_write_sol_particle_real_f('VelocityY', velocity_y, ier)
42  call cg_iric_write_sol_particle_real_f('Temperature', temperature, ier)
43  do
44    time = time + 10.0
45
46    ! (ここで計算を実行)
47
48    call iric_check_cancel_f(canceled)
49    if (canceled == 1) exit
50
51    ! 計算結果を出力
52    call iric_write_sol_start_f(condFile, ier)
53    call cg_iric_write_sol_time_f(time, ier)
54    call cg_iric_write_sol_particle_pos2d_f(numparticles, particle_x, particle_y, ier)
55    call cg_iric_write_sol_particle_real_f('VelocityX', velocity_x, ier)
56    call cg_iric_write_sol_particle_real_f('VelocityY', velocity_y, ier)
57    call cg_iric_write_sol_particle_real_f('Temperature', temperature, ier)
58    call cg_iric_flush_f(condFile, fin, ier)
59    call iric_write_sol_end_f(condFile, ier)
60
61    if (time > 1000) exit
62  end do
63
64  ! CGNS ファイルのクローズ
65  call cg_close_f(fin, ier)
66  stop
67end program SampleProgram