C
C There are a total of 93 entries in the algebraic variable array.
C There are a total of 41 entries in each of the rate and state variable arrays.
C There are a total of 148 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component model_parameters (kelvin).
C CONSTS(3) is F in component model_parameters (coulomb_per_mole).
C CONSTS(4) is Na_e in component model_parameters (millimolar).
C CONSTS(5) is Ca_e in component model_parameters (millimolar).
C CONSTS(6) is K_e in component model_parameters (millimolar).
C CONSTS(146) is Vt in component model_parameters (cm3).
C CONSTS(113) is Vd in component model_parameters (cm3).
C CONSTS(114) is Vmyo in component model_parameters (cm3).
C CONSTS(117) is Sms in component model_parameters (cm2).
C CONSTS(118) is Smt in component model_parameters (cm2).
C CONSTS(148) is Cms in component model_parameters (microF).
C CONSTS(147) is Cmt in component model_parameters (microF).
C CONSTS(107) is Vc in component model_parameters (cm3).
C CONSTS(115) is VSRup in component model_parameters (cm3).
C CONSTS(116) is VSRrel in component model_parameters (cm3).
C CONSTS(119) is pt in component model_parameters (dimensionless).
C CONSTS(144) is Lt in component model_parameters (cm).
C CONSTS(145) is Rst in component model_parameters (ohm).
C CONSTS(7) is Rot in component model_parameters (ohm_per_cm2).
C CONSTS(8) is radiust in component model_parameters (cm).
C CONSTS(9) is ptcm in component model_parameters (per_cm).
C CONSTS(10) is Smtot in component model_parameters (cm2).
C ALGBRC(2) is i_circ in component common_current (microA).
C CONSTS(11) is R_st in component common_current (kilo_ohm).
C STATES(1) is Vm_s in component membrane (millivolt).
C STATES(2) is Vm_t in component membrane (millivolt).
C ALGBRC(27) is J_CaSRrel in component SR_Ca_release_channel (millimolar_per_second).
C CONSTS(12) is kap in component SR_Ca_release_channel (per_millimolar4_per_second).
C CONSTS(13) is kam in component SR_Ca_release_channel (per_second).
C CONSTS(14) is kbp in component SR_Ca_release_channel (per_millimolar3_per_second).
C CONSTS(15) is kbm in component SR_Ca_release_channel (per_second).
C CONSTS(16) is kcp in component SR_Ca_release_channel (per_second).
C CONSTS(17) is kcm in component SR_Ca_release_channel (per_second).
C STATES(3) is F2 in component SR_Ca_release_channel (dimensionless).
C STATES(4) is F3 in component SR_Ca_release_channel (dimensionless).
C STATES(5) is F1 in component SR_Ca_release_channel (dimensionless).
C STATES(6) is F4 in component SR_Ca_release_channel (dimensionless).
C STATES(7) is Ca_ss in component ion_concentrations (millimolar).
C STATES(8) is CaSRrel in component ion_concentrations (millimolar).
C ALGBRC(50) is J_CaSRup in component SERCA2a_pump (millimolar_per_second).
C CONSTS(18) is K_fb in component SERCA2a_pump (millimolar).
C CONSTS(19) is K_rb in component SERCA2a_pump (millimolar).
C ALGBRC(48) is fb in component SERCA2a_pump (dimensionless).
C ALGBRC(49) is rb in component SERCA2a_pump (dimensionless).
C CONSTS(20) is Vmaxf in component SERCA2a_pump (millimolar_per_second).
C CONSTS(21) is Vmaxr in component SERCA2a_pump (millimolar_per_second).
C CONSTS(22) is K_SR in component SERCA2a_pump (dimensionless).
C CONSTS(23) is N_fb in component SERCA2a_pump (dimensionless).
C CONSTS(24) is N_rb in component SERCA2a_pump (dimensionless).
C STATES(9) is Ca_i in component ion_concentrations (millimolar).
C STATES(10) is CaSRup in component ion_concentrations (millimolar).
C ALGBRC(51) is J_tr in component intracellular_and_SR_Ca_fluxes (millimolar_per_second).
C ALGBRC(1) is J_xfer in component intracellular_and_SR_Ca_fluxes (millimolar_per_second).
C ALGBRC(54) is J_trpn in component intracellular_and_SR_Ca_fluxes (millimolar_per_second).
C ALGBRC(52) is J_BTRH in component intracellular_and_SR_Ca_fluxes (millimolar_per_second).
C ALGBRC(53) is J_BTRL in component intracellular_and_SR_Ca_fluxes (millimolar_per_second).
C CONSTS(25) is BTRHmax in component intracellular_and_SR_Ca_fluxes (millimolar).
C CONSTS(26) is BTRLmax in component intracellular_and_SR_Ca_fluxes (millimolar).
C CONSTS(27) is tau_tr in component intracellular_and_SR_Ca_fluxes (second).
C CONSTS(28) is tau_xfer in component intracellular_and_SR_Ca_fluxes (second).
C STATES(11) is BTRH in component intracellular_and_SR_Ca_fluxes (millimolar).
C STATES(12) is BTRL in component intracellular_and_SR_Ca_fluxes (millimolar).
C CONSTS(29) is k_htrpn_plus in component intracellular_and_SR_Ca_fluxes (per_millimolar_per_second).
C CONSTS(30) is k_htrpn_minus in component intracellular_and_SR_Ca_fluxes (per_second).
C CONSTS(31) is k_ltrpn_plus in component intracellular_and_SR_Ca_fluxes (per_millimolar_per_second).
C CONSTS(32) is k_ltrpn_minus in component intracellular_and_SR_Ca_fluxes (per_second).
C STATES(13) is Na_i in component ion_concentrations (millimolar).
C STATES(14) is K_i in component ion_concentrations (millimolar).
C STATES(15) is Na_t in component ion_concentrations (millimolar).
C STATES(16) is Ca_t in component ion_concentrations (millimolar).
C STATES(17) is K_t in component ion_concentrations (millimolar).
C ALGBRC(55) is i_Kext in component ion_concentrations (microA).
C CONSTS(33) is K_mCMDN in component ion_concentrations (millimolar).
C CONSTS(34) is K_mCSQN in component ion_concentrations (millimolar).
C CONSTS(35) is K_mEGTA in component ion_concentrations (millimolar).
C CONSTS(36) is CMDN_tot in component ion_concentrations (millimolar).
C CONSTS(37) is CSQN_tot in component ion_concentrations (millimolar).
C CONSTS(38) is EGTA_tot in component ion_concentrations (millimolar).
C CONSTS(39) is tau_d in component ion_concentrations (second).
C CONSTS(40) is STIMULATION in component ion_concentrations (dimensionless).
C CONSTS(41) is time0 in component ion_concentrations (second).
C ALGBRC(60) is i_Na_s in component sodium_current (microA).
C ALGBRC(62) is i_Ca_s in component calcium_current (microA).
C ALGBRC(69) is i_BNa_s in component background_currents (microA).
C ALGBRC(75) is i_NaCa_s in component Na_Ca_ion_exchanger_current (microA).
C ALGBRC(73) is i_NaK_s in component sodium_potassium_pump (microA).
C ALGBRC(61) is i_fNa_s in component hyperpolarisation_activated_current (microA).
C ALGBRC(67) is i_fK_s in component hyperpolarisation_activated_current (microA).
C ALGBRC(71) is i_BK_s in component background_currents (microA).
C ALGBRC(66) is i_K1_s in component inward_rectifier (microA).
C ALGBRC(64) is i_Kto_s in component Ca_independent_transient_outward_K_current (microA).
C ALGBRC(65) is i_Kss_s in component steady_state_outward_K_current (microA).
C ALGBRC(74) is i_CaP_s in component sarcolemmal_calcium_pump_current (microA).
C ALGBRC(70) is i_BCa_s in component background_currents (microA).
C ALGBRC(77) is i_Na_t in component sodium_current (microA).
C ALGBRC(79) is i_Ca_t in component calcium_current (microA).
C ALGBRC(86) is i_BNa_t in component background_currents (microA).
C ALGBRC(93) is i_NaCa_t in component Na_Ca_ion_exchanger_current (microA).
C ALGBRC(91) is i_NaK_t in component sodium_potassium_pump (microA).
C ALGBRC(78) is i_fNa_t in component hyperpolarisation_activated_current (microA).
C ALGBRC(84) is i_fK_t in component hyperpolarisation_activated_current (microA).
C ALGBRC(88) is i_BK_t in component background_currents (microA).
C ALGBRC(83) is i_K1_t in component inward_rectifier (microA).
C ALGBRC(81) is i_Kto_t in component Ca_independent_transient_outward_K_current (microA).
C ALGBRC(82) is i_Kss_t in component steady_state_outward_K_current (microA).
C ALGBRC(92) is i_CaP_t in component sarcolemmal_calcium_pump_current (microA).
C ALGBRC(87) is i_BCa_t in component background_currents (microA).
C ALGBRC(56) is Jte_Na in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(57) is Jte_Ca in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(58) is Jte_K in component t_tubular_ion_fluxes (millimolar_per_second).
C CONSTS(42) is tau_Na in component t_tubular_ion_fluxes (second).
C CONSTS(43) is tau_Ca in component t_tubular_ion_fluxes (second).
C CONSTS(44) is tau_K in component t_tubular_ion_fluxes (second).
C ALGBRC(68) is i_f_s in component hyperpolarisation_activated_current (microA).
C ALGBRC(72) is i_B_s in component background_currents (microA).
C ALGBRC(59) is E_Na in component sodium_current (millivolt).
C CONSTS(45) is g_Na in component sodium_current (milliS_per_cm2).
C CONSTS(120) is g_Nas in component sodium_current (milliS).
C CONSTS(46) is fNat in component sodium_current (dimensionless).
C STATES(18) is m in component sodium_current_m_gate (dimensionless).
C STATES(19) is h in component sodium_current_h_gate (dimensionless).
C ALGBRC(3) is alpha_m in component sodium_current_m_gate (per_second).
C ALGBRC(28) is beta_m in component sodium_current_m_gate (per_second).
C ALGBRC(4) is alpha_h in component sodium_current_h_gate (per_second).
C ALGBRC(29) is beta_h in component sodium_current_h_gate (per_second).
C STATES(20) is Co in component calcium_current (dimensionless).
C STATES(21) is C_cast in component calcium_current (dimensionless).
C STATES(22) is C1 in component calcium_current (dimensionless).
C STATES(23) is C2 in component calcium_current (dimensionless).
C ALGBRC(30) is tauCo in component calcium_current (second).
C CONSTS(47) is tauC_cast in component calcium_current (second).
C ALGBRC(31) is tauC1 in component calcium_current (second).
C ALGBRC(32) is tauC2 in component calcium_current (second).
C ALGBRC(5) is Co_inf in component calcium_current (dimensionless).
C ALGBRC(6) is C_cast_inf in component calcium_current (dimensionless).
C ALGBRC(7) is C1_inf in component calcium_current (dimensionless).
C ALGBRC(8) is C2_inf in component calcium_current (dimensionless).
C CONSTS(48) is g_Ca in component calcium_current (cm_per_second).
C CONSTS(121) is g_Cas in component calcium_current (cm3_per_second).
C CONSTS(49) is fCat in component calcium_current (dimensionless).
C ALGBRC(63) is E_K in component Ca_independent_transient_outward_K_current (millivolt).
C CONSTS(50) is a in component Ca_independent_transient_outward_K_current (dimensionless).
C CONSTS(51) is b in component Ca_independent_transient_outward_K_current (dimensionless).
C CONSTS(52) is g_Kto in component Ca_independent_transient_outward_K_current (milliS_per_cm2).
C CONSTS(122) is g_Ktos in component Ca_independent_transient_outward_K_current (milliS).
C CONSTS(53) is fKtot in component Ca_independent_transient_outward_K_current (dimensionless).
C STATES(24) is r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C STATES(25) is s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C STATES(26) is s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C ALGBRC(33) is tau_r in component Ca_independent_transient_outward_K_current_r_gate (second).
C ALGBRC(9) is r_max in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C ALGBRC(34) is tau_s in component Ca_independent_transient_outward_K_current_s_gate (second).
C ALGBRC(10) is s_max in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C ALGBRC(35) is tau_s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (second).
C ALGBRC(11) is s_slow_max in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C CONSTS(54) is g_Kss in component steady_state_outward_K_current (milliS_per_cm2).
C CONSTS(123) is g_Ksss in component steady_state_outward_K_current (milliS).
C CONSTS(55) is fKsst in component steady_state_outward_K_current (dimensionless).
C STATES(27) is r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C STATES(28) is s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C ALGBRC(36) is tau_r_ss in component steady_state_outward_K_current_r_ss_gate (second).
C ALGBRC(12) is r_ss_max in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C CONSTS(108) is tau_s_ss in component steady_state_outward_K_current_s_ss_gate (second).
C ALGBRC(13) is s_ss_max in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C CONSTS(56) is g_K1 in component inward_rectifier (milliS_per_cm2).
C CONSTS(124) is g_K1s in component inward_rectifier (milliS).
C CONSTS(57) is fK1t in component inward_rectifier (dimensionless).
C CONSTS(58) is g_f in component hyperpolarisation_activated_current (milliS_per_cm2).
C CONSTS(125) is g_fs in component hyperpolarisation_activated_current (milliS).
C CONSTS(59) is fft in component hyperpolarisation_activated_current (dimensionless).
C CONSTS(60) is f_Na in component hyperpolarisation_activated_current (dimensionless).
C CONSTS(109) is f_K in component hyperpolarisation_activated_current (dimensionless).
C STATES(29) is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(37) is tau_y in component hyperpolarisation_activated_current_y_gate (second).
C ALGBRC(14) is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
C CONSTS(61) is g_B_Ca in component background_currents (milliS_per_cm2).
C CONSTS(62) is g_B_Na in component background_currents (milliS_per_cm2).
C CONSTS(63) is g_B_K in component background_currents (milliS_per_cm2).
C CONSTS(126) is g_B_Cas in component background_currents (milliS).
C CONSTS(127) is g_B_Nas in component background_currents (milliS).
C CONSTS(128) is g_B_Ks in component background_currents (milliS).
C CONSTS(64) is fCabt in component background_currents (dimensionless).
C CONSTS(65) is fNabt in component background_currents (dimensionless).
C CONSTS(66) is fKbt in component background_currents (dimensionless).
C CONSTS(67) is E_Ca in component background_currents (millivolt).
C CONSTS(68) is i_NaK_max in component sodium_potassium_pump (microA_per_cm2).
C CONSTS(129) is i_NaK_max_s in component sodium_potassium_pump (microA).
C CONSTS(69) is fNaKt in component sodium_potassium_pump (dimensionless).
C CONSTS(70) is K_m_K in component sodium_potassium_pump (millimolar).
C CONSTS(71) is K_m_Na in component sodium_potassium_pump (millimolar).
C CONSTS(110) is sigma in component sodium_potassium_pump (dimensionless).
C CONSTS(72) is i_CaP_max in component sarcolemmal_calcium_pump_current (microA_per_cm2).
C CONSTS(130) is i_CaP_max_s in component sarcolemmal_calcium_pump_current (microA).
C CONSTS(73) is fCaPt in component sarcolemmal_calcium_pump_current (dimensionless).
C CONSTS(74) is i_NaCa_max in component Na_Ca_ion_exchanger_current (microA_per_cm2).
C CONSTS(131) is i_NaCa_max_s in component Na_Ca_ion_exchanger_current (microA).
C CONSTS(75) is fNaCat in component Na_Ca_ion_exchanger_current (dimensionless).
C ALGBRC(85) is i_f_t in component hyperpolarisation_activated_current (microA).
C ALGBRC(89) is i_B_t in component background_currents (microA).
C ALGBRC(76) is E_Na in component sodium_current (millivolt).
C CONSTS(76) is g_Na in component sodium_current (milliS_per_cm2).
C CONSTS(132) is g_Nat in component sodium_current (milliS).
C CONSTS(77) is fNat in component sodium_current (dimensionless).
C STATES(30) is m in component sodium_current_m_gate (dimensionless).
C STATES(31) is h in component sodium_current_h_gate (dimensionless).
C ALGBRC(15) is alpha_m in component sodium_current_m_gate (per_second).
C ALGBRC(38) is beta_m in component sodium_current_m_gate (per_second).
C ALGBRC(16) is alpha_h in component sodium_current_h_gate (per_second).
C ALGBRC(39) is beta_h in component sodium_current_h_gate (per_second).
C STATES(32) is Co in component calcium_current (dimensionless).
C STATES(33) is C_cast in component calcium_current (dimensionless).
C STATES(34) is C1 in component calcium_current (dimensionless).
C STATES(35) is C2 in component calcium_current (dimensionless).
C ALGBRC(40) is tauCo in component calcium_current (second).
C CONSTS(78) is tauC_cast in component calcium_current (second).
C ALGBRC(41) is tauC1 in component calcium_current (second).
C ALGBRC(42) is tauC2 in component calcium_current (second).
C ALGBRC(17) is Co_inf in component calcium_current (dimensionless).
C ALGBRC(18) is C_cast_inf in component calcium_current (dimensionless).
C ALGBRC(19) is C1_inf in component calcium_current (dimensionless).
C ALGBRC(20) is C2_inf in component calcium_current (dimensionless).
C CONSTS(79) is g_Ca in component calcium_current (cm_per_second).
C CONSTS(133) is g_Cat in component calcium_current (cm3_per_second).
C CONSTS(80) is fCat in component calcium_current (dimensionless).
C ALGBRC(80) is E_K in component Ca_independent_transient_outward_K_current (millivolt).
C CONSTS(81) is a in component Ca_independent_transient_outward_K_current (dimensionless).
C CONSTS(82) is b in component Ca_independent_transient_outward_K_current (dimensionless).
C CONSTS(83) is g_Kto in component Ca_independent_transient_outward_K_current (milliS_per_cm2).
C CONSTS(134) is g_Ktot in component Ca_independent_transient_outward_K_current (milliS).
C CONSTS(84) is fKtot in component Ca_independent_transient_outward_K_current (dimensionless).
C STATES(36) is r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C STATES(37) is s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C STATES(38) is s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C ALGBRC(43) is tau_r in component Ca_independent_transient_outward_K_current_r_gate (second).
C ALGBRC(21) is r_max in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C ALGBRC(44) is tau_t in component Ca_independent_transient_outward_K_current_s_gate (second).
C ALGBRC(22) is s_max in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C ALGBRC(45) is tau_s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (second).
C ALGBRC(23) is s_slow_max in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C CONSTS(85) is g_Kss in component steady_state_outward_K_current (milliS_per_cm2).
C CONSTS(135) is g_Ksst in component steady_state_outward_K_current (milliS).
C CONSTS(86) is fKsst in component steady_state_outward_K_current (dimensionless).
C STATES(39) is r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C STATES(40) is s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C ALGBRC(46) is tau_r_ss in component steady_state_outward_K_current_r_ss_gate (second).
C ALGBRC(24) is r_ss_max in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C CONSTS(111) is tau_s_ss in component steady_state_outward_K_current_s_ss_gate (second).
C ALGBRC(25) is s_ss_max in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C CONSTS(87) is g_K1 in component inward_rectifier (milliS_per_cm2).
C CONSTS(136) is g_K1t in component inward_rectifier (milliS).
C CONSTS(88) is fK1t in component inward_rectifier (dimensionless).
C CONSTS(89) is g_f in component hyperpolarisation_activated_current (milliS_per_cm2).
C CONSTS(137) is g_ft in component hyperpolarisation_activated_current (milliS).
C CONSTS(90) is fft in component hyperpolarisation_activated_current (dimensionless).
C CONSTS(91) is f_Na in component hyperpolarisation_activated_current (dimensionless).
C CONSTS(112) is f_K in component hyperpolarisation_activated_current (dimensionless).
C STATES(41) is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(47) is tau_y in component hyperpolarisation_activated_current_y_gate (second).
C ALGBRC(26) is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
C CONSTS(92) is g_B_Ca in component background_currents (milliS_per_cm2).
C CONSTS(93) is g_B_Na in component background_currents (milliS_per_cm2).
C CONSTS(94) is g_B_K in component background_currents (milliS_per_cm2).
C CONSTS(138) is g_B_Cat in component background_currents (milliS).
C CONSTS(139) is g_B_Nat in component background_currents (milliS).
C CONSTS(140) is g_B_Kt in component background_currents (milliS).
C CONSTS(95) is fCabt in component background_currents (dimensionless).
C CONSTS(96) is fNabt in component background_currents (dimensionless).
C CONSTS(97) is fKbt in component background_currents (dimensionless).
C CONSTS(98) is E_Ca in component background_currents (millivolt).
C CONSTS(99) is i_NaK_max in component sodium_potassium_pump (microA_per_cm2).
C CONSTS(141) is i_NaK_max_t in component sodium_potassium_pump (microA).
C CONSTS(100) is fNaKt in component sodium_potassium_pump (dimensionless).
C CONSTS(101) is K_m_K in component sodium_potassium_pump (millimolar).
C CONSTS(102) is K_m_Na in component sodium_potassium_pump (millimolar).
C ALGBRC(90) is sigma in component sodium_potassium_pump (dimensionless).
C CONSTS(103) is i_CaP_max in component sarcolemmal_calcium_pump_current (microA_per_cm2).
C CONSTS(142) is i_CaP_max_t in component sarcolemmal_calcium_pump_current (microA).
C CONSTS(104) is fCaPt in component sarcolemmal_calcium_pump_current (dimensionless).
C CONSTS(105) is i_NaCa_max in component Na_Ca_ion_exchanger_current (microA_per_cm2).
C CONSTS(143) is i_NaCa_max_t in component Na_Ca_ion_exchanger_current (microA).
C CONSTS(106) is fNaCat in component Na_Ca_ion_exchanger_current (dimensionless).
C RATES(5) is d/dt F1 in component SR_Ca_release_channel (dimensionless).
C RATES(3) is d/dt F2 in component SR_Ca_release_channel (dimensionless).
C RATES(4) is d/dt F3 in component SR_Ca_release_channel (dimensionless).
C RATES(6) is d/dt F4 in component SR_Ca_release_channel (dimensionless).
C RATES(11) is d/dt BTRH in component intracellular_and_SR_Ca_fluxes (millimolar).
C RATES(12) is d/dt BTRL in component intracellular_and_SR_Ca_fluxes (millimolar).
C RATES(13) is d/dt Na_i in component ion_concentrations (millimolar).
C RATES(15) is d/dt Na_t in component ion_concentrations (millimolar).
C RATES(14) is d/dt K_i in component ion_concentrations (millimolar).
C RATES(17) is d/dt K_t in component ion_concentrations (millimolar).
C RATES(9) is d/dt Ca_i in component ion_concentrations (millimolar).
C RATES(16) is d/dt Ca_t in component ion_concentrations (millimolar).
C RATES(7) is d/dt Ca_ss in component ion_concentrations (millimolar).
C RATES(8) is d/dt CaSRrel in component ion_concentrations (millimolar).
C RATES(10) is d/dt CaSRup in component ion_concentrations (millimolar).
C RATES(1) is d/dt Vm_s in component membrane (millivolt).
C RATES(18) is d/dt m in component sodium_current_m_gate (dimensionless).
C RATES(19) is d/dt h in component sodium_current_h_gate (dimensionless).
C RATES(20) is d/dt Co in component calcium_current (dimensionless).
C RATES(22) is d/dt C1 in component calcium_current (dimensionless).
C RATES(23) is d/dt C2 in component calcium_current (dimensionless).
C RATES(21) is d/dt C_cast in component calcium_current (dimensionless).
C RATES(24) is d/dt r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C RATES(25) is d/dt s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C RATES(26) is d/dt s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C RATES(27) is d/dt r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C RATES(28) is d/dt s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C RATES(29) is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C RATES(2) is d/dt Vm_t in component membrane (millivolt).
C RATES(30) is d/dt m in component sodium_current_m_gate (dimensionless).
C RATES(31) is d/dt h in component sodium_current_h_gate (dimensionless).
C RATES(32) is d/dt Co in component calcium_current (dimensionless).
C RATES(34) is d/dt C1 in component calcium_current (dimensionless).
C RATES(35) is d/dt C2 in component calcium_current (dimensionless).
C RATES(33) is d/dt C_cast in component calcium_current (dimensionless).
C RATES(36) is d/dt r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C RATES(37) is d/dt s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C RATES(38) is d/dt s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C RATES(39) is d/dt r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C RATES(40) is d/dt s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C RATES(41) is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 8310
      CONSTS(2) = 295
      CONSTS(3) = 96500
      CONSTS(4) = 140
      CONSTS(5) = 1.2
      CONSTS(6) = 5.4
      CONSTS(7) = 83.33
      CONSTS(8) = 150e-7
      CONSTS(9) = 30000000
      CONSTS(10) = 0.0001
      CONSTS(11) = 20.0
      STATES(1) = -7.9709547e1
      STATES(2) = -7.9709536e1
      CONSTS(12) = 12.15e12
      CONSTS(13) = 576
      CONSTS(14) = 4.05e9
      CONSTS(15) = 1930
      CONSTS(16) = 100
      CONSTS(17) = 0.8
      STATES(3) = 3.8333331e-4
      STATES(4) = 6.4516879e-11
      STATES(5) = 6.7606483e-1
      STATES(6) = 3.2355462e-1
      STATES(7) = 4.3112114e-5
      STATES(8) = 2.2638169e-1
      CONSTS(18) = 0.000168
      CONSTS(19) = 3.29
      CONSTS(20) = 0.4
      CONSTS(21) = 0.9
      CONSTS(22) = 1
      CONSTS(23) = 1.2
      CONSTS(24) = 1
      STATES(9) = 3.5486498e-5
      STATES(10) = 2.2755058e-1
      CONSTS(25) = 0.14
      CONSTS(26) = 0.07
      CONSTS(27) = 0.02
      CONSTS(28) = 0.008
      STATES(11) = 9.2804169e-1
      STATES(12) = 3.4317829e-2
      CONSTS(29) = 20000
      CONSTS(30) = 0.07
      CONSTS(31) = 40000
      CONSTS(32) = 40
      STATES(13) = 6.9435814e0
      STATES(14) = 1.3753311e2
      STATES(15) = 1.4001386e2
      STATES(16) = 1.2271399e0
      STATES(17) = 5.3627676e0
      CONSTS(33) = 0.00238
      CONSTS(34) = 0.8
      CONSTS(35) = 0.00015
      CONSTS(36) = 0.05
      CONSTS(37) = 15
      CONSTS(38) = 10
      CONSTS(39) = 0.008
      CONSTS(40) = 1
      CONSTS(41) = 0
      CONSTS(42) = 0.15
      CONSTS(43) = 0.5
      CONSTS(44) = 0.15
      CONSTS(45) = 10.0
      CONSTS(46) = 0.56
      STATES(18) = 2.1832558e-6
      STATES(19) = 5.8232781e-1
      STATES(20) = 2.5436465e-6
      STATES(21) = 9.9570063e-1
      STATES(22) = 9.9994546e-1
      STATES(23) = 9.9994546e-1
      CONSTS(47) = 0.009
      CONSTS(48) = 0.00021
      CONSTS(49) = 0.87
      CONSTS(50) = 0.886
      CONSTS(51) = 0.114
      CONSTS(52) = 0.35
      CONSTS(53) = 0.56
      STATES(24) = 2.3485072e-3
      STATES(25) = 9.8315843e-1
      STATES(26) = 6.7483370e-1
      CONSTS(54) = 0.07
      CONSTS(55) = 0.76
      STATES(27) = 3.1079213e-3
      STATES(28) = 3.0773393e-1
      CONSTS(56) = 0.24
      CONSTS(57) = 0.56
      CONSTS(58) = 0.0145
      CONSTS(59) = 0.56
      CONSTS(60) = 0.2
      STATES(29) = 3.4357356e-3
      CONSTS(61) = 0.000648
      CONSTS(62) = 0.0008015
      CONSTS(63) = 0.00138
      CONSTS(64) = 0.56
      CONSTS(65) = 0.56
      CONSTS(66) = 0.56
      CONSTS(67) = 65
      CONSTS(68) = 1
      CONSTS(69) = 0.59
      CONSTS(70) = 1.5
      CONSTS(71) = 10
      CONSTS(72) = 0.85
      CONSTS(73) = 0.56
      CONSTS(74) = 1.8e-4
      CONSTS(75) = 0.81
      CONSTS(76) = 10.0
      CONSTS(77) = 0.56
      STATES(30) = 2.1832694e-6
      STATES(31) = 5.8232717e-1
      STATES(32) = 2.5436520e-6
      STATES(33) = 9.9570063e-1
      STATES(34) = 9.9785393e-1
      STATES(35) = 9.9994546e-1
      CONSTS(78) = 0.009
      CONSTS(79) = 0.00021
      CONSTS(80) = 0.87
      CONSTS(81) = 0.886
      CONSTS(82) = 0.114
      CONSTS(83) = 0.35
      CONSTS(84) = 0.56
      STATES(36) = 2.3485094e-3
      STATES(37) = 9.8315865e-1
      STATES(38) = 6.7478532e-1
      CONSTS(85) = 0.07
      CONSTS(86) = 0.76
      STATES(39) = 3.1079241e-3
      STATES(40) = 3.0773368e-1
      CONSTS(87) = 0.24
      CONSTS(88) = 0.56
      CONSTS(89) = 0.0145
      CONSTS(90) = 0.56
      CONSTS(91) = 0.2
      STATES(41) = 3.4357175e-3
      CONSTS(92) = 0.000648
      CONSTS(93) = 0.0008015
      CONSTS(94) = 0.00138
      CONSTS(95) = 0.56
      CONSTS(96) = 0.56
      CONSTS(97) = 0.56
      CONSTS(98) = 65
      CONSTS(99) = 1
      CONSTS(100) = 0.59
      CONSTS(101) = 1.5
      CONSTS(102) = 10
      CONSTS(103) = 0.85
      CONSTS(104) = 0.56
      CONSTS(105) = 1.8e-4
      CONSTS(106) = 0.81
      CONSTS(107) = CONSTS(10)/7850.00
      CONSTS(108) = 2.10000
      CONSTS(109) = 1.00000 - CONSTS(60)
      CONSTS(110) = (EXP(CONSTS(4)/67.3000) - 1.00000)/7.00000
      CONSTS(111) = 2.10000
      CONSTS(112) = 1.00000 - CONSTS(91)
      CONSTS(113) =  CONSTS(107)*7.50000e-05
      CONSTS(114) =  CONSTS(107)*0.585000
      CONSTS(115) =  CONSTS(107)*0.0315000
      CONSTS(116) =  CONSTS(107)*0.00350000
      CONSTS(117) = CONSTS(10) -  4400.00*CONSTS(107)
      CONSTS(118) = CONSTS(10) - CONSTS(117)
      CONSTS(119) =  CONSTS(9)*CONSTS(117)
      CONSTS(120) =  CONSTS(45)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(46))
      CONSTS(121) =  CONSTS(48)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(49))
      CONSTS(122) =  CONSTS(52)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(53))
      CONSTS(123) =  CONSTS(54)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(55))
      CONSTS(124) =  CONSTS(56)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(57))
      CONSTS(125) =  CONSTS(58)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(59))
      CONSTS(126) =  CONSTS(61)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(64))
      CONSTS(127) =  CONSTS(62)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(65))
      CONSTS(128) =  CONSTS(63)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(66))
      CONSTS(129) =  CONSTS(68)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(69))
      CONSTS(130) =  CONSTS(72)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(73))
      CONSTS(131) =  CONSTS(74)*(CONSTS(117)+CONSTS(118))*(1.00000 - CONSTS(75))
      CONSTS(132) =  CONSTS(76)*(CONSTS(117)+CONSTS(118))*CONSTS(77)
      CONSTS(133) =  CONSTS(79)*(CONSTS(117)+CONSTS(118))*CONSTS(80)
      CONSTS(134) =  CONSTS(83)*(CONSTS(117)+CONSTS(118))*CONSTS(84)
      CONSTS(135) =  CONSTS(85)*(CONSTS(117)+CONSTS(118))*CONSTS(86)
      CONSTS(136) =  CONSTS(87)*(CONSTS(117)+CONSTS(118))*CONSTS(88)
      CONSTS(137) =  CONSTS(89)*(CONSTS(117)+CONSTS(118))*CONSTS(90)
      CONSTS(138) =  CONSTS(92)*(CONSTS(117)+CONSTS(118))*CONSTS(95)
      CONSTS(139) =  CONSTS(93)*(CONSTS(117)+CONSTS(118))*CONSTS(96)
      CONSTS(140) =  CONSTS(94)*(CONSTS(117)+CONSTS(118))*CONSTS(97)
      CONSTS(141) =  CONSTS(99)*(CONSTS(117)+CONSTS(118))*CONSTS(100)
      CONSTS(142) =  CONSTS(103)*(CONSTS(117)+CONSTS(118))*CONSTS(104)
      CONSTS(143) =  CONSTS(105)*(CONSTS(117)+CONSTS(118))*CONSTS(106)
      CONSTS(144) = CONSTS(118)/( 2.00000*3.14159*CONSTS(8)*CONSTS(119))
      CONSTS(145) = ( CONSTS(7)*(CONSTS(144)/2.00000))/( 3.14159*CONSTS(8) ** 2.00000*CONSTS(119))
      CONSTS(146) =  3.14159*CONSTS(8) ** 2.00000*CONSTS(144)*CONSTS(119)
      CONSTS(147) =  CONSTS(118)*1.00000
      CONSTS(148) =  CONSTS(117)*1.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(5) =  - CONSTS(12)*STATES(7) ** 4.00000*STATES(5)+ CONSTS(13)*STATES(3)
      RATES(3) = ( CONSTS(12)*STATES(7) ** 4.00000*STATES(5) - ( CONSTS(13)*STATES(3)+ CONSTS(14)*STATES(7) ** 3.00000*STATES(3)+ CONSTS(16)*STATES(3)))+ CONSTS(15)*STATES(4)+ CONSTS(17)*STATES(6)
      RATES(4) =  CONSTS(14)*STATES(7) ** 3.00000*STATES(3) -  CONSTS(15)*STATES(4)
      RATES(6) =  CONSTS(16)*STATES(3) -  CONSTS(17)*STATES(6)
      ALGBRC(6) = 1.00000/(1.00000+STATES(7)/0.0100000)
      RATES(21) = (ALGBRC(6) - STATES(21))/CONSTS(47)
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000))
      RATES(28) = (ALGBRC(13) - STATES(28))/CONSTS(108)
      ALGBRC(18) = 1.00000/(1.00000+STATES(7)/0.0100000)
      RATES(33) = (ALGBRC(18) - STATES(33))/CONSTS(78)
      ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(2)+87.5000)/10.3000))
      RATES(40) = (ALGBRC(25) - STATES(40))/CONSTS(111)
      ALGBRC(3) = ( 117.260*(STATES(1)+59.3000))/(1.00000 - EXP( - 0.550000*(STATES(1)+59.3000)))
      ALGBRC(28) =  3800.00*EXP( - 0.0720000*(STATES(1)+61.0000))
      RATES(18) =  ALGBRC(3)*(1.00000 - STATES(18)) -  ALGBRC(28)*STATES(18)
      ALGBRC(4) = 284.400/(1.00000+EXP( 0.0812000*(STATES(1)+115.900)))
      ALGBRC(29) = ( 18.7700*(STATES(1)+64.4000))/(1.00000 - EXP( - 0.220000*(STATES(1)+64.4000)))
      RATES(19) =  ALGBRC(4)*(1.00000 - STATES(19)) -  ALGBRC(29)*STATES(19)
      ALGBRC(30) =  0.00305000*EXP( - 0.00450000*STATES(1)+7.00000 ** 2.00000)+ 0.00105000*EXP( - 0.00200000*STATES(1) - 18.0000 ** 2.00000)+0.000250000
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+15.3000)/- 5.00000))
      RATES(20) = (ALGBRC(5) - STATES(20))/ALGBRC(30)
      ALGBRC(31) =  0.105000*EXP(- (STATES(1)+45.0000)/12.0000 ** 2.00000)+0.0400000/(1.00000+EXP((- STATES(1)+25.0000)/25.0000))+0.0150000/(1.00000+EXP((STATES(1)+75.0000)/25.0000))+0.00170000
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+26.7000)/5.40000))
      RATES(22) = (ALGBRC(7) - STATES(22))/ALGBRC(31)
      ALGBRC(32) =  0.0410000*EXP(- (STATES(1)+47.0000)/12.0000 ** 2.00000)+0.0800000/(1.00000+EXP((STATES(1)+55.0000)/- 5.00000))+0.0150000/(1.00000+EXP((STATES(1)+75.0000)/25.0000))+0.00170000
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+26.7000)/5.40000))
      RATES(23) = (ALGBRC(8) - STATES(23))/ALGBRC(32)
      ALGBRC(33) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200))
      RATES(24) = (ALGBRC(9) - STATES(24))/ALGBRC(33)
      ALGBRC(34) =  0.350000*EXP(-  0.0666000*(STATES(1)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
      RATES(25) = (ALGBRC(10) - STATES(25))/ALGBRC(34)
      ALGBRC(35) =  3.70000*EXP(-  0.0333000*(STATES(1)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
      RATES(26) = (ALGBRC(11) - STATES(26))/ALGBRC(35)
      ALGBRC(36) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+11.5000)/11.8200))
      RATES(27) = (ALGBRC(12) - STATES(27))/ALGBRC(36)
      ALGBRC(37) = 1.00000/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900))
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800))
      RATES(29) = (ALGBRC(14) - STATES(29))/ALGBRC(37)
      ALGBRC(15) = ( 117.260*(STATES(2)+59.3000))/(1.00000 - EXP( - 0.550000*(STATES(2)+59.3000)))
      ALGBRC(38) =  3800.00*EXP( - 0.0720000*(STATES(2)+61.0000))
      RATES(30) =  ALGBRC(15)*(1.00000 - STATES(30)) -  ALGBRC(38)*STATES(30)
      ALGBRC(16) = 284.400/(1.00000+EXP( 0.0812000*(STATES(2)+115.900)))
      ALGBRC(39) = ( 18.7700*(STATES(2)+64.4000))/(1.00000 - EXP( - 0.220000*(STATES(2)+64.4000)))
      RATES(31) =  ALGBRC(16)*(1.00000 - STATES(31)) -  ALGBRC(39)*STATES(31)
      ALGBRC(40) =  0.00305000*EXP( - 0.00450000*STATES(2)+7.00000 ** 2.00000)+ 0.00105000*EXP( - 0.00200000*STATES(2) - 18.0000 ** 2.00000)+0.000250000
      ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(2)+15.3000)/- 5.00000))
      RATES(32) = (ALGBRC(17) - STATES(32))/ALGBRC(40)
      ALGBRC(41) =  0.105000*EXP(- (STATES(2)+45.0000)/12.0000 ** 2.00000)+0.0400000/(1.00000+EXP((- STATES(2)+25.0000)/25.0000))+0.0150000/(1.00000+EXP((STATES(2)+75.0000)/25.0000))+0.00170000
      ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(2)+26.7000)/5.40000))
      RATES(34) = (ALGBRC(19) - STATES(34))/ALGBRC(41)
      ALGBRC(42) =  0.0410000*EXP(- (STATES(2)+47.0000)/12.0000 ** 2.00000)+0.0800000/(1.00000+EXP((STATES(2)+55.0000)/- 5.00000))+0.0150000/(1.00000+EXP((STATES(2)+75.0000)/25.0000))+0.00170000
      ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(2)+26.7000)/5.40000))
      RATES(35) = (ALGBRC(20) - STATES(35))/ALGBRC(42)
      ALGBRC(43) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(2)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(2)+38.0000)))
      ALGBRC(21) = 1.00000/(1.00000+EXP((STATES(2)+10.6000)/- 11.4200))
      RATES(36) = (ALGBRC(21) - STATES(36))/ALGBRC(43)
      ALGBRC(44) =  0.350000*EXP(-  0.0666000*(STATES(2)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(2)+45.3000)/6.88410))
      RATES(37) = (ALGBRC(22) - STATES(37))/ALGBRC(44)
      ALGBRC(45) =  3.70000*EXP(-  0.0333000*(STATES(2)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(2)+45.3000)/6.88410))
      RATES(38) = (ALGBRC(23) - STATES(38))/ALGBRC(45)
      ALGBRC(46) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(2)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(2)+38.0000)))
      ALGBRC(24) = 1.00000/(1.00000+EXP(- (STATES(2)+11.5000)/11.8200))
      RATES(39) = (ALGBRC(24) - STATES(39))/ALGBRC(46)
      ALGBRC(47) = 1.00000/( 0.118850*EXP((STATES(2)+80.0000)/28.3700)+ 0.562300*EXP((STATES(2)+80.0000)/- 14.1900))
      ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(2)+138.600)/10.4800))
      RATES(41) = (ALGBRC(26) - STATES(41))/ALGBRC(47)
      ALGBRC(52) =  CONSTS(29)*STATES(9)*(1.00000 - STATES(11)) -  CONSTS(30)*STATES(11)
      RATES(11) = ALGBRC(52)
      ALGBRC(27) =  1800.00*(STATES(3)+STATES(4))*(STATES(8) - STATES(7))
      ALGBRC(51) = (STATES(10) - STATES(8))/CONSTS(27)
      RATES(8) =  (1.00000/(1.00000+( CONSTS(37)*CONSTS(34))/CONSTS(34)+STATES(8) ** 2.00000))*(ALGBRC(51) - ALGBRC(27))
      ALGBRC(48) = STATES(9)/CONSTS(18) ** CONSTS(23)
      ALGBRC(49) = STATES(10)/CONSTS(19) ** CONSTS(24)
      ALGBRC(50) = ( CONSTS(22)*( CONSTS(20)*ALGBRC(48) -  CONSTS(21)*ALGBRC(49)))/(1.00000+ALGBRC(48)+ALGBRC(49))
      RATES(10) = ALGBRC(50) - ALGBRC(51)
      ALGBRC(53) =  CONSTS(31)*STATES(9)*(1.00000 - STATES(12)) -  CONSTS(32)*STATES(12)
      RATES(12) = ALGBRC(53)
      ALGBRC(2) =  1000.00*((STATES(1) - STATES(2))/CONSTS(11))
      ALGBRC(55) = TERNRY(CONSTS(40).EQ.1.00000.AND.VOI - CONSTS(41).LT.0.00100000,  52.0000*(CONSTS(117)+CONSTS(118)), 0.00000)
      ALGBRC(60) =  CONSTS(120)*STATES(18) ** 3.00000*STATES(19)*(STATES(1) -  arbitrary_log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(13)+ 0.120000*STATES(14)), 10)*(( CONSTS(1)*CONSTS(2))/CONSTS(3)))
      ALGBRC(62) =  CONSTS(121)*4.00000*(( STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*(( STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) -  0.341000*CONSTS(5))/(EXP(( 2.00000*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) - 1.00000))*STATES(20)*( (0.900000+STATES(21)/10.0000)*STATES(22)+ (0.100000+STATES(21)/10.0000)*STATES(23))
      ALGBRC(75) = ( CONSTS(131)*( EXP( 0.0187000*STATES(1))*STATES(13) ** 3.00000*CONSTS(5) -  EXP( -0.0187000*STATES(1))*CONSTS(4) ** 3.00000*STATES(9)))/(1.00000+ 0.000100000*( STATES(9)*CONSTS(4) ** 3.00000+ CONSTS(5)*STATES(13) ** 3.00000))
      ALGBRC(73) = (( (CONSTS(129)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*CONSTS(110)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*CONSTS(6))/(CONSTS(6)+CONSTS(70)))/(1.00000+CONSTS(71)/STATES(13) ** 1.50000)
      ALGBRC(63) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(6)/STATES(14))
      ALGBRC(66) =  (( CONSTS(124)*(STATES(1) - (ALGBRC(63)+1.73000)))/(1.00000+EXP( (( 1.61300*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*(STATES(1) - (ALGBRC(63)+1.73000)))))*(1.00000+EXP(- (CONSTS(6) - 0.998800)/0.124000))
      ALGBRC(64) =  CONSTS(122)*STATES(24)*( CONSTS(50)*STATES(25)+ CONSTS(51)*STATES(26))*(STATES(1) - ALGBRC(63))
      ALGBRC(65) =  CONSTS(123)*STATES(27)*STATES(28)*(STATES(1) - ALGBRC(63))
      ALGBRC(74) = ( CONSTS(130)*STATES(9))/(STATES(9)+0.000400000)
      ALGBRC(59) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(4)/STATES(13))
      ALGBRC(68) =  CONSTS(125)*STATES(29)*( 0.200000*(STATES(1) - ALGBRC(59))+ 0.800000*(STATES(1) - ALGBRC(63)))
      ALGBRC(69) =  CONSTS(127)*(STATES(1) - ALGBRC(59))
      ALGBRC(71) =  CONSTS(128)*(STATES(1) - ALGBRC(63))
      ALGBRC(70) =  CONSTS(126)*(STATES(1) - CONSTS(67))
      ALGBRC(72) = ALGBRC(69)+ALGBRC(70)+ALGBRC(71)
      RATES(1) =  1000.00*(((ALGBRC(55)+ALGBRC(2)) - (ALGBRC(60)+ALGBRC(62)+ALGBRC(64)+ALGBRC(65)+ALGBRC(68)+ALGBRC(66)+ALGBRC(72)+ALGBRC(73)+ALGBRC(75)+ALGBRC(74)))/CONSTS(148))
      ALGBRC(79) =  CONSTS(133)*4.00000*(( STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*(( STATES(7)*EXP(( 2.00000*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) -  0.341000*STATES(16))/(EXP(( 2.00000*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) - 1.00000))*STATES(32)*( (0.900000+STATES(33)/10.0000)*STATES(34)+ (0.100000+STATES(33)/10.0000)*STATES(35))
      RATES(7) =  (1.00000/(1.00000+ CONSTS(36)*(CONSTS(33)/CONSTS(33)+STATES(7) ** 2.00000)))*(- (ALGBRC(62)+ALGBRC(79))/(( 2.00000*CONSTS(113)*CONSTS(3)+ALGBRC(27)/CONSTS(113)) - (( (STATES(7) - STATES(9))*CONSTS(114))/CONSTS(39))/CONSTS(113)))
      ALGBRC(67) =  CONSTS(125)*STATES(29)*CONSTS(109)*(STATES(1) - ALGBRC(63))
      ALGBRC(90) = (EXP(STATES(15)/67.3000) - 1.00000)/7.00000
      ALGBRC(91) = (( (CONSTS(141)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*ALGBRC(90)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*STATES(17))/(STATES(17)+CONSTS(101)))/(1.00000+CONSTS(102)/STATES(13) ** 1.50000)
      ALGBRC(80) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(17)/STATES(14))
      ALGBRC(84) =  CONSTS(137)*STATES(41)*CONSTS(112)*(STATES(2) - ALGBRC(80))
      ALGBRC(88) =  CONSTS(140)*(STATES(2) - ALGBRC(80))
      ALGBRC(83) =  (( CONSTS(136)*(STATES(2) - (ALGBRC(80)+1.73000)))/(1.00000+EXP( (( 1.61300*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*(STATES(2) - (ALGBRC(80)+1.73000)))))*(1.00000+EXP(- (STATES(17) - 0.998800)/0.124000))
      ALGBRC(81) =  CONSTS(134)*STATES(36)*( CONSTS(81)*STATES(37)+ CONSTS(82)*STATES(38))*(STATES(2) - ALGBRC(80))
      ALGBRC(82) =  CONSTS(135)*STATES(39)*STATES(40)*(STATES(2) - ALGBRC(80))
      RATES(14) = - (- ALGBRC(55)+ALGBRC(67)+ALGBRC(84)+ALGBRC(65)+ALGBRC(82)+ALGBRC(66)+ALGBRC(83)+ALGBRC(71)+ALGBRC(88)+ ALGBRC(73)*- 2.00000+ ALGBRC(91)*- 2.00000+ALGBRC(64)+ALGBRC(81))/( CONSTS(114)*CONSTS(3))
      ALGBRC(58) = ( 1.00000*CONSTS(146))/( CONSTS(44)*(STATES(17) - CONSTS(6)))
      RATES(17) = ((ALGBRC(84)+ALGBRC(82)+ALGBRC(83)+ALGBRC(88)+ ALGBRC(91)*- 2.00000+ALGBRC(81))/(CONSTS(3) - ALGBRC(58)))/CONSTS(146)
      ALGBRC(77) =  CONSTS(132)*STATES(30) ** 3.00000*STATES(31)*(STATES(2) -  arbitrary_log((STATES(15)+ 0.120000*STATES(17))/(STATES(13)+ 0.120000*STATES(14)), 10)*(( CONSTS(1)*CONSTS(2))/CONSTS(3)))
      ALGBRC(76) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(15)/STATES(13))
      ALGBRC(86) =  CONSTS(139)*(STATES(2) - ALGBRC(76))
      ALGBRC(93) = ( CONSTS(143)*( EXP( 0.0187000*STATES(2))*STATES(13) ** 3.00000*STATES(16) -  EXP( -0.0187000*STATES(2))*STATES(15) ** 3.00000*STATES(9)))/(1.00000+ 0.000100000*( STATES(9)*STATES(15) ** 3.00000+ STATES(16)*STATES(13) ** 3.00000))
      RATES(13) = - (ALGBRC(60)+ALGBRC(77)+ALGBRC(69)+ALGBRC(86)+ ALGBRC(75)*3.00000+ ALGBRC(93)*3.00000+ ALGBRC(73)*3.00000+ ALGBRC(91)*3.00000)/( CONSTS(114)*CONSTS(3))
      ALGBRC(56) = ( 1.00000*CONSTS(146))/( CONSTS(42)*(STATES(15) - CONSTS(4)))
      RATES(15) = ((ALGBRC(77)+ALGBRC(86)+ ALGBRC(93)*3.00000+ ALGBRC(91)*3.00000)/(CONSTS(3) - ALGBRC(56)))/CONSTS(146)
      ALGBRC(92) = ( CONSTS(142)*STATES(9))/(STATES(9)+0.000400000)
      ALGBRC(87) =  CONSTS(138)*(STATES(2) - CONSTS(98))
      RATES(9) =  (1.00000/(1.00000+ CONSTS(36)*(CONSTS(33)/CONSTS(33)+STATES(9) ** 2.00000)))*((( -2.00000*(ALGBRC(75)+ALGBRC(93))+(ALGBRC(70)+ALGBRC(87)+ALGBRC(74)+ALGBRC(92)))/( -2.00000*CONSTS(114)*CONSTS(3))+((( (STATES(7) - STATES(9))*CONSTS(114))/CONSTS(39)+ALGBRC(27)) - ALGBRC(50))/CONSTS(114)) - ( ALGBRC(52)*CONSTS(25)+ ALGBRC(53)*CONSTS(26)))
      ALGBRC(57) = ( 1.00000*CONSTS(146))/( CONSTS(43)*(STATES(16) - CONSTS(5)))
      RATES(16) = ((ALGBRC(79)+ALGBRC(92)+ALGBRC(87)+ ALGBRC(93)*- 2.00000)/( 2.00000*CONSTS(3) - ALGBRC(57)))/CONSTS(146)
      ALGBRC(85) =  CONSTS(137)*STATES(41)*( 0.200000*(STATES(2) - ALGBRC(76))+ 0.800000*(STATES(2) - ALGBRC(80)))
      ALGBRC(89) = ALGBRC(86)+ALGBRC(87)+ALGBRC(88)
      RATES(2) =  1000.00*(- (ALGBRC(2)+ALGBRC(77)+ALGBRC(79)+ALGBRC(81)+ALGBRC(82)+ALGBRC(85)+ALGBRC(83)+ALGBRC(89)+ALGBRC(91)+ALGBRC(93)+ALGBRC(92))/CONSTS(147))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(6) = 1.00000/(1.00000+STATES(7)/0.0100000)
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000))
      ALGBRC(18) = 1.00000/(1.00000+STATES(7)/0.0100000)
      ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(2)+87.5000)/10.3000))
      ALGBRC(3) = ( 117.260*(STATES(1)+59.3000))/(1.00000 - EXP( - 0.550000*(STATES(1)+59.3000)))
      ALGBRC(28) =  3800.00*EXP( - 0.0720000*(STATES(1)+61.0000))
      ALGBRC(4) = 284.400/(1.00000+EXP( 0.0812000*(STATES(1)+115.900)))
      ALGBRC(29) = ( 18.7700*(STATES(1)+64.4000))/(1.00000 - EXP( - 0.220000*(STATES(1)+64.4000)))
      ALGBRC(30) =  0.00305000*EXP( - 0.00450000*STATES(1)+7.00000 ** 2.00000)+ 0.00105000*EXP( - 0.00200000*STATES(1) - 18.0000 ** 2.00000)+0.000250000
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+15.3000)/- 5.00000))
      ALGBRC(31) =  0.105000*EXP(- (STATES(1)+45.0000)/12.0000 ** 2.00000)+0.0400000/(1.00000+EXP((- STATES(1)+25.0000)/25.0000))+0.0150000/(1.00000+EXP((STATES(1)+75.0000)/25.0000))+0.00170000
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+26.7000)/5.40000))
      ALGBRC(32) =  0.0410000*EXP(- (STATES(1)+47.0000)/12.0000 ** 2.00000)+0.0800000/(1.00000+EXP((STATES(1)+55.0000)/- 5.00000))+0.0150000/(1.00000+EXP((STATES(1)+75.0000)/25.0000))+0.00170000
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+26.7000)/5.40000))
      ALGBRC(33) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200))
      ALGBRC(34) =  0.350000*EXP(-  0.0666000*(STATES(1)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
      ALGBRC(35) =  3.70000*EXP(-  0.0333000*(STATES(1)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
      ALGBRC(36) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+11.5000)/11.8200))
      ALGBRC(37) = 1.00000/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900))
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800))
      ALGBRC(15) = ( 117.260*(STATES(2)+59.3000))/(1.00000 - EXP( - 0.550000*(STATES(2)+59.3000)))
      ALGBRC(38) =  3800.00*EXP( - 0.0720000*(STATES(2)+61.0000))
      ALGBRC(16) = 284.400/(1.00000+EXP( 0.0812000*(STATES(2)+115.900)))
      ALGBRC(39) = ( 18.7700*(STATES(2)+64.4000))/(1.00000 - EXP( - 0.220000*(STATES(2)+64.4000)))
      ALGBRC(40) =  0.00305000*EXP( - 0.00450000*STATES(2)+7.00000 ** 2.00000)+ 0.00105000*EXP( - 0.00200000*STATES(2) - 18.0000 ** 2.00000)+0.000250000
      ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(2)+15.3000)/- 5.00000))
      ALGBRC(41) =  0.105000*EXP(- (STATES(2)+45.0000)/12.0000 ** 2.00000)+0.0400000/(1.00000+EXP((- STATES(2)+25.0000)/25.0000))+0.0150000/(1.00000+EXP((STATES(2)+75.0000)/25.0000))+0.00170000
      ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(2)+26.7000)/5.40000))
      ALGBRC(42) =  0.0410000*EXP(- (STATES(2)+47.0000)/12.0000 ** 2.00000)+0.0800000/(1.00000+EXP((STATES(2)+55.0000)/- 5.00000))+0.0150000/(1.00000+EXP((STATES(2)+75.0000)/25.0000))+0.00170000
      ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(2)+26.7000)/5.40000))
      ALGBRC(43) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(2)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(2)+38.0000)))
      ALGBRC(21) = 1.00000/(1.00000+EXP((STATES(2)+10.6000)/- 11.4200))
      ALGBRC(44) =  0.350000*EXP(-  0.0666000*(STATES(2)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(2)+45.3000)/6.88410))
      ALGBRC(45) =  3.70000*EXP(-  0.0333000*(STATES(2)+70.0000) ** 2.00000)+0.0350000
      ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(2)+45.3000)/6.88410))
      ALGBRC(46) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(2)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(2)+38.0000)))
      ALGBRC(24) = 1.00000/(1.00000+EXP(- (STATES(2)+11.5000)/11.8200))
      ALGBRC(47) = 1.00000/( 0.118850*EXP((STATES(2)+80.0000)/28.3700)+ 0.562300*EXP((STATES(2)+80.0000)/- 14.1900))
      ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(2)+138.600)/10.4800))
      ALGBRC(52) =  CONSTS(29)*STATES(9)*(1.00000 - STATES(11)) -  CONSTS(30)*STATES(11)
      ALGBRC(27) =  1800.00*(STATES(3)+STATES(4))*(STATES(8) - STATES(7))
      ALGBRC(51) = (STATES(10) - STATES(8))/CONSTS(27)
      ALGBRC(48) = STATES(9)/CONSTS(18) ** CONSTS(23)
      ALGBRC(49) = STATES(10)/CONSTS(19) ** CONSTS(24)
      ALGBRC(50) = ( CONSTS(22)*( CONSTS(20)*ALGBRC(48) -  CONSTS(21)*ALGBRC(49)))/(1.00000+ALGBRC(48)+ALGBRC(49))
      ALGBRC(53) =  CONSTS(31)*STATES(9)*(1.00000 - STATES(12)) -  CONSTS(32)*STATES(12)
      ALGBRC(2) =  1000.00*((STATES(1) - STATES(2))/CONSTS(11))
      ALGBRC(55) = TERNRY(CONSTS(40).EQ.1.00000.AND.VOI - CONSTS(41).LT.0.00100000,  52.0000*(CONSTS(117)+CONSTS(118)), 0.00000)
      ALGBRC(60) =  CONSTS(120)*STATES(18) ** 3.00000*STATES(19)*(STATES(1) -  arbitrary_log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(13)+ 0.120000*STATES(14)), 10)*(( CONSTS(1)*CONSTS(2))/CONSTS(3)))
      ALGBRC(62) =  CONSTS(121)*4.00000*(( STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*(( STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) -  0.341000*CONSTS(5))/(EXP(( 2.00000*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) - 1.00000))*STATES(20)*( (0.900000+STATES(21)/10.0000)*STATES(22)+ (0.100000+STATES(21)/10.0000)*STATES(23))
      ALGBRC(75) = ( CONSTS(131)*( EXP( 0.0187000*STATES(1))*STATES(13) ** 3.00000*CONSTS(5) -  EXP( -0.0187000*STATES(1))*CONSTS(4) ** 3.00000*STATES(9)))/(1.00000+ 0.000100000*( STATES(9)*CONSTS(4) ** 3.00000+ CONSTS(5)*STATES(13) ** 3.00000))
      ALGBRC(73) = (( (CONSTS(129)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*CONSTS(110)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*CONSTS(6))/(CONSTS(6)+CONSTS(70)))/(1.00000+CONSTS(71)/STATES(13) ** 1.50000)
      ALGBRC(63) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(6)/STATES(14))
      ALGBRC(66) =  (( CONSTS(124)*(STATES(1) - (ALGBRC(63)+1.73000)))/(1.00000+EXP( (( 1.61300*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*(STATES(1) - (ALGBRC(63)+1.73000)))))*(1.00000+EXP(- (CONSTS(6) - 0.998800)/0.124000))
      ALGBRC(64) =  CONSTS(122)*STATES(24)*( CONSTS(50)*STATES(25)+ CONSTS(51)*STATES(26))*(STATES(1) - ALGBRC(63))
      ALGBRC(65) =  CONSTS(123)*STATES(27)*STATES(28)*(STATES(1) - ALGBRC(63))
      ALGBRC(74) = ( CONSTS(130)*STATES(9))/(STATES(9)+0.000400000)
      ALGBRC(59) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(4)/STATES(13))
      ALGBRC(68) =  CONSTS(125)*STATES(29)*( 0.200000*(STATES(1) - ALGBRC(59))+ 0.800000*(STATES(1) - ALGBRC(63)))
      ALGBRC(69) =  CONSTS(127)*(STATES(1) - ALGBRC(59))
      ALGBRC(71) =  CONSTS(128)*(STATES(1) - ALGBRC(63))
      ALGBRC(70) =  CONSTS(126)*(STATES(1) - CONSTS(67))
      ALGBRC(72) = ALGBRC(69)+ALGBRC(70)+ALGBRC(71)
      ALGBRC(79) =  CONSTS(133)*4.00000*(( STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*(( STATES(7)*EXP(( 2.00000*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) -  0.341000*STATES(16))/(EXP(( 2.00000*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2))) - 1.00000))*STATES(32)*( (0.900000+STATES(33)/10.0000)*STATES(34)+ (0.100000+STATES(33)/10.0000)*STATES(35))
      ALGBRC(67) =  CONSTS(125)*STATES(29)*CONSTS(109)*(STATES(1) - ALGBRC(63))
      ALGBRC(90) = (EXP(STATES(15)/67.3000) - 1.00000)/7.00000
      ALGBRC(91) = (( (CONSTS(141)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*ALGBRC(90)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*STATES(17))/(STATES(17)+CONSTS(101)))/(1.00000+CONSTS(102)/STATES(13) ** 1.50000)
      ALGBRC(80) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(17)/STATES(14))
      ALGBRC(84) =  CONSTS(137)*STATES(41)*CONSTS(112)*(STATES(2) - ALGBRC(80))
      ALGBRC(88) =  CONSTS(140)*(STATES(2) - ALGBRC(80))
      ALGBRC(83) =  (( CONSTS(136)*(STATES(2) - (ALGBRC(80)+1.73000)))/(1.00000+EXP( (( 1.61300*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*(STATES(2) - (ALGBRC(80)+1.73000)))))*(1.00000+EXP(- (STATES(17) - 0.998800)/0.124000))
      ALGBRC(81) =  CONSTS(134)*STATES(36)*( CONSTS(81)*STATES(37)+ CONSTS(82)*STATES(38))*(STATES(2) - ALGBRC(80))
      ALGBRC(82) =  CONSTS(135)*STATES(39)*STATES(40)*(STATES(2) - ALGBRC(80))
      ALGBRC(58) = ( 1.00000*CONSTS(146))/( CONSTS(44)*(STATES(17) - CONSTS(6)))
      ALGBRC(77) =  CONSTS(132)*STATES(30) ** 3.00000*STATES(31)*(STATES(2) -  arbitrary_log((STATES(15)+ 0.120000*STATES(17))/(STATES(13)+ 0.120000*STATES(14)), 10)*(( CONSTS(1)*CONSTS(2))/CONSTS(3)))
      ALGBRC(76) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(15)/STATES(13))
      ALGBRC(86) =  CONSTS(139)*(STATES(2) - ALGBRC(76))
      ALGBRC(93) = ( CONSTS(143)*( EXP( 0.0187000*STATES(2))*STATES(13) ** 3.00000*STATES(16) -  EXP( -0.0187000*STATES(2))*STATES(15) ** 3.00000*STATES(9)))/(1.00000+ 0.000100000*( STATES(9)*STATES(15) ** 3.00000+ STATES(16)*STATES(13) ** 3.00000))
      ALGBRC(56) = ( 1.00000*CONSTS(146))/( CONSTS(42)*(STATES(15) - CONSTS(4)))
      ALGBRC(92) = ( CONSTS(142)*STATES(9))/(STATES(9)+0.000400000)
      ALGBRC(87) =  CONSTS(138)*(STATES(2) - CONSTS(98))
      ALGBRC(57) = ( 1.00000*CONSTS(146))/( CONSTS(43)*(STATES(16) - CONSTS(5)))
      ALGBRC(85) =  CONSTS(137)*STATES(41)*( 0.200000*(STATES(2) - ALGBRC(76))+ 0.800000*(STATES(2) - ALGBRC(80)))
      ALGBRC(89) = ALGBRC(86)+ALGBRC(87)+ALGBRC(88)
      ALGBRC(1) = (STATES(7) - STATES(9))/CONSTS(28)
      ALGBRC(54) = ALGBRC(52)+ALGBRC(53)
      ALGBRC(61) =  CONSTS(125)*STATES(29)*CONSTS(60)*(STATES(1) - ALGBRC(59))
      ALGBRC(78) =  CONSTS(137)*STATES(41)*CONSTS(91)*(STATES(2) - ALGBRC(76))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END