Changeset 1200 for SHX


Ignore:
Timestamp:
15.12.2016 18:01:35 (3 years ago)
Author:
klaus
Message:

changes/bugfixes for magnitude plots

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

Legend:

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

    r1193 r1200  
    9494       Does not work with <circledef> parameter. 
    9595 
     96     * /excludenets[[br]] 
     97       Exclude nets given in this comma separated list. Used only in bulk 
     98       requests. 
     99 
    96100    === example === 
    97101 
     
    141145        #"SAVE", 
    142146        "BULKLINE", 
     147        "EXCLUDENETS", 
    143148    ] 
    144149 
     
    155160        else: 
    156161            fdsnadr = Settings.config.fdsnws['default_data_source'][0] 
     162        if self.qualifiers["EXCLUDENETS"]: 
     163            excludenets = self.qualifiers["EXCLUDENETS"].split(',') 
     164        else: 
     165            excludenets = [] 
    157166         
    158167        if len(self.parameters) == 1: 
     
    200209                        print "error parsing netsta '%s'" % netsta 
    201210                        continue 
     211                    if excludenets and net in excludenets: 
     212                        continue 
    202213                    bulk += "%s\n" % self.getBulkLine( net, station, location, 
    203214                        channel, reqtime, reqlength ) 
     
    211222                    except: 
    212223                        print "error parsing netsta '%s'" % netsta 
     224                        continue 
     225                    if excludenets and net in excludenets: 
    213226                        continue 
    214227                    bulk += "%s\n" % self.getBulkLine( net, station, loc, 
  • SHX/trunk/SeismicHandler/commands/magnitude.py

    r1191 r1200  
    4747    or[[br]] 
    4848    '''command''': MAGNITUDE SET <magtype> <station> <value> 
     49    or[[br]] 
     50    '''command''': MAGNITUDE PLOT <magtype> [<fname>] 
    4951 
    5052    Magnitude management. Magnitudes are stored in a separate container, which is controlled via this command set. 
     
    201203        elif subcmd == 'plot': 
    202204            magtype = self.getMagtype( self.parameters[1] ) 
     205            if len(self.parameters) > 2: 
     206                picfile = self.parameters[2].lower() 
     207            else: 
     208                picfile = None 
    203209            meanmag = maglist.meanMagnitude( magtype ) 
    204210            if meanmag == None: 
     
    217223            dists, mags, labels = self.getMagnitudePlotData( magtype ) 
    218224            ax.scatter( dists, mags ) 
     225            if min(mags) < 0.: 
     226                ax.set_ylim( 0., 1.1*max(mags) ) 
    219227            for lab,x,y in zip(labels,dists,mags): 
    220228                plt.annotate( lab, xy=(x,y), fontsize=9 ) 
    221229            plt.draw() 
    222             picfile = self.tmpPlotfile() 
     230            showit = False 
     231            if picfile == None: 
     232                showit = True 
     233                picfile = self.tmpPlotfile() 
    223234            plt.savefig( picfile, facecolor=fig.get_facecolor() ) 
    224             try: 
    225                 dspprog = Settings.config.misc.display_prog 
    226             except: 
    227                 dspprog = 'display' 
    228             os.system( "%s %s &" % (dspprog,picfile) ) 
     235            if showit: 
     236                try: 
     237                    dspprog = Settings.config.misc.display_prog 
     238                except: 
     239                    dspprog = 'display' 
     240                os.system( "%s %s &" % (dspprog,picfile) ) 
    229241        elif subcmd == 'determine': 
    230242            ap = AnalysisPar() 
  • SHX/trunk/SeismicHandler/commands/meta.py

    r1199 r1200  
    248248            if (oldgain == None or oldgain == 1.) and newgain != 1.: 
    249249                # multiply new gain onto trace data 
    250                 trc.data *= newgain 
     250                if newgain != None: 
     251                    trc.data *= newgain 
    251252            trc.set_info( 'METASTATUS', metstat ) 
    252253     
Note: See TracChangeset for help on using the changeset viewer.