This subroutine is used for user-defined concentrated heat flux, characterized by the parameter USER on the *CFLUX card. The header and variable description is as follows:
subroutine cflux(flux,msecpt,kstep,kinc,time,node,coords,vold,
& mi)
!
! user subroutine cflux
!
!
! INPUT:
!
! msecpt number of flux values (for volume elements:1)
! kstep step number
! kinc increment number
! time(1) current step time
! time(2) current total time
! node node number
! coords(1..3) global coordinates of the node
! vold(0..4,1..nk) solution field in all nodes
! 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 freedom per node (max over all
! nodes) in fields like v(0:mi(2))...
!
! OUTPUT:
!
! flux(1..msecpt) concentrated flux in the node
!