Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 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