*HEADING
 UPRESS FOR MASS DIFFUSION ELEMENTS  [DC3D8]
 KEYWORD-PARAMETER COMBINATIONS TESTED: KAPPA,DEP; SOLUBILITY,DEP
*rest,write
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1, 0., 0., 0.
 2, 7., 0., 0.
 3, 7., 0., -7.
 4, 0., 0., -7.
 5, 0., 7., 0.
 6, 7., 7., 0.
 7, 7., 7., -7.
 8, 0., 7., -7.
 9, 3.5,  0., 0.
 10, 7.,  0., -3.5
 11, 3.5, 0., -7.
 12, 0.,  0., -3.5
 17, 0., 3.5, 0.
 18, 7., 3.5, 0.
 19, 7., 3.5, -7.
 20, 0., 3.5, -7.
 13, 3.5, 7, 0.
 14, 7.,  7, -3.5
 15, 3.5, 7, -7.
 16, 0.,  7, -3.5
*NSET,NSET=FIX1
 1,2,3,4
*NSET,NSET=FIX2
 5,6,7,8
*NSET,NSET=FIX3
 1,2,6,5
*NSET,NSET=FIX4
 2,3,7,6
*NSET,NSET=FIX5
 3,4,7,8
*NSET,NSET=FIX6
 1,4,5,8
*ELEMENT,TYPE=DC3D20, ELSET=EALL
 1, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
16,17,18,19,20
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*KAPPA,TYPE=PRESS,DEP=1
5.e-3,,,0.5
2.e-2,,,2.0
*SOLUBILITY,DEP=1
0.5,,0.5
1.5,,1.5
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=1
ALL,1.
*PHYSICAL CONST, ABSOLUTEZ=-273
*STEP
*MASS DIFFUSION,STEADY STATE
1.0,1.0
*BOUNDARY
6,11,11,100.
*PRESSURE STRESS,USER
ALL, 
*NODE PRINT
NNC11, 
*EL PRINT
CONC, 
ESOL, 
ISOL, 
SOL, 
*NODE FILE
NNC, 
*EL FILE
CONC, 
ESOL, 
ISOL, 
SOL, 
*END STEP