# History file 'prem-crt.h' to trace pPKP+pPKiKP, pPKIKP,
# PKP+PKiKP and  PKIKP rays in PREM
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
  chk.pl: "data/prem/" "prem-mod.dat"
  chk.pl: "data/prem/" "prem-crt.dat"
  chk.pl: "data/prem/" "prem-src.dat"
  chk.pl: "data/prem/" "prem-cod.dat"
  chk.pl: "data/prem/" "prem-rpa.dat"
  chk.pl: "data/prem/" "prem-rec.dat"
  chk.pl: "data/prem/" "prem-wri.dat"
  chk.pl: "forms/"     "echo.pl"

# Model
  MODEL='prem-mod.dat'

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

# Data filenames for CRT (for details see 'crt.for')
  DCRT='prem-crt.dat'
  CODE='prem-cod.dat'
  RPAR='prem-rpa.dat'
  WRIT='prem-wri.dat'
  CRTLOG='prem-log.out'

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

# Ray tracing of 4 waves
  crt: 'prem-crt.h' /

# Rewriting the travel times at the receivers into formatted files
  SRC=' ' NQ=4 KALL=0 KREC=1 KTT=0 CRTOUT='crtout.tmp'
  echo.pl:"' '    'ppkp.tmp'   'ppkp1.tmp'  ' '" ">crtout.tmp"
    PTS='ppkp.out' crtpts: 'prem-crt.h' /
  echo.pl:"' '  'ppkikp.tmp' 'ppkikp1.tmp'  ' '" ">crtout.tmp"
  PTS='ppkikp.out' crtpts: 'prem-crt.h' /
  echo.pl:"' '     'pkp.tmp'    'pkp1.tmp'  ' '" ">crtout.tmp"
     PTS='pkp.out' crtpts: 'prem-crt.h' /
  echo.pl:"' '   'pkikp.tmp'  'pkikp1.tmp'  ' '" ">crtout.tmp"
   PTS='pkikp.out' crtpts: 'prem-crt.h' /

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

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

# List of output files with travel times
# 'ppkp.out'   \
# 'ppkikp.out'  > Files with computed travel times in the model PREM
# 'pkp.out'    /
# 'pkikp.out' /