function drifters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all; startup

%
% Load model results
%

tallfigure;

gridname='../MESH/bank150';
gridstruct=loadgrid(gridname);

%

PTHFILE=['../fcast_archives/FCAST_1.0_djm/EN302_FC1/CRUISE_Q4.3/DRIFTERS']

[gridname,ndrog,ndts,tsec,pth]=read_pth(PTHFILE,4);

size(pth);
ndts=ans(1)/ndrog;
xdr=pth(:,1);
xdr=reshape(xdr,ndrog,ndts);
XO=xdr';
ydr=pth(:,2);
ydr=reshape(ydr,ndrog,ndts);
YO=ydr';

PTHFILE=['../fcast_archives/FCAST_1.0_djm/EN302_FC1a/CRUISE_Q4.1/DRIFTERS']

[gridname,ndrog,ndts,tsec,pth]=read_pth(PTHFILE,4);

size(pth);
ndts=ans(1)/ndrog;
xdr=pth(:,1);
xdr=reshape(xdr,ndrog,ndts);
xdr=xdr';
ydr=pth(:,2);
ydr=reshape(ydr,ndrog,ndts);
ydr=ydr';


subplot(2,1,1)
plotbnd(gridstruct);
axis([160000 230000 -230000 -160000]);
axis('square');set(gca,'XTick',[]);set(gca,'YTick',[]);
hold on

idr=2;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

idr=5;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

idr=6;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

title(['EN302_FC1ar7 Test: 13m Drifter Trajectories'])

subplot(2,1,2)
plotbnd(gridstruct);
axis([160000 230000 -230000 -160000]);
axis('square');set(gca,'XTick',[]);set(gca,'YTick',[]);
hold on

idr=1;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

idr=3;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

idr=4;
  hp = plot(XO(:,idr),YO(:,idr),'b-');
  set(hp,'LineWidth',4.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'r-');

title(['EN302_FC1_bc_a Test: 33m Drifter Trajectories'])
