dirname = os.path.join(output_dirname, "steady_states")
if not os.path.exists(dirname):
os.mkdir(dirname)
plt.title(title)
plt.savefig(os.path.join(dirname, title.replace(" ", "-") + ".png"))
# plt.show()
plt.close()
if __name__ == "__main__":
box = DiffusionBox()
# Rates are listed in uM, from table S1
# k_off and kcat rates are already in s^-1; k_on needs to be converted.
rates = [
box.convert_membrane_rate(0.12),
3.0,
box.convert_membrane_rate(0.11),
0.4,
box.convert_membrane_rate(0.05),
0.1,
0.038,
box.convert_membrane_rate(0.07),
1.0,
0.003,
box.convert_cytosolic_rate(1.74),
0.2,
0.1,
0.01,
1.0,
0.01,
