c anirec - program to calculate receiver-function response of <br> 
c a stack of anisotropic layers to a plane wave incident from below <br> 
c CAN BE USED IN GRIDSEARCH OVER CC AND BAZ - SEE COMMENTED LINES FOR LOOP <br> 
c cannibalized from aniprop.f 11/18/95 <br> 
c <br> 
c xf77 -o anirec -fast -native -O5 anirec.f plotlib refft.o <br> 
c <br> 
c refft.o is an object file of the C-routine refft.c <br> 
c compile this code with <br> 
c cc -c -fast -native -O5 refft.c <br> 
c <br> 
c The "-fast -native" option will generate a faster executable. Note that <br> 
c <br> 
c alias xf77 'f77 -I${OPENWINHOME}/include -L${OPENWINHOME}/lib -lxview -lolgx -lX11' <br> 
c <br> 
c plotlib is PLOTIT plotting package, used to plot results thruout the code <br> 
c PLOTIT resides at love.geology.yale.edu FTP site at pub/park/Plotxy <br> 
c <br> 
c for hexagonally symmetric media <br> 
c reads fast axis orientation, constants A,B,C,D,E from file animodel <br> 
c c reads a layered model like the following file (ignore the "c "s) <br> 
c <br> 
c K&H SoCal model, deep crust horizontal anisotropy <-- TITLE <br> 
c 3 <--- # OF LAYERS OVER HSPACE <br> 
c 45 45 <--- THETA,PHI ORIENTATION ANGLES <br> 
c 4000 5500 0.06 0.00 3175 0.03 2600 FOR SYMMETRY AXIS <br> 
c 0 0 <br> 
c 27400 6300 0.00 0.00 3637 0.00 2800 <-- DEPTH(M),VP(M/S),"B","C",VS(M/S),"E" <br> 
c 90 45 <br> 
c 32400 6800 0.04 0.00 3925 0.02 2900 <-- B,C,E ARE ANISOTROPIC PARAMETERS <br> 
c 0 0 <br> 
c 60000 7800 0.00 0.00 4500 0.00 3200 <--- HALFSPACE MUST BE ISOTROPIC <br> 
c <br> 
c calculate quadratic eigenvalue problem based on the Christoffel matrix <br> 
c see appendix of P. Shearer's thesis <br> 
c <br> 
c read model, phase velocity of incident wave, loop thru incident P, SV, SH <br> 
c <br> 
c calc the eigenvector decomps for the layers <br> 
c loop over frequency, calc reflection/transmission matrices <br> 
c calc 3-comp transfer fct response at surface <br> 
c find distortion of reference wavelet <br> 
c <br>