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

Claculation of Wood-Anderson amplitude only, if S-weight < 4.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-08-28 11:39:39 +02:00
parent 827a07a7ef
commit eb592a3426
1 changed files with 14 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -507,19 +507,20 @@ def autopickstation(wfstream, pickparam):
print 'autopickstation: S-weight: %d, SNR: %f, SNR[dB]: %f' % (
Sweight, SNRS, SNRSdB)
##################################################################
# get Wood-Anderson peak-to-peak amplitude
print "################################################"
# initialize Data object
data = Data()
# re-create stream object including both horizontal components
hdat = edat.copy()
hdat += ndat
h_copy = hdat.copy()
cordat = data.restituteWFData(invdir, h_copy)
# calculate WA-peak-to-peak amplitude
# using subclass WApp of superclass Magnitude
wapp = WApp(cordat, mpickS, 10, iplot)
Ao = wapp.getwapp()
if Sweight < 4:
# get Wood-Anderson peak-to-peak amplitude
print "################################################"
# initialize Data object
data = Data()
# re-create stream object including both horizontal components
hdat = edat.copy()
hdat += ndat
h_copy = hdat.copy()
cordat = data.restituteWFData(invdir, h_copy)
# calculate WA-peak-to-peak amplitude
# using subclass WApp of superclass Magnitude
wapp = WApp(cordat, mpickS, 10, iplot)
Ao = wapp.getwapp()
else:
print 'Bad initial (AIC) S-pick, skipping this onset!'