C C There are a total of 140 entries in the algebraic variable array. C There are a total of 26 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 CONSTS(1) is NO_on in component strain_control (dimensionless). C CONSTS(2) is pH_on in component strain_control (dimensionless). C CONSTS(3) is SAC_on in component strain_control (dimensionless). C VOI is time in component environment (ms). C CONSTS(4) is V_myo in component cell_geom (mm3). C CONSTS(5) is V_SR in component cell_geom (mm3). C CONSTS(6) is N in component cell_geom (dimensionless). C CONSTS(7) is A_cap in component cell_geom (mm2). C CONSTS(141) is rho in component cell_geom (per_mm). C STATES(1) is V in component membrane (mV). C CONSTS(8) is R in component membrane (millijoule_per_mole_kelvin). C CONSTS(9) is T in component membrane (kelvin). C CONSTS(10) is F in component membrane (coulomb_per_mole). C CONSTS(11) is CmF in component membrane (uF). C CONSTS(150) is Cm in component membrane (uF_per_mm2). C ALGBRC(25) is i_Na in component sodium_current (uA_per_mm2). C ALGBRC(31) is i_t in component Ca_independent_transient_outward_K_current (uA_per_mm2). C ALGBRC(32) is i_ss in component steady_state_outward_K_current (uA_per_mm2). C ALGBRC(37) is i_f in component hyperpolarisation_activated_current (uA_per_mm2). C ALGBRC(33) is i_K1 in component inward_rectifier (uA_per_mm2). C ALGBRC(42) is i_NaK in component sodium_potassium_pump (uA_per_mm2). C ALGBRC(10) is i_Stim in component I_stimulus (uA_per_mm2). C ALGBRC(128) is I_CaL in component I_Ca_L (uA_per_mm2). C ALGBRC(88) is I_NaCa in component NCX (uA_per_mm2). C ALGBRC(84) is I_Cab in component Cab (uA_per_mm2). C ALGBRC(82) is I_pCa in component SL_pump (uA_per_mm2). C ALGBRC(38) is i_B_Na in component background_currents (uA_per_mm2). C ALGBRC(39) is i_B_K in component background_currents (uA_per_mm2). C ALGBRC(23) is I_SAC_Na in component SAC_current (uA_per_mm2). C ALGBRC(28) is I_SAC_K in component SAC_current (uA_per_mm2). C ALGBRC(30) is I_Ko in component KSA_current (uA_per_mm2). C ALGBRC(74) is I_Cl in component intracellular_ion_concentrations (uA_per_mm2). C ALGBRC(77) is I_H in component intracellular_ion_concentrations (uA_per_mm2). C CONSTS(12) is stim_period in component I_stimulus (ms). C CONSTS(13) is stim_duration in component I_stimulus (ms). C CONSTS(14) is stim_amplitude in component I_stimulus (uA). C ALGBRC(21) is E_Na in component sodium_current (mV). C ALGBRC(26) is E_K in component Ca_independent_transient_outward_K_current (mV). C CONSTS(15) is ExtensionRatio in component Myofilaments (dimensionless). C CONSTS(16) is g_SAC in component SAC_current (mS). C ALGBRC(29) is I_SAC in component SAC_current (uA_per_mm2). C CONSTS(142) is gamma_SLSAC in component SAC_current (dimensionless). C CONSTS(17) is E_R in component SAC_current (mV). C CONSTS(143) is r in component SAC_current (dimensionless). C CONSTS(18) is g_Ko in component KSA_current (mS). C CONSTS(144) is gamma_SLKO in component KSA_current (dimensionless). C CONSTS(19) is g_Na in component sodium_current (mS). C CONSTS(145) is g_Na_endo in component sodium_current (mS). C STATES(2) is Na_i in component intracellular_ion_concentrations (mM). C CONSTS(20) is Na_o in component standard_ionic_concentrations (mM). C STATES(3) is m in component sodium_current_m_gate (dimensionless). C STATES(4) is h in component sodium_current_h_gate (dimensionless). C STATES(5) is j in component sodium_current_j_gate (dimensionless). C ALGBRC(1) is m_infinity in component sodium_current_m_gate (dimensionless). C ALGBRC(13) is tau_m in component sodium_current_m_gate (ms). C ALGBRC(2) is h_infinity in component sodium_current_h_gate (dimensionless). C ALGBRC(14) is tau_h in component sodium_current_h_gate (ms). C ALGBRC(3) is j_infinity in component sodium_current_j_gate (dimensionless). C ALGBRC(15) is tau_j in component sodium_current_j_gate (ms). C CONSTS(21) is g_t in component Ca_independent_transient_outward_K_current (mS). C CONSTS(146) is g_t_endo in component Ca_independent_transient_outward_K_current (mS). C CONSTS(22) is a_endo in component Ca_independent_transient_outward_K_current (dimensionless). C CONSTS(23) is b_endo in component Ca_independent_transient_outward_K_current (dimensionless). C CONSTS(24) is K_o in component standard_ionic_concentrations (mM). C STATES(6) is K_i in component intracellular_ion_concentrations (mM). C STATES(7) is r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless). C STATES(8) is s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless). C STATES(9) is s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless). C ALGBRC(16) is tau_r in component Ca_independent_transient_outward_K_current_r_gate (ms). C ALGBRC(4) is r_infinity in component Ca_independent_transient_outward_K_current_r_gate (dimensionless). C ALGBRC(17) is tau_s_endo in component Ca_independent_transient_outward_K_current_s_gate (ms). C ALGBRC(5) is s_infinity in component Ca_independent_transient_outward_K_current_s_gate (dimensionless). C ALGBRC(18) is tau_s_slow_endo in component Ca_independent_transient_outward_K_current_s_slow_gate (ms). C ALGBRC(6) is s_slow_infinity in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless). C CONSTS(25) is g_ss in component steady_state_outward_K_current (mS). C STATES(10) is r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless). C STATES(11) is s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless). C ALGBRC(19) is tau_r_ss in component steady_state_outward_K_current_r_ss_gate (ms). C ALGBRC(7) is r_ss_infinity in component steady_state_outward_K_current_r_ss_gate (dimensionless). C CONSTS(147) is tau_s_ss in component steady_state_outward_K_current_s_ss_gate (ms). C ALGBRC(8) is s_ss_infinity in component steady_state_outward_K_current_s_ss_gate (dimensionless). C CONSTS(26) is g_K1 in component inward_rectifier (mS). C ALGBRC(34) is i_f_Na in component hyperpolarisation_activated_current (uA_per_mm2). C ALGBRC(35) is i_f_K in component hyperpolarisation_activated_current (uA_per_mm2). C CONSTS(27) is g_f in component hyperpolarisation_activated_current (mS). C CONSTS(28) is f_Na in component hyperpolarisation_activated_current (dimensionless). C CONSTS(148) is f_K in component hyperpolarisation_activated_current (dimensionless). C STATES(12) is y in component hyperpolarisation_activated_current_y_gate (dimensionless). C ALGBRC(20) is tau_y in component hyperpolarisation_activated_current_y_gate (ms). C ALGBRC(9) is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless). C CONSTS(29) is g_B_Na in component background_currents (mS). C CONSTS(30) is g_B_K in component background_currents (mS). C CONSTS(31) is scale_Na in component background_currents (dimensionless). C CONSTS(32) is scale_K in component background_currents (dimensionless). C CONSTS(33) is i_NaK_max in component sodium_potassium_pump (uA). C CONSTS(34) is K_m_K in component sodium_potassium_pump (mM). C CONSTS(35) is K_m_Na in component sodium_potassium_pump (mM). C CONSTS(149) is sigma in component sodium_potassium_pump (dimensionless). C ALGBRC(40) is p_nai in component sodium_potassium_pump (dimensionless). C ALGBRC(41) is p_v in component sodium_potassium_pump (dimensionless). C ALGBRC(43) is jco2 in component J_CO2 (mM_per_ms). C CONSTS(36) is Pco2 in component J_CO2 (mm_per_ms). C STATES(13) is CO2i in component intracellular_ion_concentrations (mM). C CONSTS(183) is CO2e in component intracellular_ion_concentrations (mM). C STATES(14) is pH_i in component intracellular_ion_concentrations (dimensionless). C CONSTS(37) is pH_e in component intracellular_ion_concentrations (dimensionless). C ALGBRC(52) is v_nhe in component comp_v_nhe_exchanger (mM_per_ms). C CONSTS(38) is KA in component comp_v_nhe_exchanger (mM). C CONSTS(39) is KB in component comp_v_nhe_exchanger (mM). C CONSTS(40) is kp1 in component comp_v_nhe_exchanger (per_ms). C CONSTS(41) is km1 in component comp_v_nhe_exchanger (per_ms). C CONSTS(42) is kp2 in component comp_v_nhe_exchanger (per_ms). C CONSTS(164) is km2 in component comp_v_nhe_exchanger (dimensionless). C CONSTS(43) is K_Hi in component comp_v_nhe_exchanger (mM). C CONSTS(44) is n_Hi in component comp_v_nhe_exchanger (dimensionless). C CONSTS(167) is Be in component comp_v_nhe_exchanger (mM). C ALGBRC(47) is Bi in component comp_v_nhe_exchanger (mM). C ALGBRC(48) is am1 in component comp_v_nhe_exchanger (per_ms). C CONSTS(170) is ap1 in component comp_v_nhe_exchanger (per_ms). C CONSTS(172) is am2 in component comp_v_nhe_exchanger (per_ms). C ALGBRC(49) is ap2 in component comp_v_nhe_exchanger (per_ms). C ALGBRC(50) is reg in component comp_v_nhe_exchanger (dimensionless). C ALGBRC(51) is flux_nhe in component comp_v_nhe_exchanger (dimensionless). C CONSTS(45) is Q_10Scaler in component comp_v_nhe_exchanger (dimensionless). C CONSTS(157) is gamma_NHE in component comp_v_nhe_exchanger (dimensionless). C CONSTS(161) is K_Hs in component comp_v_nhe_exchanger (mM). C CONSTS(46) is Cle in component intracellular_ion_concentrations (mM). C STATES(15) is Cli in component intracellular_ion_concentrations (mM). C ALGBRC(58) is v_che in component comp_v_che_exchanger (mM_per_ms). C CONSTS(47) is K_Cl in component comp_v_che_exchanger (mM). C CONSTS(48) is K_OH in component comp_v_che_exchanger (mM). C CONSTS(49) is kp1 in component comp_v_che_exchanger (per_min). C CONSTS(50) is km1 in component comp_v_che_exchanger (per_min). C CONSTS(51) is kp2 in component comp_v_che_exchanger (per_min). C CONSTS(173) is km2 in component comp_v_che_exchanger (per_min). C CONSTS(174) is OHe in component comp_v_che_exchanger (mM). C ALGBRC(53) is OHi in component comp_v_che_exchanger (mM). C CONSTS(175) is a in component comp_v_che_exchanger (dimensionless). C CONSTS(176) is b in component comp_v_che_exchanger (dimensionless). C ALGBRC(54) is c in component comp_v_che_exchanger (dimensionless). C ALGBRC(55) is d in component comp_v_che_exchanger (dimensionless). C ALGBRC(56) is s1 in component comp_v_che_exchanger (dimensionless). C ALGBRC(57) is s6 in component comp_v_che_exchanger (dimensionless). C CONSTS(52) is Q_10Scaler in component comp_v_che_exchanger (dimensionless). C CONSTS(184) is HCO3e in component intracellular_ion_concentrations (mM). C STATES(16) is HCO3i in component intracellular_ion_concentrations (mM). C ALGBRC(65) is v_nbc in component comp_v_nbc (mM_per_ms). C CONSTS(53) is K_Na in component comp_v_nbc (mM). C CONSTS(54) is K_HCO3 in component comp_v_nbc (mM). C CONSTS(55) is kp1 in component comp_v_nbc (per_min). C CONSTS(56) is km1 in component comp_v_nbc (per_min). C CONSTS(57) is kp2 in component comp_v_nbc (per_min). C CONSTS(177) is km2 in component comp_v_nbc (per_min). C CONSTS(58) is K_Hi in component comp_v_nbc (mM). C CONSTS(59) is n_Hi in component comp_v_nbc (dimensionless). C CONSTS(60) is K_He in component comp_v_nbc (mM). C CONSTS(61) is n_He in component comp_v_nbc (dimensionless). C CONSTS(178) is He in component comp_v_nbc (mM). C ALGBRC(61) is Hi in component comp_v_nbc (mM). C CONSTS(185) is a in component comp_v_nbc (dimensionless). C CONSTS(186) is b in component comp_v_nbc (dimensionless). C ALGBRC(59) is c in component comp_v_nbc (dimensionless). C ALGBRC(60) is d in component comp_v_nbc (dimensionless). C ALGBRC(62) is s1 in component comp_v_nbc (dimensionless). C ALGBRC(63) is s6 in component comp_v_nbc (dimensionless). C ALGBRC(64) is reg in component comp_v_nbc (dimensionless). C CONSTS(62) is Q_10Scaler in component comp_v_nbc (dimensionless). C ALGBRC(72) is v_ae in component comp_v_ae (mM_per_ms). C CONSTS(63) is K_Cl in component comp_v_ae (mM). C CONSTS(64) is K_HCO3 in component comp_v_ae (mM). C CONSTS(65) is kp1 in component comp_v_ae (per_min). C CONSTS(66) is km1 in component comp_v_ae (per_min). C CONSTS(67) is kp2 in component comp_v_ae (per_min). C CONSTS(181) is km2 in component comp_v_ae (per_min). C CONSTS(68) is K_Hi in component comp_v_ae (mM). C CONSTS(69) is n_Hi in component comp_v_ae (dimensionless). C CONSTS(70) is K_He in component comp_v_ae (mM). C CONSTS(71) is n_He in component comp_v_ae (dimensionless). C CONSTS(182) is He in component comp_v_ae (mM). C ALGBRC(66) is Hi in component comp_v_ae (mM). C CONSTS(187) is a in component comp_v_ae (dimensionless). C CONSTS(188) is b in component comp_v_ae (dimensionless). C ALGBRC(67) is c in component comp_v_ae (dimensionless). C ALGBRC(68) is d in component comp_v_ae (dimensionless). C ALGBRC(69) is s1 in component comp_v_ae (dimensionless). C ALGBRC(70) is s6 in component comp_v_ae (dimensionless). C ALGBRC(71) is reg in component comp_v_ae (dimensionless). C CONSTS(72) is Q_10Scaler in component comp_v_ae (dimensionless). C CONSTS(179) is gamma_AE in component comp_v_ae (dimensionless). C CONSTS(180) is K_Hs in component comp_v_ae (mM). C CONSTS(73) is PP_co2e in component intracellular_ion_concentrations (dimensionless). C CONSTS(74) is CO_2sol in component intracellular_ion_concentrations (dimensionless). C CONSTS(75) is P_atm in component intracellular_ion_concentrations (dimensionless). C CONSTS(76) is kf_co2hyd in component intracellular_ion_concentrations (per_ms). C CONSTS(77) is kr_co2hyd in component intracellular_ion_concentrations (per_mM_per_ms). C ALGBRC(80) is v_co2hyd in component intracellular_ion_concentrations (mM_per_ms). C CONSTS(78) is pKa_ib1 in component intracellular_ion_concentrations (dimensionless). C CONSTS(79) is ib1 in component intracellular_ion_concentrations (mM). C CONSTS(80) is pKa_ib2 in component intracellular_ion_concentrations (dimensionless). C CONSTS(81) is ib2 in component intracellular_ion_concentrations (mM). C ALGBRC(78) is beta_intr in component intracellular_ion_concentrations (dimensionless). C CONSTS(82) is pH_scale in component intracellular_ion_concentrations (dimensionless). C ALGBRC(44) is nai_Nak in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(79) is nai_NHE in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(81) is nai_NBC in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(89) is nai_flux in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(36) is nai_bg in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(45) is K_Nak in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(46) is K_flux in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(73) is E_Cl in component intracellular_ion_concentrations (mV). C CONSTS(83) is g_Cl in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(76) is E_H in component intracellular_ion_concentrations (mV). C CONSTS(84) is g_H in component intracellular_ion_concentrations (mM_per_ms). C CONSTS(151) is H_o in component intracellular_ion_concentrations (mM). C ALGBRC(75) is H_i in component intracellular_ion_concentrations (mM). C CONSTS(85) is J_CO2 in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(90) is nai_NaCa in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(27) is nai_na in component intracellular_ion_concentrations (mM_per_ms). C ALGBRC(92) is nai_total in component intracellular_ion_concentrations (mM_per_ms). C CONSTS(86) is Ca_o in component standard_ionic_concentrations (mM). C CONSTS(87) is g_pCa in component SL_pump (mM_per_ms). C CONSTS(88) is K_mpCa in component SL_pump (mM). C STATES(17) is Ca_i in component ionic_concentrations (mM). C CONSTS(89) is delta in component Ca_voltage (per_mV). C ALGBRC(83) is E_Ca in component Cab (mV). C CONSTS(90) is g_Cab in component Cab (mM_per_ms). C ALGBRC(126) is J_LC in component L_flux (mM_per_ms_per_mm3). C CONSTS(91) is K_mNa in component NCX (mM). C CONSTS(92) is K_mCa in component NCX (mM). C CONSTS(93) is eta in component NCX (dimensionless). C CONSTS(94) is k_sat in component NCX (dimensionless). C CONSTS(95) is g_NCX in component NCX (mM_per_ms). C ALGBRC(85) is edv in component NCX (dimensionless). C ALGBRC(86) is edv2 in component NCX (dimensionless). C ALGBRC(87) is Nai3 in component NCX (mM3). C CONSTS(152) is Nae3 in component NCX (mM3). C CONSTS(96) is g_SERCA in component SERCA (mM_per_ms). C CONSTS(97) is K_SERCA in component SERCA (mM). C ALGBRC(91) is I_SERCA in component SERCA (mM_per_ms). C STATES(18) is Ca_SR in component ionic_concentrations (mM). C STATES(19) is TRPN in component ionic_concentrations (mM). C CONSTS(98) is B_CMDN in component ionic_concentrations (mM). C CONSTS(99) is K_CMDN in component ionic_concentrations (mM). C CONSTS(100) is g_SRl in component ionic_concentrations (per_ms). C ALGBRC(140) is J_TPRN in component troponin (dimensionless). C ALGBRC(130) is J_RY in component R_flux (mM_per_ms_per_mm3). C ALGBRC(132) is J_SR in component ionic_concentrations (mM_per_ms). C ALGBRC(11) is Ca_b in component ionic_concentrations (mM). C CONSTS(101) is B_TRPN in component troponin (mM). C CONSTS(102) is g_D in component Ca_conductances (mm3_per_ms). C CONSTS(103) is J_R in component Ca_conductances (mm3_per_ms). C CONSTS(104) is J_L in component Ca_conductances (mm3_per_ms). C ALGBRC(94) is expmdV in component Ca_voltage (dimensionless). C ALGBRC(95) is expVL in component Ca_voltage (dimensionless). C ALGBRC(93) is dV in component Ca_voltage (mV). C CONSTS(105) is V_L0 in component Ca_voltage (mV). C CONSTS(106) is delta_VL in component Ca_voltage (mV). C ALGBRC(12) is C_cc in component C_ij (mM). C ALGBRC(98) is C_oc in component C_ij (mM). C ALGBRC(97) is C_co in component C_ij (mM). C ALGBRC(96) is C_oo in component C_ij (mM). C ALGBRC(102) is J_Loo in component J_ij (mM_per_ms). C ALGBRC(101) is J_Loc in component J_ij (mM_per_ms). C ALGBRC(99) is J_Rco in component J_ij (mM_per_ms). C ALGBRC(100) is J_Roo in component J_ij (mM_per_ms). C CONSTS(153) is t_R in component Ca_tau (ms). C CONSTS(107) is t_L in component Ca_tau (ms). C ALGBRC(105) is epsilon_m in component epsilon (per_ms). C ALGBRC(103) is epsilon_pco in component epsilon (per_ms). C ALGBRC(104) is epsilon_pcc in component epsilon (per_ms). C CONSTS(108) is a in component epsilon (dimensionless). C CONSTS(109) is b in component epsilon (dimensionless). C CONSTS(110) is tau_L in component epsilon (ms). C CONSTS(111) is K_L in component epsilon (mM). C ALGBRC(106) is alpha_p in component alpha (per_ms). C CONSTS(158) is alpha_m in component alpha (per_ms). C CONSTS(112) is phi_L in component alpha (dimensionless). C CONSTS(165) is phi_R in component RyR_param (dimensionless). C CONSTS(113) is phi_R_base in component RyR_param (dimensionless). C CONSTS(114) is tau_R in component RyR_param (ms). C CONSTS(115) is theta_R in component RyR_param (dimensionless). C CONSTS(116) is K_RyR in component RyR_param (mM). C CONSTS(162) is gamma_NO in component RyR_param (dimensionless). C ALGBRC(107) is beta_poc in component beta (per_ms). C ALGBRC(108) is beta_pcc in component beta (per_ms). C CONSTS(168) is beta_m in component beta (per_ms). C ALGBRC(109) is mu_poc in component mu_ij (per_ms). C ALGBRC(110) is mu_pcc in component mu_ij (per_ms). C ALGBRC(111) is mu_moc in component mu_ij (per_ms). C ALGBRC(112) is mu_mcc in component mu_ij (per_ms). C CONSTS(117) is c in component mu_ij (dimensionless). C CONSTS(118) is d in component mu_ij (dimensionless). C ALGBRC(114) is y_oc in component y_ij (dimensionless). C ALGBRC(115) is y_co in component y_ij (dimensionless). C ALGBRC(116) is y_oo in component y_ij (dimensionless). C ALGBRC(117) is y_cc in component y_ij (dimensionless). C ALGBRC(113) is denom in component y_ij (per_ms3). C ALGBRC(119) is r_1 in component r_i (per_ms). C ALGBRC(121) is r_2 in component r_i (per_ms). C ALGBRC(123) is r_3 in component r_i (per_ms). C ALGBRC(125) is r_4 in component r_i (per_ms). C ALGBRC(127) is r_5 in component r_i (per_ms). C ALGBRC(129) is r_6 in component r_i (per_ms). C ALGBRC(131) is r_7 in component r_i (per_ms). C ALGBRC(133) is r_8 in component r_i (per_ms). C STATES(20) is z_1 in component z_i (dimensionless). C STATES(21) is z_2 in component z_i (dimensionless). C STATES(22) is z_3 in component z_i (dimensionless). C ALGBRC(134) is z_4 in component z_i (dimensionless). C ALGBRC(118) is J_R1 in component J_values (mM_per_ms). C ALGBRC(120) is J_R3 in component J_values (mM_per_ms). C ALGBRC(122) is J_L1 in component J_values (mM_per_ms). C ALGBRC(124) is J_L2 in component J_values (mM_per_ms). C ALGBRC(139) is betaCab in component troponin (per_ms). C CONSTS(119) is k_on in component troponin (per_mM_per_ms). C CONSTS(120) is k_off in component troponin (per_ms). C CONSTS(121) is gamma_trpn in component troponin (dimensionless). C CONSTS(122) is TRPN_tot in component troponin (mM). C ALGBRC(138) is Tension in component Cross_Bridges (N_per_mm2). C CONSTS(123) is T_ref in component length_independent_tension (N_per_mm2). C CONSTS(154) is lamda in component Myofilaments (dimensionless). C CONSTS(155) is dExtensionRatiodt in component Myofilaments (per_ms). C CONSTS(156) is lambda_prev in component Myofilaments (dimensionless). C STATES(23) is z in component tropomyosin (dimensionless). C CONSTS(171) is z_max in component tropomyosin (dimensionless). C CONSTS(124) is alpha_0 in component tropomyosin (per_ms). C CONSTS(125) is alpha_r1 in component tropomyosin (per_ms). C CONSTS(126) is alpha_r2 in component tropomyosin (per_ms). C CONSTS(127) is n_Rel in component tropomyosin (dimensionless). C CONSTS(128) is K_z in component tropomyosin (dimensionless). C CONSTS(129) is n_Hill in component tropomyosin (dimensionless). C CONSTS(130) is Ca_50ref in component tropomyosin (mM). C CONSTS(131) is z_p in component tropomyosin (dimensionless). C CONSTS(132) is beta_1 in component tropomyosin (dimensionless). C CONSTS(166) is Ca_50 in component tropomyosin (mM). C CONSTS(169) is Ca_TRPN_50 in component tropomyosin (mM). C CONSTS(159) is K_2 in component tropomyosin (dimensionless). C CONSTS(163) is K_1 in component tropomyosin (dimensionless). C ALGBRC(22) is alpha_Tm in component tropomyosin (per_ms). C ALGBRC(24) is beta_Tm in component tropomyosin (per_ms). C CONSTS(133) is beta_0 in component filament_overlap (dimensionless). C CONSTS(160) is overlap in component filament_overlap (dimensionless). C ALGBRC(135) is T_Base in component length_independent_tension (N_per_mm2). C ALGBRC(136) is T_0 in component isometric_tension (N_per_mm2). C ALGBRC(137) is Q in component Cross_Bridges (dimensionless). C CONSTS(134) is a in component Cross_Bridges (dimensionless). C STATES(24) is Q_1 in component Cross_Bridges (dimensionless). C STATES(25) is Q_2 in component Cross_Bridges (dimensionless). C STATES(26) is Q_3 in component Cross_Bridges (dimensionless). C CONSTS(135) is A_1 in component Cross_Bridges (dimensionless). C CONSTS(136) is A_2 in component Cross_Bridges (dimensionless). C CONSTS(137) is A_3 in component Cross_Bridges (dimensionless). C CONSTS(138) is alpha_1 in component Cross_Bridges (dimensionless). C CONSTS(139) is alpha_2 in component Cross_Bridges (dimensionless). C CONSTS(140) is alpha_3 in component Cross_Bridges (dimensionless). C RATES(1) is d/dt V in component membrane (mV). C RATES(3) is d/dt m in component sodium_current_m_gate (dimensionless). C RATES(4) is d/dt h in component sodium_current_h_gate (dimensionless). C RATES(5) is d/dt j in component sodium_current_j_gate (dimensionless). C RATES(7) is d/dt r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless). C RATES(8) is d/dt s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless). C RATES(9) is d/dt s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless). C RATES(10) is d/dt r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless). C RATES(11) is d/dt s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless). C RATES(12) is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless). C RATES(16) is d/dt HCO3i in component intracellular_ion_concentrations (mM). C RATES(13) is d/dt CO2i in component intracellular_ion_concentrations (mM). C RATES(14) is d/dt pH_i in component intracellular_ion_concentrations (dimensionless). C RATES(2) is d/dt Na_i in component intracellular_ion_concentrations (mM). C RATES(6) is d/dt K_i in component intracellular_ion_concentrations (mM). C RATES(15) is d/dt Cli in component intracellular_ion_concentrations (mM). C RATES(18) is d/dt Ca_SR in component ionic_concentrations (mM). C RATES(17) is d/dt Ca_i in component ionic_concentrations (mM). C RATES(19) is d/dt TRPN in component ionic_concentrations (mM). C RATES(20) is d/dt z_1 in component z_i (dimensionless). C RATES(21) is d/dt z_2 in component z_i (dimensionless). C RATES(22) is d/dt z_3 in component z_i (dimensionless). C RATES(23) is d/dt z in component tropomyosin (dimensionless). C RATES(24) is d/dt Q_1 in component Cross_Bridges (dimensionless). C RATES(25) is d/dt Q_2 in component Cross_Bridges (dimensionless). C RATES(26) is d/dt Q_3 in component Cross_Bridges (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 0 CONSTS(2) = 0 CONSTS(3) = 0 CONSTS(4) = 2.584e-5 CONSTS(5) = 2.098e-6 CONSTS(6) = 75000 CONSTS(7) = 0.01534 STATES(1) = -80.5561112771341 CONSTS(8) = 8314.5 CONSTS(9) = 295 CONSTS(10) = 96487 CONSTS(11) = 0.0001 CONSTS(12) = 1000 CONSTS(13) = 10 CONSTS(14) = -0.0006 CONSTS(15) = 1 CONSTS(16) = 4.1333e-7 CONSTS(17) = -10 CONSTS(18) = 1.2e-6 CONSTS(19) = 0.0008 STATES(2) = 10.9735342589175 CONSTS(20) = 140 STATES(3) = 0.00419283833975832 STATES(4) = 0.675647510010602 STATES(5) = 0.675294627571238 CONSTS(21) = 3.5e-5 CONSTS(22) = 0.583 CONSTS(23) = 0.417 CONSTS(24) = 5.4 STATES(6) = 142.021261491871 STATES(7) = 0.00218107055088884 STATES(8) = 0.922573773425153 STATES(9) = 0.530207467628341 CONSTS(25) = 7e-6 STATES(10) = 0.00289374682707489 STATES(11) = 0.308572839431787 CONSTS(26) = 2.4e-5 CONSTS(27) = 1.45e-6 CONSTS(28) = 0.2 STATES(12) = 0.00331726547133295 CONSTS(29) = 8.015e-8 CONSTS(30) = 1.38e-7 CONSTS(31) = 0 CONSTS(32) = 0 CONSTS(33) = 9.5e-5 CONSTS(34) = 1.5 CONSTS(35) = 10 CONSTS(36) = 5.8e-5 STATES(13) = 1.23610845162584 STATES(14) = 7.20905044774657 CONSTS(37) = 7.4 CONSTS(38) = 21.4935205515214 CONSTS(39) = 1.77948608115949e-7 CONSTS(40) = 0.51156085723193 CONSTS(41) = 0.14069285629105 CONSTS(42) = 0.001698 CONSTS(43) = 0.000416869383470335 CONSTS(44) = 2 CONSTS(45) = 0.35 CONSTS(46) = 126 STATES(15) = 16.2285651746901 CONSTS(47) = 17970.36 CONSTS(48) = 0.0008907 CONSTS(49) = 14975.6 CONSTS(50) = 257.4 CONSTS(51) = 4084.9 CONSTS(52) = 0.35 STATES(16) = 15.1428688811927 CONSTS(53) = 4866.11 CONSTS(54) = 0.008017 CONSTS(55) = 4677 CONSTS(56) = 777.92 CONSTS(57) = 6.656 CONSTS(58) = 0.000183 CONSTS(59) = 2.91 CONSTS(60) = 6.53e-5 CONSTS(61) = 2.18 CONSTS(62) = 0.35 CONSTS(63) = 983.5 CONSTS(64) = 110.64 CONSTS(65) = 20789.9 CONSTS(66) = 21118.2 CONSTS(67) = 21256.54 CONSTS(68) = 2.67e-5 CONSTS(69) = 5.11 CONSTS(70) = 0.000312 CONSTS(71) = 1.44 CONSTS(72) = 0.35 CONSTS(73) = 0.05 CONSTS(74) = 0.03253 CONSTS(75) = 760 CONSTS(76) = 0.000365 CONSTS(77) = 0.481 CONSTS(78) = 6.40013149605198 CONSTS(79) = 31.2044223705215 CONSTS(80) = 7.48048792987277 CONSTS(81) = 6.84513535697271 CONSTS(82) = 1 CONSTS(83) = 2e-5 CONSTS(84) = 4e-6 CONSTS(85) = 0 CONSTS(86) = 2 CONSTS(87) = 3.5e-6 CONSTS(88) = 0.0005 STATES(17) = 8.60587059887361e-5 CONSTS(89) = 0.075 CONSTS(90) = 2.4216e-8 CONSTS(91) = 87.5 CONSTS(92) = 1.38 CONSTS(93) = 0.35 CONSTS(94) = 0.1 CONSTS(95) = 0.0385 CONSTS(96) = 0.00045 CONSTS(97) = 0.0005 STATES(18) = 0.843733859307907 STATES(19) = 0.0670845085340179 CONSTS(98) = 0.05 CONSTS(99) = 0.00238 CONSTS(100) = 5.2e-6 CONSTS(101) = 0.07 CONSTS(102) = 6.5e-11 CONSTS(103) = 2e-11 CONSTS(104) = 9.13e-13 CONSTS(105) = -2 CONSTS(106) = 7 CONSTS(107) = 1 CONSTS(108) = 0.0625 CONSTS(109) = 14 CONSTS(110) = 650 CONSTS(111) = 0.00022 CONSTS(112) = 2.35 CONSTS(113) = 0.05 CONSTS(114) = 2.43 CONSTS(115) = 0.012 CONSTS(116) = 0.041 CONSTS(117) = 0.01 CONSTS(118) = 100 STATES(20) = 0.988955429362334 STATES(21) = 0.00860686643357905 STATES(22) = 0.00241692704700704 CONSTS(119) = 100 CONSTS(120) = 0.2 CONSTS(121) = 2 CONSTS(122) = 0.07 CONSTS(123) = 56.2 STATES(23) = 0.0175189833295748 CONSTS(124) = 0.008 CONSTS(125) = 0.002 CONSTS(126) = 0.00175 CONSTS(127) = 3 CONSTS(128) = 0.15 CONSTS(129) = 3 CONSTS(130) = 0.00105 CONSTS(131) = 0.85 CONSTS(132) = -4 CONSTS(133) = 4.9 CONSTS(134) = 0.35 STATES(24) = 0 STATES(25) = 0 STATES(26) = 0 CONSTS(135) = -29 CONSTS(136) = 138 CONSTS(137) = 129 CONSTS(138) = 0.03 CONSTS(139) = 0.13 CONSTS(140) = 0.625 CONSTS(141) = CONSTS(7)/CONSTS(4) CONSTS(142) = TERNRY(CONSTS(3).EQ.1.00000, (CONSTS(15) - 1.00000)*10.0000, 0.00000) CONSTS(143) = - (CONSTS(17)+85.0000)/(CONSTS(17) - 65.0000) CONSTS(144) = TERNRY(CONSTS(3).EQ.1.00000, 0.700000+ (CONSTS(15) - 1.00000)*3.00000, 0.700000) CONSTS(145) = 1.33000*CONSTS(19) CONSTS(146) = 0.464700*CONSTS(21) CONSTS(147) = 2100.00 CONSTS(148) = 1.00000 - CONSTS(28) CONSTS(149) = (EXP(CONSTS(20)/67.3000) - 1.00000)/7.00000 CONSTS(150) = CONSTS(11)/CONSTS(7) CONSTS(151) = 1000.00*10.0000 ** - CONSTS(37) CONSTS(152) = CONSTS(20) ** 3.00000 CONSTS(153) = 1.17000*CONSTS(107) CONSTS(154) = TERNRY(CONSTS(15).GT.0.800000.AND.CONSTS(15).LE.1.15000, CONSTS(15), TERNRY(CONSTS(15).GT.1.15000, 1.15000, 0.800000) CONSTS(155) = 0.00000 CONSTS(156) = CONSTS(15) CONSTS(157) = TERNRY(CONSTS(2).EQ.1.00000, 1.00000 - 2.80400*(CONSTS(15) - 1.00000), 1.00000) CONSTS(158) = CONSTS(112)/CONSTS(107) CONSTS(159) = (( CONSTS(126)*CONSTS(131) ** CONSTS(127))/(CONSTS(131) ** CONSTS(127)+CONSTS(128) ** CONSTS(127)))*(1.00000 - ( CONSTS(127)*CONSTS(128) ** CONSTS(127))/(CONSTS(131) ** CONSTS(127)+CONSTS(128) ** CONSTS(127))) CONSTS(160) = 1.00000+ CONSTS(133)*(CONSTS(154) - 1.00000) CONSTS(161) = CONSTS(43)*CONSTS(157) CONSTS(162) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000+ 22.4100*(CONSTS(15) - 1.00000), 1.00000) CONSTS(163) = ( CONSTS(126)*CONSTS(131) ** CONSTS(127) - 1.00000*CONSTS(127)*CONSTS(128) ** CONSTS(127))/CONSTS(131) ** CONSTS(127)+CONSTS(128) ** CONSTS(127) ** 2.00000 CONSTS(164) = ( CONSTS(40)*CONSTS(42))/CONSTS(41) CONSTS(165) = CONSTS(113)*CONSTS(162) CONSTS(166) = CONSTS(130)*(1.00000+ CONSTS(132)*(CONSTS(154) - 1.00000)) CONSTS(167) = 10.0000 ** - CONSTS(37)*1000.00 CONSTS(168) = CONSTS(165)/CONSTS(153) CONSTS(169) = ( CONSTS(166)*CONSTS(122))/(CONSTS(166)+ (CONSTS(120)/CONSTS(119))*(1.00000 - ( (1.00000+ CONSTS(133)*(CONSTS(154) - 1.00000))*0.500000)/CONSTS(121))) CONSTS(170) = ( CONSTS(39)*CONSTS(20)*CONSTS(40))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*CONSTS(20)+ CONSTS(20)*CONSTS(167)+ CONSTS(38)*CONSTS(167)) CONSTS(171) = (CONSTS(124)/CONSTS(169)/CONSTS(122) ** CONSTS(129) - CONSTS(159))/(CONSTS(125)+CONSTS(163)+CONSTS(124)/CONSTS(169)/CONSTS(122) ** CONSTS(129)) CONSTS(172) = ( CONSTS(38)*CONSTS(167)*CONSTS(164))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*CONSTS(20)+ CONSTS(20)*CONSTS(167)+ CONSTS(38)*CONSTS(167)) CONSTS(173) = ( CONSTS(51)*CONSTS(50))/CONSTS(49) CONSTS(174) = 1000.00*10.0000 ** - 14.0000+CONSTS(37) CONSTS(175) = 1.00000+CONSTS(48)/CONSTS(174)+( CONSTS(48)*CONSTS(46))/( CONSTS(174)*CONSTS(47)) CONSTS(176) = 1.00000+CONSTS(47)/CONSTS(46)+( CONSTS(47)*CONSTS(174))/( CONSTS(46)*CONSTS(48)) CONSTS(177) = ( CONSTS(57)*CONSTS(56))/CONSTS(55) CONSTS(178) = 1000.00*10.0000 ** - CONSTS(37) CONSTS(179) = TERNRY(CONSTS(2).EQ.1.00000, 1.00000+ 2.50000*(CONSTS(15) - 1.00000), 1.00000) CONSTS(180) = CONSTS(68)*CONSTS(179) CONSTS(181) = ( CONSTS(67)*CONSTS(66))/CONSTS(65) CONSTS(182) = 1000.00*10.0000 ** - CONSTS(37) CONSTS(183) = CONSTS(73)*CONSTS(74)*CONSTS(75) CONSTS(184) = ( (CONSTS(76)/CONSTS(77))*CONSTS(183))/( 10.0000 ** - CONSTS(37)*1000.00) CONSTS(185) = 1.00000+CONSTS(20)/CONSTS(53)+( CONSTS(20)*CONSTS(184))/( CONSTS(53)*CONSTS(54)) CONSTS(186) = 1.00000+CONSTS(54)/CONSTS(184)+( CONSTS(53)*CONSTS(54))/( CONSTS(184)*CONSTS(20)) CONSTS(187) = 1.00000+CONSTS(64)/CONSTS(184)+( CONSTS(64)*CONSTS(46))/( CONSTS(184)*CONSTS(63)) CONSTS(188) = 1.00000+CONSTS(63)/CONSTS(46)+( CONSTS(63)*CONSTS(184))/( CONSTS(46)*CONSTS(64)) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(24) = CONSTS(135)*CONSTS(155) - CONSTS(138)*STATES(24) RATES(25) = CONSTS(136)*CONSTS(155) - CONSTS(139)*STATES(25) RATES(26) = CONSTS(137)*CONSTS(155) - CONSTS(140)*STATES(26) ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000)) RATES(11) = (ALGBRC(8) - STATES(11))/CONSTS(147) ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+45.0000)/- 6.50000)) ALGBRC(13) = 1.36000/(( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))+ 0.0800000*EXP(- STATES(1)/11.0000)) RATES(3) = (ALGBRC(1) - STATES(3))/ALGBRC(13) ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000)) ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 0.453700*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 3.49000/( 0.135000*EXP(- (STATES(1)+80.0000)/6.80000)+ 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))) RATES(4) = (ALGBRC(2) - STATES(4))/ALGBRC(14) ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000)) ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 11.6300*(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))/EXP( - 2.53500e-07*STATES(1)), 3.49000/( ((STATES(1)+37.7800)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))*( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))+( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))) RATES(5) = (ALGBRC(3) - STATES(5))/ALGBRC(15) ALGBRC(16) = 1000.00/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000))) ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200)) RATES(7) = (ALGBRC(4) - STATES(7))/ALGBRC(16) ALGBRC(17) = 550.000*EXP(- (STATES(1)+70.0000)/25.0000 ** 2.00000)+49.0000 ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410)) RATES(8) = (ALGBRC(5) - STATES(8))/ALGBRC(17) ALGBRC(18) = 3300.00*EXP(- (STATES(1)+70.0000)/30.0000 ** 2.00000)+49.0000 ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410)) RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(18) ALGBRC(19) = 10000.0/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000))) ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+11.5000)/- 11.8200)) RATES(10) = (ALGBRC(7) - STATES(10))/ALGBRC(19) ALGBRC(20) = 1000.00/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900)) ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800)) RATES(12) = (ALGBRC(9) - STATES(12))/ALGBRC(20) ALGBRC(11) = CONSTS(101) - STATES(19) ALGBRC(22) = CONSTS(124)*ALGBRC(11)/CONSTS(169) ** CONSTS(129) ALGBRC(24) = CONSTS(125)+( CONSTS(126)*STATES(23) ** CONSTS(127) - 1.00000)/(STATES(23) ** CONSTS(127)+CONSTS(128) ** CONSTS(127)) RATES(23) = ALGBRC(22)*(1.00000 - STATES(23)) - ALGBRC(24)*STATES(23) ALGBRC(26) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(24)/STATES(6)) ALGBRC(31) = ( CONSTS(146)*STATES(7)*( CONSTS(22)*STATES(8)+ CONSTS(23)*STATES(9))*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(32) = ( CONSTS(25)*STATES(10)*STATES(11)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(33) = (( (0.0480000/(EXP((STATES(1)+37.0000)/25.0000)+EXP((STATES(1)+37.0000)/- 25.0000))+0.0100000)*0.00100000)/(1.00000+EXP((STATES(1) - (ALGBRC(26)+76.7700))/- 17.0000))+( CONSTS(26)*(STATES(1) - (ALGBRC(26)+1.73000)))/( (1.00000+EXP(( 1.61300*CONSTS(10)*(STATES(1) - (ALGBRC(26)+1.73000)))/( CONSTS(8)*CONSTS(9))))*(1.00000+EXP((CONSTS(24) - 0.998800)/- 0.124000))))/CONSTS(7) ALGBRC(40) = 1.00000/(1.00000+CONSTS(35)/STATES(2) ** 4.00000) ALGBRC(41) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))+ 0.0365000*CONSTS(149)*EXP(( - STATES(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))) ALGBRC(42) = (( (CONSTS(33)/CONSTS(7))*ALGBRC(41)*CONSTS(24))/(CONSTS(24)+CONSTS(34)))*ALGBRC(40) ALGBRC(10) = TERNRY(VOI - INT(VOI/CONSTS(12))*CONSTS(12).GE.0.00000.AND.VOI - INT(VOI/CONSTS(12))*CONSTS(12).LE.CONSTS(13), CONSTS(14)/CONSTS(7), 0.00000) ALGBRC(39) = ( CONSTS(32)*CONSTS(30)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(28) = ( CONSTS(16)*CONSTS(142)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(30) = ( (CONSTS(18)/(1.00000+EXP(- (10.0000+STATES(1))/45.0000)))*(STATES(1) - ALGBRC(26))*CONSTS(144))/CONSTS(7) ALGBRC(35) = ( CONSTS(27)*STATES(12)*CONSTS(148)*(STATES(1) - ALGBRC(26)))/CONSTS(7) RATES(6) = ( - (((ALGBRC(10)+ALGBRC(28)+ALGBRC(32)+ALGBRC(30)+ALGBRC(31)+ALGBRC(33)+ALGBRC(35)) - 2.00000*ALGBRC(42))+ALGBRC(39))*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(73) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(STATES(15)/CONSTS(46)) ALGBRC(74) = CONSTS(83)*(STATES(1) - ALGBRC(73)) ALGBRC(53) = 1000.00*10.0000 ** - 14.0000+STATES(14) ALGBRC(54) = 1.00000+CONSTS(47)/STATES(15)+( CONSTS(47)*ALGBRC(53))/( STATES(15)*CONSTS(48)) ALGBRC(55) = 1.00000+CONSTS(48)/ALGBRC(53)+( CONSTS(48)*STATES(15))/( ALGBRC(53)*CONSTS(47)) ALGBRC(56) = 1.00000/(CONSTS(175)+( ALGBRC(55)*(CONSTS(49)+( CONSTS(51)*CONSTS(175))/CONSTS(176)))/(CONSTS(50)+( CONSTS(173)*ALGBRC(55))/ALGBRC(54))) ALGBRC(57) = 1.00000/(ALGBRC(55)+( CONSTS(175)*(CONSTS(50)+( CONSTS(173)*ALGBRC(55))/ALGBRC(54)))/(CONSTS(49)+( CONSTS(51)*CONSTS(175))/CONSTS(176))) ALGBRC(58) = (( CONSTS(4)*( CONSTS(50)*ALGBRC(57) - CONSTS(49)*ALGBRC(56)))/( 60.0000*1000.00))*CONSTS(52) ALGBRC(67) = 1.00000+CONSTS(63)/STATES(15)+( CONSTS(63)*STATES(16))/( STATES(15)*CONSTS(64)) ALGBRC(68) = 1.00000+CONSTS(64)/STATES(16)+( CONSTS(64)*STATES(15))/( STATES(16)*CONSTS(63)) ALGBRC(69) = 1.00000/(CONSTS(187)+( ALGBRC(68)*(CONSTS(65)+( CONSTS(67)*CONSTS(187))/CONSTS(188)))/(CONSTS(66)+( CONSTS(181)*ALGBRC(68))/ALGBRC(67))) ALGBRC(70) = 1.00000/(ALGBRC(68)+( CONSTS(187)*(CONSTS(66)+( CONSTS(181)*ALGBRC(68))/ALGBRC(67)))/(CONSTS(65)+( CONSTS(67)*CONSTS(187))/CONSTS(188))) ALGBRC(66) = 1000.00*10.0000 ** - STATES(14) ALGBRC(71) = ( (CONSTS(180) ** CONSTS(69)/(ALGBRC(66) ** CONSTS(69)+CONSTS(180) ** CONSTS(69)))*CONSTS(182) ** CONSTS(71))/(CONSTS(182) ** CONSTS(71)+CONSTS(70) ** CONSTS(71)) ALGBRC(72) = (( CONSTS(4)*ALGBRC(71)*( CONSTS(66)*ALGBRC(70) - CONSTS(65)*ALGBRC(69)))/( 60.0000*1000.00))*CONSTS(72) RATES(15) = ( ALGBRC(74)*CONSTS(7))/( CONSTS(4)*CONSTS(10))+( CONSTS(82)*(ALGBRC(58)+ALGBRC(72)))/CONSTS(4) ALGBRC(59) = 1.00000+CONSTS(54)/STATES(16)+( CONSTS(53)*CONSTS(54))/( STATES(16)*STATES(2)) ALGBRC(60) = 1.00000+STATES(2)/CONSTS(53)+( STATES(2)*STATES(16))/( CONSTS(53)*CONSTS(54)) ALGBRC(62) = 1.00000/(CONSTS(185)+( ALGBRC(60)*(CONSTS(55)+( CONSTS(57)*CONSTS(185))/CONSTS(186)))/(CONSTS(56)+( CONSTS(177)*ALGBRC(60))/ALGBRC(59))) ALGBRC(63) = 1.00000/(ALGBRC(60)+( CONSTS(185)*(CONSTS(56)+( CONSTS(177)*ALGBRC(60))/ALGBRC(59)))/(CONSTS(55)+( CONSTS(57)*CONSTS(185))/CONSTS(186))) ALGBRC(61) = 1000.00*10.0000 ** - STATES(14) ALGBRC(64) = (ALGBRC(61) ** CONSTS(59)/(ALGBRC(61) ** CONSTS(59)+CONSTS(58) ** CONSTS(59)))*(1.00000 - CONSTS(178) ** CONSTS(61)/(CONSTS(178) ** CONSTS(61)+CONSTS(60) ** CONSTS(61))) ALGBRC(65) = (( CONSTS(4)*ALGBRC(64)*( CONSTS(56)*ALGBRC(63) - CONSTS(55)*ALGBRC(62)))/( 60.0000*1000.00))*CONSTS(62) ALGBRC(80) = (CONSTS(4)+CONSTS(5))*( CONSTS(76)*STATES(13) - CONSTS(77)*STATES(16)*10.0000 ** - STATES(14)*1000.00) RATES(16) = ALGBRC(80)/(CONSTS(4)+CONSTS(5))+(ALGBRC(65) - ALGBRC(72))/CONSTS(4) ALGBRC(43) = CONSTS(4)*CONSTS(141)*CONSTS(36)*(CONSTS(183) - STATES(13)) RATES(13) = (ALGBRC(43)/CONSTS(4) - ALGBRC(80)/(CONSTS(4)+CONSTS(5)))+CONSTS(85) ALGBRC(75) = 1000.00*10.0000 ** - STATES(14) ALGBRC(76) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(151)/ALGBRC(75)) ALGBRC(77) = CONSTS(84)*(STATES(1) - ALGBRC(76)) ALGBRC(47) = 10.0000 ** - STATES(14)*1000.00 ALGBRC(48) = ( CONSTS(39)*STATES(2)*CONSTS(41))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*STATES(2)+ STATES(2)*ALGBRC(47)+ CONSTS(38)*ALGBRC(47)) ALGBRC(49) = ( CONSTS(38)*ALGBRC(47)*CONSTS(42))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*STATES(2)+ STATES(2)*ALGBRC(47)+ CONSTS(38)*ALGBRC(47)) ALGBRC(50) = ALGBRC(47) ** CONSTS(44)/(ALGBRC(47) ** CONSTS(44)+CONSTS(161) ** CONSTS(44)) ALGBRC(51) = ( ALGBRC(50)*( CONSTS(170)*ALGBRC(49) - ALGBRC(48)*CONSTS(172)))/(CONSTS(170)+ALGBRC(49)+ALGBRC(48)+CONSTS(172)) ALGBRC(52) = ALGBRC(51)*CONSTS(4)*CONSTS(45) ALGBRC(78) = log(10.0000)*(10.0000 ** - STATES(14)+( 10.0000 ** STATES(14)+CONSTS(78)*CONSTS(79))/10.0000 ** STATES(14)+10.0000 ** CONSTS(78) ** 2.00000+( 10.0000 ** STATES(14)+CONSTS(80)*CONSTS(81))/10.0000 ** STATES(14)+10.0000 ** CONSTS(80) ** 2.00000) RATES(14) = (1.00000/- ALGBRC(78))*(((- ALGBRC(52)+ALGBRC(58))/CONSTS(4)+ALGBRC(80)/(CONSTS(4)+CONSTS(5))) - ( ALGBRC(77)*CONSTS(7))/( CONSTS(4)*CONSTS(10))) ALGBRC(21) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(20)/STATES(2)) ALGBRC(25) = ( CONSTS(145)*STATES(3) ** 3.00000*STATES(4)*STATES(5)*(STATES(1) - ALGBRC(21)))/CONSTS(7) ALGBRC(85) = EXP( CONSTS(89)*0.500000*STATES(1)*CONSTS(93)) ALGBRC(86) = EXP( CONSTS(89)*0.500000*STATES(1)*(CONSTS(93) - 1.00000)) ALGBRC(87) = STATES(2) ** 3.00000 ALGBRC(88) = ( (( (CONSTS(95)/( (CONSTS(152)+CONSTS(91) ** 3.00000)*(CONSTS(86)+CONSTS(92))))*( ALGBRC(85)*ALGBRC(87)*CONSTS(86) - ALGBRC(86)*CONSTS(152)*STATES(17)))/(1.00000+ CONSTS(94)*ALGBRC(86)))*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(38) = ( CONSTS(31)*CONSTS(29)*(STATES(1) - ALGBRC(21)))/CONSTS(7) ALGBRC(23) = (( CONSTS(16)*CONSTS(142)*(STATES(1) - ALGBRC(21)))/CONSTS(7))*CONSTS(143) ALGBRC(34) = ( CONSTS(27)*STATES(12)*CONSTS(28)*(STATES(1) - ALGBRC(21)))/CONSTS(7) RATES(2) = ( - (ALGBRC(38)+ALGBRC(23)+ALGBRC(25)+ ALGBRC(88)*3.00000+ ALGBRC(42)*3.00000+ALGBRC(34))*CONSTS(7))/( CONSTS(4)*CONSTS(10))+(ALGBRC(52)+ALGBRC(65))/CONSTS(4) ALGBRC(37) = ALGBRC(34)+ALGBRC(35) ALGBRC(93) = CONSTS(89)*STATES(1) ALGBRC(94) = EXP(- ALGBRC(93)) ALGBRC(102) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( (( CONSTS(104)*ALGBRC(93))/(1.00000 - ALGBRC(94)))*(( CONSTS(86)*ALGBRC(94) - STATES(17))+ (CONSTS(103)/CONSTS(102))*( CONSTS(86)*ALGBRC(94) - STATES(18))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( (( CONSTS(104)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(86)*EXP(- 1.00000e-05) - STATES(17))+ (CONSTS(103)/CONSTS(102))*( CONSTS(86)*EXP(- 1.00000e-05) - STATES(18))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(101) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( (( CONSTS(104)*ALGBRC(93))/(1.00000 - ALGBRC(94)))*( CONSTS(86)*ALGBRC(94) - STATES(17)))/(1.00000+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( (( CONSTS(104)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(86)*EXP(- 1.00000e-05) - STATES(17)))/(1.00000+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(95) = EXP((STATES(1) - CONSTS(105))/CONSTS(106)) ALGBRC(106) = ALGBRC(95)/( CONSTS(107)*(ALGBRC(95)+1.00000)) ALGBRC(108) = ( (1.00000/CONSTS(153))*STATES(17) ** 2.00000)/(STATES(17) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(98) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-09, (STATES(17)+( (CONSTS(104)/CONSTS(102))*CONSTS(86)*ALGBRC(93)*ALGBRC(94))/(1.00000 - ALGBRC(94)))/(1.00000+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), (STATES(17)+ (CONSTS(104)/CONSTS(102))*CONSTS(86))/(1.00000+CONSTS(104)/CONSTS(102))) ALGBRC(107) = ( (1.00000/CONSTS(153))*ALGBRC(98) ** 2.00000)/(ALGBRC(98) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(113) = (ALGBRC(106)+CONSTS(158))*( (CONSTS(168)+ALGBRC(107)+CONSTS(158))*(CONSTS(168)+ALGBRC(108))+ ALGBRC(106)*(CONSTS(168)+ALGBRC(107))) ALGBRC(114) = ( ALGBRC(106)*CONSTS(168)*(ALGBRC(106)+CONSTS(158)+CONSTS(168)+ALGBRC(108)))/ALGBRC(113) ALGBRC(116) = ( ALGBRC(106)*( ALGBRC(107)*(ALGBRC(106)+CONSTS(168)+ALGBRC(108))+ ALGBRC(108)*CONSTS(158)))/ALGBRC(113) ALGBRC(122) = ALGBRC(102)*ALGBRC(116)+ ALGBRC(101)*ALGBRC(114) ALGBRC(124) = ( ALGBRC(101)*ALGBRC(106))/(ALGBRC(106)+CONSTS(158)) ALGBRC(126) = ( ( STATES(20)*ALGBRC(122)+ STATES(21)*ALGBRC(124))*CONSTS(6))/CONSTS(4) ALGBRC(128) = ( - ALGBRC(126)*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(83) = arbitrary_log(CONSTS(86)/STATES(17), 10)/CONSTS(89) ALGBRC(84) = ( CONSTS(90)*(STATES(1) - ALGBRC(83))*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(82) = ( (( CONSTS(87)*STATES(17))/(CONSTS(88)+STATES(17)))*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) RATES(1) = - (ALGBRC(77)+ALGBRC(74)+ALGBRC(25)+ALGBRC(128)+ALGBRC(31)+ALGBRC(32)+ALGBRC(37)+ALGBRC(33)+ALGBRC(84)+ALGBRC(23)+ALGBRC(39)+ALGBRC(38)+ALGBRC(30)+ALGBRC(28)+ALGBRC(42)+ALGBRC(88)+ALGBRC(82)+ALGBRC(10))/CONSTS(150) ALGBRC(109) = ( (1.00000/CONSTS(114))*(ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/(ALGBRC(98) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(110) = ( (1.00000/CONSTS(114))*(STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/(STATES(17) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(117) = ( CONSTS(158)*CONSTS(168)*(CONSTS(158)+ALGBRC(106)+CONSTS(168)+ALGBRC(107)))/ALGBRC(113) ALGBRC(119) = ALGBRC(114)*ALGBRC(109)+ ALGBRC(117)*ALGBRC(110) ALGBRC(111) = ( (CONSTS(115)/CONSTS(114))*CONSTS(118)*(ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/( CONSTS(118)*ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000) ALGBRC(112) = ( (CONSTS(115)/CONSTS(114))*CONSTS(118)*(STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/( CONSTS(118)*STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000) ALGBRC(121) = ( ALGBRC(106)*ALGBRC(111)+ CONSTS(158)*ALGBRC(112))/(ALGBRC(106)+CONSTS(158)) ALGBRC(97) = (STATES(17)+ (CONSTS(103)/CONSTS(102))*STATES(18))/(1.00000+CONSTS(103)/CONSTS(102)) ALGBRC(103) = ( (( (1.00000/CONSTS(110))*ALGBRC(97))/CONSTS(111))*(ALGBRC(95)+CONSTS(108)))/(ALGBRC(95)+1.00000) ALGBRC(104) = ( (( (1.00000/CONSTS(110))*STATES(17))/CONSTS(111))*(ALGBRC(95)+CONSTS(108)))/(ALGBRC(95)+1.00000) ALGBRC(115) = ( CONSTS(158)*( ALGBRC(108)*(CONSTS(158)+CONSTS(168)+ALGBRC(107))+ ALGBRC(107)*ALGBRC(106)))/ALGBRC(113) ALGBRC(127) = ALGBRC(115)*ALGBRC(103)+ ALGBRC(117)*ALGBRC(104) ALGBRC(105) = ( (1.00000/CONSTS(110))*CONSTS(109)*(ALGBRC(95)+CONSTS(108)))/( CONSTS(109)*ALGBRC(95)+CONSTS(108)) ALGBRC(129) = ALGBRC(105) RATES(20) = - (ALGBRC(119)+ALGBRC(127))*STATES(20)+ ALGBRC(121)*STATES(21)+ ALGBRC(129)*STATES(22) ALGBRC(91) = ( CONSTS(96)*STATES(17) ** 2.00000)/(CONSTS(97) ** 2.00000+STATES(17) ** 2.00000) ALGBRC(99) = ( CONSTS(103)*(STATES(18) - STATES(17)))/(1.00000+CONSTS(103)/CONSTS(102)) ALGBRC(100) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( CONSTS(103)*((STATES(18) - STATES(17))+ (( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94)))*(STATES(18) - CONSTS(86)*ALGBRC(94))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( CONSTS(103)*((STATES(18) - STATES(17))+ (( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(18) - CONSTS(86)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(118) = ALGBRC(116)*ALGBRC(100)+ ALGBRC(99)*ALGBRC(115) ALGBRC(120) = ( ALGBRC(99)*ALGBRC(108))/(CONSTS(168)+ALGBRC(108)) ALGBRC(130) = ( ( STATES(20)*ALGBRC(118)+ STATES(22)*ALGBRC(120))*CONSTS(6))/CONSTS(4) ALGBRC(132) = (- ALGBRC(130)+ALGBRC(91)) - CONSTS(100)*(STATES(18) - STATES(17)) RATES(18) = (CONSTS(4)/CONSTS(5))*ALGBRC(132) ALGBRC(131) = ( CONSTS(158)*ALGBRC(104))/(ALGBRC(106)+CONSTS(158)) ALGBRC(133) = ALGBRC(105) ALGBRC(134) = ((1.00000 - STATES(20)) - STATES(21)) - STATES(22) RATES(21) = ( ALGBRC(119)*STATES(20) - (ALGBRC(121)+ALGBRC(131))*STATES(21))+ ALGBRC(133)*ALGBRC(134) ALGBRC(123) = ( CONSTS(168)*ALGBRC(110))/(CONSTS(168)+ALGBRC(108)) ALGBRC(125) = ALGBRC(112) RATES(22) = ( ALGBRC(127)*STATES(20) - (ALGBRC(129)+ALGBRC(123))*STATES(22))+ ALGBRC(125)*ALGBRC(134) ALGBRC(135) = ( CONSTS(123)*STATES(23))/CONSTS(171) ALGBRC(136) = ALGBRC(135)*CONSTS(160) ALGBRC(137) = STATES(24)+STATES(25)+STATES(26) ALGBRC(138) = TERNRY(ALGBRC(137).LT.0.00000, ( ALGBRC(136)*( CONSTS(134)*ALGBRC(137)+1.00000))/(1.00000 - ALGBRC(137)), ( ALGBRC(136)*(1.00000+ (CONSTS(134)+2.00000)*ALGBRC(137)))/(1.00000+ALGBRC(137))) ALGBRC(139) = TERNRY(1.00000 - ALGBRC(138)/( CONSTS(121)*CONSTS(123)).GT.0.100000, CONSTS(120)*(1.00000 - ALGBRC(138)/( CONSTS(121)*CONSTS(123))), CONSTS(120)*0.100000) ALGBRC(140) = (CONSTS(101) - STATES(19))*ALGBRC(139) - STATES(17)*STATES(19)*CONSTS(119) RATES(17) = (1.00000/(1.00000+( CONSTS(98)*CONSTS(99))/( (STATES(17)+CONSTS(99))*(STATES(17)+CONSTS(99)))))*((ALGBRC(140) - ALGBRC(132))+( ((( 2.00000*ALGBRC(88) - ALGBRC(82)) - ALGBRC(84)) - ALGBRC(128))*CONSTS(7))/( 2.00000*CONSTS(4)*CONSTS(10))) RATES(19) = ALGBRC(140) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000)) ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+45.0000)/- 6.50000)) ALGBRC(13) = 1.36000/(( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))+ 0.0800000*EXP(- STATES(1)/11.0000)) ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000)) ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 0.453700*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 3.49000/( 0.135000*EXP(- (STATES(1)+80.0000)/6.80000)+ 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))) ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000)) ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 11.6300*(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))/EXP( - 2.53500e-07*STATES(1)), 3.49000/( ((STATES(1)+37.7800)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))*( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))+( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))) ALGBRC(16) = 1000.00/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000))) ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200)) ALGBRC(17) = 550.000*EXP(- (STATES(1)+70.0000)/25.0000 ** 2.00000)+49.0000 ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410)) ALGBRC(18) = 3300.00*EXP(- (STATES(1)+70.0000)/30.0000 ** 2.00000)+49.0000 ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410)) ALGBRC(19) = 10000.0/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000))) ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+11.5000)/- 11.8200)) ALGBRC(20) = 1000.00/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900)) ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800)) ALGBRC(11) = CONSTS(101) - STATES(19) ALGBRC(22) = CONSTS(124)*ALGBRC(11)/CONSTS(169) ** CONSTS(129) ALGBRC(24) = CONSTS(125)+( CONSTS(126)*STATES(23) ** CONSTS(127) - 1.00000)/(STATES(23) ** CONSTS(127)+CONSTS(128) ** CONSTS(127)) ALGBRC(26) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(24)/STATES(6)) ALGBRC(31) = ( CONSTS(146)*STATES(7)*( CONSTS(22)*STATES(8)+ CONSTS(23)*STATES(9))*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(32) = ( CONSTS(25)*STATES(10)*STATES(11)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(33) = (( (0.0480000/(EXP((STATES(1)+37.0000)/25.0000)+EXP((STATES(1)+37.0000)/- 25.0000))+0.0100000)*0.00100000)/(1.00000+EXP((STATES(1) - (ALGBRC(26)+76.7700))/- 17.0000))+( CONSTS(26)*(STATES(1) - (ALGBRC(26)+1.73000)))/( (1.00000+EXP(( 1.61300*CONSTS(10)*(STATES(1) - (ALGBRC(26)+1.73000)))/( CONSTS(8)*CONSTS(9))))*(1.00000+EXP((CONSTS(24) - 0.998800)/- 0.124000))))/CONSTS(7) ALGBRC(40) = 1.00000/(1.00000+CONSTS(35)/STATES(2) ** 4.00000) ALGBRC(41) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))+ 0.0365000*CONSTS(149)*EXP(( - STATES(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))) ALGBRC(42) = (( (CONSTS(33)/CONSTS(7))*ALGBRC(41)*CONSTS(24))/(CONSTS(24)+CONSTS(34)))*ALGBRC(40) ALGBRC(10) = TERNRY(VOI - INT(VOI/CONSTS(12))*CONSTS(12).GE.0.00000.AND.VOI - INT(VOI/CONSTS(12))*CONSTS(12).LE.CONSTS(13), CONSTS(14)/CONSTS(7), 0.00000) ALGBRC(39) = ( CONSTS(32)*CONSTS(30)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(28) = ( CONSTS(16)*CONSTS(142)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(30) = ( (CONSTS(18)/(1.00000+EXP(- (10.0000+STATES(1))/45.0000)))*(STATES(1) - ALGBRC(26))*CONSTS(144))/CONSTS(7) ALGBRC(35) = ( CONSTS(27)*STATES(12)*CONSTS(148)*(STATES(1) - ALGBRC(26)))/CONSTS(7) ALGBRC(73) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(STATES(15)/CONSTS(46)) ALGBRC(74) = CONSTS(83)*(STATES(1) - ALGBRC(73)) ALGBRC(53) = 1000.00*10.0000 ** - 14.0000+STATES(14) ALGBRC(54) = 1.00000+CONSTS(47)/STATES(15)+( CONSTS(47)*ALGBRC(53))/( STATES(15)*CONSTS(48)) ALGBRC(55) = 1.00000+CONSTS(48)/ALGBRC(53)+( CONSTS(48)*STATES(15))/( ALGBRC(53)*CONSTS(47)) ALGBRC(56) = 1.00000/(CONSTS(175)+( ALGBRC(55)*(CONSTS(49)+( CONSTS(51)*CONSTS(175))/CONSTS(176)))/(CONSTS(50)+( CONSTS(173)*ALGBRC(55))/ALGBRC(54))) ALGBRC(57) = 1.00000/(ALGBRC(55)+( CONSTS(175)*(CONSTS(50)+( CONSTS(173)*ALGBRC(55))/ALGBRC(54)))/(CONSTS(49)+( CONSTS(51)*CONSTS(175))/CONSTS(176))) ALGBRC(58) = (( CONSTS(4)*( CONSTS(50)*ALGBRC(57) - CONSTS(49)*ALGBRC(56)))/( 60.0000*1000.00))*CONSTS(52) ALGBRC(67) = 1.00000+CONSTS(63)/STATES(15)+( CONSTS(63)*STATES(16))/( STATES(15)*CONSTS(64)) ALGBRC(68) = 1.00000+CONSTS(64)/STATES(16)+( CONSTS(64)*STATES(15))/( STATES(16)*CONSTS(63)) ALGBRC(69) = 1.00000/(CONSTS(187)+( ALGBRC(68)*(CONSTS(65)+( CONSTS(67)*CONSTS(187))/CONSTS(188)))/(CONSTS(66)+( CONSTS(181)*ALGBRC(68))/ALGBRC(67))) ALGBRC(70) = 1.00000/(ALGBRC(68)+( CONSTS(187)*(CONSTS(66)+( CONSTS(181)*ALGBRC(68))/ALGBRC(67)))/(CONSTS(65)+( CONSTS(67)*CONSTS(187))/CONSTS(188))) ALGBRC(66) = 1000.00*10.0000 ** - STATES(14) ALGBRC(71) = ( (CONSTS(180) ** CONSTS(69)/(ALGBRC(66) ** CONSTS(69)+CONSTS(180) ** CONSTS(69)))*CONSTS(182) ** CONSTS(71))/(CONSTS(182) ** CONSTS(71)+CONSTS(70) ** CONSTS(71)) ALGBRC(72) = (( CONSTS(4)*ALGBRC(71)*( CONSTS(66)*ALGBRC(70) - CONSTS(65)*ALGBRC(69)))/( 60.0000*1000.00))*CONSTS(72) ALGBRC(59) = 1.00000+CONSTS(54)/STATES(16)+( CONSTS(53)*CONSTS(54))/( STATES(16)*STATES(2)) ALGBRC(60) = 1.00000+STATES(2)/CONSTS(53)+( STATES(2)*STATES(16))/( CONSTS(53)*CONSTS(54)) ALGBRC(62) = 1.00000/(CONSTS(185)+( ALGBRC(60)*(CONSTS(55)+( CONSTS(57)*CONSTS(185))/CONSTS(186)))/(CONSTS(56)+( CONSTS(177)*ALGBRC(60))/ALGBRC(59))) ALGBRC(63) = 1.00000/(ALGBRC(60)+( CONSTS(185)*(CONSTS(56)+( CONSTS(177)*ALGBRC(60))/ALGBRC(59)))/(CONSTS(55)+( CONSTS(57)*CONSTS(185))/CONSTS(186))) ALGBRC(61) = 1000.00*10.0000 ** - STATES(14) ALGBRC(64) = (ALGBRC(61) ** CONSTS(59)/(ALGBRC(61) ** CONSTS(59)+CONSTS(58) ** CONSTS(59)))*(1.00000 - CONSTS(178) ** CONSTS(61)/(CONSTS(178) ** CONSTS(61)+CONSTS(60) ** CONSTS(61))) ALGBRC(65) = (( CONSTS(4)*ALGBRC(64)*( CONSTS(56)*ALGBRC(63) - CONSTS(55)*ALGBRC(62)))/( 60.0000*1000.00))*CONSTS(62) ALGBRC(80) = (CONSTS(4)+CONSTS(5))*( CONSTS(76)*STATES(13) - CONSTS(77)*STATES(16)*10.0000 ** - STATES(14)*1000.00) ALGBRC(43) = CONSTS(4)*CONSTS(141)*CONSTS(36)*(CONSTS(183) - STATES(13)) ALGBRC(75) = 1000.00*10.0000 ** - STATES(14) ALGBRC(76) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(151)/ALGBRC(75)) ALGBRC(77) = CONSTS(84)*(STATES(1) - ALGBRC(76)) ALGBRC(47) = 10.0000 ** - STATES(14)*1000.00 ALGBRC(48) = ( CONSTS(39)*STATES(2)*CONSTS(41))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*STATES(2)+ STATES(2)*ALGBRC(47)+ CONSTS(38)*ALGBRC(47)) ALGBRC(49) = ( CONSTS(38)*ALGBRC(47)*CONSTS(42))/( CONSTS(38)*CONSTS(39)+ CONSTS(39)*STATES(2)+ STATES(2)*ALGBRC(47)+ CONSTS(38)*ALGBRC(47)) ALGBRC(50) = ALGBRC(47) ** CONSTS(44)/(ALGBRC(47) ** CONSTS(44)+CONSTS(161) ** CONSTS(44)) ALGBRC(51) = ( ALGBRC(50)*( CONSTS(170)*ALGBRC(49) - ALGBRC(48)*CONSTS(172)))/(CONSTS(170)+ALGBRC(49)+ALGBRC(48)+CONSTS(172)) ALGBRC(52) = ALGBRC(51)*CONSTS(4)*CONSTS(45) ALGBRC(78) = log(10.0000)*(10.0000 ** - STATES(14)+( 10.0000 ** STATES(14)+CONSTS(78)*CONSTS(79))/10.0000 ** STATES(14)+10.0000 ** CONSTS(78) ** 2.00000+( 10.0000 ** STATES(14)+CONSTS(80)*CONSTS(81))/10.0000 ** STATES(14)+10.0000 ** CONSTS(80) ** 2.00000) ALGBRC(21) = (( CONSTS(8)*CONSTS(9))/CONSTS(10))*log(CONSTS(20)/STATES(2)) ALGBRC(25) = ( CONSTS(145)*STATES(3) ** 3.00000*STATES(4)*STATES(5)*(STATES(1) - ALGBRC(21)))/CONSTS(7) ALGBRC(85) = EXP( CONSTS(89)*0.500000*STATES(1)*CONSTS(93)) ALGBRC(86) = EXP( CONSTS(89)*0.500000*STATES(1)*(CONSTS(93) - 1.00000)) ALGBRC(87) = STATES(2) ** 3.00000 ALGBRC(88) = ( (( (CONSTS(95)/( (CONSTS(152)+CONSTS(91) ** 3.00000)*(CONSTS(86)+CONSTS(92))))*( ALGBRC(85)*ALGBRC(87)*CONSTS(86) - ALGBRC(86)*CONSTS(152)*STATES(17)))/(1.00000+ CONSTS(94)*ALGBRC(86)))*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(38) = ( CONSTS(31)*CONSTS(29)*(STATES(1) - ALGBRC(21)))/CONSTS(7) ALGBRC(23) = (( CONSTS(16)*CONSTS(142)*(STATES(1) - ALGBRC(21)))/CONSTS(7))*CONSTS(143) ALGBRC(34) = ( CONSTS(27)*STATES(12)*CONSTS(28)*(STATES(1) - ALGBRC(21)))/CONSTS(7) ALGBRC(37) = ALGBRC(34)+ALGBRC(35) ALGBRC(93) = CONSTS(89)*STATES(1) ALGBRC(94) = EXP(- ALGBRC(93)) ALGBRC(102) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( (( CONSTS(104)*ALGBRC(93))/(1.00000 - ALGBRC(94)))*(( CONSTS(86)*ALGBRC(94) - STATES(17))+ (CONSTS(103)/CONSTS(102))*( CONSTS(86)*ALGBRC(94) - STATES(18))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( (( CONSTS(104)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(86)*EXP(- 1.00000e-05) - STATES(17))+ (CONSTS(103)/CONSTS(102))*( CONSTS(86)*EXP(- 1.00000e-05) - STATES(18))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(101) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( (( CONSTS(104)*ALGBRC(93))/(1.00000 - ALGBRC(94)))*( CONSTS(86)*ALGBRC(94) - STATES(17)))/(1.00000+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( (( CONSTS(104)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(86)*EXP(- 1.00000e-05) - STATES(17)))/(1.00000+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(95) = EXP((STATES(1) - CONSTS(105))/CONSTS(106)) ALGBRC(106) = ALGBRC(95)/( CONSTS(107)*(ALGBRC(95)+1.00000)) ALGBRC(108) = ( (1.00000/CONSTS(153))*STATES(17) ** 2.00000)/(STATES(17) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(98) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-09, (STATES(17)+( (CONSTS(104)/CONSTS(102))*CONSTS(86)*ALGBRC(93)*ALGBRC(94))/(1.00000 - ALGBRC(94)))/(1.00000+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), (STATES(17)+ (CONSTS(104)/CONSTS(102))*CONSTS(86))/(1.00000+CONSTS(104)/CONSTS(102))) ALGBRC(107) = ( (1.00000/CONSTS(153))*ALGBRC(98) ** 2.00000)/(ALGBRC(98) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(113) = (ALGBRC(106)+CONSTS(158))*( (CONSTS(168)+ALGBRC(107)+CONSTS(158))*(CONSTS(168)+ALGBRC(108))+ ALGBRC(106)*(CONSTS(168)+ALGBRC(107))) ALGBRC(114) = ( ALGBRC(106)*CONSTS(168)*(ALGBRC(106)+CONSTS(158)+CONSTS(168)+ALGBRC(108)))/ALGBRC(113) ALGBRC(116) = ( ALGBRC(106)*( ALGBRC(107)*(ALGBRC(106)+CONSTS(168)+ALGBRC(108))+ ALGBRC(108)*CONSTS(158)))/ALGBRC(113) ALGBRC(122) = ALGBRC(102)*ALGBRC(116)+ ALGBRC(101)*ALGBRC(114) ALGBRC(124) = ( ALGBRC(101)*ALGBRC(106))/(ALGBRC(106)+CONSTS(158)) ALGBRC(126) = ( ( STATES(20)*ALGBRC(122)+ STATES(21)*ALGBRC(124))*CONSTS(6))/CONSTS(4) ALGBRC(128) = ( - ALGBRC(126)*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(83) = arbitrary_log(CONSTS(86)/STATES(17), 10)/CONSTS(89) ALGBRC(84) = ( CONSTS(90)*(STATES(1) - ALGBRC(83))*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(82) = ( (( CONSTS(87)*STATES(17))/(CONSTS(88)+STATES(17)))*2.00000*CONSTS(4)*CONSTS(10))/CONSTS(7) ALGBRC(109) = ( (1.00000/CONSTS(114))*(ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/(ALGBRC(98) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(110) = ( (1.00000/CONSTS(114))*(STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/(STATES(17) ** 2.00000+CONSTS(116) ** 2.00000) ALGBRC(117) = ( CONSTS(158)*CONSTS(168)*(CONSTS(158)+ALGBRC(106)+CONSTS(168)+ALGBRC(107)))/ALGBRC(113) ALGBRC(119) = ALGBRC(114)*ALGBRC(109)+ ALGBRC(117)*ALGBRC(110) ALGBRC(111) = ( (CONSTS(115)/CONSTS(114))*CONSTS(118)*(ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/( CONSTS(118)*ALGBRC(98) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000) ALGBRC(112) = ( (CONSTS(115)/CONSTS(114))*CONSTS(118)*(STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000))/( CONSTS(118)*STATES(17) ** 2.00000+ CONSTS(117)*CONSTS(116) ** 2.00000) ALGBRC(121) = ( ALGBRC(106)*ALGBRC(111)+ CONSTS(158)*ALGBRC(112))/(ALGBRC(106)+CONSTS(158)) ALGBRC(97) = (STATES(17)+ (CONSTS(103)/CONSTS(102))*STATES(18))/(1.00000+CONSTS(103)/CONSTS(102)) ALGBRC(103) = ( (( (1.00000/CONSTS(110))*ALGBRC(97))/CONSTS(111))*(ALGBRC(95)+CONSTS(108)))/(ALGBRC(95)+1.00000) ALGBRC(104) = ( (( (1.00000/CONSTS(110))*STATES(17))/CONSTS(111))*(ALGBRC(95)+CONSTS(108)))/(ALGBRC(95)+1.00000) ALGBRC(115) = ( CONSTS(158)*( ALGBRC(108)*(CONSTS(158)+CONSTS(168)+ALGBRC(107))+ ALGBRC(107)*ALGBRC(106)))/ALGBRC(113) ALGBRC(127) = ALGBRC(115)*ALGBRC(103)+ ALGBRC(117)*ALGBRC(104) ALGBRC(105) = ( (1.00000/CONSTS(110))*CONSTS(109)*(ALGBRC(95)+CONSTS(108)))/( CONSTS(109)*ALGBRC(95)+CONSTS(108)) ALGBRC(129) = ALGBRC(105) ALGBRC(91) = ( CONSTS(96)*STATES(17) ** 2.00000)/(CONSTS(97) ** 2.00000+STATES(17) ** 2.00000) ALGBRC(99) = ( CONSTS(103)*(STATES(18) - STATES(17)))/(1.00000+CONSTS(103)/CONSTS(102)) ALGBRC(100) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-05, ( CONSTS(103)*((STATES(18) - STATES(17))+ (( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94)))*(STATES(18) - CONSTS(86)*ALGBRC(94))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), ( CONSTS(103)*((STATES(18) - STATES(17))+ (( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(18) - CONSTS(86)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))) ALGBRC(118) = ALGBRC(116)*ALGBRC(100)+ ALGBRC(99)*ALGBRC(115) ALGBRC(120) = ( ALGBRC(99)*ALGBRC(108))/(CONSTS(168)+ALGBRC(108)) ALGBRC(130) = ( ( STATES(20)*ALGBRC(118)+ STATES(22)*ALGBRC(120))*CONSTS(6))/CONSTS(4) ALGBRC(132) = (- ALGBRC(130)+ALGBRC(91)) - CONSTS(100)*(STATES(18) - STATES(17)) ALGBRC(131) = ( CONSTS(158)*ALGBRC(104))/(ALGBRC(106)+CONSTS(158)) ALGBRC(133) = ALGBRC(105) ALGBRC(134) = ((1.00000 - STATES(20)) - STATES(21)) - STATES(22) ALGBRC(123) = ( CONSTS(168)*ALGBRC(110))/(CONSTS(168)+ALGBRC(108)) ALGBRC(125) = ALGBRC(112) ALGBRC(135) = ( CONSTS(123)*STATES(23))/CONSTS(171) ALGBRC(136) = ALGBRC(135)*CONSTS(160) ALGBRC(137) = STATES(24)+STATES(25)+STATES(26) ALGBRC(138) = TERNRY(ALGBRC(137).LT.0.00000, ( ALGBRC(136)*( CONSTS(134)*ALGBRC(137)+1.00000))/(1.00000 - ALGBRC(137)), ( ALGBRC(136)*(1.00000+ (CONSTS(134)+2.00000)*ALGBRC(137)))/(1.00000+ALGBRC(137))) ALGBRC(139) = TERNRY(1.00000 - ALGBRC(138)/( CONSTS(121)*CONSTS(123)).GT.0.100000, CONSTS(120)*(1.00000 - ALGBRC(138)/( CONSTS(121)*CONSTS(123))), CONSTS(120)*0.100000) ALGBRC(140) = (CONSTS(101) - STATES(19))*ALGBRC(139) - STATES(17)*STATES(19)*CONSTS(119) ALGBRC(12) = STATES(17) ALGBRC(27) = ( ALGBRC(25)*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(29) = ALGBRC(28)+ALGBRC(23) ALGBRC(36) = ( ALGBRC(34)*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(44) = ( ALGBRC(42)*3.00000*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(45) = ( ALGBRC(42)*- 2.00000*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(46) = ( - (ALGBRC(28)+ALGBRC(32)+ALGBRC(30)+ALGBRC(31)+ALGBRC(33)+ALGBRC(35)+ ALGBRC(42)*- 2.00000+ALGBRC(39))*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(79) = ALGBRC(52)/CONSTS(4) ALGBRC(81) = ALGBRC(65)/CONSTS(4) ALGBRC(89) = ( - (ALGBRC(38)+ALGBRC(23)+ALGBRC(25)+ ALGBRC(88)*3.00000+ ALGBRC(42)*3.00000+ALGBRC(34))*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(90) = ( ALGBRC(88)*3.00000*CONSTS(7))/( CONSTS(4)*CONSTS(10)) ALGBRC(92) = ALGBRC(89)+ALGBRC(79)+ALGBRC(81) ALGBRC(96) = TERNRY(ABS(ALGBRC(93)).GT.1.00000e-09, (STATES(17)+ (CONSTS(103)/CONSTS(102))*STATES(18)+( (CONSTS(104)/CONSTS(102))*CONSTS(86)*ALGBRC(93)*ALGBRC(94))/(1.00000 - ALGBRC(94)))/(1.00000+CONSTS(103)/CONSTS(102)+( (CONSTS(104)/CONSTS(102))*ALGBRC(93))/(1.00000 - ALGBRC(94))), (STATES(17)+ (CONSTS(103)/CONSTS(102))*STATES(18)+ (CONSTS(104)/CONSTS(102))*CONSTS(86))/(1.00000+CONSTS(103)/CONSTS(102)+CONSTS(104)/CONSTS(102))) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END