Plotting Raspberry Shake instrument response in ObsPy

This has probably been asked and answered before, but:

What is the ObsPy code for plotting the Raspberry Shake instrument response?

Hi Alan, see code below:

from obspy import read_inventory, UTCDateTime
from datetime import timedelta
sta = 'R24FA'
now = UTCDateTime.now() - timedelta(days=1)
inv = read_inventory('https://fdsnws.raspberryshakedata.com/fdsnws/station/1/query?network=AM&station=%s&level=resp&format=xml&nodata=404&starttime=%s' % (sta, str(now)))
inv.plot_response(min_freq=0.01)

Thanks, Ian!

Just what I needed :slight_smile:

2 Likes