1from obspy.taup import TauPyModel
2
3model = TauPyModel(model='iasp91')
4# 500 is from 0 to 700 km
5arrivals = model.get_ray_paths(500, 45, phase_list=['p','s','P','S'])
6arrivals.plot_rays(plot_type='spherical', phase_list=['p','s','P', 'S'],
7                   legend=True)