solar_longitude_atafter
double t = solar_longitude_atafter(lambda, rd) + zone;
t_u = solar_longitude_atafter(jq2->longitude, floor(t_u));
t_jq = solar_longitude_atafter(lambda, floor(t_jq)) + zone;
t = solar_longitude_atafter(longitude, approx);
t = solar_longitude_atafter(lambda, day_approx) + loc->zone;
double solar_longitude_atafter(double lambda, double t);