Instaseis Tutorial
Part 3: Plot a Record Section

Seismo-Live: http://seismo-live.org

Authors:

Advanced Exercise 1: Plot Record Section

Use Instaseis to calculate a record section of your choice.

image

Notes

  1. ObsPy Trace objects provide a function times() to conveniently get the time axis for plotting and the normalize() function to normalize seismograms to a common amplitude.
  2. Use high pass filtering and deep sources if you want to enhance body waves
  3. np.linspace() can help to generate equidistant stations
  4. obspy.taup.TauPyModel can be used to compare with ray theoretical arrivals

Basic lines to set up the notebook and some paths.

In [ ]:
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
import os
import obspy
plt.style.use('ggplot')
plt.rcParams['figure.figsize'] = (10, 8)

Import Instaseis and open the database:

In [ ]:
import instaseis
db = instaseis.open_db("data/database")

Solution

In [ ]: