double calcSunApparentLong(double t) { double o = calcSunTrueLong(t); double omega = 125.04 - 1934.136 * t; double lambda = o - 0.00569 - 0.00478 * sin(degToRad(omega)); return lambda; // in degrees }
/* Purpose: calculate the apparent longitude of the sun (degrees) */ double calcSunApparentLong(double t) { return calcSunTrueLong(t) - 0.00569 - 0.00478 * sin(deg2rad(125.04-1934.136*t)); }