# History file 'len-crt.h' to perform two-point ray tracing
# in the model with a LENticular inclusion
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
  chk.pl: "data/len/" "len-mod.dat"
  chk.pl: "data/len/" "len-rpa.dat"
  chk.pl: "data/len/" "len-rec.dat"
  chk.pl: "data/len/" "len-src.dat"
  chk.pl: "data/len/" "len-crt.dat"
  chk.pl: "crt/"      "codep.dat"
  chk.pl: "crt/"      "writall.dat"
  chk.pl: "forms/"    "hsv.dat"
# Output file of 'len-mod.h' used here as an input file
  chk.pl:  ""          "len-srf.wrl"  # VRML
# chk.pl:  ""          "len-srf.mx"   # GOCAD

# Model
  MODEL='len-mod.dat'

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

# Data filenames for CRT (for details see 'crt.for')
  DCRT='len-crt.dat'
  RPAR='len-rpa.dat'
  CODE='codep.dat'      # just refracted P wave
  WRIT='writall.dat'    # all rays stored in files
  CRTLOG='len-log.out'

# Kind of ray parameters (for details see 'init.for')
  INIPAR=2  # geographic-like spherical coordinates

# Ray tracing of the refracted wave
  crt:

# Generating the figures of the ray parameters
  REC='len-rec.dat' SYMBOLS='symbol.dat' COLORS='color.dat'
  ISREC=4 ICREC=1 HREC=0.4
    IRBAS=1 IRTWO=0 IRAUX=0 ITHOM=0 ISANG=1 ISUC=0
  RPPLOT='len-rp2.ps'   rpplot: 'len-crt.h' /
                    IRAUX=1                            HRAUX=0.08
  RPPLOT='len-rp4.ps'  rpplot: 'len-crt.h' /
    IRBAS=0         IRAUX=0 ITHOM=1
  RPPLOT='len-rp5.ps'  rpplot: 'len-crt.h' /
    IRBAS=1 IRTWO=1 IRAUX=1                 ISUC=1     HRTWO=0.5
  RPPLOT='len-rp6.ps'  rpplot: 'len-crt.h' /
                                    ISANG=0 ISUC=0     HRAUX=0.1
  PLIM1=0. PLIM2=60. PLIM3=-20. PLIM4=20.   SYMBOLS='symbol7.dat'
  RPPLOT='len-rp7.ps'  rpplot: 'len-crt.h' /

# Supplementing receivers with the third coordinate
  INTF='len-rec.dat' INTFOUT='len-rec.out' KSRFC=2 KOLUMN=3
  intf:
# Note  File 'len-rec.dat' contains only X1,X2 receiver coordinates.
#       If it were submitted to program CRTPTS as the receiver
#       file, the travel time could be extrapolated in the horizontal
#       plane to the receiver position, but the vertical X3 coordinate
#       would correspond to the numerical approximation of the point
#       of intersection with the reference surface (Earth surface)
#       rather than to the exact X3 coordinate of the receiver.

# Rewriting the travel times at the receivers into a formatted file
  REC='len-rec.out' SRC= PTS='len-ttr.out'  NQ=4 KALL=0 KREC=1 KTT=0
  crtpts:

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

# Adding two-point rays into the VRML or GOCAD file with surfaces
  VRML='VRML2'  WRL='len-srf.wrl'  WRLOUT='len-ray.wrl'  # VRML
# VRML='GOCAD'  WRL='len-srf.mx'   WRLOUT='len-ray.mx'   # GOCAD
  LIN='len-ray.out'
# Setting blue colour for rays
  KOLLIN=0  R=0.00  G=0.00  B=1.00
  NAME='Rays'  PROPERTIES=' '  # just for GOCAD
# Example of parameters for colouring rays according to travel time
# KOLLIN=4  COLORS='hsv.dat'  VPER=6.0  VREF=0.0  CREF=0.166667
# NAME='Rays'  PROPERTIES='T'  # just for GOCAD
  linwrl:

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

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

# Temporary files (may be deleted)
# '*.tmp'

# List of output PostScript figures
# len-rp2.ps, len-rp4.ps, ..., len-rp7.ps... Refer to paper
#     Bulant, P.  3-D two-point ray tracing in pictures.  In Seismic
#             Waves in Complex 3-D Structures, Report 6 (Department of
#             Geophysics, Charles University, Prague, September 1997)

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