7.GMSの使用イメージ --GMSの初期化

GMSではグリッド情報を下図のように定義している。

program shallow_1d
  use GMS
  implicit none
!物理境界とのりしろの大きさ
  integer, parameter :: plb = 0, pub=100
  integer, parameter :: margin = 1

!格子間隔
  real(8), parameter :: delta_x = 1.0D0
!grid情報 定義
  integer, dimension(1) :: u_grid = (/0/)
  integer, dimension(1) :: h_grid = (/1/)

!GMS変数 宣言
  type(var1d) :: u, u_a, u_b
  type(var1d) :: h, h_a, h_b


!GMSに必要な情報を教える
  call set_physical_bound_1(plb, pub)
  call set_array_margin_1(margin)
  call gms_set_interval_x(delta_x)

!--work領域の大きさ
  call init_memory_manager_1d(20)

!--GMS変数にグリッド情報を教える
  call def_var1d(u,   u_grid)
  call def_var1d(u_a, u_grid)
  call def_var1d(u_b, u_grid)

  call def_var1d(h,   h_grid)
  call def_var1d(h_a, h_grid)
  call def_var1d(h_b, h_grid)