*HEADING
EC36MFDC: DFLUX FOR MASS DIFFUSION ELEMENTS  [DC3D6]
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*RESTART,WRITE
*NODE,NSET=ALL
  1, 0.0, 0.0, 0.0
  2, 7.0, 0.0, 0.0
  3, 7.0, 7.0, 0.0
  4, 0.0, 7.0, 0.0
  5, 0.0, 0.0, 7.0
  6, 7.0, 0.0, 7.0
  7, 7.0, 7.0, 7.0
  8, 0.0, 7.0, 7.0
*ELEMENT,TYPE=DC3D6,ELSET=ALLELS
 1, 1,2,4,5,6,8
 2, 2,3,4,6,7,8
*SOLID SECTION,MATERIAL=A1, ELSET=ALLELS
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*SOLUBILITY
 1., 
*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,8,7
*NSET,NSET=FIX6
4,1,5,8
*BOUNDARY
 FIX1,11
*STEP
*MASS DIFFUSION, STEADY STATE
*NODE PRINT
 NNC,RFL
*EL PRINT,FREQ=0
*DFLUX
 ALLELS, BF, .3
*OUTPUT,FIELD
*NODE OUTPUT,NSET=ALL
 NNC,RFL
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ALL
 NNC,RFL
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX2, 11
*DFLUX, OP=NEW
 1, S1, .3
 2, S1, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX1, 11
*DFLUX,OP=NEW
 1, S2, .3
 2, S2, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX5, 11
*DFLUX,OP=NEW
 1, S3, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX6, 11
*DFLUX,OP=NEW
 2, S3, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX3, 11
*DFLUX,OP=NEW
 2, S4, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY, OP=NEW
 FIX4, 11
*DFLUX,OP=NEW
 1, S5, .3
*END STEP