*HEADING
TWO DIMENSIONAL CYCLIC DEMAND WETTABILITY, WITH GEL      
*** UNITS: M, SEC, NEWTON
*NODE,NSET=ALLN
1,0.,0.
21,.0508,0.
801,0.,.0508
821,.0508,.0508
*NGEN,NSET=BOT
1,21,1
*NGEN,NSET=TOP
801,821,1
*NFILL,NSET=ALLN
BOT,TOP,20,40
*NSET,NSET=OUTN
1,165,329,493,657,821
*NSET,NSET=QA_TEST_NODE
 1, 493, 821
*ELEMENT,TYPE=CPE8RP,ELSET=BLOCK
1,1,3,83,81,2,43,82,41
*ELGEN,ELSET=BLOCK
1,10,2,1,10,80,10
*ELSET,ELSET=OUTE
1,23,45,67,89,100
*ELSET,ELSET=QA_TEST_ELEM
 1, 67, 100
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
.02, 
*MATERIAL,NAME=CORE
*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.
*elastic
100000.0, 0.0
*POROUS BULK MODULI
 ,100000000.
*GEL
.0005,.0015,1.E8,500.
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,.05
*NSET,NSET=PORN,GENERATE
1,21,2
81,101,2
161,181,2
241,261,2
321,341,2
401,421,2
481,501,2
561,581,2
641,661,2
721,741,2
801,821,2
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN,-10000.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,5.
*RESTART,WRITE,FREQUENCY=10
*STEP,AMPLITUDE=STEP
*SOILS,CONSOLIDATION,UTOL=10000.
1.,600.,,100.0
*BOUNDARY
ALLN,PINNED
1,8,,-100.
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
.01,1.,,1.E-6
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=10,NSET=PORN 
POR,RVF,RVT
*EL PRINT,FREQUENCY=10
SAT,POR,VOIDR,GELVR
*NODE FILE,FREQUENCY=5,NSET=QA_TEST_NODE
 POR,RVF,RVT
*EL FILE,FREQUENCY=5,ELSET=QA_TEST_ELEM
 SAT,POR,VOIDR,GELVR
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=OUTN
 POR,RVF,RVT
*ELEMENT OUTPUT,ELSET=OUTE
 SAT,POR,VOIDR,GELVR
*OUTPUT,HISTORY,FREQ=5
*NODE OUTPUT,NSET=QA_TEST_NODE
 POR,RVF,RVT
*ELEMENT OUTPUT,ELSET=QA_TEST_ELEM
 SAT,POR,VOIDR,GELVR
*END STEP
*STEP,AMPLITUDE=STEP
*SOILS,CONSOLIDATION,UTOL=10000.
1.,600.,,100.0
*BOUNDARY
ALLN,PINNED
1,8,,-10000.
*END STEP
*STEP,AMPLITUDE=STEP
*SOILS,CONSOLIDATION,UTOL=4000.
.1,800.,,100.0
*BOUNDARY
ALLN,PINNED
1,8,,-100.
*END STEP