SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,KSLAY,KSPT, * COORDS,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