summaryrefslogtreecommitdiffstats
path: root/src
diff options
context:
space:
mode:
authorMattias Andrée <maandree@operamail.com>2014-02-17 01:41:14 +0100
committerMattias Andrée <maandree@operamail.com>2014-02-17 01:41:14 +0100
commitc70be6e3a3f02276f06a6da869a9ff798d393e2a (patch)
treec7a71a3272382e6a03e8b746767ea2b3814438bc /src
parentdoc (diff)
downloadblueshift-c70be6e3a3f02276f06a6da869a9ff798d393e2a.tar.gz
blueshift-c70be6e3a3f02276f06a6da869a9ff798d393e2a.tar.bz2
blueshift-c70be6e3a3f02276f06a6da869a9ff798d393e2a.tar.xz
doc
Signed-off-by: Mattias Andrée <maandree@operamail.com>
Diffstat (limited to '')
-rw-r--r--src/solar.py40
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)
-