def sunazel_lacc(JD, Lat=50.867387222, Lon=0.33612916666, H=75.357): '''Sun's azimuth and elevation at given JD and geodetic location. Low-precision algorithm''' x, y, z = sunpos_lacc(JD) az, el, r = cds.geo2azel(x, y, z, Lat, Lon, H) return az, el, r
def moonazel(JD=2448908.5, Lat=50.867387222, Lon=0.33612916666, H=75.357): x, y, z = moonpos(JD) az, el, r = cds.geo2azel(x, y, z, Lat, Lon, H) return az, el, r