program deg2dms
      
      real degdec
      integer degree, minute, sec
      
      print *, 'Enter a degree in decimals'
      read *, degdec
      
      call deg2dmssub(degdec,degree,minute,sec)
      
      print *, degdec, degree, minute, sec
      
      stop
      end
C-------------------------------------------------------------
      subroutine deg2dmssub(degdec, degree, minute, sec)
      
      real degdec
      integer degree, minute, sec, afterdec
      
      degree = int(degdec)
      
      afterdec = (degdec-degree)*1000
      
      print *, afterdec
      
      minute = afterdec * 60 / 1000
      second = afterdec * 3600 / 1000
      
      
      
      return
      end