c------------------------------------------------------------------------------
c     USDFLD subroutine
c------------------------------------------------------------------------------

      subroutine usdfld(field,statev,pnewdt,direct,t,celent,
     1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer,
     2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo,
     3 laccfla)
c
      include 'aba_param.inc'
c
      character*80 cmname,orname
      character*3  flgray(15)
      dimension field(nfield),statev(nstatv),direct(3,3),
     1 t(3,3),time(2)
      dimension array(15),jarray(15),jmac(*),jmatyp(*),
     1 coord(*)

c--------------------------------------------------------------------------

	if (kstep .eq. 1 .and. kinc .eq. 1) then
	
		x = coord(1)
		y = coord(2)
		z = coord(3)
		
		field(1) = 0.0d0 + z / 1.0d0 * 1.0d0
		
		statev = field
		
	else
		
		field = statev
		
	end if

c--------------------------------------------------------------------------

      return
      end