H=zeros(mmm,'f')
tpi=2*pi
alpha=linspace(0,tpi,mmm)
ae=sqrt(8./3.)
H=ae*0.5*(1.-cos(alpha))
return H
########################################################################
print " The input file must have two columns: time(sec) & amplitude "
a,b,num =read_two_columns()
sr,dt,mean,sd,rms,skew,kurtosis,dur=signal_stats(a,b,num)
sr,dt=sample_rate_check(a,b,num,sr,dt)
########################################################################
print " "
print " Remove mean: 1=yes 2=no "
mr_choice = GetInteger2()
print " "
print " Select Window: 1=Rectangular 2=Hanning "
h_choice = GetInteger2()
########################################################################
from scipy.stats.stats import pearsonr
import numpy as np
########################################################################
print ("Each file must have two columns: time(sec) & accel(G)")
a1,b1,num1 =read_two_columns()
sr1,dt1,mean1,sd1,rms1,skew1,kurtosis1,dur1=signal_stats(a1,b1,num1)
sr1,dt1=sample_rate_check(a1,b1,num1,sr1,dt1)
########################################################################
a2,b2,num2 =read_two_columns()
sr2,dt2,mean2,sd2,rms2,skew2,kurtosis2,dur2=signal_stats(a2,b2,num2)
sr2,dt2=sample_rate_check(a2,b2,num2,sr2,dt2)
########################################################################
a1=array(a1)
b1=array(b1)
#print ("a1 time: ", a1)
