*HEADING
 DSLOAD, DSFLUX, SFILM, AND SRADIATE FOR COUPLED TEMP-DISP [C3D10MT]
 RESULTS MATCH EC3ATKDK
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
 1,0.,0.,0.
 3,7.,0.,0.
 19,0.,7.,0.
 21,7.,7.,0.
 7,0.,0.,-7.
 9,7.,0.,-7.
 25,0.,7.,-7.
 27,7.,7.,-7.
*NGEN,NSET=ALL
 1,3
 1,7,3
 3,9,3
 7,9
 1,19,9
 3,21,9
 7,25,9
 9,27,9
 19,25,3
 21,27,3
 19,21
 25,27
 8,26,9
 2,20,9
 2,8,3
 10,16,3
 20,26,3
 12,18,3
*ELEMENT,TYPE=C3D10MPT,ELSET=ELEMS
 10, 7,  9, 25,  1,  8, 17, 16,  4,  5, 13
 20, 1,  9, 21,  3,  5, 15, 11,  2,  6, 12
 30, 9, 27, 25, 21, 18, 26, 17, 15, 24, 23
 40, 1, 21, 25, 19, 11, 23, 13, 10, 20, 22
 50, 1, 21,  9, 25, 11, 15,  5, 13, 23, 17
*SOLID SECTION,MATERIAL=ELMAT, ELSET=ELEMS
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-460.,STEFAN BOLTZMANN=5.0e-8
*MATERIAL,NAME=ELMAT
*Permeability, specific=10.
0.001,0.
*ELASTIC
 3.E6, 0.3
*CONDUCTIVITY
 3.77E-5,
*SPECIFIC HEAT
.39,
*DENSITY
 82.9, 
*EXPANSION
 .0, 
*CONDUCTIVITY, PORE FLUID
 3.77E-5,
*SPECIFIC HEAT, PORE FLUID
.39,
*DENSITY, PORE FLUID
 82.9, 
*EXPANSION, PORE FLUID
 .0, 
*NSET, NSET=NTOP
 19, 20, 21, 10, 11, 12, 1, 2, 3
*SURFACE, NAME=TOP
 20, S4
 40, S2
*NSET, NSET=NBOT
 27, 26, 25, 18, 17, 16, 7, 8, 9
*SURFACE, NAME=BOTTOM
 10, S1
 30, S1
*NSET, NSET=NSIDE_1
 25, 26, 27, 22, 23, 24, 19, 20, 21
*SURFACE, NAME=SIDE_1
 30, S3
 40, S3
*NSET, NSET=NSIDE_2
 1, 2, 3, 4, 5, 6, 7, 8, 9
*SURFACE, NAME=SIDE_2
 10, S2
 20, S2
*NSET, NSET=NSIDE_3
 25, 22, 19, 16, 13, 10, 1, 4, 7
*SURFACE, NAME=SIDE_3
 10, S4
 40, S4
*NSET, NSET=NSIDE_4
 21, 24, 27, 12, 15, 18, 3, 6, 9
*SURFACE, NAME=SIDE_4
 20, S3
 30, S2
*INITIAL CONDITIONS,TYPE=TEMPERATURE
 ALL,0.
*INITIAL CONDITIONS,TYPE=RATIO
 ALL,1.
*INITIAL CONDITIONS,TYPE=PORE
 ALL,0.
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 TOP, P, 3000.
*DSFLUX, OP=NEW
 TOP, S, .3
*BOUNDARY, OP=NEW
 ALL,1,3
 NBOT,11,11
 ALL, 8, 8, 0.0
*EL PRINT
 LOADS,FLUXS
 FILM, RAD
*NODE PRINT,TOTALS=YES
 RF,
 NT,
 RFL,
*OUTPUT,FIELD, VARIABLE=PRESELECT
*OUTPUT,FIELD
*NODE OUTPUT
 RF,
 NT,
 RFL,
*ELEMENT OUTPUT
 HFL
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 BOTTOM,P,3000.
*DSFLUX, OP=NEW
 BOTTOM,S,0.3
*BOUNDARY, OP=NEW
 ALL,1,3
 NTOP,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_1,P,3000.
*DSFLUX, OP=NEW
 SIDE_1,S,0.3
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_2,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_2,P,3000.
*DSFLUX, OP=NEW
 SIDE_2,S,0.3
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_1,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_3,P,3000.
*DSFLUX, OP=NEW
 SIDE_3,S,0.3
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_4,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_4,P,3000.
*DSFLUX, OP=NEW
 SIDE_4,S,0.3
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_3,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 BOTTOM,HP,3000.,0.,-7.
*DSFLUX, OP=NEW
*SFILM, OP=NEW
 BOTTOM,F,75.,0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NTOP,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 TOP,HP,3000.,0.,-7.
*SFILM, OP=NEW
 TOP, F, 75.0, 0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NBOT,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_2,HP,3000.,0.,-7.
*SFILM, OP=NEW
 SIDE_2,F, 75.0, 0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_1,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_4,HP,3000.,0.,-7.
*SFILM, OP=NEW
 SIDE_4,F, 75.0, 0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_3,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_1,HP,3000.,0.,-7.
*SFILM, OP=NEW
 SIDE_1,F,75.0,0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_2,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
 SIDE_3,HP,3000.,0.,-7.
*SFILM, OP=NEW
 SIDE_3,F,75.0,0.103
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_4,11,11
 ALL, 8, 8, 0.0
*END STEP
*STEP
*SOILS
1.0,1.0
*DSLOAD, OP=NEW
*SFILM, OP=NEW
*SRADIATE, OP=NEW
 BOTTOM,R,75.,1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NTOP,11,11
 ALL, 8, 8, 0.0
*NODE PRINT,TOTALS=YES
 NT,
 RFL,
*OUTPUT,FIELD
*NODE OUTPUT
 NT,
 RFL,
*END STEP
**
*STEP
*SOILS
1.0,1.0
*SRADIATE, OP=NEW
 TOP, R, 75.0, 1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NBOT,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*SRADIATE, OP=NEW
 SIDE_2,R, 75.0, 1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_1,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*SRADIATE, OP=NEW
 SIDE_4,R, 75.0, 1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_3,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*SRADIATE, OP=NEW
 SIDE_1,R,75.0,1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_2,11,11
 ALL, 8, 8, 0.0
*END STEP
**
*STEP
*SOILS
1.0,1.0
*SRADIATE, OP=NEW
 SIDE_3,R,75.0,1.e-6
*BOUNDARY, OP=NEW
 ALL,1,3
 NSIDE_4,11,11
 ALL, 8, 8, 0.0
*END STEP