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