ポリゴンもしくは折れ線ごとに値を持つ計算結果

ポリゴンもしくは折れ線ごとに値を持つ計算結果を出力する場合、 表 76 に示す関数を 使用します。

ポリゴンもしくは折れ線では、複数のグループを出力することができます。 各グループの出力の最初と最後で、 cg_iric_write_sol_polydata_groupbegin_f と cg_iric_write_sol_polydata_groupend_f を呼び出してください。

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

注釈

粒子の座標ごとに値を持つ計算結果を出力する関数では、 一度の関数の呼び出して全ての粒子の座標や値を出力しますが、 ポリゴンもしくは折れ線ごとに値を持つ計算結果では、一度の関数呼び出しでは ポリゴンもしくは折れ線一つ分のデータを出力します。

注釈

ポリゴンもしくは折れ線ごとに値を持つ計算結果は、2次元にのみ対応しています。

注釈

一つのグループの中にポリゴンと折れ線を混在させることもできます。

注釈

ポリゴンもしくは折れ線では、計算結果の値はスカラー量にのみ対応しています。

表 76 ポリゴンもしくは折れ線ごとに値を持つ計算結果の出力に利用する関数

関数

備考

cg_iric_write_sol_polydata_groupbegin_f

ポリゴンもしくは折れ線で定義された計算結果の出力を開始する

cg_iric_write_sol_polydata_groupend_f

ポリゴンもしくは折れ線で定義された計算結果の出力を終了する

cg_iric_write_sol_polydata_polygon_f

計算結果としてポリゴンの形状を出力する

cg_iric_write_sol_polydata_polyline_f

計算結果として折れ線の形状を出力する

cg_iric_write_sol_polydata_integer_f

整数のポリゴンもしくは折れ線ごとに値を持つ計算結果を出力する

cg_iric_write_sol_polydata_real_f

倍精度実数のポリゴンもしくは折れ線ごとに値を持つ計算結果を出力する
リスト 108 サンプルプログラム (ポリゴンもしくは折れ線ごとに値を持つ計算結果)
 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:: numpolygons = 10
11  integer:: numpoints = 5
12  double precision, dimension(:), allocatable:: polydata_x, polydata_y,
13  double precision:: temperature = 26
14  integer:: i
15  integer:: status = 1
16  character(len=20):: condFile
17
18  condFile = 'test.cgn'
19
20  ! CGNS ファイルのオープン
21  call cg_open_f(condFile, CG_MODE_MODIFY, fin, ier)
22  if (ier /=0) STOP "*** Open error of CGNS file ***"
23
24  ! 内部変数の初期化
25  call cg_iric_init_f(fin, ier)
26  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
27
28  ! 格子のサイズを調べる
29  call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
30  ! 格子を読み込むためのメモリを確保
31  allocate(grid_x(isize, jsize), grid_y(isize, jsize))
32  ! 計算結果を保持するメモリを確保。一つのポリゴンの点数は5点
33  allocate(polydata_x(numpoints), polydata_y(numpoints))
34
35  ! 格子を読み込む
36  call cg_iric_getgridcoord2d_f(grid_x, grid_y, ier)
37
38  ! 初期状態の情報を出力
39  time = 0
40  call cg_iric_write_sol_time_f(time, ier)
41
42  call cg_iric_write_sol_polydata_groupbegin_f('fish', ier)
43  do i = 1, numpolygons
44    ! (ここで polydata_x, polydata_y, temperature, status に適切な値を設定)
45    call cg_iric_write_sol_polydata_polygon_f(numpoints, polydata_x, polydata_y, ier)
46    call cg_iric_write_sol_polydata_real_f('Temperature', temperature, ier)
47    call cg_iric_write_sol_polydata_integer_f('Status', status, ier)
48  end do
49  call cg_iric_write_sol_polydata_groupend_f(ier)
50
51  do
52    time = time + 10.0
53
54    ! (ここで計算を実行)
55
56    call iric_check_cancel_f(canceled)
57    if (canceled == 1) exit
58
59    ! 計算結果を出力
60    call iric_write_sol_start_f(condFile, ier)
61    call cg_iric_write_sol_time_f(time, ier)
62    call cg_iric_write_sol_polydata_groupbegin_f('fish', ier)
63    do i = 1, numpolygons
64      ! (ここで polydata_x, polydata_y, temperature, status に適切な値を設定)
65      call cg_iric_write_sol_polydata_polygon_f(numpoints, polydata_x, polydata_y, ier)
66      call cg_iric_write_sol_polydata_real_f('Temperature', temperature, ier)
67      call cg_iric_write_sol_polydata_integer_f('Status', status, ier)
68    end do
69    call cg_iric_write_sol_polydata_groupend_f(ier)
70
71    if (time > 1000) exit
72  end do
73
74  ! CGNS ファイルのクローズ
75  call cg_close_f(fin, ier)
76  stop
77end program SampleProgram