#!/bin/sh
### ======================================================================
### Unix script to trace pPKP+pPKiKP,pPKIKP,PKP+PKiKP,PKIKP rays in PREM
### ======================================================================
### Ray tracing in PREM (4 elementary waves):
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "'prem-crt.dat' /" >crt.tmp
echo                   >>crt.tmp
echo                   >>crt.tmp
echo                   >>crt.tmp
echo                   >>crt.tmp
    crt                 <crt.tmp
rm                       crt.tmp
### Note: 4 blank lines in crt.tmp are required only if 'scronum.for' is
###       replaced with 'scropc.for'.  Otherwise, the above lines can be
###       simplified to a single line:
###       echo "'prem-crt.dat' /"|crt
### ......................................................................
### Rewriting the travel times at the receivers into formatted files:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "'prem-rec.dat' ' '   'ppkp.tmp'   'ppkp1.tmp'   'ppkp.out' 4 0 1 0 /"|crtpts
echo "'prem-rec.dat' ' ' 'ppkikp.tmp' 'ppkikp1.tmp' 'ppkikp.out' 4 0 1 0 /"|crtpts
echo "'prem-rec.dat' ' '    'pkp.tmp'    'pkp1.tmp'    'pkp.out' 4 0 1 0 /"|crtpts
echo "'prem-rec.dat' ' '  'pkikp.tmp'  'pkikp1.tmp'  'pkikp.out' 4 0 1 0 /"|crtpts
### ======================================================================