I'm not sure why you've added this. It's checking it doesn't blow up, but you don't check the resulting regridded surface values (values + 1 or linear_values + 1, right?). It would be good to check the extrapolation mode is correctly applied to the surface since you are using it via functools.partial (nifty!).