complex(kind(1.d0)) function zatan2 (csn, ccs)

        implicit none
  c april 1977 version.  w. fullerton, c3, los alamos scientific lab.

        complex(kind(1.d0)) csn, ccs, zatan
        real(kind(1.d0)) pi

        external zatan
        data pi / 3.1415926535 8979323846d0 /
  c
        if (abs(ccs).eq.0.) go to 10
  c
        zatan2 = zatan (csn/ccs)
        if (real(ccs).lt.0.) zatan2 = zatan2 + pi
        if (real(zatan2).gt.pi) zatan2 = zatan2 - 2.0*pi
        return
  c
   10   if (abs(csn).eq.0.) call seteru (
       1  34hzatan2  called with both args zero, 34, 1, 2)
  c

        zatan2 = cmplx (sign(0.5*pi,real(csn)), 0.0, kind(1.d0))

  c
        return
        end