SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
     1                 JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION COORDS(3)
      CHARACTER*80 SNAME
      write(*,*) COORDS(1), COORDS(2), COORDS(3) 
         F=-400/(SQRT((COORDS(1)**2)+(COORDS(2)**2)))
      RETURN
      END