座標軸設定用モジュール
real(8) :: Xmin, Xmax ! x 座標の始点・終点 real(8) :: Zmin, Zmax ! z 座標の始点・終点 real(8) :: DelX, DelZ !格子間隔 real(8) :: DelXZ !平均格子間隔(平均混合距離) integer :: NX, NZ !格子点数 integer :: DimXMin ! x 方向の配列の下限 integer :: DimXMax ! x 方向の配列の上限 integer :: DimZMin ! z 方向の配列の下限 integer :: DimZMax ! z 方向の配列の上限 integer :: RegXMin ! x 方向の物理領域の下限 integer :: RegXMax ! x 方向の物理領域の上限 integer :: RegZMin ! z 方向の物理領域の下限 integer :: RegZMax ! z 方向の物理領域の上限 integer, parameter :: Margin = 5 !境界のグリッド数 integer, parameter :: kind = 8 !精度を表す real(8), allocatable :: s_X(:) !X 座標軸(スカラー) real(8), allocatable :: f_X(:) !X 座標軸(ベクトル) real(8), allocatable :: s_Z(:) !Z 座標軸(スカラー) real(8), allocatable :: f_Z(:) !Z 座標軸(ベクトル)
NAMELIST から情報を得て, 格子点を計算する
subroutine gridset_init(cfgfile)
use dc_trace, only: BeginSub, EndSub
character(*), intent(in) :: cfgfile
NAMELIST /gridset/ & & NX, NZ, Xmin, Xmax, Zmin, Zmax