C
C There are a total of 19 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 72 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is G_o in component G_o (millimolar).
C STATES(1) is G in component G (millimolar).
C ALGBRC(4) is V_IN in component V_IN (flux).
C ALGBRC(5) is V_HK in component V_HK (flux).
C STATES(2) is G6P in component G6P (millimolar).
C ALGBRC(6) is V_PFK in component V_PFK (flux).
C ALGBRC(9) is V_G6PDH in component V_G6PDH (flux).
C STATES(3) is FDP in component FDP (millimolar).
C ALGBRC(8) is V_ALD in component V_ALD (flux).
C STATES(4) is G3P in component G3P (millimolar).
C ALGBRC(10) is V_GAPDH in component V_GAPDH (flux).
C STATES(5) is DPG in component DPG (millimolar).
C ALGBRC(11) is V_PGK in component V_PGK (flux).
C STATES(6) is PEP in component PEP (millimolar).
C ALGBRC(16) is V_PK in component V_PK (flux).
C STATES(7) is Py in component Py (millimolar).
C ALGBRC(17) is V_TCA in component V_TCA (flux).
C ALGBRC(18) is V_ADH in component V_ADH (flux).
C STATES(8) is ATP in component ATP (millimolar).
C CONSTS(2) is PO in component ATP (dimensionless).
C ALGBRC(19) is V_ATPase in component V_ATPase (flux).
C ALGBRC(1) is ADP in component ADP (millimolar).
C CONSTS(3) is Cn in component ADP (millimolar).
C CONSTS(4) is AMP in component AMP (millimolar).
C CONSTS(5) is GTP in component GTP (millimolar).
C CONSTS(6) is GDP in component GDP (millimolar).
C CONSTS(7) is H in component H (millimolar).
C CONSTS(8) is NADP in component NADP (millimolar).
C CONSTS(9) is NADH in component NADH (millimolar).
C CONSTS(10) is NAD in component NAD (millimolar).
C ALGBRC(3) is CD in component CD (millimolar).
C CONSTS(11) is CMTP in component CD (millimolar).
C STATES(9) is CT in component CT (millimolar).
C STATES(10) is CP in component CP (millimolar).
C CONSTS(12) is kpol in component CT (third_order_rate_constant).
C CONSTS(13) is kf in component CT (first_order_rate_constant).
C CONSTS(14) is kb in component CT (second_order_rate_constant).
C CONSTS(15) is kdp in component CP (first_order_rate_constant).
C STATES(11) is PKp in component PKp (millimolar).
C CONSTS(16) is kp2 in component PKp (second_order_rate_constant).
C CONSTS(17) is kp3 in component PKp (first_order_rate_constant).
C CONSTS(18) is k4 in component PKp (second_order_rate_constant).
C ALGBRC(2) is PKt in component PKt (millimolar).
C CONSTS(19) is C_PK in component PKt (millimolar).
C CONSTS(20) is Ke_in in component V_IN (millimolar).
C CONSTS(21) is KG_in in component V_IN (millimolar).
C CONSTS(22) is V_IN_max in component V_IN (flux).
C CONSTS(23) is KG_m in component V_HK (millimolar).
C CONSTS(24) is KG_s in component V_HK (millimolar).
C CONSTS(25) is KATP_m in component V_HK (millimolar).
C CONSTS(26) is V_HK_max in component V_HK (flux).
C CONSTS(27) is KG6P_r in component V_PFK (millimolar).
C CONSTS(28) is KATP_r in component V_PFK (millimolar).
C CONSTS(29) is KAMP_r in component V_PFK (millimolar).
C CONSTS(30) is cATP in component V_PFK (dimensionless).
C CONSTS(31) is cAMP in component V_PFK (dimensionless).
C CONSTS(32) is cG6P in component V_PFK (dimensionless).
C CONSTS(33) is Lo in component V_PFK (dimensionless).
C CONSTS(34) is gr in component V_PFK (dimensionless).
C CONSTS(35) is n1 in component V_PFK (dimensionless).
C CONSTS(36) is V_PFK_max in component V_PFK (flux).
C ALGBRC(7) is TUB in component V_G6PDH (millimolar).
C CONSTS(37) is KG6P in component V_G6PDH (millimolar).
C CONSTS(38) is KNADP in component V_G6PDH (millimolar).
C CONSTS(39) is KNADP_ in component V_G6PDH (millimolar).
C CONSTS(40) is KTUB in component V_G6PDH (millimolar).
C CONSTS(41) is V_G6PDH_max in component V_G6PDH (flux).
C CONSTS(42) is V_G6PDH_max_II in component V_G6PDH (flux).
C CONSTS(43) is KG3P_m in component V_ALD (millimolar).
C CONSTS(44) is KFDP_m in component V_ALD (millimolar).
C CONSTS(45) is V_ALD_max in component V_ALD (flux).
C CONSTS(46) is V_ALD_max_r in component V_ALD (flux).
C CONSTS(47) is K1 in component V_GAPDH (millimolar).
C CONSTS(48) is K2 in component V_GAPDH (millimolar).
C CONSTS(49) is K3 in component V_GAPDH (millimolar).
C CONSTS(50) is KG3P in component V_GAPDH (millimolar).
C CONSTS(51) is KNAD in component V_GAPDH (millimolar).
C CONSTS(52) is KNADH_i in component V_GAPDH (millimolar).
C CONSTS(53) is V_GAPDH_max in component V_GAPDH (flux).
C CONSTS(54) is KDPG_m in component V_PGK (millimolar).
C CONSTS(55) is V_PGK_max in component V_PGK (flux).
C ALGBRC(12) is R in component V_PK (dimensionless).
C ALGBRC(13) is T in component V_PK (dimensionless).
C CONSTS(56) is KpH in component V_PK (millimolar).
C CONSTS(57) is KPEP_r in component V_PK (millimolar).
C CONSTS(58) is KADP_r in component V_PK (millimolar).
C CONSTS(59) is KFDP_r in component V_PK (millimolar).
C CONSTS(60) is cADP in component V_PK (dimensionless).
C CONSTS(61) is cFDP in component V_PK (dimensionless).
C CONSTS(62) is cPEP in component V_PK (dimensionless).
C CONSTS(63) is Lo_PK in component V_PK (dimensionless).
C CONSTS(64) is gr_PK in component V_PK (dimensionless).
C CONSTS(65) is gt_PK in component V_PK (dimensionless).
C ALGBRC(15) is n in component V_PK (dimensionless).
C ALGBRC(14) is V_PK_max in component V_PK (flux).
C CONSTS(66) is V_PKt_max in component V_PK (flux).
C CONSTS(67) is V_PKp_max in component V_PK (flux).
C CONSTS(68) is KPy_m in component V_TCA (millimolar).
C CONSTS(69) is V_TCA_max in component V_TCA (flux).
C CONSTS(70) is KPy__m in component V_ADH (millimolar).
C CONSTS(71) is V_ADH_max in component V_ADH (flux).
C CONSTS(72) is KATP in component V_ATPase (first_order_rate_constant).
C RATES(1) is d/dt G in component G (millimolar).
C RATES(2) is d/dt G6P in component G6P (millimolar).
C RATES(3) is d/dt FDP in component FDP (millimolar).
C RATES(4) is d/dt G3P in component G3P (millimolar).
C RATES(5) is d/dt DPG in component DPG (millimolar).
C RATES(6) is d/dt PEP in component PEP (millimolar).
C RATES(7) is d/dt Py in component Py (millimolar).
C RATES(8) is d/dt ATP in component ATP (millimolar).
C RATES(9) is d/dt CT in component CT (millimolar).
C RATES(10) is d/dt CP in component CP (millimolar).
C RATES(11) is d/dt PKp in component PKp (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      STATES(1) = 0.01
      STATES(2) = 0.01
      STATES(3) = 0.01
      STATES(4) = 0.01
      STATES(5) = 0.01
      STATES(6) = 0.01
      STATES(7) = 0.01
      STATES(8) = 1.4
      CONSTS(2) = 4
      CONSTS(3) = 9
      CONSTS(4) = 0.5
      CONSTS(5) = 0.95
      CONSTS(6) = 0.05
      CONSTS(7) = 3.2e-8
      CONSTS(8) = 1
      CONSTS(9) = 0.01
      CONSTS(10) = 1
      CONSTS(11) = 0.9
      STATES(9) = 0.2
      STATES(10) = 1.2
      CONSTS(12) = 10
      CONSTS(13) = 3
      CONSTS(14) = 2.5
      CONSTS(15) = 0.0025
      STATES(11) = 0.005
      CONSTS(16) = 10
      CONSTS(17) = 0.05
      CONSTS(18) = 0.02
      CONSTS(19) = 0.01
      CONSTS(20) = 12
      CONSTS(21) = 0.001
      CONSTS(22) = 10
      CONSTS(23) = 0.11
      CONSTS(24) = 0.0062
      CONSTS(25) = 0.1
      CONSTS(26) = 13
      CONSTS(27) = 1
      CONSTS(28) = 0.06
      CONSTS(29) = 0.025
      CONSTS(30) = 1
      CONSTS(31) = 0.019
      CONSTS(32) = 0.0005
      CONSTS(33) = 25000
      CONSTS(34) = 10
      CONSTS(35) = 2
      CONSTS(36) = 30
      CONSTS(37) = 0.05
      CONSTS(38) = 0.05
      CONSTS(39) = 0.05
      CONSTS(40) = 0.4
      CONSTS(41) = 1.6
      CONSTS(42) = 1
      CONSTS(43) = 20
      CONSTS(44) = 0.5
      CONSTS(45) = 2.5
      CONSTS(46) = 1
      CONSTS(47) = 1.1
      CONSTS(48) = 1.5
      CONSTS(49) = 2.5
      CONSTS(50) = 0.0025
      CONSTS(51) = 0.18
      CONSTS(52) = 0.0003
      CONSTS(53) = 10
      CONSTS(54) = 0.002
      CONSTS(55) = 3
      CONSTS(56) = 9.5e-9
      CONSTS(57) = 1
      CONSTS(58) = 0.06
      CONSTS(59) = 0.025
      CONSTS(60) = 1
      CONSTS(61) = 0.01
      CONSTS(62) = 0.02
      CONSTS(63) = 1000
      CONSTS(64) = 0.1
      CONSTS(65) = 1
      CONSTS(66) = 25
      CONSTS(67) = 50
      CONSTS(68) = 0.329
      CONSTS(69) = 10
      CONSTS(70) = 0.169
      CONSTS(71) = 0.5
      CONSTS(72) = 5
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(10) =  CONSTS(12)*STATES(9)*STATES(10) ** 2.00000 -  CONSTS(15)*STATES(10)
      ALGBRC(2) = CONSTS(19) - STATES(11)
      RATES(11) =  0.100000*CONSTS(16)*ALGBRC(2)*STATES(10) - ( CONSTS(17)*STATES(11)+ CONSTS(18)*STATES(11)*CONSTS(5))
      ALGBRC(3) = CONSTS(11) - (STATES(9)+STATES(10))
      RATES(9) = - ( CONSTS(12)*STATES(9)*STATES(10) ** 2.00000+ CONSTS(13)*ALGBRC(3)+ CONSTS(14)*STATES(9)*CONSTS(6))
      ALGBRC(4) =  CONSTS(22)*(CONSTS(1)/( (CONSTS(21)+CONSTS(1))*(1.00000+STATES(2)/CONSTS(20))) - STATES(1)/( (CONSTS(21)+STATES(1))*(1.00000+STATES(2)/CONSTS(20))))
      ALGBRC(5) = ( CONSTS(26)*1.00000)/(1.00000+( CONSTS(24)*CONSTS(25))/( STATES(1)*STATES(8))+CONSTS(23)/STATES(1)+CONSTS(25)/STATES(8))
      RATES(1) = ALGBRC(4) - ALGBRC(5)
      ALGBRC(6) = ( (( (( CONSTS(36)*CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28))*1.00000+STATES(2)/CONSTS(27)+STATES(8)/CONSTS(28)+( (( CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28) ** CONSTS(35) - 1.00000)/(1.00000+STATES(2)/CONSTS(27)+STATES(8)/CONSTS(28)+( (( CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28) ** CONSTS(35)+ CONSTS(33)*(1.00000+( CONSTS(31)*CONSTS(4))/CONSTS(29))/(1.00000+CONSTS(4)/CONSTS(29)) ** CONSTS(35)*1.00000+( CONSTS(32)*STATES(2))/CONSTS(27)+( CONSTS(30)*STATES(8))/CONSTS(28)+( (( CONSTS(34)*CONSTS(32)*STATES(2))/CONSTS(27))*CONSTS(30)*STATES(8))/CONSTS(28) ** CONSTS(35))
      ALGBRC(8) = (( CONSTS(45)*STATES(3))/CONSTS(44) - ( CONSTS(46)*STATES(4))/CONSTS(43))/(1.00000+STATES(3)/CONSTS(44)+STATES(4)/CONSTS(43))
      RATES(3) = ALGBRC(6) - ALGBRC(8)
      ALGBRC(7) = STATES(9)+ALGBRC(3)
      ALGBRC(9) = CONSTS(41)/(( CONSTS(37)*CONSTS(38))/( STATES(2)*CONSTS(8))+CONSTS(37)/STATES(2)+CONSTS(38)/CONSTS(8)+1.00000)+CONSTS(42)/(( CONSTS(37)*CONSTS(39)*CONSTS(40))/( STATES(2)*CONSTS(8)*ALGBRC(7))+( CONSTS(37)*CONSTS(39))/( STATES(2)*CONSTS(8))+( CONSTS(39)*CONSTS(40))/( CONSTS(8)*ALGBRC(7))+( CONSTS(37)*CONSTS(40))/( STATES(2)*ALGBRC(7))+CONSTS(40)/ALGBRC(7)+CONSTS(37)/STATES(2)+CONSTS(39)/CONSTS(8)+1.00000)
      RATES(2) = ALGBRC(5) - (ALGBRC(6)+ALGBRC(9))
      ALGBRC(1) = CONSTS(3) - (STATES(8)+CONSTS(4))
      ALGBRC(10) = CONSTS(53)/(1.00000+CONSTS(50)/STATES(4)+ (CONSTS(51)/CONSTS(10))*(1.00000+CONSTS(4)/CONSTS(47)+ALGBRC(1)/CONSTS(48)+STATES(8)/CONSTS(49))+ (( CONSTS(50)*CONSTS(51))/( STATES(4)*CONSTS(10)))*(1.00000+CONSTS(9)/CONSTS(52))+1.00000+CONSTS(4)/CONSTS(47)+ALGBRC(1)/CONSTS(48)+STATES(8)/CONSTS(49))
      RATES(4) =  2.00000*ALGBRC(8) - ALGBRC(10)
      ALGBRC(11) = ( CONSTS(55)*STATES(5))/(CONSTS(54)+STATES(5))
      RATES(5) = ALGBRC(10) - ALGBRC(11)
      ALGBRC(12) = 1.00000+STATES(6)/CONSTS(57)+ALGBRC(1)/CONSTS(58)+( (( CONSTS(64)*STATES(6))/CONSTS(57))*ALGBRC(1))/CONSTS(58)
      ALGBRC(13) = 1.00000+( CONSTS(62)*STATES(6))/CONSTS(57)+( CONSTS(60)*ALGBRC(1))/CONSTS(58)+( (( CONSTS(65)*CONSTS(62)*STATES(6))/CONSTS(57))*CONSTS(60)*ALGBRC(1))/CONSTS(58)
      ALGBRC(15) = 4.00000+STATES(11)/CONSTS(19)
      ALGBRC(14) = CONSTS(66)+( (CONSTS(67) - CONSTS(66))*STATES(11))/CONSTS(19)
      ALGBRC(16) = ( (ALGBRC(14)/(1.00000+CONSTS(56)/CONSTS(7)))*( CONSTS(64)*(STATES(6)/CONSTS(57))*(ALGBRC(1)/CONSTS(58))*ALGBRC(12) ** ALGBRC(15) - 1.00000+ CONSTS(63)*(1.00000+( CONSTS(61)*STATES(3))/CONSTS(59))/(1.00000+STATES(3)/CONSTS(59)) ** ALGBRC(15)*(STATES(3)/CONSTS(59))*CONSTS(65)*(( CONSTS(62)*STATES(6))/CONSTS(57))*(( CONSTS(60)*ALGBRC(1))/CONSTS(58))*ALGBRC(13) ** ALGBRC(15) - 1.00000))/(ALGBRC(12) ** ALGBRC(15)+ CONSTS(63)*(1.00000+( CONSTS(61)*STATES(3))/CONSTS(59))/(1.00000+STATES(3)/CONSTS(59)) ** ALGBRC(15)*ALGBRC(13) ** ALGBRC(15))
      RATES(6) = ALGBRC(11) - ALGBRC(16)
      ALGBRC(17) = ( CONSTS(69)*STATES(7) ** 2.00000)/(CONSTS(68) ** 2.00000+STATES(7) ** 2.00000)
      ALGBRC(18) = ( CONSTS(71)*STATES(7))/(CONSTS(70)+STATES(7))
      RATES(7) = ALGBRC(16) - (ALGBRC(17)+ALGBRC(18))
      ALGBRC(19) =  CONSTS(72)*STATES(8)
      RATES(8) = (ALGBRC(11)+ALGBRC(16)+ CONSTS(2)*ALGBRC(17)) - (ALGBRC(5)+ALGBRC(6)+ALGBRC(19))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = CONSTS(19) - STATES(11)
      ALGBRC(3) = CONSTS(11) - (STATES(9)+STATES(10))
      ALGBRC(4) =  CONSTS(22)*(CONSTS(1)/( (CONSTS(21)+CONSTS(1))*(1.00000+STATES(2)/CONSTS(20))) - STATES(1)/( (CONSTS(21)+STATES(1))*(1.00000+STATES(2)/CONSTS(20))))
      ALGBRC(5) = ( CONSTS(26)*1.00000)/(1.00000+( CONSTS(24)*CONSTS(25))/( STATES(1)*STATES(8))+CONSTS(23)/STATES(1)+CONSTS(25)/STATES(8))
      ALGBRC(6) = ( (( (( CONSTS(36)*CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28))*1.00000+STATES(2)/CONSTS(27)+STATES(8)/CONSTS(28)+( (( CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28) ** CONSTS(35) - 1.00000)/(1.00000+STATES(2)/CONSTS(27)+STATES(8)/CONSTS(28)+( (( CONSTS(34)*STATES(2))/CONSTS(27))*STATES(8))/CONSTS(28) ** CONSTS(35)+ CONSTS(33)*(1.00000+( CONSTS(31)*CONSTS(4))/CONSTS(29))/(1.00000+CONSTS(4)/CONSTS(29)) ** CONSTS(35)*1.00000+( CONSTS(32)*STATES(2))/CONSTS(27)+( CONSTS(30)*STATES(8))/CONSTS(28)+( (( CONSTS(34)*CONSTS(32)*STATES(2))/CONSTS(27))*CONSTS(30)*STATES(8))/CONSTS(28) ** CONSTS(35))
      ALGBRC(8) = (( CONSTS(45)*STATES(3))/CONSTS(44) - ( CONSTS(46)*STATES(4))/CONSTS(43))/(1.00000+STATES(3)/CONSTS(44)+STATES(4)/CONSTS(43))
      ALGBRC(7) = STATES(9)+ALGBRC(3)
      ALGBRC(9) = CONSTS(41)/(( CONSTS(37)*CONSTS(38))/( STATES(2)*CONSTS(8))+CONSTS(37)/STATES(2)+CONSTS(38)/CONSTS(8)+1.00000)+CONSTS(42)/(( CONSTS(37)*CONSTS(39)*CONSTS(40))/( STATES(2)*CONSTS(8)*ALGBRC(7))+( CONSTS(37)*CONSTS(39))/( STATES(2)*CONSTS(8))+( CONSTS(39)*CONSTS(40))/( CONSTS(8)*ALGBRC(7))+( CONSTS(37)*CONSTS(40))/( STATES(2)*ALGBRC(7))+CONSTS(40)/ALGBRC(7)+CONSTS(37)/STATES(2)+CONSTS(39)/CONSTS(8)+1.00000)
      ALGBRC(1) = CONSTS(3) - (STATES(8)+CONSTS(4))
      ALGBRC(10) = CONSTS(53)/(1.00000+CONSTS(50)/STATES(4)+ (CONSTS(51)/CONSTS(10))*(1.00000+CONSTS(4)/CONSTS(47)+ALGBRC(1)/CONSTS(48)+STATES(8)/CONSTS(49))+ (( CONSTS(50)*CONSTS(51))/( STATES(4)*CONSTS(10)))*(1.00000+CONSTS(9)/CONSTS(52))+1.00000+CONSTS(4)/CONSTS(47)+ALGBRC(1)/CONSTS(48)+STATES(8)/CONSTS(49))
      ALGBRC(11) = ( CONSTS(55)*STATES(5))/(CONSTS(54)+STATES(5))
      ALGBRC(12) = 1.00000+STATES(6)/CONSTS(57)+ALGBRC(1)/CONSTS(58)+( (( CONSTS(64)*STATES(6))/CONSTS(57))*ALGBRC(1))/CONSTS(58)
      ALGBRC(13) = 1.00000+( CONSTS(62)*STATES(6))/CONSTS(57)+( CONSTS(60)*ALGBRC(1))/CONSTS(58)+( (( CONSTS(65)*CONSTS(62)*STATES(6))/CONSTS(57))*CONSTS(60)*ALGBRC(1))/CONSTS(58)
      ALGBRC(15) = 4.00000+STATES(11)/CONSTS(19)
      ALGBRC(14) = CONSTS(66)+( (CONSTS(67) - CONSTS(66))*STATES(11))/CONSTS(19)
      ALGBRC(16) = ( (ALGBRC(14)/(1.00000+CONSTS(56)/CONSTS(7)))*( CONSTS(64)*(STATES(6)/CONSTS(57))*(ALGBRC(1)/CONSTS(58))*ALGBRC(12) ** ALGBRC(15) - 1.00000+ CONSTS(63)*(1.00000+( CONSTS(61)*STATES(3))/CONSTS(59))/(1.00000+STATES(3)/CONSTS(59)) ** ALGBRC(15)*(STATES(3)/CONSTS(59))*CONSTS(65)*(( CONSTS(62)*STATES(6))/CONSTS(57))*(( CONSTS(60)*ALGBRC(1))/CONSTS(58))*ALGBRC(13) ** ALGBRC(15) - 1.00000))/(ALGBRC(12) ** ALGBRC(15)+ CONSTS(63)*(1.00000+( CONSTS(61)*STATES(3))/CONSTS(59))/(1.00000+STATES(3)/CONSTS(59)) ** ALGBRC(15)*ALGBRC(13) ** ALGBRC(15))
      ALGBRC(17) = ( CONSTS(69)*STATES(7) ** 2.00000)/(CONSTS(68) ** 2.00000+STATES(7) ** 2.00000)
      ALGBRC(18) = ( CONSTS(71)*STATES(7))/(CONSTS(70)+STATES(7))
      ALGBRC(19) =  CONSTS(72)*STATES(8)
      RETURN
      END