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 198 entries in the algebraic variable array.
C There are a total of 41 entries in each of the rate and state variable arrays.
C There are a total of 156 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 failtype in component environment (dimensionless).
C CONSTS(3) is nonfail_llimit_factor in component distribution_params (dimensionless).
C CONSTS(4) is nonfail_ulimit_factor in component distribution_params (dimensionless).
C CONSTS(5) is fail_downreg_llimit_factor in component distribution_params (dimensionless).
C CONSTS(6) is fail_downreg_ulimit_factor in component distribution_params (dimensionless).
C CONSTS(7) is fail_upreg_llimit_factor in component distribution_params (dimensionless).
C CONSTS(8) is fail_upreg_ulimit_factor in component distribution_params (dimensionless).
C CONSTS(109) is llimit_factor_faildep in component distribution_params (dimensionless).
C CONSTS(117) is ulimit_factor_faildep in component distribution_params (dimensionless).
C CONSTS(119) is GKr_scal_rv in component random_variation (dimensionless).
C CONSTS(101) is GKs_scal_rv in component random_variation (dimensionless).
C CONSTS(102) is GK1_scal_rv in component random_variation (dimensionless).
C CONSTS(120) is PCa_scal_rv in component random_variation (dimensionless).
C CONSTS(122) is Jup_scal_rv in component random_variation (dimensionless).
C CONSTS(126) is Gncx_b_scal_rv in component random_variation (dimensionless).
C CONSTS(121) is Gto_scal_rv in component random_variation (dimensionless).
C is x in component random_variation (dimensionless).
C CONSTS(9) is nao in component extracellular (millimolar).
C CONSTS(10) is cao in component extracellular (millimolar).
C CONSTS(11) is ko in component extracellular (millimolar).
C CONSTS(12) is R in component physical_constants (joule_per_kilomole_kelvin).
C CONSTS(13) is T in component physical_constants (kelvin).
C CONSTS(14) is F in component physical_constants (coulomb_per_mole).
C CONSTS(15) is zna in component physical_constants (dimensionless).
C CONSTS(16) is zca in component physical_constants (dimensionless).
C CONSTS(17) is zk in component physical_constants (dimensionless).
C CONSTS(18) is L in component cell_geometry (centimeter).
C CONSTS(19) is rad in component cell_geometry (centimeter).
C CONSTS(130) is vcell in component cell_geometry (microliter).
C CONSTS(133) is Ageo in component cell_geometry (centimeter_squared).
C CONSTS(134) is Acap in component cell_geometry (centimeter_squared).
C CONSTS(135) is vmyo in component cell_geometry (microliter).
C CONSTS(136) is vnsr in component cell_geometry (microliter).
C CONSTS(137) is vjsr in component cell_geometry (microliter).
C CONSTS(138) is vss in component cell_geometry (microliter).
C STATES(1) is v in component membrane (millivolt).
C ALGBRC(30) is vffrt in component membrane (coulomb_per_mole).
C ALGBRC(40) is vfrt in component membrane (dimensionless).
C ALGBRC(63) is INa in component INa (microA_per_microF).
C ALGBRC(65) is INaL in component INaL (microA_per_microF).
C ALGBRC(71) is Ito in component Ito (microA_per_microF).
C ALGBRC(82) is ICaL in component ICaL (microA_per_microF).
C ALGBRC(83) is ICaNa in component ICaL (microA_per_microF).
C ALGBRC(86) is ICaK in component ICaL (microA_per_microF).
C ALGBRC(97) is IKr in component IKr (microA_per_microF).
C ALGBRC(99) is IKs in component IKs (microA_per_microF).
C ALGBRC(101) is IK1 in component IK1 (microA_per_microF).
C ALGBRC(133) is INaCa_i in component INaCa_i (microA_per_microF).
C ALGBRC(163) is INaCa_ss in component INaCa_i (microA_per_microF).
C ALGBRC(182) is INaK in component INaK (microA_per_microF).
C ALGBRC(185) is INab in component INab (microA_per_microF).
C ALGBRC(184) is IKb in component IKb (microA_per_microF).
C ALGBRC(189) is IpCa in component IpCa (microA_per_microF).
C ALGBRC(187) is ICab in component ICab (microA_per_microF).
C ALGBRC(13) is Istim in component membrane (microA_per_microF).
C CONSTS(20) is amp in component membrane (microA_per_microF).
C CONSTS(21) is duration in component membrane (millisecond).
C CONSTS(22) is stimStart in component membrane (millisecond).
C CONSTS(23) is KmCaMK in component CaMK (millimolar).
C CONSTS(24) is aCaMK in component CaMK (per_millimolar_per_millisecond).
C CONSTS(25) is bCaMK in component CaMK (per_millisecond).
C CONSTS(26) is CaMKo in component CaMK (dimensionless).
C CONSTS(27) is KmCaM in component CaMK (millimolar).
C ALGBRC(46) is CaMKb in component CaMK (millimolar).
C ALGBRC(48) 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(28) is cmdnmax_b in component intracellular_ions (millimolar).
C CONSTS(103) is cmdnmax in component intracellular_ions (millimolar).
C CONSTS(29) is kmcmdn 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 BSRmax in component intracellular_ions (millimolar).
C CONSTS(33) is KmBSR in component intracellular_ions (millimolar).
C CONSTS(34) is BSLmax in component intracellular_ions (millimolar).
C CONSTS(35) is KmBSL in component intracellular_ions (millimolar).
C CONSTS(36) is csqnmax in component intracellular_ions (millimolar).
C CONSTS(37) 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(190) is Jdiff in component diff (millimolar_per_millisecond).
C ALGBRC(197) is Jup in component SERCA (millimolar_per_millisecond).
C ALGBRC(186) is JdiffK in component diff (millimolar_per_millisecond).
C ALGBRC(192) is Jrel in component ryr (millimolar_per_millisecond).
C ALGBRC(198) is Jtr in component trans_flux (millimolar_per_millisecond).
C ALGBRC(50) is Bcai in component intracellular_ions (dimensionless).
C ALGBRC(54) is Bcajsr in component intracellular_ions (dimensionless).
C ALGBRC(52) is Bcass in component intracellular_ions (dimensionless).
C CONSTS(38) is cm in component intracellular_ions (microF_per_centimeter_squared).
C CONSTS(39) 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 EKs in component reversal_potentials (millivolt).
C ALGBRC(1) is mss in component INa (dimensionless).
C ALGBRC(14) is tm in component INa (millisecond).
C CONSTS(40) is mssV1 in component INa (millivolt).
C CONSTS(41) is mssV2 in component INa (millivolt).
C CONSTS(42) is mtV1 in component INa (millivolt).
C CONSTS(43) is mtV2 in component INa (millivolt).
C CONSTS(44) is mtD1 in component INa (dimensionless).
C CONSTS(45) is mtD2 in component INa (dimensionless).
C CONSTS(46) is mtV3 in component INa (millivolt).
C CONSTS(47) is mtV4 in component INa (millivolt).
C STATES(11) is m in component INa (dimensionless).
C ALGBRC(2) 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(48) is hssV1 in component INa (millivolt).
C CONSTS(49) is hssV2 in component INa (millivolt).
C CONSTS(104) is Ahs in component INa (dimensionless).
C CONSTS(50) 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(60) is h in component INa (dimensionless).
C CONSTS(51) is GNa in component INa (milliS_per_microF).
C ALGBRC(17) is jss in component INa (dimensionless).
C ALGBRC(31) is tj in component INa (millisecond).
C STATES(14) is j in component INa (dimensionless).
C ALGBRC(32) is hssp in component INa (dimensionless).
C ALGBRC(41) is thsp in component INa (millisecond).
C STATES(15) is hsp in component INa (dimensionless).
C ALGBRC(61) is hp in component INa (dimensionless).
C ALGBRC(42) is tjp in component INa (millisecond).
C STATES(16) is jp in component INa (dimensionless).
C ALGBRC(62) is fINap in component INa (dimensionless).
C ALGBRC(33) is mLss in component INaL (dimensionless).
C ALGBRC(43) is tmL in component INaL (millisecond).
C STATES(17) is mL in component INaL (dimensionless).
C CONSTS(52) is thL in component INaL (millisecond).
C ALGBRC(3) is hLss in component INaL (dimensionless).
C STATES(18) is hL in component INaL (dimensionless).
C ALGBRC(4) is hLssp in component INaL (dimensionless).
C CONSTS(105) is thLp in component INaL (millisecond).
C STATES(19) is hLp in component INaL (dimensionless).
C CONSTS(53) is GNaL_b in component INaL (milliS_per_microF).
C CONSTS(106) is GNaL in component INaL (milliS_per_microF).
C ALGBRC(64) is fINaLp in component INaL (dimensionless).
C CONSTS(54) is Gto_b in component Ito (milliS_per_microF).
C ALGBRC(5) 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(6) is iss in component Ito (dimensionless).
C ALGBRC(19) is delta_epi in component Ito (dimensionless).
C ALGBRC(34) is tiF_b in component Ito (millisecond).
C ALGBRC(44) is tiS_b in component Ito (millisecond).
C ALGBRC(47) is tiF in component Ito (millisecond).
C ALGBRC(49) is tiS in component Ito (millisecond).
C ALGBRC(66) is AiF in component Ito (dimensionless).
C ALGBRC(67) 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(68) is i in component Ito (dimensionless).
C ALGBRC(35) is assp in component Ito (dimensionless).
C STATES(23) is ap in component Ito (dimensionless).
C ALGBRC(51) is dti_develop in component Ito (dimensionless).
C ALGBRC(53) is dti_recover in component Ito (dimensionless).
C ALGBRC(55) is tiFp in component Ito (millisecond).
C ALGBRC(56) 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(69) is ip in component Ito (dimensionless).
C CONSTS(123) is Gto in component Ito (milliS_per_microF).
C ALGBRC(70) is fItop in component Ito (dimensionless).
C CONSTS(55) is Kmn in component ICaL (millimolar).
C CONSTS(56) is k2n in component ICaL (per_millisecond).
C CONSTS(57) is PCa_b in component ICaL (dimensionless).
C ALGBRC(7) is dss in component ICaL (dimensionless).
C STATES(26) is d in component ICaL (dimensionless).
C ALGBRC(8) is fss in component ICaL (dimensionless).
C CONSTS(107) is Aff in component ICaL (dimensionless).
C CONSTS(114) 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(72) is f in component ICaL (dimensionless).
C ALGBRC(20) is fcass in component ICaL (dimensionless).
C ALGBRC(73) is Afcaf in component ICaL (dimensionless).
C ALGBRC(74) 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(75) 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(76) is fp in component ICaL (dimensionless).
C STATES(33) is fcafp in component ICaL (dimensionless).
C ALGBRC(77) is fcap in component ICaL (dimensionless).
C ALGBRC(9) 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(78) is PhiCaL in component ICaL (dimensionless).
C ALGBRC(79) is PhiCaNa in component ICaL (dimensionless).
C ALGBRC(80) is PhiCaK in component ICaL (dimensionless).
C CONSTS(124) is PCa in component ICaL (dimensionless).
C CONSTS(127) is PCap in component ICaL (dimensionless).
C CONSTS(128) is PCaNa in component ICaL (dimensionless).
C CONSTS(129) is PCaK in component ICaL (dimensionless).
C CONSTS(131) is PCaNap in component ICaL (dimensionless).
C CONSTS(132) is PCaKp in component ICaL (dimensionless).
C ALGBRC(81) 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(36) is tfcaf in component ICaL (millisecond).
C ALGBRC(37) is tfcas in component ICaL (millisecond).
C CONSTS(108) is tjca in component ICaL (millisecond).
C ALGBRC(38) is tffp in component ICaL (millisecond).
C ALGBRC(45) is tfcafp in component ICaL (millisecond).
C CONSTS(58) is GKr_b in component IKr (milliS_per_microF).
C CONSTS(125) is GKr in component IKr (milliS_per_microF).
C ALGBRC(10) is xrss in component IKr (dimensionless).
C ALGBRC(25) is txrf in component IKr (millisecond).
C ALGBRC(26) is txrs in component IKr (millisecond).
C ALGBRC(89) is Axrf in component IKr (dimensionless).
C ALGBRC(92) is Axrs in component IKr (dimensionless).
C STATES(35) is xrf in component IKr (dimensionless).
C STATES(36) is xrs in component IKr (dimensionless).
C ALGBRC(95) is xr in component IKr (dimensionless).
C ALGBRC(96) is rkr in component IKr (dimensionless).
C CONSTS(59) is GKs_b in component IKs (milliS_per_microF).
C CONSTS(115) is GKs in component IKs (milliS_per_microF).
C ALGBRC(11) is xs1ss in component IKs (dimensionless).
C ALGBRC(27) is xs2ss in component IKs (dimensionless).
C ALGBRC(28) is txs1 in component IKs (millisecond).
C STATES(37) is xs1 in component IKs (dimensionless).
C STATES(38) is xs2 in component IKs (dimensionless).
C ALGBRC(98) is KsCa in component IKs (dimensionless).
C ALGBRC(39) is txs2 in component IKs (millisecond).
C CONSTS(116) is GK1 in component IK1 (milliS_per_microF).
C CONSTS(60) is GK1_b in component IK1 (milliS_per_microF).
C ALGBRC(12) is xk1ss in component IK1 (dimensionless).
C ALGBRC(29) is txk1 in component IK1 (millisecond).
C STATES(39) is xk1 in component IK1 (dimensionless).
C ALGBRC(100) is rk1 in component IK1 (millisecond).
C CONSTS(61) is kna1 in component INaCa_i (per_millisecond).
C CONSTS(62) is kna2 in component INaCa_i (per_millisecond).
C CONSTS(63) is kna3 in component INaCa_i (per_millisecond).
C CONSTS(64) is kasymm in component INaCa_i (dimensionless).
C CONSTS(65) is wna in component INaCa_i (dimensionless).
C CONSTS(66) is wca in component INaCa_i (dimensionless).
C CONSTS(67) is wnaca in component INaCa_i (dimensionless).
C CONSTS(68) is kcaon in component INaCa_i (per_millisecond).
C CONSTS(69) is kcaoff in component INaCa_i (per_millisecond).
C CONSTS(70) is qna in component INaCa_i (dimensionless).
C CONSTS(71) is qca in component INaCa_i (dimensionless).
C ALGBRC(103) is hna in component INaCa_i (dimensionless).
C ALGBRC(102) is hca in component INaCa_i (dimensionless).
C CONSTS(72) is KmCaAct in component INaCa_i (millimolar).
C CONSTS(139) is Gncx_b in component INaCa_i (milliS_per_microF).
C CONSTS(146) is Gncx in component INaCa_i (milliS_per_microF).
C ALGBRC(104) is h1_i in component INaCa_i (dimensionless).
C ALGBRC(105) is h2_i in component INaCa_i (dimensionless).
C ALGBRC(106) is h3_i in component INaCa_i (dimensionless).
C ALGBRC(107) is h4_i in component INaCa_i (dimensionless).
C ALGBRC(108) is h5_i in component INaCa_i (dimensionless).
C ALGBRC(109) is h6_i in component INaCa_i (dimensionless).
C ALGBRC(110) is h7_i in component INaCa_i (dimensionless).
C ALGBRC(111) is h8_i in component INaCa_i (dimensionless).
C ALGBRC(112) is h9_i in component INaCa_i (dimensionless).
C CONSTS(140) is h10_i in component INaCa_i (dimensionless).
C CONSTS(141) is h11_i in component INaCa_i (dimensionless).
C CONSTS(142) is h12_i in component INaCa_i (dimensionless).
C CONSTS(143) is k1_i in component INaCa_i (dimensionless).
C CONSTS(144) is k2_i in component INaCa_i (dimensionless).
C ALGBRC(113) is k3p_i in component INaCa_i (dimensionless).
C ALGBRC(114) is k3pp_i in component INaCa_i (dimensionless).
C ALGBRC(115) is k3_i in component INaCa_i (dimensionless).
C ALGBRC(118) is k4_i in component INaCa_i (dimensionless).
C ALGBRC(116) is k4p_i in component INaCa_i (dimensionless).
C ALGBRC(117) is k4pp_i in component INaCa_i (dimensionless).
C CONSTS(145) is k5_i in component INaCa_i (dimensionless).
C ALGBRC(119) is k6_i in component INaCa_i (dimensionless).
C ALGBRC(120) is k7_i in component INaCa_i (dimensionless).
C ALGBRC(121) is k8_i in component INaCa_i (dimensionless).
C ALGBRC(122) is x1_i in component INaCa_i (dimensionless).
C ALGBRC(123) is x2_i in component INaCa_i (dimensionless).
C ALGBRC(124) is x3_i in component INaCa_i (dimensionless).
C ALGBRC(125) is x4_i in component INaCa_i (dimensionless).
C ALGBRC(126) is E1_i in component INaCa_i (dimensionless).
C ALGBRC(127) is E2_i in component INaCa_i (dimensionless).
C ALGBRC(128) is E3_i in component INaCa_i (dimensionless).
C ALGBRC(129) is E4_i in component INaCa_i (dimensionless).
C ALGBRC(130) is allo_i in component INaCa_i (dimensionless).
C ALGBRC(131) is JncxNa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(132) is JncxCa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(134) is h1_ss in component INaCa_i (dimensionless).
C ALGBRC(135) is h2_ss in component INaCa_i (dimensionless).
C ALGBRC(136) is h3_ss in component INaCa_i (dimensionless).
C ALGBRC(137) is h4_ss in component INaCa_i (dimensionless).
C ALGBRC(138) is h5_ss in component INaCa_i (dimensionless).
C ALGBRC(139) is h6_ss in component INaCa_i (dimensionless).
C ALGBRC(140) is h7_ss in component INaCa_i (dimensionless).
C ALGBRC(141) is h8_ss in component INaCa_i (dimensionless).
C ALGBRC(142) is h9_ss in component INaCa_i (dimensionless).
C CONSTS(147) is h10_ss in component INaCa_i (dimensionless).
C CONSTS(148) is h11_ss in component INaCa_i (dimensionless).
C CONSTS(149) is h12_ss in component INaCa_i (dimensionless).
C CONSTS(150) is k1_ss in component INaCa_i (dimensionless).
C CONSTS(151) is k2_ss in component INaCa_i (dimensionless).
C ALGBRC(143) is k3p_ss in component INaCa_i (dimensionless).
C ALGBRC(144) is k3pp_ss in component INaCa_i (dimensionless).
C ALGBRC(145) is k3_ss in component INaCa_i (dimensionless).
C ALGBRC(148) is k4_ss in component INaCa_i (dimensionless).
C ALGBRC(146) is k4p_ss in component INaCa_i (dimensionless).
C ALGBRC(147) is k4pp_ss in component INaCa_i (dimensionless).
C CONSTS(152) is k5_ss in component INaCa_i (dimensionless).
C ALGBRC(149) is k6_ss in component INaCa_i (dimensionless).
C ALGBRC(150) is k7_ss in component INaCa_i (dimensionless).
C ALGBRC(151) is k8_ss in component INaCa_i (dimensionless).
C ALGBRC(152) is x1_ss in component INaCa_i (dimensionless).
C ALGBRC(153) is x2_ss in component INaCa_i (dimensionless).
C ALGBRC(154) is x3_ss in component INaCa_i (dimensionless).
C ALGBRC(155) is x4_ss in component INaCa_i (dimensionless).
C ALGBRC(156) is E1_ss in component INaCa_i (dimensionless).
C ALGBRC(157) is E2_ss in component INaCa_i (dimensionless).
C ALGBRC(158) is E3_ss in component INaCa_i (dimensionless).
C ALGBRC(159) is E4_ss in component INaCa_i (dimensionless).
C ALGBRC(160) is allo_ss in component INaCa_i (dimensionless).
C ALGBRC(161) is JncxNa_ss in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(162) 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_b in component INaK (milliS_per_microF).
C CONSTS(156) is Pnak in component INaK (milliS_per_microF).
C ALGBRC(164) is Knai in component INaK (millimolar).
C ALGBRC(165) is Knao in component INaK (millimolar).
C ALGBRC(166) is P in component INaK (dimensionless).
C ALGBRC(167) is a1 in component INaK (dimensionless).
C CONSTS(153) is b1 in component INaK (dimensionless).
C CONSTS(154) is a2 in component INaK (dimensionless).
C ALGBRC(168) is b2 in component INaK (dimensionless).
C ALGBRC(169) is a3 in component INaK (dimensionless).
C ALGBRC(170) is b3 in component INaK (dimensionless).
C CONSTS(155) is a4 in component INaK (dimensionless).
C ALGBRC(171) is b4 in component INaK (dimensionless).
C ALGBRC(172) is x1 in component INaK (dimensionless).
C ALGBRC(173) is x2 in component INaK (dimensionless).
C ALGBRC(174) is x3 in component INaK (dimensionless).
C ALGBRC(175) is x4 in component INaK (dimensionless).
C ALGBRC(176) is E1 in component INaK (dimensionless).
C ALGBRC(177) is E2 in component INaK (dimensionless).
C ALGBRC(178) is E3 in component INaK (dimensionless).
C ALGBRC(179) is E4 in component INaK (dimensionless).
C ALGBRC(180) is JnakNa in component INaK (millimolar_per_millisecond).
C ALGBRC(181) is JnakK in component INaK (millimolar_per_millisecond).
C ALGBRC(183) is xkb in component IKb (dimensionless).
C CONSTS(95) is GKb_b in component IKb (milliS_per_microF).
C CONSTS(110) is GKb in component IKb (milliS_per_microF).
C CONSTS(96) is PNab in component INab (milliS_per_microF).
C CONSTS(97) is PCab in component ICab (milliS_per_microF).
C CONSTS(98) is GpCa in component IpCa (milliS_per_microF).
C CONSTS(99) is KmCap in component IpCa (millimolar).
C CONSTS(100) is bt in component ryr (millisecond).
C CONSTS(111) is a_rel in component ryr (millisecond).
C ALGBRC(87) is Jrel_inf in component ryr (dimensionless).
C ALGBRC(93) is tau_rel in component ryr (millisecond).
C ALGBRC(88) is Jrel_infp in component ryr (dimensionless).
C ALGBRC(85) is Jrel_temp in component ryr (dimensionless).
C ALGBRC(94) is tau_relp in component ryr (millisecond).
C STATES(40) is Jrelnp in component ryr (dimensionless).
C STATES(41) is Jrelp in component ryr (dimensionless).
C CONSTS(112) is btp in component ryr (millisecond).
C CONSTS(118) is a_relp in component ryr (millisecond).
C ALGBRC(84) is Jrel_inf_temp in component ryr (dimensionless).
C ALGBRC(191) is fJrelp in component ryr (dimensionless).
C ALGBRC(90) is tau_rel_temp in component ryr (millisecond).
C ALGBRC(91) is tau_relp_temp in component ryr (millisecond).
C CONSTS(113) is upScale in component SERCA (dimensionless).
C ALGBRC(193) is Jupnp in component SERCA (millimolar_per_millisecond).
C ALGBRC(194) is Jupp in component SERCA (millimolar_per_millisecond).
C ALGBRC(195) is fJupp in component SERCA (dimensionless).
C ALGBRC(196) is Jleak in component SERCA (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 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 xrf in component IKr (dimensionless).
C RATES(36) is d/dt xrs in component IKr (dimensionless).
C RATES(37) is d/dt xs1 in component IKs (dimensionless).
C RATES(38) is d/dt xs2 in component IKs (dimensionless).
C RATES(39) is d/dt xk1 in component IK1 (dimensionless).
C RATES(40) is d/dt Jrelnp in component ryr (dimensionless).
C RATES(41) is d/dt Jrelp in component ryr (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0
CONSTS(2) = 0
CONSTS(3) = 0.7
CONSTS(4) = 1.3
CONSTS(5) = 0.4
CONSTS(6) = 1
CONSTS(7) = 1
CONSTS(8) = 1.6
CONSTS(9) = 140
CONSTS(10) = 1.8
CONSTS(11) = 5.4
CONSTS(12) = 8314
CONSTS(13) = 310
CONSTS(14) = 96485
CONSTS(15) = 1
CONSTS(16) = 2
CONSTS(17) = 1
CONSTS(18) = 0.01
CONSTS(19) = 0.0011
STATES(1) = -87
CONSTS(20) = -80
CONSTS(21) = 0.5
CONSTS(22) = 50.0
CONSTS(23) = 0.15
CONSTS(24) = 0.05
CONSTS(25) = 0.00068
CONSTS(26) = 0.05
CONSTS(27) = 0.0015
STATES(2) = 0
STATES(3) = 1e-4
CONSTS(28) = 0.05
CONSTS(29) = 0.00238
CONSTS(30) = 0.07
CONSTS(31) = 0.0005
CONSTS(32) = 0.047
CONSTS(33) = 0.00087
CONSTS(34) = 1.124
CONSTS(35) = 0.0087
CONSTS(36) = 10
CONSTS(37) = 0.8
STATES(4) = 7
STATES(5) = 7
STATES(6) = 145
STATES(7) = 145
STATES(8) = 1.2
STATES(9) = 1.2
STATES(10) = 1e-4
CONSTS(38) = 1
CONSTS(39) = 0.01833
CONSTS(40) = 39.57
CONSTS(41) = 9.871
CONSTS(42) = 11.64
CONSTS(43) = 34.77
CONSTS(44) = 6.765
CONSTS(45) = 8.552
CONSTS(46) = 77.42
CONSTS(47) = 5.955
STATES(11) = 0
CONSTS(48) = 82.9
CONSTS(49) = 6.086
CONSTS(50) = 0.99
STATES(12) = 1
STATES(13) = 1
CONSTS(51) = 75
STATES(14) = 1
STATES(15) = 1
STATES(16) = 1
STATES(17) = 0
CONSTS(52) = 200
STATES(18) = 1
STATES(19) = 1
CONSTS(53) = 0.0075
CONSTS(54) = 0.02
STATES(20) = 0
STATES(21) = 1
STATES(22) = 1
STATES(23) = 0
STATES(24) = 1
STATES(25) = 1
CONSTS(55) = 0.002
CONSTS(56) = 1000
CONSTS(57) = 0.0001
STATES(26) = 0
STATES(27) = 1
STATES(28) = 1
STATES(29) = 1
STATES(30) = 1
STATES(31) = 1
STATES(32) = 1
STATES(33) = 1
STATES(34) = 0
CONSTS(58) = 0.046
STATES(35) = 0
STATES(36) = 0
CONSTS(59) = 0.0034
STATES(37) = 0
STATES(38) = 0
CONSTS(60) = 0.1908
STATES(39) = 1
CONSTS(61) = 15
CONSTS(62) = 5
CONSTS(63) = 88.12
CONSTS(64) = 12.5
CONSTS(65) = 6e4
CONSTS(66) = 6e4
CONSTS(67) = 5e3
CONSTS(68) = 1.5e6
CONSTS(69) = 5e3
CONSTS(70) = 0.5224
CONSTS(71) = 0.167
CONSTS(72) = 150e-6
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) = 30
CONSTS(95) = 0.003
CONSTS(96) = 3.75e-10
CONSTS(97) = 2.5e-8
CONSTS(98) = 0.0005
CONSTS(99) = 0.0005
CONSTS(100) = 4.75
STATES(40) = 0
STATES(41) = 0
CONSTS(101) = SampleUsingPDF(&pdf_0, 2, pdf_roots_0, CONSTANTS, ALGEBRAIC)
CONSTS(102) = SampleUsingPDF(&pdf_1, 2, pdf_roots_1, CONSTANTS, ALGEBRAIC)
CONSTS(103) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(28)*1.30000, CONSTS(28))
CONSTS(104) = 1.00000 - CONSTS(50)
CONSTS(105) = 3.00000*CONSTS(52)
CONSTS(106) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(53)*0.600000, CONSTS(53))
CONSTS(107) = 0.600000
CONSTS(108) = 75.0000
CONSTS(109) = TERNRY(CONSTS(2).EQ.-1.00000, CONSTS(5), TERNRY(CONSTS(2).EQ.1.00000, CONSTS(7), CONSTS(3))
CONSTS(110) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(95)*0.600000, CONSTS(95))
CONSTS(111) = 0.500000*CONSTS(100)
CONSTS(112) = 1.25000*CONSTS(100)
CONSTS(113) = TERNRY(CONSTS(1).EQ.1.00000, 1.30000, 1.00000)
CONSTS(114) = 1.00000 - CONSTS(107)
CONSTS(115) = CONSTS(101)*(TERNRY(CONSTS(1).EQ.1.00000, CONSTS(59)*1.40000, CONSTS(59)))
CONSTS(116) = CONSTS(102)*(TERNRY(CONSTS(1).EQ.1.00000, CONSTS(60)*1.20000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(60)*1.30000, CONSTS(60)))
CONSTS(117) = TERNRY(CONSTS(2).EQ.-1.00000, CONSTS(6), TERNRY(CONSTS(2).EQ.1.00000, CONSTS(8), CONSTS(4))
CONSTS(118) = 0.500000*CONSTS(112)
CONSTS(119) = SampleUsingPDF(&pdf_2, 2, pdf_roots_2, CONSTANTS, ALGEBRAIC)
CONSTS(120) = SampleUsingPDF(&pdf_3, 2, pdf_roots_3, CONSTANTS, ALGEBRAIC)
CONSTS(121) = SampleUsingPDF(&pdf_4, 2, pdf_roots_4, CONSTANTS, ALGEBRAIC)
CONSTS(122) = SampleUsingPDF(&pdf_5, 2, pdf_roots_5, CONSTANTS, ALGEBRAIC)
CONSTS(123) = CONSTS(121)*(TERNRY(CONSTS(1).EQ.1.00000, CONSTS(54)*4.00000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(54)*4.00000, CONSTS(54)))
CONSTS(124) = CONSTS(120)*(TERNRY(CONSTS(1).EQ.1.00000, CONSTS(57)*1.20000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(57)*2.50000, CONSTS(57)))
CONSTS(125) = CONSTS(119)*(TERNRY(CONSTS(1).EQ.1.00000, CONSTS(58)*1.30000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(58)*0.800000, CONSTS(58)))
CONSTS(126) = SampleUsingPDF(&pdf_6, 2, pdf_roots_6, CONSTANTS, ALGEBRAIC)
CONSTS(127) = 1.10000*CONSTS(124)
CONSTS(128) = 0.00125000*CONSTS(124)
CONSTS(129) = 0.000357400*CONSTS(124)
CONSTS(130) = 1000.00*3.14000*CONSTS(19)*CONSTS(19)*CONSTS(18)
CONSTS(131) = 0.00125000*CONSTS(127)
CONSTS(132) = 0.000357400*CONSTS(127)
CONSTS(133) = 2.00000*3.14000*CONSTS(19)*CONSTS(19)+ 2.00000*3.14000*CONSTS(19)*CONSTS(18)
CONSTS(134) = 2.00000*CONSTS(133)
CONSTS(135) = 0.680000*CONSTS(130)
CONSTS(136) = 0.0552000*CONSTS(130)
CONSTS(137) = 0.00480000*CONSTS(130)
CONSTS(138) = 0.0200000*CONSTS(130)
CONSTS(139) = CONSTS(126)*0.000800000
CONSTS(140) = CONSTS(64)+1.00000+ (CONSTS(9)/CONSTS(61))*(1.00000+CONSTS(9)/CONSTS(62))
CONSTS(141) = ( CONSTS(9)*CONSTS(9))/( CONSTS(140)*CONSTS(61)*CONSTS(62))
CONSTS(142) = 1.00000/CONSTS(140)
CONSTS(143) = CONSTS(142)*CONSTS(10)*CONSTS(68)
CONSTS(144) = CONSTS(69)
CONSTS(145) = CONSTS(69)
CONSTS(146) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(139)*1.10000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(139)*1.40000, CONSTS(139))
CONSTS(147) = CONSTS(64)+1.00000+ (CONSTS(9)/CONSTS(61))*(1.00000+CONSTS(9)/CONSTS(62))
CONSTS(148) = ( CONSTS(9)*CONSTS(9))/( CONSTS(147)*CONSTS(61)*CONSTS(62))
CONSTS(149) = 1.00000/CONSTS(147)
CONSTS(150) = CONSTS(149)*CONSTS(10)*CONSTS(68)
CONSTS(151) = CONSTS(69)
CONSTS(152) = CONSTS(69)
CONSTS(153) = CONSTS(74)*CONSTS(86)
CONSTS(154) = CONSTS(75)
CONSTS(155) = (( CONSTS(79)*CONSTS(87))/CONSTS(88))/(1.00000+CONSTS(87)/CONSTS(88))
CONSTS(156) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(94)*0.900000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(94)*0.700000, CONSTS(94))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
RATES(18) = (ALGBRC(3) - STATES(18))/CONSTS(52)
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
RATES(19) = (ALGBRC(4) - STATES(19))/CONSTS(105)
ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(40))/CONSTS(41)))
ALGBRC(14) = 1.00000/( CONSTS(44)*EXP((STATES(1)+CONSTS(42))/CONSTS(43))+ CONSTS(45)*EXP(- (STATES(1)+CONSTS(46))/CONSTS(47)))
RATES(11) = (ALGBRC(1) - STATES(11))/ALGBRC(14)
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(48))/CONSTS(49)))
ALGBRC(15) = 1.00000/( 1.43200e-05*EXP(- (STATES(1)+1.19600)/6.28500)+ 6.14900*EXP((STATES(1)+0.509600)/20.2700))
RATES(12) = (ALGBRC(2) - STATES(12))/ALGBRC(15)
ALGBRC(16) = 1.00000/( 0.00979400*EXP(- (STATES(1)+17.9500)/28.0500)+ 0.334300*EXP((STATES(1)+5.73000)/56.6600))
RATES(13) = (ALGBRC(2) - STATES(13))/ALGBRC(16)
ALGBRC(5) = 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(5) - STATES(20))/ALGBRC(18)
ALGBRC(7) = 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(7) - STATES(26))/ALGBRC(22)
ALGBRC(8) = 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(8) - 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(8) - STATES(28))/ALGBRC(24)
ALGBRC(20) = ALGBRC(8)
RATES(31) = (ALGBRC(20) - STATES(31))/CONSTS(108)
ALGBRC(9) = STATES(31)*1.00000
ALGBRC(21) = 1.00000/(CONSTS(56)/ALGBRC(9)+1.00000+CONSTS(55)/STATES(3) ** 4.00000)
RATES(34) = ALGBRC(21)*CONSTS(56) - STATES(34)*ALGBRC(9)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+8.33700)/6.78900))
ALGBRC(25) = 12.9800+1.00000/( 0.365200*EXP((STATES(1) - 31.6600)/3.86900)+ 4.12300e-05*EXP(- (STATES(1) - 47.7800)/20.3800))
RATES(35) = (ALGBRC(10) - STATES(35))/ALGBRC(25)
ALGBRC(26) = 1.86500+1.00000/( 0.0662900*EXP((STATES(1) - 34.7000)/7.35500)+ 1.12800e-05*EXP(- (STATES(1) - 29.7400)/25.9400))
RATES(36) = (ALGBRC(10) - STATES(36))/ALGBRC(26)
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
ALGBRC(28) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
RATES(37) = (ALGBRC(11) - STATES(37))/ALGBRC(28)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(11)+144.590)/( 1.56920*CONSTS(11)+3.81150)))
ALGBRC(29) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
RATES(39) = (ALGBRC(12) - STATES(39))/ALGBRC(29)
ALGBRC(17) = ALGBRC(2)
ALGBRC(31) = 2.03800+1.00000/( 0.0213600*EXP(- (STATES(1)+100.600)/8.28100)+ 0.305200*EXP((STATES(1)+0.994100)/38.4500))
RATES(14) = (ALGBRC(17) - STATES(14))/ALGBRC(31)
ALGBRC(35) = 1.00000/(1.00000+EXP(- (STATES(1) - 24.3400)/14.8200))
RATES(23) = (ALGBRC(35) - STATES(23))/ALGBRC(18)
ALGBRC(36) = 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(36)
ALGBRC(37) = 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(37)
ALGBRC(38) = 2.50000*ALGBRC(23)
RATES(32) = (ALGBRC(8) - STATES(32))/ALGBRC(38)
ALGBRC(27) = ALGBRC(11)
ALGBRC(39) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
RATES(38) = (ALGBRC(27) - STATES(38))/ALGBRC(39)
ALGBRC(46) = ( CONSTS(26)*(1.00000 - STATES(2)))/(1.00000+CONSTS(27)/STATES(3))
RATES(2) = CONSTS(24)*ALGBRC(46)*(ALGBRC(46)+STATES(2)) - CONSTS(25)*STATES(2)
ALGBRC(32) = 1.00000/(1.00000+EXP((STATES(1)+89.1000)/6.08600))
ALGBRC(41) = 3.00000*ALGBRC(16)
RATES(15) = (ALGBRC(32) - STATES(15))/ALGBRC(41)
ALGBRC(42) = 1.46000*ALGBRC(31)
RATES(16) = (ALGBRC(17) - STATES(16))/ALGBRC(42)
ALGBRC(33) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
ALGBRC(43) = ALGBRC(14)
RATES(17) = (ALGBRC(33) - STATES(17))/ALGBRC(43)
ALGBRC(45) = 2.50000*ALGBRC(36)
RATES(33) = (ALGBRC(20) - STATES(33))/ALGBRC(45)
ALGBRC(6) = 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(34) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+50.0000)/16.5900))
ALGBRC(47) = ALGBRC(34)*ALGBRC(19)
RATES(21) = (ALGBRC(6) - STATES(21))/ALGBRC(47)
ALGBRC(44) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
ALGBRC(49) = ALGBRC(44)*ALGBRC(19)
RATES(22) = (ALGBRC(6) - STATES(22))/ALGBRC(49)
ALGBRC(51) = 1.35400+0.000100000/(EXP((STATES(1) - 167.400)/15.8900)+EXP(- (STATES(1) - 12.2300)/0.215400))
ALGBRC(53) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+70.0000)/20.0000))
ALGBRC(55) = ALGBRC(51)*ALGBRC(53)*ALGBRC(47)
RATES(24) = (ALGBRC(6) - STATES(24))/ALGBRC(55)
ALGBRC(56) = ALGBRC(51)*ALGBRC(53)*ALGBRC(49)
RATES(25) = (ALGBRC(6) - STATES(25))/ALGBRC(56)
ALGBRC(72) = CONSTS(107)*STATES(27)+ CONSTS(114)*STATES(28)
ALGBRC(73) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
ALGBRC(74) = 1.00000 - ALGBRC(73)
ALGBRC(75) = ALGBRC(73)*STATES(29)+ ALGBRC(74)*STATES(30)
ALGBRC(76) = CONSTS(107)*STATES(32)+ CONSTS(114)*STATES(28)
ALGBRC(77) = ALGBRC(73)*STATES(33)+ ALGBRC(74)*STATES(30)
ALGBRC(30) = ( STATES(1)*CONSTS(14)*CONSTS(14))/( CONSTS(12)*CONSTS(13))
ALGBRC(40) = ( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))
ALGBRC(78) = ( 4.00000*ALGBRC(30)*( STATES(3)*EXP( 2.00000*ALGBRC(40)) - 0.341000*CONSTS(10)))/(EXP( 2.00000*ALGBRC(40)) - 1.00000)
ALGBRC(48) = ALGBRC(46)+STATES(2)
ALGBRC(81) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(82) = (1.00000 - ALGBRC(81))*CONSTS(124)*ALGBRC(78)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(127)*ALGBRC(78)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
ALGBRC(84) = ( CONSTS(111)*- ALGBRC(82))/(1.00000+ 1.00000*1.50000/STATES(9) ** 8.00000)
ALGBRC(87) = TERNRY(CONSTS(1).EQ.2.00000, ALGBRC(84)*1.70000, ALGBRC(84))
ALGBRC(90) = CONSTS(100)/(1.00000+0.0123000/STATES(9))
ALGBRC(93) = TERNRY(ALGBRC(90).LT.0.00100000, 0.00100000, ALGBRC(90))
RATES(40) = (ALGBRC(87) - STATES(40))/ALGBRC(93)
ALGBRC(85) = ( CONSTS(118)*- ALGBRC(82))/(1.00000+1.50000/STATES(9) ** 8.00000)
ALGBRC(88) = TERNRY(CONSTS(1).EQ.2.00000, ALGBRC(85)*1.70000, ALGBRC(85))
ALGBRC(91) = CONSTS(112)/(1.00000+0.0123000/STATES(9))
ALGBRC(94) = TERNRY(ALGBRC(91).LT.0.00100000, 0.00100000, ALGBRC(91))
RATES(41) = (ALGBRC(88) - STATES(41))/ALGBRC(94)
ALGBRC(58) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(11)/STATES(6))
ALGBRC(66) = 1.00000/(1.00000+EXP((STATES(1) - 213.600)/151.200))
ALGBRC(67) = 1.00000 - ALGBRC(66)
ALGBRC(68) = ALGBRC(66)*STATES(21)+ ALGBRC(67)*STATES(22)
ALGBRC(69) = ALGBRC(66)*STATES(24)+ ALGBRC(67)*STATES(25)
ALGBRC(70) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(71) = CONSTS(123)*(STATES(1) - ALGBRC(58))*( (1.00000 - ALGBRC(70))*STATES(20)*ALGBRC(68)+ ALGBRC(70)*STATES(23)*ALGBRC(69))
ALGBRC(89) = 1.00000/(1.00000+EXP((STATES(1)+54.8100)/38.2100))
ALGBRC(92) = 1.00000 - ALGBRC(89)
ALGBRC(95) = ALGBRC(89)*STATES(35)+ ALGBRC(92)*STATES(36)
ALGBRC(96) = ( (1.00000/(1.00000+EXP((STATES(1)+55.0000)/75.0000)))*1.00000)/(1.00000+EXP((STATES(1) - 10.0000)/30.0000))
ALGBRC(97) = CONSTS(125)* (CONSTS(11)/5.40000) ** (1.0 / 2)*ALGBRC(95)*ALGBRC(96)*(STATES(1) - ALGBRC(58))
ALGBRC(59) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(11)+ CONSTS(39)*CONSTS(9))/(STATES(6)+ CONSTS(39)*STATES(4)))
ALGBRC(98) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
ALGBRC(99) = CONSTS(115)*ALGBRC(98)*STATES(37)*STATES(38)*(STATES(1) - ALGBRC(59))
ALGBRC(100) = 1.00000/(1.00000+EXP(((STATES(1)+105.800) - 2.60000*CONSTS(11))/9.49300))
ALGBRC(101) = CONSTS(116)* CONSTS(11) ** (1.0 / 2)*ALGBRC(100)*STATES(39)*(STATES(1) - ALGBRC(58))
ALGBRC(165) = CONSTS(82)*EXP(( (1.00000 - CONSTS(83))*STATES(1)*CONSTS(14))/( 3.00000*CONSTS(12)*CONSTS(13)))
ALGBRC(169) = ( CONSTS(77)*CONSTS(11)/CONSTS(85) ** 2.00000)/((1.00000+CONSTS(9)/ALGBRC(165) ** 3.00000+1.00000+CONSTS(11)/CONSTS(85) ** 2.00000) - 1.00000)
ALGBRC(166) = CONSTS(90)/(1.00000+CONSTS(89)/CONSTS(91)+STATES(4)/CONSTS(92)+STATES(6)/CONSTS(93))
ALGBRC(170) = ( CONSTS(78)*ALGBRC(166)*CONSTS(89))/(1.00000+CONSTS(87)/CONSTS(88))
ALGBRC(164) = CONSTS(81)*EXP(( CONSTS(83)*STATES(1)*CONSTS(14))/( 3.00000*CONSTS(12)*CONSTS(13)))
ALGBRC(167) = ( CONSTS(73)*STATES(4)/ALGBRC(164) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(164) ** 3.00000+1.00000+STATES(6)/CONSTS(84) ** 2.00000) - 1.00000)
ALGBRC(168) = ( CONSTS(76)*CONSTS(9)/ALGBRC(165) ** 3.00000)/((1.00000+CONSTS(9)/ALGBRC(165) ** 3.00000+1.00000+CONSTS(11)/CONSTS(85) ** 2.00000) - 1.00000)
ALGBRC(171) = ( CONSTS(80)*STATES(6)/CONSTS(84) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(164) ** 3.00000+1.00000+STATES(6)/CONSTS(84) ** 2.00000) - 1.00000)
ALGBRC(172) = CONSTS(155)*ALGBRC(167)*CONSTS(154)+ ALGBRC(168)*ALGBRC(171)*ALGBRC(170)+ CONSTS(154)*ALGBRC(171)*ALGBRC(170)+ ALGBRC(170)*ALGBRC(167)*CONSTS(154)
ALGBRC(173) = ALGBRC(168)*CONSTS(153)*ALGBRC(171)+ ALGBRC(167)*CONSTS(154)*ALGBRC(169)+ ALGBRC(169)*CONSTS(153)*ALGBRC(171)+ CONSTS(154)*ALGBRC(169)*ALGBRC(171)
ALGBRC(174) = CONSTS(154)*ALGBRC(169)*CONSTS(155)+ ALGBRC(170)*ALGBRC(168)*CONSTS(153)+ ALGBRC(168)*CONSTS(153)*CONSTS(155)+ ALGBRC(169)*CONSTS(155)*CONSTS(153)
ALGBRC(175) = ALGBRC(171)*ALGBRC(170)*ALGBRC(168)+ ALGBRC(169)*CONSTS(155)*ALGBRC(167)+ ALGBRC(168)*CONSTS(155)*ALGBRC(167)+ ALGBRC(170)*ALGBRC(168)*ALGBRC(167)
ALGBRC(176) = ALGBRC(172)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(177) = ALGBRC(173)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(180) = 3.00000*( ALGBRC(176)*ALGBRC(169) - ALGBRC(177)*ALGBRC(170))
ALGBRC(178) = ALGBRC(174)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(179) = ALGBRC(175)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(181) = 2.00000*( ALGBRC(179)*CONSTS(153) - ALGBRC(178)*ALGBRC(167))
ALGBRC(182) = CONSTS(156)*( CONSTS(15)*ALGBRC(180)+ CONSTS(17)*ALGBRC(181))
ALGBRC(183) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
ALGBRC(184) = CONSTS(110)*ALGBRC(183)*(STATES(1) - ALGBRC(58))
ALGBRC(13) = TERNRY(VOI.GT.CONSTS(22).AND.VOI.LE.CONSTS(22)+CONSTS(21), CONSTS(20), 0.00000)
ALGBRC(186) = (STATES(7) - STATES(6))/2.00000
RATES(6) = ( - ((ALGBRC(71)+ALGBRC(97)+ALGBRC(99)+ALGBRC(101)+ALGBRC(184)+ALGBRC(13)) - 2.00000*ALGBRC(182))*CONSTS(38)*CONSTS(134))/( CONSTS(14)*CONSTS(135))+( ALGBRC(186)*CONSTS(138))/CONSTS(135)
ALGBRC(80) = ( 1.00000*ALGBRC(30)*( 0.750000*STATES(7)*EXP( 1.00000*ALGBRC(40)) - 0.750000*CONSTS(11)))/(EXP( 1.00000*ALGBRC(40)) - 1.00000)
ALGBRC(86) = (1.00000 - ALGBRC(81))*CONSTS(129)*ALGBRC(80)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(132)*ALGBRC(80)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
RATES(7) = ( - ALGBRC(86)*CONSTS(38)*CONSTS(134))/( CONSTS(14)*CONSTS(138)) - ALGBRC(186)
ALGBRC(57) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(9)/STATES(4))
ALGBRC(60) = CONSTS(50)*STATES(12)+ CONSTS(104)*STATES(13)
ALGBRC(61) = CONSTS(50)*STATES(12)+ CONSTS(104)*STATES(15)
ALGBRC(62) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(63) = CONSTS(51)*(STATES(1) - ALGBRC(57))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(62))*ALGBRC(60)*STATES(14)+ ALGBRC(62)*ALGBRC(61)*STATES(16))
ALGBRC(64) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(65) = CONSTS(106)*(STATES(1) - ALGBRC(57))*STATES(17)*( (1.00000 - ALGBRC(64))*STATES(18)+ ALGBRC(64)*STATES(19))
ALGBRC(130) = 1.00000/(1.00000+CONSTS(72)/STATES(10) ** 2.00000)
ALGBRC(103) = EXP(( CONSTS(70)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))
ALGBRC(110) = 1.00000+ (CONSTS(9)/CONSTS(63))*(1.00000+1.00000/ALGBRC(103))
ALGBRC(111) = CONSTS(9)/( CONSTS(63)*ALGBRC(103)*ALGBRC(110))
ALGBRC(114) = ALGBRC(111)*CONSTS(67)
ALGBRC(104) = 1.00000+ (STATES(4)/CONSTS(63))*(1.00000+ALGBRC(103))
ALGBRC(105) = ( STATES(4)*ALGBRC(103))/( CONSTS(63)*ALGBRC(104))
ALGBRC(117) = ALGBRC(105)*CONSTS(67)
ALGBRC(107) = 1.00000+ (STATES(4)/CONSTS(61))*(1.00000+STATES(4)/CONSTS(62))
ALGBRC(108) = ( STATES(4)*STATES(4))/( ALGBRC(107)*CONSTS(61)*CONSTS(62))
ALGBRC(120) = ALGBRC(108)*ALGBRC(105)*CONSTS(65)
ALGBRC(121) = ALGBRC(111)*CONSTS(141)*CONSTS(65)
ALGBRC(112) = 1.00000/ALGBRC(110)
ALGBRC(113) = ALGBRC(112)*CONSTS(66)
ALGBRC(115) = ALGBRC(113)+ALGBRC(114)
ALGBRC(102) = EXP(( CONSTS(71)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))
ALGBRC(106) = 1.00000/ALGBRC(104)
ALGBRC(116) = ( ALGBRC(106)*CONSTS(66))/ALGBRC(102)
ALGBRC(118) = ALGBRC(116)+ALGBRC(117)
ALGBRC(109) = 1.00000/ALGBRC(107)
ALGBRC(119) = ALGBRC(109)*STATES(10)*CONSTS(68)
ALGBRC(122) = CONSTS(144)*ALGBRC(118)*(ALGBRC(120)+ALGBRC(119))+ CONSTS(145)*ALGBRC(120)*(CONSTS(144)+ALGBRC(115))
ALGBRC(123) = CONSTS(143)*ALGBRC(120)*(ALGBRC(118)+CONSTS(145))+ ALGBRC(118)*ALGBRC(119)*(CONSTS(143)+ALGBRC(121))
ALGBRC(124) = CONSTS(143)*ALGBRC(115)*(ALGBRC(120)+ALGBRC(119))+ ALGBRC(121)*ALGBRC(119)*(CONSTS(144)+ALGBRC(115))
ALGBRC(125) = CONSTS(144)*ALGBRC(121)*(ALGBRC(118)+CONSTS(145))+ ALGBRC(115)*CONSTS(145)*(CONSTS(143)+ALGBRC(121))
ALGBRC(126) = ALGBRC(122)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(127) = ALGBRC(123)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(128) = ALGBRC(124)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(129) = ALGBRC(125)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(131) = ( 3.00000*( ALGBRC(129)*ALGBRC(120) - ALGBRC(126)*ALGBRC(121))+ ALGBRC(128)*ALGBRC(117)) - ALGBRC(127)*ALGBRC(114)
ALGBRC(132) = ALGBRC(127)*CONSTS(144) - ALGBRC(126)*CONSTS(143)
ALGBRC(133) = 0.800000*CONSTS(146)*ALGBRC(130)*( CONSTS(15)*ALGBRC(131)+ CONSTS(16)*ALGBRC(132))
ALGBRC(185) = ( CONSTS(96)*ALGBRC(30)*( STATES(4)*EXP(ALGBRC(40)) - CONSTS(9)))/(EXP(ALGBRC(40)) - 1.00000)
ALGBRC(188) = (STATES(5) - STATES(4))/2.00000
RATES(4) = ( - (ALGBRC(63)+ALGBRC(65)+ 3.00000*ALGBRC(133)+ 3.00000*ALGBRC(182)+ALGBRC(185))*CONSTS(134)*CONSTS(38))/( CONSTS(14)*CONSTS(135))+( ALGBRC(188)*CONSTS(138))/CONSTS(135)
ALGBRC(79) = ( 1.00000*ALGBRC(30)*( 0.750000*STATES(5)*EXP( 1.00000*ALGBRC(40)) - 0.750000*CONSTS(9)))/(EXP( 1.00000*ALGBRC(40)) - 1.00000)
ALGBRC(83) = (1.00000 - ALGBRC(81))*CONSTS(128)*ALGBRC(79)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(131)*ALGBRC(79)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
ALGBRC(160) = 1.00000/(1.00000+CONSTS(72)/STATES(3) ** 2.00000)
ALGBRC(140) = 1.00000+ (CONSTS(9)/CONSTS(63))*(1.00000+1.00000/ALGBRC(103))
ALGBRC(141) = CONSTS(9)/( CONSTS(63)*ALGBRC(103)*ALGBRC(140))
ALGBRC(144) = ALGBRC(141)*CONSTS(67)
ALGBRC(134) = 1.00000+ (STATES(5)/CONSTS(63))*(1.00000+ALGBRC(103))
ALGBRC(135) = ( STATES(5)*ALGBRC(103))/( CONSTS(63)*ALGBRC(134))
ALGBRC(147) = ALGBRC(135)*CONSTS(67)
ALGBRC(137) = 1.00000+ (STATES(5)/CONSTS(61))*(1.00000+STATES(5)/CONSTS(62))
ALGBRC(138) = ( STATES(5)*STATES(5))/( ALGBRC(137)*CONSTS(61)*CONSTS(62))
ALGBRC(150) = ALGBRC(138)*ALGBRC(135)*CONSTS(65)
ALGBRC(151) = ALGBRC(141)*CONSTS(148)*CONSTS(65)
ALGBRC(142) = 1.00000/ALGBRC(140)
ALGBRC(143) = ALGBRC(142)*CONSTS(66)
ALGBRC(145) = ALGBRC(143)+ALGBRC(144)
ALGBRC(136) = 1.00000/ALGBRC(134)
ALGBRC(146) = ( ALGBRC(136)*CONSTS(66))/ALGBRC(102)
ALGBRC(148) = ALGBRC(146)+ALGBRC(147)
ALGBRC(139) = 1.00000/ALGBRC(137)
ALGBRC(149) = ALGBRC(139)*STATES(3)*CONSTS(68)
ALGBRC(152) = CONSTS(151)*ALGBRC(148)*(ALGBRC(150)+ALGBRC(149))+ CONSTS(152)*ALGBRC(150)*(CONSTS(151)+ALGBRC(145))
ALGBRC(153) = CONSTS(150)*ALGBRC(150)*(ALGBRC(148)+CONSTS(152))+ ALGBRC(148)*ALGBRC(149)*(CONSTS(150)+ALGBRC(151))
ALGBRC(154) = CONSTS(150)*ALGBRC(145)*(ALGBRC(150)+ALGBRC(149))+ ALGBRC(151)*ALGBRC(149)*(CONSTS(151)+ALGBRC(145))
ALGBRC(155) = CONSTS(151)*ALGBRC(151)*(ALGBRC(148)+CONSTS(152))+ ALGBRC(145)*CONSTS(152)*(CONSTS(150)+ALGBRC(151))
ALGBRC(156) = ALGBRC(152)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(157) = ALGBRC(153)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(158) = ALGBRC(154)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(159) = ALGBRC(155)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(161) = ( 3.00000*( ALGBRC(159)*ALGBRC(150) - ALGBRC(156)*ALGBRC(151))+ ALGBRC(158)*ALGBRC(147)) - ALGBRC(157)*ALGBRC(144)
ALGBRC(162) = ALGBRC(157)*CONSTS(151) - ALGBRC(156)*CONSTS(150)
ALGBRC(163) = 0.200000*CONSTS(146)*ALGBRC(160)*( CONSTS(15)*ALGBRC(161)+ CONSTS(16)*ALGBRC(162))
RATES(5) = ( - (ALGBRC(83)+ 3.00000*ALGBRC(163))*CONSTS(38)*CONSTS(134))/( CONSTS(14)*CONSTS(138)) - ALGBRC(188)
ALGBRC(189) = ( CONSTS(98)*STATES(10))/(CONSTS(99)+STATES(10))
ALGBRC(187) = ( CONSTS(97)*4.00000*ALGBRC(30)*( STATES(10)*EXP( 2.00000*ALGBRC(40)) - 0.341000*CONSTS(10)))/(EXP( 2.00000*ALGBRC(40)) - 1.00000)
RATES(1) = - (ALGBRC(63)+ALGBRC(65)+ALGBRC(71)+ALGBRC(82)+ALGBRC(83)+ALGBRC(86)+ALGBRC(97)+ALGBRC(99)+ALGBRC(101)+ALGBRC(133)+ALGBRC(163)+ALGBRC(182)+ALGBRC(185)+ALGBRC(184)+ALGBRC(189)+ALGBRC(187)+ALGBRC(13))
ALGBRC(190) = (STATES(3) - STATES(10))/0.200000
ALGBRC(191) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(192) = (1.00000 - ALGBRC(191))*STATES(40)+ ALGBRC(191)*STATES(41)
ALGBRC(52) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(33))/CONSTS(33)+STATES(3) ** 2.00000+( CONSTS(34)*CONSTS(35))/CONSTS(35)+STATES(3) ** 2.00000)
RATES(3) = ALGBRC(52)*((( - (ALGBRC(82) - 2.00000*ALGBRC(163))*CONSTS(38)*CONSTS(134))/( 2.00000*CONSTS(14)*CONSTS(138))+( ALGBRC(192)*CONSTS(137))/CONSTS(138)) - ALGBRC(190))
ALGBRC(193) = ( CONSTS(113)*0.00437500*STATES(10))/(STATES(10)+0.000920000)
ALGBRC(194) = ( CONSTS(113)*2.75000*0.00437500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
ALGBRC(195) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(196) = ( 0.00393750*STATES(8))/15.0000
ALGBRC(197) = CONSTS(122)*(( (1.00000 - ALGBRC(195))*ALGBRC(193)+ ALGBRC(195)*ALGBRC(194)) - ALGBRC(196))
ALGBRC(50) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(29))/CONSTS(29)+STATES(10) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(10) ** 2.00000)
RATES(10) = ALGBRC(50)*((( - ((ALGBRC(189)+ALGBRC(187)) - 2.00000*ALGBRC(133))*CONSTS(38)*CONSTS(134))/( 2.00000*CONSTS(14)*CONSTS(135)) - ( ALGBRC(197)*CONSTS(136))/CONSTS(135))+( ALGBRC(190)*CONSTS(138))/CONSTS(135))
ALGBRC(198) = (STATES(8) - STATES(9))/100.000
RATES(8) = ALGBRC(197) - ( ALGBRC(198)*CONSTS(137))/CONSTS(136)
ALGBRC(54) = 1.00000/(1.00000+( CONSTS(36)*CONSTS(37))/CONSTS(37)+STATES(9) ** 2.00000)
RATES(9) = ALGBRC(54)*(ALGBRC(198) - ALGBRC(192))
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(40))/CONSTS(41)))
ALGBRC(14) = 1.00000/( CONSTS(44)*EXP((STATES(1)+CONSTS(42))/CONSTS(43))+ CONSTS(45)*EXP(- (STATES(1)+CONSTS(46))/CONSTS(47)))
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(48))/CONSTS(49)))
ALGBRC(15) = 1.00000/( 1.43200e-05*EXP(- (STATES(1)+1.19600)/6.28500)+ 6.14900*EXP((STATES(1)+0.509600)/20.2700))
ALGBRC(16) = 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(- (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(7) = 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(8) = 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(8)
ALGBRC(9) = STATES(31)*1.00000
ALGBRC(21) = 1.00000/(CONSTS(56)/ALGBRC(9)+1.00000+CONSTS(55)/STATES(3) ** 4.00000)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+8.33700)/6.78900))
ALGBRC(25) = 12.9800+1.00000/( 0.365200*EXP((STATES(1) - 31.6600)/3.86900)+ 4.12300e-05*EXP(- (STATES(1) - 47.7800)/20.3800))
ALGBRC(26) = 1.86500+1.00000/( 0.0662900*EXP((STATES(1) - 34.7000)/7.35500)+ 1.12800e-05*EXP(- (STATES(1) - 29.7400)/25.9400))
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000)/8.93200))
ALGBRC(28) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800)/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000)/230.000))
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(11)+144.590)/( 1.56920*CONSTS(11)+3.81150)))
ALGBRC(29) = 122.200/(EXP(- (STATES(1)+127.200)/20.3600)+EXP((STATES(1)+236.800)/69.3300))
ALGBRC(17) = ALGBRC(2)
ALGBRC(31) = 2.03800+1.00000/( 0.0213600*EXP(- (STATES(1)+100.600)/8.28100)+ 0.305200*EXP((STATES(1)+0.994100)/38.4500))
ALGBRC(35) = 1.00000/(1.00000+EXP(- (STATES(1) - 24.3400)/14.8200))
ALGBRC(36) = 7.00000+1.00000/( 0.0400000*EXP(- (STATES(1) - 4.00000)/7.00000)+ 0.0400000*EXP((STATES(1) - 4.00000)/7.00000))
ALGBRC(37) = 100.000+1.00000/( 0.000120000*EXP(- STATES(1)/3.00000)+ 0.000120000*EXP(STATES(1)/7.00000))
ALGBRC(38) = 2.50000*ALGBRC(23)
ALGBRC(27) = ALGBRC(11)
ALGBRC(39) = 1.00000/( 0.0100000*EXP((STATES(1) - 50.0000)/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400)/31.0000))
ALGBRC(46) = ( CONSTS(26)*(1.00000 - STATES(2)))/(1.00000+CONSTS(27)/STATES(3))
ALGBRC(32) = 1.00000/(1.00000+EXP((STATES(1)+89.1000)/6.08600))
ALGBRC(41) = 3.00000*ALGBRC(16)
ALGBRC(42) = 1.46000*ALGBRC(31)
ALGBRC(33) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
ALGBRC(43) = ALGBRC(14)
ALGBRC(45) = 2.50000*ALGBRC(36)
ALGBRC(6) = 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(34) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+50.0000)/16.5900))
ALGBRC(47) = ALGBRC(34)*ALGBRC(19)
ALGBRC(44) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+114.100)/8.07900))
ALGBRC(49) = ALGBRC(44)*ALGBRC(19)
ALGBRC(51) = 1.35400+0.000100000/(EXP((STATES(1) - 167.400)/15.8900)+EXP(- (STATES(1) - 12.2300)/0.215400))
ALGBRC(53) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+70.0000)/20.0000))
ALGBRC(55) = ALGBRC(51)*ALGBRC(53)*ALGBRC(47)
ALGBRC(56) = ALGBRC(51)*ALGBRC(53)*ALGBRC(49)
ALGBRC(72) = CONSTS(107)*STATES(27)+ CONSTS(114)*STATES(28)
ALGBRC(73) = 0.300000+0.600000/(1.00000+EXP((STATES(1) - 10.0000)/10.0000))
ALGBRC(74) = 1.00000 - ALGBRC(73)
ALGBRC(75) = ALGBRC(73)*STATES(29)+ ALGBRC(74)*STATES(30)
ALGBRC(76) = CONSTS(107)*STATES(32)+ CONSTS(114)*STATES(28)
ALGBRC(77) = ALGBRC(73)*STATES(33)+ ALGBRC(74)*STATES(30)
ALGBRC(30) = ( STATES(1)*CONSTS(14)*CONSTS(14))/( CONSTS(12)*CONSTS(13))
ALGBRC(40) = ( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))
ALGBRC(78) = ( 4.00000*ALGBRC(30)*( STATES(3)*EXP( 2.00000*ALGBRC(40)) - 0.341000*CONSTS(10)))/(EXP( 2.00000*ALGBRC(40)) - 1.00000)
ALGBRC(48) = ALGBRC(46)+STATES(2)
ALGBRC(81) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(82) = (1.00000 - ALGBRC(81))*CONSTS(124)*ALGBRC(78)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(127)*ALGBRC(78)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
ALGBRC(84) = ( CONSTS(111)*- ALGBRC(82))/(1.00000+ 1.00000*1.50000/STATES(9) ** 8.00000)
ALGBRC(87) = TERNRY(CONSTS(1).EQ.2.00000, ALGBRC(84)*1.70000, ALGBRC(84))
ALGBRC(90) = CONSTS(100)/(1.00000+0.0123000/STATES(9))
ALGBRC(93) = TERNRY(ALGBRC(90).LT.0.00100000, 0.00100000, ALGBRC(90))
ALGBRC(85) = ( CONSTS(118)*- ALGBRC(82))/(1.00000+1.50000/STATES(9) ** 8.00000)
ALGBRC(88) = TERNRY(CONSTS(1).EQ.2.00000, ALGBRC(85)*1.70000, ALGBRC(85))
ALGBRC(91) = CONSTS(112)/(1.00000+0.0123000/STATES(9))
ALGBRC(94) = TERNRY(ALGBRC(91).LT.0.00100000, 0.00100000, ALGBRC(91))
ALGBRC(58) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(11)/STATES(6))
ALGBRC(66) = 1.00000/(1.00000+EXP((STATES(1) - 213.600)/151.200))
ALGBRC(67) = 1.00000 - ALGBRC(66)
ALGBRC(68) = ALGBRC(66)*STATES(21)+ ALGBRC(67)*STATES(22)
ALGBRC(69) = ALGBRC(66)*STATES(24)+ ALGBRC(67)*STATES(25)
ALGBRC(70) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(71) = CONSTS(123)*(STATES(1) - ALGBRC(58))*( (1.00000 - ALGBRC(70))*STATES(20)*ALGBRC(68)+ ALGBRC(70)*STATES(23)*ALGBRC(69))
ALGBRC(89) = 1.00000/(1.00000+EXP((STATES(1)+54.8100)/38.2100))
ALGBRC(92) = 1.00000 - ALGBRC(89)
ALGBRC(95) = ALGBRC(89)*STATES(35)+ ALGBRC(92)*STATES(36)
ALGBRC(96) = ( (1.00000/(1.00000+EXP((STATES(1)+55.0000)/75.0000)))*1.00000)/(1.00000+EXP((STATES(1) - 10.0000)/30.0000))
ALGBRC(97) = CONSTS(125)* (CONSTS(11)/5.40000) ** (1.0 / 2)*ALGBRC(95)*ALGBRC(96)*(STATES(1) - ALGBRC(58))
ALGBRC(59) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(11)+ CONSTS(39)*CONSTS(9))/(STATES(6)+ CONSTS(39)*STATES(4)))
ALGBRC(98) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(10) ** 1.40000)
ALGBRC(99) = CONSTS(115)*ALGBRC(98)*STATES(37)*STATES(38)*(STATES(1) - ALGBRC(59))
ALGBRC(100) = 1.00000/(1.00000+EXP(((STATES(1)+105.800) - 2.60000*CONSTS(11))/9.49300))
ALGBRC(101) = CONSTS(116)* CONSTS(11) ** (1.0 / 2)*ALGBRC(100)*STATES(39)*(STATES(1) - ALGBRC(58))
ALGBRC(165) = CONSTS(82)*EXP(( (1.00000 - CONSTS(83))*STATES(1)*CONSTS(14))/( 3.00000*CONSTS(12)*CONSTS(13)))
ALGBRC(169) = ( CONSTS(77)*CONSTS(11)/CONSTS(85) ** 2.00000)/((1.00000+CONSTS(9)/ALGBRC(165) ** 3.00000+1.00000+CONSTS(11)/CONSTS(85) ** 2.00000) - 1.00000)
ALGBRC(166) = CONSTS(90)/(1.00000+CONSTS(89)/CONSTS(91)+STATES(4)/CONSTS(92)+STATES(6)/CONSTS(93))
ALGBRC(170) = ( CONSTS(78)*ALGBRC(166)*CONSTS(89))/(1.00000+CONSTS(87)/CONSTS(88))
ALGBRC(164) = CONSTS(81)*EXP(( CONSTS(83)*STATES(1)*CONSTS(14))/( 3.00000*CONSTS(12)*CONSTS(13)))
ALGBRC(167) = ( CONSTS(73)*STATES(4)/ALGBRC(164) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(164) ** 3.00000+1.00000+STATES(6)/CONSTS(84) ** 2.00000) - 1.00000)
ALGBRC(168) = ( CONSTS(76)*CONSTS(9)/ALGBRC(165) ** 3.00000)/((1.00000+CONSTS(9)/ALGBRC(165) ** 3.00000+1.00000+CONSTS(11)/CONSTS(85) ** 2.00000) - 1.00000)
ALGBRC(171) = ( CONSTS(80)*STATES(6)/CONSTS(84) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(164) ** 3.00000+1.00000+STATES(6)/CONSTS(84) ** 2.00000) - 1.00000)
ALGBRC(172) = CONSTS(155)*ALGBRC(167)*CONSTS(154)+ ALGBRC(168)*ALGBRC(171)*ALGBRC(170)+ CONSTS(154)*ALGBRC(171)*ALGBRC(170)+ ALGBRC(170)*ALGBRC(167)*CONSTS(154)
ALGBRC(173) = ALGBRC(168)*CONSTS(153)*ALGBRC(171)+ ALGBRC(167)*CONSTS(154)*ALGBRC(169)+ ALGBRC(169)*CONSTS(153)*ALGBRC(171)+ CONSTS(154)*ALGBRC(169)*ALGBRC(171)
ALGBRC(174) = CONSTS(154)*ALGBRC(169)*CONSTS(155)+ ALGBRC(170)*ALGBRC(168)*CONSTS(153)+ ALGBRC(168)*CONSTS(153)*CONSTS(155)+ ALGBRC(169)*CONSTS(155)*CONSTS(153)
ALGBRC(175) = ALGBRC(171)*ALGBRC(170)*ALGBRC(168)+ ALGBRC(169)*CONSTS(155)*ALGBRC(167)+ ALGBRC(168)*CONSTS(155)*ALGBRC(167)+ ALGBRC(170)*ALGBRC(168)*ALGBRC(167)
ALGBRC(176) = ALGBRC(172)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(177) = ALGBRC(173)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(180) = 3.00000*( ALGBRC(176)*ALGBRC(169) - ALGBRC(177)*ALGBRC(170))
ALGBRC(178) = ALGBRC(174)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(179) = ALGBRC(175)/(ALGBRC(172)+ALGBRC(173)+ALGBRC(174)+ALGBRC(175))
ALGBRC(181) = 2.00000*( ALGBRC(179)*CONSTS(153) - ALGBRC(178)*ALGBRC(167))
ALGBRC(182) = CONSTS(156)*( CONSTS(15)*ALGBRC(180)+ CONSTS(17)*ALGBRC(181))
ALGBRC(183) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
ALGBRC(184) = CONSTS(110)*ALGBRC(183)*(STATES(1) - ALGBRC(58))
ALGBRC(13) = TERNRY(VOI.GT.CONSTS(22).AND.VOI.LE.CONSTS(22)+CONSTS(21), CONSTS(20), 0.00000)
ALGBRC(186) = (STATES(7) - STATES(6))/2.00000
ALGBRC(80) = ( 1.00000*ALGBRC(30)*( 0.750000*STATES(7)*EXP( 1.00000*ALGBRC(40)) - 0.750000*CONSTS(11)))/(EXP( 1.00000*ALGBRC(40)) - 1.00000)
ALGBRC(86) = (1.00000 - ALGBRC(81))*CONSTS(129)*ALGBRC(80)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(132)*ALGBRC(80)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
ALGBRC(57) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(9)/STATES(4))
ALGBRC(60) = CONSTS(50)*STATES(12)+ CONSTS(104)*STATES(13)
ALGBRC(61) = CONSTS(50)*STATES(12)+ CONSTS(104)*STATES(15)
ALGBRC(62) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(63) = CONSTS(51)*(STATES(1) - ALGBRC(57))*STATES(11) ** 3.00000*( (1.00000 - ALGBRC(62))*ALGBRC(60)*STATES(14)+ ALGBRC(62)*ALGBRC(61)*STATES(16))
ALGBRC(64) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(65) = CONSTS(106)*(STATES(1) - ALGBRC(57))*STATES(17)*( (1.00000 - ALGBRC(64))*STATES(18)+ ALGBRC(64)*STATES(19))
ALGBRC(130) = 1.00000/(1.00000+CONSTS(72)/STATES(10) ** 2.00000)
ALGBRC(103) = EXP(( CONSTS(70)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))
ALGBRC(110) = 1.00000+ (CONSTS(9)/CONSTS(63))*(1.00000+1.00000/ALGBRC(103))
ALGBRC(111) = CONSTS(9)/( CONSTS(63)*ALGBRC(103)*ALGBRC(110))
ALGBRC(114) = ALGBRC(111)*CONSTS(67)
ALGBRC(104) = 1.00000+ (STATES(4)/CONSTS(63))*(1.00000+ALGBRC(103))
ALGBRC(105) = ( STATES(4)*ALGBRC(103))/( CONSTS(63)*ALGBRC(104))
ALGBRC(117) = ALGBRC(105)*CONSTS(67)
ALGBRC(107) = 1.00000+ (STATES(4)/CONSTS(61))*(1.00000+STATES(4)/CONSTS(62))
ALGBRC(108) = ( STATES(4)*STATES(4))/( ALGBRC(107)*CONSTS(61)*CONSTS(62))
ALGBRC(120) = ALGBRC(108)*ALGBRC(105)*CONSTS(65)
ALGBRC(121) = ALGBRC(111)*CONSTS(141)*CONSTS(65)
ALGBRC(112) = 1.00000/ALGBRC(110)
ALGBRC(113) = ALGBRC(112)*CONSTS(66)
ALGBRC(115) = ALGBRC(113)+ALGBRC(114)
ALGBRC(102) = EXP(( CONSTS(71)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))
ALGBRC(106) = 1.00000/ALGBRC(104)
ALGBRC(116) = ( ALGBRC(106)*CONSTS(66))/ALGBRC(102)
ALGBRC(118) = ALGBRC(116)+ALGBRC(117)
ALGBRC(109) = 1.00000/ALGBRC(107)
ALGBRC(119) = ALGBRC(109)*STATES(10)*CONSTS(68)
ALGBRC(122) = CONSTS(144)*ALGBRC(118)*(ALGBRC(120)+ALGBRC(119))+ CONSTS(145)*ALGBRC(120)*(CONSTS(144)+ALGBRC(115))
ALGBRC(123) = CONSTS(143)*ALGBRC(120)*(ALGBRC(118)+CONSTS(145))+ ALGBRC(118)*ALGBRC(119)*(CONSTS(143)+ALGBRC(121))
ALGBRC(124) = CONSTS(143)*ALGBRC(115)*(ALGBRC(120)+ALGBRC(119))+ ALGBRC(121)*ALGBRC(119)*(CONSTS(144)+ALGBRC(115))
ALGBRC(125) = CONSTS(144)*ALGBRC(121)*(ALGBRC(118)+CONSTS(145))+ ALGBRC(115)*CONSTS(145)*(CONSTS(143)+ALGBRC(121))
ALGBRC(126) = ALGBRC(122)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(127) = ALGBRC(123)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(128) = ALGBRC(124)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(129) = ALGBRC(125)/(ALGBRC(122)+ALGBRC(123)+ALGBRC(124)+ALGBRC(125))
ALGBRC(131) = ( 3.00000*( ALGBRC(129)*ALGBRC(120) - ALGBRC(126)*ALGBRC(121))+ ALGBRC(128)*ALGBRC(117)) - ALGBRC(127)*ALGBRC(114)
ALGBRC(132) = ALGBRC(127)*CONSTS(144) - ALGBRC(126)*CONSTS(143)
ALGBRC(133) = 0.800000*CONSTS(146)*ALGBRC(130)*( CONSTS(15)*ALGBRC(131)+ CONSTS(16)*ALGBRC(132))
ALGBRC(185) = ( CONSTS(96)*ALGBRC(30)*( STATES(4)*EXP(ALGBRC(40)) - CONSTS(9)))/(EXP(ALGBRC(40)) - 1.00000)
ALGBRC(188) = (STATES(5) - STATES(4))/2.00000
ALGBRC(79) = ( 1.00000*ALGBRC(30)*( 0.750000*STATES(5)*EXP( 1.00000*ALGBRC(40)) - 0.750000*CONSTS(9)))/(EXP( 1.00000*ALGBRC(40)) - 1.00000)
ALGBRC(83) = (1.00000 - ALGBRC(81))*CONSTS(128)*ALGBRC(79)*STATES(26)*( ALGBRC(72)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(75)*STATES(34))+ ALGBRC(81)*CONSTS(131)*ALGBRC(79)*STATES(26)*( ALGBRC(76)*(1.00000 - STATES(34))+ STATES(31)*ALGBRC(77)*STATES(34))
ALGBRC(160) = 1.00000/(1.00000+CONSTS(72)/STATES(3) ** 2.00000)
ALGBRC(140) = 1.00000+ (CONSTS(9)/CONSTS(63))*(1.00000+1.00000/ALGBRC(103))
ALGBRC(141) = CONSTS(9)/( CONSTS(63)*ALGBRC(103)*ALGBRC(140))
ALGBRC(144) = ALGBRC(141)*CONSTS(67)
ALGBRC(134) = 1.00000+ (STATES(5)/CONSTS(63))*(1.00000+ALGBRC(103))
ALGBRC(135) = ( STATES(5)*ALGBRC(103))/( CONSTS(63)*ALGBRC(134))
ALGBRC(147) = ALGBRC(135)*CONSTS(67)
ALGBRC(137) = 1.00000+ (STATES(5)/CONSTS(61))*(1.00000+STATES(5)/CONSTS(62))
ALGBRC(138) = ( STATES(5)*STATES(5))/( ALGBRC(137)*CONSTS(61)*CONSTS(62))
ALGBRC(150) = ALGBRC(138)*ALGBRC(135)*CONSTS(65)
ALGBRC(151) = ALGBRC(141)*CONSTS(148)*CONSTS(65)
ALGBRC(142) = 1.00000/ALGBRC(140)
ALGBRC(143) = ALGBRC(142)*CONSTS(66)
ALGBRC(145) = ALGBRC(143)+ALGBRC(144)
ALGBRC(136) = 1.00000/ALGBRC(134)
ALGBRC(146) = ( ALGBRC(136)*CONSTS(66))/ALGBRC(102)
ALGBRC(148) = ALGBRC(146)+ALGBRC(147)
ALGBRC(139) = 1.00000/ALGBRC(137)
ALGBRC(149) = ALGBRC(139)*STATES(3)*CONSTS(68)
ALGBRC(152) = CONSTS(151)*ALGBRC(148)*(ALGBRC(150)+ALGBRC(149))+ CONSTS(152)*ALGBRC(150)*(CONSTS(151)+ALGBRC(145))
ALGBRC(153) = CONSTS(150)*ALGBRC(150)*(ALGBRC(148)+CONSTS(152))+ ALGBRC(148)*ALGBRC(149)*(CONSTS(150)+ALGBRC(151))
ALGBRC(154) = CONSTS(150)*ALGBRC(145)*(ALGBRC(150)+ALGBRC(149))+ ALGBRC(151)*ALGBRC(149)*(CONSTS(151)+ALGBRC(145))
ALGBRC(155) = CONSTS(151)*ALGBRC(151)*(ALGBRC(148)+CONSTS(152))+ ALGBRC(145)*CONSTS(152)*(CONSTS(150)+ALGBRC(151))
ALGBRC(156) = ALGBRC(152)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(157) = ALGBRC(153)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(158) = ALGBRC(154)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(159) = ALGBRC(155)/(ALGBRC(152)+ALGBRC(153)+ALGBRC(154)+ALGBRC(155))
ALGBRC(161) = ( 3.00000*( ALGBRC(159)*ALGBRC(150) - ALGBRC(156)*ALGBRC(151))+ ALGBRC(158)*ALGBRC(147)) - ALGBRC(157)*ALGBRC(144)
ALGBRC(162) = ALGBRC(157)*CONSTS(151) - ALGBRC(156)*CONSTS(150)
ALGBRC(163) = 0.200000*CONSTS(146)*ALGBRC(160)*( CONSTS(15)*ALGBRC(161)+ CONSTS(16)*ALGBRC(162))
ALGBRC(189) = ( CONSTS(98)*STATES(10))/(CONSTS(99)+STATES(10))
ALGBRC(187) = ( CONSTS(97)*4.00000*ALGBRC(30)*( STATES(10)*EXP( 2.00000*ALGBRC(40)) - 0.341000*CONSTS(10)))/(EXP( 2.00000*ALGBRC(40)) - 1.00000)
ALGBRC(190) = (STATES(3) - STATES(10))/0.200000
ALGBRC(191) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(192) = (1.00000 - ALGBRC(191))*STATES(40)+ ALGBRC(191)*STATES(41)
ALGBRC(52) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(33))/CONSTS(33)+STATES(3) ** 2.00000+( CONSTS(34)*CONSTS(35))/CONSTS(35)+STATES(3) ** 2.00000)
ALGBRC(193) = ( CONSTS(113)*0.00437500*STATES(10))/(STATES(10)+0.000920000)
ALGBRC(194) = ( CONSTS(113)*2.75000*0.00437500*STATES(10))/((STATES(10)+0.000920000) - 0.000170000)
ALGBRC(195) = 1.00000/(1.00000+CONSTS(23)/ALGBRC(48))
ALGBRC(196) = ( 0.00393750*STATES(8))/15.0000
ALGBRC(197) = CONSTS(122)*(( (1.00000 - ALGBRC(195))*ALGBRC(193)+ ALGBRC(195)*ALGBRC(194)) - ALGBRC(196))
ALGBRC(50) = 1.00000/(1.00000+( CONSTS(103)*CONSTS(29))/CONSTS(29)+STATES(10) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(10) ** 2.00000)
ALGBRC(198) = (STATES(8) - STATES(9))/100.000
ALGBRC(54) = 1.00000/(1.00000+( CONSTS(36)*CONSTS(37))/CONSTS(37)+STATES(9) ** 2.00000)
RETURN
END
double pdf_0_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(3));
}
double pdf_0_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(4));
}
double pdf_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(3).AND.bvar.LE.CONSTS(4), 1.00000/(CONSTS(4) - CONSTS(3)), 0.00000));
}
double (*pdf_roots_0[])(double bvar, double*, double*) = {pdf_0_root_0,pdf_0_root_1,};
double pdf_1_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(3));
}
double pdf_1_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(4));
}
double pdf_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(3).AND.bvar.LE.CONSTS(4), 1.00000/(CONSTS(4) - CONSTS(3)), 0.00000));
}
double (*pdf_roots_1[])(double bvar, double*, double*) = {pdf_1_root_0,pdf_1_root_1,};
double pdf_2_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(109));
}
double pdf_2_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(117));
}
double pdf_2(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(109).AND.bvar.LE.CONSTS(117), 1.00000/(CONSTS(117) - CONSTS(109)), 0.00000));
}
double (*pdf_roots_2[])(double bvar, double*, double*) = {pdf_2_root_0,pdf_2_root_1,};
double pdf_3_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(109));
}
double pdf_3_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(117));
}
double pdf_3(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(109).AND.bvar.LE.CONSTS(117), 1.00000/(CONSTS(117) - CONSTS(109)), 0.00000));
}
double (*pdf_roots_3[])(double bvar, double*, double*) = {pdf_3_root_0,pdf_3_root_1,};
double pdf_4_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(109));
}
double pdf_4_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(117));
}
double pdf_4(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(109).AND.bvar.LE.CONSTS(117), 1.00000/(CONSTS(117) - CONSTS(109)), 0.00000));
}
double (*pdf_roots_4[])(double bvar, double*, double*) = {pdf_4_root_0,pdf_4_root_1,};
double pdf_5_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(109));
}
double pdf_5_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(117));
}
double pdf_5(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(109).AND.bvar.LE.CONSTS(117), 1.00000/(CONSTS(117) - CONSTS(109)), 0.00000));
}
double (*pdf_roots_5[])(double bvar, double*, double*) = {pdf_5_root_0,pdf_5_root_1,};
double pdf_6_root_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(109));
}
double pdf_6_root_1(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (bvar - CONSTS(117));
}
double pdf_6(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
return (TERNRY(bvar.GT.CONSTS(109).AND.bvar.LE.CONSTS(117), 1.00000/(CONSTS(117) - CONSTS(109)), 0.00000));
}
double (*pdf_roots_6[])(double bvar, double*, double*) = {pdf_6_root_0,pdf_6_root_1,};
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END