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