計算格子の出力

CGNSファイルに、計算格子を出力します。

ソルバーでは、ソルバーで計算に用いる格子を生成する場合や、 2次元格子から3次元格子を生成する場合に行います。

格子生成プログラムでは必ず行います。

ここで示す関数は、ソルバーでは計算開始時の格子を出力するために使用します。 計算中に格子形状が変化する場合の格子の出力には、計算格子の出力 (計算開始後の格子) に示す関数を使用して下さい。

表 65 利用する関数

関数

備考

cg_iric_writegridcoord1d_f

1次元構造格子を出力する

cg_iric_writegridcoord2d_f

2次元構造格子を出力する

cg_iric_writegridcoord3d_f

3次元構造格子を出力する

cg_iric_write_grid_real_node_f

格子点で定義された整数の属性を出力する

cg_iric_write_grid_integer_node_f

格子点で定義された倍精度実数の属性を出力する

cg_iric_write_grid_real_cell_f

セルで定義された整数の属性を出力する

cg_iric_write_grid_integer_cell_f

セルで定義された倍精度実数の属性を出力する

2次元格子を読み込み、それを分割して生成した3次元格子を 出力する処理の記述例を リスト 99 に示します。

リスト 99 3次元格子を出力する処理の記述例
 1program Sample7
 2  implicit none
 3  include 'cgnslib_f.h'
 4
 5  integer:: fin, ier, isize, jsize, ksize, i, j, k, aret
 6  double precision:: time
 7  double precision:: convergence
 8  double precision, dimension(:,:), allocatable::grid_x, grid_y, elevation
 9  double precision, dimension(:,:,:), allocatable::grid3d_x, grid3d_y, grid3d_z
10  double precision, dimension(:,:,:), allocatable:: velocity, density
11
12  ! CGNS ファイルのオープン
13  call cg_open_f('test3d.cgn', CG_MODE_MODIFY, fin, ier)
14  if (ier /=0) STOP "*** Open error of CGNS file ***"
15
16  ! 内部変数の初期化
17  call cg_iric_init_f(fin, ier)
18  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
19
20  ! 格子のサイズを調べる
21  call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
22  ! 格子を読み込むためのメモリを確保
23  allocate(grid_x(isize,jsize), grid_y(isize,jsize), elevation(isize,jsize))
24  ! 格子を読み込む
25  call cg_iric_getgridcoord2d_f(grid_x, grid_y, ier)
26  call cg_iric_read_grid_real_node_f('Elevation', elevation, ier)
27
28  ! 読み込んだ2次元格子を元に、3次元格子を生成。
29  ! 3次元格子は Z方向に、深さ 5 で、5分割する
30
31  ksize = 6
32  allocate(grid3d_x(isize,jsize,ksize), grid3d_y(isize,jsize,ksize), grid3d_z(isize,jsize,ksize))
33  allocate(velocity(isize,jsize,ksize), STAT = aret)
34  print *, aret
35  allocate(density(isize,jsize,ksize), STAT = aret)
36  print *, aret
37  do i = 1, isize
38    do j = 1, jsize
39      do k = 1, ksize
40        grid3d_x(i,j,k) = grid_x(i,j)
41        grid3d_y(i,j,k) = grid_y(i,j)
42        grid3d_z(i,j,k) = elevation(i,j) + (k - 1)
43        velocity(i,j,k) = 0
44        density(i,j,k) = 0
45      end do
46    end do
47  end do
48  ! 生成した3次元格子を出力
49  call cg_iric_writegridcoord3d_f(isize, jsize, ksize, grid3d_x, grid3d_y, grid3d_z, ier)
50
51  ! 初期状態の情報を出力
52  time = 0
53  convergence = 0.1
54  call cg_iric_write_sol_time_f(time, ier)
55  ! 格子を出力
56  call cg_iric_write_sol_gridcoord3d_f(grid3d_x, grid3d_y, grid3d_z, ier)
57  ! 計算結果を出力
58  call cg_iric_write_sol_real_f('Velocity', velocity, ier)
59  call cg_iric_write_sol_real_f('Density', density, ier)
60  call cg_iric_write_sol_baseiterative_real_f ('Convergence', convergence, ier)
61
62
63  do
64    time = time + 10.0
65    ! (ここで計算を実行。格子の形状も変化)
66    call cg_iric_write_sol_time_f(time, ier)
67    ! 格子を出力
68    call cg_iric_write_sol_gridcoord3d_f(grid3d_x, grid3d_y, grid3d_z, ier)
69    ! 計算結果を出力
70    call cg_iric_write_sol_real_f('Velocity', velocity, ier)
71    call cg_iric_write_sol_real_f('Density', density, ier)
72    call cg_iric_write_sol_baseiterative_real_f ('Convergence', convergence, ier)
73
74    If (time > 100) exit
75  end do
76
77  ! CGNS ファイルのクローズ
78  call cg_close_f(fin, ier)
79  stop
80end program Sample7