def test_diffracted_phases(self): """ Test of exotic diffracted phases. """ model = TauPyModel('iasp91') phs = [ "SedPdiffKP", "PdiffPdiff", "PedPdiffKKP", "PdiffKKPdiff", "PPdiff", "SKdiffP" ] dists = [155.0, 210.0, 310.0, 300.0, 220.0, 200.0] times = [1464.97, 1697.88, 2052.42, 2008.03, 1742.27, 1431.53] for ph, dist, time in zip(phs, dists, times): phase = SeismicPhase(ph, model.model) arrivals = phase.calc_time(dist) arrival = arrivals[0] phase.calc_pierce_from_arrival(arrival) phase.calc_path_from_arrival(arrival) tol = 1e-2 assert abs(arrival.time - time) < tol assert abs(arrival.pierce["time"][-1] - time) < tol assert abs(arrival.path["time"][-1] - time) < tol