OV-GSO-DC / dustem-wrapper_idl

Blame view

LabTools/IRAP/JPB/phangs_isrf_pipeline_alma1.pro 4.16 KB
Edit Raw Normal View History
bd9f60b0   Jean-Philippe Bernard   First commit 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
PRO phangs_isrf_pipeline_alma1


;PHANGS ISRF pipeline
save_data_dir=!phangs_data_dir+'/ISRF/WORK/'
str='ls '+save_data_dir+'ngc0628_*SPIRE3.sav'
print,str

str='ls '+save_data_dir+'ngc0628_*images*.sav'
spawn,str
str='ls '+save_data_dir+'ngc0628_*seds_muse_pixels*.sav'
spawn,str

;==== make sky images
;produces _jwst_images.sav and _ref_header.sav
;/data/projects/phangs/phangs_drive/Archive/PHANGS_JWST/DR1/ngc0628/f1000w_psf_matched
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop
;produces _muse_filters_images.sav,_muse_data.sav, _muse_NH.sav
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
b09a906c   Jean-Philippe Bernard   improved 
20
21
22
23
 
;produces _HI_images.sav
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop
;produces _CO_images.sav
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop
06c3ba9c   Jean-Philippe Bernard   improved 
24
25
 
;produces _herschel_images.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop
bd9f60b0   Jean-Philippe Bernard   First commit 
26
27
 
;produces _astrosat_images.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop
bd9f60b0   Jean-Philippe Bernard   First commit 
28
9c5f0ed1   Jean-Philippe Bernard   improved 
29
30
31
 
;produces _muse_data.sav,
phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop
f5f6cb5e   Jean-Philippe Bernard   improved 
32
33
 
;==== extract SEDs on MUSE pixels
;produces _jwst_seds_muse_pixels.sav, _muse_seds_muse_pixels.sav', _astrosat_seds_muse_pixels.sav
4233645a   Jean-Philippe Bernard   improved 
34
35
 
;Actually running using :
;nohup /usr/local/harris/idl87/bin/idl < {$phangs_srundir}/make_phangs_extract_seds.pro >& {$phangs_srundir}/make_phangs_extract_seds.log &
f5f6cb5e   Jean-Philippe Bernard   improved 
36
37
 
phangs_extract_seds,source_name='ngc0628'
dc862aaf   Jean-Philippe Bernard   improved 
38
39
40
41
42
 
;=== Compute ISRF prediction. This is done only at full resolution
;=== produce _isrf_min_prediction.sav, _isrf_min_prediction.fits (product)
make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save

;==== Make isrf classes
8cfdf0b6   Jean-Philippe Bernard   improved 
43
 
;produce isrf_classes_one_ratio_on_*.fits, *_isrf_classes_map.fits, *_isrf_classes_voronoi.sav
dc862aaf   Jean-Philippe Bernard   improved 
44
45
 
make_phangs_isrf_classes,bidon,source_name='ngc0628',/save
dc862aaf   Jean-Philippe Bernard   improved 
46
 
;=== Fit the SEDs with the ISRF Grids
8cfdf0b6   Jean-Philippe Bernard   improved 
47
 
;=== produce
dc862aaf   Jean-Philippe Bernard   improved 
48
49
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save
8cfdf0b6   Jean-Philippe Bernard   improved 
50
51
52
53
54
55
56
57
58
59
60
61
 
;============================================================================================================
;===== NGC0628 at PACS3 resolution
;============================================================================================================
;resolution_filter='PACS3'
;ngc0628_jwst_images_PACS3.sav, ngc0628_ref_header_PACS3.sav
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
;ngc0628_CO_images_PACS3.sav
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
;ngc0628_HI_images_PACS3.sav
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
;ngc0628_astrosat_data_PACS3.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
d175a44a   Jean-Philippe Bernard   improved 
62
 
;ngc0628_muse_filters_data_PACS3.sav, ngc0628_muse_images_PACS3.sav (not produced _muse_NH_PACS3.sav)
8cfdf0b6   Jean-Philippe Bernard   improved 
63
 
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
8cfdf0b6   Jean-Philippe Bernard   improved 
64
65
66
67
68
69
 
;ngc0628_herschel_images_PACS3.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'

phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'

phangs_extract_seds,source_name='ngc0628',resolution_filter='PACS3'
dc862aaf   Jean-Philippe Bernard   improved 
70
8cfdf0b6   Jean-Philippe Bernard   improved 
71
72
 
;ngc0628_isrf_min_prediction_PACS3.sav
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS3',/save
4233645a   Jean-Philippe Bernard   improved 
73
c7d22c08   Jean-Philippe Bernard   improved 
74
75
 
;===== j'en suis la sur alma1
eb7ea183   Jean-Philippe Bernard   included stellar ... 
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
 
;=== NGC3351

phangs_make_jwst_images,source_name='ngc3351',/save,/show,/nostop
;produces _muse_filters_data.sav,_muse_images.sav, _muse_NH.sav
phangs_make_muse_filters_images,source_name='ngc3351',/save,/show,/nostop
;produces _astrosat_data.sav
phangs_make_astrosat_images,source_name='ngc3351',/save,/show,/nostop
;produces _herschel_images.sav
phangs_make_herschel_images,source_name='ngc3351',/save,/show,/nostop
phangs_make_hi_images,source_name='ngc3351',/save,/show,/nostop
phangs_make_co_images,source_name='ngc3351',/save,/show,/nostop

phangs_extract_seds,source_name='ngc3351',/save

make_phangs_ssps_isrf_prediction,source_name='ngc3351',/save

make_phangs_isrf_classes,source_name='ngc3351',/fits

phangs_brute_force_fit_with_isrf_grid,source_name='ngc3351',/include_herschel,/normalize,/fit_G0
bd9f60b0   Jean-Philippe Bernard   First commit 
97
 
END