This subroutine is used for nonuniform heat flux, characterized by distributed load labels of the form SxNUy, cf *DFLUX. The load label can be up to 20 characters long. In particular, y can be used to distinguish different nonuniform flux patterns. The header and variable description is as follows:
! ! user subroutine dflux ! ! ! INPUT: ! ! sol current temperature value ! kstep step number ! kinc increment number ! time(1) current step time ! time(2) current total time ! noel element number ! npt integration point number ! coords(1..3) global coordinates of the integration point ! jltyp loading face kode: ! 11 = face 1 ! 12 = face 2 ! 13 = face 3 ! 14 = face 4 ! 15 = face 5 ! 16 = face 6 ! temp currently not used ! press currently not used ! loadtype load type label ! vold(0..3,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 ! ! OUTPUT: ! ! flux(1) magnitude of the flux ! flux(2) not used; please do NOT assign any value !