source: issm/oecreview/Archive/19101-20495/ISSM-20121-20122.diff

Index: ../trunk-jpl/src/m/geometry/locationtonode.m
===================================================================
--- ../trunk-jpl/src/m/geometry/locationtonode.m        (revision 20121)
+++ ../trunk-jpl/src/m/geometry/locationtonode.m        (revision 20122)
@@ -1,11 +1,26 @@
-function node=locationtonode(md,location)
+function nodes=locationtonode(md,location,varargin)
 %LOCATIONTONODE - find, given a string location (ex: 'LA', 'new york', the nearest node on a mesh3dsurface.
 %
 %   Usage:
 %      node=locationnode(md,'LA');
+%      nodes=locationnode(md,'LA',100); %option, specify a distance (km) around the location that will return nodes
 %
 %   See also: geoCode (in externalpackages), find_point
 
-coords=geoCode(location,'osm');  
-latny=coords(1); longny=coords(2);
-node=find_point(md.mesh.lat,md.mesh.long,latny,longny);
+       if nargin==3,
+               radius=varargin{1};
+       else
+               radius=0;
+       end
+
+       coords=geoCode(location,'osm');  
+       latny=coords(1); longny=coords(2);
+       node=find_point(md.mesh.lat,md.mesh.long,latny,longny);
+
+       if radius>0,
+               distance=sqrt( (md.mesh.x-md.mesh.x(node)).^2 + (md.mesh.y-md.mesh.y(node)).^2 + (md.mesh.z-md.mesh.z(node)).^2);
+               nodes=find(distance<radius*1000);
+       else
+               nodes=node;
+       end
+end