# Coupling ray theory seismograms in weakly anisotropic model QI
# Second-order perturbations of travel times
# Anisotropic common ray tracing
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file may be run directly for one of the models
# 'qi-mod.dat', 'qi-modh.dat', 'qi-mod2.dat', 'qi-mod4.dat',
# 'qi-mod8.dat', or may be called from the history file 'qi-ss4m.h'
#  for the model 'qi-mod.tmp' - edit the line containing "MODEL=" below.

# Executables required
# ~~~~~~~~~~~~~~~~~~~~
# This history file executes compiled programs 'greentc.for' and
# 'tccomp.for'.  Programs 'greentc.for' and 'tccomp.for' are located
# in directory 'data/tc'.

# Input files required
  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"

# Model QI, QIH, QI2, QI4 or QI8
# MODEL='qi-mod.dat'   chk.pl: "data/qi/" "qi-mod.dat"
# MODEL='qi-modh.dat'  chk.pl: "data/qi/" "qi-modh.dat"
# MODEL='qi-mod2.dat'  chk.pl: "data/qi/" "qi-mod2.dat"
# MODEL='qi-mod4.dat'  chk.pl: "data/qi/" "qi-mod4.dat"
# MODEL='qi-mod8.dat'  chk.pl: "data/qi/" "qi-mod8.dat"
  MODEL='qi-mod.tmp'   chk.pl: " "        "qi-mod.tmp"

# 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'
  WRIT='writ.dat'       # only two-point rays stored in files

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

# Parameter to get the same results of ray tracing as in version 5.70
  UEBMUL=1.

# Degree of the homogeneous Hamiltonian for anisotropic ray tracing
  DEGREE=-1

# 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

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

# Data to control seismogram plotting (program SP)
  KODESP=3     SPCHRH=0.25
  SPTMIN=0.40  SPTMAX=0.55  SPTLEN=6.0   SPTDIV=3   SPTSUB=5
  SPXMIN=-.03  SPXMAX=0.60  SPXLEN=15.0  SPXDIV=7   SPXSUB=9
  NORMSP=1     SPAMP=0.80
  SS1='qi-ss.gse'   KOLOR=2 KOLOR1=1

# Isotropic two-point ray tracing
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTANI=0  RPAR='qi-rpa.dat'  CRTLOG='qi-crt.out'
  crt:

# Coupling ray theory seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=0  QITT=0  QIDT=0  QIPV=0  QIRAY=0  QILST='qilst.out'
  GREEN='green.out'
  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
  GREENTC='greentc.out'
  greentc:

# Quasi-isotropic ray theory seismograms - Christoffel matrix
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=1  QITT=0  QIDT=0  QIPV=0  QIRAY=0  QILST=' '
  GREEN='greenc.out'
  SS='qi-ssc.gse'               # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Plotting seismograms
  SP1='qi-ss1c.ps'  SP2='qi-ss2c.ps'  SP3='qi-ss3c.ps'
  sp:

# Comparison of plane-wave propagator matrices
  GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errc.out'
  INORM=2
  greentc:
  tccomp:

# Quasi-isotropic ray theory seismograms - travel times
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=0  QITT=1  QIDT=1  QIPV=0  QIRAY=0
  GREEN='greent.out'
  SS='qi-sst.gse'               # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Plotting seismograms
  SP1='qi-ss1t.ps'  SP2='qi-ss2t.ps'  SP3='qi-ss3t.ps'
  sp:

# Comparison of plane-wave propagator matrices
  GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errt.out'
  INORM=2
  greentc:
  tccomp:

# Quasi-isotropic ray theory seismograms - polarization vector
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=0  QITT=0  QIDT=0  QIPV=1  QIRAY=0
  GREEN='greenp.out'
  SS='qi-ssp.gse'               # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Plotting seismograms
  SP1='qi-ss1p.ps'  SP2='qi-ss2p.ps'  SP3='qi-ss3p.ps'
  sp:

# Comparison of plane-wave propagator matrices
  GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errp.out'
  INORM=2
  greentc:
  tccomp:

# Quasi-isotropic ray theory seismograms - anisotropic common ray
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=0  QITT=0  QIDT=0  QIPV=0  QIRAY=1
  GREEN='greenr.out'
  SS='qi-ssr.gse'               # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Plotting seismograms
  SS='qi-ss.gse'  SS1='qi-ssr.gse'
  SP1='qi-ss1r.ps'  SP2='qi-ss2r.ps'  SP3='qi-ss3r.ps'
  sp:

# Comparison of plane-wave propagator matrices
  GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-errr.out'
  INORM=2
  greentc:
  tccomp:

# Plotting seismograms of 3 different QI approximations together
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  SS='qi-ssc.gse'
  SS1='qi-ssp.gse'
  SS2='qi-sst.gse'
  SS3='qi-ss.gse'
  KOLOR=2  KOLOR1=4   KOLOR2=3  KOLOR3=1
  SP1='qi-ss1.ps'  SP2='qi-ss2.ps'  SP3='qi-ss3.ps'
  sp:

# Anisotropic common ray two-point ray tracing
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTANI=1  RPAR='qi-rpaa.dat'  CRTLOG='qi-crta.out'
  crt:

# Coupling ray theory seismograms along the common anisotropic rays
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  QICHM=0  QITT=0  QIDT=0  QIPV=0  QIRAY=0  QILST='qilsta.out'
  GREEN='green.out'
  SS='qi-ssa.gse'               # GSE file with synthetic seismograms
  SIGPLOT=                      # no plotting by program SS
  green:                        # Green function
  greenss:                      # response function
  ss:                           # synthetic seismograms

# Plotting seismograms
  SS='qi-ss.gse'  SS1='qi-ssa.gse'  SS2=  SS3=
  KOLOR=2  KOLOR1=1   KOLOR2=   KOLOR3=
  SP1='qi-ss1a.ps'  SP2='qi-ss2a.ps'  SP3='qi-ss3a.ps'
  sp:
  SS='qi-ssr.gse'  SS1='qi-ssa.gse'
  SP1='qi-ss1q.ps'  SP2='qi-ss2q.ps'  SP3='qi-ss3q.ps'
  sp:

# Comparison of plane-wave propagator matrices
  GREENTC='greentcq.out' GREENTC2='greentc.out' TCCOMP='qi-erra.out'
  INORM=2
  greentc:
  tccomp:

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# 'qi-ss1?.ps','qi-ss2?.ps','qi-ss3?.ps' ... Synthetic seismograms
#             compared with coupling ray theory seismograms.
# 'qi-err?.out' ... Norm of the difference from the coupling ray theory
#             ray propagator matrix.
#     where ? identifies either the quasi-isotropic approximation being
#     compared
#             c ... Christoffel matrix
#             t ... travel times
#             p ... polarization vector
#             r ... anisotropic common ray (perturbation)
#               ... Christoffel matrix & travel times & polariz. vector
#     or the calculation along anisotropic common rays
#             a ... anisotropic common ray (real anisotropic common ray)
#             q ... anisotropic common ray real & perturbation
# 'qilst.out' ... Second-order perturbations of travel times.
# 'qilsta.out'... Second-order perturbations of travel times calculated
#                 along anisotropic common rays.