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