subroutine epem3j_z(p,z)
************************************************************************
*     John M. Campbell                                                 *
*     November, 2008.                                                  *
************************************************************************
      implicit none
      include 'constants.f'
      include 'qcdcouple.f'
      include 'scale.f'
      include 'PR_new.f'
      include 'agq.f'
      include 'colstruc.f'
      integer is
      double precision z,xl34,xl35,xl45,p(mxpart,4),dot
      double precision ff_qq,ff_gg

      xl34=dlog(+two*dot(p,3,4)/musq)
      xl35=dlog(+two*dot(p,3,5)/musq)
      xl45=dlog(+two*dot(p,4,5)/musq)

      do is=1,3
c--- leading colour (Ca) goes in (0,0)
c--- subleading colour (2*Cf-Ca) goes in (0,1)
c--- Tr terms go in (1,0)
      caonly=.true.
      nfonly=.false.
      Q1(g,g,g,is)=ason4pi*ca*(ff_qq(z,xl45,is)+0.5d0*ff_gg(z,xl45,is)
     &                        +ff_qq(z,xl35,is)+0.5d0*ff_gg(z,xl35,is))

      Q1(g,g,q,is)=ason4pi*ca*(-ff_qq(z,xl34,is)-ff_qq(z,xl34,is))/xnsq

      caonly=.false.
      nfonly=.true.
      Q1(q,q,g,is)=ason4pi*ca*(0.5d0*ff_gg(z,xl45,is)
     &                        +0.5d0*ff_gg(z,xl35,is))
      enddo
      
      return
      end