8.GMS使用イメージ--時間積分

!初期値を代入
  call input_var(u0, u_b) !u0は普通の配列
  call input_var(h0, h_b) !h0は普通の配列

・・・!u, hの計算など

   do i = 1, nloop
!leap-frog法で時間積分する
     h_a = h_b - depth * d_x(u) * 2.0D0 * dt
     u_a = u_b -  grav * d_x(h) * 2.0D0 * dt 
!従来の書き方だと
!h_a(1:imax) = h_b(1:imax) &
!     -depth*(( u(1:imax) - u(0:imax-1) ) / dx) * 2.0D0 * dt
!u_a(0:imax)=u_b(0:imax) &
!     -grav * (( h(1:imax+1) - h(0:imax) ) / dx) * 2.0D0 * dt

!配列添字の追放で可読性が大幅に向上