*Heading                                                                                            
 Multistep verification test for cgax4h elements using distributions                                
*Node                                                                                               
        1,  0.1000000E-01,  0.0000000E+00,  0.0000000E+00
        2,  0.3333333E+00,  0.0000000E+00,  0.0000000E+00
        3,  0.6666667E+00,  0.0000000E+00,  0.0000000E+00
        4,  0.1000000E+01,  0.0000000E+00,  0.0000000E+00
        5,  0.1000000E-01,  0.3333333E+00,  0.0000000E+00
        6,  0.3333333E+00,  0.3333333E+00,  0.0000000E+00
        7,  0.6666667E+00,  0.3333333E+00,  0.0000000E+00
        8,  0.1000000E+01,  0.3333333E+00,  0.0000000E+00
        9,  0.1000000E-01,  0.6666667E+00,  0.0000000E+00
       10,  0.3333333E+00,  0.6666667E+00,  0.0000000E+00
       11,  0.6666667E+00,  0.6666667E+00,  0.0000000E+00
       12,  0.1000000E+01,  0.6666667E+00,  0.0000000E+00
       13,  0.1000000E-01,  0.1000000E+01,  0.0000000E+00
       14,  0.3333333E+00,  0.1000000E+01,  0.0000000E+00
       15,  0.6666667E+00,  0.1000000E+01,  0.0000000E+00
       16,  0.1000000E+01,  0.1000000E+01,  0.0000000E+00
*Nset, nset=side1, generate                                                                         
        1,       4,       1,
*Nset, nset=side3, generate                                                                         
       13,      16,       1,
*Elset, Elset=eall, generate                                                                        
        1,       9,       1,
**------------------------------------------------------------------------------------------        
*Material, name=mat1                                                                                
*Elastic, type=engineering constants                                                                
  0.300E+07,  0.400E+07,  0.500E+07,  0.100E+00,  0.200E+00,  0.300E+00,  0.300E+06,  0.400E+06
  0.500E+06,
*Density                                                                                            
  0.750E-03,
**------------------------------------------------------------------------------------------        
*Element, type=cgax4h, elset=eall                                                                   
        1,       1,       2,       6,       5,
        2,       2,       3,       7,       6,
        3,       3,       4,       8,       7,
        4,       5,       6,      10,       9,
        5,       6,       7,      11,      10,
        6,       7,       8,      12,      11,
        7,       9,      10,      14,      13,
        8,      10,      11,      15,      14,
        9,      11,      12,      16,      15,
*Distribution table, name=OrientationTable
 coord3d, coord3d                                             
*Distribution, name=distOrient, location=element, table=orientationTable                                  
         ,   0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00,
        1,   0.1000000E+01,   0.0000000E+00,   0.0000000E+00,   0.0000000E+00,   0.1000000E+01,   0.0000000E+00,
        2,   0.9998477E+00,   0.1745241E-01,   0.0000000E+00,  -0.1745241E-01,   0.9998477E+00,   0.0000000E+00,
        3,   0.9993908E+00,   0.3489950E-01,   0.0000000E+00,  -0.3489950E-01,   0.9993908E+00,   0.0000000E+00,
        4,   0.9986295E+00,   0.5233596E-01,   0.0000000E+00,  -0.5233596E-01,   0.9986295E+00,   0.0000000E+00,
        5,   0.9975641E+00,   0.6975647E-01,   0.0000000E+00,  -0.6975647E-01,   0.9975641E+00,   0.0000000E+00,
        6,   0.9961947E+00,   0.8715574E-01,   0.0000000E+00,  -0.8715574E-01,   0.9961947E+00,   0.0000000E+00,
        7,   0.9945219E+00,   0.1045285E+00,   0.0000000E+00,  -0.1045285E+00,   0.9945219E+00,   0.0000000E+00,
        8,   0.9925462E+00,   0.1218693E+00,   0.0000000E+00,  -0.1218693E+00,   0.9925462E+00,   0.0000000E+00,
        9,   0.9902681E+00,   0.1391731E+00,   0.0000000E+00,  -0.1391731E+00,   0.9902681E+00,   0.0000000E+00,
*Distribution Table, name=AngleTable
 angle
*Distribution, name=distAngle, table=AngleTable, location=element
  , 0.0
 1, 1.0
 2, -1.0
 3, 2.0
 4, -2.0
 5, 3.0
 6, -3.0
 7, 4.0
 8, -4.0
 9, 5.0
*Solid Section, material=mat1, elset=eall, orientation=ori                                                           
1.,                                                       
*Orientation, name=ori
 distOrient
 3, distAngle
*Boundary                                                                                           
 side1, 1, 2                                                                                        
 side1, 5, 5                                                                                        
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-1
 static
*static
 .25, 1.0
*Boundary
 side3, 1, 1, 0.01
 side3, 5, 5, 0.01
*El print, freq=999                                                                                 
 s,                                                                                                 
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field, op=new                                                                                            
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-2
 static 
*static
*El print, freq=999                                                                                 
 s,      
*Node print, freq=0                                                                          
*Output, field, op=new                                                                                            
*Element output                                                                                     
 s,                                                                                                 
*load case, name=load1
*dload
 eall, br, 1.e3                                                                             
*end load case
*load case, name=load2
*dload
 eall, grav, 1.e6, 0.,1.,0.                                                                                
*end load case
*load case, name=load3
*dload
 eall, cent, 100.,0.,0.,0.,1.,0.                                                                                 
*end load case
*load case, name=load4
*dload
 eall, centrif, 1.e6,0.,0.,0.,1.,0.                                                                                 
*end load case
*end step