Changeset 1080


Ignore:
Timestamp:
06.01.2016 15:25:49 (4 years ago)
Author:
klaus
Message:

add disp/vel simulations; fix norm bug

Location:
SHX/trunk/SeismicHandler
Files:
2 added
4 edited

Legend:

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

    r1063 r1080  
    137137            col = collist[colidx] 
    138138            if normtype.upper().startswith('S'): 
    139                 norm = 2.*max(ydata) 
     139                lmaxval = max( abs(max(ydata)), abs(min(ydata)) ) 
     140                norm = 2.*lmaxval 
    140141            else: 
    141142                norm = 2.*maxval 
     
    150151        # new tick labels (station names) at y-axis instead of numbers 
    151152        plt.yticks( labloc, labtxt ) 
    152         ax.set_ylim( 0., ypos ) 
     153        ax.set_ylim( 0., ypos-0.5 ) 
     154        print "dbg: ylim 0,", ypos 
    153155        if trclabsize: 
    154156            ax.tick_params(axis='y', which='major', labelsize=trclabsize ) 
  • SHX/trunk/SeismicHandler/commands/simulate.py

    r1064 r1080  
    120120            "lrsm-lp"       :  (0.004,  10.), 
    121121            "lrsm-sp"       :  (0.01, 1000.), 
     122            "velocity"      :  (0.01,  None), 
     123            "displacement"  :  (0.01,  None), 
    122124        } 
    123125 
  • SHX/trunk/SeismicHandler/modules/wx_.py

    r1073 r1080  
    10441044        Gather information about traces. Called on redraw (at idle state). 
    10451045        """ 
    1046         min = UTCDateTime() # now 
    1047         max = UTCDateTime(0) # 1970 
     1046        mint = UTCDateTime() # now 
     1047        maxt = UTCDateTime(0) # 1970 
    10481048        minrel = 1e20 
    10491049        maxrel = 0 
     
    10581058            if e > maxrel: 
    10591059                maxrel = e 
    1060             if t.stats.starttime < min: 
    1061                 min = t.stats.starttime 
    1062             if t.stats.endtime > max: 
    1063                 max = t.stats.endtime 
     1060            if t.stats.starttime < mint: 
     1061                mint = t.stats.starttime 
     1062            if t.stats.endtime > maxt: 
     1063                maxt = t.stats.endtime 
    10641064            # needed for relative time axis 
    10651065            l = e - s 
     
    10681068            # normation 
    10691069            # no time window set or normtype forces full traces 
    1070             if self._timewindow[1] is None or self._normtype.startswith("A"): 
     1070            #if self._timewindow[1] is None or self._normtype.startswith("A"): 
     1071            if self._timewindow[1] is None: 
    10711072                data = t.data 
    10721073            else: 
    10731074                data = t.get_datawindow(*self._timewindow)  #? KS 13-Aug-15 
    1074             t_ampl = abs(data.max()) 
     1075            t_ampl = max( abs(data.max()), abs(data.min()) ) 
    10751076            if t_ampl > amplitude: 
    10761077                amplitude = t_ampl 
    10771078 
    1078         self.start = min 
     1079        self.start = mint 
    10791080        self.relstart = minrel 
    1080         self.end = max 
     1081        self.end = maxt 
    10811082        self.relend = maxrel 
    10821083        self.maxDuration = maxDuration 
  • SHX/trunk/SeismicHandler/patches/obspy_.py

    r1070 r1080  
    413413    # much faster is 
    414414    for i in range(len(self.shx._parent.traces)): 
    415         #print "dbg:", i, self.shx._parent.traces[i].stats.sh.UNIQUE_ID, self.stats.sh.UNIQUE_ID 
    416415        if self.shx._parent.traces[i].stats.sh.UNIQUE_ID == self.stats.sh.UNIQUE_ID: 
    417416            idx = i 
    418417            break 
    419418    else: 
    420         #print "dbg: trace_index: this cannot happen" 
    421419        idx = 0 
    422420    if one: 
Note: See TracChangeset for help on using the changeset viewer.