Respuesta: brújula solar
estoy modificando el programa, para la declinación solar utilizo la formula de Spencer
$declinacion = 0.006918 - (0.3999 * cos($ang)) + (0.070257 * sin($ang)) - (0.006758 * cos(2 * $ang)) + (0.000907 * sin(2 * $ang)) - (0.002697 * cos(3 * $ang)) + (0.00148 * sin(3 * $ang))
donde $ang = (2 * $pi * ($dia-1)/365)
sin embargo el valor de esta declinación es diferente del que se puede obtener en almanaques, creo que es por ser la formula "antigua" de 1971 y al pasar los años el error se ha hecho mas grande, conocéis alguna formula para la declinación solar que este mas actualizada?
Esto ya no es lo que era...