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 119 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 L_totmax in component b1_AR_Gs_parameters (uM).
C CONSTS(2) is sum_b1_AR in component b1_AR_Gs_parameters (uM).
C CONSTS(3) is Gs_tot in component b1_AR_Gs_parameters (uM).
C CONSTS(4) is Kl in component b1_AR_Gs_parameters (uM).
C CONSTS(5) is Kr in component b1_AR_Gs_parameters (uM).
C CONSTS(6) is Kc in component b1_AR_Gs_parameters (uM).
C CONSTS(7) is k_bar_kp in component b1_AR_Gs_parameters (per_sec).
C CONSTS(8) is k_bar_km in component b1_AR_Gs_parameters (per_sec).
C CONSTS(9) is k_p_kap in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(10) is k_p_kam in component b1_AR_Gs_parameters (per_sec).
C CONSTS(11) is k_g_act in component b1_AR_Gs_parameters (per_sec).
C CONSTS(12) is k_hyd in component b1_AR_Gs_parameters (per_sec).
C CONSTS(13) is k_reassoc in component b1_AR_Gs_parameters (per_uM_per_sec).
C CONSTS(14) is AC_tot in component cAMP_parameters (uM).
C CONSTS(15) is ATP in component cAMP_parameters (uM).
C CONSTS(16) is PDE_tot in component cAMP_parameters (uM).
C CONSTS(17) is IBMX_tot in component cAMP_parameters (uM).
C CONSTS(18) is Fsk_tot in component cAMP_parameters (uM).
C CONSTS(19) is k_ac_basal in component cAMP_parameters (per_sec).
C CONSTS(20) is k_ac_gsa in component cAMP_parameters (per_sec).
C CONSTS(21) is k_ac_fsk in component cAMP_parameters (per_sec).
C CONSTS(22) is k_pde in component cAMP_parameters (per_sec).
C CONSTS(23) is Km_basal in component cAMP_parameters (uM).
C CONSTS(24) is Km_gsa in component cAMP_parameters (uM).
C CONSTS(25) is Km_fsk in component cAMP_parameters (uM).
C CONSTS(26) is Km_pde in component cAMP_parameters (uM).
C CONSTS(27) is K_gsa in component cAMP_parameters (uM).
C CONSTS(28) is K_fsk in component cAMP_parameters (uM).
C CONSTS(29) is Ki_ibmx in component cAMP_parameters (uM).
C CONSTS(30) is PKAI_tot in component PKA_parameters (uM).
C CONSTS(31) is PKAII_tot in component PKA_parameters (uM).
C CONSTS(32) is PKI_tot in component PKA_parameters (uM).
C CONSTS(33) is K_a in component PKA_parameters (uM).
C CONSTS(34) is K_b in component PKA_parameters (uM).
C CONSTS(35) is K_d in component PKA_parameters (uM).
C CONSTS(36) is Ki_pki in component PKA_parameters (uM).
C CONSTS(37) is PLB_tot in component PLB_parameters (uM).
C CONSTS(38) is PP1_tot in component PLB_parameters (uM).
C CONSTS(39) is Inhib1_tot in component PLB_parameters (uM).
C CONSTS(40) is k_pka_plb in component PLB_parameters (per_sec).
C CONSTS(41) is Km_pka_plb in component PLB_parameters (uM).
C CONSTS(42) is k_pp1_plb in component PLB_parameters (per_sec).
C CONSTS(43) is Km_pp1_plb in component PLB_parameters (uM).
C CONSTS(44) is k_pka_i1 in component PLB_parameters (per_sec).
C CONSTS(45) is Km_pka_i1 in component PLB_parameters (uM).
C CONSTS(46) is Vmax_pp2a_i1 in component PLB_parameters (uM_per_sec).
C CONSTS(47) is Km_pp2a_i1 in component PLB_parameters (uM).
C CONSTS(48) is Ki_inhib1 in component PLB_parameters (uM).
C CONSTS(49) is epsilon in component LCC_parameters (dimensionless).
C CONSTS(50) is LCC_tot in component LCC_parameters (uM).
C CONSTS(51) is PP1_lcc_tot in component LCC_parameters (uM).
C CONSTS(52) is PP2A_lcc_tot in component LCC_parameters (uM).
C CONSTS(53) is k_pka_lcc in component LCC_parameters (per_sec).
C CONSTS(54) is Km_pka_lcc in component LCC_parameters (uM).
C CONSTS(55) is k_pp1_lcc in component LCC_parameters (per_sec).
C CONSTS(56) is Km_pp1_lcc in component LCC_parameters (uM).
C CONSTS(57) is k_pp2a_lcc in component LCC_parameters (per_sec).
C CONSTS(58) is Km_pp2a_lcc in component LCC_parameters (uM).
C CONSTS(59) is V_myo in component EC_Coupling_Parameters (uL).
C CONSTS(60) is V_nsr in component EC_Coupling_Parameters (uL).
C CONSTS(61) is V_jsr in component EC_Coupling_Parameters (uL).
C CONSTS(62) is A_Cap in component EC_Coupling_Parameters (cm2).
C CONSTS(63) is Temp in component EC_Coupling_Parameters (kelvin).
C CONSTS(64) is Na_ext in component EC_Coupling_Parameters (mM).
C CONSTS(65) is K_ext in component EC_Coupling_Parameters (mM).
C CONSTS(66) is Ca_ext in component EC_Coupling_Parameters (mM).
C CONSTS(67) is G_Na in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(68) is G_to in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(69) is G_ss in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(70) is G_Ki_bar in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(71) is G_Kp in component EC_Coupling_Parameters (mS_per_uF).
C CONSTS(72) is f in component EC_Coupling_Parameters (per_sec).
C CONSTS(73) is g in component EC_Coupling_Parameters (per_sec).
C CONSTS(74) is gamma_o in component EC_Coupling_Parameters (per_mM_per_sec).
C CONSTS(75) is omega in component EC_Coupling_Parameters (per_sec).
C CONSTS(76) is p_Ca in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(77) is p_K in component EC_Coupling_Parameters (cm_per_sec).
C CONSTS(78) is N_lcc in component EC_Coupling_Parameters (dimensionless).
C CONSTS(79) is I_Ca05 in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(80) is k_NaCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(81) is Km_Na in component EC_Coupling_Parameters (mM).
C CONSTS(82) is Km_Ca in component EC_Coupling_Parameters (mM).
C CONSTS(83) is k_sat in component EC_Coupling_Parameters (dimensionless).
C CONSTS(84) is eta in component EC_Coupling_Parameters (dimensionless).
C CONSTS(85) is I_bar_NaK in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(86) is Km_Nai in component EC_Coupling_Parameters (mM).
C CONSTS(87) is Km_Ko in component EC_Coupling_Parameters (mM).
C CONSTS(88) is I_bar_PCa in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(89) is Km_PCa in component EC_Coupling_Parameters (mM).
C CONSTS(90) is G_CaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(91) is G_NaB in component EC_Coupling_Parameters (uA_per_uF).
C CONSTS(92) is Pns in component EC_Coupling_Parameters (dimensionless).
C CONSTS(93) is Km_NS in component EC_Coupling_Parameters (mM).
C CONSTS(94) is I_up_bar in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(95) is Km_up0 in component EC_Coupling_Parameters (mM).
C CONSTS(96) is NSR_bar in component EC_Coupling_Parameters (mM).
C CONSTS(97) is tau_on in component EC_Coupling_Parameters (second).
C CONSTS(98) is tau_off in component EC_Coupling_Parameters (second).
C CONSTS(99) is G_max_rel in component EC_Coupling_Parameters (mM_per_sec).
C CONSTS(100) is d_Cai_th in component EC_Coupling_Parameters (mM).
C CONSTS(101) is Km_rel in component EC_Coupling_Parameters (mM).
C CONSTS(102) is CSQN_th in component EC_Coupling_Parameters (mM).
C CONSTS(103) is CSQN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(104) is Km_CSQN in component EC_Coupling_Parameters (mM).
C CONSTS(105) is tau_tr in component EC_Coupling_Parameters (second).
C CONSTS(106) is TRPN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(107) is CMDN_bar in component EC_Coupling_Parameters (mM).
C CONSTS(108) is INDO_bar in component EC_Coupling_Parameters (mM).
C CONSTS(109) is Km_TRPN in component EC_Coupling_Parameters (mM).
C CONSTS(110) is Km_CMDN in component EC_Coupling_Parameters (mM).
C CONSTS(111) is Km_INDO in component EC_Coupling_Parameters (mM).
C ALGBRC(63) is LR in component b1_AR_module (uM).
C ALGBRC(64) is LRG in component b1_AR_module (uM).
C ALGBRC(65) is RG in component b1_AR_module (uM).
C ALGBRC(81) 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(89) is PKA_DESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(24) is PKA_RESENS in component b1_AR_module (uM_per_sec).
C ALGBRC(82) is G_ACT in component b1_AR_module (uM_per_sec).
C ALGBRC(26) is HYD in component b1_AR_module (uM_per_sec).
C ALGBRC(28) is REASSOC in component b1_AR_module (per_sec).
C ALGBRC(66) is L in component b1_AR_module (uM).
C ALGBRC(67) is R in component b1_AR_module (uM).
C ALGBRC(68) 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(69) is PKAC_I in component PKA_module (uM).
C ALGBRC(70) is cAMP in component PKA_module (uM).
C ALGBRC(43) is Gsa_GTP in component cAMP_module (uM).
C ALGBRC(51) is Fsk in component cAMP_module (uM).
C ALGBRC(44) 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(45) is Gsa_GTP_AC in component cAMP_module (uM).
C ALGBRC(52) is Fsk_AC in component cAMP_module (uM).
C ALGBRC(49) is AC_ACT_GSA in component cAMP_module (uM).
C ALGBRC(47) is AC_ACT_BASAL in component cAMP_module (uM).
C ALGBRC(54) is AC_ACT_FSK in component cAMP_module (uM).
C ALGBRC(83) is PDE_ACT in component cAMP_module (uM).
C CONSTS(130) is PDE_IBMX in component cAMP_module (uM).
C ALGBRC(71) is PKI in component PKA_module (uM).
C ALGBRC(72) is A2RC_I in component PKA_module (uM).
C ALGBRC(73) is A2R_I in component PKA_module (uM).
C ALGBRC(74) is A2RC_II in component PKA_module (uM).
C ALGBRC(75) is A2R_II in component PKA_module (uM).
C ALGBRC(76) is ARC_I in component PKA_module (uM).
C ALGBRC(77) is ARC_II in component PKA_module (uM).
C ALGBRC(78) is PKA_temp in component PKA_module (uM).
C ALGBRC(79) is PKAC_II in component PKA_module (uM).
C ALGBRC(30) is PLB in component PLB_module (uM).
C ALGBRC(84) is PLB_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(61) is PLB_DEPHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(31) is Inhib1 in component PLB_module (uM).
C ALGBRC(56) is Inhib1p_PP1 in component PLB_module (uM).
C ALGBRC(85) is Inhib1_PHOSPH in component PLB_module (uM_per_sec).
C ALGBRC(33) 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(57) is Inhib1p in component PLB_module (uM).
C ALGBRC(58) is PP1 in component PLB_module (uM).
C ALGBRC(1) is frac_PLB_p in component PLB_module (dimensionless).
C ALGBRC(32) is frac_PLB in component PLB_module (dimensionless).
C CONSTS(120) is frac_PLB_o in component PLB_module (dimensionless).
C ALGBRC(35) is LCCa in component LCC_module (uM).
C ALGBRC(86) is LCCa_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(37) is LCCa_DEPHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(39) is LCCb in component LCC_module (uM).
C ALGBRC(87) is LCCb_PHOSPH in component LCC_module (uM_per_sec).
C ALGBRC(41) 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(121) is frac_LCCa_po in component LCC_module (dimensionless).
C ALGBRC(34) is frac_LCCb_p in component LCC_module (dimensionless).
C CONSTS(122) is frac_LCCb_po in component LCC_module (dimensionless).
C ALGBRC(36) is E_Na in component Nernst_Potentials (mV).
C ALGBRC(38) is E_K in component Nernst_Potentials (mV).
C ALGBRC(40) is E_Ca in component Nernst_Potentials (mV).
C CONSTS(123) is E_Cl in component Nernst_Potentials (mV).
C CONSTS(112) is R in component Nernst_Potentials (joules_per_mole_kelvin).
C CONSTS(113) is Frdy in component Nernst_Potentials (coulombs_per_mole).
C CONSTS(124) is FoRT in component Nernst_Potentials (per_mV).
C CONSTS(114) is z_Na in component Nernst_Potentials (dimensionless).
C CONSTS(115) is z_K in component Nernst_Potentials (dimensionless).
C CONSTS(116) 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(14) 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(15) is bh in component Fast_Na_Current (per_sec).
C ALGBRC(16) 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(42) 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(17) is b_lcc in component L_Type_Calcium_Current (per_sec).
C ALGBRC(19) 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(18) is tau_y_lcc in component L_Type_Calcium_Current (second).
C ALGBRC(25) is gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(27) is v_gamma in component L_Type_Calcium_Current (per_mM_per_sec).
C ALGBRC(29) 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(46) is i_bar_Ca in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(48) is i_bar_K in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(50) is f_avail in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(53) is I_Ca in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(55) is I_CaK in component L_Type_Calcium_Current (uA_per_uF).
C ALGBRC(59) 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(20) is tau_r_to in component Transient_Outward_K_Current (second).
C ALGBRC(21) is tau_s_to in component Transient_Outward_K_Current (second).
C ALGBRC(22) 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(60) 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(23) 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(125) 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(62) is I_ss in component Steady_State_K_Current (uA_per_uF).
C ALGBRC(80) is a_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(88) is b_Ki in component Time_Independent_K_Current (dimensionless).
C ALGBRC(90) is Ki_ss in component Time_Independent_K_Current (dimensionless).
C ALGBRC(91) is I_Ki in component Time_Independent_K_Current (uA_per_uF).
C ALGBRC(92) is Kp in component Plateau_K_Current (dimensionless).
C ALGBRC(93) is I_Kp in component Plateau_K_Current (uA_per_uF).
C ALGBRC(94) is s4 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(95) is s5 in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(96) is I_NCX in component Na_Ca_Exchanger_Current (uA_per_uF).
C CONSTS(126) is sigma in component Na_K_Pump_Current (dimensionless).
C ALGBRC(97) is f_NaK in component Na_K_Pump_Current (dimensionless).
C ALGBRC(98) is I_NaK in component Na_K_Pump_Current (uA_per_uF).
C ALGBRC(99) is I_PCa in component Sarcolemmal_Ca_Pump_Current (uA_per_uF).
C ALGBRC(100) is I_CaB in component Ca_Background_Current (uA_per_uF).
C ALGBRC(101) is I_NaB in component Na_Background_Current (uA_per_uF).
C ALGBRC(102) is I_Na_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(103) is I_K_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(104) is I_Ca_tot in component Total_Membrane_Currents (uA_per_uF).
C ALGBRC(105) 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(107) is ryr_off in component Calcium_Induced_Calcium_Release (dimensionless).
C ALGBRC(108) is g_rel in component Calcium_Induced_Calcium_Release (per_sec).
C ALGBRC(109) 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(110) is Km_up in component Other_SR_Fluxes_and_Concentrations (mM).
C ALGBRC(111) is I_up in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(112) is I_leak in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(113) is I_tr in component Other_SR_Fluxes_and_Concentrations (mM_per_sec).
C ALGBRC(115) 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(117) is SR_content in component Other_SR_Fluxes_and_Concentrations (uM).
C ALGBRC(114) is b_trpn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(116) is b_cmdn in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(118) is b_indo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C ALGBRC(119) is B_myo in component Cytoplasmic_Calcium_Buffering (dimensionless).
C CONSTS(127) is I_app in component Ion_Concentrations_and_Membrane_Potential (uA_per_uF).
C CONSTS(117) is V_hold in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(118) is V_test in component Ion_Concentrations_and_Membrane_Potential (mV).
C ALGBRC(13) is V_clamp in component Ion_Concentrations_and_Membrane_Potential (mV).
C CONSTS(119) is R_clamp in component Ion_Concentrations_and_Membrane_Potential (dimensionless).
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) = 1
CONSTS(2) = 0.0132
CONSTS(3) = 3.83
CONSTS(4) = 0.285
CONSTS(5) = 0.062
CONSTS(6) = 33
CONSTS(7) = 1.1e-3
CONSTS(8) = 2.2e-3
CONSTS(9) = 3.6e-3
CONSTS(10) = 2.2e-3
CONSTS(11) = 16
CONSTS(12) = 0.8
CONSTS(13) = 1.21e3
CONSTS(14) = 49.7e-3
CONSTS(15) = 5e3
CONSTS(16) = 38.9e-3
CONSTS(17) = 0
CONSTS(18) = 0
CONSTS(19) = 0.2
CONSTS(20) = 8.5
CONSTS(21) = 7.3
CONSTS(22) = 5
CONSTS(23) = 1.03e3
CONSTS(24) = 315
CONSTS(25) = 860
CONSTS(26) = 1.3
CONSTS(27) = 0.4
CONSTS(28) = 44
CONSTS(29) = 30
CONSTS(30) = 0.59
CONSTS(31) = 0.025
CONSTS(32) = 0.18
CONSTS(33) = 9.14
CONSTS(34) = 1.64
CONSTS(35) = 4.375
CONSTS(36) = 0.2e-3
CONSTS(37) = 106
CONSTS(38) = 0.89
CONSTS(39) = 0.3
CONSTS(40) = 54
CONSTS(41) = 21
CONSTS(42) = 8.5
CONSTS(43) = 7
CONSTS(44) = 60
CONSTS(45) = 1
CONSTS(46) = 14
CONSTS(47) = 1
CONSTS(48) = 1e-3
CONSTS(49) = 10
CONSTS(50) = 25e-3
CONSTS(51) = 25e-3
CONSTS(52) = 25e-3
CONSTS(53) = 54
CONSTS(54) = 21
CONSTS(55) = 8.52
CONSTS(56) = 3
CONSTS(57) = 10.1
CONSTS(58) = 3
CONSTS(59) = 20.8e-6
CONSTS(60) = 9.88e-7
CONSTS(61) = 9.3e-8
CONSTS(62) = 1.534e-4
CONSTS(63) = 310
CONSTS(64) = 140
CONSTS(65) = 5.4
CONSTS(66) = 1.8
CONSTS(67) = 8
CONSTS(68) = 0.35
CONSTS(69) = 0.07
CONSTS(70) = 0.24
CONSTS(71) = 0.008
CONSTS(72) = 300
CONSTS(73) = 2e3
CONSTS(74) = 5187.5
CONSTS(75) = 10
CONSTS(76) = 1.7469e-8
CONSTS(77) = 3.234e-11
CONSTS(78) = 3e5
CONSTS(79) = -0.458
CONSTS(80) = 1483
CONSTS(81) = 87.5
CONSTS(82) = 1.38
CONSTS(83) = 0.1
CONSTS(84) = 0.35
CONSTS(85) = 1.1
CONSTS(86) = 10
CONSTS(87) = 1.5
CONSTS(88) = 1.15
CONSTS(89) = 0.5e-3
CONSTS(90) = 2.8e-3
CONSTS(91) = 1.18e-3
CONSTS(92) = 0
CONSTS(93) = 1.2e-3
CONSTS(94) = 4.7
CONSTS(95) = 3e-4
CONSTS(96) = 15
CONSTS(97) = 2e-3
CONSTS(98) = 2e-3
CONSTS(99) = 60e3
CONSTS(100) = 0.18e-3
CONSTS(101) = 0.8e-3
CONSTS(102) = 8.75
CONSTS(103) = 15
CONSTS(104) = 0.8
CONSTS(105) = 5.7e-4
CONSTS(106) = 0.07
CONSTS(107) = 0.05
CONSTS(108) = 0.07
CONSTS(109) = 0.5128e-3
CONSTS(110) = 2.38e-3
CONSTS(111) = 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(112) = 8314
CONSTS(113) = 96485
CONSTS(114) = 1
CONSTS(115) = 1
CONSTS(116) = 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(117) = -40
CONSTS(118) = -10
CONSTS(119) = 0.02
CONSTS(120) = 0.961300
CONSTS(121) = 0.204100
CONSTS(122) = 0.233600
CONSTS(123) = - 40.0000
CONSTS(124) = (CONSTS(113)/CONSTS(112))/CONSTS(63)
CONSTS(125) = 2.10000
CONSTS(126) = (EXP(CONSTS(64)/67.3000) - 1.00000)/7.00000
CONSTS(127) = 0.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(125)
ALGBRC(3) = ( 0.320000*(STATES(18)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(18)+47.1300)))
ALGBRC(14) = 0.0800000*EXP(- STATES(18)/11.0000)
RATES(15) = 1000.00*( ALGBRC(3)*(1.00000 - STATES(15)) - ALGBRC(14)*STATES(15))
ALGBRC(4) = TERNRY(STATES(18).GE.- 40.0000, 0.00000, 0.135000*EXP((80.0000+STATES(18))/- 6.80000))
ALGBRC(15) = 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(15)*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(16) = 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(16)*STATES(17))
ALGBRC(6) = 400.000*EXP((STATES(18)+2.00000)/10.0000)
ALGBRC(17) = 50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
RATES(19) = ALGBRC(6)*(1.00000 - STATES(19)) - ALGBRC(17)*STATES(19)
RATES(20) = 2.00000*ALGBRC(6)*(1.00000 - STATES(20)) - (ALGBRC(17)/2.00000)*STATES(20)
ALGBRC(2) = STATES(10)/CONSTS(50)
ALGBRC(19) = CONSTS(72)*(( 0.375000*ALGBRC(2))/CONSTS(121)+0.625000)
RATES(21) = ALGBRC(19)*(1.00000 - STATES(21)) - CONSTS(73)*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(18) = 0.0200000+0.300000/(1.00000+EXP((STATES(18)+30.0000)/9.50000))
RATES(22) = (ALGBRC(7) - STATES(22))/ALGBRC(18)
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(18)+10.6000)/- 11.4200))
ALGBRC(20) = 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(20)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(18)+43.5000)/6.88410))
ALGBRC(21) = 0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
RATES(25) = (ALGBRC(9) - STATES(25))/ALGBRC(21)
ALGBRC(22) = 3.70000*EXP( - 1.00000*(STATES(18)+70.0000)/30.0000 ** 2.00000)+0.0350000
RATES(26) = (ALGBRC(9) - STATES(26))/ALGBRC(22)
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(18)+11.5000)/- 11.8200))
ALGBRC(23) = 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(23)
ALGBRC(26) = CONSTS(12)*STATES(4)
ALGBRC(28) = CONSTS(13)*STATES(5)*STATES(6)
RATES(5) = ALGBRC(26) - ALGBRC(28)
ALGBRC(25) = CONSTS(74)*STATES(14)
ALGBRC(27) = ALGBRC(25)*(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(19)/CONSTS(73)))
ALGBRC(29) = CONSTS(75)*(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(29)*(1.00000 - STATES(23)) - ALGBRC(27)*STATES(23)
CALL minimize(minfunc_1, CONSTS, VARIABLES, <VAR>)
ALGBRC(81) = CONSTS(7)*(ALGBRC(63)+ALGBRC(64))
ALGBRC(12) = CONSTS(8)*STATES(2)
RATES(2) = ALGBRC(81) - ALGBRC(12)
ALGBRC(82) = CONSTS(11)*(ALGBRC(65)+ALGBRC(64))
RATES(4) = ALGBRC(82) - ALGBRC(26)
RATES(6) = ALGBRC(82) - ALGBRC(28)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
ALGBRC(49) = ( CONSTS(20)*ALGBRC(45)*CONSTS(15))/(CONSTS(24)+CONSTS(15))
ALGBRC(47) = ( CONSTS(19)*ALGBRC(44)*CONSTS(15))/(CONSTS(23)+CONSTS(15))
CALL minimize(minfunc_3, CONSTS, VARIABLES, <VAR>)
ALGBRC(54) = ( CONSTS(21)*ALGBRC(52)*CONSTS(15))/(CONSTS(25)+CONSTS(15))
ALGBRC(83) = ( CONSTS(22)*CONSTS(128)*ALGBRC(70))/(CONSTS(26)+ALGBRC(70))
RATES(7) = (ALGBRC(47)+ALGBRC(49)+ALGBRC(54)) - ALGBRC(83)
ALGBRC(30) = CONSTS(37) - STATES(8)
ALGBRC(84) = ( CONSTS(40)*ALGBRC(69)*ALGBRC(30))/(CONSTS(41)+ALGBRC(30))
CALL minimize(minfunc_4, CONSTS, VARIABLES, <VAR>)
ALGBRC(61) = ( CONSTS(42)*ALGBRC(58)*STATES(8))/(CONSTS(43)+STATES(8))
RATES(8) = ALGBRC(84) - ALGBRC(61)
ALGBRC(31) = CONSTS(39) - STATES(9)
ALGBRC(85) = ( CONSTS(44)*ALGBRC(69)*ALGBRC(31))/(CONSTS(45)+ALGBRC(31))
ALGBRC(33) = ( CONSTS(46)*STATES(9))/(CONSTS(47)+STATES(9))
RATES(9) = ALGBRC(85) - ALGBRC(33)
ALGBRC(35) = CONSTS(50) - STATES(10)
ALGBRC(86) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(35))/(CONSTS(54)+ CONSTS(49)*ALGBRC(35))
ALGBRC(37) = ( CONSTS(49)*CONSTS(57)*CONSTS(52)*STATES(10))/(CONSTS(58)+ CONSTS(49)*STATES(10))
RATES(10) = ALGBRC(86) - ALGBRC(37)
ALGBRC(39) = CONSTS(50) - STATES(11)
ALGBRC(87) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(39))/(CONSTS(54)+ CONSTS(49)*ALGBRC(39))
ALGBRC(41) = ( CONSTS(49)*CONSTS(55)*CONSTS(51)*STATES(11))/(CONSTS(56)+ CONSTS(49)*STATES(11))
RATES(11) = ALGBRC(87) - ALGBRC(41)
ALGBRC(89) = CONSTS(9)*ALGBRC(69)*STATES(1)
ALGBRC(24) = CONSTS(10)*STATES(3)
RATES(1) = (ALGBRC(12) - ALGBRC(81))+(ALGBRC(24) - ALGBRC(89))
RATES(3) = ALGBRC(89) - ALGBRC(24)
ALGBRC(36) = (1.00000/( CONSTS(124)*CONSTS(114)))*log(CONSTS(64)/STATES(12))
ALGBRC(42) = CONSTS(67)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(36))
ALGBRC(94) = EXP( CONSTS(84)*STATES(18)*CONSTS(124))*STATES(12) ** 3.00000*CONSTS(66)
ALGBRC(95) = EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124))*CONSTS(64) ** 3.00000*STATES(14)
ALGBRC(96) = (CONSTS(80)/( (CONSTS(81) ** 3.00000+CONSTS(64) ** 3.00000)*(CONSTS(82)+CONSTS(66))*(1.00000+ CONSTS(83)*EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124)))))*(ALGBRC(94) - ALGBRC(95))
ALGBRC(97) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(124))+ 0.0365000*CONSTS(126)*EXP( - STATES(18)*CONSTS(124)))
ALGBRC(98) = ( (( CONSTS(85)*ALGBRC(97))/(1.00000+CONSTS(86)/STATES(12) ** 1.50000))*CONSTS(65))/(CONSTS(65)+CONSTS(87))
ALGBRC(101) = CONSTS(91)*(STATES(18) - ALGBRC(36))
ALGBRC(102) = ALGBRC(42)+ALGBRC(101)+ 3.00000*ALGBRC(96)+ 3.00000*ALGBRC(98)
RATES(12) = ( - 1000.00*ALGBRC(102)*CONSTS(62))/( CONSTS(59)*CONSTS(114)*CONSTS(113))
ALGBRC(48) = ( CONSTS(77)*STATES(18)*CONSTS(113)*CONSTS(124)*( STATES(13)*EXP( STATES(18)*CONSTS(124)) - CONSTS(65)))/(EXP( STATES(18)*CONSTS(124)) - 1.00000)
ALGBRC(34) = STATES(11)/CONSTS(50)
ALGBRC(50) = 0.500000*(( 0.400000*ALGBRC(34))/CONSTS(122)+0.600000)
ALGBRC(46) = ( CONSTS(76)*4.00000*STATES(18)*CONSTS(113)*CONSTS(124)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(124)) - 0.341000*CONSTS(66)))/(EXP( 2.00000*STATES(18)*CONSTS(124)) - 1.00000)
ALGBRC(53) = ALGBRC(46)*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(55) = (ALGBRC(48)/(1.00000+ALGBRC(53)/CONSTS(79)))*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(38) = (1.00000/( CONSTS(124)*CONSTS(115)))*log(CONSTS(65)/STATES(13))
ALGBRC(60) = CONSTS(68)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(38))
ALGBRC(62) = CONSTS(69)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(38))
ALGBRC(80) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(38)) - 59.2150)))
ALGBRC(88) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(38)))+EXP( 0.0617500*((STATES(18) - ALGBRC(38)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(38))+4.75300)))
ALGBRC(90) = ALGBRC(80)/(ALGBRC(80)+ALGBRC(88))
ALGBRC(91) = CONSTS(70)* (CONSTS(65)/5.40000) ** (1.0 / 2)*ALGBRC(90)*(STATES(18) - ALGBRC(38))
ALGBRC(92) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
ALGBRC(93) = CONSTS(71)*ALGBRC(92)*(STATES(18) - ALGBRC(38))
ALGBRC(103) = ((ALGBRC(60)+ALGBRC(62)+ALGBRC(91)+ALGBRC(93)) - 2.00000*ALGBRC(98))+ALGBRC(55)
RATES(13) = ( - 1000.00*ALGBRC(103)*CONSTS(62))/( CONSTS(59)*CONSTS(115)*CONSTS(113))
ALGBRC(99) = ( CONSTS(88)*STATES(14))/(CONSTS(89)+STATES(14))
ALGBRC(40) = (1.00000/( CONSTS(124)*CONSTS(116)))*log(CONSTS(66)/STATES(14))
ALGBRC(100) = CONSTS(90)*(STATES(18) - ALGBRC(40))
ALGBRC(104) = (ALGBRC(53)+ALGBRC(100)+ALGBRC(99)) - 2.00000*ALGBRC(96)
RATES(18) = - 1000.00*((ALGBRC(104)+ALGBRC(103)+ALGBRC(102)) - CONSTS(127))
RATES(30) = TERNRY( - 1000.00*((ALGBRC(104)+ALGBRC(103)+ALGBRC(102)) - CONSTS(127)).GT.30000.0, 1.00000 - 10000.0*STATES(30), 1.00000)
ALGBRC(32) = ALGBRC(30)/CONSTS(37)
ALGBRC(110) = ( CONSTS(95)*(1.00000+ 2.00000*ALGBRC(32)))/(1.00000+ 2.00000*CONSTS(120))
ALGBRC(111) = ( CONSTS(94)*STATES(14) ** 2.00000)/(ALGBRC(110) ** 2.00000+STATES(14) ** 2.00000)
ALGBRC(112) = ( CONSTS(94)*STATES(31))/CONSTS(96)
ALGBRC(113) = (STATES(31) - STATES(29))/CONSTS(105)
RATES(31) = (ALGBRC(111) - ALGBRC(112)) - ( ALGBRC(113)*CONSTS(61))/CONSTS(60)
ALGBRC(105) = STATES(30)+0.00200000
ALGBRC(106) = 1.00000 - EXP(- ALGBRC(105)/CONSTS(97))
ALGBRC(107) = EXP(- ALGBRC(105)/CONSTS(98))
ALGBRC(108) = CONSTS(99)/(1.00000+EXP((ALGBRC(104)+5.00000)/0.900000))
ALGBRC(109) = ALGBRC(108)*ALGBRC(106)*ALGBRC(107)*(STATES(29) - STATES(14))
ALGBRC(115) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(104))/CONSTS(104)+STATES(29) ** 2.00000)
RATES(29) = ALGBRC(115)*(ALGBRC(113) - ALGBRC(109))
ALGBRC(114) = ( CONSTS(106)*CONSTS(109))/CONSTS(109)+STATES(14) ** 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(119) = 1.00000/(1.00000+ALGBRC(116)+ALGBRC(114)+ALGBRC(114)+ALGBRC(118))
RATES(14) = - ALGBRC(119)*((( 1000.00*ALGBRC(104)*CONSTS(62))/( CONSTS(59)*CONSTS(116)*CONSTS(113))+( (ALGBRC(111) - ALGBRC(112))*CONSTS(60))/CONSTS(59)) - ( ALGBRC(109)*CONSTS(61))/CONSTS(59))
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(14) = 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(15) = 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(16) = 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(17) = 50.0000*EXP(( - 1.00000*(STATES(18)+2.00000))/13.0000)
ALGBRC(2) = STATES(10)/CONSTS(50)
ALGBRC(19) = CONSTS(72)*(( 0.375000*ALGBRC(2))/CONSTS(121)+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(18) = 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(20) = 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(21) = 0.350000*EXP( - 1.00000*(STATES(18)+70.0000)/15.0000 ** 2.00000)+0.0350000
ALGBRC(22) = 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(23) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(18)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(18)+38.0000)))
ALGBRC(26) = CONSTS(12)*STATES(4)
ALGBRC(28) = CONSTS(13)*STATES(5)*STATES(6)
ALGBRC(25) = CONSTS(74)*STATES(14)
ALGBRC(27) = ALGBRC(25)*(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(19)/CONSTS(73)))
ALGBRC(29) = CONSTS(75)*(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(81) = CONSTS(7)*(ALGBRC(63)+ALGBRC(64))
ALGBRC(12) = CONSTS(8)*STATES(2)
ALGBRC(82) = CONSTS(11)*(ALGBRC(65)+ALGBRC(64))
ALGBRC(49) = ( CONSTS(20)*ALGBRC(45)*CONSTS(15))/(CONSTS(24)+CONSTS(15))
ALGBRC(47) = ( CONSTS(19)*ALGBRC(44)*CONSTS(15))/(CONSTS(23)+CONSTS(15))
ALGBRC(54) = ( CONSTS(21)*ALGBRC(52)*CONSTS(15))/(CONSTS(25)+CONSTS(15))
ALGBRC(83) = ( CONSTS(22)*CONSTS(128)*ALGBRC(70))/(CONSTS(26)+ALGBRC(70))
ALGBRC(30) = CONSTS(37) - STATES(8)
ALGBRC(84) = ( CONSTS(40)*ALGBRC(69)*ALGBRC(30))/(CONSTS(41)+ALGBRC(30))
ALGBRC(61) = ( CONSTS(42)*ALGBRC(58)*STATES(8))/(CONSTS(43)+STATES(8))
ALGBRC(31) = CONSTS(39) - STATES(9)
ALGBRC(85) = ( CONSTS(44)*ALGBRC(69)*ALGBRC(31))/(CONSTS(45)+ALGBRC(31))
ALGBRC(33) = ( CONSTS(46)*STATES(9))/(CONSTS(47)+STATES(9))
ALGBRC(35) = CONSTS(50) - STATES(10)
ALGBRC(86) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(35))/(CONSTS(54)+ CONSTS(49)*ALGBRC(35))
ALGBRC(37) = ( CONSTS(49)*CONSTS(57)*CONSTS(52)*STATES(10))/(CONSTS(58)+ CONSTS(49)*STATES(10))
ALGBRC(39) = CONSTS(50) - STATES(11)
ALGBRC(87) = ( CONSTS(49)*CONSTS(53)*ALGBRC(79)*ALGBRC(39))/(CONSTS(54)+ CONSTS(49)*ALGBRC(39))
ALGBRC(41) = ( CONSTS(49)*CONSTS(55)*CONSTS(51)*STATES(11))/(CONSTS(56)+ CONSTS(49)*STATES(11))
ALGBRC(89) = CONSTS(9)*ALGBRC(69)*STATES(1)
ALGBRC(24) = CONSTS(10)*STATES(3)
ALGBRC(36) = (1.00000/( CONSTS(124)*CONSTS(114)))*log(CONSTS(64)/STATES(12))
ALGBRC(42) = CONSTS(67)*STATES(15) ** 3.00000*STATES(16)*STATES(17)*(STATES(18) - ALGBRC(36))
ALGBRC(94) = EXP( CONSTS(84)*STATES(18)*CONSTS(124))*STATES(12) ** 3.00000*CONSTS(66)
ALGBRC(95) = EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124))*CONSTS(64) ** 3.00000*STATES(14)
ALGBRC(96) = (CONSTS(80)/( (CONSTS(81) ** 3.00000+CONSTS(64) ** 3.00000)*(CONSTS(82)+CONSTS(66))*(1.00000+ CONSTS(83)*EXP( (CONSTS(84) - 1.00000)*STATES(18)*CONSTS(124)))))*(ALGBRC(94) - ALGBRC(95))
ALGBRC(97) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(18)*CONSTS(124))+ 0.0365000*CONSTS(126)*EXP( - STATES(18)*CONSTS(124)))
ALGBRC(98) = ( (( CONSTS(85)*ALGBRC(97))/(1.00000+CONSTS(86)/STATES(12) ** 1.50000))*CONSTS(65))/(CONSTS(65)+CONSTS(87))
ALGBRC(101) = CONSTS(91)*(STATES(18) - ALGBRC(36))
ALGBRC(102) = ALGBRC(42)+ALGBRC(101)+ 3.00000*ALGBRC(96)+ 3.00000*ALGBRC(98)
ALGBRC(48) = ( CONSTS(77)*STATES(18)*CONSTS(113)*CONSTS(124)*( STATES(13)*EXP( STATES(18)*CONSTS(124)) - CONSTS(65)))/(EXP( STATES(18)*CONSTS(124)) - 1.00000)
ALGBRC(34) = STATES(11)/CONSTS(50)
ALGBRC(50) = 0.500000*(( 0.400000*ALGBRC(34))/CONSTS(122)+0.600000)
ALGBRC(46) = ( CONSTS(76)*4.00000*STATES(18)*CONSTS(113)*CONSTS(124)*( 0.00100000*EXP( 2.00000*STATES(18)*CONSTS(124)) - 0.341000*CONSTS(66)))/(EXP( 2.00000*STATES(18)*CONSTS(124)) - 1.00000)
ALGBRC(53) = ALGBRC(46)*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(55) = (ALGBRC(48)/(1.00000+ALGBRC(53)/CONSTS(79)))*CONSTS(78)*ALGBRC(50)*STATES(19) ** 4.00000*STATES(21)*STATES(22)*STATES(23)
ALGBRC(38) = (1.00000/( CONSTS(124)*CONSTS(115)))*log(CONSTS(65)/STATES(13))
ALGBRC(60) = CONSTS(68)*STATES(24)*( 0.886000*STATES(25)+ 0.114000*STATES(26))*(STATES(18) - ALGBRC(38))
ALGBRC(62) = CONSTS(69)*STATES(27)*STATES(28)*(STATES(18) - ALGBRC(38))
ALGBRC(80) = 1.02000/(1.00000+EXP( 0.238500*((STATES(18) - ALGBRC(38)) - 59.2150)))
ALGBRC(88) = ( 0.491240*EXP( 0.0803200*((STATES(18)+5.47600) - ALGBRC(38)))+EXP( 0.0617500*((STATES(18) - ALGBRC(38)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(18) - ALGBRC(38))+4.75300)))
ALGBRC(90) = ALGBRC(80)/(ALGBRC(80)+ALGBRC(88))
ALGBRC(91) = CONSTS(70)* (CONSTS(65)/5.40000) ** (1.0 / 2)*ALGBRC(90)*(STATES(18) - ALGBRC(38))
ALGBRC(92) = 1.00000/(1.00000+EXP((7.48800 - STATES(18))/5.98000))
ALGBRC(93) = CONSTS(71)*ALGBRC(92)*(STATES(18) - ALGBRC(38))
ALGBRC(103) = ((ALGBRC(60)+ALGBRC(62)+ALGBRC(91)+ALGBRC(93)) - 2.00000*ALGBRC(98))+ALGBRC(55)
ALGBRC(99) = ( CONSTS(88)*STATES(14))/(CONSTS(89)+STATES(14))
ALGBRC(40) = (1.00000/( CONSTS(124)*CONSTS(116)))*log(CONSTS(66)/STATES(14))
ALGBRC(100) = CONSTS(90)*(STATES(18) - ALGBRC(40))
ALGBRC(104) = (ALGBRC(53)+ALGBRC(100)+ALGBRC(99)) - 2.00000*ALGBRC(96)
ALGBRC(32) = ALGBRC(30)/CONSTS(37)
ALGBRC(110) = ( CONSTS(95)*(1.00000+ 2.00000*ALGBRC(32)))/(1.00000+ 2.00000*CONSTS(120))
ALGBRC(111) = ( CONSTS(94)*STATES(14) ** 2.00000)/(ALGBRC(110) ** 2.00000+STATES(14) ** 2.00000)
ALGBRC(112) = ( CONSTS(94)*STATES(31))/CONSTS(96)
ALGBRC(113) = (STATES(31) - STATES(29))/CONSTS(105)
ALGBRC(105) = STATES(30)+0.00200000
ALGBRC(106) = 1.00000 - EXP(- ALGBRC(105)/CONSTS(97))
ALGBRC(107) = EXP(- ALGBRC(105)/CONSTS(98))
ALGBRC(108) = CONSTS(99)/(1.00000+EXP((ALGBRC(104)+5.00000)/0.900000))
ALGBRC(109) = ALGBRC(108)*ALGBRC(106)*ALGBRC(107)*(STATES(29) - STATES(14))
ALGBRC(115) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(104))/CONSTS(104)+STATES(29) ** 2.00000)
ALGBRC(114) = ( CONSTS(106)*CONSTS(109))/CONSTS(109)+STATES(14) ** 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(119) = 1.00000/(1.00000+ALGBRC(116)+ALGBRC(114)+ALGBRC(114)+ALGBRC(118))
ALGBRC(1) = STATES(8)/CONSTS(37)
ALGBRC(13) = TERNRY(VOI.GT.59.1000.AND.VOI.LT.59.5000, CONSTS(118), CONSTS(117))
ALGBRC(59) = ALGBRC(53)+ALGBRC(55)
ALGBRC(117) = 1000.00*(( (STATES(29)+STATES(29)/ALGBRC(115))*CONSTS(61))/CONSTS(59)+( STATES(31)*CONSTS(60))/CONSTS(59))
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0[1] = abs(CONSTS(130) - ( CONSTS(128)*CONSTS(129))/CONSTS(29))
minfunc_0[2] = abs(CONSTS(128) - (CONSTS(16) - CONSTS(130)))
minfunc_0[3] = abs(CONSTS(129) - (CONSTS(17) - CONSTS(130)))
RETURN
END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_1[1] = abs(ALGBRC(63) - ( ALGBRC(66)*ALGBRC(67))/CONSTS(4))
minfunc_1[2] = abs(ALGBRC(64) - ( ALGBRC(63)*ALGBRC(68))/CONSTS(5))
minfunc_1[3] = abs(ALGBRC(65) - ( ALGBRC(67)*ALGBRC(68))/CONSTS(6))
minfunc_1[4] = abs(ALGBRC(66) - ((CONSTS(1) - ALGBRC(63)) - ALGBRC(64)))
minfunc_1[5] = abs(ALGBRC(67) - (((STATES(1) - ALGBRC(63)) - ALGBRC(64)) - ALGBRC(65)))
minfunc_1[6] = abs(ALGBRC(68) - ((CONSTS(3) - ALGBRC(64)) - ALGBRC(65)))
minfunc_1[7] = abs(ALGBRC(71) - ( CONSTS(32)*CONSTS(36))/(CONSTS(36)+ALGBRC(69)+ALGBRC(79)))
minfunc_1[8] = abs(ALGBRC(72) - (ALGBRC(69)/CONSTS(35))*ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36)))
minfunc_1[9] = abs(ALGBRC(73) - ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36)))
minfunc_1[10] = abs(ALGBRC(74) - (ALGBRC(79)/CONSTS(35))*ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36)))
minfunc_1[11] = abs(ALGBRC(75) - ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36)))
minfunc_1[12] = abs(ALGBRC(76) - (CONSTS(33)/ALGBRC(70))*ALGBRC(72))
minfunc_1[13] = abs(ALGBRC(77) - (CONSTS(33)/ALGBRC(70))*ALGBRC(74))
minfunc_1[14] = abs(ALGBRC(78) - (( CONSTS(33)*CONSTS(34))/CONSTS(35)+( CONSTS(33)*ALGBRC(70))/CONSTS(35)+ALGBRC(70) ** 2.00000/CONSTS(35)))
minfunc_1[15] = abs(ALGBRC(70) - ((STATES(7) - (ALGBRC(76)+ 2.00000*ALGBRC(72)+ 2.00000*ALGBRC(73))) - (ALGBRC(77)+ 2.00000*ALGBRC(74)+ 2.00000*ALGBRC(75))))
minfunc_1[16] = abs(0.00000 - ( 2.00000*CONSTS(30)*ALGBRC(70) ** 2.00000 - ALGBRC(69)*(1.00000+ALGBRC(71)/CONSTS(36))*( ALGBRC(78)*ALGBRC(69)+ALGBRC(70) ** 2.00000)))
minfunc_1[17] = abs(0.00000 - ( 2.00000*CONSTS(31)*ALGBRC(70) ** 2.00000 - ALGBRC(79)*(1.00000+ALGBRC(71)/CONSTS(36))*( ALGBRC(78)*ALGBRC(79)+ALGBRC(70) ** 2.00000)))
RETURN
END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_2[1] = abs(ALGBRC(45) - ( ALGBRC(43)*ALGBRC(44))/CONSTS(27))
minfunc_2[2] = abs(ALGBRC(43) - (STATES(4) - ALGBRC(45)))
minfunc_2[3] = abs(ALGBRC(44) - (CONSTS(14) - ALGBRC(45)))
RETURN
END
REAL FUNCTION minfunc_3(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_3[1] = abs(ALGBRC(52) - ( ALGBRC(51)*ALGBRC(44))/CONSTS(28))
minfunc_3[2] = abs(ALGBRC(51) - (CONSTS(18) - ALGBRC(52)))
RETURN
END
REAL FUNCTION minfunc_4(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_4[1] = abs(ALGBRC(56) - ( ALGBRC(57)*ALGBRC(58))/CONSTS(48))
minfunc_4[2] = abs(ALGBRC(57) - (STATES(9) - ALGBRC(56)))
minfunc_4[3] = abs(ALGBRC(58) - (CONSTS(38) - ALGBRC(56)))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END