# Two-point ray tracing in model EX
# followed by plotting of the computed two-point rays
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file is designed to be run from another history file
# (see e.g. 'ex.h'), which runs 'srp.for' to generate source
# coordinates 'ex-src.out' and receiver coordinates 'ex-rec.out'
# and then executes this history file for the source.

# Input files required
  chk.pl: "data/ex/" "ex-mod.dat"
  chk.pl: "data/ex/" "ex-dcr.dat"
  chk.pl: "data/ex/" "ex-rpa.dat"
  chk.pl: "data/ex/" "ex-cod.dat"
  chk.pl: "crt/"      "writ.dat"
  chk.pl: ""          "ex-rec.out"
  chk.pl: ""          "ex-src.out"
# Output of 'ex-srf.h'
  chk.pl: ""          "ex-srf.wrl"
 #chk.pl: ""          "ex-srf.mx"   # Gocad

# Model
  MODEL='ex-mod.dat'

# Two-point ray tracing
# ~~~~~~~~~~~~~~~~~~~~~

# Source and receivers
  SRC='ex-src.out'  REC='ex-rec.out'

# Data filenames for CRT
  DCRT='ex-dcr.dat'
  RPAR='ex-rpa.dat'
  CODE='ex-cod.dat'
  WRIT='writ.dat'
  CRTLOG='ex-log.out'

# Kind of ray parameters
  INIPAR=3

# Screen output
  CRTSCRO='SWRG'

# Ray tracing of the reflected wave
  crt:

# Rewriting the travel times at the receivers into a formatted file
  PTS='ex-tim.out'  NQ=4
  crtpts:

# Rewriting the two-point rays into a formatted file
  RAYS='ex-ray.out'
  crtray:

# Displaying the rays in VRML or Gocad
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Form and name of input and output files
  VRML='VRML2'  WRL='ex-srf.wrl'  WRLOUT='ex-ray.wrl'  # VRML
 #VRML='GOCAD'  WRL='ex-srf.mx'   WRLOUT='ex-ray.mx'   # GOCAD

  LIN='ex-ray.out'
# Setting red colour for rays
  KOLLIN=0  R=1.00  G=0.00  B=0.00
  NAME='Rays'  PROPERTIES=' '  # just for GOCAD
  linwrl:

  VRML='VRML2'  WRL='ex-ray.wrl'  # VRML
 #VRML='GOCAD'  WRL='ex-ray.mx'   # GOCAD

# Adding source and receivers into the VRML or GOCAD file
  PTS='ex-src.out'
  KOLPTS=0  R=1.00  G=1.00  B=0.00
  NAME='Source'  PROPERTIES=' '  # just for GOCAD
  ptswrl:
  PTS='ex-rec.out'
  KOLPTS=0  R=0.00  G=0.00  B=1.00
  NAME='Receivers'  PROPERTIES=' '  # just for GOCAD
  ptswrl: