Exemplo n.º 1
0
           drawstyle='steps-mid')
    # p.errorbar(rs_n.value, sd_n.value, yerr=sd_sigma_n.value, fmt="D-",
    #            color="b", label="North")
    # p.errorbar(rs_s.value, sd_s.value, yerr=sd_sigma_s.value, fmt="o-",
    #            color="g", label="South")
    p.ylabel(r"$\Sigma$ (M$_{\odot}$ pc$^{-2}$)")
    p.xlabel(r"R (kpc)")
    p.legend(loc='best')
    p.grid("on")
    p.savefig(paper1_figures_path("M33_Sigma_profile_N_S.pdf"))
    p.savefig(paper1_figures_path("M33_Sigma_profile_N_S.png"))
    p.clf()

    # Compare to the surface density profile in Corbelli
    corbelli_filename = \
        c_hi_analysispath("rotation_curves/corbelli_rotation_curve.csv")
    corbelli = Table.read(corbelli_filename)

    p.plot(rs.value, sd.value / scale_factor,
           linestyle="-", color="b",
           label="This work", drawstyle='steps-mid')
    p.plot(corbelli["R"][corbelli["R"] <= 10.0],
           corbelli["SigmaHI"][corbelli["R"] <= 10.0], "r--",
           drawstyle='steps-mid',
           label="Corbelli et al. (2014)")
    p.ylabel(r"$\Sigma$ (M$_{\odot}$ pc$^{-2}$)")
    p.xlabel(r"R (kpc)")
    p.legend(loc='best')
    p.grid()
    p.savefig(paper1_figures_path("M33_Sigma_profile_w_Corbelli.pdf"))
    p.savefig(paper1_figures_path("M33_Sigma_profile_w_Corbelli.png"))
Exemplo n.º 2
0
p.xlabel("$\Sigma_{\mathrm{Gas}}$ (M$_{\odot}$ pc$^{-2}$)")
p.ylabel("H$_2$-to-HI Ratio $\Sigma_{\mathrm{H2}} / \Sigma_{\mathrm{HI}}$")
p.xlim([2, 22])
p.ylim([1e-4, 10])
p.legend(loc='upper left')
p.grid()

p.savefig(paper1_figures_path("ratio_totalsigma_w_krumholzmodel_dr_{}pc.pdf".format(dr.value)))
p.savefig(paper1_figures_path("ratio_totalsigma_w_krumholzmodel_dr_{}pc.png".format(dr.value)))
p.close()

# p.show()

# Also plot the total gas surface density against the stellar surface density
# from Corbelli
corbelli = Table.read(c_hi_analysispath("rotation_curves/corbelli_rotation_curve.csv"))

p.semilogy(rs.value, total_sd,
           linestyle="-", color="b",
           label="Gas", drawstyle='steps-mid')
p.semilogy(corbelli["R"][corbelli["R"] <= 6.5],
           corbelli["SigmaStellar"][corbelli["R"] <= 6.5], "g--",
           drawstyle='steps-mid',
           label="Stars")
p.ylabel(r"log $\Sigma$ / (M$_{\odot}$ pc$^{-2}$)")
p.xlabel(r"R (kpc)")
p.legend(loc='best')
p.grid()

p.savefig(paper1_figures_path("M33_Sigma_profile_gas_stars_corbelli_{}pc.pdf".format(dr.value)))
p.savefig(paper1_figures_path("M33_Sigma_profile_gas_stars_corbelli_{}pc.png".format(dr.value)))