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