c ... fillup array for 
c ... dimensional dynamic viscosity of mixture at center of cell (i,j)
c ... units  (kg/m.s)
c ... & dimensional heat conductivity of mixture at center of cell (i,j)
c ... units (J/mKs)

	subroutine getmulam(temperature,pressure,Y,ns,ospwt,amuref,
	1                   htcndref,amu,htcond,nx,ny,nz)

	IMPLICIT NONE
 
#include "parameter.com"
#include "drfm_com.com"
#include "defs.h"

	integer i, j, k, ns, nx, ny, nz

	real*8 Y(nsmx), ospwt(nsmx)
	real*8 yspl(nsmx),xspl(nsmx)
	real*8 temp(nx,ny), ys(nsmx,nx,ny)
	real*8 amu(nx,ny), htcond(nx,ny)

	double precision temperature, vtemp, pressure, cpmix, drftc
	double precision rho, wb, drfvs, amuref, htcndref

	i=1
	j=1

	temp(i,j) = temperature
	
	vtemp = temp(i,j)

	do k = 1,ns
	   ys(k,i,j) = Y(k)
	   yspl(k) = ys(k,i,j)
	end do

	call gmlfr(ns,yspl,xspl,wb,ospwt) 

	call drfstc(ns,drfm_xmol,xspl,vtemp,drftc,drfvs,drfm_xml1,
	1    drfm_sig,drfm_epsm,drfm_delm,drfm_chim,
	2    drfm_vstm,drfm_rstm,drfm_tvta,drfm_itox,
	3    drfm_error)

	amu(i,j) = drfvs*1.d-6 ! dimensional
	htcond(i,j) = drftc*1.d-3 !dimensional

c-------------------------------------------------------------------------------
	return
	end