*HEADING
ONE DIMENSIONAL COUPLED DEMAND WETTABILITY, WITH GEL, LINEAR 3D ELEMENTS 
*** UNITS: M, SEC, NEWTON
*NODE,NSET=ALLN
1,0.,0.
3,.00508,0.
101,0.,.0508
103,.00508,.0508
1001,0.,0.,0.02
1003,.00508,0.,0.02
1101,0.,.0508,0.02
1103,.00508,.0508,0.02
*NGEN,NSET=BOTB
1,3,2
*NGEN,NSET=TOPB
101,103,2
*NFILL,NSET=ALLN
BOTB,TOPB,10,10
*NGEN,NSET=BOTF
1001,1003,2
*NGEN,NSET=TOPF
1101,1103,2
*NFILL,NSET=ALLN
BOTF,TOPF,10,10
*NSET,NSET=LHS,GENERATE
1,101,10
1001,1101,10
*NSET,NSET=PORN,GENERATE
1,101,10
3,103,10
1001,1101,10
1003,1103,10
*ELEMENT,TYPE=C3D8RP,ELSET=BLOCK
1,1,3,13,11,1001,1003,1013,1011
*ELGEN,ELSET=BLOCK
1,10,10,1
*ELSET,ELSET=OUTE
1,3,5,7,9
**********************************************
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
*MATERIAL,NAME=CORE
*ELASTIC
10000.,0.
*POROUS BULK MODULI
 ,2.E9
*PERMEABILITY,SPECIFIC=10000.
3.7E-4, 
*SORPTION
-100000.,.04
-10000.,.05
-4500.,.1
-3500.,.18
-2000.,.45
-1000.,.91
0.,1.
*SORPTION,TYPE=EXSORPTION
-100000.,.09
-10000.,.1
-8000.,.11
-6000.,.18
-4500.,.33
-3000.,.79
-2000.,.91
0.,1.
*GEL
.0005,.0015,1.E8,500.
**********************************************
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,.05
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,-10000.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,5.
*BOUNDARY
BOTB,2
BOTF,2
ALLN,1
ALLN,3
*RESTART,WRITE,FREQUENCY=10
**
*STEP
*SOILS,CONSOLIDATION
1.E-7,1.E-7
*DLOAD
10,P5,-500.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=5,NSET=LHS 
U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,ELSET=OUTE
S, 
E, 
SAT,POR,GELVR,VOIDR
*NODE FILE,FREQUENCY=1,NSET=LHS
U,RF,POR,RVT
*EL FILE,FREQUENCY=1,ELSET=OUTE
S,E
SAT,POR,GELVR,VOIDR
*OUTPUT,FIELD, Frequency=100
*NODE OUTPUT,NSET=LHS
 U,POR
*ELEMENT OUTPUT,ELSET=OUTE
 S,SAT,VOIDR
*OUTPUT,HISTORY,Frequency=5
*NODE OUTPUT,NSET=LHS
 U,RF,POR,RVT
*ELEMENT OUTPUT,ELSET=OUTE
 S22,POR,GELVR
*END STEP
**********************************************
*STEP,INC=100,NLGEOM
*SOILS,CONSOLIDATION,UTOL=10000.
1.,3000.,,300.
*BOUNDARY
1,8,,-100.
3,8,,-100.
1001,8,,-100.
1003,8,,-100.
*NODE FILE,FREQUENCY=1,NSET=LHS
U,RF,POR,RVT
*EL FILE,FREQUENCY=1,ELSET=OUTE
E,
*END STEP