A 2021 publication presents a method that uses a solar azimuth formula based on the
subsolar point and the
atan2 function, as defined in
Fortran 90, that gives an unambiguous solution without the need for circumstantial treatment. The subsolar point is the point on the surface of the Earth where the Sun is overhead. The method first calculates the
declination of the Sun and
equation of time using equations from The Astronomical Almanac, then it gives the x-, y- and z-components of the unit vector pointing toward the Sun, through
vector analysis rather than
spherical trigonometry, as follows: :\begin{align} \phi_{s} &= \delta, \\ \lambda_{s} &= -15(T_\mathrm{GMT}-12+E_\mathrm{min}/60), \\ S_{x} &= \cos \phi_{s} \sin (\lambda_{s}-\lambda_{o}), \\ S_{y} &= \cos \phi_{o} \sin \phi_{s} - \sin \phi_{o} \cos \phi_{s} \cos (\lambda_{s}-\lambda_{o}), \\ S_{z} &= \sin \phi_{o} \sin \phi_{s} + \cos \phi_{o} \cos \phi_{s} \cos (\lambda_{s}-\lambda_{o}). \end{align} where • \delta is the declination of the Sun, • \phi_{s} is the latitude of the subsolar point, • \lambda_{s} is the longitude of the subsolar point, • T_\mathrm{GMT} is the Greenwich Mean Time or UTC, • E_\mathrm{min} is the
equation of time in minutes, • \phi_{o} is the latitude of the observer, • \lambda_{o} is the longitude of the observer, • S_{x}, S_{y}, S_{z} are the x-, y- and z-components, respectively, of the unit vector pointing toward the Sun. The x-, y- and z-axises of the coordinate system point to East, North and upward, respectively. It can be shown that S_{x}^{2}+S_{y}^{2}+S_{z}^{2}=1. With the above mathematical setup, the solar zenith angle and solar azimuth angle are simply :Z = \mathrm{acos}(S_{z}), :\gamma_{s} = \mathrm{atan2}(-S_{x}, -S_{y}). (South-Clockwise Convention) where • Z is the solar zenith angle, • \gamma_{s} is the solar azimuth angle following the South-Clockwise Convention. If one prefers North-Clockwise Convention, or East-Counterclockwise Convention, the formulas are :\gamma_{s} = \mathrm{atan2}(S_{x}, S_{y}), (North-Clockwise Convention) :\gamma_{s} = \mathrm{atan2}(S_{y}, S_{x}). (East-Counterclockwise Convention) Finally, the values of S_{x}, S_{y} and S_{z} at 1-hour step for an entire year can be presented in a 3D plot of "wreath of
analemmas" as a graphic depiction of all possible positions of the Sun in terms of solar zenith angle and solar azimuth angle for any given location. Refer to
sun path for similar plots for other locations. == See also ==