9.GMS使用イメージ--時間積分 後処理

!境界条件を満たすようにのりしろ部分をそろえる
  call cyclic_boundary_1(u)
  call cyclic_boundary_1(h)

!境界条件を満たすようにのりしろ部分をそろえる
     call cyclic_boundary_1(u_a)
     call cyclic_boundary_1(h_a)

!--Asselin(1972)のタイムフィルタ
     u = u + 0.1D0 * (u_a - 2.0D0 * u + u_b) 
     h = h + 0.1D0 * (h_a - 2.0D0 * h + h_b) 

!次のstepの積分のために配列のコピー
     u_b = u
     h_b = h
     h   = h_a
     u   = u_a
  end do
end program shallow_1d