Dear all, There is a problem with scipy.interpolate: In [2]: interpolate.griddata([-1.5,-1.0], [5,6],[[-1.12]]) Out[2]: array([[ 5.76]]) In [3]: interpolate.griddata([-1.0,-1.5], [5,6],[[-1.12]]) Out[3]: array([[ nan]]) It depends on the order of the input values. Cheers Wolfgang