このサブルーチンは SxNUy 形式の分布荷重ラベルで特徴づけられる非一様な熱フラックスのために使用されます。*DFLUXを参照してください。荷重ラベルの上限は20文字です。特に y は非一様なフラックスのパターンを判別するために使用できます。ヘッダーと変数の意味は以下の通りです。
subroutine dflux(flux,sol,kstep,kinc,time,noel,npt,coords,jltyp,temp,press,loadtype,area,vold,co,lakonl,konl,ipompc,nodempc,coefmpc,nmpc,ikmpc,ilmpc,iscale,mi)
!
! ユーザーサブルーチン dflux
!
!
! 入力：
!
! sol 現在の温度値
! kstep ステップ番号
! kinc インクリメント番号
! time(1) 現在のステップ時間
! time(2) 現在の総時間
! noel 要素番号
! npt 積分点番号
! coords(1..3) 積分点のグローバル座標
! jltyp 荷重面コード：
! 1 = 体積フラックス
! 11 = 面1
! 12 = 面2
! 13 = 面3
! 14 = 面4
! 15 = 面5
! 16 = 面6
! temp 現在は不使用
! press 現在は不使用
! loadtype 荷重タイプラベル
! area 表面フラックスの場合：積分点によってカバーされる面積
! 体積フラックスの場合：積分点によってカバーされる体積
! vold(0..4,1..nk) 全節点で解かれる場
! 0: 温度
! 1: グローバル X 方向の変位
! 2: グローバル Y 方向の変位
! 3: グローバル Z 方向の変位
! 4: 静圧
! co(3,1..nk) 全節点の座標
! 1: グローバル X 方向の座標
! 2: グローバル Y 方向の座標
! 3: グローバル Z 方向の座標
! lakonl 要素ラベル
! konl(1..20) 要素に属する節点
! ipompc(1..nmpc) ipompc(i) は場 nodempc の MPC i の最初の項を指す
! nodempc(1,*) MPC の項の節点番号
! nodempc(2,*) MPC の項の座標方向
! nodempc(3,*) 0以外の場合：場 nodempc の MPC の次の項を指す
! 0の場合：MPC 定義終了
! coefmpc(*) MPC の項の係数
! nmpc MPC の数
! ikmpc(1..nmpc) 順序付けられたMPCのグローバルな自由度
! グローバルな自由度は MPC の依存項の 8*(node-1)+direction（direction = 0：温度, 1-3：変位, 4：静圧, 5-7：回転）
! ilmpc(1..nmpc) ilmpc(i) は ikmpc(i) での参照番号に対応する MPC 番号
! mi(1) 要素あたりの積分点の最大数（全要素にわたる最大値）
! mi(2) 場での節点あたりの最大自由度（全節点にわたる最大値）。例 v(0:mi(2))。
!
! 出力：
!
! flux(1) フラックスの大きさ
! flux(2) 不使用：値を代入しないでください
! iscale 静計算でのステップ時間より小さなインクリメントでフラックスをスケールするかどうか
! 0：スケールしない
! 1：スケールする（デフォルト）
!