if yr==82;
   %  load /home3/ocn/jmanning/wcr/position/position_all.dat;

   load e:\wcr\position\position_all.dat;
elseif yr==97;
  load /home3/ocn/jmanning/atsee/ronsta.dat;
  position_all=[ronsta(:,3) ronsta(:,1:2) ronsta(:,5) ronsta(:,1:2)];
end;

idda=find(position_all(:,1)>=3);
ipps=position_all(idda,:);
disp('posts mooring position assuming first sta ref')
[ipr,ipc]=size(ipps);
for jk=1:ipr;
  dists(jk)=dist(q(1,1),q(1,2),ipps(jk,7),ipps(jk,8));
end;

if ref==5;
  disp('Hardcoded to compute x-distance from the 60m isobath');
  load isobath60;
  reflat=interp1(lon60,lat60,q(1,2));
  dd=dist(reflat,q(1,2),q(1,1),q(1,2));
  if reflat<q(1,1);dd=-dd;end;
  cumdist=cumdist+dd;
end;

hold on;
sizesym(dists,-1*ipps(:,5),-1*ones(length(dists),1),3,position_all(idda,1),'r')
