PRO dustem_cc_mips_vs_handbook

;Example of color correction calculations with DustemWrap
;Computes color corrections for Spitzer MIPS filters
;Results to be compared to the tables in the MIPS HandBook
;http://irsa.ipac.caltech.edu/data/SPITZER/docs/mips/mipsinstrumenthandbook/MIPS_Instrument_Handbook.pdf
;page 98

;== Initialize DustemWrap
dustem_init

;== Define filters
filters=['MIPS1','MIPS2','MIPS3']
Nfilt=n_elements(filters)
;== Define Wavelengths for spectrum definition
Nwav=1000 & wavmin=10. & wavmax=1000.
wavs=findgen(Nwav)/(Nwav-1)*(wavmax-wavmin)+wavmin

;== Compute color corrections for power laws of frequency
pows=[-3.,-2.,-1.,0.,1,2,3]
Nv=n_elements(pows)
K=fltarr(Nv,Nfilt)
frmt1='(A10,15F9.3)'
frmt2='(A10,15I9)'
FOR i=0L,Nv-1 DO BEGIN
  spec=wavs^(-1.*pows(i))
  sed=dustem_cc(wavs,spec,filters,cc=cc,fluxconv=fluxconv,help=help)
  K(i,*)=cc
ENDFOR

;== print results
print,'====== Power laws ======'
print,'Filters',pows,format=frmt1
FOR i=0L,Nfilt-1 DO print,filters(i),K(*,i),format=frmt1

;== Compute color corrections for Black Bodies
Ts=[1.e4,5.e3,1.e3,500,300,200,150,100,70,50,30,20]
Nv=n_elements(Ts)
K=fltarr(Nv,Nfilt)
FOR i=0L,Nv-1 DO BEGIN
  spec=dustem_planck_function(Ts(i),wavs)
  sed=dustem_cc(wavs,spec,filters,cc=cc,fluxconv=fluxconv,help=help)
  K(i,*)=cc
ENDFOR

;== print results
print,'====== Black Bodies ======'
print,'Filters',Ts,format=frmt2
FOR i=0L,Nfilt-1 DO print,filters(i),K(*,i),format=frmt1

END