subroutine dkqqb_ww_gs(p,msqc)
      implicit none
************************************************************************
*     Author: R.K. Ellis                                               *
*     June, 2012.                                                      *
*     calculate the subtraction term for radiation in the              *
*     W^- decay for the process                                        *
*                                                                      *
*     q(-p1)+qbar(-p2) = W^+(nu(p3)+e+(p4))                            *
*                       +W^-(s(p5)+cbar(p6)+g(p7))                     *
*                                                                      *
************************************************************************
      include 'constants.f'
      include 'ptilde.f'
      include 'qqgg.f'
      double precision msqc(maxd,-nf:nf,-nf:nf),p(mxpart,4)
      double precision 
     & msq57_6(-nf:nf,-nf:nf),msq67_5(-nf:nf,-nf:nf),
     & dummyv(-nf:nf,-nf:nf),sub57_6(4),sub67_5(4),dsubv
      integer j,k,nd
      logical incldip(0:maxd)
      common/incldip/incldip
      external qqb_ww,donothing_gvec

      ndmax=2
      
      do j=-nf,nf
      do k=-nf,nf
      do nd=1,ndmax
        msqc(nd,j,k)=0d0
        incldip(nd)=.true.
      enddo
      enddo
      enddo

      call dips(1,p,5,7,6,sub57_6,dsubv,msq57_6,dummyv,
     & qqb_ww,donothing_gvec)
      call dips(2,p,6,7,5,sub67_5,dsubv,msq67_5,dummyv,
     & qqb_ww,donothing_gvec)

      do j=-nf,nf
      do k=-nf,nf
      msqc(1,j,k)=sub57_6(qq)*msq57_6(j,k)*2d0*cf
      msqc(2,j,k)=sub67_5(qq)*msq67_5(j,k)*2d0*cf
      enddo
      enddo      
      return
      end