*HEADING
SFLOW, DSFLOW, OP=NEW  FOR COUPLED TEMP-DISPLACEMENT-PORE PRESSURE ELEMENTS [CPE4RPHT]
LOAD LABELS USED: S, Q
Results should match with flow load steps of ece4pfdl.inp (plane strain elements with pore-pressure)
**###################################################################
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
1, 1., 0.
2, 4., 0.
3, 4., 5.
4, 1., 5.
*NSET, NSET=SIDE1
1, 2
*NSET, NSET=SIDE2
2, 3
*NSET, NSET=SIDE3
3, 4
*NSET, NSET=SIDE4
4, 1
*NSET, NSET=NALL
1,2,3,4
*SURFACE, TYPE=ELEMENT, NAME=SURF1
1, S1
*SURFACE, TYPE=ELEMENT, NAME=SURF2
1, S2
*SURFACE, TYPE=ELEMENT, NAME=SURF3
1, S3
*SURFACE, TYPE=ELEMENT, NAME=SURF4
1, S4
*ELEMENT,TYPE=CPE4RPHT, ELSET=ALL_ELE
1, 1,2,3,4
*SOLID SECTION,MATERIAL=ELMAT, ELSET=ALL_ELE
*MATERIAL,NAME=ELMAT
**#################
** Properties for the solid
**#################
*ELASTIC
1.0E+08,0.0
*DENSITY
1.4142, 
*PERMEABILITY, SPECIFIC=1.0
1.0E-05
*CONDUCTIVITY
0.1
*SPECIFIC HEAT
0.39,
**####################
** Initial Conditions
**####################
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL, 0.
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*INITIAL CONDITIONS,TYPE=PORE
NALL,0.
**###################################################
** STEPS 1 TO 4:  *SFLOW   ->  Specify seepage coefficient and sink pore pressure on a surface to model pore fluid flow as a function of current pore pressure
** STEPS 5 TO 8:  *DSFLOW  ->  Specify seepage velocity on a surface
**###################################################
**############
**   STEP 1
**############
*STEP
SFLOW through face 1
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE3, 8
NALL, 11, 11
*SFLOW, OP=NEW
SURF1, Q, 14.7, 2.E-5
*EL PRINT
LOADS,
*NODE PRINT,TOTALS=YES
RF, 
*NODE PRINT
POR,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
POR
*OUTPUT,HISTORY
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP
**############
**   STEP 2
**############
*STEP
SFLOW through face 2
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE4, 8
NALL, 11, 11
*SFLOW, OP=NEW
SURF2, Q, 14.7, 2.E-5
*END STEP
**############
**   STEP 3
**############
*STEP
SFLOW through face 3
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE1, 8
NALL, 11, 11
*SFLOW, OP=NEW
SURF3, Q, 14.7, 2.E-5
*END STEP
**############
**   STEP 4
**############
*STEP
SFLOW through face 4
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY,OP=NEW
NALL, 1,2
SIDE2, 8
NALL, 11, 11
*SFLOW, OP=NEW
SURF4, Q, 14.7, 2.E-5
*END STEP
**#######################################
**#######################################
**############
**   STEP 5
**############
*STEP
DSFLOW through face 1;
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE3, 8
NALL, 11, 11
*SFLOW, OP=NEW
*DSFLOW, OP=NEW
SURF1, S, 3.E-5
*END STEP
**############
**   STEP 6
**############
*STEP
DSFLOW through face 2;
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE4, 8
NALL, 11, 11
*DSFLOW, OP=NEW
SURF2, S, 3.E-5
*END STEP
**############
**   STEP 7
**############
*STEP
DSFLOW through face 3
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE1, 8
NALL, 11, 11
*DSFLOW, OP=NEW
SURF3, S, 3.E-5
*END STEP
**############
**   STEP 8
**############
*STEP
DSFLOW through face 4
*SOILS,CONSOLIDATION
1., 1.
*BOUNDARY, OP=NEW
NALL, 1,2
SIDE2, 8
NALL, 11, 11
*DSFLOW, OP=NEW
SURF4, S, 3.E-5
*END STEP