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)