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 145 entries in the algebraic variable array.
C There are a total of 73 entries in each of the rate and state variable arrays.
C There are a total of 188 entries in the constant variable array.
C
C
C VOI is time in component Environment (second).
C CONSTS(1) is F in component Parameters (sec_A_per_mol).
C CONSTS(2) is vol_cyt in component Parameters (cubic_Meter).
C CONSTS(3) is vol_md in component Parameters (cubic_Meter).
C CONSTS(4) is vol_SR in component Parameters (cubic_Meter).
C CONSTS(5) is vol_cyt_Ca in component Parameters (cubic_Meter).
C CONSTS(6) is K_mr in component Parameters (Molar).
C CONSTS(7) is n_H in component Parameters (dimensionless).
C CONSTS(8) is v_RyR in component Parameters (per_sec).
C CONSTS(9) is K_a in component Parameters (Molar).
C CONSTS(10) is K_b in component Parameters (Molar).
C CONSTS(11) is K_c in component Parameters (Molar).
C CONSTS(12) is K_c_rate in component Parameters (per_sec).
C CONSTS(13) is v_IP3R in component Parameters (per_sec).
C CONSTS(14) is v_IP3 in component Parameters (per_sec).
C CONSTS(15) is conc_IP3_ref in component Parameters (Molar).
C CONSTS(16) is alpha_4 in component Parameters (dimensionless).
C CONSTS(17) is k4 in component Parameters (Molar).
C CONSTS(18) is I_r in component Parameters (per_sec).
C CONSTS(19) is d1 in component Parameters (Molar).
C CONSTS(20) is a4 in component Parameters (per_M_sec).
C CONSTS(21) is d4 in component Parameters (Molar).
C CONSTS(22) is a5 in component Parameters (per_M_sec).
C CONSTS(23) is d5 in component Parameters (Molar).
C CONSTS(24) is G_CaSOC_max in component Parameters (siemens2).
C CONSTS(25) is K_SOC in component Parameters (Molar).
C CONSTS(26) is delta_CT in component Parameters (volt).
C CONSTS(27) is J_NaCl_max in component Parameters (mol_per_sec).
C CONSTS(28) is J_KCl_max in component Parameters (mol_per_sec).
C CONSTS(29) is J_NKCC_max in component Parameters (mol_per_sec).
C CONSTS(30) is K_NKCC_Na in component Parameters (Molar).
C CONSTS(31) is K_NKCC_K in component Parameters (Molar).
C CONSTS(32) is K_NKCC_Cl1 in component Parameters (Molar).
C CONSTS(33) is K_NKCC_Cl2 in component Parameters (Molar).
C CONSTS(34) is k_on in component Parameters (per_M_sec).
C CONSTS(35) is k_off in component Parameters (per_sec).
C CONSTS(36) is alpha1 in component Parameters (per_sec).
C CONSTS(37) is alpha2 in component Parameters (per_sec).
C CONSTS(38) is alpha3 in component Parameters (per_sec).
C CONSTS(39) is G_ClCa_max in component Parameters (siemens2).
C CONSTS(40) is G_K_b in component Parameters (siemens2).
C CONSTS(41) is G_Na_b in component Parameters (siemens2).
C CONSTS(42) is G_Ca_b in component Parameters (siemens2).
C CONSTS(43) is G_Cl_b in component Parameters (siemens2).
C CONSTS(44) is G_Kir_max in component Parameters (siemens2).
C CONSTS(45) is conc_K_ref in component Parameters (Molar).
C CONSTS(46) is G_Kv_max in component Parameters (siemens2).
C CONSTS(47) is n_ATP in component Parameters (dimensionless).
C CONSTS(48) is G_KATP_max in component Parameters (siemens2).
C CONSTS(49) is G_KCa_max in component Parameters (siemens2).
C CONSTS(50) is Tau_PF in component Parameters (second).
C CONSTS(51) is Tau_PS in component Parameters (second).
C CONSTS(52) is K_m_K in component Parameters (Molar).
C CONSTS(53) is K_m_Na_alpha1 in component Parameters (Molar).
C CONSTS(54) is K_m_Na_alpha2 in component Parameters (Molar).
C CONSTS(55) is I_NaCa_max in component Parameters (ampere).
C CONSTS(56) is K_m_Cai in component Parameters (Molar).
C CONSTS(57) is K_m_Cao in component Parameters (Molar).
C CONSTS(58) is K_m_Nai in component Parameters (Molar).
C CONSTS(59) is K_m_Nao in component Parameters (Molar).
C CONSTS(60) is k_sat in component Parameters (dimensionless).
C CONSTS(61) is gamma in component Parameters (dimensionless).
C CONSTS(62) is K_mCa_act in component Parameters (Molar).
C CONSTS(63) is G_VONa_max in component Parameters (siemens2).
C CONSTS(64) is Tau_m in component Parameters (second).
C CONSTS(65) is Tau_h in component Parameters (second).
C CONSTS(66) is K_m_CaP in component Parameters (Molar).
C CONSTS(67) is G_CaL_max in component Parameters (siemens2).
C CONSTS(68) is conc_Bf_tot in component Parameters (Molar).
C CONSTS(69) is k_Bf_on in component Parameters (per_M_sec).
C CONSTS(70) is k_Bf_off in component Parameters (per_sec).
C CONSTS(71) is conc_Calseq_SR_tot in component Parameters (Molar).
C CONSTS(72) is k_Calseq_on in component Parameters (per_M_sec).
C CONSTS(73) is k_Calseq_off in component Parameters (per_sec).
C CONSTS(74) is f_cyt in component Parameters (dimensionless).
C CONSTS(75) is f_md in component Parameters (dimensionless).
C CONSTS(76) is R in component Parameters (J_per_mol_K).
C CONSTS(77) is T in component Parameters (kelvin).
C CONSTS(78) is Area in component Parameters (sq_meter).
C CONSTS(79) is L in component Parameters (meter).
C CONSTS(80) is hcon in component Parameters (dimensionless).
C CONSTS(81) is D_Ca in component Parameters (sq_m_per_sec).
C CONSTS(177) is D_K in component Parameters (sq_m_per_sec).
C CONSTS(178) is D_Na in component Parameters (sq_m_per_sec).
C CONSTS(82) is D_Cl in component Parameters (sq_m_per_sec).
C STATES(1) is conc_K_cyt in component Potassium_Conc_cyt (Molar).
C ALGBRC(38) is I_K_cyt_b in component Background_Currents (ampere).
C ALGBRC(88) is I_K_cyt_ir in component K_Inward_Rectifier_Current (ampere).
C ALGBRC(96) is I_K_cyt_ATP in component K_ATP_Current (ampere).
C ALGBRC(93) is I_K_cyt_v in component K_Delayed_Rectifier_Current (ampere).
C ALGBRC(99) is I_K_cyt_Ca in component K_Ca_Current (ampere).
C ALGBRC(141) is I_NaK_cyt_alpha1 in component NaK_alpha_Current (ampere).
C ALGBRC(78) is J_cyt_KCl in component Cl_flux_KCl (mol_per_sec).
C ALGBRC(80) is J_cyt_NKCC in component Cl_flux_NKCC (mol_per_sec).
C ALGBRC(116) is J_Kdiff in component Md_Cyt_Flux_K (mol_per_sec).
C STATES(2) is conc_K_md in component Potassium_Conc_md (Molar).
C ALGBRC(42) is I_K_md_b in component Background_Currents (ampere).
C ALGBRC(91) is I_K_md_ir in component K_Inward_Rectifier_Current (ampere).
C ALGBRC(97) is I_K_md_ATP in component K_ATP_Current (ampere).
C ALGBRC(95) is I_K_md_v in component K_Delayed_Rectifier_Current (ampere).
C ALGBRC(101) is I_K_md_Ca in component K_Ca_Current (ampere).
C ALGBRC(142) is I_NaK_md_alpha2 in component NaK_alpha_Current (ampere).
C ALGBRC(79) is J_md_KCl in component Cl_flux_KCl (mol_per_sec).
C ALGBRC(81) is J_md_NKCC in component Cl_flux_NKCC (mol_per_sec).
C STATES(3) is conc_Na_cyt in component Sodium_Conc_cyt (Molar).
C ALGBRC(44) is I_Na_cyt_b in component Background_Currents (ampere).
C ALGBRC(109) is I_VONa_cyt in component VONa_Current (ampere).
C ALGBRC(74) is I_Na_cyt_SOC in component SOC_Current_Na (ampere).
C ALGBRC(76) is J_cyt_NaCl in component Cl_flux_NaCl (mol_per_sec).
C ALGBRC(118) is J_Nadiff in component Md_Cyt_Flux_Na (mol_per_sec).
C STATES(4) is conc_Na_md in component Sodium_Conc_md (Molar).
C ALGBRC(46) is I_Na_md_b in component Background_Currents (ampere).
C ALGBRC(110) is I_VONa_md in component VONa_Current (ampere).
C ALGBRC(75) is I_Na_md_SOC in component SOC_Current_Na (ampere).
C ALGBRC(77) is J_md_NaCl in component Cl_flux_NaCl (mol_per_sec).
C ALGBRC(108) is I_NaCa_md in component NCX_Current (ampere).
C STATES(5) is conc_Cl_cyt in component Chloride_Conc_cyt (Molar).
C ALGBRC(48) is I_Cl_cyt_b in component Background_Currents (ampere).
C ALGBRC(82) is I_Cl_cyt_Ca in component ClCa_Current (ampere).
C ALGBRC(122) is J_Cldiff in component Md_Cyt_Flux_Cl (mol_per_sec).
C STATES(6) is conc_Cl_md in component Chloride_Conc_md (Molar).
C ALGBRC(50) is I_Cl_md_b in component Background_Currents (ampere).
C ALGBRC(83) is I_Cl_md_Ca in component ClCa_Current (ampere).
C STATES(7) is conc_Ca_cyt in component Calcium_Conc_cyt (Molar).
C ALGBRC(52) is I_Ca_cyt_b in component Background_Currents (ampere).
C ALGBRC(127) is I_Ca_cyt_P in component Ca_Pump_Current (ampere).
C ALGBRC(72) is I_Ca_cyt_SOC in component SOC_Current_Ca (ampere).
C ALGBRC(112) is I_Ca_cyt_L in component CaL_Current (ampere).
C ALGBRC(132) is I_cyt_SERCa_IP3R in component SERCA_Currents_IP3R (ampere).
C ALGBRC(56) is I_cyt_RyR in component RyR_Currents (ampere).
C ALGBRC(120) is J_Cadiff in component Md_Cyt_Flux_Ca (mol_per_sec).
C ALGBRC(64) is I_cyt_IP3R in component IP3R_Currents (ampere).
C STATES(8) is conc_CaMN__ in component MLCK_Activation_CaM (Molar).
C STATES(9) is conc_BfCa_cyt in component Bf_Kinetics (Molar).
C ALGBRC(134) is I_cyt_SERCa_RyR in component SERCA_Currents_RyR (ampere).
C STATES(10) is conc_CaM_C_ in component MLCK_Activation_CaM (Molar).
C STATES(11) is conc_CaM_CM in component MLCK_Activation_CaM (Molar).
C STATES(12) is conc_CaMNCM in component MLCK_Activation_CaM (Molar).
C STATES(13) is conc_CaM___ in component MLCK_Activation_CaM (Molar).
C STATES(14) is conc_BfCa_md in component Bf_Kinetics (Molar).
C STATES(15) is conc_Ca_md in component Calcium_Conc_md (Molar).
C ALGBRC(54) is I_Ca_md_b in component Background_Currents (ampere).
C ALGBRC(128) is I_Ca_md_P in component Ca_Pump_Current (ampere).
C ALGBRC(73) is I_Ca_md_SOC in component SOC_Current_Ca (ampere).
C ALGBRC(114) is I_Ca_md_L in component CaL_Current (ampere).
C ALGBRC(135) is I_md_SERCa_RyR in component SERCA_Currents_RyR (ampere).
C ALGBRC(58) is I_md_RyR in component RyR_Currents (ampere).
C ALGBRC(68) is I_md_IP3R in component IP3R_Currents (ampere).
C ALGBRC(133) is I_md_SERCa_IP3R in component SERCA_Currents_IP3R (ampere).
C STATES(16) is conc_CaMN__ in component MLCK_Activation_CaM_md (Molar).
C STATES(17) is conc_CaM_C_ in component MLCK_Activation_CaM_md (Molar).
C STATES(18) is conc_CaM_CM in component MLCK_Activation_CaM_md (Molar).
C STATES(19) is conc_CaMNCM in component MLCK_Activation_CaM_md (Molar).
C STATES(20) is conc_Ca_IP3R in component Calcium_Conc_SR (Molar).
C STATES(21) is conc_Ca_RyR in component Calcium_Conc_SR (Molar).
C STATES(22) is conc_CalseqCa_IP3R in component Calseq_Kinetics (Molar).
C STATES(23) is conc_CalseqCa_RyR in component Calseq_Kinetics (Molar).
C STATES(24) is V_m_cyt in component Membrane_Voltage (volt).
C STATES(25) is V_m_md in component Membrane_Voltage (volt).
C ALGBRC(35) is E_K_cyt in component Nernst_Potentials (volt).
C ALGBRC(40) is E_K_md in component Nernst_Potentials (volt).
C ALGBRC(43) is E_Na_cyt in component Nernst_Potentials (volt).
C ALGBRC(45) is E_Na_md in component Nernst_Potentials (volt).
C ALGBRC(47) is E_Cl_cyt in component Nernst_Potentials (volt).
C ALGBRC(49) is E_Cl_md in component Nernst_Potentials (volt).
C ALGBRC(51) is E_Ca_cyt in component Nernst_Potentials (volt).
C ALGBRC(53) is E_Ca_md in component Nernst_Potentials (volt).
C CONSTS(83) is conc_K_out in component Nernst_Potentials (Molar).
C CONSTS(84) is conc_Na_out in component Nernst_Potentials (Molar).
C CONSTS(85) is conc_Cl_out in component Nernst_Potentials (Molar).
C CONSTS(86) is conc_Ca_out in component Nernst_Potentials (Molar).
C CONSTS(87) is z_K in component Nernst_Potentials (dimensionless).
C CONSTS(88) is z_Na in component Nernst_Potentials (dimensionless).
C CONSTS(89) is z_Cl in component Nernst_Potentials (dimensionless).
C CONSTS(90) is z_Ca in component Nernst_Potentials (dimensionless).
C CONSTS(179) is RTF in component Nernst_Potentials (volt).
C ALGBRC(131) is I_SERCA_max in component ROS_SERCA_Interaction (ampere).
C ALGBRC(125) is K_mf in component SERCA_Activation (Molar).
C ALGBRC(55) is P_cyt_RyR in component RyR_Currents (dimensionless).
C ALGBRC(57) is P_md_RyR in component RyR_Currents (dimensionless).
C STATES(26) is w_cyt in component RyR_Currents (dimensionless).
C STATES(27) is w_md in component RyR_Currents (dimensionless).
C ALGBRC(1) is w_inf_cyt in component RyR_Currents (dimensionless).
C ALGBRC(2) is w_inf_md in component RyR_Currents (dimensionless).
C ALGBRC(63) is x110_cyt in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(67) is x110_md in component IP3R_Binding_Sites (dimensionless).
C STATES(28) is conc_IP3_cyt in component IP3R_Currents (Molar).
C STATES(29) is conc_IP3_md in component IP3R_Currents (Molar).
C STATES(30) is x000_cyt in component IP3R_Binding_Sites (dimensionless).
C STATES(31) is x000_md in component IP3R_Binding_Sites (dimensionless).
C STATES(32) is x001_cyt in component IP3R_Binding_Sites (dimensionless).
C STATES(33) is x001_md in component IP3R_Binding_Sites (dimensionless).
C STATES(34) is x010_cyt in component IP3R_Binding_Sites (dimensionless).
C STATES(35) is x010_md in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(65) is x011_cyt in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(69) is x011_md in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(61) is x101_cyt in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(62) is x101_md in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(66) is x111_cyt in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(70) is x111_md in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(59) is x100_cyt in component IP3R_Binding_Sites (dimensionless).
C ALGBRC(60) is x100_md in component IP3R_Binding_Sites (dimensionless).
C CONSTS(91) is d3 in component IP3R_Binding_Sites (Molar).
C ALGBRC(71) is conc_Ca_SR in component SOC_Current_Ca (Molar).
C CONSTS(92) is z_Na in component SOC_Current_Na (dimensionless).
C CONSTS(93) is z_Ca in component SOC_Current_Na (dimensionless).
C CONSTS(94) is P_SOC_ratio in component SOC_Current_Na (dimensionless).
C STATES(36) is y_C1_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(37) is y_C2_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(38) is y_C3_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(39) is y_O1_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(40) is y_O2_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(41) is y_O3_cyt in component Cl_Channels_cyt (dimensionless).
C ALGBRC(20) is Beta_cyt in component Cl_Channel_Rates (per_sec).
C ALGBRC(3) is y_C0_cyt in component Cl_Channels_cyt (dimensionless).
C STATES(42) is y_C1_md in component Cl_Channels_md (dimensionless).
C STATES(43) is y_C2_md in component Cl_Channels_md (dimensionless).
C STATES(44) is y_C3_md in component Cl_Channels_md (dimensionless).
C STATES(45) is y_O1_md in component Cl_Channels_md (dimensionless).
C STATES(46) is y_O2_md in component Cl_Channels_md (dimensionless).
C STATES(47) is y_O3_md in component Cl_Channels_md (dimensionless).
C ALGBRC(21) is Beta_md in component Cl_Channel_Rates (per_sec).
C ALGBRC(4) is y_C0_md in component Cl_Channels_md (dimensionless).
C CONSTS(95) is V1 in component Cl_Channel_Rates (dimensionless).
C CONSTS(96) is V2 in component Cl_Channel_Rates (per_V).
C CONSTS(97) is Lambda_Beta in component Cl_Channel_Rates (per_sec).
C ALGBRC(87) is x_Kir_cyt in component K_Inward_Rectifier_Current (dimensionless).
C ALGBRC(90) is x_Kir_md in component K_Inward_Rectifier_Current (dimensionless).
C ALGBRC(84) is alpha_Kir_cyt in component K_ir_Rates (per_sec).
C ALGBRC(85) is alpha_Kir_md in component K_ir_Rates (per_sec).
C ALGBRC(86) is Beta_Kir_cyt in component K_ir_Rates (per_sec).
C ALGBRC(89) is Beta_Kir_md in component K_ir_Rates (per_sec).
C ALGBRC(92) is P_Kv_cyt in component Kv_Activations (dimensionless).
C ALGBRC(94) is P_Kv_md in component Kv_Activations (dimensionless).
C STATES(48) is P1_cyt in component Kv_Activations (dimensionless).
C STATES(49) is P1_md in component Kv_Activations (dimensionless).
C STATES(50) is P2_cyt in component Kv_Activations (dimensionless).
C STATES(51) is P2_md in component Kv_Activations (dimensionless).
C ALGBRC(5) is P_bar_Kv_cyt in component Kv_Activations (dimensionless).
C ALGBRC(6) is P_bar_Kv_md in component Kv_Activations (dimensionless).
C ALGBRC(22) is Tau_P1_cyt in component Kv_Time_Constants (second).
C ALGBRC(23) is Tau_P1_md in component Kv_Time_Constants (second).
C ALGBRC(24) is Tau_P2_cyt in component Kv_Time_Constants (second).
C ALGBRC(25) is Tau_P2_md in component Kv_Time_Constants (second).
C ALGBRC(98) is P_KCa_cyt in component KCa_Activations (dimensionless).
C ALGBRC(100) is P_KCa_md in component KCa_Activations (dimensionless).
C STATES(52) is PF_cyt in component KCa_Activations (dimensionless).
C STATES(53) is PF_md in component KCa_Activations (dimensionless).
C STATES(54) is PS_cyt in component KCa_Activations (dimensionless).
C STATES(55) is PS_md in component KCa_Activations (dimensionless).
C ALGBRC(26) is P_bar_KCa_cyt in component KCa_Activations (dimensionless).
C ALGBRC(27) is P_bar_KCa_md in component KCa_Activations (dimensionless).
C ALGBRC(7) is V_KCa_cyt in component cGMP_Ca_Interaction (volt).
C ALGBRC(8) is V_KCa_md in component cGMP_Ca_Interaction (volt).
C ALGBRC(139) is I_NaK_alpha1_max in component ROS_NaK_Interaction (ampere).
C ALGBRC(143) is I_NaK_alpha2_max in component NaK_alpha_Current (ampere).
C ALGBRC(102) is Psi_NaK_cyt in component NaK_alpha_Waveform (dimensionless).
C ALGBRC(103) is Psi_NaK_md in component NaK_alpha_Waveform (dimensionless).
C CONSTS(180) is Sigma in component NaK_alpha_Waveform (dimensionless).
C ALGBRC(104) is K_a_NaCa in component NCX_Waveform (dimensionless).
C ALGBRC(105) is Psi_F in component NCX_Waveform (dimensionless).
C ALGBRC(106) is Psi_R in component NCX_Waveform (dimensionless).
C ALGBRC(107) is G in component NCX_Waveform (M_fourpow).
C STATES(56) is m_VONa_cyt in component VONa_Channels (dimensionless).
C STATES(57) is m_VONa_md in component VONa_Channels (dimensionless).
C STATES(58) is h_VONa_cyt in component VONa_Channels (dimensionless).
C STATES(59) is h_VONa_md in component VONa_Channels (dimensionless).
C ALGBRC(9) is m_bar_cyt in component VONa_Channels (dimensionless).
C ALGBRC(10) is m_bar_md in component VONa_Channels (dimensionless).
C ALGBRC(11) is h_bar_cyt in component VONa_Channels (dimensionless).
C ALGBRC(12) is h_bar_md in component VONa_Channels (dimensionless).
C ALGBRC(126) is I_CaP_max in component CaP_Current_Max (ampere).
C STATES(60) is d_L_cyt in component CaL_Activations (dimensionless).
C STATES(61) is d_L_md in component CaL_Activations (dimensionless).
C ALGBRC(111) is f_L_cyt in component CaL_Activations (dimensionless).
C ALGBRC(113) is f_L_md in component CaL_Activations (dimensionless).
C STATES(62) is f_F_cyt in component CaL_Activations (dimensionless).
C STATES(63) is f_F_md in component CaL_Activations (dimensionless).
C ALGBRC(13) is d_bar_L_cyt in component CaL_Activations (dimensionless).
C ALGBRC(14) is d_bar_L_md in component CaL_Activations (dimensionless).
C ALGBRC(15) is f_bar_F_cyt in component CaL_Activations (dimensionless).
C ALGBRC(16) is f_bar_F_md in component CaL_Activations (dimensionless).
C ALGBRC(28) is Tau_d_cyt in component CaL_Time_Constants (second).
C ALGBRC(29) is Tau_d_md in component CaL_Time_Constants (second).
C ALGBRC(30) is Tau_f_cyt in component CaL_Time_Constants (second).
C ALGBRC(31) is Tau_f_md in component CaL_Time_Constants (second).
C ALGBRC(115) is Xi_K in component Md_Cyt_Flux_K (dimensionless).
C CONSTS(183) is P_K_diff in component Md_Cyt_Flux_K (m_per_sec).
C CONSTS(98) is z_K in component Md_Cyt_Flux_K (dimensionless).
C CONSTS(187) is Volsec in component Md_Cyt_Flux_K (cubic_m_per_sec).
C CONSTS(99) is DiffArea in component Md_Cyt_Flux_K (sq_meter).
C ALGBRC(117) is Xi_Na in component Md_Cyt_Flux_Na (dimensionless).
C CONSTS(184) is P_Na_diff in component Md_Cyt_Flux_Na (m_per_sec).
C CONSTS(100) is z_Na in component Md_Cyt_Flux_Na (dimensionless).
C CONSTS(188) is Volsec in component Md_Cyt_Flux_Na (cubic_m_per_sec).
C CONSTS(101) is DiffArea in component Md_Cyt_Flux_Na (sq_meter).
C ALGBRC(119) is Xi_Ca in component Md_Cyt_Flux_Ca (dimensionless).
C CONSTS(181) is P_Ca_diff in component Md_Cyt_Flux_Ca (m_per_sec).
C CONSTS(102) is z_Ca in component Md_Cyt_Flux_Ca (dimensionless).
C CONSTS(185) is Volsec in component Md_Cyt_Flux_Ca (cubic_m_per_sec).
C CONSTS(103) is DiffArea in component Md_Cyt_Flux_Ca (sq_meter).
C ALGBRC(121) is Xi_Cl in component Md_Cyt_Flux_Cl (dimensionless).
C CONSTS(182) is P_Cl_diff in component Md_Cyt_Flux_Cl (m_per_sec).
C CONSTS(104) is z_Cl in component Md_Cyt_Flux_Cl (dimensionless).
C CONSTS(186) is Volsec in component Md_Cyt_Flux_Cl (cubic_m_per_sec).
C CONSTS(105) is DiffArea in component Md_Cyt_Flux_Cl (sq_meter).
C ALGBRC(144) is I_cyt_tot in component Membrane_Potential_Cyt (ampere).
C ALGBRC(136) is I_cyt_up in component Membrane_Potential_Cyt (ampere).
C ALGBRC(145) is I_md_tot in component Membrane_Potential_md (ampere).
C CONSTS(106) is C_m in component Membrane_Voltage (farad2).
C ALGBRC(138) is I_md_up in component Membrane_Potential_md (ampere).
C CONSTS(107) is k_NO_O2 in component Parameters_VR (per_M_sec).
C CONSTS(108) is k_NO_cons in component Parameters_VR (per_sec).
C CONSTS(109) is k_SOD_O2 in component Parameters_VR (per_M_sec).
C CONSTS(110) is k_cat in component Parameters_VR (per_M_sec).
C CONSTS(111) is conc_SOD in component Parameters_VR (Molar).
C CONSTS(112) is conc_CAT in component Parameters_VR (Molar).
C CONSTS(113) is J_NO in component Parameters_VR (mol_per_sq_m_sec).
C CONSTS(114) is G_O2 in component Parameters_VR (M_per_sec).
C CONSTS(115) is J_H2O2 in component Parameters_VR (mol_per_sq_m_sec).
C CONSTS(116) is P_NO in component Parameters_VR (m_per_sec).
C CONSTS(117) is P_O2 in component Parameters_VR (m_per_sec).
C CONSTS(118) is P_H2O2 in component Parameters_VR (m_per_sec).
C CONSTS(119) is k_sGC_for in component Parameters_VR (per_M_sec).
C CONSTS(120) is k_sGC_back in component Parameters_VR (per_sec).
C CONSTS(121) is k2_sGC_for in component Parameters_VR (per_sec).
C CONSTS(122) is k3_sCG_for in component Parameters_VR (per_M_sec).
C CONSTS(123) is K4_sGC_for in component Parameters_VR (per_M_sec).
C CONSTS(124) is V_sGC_max in component Parameters_VR (M_per_sec).
C CONSTS(125) is k_PDE in component Parameters_VR (per_sec).
C CONSTS(126) is K_M_PDE in component Parameters_VR (Molar).
C CONSTS(127) is K_MLCP in component Parameters_VR (Molar).
C CONSTS(128) is K_SERCA in component Parameters_VR (Molar).
C CONSTS(129) is K_PMCA in component Parameters_VR (Molar).
C CONSTS(130) is K_KCa_cGMP in component Parameters_VR (Molar).
C CONSTS(131) is V_KCa_cGMP in component Parameters_VR (volt).
C CONSTS(132) is K_KCa_NO in component Parameters_VR (Molar).
C CONSTS(133) is V_KCa_NO in component Parameters_VR (volt).
C CONSTS(134) is V_Ca in component Parameters_VR (volt).
C CONSTS(135) is V_B in component Parameters_VR (volt).
C CONSTS(136) is k_KCa in component Parameters_VR (volt).
C CONSTS(137) is K_SERCA_O2 in component Parameters_VR (Molar).
C CONSTS(138) is K_SERCA_H2O2 in component Parameters_VR (Molar).
C CONSTS(139) is K_NaK_O2 in component Parameters_VR (Molar).
C CONSTS(140) is K_NaK_H202 in component Parameters_VR (Molar).
C CONSTS(141) is K_MLCP_H2O2 in component Parameters_VR (Molar).
C CONSTS(142) is alpha_MLCP in component Parameters_VR (dimensionless).
C CONSTS(143) is k1_CAM_on in component Parameters_VR (per_sq_M_sec).
C CONSTS(144) is k1_CAM_off in component Parameters_VR (per_sec).
C CONSTS(145) is k2_CAM_on in component Parameters_VR (per_sq_M_sec).
C CONSTS(146) is k2_CAM_off in component Parameters_VR (per_sec).
C CONSTS(147) is k3_CAM_on in component Parameters_VR (per_sq_M_sec).
C CONSTS(148) is k3_CAM_off in component Parameters_VR (per_sec).
C CONSTS(149) is k4_CAM_on in component Parameters_VR (per_sq_M_sec).
C CONSTS(150) is k4_CAM_off in component Parameters_VR (per_sec).
C CONSTS(151) is k5_CAM_on in component Parameters_VR (per_M_sec).
C CONSTS(152) is k5_CAM_off in component Parameters_VR (per_sec).
C CONSTS(153) is k6_CAM_on in component Parameters_VR (per_sq_M_sec).
C CONSTS(154) is k6_CAM_off in component Parameters_VR (per_sec).
C CONSTS(155) is k7_CAM_on in component Parameters_VR (per_M_sec).
C CONSTS(156) is k7_CAM_off in component Parameters_VR (per_sec).
C CONSTS(157) is k_Myo_MLCK in component Parameters_VR (per_M_sec).
C CONSTS(158) is k_Myo_MLCP in component Parameters_VR (per_sec).
C CONSTS(159) is k3_Myo in component Parameters_VR (per_sec).
C CONSTS(160) is k4_Myo in component Parameters_VR (per_sec).
C CONSTS(161) is k7_Myo in component Parameters_VR (per_sec).
C CONSTS(162) is conc_CaM_tot in component Parameters_VR (Molar).
C CONSTS(163) is conc_MLCK_tot in component Parameters_VR (Molar).
C CONSTS(164) is conc_Myo_tot in component Parameters_VR (Molar).
C CONSTS(165) is K_mf_rest in component Parameters_VR (Molar).
C CONSTS(166) is I_PMCA_rest_max in component Parameters_VR (ampere).
C CONSTS(167) is I_SERCA_rest_max in component Parameters_VR (ampere).
C CONSTS(168) is I_NaK_alpha1_rest_max in component Parameters_VR (ampere).
C CONSTS(169) is I_NaK_alpha2_rest_max in component Parameters_VR (ampere).
C ALGBRC(123) is conc_CaMNC_ in component MLCK_Activation_CaM (Molar).
C ALGBRC(124) is conc_MLCK_free in component MLCK_Activation_CaM (Molar).
C STATES(64) is conc_Myo in component MLCK_Phospho_Myosin (Molar).
C STATES(65) is conc_MyoP in component MLCK_Phospho_Myosin (Molar).
C STATES(66) is conc_AMyo in component MLCK_Phospho_Myosin (Molar).
C ALGBRC(17) is conc_AMyoP in component MLCK_Phospho_Myosin (Molar).
C ALGBRC(32) is k1_6_Myo in component MLCP_Rate_Constants (per_sec).
C ALGBRC(39) is k2_Myo in component MLCP_Rate_Constants (per_sec).
C ALGBRC(41) is k5_Myo in component MLCP_Rate_Constants (per_sec).
C CONSTS(170) is K_M_cGMP in component MLCP_Rate_Constants (Molar).
C STATES(67) is conc_cGMP in component cGMP_Kinetics (Molar).
C ALGBRC(36) is R_MLCP in component MLCP_Rate_Constants (dimensionless).
C STATES(68) is conc_H2O2_cyt in component H2O2_Kinetics (Molar).
C STATES(69) is conc_NO_cyt in component NO_Kinetics (Molar).
C STATES(70) is conc_O2_cyt in component O2_Kinetics (Molar).
C ALGBRC(18) is zeta in component Zeta (dimensionless).
C STATES(71) is conc_Eb in component sGC_Activation (dimensionless).
C ALGBRC(19) is conc_E6c in component sGC_Activation (dimensionless).
C STATES(72) is conc_E5c in component sGC_Activation (dimensionless).
C ALGBRC(33) is k4_sGC_for in component sGC_Activation (per_sec).
C CONSTS(171) is convert in component cGMP_Ca_Interaction (dimensionless).
C CONSTS(172) is K_mod in component SERCA_Activation (dimensionless).
C ALGBRC(129) is theta_SERCA_O2 in component ROS_SERCA_Interaction (dimensionless).
C ALGBRC(130) is theta_SERCA_H2O2 in component ROS_SERCA_Interaction (dimensionless).
C STATES(73) is conc_CaM___ in component MLCK_Activation_CaM_md (Molar).
C ALGBRC(137) is conc_CaMNC_ in component MLCK_Activation_CaM_md (Molar).
C ALGBRC(140) is conc_MLCK_free in component MLCK_Activation_CaM_md (Molar).
C ALGBRC(34) is f_p in component Contractile_Force (dimensionless).
C ALGBRC(37) is F_contract in component Contractile_Force (dimensionless).
C CONSTS(173) is a in component Contractile_Force (dimensionless).
C CONSTS(174) is b in component Contractile_Force (dimensionless).
C CONSTS(175) is c in component Contractile_Force (dimensionless).
C CONSTS(176) is d in component Contractile_Force (dimensionless).
C RATES(1) is d/dt conc_K_cyt in component Potassium_Conc_cyt (Molar).
C RATES(2) is d/dt conc_K_md in component Potassium_Conc_md (Molar).
C RATES(3) is d/dt conc_Na_cyt in component Sodium_Conc_cyt (Molar).
C RATES(4) is d/dt conc_Na_md in component Sodium_Conc_md (Molar).
C RATES(5) is d/dt conc_Cl_cyt in component Chloride_Conc_cyt (Molar).
C RATES(6) is d/dt conc_Cl_md in component Chloride_Conc_md (Molar).
C RATES(7) is d/dt conc_Ca_cyt in component Calcium_Conc_cyt (Molar).
C RATES(8) is d/dt conc_CaMN__ in component MLCK_Activation_CaM (Molar).
C RATES(10) is d/dt conc_CaM_C_ in component MLCK_Activation_CaM (Molar).
C RATES(11) is d/dt conc_CaM_CM in component MLCK_Activation_CaM (Molar).
C RATES(9) is d/dt conc_BfCa_cyt in component Bf_Kinetics (Molar).
C RATES(14) is d/dt conc_BfCa_md in component Bf_Kinetics (Molar).
C RATES(15) is d/dt conc_Ca_md in component Calcium_Conc_md (Molar).
C RATES(16) is d/dt conc_CaMN__ in component MLCK_Activation_CaM_md (Molar).
C RATES(17) is d/dt conc_CaM_C_ in component MLCK_Activation_CaM_md (Molar).
C RATES(18) is d/dt conc_CaM_CM in component MLCK_Activation_CaM_md (Molar).
C RATES(20) is d/dt conc_Ca_IP3R in component Calcium_Conc_SR (Molar).
C RATES(22) is d/dt conc_CalseqCa_IP3R in component Calseq_Kinetics (Molar).
C RATES(21) is d/dt conc_Ca_RyR in component Calcium_Conc_SR (Molar).
C RATES(23) is d/dt conc_CalseqCa_RyR in component Calseq_Kinetics (Molar).
C RATES(26) is d/dt w_cyt in component RyR_Currents (dimensionless).
C RATES(27) is d/dt w_md in component RyR_Currents (dimensionless).
C RATES(28) is d/dt conc_IP3_cyt in component IP3R_Currents (Molar).
C RATES(29) is d/dt conc_IP3_md in component IP3R_Currents (Molar).
C RATES(30) is d/dt x000_cyt in component IP3R_Binding_Sites (dimensionless).
C RATES(31) is d/dt x000_md in component IP3R_Binding_Sites (dimensionless).
C RATES(32) is d/dt x001_cyt in component IP3R_Binding_Sites (dimensionless).
C RATES(33) is d/dt x001_md in component IP3R_Binding_Sites (dimensionless).
C RATES(34) is d/dt x010_cyt in component IP3R_Binding_Sites (dimensionless).
C RATES(35) is d/dt x010_md in component IP3R_Binding_Sites (dimensionless).
C RATES(36) is d/dt y_C1_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(37) is d/dt y_C2_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(38) is d/dt y_C3_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(39) is d/dt y_O1_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(40) is d/dt y_O2_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(41) is d/dt y_O3_cyt in component Cl_Channels_cyt (dimensionless).
C RATES(42) is d/dt y_C1_md in component Cl_Channels_md (dimensionless).
C RATES(43) is d/dt y_C2_md in component Cl_Channels_md (dimensionless).
C RATES(44) is d/dt y_C3_md in component Cl_Channels_md (dimensionless).
C RATES(45) is d/dt y_O1_md in component Cl_Channels_md (dimensionless).
C RATES(46) is d/dt y_O2_md in component Cl_Channels_md (dimensionless).
C RATES(47) is d/dt y_O3_md in component Cl_Channels_md (dimensionless).
C RATES(48) is d/dt P1_cyt in component Kv_Activations (dimensionless).
C RATES(49) is d/dt P1_md in component Kv_Activations (dimensionless).
C RATES(50) is d/dt P2_cyt in component Kv_Activations (dimensionless).
C RATES(51) is d/dt P2_md in component Kv_Activations (dimensionless).
C RATES(52) is d/dt PF_cyt in component KCa_Activations (dimensionless).
C RATES(53) is d/dt PF_md in component KCa_Activations (dimensionless).
C RATES(54) is d/dt PS_cyt in component KCa_Activations (dimensionless).
C RATES(55) is d/dt PS_md in component KCa_Activations (dimensionless).
C RATES(56) is d/dt m_VONa_cyt in component VONa_Channels (dimensionless).
C RATES(57) is d/dt m_VONa_md in component VONa_Channels (dimensionless).
C RATES(58) is d/dt h_VONa_cyt in component VONa_Channels (dimensionless).
C RATES(59) is d/dt h_VONa_md in component VONa_Channels (dimensionless).
C RATES(60) is d/dt d_L_cyt in component CaL_Activations (dimensionless).
C RATES(61) is d/dt d_L_md in component CaL_Activations (dimensionless).
C RATES(62) is d/dt f_F_cyt in component CaL_Activations (dimensionless).
C RATES(63) is d/dt f_F_md in component CaL_Activations (dimensionless).
C RATES(24) is d/dt V_m_cyt in component Membrane_Voltage (volt).
C RATES(25) is d/dt V_m_md in component Membrane_Voltage (volt).
C RATES(13) is d/dt conc_CaM___ in component MLCK_Activation_CaM (Molar).
C RATES(12) is d/dt conc_CaMNCM in component MLCK_Activation_CaM (Molar).
C RATES(64) is d/dt conc_Myo in component MLCK_Phospho_Myosin (Molar).
C RATES(65) is d/dt conc_MyoP in component MLCK_Phospho_Myosin (Molar).
C RATES(66) is d/dt conc_AMyo in component MLCK_Phospho_Myosin (Molar).
C RATES(69) is d/dt conc_NO_cyt in component NO_Kinetics (Molar).
C RATES(70) is d/dt conc_O2_cyt in component O2_Kinetics (Molar).
C RATES(68) is d/dt conc_H2O2_cyt in component H2O2_Kinetics (Molar).
C RATES(71) is d/dt conc_Eb in component sGC_Activation (dimensionless).
C RATES(72) is d/dt conc_E5c in component sGC_Activation (dimensionless).
C RATES(67) is d/dt conc_cGMP in component cGMP_Kinetics (Molar).
C RATES(73) is d/dt conc_CaM___ in component MLCK_Activation_CaM_md (Molar).
C RATES(19) is d/dt conc_CaMNCM in component MLCK_Activation_CaM_md (Molar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 96487
CONSTS(2) = 5E-16
CONSTS(3) = 3E-18
CONSTS(4) = 7E-17
CONSTS(5) = 3.5E-16
CONSTS(6) = 1.7
CONSTS(7) = 2
CONSTS(8) = 12
CONSTS(9) = 3.7224E-4
CONSTS(10) = 6.3601E-4
CONSTS(11) = 5.71E-5
CONSTS(12) = 0.1
CONSTS(13) = 20
CONSTS(14) = 1.85
CONSTS(15) = 2.4E-4
CONSTS(16) = 0.5
CONSTS(17) = 1.1E-3
CONSTS(18) = 1
CONSTS(19) = 1.3E-4
CONSTS(20) = 2E2
CONSTS(21) = 1.445E-4
CONSTS(22) = 2E4
CONSTS(23) = 8.234E-5
CONSTS(24) = 7.5E-11
CONSTS(25) = 0.1
CONSTS(26) = 8.78E-2
CONSTS(27) = 6E-17
CONSTS(28) = 3.5E-18
CONSTS(29) = 6E-17
CONSTS(30) = 46
CONSTS(31) = 4.0
CONSTS(32) = 78
CONSTS(33) = 78
CONSTS(34) = 3.0E3
CONSTS(35) = 50
CONSTS(36) = 10
CONSTS(37) = 30
CONSTS(38) = 300
CONSTS(39) = 1.2E-8
CONSTS(40) = 1E-11
CONSTS(41) = 1E-11
CONSTS(42) = 1.2E-11
CONSTS(43) = 1E-11
CONSTS(44) = 3E-10
CONSTS(45) = 5.4
CONSTS(46) = 9.8325E-9
CONSTS(47) = 0.463
CONSTS(48) = 1.036E-9
CONSTS(49) = 5.00E-10
CONSTS(50) = 5E-4
CONSTS(51) = 1.15E-2
CONSTS(52) = 1.5
CONSTS(53) = 12
CONSTS(54) = 22
CONSTS(55) = 1.7097E-10
CONSTS(56) = 3.59E-3
CONSTS(57) = 1.3
CONSTS(58) = 12.29
CONSTS(59) = 87.5
CONSTS(60) = 0.27
CONSTS(61) = 0.35
CONSTS(62) = 1.25E-4
CONSTS(63) = 2.860E-9
CONSTS(64) = 0.0001
CONSTS(65) = 0.001
CONSTS(66) = 1.7E-4
CONSTS(67) = 1.4151E-9
CONSTS(68) = 4E-1
CONSTS(69) = 100E3
CONSTS(70) = 300
CONSTS(71) = 1.4E-1
CONSTS(72) = 1E5
CONSTS(73) = 8E4
CONSTS(74) = 0.8580
CONSTS(75) = 0.1420
CONSTS(76) = 8.3144621
CONSTS(77) = 298
CONSTS(78) = 1.21E-9
CONSTS(79) = 5E-7
CONSTS(80) = 0.001
CONSTS(81) = 3E-10
CONSTS(82) = 3E-10
STATES(1) = 140.55
STATES(2) = 131.37
STATES(3) = 9.1964
STATES(4) = 20.742
STATES(5) = 23.257
STATES(6) = 16.737
STATES(7) = 1.0533E-4
STATES(8) = 1.1023E-5
STATES(9) = 1.3568E-2
STATES(10) = 4.1154E-5
STATES(11) = 1.9451E-3
STATES(12) = 5.3950E-5
STATES(13) = 7.9487E-3
STATES(14) = 4.4109E-2
STATES(15) = 3.7182E-4
STATES(16) = 1.2766E-4
STATES(17) = 4.7661E-4
STATES(18) = 1.4863E-3
STATES(19) = 5.1368E-4
STATES(20) = 5.6334E-1
STATES(21) = 4.9091E-1
STATES(22) = 5.7849E-2
STATES(23) = 5.3239E-2
STATES(24) = -7.8960E-2
STATES(25) = -7.6961E-2
CONSTS(83) = 5.4
CONSTS(84) = 140
CONSTS(85) = 110
CONSTS(86) = 2.0
CONSTS(87) = 1
CONSTS(88) = 1
CONSTS(89) = -1
CONSTS(90) = 2
STATES(26) = 0.11187
STATES(27) = 0.89971
STATES(28) = 2.4140E-4
STATES(29) = 2.7808E-4
STATES(30) = 0.11631
STATES(31) = 0.028019
STATES(32) = 0.084779
STATES(33) = 0.072097
STATES(34) = 0.14878
STATES(35) = 0.12652
CONSTS(91) = 9.434E-4
CONSTS(92) = 1
CONSTS(93) = 2
CONSTS(94) = 0.125
STATES(36) = 0.0062437
STATES(37) = 0.000039459
STATES(38) = 0.00000024937
STATES(39) = 0.0056570
STATES(40) = 0.00010725
STATES(41) = 0.0000067783
STATES(42) = 2.1346E-2
STATES(43) = 4.7620E-4
STATES(44) = 1.0624E-5
STATES(45) = 1.9731E-2
STATES(46) = 1.3205E-3
STATES(47) = 2.94594E-4
CONSTS(95) = -3.8
CONSTS(96) = -10.0
CONSTS(97) = 224
STATES(48) = 0.0048876
STATES(49) = 0.0056051
STATES(50) = 0.0048876
STATES(51) = 0.0056051
STATES(52) = 0.94813
STATES(53) = 0.67538
STATES(54) = 0.94813
STATES(55) = 0.67538
STATES(56) = 1.0429E-5
STATES(57) = 1.5817E-5
STATES(58) = 0.87853
STATES(59) = 0.84170
STATES(60) = 3.7841E-5
STATES(61) = 4.9277E-5
STATES(62) = 1
STATES(63) = 1
CONSTS(98) = 1
CONSTS(99) = 7.7598E-13
CONSTS(100) = 1
CONSTS(101) = 7.7598E-13
CONSTS(102) = 2
CONSTS(103) = 7.7598E-13
CONSTS(104) = -1
CONSTS(105) = 7.7598E-13
CONSTS(106) = 1.2100E-11
CONSTS(107) = 1.9E7
CONSTS(108) = 0.01
CONSTS(109) = 1.6E6
CONSTS(110) = 3.4E4
CONSTS(111) = 1E-3
CONSTS(112) = 20E-3
CONSTS(113) = 2.8E-8
CONSTS(114) = 4.5E-3
CONSTS(115) = 6E-8
CONSTS(116) = 3E-4
CONSTS(117) = 6E-4
CONSTS(118) = 6E-4
CONSTS(119) = 2E6
CONSTS(120) = 100
CONSTS(121) = 0.1
CONSTS(122) = 3E3
CONSTS(123) = 98
CONSTS(124) = 1.09E-3
CONSTS(125) = 0.032
CONSTS(126) = 2.0E-3
CONSTS(127) = 5.5E-3
CONSTS(128) = 5.0E-4
CONSTS(129) = 1.0E-3
CONSTS(130) = 5.5E-4
CONSTS(131) = 66.9
CONSTS(132) = 2E-1
CONSTS(133) = 100
CONSTS(134) = -53.7
CONSTS(135) = 283.7
CONSTS(136) = 30.8E-3
CONSTS(137) = 2.0E-2
CONSTS(138) = 6.6E-2
CONSTS(139) = 1.6E-3
CONSTS(140) = 0.56
CONSTS(141) = 1E-3
CONSTS(142) = 4.0
CONSTS(143) = 2.8E6
CONSTS(144) = 6
CONSTS(145) = 100E6
CONSTS(146) = 800
CONSTS(147) = 2.8E6
CONSTS(148) = 6.0
CONSTS(149) = 100E6
CONSTS(150) = 800
CONSTS(151) = 1000E3
CONSTS(152) = 20
CONSTS(153) = 12.5E6
CONSTS(154) = 5.0
CONSTS(155) = 1000E3
CONSTS(156) = 1
CONSTS(157) = 4E3
CONSTS(158) = 0.4
CONSTS(159) = 1.8
CONSTS(160) = 0.1
CONSTS(161) = 0.045
CONSTS(162) = 1.00E-2
CONSTS(163) = 2.0E-3
CONSTS(164) = 30E-3
CONSTS(165) = 3.0995E-4
CONSTS(166) = 2.853E-12
CONSTS(167) = 100E-12
CONSTS(168) = 2.7225E-11
CONSTS(169) = 2.7225E-11
STATES(64) = 8.0265E-3
STATES(65) = 6.3288E-4
STATES(66) = 1.8089E-2
CONSTS(170) = 5.5E-3
STATES(67) = 9.0341E-3
STATES(68) = 5.0006E-5
STATES(69) = 9.2798E-5
STATES(70) = 8.3402E-7
STATES(71) = 0.27479
STATES(72) = 0.21715
CONSTS(171) = 1E-3
CONSTS(172) = 0.5263
STATES(73) = 7.3876E-3
CONSTS(173) = 3.9796
CONSTS(174) = 97.35
CONSTS(175) = 522.25
CONSTS(176) = 329.81
CONSTS(177) = ( CONSTS(81)*1.96000)/0.790000
CONSTS(178) = ( CONSTS(81)*1.33000)/0.790000
CONSTS(179) = ( CONSTS(76)*CONSTS(77))/CONSTS(1)
CONSTS(180) = (1.00000/7.00000)*(EXP(CONSTS(84)/67.3000) - 1.00000)
CONSTS(181) = ( CONSTS(80)*CONSTS(81))/CONSTS(79)
CONSTS(182) = ( CONSTS(80)*CONSTS(82))/CONSTS(79)
CONSTS(183) = ( CONSTS(80)*CONSTS(177))/CONSTS(79)
CONSTS(184) = ( CONSTS(80)*CONSTS(178))/CONSTS(79)
CONSTS(185) = CONSTS(103)*CONSTS(181)
CONSTS(186) = CONSTS(105)*CONSTS(182)
CONSTS(187) = CONSTS(99)*CONSTS(183)
CONSTS(188) = CONSTS(101)*CONSTS(184)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(9) = CONSTS(69)*STATES(7)*(CONSTS(68) - STATES(9)) - CONSTS(70)*STATES(9)
RATES(28) = CONSTS(14)*CONSTS(15)*((STATES(7)+ (1.00000 - CONSTS(16))*CONSTS(17))/(STATES(7)+CONSTS(17))) - CONSTS(18)*STATES(28)
RATES(29) = CONSTS(14)*CONSTS(15)*((STATES(15)+ (1.00000 - CONSTS(16))*CONSTS(17))/(STATES(15)+CONSTS(17))) - CONSTS(18)*STATES(29)
RATES(30) = - CONSTS(20)*( STATES(7)*STATES(30) - CONSTS(21)*STATES(32)) - CONSTS(22)*( STATES(7)*STATES(30) - CONSTS(23)*STATES(34))
RATES(31) = - CONSTS(20)*( STATES(15)*STATES(31) - CONSTS(21)*STATES(33)) - CONSTS(22)*( STATES(15)*STATES(31) - CONSTS(23)*STATES(35))
RATES(13) = (( - CONSTS(143)*STATES(7) ** 2.00000*STATES(13)+ CONSTS(144)*STATES(10)) - CONSTS(145)*STATES(7) ** 2.00000*STATES(13))+ CONSTS(146)*STATES(8)
RATES(69) = ( (CONSTS(78)/( 2.00000*CONSTS(2)))*(CONSTS(113) - CONSTS(116)*STATES(69)) - CONSTS(107)*STATES(69)*STATES(70)) - CONSTS(108)*STATES(69)
RATES(68) = ( (CONSTS(78)/( 2.00000*CONSTS(2)))*(CONSTS(115) - CONSTS(118)*STATES(68))+ CONSTS(109)*CONSTS(111)*STATES(70)) - CONSTS(110)*STATES(68)*CONSTS(112)
RATES(67) = CONSTS(124)*STATES(72) - ( CONSTS(125)*STATES(67) ** 2.00000)/(CONSTS(126)+STATES(67))
RATES(73) = (( - CONSTS(143)*STATES(15) ** 2.00000*STATES(73)+ CONSTS(144)*STATES(17)) - CONSTS(145)*STATES(15) ** 2.00000*STATES(73))+ CONSTS(146)*STATES(16)
RATES(14) = CONSTS(69)*STATES(15)*(CONSTS(68) - STATES(14)) - CONSTS(70)*STATES(14)
ALGBRC(1) = (1.00000+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)/(1.00000+1.00000/CONSTS(11)+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)
RATES(26) = ( CONSTS(12)*(ALGBRC(1) - STATES(26)))/ALGBRC(1)
ALGBRC(2) = (1.00000+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)/(1.00000+1.00000/CONSTS(11)+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)
RATES(27) = ( CONSTS(12)*(ALGBRC(2) - STATES(27)))/ALGBRC(2)
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(24)+0.0239000)/0.00480000))
RATES(56) = (ALGBRC(9) - STATES(56))/CONSTS(64)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(25)+0.0239000)/0.00480000))
RATES(57) = (ALGBRC(10) - STATES(57))/CONSTS(64)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(24)+0.0661000)/0.00650000))
RATES(58) = (ALGBRC(11) - STATES(58))/CONSTS(65)
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(25)+0.0661000)/0.00650000))
RATES(59) = (ALGBRC(12) - STATES(59))/CONSTS(65)
ALGBRC(18) = ( - CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))
RATES(70) = ((CONSTS(114) - (CONSTS(78)/( 2.00000*CONSTS(2)))*CONSTS(117)*STATES(70)*(ALGBRC(18)/(1.00000 - EXP(- ALGBRC(18))))) - CONSTS(107)*STATES(69)*STATES(70)) - CONSTS(109)*CONSTS(111)*STATES(70)
ALGBRC(20) = CONSTS(97)*EXP(CONSTS(95)+ CONSTS(96)*STATES(24))
CALL minimize(minfunc_0, CONSTS, VARIABLES, ALGBRC(3))
RATES(36) = (( - CONSTS(34)*STATES(7)*(STATES(36) - ALGBRC(3))+ CONSTS(35)*(STATES(37) - STATES(36))) - CONSTS(36)*STATES(36))+ ALGBRC(20)*STATES(39)
RATES(37) = (( - CONSTS(34)*STATES(7)*(STATES(37) - STATES(36))+ CONSTS(35)*(STATES(38) - STATES(37))) - CONSTS(37)*STATES(37))+ ALGBRC(20)*STATES(40)
RATES(38) = (( CONSTS(34)*STATES(7)*STATES(37) - CONSTS(35)*STATES(38)) - CONSTS(38)*STATES(38))+ ALGBRC(20)*STATES(41)
RATES(39) = CONSTS(36)*STATES(36) - ALGBRC(20)*STATES(39)
RATES(40) = CONSTS(37)*STATES(37) - ALGBRC(20)*STATES(40)
RATES(41) = CONSTS(38)*STATES(38) - ALGBRC(20)*STATES(41)
ALGBRC(21) = CONSTS(97)*EXP(CONSTS(95)+ CONSTS(96)*STATES(25))
CALL minimize(minfunc_1, CONSTS, VARIABLES, ALGBRC(4))
RATES(42) = (( - CONSTS(34)*STATES(15)*(STATES(42) - ALGBRC(4))+ CONSTS(35)*(STATES(43) - STATES(42))) - CONSTS(36)*STATES(42))+ ALGBRC(21)*STATES(45)
RATES(43) = (( - CONSTS(34)*STATES(15)*(STATES(43) - STATES(42))+ CONSTS(35)*(STATES(44) - STATES(43))) - CONSTS(37)*STATES(43))+ ALGBRC(21)*STATES(46)
RATES(44) = (( CONSTS(34)*STATES(15)*STATES(43) - CONSTS(35)*STATES(44)) - CONSTS(38)*STATES(44))+ ALGBRC(21)*STATES(47)
RATES(45) = CONSTS(36)*STATES(42) - ALGBRC(21)*STATES(45)
RATES(46) = CONSTS(37)*STATES(43) - ALGBRC(21)*STATES(46)
RATES(47) = CONSTS(38)*STATES(44) - ALGBRC(21)*STATES(47)
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(24)+0.00177000)/0.0145200))
ALGBRC(22) = 0.210987*EXP(- (STATES(24)+0.214340)/0.195350 ** 2.00000) - 0.0205900
RATES(48) = (ALGBRC(5) - STATES(48))/ALGBRC(22)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(25)+0.00177000)/0.0145200))
ALGBRC(23) = 0.210987*EXP(- (STATES(25)+0.214340)/0.195350 ** 2.00000) - 0.0205900
RATES(49) = (ALGBRC(6) - STATES(49))/ALGBRC(23)
ALGBRC(24) = 0.821390*EXP(- (STATES(24)+0.0315900)/0.0274600 ** 2.00000)+0.000190000
RATES(50) = (ALGBRC(5) - STATES(50))/ALGBRC(24)
ALGBRC(25) = 0.821390*EXP(- (STATES(25)+0.0315900)/0.0274600 ** 2.00000)+0.000190000
RATES(51) = (ALGBRC(6) - STATES(51))/ALGBRC(25)
ALGBRC(7) = ((( CONSTS(134)*arbitrary_log(STATES(7)/1.00000, 10) - CONSTS(135)) - CONSTS(133)*(STATES(69)/(CONSTS(132)+STATES(69)))) - CONSTS(131)*(STATES(67) ** 2.00000/(CONSTS(130) ** 2.00000+STATES(67) ** 2.00000)))*CONSTS(171)
ALGBRC(26) = 1.00000/(1.00000+EXP(- (STATES(24) - ALGBRC(7))/0.0308000))
RATES(52) = (ALGBRC(26) - STATES(52))/CONSTS(50)
ALGBRC(8) = ((( CONSTS(134)*arbitrary_log(STATES(15)/1.00000, 10) - CONSTS(135)) - CONSTS(133)*(STATES(69)/(CONSTS(132)+STATES(69)))) - CONSTS(131)*(STATES(67) ** 2.00000/(CONSTS(130) ** 2.00000+STATES(67) ** 2.00000)))*CONSTS(171)
ALGBRC(27) = 1.00000/(1.00000+EXP(- (STATES(25) - ALGBRC(8))/0.0308000))
RATES(53) = (ALGBRC(27) - STATES(53))/CONSTS(50)
RATES(54) = (ALGBRC(26) - STATES(54))/CONSTS(51)
RATES(55) = (ALGBRC(27) - STATES(55))/CONSTS(51)
ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(24)+0.00188000)/0.00757000))
ALGBRC(28) = 0.00289000*EXP(- (STATES(24)+0.00863000)/0.0123900 ** 2.00000)+0.00243000
RATES(60) = (ALGBRC(13) - STATES(60))/ALGBRC(28)
ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(25)+0.00188000)/0.00757000))
ALGBRC(29) = 0.0289000*EXP(- (STATES(25)+0.00863000)/0.0123900 ** 2.00000)+0.00243000
RATES(61) = (ALGBRC(14) - STATES(61))/ALGBRC(29)
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(24)+0.0293200)/0.00154000))
ALGBRC(30) = 0.295590*EXP(- (STATES(24) - 0.00472000)/0.112550 ** 2.00000)+0.0231900
RATES(62) = (ALGBRC(15) - STATES(62))/ALGBRC(30)
ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(25)+0.0293200)/0.00154000))
ALGBRC(31) = 0.295590*EXP(- (STATES(25) - 0.00472000)/0.112550 ** 2.00000)+0.0231900
RATES(63) = (ALGBRC(16) - STATES(63))/ALGBRC(31)
ALGBRC(19) = 1.00000 - (STATES(71)+STATES(72))
ALGBRC(33) = CONSTS(123)*STATES(67)
RATES(71) = - CONSTS(119)*STATES(71)*STATES(69)+ CONSTS(120)*ALGBRC(19)+ ALGBRC(33)*STATES(72)
RATES(72) = ( CONSTS(122)*ALGBRC(19)*STATES(69)+ CONSTS(121)*ALGBRC(19)) - ALGBRC(33)*STATES(72)
ALGBRC(32) = CONSTS(157)*STATES(12)
ALGBRC(36) = 1.00000+ CONSTS(142)*(STATES(68) ** 2.00000/(CONSTS(141) ** 2.00000+STATES(68) ** 2.00000))
ALGBRC(39) = ( CONSTS(158)*(1.00000+ 3.65000*(STATES(67) ** 2.00000/(CONSTS(170) ** 2.00000+STATES(67) ** 2.00000))))/ALGBRC(36)
RATES(64) = - ALGBRC(32)*STATES(64)+ ALGBRC(39)*STATES(65)+ CONSTS(161)*STATES(66)
ALGBRC(17) = CONSTS(164) - (STATES(64)+STATES(65)+STATES(66))
RATES(65) = (( ALGBRC(32)*STATES(64) - ALGBRC(39)*STATES(65)) - CONSTS(159)*STATES(65))+ CONSTS(160)*ALGBRC(17)
ALGBRC(41) = ALGBRC(39)
RATES(66) = ( - CONSTS(161)*STATES(66)+ ALGBRC(41)*ALGBRC(17)) - ALGBRC(32)*STATES(66)
ALGBRC(63) = (STATES(28)/CONSTS(19))*STATES(34)
ALGBRC(61) = (STATES(28)/CONSTS(91))*STATES(32)
ALGBRC(59) = (STATES(28)/CONSTS(19))*STATES(30)
CALL minimize(minfunc_2, CONSTS, VARIABLES, <VAR>)
RATES(32) = CONSTS(20)*( STATES(7)*STATES(30) - CONSTS(21)*STATES(32)) - CONSTS(22)*( STATES(7)*STATES(32) - CONSTS(23)*ALGBRC(65))
RATES(34) = CONSTS(22)*( STATES(7)*STATES(30) - CONSTS(23)*STATES(34)) - CONSTS(20)*( STATES(7)*STATES(34) - CONSTS(21)*ALGBRC(65))
ALGBRC(67) = (STATES(29)/CONSTS(19))*STATES(35)
ALGBRC(62) = (STATES(29)/CONSTS(91))*STATES(33)
ALGBRC(60) = (STATES(29)/CONSTS(19))*STATES(31)
CALL minimize(minfunc_3, CONSTS, VARIABLES, <VAR>)
RATES(33) = CONSTS(20)*( STATES(15)*STATES(31) - CONSTS(21)*STATES(33)) - CONSTS(22)*( STATES(15)*STATES(33) - CONSTS(23)*ALGBRC(69))
RATES(35) = CONSTS(22)*( STATES(15)*STATES(31) - CONSTS(23)*STATES(35)) - CONSTS(20)*( STATES(15)*STATES(35) - CONSTS(21)*ALGBRC(69))
RATES(22) = CONSTS(72)*STATES(20)*(CONSTS(71) - STATES(22)) - CONSTS(73)*STATES(22)
RATES(23) = CONSTS(72)*STATES(21)*(CONSTS(71) - STATES(23)) - CONSTS(73)*STATES(23)
ALGBRC(35) = (CONSTS(179)/CONSTS(87))*log(CONSTS(83)/STATES(1))
ALGBRC(47) = (CONSTS(179)/CONSTS(89))*log(CONSTS(85)/STATES(5))
ALGBRC(78) = CONSTS(28)*CONSTS(74)*((ALGBRC(35) - ALGBRC(47))/((ALGBRC(35) - ALGBRC(47))+CONSTS(26)))
ALGBRC(80) = CONSTS(29)*CONSTS(74)*(( ( CONSTS(84)*CONSTS(83)*CONSTS(85) ** 2.00000 - STATES(3)*STATES(1)*STATES(5) ** 2.00000)*CONSTS(30)*CONSTS(31)*CONSTS(32)*CONSTS(33))/( (CONSTS(84)+CONSTS(30))*(CONSTS(83)+CONSTS(31))*(CONSTS(85)+CONSTS(32))*(CONSTS(85)+CONSTS(33))*(STATES(3)+CONSTS(30))*(STATES(1)+CONSTS(31))*(STATES(5)+CONSTS(32))*(STATES(5)+CONSTS(33))))
ALGBRC(43) = (CONSTS(179)/CONSTS(88))*log(CONSTS(84)/STATES(3))
ALGBRC(76) = CONSTS(27)*CONSTS(74)*(ALGBRC(43) - ALGBRC(47) ** 4.00000/(ALGBRC(43) - ALGBRC(47) ** 4.00000+CONSTS(26) ** 4.00000))
ALGBRC(48) = CONSTS(43)*CONSTS(74)*(STATES(24) - ALGBRC(47))
ALGBRC(82) = CONSTS(39)*CONSTS(74)*(STATES(39)+STATES(40)+STATES(41))*(STATES(24) - ALGBRC(47))
ALGBRC(121) = (( CONSTS(104)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(122) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(186)*(STATES(6) - STATES(5)), CONSTS(186)*ALGBRC(121)*((STATES(6) - STATES(5)*EXP(- ALGBRC(121)))/(1.00000 - EXP(- ALGBRC(121)))))
RATES(5) = (ALGBRC(48)+ALGBRC(82)+( CONSTS(1)*ALGBRC(76)+ CONSTS(1)*(ALGBRC(78)+ 2.00000*ALGBRC(80)+ALGBRC(122))))/( CONSTS(1)*CONSTS(2))
ALGBRC(40) = (CONSTS(179)/CONSTS(87))*log(CONSTS(83)/STATES(2))
ALGBRC(49) = (CONSTS(179)/CONSTS(89))*log(CONSTS(85)/STATES(6))
ALGBRC(79) = CONSTS(28)*CONSTS(75)*((ALGBRC(40) - ALGBRC(49))/((ALGBRC(40) - ALGBRC(49))+CONSTS(26)))
ALGBRC(81) = CONSTS(29)*CONSTS(75)*(( ( CONSTS(84)*CONSTS(83)*CONSTS(85) ** 2.00000 - STATES(4)*STATES(2)*STATES(6) ** 2.00000)*CONSTS(30)*CONSTS(31)*CONSTS(32)*CONSTS(33))/( (CONSTS(84)+CONSTS(30))*(CONSTS(83)+CONSTS(31))*(CONSTS(85)+CONSTS(32))*(CONSTS(85)+CONSTS(33))*(STATES(4)+CONSTS(30))*(STATES(2)+CONSTS(31))*(STATES(6)+CONSTS(32))*(STATES(6)+CONSTS(33))))
ALGBRC(45) = (CONSTS(179)/CONSTS(88))*log(CONSTS(84)/STATES(4))
ALGBRC(77) = CONSTS(27)*CONSTS(75)*(ALGBRC(45) - ALGBRC(49) ** 4.00000/(ALGBRC(45) - ALGBRC(49) ** 4.00000+CONSTS(26) ** 4.00000))
ALGBRC(50) = CONSTS(43)*CONSTS(75)*(STATES(25) - ALGBRC(49))
ALGBRC(83) = CONSTS(39)*CONSTS(75)*(STATES(45)+STATES(46)+STATES(47))*(STATES(25) - ALGBRC(49))
RATES(6) = (ALGBRC(50)+ALGBRC(83)+( CONSTS(1)*ALGBRC(77)+ CONSTS(1)*((ALGBRC(79)+ 2.00000*ALGBRC(81)) - ALGBRC(122))))/( CONSTS(1)*CONSTS(3))
ALGBRC(123) = CONSTS(162) - (STATES(13)+STATES(8)+STATES(10)+STATES(11)+STATES(12))
RATES(8) = (( CONSTS(145)*STATES(7) ** 2.00000*STATES(13) - CONSTS(146)*STATES(8)) - CONSTS(147)*STATES(7) ** 2.00000*STATES(8))+ CONSTS(148)*ALGBRC(123)
ALGBRC(124) = CONSTS(163) - (STATES(11)+STATES(12))
RATES(10) = (((( CONSTS(143)*STATES(7) ** 2.00000*STATES(13) - CONSTS(144)*STATES(10)) - CONSTS(149)*STATES(7) ** 2.00000*STATES(10))+ CONSTS(150)*ALGBRC(123)) - CONSTS(151)*ALGBRC(124)*STATES(10))+ CONSTS(152)*STATES(11)
RATES(12) = (( CONSTS(153)*STATES(7) ** 2.00000*STATES(11) - CONSTS(154)*STATES(12))+ CONSTS(155)*ALGBRC(124)*ALGBRC(123)) - CONSTS(156)*STATES(12)
RATES(11) = (( CONSTS(151)*ALGBRC(124)*STATES(10) - CONSTS(152)*STATES(11)) - CONSTS(153)*STATES(7) ** 2.00000*STATES(11))+ CONSTS(154)*STATES(12)
ALGBRC(129) = STATES(70) ** (1.0 / 2)/( CONSTS(137) ** (1.0 / 2)+ STATES(70) ** (1.0 / 2))
ALGBRC(130) = STATES(68)/(CONSTS(138)+STATES(68))
ALGBRC(131) = CONSTS(167)*(1.00000 - ALGBRC(129))*(1.00000 - ALGBRC(130))
ALGBRC(125) = (CONSTS(165)/CONSTS(172))*(1.00000 - STATES(67)/( 2.00000*(CONSTS(128)+STATES(67))))
ALGBRC(132) = ALGBRC(131)*CONSTS(74)*((STATES(7)/ALGBRC(125) ** CONSTS(7) - STATES(20)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(7)/ALGBRC(125) ** CONSTS(7)+STATES(20)/CONSTS(6) ** CONSTS(7))))
ALGBRC(64) = CONSTS(13)*ALGBRC(63) ** 3.00000*(STATES(20) - STATES(7))*2.00000*CONSTS(1)*CONSTS(4)
ALGBRC(68) = CONSTS(13)*ALGBRC(67) ** 3.00000*(STATES(20) - STATES(15))*2.00000*CONSTS(1)*CONSTS(4)
ALGBRC(133) = ALGBRC(131)*CONSTS(75)*((STATES(15)/ALGBRC(125) ** CONSTS(7) - STATES(20)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(15)/ALGBRC(125) ** CONSTS(7)+STATES(20)/CONSTS(6) ** CONSTS(7))))
RATES(20) = (((ALGBRC(132) - ALGBRC(64))+ALGBRC(133)) - ALGBRC(68))/( CONSTS(1)*CONSTS(4)) - RATES(22)
ALGBRC(51) = (CONSTS(179)/CONSTS(90))*log(CONSTS(86)/STATES(7))
ALGBRC(52) = CONSTS(42)*CONSTS(74)*(STATES(24) - ALGBRC(51))
ALGBRC(126) = CONSTS(166)*(1.00000+( 3.60000*STATES(67))/(CONSTS(129)+STATES(67)))
ALGBRC(127) = ALGBRC(126)*CONSTS(74)*(STATES(7)/(CONSTS(66)+STATES(7)))
ALGBRC(71) = (STATES(20)+STATES(21))/2.00000
ALGBRC(72) = CONSTS(24)*(CONSTS(25)/(CONSTS(25)+ALGBRC(71)))*(STATES(24) - ALGBRC(51))
ALGBRC(111) = 0.740000*STATES(62)+0.260000
ALGBRC(112) = CONSTS(67)*CONSTS(74)*STATES(60)*ALGBRC(111)*(STATES(24) - ALGBRC(51))
ALGBRC(55) = ( STATES(26)*(1.00000+STATES(7)/CONSTS(10) ** 3.00000))/(1.00000+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)
ALGBRC(56) = CONSTS(8)*ALGBRC(55)*(STATES(21) - STATES(7))*( 2.00000*CONSTS(1)*CONSTS(4))
ALGBRC(119) = (( CONSTS(102)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(120) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(185)*(STATES(15) - STATES(7)), CONSTS(185)*ALGBRC(119)*((STATES(15) - STATES(7)*EXP(- ALGBRC(119)))/(1.00000 - EXP(- ALGBRC(119)))))
ALGBRC(134) = ALGBRC(131)*CONSTS(74)*((STATES(7)/ALGBRC(125) ** CONSTS(7) - STATES(21)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(7)/ALGBRC(125) ** CONSTS(7)+STATES(21)/CONSTS(6) ** CONSTS(7))))
RATES(7) = ((- (((ALGBRC(52)+ALGBRC(127)+ALGBRC(72)+ALGBRC(112)+ALGBRC(132)+ALGBRC(134)) - ALGBRC(56)) - ALGBRC(64))+ 2.00000*CONSTS(1)*ALGBRC(120))/( 2.00000*CONSTS(1)*CONSTS(5)) - 2.00000*(RATES(8)+RATES(10)+RATES(11))) - RATES(9)
ALGBRC(135) = ALGBRC(131)*CONSTS(75)*((STATES(15)/ALGBRC(125) ** CONSTS(7) - STATES(21)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(15)/ALGBRC(125) ** CONSTS(7)+STATES(21)/CONSTS(6) ** CONSTS(7))))
ALGBRC(57) = ( STATES(27)*(1.00000+STATES(15)/CONSTS(10) ** 3.00000))/(1.00000+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)
ALGBRC(58) = CONSTS(8)*ALGBRC(57)*(STATES(21) - STATES(15))*( 2.00000*CONSTS(1)*CONSTS(4))
RATES(21) = (((ALGBRC(134) - ALGBRC(56))+ALGBRC(135)) - ALGBRC(58))/( CONSTS(1)*CONSTS(4)) - RATES(23)
ALGBRC(137) = CONSTS(162) - (STATES(73)+STATES(16)+STATES(17)+STATES(18)+STATES(19))
RATES(16) = (( CONSTS(145)*STATES(15) ** 2.00000*STATES(73) - CONSTS(146)*STATES(16)) - CONSTS(147)*STATES(15) ** 2.00000*STATES(16))+ CONSTS(148)*ALGBRC(137)
ALGBRC(140) = CONSTS(163) - (STATES(18)+STATES(19))
RATES(17) = (((( CONSTS(143)*STATES(15) ** 2.00000*STATES(73) - CONSTS(144)*STATES(17)) - CONSTS(149)*STATES(15) ** 2.00000*STATES(17))+ CONSTS(150)*ALGBRC(137)) - CONSTS(151)*ALGBRC(140)*STATES(17))+ CONSTS(152)*STATES(18)
RATES(19) = (( CONSTS(153)*STATES(15) ** 2.00000*STATES(18) - CONSTS(154)*STATES(19))+ CONSTS(155)*ALGBRC(140)*ALGBRC(137)) - CONSTS(156)*STATES(19)
ALGBRC(38) = CONSTS(40)*CONSTS(74)*(STATES(24) - ALGBRC(35))
ALGBRC(84) = 3.24100/(1.00000+EXP(((STATES(24) - ALGBRC(35)) - 185.720)/1.59600))
ALGBRC(86) = 1.00000*(( 0.0130000*EXP(((STATES(24) - ALGBRC(35))+5.79220)/0.294000)+EXP(((STATES(24) - ALGBRC(35)) - 653.733)/0.244000))/(1.00000+EXP(- ((STATES(24) - ALGBRC(35)) - 7.77750)/0.0928000)))
ALGBRC(87) = ALGBRC(84)/(ALGBRC(84)+ALGBRC(86))
ALGBRC(88) = CONSTS(44)*CONSTS(74)*ALGBRC(87)* (CONSTS(83)/CONSTS(45)) ** (1.0 / 2)*(STATES(24) - ALGBRC(35))
ALGBRC(96) = CONSTS(48)*CONSTS(74)*CONSTS(83)/CONSTS(45) ** CONSTS(47)*(STATES(24) - ALGBRC(35))
ALGBRC(92) = 0.580000*STATES(48)+ 0.420000*STATES(50)
ALGBRC(93) = CONSTS(46)*CONSTS(74)*ALGBRC(92) ** 2.00000*(STATES(24) - ALGBRC(35))
ALGBRC(98) = 0.650000*STATES(52)+ 0.350000*STATES(54)
ALGBRC(99) = CONSTS(49)*CONSTS(74)*ALGBRC(98)*(STATES(24) - ALGBRC(35))
ALGBRC(139) = CONSTS(168)*(1.00000 - STATES(70)/(CONSTS(139)+STATES(70)))*(1.00000 - STATES(68)/(CONSTS(140)+STATES(68)))
ALGBRC(102) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(24)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))+ 0.0365000*CONSTS(180)*EXP(( - STATES(24)*CONSTS(1))/( CONSTS(76)*CONSTS(77))))
ALGBRC(141) = ALGBRC(139)*ALGBRC(102)*(CONSTS(83)/(CONSTS(83)+CONSTS(52)))*(STATES(3) ** 1.50000/(STATES(3) ** 1.50000+CONSTS(53) ** 1.50000))
ALGBRC(115) = (( CONSTS(98)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(116) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(187)*(STATES(2) - STATES(1)), CONSTS(187)*ALGBRC(115)*((STATES(2) - STATES(1)*EXP(- ALGBRC(115)))/(1.00000 - EXP(- ALGBRC(115)))))
RATES(1) = (- ((ALGBRC(38)+ALGBRC(88)+ALGBRC(96)+ALGBRC(93)+ALGBRC(99)) - 2.00000*ALGBRC(141))+ CONSTS(1)*(ALGBRC(78)+ALGBRC(80)+ALGBRC(116)))/( CONSTS(1)*CONSTS(2))
ALGBRC(42) = CONSTS(40)*CONSTS(75)*(STATES(25) - ALGBRC(40))
ALGBRC(85) = 3.24100/(1.00000+EXP(((STATES(25) - ALGBRC(40)) - 185.720)/1.59600))
ALGBRC(89) = 1.00000*(( 0.0130000*EXP(((STATES(25) - ALGBRC(40))+5.79220)/0.294000)+EXP(((STATES(25) - ALGBRC(40)) - 653.733)/0.244000))/(1.00000+EXP(- (((STATES(25) - ALGBRC(40)) - 7.77750)/0.0928000))))
ALGBRC(90) = ALGBRC(85)/(ALGBRC(85)+ALGBRC(89))
ALGBRC(91) = CONSTS(44)*CONSTS(75)*ALGBRC(90)* (CONSTS(83)/CONSTS(45)) ** (1.0 / 2)*(STATES(25) - ALGBRC(40))
ALGBRC(97) = CONSTS(48)*CONSTS(75)*CONSTS(83)/CONSTS(45) ** CONSTS(47)*(STATES(25) - ALGBRC(40))
ALGBRC(94) = 0.580000*STATES(49)+ 0.420000*STATES(51)
ALGBRC(95) = CONSTS(46)*CONSTS(75)*ALGBRC(94) ** 2.00000*(STATES(25) - ALGBRC(40))
ALGBRC(100) = 0.650000*STATES(53)+ 0.350000*STATES(55)
ALGBRC(101) = CONSTS(49)*CONSTS(75)*ALGBRC(100)*(STATES(25) - ALGBRC(40))
ALGBRC(103) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))+ 0.0365000*CONSTS(180)*EXP(( - STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77))))
ALGBRC(142) = ALGBRC(139)*ALGBRC(103)*(CONSTS(83)/(CONSTS(83)+CONSTS(52)))*(STATES(4) ** 1.50000/(STATES(4) ** 1.50000+CONSTS(54) ** 1.50000))
RATES(2) = (- ((ALGBRC(42)+ALGBRC(91)+ALGBRC(97)+ALGBRC(95)+ALGBRC(101)) - 2.00000*ALGBRC(142))+ CONSTS(1)*((ALGBRC(79)+ALGBRC(81)) - ALGBRC(116)))/( CONSTS(1)*CONSTS(3))
ALGBRC(44) = CONSTS(41)*CONSTS(74)*(STATES(24) - ALGBRC(43))
ALGBRC(109) = CONSTS(63)*CONSTS(74)*STATES(56)*STATES(58)*(STATES(24) - ALGBRC(43))
ALGBRC(74) = ALGBRC(72)*( (CONSTS(92) ** 2.00000/CONSTS(93) ** 2.00000)*CONSTS(94))*((STATES(3) - CONSTS(84)*EXP(( - CONSTS(92)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))))/(STATES(7) - CONSTS(86)*EXP(( - CONSTS(93)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77)))))*((1.00000 - EXP(( - CONSTS(93)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))))/(1.00000 - EXP(( - CONSTS(92)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77)))))
ALGBRC(117) = (( CONSTS(100)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(118) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(188)*(STATES(4) - STATES(3)), CONSTS(188)*ALGBRC(117)*((STATES(4) - STATES(3)*EXP(- ALGBRC(117)))/(1.00000 - EXP(- ALGBRC(117)))))
RATES(3) = (- (ALGBRC(44)+ 3.00000*ALGBRC(141)+ALGBRC(109)+ALGBRC(74))+ CONSTS(1)*(ALGBRC(76)+ALGBRC(80)+ALGBRC(118)))/( CONSTS(1)*CONSTS(2))
ALGBRC(46) = CONSTS(41)*CONSTS(75)*(STATES(25) - ALGBRC(45))
ALGBRC(110) = CONSTS(63)*CONSTS(75)*STATES(57)*STATES(59)*(STATES(25) - ALGBRC(45))
ALGBRC(53) = (CONSTS(179)/CONSTS(90))*log(CONSTS(86)/STATES(15))
ALGBRC(73) = CONSTS(24)*(CONSTS(25)/(CONSTS(25)+ALGBRC(71)))*(STATES(25) - ALGBRC(53))
ALGBRC(75) = ALGBRC(73)*( (CONSTS(92) ** 2.00000/CONSTS(93) ** 2.00000)*CONSTS(94))*((STATES(4) - CONSTS(84)*EXP(( - CONSTS(92)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77))))/(STATES(15) - CONSTS(86)*EXP(( - CONSTS(93)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77)))))*((1.00000 - EXP(( - CONSTS(93)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77))))/(1.00000 - EXP(( - CONSTS(92)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77)))))
ALGBRC(104) = 1.00000/(1.00000+CONSTS(62)/STATES(15) ** 2.00000)
ALGBRC(105) = EXP(( CONSTS(61)*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))
ALGBRC(106) = EXP(( (CONSTS(61) - 1.00000)*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))
ALGBRC(107) = CONSTS(84) ** 3.00000*STATES(15)+ STATES(4) ** 3.00000*CONSTS(86)+ CONSTS(59) ** 3.00000*STATES(15)+ CONSTS(57)*STATES(4) ** 3.00000+ CONSTS(58) ** 3.00000*CONSTS(86)*(1.00000+STATES(15)/CONSTS(56))+ CONSTS(56)*CONSTS(84) ** 3.00000*(1.00000+STATES(4)/CONSTS(58) ** 3.00000)
ALGBRC(108) = CONSTS(55)*ALGBRC(104)*(( ALGBRC(105)*STATES(4) ** 3.00000*CONSTS(86) - ALGBRC(106)*CONSTS(84) ** 3.00000*STATES(15))/( ALGBRC(107)*(1.00000+ CONSTS(60)*ALGBRC(106))))
RATES(4) = (- (ALGBRC(46)+ 3.00000*ALGBRC(142)+ALGBRC(110)+ALGBRC(75)+ 3.00000*ALGBRC(108))+ CONSTS(1)*((ALGBRC(77)+ALGBRC(81)) - ALGBRC(118)))/( CONSTS(1)*CONSTS(3))
RATES(18) = (( CONSTS(151)*ALGBRC(140)*STATES(17) - CONSTS(152)*STATES(18)) - CONSTS(153)*STATES(15) ** 2.00000*STATES(18))+ CONSTS(154)*STATES(19)
ALGBRC(54) = CONSTS(42)*CONSTS(75)*(STATES(25) - ALGBRC(53))
ALGBRC(128) = ALGBRC(126)*CONSTS(75)*(STATES(15)/(CONSTS(66)+STATES(15)))
ALGBRC(113) = 0.740000*STATES(63)+0.260000
ALGBRC(114) = CONSTS(67)*CONSTS(75)*STATES(61)*ALGBRC(113)*(STATES(25) - ALGBRC(53))
RATES(15) = ((- (((((ALGBRC(54)+ALGBRC(128)+ALGBRC(73)+ALGBRC(114)) - 2.00000*ALGBRC(108))+ALGBRC(135)+ALGBRC(133)) - ALGBRC(58)) - ALGBRC(68)) - 2.00000*CONSTS(1)*ALGBRC(120))/( 2.00000*CONSTS(1)*CONSTS(3)) - 2.00000*(RATES(16)+RATES(17)+RATES(18))) - RATES(14)
ALGBRC(136) = ALGBRC(132)+ALGBRC(134)
ALGBRC(144) = ((((((((ALGBRC(38)+ALGBRC(88)+ALGBRC(96)+ALGBRC(93)+ALGBRC(99)) - CONSTS(1)*ALGBRC(116))+ALGBRC(141)+ALGBRC(44)+ALGBRC(109)+ALGBRC(74)) - CONSTS(1)*ALGBRC(118))+ALGBRC(48)+ALGBRC(82)+ CONSTS(1)*ALGBRC(122)+ALGBRC(52)+ALGBRC(127)+ALGBRC(112)+ALGBRC(72)) - 2.00000*CONSTS(1)*ALGBRC(120))+ALGBRC(136)) - ALGBRC(56)) - ALGBRC(64)
RATES(24) = - ALGBRC(144)/( CONSTS(74)*CONSTS(106))
ALGBRC(138) = ALGBRC(133)+ALGBRC(135)
ALGBRC(145) = ((((ALGBRC(42)+ALGBRC(91)+ALGBRC(97)+ALGBRC(95)+ALGBRC(101)+ CONSTS(1)*ALGBRC(116)+ALGBRC(142)+ALGBRC(46)+ALGBRC(110)+ALGBRC(75)+ALGBRC(108)+ CONSTS(1)*ALGBRC(118)+ALGBRC(50)+ALGBRC(83)) - CONSTS(1)*ALGBRC(122))+ALGBRC(54)+ALGBRC(128)+ALGBRC(114)+ALGBRC(73)+ 2.00000*CONSTS(1)*ALGBRC(120)+ALGBRC(138)) - ALGBRC(58)) - ALGBRC(68)
RATES(25) = - ALGBRC(145)/( CONSTS(75)*CONSTS(106))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = (1.00000+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)/(1.00000+1.00000/CONSTS(11)+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)
ALGBRC(2) = (1.00000+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)/(1.00000+1.00000/CONSTS(11)+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(24)+0.0239000)/0.00480000))
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(25)+0.0239000)/0.00480000))
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(24)+0.0661000)/0.00650000))
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(25)+0.0661000)/0.00650000))
ALGBRC(18) = ( - CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))
ALGBRC(20) = CONSTS(97)*EXP(CONSTS(95)+ CONSTS(96)*STATES(24))
ALGBRC(21) = CONSTS(97)*EXP(CONSTS(95)+ CONSTS(96)*STATES(25))
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(24)+0.00177000)/0.0145200))
ALGBRC(22) = 0.210987*EXP(- (STATES(24)+0.214340)/0.195350 ** 2.00000) - 0.0205900
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(25)+0.00177000)/0.0145200))
ALGBRC(23) = 0.210987*EXP(- (STATES(25)+0.214340)/0.195350 ** 2.00000) - 0.0205900
ALGBRC(24) = 0.821390*EXP(- (STATES(24)+0.0315900)/0.0274600 ** 2.00000)+0.000190000
ALGBRC(25) = 0.821390*EXP(- (STATES(25)+0.0315900)/0.0274600 ** 2.00000)+0.000190000
ALGBRC(7) = ((( CONSTS(134)*arbitrary_log(STATES(7)/1.00000, 10) - CONSTS(135)) - CONSTS(133)*(STATES(69)/(CONSTS(132)+STATES(69)))) - CONSTS(131)*(STATES(67) ** 2.00000/(CONSTS(130) ** 2.00000+STATES(67) ** 2.00000)))*CONSTS(171)
ALGBRC(26) = 1.00000/(1.00000+EXP(- (STATES(24) - ALGBRC(7))/0.0308000))
ALGBRC(8) = ((( CONSTS(134)*arbitrary_log(STATES(15)/1.00000, 10) - CONSTS(135)) - CONSTS(133)*(STATES(69)/(CONSTS(132)+STATES(69)))) - CONSTS(131)*(STATES(67) ** 2.00000/(CONSTS(130) ** 2.00000+STATES(67) ** 2.00000)))*CONSTS(171)
ALGBRC(27) = 1.00000/(1.00000+EXP(- (STATES(25) - ALGBRC(8))/0.0308000))
ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(24)+0.00188000)/0.00757000))
ALGBRC(28) = 0.00289000*EXP(- (STATES(24)+0.00863000)/0.0123900 ** 2.00000)+0.00243000
ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(25)+0.00188000)/0.00757000))
ALGBRC(29) = 0.0289000*EXP(- (STATES(25)+0.00863000)/0.0123900 ** 2.00000)+0.00243000
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(24)+0.0293200)/0.00154000))
ALGBRC(30) = 0.295590*EXP(- (STATES(24) - 0.00472000)/0.112550 ** 2.00000)+0.0231900
ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(25)+0.0293200)/0.00154000))
ALGBRC(31) = 0.295590*EXP(- (STATES(25) - 0.00472000)/0.112550 ** 2.00000)+0.0231900
ALGBRC(19) = 1.00000 - (STATES(71)+STATES(72))
ALGBRC(33) = CONSTS(123)*STATES(67)
ALGBRC(32) = CONSTS(157)*STATES(12)
ALGBRC(36) = 1.00000+ CONSTS(142)*(STATES(68) ** 2.00000/(CONSTS(141) ** 2.00000+STATES(68) ** 2.00000))
ALGBRC(39) = ( CONSTS(158)*(1.00000+ 3.65000*(STATES(67) ** 2.00000/(CONSTS(170) ** 2.00000+STATES(67) ** 2.00000))))/ALGBRC(36)
ALGBRC(17) = CONSTS(164) - (STATES(64)+STATES(65)+STATES(66))
ALGBRC(41) = ALGBRC(39)
ALGBRC(63) = (STATES(28)/CONSTS(19))*STATES(34)
ALGBRC(61) = (STATES(28)/CONSTS(91))*STATES(32)
ALGBRC(59) = (STATES(28)/CONSTS(19))*STATES(30)
ALGBRC(67) = (STATES(29)/CONSTS(19))*STATES(35)
ALGBRC(62) = (STATES(29)/CONSTS(91))*STATES(33)
ALGBRC(60) = (STATES(29)/CONSTS(19))*STATES(31)
ALGBRC(35) = (CONSTS(179)/CONSTS(87))*log(CONSTS(83)/STATES(1))
ALGBRC(47) = (CONSTS(179)/CONSTS(89))*log(CONSTS(85)/STATES(5))
ALGBRC(78) = CONSTS(28)*CONSTS(74)*((ALGBRC(35) - ALGBRC(47))/((ALGBRC(35) - ALGBRC(47))+CONSTS(26)))
ALGBRC(80) = CONSTS(29)*CONSTS(74)*(( ( CONSTS(84)*CONSTS(83)*CONSTS(85) ** 2.00000 - STATES(3)*STATES(1)*STATES(5) ** 2.00000)*CONSTS(30)*CONSTS(31)*CONSTS(32)*CONSTS(33))/( (CONSTS(84)+CONSTS(30))*(CONSTS(83)+CONSTS(31))*(CONSTS(85)+CONSTS(32))*(CONSTS(85)+CONSTS(33))*(STATES(3)+CONSTS(30))*(STATES(1)+CONSTS(31))*(STATES(5)+CONSTS(32))*(STATES(5)+CONSTS(33))))
ALGBRC(43) = (CONSTS(179)/CONSTS(88))*log(CONSTS(84)/STATES(3))
ALGBRC(76) = CONSTS(27)*CONSTS(74)*(ALGBRC(43) - ALGBRC(47) ** 4.00000/(ALGBRC(43) - ALGBRC(47) ** 4.00000+CONSTS(26) ** 4.00000))
ALGBRC(48) = CONSTS(43)*CONSTS(74)*(STATES(24) - ALGBRC(47))
ALGBRC(82) = CONSTS(39)*CONSTS(74)*(STATES(39)+STATES(40)+STATES(41))*(STATES(24) - ALGBRC(47))
ALGBRC(121) = (( CONSTS(104)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(122) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(186)*(STATES(6) - STATES(5)), CONSTS(186)*ALGBRC(121)*((STATES(6) - STATES(5)*EXP(- ALGBRC(121)))/(1.00000 - EXP(- ALGBRC(121)))))
ALGBRC(40) = (CONSTS(179)/CONSTS(87))*log(CONSTS(83)/STATES(2))
ALGBRC(49) = (CONSTS(179)/CONSTS(89))*log(CONSTS(85)/STATES(6))
ALGBRC(79) = CONSTS(28)*CONSTS(75)*((ALGBRC(40) - ALGBRC(49))/((ALGBRC(40) - ALGBRC(49))+CONSTS(26)))
ALGBRC(81) = CONSTS(29)*CONSTS(75)*(( ( CONSTS(84)*CONSTS(83)*CONSTS(85) ** 2.00000 - STATES(4)*STATES(2)*STATES(6) ** 2.00000)*CONSTS(30)*CONSTS(31)*CONSTS(32)*CONSTS(33))/( (CONSTS(84)+CONSTS(30))*(CONSTS(83)+CONSTS(31))*(CONSTS(85)+CONSTS(32))*(CONSTS(85)+CONSTS(33))*(STATES(4)+CONSTS(30))*(STATES(2)+CONSTS(31))*(STATES(6)+CONSTS(32))*(STATES(6)+CONSTS(33))))
ALGBRC(45) = (CONSTS(179)/CONSTS(88))*log(CONSTS(84)/STATES(4))
ALGBRC(77) = CONSTS(27)*CONSTS(75)*(ALGBRC(45) - ALGBRC(49) ** 4.00000/(ALGBRC(45) - ALGBRC(49) ** 4.00000+CONSTS(26) ** 4.00000))
ALGBRC(50) = CONSTS(43)*CONSTS(75)*(STATES(25) - ALGBRC(49))
ALGBRC(83) = CONSTS(39)*CONSTS(75)*(STATES(45)+STATES(46)+STATES(47))*(STATES(25) - ALGBRC(49))
ALGBRC(123) = CONSTS(162) - (STATES(13)+STATES(8)+STATES(10)+STATES(11)+STATES(12))
ALGBRC(124) = CONSTS(163) - (STATES(11)+STATES(12))
ALGBRC(129) = STATES(70) ** (1.0 / 2)/( CONSTS(137) ** (1.0 / 2)+ STATES(70) ** (1.0 / 2))
ALGBRC(130) = STATES(68)/(CONSTS(138)+STATES(68))
ALGBRC(131) = CONSTS(167)*(1.00000 - ALGBRC(129))*(1.00000 - ALGBRC(130))
ALGBRC(125) = (CONSTS(165)/CONSTS(172))*(1.00000 - STATES(67)/( 2.00000*(CONSTS(128)+STATES(67))))
ALGBRC(132) = ALGBRC(131)*CONSTS(74)*((STATES(7)/ALGBRC(125) ** CONSTS(7) - STATES(20)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(7)/ALGBRC(125) ** CONSTS(7)+STATES(20)/CONSTS(6) ** CONSTS(7))))
ALGBRC(64) = CONSTS(13)*ALGBRC(63) ** 3.00000*(STATES(20) - STATES(7))*2.00000*CONSTS(1)*CONSTS(4)
ALGBRC(68) = CONSTS(13)*ALGBRC(67) ** 3.00000*(STATES(20) - STATES(15))*2.00000*CONSTS(1)*CONSTS(4)
ALGBRC(133) = ALGBRC(131)*CONSTS(75)*((STATES(15)/ALGBRC(125) ** CONSTS(7) - STATES(20)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(15)/ALGBRC(125) ** CONSTS(7)+STATES(20)/CONSTS(6) ** CONSTS(7))))
ALGBRC(51) = (CONSTS(179)/CONSTS(90))*log(CONSTS(86)/STATES(7))
ALGBRC(52) = CONSTS(42)*CONSTS(74)*(STATES(24) - ALGBRC(51))
ALGBRC(126) = CONSTS(166)*(1.00000+( 3.60000*STATES(67))/(CONSTS(129)+STATES(67)))
ALGBRC(127) = ALGBRC(126)*CONSTS(74)*(STATES(7)/(CONSTS(66)+STATES(7)))
ALGBRC(71) = (STATES(20)+STATES(21))/2.00000
ALGBRC(72) = CONSTS(24)*(CONSTS(25)/(CONSTS(25)+ALGBRC(71)))*(STATES(24) - ALGBRC(51))
ALGBRC(111) = 0.740000*STATES(62)+0.260000
ALGBRC(112) = CONSTS(67)*CONSTS(74)*STATES(60)*ALGBRC(111)*(STATES(24) - ALGBRC(51))
ALGBRC(55) = ( STATES(26)*(1.00000+STATES(7)/CONSTS(10) ** 3.00000))/(1.00000+CONSTS(9)/STATES(7) ** 4.00000+STATES(7)/CONSTS(10) ** 3.00000)
ALGBRC(56) = CONSTS(8)*ALGBRC(55)*(STATES(21) - STATES(7))*( 2.00000*CONSTS(1)*CONSTS(4))
ALGBRC(119) = (( CONSTS(102)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(120) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(185)*(STATES(15) - STATES(7)), CONSTS(185)*ALGBRC(119)*((STATES(15) - STATES(7)*EXP(- ALGBRC(119)))/(1.00000 - EXP(- ALGBRC(119)))))
ALGBRC(134) = ALGBRC(131)*CONSTS(74)*((STATES(7)/ALGBRC(125) ** CONSTS(7) - STATES(21)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(7)/ALGBRC(125) ** CONSTS(7)+STATES(21)/CONSTS(6) ** CONSTS(7))))
ALGBRC(135) = ALGBRC(131)*CONSTS(75)*((STATES(15)/ALGBRC(125) ** CONSTS(7) - STATES(21)/CONSTS(6) ** CONSTS(7))/(1.00000+(STATES(15)/ALGBRC(125) ** CONSTS(7)+STATES(21)/CONSTS(6) ** CONSTS(7))))
ALGBRC(57) = ( STATES(27)*(1.00000+STATES(15)/CONSTS(10) ** 3.00000))/(1.00000+CONSTS(9)/STATES(15) ** 4.00000+STATES(15)/CONSTS(10) ** 3.00000)
ALGBRC(58) = CONSTS(8)*ALGBRC(57)*(STATES(21) - STATES(15))*( 2.00000*CONSTS(1)*CONSTS(4))
ALGBRC(137) = CONSTS(162) - (STATES(73)+STATES(16)+STATES(17)+STATES(18)+STATES(19))
ALGBRC(140) = CONSTS(163) - (STATES(18)+STATES(19))
ALGBRC(38) = CONSTS(40)*CONSTS(74)*(STATES(24) - ALGBRC(35))
ALGBRC(84) = 3.24100/(1.00000+EXP(((STATES(24) - ALGBRC(35)) - 185.720)/1.59600))
ALGBRC(86) = 1.00000*(( 0.0130000*EXP(((STATES(24) - ALGBRC(35))+5.79220)/0.294000)+EXP(((STATES(24) - ALGBRC(35)) - 653.733)/0.244000))/(1.00000+EXP(- ((STATES(24) - ALGBRC(35)) - 7.77750)/0.0928000)))
ALGBRC(87) = ALGBRC(84)/(ALGBRC(84)+ALGBRC(86))
ALGBRC(88) = CONSTS(44)*CONSTS(74)*ALGBRC(87)* (CONSTS(83)/CONSTS(45)) ** (1.0 / 2)*(STATES(24) - ALGBRC(35))
ALGBRC(96) = CONSTS(48)*CONSTS(74)*CONSTS(83)/CONSTS(45) ** CONSTS(47)*(STATES(24) - ALGBRC(35))
ALGBRC(92) = 0.580000*STATES(48)+ 0.420000*STATES(50)
ALGBRC(93) = CONSTS(46)*CONSTS(74)*ALGBRC(92) ** 2.00000*(STATES(24) - ALGBRC(35))
ALGBRC(98) = 0.650000*STATES(52)+ 0.350000*STATES(54)
ALGBRC(99) = CONSTS(49)*CONSTS(74)*ALGBRC(98)*(STATES(24) - ALGBRC(35))
ALGBRC(139) = CONSTS(168)*(1.00000 - STATES(70)/(CONSTS(139)+STATES(70)))*(1.00000 - STATES(68)/(CONSTS(140)+STATES(68)))
ALGBRC(102) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(24)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))+ 0.0365000*CONSTS(180)*EXP(( - STATES(24)*CONSTS(1))/( CONSTS(76)*CONSTS(77))))
ALGBRC(141) = ALGBRC(139)*ALGBRC(102)*(CONSTS(83)/(CONSTS(83)+CONSTS(52)))*(STATES(3) ** 1.50000/(STATES(3) ** 1.50000+CONSTS(53) ** 1.50000))
ALGBRC(115) = (( CONSTS(98)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(116) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(187)*(STATES(2) - STATES(1)), CONSTS(187)*ALGBRC(115)*((STATES(2) - STATES(1)*EXP(- ALGBRC(115)))/(1.00000 - EXP(- ALGBRC(115)))))
ALGBRC(42) = CONSTS(40)*CONSTS(75)*(STATES(25) - ALGBRC(40))
ALGBRC(85) = 3.24100/(1.00000+EXP(((STATES(25) - ALGBRC(40)) - 185.720)/1.59600))
ALGBRC(89) = 1.00000*(( 0.0130000*EXP(((STATES(25) - ALGBRC(40))+5.79220)/0.294000)+EXP(((STATES(25) - ALGBRC(40)) - 653.733)/0.244000))/(1.00000+EXP(- (((STATES(25) - ALGBRC(40)) - 7.77750)/0.0928000))))
ALGBRC(90) = ALGBRC(85)/(ALGBRC(85)+ALGBRC(89))
ALGBRC(91) = CONSTS(44)*CONSTS(75)*ALGBRC(90)* (CONSTS(83)/CONSTS(45)) ** (1.0 / 2)*(STATES(25) - ALGBRC(40))
ALGBRC(97) = CONSTS(48)*CONSTS(75)*CONSTS(83)/CONSTS(45) ** CONSTS(47)*(STATES(25) - ALGBRC(40))
ALGBRC(94) = 0.580000*STATES(49)+ 0.420000*STATES(51)
ALGBRC(95) = CONSTS(46)*CONSTS(75)*ALGBRC(94) ** 2.00000*(STATES(25) - ALGBRC(40))
ALGBRC(100) = 0.650000*STATES(53)+ 0.350000*STATES(55)
ALGBRC(101) = CONSTS(49)*CONSTS(75)*ALGBRC(100)*(STATES(25) - ALGBRC(40))
ALGBRC(103) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))+ 0.0365000*CONSTS(180)*EXP(( - STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77))))
ALGBRC(142) = ALGBRC(139)*ALGBRC(103)*(CONSTS(83)/(CONSTS(83)+CONSTS(52)))*(STATES(4) ** 1.50000/(STATES(4) ** 1.50000+CONSTS(54) ** 1.50000))
ALGBRC(44) = CONSTS(41)*CONSTS(74)*(STATES(24) - ALGBRC(43))
ALGBRC(109) = CONSTS(63)*CONSTS(74)*STATES(56)*STATES(58)*(STATES(24) - ALGBRC(43))
ALGBRC(74) = ALGBRC(72)*( (CONSTS(92) ** 2.00000/CONSTS(93) ** 2.00000)*CONSTS(94))*((STATES(3) - CONSTS(84)*EXP(( - CONSTS(92)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))))/(STATES(7) - CONSTS(86)*EXP(( - CONSTS(93)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77)))))*((1.00000 - EXP(( - CONSTS(93)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77))))/(1.00000 - EXP(( - CONSTS(92)*CONSTS(1)*STATES(24))/( CONSTS(76)*CONSTS(77)))))
ALGBRC(117) = (( CONSTS(100)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))*(STATES(25) - STATES(24))
ALGBRC(118) = TERNRY(STATES(25) - STATES(24).EQ.0.00000, CONSTS(188)*(STATES(4) - STATES(3)), CONSTS(188)*ALGBRC(117)*((STATES(4) - STATES(3)*EXP(- ALGBRC(117)))/(1.00000 - EXP(- ALGBRC(117)))))
ALGBRC(46) = CONSTS(41)*CONSTS(75)*(STATES(25) - ALGBRC(45))
ALGBRC(110) = CONSTS(63)*CONSTS(75)*STATES(57)*STATES(59)*(STATES(25) - ALGBRC(45))
ALGBRC(53) = (CONSTS(179)/CONSTS(90))*log(CONSTS(86)/STATES(15))
ALGBRC(73) = CONSTS(24)*(CONSTS(25)/(CONSTS(25)+ALGBRC(71)))*(STATES(25) - ALGBRC(53))
ALGBRC(75) = ALGBRC(73)*( (CONSTS(92) ** 2.00000/CONSTS(93) ** 2.00000)*CONSTS(94))*((STATES(4) - CONSTS(84)*EXP(( - CONSTS(92)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77))))/(STATES(15) - CONSTS(86)*EXP(( - CONSTS(93)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77)))))*((1.00000 - EXP(( - CONSTS(93)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77))))/(1.00000 - EXP(( - CONSTS(92)*CONSTS(1)*STATES(25))/( CONSTS(76)*CONSTS(77)))))
ALGBRC(104) = 1.00000/(1.00000+CONSTS(62)/STATES(15) ** 2.00000)
ALGBRC(105) = EXP(( CONSTS(61)*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))
ALGBRC(106) = EXP(( (CONSTS(61) - 1.00000)*STATES(25)*CONSTS(1))/( CONSTS(76)*CONSTS(77)))
ALGBRC(107) = CONSTS(84) ** 3.00000*STATES(15)+ STATES(4) ** 3.00000*CONSTS(86)+ CONSTS(59) ** 3.00000*STATES(15)+ CONSTS(57)*STATES(4) ** 3.00000+ CONSTS(58) ** 3.00000*CONSTS(86)*(1.00000+STATES(15)/CONSTS(56))+ CONSTS(56)*CONSTS(84) ** 3.00000*(1.00000+STATES(4)/CONSTS(58) ** 3.00000)
ALGBRC(108) = CONSTS(55)*ALGBRC(104)*(( ALGBRC(105)*STATES(4) ** 3.00000*CONSTS(86) - ALGBRC(106)*CONSTS(84) ** 3.00000*STATES(15))/( ALGBRC(107)*(1.00000+ CONSTS(60)*ALGBRC(106))))
ALGBRC(54) = CONSTS(42)*CONSTS(75)*(STATES(25) - ALGBRC(53))
ALGBRC(128) = ALGBRC(126)*CONSTS(75)*(STATES(15)/(CONSTS(66)+STATES(15)))
ALGBRC(113) = 0.740000*STATES(63)+0.260000
ALGBRC(114) = CONSTS(67)*CONSTS(75)*STATES(61)*ALGBRC(113)*(STATES(25) - ALGBRC(53))
ALGBRC(136) = ALGBRC(132)+ALGBRC(134)
ALGBRC(144) = ((((((((ALGBRC(38)+ALGBRC(88)+ALGBRC(96)+ALGBRC(93)+ALGBRC(99)) - CONSTS(1)*ALGBRC(116))+ALGBRC(141)+ALGBRC(44)+ALGBRC(109)+ALGBRC(74)) - CONSTS(1)*ALGBRC(118))+ALGBRC(48)+ALGBRC(82)+ CONSTS(1)*ALGBRC(122)+ALGBRC(52)+ALGBRC(127)+ALGBRC(112)+ALGBRC(72)) - 2.00000*CONSTS(1)*ALGBRC(120))+ALGBRC(136)) - ALGBRC(56)) - ALGBRC(64)
ALGBRC(138) = ALGBRC(133)+ALGBRC(135)
ALGBRC(145) = ((((ALGBRC(42)+ALGBRC(91)+ALGBRC(97)+ALGBRC(95)+ALGBRC(101)+ CONSTS(1)*ALGBRC(116)+ALGBRC(142)+ALGBRC(46)+ALGBRC(110)+ALGBRC(75)+ALGBRC(108)+ CONSTS(1)*ALGBRC(118)+ALGBRC(50)+ALGBRC(83)) - CONSTS(1)*ALGBRC(122))+ALGBRC(54)+ALGBRC(128)+ALGBRC(114)+ALGBRC(73)+ 2.00000*CONSTS(1)*ALGBRC(120)+ALGBRC(138)) - ALGBRC(58)) - ALGBRC(68)
ALGBRC(34) = (STATES(65)+ALGBRC(17))/CONSTS(164)
ALGBRC(37) = ((CONSTS(173) - CONSTS(174)*ALGBRC(34))+ CONSTS(175)*ALGBRC(34) ** 2.00000) - CONSTS(176)*ALGBRC(34) ** 3.00000
ALGBRC(143) = ALGBRC(139)
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0 = abs(1.00000 - ALGBRC(3)+STATES(36)+STATES(37)+STATES(38)+STATES(39)+STATES(40)+STATES(41))
RETURN
END
REAL FUNCTION minfunc_1(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_1 = abs(1.00000 - ALGBRC(4)+STATES(42)+STATES(43)+STATES(44)+STATES(45)+STATES(46)+STATES(47))
RETURN
END
REAL FUNCTION minfunc_2(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_2[1] = abs(ALGBRC(66) - (STATES(28)/CONSTS(91))*ALGBRC(65))
minfunc_2[2] = abs(ALGBRC(65) - (1.00000 - (STATES(30)+STATES(32)+STATES(34)+ALGBRC(59)+ALGBRC(61)+ALGBRC(63)+ALGBRC(66))))
RETURN
END
REAL FUNCTION minfunc_3(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_3[1] = abs(ALGBRC(70) - (STATES(29)/CONSTS(91))*ALGBRC(69))
minfunc_3[2] = abs(ALGBRC(69) - (1.00000 - (STATES(31)+STATES(33)+STATES(35)+ALGBRC(60)+ALGBRC(62)+ALGBRC(67)+ALGBRC(70))))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END