C C There are a total of 369 entries in the algebraic variable array. C There are a total of 145 entries in each of the rate and state variable arrays. C There are a total of 410 entries in the constant variable array. C C C CONSTS(325) is AC47_cyt in component ac (uM). C CONSTS(330) is AC47_eca in component ac (uM). C CONSTS(321) is AC56_cav in component ac (uM). C CONSTS(329) is AC56_cyt in component ac (uM). C CONSTS(297) is AC_tot in component ac (uM). C CONSTS(1) is ATP in component ac (uM). C CONSTS(2) is KmATP in component ac (uM). C CONSTS(3) is KmGiAC56 in component ac (uM). C CONSTS(4) is KmGsAC47 in component ac (dimensionless). C CONSTS(5) is KmGsAC56 in component ac (dimensionless). C CONSTS(6) is KmGsGiAC56 in component ac (dimensionless). C ALGBRC(1) is ac_kAC47_cyt_gsa in component ac (dimensionless). C ALGBRC(2) is ac_kAC47_eca_gsa in component ac (dimensionless). C ALGBRC(3) is ac_kAC56_cav_gsa in component ac (dimensionless). C ALGBRC(31) is ac_kAC56_cyt_gsa in component ac (dimensionless). C CONSTS(7) is afAC47 in component ac (hertz). C CONSTS(8) is afAC56 in component ac (hertz). C CONSTS(9) is basalAC47 in component ac (dimensionless). C CONSTS(10) is basalAC56 in component ac (dimensionless). C CONSTS(295) is R_b1_tot in component beta (uM). C STATES(1) is Gi_bg in component beta_cav (uM). C STATES(2) is Gs_aGTP in component beta_cav (uM). C STATES(3) is Gs_aGTP in component beta_cyt (uM). C STATES(4) is Gs_aGTP in component beta_eca (uM). C ALGBRC(69) is dcAMP_AC47_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(68) is dcAMP_AC47_eca in component ac (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(83) is dcAMP_AC56_cav in component ac (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(107) is dcAMP_AC56_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3). C CONSTS(228) is fATP in component ac (dimensionless). C CONSTS(11) is f_AC47_eca in component ac (dimensionless). C CONSTS(261) is f_AC56_AC47 in component ac (dimensionless). C CONSTS(12) is f_AC56_cav in component ac (dimensionless). C ALGBRC(33) is gsi in component ac (dimensionless). C CONSTS(13) is hGsAC47 in component ac (dimensionless). C CONSTS(14) is hGsAC56 in component ac (dimensionless). C CONSTS(15) is hGsGiAC56 in component ac (dimensionless). C ALGBRC(56) is kAC47_cyt in component ac (hertz). C ALGBRC(32) is kAC47_eca in component ac (hertz). C ALGBRC(57) is kAC56_cav in component ac (hertz). C ALGBRC(85) is kAC56_cyt in component ac (hertz). C CONSTS(16) is vGsGiAC56 in component ac (dimensionless). C CONSTS(320) is vr_cav in component cell (dimensionless). C CONSTS(322) is vr_cyt in component cell (dimensionless). C CONSTS(327) is vr_eca in component cell (dimensionless). C CONSTS(17) is ICaL_akap in component akap_sig (uM). C CONSTS(392) is ICaL_akapf in component akap_sig (uM). C CONSTS(393) is ICaL_arn in component akap_sig (uM). C CONSTS(394) is ICaL_arp in component akap_sig (uM). C CONSTS(18) is ICaL_tot in component akap_sig (uM). C CONSTS(282) is ICaLf in component akap_sig (uM). C STATES(5) is ICaLp in component akap_sig (uM). C CONSTS(19) is Ka_ical in component akap_sig (uM). C CONSTS(20) is Ka_ryr in component akap_sig (uM). C CONSTS(21) is Ki in component akap_sig (uM). C CONSTS(22) is Kp_ical in component akap_sig (uM). C CONSTS(23) is Kp_ryr in component akap_sig (uM). C CONSTS(24) is Kr in component akap_sig (uM). C CONSTS(25) is Li in component akap_sig (uM). C CONSTS(26) is Lr in component akap_sig (uM). C CONSTS(27) is Mi in component akap_sig (uM). C CONSTS(28) is Mr in component akap_sig (uM). C CONSTS(381) is PKA_cav in component pka (uM). C CONSTS(391) is PKAf in component akap_sig (uM). C CONSTS(29) is PP1_cav in component pp1 (uM). C CONSTS(292) is PP1f_cav in component akap_sig (uM). C CONSTS(30) is RyR_akap in component akap_sig (uM). C CONSTS(395) is RyR_akapf in component akap_sig (uM). C CONSTS(396) is RyR_arn in component akap_sig (uM). C CONSTS(397) is RyR_arp in component akap_sig (uM). C CONSTS(31) is RyR_tot in component akap_sig (uM). C CONSTS(294) is RyRf in component akap_sig (uM). C STATES(6) is RyRp in component akap_sig (uM). C CONSTS(274) is akap_sig_ICaLf_sum in component akap_sig (dimensionless). C ALGBRC(4) is akap_sig_ICaLp_dif in component akap_sig (uM). C CONSTS(388) is akap_sig_PKAf_arg in component akap_sig (dimensionless). C CONSTS(382) is akap_sig_PKAf_b in component akap_sig (uM). C CONSTS(383) is akap_sig_PKAf_c in component akap_sig (mM2_times_1e_minus_6). C CONSTS(384) is akap_sig_PKAf_d in component akap_sig (mM3_times_1e_minus_9). C CONSTS(389) is akap_sig_PKAf_mag in component akap_sig (uM). C CONSTS(385) is akap_sig_PKAf_rr in component akap_sig (mol6_per_m18_times_1e_minus_18). C CONSTS(390) is akap_sig_PKAf_x in component akap_sig (dimensionless). C CONSTS(386) is akap_sig_PKAf_yi in component akap_sig (mM3_times_1e_minus_9). C CONSTS(387) is akap_sig_PKAf_yr in component akap_sig (mM3_times_1e_minus_9). C CONSTS(289) is akap_sig_PP1f_cav_arg in component akap_sig (dimensionless). C CONSTS(283) is akap_sig_PP1f_cav_b in component akap_sig (uM). C CONSTS(284) is akap_sig_PP1f_cav_c in component akap_sig (mM2_times_1e_minus_6). C CONSTS(285) is akap_sig_PP1f_cav_d in component akap_sig (mM3_times_1e_minus_9). C CONSTS(290) is akap_sig_PP1f_cav_mag in component akap_sig (uM). C CONSTS(286) is akap_sig_PP1f_cav_rr in component akap_sig (mol6_per_m18_times_1e_minus_18). C CONSTS(291) is akap_sig_PP1f_cav_x in component akap_sig (dimensionless). C CONSTS(287) is akap_sig_PP1f_cav_yi in component akap_sig (mM3_times_1e_minus_9). C CONSTS(288) is akap_sig_PP1f_cav_yr in component akap_sig (mM3_times_1e_minus_9). C CONSTS(293) is akap_sig_RyRf_sum in component akap_sig (dimensionless). C ALGBRC(5) is akap_sig_RyRp_dif in component akap_sig (uM). C ALGBRC(11) is fp_ICaL in component akap_sig (dimensionless). C ALGBRC(39) is fp_RyR in component akap_sig (dimensionless). C CONSTS(32) is ka_ical in component akap_sig (hertz). C CONSTS(33) is ka_ryr in component akap_sig (hertz). C CONSTS(34) is kp_ical in component akap_sig (hertz). C CONSTS(35) is kp_ryr in component akap_sig (hertz). C STATES(7) is C in component pka_cav (uM). C VOI is time in component engine (ms). C CONSTS(36) is Gi_tot in component beta (uM). C CONSTS(296) is Gs_tot in component beta (uM). C CONSTS(298) is R_b2_tot in component beta (uM). C CONSTS(37) is f_Gi_cav in component beta (dimensionless). C CONSTS(216) is f_Gi_eca in component beta (dimensionless). C CONSTS(38) is f_Gs_cav in component beta (dimensionless). C CONSTS(217) is f_Gs_cyt in component beta (dimensionless). C CONSTS(39) is f_Gs_eca in component beta (dimensionless). C CONSTS(40) is f_Rb1_cav in component beta (dimensionless). C CONSTS(299) is f_Rb1_cyt in component beta (dimensionless). C CONSTS(41) is f_Rb1_eca in component beta (dimensionless). C CONSTS(42) is f_Rb2_cav in component beta (dimensionless). C CONSTS(300) is f_Rb2_eca in component beta (dimensionless). C CONSTS(43) is k_act1_Gi in component beta (hertz). C CONSTS(44) is k_act1_Gs in component beta (hertz). C CONSTS(45) is k_act2_Gi in component beta (hertz). C CONSTS(46) is k_act2_Gs in component beta (hertz). C CONSTS(47) is k_b1_c in component beta (uM). C CONSTS(48) is k_b1_h in component beta (uM). C CONSTS(49) is k_b1_l in component beta (uM). C CONSTS(50) is k_b2_a in component beta (uM). C CONSTS(51) is k_b2_c in component beta (uM). C CONSTS(52) is k_b2_f in component beta (uM). C CONSTS(53) is k_b2_h in component beta (uM). C CONSTS(54) is k_b2_l in component beta (uM). C CONSTS(55) is k_b2_n in component beta (uM). C CONSTS(218) is k_grk_dp in component beta (hertz). C CONSTS(219) is k_grk_p in component beta (hertz). C CONSTS(220) is k_hydr_Gi in component beta (hertz). C CONSTS(56) is k_hydr_Gs in component beta (hertz). C CONSTS(258) is k_pka_dp in component beta (hertz). C CONSTS(221) is k_pka_p in component beta (per_mM_per_ms). C CONSTS(222) is k_reas_Gi in component beta (per_mM_per_ms). C CONSTS(57) is k_reas_Gs in component beta (per_mM_per_ms). C CONSTS(58) is rate_bds in component beta (per_mM_per_ms). C CONSTS(59) is GRK in component beta_cav (dimensionless). C STATES(8) is Gi_aGDP in component beta_cav (uM). C STATES(9) is Gi_aGTP in component beta_cav (uM). C ALGBRC(6) is Gi_abg in component beta_cav (uM). C ALGBRC(87) is Gi_f in component beta_cav (uM). C STATES(10) is Gs_aGDP in component beta_cav (uM). C ALGBRC(7) is Gs_abg in component beta_cav (uM). C STATES(11) is Gs_bg in component beta_cav (uM). C ALGBRC(157) is Gs_f in component beta_cav (uM). C ALGBRC(189) is LRGs_tot in component beta_cav (uM). C ALGBRC(165) is LRb1 in component beta_cav (uM). C ALGBRC(169) is LRb1Gs in component beta_cav (uM). C ALGBRC(181) is LRb2 in component beta_cav (uM). C ALGBRC(109) is LRb2Gi in component beta_cav (uM). C ALGBRC(185) is LRb2Gs in component beta_cav (uM). C ALGBRC(197) is RGs_tot in component beta_cav (uM). C CONSTS(323) is R_b1_tot in component beta_cav (uM). C CONSTS(324) is R_b2_tot in component beta_cav (uM). C ALGBRC(173) is Rb1Gs in component beta_cav (uM). C ALGBRC(161) is Rb1_f in component beta_cav (uM). C STATES(12) is Rb1_grk_tot in component beta_cav (uM). C ALGBRC(34) is Rb1_np_tot in component beta_cav (uM). C STATES(13) is Rb1_pka_tot in component beta_cav (uM). C ALGBRC(99) is Rb2Gi in component beta_cav (uM). C ALGBRC(193) is Rb2Gs in component beta_cav (uM). C ALGBRC(177) is Rb2_f in component beta_cav (uM). C STATES(14) is Rb2_grk_tot in component beta_cav (uM). C ALGBRC(58) is Rb2_np_tot in component beta_cav (uM). C ALGBRC(71) is Rb2_pka_f in component beta_cav (uM). C STATES(15) is Rb2_pka_tot in component beta_cav (uM). C CONSTS(224) is beta_cav_Gs_f_a in component beta_cav (mol4_per_m12_times_1e_minus_12). C ALGBRC(132) is beta_cav_Gs_f_arg in component beta_cav (dimensionless). C ALGBRC(70) is beta_cav_Gs_f_b in component beta_cav (uM). C ALGBRC(86) is beta_cav_Gs_f_c in component beta_cav (mM2_times_1e_minus_6). C CONSTS(259) is beta_cav_Gs_f_c11 in component beta_cav (mol5_per_m15_times_1e_minus_15). C CONSTS(272) is beta_cav_Gs_f_c22 in component beta_cav (mol5_per_m15_times_1e_minus_15). C CONSTS(280) is beta_cav_Gs_f_c33 in component beta_cav (mol6_per_m18_times_1e_minus_18). C ALGBRC(98) is beta_cav_Gs_f_d in component beta_cav (mM3_times_1e_minus_9). C ALGBRC(148) is beta_cav_Gs_f_i in component beta_cav (uM). C ALGBRC(138) is beta_cav_Gs_f_mag in component beta_cav (uM). C ALGBRC(152) is beta_cav_Gs_f_r in component beta_cav (uM). C ALGBRC(110) is beta_cav_Gs_f_rr in component beta_cav (mol6_per_m18_times_1e_minus_18). C ALGBRC(143) is beta_cav_Gs_f_x in component beta_cav (dimensionless). C ALGBRC(118) is beta_cav_Gs_f_yi in component beta_cav (mM3_times_1e_minus_9). C ALGBRC(126) is beta_cav_Gs_f_yr in component beta_cav (mM3_times_1e_minus_9). C CONSTS(223) is beta_cav_Rb2_pka_f_a in component beta_cav (uM). C ALGBRC(35) is beta_cav_Rb2_pka_f_b in component beta_cav (mM2_times_1e_minus_6). C ALGBRC(59) is beta_cav_Rb2_pka_f_c in component beta_cav (mM3_times_1e_minus_9). C CONSTS(60) is L in component iso (uM). C CONSTS(61) is k_GsAct_b2 in component beta_cav (dimensionless). C CONSTS(62) is GRK in component beta_cyt (dimensionless). C STATES(16) is Gs_aGDP in component beta_cyt (uM). C ALGBRC(8) is Gs_abg in component beta_cyt (uM). C STATES(17) is Gs_bg in component beta_cyt (uM). C ALGBRC(100) is Gs_f in component beta_cyt (uM). C ALGBRC(111) is LRb1Gs_np in component beta_cyt (uM). C ALGBRC(119) is LRb1_np in component beta_cyt (uM). C CONSTS(326) is R_b1_tot in component beta_cyt (uM). C ALGBRC(120) is Rb1Gs_np in component beta_cyt (uM). C STATES(18) is Rb1_grk_tot in component beta_cyt (uM). C ALGBRC(88) is Rb1_np_f in component beta_cyt (uM). C ALGBRC(36) is Rb1_np_tot in component beta_cyt (uM). C STATES(19) is Rb1_pka_tot in component beta_cyt (uM). C CONSTS(225) is beta_cyt_Rb1_np_f_a in component beta_cyt (uM). C ALGBRC(60) is beta_cyt_Rb1_np_f_b in component beta_cyt (mM2_times_1e_minus_6). C ALGBRC(72) is beta_cyt_Rb1_np_f_c in component beta_cyt (mM3_times_1e_minus_9). C STATES(20) is C in component pka_cyt (uM). C CONSTS(63) is GRK in component beta_eca (dimensionless). C STATES(21) is Gi_aGDP in component beta_eca (uM). C STATES(22) is Gi_aGTP in component beta_eca (uM). C ALGBRC(9) is Gi_abg in component beta_eca (uM). C STATES(23) is Gi_bg in component beta_eca (uM). C ALGBRC(90) is Gi_f in component beta_eca (uM). C STATES(24) is Gs_aGDP in component beta_eca (uM). C ALGBRC(10) is Gs_abg in component beta_eca (uM). C STATES(25) is Gs_bg in component beta_eca (uM). C ALGBRC(158) is Gs_f in component beta_eca (uM). C ALGBRC(190) is LRGs_tot in component beta_eca (uM). C ALGBRC(166) is LRb1 in component beta_eca (uM). C ALGBRC(170) is LRb1Gs in component beta_eca (uM). C ALGBRC(182) is LRb2 in component beta_eca (uM). C ALGBRC(112) is LRb2Gi in component beta_eca (uM). C ALGBRC(186) is LRb2Gs in component beta_eca (uM). C ALGBRC(198) is RGs_tot in component beta_eca (uM). C CONSTS(331) is R_b1_tot in component beta_eca (uM). C CONSTS(332) is R_b2_tot in component beta_eca (uM). C ALGBRC(174) is Rb1Gs in component beta_eca (uM). C ALGBRC(162) is Rb1_f in component beta_eca (uM). C STATES(26) is Rb1_grk_tot in component beta_eca (uM). C ALGBRC(37) is Rb1_np_tot in component beta_eca (uM). C STATES(27) is Rb1_pka_tot in component beta_eca (uM). C ALGBRC(102) is Rb2Gi in component beta_eca (uM). C ALGBRC(194) is Rb2Gs in component beta_eca (uM). C ALGBRC(178) is Rb2_f in component beta_eca (uM). C STATES(28) is Rb2_grk_tot in component beta_eca (uM). C ALGBRC(61) is Rb2_np_tot in component beta_eca (uM). C ALGBRC(74) is Rb2_pka_f in component beta_eca (uM). C STATES(29) is Rb2_pka_tot in component beta_eca (uM). C CONSTS(227) is beta_eca_Gs_f_a in component beta_eca (mol4_per_m12_times_1e_minus_12). C ALGBRC(133) is beta_eca_Gs_f_arg in component beta_eca (dimensionless). C ALGBRC(73) is beta_eca_Gs_f_b in component beta_eca (uM). C ALGBRC(89) is beta_eca_Gs_f_c in component beta_eca (mM2_times_1e_minus_6). C CONSTS(260) is beta_eca_Gs_f_c11 in component beta_eca (mol5_per_m15_times_1e_minus_15). C CONSTS(273) is beta_eca_Gs_f_c22 in component beta_eca (mol5_per_m15_times_1e_minus_15). C CONSTS(281) is beta_eca_Gs_f_c33 in component beta_eca (mol6_per_m18_times_1e_minus_18). C ALGBRC(101) is beta_eca_Gs_f_d in component beta_eca (mM3_times_1e_minus_9). C ALGBRC(149) is beta_eca_Gs_f_i in component beta_eca (uM). C ALGBRC(139) is beta_eca_Gs_f_mag in component beta_eca (uM). C ALGBRC(153) is beta_eca_Gs_f_r in component beta_eca (uM). C ALGBRC(113) is beta_eca_Gs_f_rr in component beta_eca (mol6_per_m18_times_1e_minus_18). C ALGBRC(144) is beta_eca_Gs_f_x in component beta_eca (dimensionless). C ALGBRC(121) is beta_eca_Gs_f_yi in component beta_eca (mM3_times_1e_minus_9). C ALGBRC(127) is beta_eca_Gs_f_yr in component beta_eca (mM3_times_1e_minus_9). C CONSTS(226) is beta_eca_Rb2_pka_f_a in component beta_eca (uM). C ALGBRC(38) is beta_eca_Rb2_pka_f_b in component beta_eca (mM2_times_1e_minus_6). C ALGBRC(62) is beta_eca_Rb2_pka_f_c in component beta_eca (mM3_times_1e_minus_9). C CONSTS(64) is k_GsAct_b2 in component beta_eca (dimensionless). C STATES(30) is C in component pka_eca (uM). C CONSTS(303) is AF in component cell (m2_mol_per_s_per_A_times_1e_minus_4). C ALGBRC(187) is Ca in component calcium (mM). C ALGBRC(103) is Ca_CaL in component calcium (mM). C ALGBRC(128) is Ca_jsr in component calcium (mM). C STATES(31) is Ca_nsr in component calcium (mM). C ALGBRC(150) is Ca_sr in component calcium (mM). C ALGBRC(306) is ICa_tot in component calcium (uA_per_cm2). C ALGBRC(218) is ICab in component icab (uA_per_cm2). C ALGBRC(298) is INaCa in component inaca (uA_per_cm2). C ALGBRC(292) is INaCaSR in component inaca (uA_per_cm2). C ALGBRC(202) is Idiff_Ca in component diff (mM_per_ms). C ALGBRC(210) is Idiff_sr in component diff (mM_per_ms). C ALGBRC(305) is IpCa in component ipca (uA_per_cm2). C ALGBRC(320) is Irel in component irel (mM_per_ms). C ALGBRC(212) is Itr in component diff (mM_per_ms). C ALGBRC(329) is Iup in component iup (mM_per_ms). C CONSTS(65) is Ka_tni in component calcium (uM). C CONSTS(66) is Kp_tni in component calcium (uM). C CONSTS(67) is PP2A in component pp1 (uM). C CONSTS(229) is bar_sum in component calcium (mM). C CONSTS(68) is bsl_bar in component calcium (mM). C CONSTS(69) is bsl_km in component calcium (mM). C CONSTS(70) is bsr_bar in component calcium (mM). C CONSTS(71) is bsr_km in component calcium (mM). C ALGBRC(63) is calcium_Ca_CaL_b in component calcium (mM). C ALGBRC(75) is calcium_Ca_CaL_c in component calcium (mM2). C ALGBRC(91) is calcium_Ca_CaL_d in component calcium (mM3). C ALGBRC(179) is calcium_Ca_b in component calcium (mM). C ALGBRC(183) is calcium_Ca_c in component calcium (mM2). C ALGBRC(171) is calcium_Ca_d in component calcium (mM3). C ALGBRC(114) is calcium_Ca_jsr_b in component calcium (mM). C ALGBRC(122) is calcium_Ca_jsr_c in component calcium (mM2). C CONSTS(313) is calcium_Ca_nsr_r1 in component calcium (dimensionless). C ALGBRC(134) is calcium_Ca_sr_b in component calcium (mM). C ALGBRC(140) is calcium_Ca_sr_c in component calcium (mM2). C ALGBRC(145) is calcium_Ca_sr_d in component calcium (mM3). C ALGBRC(154) is calcium_fhat_val in component calcium (dimensionless). C CONSTS(306) is calcium_uCa_CaL_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(316) is calcium_uCa_CaL_r2 in component calcium (dimensionless). C CONSTS(309) is calcium_uCa_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(314) is calcium_uCa_r2 in component calcium (dimensionless). C CONSTS(317) is calcium_uCa_sr_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(318) is calcium_uCa_sr_r2 in component calcium (dimensionless). C CONSTS(72) is cbar in component calcium (mM). C CONSTS(73) is csqn_bar in component calcium (mM). C CONSTS(74) is csqn_km in component calcium (mM). C STATES(32) is f_tni in component calcium (dimensionless). C ALGBRC(159) is fhat in component calcium (dimensionless). C ALGBRC(249) is ICaL in component ical (uA_per_cm2). C CONSTS(75) is ka_tni in component calcium (hertz). C CONSTS(76) is kc in component calcium (mM). C CONSTS(231) is km_pro in component calcium (mM2). C CONSTS(262) is km_sum in component calcium (mM). C CONSTS(77) is kp_tni in component calcium (hertz). C ALGBRC(167) is kpro in component calcium (mM2). C ALGBRC(175) is ksum in component calcium (mM). C ALGBRC(163) is kt in component calcium (mM). C CONSTS(78) is ktn in component calcium (mM). C CONSTS(230) is ktp in component calcium (mM). C CONSTS(319) is r3 in component calcium (dimensionless). C CONSTS(263) is ss_pro in component calcium (mM2). C CONSTS(275) is ss_sum in component calcium (mM). C CONSTS(79) is tbar in component calcium (mM). C STATES(33) is uCa in component calcium (mM). C STATES(34) is uCa_CaL in component calcium (mM). C STATES(35) is uCa_jsr in component calcium (mM). C STATES(36) is uCa_sr in component calcium (mM). C CONSTS(305) is v_CaL in component cell (uL). C CONSTS(308) is v_cyt in component cell (uL). C CONSTS(311) is v_jsr in component cell (uL). C CONSTS(312) is v_nsr in component cell (uL). C CONSTS(315) is v_sr in component cell (uL). C CONSTS(80) is CaMK0 in component camk (dimensionless). C CONSTS(81) is K in component camk (dimensionless). C CONSTS(82) is Km in component camk (mM). C CONSTS(83) is PP1_eca in component pp1 (uM). C ALGBRC(207) is PP1_tot in component camk (uM). C ALGBRC(204) is PP1f_cyt in component pp1 (uM). C ALGBRC(195) is active in component camk (dimensionless). C ALGBRC(191) is bound in component camk (dimensionless). C ALGBRC(200) is c in component camk (dimensionless). C ALGBRC(201) is camk_f_ryr_d in component camk (dimensionless). C CONSTS(84) is camk_trap_alpha in component camk (mS_per_uF). C CONSTS(85) is camk_trap_beta in component camk (mS_per_uF). C STATES(37) is f_ical in component camk (dimensionless). C STATES(38) is f_ik1 in component camk (dimensionless). C STATES(39) is f_ina in component camk (dimensionless). C STATES(40) is f_ito in component camk (dimensionless). C STATES(41) is f_plb in component camk (dimensionless). C STATES(42) is f_ryr in component camk (dimensionless). C CONSTS(232) is tau_cal in component camk (ms). C CONSTS(233) is tau_ik1 in component camk (ms). C CONSTS(234) is tau_ina in component camk (ms). C CONSTS(235) is tau_ito in component camk (ms). C CONSTS(86) is tau_plb in component camk (ms). C CONSTS(87) is tau_ryr in component camk (ms). C STATES(43) is trap in component camk (dimensionless). C STATES(44) is cAMP_cav in component camp (uM). C STATES(45) is cAMP_cyt in component camp (uM). C STATES(46) is cAMP_eca in component camp (uM). C ALGBRC(96) is camp_cAMP_cav_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(106) is camp_cAMP_cav_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(136) is camp_cAMP_cav_pde in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(116) is camp_cAMP_cyt_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(124) is camp_cAMP_cyt_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(147) is camp_cAMP_cyt_pde in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(84) is camp_cAMP_eca_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(97) is camp_cAMP_eca_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(125) is camp_cAMP_eca_pde in component camp (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(115) is dcAMP_PDE2_cav in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(130) is dcAMP_PDE2_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(108) is dcAMP_PDE2_eca in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(123) is dcAMP_PDE3_cav in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(137) is dcAMP_PDE3_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(129) is dcAMP_PDE4_cav in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(142) is dcAMP_PDE4_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(117) is dcAMP_PDE4_eca in component pde (mol_per_m3_per_s_times_1e_minus_3). C CONSTS(236) is j_cav_cyt in component camp (m3_per_s_times_1e_minus_9). C CONSTS(237) is j_cav_eca in component camp (m3_per_s_times_1e_minus_9). C CONSTS(238) is j_eca_cyt in component camp (m3_per_s_times_1e_minus_9). C ALGBRC(146) is dcAMP in component pka_cav (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(155) is dcAMP in component pka_cyt (mol_per_m3_per_s_times_1e_minus_3). C ALGBRC(135) is dcAMP in component pka_eca (mol_per_m3_per_s_times_1e_minus_3). C CONSTS(307) is v_cav in component cell (uL). C CONSTS(310) is v_eca in component cell (uL). C CONSTS(88) is F in component phys (C_per_mol). C CONSTS(302) is capArea in component cell (cm2). C CONSTS(301) is geoArea in component cell (cm2). C CONSTS(89) is length in component cell (cm). C CONSTS(90) is pi in component cell (dimensionless). C CONSTS(91) is radius in component cell (cm). C CONSTS(304) is volume in component cell (uL). C ALGBRC(334) is CTKCl in component ctkcl (mM_per_ms). C ALGBRC(355) is CTNaCl in component ctnacl (mM_per_ms). C STATES(47) is Cl in component chloride (mM). C STATES(48) is Cl_sr in component chloride (mM). C ALGBRC(324) is IClCa in component iclca (uA_per_cm2). C ALGBRC(331) is ICl_tot in component chloride (uA_per_cm2). C ALGBRC(330) is IClb in component iclb (uA_per_cm2). C ALGBRC(205) is Idiff_Cl in component diff (mM_per_ms). C CONSTS(333) is chloride_Cl_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5). C CONSTS(338) is chloride_Cl_r2 in component chloride (dimensionless). C CONSTS(342) is chloride_Cl_sr_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5). C ALGBRC(328) is ECl in component nernst (mV). C ALGBRC(332) is EK in component nernst (mV). C CONSTS(92) is KClBar in component ctkcl (mM_per_ms). C ALGBRC(333) is ctkcl_CTKCl_z1 in component ctkcl (mV). C CONSTS(93) is ctkcl_CTKCl_z2 in component ctkcl (mV). C ALGBRC(353) is ENa in component nernst (mV). C CONSTS(94) is NaClBar in component ctnacl (mM_per_ms). C ALGBRC(354) is ctnacl_CTNaCl_z1 in component ctnacl (g4_m8_per_s12_per_A4). C CONSTS(347) is ctnacl_CTNaCl_z2 in component ctnacl (g4_m8_per_s12_per_A4). C ALGBRC(208) is Idiff_Na in component diff (mM_per_ms). C STATES(49) is Na in component sodium (mM). C STATES(50) is Na_sr in component sodium (mM). C CONSTS(95) is tau in component diff (ms). C CONSTS(96) is tau_sr in component diff (ms). C CONSTS(97) is tau_tr in component diff (ms). C CONSTS(98) is Cao in component extra (mM). C CONSTS(99) is Clo in component extra (mM). C CONSTS(100) is Ko in component extra (mM). C CONSTS(101) is Nao in component extra (mM). C CONSTS(264) is FRT in component phys (per_mV). C STATES(51) is V in component membrane (mV). C ALGBRC(216) is efrt in component icab (dimensionless). C CONSTS(102) is pCab in component icab (cm_per_s). C ALGBRC(214) is vfrt in component icab (dimensionless). C ALGBRC(224) is f_hat in component ical (dimensionless). C ALGBRC(236) is ICaL in component ical_camk (uA_per_cm2). C CONSTS(398) is ical_f_hat_ratio in component ical (dimensionless). C ALGBRC(221) is ical_f_hat_val in component ical (dimensionless). C ALGBRC(247) is ICaL in component ical_np (uA_per_cm2). C STATES(52) is C in component ical_camk (dimensionless). C STATES(53) is CI in component ical_camk (dimensionless). C STATES(54) is CIs in component ical_camk (dimensionless). C STATES(55) is Cs in component ical_camk (dimensionless). C CONSTS(276) is FFRT in component phys (s4_A2_per_g_per_m2_per_mol). C ALGBRC(233) is IBar in component ical_camk (uA_per_cm2). C STATES(56) is O in component ical_camk (dimensionless). C STATES(57) is OI in component ical_camk (dimensionless). C STATES(58) is OIs in component ical_camk (dimensionless). C STATES(59) is Os in component ical_camk (dimensionless). C ALGBRC(227) is PCa in component ical_camk (cm_per_s). C ALGBRC(156) is ac_inf in component ical_camk (dimensionless). C ALGBRC(196) is ac_tau in component ical_np (ms). C ALGBRC(203) is alpha in component ical_camk (mS_per_uF). C ALGBRC(206) is beta in component ical_camk (mS_per_uF). C ALGBRC(184) is delta in component ical_camk (mS_per_uF). C ALGBRC(243) is delta1 in component ical_camk (mS_per_uF). C ALGBRC(160) is delta_tau in component ical_camk (dimensionless). C ALGBRC(230) is ical_camk_IBar_vv in component ical_camk (dimensionless). C ALGBRC(231) is ical_camk_delta1_xs_cor in component ical_camk (mS_per_uF). C ALGBRC(237) is ical_camk_delta1_y_cor in component ical_camk (mS_per_uF). C ALGBRC(215) is in_a in component ical_np (dimensionless). C ALGBRC(217) is in_b in component ical_np (dimensionless). C ALGBRC(168) is in_hi_inf in component ical_camk (dimensionless). C ALGBRC(219) is in_hi_tau in component ical_camk (ms). C ALGBRC(164) is in_inf in component ical_camk (dimensionless). C ALGBRC(172) is in_lo_inf in component ical_camk (dimensionless). C ALGBRC(222) is in_lo_tau in component ical_camk (ms). C ALGBRC(188) is inca in component ical_camk (dimensionless). C ALGBRC(176) is ss_cal_10 in component ical_camk (dimensionless). C ALGBRC(180) is ss_cal_4 in component ical_camk (dimensionless). C CONSTS(103) is theta in component ical_camk (mS_per_uF). C CONSTS(104) is theta1 in component ical_camk (mS_per_uF). C ALGBRC(225) is x in component ical_camk (mS_per_uF). C ALGBRC(228) is xs in component ical_camk (mS_per_uF). C ALGBRC(234) is y in component ical_camk (mS_per_uF). C ALGBRC(240) is ys in component ical_camk (mS_per_uF). C STATES(60) is C in component ical_np (dimensionless). C STATES(61) is CI in component ical_np (dimensionless). C STATES(62) is CIs in component ical_np (dimensionless). C STATES(63) is Cs in component ical_np (dimensionless). C ALGBRC(245) is IBar in component ical_np (uA_per_cm2). C STATES(64) is O in component ical_np (dimensionless). C STATES(65) is OI in component ical_np (dimensionless). C STATES(66) is OIs in component ical_np (dimensionless). C STATES(67) is Os in component ical_np (dimensionless). C ALGBRC(239) is PCa in component ical_np (cm_per_s). C ALGBRC(192) is ac_inf in component ical_np (dimensionless). C ALGBRC(209) is alpha in component ical_np (mS_per_uF). C ALGBRC(211) is beta in component ical_np (mS_per_uF). C ALGBRC(238) is delta in component ical_np (mS_per_uF). C ALGBRC(261) is delta1 in component ical_np (mS_per_uF). C ALGBRC(213) is delta_tau in component ical_np (dimensionless). C ALGBRC(242) is ical_np_IBar_vv in component ical_np (dimensionless). C ALGBRC(250) is ical_np_delta1_xs_cor in component ical_np (mS_per_uF). C ALGBRC(255) is ical_np_delta1_y_cor in component ical_np (mS_per_uF). C ALGBRC(223) is in_hi_inf in component ical_np (dimensionless). C ALGBRC(244) is in_hi_tau in component ical_np (ms). C ALGBRC(220) is in_inf in component ical_np (dimensionless). C ALGBRC(226) is in_lo_inf in component ical_np (dimensionless). C ALGBRC(229) is in_lo_tau in component ical_np (ms). C ALGBRC(241) is inca in component ical_np (dimensionless). C ALGBRC(232) is ss_cal_10 in component ical_np (dimensionless). C ALGBRC(235) is ss_cal_4 in component ical_np (dimensionless). C CONSTS(105) is theta in component ical_np (mS_per_uF). C CONSTS(106) is theta1 in component ical_np (mS_per_uF). C ALGBRC(246) is x in component ical_np (mS_per_uF). C ALGBRC(248) is xs in component ical_np (mS_per_uF). C ALGBRC(252) is y in component ical_np (mS_per_uF). C ALGBRC(258) is ys in component ical_np (mS_per_uF). C CONSTS(107) is Gbar in component iclb (mS_per_cm2). C ALGBRC(254) is IClCa_bar in component iclca (uA_per_cm2). C ALGBRC(318) is Irel_pure in component irel (mM_per_ms). C ALGBRC(321) is KClCa in component iclca (dimensionless). C CONSTS(108) is PCl in component iclca (cm_per_s). C STATES(68) is i2 in component iclca (dimensionless). C ALGBRC(12) is iclca_i2_alpha in component iclca (dimensionless). C ALGBRC(40) is iclca_i2_beta in component iclca (dimensionless). C CONSTS(109) is kCaCl in component iclca (mM_per_ms). C CONSTS(110) is tau in component iclca (ms). C ALGBRC(251) is vexp in component iclca (dimensionless). C CONSTS(240) is Gbar in component ik1 (mS_per_cm2). C ALGBRC(340) is IK1 in component ik1 (uA_per_cm2). C ALGBRC(339) is IK1_camk in component ik1 (uA_per_cm2). C ALGBRC(338) is IK1_np in component ik1 (uA_per_cm2). C ALGBRC(336) is ik1_IK1_np_alpha in component ik1 (dimensionless). C ALGBRC(337) is ik1_IK1_np_beta in component ik1 (dimensionless). C ALGBRC(335) is ik1_IK1_np_vv in component ik1 (mV). C CONSTS(241) is GKr in component ikr (mS_per_cm2). C ALGBRC(341) is IKr in component ikr (uA_per_cm2). C STATES(69) is ac in component ikr (dimensionless). C ALGBRC(13) is ikr_ac_tau in component ikr (ms). C ALGBRC(41) is inf in component ikr (dimensionless). C ALGBRC(257) is inx in component ikr (dimensionless). C ALGBRC(349) is EKs in component nernst (mV). C ALGBRC(260) is G in component iks (mS_per_cm2). C ALGBRC(352) is IKs in component iks (uA_per_cm2). C CONSTS(404) is IKs_arn in component iks_sig (uM). C ALGBRC(350) is IKs_np in component iks (uA_per_cm2). C ALGBRC(351) is IKs_pka in component iks (uA_per_cm2). C CONSTS(111) is IKs_tot in component iks_sig (uM). C ALGBRC(267) is f_hat in component iks (dimensionless). C ALGBRC(263) is fp_iks in component iks_sig (dimensionless). C CONSTS(405) is iks_f_hat_ratio in component iks (dimensionless). C ALGBRC(265) is iks_f_hat_val in component iks (dimensionless). C STATES(70) is O1 in component iks_np (dimensionless). C STATES(71) is O2 in component iks_np (dimensionless). C STATES(72) is O1 in component iks_pka (dimensionless). C STATES(73) is O2 in component iks_pka (dimensionless). C STATES(74) is C1 in component iks_np (dimensionless). C STATES(75) is C10 in component iks_np (dimensionless). C STATES(76) is C11 in component iks_np (dimensionless). C STATES(77) is C12 in component iks_np (dimensionless). C STATES(78) is C13 in component iks_np (dimensionless). C STATES(79) is C14 in component iks_np (dimensionless). C STATES(80) is C15 in component iks_np (dimensionless). C STATES(81) is C2 in component iks_np (dimensionless). C STATES(82) is C3 in component iks_np (dimensionless). C STATES(83) is C4 in component iks_np (dimensionless). C STATES(84) is C5 in component iks_np (dimensionless). C STATES(85) is C6 in component iks_np (dimensionless). C STATES(86) is C7 in component iks_np (dimensionless). C STATES(87) is C8 in component iks_np (dimensionless). C STATES(88) is C9 in component iks_np (dimensionless). C ALGBRC(14) is a in component iks_np (mS_per_uF). C ALGBRC(42) is b in component iks_np (mS_per_uF). C ALGBRC(64) is d in component iks_np (mS_per_uF). C ALGBRC(76) is e in component iks_np (mS_per_uF). C ALGBRC(92) is g in component iks_np (mS_per_uF). C ALGBRC(93) is o in component iks_np (mS_per_uF). C ALGBRC(104) is p in component iks_np (mS_per_uF). C CONSTS(112) is t in component iks_np (mS_per_uF). C STATES(89) is C1 in component iks_pka (dimensionless). C STATES(90) is C10 in component iks_pka (dimensionless). C STATES(91) is C11 in component iks_pka (dimensionless). C STATES(92) is C12 in component iks_pka (dimensionless). C STATES(93) is C13 in component iks_pka (dimensionless). C STATES(94) is C14 in component iks_pka (dimensionless). C STATES(95) is C15 in component iks_pka (dimensionless). C STATES(96) is C2 in component iks_pka (dimensionless). C STATES(97) is C3 in component iks_pka (dimensionless). C STATES(98) is C4 in component iks_pka (dimensionless). C STATES(99) is C5 in component iks_pka (dimensionless). C STATES(100) is C6 in component iks_pka (dimensionless). C STATES(101) is C7 in component iks_pka (dimensionless). C STATES(102) is C8 in component iks_pka (dimensionless). C STATES(103) is C9 in component iks_pka (dimensionless). C ALGBRC(15) is a in component iks_pka (mS_per_uF). C ALGBRC(43) is b in component iks_pka (mS_per_uF). C ALGBRC(65) is d in component iks_pka (mS_per_uF). C ALGBRC(77) is e in component iks_pka (mS_per_uF). C ALGBRC(94) is g in component iks_pka (mS_per_uF). C ALGBRC(95) is o in component iks_pka (mS_per_uF). C ALGBRC(105) is p in component iks_pka (mS_per_uF). C CONSTS(113) is t in component iks_pka (mS_per_uF). C CONSTS(406) is IKs_arp in component iks_sig (uM). C CONSTS(357) is IKsf in component iks_sig (uM). C STATES(104) is IKsp in component iks_sig (uM). C CONSTS(114) is K in component iks_sig (uM). C CONSTS(115) is Ka_iks in component iks_sig (uM). C CONSTS(116) is Kp_iks in component iks_sig (uM). C CONSTS(117) is L in component iks_sig (uM). C CONSTS(118) is M in component iks_sig (uM). C CONSTS(400) is PKA_eca in component pka (uM). C CONSTS(402) is PKAf in component iks_sig (uM). C CONSTS(361) is PP1f_eca in component iks_sig (uM). C CONSTS(119) is Yotiao in component iks_sig (uM). C CONSTS(403) is Yotiaof in component iks_sig (uM). C CONSTS(354) is iks_sig_IKsf_sum in component iks_sig (dimensionless). C ALGBRC(16) is iks_sig_IKsp_dif in component iks_sig (uM). C CONSTS(401) is iks_sig_PKAf_sum in component iks_sig (dimensionless). C CONSTS(359) is iks_sig_PP1f_eca_sum in component iks_sig (dimensionless). C CONSTS(120) is ka_iks in component iks_sig (hertz). C CONSTS(121) is kp_iks in component iks_sig (hertz). C ALGBRC(344) is IKur in component ikur (uA_per_cm2). C ALGBRC(342) is IKur_np in component ikur (uA_per_cm2). C ALGBRC(343) is IKur_p in component ikur (uA_per_cm2). C CONSTS(122) is Ka_ikur in component ikur (uM). C CONSTS(123) is Kp_ikur in component ikur (uM). C STATES(105) is f_ikur in component ikur (dimensionless). C ALGBRC(271) is fhat in component ikur (dimensionless). C CONSTS(124) is gbar_np in component ikur (mS_per_cm2). C ALGBRC(269) is ikur_fhat_val in component ikur (dimensionless). C CONSTS(125) is ka_ikur in component ikur (hertz). C CONSTS(126) is kp_ikur in component ikur (hertz). C ALGBRC(360) is INa in component ina (uA_per_cm2). C ALGBRC(356) is INa_both in component ina (uA_per_cm2). C ALGBRC(357) is INa_camk in component ina (uA_per_cm2). C ALGBRC(358) is INa_np in component ina (uA_per_cm2). C ALGBRC(359) is INa_pka in component ina (uA_per_cm2). C CONSTS(127) is Ka_ina in component ina (uM). C CONSTS(128) is Kp_ina in component ina (uM). C ALGBRC(276) is f_both in component ina (dimensionless). C ALGBRC(277) is f_camk_only in component ina (dimensionless). C STATES(106) is f_ina in component ina (dimensionless). C ALGBRC(279) is f_np in component ina (dimensionless). C ALGBRC(275) is f_pka in component ina (dimensionless). C ALGBRC(278) is f_pka_only in component ina (dimensionless). C CONSTS(242) is gNaBar in component ina (mS_per_cm2). C STATES(107) is h in component ina_camk (dimensionless). C STATES(108) is j in component ina_camk (dimensionless). C STATES(109) is m in component ina_camk (dimensionless). C ALGBRC(274) is ina_f_pka_val in component ina (dimensionless). C STATES(110) is h in component ina_np (dimensionless). C STATES(111) is j in component ina_np (dimensionless). C STATES(112) is m in component ina_np (dimensionless). C STATES(113) is h in component ina_pka (dimensionless). C STATES(114) is j in component ina_pka (dimensionless). C STATES(115) is m in component ina_pka (dimensionless). C CONSTS(129) is ka_ina in component ina (hertz). C CONSTS(130) is kp_ina in component ina (hertz). C CONSTS(131) is dVIn in component ina_camk (mV). C ALGBRC(17) is ina_camk_h_alpha in component ina_camk (mS_per_uF). C ALGBRC(44) is ina_camk_h_beta in component ina_camk (mS_per_uF). C ALGBRC(18) is ina_camk_j_alpha in component ina_camk (mS_per_uF). C ALGBRC(45) is ina_camk_j_beta in component ina_camk (mS_per_uF). C ALGBRC(19) is ina_camk_m_alpha in component ina_camk (mS_per_uF). C ALGBRC(46) is ina_camk_m_beta in component ina_camk (mS_per_uF). C ALGBRC(20) is ina_np_h_alpha in component ina_np (mS_per_uF). C ALGBRC(47) is ina_np_h_beta in component ina_np (mS_per_uF). C ALGBRC(21) is ina_np_j_alpha in component ina_np (mS_per_uF). C ALGBRC(48) is ina_np_j_beta in component ina_np (mS_per_uF). C ALGBRC(22) is ina_np_m_alpha in component ina_np (mS_per_uF). C ALGBRC(49) is ina_np_m_beta in component ina_np (mS_per_uF). C CONSTS(132) is dVAc in component ina_pka (mV). C CONSTS(133) is dVIn in component ina_pka (mV). C ALGBRC(23) is ina_pka_h_alpha in component ina_pka (mS_per_uF). C ALGBRC(50) is ina_pka_h_beta in component ina_pka (mS_per_uF). C ALGBRC(24) is ina_pka_j_alpha in component ina_pka (mS_per_uF). C ALGBRC(51) is ina_pka_j_beta in component ina_pka (mS_per_uF). C ALGBRC(25) is ina_pka_m_alpha in component ina_pka (mS_per_uF). C ALGBRC(52) is ina_pka_m_beta in component ina_pka (mS_per_uF). C ALGBRC(282) is INab in component inab (uA_per_cm2). C CONSTS(134) is P in component inab (cm_per_s). C ALGBRC(281) is ePhi in component inab (dimensionless). C ALGBRC(280) is inab_INab_phi in component inab (dimensionless). C CONSTS(363) is KmNai3 in component inaca (mM3). C CONSTS(365) is KmNao3 in component inaca (mM3). C CONSTS(135) is Km_Ca in component inaca (mM). C CONSTS(136) is Km_Cai in component inaca (mM). C CONSTS(137) is Km_Cao in component inaca (mM). C CONSTS(138) is Km_Nai in component inaca (mM). C CONSTS(139) is Km_Nao in component inaca (mM). C ALGBRC(283) is Na_i3 in component inaca (mM3). C CONSTS(367) is Na_o3 in component inaca (mM3). C ALGBRC(284) is Na_ss3 in component inaca (mM3). C CONSTS(140) is eta in component inaca (dimensionless). C ALGBRC(285) is exp1 in component inaca (dimensionless). C ALGBRC(286) is exp2 in component inaca (dimensionless). C ALGBRC(287) is inaca_INaCaSR_denom1 in component inaca (dimensionless). C ALGBRC(288) is inaca_INaCaSR_denom2 in component inaca (dimensionless). C ALGBRC(289) is inaca_INaCaSR_denom3 in component inaca (mol4_per_m12). C ALGBRC(290) is inaca_INaCaSR_denom4 in component inaca (mol4_per_m12). C ALGBRC(291) is inaca_INaCaSR_num in component inaca (A_mol4_per_m14_times_1e_minus_2). C ALGBRC(293) is inaca_INaCa_denom1 in component inaca (dimensionless). C ALGBRC(294) is inaca_INaCa_denom2 in component inaca (dimensionless). C ALGBRC(295) is inaca_INaCa_denom3 in component inaca (mol4_per_m12). C ALGBRC(296) is inaca_INaCa_denom4 in component inaca (mol4_per_m12). C ALGBRC(297) is inaca_INaCa_num in component inaca (A_mol4_per_m14_times_1e_minus_2). C CONSTS(141) is kSat in component inaca (dimensionless). C CONSTS(142) is vMax in component inaca (uA_per_cm2). C ALGBRC(304) is INaK in component inak (uA_per_cm2). C ALGBRC(302) is INaK_np in component inak (uA_per_cm2). C ALGBRC(303) is INaK_p in component inak (uA_per_cm2). C CONSTS(143) is Ka_inak in component inak (uM). C CONSTS(144) is Kp_inak in component inak (uM). C STATES(116) is f_inak in component inak (dimensionless). C ALGBRC(300) is fhat in component inak (dimensionless). C CONSTS(145) is ibar in component inak (uA_per_cm2). C ALGBRC(299) is inak_fhat_val in component inak (dimensionless). C CONSTS(146) is ka_inak in component inak (hertz). C CONSTS(147) is km_ko in component inak (mM). C CONSTS(148) is km_np in component inak (mM). C CONSTS(149) is km_p in component inak (mM). C CONSTS(150) is kp_inak in component inak (hertz). C ALGBRC(301) is phi in component inak (uA_per_cm2). C CONSTS(243) is pk in component inak (dimensionless). C ALGBRC(364) is INaL in component inal (uA_per_cm2). C ALGBRC(362) is INaL_camk in component inal (uA_per_cm2). C ALGBRC(363) is INaL_np in component inal (uA_per_cm2). C ALGBRC(361) is conductance in component inal (mV). C STATES(117) is h in component inal (dimensionless). C ALGBRC(26) is h_inf in component inal (dimensionless). C ALGBRC(27) is inal_m_alpha in component inal (mS_per_uF). C ALGBRC(53) is inal_m_beta in component inal (mS_per_uF). C STATES(118) is m in component inal (dimensionless). C CONSTS(151) is tau_h in component inal (ms). C CONSTS(152) is IpCa_bar in component ipca (uA_per_cm2). C CONSTS(153) is Km_pCa in component ipca (mM). C ALGBRC(316) is Ileak_ryr in component irel (mM_per_ms). C ALGBRC(308) is Ileak_ryr_np in component irel (mM_per_ms). C ALGBRC(310) is Ileak_ryr_p in component irel (mM_per_ms). C STATES(119) is Irel_np in component irel (mM_per_ms). C STATES(120) is Irel_p in component irel (mM_per_ms). C CONSTS(154) is Km_ryr_leak_np in component irel (mM). C CONSTS(155) is Km_ryr_leak_p in component irel (mM). C ALGBRC(256) is alpha_np in component irel (mM_per_ms). C ALGBRC(262) is alpha_p in component irel (mM_per_ms). C CONSTS(244) is beta_0 in component irel (ms). C ALGBRC(253) is beta_np in component irel (ms). C ALGBRC(259) is beta_p in component irel (ms). C ALGBRC(314) is fhat in component irel (dimensionless). C CONSTS(399) is irel_fhat_ratio in component irel (dimensionless). C ALGBRC(312) is irel_fhat_val in component irel (dimensionless). C ALGBRC(272) is irel_inf_np in component irel (mM_per_ms). C ALGBRC(273) is irel_inf_p in component irel (mM_per_ms). C ALGBRC(266) is irel_tau_np in component irel (ms). C ALGBRC(268) is irel_tau_p in component irel (ms). C CONSTS(156) is k_ryr_leak_np in component irel (mS_per_uF). C CONSTS(157) is k_ryr_leak_p in component irel (mS_per_uF). C ALGBRC(264) is x in component irel (dimensionless). C ALGBRC(270) is y in component irel (dimensionless). C CONSTS(158) is Gbar in component ito (mS_per_cm2). C ALGBRC(348) is ITo in component ito (uA_per_cm2). C ALGBRC(346) is ITo_camk in component ito (uA_per_cm2). C ALGBRC(347) is ITo_np in component ito (uA_per_cm2). C ALGBRC(326) is R in component ito (dimensionless). C ALGBRC(28) is a_inf in component ito (dimensionless). C STATES(121) is a_np in component ito (dimensionless). C ALGBRC(78) is a_tau in component ito (ms). C ALGBRC(29) is alph_if in component ito (dimensionless). C ALGBRC(55) is alph_is in component ito (dimensionless). C ALGBRC(67) is beta_i in component ito (mS_per_uF). C STATES(122) is if_camk in component ito (dimensionless). C STATES(123) is if_np in component ito (dimensionless). C STATES(124) is is_camk in component ito (dimensionless). C STATES(125) is is_np in component ito (dimensionless). C ALGBRC(54) is ito_a_np_alpha in component ito (dimensionless). C ALGBRC(66) is ito_a_np_beta in component ito (dimensionless). C ALGBRC(79) is ito_if_camk_alpha in component ito (mS_per_uF). C ALGBRC(80) is ito_if_np_alpha in component ito (mS_per_uF). C ALGBRC(81) is ito_is_camk_alpha in component ito (mS_per_uF). C ALGBRC(82) is ito_is_np_alpha in component ito (mS_per_uF). C ALGBRC(345) is x in component ito (uA_per_cm2). C ALGBRC(309) is Imax in component iup (mM_per_ms). C CONSTS(159) is Ka_plb in component iup (uM). C CONSTS(377) is Km_both in component iup (mM). C CONSTS(370) is Km_camk in component iup (mM). C CONSTS(160) is Km_np in component iup (mM). C CONSTS(375) is Km_pka in component iup (mM). C ALGBRC(323) is Km_up in component iup (mM). C CONSTS(161) is Kp_plb in component iup (uM). C ALGBRC(307) is f_SERCA2a in component iup (dimensionless). C ALGBRC(315) is f_both in component iup (dimensionless). C ALGBRC(317) is f_camk_only in component iup (dimensionless). C ALGBRC(322) is f_np in component iup (dimensionless). C ALGBRC(313) is f_pka in component iup (dimensionless). C ALGBRC(319) is f_pka_only in component iup (dimensionless). C STATES(126) is f_plb in component iup (dimensionless). C ALGBRC(311) is iup_f_pka_val in component iup (dimensionless). C CONSTS(162) is iupmax in component iup (mM_per_ms). C CONSTS(380) is iupmaxCAMK in component iup (mM_per_ms). C CONSTS(163) is ka_plb in component iup (hertz). C CONSTS(164) is kp_plb in component iup (hertz). C ALGBRC(325) is leak in component iup (mM_per_ms). C CONSTS(165) is nsrmax in component iup (mM). C ALGBRC(327) is uptake in component iup (mM_per_ms). C ALGBRC(365) is IK_tot in component potassium (uA_per_cm2). C ALGBRC(367) is INa_tot in component sodium (uA_per_cm2). C ALGBRC(368) is i_ion in component membrane (uA_per_cm2). C ALGBRC(369) is i_stim in component stimulus (uA_per_cm2). C CONSTS(166) is PNaK in component nernst (dimensionless). C CONSTS(239) is RTF in component phys (mV). C STATES(127) is K in component potassium (mM). C CONSTS(167) is KPDEp in component pde (uM). C CONSTS(168) is KmIbmxPde2 in component pde (dimensionless). C CONSTS(169) is KmIbmxPde3 in component pde (dimensionless). C CONSTS(170) is KmIbmxPde4 in component pde (dimensionless). C CONSTS(171) is KmPDE2 in component pde (uM). C CONSTS(172) is KmPDE3 in component pde (uM). C CONSTS(173) is KmPDE4 in component pde (uM). C CONSTS(348) is PDE2_cav in component pde (uM). C CONSTS(349) is PDE2_cyt in component pde (uM). C CONSTS(350) is PDE2_eca in component pde (uM). C CONSTS(174) is PDE2_tot in component pde (uM). C STATES(128) is PDE3_P_cav in component pde (uM). C STATES(129) is PDE3_P_cyt in component pde (uM). C CONSTS(378) is PDE3_cav in component pde (uM). C CONSTS(379) is PDE3_cyt in component pde (uM). C CONSTS(364) is PDE3_tot in component pde (uM). C STATES(130) is PDE4_P_cav in component pde (uM). C STATES(131) is PDE4_P_cyt in component pde (uM). C STATES(132) is PDE4_P_eca in component pde (uM). C CONSTS(371) is PDE4_cav in component pde (uM). C CONSTS(372) is PDE4_cyt in component pde (uM). C CONSTS(373) is PDE4_eca in component pde (uM). C CONSTS(366) is PDE4_tot in component pde (uM). C CONSTS(175) is delta_k_pde34 in component pde (dimensionless). C CONSTS(176) is f_pde2_cav in component pde (dimensionless). C CONSTS(334) is f_pde2_cyt in component pde (dimensionless). C CONSTS(177) is f_pde2_eca in component pde (dimensionless). C CONSTS(335) is f_pde2_part in component pde (dimensionless). C CONSTS(368) is f_pde3_cav in component pde (dimensionless). C CONSTS(374) is f_pde3_cyt in component pde (dimensionless). C CONSTS(178) is f_pde4_cav in component pde (dimensionless). C CONSTS(340) is f_pde4_cyt in component pde (dimensionless). C CONSTS(344) is f_pde4_eca in component pde (dimensionless). C CONSTS(179) is f_pde4_part in component pde (dimensionless). C CONSTS(180) is f_pde_part in component pde (dimensionless). C CONSTS(181) is ff_pde3_cyt in component pde (dimensionless). C CONSTS(182) is h_ibmx_pde2 in component pde (dimensionless). C CONSTS(183) is h_ibmx_pde3 in component pde (dimensionless). C CONSTS(184) is h_ibmx_pde4 in component pde (dimensionless). C CONSTS(185) is ibmx in component pde (uM). C CONSTS(343) is ibmx2 in component pde (uM). C CONSTS(376) is ibmx3 in component pde (uM). C CONSTS(369) is ibmx4 in component pde (uM). C CONSTS(339) is ibmx_h2 in component pde (dimensionless). C CONSTS(351) is ibmx_h3 in component pde (dimensionless). C CONSTS(355) is ibmx_h4 in component pde (dimensionless). C CONSTS(186) is kPDE2 in component pde (hertz). C CONSTS(187) is kPDE3 in component pde (hertz). C CONSTS(188) is kPDE4 in component pde (hertz). C CONSTS(245) is kbPDEp in component pde (hertz). C CONSTS(189) is kfPDEp in component pde (per_mM_per_ms). C CONSTS(360) is pde_PDE3_tot_alpha in component pde (dimensionless). C CONSTS(362) is pde_PDE3_tot_beta in component pde (dimensionless). C CONSTS(190) is r_pde34_frac in component pde (dimensionless). C CONSTS(358) is r_pde3_cyt in component pde (dimensionless). C CONSTS(191) is R in component phys (mJ_per_mol_per_K). C CONSTS(192) is T in component phys (kelvin). C CONSTS(246) is K_pki in component pka (uM). C CONSTS(336) is PKA_cyt in component pka (uM). C CONSTS(193) is PKA_tot in component pka (uM). C CONSTS(345) is PKI_cav in component pka (uM). C CONSTS(356) is PKI_cyt in component pka (uM). C CONSTS(352) is PKI_eca in component pka (uM). C CONSTS(337) is PKI_tot in component pka (uM). C CONSTS(265) is b_pki in component pka (hertz). C CONSTS(194) is f_cav in component pka (dimensionless). C CONSTS(328) is f_cyt in component pka (dimensionless). C CONSTS(195) is f_eca in component pka (dimensionless). C CONSTS(196) is f_pki in component pka (per_mM_per_ms). C CONSTS(341) is f_pki_cav in component pka (dimensionless). C CONSTS(353) is f_pki_cyt in component pka (dimensionless). C CONSTS(346) is f_pki_eca in component pka (dimensionless). C STATES(133) is A2R in component pka_cav (uM). C STATES(134) is A2RC in component pka_cav (uM). C STATES(135) is ARC in component pka_cav (uM). C CONSTS(197) is K1 in component pka_cav (uM). C CONSTS(198) is K2 in component pka_cav (uM). C CONSTS(199) is K3 in component pka_cav (m3_per_mol_times_1e3). C STATES(136) is PKIC in component pka_cav (uM). C ALGBRC(141) is RCf in component pka_cav (uM). C CONSTS(247) is b1 in component pka_cav (hertz). C CONSTS(248) is b2 in component pka_cav (hertz). C CONSTS(249) is b3 in component pka_cav (per_mM_per_ms). C CONSTS(200) is f1 in component pka_cav (per_mM_per_ms). C CONSTS(201) is f2 in component pka_cav (per_mM_per_ms). C CONSTS(202) is f3 in component pka_cav (hertz). C STATES(137) is A2R in component pka_cyt (uM). C STATES(138) is A2RC in component pka_cyt (uM). C STATES(139) is ARC in component pka_cyt (uM). C CONSTS(203) is K1 in component pka_cyt (uM). C CONSTS(204) is K2 in component pka_cyt (uM). C CONSTS(205) is K3 in component pka_cyt (m3_per_mol_times_1e3). C STATES(140) is PKIC in component pka_cyt (uM). C ALGBRC(151) is RCf in component pka_cyt (uM). C CONSTS(266) is b1 in component pka_cyt (hertz). C CONSTS(267) is b2 in component pka_cyt (hertz). C CONSTS(268) is b3 in component pka_cyt (per_mM_per_ms). C CONSTS(250) is f1 in component pka_cyt (per_mM_per_ms). C CONSTS(251) is f2 in component pka_cyt (per_mM_per_ms). C CONSTS(252) is f3 in component pka_cyt (hertz). C STATES(141) is A2R in component pka_eca (uM). C STATES(142) is A2RC in component pka_eca (uM). C STATES(143) is ARC in component pka_eca (uM). C CONSTS(253) is K1 in component pka_eca (uM). C CONSTS(254) is K2 in component pka_eca (uM). C CONSTS(255) is K3 in component pka_eca (m3_per_mol_times_1e3). C STATES(144) is PKIC in component pka_eca (uM). C ALGBRC(131) is RCf in component pka_eca (uM). C CONSTS(277) is b1 in component pka_eca (hertz). C CONSTS(278) is b2 in component pka_eca (hertz). C CONSTS(279) is b3 in component pka_eca (per_mM_per_ms). C CONSTS(269) is f1 in component pka_eca (per_mM_per_ms). C CONSTS(270) is f2 in component pka_eca (per_mM_per_ms). C CONSTS(271) is f3 in component pka_eca (hertz). C CONSTS(407) is potassium_K_r1 in component potassium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(206) is K in component pp1 (uM). C CONSTS(207) is Kdp in component pp1 (uM). C CONSTS(208) is Kp in component pp1 (uM). C CONSTS(209) is PP1_cyt in component pp1 (uM). C ALGBRC(30) is di in component pp1 (uM). C CONSTS(210) is f in component pp1 (dimensionless). C STATES(145) is inhib1_p in component pp1 (uM). C CONSTS(256) is inhib1_tot in component pp1 (uM). C CONSTS(211) is kdp in component pp1 (hertz). C CONSTS(212) is kp in component pp1 (hertz). C ALGBRC(199) is pp1_PP1f_cyt_sum in component pp1 (uM). C ALGBRC(366) is INa_cyt in component sodium (uA_per_cm2). C CONSTS(408) is sodium_Na_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(409) is sodium_Na_r2 in component sodium (dimensionless). C CONSTS(410) is sodium_Na_sr_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5). C CONSTS(257) is amplitude in component stimulus (uA_per_cm2). C CONSTS(213) is duration in component stimulus (ms). C CONSTS(214) is offset in component stimulus (ms). C CONSTS(215) is period in component stimulus (ms). C RATES(5) is d/dt ICaLp in component akap_sig (uM). C RATES(6) is d/dt RyRp in component akap_sig (uM). C RATES(8) is d/dt Gi_aGDP in component beta_cav (uM). C RATES(9) is d/dt Gi_aGTP in component beta_cav (uM). C RATES(1) is d/dt Gi_bg in component beta_cav (uM). C RATES(10) is d/dt Gs_aGDP in component beta_cav (uM). C RATES(2) is d/dt Gs_aGTP in component beta_cav (uM). C RATES(11) is d/dt Gs_bg in component beta_cav (uM). C RATES(12) is d/dt Rb1_grk_tot in component beta_cav (uM). C RATES(13) is d/dt Rb1_pka_tot in component beta_cav (uM). C RATES(14) is d/dt Rb2_grk_tot in component beta_cav (uM). C RATES(15) is d/dt Rb2_pka_tot in component beta_cav (uM). C RATES(16) is d/dt Gs_aGDP in component beta_cyt (uM). C RATES(3) is d/dt Gs_aGTP in component beta_cyt (uM). C RATES(17) is d/dt Gs_bg in component beta_cyt (uM). C RATES(18) is d/dt Rb1_grk_tot in component beta_cyt (uM). C RATES(19) is d/dt Rb1_pka_tot in component beta_cyt (uM). C RATES(21) is d/dt Gi_aGDP in component beta_eca (uM). C RATES(22) is d/dt Gi_aGTP in component beta_eca (uM). C RATES(23) is d/dt Gi_bg in component beta_eca (uM). C RATES(24) is d/dt Gs_aGDP in component beta_eca (uM). C RATES(4) is d/dt Gs_aGTP in component beta_eca (uM). C RATES(25) is d/dt Gs_bg in component beta_eca (uM). C RATES(26) is d/dt Rb1_grk_tot in component beta_eca (uM). C RATES(27) is d/dt Rb1_pka_tot in component beta_eca (uM). C RATES(28) is d/dt Rb2_grk_tot in component beta_eca (uM). C RATES(29) is d/dt Rb2_pka_tot in component beta_eca (uM). C RATES(31) is d/dt Ca_nsr in component calcium (mM). C RATES(32) is d/dt f_tni in component calcium (dimensionless). C RATES(33) is d/dt uCa in component calcium (mM). C RATES(34) is d/dt uCa_CaL in component calcium (mM). C RATES(35) is d/dt uCa_jsr in component calcium (mM). C RATES(36) is d/dt uCa_sr in component calcium (mM). C RATES(37) is d/dt f_ical in component camk (dimensionless). C RATES(38) is d/dt f_ik1 in component camk (dimensionless). C RATES(39) is d/dt f_ina in component camk (dimensionless). C RATES(40) is d/dt f_ito in component camk (dimensionless). C RATES(41) is d/dt f_plb in component camk (dimensionless). C RATES(42) is d/dt f_ryr in component camk (dimensionless). C RATES(43) is d/dt trap in component camk (dimensionless). C RATES(44) is d/dt cAMP_cav in component camp (uM). C RATES(45) is d/dt cAMP_cyt in component camp (uM). C RATES(46) is d/dt cAMP_eca in component camp (uM). C RATES(47) is d/dt Cl in component chloride (mM). C RATES(48) is d/dt Cl_sr in component chloride (mM). C RATES(52) is d/dt C in component ical_camk (dimensionless). C RATES(53) is d/dt CI in component ical_camk (dimensionless). C RATES(54) is d/dt CIs in component ical_camk (dimensionless). C RATES(55) is d/dt Cs in component ical_camk (dimensionless). C RATES(56) is d/dt O in component ical_camk (dimensionless). C RATES(57) is d/dt OI in component ical_camk (dimensionless). C RATES(58) is d/dt OIs in component ical_camk (dimensionless). C RATES(59) is d/dt Os in component ical_camk (dimensionless). C RATES(60) is d/dt C in component ical_np (dimensionless). C RATES(61) is d/dt CI in component ical_np (dimensionless). C RATES(62) is d/dt CIs in component ical_np (dimensionless). C RATES(63) is d/dt Cs in component ical_np (dimensionless). C RATES(64) is d/dt O in component ical_np (dimensionless). C RATES(65) is d/dt OI in component ical_np (dimensionless). C RATES(66) is d/dt OIs in component ical_np (dimensionless). C RATES(67) is d/dt Os in component ical_np (dimensionless). C RATES(68) is d/dt i2 in component iclca (dimensionless). C RATES(69) is d/dt ac in component ikr (dimensionless). C RATES(74) is d/dt C1 in component iks_np (dimensionless). C RATES(75) is d/dt C10 in component iks_np (dimensionless). C RATES(76) is d/dt C11 in component iks_np (dimensionless). C RATES(77) is d/dt C12 in component iks_np (dimensionless). C RATES(78) is d/dt C13 in component iks_np (dimensionless). C RATES(79) is d/dt C14 in component iks_np (dimensionless). C RATES(80) is d/dt C15 in component iks_np (dimensionless). C RATES(81) is d/dt C2 in component iks_np (dimensionless). C RATES(82) is d/dt C3 in component iks_np (dimensionless). C RATES(83) is d/dt C4 in component iks_np (dimensionless). C RATES(84) is d/dt C5 in component iks_np (dimensionless). C RATES(85) is d/dt C6 in component iks_np (dimensionless). C RATES(86) is d/dt C7 in component iks_np (dimensionless). C RATES(87) is d/dt C8 in component iks_np (dimensionless). C RATES(88) is d/dt C9 in component iks_np (dimensionless). C RATES(70) is d/dt O1 in component iks_np (dimensionless). C RATES(71) is d/dt O2 in component iks_np (dimensionless). C RATES(89) is d/dt C1 in component iks_pka (dimensionless). C RATES(90) is d/dt C10 in component iks_pka (dimensionless). C RATES(91) is d/dt C11 in component iks_pka (dimensionless). C RATES(92) is d/dt C12 in component iks_pka (dimensionless). C RATES(93) is d/dt C13 in component iks_pka (dimensionless). C RATES(94) is d/dt C14 in component iks_pka (dimensionless). C RATES(95) is d/dt C15 in component iks_pka (dimensionless). C RATES(96) is d/dt C2 in component iks_pka (dimensionless). C RATES(97) is d/dt C3 in component iks_pka (dimensionless). C RATES(98) is d/dt C4 in component iks_pka (dimensionless). C RATES(99) is d/dt C5 in component iks_pka (dimensionless). C RATES(100) is d/dt C6 in component iks_pka (dimensionless). C RATES(101) is d/dt C7 in component iks_pka (dimensionless). C RATES(102) is d/dt C8 in component iks_pka (dimensionless). C RATES(103) is d/dt C9 in component iks_pka (dimensionless). C RATES(72) is d/dt O1 in component iks_pka (dimensionless). C RATES(73) is d/dt O2 in component iks_pka (dimensionless). C RATES(104) is d/dt IKsp in component iks_sig (uM). C RATES(105) is d/dt f_ikur in component ikur (dimensionless). C RATES(106) is d/dt f_ina in component ina (dimensionless). C RATES(107) is d/dt h in component ina_camk (dimensionless). C RATES(108) is d/dt j in component ina_camk (dimensionless). C RATES(109) is d/dt m in component ina_camk (dimensionless). C RATES(110) is d/dt h in component ina_np (dimensionless). C RATES(111) is d/dt j in component ina_np (dimensionless). C RATES(112) is d/dt m in component ina_np (dimensionless). C RATES(113) is d/dt h in component ina_pka (dimensionless). C RATES(114) is d/dt j in component ina_pka (dimensionless). C RATES(115) is d/dt m in component ina_pka (dimensionless). C RATES(116) is d/dt f_inak in component inak (dimensionless). C RATES(117) is d/dt h in component inal (dimensionless). C RATES(118) is d/dt m in component inal (dimensionless). C RATES(119) is d/dt Irel_np in component irel (mM_per_ms). C RATES(120) is d/dt Irel_p in component irel (mM_per_ms). C RATES(121) is d/dt a_np in component ito (dimensionless). C RATES(122) is d/dt if_camk in component ito (dimensionless). C RATES(123) is d/dt if_np in component ito (dimensionless). C RATES(124) is d/dt is_camk in component ito (dimensionless). C RATES(125) is d/dt is_np in component ito (dimensionless). C RATES(126) is d/dt f_plb in component iup (dimensionless). C RATES(51) is d/dt V in component membrane (mV). C RATES(128) is d/dt PDE3_P_cav in component pde (uM). C RATES(129) is d/dt PDE3_P_cyt in component pde (uM). C RATES(130) is d/dt PDE4_P_cav in component pde (uM). C RATES(131) is d/dt PDE4_P_cyt in component pde (uM). C RATES(132) is d/dt PDE4_P_eca in component pde (uM). C RATES(133) is d/dt A2R in component pka_cav (uM). C RATES(134) is d/dt A2RC in component pka_cav (uM). C RATES(135) is d/dt ARC in component pka_cav (uM). C RATES(7) is d/dt C in component pka_cav (uM). C RATES(136) is d/dt PKIC in component pka_cav (uM). C RATES(137) is d/dt A2R in component pka_cyt (uM). C RATES(138) is d/dt A2RC in component pka_cyt (uM). C RATES(139) is d/dt ARC in component pka_cyt (uM). C RATES(20) is d/dt C in component pka_cyt (uM). C RATES(140) is d/dt PKIC in component pka_cyt (uM). C RATES(141) is d/dt A2R in component pka_eca (uM). C RATES(142) is d/dt A2RC in component pka_eca (uM). C RATES(143) is d/dt ARC in component pka_eca (uM). C RATES(30) is d/dt C in component pka_eca (uM). C RATES(144) is d/dt PKIC in component pka_eca (uM). C RATES(127) is d/dt K in component potassium (mM). C RATES(145) is d/dt inhib1_p in component pp1 (uM). C RATES(49) is d/dt Na in component sodium (mM). C RATES(50) is d/dt Na_sr in component sodium (mM). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 5000.0 CONSTS(2) = 315.0 CONSTS(3) = 0.0465 CONSTS(4) = 0.031544 CONSTS(5) = 0.0852 CONSTS(6) = 0.4824 CONSTS(7) = 3.3757 CONSTS(8) = 41.32 CONSTS(9) = 0.03135 CONSTS(10) = 0.037696 STATES(1) = 2.09911481235842013e-03 STATES(2) = 6.85041638458664965e-03 STATES(3) = 7.31420577213055985e-04 STATES(4) = 1.84627603007976003e-02 CONSTS(11) = 0.16479 CONSTS(12) = 0.087459 CONSTS(13) = 1.0043 CONSTS(14) = 1.3574 CONSTS(15) = 0.6623 CONSTS(16) = 0.8569 CONSTS(17) = 0.025 CONSTS(18) = 0.025 STATES(5) = 6.73713947839316954e-04 CONSTS(19) = 1.27019999999999993e-06 CONSTS(20) = 6.62979999999999944e-05 CONSTS(21) = 0.01 CONSTS(22) = 0.0063064 CONSTS(23) = 0.043003 CONSTS(24) = 0.01 CONSTS(25) = 0.0001 CONSTS(26) = 0.0001 CONSTS(27) = 0.01 CONSTS(28) = 0.01 CONSTS(29) = 0.25 CONSTS(30) = 0.125 CONSTS(31) = 0.125 STATES(6) = 4.10693810508170991e-03 CONSTS(32) = 5.10090000000000044e-04 CONSTS(33) = 0.0025548 CONSTS(34) = 0.0006903 CONSTS(35) = 0.0038257 STATES(7) = 3.26565916584702978e-02 CONSTS(36) = 0.5 CONSTS(37) = 0.85 CONSTS(38) = 0.0011071 CONSTS(39) = 0.5664 CONSTS(40) = 0.081161 CONSTS(41) = 0.48744 CONSTS(42) = 0.85 CONSTS(43) = 4.0 CONSTS(44) = 4.9054 CONSTS(45) = 0.05 CONSTS(46) = 0.25945 CONSTS(47) = 2.449 CONSTS(48) = 0.062 CONSTS(49) = 0.567 CONSTS(50) = 1.6655 CONSTS(51) = 1.8463 CONSTS(52) = 0.1 CONSTS(53) = 0.012 CONSTS(54) = 1.053 CONSTS(55) = 1.053 CONSTS(56) = 0.8 CONSTS(57) = 1210.0 CONSTS(58) = 0.35 CONSTS(59) = 1.0 STATES(8) = 5.02792845976641014e-04 STATES(9) = 1.59632196638178009e-03 STATES(10) = 6.07316088556675987e-04 STATES(11) = 7.45773247314215036e-03 STATES(12) = 2.49592854373432000e-10 STATES(13) = 1.49041813757830998e-02 STATES(14) = 8.91799633266019011e-10 STATES(15) = 2.75455839709412009e-02 CONSTS(60) = 0.0 CONSTS(61) = 1.0 CONSTS(62) = 1.0 STATES(16) = 4.19991861054322011e-04 STATES(17) = 1.15141243826746994e-03 STATES(18) = 7.07824478944670999e-11 STATES(19) = 9.44463350378085993e-03 STATES(20) = 3.62113356111495976e-01 CONSTS(63) = 1.0 STATES(21) = 3.41341142614041016e-04 STATES(22) = 3.64315164237569004e-04 STATES(23) = 7.05656306851923029e-04 STATES(24) = 6.39038440072506948e-04 STATES(25) = 1.91017987408719017e-02 STATES(26) = 1.18055788874765002e-09 STATES(27) = 2.03016833596287999e-01 STATES(28) = 1.13428924662652000e-10 STATES(29) = 1.10248953370551007e-02 CONSTS(64) = 1.0 STATES(30) = 5.67249910261072965e-01 STATES(31) = 1.191 CONSTS(65) = 2.71430000000000008e-05 CONSTS(66) = 0.26714 CONSTS(67) = 1.0 CONSTS(68) = 1.124 CONSTS(69) = 0.0087 CONSTS(70) = 0.047 CONSTS(71) = 0.00087 CONSTS(72) = 0.05 CONSTS(73) = 10.0 CONSTS(74) = 0.8 STATES(32) = 6.73518785672381992e-01 CONSTS(75) = 0.10408 CONSTS(76) = 0.00238 CONSTS(77) = 0.052633 CONSTS(78) = 0.0005 CONSTS(79) = 0.07 STATES(33) = 0.013394 STATES(34) = 0.023413 STATES(35) = 6.8659 STATES(36) = 0.023413 CONSTS(80) = 0.05 CONSTS(81) = 0.25 CONSTS(82) = 0.0015 CONSTS(83) = 0.1 CONSTS(84) = 0.05 CONSTS(85) = 0.00068 STATES(37) = 0.0 STATES(38) = 0.0 STATES(39) = 0.0 STATES(40) = 0.0 STATES(41) = 0.0 STATES(42) = 0.0 CONSTS(86) = 100000.0 CONSTS(87) = 10000.0 STATES(43) = 0.0017546 STATES(44) = 3.47102959606005013e-01 STATES(45) = 4.74081735738210996e-01 STATES(46) = 9.62359241535767040e+00 CONSTS(88) = 96487.0 CONSTS(89) = 0.01 CONSTS(90) = 3.14159265358979312e+00 CONSTS(91) = 0.0011 STATES(47) = 20.273 STATES(48) = 20.273 CONSTS(92) = 1.77e-05 CONSTS(93) = 87.8251 CONSTS(94) = 2.46108000000000016e-05 STATES(49) = 6.8909 STATES(50) = 6.8909 CONSTS(95) = 0.2 CONSTS(96) = 0.02 CONSTS(97) = 75.0 CONSTS(98) = 1.8 CONSTS(99) = 100.0 CONSTS(100) = 5.4 CONSTS(101) = 140.0 STATES(51) = -87.491 CONSTS(102) = 1.995e-07 STATES(52) = 1.0 STATES(53) = 0.0 STATES(54) = 0.0 STATES(55) = 0.0 STATES(56) = 0.0 STATES(57) = 0.0 STATES(58) = 0.0 STATES(59) = 0.0 CONSTS(103) = 1.0 CONSTS(104) = 1e-06 STATES(60) = 1.0 STATES(61) = 0.0 STATES(62) = 0.0 STATES(63) = 0.0 STATES(64) = 0.0 STATES(65) = 0.0 STATES(66) = 0.0 STATES(67) = 0.0 CONSTS(105) = 1.0 CONSTS(106) = 1e-06 CONSTS(107) = 0.000225 CONSTS(108) = 9e-07 STATES(68) = 0.99604 CONSTS(109) = 0.4 CONSTS(110) = 8.0 STATES(69) = 1.23059999999999995e-08 CONSTS(111) = 0.025 STATES(70) = 9.37220000000000071e-16 STATES(71) = 1.65950000000000014e-17 STATES(72) = 1.12010000000000000e-15 STATES(73) = 1.66129999999999997e-18 STATES(74) = 0.91141 STATES(75) = 5.36959999999999959e-07 STATES(76) = 2.48610000000000001e-08 STATES(77) = 2.87759999999999997e-10 STATES(78) = 1.12169999999999999e-10 STATES(79) = 2.59670000000000003e-12 STATES(80) = 8.78740000000000076e-15 STATES(81) = 0.084395 STATES(82) = 0.0029306 STATES(83) = 4.52285000000000024e-05 STATES(84) = 2.61750000000000024e-07 STATES(85) = 0.0011424 STATES(86) = 7.93370000000000029e-05 STATES(87) = 1.83659999999999992e-06 STATES(88) = 1.41719999999999995e-08 CONSTS(112) = 0.0027304 STATES(89) = 0.95624 STATES(90) = 3.19559999999999984e-07 STATES(91) = 7.039e-09 STATES(92) = 3.87629999999999971e-11 STATES(93) = 5.02769999999999999e-11 STATES(94) = 5.53740000000000021e-13 STATES(95) = 2.96639999999999984e-15 STATES(96) = 0.042127 STATES(97) = 6.95969999999999965e-04 STATES(98) = 5.11010000000000042e-06 STATES(99) = 1.407e-08 STATES(100) = 9.02690000000000046e-04 STATES(101) = 2.98259999999999987e-05 STATES(102) = 3.285e-07 STATES(103) = 1.206e-09 CONSTS(113) = 0.0046171 STATES(104) = 7.65988420110534033e-04 CONSTS(114) = 0.01 CONSTS(115) = 9.97940000000000003e-05 CONSTS(116) = 1.11470000000000002e-04 CONSTS(117) = 0.0001 CONSTS(118) = 0.01 CONSTS(119) = 0.025 CONSTS(120) = 0.16305 CONSTS(121) = 1.0542 CONSTS(122) = 0.27623 CONSTS(123) = 0.002331 STATES(105) = 5.89379755147717982e-02 CONSTS(124) = 0.00384 CONSTS(125) = 0.069537 CONSTS(126) = 0.317 CONSTS(127) = 0.10988 CONSTS(128) = 7.8605 STATES(106) = 2.39479458960527997e-01 STATES(107) = 0.83805 STATES(108) = 0.99281 STATES(109) = 6.81269999999999988e-04 STATES(110) = 0.0068172 STATES(111) = 0.99709 STATES(112) = 0.90163 STATES(113) = 0.001236 STATES(114) = 0.99123 STATES(115) = 0.79472 CONSTS(129) = 0.01368 CONSTS(130) = 0.052811 CONSTS(131) = 3.25 CONSTS(132) = 3.7 CONSTS(133) = 4.9 CONSTS(134) = 3.2e-09 CONSTS(135) = 0.000125 CONSTS(136) = 0.0036 CONSTS(137) = 1.3 CONSTS(138) = 12.3 CONSTS(139) = 87.5 CONSTS(140) = 0.27 CONSTS(141) = 0.32 CONSTS(142) = 4.5 CONSTS(143) = 0.0011001 CONSTS(144) = 5.7392 STATES(116) = 1.26345311579565994e-01 CONSTS(145) = 1.4 CONSTS(146) = 0.015265 CONSTS(147) = 1.5 CONSTS(148) = 2.6 CONSTS(149) = 1.846 CONSTS(150) = 0.092455 STATES(117) = 0.36003 STATES(118) = 0.0007053 CONSTS(151) = 600.0 CONSTS(152) = 0.0575 CONSTS(153) = 0.0005 STATES(119) = 3.66750000000000000e-09 STATES(120) = 7.30739999999999981e-09 CONSTS(154) = 20.0 CONSTS(155) = 1.1 CONSTS(156) = 0.000175 CONSTS(157) = 0.0005 CONSTS(158) = 0.4975 STATES(121) = 1.76869999999999985e-05 STATES(122) = 1.0 STATES(123) = 0.99798 STATES(124) = 1.0 STATES(125) = 0.98747 CONSTS(159) = 9.88539999999999992e-04 CONSTS(160) = 0.00092 CONSTS(161) = 0.80737 STATES(126) = 5.92167467082830967e-01 CONSTS(162) = 0.004375 CONSTS(163) = 0.11348 CONSTS(164) = 0.48302 CONSTS(165) = 15.0 CONSTS(166) = 0.01833 STATES(127) = 145.62 CONSTS(167) = 0.52218 CONSTS(168) = 21.58 CONSTS(169) = 2.642 CONSTS(170) = 11.89 CONSTS(171) = 50.0 CONSTS(172) = 0.8 CONSTS(173) = 1.4 CONSTS(174) = 0.029268 STATES(128) = 2.36821659448036986e-02 STATES(129) = 1.28402905095187994e-02 STATES(130) = 6.37363047239019025e-03 STATES(131) = 9.17039986149184062e-03 STATES(132) = 4.29171113639321980e-05 CONSTS(175) = 3.0 CONSTS(176) = 0.16957 CONSTS(177) = 2.12570000000000006e-04 CONSTS(178) = 0.12481 CONSTS(179) = 0.125 CONSTS(180) = 0.2 CONSTS(181) = 0.35 CONSTS(182) = 1.167 CONSTS(183) = 0.7629 CONSTS(184) = 0.9024 CONSTS(185) = 0.0 CONSTS(186) = 20.0 CONSTS(187) = 2.5 CONSTS(188) = 4.0 CONSTS(189) = 0.0196 CONSTS(190) = 3.71 CONSTS(191) = 8314.0 CONSTS(192) = 310.0 CONSTS(193) = 0.5 CONSTS(194) = 0.0388 CONSTS(195) = 0.1 CONSTS(196) = 50.0 STATES(133) = 2.25475702283052998e-01 STATES(134) = 2.76490711096605019e-03 STATES(135) = 9.04820284659604013e-02 CONSTS(197) = 2.4984 CONSTS(198) = 11.359 CONSTS(199) = 0.3755 STATES(136) = 1.92819110624504991e-01 CONSTS(200) = 100.0 CONSTS(201) = 100.0 CONSTS(202) = 100.0 STATES(137) = 4.89063888619455989e-01 STATES(138) = 6.64997605558790977e-02 STATES(139) = 6.46928309115710060e-02 CONSTS(203) = 0.1088 CONSTS(204) = 0.4612 CONSTS(205) = 0.3755 STATES(140) = 1.26950532507959013e-01 STATES(141) = 8.17161796756963987e-01 STATES(142) = 1.74057375932567010e-01 STATES(143) = 2.05444874210056000e-01 STATES(144) = 2.49911886495889995e-01 CONSTS(206) = 0.001 CONSTS(207) = 1.95259999999999991e-05 CONSTS(208) = 0.001469 CONSTS(209) = 0.2 CONSTS(210) = 0.3 STATES(145) = 2.82662056977524001e-02 CONSTS(211) = 0.0035731 CONSTS(212) = 0.010145 CONSTS(213) = 0.5 CONSTS(214) = 100.0 CONSTS(215) = 1000.0 CONSTS(216) = 1.00000 - CONSTS(37) CONSTS(217) = (1.00000 - CONSTS(38)) - CONSTS(39) CONSTS(218) = CONSTS(58)*0.000983300 CONSTS(219) = CONSTS(58)*0.00133000 CONSTS(220) = CONSTS(56) CONSTS(221) = CONSTS(58)*0.00650000 CONSTS(222) = CONSTS(57) CONSTS(223) = ( (CONSTS(52)+CONSTS(60))*(CONSTS(55)+CONSTS(60)))/CONSTS(55) CONSTS(224) = ( ( CONSTS(49)*CONSTS(54))*(CONSTS(48)+CONSTS(60)))*(CONSTS(53)+CONSTS(60)) CONSTS(225) = ( (CONSTS(48)+CONSTS(60))*(CONSTS(49)+CONSTS(60)))/CONSTS(49) CONSTS(226) = ( (CONSTS(52)+CONSTS(60))*(CONSTS(55)+CONSTS(60)))/CONSTS(55) CONSTS(227) = ( ( CONSTS(49)*CONSTS(54))*(CONSTS(48)+CONSTS(60)))*(CONSTS(53)+CONSTS(60)) CONSTS(228) = CONSTS(1)/(CONSTS(2)+CONSTS(1)) CONSTS(229) = CONSTS(70)+CONSTS(68) CONSTS(230) = 1.50000*CONSTS(78) CONSTS(231) = CONSTS(71)*CONSTS(69) CONSTS(232) = CONSTS(87) CONSTS(233) = CONSTS(86) CONSTS(234) = CONSTS(86) CONSTS(235) = CONSTS(86) CONSTS(236) = 7.50000e-14*1.00000e+06 CONSTS(237) = 5.00000e-15*1.00000e+06 CONSTS(238) = 9.00000e-15*1.00000e+06 CONSTS(239) = ( CONSTS(191)*CONSTS(192))/CONSTS(88) CONSTS(240) = 0.500000* (CONSTS(100)/5.40000) ** (1.0 / 2) CONSTS(241) = 0.0138542* (CONSTS(100)/5.40000) ** (1.0 / 2) CONSTS(242) = ( 2.15000*8.25000)*1.10000 CONSTS(243) = CONSTS(100)/(CONSTS(100)+CONSTS(147)) CONSTS(244) = 0.666700*4.75000 CONSTS(245) = CONSTS(167)*CONSTS(189) CONSTS(246) = 0.0100000/50.0000 CONSTS(247) = CONSTS(200)*CONSTS(197) CONSTS(248) = CONSTS(201)*CONSTS(198) CONSTS(249) = CONSTS(202)*CONSTS(199) CONSTS(250) = CONSTS(200) CONSTS(251) = CONSTS(201) CONSTS(252) = CONSTS(202) CONSTS(253) = CONSTS(197) CONSTS(254) = CONSTS(198) CONSTS(255) = CONSTS(199) CONSTS(256) = (CONSTS(210)/(1.00000 - CONSTS(210)))*CONSTS(206)+ CONSTS(210)*CONSTS(209) CONSTS(257) = - 80.0000 CONSTS(258) = 0.156290*CONSTS(221) CONSTS(259) = ( ( ( CONSTS(47)*CONSTS(48))*CONSTS(54))*(CONSTS(53)+CONSTS(60)))*(CONSTS(49)+CONSTS(60)) CONSTS(260) = ( ( ( CONSTS(47)*CONSTS(48))*CONSTS(54))*(CONSTS(53)+CONSTS(60)))*(CONSTS(49)+CONSTS(60)) CONSTS(261) = 1.00000/(1.00000+0.350000) CONSTS(262) = CONSTS(71)+CONSTS(69) CONSTS(263) = ( CONSTS(70)*CONSTS(69)+ CONSTS(68)*CONSTS(71))+CONSTS(231) CONSTS(264) = 1.00000/CONSTS(239) CONSTS(265) = CONSTS(196)*CONSTS(246) CONSTS(266) = CONSTS(250)*CONSTS(203) CONSTS(267) = CONSTS(251)*CONSTS(204) CONSTS(268) = CONSTS(252)*CONSTS(205) CONSTS(269) = CONSTS(200) CONSTS(270) = CONSTS(201) CONSTS(271) = CONSTS(202) CONSTS(272) = ( ( ( CONSTS(51)*CONSTS(53))*CONSTS(49))*(CONSTS(48)+CONSTS(60)))*(CONSTS(54)+CONSTS(60)) CONSTS(273) = ( ( ( CONSTS(51)*CONSTS(53))*CONSTS(49))*(CONSTS(48)+CONSTS(60)))*(CONSTS(54)+CONSTS(60)) CONSTS(274) = 1.00000+(CONSTS(17) - CONSTS(18))/CONSTS(25) CONSTS(275) = CONSTS(229)+CONSTS(262) CONSTS(276) = CONSTS(88)*CONSTS(264) CONSTS(277) = CONSTS(269)*CONSTS(253) CONSTS(278) = CONSTS(270)*CONSTS(254) CONSTS(279) = CONSTS(271)*CONSTS(255) CONSTS(280) = ( ( ( ( CONSTS(47)*CONSTS(51))*CONSTS(48))*CONSTS(53))*(CONSTS(49)+CONSTS(60)))*(CONSTS(54)+CONSTS(60)) CONSTS(281) = ( ( ( ( CONSTS(47)*CONSTS(51))*CONSTS(48))*CONSTS(53))*(CONSTS(49)+CONSTS(60)))*(CONSTS(54)+CONSTS(60)) CONSTS(282) = (CONSTS(25)/2.00000)*( (CONSTS(274) ** 2.00000+( 4.00000*CONSTS(18))/CONSTS(25)) ** (1.0 / 2) - CONSTS(274)) CONSTS(283) = (((CONSTS(17)+CONSTS(30))+CONSTS(21))+CONSTS(24)) - CONSTS(29) CONSTS(284) = (( CONSTS(17)*CONSTS(24)+ CONSTS(30)*CONSTS(21))+ CONSTS(21)*CONSTS(24)) - CONSTS(29)*(CONSTS(21)+CONSTS(24)) CONSTS(285) = ( CONSTS(29)*CONSTS(21))*CONSTS(24) CONSTS(286) = ((( (- CONSTS(285)/27.0000)*CONSTS(283) ** 3.00000 - ( ( ( CONSTS(283)*CONSTS(283))*CONSTS(284))*CONSTS(284))/108.000)+( ( CONSTS(283)*CONSTS(284))*CONSTS(285))/6.00000)+CONSTS(284) ** 3.00000/27.0000)+( CONSTS(285)*CONSTS(285))/4.00000 CONSTS(287) = TERNRY(CONSTS(286).LT.0.00000, - CONSTS(286) ** (1.0 / 2), 0.00000) CONSTS(288) = (((TERNRY(CONSTS(286).GT.0.00000, CONSTS(286) ** (1.0 / 2), 0.00000))+CONSTS(285)/2.00000)+( CONSTS(283)*CONSTS(284))/6.00000) - CONSTS(283) ** 3.00000/27.0000 CONSTS(289) = atan(CONSTS(287)/CONSTS(288))/3.00000 CONSTS(290) = CONSTS(288)*CONSTS(288)+ CONSTS(287)*CONSTS(287) ** 1.00000/6.00000 CONSTS(291) = (CONSTS(284)/3.00000 - ( CONSTS(283)*CONSTS(283))/9.00000)/( CONSTS(290)*CONSTS(290)) CONSTS(292) = ( CONSTS(290)*cos(CONSTS(289)))*(1.00000 - CONSTS(291)) - CONSTS(283)/3.00000 CONSTS(293) = 1.00000+(CONSTS(30) - CONSTS(31))/CONSTS(26) CONSTS(294) = (CONSTS(26)/2.00000)*( (CONSTS(293) ** 2.00000+( 4.00000*CONSTS(31))/CONSTS(26)) ** (1.0 / 2) - CONSTS(293)) CONSTS(295) = 0.850000*0.0250000 CONSTS(296) = 224.000*CONSTS(295) CONSTS(297) = 3.00000*CONSTS(295) CONSTS(298) = 0.150000*0.0250000 CONSTS(299) = (1.00000 - CONSTS(40)) - CONSTS(41) CONSTS(300) = 1.00000 - CONSTS(42) CONSTS(301) = ( ( 2.00000*CONSTS(90))*CONSTS(91))*(CONSTS(91)+CONSTS(89)) CONSTS(302) = 2.00000*CONSTS(301) CONSTS(303) = CONSTS(302)/CONSTS(88) CONSTS(304) = ( ( ( 1000.00*CONSTS(90))*CONSTS(91))*CONSTS(91))*CONSTS(89) CONSTS(305) = CONSTS(304)*0.00200000 CONSTS(306) = CONSTS(303)/( 2.00000*CONSTS(305)) CONSTS(307) = 0.0200000*CONSTS(304) CONSTS(308) = CONSTS(304)*0.678000 CONSTS(309) = CONSTS(303)/( 2.00000*CONSTS(308)) CONSTS(310) = 0.0400000*CONSTS(304) CONSTS(311) = CONSTS(304)*0.00480000 CONSTS(312) = CONSTS(304)*0.0552000 CONSTS(313) = CONSTS(311)/CONSTS(312) CONSTS(314) = CONSTS(312)/CONSTS(308) CONSTS(315) = CONSTS(304)*0.0200000 CONSTS(316) = CONSTS(315)/CONSTS(305) CONSTS(317) = - CONSTS(303)/CONSTS(315) CONSTS(318) = - CONSTS(311)/CONSTS(315) CONSTS(319) = CONSTS(315)/CONSTS(308) CONSTS(320) = CONSTS(304)/CONSTS(307) CONSTS(321) = ( ( CONSTS(12)*CONSTS(261))*CONSTS(297))*CONSTS(320) CONSTS(322) = CONSTS(304)/CONSTS(308) CONSTS(323) = ( CONSTS(40)*CONSTS(295))*CONSTS(320) CONSTS(324) = ( CONSTS(42)*CONSTS(298))*CONSTS(320) CONSTS(325) = ( ( (1.00000 - CONSTS(11))*(1.00000 - CONSTS(261)))*CONSTS(297))*CONSTS(322) CONSTS(326) = ( CONSTS(299)*CONSTS(295))*CONSTS(322) CONSTS(327) = CONSTS(304)/CONSTS(310) CONSTS(328) = (1.00000 - CONSTS(194)) - CONSTS(195) CONSTS(329) = ( ( (1.00000 - CONSTS(12))*CONSTS(261))*CONSTS(297))*CONSTS(322) CONSTS(330) = ( ( CONSTS(11)*(1.00000 - CONSTS(261)))*CONSTS(297))*CONSTS(327) CONSTS(331) = ( CONSTS(41)*CONSTS(295))*CONSTS(327) CONSTS(332) = ( CONSTS(300)*CONSTS(298))*CONSTS(327) CONSTS(333) = CONSTS(303)/CONSTS(308) CONSTS(334) = (1.00000 - CONSTS(176)) - CONSTS(177) CONSTS(335) = CONSTS(176)+CONSTS(177) CONSTS(336) = ( CONSTS(328)*CONSTS(193))*CONSTS(322) CONSTS(337) = 0.200000*CONSTS(193) CONSTS(338) = CONSTS(315)/CONSTS(308) CONSTS(339) = CONSTS(185)*1.00000 ** CONSTS(182) CONSTS(340) = 1.00000 - CONSTS(179) CONSTS(341) = CONSTS(194) CONSTS(342) = CONSTS(303)/CONSTS(315) CONSTS(343) = (1.00000 - CONSTS(339)/(CONSTS(168)+CONSTS(339)))*CONSTS(174) CONSTS(344) = CONSTS(179) - CONSTS(178) CONSTS(345) = ( CONSTS(341)*CONSTS(337))*CONSTS(320) CONSTS(346) = CONSTS(195) CONSTS(347) = 87.8251 ** 4.00000 CONSTS(348) = ( CONSTS(343)*CONSTS(176))*CONSTS(320) CONSTS(349) = ( CONSTS(343)*CONSTS(334))*CONSTS(322) CONSTS(350) = ( CONSTS(343)*CONSTS(177))*CONSTS(327) CONSTS(351) = CONSTS(185)*1.00000 ** CONSTS(183) CONSTS(352) = ( CONSTS(346)*CONSTS(337))*CONSTS(327) CONSTS(353) = (1.00000 - CONSTS(341)) - CONSTS(346) CONSTS(354) = 1.00000+(CONSTS(119) - CONSTS(111))/CONSTS(117) CONSTS(355) = CONSTS(185)*1.00000 ** CONSTS(184) CONSTS(356) = ( CONSTS(353)*CONSTS(337))*CONSTS(322) CONSTS(357) = (CONSTS(117)/2.00000)*( (CONSTS(354) ** 2.00000+( 4.00000*CONSTS(111))/CONSTS(117)) ** (1.0 / 2) - CONSTS(354)) CONSTS(358) = CONSTS(181)/(1.00000 - CONSTS(181)) CONSTS(359) = 1.00000+(CONSTS(119) - CONSTS(83))/CONSTS(114) CONSTS(360) = CONSTS(358)*( CONSTS(179)*(((1.00000+CONSTS(190)) - CONSTS(190)*CONSTS(335)) - CONSTS(180))+ CONSTS(335)*(CONSTS(180) - 1.00000))+ ( CONSTS(190)*CONSTS(179))*(CONSTS(180) - CONSTS(335)) CONSTS(361) = (CONSTS(114)/2.00000)*( (CONSTS(359) ** 2.00000+( 4.00000*CONSTS(83))/CONSTS(114)) ** (1.0 / 2) - CONSTS(359)) CONSTS(362) = CONSTS(179)*((1.00000+CONSTS(190))+ CONSTS(180)*(CONSTS(358) - CONSTS(190))) - CONSTS(180)*(1.00000+CONSTS(358)) CONSTS(363) = CONSTS(138) ** 3.00000 CONSTS(364) = (CONSTS(360)/CONSTS(362))*CONSTS(174) CONSTS(365) = CONSTS(139) ** 3.00000 CONSTS(366) = ( (CONSTS(180) - CONSTS(335))*CONSTS(174)+ CONSTS(180)*CONSTS(364))/( (1.00000+CONSTS(190))*CONSTS(179) - CONSTS(180)) CONSTS(367) = CONSTS(101) ** 3.00000 CONSTS(368) = ( ( CONSTS(190)*CONSTS(179))*CONSTS(366))/CONSTS(364) CONSTS(369) = (1.00000 - CONSTS(355)/(CONSTS(170)+CONSTS(355)))*CONSTS(366) CONSTS(370) = CONSTS(160) - 0.000170000 CONSTS(371) = ( CONSTS(369)*CONSTS(178))*CONSTS(320) CONSTS(372) = ( CONSTS(369)*CONSTS(340))*CONSTS(322) CONSTS(373) = ( CONSTS(369)*CONSTS(344))*CONSTS(327) CONSTS(374) = 1.00000 - CONSTS(368) CONSTS(375) = CONSTS(160)*(1.00000 - 0.460000) CONSTS(376) = (1.00000 - CONSTS(351)/(CONSTS(169)+CONSTS(351)))*CONSTS(364) CONSTS(377) = CONSTS(375) CONSTS(378) = ( CONSTS(376)*CONSTS(368))*CONSTS(320) CONSTS(379) = ( CONSTS(376)*CONSTS(374))*CONSTS(322) CONSTS(380) = 3.25000*CONSTS(162) CONSTS(381) = ( CONSTS(194)*CONSTS(193))*CONSTS(320) CONSTS(382) = (((CONSTS(17)+CONSTS(30))+CONSTS(27))+CONSTS(28)) - CONSTS(381) CONSTS(383) = (( CONSTS(17)*CONSTS(28)+ CONSTS(30)*CONSTS(27))+ CONSTS(27)*CONSTS(28)) - CONSTS(381)*(CONSTS(27)+CONSTS(28)) CONSTS(384) = ( CONSTS(381)*CONSTS(27))*CONSTS(28) CONSTS(385) = ((( (- CONSTS(384)/27.0000)*CONSTS(382) ** 3.00000 - ( ( ( CONSTS(382)*CONSTS(382))*CONSTS(383))*CONSTS(383))/108.000)+( ( CONSTS(382)*CONSTS(383))*CONSTS(384))/6.00000)+CONSTS(383) ** 3.00000/27.0000)+( CONSTS(384)*CONSTS(384))/4.00000 CONSTS(386) = TERNRY(CONSTS(385).LT.0.00000, - CONSTS(385) ** (1.0 / 2), 0.00000) CONSTS(387) = (((TERNRY(CONSTS(385).GT.0.00000, CONSTS(385) ** (1.0 / 2), 0.00000))+CONSTS(384)/2.00000)+( CONSTS(382)*CONSTS(383))/6.00000) - CONSTS(382) ** 3.00000/27.0000 CONSTS(388) = atan(CONSTS(386)/CONSTS(387))/3.00000 CONSTS(389) = CONSTS(387)*CONSTS(387)+ CONSTS(386)*CONSTS(386) ** 1.00000/6.00000 CONSTS(390) = (CONSTS(383)/3.00000 - ( CONSTS(382)*CONSTS(382))/9.00000)/( CONSTS(389)*CONSTS(389)) CONSTS(391) = ( CONSTS(389)*cos(CONSTS(388)))*(1.00000 - CONSTS(390)) - CONSTS(382)/3.00000 CONSTS(392) = ((CONSTS(17) - CONSTS(18))+CONSTS(282))/( (CONSTS(292)/CONSTS(21)+1.00000)*(CONSTS(391)/CONSTS(27)+1.00000)) CONSTS(393) = ( ( CONSTS(282)*CONSTS(392))*CONSTS(391))/( CONSTS(25)*CONSTS(27)) CONSTS(394) = ( CONSTS(393)*CONSTS(292))/CONSTS(21) CONSTS(395) = ((CONSTS(30) - CONSTS(31))+CONSTS(294))/( (CONSTS(292)/CONSTS(24)+1.00000)*(CONSTS(391)/CONSTS(28)+1.00000)) CONSTS(396) = ( ( CONSTS(294)*CONSTS(395))*CONSTS(391))/( CONSTS(26)*CONSTS(28)) CONSTS(397) = ( CONSTS(396)*CONSTS(292))/CONSTS(24) CONSTS(398) = 0.0269000+CONSTS(393)/CONSTS(18) CONSTS(399) = 0.0329000+CONSTS(396)/CONSTS(31) CONSTS(400) = ( CONSTS(195)*CONSTS(193))*CONSTS(327) CONSTS(401) = 1.00000+(CONSTS(119) - CONSTS(400))/CONSTS(118) CONSTS(402) = (CONSTS(118)/2.00000)*( (CONSTS(401) ** 2.00000+( 4.00000*CONSTS(400))/CONSTS(118)) ** (1.0 / 2) - CONSTS(401)) CONSTS(403) = ((CONSTS(119) - CONSTS(111))+CONSTS(357))/( (1.00000+CONSTS(361)/CONSTS(114))*(1.00000+CONSTS(402)/CONSTS(118))) CONSTS(404) = ( ( CONSTS(357)*CONSTS(403))*CONSTS(402))/( CONSTS(117)*CONSTS(118)) CONSTS(405) = 0.0306000+CONSTS(404)/CONSTS(111) CONSTS(406) = ( CONSTS(404)*CONSTS(361))/CONSTS(114) CONSTS(407) = - CONSTS(303)/CONSTS(308) CONSTS(408) = - CONSTS(303)/CONSTS(308) CONSTS(409) = CONSTS(315)/CONSTS(308) CONSTS(410) = ( 3.00000*CONSTS(303))/CONSTS(315) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(8) = 0.00100000*( CONSTS(220)*STATES(9) - ( CONSTS(222)*STATES(1))*STATES(8)) RATES(10) = 0.00100000*( CONSTS(56)*STATES(2) - ( CONSTS(57)*STATES(11))*STATES(10)) RATES(16) = 0.00100000*( CONSTS(56)*STATES(3) - ( CONSTS(57)*STATES(17))*STATES(16)) RATES(21) = 0.00100000*( CONSTS(220)*STATES(22) - ( CONSTS(222)*STATES(23))*STATES(21)) RATES(24) = 0.00100000*( CONSTS(56)*STATES(4) - ( CONSTS(57)*STATES(25))*STATES(24)) RATES(32) = 0.00100000*(( ( CONSTS(75)*STATES(20))*(1.00000 - STATES(32)))/(CONSTS(65)+ (1.00000 - STATES(32))*1.00000) - ( ( CONSTS(77)*CONSTS(67))*STATES(32))/(CONSTS(66)+ STATES(32)*1.00000)) RATES(105) = 0.00100000*(( ( CONSTS(125)*STATES(30))*(1.00000 - STATES(105)))/(CONSTS(122)+ (1.00000 - STATES(105))*1.00000) - ( ( CONSTS(126)*CONSTS(83))*STATES(105))/(CONSTS(123)+ STATES(105)*1.00000)) RATES(106) = 0.00100000*(( ( CONSTS(129)*STATES(7))*(1.00000 - STATES(106)))/(CONSTS(127)+ (1.00000 - STATES(106))*1.00000) - ( ( CONSTS(130)*CONSTS(29))*STATES(106))/(CONSTS(128)+ STATES(106)*1.00000)) RATES(116) = 0.00100000*(( ( CONSTS(146)*STATES(7))*(1.00000 - STATES(116)))/(CONSTS(143)+ (1.00000 - STATES(116))*1.00000) - ( ( CONSTS(150)*CONSTS(29))*STATES(116))/(CONSTS(144)+ STATES(116)*1.00000)) RATES(128) = 0.00100000*( ( CONSTS(189)*STATES(7))*(CONSTS(378) - STATES(128)) - CONSTS(245)*STATES(128)) RATES(129) = 0.00100000*( ( CONSTS(189)*STATES(20))*(CONSTS(379) - STATES(129)) - CONSTS(245)*STATES(129)) RATES(130) = 0.00100000*( ( CONSTS(189)*STATES(7))*(CONSTS(371) - STATES(130)) - CONSTS(245)*STATES(130)) RATES(131) = 0.00100000*( ( CONSTS(189)*STATES(20))*(CONSTS(372) - STATES(131)) - CONSTS(245)*STATES(131)) RATES(132) = 0.00100000*( ( CONSTS(189)*STATES(30))*(CONSTS(373) - STATES(132)) - CONSTS(245)*STATES(132)) RATES(133) = 0.00100000*( CONSTS(202)*STATES(134) - ( CONSTS(249)*STATES(133))*STATES(7)) RATES(134) = 0.00100000*(( ( CONSTS(201)*STATES(135))*STATES(44) - (CONSTS(248)+CONSTS(202))*STATES(134))+ ( CONSTS(249)*STATES(133))*STATES(7)) RATES(7) = 0.00100000*((( CONSTS(202)*STATES(134) - ( CONSTS(249)*STATES(133))*STATES(7))+ CONSTS(265)*STATES(136)) - ( CONSTS(196)*(CONSTS(345) - STATES(136)))*STATES(7)) RATES(136) = 0.00100000*( ( CONSTS(196)*(CONSTS(345) - STATES(136)))*STATES(7) - CONSTS(265)*STATES(136)) RATES(137) = 0.00100000*( CONSTS(252)*STATES(138) - ( CONSTS(268)*STATES(137))*STATES(20)) RATES(138) = 0.00100000*(( ( CONSTS(251)*STATES(139))*STATES(45) - (CONSTS(267)+CONSTS(252))*STATES(138))+ ( CONSTS(268)*STATES(137))*STATES(20)) RATES(20) = 0.00100000*((( CONSTS(252)*STATES(138) - ( CONSTS(268)*STATES(137))*STATES(20))+ CONSTS(265)*STATES(140)) - ( CONSTS(196)*(CONSTS(356) - STATES(140)))*STATES(20)) RATES(140) = 0.00100000*( ( CONSTS(196)*(CONSTS(356) - STATES(140)))*STATES(20) - CONSTS(265)*STATES(140)) RATES(141) = 0.00100000*( CONSTS(271)*STATES(142) - ( CONSTS(279)*STATES(141))*STATES(30)) RATES(142) = 0.00100000*(( ( CONSTS(270)*STATES(143))*STATES(46) - (CONSTS(278)+CONSTS(271))*STATES(142))+ ( CONSTS(279)*STATES(141))*STATES(30)) RATES(30) = 0.00100000*((( CONSTS(271)*STATES(142) - ( CONSTS(279)*STATES(141))*STATES(30))+ CONSTS(265)*STATES(144)) - ( CONSTS(196)*(CONSTS(352) - STATES(144)))*STATES(30)) RATES(144) = 0.00100000*( ( CONSTS(196)*(CONSTS(352) - STATES(144)))*STATES(30) - CONSTS(265)*STATES(144)) ALGBRC(4) = CONSTS(394) - STATES(5) RATES(5) = 0.00100000*(( ( CONSTS(32)*STATES(7))*ALGBRC(4))/(CONSTS(19)+ALGBRC(4)) - ( ( CONSTS(34)*CONSTS(29))*STATES(5))/(CONSTS(22)+STATES(5))) ALGBRC(5) = CONSTS(397) - STATES(6) RATES(6) = 0.00100000*(( ( CONSTS(33)*STATES(7))*ALGBRC(5))/(CONSTS(20)+ALGBRC(5)) - ( ( CONSTS(35)*CONSTS(29))*STATES(6))/(CONSTS(23)+STATES(6))) ALGBRC(16) = CONSTS(406) - STATES(104) RATES(104) = 0.00100000*(( ( CONSTS(120)*STATES(30))*ALGBRC(16))/(CONSTS(115)+ALGBRC(16)) - ( ( CONSTS(121)*CONSTS(83))*STATES(104))/(CONSTS(116)+STATES(104))) ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(51)+91.0000)/6.10000)) RATES(117) = (ALGBRC(26) - STATES(117))/CONSTS(151) ALGBRC(30) = CONSTS(256) - STATES(145) RATES(145) = 0.00100000*(( ( CONSTS(212)*STATES(20))*ALGBRC(30))/(CONSTS(208)+ALGBRC(30)) - ( ( CONSTS(211)*CONSTS(67))*STATES(145))/(CONSTS(207)+STATES(145))) ALGBRC(34) = (CONSTS(323) - STATES(13)) - STATES(12) RATES(13) = 0.00100000*( ( CONSTS(221)*STATES(7))*ALGBRC(34) - CONSTS(258)*STATES(13)) ALGBRC(36) = (CONSTS(326) - STATES(19)) - STATES(18) RATES(19) = 0.00100000*( ( CONSTS(221)*STATES(20))*ALGBRC(36) - CONSTS(258)*STATES(19)) ALGBRC(37) = (CONSTS(331) - STATES(27)) - STATES(26) RATES(27) = 0.00100000*( ( CONSTS(221)*STATES(30))*ALGBRC(37) - CONSTS(258)*STATES(27)) ALGBRC(12) = 0.0250000/(1.00000+EXP((STATES(51)+58.0000)/5.00000)) ALGBRC(40) = 0.200000/(1.00000+EXP((STATES(51)+19.0000)/- 9.00000)) RATES(68) = (ALGBRC(12)/(ALGBRC(12)+ALGBRC(40)) - STATES(68))/CONSTS(110) ALGBRC(13) = 1.00000/(( 0.000600000*(STATES(51) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(51) - 1.73840))) - ( 0.000300000*(STATES(51)+38.3608))/(1.00000 - EXP( 0.152200*(STATES(51)+38.3608)))) ALGBRC(41) = 1.00000/(1.00000+EXP((STATES(51)+10.0850)/- 4.25000)) RATES(69) = (ALGBRC(41) - STATES(69))/ALGBRC(13) ALGBRC(14) = 0.00739900/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0311960))/- 0.800190)) ALGBRC(42) = 0.00569920/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0415200))/1.34890)) RATES(74) = ALGBRC(42)*STATES(81) - STATES(74)*( 4.00000*ALGBRC(14)) ALGBRC(15) = 0.00994150/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0448090))/- 0.581720)) ALGBRC(43) = 0.00332010/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0942170))/0.953640)) RATES(89) = ALGBRC(43)*STATES(96) - STATES(89)*( 4.00000*ALGBRC(15)) ALGBRC(17) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 0.00000, 0.135000*EXP(((87.0000+STATES(51))+CONSTS(131))/- 6.80000)) ALGBRC(44) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(((STATES(51)+CONSTS(131))+27.4034)/- 11.1000))), 3.56000*EXP( 0.0790000*((STATES(51)+CONSTS(131))+7.00000))+ 310000.*EXP( 0.350000*((STATES(51)+CONSTS(131))+7.00000))) RATES(107) = ALGBRC(17)*(1.00000 - STATES(107)) - ALGBRC(44)*STATES(107) ALGBRC(18) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*(STATES(51)+CONSTS(131)))+ - 6.94800e-05*EXP( - 0.0439100*(STATES(51)+CONSTS(131))))*((STATES(51)+CONSTS(131))+37.7800))/(1.00000+EXP( 0.311000*((STATES(51)+CONSTS(131))+79.2300)))) ALGBRC(45) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*(STATES(51)+CONSTS(131))))/(1.00000+EXP( - 0.100000*((STATES(51)+CONSTS(131))+32.0000))), ( 0.121200*EXP( - 0.0105200*(STATES(51)+CONSTS(131))))/(1.00000+EXP( - 0.137800*((STATES(51)+CONSTS(131))+40.1400)))) RATES(108) = ALGBRC(18)*(1.00000 - STATES(108)) - ALGBRC(45)*STATES(108) ALGBRC(19) = ( 0.320000*(STATES(51)+58.4729))/(1.00000 - EXP( - 0.100000*(STATES(51)+58.4729))) ALGBRC(46) = 0.0800000*EXP((13.7299 - STATES(51))/11.0000) RATES(109) = ALGBRC(19)*(1.00000 - STATES(109)) - ALGBRC(46)*STATES(109) ALGBRC(20) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, 0.135000*EXP((87.0000+STATES(51))/- 6.80000)) ALGBRC(47) = TERNRY(STATES(51).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(51)+27.4034)/- 11.1000))), 3.56000*EXP( 0.0790000*(STATES(51)+7.00000))+ 310000.*EXP( 0.350000*(STATES(51)+7.00000))) RATES(110) = ALGBRC(20)*(1.00000 - STATES(110)) - ALGBRC(47)*STATES(110) ALGBRC(21) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(51)) - 6.94800e-05*EXP( - 0.0439100*STATES(51)))*(STATES(51)+37.7800))/(1.00000+EXP( 0.311000*(STATES(51)+79.2300)))) ALGBRC(48) = TERNRY(STATES(51).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(51)))/(1.00000+EXP( - 0.100000*(STATES(51)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(51)))/(1.00000+EXP( - 0.137800*(STATES(51)+40.1400)))) RATES(111) = ALGBRC(21)*(1.00000 - STATES(111)) - ALGBRC(48)*STATES(111) ALGBRC(22) = ( 0.320000*(STATES(51)+58.4729))/(1.00000 - EXP( - 0.100000*(STATES(51)+58.4729))) ALGBRC(49) = 0.0800000*EXP((13.7299 - STATES(51))/11.0000) RATES(112) = ALGBRC(22)*(1.00000 - STATES(112)) - ALGBRC(49)*STATES(112) ALGBRC(23) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, 0.135000*EXP(((87.0000+STATES(51))+CONSTS(133))/- 6.80000)) ALGBRC(50) = TERNRY(STATES(51).GE.- 40.0000, (1.00000/0.130000)/(1.00000+EXP(((STATES(51)+CONSTS(133))+27.4034)/- 11.1000)), 3.56000*EXP( 0.0790000*((STATES(51)+CONSTS(133))+7.00000))+ 310000.*EXP( 0.350000*((STATES(51)+CONSTS(133))+7.00000))) RATES(113) = ALGBRC(23)*(1.00000 - STATES(113)) - ALGBRC(50)*STATES(113) ALGBRC(24) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*(STATES(51)+CONSTS(133))) - 6.94800e-05*EXP( - 0.0439100*(STATES(51)+CONSTS(133))))*((STATES(51)+CONSTS(133))+37.7800))/(1.00000+EXP( 0.311000*((STATES(51)+CONSTS(133))+79.2300)))) ALGBRC(51) = TERNRY(STATES(51).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*(STATES(51)+CONSTS(133))))/(1.00000+EXP( - 0.100000*((STATES(51)+CONSTS(133))+32.0000))), ( 0.121200*EXP( - 0.0105200*(STATES(51)+CONSTS(133))))/(1.00000+EXP( - 0.137800*((STATES(51)+CONSTS(133))+40.1400)))) RATES(114) = ALGBRC(24)*(1.00000 - STATES(114)) - ALGBRC(51)*STATES(114) ALGBRC(25) = ( 0.320000*((STATES(51)+CONSTS(132))+58.4729))/(1.00000 - EXP( - 0.100000*((STATES(51)+CONSTS(132))+58.4729))) ALGBRC(52) = 0.0800000*EXP(((STATES(51)+CONSTS(132)) - 13.7299)/- 11.0000) RATES(115) = ALGBRC(25)*(1.00000 - STATES(115)) - ALGBRC(52)*STATES(115) ALGBRC(27) = ( 0.320000*(STATES(51)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(51)+47.1300))) ALGBRC(53) = 0.0800000*EXP(STATES(51)/- 11.0000) RATES(118) = ALGBRC(27)*(1.00000 - STATES(118)) - ALGBRC(53)*STATES(118) ALGBRC(58) = (CONSTS(324) - STATES(15)) - STATES(14) RATES(15) = 0.00100000*( ( CONSTS(221)*STATES(7))*ALGBRC(58) - CONSTS(258)*STATES(15)) ALGBRC(61) = (CONSTS(332) - STATES(29)) - STATES(28) RATES(29) = 0.00100000*( ( CONSTS(221)*STATES(30))*ALGBRC(61) - CONSTS(258)*STATES(29)) ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(51)+9.43700)/- 7.13300)) ALGBRC(54) = 1.00000/(1.00000+EXP((STATES(51) - 18.4099)/- 29.3814)) ALGBRC(66) = 1.00000/(1.00000+EXP((STATES(51)+100.000)/29.3814)) ALGBRC(78) = 1.00000/(ALGBRC(54)/1.20890+ 3.50000*ALGBRC(66)) RATES(121) = (ALGBRC(28) - STATES(121))/ALGBRC(78) ALGBRC(67) = ((1.00000/(1.00000+EXP((STATES(51)+19.0000)/- 9.00000)))/0.500000)/9.79530 ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(51)+58.0000)/5.00000)) ALGBRC(79) = 0.0479600*ALGBRC(29) RATES(122) = ALGBRC(79)*(1.00000 - STATES(122)) - ALGBRC(67)*STATES(122) ALGBRC(80) = 0.0214400*ALGBRC(29) RATES(123) = ALGBRC(80)*(1.00000 - STATES(123)) - ALGBRC(67)*STATES(123) ALGBRC(55) = (1.00000/(1.00000+EXP((STATES(51)+60.0000)/5.00000)))/250.000 ALGBRC(81) = 2.46000*ALGBRC(55) RATES(124) = ALGBRC(81)*(1.00000 - STATES(124)) - ALGBRC(67)*STATES(124) ALGBRC(82) = 0.560340*ALGBRC(55) RATES(125) = ALGBRC(82)*(1.00000 - STATES(125)) - ALGBRC(67)*STATES(125) ALGBRC(64) = 0.0906540*EXP( ( - 0.111570*STATES(51))*CONSTS(264)) ALGBRC(92) = 0.388390/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.150190))/- 0.606930)) RATES(75) = ( ALGBRC(42)*STATES(76)+ ALGBRC(92)*STATES(86)) - STATES(75)*( 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(64)) RATES(76) = ((( ( 2.00000*ALGBRC(14))*STATES(75)+ ( 2.00000*ALGBRC(42))*STATES(77))+ ( 2.00000*ALGBRC(92))*STATES(87))+ ( 3.00000*ALGBRC(64))*STATES(78)) - STATES(76)*(((ALGBRC(14)+ALGBRC(42))+ALGBRC(92))+ 2.00000*ALGBRC(64)) RATES(77) = (( ALGBRC(14)*STATES(76)+ ( 3.00000*ALGBRC(92))*STATES(88))+ ( 3.00000*ALGBRC(64))*STATES(79)) - STATES(77)*(( 2.00000*ALGBRC(42)+ 2.00000*ALGBRC(92))+ 2.00000*ALGBRC(64)) RATES(78) = ( ALGBRC(42)*STATES(79)+ ALGBRC(92)*STATES(76)) - STATES(78)*(ALGBRC(14)+ 3.00000*ALGBRC(64)) RATES(79) = (( ALGBRC(14)*STATES(78)+ ( 2.00000*ALGBRC(92))*STATES(77))+ ( 4.00000*ALGBRC(64))*STATES(80)) - STATES(79)*((ALGBRC(42)+ALGBRC(92))+ 3.00000*ALGBRC(64)) ALGBRC(76) = 0.00311240+(0.0283300 - 0.00311240)/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.0516600))/1.55220)) RATES(80) = ( ALGBRC(92)*STATES(79) - STATES(80)*( 4.00000*ALGBRC(64)+CONSTS(112)))+ ALGBRC(76)*STATES(70) RATES(81) = (( ( 4.00000*ALGBRC(14))*STATES(74)+ ( 2.00000*ALGBRC(42))*STATES(82))+ ALGBRC(64)*STATES(85)) - STATES(81)*(( 3.00000*ALGBRC(14)+ALGBRC(42))+ALGBRC(92)) RATES(82) = (( ( 3.00000*ALGBRC(14))*STATES(81)+ ( 3.00000*ALGBRC(42))*STATES(83))+ ALGBRC(64)*STATES(86)) - STATES(82)*(( 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(42))+ 2.00000*ALGBRC(92)) RATES(83) = (( ( 2.00000*ALGBRC(14))*STATES(82)+ ( 4.00000*ALGBRC(42))*STATES(84))+ ALGBRC(64)*STATES(87)) - STATES(83)*((ALGBRC(14)+ 3.00000*ALGBRC(42))+ 3.00000*ALGBRC(92)) RATES(84) = ( ALGBRC(14)*STATES(83)+ ALGBRC(64)*STATES(88)) - STATES(84)*( 4.00000*ALGBRC(42)+ 4.00000*ALGBRC(92)) RATES(85) = ( ALGBRC(42)*STATES(86)+ ALGBRC(92)*STATES(81)) - STATES(85)*( 3.00000*ALGBRC(14)+ALGBRC(64)) RATES(86) = ((( ( 3.00000*ALGBRC(14))*STATES(85)+ ( 2.00000*ALGBRC(42))*STATES(87))+ ( 2.00000*ALGBRC(92))*STATES(82))+ ( 2.00000*ALGBRC(64))*STATES(75)) - STATES(86)*((( 2.00000*ALGBRC(14)+ALGBRC(42))+ALGBRC(92))+ALGBRC(64)) RATES(87) = ((( ( 2.00000*ALGBRC(14))*STATES(86)+ ( 3.00000*ALGBRC(42))*STATES(88))+ ( 3.00000*ALGBRC(92))*STATES(83))+ ( 2.00000*ALGBRC(64))*STATES(76)) - STATES(87)*(((ALGBRC(14)+ 2.00000*ALGBRC(42))+ 2.00000*ALGBRC(92))+ALGBRC(64)) RATES(88) = (( ALGBRC(14)*STATES(87)+ ( 4.00000*ALGBRC(92))*STATES(84))+ ( 2.00000*ALGBRC(64))*STATES(77)) - STATES(88)*(( 3.00000*ALGBRC(42)+ 3.00000*ALGBRC(92))+ALGBRC(64)) ALGBRC(65) = 0.0657000*EXP( ( - 0.118990*STATES(51))*CONSTS(264)) ALGBRC(94) = 0.563560/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.179860))/- 0.583810)) RATES(90) = ( ALGBRC(43)*STATES(91)+ ALGBRC(94)*STATES(101)) - STATES(90)*( 2.00000*ALGBRC(15)+ 2.00000*ALGBRC(65)) RATES(91) = ((( ( 2.00000*ALGBRC(15))*STATES(90)+ ( 2.00000*ALGBRC(43))*STATES(92))+ ( 2.00000*ALGBRC(94))*STATES(102))+ ( 3.00000*ALGBRC(65))*STATES(93)) - STATES(91)*(((ALGBRC(15)+ALGBRC(43))+ALGBRC(94))+ 2.00000*ALGBRC(65)) RATES(92) = (( ALGBRC(15)*STATES(91)+ ( 3.00000*ALGBRC(94))*STATES(103))+ ( 3.00000*ALGBRC(65))*STATES(94)) - STATES(92)*(( 2.00000*ALGBRC(43)+ 2.00000*ALGBRC(94))+ 2.00000*ALGBRC(65)) RATES(93) = ( ALGBRC(43)*STATES(94)+ ALGBRC(94)*STATES(91)) - STATES(93)*(ALGBRC(15)+ 3.00000*ALGBRC(65)) RATES(94) = (( ALGBRC(15)*STATES(93)+ ( 2.00000*ALGBRC(94))*STATES(92))+ ( 4.00000*ALGBRC(65))*STATES(95)) - STATES(94)*((ALGBRC(43)+ALGBRC(94))+ 3.00000*ALGBRC(65)) ALGBRC(77) = 0.000385250+(0.0124060 - 0.000385250)/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.0641180))/0.779920)) RATES(95) = ( ALGBRC(94)*STATES(94) - STATES(95)*( 4.00000*ALGBRC(65)+CONSTS(113)))+ ALGBRC(77)*STATES(72) RATES(96) = (( ( 4.00000*ALGBRC(15))*STATES(89)+ ( 2.00000*ALGBRC(43))*STATES(97))+ ALGBRC(65)*STATES(100)) - STATES(96)*(( 3.00000*ALGBRC(15)+ALGBRC(43))+ALGBRC(94)) RATES(97) = (( ( 3.00000*ALGBRC(15))*STATES(96)+ ( 3.00000*ALGBRC(43))*STATES(98))+ ALGBRC(65)*STATES(101)) - STATES(97)*(( 2.00000*ALGBRC(15)+ 2.00000*ALGBRC(43))+ 2.00000*ALGBRC(94)) RATES(98) = (( ( 2.00000*ALGBRC(15))*STATES(97)+ ( 4.00000*ALGBRC(43))*STATES(99))+ ALGBRC(65)*STATES(102)) - STATES(98)*((ALGBRC(15)+ 3.00000*ALGBRC(43))+ 3.00000*ALGBRC(94)) RATES(99) = ( ALGBRC(15)*STATES(98)+ ALGBRC(65)*STATES(103)) - STATES(99)*( 4.00000*ALGBRC(43)+ 4.00000*ALGBRC(94)) RATES(100) = ( ALGBRC(43)*STATES(101)+ ALGBRC(94)*STATES(96)) - STATES(100)*( 3.00000*ALGBRC(15)+ALGBRC(65)) RATES(101) = ((( ( 3.00000*ALGBRC(15))*STATES(100)+ ( 2.00000*ALGBRC(43))*STATES(102))+ ( 2.00000*ALGBRC(94))*STATES(97))+ ( 2.00000*ALGBRC(65))*STATES(90)) - STATES(101)*((( 2.00000*ALGBRC(15)+ALGBRC(43))+ALGBRC(94))+ALGBRC(65)) RATES(102) = ((( ( 2.00000*ALGBRC(15))*STATES(101)+ ( 3.00000*ALGBRC(43))*STATES(103))+ ( 3.00000*ALGBRC(94))*STATES(98))+ ( 2.00000*ALGBRC(65))*STATES(91)) - STATES(102)*(((ALGBRC(15)+ 2.00000*ALGBRC(43))+ 2.00000*ALGBRC(94))+ALGBRC(65)) RATES(103) = (( ALGBRC(15)*STATES(102)+ ( 4.00000*ALGBRC(94))*STATES(99))+ ( 2.00000*ALGBRC(65))*STATES(92)) - STATES(103)*(( 3.00000*ALGBRC(43)+ 3.00000*ALGBRC(94))+ALGBRC(65)) ALGBRC(93) = 0.000441980*EXP( ( - 1.20220*STATES(51))*CONSTS(264)) ALGBRC(104) = 0.000401730*EXP( ( 0.000208730*STATES(51))*CONSTS(264)) RATES(70) = ( - (ALGBRC(76)+ALGBRC(104))*STATES(70)+ ALGBRC(93)*STATES(71))+ CONSTS(112)*STATES(80) RATES(71) = ALGBRC(104)*STATES(70) - ALGBRC(93)*STATES(71) ALGBRC(95) = 0.000237300*EXP( ( - 1.97420*STATES(51))*CONSTS(264)) ALGBRC(105) = 0.000226520*EXP( ( 0.000246890*STATES(51))*CONSTS(264)) RATES(72) = ( - (ALGBRC(77)+ALGBRC(105))*STATES(72)+ ALGBRC(95)*STATES(73))+ CONSTS(113)*STATES(95) RATES(73) = ALGBRC(105)*STATES(72) - ALGBRC(95)*STATES(73) ALGBRC(6) = ( ( CONSTS(37)*CONSTS(36))*CONSTS(320) - STATES(9)) - STATES(8) ALGBRC(35) = ( ALGBRC(6)*(CONSTS(60)+CONSTS(52)) - STATES(15)*(CONSTS(52)+CONSTS(60)))+ ( CONSTS(50)*CONSTS(52))*(1.00000+CONSTS(60)/CONSTS(55)) ALGBRC(59) = ( - STATES(15)*CONSTS(50))*CONSTS(52) ALGBRC(71) = (- ALGBRC(35)+ ( ALGBRC(35)*ALGBRC(35) - ( 4.00000*CONSTS(223))*ALGBRC(59)) ** (1.0 / 2))/( 2.00000*CONSTS(223)) ALGBRC(87) = ALGBRC(6)/(1.00000+ (ALGBRC(71)/CONSTS(50))*(1.00000+CONSTS(60)/CONSTS(52))) ALGBRC(99) = ( ALGBRC(71)*ALGBRC(87))/CONSTS(50) ALGBRC(109) = ( ALGBRC(99)*CONSTS(60))/CONSTS(52) RATES(9) = 0.00100000*(( CONSTS(45)*ALGBRC(99)+ CONSTS(43)*ALGBRC(109)) - CONSTS(220)*STATES(9)) RATES(1) = 0.00100000*(( CONSTS(45)*ALGBRC(99)+ CONSTS(43)*ALGBRC(109)) - ( CONSTS(222)*STATES(1))*STATES(8)) ALGBRC(9) = ( ( CONSTS(216)*CONSTS(36))*CONSTS(327) - STATES(22)) - STATES(21) ALGBRC(38) = ( ALGBRC(9)*(CONSTS(60)+CONSTS(52)) - STATES(29)*(CONSTS(52)+CONSTS(60)))+ ( CONSTS(50)*CONSTS(52))*(1.00000+CONSTS(60)/CONSTS(55)) ALGBRC(62) = ( - STATES(29)*CONSTS(50))*CONSTS(52) ALGBRC(74) = (- ALGBRC(38)+ ( ALGBRC(38)*ALGBRC(38) - ( 4.00000*CONSTS(226))*ALGBRC(62)) ** (1.0 / 2))/( 2.00000*CONSTS(226)) ALGBRC(90) = ALGBRC(9)/(1.00000+ (ALGBRC(74)/CONSTS(50))*(1.00000+CONSTS(60)/CONSTS(52))) ALGBRC(102) = ( ALGBRC(74)*ALGBRC(90))/CONSTS(50) ALGBRC(112) = (CONSTS(60)/CONSTS(52))*ALGBRC(102) RATES(22) = 0.00100000*(( CONSTS(45)*ALGBRC(102)+ CONSTS(43)*ALGBRC(112)) - CONSTS(220)*STATES(22)) RATES(23) = 0.00100000*(( CONSTS(45)*ALGBRC(102)+ CONSTS(43)*ALGBRC(112)) - ( CONSTS(222)*STATES(23))*STATES(21)) ALGBRC(8) = ( ( CONSTS(217)*CONSTS(296))*CONSTS(322) - STATES(3)) - STATES(16) ALGBRC(60) = ( ALGBRC(8)*(CONSTS(48)+CONSTS(60)) - ALGBRC(36)*(CONSTS(48)+CONSTS(60)))+ ( CONSTS(47)*CONSTS(48))*(1.00000+CONSTS(60)/CONSTS(49)) ALGBRC(72) = ( - ALGBRC(36)*CONSTS(48))*CONSTS(47) ALGBRC(88) = (- ALGBRC(60)+ ( ALGBRC(60)*ALGBRC(60) - ( 4.00000*CONSTS(225))*ALGBRC(72)) ** (1.0 / 2))/( 2.00000*CONSTS(225)) ALGBRC(100) = ALGBRC(8)/(1.00000+ (ALGBRC(88)/CONSTS(47))*(1.00000+CONSTS(60)/CONSTS(48))) ALGBRC(111) = ( ( CONSTS(60)*ALGBRC(88))*ALGBRC(100))/( CONSTS(47)*CONSTS(48)) ALGBRC(120) = ( ALGBRC(100)*ALGBRC(88))/CONSTS(47) RATES(3) = 0.00100000*(( CONSTS(46)*ALGBRC(120)+ CONSTS(44)*ALGBRC(111)) - CONSTS(56)*STATES(3)) RATES(17) = 0.00100000*(( CONSTS(46)*ALGBRC(120)+ CONSTS(44)*ALGBRC(111)) - ( CONSTS(57)*STATES(17))*STATES(16)) ALGBRC(119) = ( CONSTS(60)*ALGBRC(88))/CONSTS(49) RATES(18) = 0.00100000*( ( CONSTS(219)*CONSTS(62))*(ALGBRC(119)+ALGBRC(111)) - CONSTS(218)*STATES(18)) ALGBRC(131) = ((CONSTS(400) - STATES(143)) - STATES(142)) - STATES(141) RATES(143) = 0.00100000*((( ( CONSTS(269)*ALGBRC(131))*STATES(46) - CONSTS(277)*STATES(143)) - ( CONSTS(270)*STATES(143))*STATES(46))+ CONSTS(278)*STATES(142)) ALGBRC(2) = STATES(4)*1.00000 ** CONSTS(13) ALGBRC(32) = CONSTS(7)*(CONSTS(9)+ALGBRC(2)/(CONSTS(4)+ALGBRC(2))) ALGBRC(68) = ( ALGBRC(32)*CONSTS(330))*CONSTS(228) ALGBRC(84) = ( CONSTS(237)*(STATES(44) - STATES(46)))/CONSTS(310) ALGBRC(97) = ( CONSTS(238)*(STATES(46) - STATES(45)))/CONSTS(310) ALGBRC(108) = ( CONSTS(350)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(46)) ALGBRC(117) = ( (CONSTS(373)+ (CONSTS(175) - 1.00000)*STATES(132))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(46)) ALGBRC(125) = ALGBRC(108)+ALGBRC(117) ALGBRC(135) = (( ( - CONSTS(269)*ALGBRC(131))*STATES(46)+ CONSTS(277)*STATES(143)) - ( CONSTS(270)*STATES(143))*STATES(46))+ CONSTS(278)*STATES(142) RATES(46) = 0.00100000*((((ALGBRC(135)+ALGBRC(68)) - ALGBRC(125))+ALGBRC(84)) - ALGBRC(97)) ALGBRC(141) = ((CONSTS(381) - STATES(135)) - STATES(134)) - STATES(133) RATES(135) = 0.00100000*((( ( CONSTS(200)*ALGBRC(141))*STATES(44) - CONSTS(247)*STATES(135)) - ( CONSTS(201)*STATES(135))*STATES(44))+ CONSTS(248)*STATES(134)) ALGBRC(3) = STATES(2)*1.00000 ** CONSTS(14) ALGBRC(33) = STATES(2)*1.00000 ** CONSTS(15) ALGBRC(57) = ( CONSTS(8)*(CONSTS(10)+ALGBRC(3)/(CONSTS(5)+ALGBRC(3))))*(1.00000 - ( (1.00000 - ( CONSTS(16)*ALGBRC(33))/(CONSTS(6)+ALGBRC(33)))*STATES(1))/(CONSTS(3)+STATES(1))) ALGBRC(83) = ( ALGBRC(57)*CONSTS(321))*CONSTS(228) ALGBRC(96) = ( CONSTS(237)*(STATES(44) - STATES(46)))/CONSTS(307) ALGBRC(106) = ( CONSTS(236)*(STATES(44) - STATES(45)))/CONSTS(307) ALGBRC(115) = ( CONSTS(348)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(44)) ALGBRC(123) = ( (CONSTS(378)+ (CONSTS(175) - 1.00000)*STATES(128))*CONSTS(187))/(1.00000+CONSTS(172)/STATES(44)) ALGBRC(129) = ( (CONSTS(371)+ (CONSTS(175) - 1.00000)*STATES(130))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(44)) ALGBRC(136) = (ALGBRC(115)+ALGBRC(123))+ALGBRC(129) ALGBRC(146) = (( ( - CONSTS(200)*ALGBRC(141))*STATES(44)+ CONSTS(247)*STATES(135)) - ( CONSTS(201)*STATES(135))*STATES(44))+ CONSTS(248)*STATES(134) RATES(44) = 0.00100000*((((ALGBRC(146)+ALGBRC(83)) - ALGBRC(136)) - ALGBRC(96)) - ALGBRC(106)) ALGBRC(151) = ((CONSTS(336) - STATES(139)) - STATES(138)) - STATES(137) RATES(139) = 0.00100000*((( ( CONSTS(250)*ALGBRC(151))*STATES(45) - CONSTS(266)*STATES(139)) - ( CONSTS(251)*STATES(139))*STATES(45))+ CONSTS(267)*STATES(138)) ALGBRC(1) = STATES(3)*1.00000 ** CONSTS(13) ALGBRC(56) = CONSTS(7)*(CONSTS(9)+ALGBRC(1)/(CONSTS(4)+ALGBRC(1))) ALGBRC(69) = ( ALGBRC(56)*CONSTS(325))*CONSTS(228) ALGBRC(31) = STATES(3)*1.00000 ** CONSTS(14) ALGBRC(85) = CONSTS(8)*(CONSTS(10)+ALGBRC(31)/(CONSTS(5)+ALGBRC(31))) ALGBRC(107) = ( ALGBRC(85)*CONSTS(329))*CONSTS(228) ALGBRC(116) = ( CONSTS(236)*(STATES(44) - STATES(45)))/CONSTS(308) ALGBRC(124) = ( CONSTS(238)*(STATES(46) - STATES(45)))/CONSTS(308) ALGBRC(130) = ( CONSTS(349)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(45)) ALGBRC(137) = ( (CONSTS(379)+ (CONSTS(175) - 1.00000)*STATES(129))*CONSTS(187))/(1.00000+CONSTS(172)/STATES(45)) ALGBRC(142) = ( (CONSTS(372)+ (CONSTS(175) - 1.00000)*STATES(131))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(45)) ALGBRC(147) = (ALGBRC(130)+ALGBRC(137))+ALGBRC(142) ALGBRC(155) = (( ( - CONSTS(250)*ALGBRC(151))*STATES(45)+ CONSTS(266)*STATES(139)) - ( CONSTS(251)*STATES(139))*STATES(45))+ CONSTS(267)*STATES(138) RATES(45) = 0.00100000*(((((ALGBRC(155)+ALGBRC(69))+ALGBRC(107)) - ALGBRC(147))+ALGBRC(116))+ALGBRC(124)) ALGBRC(7) = ( ( CONSTS(38)*CONSTS(296))*CONSTS(320) - STATES(2)) - STATES(10) ALGBRC(70) = (((CONSTS(259)+CONSTS(272))/CONSTS(224)+ALGBRC(34))+ALGBRC(58)) - ALGBRC(7) ALGBRC(86) = (( CONSTS(272)*(ALGBRC(34) - ALGBRC(7))+ CONSTS(259)*(ALGBRC(58) - ALGBRC(7)))+CONSTS(280))/CONSTS(224) ALGBRC(98) = ( ALGBRC(7)*CONSTS(280))/CONSTS(224) ALGBRC(110) = ((( (- ALGBRC(98)/27.0000)*ALGBRC(70) ** 3.00000 - ( ( ( ALGBRC(70)*ALGBRC(70))*ALGBRC(86))*ALGBRC(86))/108.000)+( ( ALGBRC(70)*ALGBRC(86))*ALGBRC(98))/6.00000)+ALGBRC(86) ** 3.00000/27.0000)+( ALGBRC(98)*ALGBRC(98))/4.00000 ALGBRC(118) = TERNRY(ALGBRC(110).LT.0.00000, - ALGBRC(110) ** (1.0 / 2), 0.00000) ALGBRC(126) = (((TERNRY(ALGBRC(110).GT.0.00000, ALGBRC(110) ** (1.0 / 2), 0.00000))+ALGBRC(98)/2.00000)+( ALGBRC(70)*ALGBRC(86))/6.00000) - ALGBRC(70) ** 3.00000/27.0000 ALGBRC(132) = atan(ALGBRC(118)/ALGBRC(126))/3.00000 ALGBRC(138) = ALGBRC(126)*ALGBRC(126)+ ALGBRC(118)*ALGBRC(118) ** 1.00000/6.00000 ALGBRC(143) = (ALGBRC(86)/3.00000 - ( ALGBRC(70)*ALGBRC(70))/9.00000)/( ALGBRC(138)*ALGBRC(138)) ALGBRC(148) = ( ALGBRC(138)* sin(ALGBRC(132)))*(1.00000+ALGBRC(143)) ALGBRC(152) = ( ALGBRC(138)*cos(ALGBRC(132)))*(1.00000 - ALGBRC(143)) - ALGBRC(70)/3.00000 ALGBRC(157) = ( ALGBRC(152)*ALGBRC(152)+ ALGBRC(148)*ALGBRC(148)) ** (1.0 / 2) ALGBRC(161) = ALGBRC(34)/((1.00000+CONSTS(60)/CONSTS(49))+( ALGBRC(157)*(CONSTS(48)+CONSTS(60)))/( CONSTS(47)*CONSTS(48))) ALGBRC(165) = ( CONSTS(60)*ALGBRC(161))/CONSTS(49) ALGBRC(169) = ( ( CONSTS(60)*ALGBRC(161))*ALGBRC(157))/( CONSTS(47)*CONSTS(48)) RATES(12) = 0.00100000*( ( CONSTS(219)*CONSTS(59))*(ALGBRC(165)+ALGBRC(169)) - CONSTS(218)*STATES(12)) ALGBRC(10) = ( ( CONSTS(39)*CONSTS(296))*CONSTS(327) - STATES(4)) - STATES(24) ALGBRC(73) = (((CONSTS(260)+CONSTS(273))/CONSTS(227)+ALGBRC(37))+ALGBRC(61)) - ALGBRC(10) ALGBRC(89) = (( CONSTS(273)*(ALGBRC(37) - ALGBRC(10))+ CONSTS(260)*(ALGBRC(61) - ALGBRC(10)))+CONSTS(281))/CONSTS(227) ALGBRC(101) = ( ALGBRC(10)*CONSTS(281))/CONSTS(227) ALGBRC(113) = ((( (- ALGBRC(101)/27.0000)*ALGBRC(73) ** 3.00000 - ( ( ( ALGBRC(73)*ALGBRC(73))*ALGBRC(89))*ALGBRC(89))/108.000)+( ( ALGBRC(73)*ALGBRC(89))*ALGBRC(101))/6.00000)+ALGBRC(89) ** 3.00000/27.0000)+( ALGBRC(101)*ALGBRC(101))/4.00000 ALGBRC(121) = TERNRY(ALGBRC(113).LT.0.00000, - ALGBRC(113) ** (1.0 / 2), 0.00000) ALGBRC(127) = (((TERNRY(ALGBRC(113).GT.0.00000, ALGBRC(113) ** (1.0 / 2), 0.00000))+ALGBRC(101)/2.00000)+( ALGBRC(73)*ALGBRC(89))/6.00000) - ALGBRC(73) ** 3.00000/27.0000 ALGBRC(133) = atan(ALGBRC(121)/ALGBRC(127))/3.00000 ALGBRC(139) = ALGBRC(127)*ALGBRC(127)+ ALGBRC(121)*ALGBRC(121) ** 1.00000/6.00000 ALGBRC(144) = (ALGBRC(89)/3.00000 - ( ALGBRC(73)*ALGBRC(73))/9.00000)/( ALGBRC(139)*ALGBRC(139)) ALGBRC(149) = ( ALGBRC(139)* sin(ALGBRC(133)))*(1.00000+ALGBRC(144)) ALGBRC(153) = ( ALGBRC(139)*cos(ALGBRC(133)))*(1.00000 - ALGBRC(144)) - ALGBRC(73)/3.00000 ALGBRC(158) = ( ALGBRC(153)*ALGBRC(153)+ ALGBRC(149)*ALGBRC(149)) ** (1.0 / 2) ALGBRC(162) = ALGBRC(37)/((1.00000+CONSTS(60)/CONSTS(49))+( ALGBRC(158)*(CONSTS(48)+CONSTS(60)))/( CONSTS(47)*CONSTS(48))) ALGBRC(166) = ( CONSTS(60)*ALGBRC(162))/CONSTS(49) ALGBRC(170) = ( ( CONSTS(60)*ALGBRC(162))*ALGBRC(158))/( CONSTS(47)*CONSTS(48)) RATES(26) = 0.00100000*( ( CONSTS(219)*CONSTS(63))*(ALGBRC(166)+ALGBRC(170)) - CONSTS(218)*STATES(26)) ALGBRC(177) = ALGBRC(58)/((1.00000+CONSTS(60)/CONSTS(54))+( ALGBRC(157)*(CONSTS(53)+CONSTS(60)))/( CONSTS(51)*CONSTS(53))) ALGBRC(181) = ( CONSTS(60)*ALGBRC(177))/CONSTS(54) ALGBRC(185) = ( ( CONSTS(60)*ALGBRC(177))*ALGBRC(157))/( CONSTS(51)*CONSTS(53)) RATES(14) = 0.00100000*( ( CONSTS(219)*CONSTS(59))*(ALGBRC(181)+ALGBRC(185)) - CONSTS(218)*STATES(14)) ALGBRC(178) = ALGBRC(61)/((1.00000+CONSTS(60)/CONSTS(54))+( ALGBRC(158)*(CONSTS(53)+CONSTS(60)))/( CONSTS(51)*CONSTS(53))) ALGBRC(182) = ( CONSTS(60)*ALGBRC(178))/CONSTS(54) ALGBRC(186) = ( ( CONSTS(60)*ALGBRC(178))*ALGBRC(158))/( CONSTS(51)*CONSTS(53)) RATES(28) = 0.00100000*( ( CONSTS(219)*CONSTS(63))*(ALGBRC(182)+ALGBRC(186)) - CONSTS(218)*STATES(28)) ALGBRC(189) = ALGBRC(169)+ CONSTS(61)*ALGBRC(185) ALGBRC(173) = ( ALGBRC(161)*ALGBRC(157))/CONSTS(47) ALGBRC(193) = ( ALGBRC(177)*ALGBRC(157))/CONSTS(51) ALGBRC(197) = ALGBRC(173)+ CONSTS(61)*ALGBRC(193) RATES(2) = 0.00100000*(( CONSTS(46)*ALGBRC(197)+ CONSTS(44)*ALGBRC(189)) - CONSTS(56)*STATES(2)) RATES(11) = 0.00100000*(( CONSTS(46)*ALGBRC(197)+ CONSTS(44)*ALGBRC(189)) - ( CONSTS(57)*STATES(11))*STATES(10)) ALGBRC(190) = ALGBRC(170)+ CONSTS(64)*ALGBRC(186) ALGBRC(174) = ( ALGBRC(162)*ALGBRC(158))/CONSTS(47) ALGBRC(194) = ( ALGBRC(178)*ALGBRC(158))/CONSTS(51) ALGBRC(198) = ALGBRC(174)+ CONSTS(64)*ALGBRC(194) RATES(4) = 0.00100000*(( CONSTS(46)*ALGBRC(198)+ CONSTS(44)*ALGBRC(190)) - CONSTS(56)*STATES(4)) RATES(25) = 0.00100000*(( CONSTS(46)*ALGBRC(198)+ CONSTS(44)*ALGBRC(190)) - ( CONSTS(57)*STATES(25))*STATES(24)) ALGBRC(134) = CONSTS(275) - STATES(36) ALGBRC(140) = CONSTS(263) - STATES(36)*CONSTS(262) ALGBRC(145) = - CONSTS(231)*STATES(36) ALGBRC(150) = - ALGBRC(134)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(134)*ALGBRC(134) - 3.00000*ALGBRC(140)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(134))*ALGBRC(140) - ( ( 2.00000*ALGBRC(134))*ALGBRC(134))*ALGBRC(134)) - 27.0000*ALGBRC(145))/( 2.00000* ALGBRC(134)*ALGBRC(134) - 3.00000*ALGBRC(140) ** 1.50000))/3.00000) ALGBRC(191) = ( CONSTS(80)*(1.00000 - STATES(43)))/(1.00000+CONSTS(82)/ALGBRC(150)) ALGBRC(195) = ALGBRC(191)+STATES(43) ALGBRC(200) = ALGBRC(195)/(ALGBRC(195)+CONSTS(81)) RATES(37) = (ALGBRC(200) - STATES(37))/CONSTS(232) RATES(38) = (ALGBRC(200) - STATES(38))/CONSTS(233) RATES(39) = (ALGBRC(200) - STATES(39))/CONSTS(234) RATES(40) = (ALGBRC(200) - STATES(40))/CONSTS(235) RATES(41) = (ALGBRC(200) - STATES(41))/CONSTS(86) ALGBRC(201) = 1.00000/(1.00000+CONSTS(81)/ALGBRC(195) ** 2.00000) RATES(42) = (ALGBRC(201) - STATES(42))/CONSTS(87) ALGBRC(199) = (CONSTS(206) - CONSTS(209))+STATES(145) ALGBRC(204) = 0.500000*( (ALGBRC(199) ** 2.00000+ ( 4.00000*CONSTS(206))*CONSTS(209)) ** (1.0 / 2) - ALGBRC(199)) RATES(126) = 0.00100000*(( ( CONSTS(163)*STATES(20))*(1.00000 - STATES(126)))/(CONSTS(159)+ (1.00000 - STATES(126))*1.00000) - ( ( CONSTS(164)*ALGBRC(204))*STATES(126))/(CONSTS(161)+ STATES(126)*1.00000)) ALGBRC(207) = (CONSTS(29)/CONSTS(320)+CONSTS(83)/CONSTS(327))+ALGBRC(204)/CONSTS(322) RATES(43) = ( CONSTS(84)*ALGBRC(191))*ALGBRC(195) - ( CONSTS(85)*STATES(43))*(0.100000+( 0.900000*ALGBRC(207))/0.137100) ALGBRC(156) = 1.00000/( (1.00000+EXP((4.79800+STATES(51))/- 7.56990))*(1.00000+EXP((25.0000+STATES(51))/- 5.00000))) ALGBRC(196) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(51)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(51)+13.0000))) ALGBRC(203) = ALGBRC(156)/ALGBRC(196) ALGBRC(206) = (1.00000 - ALGBRC(156))/ALGBRC(196) ALGBRC(180) = 1.00000+0.00200000/ALGBRC(150) ** 4.00000 ALGBRC(184) = 6.00000/ALGBRC(180) ALGBRC(164) = 1.00000/(1.00000+EXP((29.9790+STATES(51))/3.17750)) ALGBRC(172) = (0.100000+ALGBRC(164))/1.10000 ALGBRC(215) = 1.00000/( ( 70.0000*(1.00000 - 0.500000*STATES(37)))*(1.00000+EXP((STATES(51)+49.1000)/10.3490))) ALGBRC(217) = 1.00000/(1.00000+EXP((STATES(51)+0.213000)/- 10.8070)) ALGBRC(222) = 1.00000/(ALGBRC(215)+ALGBRC(217)/38.4940) ALGBRC(225) = ALGBRC(172)/ALGBRC(222) ALGBRC(234) = (1.00000 - ALGBRC(172))/ALGBRC(222) RATES(52) = (( - ((ALGBRC(203)+ALGBRC(184))+ALGBRC(234))*STATES(52)+ ALGBRC(206)*STATES(56))+ CONSTS(103)*STATES(55))+ ALGBRC(225)*STATES(53) RATES(56) = (( - ((ALGBRC(206)+ALGBRC(184))+ALGBRC(234))*STATES(56)+ ALGBRC(203)*STATES(52))+ CONSTS(103)*STATES(59))+ ALGBRC(225)*STATES(57) ALGBRC(168) = (0.000100000+ALGBRC(164))/1.00010 ALGBRC(160) = 0.100000*STATES(37) ALGBRC(176) = 1.00000+0.0100000/ALGBRC(150) ** 10.0000 ALGBRC(188) = (32.5000 - (18.0000 - ALGBRC(160))/ALGBRC(180)) - 10.0000/ALGBRC(176) ALGBRC(219) = 1.00000/(ALGBRC(215)+ALGBRC(217)/ALGBRC(188)) ALGBRC(228) = ALGBRC(168)/ALGBRC(219) ALGBRC(240) = (1.00000 - ALGBRC(168))/ALGBRC(219) RATES(55) = (( - ((ALGBRC(203)+CONSTS(103))+ALGBRC(240))*STATES(55)+ ALGBRC(184)*STATES(52))+ ALGBRC(206)*STATES(59))+ ALGBRC(228)*STATES(54) RATES(59) = (( - ((ALGBRC(206)+CONSTS(103))+ALGBRC(240))*STATES(59)+ ALGBRC(184)*STATES(56))+ ALGBRC(203)*STATES(55))+ ALGBRC(228)*STATES(58) ALGBRC(231) = TERNRY(ABS(ALGBRC(228)).LT.1.00000e-12, 1.00000e-12, ALGBRC(228)) ALGBRC(237) = TERNRY(ABS(ALGBRC(234)).LT.1.00000e-12, 1.00000e-12, ALGBRC(234)) ALGBRC(243) = ( CONSTS(104)*( ( ALGBRC(225)*ALGBRC(240))*ALGBRC(184)))/( ( ALGBRC(237)*ALGBRC(231))*CONSTS(103)) RATES(53) = (( - ((ALGBRC(203)+ALGBRC(243))+ALGBRC(225))*STATES(53)+ ALGBRC(234)*STATES(52))+ CONSTS(104)*STATES(54))+ ALGBRC(206)*STATES(57) RATES(54) = (( - ((ALGBRC(203)+CONSTS(104))+ALGBRC(228))*STATES(54)+ ALGBRC(240)*STATES(55))+ ALGBRC(243)*STATES(53))+ ALGBRC(206)*STATES(58) RATES(57) = (( - ((ALGBRC(206)+ALGBRC(243))+ALGBRC(225))*STATES(57)+ ALGBRC(234)*STATES(56))+ CONSTS(104)*STATES(58))+ ALGBRC(203)*STATES(53) RATES(58) = (( - ((ALGBRC(206)+CONSTS(104))+ALGBRC(228))*STATES(58)+ ALGBRC(240)*STATES(59))+ ALGBRC(243)*STATES(57))+ ALGBRC(203)*STATES(54) ALGBRC(63) = CONSTS(275) - STATES(34) ALGBRC(75) = CONSTS(263) - STATES(34)*CONSTS(262) ALGBRC(91) = - CONSTS(231)*STATES(34) ALGBRC(103) = - ALGBRC(63)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(63)*ALGBRC(63) - 3.00000*ALGBRC(75)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(63))*ALGBRC(75) - ( ( 2.00000*ALGBRC(63))*ALGBRC(63))*ALGBRC(63)) - 27.0000*ALGBRC(91))/( 2.00000* ALGBRC(63)*ALGBRC(63) - 3.00000*ALGBRC(75) ** 1.50000))/3.00000) ALGBRC(210) = (ALGBRC(150) - ALGBRC(103))/CONSTS(96) ALGBRC(11) = (STATES(5)+CONSTS(393))/CONSTS(18) ALGBRC(221) = (ALGBRC(11) - CONSTS(398))/(0.927300 - CONSTS(398)) ALGBRC(224) = TERNRY(ALGBRC(221).LT.0.00000, 0.00000, ALGBRC(221)) ALGBRC(227) = 0.000257900*(1.00000+ 0.100000*STATES(37)) ALGBRC(230) = EXP( ( 2.00000*STATES(51))*CONSTS(264)) ALGBRC(233) = ( ( ( ( ALGBRC(227)*4.00000)*STATES(51))*CONSTS(276))*( ALGBRC(150)*ALGBRC(230) - 0.341000*CONSTS(98)))/(ALGBRC(230) - 1.00000) ALGBRC(236) = ALGBRC(233)*(STATES(56)+STATES(59)) ALGBRC(239) = 0.000155200*(1.00000+ 0.400000*STATES(37)) ALGBRC(242) = EXP( ( 2.00000*STATES(51))*CONSTS(264)) ALGBRC(245) = ( ( ( ( ALGBRC(239)*4.00000)*STATES(51))*CONSTS(276))*( ALGBRC(150)*ALGBRC(242) - 0.341000*CONSTS(98)))/(ALGBRC(242) - 1.00000) ALGBRC(247) = ALGBRC(245)*(STATES(64)+STATES(67)) ALGBRC(249) = (1.00000 - ALGBRC(224))*ALGBRC(247)+ ALGBRC(224)*ALGBRC(236) RATES(34) = - CONSTS(306)*ALGBRC(249)+ CONSTS(316)*ALGBRC(210) ALGBRC(192) = 1.00000/( (1.00000+EXP((13.5600 - STATES(51))/9.45000))*(1.00000+EXP((25.0000+STATES(51))/- 5.00000))) ALGBRC(209) = ALGBRC(192)/ALGBRC(196) ALGBRC(211) = (1.00000 - ALGBRC(192))/ALGBRC(196) ALGBRC(235) = 1.00000+0.00110000/ALGBRC(150) ** 4.00000 ALGBRC(238) = 14.9186/ALGBRC(235) ALGBRC(220) = 1.00000/(1.00000+EXP((17.5000+STATES(51))/3.00000)) ALGBRC(226) = (0.247400+ALGBRC(220))/1.24740 ALGBRC(229) = 1.00000/(ALGBRC(215)+ALGBRC(217)/26.5530) ALGBRC(246) = ALGBRC(226)/ALGBRC(229) ALGBRC(252) = (1.00000 - ALGBRC(226))/ALGBRC(229) RATES(60) = (( - ((ALGBRC(209)+ALGBRC(238))+ALGBRC(252))*STATES(60)+ ALGBRC(211)*STATES(64))+ CONSTS(105)*STATES(63))+ ALGBRC(246)*STATES(61) RATES(64) = (( - ((ALGBRC(211)+ALGBRC(238))+ALGBRC(252))*STATES(64)+ ALGBRC(209)*STATES(60))+ CONSTS(105)*STATES(67))+ ALGBRC(246)*STATES(65) ALGBRC(223) = (0.00100000+ALGBRC(220))/1.00100 ALGBRC(213) = 5.00000*STATES(37) ALGBRC(232) = 1.00000+0.0120000/ALGBRC(150) ** 10.0000 ALGBRC(241) = (13.8250 - (6.38360 - ALGBRC(213))/ALGBRC(235)) - 3.36960/ALGBRC(232) ALGBRC(244) = 1.00000/(ALGBRC(215)+ALGBRC(217)/ALGBRC(241)) ALGBRC(248) = ALGBRC(223)/ALGBRC(244) ALGBRC(258) = (1.00000 - ALGBRC(223))/ALGBRC(244) RATES(63) = (( - ((ALGBRC(209)+CONSTS(105))+ALGBRC(258))*STATES(63)+ ALGBRC(238)*STATES(60))+ ALGBRC(211)*STATES(67))+ ALGBRC(248)*STATES(62) RATES(67) = (( - ((ALGBRC(211)+CONSTS(105))+ALGBRC(258))*STATES(67)+ ALGBRC(238)*STATES(64))+ ALGBRC(209)*STATES(63))+ ALGBRC(248)*STATES(66) ALGBRC(250) = TERNRY(ABS(ALGBRC(248)).LT.1.00000e-12, 1.00000e-12, ALGBRC(248)) ALGBRC(255) = TERNRY(ABS(ALGBRC(252)).LT.1.00000e-12, 1.00000e-12, ALGBRC(252)) ALGBRC(261) = ( CONSTS(106)*( ( ALGBRC(246)*ALGBRC(258))*ALGBRC(238)))/( ( ALGBRC(255)*ALGBRC(250))*CONSTS(105)) RATES(61) = (( - ((ALGBRC(209)+ALGBRC(261))+ALGBRC(246))*STATES(61)+ ALGBRC(252)*STATES(60))+ CONSTS(106)*STATES(62))+ ALGBRC(211)*STATES(65) RATES(62) = (( - ((ALGBRC(209)+CONSTS(106))+ALGBRC(248))*STATES(62)+ ALGBRC(258)*STATES(63))+ ALGBRC(261)*STATES(61))+ ALGBRC(211)*STATES(66) RATES(65) = (( - ((ALGBRC(211)+ALGBRC(261))+ALGBRC(246))*STATES(65)+ ALGBRC(252)*STATES(64))+ CONSTS(106)*STATES(66))+ ALGBRC(209)*STATES(61) RATES(66) = (( - ((ALGBRC(211)+CONSTS(106))+ALGBRC(248))*STATES(66)+ ALGBRC(258)*STATES(67))+ ALGBRC(261)*STATES(65))+ ALGBRC(209)*STATES(62) ALGBRC(253) = CONSTS(244)*(1.00000+ 2.00000*STATES(42)) ALGBRC(256) = 0.112500*ALGBRC(253) ALGBRC(114) = (CONSTS(73)+CONSTS(74)) - STATES(35) ALGBRC(122) = STATES(35)*CONSTS(74) ALGBRC(128) = ( ( ALGBRC(114)*ALGBRC(114)+ 4.00000*ALGBRC(122)) ** (1.0 / 2) - ALGBRC(114))/2.00000 ALGBRC(270) = ( ALGBRC(249)*1.00000)/(1.00000+1.00000/ALGBRC(128) ** 8.00000) ALGBRC(272) = ALGBRC(256)*ALGBRC(270) ALGBRC(264) = 1.00000+0.0123000/ALGBRC(128) ALGBRC(266) = ALGBRC(253)/ALGBRC(264) RATES(119) = - (ALGBRC(272)+STATES(119))/ALGBRC(266) ALGBRC(259) = CONSTS(244)*(1.00000+ 0.00000*STATES(42)) ALGBRC(262) = 0.112500*ALGBRC(259) ALGBRC(273) = ( 1.99250*ALGBRC(262))*ALGBRC(270) ALGBRC(268) = ( 0.535700*ALGBRC(259))/ALGBRC(264) RATES(120) = - (ALGBRC(273)+STATES(120))/ALGBRC(268) ALGBRC(287) = 1.00000+CONSTS(135)/ALGBRC(150) ** 2.00000 ALGBRC(286) = EXP( ( (CONSTS(140) - 1.00000)*STATES(51))*CONSTS(264)) ALGBRC(288) = 1.00000+ CONSTS(141)*ALGBRC(286) ALGBRC(284) = STATES(50) ** 3.00000 ALGBRC(289) = ( CONSTS(137)*ALGBRC(284)+ CONSTS(365)*ALGBRC(150))+ ( CONSTS(363)*CONSTS(98))*(1.00000+ALGBRC(150)/CONSTS(136)) ALGBRC(290) = ( ( CONSTS(136)*CONSTS(367))*(1.00000+ALGBRC(284)/CONSTS(363))+ ALGBRC(284)*CONSTS(98))+ CONSTS(367)*ALGBRC(150) ALGBRC(285) = EXP( ( CONSTS(140)*STATES(51))*CONSTS(264)) ALGBRC(291) = ( 0.200000*CONSTS(142))*( ( ALGBRC(284)*CONSTS(98))*ALGBRC(285) - ( CONSTS(367)*ALGBRC(150))*ALGBRC(286)) ALGBRC(292) = ALGBRC(291)/( ( ALGBRC(287)*ALGBRC(288))*(ALGBRC(289)+ALGBRC(290))) ALGBRC(208) = (STATES(50) - STATES(49))/CONSTS(95) RATES(50) = - ( CONSTS(410)*ALGBRC(292)+ALGBRC(208)) ALGBRC(39) = (STATES(6)+CONSTS(396))/CONSTS(31) ALGBRC(312) = (ALGBRC(39) - CONSTS(399))/(0.958600 - CONSTS(399)) ALGBRC(314) = TERNRY(ALGBRC(312).LT.0.00000, 0.00000, ALGBRC(312)) ALGBRC(318) = (1.00000 - ALGBRC(314))*STATES(119)+ ALGBRC(314)*STATES(120) ALGBRC(308) = ( CONSTS(156)*EXP(ALGBRC(128)/CONSTS(154)))*(ALGBRC(128) - ALGBRC(150)) ALGBRC(310) = ( CONSTS(157)*EXP(ALGBRC(128)/CONSTS(155)))*(ALGBRC(128) - ALGBRC(150)) ALGBRC(316) = (1.00000 - ALGBRC(314))*ALGBRC(308)+ ALGBRC(314)*ALGBRC(310) ALGBRC(320) = ALGBRC(316)+ALGBRC(318) ALGBRC(212) = (STATES(31) - ALGBRC(128))/CONSTS(97) RATES(35) = ALGBRC(212) - ALGBRC(320) ALGBRC(154) = (STATES(32) - 0.673519)/(0.999180 - 0.673519) ALGBRC(159) = TERNRY(ALGBRC(154).LT.0.00000, 0.00000, ALGBRC(154)) ALGBRC(163) = (1.00000 - ALGBRC(159))*CONSTS(78)+ ALGBRC(159)*CONSTS(230) ALGBRC(175) = ALGBRC(163)+CONSTS(76) ALGBRC(179) = ((ALGBRC(175) - STATES(33))+CONSTS(72))+CONSTS(79) ALGBRC(167) = ALGBRC(163)*CONSTS(76) ALGBRC(183) = ((ALGBRC(167) - STATES(33)*ALGBRC(175))+ CONSTS(79)*CONSTS(76))+ CONSTS(72)*ALGBRC(163) ALGBRC(171) = - ALGBRC(167)*STATES(33) ALGBRC(187) = - ALGBRC(179)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(179)*ALGBRC(179) - 3.00000*ALGBRC(183)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(179))*ALGBRC(183) - ( ( 2.00000*ALGBRC(179))*ALGBRC(179))*ALGBRC(179)) - 27.0000*ALGBRC(171))/( 2.00000* ALGBRC(179)*ALGBRC(179) - 3.00000*ALGBRC(183) ** 1.50000))/3.00000) ALGBRC(202) = (ALGBRC(150) - ALGBRC(187))/CONSTS(95) RATES(36) = - ((( CONSTS(317)*ALGBRC(292)+ CONSTS(318)*ALGBRC(320))+ALGBRC(202))+ALGBRC(210)) ALGBRC(251) = EXP( STATES(51)*CONSTS(264)) ALGBRC(254) = ( ( ( CONSTS(108)*STATES(51))*CONSTS(276))*(STATES(47) - CONSTS(99)*ALGBRC(251)))/(1.00000 - ALGBRC(251)) ALGBRC(321) = 1.00000 - 1.00000/(1.00000+ALGBRC(318)/CONSTS(109) ** 2.00000) ALGBRC(324) = ( ALGBRC(254)*ALGBRC(321))*STATES(68) ALGBRC(205) = (STATES(48) - STATES(47))/CONSTS(95) RATES(48) = CONSTS(342)*ALGBRC(324) - ALGBRC(205) ALGBRC(307) = 1.00000/(1.00000+0.0300000/ALGBRC(195) ** 2.00000) ALGBRC(309) = (1.00000 - ALGBRC(307))*CONSTS(162)+ ALGBRC(307)*CONSTS(380) ALGBRC(325) = ( ALGBRC(309)*STATES(31))/CONSTS(165) ALGBRC(311) = (STATES(126) - 0.659100)/(0.994500 - 0.659100) ALGBRC(313) = TERNRY(ALGBRC(311).LT.0.00000, 0.00000, ALGBRC(311)) ALGBRC(315) = ALGBRC(313)*STATES(41) ALGBRC(317) = STATES(41) - ALGBRC(315) ALGBRC(319) = ALGBRC(313) - ALGBRC(315) ALGBRC(322) = ((1.00000 - ALGBRC(319)) - ALGBRC(317)) - ALGBRC(315) ALGBRC(323) = (( ALGBRC(322)*CONSTS(160)+ ALGBRC(319)*CONSTS(375))+ ALGBRC(317)*CONSTS(370))+ ALGBRC(315)*CONSTS(377) ALGBRC(327) = ( ALGBRC(309)*ALGBRC(187))/(ALGBRC(187)+ALGBRC(323)) ALGBRC(329) = ALGBRC(327) - ALGBRC(325) RATES(31) = ALGBRC(329) - CONSTS(313)*ALGBRC(212) ALGBRC(214) = ( 2.00000*STATES(51))*CONSTS(264) ALGBRC(216) = EXP(ALGBRC(214)) ALGBRC(218) = ( ( ( ( CONSTS(102)*2.00000)*CONSTS(88))*ALGBRC(214))*( ALGBRC(187)*ALGBRC(216) - 0.341000*CONSTS(98)))/(ALGBRC(216) - 1.00000) ALGBRC(293) = 1.00000+CONSTS(135)/ALGBRC(187) ** 2.00000 ALGBRC(294) = 1.00000+ CONSTS(141)*ALGBRC(286) ALGBRC(283) = STATES(49) ** 3.00000 ALGBRC(295) = ( CONSTS(137)*ALGBRC(283)+ CONSTS(365)*ALGBRC(187))+ ( CONSTS(363)*CONSTS(98))*(1.00000+ALGBRC(187)/CONSTS(136)) ALGBRC(296) = ( ( CONSTS(136)*CONSTS(367))*(1.00000+ALGBRC(283)/CONSTS(363))+ ALGBRC(283)*CONSTS(98))+ CONSTS(367)*ALGBRC(187) ALGBRC(297) = ( 0.800000*CONSTS(142))*( ( ALGBRC(283)*CONSTS(98))*ALGBRC(285) - ( CONSTS(367)*ALGBRC(187))*ALGBRC(286)) ALGBRC(298) = ALGBRC(297)/( ( ALGBRC(293)*ALGBRC(294))*(ALGBRC(295)+ALGBRC(296))) ALGBRC(305) = ( CONSTS(152)*ALGBRC(187))/(CONSTS(153)+ALGBRC(187)) RATES(33) = ( - CONSTS(309)*((ALGBRC(218)+ALGBRC(305)) - 2.00000*ALGBRC(298)) - CONSTS(314)*ALGBRC(329))+ CONSTS(319)*ALGBRC(202) ALGBRC(328) = - CONSTS(239)*log(CONSTS(99)/STATES(47)) ALGBRC(332) = CONSTS(239)*log(CONSTS(100)/STATES(127)) ALGBRC(333) = ALGBRC(332) - ALGBRC(328) ALGBRC(334) = ( CONSTS(92)*ALGBRC(333))/(ALGBRC(333)+CONSTS(93)) ALGBRC(353) = CONSTS(239)*log(CONSTS(101)/STATES(49)) ALGBRC(354) = ALGBRC(353) - ALGBRC(328) ** 4.00000 ALGBRC(355) = ( CONSTS(94)*ALGBRC(354))/(ALGBRC(354)+CONSTS(347)) ALGBRC(330) = CONSTS(107)*(STATES(51) - ALGBRC(328)) RATES(47) = (( CONSTS(333)*ALGBRC(330)+ALGBRC(355))+ALGBRC(334))+ CONSTS(338)*ALGBRC(205) ALGBRC(356) = ( ( ( CONSTS(242)*STATES(115) ** 3.00000)*STATES(107))*STATES(108))*(STATES(51) - ALGBRC(353)) ALGBRC(357) = ( ( ( CONSTS(242)*STATES(109) ** 3.00000)*STATES(107))*STATES(108))*(STATES(51) - ALGBRC(353)) ALGBRC(358) = ( ( ( CONSTS(242)*STATES(112) ** 3.00000)*STATES(110))*STATES(111))*(STATES(51) - ALGBRC(353)) ALGBRC(359) = ( ( ( ( CONSTS(242)*STATES(115) ** 3.00000)*STATES(113))*STATES(114))*(STATES(51) - ALGBRC(353)))*1.25000 ALGBRC(274) = (STATES(106) - 0.239480)/(0.950143 - 0.239480) ALGBRC(275) = TERNRY(ALGBRC(274).LT.0.00000, 0.00000, ALGBRC(274)) ALGBRC(276) = ALGBRC(275)*STATES(39) ALGBRC(277) = STATES(39) - ALGBRC(276) ALGBRC(278) = ALGBRC(275) - ALGBRC(276) ALGBRC(279) = ((1.00000 - ALGBRC(278)) - ALGBRC(277)) - ALGBRC(276) ALGBRC(360) = (( ALGBRC(279)*ALGBRC(358)+ ALGBRC(278)*ALGBRC(359))+ ALGBRC(277)*ALGBRC(357))+ ALGBRC(276)*ALGBRC(356) ALGBRC(280) = STATES(51)*CONSTS(264) ALGBRC(281) = EXP(ALGBRC(280)) ALGBRC(282) = ( ( ( CONSTS(134)*CONSTS(88))*ALGBRC(280))*( STATES(49)*ALGBRC(281) - CONSTS(101)))/(ALGBRC(281) - 1.00000) ALGBRC(301) = ( CONSTS(145)*CONSTS(243))/(1.00000+EXP( - (STATES(51)+92.0000)*CONSTS(264))) ALGBRC(302) = ALGBRC(301)*STATES(49)/(STATES(49)+CONSTS(148)) ** 3.00000 ALGBRC(303) = ALGBRC(301)*STATES(49)/(STATES(49)+CONSTS(149)) ** 3.00000 ALGBRC(299) = (STATES(116) - 0.126345)/(0.998014 - 0.126345) ALGBRC(300) = TERNRY(ALGBRC(299).LT.0.00000, 0.00000, ALGBRC(299)) ALGBRC(304) = (1.00000 - ALGBRC(300))*ALGBRC(302)+ ALGBRC(300)*ALGBRC(303) ALGBRC(361) = ( STATES(118) ** 3.00000*STATES(117))*(STATES(51) - ALGBRC(353)) ALGBRC(362) = 0.0160000*ALGBRC(361) ALGBRC(363) = 0.00650000*ALGBRC(361) ALGBRC(364) = (1.00000 - STATES(39))*ALGBRC(363)+ STATES(39)*ALGBRC(362) ALGBRC(366) = (((ALGBRC(360)+ALGBRC(282))+ALGBRC(364))+ ALGBRC(304)*3.00000)+ ALGBRC(298)*3.00000 RATES(49) = ( CONSTS(408)*ALGBRC(366)+ CONSTS(409)*ALGBRC(208))+ALGBRC(355) ALGBRC(306) = (((ALGBRC(249)+ALGBRC(218))+ALGBRC(305)) - ALGBRC(298)*2.00000) - ALGBRC(292)*2.00000 ALGBRC(331) = ALGBRC(330)+ALGBRC(324) ALGBRC(335) = STATES(51) - ALGBRC(332) ALGBRC(336) = 1.02000/(1.00000+EXP( 0.238500*(ALGBRC(335) - 59.2150))) ALGBRC(337) = ( 0.491240*EXP( 0.0803200*(ALGBRC(335)+5.47600))+EXP( 0.0617500*(ALGBRC(335) - 594.310)))/(1.00000+EXP( - 0.514300*(ALGBRC(335)+4.75300))) ALGBRC(338) = ( CONSTS(240)*(ALGBRC(336)/(ALGBRC(336)+ALGBRC(337))))*ALGBRC(335) ALGBRC(339) = ALGBRC(338)*1.20000 ALGBRC(340) = (1.00000 - STATES(38))*ALGBRC(338)+ STATES(38)*ALGBRC(339) ALGBRC(257) = 1.00000/(1.00000+EXP((STATES(51)+10.0000)/15.4000)) ALGBRC(341) = ( ( CONSTS(241)*STATES(69))*ALGBRC(257))*(STATES(51) - ALGBRC(332)) ALGBRC(349) = CONSTS(239)*log((CONSTS(100)+ CONSTS(166)*CONSTS(101))/(STATES(127)+ CONSTS(166)*STATES(49))) ALGBRC(260) = 0.195610*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(187) ** 1.40000)) ALGBRC(350) = ( ALGBRC(260)*(STATES(70)+STATES(71)))*(STATES(51) - ALGBRC(349)) ALGBRC(351) = ( ALGBRC(260)*(STATES(72)+STATES(73)))*(STATES(51) - ALGBRC(349)) ALGBRC(263) = (STATES(104)+CONSTS(404))/CONSTS(111) ALGBRC(265) = (ALGBRC(263) - CONSTS(405))/(0.785000 - CONSTS(405)) ALGBRC(267) = TERNRY(ALGBRC(265).LT.0.00000, 0.00000, ALGBRC(265)) ALGBRC(352) = ALGBRC(267)*ALGBRC(351)+ (1.00000 - ALGBRC(267))*ALGBRC(350) ALGBRC(342) = ( CONSTS(124)*(STATES(51) - ALGBRC(332)))/(1.00000+EXP((15.0000 - STATES(51))/17.0000)) ALGBRC(343) = (( CONSTS(124)*(STATES(51) - ALGBRC(332)))/(1.00000+EXP((36.0000 - STATES(51))/17.0000)))*3.62000 ALGBRC(269) = (STATES(105) - 0.0589380)/(0.393747 - 0.0589380) ALGBRC(271) = TERNRY(ALGBRC(269).LT.0.00000, 0.00000, ALGBRC(269)) ALGBRC(344) = (1.00000 - ALGBRC(271))*ALGBRC(342)+ ALGBRC(271)*ALGBRC(343) ALGBRC(326) = EXP(STATES(51)/550.000) ALGBRC(345) = ( ( CONSTS(158)*STATES(121) ** 3.00000)*ALGBRC(326))*(STATES(51) - ALGBRC(332)) ALGBRC(346) = ALGBRC(345)*( 0.735600*STATES(122)+ 0.264400*STATES(124)) ALGBRC(347) = ALGBRC(345)*( 0.735600*STATES(123)+ 0.264400*STATES(125)) ALGBRC(348) = (1.00000 - STATES(40))*ALGBRC(347)+ STATES(40)*ALGBRC(346) ALGBRC(365) = ((((ALGBRC(340)+ALGBRC(341))+ALGBRC(352))+ALGBRC(344))+ALGBRC(348)) - 2.00000*ALGBRC(304) ALGBRC(367) = ALGBRC(366)+ ALGBRC(292)*3.00000 ALGBRC(368) = ((ALGBRC(367)+ALGBRC(365))+ALGBRC(306))+ALGBRC(331) ALGBRC(369) = (TERNRY((VOI - CONSTS(214)) - CONSTS(215)*INT((VOI - CONSTS(214))/CONSTS(215)).LT.CONSTS(213), 1.00000, 0.00000))*CONSTS(257) RATES(51) = - (ALGBRC(368)+ALGBRC(369))/1.00000 RATES(127) = CONSTS(407)*(ALGBRC(365)+ALGBRC(369))+ALGBRC(334) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(4) = CONSTS(394) - STATES(5) ALGBRC(5) = CONSTS(397) - STATES(6) ALGBRC(16) = CONSTS(406) - STATES(104) ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(51)+91.0000)/6.10000)) ALGBRC(30) = CONSTS(256) - STATES(145) ALGBRC(34) = (CONSTS(323) - STATES(13)) - STATES(12) ALGBRC(36) = (CONSTS(326) - STATES(19)) - STATES(18) ALGBRC(37) = (CONSTS(331) - STATES(27)) - STATES(26) ALGBRC(12) = 0.0250000/(1.00000+EXP((STATES(51)+58.0000)/5.00000)) ALGBRC(40) = 0.200000/(1.00000+EXP((STATES(51)+19.0000)/- 9.00000)) ALGBRC(13) = 1.00000/(( 0.000600000*(STATES(51) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(51) - 1.73840))) - ( 0.000300000*(STATES(51)+38.3608))/(1.00000 - EXP( 0.152200*(STATES(51)+38.3608)))) ALGBRC(41) = 1.00000/(1.00000+EXP((STATES(51)+10.0850)/- 4.25000)) ALGBRC(14) = 0.00739900/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0311960))/- 0.800190)) ALGBRC(42) = 0.00569920/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0415200))/1.34890)) ALGBRC(15) = 0.00994150/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0448090))/- 0.581720)) ALGBRC(43) = 0.00332010/(1.00000+EXP(( CONSTS(264)*(STATES(51) - 0.0942170))/0.953640)) ALGBRC(17) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 0.00000, 0.135000*EXP(((87.0000+STATES(51))+CONSTS(131))/- 6.80000)) ALGBRC(44) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP(((STATES(51)+CONSTS(131))+27.4034)/- 11.1000))), 3.56000*EXP( 0.0790000*((STATES(51)+CONSTS(131))+7.00000))+ 310000.*EXP( 0.350000*((STATES(51)+CONSTS(131))+7.00000))) ALGBRC(18) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*(STATES(51)+CONSTS(131)))+ - 6.94800e-05*EXP( - 0.0439100*(STATES(51)+CONSTS(131))))*((STATES(51)+CONSTS(131))+37.7800))/(1.00000+EXP( 0.311000*((STATES(51)+CONSTS(131))+79.2300)))) ALGBRC(45) = TERNRY(STATES(51)+CONSTS(131).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*(STATES(51)+CONSTS(131))))/(1.00000+EXP( - 0.100000*((STATES(51)+CONSTS(131))+32.0000))), ( 0.121200*EXP( - 0.0105200*(STATES(51)+CONSTS(131))))/(1.00000+EXP( - 0.137800*((STATES(51)+CONSTS(131))+40.1400)))) ALGBRC(19) = ( 0.320000*(STATES(51)+58.4729))/(1.00000 - EXP( - 0.100000*(STATES(51)+58.4729))) ALGBRC(46) = 0.0800000*EXP((13.7299 - STATES(51))/11.0000) ALGBRC(20) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, 0.135000*EXP((87.0000+STATES(51))/- 6.80000)) ALGBRC(47) = TERNRY(STATES(51).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(51)+27.4034)/- 11.1000))), 3.56000*EXP( 0.0790000*(STATES(51)+7.00000))+ 310000.*EXP( 0.350000*(STATES(51)+7.00000))) ALGBRC(21) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(51)) - 6.94800e-05*EXP( - 0.0439100*STATES(51)))*(STATES(51)+37.7800))/(1.00000+EXP( 0.311000*(STATES(51)+79.2300)))) ALGBRC(48) = TERNRY(STATES(51).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(51)))/(1.00000+EXP( - 0.100000*(STATES(51)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(51)))/(1.00000+EXP( - 0.137800*(STATES(51)+40.1400)))) ALGBRC(22) = ( 0.320000*(STATES(51)+58.4729))/(1.00000 - EXP( - 0.100000*(STATES(51)+58.4729))) ALGBRC(49) = 0.0800000*EXP((13.7299 - STATES(51))/11.0000) ALGBRC(23) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, 0.135000*EXP(((87.0000+STATES(51))+CONSTS(133))/- 6.80000)) ALGBRC(50) = TERNRY(STATES(51).GE.- 40.0000, (1.00000/0.130000)/(1.00000+EXP(((STATES(51)+CONSTS(133))+27.4034)/- 11.1000)), 3.56000*EXP( 0.0790000*((STATES(51)+CONSTS(133))+7.00000))+ 310000.*EXP( 0.350000*((STATES(51)+CONSTS(133))+7.00000))) ALGBRC(24) = TERNRY(STATES(51).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*(STATES(51)+CONSTS(133))) - 6.94800e-05*EXP( - 0.0439100*(STATES(51)+CONSTS(133))))*((STATES(51)+CONSTS(133))+37.7800))/(1.00000+EXP( 0.311000*((STATES(51)+CONSTS(133))+79.2300)))) ALGBRC(51) = TERNRY(STATES(51).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*(STATES(51)+CONSTS(133))))/(1.00000+EXP( - 0.100000*((STATES(51)+CONSTS(133))+32.0000))), ( 0.121200*EXP( - 0.0105200*(STATES(51)+CONSTS(133))))/(1.00000+EXP( - 0.137800*((STATES(51)+CONSTS(133))+40.1400)))) ALGBRC(25) = ( 0.320000*((STATES(51)+CONSTS(132))+58.4729))/(1.00000 - EXP( - 0.100000*((STATES(51)+CONSTS(132))+58.4729))) ALGBRC(52) = 0.0800000*EXP(((STATES(51)+CONSTS(132)) - 13.7299)/- 11.0000) ALGBRC(27) = ( 0.320000*(STATES(51)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(51)+47.1300))) ALGBRC(53) = 0.0800000*EXP(STATES(51)/- 11.0000) ALGBRC(58) = (CONSTS(324) - STATES(15)) - STATES(14) ALGBRC(61) = (CONSTS(332) - STATES(29)) - STATES(28) ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(51)+9.43700)/- 7.13300)) ALGBRC(54) = 1.00000/(1.00000+EXP((STATES(51) - 18.4099)/- 29.3814)) ALGBRC(66) = 1.00000/(1.00000+EXP((STATES(51)+100.000)/29.3814)) ALGBRC(78) = 1.00000/(ALGBRC(54)/1.20890+ 3.50000*ALGBRC(66)) ALGBRC(67) = ((1.00000/(1.00000+EXP((STATES(51)+19.0000)/- 9.00000)))/0.500000)/9.79530 ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(51)+58.0000)/5.00000)) ALGBRC(79) = 0.0479600*ALGBRC(29) ALGBRC(80) = 0.0214400*ALGBRC(29) ALGBRC(55) = (1.00000/(1.00000+EXP((STATES(51)+60.0000)/5.00000)))/250.000 ALGBRC(81) = 2.46000*ALGBRC(55) ALGBRC(82) = 0.560340*ALGBRC(55) ALGBRC(64) = 0.0906540*EXP( ( - 0.111570*STATES(51))*CONSTS(264)) ALGBRC(92) = 0.388390/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.150190))/- 0.606930)) ALGBRC(76) = 0.00311240+(0.0283300 - 0.00311240)/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.0516600))/1.55220)) ALGBRC(65) = 0.0657000*EXP( ( - 0.118990*STATES(51))*CONSTS(264)) ALGBRC(94) = 0.563560/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.179860))/- 0.583810)) ALGBRC(77) = 0.000385250+(0.0124060 - 0.000385250)/(1.00000+EXP(( CONSTS(264)*(STATES(51)+0.0641180))/0.779920)) ALGBRC(93) = 0.000441980*EXP( ( - 1.20220*STATES(51))*CONSTS(264)) ALGBRC(104) = 0.000401730*EXP( ( 0.000208730*STATES(51))*CONSTS(264)) ALGBRC(95) = 0.000237300*EXP( ( - 1.97420*STATES(51))*CONSTS(264)) ALGBRC(105) = 0.000226520*EXP( ( 0.000246890*STATES(51))*CONSTS(264)) ALGBRC(6) = ( ( CONSTS(37)*CONSTS(36))*CONSTS(320) - STATES(9)) - STATES(8) ALGBRC(35) = ( ALGBRC(6)*(CONSTS(60)+CONSTS(52)) - STATES(15)*(CONSTS(52)+CONSTS(60)))+ ( CONSTS(50)*CONSTS(52))*(1.00000+CONSTS(60)/CONSTS(55)) ALGBRC(59) = ( - STATES(15)*CONSTS(50))*CONSTS(52) ALGBRC(71) = (- ALGBRC(35)+ ( ALGBRC(35)*ALGBRC(35) - ( 4.00000*CONSTS(223))*ALGBRC(59)) ** (1.0 / 2))/( 2.00000*CONSTS(223)) ALGBRC(87) = ALGBRC(6)/(1.00000+ (ALGBRC(71)/CONSTS(50))*(1.00000+CONSTS(60)/CONSTS(52))) ALGBRC(99) = ( ALGBRC(71)*ALGBRC(87))/CONSTS(50) ALGBRC(109) = ( ALGBRC(99)*CONSTS(60))/CONSTS(52) ALGBRC(9) = ( ( CONSTS(216)*CONSTS(36))*CONSTS(327) - STATES(22)) - STATES(21) ALGBRC(38) = ( ALGBRC(9)*(CONSTS(60)+CONSTS(52)) - STATES(29)*(CONSTS(52)+CONSTS(60)))+ ( CONSTS(50)*CONSTS(52))*(1.00000+CONSTS(60)/CONSTS(55)) ALGBRC(62) = ( - STATES(29)*CONSTS(50))*CONSTS(52) ALGBRC(74) = (- ALGBRC(38)+ ( ALGBRC(38)*ALGBRC(38) - ( 4.00000*CONSTS(226))*ALGBRC(62)) ** (1.0 / 2))/( 2.00000*CONSTS(226)) ALGBRC(90) = ALGBRC(9)/(1.00000+ (ALGBRC(74)/CONSTS(50))*(1.00000+CONSTS(60)/CONSTS(52))) ALGBRC(102) = ( ALGBRC(74)*ALGBRC(90))/CONSTS(50) ALGBRC(112) = (CONSTS(60)/CONSTS(52))*ALGBRC(102) ALGBRC(8) = ( ( CONSTS(217)*CONSTS(296))*CONSTS(322) - STATES(3)) - STATES(16) ALGBRC(60) = ( ALGBRC(8)*(CONSTS(48)+CONSTS(60)) - ALGBRC(36)*(CONSTS(48)+CONSTS(60)))+ ( CONSTS(47)*CONSTS(48))*(1.00000+CONSTS(60)/CONSTS(49)) ALGBRC(72) = ( - ALGBRC(36)*CONSTS(48))*CONSTS(47) ALGBRC(88) = (- ALGBRC(60)+ ( ALGBRC(60)*ALGBRC(60) - ( 4.00000*CONSTS(225))*ALGBRC(72)) ** (1.0 / 2))/( 2.00000*CONSTS(225)) ALGBRC(100) = ALGBRC(8)/(1.00000+ (ALGBRC(88)/CONSTS(47))*(1.00000+CONSTS(60)/CONSTS(48))) ALGBRC(111) = ( ( CONSTS(60)*ALGBRC(88))*ALGBRC(100))/( CONSTS(47)*CONSTS(48)) ALGBRC(120) = ( ALGBRC(100)*ALGBRC(88))/CONSTS(47) ALGBRC(119) = ( CONSTS(60)*ALGBRC(88))/CONSTS(49) ALGBRC(131) = ((CONSTS(400) - STATES(143)) - STATES(142)) - STATES(141) ALGBRC(2) = STATES(4)*1.00000 ** CONSTS(13) ALGBRC(32) = CONSTS(7)*(CONSTS(9)+ALGBRC(2)/(CONSTS(4)+ALGBRC(2))) ALGBRC(68) = ( ALGBRC(32)*CONSTS(330))*CONSTS(228) ALGBRC(84) = ( CONSTS(237)*(STATES(44) - STATES(46)))/CONSTS(310) ALGBRC(97) = ( CONSTS(238)*(STATES(46) - STATES(45)))/CONSTS(310) ALGBRC(108) = ( CONSTS(350)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(46)) ALGBRC(117) = ( (CONSTS(373)+ (CONSTS(175) - 1.00000)*STATES(132))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(46)) ALGBRC(125) = ALGBRC(108)+ALGBRC(117) ALGBRC(135) = (( ( - CONSTS(269)*ALGBRC(131))*STATES(46)+ CONSTS(277)*STATES(143)) - ( CONSTS(270)*STATES(143))*STATES(46))+ CONSTS(278)*STATES(142) ALGBRC(141) = ((CONSTS(381) - STATES(135)) - STATES(134)) - STATES(133) ALGBRC(3) = STATES(2)*1.00000 ** CONSTS(14) ALGBRC(33) = STATES(2)*1.00000 ** CONSTS(15) ALGBRC(57) = ( CONSTS(8)*(CONSTS(10)+ALGBRC(3)/(CONSTS(5)+ALGBRC(3))))*(1.00000 - ( (1.00000 - ( CONSTS(16)*ALGBRC(33))/(CONSTS(6)+ALGBRC(33)))*STATES(1))/(CONSTS(3)+STATES(1))) ALGBRC(83) = ( ALGBRC(57)*CONSTS(321))*CONSTS(228) ALGBRC(96) = ( CONSTS(237)*(STATES(44) - STATES(46)))/CONSTS(307) ALGBRC(106) = ( CONSTS(236)*(STATES(44) - STATES(45)))/CONSTS(307) ALGBRC(115) = ( CONSTS(348)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(44)) ALGBRC(123) = ( (CONSTS(378)+ (CONSTS(175) - 1.00000)*STATES(128))*CONSTS(187))/(1.00000+CONSTS(172)/STATES(44)) ALGBRC(129) = ( (CONSTS(371)+ (CONSTS(175) - 1.00000)*STATES(130))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(44)) ALGBRC(136) = (ALGBRC(115)+ALGBRC(123))+ALGBRC(129) ALGBRC(146) = (( ( - CONSTS(200)*ALGBRC(141))*STATES(44)+ CONSTS(247)*STATES(135)) - ( CONSTS(201)*STATES(135))*STATES(44))+ CONSTS(248)*STATES(134) ALGBRC(151) = ((CONSTS(336) - STATES(139)) - STATES(138)) - STATES(137) ALGBRC(1) = STATES(3)*1.00000 ** CONSTS(13) ALGBRC(56) = CONSTS(7)*(CONSTS(9)+ALGBRC(1)/(CONSTS(4)+ALGBRC(1))) ALGBRC(69) = ( ALGBRC(56)*CONSTS(325))*CONSTS(228) ALGBRC(31) = STATES(3)*1.00000 ** CONSTS(14) ALGBRC(85) = CONSTS(8)*(CONSTS(10)+ALGBRC(31)/(CONSTS(5)+ALGBRC(31))) ALGBRC(107) = ( ALGBRC(85)*CONSTS(329))*CONSTS(228) ALGBRC(116) = ( CONSTS(236)*(STATES(44) - STATES(45)))/CONSTS(308) ALGBRC(124) = ( CONSTS(238)*(STATES(46) - STATES(45)))/CONSTS(308) ALGBRC(130) = ( CONSTS(349)*CONSTS(186))/(1.00000+CONSTS(171)/STATES(45)) ALGBRC(137) = ( (CONSTS(379)+ (CONSTS(175) - 1.00000)*STATES(129))*CONSTS(187))/(1.00000+CONSTS(172)/STATES(45)) ALGBRC(142) = ( (CONSTS(372)+ (CONSTS(175) - 1.00000)*STATES(131))*CONSTS(188))/(1.00000+CONSTS(173)/STATES(45)) ALGBRC(147) = (ALGBRC(130)+ALGBRC(137))+ALGBRC(142) ALGBRC(155) = (( ( - CONSTS(250)*ALGBRC(151))*STATES(45)+ CONSTS(266)*STATES(139)) - ( CONSTS(251)*STATES(139))*STATES(45))+ CONSTS(267)*STATES(138) ALGBRC(7) = ( ( CONSTS(38)*CONSTS(296))*CONSTS(320) - STATES(2)) - STATES(10) ALGBRC(70) = (((CONSTS(259)+CONSTS(272))/CONSTS(224)+ALGBRC(34))+ALGBRC(58)) - ALGBRC(7) ALGBRC(86) = (( CONSTS(272)*(ALGBRC(34) - ALGBRC(7))+ CONSTS(259)*(ALGBRC(58) - ALGBRC(7)))+CONSTS(280))/CONSTS(224) ALGBRC(98) = ( ALGBRC(7)*CONSTS(280))/CONSTS(224) ALGBRC(110) = ((( (- ALGBRC(98)/27.0000)*ALGBRC(70) ** 3.00000 - ( ( ( ALGBRC(70)*ALGBRC(70))*ALGBRC(86))*ALGBRC(86))/108.000)+( ( ALGBRC(70)*ALGBRC(86))*ALGBRC(98))/6.00000)+ALGBRC(86) ** 3.00000/27.0000)+( ALGBRC(98)*ALGBRC(98))/4.00000 ALGBRC(118) = TERNRY(ALGBRC(110).LT.0.00000, - ALGBRC(110) ** (1.0 / 2), 0.00000) ALGBRC(126) = (((TERNRY(ALGBRC(110).GT.0.00000, ALGBRC(110) ** (1.0 / 2), 0.00000))+ALGBRC(98)/2.00000)+( ALGBRC(70)*ALGBRC(86))/6.00000) - ALGBRC(70) ** 3.00000/27.0000 ALGBRC(132) = atan(ALGBRC(118)/ALGBRC(126))/3.00000 ALGBRC(138) = ALGBRC(126)*ALGBRC(126)+ ALGBRC(118)*ALGBRC(118) ** 1.00000/6.00000 ALGBRC(143) = (ALGBRC(86)/3.00000 - ( ALGBRC(70)*ALGBRC(70))/9.00000)/( ALGBRC(138)*ALGBRC(138)) ALGBRC(148) = ( ALGBRC(138)* sin(ALGBRC(132)))*(1.00000+ALGBRC(143)) ALGBRC(152) = ( ALGBRC(138)*cos(ALGBRC(132)))*(1.00000 - ALGBRC(143)) - ALGBRC(70)/3.00000 ALGBRC(157) = ( ALGBRC(152)*ALGBRC(152)+ ALGBRC(148)*ALGBRC(148)) ** (1.0 / 2) ALGBRC(161) = ALGBRC(34)/((1.00000+CONSTS(60)/CONSTS(49))+( ALGBRC(157)*(CONSTS(48)+CONSTS(60)))/( CONSTS(47)*CONSTS(48))) ALGBRC(165) = ( CONSTS(60)*ALGBRC(161))/CONSTS(49) ALGBRC(169) = ( ( CONSTS(60)*ALGBRC(161))*ALGBRC(157))/( CONSTS(47)*CONSTS(48)) ALGBRC(10) = ( ( CONSTS(39)*CONSTS(296))*CONSTS(327) - STATES(4)) - STATES(24) ALGBRC(73) = (((CONSTS(260)+CONSTS(273))/CONSTS(227)+ALGBRC(37))+ALGBRC(61)) - ALGBRC(10) ALGBRC(89) = (( CONSTS(273)*(ALGBRC(37) - ALGBRC(10))+ CONSTS(260)*(ALGBRC(61) - ALGBRC(10)))+CONSTS(281))/CONSTS(227) ALGBRC(101) = ( ALGBRC(10)*CONSTS(281))/CONSTS(227) ALGBRC(113) = ((( (- ALGBRC(101)/27.0000)*ALGBRC(73) ** 3.00000 - ( ( ( ALGBRC(73)*ALGBRC(73))*ALGBRC(89))*ALGBRC(89))/108.000)+( ( ALGBRC(73)*ALGBRC(89))*ALGBRC(101))/6.00000)+ALGBRC(89) ** 3.00000/27.0000)+( ALGBRC(101)*ALGBRC(101))/4.00000 ALGBRC(121) = TERNRY(ALGBRC(113).LT.0.00000, - ALGBRC(113) ** (1.0 / 2), 0.00000) ALGBRC(127) = (((TERNRY(ALGBRC(113).GT.0.00000, ALGBRC(113) ** (1.0 / 2), 0.00000))+ALGBRC(101)/2.00000)+( ALGBRC(73)*ALGBRC(89))/6.00000) - ALGBRC(73) ** 3.00000/27.0000 ALGBRC(133) = atan(ALGBRC(121)/ALGBRC(127))/3.00000 ALGBRC(139) = ALGBRC(127)*ALGBRC(127)+ ALGBRC(121)*ALGBRC(121) ** 1.00000/6.00000 ALGBRC(144) = (ALGBRC(89)/3.00000 - ( ALGBRC(73)*ALGBRC(73))/9.00000)/( ALGBRC(139)*ALGBRC(139)) ALGBRC(149) = ( ALGBRC(139)* sin(ALGBRC(133)))*(1.00000+ALGBRC(144)) ALGBRC(153) = ( ALGBRC(139)*cos(ALGBRC(133)))*(1.00000 - ALGBRC(144)) - ALGBRC(73)/3.00000 ALGBRC(158) = ( ALGBRC(153)*ALGBRC(153)+ ALGBRC(149)*ALGBRC(149)) ** (1.0 / 2) ALGBRC(162) = ALGBRC(37)/((1.00000+CONSTS(60)/CONSTS(49))+( ALGBRC(158)*(CONSTS(48)+CONSTS(60)))/( CONSTS(47)*CONSTS(48))) ALGBRC(166) = ( CONSTS(60)*ALGBRC(162))/CONSTS(49) ALGBRC(170) = ( ( CONSTS(60)*ALGBRC(162))*ALGBRC(158))/( CONSTS(47)*CONSTS(48)) ALGBRC(177) = ALGBRC(58)/((1.00000+CONSTS(60)/CONSTS(54))+( ALGBRC(157)*(CONSTS(53)+CONSTS(60)))/( CONSTS(51)*CONSTS(53))) ALGBRC(181) = ( CONSTS(60)*ALGBRC(177))/CONSTS(54) ALGBRC(185) = ( ( CONSTS(60)*ALGBRC(177))*ALGBRC(157))/( CONSTS(51)*CONSTS(53)) ALGBRC(178) = ALGBRC(61)/((1.00000+CONSTS(60)/CONSTS(54))+( ALGBRC(158)*(CONSTS(53)+CONSTS(60)))/( CONSTS(51)*CONSTS(53))) ALGBRC(182) = ( CONSTS(60)*ALGBRC(178))/CONSTS(54) ALGBRC(186) = ( ( CONSTS(60)*ALGBRC(178))*ALGBRC(158))/( CONSTS(51)*CONSTS(53)) ALGBRC(189) = ALGBRC(169)+ CONSTS(61)*ALGBRC(185) ALGBRC(173) = ( ALGBRC(161)*ALGBRC(157))/CONSTS(47) ALGBRC(193) = ( ALGBRC(177)*ALGBRC(157))/CONSTS(51) ALGBRC(197) = ALGBRC(173)+ CONSTS(61)*ALGBRC(193) ALGBRC(190) = ALGBRC(170)+ CONSTS(64)*ALGBRC(186) ALGBRC(174) = ( ALGBRC(162)*ALGBRC(158))/CONSTS(47) ALGBRC(194) = ( ALGBRC(178)*ALGBRC(158))/CONSTS(51) ALGBRC(198) = ALGBRC(174)+ CONSTS(64)*ALGBRC(194) ALGBRC(134) = CONSTS(275) - STATES(36) ALGBRC(140) = CONSTS(263) - STATES(36)*CONSTS(262) ALGBRC(145) = - CONSTS(231)*STATES(36) ALGBRC(150) = - ALGBRC(134)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(134)*ALGBRC(134) - 3.00000*ALGBRC(140)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(134))*ALGBRC(140) - ( ( 2.00000*ALGBRC(134))*ALGBRC(134))*ALGBRC(134)) - 27.0000*ALGBRC(145))/( 2.00000* ALGBRC(134)*ALGBRC(134) - 3.00000*ALGBRC(140) ** 1.50000))/3.00000) ALGBRC(191) = ( CONSTS(80)*(1.00000 - STATES(43)))/(1.00000+CONSTS(82)/ALGBRC(150)) ALGBRC(195) = ALGBRC(191)+STATES(43) ALGBRC(200) = ALGBRC(195)/(ALGBRC(195)+CONSTS(81)) ALGBRC(201) = 1.00000/(1.00000+CONSTS(81)/ALGBRC(195) ** 2.00000) ALGBRC(199) = (CONSTS(206) - CONSTS(209))+STATES(145) ALGBRC(204) = 0.500000*( (ALGBRC(199) ** 2.00000+ ( 4.00000*CONSTS(206))*CONSTS(209)) ** (1.0 / 2) - ALGBRC(199)) ALGBRC(207) = (CONSTS(29)/CONSTS(320)+CONSTS(83)/CONSTS(327))+ALGBRC(204)/CONSTS(322) ALGBRC(156) = 1.00000/( (1.00000+EXP((4.79800+STATES(51))/- 7.56990))*(1.00000+EXP((25.0000+STATES(51))/- 5.00000))) ALGBRC(196) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(51)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(51)+13.0000))) ALGBRC(203) = ALGBRC(156)/ALGBRC(196) ALGBRC(206) = (1.00000 - ALGBRC(156))/ALGBRC(196) ALGBRC(180) = 1.00000+0.00200000/ALGBRC(150) ** 4.00000 ALGBRC(184) = 6.00000/ALGBRC(180) ALGBRC(164) = 1.00000/(1.00000+EXP((29.9790+STATES(51))/3.17750)) ALGBRC(172) = (0.100000+ALGBRC(164))/1.10000 ALGBRC(215) = 1.00000/( ( 70.0000*(1.00000 - 0.500000*STATES(37)))*(1.00000+EXP((STATES(51)+49.1000)/10.3490))) ALGBRC(217) = 1.00000/(1.00000+EXP((STATES(51)+0.213000)/- 10.8070)) ALGBRC(222) = 1.00000/(ALGBRC(215)+ALGBRC(217)/38.4940) ALGBRC(225) = ALGBRC(172)/ALGBRC(222) ALGBRC(234) = (1.00000 - ALGBRC(172))/ALGBRC(222) ALGBRC(168) = (0.000100000+ALGBRC(164))/1.00010 ALGBRC(160) = 0.100000*STATES(37) ALGBRC(176) = 1.00000+0.0100000/ALGBRC(150) ** 10.0000 ALGBRC(188) = (32.5000 - (18.0000 - ALGBRC(160))/ALGBRC(180)) - 10.0000/ALGBRC(176) ALGBRC(219) = 1.00000/(ALGBRC(215)+ALGBRC(217)/ALGBRC(188)) ALGBRC(228) = ALGBRC(168)/ALGBRC(219) ALGBRC(240) = (1.00000 - ALGBRC(168))/ALGBRC(219) ALGBRC(231) = TERNRY(ABS(ALGBRC(228)).LT.1.00000e-12, 1.00000e-12, ALGBRC(228)) ALGBRC(237) = TERNRY(ABS(ALGBRC(234)).LT.1.00000e-12, 1.00000e-12, ALGBRC(234)) ALGBRC(243) = ( CONSTS(104)*( ( ALGBRC(225)*ALGBRC(240))*ALGBRC(184)))/( ( ALGBRC(237)*ALGBRC(231))*CONSTS(103)) ALGBRC(63) = CONSTS(275) - STATES(34) ALGBRC(75) = CONSTS(263) - STATES(34)*CONSTS(262) ALGBRC(91) = - CONSTS(231)*STATES(34) ALGBRC(103) = - ALGBRC(63)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(63)*ALGBRC(63) - 3.00000*ALGBRC(75)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(63))*ALGBRC(75) - ( ( 2.00000*ALGBRC(63))*ALGBRC(63))*ALGBRC(63)) - 27.0000*ALGBRC(91))/( 2.00000* ALGBRC(63)*ALGBRC(63) - 3.00000*ALGBRC(75) ** 1.50000))/3.00000) ALGBRC(210) = (ALGBRC(150) - ALGBRC(103))/CONSTS(96) ALGBRC(11) = (STATES(5)+CONSTS(393))/CONSTS(18) ALGBRC(221) = (ALGBRC(11) - CONSTS(398))/(0.927300 - CONSTS(398)) ALGBRC(224) = TERNRY(ALGBRC(221).LT.0.00000, 0.00000, ALGBRC(221)) ALGBRC(227) = 0.000257900*(1.00000+ 0.100000*STATES(37)) ALGBRC(230) = EXP( ( 2.00000*STATES(51))*CONSTS(264)) ALGBRC(233) = ( ( ( ( ALGBRC(227)*4.00000)*STATES(51))*CONSTS(276))*( ALGBRC(150)*ALGBRC(230) - 0.341000*CONSTS(98)))/(ALGBRC(230) - 1.00000) ALGBRC(236) = ALGBRC(233)*(STATES(56)+STATES(59)) ALGBRC(239) = 0.000155200*(1.00000+ 0.400000*STATES(37)) ALGBRC(242) = EXP( ( 2.00000*STATES(51))*CONSTS(264)) ALGBRC(245) = ( ( ( ( ALGBRC(239)*4.00000)*STATES(51))*CONSTS(276))*( ALGBRC(150)*ALGBRC(242) - 0.341000*CONSTS(98)))/(ALGBRC(242) - 1.00000) ALGBRC(247) = ALGBRC(245)*(STATES(64)+STATES(67)) ALGBRC(249) = (1.00000 - ALGBRC(224))*ALGBRC(247)+ ALGBRC(224)*ALGBRC(236) ALGBRC(192) = 1.00000/( (1.00000+EXP((13.5600 - STATES(51))/9.45000))*(1.00000+EXP((25.0000+STATES(51))/- 5.00000))) ALGBRC(209) = ALGBRC(192)/ALGBRC(196) ALGBRC(211) = (1.00000 - ALGBRC(192))/ALGBRC(196) ALGBRC(235) = 1.00000+0.00110000/ALGBRC(150) ** 4.00000 ALGBRC(238) = 14.9186/ALGBRC(235) ALGBRC(220) = 1.00000/(1.00000+EXP((17.5000+STATES(51))/3.00000)) ALGBRC(226) = (0.247400+ALGBRC(220))/1.24740 ALGBRC(229) = 1.00000/(ALGBRC(215)+ALGBRC(217)/26.5530) ALGBRC(246) = ALGBRC(226)/ALGBRC(229) ALGBRC(252) = (1.00000 - ALGBRC(226))/ALGBRC(229) ALGBRC(223) = (0.00100000+ALGBRC(220))/1.00100 ALGBRC(213) = 5.00000*STATES(37) ALGBRC(232) = 1.00000+0.0120000/ALGBRC(150) ** 10.0000 ALGBRC(241) = (13.8250 - (6.38360 - ALGBRC(213))/ALGBRC(235)) - 3.36960/ALGBRC(232) ALGBRC(244) = 1.00000/(ALGBRC(215)+ALGBRC(217)/ALGBRC(241)) ALGBRC(248) = ALGBRC(223)/ALGBRC(244) ALGBRC(258) = (1.00000 - ALGBRC(223))/ALGBRC(244) ALGBRC(250) = TERNRY(ABS(ALGBRC(248)).LT.1.00000e-12, 1.00000e-12, ALGBRC(248)) ALGBRC(255) = TERNRY(ABS(ALGBRC(252)).LT.1.00000e-12, 1.00000e-12, ALGBRC(252)) ALGBRC(261) = ( CONSTS(106)*( ( ALGBRC(246)*ALGBRC(258))*ALGBRC(238)))/( ( ALGBRC(255)*ALGBRC(250))*CONSTS(105)) ALGBRC(253) = CONSTS(244)*(1.00000+ 2.00000*STATES(42)) ALGBRC(256) = 0.112500*ALGBRC(253) ALGBRC(114) = (CONSTS(73)+CONSTS(74)) - STATES(35) ALGBRC(122) = STATES(35)*CONSTS(74) ALGBRC(128) = ( ( ALGBRC(114)*ALGBRC(114)+ 4.00000*ALGBRC(122)) ** (1.0 / 2) - ALGBRC(114))/2.00000 ALGBRC(270) = ( ALGBRC(249)*1.00000)/(1.00000+1.00000/ALGBRC(128) ** 8.00000) ALGBRC(272) = ALGBRC(256)*ALGBRC(270) ALGBRC(264) = 1.00000+0.0123000/ALGBRC(128) ALGBRC(266) = ALGBRC(253)/ALGBRC(264) ALGBRC(259) = CONSTS(244)*(1.00000+ 0.00000*STATES(42)) ALGBRC(262) = 0.112500*ALGBRC(259) ALGBRC(273) = ( 1.99250*ALGBRC(262))*ALGBRC(270) ALGBRC(268) = ( 0.535700*ALGBRC(259))/ALGBRC(264) ALGBRC(287) = 1.00000+CONSTS(135)/ALGBRC(150) ** 2.00000 ALGBRC(286) = EXP( ( (CONSTS(140) - 1.00000)*STATES(51))*CONSTS(264)) ALGBRC(288) = 1.00000+ CONSTS(141)*ALGBRC(286) ALGBRC(284) = STATES(50) ** 3.00000 ALGBRC(289) = ( CONSTS(137)*ALGBRC(284)+ CONSTS(365)*ALGBRC(150))+ ( CONSTS(363)*CONSTS(98))*(1.00000+ALGBRC(150)/CONSTS(136)) ALGBRC(290) = ( ( CONSTS(136)*CONSTS(367))*(1.00000+ALGBRC(284)/CONSTS(363))+ ALGBRC(284)*CONSTS(98))+ CONSTS(367)*ALGBRC(150) ALGBRC(285) = EXP( ( CONSTS(140)*STATES(51))*CONSTS(264)) ALGBRC(291) = ( 0.200000*CONSTS(142))*( ( ALGBRC(284)*CONSTS(98))*ALGBRC(285) - ( CONSTS(367)*ALGBRC(150))*ALGBRC(286)) ALGBRC(292) = ALGBRC(291)/( ( ALGBRC(287)*ALGBRC(288))*(ALGBRC(289)+ALGBRC(290))) ALGBRC(208) = (STATES(50) - STATES(49))/CONSTS(95) ALGBRC(39) = (STATES(6)+CONSTS(396))/CONSTS(31) ALGBRC(312) = (ALGBRC(39) - CONSTS(399))/(0.958600 - CONSTS(399)) ALGBRC(314) = TERNRY(ALGBRC(312).LT.0.00000, 0.00000, ALGBRC(312)) ALGBRC(318) = (1.00000 - ALGBRC(314))*STATES(119)+ ALGBRC(314)*STATES(120) ALGBRC(308) = ( CONSTS(156)*EXP(ALGBRC(128)/CONSTS(154)))*(ALGBRC(128) - ALGBRC(150)) ALGBRC(310) = ( CONSTS(157)*EXP(ALGBRC(128)/CONSTS(155)))*(ALGBRC(128) - ALGBRC(150)) ALGBRC(316) = (1.00000 - ALGBRC(314))*ALGBRC(308)+ ALGBRC(314)*ALGBRC(310) ALGBRC(320) = ALGBRC(316)+ALGBRC(318) ALGBRC(212) = (STATES(31) - ALGBRC(128))/CONSTS(97) ALGBRC(154) = (STATES(32) - 0.673519)/(0.999180 - 0.673519) ALGBRC(159) = TERNRY(ALGBRC(154).LT.0.00000, 0.00000, ALGBRC(154)) ALGBRC(163) = (1.00000 - ALGBRC(159))*CONSTS(78)+ ALGBRC(159)*CONSTS(230) ALGBRC(175) = ALGBRC(163)+CONSTS(76) ALGBRC(179) = ((ALGBRC(175) - STATES(33))+CONSTS(72))+CONSTS(79) ALGBRC(167) = ALGBRC(163)*CONSTS(76) ALGBRC(183) = ((ALGBRC(167) - STATES(33)*ALGBRC(175))+ CONSTS(79)*CONSTS(76))+ CONSTS(72)*ALGBRC(163) ALGBRC(171) = - ALGBRC(167)*STATES(33) ALGBRC(187) = - ALGBRC(179)/3.00000+ ( (2.00000/3.00000)* ( ALGBRC(179)*ALGBRC(179) - 3.00000*ALGBRC(183)) ** (1.0 / 2))*cos(acos((( ( 9.00000*ALGBRC(179))*ALGBRC(183) - ( ( 2.00000*ALGBRC(179))*ALGBRC(179))*ALGBRC(179)) - 27.0000*ALGBRC(171))/( 2.00000* ALGBRC(179)*ALGBRC(179) - 3.00000*ALGBRC(183) ** 1.50000))/3.00000) ALGBRC(202) = (ALGBRC(150) - ALGBRC(187))/CONSTS(95) ALGBRC(251) = EXP( STATES(51)*CONSTS(264)) ALGBRC(254) = ( ( ( CONSTS(108)*STATES(51))*CONSTS(276))*(STATES(47) - CONSTS(99)*ALGBRC(251)))/(1.00000 - ALGBRC(251)) ALGBRC(321) = 1.00000 - 1.00000/(1.00000+ALGBRC(318)/CONSTS(109) ** 2.00000) ALGBRC(324) = ( ALGBRC(254)*ALGBRC(321))*STATES(68) ALGBRC(205) = (STATES(48) - STATES(47))/CONSTS(95) ALGBRC(307) = 1.00000/(1.00000+0.0300000/ALGBRC(195) ** 2.00000) ALGBRC(309) = (1.00000 - ALGBRC(307))*CONSTS(162)+ ALGBRC(307)*CONSTS(380) ALGBRC(325) = ( ALGBRC(309)*STATES(31))/CONSTS(165) ALGBRC(311) = (STATES(126) - 0.659100)/(0.994500 - 0.659100) ALGBRC(313) = TERNRY(ALGBRC(311).LT.0.00000, 0.00000, ALGBRC(311)) ALGBRC(315) = ALGBRC(313)*STATES(41) ALGBRC(317) = STATES(41) - ALGBRC(315) ALGBRC(319) = ALGBRC(313) - ALGBRC(315) ALGBRC(322) = ((1.00000 - ALGBRC(319)) - ALGBRC(317)) - ALGBRC(315) ALGBRC(323) = (( ALGBRC(322)*CONSTS(160)+ ALGBRC(319)*CONSTS(375))+ ALGBRC(317)*CONSTS(370))+ ALGBRC(315)*CONSTS(377) ALGBRC(327) = ( ALGBRC(309)*ALGBRC(187))/(ALGBRC(187)+ALGBRC(323)) ALGBRC(329) = ALGBRC(327) - ALGBRC(325) ALGBRC(214) = ( 2.00000*STATES(51))*CONSTS(264) ALGBRC(216) = EXP(ALGBRC(214)) ALGBRC(218) = ( ( ( ( CONSTS(102)*2.00000)*CONSTS(88))*ALGBRC(214))*( ALGBRC(187)*ALGBRC(216) - 0.341000*CONSTS(98)))/(ALGBRC(216) - 1.00000) ALGBRC(293) = 1.00000+CONSTS(135)/ALGBRC(187) ** 2.00000 ALGBRC(294) = 1.00000+ CONSTS(141)*ALGBRC(286) ALGBRC(283) = STATES(49) ** 3.00000 ALGBRC(295) = ( CONSTS(137)*ALGBRC(283)+ CONSTS(365)*ALGBRC(187))+ ( CONSTS(363)*CONSTS(98))*(1.00000+ALGBRC(187)/CONSTS(136)) ALGBRC(296) = ( ( CONSTS(136)*CONSTS(367))*(1.00000+ALGBRC(283)/CONSTS(363))+ ALGBRC(283)*CONSTS(98))+ CONSTS(367)*ALGBRC(187) ALGBRC(297) = ( 0.800000*CONSTS(142))*( ( ALGBRC(283)*CONSTS(98))*ALGBRC(285) - ( CONSTS(367)*ALGBRC(187))*ALGBRC(286)) ALGBRC(298) = ALGBRC(297)/( ( ALGBRC(293)*ALGBRC(294))*(ALGBRC(295)+ALGBRC(296))) ALGBRC(305) = ( CONSTS(152)*ALGBRC(187))/(CONSTS(153)+ALGBRC(187)) ALGBRC(328) = - CONSTS(239)*log(CONSTS(99)/STATES(47)) ALGBRC(332) = CONSTS(239)*log(CONSTS(100)/STATES(127)) ALGBRC(333) = ALGBRC(332) - ALGBRC(328) ALGBRC(334) = ( CONSTS(92)*ALGBRC(333))/(ALGBRC(333)+CONSTS(93)) ALGBRC(353) = CONSTS(239)*log(CONSTS(101)/STATES(49)) ALGBRC(354) = ALGBRC(353) - ALGBRC(328) ** 4.00000 ALGBRC(355) = ( CONSTS(94)*ALGBRC(354))/(ALGBRC(354)+CONSTS(347)) ALGBRC(330) = CONSTS(107)*(STATES(51) - ALGBRC(328)) ALGBRC(356) = ( ( ( CONSTS(242)*STATES(115) ** 3.00000)*STATES(107))*STATES(108))*(STATES(51) - ALGBRC(353)) ALGBRC(357) = ( ( ( CONSTS(242)*STATES(109) ** 3.00000)*STATES(107))*STATES(108))*(STATES(51) - ALGBRC(353)) ALGBRC(358) = ( ( ( CONSTS(242)*STATES(112) ** 3.00000)*STATES(110))*STATES(111))*(STATES(51) - ALGBRC(353)) ALGBRC(359) = ( ( ( ( CONSTS(242)*STATES(115) ** 3.00000)*STATES(113))*STATES(114))*(STATES(51) - ALGBRC(353)))*1.25000 ALGBRC(274) = (STATES(106) - 0.239480)/(0.950143 - 0.239480) ALGBRC(275) = TERNRY(ALGBRC(274).LT.0.00000, 0.00000, ALGBRC(274)) ALGBRC(276) = ALGBRC(275)*STATES(39) ALGBRC(277) = STATES(39) - ALGBRC(276) ALGBRC(278) = ALGBRC(275) - ALGBRC(276) ALGBRC(279) = ((1.00000 - ALGBRC(278)) - ALGBRC(277)) - ALGBRC(276) ALGBRC(360) = (( ALGBRC(279)*ALGBRC(358)+ ALGBRC(278)*ALGBRC(359))+ ALGBRC(277)*ALGBRC(357))+ ALGBRC(276)*ALGBRC(356) ALGBRC(280) = STATES(51)*CONSTS(264) ALGBRC(281) = EXP(ALGBRC(280)) ALGBRC(282) = ( ( ( CONSTS(134)*CONSTS(88))*ALGBRC(280))*( STATES(49)*ALGBRC(281) - CONSTS(101)))/(ALGBRC(281) - 1.00000) ALGBRC(301) = ( CONSTS(145)*CONSTS(243))/(1.00000+EXP( - (STATES(51)+92.0000)*CONSTS(264))) ALGBRC(302) = ALGBRC(301)*STATES(49)/(STATES(49)+CONSTS(148)) ** 3.00000 ALGBRC(303) = ALGBRC(301)*STATES(49)/(STATES(49)+CONSTS(149)) ** 3.00000 ALGBRC(299) = (STATES(116) - 0.126345)/(0.998014 - 0.126345) ALGBRC(300) = TERNRY(ALGBRC(299).LT.0.00000, 0.00000, ALGBRC(299)) ALGBRC(304) = (1.00000 - ALGBRC(300))*ALGBRC(302)+ ALGBRC(300)*ALGBRC(303) ALGBRC(361) = ( STATES(118) ** 3.00000*STATES(117))*(STATES(51) - ALGBRC(353)) ALGBRC(362) = 0.0160000*ALGBRC(361) ALGBRC(363) = 0.00650000*ALGBRC(361) ALGBRC(364) = (1.00000 - STATES(39))*ALGBRC(363)+ STATES(39)*ALGBRC(362) ALGBRC(366) = (((ALGBRC(360)+ALGBRC(282))+ALGBRC(364))+ ALGBRC(304)*3.00000)+ ALGBRC(298)*3.00000 ALGBRC(306) = (((ALGBRC(249)+ALGBRC(218))+ALGBRC(305)) - ALGBRC(298)*2.00000) - ALGBRC(292)*2.00000 ALGBRC(331) = ALGBRC(330)+ALGBRC(324) ALGBRC(335) = STATES(51) - ALGBRC(332) ALGBRC(336) = 1.02000/(1.00000+EXP( 0.238500*(ALGBRC(335) - 59.2150))) ALGBRC(337) = ( 0.491240*EXP( 0.0803200*(ALGBRC(335)+5.47600))+EXP( 0.0617500*(ALGBRC(335) - 594.310)))/(1.00000+EXP( - 0.514300*(ALGBRC(335)+4.75300))) ALGBRC(338) = ( CONSTS(240)*(ALGBRC(336)/(ALGBRC(336)+ALGBRC(337))))*ALGBRC(335) ALGBRC(339) = ALGBRC(338)*1.20000 ALGBRC(340) = (1.00000 - STATES(38))*ALGBRC(338)+ STATES(38)*ALGBRC(339) ALGBRC(257) = 1.00000/(1.00000+EXP((STATES(51)+10.0000)/15.4000)) ALGBRC(341) = ( ( CONSTS(241)*STATES(69))*ALGBRC(257))*(STATES(51) - ALGBRC(332)) ALGBRC(349) = CONSTS(239)*log((CONSTS(100)+ CONSTS(166)*CONSTS(101))/(STATES(127)+ CONSTS(166)*STATES(49))) ALGBRC(260) = 0.195610*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(187) ** 1.40000)) ALGBRC(350) = ( ALGBRC(260)*(STATES(70)+STATES(71)))*(STATES(51) - ALGBRC(349)) ALGBRC(351) = ( ALGBRC(260)*(STATES(72)+STATES(73)))*(STATES(51) - ALGBRC(349)) ALGBRC(263) = (STATES(104)+CONSTS(404))/CONSTS(111) ALGBRC(265) = (ALGBRC(263) - CONSTS(405))/(0.785000 - CONSTS(405)) ALGBRC(267) = TERNRY(ALGBRC(265).LT.0.00000, 0.00000, ALGBRC(265)) ALGBRC(352) = ALGBRC(267)*ALGBRC(351)+ (1.00000 - ALGBRC(267))*ALGBRC(350) ALGBRC(342) = ( CONSTS(124)*(STATES(51) - ALGBRC(332)))/(1.00000+EXP((15.0000 - STATES(51))/17.0000)) ALGBRC(343) = (( CONSTS(124)*(STATES(51) - ALGBRC(332)))/(1.00000+EXP((36.0000 - STATES(51))/17.0000)))*3.62000 ALGBRC(269) = (STATES(105) - 0.0589380)/(0.393747 - 0.0589380) ALGBRC(271) = TERNRY(ALGBRC(269).LT.0.00000, 0.00000, ALGBRC(269)) ALGBRC(344) = (1.00000 - ALGBRC(271))*ALGBRC(342)+ ALGBRC(271)*ALGBRC(343) ALGBRC(326) = EXP(STATES(51)/550.000) ALGBRC(345) = ( ( CONSTS(158)*STATES(121) ** 3.00000)*ALGBRC(326))*(STATES(51) - ALGBRC(332)) ALGBRC(346) = ALGBRC(345)*( 0.735600*STATES(122)+ 0.264400*STATES(124)) ALGBRC(347) = ALGBRC(345)*( 0.735600*STATES(123)+ 0.264400*STATES(125)) ALGBRC(348) = (1.00000 - STATES(40))*ALGBRC(347)+ STATES(40)*ALGBRC(346) ALGBRC(365) = ((((ALGBRC(340)+ALGBRC(341))+ALGBRC(352))+ALGBRC(344))+ALGBRC(348)) - 2.00000*ALGBRC(304) ALGBRC(367) = ALGBRC(366)+ ALGBRC(292)*3.00000 ALGBRC(368) = ((ALGBRC(367)+ALGBRC(365))+ALGBRC(306))+ALGBRC(331) ALGBRC(369) = (TERNRY((VOI - CONSTS(214)) - CONSTS(215)*INT((VOI - CONSTS(214))/CONSTS(215)).LT.CONSTS(213), 1.00000, 0.00000))*CONSTS(257) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END