Click on each page link above to see each scanned page of code.
I have included page 1 because it may relate to coding a Klauder wavelet.
Here are some notes included on the wavelet design:
Both sides of a zero phase wavelet are stored.
The GSI type filters are the same design as in GSI's regular seismic processing
routines.
The Ricker wavelets are one sided, but not minimum phase in digital form,
although fairly close to it. (Their design equations are from an article by SSL.)
The number of points used is based on the wavelet type and frequency, and is
sufficient to include all significant amplitudes.
The Butterworth filters have cut-off slopes that are constant in db/octave,
provided the wavelet is long enough. The cut-off frequencies are the -6 db
points on the response.
This code was written in Basic by John Wardell
I'll leave them online for a while. Contact me if you need more information.
Cheers.,
Robert.