      subroutine dkqqb_hww_gs(p,msqc)
      implicit none
************************************************************************
*     Author: R.K. Ellis                                               *
*     May, 2012.                                                       *
*     calculate the subtraction term for radiation in the              *
*     W^- decay for the process                                        *
*                                                                      *
*     q(-p1)+qbar(-p2) = H-> W^+(nu(p3)+e+(p4))                        *
*                       +W^-(e-(p5)+nubar(p6)+g(p7))                   *
*                                                                      *
************************************************************************
      include 'constants.f'
      include 'ptilde.f'
      include 'qqgg.f'
      include 'incldip.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
      external qqb_hww,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_hww,donothing_gvec)
      call dips(2,p,6,7,5,sub67_5,dsubv,msq67_5,dummyv,
     & qqb_hww,donothing_gvec)

      msqc(1,0,0)=sub57_6(qq)*msq57_6(0,0)*2d0*cf
      msqc(2,0,0)=sub67_5(qq)*msq67_5(0,0)*2d0*cf


      return
      end

