main.f90

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

PROGRAM P_1DFD_DS

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

  IMPLICIT NONE

  CHARACTER(LEN=17) :: JOB_NAME

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

  OPEN  (10, FILE = 'HF_1DFD_DS', STATUS='OLD')
  READ  (10,'(A)') JOB_NAME
  CLOSE (10)

  CALL  READ_OPEN ( JOB_NAME (1:LEN_TRIM(JOB_NAME)) )
                         print *,  'read_open'
  CALL  AUX1
                         print *,  'aux1'
  CALL  ALLOC_NR
                         print *,  'alloc_nr'
  CALL  ALLOC_WAVEFIELD
                         print *,  'alloc_wavefield'
  CALL  READ_MODEL_PAR
                         print *,  'read_model_par'
  CALL  OPEN_SR_FILE
                         print *,  'open_sr_file'
  CALL  PRE_NR_COEF1
                         print *,  'pre_nr_coef1'
  CALL  TIME_LOOP


END PROGRAM P_1DFD_DS