#mulquiney1



#reactions
v_RBCvAK:
 {1.0}mgadpRBC + {1.0}adpRBC = {1.0}mgatpRBC + {1.0}ampRBC
  TsRBC * Vrbc * (K1appRBCvAK * adpRBC / Vrbc * mgadpRBC / Vrbc - K2appRBCvAK * ampRBC / Vrbc * mgatpRBC / Vrbc)
v_RBCvALD:
 {1.0}f16p2RBC = {1.0}dhapRBC + {1.0}gapRBC
  TsRBC * Vrbc * (VRBCvALD * (f16p2RBC / Vrbc * KcatfRBCvALD / Kmf16p2RBCvALD - gapRBC / Vrbc * dhapRBC / Vrbc * KcatrRBCvALD / (KidhapRBCvALD * KmgapRBCvALD)) / (1 + dhapRBC / Vrbc / KidhapRBCvALD + f16p2RBC / Vrbc / Kmf16p2RBCvALD + gapRBC / Vrbc * dhapRBC / Vrbc / (KidhapRBCvALD * KmgapRBCvALD) + f16p2RBC / Vrbc * gapRBC / Vrbc * KmdhapRBCvALD / (Kif16p2RBCvALD * KidhapRBCvALD * KmgapRBCvALD) + (b23pgRBC / Vrbc + mgb23pgRBC / Vrbc) / Kib23pgRBCvALD + gapRBC / Vrbc * KmdhapRBCvALD * (1 + (b23pgRBC / Vrbc + mgb23pgRBC / Vrbc) / Kib23pgRBCvALD) / (KidhapRBCvALD * KmgapRBCvALD)))
v_RBCvATPASE:
 {1.0}mgatpRBC = {1.0}mgadpRBC + {1.0}phosRBC
  TsRBC * Vrbc * (KRBCvATPASE * mgatpRBC / Vrbc)
v_RBCvBPGSP1:
 {1.0}b13pgRBC + {1.0}bpgspRBC = {1.0}bpgspb13pgRBC
  TsRBC * Vrbc * (K1appRBCvBPGSP1 * b13pgRBC / Vrbc * bpgspRBC / Vrbc - K2RBCvBPGSP1 * bpgspb13pgRBC / Vrbc)
v_RBCvBPGSP2:
 {1.0}bpgspb13pgRBC = {1.0}bpgsppRBC + {1.0}p3gRBC
  TsRBC * Vrbc * (K3appRBCvBPGSP2 * bpgspb13pgRBC / Vrbc)
v_RBCvBPGSP3:
 {1.0}bpgsppRBC + {1.0}p3gRBC = {1.0}bpgsppp3gRBC
  TsRBC * Vrbc * (-K5RBCvBPGSP3 * bpgsppp3gRBC / Vrbc + K4appRBCvBPGSP3 * bpgsppRBC / Vrbc * p3gRBC / Vrbc)
v_RBCvBPGSP4:
 {1.0}bpgsppRBC + {1.0}p2gRBC = {1.0}bpgsppp2gRBC
  TsRBC * Vrbc * (-K7RBCvBPGSP4 * bpgsppp2gRBC / Vrbc + K6appRBCvBPGSP4 * bpgsppRBC / Vrbc * p2gRBC / Vrbc)
v_RBCvBPGSP5:
 {1.0}bpgsppp3gRBC = {1.0}bpgspb23pgRBC
  TsRBC * Vrbc * (K8RBCvBPGSP5 * bpgsppp3gRBC / Vrbc - K9RBCvBPGSP5 * bpgspb23pgRBC / Vrbc)
v_RBCvBPGSP6:
 {1.0}bpgsppp2gRBC = {1.0}bpgspb23pgRBC
  TsRBC * Vrbc * (K10RBCvBPGSP6 * bpgsppp2gRBC / Vrbc - K11RBCvBPGSP6 * bpgspb23pgRBC / Vrbc)
v_RBCvBPGSP7:
 {1.0}bpgspb23pgRBC = {1.0}bpgspRBC + {1.0}b23pgRBC
  TsRBC * Vrbc * (-K13appRBCvBPGSP7 * b23pgRBC / Vrbc * bpgspRBC / Vrbc + K12RBCvBPGSP7 * bpgspb23pgRBC / Vrbc)
v_RBCvBPGSP8:
 {1.0}bpgsppRBC + {1.0}phosRBC = {1.0}bpgsppphosRBC
  TsRBC * Vrbc * (K14RBCvBPGSP8 * bpgsppRBC / Vrbc * phosRBC / Vrbc - K15RBCvBPGSP8 * bpgsppphosRBC / Vrbc)
v_RBCvBPGSP9:
 {1.0}bpgsppphosRBC = {1.0}bpgspRBC + {2.0}phosRBC
  TsRBC * Vrbc * (K16RBCvBPGSP9 * bpgsppphosRBC / Vrbc)
v_RBCvENO:
 {1.0}p2gRBC = {1.0}pepRBC
  TsRBC * Vrbc * (VRBCvENO * (KcatfRBCvENO * mgRBC / Vrbc * p2gRBC / Vrbc / (KimgRBCvENO * Kmp2gRBCvENO) - KcatrRBCvENO * mgRBC / Vrbc * pepRBC / Vrbc / (KipepRBCvENO * KmmgRBCvENO)) / (1 + mgRBC / Vrbc / KimgRBCvENO + p2gRBC / Vrbc / Kip2gRBCvENO + mgRBC / Vrbc * p2gRBC / Vrbc / (KimgRBCvENO * Kmp2gRBCvENO) + pepRBC / Vrbc / KipepRBCvENO + mgRBC / Vrbc * pepRBC / Vrbc / (KipepRBCvENO * KmmgRBCvENO)))
v_RBCvG6PDH:
 {1.0}g6pRBC + {1.0}nadpRBC = {1.0}nadphRBC + {1.0}p6glRBC
  TsRBC * Vrbc * (VRBCvG6PDH * (g6pRBC / Vrbc * K1RBCvG6PDH * K3RBCvG6PDH * K5RBCvG6PDH * K7RBCvG6PDH * K9RBCvG6PDH * nadpRBC / Vrbc - K2RBCvG6PDH * K4RBCvG6PDH * K6RBCvG6PDH * K8RBCvG6PDH * K10RBCvG6PDH * nadphRBC / Vrbc * p6glRBC / Vrbc) / (g6pRBC / Vrbc * K3RBCvG6PDH * K5RBCvG6PDH * K7RBCvG6PDH * K9RBCvG6PDH + K2RBCvG6PDH * (K4RBCvG6PDH * K6RBCvG6PDH + K4RBCvG6PDH * K7RBCvG6PDH + K5RBCvG6PDH * K7RBCvG6PDH) * K9RBCvG6PDH + K1RBCvG6PDH * (K4RBCvG6PDH * K6RBCvG6PDH + K4RBCvG6PDH * K7RBCvG6PDH + K5RBCvG6PDH * K7RBCvG6PDH) * K9RBCvG6PDH * nadpRBC / Vrbc + g6pRBC / Vrbc * K1RBCvG6PDH * K3RBCvG6PDH * (K5RBCvG6PDH * K7RBCvG6PDH + K5RBCvG6PDH * K9RBCvG6PDH + K6RBCvG6PDH * K9RBCvG6PDH + K7RBCvG6PDH * K9RBCvG6PDH) * nadpRBC / Vrbc + g6pRBC / Vrbc * K3RBCvG6PDH * K5RBCvG6PDH * K7RBCvG6PDH * K10RBCvG6PDH * nadphRBC / Vrbc + K2RBCvG6PDH * (K4RBCvG6PDH * K6RBCvG6PDH + K5RBCvG6PDH * K6RBCvG6PDH + K5RBCvG6PDH * K7RBCvG6PDH) * K10RBCvG6PDH * nadphRBC / Vrbc + K2RBCvG6PDH * K4RBCvG6PDH * K6RBCvG6PDH * K8RBCvG6PDH * p6glRBC / Vrbc + K1RBCvG6PDH * K4RBCvG6PDH * K6RBCvG6PDH * K8RBCvG6PDH * nadpRBC / Vrbc * p6glRBC / Vrbc + g6pRBC / Vrbc * K1RBCvG6PDH * K3RBCvG6PDH * (K5RBCvG6PDH + K6RBCvG6PDH) * K8RBCvG6PDH * nadpRBC / Vrbc * p6glRBC / Vrbc + g6pRBC / Vrbc * K3RBCvG6PDH * (K5RBCvG6PDH + K6RBCvG6PDH) * K8RBCvG6PDH * K10RBCvG6PDH * nadphRBC / Vrbc * p6glRBC / Vrbc + (K2RBCvG6PDH * K4RBCvG6PDH + K2RBCvG6PDH * K5RBCvG6PDH + K2RBCvG6PDH * K6RBCvG6PDH + K4RBCvG6PDH * K6RBCvG6PDH) * K8RBCvG6PDH * K10RBCvG6PDH * nadphRBC / Vrbc * p6glRBC / Vrbc))
v_RBCvGAPDH:
 {1.0}gapRBC + {1.0}phosRBC + {1.0}nadRBC = {1.0}nadhRBC + {1.0}b13pgRBC
  TsRBC * Vrbc * VRBCvGAPDH / (gapRBC / Vrbc / KiappgapRBCvGAPDH * (1 + gapRBC / Vrbc / KidgapRBCvGAPDH) + b13pgRBC / Vrbc / Kiappb13pgRBCvGAPDH * (1 + gapRBC / Vrbc / KidgapRBCvGAPDH) + Kmb13pgRBCvGAPDH * nadhRBC / Vrbc / (Kiappb13pgRBCvGAPDH * KmappnadhRBCvGAPDH) + KmgapRBCvGAPDH * nadRBC / Vrbc * phosRBC / Vrbc / (KiappgapRBCvGAPDH * KmnadRBCvGAPDH * KiphosRBCvGAPDH) + nadRBC / Vrbc * gapRBC / Vrbc / (KinadRBCvGAPDH * KiappgapRBCvGAPDH) + phosRBC / Vrbc * gapRBC / Vrbc / (KiphosRBCvGAPDH * KiappgapRBCvGAPDH) * (1 + gapRBC / Vrbc / KidgapRBCvGAPDH) + nadRBC / Vrbc * b13pgRBC / Vrbc / (KinadRBCvGAPDH * Kiappb13pgRBCvGAPDH) + Kmb13pgRBCvGAPDH * phosRBC / Vrbc * nadhRBC / Vrbc / (KiphosRBCvGAPDH * Kiappb13pgRBCvGAPDH * KmappnadhRBCvGAPDH) + gapRBC / Vrbc * nadhRBC / Vrbc / (KiappgapRBCvGAPDH * KiappnadhRBCvGAPDH) + b13pgRBC / Vrbc * nadhRBC / Vrbc / (Kiappb13pgRBCvGAPDH * KmappnadhRBCvGAPDH) + nadRBC / Vrbc * phosRBC / Vrbc * gapRBC / Vrbc / (KmnadRBCvGAPDH * KiphosRBCvGAPDH * KiappgapRBCvGAPDH) + KmgapRBCvGAPDH * nadRBC / Vrbc * phosRBC / Vrbc * b13pgRBC / Vrbc / (KiappgapRBCvGAPDH * KmnadRBCvGAPDH * KiphosRBCvGAPDH * Kidb13pgRBCvGAPDH) + phosRBC / Vrbc * gapRBC / Vrbc * nadhRBC / Vrbc / (KiphosRBCvGAPDH * KiappgapRBCvGAPDH * KiappnadhRBCvGAPDH) + Kmb13pgRBCvGAPDH * phosRBC / Vrbc * b13pgRBC / Vrbc * nadhRBC / Vrbc / (Kiappb13pgRBCvGAPDH * KmappnadhRBCvGAPDH * KiphosRBCvGAPDH * Kidb13pgRBCvGAPDH)) * (KcatfappRBCvGAPDH * gapRBC / Vrbc * nadRBC / Vrbc * phosRBC / Vrbc / (KmnadRBCvGAPDH * KiphosRBCvGAPDH * KiappgapRBCvGAPDH) - KcatrappRBCvGAPDH * b13pgRBC / Vrbc * nadhRBC / Vrbc / (Kiappb13pgRBCvGAPDH * KmappnadhRBCvGAPDH))
v_RBCvGSSGR:
 {1.0}nadphRBC + {1.0}gssgRBC = {2.0}gshRBC + {1.0}nadpRBC
  TsRBC * Vrbc * (VRBCvGSSGR * (-gshRBC * gshRBC / Vrbc / Vrbc * K2RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K10RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gssgRBC / Vrbc * K1RBCvGSSGR * K3RBCvGSSGR * K5RBCvGSSGR * K7RBCvGSSGR * K9RBCvGSSGR * K11RBCvGSSGR * nadphRBC / Vrbc) / (gshRBC * gshRBC / Vrbc / Vrbc * K2RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K10RBCvGSSGR + gshRBC / Vrbc * K2RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K11RBCvGSSGR + gssgRBC / Vrbc * K3RBCvGSSGR * K5RBCvGSSGR * K7RBCvGSSGR * K9RBCvGSSGR * K11RBCvGSSGR + K2RBCvGSSGR * (K4RBCvGSSGR * K6RBCvGSSGR + K4RBCvGSSGR * K7RBCvGSSGR + K5RBCvGSSGR * K7RBCvGSSGR) * K9RBCvGSSGR * K11RBCvGSSGR + gshRBC / Vrbc * K2RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gssgRBC / Vrbc * K3RBCvGSSGR * K5RBCvGSSGR * K7RBCvGSSGR * K9RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + K2RBCvGSSGR * (K4RBCvGSSGR * K6RBCvGSSGR + K4RBCvGSSGR * K7RBCvGSSGR + K5RBCvGSSGR * K7RBCvGSSGR) * K9RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gshRBC / Vrbc * gssgRBC / Vrbc * K3RBCvGSSGR * K5RBCvGSSGR * K7RBCvGSSGR * K10RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gshRBC / Vrbc * K2RBCvGSSGR * (K4RBCvGSSGR * K6RBCvGSSGR + K4RBCvGSSGR * K7RBCvGSSGR + K5RBCvGSSGR * K7RBCvGSSGR) * K10RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gshRBC * gshRBC / Vrbc / Vrbc * gssgRBC / Vrbc * K3RBCvGSSGR * (K5RBCvGSSGR + K6RBCvGSSGR) * K8RBCvGSSGR * K10RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gshRBC * gshRBC / Vrbc / Vrbc * (K2RBCvGSSGR * K4RBCvGSSGR + K2RBCvGSSGR * K5RBCvGSSGR + K2RBCvGSSGR * K6RBCvGSSGR + K4RBCvGSSGR * K6RBCvGSSGR) * K8RBCvGSSGR * K10RBCvGSSGR * K12RBCvGSSGR * nadpRBC / Vrbc + gshRBC / Vrbc * gssgRBC / Vrbc * K1RBCvGSSGR * K3RBCvGSSGR * K5RBCvGSSGR * K7RBCvGSSGR * K10RBCvGSSGR * nadphRBC / Vrbc + gshRBC * gshRBC / Vrbc / Vrbc * K1RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K10RBCvGSSGR * nadphRBC / Vrbc + gshRBC * gshRBC / Vrbc / Vrbc * gssgRBC / Vrbc * K1RBCvGSSGR * K3RBCvGSSGR * (K5RBCvGSSGR + K6RBCvGSSGR) * K8RBCvGSSGR * K10RBCvGSSGR * nadphRBC / Vrbc + gshRBC / Vrbc * K1RBCvGSSGR * K4RBCvGSSGR * K6RBCvGSSGR * K8RBCvGSSGR * K11RBCvGSSGR * nadphRBC / Vrbc + gshRBC / Vrbc * gssgRBC / Vrbc * K1RBCvGSSGR * K3RBCvGSSGR * (K5RBCvGSSGR + K6RBCvGSSGR) * K8RBCvGSSGR * K11RBCvGSSGR * nadphRBC / Vrbc + K1RBCvGSSGR * (K4RBCvGSSGR * K6RBCvGSSGR + K4RBCvGSSGR * K7RBCvGSSGR + K5RBCvGSSGR * K7RBCvGSSGR) * K9RBCvGSSGR * K11RBCvGSSGR * nadphRBC / Vrbc + gssgRBC / Vrbc * K1RBCvGSSGR * K3RBCvGSSGR * (K5RBCvGSSGR * K7RBCvGSSGR * K9RBCvGSSGR + K5RBCvGSSGR * K7RBCvGSSGR * K11RBCvGSSGR + K5RBCvGSSGR * K9RBCvGSSGR * K11RBCvGSSGR + K6RBCvGSSGR * K9RBCvGSSGR * K11RBCvGSSGR + K7RBCvGSSGR * K9RBCvGSSGR * K11RBCvGSSGR) * nadphRBC / Vrbc))
v_RBCvHBADP:
 {1.0}hbRBC + {1.0}adpRBC = {1.0}hbadpRBC
  TsRBC * Vrbc * (-hbadpRBC / Vrbc * KdRBCvHBADP + KaappRBCvHBADP * adpRBC / Vrbc * hbRBC / Vrbc)
v_RBCvHBATP:
 {1.0}hbRBC + {1.0}atpRBC = {1.0}hbatpRBC
  TsRBC * Vrbc * (-hbatpRBC / Vrbc * KdRBCvHBATP + KaappRBCvHBATP * atpRBC / Vrbc * hbRBC / Vrbc)
v_RBCvHBB13PG:
 {1.0}hbRBC + {1.0}b13pgRBC = {1.0}hbb13pgRBC
  TsRBC * Vrbc * (-hbb13pgRBC / Vrbc * KdRBCvHBB13PG + KaappRBCvHBB13PG * b13pgRBC / Vrbc * hbRBC / Vrbc)
v_RBCvHBB23PG:
 {1.0}hbRBC + {1.0}b23pgRBC = {1.0}hbb23pgRBC
  TsRBC * Vrbc * (-hbb23pgRBC / Vrbc * KdRBCvHBBPG + KaappRBCvHBBPG * b23pgRBC / Vrbc * hbRBC / Vrbc)
v_RBCvHBMGATP:
 {1.0}mgatpRBC + {1.0}hbRBC = {1.0}hbmgatpRBC
  TsRBC * Vrbc * (-hbmgatpRBC / Vrbc * KdRBCvHBMGATP + KaappRBCvHBMGATP * hbRBC / Vrbc * mgatpRBC / Vrbc)
v_RBCvHK:
 {1.0}mgatpRBC = {1.0}g6pRBC + {1.0}mgadpRBC
  TsRBC * Vrbc * (-KcatrappRBCvHK * ERBCvHK * g6pRBC / Vrbc * mgadpRBC / Vrbc / (Kig6pRBCvHK * KmmgadpRBCvHK) + KcatfappRBCvHK * ERBCvHK * glcRBC / Vrbc * mgatpRBC / Vrbc / (KiglcRBCvHK * KmmgatpRBCvHK)) / (1 + glcRBC / Vrbc / KiglcRBCvHK + glcRBC / Vrbc * b23pgRBC / Vrbc / (KdibpgRBCvHK * KiglcRBCvHK) + glcRBC / Vrbc * g16p2RBC / Vrbc / (Kdig16p2RBCvHK * KiglcRBCvHK) + glcRBC / Vrbc * g6pRBC / Vrbc / (Kdig6pRBCvHK * KiglcRBCvHK) + g6pRBC / Vrbc / Kig6pRBCvHK + glcRBC / Vrbc * gshRBC / Vrbc / (KdigshRBCvHK * KiglcRBCvHK) + mgadpRBC / Vrbc / KimgadpRBCvHK + g6pRBC / Vrbc * mgadpRBC / Vrbc / (Kig6pRBCvHK * KmmgadpRBCvHK) + mgatpRBC / Vrbc / KimgatpRBCvHK + glcRBC / Vrbc * mgatpRBC / Vrbc / (KiglcRBCvHK * KmmgatpRBCvHK))
v_RBCvLACTRANSPORT:
 {1.0}lacRBC = $pool
  TsRBC * vBld * (-KiRBCvLACTRANSPORT * lacEXT / vBld + KoRBCvLACTRANSPORT * lacRBC / Vrbc)
v_RBCvLDH:
 {1.0}nadhRBC + {1.0}pyrRBC = {1.0}lacRBC + {1.0}nadRBC
  TsRBC * Vrbc * VRBCvLDH * (-(KcatrRBCvLDH * lacRBC / Vrbc * nadRBC / Vrbc / (KinadRBCvLDH * KmapplacRBCvLDH)) + KcatfRBCvLDH * nadhRBC / Vrbc * pyrRBC / Vrbc / (KinadhRBCvLDH * KmapppyrRBCvLDH)) / (nadRBC / Vrbc / KinadRBCvLDH + lacRBC / Vrbc * nadRBC / Vrbc / (KinadRBCvLDH * KmapplacRBCvLDH) + nadhRBC / Vrbc / KinadhRBCvLDH + KmnadRBCvLDH * lacRBC / Vrbc * nadhRBC / Vrbc / (KinadRBCvLDH * KinadhRBCvLDH * KmapplacRBCvLDH) + KmnadhRBCvLDH * nadRBC / Vrbc * pyrRBC / Vrbc / (KinadRBCvLDH * KinadhRBCvLDH * KmapppyrRBCvLDH) + lacRBC / Vrbc * nadRBC / Vrbc * pyrRBC / Vrbc / (KinadRBCvLDH * KiapppyrRBCvLDH * KmapplacRBCvLDH) + nadhRBC / Vrbc * pyrRBC / Vrbc / (KinadhRBCvLDH * KmapppyrRBCvLDH) + lacRBC / Vrbc * nadhRBC / Vrbc * pyrRBC / Vrbc / (KiapplacRBCvLDH * KinadhRBCvLDH * KmapppyrRBCvLDH) + (1 + pyrRBC / Vrbc / KidpyrRBCvLDH) * (1 + KmnadRBCvLDH * lacRBC / Vrbc / (KinadRBCvLDH * KmapplacRBCvLDH) + KmnadhRBCvLDH * pyrRBC / Vrbc / (KinadhRBCvLDH * KmapppyrRBCvLDH)))
v_RBCvLDHP:
 {1.0}nadphRBC + {1.0}pyrRBC = {1.0}lacRBC + {1.0}nadpRBC
  TsRBC * Vrbc * ((-(KrRBCvLDHP * lacRBC / Vrbc * nadpRBC / Vrbc / KmlacRBCvLDHP) + KfRBCvLDHP * nadphRBC / Vrbc * pyrRBC / Vrbc / KmpyrRBCvLDHP) / (1 + lacRBC / Vrbc / KmlacRBCvLDHP + pyrRBC / Vrbc / KmpyrRBCvLDHP))
v_RBCvMGADP:
 {1.0}mgRBC + {1.0}adpRBC = {1.0}mgadpRBC
  TsRBC * Vrbc * (KaappRBCvMGADP * adpRBC / Vrbc * mgRBC / Vrbc - KdRBCvMGADP * mgadpRBC / Vrbc)
v_RBCvMGATP:
 {1.0}mgRBC + {1.0}atpRBC = {1.0}mgatpRBC
  TsRBC * Vrbc * (KaappRBCvMGATP * atpRBC / Vrbc * mgRBC / Vrbc - KdRBCvMGATP * mgatpRBC / Vrbc)
v_RBCvMGB13PG:
 {1.0}mgRBC + {1.0}b13pgRBC = {1.0}mgb13pgRBC
  TsRBC * Vrbc * (KaappRBCvMGB13PG * b13pgRBC / Vrbc * mgRBC / Vrbc - KdRBCvMGB13PG * mgb13pgRBC / Vrbc)
v_RBCvMGB23PG:
 {1.0}mgRBC + {1.0}b23pgRBC = {1.0}mgb23pgRBC
  TsRBC * Vrbc * (KaappRBCvMGB23PG * b23pgRBC / Vrbc * mgRBC / Vrbc - KdRBCvMGB23PG * mgb23pgRBC / Vrbc)
v_RBCvMGF16P2:
 {1.0}mgRBC + {1.0}f16p2RBC = {1.0}mgf16p2RBC
  TsRBC * Vrbc * (KaappRBCvMGF16P2 * f16p2RBC / Vrbc * mgRBC / Vrbc - KdRBCvMGF16P2 * mgf16p2RBC / Vrbc)
v_RBCvMGG16P2:
 {1.0}mgRBC + {1.0}g16p2RBC = {1.0}mgg16p2RBC
  TsRBC * Vrbc * (KaappRBCvMGG16P2 * g16p2RBC / Vrbc * mgRBC / Vrbc - KdRBCvMGG16P2 * mgg16p2RBC / Vrbc)
v_RBCvMGPHOS:
 {1.0}mgRBC + {1.0}phosRBC = {1.0}mgphosRBC
  TsRBC * Vrbc * (-KdRBCvMGPHOS * mgphosRBC / Vrbc + KaappRBCvMGPHOS * mgRBC / Vrbc * phosRBC / Vrbc)
v_RBCvOX:
 {2.0}gshRBC = {1.0}gssgRBC
  TsRBC * Vrbc * (KRBCvOX * gshRBC / Vrbc)
v_RBCvOXNADH:
 {1.0}nadhRBC = {1.0}nadRBC
  TsRBC * Vrbc * (KRBCvOXNADH * nadhRBC / Vrbc)
v_RBCvP6GDH:
 {1.0}p6gRBC + {1.0}nadpRBC = {1.0}nadphRBC + {1.0}ru5pRBC
  TsRBC * Vrbc * (VRBCvP6GDH * (K1RBCvP6GDH * K3RBCvP6GDH * K5RBCvP6GDH * K7RBCvP6GDH * K9RBCvP6GDH * K11RBCvP6GDH * nadpRBC / Vrbc * p6gRBC / Vrbc - co2RBC / Vrbc * K2RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K10RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * ru5pRBC / Vrbc) / (co2RBC / Vrbc * K2RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K11RBCvP6GDH + K2RBCvP6GDH * (K4RBCvP6GDH * K6RBCvP6GDH + K4RBCvP6GDH * K7RBCvP6GDH + K5RBCvP6GDH * K7RBCvP6GDH) * K9RBCvP6GDH * K11RBCvP6GDH + co2RBC / Vrbc * K1RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K11RBCvP6GDH * nadpRBC / Vrbc + K1RBCvP6GDH * (K4RBCvP6GDH * K6RBCvP6GDH + K4RBCvP6GDH * K7RBCvP6GDH + K5RBCvP6GDH * K7RBCvP6GDH) * K9RBCvP6GDH * K11RBCvP6GDH * nadpRBC / Vrbc + co2RBC / Vrbc * K2RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc + K2RBCvP6GDH * (K4RBCvP6GDH * K6RBCvP6GDH + K4RBCvP6GDH * K7RBCvP6GDH + K5RBCvP6GDH * K7RBCvP6GDH) * K9RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc + K3RBCvP6GDH * K5RBCvP6GDH * K7RBCvP6GDH * K9RBCvP6GDH * K11RBCvP6GDH * p6gRBC / Vrbc + co2RBC / Vrbc * K1RBCvP6GDH * K3RBCvP6GDH * (K5RBCvP6GDH + K6RBCvP6GDH) * K8RBCvP6GDH * K11RBCvP6GDH * nadpRBC / Vrbc * p6gRBC / Vrbc + K1RBCvP6GDH * K3RBCvP6GDH * (K5RBCvP6GDH * K7RBCvP6GDH * K9RBCvP6GDH + K5RBCvP6GDH * K7RBCvP6GDH * K11RBCvP6GDH + K5RBCvP6GDH * K9RBCvP6GDH * K11RBCvP6GDH + K6RBCvP6GDH * K9RBCvP6GDH * K11RBCvP6GDH + K7RBCvP6GDH * K9RBCvP6GDH * K11RBCvP6GDH) * nadpRBC / Vrbc * p6gRBC / Vrbc + K3RBCvP6GDH * K5RBCvP6GDH * K7RBCvP6GDH * K9RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * p6gRBC / Vrbc + co2RBC / Vrbc * K2RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K10RBCvP6GDH * ru5pRBC / Vrbc + co2RBC / Vrbc * K1RBCvP6GDH * K4RBCvP6GDH * K6RBCvP6GDH * K8RBCvP6GDH * K10RBCvP6GDH * nadpRBC / Vrbc * ru5pRBC / Vrbc + K2RBCvP6GDH * (K4RBCvP6GDH * K6RBCvP6GDH + K4RBCvP6GDH * K7RBCvP6GDH + K5RBCvP6GDH * K7RBCvP6GDH) * K10RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * ru5pRBC / Vrbc + co2RBC / Vrbc * (K2RBCvP6GDH * K4RBCvP6GDH + K2RBCvP6GDH * K5RBCvP6GDH + K2RBCvP6GDH * K6RBCvP6GDH + K4RBCvP6GDH * K6RBCvP6GDH) * K8RBCvP6GDH * K10RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * ru5pRBC / Vrbc + K1RBCvP6GDH * K3RBCvP6GDH * K5RBCvP6GDH * K7RBCvP6GDH * K10RBCvP6GDH * nadpRBC / Vrbc * p6gRBC / Vrbc * ru5pRBC / Vrbc + co2RBC / Vrbc * K1RBCvP6GDH * K3RBCvP6GDH * (K5RBCvP6GDH + K6RBCvP6GDH) * K8RBCvP6GDH * K10RBCvP6GDH * nadpRBC / Vrbc * p6gRBC / Vrbc * ru5pRBC / Vrbc + K3RBCvP6GDH * K5RBCvP6GDH * K7RBCvP6GDH * K10RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * p6gRBC / Vrbc * ru5pRBC / Vrbc + co2RBC / Vrbc * K3RBCvP6GDH * (K5RBCvP6GDH + K6RBCvP6GDH) * K8RBCvP6GDH * K10RBCvP6GDH * K12RBCvP6GDH * nadphRBC / Vrbc * p6gRBC / Vrbc * ru5pRBC / Vrbc))
v_RBCvPFK:
 {1.0}mgatpRBC + {1.0}f6pRBC = {1.0}f16p2RBC + {1.0}mgadpRBC
  TsRBC * Vrbc * (VRBCvPFK * (-(f16p2RBC / Vrbc * KcatfRBCvPFK * mgadpRBC / Vrbc / (Kmf16p2RBCvPFK * KmmgadpRBCvPFK)) + f6pRBC / Vrbc * KcatfRBCvPFK * mgatpRBC / Vrbc / (Kmf6pRBCvPFK * KmmgatpRBCvPFK)) / ((1 + f16p2RBC / Vrbc / Kmf16p2RBCvPFK + f6pRBC / Vrbc / Kmf6pRBCvPFK + mgadpRBC / Vrbc / KmmgadpRBCvPFK + f16p2RBC / Vrbc * mgadpRBC / Vrbc / (Kmf16p2RBCvPFK * KmmgadpRBCvPFK) + mgatpRBC / Vrbc / KmmgatpRBCvPFK + f6pRBC / Vrbc * mgatpRBC / Vrbc / (Kmf6pRBCvPFK * KmmgatpRBCvPFK)) * (1 + pow(hRBC / KaRBCvPFK, nRBCvPFK) * pow(1 + atpRBC / Vrbc / KtatpRBCvPFK, 4) * pow(1 + b23pgRBC / Vrbc / Ktb23pgRBCvPFK, 4) * pow(1 + mgRBC / Vrbc / KtmgRBCvPFK, 4) / (pow(1 + f16p2RBC / Vrbc / Kmf16p2RBCvPFK + f6pRBC / Vrbc / Kmf6pRBCvPFK, 4) * pow(1 + ampRBC / Vrbc / KrampRBCvPFK, 4) * pow(1 + g16p2RBC / Vrbc / Krg16p2RBCvPFK, 4) * pow(1 + phosRBC / Vrbc / KrphosRBCvPFK, 4)))))
v_RBCvPGI:
 {1.0}g6pRBC = {1.0}f6pRBC
  TsRBC * Vrbc * (ERBCvPGI * (-(f6pRBC / Vrbc * KcatrRBCvPGI / Kmf6pRBCvPGI) + g6pRBC / Vrbc * KcatfRBCvPGI / Kmg6pRBCvPGI) / (1 + f6pRBC / Vrbc / Kmf6pRBCvPGI + g6pRBC / Vrbc / Kmg6pRBCvPGI))
v_RBCvPGK:
 {1.0}b13pgRBC + {1.0}mgadpRBC = {1.0}mgatpRBC + {1.0}p3gRBC
  TsRBC * Vrbc * (VRBCvPGK * (b13pgRBC / Vrbc * KcatfRBCvPGK * mgadpRBC / Vrbc / (KimgadpRBCvPGK * Kmb13pgRBCvPGK) - KcatrRBCvPGK * mgatpRBC / Vrbc * p3gRBC / Vrbc / (KimgatpRBCvPGK * Kmp3gRBCvPGK)) / (1 + b13pgRBC / Vrbc / Kib13pgRBCvPGK + mgadpRBC / Vrbc / KimgadpRBCvPGK + b13pgRBC / Vrbc * mgadpRBC / Vrbc / (KimgadpRBCvPGK * Kmb13pgRBCvPGK) + mgatpRBC / Vrbc / KimgatpRBCvPGK + p3gRBC / Vrbc / Kip3gRBCvPGK + mgatpRBC / Vrbc * p3gRBC / Vrbc / (KimgatpRBCvPGK * Kmp3gRBCvPGK)))
v_RBCvPGLHYDROLYSIS:
 {1.0}p6glRBC = {1.0}p6gRBC
  TsRBC * Vrbc * (KhydrolRBCvPGLHYDROLYSIS * p6glRBC / Vrbc + VRBCvPGLHYDROLYSIS * Klactonase3RBCvPGLHYDROLYSIS * p6glRBC / Vrbc / ((Klactonase2RBCvPGLHYDROLYSIS + Klactonase3RBCvPGLHYDROLYSIS) / Klactonase1RBCvPGLHYDROLYSIS + p6glRBC / Vrbc))
v_RBCvPGM:
 {1.0}p3gRBC = {1.0}p2gRBC
  TsRBC * Vrbc * (VRBCvPGM * (-(KcatrRBCvPGM * p2gRBC / Vrbc / Kmp2gRBCvPGM) + KcatfRBCvPGM * p3gRBC / Vrbc / Kmp3gRBCvPGM) / (1 + p2gRBC / Vrbc / Kmp2gRBCvPGM + p3gRBC / Vrbc / Kmp3gRBCvPGM))
v_RBCvPHOSTRANSPORT:
 {1.0}phosRBC = $pool
  TsRBC * vBld * (-KiRBCvPHOSTRANSPORT * phosEXT / vBld + KoRBCvPHOSTRANSPORT * phosRBC / Vrbc)
v_RBCvPK:
 {1.0}pepRBC + {1.0}mgadpRBC = {1.0}pyrRBC + {1.0}mgatpRBC
  TsRBC * Vrbc * (VRBCvPK * (KcatfRBCvPK * mgadpRBC / Vrbc * pepRBC / Vrbc / (KrmgadpRBCvPK * KrpepRBCvPK) - KcatrRBCvPK * mgatpRBC / Vrbc * pyrRBC / Vrbc / (KrmgatpRBCvPK * KrpyrRBCvPK)) / ((1 + mgadpRBC / Vrbc / KrmgadpRBCvPK + mgatpRBC / Vrbc / KrmgatpRBCvPK + pepRBC / Vrbc / KrpepRBCvPK + mgadpRBC / Vrbc * pepRBC / Vrbc / (KrmgadpRBCvPK * KrpepRBCvPK) + pyrRBC / Vrbc / KrpyrRBCvPK + mgatpRBC / Vrbc * pyrRBC / Vrbc / (KrmgatpRBCvPK * KrpyrRBCvPK)) * (1 + LRBCvPK)))
v_RBCvPYRTRANSPORT:
 {1.0}pyrRBC = $pool
  TsRBC * vBld * (-KiRBCvPYRTRANSPORT * pyrEXT / vBld + KoRBCvPYRTRANSPORT * pyrRBC / Vrbc)
v_RBCvR5PI:
 {1.0}ru5pRBC = {1.0}rib5pRBC
  TsRBC * Vrbc * (VRBCvR5PI * (-(K2RBCvR5PI * K4RBCvR5PI * rib5pRBC / Vrbc / (K2RBCvR5PI + K3RBCvR5PI)) + K1RBCvR5PI * K3RBCvR5PI * ru5pRBC / Vrbc / (K2RBCvR5PI + K3RBCvR5PI)) / (1 + K1RBCvR5PI * ru5pRBC / Vrbc / (K2RBCvR5PI + K3RBCvR5PI) + K4RBCvR5PI * xu5pRBC / Vrbc / (K2RBCvR5PI + K3RBCvR5PI)))
v_RBCvRu5PE:
 {1.0}ru5pRBC = {1.0}xu5pRBC
  TsRBC * Vrbc * (VRBCvRu5PE * (K1RBCvRu5PE * K3RBCvRu5PE * ru5pRBC / Vrbc / (K2RBCvRu5PE + K3RBCvRu5PE) - K2RBCvRu5PE * K4RBCvRu5PE * xu5pRBC / Vrbc / (K2RBCvRu5PE + K3RBCvRu5PE)) / (1 + K1RBCvRu5PE * ru5pRBC / Vrbc / (K2RBCvRu5PE + K3RBCvRu5PE) + K4RBCvRu5PE * xu5pRBC / Vrbc / (K2RBCvRu5PE + K3RBCvRu5PE)))
v_RBCvTA:
 {1.0}sed7pRBC + {1.0}gapRBC = {1.0}ery4pRBC + {1.0}f6pRBC
  TsRBC * Vrbc * (VRBCvTA * (-ery4pRBC / Vrbc * f6pRBC / Vrbc * K2RBCvTA * K4RBCvTA * K6RBCvTA * K8RBCvTA + gapRBC / Vrbc * K1RBCvTA * K3RBCvTA * K5RBCvTA * K7RBCvTA * sed7pRBC / Vrbc) / (gapRBC / Vrbc * (K2RBCvTA + K3RBCvTA) * K5RBCvTA * K7RBCvTA + ery4pRBC / Vrbc * K2RBCvTA * K4RBCvTA * (K6RBCvTA + K7RBCvTA) + f6pRBC / Vrbc * gapRBC / Vrbc * (K2RBCvTA + K3RBCvTA) * K5RBCvTA * K8RBCvTA + f6pRBC / Vrbc * (K2RBCvTA + K3RBCvTA) * K6RBCvTA * K8RBCvTA + ery4pRBC / Vrbc * f6pRBC / Vrbc * K4RBCvTA * (K2RBCvTA + K6RBCvTA) * K8RBCvTA + gapRBC / Vrbc * K1RBCvTA * K5RBCvTA * (K3RBCvTA + K7RBCvTA) * sed7pRBC / Vrbc + K1RBCvTA * K3RBCvTA * (K6RBCvTA + K7RBCvTA) * sed7pRBC / Vrbc + ery4pRBC / Vrbc * K1RBCvTA * K4RBCvTA * (K6RBCvTA + K7RBCvTA) * sed7pRBC / Vrbc))
v_RBCvTIM:
 {1.0}gapRBC = {1.0}dhapRBC
  TsRBC * Vrbc * (VRBCvTIM * (gapRBC / Vrbc * KcatfRBCvTIM / KmgapRBCvTIM - dhapRBC / Vrbc * KcatrRBCvTIM / KmdhapRBCvTIM) / (1 + gapRBC / Vrbc / KmgapRBCvTIM + dhapRBC / Vrbc / KmdhapRBCvTIM))
v_RBCvTK1:
 {1.0}xu5pRBC + {1.0}tkRBC = {1.0}tkxu5pRBC
  TsRBC * Vrbc * (K1RBCvTK1 * tkRBC / Vrbc * xu5pRBC / Vrbc - K2RBCvTK1 * tkxu5pRBC / Vrbc)
v_RBCvTK2:
 {1.0}tkxu5pRBC = {1.0}tkgRBC + {1.0}gapRBC
  TsRBC * Vrbc * (K3RBCvTK2 * tkxu5pRBC / Vrbc - K4RBCvTK2 * tkgRBC / Vrbc * gapRBC / Vrbc)
v_RBCvTK3:
 {1.0}tkgRBC + {1.0}rib5pRBC = {1.0}tkgrib5pRBC
  TsRBC * Vrbc * (K5RBCvTK3 * tkgRBC / Vrbc * rib5pRBC / Vrbc - K6RBCvTK3 * tkgrib5pRBC / Vrbc)
v_RBCvTK4:
 {1.0}tkgrib5pRBC = {1.0}tkRBC + {1.0}sed7pRBC
  TsRBC * Vrbc * (K7RBCvTK4 * tkgrib5pRBC / Vrbc - K8RBCvTK4 * tkRBC / Vrbc * sed7pRBC / Vrbc)
v_RBCvTK5:
 {1.0}tkgRBC + {1.0}ery4pRBC = {1.0}tkgery4pRBC
  TsRBC * Vrbc * (K9RBCvTK5 * tkgRBC / Vrbc * ery4pRBC / Vrbc - K10RBCvTK5 * tkgery4pRBC / Vrbc)
v_RBCvTK6:
 {1.0}tkgery4pRBC = {1.0}tkRBC + {1.0}f6pRBC
  TsRBC * Vrbc * (K11RBCvTK6 * tkgery4pRBC / Vrbc - K12RBCvTK6 * tkRBC / Vrbc * f6pRBC / Vrbc)


#parameters
AlphaCellWaterFraction = 0.7
CRBCvMGADP = 0.00062
CRBCvMGATP = 8.4e-05
CRBCvMGB13PG = 0.0032
CRBCvMGB23PG = 0.0032
CRBCvMGF16P2 = 0.0083
CRBCvMGG16P2 = 0.0083
CRBCvMGPHOS = 40800.0
ConcADPRBC = 0.00031
ConcAMPRBC = 3e-05
ConcATPRBC = 0.0021
ConcB13PGRBC = 7e-07
ConcB23PGRBC = 0.0067
ConcBPGSPB13PGRBC = 0.0
ConcBPGSPB23PGRBC = 0.0
ConcBPGSPPP2GRBC = 0.0
ConcBPGSPPP3GRBC = 0.0
ConcBPGSPPRBC = 0.0
ConcBPGSPRBC = 3.8e-06
ConcBpgsppphosRBC = 0.0
ConcCo2RBC = 0.0012
ConcDhapRBC = 1.9e-05
ConcEry4PRBC = 1e-05
ConcF16P2RBC = 2.7e-06
ConcF6PRBC = 1.3e-05
ConcG16P2RBC = 0.000122
ConcG6PRBC = 4e-05
ConcGSHRBC = 0.0032
ConcGapRBC = 5.7e-06
ConcGlcRBC = 0.005
ConcGssgRBC = 9e-08
ConcHbADPRBC = 0.0
ConcHbATPRBC = 0.0
ConcHbB13PGRBC = 0.0
ConcHbB23PGRBC = 0.0
ConcHbRBC = 0.007
ConcHbmgatpRBC = 0.0
ConcLacEXT = 0.00182
ConcLacRBC = 0.0014
ConcMgADPRBC = 0.0
ConcMgB13PGRBC = 0.0
ConcMgB23PGRBC = 0.0
ConcMgF16P2RBC = 0.0
ConcMgRBC = 0.003
ConcMgatpRBC = 0.0
ConcMgg16p2RBC = 0.0
ConcMgphosRBC = 0.0
ConcNADPRBC = 1.25e-07
ConcNADRBC = 6e-05
ConcNadhRBC = 1.4e-07
ConcNadphRBC = 6.4e-05
ConcP2GARBC = 1e-05
ConcP3GARBC = 6.4e-05
ConcP6GLRBC = 1.4e-10
ConcP6GRBC = 1.4e-07
ConcPEPRBC = 2.3e-05
ConcPhosEXT = 0.00192
ConcPhosRBC = 0.001
ConcPyrEXT = 8.5e-05
ConcPyrRBC = 6e-05
ConcRib5PRBC = 1e-05
ConcRu5PRBC = 1e-05
ConcSed7PRBC = 1e-05
ConcTKGEry4PRBC = 0.0
ConcTKGRBC = 0.0
ConcTKGRib5PRBC = 0.0
ConcTKRBC = 3.3e-07
ConcTkxu5pRBC = 0.0
ConcXu5PRBC = 1e-06
ERBCvHK = 2.5e-08
ERBCvPGI = 2.18e-07
K10RBCvBPGSP6 = 1979.0
K10RBCvG6PDH = 1400000000.0
K10RBCvGSSGR = 50000000.0
K10RBCvP6GDH = 225000.0
K10RBCvTK5 = 175.0
K11RBCvBPGSP6 = 0.01
K11RBCvGSSGR = 7000.0
K11RBCvP6GDH = 300.0
K11RBCvTK6 = 40.0
K12RBCvBPGSP7 = 1000.0
K12RBCvGSSGR = 100000000.0
K12RBCvP6GDH = 4950000.0
K12RBCvTK6 = 21300.0
K13RBCvBPGSP7 = 1800000.0
K14RBCvBPGSP8 = 1000000000.0
K15RBCvBPGSP8 = 610000.0
K16RBCvBPGSP9 = 0.19
K1RBCvAK = 4300.0
K1RBCvBPGSP1 = 80000000.0
K1RBCvG6PDH = 110000000.0
K1RBCvGSSGR = 85000000.0
K1RBCvP6GDH = 2400000.0
K1RBCvR5PI = 60900.0
K1RBCvRu5PE = 3910000.0
K1RBCvTA = 580000.0
K1RBCvTK1 = 216000.0
K2RBCvAK = 1400.0
K2RBCvBPGSP1 = 400.0
K2RBCvG6PDH = 870.0
K2RBCvGSSGR = 510.0
K2RBCvP6GDH = 410.0
K2RBCvR5PI = 33.3
K2RBCvRu5PE = 438.0
K2RBCvTA = 45.3
K2RBCvTK1 = 38.0
K3RBCvBPGSP2 = 9.9
K3RBCvG6PDH = 26000000.0
K3RBCvGSSGR = 1000000000.0
K3RBCvP6GDH = 2000000000.0
K3RBCvR5PI = 14.2
K3RBCvRu5PE = 305.0
K3RBCvTA = 16.3
K3RBCvTK2 = 34.0
K4RBCvBPGSP3 = 185000000.0
K4RBCvG6PDH = 300.0
K4RBCvGSSGR = 72000.0
K4RBCvP6GDH = 26000.0
K4RBCvR5PI = 21600.0
K4RBCvRu5PE = 1490000.0
K4RBCvTA = 1010000.0
K4RBCvTK2 = 156000.0
K5RBCvBPGSP3 = 100000000.0
K5RBCvG6PDH = 750.0
K5RBCvGSSGR = 810.0
K5RBCvP6GDH = 48.0
K5RBCvTA = 490000.0
K5RBCvTK3 = 329000.0
K6RBCvBPGSP4 = 1000.0
K6RBCvG6PDH = 2000.0
K6RBCvGSSGR = 1000.0
K6RBCvP6GDH = 30.0
K6RBCvTA = 60.0
K6RBCvTK3 = 175.0
K7RBCvBPGSP4 = 1000.0
K7RBCvG6PDH = 220000.0
K7RBCvGSSGR = 1000000.0
K7RBCvP6GDH = 630.0
K7RBCvTA = 17.0
K7RBCvTK4 = 40.0
K8RBCvBPGSP5 = 10000.0
K8RBCvG6PDH = 1100000000.0
K8RBCvGSSGR = 50000000.0
K8RBCvP6GDH = 36000.0
K8RBCvTA = 79000.0
K8RBCvTK4 = 44800.0
K9RBCvBPGSP5 = 0.55
K9RBCvG6PDH = 10000.0
K9RBCvGSSGR = 1000000.0
K9RBCvP6GDH = 800.0
K9RBCvTK5 = 2240000.0
KRBCvATPASE = 0.000585
KRBCvMGADP = 3290.0
KRBCvMGATP = 43200.0
KRBCvMGB13PG = 7410.0
KRBCvMGB23PG = 7410.0
KRBCvMGF16P2 = 363.0
KRBCvMGG16P2 = 363.0
KRBCvOX = 3.4e-05
KRBCvOXNADH = 0.0163
KaRBCvHBADP = 300000.0
KaRBCvHBATP = 432000.0
KaRBCvHBB13PG = 380000.0
KaRBCvHBBPG = 300000.0
KaRBCvHBMGATP = 46800.0
KaRBCvMGADP = 2760000.0
KaRBCvMGATP = 31200000.0
KaRBCvMGB13PG = 228000.0
KaRBCvMGB23PG = 804000.0
KaRBCvMGF16P2 = 480000.0
KaRBCvMGG16P2 = 480000.0
KaRBCvMGPHOS = 40800.0
KaRBCvPFK = 8.91251e-08
KahbRBC = 2.51189e-07
KcatfRBCvALD = 68.0
KcatfRBCvENO = 190.0
KcatfRBCvGAPDH = 695.304
KcatfRBCvHK = 299.16
KcatfRBCvLDH = 458.0
KcatfRBCvPFK = 822.0
KcatfRBCvPGI = 1470.0
KcatfRBCvPGK = 2290.0
KcatfRBCvPGM = 795.0
KcatfRBCvPK = 1386.0
KcatfRBCvTIM = 14560.0
KcatrRBCvALD = 234.0
KcatrRBCvENO = 50.0
KcatrRBCvGAPDH = 512.487
KcatrRBCvHK = 1.92792
KcatrRBCvLDH = 115.0
KcatrRBCvPFK = 36.0
KcatrRBCvPGI = 1760.0
KcatrRBCvPGK = 917.0
KcatrRBCvPGM = 714.0
KcatrRBCvPK = 3.26
KcatrRBCvTIM = 1280.0
KdRBCvHBADP = 1200.0
KdRBCvHBATP = 1200.0
KdRBCvHBB13PG = 1200.0
KdRBCvHBBPG = 1200.0
KdRBCvHBMGATP = 1200.0
KdRBCvMGADP = 1200.0
KdRBCvMGATP = 1200.0
KdRBCvMGB13PG = 1200.0
KdRBCvMGB23PG = 1200.0
KdRBCvMGF16P2 = 1200.0
KdRBCvMGG16P2 = 1200.0
KdRBCvMGPHOS = 1200.0
KdibpgRBCvHK = 0.004
Kdig16p2RBCvHK = 3e-05
Kdig6pRBCvHK = 1e-05
KdigshRBCvHK = 0.003
KfRBCvLDHP = 0.00346
Kh2bpgRBCvMGB13PG = 4270000.0
Kh2bpgRBCvMGB23PG = 4270000.0
Kh2fRBCvMGF16P2 = 1120000.0
Kh2fRBCvMGG16P2 = 1120000.0
KhadpRBCvAK = 5420000.0
KhadpRBCvMGADP = 5420000.0
KhampRBCvAK = 3090000.0
KhatpRBCvMGATP = 9070000.0
KhbpgRBCvMGB13PG = 162000000.0
KhbpgRBCvMGB23PG = 162000000.0
KhfRBCvMGF16P2 = 7560000.0
KhfRBCvMGG16P2 = 7560000.0
KhphosRBCvMGPHOS = 5680000.0
KhydrolRBCvPGLHYDROLYSIS = 0.00071
KiRBCvLACTRANSPORT = 0.0036
KiRBCvPHOSTRANSPORT = 0.00056
KiRBCvPYRTRANSPORT = 0.018
Kib13pgRBCvGAPDH = 4.55544e-21
Kib13pgRBCvPGK = 1.6e-06
Kib23pgRBCvALD = 0.0015
Kidb13pgRBCvGAPDH = 1e-06
KidgapRBCvGAPDH = 3.1e-05
KidhapRBCvALD = 1.1e-05
KidpyrRBCvLDH = 0.000101
Kif16p2RBCvALD = 1.98e-05
Kig6pRBCvHK = 4.7e-05
KigapRBCvGAPDH = 4.76523e-19
KiglcRBCvHK = 4.7e-05
KilacRBCvLDH = 0.00733
KimgRBCvENO = 4.6e-05
KimgadpRBCvHK = 0.001
KimgadpRBCvPGK = 8e-05
KimgatpRBCvHK = 0.001
KimgatpRBCvPGK = 0.000186
KinadRBCvGAPDH = 4.5e-05
KinadRBCvLDH = 0.000503
KinadhRBCvGAPDH = 1e-05
KinadhRBCvLDH = 2.45e-06
Kip2gRBCvENO = 0.00014
Kip3gRBCvPGK = 0.000205
KipepRBCvENO = 0.0001105
KiphosRBCvGAPDH = 0.00316
KipyrRBCvLDH = 0.000228
KkadpRBCvAK = 4.8
KkadpRBCvMGADP = 4.8
KkampRBCvAK = 1.8
KkatpRBCvMGATP = 14.0
KkbpgRBCvMGB13PG = 85.1
KkbpgRBCvMGB23PG = 85.1
KkfRBCvMGF16P2 = 10.7
KkfRBCvMGG16P2 = 10.7
KkhbpgRBCvMGB13PG = 8.9
KkhbpgRBCvMGB23PG = 8.9
KkhfRBCvMGF16P2 = 3.3
KkhfRBCvMGG16P2 = 3.3
KkphosRBCvMGPHOS = 3.0
Klactonase1RBCvPGLHYDROLYSIS = 13000000.0
Klactonase2RBCvPGLHYDROLYSIS = 1000.0
Klactonase3RBCvPGLHYDROLYSIS = 29.0
Kmb13pgRBCvGAPDH = 6.71e-07
Kmb13pgRBCvPGK = 2e-06
KmdhapRBCvALD = 3.5e-05
KmdhapRBCvTIM = 0.0001624
Kmf16p2RBCvALD = 7.1e-06
Kmf16p2RBCvPFK = 0.0005
Kmf6pRBCvPFK = 7.5e-05
Kmf6pRBCvPGI = 7.1e-05
Kmg6pRBCvPGI = 0.000181
KmgapRBCvALD = 0.000189
KmgapRBCvGAPDH = 9.5e-05
KmgapRBCvTIM = 0.000446
KmghadpRBCvMGADP = 107.0
KmghatpRBCvMGATP = 748.0
KmghbpgRBCvMGB13PG = 513.0
KmghbpgRBCvMGB23PG = 513.0
KmghfRBCvMGF16P2 = 89.0
KmghfRBCvMGG16P2 = 89.0
KmlacRBCvLDH = 0.00107
KmlacRBCvLDHP = 0.000414
KmmgRBCvENO = 4.6e-05
KmmgadpRBCvHK = 0.001
KmmgadpRBCvPFK = 0.00054
KmmgadpRBCvPGK = 0.0001
KmmgatpRBCvHK = 0.001
KmmgatpRBCvPFK = 6.8e-05
KmmgatpRBCvPGK = 0.001
KmnadRBCvGAPDH = 4.5e-05
KmnadRBCvLDH = 0.000107
KmnadhRBCvGAPDH = 3.3e-06
KmnadhRBCvLDH = 8.44e-06
Kmp2gRBCvENO = 0.00014
Kmp2gRBCvPGM = 2.56e-05
Kmp3gRBCvPGK = 0.0011
Kmp3gRBCvPGM = 0.000168
KmpepRBCvENO = 0.0001105
KmphosRBCvGAPDH = 0.00316
KmpyrRBCvLDH = 0.000137
KmpyrRBCvLDHP = 0.000414
KrRBCvLDHP = 5.43e-07
KrampRBCvPFK = 0.0003
Krf16p2RBCvPK = 5e-06
Krg16p2RBCvPFK = 0.01
Krg16p2RBCvPK = 0.0001
KrmgadpRBCvPK = 0.000474
KrmgatpRBCvPK = 0.003
KrpepRBCvPK = 0.000225
KrphosRBCvPFK = 0.03
KrpyrRBCvPK = 0.002
KtatpRBCvPFK = 0.0001
KtatpRBCvPK = 0.00339
Ktb23pgRBCvPFK = 0.005
KtmgRBCvPFK = 0.004
RtvRBC = 0.69
TsRBC = 60.0
VRBCvALD = 3.7e-07
VRBCvENO = 2.2e-07
VRBCvG6PDH = 9.3e-08
VRBCvGAPDH = 7.66e-06
VRBCvGSSGR = 1.25e-07
VRBCvLDH = 3.43e-06
VRBCvP6GDH = 2.1e-06
VRBCvPFK = 1.1e-07
VRBCvPGK = 2.74e-06
VRBCvPGLHYDROLYSIS = 1.4e-05
VRBCvPGM = 4.1e-07
VRBCvPK = 8.7e-08
VRBCvR5PI = 1.42e-05
VRBCvRu5PE = 4.22e-06
VRBCvTA = 6.9e-07
VRBCvTIM = 1.14e-06
bpgxRBC = 3.8e-06
gssgxRBC = 0.00160009
hbxRBC = 0.007
kRBC = 0.15
nRBCvPFK = 5.0
nadpxRBC = 6.4125e-05
pHConversionFactor = 1.0
phRBC = 7.2
tkxRBC = 3.3e-07
xatpRBC = -0.00456
xg16p2RBC = 0.000122
xphosRBC = 0.000438
Vrbc = 90.0
vBld = 128.5


#initial values
adpRBC = 0.0279
ampRBC = 27/10000
atpRBC = 0.189
b13pgRBC = 63/1000000
b23pgRBC = 0.603
bpgspRBC = 0.00034199999999999996
bpgspb13pgRBC = 0
bpgspb23pgRBC = 0
bpgsppRBC = 0
bpgsppp2gRBC = 0
bpgsppp3gRBC = 0
bpgsppphosRBC = 0
dhapRBC = 0.0017100000000000001
ery4pRBC = 9/10000
f16p2RBC = 0.00024300000000000002
f6pRBC = 0.00117
g16p2RBC = 0.01098
g6pRBC = 9/2500
gapRBC = 0.000513
gshRBC = 0.28800000000000003
gssgRBC = 81/10000000
hbRBC = 0.63
hbadpRBC = 0
hbatpRBC = 0
hbb13pgRBC = 0
hbb23pgRBC = 0
hbmgatpRBC = 0
lacRBC = 0.126
mgRBC = 0.27
mgadpRBC = 0
mgatpRBC = 0
mgb13pgRBC = 0
mgb23pgRBC = 0
mgf16p2RBC = 0
mgg16p2RBC = 0
mgphosRBC = 0
nadRBC = 27/5000
nadhRBC = 0.000012599999999999998
nadpRBC = 0.000011249999999999999
nadphRBC = 0.00576
p2gRBC = 9/10000
p3gRBC = 0.00576
p6gRBC = 0.000012599999999999998
p6glRBC = 1.2599999999999999*^-8
pepRBC = 0.00207
phosRBC = 0.09
pyrRBC = 27/5000
rib5pRBC = 9/10000
ru5pRBC = 9/10000
sed7pRBC = 9/10000
tkRBC = 0.0000297
tkgRBC = 0
tkgery4pRBC = 0
tkgrib5pRBC = 0
tkxu5pRBC = 0
xu5pRBC = 9/100000


#assignment rules
!F KiappgapRBCvGAPDH = KigapRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
!F K13appRBCvBPGSP7 = K13RBCvBPGSP7 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
!F KaappRBCvMGB23PG = KaRBCvMGB23PG * CRBCvMGB23PG * (KRBCvMGB23PG + KhbpgRBCvMGB23PG * KmghbpgRBCvMGB23PG * hRBC) / (1.0 + KhbpgRBCvMGB23PG * hRBC + Kh2bpgRBCvMGB23PG * KhbpgRBCvMGB23PG * hRBC * hRBC + kRBC * KkbpgRBCvMGB23PG + kRBC * KhbpgRBCvMGB23PG * KkhbpgRBCvMGB23PG * hRBC)
!F KaappRBCvMGF16P2 = KaRBCvMGF16P2 * CRBCvMGF16P2 * (KRBCvMGF16P2 + KhfRBCvMGF16P2 * KmghfRBCvMGF16P2 * hRBC) / (1.0 + KhfRBCvMGF16P2 * hRBC + Kh2fRBCvMGF16P2 * KhfRBCvMGF16P2 * hRBC * hRBC + kRBC * KkfRBCvMGF16P2 + kRBC * KhfRBCvMGF16P2 * KkhfRBCvMGF16P2 * hRBC)
!F K6appRBCvBPGSP4 = K6RBCvBPGSP4 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
!F K1appRBCvAK = K1RBCvAK * (1.0 + pHConversionFactor * pow(10.0, -phRBC) * KhadpRBCvAK + kRBC * KkadpRBCvAK)
!F KaappRBCvMGATP = KaRBCvMGATP * CRBCvMGATP * ((KRBCvMGATP + KhatpRBCvMGATP * KmghatpRBCvMGATP * hRBC) / (1.0 + KhatpRBCvMGATP / pHConversionFactor / pow(10.0, phRBC) + kRBC * KkatpRBCvMGATP))
!F KaappRBCvMGG16P2 = KaRBCvMGG16P2 * CRBCvMGG16P2 * (KRBCvMGG16P2 + KhfRBCvMGG16P2 * KmghfRBCvMGG16P2 * hRBC) / (1.0 + KhfRBCvMGG16P2 * hRBC + Kh2fRBCvMGG16P2 * KhfRBCvMGG16P2 * hRBC * hRBC + kRBC * KkfRBCvMGG16P2 + kRBC * KhfRBCvMGG16P2 * KkhfRBCvMGG16P2 * hRBC)
!F K2appRBCvAK = K2RBCvAK * (1.0 + hRBC * KhampRBCvAK + kRBC * KkampRBCvAK)
!F KaappRBCvMGB13PG = KaRBCvMGB13PG * CRBCvMGB13PG * (KRBCvMGB13PG + KhbpgRBCvMGB13PG * KmghbpgRBCvMGB13PG * hRBC) / (1.0 + KhbpgRBCvMGB13PG * hRBC + Kh2bpgRBCvMGB13PG * KhbpgRBCvMGB13PG * hRBC * hRBC + kRBC * KkbpgRBCvMGB13PG + kRBC * KhbpgRBCvMGB13PG * KkhbpgRBCvMGB13PG * hRBC)
!F KaappRBCvMGPHOS = KaRBCvMGPHOS * (1.0 + pHConversionFactor * pow(10.0, -7.2) * KhphosRBCvMGPHOS + kRBC * KkphosRBCvMGPHOS) / (1.0 + hRBC * KhphosRBCvMGPHOS + kRBC * KkphosRBCvMGPHOS)
!F KoRBCvPHOSTRANSPORT = KiRBCvPHOSTRANSPORT / KeqRBCvPHOSTRANSPORT
!F KaappRBCvMGADP = KaRBCvMGADP * CRBCvMGADP * ((KRBCvMGADP + KhadpRBCvMGADP * KmghadpRBCvMGADP * hRBC) / (1.0 + KhadpRBCvMGADP * hRBC + kRBC * KkadpRBCvMGADP))
!F KiapppyrRBCvLDH = KipyrRBCvLDH * ((1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)))
!F KcatrappRBCvGAPDH = KcatrRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
!F KmapplacRBCvLDH = KmlacRBCvLDH * ((1.0 + pow(10.0, -phRBC) / pow(10.0, -6.8)) / (1.0 + pow(10.0, -7.2) / pow(10.0, -6.8)))
!F KeqRBCvLACTRANSPORT = (1.0 + pow(10.0, phRBC - 3.73)) / (1.0 + pow(10.0, phRBC - 3.73) / RtvRBC)
!F KmappnadhRBCvGAPDH = KmnadhRBCvGAPDH * (pow(10.0, -7.2) / pow(10.0, -phRBC))
!F KmapppyrRBCvLDH = KmpyrRBCvLDH * ((1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)))
!F KcatrappRBCvHK = KcatrRBCvHK / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.02) + pow(10.0, -9.55) / pow(10.0, -phRBC))
!F HbpHRBC = (1.0 + 2.0 * KahbRBC / (pHConversionFactor * pow(10.0, -7.2)) + pow(KahbRBC / (pHConversionFactor * pow(10.0, -7.2)), 2.0)) / (1.0 + 2.0 * KahbRBC / hRBC + KahbRBC / hRBC * KahbRBC / hRBC)
!F KeqRBCvPHOSTRANSPORT = (1.0 + pow(10.0, phRBC - 6.75)) / (1.0 / RtvRBC + pow(10.0, phRBC - 6.75) / (RtvRBC * RtvRBC))
!F KcatfappRBCvHK = KcatfRBCvHK / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.02) + pow(10.0, -9.55) / pow(10.0, -phRBC))
!F KaappRBCvHBB13PG = KaRBCvHBB13PG * HbpHRBC
!F KiapplacRBCvLDH = KilacRBCvLDH * ((1.0 + pow(10.0, -phRBC) / pow(10.0, -6.8)) / (1.0 + pow(10.0, -7.2) / pow(10.0, -6.8)))
!F K4appRBCvBPGSP3 = K4RBCvBPGSP3 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
!F KcatfappRBCvGAPDH = KcatfRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
!F atpRBCsum = atpRBC+adpRBC+ampRBC
!F KoRBCvPYRTRANSPORT = KiRBCvPYRTRANSPORT / RtvRBC
!F hRBC = pHConversionFactor / pow(10.0, phRBC)
!F KaappRBCvHBATP = KaRBCvHBATP * HbpHRBC
!F KaappRBCvHBADP = KaRBCvHBADP * HbpHRBC
!F KoRBCvLACTRANSPORT = KiRBCvLACTRANSPORT / KeqRBCvLACTRANSPORT
!F K3appRBCvBPGSP2 = K3RBCvBPGSP2 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
!F KaappRBCvHBBPG = KaRBCvHBBPG * HbpHRBC
!F KiappnadhRBCvGAPDH = KinadhRBCvGAPDH * (pow(10.0, -7.2) / pow(10.0, -phRBC))
!F KaappRBCvHBMGATP = KaRBCvHBMGATP * HbpHRBC
!F K1appRBCvBPGSP1 = K1RBCvBPGSP1 * ((1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)))
!F testsum = adpRBC+atpRBC+ampRBC+mgRBC +mgadpRBC + mgatpRBC +mgb13pgRBC + mgb23pgRBC + mgf16p2RBC +   mgg16p2RBC	+ mgphosRBC
!F Kiappb13pgRBCvGAPDH = Kib13pgRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
!F LRBCvPK = pHConversionFactor * pow(10.0, -6.8) / hRBC * pow(1.0 + atpRBC / Vrbc / KtatpRBCvPK, 4.0) / (pow(1.0 + f16p2RBC / Vrbc / Krf16p2RBCvPK + g16p2RBC / Vrbc / Krg16p2RBCvPK, 4.0) * pow(1.0 + pepRBC / Vrbc / KrpepRBCvPK + pyrRBC / Vrbc / KrpyrRBCvPK, 4.0))
!F lacEXT = ConcLacEXT*vBld
!F phosEXT = ConcPhosEXT*vBld
!F pyrEXT = ConcPyrEXT*vBld
!F glcRBC = ConcGlcRBC*Vrbc
!F co2RBC = ConcCo2RBC * Vrbc