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

	IMPLICIT NONE

	include 'parameter.par'
	include 'common.com'

	integer j, k, nspcon

	double precision vtemp,vpres

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,vpres,nsp,xsp,bindkj)     ! get bin Diff coeffs

	   do k=1,nsp
	      capd(j,k) = bindkj(nspcon+k)		! save dimless D_k_ns 
	   end do

	end do

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

	return
	end