#!/bin/sh
### ====================================================================
### Two-point ray tracing in the model with LENticular inclusion (Unix)
### ====================================================================
### Ray tracing (refracted wave):
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "'len-crt.dat'  /" >crt.tmp
echo                   >>crt.tmp
    crt                 <crt.tmp
rm                       crt.tmp
### Note: The blank line in crt.tmp is required only if 'scronum.for' is
###       replaced with 'scropc.for'.  Otherwise, the above lines can be
###       simplified to a single line:
###       echo "'len-crt.dat'  /"|crt
### ....................................................................
### Supplementing receivers with the third coordinate:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "'model.dat' 'len-rec.dat' 'len-rec.tmp' 2 3 /"|intf
### Note: File 'len-rec.dat' contains only X1,X2 receiver coordinates.
###       If it were submitted to program 'crtpts.for' as the receiver
###       file, the travel time could be extrapolated in the horizontal
###       plane to the receiver position, but the vertical X3 coordinate
###       would correspond to the numerical approximation of the point
###       of intersection with the reference surface (Earth surface)
###       rather than to the exact X3 coordinate of the receiver.
### ....................................................................
### Rewriting the travel times at the receivers into a formatted file:
### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "'len-rec.tmp' 3* 'len-rec.out'  4 0 1 0 /"|crtpts
### ====================================================================