This subroutine is used for user-defined concentrated load, characterized by the parameter USER on the *CLOAD card. The header and variable description is as follows:
subroutine cload(xload,kstep,kinc,time,node,idof,coords,vold,
& mi,ntrans,trab,inotr,veold)
!
! user subroutine cload
!
!
! INPUT:
!
! kstep step number
! kinc increment number
! time(1) current step time
! time(2) current total time
! node node number
! idof degree of freedom
! coords(1..3) global coordinates of the node
! vold(0..mi(2)
! ,1..nk) solution field in all nodes (for modal
! dynamics: in all nodes for which output
! was requested or a force was applied)
! 0: temperature
! 1: displacement in global x-direction
! 2: displacement in global y-direction
! 3: displacement in global z-direction
! 4: static pressure
! mi(1) max # of integration points per element (max
! over all elements)
! mi(2) max degree of freedomm per node (max over all
! nodes) in fields like v(0:mi(2))...
! veold(0..3,1..nk) derivative of the solution field w.r.t.
! time in all nodes(for modal
! dynamics: in all nodes for which output
! was requested or a force was applied)
! 0: temperature rate
! 1: velocity in global x-direction
! 2: velocity in global y-direction
! 3: velocity in global z-direction
! ntrans number of transform definitions
! trab(1..6,i) coordinates of two points defining transform i
! trab(7,i) -1: cylindrical transformation
! 1: rectangular transformation
! inotr(1,j) transformation number applied to node j
! inotr(2,j) a SPC in a node j in which a transformation
! applied corresponds to a MPC. inotr(2,j)
! contains the number of a new node generated
! for the inhomogeneous part of the MPC
!
! OUTPUT:
!
! xload concentrated load in direction idof of node
! "node" (global coordinates)
!