This subroutine should return double precision.
This subroutine should return double precision.