*HEADING
 Basic verification for COUPLING option
 Test a uniform distributed load applied through the
 COUPLING option to a surface with underlying elements C3D8
*NODE, NSET=FRONT
 1,   0.0, 0.0, 0.0
 2,   1.0, 0.0, 0.0
 3,   3.0, 0.0, 0.0
 4,   6.0, 0.0, 0.0
 5,   0.0, 1.0, 0.0
 6,   1.0, 1.0, 0.0
 7,   3.0, 1.0, 0.0
 8,   6.0, 1.0, 0.0
 9,   0.0, 3.0, 0.0
10,   1.0, 3.0, 0.0
11,   3.0, 3.0, 0.0
12,   6.0, 3.0, 0.0
13,   0.0, 6.0, 0.0
14,   1.0, 6.0, 0.0
15,   3.0, 6.0, 0.0
16,   6.0, 6.0, 0.0
17,   0.0,10.0, 0.0
18,   1.0,10.0, 0.0
19,   3.0,10.0, 0.0
20,   6.0,10.0, 0.0
*NCOPY, CHANGE NUMBER=100, OLD SET=FRONT, SHIFT
 0.0, 0.0, 1.0

*NCOPY, CHANGE NUMBER=200, OLD SET=FRONT, SHIFT
 0.0, 0.0, 3.0

*NCOPY, CHANGE NUMBER=300, OLD SET=FRONT, SHIFT
 0.0, 0.0, 6.0

*NCOPY, CHANGE NUMBER=400, OLD SET=FRONT, SHIFT
 0.0, 0.0, 10.0

*ELEMENT, TYPE=C3D8, ELSET=EALL
 1, 1, 2, 6, 5, 101, 102, 106, 105
*ELGEN, ELSET=EALL
 1, 3, 1, 1, 4, 4, 3, 4, 100, 12
*SOLID SECTION, MATERIAL=ELASTIC, ELSET=EALL
*MATERIAL, NAME=ELASTIC
*ELASTIC
 30.0e6,0.3
*ELSET, ELSET=FACEELEM, GEN
 1,46,3
*SURFACE, NAME=LSURFACE
 FACEELEM,
*NODE, NSET=NREF
 9999, -1, 5.0, 5.0
*COUPLING, CONSTRAINT NAME=LOAD_DISTRIBUTOR, SURFACE=LSURFACE, REF NODE=NREF
*DISTRIBUTING, WEIGHT=UNIFORM
*NSET,NSET=A
 4,8,12,16,20
 104,108,112,116,120
 204,208,212,216,220
 304,308,312,316,320
 404,408,412,416,420
*BOUNDARY
 A, 1, 1
 4, 2, 3
104,2, 2
204,2, 2
304,2, 2
404,2, 2
8, 3, 3
12,3, 3
16,3, 3
20,3, 3
*STEP
 Apply uniform load
*STATIC
 ,1.0
*CLOAD
 NREF, 1, 2000.0
*EL PRINT, FREQ=10000
 S
*NODE PRINT, FREQ=10000
 U, 
 RF,
*NODE FILE, FREQ=10000
 U, 
 RF,
*OUTPUT, FIELD, FREQ=1000
*NODE OUTPUT
 U, 
 RF,
*END STEP