*HEADING
 Basic verification for *COUPLING
 Test tributary area of surface with underlying element B21
*PREPRINT, MODEL=YES
*NODE, NSET=A
  1,  0.0000,  0.0,  0.0
  2,  1.2500,  0.0,  0.0
  3,  2.5000,  0.0,  0.0
  4,  3.2500,  0.0,  0.0
  5,  4.0000,  0.0,  0.0
  6,  4.5000,  0.0,  0.0
  7,  5.0000,  0.0,  0.0
  8,  5.5000,  0.0,  0.0
  9,  6.0000,  0.0,  0.0
 10,  6.7500,  0.0,  0.0
 11,  7.5000,  0.0,  0.0
 12,  8.7500,  0.0,  0.0
 13, 10.0000,  0.0,  0.0
*ELEMENT, ELSET=ELEM, TYPE=B21
 1, 1, 2
*ELGEN, ELSET=ELEM
 1, 12, 1, 1
*ELEMENT, TYPE=MASS, ELSET=MASS
1000,9999
*MASS, ELSET=MASS
1.E-6,
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY
1001,9999
*ROTARY INERTIA, ELSET=ROTARY
1.E-6,1.E-6,1.E-6
*BEAM SECTION, SECT=CIRC, ELSET=ELEM, MATERIAL=ELASTIC
 0.5,
*MATERIAL, NAME=ELASTIC
*ELASTIC
 30e6, 0.3
*DENSITY
1.0,
*SURFACE, NAME=TOP
 ELEM, SPOS
*NODE
 9999, 5.0, 5.0, 0.0
*COUPLING, CONSTRAINT NAME=LOAD, REF NODE=9999, SURFACE=TOP
*DISTRIBUTING
*NSET, NSET=NREF
9999,
*NSET, NSET=QA_TEST
6,7,8,
*ELSET, ELSET=QA_TEST
6,7,
*AMPLITUDE, NAME=AMP,DEFINITION=SMOOTH STEP
0.0,0.0,.5,1.0,1.0,1.0
*STEP
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY
 1,  1, 2, 0.0
 13, 2, 2 ,0.0
*CLOAD, AMPLITUDE=AMP
 9999, 2, -1.E+3
*OUTPUT,FIELD, NUMBER=20, VARIABLE=PRESELECT
***OUTPUT,HISTORY, TIME INTERVAL=5.E-3,VARIABLE=PRESELECT
*OUTPUT,FIELD,NUMBER=2
*NODE OUTPUT,NSET=QA_TEST
 U,
*ELEMENT OUTPUT,ELSET=QA_TEST
 S, 
*END STEP