@@ -28,7 +28,7 @@ def solar_to_UTC(time_solar: datetime, lon: float) -> datetime:
28
28
"""
29
29
return time_solar - timedelta (hours = (np .radians (lon ) / np .pi * 12 ))
30
30
31
- def UTC_offset_hours (geometry : rt .RasterGeometry ) -> rt .Raster :
31
+ def UTC_offset_hours_for_area (geometry : rt .RasterGeometry ) -> rt .Raster :
32
32
"""
33
33
Calculates the UTC offset in hours for a given raster geometry.
34
34
@@ -40,7 +40,7 @@ def UTC_offset_hours(geometry: rt.RasterGeometry) -> rt.Raster:
40
40
"""
41
41
return rt .Raster (np .radians (geometry .lon ) / np .pi * 12 , geometry = geometry )
42
42
43
- def day_of_year (time_UTC : datetime , geometry : rt .RasterGeometry ) -> rt .Raster :
43
+ def solar_day_of_year_for_area (time_UTC : datetime , geometry : rt .RasterGeometry ) -> rt .Raster :
44
44
"""
45
45
Calculates the day of the year for a given UTC time and raster geometry.
46
46
@@ -53,15 +53,15 @@ def day_of_year(time_UTC: datetime, geometry: rt.RasterGeometry) -> rt.Raster:
53
53
"""
54
54
doy_UTC = time_UTC .timetuple ().tm_yday
55
55
hour_UTC = time_UTC .hour + time_UTC .minute / 60 + time_UTC .second / 3600
56
- UTC_offset_hours = UTC_offset_hours (geometry = geometry )
56
+ UTC_offset_hours = UTC_offset_hours_for_area (geometry = geometry )
57
57
hour_of_day = hour_UTC + UTC_offset_hours
58
58
doy = doy_UTC
59
59
doy = rt .where (hour_of_day < 0 , doy - 1 , doy )
60
60
doy = rt .where (hour_of_day > 24 , doy + 1 , doy )
61
61
62
62
return doy
63
63
64
- def hour_of_day (time_UTC : datetime , geometry : rt .RasterGeometry ) -> rt .Raster :
64
+ def solar_hour_of_day_for_area (time_UTC : datetime , geometry : rt .RasterGeometry ) -> rt .Raster :
65
65
"""
66
66
Calculates the hour of the day for a given UTC time and raster geometry.
67
67
@@ -73,7 +73,7 @@ def hour_of_day(time_UTC: datetime, geometry: rt.RasterGeometry) -> rt.Raster:
73
73
rt.Raster: The hour of the day.
74
74
"""
75
75
hour_UTC = time_UTC .hour + time_UTC .minute / 60 + time_UTC .second / 3600
76
- UTC_offset_hours = UTC_offset_hours (geometry = geometry )
76
+ UTC_offset_hours = UTC_offset_hours_for_area (geometry = geometry )
77
77
hour_of_day = hour_UTC + UTC_offset_hours
78
78
hour_of_day = rt .where (hour_of_day < 0 , hour_of_day + 24 , hour_of_day )
79
79
hour_of_day = rt .where (hour_of_day > 24 , hour_of_day - 24 , hour_of_day )
0 commit comments