Index: /issm/trunk-jpl/src/m/interp/interp.py
===================================================================
--- /issm/trunk-jpl/src/m/interp/interp.py	(revision 17890)
+++ /issm/trunk-jpl/src/m/interp/interp.py	(revision 17891)
@@ -187,18 +187,17 @@
 	x,y:			data point coordinates
 	data:			data to be interpolated (same length as x,y)
-
-		function: form of radial basis function for interpolation:
-			'multiquadric': sqrt((r/self.epsilon)**2 + 1) (default)
-			'inverse': 1.0/sqrt((r/self.epsilon)**2 + 1)
-			'gaussian': exp(-(r/self.epsilon)**2)
-			'linear': r
-			'cubic': r**3
-			'quintic': r**5
-			'thin_plate': r**2 * log(r)
-		epsilon: adjustable constant for scaling radial distance.  Defaults to 
+	function:	form of radial basis function for interpolation:
+					'multiquadric': sqrt((r/self.epsilon)**2 + 1) (default)
+					'inverse': 1.0/sqrt((r/self.epsilon)**2 + 1)
+					'gaussian': exp(-(r/self.epsilon)**2)
+					'linear': r
+					'cubic': r**3
+					'quintic': r**5
+					'thin_plate': r**2 * log(r)
+	epsilon:		adjustable constant for scaling radial distance.  Defaults to 
 					approximate average distance between nodes.
-		smooth: float>0, adjusts the amount of smoothing applied.  Defaults to 0,
+	smooth:		float>0, adjusts the amount of smoothing applied.  Defaults to 0,
 					such that the function always passes through nodal points.
-		z:	coordinate array if interpolating in 3 dimensions
+	z:				coordinate array if interpolating in 3 dimensions
 
 	Usage:
@@ -211,5 +210,5 @@
 
 	# unpack kwargs
-	function=kwargs.pop('function','multiquadric')
+	function=kwargs.pop('function','gaussian')
 	if 'function' in kwargs: del kwargs['function']
 	epsilon=kwargs.pop('epsilon',None)