LunarPhase function

Uses library(oce) function moonAngle to assign lunar fraction and phase
for a given date and position. Use ?moonAngle for help