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 121 entries in the algebraic variable array.
C There are a total of 31 entries in each of the rate and state variable arrays.
C There are a total of 130 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is protocol in component environment (dimensionless).
C CONSTS(2) is L_totmax in component b1_AR_Gs_parameters (uM).
C CONSTS(3) is sum_b1_AR in component b1_AR_Gs_parameters (uM).
C CONSTS(4) is Gs_tot in component b1_AR_Gs_parameters (uM).
C CONSTS(5) is Kl in component b1_AR_Gs_parameters (uM).
C CONSTS(6) is Kr in component b1_AR_Gs_parameters (uM).
C CONSTS(7) is Kc in component b1_AR_Gs_parameters (uM).
C CONSTS(8) is k_bar_kp in component b1_AR_Gs_parameters (per_sec).
C CONSTS(9) is k_bar_km in component b1_AR_Gs_parameters (per_sec).
C CONSTS(10) is k_p_kap in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(11) is k_p_kam in component b1_AR_Gs_parameters (per_sec).
C CONSTS(12) is k_g_act in component b1_AR_Gs_parameters (per_sec).
C CONSTS(13) is k_hyd in component b1_AR_Gs_parameters (per_sec).
C CONSTS(14) is k_reassoc in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(15) is AC_tot in component cAMP_parameters (uM).
C CONSTS(16) is ATP in component cAMP_parameters (uM).
C CONSTS(17) is PDE_tot in component cAMP_parameters (uM).
C CONSTS(18) is IBMX_tot in component cAMP_parameters (uM).
C CONSTS(19) is Fsk_tot in component cAMP_parameters (uM).
C CONSTS(20) is k_ac_basal in component cAMP_parameters (per_sec).
C CONSTS(21) is k_ac_gsa in component cAMP_parameters (per_sec).
C CONSTS(22) is k_ac_fsk in component cAMP_parameters (per_sec).
C CONSTS(23) is k_pde in component cAMP_parameters (per_sec).
C CONSTS(24) is Km_basal in component cAMP_parameters (uM).
C CONSTS(25) is Km_gsa in component cAMP_parameters (uM).
C CONSTS(26) is Km_fsk in component cAMP_parameters (uM).
C CONSTS(27) is Km_pde in component cAMP_parameters (uM).
C CONSTS(28) is K_gsa in component cAMP_parameters (uM).
C CONSTS(29) is K_fsk in component cAMP_parameters (uM).
C CONSTS(30) is Ki_ibmx in component cAMP_parameters (uM).
C CONSTS(31) is PKAI_tot in component PKA_parameters (uM).
C CONSTS(32) is PKAII_tot in component PKA_parameters (uM).
C CONSTS(33) is PKI_tot in component PKA_parameters (uM).
C CONSTS(34) is K_a in component PKA_parameters (uM).
C CONSTS(35) is K_b in component PKA_parameters (uM).
C CONSTS(36) is K_d in component PKA_parameters (uM).
C CONSTS(37) is Ki_pki in component PKA_parameters (uM).
C CONSTS(38) is PLB_tot in component PLB_parameters (uM).
C CONSTS(39) is PP1_tot in component PLB_parameters (uM).
C CONSTS(40) is Inhib1_tot in component PLB_parameters (uM).
C CONSTS(41) is k_pka_plb in component PLB_parameters (per_sec).
C CONSTS(42) is Km_pka_plb in component PLB_parameters (uM).
C CONSTS(43) is k_pp1_plb in component PLB_parameters (per_sec).
C CONSTS(44) is Km_pp1_plb in component PLB_parameters (uM).
C CONSTS(45) is k_pka_i1 in component PLB_parameters (per_sec).
C CONSTS(46) is Km_pka_i1 in component PLB_parameters (uM).
C CONSTS(47) is Vmax_pp2a_i1 in component PLB_parameters (uM_per_sec).
C CONSTS(48) is Km_pp2a_i1 in component PLB_parameters (uM).
C CONSTS(49) is Ki_inhib1 in component PLB_parameters (uM).
C CONSTS(50) is epsilon in component LCC_parameters (dimensionless).
C CONSTS(51) is LCC_tot in component LCC_parameters (uM).
C CONSTS(52) is PP1_lcc_tot in component LCC_parameters (uM).
C CONSTS(53) is PP2A_lcc_tot in component LCC_parameters (uM).
C CONSTS(54) is k_pka_lcc in component LCC_parameters (per_sec).
C CONSTS(55) is Km_pka_lcc in component LCC_parameters (uM).
C CONSTS(56) is k_pp1_lcc in component LCC_parameters (per_sec).
C CONSTS(57) is Km_pp1_lcc in component LCC_parameters (uM).
C CONSTS(58) is k_pp2a_lcc in component LCC_parameters (per_sec).
C CONSTS(59) is Km_pp2a_lcc in component LCC_parameters (uM).
C CONSTS(60) is V_myo in component EC_Coupling_Parameters (uL).
C CONSTS(61) is V_nsr in component EC_Coupling_Parameters (uL).
C CONSTS(62) is V_jsr in component EC_Coupling_Parameters (uL).
C CONSTS(63) is A_Cap in component EC_Coupling_Parameters (cm2).
C CONSTS(64) is Temp in component EC_Coupling_Parameters (kelvin).
C CONSTS(65) is Na_ext in component EC_Coupling_Parameters (mM).
C CONSTS(66) is K_ext in component EC_Coupling_Parameters (mM).
C CONSTS(67) is Ca_ext in component EC_Coupling_Parameters (mM).
C CONSTS(68) is G_Na in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(69) is G_to in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(70) is G_ss in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(71) is G_Ki_bar in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(72) is G_Kp in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(73) is f in component EC_Coupling_Parameters (per_sec).
C CONSTS(74) is g in component EC_Coupling_Parameters (per_sec).
C CONSTS(75) is gamma_o in component EC_Coupling_Parameters (per_mM_per_sec).
C CONSTS(76) is omega in component EC_Coupling_Parameters (per_sec).
C CONSTS(77) is p_Ca in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(78) is p_K in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(79) is N_lcc in component EC_Coupling_Parameters (dimensionless).
C CONSTS(80) is I_Ca05 in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(81) is k_NaCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(82) is Km_Na in component EC_Coupling_Parameters (mM).
C CONSTS(83) is Km_Ca in component EC_Coupling_Parameters (mM).
C CONSTS(84) is k_sat in component EC_Coupling_Parameters (dimensionless).
C CONSTS(85) is eta in component EC_Coupling_Parameters (dimensionless).
C CONSTS(86) is I_bar_NaK in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(87) is Km_Nai in component EC_Coupling_Parameters (mM).
C CONSTS(88) is Km_Ko in component EC_Coupling_Parameters (mM).
C CONSTS(89) is I_bar_PCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(90) is Km_PCa in component EC_Coupling_Parameters (mM).
C CONSTS(91) is G_CaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(92) is G_NaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(93) is Pns in component EC_Coupling_Parameters (dimensionless).
C CONSTS(94) is Km_NS in component EC_Coupling_Parameters (mM).
C CONSTS(95) is I_up_bar in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(96) is Km_up0 in component EC_Coupling_Parameters (mM).
C CONSTS(97) is NSR_bar in component EC_Coupling_Parameters (mM).
C CONSTS(98) is tau_on in component EC_Coupling_Parameters (second).
C CONSTS(99) is tau_off in component EC_Coupling_Parameters (second).
C CONSTS(100) is G_max_rel in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(101) is d_Cai_th in component EC_Coupling_Parameters (mM).
C CONSTS(102) is Km_rel in component EC_Coupling_Parameters (mM).
C CONSTS(103) is CSQN_th in component EC_Coupling_Parameters (mM).
C CONSTS(104) is CSQN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(105) is Km_CSQN in component EC_Coupling_Parameters (mM).
C CONSTS(106) is tau_tr in component EC_Coupling_Parameters (second).
C CONSTS(107) is TRPN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(108) is CMDN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(109) is INDO_bar in component EC_Coupling_Parameters (mM).
C CONSTS(110) is Km_TRPN in component EC_Coupling_Parameters (mM).
C CONSTS(111) is Km_CMDN in component EC_Coupling_Parameters (mM).
C CONSTS(112) is Km_INDO in component EC_Coupling_Parameters (mM).
C ALGBRC(62) is LR in component b1_AR_module (uM).
C ALGBRC(63) is LRG in component b1_AR_module (uM).
C ALGBRC(64) is RG in component b1_AR_module (uM).
C ALGBRC(80) is BARK_DESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(12) is BARK_RESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(88) is PKA_DESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(23) is PKA_RESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(81) is G_ACT in component b1_AR_module (uM_per_sec).
C ALGBRC(25) is HYD in component b1_AR_module (uM_per_sec).
C ALGBRC(27) is REASSOC in component b1_AR_module (per_sec).
C ALGBRC(65) is L in component b1_AR_module (uM).
C ALGBRC(66) is R in component b1_AR_module (uM).
C ALGBRC(67) is Gs in component b1_AR_module (uM).
C STATES(1) is b1_AR_tot in component b1_AR_module (uM).
C STATES(2) is b1_AR_d in component b1_AR_module (uM).
C STATES(3) is b1_AR_p in component b1_AR_module (uM).
C STATES(4) is Gs_agtp_tot in component b1_AR_module (uM).
C STATES(5) is Gs_agdp in component b1_AR_module (uM).
C STATES(6) is Gs_bg in component b1_AR_module (uM).
C ALGBRC(68) is PKAC_I in component PKA_module (uM).
C ALGBRC(69) is cAMP in component PKA_module (uM).
C ALGBRC(42) is Gsa_GTP in component cAMP_module (uM).
C ALGBRC(50) is Fsk in component cAMP_module (uM).
C ALGBRC(43) is AC in component cAMP_module (uM).
C CONSTS(128) is PDE in component cAMP_module (uM).
C CONSTS(129) is IBMX in component cAMP_module (uM).
C STATES(7) is cAMP_tot in component cAMP_module (uM).
C ALGBRC(44) is Gsa_GTP_AC in component cAMP_module (uM).
C ALGBRC(51) is Fsk_AC in component cAMP_module (uM).
C ALGBRC(48) is AC_ACT_GSA in component cAMP_module (uM).
C ALGBRC(46) is AC_ACT_BASAL in component cAMP_module (uM).
C ALGBRC(53) is AC_ACT_FSK in component cAMP_module (uM).
C ALGBRC(82) is PDE_ACT in component cAMP_module (uM).
C CONSTS(130) is PDE_IBMX in component cAMP_module (uM).
C ALGBRC(70) is PKI in component PKA_module (uM).
C ALGBRC(71) is A2RC_I in component PKA_module (uM).
C ALGBRC(72) is A2R_I in component PKA_module (uM).
C ALGBRC(73) is A2RC_II in component PKA_module (uM).
C ALGBRC(74) is A2R_II in component PKA_module (uM).
C ALGBRC(75) is ARC_I in component PKA_module (uM).
C ALGBRC(76) is ARC_II in component PKA_module (uM).
C ALGBRC(77) is PKA_temp in component PKA_module (uM).
C ALGBRC(78) is PKAC_II in component PKA_module (uM).
C ALGBRC(29) is PLB in component PLB_module (uM).
C ALGBRC(83) is PLB_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(60) is PLB_DEPHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(30) is Inhib1 in component PLB_module (uM).
C ALGBRC(55) is Inhib1p_PP1 in component PLB_module (uM).
C ALGBRC(84) is Inhib1_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(32) is Inhib1_DEPHOSPH in component PLB_module (uM_per_sec).
C STATES(8) is PLB_p in component PLB_module (uM).
C STATES(9) is Inhib1p_tot in component PLB_module (uM).
C ALGBRC(56) is Inhib1p in component PLB_module (uM).
C ALGBRC(57) is PP1 in component PLB_module (uM).
C ALGBRC(1) is frac_PLB_p in component PLB_module (dimensionless).
C ALGBRC(31) is frac_PLB in component PLB_module (dimensionless).
C CONSTS(121) is frac_PLB_o in component PLB_module (dimensionless).
C ALGBRC(34) is LCCa in component LCC_module (uM).
C ALGBRC(85) is LCCa_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(36) is LCCa_DEPHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(38) is LCCb in component LCC_module (uM).
C ALGBRC(86) is LCCb_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(40) is LCCb_DEPHOSPH in component LCC_module (uM_per_sec).
C STATES(10) is LCCa_p in component LCC_module (uM).
C STATES(11) is LCCb_p in component LCC_module (uM).
C ALGBRC(2) is frac_LCCa_p in component LCC_module (dimensionless).
C CONSTS(122) is frac_LCCa_po in component LCC_module (dimensionless).
C ALGBRC(33) is frac_LCCb_p in component LCC_module (dimensionless).
C CONSTS(123) is frac_LCCb_po in component LCC_module (dimensionless).
C ALGBRC(35) is E_Na in component Nernst_Potentials (mV).
C ALGBRC(37) is E_K in component Nernst_Potentials (mV).
C ALGBRC(39) is E_Ca in component Nernst_Potentials (mV).
C CONSTS(124) is E_Cl in component Nernst_Potentials (mV).
C CONSTS(113) is R in component Nernst_Potentials (joules_per_mole_kelvin).
C CONSTS(114) is Frdy in component Nernst_Potentials (coulombs_per_mole).
C CONSTS(125) is FoRT in component Nernst_Potentials (per_mV).
C CONSTS(115) is z_Na in component Nernst_Potentials (dimensionless).
C CONSTS(116) is z_K in component Nernst_Potentials (dimensionless).
C CONSTS(117) is z_Ca in component Nernst_Potentials (dimensionless).
C STATES(12) is Na_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C STATES(13) is K_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C STATES(14) is Ca_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C ALGBRC(3) is am in component Fast_Na_Current (per_sec).
C ALGBRC(13) is bm in component Fast_Na_Current (per_sec).
C ALGBRC(4) is ah in component Fast_Na_Current (per_sec).
C ALGBRC(5) is aj in component Fast_Na_Current (per_sec).
C ALGBRC(14) is bh in component Fast_Na_Current (per_sec).
C ALGBRC(15) is bj in component Fast_Na_Current (per_sec).
C STATES(15) is m in component Fast_Na_Current (dimensionless).
C STATES(16) is h in component Fast_Na_Current (dimensionless).
C STATES(17) is j in component Fast_Na_Current (dimensionless).
C ALGBRC(41) is I_Na in component Fast_Na_Current (uA_per_uF).
C STATES(18) is V_m in component Ion_Concentrations_and_Membrane_Potential (mV).
C ALGBRC(6) is a_lcc in component L_Type_Calcium_Current (per_sec).
C ALGBRC(16) is b_lcc in component L_Type_Calcium_Current (per_sec).
C ALGBRC(18) is f_lcc in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(7) is y_lcc_inf in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(17) is tau_y_lcc in component L_Type_Calcium_Current (second).
C ALGBRC(24) is gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(26) is v_gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(28) is v_omega in component L_Type_Calcium_Current (dimensionless).
C STATES(19) is v in component L_Type_Calcium_Current (dimensionless).
C STATES(20) is w in component L_Type_Calcium_Current (dimensionless).
C STATES(21) is x in component L_Type_Calcium_Current (dimensionless).
C STATES(22) is y in component L_Type_Calcium_Current (dimensionless).
C STATES(23) is z in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(45) is i_bar_Ca in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(47) is i_bar_K in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(49) is f_avail in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(52) is I_Ca in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(54) is I_CaK in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(58) is I_Ca_tot in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(8) is r_toss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(9) is s_toss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(19) is tau_r_to in component Transient_Outward_K_Current (second).
C ALGBRC(20) is tau_s_to in component Transient_Outward_K_Current (second).
C ALGBRC(21) is tau_ss_to in component Transient_Outward_K_Current (second).
C STATES(24) is r_to in component Transient_Outward_K_Current (dimensionless).
C STATES(25) is s_to in component Transient_Outward_K_Current (dimensionless).
C STATES(26) is ss_to in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(59) is I_to in component Transient_Outward_K_Current (uA_per_uF).
C ALGBRC(10) is r_ss_inf in component Steady_State_K_Current (dimensionless).
C ALGBRC(22) is tau_r_ss in component Steady_State_K_Current (second).
C ALGBRC(11) is s_ss_inf in component Steady_State_K_Current (dimensionless).
C CONSTS(126) is tau_s_ss in component Steady_State_K_Current (second).
C STATES(27) is r_ss in component Steady_State_K_Current (dimensionless).
C STATES(28) is s_ss in component Steady_State_K_Current (dimensionless).
C ALGBRC(61) is I_ss in component Steady_State_K_Current (uA_per_uF).
C ALGBRC(79) is a_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(87) is b_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(89) is Ki_ss in component Time_Independent_K_Current (dimensionless).
C ALGBRC(90) is I_Ki in component Time_Independent_K_Current (uA_per_uF).
C ALGBRC(91) is Kp in component Plateau_K_Current (dimensionless).
C ALGBRC(92) is I_Kp in component Plateau_K_Current (uA_per_uF).
C ALGBRC(93) is s4 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(94) is s5 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(95) is I_NCX in component Na_Ca_Exchanger_Current (uA_per_uF).
C CONSTS(127) is sigma in component Na_K_Pump_Current (dimensionless).
C ALGBRC(96) is f_NaK in component Na_K_Pump_Current (dimensionless).
C ALGBRC(97) is I_NaK in component Na_K_Pump_Current (uA_per_uF).
C ALGBRC(98) is I_PCa in component Sarcolemmal_Ca_Pump_Current (uA_per_uF).
C ALGBRC(99) is I_CaB in component Ca_Background_Current (uA_per_uF).
C ALGBRC(100) is I_NaB in component Na_Background_Current (uA_per_uF).
C ALGBRC(101) is I_Na_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(102) is I_K_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(103) is I_Ca_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(104) is t_rel in component Calcium_Induced_Calcium_Release (second).
C ALGBRC(106) is ryr_on in component Calcium_Induced_Calcium_Release (dimensionless).
C ALGBRC(108) is ryr_off in component Calcium_Induced_Calcium_Release (dimensionless).
C ALGBRC(110) is g_rel in component Calcium_Induced_Calcium_Release (per_sec).
C ALGBRC(111) is I_rel in component Calcium_Induced_Calcium_Release (mM_per_sec).
C STATES(29) is Ca_jsr in component Other_SR_Fluxes_and_Concentrations (mM).
C STATES(30) is trel in component Ion_Concentrations_and_Membrane_Potential (second).
C ALGBRC(112) is Km_up in component Other_SR_Fluxes_and_Concentrations (mM).
C ALGBRC(113) is I_up in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(114) is I_leak in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(115) is I_tr in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(117) is B_jsr in component Other_SR_Fluxes_and_Concentrations (dimensionless).
C STATES(31) is Ca_nsr in component Other_SR_Fluxes_and_Concentrations (mM).
C ALGBRC(119) is SR_content in component Other_SR_Fluxes_and_Concentrations (uM).
C ALGBRC(116) is b_trpn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(118) is b_cmdn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(120) is b_indo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(121) is B_myo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(109) is I_app in component Ion_Concentrations_and_Membrane_Potential (uA_per_uF).
C CONSTS(118) is V_hold in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(119) is V_test in component Ion_Concentrations_and_Membrane_Potential (mV).
C ALGBRC(105) is V_clamp in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(120) is R_clamp in component Ion_Concentrations_and_Membrane_Potential (dimensionless).
C ALGBRC(107) is I_pace in component Ion_Concentrations_and_Membrane_Potential (uA_per_uF).
C RATES(1) is d/dt b1_AR_tot in component b1_AR_module (uM).
C RATES(2) is d/dt b1_AR_d in component b1_AR_module (uM).
C RATES(3) is d/dt b1_AR_p in component b1_AR_module (uM).
C RATES(4) is d/dt Gs_agtp_tot in component b1_AR_module (uM).
C RATES(5) is d/dt Gs_agdp in component b1_AR_module (uM).
C RATES(6) is d/dt Gs_bg in component b1_AR_module (uM).
C RATES(7) is d/dt cAMP_tot in component cAMP_module (uM).
C RATES(8) is d/dt PLB_p in component PLB_module (uM).
C RATES(9) is d/dt Inhib1p_tot in component PLB_module (uM).
C RATES(10) is d/dt LCCa_p in component LCC_module (uM).
C RATES(11) is d/dt LCCb_p in component LCC_module (uM).
C RATES(15) is d/dt m in component Fast_Na_Current (dimensionless).
C RATES(16) is d/dt h in component Fast_Na_Current (dimensionless).
C RATES(17) is d/dt j in component Fast_Na_Current (dimensionless).
C RATES(19) is d/dt v in component L_Type_Calcium_Current (dimensionless).
C RATES(20) is d/dt w in component L_Type_Calcium_Current (dimensionless).
C RATES(21) is d/dt x in component L_Type_Calcium_Current (dimensionless).
C RATES(22) is d/dt y in component L_Type_Calcium_Current (dimensionless).
C RATES(23) is d/dt z in component L_Type_Calcium_Current (dimensionless).
C RATES(24) is d/dt r_to in component Transient_Outward_K_Current (dimensionless).
C RATES(25) is d/dt s_to in component Transient_Outward_K_Current (dimensionless).
C RATES(26) is d/dt ss_to in component Transient_Outward_K_Current (dimensionless).
C RATES(27) is d/dt r_ss in component Steady_State_K_Current (dimensionless).
C RATES(28) is d/dt s_ss in component Steady_State_K_Current (dimensionless).
C RATES(31) is d/dt Ca_nsr in component Other_SR_Fluxes_and_Concentrations (mM).
C RATES(29) is d/dt Ca_jsr in component Other_SR_Fluxes_and_Concentrations (mM).
C RATES(12) is d/dt Na_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(13) is d/dt K_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(14) is d/dt Ca_i in component Ion_Concentrations_and_Membrane_Potential (mM).
C RATES(18) is d/dt V_m in component Ion_Concentrations_and_Membrane_Potential (mV).
C RATES(30) is d/dt trel in component Ion_Concentrations_and_Membrane_Potential (second).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0
CONSTS(2) = 1
CONSTS(3) = 0.0132
CONSTS(4) = 3.83
CONSTS(5) = 0.285
CONSTS(6) = 0.062
CONSTS(7) = 33
CONSTS(8) = 1.1e-3
CONSTS(9) = 2.2e-3
CONSTS(10) = 3.6e-3
CONSTS(11) = 2.2e-3
CONSTS(12) = 16
CONSTS(13) = 0.8
CONSTS(14) = 1.21e3
CONSTS(15) = 49.7e-3
CONSTS(16) = 5e3
CONSTS(17) = 38.9e-3
CONSTS(18) = 0
CONSTS(19) = 0
CONSTS(20) = 0.2
CONSTS(21) = 8.5
CONSTS(22) = 7.3
CONSTS(23) = 5
CONSTS(24) = 1.03e3
CONSTS(25) = 315
CONSTS(26) = 860
CONSTS(27) = 1.3
CONSTS(28) = 0.4
CONSTS(29) = 44
CONSTS(30) = 30
CONSTS(31) = 0.59
CONSTS(32) = 0.025
CONSTS(33) = 0.18
CONSTS(34) = 9.14
CONSTS(35) = 1.64
CONSTS(36) = 4.375
CONSTS(37) = 0.2e-3
CONSTS(38) = 106
CONSTS(39) = 0.89
CONSTS(40) = 0.3
CONSTS(41) = 54
CONSTS(42) = 21
CONSTS(43) = 8.5
CONSTS(44) = 7
CONSTS(45) = 60
CONSTS(46) = 1
CONSTS(47) = 14
CONSTS(48) = 1
CONSTS(49) = 1e-3
CONSTS(50) = 10
CONSTS(51) = 25e-3
CONSTS(52) = 25e-3
CONSTS(53) = 25e-3
CONSTS(54) = 54
CONSTS(55) = 21
CONSTS(56) = 8.52
CONSTS(57) = 3
CONSTS(58) = 10.1
CONSTS(59) = 3
CONSTS(60) = 20.8e-6
CONSTS(61) = 9.88e-7
CONSTS(62) = 9.3e-8
CONSTS(63) = 1.534e-4
CONSTS(64) = 310
CONSTS(65) = 140
CONSTS(66) = 5.4
CONSTS(67) = 1.8
CONSTS(68) = 8
CONSTS(69) = 0.35
CONSTS(70) = 0.07
CONSTS(71) = 0.24
CONSTS(72) = 0.008
CONSTS(73) = 300
CONSTS(74) = 2e3
CONSTS(75) = 5187.5
CONSTS(76) = 10
CONSTS(77) = 1.7469e-8
CONSTS(78) = 3.234e-11
CONSTS(79) = 3e5
CONSTS(80) = -0.458
CONSTS(81) = 1483
CONSTS(82) = 87.5
CONSTS(83) = 1.38
CONSTS(84) = 0.1
CONSTS(85) = 0.35
CONSTS(86) = 1.1
CONSTS(87) = 10
CONSTS(88) = 1.5
CONSTS(89) = 1.15
CONSTS(90) = 0.5e-3
CONSTS(91) = 2.8e-3
CONSTS(92) = 1.18e-3
CONSTS(93) = 0
CONSTS(94) = 1.2e-3
CONSTS(95) = 4.7
CONSTS(96) = 3e-4
CONSTS(97) = 15
CONSTS(98) = 2e-3
CONSTS(99) = 2e-3
CONSTS(100) = 60e3
CONSTS(101) = 0.18e-3
CONSTS(102) = 0.8e-3
CONSTS(103) = 8.75
CONSTS(104) = 15
CONSTS(105) = 0.8
CONSTS(106) = 5.7e-4
CONSTS(107) = 0.07
CONSTS(108) = 0.05
CONSTS(109) = 0.07
CONSTS(110) = 0.5128e-3
CONSTS(111) = 2.38e-3
CONSTS(112) = 8.44e-4
STATES(1) = 0.01205
STATES(2) = 0
STATES(3) = 1.154e-3
STATES(4) = 0.02505
STATES(5) = 6.446e-4
STATES(6) = 0.02569
STATES(7) = 0.8453
STATES(8) = 4.105
STATES(9) = 0.0526
STATES(10) = 5.103e-3
STATES(11) = 5.841e-3
CONSTS(113) = 8314
CONSTS(114) = 96485
CONSTS(115) = 1
CONSTS(116) = 1
CONSTS(117) = 2
STATES(12) = 16
STATES(13) = 145
STATES(14) = 1.58e-4
STATES(15) = 1.4e-3
STATES(16) = 0.99
STATES(17) = 0.99
STATES(18) = -85.66
STATES(19) = 0
STATES(20) = 0
STATES(21) = 0.13
STATES(22) = 0.96
STATES(23) = 0.92
STATES(24) = 1.4e-3
STATES(25) = 1
STATES(26) = 0.613
STATES(27) = 198e-3
STATES(28) = 0.43
STATES(29) = 1.92
STATES(30) = 0.9
STATES(31) = 1.92
CONSTS(118) = -40
CONSTS(119) = -10
CONSTS(120) = 0.02
CONSTS(121) = 0.961300
CONSTS(122) = 0.204100
CONSTS(123) = 0.233600
CONSTS(124) = - 40.0000
CONSTS(125) = (CONSTS(114)/CONSTS(113))/CONSTS(64)
CONSTS(126) = 2.10000
CONSTS(127) = (EXP(CONSTS(65)/67.3000) - 1.00000)/7.00000
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(18)+87.5000)/10.3000))
RATES(28) = (ALGBRC(11) - STATES(28))/CONSTS(126)
ALGBRC(3) = ( 0.320000*(STATES(18)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(18)+47.1300)))
ALGBRC(13) = 0.0800000*EXP(- STATES(18)/11.0000)
RATES(15) = 1000.00*( ALGBRC(3)*(1.00000 - STATES(15)) - ALGBRC(13)*STATES(15))
ALGBRC(4) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, 0.135000*EXP((80.0000+STATES(18))/- 6.80000))
ALGBRC(14) = TERNRY(STATES(18).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(- (STATES(18)+10.6600)/11.1000))), 3.56000*EXP( 0.0790000*STATES(18))+ 310000.*EXP( 0.350000*STATES(18)))
RATES(16) = 1000.00*( ALGBRC(4)*(1.00000 - STATES(16)) - ALGBRC(14)*STATES(16))
ALGBRC(5) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(18)) - 3.47400e-05*EXP( - 0.0439100*STATES(18)))*(STATES(18)+37.7800))/(1.00000+EXP( 0.311000*(STATES(18)+79.2300))))
ALGBRC(15) = TERNRY(STATES(18).GE.- 40.0000, ( 0.300000*EXP( - 2.57500e-07*STATES(18)))/(1.00000+EXP( - 0.100000*(STATES(18)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(18)))/(1.00000+EXP( - 0.137800*(STATES(18)+40.1400))))
RATES(17) = 1000.00*( ALGBRC(5)*(1.00000 - STATES(17)) - ALGBRC(15)*STATES(17))
ALGBRC(6) = 400.000*EXP((STATES(18)+2.00000)/10.0000)
ALGBRC(16) = 50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
RATES(19) = ALGBRC(6)*(1.00000 - STATES(19)) - ALGBRC(16)*STATES(19)
RATES(20) = 2.00000*ALGBRC(6)*(1.00000 - STATES(20)) - (ALGBRC(16)/2.00000)*STATES(20)
ALGBRC(2) = STATES(10)/CONSTS(51)
ALGBRC(18) = CONSTS(73)*(( 0.375000*ALGBRC(2))/CONSTS(122)+0.625000)
RATES(21) = ALGBRC(18)*(1.00000 - STATES(21)) - CONSTS(74)*STATES(21)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(18)+55.0000)/7.50000))+0.100000/(1.00000+EXP((- STATES(18)+21.0000)/6.00000))
ALGBRC(17) = 0.0200000+0.300000/(1.00000+EXP((STATES(18)+30.0000)/9.50000))
RATES(22) = (ALGBRC(7) - STATES(22))/ALGBRC(17)
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(18)+10.6000)/- 11.4200))
ALGBRC(19) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
RATES(24) = (ALGBRC(8) - STATES(24))/ALGBRC(19)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(18)+43.5000)/6.88410))
ALGBRC(20) = 0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
RATES(25) = (ALGBRC(9) - STATES(25))/ALGBRC(20)
ALGBRC(21) = 3.70000*EXP( - 1.00000*(STATES(18)+70.0000)/30.0000 ** 2.00000)+0.0350000
RATES(26) = (ALGBRC(9) - STATES(26))/ALGBRC(21)
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(18)+11.5000)/- 11.8200))
ALGBRC(22) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
RATES(27) = (ALGBRC(10) - STATES(27))/ALGBRC(22)
ALGBRC(25) = CONSTS(13)*STATES(4)
ALGBRC(27) = CONSTS(14)*STATES(5)*STATES(6)
RATES(5) = ALGBRC(25) - ALGBRC(27)
ALGBRC(24) = CONSTS(75)*STATES(14)
ALGBRC(26) = ALGBRC(24)*(1.00000 - STATES(19) ** 4.00000+ 2.00000*STATES(19)*1.00000 - STATES(19) ** 3.00000+ 4.00000*STATES(19) ** 2.00000*1.00000 - STATES(19) ** 2.00000+ 8.00000*STATES(19) ** 3.00000*(1.00000 - STATES(19))+ 16.0000*STATES(19) ** 4.00000*(1.00000 - ALGBRC(18)/CONSTS(74)))
ALGBRC(28) = CONSTS(76)*(1.00000 - STATES(20) ** 4.00000+ 0.500000*STATES(20)*1.00000 - STATES(20) ** 3.00000+ 0.250000*STATES(20) ** 2.00000*1.00000 - STATES(20) ** 2.00000+ 0.125000*STATES(20) ** 3.00000*(1.00000 - STATES(20))+ 0.0625000*STATES(20) ** 4.00000)
RATES(23) = ALGBRC(28)*(1.00000 - STATES(23)) - ALGBRC(26)*STATES(23)
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
ALGBRC(80) = CONSTS(8)*(ALGBRC(62)+ALGBRC(63))
ALGBRC(12) = CONSTS(9)*STATES(2)
RATES(2) = ALGBRC(80) - ALGBRC(12)
ALGBRC(81) = CONSTS(12)*(ALGBRC(64)+ALGBRC(63))
RATES(4) = ALGBRC(81) - ALGBRC(25)
RATES(6) = ALGBRC(81) - ALGBRC(27)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
ALGBRC(48) = ( CONSTS(21)*ALGBRC(44)*CONSTS(16))/(CONSTS(25)+CONSTS(16))
ALGBRC(46) = ( CONSTS(20)*ALGBRC(43)*CONSTS(16))/(CONSTS(24)+CONSTS(16))
CALL minimize(minfunc_3, CONSTS, VARIABLES, <VAR>)
ALGBRC(53) = ( CONSTS(22)*ALGBRC(51)*CONSTS(16))/(CONSTS(26)+CONSTS(16))
ALGBRC(82) = ( CONSTS(23)*CONSTS(128)*ALGBRC(69))/(CONSTS(27)+ALGBRC(69))
RATES(7) = (ALGBRC(46)+ALGBRC(48)+ALGBRC(53)) - ALGBRC(82)
ALGBRC(29) = CONSTS(38) - STATES(8)
ALGBRC(83) = ( CONSTS(41)*ALGBRC(68)*ALGBRC(29))/(CONSTS(42)+ALGBRC(29))
CALL minimize(minfunc_4, CONSTS, VARIABLES, <VAR>)
ALGBRC(60) = ( CONSTS(43)*ALGBRC(57)*STATES(8))/(CONSTS(44)+STATES(8))
RATES(8) = ALGBRC(83) - ALGBRC(60)
ALGBRC(30) = CONSTS(40) - STATES(9)
ALGBRC(84) = ( CONSTS(45)*ALGBRC(68)*ALGBRC(30))/(CONSTS(46)+ALGBRC(30))
ALGBRC(32) = ( CONSTS(47)*STATES(9))/(CONSTS(48)+STATES(9))
RATES(9) = ALGBRC(84) - ALGBRC(32)
ALGBRC(34) = CONSTS(51) - STATES(10)
ALGBRC(85) = ( CONSTS(50)*CONSTS(54)*ALGBRC(78)*ALGBRC(34))/(CONSTS(55)+ CONSTS(50)*ALGBRC(34))
ALGBRC(36) = ( CONSTS(50)*CONSTS(58)*CONSTS(53)*STATES(10))/(CONSTS(59)+ CONSTS(50)*STATES(10))
RATES(10) = ALGBRC(85) - ALGBRC(36)
ALGBRC(38) = CONSTS(51) - STATES(11)
ALGBRC(86) = ( CONSTS(50)*CONSTS(54)*ALGBRC(78)*ALGBRC(38))/(CONSTS(55)+ CONSTS(50)*ALGBRC(38))
ALGBRC(40) = ( CONSTS(50)*CONSTS(56)*CONSTS(52)*STATES(11))/(CONSTS(57)+ CONSTS(50)*STATES(11))
RATES(11) = ALGBRC(86) - ALGBRC(40)
ALGBRC(88) = CONSTS(10)*ALGBRC(68)*STATES(1)
ALGBRC(23) = CONSTS(11)*STATES(3)
RATES(1) = (ALGBRC(12) - ALGBRC(80))+(ALGBRC(23) - ALGBRC(88))
RATES(3) = ALGBRC(88) - ALGBRC(23)
ALGBRC(35) = (1.00000/( CONSTS(125)*CONSTS(115)))*log(CONSTS(65)/STATES(12))
ALGBRC(41) = CONSTS(68)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(35))
ALGBRC(93) = EXP( CONSTS(85)*STATES(18)*CONSTS(125))*STATES(12) ** 3.00000*CONSTS(67)
ALGBRC(94) = EXP( (CONSTS(85) - 1.00000)*STATES(18)*CONSTS(125))*CONSTS(65) ** 3.00000*STATES(14)
ALGBRC(95) = (CONSTS(81)/( (CONSTS(82) ** 3.00000+CONSTS(65) ** 3.00000)*(CONSTS(83)+CONSTS(67))*(1.00000+ CONSTS(84)*EXP( (CONSTS(85) - 1.00000)*STATES(18)*CONSTS(125)))))*(ALGBRC(93) - ALGBRC(94))
ALGBRC(96) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(125))+ 0.0365000*CONSTS(127)*EXP( - STATES(18)*CONSTS(125)))
ALGBRC(97) = ( (( CONSTS(86)*ALGBRC(96))/(1.00000+CONSTS(87)/STATES(12) ** 1.50000))*CONSTS(66))/(CONSTS(66)+CONSTS(88))
ALGBRC(100) = CONSTS(92)*(STATES(18) - ALGBRC(35))
ALGBRC(101) = ALGBRC(41)+ALGBRC(100)+ 3.00000*ALGBRC(95)+ 3.00000*ALGBRC(97)
RATES(12) = ( - 1000.00*ALGBRC(101)*CONSTS(63))/( CONSTS(60)*CONSTS(115)*CONSTS(114))
ALGBRC(47) = ( CONSTS(78)*STATES(18)*CONSTS(114)*CONSTS(125)*( STATES(13)*EXP( STATES(18)*CONSTS(125)) - CONSTS(66)))/(EXP( STATES(18)*CONSTS(125)) - 1.00000)
ALGBRC(33) = STATES(11)/CONSTS(51)
ALGBRC(49) = 0.500000*(( 0.400000*ALGBRC(33))/CONSTS(123)+0.600000)
ALGBRC(45) = ( CONSTS(77)*4.00000*STATES(18)*CONSTS(114)*CONSTS(125)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(125)) - 0.341000*CONSTS(67)))/(EXP( 2.00000*STATES(18)*CONSTS(125)) - 1.00000)
ALGBRC(52) = ALGBRC(45)*CONSTS(79)*ALGBRC(49)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(54) = (ALGBRC(47)/(1.00000+ALGBRC(52)/CONSTS(80)))*CONSTS(79)*ALGBRC(49)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(37) = (1.00000/( CONSTS(125)*CONSTS(116)))*log(CONSTS(66)/STATES(13))
ALGBRC(59) = CONSTS(69)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(37))
ALGBRC(61) = CONSTS(70)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(37))
ALGBRC(79) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(37)) - 59.2150)))
ALGBRC(87) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(37)))+EXP( 0.0617500*((STATES(18) - ALGBRC(37)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(37))+4.75300)))
ALGBRC(89) = ALGBRC(79)/(ALGBRC(79)+ALGBRC(87))
ALGBRC(90) = CONSTS(71)* (CONSTS(66)/5.40000) ** (1.0 / 2)*ALGBRC(89)*(STATES(18) - ALGBRC(37))
ALGBRC(91) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
ALGBRC(92) = CONSTS(72)*ALGBRC(91)*(STATES(18) - ALGBRC(37))
ALGBRC(102) = ((ALGBRC(59)+ALGBRC(61)+ALGBRC(90)+ALGBRC(92)) - 2.00000*ALGBRC(97))+ALGBRC(54)
RATES(13) = ( - 1000.00*ALGBRC(102)*CONSTS(63))/( CONSTS(60)*CONSTS(116)*CONSTS(114))
ALGBRC(98) = ( CONSTS(89)*STATES(14))/(CONSTS(90)+STATES(14))
ALGBRC(39) = (1.00000/( CONSTS(125)*CONSTS(117)))*log(CONSTS(67)/STATES(14))
ALGBRC(99) = CONSTS(91)*(STATES(18) - ALGBRC(39))
ALGBRC(103) = (ALGBRC(52)+ALGBRC(99)+ALGBRC(98)) - 2.00000*ALGBRC(95)
ALGBRC(105) = TERNRY(VOI.GT.59.1000.AND.VOI.LT.59.5000, CONSTS(119), CONSTS(118))
ALGBRC(107) = TERNRY( INT(MOD(VOI+0.900000, 1.00000)).LT.0.00500000, 10.0000, 0.00000)
ALGBRC(109) = TERNRY(CONSTS(1).EQ.0.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000, ALGBRC(107), (ALGBRC(105) - STATES(18))/0.0200000)
RATES(18) = - 1000.00*((ALGBRC(103)+ALGBRC(102)+ALGBRC(101)) - ALGBRC(109))
RATES(30) = TERNRY( - 1000.00*((ALGBRC(103)+ALGBRC(102)+ALGBRC(101)) - ALGBRC(109)).GT.30000.0, 1.00000 - 10000.0*STATES(30), 1.00000)
ALGBRC(31) = ALGBRC(29)/CONSTS(38)
ALGBRC(112) = ( CONSTS(96)*(1.00000+ 2.00000*ALGBRC(31)))/(1.00000+ 2.00000*CONSTS(121))
ALGBRC(113) = ( CONSTS(95)*STATES(14) ** 2.00000)/(ALGBRC(112) ** 2.00000+STATES(14) ** 2.00000)
ALGBRC(114) = ( CONSTS(95)*STATES(31))/CONSTS(97)
ALGBRC(115) = (STATES(31) - STATES(29))/CONSTS(106)
RATES(31) = (ALGBRC(113) - ALGBRC(114)) - ( ALGBRC(115)*CONSTS(62))/CONSTS(61)
ALGBRC(104) = STATES(30)+0.00200000
ALGBRC(106) = 1.00000 - EXP(- ALGBRC(104)/CONSTS(98))
ALGBRC(108) = EXP(- ALGBRC(104)/CONSTS(99))
ALGBRC(110) = CONSTS(100)/(1.00000+EXP((ALGBRC(103)+5.00000)/0.900000))
ALGBRC(111) = ALGBRC(110)*ALGBRC(106)*ALGBRC(108)*(STATES(29) - STATES(14))
ALGBRC(117) = 1.00000/(1.00000+( CONSTS(104)*CONSTS(105))/CONSTS(105)+STATES(29) ** 2.00000)
RATES(29) = ALGBRC(117)*(ALGBRC(115) - ALGBRC(111))
ALGBRC(116) = ( CONSTS(107)*CONSTS(110))/CONSTS(110)+STATES(14) ** 2.00000
ALGBRC(118) = ( CONSTS(108)*CONSTS(111))/CONSTS(111)+STATES(14) ** 2.00000
ALGBRC(120) = ( CONSTS(109)*CONSTS(112))/CONSTS(112)+STATES(14) ** 2.00000
ALGBRC(121) = 1.00000/(1.00000+ALGBRC(118)+ALGBRC(116)+ALGBRC(116)+ALGBRC(120))
RATES(14) = - ALGBRC(121)*((( 1000.00*ALGBRC(103)*CONSTS(63))/( CONSTS(60)*CONSTS(117)*CONSTS(114))+( (ALGBRC(113) - ALGBRC(114))*CONSTS(61))/CONSTS(60)) - ( ALGBRC(111)*CONSTS(62))/CONSTS(60))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(18)+87.5000)/10.3000))
ALGBRC(3) = ( 0.320000*(STATES(18)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(18)+47.1300)))
ALGBRC(13) = 0.0800000*EXP(- STATES(18)/11.0000)
ALGBRC(4) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, 0.135000*EXP((80.0000+STATES(18))/- 6.80000))
ALGBRC(14) = TERNRY(STATES(18).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(- (STATES(18)+10.6600)/11.1000))), 3.56000*EXP( 0.0790000*STATES(18))+ 310000.*EXP( 0.350000*STATES(18)))
ALGBRC(5) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(18)) - 3.47400e-05*EXP( - 0.0439100*STATES(18)))*(STATES(18)+37.7800))/(1.00000+EXP( 0.311000*(STATES(18)+79.2300))))
ALGBRC(15) = TERNRY(STATES(18).GE.- 40.0000, ( 0.300000*EXP( - 2.57500e-07*STATES(18)))/(1.00000+EXP( - 0.100000*(STATES(18)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(18)))/(1.00000+EXP( - 0.137800*(STATES(18)+40.1400))))
ALGBRC(6) = 400.000*EXP((STATES(18)+2.00000)/10.0000)
ALGBRC(16) = 50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
ALGBRC(2) = STATES(10)/CONSTS(51)
ALGBRC(18) = CONSTS(73)*(( 0.375000*ALGBRC(2))/CONSTS(122)+0.625000)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(18)+55.0000)/7.50000))+0.100000/(1.00000+EXP((- STATES(18)+21.0000)/6.00000))
ALGBRC(17) = 0.0200000+0.300000/(1.00000+EXP((STATES(18)+30.0000)/9.50000))
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(18)+10.6000)/- 11.4200))
ALGBRC(19) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(18)+43.5000)/6.88410))
ALGBRC(20) = 0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
ALGBRC(21) = 3.70000*EXP( - 1.00000*(STATES(18)+70.0000)/30.0000 ** 2.00000)+0.0350000
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(18)+11.5000)/- 11.8200))
ALGBRC(22) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
ALGBRC(25) = CONSTS(13)*STATES(4)
ALGBRC(27) = CONSTS(14)*STATES(5)*STATES(6)
ALGBRC(24) = CONSTS(75)*STATES(14)
ALGBRC(26) = ALGBRC(24)*(1.00000 - STATES(19) ** 4.00000+ 2.00000*STATES(19)*1.00000 - STATES(19) ** 3.00000+ 4.00000*STATES(19) ** 2.00000*1.00000 - STATES(19) ** 2.00000+ 8.00000*STATES(19) ** 3.00000*(1.00000 - STATES(19))+ 16.0000*STATES(19) ** 4.00000*(1.00000 - ALGBRC(18)/CONSTS(74)))
ALGBRC(28) = CONSTS(76)*(1.00000 - STATES(20) ** 4.00000+ 0.500000*STATES(20)*1.00000 - STATES(20) ** 3.00000+ 0.250000*STATES(20) ** 2.00000*1.00000 - STATES(20) ** 2.00000+ 0.125000*STATES(20) ** 3.00000*(1.00000 - STATES(20))+ 0.0625000*STATES(20) ** 4.00000)
ALGBRC(80) = CONSTS(8)*(ALGBRC(62)+ALGBRC(63))
ALGBRC(12) = CONSTS(9)*STATES(2)
ALGBRC(81) = CONSTS(12)*(ALGBRC(64)+ALGBRC(63))
ALGBRC(48) = ( CONSTS(21)*ALGBRC(44)*CONSTS(16))/(CONSTS(25)+CONSTS(16))
ALGBRC(46) = ( CONSTS(20)*ALGBRC(43)*CONSTS(16))/(CONSTS(24)+CONSTS(16))
ALGBRC(53) = ( CONSTS(22)*ALGBRC(51)*CONSTS(16))/(CONSTS(26)+CONSTS(16))
ALGBRC(82) = ( CONSTS(23)*CONSTS(128)*ALGBRC(69))/(CONSTS(27)+ALGBRC(69))
ALGBRC(29) = CONSTS(38) - STATES(8)
ALGBRC(83) = ( CONSTS(41)*ALGBRC(68)*ALGBRC(29))/(CONSTS(42)+ALGBRC(29))
ALGBRC(60) = ( CONSTS(43)*ALGBRC(57)*STATES(8))/(CONSTS(44)+STATES(8))
ALGBRC(30) = CONSTS(40) - STATES(9)
ALGBRC(84) = ( CONSTS(45)*ALGBRC(68)*ALGBRC(30))/(CONSTS(46)+ALGBRC(30))
ALGBRC(32) = ( CONSTS(47)*STATES(9))/(CONSTS(48)+STATES(9))
ALGBRC(34) = CONSTS(51) - STATES(10)
ALGBRC(85) = ( CONSTS(50)*CONSTS(54)*ALGBRC(78)*ALGBRC(34))/(CONSTS(55)+ CONSTS(50)*ALGBRC(34))
ALGBRC(36) = ( CONSTS(50)*CONSTS(58)*CONSTS(53)*STATES(10))/(CONSTS(59)+ CONSTS(50)*STATES(10))
ALGBRC(38) = CONSTS(51) - STATES(11)
ALGBRC(86) = ( CONSTS(50)*CONSTS(54)*ALGBRC(78)*ALGBRC(38))/(CONSTS(55)+ CONSTS(50)*ALGBRC(38))
ALGBRC(40) = ( CONSTS(50)*CONSTS(56)*CONSTS(52)*STATES(11))/(CONSTS(57)+ CONSTS(50)*STATES(11))
ALGBRC(88) = CONSTS(10)*ALGBRC(68)*STATES(1)
ALGBRC(23) = CONSTS(11)*STATES(3)
ALGBRC(35) = (1.00000/( CONSTS(125)*CONSTS(115)))*log(CONSTS(65)/STATES(12))
ALGBRC(41) = CONSTS(68)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(35))
ALGBRC(93) = EXP( CONSTS(85)*STATES(18)*CONSTS(125))*STATES(12) ** 3.00000*CONSTS(67)
ALGBRC(94) = EXP( (CONSTS(85) - 1.00000)*STATES(18)*CONSTS(125))*CONSTS(65) ** 3.00000*STATES(14)
ALGBRC(95) = (CONSTS(81)/( (CONSTS(82) ** 3.00000+CONSTS(65) ** 3.00000)*(CONSTS(83)+CONSTS(67))*(1.00000+ CONSTS(84)*EXP( (CONSTS(85) - 1.00000)*STATES(18)*CONSTS(125)))))*(ALGBRC(93) - ALGBRC(94))
ALGBRC(96) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(125))+ 0.0365000*CONSTS(127)*EXP( - STATES(18)*CONSTS(125)))
ALGBRC(97) = ( (( CONSTS(86)*ALGBRC(96))/(1.00000+CONSTS(87)/STATES(12) ** 1.50000))*CONSTS(66))/(CONSTS(66)+CONSTS(88))
ALGBRC(100) = CONSTS(92)*(STATES(18) - ALGBRC(35))
ALGBRC(101) = ALGBRC(41)+ALGBRC(100)+ 3.00000*ALGBRC(95)+ 3.00000*ALGBRC(97)
ALGBRC(47) = ( CONSTS(78)*STATES(18)*CONSTS(114)*CONSTS(125)*( STATES(13)*EXP( STATES(18)*CONSTS(125)) - CONSTS(66)))/(EXP( STATES(18)*CONSTS(125)) - 1.00000)
ALGBRC(33) = STATES(11)/CONSTS(51)
ALGBRC(49) = 0.500000*(( 0.400000*ALGBRC(33))/CONSTS(123)+0.600000)
ALGBRC(45) = ( CONSTS(77)*4.00000*STATES(18)*CONSTS(114)*CONSTS(125)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(125)) - 0.341000*CONSTS(67)))/(EXP( 2.00000*STATES(18)*CONSTS(125)) - 1.00000)
ALGBRC(52) = ALGBRC(45)*CONSTS(79)*ALGBRC(49)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(54) = (ALGBRC(47)/(1.00000+ALGBRC(52)/CONSTS(80)))*CONSTS(79)*ALGBRC(49)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(37) = (1.00000/( CONSTS(125)*CONSTS(116)))*log(CONSTS(66)/STATES(13))
ALGBRC(59) = CONSTS(69)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(37))
ALGBRC(61) = CONSTS(70)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(37))
ALGBRC(79) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(37)) - 59.2150)))
ALGBRC(87) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(37)))+EXP( 0.0617500*((STATES(18) - ALGBRC(37)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(37))+4.75300)))
ALGBRC(89) = ALGBRC(79)/(ALGBRC(79)+ALGBRC(87))
ALGBRC(90) = CONSTS(71)* (CONSTS(66)/5.40000) ** (1.0 / 2)*ALGBRC(89)*(STATES(18) - ALGBRC(37))
ALGBRC(91) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
ALGBRC(92) = CONSTS(72)*ALGBRC(91)*(STATES(18) - ALGBRC(37))
ALGBRC(102) = ((ALGBRC(59)+ALGBRC(61)+ALGBRC(90)+ALGBRC(92)) - 2.00000*ALGBRC(97))+ALGBRC(54)
ALGBRC(98) = ( CONSTS(89)*STATES(14))/(CONSTS(90)+STATES(14))
ALGBRC(39) = (1.00000/( CONSTS(125)*CONSTS(117)))*log(CONSTS(67)/STATES(14))
ALGBRC(99) = CONSTS(91)*(STATES(18) - ALGBRC(39))
ALGBRC(103) = (ALGBRC(52)+ALGBRC(99)+ALGBRC(98)) - 2.00000*ALGBRC(95)
ALGBRC(105) = TERNRY(VOI.GT.59.1000.AND.VOI.LT.59.5000, CONSTS(119), CONSTS(118))
ALGBRC(107) = TERNRY( INT(MOD(VOI+0.900000, 1.00000)).LT.0.00500000, 10.0000, 0.00000)
ALGBRC(109) = TERNRY(CONSTS(1).EQ.0.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000, ALGBRC(107), (ALGBRC(105) - STATES(18))/0.0200000)
ALGBRC(31) = ALGBRC(29)/CONSTS(38)
ALGBRC(112) = ( CONSTS(96)*(1.00000+ 2.00000*ALGBRC(31)))/(1.00000+ 2.00000*CONSTS(121))
ALGBRC(113) = ( CONSTS(95)*STATES(14) ** 2.00000)/(ALGBRC(112) ** 2.00000+STATES(14) ** 2.00000)
ALGBRC(114) = ( CONSTS(95)*STATES(31))/CONSTS(97)
ALGBRC(115) = (STATES(31) - STATES(29))/CONSTS(106)
ALGBRC(104) = STATES(30)+0.00200000
ALGBRC(106) = 1.00000 - EXP(- ALGBRC(104)/CONSTS(98))
ALGBRC(108) = EXP(- ALGBRC(104)/CONSTS(99))
ALGBRC(110) = CONSTS(100)/(1.00000+EXP((ALGBRC(103)+5.00000)/0.900000))
ALGBRC(111) = ALGBRC(110)*ALGBRC(106)*ALGBRC(108)*(STATES(29) - STATES(14))
ALGBRC(117) = 1.00000/(1.00000+( CONSTS(104)*CONSTS(105))/CONSTS(105)+STATES(29) ** 2.00000)
ALGBRC(116) = ( CONSTS(107)*CONSTS(110))/CONSTS(110)+STATES(14) ** 2.00000
ALGBRC(118) = ( CONSTS(108)*CONSTS(111))/CONSTS(111)+STATES(14) ** 2.00000
ALGBRC(120) = ( CONSTS(109)*CONSTS(112))/CONSTS(112)+STATES(14) ** 2.00000
ALGBRC(121) = 1.00000/(1.00000+ALGBRC(118)+ALGBRC(116)+ALGBRC(116)+ALGBRC(120))
ALGBRC(1) = STATES(8)/CONSTS(38)
ALGBRC(58) = ALGBRC(52)+ALGBRC(54)
ALGBRC(119) = 1000.00*(( (STATES(29)+STATES(29)/ALGBRC(117))*CONSTS(62))/CONSTS(60)+( STATES(31)*CONSTS(61))/CONSTS(60))
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0[1] = abs(CONSTS(130) - ( CONSTS(128)*CONSTS(129))/CONSTS(30))
minfunc_0[2] = abs(CONSTS(128) - (CONSTS(17) - CONSTS(130)))
minfunc_0[3] = abs(CONSTS(129) - (CONSTS(18) - CONSTS(130)))
RETURN
END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_1[1] = abs(ALGBRC(62) - ( ALGBRC(65)*ALGBRC(66))/CONSTS(5))
minfunc_1[2] = abs(ALGBRC(63) - ( ALGBRC(62)*ALGBRC(67))/CONSTS(6))
minfunc_1[3] = abs(ALGBRC(64) - ( ALGBRC(66)*ALGBRC(67))/CONSTS(7))
minfunc_1[4] = abs(ALGBRC(65) - ((CONSTS(2) - ALGBRC(62)) - ALGBRC(63)))
minfunc_1[5] = abs(ALGBRC(66) - (((STATES(1) - ALGBRC(62)) - ALGBRC(63)) - ALGBRC(64)))
minfunc_1[6] = abs(ALGBRC(67) - ((CONSTS(4) - ALGBRC(63)) - ALGBRC(64)))
minfunc_1[7] = abs(ALGBRC(70) - ( CONSTS(33)*CONSTS(37))/(CONSTS(37)+ALGBRC(68)+ALGBRC(78)))
minfunc_1[8] = abs(ALGBRC(71) - (ALGBRC(68)/CONSTS(36))*ALGBRC(68)*(1.00000+ALGBRC(70)/CONSTS(37)))
minfunc_1[9] = abs(ALGBRC(72) - ALGBRC(68)*(1.00000+ALGBRC(70)/CONSTS(37)))
minfunc_1[10] = abs(ALGBRC(73) - (ALGBRC(78)/CONSTS(36))*ALGBRC(78)*(1.00000+ALGBRC(70)/CONSTS(37)))
minfunc_1[11] = abs(ALGBRC(74) - ALGBRC(78)*(1.00000+ALGBRC(70)/CONSTS(37)))
minfunc_1[12] = abs(ALGBRC(75) - (CONSTS(34)/ALGBRC(69))*ALGBRC(71))
minfunc_1[13] = abs(ALGBRC(76) - (CONSTS(34)/ALGBRC(69))*ALGBRC(73))
minfunc_1[14] = abs(ALGBRC(77) - (( CONSTS(34)*CONSTS(35))/CONSTS(36)+( CONSTS(34)*ALGBRC(69))/CONSTS(36)+ALGBRC(69) ** 2.00000/CONSTS(36)))
minfunc_1[15] = abs(ALGBRC(69) - ((STATES(7) - (ALGBRC(75)+ 2.00000*ALGBRC(71)+ 2.00000*ALGBRC(72))) - (ALGBRC(76)+ 2.00000*ALGBRC(73)+ 2.00000*ALGBRC(74))))
minfunc_1[16] = abs(0.00000 - ( 2.00000*CONSTS(31)*ALGBRC(69) ** 2.00000 - ALGBRC(68)*(1.00000+ALGBRC(70)/CONSTS(37))*( ALGBRC(77)*ALGBRC(68)+ALGBRC(69) ** 2.00000)))
minfunc_1[17] = abs(0.00000 - ( 2.00000*CONSTS(32)*ALGBRC(69) ** 2.00000 - ALGBRC(78)*(1.00000+ALGBRC(70)/CONSTS(37))*( ALGBRC(77)*ALGBRC(78)+ALGBRC(69) ** 2.00000)))
RETURN
END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_2[1] = abs(ALGBRC(44) - ( ALGBRC(42)*ALGBRC(43))/CONSTS(28))
minfunc_2[2] = abs(ALGBRC(42) - (STATES(4) - ALGBRC(44)))
minfunc_2[3] = abs(ALGBRC(43) - (CONSTS(15) - ALGBRC(44)))
RETURN
END
REAL FUNCTION minfunc_3(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_3[1] = abs(ALGBRC(51) - ( ALGBRC(50)*ALGBRC(43))/CONSTS(29))
minfunc_3[2] = abs(ALGBRC(50) - (CONSTS(19) - ALGBRC(51)))
RETURN
END
REAL FUNCTION minfunc_4(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_4[1] = abs(ALGBRC(55) - ( ALGBRC(56)*ALGBRC(57))/CONSTS(49))
minfunc_4[2] = abs(ALGBRC(56) - (STATES(9) - ALGBRC(55)))
minfunc_4[3] = abs(ALGBRC(57) - (CONSTS(39) - ALGBRC(55)))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END