*HEADING
 BASIC VERIFICATION OF PBILINEAR MPC
*NODE
 1, 1.
 5, 4.
21, 1., 5.
25, 4., 5.
*NGEN,NSET=BOT
1,5
*NGEN,NSET=TOP
21,25
*NFILL,NSET=N1
BOT,TOP,4,5
*NCOPY,CHANGE NUMBER=50,OLD SET=N1,SHIFT,NEW SET=N2
0.0,0.0,-0.25
          ,           ,           ,           ,           ,  
*NCOPY,CHANGE NUMBER=100,OLD SET=N1,SHIFT,NEW SET=N3
0.0,0.0,-0.5

*NCOPY,CHANGE NUMBER=150,OLD SET=N1,SHIFT,NEW SET=N4
0.0,0.0,-0.75

*NCOPY,CHANGE NUMBER=200,OLD SET=N1,SHIFT,NEW SET=N5
0.0,0.0,-1.0

*NSET,NSET=NALL
N1,N2,N3,N4,N5
*NSET,NSET=CORNERS1
1,5,201,205
*NSET,NSET=CORNERS2
21,23,25,121,123,125,221,223,225
*NSET,NSET=CORNERS3
1,5,11,13,15,21,23,25
*NSET,NSET=CORNERS4
5,205,15,115,215,25,125,225
*NSET,NSET=CORNERS5
201,205,211,213,215,221,223,225
*NSET,NSET=CORNERS6
1,201,11,111,211,21,121,221
*MPC
P BILINEAR,113,11,15,215,211
P LINEAR,13,11,15
P LINEAR,115,15,215
P LINEAR,213,211,215
P LINEAR,111,11,211
C BIQUAD,63,11,15,115,111,13,65,113,61
C BIQUAD,114,13,15,215,213,14,115,214,113
C BIQUAD,163,111,115,215,211,113,165,213,161
C BIQUAD,112,11,13,213,211,12,113,212,111
QUADRATIC,65,15,115,215
QUADRATIC,165,15,115,215
QUADRATIC,61,11,111,211
QUADRATIC,161,11,111,211
QUADRATIC,12,11,13,15
QUADRATIC,14,11,13,15
QUADRATIC,212,211,213,215
QUADRATIC,214,211,213,215
*ELEMENT,TYPE=C3D20P,ELSET=EALL
1, 1,5,205,201, 11,15,215,211, 3,105,203,101, 13,115,213,
   111, 6,10,210,206
*ELEMENT,TYPE=C3D20P,ELSET=EALL
2, 11,13,113,111, 21,23,123,121, 12,63,112,61, 22,73,122,
   71, 16,18,118,116
*ELGEN,ELSET=EALL
2,2,2,1,2,100,2
*SOLID SECTION,MATERIAL=A1,ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 1.E8, 
*PERMEABILITY,SPECIFIC=1.0
 1.E-5, 
*INITIAL CONDITIONS,TYPE=RATIO
NALL,1.
*AMPLITUDE, NAME=RAMPP
0.,0.,  1.,1.
*STEP
*SOILS
 1. , 1.
*EL FILE
 LOADS,
*OUTPUT,FIELD
*ELEMENT OUTPUT
LOADS,
*NODE FILE
 RF, 
 POR, 
*OUTPUT,FIELD
*NODE OUTPUT
RF,
POR,
*EL PRINT
 LOADS,
*NODE PRINT
 RF, 
 POR, 
*BOUNDARY
 NALL, 1,3
 CORNERS3, 8
*DFLOW, OP=MOD, AMPLITUDE=RAMPP
 1, S5, 3.E-5
 4, S5, 3.E-5
 5, S5, 3.E-5
*END STEP
*STEP
*SOILS
 1., 1.
*BOUNDARY, OP=NEW
 NALL, 1,3
 CORNERS4, 8
*DFLOW, OP=NEW, AMPLITUDE=RAMPP
 1, S6, 3.E-5
 2, S6, 3.E-5
 4, S6, 3.E-5
*END STEP