C
C There are a total of 35 entries in the algebraic variable array.
C There are a total of 16 entries in each of the rate and state variable arrays.
C There are a total of 134 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is ADP_m in component ADP_m (millimolar).
C ALGBRC(33) is V_ANT in component calcium_dynamics (flux).
C ALGBRC(31) is V_ATPase in component oxidative_phosphorylation (flux).
C ALGBRC(12) is V_SL in component V_SL (flux).
C STATES(2) is NADH in component NADH (millimolar).
C ALGBRC(27) is V_O2 in component oxidative_phosphorylation (flux).
C ALGBRC(10) is V_IDH in component V_IDH (flux).
C ALGBRC(11) is V_KGDH in component V_KGDH (flux).
C ALGBRC(14) is V_MDH in component V_MDH (flux).
C STATES(3) is ISOC in component ISOC (millimolar).
C ALGBRC(9) is V_ACO in component V_ACO (flux).
C STATES(4) is alpha_KG in component alpha_KG (millimolar).
C ALGBRC(16) is V_AAT in component V_AAT (flux).
C STATES(5) is SCoA in component SCoA (millimolar).
C STATES(6) is Suc in component Suc (millimolar).
C ALGBRC(13) is V_SDH in component V_SDH (flux).
C STATES(7) is FUM in component FUM (millimolar).
C ALGBRC(17) is V_FH in component V_FH (flux).
C STATES(8) is MAL in component MAL (millimolar).
C STATES(9) is OAA in component OAA (millimolar).
C ALGBRC(8) is V_CS in component V_CS (flux).
C STATES(10) is ASP in component ASP (millimolar).
C ALGBRC(19) is V_C_ASP in component V_C_ASP (flux).
C STATES(11) is Ca_m in component Ca_m (micromolar).
C CONSTS(1) is f in component Ca_m (dimensionless).
C ALGBRC(34) is V_uni in component calcium_dynamics (flux).
C ALGBRC(35) is V_NaCa in component calcium_dynamics (flux).
C ALGBRC(2) is Ca_i in component Ca_i (micromolar).
C CONSTS(2) is stim_start in component Ca_i (second).
C CONSTS(3) is stim_end in component Ca_i (second).
C ALGBRC(1) is stim_period in component Ca_i (second).
C CONSTS(4) is stim_duration in component Ca_i (second).
C CONSTS(5) is pulse_value in component Ca_i (micromolar).
C CONSTS(6) is Na_i in component Na_i (millimolar).
C CONSTS(7) is ATP_i in component ATP_i (millimolar).
C ALGBRC(3) is ATP_m in component ATP_m (millimolar).
C CONSTS(8) is Cm in component ATP_m (millimolar).
C ALGBRC(4) is ADP_i in component ADP_i (millimolar).
C CONSTS(9) is pulse_value in component ADP_i (millimolar).
C CONSTS(10) is GLU in component GLU (millimolar).
C CONSTS(11) is Mg in component Mg (millimolar).
C CONSTS(12) is H in component H (millimolar).
C CONSTS(13) is Pi in component Pi (millimolar).
C CONSTS(14) is CoA in component CoA (millimolar).
C CONSTS(15) is AcCoA in component AcCoA (millimolar).
C CONSTS(16) is FAD in component FAD (millimolar).
C CONSTS(17) is FADH2 in component FADH2 (millimolar).
C ALGBRC(5) is NAD in component NAD (millimolar).
C CONSTS(18) is C_PN in component NAD (millimolar).
C CONSTS(19) is NADPH in component NADPH (millimolar).
C STATES(12) is O2_m in component O2_m (millimolar).
C CONSTS(20) is shunt in component O2_m (dimensionless).
C ALGBRC(21) is VTr_ROS in component VTr_ROS (flux).
C STATES(13) is O2_i in component O2_i (millimolar).
C ALGBRC(15) is V_SOD in component V_SOD (flux).
C STATES(14) is H2O2 in component H2O2 (millimolar).
C ALGBRC(20) is V_CAT in component V_CAT (flux).
C ALGBRC(22) is V_GPX in component V_GPX (flux).
C STATES(15) is GSH in component GSH (millimolar).
C ALGBRC(24) is V_GR in component V_GR (flux).
C ALGBRC(6) is GSSG in component GSSG (millimolar).
C CONSTS(21) is GT in component GSSG (millimolar).
C ALGBRC(7) is CIT in component CIT (millimolar).
C CONSTS(22) is C_Kint in component CIT (millimolar).
C STATES(16) is delta_psi_m in component mitochondrial_membrane (volt).
C CONSTS(23) is R in component mitochondrial_membrane (volt_coulomb_per_mole_kelvin).
C CONSTS(24) is T in component mitochondrial_membrane (kelvin).
C CONSTS(25) is F in component mitochondrial_membrane (coulomb_per_mole).
C CONSTS(26) is C_mito in component mitochondrial_membrane (millimolar_per_volt).
C ALGBRC(26) is V_He in component oxidative_phosphorylation (flux).
C ALGBRC(28) is V_He_F in component oxidative_phosphorylation (flux).
C ALGBRC(32) is V_Hu in component oxidative_phosphorylation (flux).
C ALGBRC(29) is V_Hleak in component oxidative_phosphorylation (flux).
C CONSTS(27) is Km_AcCoA in component V_CS (millimolar).
C CONSTS(28) is Km_OAA in component V_CS (millimolar).
C CONSTS(29) is Kcat_CS in component V_CS (first_order_rate_constant).
C CONSTS(30) is ET_CS in component V_CS (millimolar).
C CONSTS(31) is Kf_ACO in component V_ACO (first_order_rate_constant).
C CONSTS(32) is KE_ACO in component V_ACO (dimensionless).
C CONSTS(33) is Kh_1 in component V_IDH (millimolar).
C CONSTS(34) is Kh_2 in component V_IDH (millimolar).
C CONSTS(35) is Km_ISOC in component V_IDH (millimolar).
C CONSTS(36) is Ka_ADP in component V_IDH (millimolar).
C CONSTS(37) is Ka_Ca in component V_IDH (micromolar).
C CONSTS(38) is Km_NAD in component V_IDH (millimolar).
C CONSTS(39) is Ki_NADH in component V_IDH (millimolar).
C CONSTS(40) is Kcat_IDH in component V_IDH (first_order_rate_constant).
C CONSTS(41) is ET_IDH in component V_IDH (millimolar).
C CONSTS(42) is ni in component V_IDH (dimensionless).
C CONSTS(43) is Km_alpha_KG in component V_KGDH (millimolar).
C CONSTS(44) is Kcat_KGDH in component V_KGDH (first_order_rate_constant).
C CONSTS(45) is ET_KGDH in component V_KGDH (millimolar).
C CONSTS(46) is Kd_Mg in component V_KGDH (millimolar).
C CONSTS(47) is Kd_Ca in component V_KGDH (micromolar).
C CONSTS(48) is n_alpha_KG in component V_KGDH (dimensionless).
C CONSTS(49) is Km_NAD in component V_KGDH (millimolar).
C CONSTS(50) is kf_SL in component V_SL (second_order_rate_constant).
C CONSTS(51) is Ke_SL in component V_SL (millimolar).
C CONSTS(52) is Kisdh_OAA in component V_SDH (millimolar).
C CONSTS(53) is Kcat_SDH in component V_SDH (first_order_rate_constant).
C CONSTS(54) is ET_SDH in component V_SDH (millimolar).
C CONSTS(55) is Km_Suc in component V_SDH (millimolar).
C CONSTS(56) is Ki_FUM in component V_SDH (millimolar).
C CONSTS(57) is Km_MAL in component V_MDH (millimolar).
C CONSTS(58) is Kcat_MDH in component V_MDH (first_order_rate_constant).
C CONSTS(59) is ET_MDH in component V_MDH (millimolar).
C CONSTS(60) is Ki_OAA in component V_MDH (millimolar).
C CONSTS(133) is fh_a in component V_MDH (dimensionless).
C CONSTS(134) is fh_i in component V_MDH (dimensionless).
C CONSTS(61) is Km_NAD in component V_MDH (millimolar).
C CONSTS(62) is kh1 in component V_MDH (millimolar).
C CONSTS(63) is kh2 in component V_MDH (millimolar).
C CONSTS(64) is kh3 in component V_MDH (millimolar).
C CONSTS(65) is kh4 in component V_MDH (millimolar).
C CONSTS(66) is k_offset in component V_MDH (dimensionless).
C CONSTS(67) is Ke_FH in component V_FH (dimensionless).
C CONSTS(68) is kf_FH in component V_FH (first_order_rate_constant).
C CONSTS(69) is Ke_AAT in component V_AAT (dimensionless).
C CONSTS(70) is kf_AAT in component V_AAT (second_order_rate_constant).
C CONSTS(71) is k_C_ASP in component V_C_ASP (first_order_rate_constant).
C CONSTS(72) is k1_SOD in component V_SOD (second_order_rate_constant).
C CONSTS(73) is k3_SOD in component V_SOD (second_order_rate_constant).
C CONSTS(74) is k5_SOD in component V_SOD (first_order_rate_constant).
C CONSTS(75) is ET_SOD in component V_SOD (millimolar).
C CONSTS(76) is Ki_H2O2 in component V_SOD (millimolar).
C CONSTS(77) is k1_CAT in component V_CAT (second_order_rate_constant).
C CONSTS(78) is ET_CAT in component V_CAT (millimolar).
C CONSTS(79) is fr in component V_CAT (per_millimolar).
C CONSTS(80) is phi1_GPX in component V_GPX (millimolar_second).
C CONSTS(81) is phi2_GPX in component V_GPX (millimolar_second).
C CONSTS(82) is ET_GPX in component V_GPX (millimolar).
C CONSTS(83) is k1_GR in component V_GR (first_order_rate_constant).
C CONSTS(84) is KM_GSSG in component V_GR (millimolar).
C CONSTS(85) is KM_NADPH in component V_GR (millimolar).
C CONSTS(86) is ET_GR in component V_GR (millimolar).
C ALGBRC(18) is V_IMAC in component V_IMAC (flux).
C CONSTS(87) is Gmax in component V_IMAC (millimolar_per_second_volt).
C CONSTS(88) is GL in component V_IMAC (millimolar_per_second_volt).
C CONSTS(89) is a in component V_IMAC (dimensionless).
C CONSTS(90) is b in component V_IMAC (dimensionless).
C CONSTS(91) is kappa in component V_IMAC (per_volt).
C CONSTS(92) is delta_psi_bm in component V_IMAC (volt).
C CONSTS(93) is Kcc in component V_IMAC (millimolar).
C CONSTS(94) is j in component VTr_ROS (dimensionless).
C CONSTS(95) is rho_res in component oxidative_phosphorylation (millimolar).
C CONSTS(96) is rho_res_F in component oxidative_phosphorylation (millimolar).
C CONSTS(97) is ra in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(98) is rc1 in component oxidative_phosphorylation (first_order_rate_constant).
C ALGBRC(23) is Ares in component oxidative_phosphorylation (volt).
C CONSTS(132) is Ares_F in component oxidative_phosphorylation (volt).
C CONSTS(99) is r1 in component oxidative_phosphorylation (dimensionless).
C CONSTS(100) is r2 in component oxidative_phosphorylation (dimensionless).
C CONSTS(101) is r3 in component oxidative_phosphorylation (dimensionless).
C CONSTS(102) is rb in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(103) is rc2 in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(104) is Kres in component oxidative_phosphorylation (dimensionless).
C CONSTS(105) is Kres_F in component oxidative_phosphorylation (dimensionless).
C CONSTS(106) is gH in component oxidative_phosphorylation (millimolar_per_second_volt).
C CONSTS(107) is delta_psi_B in component oxidative_phosphorylation (volt).
C CONSTS(108) is g in component oxidative_phosphorylation (dimensionless).
C ALGBRC(25) is delta_mu_H in component oxidative_phosphorylation (volt).
C CONSTS(109) is delta_pH in component oxidative_phosphorylation (dimensionless).
C CONSTS(110) is rho_F1 in component oxidative_phosphorylation (millimolar).
C CONSTS(111) is pa in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(112) is pc1 in component oxidative_phosphorylation (first_order_rate_constant).
C ALGBRC(30) is AF1 in component oxidative_phosphorylation (volt).
C CONSTS(113) is p1 in component oxidative_phosphorylation (dimensionless).
C CONSTS(114) is p2 in component oxidative_phosphorylation (dimensionless).
C CONSTS(115) is p3 in component oxidative_phosphorylation (dimensionless).
C CONSTS(116) is pb in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(117) is pc2 in component oxidative_phosphorylation (first_order_rate_constant).
C CONSTS(118) is KF1 in component oxidative_phosphorylation (millimolar).
C CONSTS(119) is h in component calcium_dynamics (dimensionless).
C CONSTS(120) is delta_psi_0 in component calcium_dynamics (volt).
C CONSTS(121) is Vmax_ANT in component calcium_dynamics (flux).
C CONSTS(122) is L in component calcium_dynamics (dimensionless).
C CONSTS(123) is na in component calcium_dynamics (dimensionless).
C CONSTS(124) is Vmax_uni in component calcium_dynamics (flux).
C CONSTS(125) is K_act in component calcium_dynamics (micromolar).
C CONSTS(126) is K_trans in component calcium_dynamics (micromolar).
C CONSTS(127) is n in component calcium_dynamics (dimensionless).
C CONSTS(128) is Vmax_NaCa in component calcium_dynamics (flux).
C CONSTS(129) is KNa in component calcium_dynamics (millimolar).
C CONSTS(130) is KCa in component calcium_dynamics (micromolar).
C CONSTS(131) is b in component calcium_dynamics (dimensionless).
C RATES(1) is d/dt ADP_m in component ADP_m (millimolar).
C RATES(2) is d/dt NADH in component NADH (millimolar).
C RATES(3) is d/dt ISOC in component ISOC (millimolar).
C RATES(4) is d/dt alpha_KG in component alpha_KG (millimolar).
C RATES(5) is d/dt SCoA in component SCoA (millimolar).
C RATES(6) is d/dt Suc in component Suc (millimolar).
C RATES(7) is d/dt FUM in component FUM (millimolar).
C RATES(8) is d/dt MAL in component MAL (millimolar).
C RATES(9) is d/dt OAA in component OAA (millimolar).
C RATES(10) is d/dt ASP in component ASP (millimolar).
C RATES(11) is d/dt Ca_m in component Ca_m (micromolar).
C RATES(12) is d/dt O2_m in component O2_m (millimolar).
C RATES(13) is d/dt O2_i in component O2_i (millimolar).
C RATES(14) is d/dt H2O2 in component H2O2 (millimolar).
C RATES(15) is d/dt GSH in component GSH (millimolar).
C RATES(16) is d/dt delta_psi_m in component mitochondrial_membrane (volt).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.1
      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) = 0.01
      STATES(9) = 0.01
      STATES(10) = 0.01
      STATES(11) = 0.01
      CONSTS(1) = 0.0003
      CONSTS(2) = 0
      CONSTS(3) = 10000
      CONSTS(4) = 0.4
      CONSTS(5) = 1
      CONSTS(6) = 10.0
      CONSTS(7) = 6.5
      CONSTS(8) = 15.0
      CONSTS(9) = 0.15
      CONSTS(10) = 20
      CONSTS(11) = 0.4
      CONSTS(12) = 2.5E-5
      CONSTS(13) = 20.0
      CONSTS(14) = 0.02
      CONSTS(15) = 0.0002
      CONSTS(16) = 0.01
      CONSTS(17) = 1.24
      CONSTS(18) = 10.0
      CONSTS(19) = 1.0
      STATES(12) = 0.01
      CONSTS(20) = 0.05
      STATES(13) = 0.01
      STATES(14) = 0.01
      STATES(15) = 0.01
      CONSTS(21) = 1.0
      CONSTS(22) = 1.0
      STATES(16) = 0.01
      CONSTS(23) = 8.315
      CONSTS(24) = 310.16
      CONSTS(25) = 96480
      CONSTS(26) = 1.812
      CONSTS(27) = 1.26E-2
      CONSTS(28) = 6.4E-4
      CONSTS(29) = 3.2
      CONSTS(30) = 0.4
      CONSTS(31) = 12.5
      CONSTS(32) = 2.22
      CONSTS(33) = 8.1E-5
      CONSTS(34) = 5.98E-5
      CONSTS(35) = 1.52
      CONSTS(36) = 6.2E-2
      CONSTS(37) = 1.41
      CONSTS(38) = 0.923
      CONSTS(39) = 0.19
      CONSTS(40) = 1.94
      CONSTS(41) = 0.109
      CONSTS(42) = 1
      CONSTS(43) = 1.94
      CONSTS(44) = 0.15
      CONSTS(45) = 0.5
      CONSTS(46) = 0.0308
      CONSTS(47) = 1.27
      CONSTS(48) = 1.2
      CONSTS(49) = 38.7
      CONSTS(50) = 0.127
      CONSTS(51) = 3.115
      CONSTS(52) = 0.15
      CONSTS(53) = 1.0
      CONSTS(54) = 0.5
      CONSTS(55) = 3.0E-2
      CONSTS(56) = 1.3
      CONSTS(57) = 1.493
      CONSTS(58) = 2.775E1
      CONSTS(59) = 0.154
      CONSTS(60) = 3.1E-3
      CONSTS(61) = 0.2244
      CONSTS(62) = 1.13E-5
      CONSTS(63) = 26.7
      CONSTS(64) = 6.68E-9
      CONSTS(65) = 5.62E-6
      CONSTS(66) = 3.99E-2
      CONSTS(67) = 1.0
      CONSTS(68) = 0.83
      CONSTS(69) = 6.6
      CONSTS(70) = 0.644
      CONSTS(71) = 0.01
      CONSTS(72) = 2.4E6
      CONSTS(73) = 4.8E4
      CONSTS(74) = 5.0E-1
      CONSTS(75) = 0.4E-3
      CONSTS(76) = 0.5
      CONSTS(77) = 1.7E4
      CONSTS(78) = 0.001
      CONSTS(79) = 50.0
      CONSTS(80) = 0.15
      CONSTS(81) = 0.5
      CONSTS(82) = 0.00141
      CONSTS(83) = 0.0308
      CONSTS(84) = 1.94
      CONSTS(85) = 38.7
      CONSTS(86) = 1.27E-3
      CONSTS(87) = 7.82
      CONSTS(88) = 0.0782
      CONSTS(89) = 1.0E-3
      CONSTS(90) = 1.0E4
      CONSTS(91) = 70.0
      CONSTS(92) = 0.004
      CONSTS(93) = 0.01
      CONSTS(94) = 0.12
      CONSTS(95) = 0.0006
      CONSTS(96) = 0.0045
      CONSTS(97) = 6.394E-10
      CONSTS(98) = 2.656E-19
      CONSTS(99) = 2.077E-18
      CONSTS(100) = 1.728E-9
      CONSTS(101) = 1.059E-26
      CONSTS(102) = 1.762E-13
      CONSTS(103) = 8.632E-27
      CONSTS(104) = 1.35E18
      CONSTS(105) = 5.765E13
      CONSTS(106) = 0.01
      CONSTS(107) = 0.05
      CONSTS(108) = 0.85
      CONSTS(109) = -0.6
      CONSTS(110) = 0.06
      CONSTS(111) = 1.656E-5
      CONSTS(112) = 9.651E-14
      CONSTS(113) = 1.346E-8
      CONSTS(114) = 7.739E-7
      CONSTS(115) = 6.65E-15
      CONSTS(116) = 3.373E-7
      CONSTS(117) = 4.585E-14
      CONSTS(118) = 1.71E6
      CONSTS(119) = 0.5
      CONSTS(120) = 0.091
      CONSTS(121) = 0.05
      CONSTS(122) = 110.0
      CONSTS(123) = 2.8
      CONSTS(124) = 0.000625
      CONSTS(125) = 0.38
      CONSTS(126) = 19.0
      CONSTS(127) = 3.0
      CONSTS(128) = 0.005
      CONSTS(129) = 9.4
      CONSTS(130) = 3.75E-1
      CONSTS(131) = 0.5
      CONSTS(132) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(105)*CONSTS(17)/CONSTS(16) ** 0.500000)
      CONSTS(133) = 1.00000/(1.00000+CONSTS(12)/CONSTS(62)+CONSTS(12) ** 2.00000/( CONSTS(62)*CONSTS(63)))+CONSTS(66)
      CONSTS(134) = 1.00000/(1.00000+CONSTS(64)/CONSTS(12)+( CONSTS(64)*CONSTS(65))/CONSTS(12) ** 2.00000) ** 2.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(18) - STATES(2)
      ALGBRC(10) = ( CONSTS(40)*CONSTS(41))/(1.00000+CONSTS(12)/CONSTS(33)+CONSTS(34)/CONSTS(12)+CONSTS(35)/STATES(3) ** CONSTS(42)/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))+ (CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39))+( CONSTS(35)/STATES(3) ** CONSTS(42)*(CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39)))/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37))))
      ALGBRC(7) = CONSTS(22) - (STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9))
      ALGBRC(9) =  CONSTS(31)*(ALGBRC(7) - STATES(3)/CONSTS(32))
      RATES(3) = ALGBRC(9) - ALGBRC(10)
      ALGBRC(3) = CONSTS(8) - STATES(1)
      ALGBRC(12) =  CONSTS(50)*( STATES(5)*STATES(1) - ( STATES(6)*ALGBRC(3)*CONSTS(14))/CONSTS(51))
      ALGBRC(11) = ( CONSTS(44)*CONSTS(45))/(1.00000+CONSTS(43)/STATES(4) ** CONSTS(48)/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))+(CONSTS(49)/ALGBRC(5))/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47))))
      RATES(5) = ALGBRC(11) - ALGBRC(12)
      ALGBRC(13) = ( CONSTS(53)*CONSTS(54))/(1.00000+ (CONSTS(55)/STATES(6))*(1.00000+STATES(9)/CONSTS(52))*(1.00000+STATES(7)/CONSTS(56)))
      RATES(6) = ALGBRC(12) - ALGBRC(13)
      ALGBRC(16) =  CONSTS(70)*( STATES(9)*CONSTS(10) - ( STATES(4)*STATES(10))/CONSTS(69))
      RATES(4) = ALGBRC(16)+ALGBRC(10)+- ALGBRC(11)
      ALGBRC(17) =  CONSTS(68)*(STATES(7) - STATES(8)/CONSTS(67))
      RATES(7) = ALGBRC(13) - ALGBRC(17)
      ALGBRC(14) = ( CONSTS(58)*CONSTS(59)*CONSTS(133)*CONSTS(134))/(1.00000+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))+CONSTS(61)/ALGBRC(5)+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))*(CONSTS(61)/ALGBRC(5)))
      RATES(8) = ALGBRC(17) - ALGBRC(14)
      ALGBRC(8) = ( CONSTS(29)*CONSTS(30))/(1.00000+CONSTS(27)/CONSTS(15)+CONSTS(28)/STATES(9)+ (CONSTS(27)/CONSTS(15))*(CONSTS(28)/STATES(9)))
      RATES(9) = ALGBRC(14) - (ALGBRC(8)+ALGBRC(16))
      ALGBRC(19) =  CONSTS(71)*STATES(10)
      RATES(10) = ALGBRC(16) - ALGBRC(19)
      ALGBRC(18) =  (CONSTS(89)+CONSTS(90)/(1.00000+CONSTS(93)/STATES(13)))*(CONSTS(88)+CONSTS(87)/(1.00000+EXP( CONSTS(91)*(CONSTS(92) - STATES(16)))))*STATES(16)
      ALGBRC(21) =  CONSTS(94)*(ALGBRC(18)/STATES(16))*(STATES(16) -  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*arbitrary_log(STATES(12)/STATES(13), 10))
      ALGBRC(15) = ( 2.00000*CONSTS(72)*CONSTS(74)*(CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))*CONSTS(75)*STATES(13))/( CONSTS(74)*( 2.00000*CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))+ STATES(13)*CONSTS(72)*CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))
      RATES(13) = ALGBRC(21) - ALGBRC(15)
      ALGBRC(20) =  2.00000*CONSTS(77)*CONSTS(78)*STATES(14)*EXP( - CONSTS(79)*STATES(14))
      ALGBRC(22) = ( CONSTS(82)*STATES(14)*STATES(15))/( CONSTS(80)*STATES(15)+ CONSTS(81)*STATES(14))
      RATES(14) = ALGBRC(15) - (ALGBRC(20)+ALGBRC(22))
      ALGBRC(6) =  2.00000*(CONSTS(21) - STATES(15))
      ALGBRC(24) = ( CONSTS(83)*CONSTS(86))/(1.00000+CONSTS(84)/ALGBRC(6)+CONSTS(85)/CONSTS(19)+ (CONSTS(84)/ALGBRC(6))*(CONSTS(85)/CONSTS(19)))
      RATES(15) = ALGBRC(24) - ALGBRC(22)
      ALGBRC(23) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(104)*STATES(2)/ALGBRC(5) ** 0.500000)
      ALGBRC(25) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*CONSTS(109)+STATES(16)
      ALGBRC(27) =  0.500000*CONSTS(95)*((( (CONSTS(97)+ CONSTS(98)*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) -  CONSTS(97)*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))+ CONSTS(103)*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      RATES(2) = - ALGBRC(27)+ALGBRC(10)+ALGBRC(11)+ALGBRC(14)
      RATES(12) =  CONSTS(20)*ALGBRC(27) - ALGBRC(21)
      ALGBRC(4) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, CONSTS(9), 0.0500000)
      ALGBRC(33) =  CONSTS(121)*((1.00000 - ( 0.0500000*CONSTS(7)*0.450000*0.800000*STATES(1))/( 0.450000*ALGBRC(4)*0.0500000*ALGBRC(3)))/( (1.00000+ (( 0.0500000*CONSTS(7))/( 0.450000*ALGBRC(4)))*EXP(( - CONSTS(119)*CONSTS(25)*CONSTS(120))/( CONSTS(23)*CONSTS(24))))*(1.00000+( 0.450000*0.800000*STATES(1))/( 0.0500000*ALGBRC(3)))))
      ALGBRC(30) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(118)*(ALGBRC(3)/( STATES(1)*CONSTS(13))))
      ALGBRC(31) =  - CONSTS(110)*(( ( 100.000*CONSTS(111)+ CONSTS(112)*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - ( CONSTS(111)*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))+ CONSTS(117)*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      RATES(1) = ALGBRC(33) - (ALGBRC(31)+ALGBRC(12))
      ALGBRC(1) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, 0.500000, 4.00000)
      ALGBRC(2) = TERNRY(VOI.GE.CONSTS(2).AND.VOI.LE.CONSTS(3).AND.(VOI - CONSTS(2)) -  INT((VOI - CONSTS(2))/ALGBRC(1))*ALGBRC(1).LE.CONSTS(4), CONSTS(5), 0.100000)
      ALGBRC(34) =  CONSTS(124)*(( (ALGBRC(2)/CONSTS(126))*1.00000+ALGBRC(2)/CONSTS(126) ** 3.00000*(( 2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))/(1.00000+ALGBRC(2)/CONSTS(126) ** 4.00000+ (CONSTS(122)/1.00000+ALGBRC(2)/CONSTS(125) ** CONSTS(123))*(1.00000 - EXP(( -2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))))
      ALGBRC(35) =  CONSTS(128)*(( EXP(( CONSTS(131)*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))*EXP(log(ALGBRC(2)/STATES(11))))/( 1.00000+CONSTS(129)/CONSTS(6) ** CONSTS(127)*(1.00000+CONSTS(130)/STATES(11))))
      RATES(11) =  CONSTS(1)*1.00000*(ALGBRC(34) - ALGBRC(35))
      ALGBRC(26) =  6.00000*CONSTS(95)*(( CONSTS(97)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))) -  (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(28) =  6.00000*CONSTS(96)*(( CONSTS(97)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))) -  (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(32) =  -3.00000*CONSTS(110)*(( 100.000*CONSTS(111)*(1.00000+EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24)))) -  (CONSTS(111)+CONSTS(116))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(29) =  CONSTS(106)*ALGBRC(25)
      RATES(16) = (ALGBRC(26)+ALGBRC(28)+- (ALGBRC(32)+ALGBRC(33)+ALGBRC(29)+ALGBRC(35)+ 2.00000*ALGBRC(34)))/CONSTS(26)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(18) - STATES(2)
      ALGBRC(10) = ( CONSTS(40)*CONSTS(41))/(1.00000+CONSTS(12)/CONSTS(33)+CONSTS(34)/CONSTS(12)+CONSTS(35)/STATES(3) ** CONSTS(42)/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))+ (CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39))+( CONSTS(35)/STATES(3) ** CONSTS(42)*(CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39)))/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37))))
      ALGBRC(7) = CONSTS(22) - (STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9))
      ALGBRC(9) =  CONSTS(31)*(ALGBRC(7) - STATES(3)/CONSTS(32))
      ALGBRC(3) = CONSTS(8) - STATES(1)
      ALGBRC(12) =  CONSTS(50)*( STATES(5)*STATES(1) - ( STATES(6)*ALGBRC(3)*CONSTS(14))/CONSTS(51))
      ALGBRC(11) = ( CONSTS(44)*CONSTS(45))/(1.00000+CONSTS(43)/STATES(4) ** CONSTS(48)/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))+(CONSTS(49)/ALGBRC(5))/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47))))
      ALGBRC(13) = ( CONSTS(53)*CONSTS(54))/(1.00000+ (CONSTS(55)/STATES(6))*(1.00000+STATES(9)/CONSTS(52))*(1.00000+STATES(7)/CONSTS(56)))
      ALGBRC(16) =  CONSTS(70)*( STATES(9)*CONSTS(10) - ( STATES(4)*STATES(10))/CONSTS(69))
      ALGBRC(17) =  CONSTS(68)*(STATES(7) - STATES(8)/CONSTS(67))
      ALGBRC(14) = ( CONSTS(58)*CONSTS(59)*CONSTS(133)*CONSTS(134))/(1.00000+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))+CONSTS(61)/ALGBRC(5)+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))*(CONSTS(61)/ALGBRC(5)))
      ALGBRC(8) = ( CONSTS(29)*CONSTS(30))/(1.00000+CONSTS(27)/CONSTS(15)+CONSTS(28)/STATES(9)+ (CONSTS(27)/CONSTS(15))*(CONSTS(28)/STATES(9)))
      ALGBRC(19) =  CONSTS(71)*STATES(10)
      ALGBRC(18) =  (CONSTS(89)+CONSTS(90)/(1.00000+CONSTS(93)/STATES(13)))*(CONSTS(88)+CONSTS(87)/(1.00000+EXP( CONSTS(91)*(CONSTS(92) - STATES(16)))))*STATES(16)
      ALGBRC(21) =  CONSTS(94)*(ALGBRC(18)/STATES(16))*(STATES(16) -  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*arbitrary_log(STATES(12)/STATES(13), 10))
      ALGBRC(15) = ( 2.00000*CONSTS(72)*CONSTS(74)*(CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))*CONSTS(75)*STATES(13))/( CONSTS(74)*( 2.00000*CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))+ STATES(13)*CONSTS(72)*CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))
      ALGBRC(20) =  2.00000*CONSTS(77)*CONSTS(78)*STATES(14)*EXP( - CONSTS(79)*STATES(14))
      ALGBRC(22) = ( CONSTS(82)*STATES(14)*STATES(15))/( CONSTS(80)*STATES(15)+ CONSTS(81)*STATES(14))
      ALGBRC(6) =  2.00000*(CONSTS(21) - STATES(15))
      ALGBRC(24) = ( CONSTS(83)*CONSTS(86))/(1.00000+CONSTS(84)/ALGBRC(6)+CONSTS(85)/CONSTS(19)+ (CONSTS(84)/ALGBRC(6))*(CONSTS(85)/CONSTS(19)))
      ALGBRC(23) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(104)*STATES(2)/ALGBRC(5) ** 0.500000)
      ALGBRC(25) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*CONSTS(109)+STATES(16)
      ALGBRC(27) =  0.500000*CONSTS(95)*((( (CONSTS(97)+ CONSTS(98)*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) -  CONSTS(97)*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))+ CONSTS(103)*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(4) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, CONSTS(9), 0.0500000)
      ALGBRC(33) =  CONSTS(121)*((1.00000 - ( 0.0500000*CONSTS(7)*0.450000*0.800000*STATES(1))/( 0.450000*ALGBRC(4)*0.0500000*ALGBRC(3)))/( (1.00000+ (( 0.0500000*CONSTS(7))/( 0.450000*ALGBRC(4)))*EXP(( - CONSTS(119)*CONSTS(25)*CONSTS(120))/( CONSTS(23)*CONSTS(24))))*(1.00000+( 0.450000*0.800000*STATES(1))/( 0.0500000*ALGBRC(3)))))
      ALGBRC(30) =  (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(118)*(ALGBRC(3)/( STATES(1)*CONSTS(13))))
      ALGBRC(31) =  - CONSTS(110)*(( ( 100.000*CONSTS(111)+ CONSTS(112)*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - ( CONSTS(111)*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))+ CONSTS(117)*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(1) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, 0.500000, 4.00000)
      ALGBRC(2) = TERNRY(VOI.GE.CONSTS(2).AND.VOI.LE.CONSTS(3).AND.(VOI - CONSTS(2)) -  INT((VOI - CONSTS(2))/ALGBRC(1))*ALGBRC(1).LE.CONSTS(4), CONSTS(5), 0.100000)
      ALGBRC(34) =  CONSTS(124)*(( (ALGBRC(2)/CONSTS(126))*1.00000+ALGBRC(2)/CONSTS(126) ** 3.00000*(( 2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))/(1.00000+ALGBRC(2)/CONSTS(126) ** 4.00000+ (CONSTS(122)/1.00000+ALGBRC(2)/CONSTS(125) ** CONSTS(123))*(1.00000 - EXP(( -2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))))
      ALGBRC(35) =  CONSTS(128)*(( EXP(( CONSTS(131)*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))*EXP(log(ALGBRC(2)/STATES(11))))/( 1.00000+CONSTS(129)/CONSTS(6) ** CONSTS(127)*(1.00000+CONSTS(130)/STATES(11))))
      ALGBRC(26) =  6.00000*CONSTS(95)*(( CONSTS(97)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))) -  (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(28) =  6.00000*CONSTS(96)*(( CONSTS(97)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))) -  (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(32) =  -3.00000*CONSTS(110)*(( 100.000*CONSTS(111)*(1.00000+EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24)))) -  (CONSTS(111)+CONSTS(116))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))
      ALGBRC(29) =  CONSTS(106)*ALGBRC(25)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END