# Initial-value ray tracing in model Pluto 1.5
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file calculates the rays corresponding to the point
# source specified in file 'plu-src.tmp' and plots the calculated
# rays to the output file 'plu-r.ps'.
# This history file is designed to be run from another history file
# (see <A HREF="plu-test.h">'plu-test.h'</A>), which runs 'srp.for'
# to generate source coordinates 'plu-src.tmp' and then executes this
# history file for the source.

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

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

# Source
  SRC='plu-src.tmp'

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

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

# 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: "/                                         " ">>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"
  pictures:

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