function [index,distance]=nearxy(x,y,x0,y0,dist);
% NEARXY  finds the indices of (x,y) that are closest to the point (x0,y0).
%        [index,distance]=nearxy(x,y,x0,y0) finds the closest point and
%                                           the distance
%        [index,distance]=nearxy(x,y,x0,y0,dist) finds all points closer than
%                                           the value of dist.
% rsignell@crusty.er.usgs.gov
%
distance=sqrt((x-x0).^2+(y-y0).^2);
if (nargin > 4),
  index=find(distance<=dist);     %finds points closer than dist
else,
  index=find(distance==min(distance));  % finds closest point
  index=index(1);
end
distance=distance(index);