C
C There are a total of 24 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 60 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is Glc_x_0 in component Glc_x_0 (millimolar).
C STATES(1) is Glc_x in component Glc_x (millimolar).
C ALGBRC(2) is GlcTrans in component GlcTrans (flux).
C ALGBRC(1) is inGlc in component inGlc (flux).
C STATES(2) is Glc in component Glc (millimolar).
C ALGBRC(3) is HK in component HK (flux).
C STATES(3) is G6P in component G6P (millimolar).
C ALGBRC(4) is PGI in component PGI (flux).
C ALGBRC(11) is storage in component storage (flux).
C STATES(4) is F6P in component F6P (millimolar).
C ALGBRC(5) is PFK in component PFK (flux).
C STATES(5) is FBP in component FBP (millimolar).
C ALGBRC(6) is ALD in component ALD (flux).
C STATES(6) is GAP in component GAP (millimolar).
C ALGBRC(8) is GAPDH in component GAPDH (flux).
C ALGBRC(7) is TIM in component TIM (flux).
C STATES(7) is DHAP in component DHAP (millimolar).
C ALGBRC(15) is lpGlyc in component lpGlyc (flux).
C STATES(8) is BPG in component BPG (millimolar).
C ALGBRC(9) is lpPEP in component lpPEP (flux).
C STATES(9) is PEP in component PEP (millimolar).
C ALGBRC(10) is PK in component PK (flux).
C STATES(10) is Pyr in component Pyr (millimolar).
C ALGBRC(12) is PDC in component PDC (flux).
C STATES(11) is ATP in component ATP (millimolar).
C ALGBRC(13) is consum in component consum (flux).
C ALGBRC(16) is AK in component AK (flux).
C STATES(12) is ADP in component ADP (millimolar).
C STATES(13) is AMP in component AMP (millimolar).
C CONSTS(2) is CN_x_0 in component CN_x_0 (millimolar).
C STATES(14) is CN_x in component CN_x (millimolar).
C ALGBRC(22) is lacto in component lacto (flux).
C ALGBRC(24) is inCN in component inCN (flux).
C STATES(15) is ACA in component ACA (millimolar).
C ALGBRC(19) is difACA in component difACA (flux).
C ALGBRC(14) is ADH in component ADH (flux).
C STATES(16) is ACA_x in component ACA_x (millimolar).
C ALGBRC(17) is outACA in component outACA (flux).
C STATES(17) is EtOH in component EtOH (millimolar).
C ALGBRC(20) is difEtOH in component difEtOH (flux).
C STATES(18) is EtOH_x in component EtOH_x (millimolar).
C ALGBRC(18) is outEtOH in component outEtOH (flux).
C STATES(19) is Glyc in component Glyc (millimolar).
C ALGBRC(23) is difGlyc in component difGlyc (flux).
C STATES(20) is Glyc_x in component Glyc_x (millimolar).
C ALGBRC(21) is outGlyc in component outGlyc (flux).
C STATES(21) is NADH in component NADH (millimolar).
C STATES(22) is NAD in component NAD (millimolar).
C CONSTS(3) is k0 in component model_parameters (first_order_rate_constant).
C CONSTS(4) is K2_Glc in component GlcTrans (millimolar).
C CONSTS(5) is K2_IG6P in component GlcTrans (millimolar).
C CONSTS(6) is K2_IIG6P in component GlcTrans (millimolar).
C CONSTS(7) is V_2m in component GlcTrans (flux).
C CONSTS(8) is P2 in component GlcTrans (dimensionless).
C CONSTS(9) is y_vol in component model_parameters (dimensionless).
C CONSTS(10) is K3_DGlc in component HK (millimolar).
C CONSTS(11) is K3_Glc in component HK (millimolar).
C CONSTS(12) is K3_ATP in component HK (millimolar).
C CONSTS(13) is V_3m in component HK (flux).
C CONSTS(14) is V_4m in component PGI (flux).
C CONSTS(15) is K4_G6P in component PGI (millimolar).
C CONSTS(16) is K4_F6P in component PGI (millimolar).
C CONSTS(17) is K4_eq in component PGI (dimensionless).
C CONSTS(18) is K5 in component PFK (millimolar2).
C CONSTS(19) is kappa5 in component PFK (dimensionless).
C CONSTS(20) is V_5m in component PFK (flux).
C CONSTS(21) is K6_eq in component ALD (millimolar).
C CONSTS(22) is K6_FBP in component ALD (millimolar).
C CONSTS(23) is K6_DHAP in component ALD (millimolar).
C CONSTS(24) is K6_GAP in component ALD (millimolar).
C CONSTS(25) is K6_IGAP in component ALD (millimolar).
C CONSTS(26) is V_6r in component ALD (flux).
C CONSTS(27) is V_6f in component ALD (flux).
C CONSTS(28) is K7_eq in component TIM (dimensionless).
C CONSTS(29) is K7_DHAP in component TIM (millimolar).
C CONSTS(30) is K7_GAP in component TIM (millimolar).
C CONSTS(31) is V_7m in component TIM (flux).
C CONSTS(32) is K8_NAD in component GAPDH (millimolar).
C CONSTS(33) is K8_NADH in component GAPDH (millimolar).
C CONSTS(34) is K8_GAP in component GAPDH (millimolar).
C CONSTS(35) is K8_BPG in component GAPDH (millimolar).
C CONSTS(36) is K8_eq in component GAPDH (dimensionless).
C CONSTS(37) is V_8m in component GAPDH (flux).
C CONSTS(38) is k9f in component lpPEP (second_order_rate_constant).
C CONSTS(39) is k9r in component lpPEP (second_order_rate_constant).
C CONSTS(40) is K10_PEP in component PK (millimolar).
C CONSTS(41) is K10_ADP in component PK (millimolar).
C CONSTS(42) is V_10m in component PK (flux).
C CONSTS(43) is K11 in component PDC (millimolar).
C CONSTS(44) is V_11m in component PDC (flux).
C CONSTS(45) is K12_NADH in component ADH (millimolar).
C CONSTS(46) is K12_ACA in component ADH (millimolar).
C CONSTS(47) is V_12m in component ADH (flux).
C CONSTS(48) is k13 in component difEtOH (first_order_rate_constant).
C CONSTS(49) is K15_NADH in component lpGlyc (millimolar).
C CONSTS(50) is K15_INADH in component lpGlyc (millimolar).
C CONSTS(51) is K15_INAD in component lpGlyc (millimolar).
C CONSTS(52) is K15_DHAP in component lpGlyc (millimolar).
C CONSTS(53) is V_15m in component lpGlyc (flux).
C CONSTS(54) is k16 in component difGlyc (first_order_rate_constant).
C CONSTS(55) is k18 in component difACA (first_order_rate_constant).
C CONSTS(56) is k20 in component lacto (second_order_rate_constant).
C CONSTS(57) is k22 in component storage (second_order_rate_constant).
C CONSTS(58) is k23 in component consum (first_order_rate_constant).
C CONSTS(59) is k24f in component AK (second_order_rate_constant).
C CONSTS(60) is k24r in component AK (second_order_rate_constant).
C RATES(1) is d/dt Glc_x in component Glc_x (millimolar).
C RATES(2) is d/dt Glc in component Glc (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 GAP in component GAP (millimolar).
C RATES(7) is d/dt DHAP in component DHAP (millimolar).
C RATES(8) is d/dt BPG in component BPG (millimolar).
C RATES(9) is d/dt PEP in component PEP (millimolar).
C RATES(10) is d/dt Pyr in component Pyr (millimolar).
C RATES(11) is d/dt ATP in component ATP (millimolar).
C RATES(12) is d/dt ADP in component ADP (millimolar).
C RATES(13) is d/dt AMP in component AMP (millimolar).
C RATES(14) is d/dt CN_x in component CN_x (millimolar).
C RATES(15) is d/dt ACA in component ACA (millimolar).
C RATES(16) is d/dt ACA_x in component ACA_x (millimolar).
C RATES(17) is d/dt EtOH in component EtOH (millimolar).
C RATES(18) is d/dt EtOH_x in component EtOH_x (millimolar).
C RATES(19) is d/dt Glyc in component Glyc (millimolar).
C RATES(20) is d/dt Glyc_x in component Glyc_x (millimolar).
C RATES(21) is d/dt NADH in component NADH (millimolar).
C RATES(22) is d/dt NAD in component NAD (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 24.0
      STATES(1) = 6.7
      STATES(2) = 0.573074
      STATES(3) = 4.2
      STATES(4) = 0.49
      STATES(5) = 4.64
      STATES(6) = 0.115
      STATES(7) = 2.95
      STATES(8) = 0.00027
      STATES(9) = 0.04
      STATES(10) = 8.7
      STATES(11) = 2.1
      STATES(12) = 1.5
      STATES(13) = 0.33
      CONSTS(2) = 5.60
      STATES(14) = 5.20358
      STATES(15) = 1.48153
      STATES(16) = 1.28836
      STATES(17) = 19.2379
      STATES(18) = 16.4514
      STATES(19) = 4.196
      STATES(20) = 1.68478
      STATES(21) = 0.33
      STATES(22) = 0.65
      CONSTS(3) = 0.048
      CONSTS(4) = 1.7
      CONSTS(5) = 1.2
      CONSTS(6) = 7.2
      CONSTS(7) = 1014.96
      CONSTS(8) = 1.0
      CONSTS(9) = 59.0
      CONSTS(10) = 0.37
      CONSTS(11) = 0.0
      CONSTS(12) = 0.1
      CONSTS(13) = 51.7547
      CONSTS(14) = 496.042
      CONSTS(15) = 0.8
      CONSTS(16) = 0.15
      CONSTS(17) = 0.13
      CONSTS(18) = 0.021
      CONSTS(19) = 0.15
      CONSTS(20) = 45.4327
      CONSTS(21) = 0.081
      CONSTS(22) = 0.3
      CONSTS(23) = 2.0
      CONSTS(24) = 4.0
      CONSTS(25) = 10.0
      CONSTS(26) = 1.10391E4
      CONSTS(27) = 2.20782E3
      CONSTS(28) = 0.055
      CONSTS(29) = 1.23
      CONSTS(30) = 1.27
      CONSTS(31) = 1.16365E2
      CONSTS(32) = 0.1
      CONSTS(33) = 0.06
      CONSTS(34) = 0.6
      CONSTS(35) = 0.01
      CONSTS(36) = 0.0055
      CONSTS(37) = 8.33858E2
      CONSTS(38) = 4.43866E5
      CONSTS(39) = 1.52862E3
      CONSTS(40) = 0.2
      CONSTS(41) = 0.17
      CONSTS(42) = 3.43096E2
      CONSTS(43) = 0.3
      CONSTS(44) = 5.31328E1
      CONSTS(45) = 0.1
      CONSTS(46) = 0.71
      CONSTS(47) = 8.98023E1
      CONSTS(48) = 16.72
      CONSTS(49) = 0.13
      CONSTS(50) = 0.034
      CONSTS(51) = 0.13
      CONSTS(52) = 25.0
      CONSTS(53) = 8.14797E1
      CONSTS(54) = 1.9
      CONSTS(55) = 24.7
      CONSTS(56) = 2.83828E-3
      CONSTS(57) = 2.25932
      CONSTS(58) = 3.20760
      CONSTS(59) = 4.32900E2
      CONSTS(60) = 1.33333E2
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) =  (CONSTS(7)/CONSTS(9))*((STATES(1)/CONSTS(4))/(1.00000+STATES(1)/CONSTS(4)+ (( CONSTS(8)*(STATES(1)/CONSTS(4))+1.00000)/( CONSTS(8)*(STATES(2)/CONSTS(4))+1.00000))*(1.00000+STATES(2)/CONSTS(4)+STATES(3)/CONSTS(5)+( STATES(2)*STATES(3))/( CONSTS(4)*CONSTS(6))))) -  (CONSTS(7)/CONSTS(9))*((STATES(2)/CONSTS(4))/(1.00000+STATES(2)/CONSTS(4)+ (( CONSTS(8)*(STATES(2)/CONSTS(4))+1.00000)/( CONSTS(8)*(STATES(1)/CONSTS(4))+1.00000))*(1.00000+STATES(1)/CONSTS(4))+STATES(3)/CONSTS(5)+( STATES(2)*STATES(3))/( CONSTS(4)*CONSTS(6))))
      ALGBRC(1) =  CONSTS(3)*(CONSTS(1) - STATES(1))
      RATES(1) = ALGBRC(1) - ALGBRC(2)
      ALGBRC(3) = ( CONSTS(13)*STATES(11)*STATES(2))/( CONSTS(10)*CONSTS(12)+ CONSTS(11)*STATES(11)+ CONSTS(12)*STATES(2)+ STATES(2)*STATES(11))
      RATES(2) =  59.0000*ALGBRC(2) - ALGBRC(3)
      ALGBRC(4) = ( CONSTS(14)*STATES(3))/(CONSTS(15)+STATES(3)+ (CONSTS(15)/CONSTS(16))*STATES(4)) - ( CONSTS(14)*(STATES(4)/CONSTS(17)))/(CONSTS(15)+STATES(3)+ (CONSTS(15)/CONSTS(16))*STATES(4))
      ALGBRC(5) = ( CONSTS(20)*STATES(4) ** 2.00000)/( CONSTS(18)*(1.00000+ CONSTS(19)*(STATES(11)/STATES(13))*(STATES(11)/STATES(13)))+STATES(4) ** 2.00000)
      RATES(4) = ALGBRC(4) - ALGBRC(5)
      ALGBRC(6) = ( CONSTS(27)*STATES(5))/(CONSTS(22)+STATES(5)+( STATES(6)*CONSTS(23)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(7)*CONSTS(24)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(5)*STATES(6))/CONSTS(25)+( STATES(6)*STATES(7)*CONSTS(27))/( CONSTS(21)*CONSTS(26))) - ( CONSTS(27)*(( STATES(6)*STATES(7))/CONSTS(21)))/(CONSTS(22)+STATES(5)+( STATES(6)*CONSTS(23)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(7)*CONSTS(24)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(5)*STATES(6))/CONSTS(25)+( STATES(6)*STATES(7)*CONSTS(27))/( CONSTS(21)*CONSTS(26)))
      RATES(5) = ALGBRC(5) - ALGBRC(6)
      ALGBRC(8) = ( CONSTS(37)*STATES(6)*STATES(22))/( CONSTS(34)*CONSTS(32)*(1.00000+STATES(6)/CONSTS(34)+STATES(8)/CONSTS(35))*(1.00000+STATES(22)/CONSTS(32)+STATES(21)/CONSTS(33))) - ( CONSTS(37)*(( STATES(8)*STATES(21))/CONSTS(36)))/( CONSTS(34)*CONSTS(32)*(1.00000+STATES(6)/CONSTS(34)+STATES(8)/CONSTS(35))*(1.00000+STATES(22)/CONSTS(32)+STATES(21)/CONSTS(33)))
      ALGBRC(7) = ( CONSTS(31)*STATES(7))/(CONSTS(29)+STATES(7)+ (CONSTS(29)/CONSTS(30))*STATES(6)) - ( CONSTS(31)*(STATES(6)/CONSTS(28)))/(CONSTS(29)+STATES(7)+ (CONSTS(29)/CONSTS(30))*STATES(6))
      RATES(6) = (ALGBRC(6)+ALGBRC(7)) - ALGBRC(8)
      ALGBRC(9) =  CONSTS(38)*STATES(8)*STATES(12) -  CONSTS(39)*STATES(9)*STATES(11)
      RATES(8) = ALGBRC(8) - ALGBRC(9)
      ALGBRC(10) = ( CONSTS(42)*STATES(12)*STATES(9))/( (CONSTS(40)+STATES(9))*(CONSTS(41)+STATES(12)))
      RATES(9) = ALGBRC(9) - ALGBRC(10)
      ALGBRC(11) =  CONSTS(57)*STATES(11)*STATES(3)
      RATES(3) = ALGBRC(3) - (ALGBRC(4)+ALGBRC(11))
      ALGBRC(12) = ( CONSTS(44)*STATES(10))/(CONSTS(43)+STATES(10))
      RATES(10) = ALGBRC(10) - ALGBRC(12)
      ALGBRC(15) = ( CONSTS(53)*STATES(7))/( CONSTS(52)*(1.00000+ (CONSTS(50)/STATES(21))*(1.00000+STATES(22)/CONSTS(51)))+ STATES(7)*(1.00000+ (CONSTS(49)/STATES(21))*(1.00000+STATES(22)/CONSTS(51))))
      RATES(7) = ALGBRC(6) - (ALGBRC(7)+ALGBRC(15))
      ALGBRC(13) =  CONSTS(58)*STATES(11)
      ALGBRC(16) =  CONSTS(59)*STATES(13)*STATES(11) -  CONSTS(60)*STATES(12) ** 2.00000
      RATES(11) = (ALGBRC(10)+ALGBRC(9)) - (ALGBRC(5)+ALGBRC(11)+ALGBRC(3)+ALGBRC(13)+ALGBRC(16))
      RATES(12) = (ALGBRC(5)+ALGBRC(11)+ALGBRC(3)+ALGBRC(13)+ 2.00000*ALGBRC(16)) - (ALGBRC(10)+ALGBRC(9))
      RATES(13) = - ALGBRC(16)
      ALGBRC(14) = ( CONSTS(47)*STATES(15)*STATES(21))/( (CONSTS(45)+STATES(21))*(CONSTS(46)+STATES(15)))
      RATES(21) = ALGBRC(8) - (ALGBRC(14)+ALGBRC(15))
      RATES(22) = (ALGBRC(14)+ALGBRC(15)) - ALGBRC(8)
      ALGBRC(19) =  (CONSTS(55)/CONSTS(9))*(STATES(15) - STATES(16))
      RATES(15) = ALGBRC(12) - ( 59.0000*ALGBRC(19)+ALGBRC(14))
      ALGBRC(20) =  (CONSTS(48)/CONSTS(9))*(STATES(17) - STATES(18))
      RATES(17) = ALGBRC(14) -  59.0000*ALGBRC(20)
      ALGBRC(18) =  CONSTS(3)*STATES(18)
      RATES(18) = ALGBRC(20) - ALGBRC(18)
      ALGBRC(22) =  CONSTS(56)*STATES(16)*STATES(14)
      ALGBRC(17) =  CONSTS(3)*STATES(16)
      RATES(16) = ALGBRC(19) - (ALGBRC(22)+ALGBRC(17))
      ALGBRC(23) =  (CONSTS(54)/CONSTS(9))*(STATES(19) - STATES(20))
      RATES(19) = ALGBRC(15) -  59.0000*ALGBRC(23)
      ALGBRC(21) =  CONSTS(3)*STATES(20)
      RATES(20) = ALGBRC(23) - ALGBRC(21)
      ALGBRC(24) =  CONSTS(3)*(CONSTS(2) - STATES(14))
      RATES(14) = ALGBRC(24) - ALGBRC(22)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) =  (CONSTS(7)/CONSTS(9))*((STATES(1)/CONSTS(4))/(1.00000+STATES(1)/CONSTS(4)+ (( CONSTS(8)*(STATES(1)/CONSTS(4))+1.00000)/( CONSTS(8)*(STATES(2)/CONSTS(4))+1.00000))*(1.00000+STATES(2)/CONSTS(4)+STATES(3)/CONSTS(5)+( STATES(2)*STATES(3))/( CONSTS(4)*CONSTS(6))))) -  (CONSTS(7)/CONSTS(9))*((STATES(2)/CONSTS(4))/(1.00000+STATES(2)/CONSTS(4)+ (( CONSTS(8)*(STATES(2)/CONSTS(4))+1.00000)/( CONSTS(8)*(STATES(1)/CONSTS(4))+1.00000))*(1.00000+STATES(1)/CONSTS(4))+STATES(3)/CONSTS(5)+( STATES(2)*STATES(3))/( CONSTS(4)*CONSTS(6))))
      ALGBRC(1) =  CONSTS(3)*(CONSTS(1) - STATES(1))
      ALGBRC(3) = ( CONSTS(13)*STATES(11)*STATES(2))/( CONSTS(10)*CONSTS(12)+ CONSTS(11)*STATES(11)+ CONSTS(12)*STATES(2)+ STATES(2)*STATES(11))
      ALGBRC(4) = ( CONSTS(14)*STATES(3))/(CONSTS(15)+STATES(3)+ (CONSTS(15)/CONSTS(16))*STATES(4)) - ( CONSTS(14)*(STATES(4)/CONSTS(17)))/(CONSTS(15)+STATES(3)+ (CONSTS(15)/CONSTS(16))*STATES(4))
      ALGBRC(5) = ( CONSTS(20)*STATES(4) ** 2.00000)/( CONSTS(18)*(1.00000+ CONSTS(19)*(STATES(11)/STATES(13))*(STATES(11)/STATES(13)))+STATES(4) ** 2.00000)
      ALGBRC(6) = ( CONSTS(27)*STATES(5))/(CONSTS(22)+STATES(5)+( STATES(6)*CONSTS(23)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(7)*CONSTS(24)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(5)*STATES(6))/CONSTS(25)+( STATES(6)*STATES(7)*CONSTS(27))/( CONSTS(21)*CONSTS(26))) - ( CONSTS(27)*(( STATES(6)*STATES(7))/CONSTS(21)))/(CONSTS(22)+STATES(5)+( STATES(6)*CONSTS(23)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(7)*CONSTS(24)*CONSTS(27))/( CONSTS(21)*CONSTS(26))+( STATES(5)*STATES(6))/CONSTS(25)+( STATES(6)*STATES(7)*CONSTS(27))/( CONSTS(21)*CONSTS(26)))
      ALGBRC(8) = ( CONSTS(37)*STATES(6)*STATES(22))/( CONSTS(34)*CONSTS(32)*(1.00000+STATES(6)/CONSTS(34)+STATES(8)/CONSTS(35))*(1.00000+STATES(22)/CONSTS(32)+STATES(21)/CONSTS(33))) - ( CONSTS(37)*(( STATES(8)*STATES(21))/CONSTS(36)))/( CONSTS(34)*CONSTS(32)*(1.00000+STATES(6)/CONSTS(34)+STATES(8)/CONSTS(35))*(1.00000+STATES(22)/CONSTS(32)+STATES(21)/CONSTS(33)))
      ALGBRC(7) = ( CONSTS(31)*STATES(7))/(CONSTS(29)+STATES(7)+ (CONSTS(29)/CONSTS(30))*STATES(6)) - ( CONSTS(31)*(STATES(6)/CONSTS(28)))/(CONSTS(29)+STATES(7)+ (CONSTS(29)/CONSTS(30))*STATES(6))
      ALGBRC(9) =  CONSTS(38)*STATES(8)*STATES(12) -  CONSTS(39)*STATES(9)*STATES(11)
      ALGBRC(10) = ( CONSTS(42)*STATES(12)*STATES(9))/( (CONSTS(40)+STATES(9))*(CONSTS(41)+STATES(12)))
      ALGBRC(11) =  CONSTS(57)*STATES(11)*STATES(3)
      ALGBRC(12) = ( CONSTS(44)*STATES(10))/(CONSTS(43)+STATES(10))
      ALGBRC(15) = ( CONSTS(53)*STATES(7))/( CONSTS(52)*(1.00000+ (CONSTS(50)/STATES(21))*(1.00000+STATES(22)/CONSTS(51)))+ STATES(7)*(1.00000+ (CONSTS(49)/STATES(21))*(1.00000+STATES(22)/CONSTS(51))))
      ALGBRC(13) =  CONSTS(58)*STATES(11)
      ALGBRC(16) =  CONSTS(59)*STATES(13)*STATES(11) -  CONSTS(60)*STATES(12) ** 2.00000
      ALGBRC(14) = ( CONSTS(47)*STATES(15)*STATES(21))/( (CONSTS(45)+STATES(21))*(CONSTS(46)+STATES(15)))
      ALGBRC(19) =  (CONSTS(55)/CONSTS(9))*(STATES(15) - STATES(16))
      ALGBRC(20) =  (CONSTS(48)/CONSTS(9))*(STATES(17) - STATES(18))
      ALGBRC(18) =  CONSTS(3)*STATES(18)
      ALGBRC(22) =  CONSTS(56)*STATES(16)*STATES(14)
      ALGBRC(17) =  CONSTS(3)*STATES(16)
      ALGBRC(23) =  (CONSTS(54)/CONSTS(9))*(STATES(19) - STATES(20))
      ALGBRC(21) =  CONSTS(3)*STATES(20)
      ALGBRC(24) =  CONSTS(3)*(CONSTS(2) - STATES(14))
      RETURN
      END