From 85e1e45552762d1c3692aa0e2754c59c71bf36b2 Mon Sep 17 00:00:00 2001
From: Darius Arnold <Darius.Arnold@ruhr-uni-bochum.de>
Date: Sat, 26 Aug 2017 19:41:26 +0200
Subject: [PATCH] [change] wadaticheck and jackknife 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 e1ae92db..d4d64ce6 100644
--- a/pylot/core/pick/utils.py
+++ b/pylot/core/pick/utils.py
@@ -617,6 +617,7 @@ def wadaticheck(pickdic, dttolerance, iplot=0, fig_dict=None):
checkedPpicks = []
checkedSpicks = []
checkedSPtimes = []
+ badstations = []
# calculate deviations from Wadati regression
ii = 0
ibad = 0
@@ -630,6 +631,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'
@@ -641,6 +643,8 @@ def wadaticheck(pickdic, dttolerance, iplot=0, fig_dict=None):
checkedSPtimes.append(checkedSPtime)
pickdic[key]['S']['marked'] = marker
+ print("wadaticheck: the following stations failed the check:")
+ print(badstations)
if len(checkedPpicks) >= 3:
# calculate new slope
@@ -862,6 +866,7 @@ def checkPonsets(pickdic, dttolerance, jackfactor, 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'