Determining when the sun shines on different sides of a building using astropy

Question

I recently installed motorized solar blinds on various windows in our house to reduce solar gain during the summer months. I have integrated the blinds with my hubitat home automation system, and am writing a python script to control the blinds. I've got lat/lon coordinates for my house, and I've determined that my house faces 145 degrees SE. I'm trying to figure out how to use astropy to determine which side(s) of the house the sun is shining on (if any) at the current time. My house shape is basically square, so the windows are basically either parallel or perpendicular to the direction the house faces.

I've found this example for getting the solar zenith. I'm assuming I need to use something like this example, plus the azimuth, adjusted somehow for the direction my house faces, to figure out which windows are getting sun at a given time, but I'm still working out the math. Thanks in advance for your help.

Answer 1

The code below will give you the alt-az coordinates for any specified location for any specified time. From there, I think the easiest method to determine which side of the house it will hit is to subtract 145deg from the AZ, so that 0 means it's hitting the front dead on. I'll call this the "front angle".

If the front angle is between 270 and 90 (passing through 0), it's hitting the front. Between 0 and 180 the left, between 90 and 270 the back, and between 180 and 0 the right.

You'll also want to validate the sun's altitude is > 0

#!/usr/bin/python
import astropy.units as u
from astropy.time import Time
from astropy.coordinates import EarthLocation, AltAz
from astropy.coordinates import get_sun

louisville = EarthLocation.from_geodetic(-85.7636,38.2464000,0)
utcoffset = -4*u.hour  # Eastern Daylight Time
time = Time('2012-7-12 20:00:00') - utcoffset
aaframe=AltAz(obstime=time,location=louisville)

sun=get_sun(time)

sunaa=sun.transform_to(aaframe)
print(f"Sun Alt, Az = {sunaa.alt} {sunaa.az}")

frontAngle=sunaa.az-145*u.deg

if(frontAngle<0): frontAngle=frontAngle+360*u.deg
print(f"Front Angle = {frontAngle}")

if (frontAngle.is_within_bounds('270d', '360d') or frontAngle.is_within_bounds('0d', '90d')): print("front")
if (frontAngle.is_within_bounds('0d', '180d')): print("left")
if (frontAngle.is_within_bounds('90d', '270d')): print("back")
if (frontAngle.is_within_bounds('180d', '360d')): print("right")