C
C There are a total of 215 entries in the algebraic variable array.
C There are a total of 80 entries in each of the rate and state variable arrays.
C There are a total of 125 entries in the constant variable array.
C
C
C STATES(1) is CaJSR in component calcium (mM).
C STATES(2) is CaNSR in component calcium (mM).
C STATES(3) is CaSS in component calcium (mM).
C STATES(4) is Cai in component calcium (mM).
C STATES(5) is HTRPNCa in component caflux (dimensionless).
C CONSTS(1) is HTRPNtot in component caflux (mM).
C ALGBRC(1) is Jtr in component caflux (mM_per_ms).
C ALGBRC(52) is Jtrpn in component caflux (mM_per_ms).
C ALGBRC(13) is Jxfer in component caflux (mM_per_ms).
C STATES(6) is LTRPNCa in component caflux (dimensionless).
C CONSTS(2) is LTRPNtot in component caflux (mM).
C ALGBRC(26) is caflux_HTRPNCa_a1 in component caflux (mS_per_uF).
C ALGBRC(40) is caflux_LTRPNCa_a1 in component caflux (mS_per_uF).
C CONSTS(3) is khtrpn_minus in component caflux (mS_per_uF).
C CONSTS(4) is khtrpn_plus in component caflux (per_mM_per_ms).
C CONSTS(5) is kltrpn_minus in component caflux (mS_per_uF).
C CONSTS(6) is kltrpn_plus in component caflux (per_mM_per_ms).
C CONSTS(7) is tautr in component caflux (ms).
C CONSTS(8) is tauxfer in component caflux (ms).
C VOI is time in component engine (ms).
C CONSTS(9) is CMDNtot in component calcium (mM).
C CONSTS(10) is CSQNtot in component calcium (mM).
C CONSTS(11) is EGTAtot in component calcium (mM).
C ALGBRC(206) is ICa in component ical (A_per_F).
C ALGBRC(193) is INaCa in component inaca (A_per_F).
C ALGBRC(106) is IpCa in component ipca (A_per_F).
C ALGBRC(212) is Jrel in component ryr (mM_per_ms).
C ALGBRC(201) is Jup in component serca (mM_per_ms).
C CONSTS(12) is KmCMDN in component calcium (mM).
C CONSTS(13) is KmCSQN in component calcium (mM).
C CONSTS(14) is KmEGTA in component calcium (mM).
C CONSTS(15) is VJSR in component cell (uL).
C CONSTS(16) is VNSR in component cell (uL).
C CONSTS(17) is VSS in component cell (uL).
C CONSTS(18) is Vmyo in component cell (uL).
C ALGBRC(67) is beta_JSR in component calcium (dimensionless).
C ALGBRC(82) is beta_SS in component calcium (dimensionless).
C ALGBRC(97) is beta_i in component calcium (dimensionless).
C ALGBRC(214) is calcium_CaSS_a3 in component calcium (mM_per_ms).
C ALGBRC(195) is calcium_Cai_a3 in component calcium (A_per_F).
C ALGBRC(60) is calcium_beta_JSR_b1 in component calcium (dimensionless).
C ALGBRC(72) is calcium_beta_SS_b1 in component calcium (dimensionless).
C ALGBRC(77) is calcium_beta_SS_b2 in component calcium (dimensionless).
C ALGBRC(87) is calcium_beta_i_b1 in component calcium (dimensionless).
C ALGBRC(92) is calcium_beta_i_b2 in component calcium (dimensionless).
C CONSTS(104) is a1 in component cell (s3_A_mol_per_g_per_m5).
C CONSTS(114) is a2 in component cell (s3_A_mol_per_g_per_m5).
C ALGBRC(210) is I in component icat (A_per_F).
C CONSTS(19) is Acap in component cell (cm2).
C CONSTS(20) is F in component phys (C_per_mmol).
C STATES(7) is O2 in component ryr (dimensionless).
C STATES(8) is C1 in component ryr (dimensionless).
C STATES(9) is C2 in component ryr (dimensionless).
C CONSTS(21) is Cao in component extra (mM).
C CONSTS(22) is Ko in component extra (mM).
C CONSTS(23) is Nao in component extra (mM).
C STATES(10) is C0 in component ical (dimensionless).
C ALGBRC(14) is C0_to_C1 in component ical (mS_per_uF).
C ALGBRC(119) is C0_to_CCa0 in component ical (mS_per_uF).
C STATES(11) is C1 in component ical (dimensionless).
C ALGBRC(83) is C1_to_C0 in component ical (mS_per_uF).
C ALGBRC(27) is C1_to_C2 in component ical (mS_per_uF).
C ALGBRC(122) is C1_to_CCa1 in component ical (mS_per_uF).
C STATES(12) is C2 in component ical (dimensionless).
C ALGBRC(88) is C2_to_C1 in component ical (mS_per_uF).
C ALGBRC(33) is C2_to_C3 in component ical (mS_per_uF).
C ALGBRC(127) is C2_to_CCa2 in component ical (mS_per_uF).
C STATES(13) is C3 in component ical (dimensionless).
C ALGBRC(93) is C3_to_C2 in component ical (mS_per_uF).
C ALGBRC(41) is C3_to_C4 in component ical (mS_per_uF).
C ALGBRC(134) is C3_to_CCa3 in component ical (mS_per_uF).
C STATES(14) is C4 in component ical (dimensionless).
C ALGBRC(98) is C4_to_C3 in component ical (mS_per_uF).
C ALGBRC(143) is C4_to_CCa4 in component ical (mS_per_uF).
C STATES(15) is CCa0 in component ical (dimensionless).
C CONSTS(103) is CCa0_to_C0 in component ical (mS_per_uF).
C ALGBRC(53) is CCa0_to_CCa1 in component ical (mS_per_uF).
C STATES(16) is CCa1 in component ical (dimensionless).
C CONSTS(113) is CCa1_to_C1 in component ical (mS_per_uF).
C ALGBRC(104) is CCa1_to_CCa0 in component ical (mS_per_uF).
C ALGBRC(61) is CCa1_to_CCa2 in component ical (mS_per_uF).
C STATES(17) is CCa2 in component ical (dimensionless).
C CONSTS(115) is CCa2_to_C2 in component ical (mS_per_uF).
C ALGBRC(107) is CCa2_to_CCa1 in component ical (mS_per_uF).
C ALGBRC(68) is CCa2_to_CCa3 in component ical (mS_per_uF).
C STATES(18) is CCa3 in component ical (dimensionless).
C CONSTS(117) is CCa3_to_C3 in component ical (mS_per_uF).
C ALGBRC(110) is CCa3_to_CCa2 in component ical (mS_per_uF).
C ALGBRC(73) is CCa3_to_CCa4 in component ical (mS_per_uF).
C STATES(19) is CCa4 in component ical (dimensionless).
C CONSTS(119) is CCa4_to_C4 in component ical (mS_per_uF).
C ALGBRC(113) is CCa4_to_CCa3 in component ical (mS_per_uF).
C ALGBRC(209) is ICaK in component ical (A_per_F).
C CONSTS(116) is ICahalf in component ical (A_per_F).
C ALGBRC(205) is ICamax in component ical (A_per_F).
C STATES(20) is Ki in component potassium (mM).
C STATES(21) is Open in component ical (dimensionless).
C CONSTS(118) is PCa in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(120) is PK in component ical (L_per_F_per_ms_times_1e0).
C ALGBRC(208) is PKprime in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(24) is Pscale in component ical (dimensionless).
C STATES(22) is V in component membrane (mV).
C ALGBRC(204) is VFFRT in component phys (C_per_mmol).
C ALGBRC(181) is VFRT in component phys (dimensionless).
C CONSTS(25) is aL in component ical (dimensionless).
C ALGBRC(2) is alpha in component ical (mS_per_uF).
C ALGBRC(46) is alpha_prime in component ical (mS_per_uF).
C CONSTS(26) is bL in component ical (dimensionless).
C ALGBRC(78) is beta in component ical (mS_per_uF).
C ALGBRC(100) is beta_prime in component ical (mS_per_uF).
C CONSTS(27) is fL in component ical (mS_per_uF).
C CONSTS(28) is gL in component ical (mS_per_uF).
C ALGBRC(116) is gamma in component ical (mS_per_uF).
C ALGBRC(123) is ical_C0_a1 in component ical (mS_per_uF).
C ALGBRC(128) is ical_C0_a2 in component ical (mS_per_uF).
C ALGBRC(129) is ical_C1_a1 in component ical (mS_per_uF).
C ALGBRC(135) is ical_C1_a2 in component ical (mS_per_uF).
C ALGBRC(136) is ical_C2_a1 in component ical (mS_per_uF).
C ALGBRC(141) is ical_C2_a2 in component ical (mS_per_uF).
C ALGBRC(142) is ical_C3_a1 in component ical (mS_per_uF).
C ALGBRC(152) is ical_C3_a2 in component ical (mS_per_uF).
C ALGBRC(153) is ical_C4_a1 in component ical (mS_per_uF).
C ALGBRC(167) is ical_C4_a2 in component ical (mS_per_uF).
C ALGBRC(124) is ical_CCa0_a1 in component ical (mS_per_uF).
C ALGBRC(130) is ical_CCa0_a2 in component ical (mS_per_uF).
C ALGBRC(131) is ical_CCa1_a1 in component ical (mS_per_uF).
C ALGBRC(137) is ical_CCa1_a2 in component ical (mS_per_uF).
C ALGBRC(138) is ical_CCa2_a1 in component ical (mS_per_uF).
C ALGBRC(144) is ical_CCa2_a2 in component ical (mS_per_uF).
C ALGBRC(145) is ical_CCa3_a1 in component ical (mS_per_uF).
C ALGBRC(154) is ical_CCa3_a2 in component ical (mS_per_uF).
C ALGBRC(155) is ical_CCa4_a1 in component ical (mS_per_uF).
C ALGBRC(168) is ical_CCa4_a2 in component ical (mS_per_uF).
C ALGBRC(183) is ical_ICaK_a1 in component ical (mM).
C ALGBRC(185) is ical_ICaK_a2 in component ical (dimensionless).
C ALGBRC(187) is ical_ICamax_a1 in component ical (mM).
C ALGBRC(189) is ical_ICamax_a2 in component ical (dimensionless).
C CONSTS(29) is ical_yCa_yCa_inf_a1 in component ical (dimensionless).
C ALGBRC(207) is imax in component ical (A_per_F).
C CONSTS(30) is omega in component ical (mS_per_uF).
C ALGBRC(3) is tau_yCa in component ical (ms).
C STATES(23) is yCa in component ical (dimensionless).
C ALGBRC(15) is yCa_inf in component ical (dimensionless).
C CONSTS(31) is Ttypescale in component icat (L_per_F_per_ms_times_1e0).
C ALGBRC(4) is icat_l_inf in component icat (dimensionless).
C ALGBRC(16) is icat_l_tau in component icat (ms).
C ALGBRC(5) is icat_n_inf in component icat (dimensionless).
C ALGBRC(17) is icat_n_tau in component icat (ms).
C STATES(24) is l in component icat (dimensionless).
C STATES(25) is n in component icat (dimensionless).
C ALGBRC(115) is EK in component nernst (mV).
C ALGBRC(151) is ENa in component nernst (mV).
C ALGBRC(166) is IHCN in component ihcn (A_per_F).
C CONSTS(32) is IHCNmax in component ihcn (mS_per_uF).
C ALGBRC(182) is h_alpha in component ihcn (mS_per_uF).
C ALGBRC(184) is h_beta in component ihcn (mS_per_uF).
C ALGBRC(186) is h_delta in component ihcn (mS_per_uF).
C CONSTS(33) is h_f in component ihcn (dimensionless).
C ALGBRC(188) is h_gamma in component ihcn (mS_per_uF).
C STATES(26) is hcn1 in component ihcn (dimensionless).
C STATES(27) is hcn10 in component ihcn (dimensionless).
C STATES(28) is hcn2 in component ihcn (dimensionless).
C STATES(29) is hcn3 in component ihcn (dimensionless).
C STATES(30) is hcn4 in component ihcn (dimensionless).
C STATES(31) is hcn5 in component ihcn (dimensionless).
C STATES(32) is hcn6 in component ihcn (dimensionless).
C STATES(33) is hcn7 in component ihcn (dimensionless).
C STATES(34) is hcn8 in component ihcn (dimensionless).
C STATES(35) is hcn9 in component ihcn (dimensionless).
C CONSTS(34) is GK1 in component ik1 (mS_per_uF).
C ALGBRC(118) is IK1 in component ik1 (A_per_F).
C ALGBRC(101) is ik1_IK1_inf in component ik1 (dimensionless).
C CONSTS(35) is A0 in component ikr (mS_per_uF).
C CONSTS(36) is A1 in component ikr (mS_per_uF).
C CONSTS(37) is A2 in component ikr (mS_per_uF).
C CONSTS(38) is A3 in component ikr (mS_per_uF).
C CONSTS(39) is A4 in component ikr (mS_per_uF).
C CONSTS(40) is A5 in component ikr (mS_per_uF).
C CONSTS(41) is A6 in component ikr (mS_per_uF).
C CONSTS(42) is B0 in component ikr (per_mV).
C CONSTS(105) is B1 in component ikr (per_mV).
C CONSTS(43) is B2 in component ikr (per_mV).
C CONSTS(106) is B3 in component ikr (per_mV).
C CONSTS(44) is B4 in component ikr (per_mV).
C CONSTS(107) is B5 in component ikr (per_mV).
C CONSTS(45) is B6 in component ikr (per_mV).
C STATES(36) is C1 in component ikr (dimensionless).
C ALGBRC(6) is C1H_to_C2H in component ikr (mS_per_uF).
C STATES(37) is C2 in component ikr (dimensionless).
C ALGBRC(18) is C2H_to_C1H in component ikr (mS_per_uF).
C CONSTS(108) is C2H_to_C3H in component ikr (mS_per_uF).
C STATES(38) is C3 in component ikr (dimensionless).
C CONSTS(109) is C3H_to_C2H in component ikr (mS_per_uF).
C ALGBRC(7) is C3H_to_IH in component ikr (mS_per_uF).
C ALGBRC(19) is C3H_to_OH in component ikr (mS_per_uF).
C CONSTS(46) is GKr in component ikr (mS_per_uF).
C STATES(39) is I in component ikr (dimensionless).
C ALGBRC(47) is IH_to_C3H in component ikr (mS_per_uF).
C ALGBRC(28) is IH_to_OH in component ikr (mS_per_uF).
C ALGBRC(121) is IKr in component ikr (A_per_F).
C STATES(40) is O in component ikr (dimensionless).
C ALGBRC(34) is OH_to_C3H in component ikr (mS_per_uF).
C ALGBRC(42) is OH_to_IH in component ikr (mS_per_uF).
C CONSTS(47) is T_Const in component ikr (dimensionless).
C CONSTS(110) is fKo in component ikr (dimensionless).
C ALGBRC(29) is ikr_C2_a1 in component ikr (mS_per_uF).
C ALGBRC(35) is ikr_C2_a2 in component ikr (mS_per_uF).
C ALGBRC(54) is ikr_C3_a1 in component ikr (mS_per_uF).
C ALGBRC(62) is ikr_C3_a2 in component ikr (mS_per_uF).
C ALGBRC(55) is ikr_I_a1 in component ikr (mS_per_uF).
C ALGBRC(63) is ikr_I_a2 in component ikr (mS_per_uF).
C ALGBRC(48) is ikr_O_a1 in component ikr (mS_per_uF).
C ALGBRC(56) is ikr_O_a2 in component ikr (mS_per_uF).
C CONSTS(48) is GKs in component iks (mS_per_uF).
C ALGBRC(126) is IKs in component iks (A_per_F).
C ALGBRC(8) is iks_xf_wt_alpha in component iks (mS_per_uF).
C ALGBRC(20) is iks_xf_wt_beta in component iks (mS_per_uF).
C ALGBRC(9) is iks_xs_wt_alpha in component iks (mS_per_uF).
C ALGBRC(21) is iks_xs_wt_beta in component iks (mS_per_uF).
C STATES(41) is xf_wt in component iks (dimensionless).
C STATES(42) is xs_wt in component iks (dimensionless).
C CONSTS(49) is KmCa in component inaca (mM).
C CONSTS(50) is KmNa in component inaca (mM).
C STATES(43) is Nai in component sodium (mM).
C ALGBRC(190) is a1 in component inaca (mol4_per_m12).
C ALGBRC(191) is a2 in component inaca (mol4_per_m12).
C ALGBRC(192) is a3 in component inaca (dimensionless).
C CONSTS(121) is a4 in component inaca (mM).
C CONSTS(123) is a5 in component inaca (m9_per_mol3).
C CONSTS(51) is eta in component inaca (dimensionless).
C CONSTS(52) is kNaCa in component inaca (A_per_F).
C CONSTS(53) is ksat in component inaca (dimensionless).
C CONSTS(122) is nao3 in component inaca (mM3).
C ALGBRC(200) is INaK in component inak (A_per_F).
C CONSTS(54) is INaKmax in component inak (A_per_F).
C CONSTS(55) is KmKo in component inak (mM).
C CONSTS(56) is KmNai in component inak (mM).
C ALGBRC(198) is fNaK in component inak (dimensionless).
C CONSTS(124) is inak_INaK_a1 in component inak (dimensionless).
C ALGBRC(103) is inak_INaK_a2 in component inak (dimensionless).
C ALGBRC(194) is inak_fNaK_a1 in component inak (dimensionless).
C ALGBRC(196) is inak_fNaK_a2 in component inak (dimensionless).
C CONSTS(125) is sigma in component inak (dimensionless).
C CONSTS(57) is IpCamax in component ipca (A_per_F).
C CONSTS(58) is KmpCa in component ipca (mM).
C ALGBRC(133) is Isus in component isus (A_per_F).
C CONSTS(59) is Isusmax in component isus (mS_per_uF).
C STATES(44) is C0 in component ito (dimensionless).
C ALGBRC(22) is C0_to_C1 in component ito (mS_per_uF).
C ALGBRC(146) is C0_to_CI0 in component ito (mS_per_uF).
C STATES(45) is C1 in component ito (dimensionless).
C ALGBRC(105) is C1_to_C0 in component ito (mS_per_uF).
C ALGBRC(30) is C1_to_C2 in component ito (mS_per_uF).
C ALGBRC(147) is C1_to_CI1 in component ito (mS_per_uF).
C STATES(46) is C2 in component ito (dimensionless).
C ALGBRC(108) is C2_to_C1 in component ito (mS_per_uF).
C ALGBRC(36) is C2_to_C3 in component ito (mS_per_uF).
C ALGBRC(148) is C2_to_CI2 in component ito (mS_per_uF).
C STATES(47) is C3 in component ito (dimensionless).
C ALGBRC(111) is C3_to_C2 in component ito (mS_per_uF).
C ALGBRC(149) is C3_to_CI3 in component ito (mS_per_uF).
C ALGBRC(43) is C3_to_O in component ito (mS_per_uF).
C STATES(48) is CI0 in component ito (dimensionless).
C ALGBRC(79) is CI0_to_C0 in component ito (mS_per_uF).
C ALGBRC(49) is CI0_to_CI1 in component ito (mS_per_uF).
C STATES(49) is CI1 in component ito (dimensionless).
C ALGBRC(84) is CI1_to_C1 in component ito (mS_per_uF).
C ALGBRC(114) is CI1_to_CI0 in component ito (mS_per_uF).
C ALGBRC(57) is CI1_to_CI2 in component ito (mS_per_uF).
C STATES(50) is CI2 in component ito (dimensionless).
C ALGBRC(89) is CI2_to_C2 in component ito (mS_per_uF).
C ALGBRC(117) is CI2_to_CI1 in component ito (mS_per_uF).
C ALGBRC(64) is CI2_to_CI3 in component ito (mS_per_uF).
C STATES(51) is CI3 in component ito (dimensionless).
C ALGBRC(94) is CI3_to_C3 in component ito (mS_per_uF).
C ALGBRC(120) is CI3_to_CI2 in component ito (mS_per_uF).
C ALGBRC(69) is CI3_to_OI in component ito (mS_per_uF).
C CONSTS(60) is G in component ito (mS_per_uF).
C ALGBRC(140) is Ito1 in component ito (A_per_F).
C STATES(52) is O in component ito (dimensionless).
C STATES(53) is OI in component ito (dimensionless).
C ALGBRC(125) is OI_to_CI3 in component ito (mS_per_uF).
C ALGBRC(99) is OI_to_O in component ito (mS_per_uF).
C ALGBRC(132) is O_to_C3 in component ito (mS_per_uF).
C ALGBRC(150) is O_to_OI in component ito (mS_per_uF).
C CONSTS(61) is aa in component ito (per_mV).
C CONSTS(62) is ai in component ito (per_mV).
C ALGBRC(10) is alpha_act43 in component ito (mS_per_uF).
C ALGBRC(74) is alpha_inact43 in component ito (mS_per_uF).
C CONSTS(63) is alphaa0 in component ito (mS_per_uF).
C CONSTS(64) is alphai0 in component ito (mS_per_uF).
C CONSTS(65) is b1 in component ito (dimensionless).
C CONSTS(66) is b2 in component ito (dimensionless).
C CONSTS(67) is b3 in component ito (dimensionless).
C CONSTS(68) is b4 in component ito (dimensionless).
C CONSTS(69) is ba in component ito (per_mV).
C ALGBRC(102) is beta_act43 in component ito (mS_per_uF).
C ALGBRC(139) is beta_inact43 in component ito (mS_per_uF).
C CONSTS(70) is betaa0 in component ito (mS_per_uF).
C CONSTS(71) is betai0 in component ito (mS_per_uF).
C CONSTS(72) is bi in component ito (per_mV).
C CONSTS(73) is f1 in component ito (dimensionless).
C CONSTS(74) is f2 in component ito (dimensionless).
C CONSTS(75) is f3 in component ito (dimensionless).
C CONSTS(76) is f4 in component ito (dimensionless).
C ALGBRC(156) is ito_C0_a1 in component ito (mS_per_uF).
C ALGBRC(169) is ito_C0_a2 in component ito (mS_per_uF).
C ALGBRC(157) is ito_C1_a1 in component ito (mS_per_uF).
C ALGBRC(170) is ito_C1_a2 in component ito (mS_per_uF).
C ALGBRC(158) is ito_C2_a1 in component ito (mS_per_uF).
C ALGBRC(171) is ito_C2_a2 in component ito (mS_per_uF).
C ALGBRC(159) is ito_C3_a1 in component ito (mS_per_uF).
C ALGBRC(172) is ito_C3_a2 in component ito (mS_per_uF).
C ALGBRC(160) is ito_CI0_a1 in component ito (mS_per_uF).
C ALGBRC(173) is ito_CI0_a2 in component ito (mS_per_uF).
C ALGBRC(161) is ito_CI1_a1 in component ito (mS_per_uF).
C ALGBRC(174) is ito_CI1_a2 in component ito (mS_per_uF).
C ALGBRC(162) is ito_CI2_a1 in component ito (mS_per_uF).
C ALGBRC(175) is ito_CI2_a2 in component ito (mS_per_uF).
C ALGBRC(163) is ito_CI3_a1 in component ito (mS_per_uF).
C ALGBRC(176) is ito_CI3_a2 in component ito (mS_per_uF).
C ALGBRC(164) is ito_OI_a1 in component ito (mS_per_uF).
C ALGBRC(177) is ito_OI_a2 in component ito (mS_per_uF).
C ALGBRC(165) is ito_O_a1 in component ito (mS_per_uF).
C ALGBRC(178) is ito_O_a2 in component ito (mS_per_uF).
C ALGBRC(180) is INa in component nav15 (A_per_F).
C ALGBRC(179) is INa1 in component nav11 (A_per_F).
C ALGBRC(213) is a1 in component membrane (A_per_F).
C ALGBRC(202) is a2 in component membrane (A_per_F).
C ALGBRC(109) is a3 in component membrane (A_per_F).
C CONSTS(111) is amplitude in component membrane (A_per_F).
C CONSTS(77) is duration in component membrane (ms).
C CONSTS(78) is i_diff in component membrane (A_per_F).
C ALGBRC(215) is i_ion in component membrane (A_per_F).
C ALGBRC(112) is i_stim in component membrane (A_per_F).
C CONSTS(79) is offset in component membrane (ms).
C CONSTS(80) is period in component membrane (ms).
C STATES(54) is BC1 in component nav11 (dimensionless).
C STATES(55) is BC2 in component nav11 (dimensionless).
C STATES(56) is BC3 in component nav11 (dimensionless).
C STATES(57) is BO in component nav11 (dimensionless).
C STATES(58) is C1 in component nav11 (dimensionless).
C STATES(59) is C2 in component nav11 (dimensionless).
C STATES(60) is C3 in component nav11 (dimensionless).
C CONSTS(81) is GNa1 in component nav11 (mS_per_uF).
C STATES(61) is IC2 in component nav11 (dimensionless).
C STATES(62) is IC3 in component nav11 (dimensionless).
C STATES(63) is IF in component nav11 (dimensionless).
C STATES(64) is IS1 in component nav11 (dimensionless).
C STATES(65) is IS2 in component nav11 (dimensionless).
C STATES(66) is O in component nav11 (dimensionless).
C ALGBRC(11) is a11 in component nav11 (mS_per_uF).
C ALGBRC(23) is a12 in component nav11 (mS_per_uF).
C ALGBRC(31) is a13 in component nav11 (mS_per_uF).
C ALGBRC(37) is a2 in component nav11 (mS_per_uF).
C ALGBRC(44) is a3 in component nav11 (mS_per_uF).
C ALGBRC(50) is a4 in component nav11 (mS_per_uF).
C ALGBRC(58) is a5 in component nav11 (mS_per_uF).
C ALGBRC(65) is b11 in component nav11 (mS_per_uF).
C ALGBRC(70) is b12 in component nav11 (mS_per_uF).
C ALGBRC(75) is b13 in component nav11 (mS_per_uF).
C ALGBRC(85) is b2 in component nav11 (mS_per_uF).
C ALGBRC(80) is b3 in component nav11 (mS_per_uF).
C ALGBRC(90) is b4 in component nav11 (mS_per_uF).
C ALGBRC(95) is b5 in component nav11 (mS_per_uF).
C CONSTS(82) is mu1 in component nav11 (mS_per_uF).
C CONSTS(83) is mu2 in component nav11 (mS_per_uF).
C STATES(67) is BC1 in component nav15 (dimensionless).
C STATES(68) is BC2 in component nav15 (dimensionless).
C STATES(69) is BC3 in component nav15 (dimensionless).
C STATES(70) is BO in component nav15 (dimensionless).
C STATES(71) is C1 in component nav15 (dimensionless).
C STATES(72) is C2 in component nav15 (dimensionless).
C STATES(73) is C3 in component nav15 (dimensionless).
C CONSTS(84) is GNa in component nav15 (mS_per_uF).
C STATES(74) is IC2 in component nav15 (dimensionless).
C STATES(75) is IC3 in component nav15 (dimensionless).
C STATES(76) is IF in component nav15 (dimensionless).
C STATES(77) is IS1 in component nav15 (dimensionless).
C STATES(78) is IS2 in component nav15 (dimensionless).
C STATES(79) is O in component nav15 (dimensionless).
C ALGBRC(12) is a11 in component nav15 (mS_per_uF).
C ALGBRC(24) is a12 in component nav15 (mS_per_uF).
C ALGBRC(32) is a13 in component nav15 (mS_per_uF).
C ALGBRC(38) is a2 in component nav15 (mS_per_uF).
C ALGBRC(45) is a3 in component nav15 (mS_per_uF).
C ALGBRC(51) is a4 in component nav15 (mS_per_uF).
C ALGBRC(59) is a5 in component nav15 (mS_per_uF).
C ALGBRC(66) is b11 in component nav15 (mS_per_uF).
C ALGBRC(71) is b12 in component nav15 (mS_per_uF).
C ALGBRC(76) is b13 in component nav15 (mS_per_uF).
C ALGBRC(86) is b2 in component nav15 (mS_per_uF).
C ALGBRC(81) is b3 in component nav15 (mS_per_uF).
C ALGBRC(91) is b4 in component nav15 (mS_per_uF).
C ALGBRC(96) is b5 in component nav15 (mS_per_uF).
C CONSTS(85) is mu1 in component nav15 (mS_per_uF).
C CONSTS(86) is mu2 in component nav15 (mS_per_uF).
C CONSTS(112) is RTF in component phys (mV).
C CONSTS(87) is R in component phys (J_per_mol_per_K).
C CONSTS(88) is T in component phys (kelvin).
C ALGBRC(211) is IK_tot in component potassium (A_per_F).
C STATES(80) is O1 in component ryr (dimensionless).
C CONSTS(89) is kaminus in component ryr (mS_per_uF).
C CONSTS(90) is kaplus in component ryr (m12_per_s_per_mol4_times_1e15).
C CONSTS(91) is kbminus in component ryr (mS_per_uF).
C CONSTS(92) is kbplus in component ryr (m9_per_s_per_mol3_times_1e12).
C CONSTS(93) is kcminus in component ryr (mS_per_uF).
C CONSTS(94) is kcplus in component ryr (mS_per_uF).
C ALGBRC(25) is ryr_C1_a2 in component ryr (mol4_per_m12_times_1e_minus_12).
C ALGBRC(39) is ryr_O2_a1 in component ryr (mM3_times_1e_minus_9).
C CONSTS(95) is v1 in component ryr (mS_per_uF).
C CONSTS(96) is KSR in component serca (dimensionless).
C CONSTS(97) is Kfb in component serca (mM).
C CONSTS(98) is Krb in component serca (mM).
C CONSTS(99) is Nfb in component serca (dimensionless).
C CONSTS(100) is Nrb in component serca (dimensionless).
C ALGBRC(197) is fb in component serca (dimensionless).
C ALGBRC(199) is rb in component serca (dimensionless).
C CONSTS(101) is vmaxf in component serca (mM_per_ms).
C CONSTS(102) is vmaxr in component serca (mM_per_ms).
C ALGBRC(203) is INa_tot in component sodium (A_per_F).
C RATES(5) is d/dt HTRPNCa in component caflux (dimensionless).
C RATES(6) is d/dt LTRPNCa in component caflux (dimensionless).
C RATES(1) is d/dt CaJSR in component calcium (mM).
C RATES(2) is d/dt CaNSR in component calcium (mM).
C RATES(3) is d/dt CaSS in component calcium (mM).
C RATES(4) is d/dt Cai in component calcium (mM).
C RATES(10) is d/dt C0 in component ical (dimensionless).
C RATES(11) is d/dt C1 in component ical (dimensionless).
C RATES(12) is d/dt C2 in component ical (dimensionless).
C RATES(13) is d/dt C3 in component ical (dimensionless).
C RATES(14) is d/dt C4 in component ical (dimensionless).
C RATES(15) is d/dt CCa0 in component ical (dimensionless).
C RATES(16) is d/dt CCa1 in component ical (dimensionless).
C RATES(17) is d/dt CCa2 in component ical (dimensionless).
C RATES(18) is d/dt CCa3 in component ical (dimensionless).
C RATES(19) is d/dt CCa4 in component ical (dimensionless).
C RATES(21) is d/dt Open in component ical (dimensionless).
C RATES(23) is d/dt yCa in component ical (dimensionless).
C RATES(24) is d/dt l in component icat (dimensionless).
C RATES(25) is d/dt n in component icat (dimensionless).
C RATES(26) is d/dt hcn1 in component ihcn (dimensionless).
C RATES(27) is d/dt hcn10 in component ihcn (dimensionless).
C RATES(28) is d/dt hcn2 in component ihcn (dimensionless).
C RATES(29) is d/dt hcn3 in component ihcn (dimensionless).
C RATES(30) is d/dt hcn4 in component ihcn (dimensionless).
C RATES(31) is d/dt hcn5 in component ihcn (dimensionless).
C RATES(32) is d/dt hcn6 in component ihcn (dimensionless).
C RATES(33) is d/dt hcn7 in component ihcn (dimensionless).
C RATES(34) is d/dt hcn8 in component ihcn (dimensionless).
C RATES(35) is d/dt hcn9 in component ihcn (dimensionless).
C RATES(36) is d/dt C1 in component ikr (dimensionless).
C RATES(37) is d/dt C2 in component ikr (dimensionless).
C RATES(38) is d/dt C3 in component ikr (dimensionless).
C RATES(39) is d/dt I in component ikr (dimensionless).
C RATES(40) is d/dt O in component ikr (dimensionless).
C RATES(41) is d/dt xf_wt in component iks (dimensionless).
C RATES(42) is d/dt xs_wt in component iks (dimensionless).
C RATES(44) is d/dt C0 in component ito (dimensionless).
C RATES(45) is d/dt C1 in component ito (dimensionless).
C RATES(46) is d/dt C2 in component ito (dimensionless).
C RATES(47) is d/dt C3 in component ito (dimensionless).
C RATES(48) is d/dt CI0 in component ito (dimensionless).
C RATES(49) is d/dt CI1 in component ito (dimensionless).
C RATES(50) is d/dt CI2 in component ito (dimensionless).
C RATES(51) is d/dt CI3 in component ito (dimensionless).
C RATES(52) is d/dt O in component ito (dimensionless).
C RATES(53) is d/dt OI in component ito (dimensionless).
C RATES(22) is d/dt V in component membrane (mV).
C RATES(54) is d/dt BC1 in component nav11 (dimensionless).
C RATES(55) is d/dt BC2 in component nav11 (dimensionless).
C RATES(56) is d/dt BC3 in component nav11 (dimensionless).
C RATES(57) is d/dt BO in component nav11 (dimensionless).
C RATES(58) is d/dt C1 in component nav11 (dimensionless).
C RATES(59) is d/dt C2 in component nav11 (dimensionless).
C RATES(60) is d/dt C3 in component nav11 (dimensionless).
C RATES(61) is d/dt IC2 in component nav11 (dimensionless).
C RATES(62) is d/dt IC3 in component nav11 (dimensionless).
C RATES(63) is d/dt IF in component nav11 (dimensionless).
C RATES(64) is d/dt IS1 in component nav11 (dimensionless).
C RATES(65) is d/dt IS2 in component nav11 (dimensionless).
C RATES(66) is d/dt O in component nav11 (dimensionless).
C RATES(67) is d/dt BC1 in component nav15 (dimensionless).
C RATES(68) is d/dt BC2 in component nav15 (dimensionless).
C RATES(69) is d/dt BC3 in component nav15 (dimensionless).
C RATES(70) is d/dt BO in component nav15 (dimensionless).
C RATES(71) is d/dt C1 in component nav15 (dimensionless).
C RATES(72) is d/dt C2 in component nav15 (dimensionless).
C RATES(73) is d/dt C3 in component nav15 (dimensionless).
C RATES(74) is d/dt IC2 in component nav15 (dimensionless).
C RATES(75) is d/dt IC3 in component nav15 (dimensionless).
C RATES(76) is d/dt IF in component nav15 (dimensionless).
C RATES(77) is d/dt IS1 in component nav15 (dimensionless).
C RATES(78) is d/dt IS2 in component nav15 (dimensionless).
C RATES(79) is d/dt O in component nav15 (dimensionless).
C RATES(20) is d/dt Ki in component potassium (mM).
C RATES(8) is d/dt C1 in component ryr (dimensionless).
C RATES(9) is d/dt C2 in component ryr (dimensionless).
C RATES(80) is d/dt O1 in component ryr (dimensionless).
C RATES(7) is d/dt O2 in component ryr (dimensionless).
C RATES(43) is d/dt Nai in component sodium (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 2.59679515799999983e-01
      STATES(2) = 2.59898837400000027e-01
      STATES(3) = 1.24655751999999995e-04
      STATES(4) = 7.56225546699999958e-05
      STATES(5) = 9.74145534599999974e-01
      CONSTS(1) = 0.14
      STATES(6) = 7.05428108299999967e-02
      CONSTS(2) = 0.07
      CONSTS(3) = 6.6e-05
      CONSTS(4) = 20.0
      CONSTS(5) = 0.04
      CONSTS(6) = 40.0
      CONSTS(7) = 0.5747
      CONSTS(8) = 26.7
      CONSTS(9) = 0.05
      CONSTS(10) = 15.0
      CONSTS(11) = 0.0
      CONSTS(12) = 0.00238
      CONSTS(13) = 0.8
      CONSTS(14) = 0.00015
      CONSTS(15) = 1.36e-07
      CONSTS(16) = 1.785e-06
      CONSTS(17) = 1.08e-09
      CONSTS(18) = 2.196e-05
      CONSTS(19) = 0.0003912
      CONSTS(20) = 96.5
      STATES(7) = 2.58066283699999986e-09
      STATES(8) = 4.65590390899999984e-01
      STATES(9) = 5.33775261499999987e-01
      CONSTS(21) = 2.0
      CONSTS(22) = 4.0
      CONSTS(23) = 138.0
      STATES(10) = 8.74379838900000039e-01
      STATES(11) = 2.55268330899999993e-02
      STATES(12) = 2.79463949400000000e-04
      STATES(13) = 1.35979496400000004e-06
      STATES(14) = 2.48115507400000008e-09
      STATES(15) = 8.89536543300000065e-02
      STATES(16) = 1.03878563100000005e-02
      STATES(17) = 4.54899679899999998e-04
      STATES(18) = 8.85360913299999938e-06
      STATES(19) = 6.46181672899999939e-08
      STATES(20) = 1.22483450199999993e+02
      STATES(21) = 1.85837559400000008e-10
      CONSTS(24) = 1.8
      STATES(22) = -8.03864711399999976e+01
      CONSTS(25) = 2.0
      CONSTS(26) = 2.0
      CONSTS(27) = 0.3
      CONSTS(28) = 4.0
      CONSTS(29) = 0.82
      CONSTS(30) = 0.0025
      STATES(23) = 9.98998351100000015e-01
      CONSTS(31) = 0.0075614
      STATES(24) = 5.48343768800000020e-01
      STATES(25) = 2.08213191599999990e-03
      CONSTS(32) = 0.3225
      CONSTS(33) = 2.2361
      STATES(26) = 3.45377576399999997e-01
      STATES(27) = 2.64074099199999987e-03
      STATES(28) = 4.08547201099999979e-01
      STATES(29) = 1.80893426200000013e-01
      STATES(30) = 3.54868422399999967e-02
      STATES(31) = 2.59639364999999980e-03
      STATES(32) = 1.13417300900000000e-03
      STATES(33) = 5.28665430399999966e-03
      STATES(34) = 9.66413904500000066e-03
      STATES(35) = 8.13729668700000075e-03
      CONSTS(34) = 0.0226
      CONSTS(35) = 1.71476417330859998e-02
      CONSTS(36) = 3.96932838114099976e-02
      CONSTS(37) = 2.05744860597700009e-02
      CONSTS(38) = 1.34366604422999996e-03
      CONSTS(39) = 1.06663164912879999e-01
      CONSTS(40) = 6.46393910049000014e-03
      CONSTS(41) = 8.03937440300000057e-05
      CONSTS(42) = 3.30460803883500034e-02
      CONSTS(43) = 2.61741271511800010e-02
      CONSTS(44) = 5.68908859717000021e-03
      CONSTS(45) = 6.98089239999999969e-07
      STATES(36) = 7.89638103000000036e-01
      STATES(37) = 7.51317227500000002e-04
      STATES(38) = 1.31802129699999990e-04
      CONSTS(46) = 0.0383724
      STATES(39) = 7.34187676599999975e-06
      STATES(40) = 2.70143786700000017e-05
      CONSTS(47) = 5.32000000100000037e+00
      CONSTS(48) = 2.80819022457099998e-02
      STATES(41) = 7.40047161399999998e-04
      STATES(42) = 1.95694799400000008e-01
      CONSTS(49) = 1.38
      CONSTS(50) = 87.5
      STATES(43) = 1.20687050500000002e+01
      CONSTS(51) = 0.35
      CONSTS(52) = 0.44
      CONSTS(53) = 0.2
      CONSTS(54) = 1.5
      CONSTS(55) = 1.5
      CONSTS(56) = 20.0
      CONSTS(57) = 0.05
      CONSTS(58) = 0.0005
      CONSTS(59) = 0.0919908
      STATES(44) = 9.12446021800000007e-01
      STATES(45) = 5.57395224500000022e-02
      STATES(46) = 1.27681665400000005e-03
      STATES(47) = 1.29938504799999993e-05
      STATES(48) = 1.41240034599999995e-02
      STATES(49) = 1.10576340999999998e-02
      STATES(50) = 4.37517185700000023e-03
      STATES(51) = 9.10523457199999969e-04
      CONSTS(60) = 1.52138159999999995e-01
      STATES(52) = 4.95957872400000005e-08
      STATES(53) = 5.71714766299999999e-05
      CONSTS(61) = 0.028983
      CONSTS(62) = 3.73015999999999994e-04
      CONSTS(63) = 0.543708
      CONSTS(64) = 0.0498424
      CONSTS(65) = 6.77348
      CONSTS(66) = 1.56212705152000009e+01
      CONSTS(67) = 2.87532603313000017e+01
      CONSTS(68) = 5.24576206679000052e+02
      CONSTS(69) = 0.0468437
      CONSTS(70) = 0.080185
      CONSTS(71) = 8.19481999999999958e-04
      CONSTS(72) = 5.374e-08
      CONSTS(73) = 1.8936
      CONSTS(74) = 1.42246474559999996e+01
      CONSTS(75) = 1.58574378389000003e+02
      CONSTS(76) = 1.42936645351000010e+02
      CONSTS(77) = 0.5
      CONSTS(78) = 0.0
      CONSTS(79) = 0.0
      CONSTS(80) = 1000.0
      STATES(54) = 3.66022597099999997e-10
      STATES(55) = 1.72916501500000004e-07
      STATES(56) = 3.12749104299999976e-05
      STATES(57) = 6.60499340399999993e-13
      STATES(58) = 4.69401528400000025e-06
      STATES(59) = 2.21514457499999992e-03
      STATES(60) = 4.00644227700000022e-01
      CONSTS(81) = 9.0
      STATES(61) = 3.35246803199999981e-04
      STATES(62) = 5.99583493199999998e-02
      STATES(63) = 3.65674332999999992e-06
      STATES(64) = 8.87492239200000053e-03
      STATES(65) = 5.26928426500000047e-01
      STATES(66) = 9.08194443799999939e-09
      CONSTS(82) = 4.3e-08
      CONSTS(83) = 0.0003
      STATES(67) = 2.08587473400000016e-08
      STATES(68) = 1.63649800299999998e-06
      STATES(69) = 5.17784203000000017e-05
      STATES(70) = 5.01137449999999989e-11
      STATES(71) = 2.47138734400000005e-04
      STATES(72) = 1.93896494799999999e-02
      STATES(73) = 6.13484026399999993e-01
      CONSTS(84) = 35.0
      STATES(74) = 9.13117952200000020e-03
      STATES(75) = 2.88855704800000002e-01
      STATES(76) = 1.17674846899999997e-04
      STATES(77) = 1.38841821899999998e-03
      STATES(78) = 6.62152358299999966e-02
      STATES(79) = 5.92094839499999983e-07
      CONSTS(85) = 4.3e-08
      CONSTS(86) = 0.0003
      CONSTS(87) = 8.315
      CONSTS(88) = 310.0
      STATES(80) = 6.34345983199999957e-04
      CONSTS(89) = 0.576
      CONSTS(90) = 0.01215
      CONSTS(91) = 1.93
      CONSTS(92) = 0.00405
      CONSTS(93) = 0.0008
      CONSTS(94) = 0.1
      CONSTS(95) = 1.8
      CONSTS(96) = 1.2
      CONSTS(97) = 0.000168
      CONSTS(98) = 3.29
      CONSTS(99) = 1.2
      CONSTS(100) = 1.0
      CONSTS(101) = 7.48e-05
      CONSTS(102) = 0.000318
      CONSTS(103) = CONSTS(30)
      CONSTS(104) = CONSTS(19)/( ( CONSTS(18)*CONSTS(20))*1000.00)
      CONSTS(105) = - 0.0430605
      CONSTS(106) = - 0.0269139
      CONSTS(107) = - 0.0453664
      CONSTS(108) =  CONSTS(47)*0.0260836
      CONSTS(109) =  CONSTS(47)*0.148330
      CONSTS(110) =  (CONSTS(22)/4.00000) ** (1.0 / 2)
      CONSTS(111) = - 100.000
      CONSTS(112) = ( CONSTS(87)*CONSTS(88))/CONSTS(20)
      CONSTS(113) = CONSTS(103)/CONSTS(26)
      CONSTS(114) = CONSTS(19)/( ( ( 2.00000*CONSTS(17))*CONSTS(20))*1000.00)
      CONSTS(115) = CONSTS(113)/CONSTS(26)
      CONSTS(116) = - 0.265000
      CONSTS(117) = CONSTS(115)/CONSTS(26)
      CONSTS(118) =  CONSTS(24)*0.246900
      CONSTS(119) = CONSTS(117)/CONSTS(26)
      CONSTS(120) =  CONSTS(24)*0.000457400
      CONSTS(121) = CONSTS(49)+CONSTS(21)
      CONSTS(122) = CONSTS(23) ** 3.00000
      CONSTS(123) = 5000.00/(CONSTS(50) ** 3.00000+CONSTS(122))
      CONSTS(124) = CONSTS(22)/(CONSTS(22)+CONSTS(55))
      CONSTS(125) = (EXP(CONSTS(23)/67.3000) - 1.00000)/7.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(9) =  CONSTS(94)*STATES(80) -  CONSTS(93)*STATES(9)
      RATES(21) =  CONSTS(27)*STATES(14) -  CONSTS(28)*STATES(21)
      ALGBRC(3) = 1.00000/(0.00336336/(0.500000+EXP(STATES(22)/- 5.53900))+ 0.00779047*EXP(STATES(22)/- 49.5104))
      ALGBRC(15) = CONSTS(29)/(1.00000+EXP((STATES(22)+28.5000)/7.80000))+(1.00000 - CONSTS(29))
      RATES(23) = (ALGBRC(15) - STATES(23))/ALGBRC(3)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(22)+75.6000)/6.20000))
      ALGBRC(16) = TERNRY(STATES(22).LT.- 60.0000, 500.000, 18.3000+ 0.00500000*EXP(- STATES(22)/6.20000))
      RATES(24) = (ALGBRC(4) - STATES(24))/ALGBRC(16)
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(22)+48.4000)/5.20000))
      ALGBRC(17) = TERNRY(STATES(22).LT.- 56.0000,  2.44000*EXP((STATES(22)+120.000)/40.0000), 1.34000+ 0.0350000*EXP(- STATES(22)/11.8000))
      RATES(25) = (ALGBRC(5) - STATES(25))/ALGBRC(17)
      ALGBRC(6) =  ( CONSTS(47)*CONSTS(35))*EXP( CONSTS(42)*STATES(22))
      ALGBRC(18) =  ( CONSTS(47)*CONSTS(36))*EXP( CONSTS(105)*STATES(22))
      RATES(36) =  ALGBRC(18)*STATES(37) -  ALGBRC(6)*STATES(36)
      ALGBRC(8) = TERNRY(STATES(22).EQ.21.0000, 0.000146000/0.0780000, ( 0.000146000*(STATES(22) - 21.0000))/(1.00000 - EXP( - 0.0780000*(STATES(22) - 21.0000))))
      ALGBRC(20) =  0.000910000*EXP( - 0.0280000*STATES(22))
      RATES(41) =  ALGBRC(8)*(1.00000 - STATES(41)) -  ALGBRC(20)*STATES(41)
      ALGBRC(9) = TERNRY(STATES(22).EQ.- 11.0000, 3.30000e-05/0.130000, ( 3.30000e-05*(STATES(22)+11.0000))/(1.00000 - EXP( - 0.130000*(STATES(22)+11.0000))))
      ALGBRC(21) =  0.000100000*EXP( - 0.0150000*STATES(22))
      RATES(42) =  ALGBRC(9)*(1.00000 - STATES(42)) -  ALGBRC(21)*STATES(42)
      ALGBRC(25) =  STATES(3)*1000.00 ** 4.00000
      RATES(8) =  ( - CONSTS(90)*ALGBRC(25))*STATES(8)+ CONSTS(89)*STATES(80)
      ALGBRC(26) =  CONSTS(3)*STATES(5)
      RATES(5) =  ( CONSTS(4)*STATES(4))*(1.00000 - STATES(5)) - ALGBRC(26)
      ALGBRC(29) =  ALGBRC(6)*STATES(36)+ CONSTS(109)*STATES(38)
      ALGBRC(35) =  (ALGBRC(18)+CONSTS(108))*STATES(37)
      RATES(37) = ALGBRC(29) - ALGBRC(35)
      ALGBRC(39) =  STATES(3)*1000.00 ** 3.00000
      RATES(7) =  ( CONSTS(92)*ALGBRC(39))*STATES(80) -  CONSTS(91)*STATES(7)
      ALGBRC(40) =  CONSTS(5)*STATES(6)
      RATES(6) =  ( CONSTS(6)*STATES(4))*(1.00000 - STATES(6)) - ALGBRC(40)
      RATES(80) = - ((RATES(8)+RATES(7))+RATES(9))
      ALGBRC(19) =  ( CONSTS(47)*CONSTS(37))*EXP( CONSTS(43)*STATES(22))
      ALGBRC(28) =  ( CONSTS(47)*CONSTS(40))*EXP( CONSTS(107)*STATES(22))
      ALGBRC(48) =  ALGBRC(19)*STATES(38)+ ALGBRC(28)*STATES(39)
      ALGBRC(34) =  ( CONSTS(47)*CONSTS(38))*EXP( CONSTS(106)*STATES(22))
      ALGBRC(42) =  ( CONSTS(47)*CONSTS(39))*EXP( CONSTS(44)*STATES(22))
      ALGBRC(56) =  (ALGBRC(34)+ALGBRC(42))*STATES(40)
      RATES(40) = ALGBRC(48) - ALGBRC(56)
      ALGBRC(7) =  ( CONSTS(47)*CONSTS(41))*EXP( CONSTS(45)*STATES(22))
      ALGBRC(47) = ( ( ALGBRC(34)*ALGBRC(28))*ALGBRC(7))/( ALGBRC(19)*ALGBRC(42))
      ALGBRC(54) = ( CONSTS(108)*STATES(37)+ ALGBRC(34)*STATES(40))+ ALGBRC(47)*STATES(39)
      ALGBRC(62) =  ((ALGBRC(7)+ALGBRC(19))+CONSTS(109))*STATES(38)
      RATES(38) = ALGBRC(54) - ALGBRC(62)
      ALGBRC(55) =  ALGBRC(7)*STATES(38)+ ALGBRC(42)*STATES(40)
      ALGBRC(63) =  (ALGBRC(47)+ALGBRC(28))*STATES(39)
      RATES(39) = ALGBRC(55) - ALGBRC(63)
      ALGBRC(11) = 2.80200/( 0.210000*EXP(- STATES(22)/17.0000)+ 0.230000*EXP(- STATES(22)/150.000))
      ALGBRC(65) =  0.400000*EXP(- STATES(22)/20.3000)
      RATES(56) = ( STATES(55)*ALGBRC(65)+ STATES(60)*CONSTS(82)) -  STATES(56)*(ALGBRC(11)+CONSTS(83))
      ALGBRC(12) = 3.80200/( 0.102700*EXP(- STATES(22)/17.0000)+ 0.200000*EXP(- STATES(22)/150.000))
      ALGBRC(66) =  0.191700*EXP(- STATES(22)/20.3000)
      RATES(69) = ( STATES(68)*ALGBRC(66)+ STATES(73)*CONSTS(85)) -  STATES(69)*(ALGBRC(12)+CONSTS(86))
      ALGBRC(23) = 2.80200/( 0.230000*EXP(- STATES(22)/15.0000)+ 0.250000*EXP(- STATES(22)/150.000))
      ALGBRC(70) =  0.400000*EXP(- (STATES(22) - 5.00000)/20.3000)
      RATES(55) = (( STATES(54)*ALGBRC(70)+ STATES(56)*ALGBRC(11))+ STATES(59)*CONSTS(82)) -  STATES(55)*((ALGBRC(23)+ALGBRC(65))+CONSTS(83))
      ALGBRC(24) = 3.80200/( 0.102700*EXP(- STATES(22)/15.0000)+ 0.230000*EXP(- STATES(22)/150.000))
      ALGBRC(71) =  0.200000*EXP(- (STATES(22) - 5.00000)/20.3000)
      RATES(68) = (( STATES(67)*ALGBRC(71)+ STATES(69)*ALGBRC(12))+ STATES(72)*CONSTS(85)) -  STATES(68)*((ALGBRC(24)+ALGBRC(66))+CONSTS(86))
      ALGBRC(31) = 2.80200/( 0.250000*EXP(- STATES(22)/12.0000)+ 0.270000*EXP(- STATES(22)/150.000))
      ALGBRC(75) = ( 0.400000*EXP(- (STATES(22) - 10.0000)/20.3000))/4.50000
      RATES(54) = (( STATES(57)*ALGBRC(75)+ STATES(55)*ALGBRC(23))+ STATES(58)*CONSTS(82)) -  STATES(54)*((ALGBRC(31)+ALGBRC(70))+CONSTS(83))
      RATES(57) = ( STATES(54)*ALGBRC(31)+ STATES(66)*CONSTS(82)) -  STATES(57)*(ALGBRC(75)+CONSTS(83))
      ALGBRC(32) = 3.80200/( 0.102700*EXP(- STATES(22)/12.0000)+ 0.250000*EXP(- STATES(22)/150.000))
      ALGBRC(76) =  0.220000*EXP(- (STATES(22) - 10.0000)/20.3000)
      RATES(67) = (( STATES(70)*ALGBRC(76)+ STATES(68)*ALGBRC(24))+ STATES(71)*CONSTS(85)) -  STATES(67)*((ALGBRC(32)+ALGBRC(71))+CONSTS(86))
      RATES(70) = ( STATES(67)*ALGBRC(32)+ STATES(79)*CONSTS(85)) -  STATES(70)*(ALGBRC(76)+CONSTS(86))
      ALGBRC(44) =  ( 3.79330e-07*EXP(- STATES(22)/7.60000))*3.00000
      ALGBRC(80) = 0.00840000+ 2.00000e-05*STATES(22)
      RATES(58) = ((( STATES(63)*ALGBRC(44)+ STATES(54)*CONSTS(83))+ STATES(66)*ALGBRC(75))+ STATES(59)*ALGBRC(23)) -  STATES(58)*(((ALGBRC(31)+ALGBRC(80))+CONSTS(82))+ALGBRC(70))
      RATES(59) = ((( STATES(61)*ALGBRC(44)+ STATES(55)*CONSTS(83))+ STATES(58)*ALGBRC(70))+ STATES(60)*ALGBRC(11)) -  STATES(59)*(((ALGBRC(23)+ALGBRC(80))+CONSTS(82))+ALGBRC(65))
      RATES(60) = (( STATES(62)*ALGBRC(44)+ STATES(56)*CONSTS(83))+ STATES(59)*ALGBRC(65)) -  STATES(60)*((ALGBRC(11)+ALGBRC(80))+CONSTS(82))
      RATES(61) = (( STATES(63)*ALGBRC(70)+ STATES(59)*ALGBRC(80))+ STATES(62)*ALGBRC(11)) -  STATES(61)*((ALGBRC(23)+ALGBRC(44))+ALGBRC(65))
      RATES(62) = ( STATES(61)*ALGBRC(65)+ STATES(60)*ALGBRC(80)) -  STATES(62)*(ALGBRC(11)+ALGBRC(44))
      ALGBRC(45) =  3.79330e-07*EXP(- STATES(22)/7.70000)
      ALGBRC(81) = 0.00840000+ 2.00000e-05*STATES(22)
      RATES(71) = ((( STATES(76)*ALGBRC(45)+ STATES(67)*CONSTS(86))+ STATES(79)*ALGBRC(76))+ STATES(72)*ALGBRC(24)) -  STATES(71)*(((ALGBRC(32)+ALGBRC(81))+CONSTS(85))+ALGBRC(71))
      RATES(72) = ((( STATES(74)*ALGBRC(45)+ STATES(68)*CONSTS(86))+ STATES(71)*ALGBRC(71))+ STATES(73)*ALGBRC(12)) -  STATES(72)*(((ALGBRC(24)+ALGBRC(81))+CONSTS(85))+ALGBRC(66))
      RATES(73) = (( STATES(75)*ALGBRC(45)+ STATES(69)*CONSTS(86))+ STATES(72)*ALGBRC(66)) -  STATES(73)*((ALGBRC(12)+ALGBRC(81))+CONSTS(85))
      RATES(74) = (( STATES(76)*ALGBRC(71)+ STATES(72)*ALGBRC(81))+ STATES(75)*ALGBRC(12)) -  STATES(74)*((ALGBRC(24)+ALGBRC(45))+ALGBRC(66))
      RATES(75) = ( STATES(74)*ALGBRC(66)+ STATES(73)*ALGBRC(81)) -  STATES(75)*(ALGBRC(12)+ALGBRC(45))
      ALGBRC(37) = ( 9.17800*EXP(STATES(22)/29.6800))/4.50000
      ALGBRC(85) = ( ( ALGBRC(44)*ALGBRC(37))*ALGBRC(31))/( ALGBRC(80)*ALGBRC(75))
      RATES(66) = (( STATES(63)*ALGBRC(85)+ STATES(57)*CONSTS(83))+ STATES(58)*ALGBRC(31)) -  STATES(66)*((ALGBRC(37)+ALGBRC(75))+CONSTS(82))
      ALGBRC(38) =  9.17800*EXP(STATES(22)/29.6800)
      ALGBRC(86) = ( ( ALGBRC(45)*ALGBRC(38))*ALGBRC(32))/( ALGBRC(81)*ALGBRC(76))
      RATES(79) = (( STATES(76)*ALGBRC(86)+ STATES(70)*CONSTS(86))+ STATES(71)*ALGBRC(32)) -  STATES(79)*((ALGBRC(38)+ALGBRC(76))+CONSTS(85))
      ALGBRC(50) =  ( (ALGBRC(37)/100.000)*1.50000)*0.285000
      ALGBRC(90) = ALGBRC(44)/5.00000
      RATES(63) = ((( STATES(64)*ALGBRC(90)+ STATES(66)*ALGBRC(37))+ STATES(61)*ALGBRC(23))+ STATES(58)*ALGBRC(80)) -  STATES(63)*(((ALGBRC(50)+ALGBRC(85))+ALGBRC(70))+ALGBRC(44))
      ALGBRC(51) = ALGBRC(38)/100.000
      ALGBRC(91) = ALGBRC(45)
      RATES(76) = ((( STATES(77)*ALGBRC(91)+ STATES(79)*ALGBRC(38))+ STATES(74)*ALGBRC(24))+ STATES(71)*ALGBRC(81)) -  STATES(76)*(((ALGBRC(51)+ALGBRC(86))+ALGBRC(71))+ALGBRC(45))
      ALGBRC(58) =  (ALGBRC(37)/95000.0)*80.0000
      ALGBRC(95) = (ALGBRC(44)/30.0000)/10.0000
      RATES(64) = ( STATES(65)*ALGBRC(95)+ STATES(63)*ALGBRC(50)) -  STATES(64)*(ALGBRC(58)+ALGBRC(90))
      RATES(65) =  STATES(64)*ALGBRC(58) -  STATES(65)*ALGBRC(95)
      ALGBRC(59) = ALGBRC(38)/95000.0
      ALGBRC(96) = ALGBRC(45)/50.0000
      RATES(77) = ( STATES(78)*ALGBRC(96)+ STATES(76)*ALGBRC(51)) -  STATES(77)*(ALGBRC(59)+ALGBRC(91))
      RATES(78) =  STATES(77)*ALGBRC(59) -  STATES(78)*ALGBRC(96)
      ALGBRC(2) =  ( ( 4.00000*1.20000)*0.416000)*EXP( 0.0120000*(STATES(22) - 35.0000))
      ALGBRC(14) =  4.00000*ALGBRC(2)
      ALGBRC(116) =  ( 0.600000*0.0923300)*STATES(3)
      ALGBRC(119) = ALGBRC(116)
      ALGBRC(123) =  (ALGBRC(14)+ALGBRC(119))*STATES(10)
      ALGBRC(78) =  ( ( 4.00000*0.450000)*0.0490000)*EXP( - 0.0650000*(STATES(22) - 22.0000))
      ALGBRC(83) = ALGBRC(78)
      ALGBRC(128) =  ALGBRC(83)*STATES(11)+ CONSTS(103)*STATES(15)
      RATES(10) = ALGBRC(128) - ALGBRC(123)
      ALGBRC(46) =  CONSTS(25)*ALGBRC(2)
      ALGBRC(53) =  4.00000*ALGBRC(46)
      ALGBRC(124) =  (ALGBRC(53)+CONSTS(103))*STATES(15)
      ALGBRC(100) = ALGBRC(78)/CONSTS(26)
      ALGBRC(104) = ALGBRC(100)
      ALGBRC(130) =  ALGBRC(104)*STATES(16)+ ALGBRC(119)*STATES(10)
      RATES(15) = ALGBRC(130) - ALGBRC(124)
      ALGBRC(27) =  3.00000*ALGBRC(2)
      ALGBRC(122) =  CONSTS(25)*ALGBRC(119)
      ALGBRC(129) =  ((ALGBRC(83)+ALGBRC(27))+ALGBRC(122))*STATES(11)
      ALGBRC(88) =  2.00000*ALGBRC(78)
      ALGBRC(135) = ( ALGBRC(14)*STATES(10)+ ALGBRC(88)*STATES(12))+ CONSTS(113)*STATES(16)
      RATES(11) = ALGBRC(135) - ALGBRC(129)
      ALGBRC(61) =  3.00000*ALGBRC(46)
      ALGBRC(131) =  ((ALGBRC(104)+ALGBRC(61))+CONSTS(113))*STATES(16)
      ALGBRC(107) =  2.00000*ALGBRC(100)
      ALGBRC(137) = ( ALGBRC(53)*STATES(15)+ ALGBRC(107)*STATES(17))+ ALGBRC(122)*STATES(11)
      RATES(16) = ALGBRC(137) - ALGBRC(131)
      ALGBRC(33) =  2.00000*ALGBRC(2)
      ALGBRC(127) =  CONSTS(25)*ALGBRC(122)
      ALGBRC(136) =  ((ALGBRC(88)+ALGBRC(33))+ALGBRC(127))*STATES(12)
      ALGBRC(93) =  3.00000*ALGBRC(78)
      ALGBRC(141) = ( ALGBRC(27)*STATES(11)+ ALGBRC(93)*STATES(13))+ CONSTS(115)*STATES(17)
      RATES(12) = ALGBRC(141) - ALGBRC(136)
      ALGBRC(68) =  2.00000*ALGBRC(46)
      ALGBRC(138) =  ((ALGBRC(107)+ALGBRC(68))+CONSTS(115))*STATES(17)
      ALGBRC(110) =  3.00000*ALGBRC(100)
      ALGBRC(144) = ( ALGBRC(61)*STATES(16)+ ALGBRC(110)*STATES(18))+ ALGBRC(127)*STATES(12)
      RATES(17) = ALGBRC(144) - ALGBRC(138)
      ALGBRC(41) = ALGBRC(2)
      ALGBRC(134) =  CONSTS(25)*ALGBRC(127)
      ALGBRC(142) =  ((ALGBRC(93)+ALGBRC(41))+ALGBRC(134))*STATES(13)
      ALGBRC(98) =  4.00000*ALGBRC(78)
      ALGBRC(152) = ( ALGBRC(33)*STATES(12)+ ALGBRC(98)*STATES(14))+ CONSTS(117)*STATES(18)
      RATES(13) = ALGBRC(152) - ALGBRC(142)
      ALGBRC(73) = ALGBRC(46)
      ALGBRC(145) =  ((ALGBRC(110)+ALGBRC(73))+CONSTS(117))*STATES(18)
      ALGBRC(113) =  4.00000*ALGBRC(100)
      ALGBRC(154) = ( ALGBRC(68)*STATES(17)+ ALGBRC(113)*STATES(19))+ ALGBRC(134)*STATES(13)
      RATES(18) = ALGBRC(154) - ALGBRC(145)
      ALGBRC(143) =  CONSTS(25)*ALGBRC(134)
      ALGBRC(153) =  ((ALGBRC(98)+CONSTS(27))+ALGBRC(143))*STATES(14)
      ALGBRC(167) = ( ALGBRC(41)*STATES(13)+ CONSTS(28)*STATES(21))+ CONSTS(119)*STATES(19)
      RATES(14) = ALGBRC(167) - ALGBRC(153)
      ALGBRC(155) =  (ALGBRC(113)+CONSTS(119))*STATES(19)
      ALGBRC(168) =  ALGBRC(73)*STATES(18)+ ALGBRC(143)*STATES(14)
      RATES(19) = ALGBRC(168) - ALGBRC(155)
      ALGBRC(10) =  CONSTS(63)*EXP( CONSTS(61)*STATES(22))
      ALGBRC(22) =  4.00000*ALGBRC(10)
      ALGBRC(139) =  CONSTS(71)*EXP( CONSTS(72)*STATES(22))
      ALGBRC(146) = ALGBRC(139)
      ALGBRC(156) =  (ALGBRC(22)+ALGBRC(146))*STATES(44)
      ALGBRC(102) =  CONSTS(70)*EXP( - CONSTS(69)*STATES(22))
      ALGBRC(105) = ALGBRC(102)
      ALGBRC(74) =  CONSTS(64)*EXP( - CONSTS(62)*STATES(22))
      ALGBRC(79) = ALGBRC(74)
      ALGBRC(169) =  ALGBRC(105)*STATES(45)+ ALGBRC(79)*STATES(48)
      RATES(44) = ALGBRC(169) - ALGBRC(156)
      ALGBRC(30) =  3.00000*ALGBRC(10)
      ALGBRC(147) =  CONSTS(73)*ALGBRC(139)
      ALGBRC(157) =  ((ALGBRC(30)+ALGBRC(105))+ALGBRC(147))*STATES(45)
      ALGBRC(108) =  2.00000*ALGBRC(102)
      ALGBRC(84) = ALGBRC(74)/CONSTS(65)
      ALGBRC(170) = ( ALGBRC(108)*STATES(46)+ ALGBRC(84)*STATES(49))+ ALGBRC(22)*STATES(44)
      RATES(45) = ALGBRC(170) - ALGBRC(157)
      ALGBRC(36) =  2.00000*ALGBRC(10)
      ALGBRC(148) =  CONSTS(74)*ALGBRC(139)
      ALGBRC(158) =  ((ALGBRC(36)+ALGBRC(108))+ALGBRC(148))*STATES(46)
      ALGBRC(111) =  3.00000*ALGBRC(102)
      ALGBRC(89) = ALGBRC(74)/CONSTS(66)
      ALGBRC(171) = ( ALGBRC(111)*STATES(47)+ ALGBRC(89)*STATES(50))+ ALGBRC(30)*STATES(45)
      RATES(46) = ALGBRC(171) - ALGBRC(158)
      ALGBRC(149) =  CONSTS(75)*ALGBRC(139)
      ALGBRC(43) = ALGBRC(10)
      ALGBRC(159) =  ((ALGBRC(43)+ALGBRC(111))+ALGBRC(149))*STATES(47)
      ALGBRC(94) = ALGBRC(74)/CONSTS(67)
      ALGBRC(132) =  4.00000*ALGBRC(102)
      ALGBRC(172) = ( ALGBRC(132)*STATES(52)+ ALGBRC(94)*STATES(51))+ ALGBRC(36)*STATES(46)
      RATES(47) = ALGBRC(172) - ALGBRC(159)
      ALGBRC(49) =  ( 4.00000*CONSTS(65))*ALGBRC(10)
      ALGBRC(160) =  (ALGBRC(79)+ALGBRC(49))*STATES(48)
      ALGBRC(114) = ALGBRC(102)/CONSTS(73)
      ALGBRC(173) =  ALGBRC(146)*STATES(44)+ ALGBRC(114)*STATES(49)
      RATES(48) = ALGBRC(173) - ALGBRC(160)
      ALGBRC(57) = ( ( 3.00000*CONSTS(66))*ALGBRC(10))/CONSTS(65)
      ALGBRC(161) =  ((ALGBRC(57)+ALGBRC(84))+ALGBRC(114))*STATES(49)
      ALGBRC(117) = ( ( 2.00000*CONSTS(73))*ALGBRC(102))/CONSTS(74)
      ALGBRC(174) = ( ALGBRC(117)*STATES(50)+ ALGBRC(147)*STATES(45))+ ALGBRC(49)*STATES(48)
      RATES(49) = ALGBRC(174) - ALGBRC(161)
      ALGBRC(64) = ( ( 2.00000*CONSTS(67))*ALGBRC(10))/CONSTS(66)
      ALGBRC(162) =  ((ALGBRC(64)+ALGBRC(89))+ALGBRC(117))*STATES(50)
      ALGBRC(120) = ( ( 3.00000*CONSTS(74))*ALGBRC(102))/CONSTS(75)
      ALGBRC(175) = ( ALGBRC(120)*STATES(51)+ ALGBRC(148)*STATES(46))+ ALGBRC(57)*STATES(49)
      RATES(50) = ALGBRC(175) - ALGBRC(162)
      ALGBRC(69) = ( CONSTS(68)*ALGBRC(10))/CONSTS(67)
      ALGBRC(163) =  ((ALGBRC(69)+ALGBRC(94))+ALGBRC(120))*STATES(51)
      ALGBRC(125) = ( ( 4.00000*CONSTS(75))*ALGBRC(102))/CONSTS(76)
      ALGBRC(176) = ( ALGBRC(125)*STATES(53)+ ALGBRC(149)*STATES(47))+ ALGBRC(64)*STATES(50)
      RATES(51) = ALGBRC(176) - ALGBRC(163)
      ALGBRC(150) =  CONSTS(76)*ALGBRC(139)
      ALGBRC(165) =  (ALGBRC(132)+ALGBRC(150))*STATES(52)
      ALGBRC(99) = ALGBRC(74)/CONSTS(68)
      ALGBRC(178) =  ALGBRC(43)*STATES(47)+ ALGBRC(99)*STATES(53)
      RATES(52) = ALGBRC(178) - ALGBRC(165)
      ALGBRC(164) =  (ALGBRC(99)+ALGBRC(125))*STATES(53)
      ALGBRC(177) =  ALGBRC(150)*STATES(52)+ ALGBRC(69)*STATES(51)
      RATES(53) = ALGBRC(177) - ALGBRC(164)
      ALGBRC(181) = STATES(22)/CONSTS(112)
      ALGBRC(182) = ( 0.000171200*EXP( - 1.46500*ALGBRC(181)))/1000.00
      ALGBRC(184) = ( 26.1700*EXP( 1.46500*ALGBRC(181)))/1000.00
      ALGBRC(186) = ( 287.500*EXP( 1.24200*ALGBRC(181)))/1000.00
      ALGBRC(188) = ( 0.0402500*EXP( - 1.24200*ALGBRC(181)))/1000.00
      RATES(26) = ( - ( 4.00000*ALGBRC(188)+ALGBRC(182))*STATES(26)+ ALGBRC(186)*STATES(28))+ ALGBRC(184)*STATES(32)
      RATES(27) = ( - (ALGBRC(184)/CONSTS(33) ** 4.00000+( 4.00000*ALGBRC(186))/CONSTS(33))*STATES(27)+ ( ALGBRC(182)*CONSTS(33) ** 4.00000)*STATES(31))+ ( ALGBRC(188)*CONSTS(33))*STATES(35)
      RATES(28) = (( - ((ALGBRC(186)+ 3.00000*ALGBRC(188))+ ALGBRC(182)*CONSTS(33))*STATES(28)+ ( 4.00000*ALGBRC(188))*STATES(26))+ ( 2.00000*ALGBRC(186))*STATES(29))+ (ALGBRC(184)/CONSTS(33))*STATES(33)
      RATES(29) = (( - (( 2.00000*ALGBRC(186)+ 2.00000*ALGBRC(188))+ ALGBRC(182)*CONSTS(33) ** 2.00000)*STATES(29)+ ( 3.00000*ALGBRC(188))*STATES(28))+ ( 3.00000*ALGBRC(186))*STATES(30))+ (ALGBRC(184)/CONSTS(33) ** 2.00000)*STATES(34)
      RATES(30) = (( - (( 3.00000*ALGBRC(186)+ALGBRC(188))+ ALGBRC(182)*CONSTS(33) ** 3.00000)*STATES(30)+ ( 2.00000*ALGBRC(188))*STATES(29))+ ( 4.00000*ALGBRC(186))*STATES(31))+ (ALGBRC(184)/CONSTS(33) ** 3.00000)*STATES(35)
      RATES(31) = ( - ( 4.00000*ALGBRC(186)+ ALGBRC(182)*CONSTS(33) ** 4.00000)*STATES(31)+ ALGBRC(188)*STATES(30))+ (ALGBRC(184)/CONSTS(33) ** 4.00000)*STATES(27)
      RATES(32) = ( - (ALGBRC(184)+ ( 4.00000*ALGBRC(188))*CONSTS(33))*STATES(32)+ ALGBRC(182)*STATES(26))+ (ALGBRC(186)/CONSTS(33))*STATES(33)
      RATES(33) = (( - ((ALGBRC(184)/CONSTS(33)+ALGBRC(186)/CONSTS(33))+ ( 3.00000*ALGBRC(188))*CONSTS(33))*STATES(33)+ ( ALGBRC(182)*CONSTS(33))*STATES(28))+ ( ( 4.00000*ALGBRC(188))*CONSTS(33))*STATES(32))+ (( 2.00000*ALGBRC(186))/CONSTS(33))*STATES(34)
      RATES(34) = (( - ((ALGBRC(184)/CONSTS(33) ** 2.00000+( 2.00000*ALGBRC(186))/CONSTS(33))+ ( 2.00000*ALGBRC(188))*CONSTS(33))*STATES(34)+ ( ALGBRC(182)*CONSTS(33) ** 2.00000)*STATES(29))+ ( ( 3.00000*ALGBRC(188))*CONSTS(33))*STATES(33))+ (( 3.00000*ALGBRC(186))/CONSTS(33))*STATES(35)
      RATES(35) = (( - ((ALGBRC(184)/CONSTS(33) ** 3.00000+( 3.00000*ALGBRC(186))/CONSTS(33))+ ALGBRC(188)*CONSTS(33))*STATES(35)+ ( ALGBRC(182)*CONSTS(33) ** 3.00000)*STATES(30))+ ( ( 2.00000*ALGBRC(188))*CONSTS(33))*STATES(34))+ (( 4.00000*ALGBRC(186))/CONSTS(33))*STATES(27)
      ALGBRC(1) = (STATES(2) - STATES(1))/CONSTS(7)
      ALGBRC(197) = STATES(4)/CONSTS(97) ** CONSTS(99)
      ALGBRC(199) = STATES(2)/CONSTS(98) ** CONSTS(100)
      ALGBRC(201) = ( CONSTS(96)*( CONSTS(101)*ALGBRC(197) -  CONSTS(102)*ALGBRC(199)))/((1.00000+ALGBRC(197))+ALGBRC(199))
      RATES(2) = ( ALGBRC(201)*CONSTS(18))/CONSTS(16) - ( ALGBRC(1)*CONSTS(15))/CONSTS(16)
      ALGBRC(52) =  CONSTS(2)*RATES(6)+ CONSTS(1)*RATES(5)
      ALGBRC(13) = (STATES(3) - STATES(4))/CONSTS(8)
      ALGBRC(87) = ( CONSTS(9)*CONSTS(12))/STATES(4)+CONSTS(12) ** 2.00000
      ALGBRC(92) = ( CONSTS(11)*CONSTS(14))/STATES(4)+CONSTS(14) ** 2.00000
      ALGBRC(97) = 1.00000/((1.00000+ALGBRC(87))+ALGBRC(92))
      ALGBRC(190) =  ( EXP( CONSTS(51)*ALGBRC(181))*STATES(43) ** 3.00000)*CONSTS(21)
      ALGBRC(191) =  ( EXP( (CONSTS(51) - 1.00000)*ALGBRC(181))*CONSTS(122))*STATES(4)
      ALGBRC(192) = 1.00000+ CONSTS(53)*EXP( (CONSTS(51) - 1.00000)*ALGBRC(181))
      ALGBRC(193) = ( ( CONSTS(52)*CONSTS(123))*(ALGBRC(190) - ALGBRC(191)))/( CONSTS(121)*ALGBRC(192))
      ALGBRC(106) = ( CONSTS(57)*STATES(4))/(CONSTS(58)+STATES(4))
      ALGBRC(195) =  - 2.00000*ALGBRC(193)+ALGBRC(106)
      RATES(4) =  ALGBRC(97)*(((ALGBRC(13) - ALGBRC(201)) - ALGBRC(52)) -  ( ALGBRC(195)*0.500000)*CONSTS(104))
      ALGBRC(115) =  CONSTS(112)*log(CONSTS(22)/STATES(20))
      ALGBRC(151) =  CONSTS(112)*log(CONSTS(23)/STATES(43))
      ALGBRC(166) =  ( CONSTS(32)*((((STATES(32)+STATES(33))+STATES(34))+STATES(35))+STATES(27)))*(STATES(22) - (ALGBRC(151)/3.00000+( 2.00000*ALGBRC(115))/3.00000))
      ALGBRC(194) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(181))
      ALGBRC(196) =  ( 0.0365000*CONSTS(125))*EXP( - 1.33000*ALGBRC(181))
      ALGBRC(198) = 1.00000/(ALGBRC(194)+ALGBRC(196))
      ALGBRC(103) = 1.00000+CONSTS(56)/STATES(43) ** 1.50000
      ALGBRC(200) =  ( CONSTS(54)*ALGBRC(198))*(CONSTS(124)/ALGBRC(103))
      ALGBRC(180) =  ( CONSTS(84)*(STATES(79)+STATES(70)))*(STATES(22) - ALGBRC(151))
      ALGBRC(179) =  ( CONSTS(81)*(STATES(66)+STATES(57)))*(STATES(22) - ALGBRC(151))
      ALGBRC(203) = ((ALGBRC(180)+ALGBRC(179))+ALGBRC(166)/3.00000)+ 3.00000*(ALGBRC(193)+ALGBRC(200))
      RATES(43) =  - CONSTS(104)*ALGBRC(203)
      ALGBRC(204) =  CONSTS(20)*ALGBRC(181)
      ALGBRC(187) =  0.00100000*EXP( 2.00000*ALGBRC(181)) -  CONSTS(21)*0.341000
      ALGBRC(189) = EXP( 2.00000*ALGBRC(181)) - 1.00000
      ALGBRC(205) =  ( ( CONSTS(118)*4.00000)*ALGBRC(204))*(ALGBRC(187)/ALGBRC(189))
      ALGBRC(207) = TERNRY(ALGBRC(205).GT.0.00000, 0.00000, ALGBRC(205))
      ALGBRC(208) = CONSTS(120)/(1.00000+ALGBRC(207)/CONSTS(116))
      ALGBRC(183) =  STATES(20)*EXP(ALGBRC(181)) - CONSTS(22)
      ALGBRC(185) = EXP(ALGBRC(181)) - 1.00000
      ALGBRC(209) =  ( ( ( ALGBRC(208)*STATES(21))*STATES(23))*ALGBRC(204))*(ALGBRC(183)/ALGBRC(185))
      ALGBRC(101) = 1.00000/(0.968100+EXP((STATES(22)+82.1862)/15.8864))
      ALGBRC(118) =  ( ( CONSTS(34)* ( CONSTS(22)*1.00000) ** (1.0 / 2))*ALGBRC(101))*(STATES(22) - ALGBRC(115))
      ALGBRC(121) =  ( ( CONSTS(46)*CONSTS(110))*STATES(40))*(STATES(22) - ALGBRC(115))
      ALGBRC(126) =  ( ( CONSTS(48)*STATES(42))*STATES(41))*(STATES(22) - ALGBRC(115))
      ALGBRC(133) =  ( CONSTS(59)*(1.00000/(1.00000+EXP(- (STATES(22) - 20.0000)/12.0000))))*(STATES(22) - ALGBRC(115))
      ALGBRC(140) =  ( CONSTS(60)*STATES(52))*(STATES(22) - ALGBRC(115))
      ALGBRC(112) =  (TERNRY((VOI - CONSTS(79)) -  CONSTS(80)*INT((VOI - CONSTS(79))/CONSTS(80)).LT.CONSTS(77), 1.00000, 0.00000))*CONSTS(111)
      ALGBRC(211) = (((((((ALGBRC(121)+ALGBRC(126))+ALGBRC(118))+ALGBRC(209))+( 2.00000*ALGBRC(166))/3.00000)+ALGBRC(133))+ALGBRC(112)) -  2.00000*ALGBRC(200))+ALGBRC(140)
      RATES(20) =  - CONSTS(104)*ALGBRC(211)
      ALGBRC(212) =  ( CONSTS(95)*(STATES(80)+STATES(7)))*(STATES(1) - STATES(3))
      ALGBRC(60) = ( CONSTS(10)*CONSTS(13))/STATES(1)+CONSTS(13) ** 2.00000
      ALGBRC(67) = 1.00000/(1.00000+ALGBRC(60))
      RATES(1) =  ALGBRC(67)*(ALGBRC(1) - ALGBRC(212))
      ALGBRC(206) =  ( ALGBRC(205)*STATES(23))*STATES(21)
      ALGBRC(72) = ( CONSTS(9)*CONSTS(12))/STATES(3)+CONSTS(12) ** 2.00000
      ALGBRC(77) = ( CONSTS(11)*CONSTS(14))/STATES(3)+CONSTS(14) ** 2.00000
      ALGBRC(82) = 1.00000/((1.00000+ALGBRC(72))+ALGBRC(77))
      ALGBRC(214) = ( ALGBRC(212)*CONSTS(15))/CONSTS(17) - ( ALGBRC(13)*CONSTS(18))/CONSTS(17)
      ALGBRC(210) =  ( ( ( (ALGBRC(205)/CONSTS(118))*CONSTS(31))*STATES(25))*STATES(25))*STATES(24)
      RATES(3) =  ALGBRC(82)*((ALGBRC(214) -  ALGBRC(206)*CONSTS(114)) -  ALGBRC(210)*CONSTS(114))
      ALGBRC(213) = ((((((ALGBRC(180)+ALGBRC(179))+ALGBRC(206))+ALGBRC(210))+ALGBRC(209))+ALGBRC(121))+ALGBRC(126))+ALGBRC(166)
      ALGBRC(202) = (((ALGBRC(118)+ALGBRC(193))+ALGBRC(200))+ALGBRC(140))+ALGBRC(133)
      ALGBRC(109) = ALGBRC(106)
      ALGBRC(215) = (ALGBRC(213)+ALGBRC(202))+ALGBRC(109)
      RATES(22) = - ((ALGBRC(215)+CONSTS(78))+ALGBRC(112))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(0.00336336/(0.500000+EXP(STATES(22)/- 5.53900))+ 0.00779047*EXP(STATES(22)/- 49.5104))
      ALGBRC(15) = CONSTS(29)/(1.00000+EXP((STATES(22)+28.5000)/7.80000))+(1.00000 - CONSTS(29))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(22)+75.6000)/6.20000))
      ALGBRC(16) = TERNRY(STATES(22).LT.- 60.0000, 500.000, 18.3000+ 0.00500000*EXP(- STATES(22)/6.20000))
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(22)+48.4000)/5.20000))
      ALGBRC(17) = TERNRY(STATES(22).LT.- 56.0000,  2.44000*EXP((STATES(22)+120.000)/40.0000), 1.34000+ 0.0350000*EXP(- STATES(22)/11.8000))
      ALGBRC(6) =  ( CONSTS(47)*CONSTS(35))*EXP( CONSTS(42)*STATES(22))
      ALGBRC(18) =  ( CONSTS(47)*CONSTS(36))*EXP( CONSTS(105)*STATES(22))
      ALGBRC(8) = TERNRY(STATES(22).EQ.21.0000, 0.000146000/0.0780000, ( 0.000146000*(STATES(22) - 21.0000))/(1.00000 - EXP( - 0.0780000*(STATES(22) - 21.0000))))
      ALGBRC(20) =  0.000910000*EXP( - 0.0280000*STATES(22))
      ALGBRC(9) = TERNRY(STATES(22).EQ.- 11.0000, 3.30000e-05/0.130000, ( 3.30000e-05*(STATES(22)+11.0000))/(1.00000 - EXP( - 0.130000*(STATES(22)+11.0000))))
      ALGBRC(21) =  0.000100000*EXP( - 0.0150000*STATES(22))
      ALGBRC(25) =  STATES(3)*1000.00 ** 4.00000
      ALGBRC(26) =  CONSTS(3)*STATES(5)
      ALGBRC(29) =  ALGBRC(6)*STATES(36)+ CONSTS(109)*STATES(38)
      ALGBRC(35) =  (ALGBRC(18)+CONSTS(108))*STATES(37)
      ALGBRC(39) =  STATES(3)*1000.00 ** 3.00000
      ALGBRC(40) =  CONSTS(5)*STATES(6)
      ALGBRC(19) =  ( CONSTS(47)*CONSTS(37))*EXP( CONSTS(43)*STATES(22))
      ALGBRC(28) =  ( CONSTS(47)*CONSTS(40))*EXP( CONSTS(107)*STATES(22))
      ALGBRC(48) =  ALGBRC(19)*STATES(38)+ ALGBRC(28)*STATES(39)
      ALGBRC(34) =  ( CONSTS(47)*CONSTS(38))*EXP( CONSTS(106)*STATES(22))
      ALGBRC(42) =  ( CONSTS(47)*CONSTS(39))*EXP( CONSTS(44)*STATES(22))
      ALGBRC(56) =  (ALGBRC(34)+ALGBRC(42))*STATES(40)
      ALGBRC(7) =  ( CONSTS(47)*CONSTS(41))*EXP( CONSTS(45)*STATES(22))
      ALGBRC(47) = ( ( ALGBRC(34)*ALGBRC(28))*ALGBRC(7))/( ALGBRC(19)*ALGBRC(42))
      ALGBRC(54) = ( CONSTS(108)*STATES(37)+ ALGBRC(34)*STATES(40))+ ALGBRC(47)*STATES(39)
      ALGBRC(62) =  ((ALGBRC(7)+ALGBRC(19))+CONSTS(109))*STATES(38)
      ALGBRC(55) =  ALGBRC(7)*STATES(38)+ ALGBRC(42)*STATES(40)
      ALGBRC(63) =  (ALGBRC(47)+ALGBRC(28))*STATES(39)
      ALGBRC(11) = 2.80200/( 0.210000*EXP(- STATES(22)/17.0000)+ 0.230000*EXP(- STATES(22)/150.000))
      ALGBRC(65) =  0.400000*EXP(- STATES(22)/20.3000)
      ALGBRC(12) = 3.80200/( 0.102700*EXP(- STATES(22)/17.0000)+ 0.200000*EXP(- STATES(22)/150.000))
      ALGBRC(66) =  0.191700*EXP(- STATES(22)/20.3000)
      ALGBRC(23) = 2.80200/( 0.230000*EXP(- STATES(22)/15.0000)+ 0.250000*EXP(- STATES(22)/150.000))
      ALGBRC(70) =  0.400000*EXP(- (STATES(22) - 5.00000)/20.3000)
      ALGBRC(24) = 3.80200/( 0.102700*EXP(- STATES(22)/15.0000)+ 0.230000*EXP(- STATES(22)/150.000))
      ALGBRC(71) =  0.200000*EXP(- (STATES(22) - 5.00000)/20.3000)
      ALGBRC(31) = 2.80200/( 0.250000*EXP(- STATES(22)/12.0000)+ 0.270000*EXP(- STATES(22)/150.000))
      ALGBRC(75) = ( 0.400000*EXP(- (STATES(22) - 10.0000)/20.3000))/4.50000
      ALGBRC(32) = 3.80200/( 0.102700*EXP(- STATES(22)/12.0000)+ 0.250000*EXP(- STATES(22)/150.000))
      ALGBRC(76) =  0.220000*EXP(- (STATES(22) - 10.0000)/20.3000)
      ALGBRC(44) =  ( 3.79330e-07*EXP(- STATES(22)/7.60000))*3.00000
      ALGBRC(80) = 0.00840000+ 2.00000e-05*STATES(22)
      ALGBRC(45) =  3.79330e-07*EXP(- STATES(22)/7.70000)
      ALGBRC(81) = 0.00840000+ 2.00000e-05*STATES(22)
      ALGBRC(37) = ( 9.17800*EXP(STATES(22)/29.6800))/4.50000
      ALGBRC(85) = ( ( ALGBRC(44)*ALGBRC(37))*ALGBRC(31))/( ALGBRC(80)*ALGBRC(75))
      ALGBRC(38) =  9.17800*EXP(STATES(22)/29.6800)
      ALGBRC(86) = ( ( ALGBRC(45)*ALGBRC(38))*ALGBRC(32))/( ALGBRC(81)*ALGBRC(76))
      ALGBRC(50) =  ( (ALGBRC(37)/100.000)*1.50000)*0.285000
      ALGBRC(90) = ALGBRC(44)/5.00000
      ALGBRC(51) = ALGBRC(38)/100.000
      ALGBRC(91) = ALGBRC(45)
      ALGBRC(58) =  (ALGBRC(37)/95000.0)*80.0000
      ALGBRC(95) = (ALGBRC(44)/30.0000)/10.0000
      ALGBRC(59) = ALGBRC(38)/95000.0
      ALGBRC(96) = ALGBRC(45)/50.0000
      ALGBRC(2) =  ( ( 4.00000*1.20000)*0.416000)*EXP( 0.0120000*(STATES(22) - 35.0000))
      ALGBRC(14) =  4.00000*ALGBRC(2)
      ALGBRC(116) =  ( 0.600000*0.0923300)*STATES(3)
      ALGBRC(119) = ALGBRC(116)
      ALGBRC(123) =  (ALGBRC(14)+ALGBRC(119))*STATES(10)
      ALGBRC(78) =  ( ( 4.00000*0.450000)*0.0490000)*EXP( - 0.0650000*(STATES(22) - 22.0000))
      ALGBRC(83) = ALGBRC(78)
      ALGBRC(128) =  ALGBRC(83)*STATES(11)+ CONSTS(103)*STATES(15)
      ALGBRC(46) =  CONSTS(25)*ALGBRC(2)
      ALGBRC(53) =  4.00000*ALGBRC(46)
      ALGBRC(124) =  (ALGBRC(53)+CONSTS(103))*STATES(15)
      ALGBRC(100) = ALGBRC(78)/CONSTS(26)
      ALGBRC(104) = ALGBRC(100)
      ALGBRC(130) =  ALGBRC(104)*STATES(16)+ ALGBRC(119)*STATES(10)
      ALGBRC(27) =  3.00000*ALGBRC(2)
      ALGBRC(122) =  CONSTS(25)*ALGBRC(119)
      ALGBRC(129) =  ((ALGBRC(83)+ALGBRC(27))+ALGBRC(122))*STATES(11)
      ALGBRC(88) =  2.00000*ALGBRC(78)
      ALGBRC(135) = ( ALGBRC(14)*STATES(10)+ ALGBRC(88)*STATES(12))+ CONSTS(113)*STATES(16)
      ALGBRC(61) =  3.00000*ALGBRC(46)
      ALGBRC(131) =  ((ALGBRC(104)+ALGBRC(61))+CONSTS(113))*STATES(16)
      ALGBRC(107) =  2.00000*ALGBRC(100)
      ALGBRC(137) = ( ALGBRC(53)*STATES(15)+ ALGBRC(107)*STATES(17))+ ALGBRC(122)*STATES(11)
      ALGBRC(33) =  2.00000*ALGBRC(2)
      ALGBRC(127) =  CONSTS(25)*ALGBRC(122)
      ALGBRC(136) =  ((ALGBRC(88)+ALGBRC(33))+ALGBRC(127))*STATES(12)
      ALGBRC(93) =  3.00000*ALGBRC(78)
      ALGBRC(141) = ( ALGBRC(27)*STATES(11)+ ALGBRC(93)*STATES(13))+ CONSTS(115)*STATES(17)
      ALGBRC(68) =  2.00000*ALGBRC(46)
      ALGBRC(138) =  ((ALGBRC(107)+ALGBRC(68))+CONSTS(115))*STATES(17)
      ALGBRC(110) =  3.00000*ALGBRC(100)
      ALGBRC(144) = ( ALGBRC(61)*STATES(16)+ ALGBRC(110)*STATES(18))+ ALGBRC(127)*STATES(12)
      ALGBRC(41) = ALGBRC(2)
      ALGBRC(134) =  CONSTS(25)*ALGBRC(127)
      ALGBRC(142) =  ((ALGBRC(93)+ALGBRC(41))+ALGBRC(134))*STATES(13)
      ALGBRC(98) =  4.00000*ALGBRC(78)
      ALGBRC(152) = ( ALGBRC(33)*STATES(12)+ ALGBRC(98)*STATES(14))+ CONSTS(117)*STATES(18)
      ALGBRC(73) = ALGBRC(46)
      ALGBRC(145) =  ((ALGBRC(110)+ALGBRC(73))+CONSTS(117))*STATES(18)
      ALGBRC(113) =  4.00000*ALGBRC(100)
      ALGBRC(154) = ( ALGBRC(68)*STATES(17)+ ALGBRC(113)*STATES(19))+ ALGBRC(134)*STATES(13)
      ALGBRC(143) =  CONSTS(25)*ALGBRC(134)
      ALGBRC(153) =  ((ALGBRC(98)+CONSTS(27))+ALGBRC(143))*STATES(14)
      ALGBRC(167) = ( ALGBRC(41)*STATES(13)+ CONSTS(28)*STATES(21))+ CONSTS(119)*STATES(19)
      ALGBRC(155) =  (ALGBRC(113)+CONSTS(119))*STATES(19)
      ALGBRC(168) =  ALGBRC(73)*STATES(18)+ ALGBRC(143)*STATES(14)
      ALGBRC(10) =  CONSTS(63)*EXP( CONSTS(61)*STATES(22))
      ALGBRC(22) =  4.00000*ALGBRC(10)
      ALGBRC(139) =  CONSTS(71)*EXP( CONSTS(72)*STATES(22))
      ALGBRC(146) = ALGBRC(139)
      ALGBRC(156) =  (ALGBRC(22)+ALGBRC(146))*STATES(44)
      ALGBRC(102) =  CONSTS(70)*EXP( - CONSTS(69)*STATES(22))
      ALGBRC(105) = ALGBRC(102)
      ALGBRC(74) =  CONSTS(64)*EXP( - CONSTS(62)*STATES(22))
      ALGBRC(79) = ALGBRC(74)
      ALGBRC(169) =  ALGBRC(105)*STATES(45)+ ALGBRC(79)*STATES(48)
      ALGBRC(30) =  3.00000*ALGBRC(10)
      ALGBRC(147) =  CONSTS(73)*ALGBRC(139)
      ALGBRC(157) =  ((ALGBRC(30)+ALGBRC(105))+ALGBRC(147))*STATES(45)
      ALGBRC(108) =  2.00000*ALGBRC(102)
      ALGBRC(84) = ALGBRC(74)/CONSTS(65)
      ALGBRC(170) = ( ALGBRC(108)*STATES(46)+ ALGBRC(84)*STATES(49))+ ALGBRC(22)*STATES(44)
      ALGBRC(36) =  2.00000*ALGBRC(10)
      ALGBRC(148) =  CONSTS(74)*ALGBRC(139)
      ALGBRC(158) =  ((ALGBRC(36)+ALGBRC(108))+ALGBRC(148))*STATES(46)
      ALGBRC(111) =  3.00000*ALGBRC(102)
      ALGBRC(89) = ALGBRC(74)/CONSTS(66)
      ALGBRC(171) = ( ALGBRC(111)*STATES(47)+ ALGBRC(89)*STATES(50))+ ALGBRC(30)*STATES(45)
      ALGBRC(149) =  CONSTS(75)*ALGBRC(139)
      ALGBRC(43) = ALGBRC(10)
      ALGBRC(159) =  ((ALGBRC(43)+ALGBRC(111))+ALGBRC(149))*STATES(47)
      ALGBRC(94) = ALGBRC(74)/CONSTS(67)
      ALGBRC(132) =  4.00000*ALGBRC(102)
      ALGBRC(172) = ( ALGBRC(132)*STATES(52)+ ALGBRC(94)*STATES(51))+ ALGBRC(36)*STATES(46)
      ALGBRC(49) =  ( 4.00000*CONSTS(65))*ALGBRC(10)
      ALGBRC(160) =  (ALGBRC(79)+ALGBRC(49))*STATES(48)
      ALGBRC(114) = ALGBRC(102)/CONSTS(73)
      ALGBRC(173) =  ALGBRC(146)*STATES(44)+ ALGBRC(114)*STATES(49)
      ALGBRC(57) = ( ( 3.00000*CONSTS(66))*ALGBRC(10))/CONSTS(65)
      ALGBRC(161) =  ((ALGBRC(57)+ALGBRC(84))+ALGBRC(114))*STATES(49)
      ALGBRC(117) = ( ( 2.00000*CONSTS(73))*ALGBRC(102))/CONSTS(74)
      ALGBRC(174) = ( ALGBRC(117)*STATES(50)+ ALGBRC(147)*STATES(45))+ ALGBRC(49)*STATES(48)
      ALGBRC(64) = ( ( 2.00000*CONSTS(67))*ALGBRC(10))/CONSTS(66)
      ALGBRC(162) =  ((ALGBRC(64)+ALGBRC(89))+ALGBRC(117))*STATES(50)
      ALGBRC(120) = ( ( 3.00000*CONSTS(74))*ALGBRC(102))/CONSTS(75)
      ALGBRC(175) = ( ALGBRC(120)*STATES(51)+ ALGBRC(148)*STATES(46))+ ALGBRC(57)*STATES(49)
      ALGBRC(69) = ( CONSTS(68)*ALGBRC(10))/CONSTS(67)
      ALGBRC(163) =  ((ALGBRC(69)+ALGBRC(94))+ALGBRC(120))*STATES(51)
      ALGBRC(125) = ( ( 4.00000*CONSTS(75))*ALGBRC(102))/CONSTS(76)
      ALGBRC(176) = ( ALGBRC(125)*STATES(53)+ ALGBRC(149)*STATES(47))+ ALGBRC(64)*STATES(50)
      ALGBRC(150) =  CONSTS(76)*ALGBRC(139)
      ALGBRC(165) =  (ALGBRC(132)+ALGBRC(150))*STATES(52)
      ALGBRC(99) = ALGBRC(74)/CONSTS(68)
      ALGBRC(178) =  ALGBRC(43)*STATES(47)+ ALGBRC(99)*STATES(53)
      ALGBRC(164) =  (ALGBRC(99)+ALGBRC(125))*STATES(53)
      ALGBRC(177) =  ALGBRC(150)*STATES(52)+ ALGBRC(69)*STATES(51)
      ALGBRC(181) = STATES(22)/CONSTS(112)
      ALGBRC(182) = ( 0.000171200*EXP( - 1.46500*ALGBRC(181)))/1000.00
      ALGBRC(184) = ( 26.1700*EXP( 1.46500*ALGBRC(181)))/1000.00
      ALGBRC(186) = ( 287.500*EXP( 1.24200*ALGBRC(181)))/1000.00
      ALGBRC(188) = ( 0.0402500*EXP( - 1.24200*ALGBRC(181)))/1000.00
      ALGBRC(1) = (STATES(2) - STATES(1))/CONSTS(7)
      ALGBRC(197) = STATES(4)/CONSTS(97) ** CONSTS(99)
      ALGBRC(199) = STATES(2)/CONSTS(98) ** CONSTS(100)
      ALGBRC(201) = ( CONSTS(96)*( CONSTS(101)*ALGBRC(197) -  CONSTS(102)*ALGBRC(199)))/((1.00000+ALGBRC(197))+ALGBRC(199))
      ALGBRC(52) =  CONSTS(2)*RATES(6)+ CONSTS(1)*RATES(5)
      ALGBRC(13) = (STATES(3) - STATES(4))/CONSTS(8)
      ALGBRC(87) = ( CONSTS(9)*CONSTS(12))/STATES(4)+CONSTS(12) ** 2.00000
      ALGBRC(92) = ( CONSTS(11)*CONSTS(14))/STATES(4)+CONSTS(14) ** 2.00000
      ALGBRC(97) = 1.00000/((1.00000+ALGBRC(87))+ALGBRC(92))
      ALGBRC(190) =  ( EXP( CONSTS(51)*ALGBRC(181))*STATES(43) ** 3.00000)*CONSTS(21)
      ALGBRC(191) =  ( EXP( (CONSTS(51) - 1.00000)*ALGBRC(181))*CONSTS(122))*STATES(4)
      ALGBRC(192) = 1.00000+ CONSTS(53)*EXP( (CONSTS(51) - 1.00000)*ALGBRC(181))
      ALGBRC(193) = ( ( CONSTS(52)*CONSTS(123))*(ALGBRC(190) - ALGBRC(191)))/( CONSTS(121)*ALGBRC(192))
      ALGBRC(106) = ( CONSTS(57)*STATES(4))/(CONSTS(58)+STATES(4))
      ALGBRC(195) =  - 2.00000*ALGBRC(193)+ALGBRC(106)
      ALGBRC(115) =  CONSTS(112)*log(CONSTS(22)/STATES(20))
      ALGBRC(151) =  CONSTS(112)*log(CONSTS(23)/STATES(43))
      ALGBRC(166) =  ( CONSTS(32)*((((STATES(32)+STATES(33))+STATES(34))+STATES(35))+STATES(27)))*(STATES(22) - (ALGBRC(151)/3.00000+( 2.00000*ALGBRC(115))/3.00000))
      ALGBRC(194) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(181))
      ALGBRC(196) =  ( 0.0365000*CONSTS(125))*EXP( - 1.33000*ALGBRC(181))
      ALGBRC(198) = 1.00000/(ALGBRC(194)+ALGBRC(196))
      ALGBRC(103) = 1.00000+CONSTS(56)/STATES(43) ** 1.50000
      ALGBRC(200) =  ( CONSTS(54)*ALGBRC(198))*(CONSTS(124)/ALGBRC(103))
      ALGBRC(180) =  ( CONSTS(84)*(STATES(79)+STATES(70)))*(STATES(22) - ALGBRC(151))
      ALGBRC(179) =  ( CONSTS(81)*(STATES(66)+STATES(57)))*(STATES(22) - ALGBRC(151))
      ALGBRC(203) = ((ALGBRC(180)+ALGBRC(179))+ALGBRC(166)/3.00000)+ 3.00000*(ALGBRC(193)+ALGBRC(200))
      ALGBRC(204) =  CONSTS(20)*ALGBRC(181)
      ALGBRC(187) =  0.00100000*EXP( 2.00000*ALGBRC(181)) -  CONSTS(21)*0.341000
      ALGBRC(189) = EXP( 2.00000*ALGBRC(181)) - 1.00000
      ALGBRC(205) =  ( ( CONSTS(118)*4.00000)*ALGBRC(204))*(ALGBRC(187)/ALGBRC(189))
      ALGBRC(207) = TERNRY(ALGBRC(205).GT.0.00000, 0.00000, ALGBRC(205))
      ALGBRC(208) = CONSTS(120)/(1.00000+ALGBRC(207)/CONSTS(116))
      ALGBRC(183) =  STATES(20)*EXP(ALGBRC(181)) - CONSTS(22)
      ALGBRC(185) = EXP(ALGBRC(181)) - 1.00000
      ALGBRC(209) =  ( ( ( ALGBRC(208)*STATES(21))*STATES(23))*ALGBRC(204))*(ALGBRC(183)/ALGBRC(185))
      ALGBRC(101) = 1.00000/(0.968100+EXP((STATES(22)+82.1862)/15.8864))
      ALGBRC(118) =  ( ( CONSTS(34)* ( CONSTS(22)*1.00000) ** (1.0 / 2))*ALGBRC(101))*(STATES(22) - ALGBRC(115))
      ALGBRC(121) =  ( ( CONSTS(46)*CONSTS(110))*STATES(40))*(STATES(22) - ALGBRC(115))
      ALGBRC(126) =  ( ( CONSTS(48)*STATES(42))*STATES(41))*(STATES(22) - ALGBRC(115))
      ALGBRC(133) =  ( CONSTS(59)*(1.00000/(1.00000+EXP(- (STATES(22) - 20.0000)/12.0000))))*(STATES(22) - ALGBRC(115))
      ALGBRC(140) =  ( CONSTS(60)*STATES(52))*(STATES(22) - ALGBRC(115))
      ALGBRC(112) =  (TERNRY((VOI - CONSTS(79)) -  CONSTS(80)*INT((VOI - CONSTS(79))/CONSTS(80)).LT.CONSTS(77), 1.00000, 0.00000))*CONSTS(111)
      ALGBRC(211) = (((((((ALGBRC(121)+ALGBRC(126))+ALGBRC(118))+ALGBRC(209))+( 2.00000*ALGBRC(166))/3.00000)+ALGBRC(133))+ALGBRC(112)) -  2.00000*ALGBRC(200))+ALGBRC(140)
      ALGBRC(212) =  ( CONSTS(95)*(STATES(80)+STATES(7)))*(STATES(1) - STATES(3))
      ALGBRC(60) = ( CONSTS(10)*CONSTS(13))/STATES(1)+CONSTS(13) ** 2.00000
      ALGBRC(67) = 1.00000/(1.00000+ALGBRC(60))
      ALGBRC(206) =  ( ALGBRC(205)*STATES(23))*STATES(21)
      ALGBRC(72) = ( CONSTS(9)*CONSTS(12))/STATES(3)+CONSTS(12) ** 2.00000
      ALGBRC(77) = ( CONSTS(11)*CONSTS(14))/STATES(3)+CONSTS(14) ** 2.00000
      ALGBRC(82) = 1.00000/((1.00000+ALGBRC(72))+ALGBRC(77))
      ALGBRC(214) = ( ALGBRC(212)*CONSTS(15))/CONSTS(17) - ( ALGBRC(13)*CONSTS(18))/CONSTS(17)
      ALGBRC(210) =  ( ( ( (ALGBRC(205)/CONSTS(118))*CONSTS(31))*STATES(25))*STATES(25))*STATES(24)
      ALGBRC(213) = ((((((ALGBRC(180)+ALGBRC(179))+ALGBRC(206))+ALGBRC(210))+ALGBRC(209))+ALGBRC(121))+ALGBRC(126))+ALGBRC(166)
      ALGBRC(202) = (((ALGBRC(118)+ALGBRC(193))+ALGBRC(200))+ALGBRC(140))+ALGBRC(133)
      ALGBRC(109) = ALGBRC(106)
      ALGBRC(215) = (ALGBRC(213)+ALGBRC(202))+ALGBRC(109)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END