rec_dis.f90

1D有限差分波动方程模拟

!=======================================================================

SUBROUTINE  REC_DIS ( ITILE, D )

  USE PRECISION   , ONLY:                                              &
                          WP     
  USE CONTROL_DATA, ONLY:                                              &
                          DT, LREC, MR, MZ

!-----------------------------------------------------------------------

  IMPLICIT NONE

  INTEGER,  INTENT(IN) :: ITILE
  REAL(WP), DIMENSION(0:MZ), INTENT(INOUT) :: D

  INTEGER :: I

!-----------------------------------------------------------------------

  WRITE ( 21, '(250(1X,E12.6))' ) ITILE*DT, ( D( LREC(I) ), I = 1, MR )

END  SUBROUTINE  REC_DIS