%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function drifter_composite(runid,iter)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%
% Load model results
%

PTHFILE=['../fcast_archives/',runid,'/CRUISE_Q4.',...
          num2str(iter),'/DRIFTERS']

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

b=loadgrid('../MESH/bank150');
[la,lo]=xy2ll(b.x,b.y);
b.x=lo;b.y=la;

[gridname,ndrog,ndts,tsec,pth]=read_pth(PTHFILE,4);
%Convert to lo/la
[la,lo]=xy2ll(pth(:,1),pth(:,2));
pth(:,1)=lo;pth(:,2)=la;

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';

figure
hb=plotbnd(b);
axis('tight')
set(hb,'LineWidth',2)

hold on

for idr=1:ndrog
  hp = plot(xdr(:,idr),ydr(:,idr),'b-');
end

title([runid,' Iter ',num2str(iter),' Drifter Trajectories'])
hc=lcontour(b,'z',[30 45 60 75 90 120]);
grid

hold off
