c***************************************************************************
c ... fillup array for 
c ... dimensional dynamic viscosity of mixture at center of cell (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 getmu
 
	include 'parameter.par'
#include "precision.com"
	include 'common.com'

	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 getmxvis(vtemp+epso2,nsp,xsp,amuv_p)
	   call getmxvis(vtemp-epso2,nsp,xsp,amuv_m)
	   amu(j) = (amuv_p - amuv_m)/epstemp
	end do

	return
	end

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