Commit 58e1b6ca8448002367b5c9fbd561c46780a3ec57
1 parent
62eb64c6
Exists in
master
improved
Showing
2 changed files
with
12 additions
and
5 deletions
Show diff stats
LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
| ... | ... | @@ -162,7 +162,8 @@ FOR isrf_class=class_min,class_max DO BEGIN |
| 162 | 162 | IF count EQ 0 THEN BEGIN |
| 163 | 163 | params=dustem_brute_force_fit(sed,table_name,use_filters,fact=fact,chi2=chi2,/normalize,rchi2=rchi2,show_sed=show_sed $ |
| 164 | 164 | ,params_hit=params_hit,params_uncertainties=params_uncertainties,params_min=params_min,params_max=params_max $ |
| 165 | - ,fixed_parameters_description=fixed_parameters_description,fixed_parameters_values=fixed_parameters_values) | |
| 165 | + ,fixed_parameters_description=fixed_parameters_description,fixed_parameters_values=fixed_parameters_values $ | |
| 166 | + ,weighted_params=weighted_params) | |
| 166 | 167 | ;This would be if there was no fixed_parameters |
| 167 | 168 | ;G0s[vid]=params[0] |
| 168 | 169 | ;Ypahs[vid]=params[1]/fact ;extenssive quantities must be divided by normalization factor |
| ... | ... | @@ -245,16 +246,19 @@ cgplot,xv,res,psym=10,title='Ypahs histogram',xtit='Ypah',ytit='Number',/ylog,yr |
| 245 | 246 | window,win,xsize=800,ysize=900 & win=win+1 |
| 246 | 247 | |
| 247 | 248 | obp=[1.1,0.,1.15,1] |
| 248 | -image_cont20,G0_map,Href,/square,imrange=[-2,40],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='G0' | |
| 249 | +;image_cont20,G0_map,Href,/square,imrange=[-2,40],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='G0' | |
| 249 | 250 | image_cont20,la_log10(G0_map),Href,/square,imrange=[-0.2,2],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='G0' |
| 250 | 251 | |
| 252 | +window,win,xsize=800,ysize=900 & win=win+1 | |
| 251 | 253 | image_cont20,la_log10(Ypah_map),Href,/square,imrange=[-2.,6],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Ypah)' |
| 252 | 254 | |
| 255 | +window,win,xsize=800,ysize=900 & win=win+1 | |
| 253 | 256 | image_cont20,la_log10(YVSG_map),Href,/square,imrange=[-2.,6],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Yvsg)' |
| 254 | 257 | |
| 255 | -image_cont20,chi2_map,Href,/square,imrange=[0.01,100]*1.e4,image_color_table='jpbloadct',/silent | |
| 258 | +window,win,xsize=800,ysize=900 & win=win+1 | |
| 259 | +image_cont20,chi2_map,Href,/square,imrange=[0.01,100]*1.e4,image_color_table='jpbloadct',/silent,title='Chi2' | |
| 256 | 260 | |
| 257 | -;stop | |
| 261 | +stop | |
| 258 | 262 | jwst_4_coutours=fltarr((size(jwst_images))[1],(size(jwst_images))[2],2) |
| 259 | 263 | jwst_4_coutours[*,*,0]=jwst_images[*,*,0,5] |
| 260 | 264 | jwst_4_coutours[*,*,1]=jwst_images[*,*,0,5] | ... | ... |
LabTools/IRAP/JPB/phangs_compare_seds_isrf.pro
| ... | ... | @@ -52,7 +52,10 @@ FOR j=0L,Nseds-1 DO BEGIN |
| 52 | 52 | ,params_hit=params_hit,params_uncertainties=params_uncertainties,params_min=params_min,params_max=params_max $ |
| 53 | 53 | ,fixed_parameters_description=fixed_parameters_description,fixed_parameters_values=fixed_parameters_values $ |
| 54 | 54 | ,weighted_params=weighted_params) |
| 55 | - | |
| 55 | + params[1]=params[1]/fact | |
| 56 | + params[2]=params[2]/fact | |
| 57 | + weighted_params[1]=weighted_params[1]/fact | |
| 58 | + weighted_params[2]=weighted_params[2]/fact | |
| 56 | 59 | ;fitted_params[j,*]=params |
| 57 | 60 | IF keyword_set(weighted) THEN BEGIN |
| 58 | 61 | fitted_params[j,*]=weighted_params | ... | ... |