# 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-rec7.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-rec7.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

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

# 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'

# Initializing GOCAD file
# ~~~~~~~~~~~~~~~~~~~~~~~
# Form and name of temporary GOCAD file
  VRML='GOCAD'  WRL='qi-ray.mx'
  iniwrl:

# Adding model box wireframe
  MODEL='qi-mod.dat'  LIN='qi-bnd.out'
  KOLLIN=0  R=1.00  G=1.00  B=1.00  NAME='ModelBox'
  bndlin:
  linwrl:

# Adding source and ray endpoints into the VRML or GOCAD file
  PTS='qi-src.dat'
  KOLPTS=0  R=1.00  G=1.00  B=0.00  NAME='Source'
  ptswrl:
  PTS='qi-rec.dat'
  KOLPTS=0  R=1.00  G=1.00  B=0.00  NAME='Receivers'
  ptswrl:

# Temporary formatted file with rays
  RAYS='qi-ray.out'  LIN='qi-ray.out'

# Model QI
# ~~~~~~~~
  MODEL='qi-mod.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qi1-lsti.out'  PTS='ptsi.out'
  crt:
  crtpts:
  green:
  KOLLIN=0  R=1.00  G=0.00  B=0.00  NAME='Isotropic rays'
  crtray:
  linwrl:

  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qi1-lsta.out'  PTS='ptsa1.out'
  crt:
  crtpts:
  green:
  QILSTI='qi1-lsti.out'  QILSTA='qi1-lsta.out'
  watable:
  KOLLIN=0  R=0.00  G=1.00  B=1.00  NAME='QI rays'
  crtray:
  linwrl:

# Model QI2
# ~~~~~~~~~
  MODEL='qi-mod2.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qi2-lsti.out'
  crt:
  green:

  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qi2-lsta.out'  PTS='ptsa2.out'
  crt:
  crtpts:
  green:
  QILSTI='qi2-lsti.out'  QILSTA='qi2-lsta.out'
  watable:
  KOLLIN=0  R=0.00  G=1.00  B=1.00  NAME='QI2 rays'
  crtray:
  linwrl:

# Model QI4
# ~~~~~~~~~
  MODEL='qi-mod4.dat'
  CRTANI=0  RPAR='qi-rpa.dat'   QILST='qi4-lsti.out'
  crt:
  green:

  CRTANI=1  RPAR='qi-rpaa.dat'  QILST='qi4-lsta.out'  PTS='ptsa4.out'
  crt:
  crtpts:
  green:
  QILSTI='qi4-lsti.out'  QILSTA='qi4-lsta.out'
  watable:
  KOLLIN=0  R=0.00  G=1.00  B=1.00  NAME='QI4 rays'
  crtray:
  linwrl: