SunPos

MATLAB

[RECF, rS, RtAsc, Dec] = SunPos(JD);

Description

Uses the low precision formulas for the sun presented in section C of The Astronomical Almanac for determining the sun position, right ascension and declination. The precision is given as one minute of arc, but the solar mean anomaly, longitude, etc. use the IAU 2000 methods, increasing the precision. Valid for years 1950-2050 AD.

Parameters

Parameters	Description
JD [vector]	Vector (N elements) of Julian Dates (days)

Return Values

Return Values	Description
RECF [Nx3 matrix]	Solar position in ECF coordinates (m)
rS [Nx3 matrix]	Solar position in ECI coordinates (m)
RtAsc [vector]	Solar Right Ascension (deg)
Dec [vector]	Solar Declination (deg)

References

The Astronomical Almanac (2013) U.S. Government Printing Office, 2012.
D. A. Vallado, Fundamentals of astrodynamics and applications. New York: McGraw-Hill, 1997.

Examples

>> JD = Gregorian2Julian('2-Apr-1994');

>> [RECF, rS, RtAsc, Dec] = SunPos(JD);