# Controlled initial-value ray tracing in model SALT
# followed by plotting of the computed rays
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file calculates rays shot in the x1-x3 plane
# and plots the computed rays.
# This history file is designed to be run from another history file
# (see <A HREF="sal-test.h">'sal-test.h'</A>), which runs 'srp.for'
# to generate source coordinates 'sal-src.tmp' and then executes this
# history file for the source.

# Input files required
 #chk.pl: "data/sal/" "sal-dcr.dat"
 #chk.pl: "data/sal/" "sal-rp2.dat"
 #chk.pl: "crt/"      "codep.dat"
 #chk.pl: "crt/"      "writall.dat"
 #chk.pl: ""          "sal-m.out"
 #chk.pl: ""          "sal-src.tmp"

# Model
  MODEL='sal-m.out'

# Source
  SRC='sal-src.tmp'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='sal-dcr.dat'
  RPAR='sal-rp2.dat'
  CODE='codep.dat'       # just refracted P wave
  WRIT='writall.dat'     # all rays stored for interpolation
  CRTLOG='s-crtlog.out'

# Kind of ray parameters (for details see 'crtin.for')
  INIPAR=3

# Ray tracing of refracted wave
  crt:

# Plotting the rays
  RAYS='sal-ray.out' KALL=1
  crtray:
  LIN='s-bnd.out'
  bndlin:
  PICDAT='picdat.tmp' PICTURE='sal-ray.ps'
  echo.pl: "'PROJECTION'                            " "> picdat.tmp"
  echo.pl: "1.5  6.5  0.001 0.  0. 0.  0. -0.001  / " ">>picdat.tmp"
  echo.pl: "                                        " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                            " ">>picdat.tmp"
  echo.pl: "'LWIDTH' 0.35  / Linewidth scale factor." ">>picdat.tmp"
  echo.pl: "/                                       " ">>picdat.tmp"
  echo.pl: "                                        " ">>picdat.tmp"
  echo.pl: "'LINES'                                 " ">>picdat.tmp"
  echo.pl: "'s-bnd.out'                             " ">>picdat.tmp"
  echo.pl: "                                        " ">>picdat.tmp"
  echo.pl: "'LINES'                                 " ">>picdat.tmp"
  echo.pl: "'sal-ray.out'                           " ">>picdat.tmp"
  echo.pl: "'PROJECTION'                            " ">>picdat.tmp"
  echo.pl: "16.  6.5  0. 0.  0.0002 0. 0. -0.001  / " ">>picdat.tmp"
  echo.pl: "                                        " ">>picdat.tmp"
  echo.pl: "'LINES'                                 " ">>picdat.tmp"
  echo.pl: "'s-bnd.out'                             " ">>picdat.tmp"
  echo.pl: "                                        " ">>picdat.tmp"
  echo.pl: "'LINES'                                 " ">>picdat.tmp"
  echo.pl: "'sal-ray.out'                           " ">>picdat.tmp"
  pictures:

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