Changeset 1130


Ignore:
Timestamp:
20.01.2016 15:09:34 (4 years ago)
Author:
klaus
Message:

bugfix in metadata update; bugfix in menue entries using fdsn reading

Location:
SHX/trunk/SeismicHandler
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/commands/meta.py

    r1127 r1130  
    120120        incomplete = [] 
    121121        startend = {} 
     122        oldgains = [] 
    122123        for trc in traces: 
    123             metastatus = get_meta_status( trc ) 
     124            metastatus, oldgain = get_meta_status( trc, addgain=True ) 
     125            oldgains.append( oldgain ) 
    124126            if metastatus == META_STATUS_COMPLETE: 
    125127                continue 
     
    137139            startend["%s.%s.%s.%s" % (trc.stats.network,trc.stats.station, 
    138140                trc.stats.location,trc.stats.channel)] = ( 
    139                 trc.stats.starttime,trc.stats.endtime) 
     141                trc.stats.starttime,trc.stats.endtime,oldgain) 
    140142        if incomplete == []: 
    141143            return 
     
    156158                print "got unrequested info for '%s', ignored" % chan 
    157159                continue 
    158             stime, etime = startend[chan] 
     160            stime, etime, oldgain = startend[chan] 
    159161            cm = inv.get_coordinates( chan ) 
    160162            try: 
     
    194196        triggerDatabaseReload() 
    195197        #self.stations.read( clear=True ) 
    196         for trc in traces: 
    197             trc.set_info( 'METASTATUS', get_meta_status( trc ) ) 
     198        for itrc, trc in enumerate(traces): 
     199            # check for correctly applied gain 
     200            oldgain = oldgains[itrc] 
     201            metstat, newgain = get_meta_status( trc, addgain=True ) 
     202            if (oldgain == None or oldgain == 1.) and newgain != 1.: 
     203                # multiply new gain onto trace data 
     204                trc.data *= newgain 
     205            trc.set_info( 'METASTATUS', metstat ) 
    198206     
    199207    def listIncompleteTraces( self, traces, subcmd ): 
  • SHX/trunk/SeismicHandler/modules/traces.py

    r1127 r1130  
    126126        return self.netdict.get( station.upper(), None ) 
    127127     
    128     def getMetaStatus( self, trace ): 
     128    def getMetaStatus( self, trace, addgain=False ): 
    129129        "Returns 0, 1 or 2 depending on metadata status." 
    130130        stations = Stations() 
     
    136136            r = stations[(sname, metatime)] 
    137137        except: 
    138             return META_STATUS_EMPTY 
     138            if addgain: 
     139                return (META_STATUS_EMPTY,None) 
     140            else: 
     141                return META_STATUS_EMPTY 
    139142        try: 
    140143            slat = float( r.latitude ) 
    141144            slon = float( r.longitude ) 
    142145        except: 
    143             return META_STATUS_EMPTY 
     146            if addgain: 
     147                return (META_STATUS_EMPTY,None) 
     148            else: 
     149                return META_STATUS_EMPTY 
    144150        try: 
    145151            gain = float( r.gain ) 
    146152        except: 
    147             return META_STATUS_LOCATION 
     153            if addgain: 
     154                return (META_STATUS_LOCATION,None) 
     155            else: 
     156                return META_STATUS_LOCATION 
    148157        if gain == 1.0: 
    149             return META_STATUS_LOCATION 
     158            if addgain: 
     159                return (META_STATUS_LOCATION,None) 
     160            else: 
     161                return META_STATUS_LOCATION 
    150162        if r.poles and r.poles != '[]': 
    151             return META_STATUS_COMPLETE 
     163            if addgain: 
     164                return (META_STATUS_COMPLETE,gain) 
     165            else: 
     166                return META_STATUS_COMPLETE 
    152167        else: 
    153             return META_STATUS_LOCATION 
     168            if addgain: 
     169                return (META_STATUS_LOCATION,gain) 
     170            else: 
     171                return META_STATUS_LOCATION 
    154172         
    155173 
     
    386404    Traces().addreplace(stream, applyGain=gain) 
    387405 
    388 def get_meta_status( trace ): 
    389     return Traces().getMetaStatus( trace ) 
     406def get_meta_status( trace, addgain=False ): 
     407    return Traces().getMetaStatus( trace, addgain=addgain ) 
    390408 
    391409 
  • SHX/trunk/SeismicHandler/modules/wx_.py

    r1129 r1130  
    15031503    def OnReadFdsnStations( self, e ): 
    15041504        ap = AnalysisPar() 
    1505         cmd = "fdsnws %s %s %s %s %s %s" % ( 
     1505        cmd = "fdsnws %s %s %s %s %s %s /addr=%s" % ( 
    15061506            ap.getValueAsString( 'readws_net' ), 
    15071507            ap.getValueAsString( 'readws_station' ), 
     
    15101510            ap.getValueAsString( 'read_time' ), 
    15111511            ap.getValueAsString( 'read_length' ), 
     1512            ap.getValueAsString( 'readws_server' ), 
    15121513        ) 
    15131514        print "dbg: cmd", cmd 
     
    15161517    def OnReadFdsnCircle( self, e ): 
    15171518        ap = AnalysisPar() 
    1518         cmd = "fdsnws %s,%s,%s * * %s %s %s" % ( 
     1519        cmd = "fdsnws %s,%s,%s * * %s %s %s /addr=%s" % ( 
    15191520            ap.getValueAsString( 'readws_lat' ), 
    15201521            ap.getValueAsString( 'readws_lon' ), 
     
    15231524            ap.getValueAsString( 'read_time' ), 
    15241525            ap.getValueAsString( 'read_length' ), 
     1526            ap.getValueAsString( 'readws_server' ), 
    15251527        ) 
    15261528        print "dbg: cmd", cmd 
Note: See TracChangeset for help on using the changeset viewer.