*HEADING
Global analysis: 3 x 3 mesh with c3d20rp.     File name:
*NODE,NSET=ALL
*NGEN,NSET=NL1
*NGEN,NSET=NR1
*NFILL,NSET=ALL1
NL1,NR1,6,10
*NGEN,NSET=NL2
*NGEN,NSET=NR2
*NFILL,NSET=ALL2
NL2,NR2,6,10
*NGEN,NSET=NL3
*NGEN,NSET=NR3
*NFILL,NSET=ALL3
NL3,NR3,6,10
*NSET,NSET=FIX1,GEN
*NSET,NSET=FIX2,GEN
*NSET,NSET=FIX3,GEN
*NSET,NSET=FIX4,GEN
*NSET,NSET=ALL
ALL1,ALL2,ALL3
*NSET,NSET=NALL
FIX1,FIX2,FIX3,FIX4
*ELEMENT, TYPE=C3D20RP, ELSET=EALL
*ELGEN,ELSET=EALL
*ELSET,ELSET=EFLOW2
*element,type=r3d3,elset=rigid
*rigid body,elset=rigid,refnod=10000
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 1.E8,
*PERMEABILITY
 1.E-5,
*DENSITY
  1.4142,
*INITIAL CONDITIONS,TYPE=RATIO
ALL,1.
*RESTART,WRITE
*STEP
*SOILS, CONSOLIDATION
*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