diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/solar.py | 40 |
1 files changed, 34 insertions, 6 deletions
diff --git a/src/solar.py b/src/solar.py index a0620bc..c8d7677 100644 --- a/src/solar.py +++ b/src/solar.py @@ -18,6 +18,31 @@ import math import time +def sun(latitude, longitude, t = None, low = -6.0, high = 3.0): + ''' + Get the visibility of the Sun + + @param latitude:float The latitude component of your GPS coordinate + @param longitude:float The longitude component of your GPS coordinate + @param t:float? The time in Julian Centuries, `None` for current time + @param low:float The 100 % night limit elevation of the Sun (highest when not visible) + @param high:float The 100 % day limit elevation of the Sun (lowest while fully visible) + @return :float The visibilty of the sun, 0 during the night, 1 during the day, + between 0 and 1 during twilight. Other values will not occur. + ''' + t = julian_centuries() if t is None else t + e = solar_elevation(latitude, longitude, t) + e = (e - low) / (high - low) + return min(max(0, e), 1) + + + +# The following functions are used to calculate the result for `sun` +# (most of them) but could be used for anything else. There name is +# should tell you enough, `t` (and `noon`) is in Julian centuries +# except for in the convertion methods + + def julian_day_to_epoch(t): return (jd - 2440587.5) * 86400.0 @@ -37,12 +62,21 @@ def julian_centuries_to_epoch(t): return julian_day_to_epoch(julian_centuries_to_julian_day(t)) def epoch(): + ''' + Get current POSIX time + ''' return time.time() def julian_day(): + ''' + Get current Julian Day time + ''' return epoch_to_julian_day(epoch()) def julian_centuries(): + ''' + Get current Julian Centuries time + ''' return epoch_to_julian_centuries(epoch()) def radians(deg): @@ -151,9 +185,3 @@ def solar_elevation(latitude, longitude, t = None): rc = solar_elevation_from_time(rc, latitude, longitude) return degrees(rc) -def sun(latitude, longitude, t = None, low = -6.0, high = 3.0): - t = julian_centuries() if t is None else t - e = solar_elevation(latitude, longitude, t) - e = (e - low) / (high - low) - return min(max(0, e), 1) - |