*HEADING
 Basic verification for COUPLING option
 Test a uniform distributed load applied through the
 COUPLING option to a surface with underlying elements C3D20
*NODE, NSET=FRONT
  1,  0.0, 0.0
  2,  0.5, 0.0
  3,  1.0, 0.0
  4,  2.0, 0.0
  5,  3.0, 0.0
  6,  4.5, 0.0
  7,  6.0, 0.0
  8,  0.0, 0.5
  9,  0.5, 0.5
 10,  1.0, 0.5
 11,  2.0, 0.5
 12,  3.0, 0.5
 13,  4.5, 0.5
 14,  6.0, 0.5
 15,  0.0, 1.0
 16,  0.5, 1.0
 17,  1.0, 1.0
 18,  2.0, 1.0
 19,  3.0, 1.0
 20,  4.5, 1.0
 21,  6.0, 1.0
 22,  0.0, 2.0
 23,  0.5, 2.0
 24,  1.0, 2.0
 25,  2.0, 2.0
 26,  3.0, 2.0
 27,  4.5, 2.0
 28,  6.0, 2.0
 29,  0.0, 3.0
 30,  0.5, 3.0
 31,  1.0, 3.0
 32,  2.0, 3.0
 33,  3.0, 3.0
 34,  4.5, 3.0
 35,  6.0, 3.0
 36,  0.0, 4.5
 37,  0.5, 4.5
 38,  1.0, 4.5
 39,  2.0, 4.5
 40,  3.0, 4.5
 41,  4.5, 4.5
 42,  6.0, 4.5
 43,  0.0, 6.0
 44,  0.5, 6.0
 45,  1.0, 6.0
 46,  2.0, 6.0
 47,  3.0, 6.0
 48,  4.5, 6.0
 49,  6.0, 6.0
 50,  0.0, 8.0
 51,  0.5, 8.0
 52,  1.0, 8.0
 53,  2.0, 8.0
 54,  3.0, 8.0
 55,  4.5, 8.0
 56,  6.0, 8.0
 57,  0.0, 10.0
 58,  0.5, 10.0
 59,  1.0, 10.0
 60,  2.0, 10.0
 61,  3.0, 10.0
 62,  4.5, 10.0
 63,  6.0, 10.0
*NCOPY, CHANGE NUMBER=100, OLD SET=FRONT, SHIFT
 0.0, 0.0, 0.5

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

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

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

*NCOPY, CHANGE NUMBER=500, OLD SET=FRONT, SHIFT
 0.0, 0.0, 4.5

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

*NCOPY, CHANGE NUMBER=700, OLD SET=FRONT, SHIFT
 0.0, 0.0, 8.0

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

*NODE, NSET=NREF
 9999, 10, 5.0, 5.0
*ELEMENT, TYPE=C3D27R, ELSET=EALL
 1, 1, 3, 17, 15, 201, 203, 217, 215, 2, 10, 16, 8, 
    202, 210, 216, 208, 101, 103, 117, 115,
    109, 9, 209, 102, 110, 116, 108
*ELGEN, ELSET=EALL 
 1, 3, 2, 1, 4, 14, 3, 4, 200, 12
*SOLID SECTION, MATERIAL=ELASTIC, ELSET=EALL
*MATERIAL, NAME=ELASTIC
*ELASTIC
 30.0e6,0.3
*ELSET, ELSET=FACE, GEN
 3, 48, 3
*SURFACE, NAME=LSURFACE, TRIM=YES
 FACE,
*COUPLING, CONSTRAINT NAME=LOAD_DISTRIBUTOR, SURFACE=LSURFACE, REF NODE=NREF
*DISTRIBUTING, WEIGHT=UNIFORM
*NSET, NSET=A
  1, 101, 201, 301, 401, 501, 601, 701, 801
  8, 108, 208, 308, 408, 508, 608, 708, 808
 15, 115, 215, 315, 415, 515, 615, 715, 815
 22, 122, 222, 322, 422, 522, 622, 722, 822
 29, 129, 229, 329, 429, 529, 629, 729, 829
 36, 136, 236, 336, 436, 536, 636, 736, 836
 43, 143, 243, 343, 443, 543, 643, 743, 843
 50, 150, 250, 350, 450, 550, 650, 750, 850
 57, 157, 257, 357, 457, 557, 657, 757, 857
*NSET, NSET=B, GEN
 101, 801, 100
*NSET, NSET=C, GEN
 8, 57, 7
*BOUNDARY
 A, 1, 1
 1, 2, 3
 B, 2, 2
 C, 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