C
C There are a total of 117 entries in the algebraic variable array.
C There are a total of 48 entries in each of the rate and state variable arrays.
C There are a total of 95 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C CONSTS(1) is F in component model_parameters (C_per_mole).
C CONSTS(2) is T in component model_parameters (kelvin).
C CONSTS(3) is R in component model_parameters (J_per_kmole_K).
C CONSTS(4) is K_o in component model_parameters (mM).
C CONSTS(5) is Ca_o in component model_parameters (mM).
C CONSTS(6) is Na_o in component model_parameters (mM).
C CONSTS(7) is Cl_o in component model_parameters (mM).
C CONSTS(90) is Vmyo in component cell_geometry (uL).
C CONSTS(91) is Vnsr in component cell_geometry (uL).
C CONSTS(86) is Vsr in component cell_geometry (uL).
C CONSTS(92) is Vjsr in component cell_geometry (uL).
C CONSTS(93) is Vss_sr in component cell_geometry (uL).
C CONSTS(94) is Vss_CaL in component cell_geometry (uL).
C CONSTS(88) is Acap in component cell_geometry (cm2).
C CONSTS(89) is AF in component cell_geometry (uF_mole_per_C).
C CONSTS(8) is l in component cell_geometry (cm).
C CONSTS(9) is a in component cell_geometry (cm).
C CONSTS(83) is vcell in component cell_geometry (uL).
C CONSTS(87) is ageo in component cell_geometry (cm2).
C STATES(1) is Vm in component membrane (mV).
C ALGBRC(14) is i_Stim in component membrane (uA_per_uF).
C CONSTS(10) is stim_offset in component membrane (ms).
C CONSTS(11) is stim_period in component membrane (ms).
C CONSTS(12) is stim_duration in component membrane (ms).
C CONSTS(13) is stim_amplitude in component membrane (uA_per_uF).
C ALGBRC(1) is past in component membrane (ms).
C ALGBRC(71) is caiont in component membrane (uA_per_uF).
C ALGBRC(88) is naiont in component membrane (uA_per_uF).
C ALGBRC(112) is kiont in component membrane (uA_per_uF).
C ALGBRC(116) is clont in component membrane (uA_per_uF).
C ALGBRC(81) is INa in component INa (uA_per_uF).
C ALGBRC(38) is ICaL in component ICaL (uA_per_uF).
C ALGBRC(104) is IK1 in component IK1 (uA_per_uF).
C ALGBRC(94) is IKp in component IKp (uA_per_uF).
C ALGBRC(98) is IKs in component IKs (uA_per_uF).
C ALGBRC(109) is IKr in component IKr (uA_per_uF).
C ALGBRC(69) is IpCa in component IpCa (uA_per_uF).
C ALGBRC(70) is ICab in component ICab (uA_per_uF).
C ALGBRC(66) is INaCa in component INaCa (uA_per_uF).
C ALGBRC(67) is INaCa_ss_sr in component INaCa (uA_per_uF).
C ALGBRC(49) is INaK in component INaK (uA_per_uF).
C ALGBRC(111) is Ito1 in component Ito1 (uA_per_uF).
C ALGBRC(115) is Ito2 in component Ito2 (uA_per_uF).
C ALGBRC(100) is IClb in component IClb (uA_per_uF).
C ALGBRC(85) is INaL in component INaL (uA_per_uF).
C ALGBRC(75) is INab in component INab (uA_per_uF).
C ALGBRC(26) is gKs in component IKs (mS_per_uF).
C STATES(2) is C1 in component IKs (dimensionless).
C STATES(3) is C2 in component IKs (dimensionless).
C STATES(4) is C3 in component IKs (dimensionless).
C STATES(5) is C4 in component IKs (dimensionless).
C STATES(6) is C5 in component IKs (dimensionless).
C STATES(7) is C6 in component IKs (dimensionless).
C STATES(8) is C7 in component IKs (dimensionless).
C STATES(9) is C8 in component IKs (dimensionless).
C STATES(10) is C9 in component IKs (dimensionless).
C STATES(11) is C10 in component IKs (dimensionless).
C STATES(12) is C11 in component IKs (dimensionless).
C STATES(13) is C12 in component IKs (dimensionless).
C STATES(14) is C13 in component IKs (dimensionless).
C STATES(15) is C14 in component IKs (dimensionless).
C STATES(16) is C15 in component IKs (dimensionless).
C STATES(17) is O1 in component IKs (dimensionless).
C STATES(18) is O2 in component IKs (dimensionless).
C ALGBRC(31) is OKs in component IKs (dimensionless).
C ALGBRC(2) is alpha in component IKs (per_ms).
C ALGBRC(15) is beta in component IKs (per_ms).
C ALGBRC(32) is delta in component IKs (per_ms).
C ALGBRC(27) is gamma in component IKs (per_ms).
C ALGBRC(36) is eta in component IKs (per_ms).
C CONSTS(82) is theta in component IKs (per_ms).
C ALGBRC(39) is omega in component IKs (per_ms).
C ALGBRC(41) is psi in component IKs (per_ms).
C ALGBRC(96) is EKs in component reversal_potentials (mV).
C STATES(19) is Ca_i in component Ca (mM).
C ALGBRC(35) is ICaL_max in component ICaL (uA_per_uF).
C CONSTS(14) is PCa in component ICaL (L_per_F_ms).
C STATES(20) is C in component ICaL (dimensionless).
C STATES(21) is O in component ICaL (dimensionless).
C STATES(22) is C_star in component ICaL (dimensionless).
C STATES(23) is O_star in component ICaL (dimensionless).
C STATES(24) is CI in component ICaL (dimensionless).
C STATES(25) is OI in component ICaL (dimensionless).
C STATES(26) is CI_star in component ICaL (dimensionless).
C STATES(27) is OI_star in component ICaL (dimensionless).
C ALGBRC(3) is ACT_tau in component ICaL (ms).
C ALGBRC(40) is IV_tau in component ICaL (ms).
C ALGBRC(52) is IV_tau_star in component ICaL (ms).
C ALGBRC(16) is ACT_infinity in component ICaL (dimensionless).
C ALGBRC(37) is IV_infinity in component ICaL (dimensionless).
C ALGBRC(46) is IV_infinity_star in component ICaL (dimensionless).
C ALGBRC(50) is IV_beta_star in component ICaL (dimensionless).
C ALGBRC(48) is IV_beta_infinity_star in component ICaL (dimensionless).
C ALGBRC(28) is alpha in component ICaL (per_ms).
C ALGBRC(33) is beta in component ICaL (per_ms).
C ALGBRC(58) is delta in component ICaL (per_ms).
C ALGBRC(42) is x in component ICaL (per_ms).
C ALGBRC(44) is y in component ICaL (per_ms).
C CONSTS(15) is theta in component ICaL (per_ms).
C ALGBRC(60) is delta_I in component ICaL (per_ms).
C ALGBRC(54) is x_star in component ICaL (per_ms).
C ALGBRC(56) is y_star in component ICaL (per_ms).
C CONSTS(16) is theta_I in component ICaL (per_ms).
C CONSTS(17) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(18) is gamma_Cao in component ICaL (dimensionless).
C STATES(28) is Ca_ss_CaL in component Ca (mM).
C CONSTS(19) is g_Na in component INa (mS_per_uF).
C ALGBRC(78) is ENa in component reversal_potentials (mV).
C STATES(29) is m in component INa_m_gate (dimensionless).
C STATES(30) is h in component INa_h_gate (dimensionless).
C STATES(31) is j in component INa_j_gate (dimensionless).
C ALGBRC(4) is am in component INa_m_gate (per_ms).
C ALGBRC(17) is bm in component INa_m_gate (per_ms).
C ALGBRC(5) is ah in component INa_h_gate (per_ms).
C ALGBRC(18) is bh in component INa_h_gate (per_ms).
C ALGBRC(6) is aj in component INa_j_gate (per_ms).
C ALGBRC(19) is bj in component INa_j_gate (per_ms).
C CONSTS(20) is kmnai in component INaK (mM).
C CONSTS(21) is kmko in component INaK (mM).
C CONSTS(22) is ibarnak in component INaK (uA_per_uF).
C CONSTS(23) is sigma in component INaK (dimensionless).
C CONSTS(24) is V_half in component INaK (mV).
C ALGBRC(43) is phi in component INaK (dimensionless).
C ALGBRC(45) is fv in component INaK (dimensionless).
C CONSTS(95) is PK in component INaK (dimensionless).
C ALGBRC(47) is PNa in component INaK (dimensionless).
C STATES(32) is Na_i in component Na (mM).
C CONSTS(25) is KmCa_act in component INaCa (mM).
C CONSTS(26) is Vmax in component INaCa (uA_per_uF).
C CONSTS(27) is ksat in component INaCa (dimensionless).
C CONSTS(28) is eta in component INaCa (dimensionless).
C CONSTS(29) is KmNai in component INaCa (mM).
C CONSTS(30) is KmNao in component INaCa (mM).
C CONSTS(31) is KmCai in component INaCa (mM).
C CONSTS(32) is KmCao in component INaCa (mM).
C STATES(33) is Ca_ss_sr in component Ca (mM).
C STATES(34) is Na_ss_sr in component Na (mM).
C ALGBRC(51) is num in component INaCa (mM4_uA_per_uF).
C ALGBRC(61) is num_ss in component INaCa (mM4_uA_per_uF).
C ALGBRC(53) is denom_1 in component INaCa (dimensionless).
C ALGBRC(55) is denom_2 in component INaCa (dimensionless).
C ALGBRC(57) is denom_3 in component INaCa (mM4).
C ALGBRC(59) is denom_4 in component INaCa (mM4).
C ALGBRC(62) is denom_ss_1 in component INaCa (dimensionless).
C ALGBRC(63) is denom_ss_2 in component INaCa (dimensionless).
C ALGBRC(64) is denom_ss_3 in component INaCa (mM4).
C ALGBRC(65) is denom_ss_4 in component INaCa (mM4).
C CONSTS(33) is gKp in component IKp (mS_per_uF).
C ALGBRC(68) is Kp in component IKp (dimensionless).
C ALGBRC(91) is EK in component reversal_potentials (mV).
C CONSTS(34) is gpCa in component IpCa (uA_per_uF).
C CONSTS(35) is KmpCa in component IpCa (mM).
C CONSTS(36) is PCab in component ICab (L_per_F_ms).
C CONSTS(37) is gamma_Ca_i in component ICab (dimensionless).
C CONSTS(38) is gamma_Ca_o in component ICab (dimensionless).
C CONSTS(39) is PNab in component INab (L_per_F_ms).
C ALGBRC(73) is phi in component INab (dimensionless).
C CONSTS(40) is gClb in component IClb (mS_per_uF).
C ALGBRC(99) is ECl in component reversal_potentials (mV).
C CONSTS(41) is gNaL in component INaL (mS_per_uF).
C STATES(35) is mL in component INaL_mL_gate (dimensionless).
C STATES(36) is hL in component INaL_hL_gate (dimensionless).
C ALGBRC(7) is amL in component INaL_mL_gate (per_ms).
C ALGBRC(20) is bmL in component INaL_mL_gate (per_ms).
C CONSTS(42) is tau_hL in component INaL_hL_gate (ms).
C ALGBRC(8) is hL_infinity in component INaL_hL_gate (dimensionless).
C CONSTS(43) is prnak in component reversal_potentials (dimensionless).
C STATES(37) is K_i in component K (mM).
C STATES(38) is Cl_i in component Cl (mM).
C CONSTS(84) is gK1 in component IK1 (mS_per_uF).
C ALGBRC(103) is K1 in component IK1_K1_gate (dimensionless).
C ALGBRC(101) is alpha_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(102) is beta_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(105) is CT_Na_Cl in component CT_Na_Cl (mM_per_ms).
C CONSTS(44) is CT_Na_Cl_max in component CT_Na_Cl (mM_per_ms).
C ALGBRC(106) is CT_K_Cl in component CT_K_Cl (mM_per_ms).
C CONSTS(45) is CT_K_Cl_max in component CT_K_Cl (mM_per_ms).
C CONSTS(85) is gKr in component IKr (mS_per_uF).
C ALGBRC(108) is r in component IKr_r_gate (dimensionless).
C STATES(39) is xr in component IKr_xr_gate (dimensionless).
C ALGBRC(9) is xr_infinity in component IKr_xr_gate (dimensionless).
C ALGBRC(21) is tau_xr in component IKr_xr_gate (ms).
C CONSTS(46) is gto1 in component Ito1 (mS_per_uF).
C ALGBRC(110) is rto1 in component Ito1 (dimensionless).
C STATES(40) is a in component Ito1_a_gate (dimensionless).
C STATES(41) is i1f in component Ito1_i1f_gate (dimensionless).
C STATES(42) is i1s in component Ito1_i1s_gate (dimensionless).
C ALGBRC(10) is alpha_a in component Ito1_a_gate (per_ms).
C ALGBRC(22) is beta_a in component Ito1_a_gate (per_ms).
C ALGBRC(29) is tau_a in component Ito1_a_gate (ms).
C ALGBRC(34) is a_infinity in component Ito1_a_gate (dimensionless).
C ALGBRC(23) is alpha_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(11) is beta_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(24) is alpha_i1s in component Ito1_i1s_gate (per_ms).
C ALGBRC(12) is beta_i1s in component Ito1_i1s_gate (per_ms).
C CONSTS(47) is PCl in component Ito2 (L_per_F_ms).
C CONSTS(48) is zCl in component Ito2 (dimensionless).
C ALGBRC(113) is Ito2_max in component Ito2 (uA_per_uF).
C ALGBRC(114) is KCa_ito2 in component Ito2 (dimensionless).
C CONSTS(49) is kCa_ito2 in component Ito2 (mM_per_ms).
C STATES(43) is Irel in component Irel (mM_per_ms).
C STATES(44) is i2f in component Ito2_i2f_gate (dimensionless).
C ALGBRC(25) is alpha_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(13) is beta_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(30) is i2f_infinity in component Ito2_i2f_gate (dimensionless).
C CONSTS(50) is tau_i2f in component Ito2_i2f_gate (ms).
C ALGBRC(97) is rel_infinity in component Irel (mM_per_ms).
C ALGBRC(95) is tau_rel in component Irel (ms).
C CONSTS(51) is Krel_infinity in component Irel (mM).
C ALGBRC(92) is alpha_rel in component Irel (mM_per_uA_per_uF_per_ms).
C ALGBRC(89) is beta_tau in component Irel (ms).
C CONSTS(52) is beta_0 in component Irel (ms).
C CONSTS(53) is kappa in component Irel (mM_per_uA_per_uF_per_ms_per_ms).
C CONSTS(54) is h_rel in component Irel (dimensionless).
C CONSTS(55) is h_beta in component Irel (dimensionless).
C ALGBRC(86) is delta_beta_CaMK in component Irel (dimensionless).
C CONSTS(56) is delta_beta_0 in component Irel (dimensionless).
C CONSTS(57) is K_beta in component Irel (mM).
C CONSTS(58) is Krel_tau in component Irel (mM).
C STATES(45) is Ca_JSR in component Ca (mM).
C ALGBRC(82) is CaMK_active in component CaMK_active (dimensionless).
C ALGBRC(93) is Iup in component Iup (mM_per_ms).
C CONSTS(59) is kmup in component Iup (mM).
C ALGBRC(90) is delta_kmPLB in component Iup (mM).
C CONSTS(60) is delta_kmPLB_bar in component Iup (mM).
C CONSTS(61) is iupbar in component Iup (mM_per_ms).
C ALGBRC(87) is delta_iupCaMK in component Iup (dimensionless).
C CONSTS(62) is delta_iupCaMK_bar in component Iup (dimensionless).
C CONSTS(63) is kmCaMK in component Iup (mM).
C ALGBRC(72) is Ileak in component Ileak (mM_per_ms).
C CONSTS(64) is nsrbar in component Ileak (mM).
C STATES(46) is Ca_NSR in component Ca (mM).
C ALGBRC(74) is Itr in component Itr (mM_per_ms).
C CONSTS(65) is tautr in component Itr (ms).
C CONSTS(66) is tau_diff in component Ca (ms).
C CONSTS(67) is tau_diff_ss in component Ca (ms).
C ALGBRC(76) is Idiff in component Ca (mM_per_ms).
C ALGBRC(80) is Idiff_ss in component Ca (mM_per_ms).
C ALGBRC(79) is bmyo in component Ca (dimensionless).
C ALGBRC(83) is bss in component Ca (dimensionless).
C ALGBRC(84) is bss_sr in component Ca (dimensionless).
C CONSTS(68) is cmdn_bar in component Ca (mM).
C CONSTS(69) is trpn_bar in component Ca (mM).
C CONSTS(70) is km_cmdn in component Ca (mM).
C CONSTS(71) is km_trpn in component Ca (mM).
C CONSTS(72) is kmcsqn in component Ca (mM).
C CONSTS(73) is csqnbar in component Ca (mM).
C ALGBRC(77) is bcsqn in component Ca (dimensionless).
C CONSTS(74) is BSRmax in component Ca (mM).
C CONSTS(75) is KmBSR in component Ca (mM).
C CONSTS(76) is BSLmax in component Ca (mM).
C CONSTS(77) is KmBSL in component Ca (mM).
C ALGBRC(107) is Idiff_Na in component Na (mM_per_ms).
C STATES(47) is Cl_ss in component Cl (mM).
C ALGBRC(117) is Idiff_Cl in component Cl (mM_per_ms).
C CONSTS(78) is Km in component CaMK_active (mM).
C STATES(48) is CaMK_trap in component CaMK_active (dimensionless).
C CONSTS(79) is alpha_CaMK in component CaMK_active (per_ms).
C CONSTS(80) is beta_CaMK in component CaMK_active (per_ms).
C CONSTS(81) is CaMK_0 in component CaMK_active (dimensionless).
C RATES(1) is d/dt Vm in component membrane (mV).
C RATES(2) is d/dt C1 in component IKs (dimensionless).
C RATES(3) is d/dt C2 in component IKs (dimensionless).
C RATES(4) is d/dt C3 in component IKs (dimensionless).
C RATES(5) is d/dt C4 in component IKs (dimensionless).
C RATES(6) is d/dt C5 in component IKs (dimensionless).
C RATES(7) is d/dt C6 in component IKs (dimensionless).
C RATES(8) is d/dt C7 in component IKs (dimensionless).
C RATES(9) is d/dt C8 in component IKs (dimensionless).
C RATES(10) is d/dt C9 in component IKs (dimensionless).
C RATES(11) is d/dt C10 in component IKs (dimensionless).
C RATES(12) is d/dt C11 in component IKs (dimensionless).
C RATES(13) is d/dt C12 in component IKs (dimensionless).
C RATES(14) is d/dt C13 in component IKs (dimensionless).
C RATES(15) is d/dt C14 in component IKs (dimensionless).
C RATES(16) is d/dt C15 in component IKs (dimensionless).
C RATES(17) is d/dt O1 in component IKs (dimensionless).
C RATES(18) is d/dt O2 in component IKs (dimensionless).
C RATES(20) is d/dt C in component ICaL (dimensionless).
C RATES(21) is d/dt O in component ICaL (dimensionless).
C RATES(22) is d/dt C_star in component ICaL (dimensionless).
C RATES(23) is d/dt O_star in component ICaL (dimensionless).
C RATES(24) is d/dt CI in component ICaL (dimensionless).
C RATES(25) is d/dt OI in component ICaL (dimensionless).
C RATES(26) is d/dt CI_star in component ICaL (dimensionless).
C RATES(27) is d/dt OI_star in component ICaL (dimensionless).
C RATES(29) is d/dt m in component INa_m_gate (dimensionless).
C RATES(30) is d/dt h in component INa_h_gate (dimensionless).
C RATES(31) is d/dt j in component INa_j_gate (dimensionless).
C RATES(35) is d/dt mL in component INaL_mL_gate (dimensionless).
C RATES(36) is d/dt hL in component INaL_hL_gate (dimensionless).
C RATES(39) is d/dt xr in component IKr_xr_gate (dimensionless).
C RATES(40) is d/dt a in component Ito1_a_gate (dimensionless).
C RATES(41) is d/dt i1f in component Ito1_i1f_gate (dimensionless).
C RATES(42) is d/dt i1s in component Ito1_i1s_gate (dimensionless).
C RATES(44) is d/dt i2f in component Ito2_i2f_gate (dimensionless).
C RATES(43) is d/dt Irel in component Irel (mM_per_ms).
C RATES(19) is d/dt Ca_i in component Ca (mM).
C RATES(33) is d/dt Ca_ss_sr in component Ca (mM).
C RATES(28) is d/dt Ca_ss_CaL in component Ca (mM).
C RATES(46) is d/dt Ca_NSR in component Ca (mM).
C RATES(45) is d/dt Ca_JSR in component Ca (mM).
C RATES(32) is d/dt Na_i in component Na (mM).
C RATES(34) is d/dt Na_ss_sr in component Na (mM).
C RATES(38) is d/dt Cl_i in component Cl (mM).
C RATES(47) is d/dt Cl_ss in component Cl (mM).
C RATES(37) is d/dt K_i in component K (mM).
C RATES(48) is d/dt CaMK_trap in component CaMK_active (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 96487
      CONSTS(2) = 310
      CONSTS(3) = 8314
      CONSTS(4) = 5.4
      CONSTS(5) = 1.8
      CONSTS(6) = 140
      CONSTS(7) = 100
      CONSTS(8) = 0.01
      CONSTS(9) = 0.0011
      STATES(1) = -87.472528
      CONSTS(10) = 0
      CONSTS(11) = 1e3
      CONSTS(12) = 0.5
      CONSTS(13) = -80
      STATES(2) = 0.476653873938053
      STATES(3) = 0.2551771850026437
      STATES(4) = 0.05122852186842027
      STATES(5) = 0.004570874421666423
      STATES(6) = 0.0001529389211563687
      STATES(7) = 0.1328281770413747
      STATES(8) = 0.05333228495514583
      STATES(9) = 0.007137874613383642
      STATES(10) = 0.0003184390970637536
      STATES(11) = 0.01388061042363991
      STATES(12) = 0.003715500164470033
      STATES(13) = 0.0002486371465457659
      STATES(14) = 0.0006446803014431032
      STATES(15) = 8.628258026770862e-5
      STATES(16) = 1.122823524761588e-5
      STATES(17) = 1.051906959485266e-5
      STATES(18) = 2.372219883411999e-6
      STATES(19) = 7.09945e-5
      CONSTS(14) = 1.5552e-4
      STATES(20) = 0.9999613782488198
      STATES(21) = 2.27419494597377e-5
      STATES(22) = 5.936477396599124e-11
      STATES(23) = 1.350122566369778e-15
      STATES(24) = 1.587938121237188e-5
      STATES(25) = 3.611420243201921e-10
      STATES(26) = 0
      STATES(27) = 0
      CONSTS(15) = 1
      CONSTS(16) = 0.000001
      CONSTS(17) = 1
      CONSTS(18) = 0.341
      STATES(28) = 9.61545e-5
      CONSTS(19) = 9.075
      STATES(29) = 0.001022252813581797
      STATES(30) = 0.9913109066867
      STATES(31) = 0.9970771239754529
      CONSTS(20) = 2.6
      CONSTS(21) = 1.5
      CONSTS(22) = 1.4
      CONSTS(23) = 1
      CONSTS(24) = -92
      STATES(32) = 6.935575
      CONSTS(25) = 1.25e-4
      CONSTS(26) = 4.5
      CONSTS(27) = 0.32
      CONSTS(28) = 0.27
      CONSTS(29) = 12.3
      CONSTS(30) = 87.5
      CONSTS(31) = 0.0036
      CONSTS(32) = 1.3
      STATES(33) = 9.04757e-5
      STATES(34) = 6.935607
      CONSTS(33) = 0.00276
      CONSTS(34) = 0.0575
      CONSTS(35) = 0.5e-3
      CONSTS(36) = 1.995084e-7
      CONSTS(37) = 1
      CONSTS(38) = 0.341
      CONSTS(39) = 0.32e-8
      CONSTS(40) = 2.25e-4
      CONSTS(41) = 0.0065
      STATES(35) = 0.001022252813581797
      STATES(36) = 0.3593255076983574
      CONSTS(42) = 600
      CONSTS(43) = 0.01833
      STATES(37) = 145.5689
      STATES(38) = 20.26842
      CONSTS(44) = 2.46108e-5
      CONSTS(45) = 1.77e-5
      STATES(39) = 1.235996442041028e-8
      CONSTS(46) = 0.4975
      STATES(40) = 1.773318268499115e-5
      STATES(41) = 0.9979722
      STATES(42) = 0.9874435
      CONSTS(47) = 9e-7
      CONSTS(48) = -1
      CONSTS(49) = 0.4
      STATES(43) = 0.00107332
      STATES(44) = 0.9960352
      CONSTS(50) = 8
      CONSTS(51) = 1
      CONSTS(52) = 4.75
      CONSTS(53) = 0.1125
      CONSTS(54) = 8
      CONSTS(55) = 10
      CONSTS(56) = 1
      CONSTS(57) = 0.28
      CONSTS(58) = 0.0123
      STATES(45) = 0.917692717
      CONSTS(59) = 0.00092
      CONSTS(60) = 1.7e-4
      CONSTS(61) = 0.004375
      CONSTS(62) = 0.75
      CONSTS(63) = 0.15
      CONSTS(64) = 15
      STATES(46) = 0.929835335
      CONSTS(65) = 25
      CONSTS(66) = 0.2
      CONSTS(67) = 2
      CONSTS(68) = 0.05
      CONSTS(69) = 0.07
      CONSTS(70) = 0.00238
      CONSTS(71) = 0.0005
      CONSTS(72) = 0.8
      CONSTS(73) = 10
      CONSTS(74) = 0.047
      CONSTS(75) = 0.00087
      CONSTS(76) = 1.124
      CONSTS(77) = 0.0087
      STATES(47) = 20.26842
      CONSTS(78) = 0.0015
      STATES(48) = 0.0018861341
      CONSTS(79) = 0.05
      CONSTS(80) = 0.00068
      CONSTS(81) = 0.05
      CONSTS(82) = 0.0895380
      CONSTS(83) =  1000.00* 3.14159265358979*CONSTS(9)*CONSTS(9)*CONSTS(8)
      CONSTS(84) =  0.500000* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(85) =  0.0138542* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(86) =  CONSTS(83)*0.0600000
      CONSTS(87) =  2.00000* 3.14159265358979*CONSTS(9)*CONSTS(9)+ 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(8)
      CONSTS(88) =  CONSTS(87)*2.00000
      CONSTS(89) = ( 1.00000*1.00000*CONSTS(88))/( CONSTS(1)*1.00000)
      CONSTS(90) =  CONSTS(83)*0.678000
      CONSTS(91) =  CONSTS(83)*0.0552000
      CONSTS(92) =  CONSTS(83)*0.00480000
      CONSTS(93) =  CONSTS(83)*0.0200000
      CONSTS(94) =  CONSTS(83)*0.00200000
      CONSTS(95) = CONSTS(4)/(CONSTS(4)+CONSTS(21))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
      RATES(36) = (ALGBRC(8) - STATES(36))/CONSTS(42)
      ALGBRC(2) =  0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(15) =  0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(2) =  - 4.00000*ALGBRC(2)*STATES(2)+ ALGBRC(15)*STATES(3)
      ALGBRC(4) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(17) =  0.0800000*EXP(- STATES(1)/11.0000)
      RATES(29) =  ALGBRC(4)*(1.00000 - STATES(29)) -  ALGBRC(17)*STATES(29)
      ALGBRC(5) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000))
      ALGBRC(18) = TERNRY(STATES(1).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))),  3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))
      RATES(30) =  ALGBRC(5)*(1.00000 - STATES(30)) -  ALGBRC(18)*STATES(30)
      ALGBRC(6) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(1)) -  6.94800e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
      ALGBRC(19) = TERNRY(STATES(1).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
      RATES(31) =  ALGBRC(6)*(1.00000 - STATES(31)) -  ALGBRC(19)*STATES(31)
      ALGBRC(7) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(20) =  0.0800000*EXP(- STATES(1)/11.0000)
      RATES(35) =  ALGBRC(7)*(1.00000 - STATES(35)) -  ALGBRC(20)*STATES(35)
      ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
      ALGBRC(21) = 1.00000/(( 0.000600000*(STATES(1) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(1) - 1.73840)))+( 0.000300000*(STATES(1)+38.3608))/(EXP( 0.152200*(STATES(1)+38.3608)) - 1.00000))
      RATES(39) = (ALGBRC(9) - STATES(39))/ALGBRC(21)
      ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      RATES(41) =  ALGBRC(23)*(1.00000 - STATES(41)) -  ALGBRC(11)*STATES(41)
      ALGBRC(24) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
      ALGBRC(12) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      RATES(42) =  ALGBRC(24)*(1.00000 - STATES(42)) -  ALGBRC(12)*STATES(42)
      ALGBRC(25) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(13) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(30) = ALGBRC(25)/(ALGBRC(25)+ALGBRC(13))
      RATES(44) = (ALGBRC(30) - STATES(44))/CONSTS(50)
      ALGBRC(32) =  0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(27) =  0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(3) =  - ( 3.00000*ALGBRC(2)+ALGBRC(15)+ALGBRC(27))*STATES(3)+ 4.00000*ALGBRC(2)*STATES(2)+ 2.00000*ALGBRC(15)*STATES(4)+ ALGBRC(32)*STATES(7)
      RATES(4) =  - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(15)+ 2.00000*ALGBRC(27))*STATES(4)+ 3.00000*ALGBRC(2)*STATES(3)+ 3.00000*ALGBRC(15)*STATES(5)+ ALGBRC(32)*STATES(8)
      RATES(5) =  - (ALGBRC(2)+ 3.00000*ALGBRC(15)+ 3.00000*ALGBRC(27))*STATES(5)+ 2.00000*ALGBRC(2)*STATES(4)+ 4.00000*ALGBRC(15)*STATES(6)+ ALGBRC(32)*STATES(9)
      RATES(6) =  - ( 4.00000*ALGBRC(15)+ 4.00000*ALGBRC(27))*STATES(6)+ ALGBRC(2)*STATES(5)+ ALGBRC(32)*STATES(10)
      RATES(7) =  - ( 3.00000*ALGBRC(2)+ALGBRC(32))*STATES(7)+ ALGBRC(15)*STATES(8)+ ALGBRC(27)*STATES(3)
      RATES(8) =  - ( 2.00000*ALGBRC(2)+ALGBRC(15)+ALGBRC(27)+ALGBRC(32))*STATES(8)+ 3.00000*ALGBRC(2)*STATES(7)+ 2.00000*ALGBRC(15)*STATES(9)+ 2.00000*ALGBRC(27)*STATES(4)+ 2.00000*ALGBRC(32)*STATES(11)
      RATES(9) =  - (ALGBRC(2)+ 2.00000*ALGBRC(15)+ 2.00000*ALGBRC(27)+ALGBRC(32))*STATES(9)+ 2.00000*ALGBRC(2)*STATES(8)+ 3.00000*ALGBRC(15)*STATES(10)+ 3.00000*ALGBRC(27)*STATES(5)+ 2.00000*ALGBRC(32)*STATES(12)
      RATES(10) =  - ( 3.00000*ALGBRC(15)+ 3.00000*ALGBRC(27)+ALGBRC(32))*STATES(10)+ ALGBRC(2)*STATES(9)+ 4.00000*ALGBRC(27)*STATES(6)+ 2.00000*ALGBRC(32)*STATES(13)
      RATES(11) =  - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(32))*STATES(11)+ ALGBRC(15)*STATES(12)+ ALGBRC(27)*STATES(8)
      RATES(12) =  - (ALGBRC(2)+ALGBRC(15)+ALGBRC(27)+ 2.00000*ALGBRC(32))*STATES(12)+ 2.00000*ALGBRC(2)*STATES(11)+ 2.00000*ALGBRC(15)*STATES(13)+ 2.00000*ALGBRC(27)*STATES(9)+ 3.00000*ALGBRC(32)*STATES(14)
      RATES(13) =  - ( 2.00000*ALGBRC(15)+ 2.00000*ALGBRC(27)+ 2.00000*ALGBRC(32))*STATES(13)+ ALGBRC(2)*STATES(12)+ 3.00000*ALGBRC(27)*STATES(10)+ 3.00000*ALGBRC(32)*STATES(15)
      RATES(14) =  - (ALGBRC(2)+ 3.00000*ALGBRC(32))*STATES(14)+ ALGBRC(15)*STATES(15)+ ALGBRC(27)*STATES(12)
      RATES(15) =  - (ALGBRC(15)+ALGBRC(27)+ 3.00000*ALGBRC(32))*STATES(15)+ ALGBRC(2)*STATES(14)+ 2.00000*ALGBRC(27)*STATES(13)+ 4.00000*ALGBRC(32)*STATES(16)
      ALGBRC(10) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
      ALGBRC(22) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
      ALGBRC(29) = 1.00000/(ALGBRC(10)+ALGBRC(22))
      ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
      RATES(40) = (ALGBRC(34) - STATES(40))/ALGBRC(29)
      ALGBRC(36) =  0.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(16) =  - ( 4.00000*ALGBRC(32)+CONSTS(82))*STATES(16)+ ALGBRC(27)*STATES(15)+ ALGBRC(36)*STATES(17)
      ALGBRC(39) =  0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(41) =  0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(17) =  - (ALGBRC(36)+ALGBRC(41))*STATES(17)+ ALGBRC(39)*STATES(18)+ CONSTS(82)*STATES(16)
      RATES(18) =  - ALGBRC(39)*STATES(18)+ ALGBRC(41)*STATES(17)
      ALGBRC(3) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(16) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(28) = ALGBRC(16)/ALGBRC(3)
      ALGBRC(33) = (1.00000 - ALGBRC(16))/ALGBRC(3)
      ALGBRC(58) = 3.00000/(1.00000+0.00300000/STATES(28) ** 4.00000)
      ALGBRC(40) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+1.00000/( 30.5530*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070))))
      ALGBRC(37) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
      ALGBRC(42) = ALGBRC(37)/ALGBRC(40)
      ALGBRC(44) = (1.00000 - ALGBRC(37))/ALGBRC(40)
      RATES(20) =  - (ALGBRC(28)+ALGBRC(58)+ALGBRC(44))*STATES(20)+ ALGBRC(33)*STATES(21)+ CONSTS(15)*STATES(22)+ ALGBRC(42)*STATES(24)
      RATES(21) =  - (ALGBRC(33)+ALGBRC(58)+ALGBRC(44))*STATES(21)+ ALGBRC(28)*STATES(20)+ CONSTS(15)*STATES(23)+ ALGBRC(42)*STATES(25)
      ALGBRC(48) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(28) ** 4.00000)
      ALGBRC(50) = 1.00000/( ALGBRC(48)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(52) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(50))
      ALGBRC(46) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(54) = ALGBRC(46)/ALGBRC(52)
      ALGBRC(56) = (1.00000 - ALGBRC(46))/ALGBRC(52)
      RATES(22) =  - (ALGBRC(28)+CONSTS(15)+ALGBRC(56))*STATES(22)+ ALGBRC(58)*STATES(20)+ ALGBRC(33)*STATES(23)+ ALGBRC(54)*STATES(26)
      RATES(23) =  - (ALGBRC(33)+CONSTS(15)+ALGBRC(56))*STATES(23)+ ALGBRC(58)*STATES(21)+ ALGBRC(28)*STATES(22)+ ALGBRC(54)*STATES(27)
      ALGBRC(60) = ( CONSTS(16)*ALGBRC(42)*ALGBRC(56)*ALGBRC(58))/( ALGBRC(44)*ALGBRC(54)*CONSTS(15))
      RATES(24) =  - (ALGBRC(28)+ALGBRC(60)+ALGBRC(42))*STATES(24)+ ALGBRC(44)*STATES(20)+ CONSTS(16)*STATES(26)+ ALGBRC(33)*STATES(25)
      RATES(25) =  - (ALGBRC(33)+ALGBRC(60)+ALGBRC(42))*STATES(25)+ ALGBRC(44)*STATES(21)+ CONSTS(16)*STATES(27)+ ALGBRC(28)*STATES(24)
      RATES(26) =  - (ALGBRC(28)+CONSTS(16)+ALGBRC(54))*STATES(26)+ ALGBRC(60)*STATES(24)+ ALGBRC(56)*STATES(22)+ ALGBRC(33)*STATES(27)
      RATES(27) =  - (ALGBRC(33)+CONSTS(16)+ALGBRC(54))*STATES(27)+ ALGBRC(60)*STATES(25)+ ALGBRC(56)*STATES(23)+ ALGBRC(28)*STATES(26)
      ALGBRC(74) = (STATES(46) - STATES(45))/CONSTS(65)
      ALGBRC(77) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(73))/STATES(45)+CONSTS(72) ** 2.00000)
      RATES(45) =  ALGBRC(77)*(ALGBRC(74) - STATES(43))
      ALGBRC(61) =  0.200000*CONSTS(26)*( STATES(34) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(33)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(62) = 1.00000+CONSTS(25)/STATES(33) ** 2.00000
      ALGBRC(63) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(64) =  CONSTS(32)*STATES(34) ** 3.00000+ CONSTS(30) ** 3.00000*STATES(33)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+STATES(33)/CONSTS(31))
      ALGBRC(65) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(34)/CONSTS(29) ** 3.00000)+ STATES(34) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(33)
      ALGBRC(67) = ALGBRC(61)/( ALGBRC(62)*ALGBRC(63)*(ALGBRC(64)+ALGBRC(65)))
      ALGBRC(76) = (STATES(33) - STATES(19))/CONSTS(66)
      ALGBRC(80) = (STATES(33) - STATES(28))/CONSTS(67)
      ALGBRC(84) = 1.00000/(1.00000+( CONSTS(74)*CONSTS(75))/CONSTS(75)+STATES(33) ** 2.00000+( CONSTS(76)*CONSTS(77))/CONSTS(77)+STATES(33) ** 2.00000)
      RATES(33) =  ALGBRC(84)*- ((ALGBRC(76)+ALGBRC(80)) - (( 2.00000*ALGBRC(67)*CONSTS(89))/( 2.00000*CONSTS(93))+( STATES(43)*CONSTS(92))/CONSTS(93)))
      ALGBRC(35) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*STATES(28)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(18)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(38) =  ALGBRC(35)*(STATES(21)+STATES(23))
      ALGBRC(83) = 1.00000/(1.00000+( CONSTS(74)*CONSTS(75))/CONSTS(75)+STATES(28) ** 2.00000+( CONSTS(76)*CONSTS(77))/CONSTS(77)+STATES(28) ** 2.00000)
      RATES(28) =  ALGBRC(83)*- (( ALGBRC(38)*CONSTS(89))/( 2.00000*CONSTS(94)) - ( ALGBRC(80)*CONSTS(93))/CONSTS(94))
      ALGBRC(82) = ( CONSTS(81)*(1.00000 - STATES(48)))/(1.00000+CONSTS(78)/STATES(33))+STATES(48)
      RATES(48) =  CONSTS(79)*ALGBRC(82)*(ALGBRC(82) - STATES(48)) -  CONSTS(80)*STATES(48)
      ALGBRC(69) = ( CONSTS(34)*STATES(19))/(CONSTS(35)+STATES(19))
      ALGBRC(70) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*STATES(19)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(38)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(51) =  0.800000*CONSTS(26)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(19)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(53) = 1.00000+CONSTS(25)/STATES(19) ** 2.00000
      ALGBRC(55) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(57) =  CONSTS(32)*STATES(32) ** 3.00000+ CONSTS(30) ** 3.00000*STATES(19)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+STATES(19)/CONSTS(31))
      ALGBRC(59) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(29) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(19)
      ALGBRC(66) = ALGBRC(51)/( ALGBRC(53)*ALGBRC(55)*(ALGBRC(57)+ALGBRC(59)))
      ALGBRC(90) = ( CONSTS(60)*ALGBRC(82)*1.00000)/(CONSTS(63)+ ALGBRC(82)*1.00000)
      ALGBRC(87) = ( CONSTS(62)*ALGBRC(82)*1.00000)/(CONSTS(63)+ ALGBRC(82)*1.00000)
      ALGBRC(93) = ( (ALGBRC(87)+1.00000)*CONSTS(61)*STATES(19))/((STATES(19)+CONSTS(59)) - ALGBRC(90))
      ALGBRC(72) =  (0.00437500/CONSTS(64))*STATES(46)
      ALGBRC(79) = 1.00000/(1.00000+( CONSTS(68)*CONSTS(70))/STATES(19)+CONSTS(70) ** 2.00000+( CONSTS(71)*CONSTS(69))/STATES(19)+CONSTS(71) ** 2.00000)
      RATES(19) =  ALGBRC(79)*(( - ((ALGBRC(70)+ALGBRC(69)) -  2.00000*ALGBRC(66))*CONSTS(89))/( CONSTS(90)*2.00000)+( (ALGBRC(72) - ALGBRC(93))*CONSTS(91))/CONSTS(90)+( ALGBRC(76)*CONSTS(94))/CONSTS(90))
      RATES(46) = ALGBRC(93) - (ALGBRC(72)+( ALGBRC(74)*CONSTS(92))/CONSTS(91))
      ALGBRC(86) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(82)) ** CONSTS(55))
      ALGBRC(89) =  CONSTS(52)*(1.00000+ALGBRC(86))
      ALGBRC(92) =  ALGBRC(89)*CONSTS(53)
      ALGBRC(97) = ( ALGBRC(38)*ALGBRC(92))/(1.00000+CONSTS(51)/STATES(45) ** CONSTS(54))
      ALGBRC(95) = ALGBRC(89)/(1.00000+CONSTS(58)/STATES(45))
      RATES(43) = - (ALGBRC(97)+STATES(43))/ALGBRC(95)
      ALGBRC(78) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(32))
      ALGBRC(81) =  CONSTS(19)*STATES(29) ** 3.00000*STATES(30)*STATES(31)*(STATES(1) - ALGBRC(78))
      ALGBRC(43) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(43)))
      ALGBRC(47) = STATES(32)/(STATES(32)+CONSTS(20)) ** 3.00000
      ALGBRC(49) =  CONSTS(22)*ALGBRC(45)*CONSTS(95)*ALGBRC(47)
      ALGBRC(85) =  CONSTS(41)*STATES(35) ** 3.00000*STATES(36)*(STATES(1) - ALGBRC(78))
      ALGBRC(73) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(75) = ( CONSTS(1)*CONSTS(39)*ALGBRC(73)*( STATES(32)*EXP(ALGBRC(73)) - CONSTS(6)))/(EXP(ALGBRC(73)) - 1.00000)
      ALGBRC(99) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(38))
      ALGBRC(105) = ( CONSTS(44)*ALGBRC(78) - ALGBRC(99) ** 4.00000)/(ALGBRC(78) - ALGBRC(99) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(107) = (STATES(34) - STATES(32))/CONSTS(66)
      RATES(32) = - (( ( 3.00000*ALGBRC(66)+ 3.00000*ALGBRC(49)+ALGBRC(81)+ALGBRC(85)+ALGBRC(75))*CONSTS(89))/CONSTS(90) - (ALGBRC(105)+( ALGBRC(107)*CONSTS(93))/CONSTS(90)))
      RATES(34) = - (( 3.00000*ALGBRC(67)*CONSTS(89))/CONSTS(93)+ALGBRC(107))
      ALGBRC(1) =  INT(VOI/CONSTS(11))*CONSTS(11)
      ALGBRC(14) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
      ALGBRC(91) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(37))
      ALGBRC(101) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(91)) - 59.2150)))
      ALGBRC(102) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(91))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(91)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(91))+4.75300)))
      ALGBRC(103) = ALGBRC(101)/(ALGBRC(101)+ALGBRC(102))
      ALGBRC(104) =  CONSTS(84)*ALGBRC(103)*(STATES(1) - ALGBRC(91))
      ALGBRC(68) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(94) =  CONSTS(33)*ALGBRC(68)*(STATES(1) - ALGBRC(91))
      ALGBRC(26) =  0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(19) ** 1.40000))
      ALGBRC(31) = STATES(17)+STATES(18)
      ALGBRC(96) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(37)+ CONSTS(43)*STATES(32)))
      ALGBRC(98) =  ALGBRC(26)*ALGBRC(31)*(STATES(1) - ALGBRC(96))
      ALGBRC(108) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(109) =  CONSTS(85)*STATES(39)*ALGBRC(108)*(STATES(1) - ALGBRC(91))
      ALGBRC(110) = EXP(STATES(1)/550.000)
      ALGBRC(111) =  CONSTS(46)*STATES(40) ** 3.00000*STATES(41)*STATES(42)*ALGBRC(110)*(STATES(1) - ALGBRC(91))
      ALGBRC(106) = ( CONSTS(45)*(ALGBRC(91) - ALGBRC(99)))/((ALGBRC(91)+87.8251) - ALGBRC(99))
      RATES(37) = - (( ((ALGBRC(98)+ALGBRC(109)+ALGBRC(104)+ALGBRC(111)+ALGBRC(94)+ALGBRC(14)) -  2.00000*ALGBRC(49))*CONSTS(89))/CONSTS(90) - ALGBRC(106))
      ALGBRC(71) = (ALGBRC(38)+ALGBRC(70)+ALGBRC(69)) -  2.00000*(ALGBRC(66)+ALGBRC(67))
      ALGBRC(88) = ALGBRC(81)+ 3.00000*(ALGBRC(66)+ALGBRC(67))+ 3.00000*ALGBRC(49)+ALGBRC(85)+ALGBRC(75)
      ALGBRC(112) = ALGBRC(109)+ALGBRC(98)+ALGBRC(104)+ALGBRC(94)+ - 2.00000*ALGBRC(49)+ALGBRC(111)+ALGBRC(14)
      ALGBRC(113) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(38) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(114) = 1.00000 - 1.00000/(1.00000+STATES(43)/CONSTS(49) ** 2.00000)
      ALGBRC(115) =  ALGBRC(113)*STATES(44)*ALGBRC(114)
      ALGBRC(100) =  CONSTS(40)*(STATES(1) - ALGBRC(99))
      ALGBRC(116) = ALGBRC(100)+ALGBRC(115)
      RATES(1) = - (ALGBRC(88)+ALGBRC(112)+ALGBRC(71)+ALGBRC(116))
      ALGBRC(117) = (STATES(47) - STATES(38))/CONSTS(66)
      RATES(38) = - (( ALGBRC(100)*CONSTS(89))/( - 1.00000*CONSTS(90)) - (ALGBRC(105)+ALGBRC(106)+( ALGBRC(117)*CONSTS(93))/CONSTS(90)))
      RATES(47) = - (( ALGBRC(115)*CONSTS(89))/( - 1.00000*CONSTS(93))+ALGBRC(117))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
      ALGBRC(2) =  0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(15) =  0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(4) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(17) =  0.0800000*EXP(- STATES(1)/11.0000)
      ALGBRC(5) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000))
      ALGBRC(18) = TERNRY(STATES(1).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))),  3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))
      ALGBRC(6) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(1)) -  6.94800e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
      ALGBRC(19) = TERNRY(STATES(1).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
      ALGBRC(7) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(20) =  0.0800000*EXP(- STATES(1)/11.0000)
      ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
      ALGBRC(21) = 1.00000/(( 0.000600000*(STATES(1) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(1) - 1.73840)))+( 0.000300000*(STATES(1)+38.3608))/(EXP( 0.152200*(STATES(1)+38.3608)) - 1.00000))
      ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(24) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
      ALGBRC(12) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(25) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(13) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(30) = ALGBRC(25)/(ALGBRC(25)+ALGBRC(13))
      ALGBRC(32) =  0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(27) =  0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(10) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
      ALGBRC(22) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
      ALGBRC(29) = 1.00000/(ALGBRC(10)+ALGBRC(22))
      ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
      ALGBRC(36) =  0.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(39) =  0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(41) =  0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(3) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(16) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(28) = ALGBRC(16)/ALGBRC(3)
      ALGBRC(33) = (1.00000 - ALGBRC(16))/ALGBRC(3)
      ALGBRC(58) = 3.00000/(1.00000+0.00300000/STATES(28) ** 4.00000)
      ALGBRC(40) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+1.00000/( 30.5530*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070))))
      ALGBRC(37) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
      ALGBRC(42) = ALGBRC(37)/ALGBRC(40)
      ALGBRC(44) = (1.00000 - ALGBRC(37))/ALGBRC(40)
      ALGBRC(48) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(28) ** 4.00000)
      ALGBRC(50) = 1.00000/( ALGBRC(48)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(52) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(50))
      ALGBRC(46) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(54) = ALGBRC(46)/ALGBRC(52)
      ALGBRC(56) = (1.00000 - ALGBRC(46))/ALGBRC(52)
      ALGBRC(60) = ( CONSTS(16)*ALGBRC(42)*ALGBRC(56)*ALGBRC(58))/( ALGBRC(44)*ALGBRC(54)*CONSTS(15))
      ALGBRC(74) = (STATES(46) - STATES(45))/CONSTS(65)
      ALGBRC(77) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(73))/STATES(45)+CONSTS(72) ** 2.00000)
      ALGBRC(61) =  0.200000*CONSTS(26)*( STATES(34) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(33)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(62) = 1.00000+CONSTS(25)/STATES(33) ** 2.00000
      ALGBRC(63) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(64) =  CONSTS(32)*STATES(34) ** 3.00000+ CONSTS(30) ** 3.00000*STATES(33)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+STATES(33)/CONSTS(31))
      ALGBRC(65) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(34)/CONSTS(29) ** 3.00000)+ STATES(34) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(33)
      ALGBRC(67) = ALGBRC(61)/( ALGBRC(62)*ALGBRC(63)*(ALGBRC(64)+ALGBRC(65)))
      ALGBRC(76) = (STATES(33) - STATES(19))/CONSTS(66)
      ALGBRC(80) = (STATES(33) - STATES(28))/CONSTS(67)
      ALGBRC(84) = 1.00000/(1.00000+( CONSTS(74)*CONSTS(75))/CONSTS(75)+STATES(33) ** 2.00000+( CONSTS(76)*CONSTS(77))/CONSTS(77)+STATES(33) ** 2.00000)
      ALGBRC(35) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*STATES(28)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(18)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(38) =  ALGBRC(35)*(STATES(21)+STATES(23))
      ALGBRC(83) = 1.00000/(1.00000+( CONSTS(74)*CONSTS(75))/CONSTS(75)+STATES(28) ** 2.00000+( CONSTS(76)*CONSTS(77))/CONSTS(77)+STATES(28) ** 2.00000)
      ALGBRC(82) = ( CONSTS(81)*(1.00000 - STATES(48)))/(1.00000+CONSTS(78)/STATES(33))+STATES(48)
      ALGBRC(69) = ( CONSTS(34)*STATES(19))/(CONSTS(35)+STATES(19))
      ALGBRC(70) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*STATES(19)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(38)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(51) =  0.800000*CONSTS(26)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(19)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(53) = 1.00000+CONSTS(25)/STATES(19) ** 2.00000
      ALGBRC(55) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(57) =  CONSTS(32)*STATES(32) ** 3.00000+ CONSTS(30) ** 3.00000*STATES(19)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+STATES(19)/CONSTS(31))
      ALGBRC(59) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(29) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(19)
      ALGBRC(66) = ALGBRC(51)/( ALGBRC(53)*ALGBRC(55)*(ALGBRC(57)+ALGBRC(59)))
      ALGBRC(90) = ( CONSTS(60)*ALGBRC(82)*1.00000)/(CONSTS(63)+ ALGBRC(82)*1.00000)
      ALGBRC(87) = ( CONSTS(62)*ALGBRC(82)*1.00000)/(CONSTS(63)+ ALGBRC(82)*1.00000)
      ALGBRC(93) = ( (ALGBRC(87)+1.00000)*CONSTS(61)*STATES(19))/((STATES(19)+CONSTS(59)) - ALGBRC(90))
      ALGBRC(72) =  (0.00437500/CONSTS(64))*STATES(46)
      ALGBRC(79) = 1.00000/(1.00000+( CONSTS(68)*CONSTS(70))/STATES(19)+CONSTS(70) ** 2.00000+( CONSTS(71)*CONSTS(69))/STATES(19)+CONSTS(71) ** 2.00000)
      ALGBRC(86) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(82)) ** CONSTS(55))
      ALGBRC(89) =  CONSTS(52)*(1.00000+ALGBRC(86))
      ALGBRC(92) =  ALGBRC(89)*CONSTS(53)
      ALGBRC(97) = ( ALGBRC(38)*ALGBRC(92))/(1.00000+CONSTS(51)/STATES(45) ** CONSTS(54))
      ALGBRC(95) = ALGBRC(89)/(1.00000+CONSTS(58)/STATES(45))
      ALGBRC(78) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(32))
      ALGBRC(81) =  CONSTS(19)*STATES(29) ** 3.00000*STATES(30)*STATES(31)*(STATES(1) - ALGBRC(78))
      ALGBRC(43) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(43)))
      ALGBRC(47) = STATES(32)/(STATES(32)+CONSTS(20)) ** 3.00000
      ALGBRC(49) =  CONSTS(22)*ALGBRC(45)*CONSTS(95)*ALGBRC(47)
      ALGBRC(85) =  CONSTS(41)*STATES(35) ** 3.00000*STATES(36)*(STATES(1) - ALGBRC(78))
      ALGBRC(73) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(75) = ( CONSTS(1)*CONSTS(39)*ALGBRC(73)*( STATES(32)*EXP(ALGBRC(73)) - CONSTS(6)))/(EXP(ALGBRC(73)) - 1.00000)
      ALGBRC(99) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(38))
      ALGBRC(105) = ( CONSTS(44)*ALGBRC(78) - ALGBRC(99) ** 4.00000)/(ALGBRC(78) - ALGBRC(99) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(107) = (STATES(34) - STATES(32))/CONSTS(66)
      ALGBRC(1) =  INT(VOI/CONSTS(11))*CONSTS(11)
      ALGBRC(14) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
      ALGBRC(91) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(37))
      ALGBRC(101) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(91)) - 59.2150)))
      ALGBRC(102) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(91))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(91)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(91))+4.75300)))
      ALGBRC(103) = ALGBRC(101)/(ALGBRC(101)+ALGBRC(102))
      ALGBRC(104) =  CONSTS(84)*ALGBRC(103)*(STATES(1) - ALGBRC(91))
      ALGBRC(68) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(94) =  CONSTS(33)*ALGBRC(68)*(STATES(1) - ALGBRC(91))
      ALGBRC(26) =  0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(19) ** 1.40000))
      ALGBRC(31) = STATES(17)+STATES(18)
      ALGBRC(96) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(37)+ CONSTS(43)*STATES(32)))
      ALGBRC(98) =  ALGBRC(26)*ALGBRC(31)*(STATES(1) - ALGBRC(96))
      ALGBRC(108) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(109) =  CONSTS(85)*STATES(39)*ALGBRC(108)*(STATES(1) - ALGBRC(91))
      ALGBRC(110) = EXP(STATES(1)/550.000)
      ALGBRC(111) =  CONSTS(46)*STATES(40) ** 3.00000*STATES(41)*STATES(42)*ALGBRC(110)*(STATES(1) - ALGBRC(91))
      ALGBRC(106) = ( CONSTS(45)*(ALGBRC(91) - ALGBRC(99)))/((ALGBRC(91)+87.8251) - ALGBRC(99))
      ALGBRC(71) = (ALGBRC(38)+ALGBRC(70)+ALGBRC(69)) -  2.00000*(ALGBRC(66)+ALGBRC(67))
      ALGBRC(88) = ALGBRC(81)+ 3.00000*(ALGBRC(66)+ALGBRC(67))+ 3.00000*ALGBRC(49)+ALGBRC(85)+ALGBRC(75)
      ALGBRC(112) = ALGBRC(109)+ALGBRC(98)+ALGBRC(104)+ALGBRC(94)+ - 2.00000*ALGBRC(49)+ALGBRC(111)+ALGBRC(14)
      ALGBRC(113) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(38) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(114) = 1.00000 - 1.00000/(1.00000+STATES(43)/CONSTS(49) ** 2.00000)
      ALGBRC(115) =  ALGBRC(113)*STATES(44)*ALGBRC(114)
      ALGBRC(100) =  CONSTS(40)*(STATES(1) - ALGBRC(99))
      ALGBRC(116) = ALGBRC(100)+ALGBRC(115)
      ALGBRC(117) = (STATES(47) - STATES(38))/CONSTS(66)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END