Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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