function [w,wmean,t]=hodo(cdf,sta,tit)
%function [w,wmean,t]=hodo(cdf,sta,tit)
%  plots mean hodograph at a particular station 
%  sta is station number, tit is title for plot
if(nargin==1),
 sta=cdf;
 cdf='ecomts.cdf';
end
[name,size]=mcdinq(cdf,'time');
%i1=size-74;
%i2=73;
i1=0;
i2=size;
t=mcvgt(cdf,'time',i1,i2);
u=mcvgt(cdf,'u',[i1 (sta-1) 0],[i2 1 10]);
v=mcvgt(cdf,'v',[i1 (sta-1) 0],[i2 1 10]);
w=u+i*v;
%[w]=boxfilt(w,72);
%wm=mean(w);
%compass2(wm,'-');
for k=1:1:size,
  compass2(w(k,:),'-')
  text(.65,.85,['day ' num2str(t(k))],'sc')
  if nargin==3,
   title(tit)
  end
  pause
end

