if(1),
nfft = 2000;
%adv = a244_3;
fs = 25;
end

if(0),
nfft =400;
adv = a244_1;
fs = 10;
end

for bno=5,
  tt = sprintf('File: %s, burst: %d',adv.file,bno)
  ii = adv.bur_ind(bno):adv.bur_ind(bno)+adv.samp_per_burst-1;
  u = adv.vx(ii);
  v = adv.vy(ii);
  w = adv.vz(ii);
  hdg = median(adv.hdg(ii)).*ones(adv.samp_per_burst,1);
  pitch = median(adv.pitch(ii)).*ones(adv.samp_per_burst,1);
  roll = median(adv.roll(ii)).*ones(adv.samp_per_burst,1);
  [ve,vn,vu]=xyz2enu(u,v,w,adv.hdg,adv.pitch,adv.roll,adv.samp_per_burst);

  figure(1),
  clf
  [sd1,az1,sd2,az2]=pcastats(ve,vn,30,'.k');  
  [s,d]=pcoord(ve,vn);
  dr = d-az1;
  dr = dr+(d<0)*360;
  [ur,vr]=xycoord(s,dr);
  t=((1:adv.samp_per_burst)./fs)';
  
  subplot(4,1,1)
  plot(t,vr,'b')
  axis([0 900 -50 50])
  set(gca,'xticklabel','')
  title(tt)

  subplot(4,1,2)
  plot(t,ur,'r')
  axis([0 900 -50 50])
  V=cmgspectra(vr,nfft,nfft/2,hanning(nfft),fs,.9,'linear');
  U=cmgspectra(ur,nfft,nfft/2,hanning(nfft),fs,.9,'linear');
  W=cmgspectra(vu,nfft,nfft/2,hanning(nfft),fs,.9,'linear');
  [k,xdecay]=wave_decay(U.freq, 10);

  subplot(2,2,3)
  loglog(V.freq,V.spec(:,1),'b');
  hold on
  loglog(U.freq,U.spec(:,1),'r');
  loglog(W.freq,W.spec(:,1),'c');
  if(fs==10),
    P = cmgspectra(100*adv.p(ii),nfft,nfft/2,hanning(nfft),fs,.9, ...
		   'linear');
    loglog(P.freq,P.spec(:,1),'k');
  end
  axis([.01 20 1e-5 100])
  plot(U.freq,xdecay,'--k')
  plot(U.freq,0.01*exp( (-5/3)*log(U.freq) ),'--k') 
  
end