restore,/ver,'spectrumxyz.sav'
file='InputSpectrum20150511.text'
@hathaway_spectrum
Rsun=696
ell=k*Rsun
plot,el,SIM_VELOCITYSPECTRUM,/xlog,/ylog,xr=[el[1],el[4094]],ytit='!5Velocity (m/s)',xtit='!5Spherical Harmonic Degree',$
chars=1.5,thick=3
oplot,ell,sqrt((pvxm+pvym+pvzm)*k)*1.e4,lin=2,thick=3
restore,/ver,'v_k.sav'
oplot,k*6.96e10,v_k*1.e-2,lin=4,thick=4