complex(kind(1.d0)) function zgamr (z)

        implicit none
  c july 1977 edition.  w. fullerton, c3, los alamos scientific lab.
  c this version is an inaccurate preliminary one.  eventually this
  c routine should be a fundamental routine with no dependence on cgamma.
  c

        complex(kind(1.d0)) z, zlngam
        real(kind(1.d0)) x

        integer irold,ir
        external zlngam
  c
        zgamr = (0.0, 0.0)
        x = real (z)
        if (x.le.0.0 .and. aint(x).eq.x .and. 0.0.eq.aimag(z)) return
  c
        call entsrc (irold, 1)
        zgamr = zlngam(z)
        call erroff
        call entsrc (ir, irold)
        zgamr = exp (-zgamr)
  c
        return
        end