SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
     1                 JLTYP,SNAME)

      INCLUDE 'ABA_PARAM.INC'

      DIMENSION COORDS(3)
      CHARACTER*80 SNAME

      F=-400/(SQRT((COORDS(1)**2)+(COORDS(2)**2)))

      RETURN
      END