C C There are a total of 296 entries in the algebraic variable array. C There are a total of 80 entries in each of the rate and state variable arrays. C There are a total of 770 entries in the constant variable array. C C C VOI is time in component environment (second). C ALGBRC(8) is I_CaT in component I_CaT (picoampere). C CONSTS(1) is g_CaT in component I_CaT (picosiemens). C CONSTS(2) is ECaT in component I_CaT (millivolt). C CONSTS(3) is alpha_OT in component I_CaT (first_order_rate_constant). C CONSTS(4) is beta_OT in component I_CaT (first_order_rate_constant). C STATES(1) is OCaT in component I_CaT (dimensionless). C STATES(2) is Vm in component membrane (millivolt). C STATES(3) is dT in component I_CaT_d_gate (dimensionless). C STATES(4) is fT in component I_CaT_f_gate (dimensionless). C ALGBRC(1) is infinity_dT in component I_CaT_d_gate (dimensionless). C CONSTS(688) is tau_dT1 in component I_CaT_d_gate (second). C CONSTS(5) is A_dT1 in component I_CaT_d_gate (second). C CONSTS(6) is k_dT in component I_CaT_d_gate (per_millivolt). C CONSTS(7) is V_dT in component I_CaT_d_gate (millivolt). C ALGBRC(2) is infinity_fT in component I_CaT_f_gate (dimensionless). C ALGBRC(20) is tau_fT1 in component I_CaT_f_gate (second). C CONSTS(8) is A_fT1 in component I_CaT_f_gate (second). C CONSTS(9) is A_fT2 in component I_CaT_f_gate (second). C CONSTS(10) is A_fT3 in component I_CaT_f_gate (per_millivolt). C CONSTS(11) is A_fT4 in component I_CaT_f_gate (millivolt). C CONSTS(12) is k_fT in component I_CaT_f_gate (per_millivolt). C CONSTS(13) is V_fT in component I_CaT_f_gate (millivolt). C ALGBRC(26) is I_CaExt in component I_CaExt (picoampere). C CONSTS(14) is gCaExt in component I_CaExt (femtoampere). C CONSTS(15) is KCaExt in component I_CaExt (micromolar). C STATES(5) is CCy in component CCy (micromolar). C ALGBRC(30) is I_Kv1_1 in component I_Kv1_1 (picoampere). C CONSTS(16) is g_Kv1_1 in component I_Kv1_1 (picosiemens). C CONSTS(689) is EK in component Nerst_potentials (millivolt). C STATES(6) is dv1_1 in component I_Kv1_1_d_gate (dimensionless). C STATES(7) is fv1_1 in component I_Kv1_1_f_gate (dimensionless). C ALGBRC(21) is alpha_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant). C ALGBRC(3) is A_dv1_11 in component I_Kv1_1_d_gate (first_order_rate_constant). C CONSTS(17) is A_dv1_11a in component I_Kv1_1_d_gate (first_order_rate_constant). C CONSTS(18) is A_dv1_11b in component I_Kv1_1_d_gate (first_order_rate_constant). C CONSTS(19) is A_dv1_11c in component I_Kv1_1_d_gate (per_millivolt). C CONSTS(20) is A_dv1_11d in component I_Kv1_1_d_gate (millivolt). C CONSTS(21) is A_dv1_12 in component I_Kv1_1_d_gate (dimensionless). C CONSTS(22) is A_dv1_13 in component I_Kv1_1_d_gate (per_millivolt). C CONSTS(23) is A_dv1_14 in component I_Kv1_1_d_gate (millivolt). C ALGBRC(27) is beta_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant). C ALGBRC(4) is alpha_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant). C CONSTS(24) is A_fv1_11 in component I_Kv1_1_f_gate (first_order_rate_constant). C CONSTS(25) is A_fv1_12 in component I_Kv1_1_f_gate (dimensionless). C CONSTS(26) is A_fv1_13 in component I_Kv1_1_f_gate (per_millivolt). C CONSTS(27) is A_fv1_14 in component I_Kv1_1_f_gate (millivolt). C ALGBRC(22) is beta_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant). C ALGBRC(31) is I_KERG in component I_KERG (picoampere). C CONSTS(28) is g_KERG in component I_KERG (picosiemens). C STATES(8) is dERG in component I_KERG_d_gate (dimensionless). C ALGBRC(5) is alpha_dERG in component I_KERG_d_gate (first_order_rate_constant). C CONSTS(29) is A_dERG1 in component I_KERG_d_gate (first_order_rate_constant). C CONSTS(30) is A_dERG2 in component I_KERG_d_gate (dimensionless). C CONSTS(31) is A_dERG3 in component I_KERG_d_gate (per_millivolt). C CONSTS(32) is A_dERG4 in component I_KERG_d_gate (millivolt). C ALGBRC(23) is beta_dERG in component I_KERG_d_gate (first_order_rate_constant). C ALGBRC(32) is I_KB in component I_KB (picoampere). C CONSTS(33) is g_KB in component I_KB (picosiemens). C ALGBRC(33) is I_Na in component I_Na (picoampere). C CONSTS(34) is g_Na in component I_Na (picosiemens). C CONSTS(690) is ENa_Cy in component Nerst_potentials (millivolt). C STATES(9) is dNa in component I_Na_d_gate (dimensionless). C STATES(10) is fNa in component I_Na_f_gate (dimensionless). C ALGBRC(24) is alpha_dNa in component I_Na_d_gate (first_order_rate_constant). C ALGBRC(6) is A_dNa1 in component I_Na_d_gate (first_order_rate_constant). C CONSTS(35) is A_dNa1a in component I_Na_d_gate (first_order_rate_constant). C CONSTS(36) is A_dNa1b in component I_Na_d_gate (first_order_rate_constant). C CONSTS(37) is A_dNa1c in component I_Na_d_gate (per_millivolt). C CONSTS(38) is A_dNa1d in component I_Na_d_gate (millivolt). C CONSTS(39) is A_dNa2 in component I_Na_d_gate (dimensionless). C CONSTS(40) is A_dNa3 in component I_Na_d_gate (per_millivolt). C CONSTS(41) is A_dNa4 in component I_Na_d_gate (millivolt). C ALGBRC(28) is beta_dNa in component I_Na_d_gate (first_order_rate_constant). C ALGBRC(25) is alpha_fNa in component I_Na_f_gate (first_order_rate_constant). C ALGBRC(7) is A_fNa1 in component I_Na_f_gate (first_order_rate_constant). C CONSTS(42) is A_fNa1a in component I_Na_f_gate (first_order_rate_constant). C CONSTS(43) is A_fNa1b in component I_Na_f_gate (first_order_rate_constant). C CONSTS(44) is A_fNa1c in component I_Na_f_gate (per_millivolt). C CONSTS(45) is A_fNa1d in component I_Na_f_gate (millivolt). C CONSTS(46) is A_fNa2 in component I_Na_f_gate (dimensionless). C CONSTS(47) is A_fNa3 in component I_Na_f_gate (per_millivolt). C CONSTS(48) is A_fNa4 in component I_Na_f_gate (millivolt). C ALGBRC(29) is beta_fNa in component I_Na_f_gate (first_order_rate_constant). C ALGBRC(34) is I_L in component I_L (picoampere). C CONSTS(49) is g_L in component I_L (picosiemens). C CONSTS(50) is EL in component I_L (millivolt). C ALGBRC(35) is I_ion_Cy in component membrane (picoampere). C CONSTS(51) is Cm in component membrane (picofarad). C ALGBRC(73) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(97) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(121) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(145) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(169) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(193) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(217) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(241) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(265) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(288) is I_ion_PU in component PU_membrane (picoampere). C ALGBRC(9) is ECa_Cy in component Nerst_potentials (millivolt). C CONSTS(52) is T in component model_parameters (kelvin). C CONSTS(53) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(54) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(55) is CO in component model_parameters (micromolar). C CONSTS(56) is KO in component model_parameters (micromolar). C CONSTS(57) is Ki in component model_parameters (micromolar). C CONSTS(58) is NO in component model_parameters (micromolar). C CONSTS(59) is Ni in component model_parameters (micromolar). C CONSTS(60) is lambda_S2_Cy in component CCy (dimensionless). C CONSTS(61) is delta_SCy in component CCy (micromolar_per_coulomb). C ALGBRC(46) is JCy in component JCy (flux). C ALGBRC(36) is JS2Cy1 in component JS2Cy (flux). C ALGBRC(37) is JS2Cy2 in component JS2Cy (flux). C ALGBRC(38) is JS2Cy3 in component JS2Cy (flux). C ALGBRC(39) is JS2Cy4 in component JS2Cy (flux). C ALGBRC(40) is JS2Cy5 in component JS2Cy (flux). C ALGBRC(41) is JS2Cy6 in component JS2Cy (flux). C ALGBRC(42) is JS2Cy7 in component JS2Cy (flux). C ALGBRC(43) is JS2Cy8 in component JS2Cy (flux). C ALGBRC(44) is JS2Cy9 in component JS2Cy (flux). C ALGBRC(45) is JS2Cy10 in component JS2Cy (flux). C CONSTS(62) is mu_S2Cy1 in component JS2Cy (first_order_rate_constant). C CONSTS(63) is mu_S2Cy2 in component JS2Cy (first_order_rate_constant). C CONSTS(64) is mu_S2Cy3 in component JS2Cy (first_order_rate_constant). C CONSTS(65) is mu_S2Cy4 in component JS2Cy (first_order_rate_constant). C CONSTS(66) is mu_S2Cy5 in component JS2Cy (first_order_rate_constant). C CONSTS(67) is mu_S2Cy6 in component JS2Cy (first_order_rate_constant). C CONSTS(68) is mu_S2Cy7 in component JS2Cy (first_order_rate_constant). C CONSTS(69) is mu_S2Cy8 in component JS2Cy (first_order_rate_constant). C CONSTS(70) is mu_S2Cy9 in component JS2Cy (first_order_rate_constant). C CONSTS(71) is mu_S2Cy10 in component JS2Cy (first_order_rate_constant). C STATES(11) is CS2 in component CS2 (micromolar). C STATES(12) is CS2 in component CS2 (micromolar). C STATES(13) is CS2 in component CS2 (micromolar). C STATES(14) is CS2 in component CS2 (micromolar). C STATES(15) is CS2 in component CS2 (micromolar). C STATES(16) is CS2 in component CS2 (micromolar). C STATES(17) is CS2 in component CS2 (micromolar). C STATES(18) is CS2 in component CS2 (micromolar). C STATES(19) is CS2 in component CS2 (micromolar). C STATES(20) is CS2 in component CS2 (micromolar). C ALGBRC(47) is JS1S2_1 in component JS1S2 (flux). C ALGBRC(48) is JS1S2_2 in component JS1S2 (flux). C ALGBRC(49) is JS1S2_3 in component JS1S2 (flux). C ALGBRC(50) is JS1S2_4 in component JS1S2 (flux). C ALGBRC(51) is JS1S2_5 in component JS1S2 (flux). C ALGBRC(52) is JS1S2_6 in component JS1S2 (flux). C ALGBRC(53) is JS1S2_7 in component JS1S2 (flux). C ALGBRC(54) is JS1S2_8 in component JS1S2 (flux). C ALGBRC(55) is JS1S2_9 in component JS1S2 (flux). C ALGBRC(56) is JS1S2_10 in component JS1S2 (flux). C CONSTS(692) is mu_S1S2_1 in component JS1S2 (first_order_rate_constant). C CONSTS(693) is mu_S1S2_2 in component JS1S2 (first_order_rate_constant). C CONSTS(694) is mu_S1S2_3 in component JS1S2 (first_order_rate_constant). C CONSTS(695) is mu_S1S2_4 in component JS1S2 (first_order_rate_constant). C CONSTS(696) is mu_S1S2_5 in component JS1S2 (first_order_rate_constant). C CONSTS(697) is mu_S1S2_6 in component JS1S2 (first_order_rate_constant). C CONSTS(698) is mu_S1S2_7 in component JS1S2 (first_order_rate_constant). C CONSTS(699) is mu_S1S2_8 in component JS1S2 (first_order_rate_constant). C CONSTS(700) is mu_S1S2_9 in component JS1S2 (first_order_rate_constant). C CONSTS(701) is mu_S1S2_10 in component JS1S2 (first_order_rate_constant). C CONSTS(72) is mu_A in component JS1S2 (first_order_rate_constant). C CONSTS(73) is mu_B in component JS1S2 (first_order_rate_constant). C CONSTS(74) is n_PU in component model_parameters (dimensionless). C STATES(21) is CS1 in component CS1 (micromolar). C STATES(22) is CS1 in component CS1 (micromolar). C STATES(23) is CS1 in component CS1 (micromolar). C STATES(24) is CS1 in component CS1 (micromolar). C STATES(25) is CS1 in component CS1 (micromolar). C STATES(26) is CS1 in component CS1 (micromolar). C STATES(27) is CS1 in component CS1 (micromolar). C STATES(28) is CS1 in component CS1 (micromolar). C STATES(29) is CS1 in component CS1 (micromolar). C STATES(30) is CS1 in component CS1 (micromolar). C CONSTS(75) is mu_Cy in component JCy (first_order_rate_constant). C CONSTS(76) is C_infinity in component JCy (micromolar). C CONSTS(77) is P in component model_parameters (micromolar). C ALGBRC(67) is I_iCa in component PU_membrane (picoampere). C ALGBRC(71) is I_iNa in component PU_membrane (picoampere). C CONSTS(716) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(59) is I_Ca in component I_Ca (picoampere). C ALGBRC(61) is I_Na in component I_Na (picoampere). C ALGBRC(69) is I_NaP in component I_NaP (picoampere). C ALGBRC(63) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(66) is I_PM in component I_PM (picoampere). C ALGBRC(65) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(78) is gCa_ in component I_Ca (picosiemens). C CONSTS(79) is kCa in component I_Ca (per_millivolt). C CONSTS(80) is kVCa in component I_Ca (per_millivolt). C CONSTS(81) is VhCa in component I_Ca (millivolt). C ALGBRC(57) is gCa in component I_Ca (picosiemens). C ALGBRC(58) is ECa_PU in component I_Ca (millivolt). C CONSTS(82) is T in component model_parameters (kelvin). C CONSTS(83) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(84) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(85) is CO in component model_parameters (micromolar). C CONSTS(86) is gNa in component I_Na (picosiemens). C ALGBRC(60) is ENa_PU in component I_Na (millivolt). C CONSTS(87) is NO in component model_parameters (micromolar). C STATES(31) is NS1 in component NS1 (micromolar). C CONSTS(88) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(62) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(89) is hNSCC in component model_parameters (dimensionless). C CONSTS(90) is ENSCC in component model_parameters (millivolt). C CONSTS(91) is KNSCC in component model_parameters (micromolar). C CONSTS(92) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(64) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(93) is gPM in component I_PM (femtoampere). C CONSTS(94) is KPM in component I_PM (micromolar). C CONSTS(95) is gNaP in component I_NaP (picosiemens). C CONSTS(96) is KNaP in component I_NaP (micromolar). C CONSTS(97) is hNaP in component I_NaP (dimensionless). C CONSTS(98) is ENaP in component I_NaP (millivolt). C ALGBRC(68) is JSERCA in component JSERCA (flux). C CONSTS(99) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(100) is A2 in component JSERCA (dimensionless). C CONSTS(101) is A4 in component JSERCA (per_micromolar). C CONSTS(102) is A5 in component JSERCA (per_micromolar). C CONSTS(103) is A6 in component JSERCA (per_micromolar2). C STATES(32) is CER in component CER (micromolar). C ALGBRC(72) is JMCU in component JMCU (flux). C CONSTS(104) is VMCU in component JMCU (flux). C CONSTS(105) is KMCU in component JMCU (micromolar). C ALGBRC(70) is epsilon_INH in component JMCU (dimensionless). C CONSTS(106) is KINH in component JMCU (micromolar). C CONSTS(107) is hINH in component JMCU (dimensionless). C STATES(33) is CMT in component CMT (micromolar). C ALGBRC(80) is JIPR in component JIPR (flux). C CONSTS(108) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(109) is k_1 in component JIPR (flux). C CONSTS(110) is k1 in component JIPR (first_order_rate_constant). C CONSTS(111) is k2 in component JIPR (first_order_rate_constant). C CONSTS(112) is r2 in component JIPR (first_order_rate_constant). C CONSTS(113) is r_2 in component JIPR (flux). C CONSTS(114) is r4 in component JIPR (first_order_rate_constant). C CONSTS(115) is R1 in component JIPR (micromolar). C CONSTS(116) is R3 in component JIPR (micromolar). C ALGBRC(74) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(77) is phi_1 in component JIPR (flux). C ALGBRC(81) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(83) is phi3 in component JIPR (first_order_rate_constant). C STATES(34) is H in component JIPR (dimensionless). C CONSTS(117) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(118) is h_beta in component JIPR (dimensionless). C CONSTS(119) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(120) is K_beta in component JIPR (micromolar). C CONSTS(702) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(10) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(35) is zeta in component JIPR (dimensionless). C CONSTS(121) is K_phi3_act in component JIPR (micromolar). C CONSTS(122) is K_phi3_inh in component JIPR (micromolar). C CONSTS(123) is h_phi3_act in component JIPR (dimensionless). C CONSTS(124) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(125) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(126) is P in component model_parameters (micromolar). C ALGBRC(75) is JNCX in component JNCX (flux). C CONSTS(127) is VNCX in component JNCX (flux). C CONSTS(128) is KNCX in component JNCX (micromolar). C CONSTS(691) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(712) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(129) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(130) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(131) is gamma_MT in component model_parameters (dimensionless). C CONSTS(132) is gamma_ER in component model_parameters (dimensionless). C CONSTS(713) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(714) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(715) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(133) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(78) is fm in component CMT (dimensionless). C CONSTS(134) is Km in component CMT (micromolar). C CONSTS(135) is Bm in component CMT (micromolar). C CONSTS(136) is KO in component model_parameters (micromolar). C CONSTS(137) is Ki in component model_parameters (micromolar). C CONSTS(138) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(91) is I_iCa in component PU_membrane (picoampere). C ALGBRC(95) is I_iNa in component PU_membrane (picoampere). C CONSTS(722) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(82) is I_Ca in component I_Ca (picoampere). C ALGBRC(85) is I_Na in component I_Na (picoampere). C ALGBRC(93) is I_NaP in component I_NaP (picoampere). C ALGBRC(87) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(90) is I_PM in component I_PM (picoampere). C ALGBRC(89) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(139) is gCa_ in component I_Ca (picosiemens). C CONSTS(140) is kCa in component I_Ca (per_millivolt). C CONSTS(141) is kVCa in component I_Ca (per_millivolt). C CONSTS(142) is VhCa in component I_Ca (millivolt). C ALGBRC(76) is gCa in component I_Ca (picosiemens). C ALGBRC(79) is ECa_PU in component I_Ca (millivolt). C CONSTS(143) is T in component model_parameters (kelvin). C CONSTS(144) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(145) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(146) is CO in component model_parameters (micromolar). C CONSTS(147) is gNa in component I_Na (picosiemens). C ALGBRC(84) is ENa_PU in component I_Na (millivolt). C CONSTS(148) is NO in component model_parameters (micromolar). C STATES(36) is NS1 in component NS1 (micromolar). C CONSTS(149) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(86) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(150) is hNSCC in component model_parameters (dimensionless). C CONSTS(151) is ENSCC in component model_parameters (millivolt). C CONSTS(152) is KNSCC in component model_parameters (micromolar). C CONSTS(153) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(88) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(154) is gPM in component I_PM (femtoampere). C CONSTS(155) is KPM in component I_PM (micromolar). C CONSTS(156) is gNaP in component I_NaP (picosiemens). C CONSTS(157) is KNaP in component I_NaP (micromolar). C CONSTS(158) is hNaP in component I_NaP (dimensionless). C CONSTS(159) is ENaP in component I_NaP (millivolt). C ALGBRC(92) is JSERCA in component JSERCA (flux). C CONSTS(160) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(161) is A2 in component JSERCA (dimensionless). C CONSTS(162) is A4 in component JSERCA (per_micromolar). C CONSTS(163) is A5 in component JSERCA (per_micromolar). C CONSTS(164) is A6 in component JSERCA (per_micromolar2). C STATES(37) is CER in component CER (micromolar). C ALGBRC(96) is JMCU in component JMCU (flux). C CONSTS(165) is VMCU in component JMCU (flux). C CONSTS(166) is KMCU in component JMCU (micromolar). C ALGBRC(94) is epsilon_INH in component JMCU (dimensionless). C CONSTS(167) is KINH in component JMCU (micromolar). C CONSTS(168) is hINH in component JMCU (dimensionless). C STATES(38) is CMT in component CMT (micromolar). C ALGBRC(104) is JIPR in component JIPR (flux). C CONSTS(169) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(170) is k_1 in component JIPR (flux). C CONSTS(171) is k1 in component JIPR (first_order_rate_constant). C CONSTS(172) is k2 in component JIPR (first_order_rate_constant). C CONSTS(173) is r2 in component JIPR (first_order_rate_constant). C CONSTS(174) is r_2 in component JIPR (flux). C CONSTS(175) is r4 in component JIPR (first_order_rate_constant). C CONSTS(176) is R1 in component JIPR (micromolar). C CONSTS(177) is R3 in component JIPR (micromolar). C ALGBRC(98) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(101) is phi_1 in component JIPR (flux). C ALGBRC(105) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(107) is phi3 in component JIPR (first_order_rate_constant). C STATES(39) is H in component JIPR (dimensionless). C CONSTS(178) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(179) is h_beta in component JIPR (dimensionless). C CONSTS(180) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(181) is K_beta in component JIPR (micromolar). C CONSTS(703) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(11) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(40) is zeta in component JIPR (dimensionless). C CONSTS(182) is K_phi3_act in component JIPR (micromolar). C CONSTS(183) is K_phi3_inh in component JIPR (micromolar). C CONSTS(184) is h_phi3_act in component JIPR (dimensionless). C CONSTS(185) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(186) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(187) is P in component model_parameters (micromolar). C ALGBRC(99) is JNCX in component JNCX (flux). C CONSTS(188) is VNCX in component JNCX (flux). C CONSTS(189) is KNCX in component JNCX (micromolar). C CONSTS(717) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(718) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(190) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(191) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(192) is gamma_MT in component model_parameters (dimensionless). C CONSTS(193) is gamma_ER in component model_parameters (dimensionless). C CONSTS(719) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(720) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(721) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(194) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(102) is fm in component CMT (dimensionless). C CONSTS(195) is Km in component CMT (micromolar). C CONSTS(196) is Bm in component CMT (micromolar). C CONSTS(197) is KO in component model_parameters (micromolar). C CONSTS(198) is Ki in component model_parameters (micromolar). C CONSTS(199) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(115) is I_iCa in component PU_membrane (picoampere). C ALGBRC(119) is I_iNa in component PU_membrane (picoampere). C CONSTS(728) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(106) is I_Ca in component I_Ca (picoampere). C ALGBRC(109) is I_Na in component I_Na (picoampere). C ALGBRC(117) is I_NaP in component I_NaP (picoampere). C ALGBRC(111) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(114) is I_PM in component I_PM (picoampere). C ALGBRC(113) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(200) is gCa_ in component I_Ca (picosiemens). C CONSTS(201) is kCa in component I_Ca (per_millivolt). C CONSTS(202) is kVCa in component I_Ca (per_millivolt). C CONSTS(203) is VhCa in component I_Ca (millivolt). C ALGBRC(100) is gCa in component I_Ca (picosiemens). C ALGBRC(103) is ECa_PU in component I_Ca (millivolt). C CONSTS(204) is T in component model_parameters (kelvin). C CONSTS(205) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(206) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(207) is CO in component model_parameters (micromolar). C CONSTS(208) is gNa in component I_Na (picosiemens). C ALGBRC(108) is ENa_PU in component I_Na (millivolt). C CONSTS(209) is NO in component model_parameters (micromolar). C STATES(41) is NS1 in component NS1 (micromolar). C CONSTS(210) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(110) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(211) is hNSCC in component model_parameters (dimensionless). C CONSTS(212) is ENSCC in component model_parameters (millivolt). C CONSTS(213) is KNSCC in component model_parameters (micromolar). C CONSTS(214) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(112) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(215) is gPM in component I_PM (femtoampere). C CONSTS(216) is KPM in component I_PM (micromolar). C CONSTS(217) is gNaP in component I_NaP (picosiemens). C CONSTS(218) is KNaP in component I_NaP (micromolar). C CONSTS(219) is hNaP in component I_NaP (dimensionless). C CONSTS(220) is ENaP in component I_NaP (millivolt). C ALGBRC(116) is JSERCA in component JSERCA (flux). C CONSTS(221) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(222) is A2 in component JSERCA (dimensionless). C CONSTS(223) is A4 in component JSERCA (per_micromolar). C CONSTS(224) is A5 in component JSERCA (per_micromolar). C CONSTS(225) is A6 in component JSERCA (per_micromolar2). C STATES(42) is CER in component CER (micromolar). C ALGBRC(120) is JMCU in component JMCU (flux). C CONSTS(226) is VMCU in component JMCU (flux). C CONSTS(227) is KMCU in component JMCU (micromolar). C ALGBRC(118) is epsilon_INH in component JMCU (dimensionless). C CONSTS(228) is KINH in component JMCU (micromolar). C CONSTS(229) is hINH in component JMCU (dimensionless). C STATES(43) is CMT in component CMT (micromolar). C ALGBRC(128) is JIPR in component JIPR (flux). C CONSTS(230) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(231) is k_1 in component JIPR (flux). C CONSTS(232) is k1 in component JIPR (first_order_rate_constant). C CONSTS(233) is k2 in component JIPR (first_order_rate_constant). C CONSTS(234) is r2 in component JIPR (first_order_rate_constant). C CONSTS(235) is r_2 in component JIPR (flux). C CONSTS(236) is r4 in component JIPR (first_order_rate_constant). C CONSTS(237) is R1 in component JIPR (micromolar). C CONSTS(238) is R3 in component JIPR (micromolar). C ALGBRC(122) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(125) is phi_1 in component JIPR (flux). C ALGBRC(129) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(131) is phi3 in component JIPR (first_order_rate_constant). C STATES(44) is H in component JIPR (dimensionless). C CONSTS(239) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(240) is h_beta in component JIPR (dimensionless). C CONSTS(241) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(242) is K_beta in component JIPR (micromolar). C CONSTS(704) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(12) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(45) is zeta in component JIPR (dimensionless). C CONSTS(243) is K_phi3_act in component JIPR (micromolar). C CONSTS(244) is K_phi3_inh in component JIPR (micromolar). C CONSTS(245) is h_phi3_act in component JIPR (dimensionless). C CONSTS(246) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(247) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(248) is P in component model_parameters (micromolar). C ALGBRC(123) is JNCX in component JNCX (flux). C CONSTS(249) is VNCX in component JNCX (flux). C CONSTS(250) is KNCX in component JNCX (micromolar). C CONSTS(723) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(724) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(251) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(252) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(253) is gamma_MT in component model_parameters (dimensionless). C CONSTS(254) is gamma_ER in component model_parameters (dimensionless). C CONSTS(725) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(726) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(727) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(255) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(126) is fm in component CMT (dimensionless). C CONSTS(256) is Km in component CMT (micromolar). C CONSTS(257) is Bm in component CMT (micromolar). C CONSTS(258) is KO in component model_parameters (micromolar). C CONSTS(259) is Ki in component model_parameters (micromolar). C CONSTS(260) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(139) is I_iCa in component PU_membrane (picoampere). C ALGBRC(143) is I_iNa in component PU_membrane (picoampere). C CONSTS(734) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(130) is I_Ca in component I_Ca (picoampere). C ALGBRC(133) is I_Na in component I_Na (picoampere). C ALGBRC(141) is I_NaP in component I_NaP (picoampere). C ALGBRC(135) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(138) is I_PM in component I_PM (picoampere). C ALGBRC(137) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(261) is gCa_ in component I_Ca (picosiemens). C CONSTS(262) is kCa in component I_Ca (per_millivolt). C CONSTS(263) is kVCa in component I_Ca (per_millivolt). C CONSTS(264) is VhCa in component I_Ca (millivolt). C ALGBRC(124) is gCa in component I_Ca (picosiemens). C ALGBRC(127) is ECa_PU in component I_Ca (millivolt). C CONSTS(265) is T in component model_parameters (kelvin). C CONSTS(266) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(267) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(268) is CO in component model_parameters (micromolar). C CONSTS(269) is gNa in component I_Na (picosiemens). C ALGBRC(132) is ENa_PU in component I_Na (millivolt). C CONSTS(270) is NO in component model_parameters (micromolar). C STATES(46) is NS1 in component NS1 (micromolar). C CONSTS(271) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(134) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(272) is hNSCC in component model_parameters (dimensionless). C CONSTS(273) is ENSCC in component model_parameters (millivolt). C CONSTS(274) is KNSCC in component model_parameters (micromolar). C CONSTS(275) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(136) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(276) is gPM in component I_PM (femtoampere). C CONSTS(277) is KPM in component I_PM (micromolar). C CONSTS(278) is gNaP in component I_NaP (picosiemens). C CONSTS(279) is KNaP in component I_NaP (micromolar). C CONSTS(280) is hNaP in component I_NaP (dimensionless). C CONSTS(281) is ENaP in component I_NaP (millivolt). C ALGBRC(140) is JSERCA in component JSERCA (flux). C CONSTS(282) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(283) is A2 in component JSERCA (dimensionless). C CONSTS(284) is A4 in component JSERCA (per_micromolar). C CONSTS(285) is A5 in component JSERCA (per_micromolar). C CONSTS(286) is A6 in component JSERCA (per_micromolar2). C STATES(47) is CER in component CER (micromolar). C ALGBRC(144) is JMCU in component JMCU (flux). C CONSTS(287) is VMCU in component JMCU (flux). C CONSTS(288) is KMCU in component JMCU (micromolar). C ALGBRC(142) is epsilon_INH in component JMCU (dimensionless). C CONSTS(289) is KINH in component JMCU (micromolar). C CONSTS(290) is hINH in component JMCU (dimensionless). C STATES(48) is CMT in component CMT (micromolar). C ALGBRC(152) is JIPR in component JIPR (flux). C CONSTS(291) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(292) is k_1 in component JIPR (flux). C CONSTS(293) is k1 in component JIPR (first_order_rate_constant). C CONSTS(294) is k2 in component JIPR (first_order_rate_constant). C CONSTS(295) is r2 in component JIPR (first_order_rate_constant). C CONSTS(296) is r_2 in component JIPR (flux). C CONSTS(297) is r4 in component JIPR (first_order_rate_constant). C CONSTS(298) is R1 in component JIPR (micromolar). C CONSTS(299) is R3 in component JIPR (micromolar). C ALGBRC(146) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(149) is phi_1 in component JIPR (flux). C ALGBRC(153) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(155) is phi3 in component JIPR (first_order_rate_constant). C STATES(49) is H in component JIPR (dimensionless). C CONSTS(300) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(301) is h_beta in component JIPR (dimensionless). C CONSTS(302) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(303) is K_beta in component JIPR (micromolar). C CONSTS(705) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(13) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(50) is zeta in component JIPR (dimensionless). C CONSTS(304) is K_phi3_act in component JIPR (micromolar). C CONSTS(305) is K_phi3_inh in component JIPR (micromolar). C CONSTS(306) is h_phi3_act in component JIPR (dimensionless). C CONSTS(307) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(308) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(309) is P in component model_parameters (micromolar). C ALGBRC(147) is JNCX in component JNCX (flux). C CONSTS(310) is VNCX in component JNCX (flux). C CONSTS(311) is KNCX in component JNCX (micromolar). C CONSTS(729) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(730) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(312) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(313) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(314) is gamma_MT in component model_parameters (dimensionless). C CONSTS(315) is gamma_ER in component model_parameters (dimensionless). C CONSTS(731) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(732) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(733) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(316) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(150) is fm in component CMT (dimensionless). C CONSTS(317) is Km in component CMT (micromolar). C CONSTS(318) is Bm in component CMT (micromolar). C CONSTS(319) is KO in component model_parameters (micromolar). C CONSTS(320) is Ki in component model_parameters (micromolar). C CONSTS(321) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(163) is I_iCa in component PU_membrane (picoampere). C ALGBRC(167) is I_iNa in component PU_membrane (picoampere). C CONSTS(740) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(154) is I_Ca in component I_Ca (picoampere). C ALGBRC(157) is I_Na in component I_Na (picoampere). C ALGBRC(165) is I_NaP in component I_NaP (picoampere). C ALGBRC(159) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(162) is I_PM in component I_PM (picoampere). C ALGBRC(161) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(322) is gCa_ in component I_Ca (picosiemens). C CONSTS(323) is kCa in component I_Ca (per_millivolt). C CONSTS(324) is kVCa in component I_Ca (per_millivolt). C CONSTS(325) is VhCa in component I_Ca (millivolt). C ALGBRC(148) is gCa in component I_Ca (picosiemens). C ALGBRC(151) is ECa_PU in component I_Ca (millivolt). C CONSTS(326) is T in component model_parameters (kelvin). C CONSTS(327) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(328) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(329) is CO in component model_parameters (micromolar). C CONSTS(330) is gNa in component I_Na (picosiemens). C ALGBRC(156) is ENa_PU in component I_Na (millivolt). C CONSTS(331) is NO in component model_parameters (micromolar). C STATES(51) is NS1 in component NS1 (micromolar). C CONSTS(332) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(158) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(333) is hNSCC in component model_parameters (dimensionless). C CONSTS(334) is ENSCC in component model_parameters (millivolt). C CONSTS(335) is KNSCC in component model_parameters (micromolar). C CONSTS(336) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(160) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(337) is gPM in component I_PM (femtoampere). C CONSTS(338) is KPM in component I_PM (micromolar). C CONSTS(339) is gNaP in component I_NaP (picosiemens). C CONSTS(340) is KNaP in component I_NaP (micromolar). C CONSTS(341) is hNaP in component I_NaP (dimensionless). C CONSTS(342) is ENaP in component I_NaP (millivolt). C ALGBRC(164) is JSERCA in component JSERCA (flux). C CONSTS(343) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(344) is A2 in component JSERCA (dimensionless). C CONSTS(345) is A4 in component JSERCA (per_micromolar). C CONSTS(346) is A5 in component JSERCA (per_micromolar). C CONSTS(347) is A6 in component JSERCA (per_micromolar2). C STATES(52) is CER in component CER (micromolar). C ALGBRC(168) is JMCU in component JMCU (flux). C CONSTS(348) is VMCU in component JMCU (flux). C CONSTS(349) is KMCU in component JMCU (micromolar). C ALGBRC(166) is epsilon_INH in component JMCU (dimensionless). C CONSTS(350) is KINH in component JMCU (micromolar). C CONSTS(351) is hINH in component JMCU (dimensionless). C STATES(53) is CMT in component CMT (micromolar). C ALGBRC(176) is JIPR in component JIPR (flux). C CONSTS(352) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(353) is k_1 in component JIPR (flux). C CONSTS(354) is k1 in component JIPR (first_order_rate_constant). C CONSTS(355) is k2 in component JIPR (first_order_rate_constant). C CONSTS(356) is r2 in component JIPR (first_order_rate_constant). C CONSTS(357) is r_2 in component JIPR (flux). C CONSTS(358) is r4 in component JIPR (first_order_rate_constant). C CONSTS(359) is R1 in component JIPR (micromolar). C CONSTS(360) is R3 in component JIPR (micromolar). C ALGBRC(170) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(173) is phi_1 in component JIPR (flux). C ALGBRC(177) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(179) is phi3 in component JIPR (first_order_rate_constant). C STATES(54) is H in component JIPR (dimensionless). C CONSTS(361) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(362) is h_beta in component JIPR (dimensionless). C CONSTS(363) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(364) is K_beta in component JIPR (micromolar). C CONSTS(706) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(14) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(55) is zeta in component JIPR (dimensionless). C CONSTS(365) is K_phi3_act in component JIPR (micromolar). C CONSTS(366) is K_phi3_inh in component JIPR (micromolar). C CONSTS(367) is h_phi3_act in component JIPR (dimensionless). C CONSTS(368) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(369) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(370) is P in component model_parameters (micromolar). C ALGBRC(171) is JNCX in component JNCX (flux). C CONSTS(371) is VNCX in component JNCX (flux). C CONSTS(372) is KNCX in component JNCX (micromolar). C CONSTS(735) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(736) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(373) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(374) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(375) is gamma_MT in component model_parameters (dimensionless). C CONSTS(376) is gamma_ER in component model_parameters (dimensionless). C CONSTS(737) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(738) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(739) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(377) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(174) is fm in component CMT (dimensionless). C CONSTS(378) is Km in component CMT (micromolar). C CONSTS(379) is Bm in component CMT (micromolar). C CONSTS(380) is KO in component model_parameters (micromolar). C CONSTS(381) is Ki in component model_parameters (micromolar). C CONSTS(382) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(187) is I_iCa in component PU_membrane (picoampere). C ALGBRC(191) is I_iNa in component PU_membrane (picoampere). C CONSTS(746) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(178) is I_Ca in component I_Ca (picoampere). C ALGBRC(181) is I_Na in component I_Na (picoampere). C ALGBRC(189) is I_NaP in component I_NaP (picoampere). C ALGBRC(183) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(186) is I_PM in component I_PM (picoampere). C ALGBRC(185) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(383) is gCa_ in component I_Ca (picosiemens). C CONSTS(384) is kCa in component I_Ca (per_millivolt). C CONSTS(385) is kVCa in component I_Ca (per_millivolt). C CONSTS(386) is VhCa in component I_Ca (millivolt). C ALGBRC(172) is gCa in component I_Ca (picosiemens). C ALGBRC(175) is ECa_PU in component I_Ca (millivolt). C CONSTS(387) is T in component model_parameters (kelvin). C CONSTS(388) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(389) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(390) is CO in component model_parameters (micromolar). C CONSTS(391) is gNa in component I_Na (picosiemens). C ALGBRC(180) is ENa_PU in component I_Na (millivolt). C CONSTS(392) is NO in component model_parameters (micromolar). C STATES(56) is NS1 in component NS1 (micromolar). C CONSTS(393) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(182) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(394) is hNSCC in component model_parameters (dimensionless). C CONSTS(395) is ENSCC in component model_parameters (millivolt). C CONSTS(396) is KNSCC in component model_parameters (micromolar). C CONSTS(397) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(184) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(398) is gPM in component I_PM (femtoampere). C CONSTS(399) is KPM in component I_PM (micromolar). C CONSTS(400) is gNaP in component I_NaP (picosiemens). C CONSTS(401) is KNaP in component I_NaP (micromolar). C CONSTS(402) is hNaP in component I_NaP (dimensionless). C CONSTS(403) is ENaP in component I_NaP (millivolt). C ALGBRC(188) is JSERCA in component JSERCA (flux). C CONSTS(404) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(405) is A2 in component JSERCA (dimensionless). C CONSTS(406) is A4 in component JSERCA (per_micromolar). C CONSTS(407) is A5 in component JSERCA (per_micromolar). C CONSTS(408) is A6 in component JSERCA (per_micromolar2). C STATES(57) is CER in component CER (micromolar). C ALGBRC(192) is JMCU in component JMCU (flux). C CONSTS(409) is VMCU in component JMCU (flux). C CONSTS(410) is KMCU in component JMCU (micromolar). C ALGBRC(190) is epsilon_INH in component JMCU (dimensionless). C CONSTS(411) is KINH in component JMCU (micromolar). C CONSTS(412) is hINH in component JMCU (dimensionless). C STATES(58) is CMT in component CMT (micromolar). C ALGBRC(200) is JIPR in component JIPR (flux). C CONSTS(413) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(414) is k_1 in component JIPR (flux). C CONSTS(415) is k1 in component JIPR (first_order_rate_constant). C CONSTS(416) is k2 in component JIPR (first_order_rate_constant). C CONSTS(417) is r2 in component JIPR (first_order_rate_constant). C CONSTS(418) is r_2 in component JIPR (flux). C CONSTS(419) is r4 in component JIPR (first_order_rate_constant). C CONSTS(420) is R1 in component JIPR (micromolar). C CONSTS(421) is R3 in component JIPR (micromolar). C ALGBRC(194) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(197) is phi_1 in component JIPR (flux). C ALGBRC(201) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(203) is phi3 in component JIPR (first_order_rate_constant). C STATES(59) is H in component JIPR (dimensionless). C CONSTS(422) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(423) is h_beta in component JIPR (dimensionless). C CONSTS(424) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(425) is K_beta in component JIPR (micromolar). C CONSTS(707) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(15) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(60) is zeta in component JIPR (dimensionless). C CONSTS(426) is K_phi3_act in component JIPR (micromolar). C CONSTS(427) is K_phi3_inh in component JIPR (micromolar). C CONSTS(428) is h_phi3_act in component JIPR (dimensionless). C CONSTS(429) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(430) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(431) is P in component model_parameters (micromolar). C ALGBRC(195) is JNCX in component JNCX (flux). C CONSTS(432) is VNCX in component JNCX (flux). C CONSTS(433) is KNCX in component JNCX (micromolar). C CONSTS(741) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(742) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(434) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(435) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(436) is gamma_MT in component model_parameters (dimensionless). C CONSTS(437) is gamma_ER in component model_parameters (dimensionless). C CONSTS(743) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(744) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(745) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(438) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(198) is fm in component CMT (dimensionless). C CONSTS(439) is Km in component CMT (micromolar). C CONSTS(440) is Bm in component CMT (micromolar). C CONSTS(441) is KO in component model_parameters (micromolar). C CONSTS(442) is Ki in component model_parameters (micromolar). C CONSTS(443) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(211) is I_iCa in component PU_membrane (picoampere). C ALGBRC(215) is I_iNa in component PU_membrane (picoampere). C CONSTS(752) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(202) is I_Ca in component I_Ca (picoampere). C ALGBRC(205) is I_Na in component I_Na (picoampere). C ALGBRC(213) is I_NaP in component I_NaP (picoampere). C ALGBRC(207) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(210) is I_PM in component I_PM (picoampere). C ALGBRC(209) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(444) is gCa_ in component I_Ca (picosiemens). C CONSTS(445) is kCa in component I_Ca (per_millivolt). C CONSTS(446) is kVCa in component I_Ca (per_millivolt). C CONSTS(447) is VhCa in component I_Ca (millivolt). C ALGBRC(196) is gCa in component I_Ca (picosiemens). C ALGBRC(199) is ECa_PU in component I_Ca (millivolt). C CONSTS(448) is T in component model_parameters (kelvin). C CONSTS(449) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(450) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(451) is CO in component model_parameters (micromolar). C CONSTS(452) is gNa in component I_Na (picosiemens). C ALGBRC(204) is ENa_PU in component I_Na (millivolt). C CONSTS(453) is NO in component model_parameters (micromolar). C STATES(61) is NS1 in component NS1 (micromolar). C CONSTS(454) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(206) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(455) is hNSCC in component model_parameters (dimensionless). C CONSTS(456) is ENSCC in component model_parameters (millivolt). C CONSTS(457) is KNSCC in component model_parameters (micromolar). C CONSTS(458) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(208) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(459) is gPM in component I_PM (femtoampere). C CONSTS(460) is KPM in component I_PM (micromolar). C CONSTS(461) is gNaP in component I_NaP (picosiemens). C CONSTS(462) is KNaP in component I_NaP (micromolar). C CONSTS(463) is hNaP in component I_NaP (dimensionless). C CONSTS(464) is ENaP in component I_NaP (millivolt). C ALGBRC(212) is JSERCA in component JSERCA (flux). C CONSTS(465) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(466) is A2 in component JSERCA (dimensionless). C CONSTS(467) is A4 in component JSERCA (per_micromolar). C CONSTS(468) is A5 in component JSERCA (per_micromolar). C CONSTS(469) is A6 in component JSERCA (per_micromolar2). C STATES(62) is CER in component CER (micromolar). C ALGBRC(216) is JMCU in component JMCU (flux). C CONSTS(470) is VMCU in component JMCU (flux). C CONSTS(471) is KMCU in component JMCU (micromolar). C ALGBRC(214) is epsilon_INH in component JMCU (dimensionless). C CONSTS(472) is KINH in component JMCU (micromolar). C CONSTS(473) is hINH in component JMCU (dimensionless). C STATES(63) is CMT in component CMT (micromolar). C ALGBRC(224) is JIPR in component JIPR (flux). C CONSTS(474) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(475) is k_1 in component JIPR (flux). C CONSTS(476) is k1 in component JIPR (first_order_rate_constant). C CONSTS(477) is k2 in component JIPR (first_order_rate_constant). C CONSTS(478) is r2 in component JIPR (first_order_rate_constant). C CONSTS(479) is r_2 in component JIPR (flux). C CONSTS(480) is r4 in component JIPR (first_order_rate_constant). C CONSTS(481) is R1 in component JIPR (micromolar). C CONSTS(482) is R3 in component JIPR (micromolar). C ALGBRC(218) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(221) is phi_1 in component JIPR (flux). C ALGBRC(225) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(227) is phi3 in component JIPR (first_order_rate_constant). C STATES(64) is H in component JIPR (dimensionless). C CONSTS(483) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(484) is h_beta in component JIPR (dimensionless). C CONSTS(485) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(486) is K_beta in component JIPR (micromolar). C CONSTS(708) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(16) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(65) is zeta in component JIPR (dimensionless). C CONSTS(487) is K_phi3_act in component JIPR (micromolar). C CONSTS(488) is K_phi3_inh in component JIPR (micromolar). C CONSTS(489) is h_phi3_act in component JIPR (dimensionless). C CONSTS(490) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(491) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(492) is P in component model_parameters (micromolar). C ALGBRC(219) is JNCX in component JNCX (flux). C CONSTS(493) is VNCX in component JNCX (flux). C CONSTS(494) is KNCX in component JNCX (micromolar). C CONSTS(747) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(748) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(495) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(496) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(497) is gamma_MT in component model_parameters (dimensionless). C CONSTS(498) is gamma_ER in component model_parameters (dimensionless). C CONSTS(749) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(750) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(751) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(499) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(222) is fm in component CMT (dimensionless). C CONSTS(500) is Km in component CMT (micromolar). C CONSTS(501) is Bm in component CMT (micromolar). C CONSTS(502) is KO in component model_parameters (micromolar). C CONSTS(503) is Ki in component model_parameters (micromolar). C CONSTS(504) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(235) is I_iCa in component PU_membrane (picoampere). C ALGBRC(239) is I_iNa in component PU_membrane (picoampere). C CONSTS(758) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(226) is I_Ca in component I_Ca (picoampere). C ALGBRC(229) is I_Na in component I_Na (picoampere). C ALGBRC(237) is I_NaP in component I_NaP (picoampere). C ALGBRC(231) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(234) is I_PM in component I_PM (picoampere). C ALGBRC(233) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(505) is gCa_ in component I_Ca (picosiemens). C CONSTS(506) is kCa in component I_Ca (per_millivolt). C CONSTS(507) is kVCa in component I_Ca (per_millivolt). C CONSTS(508) is VhCa in component I_Ca (millivolt). C ALGBRC(220) is gCa in component I_Ca (picosiemens). C ALGBRC(223) is ECa_PU in component I_Ca (millivolt). C CONSTS(509) is T in component model_parameters (kelvin). C CONSTS(510) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(511) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(512) is CO in component model_parameters (micromolar). C CONSTS(513) is gNa in component I_Na (picosiemens). C ALGBRC(228) is ENa_PU in component I_Na (millivolt). C CONSTS(514) is NO in component model_parameters (micromolar). C STATES(66) is NS1 in component NS1 (micromolar). C CONSTS(515) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(230) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(516) is hNSCC in component model_parameters (dimensionless). C CONSTS(517) is ENSCC in component model_parameters (millivolt). C CONSTS(518) is KNSCC in component model_parameters (micromolar). C CONSTS(519) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(232) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(520) is gPM in component I_PM (femtoampere). C CONSTS(521) is KPM in component I_PM (micromolar). C CONSTS(522) is gNaP in component I_NaP (picosiemens). C CONSTS(523) is KNaP in component I_NaP (micromolar). C CONSTS(524) is hNaP in component I_NaP (dimensionless). C CONSTS(525) is ENaP in component I_NaP (millivolt). C ALGBRC(236) is JSERCA in component JSERCA (flux). C CONSTS(526) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(527) is A2 in component JSERCA (dimensionless). C CONSTS(528) is A4 in component JSERCA (per_micromolar). C CONSTS(529) is A5 in component JSERCA (per_micromolar). C CONSTS(530) is A6 in component JSERCA (per_micromolar2). C STATES(67) is CER in component CER (micromolar). C ALGBRC(240) is JMCU in component JMCU (flux). C CONSTS(531) is VMCU in component JMCU (flux). C CONSTS(532) is KMCU in component JMCU (micromolar). C ALGBRC(238) is epsilon_INH in component JMCU (dimensionless). C CONSTS(533) is KINH in component JMCU (micromolar). C CONSTS(534) is hINH in component JMCU (dimensionless). C STATES(68) is CMT in component CMT (micromolar). C ALGBRC(248) is JIPR in component JIPR (flux). C CONSTS(535) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(536) is k_1 in component JIPR (flux). C CONSTS(537) is k1 in component JIPR (first_order_rate_constant). C CONSTS(538) is k2 in component JIPR (first_order_rate_constant). C CONSTS(539) is r2 in component JIPR (first_order_rate_constant). C CONSTS(540) is r_2 in component JIPR (flux). C CONSTS(541) is r4 in component JIPR (first_order_rate_constant). C CONSTS(542) is R1 in component JIPR (micromolar). C CONSTS(543) is R3 in component JIPR (micromolar). C ALGBRC(242) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(245) is phi_1 in component JIPR (flux). C ALGBRC(249) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(251) is phi3 in component JIPR (first_order_rate_constant). C STATES(69) is H in component JIPR (dimensionless). C CONSTS(544) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(545) is h_beta in component JIPR (dimensionless). C CONSTS(546) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(547) is K_beta in component JIPR (micromolar). C CONSTS(709) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(17) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(70) is zeta in component JIPR (dimensionless). C CONSTS(548) is K_phi3_act in component JIPR (micromolar). C CONSTS(549) is K_phi3_inh in component JIPR (micromolar). C CONSTS(550) is h_phi3_act in component JIPR (dimensionless). C CONSTS(551) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(552) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(553) is P in component model_parameters (micromolar). C ALGBRC(243) is JNCX in component JNCX (flux). C CONSTS(554) is VNCX in component JNCX (flux). C CONSTS(555) is KNCX in component JNCX (micromolar). C CONSTS(753) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(754) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(556) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(557) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(558) is gamma_MT in component model_parameters (dimensionless). C CONSTS(559) is gamma_ER in component model_parameters (dimensionless). C CONSTS(755) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(756) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(757) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(560) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(246) is fm in component CMT (dimensionless). C CONSTS(561) is Km in component CMT (micromolar). C CONSTS(562) is Bm in component CMT (micromolar). C CONSTS(563) is KO in component model_parameters (micromolar). C CONSTS(564) is Ki in component model_parameters (micromolar). C CONSTS(565) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(259) is I_iCa in component PU_membrane (picoampere). C ALGBRC(263) is I_iNa in component PU_membrane (picoampere). C CONSTS(764) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(250) is I_Ca in component I_Ca (picoampere). C ALGBRC(253) is I_Na in component I_Na (picoampere). C ALGBRC(261) is I_NaP in component I_NaP (picoampere). C ALGBRC(255) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(258) is I_PM in component I_PM (picoampere). C ALGBRC(257) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(566) is gCa_ in component I_Ca (picosiemens). C CONSTS(567) is kCa in component I_Ca (per_millivolt). C CONSTS(568) is kVCa in component I_Ca (per_millivolt). C CONSTS(569) is VhCa in component I_Ca (millivolt). C ALGBRC(244) is gCa in component I_Ca (picosiemens). C ALGBRC(247) is ECa_PU in component I_Ca (millivolt). C CONSTS(570) is T in component model_parameters (kelvin). C CONSTS(571) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(572) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(573) is CO in component model_parameters (micromolar). C CONSTS(574) is gNa in component I_Na (picosiemens). C ALGBRC(252) is ENa_PU in component I_Na (millivolt). C CONSTS(575) is NO in component model_parameters (micromolar). C STATES(71) is NS1 in component NS1 (micromolar). C CONSTS(576) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(254) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(577) is hNSCC in component model_parameters (dimensionless). C CONSTS(578) is ENSCC in component model_parameters (millivolt). C CONSTS(579) is KNSCC in component model_parameters (micromolar). C CONSTS(580) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(256) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(581) is gPM in component I_PM (femtoampere). C CONSTS(582) is KPM in component I_PM (micromolar). C CONSTS(583) is gNaP in component I_NaP (picosiemens). C CONSTS(584) is KNaP in component I_NaP (micromolar). C CONSTS(585) is hNaP in component I_NaP (dimensionless). C CONSTS(586) is ENaP in component I_NaP (millivolt). C ALGBRC(260) is JSERCA in component JSERCA (flux). C CONSTS(587) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(588) is A2 in component JSERCA (dimensionless). C CONSTS(589) is A4 in component JSERCA (per_micromolar). C CONSTS(590) is A5 in component JSERCA (per_micromolar). C CONSTS(591) is A6 in component JSERCA (per_micromolar2). C STATES(72) is CER in component CER (micromolar). C ALGBRC(264) is JMCU in component JMCU (flux). C CONSTS(592) is VMCU in component JMCU (flux). C CONSTS(593) is KMCU in component JMCU (micromolar). C ALGBRC(262) is epsilon_INH in component JMCU (dimensionless). C CONSTS(594) is KINH in component JMCU (micromolar). C CONSTS(595) is hINH in component JMCU (dimensionless). C STATES(73) is CMT in component CMT (micromolar). C ALGBRC(272) is JIPR in component JIPR (flux). C CONSTS(596) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(597) is k_1 in component JIPR (flux). C CONSTS(598) is k1 in component JIPR (first_order_rate_constant). C CONSTS(599) is k2 in component JIPR (first_order_rate_constant). C CONSTS(600) is r2 in component JIPR (first_order_rate_constant). C CONSTS(601) is r_2 in component JIPR (flux). C CONSTS(602) is r4 in component JIPR (first_order_rate_constant). C CONSTS(603) is R1 in component JIPR (micromolar). C CONSTS(604) is R3 in component JIPR (micromolar). C ALGBRC(266) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(269) is phi_1 in component JIPR (flux). C ALGBRC(273) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(275) is phi3 in component JIPR (first_order_rate_constant). C STATES(74) is H in component JIPR (dimensionless). C CONSTS(605) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(606) is h_beta in component JIPR (dimensionless). C CONSTS(607) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(608) is K_beta in component JIPR (micromolar). C CONSTS(710) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(18) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(75) is zeta in component JIPR (dimensionless). C CONSTS(609) is K_phi3_act in component JIPR (micromolar). C CONSTS(610) is K_phi3_inh in component JIPR (micromolar). C CONSTS(611) is h_phi3_act in component JIPR (dimensionless). C CONSTS(612) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(613) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(614) is P in component model_parameters (micromolar). C ALGBRC(267) is JNCX in component JNCX (flux). C CONSTS(615) is VNCX in component JNCX (flux). C CONSTS(616) is KNCX in component JNCX (micromolar). C CONSTS(759) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(760) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(617) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(618) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(619) is gamma_MT in component model_parameters (dimensionless). C CONSTS(620) is gamma_ER in component model_parameters (dimensionless). C CONSTS(761) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(762) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(763) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(621) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(270) is fm in component CMT (dimensionless). C CONSTS(622) is Km in component CMT (micromolar). C CONSTS(623) is Bm in component CMT (micromolar). C CONSTS(624) is KO in component model_parameters (micromolar). C CONSTS(625) is Ki in component model_parameters (micromolar). C CONSTS(626) is n_PU_base in component model_parameters (dimensionless). C ALGBRC(283) is I_iCa in component PU_membrane (picoampere). C ALGBRC(286) is I_iNa in component PU_membrane (picoampere). C CONSTS(770) is alpha_scale in component model_parameters (dimensionless). C ALGBRC(274) is I_Ca in component I_Ca (picoampere). C ALGBRC(277) is I_Na in component I_Na (picoampere). C ALGBRC(284) is I_NaP in component I_NaP (picoampere). C ALGBRC(279) is I_NSCC_Ca in component I_NSCC_Ca (picoampere). C ALGBRC(282) is I_PM in component I_PM (picoampere). C ALGBRC(281) is I_NSCC_Na in component I_NSCC_Na (picoampere). C CONSTS(627) is gCa_ in component I_Ca (picosiemens). C CONSTS(628) is kCa in component I_Ca (per_millivolt). C CONSTS(629) is kVCa in component I_Ca (per_millivolt). C CONSTS(630) is VhCa in component I_Ca (millivolt). C ALGBRC(268) is gCa in component I_Ca (picosiemens). C ALGBRC(271) is ECa_PU in component I_Ca (millivolt). C CONSTS(631) is T in component model_parameters (kelvin). C CONSTS(632) is R in component model_parameters (attojoule_per_zeptomole_kelvin). C CONSTS(633) is F in component model_parameters (femtocoulomb_per_zeptomole). C CONSTS(634) is CO in component model_parameters (micromolar). C CONSTS(635) is gNa in component I_Na (picosiemens). C ALGBRC(276) is ENa_PU in component I_Na (millivolt). C CONSTS(636) is NO in component model_parameters (micromolar). C STATES(76) is NS1 in component NS1 (micromolar). C CONSTS(637) is gNSCC_Ca_ in component I_NSCC_Ca (picosiemens). C ALGBRC(278) is gNSCC_Ca in component I_NSCC_Ca (picosiemens). C CONSTS(638) is hNSCC in component model_parameters (dimensionless). C CONSTS(639) is ENSCC in component model_parameters (millivolt). C CONSTS(640) is KNSCC in component model_parameters (micromolar). C CONSTS(641) is gNSCC_Na_ in component I_NSCC_Na (picosiemens). C ALGBRC(280) is gNSCC_Na in component I_NSCC_Na (picosiemens). C CONSTS(642) is gPM in component I_PM (femtoampere). C CONSTS(643) is KPM in component I_PM (micromolar). C CONSTS(644) is gNaP in component I_NaP (picosiemens). C CONSTS(645) is KNaP in component I_NaP (micromolar). C CONSTS(646) is hNaP in component I_NaP (dimensionless). C CONSTS(647) is ENaP in component I_NaP (millivolt). C ALGBRC(285) is JSERCA in component JSERCA (flux). C CONSTS(648) is VSERCA in component JSERCA (first_order_rate_constant). C CONSTS(649) is A2 in component JSERCA (dimensionless). C CONSTS(650) is A4 in component JSERCA (per_micromolar). C CONSTS(651) is A5 in component JSERCA (per_micromolar). C CONSTS(652) is A6 in component JSERCA (per_micromolar2). C STATES(77) is CER in component CER (micromolar). C ALGBRC(289) is JMCU in component JMCU (flux). C CONSTS(653) is VMCU in component JMCU (flux). C CONSTS(654) is KMCU in component JMCU (micromolar). C ALGBRC(287) is epsilon_INH in component JMCU (dimensionless). C CONSTS(655) is KINH in component JMCU (micromolar). C CONSTS(656) is hINH in component JMCU (dimensionless). C STATES(78) is CMT in component CMT (micromolar). C ALGBRC(294) is JIPR in component JIPR (flux). C CONSTS(657) is kIPR in component JIPR (first_order_rate_constant). C CONSTS(658) is k_1 in component JIPR (flux). C CONSTS(659) is k1 in component JIPR (first_order_rate_constant). C CONSTS(660) is k2 in component JIPR (first_order_rate_constant). C CONSTS(661) is r2 in component JIPR (first_order_rate_constant). C CONSTS(662) is r_2 in component JIPR (flux). C CONSTS(663) is r4 in component JIPR (first_order_rate_constant). C CONSTS(664) is R1 in component JIPR (micromolar). C CONSTS(665) is R3 in component JIPR (micromolar). C ALGBRC(290) is phi1 in component JIPR (first_order_rate_constant). C ALGBRC(292) is phi_1 in component JIPR (flux). C ALGBRC(295) is phi2 in component JIPR (first_order_rate_constant). C ALGBRC(296) is phi3 in component JIPR (first_order_rate_constant). C STATES(79) is H in component JIPR (dimensionless). C CONSTS(666) is g_beta in component JIPR (first_order_rate_constant). C CONSTS(667) is h_beta in component JIPR (dimensionless). C CONSTS(668) is g_alpha in component JIPR (first_order_rate_constant). C CONSTS(669) is K_beta in component JIPR (micromolar). C CONSTS(711) is alpha_zeta in component JIPR (first_order_rate_constant). C ALGBRC(19) is beta_zeta in component JIPR (first_order_rate_constant). C STATES(80) is zeta in component JIPR (dimensionless). C CONSTS(670) is K_phi3_act in component JIPR (micromolar). C CONSTS(671) is K_phi3_inh in component JIPR (micromolar). C CONSTS(672) is h_phi3_act in component JIPR (dimensionless). C CONSTS(673) is h_phi3_inh in component JIPR (dimensionless). C CONSTS(674) is g_phi3 in component JIPR (first_order_rate_constant). C CONSTS(675) is P in component model_parameters (micromolar). C ALGBRC(291) is JNCX in component JNCX (flux). C CONSTS(676) is VNCX in component JNCX (flux). C CONSTS(677) is KNCX in component JNCX (micromolar). C CONSTS(765) is lambda_MT_S1 in component CS1 (dimensionless). C CONSTS(766) is lambda_ER_S1 in component CS1 (dimensionless). C CONSTS(678) is delta_SPU in component model_parameters (micromolar_per_coulomb). C CONSTS(679) is gamma_S1 in component model_parameters (dimensionless). C CONSTS(680) is gamma_MT in component model_parameters (dimensionless). C CONSTS(681) is gamma_ER in component model_parameters (dimensionless). C CONSTS(767) is lambda_MT_S2 in component CS2 (dimensionless). C CONSTS(768) is lambda_ER_S2 in component CS2 (dimensionless). C CONSTS(769) is lambda_S1_S2 in component CS2 (dimensionless). C CONSTS(682) is gamma_S2 in component model_parameters (dimensionless). C ALGBRC(293) is fm in component CMT (dimensionless). C CONSTS(683) is Km in component CMT (micromolar). C CONSTS(684) is Bm in component CMT (micromolar). C CONSTS(685) is KO in component model_parameters (micromolar). C CONSTS(686) is Ki in component model_parameters (micromolar). C CONSTS(687) is n_PU_base in component model_parameters (dimensionless). C RATES(1) is d/dt OCaT in component I_CaT (dimensionless). C RATES(3) is d/dt dT in component I_CaT_d_gate (dimensionless). C RATES(4) is d/dt fT in component I_CaT_f_gate (dimensionless). C RATES(6) is d/dt dv1_1 in component I_Kv1_1_d_gate (dimensionless). C RATES(7) is d/dt fv1_1 in component I_Kv1_1_f_gate (dimensionless). C RATES(8) is d/dt dERG in component I_KERG_d_gate (dimensionless). C RATES(9) is d/dt dNa in component I_Na_d_gate (dimensionless). C RATES(10) is d/dt fNa in component I_Na_f_gate (dimensionless). C RATES(2) is d/dt Vm in component membrane (millivolt). C RATES(5) is d/dt CCy in component CCy (micromolar). C RATES(34) is d/dt H in component JIPR (dimensionless). C RATES(35) is d/dt zeta in component JIPR (dimensionless). C RATES(21) is d/dt CS1 in component CS1 (micromolar). C RATES(11) is d/dt CS2 in component CS2 (micromolar). C RATES(32) is d/dt CER in component CER (micromolar). C RATES(33) is d/dt CMT in component CMT (micromolar). C RATES(31) is d/dt NS1 in component NS1 (micromolar). C RATES(39) is d/dt H in component JIPR (dimensionless). C RATES(40) is d/dt zeta in component JIPR (dimensionless). C RATES(22) is d/dt CS1 in component CS1 (micromolar). C RATES(12) is d/dt CS2 in component CS2 (micromolar). C RATES(37) is d/dt CER in component CER (micromolar). C RATES(38) is d/dt CMT in component CMT (micromolar). C RATES(36) is d/dt NS1 in component NS1 (micromolar). C RATES(44) is d/dt H in component JIPR (dimensionless). C RATES(45) is d/dt zeta in component JIPR (dimensionless). C RATES(23) is d/dt CS1 in component CS1 (micromolar). C RATES(13) is d/dt CS2 in component CS2 (micromolar). C RATES(42) is d/dt CER in component CER (micromolar). C RATES(43) is d/dt CMT in component CMT (micromolar). C RATES(41) is d/dt NS1 in component NS1 (micromolar). C RATES(49) is d/dt H in component JIPR (dimensionless). C RATES(50) is d/dt zeta in component JIPR (dimensionless). C RATES(24) is d/dt CS1 in component CS1 (micromolar). C RATES(14) is d/dt CS2 in component CS2 (micromolar). C RATES(47) is d/dt CER in component CER (micromolar). C RATES(48) is d/dt CMT in component CMT (micromolar). C RATES(46) is d/dt NS1 in component NS1 (micromolar). C RATES(54) is d/dt H in component JIPR (dimensionless). C RATES(55) is d/dt zeta in component JIPR (dimensionless). C RATES(25) is d/dt CS1 in component CS1 (micromolar). C RATES(15) is d/dt CS2 in component CS2 (micromolar). C RATES(52) is d/dt CER in component CER (micromolar). C RATES(53) is d/dt CMT in component CMT (micromolar). C RATES(51) is d/dt NS1 in component NS1 (micromolar). C RATES(59) is d/dt H in component JIPR (dimensionless). C RATES(60) is d/dt zeta in component JIPR (dimensionless). C RATES(26) is d/dt CS1 in component CS1 (micromolar). C RATES(16) is d/dt CS2 in component CS2 (micromolar). C RATES(57) is d/dt CER in component CER (micromolar). C RATES(58) is d/dt CMT in component CMT (micromolar). C RATES(56) is d/dt NS1 in component NS1 (micromolar). C RATES(64) is d/dt H in component JIPR (dimensionless). C RATES(65) is d/dt zeta in component JIPR (dimensionless). C RATES(27) is d/dt CS1 in component CS1 (micromolar). C RATES(17) is d/dt CS2 in component CS2 (micromolar). C RATES(62) is d/dt CER in component CER (micromolar). C RATES(63) is d/dt CMT in component CMT (micromolar). C RATES(61) is d/dt NS1 in component NS1 (micromolar). C RATES(69) is d/dt H in component JIPR (dimensionless). C RATES(70) is d/dt zeta in component JIPR (dimensionless). C RATES(28) is d/dt CS1 in component CS1 (micromolar). C RATES(18) is d/dt CS2 in component CS2 (micromolar). C RATES(67) is d/dt CER in component CER (micromolar). C RATES(68) is d/dt CMT in component CMT (micromolar). C RATES(66) is d/dt NS1 in component NS1 (micromolar). C RATES(74) is d/dt H in component JIPR (dimensionless). C RATES(75) is d/dt zeta in component JIPR (dimensionless). C RATES(29) is d/dt CS1 in component CS1 (micromolar). C RATES(19) is d/dt CS2 in component CS2 (micromolar). C RATES(72) is d/dt CER in component CER (micromolar). C RATES(73) is d/dt CMT in component CMT (micromolar). C RATES(71) is d/dt NS1 in component NS1 (micromolar). C RATES(79) is d/dt H in component JIPR (dimensionless). C RATES(80) is d/dt zeta in component JIPR (dimensionless). C RATES(30) is d/dt CS1 in component CS1 (micromolar). C RATES(20) is d/dt CS2 in component CS2 (micromolar). C RATES(77) is d/dt CER in component CER (micromolar). C RATES(78) is d/dt CMT in component CMT (micromolar). C RATES(76) is d/dt NS1 in component NS1 (micromolar). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 700.0 CONSTS(2) = 17.0 CONSTS(3) = 800.0 CONSTS(4) = 72 STATES(1) = 0 STATES(2) = -70.0 STATES(3) = 0.010 STATES(4) = 0.001 CONSTS(5) = 0.0025 CONSTS(6) = -0.60 CONSTS(7) = -53 CONSTS(8) = 0.019 CONSTS(9) = 6.75 CONSTS(10) = 2 CONSTS(11) = -40 CONSTS(12) = 1 CONSTS(13) = -65 CONSTS(14) = 100.0 CONSTS(15) = 1.0 STATES(5) = 0.12 CONSTS(16) = 10.0 STATES(6) = 0.000 STATES(7) = 1.0 CONSTS(17) = 0.0 CONSTS(18) = 250.0 CONSTS(19) = -0.077 CONSTS(20) = -4.8 CONSTS(21) = 0.80 CONSTS(22) = -0.13 CONSTS(23) = -25.0 CONSTS(24) = 333.0 CONSTS(25) = 0.10 CONSTS(26) = 0.23 CONSTS(27) = -44.8 CONSTS(28) = 10.0 STATES(8) = 0.000 CONSTS(29) = 1000.0 CONSTS(30) = 0.70 CONSTS(31) = -0.56 CONSTS(32) = 30 CONSTS(33) = 6.0 CONSTS(34) = 2000.0 STATES(9) = 0.000 STATES(10) = 0.100 CONSTS(35) = 250.0 CONSTS(36) = 1000.0 CONSTS(37) = -0.2 CONSTS(38) = -15.0 CONSTS(39) = 1.0 CONSTS(40) = -0.147 CONSTS(41) = -35.3 CONSTS(42) = 10.0 CONSTS(43) = 50.0 CONSTS(44) = -0.036 CONSTS(45) = -38.0 CONSTS(46) = 1.00 CONSTS(47) = 0.167 CONSTS(48) = -83.0 CONSTS(49) = 0.20 CONSTS(50) = 0.0 CONSTS(51) = 20.0 CONSTS(52) = 310.16 CONSTS(53) = 8.314E-3 CONSTS(54) = 0.09649 CONSTS(55) = 1.8E3 CONSTS(56) = 5400.0 CONSTS(57) = 145E3 CONSTS(58) = 140E3 CONSTS(59) = 10E3 CONSTS(60) = 1.0E-3 CONSTS(61) = 0.002 CONSTS(62) = 0.000 CONSTS(63) = 0.3361 CONSTS(64) = 0.3658 CONSTS(65) = 0.4357 CONSTS(66) = 0.5351 CONSTS(67) = 0.6366 CONSTS(68) = 0.7386 CONSTS(69) = 0.8361 CONSTS(70) = 0.9320 CONSTS(71) = 1.0185 STATES(11) = 0.023 STATES(12) = 0.023 STATES(13) = 0.023 STATES(14) = 0.023 STATES(15) = 0.023 STATES(16) = 0.023 STATES(17) = 0.023 STATES(18) = 0.023 STATES(19) = 0.023 STATES(20) = 0.023 CONSTS(72) = 0.30 CONSTS(73) = 0.24 CONSTS(74) = 10.0 STATES(21) = 0.120 STATES(22) = 0.120 STATES(23) = 0.120 STATES(24) = 0.120 STATES(25) = 0.120 STATES(26) = 0.120 STATES(27) = 0.120 STATES(28) = 0.120 STATES(29) = 0.120 STATES(30) = 0.120 CONSTS(75) = 1.30 CONSTS(76) = 0.12 CONSTS(77) = 1.0 CONSTS(78) = 0.074 CONSTS(79) = 0.013 CONSTS(80) = -0.20 CONSTS(81) = -85.0 CONSTS(82) = 310.16 CONSTS(83) = 8.314E-3 CONSTS(84) = 0.09649 CONSTS(85) = 1.8E3 CONSTS(86) = 13.5 CONSTS(87) = 140.0E3 STATES(31) = 1.01E4 CONSTS(88) = 0.10 CONSTS(89) = 3.0 CONSTS(90) = 0.0 CONSTS(91) = 0.12 CONSTS(92) = 160.0 CONSTS(93) = 675.0 CONSTS(94) = 1.0 CONSTS(95) = 187.5 CONSTS(96) = 1.0E4 CONSTS(97) = 4.0 CONSTS(98) = 10.0 CONSTS(99) = 1.0E5 CONSTS(100) = 6E-4 CONSTS(101) = 3.57 CONSTS(102) = 2.7E-5 CONSTS(103) = 2.31E-5 STATES(32) = 200.0 CONSTS(104) = 800.0 CONSTS(105) = 10.0 CONSTS(106) = 10.0 CONSTS(107) = 4.0 STATES(33) = 0.200 CONSTS(108) = 2000.0 CONSTS(109) = 6.4 CONSTS(110) = 0.0 CONSTS(111) = 4.0 CONSTS(112) = 250.0 CONSTS(113) = 0.0 CONSTS(114) = 750.0 CONSTS(115) = 36.0 CONSTS(116) = 300.0 STATES(34) = 0.200 CONSTS(117) = 1500.0 CONSTS(118) = 4.0 CONSTS(119) = 0.85 CONSTS(120) = 0.35 STATES(35) = 0.300 CONSTS(121) = 0.1 CONSTS(122) = 0.5 CONSTS(123) = 3.0 CONSTS(124) = 3.0 CONSTS(125) = 4.5 CONSTS(126) = 1.0 CONSTS(127) = 3.0 CONSTS(128) = 0.3 CONSTS(129) = 18.5 CONSTS(130) = 100.0 CONSTS(131) = 200.0 CONSTS(132) = 20.0 CONSTS(133) = 1.0 CONSTS(134) = 0.01 CONSTS(135) = 100.0 CONSTS(136) = 5400.0 CONSTS(137) = 145E3 CONSTS(138) = 50.0 CONSTS(139) = 0.074 CONSTS(140) = 0.013 CONSTS(141) = -0.20 CONSTS(142) = -85.0 CONSTS(143) = 310.16 CONSTS(144) = 8.314E-3 CONSTS(145) = 0.09649 CONSTS(146) = 1.8E3 CONSTS(147) = 13.5 CONSTS(148) = 140.0E3 STATES(36) = 1.01E4 CONSTS(149) = 0.10 CONSTS(150) = 3.0 CONSTS(151) = 0.0 CONSTS(152) = 0.12 CONSTS(153) = 160.0 CONSTS(154) = 675.0 CONSTS(155) = 1.0 CONSTS(156) = 187.5 CONSTS(157) = 1.0E4 CONSTS(158) = 4.0 CONSTS(159) = 10.0 CONSTS(160) = 1.0E5 CONSTS(161) = 6E-4 CONSTS(162) = 3.57 CONSTS(163) = 2.7E-5 CONSTS(164) = 2.31E-5 STATES(37) = 200.0 CONSTS(165) = 800.0 CONSTS(166) = 10.0 CONSTS(167) = 10.0 CONSTS(168) = 4.0 STATES(38) = 0.200 CONSTS(169) = 2000.0 CONSTS(170) = 6.4 CONSTS(171) = 0.0 CONSTS(172) = 4.0 CONSTS(173) = 250.0 CONSTS(174) = 0.0 CONSTS(175) = 750.0 CONSTS(176) = 36.0 CONSTS(177) = 300.0 STATES(39) = 0.200 CONSTS(178) = 1500.0 CONSTS(179) = 4.0 CONSTS(180) = 0.85 CONSTS(181) = 0.35 STATES(40) = 0.300 CONSTS(182) = 0.1 CONSTS(183) = 0.5 CONSTS(184) = 3.0 CONSTS(185) = 3.0 CONSTS(186) = 4.5 CONSTS(187) = 1.0 CONSTS(188) = 3.0 CONSTS(189) = 0.3 CONSTS(190) = 18.5 CONSTS(191) = 100.0 CONSTS(192) = 200.0 CONSTS(193) = 20.0 CONSTS(194) = 1.0 CONSTS(195) = 0.01 CONSTS(196) = 100.0 CONSTS(197) = 5400.0 CONSTS(198) = 145E3 CONSTS(199) = 50.0 CONSTS(200) = 0.074 CONSTS(201) = 0.013 CONSTS(202) = -0.20 CONSTS(203) = -85.0 CONSTS(204) = 310.16 CONSTS(205) = 8.314E-3 CONSTS(206) = 0.09649 CONSTS(207) = 1.8E3 CONSTS(208) = 13.5 CONSTS(209) = 140.0E3 STATES(41) = 1.01E4 CONSTS(210) = 0.10 CONSTS(211) = 3.0 CONSTS(212) = 0.0 CONSTS(213) = 0.12 CONSTS(214) = 160.0 CONSTS(215) = 675.0 CONSTS(216) = 1.0 CONSTS(217) = 187.5 CONSTS(218) = 1.0E4 CONSTS(219) = 4.0 CONSTS(220) = 10.0 CONSTS(221) = 1.0E5 CONSTS(222) = 6E-4 CONSTS(223) = 3.57 CONSTS(224) = 2.7E-5 CONSTS(225) = 2.31E-5 STATES(42) = 200.0 CONSTS(226) = 800.0 CONSTS(227) = 10.0 CONSTS(228) = 10.0 CONSTS(229) = 4.0 STATES(43) = 0.200 CONSTS(230) = 2000.0 CONSTS(231) = 6.4 CONSTS(232) = 0.0 CONSTS(233) = 4.0 CONSTS(234) = 250.0 CONSTS(235) = 0.0 CONSTS(236) = 750.0 CONSTS(237) = 36.0 CONSTS(238) = 300.0 STATES(44) = 0.200 CONSTS(239) = 1500.0 CONSTS(240) = 4.0 CONSTS(241) = 0.85 CONSTS(242) = 0.35 STATES(45) = 0.300 CONSTS(243) = 0.1 CONSTS(244) = 0.5 CONSTS(245) = 3.0 CONSTS(246) = 3.0 CONSTS(247) = 4.5 CONSTS(248) = 1.0 CONSTS(249) = 3.0 CONSTS(250) = 0.3 CONSTS(251) = 18.5 CONSTS(252) = 100.0 CONSTS(253) = 200.0 CONSTS(254) = 20.0 CONSTS(255) = 1.0 CONSTS(256) = 0.01 CONSTS(257) = 100.0 CONSTS(258) = 5400.0 CONSTS(259) = 145E3 CONSTS(260) = 50.0 CONSTS(261) = 0.074 CONSTS(262) = 0.013 CONSTS(263) = -0.20 CONSTS(264) = -85.0 CONSTS(265) = 310.16 CONSTS(266) = 8.314E-3 CONSTS(267) = 0.09649 CONSTS(268) = 1.8E3 CONSTS(269) = 13.5 CONSTS(270) = 140.0E3 STATES(46) = 1.01E4 CONSTS(271) = 0.10 CONSTS(272) = 3.0 CONSTS(273) = 0.0 CONSTS(274) = 0.12 CONSTS(275) = 160.0 CONSTS(276) = 675.0 CONSTS(277) = 1.0 CONSTS(278) = 187.5 CONSTS(279) = 1.0E4 CONSTS(280) = 4.0 CONSTS(281) = 10.0 CONSTS(282) = 1.0E5 CONSTS(283) = 6E-4 CONSTS(284) = 3.57 CONSTS(285) = 2.7E-5 CONSTS(286) = 2.31E-5 STATES(47) = 200.0 CONSTS(287) = 800.0 CONSTS(288) = 10.0 CONSTS(289) = 10.0 CONSTS(290) = 4.0 STATES(48) = 0.200 CONSTS(291) = 2000.0 CONSTS(292) = 6.4 CONSTS(293) = 0.0 CONSTS(294) = 4.0 CONSTS(295) = 250.0 CONSTS(296) = 0.0 CONSTS(297) = 750.0 CONSTS(298) = 36.0 CONSTS(299) = 300.0 STATES(49) = 0.200 CONSTS(300) = 1500.0 CONSTS(301) = 4.0 CONSTS(302) = 0.85 CONSTS(303) = 0.35 STATES(50) = 0.300 CONSTS(304) = 0.1 CONSTS(305) = 0.5 CONSTS(306) = 3.0 CONSTS(307) = 3.0 CONSTS(308) = 4.5 CONSTS(309) = 1.0 CONSTS(310) = 3.0 CONSTS(311) = 0.3 CONSTS(312) = 18.5 CONSTS(313) = 100.0 CONSTS(314) = 200.0 CONSTS(315) = 20.0 CONSTS(316) = 1.0 CONSTS(317) = 0.01 CONSTS(318) = 100.0 CONSTS(319) = 5400.0 CONSTS(320) = 145E3 CONSTS(321) = 50.0 CONSTS(322) = 0.074 CONSTS(323) = 0.013 CONSTS(324) = -0.20 CONSTS(325) = -85.0 CONSTS(326) = 310.16 CONSTS(327) = 8.314E-3 CONSTS(328) = 0.09649 CONSTS(329) = 1.8E3 CONSTS(330) = 13.5 CONSTS(331) = 140.0E3 STATES(51) = 1.01E4 CONSTS(332) = 0.10 CONSTS(333) = 3.0 CONSTS(334) = 0.0 CONSTS(335) = 0.12 CONSTS(336) = 160.0 CONSTS(337) = 675.0 CONSTS(338) = 1.0 CONSTS(339) = 187.5 CONSTS(340) = 1.0E4 CONSTS(341) = 4.0 CONSTS(342) = 10.0 CONSTS(343) = 1.0E5 CONSTS(344) = 6E-4 CONSTS(345) = 3.57 CONSTS(346) = 2.7E-5 CONSTS(347) = 2.31E-5 STATES(52) = 200.0 CONSTS(348) = 800.0 CONSTS(349) = 10.0 CONSTS(350) = 10.0 CONSTS(351) = 4.0 STATES(53) = 0.200 CONSTS(352) = 2000.0 CONSTS(353) = 6.4 CONSTS(354) = 0.0 CONSTS(355) = 4.0 CONSTS(356) = 250.0 CONSTS(357) = 0.0 CONSTS(358) = 750.0 CONSTS(359) = 36.0 CONSTS(360) = 300.0 STATES(54) = 0.200 CONSTS(361) = 1500.0 CONSTS(362) = 4.0 CONSTS(363) = 0.85 CONSTS(364) = 0.35 STATES(55) = 0.300 CONSTS(365) = 0.1 CONSTS(366) = 0.5 CONSTS(367) = 3.0 CONSTS(368) = 3.0 CONSTS(369) = 4.5 CONSTS(370) = 1.0 CONSTS(371) = 3.0 CONSTS(372) = 0.3 CONSTS(373) = 18.5 CONSTS(374) = 100.0 CONSTS(375) = 200.0 CONSTS(376) = 20.0 CONSTS(377) = 1.0 CONSTS(378) = 0.01 CONSTS(379) = 100.0 CONSTS(380) = 5400.0 CONSTS(381) = 145E3 CONSTS(382) = 50.0 CONSTS(383) = 0.074 CONSTS(384) = 0.013 CONSTS(385) = -0.20 CONSTS(386) = -85.0 CONSTS(387) = 310.16 CONSTS(388) = 8.314E-3 CONSTS(389) = 0.09649 CONSTS(390) = 1.8E3 CONSTS(391) = 13.5 CONSTS(392) = 140.0E3 STATES(56) = 1.01E4 CONSTS(393) = 0.10 CONSTS(394) = 3.0 CONSTS(395) = 0.0 CONSTS(396) = 0.12 CONSTS(397) = 160.0 CONSTS(398) = 675.0 CONSTS(399) = 1.0 CONSTS(400) = 187.5 CONSTS(401) = 1.0E4 CONSTS(402) = 4.0 CONSTS(403) = 10.0 CONSTS(404) = 1.0E5 CONSTS(405) = 6E-4 CONSTS(406) = 3.57 CONSTS(407) = 2.7E-5 CONSTS(408) = 2.31E-5 STATES(57) = 200.0 CONSTS(409) = 800.0 CONSTS(410) = 10.0 CONSTS(411) = 10.0 CONSTS(412) = 4.0 STATES(58) = 0.200 CONSTS(413) = 2000.0 CONSTS(414) = 6.4 CONSTS(415) = 0.0 CONSTS(416) = 4.0 CONSTS(417) = 250.0 CONSTS(418) = 0.0 CONSTS(419) = 750.0 CONSTS(420) = 36.0 CONSTS(421) = 300.0 STATES(59) = 0.200 CONSTS(422) = 1500.0 CONSTS(423) = 4.0 CONSTS(424) = 0.85 CONSTS(425) = 0.35 STATES(60) = 0.300 CONSTS(426) = 0.1 CONSTS(427) = 0.5 CONSTS(428) = 3.0 CONSTS(429) = 3.0 CONSTS(430) = 4.5 CONSTS(431) = 1.0 CONSTS(432) = 3.0 CONSTS(433) = 0.3 CONSTS(434) = 18.5 CONSTS(435) = 100.0 CONSTS(436) = 200.0 CONSTS(437) = 20.0 CONSTS(438) = 1.0 CONSTS(439) = 0.01 CONSTS(440) = 100.0 CONSTS(441) = 5400.0 CONSTS(442) = 145E3 CONSTS(443) = 50.0 CONSTS(444) = 0.074 CONSTS(445) = 0.013 CONSTS(446) = -0.20 CONSTS(447) = -85.0 CONSTS(448) = 310.16 CONSTS(449) = 8.314E-3 CONSTS(450) = 0.09649 CONSTS(451) = 1.8E3 CONSTS(452) = 13.5 CONSTS(453) = 140.0E3 STATES(61) = 1.01E4 CONSTS(454) = 0.10 CONSTS(455) = 3.0 CONSTS(456) = 0.0 CONSTS(457) = 0.12 CONSTS(458) = 160.0 CONSTS(459) = 675.0 CONSTS(460) = 1.0 CONSTS(461) = 187.5 CONSTS(462) = 1.0E4 CONSTS(463) = 4.0 CONSTS(464) = 10.0 CONSTS(465) = 1.0E5 CONSTS(466) = 6E-4 CONSTS(467) = 3.57 CONSTS(468) = 2.7E-5 CONSTS(469) = 2.31E-5 STATES(62) = 200.0 CONSTS(470) = 800.0 CONSTS(471) = 10.0 CONSTS(472) = 10.0 CONSTS(473) = 4.0 STATES(63) = 0.200 CONSTS(474) = 2000.0 CONSTS(475) = 6.4 CONSTS(476) = 0.0 CONSTS(477) = 4.0 CONSTS(478) = 250.0 CONSTS(479) = 0.0 CONSTS(480) = 750.0 CONSTS(481) = 36.0 CONSTS(482) = 300.0 STATES(64) = 0.200 CONSTS(483) = 1500.0 CONSTS(484) = 4.0 CONSTS(485) = 0.85 CONSTS(486) = 0.35 STATES(65) = 0.300 CONSTS(487) = 0.1 CONSTS(488) = 0.5 CONSTS(489) = 3.0 CONSTS(490) = 3.0 CONSTS(491) = 4.5 CONSTS(492) = 1.0 CONSTS(493) = 3.0 CONSTS(494) = 0.3 CONSTS(495) = 18.5 CONSTS(496) = 100.0 CONSTS(497) = 200.0 CONSTS(498) = 20.0 CONSTS(499) = 1.0 CONSTS(500) = 0.01 CONSTS(501) = 100.0 CONSTS(502) = 5400.0 CONSTS(503) = 145E3 CONSTS(504) = 50.0 CONSTS(505) = 0.074 CONSTS(506) = 0.013 CONSTS(507) = -0.20 CONSTS(508) = -85.0 CONSTS(509) = 310.16 CONSTS(510) = 8.314E-3 CONSTS(511) = 0.09649 CONSTS(512) = 1.8E3 CONSTS(513) = 13.5 CONSTS(514) = 140.0E3 STATES(66) = 1.01E4 CONSTS(515) = 0.10 CONSTS(516) = 3.0 CONSTS(517) = 0.0 CONSTS(518) = 0.12 CONSTS(519) = 160.0 CONSTS(520) = 675.0 CONSTS(521) = 1.0 CONSTS(522) = 187.5 CONSTS(523) = 1.0E4 CONSTS(524) = 4.0 CONSTS(525) = 10.0 CONSTS(526) = 1.0E5 CONSTS(527) = 6E-4 CONSTS(528) = 3.57 CONSTS(529) = 2.7E-5 CONSTS(530) = 2.31E-5 STATES(67) = 200.0 CONSTS(531) = 800.0 CONSTS(532) = 10.0 CONSTS(533) = 10.0 CONSTS(534) = 4.0 STATES(68) = 0.200 CONSTS(535) = 2000.0 CONSTS(536) = 6.4 CONSTS(537) = 0.0 CONSTS(538) = 4.0 CONSTS(539) = 250.0 CONSTS(540) = 0.0 CONSTS(541) = 750.0 CONSTS(542) = 36.0 CONSTS(543) = 300.0 STATES(69) = 0.200 CONSTS(544) = 1500.0 CONSTS(545) = 4.0 CONSTS(546) = 0.85 CONSTS(547) = 0.35 STATES(70) = 0.300 CONSTS(548) = 0.1 CONSTS(549) = 0.5 CONSTS(550) = 3.0 CONSTS(551) = 3.0 CONSTS(552) = 4.5 CONSTS(553) = 1.0 CONSTS(554) = 3.0 CONSTS(555) = 0.3 CONSTS(556) = 18.5 CONSTS(557) = 100.0 CONSTS(558) = 200.0 CONSTS(559) = 20.0 CONSTS(560) = 1.0 CONSTS(561) = 0.01 CONSTS(562) = 100.0 CONSTS(563) = 5400.0 CONSTS(564) = 145E3 CONSTS(565) = 50.0 CONSTS(566) = 0.074 CONSTS(567) = 0.013 CONSTS(568) = -0.20 CONSTS(569) = -85.0 CONSTS(570) = 310.16 CONSTS(571) = 8.314E-3 CONSTS(572) = 0.09649 CONSTS(573) = 1.8E3 CONSTS(574) = 13.5 CONSTS(575) = 140.0E3 STATES(71) = 1.01E4 CONSTS(576) = 0.10 CONSTS(577) = 3.0 CONSTS(578) = 0.0 CONSTS(579) = 0.12 CONSTS(580) = 160.0 CONSTS(581) = 675.0 CONSTS(582) = 1.0 CONSTS(583) = 187.5 CONSTS(584) = 1.0E4 CONSTS(585) = 4.0 CONSTS(586) = 10.0 CONSTS(587) = 1.0E5 CONSTS(588) = 6E-4 CONSTS(589) = 3.57 CONSTS(590) = 2.7E-5 CONSTS(591) = 2.31E-5 STATES(72) = 200.0 CONSTS(592) = 800.0 CONSTS(593) = 10.0 CONSTS(594) = 10.0 CONSTS(595) = 4.0 STATES(73) = 0.200 CONSTS(596) = 2000.0 CONSTS(597) = 6.4 CONSTS(598) = 0.0 CONSTS(599) = 4.0 CONSTS(600) = 250.0 CONSTS(601) = 0.0 CONSTS(602) = 750.0 CONSTS(603) = 36.0 CONSTS(604) = 300.0 STATES(74) = 0.200 CONSTS(605) = 1500.0 CONSTS(606) = 4.0 CONSTS(607) = 0.85 CONSTS(608) = 0.35 STATES(75) = 0.300 CONSTS(609) = 0.1 CONSTS(610) = 0.5 CONSTS(611) = 3.0 CONSTS(612) = 3.0 CONSTS(613) = 4.5 CONSTS(614) = 1.0 CONSTS(615) = 3.0 CONSTS(616) = 0.3 CONSTS(617) = 18.5 CONSTS(618) = 100.0 CONSTS(619) = 200.0 CONSTS(620) = 20.0 CONSTS(621) = 1.0 CONSTS(622) = 0.01 CONSTS(623) = 100.0 CONSTS(624) = 5400.0 CONSTS(625) = 145E3 CONSTS(626) = 50.0 CONSTS(627) = 0.074 CONSTS(628) = 0.013 CONSTS(629) = -0.20 CONSTS(630) = -85.0 CONSTS(631) = 310.16 CONSTS(632) = 8.314E-3 CONSTS(633) = 0.09649 CONSTS(634) = 1.8E3 CONSTS(635) = 13.5 CONSTS(636) = 140.0E3 STATES(76) = 1.01E4 CONSTS(637) = 0.10 CONSTS(638) = 3.0 CONSTS(639) = 0.0 CONSTS(640) = 0.12 CONSTS(641) = 160.0 CONSTS(642) = 675.0 CONSTS(643) = 1.0 CONSTS(644) = 187.5 CONSTS(645) = 1.0E4 CONSTS(646) = 4.0 CONSTS(647) = 10.0 CONSTS(648) = 1.0E5 CONSTS(649) = 6E-4 CONSTS(650) = 3.57 CONSTS(651) = 2.7E-5 CONSTS(652) = 2.31E-5 STATES(77) = 200.0 CONSTS(653) = 800.0 CONSTS(654) = 10.0 CONSTS(655) = 10.0 CONSTS(656) = 4.0 STATES(78) = 0.200 CONSTS(657) = 2000.0 CONSTS(658) = 6.4 CONSTS(659) = 0.0 CONSTS(660) = 4.0 CONSTS(661) = 250.0 CONSTS(662) = 0.0 CONSTS(663) = 750.0 CONSTS(664) = 36.0 CONSTS(665) = 300.0 STATES(79) = 0.200 CONSTS(666) = 1500.0 CONSTS(667) = 4.0 CONSTS(668) = 0.85 CONSTS(669) = 0.35 STATES(80) = 0.300 CONSTS(670) = 0.1 CONSTS(671) = 0.5 CONSTS(672) = 3.0 CONSTS(673) = 3.0 CONSTS(674) = 4.5 CONSTS(675) = 1.0 CONSTS(676) = 3.0 CONSTS(677) = 0.3 CONSTS(678) = 18.5 CONSTS(679) = 100.0 CONSTS(680) = 200.0 CONSTS(681) = 20.0 CONSTS(682) = 1.0 CONSTS(683) = 0.01 CONSTS(684) = 100.0 CONSTS(685) = 5400.0 CONSTS(686) = 145E3 CONSTS(687) = 50.0 CONSTS(688) = CONSTS(5) CONSTS(689) = (( CONSTS(53)*CONSTS(52))/CONSTS(54))*log(CONSTS(56)/CONSTS(57)) CONSTS(690) = (( CONSTS(53)*CONSTS(52))/CONSTS(54))*log(CONSTS(58)/CONSTS(59)) CONSTS(691) = CONSTS(131)/CONSTS(130) CONSTS(692) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((1.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(693) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((2.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(694) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((3.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(695) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((4.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(696) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((5.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(697) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((6.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(698) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((7.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(699) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((8.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(700) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((9.00000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(701) = CONSTS(72)+ (CONSTS(73) - CONSTS(72))*((10.0000 - 1.00000)/(CONSTS(74) - 1.00000)) CONSTS(702) = CONSTS(119) CONSTS(703) = CONSTS(180) CONSTS(704) = CONSTS(241) CONSTS(705) = CONSTS(302) CONSTS(706) = CONSTS(363) CONSTS(707) = CONSTS(424) CONSTS(708) = CONSTS(485) CONSTS(709) = CONSTS(546) CONSTS(710) = CONSTS(607) CONSTS(711) = CONSTS(668) CONSTS(712) = CONSTS(132)/CONSTS(130) CONSTS(713) = CONSTS(131)/CONSTS(133) CONSTS(714) = CONSTS(132)/CONSTS(133) CONSTS(715) = CONSTS(130)/CONSTS(133) CONSTS(716) = CONSTS(138)/CONSTS(74) CONSTS(717) = CONSTS(192)/CONSTS(191) CONSTS(718) = CONSTS(193)/CONSTS(191) CONSTS(719) = CONSTS(192)/CONSTS(194) CONSTS(720) = CONSTS(193)/CONSTS(194) CONSTS(721) = CONSTS(191)/CONSTS(194) CONSTS(722) = CONSTS(199)/CONSTS(74) CONSTS(723) = CONSTS(253)/CONSTS(252) CONSTS(724) = CONSTS(254)/CONSTS(252) CONSTS(725) = CONSTS(253)/CONSTS(255) CONSTS(726) = CONSTS(254)/CONSTS(255) CONSTS(727) = CONSTS(252)/CONSTS(255) CONSTS(728) = CONSTS(260)/CONSTS(74) CONSTS(729) = CONSTS(314)/CONSTS(313) CONSTS(730) = CONSTS(315)/CONSTS(313) CONSTS(731) = CONSTS(314)/CONSTS(316) CONSTS(732) = CONSTS(315)/CONSTS(316) CONSTS(733) = CONSTS(313)/CONSTS(316) CONSTS(734) = CONSTS(321)/CONSTS(74) CONSTS(735) = CONSTS(375)/CONSTS(374) CONSTS(736) = CONSTS(376)/CONSTS(374) CONSTS(737) = CONSTS(375)/CONSTS(377) CONSTS(738) = CONSTS(376)/CONSTS(377) CONSTS(739) = CONSTS(374)/CONSTS(377) CONSTS(740) = CONSTS(382)/CONSTS(74) CONSTS(741) = CONSTS(436)/CONSTS(435) CONSTS(742) = CONSTS(437)/CONSTS(435) CONSTS(743) = CONSTS(436)/CONSTS(438) CONSTS(744) = CONSTS(437)/CONSTS(438) CONSTS(745) = CONSTS(435)/CONSTS(438) CONSTS(746) = CONSTS(443)/CONSTS(74) CONSTS(747) = CONSTS(497)/CONSTS(496) CONSTS(748) = CONSTS(498)/CONSTS(496) CONSTS(749) = CONSTS(497)/CONSTS(499) CONSTS(750) = CONSTS(498)/CONSTS(499) CONSTS(751) = CONSTS(496)/CONSTS(499) CONSTS(752) = CONSTS(504)/CONSTS(74) CONSTS(753) = CONSTS(558)/CONSTS(557) CONSTS(754) = CONSTS(559)/CONSTS(557) CONSTS(755) = CONSTS(558)/CONSTS(560) CONSTS(756) = CONSTS(559)/CONSTS(560) CONSTS(757) = CONSTS(557)/CONSTS(560) CONSTS(758) = CONSTS(565)/CONSTS(74) CONSTS(759) = CONSTS(619)/CONSTS(618) CONSTS(760) = CONSTS(620)/CONSTS(618) CONSTS(761) = CONSTS(619)/CONSTS(621) CONSTS(762) = CONSTS(620)/CONSTS(621) CONSTS(763) = CONSTS(618)/CONSTS(621) CONSTS(764) = CONSTS(626)/CONSTS(74) CONSTS(765) = CONSTS(680)/CONSTS(679) CONSTS(766) = CONSTS(681)/CONSTS(679) CONSTS(767) = CONSTS(680)/CONSTS(682) CONSTS(768) = CONSTS(681)/CONSTS(682) CONSTS(769) = CONSTS(679)/CONSTS(682) CONSTS(770) = CONSTS(687)/CONSTS(74) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = CONSTS(3)*STATES(3)*STATES(4)*(1.00000 - STATES(1)) - CONSTS(4)*STATES(1) ALGBRC(1) = 1.00000/(1.00000+EXP( CONSTS(6)*(STATES(2) - CONSTS(7)))) RATES(3) = (ALGBRC(1) - STATES(3))/CONSTS(688) ALGBRC(10) = CONSTS(117)/(1.00000+CONSTS(120)/STATES(11) ** CONSTS(118)) RATES(35) = CONSTS(702)*(1.00000 - STATES(35)) - ALGBRC(10)*STATES(35) ALGBRC(11) = CONSTS(178)/(1.00000+CONSTS(181)/STATES(12) ** CONSTS(179)) RATES(40) = CONSTS(703)*(1.00000 - STATES(40)) - ALGBRC(11)*STATES(40) ALGBRC(12) = CONSTS(239)/(1.00000+CONSTS(242)/STATES(13) ** CONSTS(240)) RATES(45) = CONSTS(704)*(1.00000 - STATES(45)) - ALGBRC(12)*STATES(45) ALGBRC(13) = CONSTS(300)/(1.00000+CONSTS(303)/STATES(14) ** CONSTS(301)) RATES(50) = CONSTS(705)*(1.00000 - STATES(50)) - ALGBRC(13)*STATES(50) ALGBRC(14) = CONSTS(361)/(1.00000+CONSTS(364)/STATES(15) ** CONSTS(362)) RATES(55) = CONSTS(706)*(1.00000 - STATES(55)) - ALGBRC(14)*STATES(55) ALGBRC(15) = CONSTS(422)/(1.00000+CONSTS(425)/STATES(16) ** CONSTS(423)) RATES(60) = CONSTS(707)*(1.00000 - STATES(60)) - ALGBRC(15)*STATES(60) ALGBRC(16) = CONSTS(483)/(1.00000+CONSTS(486)/STATES(17) ** CONSTS(484)) RATES(65) = CONSTS(708)*(1.00000 - STATES(65)) - ALGBRC(16)*STATES(65) ALGBRC(17) = CONSTS(544)/(1.00000+CONSTS(547)/STATES(18) ** CONSTS(545)) RATES(70) = CONSTS(709)*(1.00000 - STATES(70)) - ALGBRC(17)*STATES(70) ALGBRC(18) = CONSTS(605)/(1.00000+CONSTS(608)/STATES(19) ** CONSTS(606)) RATES(75) = CONSTS(710)*(1.00000 - STATES(75)) - ALGBRC(18)*STATES(75) ALGBRC(19) = CONSTS(666)/(1.00000+CONSTS(669)/STATES(20) ** CONSTS(667)) RATES(80) = CONSTS(711)*(1.00000 - STATES(80)) - ALGBRC(19)*STATES(80) ALGBRC(2) = 1.00000/(1.00000+EXP( CONSTS(12)*(STATES(2) - CONSTS(13)))) ALGBRC(20) = CONSTS(8)+(CONSTS(9) - CONSTS(8))/(1.00000+EXP( CONSTS(10)*(STATES(2) - CONSTS(11)))) RATES(4) = (ALGBRC(2) - STATES(4))/ALGBRC(20) ALGBRC(4) = CONSTS(24)*(CONSTS(25)/(1.00000+EXP( CONSTS(26)*(STATES(2) - CONSTS(27))))+(1.00000 - CONSTS(25))) ALGBRC(22) = CONSTS(24)*(CONSTS(25) - CONSTS(25)/(1.00000+EXP( CONSTS(26)*(STATES(2) - CONSTS(27))))) RATES(7) = ALGBRC(4)*(1.00000 - STATES(7)) - ALGBRC(22)*STATES(7) ALGBRC(5) = CONSTS(29)*(CONSTS(30)/(1.00000+EXP( CONSTS(31)*(STATES(2) - CONSTS(32))))+(1.00000 - CONSTS(30))) ALGBRC(23) = CONSTS(29)*(CONSTS(30) - CONSTS(30)/(1.00000+EXP( CONSTS(31)*(STATES(2) - CONSTS(32))))) RATES(8) = ALGBRC(5)*(1.00000 - STATES(8)) - ALGBRC(23)*STATES(8) ALGBRC(3) = CONSTS(17)+(CONSTS(18) - CONSTS(17))/(1.00000+EXP( CONSTS(19)*(STATES(2) - CONSTS(20)))) ALGBRC(21) = ALGBRC(3)*(CONSTS(21)/(1.00000+EXP( CONSTS(22)*(STATES(2) - CONSTS(23))))+(1.00000 - CONSTS(21))) ALGBRC(27) = ALGBRC(3)*(CONSTS(21) - CONSTS(21)/(1.00000+EXP( CONSTS(22)*(STATES(2) - CONSTS(23))))) RATES(6) = ALGBRC(21)*(1.00000 - STATES(6)) - ALGBRC(27)*STATES(6) ALGBRC(6) = CONSTS(35)+(CONSTS(36) - CONSTS(35))/(1.00000+EXP( CONSTS(37)*(STATES(2) - CONSTS(38)))) ALGBRC(24) = ALGBRC(6)*(CONSTS(39)/(1.00000+EXP( CONSTS(40)*(STATES(2) - CONSTS(41))))+(1.00000 - CONSTS(39))) ALGBRC(28) = ALGBRC(6)*(CONSTS(39) - CONSTS(39)/(1.00000+EXP( CONSTS(40)*(STATES(2) - CONSTS(41))))) RATES(9) = ALGBRC(24)*(1.00000 - STATES(9)) - ALGBRC(28)*STATES(9) ALGBRC(7) = CONSTS(42)+(CONSTS(43) - CONSTS(42))/(1.00000+EXP( CONSTS(44)*(STATES(2) - CONSTS(45)))) ALGBRC(25) = ALGBRC(7)*(CONSTS(46)/(1.00000+EXP( CONSTS(47)*(STATES(2) - CONSTS(48))))+(1.00000 - CONSTS(46))) ALGBRC(29) = ALGBRC(7)*(CONSTS(46) - CONSTS(46)/(1.00000+EXP( CONSTS(47)*(STATES(2) - CONSTS(48))))) RATES(10) = ALGBRC(25)*(1.00000 - STATES(10)) - ALGBRC(29)*STATES(10) ALGBRC(8) = CONSTS(1)*STATES(1)*(STATES(2) - CONSTS(2)) ALGBRC(26) = CONSTS(14)*(STATES(5)/(CONSTS(15)+STATES(5))) ALGBRC(46) = CONSTS(75)*(CONSTS(76) - STATES(5)) ALGBRC(36) = CONSTS(62)*(STATES(5) - STATES(11)) ALGBRC(37) = CONSTS(63)*(STATES(5) - STATES(12)) ALGBRC(38) = CONSTS(64)*(STATES(5) - STATES(13)) ALGBRC(39) = CONSTS(65)*(STATES(5) - STATES(14)) ALGBRC(40) = CONSTS(66)*(STATES(5) - STATES(15)) ALGBRC(41) = CONSTS(67)*(STATES(5) - STATES(16)) ALGBRC(42) = CONSTS(68)*(STATES(5) - STATES(17)) ALGBRC(43) = CONSTS(69)*(STATES(5) - STATES(18)) ALGBRC(44) = CONSTS(70)*(STATES(5) - STATES(19)) ALGBRC(45) = CONSTS(71)*(STATES(5) - STATES(20)) RATES(5) = ALGBRC(46) - ( CONSTS(60)*ALGBRC(36)+ CONSTS(60)*ALGBRC(37)+ CONSTS(60)*ALGBRC(38)+ CONSTS(60)*ALGBRC(39)+ CONSTS(60)*ALGBRC(40)+ CONSTS(60)*ALGBRC(41)+ CONSTS(60)*ALGBRC(42)+ CONSTS(60)*ALGBRC(43)+ CONSTS(60)*ALGBRC(44)+ CONSTS(60)*ALGBRC(45)+ (CONSTS(61)/2.00000)*(ALGBRC(8)+ALGBRC(26))) ALGBRC(60) = (( CONSTS(83)*CONSTS(82))/CONSTS(84))*log(CONSTS(87)/STATES(31)) ALGBRC(61) = CONSTS(86)*(STATES(2) - ALGBRC(60)) ALGBRC(69) = CONSTS(95)*(STATES(31) ** CONSTS(97)/(CONSTS(96) ** CONSTS(97)+STATES(31) ** CONSTS(97)))*(CONSTS(98) - STATES(2)) ALGBRC(64) = CONSTS(92)*(CONSTS(91) ** CONSTS(89)/(CONSTS(91) ** CONSTS(89)+STATES(21) ** CONSTS(89))) ALGBRC(65) = ALGBRC(64)*(STATES(2) - CONSTS(90)) ALGBRC(71) = CONSTS(716)*(ALGBRC(65)+ALGBRC(61)+ALGBRC(69)) RATES(31) = - (CONSTS(129)/( 1.00000*CONSTS(716)))*ALGBRC(71) ALGBRC(47) = CONSTS(692)*(STATES(11) - STATES(21)) ALGBRC(57) = ( CONSTS(78)*EXP( CONSTS(79)*STATES(2)))/(1.00000+EXP( CONSTS(80)*(STATES(2) - CONSTS(81)))) ALGBRC(58) = (( CONSTS(83)*CONSTS(82))/( 2.00000*CONSTS(84)))*log(CONSTS(85)/STATES(21)) ALGBRC(59) = ALGBRC(57)*(STATES(2) - ALGBRC(58)) ALGBRC(62) = CONSTS(88)*(CONSTS(91) ** CONSTS(89)/(CONSTS(91) ** CONSTS(89)+STATES(21) ** CONSTS(89))) ALGBRC(63) = ALGBRC(62)*(STATES(2) - CONSTS(90)) ALGBRC(66) = CONSTS(93)*(STATES(21) ** 2.00000/(CONSTS(94) ** 2.00000+STATES(21) ** 2.00000)) ALGBRC(67) = CONSTS(716)*(ALGBRC(59)+ALGBRC(63)+ALGBRC(66)) ALGBRC(68) = ( CONSTS(99)*(STATES(21) - CONSTS(100)*STATES(32)))/(1.00000+ CONSTS(101)*STATES(21)+ CONSTS(102)*STATES(32)+ CONSTS(103)*STATES(21)*STATES(32)) ALGBRC(75) = CONSTS(127)*(STATES(33)/(STATES(33)+CONSTS(128))) RATES(21) = (ALGBRC(47)+ CONSTS(691)*ALGBRC(75)) - ( (CONSTS(129)/( CONSTS(716)*2.00000))*ALGBRC(67)+ CONSTS(712)*ALGBRC(68)) ALGBRC(70) = CONSTS(106) ** CONSTS(107)/(CONSTS(106) ** CONSTS(107)+STATES(33) ** CONSTS(107)) ALGBRC(72) = CONSTS(104)*(STATES(11) ** 2.00000/(CONSTS(105) ** 2.00000+STATES(11) ** 2.00000))*ALGBRC(70) ALGBRC(78) = 1.00000/(1.00000+( CONSTS(134)*CONSTS(135))/CONSTS(134)+STATES(33) ** 2.00000) RATES(33) = ALGBRC(78)*(ALGBRC(72) - ALGBRC(75)) ALGBRC(74) = ( CONSTS(110)*CONSTS(115)+ CONSTS(112)*STATES(11))/(CONSTS(115)+STATES(11)) ALGBRC(77) = ( (CONSTS(109)+CONSTS(113))*CONSTS(116))/(CONSTS(116)+STATES(11)) ALGBRC(80) = CONSTS(108)*( CONSTS(126)*ALGBRC(74)*STATES(34))/( CONSTS(126)*ALGBRC(74)+ALGBRC(77)) ** 4.00000*(STATES(32) - STATES(11)) RATES(11) = (ALGBRC(36)+ CONSTS(714)*ALGBRC(80)) - ( CONSTS(715)*ALGBRC(47)+ CONSTS(713)*ALGBRC(72)) RATES(32) = ALGBRC(68) - ALGBRC(80) ALGBRC(81) = ( CONSTS(111)*CONSTS(116)+ CONSTS(114)*STATES(11))/(CONSTS(116)+STATES(11)) ALGBRC(83) = ( CONSTS(125)*STATES(35))/( (1.00000+CONSTS(121)/STATES(11) ** CONSTS(123))*(1.00000+STATES(11)/CONSTS(122) ** CONSTS(124))) RATES(34) = ALGBRC(83)*(1.00000 - STATES(34)) - (( CONSTS(126)*ALGBRC(74)*ALGBRC(81))/( CONSTS(126)*ALGBRC(74)+ALGBRC(77)))*STATES(34) ALGBRC(84) = (( CONSTS(144)*CONSTS(143))/CONSTS(145))*log(CONSTS(148)/STATES(36)) ALGBRC(85) = CONSTS(147)*(STATES(2) - ALGBRC(84)) ALGBRC(93) = CONSTS(156)*(STATES(36) ** CONSTS(158)/(CONSTS(157) ** CONSTS(158)+STATES(36) ** CONSTS(158)))*(CONSTS(159) - STATES(2)) ALGBRC(88) = CONSTS(153)*(CONSTS(152) ** CONSTS(150)/(CONSTS(152) ** CONSTS(150)+STATES(22) ** CONSTS(150))) ALGBRC(89) = ALGBRC(88)*(STATES(2) - CONSTS(151)) ALGBRC(95) = CONSTS(722)*(ALGBRC(89)+ALGBRC(85)+ALGBRC(93)) RATES(36) = - (CONSTS(190)/( 1.00000*CONSTS(722)))*ALGBRC(95) ALGBRC(48) = CONSTS(693)*(STATES(12) - STATES(22)) ALGBRC(76) = ( CONSTS(139)*EXP( CONSTS(140)*STATES(2)))/(1.00000+EXP( CONSTS(141)*(STATES(2) - CONSTS(142)))) ALGBRC(79) = (( CONSTS(144)*CONSTS(143))/( 2.00000*CONSTS(145)))*log(CONSTS(146)/STATES(22)) ALGBRC(82) = ALGBRC(76)*(STATES(2) - ALGBRC(79)) ALGBRC(86) = CONSTS(149)*(CONSTS(152) ** CONSTS(150)/(CONSTS(152) ** CONSTS(150)+STATES(22) ** CONSTS(150))) ALGBRC(87) = ALGBRC(86)*(STATES(2) - CONSTS(151)) ALGBRC(90) = CONSTS(154)*(STATES(22) ** 2.00000/(CONSTS(155) ** 2.00000+STATES(22) ** 2.00000)) ALGBRC(91) = CONSTS(722)*(ALGBRC(82)+ALGBRC(87)+ALGBRC(90)) ALGBRC(92) = ( CONSTS(160)*(STATES(22) - CONSTS(161)*STATES(37)))/(1.00000+ CONSTS(162)*STATES(22)+ CONSTS(163)*STATES(37)+ CONSTS(164)*STATES(22)*STATES(37)) ALGBRC(99) = CONSTS(188)*(STATES(38)/(STATES(38)+CONSTS(189))) RATES(22) = (ALGBRC(48)+ CONSTS(717)*ALGBRC(99)) - ( (CONSTS(190)/( CONSTS(722)*2.00000))*ALGBRC(91)+ CONSTS(718)*ALGBRC(92)) ALGBRC(94) = CONSTS(167) ** CONSTS(168)/(CONSTS(167) ** CONSTS(168)+STATES(38) ** CONSTS(168)) ALGBRC(96) = CONSTS(165)*(STATES(12) ** 2.00000/(CONSTS(166) ** 2.00000+STATES(12) ** 2.00000))*ALGBRC(94) ALGBRC(102) = 1.00000/(1.00000+( CONSTS(195)*CONSTS(196))/CONSTS(195)+STATES(38) ** 2.00000) RATES(38) = ALGBRC(102)*(ALGBRC(96) - ALGBRC(99)) ALGBRC(98) = ( CONSTS(171)*CONSTS(176)+ CONSTS(173)*STATES(12))/(CONSTS(176)+STATES(12)) ALGBRC(101) = ( (CONSTS(170)+CONSTS(174))*CONSTS(177))/(CONSTS(177)+STATES(12)) ALGBRC(104) = CONSTS(169)*( CONSTS(187)*ALGBRC(98)*STATES(39))/( CONSTS(187)*ALGBRC(98)+ALGBRC(101)) ** 4.00000*(STATES(37) - STATES(12)) RATES(12) = (ALGBRC(37)+ CONSTS(720)*ALGBRC(104)) - ( CONSTS(721)*ALGBRC(48)+ CONSTS(719)*ALGBRC(96)) RATES(37) = ALGBRC(92) - ALGBRC(104) ALGBRC(105) = ( CONSTS(172)*CONSTS(177)+ CONSTS(175)*STATES(12))/(CONSTS(177)+STATES(12)) ALGBRC(107) = ( CONSTS(186)*STATES(40))/( (1.00000+CONSTS(182)/STATES(12) ** CONSTS(184))*(1.00000+STATES(12)/CONSTS(183) ** CONSTS(185))) RATES(39) = ALGBRC(107)*(1.00000 - STATES(39)) - (( CONSTS(187)*ALGBRC(98)*ALGBRC(105))/( CONSTS(187)*ALGBRC(98)+ALGBRC(101)))*STATES(39) ALGBRC(108) = (( CONSTS(205)*CONSTS(204))/CONSTS(206))*log(CONSTS(209)/STATES(41)) ALGBRC(109) = CONSTS(208)*(STATES(2) - ALGBRC(108)) ALGBRC(117) = CONSTS(217)*(STATES(41) ** CONSTS(219)/(CONSTS(218) ** CONSTS(219)+STATES(41) ** CONSTS(219)))*(CONSTS(220) - STATES(2)) ALGBRC(112) = CONSTS(214)*(CONSTS(213) ** CONSTS(211)/(CONSTS(213) ** CONSTS(211)+STATES(23) ** CONSTS(211))) ALGBRC(113) = ALGBRC(112)*(STATES(2) - CONSTS(212)) ALGBRC(119) = CONSTS(728)*(ALGBRC(113)+ALGBRC(109)+ALGBRC(117)) RATES(41) = - (CONSTS(251)/( 1.00000*CONSTS(728)))*ALGBRC(119) ALGBRC(49) = CONSTS(694)*(STATES(13) - STATES(23)) ALGBRC(100) = ( CONSTS(200)*EXP( CONSTS(201)*STATES(2)))/(1.00000+EXP( CONSTS(202)*(STATES(2) - CONSTS(203)))) ALGBRC(103) = (( CONSTS(205)*CONSTS(204))/( 2.00000*CONSTS(206)))*log(CONSTS(207)/STATES(23)) ALGBRC(106) = ALGBRC(100)*(STATES(2) - ALGBRC(103)) ALGBRC(110) = CONSTS(210)*(CONSTS(213) ** CONSTS(211)/(CONSTS(213) ** CONSTS(211)+STATES(23) ** CONSTS(211))) ALGBRC(111) = ALGBRC(110)*(STATES(2) - CONSTS(212)) ALGBRC(114) = CONSTS(215)*(STATES(23) ** 2.00000/(CONSTS(216) ** 2.00000+STATES(23) ** 2.00000)) ALGBRC(115) = CONSTS(728)*(ALGBRC(106)+ALGBRC(111)+ALGBRC(114)) ALGBRC(116) = ( CONSTS(221)*(STATES(23) - CONSTS(222)*STATES(42)))/(1.00000+ CONSTS(223)*STATES(23)+ CONSTS(224)*STATES(42)+ CONSTS(225)*STATES(23)*STATES(42)) ALGBRC(123) = CONSTS(249)*(STATES(43)/(STATES(43)+CONSTS(250))) RATES(23) = (ALGBRC(49)+ CONSTS(723)*ALGBRC(123)) - ( (CONSTS(251)/( CONSTS(728)*2.00000))*ALGBRC(115)+ CONSTS(724)*ALGBRC(116)) ALGBRC(118) = CONSTS(228) ** CONSTS(229)/(CONSTS(228) ** CONSTS(229)+STATES(43) ** CONSTS(229)) ALGBRC(120) = CONSTS(226)*(STATES(13) ** 2.00000/(CONSTS(227) ** 2.00000+STATES(13) ** 2.00000))*ALGBRC(118) ALGBRC(126) = 1.00000/(1.00000+( CONSTS(256)*CONSTS(257))/CONSTS(256)+STATES(43) ** 2.00000) RATES(43) = ALGBRC(126)*(ALGBRC(120) - ALGBRC(123)) ALGBRC(122) = ( CONSTS(232)*CONSTS(237)+ CONSTS(234)*STATES(13))/(CONSTS(237)+STATES(13)) ALGBRC(125) = ( (CONSTS(231)+CONSTS(235))*CONSTS(238))/(CONSTS(238)+STATES(13)) ALGBRC(128) = CONSTS(230)*( CONSTS(248)*ALGBRC(122)*STATES(44))/( CONSTS(248)*ALGBRC(122)+ALGBRC(125)) ** 4.00000*(STATES(42) - STATES(13)) RATES(13) = (ALGBRC(38)+ CONSTS(726)*ALGBRC(128)) - ( CONSTS(727)*ALGBRC(49)+ CONSTS(725)*ALGBRC(120)) RATES(42) = ALGBRC(116) - ALGBRC(128) ALGBRC(129) = ( CONSTS(233)*CONSTS(238)+ CONSTS(236)*STATES(13))/(CONSTS(238)+STATES(13)) ALGBRC(131) = ( CONSTS(247)*STATES(45))/( (1.00000+CONSTS(243)/STATES(13) ** CONSTS(245))*(1.00000+STATES(13)/CONSTS(244) ** CONSTS(246))) RATES(44) = ALGBRC(131)*(1.00000 - STATES(44)) - (( CONSTS(248)*ALGBRC(122)*ALGBRC(129))/( CONSTS(248)*ALGBRC(122)+ALGBRC(125)))*STATES(44) ALGBRC(132) = (( CONSTS(266)*CONSTS(265))/CONSTS(267))*log(CONSTS(270)/STATES(46)) ALGBRC(133) = CONSTS(269)*(STATES(2) - ALGBRC(132)) ALGBRC(141) = CONSTS(278)*(STATES(46) ** CONSTS(280)/(CONSTS(279) ** CONSTS(280)+STATES(46) ** CONSTS(280)))*(CONSTS(281) - STATES(2)) ALGBRC(136) = CONSTS(275)*(CONSTS(274) ** CONSTS(272)/(CONSTS(274) ** CONSTS(272)+STATES(24) ** CONSTS(272))) ALGBRC(137) = ALGBRC(136)*(STATES(2) - CONSTS(273)) ALGBRC(143) = CONSTS(734)*(ALGBRC(137)+ALGBRC(133)+ALGBRC(141)) RATES(46) = - (CONSTS(312)/( 1.00000*CONSTS(734)))*ALGBRC(143) ALGBRC(50) = CONSTS(695)*(STATES(14) - STATES(24)) ALGBRC(124) = ( CONSTS(261)*EXP( CONSTS(262)*STATES(2)))/(1.00000+EXP( CONSTS(263)*(STATES(2) - CONSTS(264)))) ALGBRC(127) = (( CONSTS(266)*CONSTS(265))/( 2.00000*CONSTS(267)))*log(CONSTS(268)/STATES(24)) ALGBRC(130) = ALGBRC(124)*(STATES(2) - ALGBRC(127)) ALGBRC(134) = CONSTS(271)*(CONSTS(274) ** CONSTS(272)/(CONSTS(274) ** CONSTS(272)+STATES(24) ** CONSTS(272))) ALGBRC(135) = ALGBRC(134)*(STATES(2) - CONSTS(273)) ALGBRC(138) = CONSTS(276)*(STATES(24) ** 2.00000/(CONSTS(277) ** 2.00000+STATES(24) ** 2.00000)) ALGBRC(139) = CONSTS(734)*(ALGBRC(130)+ALGBRC(135)+ALGBRC(138)) ALGBRC(140) = ( CONSTS(282)*(STATES(24) - CONSTS(283)*STATES(47)))/(1.00000+ CONSTS(284)*STATES(24)+ CONSTS(285)*STATES(47)+ CONSTS(286)*STATES(24)*STATES(47)) ALGBRC(147) = CONSTS(310)*(STATES(48)/(STATES(48)+CONSTS(311))) RATES(24) = (ALGBRC(50)+ CONSTS(729)*ALGBRC(147)) - ( (CONSTS(312)/( CONSTS(734)*2.00000))*ALGBRC(139)+ CONSTS(730)*ALGBRC(140)) ALGBRC(142) = CONSTS(289) ** CONSTS(290)/(CONSTS(289) ** CONSTS(290)+STATES(48) ** CONSTS(290)) ALGBRC(144) = CONSTS(287)*(STATES(14) ** 2.00000/(CONSTS(288) ** 2.00000+STATES(14) ** 2.00000))*ALGBRC(142) ALGBRC(150) = 1.00000/(1.00000+( CONSTS(317)*CONSTS(318))/CONSTS(317)+STATES(48) ** 2.00000) RATES(48) = ALGBRC(150)*(ALGBRC(144) - ALGBRC(147)) ALGBRC(146) = ( CONSTS(293)*CONSTS(298)+ CONSTS(295)*STATES(14))/(CONSTS(298)+STATES(14)) ALGBRC(149) = ( (CONSTS(292)+CONSTS(296))*CONSTS(299))/(CONSTS(299)+STATES(14)) ALGBRC(152) = CONSTS(291)*( CONSTS(309)*ALGBRC(146)*STATES(49))/( CONSTS(309)*ALGBRC(146)+ALGBRC(149)) ** 4.00000*(STATES(47) - STATES(14)) RATES(14) = (ALGBRC(39)+ CONSTS(732)*ALGBRC(152)) - ( CONSTS(733)*ALGBRC(50)+ CONSTS(731)*ALGBRC(144)) RATES(47) = ALGBRC(140) - ALGBRC(152) ALGBRC(153) = ( CONSTS(294)*CONSTS(299)+ CONSTS(297)*STATES(14))/(CONSTS(299)+STATES(14)) ALGBRC(155) = ( CONSTS(308)*STATES(50))/( (1.00000+CONSTS(304)/STATES(14) ** CONSTS(306))*(1.00000+STATES(14)/CONSTS(305) ** CONSTS(307))) RATES(49) = ALGBRC(155)*(1.00000 - STATES(49)) - (( CONSTS(309)*ALGBRC(146)*ALGBRC(153))/( CONSTS(309)*ALGBRC(146)+ALGBRC(149)))*STATES(49) ALGBRC(156) = (( CONSTS(327)*CONSTS(326))/CONSTS(328))*log(CONSTS(331)/STATES(51)) ALGBRC(157) = CONSTS(330)*(STATES(2) - ALGBRC(156)) ALGBRC(165) = CONSTS(339)*(STATES(51) ** CONSTS(341)/(CONSTS(340) ** CONSTS(341)+STATES(51) ** CONSTS(341)))*(CONSTS(342) - STATES(2)) ALGBRC(160) = CONSTS(336)*(CONSTS(335) ** CONSTS(333)/(CONSTS(335) ** CONSTS(333)+STATES(25) ** CONSTS(333))) ALGBRC(161) = ALGBRC(160)*(STATES(2) - CONSTS(334)) ALGBRC(167) = CONSTS(740)*(ALGBRC(161)+ALGBRC(157)+ALGBRC(165)) RATES(51) = - (CONSTS(373)/( 1.00000*CONSTS(740)))*ALGBRC(167) ALGBRC(51) = CONSTS(696)*(STATES(15) - STATES(25)) ALGBRC(148) = ( CONSTS(322)*EXP( CONSTS(323)*STATES(2)))/(1.00000+EXP( CONSTS(324)*(STATES(2) - CONSTS(325)))) ALGBRC(151) = (( CONSTS(327)*CONSTS(326))/( 2.00000*CONSTS(328)))*log(CONSTS(329)/STATES(25)) ALGBRC(154) = ALGBRC(148)*(STATES(2) - ALGBRC(151)) ALGBRC(158) = CONSTS(332)*(CONSTS(335) ** CONSTS(333)/(CONSTS(335) ** CONSTS(333)+STATES(25) ** CONSTS(333))) ALGBRC(159) = ALGBRC(158)*(STATES(2) - CONSTS(334)) ALGBRC(162) = CONSTS(337)*(STATES(25) ** 2.00000/(CONSTS(338) ** 2.00000+STATES(25) ** 2.00000)) ALGBRC(163) = CONSTS(740)*(ALGBRC(154)+ALGBRC(159)+ALGBRC(162)) ALGBRC(164) = ( CONSTS(343)*(STATES(25) - CONSTS(344)*STATES(52)))/(1.00000+ CONSTS(345)*STATES(25)+ CONSTS(346)*STATES(52)+ CONSTS(347)*STATES(25)*STATES(52)) ALGBRC(171) = CONSTS(371)*(STATES(53)/(STATES(53)+CONSTS(372))) RATES(25) = (ALGBRC(51)+ CONSTS(735)*ALGBRC(171)) - ( (CONSTS(373)/( CONSTS(740)*2.00000))*ALGBRC(163)+ CONSTS(736)*ALGBRC(164)) ALGBRC(166) = CONSTS(350) ** CONSTS(351)/(CONSTS(350) ** CONSTS(351)+STATES(53) ** CONSTS(351)) ALGBRC(168) = CONSTS(348)*(STATES(15) ** 2.00000/(CONSTS(349) ** 2.00000+STATES(15) ** 2.00000))*ALGBRC(166) ALGBRC(174) = 1.00000/(1.00000+( CONSTS(378)*CONSTS(379))/CONSTS(378)+STATES(53) ** 2.00000) RATES(53) = ALGBRC(174)*(ALGBRC(168) - ALGBRC(171)) ALGBRC(170) = ( CONSTS(354)*CONSTS(359)+ CONSTS(356)*STATES(15))/(CONSTS(359)+STATES(15)) ALGBRC(173) = ( (CONSTS(353)+CONSTS(357))*CONSTS(360))/(CONSTS(360)+STATES(15)) ALGBRC(176) = CONSTS(352)*( CONSTS(370)*ALGBRC(170)*STATES(54))/( CONSTS(370)*ALGBRC(170)+ALGBRC(173)) ** 4.00000*(STATES(52) - STATES(15)) RATES(15) = (ALGBRC(40)+ CONSTS(738)*ALGBRC(176)) - ( CONSTS(739)*ALGBRC(51)+ CONSTS(737)*ALGBRC(168)) RATES(52) = ALGBRC(164) - ALGBRC(176) ALGBRC(177) = ( CONSTS(355)*CONSTS(360)+ CONSTS(358)*STATES(15))/(CONSTS(360)+STATES(15)) ALGBRC(179) = ( CONSTS(369)*STATES(55))/( (1.00000+CONSTS(365)/STATES(15) ** CONSTS(367))*(1.00000+STATES(15)/CONSTS(366) ** CONSTS(368))) RATES(54) = ALGBRC(179)*(1.00000 - STATES(54)) - (( CONSTS(370)*ALGBRC(170)*ALGBRC(177))/( CONSTS(370)*ALGBRC(170)+ALGBRC(173)))*STATES(54) ALGBRC(180) = (( CONSTS(388)*CONSTS(387))/CONSTS(389))*log(CONSTS(392)/STATES(56)) ALGBRC(181) = CONSTS(391)*(STATES(2) - ALGBRC(180)) ALGBRC(189) = CONSTS(400)*(STATES(56) ** CONSTS(402)/(CONSTS(401) ** CONSTS(402)+STATES(56) ** CONSTS(402)))*(CONSTS(403) - STATES(2)) ALGBRC(184) = CONSTS(397)*(CONSTS(396) ** CONSTS(394)/(CONSTS(396) ** CONSTS(394)+STATES(26) ** CONSTS(394))) ALGBRC(185) = ALGBRC(184)*(STATES(2) - CONSTS(395)) ALGBRC(191) = CONSTS(746)*(ALGBRC(185)+ALGBRC(181)+ALGBRC(189)) RATES(56) = - (CONSTS(434)/( 1.00000*CONSTS(746)))*ALGBRC(191) ALGBRC(52) = CONSTS(697)*(STATES(16) - STATES(26)) ALGBRC(172) = ( CONSTS(383)*EXP( CONSTS(384)*STATES(2)))/(1.00000+EXP( CONSTS(385)*(STATES(2) - CONSTS(386)))) ALGBRC(175) = (( CONSTS(388)*CONSTS(387))/( 2.00000*CONSTS(389)))*log(CONSTS(390)/STATES(26)) ALGBRC(178) = ALGBRC(172)*(STATES(2) - ALGBRC(175)) ALGBRC(182) = CONSTS(393)*(CONSTS(396) ** CONSTS(394)/(CONSTS(396) ** CONSTS(394)+STATES(26) ** CONSTS(394))) ALGBRC(183) = ALGBRC(182)*(STATES(2) - CONSTS(395)) ALGBRC(186) = CONSTS(398)*(STATES(26) ** 2.00000/(CONSTS(399) ** 2.00000+STATES(26) ** 2.00000)) ALGBRC(187) = CONSTS(746)*(ALGBRC(178)+ALGBRC(183)+ALGBRC(186)) ALGBRC(188) = ( CONSTS(404)*(STATES(26) - CONSTS(405)*STATES(57)))/(1.00000+ CONSTS(406)*STATES(26)+ CONSTS(407)*STATES(57)+ CONSTS(408)*STATES(26)*STATES(57)) ALGBRC(195) = CONSTS(432)*(STATES(58)/(STATES(58)+CONSTS(433))) RATES(26) = (ALGBRC(52)+ CONSTS(741)*ALGBRC(195)) - ( (CONSTS(434)/( CONSTS(746)*2.00000))*ALGBRC(187)+ CONSTS(742)*ALGBRC(188)) ALGBRC(190) = CONSTS(411) ** CONSTS(412)/(CONSTS(411) ** CONSTS(412)+STATES(58) ** CONSTS(412)) ALGBRC(192) = CONSTS(409)*(STATES(16) ** 2.00000/(CONSTS(410) ** 2.00000+STATES(16) ** 2.00000))*ALGBRC(190) ALGBRC(198) = 1.00000/(1.00000+( CONSTS(439)*CONSTS(440))/CONSTS(439)+STATES(58) ** 2.00000) RATES(58) = ALGBRC(198)*(ALGBRC(192) - ALGBRC(195)) ALGBRC(194) = ( CONSTS(415)*CONSTS(420)+ CONSTS(417)*STATES(16))/(CONSTS(420)+STATES(16)) ALGBRC(197) = ( (CONSTS(414)+CONSTS(418))*CONSTS(421))/(CONSTS(421)+STATES(16)) ALGBRC(200) = CONSTS(413)*( CONSTS(431)*ALGBRC(194)*STATES(59))/( CONSTS(431)*ALGBRC(194)+ALGBRC(197)) ** 4.00000*(STATES(57) - STATES(16)) RATES(16) = (ALGBRC(41)+ CONSTS(744)*ALGBRC(200)) - ( CONSTS(745)*ALGBRC(52)+ CONSTS(743)*ALGBRC(192)) RATES(57) = ALGBRC(188) - ALGBRC(200) ALGBRC(201) = ( CONSTS(416)*CONSTS(421)+ CONSTS(419)*STATES(16))/(CONSTS(421)+STATES(16)) ALGBRC(203) = ( CONSTS(430)*STATES(60))/( (1.00000+CONSTS(426)/STATES(16) ** CONSTS(428))*(1.00000+STATES(16)/CONSTS(427) ** CONSTS(429))) RATES(59) = ALGBRC(203)*(1.00000 - STATES(59)) - (( CONSTS(431)*ALGBRC(194)*ALGBRC(201))/( CONSTS(431)*ALGBRC(194)+ALGBRC(197)))*STATES(59) ALGBRC(204) = (( CONSTS(449)*CONSTS(448))/CONSTS(450))*log(CONSTS(453)/STATES(61)) ALGBRC(205) = CONSTS(452)*(STATES(2) - ALGBRC(204)) ALGBRC(213) = CONSTS(461)*(STATES(61) ** CONSTS(463)/(CONSTS(462) ** CONSTS(463)+STATES(61) ** CONSTS(463)))*(CONSTS(464) - STATES(2)) ALGBRC(208) = CONSTS(458)*(CONSTS(457) ** CONSTS(455)/(CONSTS(457) ** CONSTS(455)+STATES(27) ** CONSTS(455))) ALGBRC(209) = ALGBRC(208)*(STATES(2) - CONSTS(456)) ALGBRC(215) = CONSTS(752)*(ALGBRC(209)+ALGBRC(205)+ALGBRC(213)) RATES(61) = - (CONSTS(495)/( 1.00000*CONSTS(752)))*ALGBRC(215) ALGBRC(53) = CONSTS(698)*(STATES(17) - STATES(27)) ALGBRC(196) = ( CONSTS(444)*EXP( CONSTS(445)*STATES(2)))/(1.00000+EXP( CONSTS(446)*(STATES(2) - CONSTS(447)))) ALGBRC(199) = (( CONSTS(449)*CONSTS(448))/( 2.00000*CONSTS(450)))*log(CONSTS(451)/STATES(27)) ALGBRC(202) = ALGBRC(196)*(STATES(2) - ALGBRC(199)) ALGBRC(206) = CONSTS(454)*(CONSTS(457) ** CONSTS(455)/(CONSTS(457) ** CONSTS(455)+STATES(27) ** CONSTS(455))) ALGBRC(207) = ALGBRC(206)*(STATES(2) - CONSTS(456)) ALGBRC(210) = CONSTS(459)*(STATES(27) ** 2.00000/(CONSTS(460) ** 2.00000+STATES(27) ** 2.00000)) ALGBRC(211) = CONSTS(752)*(ALGBRC(202)+ALGBRC(207)+ALGBRC(210)) ALGBRC(212) = ( CONSTS(465)*(STATES(27) - CONSTS(466)*STATES(62)))/(1.00000+ CONSTS(467)*STATES(27)+ CONSTS(468)*STATES(62)+ CONSTS(469)*STATES(27)*STATES(62)) ALGBRC(219) = CONSTS(493)*(STATES(63)/(STATES(63)+CONSTS(494))) RATES(27) = (ALGBRC(53)+ CONSTS(747)*ALGBRC(219)) - ( (CONSTS(495)/( CONSTS(752)*2.00000))*ALGBRC(211)+ CONSTS(748)*ALGBRC(212)) ALGBRC(214) = CONSTS(472) ** CONSTS(473)/(CONSTS(472) ** CONSTS(473)+STATES(63) ** CONSTS(473)) ALGBRC(216) = CONSTS(470)*(STATES(17) ** 2.00000/(CONSTS(471) ** 2.00000+STATES(17) ** 2.00000))*ALGBRC(214) ALGBRC(222) = 1.00000/(1.00000+( CONSTS(500)*CONSTS(501))/CONSTS(500)+STATES(63) ** 2.00000) RATES(63) = ALGBRC(222)*(ALGBRC(216) - ALGBRC(219)) ALGBRC(218) = ( CONSTS(476)*CONSTS(481)+ CONSTS(478)*STATES(17))/(CONSTS(481)+STATES(17)) ALGBRC(221) = ( (CONSTS(475)+CONSTS(479))*CONSTS(482))/(CONSTS(482)+STATES(17)) ALGBRC(224) = CONSTS(474)*( CONSTS(492)*ALGBRC(218)*STATES(64))/( CONSTS(492)*ALGBRC(218)+ALGBRC(221)) ** 4.00000*(STATES(62) - STATES(17)) RATES(17) = (ALGBRC(42)+ CONSTS(750)*ALGBRC(224)) - ( CONSTS(751)*ALGBRC(53)+ CONSTS(749)*ALGBRC(216)) RATES(62) = ALGBRC(212) - ALGBRC(224) ALGBRC(225) = ( CONSTS(477)*CONSTS(482)+ CONSTS(480)*STATES(17))/(CONSTS(482)+STATES(17)) ALGBRC(227) = ( CONSTS(491)*STATES(65))/( (1.00000+CONSTS(487)/STATES(17) ** CONSTS(489))*(1.00000+STATES(17)/CONSTS(488) ** CONSTS(490))) RATES(64) = ALGBRC(227)*(1.00000 - STATES(64)) - (( CONSTS(492)*ALGBRC(218)*ALGBRC(225))/( CONSTS(492)*ALGBRC(218)+ALGBRC(221)))*STATES(64) ALGBRC(228) = (( CONSTS(510)*CONSTS(509))/CONSTS(511))*log(CONSTS(514)/STATES(66)) ALGBRC(229) = CONSTS(513)*(STATES(2) - ALGBRC(228)) ALGBRC(237) = CONSTS(522)*(STATES(66) ** CONSTS(524)/(CONSTS(523) ** CONSTS(524)+STATES(66) ** CONSTS(524)))*(CONSTS(525) - STATES(2)) ALGBRC(232) = CONSTS(519)*(CONSTS(518) ** CONSTS(516)/(CONSTS(518) ** CONSTS(516)+STATES(28) ** CONSTS(516))) ALGBRC(233) = ALGBRC(232)*(STATES(2) - CONSTS(517)) ALGBRC(239) = CONSTS(758)*(ALGBRC(233)+ALGBRC(229)+ALGBRC(237)) RATES(66) = - (CONSTS(556)/( 1.00000*CONSTS(758)))*ALGBRC(239) ALGBRC(54) = CONSTS(699)*(STATES(18) - STATES(28)) ALGBRC(220) = ( CONSTS(505)*EXP( CONSTS(506)*STATES(2)))/(1.00000+EXP( CONSTS(507)*(STATES(2) - CONSTS(508)))) ALGBRC(223) = (( CONSTS(510)*CONSTS(509))/( 2.00000*CONSTS(511)))*log(CONSTS(512)/STATES(28)) ALGBRC(226) = ALGBRC(220)*(STATES(2) - ALGBRC(223)) ALGBRC(230) = CONSTS(515)*(CONSTS(518) ** CONSTS(516)/(CONSTS(518) ** CONSTS(516)+STATES(28) ** CONSTS(516))) ALGBRC(231) = ALGBRC(230)*(STATES(2) - CONSTS(517)) ALGBRC(234) = CONSTS(520)*(STATES(28) ** 2.00000/(CONSTS(521) ** 2.00000+STATES(28) ** 2.00000)) ALGBRC(235) = CONSTS(758)*(ALGBRC(226)+ALGBRC(231)+ALGBRC(234)) ALGBRC(236) = ( CONSTS(526)*(STATES(28) - CONSTS(527)*STATES(67)))/(1.00000+ CONSTS(528)*STATES(28)+ CONSTS(529)*STATES(67)+ CONSTS(530)*STATES(28)*STATES(67)) ALGBRC(243) = CONSTS(554)*(STATES(68)/(STATES(68)+CONSTS(555))) RATES(28) = (ALGBRC(54)+ CONSTS(753)*ALGBRC(243)) - ( (CONSTS(556)/( CONSTS(758)*2.00000))*ALGBRC(235)+ CONSTS(754)*ALGBRC(236)) ALGBRC(238) = CONSTS(533) ** CONSTS(534)/(CONSTS(533) ** CONSTS(534)+STATES(68) ** CONSTS(534)) ALGBRC(240) = CONSTS(531)*(STATES(18) ** 2.00000/(CONSTS(532) ** 2.00000+STATES(18) ** 2.00000))*ALGBRC(238) ALGBRC(246) = 1.00000/(1.00000+( CONSTS(561)*CONSTS(562))/CONSTS(561)+STATES(68) ** 2.00000) RATES(68) = ALGBRC(246)*(ALGBRC(240) - ALGBRC(243)) ALGBRC(242) = ( CONSTS(537)*CONSTS(542)+ CONSTS(539)*STATES(18))/(CONSTS(542)+STATES(18)) ALGBRC(245) = ( (CONSTS(536)+CONSTS(540))*CONSTS(543))/(CONSTS(543)+STATES(18)) ALGBRC(248) = CONSTS(535)*( CONSTS(553)*ALGBRC(242)*STATES(69))/( CONSTS(553)*ALGBRC(242)+ALGBRC(245)) ** 4.00000*(STATES(67) - STATES(18)) RATES(18) = (ALGBRC(43)+ CONSTS(756)*ALGBRC(248)) - ( CONSTS(757)*ALGBRC(54)+ CONSTS(755)*ALGBRC(240)) RATES(67) = ALGBRC(236) - ALGBRC(248) ALGBRC(249) = ( CONSTS(538)*CONSTS(543)+ CONSTS(541)*STATES(18))/(CONSTS(543)+STATES(18)) ALGBRC(251) = ( CONSTS(552)*STATES(70))/( (1.00000+CONSTS(548)/STATES(18) ** CONSTS(550))*(1.00000+STATES(18)/CONSTS(549) ** CONSTS(551))) RATES(69) = ALGBRC(251)*(1.00000 - STATES(69)) - (( CONSTS(553)*ALGBRC(242)*ALGBRC(249))/( CONSTS(553)*ALGBRC(242)+ALGBRC(245)))*STATES(69) ALGBRC(252) = (( CONSTS(571)*CONSTS(570))/CONSTS(572))*log(CONSTS(575)/STATES(71)) ALGBRC(253) = CONSTS(574)*(STATES(2) - ALGBRC(252)) ALGBRC(261) = CONSTS(583)*(STATES(71) ** CONSTS(585)/(CONSTS(584) ** CONSTS(585)+STATES(71) ** CONSTS(585)))*(CONSTS(586) - STATES(2)) ALGBRC(256) = CONSTS(580)*(CONSTS(579) ** CONSTS(577)/(CONSTS(579) ** CONSTS(577)+STATES(29) ** CONSTS(577))) ALGBRC(257) = ALGBRC(256)*(STATES(2) - CONSTS(578)) ALGBRC(263) = CONSTS(764)*(ALGBRC(257)+ALGBRC(253)+ALGBRC(261)) RATES(71) = - (CONSTS(617)/( 1.00000*CONSTS(764)))*ALGBRC(263) ALGBRC(55) = CONSTS(700)*(STATES(19) - STATES(29)) ALGBRC(244) = ( CONSTS(566)*EXP( CONSTS(567)*STATES(2)))/(1.00000+EXP( CONSTS(568)*(STATES(2) - CONSTS(569)))) ALGBRC(247) = (( CONSTS(571)*CONSTS(570))/( 2.00000*CONSTS(572)))*log(CONSTS(573)/STATES(29)) ALGBRC(250) = ALGBRC(244)*(STATES(2) - ALGBRC(247)) ALGBRC(254) = CONSTS(576)*(CONSTS(579) ** CONSTS(577)/(CONSTS(579) ** CONSTS(577)+STATES(29) ** CONSTS(577))) ALGBRC(255) = ALGBRC(254)*(STATES(2) - CONSTS(578)) ALGBRC(258) = CONSTS(581)*(STATES(29) ** 2.00000/(CONSTS(582) ** 2.00000+STATES(29) ** 2.00000)) ALGBRC(259) = CONSTS(764)*(ALGBRC(250)+ALGBRC(255)+ALGBRC(258)) ALGBRC(260) = ( CONSTS(587)*(STATES(29) - CONSTS(588)*STATES(72)))/(1.00000+ CONSTS(589)*STATES(29)+ CONSTS(590)*STATES(72)+ CONSTS(591)*STATES(29)*STATES(72)) ALGBRC(267) = CONSTS(615)*(STATES(73)/(STATES(73)+CONSTS(616))) RATES(29) = (ALGBRC(55)+ CONSTS(759)*ALGBRC(267)) - ( (CONSTS(617)/( CONSTS(764)*2.00000))*ALGBRC(259)+ CONSTS(760)*ALGBRC(260)) ALGBRC(262) = CONSTS(594) ** CONSTS(595)/(CONSTS(594) ** CONSTS(595)+STATES(73) ** CONSTS(595)) ALGBRC(264) = CONSTS(592)*(STATES(19) ** 2.00000/(CONSTS(593) ** 2.00000+STATES(19) ** 2.00000))*ALGBRC(262) ALGBRC(270) = 1.00000/(1.00000+( CONSTS(622)*CONSTS(623))/CONSTS(622)+STATES(73) ** 2.00000) RATES(73) = ALGBRC(270)*(ALGBRC(264) - ALGBRC(267)) ALGBRC(266) = ( CONSTS(598)*CONSTS(603)+ CONSTS(600)*STATES(19))/(CONSTS(603)+STATES(19)) ALGBRC(269) = ( (CONSTS(597)+CONSTS(601))*CONSTS(604))/(CONSTS(604)+STATES(19)) ALGBRC(272) = CONSTS(596)*( CONSTS(614)*ALGBRC(266)*STATES(74))/( CONSTS(614)*ALGBRC(266)+ALGBRC(269)) ** 4.00000*(STATES(72) - STATES(19)) RATES(19) = (ALGBRC(44)+ CONSTS(762)*ALGBRC(272)) - ( CONSTS(763)*ALGBRC(55)+ CONSTS(761)*ALGBRC(264)) RATES(72) = ALGBRC(260) - ALGBRC(272) ALGBRC(273) = ( CONSTS(599)*CONSTS(604)+ CONSTS(602)*STATES(19))/(CONSTS(604)+STATES(19)) ALGBRC(275) = ( CONSTS(613)*STATES(75))/( (1.00000+CONSTS(609)/STATES(19) ** CONSTS(611))*(1.00000+STATES(19)/CONSTS(610) ** CONSTS(612))) RATES(74) = ALGBRC(275)*(1.00000 - STATES(74)) - (( CONSTS(614)*ALGBRC(266)*ALGBRC(273))/( CONSTS(614)*ALGBRC(266)+ALGBRC(269)))*STATES(74) ALGBRC(276) = (( CONSTS(632)*CONSTS(631))/CONSTS(633))*log(CONSTS(636)/STATES(76)) ALGBRC(277) = CONSTS(635)*(STATES(2) - ALGBRC(276)) ALGBRC(284) = CONSTS(644)*(STATES(76) ** CONSTS(646)/(CONSTS(645) ** CONSTS(646)+STATES(76) ** CONSTS(646)))*(CONSTS(647) - STATES(2)) ALGBRC(280) = CONSTS(641)*(CONSTS(640) ** CONSTS(638)/(CONSTS(640) ** CONSTS(638)+STATES(30) ** CONSTS(638))) ALGBRC(281) = ALGBRC(280)*(STATES(2) - CONSTS(639)) ALGBRC(286) = CONSTS(770)*(ALGBRC(281)+ALGBRC(277)+ALGBRC(284)) RATES(76) = - (CONSTS(678)/( 1.00000*CONSTS(770)))*ALGBRC(286) ALGBRC(30) = CONSTS(16)*STATES(6)*STATES(7)*(STATES(2) - CONSTS(689)) ALGBRC(31) = CONSTS(28)*STATES(8)*(STATES(2) - CONSTS(689)) ALGBRC(32) = CONSTS(33)*(STATES(2) - CONSTS(689)) ALGBRC(33) = CONSTS(34)*STATES(9)*STATES(10)*(STATES(2) - CONSTS(690)) ALGBRC(34) = CONSTS(49)*(STATES(2) - CONSTS(50)) ALGBRC(35) = ALGBRC(8)+ALGBRC(26)+ALGBRC(30)+ALGBRC(31)+ALGBRC(32)+ALGBRC(33)+ALGBRC(34) ALGBRC(73) = ALGBRC(67)+ALGBRC(71) ALGBRC(97) = ALGBRC(91)+ALGBRC(95) ALGBRC(121) = ALGBRC(115)+ALGBRC(119) ALGBRC(145) = ALGBRC(139)+ALGBRC(143) ALGBRC(169) = ALGBRC(163)+ALGBRC(167) ALGBRC(193) = ALGBRC(187)+ALGBRC(191) ALGBRC(217) = ALGBRC(211)+ALGBRC(215) ALGBRC(241) = ALGBRC(235)+ALGBRC(239) ALGBRC(265) = ALGBRC(259)+ALGBRC(263) ALGBRC(268) = ( CONSTS(627)*EXP( CONSTS(628)*STATES(2)))/(1.00000+EXP( CONSTS(629)*(STATES(2) - CONSTS(630)))) ALGBRC(271) = (( CONSTS(632)*CONSTS(631))/( 2.00000*CONSTS(633)))*log(CONSTS(634)/STATES(30)) ALGBRC(274) = ALGBRC(268)*(STATES(2) - ALGBRC(271)) ALGBRC(278) = CONSTS(637)*(CONSTS(640) ** CONSTS(638)/(CONSTS(640) ** CONSTS(638)+STATES(30) ** CONSTS(638))) ALGBRC(279) = ALGBRC(278)*(STATES(2) - CONSTS(639)) ALGBRC(282) = CONSTS(642)*(STATES(30) ** 2.00000/(CONSTS(643) ** 2.00000+STATES(30) ** 2.00000)) ALGBRC(283) = CONSTS(770)*(ALGBRC(274)+ALGBRC(279)+ALGBRC(282)) ALGBRC(288) = ALGBRC(283)+ALGBRC(286) RATES(2) = - (1.00000/CONSTS(51))*(ALGBRC(35)+ALGBRC(73)+ALGBRC(97)+ALGBRC(121)+ALGBRC(145)+ALGBRC(169)+ALGBRC(193)+ALGBRC(217)+ALGBRC(241)+ALGBRC(265)+ALGBRC(288)) ALGBRC(56) = CONSTS(701)*(STATES(20) - STATES(30)) ALGBRC(285) = ( CONSTS(648)*(STATES(30) - CONSTS(649)*STATES(77)))/(1.00000+ CONSTS(650)*STATES(30)+ CONSTS(651)*STATES(77)+ CONSTS(652)*STATES(30)*STATES(77)) ALGBRC(291) = CONSTS(676)*(STATES(78)/(STATES(78)+CONSTS(677))) RATES(30) = (ALGBRC(56)+ CONSTS(765)*ALGBRC(291)) - ( (CONSTS(678)/( CONSTS(770)*2.00000))*ALGBRC(283)+ CONSTS(766)*ALGBRC(285)) ALGBRC(287) = CONSTS(655) ** CONSTS(656)/(CONSTS(655) ** CONSTS(656)+STATES(78) ** CONSTS(656)) ALGBRC(289) = CONSTS(653)*(STATES(20) ** 2.00000/(CONSTS(654) ** 2.00000+STATES(20) ** 2.00000))*ALGBRC(287) ALGBRC(293) = 1.00000/(1.00000+( CONSTS(683)*CONSTS(684))/CONSTS(683)+STATES(78) ** 2.00000) RATES(78) = ALGBRC(293)*(ALGBRC(289) - ALGBRC(291)) ALGBRC(290) = ( CONSTS(659)*CONSTS(664)+ CONSTS(661)*STATES(20))/(CONSTS(664)+STATES(20)) ALGBRC(292) = ( (CONSTS(658)+CONSTS(662))*CONSTS(665))/(CONSTS(665)+STATES(20)) ALGBRC(294) = CONSTS(657)*( CONSTS(675)*ALGBRC(290)*STATES(79))/( CONSTS(675)*ALGBRC(290)+ALGBRC(292)) ** 4.00000*(STATES(77) - STATES(20)) RATES(20) = (ALGBRC(45)+ CONSTS(768)*ALGBRC(294)) - ( CONSTS(769)*ALGBRC(56)+ CONSTS(767)*ALGBRC(289)) RATES(77) = ALGBRC(285) - ALGBRC(294) ALGBRC(295) = ( CONSTS(660)*CONSTS(665)+ CONSTS(663)*STATES(20))/(CONSTS(665)+STATES(20)) ALGBRC(296) = ( CONSTS(674)*STATES(80))/( (1.00000+CONSTS(670)/STATES(20) ** CONSTS(672))*(1.00000+STATES(20)/CONSTS(671) ** CONSTS(673))) RATES(79) = ALGBRC(296)*(1.00000 - STATES(79)) - (( CONSTS(675)*ALGBRC(290)*ALGBRC(295))/( CONSTS(675)*ALGBRC(290)+ALGBRC(292)))*STATES(79) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = 1.00000/(1.00000+EXP( CONSTS(6)*(STATES(2) - CONSTS(7)))) ALGBRC(10) = CONSTS(117)/(1.00000+CONSTS(120)/STATES(11) ** CONSTS(118)) ALGBRC(11) = CONSTS(178)/(1.00000+CONSTS(181)/STATES(12) ** CONSTS(179)) ALGBRC(12) = CONSTS(239)/(1.00000+CONSTS(242)/STATES(13) ** CONSTS(240)) ALGBRC(13) = CONSTS(300)/(1.00000+CONSTS(303)/STATES(14) ** CONSTS(301)) ALGBRC(14) = CONSTS(361)/(1.00000+CONSTS(364)/STATES(15) ** CONSTS(362)) ALGBRC(15) = CONSTS(422)/(1.00000+CONSTS(425)/STATES(16) ** CONSTS(423)) ALGBRC(16) = CONSTS(483)/(1.00000+CONSTS(486)/STATES(17) ** CONSTS(484)) ALGBRC(17) = CONSTS(544)/(1.00000+CONSTS(547)/STATES(18) ** CONSTS(545)) ALGBRC(18) = CONSTS(605)/(1.00000+CONSTS(608)/STATES(19) ** CONSTS(606)) ALGBRC(19) = CONSTS(666)/(1.00000+CONSTS(669)/STATES(20) ** CONSTS(667)) ALGBRC(2) = 1.00000/(1.00000+EXP( CONSTS(12)*(STATES(2) - CONSTS(13)))) ALGBRC(20) = CONSTS(8)+(CONSTS(9) - CONSTS(8))/(1.00000+EXP( CONSTS(10)*(STATES(2) - CONSTS(11)))) ALGBRC(4) = CONSTS(24)*(CONSTS(25)/(1.00000+EXP( CONSTS(26)*(STATES(2) - CONSTS(27))))+(1.00000 - CONSTS(25))) ALGBRC(22) = CONSTS(24)*(CONSTS(25) - CONSTS(25)/(1.00000+EXP( CONSTS(26)*(STATES(2) - CONSTS(27))))) ALGBRC(5) = CONSTS(29)*(CONSTS(30)/(1.00000+EXP( CONSTS(31)*(STATES(2) - CONSTS(32))))+(1.00000 - CONSTS(30))) ALGBRC(23) = CONSTS(29)*(CONSTS(30) - CONSTS(30)/(1.00000+EXP( CONSTS(31)*(STATES(2) - CONSTS(32))))) ALGBRC(3) = CONSTS(17)+(CONSTS(18) - CONSTS(17))/(1.00000+EXP( CONSTS(19)*(STATES(2) - CONSTS(20)))) ALGBRC(21) = ALGBRC(3)*(CONSTS(21)/(1.00000+EXP( CONSTS(22)*(STATES(2) - CONSTS(23))))+(1.00000 - CONSTS(21))) ALGBRC(27) = ALGBRC(3)*(CONSTS(21) - CONSTS(21)/(1.00000+EXP( CONSTS(22)*(STATES(2) - CONSTS(23))))) ALGBRC(6) = CONSTS(35)+(CONSTS(36) - CONSTS(35))/(1.00000+EXP( CONSTS(37)*(STATES(2) - CONSTS(38)))) ALGBRC(24) = ALGBRC(6)*(CONSTS(39)/(1.00000+EXP( CONSTS(40)*(STATES(2) - CONSTS(41))))+(1.00000 - CONSTS(39))) ALGBRC(28) = ALGBRC(6)*(CONSTS(39) - CONSTS(39)/(1.00000+EXP( CONSTS(40)*(STATES(2) - CONSTS(41))))) ALGBRC(7) = CONSTS(42)+(CONSTS(43) - CONSTS(42))/(1.00000+EXP( CONSTS(44)*(STATES(2) - CONSTS(45)))) ALGBRC(25) = ALGBRC(7)*(CONSTS(46)/(1.00000+EXP( CONSTS(47)*(STATES(2) - CONSTS(48))))+(1.00000 - CONSTS(46))) ALGBRC(29) = ALGBRC(7)*(CONSTS(46) - CONSTS(46)/(1.00000+EXP( CONSTS(47)*(STATES(2) - CONSTS(48))))) ALGBRC(8) = CONSTS(1)*STATES(1)*(STATES(2) - CONSTS(2)) ALGBRC(26) = CONSTS(14)*(STATES(5)/(CONSTS(15)+STATES(5))) ALGBRC(46) = CONSTS(75)*(CONSTS(76) - STATES(5)) ALGBRC(36) = CONSTS(62)*(STATES(5) - STATES(11)) ALGBRC(37) = CONSTS(63)*(STATES(5) - STATES(12)) ALGBRC(38) = CONSTS(64)*(STATES(5) - STATES(13)) ALGBRC(39) = CONSTS(65)*(STATES(5) - STATES(14)) ALGBRC(40) = CONSTS(66)*(STATES(5) - STATES(15)) ALGBRC(41) = CONSTS(67)*(STATES(5) - STATES(16)) ALGBRC(42) = CONSTS(68)*(STATES(5) - STATES(17)) ALGBRC(43) = CONSTS(69)*(STATES(5) - STATES(18)) ALGBRC(44) = CONSTS(70)*(STATES(5) - STATES(19)) ALGBRC(45) = CONSTS(71)*(STATES(5) - STATES(20)) ALGBRC(60) = (( CONSTS(83)*CONSTS(82))/CONSTS(84))*log(CONSTS(87)/STATES(31)) ALGBRC(61) = CONSTS(86)*(STATES(2) - ALGBRC(60)) ALGBRC(69) = CONSTS(95)*(STATES(31) ** CONSTS(97)/(CONSTS(96) ** CONSTS(97)+STATES(31) ** CONSTS(97)))*(CONSTS(98) - STATES(2)) ALGBRC(64) = CONSTS(92)*(CONSTS(91) ** CONSTS(89)/(CONSTS(91) ** CONSTS(89)+STATES(21) ** CONSTS(89))) ALGBRC(65) = ALGBRC(64)*(STATES(2) - CONSTS(90)) ALGBRC(71) = CONSTS(716)*(ALGBRC(65)+ALGBRC(61)+ALGBRC(69)) ALGBRC(47) = CONSTS(692)*(STATES(11) - STATES(21)) ALGBRC(57) = ( CONSTS(78)*EXP( CONSTS(79)*STATES(2)))/(1.00000+EXP( CONSTS(80)*(STATES(2) - CONSTS(81)))) ALGBRC(58) = (( CONSTS(83)*CONSTS(82))/( 2.00000*CONSTS(84)))*log(CONSTS(85)/STATES(21)) ALGBRC(59) = ALGBRC(57)*(STATES(2) - ALGBRC(58)) ALGBRC(62) = CONSTS(88)*(CONSTS(91) ** CONSTS(89)/(CONSTS(91) ** CONSTS(89)+STATES(21) ** CONSTS(89))) ALGBRC(63) = ALGBRC(62)*(STATES(2) - CONSTS(90)) ALGBRC(66) = CONSTS(93)*(STATES(21) ** 2.00000/(CONSTS(94) ** 2.00000+STATES(21) ** 2.00000)) ALGBRC(67) = CONSTS(716)*(ALGBRC(59)+ALGBRC(63)+ALGBRC(66)) ALGBRC(68) = ( CONSTS(99)*(STATES(21) - CONSTS(100)*STATES(32)))/(1.00000+ CONSTS(101)*STATES(21)+ CONSTS(102)*STATES(32)+ CONSTS(103)*STATES(21)*STATES(32)) ALGBRC(75) = CONSTS(127)*(STATES(33)/(STATES(33)+CONSTS(128))) ALGBRC(70) = CONSTS(106) ** CONSTS(107)/(CONSTS(106) ** CONSTS(107)+STATES(33) ** CONSTS(107)) ALGBRC(72) = CONSTS(104)*(STATES(11) ** 2.00000/(CONSTS(105) ** 2.00000+STATES(11) ** 2.00000))*ALGBRC(70) ALGBRC(78) = 1.00000/(1.00000+( CONSTS(134)*CONSTS(135))/CONSTS(134)+STATES(33) ** 2.00000) ALGBRC(74) = ( CONSTS(110)*CONSTS(115)+ CONSTS(112)*STATES(11))/(CONSTS(115)+STATES(11)) ALGBRC(77) = ( (CONSTS(109)+CONSTS(113))*CONSTS(116))/(CONSTS(116)+STATES(11)) ALGBRC(80) = CONSTS(108)*( CONSTS(126)*ALGBRC(74)*STATES(34))/( CONSTS(126)*ALGBRC(74)+ALGBRC(77)) ** 4.00000*(STATES(32) - STATES(11)) ALGBRC(81) = ( CONSTS(111)*CONSTS(116)+ CONSTS(114)*STATES(11))/(CONSTS(116)+STATES(11)) ALGBRC(83) = ( CONSTS(125)*STATES(35))/( (1.00000+CONSTS(121)/STATES(11) ** CONSTS(123))*(1.00000+STATES(11)/CONSTS(122) ** CONSTS(124))) ALGBRC(84) = (( CONSTS(144)*CONSTS(143))/CONSTS(145))*log(CONSTS(148)/STATES(36)) ALGBRC(85) = CONSTS(147)*(STATES(2) - ALGBRC(84)) ALGBRC(93) = CONSTS(156)*(STATES(36) ** CONSTS(158)/(CONSTS(157) ** CONSTS(158)+STATES(36) ** CONSTS(158)))*(CONSTS(159) - STATES(2)) ALGBRC(88) = CONSTS(153)*(CONSTS(152) ** CONSTS(150)/(CONSTS(152) ** CONSTS(150)+STATES(22) ** CONSTS(150))) ALGBRC(89) = ALGBRC(88)*(STATES(2) - CONSTS(151)) ALGBRC(95) = CONSTS(722)*(ALGBRC(89)+ALGBRC(85)+ALGBRC(93)) ALGBRC(48) = CONSTS(693)*(STATES(12) - STATES(22)) ALGBRC(76) = ( CONSTS(139)*EXP( CONSTS(140)*STATES(2)))/(1.00000+EXP( CONSTS(141)*(STATES(2) - CONSTS(142)))) ALGBRC(79) = (( CONSTS(144)*CONSTS(143))/( 2.00000*CONSTS(145)))*log(CONSTS(146)/STATES(22)) ALGBRC(82) = ALGBRC(76)*(STATES(2) - ALGBRC(79)) ALGBRC(86) = CONSTS(149)*(CONSTS(152) ** CONSTS(150)/(CONSTS(152) ** CONSTS(150)+STATES(22) ** CONSTS(150))) ALGBRC(87) = ALGBRC(86)*(STATES(2) - CONSTS(151)) ALGBRC(90) = CONSTS(154)*(STATES(22) ** 2.00000/(CONSTS(155) ** 2.00000+STATES(22) ** 2.00000)) ALGBRC(91) = CONSTS(722)*(ALGBRC(82)+ALGBRC(87)+ALGBRC(90)) ALGBRC(92) = ( CONSTS(160)*(STATES(22) - CONSTS(161)*STATES(37)))/(1.00000+ CONSTS(162)*STATES(22)+ CONSTS(163)*STATES(37)+ CONSTS(164)*STATES(22)*STATES(37)) ALGBRC(99) = CONSTS(188)*(STATES(38)/(STATES(38)+CONSTS(189))) ALGBRC(94) = CONSTS(167) ** CONSTS(168)/(CONSTS(167) ** CONSTS(168)+STATES(38) ** CONSTS(168)) ALGBRC(96) = CONSTS(165)*(STATES(12) ** 2.00000/(CONSTS(166) ** 2.00000+STATES(12) ** 2.00000))*ALGBRC(94) ALGBRC(102) = 1.00000/(1.00000+( CONSTS(195)*CONSTS(196))/CONSTS(195)+STATES(38) ** 2.00000) ALGBRC(98) = ( CONSTS(171)*CONSTS(176)+ CONSTS(173)*STATES(12))/(CONSTS(176)+STATES(12)) ALGBRC(101) = ( (CONSTS(170)+CONSTS(174))*CONSTS(177))/(CONSTS(177)+STATES(12)) ALGBRC(104) = CONSTS(169)*( CONSTS(187)*ALGBRC(98)*STATES(39))/( CONSTS(187)*ALGBRC(98)+ALGBRC(101)) ** 4.00000*(STATES(37) - STATES(12)) ALGBRC(105) = ( CONSTS(172)*CONSTS(177)+ CONSTS(175)*STATES(12))/(CONSTS(177)+STATES(12)) ALGBRC(107) = ( CONSTS(186)*STATES(40))/( (1.00000+CONSTS(182)/STATES(12) ** CONSTS(184))*(1.00000+STATES(12)/CONSTS(183) ** CONSTS(185))) ALGBRC(108) = (( CONSTS(205)*CONSTS(204))/CONSTS(206))*log(CONSTS(209)/STATES(41)) ALGBRC(109) = CONSTS(208)*(STATES(2) - ALGBRC(108)) ALGBRC(117) = CONSTS(217)*(STATES(41) ** CONSTS(219)/(CONSTS(218) ** CONSTS(219)+STATES(41) ** CONSTS(219)))*(CONSTS(220) - STATES(2)) ALGBRC(112) = CONSTS(214)*(CONSTS(213) ** CONSTS(211)/(CONSTS(213) ** CONSTS(211)+STATES(23) ** CONSTS(211))) ALGBRC(113) = ALGBRC(112)*(STATES(2) - CONSTS(212)) ALGBRC(119) = CONSTS(728)*(ALGBRC(113)+ALGBRC(109)+ALGBRC(117)) ALGBRC(49) = CONSTS(694)*(STATES(13) - STATES(23)) ALGBRC(100) = ( CONSTS(200)*EXP( CONSTS(201)*STATES(2)))/(1.00000+EXP( CONSTS(202)*(STATES(2) - CONSTS(203)))) ALGBRC(103) = (( CONSTS(205)*CONSTS(204))/( 2.00000*CONSTS(206)))*log(CONSTS(207)/STATES(23)) ALGBRC(106) = ALGBRC(100)*(STATES(2) - ALGBRC(103)) ALGBRC(110) = CONSTS(210)*(CONSTS(213) ** CONSTS(211)/(CONSTS(213) ** CONSTS(211)+STATES(23) ** CONSTS(211))) ALGBRC(111) = ALGBRC(110)*(STATES(2) - CONSTS(212)) ALGBRC(114) = CONSTS(215)*(STATES(23) ** 2.00000/(CONSTS(216) ** 2.00000+STATES(23) ** 2.00000)) ALGBRC(115) = CONSTS(728)*(ALGBRC(106)+ALGBRC(111)+ALGBRC(114)) ALGBRC(116) = ( CONSTS(221)*(STATES(23) - CONSTS(222)*STATES(42)))/(1.00000+ CONSTS(223)*STATES(23)+ CONSTS(224)*STATES(42)+ CONSTS(225)*STATES(23)*STATES(42)) ALGBRC(123) = CONSTS(249)*(STATES(43)/(STATES(43)+CONSTS(250))) ALGBRC(118) = CONSTS(228) ** CONSTS(229)/(CONSTS(228) ** CONSTS(229)+STATES(43) ** CONSTS(229)) ALGBRC(120) = CONSTS(226)*(STATES(13) ** 2.00000/(CONSTS(227) ** 2.00000+STATES(13) ** 2.00000))*ALGBRC(118) ALGBRC(126) = 1.00000/(1.00000+( CONSTS(256)*CONSTS(257))/CONSTS(256)+STATES(43) ** 2.00000) ALGBRC(122) = ( CONSTS(232)*CONSTS(237)+ CONSTS(234)*STATES(13))/(CONSTS(237)+STATES(13)) ALGBRC(125) = ( (CONSTS(231)+CONSTS(235))*CONSTS(238))/(CONSTS(238)+STATES(13)) ALGBRC(128) = CONSTS(230)*( CONSTS(248)*ALGBRC(122)*STATES(44))/( CONSTS(248)*ALGBRC(122)+ALGBRC(125)) ** 4.00000*(STATES(42) - STATES(13)) ALGBRC(129) = ( CONSTS(233)*CONSTS(238)+ CONSTS(236)*STATES(13))/(CONSTS(238)+STATES(13)) ALGBRC(131) = ( CONSTS(247)*STATES(45))/( (1.00000+CONSTS(243)/STATES(13) ** CONSTS(245))*(1.00000+STATES(13)/CONSTS(244) ** CONSTS(246))) ALGBRC(132) = (( CONSTS(266)*CONSTS(265))/CONSTS(267))*log(CONSTS(270)/STATES(46)) ALGBRC(133) = CONSTS(269)*(STATES(2) - ALGBRC(132)) ALGBRC(141) = CONSTS(278)*(STATES(46) ** CONSTS(280)/(CONSTS(279) ** CONSTS(280)+STATES(46) ** CONSTS(280)))*(CONSTS(281) - STATES(2)) ALGBRC(136) = CONSTS(275)*(CONSTS(274) ** CONSTS(272)/(CONSTS(274) ** CONSTS(272)+STATES(24) ** CONSTS(272))) ALGBRC(137) = ALGBRC(136)*(STATES(2) - CONSTS(273)) ALGBRC(143) = CONSTS(734)*(ALGBRC(137)+ALGBRC(133)+ALGBRC(141)) ALGBRC(50) = CONSTS(695)*(STATES(14) - STATES(24)) ALGBRC(124) = ( CONSTS(261)*EXP( CONSTS(262)*STATES(2)))/(1.00000+EXP( CONSTS(263)*(STATES(2) - CONSTS(264)))) ALGBRC(127) = (( CONSTS(266)*CONSTS(265))/( 2.00000*CONSTS(267)))*log(CONSTS(268)/STATES(24)) ALGBRC(130) = ALGBRC(124)*(STATES(2) - ALGBRC(127)) ALGBRC(134) = CONSTS(271)*(CONSTS(274) ** CONSTS(272)/(CONSTS(274) ** CONSTS(272)+STATES(24) ** CONSTS(272))) ALGBRC(135) = ALGBRC(134)*(STATES(2) - CONSTS(273)) ALGBRC(138) = CONSTS(276)*(STATES(24) ** 2.00000/(CONSTS(277) ** 2.00000+STATES(24) ** 2.00000)) ALGBRC(139) = CONSTS(734)*(ALGBRC(130)+ALGBRC(135)+ALGBRC(138)) ALGBRC(140) = ( CONSTS(282)*(STATES(24) - CONSTS(283)*STATES(47)))/(1.00000+ CONSTS(284)*STATES(24)+ CONSTS(285)*STATES(47)+ CONSTS(286)*STATES(24)*STATES(47)) ALGBRC(147) = CONSTS(310)*(STATES(48)/(STATES(48)+CONSTS(311))) ALGBRC(142) = CONSTS(289) ** CONSTS(290)/(CONSTS(289) ** CONSTS(290)+STATES(48) ** CONSTS(290)) ALGBRC(144) = CONSTS(287)*(STATES(14) ** 2.00000/(CONSTS(288) ** 2.00000+STATES(14) ** 2.00000))*ALGBRC(142) ALGBRC(150) = 1.00000/(1.00000+( CONSTS(317)*CONSTS(318))/CONSTS(317)+STATES(48) ** 2.00000) ALGBRC(146) = ( CONSTS(293)*CONSTS(298)+ CONSTS(295)*STATES(14))/(CONSTS(298)+STATES(14)) ALGBRC(149) = ( (CONSTS(292)+CONSTS(296))*CONSTS(299))/(CONSTS(299)+STATES(14)) ALGBRC(152) = CONSTS(291)*( CONSTS(309)*ALGBRC(146)*STATES(49))/( CONSTS(309)*ALGBRC(146)+ALGBRC(149)) ** 4.00000*(STATES(47) - STATES(14)) ALGBRC(153) = ( CONSTS(294)*CONSTS(299)+ CONSTS(297)*STATES(14))/(CONSTS(299)+STATES(14)) ALGBRC(155) = ( CONSTS(308)*STATES(50))/( (1.00000+CONSTS(304)/STATES(14) ** CONSTS(306))*(1.00000+STATES(14)/CONSTS(305) ** CONSTS(307))) ALGBRC(156) = (( CONSTS(327)*CONSTS(326))/CONSTS(328))*log(CONSTS(331)/STATES(51)) ALGBRC(157) = CONSTS(330)*(STATES(2) - ALGBRC(156)) ALGBRC(165) = CONSTS(339)*(STATES(51) ** CONSTS(341)/(CONSTS(340) ** CONSTS(341)+STATES(51) ** CONSTS(341)))*(CONSTS(342) - STATES(2)) ALGBRC(160) = CONSTS(336)*(CONSTS(335) ** CONSTS(333)/(CONSTS(335) ** CONSTS(333)+STATES(25) ** CONSTS(333))) ALGBRC(161) = ALGBRC(160)*(STATES(2) - CONSTS(334)) ALGBRC(167) = CONSTS(740)*(ALGBRC(161)+ALGBRC(157)+ALGBRC(165)) ALGBRC(51) = CONSTS(696)*(STATES(15) - STATES(25)) ALGBRC(148) = ( CONSTS(322)*EXP( CONSTS(323)*STATES(2)))/(1.00000+EXP( CONSTS(324)*(STATES(2) - CONSTS(325)))) ALGBRC(151) = (( CONSTS(327)*CONSTS(326))/( 2.00000*CONSTS(328)))*log(CONSTS(329)/STATES(25)) ALGBRC(154) = ALGBRC(148)*(STATES(2) - ALGBRC(151)) ALGBRC(158) = CONSTS(332)*(CONSTS(335) ** CONSTS(333)/(CONSTS(335) ** CONSTS(333)+STATES(25) ** CONSTS(333))) ALGBRC(159) = ALGBRC(158)*(STATES(2) - CONSTS(334)) ALGBRC(162) = CONSTS(337)*(STATES(25) ** 2.00000/(CONSTS(338) ** 2.00000+STATES(25) ** 2.00000)) ALGBRC(163) = CONSTS(740)*(ALGBRC(154)+ALGBRC(159)+ALGBRC(162)) ALGBRC(164) = ( CONSTS(343)*(STATES(25) - CONSTS(344)*STATES(52)))/(1.00000+ CONSTS(345)*STATES(25)+ CONSTS(346)*STATES(52)+ CONSTS(347)*STATES(25)*STATES(52)) ALGBRC(171) = CONSTS(371)*(STATES(53)/(STATES(53)+CONSTS(372))) ALGBRC(166) = CONSTS(350) ** CONSTS(351)/(CONSTS(350) ** CONSTS(351)+STATES(53) ** CONSTS(351)) ALGBRC(168) = CONSTS(348)*(STATES(15) ** 2.00000/(CONSTS(349) ** 2.00000+STATES(15) ** 2.00000))*ALGBRC(166) ALGBRC(174) = 1.00000/(1.00000+( CONSTS(378)*CONSTS(379))/CONSTS(378)+STATES(53) ** 2.00000) ALGBRC(170) = ( CONSTS(354)*CONSTS(359)+ CONSTS(356)*STATES(15))/(CONSTS(359)+STATES(15)) ALGBRC(173) = ( (CONSTS(353)+CONSTS(357))*CONSTS(360))/(CONSTS(360)+STATES(15)) ALGBRC(176) = CONSTS(352)*( CONSTS(370)*ALGBRC(170)*STATES(54))/( CONSTS(370)*ALGBRC(170)+ALGBRC(173)) ** 4.00000*(STATES(52) - STATES(15)) ALGBRC(177) = ( CONSTS(355)*CONSTS(360)+ CONSTS(358)*STATES(15))/(CONSTS(360)+STATES(15)) ALGBRC(179) = ( CONSTS(369)*STATES(55))/( (1.00000+CONSTS(365)/STATES(15) ** CONSTS(367))*(1.00000+STATES(15)/CONSTS(366) ** CONSTS(368))) ALGBRC(180) = (( CONSTS(388)*CONSTS(387))/CONSTS(389))*log(CONSTS(392)/STATES(56)) ALGBRC(181) = CONSTS(391)*(STATES(2) - ALGBRC(180)) ALGBRC(189) = CONSTS(400)*(STATES(56) ** CONSTS(402)/(CONSTS(401) ** CONSTS(402)+STATES(56) ** CONSTS(402)))*(CONSTS(403) - STATES(2)) ALGBRC(184) = CONSTS(397)*(CONSTS(396) ** CONSTS(394)/(CONSTS(396) ** CONSTS(394)+STATES(26) ** CONSTS(394))) ALGBRC(185) = ALGBRC(184)*(STATES(2) - CONSTS(395)) ALGBRC(191) = CONSTS(746)*(ALGBRC(185)+ALGBRC(181)+ALGBRC(189)) ALGBRC(52) = CONSTS(697)*(STATES(16) - STATES(26)) ALGBRC(172) = ( CONSTS(383)*EXP( CONSTS(384)*STATES(2)))/(1.00000+EXP( CONSTS(385)*(STATES(2) - CONSTS(386)))) ALGBRC(175) = (( CONSTS(388)*CONSTS(387))/( 2.00000*CONSTS(389)))*log(CONSTS(390)/STATES(26)) ALGBRC(178) = ALGBRC(172)*(STATES(2) - ALGBRC(175)) ALGBRC(182) = CONSTS(393)*(CONSTS(396) ** CONSTS(394)/(CONSTS(396) ** CONSTS(394)+STATES(26) ** CONSTS(394))) ALGBRC(183) = ALGBRC(182)*(STATES(2) - CONSTS(395)) ALGBRC(186) = CONSTS(398)*(STATES(26) ** 2.00000/(CONSTS(399) ** 2.00000+STATES(26) ** 2.00000)) ALGBRC(187) = CONSTS(746)*(ALGBRC(178)+ALGBRC(183)+ALGBRC(186)) ALGBRC(188) = ( CONSTS(404)*(STATES(26) - CONSTS(405)*STATES(57)))/(1.00000+ CONSTS(406)*STATES(26)+ CONSTS(407)*STATES(57)+ CONSTS(408)*STATES(26)*STATES(57)) ALGBRC(195) = CONSTS(432)*(STATES(58)/(STATES(58)+CONSTS(433))) ALGBRC(190) = CONSTS(411) ** CONSTS(412)/(CONSTS(411) ** CONSTS(412)+STATES(58) ** CONSTS(412)) ALGBRC(192) = CONSTS(409)*(STATES(16) ** 2.00000/(CONSTS(410) ** 2.00000+STATES(16) ** 2.00000))*ALGBRC(190) ALGBRC(198) = 1.00000/(1.00000+( CONSTS(439)*CONSTS(440))/CONSTS(439)+STATES(58) ** 2.00000) ALGBRC(194) = ( CONSTS(415)*CONSTS(420)+ CONSTS(417)*STATES(16))/(CONSTS(420)+STATES(16)) ALGBRC(197) = ( (CONSTS(414)+CONSTS(418))*CONSTS(421))/(CONSTS(421)+STATES(16)) ALGBRC(200) = CONSTS(413)*( CONSTS(431)*ALGBRC(194)*STATES(59))/( CONSTS(431)*ALGBRC(194)+ALGBRC(197)) ** 4.00000*(STATES(57) - STATES(16)) ALGBRC(201) = ( CONSTS(416)*CONSTS(421)+ CONSTS(419)*STATES(16))/(CONSTS(421)+STATES(16)) ALGBRC(203) = ( CONSTS(430)*STATES(60))/( (1.00000+CONSTS(426)/STATES(16) ** CONSTS(428))*(1.00000+STATES(16)/CONSTS(427) ** CONSTS(429))) ALGBRC(204) = (( CONSTS(449)*CONSTS(448))/CONSTS(450))*log(CONSTS(453)/STATES(61)) ALGBRC(205) = CONSTS(452)*(STATES(2) - ALGBRC(204)) ALGBRC(213) = CONSTS(461)*(STATES(61) ** CONSTS(463)/(CONSTS(462) ** CONSTS(463)+STATES(61) ** CONSTS(463)))*(CONSTS(464) - STATES(2)) ALGBRC(208) = CONSTS(458)*(CONSTS(457) ** CONSTS(455)/(CONSTS(457) ** CONSTS(455)+STATES(27) ** CONSTS(455))) ALGBRC(209) = ALGBRC(208)*(STATES(2) - CONSTS(456)) ALGBRC(215) = CONSTS(752)*(ALGBRC(209)+ALGBRC(205)+ALGBRC(213)) ALGBRC(53) = CONSTS(698)*(STATES(17) - STATES(27)) ALGBRC(196) = ( CONSTS(444)*EXP( CONSTS(445)*STATES(2)))/(1.00000+EXP( CONSTS(446)*(STATES(2) - CONSTS(447)))) ALGBRC(199) = (( CONSTS(449)*CONSTS(448))/( 2.00000*CONSTS(450)))*log(CONSTS(451)/STATES(27)) ALGBRC(202) = ALGBRC(196)*(STATES(2) - ALGBRC(199)) ALGBRC(206) = CONSTS(454)*(CONSTS(457) ** CONSTS(455)/(CONSTS(457) ** CONSTS(455)+STATES(27) ** CONSTS(455))) ALGBRC(207) = ALGBRC(206)*(STATES(2) - CONSTS(456)) ALGBRC(210) = CONSTS(459)*(STATES(27) ** 2.00000/(CONSTS(460) ** 2.00000+STATES(27) ** 2.00000)) ALGBRC(211) = CONSTS(752)*(ALGBRC(202)+ALGBRC(207)+ALGBRC(210)) ALGBRC(212) = ( CONSTS(465)*(STATES(27) - CONSTS(466)*STATES(62)))/(1.00000+ CONSTS(467)*STATES(27)+ CONSTS(468)*STATES(62)+ CONSTS(469)*STATES(27)*STATES(62)) ALGBRC(219) = CONSTS(493)*(STATES(63)/(STATES(63)+CONSTS(494))) ALGBRC(214) = CONSTS(472) ** CONSTS(473)/(CONSTS(472) ** CONSTS(473)+STATES(63) ** CONSTS(473)) ALGBRC(216) = CONSTS(470)*(STATES(17) ** 2.00000/(CONSTS(471) ** 2.00000+STATES(17) ** 2.00000))*ALGBRC(214) ALGBRC(222) = 1.00000/(1.00000+( CONSTS(500)*CONSTS(501))/CONSTS(500)+STATES(63) ** 2.00000) ALGBRC(218) = ( CONSTS(476)*CONSTS(481)+ CONSTS(478)*STATES(17))/(CONSTS(481)+STATES(17)) ALGBRC(221) = ( (CONSTS(475)+CONSTS(479))*CONSTS(482))/(CONSTS(482)+STATES(17)) ALGBRC(224) = CONSTS(474)*( CONSTS(492)*ALGBRC(218)*STATES(64))/( CONSTS(492)*ALGBRC(218)+ALGBRC(221)) ** 4.00000*(STATES(62) - STATES(17)) ALGBRC(225) = ( CONSTS(477)*CONSTS(482)+ CONSTS(480)*STATES(17))/(CONSTS(482)+STATES(17)) ALGBRC(227) = ( CONSTS(491)*STATES(65))/( (1.00000+CONSTS(487)/STATES(17) ** CONSTS(489))*(1.00000+STATES(17)/CONSTS(488) ** CONSTS(490))) ALGBRC(228) = (( CONSTS(510)*CONSTS(509))/CONSTS(511))*log(CONSTS(514)/STATES(66)) ALGBRC(229) = CONSTS(513)*(STATES(2) - ALGBRC(228)) ALGBRC(237) = CONSTS(522)*(STATES(66) ** CONSTS(524)/(CONSTS(523) ** CONSTS(524)+STATES(66) ** CONSTS(524)))*(CONSTS(525) - STATES(2)) ALGBRC(232) = CONSTS(519)*(CONSTS(518) ** CONSTS(516)/(CONSTS(518) ** CONSTS(516)+STATES(28) ** CONSTS(516))) ALGBRC(233) = ALGBRC(232)*(STATES(2) - CONSTS(517)) ALGBRC(239) = CONSTS(758)*(ALGBRC(233)+ALGBRC(229)+ALGBRC(237)) ALGBRC(54) = CONSTS(699)*(STATES(18) - STATES(28)) ALGBRC(220) = ( CONSTS(505)*EXP( CONSTS(506)*STATES(2)))/(1.00000+EXP( CONSTS(507)*(STATES(2) - CONSTS(508)))) ALGBRC(223) = (( CONSTS(510)*CONSTS(509))/( 2.00000*CONSTS(511)))*log(CONSTS(512)/STATES(28)) ALGBRC(226) = ALGBRC(220)*(STATES(2) - ALGBRC(223)) ALGBRC(230) = CONSTS(515)*(CONSTS(518) ** CONSTS(516)/(CONSTS(518) ** CONSTS(516)+STATES(28) ** CONSTS(516))) ALGBRC(231) = ALGBRC(230)*(STATES(2) - CONSTS(517)) ALGBRC(234) = CONSTS(520)*(STATES(28) ** 2.00000/(CONSTS(521) ** 2.00000+STATES(28) ** 2.00000)) ALGBRC(235) = CONSTS(758)*(ALGBRC(226)+ALGBRC(231)+ALGBRC(234)) ALGBRC(236) = ( CONSTS(526)*(STATES(28) - CONSTS(527)*STATES(67)))/(1.00000+ CONSTS(528)*STATES(28)+ CONSTS(529)*STATES(67)+ CONSTS(530)*STATES(28)*STATES(67)) ALGBRC(243) = CONSTS(554)*(STATES(68)/(STATES(68)+CONSTS(555))) ALGBRC(238) = CONSTS(533) ** CONSTS(534)/(CONSTS(533) ** CONSTS(534)+STATES(68) ** CONSTS(534)) ALGBRC(240) = CONSTS(531)*(STATES(18) ** 2.00000/(CONSTS(532) ** 2.00000+STATES(18) ** 2.00000))*ALGBRC(238) ALGBRC(246) = 1.00000/(1.00000+( CONSTS(561)*CONSTS(562))/CONSTS(561)+STATES(68) ** 2.00000) ALGBRC(242) = ( CONSTS(537)*CONSTS(542)+ CONSTS(539)*STATES(18))/(CONSTS(542)+STATES(18)) ALGBRC(245) = ( (CONSTS(536)+CONSTS(540))*CONSTS(543))/(CONSTS(543)+STATES(18)) ALGBRC(248) = CONSTS(535)*( CONSTS(553)*ALGBRC(242)*STATES(69))/( CONSTS(553)*ALGBRC(242)+ALGBRC(245)) ** 4.00000*(STATES(67) - STATES(18)) ALGBRC(249) = ( CONSTS(538)*CONSTS(543)+ CONSTS(541)*STATES(18))/(CONSTS(543)+STATES(18)) ALGBRC(251) = ( CONSTS(552)*STATES(70))/( (1.00000+CONSTS(548)/STATES(18) ** CONSTS(550))*(1.00000+STATES(18)/CONSTS(549) ** CONSTS(551))) ALGBRC(252) = (( CONSTS(571)*CONSTS(570))/CONSTS(572))*log(CONSTS(575)/STATES(71)) ALGBRC(253) = CONSTS(574)*(STATES(2) - ALGBRC(252)) ALGBRC(261) = CONSTS(583)*(STATES(71) ** CONSTS(585)/(CONSTS(584) ** CONSTS(585)+STATES(71) ** CONSTS(585)))*(CONSTS(586) - STATES(2)) ALGBRC(256) = CONSTS(580)*(CONSTS(579) ** CONSTS(577)/(CONSTS(579) ** CONSTS(577)+STATES(29) ** CONSTS(577))) ALGBRC(257) = ALGBRC(256)*(STATES(2) - CONSTS(578)) ALGBRC(263) = CONSTS(764)*(ALGBRC(257)+ALGBRC(253)+ALGBRC(261)) ALGBRC(55) = CONSTS(700)*(STATES(19) - STATES(29)) ALGBRC(244) = ( CONSTS(566)*EXP( CONSTS(567)*STATES(2)))/(1.00000+EXP( CONSTS(568)*(STATES(2) - CONSTS(569)))) ALGBRC(247) = (( CONSTS(571)*CONSTS(570))/( 2.00000*CONSTS(572)))*log(CONSTS(573)/STATES(29)) ALGBRC(250) = ALGBRC(244)*(STATES(2) - ALGBRC(247)) ALGBRC(254) = CONSTS(576)*(CONSTS(579) ** CONSTS(577)/(CONSTS(579) ** CONSTS(577)+STATES(29) ** CONSTS(577))) ALGBRC(255) = ALGBRC(254)*(STATES(2) - CONSTS(578)) ALGBRC(258) = CONSTS(581)*(STATES(29) ** 2.00000/(CONSTS(582) ** 2.00000+STATES(29) ** 2.00000)) ALGBRC(259) = CONSTS(764)*(ALGBRC(250)+ALGBRC(255)+ALGBRC(258)) ALGBRC(260) = ( CONSTS(587)*(STATES(29) - CONSTS(588)*STATES(72)))/(1.00000+ CONSTS(589)*STATES(29)+ CONSTS(590)*STATES(72)+ CONSTS(591)*STATES(29)*STATES(72)) ALGBRC(267) = CONSTS(615)*(STATES(73)/(STATES(73)+CONSTS(616))) ALGBRC(262) = CONSTS(594) ** CONSTS(595)/(CONSTS(594) ** CONSTS(595)+STATES(73) ** CONSTS(595)) ALGBRC(264) = CONSTS(592)*(STATES(19) ** 2.00000/(CONSTS(593) ** 2.00000+STATES(19) ** 2.00000))*ALGBRC(262) ALGBRC(270) = 1.00000/(1.00000+( CONSTS(622)*CONSTS(623))/CONSTS(622)+STATES(73) ** 2.00000) ALGBRC(266) = ( CONSTS(598)*CONSTS(603)+ CONSTS(600)*STATES(19))/(CONSTS(603)+STATES(19)) ALGBRC(269) = ( (CONSTS(597)+CONSTS(601))*CONSTS(604))/(CONSTS(604)+STATES(19)) ALGBRC(272) = CONSTS(596)*( CONSTS(614)*ALGBRC(266)*STATES(74))/( CONSTS(614)*ALGBRC(266)+ALGBRC(269)) ** 4.00000*(STATES(72) - STATES(19)) ALGBRC(273) = ( CONSTS(599)*CONSTS(604)+ CONSTS(602)*STATES(19))/(CONSTS(604)+STATES(19)) ALGBRC(275) = ( CONSTS(613)*STATES(75))/( (1.00000+CONSTS(609)/STATES(19) ** CONSTS(611))*(1.00000+STATES(19)/CONSTS(610) ** CONSTS(612))) ALGBRC(276) = (( CONSTS(632)*CONSTS(631))/CONSTS(633))*log(CONSTS(636)/STATES(76)) ALGBRC(277) = CONSTS(635)*(STATES(2) - ALGBRC(276)) ALGBRC(284) = CONSTS(644)*(STATES(76) ** CONSTS(646)/(CONSTS(645) ** CONSTS(646)+STATES(76) ** CONSTS(646)))*(CONSTS(647) - STATES(2)) ALGBRC(280) = CONSTS(641)*(CONSTS(640) ** CONSTS(638)/(CONSTS(640) ** CONSTS(638)+STATES(30) ** CONSTS(638))) ALGBRC(281) = ALGBRC(280)*(STATES(2) - CONSTS(639)) ALGBRC(286) = CONSTS(770)*(ALGBRC(281)+ALGBRC(277)+ALGBRC(284)) ALGBRC(30) = CONSTS(16)*STATES(6)*STATES(7)*(STATES(2) - CONSTS(689)) ALGBRC(31) = CONSTS(28)*STATES(8)*(STATES(2) - CONSTS(689)) ALGBRC(32) = CONSTS(33)*(STATES(2) - CONSTS(689)) ALGBRC(33) = CONSTS(34)*STATES(9)*STATES(10)*(STATES(2) - CONSTS(690)) ALGBRC(34) = CONSTS(49)*(STATES(2) - CONSTS(50)) ALGBRC(35) = ALGBRC(8)+ALGBRC(26)+ALGBRC(30)+ALGBRC(31)+ALGBRC(32)+ALGBRC(33)+ALGBRC(34) ALGBRC(73) = ALGBRC(67)+ALGBRC(71) ALGBRC(97) = ALGBRC(91)+ALGBRC(95) ALGBRC(121) = ALGBRC(115)+ALGBRC(119) ALGBRC(145) = ALGBRC(139)+ALGBRC(143) ALGBRC(169) = ALGBRC(163)+ALGBRC(167) ALGBRC(193) = ALGBRC(187)+ALGBRC(191) ALGBRC(217) = ALGBRC(211)+ALGBRC(215) ALGBRC(241) = ALGBRC(235)+ALGBRC(239) ALGBRC(265) = ALGBRC(259)+ALGBRC(263) ALGBRC(268) = ( CONSTS(627)*EXP( CONSTS(628)*STATES(2)))/(1.00000+EXP( CONSTS(629)*(STATES(2) - CONSTS(630)))) ALGBRC(271) = (( CONSTS(632)*CONSTS(631))/( 2.00000*CONSTS(633)))*log(CONSTS(634)/STATES(30)) ALGBRC(274) = ALGBRC(268)*(STATES(2) - ALGBRC(271)) ALGBRC(278) = CONSTS(637)*(CONSTS(640) ** CONSTS(638)/(CONSTS(640) ** CONSTS(638)+STATES(30) ** CONSTS(638))) ALGBRC(279) = ALGBRC(278)*(STATES(2) - CONSTS(639)) ALGBRC(282) = CONSTS(642)*(STATES(30) ** 2.00000/(CONSTS(643) ** 2.00000+STATES(30) ** 2.00000)) ALGBRC(283) = CONSTS(770)*(ALGBRC(274)+ALGBRC(279)+ALGBRC(282)) ALGBRC(288) = ALGBRC(283)+ALGBRC(286) ALGBRC(56) = CONSTS(701)*(STATES(20) - STATES(30)) ALGBRC(285) = ( CONSTS(648)*(STATES(30) - CONSTS(649)*STATES(77)))/(1.00000+ CONSTS(650)*STATES(30)+ CONSTS(651)*STATES(77)+ CONSTS(652)*STATES(30)*STATES(77)) ALGBRC(291) = CONSTS(676)*(STATES(78)/(STATES(78)+CONSTS(677))) ALGBRC(287) = CONSTS(655) ** CONSTS(656)/(CONSTS(655) ** CONSTS(656)+STATES(78) ** CONSTS(656)) ALGBRC(289) = CONSTS(653)*(STATES(20) ** 2.00000/(CONSTS(654) ** 2.00000+STATES(20) ** 2.00000))*ALGBRC(287) ALGBRC(293) = 1.00000/(1.00000+( CONSTS(683)*CONSTS(684))/CONSTS(683)+STATES(78) ** 2.00000) ALGBRC(290) = ( CONSTS(659)*CONSTS(664)+ CONSTS(661)*STATES(20))/(CONSTS(664)+STATES(20)) ALGBRC(292) = ( (CONSTS(658)+CONSTS(662))*CONSTS(665))/(CONSTS(665)+STATES(20)) ALGBRC(294) = CONSTS(657)*( CONSTS(675)*ALGBRC(290)*STATES(79))/( CONSTS(675)*ALGBRC(290)+ALGBRC(292)) ** 4.00000*(STATES(77) - STATES(20)) ALGBRC(295) = ( CONSTS(660)*CONSTS(665)+ CONSTS(663)*STATES(20))/(CONSTS(665)+STATES(20)) ALGBRC(296) = ( CONSTS(674)*STATES(80))/( (1.00000+CONSTS(670)/STATES(20) ** CONSTS(672))*(1.00000+STATES(20)/CONSTS(671) ** CONSTS(673))) ALGBRC(9) = (( CONSTS(53)*CONSTS(52))/( 2.00000*CONSTS(54)))*log(CONSTS(55)/STATES(5)) RETURN END