% routine for extracting/subsampling bathy data

%ph=menu('Which database?','oreilly','candela');
ph=2;

if ph==1;

  load /data5/jmanning/bathy/necs_60m
  for i=1:100;
    x=eval(['s',num2str(i),'(:,1)']);
    y=eval(['s',num2str(i),'(:,2)']);
    id=find(x>-69&x<-67&y>40&y<42);
    if id~=[];
      plot(x,y);
      title(num2str(i));
      pause(1);
    end;
%   hold on;
  end;

elseif ph==2;
% try something else
  load map
  for i=1:5;
    eval(['x=isobath60m',num2str(i),'(:,4);']);
    eval(['y=isobath60m',num2str(i),'(:,3);']);
   id=find(x>-69&x<-66&y>40&y<41.5);
   if id~=[];
     plot(x(id),y(id));
     title(num2str(i));
     lat=y(id);lon=x(id);
     save isobath60 lat lon
    
   end;
%   hold on?;
  end;
end;

