c
c This subroutine calculates total internal coordinate gradient.
c Depending on the number of function evaluations, either whole energy
c
c This subroutine calculates total internal coordinate gradient.
c Depending on the number of function evaluations, either whole energy