*HEADING
THERMAL LOADING CONTOUR INTEGRAL CALC., MESH 2, 3-D HEAT TRANSFER    
*NODE
1, 
25, 
361,10.,0.
367,10.,10.
379,-10.,10.
385,-10.,0.
647,10.,40.
659,-10.,40.
**
1001,0.,0.,.5
1025,0.,0.,.5
1361,10.,0.,.5
1367,10.,10.,.5
1379,-10.,10.,.5
1385,-10.,0.,.5
1647,10.,40.,.5
1659,-10.,40.,.5
**
2001,0.,0.,1.
2025,0.,0.,1.
2361,10.,0.,1.
2367,10.,10.,1.
2379,-10.,10.,1.
2385,-10.,0.,1.
2647,10.,40.,1.
2659,-10.,40.,1.
*NSET,NSET=RIGHT,GENERATE
361,367
1361,1367
2361,2367
367,647,20
1367,1647,20
2367,2647,20
*NSET,NSET=YAXIS,GENERATE
379,385
1379,1385
2379,2385
379,659,20
1379,1659,20
2379,2659,20
*NGEN,NSET=CRFR
1,25
*NGEN,NSET=CRMD
1001,1025
*NGEN,NSET=CRBK
2001,2025
*NSET,NSET=CRACK
CRFR,CRMD,CRBK
*NGEN,NSET=OUTER
361,367
367,379
379,385
1361,1367
1367,1379
1379,1385
2361,2367
2367,2379
2379,2385
*NSET,NSET=LOW,GENERATE
367,379
1367,1379
2367,2379
*NGEN,NSET=TOP
647,659
1647,1659
2647,2659
*NFILL,BIAS=.9
LOW,TOP,14,20
*NFILL,SINGULAR=1
CRACK,OUTER,12,30
*NSET,NSET=FRONT,GENERATE
1,385
*NSET,NSET=BACK,GENERATE
2001,2385
*ELEMENT,TYPE=DC3D20
1,1,61,63,3,2001,2061,2063,2003,31,62,33,2,2031,2062,2033,
2002,1001,1061,1063,1003
***
73,367,407,409,369,2367,2407,2409,2369,387,408,389,368,2387,2408,2389,
2368,1367,1407,1409,1369
*ELGEN,ELSET=ELS
1,6,60,1,12,2,6
73,6,2,1,7,40,6
*ELSET,ELSET=PRINT,GENERATE
1,6
*NSET,NSET=XAXIS,GENERATE
25,385,30
1025,1385,30
2025,2385,30
*SOLID SECTION,ELSET=ELS,MATERIAL=HEAT
*MATERIAL,NAME=HEAT
*SPECIFIC HEAT
.1431, 
*DENSITY
.2829, 
*CONDUCTIVITY
7.872E-4, 
*STEP
*HEAT TRANSFER,STEADY STATE
1.,1.
*BOUNDARY
YAXIS,11,,1000.
RIGHT,11,,-1000.
*NODE PRINT
NT, 
*NODE FILE
NT, 
*OUTPUT,FIELD
*NODE OUTPUT
NT,
*END STEP