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 199 entries in the algebraic variable array.
C There are a total of 46 entries in each of the rate and state variable arrays.
C There are a total of 146 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is nao in component extracellular (millimolar).
C CONSTS(2) is cao in component extracellular (millimolar).
C CONSTS(3) is ko in component extracellular (millimolar).
C CONSTS(4) is R in component physical_constants (joule_per_kilomole_kelvin).
C CONSTS(5) is T in component physical_constants (kelvin).
C CONSTS(6) is F in component physical_constants (coulomb_per_mole).
C CONSTS(7) is zna in component physical_constants (dimensionless).
C CONSTS(8) is zca in component physical_constants (dimensionless).
C CONSTS(9) is zk in component physical_constants (dimensionless).
C CONSTS(10) is L in component cell_geometry (centimeter).
C CONSTS(11) is rad in component cell_geometry (centimeter).
C CONSTS(12) is greekpi in component cell_geometry (dimensionless).
C CONSTS(120) is vcell in component cell_geometry (microliter).
C CONSTS(124) is Ageo in component cell_geometry (centimeter_squared).
C CONSTS(125) is Acap in component cell_geometry (centimeter_squared).
C CONSTS(126) is vmyo in component cell_geometry (microliter).
C CONSTS(127) is vnsr in component cell_geometry (microliter).
C CONSTS(128) is vjsr in component cell_geometry (microliter).
C CONSTS(129) is vcsr in component cell_geometry (microliter).
C CONSTS(131) is vsl in component cell_geometry (microliter).
C CONSTS(130) is vss in component cell_geometry (microliter).
C STATES(1) is v in component membrane (millivolt).
C ALGBRC(37) is vffrt in component membrane (coulomb_per_mole).
C ALGBRC(45) is vfrt in component membrane (dimensionless).
C ALGBRC(64) is INa in component INa (microA_per_microF).
C ALGBRC(66) is INaL in component INaL (microA_per_microF).
C ALGBRC(67) is Ito in component Ito (microA_per_microF).
C ALGBRC(69) is Isus in component Isus (microA_per_microF).
C ALGBRC(80) is ICaL in component ICaL (microA_per_microF).
C ALGBRC(83) is ICaT in component ICaT (microA_per_microF).
C ALGBRC(81) is ICaNa in component ICaL (microA_per_microF).
C ALGBRC(82) is ICaK in component ICaL (microA_per_microF).
C ALGBRC(88) is IKr in component IKr (microA_per_microF).
C ALGBRC(90) is IKs in component IKs (microA_per_microF).
C ALGBRC(93) is If in component If (microA_per_microF).
C ALGBRC(95) is IK1 in component IK1 (microA_per_microF).
C ALGBRC(127) is INaCa_i in component INaCa_i (microA_per_microF).
C ALGBRC(157) is INaCa_ss in component INaCa_i (microA_per_microF).
C ALGBRC(176) is INaK in component INaK (microA_per_microF).
C ALGBRC(177) is INab in component INab (microA_per_microF).
C ALGBRC(183) is IpCa in component IpCa (microA_per_microF).
C ALGBRC(180) is ICab in component ICab (microA_per_microF).
C ALGBRC(16) is Istim in component membrane (microA_per_microF).
C CONSTS(13) is amp in component membrane (microA_per_microF).
C CONSTS(14) is duration in component membrane (millisecond).
C CONSTS(15) is KmCaMK in component CaMK (millimolar).
C CONSTS(16) is aCaMK in component CaMK (per_millimolar_per_millisecond).
C CONSTS(17) is bCaMK in component CaMK (per_millisecond).
C CONSTS(18) is CaMKo in component CaMK (dimensionless).
C CONSTS(19) is KmCaM in component CaMK (millimolar).
C ALGBRC(50) is CaMKb in component CaMK (millimolar).
C ALGBRC(51) 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(20) is BSRmax in component intracellular_ions (millimolar).
C CONSTS(21) is KmBSR in component intracellular_ions (millimolar).
C CONSTS(22) is BSLmax in component intracellular_ions (millimolar).
C CONSTS(23) is KmBSL in component intracellular_ions (millimolar).
C CONSTS(24) is csqnmax in component intracellular_ions (millimolar).
C CONSTS(25) is kmcsqn in component intracellular_ions (millimolar).
C CONSTS(26) is csqnmaxsl in component intracellular_ions (millimolar).
C CONSTS(27) is cmdnmax in component intracellular_ions (millimolar).
C CONSTS(28) is kmcmdn in component intracellular_ions (millimolar).
C CONSTS(29) is cmdnmaxsl in component intracellular_ions (millimolar).
C CONSTS(30) is trpnmax in component intracellular_ions (millimolar).
C CONSTS(31) is kmtrpn in component intracellular_ions (millimolar).
C CONSTS(32) is trpnmaxsl in component intracellular_ions (millimolar).
C STATES(4) is nai in component intracellular_ions (millimolar).
C STATES(5) is nasl in component intracellular_ions (millimolar).
C STATES(6) is nass in component intracellular_ions (millimolar).
C STATES(7) is ki in component intracellular_ions (millimolar).
C STATES(8) is kss in component intracellular_ions (millimolar).
C STATES(9) is ksl in component intracellular_ions (millimolar).
C STATES(10) is cai in component intracellular_ions (millimolar).
C STATES(11) is casl in component intracellular_ions (millimolar).
C STATES(12) is cansr in component intracellular_ions (millimolar).
C STATES(13) is cajsr in component intracellular_ions (millimolar).
C STATES(14) is cacsr in component intracellular_ions (millimolar).
C ALGBRC(91) is IfNa in component If (microA_per_microF).
C ALGBRC(92) is IfK in component If (microA_per_microF).
C ALGBRC(179) is JdiffNa in component diff (millimolar_per_millisecond).
C ALGBRC(182) is JgapNa in component diff (millimolar_per_millisecond).
C ALGBRC(184) is Jdiff in component diff (millimolar_per_millisecond).
C ALGBRC(185) is Jgap in component diff (millimolar_per_millisecond).
C ALGBRC(178) is JdiffK in component diff (millimolar_per_millisecond).
C ALGBRC(181) is JgapK in component diff (millimolar_per_millisecond).
C ALGBRC(190) is Jup1 in component SERCA (millimolar_per_millisecond).
C ALGBRC(191) is Jup2 in component SERCA (millimolar_per_millisecond).
C STATES(15) is Jrel1 in component ryr (millimolar_per_millisecond).
C STATES(16) is Jrel2 in component ryr (millimolar_per_millisecond).
C ALGBRC(194) is Jip3 in component IP3 (millimolar_per_millisecond).
C ALGBRC(196) is Jtr1 in component trans_flux (millimolar_per_millisecond).
C ALGBRC(198) is Jtr2 in component trans_flux (millimolar_per_millisecond).
C ALGBRC(52) is Bcai in component intracellular_ions (dimensionless).
C ALGBRC(55) is Bcajsr in component intracellular_ions (dimensionless).
C ALGBRC(56) is Bcacsr in component intracellular_ions (dimensionless).
C ALGBRC(53) is Bcass in component intracellular_ions (dimensionless).
C ALGBRC(54) is Bcasl in component intracellular_ions (dimensionless).
C CONSTS(33) is cm in component intracellular_ions (microF_per_centimeter_squared).
C CONSTS(34) is PKNa in component reversal_potentials (dimensionless).
C ALGBRC(57) is ENa in component reversal_potentials (millivolt).
C ALGBRC(58) is EK in component reversal_potentials (millivolt).
C ALGBRC(59) is ECa in component reversal_potentials (millivolt).
C ALGBRC(60) is EKs in component reversal_potentials (millivolt).
C ALGBRC(1) is mss in component INa (dimensionless).
C ALGBRC(17) is tm in component INa (millisecond).
C CONSTS(35) is mssV1 in component INa (millivolt).
C CONSTS(36) is mssV2 in component INa (millivolt).
C CONSTS(37) is mtV1 in component INa (millivolt).
C CONSTS(38) is mtV2 in component INa (millivolt).
C CONSTS(39) is mtD1 in component INa (dimensionless).
C CONSTS(40) is mtD2 in component INa (dimensionless).
C CONSTS(41) is mtV3 in component INa (millivolt).
C CONSTS(42) is mtV4 in component INa (millivolt).
C STATES(17) is m in component INa (dimensionless).
C ALGBRC(2) is hss in component INa (dimensionless).
C ALGBRC(18) is thf in component INa (millisecond).
C ALGBRC(19) is ths in component INa (millisecond).
C CONSTS(43) is hssV1 in component INa (millivolt).
C CONSTS(44) is hssV2 in component INa (millivolt).
C CONSTS(113) is Ahs in component INa (dimensionless).
C CONSTS(45) is Ahf in component INa (dimensionless).
C STATES(18) is hf in component INa (dimensionless).
C STATES(19) is hs in component INa (dimensionless).
C ALGBRC(61) is h in component INa (dimensionless).
C CONSTS(46) is GNa in component INa (milliS_per_microF).
C ALGBRC(20) is jss in component INa (dimensionless).
C ALGBRC(38) is tj in component INa (millisecond).
C STATES(20) is j in component INa (dimensionless).
C ALGBRC(39) is hssp in component INa (dimensionless).
C ALGBRC(46) is thsp in component INa (millisecond).
C STATES(21) is hsp in component INa (dimensionless).
C ALGBRC(62) is hp in component INa (dimensionless).
C ALGBRC(47) is tjp in component INa (millisecond).
C STATES(22) is jp in component INa (dimensionless).
C ALGBRC(63) is fINap in component INa (dimensionless).
C ALGBRC(40) is mLss in component INaL (dimensionless).
C ALGBRC(48) is tmL in component INaL (millisecond).
C STATES(23) is mL in component INaL (dimensionless).
C CONSTS(47) is thL in component INaL (millisecond).
C ALGBRC(3) is hLss in component INaL (dimensionless).
C STATES(24) is hL in component INaL (dimensionless).
C ALGBRC(4) is hLssp in component INaL (dimensionless).
C CONSTS(114) is thLp in component INaL (millisecond).
C STATES(25) is hLp in component INaL (dimensionless).
C CONSTS(48) is GNaL in component INaL (milliS_per_microF).
C ALGBRC(65) is fINaLp in component INaL (dimensionless).
C ALGBRC(5) is ass in component Ito (dimensionless).
C ALGBRC(21) is taua in component Ito (millisecond).
C ALGBRC(6) is iss in component Ito (dimensionless).
C ALGBRC(22) is tauis in component Ito (millisecond).
C ALGBRC(23) is tauif in component Ito (millisecond).
C STATES(26) is a in component Ito (dimensionless).
C STATES(27) is i1 in component Ito (dimensionless).
C STATES(28) is i2 in component Ito (dimensionless).
C CONSTS(49) is Gto in component Ito (milliS_per_microF).
C CONSTS(50) is Gsus in component Isus (milliS_per_microF).
C ALGBRC(68) is asus in component Isus (dimensionless).
C CONSTS(51) is Kmn in component ICaL (millimolar).
C CONSTS(52) is k2n in component ICaL (per_millisecond).
C ALGBRC(7) is dss in component ICaL (dimensionless).
C STATES(29) is d in component ICaL (dimensionless).
C ALGBRC(8) is fss in component ICaL (dimensionless).
C CONSTS(115) is Aff in component ICaL (dimensionless).
C CONSTS(121) is Afs in component ICaL (dimensionless).
C STATES(30) is ff in component ICaL (dimensionless).
C STATES(31) is fs in component ICaL (dimensionless).
C ALGBRC(70) is f in component ICaL (dimensionless).
C ALGBRC(24) is fcass in component ICaL (dimensionless).
C ALGBRC(71) is Afcaf in component ICaL (dimensionless).
C ALGBRC(72) is Afcas in component ICaL (dimensionless).
C STATES(32) is fcaf in component ICaL (dimensionless).
C STATES(33) is fcas in component ICaL (dimensionless).
C ALGBRC(73) is fca in component ICaL (dimensionless).
C STATES(34) is jca in component ICaL (dimensionless).
C STATES(35) is ffp in component ICaL (dimensionless).
C ALGBRC(74) is fp in component ICaL (dimensionless).
C STATES(36) is fcafp in component ICaL (dimensionless).
C ALGBRC(75) is fcap in component ICaL (dimensionless).
C ALGBRC(9) is km2n in component ICaL (per_millisecond).
C ALGBRC(25) is anca in component ICaL (dimensionless).
C STATES(37) is nca in component ICaL (dimensionless).
C ALGBRC(76) is PhiCaL in component ICaL (dimensionless).
C ALGBRC(77) is PhiCaNa in component ICaL (dimensionless).
C ALGBRC(78) is PhiCaK in component ICaL (dimensionless).
C CONSTS(53) is PCa in component ICaL (dimensionless).
C CONSTS(116) is PCap in component ICaL (dimensionless).
C CONSTS(117) is PCaNa in component ICaL (dimensionless).
C CONSTS(118) is PCaK in component ICaL (dimensionless).
C CONSTS(122) is PCaNap in component ICaL (dimensionless).
C CONSTS(123) is PCaKp in component ICaL (dimensionless).
C ALGBRC(79) is fICaLp in component ICaL (dimensionless).
C ALGBRC(26) is td in component ICaL (millisecond).
C ALGBRC(27) is tff in component ICaL (millisecond).
C ALGBRC(28) is tfs in component ICaL (millisecond).
C ALGBRC(41) is tfcaf in component ICaL (millisecond).
C ALGBRC(42) is tfcas in component ICaL (millisecond).
C CONSTS(119) is tjca in component ICaL (millisecond).
C ALGBRC(43) is tffp in component ICaL (millisecond).
C ALGBRC(49) is tfcafp in component ICaL (millisecond).
C CONSTS(54) is GCaT in component ICaT (milliS_per_microF).
C ALGBRC(10) is bss in component ICaT (dimensionless).
C ALGBRC(11) is gss in component ICaT (dimensionless).
C ALGBRC(29) is taub in component ICaT (millisecond).
C ALGBRC(30) is taug in component ICaT (millisecond).
C STATES(38) is b in component ICaT (dimensionless).
C STATES(39) is g in component ICaT (dimensionless).
C CONSTS(55) is GKr in component IKr (milliS_per_microF).
C ALGBRC(12) is xrss in component IKr (dimensionless).
C ALGBRC(31) is txrf in component IKr (millisecond).
C ALGBRC(32) is txrs in component IKr (millisecond).
C ALGBRC(84) is Axrf in component IKr (dimensionless).
C ALGBRC(85) is Axrs in component IKr (dimensionless).
C STATES(40) is xrf in component IKr (dimensionless).
C STATES(41) is xrs in component IKr (dimensionless).
C ALGBRC(86) is xr in component IKr (dimensionless).
C ALGBRC(87) is rkr in component IKr (dimensionless).
C CONSTS(56) is GKs in component IKs (milliS_per_microF).
C ALGBRC(13) is xs1ss in component IKs (dimensionless).
C ALGBRC(33) is xs2ss in component IKs (dimensionless).
C ALGBRC(34) is txs1 in component IKs (millisecond).
C STATES(42) is xs1 in component IKs (dimensionless).
C STATES(43) is xs2 in component IKs (dimensionless).
C ALGBRC(89) is KsCa in component IKs (dimensionless).
C ALGBRC(44) is txs2 in component IKs (millisecond).
C CONSTS(57) is GfNa in component If (milliS_per_microF).
C CONSTS(58) is GfK in component If (milliS_per_microF).
C ALGBRC(14) is yss in component If (dimensionless).
C ALGBRC(35) is tauy in component If (millisecond).
C STATES(44) is y in component If (dimensionless).
C CONSTS(59) is GK1 in component IK1 (milliS_per_microF).
C ALGBRC(15) is xk1ss in component IK1 (dimensionless).
C ALGBRC(36) is txk1 in component IK1 (millisecond).
C STATES(45) is xk1 in component IK1 (dimensionless).
C ALGBRC(94) is rk1 in component IK1 (millisecond).
C CONSTS(60) is kna1 in component INaCa_i (per_millisecond).
C CONSTS(61) is kna2 in component INaCa_i (per_millisecond).
C CONSTS(62) is kna3 in component INaCa_i (per_millisecond).
C CONSTS(63) is kasymm in component INaCa_i (dimensionless).
C CONSTS(64) is wna in component INaCa_i (dimensionless).
C CONSTS(65) is wca in component INaCa_i (dimensionless).
C CONSTS(66) is wnaca in component INaCa_i (dimensionless).
C CONSTS(67) is kcaon in component INaCa_i (per_millisecond).
C CONSTS(68) is kcaoff in component INaCa_i (per_millisecond).
C CONSTS(69) is qna in component INaCa_i (dimensionless).
C CONSTS(70) is qca in component INaCa_i (dimensionless).
C ALGBRC(97) is hna in component INaCa_i (dimensionless).
C ALGBRC(96) is hca in component INaCa_i (dimensionless).
C CONSTS(71) is KmCaAct in component INaCa_i (millimolar).
C CONSTS(72) is Gncx in component INaCa_i (milliS_per_microF).
C ALGBRC(98) is h1_i in component INaCa_i (dimensionless).
C ALGBRC(99) is h2_i in component INaCa_i (dimensionless).
C ALGBRC(100) is h3_i in component INaCa_i (dimensionless).
C ALGBRC(101) is h4_i in component INaCa_i (dimensionless).
C ALGBRC(102) is h5_i in component INaCa_i (dimensionless).
C ALGBRC(103) is h6_i in component INaCa_i (dimensionless).
C ALGBRC(104) is h7_i in component INaCa_i (dimensionless).
C ALGBRC(105) is h8_i in component INaCa_i (dimensionless).
C ALGBRC(106) is h9_i in component INaCa_i (dimensionless).
C CONSTS(132) is h10_i in component INaCa_i (dimensionless).
C CONSTS(133) is h11_i in component INaCa_i (dimensionless).
C CONSTS(134) is h12_i in component INaCa_i (dimensionless).
C CONSTS(135) is k1_i in component INaCa_i (dimensionless).
C CONSTS(136) is k2_i in component INaCa_i (dimensionless).
C ALGBRC(107) is k3p_i in component INaCa_i (dimensionless).
C ALGBRC(108) is k3pp_i in component INaCa_i (dimensionless).
C ALGBRC(109) is k3_i in component INaCa_i (dimensionless).
C ALGBRC(112) is k4_i in component INaCa_i (dimensionless).
C ALGBRC(110) is k4p_i in component INaCa_i (dimensionless).
C ALGBRC(111) is k4pp_i in component INaCa_i (dimensionless).
C CONSTS(137) is k5_i in component INaCa_i (dimensionless).
C ALGBRC(113) is k6_i in component INaCa_i (dimensionless).
C ALGBRC(114) is k7_i in component INaCa_i (dimensionless).
C ALGBRC(115) is k8_i in component INaCa_i (dimensionless).
C ALGBRC(116) is x1_i in component INaCa_i (dimensionless).
C ALGBRC(117) is x2_i in component INaCa_i (dimensionless).
C ALGBRC(118) is x3_i in component INaCa_i (dimensionless).
C ALGBRC(119) is x4_i in component INaCa_i (dimensionless).
C ALGBRC(120) is E1_i in component INaCa_i (dimensionless).
C ALGBRC(121) is E2_i in component INaCa_i (dimensionless).
C ALGBRC(122) is E3_i in component INaCa_i (dimensionless).
C ALGBRC(123) is E4_i in component INaCa_i (dimensionless).
C ALGBRC(124) is allo_i in component INaCa_i (dimensionless).
C ALGBRC(125) is JncxNa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(126) is JncxCa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(128) is h1_ss in component INaCa_i (dimensionless).
C ALGBRC(129) is h2_ss in component INaCa_i (dimensionless).
C ALGBRC(130) is h3_ss in component INaCa_i (dimensionless).
C ALGBRC(131) is h4_ss in component INaCa_i (dimensionless).
C ALGBRC(132) is h5_ss in component INaCa_i (dimensionless).
C ALGBRC(133) is h6_ss in component INaCa_i (dimensionless).
C ALGBRC(134) is h7_ss in component INaCa_i (dimensionless).
C ALGBRC(135) is h8_ss in component INaCa_i (dimensionless).
C ALGBRC(136) is h9_ss in component INaCa_i (dimensionless).
C CONSTS(138) is h10_ss in component INaCa_i (dimensionless).
C CONSTS(139) is h11_ss in component INaCa_i (dimensionless).
C CONSTS(140) is h12_ss in component INaCa_i (dimensionless).
C CONSTS(141) is k1_ss in component INaCa_i (dimensionless).
C CONSTS(142) is k2_ss in component INaCa_i (dimensionless).
C ALGBRC(137) is k3p_ss in component INaCa_i (dimensionless).
C ALGBRC(138) is k3pp_ss in component INaCa_i (dimensionless).
C ALGBRC(139) is k3_ss in component INaCa_i (dimensionless).
C ALGBRC(142) is k4_ss in component INaCa_i (dimensionless).
C ALGBRC(140) is k4p_ss in component INaCa_i (dimensionless).
C ALGBRC(141) is k4pp_ss in component INaCa_i (dimensionless).
C CONSTS(143) is k5_ss in component INaCa_i (dimensionless).
C ALGBRC(143) is k6_ss in component INaCa_i (dimensionless).
C ALGBRC(144) is k7_ss in component INaCa_i (dimensionless).
C ALGBRC(145) is k8_ss in component INaCa_i (dimensionless).
C ALGBRC(146) is x1_ss in component INaCa_i (dimensionless).
C ALGBRC(147) is x2_ss in component INaCa_i (dimensionless).
C ALGBRC(148) is x3_ss in component INaCa_i (dimensionless).
C ALGBRC(149) is x4_ss in component INaCa_i (dimensionless).
C ALGBRC(150) is E1_ss in component INaCa_i (dimensionless).
C ALGBRC(151) is E2_ss in component INaCa_i (dimensionless).
C ALGBRC(152) is E3_ss in component INaCa_i (dimensionless).
C ALGBRC(153) is E4_ss in component INaCa_i (dimensionless).
C ALGBRC(154) is allo_ss in component INaCa_i (dimensionless).
C ALGBRC(155) is JncxNa_ss in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(156) is JncxCa_ss in component INaCa_i (millimolar_per_millisecond).
C CONSTS(73) is k1p in component INaK (per_millisecond).
C CONSTS(74) is k1m in component INaK (per_millisecond).
C CONSTS(75) is k2p in component INaK (per_millisecond).
C CONSTS(76) is k2m in component INaK (per_millisecond).
C CONSTS(77) is k3p in component INaK (per_millisecond).
C CONSTS(78) is k3m in component INaK (per_millisecond).
C CONSTS(79) is k4p in component INaK (per_millisecond).
C CONSTS(80) is k4m in component INaK (per_millisecond).
C CONSTS(81) is Knai0 in component INaK (millimolar).
C CONSTS(82) is Knao0 in component INaK (millimolar).
C CONSTS(83) is delta in component INaK (millivolt).
C CONSTS(84) is Kki in component INaK (per_millisecond).
C CONSTS(85) is Kko in component INaK (per_millisecond).
C CONSTS(86) is MgADP in component INaK (millimolar).
C CONSTS(87) is MgATP in component INaK (millimolar).
C CONSTS(88) is Kmgatp in component INaK (millimolar).
C CONSTS(89) is H in component INaK (millimolar).
C CONSTS(90) is eP in component INaK (dimensionless).
C CONSTS(91) is Khp in component INaK (millimolar).
C CONSTS(92) is Knap in component INaK (millimolar).
C CONSTS(93) is Kxkur in component INaK (millimolar).
C CONSTS(94) is Pnak in component INaK (milliS_per_microF).
C ALGBRC(158) is Knai in component INaK (millimolar).
C ALGBRC(159) is Knao in component INaK (millimolar).
C ALGBRC(160) is P in component INaK (dimensionless).
C ALGBRC(161) is a1 in component INaK (dimensionless).
C CONSTS(144) is b1 in component INaK (dimensionless).
C CONSTS(145) is a2 in component INaK (dimensionless).
C ALGBRC(162) is b2 in component INaK (dimensionless).
C ALGBRC(163) is a3 in component INaK (dimensionless).
C ALGBRC(164) is b3 in component INaK (dimensionless).
C CONSTS(146) is a4 in component INaK (dimensionless).
C ALGBRC(165) is b4 in component INaK (dimensionless).
C ALGBRC(166) is x1 in component INaK (dimensionless).
C ALGBRC(167) is x2 in component INaK (dimensionless).
C ALGBRC(168) is x3 in component INaK (dimensionless).
C ALGBRC(169) is x4 in component INaK (dimensionless).
C ALGBRC(170) is E1 in component INaK (dimensionless).
C ALGBRC(171) is E2 in component INaK (dimensionless).
C ALGBRC(172) is E3 in component INaK (dimensionless).
C ALGBRC(173) is E4 in component INaK (dimensionless).
C ALGBRC(174) is JnakNa in component INaK (millimolar_per_millisecond).
C ALGBRC(175) is JnakK in component INaK (millimolar_per_millisecond).
C CONSTS(95) is PNab in component INab (milliS_per_microF).
C CONSTS(96) is PCab in component ICab (milliS_per_microF).
C CONSTS(97) is GpCa in component IpCa (milliS_per_microF).
C CONSTS(98) is KmCap in component IpCa (millimolar).
C CONSTS(99) is sstau in component diff (millisecond).
C CONSTS(100) is gaptau in component diff (millisecond).
C ALGBRC(197) is REL in component ryr (millimolar_per_millisecond).
C ALGBRC(199) is irelss in component ryr (millimolar_per_millisecond).
C ALGBRC(186) is ireltau in component ryr (dimensionless).
C ALGBRC(193) is REL2 in component ryr (millimolar_per_millisecond).
C ALGBRC(195) is irelss2 in component ryr (millimolar_per_millisecond).
C ALGBRC(187) is ireltau2 in component ryr (dimensionless).
C CONSTS(101) is dqupcamkbar in component SERCA (dimensionless).
C CONSTS(102) is dkmplbbar in component SERCA (dimensionless).
C CONSTS(103) is kmup in component SERCA (dimensionless).
C CONSTS(104) is nsrbar in component SERCA (dimensionless).
C ALGBRC(188) is dkmplb in component SERCA (dimensionless).
C ALGBRC(189) is dqupcamk in component SERCA (dimensionless).
C CONSTS(105) is IP3 in component IP3 (dimensionless).
C CONSTS(106) is k1 in component IP3 (dimensionless).
C CONSTS(107) is k1a in component IP3 (dimensionless).
C CONSTS(108) is k0 in component IP3 (dimensionless).
C CONSTS(109) is k0a in component IP3 (dimensionless).
C CONSTS(110) is k2 in component IP3 (dimensionless).
C CONSTS(111) is k2a in component IP3 (dimensionless).
C CONSTS(112) is tauip3 in component IP3 (millisecond).
C ALGBRC(192) is POip3 in component IP3 (dimensionless).
C STATES(46) is u in component IP3 (millimolar_per_millisecond).
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 nasl in component intracellular_ions (millimolar).
C RATES(6) is d/dt nass in component intracellular_ions (millimolar).
C RATES(7) is d/dt ki in component intracellular_ions (millimolar).
C RATES(9) is d/dt ksl in component intracellular_ions (millimolar).
C RATES(8) 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(11) is d/dt casl in component intracellular_ions (millimolar).
C RATES(12) is d/dt cansr in component intracellular_ions (millimolar).
C RATES(13) is d/dt cajsr in component intracellular_ions (millimolar).
C RATES(14) is d/dt cacsr in component intracellular_ions (millimolar).
C RATES(17) is d/dt m in component INa (dimensionless).
C RATES(18) is d/dt hf in component INa (dimensionless).
C RATES(19) is d/dt hs in component INa (dimensionless).
C RATES(20) is d/dt j in component INa (dimensionless).
C RATES(21) is d/dt hsp in component INa (dimensionless).
C RATES(22) is d/dt jp in component INa (dimensionless).
C RATES(23) is d/dt mL in component INaL (dimensionless).
C RATES(24) is d/dt hL in component INaL (dimensionless).
C RATES(25) is d/dt hLp in component INaL (dimensionless).
C RATES(26) is d/dt a in component Ito (dimensionless).
C RATES(27) is d/dt i1 in component Ito (dimensionless).
C RATES(28) is d/dt i2 in component Ito (dimensionless).
C RATES(29) is d/dt d in component ICaL (dimensionless).
C RATES(30) is d/dt ff in component ICaL (dimensionless).
C RATES(31) is d/dt fs in component ICaL (dimensionless).
C RATES(32) is d/dt fcaf in component ICaL (dimensionless).
C RATES(33) is d/dt fcas in component ICaL (dimensionless).
C RATES(34) is d/dt jca in component ICaL (dimensionless).
C RATES(35) is d/dt ffp in component ICaL (dimensionless).
C RATES(36) is d/dt fcafp in component ICaL (dimensionless).
C RATES(37) is d/dt nca in component ICaL (dimensionless).
C RATES(38) is d/dt b in component ICaT (dimensionless).
C RATES(39) is d/dt g in component ICaT (dimensionless).
C RATES(40) is d/dt xrf in component IKr (dimensionless).
C RATES(41) is d/dt xrs in component IKr (dimensionless).
C RATES(42) is d/dt xs1 in component IKs (dimensionless).
C RATES(43) is d/dt xs2 in component IKs (dimensionless).
C RATES(44) is d/dt y in component If (dimensionless).
C RATES(45) is d/dt xk1 in component IK1 (dimensionless).
C RATES(15) is d/dt Jrel1 in component ryr (millimolar_per_millisecond).
C RATES(16) is d/dt Jrel2 in component ryr (millimolar_per_millisecond).
C RATES(46) is d/dt u in component IP3 (millimolar_per_millisecond).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 140
      CONSTS(2) = 1.8
      CONSTS(3) = 5.4
      CONSTS(4) = 8314
      CONSTS(5) = 310
      CONSTS(6) = 96485
      CONSTS(7) = 1
      CONSTS(8) = 2
      CONSTS(9) = 1
      CONSTS(10) = 0.0164
      CONSTS(11) = 0.00175
      CONSTS(12) = 3.14159265
      STATES(1) = -86.6814002878592
      CONSTS(13) = -40
      CONSTS(14) = 1
      CONSTS(15) = 0.15
      CONSTS(16) = 0.05
      CONSTS(17) = 0.00068
      CONSTS(18) = 0.05
      CONSTS(19) = 0.0015
      STATES(2) = 0.00505983330678751
      STATES(3) = 0.000101777993438818
      CONSTS(20) = 0.019975
      CONSTS(21) = 0.00087
      CONSTS(22) = 0.4777
      CONSTS(23) = 0.0087
      CONSTS(24) = 2.88
      CONSTS(25) = 0.8
      CONSTS(26) = 1.2
      CONSTS(27) = 0.1125
      CONSTS(28) = 0.00238
      CONSTS(29) = 0.0125
      CONSTS(30) = 0.0315
      CONSTS(31) = 0.0005
      CONSTS(32) = 0.0035
      STATES(4) = 8.23183964616932
      STATES(5) = 8.23153516580562
      STATES(6) = 8.23154325237268
      STATES(7) = 143.767359809132
      STATES(8) = 143.767768218104
      STATES(9) = 143.767769906216
      STATES(10) = 4.36004404734282e-5
      STATES(11) = 0.000102004317781147
      STATES(12) = 1.26350902016858
      STATES(13) = 1.24811940209535
      STATES(14) = 1.26516959198518
      STATES(15) = 0.000108240945806962
      STATES(16) = 1.25045800437317e-69
      CONSTS(33) = 1
      CONSTS(34) = 0.01833
      CONSTS(35) = 48.4264
      CONSTS(36) = 7.5653
      CONSTS(37) = 11.64
      CONSTS(38) = 34.77
      CONSTS(39) = 6.765
      CONSTS(40) = 8.552
      CONSTS(41) = 77.42
      CONSTS(42) = 5.955
      STATES(17) = 0.00632661703915808
      CONSTS(43) = 78.5
      CONSTS(44) = 6.22
      CONSTS(45) = 0.99
      STATES(18) = 0.788611739889677
      STATES(19) = 0.788545979951331
      CONSTS(46) = 39.4572
      STATES(20) = 0.790474358603666
      STATES(21) = 0.579693514309867
      STATES(22) = 0.790947058236417
      STATES(23) = 0.000241925773627233
      CONSTS(47) = 200
      STATES(24) = 0.463574582508218
      STATES(25) = 0.240216198686475
      CONSTS(48) = 0.0189
      STATES(26) = 0.000272851144435704
      STATES(27) = 0.649604795721571
      STATES(28) = 0.989965695822495
      CONSTS(49) = 0.192
      CONSTS(50) = 0.0301
      CONSTS(51) = 0.002
      CONSTS(52) = 1000
      STATES(29) = 6.97735089296892e-9
      STATES(30) = 0.999999968230738
      STATES(31) = 0.926692153319136
      STATES(32) = 0.99999996819573
      STATES(33) = 0.999999905741936
      STATES(34) = 0.999978907334662
      STATES(35) = 0.999999968365903
      STATES(36) = 0.999999968278239
      STATES(37) = 0.00547252500964926
      CONSTS(53) = 7.7677e-5
      CONSTS(54) = 0.0754
      STATES(38) = 0.000304250912559619
      STATES(39) = 0.994214357917907
      CONSTS(55) = 0.0342
      STATES(40) = 0.000331691184084272
      STATES(41) = 0.568716473334161
      CONSTS(56) = 0.0029
      STATES(42) = 0.191165248085394
      STATES(43) = 0.000222677365291219
      CONSTS(57) = 0.0116
      CONSTS(58) = 0.0232
      STATES(44) = 0.233119011214908
      CONSTS(59) = 0.0455
      STATES(45) = 0.997084813729909
      CONSTS(60) = 15
      CONSTS(61) = 5
      CONSTS(62) = 88.12
      CONSTS(63) = 12.5
      CONSTS(64) = 60000
      CONSTS(65) = 60000
      CONSTS(66) = 5000
      CONSTS(67) = 1500000
      CONSTS(68) = 5000
      CONSTS(69) = 0.5224
      CONSTS(70) = 0.167
      CONSTS(71) = 0.00015
      CONSTS(72) = 0.00095709
      CONSTS(73) = 949.5
      CONSTS(74) = 182.4
      CONSTS(75) = 687.2
      CONSTS(76) = 39.4
      CONSTS(77) = 1899
      CONSTS(78) = 79300
      CONSTS(79) = 639
      CONSTS(80) = 40
      CONSTS(81) = 9.073
      CONSTS(82) = 27.78
      CONSTS(83) = -0.155
      CONSTS(84) = 0.5
      CONSTS(85) = 0.3582
      CONSTS(86) = 0.05
      CONSTS(87) = 9.8
      CONSTS(88) = 1.698e-7
      CONSTS(89) = 1e-7
      CONSTS(90) = 4.2
      CONSTS(91) = 1.698e-7
      CONSTS(92) = 224
      CONSTS(93) = 292
      CONSTS(94) = 32.4872
      CONSTS(95) = 9.375e-10
      CONSTS(96) = 2.5e-8
      CONSTS(97) = 0.0005
      CONSTS(98) = 0.0005
      CONSTS(99) = 0.2
      CONSTS(100) = 12
      CONSTS(101) = 0.75
      CONSTS(102) = 0.00017
      CONSTS(103) = 0.00028
      CONSTS(104) = 15
      CONSTS(105) = 0.0001
      CONSTS(106) = 150000
      CONSTS(107) = 16.5
      CONSTS(108) = 96000
      CONSTS(109) = 9.6
      CONSTS(110) = 1800
      CONSTS(111) = 0.21
      CONSTS(112) = 3.7
      STATES(46) = 0.466236137183558
      CONSTS(113) = 1.00000 - CONSTS(45)
      CONSTS(114) =  3.00000*CONSTS(47)
      CONSTS(115) = 0.600000
      CONSTS(116) =  1.10000*CONSTS(53)
      CONSTS(117) =  0.00125000*CONSTS(53)
      CONSTS(118) =  0.000357400*CONSTS(53)
      CONSTS(119) = 75.0000
      CONSTS(120) =  1000.00*3.14159*CONSTS(11)*CONSTS(11)*CONSTS(10)
      CONSTS(121) = 1.00000 - CONSTS(115)
      CONSTS(122) =  0.00125000*CONSTS(116)
      CONSTS(123) =  0.000357400*CONSTS(116)
      CONSTS(124) =  2.00000*CONSTS(12)*CONSTS(11)*CONSTS(11)+ 2.00000*CONSTS(12)*CONSTS(11)*CONSTS(10)
      CONSTS(125) =  2.00000*CONSTS(124)
      CONSTS(126) =  0.600000*CONSTS(120)
      CONSTS(127) =  0.0400000*CONSTS(120)
      CONSTS(128) =  0.00200000*CONSTS(120)
      CONSTS(129) =  0.00800000*CONSTS(120)
      CONSTS(130) =  0.0200000*CONSTS(120)
      CONSTS(131) =  0.150000*CONSTS(120)
      CONSTS(132) = CONSTS(63)+1.00000+ (CONSTS(1)/CONSTS(60))*(1.00000+CONSTS(1)/CONSTS(61))
      CONSTS(133) = ( CONSTS(1)*CONSTS(1))/( CONSTS(132)*CONSTS(60)*CONSTS(61))
      CONSTS(134) = 1.00000/CONSTS(132)
      CONSTS(135) =  CONSTS(134)*CONSTS(2)*CONSTS(67)
      CONSTS(136) = CONSTS(68)
      CONSTS(137) = CONSTS(68)
      CONSTS(138) = CONSTS(63)+1.00000+ (CONSTS(1)/CONSTS(60))*(1.00000+CONSTS(1)/CONSTS(61))
      CONSTS(139) = ( CONSTS(1)*CONSTS(1))/( CONSTS(138)*CONSTS(60)*CONSTS(61))
      CONSTS(140) = 1.00000/CONSTS(138)
      CONSTS(141) =  CONSTS(140)*CONSTS(2)*CONSTS(67)
      CONSTS(142) = CONSTS(68)
      CONSTS(143) = CONSTS(68)
      CONSTS(144) =  CONSTS(74)*CONSTS(86)
      CONSTS(145) = CONSTS(75)
      CONSTS(146) = (( CONSTS(79)*CONSTS(87))/CONSTS(88))/(1.00000+CONSTS(87)/CONSTS(88))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(46) =  STATES(3)*CONSTS(110)*(1.00000 - STATES(46)) -  CONSTS(111)*STATES(46)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
      RATES(24) = (ALGBRC(3) - STATES(24))/CONSTS(47)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
      RATES(25) = (ALGBRC(4) - STATES(25))/CONSTS(114)
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(35))/CONSTS(36)))
      ALGBRC(17) = 1.00000/( CONSTS(39)*EXP((STATES(1)+CONSTS(37))/CONSTS(38))+ CONSTS(40)*EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
      RATES(17) = (ALGBRC(1) - STATES(17))/ALGBRC(17)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(43))/CONSTS(44)))
      ALGBRC(18) = 1.00000/( 3.68600e-06*EXP(- (STATES(1)+3.88750)/7.85790)+ 16.0000*EXP((STATES(1) - 0.496300)/9.18430))
      RATES(18) = (ALGBRC(2) - STATES(18))/ALGBRC(18)
      ALGBRC(19) = 1.00000/( 0.00979400*EXP(- (STATES(1)+17.9500)/28.0500)+ 0.334300*EXP((STATES(1)+5.73000)/56.6600))
      RATES(19) = (ALGBRC(2) - STATES(19))/ALGBRC(19)
      ALGBRC(5) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/13.0000))
      ALGBRC(21) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- (STATES(1) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814)))
      RATES(26) = (ALGBRC(5) - STATES(26))/ALGBRC(21)
      ALGBRC(6) = 1.00000/(1.00000+EXP((27.0000+STATES(1))/13.0000))
      ALGBRC(22) = 43.0000+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
      RATES(27) = (ALGBRC(6) - STATES(27))/ALGBRC(22)
      ALGBRC(23) = 6.16200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1) - 8.00000)/8.59000))
      RATES(28) = (ALGBRC(6) - STATES(28))/ALGBRC(23)
      ALGBRC(7) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000+3.30000)/4.23000))
      ALGBRC(26) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
      RATES(29) = (ALGBRC(7) - STATES(29))/ALGBRC(26)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+19.5800+3.30000)/3.69600))
      ALGBRC(27) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000+15.1900)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000+15.1900)/10.0000))
      RATES(30) = (ALGBRC(8) - STATES(30))/ALGBRC(27)
      ALGBRC(28) = 1000.00+1.00000/( 3.50000e-05*EXP(- (STATES(1)+5.00000+15.1900)/4.00000)+ 3.50000e-05*EXP((STATES(1)+5.00000+15.1900)/6.00000))
      RATES(31) = (ALGBRC(8) - STATES(31))/ALGBRC(28)
      ALGBRC(24) = ALGBRC(8)
      RATES(34) = (ALGBRC(24) - STATES(34))/CONSTS(119)
      ALGBRC(9) =  STATES(34)*1.00000
      ALGBRC(25) = 1.00000/(CONSTS(52)/ALGBRC(9)+1.00000+CONSTS(51)/STATES(3) ** 4.00000)
      RATES(37) =  ALGBRC(25)*CONSTS(52) -  STATES(37)*ALGBRC(9)
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+30.0000)/7.00000))
      ALGBRC(29) = 1.00000/( 1.06800*EXP((STATES(1)+16.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+16.3000)/30.0000))
      RATES(38) = (ALGBRC(10) - STATES(38))/ALGBRC(29)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+61.0000)/5.00000))
      ALGBRC(30) = 1.00000/( 0.0150000*EXP((STATES(1)+71.7000)/15.4000)+ 0.0150000*EXP(- (STATES(1)+71.7000)/83.3000))
      RATES(39) = (ALGBRC(11) - STATES(39))/ALGBRC(30)
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+8.33700)/6.78900))
      ALGBRC(31) = 12.9800+1.00000/( 0.365200*EXP(((STATES(1)+17.6000) - 31.6600)/3.86900)+ 4.12300e-05*EXP(- ((STATES(1)+17.6000) - 47.7800)/20.3800))
      RATES(40) = (ALGBRC(12) - STATES(40))/ALGBRC(31)
      ALGBRC(32) = 1.86500+1.00000/( 0.0662900*EXP(((STATES(1)+17.2000) - 34.7000)/7.35500)+ 1.12800e-05*EXP(- ((STATES(1)+17.2000) - 29.7400)/25.9400))
      RATES(41) = (ALGBRC(12) - STATES(41))/ALGBRC(32)
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(34) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      RATES(42) = (ALGBRC(13) - STATES(42))/ALGBRC(34)
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+87.0000)/9.50000))
      ALGBRC(35) = 2000.00/(EXP((STATES(1)+57.0000)/60.0000)+EXP(- (STATES(1)+132.000)/10.0000))
      RATES(44) = (ALGBRC(14) - STATES(44))/ALGBRC(35)
      ALGBRC(15) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(3)+144.590)/( 1.56920*CONSTS(3)+3.81150)))
      ALGBRC(36) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
      RATES(45) = (ALGBRC(15) - STATES(45))/ALGBRC(36)
      ALGBRC(20) = ALGBRC(2)
      ALGBRC(38) = 4.85900+1.00000/( 0.862800*EXP(- (STATES(1)+116.726)/7.60050)+ 1.10960*EXP((STATES(1)+6.27190)/9.03580))
      RATES(20) = (ALGBRC(20) - STATES(20))/ALGBRC(38)
      ALGBRC(41) =  0.720000*(7.00000+1.00000/( 0.0400000*EXP(- ((STATES(1)+15.1900) - 4.00000)/7.00000)+ 0.0400000*EXP(((STATES(1)+15.1900) - 4.00000)/7.00000)))
      RATES(32) = (ALGBRC(24) - STATES(32))/ALGBRC(41)
      ALGBRC(42) =  0.490000*(100.000+1.00000/( 0.000120000*EXP(- (STATES(1)+15.1900)/3.00000)+ 0.000120000*EXP((STATES(1)+15.1900)/7.00000)))
      RATES(33) = (ALGBRC(24) - STATES(33))/ALGBRC(42)
      ALGBRC(43) =  2.50000*ALGBRC(27)
      RATES(35) = (ALGBRC(8) - STATES(35))/ALGBRC(43)
      ALGBRC(33) = ALGBRC(13)
      ALGBRC(44) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      RATES(43) = (ALGBRC(33) - STATES(43))/ALGBRC(44)
      ALGBRC(50) = ( CONSTS(18)*(1.00000 - STATES(2)))/(1.00000+CONSTS(19)/STATES(3))
      RATES(2) =  CONSTS(16)*ALGBRC(50)*(ALGBRC(50)+STATES(2)) -  CONSTS(17)*STATES(2)
      ALGBRC(39) = 1.00000/(1.00000+EXP((STATES(1)+84.7000)/6.22000))
      ALGBRC(46) =  3.00000*ALGBRC(19)
      RATES(21) = (ALGBRC(39) - STATES(21))/ALGBRC(46)
      ALGBRC(47) =  1.46000*ALGBRC(38)
      RATES(22) = (ALGBRC(20) - STATES(22))/ALGBRC(47)
      ALGBRC(40) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(48) = ALGBRC(17)
      RATES(23) = (ALGBRC(40) - STATES(23))/ALGBRC(48)
      ALGBRC(49) =  2.50000*ALGBRC(41)
      RATES(36) = (ALGBRC(24) - STATES(36))/ALGBRC(49)
      ALGBRC(70) =  CONSTS(115)*STATES(30)+ CONSTS(121)*STATES(31)
      ALGBRC(71) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(72) = 1.00000 - ALGBRC(71)
      ALGBRC(73) =  ALGBRC(71)*STATES(32)+ ALGBRC(72)*STATES(33)
      ALGBRC(74) =  CONSTS(115)*STATES(35)+ CONSTS(121)*STATES(31)
      ALGBRC(75) =  ALGBRC(71)*STATES(36)+ ALGBRC(72)*STATES(33)
      ALGBRC(37) = ( STATES(1)*CONSTS(6)*CONSTS(6))/( CONSTS(4)*CONSTS(5))
      ALGBRC(45) = ( STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5))
      ALGBRC(78) = ( 1.00000*ALGBRC(37)*( 0.750000*STATES(8)*EXP( 1.00000*ALGBRC(45)) -  0.750000*CONSTS(3)))/(EXP( 1.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(51) = ALGBRC(50)+STATES(2)
      ALGBRC(79) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(82) =  (1.00000 - ALGBRC(79))*CONSTS(118)*ALGBRC(78)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(123)*ALGBRC(78)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(178) = (STATES(8) - STATES(9))/CONSTS(99)
      RATES(8) = ( - ALGBRC(82)*CONSTS(33)*CONSTS(125))/( CONSTS(6)*CONSTS(130)) - ALGBRC(178)
      ALGBRC(77) = ( 1.00000*ALGBRC(37)*( 0.750000*STATES(6)*EXP( 1.00000*ALGBRC(45)) -  0.750000*CONSTS(1)))/(EXP( 1.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(81) =  (1.00000 - ALGBRC(79))*CONSTS(117)*ALGBRC(77)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(122)*ALGBRC(77)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(154) = 1.00000/(1.00000+CONSTS(71)/STATES(3) ** 2.00000)
      ALGBRC(97) = EXP(( CONSTS(69)*STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5)))
      ALGBRC(134) = 1.00000+ (CONSTS(1)/CONSTS(62))*(1.00000+1.00000/ALGBRC(97))
      ALGBRC(135) = CONSTS(1)/( CONSTS(62)*ALGBRC(97)*ALGBRC(134))
      ALGBRC(138) =  ALGBRC(135)*CONSTS(66)
      ALGBRC(128) = 1.00000+ (STATES(6)/CONSTS(62))*(1.00000+ALGBRC(97))
      ALGBRC(129) = ( STATES(6)*ALGBRC(97))/( CONSTS(62)*ALGBRC(128))
      ALGBRC(141) =  ALGBRC(129)*CONSTS(66)
      ALGBRC(131) = 1.00000+ (STATES(6)/CONSTS(60))*(1.00000+STATES(6)/CONSTS(61))
      ALGBRC(132) = ( STATES(6)*STATES(6))/( ALGBRC(131)*CONSTS(60)*CONSTS(61))
      ALGBRC(144) =  ALGBRC(132)*ALGBRC(129)*CONSTS(64)
      ALGBRC(145) =  ALGBRC(135)*CONSTS(139)*CONSTS(64)
      ALGBRC(136) = 1.00000/ALGBRC(134)
      ALGBRC(137) =  ALGBRC(136)*CONSTS(65)
      ALGBRC(139) = ALGBRC(137)+ALGBRC(138)
      ALGBRC(96) = EXP(( CONSTS(70)*STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5)))
      ALGBRC(130) = 1.00000/ALGBRC(128)
      ALGBRC(140) = ( ALGBRC(130)*CONSTS(65))/ALGBRC(96)
      ALGBRC(142) = ALGBRC(140)+ALGBRC(141)
      ALGBRC(133) = 1.00000/ALGBRC(131)
      ALGBRC(143) =  ALGBRC(133)*STATES(3)*CONSTS(67)
      ALGBRC(146) =  CONSTS(142)*ALGBRC(142)*(ALGBRC(144)+ALGBRC(143))+ CONSTS(143)*ALGBRC(144)*(CONSTS(142)+ALGBRC(139))
      ALGBRC(147) =  CONSTS(141)*ALGBRC(144)*(ALGBRC(142)+CONSTS(143))+ ALGBRC(142)*ALGBRC(143)*(CONSTS(141)+ALGBRC(145))
      ALGBRC(148) =  CONSTS(141)*ALGBRC(139)*(ALGBRC(144)+ALGBRC(143))+ ALGBRC(145)*ALGBRC(143)*(CONSTS(142)+ALGBRC(139))
      ALGBRC(149) =  CONSTS(142)*ALGBRC(145)*(ALGBRC(142)+CONSTS(143))+ ALGBRC(139)*CONSTS(143)*(CONSTS(141)+ALGBRC(145))
      ALGBRC(150) = ALGBRC(146)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(151) = ALGBRC(147)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(152) = ALGBRC(148)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(153) = ALGBRC(149)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(155) = ( 3.00000*( ALGBRC(153)*ALGBRC(144) -  ALGBRC(150)*ALGBRC(145))+ ALGBRC(152)*ALGBRC(141)) -  ALGBRC(151)*ALGBRC(138)
      ALGBRC(156) =  ALGBRC(151)*CONSTS(142) -  ALGBRC(150)*CONSTS(141)
      ALGBRC(157) =  0.200000*CONSTS(72)*ALGBRC(154)*( CONSTS(7)*ALGBRC(155)+ CONSTS(8)*ALGBRC(156))
      ALGBRC(179) = (STATES(6) - STATES(5))/CONSTS(99)
      RATES(6) = ( - (ALGBRC(81)+ 3.00000*ALGBRC(157))*CONSTS(33)*CONSTS(125))/( CONSTS(6)*CONSTS(130)) - ALGBRC(179)
      ALGBRC(181) = (STATES(9) - STATES(7))/CONSTS(100)
      RATES(7) = ( ALGBRC(181)*CONSTS(131))/CONSTS(126)
      ALGBRC(58) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log(CONSTS(3)/STATES(9))
      ALGBRC(67) =  CONSTS(49)*STATES(26)*STATES(27)*STATES(28)*(STATES(1) - ALGBRC(58))
      ALGBRC(68) = 1.00000/(1.00000+EXP(- (STATES(1) - 12.0000)/16.0000))
      ALGBRC(69) =  CONSTS(50)*ALGBRC(68)*(STATES(1) - ALGBRC(58))
      ALGBRC(84) = 1.00000/(1.00000+EXP((STATES(1)+54.8100)/38.2100))
      ALGBRC(85) = 1.00000 - ALGBRC(84)
      ALGBRC(86) =  ALGBRC(84)*STATES(40)+ ALGBRC(85)*STATES(41)
      ALGBRC(87) = ( (1.00000/(1.00000+EXP((STATES(1)+55.0000)/( 0.320000*75.0000))))*1.00000)/(1.00000+EXP((STATES(1) - 10.0000)/( 0.320000*30.0000)))
      ALGBRC(88) =  CONSTS(55)* (CONSTS(3)/5.40000) ** (1.0 / 2)*ALGBRC(86)*ALGBRC(87)*(STATES(1) - ALGBRC(58))
      ALGBRC(60) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log((CONSTS(3)+ CONSTS(34)*CONSTS(1))/(STATES(9)+ CONSTS(34)*STATES(5)))
      ALGBRC(89) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(11) ** 1.40000)
      ALGBRC(90) =  CONSTS(56)*ALGBRC(89)*STATES(42)*STATES(43)*(STATES(1) - ALGBRC(60))
      ALGBRC(94) = 1.00000/(1.00000+EXP(((STATES(1)+116.000) -  5.50000*CONSTS(3))/11.0000))
      ALGBRC(95) =  CONSTS(59)*2.32380* (CONSTS(3)/5.40000) ** (1.0 / 2)*ALGBRC(94)*STATES(45)*(STATES(1) - ALGBRC(58))
      ALGBRC(159) =  CONSTS(82)*EXP(( (1.00000 - CONSTS(83))*STATES(1)*CONSTS(6))/( 3.00000*CONSTS(4)*CONSTS(5)))
      ALGBRC(163) = ( CONSTS(77)*CONSTS(3)/CONSTS(85) ** 2.00000)/((1.00000+CONSTS(1)/ALGBRC(159) ** 3.00000+1.00000+CONSTS(3)/CONSTS(85) ** 2.00000) - 1.00000)
      ALGBRC(160) = CONSTS(90)/(1.00000+CONSTS(89)/CONSTS(91)+STATES(5)/CONSTS(92)+STATES(9)/CONSTS(93))
      ALGBRC(164) = ( CONSTS(78)*ALGBRC(160)*CONSTS(89))/(1.00000+CONSTS(87)/CONSTS(88))
      ALGBRC(158) =  CONSTS(81)*EXP(( CONSTS(83)*STATES(1)*CONSTS(6))/( 3.00000*CONSTS(4)*CONSTS(5)))
      ALGBRC(161) = ( CONSTS(73)*STATES(5)/ALGBRC(158) ** 3.00000)/((1.00000+STATES(5)/ALGBRC(158) ** 3.00000+1.00000+STATES(9)/CONSTS(84) ** 2.00000) - 1.00000)
      ALGBRC(162) = ( CONSTS(76)*CONSTS(1)/ALGBRC(159) ** 3.00000)/((1.00000+CONSTS(1)/ALGBRC(159) ** 3.00000+1.00000+CONSTS(3)/CONSTS(85) ** 2.00000) - 1.00000)
      ALGBRC(165) = ( CONSTS(80)*STATES(9)/CONSTS(84) ** 2.00000)/((1.00000+STATES(5)/ALGBRC(158) ** 3.00000+1.00000+STATES(9)/CONSTS(84) ** 2.00000) - 1.00000)
      ALGBRC(166) =  CONSTS(146)*ALGBRC(161)*CONSTS(145)+ ALGBRC(162)*ALGBRC(165)*ALGBRC(164)+ CONSTS(145)*ALGBRC(165)*ALGBRC(164)+ ALGBRC(164)*ALGBRC(161)*CONSTS(145)
      ALGBRC(167) =  ALGBRC(162)*CONSTS(144)*ALGBRC(165)+ ALGBRC(161)*CONSTS(145)*ALGBRC(163)+ ALGBRC(163)*CONSTS(144)*ALGBRC(165)+ CONSTS(145)*ALGBRC(163)*ALGBRC(165)
      ALGBRC(168) =  CONSTS(145)*ALGBRC(163)*CONSTS(146)+ ALGBRC(164)*ALGBRC(162)*CONSTS(144)+ ALGBRC(162)*CONSTS(144)*CONSTS(146)+ ALGBRC(163)*CONSTS(146)*CONSTS(144)
      ALGBRC(169) =  ALGBRC(165)*ALGBRC(164)*ALGBRC(162)+ ALGBRC(163)*CONSTS(146)*ALGBRC(161)+ ALGBRC(162)*CONSTS(146)*ALGBRC(161)+ ALGBRC(164)*ALGBRC(162)*ALGBRC(161)
      ALGBRC(170) = ALGBRC(166)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(171) = ALGBRC(167)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(174) =  3.00000*( ALGBRC(170)*ALGBRC(163) -  ALGBRC(171)*ALGBRC(164))
      ALGBRC(172) = ALGBRC(168)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(173) = ALGBRC(169)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(175) =  2.00000*( ALGBRC(173)*CONSTS(144) -  ALGBRC(172)*ALGBRC(161))
      ALGBRC(176) =  CONSTS(94)*( CONSTS(7)*ALGBRC(174)+ CONSTS(9)*ALGBRC(175))
      ALGBRC(16) = TERNRY(VOI.LE.CONSTS(14), CONSTS(13), 0.00000)
      ALGBRC(92) =  CONSTS(58)*STATES(44)*STATES(44)*(STATES(1) - ALGBRC(58))
      RATES(9) = ( - ((ALGBRC(67)+ALGBRC(69)+ALGBRC(88)+ALGBRC(90)+ALGBRC(92)+ALGBRC(95)+ALGBRC(16)) -  2.00000*ALGBRC(176))*CONSTS(33)*CONSTS(125))/( CONSTS(6)*CONSTS(131))+( ALGBRC(178)*CONSTS(130))/CONSTS(131)+- ALGBRC(181)
      ALGBRC(57) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log(CONSTS(1)/STATES(5))
      ALGBRC(61) =  CONSTS(45)*STATES(18)+ CONSTS(113)*STATES(19)
      ALGBRC(62) =  CONSTS(45)*STATES(18)+ CONSTS(113)*STATES(21)
      ALGBRC(63) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(64) =  CONSTS(46)*(STATES(1) - ALGBRC(57))*STATES(17) ** 3.00000*( (1.00000 - ALGBRC(63))*ALGBRC(61)*STATES(20)+ ALGBRC(63)*ALGBRC(62)*STATES(22))
      ALGBRC(65) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(66) =  CONSTS(48)*(STATES(1) - ALGBRC(57))*STATES(23)*( (1.00000 - ALGBRC(65))*STATES(24)+ ALGBRC(65)*STATES(25))
      ALGBRC(76) = ( 4.00000*ALGBRC(37)*( STATES(3)*EXP( 2.00000*ALGBRC(45)) -  0.341000*CONSTS(2)))/(EXP( 2.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(80) =  (1.00000 - ALGBRC(79))*CONSTS(53)*ALGBRC(76)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(116)*ALGBRC(76)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(59) =  (( CONSTS(4)*CONSTS(5))/( 2.00000*CONSTS(6)))*log(CONSTS(2)/STATES(11))
      ALGBRC(83) =  CONSTS(54)*STATES(38)*STATES(39)*(STATES(1) - ALGBRC(59))
      ALGBRC(91) =  CONSTS(57)*STATES(44)*STATES(44)*(STATES(1) - ALGBRC(57))
      ALGBRC(93) = ALGBRC(91)+ALGBRC(92)
      ALGBRC(124) = 1.00000/(1.00000+CONSTS(71)/STATES(11) ** 2.00000)
      ALGBRC(104) = 1.00000+ (CONSTS(1)/CONSTS(62))*(1.00000+1.00000/ALGBRC(97))
      ALGBRC(105) = CONSTS(1)/( CONSTS(62)*ALGBRC(97)*ALGBRC(104))
      ALGBRC(108) =  ALGBRC(105)*CONSTS(66)
      ALGBRC(98) = 1.00000+ (STATES(5)/CONSTS(62))*(1.00000+ALGBRC(97))
      ALGBRC(99) = ( STATES(5)*ALGBRC(97))/( CONSTS(62)*ALGBRC(98))
      ALGBRC(111) =  ALGBRC(99)*CONSTS(66)
      ALGBRC(101) = 1.00000+ (STATES(5)/CONSTS(60))*(1.00000+STATES(5)/CONSTS(61))
      ALGBRC(102) = ( STATES(5)*STATES(5))/( ALGBRC(101)*CONSTS(60)*CONSTS(61))
      ALGBRC(114) =  ALGBRC(102)*ALGBRC(99)*CONSTS(64)
      ALGBRC(115) =  ALGBRC(105)*CONSTS(133)*CONSTS(64)
      ALGBRC(106) = 1.00000/ALGBRC(104)
      ALGBRC(107) =  ALGBRC(106)*CONSTS(65)
      ALGBRC(109) = ALGBRC(107)+ALGBRC(108)
      ALGBRC(100) = 1.00000/ALGBRC(98)
      ALGBRC(110) = ( ALGBRC(100)*CONSTS(65))/ALGBRC(96)
      ALGBRC(112) = ALGBRC(110)+ALGBRC(111)
      ALGBRC(103) = 1.00000/ALGBRC(101)
      ALGBRC(113) =  ALGBRC(103)*STATES(11)*CONSTS(67)
      ALGBRC(116) =  CONSTS(136)*ALGBRC(112)*(ALGBRC(114)+ALGBRC(113))+ CONSTS(137)*ALGBRC(114)*(CONSTS(136)+ALGBRC(109))
      ALGBRC(117) =  CONSTS(135)*ALGBRC(114)*(ALGBRC(112)+CONSTS(137))+ ALGBRC(112)*ALGBRC(113)*(CONSTS(135)+ALGBRC(115))
      ALGBRC(118) =  CONSTS(135)*ALGBRC(109)*(ALGBRC(114)+ALGBRC(113))+ ALGBRC(115)*ALGBRC(113)*(CONSTS(136)+ALGBRC(109))
      ALGBRC(119) =  CONSTS(136)*ALGBRC(115)*(ALGBRC(112)+CONSTS(137))+ ALGBRC(109)*CONSTS(137)*(CONSTS(135)+ALGBRC(115))
      ALGBRC(120) = ALGBRC(116)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(121) = ALGBRC(117)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(122) = ALGBRC(118)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(123) = ALGBRC(119)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(125) = ( 3.00000*( ALGBRC(123)*ALGBRC(114) -  ALGBRC(120)*ALGBRC(115))+ ALGBRC(122)*ALGBRC(111)) -  ALGBRC(121)*ALGBRC(108)
      ALGBRC(126) =  ALGBRC(121)*CONSTS(136) -  ALGBRC(120)*CONSTS(135)
      ALGBRC(127) =  0.800000*CONSTS(72)*ALGBRC(124)*( CONSTS(7)*ALGBRC(125)+ CONSTS(8)*ALGBRC(126))
      ALGBRC(177) = ( CONSTS(95)*ALGBRC(37)*( STATES(5)*EXP(ALGBRC(45)) - CONSTS(1)))/(EXP(ALGBRC(45)) - 1.00000)
      ALGBRC(183) = ( CONSTS(97)*STATES(11))/(CONSTS(98)+STATES(11))
      ALGBRC(180) = ( CONSTS(96)*4.00000*ALGBRC(37)*( STATES(11)*EXP( 2.00000*ALGBRC(45)) -  0.341000*CONSTS(2)))/(EXP( 2.00000*ALGBRC(45)) - 1.00000)
      RATES(1) = - (ALGBRC(64)+ALGBRC(66)+ALGBRC(67)+ALGBRC(69)+ALGBRC(80)+ALGBRC(83)+ALGBRC(81)+ALGBRC(82)+ALGBRC(88)+ALGBRC(90)+ALGBRC(93)+ALGBRC(95)+ALGBRC(127)+ALGBRC(157)+ALGBRC(176)+ALGBRC(177)+ALGBRC(183)+ALGBRC(180)+ALGBRC(16))
      ALGBRC(182) = (STATES(5) - STATES(4))/CONSTS(100)
      RATES(4) = ( ALGBRC(182)*CONSTS(131))/CONSTS(126)
      RATES(5) = ( - (ALGBRC(64)+ALGBRC(66)+ 3.00000*ALGBRC(127)+ 3.00000*ALGBRC(176)+ALGBRC(91)+ALGBRC(177))*CONSTS(125)*CONSTS(33))/( CONSTS(6)*CONSTS(131))+( ALGBRC(179)*CONSTS(130))/CONSTS(131)+- ALGBRC(182)
      ALGBRC(184) = (STATES(3) - STATES(11))/CONSTS(99)
      ALGBRC(185) = (STATES(11) - STATES(10))/CONSTS(100)
      ALGBRC(188) = ( CONSTS(102)*ALGBRC(51))/(ALGBRC(51)+CONSTS(15))
      ALGBRC(189) = ( CONSTS(101)*ALGBRC(51))/(ALGBRC(51)+CONSTS(15))
      ALGBRC(190) = ( 0.000200000*(ALGBRC(189)+1.00000))/(1.00000+(CONSTS(103) - ALGBRC(188))/STATES(11)) - ( 0.00105000*STATES(12))/CONSTS(104)
      ALGBRC(54) = 1.00000/(1.00000+( CONSTS(29)*CONSTS(28))/CONSTS(28)+STATES(11) ** 2.00000+( CONSTS(32)*CONSTS(31))/CONSTS(31)+STATES(11) ** 2.00000)
      RATES(11) =  ALGBRC(54)*((( - ((ALGBRC(183)+ALGBRC(180)+ALGBRC(83)) -  2.00000*ALGBRC(127))*CONSTS(33)*CONSTS(125))/( 2.00000*CONSTS(6)*CONSTS(131)) - ( ALGBRC(190)*CONSTS(127))/CONSTS(131))+- ALGBRC(185)+( ALGBRC(184)*CONSTS(130))/CONSTS(131))
      ALGBRC(191) = ( 0.00260000*(ALGBRC(189)+1.00000))/(1.00000+(CONSTS(103) - ALGBRC(188))/STATES(10)) - ( 0.00420000*STATES(12))/CONSTS(104)
      ALGBRC(52) = 1.00000/(1.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(10) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(10) ** 2.00000)
      RATES(10) =  ALGBRC(52)*((( ALGBRC(185)*CONSTS(131))/CONSTS(126)+( STATES(16)*CONSTS(129))/CONSTS(126)) - ( ALGBRC(191)*CONSTS(127))/CONSTS(126))
      ALGBRC(192) = ( CONSTS(112)*CONSTS(105)*STATES(3)*(1.00000 - STATES(46)))/( (1.00000+( CONSTS(105)*CONSTS(108))/CONSTS(109))*(1.00000+( STATES(3)*CONSTS(106))/CONSTS(107)))
      ALGBRC(194) =  10.9200*ALGBRC(192)*(STATES(13) - STATES(3))
      ALGBRC(53) = 1.00000/(1.00000+( CONSTS(20)*CONSTS(21))/CONSTS(21)+STATES(3) ** 2.00000+( CONSTS(22)*CONSTS(23))/CONSTS(23)+STATES(3) ** 2.00000)
      RATES(3) =  ALGBRC(53)*((( - (ALGBRC(80) -  2.00000*ALGBRC(157))*CONSTS(33)*CONSTS(125))/( 2.00000*CONSTS(6)*CONSTS(130))+( (STATES(15)+ALGBRC(194))*CONSTS(128))/CONSTS(130)) - ALGBRC(184))
      ALGBRC(193) = ( ALGBRC(185)*CONSTS(131))/CONSTS(126)+( - ALGBRC(191)*CONSTS(127))/CONSTS(126)+( STATES(16)*CONSTS(129))/CONSTS(126)
      ALGBRC(195) = TERNRY(ALGBRC(193).GT.0.00000, ( 91.0000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000))*ALGBRC(193))/(1.00000+1.00000/STATES(14) ** 8.00000), 0.00000)
      ALGBRC(187) = ( 6.00000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000)))/(1.00000+0.0123000/STATES(14))
      RATES(16) = (ALGBRC(195) - STATES(16))/ALGBRC(187)
      ALGBRC(196) = (STATES(12) - STATES(13))/120.000
      ALGBRC(55) = 1.00000/(1.00000+( CONSTS(26)*CONSTS(25))/CONSTS(25)+STATES(13) ** 2.00000)
      RATES(13) =  ALGBRC(55)*(ALGBRC(196) - (STATES(15)+ALGBRC(194)))
      ALGBRC(198) = (STATES(12) - STATES(14))/120.000
      RATES(12) = (ALGBRC(190)+ALGBRC(191)) - (( ALGBRC(196)*CONSTS(128))/CONSTS(127)+( ALGBRC(198)*CONSTS(129))/CONSTS(127))
      ALGBRC(56) = 1.00000/(1.00000+( CONSTS(24)*CONSTS(25))/CONSTS(25)+STATES(14) ** 2.00000)
      RATES(14) =  ALGBRC(56)*(ALGBRC(198) - STATES(16))
      ALGBRC(197) = - (( ALGBRC(80)*CONSTS(125))/( CONSTS(130)*CONSTS(6)*2.00000)+( - (STATES(15)+ALGBRC(194))*CONSTS(128))/CONSTS(130)+ALGBRC(184))
      ALGBRC(199) = TERNRY(ALGBRC(197).GT.0.00000, ( 15.0000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000))*ALGBRC(197))/(1.00000+1.00000/STATES(13) ** 8.00000), 0.00000)
      ALGBRC(186) = ( 2.00000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000)))/(1.00000+0.0123000/STATES(13))
      RATES(15) = (ALGBRC(199) - STATES(15))/ALGBRC(186)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(35))/CONSTS(36)))
      ALGBRC(17) = 1.00000/( CONSTS(39)*EXP((STATES(1)+CONSTS(37))/CONSTS(38))+ CONSTS(40)*EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(43))/CONSTS(44)))
      ALGBRC(18) = 1.00000/( 3.68600e-06*EXP(- (STATES(1)+3.88750)/7.85790)+ 16.0000*EXP((STATES(1) - 0.496300)/9.18430))
      ALGBRC(19) = 1.00000/( 0.00979400*EXP(- (STATES(1)+17.9500)/28.0500)+ 0.334300*EXP((STATES(1)+5.73000)/56.6600))
      ALGBRC(5) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/13.0000))
      ALGBRC(21) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- (STATES(1) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814)))
      ALGBRC(6) = 1.00000/(1.00000+EXP((27.0000+STATES(1))/13.0000))
      ALGBRC(22) = 43.0000+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
      ALGBRC(23) = 6.16200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1) - 8.00000)/8.59000))
      ALGBRC(7) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000+3.30000)/4.23000))
      ALGBRC(26) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+19.5800+3.30000)/3.69600))
      ALGBRC(27) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000+15.1900)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000+15.1900)/10.0000))
      ALGBRC(28) = 1000.00+1.00000/( 3.50000e-05*EXP(- (STATES(1)+5.00000+15.1900)/4.00000)+ 3.50000e-05*EXP((STATES(1)+5.00000+15.1900)/6.00000))
      ALGBRC(24) = ALGBRC(8)
      ALGBRC(9) =  STATES(34)*1.00000
      ALGBRC(25) = 1.00000/(CONSTS(52)/ALGBRC(9)+1.00000+CONSTS(51)/STATES(3) ** 4.00000)
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+30.0000)/7.00000))
      ALGBRC(29) = 1.00000/( 1.06800*EXP((STATES(1)+16.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+16.3000)/30.0000))
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+61.0000)/5.00000))
      ALGBRC(30) = 1.00000/( 0.0150000*EXP((STATES(1)+71.7000)/15.4000)+ 0.0150000*EXP(- (STATES(1)+71.7000)/83.3000))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+8.33700)/6.78900))
      ALGBRC(31) = 12.9800+1.00000/( 0.365200*EXP(((STATES(1)+17.6000) - 31.6600)/3.86900)+ 4.12300e-05*EXP(- ((STATES(1)+17.6000) - 47.7800)/20.3800))
      ALGBRC(32) = 1.86500+1.00000/( 0.0662900*EXP(((STATES(1)+17.2000) - 34.7000)/7.35500)+ 1.12800e-05*EXP(- ((STATES(1)+17.2000) - 29.7400)/25.9400))
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(34) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+87.0000)/9.50000))
      ALGBRC(35) = 2000.00/(EXP((STATES(1)+57.0000)/60.0000)+EXP(- (STATES(1)+132.000)/10.0000))
      ALGBRC(15) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(3)+144.590)/( 1.56920*CONSTS(3)+3.81150)))
      ALGBRC(36) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
      ALGBRC(20) = ALGBRC(2)
      ALGBRC(38) = 4.85900+1.00000/( 0.862800*EXP(- (STATES(1)+116.726)/7.60050)+ 1.10960*EXP((STATES(1)+6.27190)/9.03580))
      ALGBRC(41) =  0.720000*(7.00000+1.00000/( 0.0400000*EXP(- ((STATES(1)+15.1900) - 4.00000)/7.00000)+ 0.0400000*EXP(((STATES(1)+15.1900) - 4.00000)/7.00000)))
      ALGBRC(42) =  0.490000*(100.000+1.00000/( 0.000120000*EXP(- (STATES(1)+15.1900)/3.00000)+ 0.000120000*EXP((STATES(1)+15.1900)/7.00000)))
      ALGBRC(43) =  2.50000*ALGBRC(27)
      ALGBRC(33) = ALGBRC(13)
      ALGBRC(44) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      ALGBRC(50) = ( CONSTS(18)*(1.00000 - STATES(2)))/(1.00000+CONSTS(19)/STATES(3))
      ALGBRC(39) = 1.00000/(1.00000+EXP((STATES(1)+84.7000)/6.22000))
      ALGBRC(46) =  3.00000*ALGBRC(19)
      ALGBRC(47) =  1.46000*ALGBRC(38)
      ALGBRC(40) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(48) = ALGBRC(17)
      ALGBRC(49) =  2.50000*ALGBRC(41)
      ALGBRC(70) =  CONSTS(115)*STATES(30)+ CONSTS(121)*STATES(31)
      ALGBRC(71) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(72) = 1.00000 - ALGBRC(71)
      ALGBRC(73) =  ALGBRC(71)*STATES(32)+ ALGBRC(72)*STATES(33)
      ALGBRC(74) =  CONSTS(115)*STATES(35)+ CONSTS(121)*STATES(31)
      ALGBRC(75) =  ALGBRC(71)*STATES(36)+ ALGBRC(72)*STATES(33)
      ALGBRC(37) = ( STATES(1)*CONSTS(6)*CONSTS(6))/( CONSTS(4)*CONSTS(5))
      ALGBRC(45) = ( STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5))
      ALGBRC(78) = ( 1.00000*ALGBRC(37)*( 0.750000*STATES(8)*EXP( 1.00000*ALGBRC(45)) -  0.750000*CONSTS(3)))/(EXP( 1.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(51) = ALGBRC(50)+STATES(2)
      ALGBRC(79) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(82) =  (1.00000 - ALGBRC(79))*CONSTS(118)*ALGBRC(78)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(123)*ALGBRC(78)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(178) = (STATES(8) - STATES(9))/CONSTS(99)
      ALGBRC(77) = ( 1.00000*ALGBRC(37)*( 0.750000*STATES(6)*EXP( 1.00000*ALGBRC(45)) -  0.750000*CONSTS(1)))/(EXP( 1.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(81) =  (1.00000 - ALGBRC(79))*CONSTS(117)*ALGBRC(77)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(122)*ALGBRC(77)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(154) = 1.00000/(1.00000+CONSTS(71)/STATES(3) ** 2.00000)
      ALGBRC(97) = EXP(( CONSTS(69)*STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5)))
      ALGBRC(134) = 1.00000+ (CONSTS(1)/CONSTS(62))*(1.00000+1.00000/ALGBRC(97))
      ALGBRC(135) = CONSTS(1)/( CONSTS(62)*ALGBRC(97)*ALGBRC(134))
      ALGBRC(138) =  ALGBRC(135)*CONSTS(66)
      ALGBRC(128) = 1.00000+ (STATES(6)/CONSTS(62))*(1.00000+ALGBRC(97))
      ALGBRC(129) = ( STATES(6)*ALGBRC(97))/( CONSTS(62)*ALGBRC(128))
      ALGBRC(141) =  ALGBRC(129)*CONSTS(66)
      ALGBRC(131) = 1.00000+ (STATES(6)/CONSTS(60))*(1.00000+STATES(6)/CONSTS(61))
      ALGBRC(132) = ( STATES(6)*STATES(6))/( ALGBRC(131)*CONSTS(60)*CONSTS(61))
      ALGBRC(144) =  ALGBRC(132)*ALGBRC(129)*CONSTS(64)
      ALGBRC(145) =  ALGBRC(135)*CONSTS(139)*CONSTS(64)
      ALGBRC(136) = 1.00000/ALGBRC(134)
      ALGBRC(137) =  ALGBRC(136)*CONSTS(65)
      ALGBRC(139) = ALGBRC(137)+ALGBRC(138)
      ALGBRC(96) = EXP(( CONSTS(70)*STATES(1)*CONSTS(6))/( CONSTS(4)*CONSTS(5)))
      ALGBRC(130) = 1.00000/ALGBRC(128)
      ALGBRC(140) = ( ALGBRC(130)*CONSTS(65))/ALGBRC(96)
      ALGBRC(142) = ALGBRC(140)+ALGBRC(141)
      ALGBRC(133) = 1.00000/ALGBRC(131)
      ALGBRC(143) =  ALGBRC(133)*STATES(3)*CONSTS(67)
      ALGBRC(146) =  CONSTS(142)*ALGBRC(142)*(ALGBRC(144)+ALGBRC(143))+ CONSTS(143)*ALGBRC(144)*(CONSTS(142)+ALGBRC(139))
      ALGBRC(147) =  CONSTS(141)*ALGBRC(144)*(ALGBRC(142)+CONSTS(143))+ ALGBRC(142)*ALGBRC(143)*(CONSTS(141)+ALGBRC(145))
      ALGBRC(148) =  CONSTS(141)*ALGBRC(139)*(ALGBRC(144)+ALGBRC(143))+ ALGBRC(145)*ALGBRC(143)*(CONSTS(142)+ALGBRC(139))
      ALGBRC(149) =  CONSTS(142)*ALGBRC(145)*(ALGBRC(142)+CONSTS(143))+ ALGBRC(139)*CONSTS(143)*(CONSTS(141)+ALGBRC(145))
      ALGBRC(150) = ALGBRC(146)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(151) = ALGBRC(147)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(152) = ALGBRC(148)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(153) = ALGBRC(149)/(ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(149))
      ALGBRC(155) = ( 3.00000*( ALGBRC(153)*ALGBRC(144) -  ALGBRC(150)*ALGBRC(145))+ ALGBRC(152)*ALGBRC(141)) -  ALGBRC(151)*ALGBRC(138)
      ALGBRC(156) =  ALGBRC(151)*CONSTS(142) -  ALGBRC(150)*CONSTS(141)
      ALGBRC(157) =  0.200000*CONSTS(72)*ALGBRC(154)*( CONSTS(7)*ALGBRC(155)+ CONSTS(8)*ALGBRC(156))
      ALGBRC(179) = (STATES(6) - STATES(5))/CONSTS(99)
      ALGBRC(181) = (STATES(9) - STATES(7))/CONSTS(100)
      ALGBRC(58) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log(CONSTS(3)/STATES(9))
      ALGBRC(67) =  CONSTS(49)*STATES(26)*STATES(27)*STATES(28)*(STATES(1) - ALGBRC(58))
      ALGBRC(68) = 1.00000/(1.00000+EXP(- (STATES(1) - 12.0000)/16.0000))
      ALGBRC(69) =  CONSTS(50)*ALGBRC(68)*(STATES(1) - ALGBRC(58))
      ALGBRC(84) = 1.00000/(1.00000+EXP((STATES(1)+54.8100)/38.2100))
      ALGBRC(85) = 1.00000 - ALGBRC(84)
      ALGBRC(86) =  ALGBRC(84)*STATES(40)+ ALGBRC(85)*STATES(41)
      ALGBRC(87) = ( (1.00000/(1.00000+EXP((STATES(1)+55.0000)/( 0.320000*75.0000))))*1.00000)/(1.00000+EXP((STATES(1) - 10.0000)/( 0.320000*30.0000)))
      ALGBRC(88) =  CONSTS(55)* (CONSTS(3)/5.40000) ** (1.0 / 2)*ALGBRC(86)*ALGBRC(87)*(STATES(1) - ALGBRC(58))
      ALGBRC(60) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log((CONSTS(3)+ CONSTS(34)*CONSTS(1))/(STATES(9)+ CONSTS(34)*STATES(5)))
      ALGBRC(89) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(11) ** 1.40000)
      ALGBRC(90) =  CONSTS(56)*ALGBRC(89)*STATES(42)*STATES(43)*(STATES(1) - ALGBRC(60))
      ALGBRC(94) = 1.00000/(1.00000+EXP(((STATES(1)+116.000) -  5.50000*CONSTS(3))/11.0000))
      ALGBRC(95) =  CONSTS(59)*2.32380* (CONSTS(3)/5.40000) ** (1.0 / 2)*ALGBRC(94)*STATES(45)*(STATES(1) - ALGBRC(58))
      ALGBRC(159) =  CONSTS(82)*EXP(( (1.00000 - CONSTS(83))*STATES(1)*CONSTS(6))/( 3.00000*CONSTS(4)*CONSTS(5)))
      ALGBRC(163) = ( CONSTS(77)*CONSTS(3)/CONSTS(85) ** 2.00000)/((1.00000+CONSTS(1)/ALGBRC(159) ** 3.00000+1.00000+CONSTS(3)/CONSTS(85) ** 2.00000) - 1.00000)
      ALGBRC(160) = CONSTS(90)/(1.00000+CONSTS(89)/CONSTS(91)+STATES(5)/CONSTS(92)+STATES(9)/CONSTS(93))
      ALGBRC(164) = ( CONSTS(78)*ALGBRC(160)*CONSTS(89))/(1.00000+CONSTS(87)/CONSTS(88))
      ALGBRC(158) =  CONSTS(81)*EXP(( CONSTS(83)*STATES(1)*CONSTS(6))/( 3.00000*CONSTS(4)*CONSTS(5)))
      ALGBRC(161) = ( CONSTS(73)*STATES(5)/ALGBRC(158) ** 3.00000)/((1.00000+STATES(5)/ALGBRC(158) ** 3.00000+1.00000+STATES(9)/CONSTS(84) ** 2.00000) - 1.00000)
      ALGBRC(162) = ( CONSTS(76)*CONSTS(1)/ALGBRC(159) ** 3.00000)/((1.00000+CONSTS(1)/ALGBRC(159) ** 3.00000+1.00000+CONSTS(3)/CONSTS(85) ** 2.00000) - 1.00000)
      ALGBRC(165) = ( CONSTS(80)*STATES(9)/CONSTS(84) ** 2.00000)/((1.00000+STATES(5)/ALGBRC(158) ** 3.00000+1.00000+STATES(9)/CONSTS(84) ** 2.00000) - 1.00000)
      ALGBRC(166) =  CONSTS(146)*ALGBRC(161)*CONSTS(145)+ ALGBRC(162)*ALGBRC(165)*ALGBRC(164)+ CONSTS(145)*ALGBRC(165)*ALGBRC(164)+ ALGBRC(164)*ALGBRC(161)*CONSTS(145)
      ALGBRC(167) =  ALGBRC(162)*CONSTS(144)*ALGBRC(165)+ ALGBRC(161)*CONSTS(145)*ALGBRC(163)+ ALGBRC(163)*CONSTS(144)*ALGBRC(165)+ CONSTS(145)*ALGBRC(163)*ALGBRC(165)
      ALGBRC(168) =  CONSTS(145)*ALGBRC(163)*CONSTS(146)+ ALGBRC(164)*ALGBRC(162)*CONSTS(144)+ ALGBRC(162)*CONSTS(144)*CONSTS(146)+ ALGBRC(163)*CONSTS(146)*CONSTS(144)
      ALGBRC(169) =  ALGBRC(165)*ALGBRC(164)*ALGBRC(162)+ ALGBRC(163)*CONSTS(146)*ALGBRC(161)+ ALGBRC(162)*CONSTS(146)*ALGBRC(161)+ ALGBRC(164)*ALGBRC(162)*ALGBRC(161)
      ALGBRC(170) = ALGBRC(166)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(171) = ALGBRC(167)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(174) =  3.00000*( ALGBRC(170)*ALGBRC(163) -  ALGBRC(171)*ALGBRC(164))
      ALGBRC(172) = ALGBRC(168)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(173) = ALGBRC(169)/(ALGBRC(166)+ALGBRC(167)+ALGBRC(168)+ALGBRC(169))
      ALGBRC(175) =  2.00000*( ALGBRC(173)*CONSTS(144) -  ALGBRC(172)*ALGBRC(161))
      ALGBRC(176) =  CONSTS(94)*( CONSTS(7)*ALGBRC(174)+ CONSTS(9)*ALGBRC(175))
      ALGBRC(16) = TERNRY(VOI.LE.CONSTS(14), CONSTS(13), 0.00000)
      ALGBRC(92) =  CONSTS(58)*STATES(44)*STATES(44)*(STATES(1) - ALGBRC(58))
      ALGBRC(57) =  (( CONSTS(4)*CONSTS(5))/CONSTS(6))*log(CONSTS(1)/STATES(5))
      ALGBRC(61) =  CONSTS(45)*STATES(18)+ CONSTS(113)*STATES(19)
      ALGBRC(62) =  CONSTS(45)*STATES(18)+ CONSTS(113)*STATES(21)
      ALGBRC(63) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(64) =  CONSTS(46)*(STATES(1) - ALGBRC(57))*STATES(17) ** 3.00000*( (1.00000 - ALGBRC(63))*ALGBRC(61)*STATES(20)+ ALGBRC(63)*ALGBRC(62)*STATES(22))
      ALGBRC(65) = 1.00000/(1.00000+CONSTS(15)/ALGBRC(51))
      ALGBRC(66) =  CONSTS(48)*(STATES(1) - ALGBRC(57))*STATES(23)*( (1.00000 - ALGBRC(65))*STATES(24)+ ALGBRC(65)*STATES(25))
      ALGBRC(76) = ( 4.00000*ALGBRC(37)*( STATES(3)*EXP( 2.00000*ALGBRC(45)) -  0.341000*CONSTS(2)))/(EXP( 2.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(80) =  (1.00000 - ALGBRC(79))*CONSTS(53)*ALGBRC(76)*STATES(29)*( ALGBRC(70)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(73)*STATES(37))+ ALGBRC(79)*CONSTS(116)*ALGBRC(76)*STATES(29)*( ALGBRC(74)*(1.00000 - STATES(37))+ STATES(34)*ALGBRC(75)*STATES(37))
      ALGBRC(59) =  (( CONSTS(4)*CONSTS(5))/( 2.00000*CONSTS(6)))*log(CONSTS(2)/STATES(11))
      ALGBRC(83) =  CONSTS(54)*STATES(38)*STATES(39)*(STATES(1) - ALGBRC(59))
      ALGBRC(91) =  CONSTS(57)*STATES(44)*STATES(44)*(STATES(1) - ALGBRC(57))
      ALGBRC(93) = ALGBRC(91)+ALGBRC(92)
      ALGBRC(124) = 1.00000/(1.00000+CONSTS(71)/STATES(11) ** 2.00000)
      ALGBRC(104) = 1.00000+ (CONSTS(1)/CONSTS(62))*(1.00000+1.00000/ALGBRC(97))
      ALGBRC(105) = CONSTS(1)/( CONSTS(62)*ALGBRC(97)*ALGBRC(104))
      ALGBRC(108) =  ALGBRC(105)*CONSTS(66)
      ALGBRC(98) = 1.00000+ (STATES(5)/CONSTS(62))*(1.00000+ALGBRC(97))
      ALGBRC(99) = ( STATES(5)*ALGBRC(97))/( CONSTS(62)*ALGBRC(98))
      ALGBRC(111) =  ALGBRC(99)*CONSTS(66)
      ALGBRC(101) = 1.00000+ (STATES(5)/CONSTS(60))*(1.00000+STATES(5)/CONSTS(61))
      ALGBRC(102) = ( STATES(5)*STATES(5))/( ALGBRC(101)*CONSTS(60)*CONSTS(61))
      ALGBRC(114) =  ALGBRC(102)*ALGBRC(99)*CONSTS(64)
      ALGBRC(115) =  ALGBRC(105)*CONSTS(133)*CONSTS(64)
      ALGBRC(106) = 1.00000/ALGBRC(104)
      ALGBRC(107) =  ALGBRC(106)*CONSTS(65)
      ALGBRC(109) = ALGBRC(107)+ALGBRC(108)
      ALGBRC(100) = 1.00000/ALGBRC(98)
      ALGBRC(110) = ( ALGBRC(100)*CONSTS(65))/ALGBRC(96)
      ALGBRC(112) = ALGBRC(110)+ALGBRC(111)
      ALGBRC(103) = 1.00000/ALGBRC(101)
      ALGBRC(113) =  ALGBRC(103)*STATES(11)*CONSTS(67)
      ALGBRC(116) =  CONSTS(136)*ALGBRC(112)*(ALGBRC(114)+ALGBRC(113))+ CONSTS(137)*ALGBRC(114)*(CONSTS(136)+ALGBRC(109))
      ALGBRC(117) =  CONSTS(135)*ALGBRC(114)*(ALGBRC(112)+CONSTS(137))+ ALGBRC(112)*ALGBRC(113)*(CONSTS(135)+ALGBRC(115))
      ALGBRC(118) =  CONSTS(135)*ALGBRC(109)*(ALGBRC(114)+ALGBRC(113))+ ALGBRC(115)*ALGBRC(113)*(CONSTS(136)+ALGBRC(109))
      ALGBRC(119) =  CONSTS(136)*ALGBRC(115)*(ALGBRC(112)+CONSTS(137))+ ALGBRC(109)*CONSTS(137)*(CONSTS(135)+ALGBRC(115))
      ALGBRC(120) = ALGBRC(116)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(121) = ALGBRC(117)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(122) = ALGBRC(118)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(123) = ALGBRC(119)/(ALGBRC(116)+ALGBRC(117)+ALGBRC(118)+ALGBRC(119))
      ALGBRC(125) = ( 3.00000*( ALGBRC(123)*ALGBRC(114) -  ALGBRC(120)*ALGBRC(115))+ ALGBRC(122)*ALGBRC(111)) -  ALGBRC(121)*ALGBRC(108)
      ALGBRC(126) =  ALGBRC(121)*CONSTS(136) -  ALGBRC(120)*CONSTS(135)
      ALGBRC(127) =  0.800000*CONSTS(72)*ALGBRC(124)*( CONSTS(7)*ALGBRC(125)+ CONSTS(8)*ALGBRC(126))
      ALGBRC(177) = ( CONSTS(95)*ALGBRC(37)*( STATES(5)*EXP(ALGBRC(45)) - CONSTS(1)))/(EXP(ALGBRC(45)) - 1.00000)
      ALGBRC(183) = ( CONSTS(97)*STATES(11))/(CONSTS(98)+STATES(11))
      ALGBRC(180) = ( CONSTS(96)*4.00000*ALGBRC(37)*( STATES(11)*EXP( 2.00000*ALGBRC(45)) -  0.341000*CONSTS(2)))/(EXP( 2.00000*ALGBRC(45)) - 1.00000)
      ALGBRC(182) = (STATES(5) - STATES(4))/CONSTS(100)
      ALGBRC(184) = (STATES(3) - STATES(11))/CONSTS(99)
      ALGBRC(185) = (STATES(11) - STATES(10))/CONSTS(100)
      ALGBRC(188) = ( CONSTS(102)*ALGBRC(51))/(ALGBRC(51)+CONSTS(15))
      ALGBRC(189) = ( CONSTS(101)*ALGBRC(51))/(ALGBRC(51)+CONSTS(15))
      ALGBRC(190) = ( 0.000200000*(ALGBRC(189)+1.00000))/(1.00000+(CONSTS(103) - ALGBRC(188))/STATES(11)) - ( 0.00105000*STATES(12))/CONSTS(104)
      ALGBRC(54) = 1.00000/(1.00000+( CONSTS(29)*CONSTS(28))/CONSTS(28)+STATES(11) ** 2.00000+( CONSTS(32)*CONSTS(31))/CONSTS(31)+STATES(11) ** 2.00000)
      ALGBRC(191) = ( 0.00260000*(ALGBRC(189)+1.00000))/(1.00000+(CONSTS(103) - ALGBRC(188))/STATES(10)) - ( 0.00420000*STATES(12))/CONSTS(104)
      ALGBRC(52) = 1.00000/(1.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(10) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(10) ** 2.00000)
      ALGBRC(192) = ( CONSTS(112)*CONSTS(105)*STATES(3)*(1.00000 - STATES(46)))/( (1.00000+( CONSTS(105)*CONSTS(108))/CONSTS(109))*(1.00000+( STATES(3)*CONSTS(106))/CONSTS(107)))
      ALGBRC(194) =  10.9200*ALGBRC(192)*(STATES(13) - STATES(3))
      ALGBRC(53) = 1.00000/(1.00000+( CONSTS(20)*CONSTS(21))/CONSTS(21)+STATES(3) ** 2.00000+( CONSTS(22)*CONSTS(23))/CONSTS(23)+STATES(3) ** 2.00000)
      ALGBRC(193) = ( ALGBRC(185)*CONSTS(131))/CONSTS(126)+( - ALGBRC(191)*CONSTS(127))/CONSTS(126)+( STATES(16)*CONSTS(129))/CONSTS(126)
      ALGBRC(195) = TERNRY(ALGBRC(193).GT.0.00000, ( 91.0000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000))*ALGBRC(193))/(1.00000+1.00000/STATES(14) ** 8.00000), 0.00000)
      ALGBRC(187) = ( 6.00000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000)))/(1.00000+0.0123000/STATES(14))
      ALGBRC(196) = (STATES(12) - STATES(13))/120.000
      ALGBRC(55) = 1.00000/(1.00000+( CONSTS(26)*CONSTS(25))/CONSTS(25)+STATES(13) ** 2.00000)
      ALGBRC(198) = (STATES(12) - STATES(14))/120.000
      ALGBRC(56) = 1.00000/(1.00000+( CONSTS(24)*CONSTS(25))/CONSTS(25)+STATES(14) ** 2.00000)
      ALGBRC(197) = - (( ALGBRC(80)*CONSTS(125))/( CONSTS(130)*CONSTS(6)*2.00000)+( - (STATES(15)+ALGBRC(194))*CONSTS(128))/CONSTS(130)+ALGBRC(184))
      ALGBRC(199) = TERNRY(ALGBRC(197).GT.0.00000, ( 15.0000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000))*ALGBRC(197))/(1.00000+1.00000/STATES(13) ** 8.00000), 0.00000)
      ALGBRC(186) = ( 2.00000*(1.00000+( 1.00000*1.00000)/(1.00000+0.280000/ALGBRC(51) ** 8.00000)))/(1.00000+0.0123000/STATES(13))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END