quick-ish calculation of sun and moon location

sunmoon.py

from astropy import coordinates as coord
from astropy.time import Time
from astropy import units as u

lapalma = coord.EarthLocation.of_site('lapalma')
times = Time.now() + u.minute*np.linspace(-720, 720, 1440)
altaz_frame = coord.AltAz(location=lapalma, obstime=times)
sun = coord.get_body('sun', location=lapalma, time=times)
moon = coord.get_body('moon', location=lapalma, time=times)

sun_altaz = sun.transform_to(altaz_frame)
moon_altaz = moon.transform_to(altaz_frame)