Changeset 1236


Ignore:
Timestamp:
04/07/20 16:56:53 (3 months ago)
Author:
klaus
Message:

bugfixes

Location:
SHX/trunk/SeismicHandler
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/cmdscripts/READ_GERMAN_STATIONS.SHC

    r1213 r1236  
    77sdef lon 11.6156 
    88 
     9echo perf__ reading starts at $systime 
    910fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 
    1011fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=eth /excludenets=z3,s 
    11 fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=lmu 
     12fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=lmu /excludenets=z3 
    1213fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=resif /excludenets=z3 
    13 fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=orfeus 
     14fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=orfeus /excludenets=z3,nl 
    1415fdsnws |"lat|,|"lon|,|6.| * * ?h? #1 #2 /addr=gfz /excludenets=z3 
     16echo perf__ reading ended at $systime 
     17echo perf__ $tottrcs traces read 
    1518chanselect hh,eh,bh ,00,10,11 
     19echo perf__ meta started at $systime 
    1620meta all complete /addr=gfz 
    1721meta all complete /addr=orfeus 
     
    2024meta all complete /addr=eth 
    2125meta all complete /del_incomplete 
     26echo perf__ meta ended at $systime 
     27echo perf__ $tottrcs traces remained 
    2228 
    2329 
  • SHX/trunk/SeismicHandler/commands/magnitude.py

    r1234 r1236  
    428428     
    429429    def determineMagnitudeMl( self, trc, stime, etime, degdist ): 
     430        if stime >= etime: 
     431            return None 
    430432        # find horizontals 
    431433        horizcomp = "NE12XY" 
     
    450452            sl = trc.slice_relative( stime, etime ) 
    451453            #ampl.append( max( abs(max(sl.data)), abs(min(sl.data)) ) ) 
     454            if len(sl.data) < 10: 
     455                print "dbg: short data on", trc.stats.station, trc.stats.channel 
     456                return None 
    452457            mean = sl.data.mean() 
    453458            ampl.append( max(sl.data.max()-mean, mean-sl.data.min()) ) 
  • SHX/trunk/SeismicHandler/patches/obspy_.py

    r1191 r1236  
    482482    if end > hi: 
    483483        end = hi 
    484  
     484     
    485485    start = self.stats.starttime + start - torigin 
    486486    end = self.stats.starttime + end - torigin 
     487     
     488    if start > end: 
     489        end = start 
    487490 
    488491    return self.slice(start, end) 
  • SHX/trunk/SeismicHandler/tools/netstodb.s

    r1233 r1236  
     1# update station metadata for SHX doing magnitude estimations: 
     2# run this script on smpc52: 'source netstodb.s' 
     3# copy DB file to shx1f: 'scp ~/.shx/inventory.db detector@shx1f:.shx/' 
    14~/svn/SHX/SeismicHandler/tools/FdsnMeta2Db.py SX,TH,GR BGR 
    25~/svn/SHX/SeismicHandler/tools/FdsnMeta2Db.py CH ETH 
Note: See TracChangeset for help on using the changeset viewer.