Adding in extra interpolation lines for testing

Comparing delaunay, map_coordinates order 1, and map_coordinates order 3 interpolation times for x and y.
Playing around with interpolation schemes

Interpolating with Delaunay and map_coordinates, to find either xp,yp or xp,yp,zp. Fixed the indices going into the interpolation functions.
Added interpolate function to tools.py

Have moved interpolation functions into a separate function, interpolate, in tools.py which contains options for moving between grid space, projected xy space, and lon/lat. It also has functions for moving from grid space to xy or lon/lat (but not the other way since it depends on the uniformity of the grid).
Separated 2d and 3d interpolation functions

There are now tools.interpolate2d() and tools.interpolate3d() for transforming between coordinate systems.
Changed map_coordinates inputs to get fix problem

The horizontal interpolation through map_coordinates with mode='nearest' gave weird corner behavior, but setting mode='constant' with cval=np.nan is helping somewhat at least.
Cleaned up code

Cleaned up old things I tried in this branch in run.py