def processGeo(fname1, fname2): #Load the relevent files dat1 = np.loadtxt(str(fname1), skiprows=17) head1 = np.loadtxt(str(fname1) + '.h') dat2 = np.loadtxt(str(fname2), skiprows=17) head2 = np.loadtxt(str(fname2) + '.h') #assign real units to the seismograms - reshape dat1 = preprocessdata(dat1, head1) dat2 = preprocessdata(dat2, head2) #calculate the geometric mean for the two horizontal axes FAS1 = calcFAS(dat1, head1) FAS2 = calcFAS(dat2, head2) FAS = geoMean(FAS1, FAS2) return FAS #header just for sampling frequency.
def processNorm(fname): #Load the revelent files data = np.loadtxt(str(fname), skiprows=17) header = np.loadtxt(str(fname)+'.h') #assign real units to the seismograms - reshape data = preprocessdata(data, header) FAS = calcFAS(data, head1) return FAS