Geodesic routines for modern Fortran
