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 !