C
C  User subroutine to define temperature

      subroutine upress(press,kstep,kinc,time,node,coords)
c
C
      include 'aba_param.inc'
c
      parameter(const=1D4,rlen=7.0,div=3.0)
      dimension coords(3),time(2)
c
      pos=(coords(1)+coords(2)+coords(3))/rlen
      press=const*pos/div
c   
c
      return
      end