c*******************************************************************************
c ... capd(j,k) = bin_D_k_nsp	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
c-------------------------------------------------------------------------------

	subroutine gcapd

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

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

	do j=1,n

	   nspcon = (nsp-1)*nsp

	   vtemp = temp(j)		! dimensional temp
	   vpres = p0v    		! dimensional pressure

	   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 getbidij(vtemp+epso2,vpres,nsp,xsp,bindkj_p) ! get bin Diff coeffs
	   call getbidij(vtemp-epso2,vpres,nsp,xsp,bindkj_m) ! get bin Diff coeffs

	   do k=1,nsp
	      capd(j,k)=(bindkj_p(nspcon+k)-bindkj_m(nspcon+k))/epstemp
	   end do

	end do

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

	return
	end