def pspier(self, depth, slat, slon, slowness, azi, phase='P'):
""""
Calculate piercing points of events.
:param depth: depth of pp
:param slat: station latitude
:param slat: station longitude
:param slowness: slowness of wave
:param azi: azimuth of wave
:param model: use this model (instance of seis.model.EarthModel1D)
:return: (distance, latitude of pp, longitude of pp"""
phi = self.trace3(slowness, phase=phase, till_turn=True)[2] #@UnusedVariable
x = phi[self.layer(depth)] * r_earth
return (x,) + dist2gps(x, azi, slat, slon)
def pspier(depth, slat, slon, slowness, azi, model=iasp91):
""""
Calculate piercing points of events.
:param depth: depth of pp
:param slat: station latitude
:param slat: station longitude
:param slowness: slowness of wave
:param azi: azimuth of wave
:param model: use this model (instance of seis.model.EarthModel1D)
:return: (distance, latitude of pp, longitude of pp"""
tp, ts, xp, xs = model.trace(slowness) #@UnusedVariable
xs = xs[model.layer(depth)]
return (xs,) + dist2gps(xs, azi, slat, slon)
def pspier(depth, slat, slon, slowness, azi, model=iasp91):
""""
Calculate piercing points of events.
:param depth: depth of pp
:param slat: station latitude
:param slat: station longitude
:param slowness: slowness of wave
:param azi: azimuth of wave
:param model: use this model (instance of seis.model.EarthModel1D)
:return: (distance, latitude of pp, longitude of pp"""
tp, ts, xp, xs = model.trace(slowness) #@UnusedVariable
xs = xs[model.layer(depth)]
return (xs, ) + dist2gps(xs, azi, slat, slon)
def pspier(self, depth, slat, slon, slowness, azi, phase='P'):
""""
Calculate piercing points of events.
:param depth: depth of pp
:param slat: station latitude
:param slat: station longitude
:param slowness: slowness of wave
:param azi: azimuth of wave
:param model: use this model (instance of seis.model.EarthModel1D)
:return: (distance, latitude of pp, longitude of pp"""
phi = self.trace3(slowness, phase=phase,
till_turn=True)[2] #@UnusedVariable
x = phi[self.layer(depth)] * r_earth
return (x, ) + dist2gps(x, azi, slat, slon)
