solar_altitude
double alt = solar_altitude(t_u, loc->latitude, loc->longitude);
double solar_altitude(double t, double latitude, double longitude);