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

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

%load 10-m data
load gb10mallA
t=z(:,1);
% load 40-m data
load gb40m99
td=zd(:,1);

load windstress99

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

V=[-71 -65 39 43];

% start loop
err=0.001;
for n=1:N
   driftID=[];
   speed=ones(1,4);
   dspeed=speed;
   clf
   gbmapss2(V)
   aspect5(41);
   axis(V);
   hold on
   % extract all drifter data for time window
   k=find(t>=t1(n)-err & t<=t2(n)+err);
   zn=z(k,:);
   kd=find(td>=t1(n)-err & td<=t2(n)+err);
   zdn=zd(kd,:);
   znew=reconfig1(zn);
   zdnew=reconfig1(zdn);
   % plot all 10-m 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*')
   % denote each day on track
   plotday=1;
   if plotday==1
      [xp,yp,tp]=get_day_xyt(znew);
      plot(xp,yp,'r+')
      hold on
   end
   text(znew(kk-1,2),znew(kk-1,3)+.05,int2str(znew(kk-1,8)))
   title([int2str(dd2yd(t1(n))), ' - ',int2str(dd2yd(t2(n)))])
   % get drifter id and speed
   did=sort(znew(kk-1,8));
   knan=find(isnan(did));
   did(knan)=[];
   Ndid=length(did);
   GSP=0;
   if GSP==1,
      for ndid=1:Ndid;
      m=find(znew(:,8)==did(ndid));
      spd=sqrt(znew(m,4).^2+znew(m,5).^2);
      spmax=max(spd);
      spmean=mean(spd);
      maxtim=find(spd==spmax);
      timmax=dd2yd(znew(m(maxtim),1));
      speed=[speed;did(ndid) spmean spmax timmax];
      end
   end
   % disp(speed)
      
   % plot all 40-m drifters   
   kk=find(isnan(zdnew(:,1)));
   h=plot(zdnew(:,2),zdnew(:,3),'r');
   set(h,'linewidth',2);
   hold on
   plot(zdnew(kk-1,2),zdnew(kk-1,3),'ro')
   % denote each day on track
   plotday=1;
   if plotday==1
      [xp,yp,tp]=get_day_xyt(zdnew);
      plot(xp,yp,'r+')
   end
   text(zdnew(kk-1,2),zdnew(kk-1,3)+.05,int2str(zdnew(kk-1,8)))
   % get drifter id and speed
   ddid=sort(zdnew(kk-1,8));
   dknan=find(isnan(ddid));
   ddid(dknan)=[];
   Nddid=length(ddid);
   GSP=0;
   if GSP==1,
      for nddid=1:Nddid;
      m=find(zdnew(:,8)==ddid(nddid));
      spd=sqrt(zdnew(m,4).^2+zdnew(m,5).^2);
      spmax=max(spd);
      spmean=mean(spd);
      maxtim=find(spd==spmax);
      timmax=dd2yd(zdnew(m(maxtim),1));
      dspeed=[dspeed;ddid(nddid) spmean spmax timmax];
      end
   end
   % disp(dspeed)
      
   % plot mean wind stress
   err=0.01;
   K=find(dd2yd(t1(n))-err<=t99 & t99<=dd2yd(t2(n))+err);
   txm=mean(tx99(K));
   tym=mean(ty99(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);
   pause
   % plot wind stress scale
   xst1=-66.5;yst1=40.2;
   KK=2;
   KXY=.755;
   yst2=yst1;
   xst2=xst1+KK.*.2./KXY;
   h=plot([xst1 xst2],[yst1 yst2],'g',xst1,yst1,'g.',xst2,yst2,'g*');set(h,'linewidth',2);
end

