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

;=== 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