* a bit of rearanging routines to get a more logical order
  * new intermediate function: ckdtree_points to calculate grid points
  * new function: calc_north_pole_interp_grid_points for a more
    appropriate grid for polar projections
  * new functions: spherical_to_cartesian and cartesian_to_spherical
  * splitting: mask_big_triangles