With this subroutine the user can define a temperature field. It is triggered bythe parameter USER on the *TEMPERATURE card. The header and variable description is as follows:
subroutine utemp(temp,msecpt,kstep,kinc,time,node,coords,vold,
& mi)
!
! user subroutine utemp
!
!
! INPUT:
!
! msecpt number of temperature 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 freedomm per node (max over all
! nodes) in fields like v(0:mi(2))...
!
! OUTPUT:
!
! temp(1..msecpt) temperature in the node
!