*HEADING
 Test of multistep analysis with mass scaling
*RESTART,WRITE	
*NODE,NSET=NALL
  1, 0.0, 0.
  2, 2.0, 0.
  3, 2.0, 1.
  4, 0.0, 1.
 11, 0.0, 3.
 12, 2.0, 3.
 13, 2.0, 4.
 14, 0.0, 4.
 21, 0.0, 6.
 22, 2.0, 6.
 23, 2.0, 7.
 24, 0.0, 7.
101, 4.0, 0.
102, 6.0, 0.
103, 6.0, 1.
104, 4.0, 1.
111, 4.0, 3.
112, 6.0, 3.
113, 6.0, 4.
114, 4.0, 4.
121, 4.0, 6.
122, 6.0, 6.
123, 6.0, 7.
124, 4.0, 7.
*NSET,NSET=NOUT
1,11,21,101,111,121
*NSET,NSET=LEFT
1,4,11,14,21,24, 101,104,111,114,121,124
*NSET,NSET=BOT
1,2,11,12,21,22, 101,102,111,112,121,122
*NSET,NSET=RIGHT
2,3,12,13,22,23, 102,103,112,113,122,123 
*NSET,NSET=QA_TEST
RIGHT,
*ELEMENT,TYPE=M3D4R,ELSET=REFER
1,  1,2,3,4
11, 11,12,13,14
21, 21,22,23,24
*ELEMENT,TYPE=M3D4R,ELSET=TEST1
101,  101,102,103,104
*ELEMENT,TYPE=M3D4R,ELSET=TEST2
111,  111,112,113,114
*ELEMENT,TYPE=M3D4R,ELSET=TEST3
121,  121,122,123,124
*ELSET,ELSET=EOUT
REFER,TEST1,TEST2,TEST3
*ELSET,ELSET=QA_TEST
EOUT,
*MEMBRANE SECTION,MATERIAL=ELAS,ELSET=REFER, POISSON=0.25
.1, 
*MEMBRANE SECTION,MATERIAL=ELAS1,ELSET=TEST1, POISSON=0.25
.1, 
*MEMBRANE SECTION,MATERIAL=ELAS2,ELSET=TEST2, POISSON=0.25
.1, 
*MEMBRANE SECTION,MATERIAL=ELAS3,ELSET=TEST3, POISSON=0.25
.1, 
*MATERIAL,NAME=ELAS
*ELASTIC
1.E2, .25
*DENSITY
10.E-5, 
*MATERIAL,NAME=ELAS1
*ELASTIC
1.E2, .25
*DENSITY
5.E-5, 
*MATERIAL,NAME=ELAS2
*ELASTIC
1.E2, .25
*DENSITY
2.E-5, 
*MATERIAL,NAME=ELAS3
*ELASTIC
1.E2, .25
*DENSITY
1.E-5, 
*boundary
nall,3
BOT,2
LEFT,1
*ELSET,ELSET=EALL
REFER,TEST1,TEST2,TEST3
*ELSET,ELSET=SCALE
TEST1,TEST2,TEST3
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.,0.,  .01, 1.0
*STEP
 Step 1: Test of local mass scaling
*DYNAMIC,EXPLICIT
,.01
*FIXED MASS SCALING,FACTOR=2., ELSET=TEST1
*FIXED MASS SCALING,FACTOR=5., ELSET=TEST2
*FIXED MASS SCALING,FACTOR=10.,ELSET=TEST3
*RESTART,WRITE,NUMBER INTERVAL=10
*BOUNDARY,AMP=SMOOTH
RIGHT,1,1,.1
*************************
** RF not expected to be noisy
** EDT not expected to be noisy
********************************
*OUTPUT,FIELD, NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
RF,
*ELEMENT OUTPUT,ELSET=QA_TEST
EDT,
*OUTPUT,HISTORY,TIME INTERVAL=0.01
*ENERGY OUTPUT, VAR=ALL
**********************
***OUTPUT,HISTORY,TIME INTERVAL=0.0005
***NODE OUTPUT,NSET=NALL
**RF,U
***ELEMENT OUTPUT,ELSET=EOUT
**EDT,EMSF
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,
*end step
*STEP
 Step 2: Test of propagation of mass matrix
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,15.
*end step
*STEP
 Step 3: Removing all fixed mass scaling
         Apply new variable mass scaling
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,-10.
*FIXED MASS SCALING
*VARIABLE MASS SCALING,DT=.001,TYPE=SET EQUAL DT,FREQUENCY=1
*end step
*STEP
 Step 4: Test of propagation of variable mass scaling
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,-20.
*end step
*STEP
 Step 5: Remove all variable mass scaling defintions
         Respecify fixed mass scaling
*DYNAMIC,EXPLICIT
,.01
*VARIABLE MASS SCALING
*FIXED MASS SCALING,FACTOR=2.0
*FIXED MASS SCALING,ELSET=TEST1,FACTOR= 4.0
*FIXED MASS SCALING,ELSET=TEST2,FACTOR=10.0
*FIXED MASS SCALING,ELSET=TEST3,FACTOR=20.0
*BOUNDARY,TYPE=VEL
RIGHT,1,1,10.
*END STEP
*STEP
 Step 6: Test of propagation of mass matrix
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,15.
*end step
*STEP
 Step 7: Set fixed mass scaling
*DYNAMIC,EXPLICIT
,.01
*FIXED MASS SCALING,ELSET=REFER,FACTOR=2.0
*FIXED MASS SCALING,ELSET=TEST1,FACTOR= 4.0
*FIXED MASS SCALING,ELSET=TEST2,FACTOR=10.0
*FIXED MASS SCALING,ELSET=TEST3,FACTOR=20.0
*BOUNDARY,TYPE=VEL
RIGHT,1,1,20.
*end step
*STEP
 Step 8: Test propagation of local and global fixed m.s.
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,25.
*end step
*STEP
 Step 9: Set variable mass scaling
*DYNAMIC,EXPLICIT
,.01
*VARIABLE MASS SCALING,DT=.001,TYPE=SET EQUAL DT,FREQ=2
*VARIABLE MASS SCALING,DT=.001,TYPE=SET EQUAL DT,ELSET=SCALE,FREQ=1
*BOUNDARY,TYPE=VEL
RIGHT,1,1,30.
*end step
*STEP
 Step 10: Test propagation of local and global variable m.s.
*DYNAMIC,EXPLICIT
,.01
*BOUNDARY,TYPE=VEL
RIGHT,1,1,35.
*end step