dustem_cc_mips_vs_handbook.pro
1.42 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
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