function drifters(runid,isim)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Load observational data
%
OPATH='../OBS/'
UVOF=[OPATH,'DRIFTERS.m3d']
%
[UVO,gridname,year,ncol,gmt]=read_obs(UVOF);
no=length(UVO(:,1));
ndr=6;ndts=no/ndr;
TO=(reshape(UVO(:,1),ndr,ndts))';
XO=(reshape(UVO(:,2),ndr,ndts))';
YO=(reshape(UVO(:,3),ndr,ndts))';
UO=(reshape(UVO(:,5),ndr,ndts))';
VO=(reshape(UVO(:,6),ndr,ndts))';

%
% Load model results
%

tallfigure;

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

%

PTHFILE=['../fcast_archives/',runid,'/CRUISE_Q4.',num2str(isim),'/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',2.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'b-');

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

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

title([runid,' 13m Drifter Trajectories ',isim])

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',2.0)
  hp = plot(xdr(:,idr),ydr(:,idr),'b-');

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

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

title([runid,' 33m Drifter Trajectories ',isim])
