moon_routines: added subroutine moon_phase_next()

* moon_routines: added subroutine moon_phase_next():
  * compute the next lunar phase after a given JD