And do recall that I meant this comment to be for line 655:

    new_values = scipy.interpolate.griddata((lons, lats), values, zip(xx.flat, yy.flat))