*HEADING
 DSLOAD, DSFLUX, SFILM, AND SRADIATE FOR COUPLED TEMP-DISP [C3D4PT]
*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=C3D4PT,ELSET=ELEMS
 10, 7,  9, 25,  1
 20, 1,  9, 21,  3
 30, 9, 27, 25, 21
 40, 1, 21, 25, 19
 50, 1, 21,  9, 25
*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