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