coordinates.f90: new subroutine horiz2eq() to convert coordinates from horizontal to equatorial