MoonPhase

new MoonPhase()
Instance Members
flmoon(n, nph)
nextphase(year, mon, day, hr, min, nph, jd)
isDST(date)
Julian2Calendar(jd, mon, day, year)
Calendar2Julian(year, month, day, mon)

sun_anomaly

Sun's mean anomaly

sun_anomaly

moon_anomaly

Moon's mean anomaly

moon_anomaly

xtra

Not sure of what's up here, but two of the numbers very interesting Notice that 2415020.75933 is used a referce epochs for mean time of lunar cycle which is 1900 Jan 1 at 6:13 AM And 29.53058868 is the "Synodic month" or the mean time in days between new moons.

xtra

xtra

Looks like this is all being done to adjust the variations in the maen lunar cycle of 29.53058868 Which from what I understand is due the moons orbit in relationship to the earths orbit around the sun.

xtra

jd

Ok 2415020 is Julian date 1899 at noon + .75933 is the Julian date of a new moon 28 + 1.53058868 is mean time of a lunar cycle , and 7 close to 7.38264717 a lunar phase but unlike the first two the tractional part dose not seem to be accounted for.

jd

int_part

This parts easy just put the integer and fractional parts right

int_part

n

12.37 is the number of lunar cycles per year n then is the number of lunar cycles since 1900 and the year month entered

n