load test_data
[P,F]=spectrum(u(:,1),v(:,1),1024,512,1024,1/min(diff(jd)));
flim=F([[1:2:20] [25:5:80] [90:20:513]]);
ss=specav(P,F,flim,0,.9);
loglog(ss.Freq.central,[ss.Pxx.est ss.Pyy.est])
