*HEADING
Global analysis: 3 x 3 mesh with c3d8pt 
*NODE,NSET=ALL
 1, 1.
61, 4.
67, 4.,5.
 7, 1.,5.
101, 1.,, .5
161, 4.,, .5
167, 4.,5., .5
107, 1.,5., .5
201, 1.,,1.0
261, 4.,,1.0
267, 4.,5.,1.0
207, 1.,5.,1.0
*NGEN,NSET=NL1
1,7,1
*NGEN,NSET=NR1
61,67,1
*NFILL,NSET=ALL1
NL1,NR1,6,10
*NGEN,NSET=NL2
101,107,1
*NGEN,NSET=NR2
161,167,1
*NFILL,NSET=ALL2
NL2,NR2,6,10
*NGEN,NSET=NL3
201,207,1
*NGEN,NSET=NR3
261,267,1
*NFILL,NSET=ALL3
NL3,NR3,6,10
*NSET,NSET=FIX1,GEN
1,61,10
101,161,10
201,261,10
*NSET,NSET=FIX2,GEN
61,67
161,167
261,267
*NSET,NSET=FIX3,GEN
7,67,10
107,167,10
207,267,10
*NSET,NSET=FIX4,GEN
1,7
101,107
201,207
*NSET,NSET=ALL
ALL1,ALL2,ALL3
*NSET,NSET=NALL
FIX1,FIX2,FIX3,FIX4
*ELEMENT, TYPE=C3D8PT, ELSET=EALL
1,1,21,23,3, 201,221,223,203
*ELGEN,ELSET=EALL
1,3,2,1,3,20,3
*ELSET,ELSET=EFLOW2
7,8,9
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 1.E8, 
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 
*DENSITY
  1.4142, 
*Specific heat
 10.0
*Conductivity
 1.0
*expansion
1e-06
*DENSITY, pore fluid
 1.4142, 
*Specific heat, pore fluid
 10.0
*Conductivity , pore fluid
 1.0
*expansion, pore fluid
 1e-06
*INITIAL CONDITIONS,TYPE=RATIO
ALL,1.
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL,0.0
*RESTART,WRITE
*STEP,nlgeom=yes,unsymm=yes
*SOILS, CONSOLIDATION
 1. , 1.
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*NODE OUTPUT
U,RF,POR,RVT,NT
*Element Output
S
*EL PRINT,TOTALS=YES
 LOADS,
 IVOL,
 EVOL,
*NODE PRINT
 RF, 
 POR, 
*BOUNDARY,OP=NEW
 ALL1,  3
 ALL3,  3
 FIX1, 2
 FIX4, 1
 FIX4, 8
 ALL, 11, 11, 0.0
*FLOW, OP=NEW
 EFLOW2, Q4, 14.7, 2.E-3
*END STEP