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 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