factor out main loop

Now it will be relatively straightforward to factor out
_interpolate_field to Cython or C.