ax_mg.minorticks_on()
mb_ylims = [2.37, 2.84]
mg_ylims = [2.09, 2.55]
xlims = [5, 11]
if eosName == "HShenEOS":
eosName = "HShen"
mb_ylims = [2.401, 3.1]
mg_ylims = [2.21, 2.9]
xlims = [2, 9]
ax_mb.set_ylim(mb_ylims)
ax_mb.set_xlim(xlims)
ax_mg.set_ylim(mg_ylims)
textPos = (0.45, 0.06)
plt.annotate(
"$\\rho_\mathrm{b, max}=\,$%s $\,$ g cm$^{-3}$" % fixScientificNotation(rhob),
textPos,
xytext=textPos,
xycoords="axes fraction",
textcoords="axes fraction",
fontsize=22,
)
textPos = (0.66, 0.16)
plt.annotate(
eosName + " $\,\,\\tilde{A}=%s\,\,$" % a,
textPos,
xytext=textPos,
xycoords="axes fraction",
textcoords="axes fraction",
fontsize=22,
)
c=colors[script],
marker=symbols[script],
ms=6,
lw=lineWidths[script],
dashes=dashList[i],
markeredgecolor=colors[script],
**labelKwarg)
# plt.plot(*mbToroid, c=colors[script], marker=symbols[script], ms=6, lw=lineWidths[script],
# dashes=(20, 5), markeredgecolor=colors[script]
# )
del thisSet
textPos = (0.7, 0.9 - 0.08 * i)
dashPlotLeg.append(plert)
dashLeg.append(r"$\rho_\mathrm{b, max}=\,$%s$\,$g cm$^{-3}$" %
fixScientificNotation(rhob))
# plt.annotate(r"$\rho_\mathrm{b, max}=$%s" % fixScientificNotation(rhob), textPos, xytext=textPos,
# xycoords='axes fraction', textcoords='axes fraction', fontsize=20)
plt.xlabel(r"$\Omega_c$ [$10^3$ rad s$^{-1}$]", labelpad=10)
#plt.axes().yaxis.set_minor_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
#plt.axes().yaxis.set_major_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
plt.ylabel("$\sim \\bar{E}$ [$10^{15}$ g cm$^{-3}$]", labelpad=5)
plt.ylabel("$r_e$ [km]", labelpad=5)
#plt.ylabel("$M_\mathrm{g}$ [$M_\odot$$]")
#removeExponentialNotationOnAxis('y')
legend1 = plt.legend(loc=2, handletextpad=0.2)
plt.minorticks_on()
if eosName == "HShenEOS":
eosName = "HShen"
thisSeq = cstSequence(thisSet, theSlice, filters)
#mbToroid = thisSeq.getSeqPlot([xVar], ['baryMass'], ('RedMax>0.0',), xcolFunc=lambda x: x/1000.0)
mgToroid = thisSeq.getSeqPlot([xVar], yVar, ('RedMax>0.0',), xcolFunc=xFunc, ycolFunc=yFunc)
#if script == 'c30p10':
plt.plot(*mgToroid, c=colors[script], marker=symbols[script], ms=6, lw=lineWidths[script],
dashes=dashList[i], markeredgecolor=colors[script],
**labelKwarg )
# plt.plot(*mbToroid, c=colors[script], marker=symbols[script], ms=6, lw=lineWidths[script],
# dashes=(20, 5), markeredgecolor=colors[script]
# )
del thisSet
textPos = (0.7, 0.9 - 0.08 * i)
dashPlotLeg.append(plert)
dashLeg.append(r"$\rho_\mathrm{b, max}=\,$%s$\,$g cm$^{-3}$" % fixScientificNotation(rhob))
# plt.annotate(r"$\rho_\mathrm{b, max}=$%s" % fixScientificNotation(rhob), textPos, xytext=textPos,
# xycoords='axes fraction', textcoords='axes fraction', fontsize=20)
plt.xlabel(r"$\Omega_c$ [$10^3$ rad s$^{-1}$]", labelpad=10)
#plt.axes().yaxis.set_minor_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
#plt.axes().yaxis.set_major_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
plt.ylabel("$\sim \\bar{E}$ [$10^{15}$ g cm$^{-3}$]", labelpad=5)
plt.ylabel("$r_e$ [km]", labelpad=5)
#plt.ylabel("$M_\mathrm{g}$ [$M_\odot$$]")
#removeExponentialNotationOnAxis('y')
legend1 = plt.legend(loc=2, handletextpad=0.2)
plt.minorticks_on()
if eosName == "HShenEOS":
eosName = "HShen"
del thisSet
plt.xlabel(r"$\Omega_c$ [$10^3$ rad s$^{-1}$]", labelpad=10)
#plt.axes().yaxis.set_minor_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
#plt.axes().yaxis.set_major_formatter(matplotlib.pyplot.FormatStrFormatter('%.0f'))
plt.ylabel("$M_\mathrm{b,g} \,\, [M_\odot]$", labelpad=5)
#removeExponentialNotationOnAxis('y')
plt.legend(loc=2)
plt.minorticks_on()
if eosName == "HShenEOS":
eosName = "HShen"
textPos = (0.32, 0.85)
plt.annotate(eosName, textPos, xytext=textPos, xycoords='axes fraction', textcoords='axes fraction',
fontsize=28)
textPos = (0.32, 0.74)
plt.annotate("$\\tilde{A}=%s$" % a, textPos, xytext=textPos, xycoords='axes fraction', textcoords='axes fraction',
fontsize=28)
textPos = (0.45, 0.075)
plt.annotate(r"$\rho_\mathrm{b, max}=\,$%s$\,$g cm$^{-3}$" % fixScientificNotation(rhob), textPos, xytext=textPos, xycoords='axes fraction', textcoords='axes fraction',
fontsize=25)
textPos = (0.7, 0.3)
plt.annotate("$M_\mathrm{g}$",
textPos, xytext=textPos, xycoords='axes fraction', textcoords='axes fraction',
fontsize=26)
textPos = (0.4, 0.55)
plt.annotate("$M_\mathrm{b}$",
textPos, xytext=textPos, xycoords='axes fraction', textcoords='axes fraction',
fontsize=26)
#plt.xlim(2, 9)
plt.show()
if eosName == "HShenEOS":
eosName = "HShen"
#plt.gca().xaxis.set_major_formatter(locator)
#plt.text(1.2e15, 1.5, eosName, fontsize=26) # Mg LS220
#plt.text(1.2e15, 1.8, eosName, fontsize=26) # Mb LS220
plt.text(10**15, 2.0, eosName, fontsize=26) # Mb Shen
#plt.text(10 **15, 1.8, eosName, fontsize=26) # Mg Shen
matplotlib.rc('xtick', labelsize=20)
matplotlib.rc('xtick.major', pad=6)
matplotlib.rc('ytick', labelsize=20)
inset = plt.axes([0.53, 0.22, 0.41, 0.32]) # OTHER
#inset = plt.axes([0.52, 0.22, 0.42, 0.31]) # Mg Shen
for thePlot in plotList:
#print thePlot
plt.plot(*thePlot[0], c=thePlot[1], ls=thePlot[2], dashes=thePlot[3])
plt.minorticks_on()
locator = matplotlib.ticker.MaxNLocator(3)
locator = matplotlib.ticker.MultipleLocator(3e14)
formatter = matplotlib.ticker.FuncFormatter(
lambda x, y: fixScientificNotation(x, 2))
inset.xaxis.set_major_locator(locator)
inset.xaxis.set_major_formatter(formatter)
#inset.xaxis.set_label_text(fontsize=20)
plt.xlim([10e14, 2.0e15]) # Mb/Mg LS220
plt.xlim([8e14, 1.6e15]) # Mg/Mb Shen
#plt.ylim([2.25, 2.43]) # Mg LS220
plt.ylim([2.25, 2.79]) # Mb LS220
plt.ylim([2.51, 3.05]) # Mb Shen
#plt.ylim([2.53, 2.675]) # Mg Shen
plt.show()
