solar.solar_depression = 'civil'
#l = Location()
#l.name = 'Wohnwagen'
#l.latitude = latitude
#l.longitude = longitude
#l.timezone = 'Europe/Berlin'
#l.elevation = 0
#sun = l.sun()
#data['sunrise'] = str(sun['sunrise'])
#data['sunset'] = str(sun['sunset'])
data['azimuth'] = solar.solar_azimuth(dateandtime, latitude, longitude)
data['elevation'] = solar.solar_elevation(dateandtime, latitude, longitude)
if data['elevation'] < 0:
data['elevation'] = 0
data['zenith'] = solar.solar_zenith(dateandtime, latitude, longitude)
# Set up a client for InfluxDB
dbclient = InfluxDBClient(config[influx_server]['host'], config[influx_server]['port'], config[influx_server]['username'], config[influx_server]['password'], config[influx_server]['database'])
json_body = [{
"measurement": "power",
"fields": {
"pv_w": float(data['/Pv/W']),
"pv_v": float(data['/Pv/V']),
"pv_a": float(data['/Pv/I']),
"pv_state": float(data['/State']),
"elev": float(data['elevation']),
"bat_a": float(resp[battery]['/Dc/0/Current']),
"bat_w": float(resp[battery]['/Dc/0/Power']),
"bat_v": float(resp[battery]['/Dc/0/Voltage']),
def test_Astral_SolarZenith():
dd = Astral()
dt = datetime.datetime(2015, 2, 3, 9, 0, 0, tzinfo=pytz.UTC)
zenith = dd.solar_zenith(dt, 51.5, -0.12)
assert float_almost_equal(zenith, 90.0 - 9.97, 0.1)
azRad, elevRad = (360. - az + 90.)*np.pi/180., (90.-elev)*np.pi/180.
Sx, Sy = calcFiniteSlopes(elevGrid, dx) # Calculate slope in X and Y directions
AspectRad = np.arctan2(Sy, Sx) # Angle of aspect
SmagRad = np.arctan(np.sqrt(Sx**2. + Sy**2.)) # magnitude of slope in radians
return (((np.cos(elevRad) * np.cos(SmagRad)) + (np.sin(elevRad)* np.sin(SmagRad) * np.cos(azRad - AspectRad))))*255
a = Astral()
Latitude = 47.074531
longitude = 12.846210
b = datetime.datetime(2015,10,3,7,30)
c= a.solar_azimuth(b,Latitude,longitude)
d= a.solar_elevation(b,Latitude,longitude)
print d, a.solar_zenith(b, Latitude, longitude)
print c
headinfo,dtm = ascii.read_ascii("C:\Master\settings/vernagtferner14-16/dgm_vernagtferner.txt")
test = calcHillshade(dtm,headinfo[-2],c,d)
print np.min(test), np.max(test)
plt.imshow(test, cmap="Greys")
plt.show()
ascii.write_ascii("test.asc",headinfo,test,format="%i")
