c*******************************************************************************
c ... evaluate prod rate of sp ksp

	subroutine wchem(tc,tempref,cspl,wspl)

	IMPLICIT NONE

        include 'parameter.par' 
	include 'chem.com'

	real*8 wspl(nsmx),cspl(nsmx)

	integer ksp,k,ir,ip3
	
	double precision ropl,rop,tc,tempref



c-------------------------------------------------------------------------------
c ... init

        do ksp=1,ns
           wspl(ksp) = 0.0d0
        end do

c-------------------------------------------------------------------------------
c ... evaluate contributions to reactants and products reaction rates

        do k=1,numreac

           ropl = rop(k,tc,tempref,cspl)

           do ir = 1, nreac(k)
              ksp    = nunk(ir,k)
              wspl(ksp) = wspl(ksp) + anu(ir,k)*ropl         ! reactants
           end do

           do ip3 = 4, nprod(k)+3
              ksp    = nunk(ip3,k)
              wspl(ksp) = wspl(ksp) + anu(ip3,k)*ropl        ! products
           end do

        end do

c-------------------------------------------------------------------------------
c ... final scaling

        do ksp=1,ns
           wspl(ksp) = wspl(ksp)*spwt(ksp)
        end do

c-------------------------------------------------------------------------------

        return
	end

c*******************************************************************************