Documentation

PyRaysum is a Python wrapper around the Fortran software Raysum, originally developed by Andrew Frederiksen in collaboration with Michael Bostock. This program generates sets of ray-theoretical seismograms for an incident plane wave (teleseismic approximation) for models consisting of a stack of layers with planar (but possibly non-parallel and dipping) interfaces, allowing the flexibility of adding seismic anisotropy in the layers. Incident P and S waves are supported. The code includes methods to process the synthetic seismograms for receiver function calculation and for plotting.

PyRaysum is bundled with a trimmed and streamlined version of the original Fortran software, and provides a module to interact with Raysum, as well as functions to post-process the data and generate plots of models and seismograms. Common computational workflows are covered in the Jupyter notebooks bundled with this package.

Quick Links

• GitHub Repository
• References

Getting Started

• Licence
• Installation
• Usage

API

• Classes
• Modules

Tutorials

• Example 1: Compare PyRaysum with other synthetic and real data
• Example 2: Reproducing Receiver Functions for Dipping and Anisotropic Subsurface Structure
• Example 3: Invert a Receiver Function

Jupyter Notebooks

• Example 1: Compare PyRaysum with other synthetic and real data
• Example 2: Reproducing Figure 2 in Porter et al. (2011)
• Example 3: Invert a receiver function