# Test of anisotropic common S wave ray tracing in models QI, QI2, QI4
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
 #chk.pl: "data/qi/" "qi-mod.dat"
 #chk.pl: "data/qi/" "qi-mod2.dat"
 #chk.pl: "data/qi/" "qi-mod4.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"

# 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
  CRTLOG='qi-crt.out'

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

# Degree of averaged homogeneous S-wave Hamiltonians for ray tracing
  DEGREE=-1

# Frequency domain and other data for program GREEN
  OF=0.0  NF=101  DF=1.000576332
  QICHM=0  QITT=0  QIDT=0  QIPV=0  QIRAY=0
  GREEN='green.out'

# Model QI
# ~~~~~~~~
  MODEL='qi-mod.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qilsti.out'  PTS='ptsi1.out'
  crt:
  crtpts:
  green:
  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qilsta.out'  PTS='ptsa1.out'
  crt:
  crtpts:
  green:
  QILSTI='qilsti.out'  QILSTA='qilsta.out'  QILSTC='qilst1.out'
  qilst:

# Model QI2
# ~~~~~~~~~
  MODEL='qi-mod2.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qilsti.out'  PTS='ptsi2.out'
  crt:
  crtpts:
  green:
  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qilsta.out'  PTS='ptsa2.out'
  crt:
  crtpts:
  green:
  QILSTI='qilsti.out'  QILSTA='qilsta.out'  QILSTC='qilst2.out'
  qilst:

# Model QI4
# ~~~~~~~~~
  MODEL='qi-mod4.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qilsti.out'  PTS='ptsi4.out'
  crt:
  crtpts:
  green:
  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qilsta.out'  PTS='ptsa4.out'
  crt:
  crtpts:
  green:
  QILSTI='qilsti.out'  QILSTA='qilsta.out'  QILSTC='qilst4.out'
  qilst: