From a96f8dcb7dce48d3b6ab4dd225dc6fb1e99aea96 Mon Sep 17 00:00:00 2001 From: Darius Arnold Date: Fri, 15 Sep 2017 13:24:13 +0200 Subject: [PATCH] [add] wadati, jackknife and median test print excluded stations --- pylot/core/pick/utils.py | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/pylot/core/pick/utils.py b/pylot/core/pick/utils.py index 92b4f23f..dd5a4567 100644 --- a/pylot/core/pick/utils.py +++ b/pylot/core/pick/utils.py @@ -616,6 +616,7 @@ def wadaticheck(pickdic, dttolerance, iplot=0, fig_dict=None): checkedPpicks = [] checkedSpicks = [] checkedSPtimes = [] + badstations = [] # calculate deviations from Wadati regression ii = 0 ibad = 0 @@ -631,6 +632,7 @@ def wadaticheck(pickdic, dttolerance, iplot=0, fig_dict=None): # # (not used anymore) # marker = 'badWadatiCheck' # pickdic[key]['S']['weight'] = 9 + badstations.append(key) ibad += 1 else: marker = 'goodWadatiCheck' @@ -643,6 +645,8 @@ def wadaticheck(pickdic, dttolerance, iplot=0, fig_dict=None): pickdic[key]['S']['marked'] = marker #pickdic[key]['S']['marked'] = marker + print("wadaticheck: the following stations failed the check:") + print(badstations) if len(checkedPpicks) >= 3: # calculate new slope @@ -868,6 +872,7 @@ def checkPonsets(pickdic, dttolerance, jackfactor=5, iplot=0, fig_dict=None): print("checkPonsets: %d pick(s) deviate too much from median!" % len(ibad)) print("checkPonsets: Skipped %d P pick(s) out of %d" % (len(badstations) \ + len(badjkstations), len(stations))) + print(badstations) goodmarker = 'goodPonsetcheck' badmarker = 'badPonsetcheck'