pro groups_np,file,wave,ray ;data from HST archive ;file='b202t' root='y0ri'+file dir='./data/' wavefile=dir+root+'_c0f.fits' fluxfile=dir+root+'_c1f.fits' w=readfits(wavefile,hw1) d=readfits(fluxfile,hd1) g=sxpar(hw1,'sdasmgnu') gg=g-1 gr=string(g,'$(i2)') rw=rebin(w,516,g) rd=rebin(d,516,g) wave=rw(*,0) ray=fltarr(516) omega=1.45e-10; pre-costar ray=rd(*,gg)*634.7*wave/omega end