!***********************************************************************************************************************************
!  FRAC
!
!  Return the fractional part of a double-precision number.
!***********************************************************************************************************************************

      FUNCTION FRAC (X) RESULT (Y)

      IMPLICIT NONE

      DOUBLE PRECISION, INTENT(IN) :: X
      DOUBLE PRECISION :: Y, Z

      Z = ABS(X)
      Y = Z - INT(Z)
      Y = SIGN(Y,X)

      RETURN

      END FUNCTION FRAC