From 7281220b5571df1edae5c006912463897b9b3675 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Wed, 24 Jun 2015 16:42:50 +0200 Subject: [PATCH] wadaticheck: bug fixed, initialization of running index at wrong place. --- pylot/core/pick/utils.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/pylot/core/pick/utils.py b/pylot/core/pick/utils.py index c869ef69..321348b9 100644 --- a/pylot/core/pick/utils.py +++ b/pylot/core/pick/utils.py @@ -449,11 +449,11 @@ def wadaticheck(pickdic, dttolerance, iplot): checkedSpicks = [] checkedSPtimes = [] # calculate deviations from Wadati regression + ii = 0 for key in pickdic: if pickdic[key].has_key('SPt'): - ii = 0 wddiff = abs(pickdic[key]['SPt'] - wdfit[ii]) - ii += 1 + ii += 1 # check, if deviation is larger than adjusted if wddiff >= dttolerance: # mark onset and downgrade S-weight to 9 @@ -489,7 +489,7 @@ def wadaticheck(pickdic, dttolerance, iplot): print 'wadaticheck: Not enough S-P times available for reliable regression!' print 'Skip wadati check!' wfitflag = 1 - + iplot=2 # plot results if iplot > 1: plt.figure(iplot) @@ -498,8 +498,8 @@ def wadaticheck(pickdic, dttolerance, iplot): f2, = plt.plot(Ppicks, wdfit, 'k') f3, = plt.plot(checkedPpicks, checkedSPtimes, 'ko') f4, = plt.plot(checkedPpicks, wdfit2, 'g') - plt.title('Wadati-Diagram, %d S-P Times, Vp/Vs(raw)=%5.2f, \ - Vp/Vs(checked)=%5.2f' % (len(SPtimes), vpvsr, cvpvsr)) + plt.title('Wadati-Diagram, %d S-P Times, Vp/Vs(raw)=%5.2f,' \ + 'Vp/Vs(checked)=%5.2f' % (len(SPtimes), vpvsr, cvpvsr)) plt.legend([f1, f2, f3, f4], ['Skipped S-Picks', 'Wadati 1', \ 'Reliable S-Picks', 'Wadati 2'], loc='best') else: