(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
adpRBC[t], ampRBC[t], atpRBC[t], b13pgRBC[t], b23pgRBC[t], bpgspRBC[t], bpgspb13pgRBC[t], bpgspb23pgRBC[t], bpgsppRBC[t], bpgsppp2gRBC[t], bpgsppp3gRBC[t], bpgsppphosRBC[t], dhapRBC[t], ery4pRBC[t], f16p2RBC[t], f6pRBC[t], g16p2RBC[t], g6pRBC[t], gapRBC[t], glcRBC[t], gshRBC[t], gssgRBC[t], hbRBC[t], hbadpRBC[t], hbatpRBC[t], hbb13pgRBC[t], hbb23pgRBC[t], hbmgatpRBC[t], lacRBC[t], mgRBC[t], mgadpRBC[t], mgatpRBC[t], mgb13pgRBC[t], mgb23pgRBC[t], mgf16p2RBC[t], mgg16p2RBC[t], mgphosRBC[t], nadRBC[t], nadhRBC[t], nadpRBC[t], nadphRBC[t], p2gRBC[t], p3gRBC[t], p6gRBC[t], p6glRBC[t], pepRBC[t], phosRBC[t], pyrRBC[t], rib5pRBC[t], ru5pRBC[t], sed7pRBC[t], tkRBC[t], tkgRBC[t], tkgery4pRBC[t], tkgrib5pRBC[t], tkxu5pRBC[t], xu5pRBC[t]
};

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

rates = {
vRBCvAK, vRBCvALD, vRBCvATPASE, vRBCvBPGSP1, vRBCvBPGSP2, vRBCvBPGSP3, vRBCvBPGSP4, vRBCvBPGSP5, vRBCvBPGSP6, vRBCvBPGSP7, vRBCvBPGSP8, vRBCvBPGSP9, vRBCvENO, vRBCvG6PDH, vRBCvGAPDH, vRBCvGLCTRANSPORT, vRBCvGSSGR, vRBCvHBADP, vRBCvHBATP, vRBCvHBB13PG, vRBCvHBB23PG, vRBCvHBMGATP, vRBCvHK, vRBCvLACTRANSPORT, vRBCvLDH, vRBCvLDHP, vRBCvMGADP, vRBCvMGATP, vRBCvMGB13PG, vRBCvMGB23PG, vRBCvMGF16P2, vRBCvMGG16P2, vRBCvMGPHOS, vRBCvOX, vRBCvOXNADH, vRBCvP6GDH, vRBCvPFK, vRBCvPGI, vRBCvPGK, vRBCvPGLHYDROLYSIS, vRBCvPGM, vRBCvPHOSTRANSPORT, vRBCvPK, vRBCvPYRTRANSPORT, vRBCvR5PI, vRBCvRu5PE, vRBCvTA, vRBCvTIM, vRBCvTK1, vRBCvTK2, vRBCvTK3, vRBCvTK4, vRBCvTK5, vRBCvTK6
};

rateEquations = {
vRBCvAK -> TsRBC*Vrbc*((K1appRBCvAK*adpRBC[t]*mgadpRBC[t])/Vrbc^2 - (K2appRBCvAK*ampRBC[t]*mgatpRBC[t])/Vrbc^2), vRBCvALD -> (TsRBC*Vrbc*VRBCvALD*((KcatfRBCvALD*f16p2RBC[t])/(Kmf16p2RBCvALD*Vrbc) - (KcatrRBCvALD*dhapRBC[t]*gapRBC[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbc^2)))/(1 + dhapRBC[t]/(KidhapRBCvALD*Vrbc) + f16p2RBC[t]/(Kmf16p2RBCvALD*Vrbc) + (dhapRBC[t]*gapRBC[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbc^2) + (KmdhapRBCvALD*f16p2RBC[t]*gapRBC[t])/(KidhapRBCvALD*Kif16p2RBCvALD*KmgapRBCvALD*Vrbc^2) + (b23pgRBC[t]/Vrbc + mgb23pgRBC[t]/Vrbc)/Kib23pgRBCvALD + (KmdhapRBCvALD*gapRBC[t]*(1 + (b23pgRBC[t]/Vrbc + mgb23pgRBC[t]/Vrbc)/Kib23pgRBCvALD))/(KidhapRBCvALD*KmgapRBCvALD*Vrbc)), vRBCvATPASE -> KRBCvATPASE*TsRBC*mgatpRBC[t], vRBCvBPGSP1 -> TsRBC*Vrbc*(-((K2RBCvBPGSP1*bpgspb13pgRBC[t])/Vrbc) + (K1appRBCvBPGSP1*b13pgRBC[t]*bpgspRBC[t])/Vrbc^2), vRBCvBPGSP2 -> K3appRBCvBPGSP2*TsRBC*bpgspb13pgRBC[t], vRBCvBPGSP3 -> TsRBC*Vrbc*(-((K5RBCvBPGSP3*bpgsppp3gRBC[t])/Vrbc) + (K4appRBCvBPGSP3*bpgsppRBC[t]*p3gRBC[t])/Vrbc^2), vRBCvBPGSP4 -> TsRBC*Vrbc*(-((K7RBCvBPGSP4*bpgsppp2gRBC[t])/Vrbc) + (K6appRBCvBPGSP4*bpgsppRBC[t]*p2gRBC[t])/Vrbc^2), vRBCvBPGSP5 -> TsRBC*Vrbc*(-((K9RBCvBPGSP5*bpgspb23pgRBC[t])/Vrbc) + (K8RBCvBPGSP5*bpgsppp3gRBC[t])/Vrbc), vRBCvBPGSP6 -> TsRBC*Vrbc*(-((K11RBCvBPGSP6*bpgspb23pgRBC[t])/Vrbc) + (K10RBCvBPGSP6*bpgsppp2gRBC[t])/Vrbc), vRBCvBPGSP7 -> TsRBC*Vrbc*((K12RBCvBPGSP7*bpgspb23pgRBC[t])/Vrbc - (K13appRBCvBPGSP7*b23pgRBC[t]*bpgspRBC[t])/Vrbc^2), vRBCvBPGSP8 -> TsRBC*Vrbc*(-((K15RBCvBPGSP8*bpgsppphosRBC[t])/Vrbc) + (K14RBCvBPGSP8*bpgsppRBC[t]*phosRBC[t])/Vrbc^2), vRBCvBPGSP9 -> K16RBCvBPGSP9*TsRBC*bpgsppphosRBC[t], vRBCvENO -> (TsRBC*Vrbc*VRBCvENO*((KcatfRBCvENO*mgRBC[t]*p2gRBC[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbc^2) - (KcatrRBCvENO*mgRBC[t]*pepRBC[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbc^2)))/(1 + mgRBC[t]/(KimgRBCvENO*Vrbc) + p2gRBC[t]/(Kip2gRBCvENO*Vrbc) + (mgRBC[t]*p2gRBC[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbc^2) + pepRBC[t]/(KipepRBCvENO*Vrbc) + (mgRBC[t]*pepRBC[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbc^2)), vRBCvG6PDH -> (TsRBC*Vrbc*VRBCvG6PDH*((K1RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBC[t]*nadpRBC[t])/Vrbc^2 - (K10RBCvG6PDH*K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadphRBC[t]*p6glRBC[t])/Vrbc^2))/(K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH + (K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBC[t])/Vrbc + (K10RBCvG6PDH*K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*nadphRBC[t])/Vrbc + (K10RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*g6pRBC[t]*nadphRBC[t])/Vrbc^2 + (K1RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH*nadpRBC[t])/Vrbc + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K9RBCvG6PDH + K6RBCvG6PDH*K9RBCvG6PDH + K7RBCvG6PDH*K9RBCvG6PDH)*g6pRBC[t]*nadpRBC[t])/Vrbc^2 + (K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*p6glRBC[t])/Vrbc + (K10RBCvG6PDH*(K2RBCvG6PDH*K4RBCvG6PDH + K2RBCvG6PDH*K5RBCvG6PDH + K2RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K6RBCvG6PDH)*K8RBCvG6PDH*nadphRBC[t]*p6glRBC[t])/Vrbc^2 + (K10RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBC[t]*nadphRBC[t]*p6glRBC[t])/Vrbc^3 + (K1RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadpRBC[t]*p6glRBC[t])/Vrbc^2 + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBC[t]*nadpRBC[t]*p6glRBC[t])/Vrbc^3), vRBCvGAPDH -> (TsRBC*Vrbc*VRBCvGAPDH*(-((KcatrappRBCvGAPDH*b13pgRBC[t]*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbc^2)) + (KcatfappRBCvGAPDH*gapRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbc^3)))/((b13pgRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbc)))/(Kiappb13pgRBCvGAPDH*Vrbc) + (gapRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbc)))/(KiappgapRBCvGAPDH*Vrbc) + (Kmb13pgRBCvGAPDH*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbc) + (b13pgRBC[t]*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbc^2) + (gapRBC[t]*nadhRBC[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*Vrbc^2) + (b13pgRBC[t]*nadRBC[t])/(Kiappb13pgRBCvGAPDH*KinadRBCvGAPDH*Vrbc^2) + (gapRBC[t]*nadRBC[t])/(KiappgapRBCvGAPDH*KinadRBCvGAPDH*Vrbc^2) + (gapRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbc))*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*Vrbc^2) + (Kmb13pgRBCvGAPDH*nadhRBC[t]*phosRBC[t])/(Kiappb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbc^2) + (Kmb13pgRBCvGAPDH*b13pgRBC[t]*nadhRBC[t]*phosRBC[t])/(Kiappb13pgRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbc^3) + (gapRBC[t]*nadhRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*KiphosRBCvGAPDH*Vrbc^3) + (KmgapRBCvGAPDH*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbc^2) + (KmgapRBCvGAPDH*b13pgRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbc^3) + (gapRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbc^3)), vRBCvGLCTRANSPORT -> TsRBC*Vrbc*((glcAlpha*glcEXT*VRBCvGLCTRANSPORT)/((kAlpha + (glcAlpha*glcEXT)/vBld + ((1 - glcAlpha)*glcEXT*kAlpha)/(kBeta*vBld))*vBld) + ((1 - glcAlpha)*glcEXT*VRBCvGLCTRANSPORT)/((kBeta + ((1 - glcAlpha)*glcEXT)/vBld + (glcAlpha*glcEXT*kBeta)/(kAlpha*vBld))*vBld) - (glcAlpha*VRBCvGLCTRANSPORT*glcRBC[t])/(Vrbc*(kAlpha + (glcAlpha*glcRBC[t])/Vrbc + ((1 - glcAlpha)*kAlpha*glcRBC[t])/(kBeta*Vrbc))) - ((1 - glcAlpha)*VRBCvGLCTRANSPORT*glcRBC[t])/(Vrbc*(kBeta + ((1 - glcAlpha)*glcRBC[t])/Vrbc + (glcAlpha*kBeta*glcRBC[t])/(kAlpha*Vrbc)))), vRBCvGSSGR -> (TsRBC*Vrbc*VRBCvGSSGR*((K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t]*nadphRBC[t])/Vrbc^2 - (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2*nadpRBC[t])/Vrbc^3))/(K11RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR + (K11RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t])/Vrbc + (K10RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2)/Vrbc^2 + (K11RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t])/Vrbc + (K11RBCvGSSGR*K1RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadphRBC[t])/Vrbc + (K11RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]*nadphRBC[t])/Vrbc^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2*nadphRBC[t])/Vrbc^3 + (K1RBCvGSSGR*K3RBCvGSSGR*(K11RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR + K11RBCvGSSGR*K5RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K6RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR)*gssgRBC[t]*nadphRBC[t])/Vrbc^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadphRBC[t])/Vrbc^3 + (K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadphRBC[t])/Vrbc^3 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*gssgRBC[t]*nadphRBC[t])/Vrbc^4 + (K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadpRBC[t])/Vrbc + (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*gshRBC[t]*nadpRBC[t])/Vrbc^2 + (K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]*nadpRBC[t])/Vrbc^2 + (K10RBCvGSSGR*K12RBCvGSSGR*(K2RBCvGSSGR*K4RBCvGSSGR + K2RBCvGSSGR*K5RBCvGSSGR + K2RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*nadpRBC[t])/Vrbc^3 + (K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t]*nadpRBC[t])/Vrbc^2 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadpRBC[t])/Vrbc^3 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*gssgRBC[t]*nadpRBC[t])/Vrbc^4), vRBCvHBADP -> TsRBC*Vrbc*(-((KdRBCvHBADP*hbadpRBC[t])/Vrbc) + (KaappRBCvHBADP*adpRBC[t]*hbRBC[t])/Vrbc^2), vRBCvHBATP -> TsRBC*Vrbc*(-((KdRBCvHBATP*hbatpRBC[t])/Vrbc) + (KaappRBCvHBATP*atpRBC[t]*hbRBC[t])/Vrbc^2), vRBCvHBB13PG -> TsRBC*Vrbc*(-((KdRBCvHBB13PG*hbb13pgRBC[t])/Vrbc) + (KaappRBCvHBB13PG*b13pgRBC[t]*hbRBC[t])/Vrbc^2), vRBCvHBB23PG -> TsRBC*Vrbc*(-((KdRBCvHBBPG*hbb23pgRBC[t])/Vrbc) + (KaappRBCvHBBPG*b23pgRBC[t]*hbRBC[t])/Vrbc^2), vRBCvHBMGATP -> TsRBC*Vrbc*(-((KdRBCvHBMGATP*hbmgatpRBC[t])/Vrbc) + (KaappRBCvHBMGATP*hbRBC[t]*mgatpRBC[t])/Vrbc^2), vRBCvHK -> (TsRBC*Vrbc*(-((ERBCvHK*KcatrappRBCvHK*g6pRBC[t]*mgadpRBC[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbc^2)) + (ERBCvHK*KcatfappRBCvHK*glcRBC[t]*mgatpRBC[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbc^2)))/(1 + g6pRBC[t]/(Kig6pRBCvHK*Vrbc) + glcRBC[t]/(KiglcRBCvHK*Vrbc) + (b23pgRBC[t]*glcRBC[t])/(KdibpgRBCvHK*KiglcRBCvHK*Vrbc^2) + (g16p2RBC[t]*glcRBC[t])/(Kdig16p2RBCvHK*KiglcRBCvHK*Vrbc^2) + (g6pRBC[t]*glcRBC[t])/(Kdig6pRBCvHK*KiglcRBCvHK*Vrbc^2) + (glcRBC[t]*gshRBC[t])/(KdigshRBCvHK*KiglcRBCvHK*Vrbc^2) + mgadpRBC[t]/(KimgadpRBCvHK*Vrbc) + (g6pRBC[t]*mgadpRBC[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbc^2) + mgatpRBC[t]/(KimgatpRBCvHK*Vrbc) + (glcRBC[t]*mgatpRBC[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbc^2)), vRBCvLACTRANSPORT -> TsRBC*vBld*(-((KiRBCvLACTRANSPORT*lacEXT)/vBld) + (KoRBCvLACTRANSPORT*lacRBC[t])/Vrbc), vRBCvLDH -> (TsRBC*Vrbc*VRBCvLDH*(-((KcatrRBCvLDH*lacRBC[t]*nadRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbc^2)) + (KcatfRBCvLDH*nadhRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbc^2)))/(nadhRBC[t]/(KinadhRBCvLDH*Vrbc) + (KmnadRBCvLDH*lacRBC[t]*nadhRBC[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbc^2) + nadRBC[t]/(KinadRBCvLDH*Vrbc) + (lacRBC[t]*nadRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbc^2) + (nadhRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbc^2) + (lacRBC[t]*nadhRBC[t]*pyrRBC[t])/(KiapplacRBCvLDH*KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbc^3) + (KmnadhRBCvLDH*nadRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapppyrRBCvLDH*Vrbc^2) + (lacRBC[t]*nadRBC[t]*pyrRBC[t])/(KiapppyrRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbc^3) + (1 + pyrRBC[t]/(KidpyrRBCvLDH*Vrbc))*(1 + (KmnadRBCvLDH*lacRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbc) + (KmnadhRBCvLDH*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbc))), vRBCvLDHP -> (TsRBC*Vrbc*(-((KrRBCvLDHP*lacRBC[t]*nadpRBC[t])/(KmlacRBCvLDHP*Vrbc^2)) + (KfRBCvLDHP*nadphRBC[t]*pyrRBC[t])/(KmpyrRBCvLDHP*Vrbc^2)))/(1 + lacRBC[t]/(KmlacRBCvLDHP*Vrbc) + pyrRBC[t]/(KmpyrRBCvLDHP*Vrbc)), vRBCvMGADP -> TsRBC*Vrbc*(-((KdRBCvMGADP*mgadpRBC[t])/Vrbc) + (KaappRBCvMGADP*adpRBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGATP -> TsRBC*Vrbc*(-((KdRBCvMGATP*mgatpRBC[t])/Vrbc) + (KaappRBCvMGATP*atpRBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGB13PG -> TsRBC*Vrbc*(-((KdRBCvMGB13PG*mgb13pgRBC[t])/Vrbc) + (KaappRBCvMGB13PG*b13pgRBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGB23PG -> TsRBC*Vrbc*(-((KdRBCvMGB23PG*mgb23pgRBC[t])/Vrbc) + (KaappRBCvMGB23PG*b23pgRBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGF16P2 -> TsRBC*Vrbc*(-((KdRBCvMGF16P2*mgf16p2RBC[t])/Vrbc) + (KaappRBCvMGF16P2*f16p2RBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGG16P2 -> TsRBC*Vrbc*(-((KdRBCvMGG16P2*mgg16p2RBC[t])/Vrbc) + (KaappRBCvMGG16P2*g16p2RBC[t]*mgRBC[t])/Vrbc^2), vRBCvMGPHOS -> TsRBC*Vrbc*(-((KdRBCvMGPHOS*mgphosRBC[t])/Vrbc) + (KaappRBCvMGPHOS*mgRBC[t]*phosRBC[t])/Vrbc^2), vRBCvOX -> KRBCvOX*TsRBC*gshRBC[t], vRBCvOXNADH -> KRBCvOXNADH*TsRBC*nadhRBC[t], vRBCvP6GDH -> (TsRBC*Vrbc*VRBCvP6GDH*((K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadpRBC[t]*p6gRBC[t])/Vrbc^2 - (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBC[t]*ru5pRBC[t])/Vrbc^3))/(K11RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH + (co2RBC*K11RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH)/Vrbc + (co2RBC*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBC[t])/Vrbc^2 + (K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadphRBC[t])/Vrbc + (co2RBC*K11RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBC[t])/Vrbc^2 + (K11RBCvP6GDH*K1RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadpRBC[t])/Vrbc + (K11RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*p6gRBC[t])/Vrbc + (K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadphRBC[t]*p6gRBC[t])/Vrbc^2 + (co2RBC*K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBC[t]*p6gRBC[t])/Vrbc^3 + (K1RBCvP6GDH*K3RBCvP6GDH*(K11RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH + K11RBCvP6GDH*K5RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K6RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH)*nadpRBC[t]*p6gRBC[t])/Vrbc^2 + (co2RBC*K10RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*ru5pRBC[t])/Vrbc^2 + (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*(K2RBCvP6GDH*K4RBCvP6GDH + K2RBCvP6GDH*K5RBCvP6GDH + K2RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBC[t]*ru5pRBC[t])/Vrbc^3 + (K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*nadphRBC[t]*ru5pRBC[t])/Vrbc^2 + (co2RBC*K10RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBC[t]*ru5pRBC[t])/Vrbc^3 + (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbc^4 + (K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadphRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbc^3 + (co2RBC*K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbc^4 + (K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadpRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbc^3), vRBCvPFK -> (TsRBC*Vrbc*VRBCvPFK*(-((KcatfRBCvPFK*f16p2RBC[t]*mgadpRBC[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbc^2)) + (KcatfRBCvPFK*f6pRBC[t]*mgatpRBC[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbc^2)))/((1 + f16p2RBC[t]/(Kmf16p2RBCvPFK*Vrbc) + f6pRBC[t]/(Kmf6pRBCvPFK*Vrbc) + mgadpRBC[t]/(KmmgadpRBCvPFK*Vrbc) + (f16p2RBC[t]*mgadpRBC[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbc^2) + mgatpRBC[t]/(KmmgatpRBCvPFK*Vrbc) + (f6pRBC[t]*mgatpRBC[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbc^2))*(1 + ((hRBC/KaRBCvPFK)^nRBCvPFK*(1 + atpRBC[t]/(KtatpRBCvPFK*Vrbc))^4*(1 + b23pgRBC[t]/(Ktb23pgRBCvPFK*Vrbc))^4*(1 + mgRBC[t]/(KtmgRBCvPFK*Vrbc))^4)/((1 + ampRBC[t]/(KrampRBCvPFK*Vrbc))^4*(1 + f16p2RBC[t]/(Kmf16p2RBCvPFK*Vrbc) + f6pRBC[t]/(Kmf6pRBCvPFK*Vrbc))^4*(1 + g16p2RBC[t]/(Krg16p2RBCvPFK*Vrbc))^4*(1 + phosRBC[t]/(KrphosRBCvPFK*Vrbc))^4))), vRBCvPGI -> (ERBCvPGI*TsRBC*Vrbc*(-((KcatrRBCvPGI*f6pRBC[t])/(Kmf6pRBCvPGI*Vrbc)) + (KcatfRBCvPGI*g6pRBC[t])/(Kmg6pRBCvPGI*Vrbc)))/(1 + f6pRBC[t]/(Kmf6pRBCvPGI*Vrbc) + g6pRBC[t]/(Kmg6pRBCvPGI*Vrbc)), vRBCvPGK -> (TsRBC*Vrbc*VRBCvPGK*((KcatfRBCvPGK*b13pgRBC[t]*mgadpRBC[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbc^2) - (KcatrRBCvPGK*mgatpRBC[t]*p3gRBC[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbc^2)))/(1 + b13pgRBC[t]/(Kib13pgRBCvPGK*Vrbc) + mgadpRBC[t]/(KimgadpRBCvPGK*Vrbc) + (b13pgRBC[t]*mgadpRBC[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbc^2) + mgatpRBC[t]/(KimgatpRBCvPGK*Vrbc) + p3gRBC[t]/(Kip3gRBCvPGK*Vrbc) + (mgatpRBC[t]*p3gRBC[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbc^2)), vRBCvPGLHYDROLYSIS -> TsRBC*Vrbc*((KhydrolRBCvPGLHYDROLYSIS*p6glRBC[t])/Vrbc + (Klactonase3RBCvPGLHYDROLYSIS*VRBCvPGLHYDROLYSIS*p6glRBC[t])/(Vrbc*((Klactonase2RBCvPGLHYDROLYSIS + Klactonase3RBCvPGLHYDROLYSIS)/Klactonase1RBCvPGLHYDROLYSIS + p6glRBC[t]/Vrbc))), vRBCvPGM -> (TsRBC*Vrbc*VRBCvPGM*(-((KcatrRBCvPGM*p2gRBC[t])/(Kmp2gRBCvPGM*Vrbc)) + (KcatfRBCvPGM*p3gRBC[t])/(Kmp3gRBCvPGM*Vrbc)))/(1 + p2gRBC[t]/(Kmp2gRBCvPGM*Vrbc) + p3gRBC[t]/(Kmp3gRBCvPGM*Vrbc)), vRBCvPHOSTRANSPORT -> TsRBC*vBld*(-((KiRBCvPHOSTRANSPORT*phosEXT)/vBld) + (KoRBCvPHOSTRANSPORT*phosRBC[t])/Vrbc), vRBCvPK -> (TsRBC*Vrbc*VRBCvPK*((KcatfRBCvPK*mgadpRBC[t]*pepRBC[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbc^2) - (KcatrRBCvPK*mgatpRBC[t]*pyrRBC[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbc^2)))/((1 + LRBCvPK)*(1 + mgadpRBC[t]/(KrmgadpRBCvPK*Vrbc) + mgatpRBC[t]/(KrmgatpRBCvPK*Vrbc) + pepRBC[t]/(KrpepRBCvPK*Vrbc) + (mgadpRBC[t]*pepRBC[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbc^2) + pyrRBC[t]/(KrpyrRBCvPK*Vrbc) + (mgatpRBC[t]*pyrRBC[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbc^2))), vRBCvPYRTRANSPORT -> TsRBC*vBld*(-((KiRBCvPYRTRANSPORT*pyrEXT)/vBld) + (KoRBCvPYRTRANSPORT*pyrRBC[t])/Vrbc), vRBCvR5PI -> (TsRBC*Vrbc*VRBCvR5PI*(-((K2RBCvR5PI*K4RBCvR5PI*rib5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbc)) + (K1RBCvR5PI*K3RBCvR5PI*ru5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbc)))/(1 + (K1RBCvR5PI*ru5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbc) + (K4RBCvR5PI*xu5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbc)), vRBCvRu5PE -> (TsRBC*Vrbc*VRBCvRu5PE*((K1RBCvRu5PE*K3RBCvRu5PE*ru5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbc) - (K2RBCvRu5PE*K4RBCvRu5PE*xu5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbc)))/(1 + (K1RBCvRu5PE*ru5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbc) + (K4RBCvRu5PE*xu5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbc)), vRBCvTA -> (TsRBC*Vrbc*VRBCvTA*(-((K2RBCvTA*K4RBCvTA*K6RBCvTA*K8RBCvTA*ery4pRBC[t]*f6pRBC[t])/Vrbc^2) + (K1RBCvTA*K3RBCvTA*K5RBCvTA*K7RBCvTA*gapRBC[t]*sed7pRBC[t])/Vrbc^2))/((K2RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBC[t])/Vrbc + ((K2RBCvTA + K3RBCvTA)*K6RBCvTA*K8RBCvTA*f6pRBC[t])/Vrbc + (K4RBCvTA*(K2RBCvTA + K6RBCvTA)*K8RBCvTA*ery4pRBC[t]*f6pRBC[t])/Vrbc^2 + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K7RBCvTA*gapRBC[t])/Vrbc + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K8RBCvTA*f6pRBC[t]*gapRBC[t])/Vrbc^2 + (K1RBCvTA*K3RBCvTA*(K6RBCvTA + K7RBCvTA)*sed7pRBC[t])/Vrbc + (K1RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBC[t]*sed7pRBC[t])/Vrbc^2 + (K1RBCvTA*K5RBCvTA*(K3RBCvTA + K7RBCvTA)*gapRBC[t]*sed7pRBC[t])/Vrbc^2), vRBCvTIM -> (TsRBC*Vrbc*VRBCvTIM*(-((KcatrRBCvTIM*dhapRBC[t])/(KmdhapRBCvTIM*Vrbc)) + (KcatfRBCvTIM*gapRBC[t])/(KmgapRBCvTIM*Vrbc)))/(1 + dhapRBC[t]/(KmdhapRBCvTIM*Vrbc) + gapRBC[t]/(KmgapRBCvTIM*Vrbc)), vRBCvTK1 -> TsRBC*Vrbc*(-((K2RBCvTK1*tkxu5pRBC[t])/Vrbc) + (K1RBCvTK1*tkRBC[t]*xu5pRBC[t])/Vrbc^2), vRBCvTK2 -> TsRBC*Vrbc*(-((K4RBCvTK2*gapRBC[t]*tkgRBC[t])/Vrbc^2) + (K3RBCvTK2*tkxu5pRBC[t])/Vrbc), vRBCvTK3 -> TsRBC*Vrbc*((K5RBCvTK3*rib5pRBC[t]*tkgRBC[t])/Vrbc^2 - (K6RBCvTK3*tkgrib5pRBC[t])/Vrbc), vRBCvTK4 -> TsRBC*Vrbc*((K7RBCvTK4*tkgrib5pRBC[t])/Vrbc - (K8RBCvTK4*sed7pRBC[t]*tkRBC[t])/Vrbc^2), vRBCvTK5 -> TsRBC*Vrbc*(-((K10RBCvTK5*tkgery4pRBC[t])/Vrbc) + (K9RBCvTK5*ery4pRBC[t]*tkgRBC[t])/Vrbc^2), vRBCvTK6 -> TsRBC*Vrbc*((K11RBCvTK6*tkgery4pRBC[t])/Vrbc - (K12RBCvTK6*f6pRBC[t]*tkRBC[t])/Vrbc^2)
};

parameters = {
AlphaCellWaterFraction -> 0.7, CRBCvMGADP -> 0.00062, CRBCvMGATP -> 8.4*^-05, CRBCvMGB13PG -> 0.0032, CRBCvMGB23PG -> 0.0032, CRBCvMGF16P2 -> 0.0083, CRBCvMGG16P2 -> 0.0083, CRBCvMGPHOS -> 40800.0, ConcADPRBC -> 0.00031, ConcAMPRBC -> 3*^-05, ConcATPRBC -> 0.0021, ConcB13PGRBC -> 7*^-07, ConcB23PGRBC -> 0.0067, ConcBPGSPB13PGRBC -> 0.0, ConcBPGSPB23PGRBC -> 0.0, ConcBPGSPPP2GRBC -> 0.0, ConcBPGSPPP3GRBC -> 0.0, ConcBPGSPPRBC -> 0.0, ConcBPGSPRBC -> 3.8*^-06, ConcBpgsppphosRBC -> 0.0, ConcCo2RBC -> 0.0012, ConcDhapRBC -> 1.9*^-05, ConcEry4PRBC -> 1*^-05, ConcF16P2RBC -> 2.7*^-06, ConcF6PRBC -> 1.3*^-05, ConcG16P2RBC -> 0.000122, ConcG6PRBC -> 4*^-05, ConcGSHRBC -> 0.0032, ConcGapRBC -> 5.7*^-06, ConcGlcEXT -> 0.005, ConcGlcRBC -> 0.005, ConcGssgRBC -> 9*^-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.25*^-07, ConcNADRBC -> 6*^-05, ConcNadhRBC -> 1.4*^-07, ConcNadphRBC -> 6.4*^-05, ConcP2GARBC -> 1*^-05, ConcP3GARBC -> 6.4*^-05, ConcP6GLRBC -> 1.4*^-10, ConcP6GRBC -> 1.4*^-07, ConcPEPRBC -> 2.3*^-05, ConcPhosEXT -> 0.00192, ConcPhosRBC -> 0.001, ConcPyrEXT -> 8.5*^-05, ConcPyrRBC -> 6*^-05, ConcRib5PRBC -> 1*^-05, ConcRu5PRBC -> 1*^-05, ConcSed7PRBC -> 1*^-05, ConcTKGEry4PRBC -> 0.0, ConcTKGRBC -> 0.0, ConcTKGRib5PRBC -> 0.0, ConcTKRBC -> 3.3*^-07, ConcTkxu5pRBC -> 0.0, ConcXu5PRBC -> 1*^-06, ERBCvHK -> 2.5*^-08, ERBCvPGI -> 2.18*^-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.4*^-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.91251*^-08, KahbRBC -> 2.51189*^-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 -> 3*^-05, Kdig6pRBCvHK -> 1*^-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.55544*^-21, Kib13pgRBCvPGK -> 1.6*^-06, Kib23pgRBCvALD -> 0.0015, Kidb13pgRBCvGAPDH -> 1*^-06, KidgapRBCvGAPDH -> 3.1*^-05, KidhapRBCvALD -> 1.1*^-05, KidpyrRBCvLDH -> 0.000101, Kif16p2RBCvALD -> 1.98*^-05, Kig6pRBCvHK -> 4.7*^-05, KigapRBCvGAPDH -> 4.76523*^-19, KiglcRBCvHK -> 4.7*^-05, KilacRBCvLDH -> 0.00733, KimgRBCvENO -> 4.6*^-05, KimgadpRBCvHK -> 0.001, KimgadpRBCvPGK -> 8*^-05, KimgatpRBCvHK -> 0.001, KimgatpRBCvPGK -> 0.000186, KinadRBCvGAPDH -> 4.5*^-05, KinadRBCvLDH -> 0.000503, KinadhRBCvGAPDH -> 1*^-05, KinadhRBCvLDH -> 2.45*^-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.71*^-07, Kmb13pgRBCvPGK -> 2*^-06, KmdhapRBCvALD -> 3.5*^-05, KmdhapRBCvTIM -> 0.0001624, Kmf16p2RBCvALD -> 7.1*^-06, Kmf16p2RBCvPFK -> 0.0005, Kmf6pRBCvPFK -> 7.5*^-05, Kmf6pRBCvPGI -> 7.1*^-05, Kmg6pRBCvPGI -> 0.000181, KmgapRBCvALD -> 0.000189, KmgapRBCvGAPDH -> 9.5*^-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.6*^-05, KmmgadpRBCvHK -> 0.001, KmmgadpRBCvPFK -> 0.00054, KmmgadpRBCvPGK -> 0.0001, KmmgatpRBCvHK -> 0.001, KmmgatpRBCvPFK -> 6.8*^-05, KmmgatpRBCvPGK -> 0.001, KmnadRBCvGAPDH -> 4.5*^-05, KmnadRBCvLDH -> 0.000107, KmnadhRBCvGAPDH -> 3.3*^-06, KmnadhRBCvLDH -> 8.44*^-06, Kmp2gRBCvENO -> 0.00014, Kmp2gRBCvPGM -> 2.56*^-05, Kmp3gRBCvPGK -> 0.0011, Kmp3gRBCvPGM -> 0.000168, KmpepRBCvENO -> 0.0001105, KmphosRBCvGAPDH -> 0.00316, KmpyrRBCvLDH -> 0.000137, KmpyrRBCvLDHP -> 0.000414, KrRBCvLDHP -> 5.43*^-07, KrampRBCvPFK -> 0.0003, Krf16p2RBCvPK -> 5*^-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.7*^-07, VRBCvENO -> 2.2*^-07, VRBCvG6PDH -> 9.3*^-08, VRBCvGAPDH -> 7.66*^-06, VRBCvGLCTRANSPORT -> 0.047, VRBCvGSSGR -> 1.25*^-07, VRBCvLDH -> 3.43*^-06, VRBCvP6GDH -> 2.1*^-06, VRBCvPFK -> 1.1*^-07, VRBCvPGK -> 2.74*^-06, VRBCvPGLHYDROLYSIS -> 1.4*^-05, VRBCvPGM -> 4.1*^-07, VRBCvPK -> 8.7*^-08, VRBCvR5PI -> 1.42*^-05, VRBCvRu5PE -> 4.22*^-06, VRBCvTA -> 6.9*^-07, VRBCvTIM -> 1.14*^-06, bpgxRBC -> 3.8*^-06, glcAlpha -> 0.36, gssgxRBC -> 0.00160009, hbxRBC -> 0.007, kAlpha -> 0.0083, kBeta -> 0.007, kRBC -> 0.15, nRBCvPFK -> 5.0, nadpxRBC -> 6.4125*^-05, pHConversionFactor -> 1.0, phRBC -> 7.2, tkxRBC -> 3.3*^-07, xatpRBC -> -0.00456, xg16p2RBC -> 0.000122, xphosRBC -> 0.000438, Vrbc -> 90.0, vBld -> 128.5
};

assignments = {
pyrEXT -> ConcPyrEXT*vBld, KcatrappRBCvHK -> KcatrRBCvHK/(1 + 1.0471285480508985*^7/10^phRBC + 2.8183829312644494*^-10*10^phRBC), KmapppyrRBCvLDH -> 0.2847472489508012*(1 + 1.584893192461114*^-7*10^phRBC)*KmpyrRBCvLDH, KmappnadhRBCvGAPDH -> 6.30957344480193*^-8*10^phRBC*KmnadhRBCvGAPDH, KeqRBCvLACTRANSPORT -> (1 + 10^(-3.73 + phRBC))/(1 + 10^(-3.73 + phRBC)/RtvRBC), KcatrappRBCvGAPDH -> KcatrRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.1622776601683795*^7/10^phRBC), co2RBC -> ConcCo2RBC*Vrbc, KiapppyrRBCvLDH -> 0.2847472489508012*(1 + 1.584893192461114*^-7*10^phRBC)*KipyrRBCvLDH, KoRBCvPHOSTRANSPORT -> KiRBCvPHOSTRANSPORT/KeqRBCvPHOSTRANSPORT, KaappRBCvHBMGATP -> HbpHRBC*KaRBCvHBMGATP, K2appRBCvAK -> K2RBCvAK*(1 + hRBC*KhampRBCvAK + KkampRBCvAK*kRBC), KiappnadhRBCvGAPDH -> 6.30957344480193*^-8*10^phRBC*KinadhRBCvGAPDH, KaappRBCvHBBPG -> HbpHRBC*KaRBCvHBBPG, K6appRBCvBPGSP4 -> (1.758577575029182*K6RBCvBPGSP4)/(1 + 4.786300923226383*^28/10^(4*phRBC)), KaappRBCvMGG16P2 -> (CRBCvMGG16P2*KaRBCvMGG16P2*(hRBC*KhfRBCvMGG16P2*KmghfRBCvMGG16P2 + KRBCvMGG16P2))/(1 + hRBC*KhfRBCvMGG16P2 + hRBC^2*Kh2fRBCvMGG16P2*KhfRBCvMGG16P2 + KkfRBCvMGG16P2*kRBC + hRBC*KhfRBCvMGG16P2*KkhfRBCvMGG16P2*kRBC), KaappRBCvMGATP -> (CRBCvMGATP*KaRBCvMGATP*(hRBC*KhatpRBCvMGATP*KmghatpRBCvMGATP + KRBCvMGATP))/(1 + KkatpRBCvMGATP*kRBC + KhatpRBCvMGATP/(10^phRBC*pHConversionFactor)), K13appRBCvBPGSP7 -> (1.758577575029182*K13RBCvBPGSP7)/(1 + 4.786300923226383*^28/10^(4*phRBC)), KiappgapRBCvGAPDH -> KigapRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.1622776601683795*^7/10^phRBC), lacEXT -> ConcLacEXT*vBld, phosEXT -> ConcPhosEXT*vBld, KoRBCvLACTRANSPORT -> KiRBCvLACTRANSPORT/KeqRBCvLACTRANSPORT, KaappRBCvHBADP -> HbpHRBC*KaRBCvHBADP, KaappRBCvHBATP -> HbpHRBC*KaRBCvHBATP, hRBC -> pHConversionFactor/10^phRBC, KoRBCvPYRTRANSPORT -> KiRBCvPYRTRANSPORT/RtvRBC, K4appRBCvBPGSP3 -> (1.758577575029182*K4RBCvBPGSP3)/(1 + 4.786300923226383*^28/10^(4*phRBC)), KiapplacRBCvLDH -> 0.7152527510491988*(1 + 6.30957344480193*^6/10^phRBC)*KilacRBCvLDH, KaappRBCvHBB13PG -> HbpHRBC*KaRBCvHBB13PG, KeqRBCvPHOSTRANSPORT -> (1 + 10^(-6.75 + phRBC))/(10^(-6.75 + phRBC)/RtvRBC^2 + RtvRBC^(-1)), testsum -> adpRBC[t] + ampRBC[t] + atpRBC[t] + mgadpRBC[t] + mgatpRBC[t] + mgb13pgRBC[t] + mgb23pgRBC[t] + mgf16p2RBC[t] + mgg16p2RBC[t] + mgphosRBC[t] + mgRBC[t], Kiappb13pgRBCvGAPDH -> Kib13pgRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.1622776601683795*^7/10^phRBC), LRBCvPK -> (1.584893192461114*^-7*pHConversionFactor*(1 + atpRBC[t]/(KtatpRBCvPK*Vrbc))^4)/(hRBC*(1 + f16p2RBC[t]/(Krf16p2RBCvPK*Vrbc) + g16p2RBC[t]/(Krg16p2RBCvPK*Vrbc))^4*(1 + pepRBC[t]/(KrpepRBCvPK*Vrbc) + pyrRBC[t]/(KrpyrRBCvPK*Vrbc))^4), KaappRBCvMGPHOS -> (KaRBCvMGPHOS*(1 + KkphosRBCvMGPHOS*kRBC + 6.30957344480193*^-8*KhphosRBCvMGPHOS*pHConversionFactor))/(1 + hRBC*KhphosRBCvMGPHOS + KkphosRBCvMGPHOS*kRBC), KaappRBCvMGADP -> (CRBCvMGADP*KaRBCvMGADP*(hRBC*KhadpRBCvMGADP*KmghadpRBCvMGADP + KRBCvMGADP))/(1 + hRBC*KhadpRBCvMGADP + KkadpRBCvMGADP*kRBC), KmapplacRBCvLDH -> 0.7152527510491988*(1 + 6.30957344480193*^6/10^phRBC)*KmlacRBCvLDH, KcatfappRBCvHK -> KcatfRBCvHK/(1 + 1.0471285480508985*^7/10^phRBC + 2.8183829312644494*^-10*10^phRBC), KcatfappRBCvGAPDH -> KcatfRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.1622776601683795*^7/10^phRBC), atpRBCsum -> adpRBC[t] + ampRBC[t] + atpRBC[t], K3appRBCvBPGSP2 -> (1.758577575029182*K3RBCvBPGSP2)/(1 + 4.786300923226383*^28/10^(4*phRBC)), K1appRBCvBPGSP1 -> (3.5118864315095824*K1RBCvBPGSP1)/(1 + 1.584893192461114*^-7*10^phRBC), KaappRBCvMGF16P2 -> (CRBCvMGF16P2*KaRBCvMGF16P2*(hRBC*KhfRBCvMGF16P2*KmghfRBCvMGF16P2 + KRBCvMGF16P2))/(1 + hRBC*KhfRBCvMGF16P2 + hRBC^2*Kh2fRBCvMGF16P2*KhfRBCvMGF16P2 + KkfRBCvMGF16P2*kRBC + hRBC*KhfRBCvMGF16P2*KkhfRBCvMGF16P2*kRBC), K1appRBCvAK -> K1RBCvAK*(1 + KkadpRBCvAK*kRBC + (KhadpRBCvAK*pHConversionFactor)/10^phRBC), KaappRBCvMGB13PG -> (CRBCvMGB13PG*KaRBCvMGB13PG*(hRBC*KhbpgRBCvMGB13PG*KmghbpgRBCvMGB13PG + KRBCvMGB13PG))/(1 + hRBC*KhbpgRBCvMGB13PG + hRBC^2*Kh2bpgRBCvMGB13PG*KhbpgRBCvMGB13PG + KkbpgRBCvMGB13PG*kRBC + hRBC*KhbpgRBCvMGB13PG*KkhbpgRBCvMGB13PG*kRBC), HbpHRBC -> (1 + (2.5118864315095825*^14*KahbRBC^2)/pHConversionFactor^2 + (3.1697863849222284*^7*KahbRBC)/pHConversionFactor)/(1 + (2*KahbRBC)/hRBC + KahbRBC^2/hRBC^2), glcEXT -> ConcGlcEXT*vBld, KaappRBCvMGB23PG -> (CRBCvMGB23PG*KaRBCvMGB23PG*(hRBC*KhbpgRBCvMGB23PG*KmghbpgRBCvMGB23PG + KRBCvMGB23PG))/(1 + hRBC*KhbpgRBCvMGB23PG + hRBC^2*Kh2bpgRBCvMGB23PG*KhbpgRBCvMGB23PG + KkbpgRBCvMGB23PG*kRBC + hRBC*KhbpgRBCvMGB23PG*KkhbpgRBCvMGB23PG*kRBC)
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "min", "metabolite" -> "fmol", "extent" -> "fmol"}
};

(* Time evolution *)
odes = {
adpRBC'[t] ==  -1.0*vRBCvHBADP -1.0*vRBCvMGADP -1.0*vRBCvAK, ampRBC'[t] == 1.0*vRBCvAK , atpRBC'[t] ==  -1.0*vRBCvMGATP -1.0*vRBCvHBATP, b13pgRBC'[t] == 1.0*vRBCvGAPDH -1.0*vRBCvMGB13PG -1.0*vRBCvPGK -1.0*vRBCvBPGSP1 -1.0*vRBCvHBB13PG, b23pgRBC'[t] == 1.0*vRBCvBPGSP7 -1.0*vRBCvMGB23PG -1.0*vRBCvHBB23PG, bpgspRBC'[t] == 1.0*vRBCvBPGSP7 +1.0*vRBCvBPGSP9 -1.0*vRBCvBPGSP1, bpgspb13pgRBC'[t] == 1.0*vRBCvBPGSP1 -1.0*vRBCvBPGSP2, bpgspb23pgRBC'[t] == 1.0*vRBCvBPGSP5 +1.0*vRBCvBPGSP6 -1.0*vRBCvBPGSP7, bpgsppRBC'[t] == 1.0*vRBCvBPGSP2 -1.0*vRBCvBPGSP8 -1.0*vRBCvBPGSP3 -1.0*vRBCvBPGSP4, bpgsppp2gRBC'[t] == 1.0*vRBCvBPGSP4 -1.0*vRBCvBPGSP6, bpgsppp3gRBC'[t] == 1.0*vRBCvBPGSP3 -1.0*vRBCvBPGSP5, bpgsppphosRBC'[t] == 1.0*vRBCvBPGSP8 -1.0*vRBCvBPGSP9, dhapRBC'[t] == 1.0*vRBCvALD +1.0*vRBCvTIM , ery4pRBC'[t] == 1.0*vRBCvTA -1.0*vRBCvTK5, f16p2RBC'[t] == 1.0*vRBCvPFK -1.0*vRBCvMGF16P2 -1.0*vRBCvALD, f6pRBC'[t] == 1.0*vRBCvPGI +1.0*vRBCvTA +1.0*vRBCvTK6 -1.0*vRBCvPFK, g16p2RBC'[t] ==  -1.0*vRBCvMGG16P2, g6pRBC'[t] == 1.0*vRBCvHK -1.0*vRBCvG6PDH -1.0*vRBCvPGI, gapRBC'[t] == 1.0*vRBCvTK2 +1.0*vRBCvALD -1.0*vRBCvGAPDH -1.0*vRBCvTA -1.0*vRBCvTIM, glcRBC'[t] == 1.0*vRBCvGLCTRANSPORT -1.0*vRBCvHK, gshRBC'[t] == 2.0*vRBCvGSSGR -2.0*vRBCvOX, gssgRBC'[t] == 1.0*vRBCvOX -1.0*vRBCvGSSGR, hbRBC'[t] ==  -1.0*vRBCvHBB23PG -1.0*vRBCvHBMGATP -1.0*vRBCvHBADP -1.0*vRBCvHBATP -1.0*vRBCvHBB13PG, hbadpRBC'[t] == 1.0*vRBCvHBADP , hbatpRBC'[t] == 1.0*vRBCvHBATP , hbb13pgRBC'[t] == 1.0*vRBCvHBB13PG , hbb23pgRBC'[t] == 1.0*vRBCvHBB23PG , hbmgatpRBC'[t] == 1.0*vRBCvHBMGATP , lacRBC'[t] == 1.0*vRBCvLDH +1.0*vRBCvLDHP -1.0*vRBCvLACTRANSPORT, mgRBC'[t] ==  -1.0*vRBCvMGADP -1.0*vRBCvMGATP -1.0*vRBCvMGB13PG -1.0*vRBCvMGB23PG -1.0*vRBCvMGF16P2 -1.0*vRBCvMGG16P2 -1.0*vRBCvMGPHOS, mgadpRBC'[t] == 1.0*vRBCvHK +1.0*vRBCvATPASE +1.0*vRBCvMGADP +1.0*vRBCvPFK -1.0*vRBCvPGK -1.0*vRBCvPK -1.0*vRBCvAK, mgatpRBC'[t] == 1.0*vRBCvMGATP +1.0*vRBCvPGK +1.0*vRBCvPK +1.0*vRBCvAK -1.0*vRBCvHBMGATP -1.0*vRBCvHK -1.0*vRBCvPFK -1.0*vRBCvATPASE, mgb13pgRBC'[t] == 1.0*vRBCvMGB13PG , mgb23pgRBC'[t] == 1.0*vRBCvMGB23PG , mgf16p2RBC'[t] == 1.0*vRBCvMGF16P2 , mgg16p2RBC'[t] == 1.0*vRBCvMGG16P2 , mgphosRBC'[t] == 1.0*vRBCvMGPHOS , nadRBC'[t] == 1.0*vRBCvLDH +1.0*vRBCvOXNADH -1.0*vRBCvGAPDH, nadhRBC'[t] == 1.0*vRBCvGAPDH -1.0*vRBCvLDH -1.0*vRBCvOXNADH, nadpRBC'[t] == 1.0*vRBCvLDHP +1.0*vRBCvGSSGR -1.0*vRBCvG6PDH -1.0*vRBCvP6GDH, nadphRBC'[t] == 1.0*vRBCvP6GDH +1.0*vRBCvG6PDH -1.0*vRBCvGSSGR -1.0*vRBCvLDHP, p2gRBC'[t] == 1.0*vRBCvPGM -1.0*vRBCvENO -1.0*vRBCvBPGSP4, p3gRBC'[t] == 1.0*vRBCvPGK +1.0*vRBCvBPGSP2 -1.0*vRBCvPGM -1.0*vRBCvBPGSP3, p6gRBC'[t] == 1.0*vRBCvPGLHYDROLYSIS -1.0*vRBCvP6GDH, p6glRBC'[t] == 1.0*vRBCvG6PDH -1.0*vRBCvPGLHYDROLYSIS, pepRBC'[t] == 1.0*vRBCvENO -1.0*vRBCvPK, phosRBC'[t] == 1.0*vRBCvATPASE +2.0*vRBCvBPGSP9 -1.0*vRBCvBPGSP8 -1.0*vRBCvGAPDH -1.0*vRBCvMGPHOS -1.0*vRBCvPHOSTRANSPORT, pyrRBC'[t] == 1.0*vRBCvPK -1.0*vRBCvLDH -1.0*vRBCvLDHP -1.0*vRBCvPYRTRANSPORT, rib5pRBC'[t] == 1.0*vRBCvR5PI -1.0*vRBCvTK3, ru5pRBC'[t] == 1.0*vRBCvP6GDH -1.0*vRBCvR5PI -1.0*vRBCvRu5PE, sed7pRBC'[t] == 1.0*vRBCvTK4 -1.0*vRBCvTA, tkRBC'[t] == 1.0*vRBCvTK4 +1.0*vRBCvTK6 -1.0*vRBCvTK1, tkgRBC'[t] == 1.0*vRBCvTK2 -1.0*vRBCvTK3 -1.0*vRBCvTK5, tkgery4pRBC'[t] == 1.0*vRBCvTK5 -1.0*vRBCvTK6, tkgrib5pRBC'[t] == 1.0*vRBCvTK3 -1.0*vRBCvTK4, tkxu5pRBC'[t] == 1.0*vRBCvTK1 -1.0*vRBCvTK2, xu5pRBC'[t] == 1.0*vRBCvRu5PE -1.0*vRBCvTK1
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]