計算格子の出力¶
CGNSファイルに、計算格子を出力します。
ソルバーでは、ソルバーで計算に用いる格子を生成する場合や、 2次元格子から3次元格子を生成する場合に行います。
格子生成プログラムでは必ず行います。
ここで示す関数は、ソルバーでは計算開始時の格子を出力するために使用します。 計算中に格子形状が変化する場合の格子の出力には、計算格子の出力 (計算開始後の格子) に示す関数を使用して下さい。
関数 |
備考 |
---|---|
cg_iric_write_grid1d_coords |
1次元構造格子を出力する |
cg_iric_write_grid2d_coords |
2次元構造格子を出力する |
cg_iric_write_grid3d_coords |
3次元構造格子を出力する |
cg_iric_write_grid_real_node |
格子点で定義された整数の属性を出力する |
cg_iric_write_grid_integer_node |
格子点で定義された倍精度実数の属性を出力する |
cg_iric_write_grid_real_cell |
セルで定義された整数の属性を出力する |
cg_iric_write_grid_integer_cell |
セルで定義された倍精度実数の属性を出力する |
2次元格子を読み込み、それを分割して生成した3次元格子を 出力する処理の記述例を リスト 166 に示します。
1program Sample7
2 use iric
3 implicit none
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_iric_open('test3d.cgn', IRIC_MODE_MODIFY, fin, ier)
14 if (ier /=0) STOP "*** Open error of CGNS file ***"
15
16 ! 格子のサイズを調べる
17 call cg_iric_read_grid2d_str_size(fin, isize, jsize, ier)
18 ! 格子を読み込むためのメモリを確保
19 allocate(grid_x(isize,jsize), grid_y(isize,jsize), elevation(isize,jsize))
20 ! 格子を読み込む
21 call cg_iric_read_grid2d_coords(fin, grid_x, grid_y, ier)
22 call cg_iric_read_grid_real_node(fin, 'Elevation', elevation, ier)
23
24 ! 読み込んだ2次元格子を元に、3次元格子を生成。
25 ! 3次元格子は Z方向に、深さ 5 で、5分割する
26
27 ksize = 6
28 allocate(grid3d_x(isize,jsize,ksize), grid3d_y(isize,jsize,ksize), grid3d_z(isize,jsize,ksize))
29 allocate(velocity(isize,jsize,ksize), STAT = aret)
30 print *, aret
31 allocate(density(isize,jsize,ksize), STAT = aret)
32 print *, aret
33 do i = 1, isize
34 do j = 1, jsize
35 do k = 1, ksize
36 grid3d_x(i,j,k) = grid_x(i,j)
37 grid3d_y(i,j,k) = grid_y(i,j)
38 grid3d_z(i,j,k) = elevation(i,j) + (k - 1)
39 velocity(i,j,k) = 0
40 density(i,j,k) = 0
41 end do
42 end do
43 end do
44 ! 生成した3次元格子を出力
45 call cg_iric_write_grid3d_coords(fin, isize, jsize, ksize, grid3d_x, grid3d_y, grid3d_z, ier)
46
47 ! 初期状態の情報を出力
48 time = 0
49 convergence = 0.1
50 call cg_iric_write_sol_start(fin, ier)
51 call cg_iric_write_sol_time(fin, time, ier)
52 ! 格子を出力
53 call cg_iric_write_sol_grid3d_coords(fin, grid3d_x, grid3d_y, grid3d_z, ier)
54 ! 計算結果を出力
55 call cg_iric_write_sol_node_real(fin, 'Velocity', velocity, ier)
56 call cg_iric_write_sol_node_real(fin, 'Density', density, ier)
57 call cg_iric_write_sol_baseiterative_real(fin, 'Convergence', convergence, ier)
58 call cg_iric_write_sol_end(fin, ier)
59
60 do
61 time = time + 10.0
62 ! (ここで計算を実行。格子の形状も変化)
63 call cg_iric_write_sol_start(fin, ier)
64 call cg_iric_write_sol_time(fin, time, ier)
65 ! 格子を出力
66 call cg_iric_write_sol_grid3d_coords(fin, grid3d_x, grid3d_y, grid3d_z, ier)
67 ! 計算結果を出力
68 call cg_iric_write_sol_node_real(fin, 'Velocity', velocity, ier)
69 call cg_iric_write_sol_node_real(fin, 'Density', density, ier)
70 call cg_iric_write_sol_baseiterative_real(fin, 'Convergence', convergence, ier)
71 call cg_iric_write_sol_end(fin, ier)
72
73 If (time > 100) exit
74 end do
75
76 ! CGNS ファイルのクローズ
77 call cg_iric_close(fin, ier)
78 stop
79end program Sample7