spec=mrdfits('spSpec-54592-2974-004.fit',0,header)
IH0=where(strmid(header,0,6) eq 'COEFF0')  
IH1=where(strmid(header,0,6) eq 'COEFF1')
reads,strmid(header[IH0],9,22), COEFF0          
reads,strmid(header[IH1],9,22), COEFF1


Imagg=where(strmid(header,0,5) eq 'MAG_G')
Imagr=where(strmid(header,0,5) eq 'MAG_R')
Imagi=where(strmid(header,0,5) eq 'MAG_I')

reads,strmid(header[Imagg],9,22), magg          
reads,strmid(header[Imagr],9,22), magr
reads,strmid(header[Imagi],9,22), magi

print,magg,magr,magi

flux=spec[*,0]*1.d-17    
Nw=N_elements(flux)  
wav=10^(indgen(Nw)*COEFF1+COEFF0)

fmagg=im_filtermag(wav,flux,filterlist='sdss_g0.par')
fmagr=im_filtermag(wav,flux,filterlist='sdss_r0.par')
fmagi=im_filtermag(wav,flux,filterlist='sdss_i0.par')
fmagg=-2.5*alog10(fmagg)
fmagr=-2.5*alog10(fmagr)
fmagi=-2.5*alog10(fmagi)

print,fmagg,fmagr,fmagi

end