*HEADING
Testing IncState data transfer between steps
Expanded domain in the second step.
In the first step, edge-edge overclosure is maintained by appropriate
boundary conditions. In the second step, bc are removed
**
*NODE,NSET=N_EBODY1
      1,           100.,           -30.,             0.1
      2,           100.,           -30.,            -30.
      3,           100.,            30.,            -30.
      4,           100.,            30.,             0.1
**
*ELEMENT,TYPE=S4R,ELSET=E_EBODY1
  1,   1,   2,   3,  4
*SHELL SECTION,ELSET=E_EBODY1,MATERIAL=EMAT1
 5.0,
*MATERIAL, NAME=EMAT1
*DENSITY
 7e-09,
*ELASTIC
 2e+05, .3
**
*NODE,NSET=N_EBODY2
      5,           70.,            0.,            0.
      6,           70.,            0.,           30.
      7,          130.,            0.,           30.
      8,          130.,            0.,            0.
**
*ELEMENT,TYPE=S4R,ELSET=E_EBODY2
  2,   5,   6,   7,  8
*SHELL SECTION,ELSET=E_EBODY2,MATERIAL=EMAT2
 5.0,
*MATERIAL, NAME=EMAT2
*DENSITY
 7e-09,
*ELASTIC
 2e+05, .3
**
*ELSET,ELSET=E_EALL
E_EBODY1,E_EBODY2 
** 
*SURFACE,TYPE=ELEMENT, NAME=EAUTO
 E_EALL,
**
*NODE,NSET=N_FBODY1
      101,         -25.,          25.,           0.1
      102,         -25.,         -25.,           0.1
      103,          25.,         -25.,           0.1
      104,          25.,          25.,           0.1
**
*ELEMENT,TYPE=S4R,ELSET=E_FBODY1
  101,   101,   102,   103,  104
*SHELL SECTION,ELSET=E_FBODY1,MATERIAL=FMAT1
 5.0,
*MATERIAL, NAME=FMAT1
*DENSITY
 7e-09,
*ELASTIC
 2e+05, .3
**
*NODE,NSET=N_FBODY2
      105,         -30.,          30.,           5.
      106,         -30.,         -30.,           5.
      107,          30.,         -30.,           5.
      108,          30.,          30.,           5.
**
*ELEMENT,TYPE=S4R,ELSET=E_FBODY2
  102,   105,   106,   107,  108
*SHELL SECTION,ELSET=E_FBODY2,MATERIAL=FMAT2
 5.0,
*MATERIAL, NAME=FMAT2
*DENSITY
 7e-09,
*ELASTIC
 2e+05, .3
**
*ELSET,ELSET=E_FALL
E_FBODY1,E_FBODY2
**
*SURFACE,TYPE=ELEMENT, NAME=FAUTO
 E_FALL,
**
*ELSET,ELSET=E_ALL
E_EBODY1,E_EBODY2,E_FBODY1,E_FBODY2
**
*SURFACE,TYPE=ELEMENT, NAME=ALL_AUTO
 E_ALL, 
**
*NSET, NSET=QA_TEST_ALL_NODES
N_FBODY1, N_FBODY2, N_EBODY1, N_EBODY2
**
*STEP
*DYNAMIC, EXPLICIT
 , 5.0e-4
** Only edge case is processed
*FILE OUTPUT, NUM=1
*NODE FILE
 U,
*CONTACT
*CONTACT INCLUSIONS
 EAUTO,
**
**Output requests for post-processing
*OUTPUT, HISTORY, VARIABLE=PRESELECT, TIME INTERVAL=10
**
*OUTPUT, FIELD, VARIABLE=PRESELECT
**
*BOUNDARY, TYPE=velocity
N_EBODY1, 3, 3, 1000
*BOUNDARY, TYPE=velocity
N_EBODY2, 3, 3, -1000
**
**Output requests for QA
*OUTPUT, FIELD, NUMBER INTERVAL=10
*NODE OUTPUT, NSET=QA_TEST_ALL_NODES
 U,
**
*END STEP
**
*STEP
*DYNAMIC, EXPLICIT
 , 5.0e-4
**
*CONTACT, OP = NEW
** Both pairs are processed in this step
*CONTACT INCLUSIONS,  ALL EXTERIOR
**
**Output requests for post-processing
*OUTPUT, HISTORY, VARIABLE=PRESELECT, TIME INTERVAL=10
*OUTPUT, FIELD, VARIABLE=PRESELECT
**
**BOUNDARY, OP=NEW
*BOUNDARY, TYPE=velocity
N_EBODY2, 3, 3, 0
**
**Output requests for QA
*OUTPUT, FIELD, NUMBER INTERVAL=10
*NODE OUTPUT, NSET=QA_TEST_ALL_NODES
 U,
**
*END STEP