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

	IMPLICIT NONE

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

	integer j, k
	
	double precision amuv,vtemp

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 getmxvis(vtemp,nsp,xsp,amuv)
	   amu(j) = amuv
	end do

	return
	end

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