#include <boost/python.hpp>
#include <GeographicLib/Geoid.hpp>

using namespace boost::python;
using namespace GeographicLib;

double EllipsoidHeight(Geoid& geoid,
 double lat, double lon, double hmsl) {
 return hmsl + Geoid::GEOIDTOELLIPSOID * geoid(lat, lon);
}

BOOST_PYTHON_MODULE(PyGeographicLib) {

 class_<Geoid, boost::noncopyable>("Geoid", init<std::string>())
 .def("EllipsoidHeight", &EllipsoidHeight,
 "Return geoid height:\n\
 input: lat, lon, height_above_geoid\n\
 output: height_above_ellipsoid")
 ;

}