C
C There are a total of 120 entries in the algebraic variable array.
C There are a total of 46 entries in each of the rate and state variable arrays.
C There are a total of 99 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(94) is Vmyo in component cell_geometry (uL).
C CONSTS(95) is Vnsr in component cell_geometry (uL).
C CONSTS(90) is Vsr in component cell_geometry (uL).
C CONSTS(96) is Vjsr in component cell_geometry (uL).
C CONSTS(97) is Vss_sr in component cell_geometry (uL).
C CONSTS(98) is Vss_CaL in component cell_geometry (uL).
C CONSTS(92) is Acap in component cell_geometry (cm2).
C CONSTS(93) 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(87) is vcell in component cell_geometry (uL).
C CONSTS(91) is ageo in component cell_geometry (cm2).
C STATES(1) is Vm in component membrane (mV).
C ALGBRC(13) 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(74) is caiont in component membrane (uA_per_uF).
C ALGBRC(91) is naiont in component membrane (uA_per_uF).
C ALGBRC(115) is kiont in component membrane (uA_per_uF).
C ALGBRC(119) is clont in component membrane (uA_per_uF).
C ALGBRC(84) is INa in component INa (uA_per_uF).
C ALGBRC(37) is ICaL in component ICaL (uA_per_uF).
C ALGBRC(107) is IK1 in component IK1 (uA_per_uF).
C ALGBRC(97) is IKp in component IKp (uA_per_uF).
C ALGBRC(101) is IKs in component IKs (uA_per_uF).
C ALGBRC(112) is IKr in component IKr (uA_per_uF).
C ALGBRC(72) is IpCa in component IpCa (uA_per_uF).
C ALGBRC(73) is ICab in component ICab (uA_per_uF).
C ALGBRC(70) is INaCa in component INaCa (uA_per_uF).
C ALGBRC(68) is INaCa_ss_sr in component INaCa (uA_per_uF).
C ALGBRC(51) is INaK in component INaK (uA_per_uF).
C ALGBRC(114) is Ito1 in component Ito1 (uA_per_uF).
C ALGBRC(118) is Ito2 in component Ito2 (uA_per_uF).
C ALGBRC(103) is IClb in component IClb (uA_per_uF).
C ALGBRC(88) is INaL in component INaL (uA_per_uF).
C ALGBRC(78) is INab in component INab (uA_per_uF).
C ALGBRC(24) is gKs in component IKs (mS_per_uF).
C CONSTS(14) is gKs_max 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 ALGBRC(28) is O2 in component IKs (dimensionless).
C ALGBRC(31) is OKs in component IKs (dimensionless).
C ALGBRC(32) is alpha in component IKs (per_ms).
C ALGBRC(35) is beta in component IKs (per_ms).
C ALGBRC(40) is delta in component IKs (per_ms).
C ALGBRC(38) is gamma in component IKs (per_ms).
C ALGBRC(43) is eta in component IKs (per_ms).
C CONSTS(86) is theta in component IKs (per_ms).
C ALGBRC(46) is omega in component IKs (per_ms).
C ALGBRC(49) is psi in component IKs (per_ms).
C ALGBRC(99) is EKs in component reversal_potentials (mV).
C STATES(18) is Ca_i in component Ca (mM).
C ALGBRC(34) is ICaL_max in component ICaL (uA_per_uF).
C CONSTS(15) is PCa in component ICaL (L_per_F_ms).
C STATES(19) is C in component ICaL (dimensionless).
C STATES(20) is O in component ICaL (dimensionless).
C STATES(21) is C_star in component ICaL (dimensionless).
C STATES(22) is O_star in component ICaL (dimensionless).
C STATES(23) is CI in component ICaL (dimensionless).
C STATES(24) is OI in component ICaL (dimensionless).
C STATES(25) is CI_star in component ICaL (dimensionless).
C ALGBRC(2) is OI_star in component ICaL (dimensionless).
C ALGBRC(14) is ACT_tau in component ICaL (ms).
C ALGBRC(39) is IV_tau in component ICaL (ms).
C ALGBRC(54) is IV_tau_star in component ICaL (ms).
C ALGBRC(25) is ACT_infinity in component ICaL (dimensionless).
C ALGBRC(36) is IV_infinity in component ICaL (dimensionless).
C ALGBRC(47) is IV_infinity_star in component ICaL (dimensionless).
C ALGBRC(52) is IV_beta_star in component ICaL (dimensionless).
C ALGBRC(50) is IV_beta_infinity_star in component ICaL (dimensionless).
C ALGBRC(29) is alpha in component ICaL (per_ms).
C ALGBRC(33) is beta in component ICaL (per_ms).
C ALGBRC(60) is delta in component ICaL (per_ms).
C ALGBRC(41) is x in component ICaL (per_ms).
C ALGBRC(44) is y in component ICaL (per_ms).
C CONSTS(16) is theta in component ICaL (per_ms).
C ALGBRC(62) is delta_I in component ICaL (per_ms).
C ALGBRC(56) is x_star in component ICaL (per_ms).
C ALGBRC(58) is y_star in component ICaL (per_ms).
C CONSTS(17) is theta_I in component ICaL (per_ms).
C CONSTS(18) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(19) is gamma_Cao in component ICaL (dimensionless).
C STATES(26) is Ca_ss_CaL in component Ca (mM).
C CONSTS(20) is g_Na in component INa (mS_per_uF).
C ALGBRC(81) is ENa in component reversal_potentials (mV).
C STATES(27) is m in component INa_m_gate (dimensionless).
C STATES(28) is h in component INa_h_gate (dimensionless).
C STATES(29) is j in component INa_j_gate (dimensionless).
C ALGBRC(3) is am in component INa_m_gate (per_ms).
C ALGBRC(15) is bm in component INa_m_gate (per_ms).
C ALGBRC(4) is ah in component INa_h_gate (per_ms).
C ALGBRC(16) is bh in component INa_h_gate (per_ms).
C ALGBRC(5) is aj in component INa_j_gate (per_ms).
C ALGBRC(17) is bj in component INa_j_gate (per_ms).
C CONSTS(21) is kmnai in component INaK (mM).
C CONSTS(22) is kmko in component INaK (mM).
C CONSTS(23) is ibarnak in component INaK (uA_per_uF).
C CONSTS(24) is sigma in component INaK (dimensionless).
C CONSTS(25) is V_half in component INaK (mV).
C ALGBRC(42) is phi in component INaK (dimensionless).
C ALGBRC(45) is fv in component INaK (dimensionless).
C CONSTS(99) is PK in component INaK (dimensionless).
C ALGBRC(48) is PNa in component INaK (dimensionless).
C STATES(30) is Na_i in component Na (mM).
C ALGBRC(69) is INaCa_cai in component INaCa (uA_per_uF).
C CONSTS(26) is KmCa_act in component INaCa (mM).
C CONSTS(27) is Vmax in component INaCa (uA_per_uF).
C CONSTS(28) is ksat in component INaCa (dimensionless).
C CONSTS(29) is eta in component INaCa (dimensionless).
C CONSTS(30) is KmNai in component INaCa (mM).
C CONSTS(31) is KmNao in component INaCa (mM).
C CONSTS(32) is KmCai in component INaCa (mM).
C CONSTS(33) is KmCao in component INaCa (mM).
C STATES(31) is Ca_ss_sr in component Ca (mM).
C STATES(32) is Na_ss_sr in component Na (mM).
C ALGBRC(53) is numerator in component INaCa (mM4_uA_per_uF).
C ALGBRC(63) is num_ss in component INaCa (mM4_uA_per_uF).
C ALGBRC(55) is denom_1 in component INaCa (dimensionless).
C ALGBRC(57) is denom_2 in component INaCa (dimensionless).
C ALGBRC(59) is denom_3 in component INaCa (mM4).
C ALGBRC(61) is denom_4 in component INaCa (mM4).
C ALGBRC(64) is denom_ss_1 in component INaCa (dimensionless).
C ALGBRC(65) is denom_ss_2 in component INaCa (dimensionless).
C ALGBRC(66) is denom_ss_3 in component INaCa (mM4).
C ALGBRC(67) is denom_ss_4 in component INaCa (mM4).
C CONSTS(34) is gKp in component IKp (mS_per_uF).
C ALGBRC(71) is Kp in component IKp (dimensionless).
C ALGBRC(94) is EK in component reversal_potentials (mV).
C CONSTS(35) is gpCa in component IpCa (uA_per_uF).
C CONSTS(36) is KmpCa in component IpCa (mM).
C CONSTS(37) is PCab in component ICab (L_per_F_ms).
C CONSTS(38) is gamma_Ca_i in component ICab (dimensionless).
C CONSTS(39) is gamma_Ca_o in component ICab (dimensionless).
C CONSTS(40) is PNab in component INab (L_per_F_ms).
C ALGBRC(76) is phi in component INab (dimensionless).
C CONSTS(41) is gClb in component IClb (mS_per_uF).
C ALGBRC(102) is ECl in component reversal_potentials (mV).
C CONSTS(42) is gNaL in component INaL (mS_per_uF).
C STATES(33) is mL in component INaL_mL_gate (dimensionless).
C STATES(34) is hL in component INaL_hL_gate (dimensionless).
C ALGBRC(6) is amL in component INaL_mL_gate (per_ms).
C ALGBRC(18) is bmL in component INaL_mL_gate (per_ms).
C CONSTS(43) is tau_hL in component INaL_hL_gate (ms).
C ALGBRC(7) is hL_infinity in component INaL_hL_gate (dimensionless).
C CONSTS(44) is prnak in component reversal_potentials (dimensionless).
C STATES(35) is K_i in component K (mM).
C STATES(36) is Cl_i in component Cl (mM).
C CONSTS(45) is g_K1_max in component IK1 (mS_per_uF).
C CONSTS(88) is gK1 in component IK1 (mS_per_uF).
C ALGBRC(106) is K1 in component IK1_K1_gate (dimensionless).
C ALGBRC(104) is alpha_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(105) is beta_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(108) is CT_Na_Cl in component CT_Na_Cl (mM_per_ms).
C CONSTS(46) is CT_Na_Cl_max in component CT_Na_Cl (mM_per_ms).
C ALGBRC(109) is CT_K_Cl in component CT_K_Cl (mM_per_ms).
C CONSTS(47) is CT_K_Cl_max in component CT_K_Cl (mM_per_ms).
C CONSTS(89) is gKr in component IKr (mS_per_uF).
C CONSTS(48) is gKr_max in component IKr (mS_per_uF).
C ALGBRC(111) is r in component IKr_r_gate (dimensionless).
C STATES(37) is xr in component IKr_xr_gate (dimensionless).
C ALGBRC(8) is xr_infinity in component IKr_xr_gate (dimensionless).
C ALGBRC(19) is tau_xr in component IKr_xr_gate (ms).
C CONSTS(49) is gto1 in component Ito1 (mS_per_uF).
C ALGBRC(113) is rto1 in component Ito1 (dimensionless).
C STATES(38) is a in component Ito1_a_gate (dimensionless).
C STATES(39) is i1f in component Ito1_i1f_gate (dimensionless).
C STATES(40) is i1s in component Ito1_i1s_gate (dimensionless).
C ALGBRC(9) is alpha_a in component Ito1_a_gate (per_ms).
C ALGBRC(20) is beta_a in component Ito1_a_gate (per_ms).
C ALGBRC(26) is tau_a in component Ito1_a_gate (ms).
C ALGBRC(30) is a_infinity in component Ito1_a_gate (dimensionless).
C ALGBRC(21) is alpha_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(10) is beta_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(22) is alpha_i1s in component Ito1_i1s_gate (per_ms).
C ALGBRC(11) is beta_i1s in component Ito1_i1s_gate (per_ms).
C CONSTS(50) is PCl in component Ito2 (L_per_F_ms).
C CONSTS(51) is zCl in component Ito2 (dimensionless).
C ALGBRC(116) is Ito2_max in component Ito2 (uA_per_uF).
C CONSTS(52) is Ito2_max_scaling_factor in component Ito2 (dimensionless).
C ALGBRC(117) is KCa_ito2 in component Ito2 (dimensionless).
C CONSTS(53) is kCa_ito2 in component Ito2 (mM_per_ms).
C STATES(41) is Irel in component Irel (mM_per_ms).
C STATES(42) is i2f in component Ito2_i2f_gate (dimensionless).
C ALGBRC(23) is alpha_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(12) is beta_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(27) is i2f_infinity in component Ito2_i2f_gate (dimensionless).
C CONSTS(54) is tau_i2f in component Ito2_i2f_gate (ms).
C ALGBRC(100) is rel_infinity in component Irel (mM_per_ms).
C ALGBRC(98) is tau_rel in component Irel (ms).
C CONSTS(55) is Krel_infinity in component Irel (mM).
C ALGBRC(95) is alpha_rel in component Irel (mM_per_uA_per_uF_per_ms).
C ALGBRC(92) is beta_tau in component Irel (ms).
C CONSTS(56) is beta_0 in component Irel (ms).
C CONSTS(57) is kappa in component Irel (mM_per_uA_per_uF_per_ms_per_ms).
C CONSTS(58) is h_rel in component Irel (dimensionless).
C CONSTS(59) is h_beta in component Irel (dimensionless).
C ALGBRC(89) is delta_beta_CaMK in component Irel (dimensionless).
C CONSTS(60) is delta_beta_0 in component Irel (dimensionless).
C CONSTS(61) is K_beta in component Irel (mM).
C CONSTS(62) is Krel_tau in component Irel (mM).
C STATES(43) is Ca_JSR in component Ca (mM).
C ALGBRC(85) is CaMK_active in component CaMK_active (dimensionless).
C ALGBRC(96) is Iup in component Iup (mM_per_ms).
C CONSTS(63) is kmup in component Iup (mM).
C ALGBRC(93) is delta_kmPLB in component Iup (mM).
C CONSTS(64) is delta_kmPLB_bar in component Iup (mM).
C CONSTS(65) is iupbar in component Iup (mM_per_ms).
C ALGBRC(90) is delta_iupCaMK in component Iup (dimensionless).
C CONSTS(66) is delta_iupCaMK_bar in component Iup (dimensionless).
C CONSTS(67) is kmCaMK in component Iup (mM).
C ALGBRC(75) is Ileak in component Ileak (mM_per_ms).
C CONSTS(68) is nsrbar in component Ileak (mM).
C STATES(44) is Ca_NSR in component Ca (mM).
C ALGBRC(77) is Itr in component Itr (mM_per_ms).
C CONSTS(69) is tautr in component Itr (ms).
C CONSTS(70) is tau_diff in component Ca (ms).
C CONSTS(71) is tau_diff_ss in component Ca (ms).
C ALGBRC(79) is Idiff in component Ca (mM_per_ms).
C ALGBRC(83) is Idiff_ss in component Ca (mM_per_ms).
C ALGBRC(82) is bmyo in component Ca (dimensionless).
C ALGBRC(86) is bss_sr in component Ca (dimensionless).
C ALGBRC(87) is bss_cal in component Ca (dimensionless).
C ALGBRC(80) is bcsqn in component Ca (dimensionless).
C CONSTS(72) is cmdn_bar in component Ca (mM).
C CONSTS(73) is trpn_bar in component Ca (mM).
C CONSTS(74) is km_cmdn in component Ca (mM).
C CONSTS(75) is km_trpn in component Ca (mM).
C CONSTS(76) is kmcsqn in component Ca (mM).
C CONSTS(77) is csqnbar in component Ca (mM).
C CONSTS(78) is BSRmax in component Ca (mM).
C CONSTS(79) is KmBSR in component Ca (mM).
C CONSTS(80) is BSLmax in component Ca (mM).
C CONSTS(81) is KmBSL in component Ca (mM).
C ALGBRC(110) is Idiff_Na in component Na (mM_per_ms).
C STATES(45) is Cl_ss in component Cl (mM).
C ALGBRC(120) is Idiff_Cl in component Cl (mM_per_ms).
C CONSTS(82) is Km in component CaMK_active (mM).
C STATES(46) is CaMK_trap in component CaMK_active (dimensionless).
C CONSTS(83) is alpha_CaMK in component CaMK_active (per_ms).
C CONSTS(84) is beta_CaMK in component CaMK_active (per_ms).
C CONSTS(85) 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(19) is d/dt C in component ICaL (dimensionless).
C RATES(20) is d/dt O in component ICaL (dimensionless).
C RATES(21) is d/dt C_star in component ICaL (dimensionless).
C RATES(22) is d/dt O_star in component ICaL (dimensionless).
C RATES(23) is d/dt CI in component ICaL (dimensionless).
C RATES(24) is d/dt OI in component ICaL (dimensionless).
C RATES(25) is d/dt CI_star in component ICaL (dimensionless).
C RATES(27) is d/dt m in component INa_m_gate (dimensionless).
C RATES(28) is d/dt h in component INa_h_gate (dimensionless).
C RATES(29) is d/dt j in component INa_j_gate (dimensionless).
C RATES(33) is d/dt mL in component INaL_mL_gate (dimensionless).
C RATES(34) is d/dt hL in component INaL_hL_gate (dimensionless).
C RATES(37) is d/dt xr in component IKr_xr_gate (dimensionless).
C RATES(38) is d/dt a in component Ito1_a_gate (dimensionless).
C RATES(39) is d/dt i1f in component Ito1_i1f_gate (dimensionless).
C RATES(40) is d/dt i1s in component Ito1_i1s_gate (dimensionless).
C RATES(42) is d/dt i2f in component Ito2_i2f_gate (dimensionless).
C RATES(41) is d/dt Irel in component Irel (mM_per_ms).
C RATES(18) is d/dt Ca_i in component Ca (mM).
C RATES(31) is d/dt Ca_ss_sr in component Ca (mM).
C RATES(26) is d/dt Ca_ss_CaL in component Ca (mM).
C RATES(44) is d/dt Ca_NSR in component Ca (mM).
C RATES(43) is d/dt Ca_JSR in component Ca (mM).
C RATES(30) is d/dt Na_i in component Na (mM).
C RATES(32) is d/dt Na_ss_sr in component Na (mM).
C RATES(36) is d/dt Cl_i in component Cl (mM).
C RATES(45) is d/dt Cl_ss in component Cl (mM).
C RATES(35) is d/dt K_i in component K (mM).
C RATES(46) is d/dt CaMK_trap in component CaMK_active (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 96485
      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.4947322902196
      CONSTS(10) = 0
      CONSTS(11) = 1000
      CONSTS(12) = 0.5
      CONSTS(13) = -80
      CONSTS(14) = 0.0826
      STATES(2) = 0.476679784307081
      STATES(3) = 0.255181042718489
      STATES(4) = 0.0512272866698557
      STATES(5) = 0.0045705849736195
      STATES(6) = 0.000152923243503507
      STATES(7) = 0.132814156809568
      STATES(8) = 0.053324572164618
      STATES(9) = 0.00713656414955778
      STATES(10) = 0.000318368292970148
      STATES(11) = 0.0138769647384953
      STATES(12) = 0.00371438476499172
      STATES(13) = 0.000248553788631871
      STATES(14) = 0.000644426171899078
      STATES(15) = 8.62495047666861e-5
      STATES(16) = 1.12344645485893e-5
      STATES(17) = 1.05292905767805e-5
      STATES(18) = 8.38749860095482e-5
      CONSTS(15) = 0.00015552
      STATES(19) = 0.99996982093287
      STATES(20) = 2.26882961412491e-5
      STATES(21) = 7.49117189549873e-6
      STATES(22) = 1.23556404128518e-10
      STATES(23) = 5.89853003114132e-11
      STATES(24) = 8.01952813127322e-17
      STATES(25) = 1.9237001017021e-9
      CONSTS(16) = 1
      CONSTS(17) = 1e-6
      CONSTS(18) = 1
      CONSTS(19) = 0.341
      STATES(26) = 0.000119239795531548
      CONSTS(20) = 9.075
      STATES(27) = 0.00101841101248412
      STATES(28) = 0.991354143386887
      STATES(29) = 0.99709036251243
      CONSTS(21) = 2.6
      CONSTS(22) = 1.5
      CONSTS(23) = 1.4
      CONSTS(24) = 1
      CONSTS(25) = -92
      STATES(30) = 8.01031392415197
      CONSTS(26) = 0.000125
      CONSTS(27) = 4.5
      CONSTS(28) = 0.32
      CONSTS(29) = 0.27
      CONSTS(30) = 12.3
      CONSTS(31) = 87.5
      CONSTS(32) = 0.0036
      CONSTS(33) = 1.3
      STATES(31) = 0.00011296707102434
      STATES(32) = 8.01033434859049
      CONSTS(34) = 0.00276
      CONSTS(35) = 0.0575
      CONSTS(36) = 0.0005
      CONSTS(37) = 1.99508e-7
      CONSTS(38) = 1
      CONSTS(39) = 0.341
      CONSTS(40) = 3.2e-9
      CONSTS(41) = 0.000225
      CONSTS(42) = 0.0065
      STATES(33) = 0.00101841101248412
      STATES(34) = 0.319611809402137
      CONSTS(43) = 600
      CONSTS(44) = 0.01833
      STATES(35) = 144.907134343314
      STATES(36) = 20.6683603724432
      CONSTS(45) = 0.5
      CONSTS(46) = 2.46108e-5
      CONSTS(47) = 1.77e-5
      CONSTS(48) = 0.0138542
      STATES(37) = 4.06959064407486e-6
      CONSTS(49) = 0.497458
      STATES(38) = 1.76776038804228e-5
      STATES(39) = 0.997976503910799
      STATES(40) = 0.947581184798815
      CONSTS(50) = 9e-7
      CONSTS(51) = -1
      CONSTS(52) = 1
      CONSTS(53) = 0.4
      STATES(41) = 0.000714227102027745
      STATES(42) = 0.996044894628276
      CONSTS(54) = 8
      CONSTS(55) = 1
      CONSTS(56) = 4.75
      CONSTS(57) = 0.1125
      CONSTS(58) = 8
      CONSTS(59) = 10
      CONSTS(60) = 1
      CONSTS(61) = 0.28
      CONSTS(62) = 0.0123
      STATES(43) = 0.993914988616979
      CONSTS(63) = 0.00092
      CONSTS(64) = 0.00017
      CONSTS(65) = 0.004375
      CONSTS(66) = 0.75
      CONSTS(67) = 0.15
      CONSTS(68) = 15
      STATES(44) = 1.01475649943057
      CONSTS(69) = 25
      CONSTS(70) = 0.2
      CONSTS(71) = 2
      CONSTS(72) = 0.05
      CONSTS(73) = 0.07
      CONSTS(74) = 0.00238
      CONSTS(75) = 0.0005
      CONSTS(76) = 0.8
      CONSTS(77) = 10
      CONSTS(78) = 0.047
      CONSTS(79) = 0.00087
      CONSTS(80) = 1.124
      CONSTS(81) = 0.0087
      STATES(45) = 20.6683605553303
      CONSTS(82) = 0.0015
      STATES(46) = 0.0114911034744777
      CONSTS(83) = 0.05
      CONSTS(84) = 0.00068
      CONSTS(85) = 0.05
      CONSTS(86) = 0.0895383
      CONSTS(87) =  1000.00* 3.14159265358979*CONSTS(9)*CONSTS(9)*CONSTS(8)
      CONSTS(88) =  CONSTS(45)* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(89) =  CONSTS(48)* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(90) =  CONSTS(87)*0.0600000
      CONSTS(91) =  2.00000* 3.14159265358979*CONSTS(9)*CONSTS(9)+ 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(8)
      CONSTS(92) =  CONSTS(91)*2.00000
      CONSTS(93) = ( 1.00000*1.00000*CONSTS(92))/( CONSTS(1)*1.00000)
      CONSTS(94) =  CONSTS(87)*0.678000
      CONSTS(95) =  CONSTS(87)*0.0552000
      CONSTS(96) =  CONSTS(87)*0.00480000
      CONSTS(97) =  CONSTS(87)*0.0200000
      CONSTS(98) =  CONSTS(87)*0.00200000
      CONSTS(99) = CONSTS(4)/(CONSTS(4)+CONSTS(22))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
      RATES(34) = (ALGBRC(7) - STATES(34))/CONSTS(43)
      ALGBRC(3) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(15) =  0.0800000*EXP(- STATES(1)/11.0000)
      RATES(27) =  ALGBRC(3)*(1.00000 - STATES(27)) -  ALGBRC(15)*STATES(27)
      ALGBRC(4) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000))
      ALGBRC(16) = 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(28) =  ALGBRC(4)*(1.00000 - STATES(28)) -  ALGBRC(16)*STATES(28)
      ALGBRC(5) = 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(17) = 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(29) =  ALGBRC(5)*(1.00000 - STATES(29)) -  ALGBRC(17)*STATES(29)
      ALGBRC(6) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(18) =  0.0800000*EXP(- STATES(1)/11.0000)
      RATES(33) =  ALGBRC(6)*(1.00000 - STATES(33)) -  ALGBRC(18)*STATES(33)
      ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
      ALGBRC(19) = 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(37) = (ALGBRC(8) - STATES(37))/ALGBRC(19)
      ALGBRC(21) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(10) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      RATES(39) =  ALGBRC(21)*(1.00000 - STATES(39)) -  ALGBRC(10)*STATES(39)
      ALGBRC(22) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
      ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      RATES(40) =  ALGBRC(22)*(1.00000 - STATES(40)) -  ALGBRC(11)*STATES(40)
      ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(12) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(27) = ALGBRC(23)/(ALGBRC(23)+ALGBRC(12))
      RATES(42) = (ALGBRC(27) - STATES(42))/CONSTS(54)
      ALGBRC(9) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
      ALGBRC(20) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
      ALGBRC(26) = 1.00000/(ALGBRC(9)+ALGBRC(20))
      ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
      RATES(38) = (ALGBRC(30) - STATES(38))/ALGBRC(26)
      ALGBRC(32) =  0.0148646*EXP(( 0.0298773*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(35) =  0.0839863*EXP(( - 0.0554611*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(2) =  - 4.00000*ALGBRC(32)*STATES(2)+ ALGBRC(35)*STATES(3)
      ALGBRC(40) =  0.00311733*EXP(( - 0.426255*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(38) =  0.0146007*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(3) =  - ( 3.00000*ALGBRC(32)+ALGBRC(35)+ALGBRC(38))*STATES(3)+ 4.00000*ALGBRC(32)*STATES(2)+ 2.00000*ALGBRC(35)*STATES(4)+ ALGBRC(40)*STATES(7)
      RATES(4) =  - ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38))*STATES(4)+ 3.00000*ALGBRC(32)*STATES(3)+ 3.00000*ALGBRC(35)*STATES(5)+ ALGBRC(40)*STATES(8)
      RATES(5) =  - (ALGBRC(32)+ 3.00000*ALGBRC(35)+ 3.00000*ALGBRC(38))*STATES(5)+ 2.00000*ALGBRC(32)*STATES(4)+ 4.00000*ALGBRC(35)*STATES(6)+ ALGBRC(40)*STATES(9)
      RATES(6) =  - ( 4.00000*ALGBRC(35)+ 4.00000*ALGBRC(38))*STATES(6)+ ALGBRC(32)*STATES(5)+ ALGBRC(40)*STATES(10)
      RATES(7) =  - ( 3.00000*ALGBRC(32)+ALGBRC(40))*STATES(7)+ ALGBRC(35)*STATES(8)+ ALGBRC(38)*STATES(3)
      RATES(8) =  - ( 2.00000*ALGBRC(32)+ALGBRC(35)+ALGBRC(38)+ALGBRC(40))*STATES(8)+ 3.00000*ALGBRC(32)*STATES(7)+ 2.00000*ALGBRC(35)*STATES(9)+ 2.00000*ALGBRC(38)*STATES(4)+ 2.00000*ALGBRC(40)*STATES(11)
      RATES(9) =  - (ALGBRC(32)+ 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38)+ALGBRC(40))*STATES(9)+ 2.00000*ALGBRC(32)*STATES(8)+ 3.00000*ALGBRC(35)*STATES(10)+ 3.00000*ALGBRC(38)*STATES(5)+ 2.00000*ALGBRC(40)*STATES(12)
      RATES(10) =  - ( 3.00000*ALGBRC(35)+ 3.00000*ALGBRC(38)+ALGBRC(40))*STATES(10)+ ALGBRC(32)*STATES(9)+ 4.00000*ALGBRC(38)*STATES(6)+ 2.00000*ALGBRC(40)*STATES(13)
      RATES(11) =  - ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(40))*STATES(11)+ ALGBRC(35)*STATES(12)+ ALGBRC(38)*STATES(8)
      RATES(12) =  - (ALGBRC(32)+ALGBRC(35)+ALGBRC(38)+ 2.00000*ALGBRC(40))*STATES(12)+ 2.00000*ALGBRC(32)*STATES(11)+ 2.00000*ALGBRC(35)*STATES(13)+ 2.00000*ALGBRC(38)*STATES(9)+ 3.00000*ALGBRC(40)*STATES(14)
      RATES(13) =  - ( 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38)+ 2.00000*ALGBRC(40))*STATES(13)+ ALGBRC(32)*STATES(12)+ 3.00000*ALGBRC(38)*STATES(10)+ 3.00000*ALGBRC(40)*STATES(15)
      RATES(14) =  - (ALGBRC(32)+ 3.00000*ALGBRC(40))*STATES(14)+ ALGBRC(35)*STATES(15)+ ALGBRC(38)*STATES(12)
      RATES(15) =  - (ALGBRC(35)+ALGBRC(38)+ 3.00000*ALGBRC(40))*STATES(15)+ ALGBRC(32)*STATES(14)+ 2.00000*ALGBRC(38)*STATES(13)+ 4.00000*ALGBRC(40)*STATES(16)
      ALGBRC(43) =  0.0773199*EXP(( - 0.0647261*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(16) =  - ( 4.00000*ALGBRC(40)+CONSTS(86))*STATES(16)+ ALGBRC(38)*STATES(15)+ ALGBRC(43)*STATES(17)
      ALGBRC(28) = 1.00000 - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17))
      ALGBRC(46) =  0.794055*EXP(( - 0.0801738*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(49) =  0.586382*EXP(( 0.282056*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(17) =  - (ALGBRC(43)+ALGBRC(49))*STATES(17)+ ALGBRC(46)*ALGBRC(28)+ CONSTS(86)*STATES(16)
      ALGBRC(14) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(25) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(29) = ALGBRC(25)/ALGBRC(14)
      ALGBRC(33) = (1.00000 - ALGBRC(25))/ALGBRC(14)
      ALGBRC(60) = 3.00000/(1.00000+0.00300000/STATES(26) ** 4.00000)
      ALGBRC(39) = 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(36) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
      ALGBRC(41) = ALGBRC(36)/ALGBRC(39)
      ALGBRC(44) = (1.00000 - ALGBRC(36))/ALGBRC(39)
      RATES(19) =  - (ALGBRC(29)+ALGBRC(60)+ALGBRC(44))*STATES(19)+ ALGBRC(33)*STATES(20)+ CONSTS(16)*STATES(21)+ ALGBRC(41)*STATES(23)
      RATES(20) =  - (ALGBRC(33)+ALGBRC(60)+ALGBRC(44))*STATES(20)+ ALGBRC(29)*STATES(19)+ CONSTS(16)*STATES(22)+ ALGBRC(41)*STATES(24)
      ALGBRC(50) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(26) ** 4.00000)
      ALGBRC(52) = 1.00000/( ALGBRC(50)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(54) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(52))
      ALGBRC(47) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(56) = ALGBRC(47)/ALGBRC(54)
      ALGBRC(58) = (1.00000 - ALGBRC(47))/ALGBRC(54)
      RATES(21) =  - (ALGBRC(29)+CONSTS(16)+ALGBRC(58))*STATES(21)+ ALGBRC(60)*STATES(19)+ ALGBRC(33)*STATES(22)+ ALGBRC(56)*STATES(25)
      ALGBRC(2) = 1.00000 - (STATES(19)+STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25))
      RATES(22) =  - (ALGBRC(33)+CONSTS(16)+ALGBRC(58))*STATES(22)+ ALGBRC(60)*STATES(20)+ ALGBRC(29)*STATES(21)+ ALGBRC(56)*ALGBRC(2)
      ALGBRC(62) = ( CONSTS(17)*ALGBRC(41)*ALGBRC(58)*ALGBRC(60))/( ALGBRC(44)*ALGBRC(56)*CONSTS(16))
      RATES(23) =  - (ALGBRC(29)+ALGBRC(62)+ALGBRC(41))*STATES(23)+ ALGBRC(44)*STATES(19)+ CONSTS(17)*STATES(25)+ ALGBRC(33)*STATES(24)
      RATES(24) =  - (ALGBRC(33)+ALGBRC(62)+ALGBRC(41))*STATES(24)+ ALGBRC(44)*STATES(20)+ CONSTS(17)*ALGBRC(2)+ ALGBRC(29)*STATES(23)
      RATES(25) =  - (ALGBRC(29)+CONSTS(17)+ALGBRC(56))*STATES(25)+ ALGBRC(62)*STATES(23)+ ALGBRC(58)*STATES(21)+ ALGBRC(33)*ALGBRC(2)
      ALGBRC(77) = (STATES(44) - STATES(43))/CONSTS(69)
      ALGBRC(80) = 1.00000/(1.00000+( CONSTS(76)*CONSTS(77))/STATES(43)+CONSTS(76) ** 2.00000)
      RATES(43) =  ALGBRC(80)*(ALGBRC(77) - STATES(41))
      ALGBRC(63) =  0.200000*CONSTS(27)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(31)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(64) = 1.00000+CONSTS(26)/STATES(31) ** 2.00000
      ALGBRC(65) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(66) =  CONSTS(33)*STATES(32) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(31)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(31)/CONSTS(32))
      ALGBRC(67) =  CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(30) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(31)
      ALGBRC(68) = ALGBRC(63)/( ALGBRC(64)*ALGBRC(65)*(ALGBRC(66)+ALGBRC(67)))
      ALGBRC(79) = (STATES(31) - STATES(18))/CONSTS(70)
      ALGBRC(83) = (STATES(31) - STATES(26))/CONSTS(71)
      ALGBRC(86) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(31) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(31) ** 2.00000)
      RATES(31) =  - ALGBRC(86)*((ALGBRC(79)+ALGBRC(83)) - (( 2.00000*ALGBRC(68)*CONSTS(93))/( 2.00000*CONSTS(97))+( STATES(41)*CONSTS(96))/CONSTS(97)))
      ALGBRC(34) = ( (( CONSTS(15)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(18)*STATES(26)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(19)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(37) =  ALGBRC(34)*(STATES(20)+STATES(22))
      ALGBRC(87) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(26) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(26) ** 2.00000)
      RATES(26) =  - ALGBRC(87)*(( ALGBRC(37)*CONSTS(93))/( 2.00000*CONSTS(98)) - ( ALGBRC(83)*CONSTS(97))/CONSTS(98))
      ALGBRC(85) = ( CONSTS(85)*(1.00000 - STATES(46)))/(1.00000+CONSTS(82)/STATES(31))+STATES(46)
      RATES(46) =  CONSTS(83)*ALGBRC(85)*(ALGBRC(85) - STATES(46)) -  CONSTS(84)*STATES(46)
      ALGBRC(72) = ( CONSTS(35)*STATES(18))/(CONSTS(36)+STATES(18))
      ALGBRC(73) = ( (( CONSTS(37)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(38)*STATES(18)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(39)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(53) =  0.800000*CONSTS(27)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(18)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(55) = 1.00000+CONSTS(26)/STATES(18) ** 2.00000
      ALGBRC(57) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(59) =  CONSTS(33)*STATES(30) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(18)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(18)/CONSTS(32))
      ALGBRC(61) =  CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(30) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(18)
      ALGBRC(69) = ALGBRC(53)/( ALGBRC(55)*ALGBRC(57)*(ALGBRC(59)+ALGBRC(61)))
      ALGBRC(70) = ALGBRC(69)+ALGBRC(68)
      ALGBRC(93) = ( CONSTS(64)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
      ALGBRC(90) = ( CONSTS(66)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
      ALGBRC(96) = ( (ALGBRC(90)+1.00000)*CONSTS(65)*STATES(18))/((STATES(18)+CONSTS(63)) - ALGBRC(93))
      ALGBRC(75) =  (0.00437500/CONSTS(68))*STATES(44)
      ALGBRC(82) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(74))/STATES(18)+CONSTS(74) ** 2.00000+( CONSTS(75)*CONSTS(73))/STATES(18)+CONSTS(75) ** 2.00000)
      RATES(18) =  ALGBRC(82)*(( - ((ALGBRC(73)+ALGBRC(72)) -  2.00000*ALGBRC(70))*CONSTS(93))/( CONSTS(94)*2.00000)+( (ALGBRC(75) - ALGBRC(96))*CONSTS(95))/CONSTS(94)+( ALGBRC(79)*CONSTS(97))/CONSTS(94))
      RATES(44) = ALGBRC(96) - (ALGBRC(75)+( ALGBRC(77)*CONSTS(96))/CONSTS(95))
      ALGBRC(89) = CONSTS(60)/(1.00000+CONSTS(61)/( 1.00000*ALGBRC(85)) ** CONSTS(59))
      ALGBRC(92) =  CONSTS(56)*(1.00000+ALGBRC(89))
      ALGBRC(95) =  ALGBRC(92)*CONSTS(57)
      ALGBRC(100) = ( ALGBRC(37)*ALGBRC(95))/(1.00000+CONSTS(55)/STATES(43) ** CONSTS(58))
      ALGBRC(98) = ALGBRC(92)/(1.00000+CONSTS(62)/STATES(43))
      RATES(41) = - (ALGBRC(100)+STATES(41))/ALGBRC(98)
      ALGBRC(81) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
      ALGBRC(84) =  CONSTS(20)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(81))
      ALGBRC(42) = ( CONSTS(24)*(STATES(1) - CONSTS(25))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
      ALGBRC(48) = STATES(30)/(STATES(30)+CONSTS(21)) ** 3.00000
      ALGBRC(51) =  CONSTS(23)*ALGBRC(45)*CONSTS(99)*ALGBRC(48)
      ALGBRC(88) =  CONSTS(42)*STATES(33) ** 3.00000*STATES(34)*(STATES(1) - ALGBRC(81))
      ALGBRC(76) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(78) = ( CONSTS(1)*CONSTS(40)*ALGBRC(76)*( STATES(30)*EXP(ALGBRC(76)) - CONSTS(6)))/(EXP(ALGBRC(76)) - 1.00000)
      ALGBRC(102) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(36))
      ALGBRC(108) = ( CONSTS(46)*ALGBRC(81) - ALGBRC(102) ** 4.00000)/(ALGBRC(81) - ALGBRC(102) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(110) = (STATES(32) - STATES(30))/CONSTS(70)
      RATES(30) = - (( ( 3.00000*ALGBRC(70)+ 3.00000*ALGBRC(51)+ALGBRC(84)+ALGBRC(88)+ALGBRC(78))*CONSTS(93))/CONSTS(94) - (ALGBRC(108)+( ALGBRC(110)*CONSTS(97))/CONSTS(94)))
      RATES(32) = - (( 3.00000*ALGBRC(68)*CONSTS(93))/CONSTS(97)+ALGBRC(110))
      ALGBRC(1) =  INT(VOI/CONSTS(11))*CONSTS(11)
      ALGBRC(13) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
      ALGBRC(94) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(35))
      ALGBRC(104) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(94)) - 59.2150)))
      ALGBRC(105) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(94))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(94)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(94))+4.75300)))
      ALGBRC(106) = ALGBRC(104)/(ALGBRC(104)+ALGBRC(105))
      ALGBRC(107) =  CONSTS(88)*ALGBRC(106)*(STATES(1) - ALGBRC(94))
      ALGBRC(71) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(97) =  CONSTS(34)*ALGBRC(71)*(STATES(1) - ALGBRC(94))
      ALGBRC(24) =  CONSTS(14)*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(18) ** 1.40000))
      ALGBRC(31) = STATES(17)+ALGBRC(28)
      ALGBRC(99) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(44)*CONSTS(6))/(STATES(35)+ CONSTS(44)*STATES(30)))
      ALGBRC(101) =  ALGBRC(24)*ALGBRC(31)*(STATES(1) - ALGBRC(99))
      ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(112) =  CONSTS(89)*STATES(37)*ALGBRC(111)*(STATES(1) - ALGBRC(94))
      ALGBRC(113) = EXP(STATES(1)/550.000)
      ALGBRC(114) =  CONSTS(49)*STATES(38) ** 3.00000*STATES(39)*STATES(40)*ALGBRC(113)*(STATES(1) - ALGBRC(94))
      ALGBRC(109) = ( CONSTS(47)*(ALGBRC(94) - ALGBRC(102)))/((ALGBRC(94)+87.8251) - ALGBRC(102))
      RATES(35) = - (( ((ALGBRC(101)+ALGBRC(112)+ALGBRC(107)+ALGBRC(114)+ALGBRC(97)+ALGBRC(13)) -  2.00000*ALGBRC(51))*CONSTS(93))/CONSTS(94) - ALGBRC(109))
      ALGBRC(74) = (ALGBRC(37)+ALGBRC(73)+ALGBRC(72)) -  2.00000*(ALGBRC(70)+ALGBRC(68))
      ALGBRC(91) = ALGBRC(84)+ 3.00000*(ALGBRC(70)+ALGBRC(68))+ 3.00000*ALGBRC(51)+ALGBRC(88)+ALGBRC(78)
      ALGBRC(115) = ALGBRC(112)+ALGBRC(101)+ALGBRC(107)+ALGBRC(97)+ - 2.00000*ALGBRC(51)+ALGBRC(114)+ALGBRC(13)
      ALGBRC(116) = ( (( CONSTS(50)*CONSTS(51) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(36) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(117) = 1.00000 - 1.00000/(1.00000+STATES(41)/CONSTS(53) ** 2.00000)
      ALGBRC(118) =  ALGBRC(116)*CONSTS(52)*STATES(42)*ALGBRC(117)
      ALGBRC(103) =  CONSTS(41)*(STATES(1) - ALGBRC(102))
      ALGBRC(119) = ALGBRC(103)+ALGBRC(118)
      RATES(1) = - (ALGBRC(91)+ALGBRC(115)+ALGBRC(74)+ALGBRC(119))
      ALGBRC(120) = (STATES(45) - STATES(36))/CONSTS(70)
      RATES(36) = - (( ALGBRC(103)*CONSTS(93))/( - 1.00000*CONSTS(94)) - (ALGBRC(108)+ALGBRC(109)+( ALGBRC(120)*CONSTS(97))/CONSTS(94)))
      RATES(45) = - (( ALGBRC(118)*CONSTS(93))/( - 1.00000*CONSTS(97))+ALGBRC(120))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
      ALGBRC(3) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(15) =  0.0800000*EXP(- STATES(1)/11.0000)
      ALGBRC(4) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.135000*EXP((80.0000+STATES(1))/- 6.80000))
      ALGBRC(16) = 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(5) = 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(17) = 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(6) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
      ALGBRC(18) =  0.0800000*EXP(- STATES(1)/11.0000)
      ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
      ALGBRC(19) = 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(21) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(10) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(22) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
      ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
      ALGBRC(12) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
      ALGBRC(27) = ALGBRC(23)/(ALGBRC(23)+ALGBRC(12))
      ALGBRC(9) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
      ALGBRC(20) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
      ALGBRC(26) = 1.00000/(ALGBRC(9)+ALGBRC(20))
      ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
      ALGBRC(32) =  0.0148646*EXP(( 0.0298773*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(35) =  0.0839863*EXP(( - 0.0554611*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(40) =  0.00311733*EXP(( - 0.426255*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(38) =  0.0146007*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(43) =  0.0773199*EXP(( - 0.0647261*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(28) = 1.00000 - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17))
      ALGBRC(46) =  0.794055*EXP(( - 0.0801738*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(49) =  0.586382*EXP(( 0.282056*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(14) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(25) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(29) = ALGBRC(25)/ALGBRC(14)
      ALGBRC(33) = (1.00000 - ALGBRC(25))/ALGBRC(14)
      ALGBRC(60) = 3.00000/(1.00000+0.00300000/STATES(26) ** 4.00000)
      ALGBRC(39) = 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(36) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
      ALGBRC(41) = ALGBRC(36)/ALGBRC(39)
      ALGBRC(44) = (1.00000 - ALGBRC(36))/ALGBRC(39)
      ALGBRC(50) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(26) ** 4.00000)
      ALGBRC(52) = 1.00000/( ALGBRC(50)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(54) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(52))
      ALGBRC(47) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(56) = ALGBRC(47)/ALGBRC(54)
      ALGBRC(58) = (1.00000 - ALGBRC(47))/ALGBRC(54)
      ALGBRC(2) = 1.00000 - (STATES(19)+STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25))
      ALGBRC(62) = ( CONSTS(17)*ALGBRC(41)*ALGBRC(58)*ALGBRC(60))/( ALGBRC(44)*ALGBRC(56)*CONSTS(16))
      ALGBRC(77) = (STATES(44) - STATES(43))/CONSTS(69)
      ALGBRC(80) = 1.00000/(1.00000+( CONSTS(76)*CONSTS(77))/STATES(43)+CONSTS(76) ** 2.00000)
      ALGBRC(63) =  0.200000*CONSTS(27)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(31)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(64) = 1.00000+CONSTS(26)/STATES(31) ** 2.00000
      ALGBRC(65) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(66) =  CONSTS(33)*STATES(32) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(31)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(31)/CONSTS(32))
      ALGBRC(67) =  CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(30) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(31)
      ALGBRC(68) = ALGBRC(63)/( ALGBRC(64)*ALGBRC(65)*(ALGBRC(66)+ALGBRC(67)))
      ALGBRC(79) = (STATES(31) - STATES(18))/CONSTS(70)
      ALGBRC(83) = (STATES(31) - STATES(26))/CONSTS(71)
      ALGBRC(86) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(31) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(31) ** 2.00000)
      ALGBRC(34) = ( (( CONSTS(15)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(18)*STATES(26)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(19)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(37) =  ALGBRC(34)*(STATES(20)+STATES(22))
      ALGBRC(87) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(26) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(26) ** 2.00000)
      ALGBRC(85) = ( CONSTS(85)*(1.00000 - STATES(46)))/(1.00000+CONSTS(82)/STATES(31))+STATES(46)
      ALGBRC(72) = ( CONSTS(35)*STATES(18))/(CONSTS(36)+STATES(18))
      ALGBRC(73) = ( (( CONSTS(37)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(38)*STATES(18)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(39)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
      ALGBRC(53) =  0.800000*CONSTS(27)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*STATES(18)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(55) = 1.00000+CONSTS(26)/STATES(18) ** 2.00000
      ALGBRC(57) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(59) =  CONSTS(33)*STATES(30) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(18)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(18)/CONSTS(32))
      ALGBRC(61) =  CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(30) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(18)
      ALGBRC(69) = ALGBRC(53)/( ALGBRC(55)*ALGBRC(57)*(ALGBRC(59)+ALGBRC(61)))
      ALGBRC(70) = ALGBRC(69)+ALGBRC(68)
      ALGBRC(93) = ( CONSTS(64)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
      ALGBRC(90) = ( CONSTS(66)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
      ALGBRC(96) = ( (ALGBRC(90)+1.00000)*CONSTS(65)*STATES(18))/((STATES(18)+CONSTS(63)) - ALGBRC(93))
      ALGBRC(75) =  (0.00437500/CONSTS(68))*STATES(44)
      ALGBRC(82) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(74))/STATES(18)+CONSTS(74) ** 2.00000+( CONSTS(75)*CONSTS(73))/STATES(18)+CONSTS(75) ** 2.00000)
      ALGBRC(89) = CONSTS(60)/(1.00000+CONSTS(61)/( 1.00000*ALGBRC(85)) ** CONSTS(59))
      ALGBRC(92) =  CONSTS(56)*(1.00000+ALGBRC(89))
      ALGBRC(95) =  ALGBRC(92)*CONSTS(57)
      ALGBRC(100) = ( ALGBRC(37)*ALGBRC(95))/(1.00000+CONSTS(55)/STATES(43) ** CONSTS(58))
      ALGBRC(98) = ALGBRC(92)/(1.00000+CONSTS(62)/STATES(43))
      ALGBRC(81) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
      ALGBRC(84) =  CONSTS(20)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(81))
      ALGBRC(42) = ( CONSTS(24)*(STATES(1) - CONSTS(25))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
      ALGBRC(48) = STATES(30)/(STATES(30)+CONSTS(21)) ** 3.00000
      ALGBRC(51) =  CONSTS(23)*ALGBRC(45)*CONSTS(99)*ALGBRC(48)
      ALGBRC(88) =  CONSTS(42)*STATES(33) ** 3.00000*STATES(34)*(STATES(1) - ALGBRC(81))
      ALGBRC(76) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(78) = ( CONSTS(1)*CONSTS(40)*ALGBRC(76)*( STATES(30)*EXP(ALGBRC(76)) - CONSTS(6)))/(EXP(ALGBRC(76)) - 1.00000)
      ALGBRC(102) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(36))
      ALGBRC(108) = ( CONSTS(46)*ALGBRC(81) - ALGBRC(102) ** 4.00000)/(ALGBRC(81) - ALGBRC(102) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(110) = (STATES(32) - STATES(30))/CONSTS(70)
      ALGBRC(1) =  INT(VOI/CONSTS(11))*CONSTS(11)
      ALGBRC(13) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
      ALGBRC(94) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(35))
      ALGBRC(104) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(94)) - 59.2150)))
      ALGBRC(105) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(94))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(94)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(94))+4.75300)))
      ALGBRC(106) = ALGBRC(104)/(ALGBRC(104)+ALGBRC(105))
      ALGBRC(107) =  CONSTS(88)*ALGBRC(106)*(STATES(1) - ALGBRC(94))
      ALGBRC(71) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(97) =  CONSTS(34)*ALGBRC(71)*(STATES(1) - ALGBRC(94))
      ALGBRC(24) =  CONSTS(14)*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(18) ** 1.40000))
      ALGBRC(31) = STATES(17)+ALGBRC(28)
      ALGBRC(99) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(44)*CONSTS(6))/(STATES(35)+ CONSTS(44)*STATES(30)))
      ALGBRC(101) =  ALGBRC(24)*ALGBRC(31)*(STATES(1) - ALGBRC(99))
      ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(112) =  CONSTS(89)*STATES(37)*ALGBRC(111)*(STATES(1) - ALGBRC(94))
      ALGBRC(113) = EXP(STATES(1)/550.000)
      ALGBRC(114) =  CONSTS(49)*STATES(38) ** 3.00000*STATES(39)*STATES(40)*ALGBRC(113)*(STATES(1) - ALGBRC(94))
      ALGBRC(109) = ( CONSTS(47)*(ALGBRC(94) - ALGBRC(102)))/((ALGBRC(94)+87.8251) - ALGBRC(102))
      ALGBRC(74) = (ALGBRC(37)+ALGBRC(73)+ALGBRC(72)) -  2.00000*(ALGBRC(70)+ALGBRC(68))
      ALGBRC(91) = ALGBRC(84)+ 3.00000*(ALGBRC(70)+ALGBRC(68))+ 3.00000*ALGBRC(51)+ALGBRC(88)+ALGBRC(78)
      ALGBRC(115) = ALGBRC(112)+ALGBRC(101)+ALGBRC(107)+ALGBRC(97)+ - 2.00000*ALGBRC(51)+ALGBRC(114)+ALGBRC(13)
      ALGBRC(116) = ( (( CONSTS(50)*CONSTS(51) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(36) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(117) = 1.00000 - 1.00000/(1.00000+STATES(41)/CONSTS(53) ** 2.00000)
      ALGBRC(118) =  ALGBRC(116)*CONSTS(52)*STATES(42)*ALGBRC(117)
      ALGBRC(103) =  CONSTS(41)*(STATES(1) - ALGBRC(102))
      ALGBRC(119) = ALGBRC(103)+ALGBRC(118)
      ALGBRC(120) = (STATES(45) - STATES(36))/CONSTS(70)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END