格子の出力処理の記述

格子の出力処理を記述します。

まずは、iRIC との連携が正しく行えることを確認するため、 単純な格子を生成して出力する処理を記述します。

格子を出力する処理を追記したソースコードを リスト 20 に示します。追記した部分を強調して示します。

リスト 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 に示します。

../_images/gridgen_cond_dialog_for_testing.png

図 36 格子生成条件設定ダイアログ 表示例

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

../_images/preprocessor_after_gridgen.png

図 37 プリプロセッサ表示例

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