Commit d3b1715c134e5b2a8741e24f27fb04eaf51da75b

Authored by Jean-Philippe Bernard
1 parent c5ee1d5e
Exists in master

improved

LabTools/IRAP/JPB/read_muse_phangs_weights.pro
... ... @@ -13,7 +13,10 @@ IF keyword_set(young) THEN BEGIN
13 13 dir=!phangs_data_dir+'/phangs_drive/Archive/PHANGS_MUSE/DR2.2/copt/MUSEDAP/young'
14 14 ENDIF
15 15  
16   -file=dir+'/'+use_object+'-0.92asec_ppxf_SFH.fits'
  16 +reso_str=phangs_objectname2reso_str(object)
  17 +
  18 +;file=dir+'/'+use_object+'-0.92asec_ppxf_SFH.fits'
  19 +file=dir+'/'+use_object+'-'+reso_str+'_ppxf_SFH.fits'
17 20 st_info=file_info(file)
18 21 IF st_info.exists NE 1 THEN BEGIN
19 22 message,'Could not find '+file,/continue
... ...
LabTools/IRAP/RM/phangs_make_herschel_images.pro
... ... @@ -82,7 +82,7 @@ data_dir=!phangs_data_dir+'/phangs_drive/Archive/Ancillary/herschel/PHANGS_Hersc
82 82 save_data_dir=!phangs_data_dir+'/ISRF/WORK/'
83 83  
84 84 ;This is to get the reference header
85   -file=save_data_dir+'ngc0628_ref_header'+reso_str+'.sav'
  85 +file=save_data_dir+use_source_name+'_ref_header'+reso_str+'.sav'
86 86 st_info=file_info(file)
87 87 IF st_info.exists NE 1 THEN BEGIN
88 88 message,'Could not find '+file,/continue
... ... @@ -98,6 +98,8 @@ herschel_filters=dustem_filter_names2filters(name_loc_filters)
98 98 Nherschel_filter=n_elements(herschel_filters)
99 99 herschel_images=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'),2,Nherschel_filter)
100 100  
  101 +;stop
  102 +
101 103 ;JPB: added missing SPIRE filters
102 104 name_filters = ['pacs70','pacs100','pacs160','spire250','spire350','spire500']
103 105  
... ... @@ -155,14 +157,18 @@ FOR i=0L,Nherschel_filter-1 DO BEGIN
155 157  
156 158 IF keyword_set(resolution_filter) or keyword_set(resolution_value) THEN BEGIN
157 159 data_reso=dustem_filter2reso(herschel_filters[i]) ;This is undefined in instrument_descritpion.xcat
158   - IF final_reso GT data_reso THEN BEGIN
  160 + IF final_reso GE data_reso THEN BEGIN
159 161 d=degrade_res(d,h,data_reso,final_reso,hout)
160 162 d_err=degrade_res(d_err,h,data_reso,final_reso,hout)
161 163 h=hout
162 164 ENDIF ELSE BEGIN
  165 + message,'Not convolving image as data_reso='+strtrim(data_reso,2)+' > finale_reso ='+strtrim(final_reso,2),/continue
  166 + stop
163 167 ;left blank intentionnally
164 168 ENDELSE
165 169 ENDIF
  170 + print,la_min(d),la_max(d)
  171 + ;stop
166 172  
167 173 IF sxpar(h,'NAXIS1') NE sxpar(href,'NAXIS1') OR sxpar(h,'NAXIS2') NE sxpar(href,'NAXIS2') THEN BEGIN
168 174 d=project2(h,d,href,/silent,save_it='/tmp/saved_proj_info.sav')
... ...
src/idl/dustem_plugin_phangs_stellar_isrf.pro
... ... @@ -6,7 +6,7 @@ FUNCTION dustem_plugin_phangs_stellar_isrf,key=key $
6 6 ,metalicity_values=metalicity_values $
7 7 ,paramdefault=paramdefault $
8 8 ; ,object_distance=object_distance $
9   - ,object_thickness=object_thickness $
  9 +; ,object_thickness=object_thickness $
10 10 ,Voronoi_Npix=Voronoi_Npix $
11 11 ,help=help
12 12  
... ...