*HEADING
EARTH DAM - STEADY STATE FREE SURFACE SEEPAGE      
*** UNITS: M, KG, SEC, NEWTON
*NODE,NSET=ALLN
1,0.,0.
39,4.8768,0.
601,1.8288,1.8288
639,3.048,1.8288
*NGEN,NSET=BOT
1,39,1
*NGEN,NSET=TOP
601,639,1
*NFILL,NSET=ALLN
BOT,TOP,12,50
*NSET,NSET=PORN0,GENERATE
1,39,2
*NSET,NSET=PORN1,GENERATE
101,139,2
*NSET,NSET=PORN2,GENERATE
201,239,2
*NSET,NSET=PORN3,GENERATE
301,339,2
*NSET,NSET=PORN4,GENERATE
401,439,2
*NSET,NSET=PORN5,GENERATE
501,539,2
*NSET,NSET=PORN6,GENERATE
601,639,2
*NSET,NSET=PORN
 PORN0, PORN1, PORN2, PORN3
*NSET,NSET=OUTN,GENERATE
1,601,100
21,621,100
*ELEMENT,TYPE=CPE4P,ELSET=DAM
1,1,3,103,101
*ELGEN,ELSET=DAM
1,19,2,1,6,100,20
*ELSET,ELSET=FSIDE,GENERATE
19,119,20
*ELSET,ELSET=FBOT,GENERATE
16,19
*ELSET,ELSET=OUTE,GENERATE
11,111,20
*SOLID SECTION,ELSET=DAM,MATERIAL=FILL
*MATERIAL,NAME=FILL
*ELASTIC
1000., 
*DENSITY
2000., 
*PERMEABILITY,SPECIFIC=10000.
2.1167E-4, 
*SORPTION
-100000.,.04
-10000.,.05
0.,1.
*** DUMMY CONTACT DEFINITION
*ELEMENT,TYPE=CPE4P,ELSET=DAM
6,3011,13,113,3111
26,3111,113,213,3211
46,3211,213,3313,3311
47,213,215,3315,3313
48,215,217,3317,3315
49,217,219,3319,3317
50,219,221,3321,3319
51,221,223,3323,3321
52,223,225,3325,3323
53,225,227,3327,3325
54,227,3229,3329,3327
34,127,3129,3229,227
14,27,3029,3129,127
*NSET, NSET=CONTACT, GEN
11,311,100
311,329,2
329,29,-100
*NCOPY, change number=3000,oldset=CONTACT,shift,new set=CONTACT_NEW
0,0,0
0,0,0,0,0,1
*NSET,NSET=ALLN
ALLN,CONTACT_NEW
*NSET,NSET=PORN
PORN,CONTACT_NEW
*ELSET,ELSET=sla1
6,26,46
*ELSET,ELSET=sla2,gen
46,54,1
*ELSET,ELSET=sla3
14,34,54
*ELSET,ELSET=mas1
5,25,45
*ELSET,ELSET=mas2,gen
66,74,1
*ELSET,ELSET=mas3
15,35,55
*surface,NAME=MAIN
mas1,S2
mas2,S1
mas3,S4
*surface, NAME=SECONDARY
sla1,S4
sla2,S3
sla3,S2
*CONTACT PAIR, INTERACTION=DUMMY, TIED, ADJUST=CONTACT_NEW
SECONDARY,MAIN
*SURFACE INTERACTION,NAME=DUMMY
1.0,
**
**
*MAP SOLUTION
*INITIAL CONDITIONS,TYPE=SATURATION
ALLN,1.
*INITIAL CONDITIONS,TYPE=PORE PRESSURE
PORN, 12192.0, 0.0, 0.0, 1.2192
PORN4,0.
PORN5,0.
PORN6,0.
*INITIAL CONDITIONS,TYPE=RATIO
ALLN,1.
*BOUNDARY
ALLN,1
ALLN,2
*RESTART,WRITE,FREQUENCY=10
*STEP,INC=5
*SOILS
.2,1.
*DLOAD
DAM,GRAV,10.,0.,-1.,0.
*BOUNDARY
1,8,,12192.
101,8,,9144.
201,8,,6096.
301,8,,3048.
401,8,,0.
*FLOW
FSIDE,Q2D,0.1
FBOT,Q1D,0.1
*CONTROLS,ANALYSIS=DISCONTINUOUS
*NODE PRINT,FREQUENCY=10,NSET=OUTN
POR, 
*EL PRINT,FREQUENCY=10,ELSET=OUTE
SAT,POR
*NODE FILE,FREQUENCY=10,NSET=OUTN
POR, 
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
POR,
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
SAT,POR
*END STEP