function is=seasons(jd);
%function is=seasons(jd);
% returns the indices of the Mass Bay seasons
season(1,:)=[1987  2 9 6  0 0];
season(2,:)=[1987  3 1 0  0 0];
season(3,:)=[1987  5 11 0  0 0];

for i=1:3,
  is(i)=tind(jd,season(i,:));
end
