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