1from obspy.taup import TauPyModel
2
3model = TauPyModel(model='iasp91')
4arrivals = model.get_ray_paths(150, 40, phase_list=['P', 'S'])
5arrivals.plot_rays(plot_type='cartesian', phase_list=['P', 'S'],
6                   plot_all=False, legend=True)
7