phangs_isrf_pipeline_alma1.pro
4.16 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
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
;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
;produces _herschel_images.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop
;produces _astrosat_images.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop
;produces _muse_data.sav,
phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop
;==== extract SEDs on MUSE pixels
;produces _jwst_seds_muse_pixels.sav, _muse_seds_muse_pixels.sav', _astrosat_seds_muse_pixels.sav
;Actually running using :
;nohup /usr/local/harris/idl87/bin/idl < {$phangs_srundir}/make_phangs_extract_seds.pro >& {$phangs_srundir}/make_phangs_extract_seds.log &
phangs_extract_seds,source_name='ngc0628'
;=== 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
;produce isrf_classes_one_ratio_on_*.fits, *_isrf_classes_map.fits, *_isrf_classes_voronoi.sav
make_phangs_isrf_classes,bidon,source_name='ngc0628',/save
;=== Fit the SEDs with the ISRF Grids
;=== produce
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save
;============================================================================================================
;===== 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'
;ngc0628_muse_filters_data_PACS3.sav, ngc0628_muse_images_PACS3.sav (not produced _muse_NH_PACS3.sav)
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
;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'
;ngc0628_isrf_min_prediction_PACS3.sav
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS3',/save
;===== j'en suis la sur alma1
;=== 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
END