added interp_polar function

Function for angular interpolation accepting as input any range of angles.
You have indices sorting the array from the previous line. Sorting is O(n*log n), reordering with an index array is O(n), so go with the second, I think something like `sdist = dist[np.arange(dist.shape[0])[:, None], asort]` should get you there.