![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_01.png\")
<\/p>\nHappy Crepusculus! Never heard of Crepusculus? I’ll come back to that later. Instead, let me begin with a fun fact: the shortest day of the year is December 21st. <\/p>\n
<\/p>\n
What? You already knew that fun fact? Okay smarty-pants, here’s a trick follow-up question: when is the earliest sunset? It turns out it comes well before the shortest day. <\/p>\n
Sounds crazy doesn’t it? But I’ll prove it to you, courtesy of a lovely File Exchange submission called SUNRISE<\/a> by François Beauducel<\/a>. As François says, it “computes sunrise and sunset times from any geographical location on Earth.” He works at the Institut de Physique du Globe de Paris<\/a>, so I’m pretty sure he knows what he’s doing here. <\/p>\n Watch this. Here I am in Natick, Massachusetts.<\/p>\n I’m going to use our lovely DATETIME<\/a> to create a year’s worth of dates. I’m going to start my calendar in July because I want December to be in the middle. <\/p>\n See what I did there? I created a vector of dates from July 1, 2019 to July 365, 2019! DATETIME handles the part that says “Silly! There is no July 200, 2019, so I’ll convert it to January 16, 2020 for you.” It’s a very convenient way to make date vectors. <\/p>\n The SUNRISE function expects an old-school DATENUM<\/a>, but that’s easily managed. <\/p>\n Merci François! Now we just calculate the hour and plot. Voilà!<\/p>\n You can see the two curves are somewhat offset. The earliest sunset advances, but the latest sunrise retreats, leaving the shortest day safely in place between them. <\/p>\n Where I live, the earliest sunset is on December 10th this year, fully eleven days before the solstice.<\/p>\n I like to call this date Crepusculus in honor of the word “crepuscular”, which means “of or relating to twilight”. I give it a fancy name because 1. it sounds cool and 2. it’s worth celebrating! After this date, the sun will set later every day until June. That’s almost as good as the solstice itself, especially for those of us that don’t often see the sun rise. And since it happens before the solstice, it gives me a sorely needed head start on solstitial merriment. <\/p>\n But it’s not the same for everybody! The earliest sunset changes by latitude. How much? Let’s look! We’ll loop across a whole bunch of latitudes and make a plot. <\/p>\n The NaT<\/a> function (NaT stands for “Not a Time”) is a nice shorthand for allocating a vector of DATETIMEs. It functions like ONES or ZEROS, only for dates. <\/p>\n Use the DIFF<\/a> function to find the leading edge of each step change. <\/p>\n Since tables are nice, let’s make a TABLE<\/a> of the mapping from latitude to date. <\/p>\nlat = 42.3;\r\nlon = -71.3;\r\nalt = 0;\r\ntzone = -5;\r\n<\/pre>\n
d = datetime(2019,7,1:365);\r\n<\/pre>\n
[srise,sset] = sunrise(lat,lon,alt,tzone,datenum(d));\r\n<\/pre>\n
sriseHour = 24*(srise - floor(srise));\r\nssetHour = 24*(sset - floor(sset));\r\nplot(d, sriseHour, d, ssetHour, 'LineWidth'<\/span>, 3)\r\nset(gca, ...<\/span>\r\n 'YLim'<\/span>,[0 24], ...<\/span>\r\n 'YTick'<\/span>,[0 6 12 18 24], ...<\/span>\r\n 'YTickLabel'<\/span>,{'Midnight'<\/span>,'6:00 AM'<\/span>,'Noon'<\/span>,'6:00 PM'<\/span>,'Midnight'<\/span>})\r\ntitle('Sunrise and Sunset'<\/span>)\r\nxlabel('Date'<\/span>)\r\nylabel('Time of Day'<\/span>)\r\nlegend({'Sunrise'<\/span>,'Sunset'<\/span>},'Location'<\/span>,'best'<\/span>)\r\nset(gca,'YDir'<\/span>,'reverse'<\/span>)\r\ngrid on<\/span>\r\nbox on<\/span>\r\n<\/pre>\n![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_02.png\")
<\/p>\n[~,ix] = min(ssetHour);\r\nd.Format = 'dd MMM yyyy HH:mm:ss'<\/span>;\r\nearliestSunset = datetime(sset(ix),'ConvertFrom'<\/span>,'datenum'<\/span>);\r\nlegend('off'<\/span>)\r\nline(d(ix),ssetHour(ix), ...<\/span>\r\n 'Marker'<\/span>,'o'<\/span>,'Color'<\/span>,[0.85 0.325 0.098], ...<\/span>\r\n 'LineWidth'<\/span>,3,'MarkerSize'<\/span>,18)\r\nxlim([datetime(2019,7,1) datetime(2020,7,1)])\r\nylim([12.0 24.0])\r\ntext(d(ix),ssetHour(ix), ...<\/span>\r\n sprintf('Crepusculus!\\n%s\\n\\n\\n\\n'<\/span>,string(earliestSunset)), ...<\/span>\r\n 'HorizontalAlignment'<\/span>,'center'<\/span>)\r\n<\/pre>\n![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_03.png\")
<\/p>\nlatitudeList = 1:0.2:65;\r\nearliestSunset = NaT(size(latitudeList));\r\n<\/pre>\n
for<\/span> i = 1:length(latitudeList)\r\n lat = latitudeList(i);\r\n \r\n [srise,sset] = sunrise(lat,lon,alt,tzone,datenum(d));\r\n ssetHour = 24*(sset - floor(sset));\r\n \r\n [~,ix] = min(ssetHour);\r\n earliestSunset(i) = d(ix);\r\nend<\/span>\r\nplot(latitudeList,earliestSunset,'.'<\/span>)\r\ngrid\r\nxlabel('Latitude (deg)'<\/span>)\r\nylabel('Date of the Earliest Sunset'<\/span>)\r\n<\/pre>\n![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_04.png\")
<\/p>\nxlim([25 50])\r\nix = find(diff(datenum(earliestSunset)))+1;\r\nhold on<\/span>\r\nplot(latitudeList(ix),earliestSunset(ix),'ro'<\/span>)\r\nhold off<\/span>\r\n<\/pre>\n![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_05.png\")
<\/p>\nlatTable = table(latitudeList(ix)',earliestSunset(ix)', ...<\/span>\r\n 'VariableNames'<\/span>,{'Latitude'<\/span>,'EarliestSunset'<\/span>})\r\n<\/pre>\nlatTable=45×2 table<\/em>\r\n Latitude EarliestSunset\r\n ________ ______________\r\n\r\n 1.6 05-Nov-2019 \r\n 2.2 06-Nov-2019 \r\n 2.8 07-Nov-2019 \r\n 3.4 08-Nov-2019 \r\n 4.2 09-Nov-2019 \r\n 4.8 10-Nov-2019 \r\n 5.6 11-Nov-2019 \r\n 6.2 12-Nov-2019 \r\n 7 13-Nov-2019 \r\n 7.8 14-Nov-2019 \r\n 8.6 15-Nov-2019 \r\n 9.4 16-Nov-2019 \r\n 10.2 17-Nov-2019 \r\n 11.2 18-Nov-2019 \r\n 12 19-Nov-2019 \r\n 13 20-Nov-2019 \r\n ⋮\r\n\r\n<\/pre>\n
![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_06.png\")
<\/p>\n![\"\"](\"https:\/\/blogs.mathworks.com\/community\/files\/crepusculus_07.png\")
<\/p>\n