;
;  Doppler Data Simulation spectrum
;
  lmax=4096
  amp1 = fltarr(lmax)
  amp2 = fltarr(lmax)
  amp3 = fltarr(lmax)
  SIM_wavenumber=findgen(lmax)
;
; Read in old spectral coefficients
;
; openr,1,'path/InputSpectrum20150511.text'
  openr,1,file
  for i=0,lmax-1 do begin $
    readf,1,l,a1,a2,a3 &$
    amp1(l-1)=a1 &$
    amp2(l-1)=a2 &$
    amp3(l-1)=a3 &$
  endfor
  close,1
  
  SIM_velocitySpectrum=fltarr(lmax)
  KE=fltarr(lmax)
  el=fltarr(lmax)
  for l=0,lmax-1 do begin $
    el[l]=float(l) &$
    KE[l]=(2.*el[l] + 1. )*(amp1(l)*amp1(l) + el[l]*(el[l] + 1.)*(amp2(l)*amp2(l) + amp3(l)*amp3(l)))/2. &$
    SIM_velocitySpectrum[l]=sqrt(el[l]*KE[l]) &$
  endfor

  save,/ver,file='hmi_spec.save',SIM_velocitySpectrum,KE,el