*HEADING
ECA4MFDC: DFLUX FOR MASS DIFFUSION ELEMENTS  [DCAX4]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1, 1.
 2, 8.
 3, 8.,7.
 4, 1.,7.
*NSET,NSET=FIX1
 1,2
*NSET,NSET=FIX2
 2,3
*NSET,NSET=FIX3
 3,4
*NSET,NSET=FIX4
 1,4
*ELEMENT,TYPE=DCAX4, ELSET=EALL
 1,1,2,3,4
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*SOLUBILITY
 1., 
*STEP
*MASS DIFFUSION, STEADY STATE
*NODE PRINT
 NNC,RFL
*NODE FILE
 NNC,RFL
*EL PRINT,FREQ=0
*EL FILE
*DFLUX,OP=NEW
 1, BF, .3
*BOUNDARY, OP=NEW
 FIX1, 11
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX3,11
*DFLUX,OP=NEW
 1, S1, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX4,11
*DFLUX,OP=NEW
 1, S2, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX1,11
*DFLUX,OP=NEW
 1, S3, .3
*END STEP
*STEP
*MASS DIFFUSION, STEADY STATE
*BOUNDARY,OP=NEW
 FIX2,11
*DFLUX,OP=NEW
 1, S4, .3
*END STEP