Summary: Computes a receiver-function for a stack of anisotropic layers over an isotropic halfspace, via a reflectivity algorithm, assuming a plane wave incident from below

Author(s): Jeffrey Park, Vadim Levin

Description: The code reads a 1D model of velocity, density and anisotropy layers with depth, computes 3 components of motion for the response of this
model to an incident compressional pulse. On output, the code computes receiver functions from the division of smoothed spectra of respective
components. Adjusts for the phase velocity and BAZ of the incident compressional phase. Computes waveforms up to a user-specified frequency. README