function gb_plot97(tstart,tstop)
% gb_plot97(tstart,tstop) plots tracks of all drifters 
% on bank during the time between tstart and tstop
% 3/25/01; 4/1/01

%load data
load gb10mallA
t=z(:,1);
load windstress97
tau=sqrt(tx97.^2+ty97.^2);
mws=find(t97>=tstart & t97<=tstop);
PS=1;
% plot stress
if PS==1
   figure(2)
   clf
   plot(t97(mws),tau(mws));
   hold on
end

% convert to 1995 drifter days
tstart=yd2dd(tstart,1997);
tstop=yd2dd(tstop,1997);

figure(1)
clf

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

% V=[-69 -65.5 40 42.5];
V=[-70 -65.5 39.75 42.5];

% start loop
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,:);
   % find all drifters inside map domain
   k=find_inside(zn,V);
   znew=reconfig1(zn(k,:));
   
   % specify which drifters NOT to plot
   RD=0;
   if RD==1;
      DP=[72,76,86,74,84];
      nDP=length(DP);
      for ndp=1:nDP
         nnn=find(znew(:,8)==DP(ndp));
         znew(nnn,:)=[];
      end
   end
   
   % plot all drifters
   kk=find(isnan(znew(:,1)));
   h=plot(znew(:,2),znew(:,3));set(h,'linewidth',2);
   hold on
   plot(znew(kk-1,2),znew(kk-1,3),'r*')
   text(znew(kk-1,2),znew(kk-1,3)+.1,int2str(znew(kk-1,8)))
   title([int2str(dd2yd(t1(n))), ' - ',int2str(dd2yd(t2(n)))])
      
   % plot mean wind stress
   err=0.01;
   K=find(dd2yd(t1(n))-err<=t97 & t97<=dd2yd(t2(n))+err);
   txm=mean(tx97(K));
   tym=mean(ty97(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],'g',xst2,yst2,'g*');set(h,'linewidth',2);
   taum=sqrt(txm.^2+tym.^2);
   taudr=180.*atan2(tym,txm)./pi;
   % disp([taum  taudr])
   pause
end
