subroutine utemp (temp,nsecpt,kstep,kinc,time,node,coords)
c
      include 'aba_param.inc'
c
      dimension temp(nsecpt),time(2),coords(3)
c Define quadratic temperature gradient in x
      temp(1)=0.01*coords(1)**2
c
      return
      end