*HEADING
 HEAT, GAP CONDUCTANCE, MF DEP USER SUBROUTINE - UTFVMHM
*NODE,NSET=PLANEA
1, 
3,,,1.
5,,1.,1.
7,,1.
*NCOPY,CHANGE NUMBER=100,OLD SET=PLANEA,SHIFT,NEW SET=PLANEB
1.,0.,0.

*NCOPY,CHANGE NUMBER=200,OLD SET=PLANEA,SHIFT,NEW SET=PLANEC
1.,0.,0.

*NCOPY,CHANGE NUMBER=300,OLD SET=PLANEA,SHIFT,NEW SET=PLANED
2.,0.,0.

*NSET,NSET=BCPLANE
 PLANEB, PLANEC
*ELEMENT,TYPE=DCC3D8,ELSET=ONE
1,101,103,105,107,1,3,5,7
*ELEMENT,TYPE=DCC3D8,ELSET=TWO
2,301,303,305,307,201,203,205,207
*surface, NAME=SURF1
ONE, S1
*surface, NAME=SURF2
TWO, S2
*CONTACT PAIR, INTERACTION=GAP
SURF1, SURF2
*ELSET,ELSET=BOTH
ONE,TWO
*SOLID SECTION,ELSET=BOTH,MATERIAL=BLOCK
*MATERIAL,NAME=BLOCK
*CONDUCTIVITY
150., 
*DENSITY
1., 
*SPECIFIC HEAT
1., 
*SURFACE INTERACTION, NAME=GAP
*GAP CONDUCTANCE
15.,,,25.
15.,1.,,25.
150.,,,50.
150.,1.,,50.
300.,,,75.
300.,1.,,75.
3000.,,,100.
3000.,1.,,100.
*NSET,NSET=EDGE1,GEN
1,301,100
*NSET,NSET=EDGE3,GEN
3,303,100
*NSET,NSET=EDGE5,GEN
5,305,100
*NSET,NSET=EDGE7,GEN
7,307,100
*STEP,AMP=STEP
*HEAT TRANSFER,STEADY STATE
.25,1.,.25,.25
*BOUNDARY
PLANEA,11,11,1000.
PLANED,11,11
*MASS FLOW RATE,USER
EDGE1,
EDGE3,
EDGE5,
EDGE7,
*EL PRINT,F=0
*NODE PRINT
NT, 
*NODE FILE
NT, 
*OUTPUT,FIELD
*NODE OUTPUT
NT,
*CONTACT PRINT
*CONTACT FILE
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
*END STEP