Index: /issm/trunk-jpl/src/m/interp/interp.py
===================================================================
--- /issm/trunk-jpl/src/m/interp/interp.py	(revision 17891)
+++ /issm/trunk-jpl/src/m/interp/interp.py	(revision 17892)
@@ -178,5 +178,5 @@
 	return interpdata
 #}}}
-def RadialInterp(x,y,data,**kwargs):#{{{
+def RadialInterp(x,y,data,xi,yi,**kwargs):#{{{
 	'''
 	Interpolation using a radial basis function in 2 or 3 dimensions.
@@ -187,4 +187,5 @@
 	x,y:			data point coordinates
 	data:			data to be interpolated (same length as x,y)
+	xi,yi:		coordinates to interpolate onto
 	function:	form of radial basis function for interpolation:
 					'multiquadric': sqrt((r/self.epsilon)**2 + 1) (default)
@@ -200,4 +201,5 @@
 					such that the function always passes through nodal points.
 	z:				coordinate array if interpolating in 3 dimensions
+	zi:			coordinate array if interpolating in 3 dimensions
 
 	Usage:
@@ -225,5 +227,5 @@
 		else:
 			rbfi=Rbf(x,y,z,data,function=function,smooth=smooth)
-		interpdata=rbfi(x,y,z)
+		interpdata=rbfi(xi,yi,zi)
 	else:
 		if epsilon:
@@ -231,5 +233,5 @@
 		else:
 			rbfi=Rbf(x,y,data,function=function,smooth=smooth)
-		interpdata=rbfi(x,y)
+		interpdata=rbfi(xi,yi)
 	
 	return interpdata