city=;yearrange=Range[2003,2020];dailyresults={};dailyresults=Reap[Do[dates=DateRange[{year,4,16},{year,5,12}];Do[grahaRAs={};fixeddate=DateObject[date,TimeZone0];T=DateDifference[DateObject[{2000,1,1,11,58.816},TimeZone"UTC"],fixeddate]/Quantity[36525,"Days"];delta=UnitConvert[Quantity[5028.796195*T+1.1054348*T^2,"ArcSeconds"],"Degrees"];ayanamsha=QuantityMagnitude[Quantity[23.439291111,"Degrees"]+delta];grahaRAs={Mod[QuantityMagnitude[UnitConvert[StarData[Entity["Star","Sun"],EntityProperty["Star","RightAscension",{"Date"fixeddate,"Location"city}]],"Degrees"]]-ayanamsha,360],Mod[QuantityMagnitude[UnitConvert[PlanetData[Entity["Planet","Jupiter"],EntityProperty["Planet","RightAscension",{"Date"fixeddate,"Location"city}]],"Degrees"]]-ayanamsha,360],Mod[QuantityMagnitude[UnitConvert[PlanetData[Entity["Planet","Mars"],EntityProperty["Planet","RightAscension",{"Date"fixeddate,"Location"city}]],"Degrees"]]-ayanamsha,360]};Sow[{grahaRAs,Mean[Table[timeseries[[loc]][DateObject[date]],{loc,1,Length[timeseries]}]]}],{date,dates}];,{year,yearrange}]];