Commit 6ccbd4f5 authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Centered SF cores positions on peak

parent d7fb9166

etl/mdl/mdl_extract_ch3oh_spectra.py

+1 −1
Changes for etl/mdl/mdl_extract_ch3oh_spectra.py: 1 added line, 1 removed line.
Original line number Diff line number Diff line
@@ -751,7 +751,7 @@ def overlay_regions_density(source, regions, densities, moment_map, wcs): 
        pc.set_facecolor('none')

        pc.set_clim(4.8, 5.8)

        ax.add_collection(pc)

    coord = SkyCoord(core_data_filtered['RA_1'], core_data_filtered['DEC_1'], frame='icrs')

    coord = SkyCoord(core_data_filtered['RA_peak'], core_data_filtered['DEC_peak'], frame='icrs')

    ax.scatter(coord.ra.deg, coord.dec.deg, transform=ax.get_transform('world'), s=100, marker='*',

               edgecolor='lightcyan', facecolor='none', zorder=1, linewidths=1)

    plt.ylabel('Dec')

etl/prs/prs_mom0_dust.py

+1 −1
Changes for etl/prs/prs_mom0_dust.py: 1 added line, 1 removed line.
Original line number Diff line number Diff line
@@ -85,7 +85,7 @@ def save_fits(data, header, filename): 
def plot_swapped(moment0, moment_hdr, continuum, cont_wcs, source, species, cores_data):

    # Background: continuum; Contours: moment0

    wcs2d = WCS(moment_hdr)

    coord = SkyCoord(cores_data['RA_1'], cores_data['DEC_1'], frame='icrs')

    coord = SkyCoord(cores_data['RA_peak'], cores_data['DEC_peak'], frame='icrs')



    plt.figure()

    ax = plt.subplot(projection=cont_wcs)