C C There are a total of 134 entries in the algebraic variable array. C There are a total of 51 entries in each of the rate and state variable arrays. C There are a total of 388 entries in the constant variable array. C C C VOI is time in component time (dimensionless). C ALGBRC(2) is Istim in component stimulus_current (dimensionless). C CONSTS(1) is stimPeriod in component stimulus_current (dimensionless). C CONSTS(2) is stimDuration in component stimulus_current (dimensionless). C CONSTS(3) is stimCurrent in component stimulus_current (dimensionless). C CONSTS(4) is clamp_ATPi in component parameters (dimensionless). C CONSTS(5) is clamp_ATPi_cyto in component parameters (dimensionless). C CONSTS(6) is Faraday in component parameters (C_per_mmol). C CONSTS(7) is R in component parameters (J_per_mol_per_K). C CONSTS(8) is T in component parameters (K). C CONSTS(9) is LHospitalThreshold in component parameters (mV). C CONSTS(10) is kt_2 in component parameters (dimensionless). C CONSTS(11) is kf_2 in component parameters (dimensionless). C CONSTS(12) is kf_3 in component parameters (dimensionless). C CONSTS(13) is keq in component parameters (dimensionless). C CONSTS(14) is CRT_cyto in component parameters (mM). C CONSTS(15) is CRT_mito in component parameters (mM). C CONSTS(16) is VATPase_cyto in component parameters (dimensionless). C CONSTS(17) is Acap in component parameters (cm2). C CONSTS(18) is AcCoA in component parameters (mM). C CONSTS(19) is aL in component parameters (dimensionless). C CONSTS(20) is b in component parameters (dimensionless). C CONSTS(21) is bL in component parameters (dimensionless). C CONSTS(22) is C_m in component parameters (uF_per_cm2). C CONSTS(23) is Cao in component parameters (mM). C CONSTS(24) is CIK in component parameters (mM). C CONSTS(25) is Cm in component parameters (mM). C CONSTS(26) is CMDNtot in component parameters (mM). C CONSTS(27) is Cmito in component parameters (mM_per_mV). C CONSTS(28) is CPN in component parameters (mM). C CONSTS(29) is CoA in component parameters (mM). C CONSTS(30) is CSQNtot in component parameters (mM). C CONSTS(31) is DpH in component parameters (pH_units). C CONSTS(32) is Dpsio in component parameters (mV). C CONSTS(33) is eta in component parameters (dimensionless). C CONSTS(34) is EtCS in component parameters (mM). C CONSTS(35) is EtID in component parameters (mM). C CONSTS(36) is EtKG in component parameters (mM). C CONSTS(37) is EtMD in component parameters (mM). C CONSTS(38) is EtSDH in component parameters (mM). C CONSTS(39) is FAD in component parameters (mM). C CONSTS(40) is FADH2 in component parameters (mM). C CONSTS(41) is fL in component parameters (per_ms). C CONSTS(42) is fm in component parameters (dimensionless). C CONSTS(43) is fprime in component parameters (per_ms). C CONSTS(44) is g in component parameters (dimensionless). C CONSTS(45) is G_Cab in component parameters (mS_per_cm2). C CONSTS(46) is G_Kp in component parameters (mS_per_cm2). C CONSTS(47) is G_Na in component parameters (mS_per_cm2). C CONSTS(48) is G_Nab in component parameters (mS_per_cm2). C CONSTS(49) is gh in component parameters (mM_per_ms_per_mV). C CONSTS(50) is gL in component parameters (per_ms). C CONSTS(51) is GLU in component parameters (mM). C CONSTS(52) is gprime in component parameters (per_ms). C CONSTS(53) is H in component parameters (mM). C CONSTS(54) is hm in component parameters (dimensionless). C CONSTS(55) is HTRPNtot in component parameters (mM). C CONSTS(56) is ICahalf in component parameters (uA_per_cm2). C CONSTS(57) is INaKmax in component parameters (uA_per_cm2). C CONSTS(58) is IpCamax in component parameters (uA_per_cm2). C CONSTS(59) is KAATeq in component parameters (dimensionless). C CONSTS(60) is KaCa in component parameters (mM). C CONSTS(61) is KACOeq in component parameters (dimensionless). C CONSTS(62) is kact in component parameters (mM). C CONSTS(63) is KADP in component parameters (mM). C CONSTS(64) is kaminus in component parameters (per_ms). C CONSTS(65) is kaplus in component parameters (per_mM4_per_ms). C CONSTS(66) is kbminus in component parameters (per_ms). C CONSTS(67) is kbplus in component parameters (per_mM3_per_ms). C CONSTS(68) is Kca in component parameters (mM). C CONSTS(69) is kcminus in component parameters (per_ms). C CONSTS(70) is kcnsASP in component parameters (per_ms). C CONSTS(71) is kcplus in component parameters (per_ms). C CONSTS(72) is KCS in component parameters (per_ms). C CONSTS(73) is kf1 in component parameters (dimensionless). C CONSTS(74) is kfAAT in component parameters (per_ms). C CONSTS(75) is kfACO in component parameters (per_ms). C CONSTS(76) is Kfb in component parameters (mM). C CONSTS(77) is kfFH in component parameters (per_ms). C CONSTS(78) is KFHeq in component parameters (dimensionless). C CONSTS(79) is kfSL in component parameters (per_mM_per_ms). C CONSTS(80) is kh_1 in component parameters (mM). C CONSTS(81) is kh_2 in component parameters (mM). C CONSTS(82) is Kh1 in component parameters (mM). C CONSTS(83) is Kh2 in component parameters (mM). C CONSTS(84) is Kh3 in component parameters (mM). C CONSTS(85) is Kh4 in component parameters (mM). C CONSTS(86) is khtrpn_minus in component parameters (per_ms). C CONSTS(87) is khtrpn_plus in component parameters (per_mM_per_ms). C CONSTS(88) is Ki_AM in component parameters (mM). C CONSTS(89) is Ki_prime_SR in component parameters (mM). C CONSTS(90) is Ki_SR in component parameters (mM). C CONSTS(91) is Ki1AD_NaK in component parameters (mM). C CONSTS(92) is KiADP_CaP in component parameters (mM). C CONSTS(93) is kIDH in component parameters (per_ms). C CONSTS(94) is KidhNADH in component parameters (mM). C CONSTS(95) is KiFUM in component parameters (mM). C CONSTS(96) is Kioaa in component parameters (mM). C CONSTS(97) is KiOxaa in component parameters (mM). C CONSTS(98) is kKGDH in component parameters (per_ms). C CONSTS(99) is kltrpn_plus in component parameters (per_mM_per_ms). C CONSTS(100) is kltrpn_minus in component parameters (per_ms). C CONSTS(101) is Km1AT_NaK in component parameters (mM). C CONSTS(102) is Km1ATP_CaP in component parameters (mM). C CONSTS(103) is Km2ATP_CaP in component parameters (mM). C CONSTS(104) is KmAcCoA in component parameters (mM). C CONSTS(105) is Kmal in component parameters (mM). C CONSTS(106) is KmATP_AM in component parameters (mM). C CONSTS(107) is KmATP_SR in component parameters (mM). C CONSTS(108) is KmCMDN in component parameters (mM). C CONSTS(109) is KmCa in component parameters (mM). C CONSTS(110) is KmCSQN in component parameters (mM). C CONSTS(111) is kMDH in component parameters (per_ms). C CONSTS(112) is Kmg in component parameters (mM). C CONSTS(113) is KmIDNAD in component parameters (mM). C CONSTS(114) is Kmiso in component parameters (mM). C CONSTS(115) is KmKG in component parameters (mM). C CONSTS(116) is KmKGNAD in component parameters (mM). C CONSTS(117) is KmKo in component parameters (mM). C CONSTS(118) is KmmNAD in component parameters (mM). C CONSTS(119) is KmNa in component parameters (mM). C CONSTS(120) is KmNai in component parameters (mM). C CONSTS(121) is KmnsCa in component parameters (mM). C CONSTS(122) is KmOaa in component parameters (mM). C CONSTS(123) is KmpCa in component parameters (mM). C CONSTS(124) is KmSucc in component parameters (mM). C CONSTS(125) is Kna in component parameters (mM). C CONSTS(126) is kNaCa in component parameters (uA_per_cm2). C CONSTS(127) is Knca in component parameters (mM). C CONSTS(128) is Ko in component parameters (mM). C CONSTS(129) is Koff in component parameters (dimensionless). C CONSTS(130) is Krb in component parameters (mM). C CONSTS(131) is kres in component parameters (dimensionless). C CONSTS(132) is kresf in component parameters (dimensionless). C CONSTS(133) is ksat in component parameters (dimensionless). C CONSTS(134) is kSDH in component parameters (per_ms). C CONSTS(135) is KSLeq in component parameters (dimensionless). C CONSTS(136) is KSR in component parameters (dimensionless). C CONSTS(137) is ktrans in component parameters (mM). C CONSTS(138) is kTrop_pn in component parameters (per_ms). C CONSTS(139) is L in component parameters (dimensionless). C CONSTS(140) is LTRPNtot in component parameters (mM). C CONSTS(141) is mcoop in component parameters (dimensionless). C CONSTS(142) is Mg in component parameters (mM). C CONSTS(143) is n in component parameters (dimensionless). C CONSTS(144) is na in component parameters (dimensionless). C CONSTS(145) is Nao in component parameters (mM). C CONSTS(146) is ncoop in component parameters (dimensionless). C CONSTS(147) is Nfb in component parameters (dimensionless). C CONSTS(148) is nID in component parameters (dimensionless). C CONSTS(149) is nKG in component parameters (dimensionless). C CONSTS(150) is Nrb in component parameters (dimensionless). C CONSTS(151) is omega in component parameters (per_ms). C CONSTS(152) is p1 in component parameters (dimensionless). C CONSTS(153) is p2 in component parameters (dimensionless). C CONSTS(154) is p3 in component parameters (dimensionless). C CONSTS(155) is pa in component parameters (per_ms). C CONSTS(156) is pb in component parameters (per_ms). C CONSTS(157) is pc1 in component parameters (per_ms). C CONSTS(158) is pc2 in component parameters (per_ms). C CONSTS(159) is PCa in component parameters (cm_per_s). C CONSTS(160) is Pi in component parameters (mM). C CONSTS(161) is PK in component parameters (cm_per_s). C CONSTS(162) is PnsK in component parameters (cm_per_s). C CONSTS(163) is PnsNa in component parameters (cm_per_s). C CONSTS(164) is r1 in component parameters (dimensionless). C CONSTS(165) is r2 in component parameters (dimensionless). C CONSTS(166) is r3 in component parameters (dimensionless). C CONSTS(167) is ra in component parameters (per_ms). C CONSTS(168) is rb in component parameters (per_ms). C CONSTS(169) is rc1 in component parameters (per_ms). C CONSTS(170) is rc2 in component parameters (per_ms). C CONSTS(171) is rhoF1 in component parameters (mM). C CONSTS(172) is rhoREF in component parameters (mM). C CONSTS(173) is rhoREN in component parameters (mM). C CONSTS(174) is tautr in component parameters (ms). C CONSTS(175) is tauxfer in component parameters (ms). C CONSTS(176) is V_AM_scaler in component parameters (dimensionless). C CONSTS(177) is V_AM_max in component parameters (mM_per_ms). C CONSTS(178) is v1 in component parameters (per_ms). C CONSTS(179) is VJSR in component parameters (uL). C CONSTS(180) is vmaxf in component parameters (per_ms). C CONSTS(181) is vmaxr in component parameters (per_ms). C CONSTS(182) is VmDT in component parameters (mM_per_ms). C CONSTS(183) is VmNC in component parameters (mM_per_ms). C CONSTS(184) is Vmuni in component parameters (mM_per_ms). C CONSTS(185) is Vmyo in component parameters (uL). C CONSTS(186) is VNSR in component parameters (uL). C CONSTS(187) is VSS in component parameters (uL). C CONSTS(188) is zeta in component parameters (N_per_mm2). C CONSTS(189) is f_xb in component parameters (per_ms). C CONSTS(190) is SL in component parameters (um). C CONSTS(191) is gmin_xb in component parameters (per_ms). C CONSTS(192) is zNa in component parameters (dimensionless). C CONSTS(193) is zK in component parameters (dimensionless). C CONSTS(194) is zCa in component parameters (dimensionless). C CONSTS(195) is one in component parameters (dimensionless). C CONSTS(196) is p_KNa in component parameters (dimensionless). C CONSTS(197) is ATPi_initial in component initial_values (dimensionless). C CONSTS(198) is V_initial in component initial_values (dimensionless). C CONSTS(199) is C1_RyR_initial in component initial_values (dimensionless). C CONSTS(200) is C2_RyR_initial in component initial_values (dimensionless). C CONSTS(201) is O2_RyR_initial in component initial_values (dimensionless). C CONSTS(202) is mNa_initial in component initial_values (dimensionless). C CONSTS(203) is hNa_initial in component initial_values (dimensionless). C CONSTS(204) is jNa_initial in component initial_values (dimensionless). C CONSTS(205) is xKs_initial in component initial_values (dimensionless). C CONSTS(206) is Open_initial in component initial_values (dimensionless). C CONSTS(207) is OCa_initial in component initial_values (dimensionless). C CONSTS(208) is yCa_initial in component initial_values (dimensionless). C CONSTS(209) is Ki_initial in component initial_values (dimensionless). C CONSTS(210) is Nai_initial in component initial_values (dimensionless). C CONSTS(211) is Cai_initial in component initial_values (dimensionless). C CONSTS(212) is ADPm_initial in component initial_values (dimensionless). C CONSTS(213) is Dpsi_initial in component initial_values (dimensionless). C CONSTS(214) is NADH_initial in component initial_values (dimensionless). C CONSTS(215) is Oaa_initial in component initial_values (dimensionless). C CONSTS(216) is AKG_initial in component initial_values (dimensionless). C CONSTS(217) is SCoA_initial in component initial_values (dimensionless). C CONSTS(218) is Succ_initial in component initial_values (dimensionless). C CONSTS(219) is FUM_initial in component initial_values (dimensionless). C CONSTS(220) is MAL_initial in component initial_values (dimensionless). C CONSTS(221) is ISOC_initial in component initial_values (dimensionless). C CONSTS(222) is Cam_initial in component initial_values (dimensionless). C CONSTS(223) is CaNSR_initial in component initial_values (dimensionless). C CONSTS(224) is CaJSR_initial in component initial_values (dimensionless). C CONSTS(225) is CaSS_initial in component initial_values (dimensionless). C CONSTS(226) is N1_initial in component initial_values (dimensionless). C CONSTS(227) is P0_initial in component initial_values (dimensionless). C CONSTS(228) is P1_initial in component initial_values (dimensionless). C CONSTS(229) is P2_initial in component initial_values (dimensionless). C CONSTS(230) is P3_initial in component initial_values (dimensionless). C CONSTS(231) is LTRPNCa_initial in component initial_values (dimensionless). C CONSTS(232) is HTRPNCa_initial in component initial_values (dimensionless). C CONSTS(233) is C1_initial in component initial_values (dimensionless). C CONSTS(234) is CCa0_initial in component initial_values (dimensionless). C CONSTS(235) is C0_initial in component initial_values (dimensionless). C CONSTS(236) is C2_initial in component initial_values (dimensionless). C CONSTS(237) is CCa1_initial in component initial_values (dimensionless). C CONSTS(238) is C3_initial in component initial_values (dimensionless). C CONSTS(239) is CCa2_initial in component initial_values (dimensionless). C CONSTS(240) is C4_initial in component initial_values (dimensionless). C CONSTS(241) is CCa3_initial in component initial_values (dimensionless). C CONSTS(242) is CCa4_initial in component initial_values (dimensionless). C CONSTS(243) is CrPi_mito_initial in component initial_values (dimensionless). C CONSTS(244) is CrPi_cyto_initial in component initial_values (dimensionless). C CONSTS(245) is ATPi_cyto_initial in component initial_values (dimensionless). C CONSTS(246) is ASP_initial in component initial_values (dimensionless). C CONSTS(247) is N0_initial in component initial_values (dimensionless). C CONSTS(248) is Vtotal in component computed_constants (uL). C CONSTS(304) is Vmito in component computed_constants (uL). C CONSTS(249) is f_01 in component computed_constants (per_ms). C CONSTS(305) is f_12 in component computed_constants (per_ms). C CONSTS(313) is f_23 in component computed_constants (per_ms). C CONSTS(317) is g0_01 in component computed_constants (per_ms). C CONSTS(322) is g0_12 in component computed_constants (per_ms). C CONSTS(327) is g0_23 in component computed_constants (per_ms). C CONSTS(330) is paths in component computed_constants (per_ms3). C CONSTS(333) is P1max in component computed_constants (dimensionless). C CONSTS(339) is P2max in component computed_constants (dimensionless). C CONSTS(344) is P3max in component computed_constants (dimensionless). C CONSTS(346) is Fmax in component computed_constants (dimensionless). C CONSTS(348) is fnormmax in component computed_constants (dimensionless). C CONSTS(250) is SLnorm in component computed_constants (dimensionless). C CONSTS(251) is Ktrop_Ca in component computed_constants (mM). C CONSTS(306) is Ktrop_half in component computed_constants (dimensionless). C CONSTS(252) is Ntrop in component computed_constants (dimensionless). C CONSTS(350) is fnormmax2 in component computed_constants (dimensionless). C CONSTS(352) is La in component computed_constants (um). C CONSTS(354) is Lm_prime in component computed_constants (um). C CONSTS(356) is Lz in component computed_constants (um). C CONSTS(358) is Lb in component computed_constants (um). C CONSTS(360) is Lm in component computed_constants (um). C CONSTS(331) is mod_factor in component computed_constants (dimensionless). C CONSTS(334) is g_01_mod in component computed_constants (per_ms). C CONSTS(335) is g_12_mod in component computed_constants (per_ms). C CONSTS(336) is g_23_mod in component computed_constants (per_ms). C CONSTS(337) is g_01_off in component computed_constants (per_ms). C CONSTS(340) is g_01_off_mod in component computed_constants (per_ms). C CONSTS(302) is RT_over_F in component computed_constants (mV). C CONSTS(253) is G_Ks in component computed_constants (mS_per_cm2). C CONSTS(311) is G_K1 in component computed_constants (mS_per_cm2). C CONSTS(254) is inv_5p98 in component computed_constants (dimensionless). C CONSTS(315) is FaradayE3 in component computed_constants (C_per_mol). C CONSTS(255) is Cao_341 in component computed_constants (mM). C CONSTS(256) is ICamax_LHospital in component computed_constants (uA_per_cm2). C CONSTS(257) is Pca_4En3 in component computed_constants (cm_per_ms). C CONSTS(320) is F_over_RT in component computed_constants (per_mV). C CONSTS(258) is inv_ICahalf in component computed_constants (cm2_per_uA). C CONSTS(318) is PKFe3 in component computed_constants (C_cm_per_s_per_mol). C CONSTS(325) is sigma in component computed_constants (dimensionless). C CONSTS(259) is inv_KmNai in component computed_constants (per_mM). C CONSTS(328) is INaKmax_Ko_Ko_KmKo in component computed_constants (uA_per_cm2). C CONSTS(332) is inv_Ki1AD_NaK in component computed_constants (per_mM). C CONSTS(260) is eta_1 in component computed_constants (dimensionless). C CONSTS(261) is Nao_p3 in component computed_constants (mM2). C CONSTS(262) is KmCa_Cao in component computed_constants (mM3_cm2_per_uA). C CONSTS(307) is KmCa_Cao_ksat in component computed_constants (mM3_cm2_per_uA). C CONSTS(263) is inv_KiADP_CaP in component computed_constants (mM). C CONSTS(264) is KmnsCa_p3 in component computed_constants (mM3). C CONSTS(319) is V_AM_scaler_max_1_f_01_12_23 in component computed_constants (mM_per_ms4). C CONSTS(265) is KmATP_AM_Ki_AM in component computed_constants (dimensionless). C CONSTS(308) is DmuH_Constant in component computed_constants (mV). C CONSTS(266) is VCS_C1 in component computed_constants (mM_per_ms). C CONSTS(267) is one_inv_KACOeq in component computed_constants (dimensionless). C CONSTS(268) is VIDH_Constant in component computed_constants (dimensionless). C CONSTS(269) is kIDH_EtID in component computed_constants (mM_per_ms). C CONSTS(270) is inv_KADP in component computed_constants (per_mM). C CONSTS(271) is inv_KaCa in component computed_constants (per_mM). C CONSTS(272) is inv_KidhNADH in component computed_constants (per_mM). C CONSTS(273) is KmKGNAD_KmIDNAD in component computed_constants (dimensionless). C CONSTS(274) is Mg_Kmg_1 in component computed_constants (dimensionless). C CONSTS(309) is Mg_Kmg_1_Kca in component computed_constants (mM). C CONSTS(275) is kKGDH_EtKG in component computed_constants (mM_per_ms). C CONSTS(276) is CoA_KSLeq in component computed_constants (mM). C CONSTS(277) is kSDH_EtSDH in component computed_constants (mM_per_ms). C CONSTS(278) is KmSucc_KiFUM in component computed_constants (dimensionless). C CONSTS(279) is inv_KiOxaa in component computed_constants (per_mM). C CONSTS(280) is kfFH_KFHeq in component computed_constants (per_ms). C CONSTS(281) is kMDH_Fh_EtMD in component computed_constants (mM_per_ms). C CONSTS(282) is Kmal_Kioaa in component computed_constants (dimensionless). C CONSTS(283) is VAAT_Constant in component computed_constants (mM_per_ms). C CONSTS(284) is kcnsASP_KAATeq_kfAAT in component computed_constants (dimensionless). C CONSTS(285) is KfAAT_GLU in component computed_constants (per_ms_per_mM). C CONSTS(286) is KfAAT_KAATeq in component computed_constants (per_ms). C CONSTS(287) is kres_sq_KmIDNAD in component computed_constants (per_mM). C CONSTS(338) is exp_6_FRT_Dpsio in component computed_constants (dimensionless). C CONSTS(343) is FRT_6_g in component computed_constants (per_mV). C CONSTS(345) is ra_rc1_exp_6_FRT_Dpsio in component computed_constants (per_ms). C CONSTS(347) is r1_exp_6_FRT_Dpsio in component computed_constants (dimensionless). C CONSTS(349) is rhoREN_ra_rc1_exp_6_FRT_Dpsio in component computed_constants (mM_per_ms). C CONSTS(351) is rhoREN_rc2 in component computed_constants (mM_per_ms). C CONSTS(353) is rhoREN_ra in component computed_constants (mM_per_ms). C CONSTS(355) is rhoRen_6_ra in component computed_constants (mM_per_ms). C CONSTS(357) is rhoRen_6_ra_rb in component computed_constants (mM_per_ms). C CONSTS(359) is AREF in component computed_constants (mV). C CONSTS(361) is exp_AREF_FRT in component computed_constants (dimensionless). C CONSTS(362) is ra_rc2_exp_AREF_FRT in component computed_constants (per_ms). C CONSTS(363) is VFO_C1 in component computed_constants (per_ms). C CONSTS(365) is ra_exp_AREF_FRT in component computed_constants (per_ms). C CONSTS(368) is ra_rb in component computed_constants (per_ms). C CONSTS(370) is VFO_VHFe_C1 in component computed_constants (per_ms). C CONSTS(371) is r2_r3_exp_AREF_FRT in component computed_constants (dimensionless). C CONSTS(372) is exp_3_FRT_Dpsio in component computed_constants (dimensionless). C CONSTS(373) is FRT_3 in component computed_constants (per_mV). C CONSTS(288) is kf1_Pi in component computed_constants (mM). C CONSTS(374) is VATPase_C1 in component computed_constants (per_ms). C CONSTS(375) is pa_pb_3 in component computed_constants (per_ms). C CONSTS(376) is pa_300 in component computed_constants (per_ms). C CONSTS(377) is p1_exp_3_FRT_Dpsio in component computed_constants (dimensionless). C CONSTS(323) is hm_F_over_RT in component computed_constants (per_mV). C CONSTS(289) is VmDT_75 in component computed_constants (mM_per_ms). C CONSTS(290) is VmDT_20 in component computed_constants (mM_per_ms). C CONSTS(291) is inv_Kfb in component computed_constants (per_mM). C CONSTS(292) is inv_Krb in component computed_constants (per_mM). C CONSTS(293) is inv_tautr in component computed_constants (per_ms). C CONSTS(294) is inv_tauxfer in component computed_constants (per_ms). C CONSTS(378) is KmATP_SR_Ki_SR in component computed_constants (dimensionless). C CONSTS(379) is inv_Ki_prime_SR in component computed_constants (per_mM). C CONSTS(366) is alpha_SL_fnormmax2 in component computed_constants (dimensionless). C CONSTS(367) is alpha_SL_fnormmax in component computed_constants (dimensionless). C CONSTS(314) is inv_LTRPNtot_Ktrop_half in component computed_constants (per_mM). C CONSTS(310) is kTrop_pn_f_01 in component computed_constants (per_ms). C CONSTS(341) is kTrop_pn_f_12_g_01_mod in component computed_constants (per_ms). C CONSTS(342) is f_23_g_12_mod in component computed_constants (per_ms). C CONSTS(295) is CMDNtot_KmCMDN in component computed_constants (mM2). C CONSTS(296) is CSQNtot_KmCSQN in component computed_constants (mM2). C CONSTS(297) is inv_ktrans in component computed_constants (per_mM). C CONSTS(298) is inv_kact in component computed_constants (per_mM). C CONSTS(299) is Vmuni_ktrans in component computed_constants (per_ms). C CONSTS(324) is FRT2 in component computed_constants (per_mV). C CONSTS(300) is b_05 in component computed_constants (dimensionless). C CONSTS(380) is Acap_Vmyo_F in component computed_constants (mmol_per_C_per_cm). C CONSTS(381) is Acap_VSS_F in component computed_constants (mmol_per_C_per_cm). C CONSTS(382) is VJSR_VSS in component computed_constants (dimensionless). C CONSTS(383) is Vmyo_VSS in component computed_constants (dimensionless). C CONSTS(384) is Vmyo_VNSR in component computed_constants (dimensionless). C CONSTS(385) is VJSR_VNSR in component computed_constants (dimensionless). C CONSTS(386) is inv_C_m in component computed_constants (cm2_per_uF). C CONSTS(303) is inv_bL in component computed_constants (dimensionless). C CONSTS(387) is inv_Cmito in component computed_constants (mV_per_mM). C CONSTS(388) is two_b in component computed_constants (dimensionless). C CONSTS(301) is inv_keq in component computed_constants (dimensionless). C CONSTS(369) is zeta_alpha_SL_fnormmax in component computed_constants (N_per_mm2). C CONSTS(364) is alpha_SL in component computed_constants (dimensionless). C ALGBRC(8) is ADP in component model_2 (dimensionless). C STATES(1) is ATPi in component model_2 (dimensionless). C ALGBRC(16) is inv_ATPi in component model_2 (dimensionless). C ALGBRC(20) is VF_over_RT in component model_2 (dimensionless). C STATES(2) is V in component model_2 (dimensionless). C ALGBRC(23) is exp_VF_over_RT in component model_2 (dimensionless). C ALGBRC(26) is VFsq_over_RT in component model_2 (dimensionless). C ALGBRC(28) is exp2VFRT in component model_2 (dimensionless). C ALGBRC(30) is O1_RyR in component model_2 (dimensionless). C STATES(3) is C1_RyR in component model_2 (dimensionless). C STATES(4) is C2_RyR in component model_2 (dimensionless). C STATES(5) is O2_RyR in component model_2 (dimensionless). C ALGBRC(1) is V_30 in component model_2 (dimensionless). C ALGBRC(125) is V_E_K in component model_2 (dimensionless). C ALGBRC(122) is INa in component model_2 (dimensionless). C STATES(6) is mNa in component model_2 (dimensionless). C STATES(7) is hNa in component model_2 (dimensionless). C STATES(8) is jNa in component model_2 (dimensionless). C ALGBRC(134) is IKs in component model_2 (dimensionless). C STATES(9) is xKs in component model_2 (dimensionless). C ALGBRC(126) is K1Alpha in component model_2 (dimensionless). C ALGBRC(127) is K1Beta in component model_2 (dimensionless). C ALGBRC(128) is K1_inf in component model_2 (dimensionless). C ALGBRC(129) is IK1 in component model_2 (dimensionless). C ALGBRC(123) is INab in component model_2 (dimensionless). C ALGBRC(130) is IKp in component model_2 (dimensionless). C ALGBRC(32) is ICamax in component model_2 (dimensionless). C ALGBRC(34) is ICaK in component model_2 (dimensionless). C STATES(10) is Open in component model_2 (dimensionless). C STATES(11) is OCa in component model_2 (dimensionless). C STATES(12) is yCa in component model_2 (dimensionless). C STATES(13) is Ki in component model_2 (dimensionless). C ALGBRC(36) is ICa in component model_2 (dimensionless). C ALGBRC(38) is NaiP1p5 in component model_2 (dimensionless). C STATES(14) is Nai in component model_2 (dimensionless). C ALGBRC(40) is INaK in component model_2 (dimensionless). C ALGBRC(42) is exp_eta_VF_over_RT in component model_2 (dimensionless). C ALGBRC(44) is exp_eta1_VF_over_RT in component model_2 (dimensionless). C ALGBRC(46) is INaCa in component model_2 (dimensionless). C STATES(15) is Cai in component model_2 (dimensionless). C ALGBRC(132) is ICab in component model_2 (dimensionless). C ALGBRC(48) is IpCa in component model_2 (dimensionless). C ALGBRC(50) is CaiP3 in component model_2 (dimensionless). C ALGBRC(52) is common in component model_2 (dimensionless). C ALGBRC(54) is InsNa in component model_2 (dimensionless). C ALGBRC(56) is InsK in component model_2 (dimensionless). C ALGBRC(58) is InsCa in component model_2 (dimensionless). C ALGBRC(60) is V_AM in component model_2 (dimensionless). C ALGBRC(62) is ATPm in component model_2 (dimensionless). C STATES(16) is ADPm in component model_2 (dimensionless). C ALGBRC(64) is DmuH in component model_2 (dimensionless). C STATES(17) is Dpsi in component model_2 (dimensionless). C ALGBRC(66) is NAD in component model_2 (dimensionless). C STATES(18) is NADH in component model_2 (dimensionless). C ALGBRC(68) is KmIDNAD_NAD in component model_2 (dimensionless). C ALGBRC(69) is exp_FRT_6_g_DmuH in component model_2 (dimensionless). C ALGBRC(70) is FRT2_Dpsi in component model_2 (dimensionless). C ALGBRC(71) is VCS in component model_2 (dimensionless). C STATES(19) is Oaa in component model_2 (dimensionless). C ALGBRC(72) is VACO in component model_2 (dimensionless). C STATES(20) is AKG in component model_2 (dimensionless). C STATES(21) is SCoA in component model_2 (dimensionless). C STATES(22) is Succ in component model_2 (dimensionless). C STATES(23) is FUM in component model_2 (dimensionless). C STATES(24) is MAL in component model_2 (dimensionless). C STATES(25) is ISOC in component model_2 (dimensionless). C ALGBRC(73) is Fa in component model_2 (dimensionless). C STATES(26) is Cam in component model_2 (dimensionless). C ALGBRC(74) is Fi in component model_2 (dimensionless). C ALGBRC(75) is VIDH in component model_2 (dimensionless). C ALGBRC(76) is a__1 in component model_2 (dimensionless). C ALGBRC(77) is VKGDH in component model_2 (dimensionless). C ALGBRC(78) is VSL in component model_2 (dimensionless). C ALGBRC(79) is VSDH in component model_2 (dimensionless). C ALGBRC(80) is VFH in component model_2 (dimensionless). C ALGBRC(81) is VMDH in component model_2 (dimensionless). C ALGBRC(82) is VAAT in component model_2 (dimensionless). C ALGBRC(83) is AREN in component model_2 (dimensionless). C ALGBRC(84) is denominator1 in component model_2 (dimensionless). C ALGBRC(87) is denominator2 in component model_2 (dimensionless). C ALGBRC(91) is denominator3 in component model_2 (dimensionless). C ALGBRC(85) is VNO in component model_2 (dimensionless). C ALGBRC(86) is VHNe in component model_2 (dimensionless). C ALGBRC(88) is VHFe in component model_2 (dimensionless). C ALGBRC(89) is exp_3FRT_DmuH in component model_2 (dimensionless). C ALGBRC(90) is AF1 in component model_2 (dimensionless). C ALGBRC(92) is VATPase in component model_2 (dimensionless). C ALGBRC(93) is Vhu in component model_2 (dimensionless). C ALGBRC(94) is ATPi_ADP in component model_2 (dimensionless). C ALGBRC(95) is ADPm_ATPm in component model_2 (dimensionless). C ALGBRC(96) is VANT in component model_2 (dimensionless). C ALGBRC(97) is Vhleak in component model_2 (dimensionless). C ALGBRC(3) is MAlpha in component model_2 (dimensionless). C ALGBRC(9) is MBeta in component model_2 (dimensionless). C ALGBRC(14) is inv_MBeta_MAlpha in component model_2 (dimensionless). C ALGBRC(18) is tmNa in component model_2 (dimensionless). C ALGBRC(4) is HAlpha in component model_2 (dimensionless). C ALGBRC(10) is HBeta in component model_2 (dimensionless). C ALGBRC(5) is JAlpha in component model_2 (dimensionless). C ALGBRC(11) is JBeta in component model_2 (dimensionless). C ALGBRC(98) is fb in component model_2 (dimensionless). C ALGBRC(99) is rb in component model_2 (per_ms). C STATES(27) is CaNSR in component model_2 (dimensionless). C ALGBRC(100) is Jup in component model_2 (dimensionless). C ALGBRC(102) is Jrel in component model_2 (dimensionless). C STATES(28) is CaJSR in component model_2 (dimensionless). C STATES(29) is CaSS in component model_2 (dimensionless). C ALGBRC(103) is Jtr in component model_2 (dimensionless). C ALGBRC(106) is Jxfer in component model_2 (dimensionless). C ALGBRC(109) is P1_N1_P2_P3 in component model_2 (dimensionless). C STATES(30) is N1 in component model_2 (dimensionless). C STATES(31) is P0 in component model_2 (dimensionless). C STATES(32) is P1 in component model_2 (dimensionless). C STATES(33) is P2 in component model_2 (dimensionless). C STATES(34) is P3 in component model_2 (dimensionless). C ALGBRC(25) is N1_differential in component model_2 (dimensionless). C ALGBRC(12) is P0_differential in component model_2 (dimensionless). C ALGBRC(15) is P1_differential in component model_2 (dimensionless). C ALGBRC(19) is P2_differential in component model_2 (dimensionless). C ALGBRC(22) is P3_differential in component model_2 (dimensionless). C ALGBRC(112) is FN_Ca in component model_2 (dimensionless). C ALGBRC(110) is force_norm in component model_2 (dimensionless). C ALGBRC(111) is force in component model_2 (dimensionless). C ALGBRC(6) is kTrop_np in component model_2 (dimensionless). C STATES(35) is LTRPNCa in component model_2 (dimensionless). C STATES(36) is HTRPNCa in component model_2 (dimensionless). C ALGBRC(113) is LTRPNCa_differential in component model_2 (dimensionless). C ALGBRC(114) is HTRPNCa_differential in component model_2 (dimensionless). C ALGBRC(115) is Jtrpn in component model_2 (dimensionless). C ALGBRC(107) is beta_SS in component model_2 (dimensionless). C ALGBRC(108) is beta_JSR in component model_2 (dimensionless). C ALGBRC(116) is beta_i in component model_2 (dimensionless). C ALGBRC(117) is Cai_ktrans_plus1 in component model_2 (dimensionless). C ALGBRC(118) is Cai_ktrans_plus1_p3 in component model_2 (dimensionless). C ALGBRC(119) is Vuni in component model_2 (dimensionless). C ALGBRC(120) is VnaCa in component model_2 (dimensionless). C ALGBRC(7) is alpha in component model_2 (dimensionless). C ALGBRC(13) is beta in component model_2 (dimensionless). C ALGBRC(17) is alpha_prime in component model_2 (dimensionless). C ALGBRC(21) is beta_prime in component model_2 (dimensionless). C ALGBRC(24) is C0_to_C1 in component model_2 (dimensionless). C ALGBRC(27) is C1_to_C2 in component model_2 (dimensionless). C ALGBRC(29) is C2_to_C3 in component model_2 (dimensionless). C ALGBRC(31) is C3_to_C4 in component model_2 (dimensionless). C ALGBRC(33) is CCa0_to_CCa1 in component model_2 (dimensionless). C ALGBRC(35) is CCa1_to_CCa2 in component model_2 (dimensionless). C ALGBRC(37) is CCa2_to_CCa3 in component model_2 (dimensionless). C ALGBRC(39) is CCa3_to_CCa4 in component model_2 (dimensionless). C ALGBRC(41) is C1_to_C0 in component model_2 (dimensionless). C ALGBRC(43) is C2_to_C1 in component model_2 (dimensionless). C ALGBRC(45) is C3_to_C2 in component model_2 (dimensionless). C ALGBRC(47) is C4_to_C3 in component model_2 (dimensionless). C ALGBRC(49) is CCa1_to_CCa0 in component model_2 (dimensionless). C ALGBRC(51) is CCa2_to_CCa1 in component model_2 (dimensionless). C ALGBRC(53) is CCa3_to_CCa2 in component model_2 (dimensionless). C ALGBRC(55) is CCa4_to_CCa3 in component model_2 (dimensionless). C ALGBRC(57) is gamma in component model_2 (dimensionless). C ALGBRC(59) is C0_to_CCa0 in component model_2 (dimensionless). C ALGBRC(61) is C1_to_CCa1 in component model_2 (dimensionless). C ALGBRC(63) is C2_to_CCa2 in component model_2 (dimensionless). C ALGBRC(65) is C3_to_CCa3 in component model_2 (dimensionless). C ALGBRC(67) is C4_to_CCa4 in component model_2 (dimensionless). C CONSTS(312) is CCa0_to_C0 in component model_2 (dimensionless). C CONSTS(316) is CCa1_to_C1 in component model_2 (dimensionless). C CONSTS(321) is CCa2_to_C2 in component model_2 (dimensionless). C CONSTS(326) is CCa3_to_C3 in component model_2 (dimensionless). C CONSTS(329) is CCa4_to_C4 in component model_2 (dimensionless). C STATES(37) is C1 in component model_2 (dimensionless). C STATES(38) is CCa0 in component model_2 (dimensionless). C STATES(39) is C0 in component model_2 (dimensionless). C STATES(40) is C2 in component model_2 (dimensionless). C STATES(41) is CCa1 in component model_2 (dimensionless). C STATES(42) is C3 in component model_2 (dimensionless). C STATES(43) is CCa2 in component model_2 (dimensionless). C STATES(44) is C4 in component model_2 (dimensionless). C STATES(45) is CCa3 in component model_2 (dimensionless). C STATES(46) is CCa4 in component model_2 (dimensionless). C ALGBRC(101) is Vt_CRP2 in component model_2 (dimensionless). C STATES(47) is CrPi_mito in component model_2 (dimensionless). C STATES(48) is CrPi_cyto in component model_2 (dimensionless). C ALGBRC(104) is VCK_cyto in component model_2 (dimensionless). C STATES(49) is ATPi_cyto in component model_2 (dimensionless). C ALGBRC(105) is VCK_mito in component model_2 (dimensionless). C STATES(50) is ASP in component model_2 (dimensionless). C STATES(51) is N0 in component model_2 (dimensionless). C ALGBRC(124) is reversal_potential in component EK (dimensionless). C ALGBRC(121) is reversal_potential in component ENa (dimensionless). C ALGBRC(133) is reversal_potential in component EKs (dimensionless). C ALGBRC(131) is reversal_potential in component ECa (dimensionless). C RATES(6) is d/dt mNa in component model_2 (dimensionless). C RATES(7) is d/dt hNa in component model_2 (dimensionless). C RATES(8) is d/dt jNa in component model_2 (dimensionless). C RATES(9) is d/dt xKs in component model_2 (dimensionless). C RATES(31) is d/dt P0 in component model_2 (dimensionless). C RATES(32) is d/dt P1 in component model_2 (dimensionless). C RATES(33) is d/dt P2 in component model_2 (dimensionless). C RATES(34) is d/dt P3 in component model_2 (dimensionless). C RATES(30) is d/dt N1 in component model_2 (dimensionless). C RATES(51) is d/dt N0 in component model_2 (dimensionless). C RATES(35) is d/dt LTRPNCa in component model_2 (dimensionless). C RATES(36) is d/dt HTRPNCa in component model_2 (dimensionless). C RATES(14) is d/dt Nai in component model_2 (dimensionless). C RATES(13) is d/dt Ki in component model_2 (dimensionless). C RATES(15) is d/dt Cai in component model_2 (dimensionless). C RATES(29) is d/dt CaSS in component model_2 (dimensionless). C RATES(28) is d/dt CaJSR in component model_2 (dimensionless). C RATES(27) is d/dt CaNSR in component model_2 (dimensionless). C RATES(2) is d/dt V in component model_2 (dimensionless). C RATES(3) is d/dt C1_RyR in component model_2 (dimensionless). C RATES(5) is d/dt O2_RyR in component model_2 (dimensionless). C RATES(4) is d/dt C2_RyR in component model_2 (dimensionless). C RATES(39) is d/dt C0 in component model_2 (dimensionless). C RATES(37) is d/dt C1 in component model_2 (dimensionless). C RATES(40) is d/dt C2 in component model_2 (dimensionless). C RATES(42) is d/dt C3 in component model_2 (dimensionless). C RATES(44) is d/dt C4 in component model_2 (dimensionless). C RATES(10) is d/dt Open in component model_2 (dimensionless). C RATES(38) is d/dt CCa0 in component model_2 (dimensionless). C RATES(41) is d/dt CCa1 in component model_2 (dimensionless). C RATES(43) is d/dt CCa2 in component model_2 (dimensionless). C RATES(45) is d/dt CCa3 in component model_2 (dimensionless). C RATES(46) is d/dt CCa4 in component model_2 (dimensionless). C RATES(12) is d/dt yCa in component model_2 (dimensionless). C RATES(11) is d/dt OCa in component model_2 (dimensionless). C RATES(47) is d/dt CrPi_mito in component model_2 (dimensionless). C RATES(48) is d/dt CrPi_cyto in component model_2 (dimensionless). C RATES(1) is d/dt ATPi in component model_2 (dimensionless). C RATES(49) is d/dt ATPi_cyto in component model_2 (dimensionless). C RATES(26) is d/dt Cam in component model_2 (dimensionless). C RATES(16) is d/dt ADPm in component model_2 (dimensionless). C RATES(17) is d/dt Dpsi in component model_2 (dimensionless). C RATES(18) is d/dt NADH in component model_2 (dimensionless). C RATES(25) is d/dt ISOC in component model_2 (dimensionless). C RATES(20) is d/dt AKG in component model_2 (dimensionless). C RATES(21) is d/dt SCoA in component model_2 (dimensionless). C RATES(22) is d/dt Succ in component model_2 (dimensionless). C RATES(23) is d/dt FUM in component model_2 (dimensionless). C RATES(24) is d/dt MAL in component model_2 (dimensionless). C RATES(19) is d/dt Oaa in component model_2 (dimensionless). C RATES(50) is d/dt ASP in component model_2 (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 4e3 CONSTS(2) = 0.5 CONSTS(3) = -100.0 CONSTS(4) = 0.0 CONSTS(5) = 0.0 CONSTS(6) = 96.5 CONSTS(7) = 8.314 CONSTS(8) = 310.0 CONSTS(9) = 1e-7 CONSTS(10) = 2.0e-3 CONSTS(11) = 1.4e-4 CONSTS(12) = 1.33e-6 CONSTS(13) = 0.0095 CONSTS(14) = 25.0 CONSTS(15) = 25.0 CONSTS(16) = 1.0E-5 CONSTS(17) = 1.534E-4 CONSTS(18) = 1.0 CONSTS(19) = 2.0 CONSTS(20) = 0.5 CONSTS(21) = 2.0 CONSTS(22) = 1.0 CONSTS(23) = 2.0 CONSTS(24) = 1.0 CONSTS(25) = 1.5 CONSTS(26) = 50E-3 CONSTS(27) = 1.812E-3 CONSTS(28) = 10.0 CONSTS(29) = 0.02 CONSTS(30) = 35.0 CONSTS(31) = -0.6 CONSTS(32) = 50.0 CONSTS(33) = 0.35 CONSTS(34) = 0.4 CONSTS(35) = 0.109 CONSTS(36) = 0.5 CONSTS(37) = 0.154 CONSTS(38) = 0.5 CONSTS(39) = 0.01 CONSTS(40) = 1.24 CONSTS(41) = 0.3 CONSTS(42) = 0.0003 CONSTS(43) = 0.0 CONSTS(44) = 0.85 CONSTS(45) = 0.003217 CONSTS(46) = 0.00828 CONSTS(47) = 12.8 CONSTS(48) = 0.000545 CONSTS(49) = 2.0E-7 CONSTS(50) = 2.0 CONSTS(51) = 10.0 CONSTS(52) = 0.0 CONSTS(53) = 2.5E-5 CONSTS(54) = 0.5 CONSTS(55) = 140E-3 CONSTS(56) = -0.4583 CONSTS(57) = 3.147 CONSTS(58) = 0.575 CONSTS(59) = 6.6 CONSTS(60) = 0.0005 CONSTS(61) = 2.22 CONSTS(62) = 0.00038 CONSTS(63) = 0.62 CONSTS(64) = 0.576 CONSTS(65) = 1.215E10 CONSTS(66) = 1.930 CONSTS(67) = 4.05E6 CONSTS(68) = 0.00127 CONSTS(69) = 0.0008 CONSTS(70) = 1.5E-6 CONSTS(71) = 0.10 CONSTS(72) = 0.5 CONSTS(73) = 1.71E6 CONSTS(74) = 6.44E-4 CONSTS(75) = 0.0125 CONSTS(76) = 0.00024 CONSTS(77) = 0.00332 CONSTS(78) = 1.0 CONSTS(79) = 0.005 CONSTS(80) = 8.1E-5 CONSTS(81) = 5.98E-5 CONSTS(82) = 1.131E-5 CONSTS(83) = 26.7 CONSTS(84) = 6.68E-9 CONSTS(85) = 5.62E-6 CONSTS(86) = 3.3E-4 CONSTS(87) = 100.0 CONSTS(88) = 0.26 CONSTS(89) = 5.1 CONSTS(90) = 0.14 CONSTS(91) = 0.1 CONSTS(92) = 1.0 CONSTS(93) = 0.05 CONSTS(94) = 0.19 CONSTS(95) = 1.3 CONSTS(96) = 0.0031 CONSTS(97) = 0.15 CONSTS(98) = 7.5E-2 CONSTS(99) = 100.0 CONSTS(100) = 40E-3 CONSTS(101) = 0.008 CONSTS(102) = 0.012 CONSTS(103) = 0.23 CONSTS(104) = 0.0126 CONSTS(105) = 1.493 CONSTS(106) = 0.03 CONSTS(107) = 0.01 CONSTS(108) = 2.38E-3 CONSTS(109) = 1.38 CONSTS(110) = 0.8 CONSTS(111) = 0.111 CONSTS(112) = 0.0308 CONSTS(113) = 0.923 CONSTS(114) = 1.52 CONSTS(115) = 1.94 CONSTS(116) = 38.7 CONSTS(117) = 1.5 CONSTS(118) = 0.2244 CONSTS(119) = 87.5 CONSTS(120) = 10.0 CONSTS(121) = 0.0012 CONSTS(122) = 0.00064 CONSTS(123) = 0.0005 CONSTS(124) = 0.03 CONSTS(125) = 9.4 CONSTS(126) = 9000.0 CONSTS(127) = 3.75E-4 CONSTS(128) = 5.4 CONSTS(129) = 0.0399 CONSTS(130) = 1.64269 CONSTS(131) = 1.35E18 CONSTS(132) = 5.765E13 CONSTS(133) = 0.1 CONSTS(134) = 0.005 CONSTS(135) = 3.115 CONSTS(136) = 2.5 CONSTS(137) = 0.019 CONSTS(138) = 0.04 CONSTS(139) = 110.0 CONSTS(140) = 70E-3 CONSTS(141) = 3.0 CONSTS(142) = 0.4 CONSTS(143) = 3.0 CONSTS(144) = 2.8 CONSTS(145) = 140.0 CONSTS(146) = 4.0 CONSTS(147) = 1.4 CONSTS(148) = 2.0 CONSTS(149) = 1.2 CONSTS(150) = 1.0 CONSTS(151) = 0.01 CONSTS(152) = 1.346E-8 CONSTS(153) = 7.739E-7 CONSTS(154) = 6.65E-15 CONSTS(155) = 1.656E-8 CONSTS(156) = 3.373E-10 CONSTS(157) = 9.651E-17 CONSTS(158) = 4.585E-17 CONSTS(159) = 0.00124 CONSTS(160) = 2.0 CONSTS(161) = 1.11E-11 CONSTS(162) = 0.0 CONSTS(163) = 1.75E-7 CONSTS(164) = 2.077E-18 CONSTS(165) = 1.728E-9 CONSTS(166) = 1.059E-26 CONSTS(167) = 6.394E-13 CONSTS(168) = 1.762E-16 CONSTS(169) = 2.656E-22 CONSTS(170) = 8.632E-30 CONSTS(171) = 0.05 CONSTS(172) = 3.75E-4 CONSTS(173) = 1.0E-1 CONSTS(174) = 0.574713 CONSTS(175) = 9.09091 CONSTS(176) = 15.0 CONSTS(177) = 0.00048 CONSTS(178) = 3.6 CONSTS(179) = 0.16E-6 CONSTS(180) = 0.29892E-3 CONSTS(181) = 0.31793E-3 CONSTS(182) = 0.015 CONSTS(183) = 0.625E-4 CONSTS(184) = 0.0275 CONSTS(185) = 25.84E-6 CONSTS(186) = 1.4E-6 CONSTS(187) = 0.495E-9 CONSTS(188) = 0.1 CONSTS(189) = 0.05 CONSTS(190) = 2.15 CONSTS(191) = 0.1 CONSTS(192) = 1 CONSTS(193) = 1 CONSTS(194) = 2 CONSTS(195) = 1 CONSTS(196) = 0.01833 CONSTS(197) = 7.98691 CONSTS(198) = -86.4154 CONSTS(199) = 0.748798 CONSTS(200) = 0.250904 CONSTS(201) = 3.5175e-10 CONSTS(202) = 0.0327969 CONSTS(203) = 0.988991 CONSTS(204) = 0.992918 CONSTS(205) = 0.000127552 CONSTS(206) = 6.96577e-24 CONSTS(207) = 0 CONSTS(208) = 0.941505 CONSTS(209) = 148.065 CONSTS(210) = 7.51215 CONSTS(211) = 6.71512e-05 CONSTS(212) = 0.0258454 CONSTS(213) = 164.073 CONSTS(214) = 7.76311 CONSTS(215) = 1.08962e-07 CONSTS(216) = 0.000610554 CONSTS(217) = 0.247364 CONSTS(218) = 0.000394005 CONSTS(219) = 0.0191924 CONSTS(220) = 0.00925885 CONSTS(221) = 0.496682 CONSTS(222) = 0.000137836 CONSTS(223) = 0.255288 CONSTS(224) = 0.255138 CONSTS(225) = 8.25461e-05 CONSTS(226) = 5.71121e-05 CONSTS(227) = 6.63724e-05 CONSTS(228) = 5.73114e-05 CONSTS(229) = 0.000106992 CONSTS(230) = 9.31884e-05 CONSTS(231) = 0.0100663 CONSTS(232) = 0.133453 CONSTS(233) = 1.04295e-05 CONSTS(234) = 0.00156783 CONSTS(235) = 0.998424 CONSTS(236) = 4.08551e-11 CONSTS(237) = 6.55103e-08 CONSTS(238) = 7.11287e-17 CONSTS(239) = 1.02648e-12 CONSTS(240) = 4.64382e-23 CONSTS(241) = 7.14846e-18 CONSTS(242) = 1.86682e-23 CONSTS(243) = 19.4385 CONSTS(244) = 19.4302 CONSTS(245) = 7.97824 CONSTS(246) = 0.0502614 CONSTS(247) = 0.999618 STATES(1) = 7.98691 STATES(2) = -86.4154 STATES(3) = 0.748798 STATES(4) = 0.250904 STATES(5) = 3.5175e-10 STATES(6) = 0.0327969 STATES(7) = 0.988991 STATES(8) = 0.992918 STATES(9) = 0.000127552 STATES(10) = 6.96577e-24 STATES(11) = 0 STATES(12) = 0.941505 STATES(13) = 148.065 STATES(14) = 7.51215 STATES(15) = 6.71512e-05 STATES(16) = 0.0258454 STATES(17) = 164.073 STATES(18) = 7.76311 STATES(19) = 1.08962e-07 STATES(20) = 0.000610554 STATES(21) = 0.247364 STATES(22) = 0.000394005 STATES(23) = 0.0191924 STATES(24) = 0.00925885 STATES(25) = 0.496682 STATES(26) = 0.000137836 STATES(27) = 0.255288 STATES(28) = 0.255138 STATES(29) = 8.25461e-05 STATES(30) = 5.71121e-05 STATES(31) = 6.63724e-05 STATES(32) = 5.73114e-05 STATES(33) = 0.000106992 STATES(34) = 9.31884e-05 STATES(35) = 0.0100663 STATES(36) = 0.133453 STATES(37) = 1.04295e-05 STATES(38) = 0.00156783 STATES(39) = 0.998424 STATES(40) = 4.08551e-11 STATES(41) = 6.55103e-08 STATES(42) = 7.11287e-17 STATES(43) = 1.02648e-12 STATES(44) = 4.64382e-23 STATES(45) = 7.14846e-18 STATES(46) = 1.86682e-23 STATES(47) = 19.4385 STATES(48) = 19.4302 STATES(49) = 7.97824 STATES(50) = 0.0502614 STATES(51) = 0.999618 CONSTS(248) = (CONSTS(185)+CONSTS(179)+CONSTS(186)+CONSTS(187))/0.640000 CONSTS(249) = 3.00000*CONSTS(189) CONSTS(250) = (CONSTS(190) - 1.70000)/0.600000 CONSTS(251) = CONSTS(100)/CONSTS(99) CONSTS(252) = 3.50000*CONSTS(190) - 2.00000 CONSTS(253) = 0.282000* (CONSTS(128)/5.40000) ** (1.0 / 2) CONSTS(254) = 1.00000/5.98000 CONSTS(255) = CONSTS(23)*341.000 CONSTS(256) = 2.00000*CONSTS(159)*1000.00*CONSTS(6)*(1.00000 - 341.000*CONSTS(23)) CONSTS(257) = 4.00000*0.00100000*CONSTS(159) CONSTS(258) = 1.00000/CONSTS(56) CONSTS(259) = 1.00000/CONSTS(120) CONSTS(260) = CONSTS(33) - 1.00000 CONSTS(261) = CONSTS(145) ** 3.00000/CONSTS(23) CONSTS(262) = (CONSTS(109)+CONSTS(23))*((CONSTS(119) ** 3.00000+CONSTS(145) ** 3.00000)/( CONSTS(126)*CONSTS(23))) CONSTS(263) = 1.00000/CONSTS(92) CONSTS(264) = CONSTS(121) ** 3.00000 CONSTS(265) = CONSTS(106)/CONSTS(88) CONSTS(266) = ( CONSTS(72)*CONSTS(34)*CONSTS(18))/(CONSTS(104)+CONSTS(18)) CONSTS(267) = 1.00000+1.00000/CONSTS(61) CONSTS(268) = 1.00000+CONSTS(53)/CONSTS(80)+CONSTS(81)/CONSTS(53) CONSTS(269) = CONSTS(93)*CONSTS(35) CONSTS(270) = 1.00000/CONSTS(63) CONSTS(271) = 1.00000/CONSTS(60) CONSTS(272) = 1.00000/CONSTS(94) CONSTS(273) = CONSTS(116)/CONSTS(113) CONSTS(274) = CONSTS(142)/CONSTS(112)+1.00000 CONSTS(275) = CONSTS(98)*CONSTS(36) CONSTS(276) = CONSTS(29)/CONSTS(135) CONSTS(277) = CONSTS(134)*CONSTS(38) CONSTS(278) = CONSTS(124)/CONSTS(95) CONSTS(279) = 1.00000/CONSTS(97) CONSTS(280) = CONSTS(77)/CONSTS(78) CONSTS(281) = 1.00000/(1.00000+CONSTS(84)/CONSTS(53)+ CONSTS(84)*(CONSTS(85)/CONSTS(53) ** 2.00000)) ** 2.00000*(1.00000/(1.00000+CONSTS(53)/CONSTS(82)+CONSTS(53) ** 2.00000/( CONSTS(82)*CONSTS(83)))+CONSTS(129))*CONSTS(111)*CONSTS(37) CONSTS(282) = CONSTS(105)/CONSTS(96) CONSTS(283) = CONSTS(74)*CONSTS(51)*CONSTS(70)*(CONSTS(59)/CONSTS(74)) CONSTS(284) = CONSTS(70)*(CONSTS(59)/CONSTS(74)) CONSTS(285) = CONSTS(74)*CONSTS(51) CONSTS(286) = CONSTS(74)/CONSTS(59) CONSTS(287) = CONSTS(131)*(CONSTS(131)/CONSTS(113)) CONSTS(288) = CONSTS(73)/CONSTS(160) CONSTS(289) = 0.750000*CONSTS(182) CONSTS(290) = 20.0000*CONSTS(182) CONSTS(291) = 1.00000/CONSTS(76) CONSTS(292) = 1.00000/CONSTS(130) CONSTS(293) = 1.00000/CONSTS(174) CONSTS(294) = 1.00000/CONSTS(175) CONSTS(295) = CONSTS(26)*CONSTS(108) CONSTS(296) = CONSTS(30)*CONSTS(110) CONSTS(297) = 1.00000/CONSTS(137) CONSTS(298) = 1.00000/CONSTS(62) CONSTS(299) = CONSTS(184)/CONSTS(137) CONSTS(300) = CONSTS(20)*0.500000 CONSTS(301) = 1.00000/CONSTS(13) CONSTS(302) = ( CONSTS(7)*CONSTS(8))/CONSTS(6) CONSTS(303) = 1.00000/CONSTS(21) CONSTS(304) = CONSTS(248)*0.360000 CONSTS(305) = 10.0000*CONSTS(189) CONSTS(306) = 1.00000/(1.00000+CONSTS(251)/(1.70000/1000.00+ ((0.900000/1000.00 - 1.70000/1000.00)/(2.30000 - 1.70000))*(CONSTS(190) - 1.70000))) CONSTS(307) = CONSTS(262)*CONSTS(133) CONSTS(308) = -2.30300*CONSTS(302)*CONSTS(31) CONSTS(309) = CONSTS(274)/CONSTS(68) CONSTS(310) = - CONSTS(138) - CONSTS(249) CONSTS(311) = 0.750000* (CONSTS(128)/5.40000) ** (1.0 / 2) CONSTS(312) = CONSTS(151) CONSTS(313) = 7.00000*CONSTS(189) CONSTS(314) = 1.00000/( CONSTS(140)*CONSTS(306)) CONSTS(315) = 1000.00*CONSTS(6) CONSTS(316) = CONSTS(312)*CONSTS(303) CONSTS(317) = 1.00000*CONSTS(191) CONSTS(318) = CONSTS(315)*CONSTS(161) CONSTS(319) = CONSTS(176)*(CONSTS(177)/(CONSTS(249)+CONSTS(305)+CONSTS(313))) CONSTS(320) = 1.00000/CONSTS(302) CONSTS(321) = CONSTS(316)*CONSTS(303) CONSTS(322) = 2.00000*CONSTS(191) CONSTS(323) = CONSTS(54)*CONSTS(320) CONSTS(324) = 2.00000*CONSTS(320) CONSTS(325) = 0.0365000*((EXP(CONSTS(145)/67.3000) - 1.00000)/7.00000) CONSTS(326) = CONSTS(321)*CONSTS(303) CONSTS(327) = 3.00000*CONSTS(191) CONSTS(328) = CONSTS(57)*(CONSTS(128)/(CONSTS(128)+CONSTS(117))) CONSTS(329) = CONSTS(326)*CONSTS(303) CONSTS(330) = CONSTS(317)*CONSTS(322)*CONSTS(327)+ CONSTS(249)*CONSTS(322)*CONSTS(327)+ CONSTS(249)*CONSTS(305)*CONSTS(327)+ CONSTS(249)*CONSTS(305)*CONSTS(313) CONSTS(331) = 1.00000+(2.30000 - CONSTS(190))/(2.30000 - 1.70000)/1.00000 ** 1.60000 CONSTS(332) = 1.00000/CONSTS(91) CONSTS(333) = ( CONSTS(249)*( 2.00000*CONSTS(191))*( 3.00000*CONSTS(191)))/CONSTS(330) CONSTS(334) = CONSTS(317)*CONSTS(331) CONSTS(335) = CONSTS(322)*CONSTS(331) CONSTS(336) = CONSTS(327)*CONSTS(331) CONSTS(337) = 30.0000/1000.00 CONSTS(338) = EXP( 6.00000*CONSTS(32)*CONSTS(320)) CONSTS(339) = ( CONSTS(249)*CONSTS(305)*( 3.00000*CONSTS(191)))/CONSTS(330) CONSTS(340) = CONSTS(337)*CONSTS(331) CONSTS(341) = - (CONSTS(138)+CONSTS(305)+CONSTS(334)) CONSTS(342) = - (CONSTS(313)+CONSTS(335)) CONSTS(343) = 6.00000*CONSTS(44)*CONSTS(320) CONSTS(344) = ( CONSTS(249)*CONSTS(305)*CONSTS(313))/CONSTS(330) CONSTS(345) = CONSTS(167)+ CONSTS(169)*CONSTS(338) CONSTS(346) = CONSTS(333)+ 2.00000*CONSTS(339)+ 3.00000*CONSTS(344) CONSTS(347) = CONSTS(164)*CONSTS(338) CONSTS(348) = CONSTS(346)/3.00000 CONSTS(349) = 0.500000*CONSTS(173)*CONSTS(345) CONSTS(350) = CONSTS(333)+CONSTS(339)+CONSTS(344) CONSTS(351) = 0.500000*CONSTS(173)*CONSTS(170) CONSTS(352) = 1.00000 CONSTS(353) = 0.500000*CONSTS(173)*CONSTS(167) CONSTS(354) = 1.50000 CONSTS(355) = 6.00000*CONSTS(173)*CONSTS(167) CONSTS(356) = 0.100000 CONSTS(357) = 6.00000*CONSTS(173)*(CONSTS(167)+CONSTS(168)) CONSTS(358) = 0.100000 CONSTS(359) = CONSTS(302)*arbitrary_log( CONSTS(132)* (CONSTS(40)/CONSTS(39)) ** (1.0 / 2), 10) CONSTS(360) = CONSTS(354) - CONSTS(358) CONSTS(361) = EXP( CONSTS(359)*CONSTS(320)) CONSTS(362) = 0.500000*(CONSTS(167)+ CONSTS(170)*CONSTS(361)) CONSTS(363) = (CONSTS(167)+ CONSTS(169)*CONSTS(338))*CONSTS(361)*0.500000 CONSTS(364) = TERNRY(CONSTS(190).LT.2.20000, MIN(2, 1.00000, ((CONSTS(190) - 2.00000*CONSTS(352))+(CONSTS(354) - CONSTS(356)))/CONSTS(360)), 1.00000 - (CONSTS(190) - 2.20000)/CONSTS(360)) CONSTS(365) = 4.00000*CONSTS(167)*CONSTS(361) CONSTS(366) = CONSTS(364)/CONSTS(350) CONSTS(367) = CONSTS(364)/( CONSTS(348)*3.00000) CONSTS(368) = 4.00000*(CONSTS(167)+CONSTS(168)) CONSTS(369) = CONSTS(188)*CONSTS(367) CONSTS(370) = (1.00000+ CONSTS(164)*CONSTS(361))*CONSTS(338) CONSTS(371) = CONSTS(165)+ CONSTS(166)*CONSTS(361) CONSTS(372) = EXP( 3.00000*CONSTS(32)*CONSTS(320)) CONSTS(373) = 3.00000*CONSTS(320) CONSTS(374) = 100.000*CONSTS(155)+ CONSTS(157)*CONSTS(372) CONSTS(375) = 3.00000*(CONSTS(155)+CONSTS(156)) CONSTS(376) = 300.000*CONSTS(155) CONSTS(377) = CONSTS(152)*CONSTS(372) CONSTS(378) = CONSTS(107)/CONSTS(90) CONSTS(379) = 1.00000/CONSTS(89) CONSTS(380) = CONSTS(17)/( CONSTS(185)*CONSTS(6)*1000.00) CONSTS(381) = CONSTS(17)/( 2.00000*CONSTS(187)*CONSTS(6)*1000.00) CONSTS(382) = CONSTS(179)/CONSTS(187) CONSTS(383) = CONSTS(185)/CONSTS(187) CONSTS(384) = CONSTS(185)/CONSTS(186) CONSTS(385) = CONSTS(179)/CONSTS(186) CONSTS(386) = 1.00000/CONSTS(22) CONSTS(387) = 1.00000/CONSTS(27) CONSTS(388) = 2.00000*CONSTS(20) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(10) = CONSTS(41)*STATES(44) - CONSTS(50)*STATES(10) RATES(12) = (1.00000/(1.00000+EXP( (STATES(2)+55.0000)*(1.00000/7.50000)))+(0.500000/(1.00000+EXP( (21.0000 - STATES(2))*(1.00000/6.00000))) - STATES(12)))/(20.0000+600.000/(1.00000+EXP( (STATES(2)+30.0000)*(1.00000/9.50000)))) RATES(11) = CONSTS(43)*STATES(46) - CONSTS(52)*STATES(11) ALGBRC(1) = STATES(2)+30.0000 RATES(9) = 7.19000e-05*(ALGBRC(1)/(1.00000 - EXP( -0.148000*ALGBRC(1))))*(1.00000 - STATES(9)) - 0.000131000*(ALGBRC(1)/(EXP( 0.0687000*ALGBRC(1)) - 1.00000))*STATES(9) ALGBRC(4) = TERNRY(STATES(2).LT.-40.0000, ( 0.135000*EXP(-80.0000/6.80000))*EXP( (-1.00000/6.80000)*STATES(2)), 0.00000) ALGBRC(10) = TERNRY(STATES(2).LT.-40.0000, 3.56000*EXP( 0.0790000*STATES(2))+ 310000.*EXP( 0.350000*STATES(2)), 1.00000/(0.130000+ ( 0.130000*EXP(- (10.6600/11.1000)))*EXP( STATES(2)*(-1.00000/11.1000)))) RATES(7) = ALGBRC(4)*(1.00000 - STATES(7)) - ALGBRC(10)*STATES(7) ALGBRC(5) = TERNRY(STATES(2).LT.-40.0000, ( -127140.*EXP( 0.244400*STATES(2)) - 3.47400e-05*EXP( -0.0439100*STATES(2)))*((STATES(2)+37.7800)/(1.00000+EXP( 0.311000*(STATES(2)+79.2300)))), 0.00000) ALGBRC(11) = TERNRY(STATES(2).LT.-40.0000, 0.121200*(EXP( -0.0105200*STATES(2))/(1.00000+EXP( -0.137800*(STATES(2)+40.1400)))), 0.300000*(EXP( -2.53500e-07*STATES(2))/(1.00000+EXP( -0.100000*STATES(2) - 3.20000)))) RATES(8) = ALGBRC(5)*(1.00000 - STATES(8)) - ALGBRC(11)*STATES(8) ALGBRC(6) = CONSTS(138)* STATES(35)*CONSTS(314) ** CONSTS(252) ALGBRC(12) = CONSTS(310)*STATES(31)+ ALGBRC(6)*STATES(51)+ CONSTS(334)*STATES(32) RATES(31) = ALGBRC(12) ALGBRC(15) = CONSTS(341)*STATES(32)+ ALGBRC(6)*STATES(30)+ CONSTS(249)*STATES(31)+ CONSTS(335)*STATES(33) RATES(32) = ALGBRC(15) ALGBRC(3) = TERNRY(STATES(2).EQ.-47.1300, 3.20000, 0.320000*((STATES(2)+47.1300)/(1.00000 - EXP( -0.100000*(STATES(2)+47.1300))))) ALGBRC(9) = 0.0800000*EXP(- ( STATES(2)*(1.00000/11.0000))) ALGBRC(14) = 1.00000/(ALGBRC(9)+ALGBRC(3)) ALGBRC(18) = TERNRY(ALGBRC(14).LT.0.0300000, ALGBRC(3)*ALGBRC(14), STATES(6)) RATES(6) = ALGBRC(3)*(1.00000 - ALGBRC(18)) - ALGBRC(9)*ALGBRC(18) ALGBRC(19) = CONSTS(342)*STATES(33)+ CONSTS(305)*STATES(32)+ CONSTS(336)*STATES(34) RATES(33) = ALGBRC(19) ALGBRC(22) = - ( CONSTS(336)*STATES(34))+ CONSTS(313)*STATES(33) RATES(34) = ALGBRC(22) ALGBRC(25) = CONSTS(138)*STATES(32) - (ALGBRC(6)+CONSTS(340))*STATES(30) RATES(30) = ALGBRC(25) RATES(51) = - ALGBRC(12) - (ALGBRC(15)+ALGBRC(19)+ALGBRC(22)+ALGBRC(25)) ALGBRC(30) = 1.00000 - (STATES(3)+STATES(4)+STATES(5)) RATES(3) = - ( CONSTS(65)*STATES(29) ** CONSTS(146)*STATES(3))+ CONSTS(64)*ALGBRC(30) RATES(5) = CONSTS(67)*STATES(29) ** CONSTS(141)*ALGBRC(30) - CONSTS(66)*STATES(5) RATES(4) = CONSTS(71)*ALGBRC(30) - CONSTS(69)*STATES(4) ALGBRC(7) = 0.400000*EXP( (STATES(2)+2.00000)*0.100000) ALGBRC(24) = 4.00000*ALGBRC(7) ALGBRC(13) = 0.0500000*EXP( (STATES(2)+2.00000)*(-1.00000/13.0000)) ALGBRC(41) = ALGBRC(13) ALGBRC(57) = 0.187500*STATES(29) ALGBRC(59) = ALGBRC(57) RATES(39) = ALGBRC(41)*STATES(37)+( CONSTS(312)*STATES(38) - (ALGBRC(24)+ALGBRC(59))*STATES(39)) ALGBRC(17) = CONSTS(19)*ALGBRC(7) ALGBRC(33) = 4.00000*ALGBRC(17) ALGBRC(21) = ALGBRC(13)*CONSTS(303) ALGBRC(49) = ALGBRC(21) RATES(38) = ALGBRC(49)*STATES(41)+( ALGBRC(59)*STATES(39) - (ALGBRC(33)+CONSTS(312))*STATES(38)) ALGBRC(27) = 3.00000*ALGBRC(7) ALGBRC(43) = 2.00000*ALGBRC(13) ALGBRC(61) = CONSTS(19)*ALGBRC(59) RATES(37) = ALGBRC(24)*STATES(39)+ ALGBRC(43)*STATES(40)+( CONSTS(316)*STATES(41) - (ALGBRC(41)+ALGBRC(27)+ALGBRC(61))*STATES(37)) ALGBRC(35) = 3.00000*ALGBRC(17) ALGBRC(51) = 2.00000*ALGBRC(21) RATES(41) = ALGBRC(33)*STATES(38)+ ALGBRC(51)*STATES(43)+( ALGBRC(61)*STATES(37) - (ALGBRC(49)+ALGBRC(35)+CONSTS(316))*STATES(41)) ALGBRC(29) = 2.00000*ALGBRC(7) ALGBRC(45) = 3.00000*ALGBRC(13) ALGBRC(63) = CONSTS(19)*ALGBRC(61) RATES(40) = ALGBRC(27)*STATES(37)+ ALGBRC(45)*STATES(42)+( CONSTS(321)*STATES(43) - (ALGBRC(43)+ALGBRC(29)+ALGBRC(63))*STATES(40)) ALGBRC(37) = 2.00000*ALGBRC(17) ALGBRC(53) = 3.00000*ALGBRC(21) RATES(43) = ALGBRC(35)*STATES(41)+ ALGBRC(53)*STATES(45)+( ALGBRC(63)*STATES(40) - (ALGBRC(51)+ALGBRC(37)+CONSTS(321))*STATES(43)) ALGBRC(31) = ALGBRC(7) ALGBRC(47) = 4.00000*ALGBRC(13) ALGBRC(65) = CONSTS(19)*ALGBRC(63) RATES(42) = ALGBRC(29)*STATES(40)+ ALGBRC(47)*STATES(44)+( CONSTS(326)*STATES(45) - (ALGBRC(45)+ALGBRC(31)+ALGBRC(65))*STATES(42)) ALGBRC(39) = ALGBRC(17) ALGBRC(55) = 4.00000*ALGBRC(21) RATES(45) = ALGBRC(37)*STATES(43)+ ALGBRC(55)*STATES(46)+( ALGBRC(65)*STATES(42) - (ALGBRC(53)+ALGBRC(39)+CONSTS(326))*STATES(45)) ALGBRC(67) = CONSTS(19)*ALGBRC(65) RATES(44) = ALGBRC(31)*STATES(42)+ CONSTS(50)*STATES(10)+( CONSTS(329)*STATES(46) - (ALGBRC(47)+CONSTS(41)+ALGBRC(67))*STATES(44)) RATES(46) = ALGBRC(39)*STATES(45)+ CONSTS(52)*STATES(11)+( ALGBRC(67)*STATES(44) - (ALGBRC(55)+CONSTS(43)+CONSTS(329))*STATES(46)) ALGBRC(72) = CONSTS(75)*(CONSTS(24) - (STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(19)+ STATES(25)*CONSTS(267))) ALGBRC(66) = CONSTS(28) - STATES(18) ALGBRC(68) = CONSTS(113)/ALGBRC(66) ALGBRC(73) = 1.00000/( (1.00000+ STATES(16)*CONSTS(270))*(1.00000+ STATES(26)*CONSTS(271))) ALGBRC(74) = 1.00000+ STATES(18)*CONSTS(272) ALGBRC(75) = CONSTS(269)/(CONSTS(268)+ ALGBRC(68)*ALGBRC(74)+ CONSTS(114)/STATES(25) ** CONSTS(148)*ALGBRC(73)*(1.00000+ ALGBRC(68)*ALGBRC(74))) RATES(25) = ALGBRC(72) - ALGBRC(75) ALGBRC(76) = CONSTS(274)+ CONSTS(309)*STATES(26) ALGBRC(77) = CONSTS(275)*(ALGBRC(76)/(ALGBRC(76)+CONSTS(115)/STATES(20) ** CONSTS(149)+ CONSTS(273)*ALGBRC(68))) ALGBRC(62) = CONSTS(25) - STATES(16) ALGBRC(78) = CONSTS(79)*( STATES(21)*STATES(16) - CONSTS(276)*STATES(22)*ALGBRC(62)) RATES(21) = ALGBRC(77) - ALGBRC(78) ALGBRC(79) = CONSTS(277)*(STATES(22)/(STATES(22)+ (CONSTS(124)+ CONSTS(278)*STATES(23))*(1.00000+ CONSTS(279)*STATES(19)))) RATES(22) = ALGBRC(78) - ALGBRC(79) ALGBRC(80) = CONSTS(77)*STATES(23) - CONSTS(280)*STATES(24) RATES(23) = ALGBRC(79) - ALGBRC(80) ALGBRC(81) = CONSTS(281)*STATES(24)*(ALGBRC(66)/( (STATES(24)+CONSTS(105)+ STATES(19)*CONSTS(282))*(CONSTS(118)+ALGBRC(66)))) RATES(24) = ALGBRC(80) - ALGBRC(81) ALGBRC(82) = CONSTS(283)*(STATES(19)/(CONSTS(284)+STATES(20))) RATES(20) = ALGBRC(75)+(ALGBRC(82) - ALGBRC(77)) ALGBRC(71) = CONSTS(266)*(STATES(19)/(STATES(19)+CONSTS(122))) RATES(19) = ALGBRC(81) - (ALGBRC(71)+ALGBRC(82)) RATES(50) = ALGBRC(82) - CONSTS(70)*STATES(50) ALGBRC(64) = CONSTS(308)+STATES(17) ALGBRC(69) = EXP( CONSTS(343)*ALGBRC(64)) ALGBRC(83) = ( STATES(18)*CONSTS(287)*ALGBRC(68)) ** (1.0 / 2) ALGBRC(84) = 1.00000/((CONSTS(338)+ CONSTS(347)*ALGBRC(83))+ (CONSTS(165)+ CONSTS(166)*ALGBRC(83))*ALGBRC(69)) ALGBRC(85) = ( (CONSTS(349)+ CONSTS(351)*ALGBRC(69))*ALGBRC(83) - CONSTS(353)*ALGBRC(69))*ALGBRC(84) RATES(18) = - ALGBRC(85)+ALGBRC(75)+ALGBRC(77)+ALGBRC(81) ALGBRC(89) = EXP( CONSTS(373)*ALGBRC(64)) ALGBRC(90) = CONSTS(288)*(ALGBRC(62)/STATES(16)) ALGBRC(91) = - (CONSTS(171)/(CONSTS(372)+ CONSTS(377)*ALGBRC(90)+ (CONSTS(153)+ CONSTS(154)*ALGBRC(90))*ALGBRC(89))) ALGBRC(92) = ( (CONSTS(374)+ CONSTS(158)*ALGBRC(89))*ALGBRC(90) - CONSTS(155)*ALGBRC(89))*ALGBRC(91) ALGBRC(8) = 8.00000 - STATES(1) ALGBRC(94) = STATES(1)/ALGBRC(8) ALGBRC(95) = STATES(16)/ALGBRC(62) ALGBRC(96) = (CONSTS(289) - CONSTS(290)*ALGBRC(94)*ALGBRC(95)*EXP(- ( CONSTS(320)*STATES(17))))/( (1.00000+ (10.0000/9.00000)*ALGBRC(94)*EXP(- ( CONSTS(323)*STATES(17))))*(1.00000+ 18.0000*ALGBRC(95))) RATES(16) = ALGBRC(96) - (ALGBRC(92)+ALGBRC(78)) ALGBRC(16) = 1.00000/STATES(1) ALGBRC(98) = STATES(15)*CONSTS(291) ** CONSTS(147) ALGBRC(99) = STATES(27)*CONSTS(292) ** CONSTS(150) ALGBRC(100) = CONSTS(136)*(( CONSTS(180)*ALGBRC(98) - CONSTS(181)*ALGBRC(99))/( (1.00000+ALGBRC(98)+ALGBRC(99))*( ALGBRC(16)*(CONSTS(107)+ ALGBRC(8)*CONSTS(378))+(1.00000+ ALGBRC(8)*CONSTS(379))))) ALGBRC(103) = (STATES(27) - STATES(28))*CONSTS(293) RATES(27) = ALGBRC(100)*CONSTS(384) - ALGBRC(103)*CONSTS(385) ALGBRC(101) = CONSTS(10)*(STATES(47) - STATES(48)) ALGBRC(105) = CONSTS(12)*( (CONSTS(15) - STATES(47))*STATES(1) - STATES(47)*ALGBRC(8)*CONSTS(301)) RATES(47) = ALGBRC(105) - ALGBRC(101) ALGBRC(104) = CONSTS(11)*( (CONSTS(14) - STATES(48))*STATES(49) - STATES(48)*(8.00000 - STATES(49))*CONSTS(301)) RATES(48) = ALGBRC(101)+ALGBRC(104) ALGBRC(20) = STATES(2)*CONSTS(320) ALGBRC(23) = EXP(ALGBRC(20)) ALGBRC(38) = ( STATES(14)*STATES(14)*STATES(14)) ** (1.0 / 2) ALGBRC(40) = CONSTS(328)*(ALGBRC(38)/( (ALGBRC(38)+ ( CONSTS(120)*CONSTS(120)*CONSTS(120)) ** (1.0 / 2))*(1.00000+ 0.124500*EXP( -0.100000*ALGBRC(20))+CONSTS(325)/ALGBRC(23))*(1.00000+ ( CONSTS(101)*ALGBRC(16))*(1.00000+ ALGBRC(8)*CONSTS(332))))) ALGBRC(48) = CONSTS(58)*(STATES(15)/(CONSTS(123)+STATES(15)))*(1.00000/(1.00000+ ( CONSTS(102)*ALGBRC(16))*(1.00000+ ALGBRC(8)*CONSTS(263)))+1.00000/(1.00000+ CONSTS(103)*ALGBRC(16))) ALGBRC(60) = CONSTS(319)*(( CONSTS(249)*STATES(31)+ CONSTS(305)*STATES(32)+ CONSTS(313)*STATES(33))/(1.00000+ ALGBRC(16)*(CONSTS(106)+ CONSTS(265)*ALGBRC(8)))) RATES(1) = (1.00000 - CONSTS(4))*( 0.615000*ALGBRC(96) - (ALGBRC(60)+ 0.500000*ALGBRC(100)+ 6.37100e-05*(ALGBRC(40)+ALGBRC(48))+ALGBRC(105))) RATES(49) = (1.00000 - CONSTS(5))*(- ALGBRC(104) - CONSTS(16)) ALGBRC(26) = CONSTS(315)*ALGBRC(20) ALGBRC(28) = ALGBRC(23)*ALGBRC(23) ALGBRC(32) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), CONSTS(257)*(ALGBRC(28) - CONSTS(255))*(0.500000 - 0.0200000*STATES(2)), CONSTS(257)*ALGBRC(26)*((ALGBRC(28) - CONSTS(255))/(ALGBRC(28) - 1.00000))) ALGBRC(36) = 6.00000*ALGBRC(32)*STATES(12)*STATES(10) ALGBRC(102) = CONSTS(178)*(ALGBRC(30)+STATES(5))*(STATES(28) - STATES(29)) ALGBRC(106) = (STATES(29) - STATES(15))*CONSTS(294) ALGBRC(107) = 1.00000/(1.00000+CONSTS(295)/( (STATES(29)+CONSTS(108))*(STATES(29)+CONSTS(108)))) RATES(29) = ALGBRC(107)*( ALGBRC(102)*CONSTS(382) - ( ALGBRC(106)*CONSTS(383)+ ALGBRC(36)*CONSTS(381))) ALGBRC(108) = 1.00000/(1.00000+CONSTS(296)/( (STATES(28)+CONSTS(110))*(STATES(28)+CONSTS(110)))) RATES(28) = ALGBRC(108)*(ALGBRC(103) - ALGBRC(102)) ALGBRC(109) = STATES(32)+STATES(30)+STATES(33)+STATES(34) ALGBRC(112) = CONSTS(366)*ALGBRC(109) ALGBRC(113) = CONSTS(99)*STATES(15)*(CONSTS(140) - STATES(35)) - ( CONSTS(100)*STATES(35))*(1.00000 - (2.00000/3.00000)*ALGBRC(112)) RATES(35) = ALGBRC(113) ALGBRC(114) = CONSTS(87)*STATES(15)*(CONSTS(55) - STATES(36)) - CONSTS(86)*STATES(36) RATES(36) = ALGBRC(114) ALGBRC(70) = CONSTS(324)*(STATES(17) - 91.0000) ALGBRC(117) = 1.00000+ STATES(15)*CONSTS(297) ALGBRC(118) = ALGBRC(117)*ALGBRC(117)*ALGBRC(117) ALGBRC(119) = CONSTS(299)*STATES(15)*ALGBRC(70)*(ALGBRC(118)/( ( ALGBRC(118)*ALGBRC(117)+CONSTS(139)/1.00000+ STATES(15)*CONSTS(298) ** CONSTS(144))*(1.00000 - EXP(- ALGBRC(70))))) ALGBRC(120) = CONSTS(183)*EXP( CONSTS(300)*ALGBRC(70))*(STATES(26)/( STATES(15)*1.00000+CONSTS(125)/STATES(14) ** CONSTS(143)*(1.00000+CONSTS(127)/STATES(26)))) RATES(26) = CONSTS(42)*(ALGBRC(119) - ALGBRC(120)) ALGBRC(86) = ( CONSTS(355)*ALGBRC(83) - CONSTS(357)*ALGBRC(69))*ALGBRC(84) ALGBRC(87) = CONSTS(172)/(CONSTS(370)+ CONSTS(371)*ALGBRC(69)) ALGBRC(88) = (CONSTS(365) - CONSTS(368)*ALGBRC(69))*ALGBRC(87) ALGBRC(93) = (CONSTS(376)+( CONSTS(376)*ALGBRC(90) - CONSTS(375)*ALGBRC(89)))*ALGBRC(91) ALGBRC(97) = CONSTS(49)*ALGBRC(64) RATES(17) = - ( ((- ALGBRC(86) - ALGBRC(88))+ALGBRC(93)+ALGBRC(96)+ALGBRC(97)+ CONSTS(388)*ALGBRC(120)+ 2.00000*ALGBRC(119))*CONSTS(387)) ALGBRC(121) = (( CONSTS(7)*CONSTS(8))/( CONSTS(192)*CONSTS(6)))*log(CONSTS(145)/STATES(14)) ALGBRC(122) = CONSTS(47)*STATES(6)*STATES(6)*STATES(6)*STATES(7)*STATES(8)*(STATES(2) - ALGBRC(121)) ALGBRC(123) = CONSTS(48)*(STATES(2) - ALGBRC(121)) ALGBRC(42) = EXP( CONSTS(33)*ALGBRC(20)) ALGBRC(44) = ALGBRC(42)/ALGBRC(23) ALGBRC(46) = ( ALGBRC(42)*STATES(14)*STATES(14)*STATES(14) - ALGBRC(44)*CONSTS(261)*STATES(15))/(CONSTS(262)+ CONSTS(307)*ALGBRC(44)) ALGBRC(50) = STATES(15)*STATES(15)*STATES(15) ALGBRC(52) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), 0.750000*ALGBRC(50)*((1.00000 - 0.0200000*STATES(2))/(ALGBRC(50)+CONSTS(264))), 0.750000*ALGBRC(50)*(ALGBRC(26)/( (ALGBRC(23) - 1.00000)*(ALGBRC(50)+CONSTS(264))))) ALGBRC(54) = CONSTS(163)*ALGBRC(52)*( STATES(14)*ALGBRC(23) - CONSTS(145)) RATES(14) = - ( (ALGBRC(122)+ALGBRC(123)+ALGBRC(54)+ 3.00000*(ALGBRC(46)+ALGBRC(40)))*CONSTS(380)) - ALGBRC(120)*0.615000 ALGBRC(131) = (( CONSTS(7)*CONSTS(8))/( CONSTS(194)*CONSTS(6)))*log(CONSTS(23)/STATES(15)) ALGBRC(132) = CONSTS(45)*(STATES(2) - ALGBRC(131)) ALGBRC(115) = ALGBRC(113)+ALGBRC(114) ALGBRC(116) = 1.00000/(1.00000+CONSTS(295)/( (STATES(15)+CONSTS(108))*(STATES(15)+CONSTS(108)))) RATES(15) = ALGBRC(116)*((ALGBRC(106) - (ALGBRC(100)+ALGBRC(115)+ 0.250000*CONSTS(380)*((ALGBRC(132) - 2.00000*ALGBRC(46))+ALGBRC(48))))+ (ALGBRC(120) - ALGBRC(119))*0.615000) ALGBRC(2) = TERNRY( INT(MOD(VOI, CONSTS(1))).LT.CONSTS(2), CONSTS(3), 0.00000) ALGBRC(133) = (( CONSTS(7)*CONSTS(8))/( CONSTS(193)*CONSTS(6)))*log(( CONSTS(195)*CONSTS(128)+ CONSTS(196)*CONSTS(145))/( CONSTS(195)*STATES(13)+ CONSTS(196)*STATES(14))) ALGBRC(134) = CONSTS(253)*STATES(9)*STATES(9)*((STATES(2) - ALGBRC(133))/(1.00000+EXP((STATES(2) - 40.0000)/40.0000))) ALGBRC(124) = (( CONSTS(7)*CONSTS(8))/( CONSTS(193)*CONSTS(6)))*log(CONSTS(128)/STATES(13)) ALGBRC(125) = STATES(2) - ALGBRC(124) ALGBRC(126) = 1.02000/(1.00000+EXP( 0.238500*(ALGBRC(125) - 59.2150))) ALGBRC(127) = ( 0.491200*EXP( 0.0803200*(ALGBRC(125)+5.47600))+EXP( 0.0617500*(ALGBRC(125) - 594.310)))/(1.00000+EXP( -0.514300*((STATES(2) - ALGBRC(124))+4.75300))) ALGBRC(128) = ALGBRC(126)/(ALGBRC(126)+ALGBRC(127)) ALGBRC(129) = CONSTS(311)*ALGBRC(128)*ALGBRC(125) ALGBRC(130) = CONSTS(46)*(ALGBRC(125)/(1.00000+EXP( (7.48800 - STATES(2))*CONSTS(254)))) ALGBRC(34) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), CONSTS(318)*(STATES(10)+STATES(11))*STATES(12)*( STATES(13)*ALGBRC(28) - CONSTS(128))*((0.500000 - 0.0200000*STATES(2))/(1.00000+ ALGBRC(32)*CONSTS(258))), CONSTS(318)*(STATES(10)+STATES(11))*STATES(12)*( STATES(13)*ALGBRC(28) - CONSTS(128))*(ALGBRC(20)/( (ALGBRC(28) - 1.00000)*(1.00000+ ALGBRC(32)*CONSTS(258))))) ALGBRC(56) = CONSTS(162)*ALGBRC(52)*( STATES(13)*ALGBRC(23) - CONSTS(128)) RATES(13) = - ( (ALGBRC(56)+ALGBRC(134)+ALGBRC(129)+ALGBRC(130)+ALGBRC(34)+(ALGBRC(2) - 2.00000*ALGBRC(40)))*CONSTS(380)) ALGBRC(58) = ALGBRC(54)+ALGBRC(56) RATES(2) = - ( CONSTS(386)*(ALGBRC(122)+ALGBRC(36)+ALGBRC(34)+ALGBRC(134)+ALGBRC(129)+ALGBRC(130)+ALGBRC(46)+ALGBRC(40)+ALGBRC(58)+ALGBRC(48)+ALGBRC(132)+ALGBRC(123)+ALGBRC(2))) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = STATES(2)+30.0000 ALGBRC(4) = TERNRY(STATES(2).LT.-40.0000, ( 0.135000*EXP(-80.0000/6.80000))*EXP( (-1.00000/6.80000)*STATES(2)), 0.00000) ALGBRC(10) = TERNRY(STATES(2).LT.-40.0000, 3.56000*EXP( 0.0790000*STATES(2))+ 310000.*EXP( 0.350000*STATES(2)), 1.00000/(0.130000+ ( 0.130000*EXP(- (10.6600/11.1000)))*EXP( STATES(2)*(-1.00000/11.1000)))) ALGBRC(5) = TERNRY(STATES(2).LT.-40.0000, ( -127140.*EXP( 0.244400*STATES(2)) - 3.47400e-05*EXP( -0.0439100*STATES(2)))*((STATES(2)+37.7800)/(1.00000+EXP( 0.311000*(STATES(2)+79.2300)))), 0.00000) ALGBRC(11) = TERNRY(STATES(2).LT.-40.0000, 0.121200*(EXP( -0.0105200*STATES(2))/(1.00000+EXP( -0.137800*(STATES(2)+40.1400)))), 0.300000*(EXP( -2.53500e-07*STATES(2))/(1.00000+EXP( -0.100000*STATES(2) - 3.20000)))) ALGBRC(6) = CONSTS(138)* STATES(35)*CONSTS(314) ** CONSTS(252) ALGBRC(12) = CONSTS(310)*STATES(31)+ ALGBRC(6)*STATES(51)+ CONSTS(334)*STATES(32) ALGBRC(15) = CONSTS(341)*STATES(32)+ ALGBRC(6)*STATES(30)+ CONSTS(249)*STATES(31)+ CONSTS(335)*STATES(33) ALGBRC(3) = TERNRY(STATES(2).EQ.-47.1300, 3.20000, 0.320000*((STATES(2)+47.1300)/(1.00000 - EXP( -0.100000*(STATES(2)+47.1300))))) ALGBRC(9) = 0.0800000*EXP(- ( STATES(2)*(1.00000/11.0000))) ALGBRC(14) = 1.00000/(ALGBRC(9)+ALGBRC(3)) ALGBRC(18) = TERNRY(ALGBRC(14).LT.0.0300000, ALGBRC(3)*ALGBRC(14), STATES(6)) ALGBRC(19) = CONSTS(342)*STATES(33)+ CONSTS(305)*STATES(32)+ CONSTS(336)*STATES(34) ALGBRC(22) = - ( CONSTS(336)*STATES(34))+ CONSTS(313)*STATES(33) ALGBRC(25) = CONSTS(138)*STATES(32) - (ALGBRC(6)+CONSTS(340))*STATES(30) ALGBRC(30) = 1.00000 - (STATES(3)+STATES(4)+STATES(5)) ALGBRC(7) = 0.400000*EXP( (STATES(2)+2.00000)*0.100000) ALGBRC(24) = 4.00000*ALGBRC(7) ALGBRC(13) = 0.0500000*EXP( (STATES(2)+2.00000)*(-1.00000/13.0000)) ALGBRC(41) = ALGBRC(13) ALGBRC(57) = 0.187500*STATES(29) ALGBRC(59) = ALGBRC(57) ALGBRC(17) = CONSTS(19)*ALGBRC(7) ALGBRC(33) = 4.00000*ALGBRC(17) ALGBRC(21) = ALGBRC(13)*CONSTS(303) ALGBRC(49) = ALGBRC(21) ALGBRC(27) = 3.00000*ALGBRC(7) ALGBRC(43) = 2.00000*ALGBRC(13) ALGBRC(61) = CONSTS(19)*ALGBRC(59) ALGBRC(35) = 3.00000*ALGBRC(17) ALGBRC(51) = 2.00000*ALGBRC(21) ALGBRC(29) = 2.00000*ALGBRC(7) ALGBRC(45) = 3.00000*ALGBRC(13) ALGBRC(63) = CONSTS(19)*ALGBRC(61) ALGBRC(37) = 2.00000*ALGBRC(17) ALGBRC(53) = 3.00000*ALGBRC(21) ALGBRC(31) = ALGBRC(7) ALGBRC(47) = 4.00000*ALGBRC(13) ALGBRC(65) = CONSTS(19)*ALGBRC(63) ALGBRC(39) = ALGBRC(17) ALGBRC(55) = 4.00000*ALGBRC(21) ALGBRC(67) = CONSTS(19)*ALGBRC(65) ALGBRC(72) = CONSTS(75)*(CONSTS(24) - (STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(19)+ STATES(25)*CONSTS(267))) ALGBRC(66) = CONSTS(28) - STATES(18) ALGBRC(68) = CONSTS(113)/ALGBRC(66) ALGBRC(73) = 1.00000/( (1.00000+ STATES(16)*CONSTS(270))*(1.00000+ STATES(26)*CONSTS(271))) ALGBRC(74) = 1.00000+ STATES(18)*CONSTS(272) ALGBRC(75) = CONSTS(269)/(CONSTS(268)+ ALGBRC(68)*ALGBRC(74)+ CONSTS(114)/STATES(25) ** CONSTS(148)*ALGBRC(73)*(1.00000+ ALGBRC(68)*ALGBRC(74))) ALGBRC(76) = CONSTS(274)+ CONSTS(309)*STATES(26) ALGBRC(77) = CONSTS(275)*(ALGBRC(76)/(ALGBRC(76)+CONSTS(115)/STATES(20) ** CONSTS(149)+ CONSTS(273)*ALGBRC(68))) ALGBRC(62) = CONSTS(25) - STATES(16) ALGBRC(78) = CONSTS(79)*( STATES(21)*STATES(16) - CONSTS(276)*STATES(22)*ALGBRC(62)) ALGBRC(79) = CONSTS(277)*(STATES(22)/(STATES(22)+ (CONSTS(124)+ CONSTS(278)*STATES(23))*(1.00000+ CONSTS(279)*STATES(19)))) ALGBRC(80) = CONSTS(77)*STATES(23) - CONSTS(280)*STATES(24) ALGBRC(81) = CONSTS(281)*STATES(24)*(ALGBRC(66)/( (STATES(24)+CONSTS(105)+ STATES(19)*CONSTS(282))*(CONSTS(118)+ALGBRC(66)))) ALGBRC(82) = CONSTS(283)*(STATES(19)/(CONSTS(284)+STATES(20))) ALGBRC(71) = CONSTS(266)*(STATES(19)/(STATES(19)+CONSTS(122))) ALGBRC(64) = CONSTS(308)+STATES(17) ALGBRC(69) = EXP( CONSTS(343)*ALGBRC(64)) ALGBRC(83) = ( STATES(18)*CONSTS(287)*ALGBRC(68)) ** (1.0 / 2) ALGBRC(84) = 1.00000/((CONSTS(338)+ CONSTS(347)*ALGBRC(83))+ (CONSTS(165)+ CONSTS(166)*ALGBRC(83))*ALGBRC(69)) ALGBRC(85) = ( (CONSTS(349)+ CONSTS(351)*ALGBRC(69))*ALGBRC(83) - CONSTS(353)*ALGBRC(69))*ALGBRC(84) ALGBRC(89) = EXP( CONSTS(373)*ALGBRC(64)) ALGBRC(90) = CONSTS(288)*(ALGBRC(62)/STATES(16)) ALGBRC(91) = - (CONSTS(171)/(CONSTS(372)+ CONSTS(377)*ALGBRC(90)+ (CONSTS(153)+ CONSTS(154)*ALGBRC(90))*ALGBRC(89))) ALGBRC(92) = ( (CONSTS(374)+ CONSTS(158)*ALGBRC(89))*ALGBRC(90) - CONSTS(155)*ALGBRC(89))*ALGBRC(91) ALGBRC(8) = 8.00000 - STATES(1) ALGBRC(94) = STATES(1)/ALGBRC(8) ALGBRC(95) = STATES(16)/ALGBRC(62) ALGBRC(96) = (CONSTS(289) - CONSTS(290)*ALGBRC(94)*ALGBRC(95)*EXP(- ( CONSTS(320)*STATES(17))))/( (1.00000+ (10.0000/9.00000)*ALGBRC(94)*EXP(- ( CONSTS(323)*STATES(17))))*(1.00000+ 18.0000*ALGBRC(95))) ALGBRC(16) = 1.00000/STATES(1) ALGBRC(98) = STATES(15)*CONSTS(291) ** CONSTS(147) ALGBRC(99) = STATES(27)*CONSTS(292) ** CONSTS(150) ALGBRC(100) = CONSTS(136)*(( CONSTS(180)*ALGBRC(98) - CONSTS(181)*ALGBRC(99))/( (1.00000+ALGBRC(98)+ALGBRC(99))*( ALGBRC(16)*(CONSTS(107)+ ALGBRC(8)*CONSTS(378))+(1.00000+ ALGBRC(8)*CONSTS(379))))) ALGBRC(103) = (STATES(27) - STATES(28))*CONSTS(293) ALGBRC(101) = CONSTS(10)*(STATES(47) - STATES(48)) ALGBRC(105) = CONSTS(12)*( (CONSTS(15) - STATES(47))*STATES(1) - STATES(47)*ALGBRC(8)*CONSTS(301)) ALGBRC(104) = CONSTS(11)*( (CONSTS(14) - STATES(48))*STATES(49) - STATES(48)*(8.00000 - STATES(49))*CONSTS(301)) ALGBRC(20) = STATES(2)*CONSTS(320) ALGBRC(23) = EXP(ALGBRC(20)) ALGBRC(38) = ( STATES(14)*STATES(14)*STATES(14)) ** (1.0 / 2) ALGBRC(40) = CONSTS(328)*(ALGBRC(38)/( (ALGBRC(38)+ ( CONSTS(120)*CONSTS(120)*CONSTS(120)) ** (1.0 / 2))*(1.00000+ 0.124500*EXP( -0.100000*ALGBRC(20))+CONSTS(325)/ALGBRC(23))*(1.00000+ ( CONSTS(101)*ALGBRC(16))*(1.00000+ ALGBRC(8)*CONSTS(332))))) ALGBRC(48) = CONSTS(58)*(STATES(15)/(CONSTS(123)+STATES(15)))*(1.00000/(1.00000+ ( CONSTS(102)*ALGBRC(16))*(1.00000+ ALGBRC(8)*CONSTS(263)))+1.00000/(1.00000+ CONSTS(103)*ALGBRC(16))) ALGBRC(60) = CONSTS(319)*(( CONSTS(249)*STATES(31)+ CONSTS(305)*STATES(32)+ CONSTS(313)*STATES(33))/(1.00000+ ALGBRC(16)*(CONSTS(106)+ CONSTS(265)*ALGBRC(8)))) ALGBRC(26) = CONSTS(315)*ALGBRC(20) ALGBRC(28) = ALGBRC(23)*ALGBRC(23) ALGBRC(32) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), CONSTS(257)*(ALGBRC(28) - CONSTS(255))*(0.500000 - 0.0200000*STATES(2)), CONSTS(257)*ALGBRC(26)*((ALGBRC(28) - CONSTS(255))/(ALGBRC(28) - 1.00000))) ALGBRC(36) = 6.00000*ALGBRC(32)*STATES(12)*STATES(10) ALGBRC(102) = CONSTS(178)*(ALGBRC(30)+STATES(5))*(STATES(28) - STATES(29)) ALGBRC(106) = (STATES(29) - STATES(15))*CONSTS(294) ALGBRC(107) = 1.00000/(1.00000+CONSTS(295)/( (STATES(29)+CONSTS(108))*(STATES(29)+CONSTS(108)))) ALGBRC(108) = 1.00000/(1.00000+CONSTS(296)/( (STATES(28)+CONSTS(110))*(STATES(28)+CONSTS(110)))) ALGBRC(109) = STATES(32)+STATES(30)+STATES(33)+STATES(34) ALGBRC(112) = CONSTS(366)*ALGBRC(109) ALGBRC(113) = CONSTS(99)*STATES(15)*(CONSTS(140) - STATES(35)) - ( CONSTS(100)*STATES(35))*(1.00000 - (2.00000/3.00000)*ALGBRC(112)) ALGBRC(114) = CONSTS(87)*STATES(15)*(CONSTS(55) - STATES(36)) - CONSTS(86)*STATES(36) ALGBRC(70) = CONSTS(324)*(STATES(17) - 91.0000) ALGBRC(117) = 1.00000+ STATES(15)*CONSTS(297) ALGBRC(118) = ALGBRC(117)*ALGBRC(117)*ALGBRC(117) ALGBRC(119) = CONSTS(299)*STATES(15)*ALGBRC(70)*(ALGBRC(118)/( ( ALGBRC(118)*ALGBRC(117)+CONSTS(139)/1.00000+ STATES(15)*CONSTS(298) ** CONSTS(144))*(1.00000 - EXP(- ALGBRC(70))))) ALGBRC(120) = CONSTS(183)*EXP( CONSTS(300)*ALGBRC(70))*(STATES(26)/( STATES(15)*1.00000+CONSTS(125)/STATES(14) ** CONSTS(143)*(1.00000+CONSTS(127)/STATES(26)))) ALGBRC(86) = ( CONSTS(355)*ALGBRC(83) - CONSTS(357)*ALGBRC(69))*ALGBRC(84) ALGBRC(87) = CONSTS(172)/(CONSTS(370)+ CONSTS(371)*ALGBRC(69)) ALGBRC(88) = (CONSTS(365) - CONSTS(368)*ALGBRC(69))*ALGBRC(87) ALGBRC(93) = (CONSTS(376)+( CONSTS(376)*ALGBRC(90) - CONSTS(375)*ALGBRC(89)))*ALGBRC(91) ALGBRC(97) = CONSTS(49)*ALGBRC(64) ALGBRC(121) = (( CONSTS(7)*CONSTS(8))/( CONSTS(192)*CONSTS(6)))*log(CONSTS(145)/STATES(14)) ALGBRC(122) = CONSTS(47)*STATES(6)*STATES(6)*STATES(6)*STATES(7)*STATES(8)*(STATES(2) - ALGBRC(121)) ALGBRC(123) = CONSTS(48)*(STATES(2) - ALGBRC(121)) ALGBRC(42) = EXP( CONSTS(33)*ALGBRC(20)) ALGBRC(44) = ALGBRC(42)/ALGBRC(23) ALGBRC(46) = ( ALGBRC(42)*STATES(14)*STATES(14)*STATES(14) - ALGBRC(44)*CONSTS(261)*STATES(15))/(CONSTS(262)+ CONSTS(307)*ALGBRC(44)) ALGBRC(50) = STATES(15)*STATES(15)*STATES(15) ALGBRC(52) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), 0.750000*ALGBRC(50)*((1.00000 - 0.0200000*STATES(2))/(ALGBRC(50)+CONSTS(264))), 0.750000*ALGBRC(50)*(ALGBRC(26)/( (ALGBRC(23) - 1.00000)*(ALGBRC(50)+CONSTS(264))))) ALGBRC(54) = CONSTS(163)*ALGBRC(52)*( STATES(14)*ALGBRC(23) - CONSTS(145)) ALGBRC(131) = (( CONSTS(7)*CONSTS(8))/( CONSTS(194)*CONSTS(6)))*log(CONSTS(23)/STATES(15)) ALGBRC(132) = CONSTS(45)*(STATES(2) - ALGBRC(131)) ALGBRC(115) = ALGBRC(113)+ALGBRC(114) ALGBRC(116) = 1.00000/(1.00000+CONSTS(295)/( (STATES(15)+CONSTS(108))*(STATES(15)+CONSTS(108)))) ALGBRC(2) = TERNRY( INT(MOD(VOI, CONSTS(1))).LT.CONSTS(2), CONSTS(3), 0.00000) ALGBRC(133) = (( CONSTS(7)*CONSTS(8))/( CONSTS(193)*CONSTS(6)))*log(( CONSTS(195)*CONSTS(128)+ CONSTS(196)*CONSTS(145))/( CONSTS(195)*STATES(13)+ CONSTS(196)*STATES(14))) ALGBRC(134) = CONSTS(253)*STATES(9)*STATES(9)*((STATES(2) - ALGBRC(133))/(1.00000+EXP((STATES(2) - 40.0000)/40.0000))) ALGBRC(124) = (( CONSTS(7)*CONSTS(8))/( CONSTS(193)*CONSTS(6)))*log(CONSTS(128)/STATES(13)) ALGBRC(125) = STATES(2) - ALGBRC(124) ALGBRC(126) = 1.02000/(1.00000+EXP( 0.238500*(ALGBRC(125) - 59.2150))) ALGBRC(127) = ( 0.491200*EXP( 0.0803200*(ALGBRC(125)+5.47600))+EXP( 0.0617500*(ALGBRC(125) - 594.310)))/(1.00000+EXP( -0.514300*((STATES(2) - ALGBRC(124))+4.75300))) ALGBRC(128) = ALGBRC(126)/(ALGBRC(126)+ALGBRC(127)) ALGBRC(129) = CONSTS(311)*ALGBRC(128)*ALGBRC(125) ALGBRC(130) = CONSTS(46)*(ALGBRC(125)/(1.00000+EXP( (7.48800 - STATES(2))*CONSTS(254)))) ALGBRC(34) = TERNRY(ABS(STATES(2)).LT.CONSTS(9), CONSTS(318)*(STATES(10)+STATES(11))*STATES(12)*( STATES(13)*ALGBRC(28) - CONSTS(128))*((0.500000 - 0.0200000*STATES(2))/(1.00000+ ALGBRC(32)*CONSTS(258))), CONSTS(318)*(STATES(10)+STATES(11))*STATES(12)*( STATES(13)*ALGBRC(28) - CONSTS(128))*(ALGBRC(20)/( (ALGBRC(28) - 1.00000)*(1.00000+ ALGBRC(32)*CONSTS(258))))) ALGBRC(56) = CONSTS(162)*ALGBRC(52)*( STATES(13)*ALGBRC(23) - CONSTS(128)) ALGBRC(58) = ALGBRC(54)+ALGBRC(56) ALGBRC(110) = CONSTS(367)*(ALGBRC(109)+STATES(33)+STATES(34)+STATES(34)) ALGBRC(111) = CONSTS(369)*(ALGBRC(109)+STATES(33)+STATES(34)+STATES(34)) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END