  function S = StartTimes()
% S is an 80-by-1 vector that houses season start times.
% To specify its entries, let
%    Spring(k) = the time inspring starts in year k
%    Summer(k) = the time summer starts in year k
%    Autumn(k) = the time autumn starts in year k
%    Winter(k) = the time winter starts in year k
%
% Time is measured in days from the 2000 vernal equinox, 
% i.e., Spring(2000) = 0;
% 
%  S = [Summer(2000) ; Autumn(2000) ; Winter(2000)  Spring(2001) ;...
%       Summer(2001) ; Autumn(2001) ; Winter(2001)  Spring(2002) ;...
%       Summer(2002) ; Autumn(2002) ; Winter(2002)  Spring(2003) ;...
%                              etc
%       Summer(2018) ; Autumn(2018) ; Winter(2018)  Spring(2019) ;...
%       Summer(2019) ; Autumn(2019) ; Winter(2019)  Spring(2020) ]
%  
% Terminology: "Predict the first k season events" means  
%              "Predict the values in s(1:k)"
%              
SeasonLengths = [...
  % 2000-2004 
  92.75902778   93.65277778   89.83958333    88.99583333;...
  92.75486111   93.64305556   89.84513889    88.99652778;...
  92.75555556   93.64652778   89.84652778    88.99027778;...
  92.75694444   93.65069444   89.84513889    88.98958333;...
  92.75555556   93.64791667   89.84166667    88.99375000;...
  % 2005-2009
  92.75902778   93.65069444   89.84166667    88.99375000;...
  92.75000000   93.65069444   89.84652778    88.98958333;...
  92.74930556   93.65625000   89.84513889    88.98611111;...
  92.75763889   93.65625000   89.84722222    88.98611111;...
  92.75138889   93.64791667   89.85277778    88.98958333;...
  % 2010-2014
  92.74722222   93.65347222   89.85347222    88.98819444;...
  92.74652778   93.65902778   89.85069444    88.98888889;...
  92.74652778   93.65277778   89.84930556    88.99305556;...
  92.75138889   93.65277778   89.85208333    88.99027778;...
  92.74583333   93.65138889   89.85694444    88.98750000;...
  % 2015-2019
  92.74513889   93.65486111   89.85208333    88.98750000;...
  92.75277778   93.65763889   89.84930556    88.98958333;...
  92.74652778   93.65138889   89.85138889    88.99097222;...
  92.74444444   93.65763889   89.85347222    88.98263889;...
  92.74722222   93.66388889   89.85347222    88.97986111;...
  ];

S = cumsum(reshape(SeasonLengths',80,1));