nt=n_elements(t)
it2=nt-1
plot_oo,k[1:*],grav1[1:*,it2]*k^2.0/total(spec1[1:*,it2])^2,yr=[.01,500]
END