% routine to compare ncep & rick wind
load c:\gbs99\model\taubuoy11.dat
r=taubuoy11;

load c:\gbs99\model\apr20.dat
j=apr20;

load c:\gbs99\model\EL9904at.dat
e=EL9904at;

%subsample ship
e=e(1:30:length(e),:);

%get u&v in oceanographic
[eu,ev]=sd2uv(e(:,5),e(:,6));

%filter ship?
feu=plfilt(eu);fev=plfilt(ev);

disp('Calculating stress...');
for i=1:length(eu);
   [stx(i),sty(i)]=wstr(-1*feu(i),-1*fev(i),15);
end;   

subplot(2,1,1);
plot(r(:,1),r(:,2),'r',j(:,1),j(:,17),'g',...
   e(:,4),stx,'b','linewidth',2);
legend('rick','jim','ship')
ylabel('taux')
title('BUOY 11');
axis([104 113 -.2 .4])
set(gca,'xticklabels',[]);
subplot(2,1,2);
plot(r(:,1),r(:,3),'r',j(:,1),j(:,18),'g',...
   e(:,4),sty,'b','linewidth',2);
ylabel('tauy')
axis([104 113 -.2 .1])
yd2calab(99)