; j is the 
restore,'spec_range.sav'
!p.multi=[0,1,2]
mu0=1e4
yr1=[1e9,1e10] & yr2=[1e9,1e10] & gam2=.5e5 & i=50 & it=0 ;(B10)
yr1=[1e2,1e4] & yr2=[1e11,1e14] & gam2=.5e2 & i=50 & it=0 ;(B1)
yr1=[1e5,2e7] & yr2=[1e4,1e5] & gam2=1e3  & i=6 & it=10 ;(M512sig1_k6_ramp1b)
;
plot_io,grav1_range[i,*]/grah1_range[i,*],yr=yr1
print,grav1_range[i,it]/grah1_range[i,it]/gam2^2
;
plot_io,stre1_range[i,*]/grav1_range[i,*],yr=yr2
print,stre1_range[i,it]/grav1_range[i,it]/gam2^2
print,stre1_range[i,it]/grav1_range[i,it]/(gam2*mu0)^2
;
!p.multi=0
END