Commit 38087e96a2b7c6168eba0557745f17e27b61afee
1 parent
3a667b1d
Exists in
master
first commit
Showing
1 changed file
with
9 additions
and
5 deletions
Show diff stats
LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro
1 | 1 | |
2 | -PRO phangs_compare_seds_isrf | |
2 | +PRO phangs_compare_seds_isrf,isrf_class | |
3 | + | |
4 | +;phangs_compare_seds_isrf,15 | |
3 | 5 | |
4 | 6 | dir=!phangs_data_dir+'/ISRF/GRIDS/' |
5 | -isrf_class_str='_isrfclass15' | |
7 | +isrf_class_str='_isrfclass'+strtrim(isrf_class,2) | |
6 | 8 | model='DBP90' |
7 | 9 | table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs'+isrf_class_str+'.fits' |
8 | 10 | st=mrdfits(table_name,1,h) |
... | ... | @@ -14,7 +16,7 @@ Nfilters=n_elements(filters) |
14 | 16 | |
15 | 17 | ii=100L |
16 | 18 | |
17 | -Nseds=1000L | |
19 | +Nseds=100L | |
18 | 20 | Nparams=3 |
19 | 21 | true_params=fltarr(Nseds,Nparams) |
20 | 22 | fitted_params=fltarr(Nseds,Nparams) |
... | ... | @@ -24,7 +26,7 @@ FOR j=0L,Nseds-1 DO BEGIN |
24 | 26 | ;seed=1 |
25 | 27 | inn=abs(randomu(seed)) |
26 | 28 | ii=inn*Nst |
27 | - print,inn,ii | |
29 | + ;print,inn,ii | |
28 | 30 | sed=dustem_initialize_sed(Nfilters) |
29 | 31 | ;sed=fltarr(Nfilters) |
30 | 32 | tags=tag_names(st) |
... | ... | @@ -53,15 +55,17 @@ ENDFOR |
53 | 55 | stop |
54 | 56 | |
55 | 57 | win=0L |
56 | -xtit='True parameter (class 15 ISRF)' | |
58 | +xtit='True parameter ('+isrf_class_str+')' | |
57 | 59 | ytit='recovered parameter (Mathis ISRF)' |
58 | 60 | |
59 | 61 | window,win & win=win+1 |
60 | 62 | cgplot,true_params[*,0],fitted_params[*,0],tit='G0',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit |
63 | +cgoplot,!x.crange,!y.crange,linestyle=1 | |
61 | 64 | window,win & win=win+1 |
62 | 65 | cgplot,true_params[*,1],fitted_params[*,1],tit='Ypah',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit |
63 | 66 | window,win & win=win+1 |
64 | 67 | cgplot,true_params[*,2],fitted_params[*,2],tit='Yvsg',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit |
68 | +cgoplot,10^!x.crange,10^!y.crange,linestyle=2,color='red',thickness=3 | |
65 | 69 | |
66 | 70 | END |
67 | 71 | ... | ... |