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, ) 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, ) 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