subroutine usdfld(field,statev,pnewdt,direct,t,celent,time,dtime,
     1 cmname,orname,nfield,nstatv,noel,npt,layer,kspt,kstep,kinc,
     2 ndi,nshr,coord,jmac,jmtyp,matlayo,laccflg)
c
      include 'aba_param.inc'
c
      character*80 cmname,orname
      character*3  flgray(15)
      dimension field(nfield),statev(nstatv),direct(3,3),t(3,3),time(2),
     * coord(*),jmac(*),jmtyp(*)
      dimension array(15),jarray(15)
c
c Get temperatures from previous increment
      call getvrm('TEMP',array,jarray,flgray,jrcd,
     $     jmac, jmtyp, matlayo, laccflg)
      call getvrm('ER',array,jarray,flgray,jrcd,
     $     jmac, jmtyp, matlayo, laccflg)
      call getvrm('S',array,jarray,flgray,jrcd,
     $     jmac, jmtyp, matlayo, laccflg)
      temp = array(1)
      epsdot11=array(2)
      epsdot22=array(3)
      epsdot33=array(4)
      epsdot12=array(5)
      epsdot13=array(6)
      epsdot23=array(7)
      stress11=array(8)
      stress22=array(9)
      stress33=array(10)
      stress12=array(11)
      stress13=array(12)
      stress23=array(13)

      field(1) = temp
      field(2) = epsdot11
      field(3) = epsdot22
      field(4) = epsdot33
      field(5) = epsdot12
      field(6) = epsdot13
      field(7) = epsdot23
      field(8) = stress11
      field(9) = stress22
      field(10) = stress33
      field(11) = stress12
      field(12) = stress13
      field(13) = stress23
      
      statev(1) = epsdot11*stress11 + epsdot22*stress22 + epsdot33*stress33


      RETURN
      END