From 38087e96a2b7c6168eba0557745f17e27b61afee Mon Sep 17 00:00:00 2001 From: Jean-Philippe Bernard Date: Mon, 26 Feb 2024 16:56:56 +0100 Subject: [PATCH] first commit --- LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro | 14 +++++++++----- 1 file changed, 9 insertions(+), 5 deletions(-) diff --git a/LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro b/LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro index 82049a6..afa8110 100644 --- a/LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro +++ b/LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro @@ -1,8 +1,10 @@ -PRO phangs_compare_seds_isrf +PRO phangs_compare_seds_isrf,isrf_class + +;phangs_compare_seds_isrf,15 dir=!phangs_data_dir+'/ISRF/GRIDS/' -isrf_class_str='_isrfclass15' +isrf_class_str='_isrfclass'+strtrim(isrf_class,2) model='DBP90' table_name=dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs'+isrf_class_str+'.fits' st=mrdfits(table_name,1,h) @@ -14,7 +16,7 @@ Nfilters=n_elements(filters) ii=100L -Nseds=1000L +Nseds=100L Nparams=3 true_params=fltarr(Nseds,Nparams) fitted_params=fltarr(Nseds,Nparams) @@ -24,7 +26,7 @@ FOR j=0L,Nseds-1 DO BEGIN ;seed=1 inn=abs(randomu(seed)) ii=inn*Nst - print,inn,ii + ;print,inn,ii sed=dustem_initialize_sed(Nfilters) ;sed=fltarr(Nfilters) tags=tag_names(st) @@ -53,15 +55,17 @@ ENDFOR stop win=0L -xtit='True parameter (class 15 ISRF)' +xtit='True parameter ('+isrf_class_str+')' ytit='recovered parameter (Mathis ISRF)' window,win & win=win+1 cgplot,true_params[*,0],fitted_params[*,0],tit='G0',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit +cgoplot,!x.crange,!y.crange,linestyle=1 window,win & win=win+1 cgplot,true_params[*,1],fitted_params[*,1],tit='Ypah',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit window,win & win=win+1 cgplot,true_params[*,2],fitted_params[*,2],tit='Yvsg',psym=4,/xlog,/ylog,xtit=xtit,ytit=ytit +cgoplot,10^!x.crange,10^!y.crange,linestyle=2,color='red',thickness=3 END -- libgit2 0.21.2