格子の出力処理の記述¶
格子の出力処理を記述します。
まずは、iRIC との連携が正しく行えることを確認するため、 単純な格子を生成して出力する処理を記述します。
格子を出力する処理を追記したソースコードを リスト 20 に示します。追記した部分を強調して示します。
1program SampleProgram
2 use iric
3 implicit none
4
5 integer:: fin, ier
6 integer:: icount, istatus
7 integer:: imax, jmax
8 double precision, dimension(:,:), allocatable::grid_x, grid_y
9 character(200)::condFile
10
11 icount = nargs()
12 if ( icount.eq.2 ) then
13 call getarg(1, condFile, istatus)
14 else
15 stop "Input File not specified."
16 endif
17
18 ! 格子生成データファイルを開く
19 call cg_iric_open(condFile, IRIC_MODE_MODIFY, fin, ier)
20 if (ier /=0) stop "*** Open error of CGNS file ***"
21
22 imax = 10
23 jmax = 10
24
25 ! 格子生成用のメモリを確保
26 allocate(grid_x(imax,jmax), grid_y(imax,jmax)
27
28 ! 格子を生成
29 do i = 1, imax
30 do j = 1, jmax
31 grid_x(i, j) = i
32 grid_y(i, j) = j
33 end do
34 end do
35
36 ! 格子を出力
37 cg_iric_write_grid2d_coords(fin, imax, jmax, grid_x, grid_y, ier)
38
39 ! 格子生成データファイルを閉じる
40 call cg_iric_close(fin, ier)
41end program SampleProgram
コンパイルしたら、できた実行プログラムを フォルダの作成 で作成したフォルダにコピーし、名前を 基本情報の作成 で executable 属性に指定した名前 (この例なら "generator.exe" ) に変更してください。またこの時、格子生成プログラムの実行に必要な DLLなども同じフォルダにコピーしてください。
この段階で、iRIC から格子生成プログラムが正しく起動できるか確認します。
ソルバーに "Nays2DH" を指定して、新しいプロジェクトを開始し、 基本情報の作成 で行ったのと同じ操作で 格子生成アルゴリズムに "Sample Grid Creator" を選択し、格子生成ダイアログを表示します。表示されるダイアログを 図 36 に示します。

図 36 格子生成条件設定ダイアログ 表示例¶
"格子生成" ボタンを押します。すると、格子生成プログラムが 10 x 10 の 格子を生成し、それがiRIC 上に読み込まれるのが確認できます。 "格子生成" ボタンを押した後のプリプロセッサの表示画面を 図 37 に示します。

図 37 プリプロセッサ表示例¶
なお、この節で追加した格子出力用の関数の詳細については、 計算格子の出力 を参照してください。ただし、 計算格子の出力 では3次元格子の出力用関数についても 解説していますが、格子生成プログラムで利用できるのは、2次元格子の出力用関数だけです。