marcel/pylot
3
0
Fork 0
Code Issues 7 Pull Requests 1 Releases Wiki Activity

Debuged: Checks additionally for component 3 if component Z not available.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-11-11 13:58:25 +01:00
parent 02dbca06f2
commit 4f67e3e31a
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
pylot/core/pick/autopick.py
View File

@ -11,6 +11,7 @@ function conglomerate utils.
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import numpy as np import numpy as np
import pdb
from scipy import integrate from scipy import integrate
from pylot.core.pick.Picker import AICPicker, PragPicker from pylot.core.pick.Picker import AICPicker, PragPicker
from pylot.core.pick.CharFuns import HOScf, AICcf, ARZcf, ARHcf, AR3Ccf from pylot.core.pick.CharFuns import HOScf, AICcf, ARZcf, ARHcf, AR3Ccf
@ -28,7 +29,6 @@ def autopickevent(data, param):
wdttolerance = param.getParam('wdttolerance') wdttolerance = param.getParam('wdttolerance')
mdttolerance = param.getParam('mdttolerance') mdttolerance = param.getParam('mdttolerance')
iplot = param.getParam('iplot') iplot = param.getParam('iplot')
for n in range(len(data)): for n in range(len(data)):
station = data[n].stats.station station = data[n].stats.station
if station not in stations: if station not in stations:
@ -141,6 +141,8 @@ def autopickstation(wfstream, pickparam):
# split components # split components
zdat = wfstream.select(component="Z") zdat = wfstream.select(component="Z")
if len(zdat) == 0: # check for other components
zdat = wfstream.select(component="3")
edat = wfstream.select(component="E") edat = wfstream.select(component="E")
if len(edat) == 0: # check for other components if len(edat) == 0: # check for other components
edat = wfstream.select(component="2") edat = wfstream.select(component="2")