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"))
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)))