# History file 'qih-ss.h' for computation of anisotropic, isotropic, and
# coupling-ray theory S-wave synthetic seismograms in model QIH.

# Executables required
# ~~~~~~~~~~~~~~~~~~~~
# This history file requires programs FRESAN, SYNFAN and ANRAYGSE
# of package ANRAY, and modified version of program ANRAY, please link
# program 'anray_44.for' located in directory 'data/qi'.

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: "data/qi/" "anrayh.qi"
  chk.pl: "data/qi/" "fresan.anv"
  chk.pl: "data/qi/" "fresan.ant"
  chk.pl: "data/qi/" "fresan.anr"
  chk.pl: "data/qi/" "synfan.an"
  chk.pl: "data/qi/" "qi-modh.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/" "qi-rpaa.dat"
  chk.pl: "data/qi/" "force.dat"
  chk.pl: "crt/"     "writ.dat"

# Data to control seismogram plotting (program SP)
  KODESP=3     SPCHRH=0.25
  SPTMIN=0.40  SPTMAX=0.50  SPTLEN=6.0   SPTDIV=1   SPTSUB=10
  SPXMIN=-.03  SPXMAX=0.60  SPXLEN=15.0  SPXDIV=7   SPXSUB=9
  NORMSP=1     SPAMP=0.80
  GSEWIDTH=8 (writing GSE files in greenss.for)

# Anisotropic ray theory synthetic seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  anray_44:   'anrayh.qi' 1* 'lu1.out' 'lu2.out' /

# Vertical component
  fresan:  'fresan.anv' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.an'  1* 'lu7.out' 'lu8.out' /
  LU8='lu8.out'
  SS='qih-ss3a.gse'
  anraygse:
  SP1=' '           SP2=' '           SP3='qih-ss3a.ps'
  sp:

# Transverse component
  fresan:  'fresan.ant' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.an'  1* 'lu7.out' 'lu8.out' /
  SS='qih-ss2a.gse'
  anraygse:
  SP1=' '           SP2='qih-ss2a.ps' SP3=' '
  sp:

# Radial component
  fresan:  'fresan.anr' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.an'  1* 'lu7.out' 'lu8.out' /
  SS='qih-ss1a.gse'
  anraygse:
  SP1='qih-ss1a.ps' SP2=' '           SP3=' '
  sp:

# Isotropic ray theory synthetic seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  SRC='qi-src.dat'      REC='qi-rec.dat'      SOURCE='force.dat'
  DCRT='qi-crt.dat'     CODE='qi-crt.dat'     WRIT='writ.dat'
  INIPAR=3

# 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 and frequency step and interval for the Fast Fourier Transform
  DT=0.000488  NFFT=2048
  OF=0.0    NF=101
  DF=1.000576332

# Ray tracing and synthetic seismograms
  MODEL='qi-modh.dat' CRTANI=0  RPAR='qi-rpa.dat'  CRTLOG='qi-crt.out'
  crt:
  MODEL=' '  GREEN='green.out'  SS='qih-ssi.gse'
  green:
  greenss:
  ss:
  SP1='qih-ss1i.ps'  SP2='qih-ss2i.ps'  SP3='qih-ss3i.ps'
  sp:

# Coupling ray theory synthetic seismograms - anisotropic common ray
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='qi-modh.dat' CRTANI=1  RPAR='qi-rpaa.dat'  CRTLOG='qi-crta.out'
  SS='qih-ssc.gse'
  crt:
  green:
  greenss:
  ss:
  SP1='qih-ss1c.ps'  SP2='qih-ss2c.ps'  SP3='qih-ss3c.ps'
  sp:

# Coupling ray theory synthetic seismograms - anisotropic common ray,
# second-order perturbations of travel time
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='qi-modh.dat' CRTANI=1 RPAR='qi-rpaa.dat'  CRTLOG='qi-crtas.out'
  QIRAY=1
  SS='qih-sscs.gse'
  crt:
  green:
  greenss:
  ss:
  SP1='qih-ss1s.ps'  SP2='qih-ss2s.ps'  SP3='qih-ss3s.ps'
  sp:
# Prevailing-frequency approximation
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTLOG='qi-crtaf.out'
  SS='qih-sscf.gse'
  crt:
  SINGLF=50.
  green:
  SINGLF=
  greenss:
  ss:
  SP1='qih-ss1f.ps'  SP2='qih-ss2f.ps'  SP3='qih-ss3f.ps'
  sp:

# Coupling, isotropic and anisotropic seismograms in a single plot
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  SS='qih-ssc.gse'   SRC=   REC=
  SS2='qih-ssi.gse'
  SS1='qih-ss1a.gse' SP1='qih-ss1.ps'  SP2=' '           SP3=' '
  sp:
  SS1='qih-ss2a.gse' SP1=' '           SP2='qih-ss2.ps'  SP3=' '
  sp:
  SS1='qih-ss3a.gse' SP1=' '           SP2=' '           SP3='qih-ss3.ps'
  sp:
  SS=    SS1=   SS2=

# Both coupling ray theory seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Data to control seismogram plotting (program SP)
  SPXMIN=0.0   SPXMAX=0.60  SPXLEN=15.0  SPXDIV=6   SPXSUB=5
  SPTTEXT='Time (s)'  SPXTEXT='Depth (km)'
  NORMSP=1     SPAMP=1.
  CALCOPS='0.50 setlinewidth'
  KOLOR=2 red   KOLOR1=1 black
  SS='qih-sscf.gse'   SRC=   REC=
  SS1='qih-sscs.gse'
  SP1='qih-sf1.ps'   SP2='qih-sf2.ps'   SP3='qih-sf3.ps'
  sp:
  KOLOR= KOLOR1=
  SS=    SS1=
  SPAMP= SPAMP1=

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# 'qih-ss1?.ps','qih-ss2?.ps','qih-ss3?.ps' ... Synthetic seismograms
#       where ? identifies the calculation
#             a ... anisotropic ray theory
#             i ... isotropic ray theory
#             c ... coupling ray theory along anisotropic common ray
#             f ... prevailing-frequency approximation of the coupling
#                   ray theory
# 'qih-sf1.ps','qih-sf2.ps','qih-sf3.ps' ... Prevailing-frequency
#       approximation seismograms compared with standard coupling ray
#       theory seismograms.
# 'qih-ss1.ps','qih-ss2.ps','qih-ss3.ps' ... Coupling, isotropic and
#                   anisotropic seismograms plotted together.