Hello Jimeno,I'm not sure exactly what you're looking for, but you should be able to use the solar functions, perhaps in conjunction with other packages, to answer your question. I see that as you mention, when you run Points Solar Radiation on a single day, you get back the time steps but not the actual time that it corresponds to -- one way to get back the corresponding actual hours would be to use a script such as this one calculating sunrise and sunset or use a package like PyEphem to get more precise calculations.From there, you should be able to combine the solar observations from ArcGIS with 'actual time', so you can correspond the output T0 with the sunrise time and the final Tn with the sunset time. So my psuedo-code for this on a single day would be something like:
import arcpy
import ephem
julian_day = 165
# run solar calculation for all sites on this day; defaults used here
arcpy.PointsSolarRadiation("elevation", "my_sites",
"output_global_radiation", "", 35, 200,
arcpy.sa.TimeWithinDay(julian_day, 0, 24))
# figure out sunrise and sunset. if the points are close together, do this
# once for all sites, if they're geographically spread over a large area, you'll need
# to make these calculations for each site.
obs = ephem.Observer()
obs.lat = ' 34' # obtained from my_sites
obs.long = '-120.0' # obtained from my_sites
obs.elev = 15 # extract from elevation raster at site location
sunrise = obs.previous_rising(ephem.Sun()) # sunrise
noon = obs.next_transit(ephem.Sun(), start=sunrise) # solar noon
sunset = obs.next_setting(ephem.Sun()) # sunset
# then, map those times to the T0...Tn time steps, which are by default in half-hour intervals.
# output this result to a CSV or the like.
Does that help get you started? I think something like this is easier than trying to dig up the internals to figure out the corresponding times.cheers,Shaun
