このルーチンは各ステップの開始時に呼び出され、荷重が適用される面の領域を決定するのに使用することができます。この時、面を通過するフラックスを計算し、別のファイルに保存することもできます。熱計算で対流と輻射による熱フラックスをチェックする場合にはこれが便利なことがあります。
subroutine ufaceload(co,ipkon,kon,lakon,nelemload,sideload,nload)
!
!
! 入力：
!
! co(0..3,1..nk) 節点の座標
! ipkon(*) 場 kon への要素トポロジーポインター
! kon(*) 全要素のトポロジーベクター
! lakon(*) 要素ラベルのベクター
! nelemload(1..2,*) 1: 荷重が適用される要素面
! 2: 環境温度のための節点
! sideload(*) 荷重ラベル
! nload 面分布荷重の数
!
! 各ステップの開始時に呼ばれるユーザールーチン；
! 要素セットの面積を計算し、後で膜や輻射係数の計算に使用することができます。
! 共有ブロックを使うと面積を共有することができます。
!