# Model with a LENticular inclusion, boundary-value rays corresponding
# to the plane wave incident at the Bottom of the model
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
 #chk.pl: "data/len/" "len-mod.dat"
 #chk.pl: "data/len/" "len-rec.dat"
 #chk.pl: "data/len/" "lenb-rpa.dat"
 #chk.pl: "data/len/" "len-crt.dat"
 #chk.pl: "crt/"      "codep.dat"
 #chk.pl: "crt/"      "writ.dat"

# Model and receivers
  MODEL='len-mod.dat'
  REC='len-rec.dat'

# Initial conditions for rays (for details see 'init.for')
  INIDIM=2             # initial surface
  INIPAR=0             # ray parameters = Cartesian coordinates
  INIT='lenb-rpa.dat'  # plane wave incident at the bottom of the model
  RPAR='lenb-rpa.dat'  # data for two-point ray tracing in the same file

# Other data filenames for CRT (for details see 'crt.for')
  DCRT='len-crt.dat'
  CODE='codep.dat'     # just refracted P wave
  WRIT='writ.dat'      # only two-point rays stored in files
  CRTLOG='lenb-log.out'

# Ray tracing of the refracted wave
  crt: 'lenb-crt.h' /

# Supplementing receivers with the third coordinate
  INTF='len-rec.dat' INTFOUT='len-rec.tmp' KSRFC=2 KOLUMN=3
  intf: 'lenb-crt.h' /

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

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

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

# List of output files with travel times
# 'lenb-rec.out' File with computed travel times in the model
#                with a LENticular inclusion