APPENDIX D 



Computer Program for Wave Refraction (Linear-Interpolation) 



Using the IBM I62O 



PROGRAM TITLE : MAIN I62O. 



Input VariaMes ; 



S,'M Matrices used in determining coefficients for the plane of 



best-fit to four data points. 



XLABL. Arbitrary title used to designate each batch of input. 



MM, M Maximum values of I and J, respectively, for the CMAT grid 



of velocity values; I=X+1, J=Y+1, with the grid origin 

 located at X = 0, Y = 0. 



CHECK Allows either a two-dimensional CMAT to be input, or a single- 

 rowed CMAT (extending from shore to deep water) to be input 

 and all other rows to be made identical to the first. (This 

 procedure creates a set of grid values which can be character- 

 ized by straight and parallel contours . ) 



TT Wave period (seconds) . 



NOJ Number of wave rays to be rim in each batch. 



D Distance incremented between successive points along a ray 



path. 



CMAT(I,J) Wave velocity (grid units/second) at grid position I, J. 



A Angle (degrees) measured from the direction of increasing X 



along the X-axis moving counter-clockwise to the direction of 

 travel of a wave ray. 



X, Y Grid origin position for a wave ray. 



40 



