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.D0
      P=0.
      IF(JLTYP.EQ.22) P=10.D3
      IF(JLTYP.EQ.24) P=30.D3
      F=P* 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