Changeset 17892

Timestamp:

04/30/14 12:35:26 (11 years ago)

Author:

cborstad

Message:

CHG: generalize to allow interpolation onto arbitrary x,y coords

File:

• issm/trunk-jpl/src/m/interp/interp.py (modified) (5 diffs)

issm/trunk-jpl/src/m/interp/interp.py

@@ -178 +178 @@
         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 +187 @@
         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 +201 @@
                                         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 +227 @@
                 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 +233 @@
                 else:
                         rbfi=Rbf(x,y,data,function=function,smooth=smooth)
-                interpdata=rbfi(x,y)
+                interpdata=rbfi(xi,yi)
         
         return interpdata