c
c
      subroutine umatht(u,dudt,dudg,flux,dfdt,dfdg,statev,temp,
     $     dtemp,dtemdx,time,dtime,predef,dpred,cmname,ntgrd,nstatv,
     $     props,nprops,coords,pnewdt,noel,npt,layer,kspt,kstep,kinc)
c
      include 'aba_param.inc'
c
      character*80 cmname
c
      dimension dudg(ntgrd),flux(ntgrd),dfdt(ntgrd),
     $     dfdg(ntgrd,ntgrd),statev(nstatv),dtemdx(ntgrd),time(2),
     $     predef(1),dpred(1),props(nprops),coords(3)
c
c
      cond = props(1)
      specht = props(2)
c
c                    input specific heat
      dudt = specht
      du = dudt*dtemp
      u = u+du
c
c                    input flux = -[k]*{dtemdx}
      do i=1, ntgrd
         flux(i) = -cond*dtemdx(i)
      end do
c              
c                    input isotropic conductivity
c
      do i=1, ntgrd
         dfdg(i,i) = -cond
      end do
c
      return
      end