moon_phase.py

import calendar
import ephem


YEAR = 2018


observer = ephem.Observer()
observer.lat = -3222.78
observer.lon = 2048.60
observer.elevation = 1822


moon = ephem.Moon()


for month in range(1, 13):
    weekday, days_in_month = calendar.monthrange(YEAR, month)

    for day in range(1, days_in_month+1):
        date = '{}/{}/{} 00:00:00'.format(YEAR, month, day)
        observer.date = date
        moon.compute(observer)
        print(date, int(round(moon.moon_phase * 100)))