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
