[pycalcal] push by enrico.spinielli - still converting...wrongly on 2010-02-02 14:19 GMT

[pyca...@googlecode.com]

Revision: 18113e0061
Author: Enrico Spinielli <enrico.s...@gmail.com>
Date: Tue Feb 2 06:19:37 2010
Log: still converting...wrongly
http://code.google.com/p/pycalcal/source/detail?r=18113e0061

Modified:
/pycalcal.nw

=======================================
--- /pycalcal.nw Sun Jan 31 05:16:34 2010
+++ /pycalcal.nw Tue Feb 2 06:19:37 2010
@@ -3361,15 +3361,31 @@
'latitude' is the observer's geographic latitude.
NOTE: azimuth is measured westward from the South.
"""
- ca = cos_degrees(azimuth)
- sl = sin_degrees(latitude)
- cl = cos_degrees(latitude)
- H = degrees_from_radians(
- atan2(sin_degrees(azimuth),
- ca * sl + tan_degrees(altitude) * cl))
- d = arcsin_degrees(sl * sin_degrees(altitude) -
- cl * cos_degrees(altitude) * ca)
- return [H, d]
+ # ca = cos_degrees(azimuth)
+ # sl = sin_degrees(latitude)
+ # cl = cos_degrees(latitude)
+ # H = degrees_from_radians(
+ # atan2(sin_degrees(azimuth),
+ # ca * sl + tan_degrees(altitude) * cl))
+ # d = arcsin_degrees(sl * sin_degrees(altitude) -
+ # cl * cos_degrees(altitude) * ca)
+ a=altitude #altitude
+ A=azimuth #azimuth
+ phi=latitude #latitude
+
+ #cos delta * sin H = cos a * sin A
+ # sin delta = sin phi * sin a + cos phi * cos a * cos A
+ #cos delta * cos H = cos phi * sin a + sin phi * cos a * cos A
+
+ delta = arcsin_degrees(sin_degrees(phi) * sin_degrees(a) +
+ cos_degrees(phi) * cos_degrees(a) *
cos_degrees(A))
+ H = arcsin_degrees((cos_degrees(a) * sin_degrees(A)) /
cos_degrees(delta))
+
+
+ #alpha= #right ascension
+ #delta= #declination
+ #H= #local hour angle
+ return [H, delta]