*HEADING
 USER SUBROUTINE GAPCON, MF DEP
*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
** tabulated clearance dependent gap conductance not present
*GAP CONDUCTANCE,USER
*AMPLITUDE,NAME=RISE
0.,0.,1.,1.
*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,AMP=RISE
EDGE1,0.,70.71068,70.71068
EDGE3,0.,70.71068,-70.71068
EDGE5,0.,-70.71068,-70.71068
EDGE7,0.,-70.71068,70.71068
*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