c*******************************************************************************
c ... fillup array of thermal conductivity htcond(j)
c ... dimensional thermal conductivity at (j)
c ... if in_n2 is .true. then use a pure N2 gas
c ... otherwise, use the local 1D flame mixture fractions at each point
c ... assumes species(nsp) = N2

	subroutine gthtcond

        include 'parameter.par'
#include "precision.com"
	include 'common.com'

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

	do j=1,n
	   vtemp = temp(j)
	   if(in_n2)then
	      do k=1,nsp-1
		 xsp(k) = 0.0
	      end do
	      xsp(nsp)=1.0
	   else
	      do k=1,nsp
	         xsp(k) = xs(j,k)
	      end do
	   endif
	   call getmxcon(vtemp+epso2,nsp,xsp,htc_p)
	   call getmxcon(vtemp-epso2,nsp,xsp,htc_m)
	   htcond(j) = ( htc_p - htc_m ) / epstemp
	end do

	return
	end

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