Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 21 entries in the algebraic variable array.
C There are a total of 21 entries in each of the rate and state variable arrays.
C There are a total of 116 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is GLY in component GLY (millimolar).
C ALGBRC(4) is flux_GP in component flux_GP (flux).
C STATES(2) is G1P in component G1P (millimolar).
C ALGBRC(5) is V_PGLM in component V_PGLM (flux).
C STATES(3) is G6P in component G6P (millimolar).
C ALGBRC(6) is V_PGI in component V_PGI (flux).
C STATES(4) is F6P in component F6P (millimolar).
C ALGBRC(10) is V_PFK in component V_PFK (flux).
C STATES(5) is FBP in component FBP (millimolar).
C ALGBRC(11) is V_ALD in component V_ALD (flux).
C STATES(6) is DHAP in component DHAP (millimolar).
C ALGBRC(12) is V_TPI in component V_TPI (flux).
C STATES(7) is GAP in component GAP (millimolar).
C ALGBRC(14) is V_GAPDH in component V_GAPDH (flux).
C STATES(8) is Thirteen_BPG in component Thirteen_BPG (millimolar).
C ALGBRC(15) is V_PGK in component V_PGK (flux).
C STATES(9) is three_PG in component three_PG (millimolar).
C ALGBRC(16) is V_PGM in component V_PGM (flux).
C STATES(10) is two_PG in component two_PG (millimolar).
C ALGBRC(17) is V_ENOL in component V_ENOL (flux).
C STATES(11) is PEP in component PEP (millimolar).
C ALGBRC(18) is V_PK in component V_PK (flux).
C STATES(12) is PYR in component PYR (millimolar).
C ALGBRC(19) is V_LDH in component V_LDH (flux).
C STATES(13) is LAC in component LAC (millimolar).
C ALGBRC(1) is output in component LAC (flux).
C STATES(14) is Pi in component Pi (millimolar).
C CONSTS(1) is V_ATPase in component V_ATPase (flux).
C STATES(15) is ADP in component ADP (millimolar).
C ALGBRC(21) is V_ADK in component V_ADK (flux).
C ALGBRC(20) is V_CK in component V_CK (flux).
C STATES(16) is ATP in component ATP (millimolar).
C STATES(17) is AMP in component AMP (millimolar).
C STATES(18) is PCr in component PCr (millimolar).
C STATES(19) is Cr in component Cr (millimolar).
C STATES(20) is NADH in component NADH (millimolar).
C STATES(21) is NAD in component NAD (millimolar).
C CONSTS(2) is frac_a in component flux_GP (dimensionless).
C CONSTS(3) is frac_b in component flux_GP (dimensionless).
C ALGBRC(2) is V_GPa in component V_GPa (flux).
C ALGBRC(3) is V_GPb in component V_GPb (flux).
C CONSTS(4) is Ki_GLY_f in component V_GPa (millimolar).
C CONSTS(5) is Ki_GLY_b in component V_GPa (millimolar).
C CONSTS(6) is K_GLY_f in component V_GPa (millimolar).
C CONSTS(7) is K_GLY_b in component V_GPa (millimolar).
C CONSTS(8) is Keq_GP in component V_GPa (dimensionless).
C CONSTS(9) is Ki_Pi in component V_GPa (millimolar).
C CONSTS(10) is K_Pi in component V_GPa (millimolar).
C CONSTS(11) is Ki_G1P in component V_GPa (millimolar).
C CONSTS(101) is V_max_r in component V_GPa (flux).
C CONSTS(12) is V_max_f in component V_GPa (flux).
C CONSTS(13) is Ki_GLY_f in component V_GPb (millimolar).
C CONSTS(14) is Ki_GLY_b in component V_GPb (millimolar).
C CONSTS(15) is K_GLY_f in component V_GPb (millimolar).
C CONSTS(16) is Keq_GP in component V_GPb (dimensionless).
C CONSTS(17) is nH in component V_GPb (dimensionless).
C CONSTS(18) is K_AMP in component V_GPb (millimolar_1_75).
C CONSTS(19) is K_GLY_b in component V_GPb (millimolar).
C CONSTS(20) is Ki_Pi in component V_GPb (millimolar).
C CONSTS(21) is K_Pi in component V_GPb (millimolar).
C CONSTS(22) is Ki_G1P in component V_GPb (millimolar).
C CONSTS(23) is K_G1P in component V_GPb (millimolar).
C CONSTS(102) is V_max_r in component V_GPb (flux).
C CONSTS(24) is V_max_f in component V_GPb (flux).
C CONSTS(25) is K_G1P in component V_PGLM (millimolar).
C CONSTS(26) is K_G6P in component V_PGLM (millimolar).
C CONSTS(27) is Keq_PGLM in component V_PGLM (dimensionless).
C CONSTS(103) is V_max_r in component V_PGLM (flux).
C CONSTS(28) is V_max_f in component V_PGLM (flux).
C CONSTS(29) is K_G6P in component V_PGI (millimolar).
C CONSTS(30) is K_F6P in component V_PGI (millimolar).
C CONSTS(31) is Keq_PGI in component V_PGI (dimensionless).
C CONSTS(32) is V_max_r in component V_PGI (flux).
C CONSTS(104) is V_max_f in component V_PGI (flux).
C CONSTS(33) is K_FBP in component V_PFK (millimolar).
C CONSTS(34) is K_FBP_ in component V_PFK (millimolar).
C CONSTS(35) is K_F6P in component V_PFK (millimolar).
C CONSTS(36) is K_ATP in component V_PFK (millimolar).
C CONSTS(37) is K_F6P_ in component V_PFK (millimolar).
C CONSTS(38) is K_ATP_ in component V_PFK (millimolar).
C CONSTS(39) is K_ADP in component V_PFK (millimolar).
C CONSTS(40) is K_ADP_ in component V_PFK (millimolar).
C CONSTS(41) is Ki_ATP in component V_PFK (millimolar).
C CONSTS(42) is Ka_AMP in component V_PFK (millimolar).
C CONSTS(105) is V_max_r in component V_PFK (flux).
C CONSTS(43) is V_max_f in component V_PFK (flux).
C CONSTS(106) is alpha in component V_PFK (dimensionless).
C ALGBRC(9) is L in component V_PFK (dimensionless).
C CONSTS(44) is Lo in component V_PFK (dimensionless).
C ALGBRC(7) is delta in component V_PFK (dimensionless).
C ALGBRC(8) is delta_ in component V_PFK (dimensionless).
C CONSTS(45) is d in component V_PFK (dimensionless).
C CONSTS(46) is var_e in component V_PFK (dimensionless).
C CONSTS(47) is K_GAP in component V_ALD (millimolar).
C CONSTS(48) is K_FBP in component V_ALD (millimolar).
C CONSTS(49) is K_DHAP in component V_ALD (millimolar).
C CONSTS(50) is Keq_ALD in component V_ALD (millimolar).
C CONSTS(107) is V_max_r in component V_ALD (flux).
C CONSTS(51) is V_max_f in component V_ALD (flux).
C CONSTS(52) is K_GAP in component V_TPI (millimolar).
C CONSTS(53) is K_DHAP in component V_TPI (millimolar).
C CONSTS(54) is Keq_TPI in component V_TPI (dimensionless).
C CONSTS(108) is V_max_r in component V_TPI (flux).
C CONSTS(55) is V_max_f in component V_TPI (flux).
C ALGBRC(13) is D_GAPDH in component V_GAPDH (dimensionless).
C CONSTS(56) is K_GAP in component V_GAPDH (millimolar).
C CONSTS(57) is K_Pi in component V_GAPDH (millimolar).
C CONSTS(58) is K_NAD in component V_GAPDH (millimolar).
C CONSTS(59) is K_NADH in component V_GAPDH (millimolar).
C CONSTS(60) is K_Thirteen_BPG in component V_GAPDH (millimolar).
C CONSTS(61) is Keq_GAPDH in component V_GAPDH (dimensionless).
C CONSTS(109) is V_max_r in component V_GAPDH (flux).
C CONSTS(62) is V_max_f in component V_GAPDH (flux).
C CONSTS(63) is K_three_PG in component V_PGK (millimolar).
C CONSTS(64) is K_ADP in component V_PGK (millimolar).
C CONSTS(65) is K_ATP in component V_PGK (millimolar).
C CONSTS(66) is K_Thirteen_BPG in component V_PGK (millimolar).
C CONSTS(67) is Keq_PGK in component V_PGK (dimensionless).
C CONSTS(68) is V_max_r in component V_PGK (flux).
C CONSTS(110) is V_max_f in component V_PGK (flux).
C CONSTS(69) is K_two_PG in component V_PGM (millimolar).
C CONSTS(70) is K_three_PG in component V_PGM (millimolar).
C CONSTS(71) is Keq_PGM in component V_PGM (dimensionless).
C CONSTS(111) is V_max_r in component V_PGM (flux).
C CONSTS(72) is V_max_f in component V_PGM (flux).
C CONSTS(73) is K_PEP in component V_ENOL (millimolar).
C CONSTS(74) is K_two_PG in component V_ENOL (millimolar).
C CONSTS(75) is Keq_ENOL in component V_ENOL (dimensionless).
C CONSTS(112) is V_max_r in component V_ENOL (flux).
C CONSTS(76) is V_max_f in component V_ENOL (flux).
C CONSTS(77) is K_PYR in component V_PK (millimolar).
C CONSTS(78) is K_ADP in component V_PK (millimolar).
C CONSTS(79) is K_ATP in component V_PK (millimolar).
C CONSTS(80) is K_PEP in component V_PK (millimolar).
C CONSTS(81) is Keq_PK in component V_PK (dimensionless).
C CONSTS(113) is V_max_r in component V_PK (flux).
C CONSTS(82) is V_max_f in component V_PK (flux).
C CONSTS(83) is K_LAC in component V_LDH (millimolar).
C CONSTS(84) is K_NADH in component V_LDH (millimolar).
C CONSTS(85) is K_NAD in component V_LDH (millimolar).
C CONSTS(86) is K_PYR in component V_LDH (millimolar).
C CONSTS(87) is Keq_LDH in component V_LDH (dimensionless).
C CONSTS(114) is V_max_r in component V_LDH (flux).
C CONSTS(88) is V_max_f in component V_LDH (flux).
C CONSTS(89) is Ki_ADP in component V_CK (millimolar).
C CONSTS(90) is K_Cr in component V_CK (millimolar).
C CONSTS(91) is K_PCr in component V_CK (millimolar).
C CONSTS(92) is Ki_PCr in component V_CK (millimolar).
C CONSTS(93) is Ki_ATP in component V_CK (millimolar).
C CONSTS(94) is Keq_CK in component V_CK (dimensionless).
C CONSTS(95) is V_max_r in component V_CK (flux).
C CONSTS(115) is V_max_f in component V_CK (flux).
C CONSTS(96) is K_ADP in component V_ADK (millimolar).
C CONSTS(97) is K_AMP in component V_ADK (millimolar).
C CONSTS(98) is K_ATP in component V_ADK (millimolar).
C CONSTS(99) is Keq_ADK in component V_ADK (dimensionless).
C CONSTS(116) is V_max_r in component V_ADK (flux).
C CONSTS(100) is V_max_f in component V_ADK (flux).
C RATES(1) is d/dt GLY in component GLY (millimolar).
C RATES(2) is d/dt G1P in component G1P (millimolar).
C RATES(3) is d/dt G6P in component G6P (millimolar).
C RATES(4) is d/dt F6P in component F6P (millimolar).
C RATES(5) is d/dt FBP in component FBP (millimolar).
C RATES(6) is d/dt DHAP in component DHAP (millimolar).
C RATES(7) is d/dt GAP in component GAP (millimolar).
C RATES(8) is d/dt Thirteen_BPG in component Thirteen_BPG (millimolar).
C RATES(9) is d/dt three_PG in component three_PG (millimolar).
C RATES(10) is d/dt two_PG in component two_PG (millimolar).
C RATES(11) is d/dt PEP in component PEP (millimolar).
C RATES(12) is d/dt PYR in component PYR (millimolar).
C RATES(13) is d/dt LAC in component LAC (millimolar).
C RATES(14) is d/dt Pi in component Pi (millimolar).
C RATES(15) is d/dt ADP in component ADP (millimolar).
C RATES(16) is d/dt ATP in component ATP (millimolar).
C RATES(17) is d/dt AMP in component AMP (millimolar).
C RATES(18) is d/dt PCr in component PCr (millimolar).
C RATES(19) is d/dt Cr in component Cr (millimolar).
C RATES(20) is d/dt NADH in component NADH (millimolar).
C RATES(21) is d/dt NAD in component NAD (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 112.0
STATES(2) = 0.0589
STATES(3) = 0.75
STATES(4) = 0.228
STATES(5) = 0.0723
STATES(6) = 0.0764
STATES(7) = 0.0355
STATES(8) = 0.065
STATES(9) = 0.052
STATES(10) = 0.005
STATES(11) = 0.0194
STATES(12) = 0.0994
STATES(13) = 1.3
STATES(14) = 4.1
CONSTS(1) = 600.0
STATES(15) = 0.013
STATES(16) = 8.2
STATES(17) = 0.99
STATES(18) = 34.67
STATES(19) = 35.0
STATES(20) = 0.001
STATES(21) = 0.5
CONSTS(2) = 0.4
CONSTS(3) = 0.6
CONSTS(4) = 2.0
CONSTS(5) = 2.0
CONSTS(6) = 1.7
CONSTS(7) = 0.15
CONSTS(8) = 0.42
CONSTS(9) = 4.7
CONSTS(10) = 4.0
CONSTS(11) = 10.1
CONSTS(12) = 20.0
CONSTS(13) = 15.0
CONSTS(14) = 4.4
CONSTS(15) = 1.7
CONSTS(16) = 0.42
CONSTS(17) = 1.75
CONSTS(18) = 1.9E-6
CONSTS(19) = 0.15
CONSTS(20) = 4.6
CONSTS(21) = 0.2
CONSTS(22) = 7.4
CONSTS(23) = 1.5
CONSTS(24) = 30.0
CONSTS(25) = 0.063
CONSTS(26) = 0.03
CONSTS(27) = 16.62
CONSTS(28) = 480.0
CONSTS(29) = 0.18
CONSTS(30) = 0.119
CONSTS(31) = 0.45
CONSTS(32) = 880.0
CONSTS(33) = 4.02
CONSTS(34) = 4.02
CONSTS(35) = 0.18
CONSTS(36) = 0.08
CONSTS(37) = 20.0
CONSTS(38) = 0.25
CONSTS(39) = 2.7
CONSTS(40) = 2.7
CONSTS(41) = 0.87
CONSTS(42) = 0.06
CONSTS(43) = 56.0
CONSTS(44) = 13.0
CONSTS(45) = 0.01
CONSTS(46) = 0.01
CONSTS(47) = 1.0
CONSTS(48) = 0.05
CONSTS(49) = 2.0
CONSTS(50) = 9.5E-5
CONSTS(51) = 104.0
CONSTS(52) = 0.32
CONSTS(53) = 0.61
CONSTS(54) = 0.052
CONSTS(55) = 12000.0
CONSTS(56) = 0.0025
CONSTS(57) = 0.29
CONSTS(58) = 0.09
CONSTS(59) = 0.0033
CONSTS(60) = 0.0008
CONSTS(61) = 0.089
CONSTS(62) = 1265.0
CONSTS(63) = 1.2
CONSTS(64) = 0.008
CONSTS(65) = 0.35
CONSTS(66) = 0.002
CONSTS(67) = 57109.0
CONSTS(68) = 1120.0
CONSTS(69) = 0.014
CONSTS(70) = 0.2
CONSTS(71) = 0.18
CONSTS(72) = 1120.0
CONSTS(73) = 0.37
CONSTS(74) = 0.1
CONSTS(75) = 0.49
CONSTS(76) = 192.0
CONSTS(77) = 7.05
CONSTS(78) = 0.3
CONSTS(79) = 1.13
CONSTS(80) = 0.08
CONSTS(81) = 10304.0
CONSTS(82) = 1440.0
CONSTS(83) = 17.0
CONSTS(84) = 0.002
CONSTS(85) = 0.849
CONSTS(86) = 0.335
CONSTS(87) = 16198.0
CONSTS(88) = 1920.0
CONSTS(89) = 0.135
CONSTS(90) = 3.8
CONSTS(91) = 1.11
CONSTS(92) = 3.9
CONSTS(93) = 3.5
CONSTS(94) = 233.0
CONSTS(95) = 500.0
CONSTS(96) = 0.35
CONSTS(97) = 0.32
CONSTS(98) = 0.27
CONSTS(99) = 2.21
CONSTS(100) = 880.0
CONSTS(101) = ( CONSTS(12)*CONSTS(7)*CONSTS(11))/( CONSTS(4)*CONSTS(10)*CONSTS(8))
CONSTS(102) = ( CONSTS(24)*CONSTS(14)*CONSTS(23))/( CONSTS(15)*CONSTS(20)*CONSTS(16))
CONSTS(103) = ( CONSTS(28)*CONSTS(26))/( CONSTS(25)*CONSTS(27))
CONSTS(104) = ( CONSTS(32)*CONSTS(29)*CONSTS(31))/CONSTS(30)
CONSTS(105) = ( CONSTS(43)*CONSTS(39)*CONSTS(33))/( CONSTS(36)*CONSTS(35))
CONSTS(106) = ( CONSTS(35)*CONSTS(36))/( CONSTS(37)*CONSTS(38))
CONSTS(107) = ( CONSTS(51)*CONSTS(49)*CONSTS(47))/( CONSTS(48)*CONSTS(50))
CONSTS(108) = ( CONSTS(55)*CONSTS(53))/( CONSTS(52)*CONSTS(54))
CONSTS(109) = ( CONSTS(62)*1.00000*CONSTS(60)*CONSTS(59))/( CONSTS(56)*CONSTS(58)*CONSTS(57)*CONSTS(61))
CONSTS(110) = ( CONSTS(68)*CONSTS(66)*CONSTS(64)*CONSTS(67))/( CONSTS(63)*CONSTS(65))
CONSTS(111) = ( CONSTS(72)*CONSTS(69))/( CONSTS(70)*CONSTS(71))
CONSTS(112) = ( CONSTS(76)*CONSTS(73))/( CONSTS(74)*CONSTS(75))
CONSTS(113) = ( CONSTS(82)*CONSTS(79)*CONSTS(77))/( CONSTS(80)*CONSTS(78)*CONSTS(81))
CONSTS(114) = ( CONSTS(88)*CONSTS(83)*CONSTS(85))/( CONSTS(86)*CONSTS(84)*CONSTS(87))
CONSTS(115) = ( CONSTS(95)*CONSTS(93)*CONSTS(90)*CONSTS(94))/( CONSTS(89)*CONSTS(91))
CONSTS(116) = ( CONSTS(100)*CONSTS(96) ** 2.00000)/( CONSTS(98)*CONSTS(97)*CONSTS(99))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = ( CONSTS(12)*(( STATES(1)*STATES(14))/( CONSTS(4)*CONSTS(10))) - CONSTS(101)*(( STATES(1)*STATES(2))/( CONSTS(5)*CONSTS(11))))/(1.00000+STATES(1)/CONSTS(4)+STATES(14)/CONSTS(9)+STATES(1)/CONSTS(5)+STATES(2)/CONSTS(11)+( STATES(1)*STATES(14))/( CONSTS(6)*CONSTS(9))+( STATES(1)*STATES(2))/( CONSTS(7)*CONSTS(11)))
ALGBRC(3) = (( CONSTS(24)*(( STATES(1)*STATES(14))/( CONSTS(13)*CONSTS(21))) - CONSTS(102)*(( STATES(1)*STATES(2))/( CONSTS(14)*CONSTS(22))))/(1.00000+STATES(1)/CONSTS(13)+STATES(14)/CONSTS(20)+STATES(1)/CONSTS(14)+STATES(2)/CONSTS(22)+( STATES(1)*STATES(14))/( CONSTS(15)*CONSTS(20))+( STATES(1)*STATES(2))/( CONSTS(19)*CONSTS(22))))*((STATES(17) ** CONSTS(17)/CONSTS(18))/(1.00000+STATES(17) ** CONSTS(17)/CONSTS(18)))
ALGBRC(4) = CONSTS(2)*ALGBRC(2)+ CONSTS(3)*ALGBRC(3)
RATES(1) = - ALGBRC(4)
ALGBRC(5) = ( CONSTS(28)*(STATES(2)/CONSTS(25)) - CONSTS(103)*(STATES(3)/CONSTS(26)))/(1.00000+STATES(2)/CONSTS(25)+STATES(3)/CONSTS(26))
RATES(2) = ALGBRC(4) - ALGBRC(5)
ALGBRC(6) = ( CONSTS(104)*(STATES(3)/CONSTS(29)) - CONSTS(32)*(STATES(4)/CONSTS(30)))/(1.00000+STATES(3)/CONSTS(29)+STATES(4)/CONSTS(30))
RATES(3) = ALGBRC(5) - ALGBRC(6)
ALGBRC(9) = CONSTS(44)* ((1.00000+STATES(16)/CONSTS(41))/(1.00000+ CONSTS(45)*(STATES(16)/CONSTS(41))))*((1.00000+ CONSTS(46)*(STATES(17)/CONSTS(42)))/(1.00000+STATES(17)/CONSTS(42))) ** 4.00000
ALGBRC(7) = (1.00000+STATES(4)/CONSTS(35))*(1.00000+STATES(16)/CONSTS(36))+STATES(15)/CONSTS(39)+ (STATES(5)/CONSTS(33))*(1.00000+STATES(15)/CONSTS(39))
ALGBRC(8) = (1.00000+STATES(4)/CONSTS(37))*(1.00000+STATES(16)/CONSTS(38))+STATES(15)/CONSTS(40)+ (STATES(5)/CONSTS(34))*(1.00000+STATES(15)/CONSTS(40))
ALGBRC(10) = (( CONSTS(43)*(( STATES(16)*STATES(4))/( CONSTS(36)*CONSTS(35))) - CONSTS(105)*(( STATES(15)*STATES(5))/( CONSTS(39)*CONSTS(33))))/ALGBRC(7))*((1.00000+ CONSTS(106)*ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 3.00000)/(1.00000+ ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 4.00000))
RATES(4) = ALGBRC(6) - ALGBRC(10)
ALGBRC(11) = ( CONSTS(51)*(STATES(5)/CONSTS(48)) - CONSTS(107)*(( STATES(6)*STATES(7))/( CONSTS(49)*CONSTS(47))))/(1.00000+STATES(5)/CONSTS(48)+STATES(6)/CONSTS(49)+STATES(7)/CONSTS(47))
RATES(5) = ALGBRC(10) - ALGBRC(11)
ALGBRC(12) = ( CONSTS(55)*(STATES(7)/CONSTS(52)) - CONSTS(108)*(STATES(6)/CONSTS(53)))/(1.00000+STATES(6)/CONSTS(53)+STATES(7)/CONSTS(52))
RATES(6) = ALGBRC(11)+ALGBRC(12)
ALGBRC(13) = 1.00000+STATES(7)/CONSTS(56)+STATES(21)/CONSTS(58)+STATES(14)/CONSTS(57)+( STATES(7)*STATES(21))/( CONSTS(56)*CONSTS(58))+( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))+STATES(8)/CONSTS(60)+STATES(20)/CONSTS(59)+( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))
ALGBRC(14) = ( CONSTS(62)*(( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))) - CONSTS(109)*(( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))))/ALGBRC(13)
RATES(7) = ALGBRC(11) - (ALGBRC(12)+ALGBRC(14))
RATES(14) = CONSTS(1) - (ALGBRC(4)+ALGBRC(14))
ALGBRC(15) = ( CONSTS(110)*(( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))) - CONSTS(68)*(( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65))))/(1.00000+STATES(8)/CONSTS(66)+STATES(15)/CONSTS(64)+( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))+STATES(9)/CONSTS(63)+STATES(16)/CONSTS(65)+( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65)))
RATES(8) = ALGBRC(14) - ALGBRC(15)
ALGBRC(16) = ( CONSTS(72)*(STATES(9)/CONSTS(70)) - CONSTS(111)*(STATES(10)/CONSTS(69)))/(1.00000+STATES(9)/CONSTS(70)+STATES(10)/CONSTS(69))
RATES(9) = ALGBRC(15) - ALGBRC(16)
ALGBRC(17) = ( CONSTS(76)*(STATES(10)/CONSTS(74)) - CONSTS(112)*(STATES(11)/CONSTS(73)))/(1.00000+STATES(10)/CONSTS(74)+STATES(11)/CONSTS(73))
RATES(10) = ALGBRC(16) - ALGBRC(17)
ALGBRC(18) = ( CONSTS(82)*(( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))) - CONSTS(113)*(( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79))))/(1.00000+STATES(11)/CONSTS(80)+STATES(15)/CONSTS(78)+( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))+STATES(12)/CONSTS(77)+STATES(16)/CONSTS(79)+( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79)))
RATES(11) = ALGBRC(17) - ALGBRC(18)
ALGBRC(19) = ( CONSTS(88)*(( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))) - CONSTS(114)*(( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85))))/(1.00000+STATES(12)/CONSTS(86)+STATES(20)/CONSTS(84)+( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))+STATES(13)/CONSTS(83)+STATES(21)/CONSTS(85)+( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85)))
RATES(12) = ALGBRC(18) - ALGBRC(19)
ALGBRC(1) = 0.200000*STATES(13)
RATES(13) = ALGBRC(19) - ALGBRC(1)
ALGBRC(20) = ( CONSTS(115)*(( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90))) - CONSTS(95)*(( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))))/(1.00000+STATES(15)/CONSTS(89)+STATES(18)/CONSTS(92)+( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))+STATES(16)/CONSTS(93)+( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90)))
RATES(18) = ALGBRC(20)
RATES(19) = - ALGBRC(20)
RATES(20) = ALGBRC(14) - ALGBRC(19)
RATES(21) = ALGBRC(19) - ALGBRC(14)
ALGBRC(21) = ( CONSTS(100)*(( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))) - CONSTS(116)*(STATES(15) ** 2.00000/CONSTS(96) ** 2.00000))/(1.00000+STATES(16)/CONSTS(98)+STATES(17)/CONSTS(97)+( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))+( 2.00000*STATES(15))/CONSTS(96)+STATES(15) ** 2.00000/CONSTS(96) ** 2.00000)
RATES(15) = ALGBRC(10)+ 2.00000*ALGBRC(21)+ALGBRC(20)+CONSTS(1)+- (ALGBRC(15)+ALGBRC(18))
RATES(16) = ALGBRC(15)+ALGBRC(18)+- (ALGBRC(10)+ALGBRC(21)+ALGBRC(20)+CONSTS(1))
RATES(17) = - ALGBRC(21)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = ( CONSTS(12)*(( STATES(1)*STATES(14))/( CONSTS(4)*CONSTS(10))) - CONSTS(101)*(( STATES(1)*STATES(2))/( CONSTS(5)*CONSTS(11))))/(1.00000+STATES(1)/CONSTS(4)+STATES(14)/CONSTS(9)+STATES(1)/CONSTS(5)+STATES(2)/CONSTS(11)+( STATES(1)*STATES(14))/( CONSTS(6)*CONSTS(9))+( STATES(1)*STATES(2))/( CONSTS(7)*CONSTS(11)))
ALGBRC(3) = (( CONSTS(24)*(( STATES(1)*STATES(14))/( CONSTS(13)*CONSTS(21))) - CONSTS(102)*(( STATES(1)*STATES(2))/( CONSTS(14)*CONSTS(22))))/(1.00000+STATES(1)/CONSTS(13)+STATES(14)/CONSTS(20)+STATES(1)/CONSTS(14)+STATES(2)/CONSTS(22)+( STATES(1)*STATES(14))/( CONSTS(15)*CONSTS(20))+( STATES(1)*STATES(2))/( CONSTS(19)*CONSTS(22))))*((STATES(17) ** CONSTS(17)/CONSTS(18))/(1.00000+STATES(17) ** CONSTS(17)/CONSTS(18)))
ALGBRC(4) = CONSTS(2)*ALGBRC(2)+ CONSTS(3)*ALGBRC(3)
ALGBRC(5) = ( CONSTS(28)*(STATES(2)/CONSTS(25)) - CONSTS(103)*(STATES(3)/CONSTS(26)))/(1.00000+STATES(2)/CONSTS(25)+STATES(3)/CONSTS(26))
ALGBRC(6) = ( CONSTS(104)*(STATES(3)/CONSTS(29)) - CONSTS(32)*(STATES(4)/CONSTS(30)))/(1.00000+STATES(3)/CONSTS(29)+STATES(4)/CONSTS(30))
ALGBRC(9) = CONSTS(44)* ((1.00000+STATES(16)/CONSTS(41))/(1.00000+ CONSTS(45)*(STATES(16)/CONSTS(41))))*((1.00000+ CONSTS(46)*(STATES(17)/CONSTS(42)))/(1.00000+STATES(17)/CONSTS(42))) ** 4.00000
ALGBRC(7) = (1.00000+STATES(4)/CONSTS(35))*(1.00000+STATES(16)/CONSTS(36))+STATES(15)/CONSTS(39)+ (STATES(5)/CONSTS(33))*(1.00000+STATES(15)/CONSTS(39))
ALGBRC(8) = (1.00000+STATES(4)/CONSTS(37))*(1.00000+STATES(16)/CONSTS(38))+STATES(15)/CONSTS(40)+ (STATES(5)/CONSTS(34))*(1.00000+STATES(15)/CONSTS(40))
ALGBRC(10) = (( CONSTS(43)*(( STATES(16)*STATES(4))/( CONSTS(36)*CONSTS(35))) - CONSTS(105)*(( STATES(15)*STATES(5))/( CONSTS(39)*CONSTS(33))))/ALGBRC(7))*((1.00000+ CONSTS(106)*ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 3.00000)/(1.00000+ ALGBRC(9)*ALGBRC(8)/ALGBRC(7) ** 4.00000))
ALGBRC(11) = ( CONSTS(51)*(STATES(5)/CONSTS(48)) - CONSTS(107)*(( STATES(6)*STATES(7))/( CONSTS(49)*CONSTS(47))))/(1.00000+STATES(5)/CONSTS(48)+STATES(6)/CONSTS(49)+STATES(7)/CONSTS(47))
ALGBRC(12) = ( CONSTS(55)*(STATES(7)/CONSTS(52)) - CONSTS(108)*(STATES(6)/CONSTS(53)))/(1.00000+STATES(6)/CONSTS(53)+STATES(7)/CONSTS(52))
ALGBRC(13) = 1.00000+STATES(7)/CONSTS(56)+STATES(21)/CONSTS(58)+STATES(14)/CONSTS(57)+( STATES(7)*STATES(21))/( CONSTS(56)*CONSTS(58))+( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))+STATES(8)/CONSTS(60)+STATES(20)/CONSTS(59)+( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))
ALGBRC(14) = ( CONSTS(62)*(( STATES(7)*STATES(21)*STATES(14))/( CONSTS(56)*CONSTS(58)*CONSTS(57))) - CONSTS(109)*(( STATES(8)*STATES(20))/( CONSTS(60)*CONSTS(59))))/ALGBRC(13)
ALGBRC(15) = ( CONSTS(110)*(( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))) - CONSTS(68)*(( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65))))/(1.00000+STATES(8)/CONSTS(66)+STATES(15)/CONSTS(64)+( STATES(8)*STATES(15))/( CONSTS(66)*CONSTS(64))+STATES(9)/CONSTS(63)+STATES(16)/CONSTS(65)+( STATES(9)*STATES(16))/( CONSTS(63)*CONSTS(65)))
ALGBRC(16) = ( CONSTS(72)*(STATES(9)/CONSTS(70)) - CONSTS(111)*(STATES(10)/CONSTS(69)))/(1.00000+STATES(9)/CONSTS(70)+STATES(10)/CONSTS(69))
ALGBRC(17) = ( CONSTS(76)*(STATES(10)/CONSTS(74)) - CONSTS(112)*(STATES(11)/CONSTS(73)))/(1.00000+STATES(10)/CONSTS(74)+STATES(11)/CONSTS(73))
ALGBRC(18) = ( CONSTS(82)*(( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))) - CONSTS(113)*(( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79))))/(1.00000+STATES(11)/CONSTS(80)+STATES(15)/CONSTS(78)+( STATES(11)*STATES(15))/( CONSTS(80)*CONSTS(78))+STATES(12)/CONSTS(77)+STATES(16)/CONSTS(79)+( STATES(12)*STATES(16))/( CONSTS(77)*CONSTS(79)))
ALGBRC(19) = ( CONSTS(88)*(( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))) - CONSTS(114)*(( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85))))/(1.00000+STATES(12)/CONSTS(86)+STATES(20)/CONSTS(84)+( STATES(12)*STATES(20))/( CONSTS(86)*CONSTS(84))+STATES(13)/CONSTS(83)+STATES(21)/CONSTS(85)+( STATES(13)*STATES(21))/( CONSTS(83)*CONSTS(85)))
ALGBRC(1) = 0.200000*STATES(13)
ALGBRC(20) = ( CONSTS(115)*(( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90))) - CONSTS(95)*(( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))))/(1.00000+STATES(15)/CONSTS(89)+STATES(18)/CONSTS(92)+( STATES(15)*STATES(18))/( CONSTS(89)*CONSTS(91))+STATES(16)/CONSTS(93)+( STATES(16)*STATES(19))/( CONSTS(93)*CONSTS(90)))
ALGBRC(21) = ( CONSTS(100)*(( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))) - CONSTS(116)*(STATES(15) ** 2.00000/CONSTS(96) ** 2.00000))/(1.00000+STATES(16)/CONSTS(98)+STATES(17)/CONSTS(97)+( STATES(16)*STATES(17))/( CONSTS(98)*CONSTS(97))+( 2.00000*STATES(15))/CONSTS(96)+STATES(15) ** 2.00000/CONSTS(96) ** 2.00000)
RETURN
END