seismograms.py 文件源码

def read(self, station_name, plot=True):

          seismogram = np.loadtxt(station_name) 
          self.t = seismogram[:,0]
          self.trace_x = seismogram[:,1]
          self.nt = len(self.t)
          self.dt = np.max(self.t)/len(self.t)

          #if (plot == true):

          plt.plot(self.t, self.trace_x)

 ################################################################
 #              cross correlate travel times
 ################################################################