def test_degrees2kilometers(self):
"""
"""
# Test if it works.
self.assertEqual(degrees2kilometers(1.0, radius=6371),
111.19492664455873)
# Test if setting the radius actually does something. Round to avoid
# some precision problems on different machines.
self.assertEqual(round(degrees2kilometers(1.0, radius=6381), 5),
round(111.36945956975816, 5))
def test_degrees2kilometers(self):
"""
"""
# Test if it works.
self.assertEqual(degrees2kilometers(1.0, radius=6371),
111.19492664455873)
# Test if setting the radius actually does something. Round to avoid
# some precision problems on different machines.
self.assertEqual(round(degrees2kilometers(1.0, radius=6381), 5),
round(111.36945956975816, 5))
def test_ppointvsobspytaup_S2P(self):
slowness = 12.33
evdep = 12.4
evdist = 67.7
pp1 = self.model.ppoint_distance(200, slowness, phase="P")
model = TauPyModel(model="iasp91")
arrivals = model.get_ray_paths(evdep, evdist, ("S250p",))
arrival = arrivals[0]
index = np.searchsorted(arrival.path["depth"][::-1], 200)
pdist = arrival.path["dist"]
pp2 = degrees2kilometers((pdist[-1] - pdist[-index - 1]) * 180 / np.pi)
self.assertLess(abs(pp1 - pp2) / pp2, 0.2)
def test_ppointvsobspytaup_S2P(self):
slowness = 12.33
evdep = 12.4
evdist = 67.7
pp1 = self.model.ppoint_distance(200, slowness, phase='P')
model = TauPyModel(model='iasp91')
arrivals = model.get_ray_paths(evdep, evdist, ('S250p', ))
arrival = arrivals[0]
index = np.searchsorted(arrival.path['depth'][::-1], 200)
pdist = arrival.path['dist']
pp2 = degrees2kilometers((pdist[-1] - pdist[-index - 1]) * 180 / np.pi)
self.assertLess(abs(pp1 - pp2) / pp2, 0.2)
def test_ppointvsobspytaup_P2S(self):
slowness = 6.28
evdep = 12.4
evdist = 67.7
depth = 200
pp1 = self.model.ppoint_distance(depth, slowness)
model = TauPyModel(model="iasp91")
arrivals = model.get_ray_paths(evdep, evdist, ("P250s",))
arrival = arrivals[0]
index = np.searchsorted(arrival.path["depth"][::-1], depth)
pdist = arrival.path["dist"]
pp2 = degrees2kilometers((pdist[-1] - pdist[-index - 1]) * 180 / np.pi)
self.assertLess(abs(pp1 - pp2) / pp2, 0.1)
def test_ppointvsobspytaup_P2S(self):
slowness = 6.28
evdep = 12.4
evdist = 67.7
depth = 200
pp1 = self.model.ppoint_distance(depth, slowness)
model = TauPyModel(model='iasp91')
arrivals = model.get_ray_paths(evdep, evdist, ('P250s', ))
arrival = arrivals[0]
index = np.searchsorted(arrival.path['depth'][::-1], depth)
pdist = arrival.path['dist']
pp2 = degrees2kilometers((pdist[-1] - pdist[-index - 1]) * 180 / np.pi)
self.assertLess(abs(pp1 - pp2) / pp2, 0.1)
