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

      subroutine getrops(pressure,temperature,nspecies,
     1                   ysp,data_out)

      IMPLICIT NONE
      
      include 'parameter.par' 
      include 'chem.com'
      include 'chem_ref.com'

      double precision pressure, press, temperature
      double precision temp, rop
      double precision sum,V,rho_k,rho

      integer k, nspecies

      real*8 cspl(nspecies)
      real*8 ysp(nspecies)  !Y(species)
      real*8 data_out(numreac)  !Return the reactions rops


c     press is dimensional pressure [N/m^2] and temp is non-dimensional temperature
#ifdef NON_DIM
      temp = temperature
      press = pressure * pref 
#else
      temp = temperature / tempref
      press = pressure
#endif

c-------------------------------------------------------------------------------
c	Evaluate the cell density of the mixture:

c       1. rho(mix) = m(mix) / V(mix)
c	2. m(mix) = 1 kgr
c       3. V(mix) = V(1) + V(2) +...+ V(ns), for species 1..ns
c       4. V(k) = m(k) / rho(k)
c       5. m(k) = ys(k,i,j) * m(mix) = ys(k,i,j)
c       6. rho(k) = pressure / ( rg(k)*temp(i,j)*tempref )
c       7. rg(k) = r0 / wt

      V = 0.0
      sum = 0.0

      do k = 1,nspecies
	      
         rg(k) = r0 / (1000.0 * wt(k))
         rho_k= press / (1000.0 * rg(k) * temp * tempref ) 
         V = ysp(k) / rho_k

         sum = sum + V
      end do

c     non-dimensionalize rho
      rho = (1.d0 / sum) / rhoref 

c-------------------------------------------------------------------------------
c	find the cspl's

      call y2c(rho,ysp,cspl)  

c-------------------------------------------------------------------------------
c ... evaluate total reaction rates
      do k=1,numreac
         data_out(k) = rop(k,temp,tempref,cspl) 
      end do
c-------------------------------------------------------------------------------
      return
      end

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