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 223 entries in the algebraic variable array.
C There are a total of 43 entries in each of the rate and state variable arrays.
C There are a total of 163 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is celltype in component environment (dimensionless).
C CONSTS(2) is nao in component extracellular (millimolar).
C CONSTS(3) is cao in component extracellular (millimolar).
C CONSTS(4) is ko in component extracellular (millimolar).
C CONSTS(5) is clo in component extracellular (millimolar).
C CONSTS(6) is R in component physical_constants (joule_per_kilomole_kelvin).
C CONSTS(7) is T in component physical_constants (kelvin).
C CONSTS(8) is F in component physical_constants (coulomb_per_mole).
C CONSTS(9) is zna in component physical_constants (dimensionless).
C CONSTS(10) is zca in component physical_constants (dimensionless).
C CONSTS(11) is zk in component physical_constants (dimensionless).
C CONSTS(12) is zcl in component physical_constants (dimensionless).
C CONSTS(13) is L in component cell_geometry (centimeter).
C CONSTS(14) is rad in component cell_geometry (centimeter).
C CONSTS(113) is vcell in component cell_geometry (microliter).
C CONSTS(131) is Ageo in component cell_geometry (centimeter_squared).
C CONSTS(137) is Acap in component cell_geometry (centimeter_squared).
C CONSTS(143) is vmyo in component cell_geometry (microliter).
C CONSTS(144) is vnsr in component cell_geometry (microliter).
C CONSTS(145) is vjsr in component cell_geometry (microliter).
C CONSTS(146) is vss in component cell_geometry (microliter).
C STATES(1) is v in component membrane (millivolt).
C ALGBRC(26) is vffrt in component membrane (coulomb_per_mole).
C ALGBRC(29) is vfrt in component membrane (dimensionless).
C ALGBRC(69) is INa in component INa (microA_per_microF).
C ALGBRC(71) is INaL in component INaL (microA_per_microF).
C ALGBRC(77) is Ito in component Ito (microA_per_microF).
C ALGBRC(113) is ICaL in component ICaL (microA_per_microF).
C ALGBRC(114) is ICaNa in component ICaL (microA_per_microF).
C ALGBRC(115) is ICaK in component ICaL (microA_per_microF).
C ALGBRC(116) is IKr in component IKr (microA_per_microF).
C ALGBRC(118) is IKs in component IKs (microA_per_microF).
C ALGBRC(122) is IK1 in component IK1 (microA_per_microF).
C ALGBRC(154) is INaCa_i in component INaCa (microA_per_microF).
C ALGBRC(184) is INaCa_ss in component INaCa (microA_per_microF).
C ALGBRC(203) is INaK in component INaK (microA_per_microF).
C ALGBRC(206) is INab in component INab (microA_per_microF).
C ALGBRC(205) is IKb in component IKb (microA_per_microF).
C ALGBRC(210) is IpCa in component IpCa (microA_per_microF).
C ALGBRC(208) is ICab in component ICab (microA_per_microF).
C ALGBRC(215) is IClCa in component ICl (microA_per_microF).
C ALGBRC(217) is IClb in component ICl (microA_per_microF).
C ALGBRC(67) is I_katp in component I_katp (microA_per_microF).
C ALGBRC(12) is Istim in component membrane (microA_per_microF).
C CONSTS(15) is i_Stim_Start in component membrane (millisecond).
C CONSTS(16) is i_Stim_End in component membrane (millisecond).
C CONSTS(17) is i_Stim_Amplitude in component membrane (microA_per_microF).
C CONSTS(18) is i_Stim_Period in component membrane (millisecond).
C CONSTS(19) is i_Stim_PulseDuration in component membrane (millisecond).
C CONSTS(20) is KmCaMK in component CaMK (millimolar).
C CONSTS(21) is aCaMK in component CaMK (per_millimolar_per_millisecond).
C CONSTS(22) is bCaMK in component CaMK (per_millisecond).
C CONSTS(23) is CaMKo in component CaMK (dimensionless).
C CONSTS(24) is KmCaM in component CaMK (millimolar).
C ALGBRC(44) is CaMKb in component CaMK (millimolar).
C ALGBRC(50) is CaMKa in component CaMK (millimolar).
C STATES(2) is CaMKt in component CaMK (millimolar).
C STATES(3) is cass in component intracellular_ions (millimolar).
C CONSTS(25) is cmdnmax_b in component intracellular_ions (millimolar).
C CONSTS(114) is cmdnmax in component intracellular_ions (millimolar).
C CONSTS(26) is kmcmdn in component intracellular_ions (millimolar).
C CONSTS(27) is trpnmax in component intracellular_ions (millimolar).
C CONSTS(28) is kmtrpn in component intracellular_ions (millimolar).
C CONSTS(29) is BSRmax in component intracellular_ions (millimolar).
C CONSTS(30) is KmBSR in component intracellular_ions (millimolar).
C CONSTS(31) is BSLmax in component intracellular_ions (millimolar).
C CONSTS(32) is KmBSL in component intracellular_ions (millimolar).
C CONSTS(33) is csqnmax in component intracellular_ions (millimolar).
C CONSTS(34) is kmcsqn in component intracellular_ions (millimolar).
C STATES(4) is nai in component intracellular_ions (millimolar).
C STATES(5) is nass in component intracellular_ions (millimolar).
C STATES(6) is ki in component intracellular_ions (millimolar).
C STATES(7) is kss in component intracellular_ions (millimolar).
C STATES(8) is cansr in component intracellular_ions (millimolar).
C STATES(9) is cajsr in component intracellular_ions (millimolar).
C STATES(10) is cai in component intracellular_ions (millimolar).
C CONSTS(35) is cli in component intracellular_ions (millimolar).
C ALGBRC(92) is ICaL_ss in component ICaL (microA_per_microF).
C ALGBRC(93) is ICaNa_ss in component ICaL (microA_per_microF).
C ALGBRC(96) is ICaK_ss in component ICaL (microA_per_microF).
C ALGBRC(110) is ICaL_i in component ICaL (microA_per_microF).
C ALGBRC(111) is ICaNa_i in component ICaL (microA_per_microF).
C ALGBRC(112) is ICaK_i in component ICaL (microA_per_microF).
C ALGBRC(209) is JdiffNa in component diff (millimolar_per_millisecond).
C ALGBRC(212) is Jdiff in component diff (millimolar_per_millisecond).
C ALGBRC(222) is Jup in component SERCA (millimolar_per_millisecond).
C ALGBRC(207) is JdiffK in component diff (millimolar_per_millisecond).
C ALGBRC(216) is Jrel in component ryr (millimolar_per_millisecond).
C ALGBRC(223) is Jtr in component trans_flux (millimolar_per_millisecond).
C ALGBRC(54) is Bcai in component intracellular_ions (dimensionless).
C ALGBRC(60) is Bcajsr in component intracellular_ions (dimensionless).
C ALGBRC(57) is Bcass in component intracellular_ions (dimensionless).
C CONSTS(36) is PKNa in component reversal_potentials (dimensionless).
C ALGBRC(64) is ENa in component reversal_potentials (millivolt).
C ALGBRC(65) is EK in component reversal_potentials (millivolt).
C ALGBRC(66) is EKs in component reversal_potentials (millivolt).
C CONSTS(115) is ECl in component reversal_potentials (millivolt).
C CONSTS(37) is gkatp in component I_katp (milliS_per_microF).
C CONSTS(38) is fkatp in component I_katp (dimensionless).
C CONSTS(39) is K_o_n in component I_katp (millimolar).
C CONSTS(40) is A_atp in component I_katp (millimolar).
C CONSTS(41) is K_atp in component I_katp (millimolar).
C CONSTS(116) is akik in component I_katp (dimensionless).
C CONSTS(117) is bkik in component I_katp (dimensionless).
C ALGBRC(1) is mss in component INa (dimensionless).
C ALGBRC(14) is tm in component INa (millisecond).
C STATES(11) is m in component INa (dimensionless).
C ALGBRC(2) is hss in component INa (dimensionless).
C ALGBRC(15) is ah in component INa (dimensionless).
C ALGBRC(30) is bh in component INa (dimensionless).
C ALGBRC(38) is th in component INa (millisecond).
C STATES(12) is h in component INa (dimensionless).
C ALGBRC(39) is jss in component INa (dimensionless).
C ALGBRC(16) is aj in component INa (dimensionless).
C ALGBRC(31) is bj in component INa (dimensionless).
C ALGBRC(45) is tj in component INa (millisecond).
C STATES(13) is j in component INa (dimensionless).
C ALGBRC(46) is hssp in component INa (dimensionless).
C STATES(14) is hp in component INa (dimensionless).
C ALGBRC(51) is tjp in component INa (millisecond).
C STATES(15) is jp in component INa (dimensionless).
C ALGBRC(68) is fINap in component INa (dimensionless).
C CONSTS(42) is GNa in component INa (milliS_per_microF).
C ALGBRC(3) is mLss in component INaL (dimensionless).
C ALGBRC(17) is tmL in component INaL (millisecond).
C STATES(16) is mL in component INaL (dimensionless).
C CONSTS(43) is thL in component INaL (millisecond).
C ALGBRC(4) is hLss in component INaL (dimensionless).
C STATES(17) is hL in component INaL (dimensionless).
C ALGBRC(5) is hLssp in component INaL (dimensionless).
C CONSTS(118) is thLp in component INaL (millisecond).
C STATES(18) is hLp in component INaL (dimensionless).
C CONSTS(44) is GNaL_b in component INaL (milliS_per_microF).
C CONSTS(119) is GNaL in component INaL (milliS_per_microF).
C ALGBRC(70) is fINaLp in component INaL (dimensionless).
C CONSTS(45) is Gto_b in component Ito (milliS_per_microF).
C ALGBRC(6) is ass in component Ito (dimensionless).
C ALGBRC(18) is ta in component Ito (millisecond).
C STATES(19) is a in component Ito (dimensionless).
C CONSTS(46) is EKshift in component Ito (millivolt).
C ALGBRC(7) is iss in component Ito (dimensionless).
C ALGBRC(19) is delta_epi in component Ito (dimensionless).
C ALGBRC(32) is tiF_b in component Ito (millisecond).
C ALGBRC(40) is tiS_b in component Ito (millisecond).
C ALGBRC(47) is tiF in component Ito (millisecond).
C ALGBRC(52) is tiS in component Ito (millisecond).
C ALGBRC(72) is AiF in component Ito (dimensionless).
C ALGBRC(73) is AiS in component Ito (dimensionless).
C STATES(20) is iF in component Ito (dimensionless).
C STATES(21) is iS in component Ito (dimensionless).
C ALGBRC(74) is i in component Ito (dimensionless).
C ALGBRC(33) is assp in component Ito (dimensionless).
C STATES(22) is ap in component Ito (dimensionless).
C ALGBRC(55) is dti_develop in component Ito (dimensionless).
C ALGBRC(58) is dti_recover in component Ito (dimensionless).
C ALGBRC(61) is tiFp in component Ito (millisecond).
C ALGBRC(62) is tiSp in component Ito (millisecond).
C STATES(23) is iFp in component Ito (dimensionless).
C STATES(24) is iSp in component Ito (dimensionless).
C ALGBRC(75) is ip in component Ito (dimensionless).
C CONSTS(120) is Gto in component Ito (milliS_per_microF).
C ALGBRC(76) is fItop in component Ito (dimensionless).
C CONSTS(47) is Kmn in component ICaL (millimolar).
C CONSTS(48) is k2n in component ICaL (per_millisecond).
C CONSTS(49) is PCa_b in component ICaL (dimensionless).
C ALGBRC(8) is dss in component ICaL (dimensionless).
C STATES(25) is d in component ICaL (dimensionless).
C ALGBRC(9) is fss in component ICaL (dimensionless).
C CONSTS(50) is Aff in component ICaL (dimensionless).
C CONSTS(121) is Afs in component ICaL (dimensionless).
C STATES(26) is ff in component ICaL (dimensionless).
C STATES(27) is fs in component ICaL (dimensionless).
C ALGBRC(78) is f in component ICaL (dimensionless).
C ALGBRC(20) is fcass in component ICaL (dimensionless).
C ALGBRC(10) is jcass in component ICaL (dimensionless).
C ALGBRC(79) is Afcaf in component ICaL (dimensionless).
C ALGBRC(80) is Afcas in component ICaL (dimensionless).
C STATES(28) is fcaf in component ICaL (dimensionless).
C STATES(29) is fcas in component ICaL (dimensionless).
C ALGBRC(81) is fca in component ICaL (dimensionless).
C STATES(30) is jca in component ICaL (dimensionless).
C STATES(31) is ffp in component ICaL (dimensionless).
C ALGBRC(82) is fp in component ICaL (dimensionless).
C STATES(32) is fcafp in component ICaL (dimensionless).
C ALGBRC(83) is fcap in component ICaL (dimensionless).
C ALGBRC(11) is km2n in component ICaL (per_millisecond).
C ALGBRC(21) is anca_ss in component ICaL (dimensionless).
C STATES(33) is nca_ss in component ICaL (dimensionless).
C ALGBRC(22) is anca_i in component ICaL (dimensionless).
C STATES(34) is nca_i in component ICaL (dimensionless).
C ALGBRC(88) is PhiCaL_ss in component ICaL (dimensionless).
C ALGBRC(89) is PhiCaNa_ss in component ICaL (dimensionless).
C ALGBRC(90) is PhiCaK_ss in component ICaL (dimensionless).
C ALGBRC(107) is PhiCaL_i in component ICaL (dimensionless).
C ALGBRC(108) is PhiCaNa_i in component ICaL (dimensionless).
C ALGBRC(109) is PhiCaK_i in component ICaL (dimensionless).
C CONSTS(122) is PCa in component ICaL (dimensionless).
C CONSTS(132) is PCap in component ICaL (dimensionless).
C CONSTS(133) is PCaNa in component ICaL (dimensionless).
C CONSTS(134) is PCaK in component ICaL (dimensionless).
C CONSTS(138) is PCaNap in component ICaL (dimensionless).
C CONSTS(139) is PCaKp in component ICaL (dimensionless).
C ALGBRC(91) is fICaLp in component ICaL (dimensionless).
C ALGBRC(23) is td in component ICaL (millisecond).
C ALGBRC(24) is tff in component ICaL (millisecond).
C ALGBRC(25) is tfs in component ICaL (millisecond).
C ALGBRC(34) is tfcaf in component ICaL (millisecond).
C ALGBRC(35) is tfcas in component ICaL (millisecond).
C CONSTS(51) is tjca in component ICaL (millisecond).
C ALGBRC(36) is tffp in component ICaL (millisecond).
C ALGBRC(41) is tfcafp in component ICaL (millisecond).
C CONSTS(52) is vShift in component ICaL (millivolt).
C CONSTS(53) is offset in component ICaL (millisecond).
C CONSTS(123) is Io in component ICaL (dimensionless).
C ALGBRC(84) is Iss in component ICaL (dimensionless).
C ALGBRC(99) is Ii in component ICaL (dimensionless).
C CONSTS(54) is dielConstant in component ICaL (per_kelvin).
C CONSTS(135) is constA in component ICaL (dimensionless).
C CONSTS(140) is gamma_cao in component ICaL (dimensionless).
C ALGBRC(85) is gamma_cass in component ICaL (dimensionless).
C ALGBRC(102) is gamma_cai in component ICaL (dimensionless).
C CONSTS(141) is gamma_nao in component ICaL (dimensionless).
C ALGBRC(86) is gamma_nass in component ICaL (dimensionless).
C ALGBRC(105) is gamma_nai in component ICaL (dimensionless).
C CONSTS(142) is gamma_ko in component ICaL (dimensionless).
C ALGBRC(87) is gamma_kss in component ICaL (dimensionless).
C ALGBRC(106) is gamma_ki in component ICaL (dimensionless).
C CONSTS(55) is ICaL_fractionSS in component ICaL (dimensionless).
C CONSTS(56) is GKr_b in component IKr (milliS_per_microF).
C STATES(35) is C1 in component IKr (dimensionless).
C STATES(36) is C2 in component IKr (dimensionless).
C STATES(37) is C3 in component IKr (dimensionless).
C STATES(38) is I in component IKr (dimensionless).
C STATES(39) is O in component IKr (dimensionless).
C ALGBRC(42) is alpha in component IKr (per_millisecond).
C ALGBRC(48) is beta in component IKr (per_millisecond).
C CONSTS(57) is alpha_1 in component IKr (per_millisecond).
C CONSTS(58) is beta_1 in component IKr (per_millisecond).
C ALGBRC(43) is alpha_2 in component IKr (per_millisecond).
C ALGBRC(49) is beta_2 in component IKr (per_millisecond).
C ALGBRC(53) is alpha_i in component IKr (per_millisecond).
C ALGBRC(56) is beta_i in component IKr (per_millisecond).
C ALGBRC(59) is alpha_C2ToI in component IKr (per_millisecond).
C ALGBRC(63) is beta_ItoC2 in component IKr (per_millisecond).
C CONSTS(124) is GKr in component IKr (milliS_per_microF).
C CONSTS(59) is GKs_b in component IKs (milliS_per_microF).
C CONSTS(125) is GKs in component IKs (milliS_per_microF).
C ALGBRC(13) is xs1ss in component IKs (dimensionless).
C ALGBRC(27) is xs2ss in component IKs (dimensionless).
C ALGBRC(28) is txs1 in component IKs (millisecond).
C STATES(40) is xs1 in component IKs (dimensionless).
C STATES(41) is xs2 in component IKs (dimensionless).
C ALGBRC(117) is KsCa in component IKs (dimensionless).
C ALGBRC(37) is txs2 in component IKs (millisecond).
C CONSTS(126) is GK1 in component IK1 (milliS_per_microF).
C CONSTS(60) is GK1_b in component IK1 (milliS_per_microF).
C ALGBRC(119) is aK1 in component IK1 (dimensionless).
C ALGBRC(120) is bK1 in component IK1 (dimensionless).
C ALGBRC(121) is K1ss in component IK1 (dimensionless).
C CONSTS(61) is INaCa_fractionSS in component INaCa (dimensionless).
C CONSTS(62) is kna1 in component INaCa (per_millisecond).
C CONSTS(63) is kna2 in component INaCa (per_millisecond).
C CONSTS(64) is kna3 in component INaCa (per_millisecond).
C CONSTS(65) is kasymm in component INaCa (dimensionless).
C CONSTS(66) is wna in component INaCa (dimensionless).
C CONSTS(67) is wca in component INaCa (dimensionless).
C CONSTS(68) is wnaca in component INaCa (dimensionless).
C CONSTS(69) is kcaon in component INaCa (per_millisecond).
C CONSTS(70) is kcaoff in component INaCa (per_millisecond).
C CONSTS(71) is qna in component INaCa (dimensionless).
C CONSTS(72) is qca in component INaCa (dimensionless).
C ALGBRC(124) is hna in component INaCa (dimensionless).
C ALGBRC(123) is hca in component INaCa (dimensionless).
C CONSTS(73) is KmCaAct in component INaCa (millimolar).
C CONSTS(74) is Gncx_b in component INaCa (milliS_per_microF).
C CONSTS(153) is Gncx in component INaCa (milliS_per_microF).
C ALGBRC(125) is h1_i in component INaCa (dimensionless).
C ALGBRC(126) is h2_i in component INaCa (dimensionless).
C ALGBRC(127) is h3_i in component INaCa (dimensionless).
C ALGBRC(128) is h4_i in component INaCa (dimensionless).
C ALGBRC(129) is h5_i in component INaCa (dimensionless).
C ALGBRC(130) is h6_i in component INaCa (dimensionless).
C ALGBRC(131) is h7_i in component INaCa (dimensionless).
C ALGBRC(132) is h8_i in component INaCa (dimensionless).
C ALGBRC(133) is h9_i in component INaCa (dimensionless).
C CONSTS(147) is h10_i in component INaCa (dimensionless).
C CONSTS(148) is h11_i in component INaCa (dimensionless).
C CONSTS(149) is h12_i in component INaCa (dimensionless).
C CONSTS(150) is k1_i in component INaCa (dimensionless).
C CONSTS(151) is k2_i in component INaCa (dimensionless).
C ALGBRC(134) is k3p_i in component INaCa (dimensionless).
C ALGBRC(135) is k3pp_i in component INaCa (dimensionless).
C ALGBRC(136) is k3_i in component INaCa (dimensionless).
C ALGBRC(139) is k4_i in component INaCa (dimensionless).
C ALGBRC(137) is k4p_i in component INaCa (dimensionless).
C ALGBRC(138) is k4pp_i in component INaCa (dimensionless).
C CONSTS(152) is k5_i in component INaCa (dimensionless).
C ALGBRC(140) is k6_i in component INaCa (dimensionless).
C ALGBRC(141) is k7_i in component INaCa (dimensionless).
C ALGBRC(142) is k8_i in component INaCa (dimensionless).
C ALGBRC(143) is x1_i in component INaCa (dimensionless).
C ALGBRC(144) is x2_i in component INaCa (dimensionless).
C ALGBRC(145) is x3_i in component INaCa (dimensionless).
C ALGBRC(146) is x4_i in component INaCa (dimensionless).
C ALGBRC(147) is E1_i in component INaCa (dimensionless).
C ALGBRC(148) is E2_i in component INaCa (dimensionless).
C ALGBRC(149) is E3_i in component INaCa (dimensionless).
C ALGBRC(150) is E4_i in component INaCa (dimensionless).
C ALGBRC(151) is allo_i in component INaCa (dimensionless).
C ALGBRC(152) is JncxNa_i in component INaCa (millimolar_per_millisecond).
C ALGBRC(153) is JncxCa_i in component INaCa (millimolar_per_millisecond).
C ALGBRC(155) is h1_ss in component INaCa (dimensionless).
C ALGBRC(156) is h2_ss in component INaCa (dimensionless).
C ALGBRC(157) is h3_ss in component INaCa (dimensionless).
C ALGBRC(158) is h4_ss in component INaCa (dimensionless).
C ALGBRC(159) is h5_ss in component INaCa (dimensionless).
C ALGBRC(160) is h6_ss in component INaCa (dimensionless).
C ALGBRC(161) is h7_ss in component INaCa (dimensionless).
C ALGBRC(162) is h8_ss in component INaCa (dimensionless).
C ALGBRC(163) is h9_ss in component INaCa (dimensionless).
C CONSTS(154) is h10_ss in component INaCa (dimensionless).
C CONSTS(155) is h11_ss in component INaCa (dimensionless).
C CONSTS(156) is h12_ss in component INaCa (dimensionless).
C CONSTS(157) is k1_ss in component INaCa (dimensionless).
C CONSTS(158) is k2_ss in component INaCa (dimensionless).
C ALGBRC(164) is k3p_ss in component INaCa (dimensionless).
C ALGBRC(165) is k3pp_ss in component INaCa (dimensionless).
C ALGBRC(166) is k3_ss in component INaCa (dimensionless).
C ALGBRC(169) is k4_ss in component INaCa (dimensionless).
C ALGBRC(167) is k4p_ss in component INaCa (dimensionless).
C ALGBRC(168) is k4pp_ss in component INaCa (dimensionless).
C CONSTS(159) is k5_ss in component INaCa (dimensionless).
C ALGBRC(170) is k6_ss in component INaCa (dimensionless).
C ALGBRC(171) is k7_ss in component INaCa (dimensionless).
C ALGBRC(172) is k8_ss in component INaCa (dimensionless).
C ALGBRC(173) is x1_ss in component INaCa (dimensionless).
C ALGBRC(174) is x2_ss in component INaCa (dimensionless).
C ALGBRC(175) is x3_ss in component INaCa (dimensionless).
C ALGBRC(176) is x4_ss in component INaCa (dimensionless).
C ALGBRC(177) is E1_ss in component INaCa (dimensionless).
C ALGBRC(178) is E2_ss in component INaCa (dimensionless).
C ALGBRC(179) is E3_ss in component INaCa (dimensionless).
C ALGBRC(180) is E4_ss in component INaCa (dimensionless).
C ALGBRC(181) is allo_ss in component INaCa (dimensionless).
C ALGBRC(182) is JncxNa_ss in component INaCa (millimolar_per_millisecond).
C ALGBRC(183) is JncxCa_ss in component INaCa (millimolar_per_millisecond).
C CONSTS(75) is k1p in component INaK (per_millisecond).
C CONSTS(76) is k1m in component INaK (per_millisecond).
C CONSTS(77) is k2p in component INaK (per_millisecond).
C CONSTS(78) is k2m in component INaK (per_millisecond).
C CONSTS(79) is k3p in component INaK (per_millisecond).
C CONSTS(80) is k3m in component INaK (per_millisecond).
C CONSTS(81) is k4p in component INaK (per_millisecond).
C CONSTS(82) is k4m in component INaK (per_millisecond).
C CONSTS(83) is Knai0 in component INaK (millimolar).
C CONSTS(84) is Knao0 in component INaK (millimolar).
C CONSTS(85) is delta in component INaK (millivolt).
C CONSTS(86) is Kki in component INaK (per_millisecond).
C CONSTS(87) is Kko in component INaK (per_millisecond).
C CONSTS(88) is MgADP in component INaK (millimolar).
C CONSTS(89) is MgATP in component INaK (millimolar).
C CONSTS(90) is Kmgatp in component INaK (millimolar).
C CONSTS(91) is H in component INaK (millimolar).
C CONSTS(92) is eP in component INaK (dimensionless).
C CONSTS(93) is Khp in component INaK (millimolar).
C CONSTS(94) is Knap in component INaK (millimolar).
C CONSTS(95) is Kxkur in component INaK (millimolar).
C CONSTS(96) is Pnak_b in component INaK (milliS_per_microF).
C CONSTS(163) is Pnak in component INaK (milliS_per_microF).
C ALGBRC(185) is Knai in component INaK (millimolar).
C ALGBRC(186) is Knao in component INaK (millimolar).
C ALGBRC(187) is P in component INaK (dimensionless).
C ALGBRC(188) is a1 in component INaK (dimensionless).
C CONSTS(160) is b1 in component INaK (dimensionless).
C CONSTS(161) is a2 in component INaK (dimensionless).
C ALGBRC(189) is b2 in component INaK (dimensionless).
C ALGBRC(190) is a3 in component INaK (dimensionless).
C ALGBRC(191) is b3 in component INaK (dimensionless).
C CONSTS(162) is a4 in component INaK (dimensionless).
C ALGBRC(192) is b4 in component INaK (dimensionless).
C ALGBRC(193) is x1 in component INaK (dimensionless).
C ALGBRC(194) is x2 in component INaK (dimensionless).
C ALGBRC(195) is x3 in component INaK (dimensionless).
C ALGBRC(196) is x4 in component INaK (dimensionless).
C ALGBRC(197) is E1 in component INaK (dimensionless).
C ALGBRC(198) is E2 in component INaK (dimensionless).
C ALGBRC(199) is E3 in component INaK (dimensionless).
C ALGBRC(200) is E4 in component INaK (dimensionless).
C ALGBRC(201) is JnakNa in component INaK (millimolar_per_millisecond).
C ALGBRC(202) is JnakK in component INaK (millimolar_per_millisecond).
C ALGBRC(204) is xkb in component IKb (dimensionless).
C CONSTS(97) is GKb_b in component IKb (milliS_per_microF).
C CONSTS(127) is GKb in component IKb (milliS_per_microF).
C CONSTS(98) is PNab in component INab (milliS_per_microF).
C CONSTS(99) is PCab in component ICab (milliS_per_microF).
C CONSTS(100) is GpCa in component IpCa (milliS_per_microF).
C CONSTS(101) is KmCap in component IpCa (millimolar).
C CONSTS(102) is GClCa in component ICl (milliS_per_microF).
C CONSTS(103) is GClb in component ICl (milliS_per_microF).
C CONSTS(104) is KdClCa in component ICl (millimolar).
C CONSTS(105) is Fjunc in component ICl (dimensionless).
C ALGBRC(211) is IClCa_junc in component ICl (microA_per_microF).
C ALGBRC(213) is IClCa_sl in component ICl (microA_per_microF).
C CONSTS(106) is tauNa in component diff (millisecond).
C CONSTS(107) is tauK in component diff (millisecond).
C CONSTS(108) is tauCa in component diff (millisecond).
C CONSTS(109) is bt in component ryr (millisecond).
C CONSTS(128) is a_rel in component ryr (millimolar_per_millisecond).
C ALGBRC(94) is Jrel_inf_b in component ryr (millimolar_per_millisecond).
C ALGBRC(97) is Jrel_inf in component ryr (millimolar_per_millisecond).
C ALGBRC(100) is tau_rel_b in component ryr (millisecond).
C ALGBRC(103) is tau_rel in component ryr (millisecond).
C STATES(42) is Jrel_np in component ryr (millimolar_per_millisecond).
C CONSTS(129) is btp in component ryr (millisecond).
C CONSTS(136) is a_relp in component ryr (millimolar_per_millisecond).
C ALGBRC(95) is Jrel_infp_b in component ryr (millimolar_per_millisecond).
C ALGBRC(98) is Jrel_infp in component ryr (millimolar_per_millisecond).
C ALGBRC(101) is tau_relp_b in component ryr (millisecond).
C ALGBRC(104) is tau_relp in component ryr (millisecond).
C STATES(43) is Jrel_p in component ryr (millimolar_per_millisecond).
C CONSTS(110) is cajsr_half in component ryr (millimolar).
C ALGBRC(214) is fJrelp in component ryr (dimensionless).
C CONSTS(111) is Jrel_b in component ryr (dimensionless).
C CONSTS(130) is upScale in component SERCA (dimensionless).
C ALGBRC(218) is Jupnp in component SERCA (millimolar_per_millisecond).
C ALGBRC(219) is Jupp in component SERCA (millimolar_per_millisecond).
C ALGBRC(220) is fJupp in component SERCA (dimensionless).
C ALGBRC(221) is Jleak in component SERCA (millimolar_per_millisecond).
C CONSTS(112) is Jup_b in component SERCA (dimensionless).
C RATES(1) is d/dt v in component membrane (millivolt).
C RATES(2) is d/dt CaMKt in component CaMK (millimolar).
C RATES(4) is d/dt nai in component intracellular_ions (millimolar).
C RATES(5) is d/dt nass in component intracellular_ions (millimolar).
C RATES(6) is d/dt ki in component intracellular_ions (millimolar).
C RATES(7) is d/dt kss in component intracellular_ions (millimolar).
C RATES(10) is d/dt cai in component intracellular_ions (millimolar).
C RATES(3) is d/dt cass in component intracellular_ions (millimolar).
C RATES(8) is d/dt cansr in component intracellular_ions (millimolar).
C RATES(9) is d/dt cajsr in component intracellular_ions (millimolar).
C RATES(11) is d/dt m in component INa (dimensionless).
C RATES(12) is d/dt h in component INa (dimensionless).
C RATES(13) is d/dt j in component INa (dimensionless).
C RATES(14) is d/dt hp in component INa (dimensionless).
C RATES(15) is d/dt jp in component INa (dimensionless).
C RATES(16) is d/dt mL in component INaL (dimensionless).
C RATES(17) is d/dt hL in component INaL (dimensionless).
C RATES(18) is d/dt hLp in component INaL (dimensionless).
C RATES(19) is d/dt a in component Ito (dimensionless).
C RATES(20) is d/dt iF in component Ito (dimensionless).
C RATES(21) is d/dt iS in component Ito (dimensionless).
C RATES(22) is d/dt ap in component Ito (dimensionless).
C RATES(23) is d/dt iFp in component Ito (dimensionless).
C RATES(24) is d/dt iSp in component Ito (dimensionless).
C RATES(25) is d/dt d in component ICaL (dimensionless).
C RATES(26) is d/dt ff in component ICaL (dimensionless).
C RATES(27) is d/dt fs in component ICaL (dimensionless).
C RATES(28) is d/dt fcaf in component ICaL (dimensionless).
C RATES(29) is d/dt fcas in component ICaL (dimensionless).
C RATES(30) is d/dt jca in component ICaL (dimensionless).
C RATES(31) is d/dt ffp in component ICaL (dimensionless).
C RATES(32) is d/dt fcafp in component ICaL (dimensionless).
C RATES(33) is d/dt nca_ss in component ICaL (dimensionless).
C RATES(34) is d/dt nca_i in component ICaL (dimensionless).
C RATES(37) is d/dt C3 in component IKr (dimensionless).
C RATES(36) is d/dt C2 in component IKr (dimensionless).
C RATES(35) is d/dt C1 in component IKr (dimensionless).
C RATES(39) is d/dt O in component IKr (dimensionless).
C RATES(38) is d/dt I in component IKr (dimensionless).
C RATES(40) is d/dt xs1 in component IKs (dimensionless).
C RATES(41) is d/dt xs2 in component IKs (dimensionless).
C RATES(42) is d/dt Jrel_np in component ryr (millimolar_per_millisecond).
C RATES(43) is d/dt Jrel_p in component ryr (millimolar_per_millisecond).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 140.0
      CONSTS(3) = 1.8
      CONSTS(4) = 5.0
      CONSTS(5) = 150.0
      CONSTS(6) = 8314
      CONSTS(7) = 310
      CONSTS(8) = 96485
      CONSTS(9) = 1
      CONSTS(10) = 2
      CONSTS(11) = 1
      CONSTS(12) = -1
      CONSTS(13) = 0.01
      CONSTS(14) = 0.0011
      STATES(1) = -89.1400
      CONSTS(15) = 0
      CONSTS(16) = 100000000000000000
      CONSTS(17) = -53
      CONSTS(18) = 1000
      CONSTS(19) = 1.0
      CONSTS(20) = 0.15
      CONSTS(21) = 0.05
      CONSTS(22) = 0.00068
      CONSTS(23) = 0.05
      CONSTS(24) = 0.0015
      STATES(2) = 0.0129
      STATES(3) = 5.7672e-05
      CONSTS(25) = 0.05
      CONSTS(26) = 0.00238
      CONSTS(27) = 0.07
      CONSTS(28) = 0.0005
      CONSTS(29) = 0.047
      CONSTS(30) = 0.00087
      CONSTS(31) = 1.124
      CONSTS(32) = 0.0087
      CONSTS(33) = 10
      CONSTS(34) = 0.8
      STATES(4) = 12.8363
      STATES(5) = 12.8366
      STATES(6) = 142.6951
      STATES(7) = 142.6951
      STATES(8) = 1.8119
      STATES(9) = 1.8102
      STATES(10) = 6.6309e-05
      CONSTS(35) = 24.0
      CONSTS(36) = 0.01833
      CONSTS(37) = 4.3195
      CONSTS(38) = 0.0
      CONSTS(39) = 5
      CONSTS(40) = 2
      CONSTS(41) = 0.25
      STATES(11) = 0.00074303
      STATES(12) = 0.8360
      STATES(13) = 0.8359
      STATES(14) = 0.6828
      STATES(15) = 0.8357
      CONSTS(42) = 11.7802
      STATES(16) = 0.00015166
      CONSTS(43) = 200
      STATES(17) = 0.5401
      STATES(18) = 0.3034
      CONSTS(44) = 0.0279
      CONSTS(45) = 0.16
      STATES(19) = 0.00092716
      CONSTS(46) = 0
      STATES(20) = 0.9996
      STATES(21) = 0.9996
      STATES(22) = 0.0004724
      STATES(23) = 0.9996
      STATES(24) = 0.9996
      CONSTS(47) = 0.002
      CONSTS(48) = 500
      CONSTS(49) = 8.3757e-05
      STATES(25) = 0
      CONSTS(50) = 0.6
      STATES(26) = 1.0
      STATES(27) = 0.9485
      STATES(28) = 1.0
      STATES(29) = 0.9999
      STATES(30) = 1.0
      STATES(31) = 1.0
      STATES(32) = 1.0
      STATES(33) = 0.00030853
      STATES(34) = 0.00053006
      CONSTS(51) = 75
      CONSTS(52) = 0
      CONSTS(53) = 0
      CONSTS(54) = 74
      CONSTS(55) = 0.8
      CONSTS(56) = 0.0321
      STATES(35) = 0.00067941
      STATES(36) = 0.00082869
      STATES(37) = 0.9982
      STATES(38) = 9.5416e-06
      STATES(39) = 0.00027561
      CONSTS(57) = 0.154375
      CONSTS(58) = 0.1911
      CONSTS(59) = 0.0011
      STATES(40) = 0.2309
      STATES(41) = 0.00016975
      CONSTS(60) = 0.6992
      CONSTS(61) = 0.35
      CONSTS(62) = 15
      CONSTS(63) = 5
      CONSTS(64) = 88.12
      CONSTS(65) = 12.5
      CONSTS(66) = 6e4
      CONSTS(67) = 6e4
      CONSTS(68) = 5e3
      CONSTS(69) = 1.5e6
      CONSTS(70) = 5e3
      CONSTS(71) = 0.5224
      CONSTS(72) = 0.167
      CONSTS(73) = 150e-6
      CONSTS(74) = 0.0034
      CONSTS(75) = 949.5
      CONSTS(76) = 182.4
      CONSTS(77) = 687.2
      CONSTS(78) = 39.4
      CONSTS(79) = 1899
      CONSTS(80) = 79300
      CONSTS(81) = 639
      CONSTS(82) = 40
      CONSTS(83) = 9.073
      CONSTS(84) = 27.78
      CONSTS(85) = -0.155
      CONSTS(86) = 0.5
      CONSTS(87) = 0.3582
      CONSTS(88) = 0.05
      CONSTS(89) = 9.8
      CONSTS(90) = 1.698e-7
      CONSTS(91) = 1e-7
      CONSTS(92) = 4.2
      CONSTS(93) = 1.698e-7
      CONSTS(94) = 224
      CONSTS(95) = 292
      CONSTS(96) = 15.4509
      CONSTS(97) = 0.0189
      CONSTS(98) = 1.9239e-09
      CONSTS(99) = 5.9194e-08
      CONSTS(100) = 5e-04
      CONSTS(101) = 0.0005
      CONSTS(102) = 0.2843
      CONSTS(103) = 1.98e-3
      CONSTS(104) = 0.1
      CONSTS(105) = 1
      CONSTS(106) = 2.0
      CONSTS(107) = 2.0
      CONSTS(108) = 0.2
      CONSTS(109) = 4.75
      STATES(42) = 2.8189e-24
      STATES(43) = 0
      CONSTS(110) = 1.7
      CONSTS(111) = 1.5378
      CONSTS(112) = 1.0
      CONSTS(113) =  1000.00*3.14000*CONSTS(14)*CONSTS(14)*CONSTS(13)
      CONSTS(114) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(25)*1.30000, CONSTS(25))
      CONSTS(115) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(12)*CONSTS(8)))*log(CONSTS(5)/CONSTS(35))
      CONSTS(116) = CONSTS(4)/CONSTS(39) ** 0.240000
      CONSTS(117) = 1.00000/(1.00000+CONSTS(40)/CONSTS(41) ** 2.00000)
      CONSTS(118) =  3.00000*CONSTS(43)
      CONSTS(119) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(44)*0.600000, CONSTS(44))
      CONSTS(120) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(45)*2.00000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(45)*2.00000, CONSTS(45))
      CONSTS(121) = 1.00000 - CONSTS(50)
      CONSTS(122) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(49)*1.20000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(49)*2.00000, CONSTS(49))
      CONSTS(123) = ( 0.500000*(CONSTS(2)+CONSTS(4)+CONSTS(5)+ 4.00000*CONSTS(3)))/1000.00
      CONSTS(124) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(56)*1.30000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(56)*0.800000, CONSTS(56))
      CONSTS(125) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(59)*1.40000, CONSTS(59))
      CONSTS(126) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(60)*1.20000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(60)*1.30000, CONSTS(60))
      CONSTS(127) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(97)*0.600000, CONSTS(97))
      CONSTS(128) = ( 0.500000*CONSTS(109))/1.00000
      CONSTS(129) =  1.25000*CONSTS(109)
      CONSTS(130) = TERNRY(CONSTS(1).EQ.1.00000, 1.30000, 1.00000)
      CONSTS(131) =  2.00000*3.14000*CONSTS(14)*CONSTS(14)+ 2.00000*3.14000*CONSTS(14)*CONSTS(13)
      CONSTS(132) =  1.10000*CONSTS(122)
      CONSTS(133) =  0.00125000*CONSTS(122)
      CONSTS(134) =  0.000357400*CONSTS(122)
      CONSTS(135) =  1.82000e+06* CONSTS(54)*CONSTS(7) ** - 1.50000
      CONSTS(136) = ( 0.500000*CONSTS(129))/1.00000
      CONSTS(137) =  2.00000*CONSTS(131)
      CONSTS(138) =  0.00125000*CONSTS(132)
      CONSTS(139) =  0.000357400*CONSTS(132)
      CONSTS(140) = EXP( - CONSTS(135)*4.00000*( CONSTS(123) ** (1.0 / 2)/(1.00000+ CONSTS(123) ** (1.0 / 2)) -  0.300000*CONSTS(123)))
      CONSTS(141) = EXP( - CONSTS(135)*1.00000*( CONSTS(123) ** (1.0 / 2)/(1.00000+ CONSTS(123) ** (1.0 / 2)) -  0.300000*CONSTS(123)))
      CONSTS(142) = EXP( - CONSTS(135)*1.00000*( CONSTS(123) ** (1.0 / 2)/(1.00000+ CONSTS(123) ** (1.0 / 2)) -  0.300000*CONSTS(123)))
      CONSTS(143) =  0.680000*CONSTS(113)
      CONSTS(144) =  0.0552000*CONSTS(113)
      CONSTS(145) =  0.00480000*CONSTS(113)
      CONSTS(146) =  0.0200000*CONSTS(113)
      CONSTS(147) = CONSTS(65)+1.00000+ (CONSTS(2)/CONSTS(62))*(1.00000+CONSTS(2)/CONSTS(63))
      CONSTS(148) = ( CONSTS(2)*CONSTS(2))/( CONSTS(147)*CONSTS(62)*CONSTS(63))
      CONSTS(149) = 1.00000/CONSTS(147)
      CONSTS(150) =  CONSTS(149)*CONSTS(3)*CONSTS(69)
      CONSTS(151) = CONSTS(70)
      CONSTS(152) = CONSTS(70)
      CONSTS(153) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(74)*1.10000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(74)*1.40000, CONSTS(74))
      CONSTS(154) = CONSTS(65)+1.00000+ (CONSTS(2)/CONSTS(62))*(1.00000+CONSTS(2)/CONSTS(63))
      CONSTS(155) = ( CONSTS(2)*CONSTS(2))/( CONSTS(154)*CONSTS(62)*CONSTS(63))
      CONSTS(156) = 1.00000/CONSTS(154)
      CONSTS(157) =  CONSTS(156)*CONSTS(3)*CONSTS(69)
      CONSTS(158) = CONSTS(70)
      CONSTS(159) = CONSTS(70)
      CONSTS(160) =  CONSTS(76)*CONSTS(88)
      CONSTS(161) = CONSTS(77)
      CONSTS(162) = (( CONSTS(81)*CONSTS(89))/CONSTS(90))/(1.00000+CONSTS(89)/CONSTS(90))
      CONSTS(163) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(96)*0.900000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(96)*0.700000, CONSTS(96))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
      RATES(17) = (ALGBRC(4) - STATES(17))/CONSTS(43)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
      RATES(18) = (ALGBRC(5) - STATES(18))/CONSTS(118)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+18.0800)/2.79160))
      RATES(30) = (ALGBRC(10) - STATES(30))/CONSTS(51)
      ALGBRC(1) = 1.00000/1.00000+EXP(- (STATES(1)+56.8600)/9.03000) ** 2.00000
      ALGBRC(14) =  0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
      RATES(11) = (ALGBRC(1) - STATES(11))/ALGBRC(14)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(17) =  0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
      RATES(16) = (ALGBRC(3) - STATES(16))/ALGBRC(17)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 14.3400)/14.8200))
      ALGBRC(18) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+CONSTS(46)+100.000)/29.3814)))
      RATES(19) = (ALGBRC(6) - STATES(19))/ALGBRC(18)
      ALGBRC(8) = TERNRY(STATES(1).GE.31.4978, 1.00000,  1.07630*EXP( - 1.00700*EXP( - 0.0829000*STATES(1))))
      ALGBRC(23) = CONSTS(53)+0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+CONSTS(52)+6.00000))+EXP( 0.0900000*(STATES(1)+CONSTS(52)+14.0000)))
      RATES(25) = (ALGBRC(8) - STATES(25))/ALGBRC(23)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
      ALGBRC(24) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000)/10.0000))
      RATES(26) = (ALGBRC(9) - STATES(26))/ALGBRC(24)
      ALGBRC(25) = 1000.00+1.00000/( 3.50000e-05*EXP(- (STATES(1)+5.00000)/4.00000)+ 3.50000e-05*EXP((STATES(1)+5.00000)/6.00000))
      RATES(27) = (ALGBRC(9) - STATES(27))/ALGBRC(25)
      ALGBRC(11) =  STATES(30)*1.00000
      ALGBRC(21) = 1.00000/(CONSTS(48)/ALGBRC(11)+1.00000+CONSTS(47)/STATES(3) ** 4.00000)
      RATES(33) =  ALGBRC(21)*CONSTS(48) -  STATES(33)*ALGBRC(11)
      ALGBRC(22) = 1.00000/(CONSTS(48)/ALGBRC(11)+1.00000+CONSTS(47)/STATES(10) ** 4.00000)
      RATES(34) =  ALGBRC(22)*CONSTS(48) -  STATES(34)*ALGBRC(11)
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(28) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      RATES(40) = (ALGBRC(13) - STATES(40))/ALGBRC(28)
      ALGBRC(33) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 24.3400)/14.8200))
      RATES(22) = (ALGBRC(33) - STATES(22))/ALGBRC(18)
      ALGBRC(20) = ALGBRC(9)
      ALGBRC(34) = 7.00000+1.00000/( 0.0400000*EXP(- (STATES(1) - 4.00000)/7.00000)+ 0.0400000*EXP((STATES(1) - 4.00000)/7.00000))
      RATES(28) = (ALGBRC(20) - STATES(28))/ALGBRC(34)
      ALGBRC(35) = 100.000+1.00000/( 0.000120000*EXP(- STATES(1)/3.00000)+ 0.000120000*EXP(STATES(1)/7.00000))
      RATES(29) = (ALGBRC(20) - STATES(29))/ALGBRC(35)
      ALGBRC(36) =  2.50000*ALGBRC(24)
      RATES(31) = (ALGBRC(9) - STATES(31))/ALGBRC(36)
      ALGBRC(27) = ALGBRC(13)
      ALGBRC(37) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      RATES(41) = (ALGBRC(27) - STATES(41))/ALGBRC(37)
      ALGBRC(44) = ( CONSTS(23)*(1.00000 - STATES(2)))/(1.00000+CONSTS(24)/STATES(3))
      RATES(2) =  CONSTS(21)*ALGBRC(44)*(ALGBRC(44)+STATES(2)) -  CONSTS(22)*STATES(2)
      ALGBRC(2) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
      ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.0570000*EXP(- (STATES(1)+80.0000)/6.80000))
      ALGBRC(30) = TERNRY(STATES(1).GE.- 40.0000, 0.770000/( 0.130000*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000))),  2.70000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.348500*STATES(1)))
      ALGBRC(38) = 1.00000/(ALGBRC(15)+ALGBRC(30))
      RATES(12) = (ALGBRC(2) - STATES(12))/ALGBRC(38)
      ALGBRC(41) =  2.50000*ALGBRC(34)
      RATES(32) = (ALGBRC(20) - STATES(32))/ALGBRC(41)
      ALGBRC(39) = ALGBRC(2)
      ALGBRC(16) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 25428.0*EXP( 0.244400*STATES(1)) -  6.94800e-06*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
      ALGBRC(31) = TERNRY(STATES(1).GE.- 40.0000, ( 0.600000*EXP( 0.0570000*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.0242400*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
      ALGBRC(45) = 1.00000/(ALGBRC(16)+ALGBRC(31))
      RATES(13) = (ALGBRC(39) - STATES(13))/ALGBRC(45)
      ALGBRC(46) = 1.00000/1.00000+EXP((STATES(1)+77.5500)/7.43000) ** 2.00000
      RATES(14) = (ALGBRC(46) - STATES(14))/ALGBRC(38)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(46)+43.9400)/5.71100))
      ALGBRC(19) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+CONSTS(46)+70.0000)/5.00000)), 1.00000)
      ALGBRC(32) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+CONSTS(46)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+CONSTS(46)+50.0000)/16.5900))
      ALGBRC(47) =  ALGBRC(32)*ALGBRC(19)
      RATES(20) = (ALGBRC(7) - STATES(20))/ALGBRC(47)
      ALGBRC(29) = ( STATES(1)*CONSTS(8))/( CONSTS(6)*CONSTS(7))
      ALGBRC(42) =  0.116100*EXP( 0.299000*ALGBRC(29))
      ALGBRC(48) =  0.244200*EXP( - 1.60400*ALGBRC(29))
      RATES(37) =  ALGBRC(48)*STATES(36) -  ALGBRC(42)*STATES(37)
      RATES(36) = ( ALGBRC(42)*STATES(37)+ CONSTS(58)*STATES(35)) -  (ALGBRC(48)+CONSTS(57))*STATES(36)
      ALGBRC(51) =  1.46000*ALGBRC(45)
      RATES(15) = (ALGBRC(39) - STATES(15))/ALGBRC(51)
      ALGBRC(40) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+CONSTS(46)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+CONSTS(46)+114.100)/8.07900))
      ALGBRC(52) =  ALGBRC(40)*ALGBRC(19)
      RATES(21) = (ALGBRC(7) - STATES(21))/ALGBRC(52)
      ALGBRC(43) =  0.0578000*EXP( 0.971000*ALGBRC(29))
      ALGBRC(49) =  0.000349000*EXP( - 1.06200*ALGBRC(29))
      ALGBRC(53) =  0.253300*EXP( 0.595300*ALGBRC(29))
      ALGBRC(56) =  0.0652500*EXP( - 0.820900*ALGBRC(29))
      RATES(39) = ( ALGBRC(43)*STATES(35)+ ALGBRC(56)*STATES(38)) -  (ALGBRC(49)+ALGBRC(53))*STATES(39)
      ALGBRC(55) = 1.35400+0.000100000/(EXP(((STATES(1)+CONSTS(46)) - 167.400)/15.8900)+EXP(- ((STATES(1)+CONSTS(46)) - 12.2300)/0.215400))
      ALGBRC(58) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+CONSTS(46)+70.0000)/20.0000))
      ALGBRC(61) =  ALGBRC(55)*ALGBRC(58)*ALGBRC(47)
      RATES(23) = (ALGBRC(7) - STATES(23))/ALGBRC(61)
      ALGBRC(62) =  ALGBRC(55)*ALGBRC(58)*ALGBRC(52)
      RATES(24) = (ALGBRC(7) - STATES(24))/ALGBRC(62)
      ALGBRC(59) =  5.20000e-05*EXP( 1.52500*ALGBRC(29))
      ALGBRC(63) = ( ALGBRC(49)*ALGBRC(56)*ALGBRC(59))/( ALGBRC(43)*ALGBRC(53))
      RATES(35) = ( CONSTS(57)*STATES(36)+ ALGBRC(49)*STATES(39)+ ALGBRC(63)*STATES(38)) -  (CONSTS(58)+ALGBRC(43)+ALGBRC(59))*STATES(35)
      RATES(38) = ( ALGBRC(59)*STATES(35)+ ALGBRC(53)*STATES(39)) -  (ALGBRC(63)+ALGBRC(56))*STATES(38)
      ALGBRC(78) =  CONSTS(50)*STATES(26)+ CONSTS(121)*STATES(27)
      ALGBRC(79) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(80) = 1.00000 - ALGBRC(79)
      ALGBRC(81) =  ALGBRC(79)*STATES(28)+ ALGBRC(80)*STATES(29)
      ALGBRC(82) =  CONSTS(50)*STATES(31)+ CONSTS(121)*STATES(27)
      ALGBRC(83) =  ALGBRC(79)*STATES(32)+ ALGBRC(80)*STATES(29)
      ALGBRC(26) = ( STATES(1)*CONSTS(8)*CONSTS(8))/( CONSTS(6)*CONSTS(7))
      ALGBRC(84) = ( 0.500000*(STATES(5)+STATES(7)+CONSTS(35)+ 4.00000*STATES(3)))/1000.00
      ALGBRC(85) = EXP( - CONSTS(135)*4.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(88) = ( 4.00000*ALGBRC(26)*( ALGBRC(85)*STATES(3)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(50) = ALGBRC(44)+STATES(2)
      ALGBRC(91) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(92) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(122)*ALGBRC(88)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(132)*ALGBRC(88)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      ALGBRC(94) = (( - CONSTS(128)*ALGBRC(92))/1.00000)/(1.00000+CONSTS(110)/STATES(9) ** 8.00000)
      ALGBRC(97) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(94)*1.70000, ALGBRC(94))
      ALGBRC(100) = CONSTS(109)/(1.00000+0.0123000/STATES(9))
      ALGBRC(103) = TERNRY(ALGBRC(100).LT.0.00100000, 0.00100000, ALGBRC(100))
      RATES(42) = (ALGBRC(97) - STATES(42))/ALGBRC(103)
      ALGBRC(95) = (( - CONSTS(136)*ALGBRC(92))/1.00000)/(1.00000+CONSTS(110)/STATES(9) ** 8.00000)
      ALGBRC(98) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(95)*1.70000, ALGBRC(95))
      ALGBRC(101) = CONSTS(129)/(1.00000+0.0123000/STATES(9))
      ALGBRC(104) = TERNRY(ALGBRC(101).LT.0.00100000, 0.00100000, ALGBRC(101))
      RATES(43) = (ALGBRC(98) - STATES(43))/ALGBRC(104)
      ALGBRC(65) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(11)*CONSTS(8)))*log(CONSTS(4)/STATES(6))
      ALGBRC(72) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(46)) - 213.600)/151.200))
      ALGBRC(73) = 1.00000 - ALGBRC(72)
      ALGBRC(74) =  ALGBRC(72)*STATES(20)+ ALGBRC(73)*STATES(21)
      ALGBRC(75) =  ALGBRC(72)*STATES(23)+ ALGBRC(73)*STATES(24)
      ALGBRC(76) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(77) =  CONSTS(120)*(STATES(1) - ALGBRC(65))*( (1.00000 - ALGBRC(76))*STATES(19)*ALGBRC(74)+ ALGBRC(76)*STATES(22)*ALGBRC(75))
      ALGBRC(116) =  CONSTS(124)* (CONSTS(4)/5.00000) ** (1.0 / 2)*STATES(39)*(STATES(1) - ALGBRC(65))
      ALGBRC(66) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(11)*CONSTS(8)))*log((CONSTS(4)+ CONSTS(36)*CONSTS(2))/(STATES(6)+ CONSTS(36)*STATES(4)))
      ALGBRC(117) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
      ALGBRC(118) =  CONSTS(125)*ALGBRC(117)*STATES(40)*STATES(41)*(STATES(1) - ALGBRC(66))
      ALGBRC(119) = 4.09400/(1.00000+EXP( 0.121700*((STATES(1) - ALGBRC(65)) - 49.9340)))
      ALGBRC(120) = ( 15.7200*EXP( 0.0674000*((STATES(1) - ALGBRC(65)) - 3.25700))+EXP( 0.0618000*((STATES(1) - ALGBRC(65)) - 594.310)))/(1.00000+EXP( - 0.162900*((STATES(1) - ALGBRC(65))+14.2070)))
      ALGBRC(121) = ALGBRC(119)/(ALGBRC(119)+ALGBRC(120))
      ALGBRC(122) =  CONSTS(126)* (CONSTS(4)/5.00000) ** (1.0 / 2)*ALGBRC(121)*(STATES(1) - ALGBRC(65))
      ALGBRC(186) =  CONSTS(84)*EXP(( (1.00000 - CONSTS(85))*ALGBRC(29))/3.00000)
      ALGBRC(190) = ( CONSTS(79)*CONSTS(4)/CONSTS(87) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(186) ** 3.00000+1.00000+CONSTS(4)/CONSTS(87) ** 2.00000) - 1.00000)
      ALGBRC(187) = CONSTS(92)/(1.00000+CONSTS(91)/CONSTS(93)+STATES(4)/CONSTS(94)+STATES(6)/CONSTS(95))
      ALGBRC(191) = ( CONSTS(80)*ALGBRC(187)*CONSTS(91))/(1.00000+CONSTS(89)/CONSTS(90))
      ALGBRC(185) =  CONSTS(83)*EXP(( CONSTS(85)*ALGBRC(29))/3.00000)
      ALGBRC(188) = ( CONSTS(75)*STATES(4)/ALGBRC(185) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(185) ** 3.00000+1.00000+STATES(6)/CONSTS(86) ** 2.00000) - 1.00000)
      ALGBRC(189) = ( CONSTS(78)*CONSTS(2)/ALGBRC(186) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(186) ** 3.00000+1.00000+CONSTS(4)/CONSTS(87) ** 2.00000) - 1.00000)
      ALGBRC(192) = ( CONSTS(82)*STATES(6)/CONSTS(86) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(185) ** 3.00000+1.00000+STATES(6)/CONSTS(86) ** 2.00000) - 1.00000)
      ALGBRC(193) =  CONSTS(162)*ALGBRC(188)*CONSTS(161)+ ALGBRC(189)*ALGBRC(192)*ALGBRC(191)+ CONSTS(161)*ALGBRC(192)*ALGBRC(191)+ ALGBRC(191)*ALGBRC(188)*CONSTS(161)
      ALGBRC(194) =  ALGBRC(189)*CONSTS(160)*ALGBRC(192)+ ALGBRC(188)*CONSTS(161)*ALGBRC(190)+ ALGBRC(190)*CONSTS(160)*ALGBRC(192)+ CONSTS(161)*ALGBRC(190)*ALGBRC(192)
      ALGBRC(195) =  CONSTS(161)*ALGBRC(190)*CONSTS(162)+ ALGBRC(191)*ALGBRC(189)*CONSTS(160)+ ALGBRC(189)*CONSTS(160)*CONSTS(162)+ ALGBRC(190)*CONSTS(162)*CONSTS(160)
      ALGBRC(196) =  ALGBRC(192)*ALGBRC(191)*ALGBRC(189)+ ALGBRC(190)*CONSTS(162)*ALGBRC(188)+ ALGBRC(189)*CONSTS(162)*ALGBRC(188)+ ALGBRC(191)*ALGBRC(189)*ALGBRC(188)
      ALGBRC(197) = ALGBRC(193)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(198) = ALGBRC(194)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(201) =  3.00000*( ALGBRC(197)*ALGBRC(190) -  ALGBRC(198)*ALGBRC(191))
      ALGBRC(199) = ALGBRC(195)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(200) = ALGBRC(196)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(202) =  2.00000*( ALGBRC(200)*CONSTS(160) -  ALGBRC(199)*ALGBRC(188))
      ALGBRC(203) =  CONSTS(163)*( CONSTS(9)*ALGBRC(201)+ CONSTS(11)*ALGBRC(202))
      ALGBRC(204) = 1.00000/(1.00000+EXP(- (STATES(1) - 10.8968)/23.9871))
      ALGBRC(205) =  CONSTS(127)*ALGBRC(204)*(STATES(1) - ALGBRC(65))
      ALGBRC(67) =  CONSTS(38)*CONSTS(37)*CONSTS(116)*CONSTS(117)*(STATES(1) - ALGBRC(65))
      ALGBRC(12) = TERNRY(VOI.GE.CONSTS(15).AND.VOI.LE.CONSTS(16).AND.(VOI - CONSTS(15)) -  INT((VOI - CONSTS(15))/CONSTS(18))*CONSTS(18).LE.CONSTS(19), CONSTS(17), 0.00000)
      ALGBRC(99) = ( 0.500000*(STATES(4)+STATES(6)+CONSTS(35)+ 4.00000*STATES(10)))/1000.00
      ALGBRC(106) = EXP( - CONSTS(135)*1.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(109) = ( 1.00000*ALGBRC(26)*( ALGBRC(106)*STATES(6)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(142)*CONSTS(4)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(112) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(134)*ALGBRC(109)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(139)*ALGBRC(109)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(207) = (STATES(7) - STATES(6))/CONSTS(107)
      RATES(6) = ( - (((ALGBRC(77)+ALGBRC(116)+ALGBRC(118)+ALGBRC(122)+ALGBRC(205)+ALGBRC(67)+ALGBRC(12)) -  2.00000*ALGBRC(203))+ALGBRC(112))*CONSTS(137))/( CONSTS(8)*CONSTS(143))+( ALGBRC(207)*CONSTS(146))/CONSTS(143)
      ALGBRC(87) = EXP( - CONSTS(135)*1.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(90) = ( 1.00000*ALGBRC(26)*( ALGBRC(87)*STATES(7)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(142)*CONSTS(4)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(96) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(134)*ALGBRC(90)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(139)*ALGBRC(90)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      RATES(7) = ( - ALGBRC(96)*CONSTS(137))/( CONSTS(8)*CONSTS(146)) - ALGBRC(207)
      ALGBRC(64) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(9)*CONSTS(8)))*log(CONSTS(2)/STATES(4))
      ALGBRC(68) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(69) =  CONSTS(42)*(STATES(1) - ALGBRC(64))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(68))*STATES(12)*STATES(13)+ ALGBRC(68)*STATES(14)*STATES(15))
      ALGBRC(70) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(71) =  CONSTS(119)*(STATES(1) - ALGBRC(64))*STATES(16)*( (1.00000 - ALGBRC(70))*STATES(17)+ ALGBRC(70)*STATES(18))
      ALGBRC(151) = 1.00000/(1.00000+CONSTS(73)/STATES(10) ** 2.00000)
      ALGBRC(124) = EXP( CONSTS(71)*ALGBRC(29))
      ALGBRC(131) = 1.00000+ (CONSTS(2)/CONSTS(64))*(1.00000+1.00000/ALGBRC(124))
      ALGBRC(132) = CONSTS(2)/( CONSTS(64)*ALGBRC(124)*ALGBRC(131))
      ALGBRC(135) =  ALGBRC(132)*CONSTS(68)
      ALGBRC(125) = 1.00000+ (STATES(4)/CONSTS(64))*(1.00000+ALGBRC(124))
      ALGBRC(126) = ( STATES(4)*ALGBRC(124))/( CONSTS(64)*ALGBRC(125))
      ALGBRC(138) =  ALGBRC(126)*CONSTS(68)
      ALGBRC(128) = 1.00000+ (STATES(4)/CONSTS(62))*(1.00000+STATES(4)/CONSTS(63))
      ALGBRC(129) = ( STATES(4)*STATES(4))/( ALGBRC(128)*CONSTS(62)*CONSTS(63))
      ALGBRC(141) =  ALGBRC(129)*ALGBRC(126)*CONSTS(66)
      ALGBRC(142) =  ALGBRC(132)*CONSTS(148)*CONSTS(66)
      ALGBRC(133) = 1.00000/ALGBRC(131)
      ALGBRC(134) =  ALGBRC(133)*CONSTS(67)
      ALGBRC(136) = ALGBRC(134)+ALGBRC(135)
      ALGBRC(123) = EXP( CONSTS(72)*ALGBRC(29))
      ALGBRC(127) = 1.00000/ALGBRC(125)
      ALGBRC(137) = ( ALGBRC(127)*CONSTS(67))/ALGBRC(123)
      ALGBRC(139) = ALGBRC(137)+ALGBRC(138)
      ALGBRC(130) = 1.00000/ALGBRC(128)
      ALGBRC(140) =  ALGBRC(130)*STATES(10)*CONSTS(69)
      ALGBRC(143) =  CONSTS(151)*ALGBRC(139)*(ALGBRC(141)+ALGBRC(140))+ CONSTS(152)*ALGBRC(141)*(CONSTS(151)+ALGBRC(136))
      ALGBRC(144) =  CONSTS(150)*ALGBRC(141)*(ALGBRC(139)+CONSTS(152))+ ALGBRC(139)*ALGBRC(140)*(CONSTS(150)+ALGBRC(142))
      ALGBRC(145) =  CONSTS(150)*ALGBRC(136)*(ALGBRC(141)+ALGBRC(140))+ ALGBRC(142)*ALGBRC(140)*(CONSTS(151)+ALGBRC(136))
      ALGBRC(146) =  CONSTS(151)*ALGBRC(142)*(ALGBRC(139)+CONSTS(152))+ ALGBRC(136)*CONSTS(152)*(CONSTS(150)+ALGBRC(142))
      ALGBRC(147) = ALGBRC(143)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(148) = ALGBRC(144)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(149) = ALGBRC(145)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(150) = ALGBRC(146)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(152) = ( 3.00000*( ALGBRC(150)*ALGBRC(141) -  ALGBRC(147)*ALGBRC(142))+ ALGBRC(149)*ALGBRC(138)) -  ALGBRC(148)*ALGBRC(135)
      ALGBRC(153) =  ALGBRC(148)*CONSTS(151) -  ALGBRC(147)*CONSTS(150)
      ALGBRC(154) =  (1.00000 - CONSTS(61))*CONSTS(153)*ALGBRC(151)*( CONSTS(9)*ALGBRC(152)+ CONSTS(10)*ALGBRC(153))
      ALGBRC(206) = ( CONSTS(98)*ALGBRC(26)*( STATES(4)*EXP(ALGBRC(29)) - CONSTS(2)))/(EXP(ALGBRC(29)) - 1.00000)
      ALGBRC(105) = EXP( - CONSTS(135)*1.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(108) = ( 1.00000*ALGBRC(26)*( ALGBRC(105)*STATES(4)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(141)*CONSTS(2)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(111) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(133)*ALGBRC(108)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(138)*ALGBRC(108)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(209) = (STATES(5) - STATES(4))/CONSTS(106)
      RATES(4) = ( - (ALGBRC(69)+ALGBRC(71)+ 3.00000*ALGBRC(154)+ALGBRC(111)+ 3.00000*ALGBRC(203)+ALGBRC(206))*CONSTS(137))/( CONSTS(8)*CONSTS(143))+( ALGBRC(209)*CONSTS(146))/CONSTS(143)
      ALGBRC(181) = 1.00000/(1.00000+CONSTS(73)/STATES(3) ** 2.00000)
      ALGBRC(161) = 1.00000+ (CONSTS(2)/CONSTS(64))*(1.00000+1.00000/ALGBRC(124))
      ALGBRC(162) = CONSTS(2)/( CONSTS(64)*ALGBRC(124)*ALGBRC(161))
      ALGBRC(165) =  ALGBRC(162)*CONSTS(68)
      ALGBRC(155) = 1.00000+ (STATES(5)/CONSTS(64))*(1.00000+ALGBRC(124))
      ALGBRC(156) = ( STATES(5)*ALGBRC(124))/( CONSTS(64)*ALGBRC(155))
      ALGBRC(168) =  ALGBRC(156)*CONSTS(68)
      ALGBRC(158) = 1.00000+ (STATES(5)/CONSTS(62))*(1.00000+STATES(5)/CONSTS(63))
      ALGBRC(159) = ( STATES(5)*STATES(5))/( ALGBRC(158)*CONSTS(62)*CONSTS(63))
      ALGBRC(171) =  ALGBRC(159)*ALGBRC(156)*CONSTS(66)
      ALGBRC(172) =  ALGBRC(162)*CONSTS(155)*CONSTS(66)
      ALGBRC(163) = 1.00000/ALGBRC(161)
      ALGBRC(164) =  ALGBRC(163)*CONSTS(67)
      ALGBRC(166) = ALGBRC(164)+ALGBRC(165)
      ALGBRC(157) = 1.00000/ALGBRC(155)
      ALGBRC(167) = ( ALGBRC(157)*CONSTS(67))/ALGBRC(123)
      ALGBRC(169) = ALGBRC(167)+ALGBRC(168)
      ALGBRC(160) = 1.00000/ALGBRC(158)
      ALGBRC(170) =  ALGBRC(160)*STATES(3)*CONSTS(69)
      ALGBRC(173) =  CONSTS(158)*ALGBRC(169)*(ALGBRC(171)+ALGBRC(170))+ CONSTS(159)*ALGBRC(171)*(CONSTS(158)+ALGBRC(166))
      ALGBRC(174) =  CONSTS(157)*ALGBRC(171)*(ALGBRC(169)+CONSTS(159))+ ALGBRC(169)*ALGBRC(170)*(CONSTS(157)+ALGBRC(172))
      ALGBRC(175) =  CONSTS(157)*ALGBRC(166)*(ALGBRC(171)+ALGBRC(170))+ ALGBRC(172)*ALGBRC(170)*(CONSTS(158)+ALGBRC(166))
      ALGBRC(176) =  CONSTS(158)*ALGBRC(172)*(ALGBRC(169)+CONSTS(159))+ ALGBRC(166)*CONSTS(159)*(CONSTS(157)+ALGBRC(172))
      ALGBRC(177) = ALGBRC(173)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(178) = ALGBRC(174)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(179) = ALGBRC(175)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(180) = ALGBRC(176)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(182) = ( 3.00000*( ALGBRC(180)*ALGBRC(171) -  ALGBRC(177)*ALGBRC(172))+ ALGBRC(179)*ALGBRC(168)) -  ALGBRC(178)*ALGBRC(165)
      ALGBRC(183) =  ALGBRC(178)*CONSTS(158) -  ALGBRC(177)*CONSTS(157)
      ALGBRC(184) =  CONSTS(61)*CONSTS(153)*ALGBRC(181)*( CONSTS(9)*ALGBRC(182)+ CONSTS(10)*ALGBRC(183))
      ALGBRC(86) = EXP( - CONSTS(135)*1.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(89) = ( 1.00000*ALGBRC(26)*( ALGBRC(86)*STATES(5)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(141)*CONSTS(2)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(93) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(133)*ALGBRC(89)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(138)*ALGBRC(89)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      RATES(5) = ( - (ALGBRC(93)+ 3.00000*ALGBRC(184))*CONSTS(137))/( CONSTS(8)*CONSTS(146)) - ALGBRC(209)
      ALGBRC(212) = (STATES(3) - STATES(10))/CONSTS(108)
      ALGBRC(214) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(216) =  CONSTS(111)*( (1.00000 - ALGBRC(214))*STATES(42)+ ALGBRC(214)*STATES(43))
      ALGBRC(57) = 1.00000/(1.00000+( CONSTS(29)*CONSTS(30))/CONSTS(30)+STATES(3) ** 2.00000+( CONSTS(31)*CONSTS(32))/CONSTS(32)+STATES(3) ** 2.00000)
      RATES(3) =  ALGBRC(57)*((( - (ALGBRC(92) -  2.00000*ALGBRC(184))*CONSTS(137))/( 2.00000*CONSTS(8)*CONSTS(146))+( ALGBRC(216)*CONSTS(145))/CONSTS(146)) - ALGBRC(212))
      ALGBRC(102) = EXP( - CONSTS(135)*4.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(107) = ( 4.00000*ALGBRC(26)*( ALGBRC(102)*STATES(10)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(110) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(122)*ALGBRC(107)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(132)*ALGBRC(107)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(113) = ALGBRC(92)+ALGBRC(110)
      ALGBRC(114) = ALGBRC(93)+ALGBRC(111)
      ALGBRC(115) = ALGBRC(96)+ALGBRC(112)
      ALGBRC(210) = ( CONSTS(100)*STATES(10))/(CONSTS(101)+STATES(10))
      ALGBRC(208) = ( CONSTS(99)*4.00000*ALGBRC(26)*( ALGBRC(102)*STATES(10)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(211) =  (( CONSTS(105)*CONSTS(102))/(1.00000+CONSTS(104)/STATES(3)))*(STATES(1) - CONSTS(115))
      ALGBRC(213) =  (( (1.00000 - CONSTS(105))*CONSTS(102))/(1.00000+CONSTS(104)/STATES(10)))*(STATES(1) - CONSTS(115))
      ALGBRC(215) = ALGBRC(211)+ALGBRC(213)
      ALGBRC(217) =  CONSTS(103)*(STATES(1) - CONSTS(115))
      RATES(1) = - (ALGBRC(69)+ALGBRC(71)+ALGBRC(77)+ALGBRC(113)+ALGBRC(114)+ALGBRC(115)+ALGBRC(116)+ALGBRC(118)+ALGBRC(122)+ALGBRC(154)+ALGBRC(184)+ALGBRC(203)+ALGBRC(206)+ALGBRC(205)+ALGBRC(210)+ALGBRC(208)+ALGBRC(215)+ALGBRC(217)+ALGBRC(67)+ALGBRC(12))
      ALGBRC(218) = ( CONSTS(130)*0.00542500*STATES(10))/(STATES(10)+0.000920000)
      ALGBRC(219) = ( CONSTS(130)*2.75000*0.00542500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
      ALGBRC(220) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(221) = ( 0.00488250*STATES(8))/15.0000
      ALGBRC(222) =  CONSTS(112)*(( (1.00000 - ALGBRC(220))*ALGBRC(218)+ ALGBRC(220)*ALGBRC(219)) - ALGBRC(221))
      ALGBRC(54) = 1.00000/(1.00000+( CONSTS(114)*CONSTS(26))/CONSTS(26)+STATES(10) ** 2.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(10) ** 2.00000)
      RATES(10) =  ALGBRC(54)*((( - ((ALGBRC(110)+ALGBRC(210)+ALGBRC(208)) -  2.00000*ALGBRC(154))*CONSTS(137))/( 2.00000*CONSTS(8)*CONSTS(143)) - ( ALGBRC(222)*CONSTS(144))/CONSTS(143))+( ALGBRC(212)*CONSTS(146))/CONSTS(143))
      ALGBRC(223) = (STATES(8) - STATES(9))/60.0000
      RATES(8) = ALGBRC(222) - ( ALGBRC(223)*CONSTS(145))/CONSTS(144)
      ALGBRC(60) = 1.00000/(1.00000+( CONSTS(33)*CONSTS(34))/CONSTS(34)+STATES(9) ** 2.00000)
      RATES(9) =  ALGBRC(60)*(ALGBRC(223) - ALGBRC(216))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+18.0800)/2.79160))
      ALGBRC(1) = 1.00000/1.00000+EXP(- (STATES(1)+56.8600)/9.03000) ** 2.00000
      ALGBRC(14) =  0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(17) =  0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 14.3400)/14.8200))
      ALGBRC(18) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+CONSTS(46)+100.000)/29.3814)))
      ALGBRC(8) = TERNRY(STATES(1).GE.31.4978, 1.00000,  1.07630*EXP( - 1.00700*EXP( - 0.0829000*STATES(1))))
      ALGBRC(23) = CONSTS(53)+0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+CONSTS(52)+6.00000))+EXP( 0.0900000*(STATES(1)+CONSTS(52)+14.0000)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
      ALGBRC(24) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000)/10.0000))
      ALGBRC(25) = 1000.00+1.00000/( 3.50000e-05*EXP(- (STATES(1)+5.00000)/4.00000)+ 3.50000e-05*EXP((STATES(1)+5.00000)/6.00000))
      ALGBRC(11) =  STATES(30)*1.00000
      ALGBRC(21) = 1.00000/(CONSTS(48)/ALGBRC(11)+1.00000+CONSTS(47)/STATES(3) ** 4.00000)
      ALGBRC(22) = 1.00000/(CONSTS(48)/ALGBRC(11)+1.00000+CONSTS(47)/STATES(10) ** 4.00000)
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(28) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      ALGBRC(33) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(46)) - 24.3400)/14.8200))
      ALGBRC(20) = ALGBRC(9)
      ALGBRC(34) = 7.00000+1.00000/( 0.0400000*EXP(- (STATES(1) - 4.00000)/7.00000)+ 0.0400000*EXP((STATES(1) - 4.00000)/7.00000))
      ALGBRC(35) = 100.000+1.00000/( 0.000120000*EXP(- STATES(1)/3.00000)+ 0.000120000*EXP(STATES(1)/7.00000))
      ALGBRC(36) =  2.50000*ALGBRC(24)
      ALGBRC(27) = ALGBRC(13)
      ALGBRC(37) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      ALGBRC(44) = ( CONSTS(23)*(1.00000 - STATES(2)))/(1.00000+CONSTS(24)/STATES(3))
      ALGBRC(2) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
      ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, 0.00000,  0.0570000*EXP(- (STATES(1)+80.0000)/6.80000))
      ALGBRC(30) = TERNRY(STATES(1).GE.- 40.0000, 0.770000/( 0.130000*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000))),  2.70000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.348500*STATES(1)))
      ALGBRC(38) = 1.00000/(ALGBRC(15)+ALGBRC(30))
      ALGBRC(41) =  2.50000*ALGBRC(34)
      ALGBRC(39) = ALGBRC(2)
      ALGBRC(16) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 25428.0*EXP( 0.244400*STATES(1)) -  6.94800e-06*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
      ALGBRC(31) = TERNRY(STATES(1).GE.- 40.0000, ( 0.600000*EXP( 0.0570000*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.0242400*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
      ALGBRC(45) = 1.00000/(ALGBRC(16)+ALGBRC(31))
      ALGBRC(46) = 1.00000/1.00000+EXP((STATES(1)+77.5500)/7.43000) ** 2.00000
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(46)+43.9400)/5.71100))
      ALGBRC(19) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+CONSTS(46)+70.0000)/5.00000)), 1.00000)
      ALGBRC(32) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+CONSTS(46)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+CONSTS(46)+50.0000)/16.5900))
      ALGBRC(47) =  ALGBRC(32)*ALGBRC(19)
      ALGBRC(29) = ( STATES(1)*CONSTS(8))/( CONSTS(6)*CONSTS(7))
      ALGBRC(42) =  0.116100*EXP( 0.299000*ALGBRC(29))
      ALGBRC(48) =  0.244200*EXP( - 1.60400*ALGBRC(29))
      ALGBRC(51) =  1.46000*ALGBRC(45)
      ALGBRC(40) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+CONSTS(46)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+CONSTS(46)+114.100)/8.07900))
      ALGBRC(52) =  ALGBRC(40)*ALGBRC(19)
      ALGBRC(43) =  0.0578000*EXP( 0.971000*ALGBRC(29))
      ALGBRC(49) =  0.000349000*EXP( - 1.06200*ALGBRC(29))
      ALGBRC(53) =  0.253300*EXP( 0.595300*ALGBRC(29))
      ALGBRC(56) =  0.0652500*EXP( - 0.820900*ALGBRC(29))
      ALGBRC(55) = 1.35400+0.000100000/(EXP(((STATES(1)+CONSTS(46)) - 167.400)/15.8900)+EXP(- ((STATES(1)+CONSTS(46)) - 12.2300)/0.215400))
      ALGBRC(58) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+CONSTS(46)+70.0000)/20.0000))
      ALGBRC(61) =  ALGBRC(55)*ALGBRC(58)*ALGBRC(47)
      ALGBRC(62) =  ALGBRC(55)*ALGBRC(58)*ALGBRC(52)
      ALGBRC(59) =  5.20000e-05*EXP( 1.52500*ALGBRC(29))
      ALGBRC(63) = ( ALGBRC(49)*ALGBRC(56)*ALGBRC(59))/( ALGBRC(43)*ALGBRC(53))
      ALGBRC(78) =  CONSTS(50)*STATES(26)+ CONSTS(121)*STATES(27)
      ALGBRC(79) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(80) = 1.00000 - ALGBRC(79)
      ALGBRC(81) =  ALGBRC(79)*STATES(28)+ ALGBRC(80)*STATES(29)
      ALGBRC(82) =  CONSTS(50)*STATES(31)+ CONSTS(121)*STATES(27)
      ALGBRC(83) =  ALGBRC(79)*STATES(32)+ ALGBRC(80)*STATES(29)
      ALGBRC(26) = ( STATES(1)*CONSTS(8)*CONSTS(8))/( CONSTS(6)*CONSTS(7))
      ALGBRC(84) = ( 0.500000*(STATES(5)+STATES(7)+CONSTS(35)+ 4.00000*STATES(3)))/1000.00
      ALGBRC(85) = EXP( - CONSTS(135)*4.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(88) = ( 4.00000*ALGBRC(26)*( ALGBRC(85)*STATES(3)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(50) = ALGBRC(44)+STATES(2)
      ALGBRC(91) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(92) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(122)*ALGBRC(88)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(132)*ALGBRC(88)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      ALGBRC(94) = (( - CONSTS(128)*ALGBRC(92))/1.00000)/(1.00000+CONSTS(110)/STATES(9) ** 8.00000)
      ALGBRC(97) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(94)*1.70000, ALGBRC(94))
      ALGBRC(100) = CONSTS(109)/(1.00000+0.0123000/STATES(9))
      ALGBRC(103) = TERNRY(ALGBRC(100).LT.0.00100000, 0.00100000, ALGBRC(100))
      ALGBRC(95) = (( - CONSTS(136)*ALGBRC(92))/1.00000)/(1.00000+CONSTS(110)/STATES(9) ** 8.00000)
      ALGBRC(98) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(95)*1.70000, ALGBRC(95))
      ALGBRC(101) = CONSTS(129)/(1.00000+0.0123000/STATES(9))
      ALGBRC(104) = TERNRY(ALGBRC(101).LT.0.00100000, 0.00100000, ALGBRC(101))
      ALGBRC(65) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(11)*CONSTS(8)))*log(CONSTS(4)/STATES(6))
      ALGBRC(72) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(46)) - 213.600)/151.200))
      ALGBRC(73) = 1.00000 - ALGBRC(72)
      ALGBRC(74) =  ALGBRC(72)*STATES(20)+ ALGBRC(73)*STATES(21)
      ALGBRC(75) =  ALGBRC(72)*STATES(23)+ ALGBRC(73)*STATES(24)
      ALGBRC(76) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(77) =  CONSTS(120)*(STATES(1) - ALGBRC(65))*( (1.00000 - ALGBRC(76))*STATES(19)*ALGBRC(74)+ ALGBRC(76)*STATES(22)*ALGBRC(75))
      ALGBRC(116) =  CONSTS(124)* (CONSTS(4)/5.00000) ** (1.0 / 2)*STATES(39)*(STATES(1) - ALGBRC(65))
      ALGBRC(66) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(11)*CONSTS(8)))*log((CONSTS(4)+ CONSTS(36)*CONSTS(2))/(STATES(6)+ CONSTS(36)*STATES(4)))
      ALGBRC(117) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
      ALGBRC(118) =  CONSTS(125)*ALGBRC(117)*STATES(40)*STATES(41)*(STATES(1) - ALGBRC(66))
      ALGBRC(119) = 4.09400/(1.00000+EXP( 0.121700*((STATES(1) - ALGBRC(65)) - 49.9340)))
      ALGBRC(120) = ( 15.7200*EXP( 0.0674000*((STATES(1) - ALGBRC(65)) - 3.25700))+EXP( 0.0618000*((STATES(1) - ALGBRC(65)) - 594.310)))/(1.00000+EXP( - 0.162900*((STATES(1) - ALGBRC(65))+14.2070)))
      ALGBRC(121) = ALGBRC(119)/(ALGBRC(119)+ALGBRC(120))
      ALGBRC(122) =  CONSTS(126)* (CONSTS(4)/5.00000) ** (1.0 / 2)*ALGBRC(121)*(STATES(1) - ALGBRC(65))
      ALGBRC(186) =  CONSTS(84)*EXP(( (1.00000 - CONSTS(85))*ALGBRC(29))/3.00000)
      ALGBRC(190) = ( CONSTS(79)*CONSTS(4)/CONSTS(87) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(186) ** 3.00000+1.00000+CONSTS(4)/CONSTS(87) ** 2.00000) - 1.00000)
      ALGBRC(187) = CONSTS(92)/(1.00000+CONSTS(91)/CONSTS(93)+STATES(4)/CONSTS(94)+STATES(6)/CONSTS(95))
      ALGBRC(191) = ( CONSTS(80)*ALGBRC(187)*CONSTS(91))/(1.00000+CONSTS(89)/CONSTS(90))
      ALGBRC(185) =  CONSTS(83)*EXP(( CONSTS(85)*ALGBRC(29))/3.00000)
      ALGBRC(188) = ( CONSTS(75)*STATES(4)/ALGBRC(185) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(185) ** 3.00000+1.00000+STATES(6)/CONSTS(86) ** 2.00000) - 1.00000)
      ALGBRC(189) = ( CONSTS(78)*CONSTS(2)/ALGBRC(186) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(186) ** 3.00000+1.00000+CONSTS(4)/CONSTS(87) ** 2.00000) - 1.00000)
      ALGBRC(192) = ( CONSTS(82)*STATES(6)/CONSTS(86) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(185) ** 3.00000+1.00000+STATES(6)/CONSTS(86) ** 2.00000) - 1.00000)
      ALGBRC(193) =  CONSTS(162)*ALGBRC(188)*CONSTS(161)+ ALGBRC(189)*ALGBRC(192)*ALGBRC(191)+ CONSTS(161)*ALGBRC(192)*ALGBRC(191)+ ALGBRC(191)*ALGBRC(188)*CONSTS(161)
      ALGBRC(194) =  ALGBRC(189)*CONSTS(160)*ALGBRC(192)+ ALGBRC(188)*CONSTS(161)*ALGBRC(190)+ ALGBRC(190)*CONSTS(160)*ALGBRC(192)+ CONSTS(161)*ALGBRC(190)*ALGBRC(192)
      ALGBRC(195) =  CONSTS(161)*ALGBRC(190)*CONSTS(162)+ ALGBRC(191)*ALGBRC(189)*CONSTS(160)+ ALGBRC(189)*CONSTS(160)*CONSTS(162)+ ALGBRC(190)*CONSTS(162)*CONSTS(160)
      ALGBRC(196) =  ALGBRC(192)*ALGBRC(191)*ALGBRC(189)+ ALGBRC(190)*CONSTS(162)*ALGBRC(188)+ ALGBRC(189)*CONSTS(162)*ALGBRC(188)+ ALGBRC(191)*ALGBRC(189)*ALGBRC(188)
      ALGBRC(197) = ALGBRC(193)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(198) = ALGBRC(194)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(201) =  3.00000*( ALGBRC(197)*ALGBRC(190) -  ALGBRC(198)*ALGBRC(191))
      ALGBRC(199) = ALGBRC(195)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(200) = ALGBRC(196)/(ALGBRC(193)+ALGBRC(194)+ALGBRC(195)+ALGBRC(196))
      ALGBRC(202) =  2.00000*( ALGBRC(200)*CONSTS(160) -  ALGBRC(199)*ALGBRC(188))
      ALGBRC(203) =  CONSTS(163)*( CONSTS(9)*ALGBRC(201)+ CONSTS(11)*ALGBRC(202))
      ALGBRC(204) = 1.00000/(1.00000+EXP(- (STATES(1) - 10.8968)/23.9871))
      ALGBRC(205) =  CONSTS(127)*ALGBRC(204)*(STATES(1) - ALGBRC(65))
      ALGBRC(67) =  CONSTS(38)*CONSTS(37)*CONSTS(116)*CONSTS(117)*(STATES(1) - ALGBRC(65))
      ALGBRC(12) = TERNRY(VOI.GE.CONSTS(15).AND.VOI.LE.CONSTS(16).AND.(VOI - CONSTS(15)) -  INT((VOI - CONSTS(15))/CONSTS(18))*CONSTS(18).LE.CONSTS(19), CONSTS(17), 0.00000)
      ALGBRC(99) = ( 0.500000*(STATES(4)+STATES(6)+CONSTS(35)+ 4.00000*STATES(10)))/1000.00
      ALGBRC(106) = EXP( - CONSTS(135)*1.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(109) = ( 1.00000*ALGBRC(26)*( ALGBRC(106)*STATES(6)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(142)*CONSTS(4)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(112) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(134)*ALGBRC(109)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(139)*ALGBRC(109)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(207) = (STATES(7) - STATES(6))/CONSTS(107)
      ALGBRC(87) = EXP( - CONSTS(135)*1.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(90) = ( 1.00000*ALGBRC(26)*( ALGBRC(87)*STATES(7)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(142)*CONSTS(4)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(96) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(134)*ALGBRC(90)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(139)*ALGBRC(90)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      ALGBRC(64) =  (( CONSTS(6)*CONSTS(7))/( CONSTS(9)*CONSTS(8)))*log(CONSTS(2)/STATES(4))
      ALGBRC(68) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(69) =  CONSTS(42)*(STATES(1) - ALGBRC(64))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(68))*STATES(12)*STATES(13)+ ALGBRC(68)*STATES(14)*STATES(15))
      ALGBRC(70) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(71) =  CONSTS(119)*(STATES(1) - ALGBRC(64))*STATES(16)*( (1.00000 - ALGBRC(70))*STATES(17)+ ALGBRC(70)*STATES(18))
      ALGBRC(151) = 1.00000/(1.00000+CONSTS(73)/STATES(10) ** 2.00000)
      ALGBRC(124) = EXP( CONSTS(71)*ALGBRC(29))
      ALGBRC(131) = 1.00000+ (CONSTS(2)/CONSTS(64))*(1.00000+1.00000/ALGBRC(124))
      ALGBRC(132) = CONSTS(2)/( CONSTS(64)*ALGBRC(124)*ALGBRC(131))
      ALGBRC(135) =  ALGBRC(132)*CONSTS(68)
      ALGBRC(125) = 1.00000+ (STATES(4)/CONSTS(64))*(1.00000+ALGBRC(124))
      ALGBRC(126) = ( STATES(4)*ALGBRC(124))/( CONSTS(64)*ALGBRC(125))
      ALGBRC(138) =  ALGBRC(126)*CONSTS(68)
      ALGBRC(128) = 1.00000+ (STATES(4)/CONSTS(62))*(1.00000+STATES(4)/CONSTS(63))
      ALGBRC(129) = ( STATES(4)*STATES(4))/( ALGBRC(128)*CONSTS(62)*CONSTS(63))
      ALGBRC(141) =  ALGBRC(129)*ALGBRC(126)*CONSTS(66)
      ALGBRC(142) =  ALGBRC(132)*CONSTS(148)*CONSTS(66)
      ALGBRC(133) = 1.00000/ALGBRC(131)
      ALGBRC(134) =  ALGBRC(133)*CONSTS(67)
      ALGBRC(136) = ALGBRC(134)+ALGBRC(135)
      ALGBRC(123) = EXP( CONSTS(72)*ALGBRC(29))
      ALGBRC(127) = 1.00000/ALGBRC(125)
      ALGBRC(137) = ( ALGBRC(127)*CONSTS(67))/ALGBRC(123)
      ALGBRC(139) = ALGBRC(137)+ALGBRC(138)
      ALGBRC(130) = 1.00000/ALGBRC(128)
      ALGBRC(140) =  ALGBRC(130)*STATES(10)*CONSTS(69)
      ALGBRC(143) =  CONSTS(151)*ALGBRC(139)*(ALGBRC(141)+ALGBRC(140))+ CONSTS(152)*ALGBRC(141)*(CONSTS(151)+ALGBRC(136))
      ALGBRC(144) =  CONSTS(150)*ALGBRC(141)*(ALGBRC(139)+CONSTS(152))+ ALGBRC(139)*ALGBRC(140)*(CONSTS(150)+ALGBRC(142))
      ALGBRC(145) =  CONSTS(150)*ALGBRC(136)*(ALGBRC(141)+ALGBRC(140))+ ALGBRC(142)*ALGBRC(140)*(CONSTS(151)+ALGBRC(136))
      ALGBRC(146) =  CONSTS(151)*ALGBRC(142)*(ALGBRC(139)+CONSTS(152))+ ALGBRC(136)*CONSTS(152)*(CONSTS(150)+ALGBRC(142))
      ALGBRC(147) = ALGBRC(143)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(148) = ALGBRC(144)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(149) = ALGBRC(145)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(150) = ALGBRC(146)/(ALGBRC(143)+ALGBRC(144)+ALGBRC(145)+ALGBRC(146))
      ALGBRC(152) = ( 3.00000*( ALGBRC(150)*ALGBRC(141) -  ALGBRC(147)*ALGBRC(142))+ ALGBRC(149)*ALGBRC(138)) -  ALGBRC(148)*ALGBRC(135)
      ALGBRC(153) =  ALGBRC(148)*CONSTS(151) -  ALGBRC(147)*CONSTS(150)
      ALGBRC(154) =  (1.00000 - CONSTS(61))*CONSTS(153)*ALGBRC(151)*( CONSTS(9)*ALGBRC(152)+ CONSTS(10)*ALGBRC(153))
      ALGBRC(206) = ( CONSTS(98)*ALGBRC(26)*( STATES(4)*EXP(ALGBRC(29)) - CONSTS(2)))/(EXP(ALGBRC(29)) - 1.00000)
      ALGBRC(105) = EXP( - CONSTS(135)*1.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(108) = ( 1.00000*ALGBRC(26)*( ALGBRC(105)*STATES(4)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(141)*CONSTS(2)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(111) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(133)*ALGBRC(108)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(138)*ALGBRC(108)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(209) = (STATES(5) - STATES(4))/CONSTS(106)
      ALGBRC(181) = 1.00000/(1.00000+CONSTS(73)/STATES(3) ** 2.00000)
      ALGBRC(161) = 1.00000+ (CONSTS(2)/CONSTS(64))*(1.00000+1.00000/ALGBRC(124))
      ALGBRC(162) = CONSTS(2)/( CONSTS(64)*ALGBRC(124)*ALGBRC(161))
      ALGBRC(165) =  ALGBRC(162)*CONSTS(68)
      ALGBRC(155) = 1.00000+ (STATES(5)/CONSTS(64))*(1.00000+ALGBRC(124))
      ALGBRC(156) = ( STATES(5)*ALGBRC(124))/( CONSTS(64)*ALGBRC(155))
      ALGBRC(168) =  ALGBRC(156)*CONSTS(68)
      ALGBRC(158) = 1.00000+ (STATES(5)/CONSTS(62))*(1.00000+STATES(5)/CONSTS(63))
      ALGBRC(159) = ( STATES(5)*STATES(5))/( ALGBRC(158)*CONSTS(62)*CONSTS(63))
      ALGBRC(171) =  ALGBRC(159)*ALGBRC(156)*CONSTS(66)
      ALGBRC(172) =  ALGBRC(162)*CONSTS(155)*CONSTS(66)
      ALGBRC(163) = 1.00000/ALGBRC(161)
      ALGBRC(164) =  ALGBRC(163)*CONSTS(67)
      ALGBRC(166) = ALGBRC(164)+ALGBRC(165)
      ALGBRC(157) = 1.00000/ALGBRC(155)
      ALGBRC(167) = ( ALGBRC(157)*CONSTS(67))/ALGBRC(123)
      ALGBRC(169) = ALGBRC(167)+ALGBRC(168)
      ALGBRC(160) = 1.00000/ALGBRC(158)
      ALGBRC(170) =  ALGBRC(160)*STATES(3)*CONSTS(69)
      ALGBRC(173) =  CONSTS(158)*ALGBRC(169)*(ALGBRC(171)+ALGBRC(170))+ CONSTS(159)*ALGBRC(171)*(CONSTS(158)+ALGBRC(166))
      ALGBRC(174) =  CONSTS(157)*ALGBRC(171)*(ALGBRC(169)+CONSTS(159))+ ALGBRC(169)*ALGBRC(170)*(CONSTS(157)+ALGBRC(172))
      ALGBRC(175) =  CONSTS(157)*ALGBRC(166)*(ALGBRC(171)+ALGBRC(170))+ ALGBRC(172)*ALGBRC(170)*(CONSTS(158)+ALGBRC(166))
      ALGBRC(176) =  CONSTS(158)*ALGBRC(172)*(ALGBRC(169)+CONSTS(159))+ ALGBRC(166)*CONSTS(159)*(CONSTS(157)+ALGBRC(172))
      ALGBRC(177) = ALGBRC(173)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(178) = ALGBRC(174)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(179) = ALGBRC(175)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(180) = ALGBRC(176)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
      ALGBRC(182) = ( 3.00000*( ALGBRC(180)*ALGBRC(171) -  ALGBRC(177)*ALGBRC(172))+ ALGBRC(179)*ALGBRC(168)) -  ALGBRC(178)*ALGBRC(165)
      ALGBRC(183) =  ALGBRC(178)*CONSTS(158) -  ALGBRC(177)*CONSTS(157)
      ALGBRC(184) =  CONSTS(61)*CONSTS(153)*ALGBRC(181)*( CONSTS(9)*ALGBRC(182)+ CONSTS(10)*ALGBRC(183))
      ALGBRC(86) = EXP( - CONSTS(135)*1.00000*( ALGBRC(84) ** (1.0 / 2)/(1.00000+ ALGBRC(84) ** (1.0 / 2)) -  0.300000*ALGBRC(84)))
      ALGBRC(89) = ( 1.00000*ALGBRC(26)*( ALGBRC(86)*STATES(5)*EXP( 1.00000*ALGBRC(29)) -  CONSTS(141)*CONSTS(2)))/(EXP( 1.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(93) =  CONSTS(55)*( (1.00000 - ALGBRC(91))*CONSTS(133)*ALGBRC(89)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(81)*STATES(33))+ ALGBRC(91)*CONSTS(138)*ALGBRC(89)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(33))+ STATES(30)*ALGBRC(83)*STATES(33)))
      ALGBRC(212) = (STATES(3) - STATES(10))/CONSTS(108)
      ALGBRC(214) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(216) =  CONSTS(111)*( (1.00000 - ALGBRC(214))*STATES(42)+ ALGBRC(214)*STATES(43))
      ALGBRC(57) = 1.00000/(1.00000+( CONSTS(29)*CONSTS(30))/CONSTS(30)+STATES(3) ** 2.00000+( CONSTS(31)*CONSTS(32))/CONSTS(32)+STATES(3) ** 2.00000)
      ALGBRC(102) = EXP( - CONSTS(135)*4.00000*( ALGBRC(99) ** (1.0 / 2)/(1.00000+ ALGBRC(99) ** (1.0 / 2)) -  0.300000*ALGBRC(99)))
      ALGBRC(107) = ( 4.00000*ALGBRC(26)*( ALGBRC(102)*STATES(10)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(110) =  (1.00000 - CONSTS(55))*( (1.00000 - ALGBRC(91))*CONSTS(122)*ALGBRC(107)*STATES(25)*( ALGBRC(78)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(81)*STATES(34))+ ALGBRC(91)*CONSTS(132)*ALGBRC(107)*STATES(25)*( ALGBRC(82)*(1.00000 - STATES(34))+ STATES(30)*ALGBRC(83)*STATES(34)))
      ALGBRC(113) = ALGBRC(92)+ALGBRC(110)
      ALGBRC(114) = ALGBRC(93)+ALGBRC(111)
      ALGBRC(115) = ALGBRC(96)+ALGBRC(112)
      ALGBRC(210) = ( CONSTS(100)*STATES(10))/(CONSTS(101)+STATES(10))
      ALGBRC(208) = ( CONSTS(99)*4.00000*ALGBRC(26)*( ALGBRC(102)*STATES(10)*EXP( 2.00000*ALGBRC(29)) -  CONSTS(140)*CONSTS(3)))/(EXP( 2.00000*ALGBRC(29)) - 1.00000)
      ALGBRC(211) =  (( CONSTS(105)*CONSTS(102))/(1.00000+CONSTS(104)/STATES(3)))*(STATES(1) - CONSTS(115))
      ALGBRC(213) =  (( (1.00000 - CONSTS(105))*CONSTS(102))/(1.00000+CONSTS(104)/STATES(10)))*(STATES(1) - CONSTS(115))
      ALGBRC(215) = ALGBRC(211)+ALGBRC(213)
      ALGBRC(217) =  CONSTS(103)*(STATES(1) - CONSTS(115))
      ALGBRC(218) = ( CONSTS(130)*0.00542500*STATES(10))/(STATES(10)+0.000920000)
      ALGBRC(219) = ( CONSTS(130)*2.75000*0.00542500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
      ALGBRC(220) = 1.00000/(1.00000+CONSTS(20)/ALGBRC(50))
      ALGBRC(221) = ( 0.00488250*STATES(8))/15.0000
      ALGBRC(222) =  CONSTS(112)*(( (1.00000 - ALGBRC(220))*ALGBRC(218)+ ALGBRC(220)*ALGBRC(219)) - ALGBRC(221))
      ALGBRC(54) = 1.00000/(1.00000+( CONSTS(114)*CONSTS(26))/CONSTS(26)+STATES(10) ** 2.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(10) ** 2.00000)
      ALGBRC(223) = (STATES(8) - STATES(9))/60.0000
      ALGBRC(60) = 1.00000/(1.00000+( CONSTS(33)*CONSTS(34))/CONSTS(34)+STATES(9) ** 2.00000)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END