C
C There are a total of 30 entries in the algebraic variable array.
C There are a total of 8 entries in each of the rate and state variable arrays.
C There are a total of 81 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is PYR in component PYR (millimolar).
C ALGBRC(5) is V_GLYC in component V_GLYC (flux).
C ALGBRC(7) is V_LDH in component V_LDH (flux).
C ALGBRC(9) is V_PDH in component V_PDH (flux).
C ALGBRC(21) is V_ALS in component V_ALS (flux).
C STATES(2) is ACP in component ACP (millimolar).
C ALGBRC(11) is V_PTA in component V_PTA (flux).
C ALGBRC(15) is V_ACK in component V_ACK (flux).
C STATES(3) is ACAL in component ACAL (millimolar).
C ALGBRC(14) is V_ACALDH in component V_ACALDH (flux).
C ALGBRC(18) is V_ADH in component V_ADH (flux).
C STATES(4) is ACLAC in component ACLAC (millimolar).
C ALGBRC(23) is V_ALDC in component V_ALDC (flux).
C ALGBRC(28) is V_NEALC in component V_NEALC (flux).
C STATES(5) is ACET in component ACET (millimolar).
C ALGBRC(27) is V_ACETDH in component V_ACETDH (flux).
C ALGBRC(25) is V_ACETEFF in component V_ACETEFF (flux).
C STATES(6) is ATP in component ATP (millimolar).
C ALGBRC(19) is V_ATPase in component V_ATPase (flux).
C ALGBRC(1) is ADP in component ADP (millimolar).
C CONSTS(1) is A_tot in component ADP (millimolar).
C STATES(7) is NADH in component NADH (millimolar).
C ALGBRC(30) is V_NOX in component V_NOX (flux).
C ALGBRC(2) is NAD in component NAD (millimolar).
C CONSTS(2) is NAD_tot in component NAD (millimolar).
C STATES(8) is ACCOA in component ACCOA (millimolar).
C ALGBRC(3) is COA in component COA (millimolar).
C CONSTS(3) is C_tot in component COA (millimolar).
C CONSTS(4) is AC in component AC (millimolar).
C CONSTS(5) is BUT in component BUT (millimolar).
C CONSTS(6) is ETOH in component ETOH (millimolar).
C CONSTS(7) is GLC in component GLC (millimolar).
C CONSTS(8) is LAC in component LAC (millimolar).
C CONSTS(9) is O in component O (millimolar).
C CONSTS(10) is P in component P (millimolar).
C ALGBRC(4) is V_GLYC_temp in component V_GLYC (flux).
C CONSTS(11) is Km_GLC in component V_GLYC (millimolar).
C CONSTS(12) is Km_NAD in component V_GLYC (millimolar).
C CONSTS(13) is Km_ADP in component V_GLYC (millimolar).
C CONSTS(14) is Km_PYR in component V_GLYC (millimolar).
C CONSTS(15) is Km_NADH in component V_GLYC (millimolar).
C CONSTS(16) is Km_ATP in component V_GLYC (millimolar).
C CONSTS(17) is V_GLYC_max in component V_GLYC (flux).
C ALGBRC(6) is V_LDH_temp in component V_LDH (flux).
C CONSTS(18) is Keq in component V_LDH (dimensionless).
C CONSTS(19) is Km_LAC in component V_LDH (millimolar).
C CONSTS(20) is Km_NAD in component V_LDH (millimolar).
C CONSTS(21) is Km_PYR in component V_LDH (millimolar).
C CONSTS(22) is Km_NADH in component V_LDH (millimolar).
C CONSTS(23) is V_LDH_max in component V_LDH (flux).
C ALGBRC(8) is V_PDH_temp in component V_PDH (flux).
C CONSTS(24) is Ki in component V_PDH (dimensionless).
C CONSTS(25) is Km_NAD in component V_PDH (millimolar).
C CONSTS(26) is Km_COA in component V_PDH (millimolar).
C CONSTS(27) is Km_PYR in component V_PDH (millimolar).
C CONSTS(28) is Km_NADH in component V_PDH (millimolar).
C CONSTS(29) is Km_ACCOA in component V_PDH (millimolar).
C CONSTS(30) is V_PDH_max in component V_PDH (flux).
C ALGBRC(10) is V_PTA_temp in component V_PTA (flux).
C CONSTS(31) is Keq in component V_PTA (dimensionless).
C CONSTS(32) is Km_P in component V_PTA (millimolar).
C CONSTS(33) is Ki_P in component V_PTA (millimolar).
C CONSTS(34) is Ki_COA in component V_PTA (millimolar).
C CONSTS(35) is Km_ACP in component V_PTA (millimolar).
C CONSTS(36) is Ki_ACP in component V_PTA (millimolar).
C CONSTS(37) is Ki_ACCOA in component V_PTA (millimolar).
C CONSTS(38) is V_PTA_max in component V_PTA (flux).
C ALGBRC(13) is V_ACK_temp in component V_ACK (flux).
C CONSTS(39) is Keq in component V_ACK (dimensionless).
C CONSTS(40) is Km_AC in component V_ACK (millimolar).
C CONSTS(41) is Km_ATP in component V_ACK (millimolar).
C CONSTS(42) is Km_ADP in component V_ACK (millimolar).
C CONSTS(43) is Km_ACP in component V_ACK (millimolar).
C CONSTS(44) is V_ACK_max in component V_ACK (flux).
C ALGBRC(12) is V_ACALDH_temp in component V_ACALDH (flux).
C CONSTS(45) is Keq in component V_ACALDH (millimolar).
C CONSTS(46) is Km_NAD in component V_ACALDH (millimolar).
C CONSTS(47) is Km_NADH in component V_ACALDH (millimolar).
C CONSTS(48) is Km_COA in component V_ACALDH (millimolar).
C CONSTS(49) is Km_ACCOA in component V_ACALDH (millimolar).
C CONSTS(50) is Km_ACAL in component V_ACALDH (millimolar).
C CONSTS(51) is V_ACALDH_max in component V_ACALDH (flux).
C ALGBRC(16) is V_ADH_temp in component V_ADH (flux).
C CONSTS(52) is Keq in component V_ADH (dimensionless).
C CONSTS(53) is Km_NAD in component V_ADH (millimolar).
C CONSTS(54) is Km_NADH in component V_ADH (millimolar).
C CONSTS(55) is Km_ETOH in component V_ADH (millimolar).
C CONSTS(56) is Km_ACAL in component V_ADH (millimolar).
C CONSTS(57) is V_ADH_max in component V_ADH (flux).
C ALGBRC(20) is V_ALS_temp in component V_ALS (flux).
C CONSTS(58) is N in component V_ALS (dimensionless).
C CONSTS(59) is Keq in component V_ALS (dimensionless).
C CONSTS(60) is Km_ACLAC in component V_ALS (millimolar).
C CONSTS(61) is Km_PYR in component V_ALS (millimolar).
C CONSTS(62) is V_ALS_max in component V_ALS (flux).
C ALGBRC(22) is V_ALDC_temp in component V_ALDC (flux).
C CONSTS(63) is Km_ACLAC in component V_ALDC (millimolar).
C CONSTS(64) is Km_ACET in component V_ALDC (millimolar).
C CONSTS(65) is V_ALDC_max in component V_ALDC (flux).
C ALGBRC(24) is V_ACETEFF_temp in component V_ACETEFF (flux).
C CONSTS(66) is Km_ACET in component V_ACETEFF (millimolar).
C CONSTS(67) is V_ACETEFF_max in component V_ACETEFF (flux).
C ALGBRC(26) is V_ACETDH_temp in component V_ACETDH (flux).
C CONSTS(68) is Keq in component V_ACETDH (dimensionless).
C CONSTS(69) is Km_NAD in component V_ACETDH (millimolar).
C CONSTS(70) is Km_NADH in component V_ACETDH (millimolar).
C CONSTS(71) is Km_BUT in component V_ACETDH (millimolar).
C CONSTS(72) is Km_ACET in component V_ACETDH (millimolar).
C CONSTS(73) is V_ACETDH_max in component V_ACETDH (flux).
C ALGBRC(17) is V_ATPase_temp in component V_ATPase (flux).
C CONSTS(74) is N in component V_ATPase (dimensionless).
C CONSTS(75) is Km_ATP in component V_ATPase (dimensionless).
C CONSTS(76) is V_ATPase_max in component V_ATPase (flux).
C ALGBRC(29) is V_NOX_temp in component V_NOX (flux).
C CONSTS(77) is Km_NAD in component V_NOX (millimolar).
C CONSTS(78) is Km_NADH in component V_NOX (millimolar).
C CONSTS(79) is Km_O in component V_NOX (millimolar).
C CONSTS(80) is V_NOX_max in component V_NOX (flux).
C CONSTS(81) is k in component V_NEALC (first_order_rate_constant).
C RATES(1) is d/dt PYR in component PYR (millimolar).
C RATES(2) is d/dt ACP in component ACP (millimolar).
C RATES(3) is d/dt ACAL in component ACAL (millimolar).
C RATES(4) is d/dt ACLAC in component ACLAC (millimolar).
C RATES(5) is d/dt ACET in component ACET (millimolar).
C RATES(6) is d/dt ATP in component ATP (millimolar).
C RATES(7) is d/dt NADH in component NADH (millimolar).
C RATES(8) is d/dt ACCOA in component ACCOA (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 1
      STATES(2) = 0.03145
      STATES(3) = 0.11
      STATES(4) = 1e-5
      STATES(5) = 1e-5
      STATES(6) = 0.1
      CONSTS(1) = 5
      STATES(7) = 3.67
      CONSTS(2) = 10
      STATES(8) = 0.11
      CONSTS(3) = 1
      CONSTS(4) = 0.01
      CONSTS(5) = 0.01
      CONSTS(6) = 0.1
      CONSTS(7) = 15
      CONSTS(8) = 0.1
      CONSTS(9) = 0.2
      CONSTS(10) = 10
      CONSTS(11) = 0.1
      CONSTS(12) = 0.1412
      CONSTS(13) = 0.04699
      CONSTS(14) = 2.5
      CONSTS(15) = 0.08999
      CONSTS(16) = 0.01867
      CONSTS(17) = 2397
      CONSTS(18) = 21120.69
      CONSTS(19) = 100
      CONSTS(20) = 2.4
      CONSTS(21) = 1.5
      CONSTS(22) = 0.08
      CONSTS(23) = 5118
      CONSTS(24) = 46.4159
      CONSTS(25) = 0.4
      CONSTS(26) = 0.014
      CONSTS(27) = 1
      CONSTS(28) = 0.1
      CONSTS(29) = 0.008
      CONSTS(30) = 259
      CONSTS(31) = 0.0065
      CONSTS(32) = 2.6
      CONSTS(33) = 2.6
      CONSTS(34) = 0.029
      CONSTS(35) = 0.7
      CONSTS(36) = 0.2
      CONSTS(37) = 0.2
      CONSTS(38) = 42
      CONSTS(39) = 174.217
      CONSTS(40) = 7
      CONSTS(41) = 0.07
      CONSTS(42) = 0.5
      CONSTS(43) = 0.16
      CONSTS(44) = 2700
      CONSTS(45) = 1
      CONSTS(46) = 0.08
      CONSTS(47) = 0.025
      CONSTS(48) = 0.008
      CONSTS(49) = 0.007
      CONSTS(50) = 10
      CONSTS(51) = 97
      CONSTS(52) = 12354.9
      CONSTS(53) = 0.08
      CONSTS(54) = 0.05
      CONSTS(55) = 1
      CONSTS(56) = 0.03
      CONSTS(57) = 162
      CONSTS(58) = 2.4
      CONSTS(59) = 9e12
      CONSTS(60) = 100
      CONSTS(61) = 50
      CONSTS(62) = 600
      CONSTS(63) = 10
      CONSTS(64) = 100
      CONSTS(65) = 106
      CONSTS(66) = 5
      CONSTS(67) = 200
      CONSTS(68) = 1400
      CONSTS(69) = 0.16
      CONSTS(70) = 0.02
      CONSTS(71) = 2.6
      CONSTS(72) = 0.06
      CONSTS(73) = 105
      CONSTS(74) = 2.58
      CONSTS(75) = 6.196
      CONSTS(76) = 900
      CONSTS(77) = 1
      CONSTS(78) = 0.041
      CONSTS(79) = 0.2
      CONSTS(80) = 118
      CONSTS(81) = 0.0003
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = CONSTS(3) - STATES(8)
      ALGBRC(10) = ( (CONSTS(38)/( CONSTS(37)*CONSTS(32)))*( STATES(8)*CONSTS(10) - ( STATES(2)*ALGBRC(3))/CONSTS(31)))/(1.00000+STATES(8)/CONSTS(37)+CONSTS(10)/CONSTS(33)+STATES(2)/CONSTS(36)+ALGBRC(3)/CONSTS(34)+( STATES(8)*CONSTS(10))/( CONSTS(37)*CONSTS(32))+( STATES(2)*ALGBRC(3))/( CONSTS(35)*CONSTS(34)))
      ALGBRC(11) = TERNRY(ALGBRC(10).GE.0.00000, ALGBRC(10), ALGBRC(10))
      ALGBRC(1) = CONSTS(1) - STATES(6)
      ALGBRC(13) = ( (CONSTS(44)/( CONSTS(42)*CONSTS(43)))*( STATES(2)*ALGBRC(1) - ( CONSTS(4)*STATES(6))/CONSTS(39)))/( (1.00000+STATES(2)/CONSTS(43)+CONSTS(4)/CONSTS(40))*(1.00000+ALGBRC(1)/CONSTS(42)+STATES(6)/CONSTS(41)))
      ALGBRC(15) = TERNRY(ALGBRC(13).GE.0.00000, ALGBRC(13), ALGBRC(13))
      RATES(2) = ALGBRC(11) - ALGBRC(15)
      ALGBRC(2) = CONSTS(2) - STATES(7)
      ALGBRC(8) = (( (( (( (CONSTS(30)/(1.00000+( CONSTS(24)*STATES(7))/ALGBRC(2)))*STATES(1))/CONSTS(27))*ALGBRC(2))/CONSTS(25))*ALGBRC(3))/CONSTS(26))/( (1.00000+STATES(1)/CONSTS(27))*(1.00000+ALGBRC(2)/CONSTS(25)+STATES(7)/CONSTS(28))*(1.00000+ALGBRC(3)/CONSTS(26)+STATES(8)/CONSTS(29)))
      ALGBRC(9) = TERNRY(ALGBRC(8).GE.0.00000, ALGBRC(8), ALGBRC(8))
      ALGBRC(12) = ( (CONSTS(51)/( CONSTS(49)*CONSTS(47)))*( STATES(8)*STATES(7) - ( ALGBRC(2)*ALGBRC(3)*STATES(3))/CONSTS(45)))/( (1.00000+ALGBRC(2)/CONSTS(46)+STATES(7)/CONSTS(47))*(1.00000+STATES(8)/CONSTS(49)+ALGBRC(3)/CONSTS(48))*(1.00000+STATES(3)/CONSTS(50)))
      ALGBRC(14) = TERNRY(ALGBRC(12).GE.0.00000, ALGBRC(12), ALGBRC(12))
      RATES(8) = ALGBRC(9) - (ALGBRC(14)+ALGBRC(11))
      ALGBRC(16) = ( (CONSTS(57)/( CONSTS(56)*CONSTS(54)))*( STATES(3)*STATES(7) - ( CONSTS(6)*ALGBRC(2))/CONSTS(52)))/( (1.00000+ALGBRC(2)/CONSTS(53)+STATES(7)/CONSTS(54))*(1.00000+STATES(3)/CONSTS(56)+CONSTS(6)/CONSTS(55)))
      ALGBRC(18) = TERNRY(ALGBRC(16).GE.0.00000, ALGBRC(16), ALGBRC(16))
      RATES(3) = ALGBRC(14) - ALGBRC(18)
      ALGBRC(4) = (( (( (( CONSTS(17)*CONSTS(7))/CONSTS(11))*ALGBRC(2))/CONSTS(12))*ALGBRC(1))/CONSTS(13))/( (1.00000+CONSTS(7)/CONSTS(11)+STATES(1)/CONSTS(14))*(1.00000+ALGBRC(2)/CONSTS(12)+STATES(7)/CONSTS(15))*(1.00000+ALGBRC(1)/CONSTS(13)+STATES(6)/CONSTS(16)))
      ALGBRC(5) = TERNRY(ALGBRC(4).GE.0.00000, ALGBRC(4), ALGBRC(4))
      ALGBRC(17) = ( CONSTS(76)*STATES(6)/ALGBRC(1) ** CONSTS(74))/(CONSTS(75) ** CONSTS(74)+STATES(6)/ALGBRC(1) ** CONSTS(74))
      ALGBRC(19) = TERNRY(ALGBRC(17).GE.0.00000, ALGBRC(17), ALGBRC(17))
      RATES(6) = (ALGBRC(5)+ALGBRC(15)) - ALGBRC(19)
      ALGBRC(6) = ( (CONSTS(23)/( CONSTS(21)*CONSTS(22)))*( STATES(1)*STATES(7) - ( CONSTS(8)*ALGBRC(2))/CONSTS(18)))/( (1.00000+STATES(1)/CONSTS(21)+CONSTS(8)/CONSTS(19))*(1.00000+STATES(7)/CONSTS(22)+ALGBRC(2)/CONSTS(20)))
      ALGBRC(7) = TERNRY(ALGBRC(6).GE.0.00000, ALGBRC(6), ALGBRC(6))
      ALGBRC(20) = ( (( CONSTS(62)*STATES(1))/CONSTS(61))*(1.00000 - STATES(4)/( STATES(1)*CONSTS(59)))*STATES(1)/CONSTS(61)+STATES(4)/CONSTS(60) ** CONSTS(58) - 1.00000)/(1.00000+STATES(1)/CONSTS(61)+STATES(4)/CONSTS(60) ** CONSTS(58))
      ALGBRC(21) = TERNRY(ALGBRC(20).GE.0.00000, ALGBRC(20), ALGBRC(20))
      RATES(1) = ALGBRC(5) - (ALGBRC(7)+ALGBRC(9)+ALGBRC(21))
      ALGBRC(22) = (( CONSTS(65)*STATES(4))/CONSTS(63))/(1.00000+STATES(4)/CONSTS(63)+STATES(5)/CONSTS(64))
      ALGBRC(23) = TERNRY(ALGBRC(22).GE.0.00000, ALGBRC(22), ALGBRC(22))
      ALGBRC(28) =  CONSTS(81)*STATES(4)
      RATES(4) =  0.500000*ALGBRC(21) - (ALGBRC(23)+ALGBRC(28))
      ALGBRC(26) = ( (CONSTS(73)/( CONSTS(72)*CONSTS(70)))*( STATES(5)*STATES(7) - ( CONSTS(5)*ALGBRC(2))/CONSTS(68)))/( (1.00000+STATES(5)/CONSTS(72)+CONSTS(5)/CONSTS(71))*(1.00000+STATES(7)/CONSTS(70)+ALGBRC(2)/CONSTS(69)))
      ALGBRC(27) = TERNRY(ALGBRC(26).GE.0.00000, ALGBRC(26), ALGBRC(26))
      ALGBRC(24) = (( CONSTS(67)*STATES(5))/CONSTS(66))/(1.00000+STATES(5)/CONSTS(66))
      ALGBRC(25) = TERNRY(ALGBRC(24).GE.0.00000, ALGBRC(24), ALGBRC(24))
      RATES(5) = (ALGBRC(23)+ALGBRC(28)) - (ALGBRC(27)+ALGBRC(25))
      ALGBRC(29) = (( CONSTS(80)*STATES(7)*CONSTS(9))/( CONSTS(78)*CONSTS(79)))/( (1.00000+STATES(7)/CONSTS(78)+ALGBRC(2)/CONSTS(77))*(1.00000+CONSTS(9)/CONSTS(79)))
      ALGBRC(30) = TERNRY(ALGBRC(29).GE.0.00000, ALGBRC(29), ALGBRC(29))
      RATES(7) = (ALGBRC(5)+ALGBRC(9)) - (ALGBRC(7)+ALGBRC(14)+ALGBRC(18)+ALGBRC(27)+ALGBRC(30))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = CONSTS(3) - STATES(8)
      ALGBRC(10) = ( (CONSTS(38)/( CONSTS(37)*CONSTS(32)))*( STATES(8)*CONSTS(10) - ( STATES(2)*ALGBRC(3))/CONSTS(31)))/(1.00000+STATES(8)/CONSTS(37)+CONSTS(10)/CONSTS(33)+STATES(2)/CONSTS(36)+ALGBRC(3)/CONSTS(34)+( STATES(8)*CONSTS(10))/( CONSTS(37)*CONSTS(32))+( STATES(2)*ALGBRC(3))/( CONSTS(35)*CONSTS(34)))
      ALGBRC(11) = TERNRY(ALGBRC(10).GE.0.00000, ALGBRC(10), ALGBRC(10))
      ALGBRC(1) = CONSTS(1) - STATES(6)
      ALGBRC(13) = ( (CONSTS(44)/( CONSTS(42)*CONSTS(43)))*( STATES(2)*ALGBRC(1) - ( CONSTS(4)*STATES(6))/CONSTS(39)))/( (1.00000+STATES(2)/CONSTS(43)+CONSTS(4)/CONSTS(40))*(1.00000+ALGBRC(1)/CONSTS(42)+STATES(6)/CONSTS(41)))
      ALGBRC(15) = TERNRY(ALGBRC(13).GE.0.00000, ALGBRC(13), ALGBRC(13))
      ALGBRC(2) = CONSTS(2) - STATES(7)
      ALGBRC(8) = (( (( (( (CONSTS(30)/(1.00000+( CONSTS(24)*STATES(7))/ALGBRC(2)))*STATES(1))/CONSTS(27))*ALGBRC(2))/CONSTS(25))*ALGBRC(3))/CONSTS(26))/( (1.00000+STATES(1)/CONSTS(27))*(1.00000+ALGBRC(2)/CONSTS(25)+STATES(7)/CONSTS(28))*(1.00000+ALGBRC(3)/CONSTS(26)+STATES(8)/CONSTS(29)))
      ALGBRC(9) = TERNRY(ALGBRC(8).GE.0.00000, ALGBRC(8), ALGBRC(8))
      ALGBRC(12) = ( (CONSTS(51)/( CONSTS(49)*CONSTS(47)))*( STATES(8)*STATES(7) - ( ALGBRC(2)*ALGBRC(3)*STATES(3))/CONSTS(45)))/( (1.00000+ALGBRC(2)/CONSTS(46)+STATES(7)/CONSTS(47))*(1.00000+STATES(8)/CONSTS(49)+ALGBRC(3)/CONSTS(48))*(1.00000+STATES(3)/CONSTS(50)))
      ALGBRC(14) = TERNRY(ALGBRC(12).GE.0.00000, ALGBRC(12), ALGBRC(12))
      ALGBRC(16) = ( (CONSTS(57)/( CONSTS(56)*CONSTS(54)))*( STATES(3)*STATES(7) - ( CONSTS(6)*ALGBRC(2))/CONSTS(52)))/( (1.00000+ALGBRC(2)/CONSTS(53)+STATES(7)/CONSTS(54))*(1.00000+STATES(3)/CONSTS(56)+CONSTS(6)/CONSTS(55)))
      ALGBRC(18) = TERNRY(ALGBRC(16).GE.0.00000, ALGBRC(16), ALGBRC(16))
      ALGBRC(4) = (( (( (( CONSTS(17)*CONSTS(7))/CONSTS(11))*ALGBRC(2))/CONSTS(12))*ALGBRC(1))/CONSTS(13))/( (1.00000+CONSTS(7)/CONSTS(11)+STATES(1)/CONSTS(14))*(1.00000+ALGBRC(2)/CONSTS(12)+STATES(7)/CONSTS(15))*(1.00000+ALGBRC(1)/CONSTS(13)+STATES(6)/CONSTS(16)))
      ALGBRC(5) = TERNRY(ALGBRC(4).GE.0.00000, ALGBRC(4), ALGBRC(4))
      ALGBRC(17) = ( CONSTS(76)*STATES(6)/ALGBRC(1) ** CONSTS(74))/(CONSTS(75) ** CONSTS(74)+STATES(6)/ALGBRC(1) ** CONSTS(74))
      ALGBRC(19) = TERNRY(ALGBRC(17).GE.0.00000, ALGBRC(17), ALGBRC(17))
      ALGBRC(6) = ( (CONSTS(23)/( CONSTS(21)*CONSTS(22)))*( STATES(1)*STATES(7) - ( CONSTS(8)*ALGBRC(2))/CONSTS(18)))/( (1.00000+STATES(1)/CONSTS(21)+CONSTS(8)/CONSTS(19))*(1.00000+STATES(7)/CONSTS(22)+ALGBRC(2)/CONSTS(20)))
      ALGBRC(7) = TERNRY(ALGBRC(6).GE.0.00000, ALGBRC(6), ALGBRC(6))
      ALGBRC(20) = ( (( CONSTS(62)*STATES(1))/CONSTS(61))*(1.00000 - STATES(4)/( STATES(1)*CONSTS(59)))*STATES(1)/CONSTS(61)+STATES(4)/CONSTS(60) ** CONSTS(58) - 1.00000)/(1.00000+STATES(1)/CONSTS(61)+STATES(4)/CONSTS(60) ** CONSTS(58))
      ALGBRC(21) = TERNRY(ALGBRC(20).GE.0.00000, ALGBRC(20), ALGBRC(20))
      ALGBRC(22) = (( CONSTS(65)*STATES(4))/CONSTS(63))/(1.00000+STATES(4)/CONSTS(63)+STATES(5)/CONSTS(64))
      ALGBRC(23) = TERNRY(ALGBRC(22).GE.0.00000, ALGBRC(22), ALGBRC(22))
      ALGBRC(28) =  CONSTS(81)*STATES(4)
      ALGBRC(26) = ( (CONSTS(73)/( CONSTS(72)*CONSTS(70)))*( STATES(5)*STATES(7) - ( CONSTS(5)*ALGBRC(2))/CONSTS(68)))/( (1.00000+STATES(5)/CONSTS(72)+CONSTS(5)/CONSTS(71))*(1.00000+STATES(7)/CONSTS(70)+ALGBRC(2)/CONSTS(69)))
      ALGBRC(27) = TERNRY(ALGBRC(26).GE.0.00000, ALGBRC(26), ALGBRC(26))
      ALGBRC(24) = (( CONSTS(67)*STATES(5))/CONSTS(66))/(1.00000+STATES(5)/CONSTS(66))
      ALGBRC(25) = TERNRY(ALGBRC(24).GE.0.00000, ALGBRC(24), ALGBRC(24))
      ALGBRC(29) = (( CONSTS(80)*STATES(7)*CONSTS(9))/( CONSTS(78)*CONSTS(79)))/( (1.00000+STATES(7)/CONSTS(78)+ALGBRC(2)/CONSTS(77))*(1.00000+CONSTS(9)/CONSTS(79)))
      ALGBRC(30) = TERNRY(ALGBRC(29).GE.0.00000, ALGBRC(29), ALGBRC(29))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END