program xdeg2dms
!
!     Unit test for deg2dms.for
!
!     B. Knapp, 94.12.21
!
C
C     RCS DATA
C     
C     $Header$
C     
C     $Log$
C
C
      implicit none
      real*8 deg,deg2,sec
      integer*4 ideg,min

    1 write(*,'(/A$)') '  Degrees: '
      read(*,*) deg
      call deg2dms(deg,ideg,min,sec)
      write(*,10) deg,ideg,min,sec
   10 format('  ',f12.6,' = ',i4,i4,f16.12)
      call dms2deg(ideg,min,sec,deg2)
      write(*,20) deg, deg2, deg-deg2
   20 format(2f20.12,1pe12.4)
      goto 1
      end