計算条件、計算格子、境界条件の読み込み処理の記述

計算条件、計算格子、境界条件の読み込み処理を記述します。

iRIC は、 ソルバー定義ファイルの作成 で作成したソルバー定義ファイルに従って、 計算条件、格子、格子属性、境界条件を計算データファイルに出力しますので、 ソルバー定義ファイルでの記述に対応するように、計算条件、計算格子、境界条件の 読み込み処理を記述します。

計算条件、計算格子の読み込み処理を追記したソースコードを リスト 8 に示します。強調して示したのが追記した部分です。

リスト 8 計算データファイルを開く処理、閉じる処理を追記したソースコード
  1program SampleProgram
  2  implicit none
  3  include 'cgnslib_f.h'
  4  include 'iriclib_f.h'
  5  integer:: fin, ier
  6  integer:: icount, istatus
  7  character(200)::condFile
  8  integer:: maxiterations
  9  double precision:: timestep
 10  integer:: surfacetype
 11  double precision:: constantsurface
 12  integer:: variable_surface_size
 13  double precision, dimension(:), allocatable:: variable_surface_time
 14  double precision, dimension(:), allocatable:: variable_surface_elevation
 15
 16  integer:: isize, jsize
 17  double precision, dimension(:,:), allocatable:: grid_x, grid_y
 18  double precision, dimension(:,:), allocatable:: elevation
 19  integer, dimension(:,:), allocatable:: obstacle
 20
 21  integer:: inflowid
 22  integer:: inflow_count
 23  integer:: inflow_element_max
 24  integer:: discharge_variable_sizemax
 25  integer, dimension(:), allocatable:: inflow_element_count
 26  integer, dimension(:,:,:), allocatable:: inflow_element
 27  integer, dimension(:), allocatable:: discharge_type
 28  double precision, dimension(:), allocatable:: discharge_constant
 29  integer, dimension(:), allocatable:: discharge_variable_size
 30  double precision, dimension(:,:), allocatable:: discharge_variable_time
 31  double precision, dimension(:,:), allocatable:: discharge_variable_value
 32
 33  write(*,*) "Sample Program"
 34
 35  ! (略)
 36
 37  ! 内部変数の初期化
 38  call cg_iric_init_f(fin, ier)
 39  if (ier /=0) STOP "*** Initialize error of CGNS file ***"
 40  ! オプションの設定
 41  call iric_initoption_f(IRIC_OPTION_CANCEL, ier)
 42  if (ier /=0) STOP "*** Initialize option error***"
 43
 44  ! 計算条件の読み込み
 45  call cg_iric_read_integer_f("maxIteretions", maxiterations, ier)
 46  call cg_iric_read_real_f("timeStep", timestep, ier)
 47  call cg_iric_read_integer_f("surfaceType", surfacetype, ier)
 48  call cg_iric_read_real_f("constantSurface", constantsurface, ier)
 49
 50  call cg_iric_read_functionalsize_f("variableSurface", variable_surface_size, ier)
 51  allocate(variable_surface_time(variable_surface_size))
 52  allocate(variable_surface_elevation(variable_surface_size))
 53  call cg_iric_read_functional_f("variableSurface", variable_surface_time, variable_surface_elevation, ier)
 54
 55  ! 格子のサイズを調べる
 56  call cg_iric_gotogridcoord2d_f(isize, jsize, ier)
 57
 58  ! 格子を読み込むためのメモリを確保
 59  allocate(grid_x(isize,jsize), grid_y(isize,jsize))
 60  ! 格子を読み込む
 61  call cg_iric_getgridcoord2d_f(grid_x, grid_y, ier)
 62
 63  ! 格子点で定義された属性 のメモリを確保
 64  allocate(elevation(isize, jsize))
 65  allocate(obstacle(isize - 1, jsize - 1))
 66
 67  ! 属性を読み込む
 68  call cg_iric_read_grid_real_node_f("Elevation", elevation, ier)
 69  call cg_iric_read_grid_integer_cell_f("Obstacle", obstacle, ier)
 70
 71  ! 流入口の数に従って、境界条件を保持するメモリを確保。
 72  allocate(inflow_element_count(inflow_count))
 73  allocate(discharge_type(inflow_count), discharge_constant(inflow_count))
 74  allocate(discharge_variable_size(inflow_count))
 75
 76  ! 流入口に指定された格子点の数と、時間依存の流入量のサイズを調べる
 77  inflow_element_max = 0
 78  do inflowid = 1, inflow_count
 79    ! 流入口に指定された格子点の数
 80    call cg_iric_read_bc_indicessize_f('inflow', inflowid, inflow_element_count(inflowid))
 81    if (inflow_element_max < inflow_element_count(inflowid)) then
 82      inflow_element_max = inflow_element_count(inflowid)
 83    end if
 84    ! 流入口の時間依存の流入量のデータの数
 85    call cg_iric_read_bc_functionalsize_f('inflow', inflowid, 'FunctionalDischarge', discharge_variable_size(inflowid), ier);
 86    if (discharge_variable_sizemax < discharge_variable_size(inflowid)) then
 87      discharge_variable_sizemax = discharge_variable_size(inflowid)
 88    end if
 89  end do
 90
 91  ! 流入口に指定された格子点と、時間依存の流入量を保持するメモリを確保。
 92  allocate(inflow_element(inflow_count, 2, inflow_element_max))
 93  allocate(discharge_variable_time(inflow_count, discharge_variable_sizemax))
 94  allocate(discharge_variable_value(inflow_count, discharge_variable_sizemax))
 95
 96  ! 境界条件の読み込み
 97  do inflowid = 1, inflow_count
 98    ! 流入口に指定された格子点
 99    call cg_iric_read_bc_indices_f('inflow', inflowid, inflow_element(inflowid:inflowid,:,:), ier)
100    ! 流入量の種類 (0 = 一定、1 = 時間依存)
101    call cg_iric_read_bc_integer_f('inflow', inflowid, 'Type', discharge_type(inflowid:inflowid), ier)
102    ! 流入量 (一定)
103    call cg_iric_read_bc_real_f('inflow', inflowid, 'ConstantDischarge', discharge_constant(inflowid:inflowid), ier)
104    ! 流入量 (時間依存)
105    call cg_iric_read_bc_functional_f('inflow', inflowid, 'FunctionalDischarge', discharge_variable_time(inflowid:inflowid,:), discharge_variable_value(inflowid:inflowid,:), ier)
106  end do
107
108  ! 計算データファイルを閉じる
109  call cg_close_f(fin, ier)
110  stop
111end program SampleProgram

計算条件などを読み込む関数に渡す引数が、 計算条件の定義, 格子属性の定義 でソルバー定義ファイルに定義した Item 要素の name 属性と一致していることに注目してください。

なお、ソルバー定義ファイルで定義する計算条件、格子、格子属性と、それを読み込むための iRIClib の関数の対応関係については、 計算条件・境界条件・格子生成条件の項目の定義と読み込み処理の例 を参照してください。

また、計算条件、計算格子、境界条件の読み込みに使う関数の詳細については、 計算条件 (もしくは格子生成条件) の読み込み, 計算格子の読み込み, 境界条件の読み込み を参照してください。