# History file 'l7-crt.h' to perform two-point ray tracing
# in the model L7
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Under MS-DOS Lahey Fortran use the option "assume infinite NDP stack"
#
# Input files required
  chk.pl: "data/l7/" "l7-mod.dat"
  #chk.pl: "data/l7/" "l7s-mod.dat"  # alternative file with the model
  chk.pl: "data/l7/" "l7-rpa.dat"
  chk.pl: "data/l7/" "l7-rec.dat"
  chk.pl: "data/l7/" "l7-src.dat"
  chk.pl: "data/l7/" "l7-dcr.dat"
  chk.pl: "data/l7/" "l7-cod.dat"
  chk.pl: "crt/"      "writ.dat"
  chk.pl: "forms/"    "hsv.dat"
# Output file of 'l7-mod.h' used here as an input file
  chk.pl:  ""          "l7-srf.wrl"  # VRML
  #chk.pl: ""          "l7-srf.mx"   # GOCAD

# Model
  MODEL='l7-mod.dat'
 #MODEL='l7s-mod.dat' # alternative file with the model

# Source and receivers
  SRC='l7-src.dat'  REC='l7-rec.dat'

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

# Kind of ray parameters
  INIPAR=3

# Screen output
  CRTSCRO='SWR'

# Ray tracing of the reflected wave
# (24757 rays for l7-mod.dat and 45600 rays for l7s-mod.dat)
  crt:

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

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

# Adding two-point rays into the VRML or GOCAD file with surfaces
  VRML='VRML2'  WRL='l7-srf.wrl'  WRLOUT='l7-ray.wrl'  # VRML
 #VRML='GOCAD'  WRL='l7-srf.mx'   WRLOUT='l7-ray.mx'   # GOCAD
  LIN='l7-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:

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

# ======================================================================

# Output VRML or GOCAD file containing 3-D sculptures of rays
# and interfaces
# l7-ray.wrl or l7-ray.mx