axisA.text(5, 0.15, "3", color=colorsA[3], fontsize=fs)
axisA.text(8, 0, "4", color=colorsA[3], fontsize=fs)
axisA.text(5.75, 0.225, "3", color="#8B4513", fontsize=fs)
axisA.text(9, 0.225, "4", color="#8B4513", fontsize=fs)
axisA.text(12.5, 0.125, "5", color=colorsA[4], fontsize=fs)
axisA.text(19.5, 0.08, "6", color=colorsA[4], fontsize=fs)
axisA.text(22, 0.04, "7", color=colorsA[4], fontsize=fs)
axisA.text(12, 0.235, "5", color=colorsA[5], fontsize=fs)
axisA.text(15, 0.235, "6", color=colorsA[5], fontsize=fs)
axisA.text(17.5, 0.235, "7", color=colorsA[5], fontsize=fs)
axisA.text(20, 0.235, "8", color=colorsA[5], fontsize=fs)
axisA.text(23, 0.235, "9", color=colorsA[5], fontsize=fs)
aging_extras.draw_cell(
    axisA, 0.59 * xlim, 0.69 * xlim, 0.9 * ylim, 0.05 * xlim, y2xscale=ylim / xlim, toxic=True, arrow=0.01 * xlim
)
aging_extras.draw_const_env(axisA, 0.79 * xlim, 0.8 * ylim, 0.99 * xlim, ylim, y2xscale=ylim / xlim)
axisA.set_xlim([0, xlim])
axisA.set_xticks([0, 4, 8, 12, 16, 20, 24])
# axisA.set_ylim([0, ylim])
axisA.set_xlabel(r"Time (h)")
axisA.set_ylabel(r"Cellular specific growth rate (h$^{-1}$)")
axisA.set_title("Following Old Pole Cell Over Time")


axisB = fig.add_subplot("B", 222)
paths = ["T010AS12NR100_results.xml", "T010AS12007100_rerun_results.xml", "T010AS12OB100_results.xml"]
attributes = {"name": "specific growth rate population structure", "starting_time": "2400", "header": "bin,frequency"}
colorsB = ["red", "#F781F3", "#fbbf07"]

xlim, ylim = 0.082, 0.6
axisA.text(1,0.56,'Asymmetric, No repair', color=colors[0], va='center',
                                                                   fontsize=fs)
axisA.text(2,0.48,'Asymmetric, Optimal repair', color=colors[1], va='center',
                                                                   fontsize=fs)
axisA.text(1.5,0.51,'0', fontsize=fs)
axisA.text(3.2, 0.395, '1', fontsize=fs)
axisA.text(3,0.25,'2',color=colors[0], fontsize=fs)
axisA.text(6,0.1,'3',color=colors[0], fontsize=fs)
axisA.text(4.5,0.32,'2',color=colors[1], fontsize=fs)
axisA.text(7,0.26,'3',color=colors[1], fontsize=fs)
axisA.text(10,0.21,'4',color=colors[1], fontsize=fs)
axisA.text(14,0.16,'5',color=colors[1], fontsize=fs)
axisA.text(18,0.11,'6',color=colors[1], fontsize=fs)
axisA.text(22,0.03,'7',color=colors[1], fontsize=fs)
aging_extras.draw_cell(axisA, 0.59*xlim, 0.69*xlim, 0.9*ylim, 0.05*xlim, y2xscale=ylim/xlim, toxic=True, arrow=0.01*xlim)
aging_extras.draw_const_env(axisA, 0.79*xlim, 0.8*ylim, 0.99*xlim, ylim, y2xscale=ylim/xlim)
axisA.set_xlim([0, xlim])
axisA.set_xticks([0, 4, 8, 12, 16, 20])#, 24])
axisA.set_ylim([0, ylim])
axisA.set_xlabel(r'Time (h)')
axisA.set_ylabel(r'Cellular specific growth rate (h$^{-1}$)')
axisA.set_title('Following Old Pole Cell Over Time')


axisB = fig.add_subplot('B', 222)
paths = ['an10_const_steady_results.xml', 'ao10_const_steady_results.xml']
attributes = {'name':'specific growth rate population structure',
                    'starting_time':'2400', 'header':'bin,frequency'}

xlim, ylim = 0.082, 0.6