Computes the declination of the Sun for a given Julian Day.

declination(jd)

Arguments

jd

Julian Day.

Value

Declination in degrees and decimal fraction.

References

https://gml.noaa.gov/gmd/grad/solcalc/calcdetails.html

Meeus, J. 1999. Astronomical Algorithms. Willmann-Bell, Richmond, Virginia, USA.

Reda, I. and Andreas, A. 2003. Solar Position Algorithm for Solar Radiation Applications. 55 pp.; NREL Report No. TP-560-34302, Revised January 2008. https://www.nrel.gov/docs/fy08osti/34302.pdf

Author

Javier G. Corripio jgc@meteoexploration.com

Examples


if (FALSE) {
declination(120)

## Plot daily changes in declination
jdays=c(1:365)
plot(declination(jdays), xlab="days of the year",
     ylab="declination", type="l")
}