Index: /issm/trunk-jpl/src/m/interp/interp.py
===================================================================
--- /issm/trunk-jpl/src/m/interp/interp.py	(revision 18108)
+++ /issm/trunk-jpl/src/m/interp/interp.py	(revision 18109)
@@ -45,6 +45,8 @@
 	# create sub-vectors that just cover the limits of xi and yi
 	# TODO x,y not necessarily a grid, so need a better definition of dx,dy (e.g. average element size)
-	dx=x[1]-x[0]
-	dy=y[1]-y[0]
+	dx=500
+	dy=500
+	#dx=x[1]-x[0]
+	#dy=y[1]-y[0]
 	xlim=[min(xi)-dx,max(xi)+dx]
 	ylim=[min(yi)-dy,max(yi)+dy]
@@ -76,5 +78,5 @@
 		xyinterp=npy.dstack([xi,yi])[0]
 		xg,yg=npy.meshgrid(subx,suby)
-		xydata=npy.dstack([xg.ravel(),yg.ravel()])[0]
+		xydata=npy.dstack([subx,suby])[0]
 		tree=cKDTree(xydata)
 		nearest=tree.query(xyinterp)[1]