mod_auxil.f90

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

!_______________________________________________________________________
 
MODULE AUXIL
  USE PRECISION, ONLY: WP

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

  IMPLICIT NONE

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

  INTEGER :: NSTFS
  
  REAL (KIND=WP), PARAMETER :: PI=3.141592653589793,                   &
                               A = -1./24., B = 9./8.

  REAL (KIND=WP)::  AH, BH, CH

  REAL (KIND=WP), DIMENSION (:), ALLOCATABLE ::  SOURTF

  REAL (KIND=WP), DIMENSION (:), ALLOCATABLE ::  OTA, OTB, G1, G2

END MODULE AUXIL