*HEADING
ONE DIMENSIONAL WICKING PROBLEM, COUPLED      
*** UNITS: M, TON, SEC, KN
*NODE,NSET=ALLN
 1,0.,0.,0.
 3,.1,0.,0.
 101,0.,0.,1.
 103,.1,0.,1.
 1001,0.,-1.,0.
 1003,.1,-1.,0.
 1101,0.,-1.,1.
 1103,.1,-1.,1.
*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=POREP,GENERATE
 1, 101, 10
 3, 103, 10
 1001, 1101, 10
 1003, 1103, 10
*ELEMENT,TYPE=C3D6Ph,ELSET=BLOCK
  1, 1013,1003,1001,  13,   3,   1
  3, 1023,1013,1011,  23,  13,  11
  5, 1033,1023,1021,  33,  23,  21
  7, 1043,1033,1031,  43,  33,  31
  9, 1053,1043,1041,  53,  43,  41
 11, 1063,1053,1051,  63,  53,  51
 13, 1073,1063,1061,  73,  63,  61
 15, 1083,1073,1071,  83,  73,  71
 17, 1093,1083,1081,  93,  83,  81
 19, 1103,1093,1091, 103,  93,  91
  2, 1001,1011,1013,   1,  11,  13
  4, 1011,1021,1023,  11,  21,  23
  6, 1021,1031,1033,  21,  31,  33
  8, 1031,1041,1043,  31,  41,  43
 10, 1041,1051,1053,  41,  51,  53
 12, 1051,1061,1063,  51,  61,  63
 14, 1061,1071,1073,  61,  71,  73
 16, 1071,1081,1083,  71,  81,  83
 18, 1081,1091,1093,  81,  91,  93
 20, 1091,1101,1103,  91, 101, 103
*SOLID SECTION,ELSET=BLOCK,MATERIAL=CORE
*MATERIAL,NAME=CORE
*ELASTIC
 50.,0.
*DENSITY
 .1, 
*POROUS BULK MODULI
 ,2.E6
*PERMEABILITY,SPECIFIC=10.
 3.7E-4, 
*SORPTION
 -100.,.04
 -10.,.05
 -4.5,.1
 -3.5,.18
 -2.,.45
 -1.,.91
 0.,1.
*SORPTION,TYPE=EXSORPTION
 -100.,.09
 -10.,.1
 -8.,.11
 -6.,.18
 -4.5,.33
 -3.,.79
 -2.,.91
 0.,1.
*INITIAL CONDITIONS,TYPE=SATURATION
 ALLN,.05
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
 POREP, -22.0, 1.0, -12.0, 0.0
*INITIAL CONDITIONS,TYPE=RATIO
 ALLN,5.
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
 BLOCK,-1.1,1.,-2.016666667,0.,0.,0.
*EQUATION
 2, 
 3,8,1.,1,8,-1.
 2, 
 1003,8,1.,1,8,-1.
 2, 
 1001,8,1.,1,8,-1.
*BOUNDARY
 ALLN,1
 ALLN,2
 BOTB,3
 BOTF,3
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=1
*GEOSTATIC
 1.E-6,1.E-6
*BOUNDARY
 1,8,,-12.
*DLOAD
 BLOCK,GRAV,10.,0.,0.,-1.
*NODE PRINT,FREQUENCY=5,NSET=LHS 
 U,RF,POR,RVT
*EL PRINT,FREQUENCY=5,POSITION=AVERAGED AT NODES
 S, 
 E, 
 SAT,POR,VOIDR
*OUTPUT, FIELD, frequency=10
*ELEMENT OUTPUT
 S,POR,SAT
*NODE OUTPUT
 U,POR
*END STEP
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=20.
 1.,1000000.
*BOUNDARY
 1,8,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT 
,1.,
*CONTROLS,PARAMETERS=FIELD,FIELD=PORE FLUID PRESSURE
,1.,
*END STEP