*HEADING
Submodel analysis: 3 x 3 mesh with c3d8ph.     File name: 
*NODE,NSET=ALL
 1, 1.95
61, 4.
67, 4.,3.45
 7, 1.95,3.45
101, 1.95,, .5
161, 4.,, .5
167, 4.,3.45, .5
107, 1.95,3.45, .5
201, 1.95,,1.0
261, 4.,,1.0
267, 4.,3.45,1.0
207, 1.95,3.45,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
11,61,10
121,161,20
211,261,10
**1,61,10
**101,161,20
**201,261,10
*NSET,NSET=FIX2,GEN
61,67
161,167,2
261,267
*NSET,NSET=FIX3,GEN
7,67,10
107,167,20
207,267,10
*NSET,NSET=FIX4,GEN
1,7
101,107,2
201,207
*NSET,NSET=ALL
ALL1,ALL2,ALL3
*NSET,NSET=NALL
FIX1,FIX2,FIX3,FIX4
*ELEMENT, TYPE=C3D8PH, 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, 
*INITIAL CONDITIONS,TYPE=RATIO
ALL,1.
*SUBMODEL
NALL,
*RESTART,WRITE
*STEP
*SOILS, CONSOLIDATION
 1. , 1.
*ENERGY FILE
*EL FILE
 LOADS,
 IVOL,
 EVOL,
*NODE FILE
U,RF,POR,RVT
*EL PRINT,TOTALS=YES
 LOADS,
 IVOL,
 EVOL,
*NODE PRINT
 RF, 
 POR, 
*BOUNDARY
 ALL1,  3
 ALL3,  3
FIX1, 2
*BOUNDARY, SUBMODEL, STEP=1
FIX3,1,2
FIX3,8
FIX4,1,2
FIX4,8
*FLOW, OP=NEW
EFLOW2, Q4, 14.7, 2.E-3
*END STEP