*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=C3D6PhT,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
1.3E3,0.
*DENSITY
1.5, 
*EXPANSION
 .0001, 
*CONDUCTIVITY
 0.1
*SPECIFIC HEAT
 10.0
*POROUS BULK MODULI
 ,2.E6
*PERMEABILITY,SPECIFIC=10.
4.5E-6, 
*PERMEABILITY,TYPE=SATURATION
.666666,.85
1.,1.
*SORPTION
-100.,.8
-10.,.85
-9.5,.9
-9.,.922
-8.,.947
-7.,.961
-6.,.973
-4.,.988
-2.,.999
0.,1.
*SORPTION,TYPE=EXSORPTION
-100.,.8
-10.,.85
-9.5,.9
-9.,.922
-8.,.947
-7.,.961
-6.,.973
-4.,.988
-2.,.999
0.,1.
*INITIAL CONDITIONS,TYPE=SATURATION
 ALLN,1.
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
 POREP, 0.
*INITIAL CONDITIONS,TYPE=RATIO
 ALLN,.4235
*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC
 BLOCK,0.,1.,-17.9750615,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,2
 BOTB,3
 BOTF,3
 BOTF,11
*STEP,INC=1
*GEOSTATIC
 1.E-6,1.E-6
*DLOAD
 BLOCK,GRAV,10.,0.,0.,-1.
*BOUNDARY
 ALLN,8,,0.
*OUTPUT,FIELD, frequency=10
*NODE OUTPUT
 U,POR
*ELEMENT OUTPUT
 SAT,S
*END STEP
*STEP,INC=100
*SOILS,CONSOLIDATION,UTOL=10.
 20.,50000.
*BOUNDARY,OP=NEW
 1,8,,0.
 ALLN,1,2,0.
 BOTB,3,,0.
 BOTF,3,,0.
*CONTROLS,ANALYSIS=DISCONTINUOUS
*END STEP