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

Marginal changes.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-06-26 16:00:20 +02:00
parent 99adb5ce9c
commit 27155ad816
1 changed files with 14 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
autoPyLoT.py
View File

@ -12,7 +12,7 @@ from pylot.core.util import _getVersionString
from pylot.core.read import Data, AutoPickParameter
from pylot.core.pick.run_autopicking import run_autopicking
from pylot.core.util.structure import DATASTRUCTURE
from pylot.core.pick.utils import wadaticheck
from pylot.core.pick.utils import wadaticheck, checkPonsets
import pdb
__version__ = _getVersionString()
@ -51,6 +51,7 @@ def autoPyLoT(inputfile):
# get some parameters for quality control from
# parameter input file (usually autoPyLoT.in).
wdttolerance = parameter.getParam('wdttolerance')
mdttolerance = parameter.getParam('mdttolerance')
iplot = parameter.getParam('iplot')
data = Data()
@ -105,10 +106,11 @@ def autoPyLoT(inputfile):
allonsets[station] = picks
# quality control
# jackknife on P onset times
# median check and jackknife on P onset times
checkedonsetsjk = checkPonsets(allonsets, mdttolerance, iplot)
# check S-P times (Wadati)
checkedonsets = wadaticheck(allonsets, wdttolerance, iplot)
# jackknife on S onset times
checkedonsetwd = wadaticheck(checkedonsetsjk, wdttolerance, iplot)
print '------------------------------------------'
print '-----Finished event %s!-----' % event
print '------------------------------------------'
@ -128,6 +130,7 @@ def autoPyLoT(inputfile):
station = wfdat[0].stats.station
allonsets = {station: picks}
for i in range(len(wfdat)):
#for i in range(0,10):
stationID = wfdat[i].stats.station
#check if station has already been processed
if stationID not in procstats:
@ -143,10 +146,11 @@ def autoPyLoT(inputfile):
allonsets[station] = picks
# quality control
#jackknife on P onset times
# median check and jackknife on P onset times
checkedonsetsjk = checkPonsets(allonsets, mdttolerance, iplot)
# check S-P times (Wadati)
checkedonsets = wadaticheck(allonsets, wdttolerance, iplot)
#jackknife on S onset times
checkedonsetswd = wadaticheck(checkedonsetsjk, wdttolerance, iplot)
print '------------------------------------------'
print '-------Finished event %s!-------' % parameter.getParam('eventID')
print '------------------------------------------'