sfricker1 implements 1-D convolution with the Ricker wavelet.

The following example from rsf/rsf/wedge shows convolution modeling with a wedge model using sfricker1.

The convolution is implemented in the frequency domain, where the Ricker wavelet takes the form
$$F(\omega) = \frac{\omega^2}{\omega_0^2}\,\exp\left(-\frac{\omega^2}{\omega_0^2}\right)$$
The peak frequency $\omega_0$ can be specified either in hertz (with frequency= parameter) as a fraction of the Nyquist frequency (with freq= parameter). Optionally, it is possible to combine the Ricker wavelet with a half-order derivative filter using deriv=y.

10 previous programs of the month