*HEADING
  Restart test; restart from pmcp_2drebar, and duplicate steps 4 onwards
*RESTART,READ,step=3,endstep
*FILE FORMAT
**
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
** set initial time averaged force to what it was at the beginning of
** step 4 from original analysis (pmcp_2drebar)
*CONTROLS,PARAMETERS=FIELD,FIELD=DISPL
 , , 4.277E+07
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
6003,P4,-0.25
6004,P4,-0.25
6001,P4,-0.25
6002,P4,-0.25
10001,P4,-0.25
10002,P4,-0.25
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
 S,E
*ELEMENT OUTPUT, REBAR
S,
E,
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*ENERGY FILE,F=0
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
RTEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
RTEST,
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #7: APPLY THERMAL LOAD
*STATIC
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*TEMPERATURE
NALL,100.
*EL PRINT,FREQ=999
S,THE,TEMP
*EL PRINT,REBAR,FREQ=999
S,THE,TEMP
*EL FILE,FREQ=999
S,THE,TEMP
*ENERGY FILE,FREQ=999
*END STEP
**
**
*STEP
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
32, , 
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0004,.004
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E,S
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
1000., 10000., 5
*BOUNDARY,OP=NEW
LEFT,1,2
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*END STEP
**
**
*STEP
 STEP #12: DYNAMICS
*DYNAMIC,DIRECT
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,0.2
*ENERGY FILE,FREQ=999
*END STEP