This routine is called at the beginning of each step and can be used to determine the area of faces on which loading is applied. In that way the flux through the face can be calculated and stored in an extra file. This can be beneficial for thermal calculations to check the heat flux due to convection and radiation.

      subroutine ufaceload(co,ipkon,kon,lakon,
!
!
!     INPUT:
!
!     co(0..3,1..nk)     coordinates of the nodes
!     ipkon(*)           element topology pointer into field kon
!     kon(*)             topology vector of all elements
!     lakon(*)           vector with elements labels