計算格子の出力 (計算開始後の格子)

CGNSファイルに、計算開始後の計算格子を出力します。計算中に格子形状が 変化するソルバーでのみ行います。

特定の時間での計算格子を出力する前に、必ず 時刻 (もしくはループ回数) の出力 で示した時刻 (もしくはループ回数) の出力を行ってください。

以下に示す場合の格子の出力については、計算格子の出力 で示した関数を利用してください。

  • ソルバーで新たに格子を生成した
  • ソルバーで格子を再分割するなどして、次元や格子点数が異なる格子を生成した
  • 格子生成プログラム内で格子を生成した
表 68 利用する関数
関数 備考
cg_iric_write_sol_grid2d_coords 2次元構造格子を出力する
cg_iric_write_sol_grid3d_coords 3次元構造格子を出力する

2次元構造格子を出力する処理の例を リスト 101 に示します。

リスト 101 2次元構造格子を出力する処理の記述例
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
program Sample5
  use iric
  implicit none

  integer:: fin, ier, isize, jsize
  double precision:: time
  double precision, dimension(:,:), allocatable:: grid_x, grid_y

  ! CGNS ファイルのオープン
  call cg_iric_open('test.cgn', IRIC_MODE_MODIFY, fin, ier)
  if (ier /=0) STOP "*** Open error of CGNS file ***"

  ! 格子のサイズを調べる
  call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
  ! 格子を読み込むためのメモリを確保
  allocate(grid_x(isize,jsize), grid_y(isize,jsize))
  ! 格子を読み込む
  call cg_iric_read_grid2d_coords(fin, grid_x, grid_y, ier)

  ! 初期状態の情報を出力
  time = 0

  call cg_iric_write_sol_start(fin, ier)
  call cg_iric_write_sol_time(fin, time, ier)
  ! 格子を出力
  call cg_iric_write_sol_grid2d_coords(fin, grid_x, grid_y, ier)
  call cg_iric_write_sol_end(fin, ier)

  do
    time = time + 10.0
    ! (ここで計算を実行)
    call cg_iric_write_sol_start(fin, ier)
    call cg_iric_write_sol_time(fin, time, ier)
    call cg_iric_write_sol_grid2d_coords(fin, grid_x, grid_y, ier)
    call cg_iric_write_sol_end(fin, ier)
    If (time > 1000) exit
  end do

  ! CGNS ファイルのクローズ
  call cg_iric_close(fin, ier)
  stop
end program Sample5