load BP;
H=H-mean(H);
om=24*2*pi/12.4206;
j=sqrt(-1);
X=2*H.*exp(j*om*tm);
om=2*pi/48;
ns=71;
[tlow dmod xx]=lowpass(tm,X,om,ns);
amp=abs(dmod);
hold on;
plot(tm,H,'b', tlow,amp,'r');

om=2*pi/32
ns=71;
[tlow,hlow,hhi]=lowpass(tm,H,om,ns);
tm=tlow;
save eta_BP tm  amp hlow;