時間積分・ファイル出力

do i = 0, nloop

  time = i * dt

  u_a = u_b -  grav * d_x(h) * 2.0D0 * dt !Leap-Frog法
...
  call cyclic_boundary_x(u_a)
  u = u + c_asselin * (u_a - 2.0D0 * u + u_b) !Asselin time filter

...
!---file出力
     call HistoryPut ("t", real(time)  , hst_u)
     call HistoryPut ("u", get(u), hst_u)
...

  u_b = u
  u = u_a
...
end do