(* 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 = {
ADPc[t], ADPg[t], AMPc[t], AMPg[t], ATPc[t], ATPg[t], BPGA13c[t], BPGA13g[t], DHAPc[t], DHAPg[t], Fru16BPc[t], Fru16BPg[t], Fru6Pc[t], Fru6Pg[t], GA3Pc[t], GA3Pg[t], Glc6Pc[t], Glc6Pg[t], Glcc[t], Glcg[t], Gly3Pc[t], Gly3Pg[t], Glyc[t], Glyg[t], NADHg[t], NADg[t], PEPc[t], PGA2c[t], PGA3c[t], PGA3g[t], Pyrc[t]
};

initialValues = {
ADPc[0] == 0.77, ADPg[0] == 0.77, AMPc[0] == 0.83, AMPg[0] == 0.83, ATPc[0] == 2.25, ATPg[0] == 2.25, BPGA13c[0] == 0.01, BPGA13g[0] == 0.4, DHAPc[0] == 4.6, DHAPg[0] == 4.6, Fru16BPc[0] == 0.01, Fru16BPg[0] == 13.0, Fru6Pc[0] == 0.2, Fru6Pg[0] == 0.2, GA3Pc[0] == 0.01, GA3Pg[0] == 0.14, Glc6Pc[0] == 0.8, Glc6Pg[0] == 0.8, Glcc[0] == 0.01, Glcg[0] == 0.01, Gly3Pc[0] == 1.5, Gly3Pg[0] == 1.5, Glyc[0] == 1*^-06, Glyg[0] == 1*^-06, NADHg[0] == 0.1, NADg[0] == 3.9, PEPc[0] == 1.4, PGA2c[0] == 0.5, PGA3c[0] == 17.0, PGA3g[0] == 17.0, Pyrc[0] == 14.0
};

rates = {
ADPtg, AKc, AKg, ALDc, ALDg, AMPtg, ATPtg, ATPuc, BPGA13tg, DHAPtg, ENOc, F16BPtg, F6Ptg, G3PDHc, G3PDHg, G6PPc, G6Ptg, GA3Ptg, GAPDHc, GAPDHg, GKc, GKg, GPOc, GlcTc, GlcTg, Gly3Ptg, GlyTc, Glytg, HXKc, HXKg, PFKc, PFKg, PGA3Tg, PGAMc, PGIc, PGIg, PGKc, PGKg, PYKc, PyrTc, TPIc, TPIg
};

rateEquations = {
ADPtg -> (ADPtgADPtgk*(-ADPc[t] + ADPg[t]))/cytosol, AKc -> (AKcAKck1*ADPc[t]^2 - AKcAKck2*AMPc[t]*ATPc[t])/cytosol, AKg -> (AKgAKgk1*ADPg[t]^2 - AKgAKgk2*AMPg[t]*ATPg[t])/cytosol, ALDc -> (ALDcALDcVmax*Fru16BPc[t]*(1 - (DHAPc[t]*GA3Pc[t])/(ALDcALDgKeq*Fru16BPc[t])))/(ALDcALDgKmFru16BP*cytosol*(1 + ADPc[t]/ALDcALDgKiADP + AMPc[t]/ALDcALDgKiAMP + ATPc[t]/ALDcALDgKiATP)*(1 + DHAPc[t]/ALDcALDgKmDHAP + Fru16BPc[t]/(ALDcALDgKmFru16BP*(1 + ADPc[t]/ALDcALDgKiADP + AMPc[t]/ALDcALDgKiAMP + ATPc[t]/ALDcALDgKiATP)) + GA3Pc[t]/ALDcALDgKmGA3P + (DHAPc[t]*GA3Pc[t])/(ALDcALDgKmDHAP*ALDcALDgKmGA3P) + (Fru16BPc[t]*GA3Pc[t])/(ALDcALDgKiGA3P*ALDcALDgKmFru16BP*(1 + ADPc[t]/ALDcALDgKiADP + AMPc[t]/ALDcALDgKiAMP + ATPc[t]/ALDcALDgKiATP)))), ALDg -> (ALDgALDgVmax*Fru16BPg[t]*(1 - (DHAPg[t]*GA3Pg[t])/(ALDgALDgKeq*Fru16BPg[t])))/(ALDgALDgKmFru16BP*cytosol*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)*(1 + DHAPg[t]/ALDgALDgKmDHAP + Fru16BPg[t]/(ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)) + GA3Pg[t]/ALDgALDgKmGA3P + (DHAPg[t]*GA3Pg[t])/(ALDgALDgKmDHAP*ALDgALDgKmGA3P) + (Fru16BPg[t]*GA3Pg[t])/(ALDgALDgKiGA3P*ALDgALDgKmFru16BP*(1 + ADPg[t]/ALDgALDgKiADP + AMPg[t]/ALDgALDgKiAMP + ATPg[t]/ALDgALDgKiATP)))), AMPtg -> (AMPtgAMPtgk*(-AMPc[t] + AMPg[t]))/cytosol, ATPtg -> (ATPtgATPtgk*(-ATPc[t] + ATPg[t]))/cytosol, ATPuc -> (ATPucATPuck*ATPc[t])/(cytosol*ADPc[t]), BPGA13tg -> (BPGA13tgBPGA13tgk*(-BPGA13c[t] + BPGA13g[t]))/cytosol, DHAPtg -> (DHAPtgDHAPtgk*(DHAPc[t] - DHAPg[t]))/cytosol, ENOc -> (ENOcENOcVmax*(1 - PEPc[t]/(ENOcENOcKeq*PGA2c[t]))*PGA2c[t])/(cytosol*ENOcENOcKmPGA2*(1 + PEPc[t]/ENOcENOcKmPEP + PGA2c[t]/ENOcENOcKmPGA2)), F16BPtg -> (F16BPtgF16BPtgk*(-Fru16BPc[t] + Fru16BPg[t]))/cytosol, F6Ptg -> (F6PtgF6Ptgk*(-Fru6Pc[t] + Fru6Pg[t]))/cytosol, G3PDHc -> (G3PDHcG3PDHcVmax*NADHc*DHAPc[t]*(1 - (NADc*Gly3Pc[t])/(G3PDHcG3PDHgKeq*NADHc*DHAPc[t])))/(cytosol*G3PDHcG3PDHgKmDHAP*G3PDHcG3PDHgKmNADH*(1 + NADc/G3PDHcG3PDHgKmNAD + NADHc/G3PDHcG3PDHgKmNADH)*(1 + DHAPc[t]/G3PDHcG3PDHgKmDHAP + Gly3Pc[t]/G3PDHcG3PDHgKmGly3P)), G3PDHg -> (G3PDHgG3PDHgVmax*DHAPg[t]*(1 - (Gly3Pg[t]*NADg[t])/(G3PDHgG3PDHgKeq*DHAPg[t]*NADHg[t]))*NADHg[t])/(cytosol*G3PDHgG3PDHgKmDHAP*G3PDHgG3PDHgKmNADH*(1 + DHAPg[t]/G3PDHgG3PDHgKmDHAP + Gly3Pg[t]/G3PDHgG3PDHgKmGly3P)*(1 + NADg[t]/G3PDHgG3PDHgKmNAD + NADHg[t]/G3PDHgG3PDHgKmNADH)), G6PPc -> (G6PPcG6PPcVmax*Glc6Pc[t]*(1 - Glcc[t]/(G6PPcG6PPcKeq*Glc6Pc[t])))/(cytosol*G6PPcG6PPcKmGlc6P*(1 + Glc6Pc[t]/G6PPcG6PPcKmGlc6P + Glcc[t]/G6PPcG6PPcKmGlc)), G6Ptg -> (G6PtgG6Ptgk*(-Glc6Pc[t] + Glc6Pg[t]))/cytosol, GA3Ptg -> (GA3PtgGA3Ptgk*(-GA3Pc[t] + GA3Pg[t]))/cytosol, GAPDHc -> (GAPDHcGAPDHcVmax*NADc*(1 - (NADHc*BPGA13c[t])/(GAPDHcGAPDHcKeq*NADc*GA3Pc[t]))*GA3Pc[t])/(cytosol*GAPDHcGAPDHcKmGA3P*GAPDHcGAPDHcKmNAD*(1 + NADc/GAPDHcGAPDHcKmNAD + NADHc/GAPDHcGAPDHcKmNADH)*(1 + BPGA13c[t]/GAPDHcGAPDHcKmBPGA13 + GA3Pc[t]/GAPDHcGAPDHcKmGA3P)), GAPDHg -> (GAPDHgGAPDHgVmax*GA3Pg[t]*NADg[t]*(1 - (BPGA13g[t]*NADHg[t])/(GAPDHgGAPDHgKeq*GA3Pg[t]*NADg[t])))/(cytosol*GAPDHgGAPDHgKmGA3P*GAPDHgGAPDHgKmNAD*(1 + BPGA13g[t]/GAPDHgGAPDHgKmBPGA13 + GA3Pg[t]/GAPDHgGAPDHgKmGA3P)*(1 + NADg[t]/GAPDHgGAPDHgKmNAD + NADHg[t]/GAPDHgGAPDHgKmNADH)), GKc -> (GKcGKcVmax*ADPc[t]*Gly3Pc[t]*(1 - (ATPc[t]*Glyc[t])/(GKcGKgKeq*ADPc[t]*Gly3Pc[t])))/(cytosol*GKcGKgKmADP*GKcGKgKmGly3P*(1 + ADPc[t]/GKcGKgKmADP + ATPc[t]/GKcGKgKmATP)*(1 + Gly3Pc[t]/GKcGKgKmGly3P + Glyc[t]/GKcGKgKmGly)), GKg -> (GKgGKgVmax*ADPg[t]*Gly3Pg[t]*(1 - (ATPg[t]*Glyg[t])/(GKgGKgKeq*ADPg[t]*Gly3Pg[t])))/(cytosol*GKgGKgKmADP*GKgGKgKmGly3P*(1 + ADPg[t]/GKgGKgKmADP + ATPg[t]/GKgGKgKmATP)*(1 + Gly3Pg[t]/GKgGKgKmGly3P + Glyg[t]/GKgGKgKmGly)), GPOc -> (GPOcGPOcVmax*Gly3Pc[t])/(cytosol*GPOcGPOcKmGly3P*(1 + ((1 + SHAM/GPOcGPOcKiSHAM)*Gly3Pc[t])/GPOcGPOcKmGly3P)), GlcTc -> (GlcTcGlcTcVmax*(Glce - Glcc[t]))/(cytosol*(Glce + GlcTcGlcTcKmGlc + Glcc[t] + (Glce*GlcTcGlcTcalpha*Glcc[t])/GlcTcGlcTcKmGlc)), GlcTg -> (GlcTgGlcTgk*Glcc[t] - GlcTgGlcTgk*Glcg[t])/cytosol, Gly3Ptg -> (Gly3PtgGly3Ptgk*(-Gly3Pc[t] + Gly3Pg[t]))/cytosol, GlyTc -> (GlyTcGlyTck*(-Glye + Glyc[t]) + (GlyTcGlyTcVrmax*GlyTcGlyTcy*(-Glye + Glyc[t]))/(Glye*GlyTcGlyTcy + GlyTcGlyTcKmGlye*GlyTcGlyTcy + Glyc[t] + (Glye*GlyTcGlyTcx*(1 - GlyTcGlyTcx + GlyTcGlyTcy)*Glyc[t])/(GlyTcGlyTcKmGlye*GlyTcGlyTcy)))/cytosol, Glytg -> (GlytgGlytgk*(-Glyc[t] + Glyg[t]))/cytosol, HXKc -> (HXKcHXKcVmax*ATPc[t]*(1 - (ADPc[t]*Glc6Pc[t])/(HXKcHXKgKeq*ATPc[t]*Glcc[t]))*Glcc[t])/(cytosol*HXKcHXKgKmATP*HXKcHXKgKmGlc*(1 + ADPc[t]/HXKcHXKgKmADP + ATPc[t]/HXKcHXKgKmATP)*(1 + Glc6Pc[t]/HXKcHXKgKmGlc6P + Glcc[t]/HXKcHXKgKmGlc)), HXKg -> (HXKgHXKgVmax*ATPg[t]*(1 - (ADPg[t]*Glc6Pg[t])/(HXKgHXKgKeq*ATPg[t]*Glcg[t]))*Glcg[t])/(cytosol*HXKgHXKgKmATP*HXKgHXKgKmGlc*(1 + ADPg[t]/HXKgHXKgKmADP + ATPg[t]/HXKgHXKgKmATP)*(1 + Glc6Pg[t]/HXKgHXKgKmGlc6P + Glcg[t]/HXKgHXKgKmGlc)), PFKc -> (PFKcPFKcVmax*PFKcPFKgKi1*ATPc[t]*(1 - (ADPc[t]*Fru16BPc[t])/(PFKcPFKgKeq*ATPc[t]*Fru6Pc[t]))*Fru6Pc[t])/(cytosol*PFKcPFKgKmATP*PFKcPFKgKmFru6P*(PFKcPFKgKi1 + Fru16BPc[t])*(PFKcPFKgKsATP/PFKcPFKgKmATP + ADPc[t]/PFKcPFKgKmADP + ATPc[t]/PFKcPFKgKmATP + (ADPc[t]*Fru16BPc[t])/(PFKcPFKgKi2*PFKcPFKgKmADP) + Fru6Pc[t]/PFKcPFKgKmFru6P + (ATPc[t]*Fru6Pc[t])/(PFKcPFKgKmATP*PFKcPFKgKmFru6P))), PFKg -> (PFKgPFKgKi1*PFKgPFKgVmax*ATPg[t]*(1 - (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKeq*ATPg[t]*Fru6Pg[t]))*Fru6Pg[t])/(cytosol*PFKgPFKgKmATP*PFKgPFKgKmFru6P*(PFKgPFKgKi1 + Fru16BPg[t])*(PFKgPFKgKsATP/PFKgPFKgKmATP + ADPg[t]/PFKgPFKgKmADP + ATPg[t]/PFKgPFKgKmATP + (ADPg[t]*Fru16BPg[t])/(PFKgPFKgKi2*PFKgPFKgKmADP) + Fru6Pg[t]/PFKgPFKgKmFru6P + (ATPg[t]*Fru6Pg[t])/(PFKgPFKgKmATP*PFKgPFKgKmFru6P))), PGA3Tg -> (-(PGA3TgPGA3Tgk*PGA3c[t]) + PGA3TgPGA3Tgk*PGA3g[t])/cytosol, PGAMc -> (PGAMcPGAMcVmax*(1 - PGA2c[t]/(PGAMcPGAMcKeq*PGA3c[t]))*PGA3c[t])/(cytosol*PGAMcPGAMcKmPGA3*(1 + PGA2c[t]/PGAMcPGAMcKmPGA2 + PGA3c[t]/PGAMcPGAMcKmPGA3)), PGIc -> (PGIcPGIcVmax*(1 - Fru6Pc[t]/(PGIcPGIgKeq*Glc6Pc[t]))*Glc6Pc[t])/(cytosol*PGIcPGIgKmGlc6P*(1 + Fru6Pc[t]/PGIcPGIgKmFru6P + Glc6Pc[t]/PGIcPGIgKmGlc6P)), PGIg -> (PGIgPGIgVmax*(1 - Fru6Pg[t]/(PGIgPGIgKeq*Glc6Pg[t]))*Glc6Pg[t])/(cytosol*PGIgPGIgKmGlc6P*(1 + Fru6Pg[t]/PGIgPGIgKmFru6P + Glc6Pg[t]/PGIgPGIgKmGlc6P)), PGKc -> (PGKcPGKcVmax*ADPc[t]*BPGA13c[t]*(1 - (ATPc[t]*PGA3c[t])/(PGKcPGKgKeq*ADPc[t]*BPGA13c[t])))/(cytosol*PGKcPGKgKmADP*PGKcPGKgKmBPGA13*(1 + ADPc[t]/PGKcPGKgKmADP + ATPc[t]/PGKcPGKgKmATP)*(1 + BPGA13c[t]/PGKcPGKgKmBPGA13 + PGA3c[t]/PGKcPGKgKmPGA3)), PGKg -> (PGKgPGKgVmax*ADPg[t]*BPGA13g[t]*(1 - (ATPg[t]*PGA3g[t])/(PGKgPGKgKeq*ADPg[t]*BPGA13g[t])))/(cytosol*PGKgPGKgKmADP*PGKgPGKgKmBPGA13*(1 + ADPg[t]/PGKgPGKgKmADP + ATPg[t]/PGKgPGKgKmATP)*(1 + BPGA13g[t]/PGKgPGKgKmBPGA13 + PGA3g[t]/PGKgPGKgKmPGA3)), PYKc -> (PYKcPYKcVmax*ADPc[t]*(PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn*(1 - (ATPc[t]*Pyrc[t])/(PYKcPYKcKeq*ADPc[t]*PEPc[t])))/(cytosol*PYKcPYKcKmADP*(1 + ADPc[t]/PYKcPYKcKmADP + ATPc[t]/PYKcPYKcKmATP)*(1 + (PEPc[t]/(PYKcPYKcKmPEP*(1 + ADPc[t]/PYKcPYKcKiADP + ATPc[t]/PYKcPYKcKiATP)))^PYKcPYKcn + Pyrc[t]/PYKcPYKcKmPyr)), PyrTc -> (PyrTcPyrTcVmax*Pyrc[t])/(cytosol*PyrTcPyrTcKmPyr*(1 + Pyrc[t]/PyrTcPyrTcKmPyr)), TPIc -> (TPIcTPIcVmax*DHAPc[t]*(1 - GA3Pc[t]/(TPIcTPIgKeq*DHAPc[t])))/(cytosol*TPIcTPIgKmDHAP*(1 + DHAPc[t]/TPIcTPIgKmDHAP + GA3Pc[t]/TPIcTPIgKmGA3P)), TPIg -> (TPIgTPIgVmax*DHAPg[t]*(1 - GA3Pg[t]/(TPIgTPIgKeq*DHAPg[t])))/(cytosol*TPIgTPIgKmDHAP*(1 + DHAPg[t]/TPIgTPIgKmDHAP + GA3Pg[t]/TPIgTPIgKmGA3P))
};

parameters = {
ADPtgADPtgk -> 1000.0, AKcAKck1 -> 442.0, AKcAKck2 -> 1000.0, AKgAKgk1 -> 442.0, AKgAKgk2 -> 1000.0, ALDcALDcVmax -> 0.056, ALDcALDgKeq -> 0.084, ALDcALDgKiADP -> 1.51, ALDcALDgKiAMP -> 3.65, ALDcALDgKiATP -> 0.68, ALDcALDgKiGA3P -> 0.098, ALDcALDgKmDHAP -> 0.015, ALDcALDgKmFru16BP -> 0.009, ALDcALDgKmGA3P -> 0.067, ALDgALDgKeq -> 0.084, ALDgALDgKiADP -> 1.51, ALDgALDgKiAMP -> 3.65, ALDgALDgKiATP -> 0.68, ALDgALDgKiGA3P -> 0.098, ALDgALDgKmDHAP -> 0.015, ALDgALDgKmFru16BP -> 0.009, ALDgALDgKmGA3P -> 0.067, ALDgALDgVmax -> 560.0, AMPtgAMPtgk -> 1000.0, ATPtgATPtgk -> 1000.0, ATPucATPuck -> 50.0, BPGA13tgBPGA13tgk -> 1000.0, DHAPtgDHAPtgk -> 1000.0, ENOcENOcKeq -> 4.12, ENOcENOcKmPEP -> 0.24, ENOcENOcKmPGA2 -> 0.054, ENOcENOcVmax -> 598.0, EXTERNAL -> 0.0, F16BPtgF16BPtgk -> 1000.0, F6PtgF6Ptgk -> 1000.0, G3PDHcG3PDHcVmax -> 0.047, G3PDHcG3PDHgKeq -> 17085.0, G3PDHcG3PDHgKmDHAP -> 0.1, G3PDHcG3PDHgKmGly3P -> 2.0, G3PDHcG3PDHgKmNAD -> 0.4, G3PDHcG3PDHgKmNADH -> 0.01, G3PDHgG3PDHgKeq -> 17085.0, G3PDHgG3PDHgKmDHAP -> 0.1, G3PDHgG3PDHgKmGly3P -> 2.0, G3PDHgG3PDHgKmNAD -> 0.4, G3PDHgG3PDHgKmNADH -> 0.01, G3PDHgG3PDHgVmax -> 465.0, G6PPcG6PPcKeq -> 268.0, G6PPcG6PPcKmGlc -> 5.0, G6PPcG6PPcKmGlc6P -> 5.0, G6PPcG6PPcVmax -> 1.0, G6PtgG6Ptgk -> 1000.0, GA3PtgGA3Ptgk -> 1000.0, GAPDHcGAPDHcKeq -> 0.066, GAPDHcGAPDHcKmBPGA13 -> 0.14, GAPDHcGAPDHcKmGA3P -> 0.17, GAPDHcGAPDHcKmNAD -> 0.04, GAPDHcGAPDHcKmNADH -> 0.007, GAPDHcGAPDHcVmax -> 144.18, GAPDHgGAPDHgKeq -> 0.066, GAPDHgGAPDHgKmBPGA13 -> 0.1, GAPDHgGAPDHgKmGA3P -> 0.15, GAPDHgGAPDHgKmNAD -> 0.45, GAPDHgGAPDHgKmNADH -> 0.02, GAPDHgGAPDHgVmax -> 576.72, GKcGKcVmax -> 0.02, GKcGKgKeq -> 0.000837, GKcGKgKmADP -> 0.56, GKcGKgKmATP -> 0.24, GKcGKgKmGly -> 0.44, GKcGKgKmGly3P -> 3.83, GKgGKgKeq -> 0.000837, GKgGKgKmADP -> 0.56, GKgGKgKmATP -> 0.24, GKgGKgKmGly -> 0.44, GKgGKgKmGly3P -> 3.83, GKgGKgVmax -> 200.0, GPOcGPOcKiSHAM -> 0.021, GPOcGPOcKmGly3P -> 1.7, GPOcGPOcVmax -> 368.0, GlcTcGlcTcKmGlc -> 1.1, GlcTcGlcTcVmax -> 111.7, GlcTcGlcTcalpha -> 0.75, GlcTgGlcTgk -> 250000.0, Gly3PtgGly3Ptgk -> 1000.0, GlyTcGlyTcKmGlye -> 0.17, GlyTcGlyTcVrmax -> 85.0, GlyTcGlyTck -> 8.5, GlyTcGlyTcx -> 1.0, GlyTcGlyTcy -> 1.0, GlytgGlytgk -> 1000.0, HXKcHXKcVmax -> 0.1929, HXKcHXKgKeq -> 759.0, HXKcHXKgKmADP -> 0.126, HXKcHXKgKmATP -> 0.116, HXKcHXKgKmGlc -> 0.1, HXKcHXKgKmGlc6P -> 12.0, HXKgHXKgKeq -> 759.0, HXKgHXKgKmADP -> 0.126, HXKgHXKgKmATP -> 0.116, HXKgHXKgKmGlc -> 0.1, HXKgHXKgKmGlc6P -> 12.0, HXKgHXKgVmax -> 1929.0, O2c -> 1.0, PFKcPFKcVmax -> 0.1708, PFKcPFKgKeq -> 1035.0, PFKcPFKgKi1 -> 15.8, PFKcPFKgKi2 -> 10.7, PFKcPFKgKmADP -> 1.0, PFKcPFKgKmATP -> 0.0648, PFKcPFKgKmFru6P -> 0.999, PFKcPFKgKsATP -> 0.0393, PFKgPFKgKeq -> 1035.0, PFKgPFKgKi1 -> 15.8, PFKgPFKgKi2 -> 10.7, PFKgPFKgKmADP -> 1.0, PFKgPFKgKmATP -> 0.0648, PFKgPFKgKmFru6P -> 0.999, PFKgPFKgKsATP -> 0.0393, PFKgPFKgVmax -> 1708.0, PGA3TgPGA3Tgk -> 250.0, PGAMcPGAMcKeq -> 0.17, PGAMcPGAMcKmPGA2 -> 0.16, PGAMcPGAMcKmPGA3 -> 0.15, PGAMcPGAMcVmax -> 225.0, PGIcPGIcVmax -> 0.1305, PGIcPGIgKeq -> 0.46, PGIcPGIgKmFru6P -> 0.12, PGIcPGIgKmGlc6P -> 0.4, PGIgPGIgKeq -> 0.46, PGIgPGIgKmFru6P -> 0.12, PGIgPGIgKmGlc6P -> 0.4, PGIgPGIgVmax -> 1305.0, PGKcPGKcVmax -> 0.0286, PGKcPGKgKeq -> 3376.0, PGKcPGKgKmADP -> 0.1, PGKcPGKgKmATP -> 0.29, PGKcPGKgKmBPGA13 -> 0.003, PGKcPGKgKmPGA3 -> 1.62, PGKgPGKgKeq -> 3376.0, PGKgPGKgKmADP -> 0.1, PGKgPGKgKmATP -> 0.29, PGKgPGKgKmBPGA13 -> 0.003, PGKgPGKgKmPGA3 -> 1.62, PGKgPGKgVmax -> 2862.0, PYKcPYKcKeq -> 10800.0, PYKcPYKcKiADP -> 0.64, PYKcPYKcKiATP -> 0.57, PYKcPYKcKmADP -> 0.114, PYKcPYKcKmATP -> 15.0, PYKcPYKcKmPEP -> 0.34, PYKcPYKcKmPyr -> 50.0, PYKcPYKcVmax -> 1020.0, PYKcPYKcn -> 2.5, PyrTcPyrTcKmPyr -> 1.96, PyrTcPyrTcVmax -> 230.0, SHAM -> 0.0, TPIcTPIcVmax -> 0.09993, TPIcTPIgKeq -> 0.046, TPIcTPIgKmDHAP -> 1.2, TPIcTPIgKmGA3P -> 0.25, TPIgTPIgKeq -> 0.046, TPIgTPIgKmDHAP -> 1.2, TPIgTPIgKmGA3P -> 0.25, TPIgTPIgVmax -> 999.3, cytosol -> 5.4549, default -> 1.0, glycosome -> 0.2451, Glce -> 5.0, Glye -> 1*^-06, NADHc -> 0.0005, NADc -> 0.36, Pic -> 0.0, Pig -> 0.0, Pyre -> 0.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
ADPc'[t] == 1.0*ADPtg +1.0*ATPuc +1.0*HXKc +1.0*PFKc -1.0*GKc -2.0*AKc -1.0*PYKc -1.0*PGKc, ADPg'[t] == 22.2558*PFKg +22.2558*HXKg -22.2558*PGKg -44.5116*AKg -22.2558*GKg -22.2558*ADPtg, AMPc'[t] == 1.0*AKc +1.0*AMPtg , AMPg'[t] == 22.2558*AKg -22.2558*AMPtg, ATPc'[t] == 1.0*GKc +1.0*AKc +1.0*PYKc +1.0*PGKc +1.0*ATPtg -1.0*ATPuc -1.0*HXKc -1.0*PFKc, ATPg'[t] == 22.2558*PGKg +22.2558*AKg +22.2558*GKg -22.2558*PFKg -22.2558*HXKg -22.2558*ATPtg, BPGA13c'[t] == 1.0*GAPDHc +1.0*BPGA13tg -1.0*PGKc, BPGA13g'[t] == 22.2558*GAPDHg -22.2558*PGKg -22.2558*BPGA13tg, DHAPc'[t] == 1.0*GPOc +1.0*ALDc -1.0*G3PDHc -1.0*TPIc -1.0*DHAPtg, DHAPg'[t] == 22.2558*ALDg +22.2558*DHAPtg -22.2558*G3PDHg -22.2558*TPIg, Fru16BPc'[t] == 1.0*F16BPtg +1.0*PFKc -1.0*ALDc, Fru16BPg'[t] == 22.2558*PFKg -22.2558*F16BPtg -22.2558*ALDg, Fru6Pc'[t] == 1.0*F6Ptg +1.0*PGIc -1.0*PFKc, Fru6Pg'[t] == 22.2558*PGIg -22.2558*F6Ptg -22.2558*PFKg, GA3Pc'[t] == 1.0*ALDc +1.0*TPIc +1.0*GA3Ptg -1.0*GAPDHc, GA3Pg'[t] == 22.2558*ALDg +22.2558*TPIg -22.2558*GAPDHg -22.2558*GA3Ptg, Glc6Pc'[t] == 1.0*G6Ptg +1.0*HXKc -1.0*PGIc -1.0*G6PPc, Glc6Pg'[t] == 22.2558*HXKg -22.2558*G6Ptg -22.2558*PGIg, Glcc'[t] == 1.0*GlcTc +1.0*G6PPc -1.0*GlcTg -1.0*HXKc, Glcg'[t] == 22.2558*GlcTg -22.2558*HXKg, Gly3Pc'[t] == 1.0*G3PDHc +1.0*Gly3Ptg -1.0*GPOc -1.0*GKc, Gly3Pg'[t] == 22.2558*G3PDHg -22.2558*GKg -22.2558*Gly3Ptg, Glyc'[t] == 1.0*GKc +1.0*Glytg -1.0*GlyTc, Glyg'[t] == 22.2558*GKg -22.2558*Glytg, NADHg'[t] == 22.2558*GAPDHg -22.2558*G3PDHg, NADg'[t] == 22.2558*G3PDHg -22.2558*GAPDHg, PEPc'[t] == 1.0*ENOc -1.0*PYKc, PGA2c'[t] == 1.0*PGAMc -1.0*ENOc, PGA3c'[t] == 1.0*PGA3Tg +1.0*PGKc -1.0*PGAMc, PGA3g'[t] == 22.2558*PGKg -22.2558*PGA3Tg, Pyrc'[t] == 1.0*PYKc -1.0*PyrTc
};
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]}]