C
C There are a total of 34 entries in the algebraic variable array.
C There are a total of 19 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 t in component Environment (second).
C CONSTS(1) is Mg_tot in component Environment (molar).
C CONSTS(2) is Pi_e in component Environment (molar).
C CONSTS(3) is ADP_e in component Environment (molar).
C CONSTS(4) is RT in component Fixed_parameters (kilojoule_per_mole).
C CONSTS(5) is F in component Fixed_parameters (kilojoule_per_mole_per_millivolt).
C CONSTS(6) is n_A in component Fixed_parameters (dimensionless).
C CONSTS(7) is dG_C1o in component Fixed_parameters (kilojoule_per_mole).
C CONSTS(8) is dG_C3o in component Fixed_parameters (kilojoule_per_mole).
C CONSTS(9) is dG_C4o in component Fixed_parameters (kilojoule_per_mole).
C CONSTS(10) is dG_F1o in component Fixed_parameters (kilojoule_per_mole).
C CONSTS(11) is pH_e in component Fixed_parameters (dimensionless).
C CONSTS(48) is H_e in component Fixed_parameters (molar).
C CONSTS(12) is K_e in component Fixed_parameters (molar).
C CONSTS(13) is ATP_e in component Fixed_parameters (molar).
C CONSTS(14) is AMP_e in component Fixed_parameters (molar).
C CONSTS(49) is k_dHPi in component Fixed_parameters (molar).
C CONSTS(50) is k_dHatp in component Fixed_parameters (molar).
C CONSTS(51) is k_dHadp in component Fixed_parameters (molar).
C CONSTS(15) is K_DT in component Fixed_parameters (molar).
C CONSTS(16) is K_DD in component Fixed_parameters (molar).
C CONSTS(17) is K_AK in component Fixed_parameters (dimensionless).
C CONSTS(18) is W_m in component Fixed_parameters (l_water_per_l_mito).
C CONSTS(53) is W_x in component Fixed_parameters (l_water_per_l_mito).
C CONSTS(56) is W_i in component Fixed_parameters (l_water_per_l_mito).
C CONSTS(19) is gamma in component Fixed_parameters (per_micron).
C CONSTS(20) is Ctot in component Fixed_parameters (molar).
C CONSTS(21) is Qtot in component Fixed_parameters (molar).
C CONSTS(22) is NADtot in component Fixed_parameters (molar).
C CONSTS(55) is H_i in component Fixed_parameters (molar).
C CONSTS(52) is K_i in component Fixed_parameters (molar).
C CONSTS(23) is k_Pi1 in component Adjustable_parameters (molar).
C CONSTS(24) is k_Pi2 in component Adjustable_parameters (molar).
C CONSTS(25) is k_Pi3 in component Adjustable_parameters (molar).
C CONSTS(26) is k_Pi4 in component Adjustable_parameters (molar).
C CONSTS(27) is k_PiH in component Adjustable_parameters (molar).
C CONSTS(28) is r in component Adjustable_parameters (dimensionless).
C CONSTS(29) is x_DH in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar).
C CONSTS(30) is x_C1 in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar).
C CONSTS(31) is x_C3 in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar).
C CONSTS(32) is x_C4 in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar).
C CONSTS(33) is x_F1 in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
C CONSTS(34) is x_ANT in component Adjustable_parameters (mole_per_second_per_l_mito).
C CONSTS(35) is x_Pi1 in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar).
C CONSTS(36) is x_KH in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
C CONSTS(37) is x_Hle in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
C CONSTS(38) is x_K in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_millivolt).
C CONSTS(39) is k_mADP in component Adjustable_parameters (molar).
C CONSTS(40) is x_AK in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
C CONSTS(41) is p_A in component Adjustable_parameters (micron_per_second).
C CONSTS(42) is k_O2 in component Adjustable_parameters (molar).
C CONSTS(43) is x_buff in component Adjustable_parameters (per_molar).
C CONSTS(44) is x_MgA in component Adjustable_parameters (mole_per_second_per_l_mito_per_molar_per_molar).
C CONSTS(45) is x_Pi2 in component Adjustable_parameters (micron_per_second).
C ALGBRC(1) is dG_H in component Proton_motive_force (kilojoule_per_mole).
C STATES(1) is dPsi in component dPsi_dt (millivolt).
C STATES(2) is H_x in component dH_x_dt (molar).
C ALGBRC(27) is J_DH in component Dehydrogenase_flux (mole_per_second_per_l_mito).
C ALGBRC(26) is NAD_x in component NAD_x_concentration (molar).
C STATES(3) is NADH_x in component dNADH_x_dt (molar).
C STATES(4) is Pi_x in component dPi_x_dt (molar).
C ALGBRC(30) is J_C1 in component Electron_flux_complex_I (mole_per_second_per_l_mito).
C ALGBRC(29) is dG_C1op in component Electron_flux_complex_I (kilojoule_per_mole).
C ALGBRC(28) is Q in component Q_concentration (molar).
C STATES(5) is QH2 in component dQH2_dt (molar).
C ALGBRC(33) is J_C3 in component Electron_flux_complex_III (mole_per_second_per_l_mito).
C ALGBRC(31) is dG_C3op in component Electron_flux_complex_III (kilojoule_per_mole).
C ALGBRC(32) is Cox in component Cox_concentration (molar).
C STATES(6) is Cred in component dCred_dt (molar).
C ALGBRC(34) is J_C4 in component Electron_flux_complex_IV (mole_per_second_per_l_mito).
C ALGBRC(2) is dG_C4op in component Electron_flux_complex_IV (kilojoule_per_mole).
C STATES(7) is O2 in component dO2_dt (molar).
C ALGBRC(3) is J_F1 in component ATP_synthesis_flux (mole_per_second_per_l_mito).
C STATES(8) is ADP_mx in component dADP_mx_dt (molar).
C STATES(9) is ATP_mx in component dATP_mx_dt (molar).
C ALGBRC(13) is J_ANT in component ANT_flux (mole_per_second_per_l_mito).
C ALGBRC(4) is Psi_x in component ANT_flux (millivolt).
C ALGBRC(5) is Psi_i in component ANT_flux (millivolt).
C ALGBRC(11) is ADP_fi in component MgADPi_binding_flux (molar).
C ALGBRC(9) is ATP_fi in component MgATPi_binding_flux (molar).
C ALGBRC(8) is ADP_fx in component MgADPx_binding_flux (molar).
C ALGBRC(6) is ATP_fx in component MgATPx_binding_flux (molar).
C CONSTS(46) is mincond in component ANT_flux (molar).
C ALGBRC(7) is J_MgATPx in component MgATPx_binding_flux (mole_per_second_per_l_mito).
C STATES(10) is ATP_x in component dATP_x_dt (molar).
C STATES(11) is Mg_x in component dMg_x_dt (molar).
C ALGBRC(10) is J_MgADPx in component MgADPx_binding_flux (mole_per_second_per_l_mito).
C STATES(12) is ADP_x in component dADP_x_dt (molar).
C ALGBRC(12) is J_MgATPi in component MgATPi_binding_flux (mole_per_second_per_l_mito).
C STATES(13) is ATP_i in component dATP_i_dt (molar).
C STATES(14) is ATP_mi in component dATP_mi_dt (molar).
C CONSTS(59) is Mg_i in component Mg_binding (molar).
C ALGBRC(14) is J_MgADPi in component MgADPi_binding_flux (mole_per_second_per_l_mito).
C STATES(15) is ADP_i in component dADP_i_dt (molar).
C STATES(16) is ADP_mi in component dADP_mi_dt (molar).
C ALGBRC(16) is J_ATP in component ATP_substrate_flux (mole_per_second_per_l_mito).
C ALGBRC(18) is J_ADP in component ADP_substrate_flux (mole_per_second_per_l_mito).
C ALGBRC(20) is J_AMP in component AMP_substrate_flux (mole_per_second_per_l_mito).
C STATES(17) is AMP_i in component dAMP_i_dt (molar).
C ALGBRC(15) is J_Pi2 in component Pi_substrate_flux (mole_per_second_per_l_mito).
C STATES(18) is Pi_i in component dPi_i_dt (molar).
C ALGBRC(21) is J_Pi1 in component Phosphate_hydrogen_cotransporter_flux (mole_per_second_per_l_mito).
C ALGBRC(17) is H2PIi in component Phosphate_hydrogen_cotransporter_flux (molar).
C ALGBRC(19) is H2PIx in component Phosphate_hydrogen_cotransporter_flux (molar).
C ALGBRC(22) is J_AKi in component Adenylate_kinase_flux (mole_per_second_per_l_mito).
C ALGBRC(23) is J_Hle in component Hydrogen_leak_flux (mole_per_second_per_l_mito).
C ALGBRC(24) is J_K in component Passive_potassium_flux (mole_per_second_per_l_mito).
C STATES(19) is K_x in component dK_x_dt (molar).
C ALGBRC(25) is J_KH in component Potassium_hydrogen_flux (mole_per_second_per_l_mito).
C CONSTS(54) is ADP_me in component ADP_binding (molar).
C CONSTS(57) is ADP_fe in component ADP_binding (molar).
C CONSTS(58) is Mg_e in component Mg_binding (molar).
C CONSTS(47) is C_im in component dPsi_dt (mole_per_l_mito_per_millivolt).
C RATES(2) is d/dt H_x in component dH_x_dt (molar).
C RATES(19) is d/dt K_x in component dK_x_dt (molar).
C RATES(11) is d/dt Mg_x in component dMg_x_dt (molar).
C RATES(3) is d/dt NADH_x in component dNADH_x_dt (molar).
C RATES(5) is d/dt QH2 in component dQH2_dt (molar).
C RATES(6) is d/dt Cred in component dCred_dt (molar).
C RATES(10) is d/dt ATP_x in component dATP_x_dt (molar).
C RATES(12) is d/dt ADP_x in component dADP_x_dt (molar).
C RATES(9) is d/dt ATP_mx in component dATP_mx_dt (molar).
C RATES(8) is d/dt ADP_mx in component dADP_mx_dt (molar).
C RATES(4) is d/dt Pi_x in component dPi_x_dt (molar).
C RATES(13) is d/dt ATP_i in component dATP_i_dt (molar).
C RATES(15) is d/dt ADP_i in component dADP_i_dt (molar).
C RATES(17) is d/dt AMP_i in component dAMP_i_dt (molar).
C RATES(14) is d/dt ATP_mi in component dATP_mi_dt (molar).
C RATES(16) is d/dt ADP_mi in component dADP_mi_dt (molar).
C RATES(18) is d/dt Pi_i in component dPi_i_dt (molar).
C RATES(1) is d/dt dPsi in component dPsi_dt (millivolt).
C RATES(7) is d/dt O2 in component dO2_dt (molar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.005
      CONSTS(2) = 0.000125
      CONSTS(3) = 0
      CONSTS(4) = 2.4734
      CONSTS(5) = 0.096484
      CONSTS(6) = 3
      CONSTS(7) = -69.37
      CONSTS(8) = -32.53
      CONSTS(9) = -122.94
      CONSTS(10) = 36.03
      CONSTS(11) = 7.1
      CONSTS(12) = 0.15
      CONSTS(13) = 0
      CONSTS(14) = 0
      CONSTS(15) = 2.4e-5
      CONSTS(16) = 3.47e-4
      CONSTS(17) = 0.4331
      CONSTS(18) = 0.72376
      CONSTS(19) = 5.99
      CONSTS(20) = 0.0027
      CONSTS(21) = 0.00135
      CONSTS(22) = 0.00297
      CONSTS(23) = 1.3413e-4
      CONSTS(24) = 6.7668e-4
      CONSTS(25) = 1.9172e-4
      CONSTS(26) = 0.02531
      CONSTS(27) = 4.5082e-4
      CONSTS(28) = 4.5807
      CONSTS(29) = 0.09183
      CONSTS(30) = 0.36923
      CONSTS(31) = 0.091737
      CONSTS(32) = 3.2562e-5
      CONSTS(33) = 150.93
      CONSTS(34) = 0.0079204
      CONSTS(35) = 339430
      CONSTS(36) = 2.9802e7
      CONSTS(37) = 250
      CONSTS(38) = 0
      CONSTS(39) = 3.5e-6
      CONSTS(40) = 0
      CONSTS(41) = 85
      CONSTS(42) = 1.2e-4
      CONSTS(43) = 100
      CONSTS(44) = 1000000
      CONSTS(45) = 327
      STATES(1) = 160
      STATES(2) = 6.30957344480193e-8
      STATES(3) = 0.0015
      STATES(4) = 0.001
      STATES(5) = 8e-4
      STATES(6) = 0.001
      STATES(7) = 2.6e-5
      STATES(8) = 0
      STATES(9) = 0
      CONSTS(46) = 1e-12
      STATES(10) = 0
      STATES(11) = 0.005
      STATES(12) = 0.01
      STATES(13) = 0
      STATES(14) = 0
      STATES(15) = 0
      STATES(16) = 0
      STATES(17) = 0
      STATES(18) = 0.001
      STATES(19) = 0.14
      CONSTS(47) = 6.756756756756757e-6
      CONSTS(48) =  1.00000*10.0000 ** - CONSTS(11)
      CONSTS(49) =  1.00000*10.0000 ** - 6.75000
      CONSTS(50) =  1.00000*10.0000 ** - 6.48000
      CONSTS(51) =  1.00000*10.0000 ** - 6.29000
      CONSTS(52) = CONSTS(12)
      CONSTS(53) =  0.900000*CONSTS(18)
      CONSTS(54) = ((CONSTS(16)+CONSTS(3)+CONSTS(1)) -  (CONSTS(16)+CONSTS(3)+CONSTS(1) ** 2.00000 -  4.00000*CONSTS(1)*CONSTS(3)) ** (1.0 / 2))/2.00000
      CONSTS(59) = 0.00000
      CONSTS(55) = CONSTS(48)
      CONSTS(56) =  0.100000*CONSTS(18)
      CONSTS(57) = CONSTS(3) - CONSTS(54)
      CONSTS(58) = CONSTS(1) - CONSTS(54)
      CONSTS(59) = CONSTS(58)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(7) = CONSTS(59)
      ALGBRC(6) = STATES(10) - STATES(9)
      ALGBRC(7) =  CONSTS(44)*( ALGBRC(6)*STATES(11) -  CONSTS(15)*STATES(9))
      RATES(9) = ALGBRC(7)/CONSTS(53)
      ALGBRC(8) = STATES(12) - STATES(8)
      ALGBRC(10) =  CONSTS(44)*( ALGBRC(8)*STATES(11) -  CONSTS(16)*STATES(8))
      RATES(11) = (- ALGBRC(7) - ALGBRC(10))/CONSTS(53)
      RATES(8) = ALGBRC(10)/CONSTS(53)
      ALGBRC(9) = STATES(13) - STATES(14)
      ALGBRC(12) =  CONSTS(44)*( ALGBRC(9)*CONSTS(59) -  CONSTS(15)*STATES(14))
      RATES(14) = ALGBRC(12)/CONSTS(56)
      ALGBRC(1) =  CONSTS(5)*STATES(1)+ CONSTS(4)*log(CONSTS(55)/STATES(2))
      ALGBRC(3) =  CONSTS(33)*( (( EXP(- (CONSTS(10) -  CONSTS(6)*ALGBRC(1))/CONSTS(4))*CONSTS(16))/CONSTS(15))*STATES(8)*STATES(4) -  STATES(9)*1.00000)
      ALGBRC(4) =  - 0.650000*STATES(1)
      ALGBRC(5) =  0.350000*STATES(1)
      ALGBRC(11) = STATES(15) - STATES(16)
      ALGBRC(13) = TERNRY(ALGBRC(11).GT.CONSTS(46).OR.ALGBRC(9).GT.CONSTS(46), ( CONSTS(34)*(ALGBRC(11)/(ALGBRC(11)+ ALGBRC(9)*EXP(( - CONSTS(5)*ALGBRC(5))/CONSTS(4))) - ALGBRC(8)/(ALGBRC(8)+ ALGBRC(6)*EXP(( - CONSTS(5)*ALGBRC(4))/CONSTS(4))))*ALGBRC(11))/(CONSTS(39)+ALGBRC(11)), 0.00000)
      RATES(10) = (ALGBRC(3) - ALGBRC(13))/CONSTS(53)
      RATES(12) = (- ALGBRC(3)+ALGBRC(13))/CONSTS(53)
      ALGBRC(14) =  CONSTS(44)*( ALGBRC(11)*CONSTS(59) -  CONSTS(16)*STATES(16))
      RATES(16) = ALGBRC(14)/CONSTS(56)
      ALGBRC(17) = ( STATES(18)*CONSTS(55))/(CONSTS(55)+CONSTS(49))
      ALGBRC(19) = ( STATES(4)*STATES(2))/(STATES(2)+CONSTS(49))
      ALGBRC(21) = ( CONSTS(35)*( STATES(2)*ALGBRC(17) -  CONSTS(55)*ALGBRC(19)))/(ALGBRC(17)+CONSTS(27))
      RATES(4) = (- ALGBRC(3)+ALGBRC(21))/CONSTS(53)
      ALGBRC(16) =  CONSTS(19)*CONSTS(41)*(CONSTS(13) - STATES(13))
      ALGBRC(22) =  CONSTS(40)*( CONSTS(17)*STATES(15)*STATES(15) -  STATES(17)*STATES(13))
      RATES(13) = (ALGBRC(16)+ALGBRC(13)+ALGBRC(22))/CONSTS(56)
      ALGBRC(18) =  CONSTS(19)*CONSTS(41)*(CONSTS(3) - STATES(15))
      RATES(15) = ((ALGBRC(18) - ALGBRC(13)) -  2.00000*ALGBRC(22))/CONSTS(56)
      ALGBRC(20) =  CONSTS(19)*CONSTS(41)*(CONSTS(14) - STATES(17))
      RATES(17) = (ALGBRC(20)+ALGBRC(22))/CONSTS(56)
      ALGBRC(15) =  CONSTS(19)*CONSTS(45)*(CONSTS(2) - STATES(18))
      RATES(18) = (- ALGBRC(21)+ALGBRC(15))/CONSTS(56)
      ALGBRC(24) = ( CONSTS(38)*STATES(1)*( CONSTS(52)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - STATES(19)))/(EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - 1.00000)
      ALGBRC(25) =  CONSTS(36)*( CONSTS(52)*STATES(2) -  STATES(19)*CONSTS(55))
      RATES(19) = (ALGBRC(25)+ALGBRC(24))/CONSTS(53)
      ALGBRC(26) = CONSTS(22) - STATES(3)
      ALGBRC(27) = ( CONSTS(29)*( CONSTS(28)*ALGBRC(26) - STATES(3))*(1.00000+STATES(4)/CONSTS(23)))/(1.00000+STATES(4)/CONSTS(24))
      ALGBRC(28) = CONSTS(21) - STATES(5)
      ALGBRC(29) = (CONSTS(7) -  CONSTS(4)*log(STATES(2)/1.00000e-07)) -  CONSTS(4)*log(ALGBRC(28)/STATES(5))
      ALGBRC(30) =  CONSTS(30)*( EXP(- (ALGBRC(29)+ 4.00000*ALGBRC(1))/CONSTS(4))*STATES(3) - ALGBRC(26))
      RATES(3) = (ALGBRC(27) - ALGBRC(30))/CONSTS(53)
      ALGBRC(31) = (CONSTS(8)+ 2.00000*CONSTS(4)*log(STATES(2)/1.00000e-07)) -  CONSTS(4)*log(STATES(5)/ALGBRC(28))
      ALGBRC(32) = CONSTS(20) - STATES(6)
      ALGBRC(33) =  (( CONSTS(31)*(1.00000+STATES(4)/CONSTS(25)))/(1.00000+STATES(4)/CONSTS(26)))*( EXP(- ((ALGBRC(31)+ 4.00000*ALGBRC(1)) -  2.00000*CONSTS(5)*STATES(1))/( 2.00000*CONSTS(4)))*ALGBRC(32) - STATES(6))
      RATES(5) = (ALGBRC(30) - ALGBRC(33))/CONSTS(53)
      ALGBRC(2) = (CONSTS(9) -  2.00000*CONSTS(4)*log(STATES(2)/1.00000e-07)) -  (CONSTS(4)/2.00000)*log(STATES(7)/1.00000)
      ALGBRC(34) =  (( (( CONSTS(32)*1.00000)/(1.00000+CONSTS(42)/STATES(7)))*STATES(6))/CONSTS(20))*( EXP(- (ALGBRC(2)+ 2.00000*ALGBRC(1))/( 2.00000*CONSTS(4)))*STATES(6) -  ALGBRC(32)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)))
      ALGBRC(23) = ( CONSTS(37)*STATES(1)*( CONSTS(55)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - STATES(2)))/(EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - 1.00000)
      RATES(2) = ( CONSTS(43)*STATES(2)*(((((ALGBRC(27) -  5.00000*ALGBRC(30)) -  2.00000*ALGBRC(33)) -  4.00000*ALGBRC(34))+ (CONSTS(6) - 1.00000)*ALGBRC(3)+ 2.00000*ALGBRC(21)+ALGBRC(23)) - ALGBRC(25)))/CONSTS(53)
      RATES(6) = ( 2.00000*ALGBRC(33) -  2.00000*ALGBRC(34))/CONSTS(56)
      RATES(1) = ((((( 4.00000*ALGBRC(30)+ 2.00000*ALGBRC(33)+ 4.00000*ALGBRC(34)) -  CONSTS(6)*ALGBRC(3)) - ALGBRC(13)) - ALGBRC(23)) - ALGBRC(24))/CONSTS(47)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(6) = STATES(10) - STATES(9)
      ALGBRC(7) =  CONSTS(44)*( ALGBRC(6)*STATES(11) -  CONSTS(15)*STATES(9))
      ALGBRC(8) = STATES(12) - STATES(8)
      ALGBRC(10) =  CONSTS(44)*( ALGBRC(8)*STATES(11) -  CONSTS(16)*STATES(8))
      ALGBRC(9) = STATES(13) - STATES(14)
      ALGBRC(12) =  CONSTS(44)*( ALGBRC(9)*CONSTS(59) -  CONSTS(15)*STATES(14))
      ALGBRC(1) =  CONSTS(5)*STATES(1)+ CONSTS(4)*log(CONSTS(55)/STATES(2))
      ALGBRC(3) =  CONSTS(33)*( (( EXP(- (CONSTS(10) -  CONSTS(6)*ALGBRC(1))/CONSTS(4))*CONSTS(16))/CONSTS(15))*STATES(8)*STATES(4) -  STATES(9)*1.00000)
      ALGBRC(4) =  - 0.650000*STATES(1)
      ALGBRC(5) =  0.350000*STATES(1)
      ALGBRC(11) = STATES(15) - STATES(16)
      ALGBRC(13) = TERNRY(ALGBRC(11).GT.CONSTS(46).OR.ALGBRC(9).GT.CONSTS(46), ( CONSTS(34)*(ALGBRC(11)/(ALGBRC(11)+ ALGBRC(9)*EXP(( - CONSTS(5)*ALGBRC(5))/CONSTS(4))) - ALGBRC(8)/(ALGBRC(8)+ ALGBRC(6)*EXP(( - CONSTS(5)*ALGBRC(4))/CONSTS(4))))*ALGBRC(11))/(CONSTS(39)+ALGBRC(11)), 0.00000)
      ALGBRC(14) =  CONSTS(44)*( ALGBRC(11)*CONSTS(59) -  CONSTS(16)*STATES(16))
      ALGBRC(17) = ( STATES(18)*CONSTS(55))/(CONSTS(55)+CONSTS(49))
      ALGBRC(19) = ( STATES(4)*STATES(2))/(STATES(2)+CONSTS(49))
      ALGBRC(21) = ( CONSTS(35)*( STATES(2)*ALGBRC(17) -  CONSTS(55)*ALGBRC(19)))/(ALGBRC(17)+CONSTS(27))
      ALGBRC(16) =  CONSTS(19)*CONSTS(41)*(CONSTS(13) - STATES(13))
      ALGBRC(22) =  CONSTS(40)*( CONSTS(17)*STATES(15)*STATES(15) -  STATES(17)*STATES(13))
      ALGBRC(18) =  CONSTS(19)*CONSTS(41)*(CONSTS(3) - STATES(15))
      ALGBRC(20) =  CONSTS(19)*CONSTS(41)*(CONSTS(14) - STATES(17))
      ALGBRC(15) =  CONSTS(19)*CONSTS(45)*(CONSTS(2) - STATES(18))
      ALGBRC(24) = ( CONSTS(38)*STATES(1)*( CONSTS(52)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - STATES(19)))/(EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - 1.00000)
      ALGBRC(25) =  CONSTS(36)*( CONSTS(52)*STATES(2) -  STATES(19)*CONSTS(55))
      ALGBRC(26) = CONSTS(22) - STATES(3)
      ALGBRC(27) = ( CONSTS(29)*( CONSTS(28)*ALGBRC(26) - STATES(3))*(1.00000+STATES(4)/CONSTS(23)))/(1.00000+STATES(4)/CONSTS(24))
      ALGBRC(28) = CONSTS(21) - STATES(5)
      ALGBRC(29) = (CONSTS(7) -  CONSTS(4)*log(STATES(2)/1.00000e-07)) -  CONSTS(4)*log(ALGBRC(28)/STATES(5))
      ALGBRC(30) =  CONSTS(30)*( EXP(- (ALGBRC(29)+ 4.00000*ALGBRC(1))/CONSTS(4))*STATES(3) - ALGBRC(26))
      ALGBRC(31) = (CONSTS(8)+ 2.00000*CONSTS(4)*log(STATES(2)/1.00000e-07)) -  CONSTS(4)*log(STATES(5)/ALGBRC(28))
      ALGBRC(32) = CONSTS(20) - STATES(6)
      ALGBRC(33) =  (( CONSTS(31)*(1.00000+STATES(4)/CONSTS(25)))/(1.00000+STATES(4)/CONSTS(26)))*( EXP(- ((ALGBRC(31)+ 4.00000*ALGBRC(1)) -  2.00000*CONSTS(5)*STATES(1))/( 2.00000*CONSTS(4)))*ALGBRC(32) - STATES(6))
      ALGBRC(2) = (CONSTS(9) -  2.00000*CONSTS(4)*log(STATES(2)/1.00000e-07)) -  (CONSTS(4)/2.00000)*log(STATES(7)/1.00000)
      ALGBRC(34) =  (( (( CONSTS(32)*1.00000)/(1.00000+CONSTS(42)/STATES(7)))*STATES(6))/CONSTS(20))*( EXP(- (ALGBRC(2)+ 2.00000*ALGBRC(1))/( 2.00000*CONSTS(4)))*STATES(6) -  ALGBRC(32)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)))
      ALGBRC(23) = ( CONSTS(37)*STATES(1)*( CONSTS(55)*EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - STATES(2)))/(EXP(( CONSTS(5)*STATES(1))/CONSTS(4)) - 1.00000)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END