lunar_phase function

lunar_phase