Next: Temperature (utemp.f) Up: User-defined loading Previous: Boundary conditions(uboun.f) Contents

###
Heat radiation (radiate.f)

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.
    !