      PROGRAM XZENITH_ANGLE
C
C     EXERCISE DRIVER FOR ZENITH_ANGLE.FOR
C
C     B. KNAPP, 95.03.06
C
C
C     RCS DATA
C     
C     $Header$
C     
C     $Log$
C
C
      IMPLICIT NONE
C
      REAL*8 JD,LAT,LON,Z
      INTEGER*4 P, STATUS
C
C     EXTERNALS
      EXTERNAL ZENITH_ANGLE
C
C
    1 WRITE(*,'(/A$)') ' JD, P, LAT, LON? '
      READ(*,*) JD, P, LAT, LON
C
      CALL ZENITH_ANGLE( JD, LAT, LON, P, Z, STATUS )
      WRITE(*,10) Z
   10 FORMAT(' ZENITH ANGLE (DEG): ',F24.12)
      GOTO 1
      END
