New subroutine planet_visibility_tonight() computes when a planet is visible at night

* New subroutine planet_visibility_tonight() computes when a planet is visible at night:
  * input: JD, planet ID, max Sun altitude, min planet altitude
  * output: twilight times, rise/set times for the Sun and the planet, transit altitudes
    	    begin and end times for the planet visibility