*HEADING
Global analysis: 3 x 3 mesh with c3d20rp.     File name: 
*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
10000,
1000,
1001, 0.1
1002, 0.1, 0.1
*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=C3D20RP, ELSET=EALL
1,1,21,23,3, 201,221,223,203, 11,22,13,2, 211,222,213,
  202, 101,121,123,103
*ELGEN,ELSET=EALL
1,3,2,1,3,20,3
*ELSET,ELSET=EFLOW2
7,8,9
*element,type=r3d3,elset=rigid
1001,1000,1001,1002
*rigid body,elset=rigid,refnod=10000
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 1.E8, 
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 
*DENSITY
  1.4142, 
*INITIAL CONDITIONS,TYPE=RATIO
ALL,1.
*RESTART,WRITE
*STEP
*SOILS, CONSOLIDATION
 1. , 1.
*ENERGY FILE
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*EL FILE
 LOADS,
 IVOL,
 EVOL,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
IVOL,
EVOL,
*NODE FILE
U,RF,POR,RVT
*OUTPUT,FIELD
*NODE OUTPUT
U,RF,POR,RVT
*EL PRINT,TOTALS=YES
 LOADS,
 IVOL,
 EVOL,
*NODE PRINT
 RF, 
 POR, 
*BOUNDARY,OP=NEW
 ALL1,  3
 ALL3,  3
FIX1, 2
FIX4, 1
FIX4, 8
10000,1,6
*FLOW, OP=NEW
EFLOW2, Q4, 14.7, 2.E-3
*END STEP