#intervals += [TimeInterval(datetime(years,month,1), datetime(years,month,calendar.monthrange(years,month)[1])+day-wedge) for month in range(1,3) ]
#intervals += [TimeInterval(datetime(years,month,1), datetime(years,month,calendar.monthrange(years,month)[1])+day-wedge) for month in range(7,9) ]
#intervals = [TimeInterval(datetime.strptime('{}{:03d}'.format(years, jday),'%Y%j'), datetime.strptime('{}{:03d}'.format(years, jday),'%Y%j')+day-wedge) for jday in range(1,366)]
#intervals = [TimeInterval(datetime.strptime('{}{:03d}'.format(years, jday),'%Y%j'), datetime.strptime('{}{:03d}'.format(years, jday),'%Y%j')+day-wedge) for jday in jdays]