phangs_isrf_pipeline.pro
3.09 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
PRO phangs_isrf_pipeline
;PHANGS ISRF pipeline
;==== make sky images
;produces _jwst_images.sav and _ref_header.sav
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop
;produces _muse_filters_data.sav,_muse_images.sav, _muse_NH.sav
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
;produces _astrosat_data.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop
;produces _herschel_images.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop
phangs_make_co_images,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
phangs_extract_seds,source_name='ngc0628'
;fit_phangs_ngc0628_nir_continuum,/from_restore,/show_plots ;fit seds of some voronoi bins
;fit_phangs_ngc0628_nir_continuum,/from_restore,/brute_force
;fit_phangs_ngc0628_nir_continuum,/from_restore,/grid_brute_force
;==== Check the classes and do a plot
check_phangs_ssps_isrf_prediction,/show_map
plot_phangs_ssps_against_uv_fast
make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save
check_phangs_ssps_isrf_prediction
;==== Make isrf classes
make_phangs_isrf_classes,bidon,source_name='ngc0628',/fits
;make_sed_phangs_tables,/test,/show_seds ;test of plugin dustem_plugin_phangs_class_isrf.pro
;Generally very long. Do it using nohup (see make_phangs_grid.pro)
make_sed_phangs_tables,grid_type=2,isrf_class=15
make_sed_phangs_tables,grid_type=2,isrf_class=0
;...
;=== had to convert previous grids like this:
convert_all_grid_tables_from_v1_to_v2
;=== Fit the SEDs with the ISRF Grids
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/from_restore
resolution_filter='SPIRE3'
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter=resolution_filter
phangs_smooth_muse_isrf,'ngc0628',reso_filter='SPIRE3',/save
phangs_smooth_muse_isrf,'ngc0628',reso_filter='SPIRE1',/save
;=== NGC3351
;Note: faire un chmod a+rwx * apres transfert sur alma1
source_name='ngc3351'
;produces _jwst_images.sav and _ref_header.sav
phangs_make_jwst_images,source_name=source_name,/save,/show,/nostop
;produces _muse_filters_data.sav,_muse_images.sav, _muse_NH.sav
phangs_make_muse_filters_images,source_name=source_name,/save,/show,/nostop
;produces _astrosat_data.sav
phangs_make_astrosat_images,source_name=source_name,/save,/show,/nostop
;produces _herschel_images.sav
phangs_make_herschel_images,source_name=source_name,/save,/show,/nostop
phangs_make_hi_images,source_name=source_name,/save,/show,/nostop
phangs_make_co_images,source_name=source_name,/save,/show,/nostop
phangs_extract_seds,source_name=source_name,/save
make_phangs_ssps_isrf_prediction,source_name=source_name,/save
make_phangs_isrf_classes,source_name=source_name,/fits
phangs_brute_force_fit_with_isrf_grid,source_name=source_name,/include_herschel,/normalize,/fit_G0
END