Pre-azimuthal fn generation

Current functionality works up to the interpolation of lat/lon data
points.

to do: get the azeqd function and inverse function working