# Zero-offset rays in model Pluto 1.5
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file calculates zero-offset rays
# and plots the calculated rays to the output file 'plu-ray.ps'.

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: ""          "plu-mod.out" # output of 'plu-inv.h'
  chk.pl: "data/plu/" "plu-crt.dat"
  chk.pl: "data/plu/" "plu-rpaz.dat"
  chk.pl: "crt/"      "codep.dat"
  chk.pl: "crt/"      "writall.dat"
  chk.pl: "forms/"    "echo.pl"

# Calculating rays
# ~~~~~~~~~~~~~~~~
# Model
  MODEL='plu-mod.out'

# Initial conditions for rays
  INIDIM=1             # initial line
  INIPAR=1             # ray parameters = X1
  INIT='plu-rpaz.dat'  # zero-offset rays (exploding reflector)
  RPAR='plu-rpaz.dat'  # data for initial-value ray tracing

# Data filenames for CRT
  DCRT='plu-crt.dat'
  CODE='codep.dat'       # just refracted P wave
  WRIT='writall.dat'     # all rays stored
  CRTLOG='plu-logz.out'

# Ray tracing of refracted wave
  crt:

# Plotting rays
# ~~~~~~~~~~~~~
  RAYS='plu-rays.out' KALL=1
  crtray:
  LIN='plu-bnd.out'
  bndlin:
  PICDAT='picdat.tmp' PICTURE='plu-ray.ps'
  echo.pl: "'PROJECTION'                              " "> picdat.tmp"
  echo.pl: "6. 2.35 0.0001438 0. 0. 0.0001438 0. 0. / " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                              " ">>picdat.tmp"
  echo.pl: "'LWIDTH' 0.1   / Linewidth scale factor.  " ">>picdat.tmp"
  echo.pl: "'ITC'    0     / No writing of texts.     " ">>picdat.tmp"
  echo.pl: "/                                         " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'plu-bnd.out'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'plu-rays.out'                            " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  pictures:

# Important output files
# ~~~~~~~~~~~~~~~~~~~~~~
# 'plu-rays.out'... Computed rays.
# 'plu-ray.ps'  ... Picture of the rays.