Receiver function inversion using wavefield decomposition and simultaneous frequency-domain deconvolution of multichannel measurements, based on the method of Bostock, 1998. These codes were for calculating Ps receiver function for stations in Ontario and Quebec part of the Superior geological Province as part of a research project with Metal Earth in 2018.
This code is not currently properly packaged to be distributed.