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 200 entries in the algebraic variable array.
C There are a total of 49 entries in each of the rate and state variable arrays.
C There are a total of 206 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 R in component physical_constants (joule_per_kilomole_kelvin).
C CONSTS(6) is T in component physical_constants (kelvin).
C CONSTS(7) is F in component physical_constants (coulomb_per_mole).
C CONSTS(8) is zna in component physical_constants (dimensionless).
C CONSTS(9) is zca in component physical_constants (dimensionless).
C CONSTS(10) is zk in component physical_constants (dimensionless).
C CONSTS(11) is L in component cell_geometry (centimeter).
C CONSTS(12) is rad in component cell_geometry (centimeter).
C CONSTS(163) is vcell in component cell_geometry (microliter).
C CONSTS(178) is Ageo in component cell_geometry (centimeter_squared).
C CONSTS(184) is Acap in component cell_geometry (centimeter_squared).
C CONSTS(185) is vmyo in component cell_geometry (microliter).
C CONSTS(186) is vnsr in component cell_geometry (microliter).
C CONSTS(187) is vjsr in component cell_geometry (microliter).
C CONSTS(188) is vss in component cell_geometry (microliter).
C STATES(1) is v in component membrane (millivolt).
C ALGBRC(13) is vfrt in component membrane (dimensionless).
C CONSTS(170) is ffrt in component membrane (coulomb_per_mole_millivolt).
C CONSTS(150) is frt in component membrane (per_millivolt).
C ALGBRC(59) is INa in component INa (microA_per_microF).
C ALGBRC(61) is INaL in component INaL (microA_per_microF).
C ALGBRC(67) is Ito in component Ito (microA_per_microF).
C ALGBRC(84) is ICaL in component ICaL (microA_per_microF).
C ALGBRC(85) is ICaNa in component ICaL (microA_per_microF).
C ALGBRC(88) is ICaK in component ICaL (microA_per_microF).
C ALGBRC(91) is IKr in component IKr (microA_per_microF).
C ALGBRC(97) is IKs in component IKs (microA_per_microF).
C ALGBRC(99) is IK1 in component IK1 (microA_per_microF).
C ALGBRC(131) is INaCa_i in component INaCa_i (microA_per_microF).
C ALGBRC(161) is INaCa_ss in component INaCa_i (microA_per_microF).
C ALGBRC(180) is INaK in component INaK (microA_per_microF).
C ALGBRC(186) is INab in component INab (microA_per_microF).
C ALGBRC(182) is IKb in component IKb (microA_per_microF).
C ALGBRC(191) is IpCa in component IpCa (microA_per_microF).
C ALGBRC(190) is ICab in component ICab (microA_per_microF).
C ALGBRC(1) is Istim in component membrane (microA_per_microF).
C CONSTS(13) is i_Stim_Start in component membrane (millisecond).
C CONSTS(14) is i_Stim_End in component membrane (millisecond).
C CONSTS(15) is i_Stim_Amplitude in component membrane (microA_per_microF).
C CONSTS(16) is i_Stim_Period in component membrane (millisecond).
C CONSTS(17) is i_Stim_PulseDuration in component membrane (millisecond).
C CONSTS(18) is KmCaMK in component CaMK (millimolar).
C CONSTS(19) is aCaMK in component CaMK (per_millimolar_per_millisecond).
C CONSTS(20) is bCaMK in component CaMK (per_millisecond).
C CONSTS(21) is CaMKo in component CaMK (dimensionless).
C CONSTS(22) is KmCaM in component CaMK (millimolar).
C ALGBRC(37) is CaMKb in component CaMK (millimolar).
C ALGBRC(43) 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(23) is cmdnmax_b in component intracellular_ions (millimolar).
C CONSTS(151) is cmdnmax in component intracellular_ions (millimolar).
C CONSTS(24) is kmcmdn in component intracellular_ions (millimolar).
C CONSTS(25) is trpnmax in component intracellular_ions (millimolar).
C CONSTS(26) is kmtrpn in component intracellular_ions (millimolar).
C CONSTS(27) is BSRmax in component intracellular_ions (millimolar).
C CONSTS(28) is KmBSR in component intracellular_ions (millimolar).
C CONSTS(29) is BSLmax in component intracellular_ions (millimolar).
C CONSTS(30) is KmBSL in component intracellular_ions (millimolar).
C CONSTS(31) is csqnmax in component intracellular_ions (millimolar).
C CONSTS(32) 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 ALGBRC(188) is JdiffNa in component diff (millimolar_per_millisecond).
C ALGBRC(192) is Jdiff in component diff (millimolar_per_millisecond).
C ALGBRC(199) is Jup in component SERCA (millimolar_per_millisecond).
C ALGBRC(184) is JdiffK in component diff (millimolar_per_millisecond).
C ALGBRC(194) is Jrel in component ryr (millimolar_per_millisecond).
C ALGBRC(200) is Jtr in component trans_flux (millimolar_per_millisecond).
C ALGBRC(45) is Bcai in component intracellular_ions (dimensionless).
C ALGBRC(49) is Bcajsr in component intracellular_ions (dimensionless).
C ALGBRC(47) is Bcass 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(51) is ENa in component reversal_potentials (millivolt).
C ALGBRC(54) is EK in component reversal_potentials (millivolt).
C ALGBRC(55) is EKs in component reversal_potentials (millivolt).
C ALGBRC(2) is mss in component INa (dimensionless).
C ALGBRC(14) 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(11) is m in component INa (dimensionless).
C ALGBRC(3) is hss in component INa (dimensionless).
C ALGBRC(15) is thf in component INa (millisecond).
C ALGBRC(16) 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(152) is Ahs in component INa (dimensionless).
C CONSTS(45) is Ahf in component INa (dimensionless).
C STATES(12) is hf in component INa (dimensionless).
C STATES(13) is hs in component INa (dimensionless).
C ALGBRC(56) is h in component INa (dimensionless).
C CONSTS(46) is GNa in component INa (milliS_per_microF).
C CONSTS(47) is shift_INa_inact in component INa (millivolt).
C ALGBRC(17) is jss in component INa (dimensionless).
C ALGBRC(28) is tj in component INa (millisecond).
C STATES(14) is j in component INa (dimensionless).
C ALGBRC(29) is hssp in component INa (dimensionless).
C ALGBRC(38) is thsp in component INa (millisecond).
C STATES(15) is hsp in component INa (dimensionless).
C ALGBRC(57) is hp in component INa (dimensionless).
C ALGBRC(39) is tjp in component INa (millisecond).
C STATES(16) is jp in component INa (dimensionless).
C ALGBRC(58) is fINap in component INa (dimensionless).
C ALGBRC(30) is mLss in component INaL (dimensionless).
C ALGBRC(40) is tmL in component INaL (millisecond).
C STATES(17) is mL in component INaL (dimensionless).
C CONSTS(48) is thL in component INaL (millisecond).
C ALGBRC(4) is hLss in component INaL (dimensionless).
C STATES(18) is hL in component INaL (dimensionless).
C ALGBRC(5) is hLssp in component INaL (dimensionless).
C CONSTS(153) is thLp in component INaL (millisecond).
C STATES(19) is hLp in component INaL (dimensionless).
C CONSTS(49) is GNaL_b in component INaL (milliS_per_microF).
C CONSTS(154) is GNaL in component INaL (milliS_per_microF).
C ALGBRC(60) is fINaLp in component INaL (dimensionless).
C CONSTS(50) 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(20) is a in component Ito (dimensionless).
C ALGBRC(7) is iss in component Ito (dimensionless).
C ALGBRC(19) is delta_epi in component Ito (dimensionless).
C ALGBRC(31) is tiF_b in component Ito (millisecond).
C ALGBRC(41) is tiS_b in component Ito (millisecond).
C ALGBRC(44) is tiF in component Ito (millisecond).
C ALGBRC(46) is tiS in component Ito (millisecond).
C ALGBRC(62) is AiF in component Ito (dimensionless).
C ALGBRC(63) is AiS in component Ito (dimensionless).
C STATES(21) is iF in component Ito (dimensionless).
C STATES(22) is iS in component Ito (dimensionless).
C ALGBRC(64) is i in component Ito (dimensionless).
C ALGBRC(32) is assp in component Ito (dimensionless).
C STATES(23) is ap in component Ito (dimensionless).
C ALGBRC(48) is dti_develop in component Ito (dimensionless).
C ALGBRC(50) is dti_recover in component Ito (dimensionless).
C ALGBRC(52) is tiFp in component Ito (millisecond).
C ALGBRC(53) is tiSp in component Ito (millisecond).
C STATES(24) is iFp in component Ito (dimensionless).
C STATES(25) is iSp in component Ito (dimensionless).
C ALGBRC(65) is ip in component Ito (dimensionless).
C CONSTS(155) is Gto in component Ito (milliS_per_microF).
C ALGBRC(66) is fItop in component Ito (dimensionless).
C CONSTS(51) is Kmn in component ICaL (millimolar).
C CONSTS(52) is k2n in component ICaL (per_millisecond).
C CONSTS(53) is PCa_b in component ICaL (dimensionless).
C ALGBRC(8) is dss in component ICaL (dimensionless).
C STATES(26) is d in component ICaL (dimensionless).
C ALGBRC(9) is fss in component ICaL (dimensionless).
C CONSTS(156) is Aff in component ICaL (dimensionless).
C CONSTS(171) is Afs in component ICaL (dimensionless).
C STATES(27) is ff in component ICaL (dimensionless).
C STATES(28) is fs in component ICaL (dimensionless).
C ALGBRC(68) is f in component ICaL (dimensionless).
C ALGBRC(20) is fcass in component ICaL (dimensionless).
C ALGBRC(69) is Afcaf in component ICaL (dimensionless).
C ALGBRC(70) is Afcas in component ICaL (dimensionless).
C STATES(29) is fcaf in component ICaL (dimensionless).
C STATES(30) is fcas in component ICaL (dimensionless).
C ALGBRC(71) is fca in component ICaL (dimensionless).
C STATES(31) is jca in component ICaL (dimensionless).
C STATES(32) is ffp in component ICaL (dimensionless).
C ALGBRC(72) is fp in component ICaL (dimensionless).
C STATES(33) is fcafp in component ICaL (dimensionless).
C ALGBRC(73) is fcap in component ICaL (dimensionless).
C ALGBRC(10) is km2n in component ICaL (per_millisecond).
C ALGBRC(21) is anca in component ICaL (dimensionless).
C STATES(34) is nca in component ICaL (dimensionless).
C ALGBRC(76) is PhiCaL in component ICaL (dimensionless).
C ALGBRC(79) is PhiCaNa in component ICaL (dimensionless).
C ALGBRC(82) is PhiCaK in component ICaL (dimensionless).
C CONSTS(157) is PCa in component ICaL (dimensionless).
C CONSTS(172) is PCap in component ICaL (dimensionless).
C CONSTS(173) is PCaNa in component ICaL (dimensionless).
C CONSTS(174) is PCaK in component ICaL (dimensionless).
C CONSTS(182) is PCaNap in component ICaL (dimensionless).
C CONSTS(183) is PCaKp in component ICaL (dimensionless).
C ALGBRC(83) is fICaLp in component ICaL (dimensionless).
C ALGBRC(22) is td in component ICaL (millisecond).
C ALGBRC(23) is tff in component ICaL (millisecond).
C ALGBRC(24) is tfs in component ICaL (millisecond).
C ALGBRC(33) is tfcaf in component ICaL (millisecond).
C ALGBRC(34) is tfcas in component ICaL (millisecond).
C CONSTS(158) is tjca in component ICaL (millisecond).
C ALGBRC(35) is tffp in component ICaL (millisecond).
C ALGBRC(42) is tfcafp in component ICaL (millisecond).
C CONSTS(159) is v0 in component ICaL (millivolt).
C ALGBRC(74) is A_1 in component ICaL (dimensionless).
C CONSTS(175) is B_1 in component ICaL (per_millivolt).
C ALGBRC(75) is U_1 in component ICaL (dimensionless).
C ALGBRC(77) is A_2 in component ICaL (dimensionless).
C CONSTS(176) is B_2 in component ICaL (per_millivolt).
C ALGBRC(78) is U_2 in component ICaL (dimensionless).
C ALGBRC(80) is A_3 in component ICaL (dimensionless).
C CONSTS(177) is B_3 in component ICaL (per_millivolt).
C ALGBRC(81) is U_3 in component ICaL (dimensionless).
C CONSTS(54) is GKr_b in component IKr (milliS_per_microF).
C STATES(35) is IC1 in component IKr (dimensionless).
C STATES(36) is IC2 in component IKr (dimensionless).
C STATES(37) is C1 in component IKr (dimensionless).
C STATES(38) is C2 in component IKr (dimensionless).
C STATES(39) is O in component IKr (dimensionless).
C STATES(40) is IO in component IKr (dimensionless).
C STATES(41) is IObound in component IKr (dimensionless).
C STATES(42) is Obound in component IKr (dimensionless).
C STATES(43) is Cbound in component IKr (dimensionless).
C STATES(44) is D in component IKr (dimensionless).
C CONSTS(160) is GKr in component IKr (milliS_per_microF).
C CONSTS(55) is A1 in component IKr (per_millisecond).
C CONSTS(56) is B1 in component IKr (per_millivolt).
C CONSTS(57) is q1 in component IKr (dimensionless).
C CONSTS(58) is A2 in component IKr (per_millisecond).
C CONSTS(59) is B2 in component IKr (per_millivolt).
C CONSTS(60) is q2 in component IKr (dimensionless).
C CONSTS(61) is A3 in component IKr (per_millisecond).
C CONSTS(62) is B3 in component IKr (per_millivolt).
C CONSTS(63) is q3 in component IKr (dimensionless).
C CONSTS(64) is A4 in component IKr (per_millisecond).
C CONSTS(65) is B4 in component IKr (per_millivolt).
C CONSTS(66) is q4 in component IKr (dimensionless).
C CONSTS(67) is A11 in component IKr (per_millisecond).
C CONSTS(68) is B11 in component IKr (per_millivolt).
C CONSTS(69) is q11 in component IKr (dimensionless).
C CONSTS(70) is A21 in component IKr (per_millisecond).
C CONSTS(71) is B21 in component IKr (per_millivolt).
C CONSTS(72) is q21 in component IKr (dimensionless).
C CONSTS(73) is A31 in component IKr (per_millisecond).
C CONSTS(74) is B31 in component IKr (per_millivolt).
C CONSTS(75) is q31 in component IKr (dimensionless).
C CONSTS(76) is A41 in component IKr (per_millisecond).
C CONSTS(77) is B41 in component IKr (per_millivolt).
C CONSTS(78) is q41 in component IKr (dimensionless).
C CONSTS(79) is A51 in component IKr (per_millisecond).
C CONSTS(80) is B51 in component IKr (per_millivolt).
C CONSTS(81) is q51 in component IKr (dimensionless).
C CONSTS(82) is A52 in component IKr (per_millisecond).
C CONSTS(83) is B52 in component IKr (per_millivolt).
C CONSTS(84) is q52 in component IKr (dimensionless).
C CONSTS(85) is A53 in component IKr (per_millisecond).
C CONSTS(86) is B53 in component IKr (per_millivolt).
C CONSTS(87) is q53 in component IKr (dimensionless).
C CONSTS(88) is A61 in component IKr (per_millisecond).
C CONSTS(89) is B61 in component IKr (per_millivolt).
C CONSTS(90) is q61 in component IKr (dimensionless).
C CONSTS(91) is A62 in component IKr (per_millisecond).
C CONSTS(92) is B62 in component IKr (per_millivolt).
C CONSTS(93) is q62 in component IKr (dimensionless).
C CONSTS(94) is A63 in component IKr (per_millisecond).
C CONSTS(95) is B63 in component IKr (per_millivolt).
C CONSTS(96) is q63 in component IKr (dimensionless).
C CONSTS(97) is Kmax in component IKr (dimensionless).
C CONSTS(98) is Ku in component IKr (per_millisecond).
C CONSTS(99) is n in component IKr (dimensionless).
C CONSTS(100) is halfmax in component IKr (dimensionless).
C CONSTS(101) is Kt in component IKr (per_millisecond).
C CONSTS(102) is Vhalf in component IKr (millivolt).
C CONSTS(103) is Temp in component IKr (dimensionless).
C CONSTS(104) is GKs_b in component IKs (milliS_per_microF).
C CONSTS(161) is GKs in component IKs (milliS_per_microF).
C ALGBRC(11) is xs1ss in component IKs (dimensionless).
C ALGBRC(25) is xs2ss in component IKs (dimensionless).
C ALGBRC(26) is txs1 in component IKs (millisecond).
C CONSTS(105) is txs1_max in component IKs (millisecond).
C STATES(45) is xs1 in component IKs (dimensionless).
C STATES(46) is xs2 in component IKs (dimensionless).
C ALGBRC(94) is KsCa in component IKs (dimensionless).
C ALGBRC(36) is txs2 in component IKs (millisecond).
C CONSTS(162) is GK1 in component IK1 (milliS_per_microF).
C CONSTS(106) is GK1_b in component IK1 (milliS_per_microF).
C ALGBRC(12) is xk1ss in component IK1 (dimensionless).
C ALGBRC(27) is txk1 in component IK1 (millisecond).
C STATES(47) is xk1 in component IK1 (dimensionless).
C ALGBRC(98) is rk1 in component IK1 (millisecond).
C CONSTS(107) is kna1 in component INaCa_i (per_millisecond).
C CONSTS(108) is kna2 in component INaCa_i (per_millisecond).
C CONSTS(109) is kna3 in component INaCa_i (per_millisecond).
C CONSTS(110) is kasymm in component INaCa_i (dimensionless).
C CONSTS(111) is wna in component INaCa_i (dimensionless).
C CONSTS(112) is wca in component INaCa_i (dimensionless).
C CONSTS(113) is wnaca in component INaCa_i (dimensionless).
C CONSTS(114) is kcaon in component INaCa_i (per_millisecond).
C CONSTS(115) is kcaoff in component INaCa_i (per_millisecond).
C CONSTS(116) is qna in component INaCa_i (dimensionless).
C CONSTS(117) is qca in component INaCa_i (dimensionless).
C ALGBRC(101) is hna in component INaCa_i (dimensionless).
C ALGBRC(100) is hca in component INaCa_i (dimensionless).
C CONSTS(118) is KmCaAct in component INaCa_i (millimolar).
C CONSTS(119) is Gncx_b in component INaCa_i (milliS_per_microF).
C CONSTS(195) is Gncx in component INaCa_i (milliS_per_microF).
C ALGBRC(102) is h1_i in component INaCa_i (dimensionless).
C ALGBRC(103) is h2_i in component INaCa_i (dimensionless).
C ALGBRC(104) is h3_i in component INaCa_i (dimensionless).
C ALGBRC(105) is h4_i in component INaCa_i (dimensionless).
C ALGBRC(106) is h5_i in component INaCa_i (dimensionless).
C ALGBRC(107) is h6_i in component INaCa_i (dimensionless).
C ALGBRC(108) is h7_i in component INaCa_i (dimensionless).
C ALGBRC(109) is h8_i in component INaCa_i (dimensionless).
C ALGBRC(110) is h9_i in component INaCa_i (dimensionless).
C CONSTS(189) is h10_i in component INaCa_i (dimensionless).
C CONSTS(190) is h11_i in component INaCa_i (dimensionless).
C CONSTS(191) is h12_i in component INaCa_i (dimensionless).
C CONSTS(192) is k1_i in component INaCa_i (dimensionless).
C CONSTS(193) is k2_i in component INaCa_i (dimensionless).
C ALGBRC(111) is k3p_i in component INaCa_i (dimensionless).
C ALGBRC(112) is k3pp_i in component INaCa_i (dimensionless).
C ALGBRC(113) is k3_i in component INaCa_i (dimensionless).
C ALGBRC(116) is k4_i in component INaCa_i (dimensionless).
C ALGBRC(114) is k4p_i in component INaCa_i (dimensionless).
C ALGBRC(115) is k4pp_i in component INaCa_i (dimensionless).
C CONSTS(194) is k5_i in component INaCa_i (dimensionless).
C ALGBRC(117) is k6_i in component INaCa_i (dimensionless).
C ALGBRC(118) is k7_i in component INaCa_i (dimensionless).
C ALGBRC(119) is k8_i in component INaCa_i (dimensionless).
C ALGBRC(120) is x1_i in component INaCa_i (dimensionless).
C ALGBRC(121) is x2_i in component INaCa_i (dimensionless).
C ALGBRC(122) is x3_i in component INaCa_i (dimensionless).
C ALGBRC(123) is x4_i in component INaCa_i (dimensionless).
C ALGBRC(124) is E1_i in component INaCa_i (dimensionless).
C ALGBRC(125) is E2_i in component INaCa_i (dimensionless).
C ALGBRC(126) is E3_i in component INaCa_i (dimensionless).
C ALGBRC(127) is E4_i in component INaCa_i (dimensionless).
C ALGBRC(128) is allo_i in component INaCa_i (dimensionless).
C ALGBRC(129) is JncxNa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(130) is JncxCa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(132) is h1_ss in component INaCa_i (dimensionless).
C ALGBRC(133) is h2_ss in component INaCa_i (dimensionless).
C ALGBRC(134) is h3_ss in component INaCa_i (dimensionless).
C ALGBRC(135) is h4_ss in component INaCa_i (dimensionless).
C ALGBRC(136) is h5_ss in component INaCa_i (dimensionless).
C ALGBRC(137) is h6_ss in component INaCa_i (dimensionless).
C ALGBRC(138) is h7_ss in component INaCa_i (dimensionless).
C ALGBRC(139) is h8_ss in component INaCa_i (dimensionless).
C ALGBRC(140) is h9_ss in component INaCa_i (dimensionless).
C CONSTS(196) is h10_ss in component INaCa_i (dimensionless).
C CONSTS(197) is h11_ss in component INaCa_i (dimensionless).
C CONSTS(198) is h12_ss in component INaCa_i (dimensionless).
C CONSTS(199) is k1_ss in component INaCa_i (dimensionless).
C CONSTS(200) is k2_ss in component INaCa_i (dimensionless).
C ALGBRC(141) is k3p_ss in component INaCa_i (dimensionless).
C ALGBRC(142) is k3pp_ss in component INaCa_i (dimensionless).
C ALGBRC(143) is k3_ss in component INaCa_i (dimensionless).
C ALGBRC(146) is k4_ss in component INaCa_i (dimensionless).
C ALGBRC(144) is k4p_ss in component INaCa_i (dimensionless).
C ALGBRC(145) is k4pp_ss in component INaCa_i (dimensionless).
C CONSTS(201) is k5_ss in component INaCa_i (dimensionless).
C ALGBRC(147) is k6_ss in component INaCa_i (dimensionless).
C ALGBRC(148) is k7_ss in component INaCa_i (dimensionless).
C ALGBRC(149) is k8_ss in component INaCa_i (dimensionless).
C ALGBRC(150) is x1_ss in component INaCa_i (dimensionless).
C ALGBRC(151) is x2_ss in component INaCa_i (dimensionless).
C ALGBRC(152) is x3_ss in component INaCa_i (dimensionless).
C ALGBRC(153) is x4_ss in component INaCa_i (dimensionless).
C ALGBRC(154) is E1_ss in component INaCa_i (dimensionless).
C ALGBRC(155) is E2_ss in component INaCa_i (dimensionless).
C ALGBRC(156) is E3_ss in component INaCa_i (dimensionless).
C ALGBRC(157) is E4_ss in component INaCa_i (dimensionless).
C ALGBRC(158) is allo_ss in component INaCa_i (dimensionless).
C ALGBRC(159) is JncxNa_ss in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(160) is JncxCa_ss in component INaCa_i (millimolar_per_millisecond).
C CONSTS(120) is k1p in component INaK (per_millisecond).
C CONSTS(121) is k1m in component INaK (per_millisecond).
C CONSTS(122) is k2p in component INaK (per_millisecond).
C CONSTS(123) is k2m in component INaK (per_millisecond).
C CONSTS(124) is k3p in component INaK (per_millisecond).
C CONSTS(125) is k3m in component INaK (per_millisecond).
C CONSTS(126) is k4p in component INaK (per_millisecond).
C CONSTS(127) is k4m in component INaK (per_millisecond).
C CONSTS(128) is Knai0 in component INaK (millimolar).
C CONSTS(129) is Knao0 in component INaK (millimolar).
C CONSTS(130) is delta in component INaK (millivolt).
C CONSTS(131) is Kki in component INaK (per_millisecond).
C CONSTS(132) is Kko in component INaK (per_millisecond).
C CONSTS(133) is MgADP in component INaK (millimolar).
C CONSTS(134) is MgATP in component INaK (millimolar).
C CONSTS(135) is Kmgatp in component INaK (millimolar).
C CONSTS(136) is H in component INaK (millimolar).
C CONSTS(137) is eP in component INaK (dimensionless).
C CONSTS(138) is Khp in component INaK (millimolar).
C CONSTS(139) is Knap in component INaK (millimolar).
C CONSTS(140) is Kxkur in component INaK (millimolar).
C CONSTS(141) is Pnak_b in component INaK (milliS_per_microF).
C CONSTS(205) is Pnak in component INaK (milliS_per_microF).
C ALGBRC(162) is Knai in component INaK (millimolar).
C ALGBRC(163) is Knao in component INaK (millimolar).
C ALGBRC(164) is P in component INaK (dimensionless).
C ALGBRC(165) is a1 in component INaK (dimensionless).
C CONSTS(202) is b1 in component INaK (dimensionless).
C CONSTS(203) is a2 in component INaK (dimensionless).
C ALGBRC(166) is b2 in component INaK (dimensionless).
C ALGBRC(167) is a3 in component INaK (dimensionless).
C ALGBRC(168) is b3 in component INaK (dimensionless).
C CONSTS(204) is a4 in component INaK (dimensionless).
C ALGBRC(169) is b4 in component INaK (dimensionless).
C ALGBRC(170) is x1 in component INaK (dimensionless).
C ALGBRC(171) is x2 in component INaK (dimensionless).
C ALGBRC(172) is x3 in component INaK (dimensionless).
C ALGBRC(173) is x4 in component INaK (dimensionless).
C ALGBRC(174) is E1 in component INaK (dimensionless).
C ALGBRC(175) is E2 in component INaK (dimensionless).
C ALGBRC(176) is E3 in component INaK (dimensionless).
C ALGBRC(177) is E4 in component INaK (dimensionless).
C ALGBRC(178) is JnakNa in component INaK (millimolar_per_millisecond).
C ALGBRC(179) is JnakK in component INaK (millimolar_per_millisecond).
C ALGBRC(181) is xkb in component IKb (dimensionless).
C CONSTS(142) is GKb_b in component IKb (milliS_per_microF).
C CONSTS(164) is GKb in component IKb (milliS_per_microF).
C CONSTS(143) is PNab in component INab (milliS_per_microF).
C ALGBRC(183) is A in component INab (microA_per_microF).
C CONSTS(179) is B in component INab (per_millivolt).
C CONSTS(165) is v0 in component INab (millivolt).
C ALGBRC(185) is U in component INab (dimensionless).
C CONSTS(144) is PCab in component ICab (milliS_per_microF).
C ALGBRC(187) is A in component ICab (microA_per_microF).
C CONSTS(180) is B in component ICab (per_millivolt).
C CONSTS(166) is v0 in component ICab (millivolt).
C ALGBRC(189) is U in component ICab (dimensionless).
C CONSTS(145) is GpCa in component IpCa (milliS_per_microF).
C CONSTS(146) is KmCap in component IpCa (millimolar).
C CONSTS(147) is bt in component ryr (millisecond).
C CONSTS(167) is a_rel in component ryr (millisecond).
C ALGBRC(89) is Jrel_inf in component ryr (dimensionless).
C ALGBRC(95) is tau_rel in component ryr (millisecond).
C ALGBRC(90) is Jrel_infp in component ryr (dimensionless).
C ALGBRC(87) is Jrel_temp in component ryr (dimensionless).
C ALGBRC(96) is tau_relp in component ryr (millisecond).
C STATES(48) is Jrelnp in component ryr (dimensionless).
C STATES(49) is Jrelp in component ryr (dimensionless).
C CONSTS(168) is btp in component ryr (millisecond).
C CONSTS(181) is a_relp in component ryr (millisecond).
C ALGBRC(86) is Jrel_inf_temp in component ryr (dimensionless).
C ALGBRC(193) is fJrelp in component ryr (dimensionless).
C CONSTS(148) is Jrel_scaling_factor in component ryr (dimensionless).
C ALGBRC(92) is tau_rel_temp in component ryr (millisecond).
C ALGBRC(93) is tau_relp_temp in component ryr (millisecond).
C CONSTS(169) is upScale in component SERCA (dimensionless).
C ALGBRC(195) is Jupnp in component SERCA (millimolar_per_millisecond).
C ALGBRC(196) is Jupp in component SERCA (millimolar_per_millisecond).
C ALGBRC(197) is fJupp in component SERCA (dimensionless).
C ALGBRC(198) is Jleak in component SERCA (millimolar_per_millisecond).
C CONSTS(149) 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 hf in component INa (dimensionless).
C RATES(13) is d/dt hs in component INa (dimensionless).
C RATES(14) is d/dt j in component INa (dimensionless).
C RATES(15) is d/dt hsp in component INa (dimensionless).
C RATES(16) is d/dt jp in component INa (dimensionless).
C RATES(17) is d/dt mL in component INaL (dimensionless).
C RATES(18) is d/dt hL in component INaL (dimensionless).
C RATES(19) is d/dt hLp in component INaL (dimensionless).
C RATES(20) is d/dt a in component Ito (dimensionless).
C RATES(21) is d/dt iF in component Ito (dimensionless).
C RATES(22) is d/dt iS in component Ito (dimensionless).
C RATES(23) is d/dt ap in component Ito (dimensionless).
C RATES(24) is d/dt iFp in component Ito (dimensionless).
C RATES(25) is d/dt iSp in component Ito (dimensionless).
C RATES(26) is d/dt d in component ICaL (dimensionless).
C RATES(27) is d/dt ff in component ICaL (dimensionless).
C RATES(28) is d/dt fs in component ICaL (dimensionless).
C RATES(29) is d/dt fcaf in component ICaL (dimensionless).
C RATES(30) is d/dt fcas in component ICaL (dimensionless).
C RATES(31) is d/dt jca in component ICaL (dimensionless).
C RATES(32) is d/dt ffp in component ICaL (dimensionless).
C RATES(33) is d/dt fcafp in component ICaL (dimensionless).
C RATES(34) is d/dt nca in component ICaL (dimensionless).
C RATES(35) is d/dt IC1 in component IKr (dimensionless).
C RATES(36) is d/dt IC2 in component IKr (dimensionless).
C RATES(37) is d/dt C1 in component IKr (dimensionless).
C RATES(38) is d/dt C2 in component IKr (dimensionless).
C RATES(39) is d/dt O in component IKr (dimensionless).
C RATES(40) is d/dt IO in component IKr (dimensionless).
C RATES(41) is d/dt IObound in component IKr (dimensionless).
C RATES(42) is d/dt Obound in component IKr (dimensionless).
C RATES(43) is d/dt Cbound in component IKr (dimensionless).
C RATES(44) is d/dt D in component IKr (dimensionless).
C RATES(45) is d/dt xs1 in component IKs (dimensionless).
C RATES(46) is d/dt xs2 in component IKs (dimensionless).
C RATES(47) is d/dt xk1 in component IK1 (dimensionless).
C RATES(48) is d/dt Jrelnp in component ryr (dimensionless).
C RATES(49) is d/dt Jrelp in component ryr (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0
      CONSTS(2) = 140
      CONSTS(3) = 1.8
      CONSTS(4) = 5.4
      CONSTS(5) = 8314
      CONSTS(6) = 310
      CONSTS(7) = 96485
      CONSTS(8) = 1
      CONSTS(9) = 2
      CONSTS(10) = 1
      CONSTS(11) = 0.01
      CONSTS(12) = 0.0011
      STATES(1) = -88.00190465
      CONSTS(13) = 10
      CONSTS(14) = 100000000000000000
      CONSTS(15) = -80
      CONSTS(16) = 1000
      CONSTS(17) = 0.5
      CONSTS(18) = 0.15
      CONSTS(19) = 0.05
      CONSTS(20) = 0.00068
      CONSTS(21) = 0.05
      CONSTS(22) = 0.0015
      STATES(2) = 0.0125840447
      STATES(3) = 8.49e-05
      CONSTS(23) = 0.05
      CONSTS(24) = 0.00238
      CONSTS(25) = 0.07
      CONSTS(26) = 0.0005
      CONSTS(27) = 0.047
      CONSTS(28) = 0.00087
      CONSTS(29) = 1.124
      CONSTS(30) = 0.0087
      CONSTS(31) = 10
      CONSTS(32) = 0.8
      STATES(4) = 7.268004498
      STATES(5) = 7.268089977
      STATES(6) = 144.6555918
      STATES(7) = 144.6555651
      STATES(8) = 1.619574538
      STATES(9) = 1.571234014
      STATES(10) = 8.6e-05
      CONSTS(33) = 1
      CONSTS(34) = 0.01833
      CONSTS(35) = 39.57
      CONSTS(36) = 9.871
      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(11) = 0.007344121102
      CONSTS(43) = 82.9
      CONSTS(44) = 6.086
      CONSTS(45) = 0.99
      STATES(12) = 0.6981071913
      STATES(13) = 0.6980895801
      CONSTS(46) = 75
      CONSTS(47) = 0
      STATES(14) = 0.6979908432
      STATES(15) = 0.4549485525
      STATES(16) = 0.6979245865
      STATES(17) = 0.0001882617273
      CONSTS(48) = 200
      STATES(18) = 0.5008548855
      STATES(19) = 0.2693065357
      CONSTS(49) = 0.019957499999999975
      CONSTS(50) = 0.02
      STATES(20) = 0.001001097687
      STATES(21) = 0.9995541745
      STATES(22) = 0.5865061736
      STATES(23) = 0.0005100862934
      STATES(24) = 0.9995541823
      STATES(25) = 0.6393399482
      CONSTS(51) = 0.002
      CONSTS(52) = 1000
      CONSTS(53) = 0.0001007
      STATES(26) = 2.34e-9
      STATES(27) = 0.9999999909
      STATES(28) = 0.9102412777
      STATES(29) = 0.9999999909
      STATES(30) = 0.9998046777
      STATES(31) = 0.9999738312
      STATES(32) = 0.9999999909
      STATES(33) = 0.9999999909
      STATES(34) = 0.002749414044
      CONSTS(54) = 0.04658545454545456
      STATES(35) = 0.999637
      STATES(36) = 6.83208e-05
      STATES(37) = 1.80145e-08
      STATES(38) = 8.26619e-05
      STATES(39) = 0.00015551
      STATES(40) = 5.67623e-05
      STATES(41) = 0
      STATES(42) = 0
      STATES(43) = 0
      STATES(44) = 0
      CONSTS(55) = 0.0264
      CONSTS(56) = 4.631E-05
      CONSTS(57) = 4.843
      CONSTS(58) = 4.986E-06
      CONSTS(59) = -0.004226
      CONSTS(60) = 4.23
      CONSTS(61) = 0.001214
      CONSTS(62) = 0.008516
      CONSTS(63) = 4.962
      CONSTS(64) = 1.854E-05
      CONSTS(65) = -0.04641
      CONSTS(66) = 3.769
      CONSTS(67) = 0.0007868
      CONSTS(68) = 1.535E-08
      CONSTS(69) = 4.942
      CONSTS(70) = 5.455E-06
      CONSTS(71) = -0.1688
      CONSTS(72) = 4.156
      CONSTS(73) = 0.005509
      CONSTS(74) = 7.771E-09
      CONSTS(75) = 4.22
      CONSTS(76) = 0.001416
      CONSTS(77) = -0.02877
      CONSTS(78) = 1.459
      CONSTS(79) = 0.4492
      CONSTS(80) = 0.008595
      CONSTS(81) = 5
      CONSTS(82) = 0.3181
      CONSTS(83) = 3.613E-08
      CONSTS(84) = 4.663
      CONSTS(85) = 0.149
      CONSTS(86) = 0.004668
      CONSTS(87) = 2.412
      CONSTS(88) = 0.01241
      CONSTS(89) = 0.1725
      CONSTS(90) = 5.568
      CONSTS(91) = 0.3226
      CONSTS(92) = -0.0006575
      CONSTS(93) = 5
      CONSTS(94) = 0.008978
      CONSTS(95) = -0.02215
      CONSTS(96) = 5.682
      CONSTS(97) = 0
      CONSTS(98) = 0
      CONSTS(99) = 1
      CONSTS(100) = 1
      CONSTS(101) = 0
      CONSTS(102) = 1
      CONSTS(103) = 37
      CONSTS(104) = 0.006358000000000001
      CONSTS(105) = 817.3
      STATES(45) = 0.2707758025
      STATES(46) = 0.0001928503426
      CONSTS(106) = 0.3239783999999998
      STATES(47) = 0.9967597594
      CONSTS(107) = 15
      CONSTS(108) = 5
      CONSTS(109) = 88.12
      CONSTS(110) = 12.5
      CONSTS(111) = 6e4
      CONSTS(112) = 6e4
      CONSTS(113) = 5e3
      CONSTS(114) = 1.5e6
      CONSTS(115) = 5e3
      CONSTS(116) = 0.5224
      CONSTS(117) = 0.167
      CONSTS(118) = 150e-6
      CONSTS(119) = 0.0008
      CONSTS(120) = 949.5
      CONSTS(121) = 182.4
      CONSTS(122) = 687.2
      CONSTS(123) = 39.4
      CONSTS(124) = 1899
      CONSTS(125) = 79300
      CONSTS(126) = 639
      CONSTS(127) = 40
      CONSTS(128) = 9.073
      CONSTS(129) = 27.78
      CONSTS(130) = -0.155
      CONSTS(131) = 0.5
      CONSTS(132) = 0.3582
      CONSTS(133) = 0.05
      CONSTS(134) = 9.8
      CONSTS(135) = 1.698e-7
      CONSTS(136) = 1e-7
      CONSTS(137) = 4.2
      CONSTS(138) = 1.698e-7
      CONSTS(139) = 224
      CONSTS(140) = 292
      CONSTS(141) = 30
      CONSTS(142) = 0.003
      CONSTS(143) = 3.75e-10
      CONSTS(144) = 2.5e-8
      CONSTS(145) = 0.0005
      CONSTS(146) = 0.0005
      CONSTS(147) = 4.75
      STATES(48) = 2.5e-7
      STATES(49) = 3.12e-7
      CONSTS(148) = 1.0
      CONSTS(149) = 1.0
      CONSTS(150) = CONSTS(7)/( CONSTS(5)*CONSTS(6))
      CONSTS(151) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(23)*1.30000, CONSTS(23))
      CONSTS(152) = 1.00000 - CONSTS(45)
      CONSTS(153) =  3.00000*CONSTS(48)
      CONSTS(154) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(49)*0.600000, CONSTS(49))
      CONSTS(155) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(50)*4.00000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(50)*4.00000, CONSTS(50))
      CONSTS(156) = 0.600000
      CONSTS(157) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(53)*1.20000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(53)*2.50000, CONSTS(53))
      CONSTS(158) = 75.0000
      CONSTS(159) = 0.00000
      CONSTS(160) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(54)*1.30000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(54)*0.800000, CONSTS(54))
      CONSTS(161) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(104)*1.40000, CONSTS(104))
      CONSTS(162) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(106)*1.20000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(106)*1.30000, CONSTS(106))
      CONSTS(163) =  1000.00*3.14000*CONSTS(12)*CONSTS(12)*CONSTS(11)
      CONSTS(164) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(142)*0.600000, CONSTS(142))
      CONSTS(165) = 0.00000
      CONSTS(166) = 0.00000
      CONSTS(167) =  0.500000*CONSTS(147)
      CONSTS(168) =  1.25000*CONSTS(147)
      CONSTS(169) = TERNRY(CONSTS(1).EQ.1.00000, 1.30000, 1.00000)
      CONSTS(205) = 0.00000
      CONSTS(170) =  CONSTS(7)*CONSTS(150)
      CONSTS(171) = 1.00000 - CONSTS(156)
      CONSTS(172) =  1.10000*CONSTS(157)
      CONSTS(173) =  0.00125000*CONSTS(157)
      CONSTS(174) =  0.000357400*CONSTS(157)
      CONSTS(175) =  2.00000*CONSTS(150)
      CONSTS(176) = CONSTS(150)
      CONSTS(177) = CONSTS(150)
      CONSTS(178) =  2.00000*3.14000*CONSTS(12)*CONSTS(12)+ 2.00000*3.14000*CONSTS(12)*CONSTS(11)
      CONSTS(179) = CONSTS(150)
      CONSTS(180) =  2.00000*CONSTS(150)
      CONSTS(181) =  0.500000*CONSTS(168)
      CONSTS(182) =  0.00125000*CONSTS(172)
      CONSTS(183) =  0.000357400*CONSTS(172)
      CONSTS(184) =  2.00000*CONSTS(178)
      CONSTS(185) =  0.680000*CONSTS(163)
      CONSTS(186) =  0.0552000*CONSTS(163)
      CONSTS(187) =  0.00480000*CONSTS(163)
      CONSTS(188) =  0.0200000*CONSTS(163)
      CONSTS(189) = CONSTS(110)+1.00000+ (CONSTS(2)/CONSTS(107))*(1.00000+CONSTS(2)/CONSTS(108))
      CONSTS(190) = ( CONSTS(2)*CONSTS(2))/( CONSTS(189)*CONSTS(107)*CONSTS(108))
      CONSTS(191) = 1.00000/CONSTS(189)
      CONSTS(192) =  CONSTS(191)*CONSTS(3)*CONSTS(114)
      CONSTS(193) = CONSTS(115)
      CONSTS(194) = CONSTS(115)
      CONSTS(195) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(119)*1.10000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(119)*1.40000, CONSTS(119))
      CONSTS(196) = CONSTS(110)+1.00000+ (CONSTS(2)/CONSTS(107))*(1.00000+CONSTS(2)/CONSTS(108))
      CONSTS(197) = ( CONSTS(2)*CONSTS(2))/( CONSTS(196)*CONSTS(107)*CONSTS(108))
      CONSTS(198) = 1.00000/CONSTS(196)
      CONSTS(199) =  CONSTS(198)*CONSTS(3)*CONSTS(114)
      CONSTS(200) = CONSTS(115)
      CONSTS(201) = CONSTS(115)
      CONSTS(202) =  CONSTS(121)*CONSTS(133)
      CONSTS(203) = CONSTS(122)
      CONSTS(204) = (( CONSTS(126)*CONSTS(134))/CONSTS(135))/(1.00000+CONSTS(134)/CONSTS(135))
      CONSTS(205) = TERNRY(CONSTS(1).EQ.1.00000,  CONSTS(141)*0.900000, TERNRY(CONSTS(1).EQ.2.00000,  CONSTS(141)*0.700000, CONSTS(141))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(44) = CONSTS(205)
      RATES(35) = (- ( CONSTS(67)*EXP( CONSTS(68)*STATES(1))*STATES(35)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(69)))/10.0000) -  CONSTS(70)*EXP( CONSTS(71)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(72)))/10.0000))+ CONSTS(79)*EXP( CONSTS(80)*STATES(1))*STATES(37)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(81)))/10.0000)) -  CONSTS(88)*EXP( CONSTS(89)*STATES(1))*STATES(35)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(90)))/10.0000)
      RATES(36) = ((( CONSTS(67)*EXP( CONSTS(68)*STATES(1))*STATES(35)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(69)))/10.0000) -  CONSTS(70)*EXP( CONSTS(71)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(72)))/10.0000)) - ( CONSTS(61)*EXP( CONSTS(62)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(63)))/10.0000) -  CONSTS(64)*EXP( CONSTS(65)*STATES(1))*STATES(40)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(66)))/10.0000)))+ CONSTS(82)*EXP( CONSTS(83)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(84)))/10.0000)) -  CONSTS(91)*EXP( CONSTS(92)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(93)))/10.0000)
      RATES(37) = - ( CONSTS(55)*EXP( CONSTS(56)*STATES(1))*STATES(37)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(57)))/10.0000) -  CONSTS(58)*EXP( CONSTS(59)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(60)))/10.0000)) - ( CONSTS(79)*EXP( CONSTS(80)*STATES(1))*STATES(37)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(81)))/10.0000) -  CONSTS(88)*EXP( CONSTS(89)*STATES(1))*STATES(35)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(90)))/10.0000))
      RATES(38) = (( CONSTS(55)*EXP( CONSTS(56)*STATES(1))*STATES(37)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(57)))/10.0000) -  CONSTS(58)*EXP( CONSTS(59)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(60)))/10.0000)) - ( CONSTS(73)*EXP( CONSTS(74)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(75)))/10.0000) -  CONSTS(76)*EXP( CONSTS(77)*STATES(1))*STATES(39)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(78)))/10.0000))) - ( CONSTS(82)*EXP( CONSTS(83)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(84)))/10.0000) -  CONSTS(91)*EXP( CONSTS(92)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(93)))/10.0000))
      RATES(39) = (( CONSTS(73)*EXP( CONSTS(74)*STATES(1))*STATES(38)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(75)))/10.0000) -  CONSTS(76)*EXP( CONSTS(77)*STATES(1))*STATES(39)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(78)))/10.0000)) - ( CONSTS(85)*EXP( CONSTS(86)*STATES(1))*STATES(39)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(87)))/10.0000) -  CONSTS(94)*EXP( CONSTS(95)*STATES(1))*STATES(40)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(96)))/10.0000))) - ( (( CONSTS(97)*CONSTS(98)*EXP( CONSTS(99)*log(STATES(44))))/(EXP( CONSTS(99)*log(STATES(44)))+CONSTS(100)))*STATES(39) -  CONSTS(98)*STATES(42))
      RATES(40) = ((( CONSTS(61)*EXP( CONSTS(62)*STATES(1))*STATES(36)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(63)))/10.0000) -  CONSTS(64)*EXP( CONSTS(65)*STATES(1))*STATES(40)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(66)))/10.0000))+ CONSTS(85)*EXP( CONSTS(86)*STATES(1))*STATES(39)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(87)))/10.0000)) -  CONSTS(94)*EXP( CONSTS(95)*STATES(1))*STATES(40)*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(96)))/10.0000)) - ( (( CONSTS(97)*CONSTS(98)*EXP( CONSTS(99)*log(STATES(44))))/(EXP( CONSTS(99)*log(STATES(44)))+CONSTS(100)))*STATES(40) -  (( CONSTS(98)*CONSTS(85)*EXP( CONSTS(86)*STATES(1))*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(87)))/10.0000))/( CONSTS(94)*EXP( CONSTS(95)*STATES(1))*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(96)))/10.0000)))*STATES(41))
      RATES(41) = (( (( CONSTS(97)*CONSTS(98)*EXP( CONSTS(99)*log(STATES(44))))/(EXP( CONSTS(99)*log(STATES(44)))+CONSTS(100)))*STATES(40) -  (( CONSTS(98)*CONSTS(85)*EXP( CONSTS(86)*STATES(1))*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(87)))/10.0000))/( CONSTS(94)*EXP( CONSTS(95)*STATES(1))*EXP(( (CONSTS(103) - 20.0000)*log(CONSTS(96)))/10.0000)))*STATES(41))+ (CONSTS(101)/(1.00000+EXP(- (STATES(1) - CONSTS(102))/6.78900)))*STATES(43)) -  CONSTS(101)*STATES(41)
      RATES(42) = (( (( CONSTS(97)*CONSTS(98)*EXP( CONSTS(99)*log(STATES(44))))/(EXP( CONSTS(99)*log(STATES(44)))+CONSTS(100)))*STATES(39) -  CONSTS(98)*STATES(42))+ (CONSTS(101)/(1.00000+EXP(- (STATES(1) - CONSTS(102))/6.78900)))*STATES(43)) -  CONSTS(101)*STATES(42)
      RATES(43) = - ( (CONSTS(101)/(1.00000+EXP(- (STATES(1) - CONSTS(102))/6.78900)))*STATES(43) -  CONSTS(101)*STATES(42)) - ( (CONSTS(101)/(1.00000+EXP(- (STATES(1) - CONSTS(102))/6.78900)))*STATES(43) -  CONSTS(101)*STATES(41))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
      RATES(18) = (ALGBRC(4) - STATES(18))/CONSTS(48)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
      RATES(19) = (ALGBRC(5) - STATES(19))/CONSTS(153)
      ALGBRC(2) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(35))/CONSTS(36)))
      ALGBRC(14) = 1.00000/( CONSTS(39)*EXP((STATES(1)+CONSTS(37))/CONSTS(38))+ CONSTS(40)*EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
      RATES(11) = (ALGBRC(2) - STATES(11))/ALGBRC(14)
      ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(43)) - CONSTS(47))/CONSTS(44)))
      ALGBRC(15) = 1.00000/( 1.43200e-05*EXP(- ((STATES(1)+1.19600) - CONSTS(47))/6.28500)+ 6.14900*EXP(((STATES(1)+0.509600) - CONSTS(47))/20.2700))
      RATES(12) = (ALGBRC(3) - STATES(12))/ALGBRC(15)
      ALGBRC(16) = 1.00000/( 0.00979400*EXP(- ((STATES(1)+17.9500) - CONSTS(47))/28.0500)+ 0.334300*EXP(((STATES(1)+5.73000) - CONSTS(47))/56.6600))
      RATES(13) = (ALGBRC(3) - STATES(13))/ALGBRC(16)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.3400)/14.8200))
      ALGBRC(18) = 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(20) = (ALGBRC(6) - STATES(20))/ALGBRC(18)
      ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000)/4.23000))
      ALGBRC(22) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
      RATES(26) = (ALGBRC(8) - STATES(26))/ALGBRC(22)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
      ALGBRC(23) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000)/10.0000))
      RATES(27) = (ALGBRC(9) - STATES(27))/ALGBRC(23)
      ALGBRC(24) = 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(28) = (ALGBRC(9) - STATES(28))/ALGBRC(24)
      ALGBRC(20) = ALGBRC(9)
      RATES(31) = (ALGBRC(20) - STATES(31))/CONSTS(158)
      ALGBRC(10) =  STATES(31)*1.00000
      ALGBRC(21) = 1.00000/(CONSTS(52)/ALGBRC(10)+1.00000+CONSTS(51)/STATES(3) ** 4.00000)
      RATES(34) =  ALGBRC(21)*CONSTS(52) -  STATES(34)*ALGBRC(10)
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(26) = CONSTS(105)+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      RATES(45) = (ALGBRC(11) - STATES(45))/ALGBRC(26)
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(4)+144.590)/( 1.56920*CONSTS(4)+3.81150)))
      ALGBRC(27) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
      RATES(47) = (ALGBRC(12) - STATES(47))/ALGBRC(27)
      ALGBRC(37) = ( CONSTS(21)*(1.00000 - STATES(2)))/(1.00000+CONSTS(22)/STATES(3))
      RATES(2) =  CONSTS(19)*ALGBRC(37)*(ALGBRC(37)+STATES(2)) -  CONSTS(20)*STATES(2)
      ALGBRC(17) = ALGBRC(3)
      ALGBRC(28) = 2.03800+1.00000/( 0.0213600*EXP(- ((STATES(1)+100.600) - CONSTS(47))/8.28100)+ 0.305200*EXP(((STATES(1)+0.994100) - CONSTS(47))/38.4500))
      RATES(14) = (ALGBRC(17) - STATES(14))/ALGBRC(28)
      ALGBRC(32) = 1.00000/(1.00000+EXP(- (STATES(1) - 24.3400)/14.8200))
      RATES(23) = (ALGBRC(32) - STATES(23))/ALGBRC(18)
      ALGBRC(33) = 7.00000+1.00000/( 0.0400000*EXP(- (STATES(1) - 4.00000)/7.00000)+ 0.0400000*EXP((STATES(1) - 4.00000)/7.00000))
      RATES(29) = (ALGBRC(20) - STATES(29))/ALGBRC(33)
      ALGBRC(34) = 100.000+1.00000/( 0.000120000*EXP(- STATES(1)/3.00000)+ 0.000120000*EXP(STATES(1)/7.00000))
      RATES(30) = (ALGBRC(20) - STATES(30))/ALGBRC(34)
      ALGBRC(35) =  2.50000*ALGBRC(23)
      RATES(32) = (ALGBRC(9) - STATES(32))/ALGBRC(35)
      ALGBRC(25) = ALGBRC(11)
      ALGBRC(36) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      RATES(46) = (ALGBRC(25) - STATES(46))/ALGBRC(36)
      ALGBRC(29) = 1.00000/(1.00000+EXP(((STATES(1)+89.1000) - CONSTS(47))/6.08600))
      ALGBRC(38) =  3.00000*ALGBRC(16)
      RATES(15) = (ALGBRC(29) - STATES(15))/ALGBRC(38)
      ALGBRC(39) =  1.46000*ALGBRC(28)
      RATES(16) = (ALGBRC(17) - STATES(16))/ALGBRC(39)
      ALGBRC(30) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(40) = ALGBRC(14)
      RATES(17) = (ALGBRC(30) - STATES(17))/ALGBRC(40)
      ALGBRC(42) =  2.50000*ALGBRC(33)
      RATES(33) = (ALGBRC(20) - STATES(33))/ALGBRC(42)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+43.9400)/5.71100))
      ALGBRC(19) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+70.0000)/5.00000)), 1.00000)
      ALGBRC(31) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+50.0000)/16.5900))
      ALGBRC(44) =  ALGBRC(31)*ALGBRC(19)
      RATES(21) = (ALGBRC(7) - STATES(21))/ALGBRC(44)
      ALGBRC(41) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
      ALGBRC(46) =  ALGBRC(41)*ALGBRC(19)
      RATES(22) = (ALGBRC(7) - STATES(22))/ALGBRC(46)
      ALGBRC(48) = 1.35400+0.000100000/(EXP((STATES(1) - 167.400)/15.8900)+EXP(- (STATES(1) - 12.2300)/0.215400))
      ALGBRC(50) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+70.0000)/20.0000))
      ALGBRC(52) =  ALGBRC(48)*ALGBRC(50)*ALGBRC(44)
      RATES(24) = (ALGBRC(7) - STATES(24))/ALGBRC(52)
      ALGBRC(53) =  ALGBRC(48)*ALGBRC(50)*ALGBRC(46)
      RATES(25) = (ALGBRC(7) - STATES(25))/ALGBRC(53)
      ALGBRC(68) =  CONSTS(156)*STATES(27)+ CONSTS(171)*STATES(28)
      ALGBRC(69) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(70) = 1.00000 - ALGBRC(69)
      ALGBRC(71) =  ALGBRC(69)*STATES(29)+ ALGBRC(70)*STATES(30)
      ALGBRC(72) =  CONSTS(156)*STATES(32)+ CONSTS(171)*STATES(28)
      ALGBRC(73) =  ALGBRC(69)*STATES(33)+ ALGBRC(70)*STATES(30)
      ALGBRC(13) =  STATES(1)*CONSTS(150)
      ALGBRC(74) = ( 4.00000*CONSTS(170)*( STATES(3)*EXP( 2.00000*ALGBRC(13)) -  0.341000*CONSTS(3)))/CONSTS(175)
      ALGBRC(75) =  CONSTS(175)*(STATES(1) - CONSTS(159))
      ALGBRC(76) = TERNRY(- 1.00000e-07.LE.ALGBRC(75).AND.ALGBRC(75).LE.1.00000e-07,  ALGBRC(74)*(1.00000 -  0.500000*ALGBRC(75)), ( ALGBRC(74)*ALGBRC(75))/(EXP(ALGBRC(75)) - 1.00000))
      ALGBRC(43) = ALGBRC(37)+STATES(2)
      ALGBRC(83) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(84) =  (1.00000 - ALGBRC(83))*CONSTS(157)*ALGBRC(76)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(172)*ALGBRC(76)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      ALGBRC(86) = ( CONSTS(167)*- ALGBRC(84))/(1.00000+ 1.00000*1.50000/STATES(9) ** 8.00000)
      ALGBRC(89) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(86)*1.70000, ALGBRC(86))
      ALGBRC(92) = CONSTS(147)/(1.00000+0.0123000/STATES(9))
      ALGBRC(95) = TERNRY(ALGBRC(92).LT.0.00100000, 0.00100000, ALGBRC(92))
      RATES(48) = (ALGBRC(89) - STATES(48))/ALGBRC(95)
      ALGBRC(87) = ( CONSTS(181)*- ALGBRC(84))/(1.00000+1.50000/STATES(9) ** 8.00000)
      ALGBRC(90) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(87)*1.70000, ALGBRC(87))
      ALGBRC(93) = CONSTS(168)/(1.00000+0.0123000/STATES(9))
      ALGBRC(96) = TERNRY(ALGBRC(93).LT.0.00100000, 0.00100000, ALGBRC(93))
      RATES(49) = (ALGBRC(90) - STATES(49))/ALGBRC(96)
      ALGBRC(54) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(4)/STATES(6))
      ALGBRC(62) = 1.00000/(1.00000+EXP((STATES(1) - 213.600)/151.200))
      ALGBRC(63) = 1.00000 - ALGBRC(62)
      ALGBRC(64) =  ALGBRC(62)*STATES(21)+ ALGBRC(63)*STATES(22)
      ALGBRC(65) =  ALGBRC(62)*STATES(24)+ ALGBRC(63)*STATES(25)
      ALGBRC(66) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(67) =  CONSTS(155)*(STATES(1) - ALGBRC(54))*( (1.00000 - ALGBRC(66))*STATES(20)*ALGBRC(64)+ ALGBRC(66)*STATES(23)*ALGBRC(65))
      ALGBRC(91) =  CONSTS(160)* (CONSTS(4)/5.40000) ** (1.0 / 2)*STATES(39)*(STATES(1) - ALGBRC(54))
      ALGBRC(55) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log((CONSTS(4)+ CONSTS(34)*CONSTS(2))/(STATES(6)+ CONSTS(34)*STATES(4)))
      ALGBRC(94) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
      ALGBRC(97) =  CONSTS(161)*ALGBRC(94)*STATES(45)*STATES(46)*(STATES(1) - ALGBRC(55))
      ALGBRC(98) = 1.00000/(1.00000+EXP(((STATES(1)+105.800) -  2.60000*CONSTS(4))/9.49300))
      ALGBRC(99) =  CONSTS(162)* CONSTS(4) ** (1.0 / 2)*ALGBRC(98)*STATES(47)*(STATES(1) - ALGBRC(54))
      ALGBRC(163) =  CONSTS(129)*EXP(( (1.00000 - CONSTS(130))*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
      ALGBRC(167) = ( CONSTS(124)*CONSTS(4)/CONSTS(132) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(163) ** 3.00000+1.00000+CONSTS(4)/CONSTS(132) ** 2.00000) - 1.00000)
      ALGBRC(164) = CONSTS(137)/(1.00000+CONSTS(136)/CONSTS(138)+STATES(4)/CONSTS(139)+STATES(6)/CONSTS(140))
      ALGBRC(168) = ( CONSTS(125)*ALGBRC(164)*CONSTS(136))/(1.00000+CONSTS(134)/CONSTS(135))
      ALGBRC(162) =  CONSTS(128)*EXP(( CONSTS(130)*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
      ALGBRC(165) = ( CONSTS(120)*STATES(4)/ALGBRC(162) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(162) ** 3.00000+1.00000+STATES(6)/CONSTS(131) ** 2.00000) - 1.00000)
      ALGBRC(166) = ( CONSTS(123)*CONSTS(2)/ALGBRC(163) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(163) ** 3.00000+1.00000+CONSTS(4)/CONSTS(132) ** 2.00000) - 1.00000)
      ALGBRC(169) = ( CONSTS(127)*STATES(6)/CONSTS(131) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(162) ** 3.00000+1.00000+STATES(6)/CONSTS(131) ** 2.00000) - 1.00000)
      ALGBRC(170) =  CONSTS(204)*ALGBRC(165)*CONSTS(203)+ ALGBRC(166)*ALGBRC(169)*ALGBRC(168)+ CONSTS(203)*ALGBRC(169)*ALGBRC(168)+ ALGBRC(168)*ALGBRC(165)*CONSTS(203)
      ALGBRC(171) =  ALGBRC(166)*CONSTS(202)*ALGBRC(169)+ ALGBRC(165)*CONSTS(203)*ALGBRC(167)+ ALGBRC(167)*CONSTS(202)*ALGBRC(169)+ CONSTS(203)*ALGBRC(167)*ALGBRC(169)
      ALGBRC(172) =  CONSTS(203)*ALGBRC(167)*CONSTS(204)+ ALGBRC(168)*ALGBRC(166)*CONSTS(202)+ ALGBRC(166)*CONSTS(202)*CONSTS(204)+ ALGBRC(167)*CONSTS(204)*CONSTS(202)
      ALGBRC(173) =  ALGBRC(169)*ALGBRC(168)*ALGBRC(166)+ ALGBRC(167)*CONSTS(204)*ALGBRC(165)+ ALGBRC(166)*CONSTS(204)*ALGBRC(165)+ ALGBRC(168)*ALGBRC(166)*ALGBRC(165)
      ALGBRC(174) = ALGBRC(170)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(175) = ALGBRC(171)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(178) =  3.00000*( ALGBRC(174)*ALGBRC(167) -  ALGBRC(175)*ALGBRC(168))
      ALGBRC(176) = ALGBRC(172)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(177) = ALGBRC(173)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(179) =  2.00000*( ALGBRC(177)*CONSTS(202) -  ALGBRC(176)*ALGBRC(165))
      ALGBRC(180) =  CONSTS(205)*( CONSTS(8)*ALGBRC(178)+ CONSTS(10)*ALGBRC(179))
      ALGBRC(181) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
      ALGBRC(182) =  CONSTS(164)*ALGBRC(181)*(STATES(1) - ALGBRC(54))
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(13).AND.VOI.LE.CONSTS(14).AND.(VOI - CONSTS(13)) -  INT((VOI - CONSTS(13))/CONSTS(16))*CONSTS(16).LE.CONSTS(17), CONSTS(15), 0.00000)
      ALGBRC(184) = (STATES(7) - STATES(6))/2.00000
      RATES(6) = ( - ((ALGBRC(67)+ALGBRC(91)+ALGBRC(97)+ALGBRC(99)+ALGBRC(182)+ALGBRC(1)) -  2.00000*ALGBRC(180))*CONSTS(33)*CONSTS(184))/( CONSTS(7)*CONSTS(185))+( ALGBRC(184)*CONSTS(188))/CONSTS(185)
      ALGBRC(80) = ( 0.750000*CONSTS(170)*( STATES(7)*EXP(ALGBRC(13)) - CONSTS(4)))/CONSTS(177)
      ALGBRC(81) =  CONSTS(177)*(STATES(1) - CONSTS(159))
      ALGBRC(82) = TERNRY(- 1.00000e-07.LE.ALGBRC(81).AND.ALGBRC(81).LE.1.00000e-07,  ALGBRC(80)*(1.00000 -  0.500000*ALGBRC(81)), ( ALGBRC(80)*ALGBRC(81))/(EXP(ALGBRC(81)) - 1.00000))
      ALGBRC(88) =  (1.00000 - ALGBRC(83))*CONSTS(174)*ALGBRC(82)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(183)*ALGBRC(82)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      RATES(7) = ( - ALGBRC(88)*CONSTS(33)*CONSTS(184))/( CONSTS(7)*CONSTS(188)) - ALGBRC(184)
      ALGBRC(51) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(2)/STATES(4))
      ALGBRC(56) =  CONSTS(45)*STATES(12)+ CONSTS(152)*STATES(13)
      ALGBRC(57) =  CONSTS(45)*STATES(12)+ CONSTS(152)*STATES(15)
      ALGBRC(58) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(59) =  CONSTS(46)*(STATES(1) - ALGBRC(51))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(58))*ALGBRC(56)*STATES(14)+ ALGBRC(58)*ALGBRC(57)*STATES(16))
      ALGBRC(60) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(61) =  CONSTS(154)*(STATES(1) - ALGBRC(51))*STATES(17)*( (1.00000 - ALGBRC(60))*STATES(18)+ ALGBRC(60)*STATES(19))
      ALGBRC(128) = 1.00000/(1.00000+CONSTS(118)/STATES(10) ** 2.00000)
      ALGBRC(101) = EXP(( CONSTS(116)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
      ALGBRC(108) = 1.00000+ (CONSTS(2)/CONSTS(109))*(1.00000+1.00000/ALGBRC(101))
      ALGBRC(109) = CONSTS(2)/( CONSTS(109)*ALGBRC(101)*ALGBRC(108))
      ALGBRC(112) =  ALGBRC(109)*CONSTS(113)
      ALGBRC(102) = 1.00000+ (STATES(4)/CONSTS(109))*(1.00000+ALGBRC(101))
      ALGBRC(103) = ( STATES(4)*ALGBRC(101))/( CONSTS(109)*ALGBRC(102))
      ALGBRC(115) =  ALGBRC(103)*CONSTS(113)
      ALGBRC(105) = 1.00000+ (STATES(4)/CONSTS(107))*(1.00000+STATES(4)/CONSTS(108))
      ALGBRC(106) = ( STATES(4)*STATES(4))/( ALGBRC(105)*CONSTS(107)*CONSTS(108))
      ALGBRC(118) =  ALGBRC(106)*ALGBRC(103)*CONSTS(111)
      ALGBRC(119) =  ALGBRC(109)*CONSTS(190)*CONSTS(111)
      ALGBRC(110) = 1.00000/ALGBRC(108)
      ALGBRC(111) =  ALGBRC(110)*CONSTS(112)
      ALGBRC(113) = ALGBRC(111)+ALGBRC(112)
      ALGBRC(100) = EXP(( CONSTS(117)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
      ALGBRC(104) = 1.00000/ALGBRC(102)
      ALGBRC(114) = ( ALGBRC(104)*CONSTS(112))/ALGBRC(100)
      ALGBRC(116) = ALGBRC(114)+ALGBRC(115)
      ALGBRC(107) = 1.00000/ALGBRC(105)
      ALGBRC(117) =  ALGBRC(107)*STATES(10)*CONSTS(114)
      ALGBRC(120) =  CONSTS(193)*ALGBRC(116)*(ALGBRC(118)+ALGBRC(117))+ CONSTS(194)*ALGBRC(118)*(CONSTS(193)+ALGBRC(113))
      ALGBRC(121) =  CONSTS(192)*ALGBRC(118)*(ALGBRC(116)+CONSTS(194))+ ALGBRC(116)*ALGBRC(117)*(CONSTS(192)+ALGBRC(119))
      ALGBRC(122) =  CONSTS(192)*ALGBRC(113)*(ALGBRC(118)+ALGBRC(117))+ ALGBRC(119)*ALGBRC(117)*(CONSTS(193)+ALGBRC(113))
      ALGBRC(123) =  CONSTS(193)*ALGBRC(119)*(ALGBRC(116)+CONSTS(194))+ ALGBRC(113)*CONSTS(194)*(CONSTS(192)+ALGBRC(119))
      ALGBRC(124) = ALGBRC(120)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(125) = ALGBRC(121)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(126) = ALGBRC(122)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(127) = ALGBRC(123)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(129) = ( 3.00000*( ALGBRC(127)*ALGBRC(118) -  ALGBRC(124)*ALGBRC(119))+ ALGBRC(126)*ALGBRC(115)) -  ALGBRC(125)*ALGBRC(112)
      ALGBRC(130) =  ALGBRC(125)*CONSTS(193) -  ALGBRC(124)*CONSTS(192)
      ALGBRC(131) =  0.800000*CONSTS(195)*ALGBRC(128)*( CONSTS(8)*ALGBRC(129)+ CONSTS(9)*ALGBRC(130))
      ALGBRC(183) = ( CONSTS(143)*CONSTS(170)*( STATES(4)*EXP(ALGBRC(13)) - CONSTS(2)))/CONSTS(179)
      ALGBRC(185) =  CONSTS(179)*(STATES(1) - CONSTS(165))
      ALGBRC(186) = TERNRY(- 1.00000e-07.LE.ALGBRC(185).AND.ALGBRC(185).LE.1.00000e-07,  ALGBRC(183)*(1.00000 -  0.500000*ALGBRC(185)), ( ALGBRC(183)*ALGBRC(185))/(EXP(ALGBRC(185)) - 1.00000))
      ALGBRC(188) = (STATES(5) - STATES(4))/2.00000
      RATES(4) = ( - (ALGBRC(59)+ALGBRC(61)+ 3.00000*ALGBRC(131)+ 3.00000*ALGBRC(180)+ALGBRC(186))*CONSTS(184)*CONSTS(33))/( CONSTS(7)*CONSTS(185))+( ALGBRC(188)*CONSTS(188))/CONSTS(185)
      ALGBRC(77) = ( 0.750000*CONSTS(170)*( STATES(5)*EXP(ALGBRC(13)) - CONSTS(2)))/CONSTS(176)
      ALGBRC(78) =  CONSTS(176)*(STATES(1) - CONSTS(159))
      ALGBRC(79) = TERNRY(- 1.00000e-07.LE.ALGBRC(78).AND.ALGBRC(78).LE.1.00000e-07,  ALGBRC(77)*(1.00000 -  0.500000*ALGBRC(78)), ( ALGBRC(77)*ALGBRC(78))/(EXP(ALGBRC(78)) - 1.00000))
      ALGBRC(85) =  (1.00000 - ALGBRC(83))*CONSTS(173)*ALGBRC(79)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(182)*ALGBRC(79)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      ALGBRC(158) = 1.00000/(1.00000+CONSTS(118)/STATES(3) ** 2.00000)
      ALGBRC(138) = 1.00000+ (CONSTS(2)/CONSTS(109))*(1.00000+1.00000/ALGBRC(101))
      ALGBRC(139) = CONSTS(2)/( CONSTS(109)*ALGBRC(101)*ALGBRC(138))
      ALGBRC(142) =  ALGBRC(139)*CONSTS(113)
      ALGBRC(132) = 1.00000+ (STATES(5)/CONSTS(109))*(1.00000+ALGBRC(101))
      ALGBRC(133) = ( STATES(5)*ALGBRC(101))/( CONSTS(109)*ALGBRC(132))
      ALGBRC(145) =  ALGBRC(133)*CONSTS(113)
      ALGBRC(135) = 1.00000+ (STATES(5)/CONSTS(107))*(1.00000+STATES(5)/CONSTS(108))
      ALGBRC(136) = ( STATES(5)*STATES(5))/( ALGBRC(135)*CONSTS(107)*CONSTS(108))
      ALGBRC(148) =  ALGBRC(136)*ALGBRC(133)*CONSTS(111)
      ALGBRC(149) =  ALGBRC(139)*CONSTS(197)*CONSTS(111)
      ALGBRC(140) = 1.00000/ALGBRC(138)
      ALGBRC(141) =  ALGBRC(140)*CONSTS(112)
      ALGBRC(143) = ALGBRC(141)+ALGBRC(142)
      ALGBRC(134) = 1.00000/ALGBRC(132)
      ALGBRC(144) = ( ALGBRC(134)*CONSTS(112))/ALGBRC(100)
      ALGBRC(146) = ALGBRC(144)+ALGBRC(145)
      ALGBRC(137) = 1.00000/ALGBRC(135)
      ALGBRC(147) =  ALGBRC(137)*STATES(3)*CONSTS(114)
      ALGBRC(150) =  CONSTS(200)*ALGBRC(146)*(ALGBRC(148)+ALGBRC(147))+ CONSTS(201)*ALGBRC(148)*(CONSTS(200)+ALGBRC(143))
      ALGBRC(151) =  CONSTS(199)*ALGBRC(148)*(ALGBRC(146)+CONSTS(201))+ ALGBRC(146)*ALGBRC(147)*(CONSTS(199)+ALGBRC(149))
      ALGBRC(152) =  CONSTS(199)*ALGBRC(143)*(ALGBRC(148)+ALGBRC(147))+ ALGBRC(149)*ALGBRC(147)*(CONSTS(200)+ALGBRC(143))
      ALGBRC(153) =  CONSTS(200)*ALGBRC(149)*(ALGBRC(146)+CONSTS(201))+ ALGBRC(143)*CONSTS(201)*(CONSTS(199)+ALGBRC(149))
      ALGBRC(154) = ALGBRC(150)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(155) = ALGBRC(151)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(156) = ALGBRC(152)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(157) = ALGBRC(153)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(159) = ( 3.00000*( ALGBRC(157)*ALGBRC(148) -  ALGBRC(154)*ALGBRC(149))+ ALGBRC(156)*ALGBRC(145)) -  ALGBRC(155)*ALGBRC(142)
      ALGBRC(160) =  ALGBRC(155)*CONSTS(200) -  ALGBRC(154)*CONSTS(199)
      ALGBRC(161) =  0.200000*CONSTS(195)*ALGBRC(158)*( CONSTS(8)*ALGBRC(159)+ CONSTS(9)*ALGBRC(160))
      RATES(5) = ( - (ALGBRC(85)+ 3.00000*ALGBRC(161))*CONSTS(33)*CONSTS(184))/( CONSTS(7)*CONSTS(188)) - ALGBRC(188)
      ALGBRC(191) = ( CONSTS(145)*STATES(10))/(CONSTS(146)+STATES(10))
      ALGBRC(187) = ( CONSTS(144)*4.00000*CONSTS(170)*( STATES(10)*EXP( 2.00000*ALGBRC(13)) -  0.341000*CONSTS(3)))/CONSTS(180)
      ALGBRC(189) =  CONSTS(180)*(STATES(1) - CONSTS(166))
      ALGBRC(190) = TERNRY(- 1.00000e-07.LE.ALGBRC(189).AND.ALGBRC(189).LE.1.00000e-07,  ALGBRC(187)*(1.00000 -  0.500000*ALGBRC(189)), ( ALGBRC(187)*ALGBRC(189))/(EXP(ALGBRC(189)) - 1.00000))
      RATES(1) = - (ALGBRC(59)+ALGBRC(61)+ALGBRC(67)+ALGBRC(84)+ALGBRC(85)+ALGBRC(88)+ALGBRC(91)+ALGBRC(97)+ALGBRC(99)+ALGBRC(131)+ALGBRC(161)+ALGBRC(180)+ALGBRC(186)+ALGBRC(182)+ALGBRC(191)+ALGBRC(190)+ALGBRC(1))
      ALGBRC(192) = (STATES(3) - STATES(10))/0.200000
      ALGBRC(193) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(194) =  CONSTS(148)*( (1.00000 - ALGBRC(193))*STATES(48)+ ALGBRC(193)*STATES(49))
      ALGBRC(47) = 1.00000/(1.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(3) ** 2.00000+( CONSTS(29)*CONSTS(30))/CONSTS(30)+STATES(3) ** 2.00000)
      RATES(3) =  ALGBRC(47)*((( - (ALGBRC(84) -  2.00000*ALGBRC(161))*CONSTS(33)*CONSTS(184))/( 2.00000*CONSTS(7)*CONSTS(188))+( ALGBRC(194)*CONSTS(187))/CONSTS(188)) - ALGBRC(192))
      ALGBRC(195) = ( CONSTS(169)*0.00437500*STATES(10))/(STATES(10)+0.000920000)
      ALGBRC(196) = ( CONSTS(169)*2.75000*0.00437500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
      ALGBRC(197) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(198) = ( 0.00393750*STATES(8))/15.0000
      ALGBRC(199) =  CONSTS(149)*(( (1.00000 - ALGBRC(197))*ALGBRC(195)+ ALGBRC(197)*ALGBRC(196)) - ALGBRC(198))
      ALGBRC(45) = 1.00000/(1.00000+( CONSTS(151)*CONSTS(24))/CONSTS(24)+STATES(10) ** 2.00000+( CONSTS(25)*CONSTS(26))/CONSTS(26)+STATES(10) ** 2.00000)
      RATES(10) =  ALGBRC(45)*((( - ((ALGBRC(191)+ALGBRC(190)) -  2.00000*ALGBRC(131))*CONSTS(33)*CONSTS(184))/( 2.00000*CONSTS(7)*CONSTS(185)) - ( ALGBRC(199)*CONSTS(186))/CONSTS(185))+( ALGBRC(192)*CONSTS(188))/CONSTS(185))
      ALGBRC(200) = (STATES(8) - STATES(9))/100.000
      RATES(8) = ALGBRC(199) - ( ALGBRC(200)*CONSTS(187))/CONSTS(186)
      ALGBRC(49) = 1.00000/(1.00000+( CONSTS(31)*CONSTS(32))/CONSTS(32)+STATES(9) ** 2.00000)
      RATES(9) =  ALGBRC(49)*(ALGBRC(200) - ALGBRC(194))
      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(2) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(35))/CONSTS(36)))
      ALGBRC(14) = 1.00000/( CONSTS(39)*EXP((STATES(1)+CONSTS(37))/CONSTS(38))+ CONSTS(40)*EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
      ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(43)) - CONSTS(47))/CONSTS(44)))
      ALGBRC(15) = 1.00000/( 1.43200e-05*EXP(- ((STATES(1)+1.19600) - CONSTS(47))/6.28500)+ 6.14900*EXP(((STATES(1)+0.509600) - CONSTS(47))/20.2700))
      ALGBRC(16) = 1.00000/( 0.00979400*EXP(- ((STATES(1)+17.9500) - CONSTS(47))/28.0500)+ 0.334300*EXP(((STATES(1)+5.73000) - CONSTS(47))/56.6600))
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.3400)/14.8200))
      ALGBRC(18) = 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(8) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000)/4.23000))
      ALGBRC(22) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
      ALGBRC(23) = 7.00000+1.00000/( 0.00450000*EXP(- (STATES(1)+20.0000)/10.0000)+ 0.00450000*EXP((STATES(1)+20.0000)/10.0000))
      ALGBRC(24) = 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(20) = ALGBRC(9)
      ALGBRC(10) =  STATES(31)*1.00000
      ALGBRC(21) = 1.00000/(CONSTS(52)/ALGBRC(10)+1.00000+CONSTS(51)/STATES(3) ** 4.00000)
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
      ALGBRC(26) = CONSTS(105)+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(4)+144.590)/( 1.56920*CONSTS(4)+3.81150)))
      ALGBRC(27) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
      ALGBRC(37) = ( CONSTS(21)*(1.00000 - STATES(2)))/(1.00000+CONSTS(22)/STATES(3))
      ALGBRC(17) = ALGBRC(3)
      ALGBRC(28) = 2.03800+1.00000/( 0.0213600*EXP(- ((STATES(1)+100.600) - CONSTS(47))/8.28100)+ 0.305200*EXP(((STATES(1)+0.994100) - CONSTS(47))/38.4500))
      ALGBRC(32) = 1.00000/(1.00000+EXP(- (STATES(1) - 24.3400)/14.8200))
      ALGBRC(33) = 7.00000+1.00000/( 0.0400000*EXP(- (STATES(1) - 4.00000)/7.00000)+ 0.0400000*EXP((STATES(1) - 4.00000)/7.00000))
      ALGBRC(34) = 100.000+1.00000/( 0.000120000*EXP(- STATES(1)/3.00000)+ 0.000120000*EXP(STATES(1)/7.00000))
      ALGBRC(35) =  2.50000*ALGBRC(23)
      ALGBRC(25) = ALGBRC(11)
      ALGBRC(36) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
      ALGBRC(29) = 1.00000/(1.00000+EXP(((STATES(1)+89.1000) - CONSTS(47))/6.08600))
      ALGBRC(38) =  3.00000*ALGBRC(16)
      ALGBRC(39) =  1.46000*ALGBRC(28)
      ALGBRC(30) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
      ALGBRC(40) = ALGBRC(14)
      ALGBRC(42) =  2.50000*ALGBRC(33)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+43.9400)/5.71100))
      ALGBRC(19) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+70.0000)/5.00000)), 1.00000)
      ALGBRC(31) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+50.0000)/16.5900))
      ALGBRC(44) =  ALGBRC(31)*ALGBRC(19)
      ALGBRC(41) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
      ALGBRC(46) =  ALGBRC(41)*ALGBRC(19)
      ALGBRC(48) = 1.35400+0.000100000/(EXP((STATES(1) - 167.400)/15.8900)+EXP(- (STATES(1) - 12.2300)/0.215400))
      ALGBRC(50) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+70.0000)/20.0000))
      ALGBRC(52) =  ALGBRC(48)*ALGBRC(50)*ALGBRC(44)
      ALGBRC(53) =  ALGBRC(48)*ALGBRC(50)*ALGBRC(46)
      ALGBRC(68) =  CONSTS(156)*STATES(27)+ CONSTS(171)*STATES(28)
      ALGBRC(69) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
      ALGBRC(70) = 1.00000 - ALGBRC(69)
      ALGBRC(71) =  ALGBRC(69)*STATES(29)+ ALGBRC(70)*STATES(30)
      ALGBRC(72) =  CONSTS(156)*STATES(32)+ CONSTS(171)*STATES(28)
      ALGBRC(73) =  ALGBRC(69)*STATES(33)+ ALGBRC(70)*STATES(30)
      ALGBRC(13) =  STATES(1)*CONSTS(150)
      ALGBRC(74) = ( 4.00000*CONSTS(170)*( STATES(3)*EXP( 2.00000*ALGBRC(13)) -  0.341000*CONSTS(3)))/CONSTS(175)
      ALGBRC(75) =  CONSTS(175)*(STATES(1) - CONSTS(159))
      ALGBRC(76) = TERNRY(- 1.00000e-07.LE.ALGBRC(75).AND.ALGBRC(75).LE.1.00000e-07,  ALGBRC(74)*(1.00000 -  0.500000*ALGBRC(75)), ( ALGBRC(74)*ALGBRC(75))/(EXP(ALGBRC(75)) - 1.00000))
      ALGBRC(43) = ALGBRC(37)+STATES(2)
      ALGBRC(83) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(84) =  (1.00000 - ALGBRC(83))*CONSTS(157)*ALGBRC(76)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(172)*ALGBRC(76)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      ALGBRC(86) = ( CONSTS(167)*- ALGBRC(84))/(1.00000+ 1.00000*1.50000/STATES(9) ** 8.00000)
      ALGBRC(89) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(86)*1.70000, ALGBRC(86))
      ALGBRC(92) = CONSTS(147)/(1.00000+0.0123000/STATES(9))
      ALGBRC(95) = TERNRY(ALGBRC(92).LT.0.00100000, 0.00100000, ALGBRC(92))
      ALGBRC(87) = ( CONSTS(181)*- ALGBRC(84))/(1.00000+1.50000/STATES(9) ** 8.00000)
      ALGBRC(90) = TERNRY(CONSTS(1).EQ.2.00000,  ALGBRC(87)*1.70000, ALGBRC(87))
      ALGBRC(93) = CONSTS(168)/(1.00000+0.0123000/STATES(9))
      ALGBRC(96) = TERNRY(ALGBRC(93).LT.0.00100000, 0.00100000, ALGBRC(93))
      ALGBRC(54) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(4)/STATES(6))
      ALGBRC(62) = 1.00000/(1.00000+EXP((STATES(1) - 213.600)/151.200))
      ALGBRC(63) = 1.00000 - ALGBRC(62)
      ALGBRC(64) =  ALGBRC(62)*STATES(21)+ ALGBRC(63)*STATES(22)
      ALGBRC(65) =  ALGBRC(62)*STATES(24)+ ALGBRC(63)*STATES(25)
      ALGBRC(66) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(67) =  CONSTS(155)*(STATES(1) - ALGBRC(54))*( (1.00000 - ALGBRC(66))*STATES(20)*ALGBRC(64)+ ALGBRC(66)*STATES(23)*ALGBRC(65))
      ALGBRC(91) =  CONSTS(160)* (CONSTS(4)/5.40000) ** (1.0 / 2)*STATES(39)*(STATES(1) - ALGBRC(54))
      ALGBRC(55) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log((CONSTS(4)+ CONSTS(34)*CONSTS(2))/(STATES(6)+ CONSTS(34)*STATES(4)))
      ALGBRC(94) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
      ALGBRC(97) =  CONSTS(161)*ALGBRC(94)*STATES(45)*STATES(46)*(STATES(1) - ALGBRC(55))
      ALGBRC(98) = 1.00000/(1.00000+EXP(((STATES(1)+105.800) -  2.60000*CONSTS(4))/9.49300))
      ALGBRC(99) =  CONSTS(162)* CONSTS(4) ** (1.0 / 2)*ALGBRC(98)*STATES(47)*(STATES(1) - ALGBRC(54))
      ALGBRC(163) =  CONSTS(129)*EXP(( (1.00000 - CONSTS(130))*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
      ALGBRC(167) = ( CONSTS(124)*CONSTS(4)/CONSTS(132) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(163) ** 3.00000+1.00000+CONSTS(4)/CONSTS(132) ** 2.00000) - 1.00000)
      ALGBRC(164) = CONSTS(137)/(1.00000+CONSTS(136)/CONSTS(138)+STATES(4)/CONSTS(139)+STATES(6)/CONSTS(140))
      ALGBRC(168) = ( CONSTS(125)*ALGBRC(164)*CONSTS(136))/(1.00000+CONSTS(134)/CONSTS(135))
      ALGBRC(162) =  CONSTS(128)*EXP(( CONSTS(130)*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
      ALGBRC(165) = ( CONSTS(120)*STATES(4)/ALGBRC(162) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(162) ** 3.00000+1.00000+STATES(6)/CONSTS(131) ** 2.00000) - 1.00000)
      ALGBRC(166) = ( CONSTS(123)*CONSTS(2)/ALGBRC(163) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(163) ** 3.00000+1.00000+CONSTS(4)/CONSTS(132) ** 2.00000) - 1.00000)
      ALGBRC(169) = ( CONSTS(127)*STATES(6)/CONSTS(131) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(162) ** 3.00000+1.00000+STATES(6)/CONSTS(131) ** 2.00000) - 1.00000)
      ALGBRC(170) =  CONSTS(204)*ALGBRC(165)*CONSTS(203)+ ALGBRC(166)*ALGBRC(169)*ALGBRC(168)+ CONSTS(203)*ALGBRC(169)*ALGBRC(168)+ ALGBRC(168)*ALGBRC(165)*CONSTS(203)
      ALGBRC(171) =  ALGBRC(166)*CONSTS(202)*ALGBRC(169)+ ALGBRC(165)*CONSTS(203)*ALGBRC(167)+ ALGBRC(167)*CONSTS(202)*ALGBRC(169)+ CONSTS(203)*ALGBRC(167)*ALGBRC(169)
      ALGBRC(172) =  CONSTS(203)*ALGBRC(167)*CONSTS(204)+ ALGBRC(168)*ALGBRC(166)*CONSTS(202)+ ALGBRC(166)*CONSTS(202)*CONSTS(204)+ ALGBRC(167)*CONSTS(204)*CONSTS(202)
      ALGBRC(173) =  ALGBRC(169)*ALGBRC(168)*ALGBRC(166)+ ALGBRC(167)*CONSTS(204)*ALGBRC(165)+ ALGBRC(166)*CONSTS(204)*ALGBRC(165)+ ALGBRC(168)*ALGBRC(166)*ALGBRC(165)
      ALGBRC(174) = ALGBRC(170)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(175) = ALGBRC(171)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(178) =  3.00000*( ALGBRC(174)*ALGBRC(167) -  ALGBRC(175)*ALGBRC(168))
      ALGBRC(176) = ALGBRC(172)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(177) = ALGBRC(173)/(ALGBRC(170)+ALGBRC(171)+ALGBRC(172)+ALGBRC(173))
      ALGBRC(179) =  2.00000*( ALGBRC(177)*CONSTS(202) -  ALGBRC(176)*ALGBRC(165))
      ALGBRC(180) =  CONSTS(205)*( CONSTS(8)*ALGBRC(178)+ CONSTS(10)*ALGBRC(179))
      ALGBRC(181) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
      ALGBRC(182) =  CONSTS(164)*ALGBRC(181)*(STATES(1) - ALGBRC(54))
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(13).AND.VOI.LE.CONSTS(14).AND.(VOI - CONSTS(13)) -  INT((VOI - CONSTS(13))/CONSTS(16))*CONSTS(16).LE.CONSTS(17), CONSTS(15), 0.00000)
      ALGBRC(184) = (STATES(7) - STATES(6))/2.00000
      ALGBRC(80) = ( 0.750000*CONSTS(170)*( STATES(7)*EXP(ALGBRC(13)) - CONSTS(4)))/CONSTS(177)
      ALGBRC(81) =  CONSTS(177)*(STATES(1) - CONSTS(159))
      ALGBRC(82) = TERNRY(- 1.00000e-07.LE.ALGBRC(81).AND.ALGBRC(81).LE.1.00000e-07,  ALGBRC(80)*(1.00000 -  0.500000*ALGBRC(81)), ( ALGBRC(80)*ALGBRC(81))/(EXP(ALGBRC(81)) - 1.00000))
      ALGBRC(88) =  (1.00000 - ALGBRC(83))*CONSTS(174)*ALGBRC(82)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(183)*ALGBRC(82)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      ALGBRC(51) =  (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(2)/STATES(4))
      ALGBRC(56) =  CONSTS(45)*STATES(12)+ CONSTS(152)*STATES(13)
      ALGBRC(57) =  CONSTS(45)*STATES(12)+ CONSTS(152)*STATES(15)
      ALGBRC(58) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(59) =  CONSTS(46)*(STATES(1) - ALGBRC(51))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(58))*ALGBRC(56)*STATES(14)+ ALGBRC(58)*ALGBRC(57)*STATES(16))
      ALGBRC(60) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(61) =  CONSTS(154)*(STATES(1) - ALGBRC(51))*STATES(17)*( (1.00000 - ALGBRC(60))*STATES(18)+ ALGBRC(60)*STATES(19))
      ALGBRC(128) = 1.00000/(1.00000+CONSTS(118)/STATES(10) ** 2.00000)
      ALGBRC(101) = EXP(( CONSTS(116)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
      ALGBRC(108) = 1.00000+ (CONSTS(2)/CONSTS(109))*(1.00000+1.00000/ALGBRC(101))
      ALGBRC(109) = CONSTS(2)/( CONSTS(109)*ALGBRC(101)*ALGBRC(108))
      ALGBRC(112) =  ALGBRC(109)*CONSTS(113)
      ALGBRC(102) = 1.00000+ (STATES(4)/CONSTS(109))*(1.00000+ALGBRC(101))
      ALGBRC(103) = ( STATES(4)*ALGBRC(101))/( CONSTS(109)*ALGBRC(102))
      ALGBRC(115) =  ALGBRC(103)*CONSTS(113)
      ALGBRC(105) = 1.00000+ (STATES(4)/CONSTS(107))*(1.00000+STATES(4)/CONSTS(108))
      ALGBRC(106) = ( STATES(4)*STATES(4))/( ALGBRC(105)*CONSTS(107)*CONSTS(108))
      ALGBRC(118) =  ALGBRC(106)*ALGBRC(103)*CONSTS(111)
      ALGBRC(119) =  ALGBRC(109)*CONSTS(190)*CONSTS(111)
      ALGBRC(110) = 1.00000/ALGBRC(108)
      ALGBRC(111) =  ALGBRC(110)*CONSTS(112)
      ALGBRC(113) = ALGBRC(111)+ALGBRC(112)
      ALGBRC(100) = EXP(( CONSTS(117)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
      ALGBRC(104) = 1.00000/ALGBRC(102)
      ALGBRC(114) = ( ALGBRC(104)*CONSTS(112))/ALGBRC(100)
      ALGBRC(116) = ALGBRC(114)+ALGBRC(115)
      ALGBRC(107) = 1.00000/ALGBRC(105)
      ALGBRC(117) =  ALGBRC(107)*STATES(10)*CONSTS(114)
      ALGBRC(120) =  CONSTS(193)*ALGBRC(116)*(ALGBRC(118)+ALGBRC(117))+ CONSTS(194)*ALGBRC(118)*(CONSTS(193)+ALGBRC(113))
      ALGBRC(121) =  CONSTS(192)*ALGBRC(118)*(ALGBRC(116)+CONSTS(194))+ ALGBRC(116)*ALGBRC(117)*(CONSTS(192)+ALGBRC(119))
      ALGBRC(122) =  CONSTS(192)*ALGBRC(113)*(ALGBRC(118)+ALGBRC(117))+ ALGBRC(119)*ALGBRC(117)*(CONSTS(193)+ALGBRC(113))
      ALGBRC(123) =  CONSTS(193)*ALGBRC(119)*(ALGBRC(116)+CONSTS(194))+ ALGBRC(113)*CONSTS(194)*(CONSTS(192)+ALGBRC(119))
      ALGBRC(124) = ALGBRC(120)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(125) = ALGBRC(121)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(126) = ALGBRC(122)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(127) = ALGBRC(123)/(ALGBRC(120)+ALGBRC(121)+ALGBRC(122)+ALGBRC(123))
      ALGBRC(129) = ( 3.00000*( ALGBRC(127)*ALGBRC(118) -  ALGBRC(124)*ALGBRC(119))+ ALGBRC(126)*ALGBRC(115)) -  ALGBRC(125)*ALGBRC(112)
      ALGBRC(130) =  ALGBRC(125)*CONSTS(193) -  ALGBRC(124)*CONSTS(192)
      ALGBRC(131) =  0.800000*CONSTS(195)*ALGBRC(128)*( CONSTS(8)*ALGBRC(129)+ CONSTS(9)*ALGBRC(130))
      ALGBRC(183) = ( CONSTS(143)*CONSTS(170)*( STATES(4)*EXP(ALGBRC(13)) - CONSTS(2)))/CONSTS(179)
      ALGBRC(185) =  CONSTS(179)*(STATES(1) - CONSTS(165))
      ALGBRC(186) = TERNRY(- 1.00000e-07.LE.ALGBRC(185).AND.ALGBRC(185).LE.1.00000e-07,  ALGBRC(183)*(1.00000 -  0.500000*ALGBRC(185)), ( ALGBRC(183)*ALGBRC(185))/(EXP(ALGBRC(185)) - 1.00000))
      ALGBRC(188) = (STATES(5) - STATES(4))/2.00000
      ALGBRC(77) = ( 0.750000*CONSTS(170)*( STATES(5)*EXP(ALGBRC(13)) - CONSTS(2)))/CONSTS(176)
      ALGBRC(78) =  CONSTS(176)*(STATES(1) - CONSTS(159))
      ALGBRC(79) = TERNRY(- 1.00000e-07.LE.ALGBRC(78).AND.ALGBRC(78).LE.1.00000e-07,  ALGBRC(77)*(1.00000 -  0.500000*ALGBRC(78)), ( ALGBRC(77)*ALGBRC(78))/(EXP(ALGBRC(78)) - 1.00000))
      ALGBRC(85) =  (1.00000 - ALGBRC(83))*CONSTS(173)*ALGBRC(79)*STATES(26)*( ALGBRC(68)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(71)*STATES(34))+ ALGBRC(83)*CONSTS(182)*ALGBRC(79)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(73)*STATES(34))
      ALGBRC(158) = 1.00000/(1.00000+CONSTS(118)/STATES(3) ** 2.00000)
      ALGBRC(138) = 1.00000+ (CONSTS(2)/CONSTS(109))*(1.00000+1.00000/ALGBRC(101))
      ALGBRC(139) = CONSTS(2)/( CONSTS(109)*ALGBRC(101)*ALGBRC(138))
      ALGBRC(142) =  ALGBRC(139)*CONSTS(113)
      ALGBRC(132) = 1.00000+ (STATES(5)/CONSTS(109))*(1.00000+ALGBRC(101))
      ALGBRC(133) = ( STATES(5)*ALGBRC(101))/( CONSTS(109)*ALGBRC(132))
      ALGBRC(145) =  ALGBRC(133)*CONSTS(113)
      ALGBRC(135) = 1.00000+ (STATES(5)/CONSTS(107))*(1.00000+STATES(5)/CONSTS(108))
      ALGBRC(136) = ( STATES(5)*STATES(5))/( ALGBRC(135)*CONSTS(107)*CONSTS(108))
      ALGBRC(148) =  ALGBRC(136)*ALGBRC(133)*CONSTS(111)
      ALGBRC(149) =  ALGBRC(139)*CONSTS(197)*CONSTS(111)
      ALGBRC(140) = 1.00000/ALGBRC(138)
      ALGBRC(141) =  ALGBRC(140)*CONSTS(112)
      ALGBRC(143) = ALGBRC(141)+ALGBRC(142)
      ALGBRC(134) = 1.00000/ALGBRC(132)
      ALGBRC(144) = ( ALGBRC(134)*CONSTS(112))/ALGBRC(100)
      ALGBRC(146) = ALGBRC(144)+ALGBRC(145)
      ALGBRC(137) = 1.00000/ALGBRC(135)
      ALGBRC(147) =  ALGBRC(137)*STATES(3)*CONSTS(114)
      ALGBRC(150) =  CONSTS(200)*ALGBRC(146)*(ALGBRC(148)+ALGBRC(147))+ CONSTS(201)*ALGBRC(148)*(CONSTS(200)+ALGBRC(143))
      ALGBRC(151) =  CONSTS(199)*ALGBRC(148)*(ALGBRC(146)+CONSTS(201))+ ALGBRC(146)*ALGBRC(147)*(CONSTS(199)+ALGBRC(149))
      ALGBRC(152) =  CONSTS(199)*ALGBRC(143)*(ALGBRC(148)+ALGBRC(147))+ ALGBRC(149)*ALGBRC(147)*(CONSTS(200)+ALGBRC(143))
      ALGBRC(153) =  CONSTS(200)*ALGBRC(149)*(ALGBRC(146)+CONSTS(201))+ ALGBRC(143)*CONSTS(201)*(CONSTS(199)+ALGBRC(149))
      ALGBRC(154) = ALGBRC(150)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(155) = ALGBRC(151)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(156) = ALGBRC(152)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(157) = ALGBRC(153)/(ALGBRC(150)+ALGBRC(151)+ALGBRC(152)+ALGBRC(153))
      ALGBRC(159) = ( 3.00000*( ALGBRC(157)*ALGBRC(148) -  ALGBRC(154)*ALGBRC(149))+ ALGBRC(156)*ALGBRC(145)) -  ALGBRC(155)*ALGBRC(142)
      ALGBRC(160) =  ALGBRC(155)*CONSTS(200) -  ALGBRC(154)*CONSTS(199)
      ALGBRC(161) =  0.200000*CONSTS(195)*ALGBRC(158)*( CONSTS(8)*ALGBRC(159)+ CONSTS(9)*ALGBRC(160))
      ALGBRC(191) = ( CONSTS(145)*STATES(10))/(CONSTS(146)+STATES(10))
      ALGBRC(187) = ( CONSTS(144)*4.00000*CONSTS(170)*( STATES(10)*EXP( 2.00000*ALGBRC(13)) -  0.341000*CONSTS(3)))/CONSTS(180)
      ALGBRC(189) =  CONSTS(180)*(STATES(1) - CONSTS(166))
      ALGBRC(190) = TERNRY(- 1.00000e-07.LE.ALGBRC(189).AND.ALGBRC(189).LE.1.00000e-07,  ALGBRC(187)*(1.00000 -  0.500000*ALGBRC(189)), ( ALGBRC(187)*ALGBRC(189))/(EXP(ALGBRC(189)) - 1.00000))
      ALGBRC(192) = (STATES(3) - STATES(10))/0.200000
      ALGBRC(193) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(194) =  CONSTS(148)*( (1.00000 - ALGBRC(193))*STATES(48)+ ALGBRC(193)*STATES(49))
      ALGBRC(47) = 1.00000/(1.00000+( CONSTS(27)*CONSTS(28))/CONSTS(28)+STATES(3) ** 2.00000+( CONSTS(29)*CONSTS(30))/CONSTS(30)+STATES(3) ** 2.00000)
      ALGBRC(195) = ( CONSTS(169)*0.00437500*STATES(10))/(STATES(10)+0.000920000)
      ALGBRC(196) = ( CONSTS(169)*2.75000*0.00437500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
      ALGBRC(197) = 1.00000/(1.00000+CONSTS(18)/ALGBRC(43))
      ALGBRC(198) = ( 0.00393750*STATES(8))/15.0000
      ALGBRC(199) =  CONSTS(149)*(( (1.00000 - ALGBRC(197))*ALGBRC(195)+ ALGBRC(197)*ALGBRC(196)) - ALGBRC(198))
      ALGBRC(45) = 1.00000/(1.00000+( CONSTS(151)*CONSTS(24))/CONSTS(24)+STATES(10) ** 2.00000+( CONSTS(25)*CONSTS(26))/CONSTS(26)+STATES(10) ** 2.00000)
      ALGBRC(200) = (STATES(8) - STATES(9))/100.000
      ALGBRC(49) = 1.00000/(1.00000+( CONSTS(31)*CONSTS(32))/CONSTS(32)+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
Source
Derived from workspace O'Hara-Rudy-CiPA-v1.0 (2017) at changeset 2c2629a8e2b6.
Collaboration
To begin collaborating on this work, please use your git client and issue this command:
License
The terms of use/license for this work is unspecified.