function gb_plot96(tstart,tstop)
% gb_plot96(tstart,tstop) plots tracks of all drifters 
% on bank during the time between tstart and tstop
% 4/2/01

% convert time to drifter time
tstart=yd2dd(tstart,1996);
tstop=yd2dd(tstop,1996);

%load data
load gb10mallA
t=z(:,1);
load windstress96
tws=yd2dd(t96,1996);

figure(2)
clf
plot(t96,tau96)
hold on
figure(1)

% set start, stop times
dt=5;
t1=[tstart:dt:tstop-dt];
t2=t1+dt;
N=length(t1);

V=[-69 -65.5 40 42.5];

% start loop
clf
for n=1:N
   clf
   gbmapss2(V);
   aspect5(41);
   axis(V);
   hold on
   % extract all drifter data for time window
   k=find(t>=t1(n) & t<=t2(n));
   zn=z(k,:);
   znew=reconfig1(zn);
   % plot all drifters
   kk=find(isnan(znew(:,1)));
   h=plot(znew(:,2),znew(:,3));set(h,'linewidth',2);
   plot(znew(kk-1,2),znew(kk-1,3),'r*')
   % label period
   pdate=[dd2cd1(t1(n)),' - ',dd2cd1(t2(n))];
   title(pdate);
      
   % plot mean wind stress
   err=0.01;
   K=find(dd2yd(t1(n))-err<=t96 & t96<=dd2yd(t2(n))+err);
   txm=mean(tx96(K));
   tym=mean(ty96(K));
   xst1=-66.75;yst1=41.75;
   KK=2;
   KXY=.755;
   yst2=yst1+KK.*tym;
   xst2=xst1+KK.*txm./KXY;
   h=plot([xst1 xst2],[yst1 yst2],'k',xst2,yst2,'k*');
   set(h,'linewidth',2);
   hold on
   pause
end

