This subroutine is used for nonuniform radiation heat flux, characterized by distributed load labels of the form RxNUy, cf *RADIATE. The load label can be up to 20 characters long. In particular, y can be used to distinguish different nonuniform radiation patterns. The header and variable description is as follows:
!
subroutine radiate(e,sink,temp,kstep,kinc,time,noel,npt,
& coords,jltyp,field,nfield,loadtype,node,area,vold,mi,
& iemchange)
!
! user subroutine radiate
!
!
! INPUT:
!
! sink present sink temperature
! temp 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
! field currently not used
! nfield currently not used (value = 1)
! loadtype load type label
! node currently not used
! area area covered by the integration point
! 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 freedomm per node (max over all
! nodes) in fields like v(0:mi(2))...
!
! OUTPUT:
!
! e(1) magnitude of the emissivity
! e(2) not used; please do NOT assign any value
! sink sink temperature (need not be defined
! for cavity radiation)
! iemchange = 1 if the emissivity is changed during
! a step, else zero.
!