*Heading                                                                                            
 restary verification test for s4 elements using distributions   
*Restart, read, step=8  
*Boundary                                                                                           
 side1, 1, 6                                   
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-1
 static
*static
 0.5, 1.0
*controls, parameters=field
 1.e-3, 1.e-4
*Boundary, op=new
 side1, 1,6
*end step
**
** --------------------------------------------------
**
*step, nlgeom=yes, name=step-2
 static
*static
 1.0, 1.0
*Boundary
 side3, 2, 2, 0.025
 side3, 3, 3, 0.050
*end step
**
** --------------------------------------------------
**
*step, perturbation, name=step-3
 static 
*static
*El print, freq=999                                                                                 
 s,                                                                                                 
 sf,
 ssavg,
 se,
 sth,      
*El print, freq=999                                                                                 
 3,                                                                                             
 tshr,                                                                                              
*Node print, freq=999                                                                               
 u,                                                                                                 
*Output, field                                                                                      
*Node output                                                                                        
 u,                                                                                                 
*Element output                                                                                     
 s,                                                                                                 
 sf,
 ssavg,
** se,
 sth,
*Element output                                                                                     
 3,                                                                                              
 tshr,            
*load case, name=load1
*dload
 eall, bz, 1.e3                                                                             
*end load case
*load case, name=load2
*dload
 eall, grav, 1.e6, 0.,0.,1.                                                                                
*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