# Coupling ray theory seismograms in weakly anisotropic model QI
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
 #chk.pl: "data/qi/" "qi-mod.dat"
 #chk.pl: "data/qi/" "qi-src.dat"
 #chk.pl: "data/qi/" "qi-rec.dat"
 #chk.pl: "data/qi/" "qi-crt.dat"
 #chk.pl: "data/qi/" "qi-rpa.dat"
 #chk.pl: "data/qi/" "force.dat"
 #chk.pl: "crt/"     "writ.dat"

# Model QI for comparison with package ANRAY
  MODEL='qi-mod.dat'

# Source and receivers
  SRC='qi-src.dat'  REC='qi-rec.dat'  SOURCE='force.dat'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='qi-crt.dat'
  CODE='qi-crt.dat'
  RPAR='qi-rpa.dat'
  WRIT='writ.dat'       # only two-point rays stored in files
  CRTLOG='qi-crt.out'

# Kind of ray parameters (for details see 'init.for')
  INIPAR=3  # azimuthal equidistant projection

# Two-point ray tracing
  crt:

# Data describing the source time function
  KSIG=1 # (Gabor signal)
  SIGT=0.  SIGF=50.  SIGW=4.  SIGPH=0.  SIGA=100.

# Data describing the filtration of the source time function
  FMIN=0.0  FLOW=5.0  FHIGH=60.  FMAX=100.  # (cosine filter)

# Time step and time interval for the Fast Fourier Transform
  DT=0.000488  NFFT=2048

# Data to control optional plotting of program SS
  MPTS=200

# Coupling ray theory seismograms
  QICHM=0  QITT=0  QIDT=0
  SS='qi-ss.gse'                # GSE file with synthetic seismograms
  SIGPLOT='qi-sig.ps'           # Optional plot of program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Quasi-isotropic ray theory seismograms
  QICHM=1  QITT=1  QIDT=1
  SS='qi-ssqi.gse'              # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Data to control seismogram plotting (program SP)
  KODESP=3     SPCHRH=0.25
  SPTMIN=0.40  SPTMAX=0.60  SPTLEN=12.0  SPTDIV=2   SPTSUB=10
  SPXMIN=-.03  SPXMAX=0.60  SPXLEN=15.0  SPXDIV=7   SPXSUB=9
  NORMSP=1     SPAMP=0.80

# Plotting seismograms
  SS='qi-ssqi.gse'  KOLOR=2     # first plotting QI seismograms in red
  SS1='qi-ss.gse'   KOLOR1=1    # overlaying CRT seismograms in black
  SP1='qi-ss1.ps'  SP2='qi-ss2.ps'  SP3='qi-ss3.ps'
  sp:                           # seismogram plotting