[Tutor] Angles
Carlos
carloslara at web.de
Tue Nov 28 22:59:06 CET 2006
Hello to All:
I'm trying to write a sun positioning system with the help of python.
The idea is that you give the program a location, date and hour and it
returns the position of the sun.
I found a webpage that details the math behind this
(http://www.usc.edu/dept/architecture/mbs/tools/vrsolar/Help/solar_concepts.html)
, it was fairly trivial to translate this to python until I got to the
azimuth equation. It looks like this:
x_azm = sin(hour_R) * cos(decl_R)
y_azm = (-(cos(hour_R))*cos(decl_R)*sin(lat_R))+(cos(lat_R)* sin(decl_R))
azimuth = atan(x_azm/y_azm)*TODEGREE
where:
Alt = Altitude
Azm = Azimuth
Decl = Declination
HAngle = Hour angle
alt_R = Altitude in radians
azm_R = Azimuth in radians
lat_R = Latitude in radians
hour_R = Hour angle in radians
x_azm = x component of azimuth
y_azm = y component of azimuth
TODEGREE = Constant equal to 180/p
My python code for this particular equation:
from math import *
def Az(Lat, Dec, H_Ang):
lat_r = radians(Lat)
decl_r = radians(Dec)
hour_r = radians(H_Ang)
x_azm = sin(hour_r) * cos(decl_r)
y_azm = (-(cos(hour_r)) * cos(decl_r) * sin(lat_r)) +
(cos(lat_r) * sin(decl_r))
Azimuth = degrees(atan(x_azm/y_azm))
Azimuth = Azimuth *(-1)
return Azimuth
I was never very good at trigonometry, but looks like my translation of
the equation is ok and the problem is some kind of python behavior,
because whenever the results exceed 100° (deg) Python returns the
complementary angle, it is possible to avoid this? Or I'm overlooking
somethig?
You can check this by yourself:
If you use Az(34, -21.67, -150) you will get 72.79, the result in the
webpage tool will be 107.21. And 72.79 + 107.21 = 180
but if you use Az(34, -21.67, -150) you will get 54.88 in both, my
script and the webpage tool.
If you want to take a look at the webpage tool this is the link:
http://www.usc.edu/dept/architecture/mbs/tools/vrsolar/frameset.html
Thanks in advance for your help,
Carlos
More information about the Tutor
mailing list