Sinusodial sunset

To within 2 minutes of accuracy, Albuquerque sunset for 2012 is given by:

where x is the day of the year (1 = January 1st, 366 = December 31st, since this is a leap year), and the result is the number of decimal hours from midnight (ie, 17 means 5pm)

This isn’t entirely surprising, since sunset just depends on the time of true noon, which is given by the equation of time, which is itself the sum of two sine waves, and by the length of day, which also appears to be sinusoidal.

I’m guessing this works for any city, even those north of the Arctic circle, providing you allow values > 24 (meaning no sunset) and < 0 (meaning no sunrise), but do the Fourier analysis only for days where sunset actually occurs.

This entry was posted in Programming. Bookmark the permalink.

Leave a Reply

Your email address will not be published. Required fields are marked *