function [index,distance]=nearxy(x,y,x0,y0,dist);
%FUNCTION NEARXY  finds the indices of the array 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(min(min(min(min(min(distance)))))));  % finds closest point
  index=index(1);
end
distance=distance(index);
