ERZSOL3 (Reflectivity Method)

Code for calculating the response of a stratified set of uniform solid layers to excitation by a point moment tensor source using the Reflectivity Method (RM).

Project Description

Author: Brian Kennett

The routine ERZSOL3 is designed for calculating the response of a stratified set of uniform solid layers to excitation by a point moment tensor source. The calculations using slowness-frequency domain synthesis include all near field terms.

The model includes attenuation - constant Q (without causal corrections) with up to 500 layers in the standard settings. Integrations are carried out over slowness (2500 slowness points currently allowed) for each distance and then converted to time using an FFT.

The seismograms are written in a very simple format (ZST) that can be plotted using the accompanying zst routine or readily converted to other formats.

The algorithm is described in

• "Seismic Waves Propagation in Stratified Media" by B.L.N. Kennett first published by Cambridge University press, 1993.