*HEADING
 USER SUBROUTINE GAPCON, FV 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=DC3D8,ELSET=ONE
1,101,103,105,107,1,3,5,7
*ELEMENT,TYPE=DC3D8,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., 
*SURFACE INTERACTION, NAME=GAP
** tabulated clearance dependent gap conductance not present
*GAP CONDUCTANCE,USER
*AMPLITUDE,NAME=RISE
0.,0.,1.,1.
*STEP,AMP=STEP
*HEAT TRANSFER,STEADY STATE
.25,1.,.25,.25
*BOUNDARY
PLANEA,11,11,1000.
PLANED,11,11
*FIELD,VARIABLE=1,AMP=RISE
PLANEB,100.
PLANEC,100.
*EL PRINT
COORD,FV
*EL FILE
COORD,FV
*OUTPUT,FIELD
*ELEMENT OUTPUT
COORD,FV
*NODE PRINT
NT, 
*NODE FILE
NT, 
*OUTPUT,FIELD
*NODE OUTPUT
NT,
*CONTACT PRINT
*CONTACT FILE
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
*END STEP