Generated Code

The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

C
C There are a total of 253 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 174 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C CONSTS(144) is a1 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
C CONSTS(149) is a2 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
C CONSTS(1) is Faraday in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (coulomb_per_millimole).
C CONSTS(2) is Temp in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (kelvin).
C CONSTS(3) is Rgas in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (joule_per_mole_kelvin).
C CONSTS(116) is RT_over_F in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV).
C CONSTS(4) is Acap in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (cm2).
C CONSTS(124) is C in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mF).
C CONSTS(5) is Vmyo in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(6) is VJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(7) is VNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
C CONSTS(8) 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(129) is INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(226) is INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(241) is INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(237) is INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(204) is IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(133) is IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(137) is IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(224) is IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(252) is ICaK in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(247) is ICa in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(218) 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(142) is IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(186) is IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C ALGBRC(228) is ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(230) is IpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C ALGBRC(59) 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(89) is Jtrpn in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM_per_ms).
C ALGBRC(51) is Jrel in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(60) is Jtr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
C ALGBRC(99) is beta_SS in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(69) is beta_JSR in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(109) is beta_i in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
C ALGBRC(253) is i_tot in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uA_per_uF).
C CONSTS(9) is Ko in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(10) is Nao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(11) is Cao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
C CONSTS(12) is stim_period in component I_stimulus (ms).
C CONSTS(13) is stim_duration in component I_stimulus (ms).
C CONSTS(14) is stim_amplitude in component I_stimulus (uA_per_uF).
C CONSTS(15) 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(16) is Kfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(17) is Nfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C ALGBRC(34) is rb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C CONSTS(18) is Krb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(19) is Nrb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
C CONSTS(20) is KSR in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
C CONSTS(21) is vmaxf in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
C CONSTS(22) is vmaxr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
C CONSTS(23) 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(24) is tautr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
C CONSTS(25) is tauxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
C CONSTS(26) is LTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(27) is HTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(28) is khtrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
C CONSTS(29) is khtrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C CONSTS(30) is kltrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
C CONSTS(31) is kltrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C CONSTS(32) is CMDNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(33) is CSQNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(34) is EGTAtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(35) is KmCMDN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(36) is KmCSQN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C CONSTS(37) is KmEGTA in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
C ALGBRC(68) is dLTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
C ALGBRC(79) 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(38) is kaplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM4_per_ms).
C CONSTS(39) is kaminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(40) is kbplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM3_per_ms).
C CONSTS(41) is kbminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(42) is kcplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(43) is kcminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C CONSTS(44) is ncoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C CONSTS(45) is mcoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
C ALGBRC(2) is dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(13) is dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(24) is dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
C ALGBRC(35) 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(46) is fL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(47) is gL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(48) is bL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C CONSTS(49) is aL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
C ALGBRC(52) is C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(61) is C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(70) is C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(80) is C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(122) is C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(126) is C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(130) is C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(134) is C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(90) is CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(100) is CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(110) is CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(118) is CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(138) is CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(143) is CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(149) is CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(157) is CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(167) is C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(176) is C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(187) is C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(198) is C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(205) is C4_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(132) is CCa0_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(148) is CCa1_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(151) is CCa2_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(157) is CCa3_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(162) 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(44) is gamma in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C CONSTS(117) is omega in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(177) is a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(188) is a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(189) is a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(199) is a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(200) is a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(206) is a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(207) is a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(212) is a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(213) is a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(219) is a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(178) is a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(190) is a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(191) is a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(201) is a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(202) is a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(208) is a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(209) is a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(214) is a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(215) is a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
C ALGBRC(220) 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(118) 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(45) is C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(53) is C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(62) is C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(71) is C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(81) is CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(91) is CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(101) is CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(111) is CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(119) is C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(123) is C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(127) is C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(131) is OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(135) is CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(140) is CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(145) is CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(151) is OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(139) is C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(144) is C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(150) is C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(158) is C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(159) is OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(146) is CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(152) is CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(160) is CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(168) is CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(169) 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(153) is a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(161) is a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(162) is a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(170) is a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(171) is a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(179) is a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(180) is a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(192) is a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(181) is a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(193) is a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(154) is a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(163) is a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(164) is a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(172) is a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(173) is a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(182) is a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(183) is a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(194) is a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(184) is a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C ALGBRC(195) is a2_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(50) is alphaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(51) is aaKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(52) is betaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(53) is baKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(54) is alphai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(55) is aiKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(56) is betai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
C CONSTS(57) is biKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
C CONSTS(58) is f1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(59) is f2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(60) is f3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(61) is f4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(62) is b1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(63) is b2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(64) is b3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(65) is b4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
C CONSTS(66) is f1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(67) is f2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(68) is f3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(69) is f4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(70) is b1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(71) is b2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(72) is b3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(73) is b4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
C CONSTS(74) is alphaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(75) is aaKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
C CONSTS(76) is betaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(77) is baKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
C CONSTS(78) is alphai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(79) 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(119) is alpha_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(120) 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(46) is C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(54) 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(47) is CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(55) is CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(63) is C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(72) is C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(82) is C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(92) is OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(64) is CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(73) is CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(83) is CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(93) is OIKv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(133) is C0Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(134) is C1Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(135) is C2Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(136) is C3Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(137) is OKv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(138) is CI0Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(139) is CI1Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(140) is CI2Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(141) is CI3Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C CONSTS(142) is OIKv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(74) is a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(84) 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(85) is a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(94) 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(95) is a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(102) 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(103) is a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(112) 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(104) is a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(113) 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(75) is a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(86) 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(87) is a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(96) 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(97) is a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(105) 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(106) is a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(114) 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(107) is a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
C ALGBRC(115) 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(117) is ENa in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(121) is EK in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(18) is EKs in component COMPUTE_REVERSAL_POTENTIALS (mV).
C ALGBRC(125) is ECa in component COMPUTE_REVERSAL_POTENTIALS (mV).
C CONSTS(143) is a1 in component COMPUTE_REVERSAL_POTENTIALS (mM).
C ALGBRC(7) is a2 in component COMPUTE_REVERSAL_POTENTIALS (mM).
C CONSTS(80) is GKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(81) is GKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(82) is GK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(83) is GNa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(84) is GNab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(85) is KvScale in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C CONSTS(86) is Kv43Frac in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C CONSTS(121) is GKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
C CONSTS(122) 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(123) is fKo in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(148) is VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(156) is VFsq_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (coulomb_per_millimole).
C ALGBRC(166) is a1_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
C ALGBRC(175) is a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(197) is a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
C ALGBRC(222) is K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
C ALGBRC(211) is IKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
C CONSTS(87) is kNaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(88) is KmNa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(89) is KmCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(90) is ksat in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(91) is eta in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(92) is INaKmax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(93) is KmNai in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(94) is KmKo in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(95) is IpCamax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
C CONSTS(96) is KmpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(97) is GCab in component COMPUTE_INaK_INaCa_ICab_IpCa (mS_per_uF).
C ALGBRC(232) is VF_over_RT in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(152) is sigma in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(233) is a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(234) is a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(235) is fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(156) is a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(236) is a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C ALGBRC(238) is a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
C ALGBRC(239) is a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
C ALGBRC(240) is a3_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
C CONSTS(161) is a4_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
C CONSTS(164) is a5_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM3).
C ALGBRC(249) is PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C ALGBRC(242) is VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(243) is VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole).
C ALGBRC(244) is a1_Ca in component COMPUTE_ICa_ICaK (mM).
C ALGBRC(245) is a2_Ca in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(250) is a1_K in component COMPUTE_ICa_ICaK (mM).
C ALGBRC(251) is a2_K in component COMPUTE_ICa_ICaK (dimensionless).
C ALGBRC(246) is ICamax in component COMPUTE_ICa_ICaK (uA_per_uF).
C ALGBRC(248) is Icabar in component COMPUTE_ICa_ICaK (uA_per_uF).
C CONSTS(145) is PCa in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C CONSTS(125) is PK in component COMPUTE_ICa_ICaK (litre_per_farad_second).
C CONSTS(98) is ICahalf in component COMPUTE_ICa_ICaK (uA_per_uF).
C CONSTS(99) is Pscale in component COMPUTE_ICa_ICaK (dimensionless).
C CONSTS(150) 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(48) is On in component INa (per_ms).
C ALGBRC(56) is Of in component INa (per_ms).
C ALGBRC(65) is GammaGamma in component INa (per_ms).
C ALGBRC(76) is DeltaDelta in component INa (per_ms).
C CONSTS(154) is epsilon in component INa (per_ms).
C CONSTS(155) is omega_na in component INa (per_ms).
C ALGBRC(88) is rho in component INa (per_ms).
C ALGBRC(98) is mu in component INa (per_ms).
C CONSTS(153) is Cn in component INa (per_ms).
C CONSTS(158) is Cf in component INa (per_ms).
C CONSTS(163) is parameter_a in component INa (dimensionless).
C ALGBRC(108) is k12 in component INa (per_ms).
C ALGBRC(116) is k23 in component INa (per_ms).
C ALGBRC(120) is k34 in component INa (per_ms).
C ALGBRC(124) is k45 in component INa (per_ms).
C ALGBRC(128) is k56 in component INa (per_ms).
C CONSTS(159) is k67 in component INa (per_ms).
C ALGBRC(132) is k89 in component INa (per_ms).
C ALGBRC(136) is k910 in component INa (per_ms).
C ALGBRC(141) is k1011 in component INa (per_ms).
C ALGBRC(147) is k1112 in component INa (per_ms).
C ALGBRC(155) is k1213 in component INa (per_ms).
C ALGBRC(165) is k57 in component INa (per_ms).
C ALGBRC(174) is k21 in component INa (per_ms).
C ALGBRC(185) is k32 in component INa (per_ms).
C ALGBRC(196) is k43 in component INa (per_ms).
C ALGBRC(203) is k54 in component INa (per_ms).
C ALGBRC(210) is k65 in component INa (per_ms).
C CONSTS(160) is k76 in component INa (per_ms).
C ALGBRC(216) is k98 in component INa (per_ms).
C ALGBRC(221) is k109 in component INa (per_ms).
C ALGBRC(223) is k1110 in component INa (per_ms).
C ALGBRC(225) is k1211 in component INa (per_ms).
C ALGBRC(227) is k1312 in component INa (per_ms).
C ALGBRC(217) is k75 in component INa (per_ms).
C CONSTS(165) is k81 in component INa (per_ms).
C CONSTS(167) is k92 in component INa (per_ms).
C CONSTS(169) is k103 in component INa (per_ms).
C CONSTS(171) is k114 in component INa (per_ms).
C CONSTS(173) is k125 in component INa (per_ms).
C ALGBRC(229) is k136 in component INa (per_ms).
C CONSTS(166) is k18 in component INa (per_ms).
C CONSTS(168) is k29 in component INa (per_ms).
C CONSTS(170) is k310 in component INa (per_ms).
C CONSTS(172) is k411 in component INa (per_ms).
C CONSTS(174) is k512 in component INa (per_ms).
C ALGBRC(231) 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(100) is TNa in component INa (kelvin).
C CONSTS(126) is KToverH in component INa (per_ms).
C CONSTS(147) is FoverRT in component INa (per_mV).
C CONSTS(146) is RTNa in component INa (joule_per_mole).
C CONSTS(127) is RTNaF in component INa (mV).
C CONSTS(101) is T_Const_HERG in component IKr (dimensionless).
C CONSTS(102) is A0_HERG in component IKr (per_ms).
C CONSTS(103) is B0_HERG in component IKr (per_mV).
C CONSTS(104) is A1_HERG in component IKr (per_ms).
C CONSTS(105) is B1_HERG in component IKr (per_mV).
C CONSTS(106) is A2_HERG in component IKr (per_ms).
C CONSTS(107) is B2_HERG in component IKr (per_mV).
C CONSTS(108) is A3_HERG in component IKr (per_ms).
C CONSTS(109) is B3_HERG in component IKr (per_mV).
C CONSTS(110) is A4_HERG in component IKr (per_ms).
C CONSTS(111) is B4_HERG in component IKr (per_mV).
C CONSTS(112) is A5_HERG in component IKr (per_ms).
C CONSTS(113) is B5_HERG in component IKr (per_mV).
C CONSTS(114) is A6_HERG in component IKr (per_ms).
C CONSTS(115) is B6_HERG in component IKr (per_mV).
C CONSTS(128) is C2H_to_C3H in component IKr (per_ms).
C CONSTS(129) 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(49) is C3H_to_IH in component IKr (per_ms).
C ALGBRC(57) 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(66) is a1_C3 in component IKr (per_ms).
C ALGBRC(77) is a2_C3 in component IKr (per_ms).
C STATES(64) is C3Herg in component IKr (dimensionless).
C ALGBRC(50) is a1_O in component IKr (per_ms).
C ALGBRC(58) is a2_O in component IKr (per_ms).
C ALGBRC(67) is a1_I in component IKr (per_ms).
C ALGBRC(78) is a2_I in component IKr (per_ms).
C STATES(65) is IHerg in component IKr (dimensionless).
C CONSTS(130) is C0ks_C1ks in component IKs (per_ms).
C CONSTS(131) 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) = 96.5
      CONSTS(2) = 310
      CONSTS(3) = 8.315
      CONSTS(4) = 1.534e-4
      CONSTS(5) = 25.84e-6
      CONSTS(6) = 0.16e-6
      CONSTS(7) = 2.1e-6
      CONSTS(8) = 1.2e-9
      STATES(1) = 0.9798304162e1
      STATES(2) = 0.1255589432e3
      STATES(3) = 0.8601192016e-4
      STATES(4) = 0.1420215245e-3
      STATES(5) = 0.2852239446
      STATES(6) = 0.2855294915
      STATES(7) = -0.9065755929e2
      CONSTS(9) = 4
      CONSTS(10) = 138
      CONSTS(11) = 2
      CONSTS(12) = 1e3
      CONSTS(13) = 3
      CONSTS(14) = -15
      CONSTS(15) = 0
      CONSTS(16) = 0.000168
      CONSTS(17) = 1.2
      CONSTS(18) = 3.29
      CONSTS(19) = 1
      CONSTS(20) = 1.2
      CONSTS(21) = 0.0748e-3
      CONSTS(22) = 0.318e-3
      CONSTS(23) = 1.8
      STATES(8) = 0.6601783287e-3
      STATES(9) = 0.397392006e-8
      CONSTS(24) = 0.5747
      CONSTS(25) = 26.7
      CONSTS(26) = 70e-3
      CONSTS(27) = 140e-3
      CONSTS(28) = 20e0
      CONSTS(29) = 0.066e-3
      CONSTS(30) = 40e0
      CONSTS(31) = 40e-3
      CONSTS(32) = 50e-3
      CONSTS(33) = 15e0
      CONSTS(34) = 0e0
      CONSTS(35) = 2.38e-3
      CONSTS(36) = 0.8e0
      CONSTS(37) = 1.5e-4
      STATES(10) = 0.8046584973e-1
      STATES(11) = 0.9772152481
      CONSTS(38) = 0.01215
      CONSTS(39) = 0.576
      CONSTS(40) = 0.00405
      CONSTS(41) = 1.93
      CONSTS(42) = 0.1
      CONSTS(43) = 0.0008
      CONSTS(44) = 4
      CONSTS(45) = 3
      STATES(12) = 0.4448038946
      STATES(13) = 0.554535924
      CONSTS(46) = 0.3
      CONSTS(47) = 4
      CONSTS(48) = 2
      CONSTS(49) = 2
      STATES(14) = 0.8621934054
      STATES(15) = 0.01141398211
      STATES(16) = 5.666325225e-5
      STATES(17) = 1.25021357e-7
      STATES(18) = 1.034426748e-10
      STATES(19) = 0.1198547081
      STATES(20) = 0.006346794302
      STATES(21) = 1.260326488e-4
      STATES(22) = 1.112315238e-6
      STATES(23) = 3.681315892e-9
      STATES(24) = 7.757981563e-12
      STATES(25) = 0.9997157074
      STATES(26) = 0.9513721351
      STATES(27) = 0.02668288089
      STATES(28) = 2.806380358e-4
      STATES(29) = 1.311837579e-6
      STATES(30) = 2.29955624e-9
      STATES(31) = 0.01513025204
      STATES(32) = 0.005442964601
      STATES(33) = 9.918373359000001e-4
      STATES(34) = 9.514386057000001e-5
      STATES(35) = 2.742677382e-6
      CONSTS(50) = 0.543708
      CONSTS(51) = 0.028983
      CONSTS(52) = 0.080185
      CONSTS(53) = 0.0468437
      CONSTS(54) = 0.0498424
      CONSTS(55) = 0.000373016
      CONSTS(56) = 0.000819482
      CONSTS(57) = 0.00000005374
      CONSTS(58) = 1.8936
      CONSTS(59) = 14.224647456
      CONSTS(60) = 158.574378389
      CONSTS(61) = 142.936645351
      CONSTS(62) = 6.77348
      CONSTS(63) = 15.6212705152
      CONSTS(64) = 28.7532603313
      CONSTS(65) = 524.576206679
      CONSTS(66) = 0.52465073996
      CONSTS(67) = 17.51885408639
      CONSTS(68) = 938.58764534556
      CONSTS(69) = 54749.19473332601
      CONSTS(70) = 1.00947847105
      CONSTS(71) = 1.17100540567
      CONSTS(72) = 0.63902768758
      CONSTS(73) = 2.12035379095
      CONSTS(74) = 1.84002414554
      CONSTS(75) = 0.00768548031
      CONSTS(76) = 0.0108174834
      CONSTS(77) = 0.07793378174
      CONSTS(78) = 0.00305767916
      CONSTS(79) = 0.00000244936
      STATES(36) = 0.5977099765
      STATES(37) = 0.1730990528
      STATES(38) = 0.01881072386
      STATES(39) = 9.16070135e-4
      STATES(40) = 1.975541357e-5
      STATES(41) = 0.03539084346
      STATES(42) = 0.005428824353
      STATES(43) = 0.02287858869
      STATES(44) = 0.03233800003
      STATES(45) = 0.1134082058
      CONSTS(80) = 0.0186e0
      CONSTS(81) = 0.0035e0
      CONSTS(82) = 0.1253051261188079722568586886469e0
      CONSTS(83) = 56.32e0
      CONSTS(84) = 0.001e0
      CONSTS(85) = 0.872
      CONSTS(86) = 0.889
      STATES(46) = 1.57442749e-8
      STATES(47) = 2.856607179e-9
      STATES(48) = 9.511789113e-6
      STATES(49) = 2.491710696e-7
      STATES(50) = 1.298547822e-5
      CONSTS(87) = 0.44
      CONSTS(88) = 87.5
      CONSTS(89) = 1.38
      CONSTS(90) = 0.2
      CONSTS(91) = 0.35
      CONSTS(92) = 2.387
      CONSTS(93) = 20
      CONSTS(94) = 1.5
      CONSTS(95) = 0.05
      CONSTS(96) = 0.0005
      CONSTS(97) = 7.684e-5
      CONSTS(98) = -0.265
      CONSTS(99) = 7
      STATES(51) = 0.1437575649
      STATES(52) = 0.0417776208
      STATES(53) = 0.004552898364
      STATES(54) = 2.20520743e-4
      STATES(55) = 4.005266484e-6
      STATES(56) = 0.4749950008
      STATES(57) = 0.2707214097
      STATES(58) = 0.05786120057
      STATES(59) = 0.005496292279
      STATES(60) = 1.957874519e-4
      STATES(61) = 4.17679396e-4
      CONSTS(100) = 294.16
      CONSTS(101) = 5.320000001
      CONSTS(102) = 0.017147641733086
      CONSTS(103) = 0.03304608038835
      CONSTS(104) = 0.03969328381141
      CONSTS(105) = -0.0430605416398
      CONSTS(106) = 0.02057448605977
      CONSTS(107) = 0.02617412715118
      CONSTS(108) = 0.00134366604423
      CONSTS(109) = -0.02691385498399
      CONSTS(110) = 0.10666316491288
      CONSTS(111) = 0.00568908859717
      CONSTS(112) = 0.00646393910049
      CONSTS(113) = -0.04536642959543
      CONSTS(114) = 0.00008039374403
      CONSTS(115) = 0.00000069808924
      STATES(62) = 0.996697338
      STATES(63) = 4.340879648e-4
      STATES(64) = 7.634099755000001e-5
      STATES(65) = 1.533347007e-6
      STATES(66) = 0.9645606295
      STATES(67) = 0.03542613568
      CONSTS(116) = ( CONSTS(3)*CONSTS(2))/CONSTS(1)
      CONSTS(117) =  0.250000*0.0100000
      CONSTS(118) = 0.820000
      CONSTS(119) = CONSTS(78)
      CONSTS(120) = CONSTS(79)
      CONSTS(121) =  CONSTS(86)*CONSTS(85)*0.100000
      CONSTS(122) =  (1.00000 - CONSTS(86))*CONSTS(85)*4.29860e-07
      CONSTS(123) =  (CONSTS(9)/4.00000) ** (1.0 / 2)
      CONSTS(124) =  CONSTS(4)*0.00100000
      CONSTS(125) =  CONSTS(99)*4.57400e-07
      CONSTS(126) = ( 1.38100e-23*CONSTS(100))/6.62600e-31
      CONSTS(127) = ( CONSTS(3)*CONSTS(100))/CONSTS(1)
      CONSTS(128) =  CONSTS(101)*0.0260836
      CONSTS(129) =  CONSTS(101)*0.148330
      CONSTS(130) = 0.00795601
      CONSTS(131) = 0.0396672
      CONSTS(132) = CONSTS(117)
      CONSTS(133) = CONSTS(120)
      CONSTS(134) =  CONSTS(66)*CONSTS(120)
      CONSTS(135) =  CONSTS(67)*CONSTS(120)
      CONSTS(136) =  CONSTS(68)*CONSTS(120)
      CONSTS(137) =  CONSTS(69)*CONSTS(120)
      CONSTS(138) = CONSTS(119)
      CONSTS(139) = CONSTS(119)/CONSTS(70)
      CONSTS(140) = CONSTS(119)/CONSTS(71)
      CONSTS(141) = CONSTS(119)/CONSTS(72)
      CONSTS(142) = CONSTS(119)/CONSTS(73)
      CONSTS(143) = CONSTS(9)+ 0.0183300*CONSTS(10)
      CONSTS(144) = CONSTS(124)/( CONSTS(5)*CONSTS(1))
      CONSTS(145) =  CONSTS(99)*0.000246900
      CONSTS(146) =  CONSTS(3)*CONSTS(100)
      CONSTS(147) = 1.00000/CONSTS(116)
      CONSTS(148) = CONSTS(132)/CONSTS(48)
      CONSTS(149) = CONSTS(124)/( 2.00000*CONSTS(8)*CONSTS(1))
      CONSTS(150) = 1.38862
      CONSTS(151) = CONSTS(148)/CONSTS(48)
      CONSTS(152) = (EXP(CONSTS(10)/67.3000) - 1.00000)/7.00000
      CONSTS(153) =  CONSTS(150)*CONSTS(126)*EXP(- 287913./CONSTS(146)+786.217/CONSTS(3))
      CONSTS(154) =  CONSTS(150)*CONSTS(126)*EXP(- 85800.4/CONSTS(146)+70.0780/CONSTS(3))
      CONSTS(155) =  CONSTS(150)*CONSTS(126)*EXP(- 121955./CONSTS(146)+225.175/CONSTS(3))
      CONSTS(156) = CONSTS(9)/(CONSTS(9)+CONSTS(94))
      CONSTS(157) = CONSTS(151)/CONSTS(48)
      CONSTS(158) =  CONSTS(150)*CONSTS(126)*EXP(- 59565.2/CONSTS(146)+0.00711000/CONSTS(3))
      CONSTS(159) = CONSTS(154)
      CONSTS(160) = CONSTS(155)
      CONSTS(161) = CONSTS(89)+CONSTS(11)
      CONSTS(162) = CONSTS(157)/CONSTS(48)
      CONSTS(163) = 1.40043
      CONSTS(164) = (CONSTS(88) ** 3.00000+CONSTS(10) ** 3.00000)/5000.00
      CONSTS(165) = CONSTS(158)
      CONSTS(166) = CONSTS(153)
      CONSTS(167) = CONSTS(165)/CONSTS(163)
      CONSTS(168) =  CONSTS(166)*CONSTS(163)
      CONSTS(169) = CONSTS(167)/CONSTS(163)
      CONSTS(170) =  CONSTS(168)*CONSTS(163)
      CONSTS(171) = CONSTS(169)/CONSTS(163)
      CONSTS(172) =  CONSTS(170)*CONSTS(163)
      CONSTS(173) = CONSTS(171)/CONSTS(163)
      CONSTS(174) =  CONSTS(172)*CONSTS(163)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(24) =  CONSTS(46)*STATES(18) -  CONSTS(47)*STATES(24)
      ALGBRC(2) =  - CONSTS(38)* STATES(4)*1000.00 ** CONSTS(44)*STATES(12)+ CONSTS(39)*STATES(8)
      RATES(12) = ALGBRC(2)
      ALGBRC(13) =  CONSTS(40)* STATES(4)*1000.00 ** CONSTS(45)*STATES(8) -  CONSTS(41)*STATES(9)
      RATES(9) = ALGBRC(13)
      ALGBRC(4) = (CONSTS(118)/(1.00000+EXP((STATES(7)+28.5000)/7.80000))+1.00000) - CONSTS(118)
      ALGBRC(15) = 1.00000/(0.00653000/(0.500000+EXP(STATES(7)/- 7.10000))+ 0.00512000*EXP(STATES(7)/- 39.8000))
      RATES(25) = (ALGBRC(4) - STATES(25))/ALGBRC(15)
      ALGBRC(9) =  CONSTS(101)*CONSTS(102)*EXP( CONSTS(103)*STATES(7))
      ALGBRC(20) =  CONSTS(101)*CONSTS(104)*EXP( CONSTS(105)*STATES(7))
      RATES(62) =  ALGBRC(20)*STATES(63) -  ALGBRC(9)*STATES(62)
      ALGBRC(24) =  CONSTS(42)*STATES(8) -  CONSTS(43)*STATES(13)
      RATES(13) = ALGBRC(24)
      ALGBRC(33) =  0.216256*EXP( - 1.88912e-05*STATES(7))
      RATES(66) =  - CONSTS(130)*STATES(66)+ ALGBRC(33)*STATES(67)
      ALGBRC(22) =  0.00700807*EXP( - 0.149998*STATES(7))
      RATES(67) = ( CONSTS(130)*STATES(66) -  (ALGBRC(33)+CONSTS(131))*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(131)*STATES(67) -  (ALGBRC(22)+ALGBRC(11))*STATES(49))+ ALGBRC(32)*STATES(50)
      RATES(50) =  ALGBRC(11)*STATES(49) -  ALGBRC(32)*STATES(50)
      ALGBRC(35) = - (ALGBRC(2)+ALGBRC(13)+ALGBRC(24))
      RATES(8) = ALGBRC(35)
      ALGBRC(31) =  ALGBRC(9)*STATES(62)+ CONSTS(129)*STATES(64)
      ALGBRC(42) =  (ALGBRC(20)+CONSTS(128))*STATES(63)
      RATES(63) = ALGBRC(31) - ALGBRC(42)
      ALGBRC(10) =  CONSTS(101)*CONSTS(106)*EXP( CONSTS(107)*STATES(7))
      ALGBRC(41) =  CONSTS(101)*CONSTS(112)*EXP( CONSTS(113)*STATES(7))
      ALGBRC(50) =  ALGBRC(10)*STATES(64)+ ALGBRC(41)*STATES(65)
      ALGBRC(21) =  CONSTS(101)*CONSTS(108)*EXP( CONSTS(109)*STATES(7))
      ALGBRC(30) =  CONSTS(101)*CONSTS(110)*EXP( CONSTS(111)*STATES(7))
      ALGBRC(58) =  (ALGBRC(21)+ALGBRC(30))*STATES(48)
      RATES(48) = ALGBRC(50) - ALGBRC(58)
      ALGBRC(23) = STATES(3)/CONSTS(16) ** CONSTS(17)
      ALGBRC(34) = STATES(6)/CONSTS(18) ** CONSTS(19)
      ALGBRC(43) = ( CONSTS(20)*( CONSTS(21)*ALGBRC(23) -  CONSTS(22)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34))
      ALGBRC(60) = (STATES(6) - STATES(5))/CONSTS(24)
      RATES(6) = ( ALGBRC(43)*CONSTS(5))/CONSTS(7) - ( ALGBRC(60)*CONSTS(6))/CONSTS(7)
      ALGBRC(51) =  CONSTS(23)*(STATES(8)+STATES(9))*(STATES(5) - STATES(4))
      ALGBRC(69) = 1.00000/(1.00000+( CONSTS(33)*CONSTS(36))/STATES(5)+CONSTS(36) ** 2.00000)
      RATES(5) =  ALGBRC(69)*(ALGBRC(60) - ALGBRC(51))
      ALGBRC(68) =  CONSTS(30)*STATES(3)*(1.00000 - STATES(10)) -  CONSTS(31)*STATES(10)
      RATES(10) = ALGBRC(68)
      ALGBRC(49) =  CONSTS(101)*CONSTS(114)*EXP( CONSTS(115)*STATES(7))
      ALGBRC(57) = ( ALGBRC(21)*ALGBRC(41)*ALGBRC(49))/( ALGBRC(10)*ALGBRC(30))
      ALGBRC(66) =  CONSTS(128)*STATES(63)+ ALGBRC(21)*STATES(48)+ ALGBRC(57)*STATES(65)
      ALGBRC(77) =  (ALGBRC(49)+ALGBRC(10)+CONSTS(129))*STATES(64)
      RATES(64) = ALGBRC(66) - ALGBRC(77)
      ALGBRC(67) =  ALGBRC(49)*STATES(64)+ ALGBRC(30)*STATES(48)
      ALGBRC(78) =  (ALGBRC(57)+ALGBRC(41))*STATES(65)
      RATES(65) = ALGBRC(67) - ALGBRC(78)
      ALGBRC(79) =  CONSTS(28)*STATES(3)*(1.00000 - STATES(11)) -  CONSTS(29)*STATES(11)
      RATES(11) = ALGBRC(79)
      ALGBRC(6) =  CONSTS(74)*EXP( CONSTS(75)*STATES(7))
      ALGBRC(27) =  4.00000*ALGBRC(6)
      ALGBRC(74) =  (ALGBRC(27)+CONSTS(133))*STATES(36)
      ALGBRC(17) =  CONSTS(76)*EXP( - CONSTS(77)*STATES(7))
      ALGBRC(63) = ALGBRC(17)
      ALGBRC(84) =  ALGBRC(63)*STATES(37)+ CONSTS(138)*STATES(41)
      RATES(36) = ALGBRC(84) - ALGBRC(74)
      ALGBRC(28) =  4.00000*CONSTS(70)*ALGBRC(6)
      ALGBRC(75) =  (CONSTS(138)+ALGBRC(28))*STATES(41)
      ALGBRC(64) = ALGBRC(17)/CONSTS(66)
      ALGBRC(86) =  CONSTS(133)*STATES(36)+ ALGBRC(64)*STATES(42)
      RATES(41) = ALGBRC(86) - ALGBRC(75)
      ALGBRC(38) =  3.00000*ALGBRC(6)
      ALGBRC(85) =  (ALGBRC(38)+ALGBRC(63)+CONSTS(134))*STATES(37)
      ALGBRC(72) =  2.00000*ALGBRC(17)
      ALGBRC(94) =  ALGBRC(72)*STATES(38)+ CONSTS(139)*STATES(42)+ ALGBRC(27)*STATES(36)
      RATES(37) = ALGBRC(94) - ALGBRC(85)
      ALGBRC(39) = ( 3.00000*CONSTS(71)*ALGBRC(6))/CONSTS(70)
      ALGBRC(87) =  (ALGBRC(39)+CONSTS(139)+ALGBRC(64))*STATES(42)
      ALGBRC(73) = ( 2.00000*CONSTS(66)*ALGBRC(17))/CONSTS(67)
      ALGBRC(96) =  ALGBRC(73)*STATES(43)+ CONSTS(134)*STATES(37)+ ALGBRC(28)*STATES(41)
      RATES(42) = ALGBRC(96) - ALGBRC(87)
      ALGBRC(46) =  2.00000*ALGBRC(6)
      ALGBRC(95) =  (ALGBRC(46)+ALGBRC(72)+CONSTS(135))*STATES(38)
      ALGBRC(82) =  3.00000*ALGBRC(17)
      ALGBRC(102) =  ALGBRC(82)*STATES(39)+ CONSTS(140)*STATES(43)+ ALGBRC(38)*STATES(37)
      RATES(38) = ALGBRC(102) - ALGBRC(95)
      ALGBRC(47) = ( 2.00000*CONSTS(72)*ALGBRC(6))/CONSTS(71)
      ALGBRC(97) =  (ALGBRC(47)+CONSTS(140)+ALGBRC(73))*STATES(43)
      ALGBRC(83) = ( 3.00000*CONSTS(67)*ALGBRC(17))/CONSTS(68)
      ALGBRC(105) =  ALGBRC(83)*STATES(44)+ CONSTS(135)*STATES(38)+ ALGBRC(39)*STATES(42)
      RATES(43) = ALGBRC(105) - ALGBRC(97)
      ALGBRC(54) = ALGBRC(6)
      ALGBRC(103) =  (ALGBRC(54)+ALGBRC(82)+CONSTS(136))*STATES(39)
      ALGBRC(92) =  4.00000*ALGBRC(17)
      ALGBRC(112) =  ALGBRC(92)*STATES(40)+ CONSTS(141)*STATES(44)+ ALGBRC(46)*STATES(38)
      RATES(39) = ALGBRC(112) - ALGBRC(103)
      ALGBRC(104) =  (ALGBRC(92)+CONSTS(137))*STATES(40)
      ALGBRC(113) =  ALGBRC(54)*STATES(39)+ CONSTS(142)*STATES(45)
      RATES(40) = ALGBRC(113) - ALGBRC(104)
      ALGBRC(55) = ( CONSTS(73)*ALGBRC(6))/CONSTS(72)
      ALGBRC(106) =  (ALGBRC(55)+CONSTS(141)+ALGBRC(83))*STATES(44)
      ALGBRC(93) = ( 4.00000*CONSTS(68)*ALGBRC(17))/CONSTS(69)
      ALGBRC(114) =  ALGBRC(93)*STATES(45)+ CONSTS(136)*STATES(39)+ ALGBRC(47)*STATES(43)
      RATES(44) = ALGBRC(114) - ALGBRC(106)
      ALGBRC(107) =  (CONSTS(142)+ALGBRC(93))*STATES(45)
      ALGBRC(115) =  CONSTS(137)*STATES(40)+ ALGBRC(55)*STATES(44)
      RATES(45) = ALGBRC(115) - ALGBRC(107)
      ALGBRC(5) =  CONSTS(50)*EXP( CONSTS(51)*STATES(7))
      ALGBRC(45) =  4.00000*ALGBRC(5)
      ALGBRC(37) =  CONSTS(56)*EXP( CONSTS(57)*STATES(7))
      ALGBRC(139) = ALGBRC(37)
      ALGBRC(153) =  (ALGBRC(45)+ALGBRC(139))*STATES(26)
      ALGBRC(16) =  CONSTS(52)*EXP( - CONSTS(53)*STATES(7))
      ALGBRC(119) = ALGBRC(16)
      ALGBRC(26) =  CONSTS(54)*EXP( - CONSTS(55)*STATES(7))
      ALGBRC(146) = ALGBRC(26)
      ALGBRC(161) =  ALGBRC(119)*STATES(27)+ ALGBRC(146)*STATES(31)
      RATES(26) = ALGBRC(161) - ALGBRC(153)
      ALGBRC(81) =  4.00000*CONSTS(62)*ALGBRC(5)
      ALGBRC(154) =  (ALGBRC(146)+ALGBRC(81))*STATES(31)
      ALGBRC(135) = ALGBRC(16)/CONSTS(58)
      ALGBRC(163) =  ALGBRC(139)*STATES(26)+ ALGBRC(135)*STATES(32)
      RATES(31) = ALGBRC(163) - ALGBRC(154)
      ALGBRC(53) =  3.00000*ALGBRC(5)
      ALGBRC(144) =  CONSTS(58)*ALGBRC(37)
      ALGBRC(162) =  (ALGBRC(53)+ALGBRC(119)+ALGBRC(144))*STATES(27)
      ALGBRC(123) =  2.00000*ALGBRC(16)
      ALGBRC(152) = ALGBRC(26)/CONSTS(62)
      ALGBRC(170) =  ALGBRC(123)*STATES(28)+ ALGBRC(152)*STATES(32)+ ALGBRC(45)*STATES(26)
      RATES(27) = ALGBRC(170) - ALGBRC(162)
      ALGBRC(91) = ( 3.00000*CONSTS(63)*ALGBRC(5))/CONSTS(62)
      ALGBRC(164) =  (ALGBRC(91)+ALGBRC(152)+ALGBRC(135))*STATES(32)
      ALGBRC(140) = ( 2.00000*CONSTS(58)*ALGBRC(16))/CONSTS(59)
      ALGBRC(172) =  ALGBRC(140)*STATES(33)+ ALGBRC(144)*STATES(27)+ ALGBRC(81)*STATES(31)
      RATES(32) = ALGBRC(172) - ALGBRC(164)
      ALGBRC(8) =  CONSTS(150)*CONSTS(126)*EXP(- 114007./CONSTS(146)+224.114/CONSTS(3)+( 0.286374*STATES(7))/CONSTS(127))
      ALGBRC(108) =  4.00000*ALGBRC(8)
      ALGBRC(19) =  CONSTS(150)*CONSTS(126)*EXP(- 272470./CONSTS(146)+708.146/CONSTS(3)+( - 2.28528*STATES(7))/CONSTS(127))
      ALGBRC(174) = ALGBRC(19)
      RATES(51) =  - (CONSTS(166)+ALGBRC(108))*STATES(51)+ ALGBRC(174)*STATES(52)+ CONSTS(165)*STATES(56)
      ALGBRC(62) =  2.00000*ALGBRC(5)
      ALGBRC(150) =  CONSTS(59)*ALGBRC(37)
      ALGBRC(171) =  (ALGBRC(62)+ALGBRC(123)+ALGBRC(150))*STATES(28)
      ALGBRC(127) =  3.00000*ALGBRC(16)
      ALGBRC(160) = ALGBRC(26)/CONSTS(63)
      ALGBRC(179) =  ALGBRC(127)*STATES(29)+ ALGBRC(160)*STATES(33)+ ALGBRC(53)*STATES(27)
      RATES(28) = ALGBRC(179) - ALGBRC(171)
      ALGBRC(101) = ( 2.00000*CONSTS(64)*ALGBRC(5))/CONSTS(63)
      ALGBRC(173) =  (ALGBRC(101)+ALGBRC(160)+ALGBRC(140))*STATES(33)
      ALGBRC(145) = ( 3.00000*CONSTS(59)*ALGBRC(16))/CONSTS(60)
      ALGBRC(182) =  ALGBRC(145)*STATES(34)+ ALGBRC(150)*STATES(28)+ ALGBRC(91)*STATES(32)
      RATES(33) = ALGBRC(182) - ALGBRC(173)
      ALGBRC(116) =  3.00000*ALGBRC(8)
      ALGBRC(185) =  2.00000*ALGBRC(19)
      RATES(52) = ( ALGBRC(108)*STATES(51) -  (ALGBRC(174)+ALGBRC(116)+CONSTS(168))*STATES(52))+ ALGBRC(185)*STATES(53)+ CONSTS(167)*STATES(57)
      ALGBRC(3) =  4.00000*1.20000*0.416000*EXP( 0.0120000*(STATES(7) - 35.0000))
      ALGBRC(52) =  4.00000*ALGBRC(3)
      ALGBRC(44) =  0.600000*0.0923300*STATES(4)
      ALGBRC(167) = ALGBRC(44)
      ALGBRC(178) =  (ALGBRC(52)+ALGBRC(167))*STATES(14)
      ALGBRC(14) =  4.00000*0.450000*0.0490000*EXP( - 0.0650000*(STATES(7) - 22.0000))
      ALGBRC(122) = ALGBRC(14)
      ALGBRC(190) =  ALGBRC(122)*STATES(15)+ CONSTS(132)*STATES(19)
      RATES(14) = ALGBRC(190) - ALGBRC(178)
      ALGBRC(25) =  CONSTS(49)*ALGBRC(3)
      ALGBRC(90) =  4.00000*ALGBRC(25)
      ALGBRC(177) =  (ALGBRC(90)+CONSTS(132))*STATES(19)
      ALGBRC(36) = ALGBRC(14)/CONSTS(48)
      ALGBRC(138) = ALGBRC(36)
      ALGBRC(188) =  ALGBRC(138)*STATES(20)+ ALGBRC(167)*STATES(14)
      RATES(19) = ALGBRC(188) - ALGBRC(177)
      ALGBRC(71) = ALGBRC(5)
      ALGBRC(158) =  CONSTS(60)*ALGBRC(37)
      ALGBRC(180) =  (ALGBRC(71)+ALGBRC(127)+ALGBRC(158))*STATES(29)
      ALGBRC(131) =  4.00000*ALGBRC(16)
      ALGBRC(168) = ALGBRC(26)/CONSTS(64)
      ALGBRC(192) =  ALGBRC(131)*STATES(30)+ ALGBRC(168)*STATES(34)+ ALGBRC(62)*STATES(28)
      RATES(29) = ALGBRC(192) - ALGBRC(180)
      ALGBRC(159) =  CONSTS(61)*ALGBRC(37)
      ALGBRC(181) =  (ALGBRC(131)+ALGBRC(159))*STATES(30)
      ALGBRC(169) = ALGBRC(26)/CONSTS(65)
      ALGBRC(193) =  ALGBRC(71)*STATES(29)+ ALGBRC(169)*STATES(35)
      RATES(30) = ALGBRC(193) - ALGBRC(181)
      ALGBRC(111) = ( CONSTS(65)*ALGBRC(5))/CONSTS(64)
      ALGBRC(183) =  (ALGBRC(111)+ALGBRC(168)+ALGBRC(145))*STATES(34)
      ALGBRC(151) = ( 4.00000*CONSTS(60)*ALGBRC(16))/CONSTS(61)
      ALGBRC(194) =  ALGBRC(151)*STATES(35)+ ALGBRC(158)*STATES(29)+ ALGBRC(101)*STATES(33)
      RATES(34) = ALGBRC(194) - ALGBRC(183)
      ALGBRC(184) =  (ALGBRC(169)+ALGBRC(151))*STATES(35)
      ALGBRC(195) =  ALGBRC(159)*STATES(30)+ ALGBRC(111)*STATES(34)
      RATES(35) = ALGBRC(195) - ALGBRC(184)
      ALGBRC(120) =  2.00000*ALGBRC(8)
      ALGBRC(196) =  3.00000*ALGBRC(19)
      RATES(53) = ( ALGBRC(116)*STATES(52) -  (ALGBRC(185)+ALGBRC(120)+CONSTS(170))*STATES(53))+ ALGBRC(196)*STATES(54)+ CONSTS(169)*STATES(58)
      ALGBRC(61) =  3.00000*ALGBRC(3)
      ALGBRC(176) =  CONSTS(49)*ALGBRC(167)
      ALGBRC(191) =  (ALGBRC(122)+ALGBRC(61)+ALGBRC(176))*STATES(15)
      ALGBRC(126) =  2.00000*ALGBRC(14)
      ALGBRC(201) =  ALGBRC(52)*STATES(14)+ ALGBRC(126)*STATES(16)+ CONSTS(148)*STATES(20)
      RATES(15) = ALGBRC(201) - ALGBRC(191)
      ALGBRC(100) =  3.00000*ALGBRC(25)
      ALGBRC(189) =  (ALGBRC(138)+ALGBRC(100)+CONSTS(148))*STATES(20)
      ALGBRC(143) =  2.00000*ALGBRC(36)
      ALGBRC(199) =  ALGBRC(90)*STATES(19)+ ALGBRC(143)*STATES(21)+ ALGBRC(176)*STATES(15)
      RATES(20) = ALGBRC(199) - ALGBRC(189)
      ALGBRC(124) = ALGBRC(8)
      ALGBRC(203) =  4.00000*ALGBRC(19)
      RATES(54) = ( ALGBRC(120)*STATES(53) -  (ALGBRC(196)+ALGBRC(124)+CONSTS(172))*STATES(54))+ ALGBRC(203)*STATES(55)+ CONSTS(171)*STATES(59)
      ALGBRC(70) =  2.00000*ALGBRC(3)
      ALGBRC(187) =  CONSTS(49)*ALGBRC(176)
      ALGBRC(202) =  (ALGBRC(126)+ALGBRC(70)+ALGBRC(187))*STATES(16)
      ALGBRC(130) =  3.00000*ALGBRC(14)
      ALGBRC(208) =  ALGBRC(61)*STATES(15)+ ALGBRC(130)*STATES(17)+ CONSTS(151)*STATES(21)
      RATES(16) = ALGBRC(208) - ALGBRC(202)
      ALGBRC(110) =  2.00000*ALGBRC(25)
      ALGBRC(200) =  (ALGBRC(143)+ALGBRC(110)+CONSTS(151))*STATES(21)
      ALGBRC(149) =  3.00000*ALGBRC(36)
      ALGBRC(206) =  ALGBRC(100)*STATES(20)+ ALGBRC(149)*STATES(22)+ ALGBRC(187)*STATES(16)
      RATES(21) = ALGBRC(206) - ALGBRC(200)
      ALGBRC(80) = ALGBRC(3)
      ALGBRC(198) =  CONSTS(49)*ALGBRC(187)
      ALGBRC(209) =  (ALGBRC(130)+ALGBRC(80)+ALGBRC(198))*STATES(17)
      ALGBRC(134) =  4.00000*ALGBRC(14)
      ALGBRC(214) =  ALGBRC(70)*STATES(16)+ ALGBRC(134)*STATES(18)+ CONSTS(157)*STATES(22)
      RATES(17) = ALGBRC(214) - ALGBRC(209)
      ALGBRC(118) = ALGBRC(25)
      ALGBRC(207) =  (ALGBRC(149)+ALGBRC(118)+CONSTS(157))*STATES(22)
      ALGBRC(157) =  4.00000*ALGBRC(36)
      ALGBRC(212) =  ALGBRC(110)*STATES(21)+ ALGBRC(157)*STATES(23)+ ALGBRC(198)*STATES(17)
      RATES(22) = ALGBRC(212) - ALGBRC(207)
      ALGBRC(29) =  CONSTS(150)*CONSTS(126)*EXP(- 196337./CONSTS(146)+529.952/CONSTS(3)+( 2.78085*STATES(7))/CONSTS(127))
      ALGBRC(128) = ALGBRC(29)
      ALGBRC(88) =  CONSTS(150)*CONSTS(126)*EXP(- 147814./CONSTS(146)+338.915/CONSTS(3)+( 2.13600*STATES(7))/CONSTS(127))
      ALGBRC(165) = ALGBRC(88)
      ALGBRC(40) =  CONSTS(150)*CONSTS(126)*EXP(- 133690./CONSTS(146)+229.205/CONSTS(3)+( - 1.55804*STATES(7))/CONSTS(127))
      ALGBRC(210) = ALGBRC(40)
      ALGBRC(98) =  CONSTS(150)*CONSTS(126)*EXP(- 121322./CONSTS(146)+193.265/CONSTS(3)+( - 1.74290*STATES(7))/CONSTS(127))
      ALGBRC(217) = ALGBRC(98)
      RATES(55) = ( ALGBRC(124)*STATES(54) -  (ALGBRC(203)+ALGBRC(128)+ALGBRC(165)+CONSTS(174))*STATES(55))+ ALGBRC(210)*STATES(46)+ ALGBRC(217)*STATES(47)+ CONSTS(173)*STATES(60)
      RATES(47) = ( ALGBRC(165)*STATES(55)+ CONSTS(159)*STATES(46)) -  (ALGBRC(217)+CONSTS(160))*STATES(47)
      ALGBRC(132) =  ALGBRC(108)*CONSTS(163)
      ALGBRC(216) = ALGBRC(174)/CONSTS(163)
      RATES(56) = ( CONSTS(166)*STATES(51) -  (CONSTS(165)+ALGBRC(132))*STATES(56))+ ALGBRC(216)*STATES(57)
      ALGBRC(205) =  CONSTS(49)*ALGBRC(198)
      ALGBRC(215) =  (ALGBRC(134)+CONSTS(46)+ALGBRC(205))*STATES(18)
      ALGBRC(220) =  ALGBRC(80)*STATES(17)+ CONSTS(47)*STATES(24)+ CONSTS(162)*STATES(23)
      RATES(18) = ALGBRC(220) - ALGBRC(215)
      ALGBRC(213) =  (ALGBRC(157)+CONSTS(162))*STATES(23)
      ALGBRC(219) =  ALGBRC(118)*STATES(22)+ ALGBRC(205)*STATES(18)
      RATES(23) = ALGBRC(219) - ALGBRC(213)
      ALGBRC(136) =  ALGBRC(116)*CONSTS(163)
      ALGBRC(221) = ALGBRC(185)/CONSTS(163)
      RATES(57) = (( CONSTS(168)*STATES(52)+ ALGBRC(132)*STATES(56)) -  (ALGBRC(216)+CONSTS(167)+ALGBRC(136))*STATES(57))+ ALGBRC(221)*STATES(58)
      ALGBRC(141) =  ALGBRC(120)*CONSTS(163)
      ALGBRC(223) = ALGBRC(196)/CONSTS(163)
      RATES(58) = (( CONSTS(170)*STATES(53)+ ALGBRC(136)*STATES(57)) -  (ALGBRC(141)+CONSTS(169)+ALGBRC(221))*STATES(58))+ ALGBRC(223)*STATES(59)
      ALGBRC(147) =  ALGBRC(124)*CONSTS(163)
      ALGBRC(225) = ALGBRC(203)/CONSTS(163)
      RATES(59) = (( CONSTS(172)*STATES(54)+ ALGBRC(141)*STATES(58)) -  (ALGBRC(223)+CONSTS(171)+ALGBRC(147))*STATES(59))+ ALGBRC(225)*STATES(60)
      ALGBRC(65) =  CONSTS(150)*CONSTS(126)*EXP(116431./CONSTS(146)+- 578.317/CONSTS(3)+( 0.764126*STATES(7))/CONSTS(127))
      ALGBRC(155) = ALGBRC(65)
      ALGBRC(76) =  CONSTS(150)*CONSTS(126)*EXP(- 55700.7/CONSTS(146)+- 130.639/CONSTS(3)+( - 3.64982*STATES(7))/CONSTS(127))
      ALGBRC(227) = ALGBRC(76)
      RATES(60) = (( CONSTS(174)*STATES(55)+ ALGBRC(147)*STATES(59)) -  (ALGBRC(225)+CONSTS(173)+ALGBRC(155))*STATES(60))+ ALGBRC(227)*STATES(61)
      ALGBRC(56) =  CONSTS(150)*CONSTS(126)*EXP(- 97657.8/CONSTS(146)+1.51000/CONSTS(3)+( 0.0684862*STATES(7))/CONSTS(127))
      ALGBRC(229) = ALGBRC(56)
      ALGBRC(48) =  CONSTS(150)*CONSTS(126)*EXP(- 62123.1/CONSTS(146)+39.2950/CONSTS(3)+( 0.288816*STATES(7))/CONSTS(127))
      ALGBRC(231) = ALGBRC(48)
      RATES(46) = ( ALGBRC(128)*STATES(55) -  (ALGBRC(210)+CONSTS(159)+ALGBRC(231))*STATES(46))+ CONSTS(160)*STATES(47)+ ALGBRC(229)*STATES(61)
      RATES(61) = ( ALGBRC(231)*STATES(46)+ ALGBRC(155)*STATES(60)) -  (ALGBRC(227)+ALGBRC(229))*STATES(61)
      ALGBRC(117) =  CONSTS(116)*log(CONSTS(10)/STATES(1))
      ALGBRC(129) =  CONSTS(83)*(STATES(46)+STATES(47))*(STATES(7) - ALGBRC(117))
      ALGBRC(226) =  CONSTS(84)*(STATES(7) - ALGBRC(117))
      ALGBRC(232) = STATES(7)/CONSTS(116)
      ALGBRC(238) =  EXP( CONSTS(91)*ALGBRC(232))*STATES(1) ** 3.00000*CONSTS(11)
      ALGBRC(239) =  EXP( (CONSTS(91) - 1.00000)*ALGBRC(232))*CONSTS(10) ** 3.00000*STATES(3)
      ALGBRC(240) = 1.00000+ CONSTS(90)*EXP( (CONSTS(91) - 1.00000)*ALGBRC(232))
      ALGBRC(241) = ( CONSTS(87)*(ALGBRC(238) - ALGBRC(239)))/( CONSTS(161)*ALGBRC(240)*CONSTS(164))
      ALGBRC(233) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(232))
      ALGBRC(234) =  0.0365000*CONSTS(152)*EXP( - 1.33000*ALGBRC(232))
      ALGBRC(235) = 1.00000/(ALGBRC(233)+ALGBRC(234))
      ALGBRC(236) = 1.00000+CONSTS(93)/STATES(1) ** 1.50000
      ALGBRC(237) = ( CONSTS(92)*ALGBRC(235)*CONSTS(156))/ALGBRC(236)
      ALGBRC(148) = STATES(7)/CONSTS(116)
      ALGBRC(156) =  1000.00*CONSTS(1)*ALGBRC(148)
      ALGBRC(175) = EXP(ALGBRC(148)) - 1.00000
      ALGBRC(197) =  STATES(1)*EXP(ALGBRC(148)) - CONSTS(10)
      ALGBRC(204) = ( 0.0200000*CONSTS(122)*STATES(40)*ALGBRC(156)*ALGBRC(197))/ALGBRC(175)
      RATES(1) =  - (ALGBRC(129)+ALGBRC(226)+ 3.00000*(ALGBRC(241)+ALGBRC(237))+ALGBRC(204))*CONSTS(144)
      ALGBRC(125) =  0.500000*CONSTS(116)*log(CONSTS(11)/STATES(3))
      ALGBRC(228) =  CONSTS(97)*(STATES(7) - ALGBRC(125))
      ALGBRC(230) = ( CONSTS(95)*STATES(3))/(CONSTS(96)+STATES(3))
      ALGBRC(59) = (STATES(4) - STATES(3))/CONSTS(25)
      ALGBRC(89) =  CONSTS(26)*ALGBRC(68)+ CONSTS(27)*ALGBRC(79)
      ALGBRC(109) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(35))/STATES(3)+CONSTS(35) ** 2.00000+( CONSTS(34)*CONSTS(37))/STATES(3)+CONSTS(37) ** 2.00000)
      RATES(3) =  ALGBRC(109)*(((ALGBRC(59) - ALGBRC(43)) - ALGBRC(89)) -  ((ALGBRC(228) -  2.00000*ALGBRC(241))+ALGBRC(230))*0.500000*CONSTS(144))
      ALGBRC(242) = STATES(7)/CONSTS(116)
      ALGBRC(243) =  1000.00*CONSTS(1)*ALGBRC(242)
      ALGBRC(244) =  0.00100000*EXP( 2.00000*ALGBRC(242)) -  CONSTS(11)*0.341000
      ALGBRC(245) = EXP( 2.00000*ALGBRC(242)) - 1.00000
      ALGBRC(246) = ( CONSTS(145)*4.00000*ALGBRC(243)*ALGBRC(244))/ALGBRC(245)
      ALGBRC(247) =  ALGBRC(246)*STATES(25)*STATES(24)
      ALGBRC(99) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(35))/STATES(4)+CONSTS(35) ** 2.00000+( CONSTS(34)*CONSTS(37))/STATES(4)+CONSTS(37) ** 2.00000)
      RATES(4) =  ALGBRC(99)*((( ALGBRC(51)*CONSTS(6))/CONSTS(8) - ( ALGBRC(59)*CONSTS(5))/CONSTS(8)) -  ALGBRC(247)*CONSTS(149))
      ALGBRC(121) =  CONSTS(116)*log(CONSTS(9)/STATES(2))
      ALGBRC(133) =  CONSTS(80)*CONSTS(123)*STATES(48)*(STATES(7) - ALGBRC(121))
      ALGBRC(137) =  CONSTS(81)*(STATES(49)+STATES(50))*(STATES(7) - ALGBRC(121))
      ALGBRC(222) = 1.00000/(0.940000+EXP( (1.26000/CONSTS(116))*(STATES(7) - ALGBRC(121))))
      ALGBRC(224) =  CONSTS(82)* (CONSTS(9)/1.00000) ** (1.0 / 2)*ALGBRC(222)*(STATES(7) - ALGBRC(121))
      ALGBRC(248) = TERNRY(ALGBRC(246).GE.0.00000, 0.00000, ALGBRC(246))
      ALGBRC(249) = CONSTS(125)/(1.00000+ALGBRC(248)/CONSTS(98))
      ALGBRC(250) =  STATES(2)*EXP(ALGBRC(242)) - CONSTS(9)
      ALGBRC(251) = EXP(ALGBRC(242)) - 1.00000
      ALGBRC(252) = ( ALGBRC(249)*STATES(24)*STATES(25)*ALGBRC(243)*ALGBRC(250))/ALGBRC(251)
      ALGBRC(1) =  INT(VOI/CONSTS(12))*CONSTS(12)
      ALGBRC(12) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(15).AND.VOI - ALGBRC(1).LE.CONSTS(15)+CONSTS(13), CONSTS(14), 0.00000)
      ALGBRC(142) =  CONSTS(121)*STATES(30)*(STATES(7) - ALGBRC(121))
      ALGBRC(166) =  STATES(2)*EXP(ALGBRC(148)) - CONSTS(9)
      ALGBRC(186) = ( CONSTS(122)*STATES(40)*ALGBRC(156)*ALGBRC(166))/ALGBRC(175)
      RATES(2) =  - (((ALGBRC(133)+ALGBRC(137)+ALGBRC(224)+ALGBRC(252)+ALGBRC(12)) -  2.00000*ALGBRC(237))+ALGBRC(142)+ALGBRC(186))*CONSTS(144)
      ALGBRC(211) = ALGBRC(186)+ALGBRC(204)
      ALGBRC(218) = ALGBRC(142)+ALGBRC(211)
      ALGBRC(253) = ALGBRC(129)+ALGBRC(247)+ALGBRC(252)+ALGBRC(133)+ALGBRC(137)+ALGBRC(224)+ALGBRC(241)+ALGBRC(237)+ALGBRC(218)+ALGBRC(230)+ALGBRC(228)+ALGBRC(226)+ALGBRC(12)
      RATES(7) = - ALGBRC(253)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) =  - CONSTS(38)* STATES(4)*1000.00 ** CONSTS(44)*STATES(12)+ CONSTS(39)*STATES(8)
      ALGBRC(13) =  CONSTS(40)* STATES(4)*1000.00 ** CONSTS(45)*STATES(8) -  CONSTS(41)*STATES(9)
      ALGBRC(4) = (CONSTS(118)/(1.00000+EXP((STATES(7)+28.5000)/7.80000))+1.00000) - CONSTS(118)
      ALGBRC(15) = 1.00000/(0.00653000/(0.500000+EXP(STATES(7)/- 7.10000))+ 0.00512000*EXP(STATES(7)/- 39.8000))
      ALGBRC(9) =  CONSTS(101)*CONSTS(102)*EXP( CONSTS(103)*STATES(7))
      ALGBRC(20) =  CONSTS(101)*CONSTS(104)*EXP( CONSTS(105)*STATES(7))
      ALGBRC(24) =  CONSTS(42)*STATES(8) -  CONSTS(43)*STATES(13)
      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(35) = - (ALGBRC(2)+ALGBRC(13)+ALGBRC(24))
      ALGBRC(31) =  ALGBRC(9)*STATES(62)+ CONSTS(129)*STATES(64)
      ALGBRC(42) =  (ALGBRC(20)+CONSTS(128))*STATES(63)
      ALGBRC(10) =  CONSTS(101)*CONSTS(106)*EXP( CONSTS(107)*STATES(7))
      ALGBRC(41) =  CONSTS(101)*CONSTS(112)*EXP( CONSTS(113)*STATES(7))
      ALGBRC(50) =  ALGBRC(10)*STATES(64)+ ALGBRC(41)*STATES(65)
      ALGBRC(21) =  CONSTS(101)*CONSTS(108)*EXP( CONSTS(109)*STATES(7))
      ALGBRC(30) =  CONSTS(101)*CONSTS(110)*EXP( CONSTS(111)*STATES(7))
      ALGBRC(58) =  (ALGBRC(21)+ALGBRC(30))*STATES(48)
      ALGBRC(23) = STATES(3)/CONSTS(16) ** CONSTS(17)
      ALGBRC(34) = STATES(6)/CONSTS(18) ** CONSTS(19)
      ALGBRC(43) = ( CONSTS(20)*( CONSTS(21)*ALGBRC(23) -  CONSTS(22)*ALGBRC(34)))/(1.00000+ALGBRC(23)+ALGBRC(34))
      ALGBRC(60) = (STATES(6) - STATES(5))/CONSTS(24)
      ALGBRC(51) =  CONSTS(23)*(STATES(8)+STATES(9))*(STATES(5) - STATES(4))
      ALGBRC(69) = 1.00000/(1.00000+( CONSTS(33)*CONSTS(36))/STATES(5)+CONSTS(36) ** 2.00000)
      ALGBRC(68) =  CONSTS(30)*STATES(3)*(1.00000 - STATES(10)) -  CONSTS(31)*STATES(10)
      ALGBRC(49) =  CONSTS(101)*CONSTS(114)*EXP( CONSTS(115)*STATES(7))
      ALGBRC(57) = ( ALGBRC(21)*ALGBRC(41)*ALGBRC(49))/( ALGBRC(10)*ALGBRC(30))
      ALGBRC(66) =  CONSTS(128)*STATES(63)+ ALGBRC(21)*STATES(48)+ ALGBRC(57)*STATES(65)
      ALGBRC(77) =  (ALGBRC(49)+ALGBRC(10)+CONSTS(129))*STATES(64)
      ALGBRC(67) =  ALGBRC(49)*STATES(64)+ ALGBRC(30)*STATES(48)
      ALGBRC(78) =  (ALGBRC(57)+ALGBRC(41))*STATES(65)
      ALGBRC(79) =  CONSTS(28)*STATES(3)*(1.00000 - STATES(11)) -  CONSTS(29)*STATES(11)
      ALGBRC(6) =  CONSTS(74)*EXP( CONSTS(75)*STATES(7))
      ALGBRC(27) =  4.00000*ALGBRC(6)
      ALGBRC(74) =  (ALGBRC(27)+CONSTS(133))*STATES(36)
      ALGBRC(17) =  CONSTS(76)*EXP( - CONSTS(77)*STATES(7))
      ALGBRC(63) = ALGBRC(17)
      ALGBRC(84) =  ALGBRC(63)*STATES(37)+ CONSTS(138)*STATES(41)
      ALGBRC(28) =  4.00000*CONSTS(70)*ALGBRC(6)
      ALGBRC(75) =  (CONSTS(138)+ALGBRC(28))*STATES(41)
      ALGBRC(64) = ALGBRC(17)/CONSTS(66)
      ALGBRC(86) =  CONSTS(133)*STATES(36)+ ALGBRC(64)*STATES(42)
      ALGBRC(38) =  3.00000*ALGBRC(6)
      ALGBRC(85) =  (ALGBRC(38)+ALGBRC(63)+CONSTS(134))*STATES(37)
      ALGBRC(72) =  2.00000*ALGBRC(17)
      ALGBRC(94) =  ALGBRC(72)*STATES(38)+ CONSTS(139)*STATES(42)+ ALGBRC(27)*STATES(36)
      ALGBRC(39) = ( 3.00000*CONSTS(71)*ALGBRC(6))/CONSTS(70)
      ALGBRC(87) =  (ALGBRC(39)+CONSTS(139)+ALGBRC(64))*STATES(42)
      ALGBRC(73) = ( 2.00000*CONSTS(66)*ALGBRC(17))/CONSTS(67)
      ALGBRC(96) =  ALGBRC(73)*STATES(43)+ CONSTS(134)*STATES(37)+ ALGBRC(28)*STATES(41)
      ALGBRC(46) =  2.00000*ALGBRC(6)
      ALGBRC(95) =  (ALGBRC(46)+ALGBRC(72)+CONSTS(135))*STATES(38)
      ALGBRC(82) =  3.00000*ALGBRC(17)
      ALGBRC(102) =  ALGBRC(82)*STATES(39)+ CONSTS(140)*STATES(43)+ ALGBRC(38)*STATES(37)
      ALGBRC(47) = ( 2.00000*CONSTS(72)*ALGBRC(6))/CONSTS(71)
      ALGBRC(97) =  (ALGBRC(47)+CONSTS(140)+ALGBRC(73))*STATES(43)
      ALGBRC(83) = ( 3.00000*CONSTS(67)*ALGBRC(17))/CONSTS(68)
      ALGBRC(105) =  ALGBRC(83)*STATES(44)+ CONSTS(135)*STATES(38)+ ALGBRC(39)*STATES(42)
      ALGBRC(54) = ALGBRC(6)
      ALGBRC(103) =  (ALGBRC(54)+ALGBRC(82)+CONSTS(136))*STATES(39)
      ALGBRC(92) =  4.00000*ALGBRC(17)
      ALGBRC(112) =  ALGBRC(92)*STATES(40)+ CONSTS(141)*STATES(44)+ ALGBRC(46)*STATES(38)
      ALGBRC(104) =  (ALGBRC(92)+CONSTS(137))*STATES(40)
      ALGBRC(113) =  ALGBRC(54)*STATES(39)+ CONSTS(142)*STATES(45)
      ALGBRC(55) = ( CONSTS(73)*ALGBRC(6))/CONSTS(72)
      ALGBRC(106) =  (ALGBRC(55)+CONSTS(141)+ALGBRC(83))*STATES(44)
      ALGBRC(93) = ( 4.00000*CONSTS(68)*ALGBRC(17))/CONSTS(69)
      ALGBRC(114) =  ALGBRC(93)*STATES(45)+ CONSTS(136)*STATES(39)+ ALGBRC(47)*STATES(43)
      ALGBRC(107) =  (CONSTS(142)+ALGBRC(93))*STATES(45)
      ALGBRC(115) =  CONSTS(137)*STATES(40)+ ALGBRC(55)*STATES(44)
      ALGBRC(5) =  CONSTS(50)*EXP( CONSTS(51)*STATES(7))
      ALGBRC(45) =  4.00000*ALGBRC(5)
      ALGBRC(37) =  CONSTS(56)*EXP( CONSTS(57)*STATES(7))
      ALGBRC(139) = ALGBRC(37)
      ALGBRC(153) =  (ALGBRC(45)+ALGBRC(139))*STATES(26)
      ALGBRC(16) =  CONSTS(52)*EXP( - CONSTS(53)*STATES(7))
      ALGBRC(119) = ALGBRC(16)
      ALGBRC(26) =  CONSTS(54)*EXP( - CONSTS(55)*STATES(7))
      ALGBRC(146) = ALGBRC(26)
      ALGBRC(161) =  ALGBRC(119)*STATES(27)+ ALGBRC(146)*STATES(31)
      ALGBRC(81) =  4.00000*CONSTS(62)*ALGBRC(5)
      ALGBRC(154) =  (ALGBRC(146)+ALGBRC(81))*STATES(31)
      ALGBRC(135) = ALGBRC(16)/CONSTS(58)
      ALGBRC(163) =  ALGBRC(139)*STATES(26)+ ALGBRC(135)*STATES(32)
      ALGBRC(53) =  3.00000*ALGBRC(5)
      ALGBRC(144) =  CONSTS(58)*ALGBRC(37)
      ALGBRC(162) =  (ALGBRC(53)+ALGBRC(119)+ALGBRC(144))*STATES(27)
      ALGBRC(123) =  2.00000*ALGBRC(16)
      ALGBRC(152) = ALGBRC(26)/CONSTS(62)
      ALGBRC(170) =  ALGBRC(123)*STATES(28)+ ALGBRC(152)*STATES(32)+ ALGBRC(45)*STATES(26)
      ALGBRC(91) = ( 3.00000*CONSTS(63)*ALGBRC(5))/CONSTS(62)
      ALGBRC(164) =  (ALGBRC(91)+ALGBRC(152)+ALGBRC(135))*STATES(32)
      ALGBRC(140) = ( 2.00000*CONSTS(58)*ALGBRC(16))/CONSTS(59)
      ALGBRC(172) =  ALGBRC(140)*STATES(33)+ ALGBRC(144)*STATES(27)+ ALGBRC(81)*STATES(31)
      ALGBRC(8) =  CONSTS(150)*CONSTS(126)*EXP(- 114007./CONSTS(146)+224.114/CONSTS(3)+( 0.286374*STATES(7))/CONSTS(127))
      ALGBRC(108) =  4.00000*ALGBRC(8)
      ALGBRC(19) =  CONSTS(150)*CONSTS(126)*EXP(- 272470./CONSTS(146)+708.146/CONSTS(3)+( - 2.28528*STATES(7))/CONSTS(127))
      ALGBRC(174) = ALGBRC(19)
      ALGBRC(62) =  2.00000*ALGBRC(5)
      ALGBRC(150) =  CONSTS(59)*ALGBRC(37)
      ALGBRC(171) =  (ALGBRC(62)+ALGBRC(123)+ALGBRC(150))*STATES(28)
      ALGBRC(127) =  3.00000*ALGBRC(16)
      ALGBRC(160) = ALGBRC(26)/CONSTS(63)
      ALGBRC(179) =  ALGBRC(127)*STATES(29)+ ALGBRC(160)*STATES(33)+ ALGBRC(53)*STATES(27)
      ALGBRC(101) = ( 2.00000*CONSTS(64)*ALGBRC(5))/CONSTS(63)
      ALGBRC(173) =  (ALGBRC(101)+ALGBRC(160)+ALGBRC(140))*STATES(33)
      ALGBRC(145) = ( 3.00000*CONSTS(59)*ALGBRC(16))/CONSTS(60)
      ALGBRC(182) =  ALGBRC(145)*STATES(34)+ ALGBRC(150)*STATES(28)+ ALGBRC(91)*STATES(32)
      ALGBRC(116) =  3.00000*ALGBRC(8)
      ALGBRC(185) =  2.00000*ALGBRC(19)
      ALGBRC(3) =  4.00000*1.20000*0.416000*EXP( 0.0120000*(STATES(7) - 35.0000))
      ALGBRC(52) =  4.00000*ALGBRC(3)
      ALGBRC(44) =  0.600000*0.0923300*STATES(4)
      ALGBRC(167) = ALGBRC(44)
      ALGBRC(178) =  (ALGBRC(52)+ALGBRC(167))*STATES(14)
      ALGBRC(14) =  4.00000*0.450000*0.0490000*EXP( - 0.0650000*(STATES(7) - 22.0000))
      ALGBRC(122) = ALGBRC(14)
      ALGBRC(190) =  ALGBRC(122)*STATES(15)+ CONSTS(132)*STATES(19)
      ALGBRC(25) =  CONSTS(49)*ALGBRC(3)
      ALGBRC(90) =  4.00000*ALGBRC(25)
      ALGBRC(177) =  (ALGBRC(90)+CONSTS(132))*STATES(19)
      ALGBRC(36) = ALGBRC(14)/CONSTS(48)
      ALGBRC(138) = ALGBRC(36)
      ALGBRC(188) =  ALGBRC(138)*STATES(20)+ ALGBRC(167)*STATES(14)
      ALGBRC(71) = ALGBRC(5)
      ALGBRC(158) =  CONSTS(60)*ALGBRC(37)
      ALGBRC(180) =  (ALGBRC(71)+ALGBRC(127)+ALGBRC(158))*STATES(29)
      ALGBRC(131) =  4.00000*ALGBRC(16)
      ALGBRC(168) = ALGBRC(26)/CONSTS(64)
      ALGBRC(192) =  ALGBRC(131)*STATES(30)+ ALGBRC(168)*STATES(34)+ ALGBRC(62)*STATES(28)
      ALGBRC(159) =  CONSTS(61)*ALGBRC(37)
      ALGBRC(181) =  (ALGBRC(131)+ALGBRC(159))*STATES(30)
      ALGBRC(169) = ALGBRC(26)/CONSTS(65)
      ALGBRC(193) =  ALGBRC(71)*STATES(29)+ ALGBRC(169)*STATES(35)
      ALGBRC(111) = ( CONSTS(65)*ALGBRC(5))/CONSTS(64)
      ALGBRC(183) =  (ALGBRC(111)+ALGBRC(168)+ALGBRC(145))*STATES(34)
      ALGBRC(151) = ( 4.00000*CONSTS(60)*ALGBRC(16))/CONSTS(61)
      ALGBRC(194) =  ALGBRC(151)*STATES(35)+ ALGBRC(158)*STATES(29)+ ALGBRC(101)*STATES(33)
      ALGBRC(184) =  (ALGBRC(169)+ALGBRC(151))*STATES(35)
      ALGBRC(195) =  ALGBRC(159)*STATES(30)+ ALGBRC(111)*STATES(34)
      ALGBRC(120) =  2.00000*ALGBRC(8)
      ALGBRC(196) =  3.00000*ALGBRC(19)
      ALGBRC(61) =  3.00000*ALGBRC(3)
      ALGBRC(176) =  CONSTS(49)*ALGBRC(167)
      ALGBRC(191) =  (ALGBRC(122)+ALGBRC(61)+ALGBRC(176))*STATES(15)
      ALGBRC(126) =  2.00000*ALGBRC(14)
      ALGBRC(201) =  ALGBRC(52)*STATES(14)+ ALGBRC(126)*STATES(16)+ CONSTS(148)*STATES(20)
      ALGBRC(100) =  3.00000*ALGBRC(25)
      ALGBRC(189) =  (ALGBRC(138)+ALGBRC(100)+CONSTS(148))*STATES(20)
      ALGBRC(143) =  2.00000*ALGBRC(36)
      ALGBRC(199) =  ALGBRC(90)*STATES(19)+ ALGBRC(143)*STATES(21)+ ALGBRC(176)*STATES(15)
      ALGBRC(124) = ALGBRC(8)
      ALGBRC(203) =  4.00000*ALGBRC(19)
      ALGBRC(70) =  2.00000*ALGBRC(3)
      ALGBRC(187) =  CONSTS(49)*ALGBRC(176)
      ALGBRC(202) =  (ALGBRC(126)+ALGBRC(70)+ALGBRC(187))*STATES(16)
      ALGBRC(130) =  3.00000*ALGBRC(14)
      ALGBRC(208) =  ALGBRC(61)*STATES(15)+ ALGBRC(130)*STATES(17)+ CONSTS(151)*STATES(21)
      ALGBRC(110) =  2.00000*ALGBRC(25)
      ALGBRC(200) =  (ALGBRC(143)+ALGBRC(110)+CONSTS(151))*STATES(21)
      ALGBRC(149) =  3.00000*ALGBRC(36)
      ALGBRC(206) =  ALGBRC(100)*STATES(20)+ ALGBRC(149)*STATES(22)+ ALGBRC(187)*STATES(16)
      ALGBRC(80) = ALGBRC(3)
      ALGBRC(198) =  CONSTS(49)*ALGBRC(187)
      ALGBRC(209) =  (ALGBRC(130)+ALGBRC(80)+ALGBRC(198))*STATES(17)
      ALGBRC(134) =  4.00000*ALGBRC(14)
      ALGBRC(214) =  ALGBRC(70)*STATES(16)+ ALGBRC(134)*STATES(18)+ CONSTS(157)*STATES(22)
      ALGBRC(118) = ALGBRC(25)
      ALGBRC(207) =  (ALGBRC(149)+ALGBRC(118)+CONSTS(157))*STATES(22)
      ALGBRC(157) =  4.00000*ALGBRC(36)
      ALGBRC(212) =  ALGBRC(110)*STATES(21)+ ALGBRC(157)*STATES(23)+ ALGBRC(198)*STATES(17)
      ALGBRC(29) =  CONSTS(150)*CONSTS(126)*EXP(- 196337./CONSTS(146)+529.952/CONSTS(3)+( 2.78085*STATES(7))/CONSTS(127))
      ALGBRC(128) = ALGBRC(29)
      ALGBRC(88) =  CONSTS(150)*CONSTS(126)*EXP(- 147814./CONSTS(146)+338.915/CONSTS(3)+( 2.13600*STATES(7))/CONSTS(127))
      ALGBRC(165) = ALGBRC(88)
      ALGBRC(40) =  CONSTS(150)*CONSTS(126)*EXP(- 133690./CONSTS(146)+229.205/CONSTS(3)+( - 1.55804*STATES(7))/CONSTS(127))
      ALGBRC(210) = ALGBRC(40)
      ALGBRC(98) =  CONSTS(150)*CONSTS(126)*EXP(- 121322./CONSTS(146)+193.265/CONSTS(3)+( - 1.74290*STATES(7))/CONSTS(127))
      ALGBRC(217) = ALGBRC(98)
      ALGBRC(132) =  ALGBRC(108)*CONSTS(163)
      ALGBRC(216) = ALGBRC(174)/CONSTS(163)
      ALGBRC(205) =  CONSTS(49)*ALGBRC(198)
      ALGBRC(215) =  (ALGBRC(134)+CONSTS(46)+ALGBRC(205))*STATES(18)
      ALGBRC(220) =  ALGBRC(80)*STATES(17)+ CONSTS(47)*STATES(24)+ CONSTS(162)*STATES(23)
      ALGBRC(213) =  (ALGBRC(157)+CONSTS(162))*STATES(23)
      ALGBRC(219) =  ALGBRC(118)*STATES(22)+ ALGBRC(205)*STATES(18)
      ALGBRC(136) =  ALGBRC(116)*CONSTS(163)
      ALGBRC(221) = ALGBRC(185)/CONSTS(163)
      ALGBRC(141) =  ALGBRC(120)*CONSTS(163)
      ALGBRC(223) = ALGBRC(196)/CONSTS(163)
      ALGBRC(147) =  ALGBRC(124)*CONSTS(163)
      ALGBRC(225) = ALGBRC(203)/CONSTS(163)
      ALGBRC(65) =  CONSTS(150)*CONSTS(126)*EXP(116431./CONSTS(146)+- 578.317/CONSTS(3)+( 0.764126*STATES(7))/CONSTS(127))
      ALGBRC(155) = ALGBRC(65)
      ALGBRC(76) =  CONSTS(150)*CONSTS(126)*EXP(- 55700.7/CONSTS(146)+- 130.639/CONSTS(3)+( - 3.64982*STATES(7))/CONSTS(127))
      ALGBRC(227) = ALGBRC(76)
      ALGBRC(56) =  CONSTS(150)*CONSTS(126)*EXP(- 97657.8/CONSTS(146)+1.51000/CONSTS(3)+( 0.0684862*STATES(7))/CONSTS(127))
      ALGBRC(229) = ALGBRC(56)
      ALGBRC(48) =  CONSTS(150)*CONSTS(126)*EXP(- 62123.1/CONSTS(146)+39.2950/CONSTS(3)+( 0.288816*STATES(7))/CONSTS(127))
      ALGBRC(231) = ALGBRC(48)
      ALGBRC(117) =  CONSTS(116)*log(CONSTS(10)/STATES(1))
      ALGBRC(129) =  CONSTS(83)*(STATES(46)+STATES(47))*(STATES(7) - ALGBRC(117))
      ALGBRC(226) =  CONSTS(84)*(STATES(7) - ALGBRC(117))
      ALGBRC(232) = STATES(7)/CONSTS(116)
      ALGBRC(238) =  EXP( CONSTS(91)*ALGBRC(232))*STATES(1) ** 3.00000*CONSTS(11)
      ALGBRC(239) =  EXP( (CONSTS(91) - 1.00000)*ALGBRC(232))*CONSTS(10) ** 3.00000*STATES(3)
      ALGBRC(240) = 1.00000+ CONSTS(90)*EXP( (CONSTS(91) - 1.00000)*ALGBRC(232))
      ALGBRC(241) = ( CONSTS(87)*(ALGBRC(238) - ALGBRC(239)))/( CONSTS(161)*ALGBRC(240)*CONSTS(164))
      ALGBRC(233) = 1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(232))
      ALGBRC(234) =  0.0365000*CONSTS(152)*EXP( - 1.33000*ALGBRC(232))
      ALGBRC(235) = 1.00000/(ALGBRC(233)+ALGBRC(234))
      ALGBRC(236) = 1.00000+CONSTS(93)/STATES(1) ** 1.50000
      ALGBRC(237) = ( CONSTS(92)*ALGBRC(235)*CONSTS(156))/ALGBRC(236)
      ALGBRC(148) = STATES(7)/CONSTS(116)
      ALGBRC(156) =  1000.00*CONSTS(1)*ALGBRC(148)
      ALGBRC(175) = EXP(ALGBRC(148)) - 1.00000
      ALGBRC(197) =  STATES(1)*EXP(ALGBRC(148)) - CONSTS(10)
      ALGBRC(204) = ( 0.0200000*CONSTS(122)*STATES(40)*ALGBRC(156)*ALGBRC(197))/ALGBRC(175)
      ALGBRC(125) =  0.500000*CONSTS(116)*log(CONSTS(11)/STATES(3))
      ALGBRC(228) =  CONSTS(97)*(STATES(7) - ALGBRC(125))
      ALGBRC(230) = ( CONSTS(95)*STATES(3))/(CONSTS(96)+STATES(3))
      ALGBRC(59) = (STATES(4) - STATES(3))/CONSTS(25)
      ALGBRC(89) =  CONSTS(26)*ALGBRC(68)+ CONSTS(27)*ALGBRC(79)
      ALGBRC(109) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(35))/STATES(3)+CONSTS(35) ** 2.00000+( CONSTS(34)*CONSTS(37))/STATES(3)+CONSTS(37) ** 2.00000)
      ALGBRC(242) = STATES(7)/CONSTS(116)
      ALGBRC(243) =  1000.00*CONSTS(1)*ALGBRC(242)
      ALGBRC(244) =  0.00100000*EXP( 2.00000*ALGBRC(242)) -  CONSTS(11)*0.341000
      ALGBRC(245) = EXP( 2.00000*ALGBRC(242)) - 1.00000
      ALGBRC(246) = ( CONSTS(145)*4.00000*ALGBRC(243)*ALGBRC(244))/ALGBRC(245)
      ALGBRC(247) =  ALGBRC(246)*STATES(25)*STATES(24)
      ALGBRC(99) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(35))/STATES(4)+CONSTS(35) ** 2.00000+( CONSTS(34)*CONSTS(37))/STATES(4)+CONSTS(37) ** 2.00000)
      ALGBRC(121) =  CONSTS(116)*log(CONSTS(9)/STATES(2))
      ALGBRC(133) =  CONSTS(80)*CONSTS(123)*STATES(48)*(STATES(7) - ALGBRC(121))
      ALGBRC(137) =  CONSTS(81)*(STATES(49)+STATES(50))*(STATES(7) - ALGBRC(121))
      ALGBRC(222) = 1.00000/(0.940000+EXP( (1.26000/CONSTS(116))*(STATES(7) - ALGBRC(121))))
      ALGBRC(224) =  CONSTS(82)* (CONSTS(9)/1.00000) ** (1.0 / 2)*ALGBRC(222)*(STATES(7) - ALGBRC(121))
      ALGBRC(248) = TERNRY(ALGBRC(246).GE.0.00000, 0.00000, ALGBRC(246))
      ALGBRC(249) = CONSTS(125)/(1.00000+ALGBRC(248)/CONSTS(98))
      ALGBRC(250) =  STATES(2)*EXP(ALGBRC(242)) - CONSTS(9)
      ALGBRC(251) = EXP(ALGBRC(242)) - 1.00000
      ALGBRC(252) = ( ALGBRC(249)*STATES(24)*STATES(25)*ALGBRC(243)*ALGBRC(250))/ALGBRC(251)
      ALGBRC(1) =  INT(VOI/CONSTS(12))*CONSTS(12)
      ALGBRC(12) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(15).AND.VOI - ALGBRC(1).LE.CONSTS(15)+CONSTS(13), CONSTS(14), 0.00000)
      ALGBRC(142) =  CONSTS(121)*STATES(30)*(STATES(7) - ALGBRC(121))
      ALGBRC(166) =  STATES(2)*EXP(ALGBRC(148)) - CONSTS(9)
      ALGBRC(186) = ( CONSTS(122)*STATES(40)*ALGBRC(156)*ALGBRC(166))/ALGBRC(175)
      ALGBRC(211) = ALGBRC(186)+ALGBRC(204)
      ALGBRC(218) = ALGBRC(142)+ALGBRC(211)
      ALGBRC(253) = ALGBRC(129)+ALGBRC(247)+ALGBRC(252)+ALGBRC(133)+ALGBRC(137)+ALGBRC(224)+ALGBRC(241)+ALGBRC(237)+ALGBRC(218)+ALGBRC(230)+ALGBRC(228)+ALGBRC(226)+ALGBRC(12)
      ALGBRC(7) = STATES(2)+ 0.0183300*STATES(1)
      ALGBRC(18) =  CONSTS(116)*log(CONSTS(143)/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