SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
     1                 JLTYP,SNAME)
      INCLUDE 'ABA_PARAM.INC'
      DIMENSION COORDS(3)
      CHARACTER*80 SNAME
C NOTE THAT COORDS(3) IS THE ANGULAR COORD IN DEGREES
      PI=2.*ASIN(1.D0)
      THETA=PI*COORDS(3)/180.
      F=-COORDS(1)*COS(THETA)
C      WRITE(6,9991)NOEL,NPT,JLTYP,F,PI,THETA,COORDS
C9991  FORMAT(' NOEL,NPT,JLTYP,F,PI,THETA,COORDS',3I5,1P6G10.3)
      RETURN
      END