C
C There are a total of 128 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 94 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(89) is Vmyo in component cell_geometry (uL).
C CONSTS(90) is Vnsr in component cell_geometry (uL).
C CONSTS(85) is Vsr in component cell_geometry (uL).
C CONSTS(91) is Vjsr in component cell_geometry (uL).
C CONSTS(92) is Vss_sr in component cell_geometry (uL).
C CONSTS(93) is Vss_CaL in component cell_geometry (uL).
C CONSTS(87) is Acap in component cell_geometry (cm2).
C CONSTS(88) 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(82) is vcell in component cell_geometry (uL).
C CONSTS(86) 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(111) is caiont in component membrane (uA_per_uF).
C ALGBRC(114) is naiont in component membrane (uA_per_uF).
C ALGBRC(81) is kiont in component membrane (uA_per_uF).
C ALGBRC(72) is clont in component membrane (uA_per_uF).
C ALGBRC(52) is INa in component INa (uA_per_uF).
C ALGBRC(99) is ICaL in component ICaL (uA_per_uF).
C ALGBRC(62) is IK1 in component IK1 (uA_per_uF).
C ALGBRC(55) is IKp in component IKp (uA_per_uF).
C ALGBRC(80) is IKs in component IKs (uA_per_uF).
C ALGBRC(66) is IKr in component IKr (uA_per_uF).
C ALGBRC(87) is IpCa in component IpCa (uA_per_uF).
C ALGBRC(88) is ICab in component ICab (uA_per_uF).
C ALGBRC(86) is INaCa in component INaCa (uA_per_uF).
C ALGBRC(110) is INaCa_ss_sr in component INaCa (uA_per_uF).
C ALGBRC(45) is INaK in component INaK (uA_per_uF).
C ALGBRC(68) is Ito1 in component Ito1 (uA_per_uF).
C ALGBRC(71) is Ito2 in component Ito2 (uA_per_uF).
C ALGBRC(58) is IClb in component IClb (uA_per_uF).
C ALGBRC(53) is INaL in component INaL (uA_per_uF).
C ALGBRC(50) is INab in component INab (uA_per_uF).
C ALGBRC(79) 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(24) is OKs in component IKs (dimensionless).
C ALGBRC(2) is alpha in component IKs (per_ms).
C ALGBRC(14) is beta in component IKs (per_ms).
C ALGBRC(29) is delta in component IKs (per_ms).
C ALGBRC(25) is gamma in component IKs (per_ms).
C ALGBRC(32) is eta in component IKs (per_ms).
C CONSTS(81) is theta in component IKs (per_ms).
C ALGBRC(34) is omega in component IKs (per_ms).
C ALGBRC(36) is psi in component IKs (per_ms).
C ALGBRC(56) is EKs in component reversal_potentials (mV).
C ALGBRC(78) is Ca_i in component Ca_buffering (mM).
C ALGBRC(97) is ICaL_max in component ICaL (uA_per_uF).
C CONSTS(14) 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 STATES(26) is OI_star in component ICaL (dimensionless).
C ALGBRC(28) is ACT_tau in component ICaL (ms).
C ALGBRC(38) is IV_tau in component ICaL (ms).
C ALGBRC(100) is IV_tau_star in component ICaL (ms).
C ALGBRC(31) is ACT_infinity in component ICaL (dimensionless).
C ALGBRC(37) is IV_infinity in component ICaL (dimensionless).
C ALGBRC(41) is IV_infinity_star in component ICaL (dimensionless).
C ALGBRC(98) is IV_beta_star in component ICaL (dimensionless).
C ALGBRC(96) is IV_beta_infinity_star in component ICaL (dimensionless).
C ALGBRC(33) is alpha in component ICaL (per_ms).
C ALGBRC(35) is beta in component ICaL (per_ms).
C ALGBRC(106) is delta in component ICaL (per_ms).
C ALGBRC(39) is x in component ICaL (per_ms).
C ALGBRC(40) is y in component ICaL (per_ms).
C CONSTS(15) is theta in component ICaL (per_ms).
C ALGBRC(108) is delta_I in component ICaL (per_ms).
C ALGBRC(102) is x_star in component ICaL (per_ms).
C ALGBRC(104) 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 ALGBRC(95) is Ca_ss_CaL in component Ca_ss_buffering (mM).
C CONSTS(19) is g_Na in component INa (mS_per_uF).
C ALGBRC(51) 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(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(42) is phi in component INaK (dimensionless).
C ALGBRC(43) is fv in component INaK (dimensionless).
C CONSTS(94) is PK in component INaK (dimensionless).
C ALGBRC(44) is PNa in component INaK (dimensionless).
C STATES(30) 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 ALGBRC(101) is Ca_ss_sr in component Ca_ss_buffering (mM).
C STATES(31) is Na_ss_sr in component Na (mM).
C ALGBRC(82) is num in component INaCa (mM4_uA_per_uF).
C ALGBRC(103) is num_ss in component INaCa (mM4_uA_per_uF).
C ALGBRC(83) is denom_1 in component INaCa (dimensionless).
C ALGBRC(46) is denom_2 in component INaCa (dimensionless).
C ALGBRC(84) is denom_3 in component INaCa (mM4).
C ALGBRC(85) is denom_4 in component INaCa (mM4).
C ALGBRC(105) is denom_ss_1 in component INaCa (dimensionless).
C ALGBRC(47) is denom_ss_2 in component INaCa (dimensionless).
C ALGBRC(107) is denom_ss_3 in component INaCa (mM4).
C ALGBRC(109) is denom_ss_4 in component INaCa (mM4).
C CONSTS(33) is gKp in component IKp (mS_per_uF).
C ALGBRC(48) is Kp in component IKp (dimensionless).
C ALGBRC(54) 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(49) is phi in component INab (dimensionless).
C CONSTS(40) is gClb in component IClb (mS_per_uF).
C ALGBRC(57) is ECl in component reversal_potentials (mV).
C CONSTS(41) is gNaL in component INaL (mS_per_uF).
C STATES(32) is mL in component INaL_mL_gate (dimensionless).
C STATES(33) 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(42) is tau_hL in component INaL_hL_gate (ms).
C ALGBRC(7) is hL_infinity in component INaL_hL_gate (dimensionless).
C CONSTS(43) is prnak in component reversal_potentials (dimensionless).
C STATES(34) is K_i in component K (mM).
C STATES(35) is Cl_i in component Cl (mM).
C CONSTS(83) is gK1 in component IK1 (mS_per_uF).
C ALGBRC(61) is K1 in component IK1_K1_gate (dimensionless).
C ALGBRC(59) is alpha_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(60) is beta_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(63) 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(64) 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(84) is gKr in component IKr (mS_per_uF).
C ALGBRC(65) is r in component IKr_r_gate (dimensionless).
C STATES(36) 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(46) is gto1 in component Ito1 (mS_per_uF).
C ALGBRC(67) is rto1 in component Ito1 (dimensionless).
C STATES(37) is a in component Ito1_a_gate (dimensionless).
C STATES(38) is i1f in component Ito1_i1f_gate (dimensionless).
C STATES(39) 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(47) is PCl in component Ito2 (L_per_F_ms).
C CONSTS(48) is zCl in component Ito2 (dimensionless).
C ALGBRC(69) is Ito2_max in component Ito2 (uA_per_uF).
C ALGBRC(70) is KCa_ito2 in component Ito2 (dimensionless).
C CONSTS(49) is kCa_ito2 in component Ito2 (mM_per_ms).
C STATES(40) is Irel in component Irel (mM_per_ms).
C STATES(41) 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(50) is tau_i2f in component Ito2_i2f_gate (ms).
C ALGBRC(128) is rel_infinity in component Irel (mM_per_ms).
C ALGBRC(127) is tau_rel in component Irel (ms).
C CONSTS(51) is Krel_infinity in component Irel (mM).
C ALGBRC(125) is alpha_rel in component Irel (mM_per_uA_per_uF_per_ms).
C ALGBRC(123) 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(121) 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 ALGBRC(118) is Ca_JSR in component Ca_JSR_buffering (mM).
C ALGBRC(120) is CaMK_active in component CaMK_active (dimensionless).
C ALGBRC(126) is Iup in component Iup (mM_per_ms).
C CONSTS(59) is kmup in component Iup (mM).
C ALGBRC(124) 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(122) 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 CaMK_active (mM).
C ALGBRC(74) is Ileak in component Ileak (mM_per_ms).
C CONSTS(64) is nsrbar in component Ileak (mM).
C STATES(42) is Ca_NSR in component Ca (mM).
C ALGBRC(119) is Itr in component Itr (mM_per_ms).
C CONSTS(65) is tautr in component Itr (ms).
C STATES(43) is Ca_JSR_0 in component Ca (mM).
C STATES(44) is Ca_i_0 in component Ca (mM).
C STATES(45) is Ca_ss_sr_0 in component Ca (mM).
C STATES(46) is Ca_ss_CaL_0 in component Ca (mM).
C CONSTS(66) is tau_diff_ss in component Ca (ms).
C CONSTS(67) is tau_diff in component Ca (ms).
C ALGBRC(112) is Idiff in component Ca (mM_per_ms).
C ALGBRC(116) is Idiff_ss in component Ca (mM_per_ms).
C ALGBRC(75) is b_myo in component Ca_buffering (mM).
C ALGBRC(76) is c_myo in component Ca_buffering (mM2).
C ALGBRC(77) is d_myo in component Ca_buffering (mM3).
C CONSTS(68) is cmdn_bar in component Ca_buffering (mM).
C CONSTS(69) is trpn_bar in component Ca_buffering (mM).
C CONSTS(70) is km_cmdn in component Ca_buffering (mM).
C CONSTS(71) is km_trpn in component Ca_buffering (mM).
C ALGBRC(89) is b_myo1 in component Ca_ss_buffering (mM).
C ALGBRC(90) is c_myo1 in component Ca_ss_buffering (mM2).
C ALGBRC(91) is d_myo1 in component Ca_ss_buffering (mM3).
C ALGBRC(92) is b_myo2 in component Ca_ss_buffering (mM).
C ALGBRC(93) is c_myo2 in component Ca_ss_buffering (mM2).
C ALGBRC(94) is d_myo2 in component Ca_ss_buffering (mM3).
C CONSTS(72) is BSL_max in component Ca_ss_buffering (mM).
C CONSTS(73) is BSR_max in component Ca_ss_buffering (mM).
C CONSTS(74) is Km_BSR in component Ca_ss_buffering (mM).
C CONSTS(75) is Km_BSL in component Ca_ss_buffering (mM).
C ALGBRC(115) is b in component Ca_JSR_buffering (mM).
C ALGBRC(117) is c in component Ca_JSR_buffering (mM2).
C CONSTS(76) is csqn_bar in component Ca_JSR_buffering (mM).
C CONSTS(77) is km_csqn in component Ca_JSR_buffering (mM).
C ALGBRC(113) is Idiff_Na in component Na (mM_per_ms).
C STATES(47) is Cl_ss in component Cl (mM).
C ALGBRC(73) is Idiff_Cl in component Cl (mM_per_ms).
C STATES(48) is CaMK_trap in component CaMK_active (dimensionless).
C CONSTS(78) is alpha_CaMK in component CaMK_active (per_ms).
C CONSTS(79) is beta_CaMK in component CaMK_active (per_ms).
C CONSTS(80) 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(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(26) is d/dt OI_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(32) is d/dt mL in component INaL_mL_gate (dimensionless).
C RATES(33) is d/dt hL in component INaL_hL_gate (dimensionless).
C RATES(36) is d/dt xr in component IKr_xr_gate (dimensionless).
C RATES(37) is d/dt a in component Ito1_a_gate (dimensionless).
C RATES(38) is d/dt i1f in component Ito1_i1f_gate (dimensionless).
C RATES(39) is d/dt i1s in component Ito1_i1s_gate (dimensionless).
C RATES(41) is d/dt i2f in component Ito2_i2f_gate (dimensionless).
C RATES(40) is d/dt Irel in component Irel (mM_per_ms).
C RATES(44) is d/dt Ca_i_0 in component Ca (mM).
C RATES(45) is d/dt Ca_ss_sr_0 in component Ca (mM).
C RATES(46) is d/dt Ca_ss_CaL_0 in component Ca (mM).
C RATES(42) is d/dt Ca_NSR in component Ca (mM).
C RATES(43) is d/dt Ca_JSR_0 in component Ca (mM).
C RATES(30) is d/dt Na_i in component Na (mM).
C RATES(31) is d/dt Na_ss_sr in component Na (mM).
C RATES(35) is d/dt Cl_i in component Cl (mM).
C RATES(47) is d/dt Cl_ss in component Cl (mM).
C RATES(34) 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
      CONSTS(14) = 1.5552e-4
      STATES(19) = 0.9999613782488198
      STATES(20) = 2.27419494597377e-5
      STATES(21) = 5.936477396599124e-11
      STATES(22) = 1.350122566369778e-15
      STATES(23) = 1.587938121237188e-5
      STATES(24) = 3.611420243201921e-10
      STATES(25) = 0
      STATES(26) = 0
      CONSTS(15) = 1
      CONSTS(16) = 0.000001
      CONSTS(17) = 1
      CONSTS(18) = 0.341
      CONSTS(19) = 9.075
      STATES(27) = 0.001022252813581797
      STATES(28) = 0.9913109066867
      STATES(29) = 0.9970771239754529
      CONSTS(20) = 2.6
      CONSTS(21) = 1.5
      CONSTS(22) = 1.4
      CONSTS(23) = 1
      CONSTS(24) = -92
      STATES(30) = 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(31) = 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(32) = 0.001022252813581797
      STATES(33) = 0.3593255076983574
      CONSTS(42) = 600
      CONSTS(43) = 0.01833
      STATES(34) = 145.5689
      STATES(35) = 20.26842
      CONSTS(44) = 2.46108e-5
      CONSTS(45) = 1.77e-5
      STATES(36) = 1.235996442041028e-8
      CONSTS(46) = 0.4975
      STATES(37) = 1.773318268499115e-5
      STATES(38) = 0.9979722
      STATES(39) = 0.9874435
      CONSTS(47) = 9e-7
      CONSTS(48) = -1
      CONSTS(49) = 0.4
      STATES(40) = 0.00107332
      STATES(41) = 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
      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(42) = 1.159939
      CONSTS(65) = 25
      STATES(43) = 6.99469048
      STATES(44) = 0.013408
      STATES(45) = 0.02414815
      STATES(46) = 0.02510283
      CONSTS(66) = 2
      CONSTS(67) = 0.2
      CONSTS(68) = 0.05
      CONSTS(69) = 0.07
      CONSTS(70) = 0.00238
      CONSTS(71) = 0.0005
      CONSTS(72) = 1.124
      CONSTS(73) = 0.047
      CONSTS(74) = 0.00087
      CONSTS(75) = 0.0087
      CONSTS(76) = 10
      CONSTS(77) = 0.8
      STATES(47) = 20.26842
      STATES(48) = 0.0018861341
      CONSTS(78) = 0.05
      CONSTS(79) = 0.00068
      CONSTS(80) = 0.05
      CONSTS(81) = 0.0895380
      CONSTS(82) =  1000.00* 3.14159265358979*CONSTS(9)*CONSTS(9)*CONSTS(8)
      CONSTS(83) =  0.500000* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(84) =  0.0138542* (CONSTS(4)/5.40000) ** (1.0 / 2)
      CONSTS(85) =  CONSTS(82)*0.0600000
      CONSTS(86) =  2.00000* 3.14159265358979*CONSTS(9)*CONSTS(9)+ 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(8)
      CONSTS(87) =  CONSTS(86)*2.00000
      CONSTS(88) = ( 1.00000*1.00000*CONSTS(87))/( CONSTS(1)*1.00000)
      CONSTS(89) =  CONSTS(82)*0.678000
      CONSTS(90) =  CONSTS(82)*0.0552000
      CONSTS(91) =  CONSTS(82)*0.00480000
      CONSTS(92) =  CONSTS(82)*0.0200000
      CONSTS(93) =  CONSTS(82)*0.00200000
      CONSTS(94) = CONSTS(4)/(CONSTS(4)+CONSTS(21))
      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(33) = (ALGBRC(7) - STATES(33))/CONSTS(42)
      ALGBRC(2) =  0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(14) =  0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(2) =  - 4.00000*ALGBRC(2)*STATES(2)+ ALGBRC(14)*STATES(3)
      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(32) =  ALGBRC(6)*(1.00000 - STATES(32)) -  ALGBRC(18)*STATES(32)
      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(36) = (ALGBRC(8) - STATES(36))/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(38) =  ALGBRC(21)*(1.00000 - STATES(38)) -  ALGBRC(10)*STATES(38)
      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(39) =  ALGBRC(22)*(1.00000 - STATES(39)) -  ALGBRC(11)*STATES(39)
      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(41) = (ALGBRC(27) - STATES(41))/CONSTS(50)
      ALGBRC(29) =  0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(25) =  0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(3) =  - ( 3.00000*ALGBRC(2)+ALGBRC(14)+ALGBRC(25))*STATES(3)+ 4.00000*ALGBRC(2)*STATES(2)+ 2.00000*ALGBRC(14)*STATES(4)+ ALGBRC(29)*STATES(7)
      RATES(4) =  - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25))*STATES(4)+ 3.00000*ALGBRC(2)*STATES(3)+ 3.00000*ALGBRC(14)*STATES(5)+ ALGBRC(29)*STATES(8)
      RATES(5) =  - (ALGBRC(2)+ 3.00000*ALGBRC(14)+ 3.00000*ALGBRC(25))*STATES(5)+ 2.00000*ALGBRC(2)*STATES(4)+ 4.00000*ALGBRC(14)*STATES(6)+ ALGBRC(29)*STATES(9)
      RATES(6) =  - ( 4.00000*ALGBRC(14)+ 4.00000*ALGBRC(25))*STATES(6)+ ALGBRC(2)*STATES(5)+ ALGBRC(29)*STATES(10)
      RATES(7) =  - ( 3.00000*ALGBRC(2)+ALGBRC(29))*STATES(7)+ ALGBRC(14)*STATES(8)+ ALGBRC(25)*STATES(3)
      RATES(8) =  - ( 2.00000*ALGBRC(2)+ALGBRC(14)+ALGBRC(25)+ALGBRC(29))*STATES(8)+ 3.00000*ALGBRC(2)*STATES(7)+ 2.00000*ALGBRC(14)*STATES(9)+ 2.00000*ALGBRC(25)*STATES(4)+ 2.00000*ALGBRC(29)*STATES(11)
      RATES(9) =  - (ALGBRC(2)+ 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25)+ALGBRC(29))*STATES(9)+ 2.00000*ALGBRC(2)*STATES(8)+ 3.00000*ALGBRC(14)*STATES(10)+ 3.00000*ALGBRC(25)*STATES(5)+ 2.00000*ALGBRC(29)*STATES(12)
      RATES(10) =  - ( 3.00000*ALGBRC(14)+ 3.00000*ALGBRC(25)+ALGBRC(29))*STATES(10)+ ALGBRC(2)*STATES(9)+ 4.00000*ALGBRC(25)*STATES(6)+ 2.00000*ALGBRC(29)*STATES(13)
      RATES(11) =  - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(29))*STATES(11)+ ALGBRC(14)*STATES(12)+ ALGBRC(25)*STATES(8)
      RATES(12) =  - (ALGBRC(2)+ALGBRC(14)+ALGBRC(25)+ 2.00000*ALGBRC(29))*STATES(12)+ 2.00000*ALGBRC(2)*STATES(11)+ 2.00000*ALGBRC(14)*STATES(13)+ 2.00000*ALGBRC(25)*STATES(9)+ 3.00000*ALGBRC(29)*STATES(14)
      RATES(13) =  - ( 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25)+ 2.00000*ALGBRC(29))*STATES(13)+ ALGBRC(2)*STATES(12)+ 3.00000*ALGBRC(25)*STATES(10)+ 3.00000*ALGBRC(29)*STATES(15)
      RATES(14) =  - (ALGBRC(2)+ 3.00000*ALGBRC(29))*STATES(14)+ ALGBRC(14)*STATES(15)+ ALGBRC(25)*STATES(12)
      RATES(15) =  - (ALGBRC(14)+ALGBRC(25)+ 3.00000*ALGBRC(29))*STATES(15)+ ALGBRC(2)*STATES(14)+ 2.00000*ALGBRC(25)*STATES(13)+ 4.00000*ALGBRC(29)*STATES(16)
      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(37) = (ALGBRC(30) - STATES(37))/ALGBRC(26)
      ALGBRC(32) =  0.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(16) =  - ( 4.00000*ALGBRC(29)+CONSTS(81))*STATES(16)+ ALGBRC(25)*STATES(15)+ ALGBRC(32)*STATES(17)
      ALGBRC(34) =  0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(36) =  0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      RATES(17) =  - (ALGBRC(32)+ALGBRC(36))*STATES(17)+ ALGBRC(34)*STATES(18)+ CONSTS(81)*STATES(16)
      RATES(18) =  - ALGBRC(34)*STATES(18)+ ALGBRC(36)*STATES(17)
      ALGBRC(57) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(35))
      ALGBRC(58) =  CONSTS(40)*(STATES(1) - ALGBRC(57))
      ALGBRC(51) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
      ALGBRC(63) = ( CONSTS(44)*ALGBRC(51) - ALGBRC(57) ** 4.00000)/(ALGBRC(51) - ALGBRC(57) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(54) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(34))
      ALGBRC(64) = ( CONSTS(45)*(ALGBRC(54) - ALGBRC(57)))/((ALGBRC(54)+87.8251) - ALGBRC(57))
      ALGBRC(73) = (STATES(47) - STATES(35))/CONSTS(67)
      RATES(35) = - (( ALGBRC(58)*CONSTS(88))/( - 1.00000*CONSTS(89)) - (ALGBRC(63)+ALGBRC(64)+( ALGBRC(73)*CONSTS(92))/CONSTS(89)))
      ALGBRC(69) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(35) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(70) = 1.00000 - 1.00000/(1.00000+STATES(40)/CONSTS(49) ** 2.00000)
      ALGBRC(71) =  ALGBRC(69)*STATES(41)*ALGBRC(70)
      RATES(47) = - (( ALGBRC(71)*CONSTS(88))/( - 1.00000*CONSTS(92))+ALGBRC(73))
      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(59) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(54)) - 59.2150)))
      ALGBRC(60) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(54))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(54)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(54))+4.75300)))
      ALGBRC(61) = ALGBRC(59)/(ALGBRC(59)+ALGBRC(60))
      ALGBRC(62) =  CONSTS(83)*ALGBRC(61)*(STATES(1) - ALGBRC(54))
      ALGBRC(48) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(55) =  CONSTS(33)*ALGBRC(48)*(STATES(1) - ALGBRC(54))
      ALGBRC(75) = ((CONSTS(68)+CONSTS(69)) - STATES(44))+CONSTS(71)+CONSTS(70)
      ALGBRC(76) = ( CONSTS(71)*CONSTS(70) -  STATES(44)*(CONSTS(71)+CONSTS(70)))+ CONSTS(69)*CONSTS(70)+ CONSTS(68)*CONSTS(71)
      ALGBRC(77) =  - CONSTS(71)*CONSTS(70)*STATES(44)
      ALGBRC(78) =  (( 2.00000* (ALGBRC(75) ** 2.00000 -  3.00000*ALGBRC(76)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(75)*ALGBRC(76) -  2.00000*ALGBRC(75) ** 3.00000) -  27.0000*ALGBRC(77))/( 2.00000*ALGBRC(75) ** 2.00000 -  3.00000*ALGBRC(76) ** 1.50000))/3.00000) - ALGBRC(75)/3.00000
      ALGBRC(79) =  0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(78) ** 1.40000))
      ALGBRC(24) = STATES(17)+STATES(18)
      ALGBRC(56) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(34)+ CONSTS(43)*STATES(30)))
      ALGBRC(80) =  ALGBRC(79)*ALGBRC(24)*(STATES(1) - ALGBRC(56))
      ALGBRC(65) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(66) =  CONSTS(84)*STATES(36)*ALGBRC(65)*(STATES(1) - ALGBRC(54))
      ALGBRC(42) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(43) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
      ALGBRC(44) = STATES(30)/(STATES(30)+CONSTS(20)) ** 3.00000
      ALGBRC(45) =  CONSTS(22)*ALGBRC(43)*CONSTS(94)*ALGBRC(44)
      ALGBRC(67) = EXP(STATES(1)/550.000)
      ALGBRC(68) =  CONSTS(46)*STATES(37) ** 3.00000*STATES(38)*STATES(39)*ALGBRC(67)*(STATES(1) - ALGBRC(54))
      RATES(34) = - (( ((ALGBRC(80)+ALGBRC(66)+ALGBRC(62)+ALGBRC(68)+ALGBRC(55)+ALGBRC(13)) -  2.00000*ALGBRC(45))*CONSTS(88))/CONSTS(89) - ALGBRC(64))
      ALGBRC(28) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(33) = ALGBRC(31)/ALGBRC(28)
      ALGBRC(35) = (1.00000 - ALGBRC(31))/ALGBRC(28)
      ALGBRC(89) = ((CONSTS(72)+CONSTS(73)) - STATES(46))+CONSTS(74)+CONSTS(75)
      ALGBRC(90) = ( CONSTS(75)*CONSTS(74) -  STATES(46)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
      ALGBRC(91) =  - CONSTS(74)*CONSTS(75)*STATES(46)
      ALGBRC(95) =  (( 2.00000* (ALGBRC(89) ** 2.00000 -  3.00000*ALGBRC(90)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(89)*ALGBRC(90) -  2.00000*ALGBRC(89) ** 3.00000) -  27.0000*ALGBRC(91))/( 2.00000*ALGBRC(89) ** 2.00000 -  3.00000*ALGBRC(90) ** 1.50000))/3.00000) - ALGBRC(89)/3.00000
      ALGBRC(106) = 3.00000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
      ALGBRC(38) = 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(39) = ALGBRC(37)/ALGBRC(38)
      ALGBRC(40) = (1.00000 - ALGBRC(37))/ALGBRC(38)
      RATES(19) =  - (ALGBRC(33)+ALGBRC(106)+ALGBRC(40))*STATES(19)+ ALGBRC(35)*STATES(20)+ CONSTS(15)*STATES(21)+ ALGBRC(39)*STATES(23)
      RATES(20) =  - (ALGBRC(35)+ALGBRC(106)+ALGBRC(40))*STATES(20)+ ALGBRC(33)*STATES(19)+ CONSTS(15)*STATES(22)+ ALGBRC(39)*STATES(24)
      ALGBRC(96) = 25.0000 - 17.5000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
      ALGBRC(98) = 1.00000/( ALGBRC(96)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(100) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(98))
      ALGBRC(41) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(102) = ALGBRC(41)/ALGBRC(100)
      ALGBRC(104) = (1.00000 - ALGBRC(41))/ALGBRC(100)
      RATES(21) =  - (ALGBRC(33)+CONSTS(15)+ALGBRC(104))*STATES(21)+ ALGBRC(106)*STATES(19)+ ALGBRC(35)*STATES(22)+ ALGBRC(102)*STATES(25)
      RATES(22) =  - (ALGBRC(35)+CONSTS(15)+ALGBRC(104))*STATES(22)+ ALGBRC(106)*STATES(20)+ ALGBRC(33)*STATES(21)+ ALGBRC(102)*STATES(26)
      ALGBRC(108) = ( CONSTS(16)*ALGBRC(39)*ALGBRC(104)*ALGBRC(106))/( ALGBRC(40)*ALGBRC(102)*CONSTS(15))
      RATES(23) =  - (ALGBRC(33)+ALGBRC(108)+ALGBRC(39))*STATES(23)+ ALGBRC(40)*STATES(19)+ CONSTS(16)*STATES(25)+ ALGBRC(35)*STATES(24)
      RATES(24) =  - (ALGBRC(35)+ALGBRC(108)+ALGBRC(39))*STATES(24)+ ALGBRC(40)*STATES(20)+ CONSTS(16)*STATES(26)+ ALGBRC(33)*STATES(23)
      RATES(25) =  - (ALGBRC(33)+CONSTS(16)+ALGBRC(102))*STATES(25)+ ALGBRC(108)*STATES(23)+ ALGBRC(104)*STATES(21)+ ALGBRC(35)*STATES(26)
      RATES(26) =  - (ALGBRC(35)+CONSTS(16)+ALGBRC(102))*STATES(26)+ ALGBRC(108)*STATES(24)+ ALGBRC(104)*STATES(22)+ ALGBRC(33)*STATES(25)
      ALGBRC(52) =  CONSTS(19)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(51))
      ALGBRC(82) =  0.800000*CONSTS(26)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*ALGBRC(78)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(83) = 1.00000+CONSTS(25)/ALGBRC(78) ** 2.00000
      ALGBRC(46) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(84) =  CONSTS(32)*STATES(30) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(78)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(78)/CONSTS(31))
      ALGBRC(85) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(29) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(78)
      ALGBRC(86) = ALGBRC(82)/( ALGBRC(83)*ALGBRC(46)*(ALGBRC(84)+ALGBRC(85)))
      ALGBRC(53) =  CONSTS(41)*STATES(32) ** 3.00000*STATES(33)*(STATES(1) - ALGBRC(51))
      ALGBRC(49) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(50) = ( CONSTS(1)*CONSTS(39)*ALGBRC(49)*( STATES(30)*EXP(ALGBRC(49)) - CONSTS(6)))/(EXP(ALGBRC(49)) - 1.00000)
      ALGBRC(113) = (STATES(31) - STATES(30))/CONSTS(67)
      RATES(30) = - (( ( 3.00000*ALGBRC(86)+ 3.00000*ALGBRC(45)+ALGBRC(52)+ALGBRC(53)+ALGBRC(50))*CONSTS(88))/CONSTS(89) - (ALGBRC(63)+( ALGBRC(113)*CONSTS(92))/CONSTS(89)))
      ALGBRC(92) = ((CONSTS(72)+CONSTS(73)) - STATES(45))+CONSTS(74)+CONSTS(75)
      ALGBRC(93) = ( CONSTS(75)*CONSTS(74) -  STATES(45)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
      ALGBRC(94) =  - CONSTS(74)*CONSTS(75)*STATES(45)
      ALGBRC(101) =  (( 2.00000* (ALGBRC(92) ** 2.00000 -  3.00000*ALGBRC(93)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(92)*ALGBRC(93) -  2.00000*ALGBRC(92) ** 3.00000) -  27.0000*ALGBRC(94))/( 2.00000*ALGBRC(92) ** 2.00000 -  3.00000*ALGBRC(93) ** 1.50000))/3.00000) - ALGBRC(92)/3.00000
      ALGBRC(103) =  0.200000*CONSTS(26)*( STATES(31) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*ALGBRC(101)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(105) = 1.00000+CONSTS(25)/ALGBRC(101) ** 2.00000
      ALGBRC(47) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(107) =  CONSTS(32)*STATES(31) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(101)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(101)/CONSTS(31))
      ALGBRC(109) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(31)/CONSTS(29) ** 3.00000)+ STATES(31) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(101)
      ALGBRC(110) = ALGBRC(103)/( ALGBRC(105)*ALGBRC(47)*(ALGBRC(107)+ALGBRC(109)))
      RATES(31) = - (( 3.00000*ALGBRC(110)*CONSTS(88))/CONSTS(92)+ALGBRC(113))
      ALGBRC(97) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*ALGBRC(95)*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(99) =  ALGBRC(97)*(STATES(20)+STATES(22))
      ALGBRC(87) = ( CONSTS(34)*ALGBRC(78))/(CONSTS(35)+ALGBRC(78))
      ALGBRC(88) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*ALGBRC(78)*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(111) = (ALGBRC(99)+ALGBRC(88)+ALGBRC(87)) -  2.00000*(ALGBRC(86)+ALGBRC(110))
      ALGBRC(114) = ALGBRC(52)+ 3.00000*(ALGBRC(86)+ALGBRC(110))+ 3.00000*ALGBRC(45)+ALGBRC(53)+ALGBRC(50)
      ALGBRC(81) = ALGBRC(66)+ALGBRC(80)+ALGBRC(62)+ALGBRC(55)+ - 2.00000*ALGBRC(45)+ALGBRC(68)+ALGBRC(13)
      ALGBRC(72) = ALGBRC(58)+ALGBRC(71)
      RATES(1) = - (ALGBRC(114)+ALGBRC(81)+ALGBRC(111)+ALGBRC(72))
      ALGBRC(112) = (ALGBRC(101) - ALGBRC(78))/CONSTS(67)
      ALGBRC(116) = (ALGBRC(101) - ALGBRC(95))/CONSTS(66)
      RATES(45) = - ((ALGBRC(112)+ALGBRC(116)) - (( 2.00000*ALGBRC(110)*CONSTS(88))/( 2.00000*CONSTS(92))+( STATES(40)*CONSTS(91))/CONSTS(92)))
      RATES(46) = - (( ALGBRC(99)*CONSTS(88))/( 2.00000*CONSTS(93)) - ( ALGBRC(116)*CONSTS(92))/CONSTS(93))
      ALGBRC(115) = (CONSTS(76)+CONSTS(77)) - STATES(43)
      ALGBRC(117) =  STATES(43)*CONSTS(77)
      ALGBRC(118) = - ALGBRC(115)/2.00000+ (ALGBRC(115) ** 2.00000+ 4.00000*ALGBRC(117)) ** (1.0 / 2)/2.00000
      ALGBRC(119) = (STATES(42) - ALGBRC(118))/CONSTS(65)
      RATES(43) = ALGBRC(119) - STATES(40)
      ALGBRC(120) = ( CONSTS(80)*(1.00000 - STATES(48))*1.00000)/(1.00000+CONSTS(63)/ALGBRC(101))
      RATES(48) =  CONSTS(78)*ALGBRC(120)*(ALGBRC(120) - STATES(48)) -  CONSTS(79)*STATES(48)
      ALGBRC(124) = ( CONSTS(60)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
      ALGBRC(122) = ( CONSTS(62)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
      ALGBRC(126) = ( (ALGBRC(122)+1.00000)*CONSTS(61)*ALGBRC(78))/((ALGBRC(78)+CONSTS(59)) - ALGBRC(124))
      ALGBRC(74) =  (0.00437500/CONSTS(64))*STATES(42)
      RATES(44) = (( - ((ALGBRC(88)+ALGBRC(87)) -  2.00000*ALGBRC(86))*CONSTS(88))/( 2.00000*CONSTS(89))+( (ALGBRC(126) - ALGBRC(74))*CONSTS(90))/CONSTS(89)) - ( ALGBRC(112)*CONSTS(92))/CONSTS(89)
      RATES(42) = ALGBRC(126) - (ALGBRC(74)+( ALGBRC(119)*CONSTS(91))/CONSTS(90))
      ALGBRC(121) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(120)) ** CONSTS(55))
      ALGBRC(123) =  CONSTS(52)*(1.00000+ALGBRC(121))
      ALGBRC(125) =  ALGBRC(123)*CONSTS(53)
      ALGBRC(128) = ( ALGBRC(99)*ALGBRC(125))/(1.00000+CONSTS(51)/ALGBRC(118) ** CONSTS(54))
      ALGBRC(127) = ALGBRC(123)/(1.00000+CONSTS(58)/ALGBRC(118))
      RATES(40) = - (ALGBRC(128)+STATES(40))/ALGBRC(127)
      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(2) =  0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(14) =  0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      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(29) =  0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(25) =  0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      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.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(34) =  0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(36) =  0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(57) =  (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(35))
      ALGBRC(58) =  CONSTS(40)*(STATES(1) - ALGBRC(57))
      ALGBRC(51) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
      ALGBRC(63) = ( CONSTS(44)*ALGBRC(51) - ALGBRC(57) ** 4.00000)/(ALGBRC(51) - ALGBRC(57) ** 4.00000+87.8251 ** 4.00000)
      ALGBRC(54) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(34))
      ALGBRC(64) = ( CONSTS(45)*(ALGBRC(54) - ALGBRC(57)))/((ALGBRC(54)+87.8251) - ALGBRC(57))
      ALGBRC(73) = (STATES(47) - STATES(35))/CONSTS(67)
      ALGBRC(69) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(35) -  CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(70) = 1.00000 - 1.00000/(1.00000+STATES(40)/CONSTS(49) ** 2.00000)
      ALGBRC(71) =  ALGBRC(69)*STATES(41)*ALGBRC(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(59) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(54)) - 59.2150)))
      ALGBRC(60) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(54))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(54)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(54))+4.75300)))
      ALGBRC(61) = ALGBRC(59)/(ALGBRC(59)+ALGBRC(60))
      ALGBRC(62) =  CONSTS(83)*ALGBRC(61)*(STATES(1) - ALGBRC(54))
      ALGBRC(48) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
      ALGBRC(55) =  CONSTS(33)*ALGBRC(48)*(STATES(1) - ALGBRC(54))
      ALGBRC(75) = ((CONSTS(68)+CONSTS(69)) - STATES(44))+CONSTS(71)+CONSTS(70)
      ALGBRC(76) = ( CONSTS(71)*CONSTS(70) -  STATES(44)*(CONSTS(71)+CONSTS(70)))+ CONSTS(69)*CONSTS(70)+ CONSTS(68)*CONSTS(71)
      ALGBRC(77) =  - CONSTS(71)*CONSTS(70)*STATES(44)
      ALGBRC(78) =  (( 2.00000* (ALGBRC(75) ** 2.00000 -  3.00000*ALGBRC(76)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(75)*ALGBRC(76) -  2.00000*ALGBRC(75) ** 3.00000) -  27.0000*ALGBRC(77))/( 2.00000*ALGBRC(75) ** 2.00000 -  3.00000*ALGBRC(76) ** 1.50000))/3.00000) - ALGBRC(75)/3.00000
      ALGBRC(79) =  0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(78) ** 1.40000))
      ALGBRC(24) = STATES(17)+STATES(18)
      ALGBRC(56) =  (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(34)+ CONSTS(43)*STATES(30)))
      ALGBRC(80) =  ALGBRC(79)*ALGBRC(24)*(STATES(1) - ALGBRC(56))
      ALGBRC(65) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
      ALGBRC(66) =  CONSTS(84)*STATES(36)*ALGBRC(65)*(STATES(1) - ALGBRC(54))
      ALGBRC(42) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(43) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
      ALGBRC(44) = STATES(30)/(STATES(30)+CONSTS(20)) ** 3.00000
      ALGBRC(45) =  CONSTS(22)*ALGBRC(43)*CONSTS(94)*ALGBRC(44)
      ALGBRC(67) = EXP(STATES(1)/550.000)
      ALGBRC(68) =  CONSTS(46)*STATES(37) ** 3.00000*STATES(38)*STATES(39)*ALGBRC(67)*(STATES(1) - ALGBRC(54))
      ALGBRC(28) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
      ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
      ALGBRC(33) = ALGBRC(31)/ALGBRC(28)
      ALGBRC(35) = (1.00000 - ALGBRC(31))/ALGBRC(28)
      ALGBRC(89) = ((CONSTS(72)+CONSTS(73)) - STATES(46))+CONSTS(74)+CONSTS(75)
      ALGBRC(90) = ( CONSTS(75)*CONSTS(74) -  STATES(46)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
      ALGBRC(91) =  - CONSTS(74)*CONSTS(75)*STATES(46)
      ALGBRC(95) =  (( 2.00000* (ALGBRC(89) ** 2.00000 -  3.00000*ALGBRC(90)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(89)*ALGBRC(90) -  2.00000*ALGBRC(89) ** 3.00000) -  27.0000*ALGBRC(91))/( 2.00000*ALGBRC(89) ** 2.00000 -  3.00000*ALGBRC(90) ** 1.50000))/3.00000) - ALGBRC(89)/3.00000
      ALGBRC(106) = 3.00000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
      ALGBRC(38) = 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(39) = ALGBRC(37)/ALGBRC(38)
      ALGBRC(40) = (1.00000 - ALGBRC(37))/ALGBRC(38)
      ALGBRC(96) = 25.0000 - 17.5000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
      ALGBRC(98) = 1.00000/( ALGBRC(96)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
      ALGBRC(100) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(98))
      ALGBRC(41) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
      ALGBRC(102) = ALGBRC(41)/ALGBRC(100)
      ALGBRC(104) = (1.00000 - ALGBRC(41))/ALGBRC(100)
      ALGBRC(108) = ( CONSTS(16)*ALGBRC(39)*ALGBRC(104)*ALGBRC(106))/( ALGBRC(40)*ALGBRC(102)*CONSTS(15))
      ALGBRC(52) =  CONSTS(19)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(51))
      ALGBRC(82) =  0.800000*CONSTS(26)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*ALGBRC(78)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(83) = 1.00000+CONSTS(25)/ALGBRC(78) ** 2.00000
      ALGBRC(46) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(84) =  CONSTS(32)*STATES(30) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(78)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(78)/CONSTS(31))
      ALGBRC(85) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(29) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(78)
      ALGBRC(86) = ALGBRC(82)/( ALGBRC(83)*ALGBRC(46)*(ALGBRC(84)+ALGBRC(85)))
      ALGBRC(53) =  CONSTS(41)*STATES(32) ** 3.00000*STATES(33)*(STATES(1) - ALGBRC(51))
      ALGBRC(49) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
      ALGBRC(50) = ( CONSTS(1)*CONSTS(39)*ALGBRC(49)*( STATES(30)*EXP(ALGBRC(49)) - CONSTS(6)))/(EXP(ALGBRC(49)) - 1.00000)
      ALGBRC(113) = (STATES(31) - STATES(30))/CONSTS(67)
      ALGBRC(92) = ((CONSTS(72)+CONSTS(73)) - STATES(45))+CONSTS(74)+CONSTS(75)
      ALGBRC(93) = ( CONSTS(75)*CONSTS(74) -  STATES(45)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
      ALGBRC(94) =  - CONSTS(74)*CONSTS(75)*STATES(45)
      ALGBRC(101) =  (( 2.00000* (ALGBRC(92) ** 2.00000 -  3.00000*ALGBRC(93)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(92)*ALGBRC(93) -  2.00000*ALGBRC(92) ** 3.00000) -  27.0000*ALGBRC(94))/( 2.00000*ALGBRC(92) ** 2.00000 -  3.00000*ALGBRC(93) ** 1.50000))/3.00000) - ALGBRC(92)/3.00000
      ALGBRC(103) =  0.200000*CONSTS(26)*( STATES(31) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) -  CONSTS(6) ** 3.00000*ALGBRC(101)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
      ALGBRC(105) = 1.00000+CONSTS(25)/ALGBRC(101) ** 2.00000
      ALGBRC(47) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
      ALGBRC(107) =  CONSTS(32)*STATES(31) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(101)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(101)/CONSTS(31))
      ALGBRC(109) =  CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(31)/CONSTS(29) ** 3.00000)+ STATES(31) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(101)
      ALGBRC(110) = ALGBRC(103)/( ALGBRC(105)*ALGBRC(47)*(ALGBRC(107)+ALGBRC(109)))
      ALGBRC(97) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*ALGBRC(95)*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(99) =  ALGBRC(97)*(STATES(20)+STATES(22))
      ALGBRC(87) = ( CONSTS(34)*ALGBRC(78))/(CONSTS(35)+ALGBRC(78))
      ALGBRC(88) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*ALGBRC(78)*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(111) = (ALGBRC(99)+ALGBRC(88)+ALGBRC(87)) -  2.00000*(ALGBRC(86)+ALGBRC(110))
      ALGBRC(114) = ALGBRC(52)+ 3.00000*(ALGBRC(86)+ALGBRC(110))+ 3.00000*ALGBRC(45)+ALGBRC(53)+ALGBRC(50)
      ALGBRC(81) = ALGBRC(66)+ALGBRC(80)+ALGBRC(62)+ALGBRC(55)+ - 2.00000*ALGBRC(45)+ALGBRC(68)+ALGBRC(13)
      ALGBRC(72) = ALGBRC(58)+ALGBRC(71)
      ALGBRC(112) = (ALGBRC(101) - ALGBRC(78))/CONSTS(67)
      ALGBRC(116) = (ALGBRC(101) - ALGBRC(95))/CONSTS(66)
      ALGBRC(115) = (CONSTS(76)+CONSTS(77)) - STATES(43)
      ALGBRC(117) =  STATES(43)*CONSTS(77)
      ALGBRC(118) = - ALGBRC(115)/2.00000+ (ALGBRC(115) ** 2.00000+ 4.00000*ALGBRC(117)) ** (1.0 / 2)/2.00000
      ALGBRC(119) = (STATES(42) - ALGBRC(118))/CONSTS(65)
      ALGBRC(120) = ( CONSTS(80)*(1.00000 - STATES(48))*1.00000)/(1.00000+CONSTS(63)/ALGBRC(101))
      ALGBRC(124) = ( CONSTS(60)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
      ALGBRC(122) = ( CONSTS(62)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
      ALGBRC(126) = ( (ALGBRC(122)+1.00000)*CONSTS(61)*ALGBRC(78))/((ALGBRC(78)+CONSTS(59)) - ALGBRC(124))
      ALGBRC(74) =  (0.00437500/CONSTS(64))*STATES(42)
      ALGBRC(121) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(120)) ** CONSTS(55))
      ALGBRC(123) =  CONSTS(52)*(1.00000+ALGBRC(121))
      ALGBRC(125) =  ALGBRC(123)*CONSTS(53)
      ALGBRC(128) = ( ALGBRC(99)*ALGBRC(125))/(1.00000+CONSTS(51)/ALGBRC(118) ** CONSTS(54))
      ALGBRC(127) = ALGBRC(123)/(1.00000+CONSTS(58)/ALGBRC(118))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END