refraction
double refraction(double elevation);
refraction(loc->elevation) + moon_radius);
double alpha = refraction(loc->elevation) + sun_radius;
double alpha = refraction(loc->elevation) + sun_radius;