Generated Code

C
C There are a total of 256 entries in the algebraic variable array.
C There are a total of 67 entries in each of the rate and state variable arrays.
C There are a total of 186 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C CONSTS(1) is iso in component environment (dimensionless).
C CONSTS(2) is CSQN2 in component environment (dimensionless).
C CONSTS(3) is RyR2 in component environment (dimensionless).
C CONSTS(156) is a1 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
C CONSTS(161) is a2 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
C CONSTS(4) is Faraday in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (coulomb_per_millimole).
C CONSTS(5) is Temp in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (kelvin).
C CONSTS(6) is Rgas in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (joule_per_mole_kelvin).
C CONSTS(125) is RT_over_F in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV).
C CONSTS(7) is Acap in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (cm2).
C CONSTS(136) is C in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mF).
C CONSTS(8) is Vmyo in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(9) is VJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(10) is VNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(11) is VSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C STATES(1) is Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(2) is Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(3) is Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(4) is CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(5) is CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(6) is CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C STATES(7) is V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV).
C ALGBRC(132) is INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(229) is INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(244) is INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(240) is INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(207) is IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(136) is IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(140) is IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(227) is IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(255) is ICaK in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(250) is ICa in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(221) is Ito1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(12) is i_Stim in component I_stimulus (uA_per_uF).
C ALGBRC(145) is IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(189) is IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(231) is ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(233) is IpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(61) is Jxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(43) is Jup in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(92) is Jtrpn in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM_per_ms).
C ALGBRC(52) is Jrel in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(62) is Jtr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(102) is beta_SS in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(72) is beta_JSR in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(112) is beta_i in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(256) is i_tot in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uA_per_uF).
C CONSTS(12) is Ko in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(13) is Nao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(14) is Cao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(15) is stim_period in component I_stimulus (ms).
C CONSTS(16) is stim_duration in component I_stimulus (ms).
C CONSTS(17) is stim_amplitude in component I_stimulus (uA_per_uF).
C CONSTS(18) is stim_offset in component I_stimulus (ms).
C ALGBRC(1) is past in component I_stimulus (ms).
C ALGBRC(23) is fb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C CONSTS(19) is Kfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(20) is Nfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C ALGBRC(34) is rb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C CONSTS(21) is Krb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(22) is Nrb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C CONSTS(23) is KSR in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(24) is vmaxf in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
C CONSTS(25) is vmaxr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
C CONSTS(26) is v1 in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
C STATES(8) is O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C STATES(9) is O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(27) is tautr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
C CONSTS(28) is tauxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
C CONSTS(29) is LTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(30) is HTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(31) is khtrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
C CONSTS(32) is khtrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C CONSTS(33) is kltrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
C CONSTS(34) is kltrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C CONSTS(35) is CMDNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(126) is CSQNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(36) is EGTAtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(37) is KmCMDN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(38) is KmCSQN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(39) is KmEGTA in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C ALGBRC(71) is dLTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C ALGBRC(82) is dHTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C STATES(10) is LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C STATES(11) is HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C CONSTS(40) is kaplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM4_per_ms).
C CONSTS(41) is kaminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(42) is kbplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM3_per_ms).
C CONSTS(43) is kbminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(44) is kcplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(45) is kcminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(46) is ncoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(47) is mcoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(48) is kryr in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C ALGBRC(2) is klumenC1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C ALGBRC(13) is klumenO1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C ALGBRC(24) is klumenC2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(49) is HmaxC1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(127) is HmaxO1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(50) is HmaxC2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(51) is Hmin in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(52) is H50C1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM).
C CONSTS(53) is H50O1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM).
C CONSTS(128) is H50C2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM).
C CONSTS(54) is HN in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C ALGBRC(35) is dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(44) is dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(53) is dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(63) is dO1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C STATES(12) is C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C STATES(13) is C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(55) is fL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(56) is gL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(57) is bL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C CONSTS(58) is aL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C ALGBRC(54) is C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(64) is C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(73) is C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(83) is C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(125) is C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(129) is C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(133) is C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(137) is C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(93) is CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(103) is CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(113) is CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(121) is CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(141) is CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(146) is CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(152) is CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(160) is CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(170) is C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(179) is C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(190) is C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(201) is C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(208) is C4_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(144) is CCa0_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(160) is CCa1_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(163) is CCa2_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(169) is CCa3_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(174) is CCa4_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(3) is alpha in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(14) is beta in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(25) is alpha_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(36) is beta_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(45) is gamma in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(129) is omega in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(180) is a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(191) is a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(192) is a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(202) is a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(203) is a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(209) is a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(210) is a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(215) is a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(216) is a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(222) is a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(181) is a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(193) is a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(194) is a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(204) is a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(205) is a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(211) is a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(212) is a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(217) is a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(218) is a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(223) is a2_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C STATES(14) is C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(15) is C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(16) is C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(17) is C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(18) is C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(19) is CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(20) is CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(21) is CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(22) is CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(23) is CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(24) is Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C STATES(25) is yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C CONSTS(130) is a1_Cainf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C ALGBRC(4) is yCa_inf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C ALGBRC(15) is tau_yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (ms).
C ALGBRC(5) is alpha_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(16) is beta_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(26) is alpha_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(37) is beta_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(46) is C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(55) is C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(65) is C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(74) is C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(84) is CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(94) is CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(104) is CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(114) is CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(122) is C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(126) is C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(130) is C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(134) is OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(138) is CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(143) is CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(148) is CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(154) is OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(142) is C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(147) is C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(153) is C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(161) is C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(162) is OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(149) is CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(155) is CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(163) is CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(171) is CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(172) is OIKv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C STATES(26) is C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(27) is C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(28) is C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(29) is C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(30) is OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(31) is CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(32) is CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(33) is CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(34) is CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C STATES(35) is OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C ALGBRC(156) is a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(164) is a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(165) is a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(173) is a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(174) is a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(182) is a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(183) is a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(195) is a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(184) is a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(196) is a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(157) is a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(166) is a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(167) is a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(175) is a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(176) is a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(185) is a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(186) is a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(197) is a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(187) is a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(198) is a2_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(59) is alphaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(60) is aaKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(61) is betaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(62) is baKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(63) is alphai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(64) is aiKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(65) is betai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(66) is biKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(67) is f1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(68) is f2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(69) is f3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(70) is f4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(71) is b1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(72) is b2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(73) is b3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(74) is b4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(75) is f1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(76) is f2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(77) is f3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(78) is f4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(79) is b1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(80) is b2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(81) is b3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(82) is b4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(83) is alphaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(84) is aaKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
C CONSTS(85) is betaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(86) is baKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
C CONSTS(87) is alphai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(88) is betai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(6) is alpha_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(17) is beta_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(131) is alpha_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(132) is beta_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(27) is C0Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(38) is C1Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(47) is C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(56) is C3Kv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(28) is CI0Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(39) is CI1Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(48) is CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(57) is CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(66) is C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(75) is C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(85) is C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(95) is OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(67) is CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(76) is CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(86) is CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(96) is OIKv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(145) is C0Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(146) is C1Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(147) is C2Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(148) is C3Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(149) is OKv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(150) is CI0Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(151) is CI1Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(152) is CI2Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(153) is CI3Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(154) is OIKv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(77) is a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(87) is a2_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(36) is C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(88) is a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(97) is a2_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(37) is C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(98) is a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(105) is a2_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(38) is C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(106) is a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(115) is a2_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(39) is C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(107) is a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(116) is a2_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(40) is OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(78) is a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(89) is a2_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(41) is CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(90) is a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(99) is a2_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(42) is CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(100) is a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(108) is a2_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(43) is CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(109) is a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(117) is a2_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(44) is CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(110) is a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(118) is a2_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C STATES(45) is OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C ALGBRC(120) is ENa in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(124) is EK in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(18) is EKs in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(128) is ECa in component COMPUTE_REVERSAL_POTENTIALS (mV).
C CONSTS(155) is a1 in component COMPUTE_REVERSAL_POTENTIALS (mM).
C ALGBRC(7) is a2 in component COMPUTE_REVERSAL_POTENTIALS (mM).
C CONSTS(89) is GKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(90) is GKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(91) is GK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(92) is GNa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(93) is GNab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(94) is KvScale in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C CONSTS(95) is Kv43Frac in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C CONSTS(133) is GKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(134) is PKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (litre_per_farad_second).
C STATES(46) is na6 in component INa (dimensionless).
C STATES(47) is na7 in component INa (dimensionless).
C STATES(48) is OHerg in component IKr (dimensionless).
C STATES(49) is O1ks in component IKs (dimensionless).
C STATES(50) is O2ks in component IKs (dimensionless).
C CONSTS(135) is fKo in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(151) is VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(159) is VFsq_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (coulomb_per_millimole).
C ALGBRC(169) is a1_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
C ALGBRC(178) is a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(200) is a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
C ALGBRC(225) is K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(214) is IKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C CONSTS(96) is kNaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(97) is KmNa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(98) is KmCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(99) is ksat in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(100) is eta in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(101) is INaKmax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(102) is KmNai in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(103) is KmKo in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(104) is IpCamax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(105) is KmpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(106) is GCab in component COMPUTE_INaK_INaCa_ICab_IpCa (mS_per_uF).
C ALGBRC(235) is VF_over_RT in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(164) is sigma in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(236) is a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(237) is a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(238) is fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(168) is a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(239) is a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(241) is a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
C ALGBRC(242) is a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
C ALGBRC(243) is a3_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(173) is a4_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(176) is a5_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM3).
C ALGBRC(252) is PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C ALGBRC(245) is VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(246) is VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole).
C ALGBRC(247) is a1_Ca in component COMPUTE_ICa_ICaK (mM).
C ALGBRC(248) is a2_Ca in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(253) is a1_K in component COMPUTE_ICa_ICaK (mM).
C ALGBRC(254) is a2_K in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(249) is ICamax in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(251) is Icabar in component COMPUTE_ICa_ICaK (uA_per_uF).
C CONSTS(157) is PCa in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C CONSTS(137) is PK in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C CONSTS(107) is ICahalf in component COMPUTE_ICa_ICaK (uA_per_uF).
C CONSTS(108) is Pscale in component COMPUTE_ICa_ICaK (dimensionless).
C CONSTS(162) is Temp_Scale in component INa (dimensionless).
C ALGBRC(8) is alpha1 in component INa (per_ms).
C ALGBRC(19) is beta1 in component INa (per_ms).
C ALGBRC(29) is gamma1 in component INa (per_ms).
C ALGBRC(40) is Delta1 in component INa (per_ms).
C ALGBRC(49) is On in component INa (per_ms).
C ALGBRC(58) is Of in component INa (per_ms).
C ALGBRC(68) is GammaGamma in component INa (per_ms).
C ALGBRC(79) is DeltaDelta in component INa (per_ms).
C CONSTS(166) is epsilon in component INa (per_ms).
C CONSTS(167) is omega_na in component INa (per_ms).
C ALGBRC(91) is rho in component INa (per_ms).
C ALGBRC(101) is mu in component INa (per_ms).
C CONSTS(165) is Cn in component INa (per_ms).
C CONSTS(170) is Cf in component INa (per_ms).
C CONSTS(175) is parameter_a in component INa (dimensionless).
C ALGBRC(111) is k12 in component INa (per_ms).
C ALGBRC(119) is k23 in component INa (per_ms).
C ALGBRC(123) is k34 in component INa (per_ms).
C ALGBRC(127) is k45 in component INa (per_ms).
C ALGBRC(131) is k56 in component INa (per_ms).
C CONSTS(171) is k67 in component INa (per_ms).
C ALGBRC(135) is k89 in component INa (per_ms).
C ALGBRC(139) is k910 in component INa (per_ms).
C ALGBRC(144) is k1011 in component INa (per_ms).
C ALGBRC(150) is k1112 in component INa (per_ms).
C ALGBRC(158) is k1213 in component INa (per_ms).
C ALGBRC(168) is k57 in component INa (per_ms).
C ALGBRC(177) is k21 in component INa (per_ms).
C ALGBRC(188) is k32 in component INa (per_ms).
C ALGBRC(199) is k43 in component INa (per_ms).
C ALGBRC(206) is k54 in component INa (per_ms).
C ALGBRC(213) is k65 in component INa (per_ms).
C CONSTS(172) is k76 in component INa (per_ms).
C ALGBRC(219) is k98 in component INa (per_ms).
C ALGBRC(224) is k109 in component INa (per_ms).
C ALGBRC(226) is k1110 in component INa (per_ms).
C ALGBRC(228) is k1211 in component INa (per_ms).
C ALGBRC(230) is k1312 in component INa (per_ms).
C ALGBRC(220) is k75 in component INa (per_ms).
C CONSTS(177) is k81 in component INa (per_ms).
C CONSTS(179) is k92 in component INa (per_ms).
C CONSTS(181) is k103 in component INa (per_ms).
C CONSTS(183) is k114 in component INa (per_ms).
C CONSTS(185) is k125 in component INa (per_ms).
C ALGBRC(232) is k136 in component INa (per_ms).
C CONSTS(178) is k18 in component INa (per_ms).
C CONSTS(180) is k29 in component INa (per_ms).
C CONSTS(182) is k310 in component INa (per_ms).
C CONSTS(184) is k411 in component INa (per_ms).
C CONSTS(186) is k512 in component INa (per_ms).
C ALGBRC(234) is k613 in component INa (per_ms).
C STATES(51) is na1 in component INa (dimensionless).
C STATES(52) is na2 in component INa (dimensionless).
C STATES(53) is na3 in component INa (dimensionless).
C STATES(54) is na4 in component INa (dimensionless).
C STATES(55) is na5 in component INa (dimensionless).
C STATES(56) is na8 in component INa (dimensionless).
C STATES(57) is na9 in component INa (dimensionless).
C STATES(58) is na10 in component INa (dimensionless).
C STATES(59) is na11 in component INa (dimensionless).
C STATES(60) is na12 in component INa (dimensionless).
C STATES(61) is na13 in component INa (dimensionless).
C CONSTS(109) is TNa in component INa (kelvin).
C CONSTS(138) is KToverH in component INa (per_ms).
C CONSTS(159) is FoverRT in component INa (per_mV).
C CONSTS(158) is RTNa in component INa (joule_per_mole).
C CONSTS(139) is RTNaF in component INa (mV).
C CONSTS(110) is T_Const_HERG in component IKr (dimensionless).
C CONSTS(111) is A0_HERG in component IKr (per_ms).
C CONSTS(112) is B0_HERG in component IKr (per_mV).
C CONSTS(113) is A1_HERG in component IKr (per_ms).
C CONSTS(114) is B1_HERG in component IKr (per_mV).
C CONSTS(115) is A2_HERG in component IKr (per_ms).
C CONSTS(116) is B2_HERG in component IKr (per_mV).
C CONSTS(117) is A3_HERG in component IKr (per_ms).
C CONSTS(118) is B3_HERG in component IKr (per_mV).
C CONSTS(119) is A4_HERG in component IKr (per_ms).
C CONSTS(120) is B4_HERG in component IKr (per_mV).
C CONSTS(121) is A5_HERG in component IKr (per_ms).
C CONSTS(122) is B5_HERG in component IKr (per_mV).
C CONSTS(123) is A6_HERG in component IKr (per_ms).
C CONSTS(124) is B6_HERG in component IKr (per_mV).
C CONSTS(140) is C2H_to_C3H in component IKr (per_ms).
C CONSTS(141) is C3H_to_C2H in component IKr (per_ms).
C ALGBRC(9) is C1H_to_C2H in component IKr (per_ms).
C ALGBRC(20) is C2H_to_C1H in component IKr (per_ms).
C ALGBRC(10) is C3H_to_OH in component IKr (per_ms).
C ALGBRC(21) is OH_to_C3H in component IKr (per_ms).
C ALGBRC(30) is OH_to_IH in component IKr (per_ms).
C ALGBRC(41) is IH_to_OH in component IKr (per_ms).
C ALGBRC(50) is C3H_to_IH in component IKr (per_ms).
C ALGBRC(59) is IH_to_C3H in component IKr (per_ms).
C STATES(62) is C1Herg in component IKr (dimensionless).
C ALGBRC(31) is a1_C2 in component IKr (per_ms).
C ALGBRC(42) is a2_C2 in component IKr (per_ms).
C STATES(63) is C2Herg in component IKr (dimensionless).
C ALGBRC(69) is a1_C3 in component IKr (per_ms).
C ALGBRC(80) is a2_C3 in component IKr (per_ms).
C STATES(64) is C3Herg in component IKr (dimensionless).
C ALGBRC(51) is a1_O in component IKr (per_ms).
C ALGBRC(60) is a2_O in component IKr (per_ms).
C ALGBRC(70) is a1_I in component IKr (per_ms).
C ALGBRC(81) is a2_I in component IKr (per_ms).
C STATES(65) is IHerg in component IKr (dimensionless).
C CONSTS(142) is C0ks_C1ks in component IKs (per_ms).
C CONSTS(143) is C1ks_O1ks in component IKs (per_ms).
C ALGBRC(11) is O1ks_O2ks in component IKs (per_ms).
C ALGBRC(22) is O1ks_C1ks in component IKs (per_ms).
C ALGBRC(32) is O2ks_O1ks in component IKs (per_ms).
C ALGBRC(33) is C1ks_C0ks in component IKs (per_ms).
C STATES(66) is C0ks in component IKs (dimensionless).
C STATES(67) is C1ks in component IKs (dimensionless).
C RATES(1) is d/dt Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(2) is d/dt Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(3) is d/dt Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(4) is d/dt CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(5) is d/dt CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(6) is d/dt CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C RATES(7) is d/dt V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV).
C RATES(10) is d/dt LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C RATES(11) is d/dt HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C RATES(12) is d/dt C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C RATES(9) is d/dt O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C RATES(13) is d/dt C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C RATES(8) is d/dt O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C RATES(14) is d/dt C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(15) is d/dt C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(16) is d/dt C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(17) is d/dt C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(18) is d/dt C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(24) is d/dt Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(19) is d/dt CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(20) is d/dt CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(21) is d/dt CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(22) is d/dt CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(23) is d/dt CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(25) is d/dt yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C RATES(26) is d/dt C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(27) is d/dt C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(28) is d/dt C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(29) is d/dt C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(30) is d/dt OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(31) is d/dt CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(32) is d/dt CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(33) is d/dt CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(34) is d/dt CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(35) is d/dt OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C RATES(36) is d/dt C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(37) is d/dt C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(38) is d/dt C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(39) is d/dt C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(40) is d/dt OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(41) is d/dt CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(42) is d/dt CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(43) is d/dt CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(44) is d/dt CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(45) is d/dt OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C RATES(51) is d/dt na1 in component INa (dimensionless).
C RATES(52) is d/dt na2 in component INa (dimensionless).
C RATES(53) is d/dt na3 in component INa (dimensionless).
C RATES(54) is d/dt na4 in component INa (dimensionless).
C RATES(55) is d/dt na5 in component INa (dimensionless).
C RATES(46) is d/dt na6 in component INa (dimensionless).
C RATES(47) is d/dt na7 in component INa (dimensionless).
C RATES(56) is d/dt na8 in component INa (dimensionless).
C RATES(57) is d/dt na9 in component INa (dimensionless).
C RATES(58) is d/dt na10 in component INa (dimensionless).
C RATES(59) is d/dt na11 in component INa (dimensionless).
C RATES(60) is d/dt na12 in component INa (dimensionless).
C RATES(61) is d/dt na13 in component INa (dimensionless).
C RATES(62) is d/dt C1Herg in component IKr (dimensionless).
C RATES(63) is d/dt C2Herg in component IKr (dimensionless).
C RATES(64) is d/dt C3Herg in component IKr (dimensionless).
C RATES(48) is d/dt OHerg in component IKr (dimensionless).
C RATES(65) is d/dt IHerg in component IKr (dimensionless).
C RATES(66) is d/dt C0ks in component IKs (dimensionless).
C RATES(67) is d/dt C1ks in component IKs (dimensionless).
C RATES(49) is d/dt O1ks in component IKs (dimensionless).
C RATES(50) is d/dt O2ks in component IKs (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0
      CONSTS(2) = 0
      CONSTS(3) = 0
      CONSTS(4) = 96.5
      CONSTS(5) = 310
      CONSTS(6) = 8.315
      CONSTS(7) = 0.0001534
      CONSTS(8) = 2.584e-5
      CONSTS(9) = 1.6e-7
      CONSTS(10) = 2.1e-6
      CONSTS(11) = 1.2e-9
      STATES(1) = 9.85573275838928
      STATES(2) = 125.427082712469
      STATES(3) = 0.000363968672182656
      STATES(4) = 0.000506604278037024
      STATES(5) = 0.421936980515042
      STATES(6) = 0.423551621440241
      STATES(7) = -86.7261544519706
      CONSTS(12) = 4
      CONSTS(13) = 138
      CONSTS(14) = 2
      CONSTS(15) = 500
      CONSTS(16) = 3
      CONSTS(17) = -15
      CONSTS(18) = 100
      CONSTS(19) = 0.000168
      CONSTS(20) = 1.2
      CONSTS(21) = 3.29
      CONSTS(22) = 1
      CONSTS(23) = 1.2
      CONSTS(24) = 7.48e-5
      CONSTS(25) = 0.000318
      CONSTS(26) = 1.8
      STATES(8) = 0.00113684728532807
      STATES(9) = 3.11350788541838e-7
      CONSTS(27) = 0.5747
      CONSTS(28) = 26.7
      CONSTS(29) = 0.07
      CONSTS(30) = 0.14
      CONSTS(31) = 20
      CONSTS(32) = 6.6e-5
      CONSTS(33) = 40
      CONSTS(34) = 0.04
      CONSTS(35) = 0.05
      CONSTS(36) = 0
      CONSTS(37) = 0.00238
      CONSTS(38) = 0.8
      CONSTS(39) = 0.00015
      STATES(10) = 0.280466039150394
      STATES(11) = 0.99347761599363
      CONSTS(40) = 0.01215
      CONSTS(41) = 0.576
      CONSTS(42) = 0.00405
      CONSTS(43) = 1.93
      CONSTS(44) = 0.1
      CONSTS(45) = 0.0008
      CONSTS(46) = 4
      CONSTS(47) = 3
      CONSTS(48) = 1
      CONSTS(49) = 0.5
      CONSTS(50) = 0.5
      CONSTS(51) = 5
      CONSTS(52) = 1
      CONSTS(53) = 1
      CONSTS(54) = 2.5
      STATES(12) = 0.132070890861418
      STATES(13) = 0.866791951404883
      CONSTS(55) = 0.3
      CONSTS(56) = 4
      CONSTS(57) = 2
      CONSTS(58) = 2
      STATES(14) = 0.465679150104636
      STATES(15) = 0.00834457719966281
      STATES(16) = 5.60736209083906e-5
      STATES(17) = 1.6747092904465e-7
      STATES(18) = 1.87571666668874e-10
      STATES(19) = 0.489846779190386
      STATES(20) = 0.035111008610982
      STATES(21) = 0.000943745917866092
      STATES(22) = 1.12741202215634e-5
      STATES(23) = 5.05056944609524e-8
      STATES(24) = 1.40806027419488e-11
      STATES(25) = 0.995434385054729
      STATES(26) = 0.908189132330738
      STATES(27) = 0.0343385704915328
      STATES(28) = 0.000487654173162347
      STATES(29) = 3.11550715247964e-6
      STATES(30) = 7.42911977991342e-9
      STATES(31) = 0.0349937004781012
      STATES(32) = 0.0171163265867255
      STATES(33) = 0.00428471710061031
      STATES(34) = 0.000564724236640674
      STATES(35) = 2.19603439704397e-5
      CONSTS(59) = 0.543708
      CONSTS(60) = 0.028983
      CONSTS(61) = 0.080185
      CONSTS(62) = 0.0468437
      CONSTS(63) = 0.0498424
      CONSTS(64) = 0.000373016
      CONSTS(65) = 0.000819482
      CONSTS(66) = 5.374e-8
      CONSTS(67) = 1.8936
      CONSTS(68) = 14.224647456
      CONSTS(69) = 158.574378389
      CONSTS(70) = 142.936645351
      CONSTS(71) = 6.77348
      CONSTS(72) = 15.6212705152
      CONSTS(73) = 28.7532603313
      CONSTS(74) = 524.576206679
      CONSTS(75) = 0.52465073996
      CONSTS(76) = 17.51885408639
      CONSTS(77) = 938.58764534556
      CONSTS(78) = 54749.194733326
      CONSTS(79) = 1.00947847105
      CONSTS(80) = 1.17100540567
      CONSTS(81) = 0.63902768758
      CONSTS(82) = 2.12035379095
      CONSTS(83) = 1.84002414554
      CONSTS(84) = 0.00768548031
      CONSTS(85) = 0.0108174834
      CONSTS(86) = 0.07793378174
      CONSTS(87) = 0.00305767916
      CONSTS(88) = 2.44936e-6
      STATES(36) = 0.149374350989705
      STATES(37) = 0.0606794865684932
      STATES(38) = 0.00930314185504921
      STATES(39) = 0.000676403999474111
      STATES(40) = 3.85187206387239e-5
      STATES(41) = 0.0442722560882536
      STATES(42) = 0.00952432663172288
      STATES(43) = 0.0567396669678271
      STATES(44) = 0.113122845136053
      STATES(45) = 0.556269044084734
      CONSTS(89) = 0.0186
      CONSTS(90) = 0.0035
      CONSTS(91) = 0.125305126118808
      CONSTS(92) = 56.32
      CONSTS(93) = 0.001
      CONSTS(94) = 0.872
      CONSTS(95) = 0.889
      STATES(46) = 1.02118700961583e-7
      STATES(47) = 1.93499158844817e-8
      STATES(48) = 0.00120284688677794
      STATES(49) = 5.65460174551007e-7
      STATES(50) = 0.0258818770122187
      CONSTS(96) = 0.44
      CONSTS(97) = 87.5
      CONSTS(98) = 1.38
      CONSTS(99) = 0.2
      CONSTS(100) = 0.35
      CONSTS(101) = 2.387
      CONSTS(102) = 20
      CONSTS(103) = 1.5
      CONSTS(104) = 0.05
      CONSTS(105) = 0.0005
      CONSTS(106) = 7.684e-5
      CONSTS(107) = -0.265
      CONSTS(108) = 7
      STATES(51) = 0.111284526171411
      STATES(52) = 0.0481019786429977
      STATES(53) = 0.00779692701457915
      STATES(54) = 0.000561699600929369
      STATES(55) = 1.51746424723121e-5
      STATES(56) = 0.368582741846592
      STATES(57) = 0.312463212648791
      STATES(58) = 0.0993398770493615
      STATES(59) = 0.0140431688972267
      STATES(60) = 0.000750073829883749
      STATES(61) = 0.0370604970714329
      CONSTS(109) = 294.16
      CONSTS(110) = 5.320000001
      CONSTS(111) = 0.017147641733086
      CONSTS(112) = 0.03304608038835
      CONSTS(113) = 0.03969328381141
      CONSTS(114) = -0.0430605416398
      CONSTS(115) = 0.02057448605977
      CONSTS(116) = 0.02617412715118
      CONSTS(117) = 0.00134366604423
      CONSTS(118) = -0.02691385498399
      CONSTS(119) = 0.10666316491288
      CONSTS(120) = 0.00568908859717
      CONSTS(121) = 0.00646393910049
      CONSTS(122) = -0.04536642959543
      CONSTS(123) = 8.039374403e-5
      CONSTS(124) = 6.9808924e-7
      STATES(62) = 0.994948338598163
      STATES(63) = 0.000595653663190548
      STATES(64) = 0.000228183228829573
      STATES(65) = 0.000243789721526602
      STATES(66) = 0.938064990549233
      STATES(67) = 0.0360525668093578
      CONSTS(125) = ( CONSTS(6)*CONSTS(5))/CONSTS(4)
      CONSTS(126) = TERNRY(CONSTS(3).EQ.1.00000, 15.0000,  0.500000*15.0000)
      CONSTS(127) = TERNRY(CONSTS(2).EQ.1.00000, 0.200000, 0.500000)
      CONSTS(128) = TERNRY(CONSTS(2).EQ.1.00000, 0.500000, 1.00000)
      CONSTS(129) =  0.250000*0.0100000
      CONSTS(130) = 0.820000
      CONSTS(131) = CONSTS(87)
      CONSTS(132) = CONSTS(88)
      CONSTS(133) =  CONSTS(95)*CONSTS(94)*0.100000
      CONSTS(134) =  (1.00000 - CONSTS(95))*CONSTS(94)*4.29860e-07
      CONSTS(135) =  (CONSTS(12)/4.00000) ** (1.0 / 2)
      CONSTS(136) =  CONSTS(7)*0.00100000
      CONSTS(137) =  CONSTS(108)*4.57400e-07
      CONSTS(138) = ( 1.38100e-23*CONSTS(109))/6.62600e-31
      CONSTS(139) = ( CONSTS(6)*CONSTS(109))/CONSTS(4)
      CONSTS(140) =  CONSTS(110)*0.0260836
      CONSTS(141) =  CONSTS(110)*0.148330
      CONSTS(142) = 0.00795601
      CONSTS(143) = 0.0396672
      CONSTS(144) = CONSTS(129)
      CONSTS(145) = CONSTS(132)
      CONSTS(146) =  CONSTS(75)*CONSTS(132)
      CONSTS(147) =  CONSTS(76)*CONSTS(132)
      CONSTS(148) =  CONSTS(77)*CONSTS(132)
      CONSTS(149) =  CONSTS(78)*CONSTS(132)
      CONSTS(150) = CONSTS(131)
      CONSTS(151) = CONSTS(131)/CONSTS(79)
      CONSTS(152) = CONSTS(131)/CONSTS(80)
      CONSTS(153) = CONSTS(131)/CONSTS(81)
      CONSTS(154) = CONSTS(131)/CONSTS(82)
      CONSTS(155) = CONSTS(12)+ 0.0183300*CONSTS(13)
      CONSTS(156) = CONSTS(136)/( CONSTS(8)*CONSTS(4))
      CONSTS(157) = TERNRY(CONSTS(1).EQ.0.00000,  CONSTS(108)*0.000246900,  1.50000*CONSTS(108)*0.000246900)
      CONSTS(158) =  CONSTS(6)*CONSTS(109)
      CONSTS(159) = 1.00000/CONSTS(125)
      CONSTS(160) = CONSTS(144)/CONSTS(57)
      CONSTS(161) = CONSTS(136)/( 2.00000*CONSTS(11)*CONSTS(4))
      CONSTS(162) = 1.38862
      CONSTS(163) = CONSTS(160)/CONSTS(57)
      CONSTS(164) = (EXP(CONSTS(13)/67.3000) - 1.00000)/7.00000
      CONSTS(165) =  CONSTS(162)*CONSTS(138)*EXP(- 287913./CONSTS(158)+786.217/CONSTS(6))
      CONSTS(166) =  CONSTS(162)*CONSTS(138)*EXP(- 85800.4/CONSTS(158)+70.0780/CONSTS(6))
      CONSTS(167) =  CONSTS(162)*CONSTS(138)*EXP(- 121955./CONSTS(158)+225.175/CONSTS(6))
      CONSTS(168) = CONSTS(12)/(CONSTS(12)+CONSTS(103))
      CONSTS(169) = CONSTS(163)/CONSTS(57)
      CONSTS(170) =  CONSTS(162)*CONSTS(138)*EXP(- 59565.2/CONSTS(158)+0.00711000/CONSTS(6))
      CONSTS(171) = CONSTS(166)
      CONSTS(172) = CONSTS(167)
      CONSTS(173) = CONSTS(98)+CONSTS(14)
      CONSTS(174) = CONSTS(169)/CONSTS(57)
      CONSTS(175) = 1.40043
      CONSTS(176) = (CONSTS(97) ** 3.00000+CONSTS(13) ** 3.00000)/5000.00
      CONSTS(177) = CONSTS(170)
      CONSTS(178) = CONSTS(165)
      CONSTS(179) = CONSTS(177)/CONSTS(175)
      CONSTS(180) =  CONSTS(178)*CONSTS(175)
      CONSTS(181) = CONSTS(179)/CONSTS(175)
      CONSTS(182) =  CONSTS(180)*CONSTS(175)
      CONSTS(183) = CONSTS(181)/CONSTS(175)
      CONSTS(184) =  CONSTS(182)*CONSTS(175)
      CONSTS(185) = CONSTS(183)/CONSTS(175)
      CONSTS(186) =  CONSTS(184)*CONSTS(175)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(24) =  CONSTS(55)*STATES(18) -  CONSTS(56)*STATES(24)
      ALGBRC(4) = (CONSTS(130)/(1.00000+EXP((STATES(7)+28.5000)/7.80000))+1.00000) - CONSTS(130)
      ALGBRC(15) = 1.00000/(0.00336336/(0.500000+EXP(STATES(7)/- 5.53900))+ 0.00779047*EXP(STATES(7)/- 49.5104))
      RATES(25) = (ALGBRC(4) - STATES(25))/ALGBRC(15)
      ALGBRC(9) =  CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*STATES(7))
      ALGBRC(20) =  CONSTS(110)*CONSTS(113)*EXP( CONSTS(114)*STATES(7))
      RATES(62) =  ALGBRC(20)*STATES(63) -  ALGBRC(9)*STATES(62)
      ALGBRC(33) =  0.216256*EXP( - 1.88912e-05*STATES(7))
      RATES(66) =  - CONSTS(142)*STATES(66)+ ALGBRC(33)*STATES(67)
      ALGBRC(22) =  0.00700807*EXP( - 0.149998*STATES(7))
      RATES(67) = ( CONSTS(142)*STATES(66) -  (ALGBRC(33)+CONSTS(143))*STATES(67))+ ALGBRC(22)*STATES(49)
      ALGBRC(11) =  0.00767254*EXP( 0.0866295*STATES(7))
      ALGBRC(32) =  0.00379738*EXP( - 0.0142567*STATES(7))
      RATES(49) = ( CONSTS(143)*STATES(67) -  (ALGBRC(22)+ALGBRC(11))*STATES(49))+ ALGBRC(32)*STATES(50)
      RATES(50) =  ALGBRC(11)*STATES(49) -  ALGBRC(32)*STATES(50)
      ALGBRC(2) = CONSTS(49) - (CONSTS(49) - CONSTS(51))/(1.00000+CONSTS(52)/STATES(5) ** CONSTS(54))
      ALGBRC(35) =  - CONSTS(40)*ALGBRC(2)* STATES(4)*1000.00 **  CONSTS(46)*CONSTS(48)*STATES(12)+ CONSTS(41)*STATES(8)
      RATES(12) = ALGBRC(35)
      ALGBRC(31) =  ALGBRC(9)*STATES(62)+ CONSTS(141)*STATES(64)
      ALGBRC(42) =  (ALGBRC(20)+CONSTS(140))*STATES(63)
      RATES(63) = ALGBRC(31) - ALGBRC(42)
      ALGBRC(44) =  CONSTS(42)* STATES(4)*1000.00 **  CONSTS(47)*CONSTS(48)*STATES(8) -  CONSTS(43)*STATES(9)
      RATES(9) = ALGBRC(44)
      ALGBRC(13) = CONSTS(127) - (CONSTS(127) - CONSTS(51))/(1.00000+CONSTS(53)/STATES(5) ** CONSTS(54))
      ALGBRC(24) = CONSTS(50) - (CONSTS(50) - CONSTS(51))/(1.00000+CONSTS(128)/STATES(5) ** CONSTS(54))
      ALGBRC(53) =  (CONSTS(44)/ALGBRC(13))*STATES(8) -  CONSTS(45)*ALGBRC(24)*STATES(13)
      RATES(13) = ALGBRC(53)
      ALGBRC(10) =  CONSTS(110)*CONSTS(115)*EXP( CONSTS(116)*STATES(7))
      ALGBRC(41) =  CONSTS(110)*CONSTS(121)*EXP( CONSTS(122)*STATES(7))
      ALGBRC(51) =  ALGBRC(10)*STATES(64)+ ALGBRC(41)*STATES(65)
      ALGBRC(21) =  CONSTS(110)*CONSTS(117)*EXP( CONSTS(118)*STATES(7))
      ALGBRC(30) =  CONSTS(110)*CONSTS(119)*EXP( CONSTS(120)*STATES(7))
      ALGBRC(60) =  (ALGBRC(21)+ALGBRC(30))*STATES(48)
      RATES(48) = ALGBRC(51) - ALGBRC(60)
      ALGBRC(23) = STATES(3)/CONSTS(19) ** CONSTS(20)
      ALGBRC(34) = STATES(6)/CONSTS(21) ** CONSTS(22)
      ALGBRC(43) = TERNRY(CONSTS(1).EQ.0.00000, ( CONSTS(23)*( CONSTS(24)*ALGBRC(23) -  CONSTS(25)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34)), ( 1.50000*CONSTS(23)*( CONSTS(24)*ALGBRC(23) -  CONSTS(25)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34)))
      ALGBRC(62) = (STATES(6) - STATES(5))/CONSTS(27)
      RATES(6) = ( ALGBRC(43)*CONSTS(8))/CONSTS(10) - ( ALGBRC(62)*CONSTS(9))/CONSTS(10)
      ALGBRC(63) = - (ALGBRC(35)+ALGBRC(44)+ALGBRC(53))
      RATES(8) = ALGBRC(63)
      ALGBRC(52) =  CONSTS(26)*(STATES(8)+STATES(9))*(STATES(5) - STATES(4))
      ALGBRC(72) = 1.00000/(1.00000+( CONSTS(126)*CONSTS(38))/STATES(5)+CONSTS(38) ** 2.00000)
      RATES(5) =  ALGBRC(72)*(ALGBRC(62) - ALGBRC(52))
      ALGBRC(71) =  CONSTS(33)*STATES(3)*(1.00000 - STATES(10)) -  CONSTS(34)*STATES(10)
      RATES(10) = ALGBRC(71)
      ALGBRC(50) =  CONSTS(110)*CONSTS(123)*EXP( CONSTS(124)*STATES(7))
      ALGBRC(59) = ( ALGBRC(21)*ALGBRC(41)*ALGBRC(50))/( ALGBRC(10)*ALGBRC(30))
      ALGBRC(69) =  CONSTS(140)*STATES(63)+ ALGBRC(21)*STATES(48)+ ALGBRC(59)*STATES(65)
      ALGBRC(80) =  (ALGBRC(50)+ALGBRC(10)+CONSTS(141))*STATES(64)
      RATES(64) = ALGBRC(69) - ALGBRC(80)
      ALGBRC(70) =  ALGBRC(50)*STATES(64)+ ALGBRC(30)*STATES(48)
      ALGBRC(81) =  (ALGBRC(59)+ALGBRC(41))*STATES(65)
      RATES(65) = ALGBRC(70) - ALGBRC(81)
      ALGBRC(82) =  CONSTS(31)*STATES(3)*(1.00000 - STATES(11)) -  CONSTS(32)*STATES(11)
      RATES(11) = ALGBRC(82)
      ALGBRC(6) =  CONSTS(83)*EXP( CONSTS(84)*STATES(7))
      ALGBRC(27) =  4.00000*ALGBRC(6)
      ALGBRC(77) =  (ALGBRC(27)+CONSTS(145))*STATES(36)
      ALGBRC(17) =  CONSTS(85)*EXP( - CONSTS(86)*STATES(7))
      ALGBRC(66) = ALGBRC(17)
      ALGBRC(87) =  ALGBRC(66)*STATES(37)+ CONSTS(150)*STATES(41)
      RATES(36) = ALGBRC(87) - ALGBRC(77)
      ALGBRC(28) =  4.00000*CONSTS(79)*ALGBRC(6)
      ALGBRC(78) =  (CONSTS(150)+ALGBRC(28))*STATES(41)
      ALGBRC(67) = ALGBRC(17)/CONSTS(75)
      ALGBRC(89) =  CONSTS(145)*STATES(36)+ ALGBRC(67)*STATES(42)
      RATES(41) = ALGBRC(89) - ALGBRC(78)
      ALGBRC(38) =  3.00000*ALGBRC(6)
      ALGBRC(88) =  (ALGBRC(38)+ALGBRC(66)+CONSTS(146))*STATES(37)
      ALGBRC(75) =  2.00000*ALGBRC(17)
      ALGBRC(97) =  ALGBRC(75)*STATES(38)+ CONSTS(151)*STATES(42)+ ALGBRC(27)*STATES(36)
      RATES(37) = ALGBRC(97) - ALGBRC(88)
      ALGBRC(39) = ( 3.00000*CONSTS(80)*ALGBRC(6))/CONSTS(79)
      ALGBRC(90) =  (ALGBRC(39)+CONSTS(151)+ALGBRC(67))*STATES(42)
      ALGBRC(76) = ( 2.00000*CONSTS(75)*ALGBRC(17))/CONSTS(76)
      ALGBRC(99) =  ALGBRC(76)*STATES(43)+ CONSTS(146)*STATES(37)+ ALGBRC(28)*STATES(41)
      RATES(42) = ALGBRC(99) - ALGBRC(90)
      ALGBRC(47) =  2.00000*ALGBRC(6)
      ALGBRC(98) =  (ALGBRC(47)+ALGBRC(75)+CONSTS(147))*STATES(38)
      ALGBRC(85) =  3.00000*ALGBRC(17)
      ALGBRC(105) =  ALGBRC(85)*STATES(39)+ CONSTS(152)*STATES(43)+ ALGBRC(38)*STATES(37)
      RATES(38) = ALGBRC(105) - ALGBRC(98)
      ALGBRC(48) = ( 2.00000*CONSTS(81)*ALGBRC(6))/CONSTS(80)
      ALGBRC(100) =  (ALGBRC(48)+CONSTS(152)+ALGBRC(76))*STATES(43)
      ALGBRC(86) = ( 3.00000*CONSTS(76)*ALGBRC(17))/CONSTS(77)
      ALGBRC(108) =  ALGBRC(86)*STATES(44)+ CONSTS(147)*STATES(38)+ ALGBRC(39)*STATES(42)
      RATES(43) = ALGBRC(108) - ALGBRC(100)
      ALGBRC(56) = ALGBRC(6)
      ALGBRC(106) =  (ALGBRC(56)+ALGBRC(85)+CONSTS(148))*STATES(39)
      ALGBRC(95) =  4.00000*ALGBRC(17)
      ALGBRC(115) =  ALGBRC(95)*STATES(40)+ CONSTS(153)*STATES(44)+ ALGBRC(47)*STATES(38)
      RATES(39) = ALGBRC(115) - ALGBRC(106)
      ALGBRC(107) =  (ALGBRC(95)+CONSTS(149))*STATES(40)
      ALGBRC(116) =  ALGBRC(56)*STATES(39)+ CONSTS(154)*STATES(45)
      RATES(40) = ALGBRC(116) - ALGBRC(107)
      ALGBRC(57) = ( CONSTS(82)*ALGBRC(6))/CONSTS(81)
      ALGBRC(109) =  (ALGBRC(57)+CONSTS(153)+ALGBRC(86))*STATES(44)
      ALGBRC(96) = ( 4.00000*CONSTS(77)*ALGBRC(17))/CONSTS(78)
      ALGBRC(117) =  ALGBRC(96)*STATES(45)+ CONSTS(148)*STATES(39)+ ALGBRC(48)*STATES(43)
      RATES(44) = ALGBRC(117) - ALGBRC(109)
      ALGBRC(110) =  (CONSTS(154)+ALGBRC(96))*STATES(45)
      ALGBRC(118) =  CONSTS(149)*STATES(40)+ ALGBRC(57)*STATES(44)
      RATES(45) = ALGBRC(118) - ALGBRC(110)
      ALGBRC(5) =  CONSTS(59)*EXP( CONSTS(60)*STATES(7))
      ALGBRC(46) =  4.00000*ALGBRC(5)
      ALGBRC(37) =  CONSTS(65)*EXP( CONSTS(66)*STATES(7))
      ALGBRC(142) = ALGBRC(37)
      ALGBRC(156) =  (ALGBRC(46)+ALGBRC(142))*STATES(26)
      ALGBRC(16) =  CONSTS(61)*EXP( - CONSTS(62)*STATES(7))
      ALGBRC(122) = ALGBRC(16)
      ALGBRC(26) =  CONSTS(63)*EXP( - CONSTS(64)*STATES(7))
      ALGBRC(149) = ALGBRC(26)
      ALGBRC(164) =  ALGBRC(122)*STATES(27)+ ALGBRC(149)*STATES(31)
      RATES(26) = ALGBRC(164) - ALGBRC(156)
      ALGBRC(84) =  4.00000*CONSTS(71)*ALGBRC(5)
      ALGBRC(157) =  (ALGBRC(149)+ALGBRC(84))*STATES(31)
      ALGBRC(138) = ALGBRC(16)/CONSTS(67)
      ALGBRC(166) =  ALGBRC(142)*STATES(26)+ ALGBRC(138)*STATES(32)
      RATES(31) = ALGBRC(166) - ALGBRC(157)
      ALGBRC(55) =  3.00000*ALGBRC(5)
      ALGBRC(147) =  CONSTS(67)*ALGBRC(37)
      ALGBRC(165) =  (ALGBRC(55)+ALGBRC(122)+ALGBRC(147))*STATES(27)
      ALGBRC(126) =  2.00000*ALGBRC(16)
      ALGBRC(155) = ALGBRC(26)/CONSTS(71)
      ALGBRC(173) =  ALGBRC(126)*STATES(28)+ ALGBRC(155)*STATES(32)+ ALGBRC(46)*STATES(26)
      RATES(27) = ALGBRC(173) - ALGBRC(165)
      ALGBRC(94) = ( 3.00000*CONSTS(72)*ALGBRC(5))/CONSTS(71)
      ALGBRC(167) =  (ALGBRC(94)+ALGBRC(155)+ALGBRC(138))*STATES(32)
      ALGBRC(143) = ( 2.00000*CONSTS(67)*ALGBRC(16))/CONSTS(68)
      ALGBRC(175) =  ALGBRC(143)*STATES(33)+ ALGBRC(147)*STATES(27)+ ALGBRC(84)*STATES(31)
      RATES(32) = ALGBRC(175) - ALGBRC(167)
      ALGBRC(8) =  CONSTS(162)*CONSTS(138)*EXP(- 114007./CONSTS(158)+224.114/CONSTS(6)+( 0.286374*STATES(7))/CONSTS(139))
      ALGBRC(111) =  4.00000*ALGBRC(8)
      ALGBRC(19) =  CONSTS(162)*CONSTS(138)*EXP(- 272470./CONSTS(158)+708.146/CONSTS(6)+( - 2.28528*STATES(7))/CONSTS(139))
      ALGBRC(177) = ALGBRC(19)
      RATES(51) =  - (CONSTS(178)+ALGBRC(111))*STATES(51)+ ALGBRC(177)*STATES(52)+ CONSTS(177)*STATES(56)
      ALGBRC(65) =  2.00000*ALGBRC(5)
      ALGBRC(153) =  CONSTS(68)*ALGBRC(37)
      ALGBRC(174) =  (ALGBRC(65)+ALGBRC(126)+ALGBRC(153))*STATES(28)
      ALGBRC(130) =  3.00000*ALGBRC(16)
      ALGBRC(163) = ALGBRC(26)/CONSTS(72)
      ALGBRC(182) =  ALGBRC(130)*STATES(29)+ ALGBRC(163)*STATES(33)+ ALGBRC(55)*STATES(27)
      RATES(28) = ALGBRC(182) - ALGBRC(174)
      ALGBRC(104) = ( 2.00000*CONSTS(73)*ALGBRC(5))/CONSTS(72)
      ALGBRC(176) =  (ALGBRC(104)+ALGBRC(163)+ALGBRC(143))*STATES(33)
      ALGBRC(148) = ( 3.00000*CONSTS(68)*ALGBRC(16))/CONSTS(69)
      ALGBRC(185) =  ALGBRC(148)*STATES(34)+ ALGBRC(153)*STATES(28)+ ALGBRC(94)*STATES(32)
      RATES(33) = ALGBRC(185) - ALGBRC(176)
      ALGBRC(119) =  3.00000*ALGBRC(8)
      ALGBRC(188) =  2.00000*ALGBRC(19)
      RATES(52) = ( ALGBRC(111)*STATES(51) -  (ALGBRC(177)+ALGBRC(119)+CONSTS(180))*STATES(52))+ ALGBRC(188)*STATES(53)+ CONSTS(179)*STATES(57)
      ALGBRC(3) =  4.00000*1.20000*0.416000*EXP( 0.0120000*(STATES(7) - 35.0000))
      ALGBRC(54) =  4.00000*ALGBRC(3)
      ALGBRC(45) =  0.600000*0.0923300*STATES(4)
      ALGBRC(170) = ALGBRC(45)
      ALGBRC(181) =  (ALGBRC(54)+ALGBRC(170))*STATES(14)
      ALGBRC(14) =  4.00000*0.450000*0.0490000*EXP( - 0.0650000*(STATES(7) - 22.0000))
      ALGBRC(125) = ALGBRC(14)
      ALGBRC(193) =  ALGBRC(125)*STATES(15)+ CONSTS(144)*STATES(19)
      RATES(14) = ALGBRC(193) - ALGBRC(181)
      ALGBRC(25) =  CONSTS(58)*ALGBRC(3)
      ALGBRC(93) =  4.00000*ALGBRC(25)
      ALGBRC(180) =  (ALGBRC(93)+CONSTS(144))*STATES(19)
      ALGBRC(36) = ALGBRC(14)/CONSTS(57)
      ALGBRC(141) = ALGBRC(36)
      ALGBRC(191) =  ALGBRC(141)*STATES(20)+ ALGBRC(170)*STATES(14)
      RATES(19) = ALGBRC(191) - ALGBRC(180)
      ALGBRC(74) = ALGBRC(5)
      ALGBRC(161) =  CONSTS(69)*ALGBRC(37)
      ALGBRC(183) =  (ALGBRC(74)+ALGBRC(130)+ALGBRC(161))*STATES(29)
      ALGBRC(134) =  4.00000*ALGBRC(16)
      ALGBRC(171) = ALGBRC(26)/CONSTS(73)
      ALGBRC(195) =  ALGBRC(134)*STATES(30)+ ALGBRC(171)*STATES(34)+ ALGBRC(65)*STATES(28)
      RATES(29) = ALGBRC(195) - ALGBRC(183)
      ALGBRC(162) =  CONSTS(70)*ALGBRC(37)
      ALGBRC(184) =  (ALGBRC(134)+ALGBRC(162))*STATES(30)
      ALGBRC(172) = ALGBRC(26)/CONSTS(74)
      ALGBRC(196) =  ALGBRC(74)*STATES(29)+ ALGBRC(172)*STATES(35)
      RATES(30) = ALGBRC(196) - ALGBRC(184)
      ALGBRC(114) = ( CONSTS(74)*ALGBRC(5))/CONSTS(73)
      ALGBRC(186) =  (ALGBRC(114)+ALGBRC(171)+ALGBRC(148))*STATES(34)
      ALGBRC(154) = ( 4.00000*CONSTS(69)*ALGBRC(16))/CONSTS(70)
      ALGBRC(197) =  ALGBRC(154)*STATES(35)+ ALGBRC(161)*STATES(29)+ ALGBRC(104)*STATES(33)
      RATES(34) = ALGBRC(197) - ALGBRC(186)
      ALGBRC(187) =  (ALGBRC(172)+ALGBRC(154))*STATES(35)
      ALGBRC(198) =  ALGBRC(162)*STATES(30)+ ALGBRC(114)*STATES(34)
      RATES(35) = ALGBRC(198) - ALGBRC(187)
      ALGBRC(123) =  2.00000*ALGBRC(8)
      ALGBRC(199) =  3.00000*ALGBRC(19)
      RATES(53) = ( ALGBRC(119)*STATES(52) -  (ALGBRC(188)+ALGBRC(123)+CONSTS(182))*STATES(53))+ ALGBRC(199)*STATES(54)+ CONSTS(181)*STATES(58)
      ALGBRC(64) =  3.00000*ALGBRC(3)
      ALGBRC(179) =  CONSTS(58)*ALGBRC(170)
      ALGBRC(194) =  (ALGBRC(125)+ALGBRC(64)+ALGBRC(179))*STATES(15)
      ALGBRC(129) =  2.00000*ALGBRC(14)
      ALGBRC(204) =  ALGBRC(54)*STATES(14)+ ALGBRC(129)*STATES(16)+ CONSTS(160)*STATES(20)
      RATES(15) = ALGBRC(204) - ALGBRC(194)
      ALGBRC(103) =  3.00000*ALGBRC(25)
      ALGBRC(192) =  (ALGBRC(141)+ALGBRC(103)+CONSTS(160))*STATES(20)
      ALGBRC(146) =  2.00000*ALGBRC(36)
      ALGBRC(202) =  ALGBRC(93)*STATES(19)+ ALGBRC(146)*STATES(21)+ ALGBRC(179)*STATES(15)
      RATES(20) = ALGBRC(202) - ALGBRC(192)
      ALGBRC(127) = ALGBRC(8)
      ALGBRC(206) =  4.00000*ALGBRC(19)
      RATES(54) = ( ALGBRC(123)*STATES(53) -  (ALGBRC(199)+ALGBRC(127)+CONSTS(184))*STATES(54))+ ALGBRC(206)*STATES(55)+ CONSTS(183)*STATES(59)
      ALGBRC(73) =  2.00000*ALGBRC(3)
      ALGBRC(190) =  CONSTS(58)*ALGBRC(179)
      ALGBRC(205) =  (ALGBRC(129)+ALGBRC(73)+ALGBRC(190))*STATES(16)
      ALGBRC(133) =  3.00000*ALGBRC(14)
      ALGBRC(211) =  ALGBRC(64)*STATES(15)+ ALGBRC(133)*STATES(17)+ CONSTS(163)*STATES(21)
      RATES(16) = ALGBRC(211) - ALGBRC(205)
      ALGBRC(113) =  2.00000*ALGBRC(25)
      ALGBRC(203) =  (ALGBRC(146)+ALGBRC(113)+CONSTS(163))*STATES(21)
      ALGBRC(152) =  3.00000*ALGBRC(36)
      ALGBRC(209) =  ALGBRC(103)*STATES(20)+ ALGBRC(152)*STATES(22)+ ALGBRC(190)*STATES(16)
      RATES(21) = ALGBRC(209) - ALGBRC(203)
      ALGBRC(83) = ALGBRC(3)
      ALGBRC(201) =  CONSTS(58)*ALGBRC(190)
      ALGBRC(212) =  (ALGBRC(133)+ALGBRC(83)+ALGBRC(201))*STATES(17)
      ALGBRC(137) =  4.00000*ALGBRC(14)
      ALGBRC(217) =  ALGBRC(73)*STATES(16)+ ALGBRC(137)*STATES(18)+ CONSTS(169)*STATES(22)
      RATES(17) = ALGBRC(217) - ALGBRC(212)
      ALGBRC(121) = ALGBRC(25)
      ALGBRC(210) =  (ALGBRC(152)+ALGBRC(121)+CONSTS(169))*STATES(22)
      ALGBRC(160) =  4.00000*ALGBRC(36)
      ALGBRC(215) =  ALGBRC(113)*STATES(21)+ ALGBRC(160)*STATES(23)+ ALGBRC(201)*STATES(17)
      RATES(22) = ALGBRC(215) - ALGBRC(210)
      ALGBRC(29) =  CONSTS(162)*CONSTS(138)*EXP(- 196337./CONSTS(158)+529.952/CONSTS(6)+( 2.78085*STATES(7))/CONSTS(139))
      ALGBRC(131) = ALGBRC(29)
      ALGBRC(91) =  CONSTS(162)*CONSTS(138)*EXP(- 147814./CONSTS(158)+338.915/CONSTS(6)+( 2.13600*STATES(7))/CONSTS(139))
      ALGBRC(168) = ALGBRC(91)
      ALGBRC(40) =  CONSTS(162)*CONSTS(138)*EXP(- 133690./CONSTS(158)+229.205/CONSTS(6)+( - 1.55804*STATES(7))/CONSTS(139))
      ALGBRC(213) = ALGBRC(40)
      ALGBRC(101) =  CONSTS(162)*CONSTS(138)*EXP(- 121322./CONSTS(158)+193.265/CONSTS(6)+( - 1.74290*STATES(7))/CONSTS(139))
      ALGBRC(220) = ALGBRC(101)
      RATES(55) = ( ALGBRC(127)*STATES(54) -  (ALGBRC(206)+ALGBRC(131)+ALGBRC(168)+CONSTS(186))*STATES(55))+ ALGBRC(213)*STATES(46)+ ALGBRC(220)*STATES(47)+ CONSTS(185)*STATES(60)
      RATES(47) = ( ALGBRC(168)*STATES(55)+ CONSTS(171)*STATES(46)) -  (ALGBRC(220)+CONSTS(172))*STATES(47)
      ALGBRC(135) =  ALGBRC(111)*CONSTS(175)
      ALGBRC(219) = ALGBRC(177)/CONSTS(175)
      RATES(56) = ( CONSTS(178)*STATES(51) -  (CONSTS(177)+ALGBRC(135))*STATES(56))+ ALGBRC(219)*STATES(57)
      ALGBRC(208) =  CONSTS(58)*ALGBRC(201)
      ALGBRC(218) =  (ALGBRC(137)+CONSTS(55)+ALGBRC(208))*STATES(18)
      ALGBRC(223) =  ALGBRC(83)*STATES(17)+ CONSTS(56)*STATES(24)+ CONSTS(174)*STATES(23)
      RATES(18) = ALGBRC(223) - ALGBRC(218)
      ALGBRC(216) =  (ALGBRC(160)+CONSTS(174))*STATES(23)
      ALGBRC(222) =  ALGBRC(121)*STATES(22)+ ALGBRC(208)*STATES(18)
      RATES(23) = ALGBRC(222) - ALGBRC(216)
      ALGBRC(139) =  ALGBRC(119)*CONSTS(175)
      ALGBRC(224) = ALGBRC(188)/CONSTS(175)
      RATES(57) = (( CONSTS(180)*STATES(52)+ ALGBRC(135)*STATES(56)) -  (ALGBRC(219)+CONSTS(179)+ALGBRC(139))*STATES(57))+ ALGBRC(224)*STATES(58)
      ALGBRC(144) =  ALGBRC(123)*CONSTS(175)
      ALGBRC(226) = ALGBRC(199)/CONSTS(175)
      RATES(58) = (( CONSTS(182)*STATES(53)+ ALGBRC(139)*STATES(57)) -  (ALGBRC(144)+CONSTS(181)+ALGBRC(224))*STATES(58))+ ALGBRC(226)*STATES(59)
      ALGBRC(150) =  ALGBRC(127)*CONSTS(175)
      ALGBRC(228) = ALGBRC(206)/CONSTS(175)
      RATES(59) = (( CONSTS(184)*STATES(54)+ ALGBRC(144)*STATES(58)) -  (ALGBRC(226)+CONSTS(183)+ALGBRC(150))*STATES(59))+ ALGBRC(228)*STATES(60)
      ALGBRC(68) =  CONSTS(162)*CONSTS(138)*EXP(116431./CONSTS(158)+- 578.317/CONSTS(6)+( 0.764126*STATES(7))/CONSTS(139))
      ALGBRC(158) = ALGBRC(68)
      ALGBRC(79) =  CONSTS(162)*CONSTS(138)*EXP(- 55700.7/CONSTS(158)+- 130.639/CONSTS(6)+( - 3.64982*STATES(7))/CONSTS(139))
      ALGBRC(230) = ALGBRC(79)
      RATES(60) = (( CONSTS(186)*STATES(55)+ ALGBRC(150)*STATES(59)) -  (ALGBRC(228)+CONSTS(185)+ALGBRC(158))*STATES(60))+ ALGBRC(230)*STATES(61)
      ALGBRC(58) =  CONSTS(162)*CONSTS(138)*EXP(- 97657.8/CONSTS(158)+1.51000/CONSTS(6)+( 0.0684862*STATES(7))/CONSTS(139))
      ALGBRC(232) = ALGBRC(58)
      ALGBRC(49) =  CONSTS(162)*CONSTS(138)*EXP(- 62123.1/CONSTS(158)+39.2950/CONSTS(6)+( 0.288816*STATES(7))/CONSTS(139))
      ALGBRC(234) = ALGBRC(49)
      RATES(46) = ( ALGBRC(131)*STATES(55) -  (ALGBRC(213)+CONSTS(171)+ALGBRC(234))*STATES(46))+ CONSTS(172)*STATES(47)+ ALGBRC(232)*STATES(61)
      RATES(61) = ( ALGBRC(234)*STATES(46)+ ALGBRC(158)*STATES(60)) -  (ALGBRC(230)+ALGBRC(232))*STATES(61)
      ALGBRC(120) =  CONSTS(125)*log(CONSTS(13)/STATES(1))
      ALGBRC(132) =  CONSTS(92)*(STATES(46)+STATES(47))*(STATES(7) - ALGBRC(120))
      ALGBRC(229) =  CONSTS(93)*(STATES(7) - ALGBRC(120))
      ALGBRC(235) = STATES(7)/CONSTS(125)
      ALGBRC(241) =  EXP( CONSTS(100)*ALGBRC(235))*STATES(1) ** 3.00000*CONSTS(14)
      ALGBRC(242) =  EXP( (CONSTS(100) - 1.00000)*ALGBRC(235))*CONSTS(13) ** 3.00000*STATES(3)
      ALGBRC(243) = 1.00000+ CONSTS(99)*EXP( (CONSTS(100) - 1.00000)*ALGBRC(235))
      ALGBRC(244) = ( CONSTS(96)*(ALGBRC(241) - ALGBRC(242)))/( CONSTS(173)*ALGBRC(243)*CONSTS(176))
      ALGBRC(236) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(235))
      ALGBRC(237) =  0.0365000*CONSTS(164)*EXP( - 1.33000*ALGBRC(235))
      ALGBRC(238) = 1.00000/(ALGBRC(236)+ALGBRC(237))
      ALGBRC(239) = 1.00000+CONSTS(102)/STATES(1) ** 1.50000
      ALGBRC(240) = ( CONSTS(101)*ALGBRC(238)*CONSTS(168))/ALGBRC(239)
      ALGBRC(151) = STATES(7)/CONSTS(125)
      ALGBRC(159) =  1000.00*CONSTS(4)*ALGBRC(151)
      ALGBRC(178) = EXP(ALGBRC(151)) - 1.00000
      ALGBRC(200) =  STATES(1)*EXP(ALGBRC(151)) - CONSTS(13)
      ALGBRC(207) = ( 0.0200000*CONSTS(134)*STATES(40)*ALGBRC(159)*ALGBRC(200))/ALGBRC(178)
      RATES(1) =  - (ALGBRC(132)+ALGBRC(229)+ 3.00000*(ALGBRC(244)+ALGBRC(240))+ALGBRC(207))*CONSTS(156)
      ALGBRC(128) =  0.500000*CONSTS(125)*log(CONSTS(14)/STATES(3))
      ALGBRC(231) =  CONSTS(106)*(STATES(7) - ALGBRC(128))
      ALGBRC(233) = ( CONSTS(104)*STATES(3))/(CONSTS(105)+STATES(3))
      ALGBRC(61) = (STATES(4) - STATES(3))/CONSTS(28)
      ALGBRC(92) =  CONSTS(29)*ALGBRC(71)+ CONSTS(30)*ALGBRC(82)
      ALGBRC(112) = 1.00000/(1.00000+( CONSTS(35)*CONSTS(37))/STATES(3)+CONSTS(37) ** 2.00000+( CONSTS(36)*CONSTS(39))/STATES(3)+CONSTS(39) ** 2.00000)
      RATES(3) =  ALGBRC(112)*(((ALGBRC(61) - ALGBRC(43)) - ALGBRC(92)) -  ((ALGBRC(231) -  2.00000*ALGBRC(244))+ALGBRC(233))*0.500000*CONSTS(156))
      ALGBRC(245) = STATES(7)/CONSTS(125)
      ALGBRC(246) =  1000.00*CONSTS(4)*ALGBRC(245)
      ALGBRC(247) =  0.00100000*EXP( 2.00000*ALGBRC(245)) -  CONSTS(14)*0.341000
      ALGBRC(248) = EXP( 2.00000*ALGBRC(245)) - 1.00000
      ALGBRC(249) = ( CONSTS(157)*4.00000*ALGBRC(246)*ALGBRC(247))/ALGBRC(248)
      ALGBRC(250) =  ALGBRC(249)*STATES(25)*STATES(24)
      ALGBRC(102) = 1.00000/(1.00000+( CONSTS(35)*CONSTS(37))/STATES(4)+CONSTS(37) ** 2.00000+( CONSTS(36)*CONSTS(39))/STATES(4)+CONSTS(39) ** 2.00000)
      RATES(4) =  ALGBRC(102)*((( ALGBRC(52)*CONSTS(9))/CONSTS(11) - ( ALGBRC(61)*CONSTS(8))/CONSTS(11)) -  ALGBRC(250)*CONSTS(161))
      ALGBRC(124) =  CONSTS(125)*log(CONSTS(12)/STATES(2))
      ALGBRC(136) =  CONSTS(89)*CONSTS(135)*STATES(48)*(STATES(7) - ALGBRC(124))
      ALGBRC(140) =  CONSTS(90)*(STATES(49)+STATES(50))*(STATES(7) - ALGBRC(124))
      ALGBRC(225) = 1.00000/(0.940000+EXP( (1.26000/CONSTS(125))*(STATES(7) - ALGBRC(124))))
      ALGBRC(227) =  CONSTS(91)* (CONSTS(12)/1.00000) ** (1.0 / 2)*ALGBRC(225)*(STATES(7) - ALGBRC(124))
      ALGBRC(251) = TERNRY(ALGBRC(249).GE.0.00000, 0.00000, ALGBRC(249))
      ALGBRC(252) = CONSTS(137)/(1.00000+ALGBRC(251)/CONSTS(107))
      ALGBRC(253) =  STATES(2)*EXP(ALGBRC(245)) - CONSTS(12)
      ALGBRC(254) = EXP(ALGBRC(245)) - 1.00000
      ALGBRC(255) = ( ALGBRC(252)*STATES(24)*STATES(25)*ALGBRC(246)*ALGBRC(253))/ALGBRC(254)
      ALGBRC(1) =  INT(VOI/CONSTS(15))*CONSTS(15)
      ALGBRC(12) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(18).AND.VOI - ALGBRC(1).LE.CONSTS(18)+CONSTS(16), CONSTS(17), 0.00000)
      ALGBRC(145) =  CONSTS(133)*STATES(30)*(STATES(7) - ALGBRC(124))
      ALGBRC(169) =  STATES(2)*EXP(ALGBRC(151)) - CONSTS(12)
      ALGBRC(189) = ( CONSTS(134)*STATES(40)*ALGBRC(159)*ALGBRC(169))/ALGBRC(178)
      RATES(2) =  - (((ALGBRC(136)+ALGBRC(140)+ALGBRC(227)+ALGBRC(255)+ALGBRC(12)) -  2.00000*ALGBRC(240))+ALGBRC(145)+ALGBRC(189))*CONSTS(156)
      ALGBRC(214) = ALGBRC(189)+ALGBRC(207)
      ALGBRC(221) = ALGBRC(145)+ALGBRC(214)
      ALGBRC(256) = ALGBRC(132)+ALGBRC(250)+ALGBRC(255)+ALGBRC(136)+ALGBRC(140)+ALGBRC(227)+ALGBRC(244)+ALGBRC(240)+ALGBRC(221)+ALGBRC(233)+ALGBRC(231)+ALGBRC(229)+ALGBRC(12)
      RATES(7) = - ALGBRC(256)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = (CONSTS(130)/(1.00000+EXP((STATES(7)+28.5000)/7.80000))+1.00000) - CONSTS(130)
      ALGBRC(15) = 1.00000/(0.00336336/(0.500000+EXP(STATES(7)/- 5.53900))+ 0.00779047*EXP(STATES(7)/- 49.5104))
      ALGBRC(9) =  CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*STATES(7))
      ALGBRC(20) =  CONSTS(110)*CONSTS(113)*EXP( CONSTS(114)*STATES(7))
      ALGBRC(33) =  0.216256*EXP( - 1.88912e-05*STATES(7))
      ALGBRC(22) =  0.00700807*EXP( - 0.149998*STATES(7))
      ALGBRC(11) =  0.00767254*EXP( 0.0866295*STATES(7))
      ALGBRC(32) =  0.00379738*EXP( - 0.0142567*STATES(7))
      ALGBRC(2) = CONSTS(49) - (CONSTS(49) - CONSTS(51))/(1.00000+CONSTS(52)/STATES(5) ** CONSTS(54))
      ALGBRC(35) =  - CONSTS(40)*ALGBRC(2)* STATES(4)*1000.00 **  CONSTS(46)*CONSTS(48)*STATES(12)+ CONSTS(41)*STATES(8)
      ALGBRC(31) =  ALGBRC(9)*STATES(62)+ CONSTS(141)*STATES(64)
      ALGBRC(42) =  (ALGBRC(20)+CONSTS(140))*STATES(63)
      ALGBRC(44) =  CONSTS(42)* STATES(4)*1000.00 **  CONSTS(47)*CONSTS(48)*STATES(8) -  CONSTS(43)*STATES(9)
      ALGBRC(13) = CONSTS(127) - (CONSTS(127) - CONSTS(51))/(1.00000+CONSTS(53)/STATES(5) ** CONSTS(54))
      ALGBRC(24) = CONSTS(50) - (CONSTS(50) - CONSTS(51))/(1.00000+CONSTS(128)/STATES(5) ** CONSTS(54))
      ALGBRC(53) =  (CONSTS(44)/ALGBRC(13))*STATES(8) -  CONSTS(45)*ALGBRC(24)*STATES(13)
      ALGBRC(10) =  CONSTS(110)*CONSTS(115)*EXP( CONSTS(116)*STATES(7))
      ALGBRC(41) =  CONSTS(110)*CONSTS(121)*EXP( CONSTS(122)*STATES(7))
      ALGBRC(51) =  ALGBRC(10)*STATES(64)+ ALGBRC(41)*STATES(65)
      ALGBRC(21) =  CONSTS(110)*CONSTS(117)*EXP( CONSTS(118)*STATES(7))
      ALGBRC(30) =  CONSTS(110)*CONSTS(119)*EXP( CONSTS(120)*STATES(7))
      ALGBRC(60) =  (ALGBRC(21)+ALGBRC(30))*STATES(48)
      ALGBRC(23) = STATES(3)/CONSTS(19) ** CONSTS(20)
      ALGBRC(34) = STATES(6)/CONSTS(21) ** CONSTS(22)
      ALGBRC(43) = TERNRY(CONSTS(1).EQ.0.00000, ( CONSTS(23)*( CONSTS(24)*ALGBRC(23) -  CONSTS(25)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34)), ( 1.50000*CONSTS(23)*( CONSTS(24)*ALGBRC(23) -  CONSTS(25)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34)))
      ALGBRC(62) = (STATES(6) - STATES(5))/CONSTS(27)
      ALGBRC(63) = - (ALGBRC(35)+ALGBRC(44)+ALGBRC(53))
      ALGBRC(52) =  CONSTS(26)*(STATES(8)+STATES(9))*(STATES(5) - STATES(4))
      ALGBRC(72) = 1.00000/(1.00000+( CONSTS(126)*CONSTS(38))/STATES(5)+CONSTS(38) ** 2.00000)
      ALGBRC(71) =  CONSTS(33)*STATES(3)*(1.00000 - STATES(10)) -  CONSTS(34)*STATES(10)
      ALGBRC(50) =  CONSTS(110)*CONSTS(123)*EXP( CONSTS(124)*STATES(7))
      ALGBRC(59) = ( ALGBRC(21)*ALGBRC(41)*ALGBRC(50))/( ALGBRC(10)*ALGBRC(30))
      ALGBRC(69) =  CONSTS(140)*STATES(63)+ ALGBRC(21)*STATES(48)+ ALGBRC(59)*STATES(65)
      ALGBRC(80) =  (ALGBRC(50)+ALGBRC(10)+CONSTS(141))*STATES(64)
      ALGBRC(70) =  ALGBRC(50)*STATES(64)+ ALGBRC(30)*STATES(48)
      ALGBRC(81) =  (ALGBRC(59)+ALGBRC(41))*STATES(65)
      ALGBRC(82) =  CONSTS(31)*STATES(3)*(1.00000 - STATES(11)) -  CONSTS(32)*STATES(11)
      ALGBRC(6) =  CONSTS(83)*EXP( CONSTS(84)*STATES(7))
      ALGBRC(27) =  4.00000*ALGBRC(6)
      ALGBRC(77) =  (ALGBRC(27)+CONSTS(145))*STATES(36)
      ALGBRC(17) =  CONSTS(85)*EXP( - CONSTS(86)*STATES(7))
      ALGBRC(66) = ALGBRC(17)
      ALGBRC(87) =  ALGBRC(66)*STATES(37)+ CONSTS(150)*STATES(41)
      ALGBRC(28) =  4.00000*CONSTS(79)*ALGBRC(6)
      ALGBRC(78) =  (CONSTS(150)+ALGBRC(28))*STATES(41)
      ALGBRC(67) = ALGBRC(17)/CONSTS(75)
      ALGBRC(89) =  CONSTS(145)*STATES(36)+ ALGBRC(67)*STATES(42)
      ALGBRC(38) =  3.00000*ALGBRC(6)
      ALGBRC(88) =  (ALGBRC(38)+ALGBRC(66)+CONSTS(146))*STATES(37)
      ALGBRC(75) =  2.00000*ALGBRC(17)
      ALGBRC(97) =  ALGBRC(75)*STATES(38)+ CONSTS(151)*STATES(42)+ ALGBRC(27)*STATES(36)
      ALGBRC(39) = ( 3.00000*CONSTS(80)*ALGBRC(6))/CONSTS(79)
      ALGBRC(90) =  (ALGBRC(39)+CONSTS(151)+ALGBRC(67))*STATES(42)
      ALGBRC(76) = ( 2.00000*CONSTS(75)*ALGBRC(17))/CONSTS(76)
      ALGBRC(99) =  ALGBRC(76)*STATES(43)+ CONSTS(146)*STATES(37)+ ALGBRC(28)*STATES(41)
      ALGBRC(47) =  2.00000*ALGBRC(6)
      ALGBRC(98) =  (ALGBRC(47)+ALGBRC(75)+CONSTS(147))*STATES(38)
      ALGBRC(85) =  3.00000*ALGBRC(17)
      ALGBRC(105) =  ALGBRC(85)*STATES(39)+ CONSTS(152)*STATES(43)+ ALGBRC(38)*STATES(37)
      ALGBRC(48) = ( 2.00000*CONSTS(81)*ALGBRC(6))/CONSTS(80)
      ALGBRC(100) =  (ALGBRC(48)+CONSTS(152)+ALGBRC(76))*STATES(43)
      ALGBRC(86) = ( 3.00000*CONSTS(76)*ALGBRC(17))/CONSTS(77)
      ALGBRC(108) =  ALGBRC(86)*STATES(44)+ CONSTS(147)*STATES(38)+ ALGBRC(39)*STATES(42)
      ALGBRC(56) = ALGBRC(6)
      ALGBRC(106) =  (ALGBRC(56)+ALGBRC(85)+CONSTS(148))*STATES(39)
      ALGBRC(95) =  4.00000*ALGBRC(17)
      ALGBRC(115) =  ALGBRC(95)*STATES(40)+ CONSTS(153)*STATES(44)+ ALGBRC(47)*STATES(38)
      ALGBRC(107) =  (ALGBRC(95)+CONSTS(149))*STATES(40)
      ALGBRC(116) =  ALGBRC(56)*STATES(39)+ CONSTS(154)*STATES(45)
      ALGBRC(57) = ( CONSTS(82)*ALGBRC(6))/CONSTS(81)
      ALGBRC(109) =  (ALGBRC(57)+CONSTS(153)+ALGBRC(86))*STATES(44)
      ALGBRC(96) = ( 4.00000*CONSTS(77)*ALGBRC(17))/CONSTS(78)
      ALGBRC(117) =  ALGBRC(96)*STATES(45)+ CONSTS(148)*STATES(39)+ ALGBRC(48)*STATES(43)
      ALGBRC(110) =  (CONSTS(154)+ALGBRC(96))*STATES(45)
      ALGBRC(118) =  CONSTS(149)*STATES(40)+ ALGBRC(57)*STATES(44)
      ALGBRC(5) =  CONSTS(59)*EXP( CONSTS(60)*STATES(7))
      ALGBRC(46) =  4.00000*ALGBRC(5)
      ALGBRC(37) =  CONSTS(65)*EXP( CONSTS(66)*STATES(7))
      ALGBRC(142) = ALGBRC(37)
      ALGBRC(156) =  (ALGBRC(46)+ALGBRC(142))*STATES(26)
      ALGBRC(16) =  CONSTS(61)*EXP( - CONSTS(62)*STATES(7))
      ALGBRC(122) = ALGBRC(16)
      ALGBRC(26) =  CONSTS(63)*EXP( - CONSTS(64)*STATES(7))
      ALGBRC(149) = ALGBRC(26)
      ALGBRC(164) =  ALGBRC(122)*STATES(27)+ ALGBRC(149)*STATES(31)
      ALGBRC(84) =  4.00000*CONSTS(71)*ALGBRC(5)
      ALGBRC(157) =  (ALGBRC(149)+ALGBRC(84))*STATES(31)
      ALGBRC(138) = ALGBRC(16)/CONSTS(67)
      ALGBRC(166) =  ALGBRC(142)*STATES(26)+ ALGBRC(138)*STATES(32)
      ALGBRC(55) =  3.00000*ALGBRC(5)
      ALGBRC(147) =  CONSTS(67)*ALGBRC(37)
      ALGBRC(165) =  (ALGBRC(55)+ALGBRC(122)+ALGBRC(147))*STATES(27)
      ALGBRC(126) =  2.00000*ALGBRC(16)
      ALGBRC(155) = ALGBRC(26)/CONSTS(71)
      ALGBRC(173) =  ALGBRC(126)*STATES(28)+ ALGBRC(155)*STATES(32)+ ALGBRC(46)*STATES(26)
      ALGBRC(94) = ( 3.00000*CONSTS(72)*ALGBRC(5))/CONSTS(71)
      ALGBRC(167) =  (ALGBRC(94)+ALGBRC(155)+ALGBRC(138))*STATES(32)
      ALGBRC(143) = ( 2.00000*CONSTS(67)*ALGBRC(16))/CONSTS(68)
      ALGBRC(175) =  ALGBRC(143)*STATES(33)+ ALGBRC(147)*STATES(27)+ ALGBRC(84)*STATES(31)
      ALGBRC(8) =  CONSTS(162)*CONSTS(138)*EXP(- 114007./CONSTS(158)+224.114/CONSTS(6)+( 0.286374*STATES(7))/CONSTS(139))
      ALGBRC(111) =  4.00000*ALGBRC(8)
      ALGBRC(19) =  CONSTS(162)*CONSTS(138)*EXP(- 272470./CONSTS(158)+708.146/CONSTS(6)+( - 2.28528*STATES(7))/CONSTS(139))
      ALGBRC(177) = ALGBRC(19)
      ALGBRC(65) =  2.00000*ALGBRC(5)
      ALGBRC(153) =  CONSTS(68)*ALGBRC(37)
      ALGBRC(174) =  (ALGBRC(65)+ALGBRC(126)+ALGBRC(153))*STATES(28)
      ALGBRC(130) =  3.00000*ALGBRC(16)
      ALGBRC(163) = ALGBRC(26)/CONSTS(72)
      ALGBRC(182) =  ALGBRC(130)*STATES(29)+ ALGBRC(163)*STATES(33)+ ALGBRC(55)*STATES(27)
      ALGBRC(104) = ( 2.00000*CONSTS(73)*ALGBRC(5))/CONSTS(72)
      ALGBRC(176) =  (ALGBRC(104)+ALGBRC(163)+ALGBRC(143))*STATES(33)
      ALGBRC(148) = ( 3.00000*CONSTS(68)*ALGBRC(16))/CONSTS(69)
      ALGBRC(185) =  ALGBRC(148)*STATES(34)+ ALGBRC(153)*STATES(28)+ ALGBRC(94)*STATES(32)
      ALGBRC(119) =  3.00000*ALGBRC(8)
      ALGBRC(188) =  2.00000*ALGBRC(19)
      ALGBRC(3) =  4.00000*1.20000*0.416000*EXP( 0.0120000*(STATES(7) - 35.0000))
      ALGBRC(54) =  4.00000*ALGBRC(3)
      ALGBRC(45) =  0.600000*0.0923300*STATES(4)
      ALGBRC(170) = ALGBRC(45)
      ALGBRC(181) =  (ALGBRC(54)+ALGBRC(170))*STATES(14)
      ALGBRC(14) =  4.00000*0.450000*0.0490000*EXP( - 0.0650000*(STATES(7) - 22.0000))
      ALGBRC(125) = ALGBRC(14)
      ALGBRC(193) =  ALGBRC(125)*STATES(15)+ CONSTS(144)*STATES(19)
      ALGBRC(25) =  CONSTS(58)*ALGBRC(3)
      ALGBRC(93) =  4.00000*ALGBRC(25)
      ALGBRC(180) =  (ALGBRC(93)+CONSTS(144))*STATES(19)
      ALGBRC(36) = ALGBRC(14)/CONSTS(57)
      ALGBRC(141) = ALGBRC(36)
      ALGBRC(191) =  ALGBRC(141)*STATES(20)+ ALGBRC(170)*STATES(14)
      ALGBRC(74) = ALGBRC(5)
      ALGBRC(161) =  CONSTS(69)*ALGBRC(37)
      ALGBRC(183) =  (ALGBRC(74)+ALGBRC(130)+ALGBRC(161))*STATES(29)
      ALGBRC(134) =  4.00000*ALGBRC(16)
      ALGBRC(171) = ALGBRC(26)/CONSTS(73)
      ALGBRC(195) =  ALGBRC(134)*STATES(30)+ ALGBRC(171)*STATES(34)+ ALGBRC(65)*STATES(28)
      ALGBRC(162) =  CONSTS(70)*ALGBRC(37)
      ALGBRC(184) =  (ALGBRC(134)+ALGBRC(162))*STATES(30)
      ALGBRC(172) = ALGBRC(26)/CONSTS(74)
      ALGBRC(196) =  ALGBRC(74)*STATES(29)+ ALGBRC(172)*STATES(35)
      ALGBRC(114) = ( CONSTS(74)*ALGBRC(5))/CONSTS(73)
      ALGBRC(186) =  (ALGBRC(114)+ALGBRC(171)+ALGBRC(148))*STATES(34)
      ALGBRC(154) = ( 4.00000*CONSTS(69)*ALGBRC(16))/CONSTS(70)
      ALGBRC(197) =  ALGBRC(154)*STATES(35)+ ALGBRC(161)*STATES(29)+ ALGBRC(104)*STATES(33)
      ALGBRC(187) =  (ALGBRC(172)+ALGBRC(154))*STATES(35)
      ALGBRC(198) =  ALGBRC(162)*STATES(30)+ ALGBRC(114)*STATES(34)
      ALGBRC(123) =  2.00000*ALGBRC(8)
      ALGBRC(199) =  3.00000*ALGBRC(19)
      ALGBRC(64) =  3.00000*ALGBRC(3)
      ALGBRC(179) =  CONSTS(58)*ALGBRC(170)
      ALGBRC(194) =  (ALGBRC(125)+ALGBRC(64)+ALGBRC(179))*STATES(15)
      ALGBRC(129) =  2.00000*ALGBRC(14)
      ALGBRC(204) =  ALGBRC(54)*STATES(14)+ ALGBRC(129)*STATES(16)+ CONSTS(160)*STATES(20)
      ALGBRC(103) =  3.00000*ALGBRC(25)
      ALGBRC(192) =  (ALGBRC(141)+ALGBRC(103)+CONSTS(160))*STATES(20)
      ALGBRC(146) =  2.00000*ALGBRC(36)
      ALGBRC(202) =  ALGBRC(93)*STATES(19)+ ALGBRC(146)*STATES(21)+ ALGBRC(179)*STATES(15)
      ALGBRC(127) = ALGBRC(8)
      ALGBRC(206) =  4.00000*ALGBRC(19)
      ALGBRC(73) =  2.00000*ALGBRC(3)
      ALGBRC(190) =  CONSTS(58)*ALGBRC(179)
      ALGBRC(205) =  (ALGBRC(129)+ALGBRC(73)+ALGBRC(190))*STATES(16)
      ALGBRC(133) =  3.00000*ALGBRC(14)
      ALGBRC(211) =  ALGBRC(64)*STATES(15)+ ALGBRC(133)*STATES(17)+ CONSTS(163)*STATES(21)
      ALGBRC(113) =  2.00000*ALGBRC(25)
      ALGBRC(203) =  (ALGBRC(146)+ALGBRC(113)+CONSTS(163))*STATES(21)
      ALGBRC(152) =  3.00000*ALGBRC(36)
      ALGBRC(209) =  ALGBRC(103)*STATES(20)+ ALGBRC(152)*STATES(22)+ ALGBRC(190)*STATES(16)
      ALGBRC(83) = ALGBRC(3)
      ALGBRC(201) =  CONSTS(58)*ALGBRC(190)
      ALGBRC(212) =  (ALGBRC(133)+ALGBRC(83)+ALGBRC(201))*STATES(17)
      ALGBRC(137) =  4.00000*ALGBRC(14)
      ALGBRC(217) =  ALGBRC(73)*STATES(16)+ ALGBRC(137)*STATES(18)+ CONSTS(169)*STATES(22)
      ALGBRC(121) = ALGBRC(25)
      ALGBRC(210) =  (ALGBRC(152)+ALGBRC(121)+CONSTS(169))*STATES(22)
      ALGBRC(160) =  4.00000*ALGBRC(36)
      ALGBRC(215) =  ALGBRC(113)*STATES(21)+ ALGBRC(160)*STATES(23)+ ALGBRC(201)*STATES(17)
      ALGBRC(29) =  CONSTS(162)*CONSTS(138)*EXP(- 196337./CONSTS(158)+529.952/CONSTS(6)+( 2.78085*STATES(7))/CONSTS(139))
      ALGBRC(131) = ALGBRC(29)
      ALGBRC(91) =  CONSTS(162)*CONSTS(138)*EXP(- 147814./CONSTS(158)+338.915/CONSTS(6)+( 2.13600*STATES(7))/CONSTS(139))
      ALGBRC(168) = ALGBRC(91)
      ALGBRC(40) =  CONSTS(162)*CONSTS(138)*EXP(- 133690./CONSTS(158)+229.205/CONSTS(6)+( - 1.55804*STATES(7))/CONSTS(139))
      ALGBRC(213) = ALGBRC(40)
      ALGBRC(101) =  CONSTS(162)*CONSTS(138)*EXP(- 121322./CONSTS(158)+193.265/CONSTS(6)+( - 1.74290*STATES(7))/CONSTS(139))
      ALGBRC(220) = ALGBRC(101)
      ALGBRC(135) =  ALGBRC(111)*CONSTS(175)
      ALGBRC(219) = ALGBRC(177)/CONSTS(175)
      ALGBRC(208) =  CONSTS(58)*ALGBRC(201)
      ALGBRC(218) =  (ALGBRC(137)+CONSTS(55)+ALGBRC(208))*STATES(18)
      ALGBRC(223) =  ALGBRC(83)*STATES(17)+ CONSTS(56)*STATES(24)+ CONSTS(174)*STATES(23)
      ALGBRC(216) =  (ALGBRC(160)+CONSTS(174))*STATES(23)
      ALGBRC(222) =  ALGBRC(121)*STATES(22)+ ALGBRC(208)*STATES(18)
      ALGBRC(139) =  ALGBRC(119)*CONSTS(175)
      ALGBRC(224) = ALGBRC(188)/CONSTS(175)
      ALGBRC(144) =  ALGBRC(123)*CONSTS(175)
      ALGBRC(226) = ALGBRC(199)/CONSTS(175)
      ALGBRC(150) =  ALGBRC(127)*CONSTS(175)
      ALGBRC(228) = ALGBRC(206)/CONSTS(175)
      ALGBRC(68) =  CONSTS(162)*CONSTS(138)*EXP(116431./CONSTS(158)+- 578.317/CONSTS(6)+( 0.764126*STATES(7))/CONSTS(139))
      ALGBRC(158) = ALGBRC(68)
      ALGBRC(79) =  CONSTS(162)*CONSTS(138)*EXP(- 55700.7/CONSTS(158)+- 130.639/CONSTS(6)+( - 3.64982*STATES(7))/CONSTS(139))
      ALGBRC(230) = ALGBRC(79)
      ALGBRC(58) =  CONSTS(162)*CONSTS(138)*EXP(- 97657.8/CONSTS(158)+1.51000/CONSTS(6)+( 0.0684862*STATES(7))/CONSTS(139))
      ALGBRC(232) = ALGBRC(58)
      ALGBRC(49) =  CONSTS(162)*CONSTS(138)*EXP(- 62123.1/CONSTS(158)+39.2950/CONSTS(6)+( 0.288816*STATES(7))/CONSTS(139))
      ALGBRC(234) = ALGBRC(49)
      ALGBRC(120) =  CONSTS(125)*log(CONSTS(13)/STATES(1))
      ALGBRC(132) =  CONSTS(92)*(STATES(46)+STATES(47))*(STATES(7) - ALGBRC(120))
      ALGBRC(229) =  CONSTS(93)*(STATES(7) - ALGBRC(120))
      ALGBRC(235) = STATES(7)/CONSTS(125)
      ALGBRC(241) =  EXP( CONSTS(100)*ALGBRC(235))*STATES(1) ** 3.00000*CONSTS(14)
      ALGBRC(242) =  EXP( (CONSTS(100) - 1.00000)*ALGBRC(235))*CONSTS(13) ** 3.00000*STATES(3)
      ALGBRC(243) = 1.00000+ CONSTS(99)*EXP( (CONSTS(100) - 1.00000)*ALGBRC(235))
      ALGBRC(244) = ( CONSTS(96)*(ALGBRC(241) - ALGBRC(242)))/( CONSTS(173)*ALGBRC(243)*CONSTS(176))
      ALGBRC(236) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(235))
      ALGBRC(237) =  0.0365000*CONSTS(164)*EXP( - 1.33000*ALGBRC(235))
      ALGBRC(238) = 1.00000/(ALGBRC(236)+ALGBRC(237))
      ALGBRC(239) = 1.00000+CONSTS(102)/STATES(1) ** 1.50000
      ALGBRC(240) = ( CONSTS(101)*ALGBRC(238)*CONSTS(168))/ALGBRC(239)
      ALGBRC(151) = STATES(7)/CONSTS(125)
      ALGBRC(159) =  1000.00*CONSTS(4)*ALGBRC(151)
      ALGBRC(178) = EXP(ALGBRC(151)) - 1.00000
      ALGBRC(200) =  STATES(1)*EXP(ALGBRC(151)) - CONSTS(13)
      ALGBRC(207) = ( 0.0200000*CONSTS(134)*STATES(40)*ALGBRC(159)*ALGBRC(200))/ALGBRC(178)
      ALGBRC(128) =  0.500000*CONSTS(125)*log(CONSTS(14)/STATES(3))
      ALGBRC(231) =  CONSTS(106)*(STATES(7) - ALGBRC(128))
      ALGBRC(233) = ( CONSTS(104)*STATES(3))/(CONSTS(105)+STATES(3))
      ALGBRC(61) = (STATES(4) - STATES(3))/CONSTS(28)
      ALGBRC(92) =  CONSTS(29)*ALGBRC(71)+ CONSTS(30)*ALGBRC(82)
      ALGBRC(112) = 1.00000/(1.00000+( CONSTS(35)*CONSTS(37))/STATES(3)+CONSTS(37) ** 2.00000+( CONSTS(36)*CONSTS(39))/STATES(3)+CONSTS(39) ** 2.00000)
      ALGBRC(245) = STATES(7)/CONSTS(125)
      ALGBRC(246) =  1000.00*CONSTS(4)*ALGBRC(245)
      ALGBRC(247) =  0.00100000*EXP( 2.00000*ALGBRC(245)) -  CONSTS(14)*0.341000
      ALGBRC(248) = EXP( 2.00000*ALGBRC(245)) - 1.00000
      ALGBRC(249) = ( CONSTS(157)*4.00000*ALGBRC(246)*ALGBRC(247))/ALGBRC(248)
      ALGBRC(250) =  ALGBRC(249)*STATES(25)*STATES(24)
      ALGBRC(102) = 1.00000/(1.00000+( CONSTS(35)*CONSTS(37))/STATES(4)+CONSTS(37) ** 2.00000+( CONSTS(36)*CONSTS(39))/STATES(4)+CONSTS(39) ** 2.00000)
      ALGBRC(124) =  CONSTS(125)*log(CONSTS(12)/STATES(2))
      ALGBRC(136) =  CONSTS(89)*CONSTS(135)*STATES(48)*(STATES(7) - ALGBRC(124))
      ALGBRC(140) =  CONSTS(90)*(STATES(49)+STATES(50))*(STATES(7) - ALGBRC(124))
      ALGBRC(225) = 1.00000/(0.940000+EXP( (1.26000/CONSTS(125))*(STATES(7) - ALGBRC(124))))
      ALGBRC(227) =  CONSTS(91)* (CONSTS(12)/1.00000) ** (1.0 / 2)*ALGBRC(225)*(STATES(7) - ALGBRC(124))
      ALGBRC(251) = TERNRY(ALGBRC(249).GE.0.00000, 0.00000, ALGBRC(249))
      ALGBRC(252) = CONSTS(137)/(1.00000+ALGBRC(251)/CONSTS(107))
      ALGBRC(253) =  STATES(2)*EXP(ALGBRC(245)) - CONSTS(12)
      ALGBRC(254) = EXP(ALGBRC(245)) - 1.00000
      ALGBRC(255) = ( ALGBRC(252)*STATES(24)*STATES(25)*ALGBRC(246)*ALGBRC(253))/ALGBRC(254)
      ALGBRC(1) =  INT(VOI/CONSTS(15))*CONSTS(15)
      ALGBRC(12) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(18).AND.VOI - ALGBRC(1).LE.CONSTS(18)+CONSTS(16), CONSTS(17), 0.00000)
      ALGBRC(145) =  CONSTS(133)*STATES(30)*(STATES(7) - ALGBRC(124))
      ALGBRC(169) =  STATES(2)*EXP(ALGBRC(151)) - CONSTS(12)
      ALGBRC(189) = ( CONSTS(134)*STATES(40)*ALGBRC(159)*ALGBRC(169))/ALGBRC(178)
      ALGBRC(214) = ALGBRC(189)+ALGBRC(207)
      ALGBRC(221) = ALGBRC(145)+ALGBRC(214)
      ALGBRC(256) = ALGBRC(132)+ALGBRC(250)+ALGBRC(255)+ALGBRC(136)+ALGBRC(140)+ALGBRC(227)+ALGBRC(244)+ALGBRC(240)+ALGBRC(221)+ALGBRC(233)+ALGBRC(231)+ALGBRC(229)+ALGBRC(12)
      ALGBRC(7) = STATES(2)+ 0.0183300*STATES(1)
      ALGBRC(18) =  CONSTS(125)*log(CONSTS(155)/ALGBRC(7))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END