Index: /issm/trunk-jpl/src/m/mesh/meshintersect.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/meshintersect.m	(revision 20126)
+++ /issm/trunk-jpl/src/m/mesh/meshintersect.m	(revision 20126)
@@ -0,0 +1,22 @@
+function indices=meshintersect(lat,long,lats,longs,varargin)
+%MESHINTERSECT - return indices (into lat and long) of common values between (lat,long) and (lats,longs). 
+%  i.e: lat(index)=lats; long(index)=longs;
+%
+%   Usage:
+%      index=meshintersect(md.mesh.lat,md.mesh.long,mdsmaller.mesh.lat,mdsmaller.mesh.long);
+%      index=meshintersect(md.mesh.lat,md.mesh.long,mdsmaller.mesh.lat,mdsmaller.mesh.long,'tolerance',1e-10); %within a certain tolerance.
+
+
+	%process options: 
+	options=pairoptions(varargin{:});
+
+	%retrieve tolerance: 
+	tolerance=getfieldvalue(options,'tolerance',1e-10);
+
+	%go through lats,longs and find within tolerance, the index of the corresponding value in lat,long: 
+	indices=zeros(length(lats),1);
+	
+	for i=1:length(lats),
+		distance=sqrt((lat-lats(i)).^2+(long-longs(i)).^2);
+		indices(i)=find(distance<tolerance);
+	end