lunar_altitude
double geo = lunar_altitude(t, latitude, longitude);
return (lunar_altitude(t, latitude, longitude) -
double lunar_altitude(double t, double latitude, double longitude);