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 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 (dimensionless).
C STATES(1) is ATPi in component model (dimensionless).
C ALGBRC(16) is inv_ATPi in component model (dimensionless).
C ALGBRC(20) is VF_over_RT in component model (dimensionless).
C STATES(2) is V in component model (dimensionless).
C ALGBRC(23) is exp_VF_over_RT in component model (dimensionless).
C ALGBRC(26) is VFsq_over_RT in component model (dimensionless).
C ALGBRC(28) is exp2VFRT in component model (dimensionless).
C ALGBRC(30) is O1_RyR in component model (dimensionless).
C STATES(3) is C1_RyR in component model (dimensionless).
C STATES(4) is C2_RyR in component model (dimensionless).
C STATES(5) is O2_RyR in component model (dimensionless).
C ALGBRC(1) is V_30 in component model (dimensionless).
C ALGBRC(125) is V_E_K in component model (dimensionless).
C ALGBRC(122) is INa in component model (dimensionless).
C STATES(6) is mNa in component model (dimensionless).
C STATES(7) is hNa in component model (dimensionless).
C STATES(8) is jNa in component model (dimensionless).
C ALGBRC(134) is IKs in component model (dimensionless).
C STATES(9) is xKs in component model (dimensionless).
C ALGBRC(126) is K1Alpha in component model (dimensionless).
C ALGBRC(127) is K1Beta in component model (dimensionless).
C ALGBRC(128) is K1_inf in component model (dimensionless).
C ALGBRC(129) is IK1 in component model (dimensionless).
C ALGBRC(123) is INab in component model (dimensionless).
C ALGBRC(130) is IKp in component model (dimensionless).
C ALGBRC(32) is ICamax in component model (dimensionless).
C ALGBRC(34) is ICaK in component model (dimensionless).
C STATES(10) is Open in component model (dimensionless).
C STATES(11) is OCa in component model (dimensionless).
C STATES(12) is yCa in component model (dimensionless).
C STATES(13) is Ki in component model (dimensionless).
C ALGBRC(36) is ICa in component model (dimensionless).
C ALGBRC(38) is NaiP1p5 in component model (dimensionless).
C STATES(14) is Nai in component model (dimensionless).
C ALGBRC(40) is INaK in component model (dimensionless).
C ALGBRC(42) is exp_eta_VF_over_RT in component model (dimensionless).
C ALGBRC(44) is exp_eta1_VF_over_RT in component model (dimensionless).
C ALGBRC(46) is INaCa in component model (dimensionless).
C STATES(15) is Cai in component model (dimensionless).
C ALGBRC(132) is ICab in component model (dimensionless).
C ALGBRC(48) is IpCa in component model (dimensionless).
C ALGBRC(50) is CaiP3 in component model (dimensionless).
C ALGBRC(52) is common in component model (dimensionless).
C ALGBRC(54) is InsNa in component model (dimensionless).
C ALGBRC(56) is InsK in component model (dimensionless).
C ALGBRC(58) is InsCa in component model (dimensionless).
C ALGBRC(60) is V_AM in component model (dimensionless).
C ALGBRC(62) is ATPm in component model (dimensionless).
C STATES(16) is ADPm in component model (dimensionless).
C ALGBRC(64) is DmuH in component model (dimensionless).
C STATES(17) is Dpsi in component model (dimensionless).
C ALGBRC(66) is NAD in component model (dimensionless).
C STATES(18) is NADH in component model (dimensionless).
C ALGBRC(68) is KmIDNAD_NAD in component model (dimensionless).
C ALGBRC(69) is exp_FRT_6_g_DmuH in component model (dimensionless).
C ALGBRC(70) is FRT2_Dpsi in component model (dimensionless).
C ALGBRC(71) is VCS in component model (dimensionless).
C STATES(19) is Oaa in component model (dimensionless).
C ALGBRC(72) is VACO in component model (dimensionless).
C STATES(20) is AKG in component model (dimensionless).
C STATES(21) is SCoA in component model (dimensionless).
C STATES(22) is Succ in component model (dimensionless).
C STATES(23) is FUM in component model (dimensionless).
C STATES(24) is MAL in component model (dimensionless).
C STATES(25) is ISOC in component model (dimensionless).
C ALGBRC(73) is Fa in component model (dimensionless).
C STATES(26) is Cam in component model (dimensionless).
C ALGBRC(74) is Fi in component model (dimensionless).
C ALGBRC(75) is VIDH in component model (dimensionless).
C ALGBRC(76) is a__1 in component model (dimensionless).
C ALGBRC(77) is VKGDH in component model (dimensionless).
C ALGBRC(78) is VSL in component model (dimensionless).
C ALGBRC(79) is VSDH in component model (dimensionless).
C ALGBRC(80) is VFH in component model (dimensionless).
C ALGBRC(81) is VMDH in component model (dimensionless).
C ALGBRC(82) is VAAT in component model (dimensionless).
C ALGBRC(83) is AREN in component model (dimensionless).
C ALGBRC(84) is denominator1 in component model (dimensionless).
C ALGBRC(87) is denominator2 in component model (dimensionless).
C ALGBRC(91) is denominator3 in component model (dimensionless).
C ALGBRC(85) is VNO in component model (dimensionless).
C ALGBRC(86) is VHNe in component model (dimensionless).
C ALGBRC(88) is VHFe in component model (dimensionless).
C ALGBRC(89) is exp_3FRT_DmuH in component model (dimensionless).
C ALGBRC(90) is AF1 in component model (dimensionless).
C ALGBRC(92) is VATPase in component model (dimensionless).
C ALGBRC(93) is Vhu in component model (dimensionless).
C ALGBRC(94) is ATPi_ADP in component model (dimensionless).
C ALGBRC(95) is ADPm_ATPm in component model (dimensionless).
C ALGBRC(96) is VANT in component model (dimensionless).
C ALGBRC(97) is Vhleak in component model (dimensionless).
C ALGBRC(3) is MAlpha in component model (dimensionless).
C ALGBRC(9) is MBeta in component model (dimensionless).
C ALGBRC(14) is inv_MBeta_MAlpha in component model (dimensionless).
C ALGBRC(18) is tmNa in component model (dimensionless).
C ALGBRC(4) is HAlpha in component model (dimensionless).
C ALGBRC(10) is HBeta in component model (dimensionless).
C ALGBRC(5) is JAlpha in component model (dimensionless).
C ALGBRC(11) is JBeta in component model (dimensionless).
C ALGBRC(98) is fb in component model (dimensionless).
C ALGBRC(99) is rb in component model (per_ms).
C STATES(27) is CaNSR in component model (dimensionless).
C ALGBRC(100) is Jup in component model (dimensionless).
C ALGBRC(102) is Jrel in component model (dimensionless).
C STATES(28) is CaJSR in component model (dimensionless).
C STATES(29) is CaSS in component model (dimensionless).
C ALGBRC(103) is Jtr in component model (dimensionless).
C ALGBRC(106) is Jxfer in component model (dimensionless).
C ALGBRC(109) is P1_N1_P2_P3 in component model (dimensionless).
C STATES(30) is N1 in component model (dimensionless).
C STATES(31) is P0 in component model (dimensionless).
C STATES(32) is P1 in component model (dimensionless).
C STATES(33) is P2 in component model (dimensionless).
C STATES(34) is P3 in component model (dimensionless).
C ALGBRC(25) is N1_differential in component model (dimensionless).
C ALGBRC(12) is P0_differential in component model (dimensionless).
C ALGBRC(15) is P1_differential in component model (dimensionless).
C ALGBRC(19) is P2_differential in component model (dimensionless).
C ALGBRC(22) is P3_differential in component model (dimensionless).
C ALGBRC(112) is FN_Ca in component model (dimensionless).
C ALGBRC(110) is force_norm in component model (dimensionless).
C ALGBRC(111) is force in component model (dimensionless).
C ALGBRC(6) is kTrop_np in component model (dimensionless).
C STATES(35) is LTRPNCa in component model (dimensionless).
C STATES(36) is HTRPNCa in component model (dimensionless).
C ALGBRC(113) is LTRPNCa_differential in component model (dimensionless).
C ALGBRC(114) is HTRPNCa_differential in component model (dimensionless).
C ALGBRC(115) is Jtrpn in component model (dimensionless).
C ALGBRC(107) is beta_SS in component model (dimensionless).
C ALGBRC(108) is beta_JSR in component model (dimensionless).
C ALGBRC(116) is beta_i in component model (dimensionless).
C ALGBRC(117) is Cai_ktrans_plus1 in component model (dimensionless).
C ALGBRC(118) is Cai_ktrans_plus1_p3 in component model (dimensionless).
C ALGBRC(119) is Vuni in component model (dimensionless).
C ALGBRC(120) is VnaCa in component model (dimensionless).
C ALGBRC(7) is alpha in component model (dimensionless).
C ALGBRC(13) is beta in component model (dimensionless).
C ALGBRC(17) is alpha_prime in component model (dimensionless).
C ALGBRC(21) is beta_prime in component model (dimensionless).
C ALGBRC(24) is C0_to_C1 in component model (dimensionless).
C ALGBRC(27) is C1_to_C2 in component model (dimensionless).
C ALGBRC(29) is C2_to_C3 in component model (dimensionless).
C ALGBRC(31) is C3_to_C4 in component model (dimensionless).
C ALGBRC(33) is CCa0_to_CCa1 in component model (dimensionless).
C ALGBRC(35) is CCa1_to_CCa2 in component model (dimensionless).
C ALGBRC(37) is CCa2_to_CCa3 in component model (dimensionless).
C ALGBRC(39) is CCa3_to_CCa4 in component model (dimensionless).
C ALGBRC(41) is C1_to_C0 in component model (dimensionless).
C ALGBRC(43) is C2_to_C1 in component model (dimensionless).
C ALGBRC(45) is C3_to_C2 in component model (dimensionless).
C ALGBRC(47) is C4_to_C3 in component model (dimensionless).
C ALGBRC(49) is CCa1_to_CCa0 in component model (dimensionless).
C ALGBRC(51) is CCa2_to_CCa1 in component model (dimensionless).
C ALGBRC(53) is CCa3_to_CCa2 in component model (dimensionless).
C ALGBRC(55) is CCa4_to_CCa3 in component model (dimensionless).
C ALGBRC(57) is gamma in component model (dimensionless).
C ALGBRC(59) is C0_to_CCa0 in component model (dimensionless).
C ALGBRC(61) is C1_to_CCa1 in component model (dimensionless).
C ALGBRC(63) is C2_to_CCa2 in component model (dimensionless).
C ALGBRC(65) is C3_to_CCa3 in component model (dimensionless).
C ALGBRC(67) is C4_to_CCa4 in component model (dimensionless).
C CONSTS(312) is CCa0_to_C0 in component model (dimensionless).
C CONSTS(316) is CCa1_to_C1 in component model (dimensionless).
C CONSTS(321) is CCa2_to_C2 in component model (dimensionless).
C CONSTS(326) is CCa3_to_C3 in component model (dimensionless).
C CONSTS(329) is CCa4_to_C4 in component model (dimensionless).
C STATES(37) is C1 in component model (dimensionless).
C STATES(38) is CCa0 in component model (dimensionless).
C STATES(39) is C0 in component model (dimensionless).
C STATES(40) is C2 in component model (dimensionless).
C STATES(41) is CCa1 in component model (dimensionless).
C STATES(42) is C3 in component model (dimensionless).
C STATES(43) is CCa2 in component model (dimensionless).
C STATES(44) is C4 in component model (dimensionless).
C STATES(45) is CCa3 in component model (dimensionless).
C STATES(46) is CCa4 in component model (dimensionless).
C ALGBRC(101) is Vt_CRP2 in component model (dimensionless).
C STATES(47) is CrPi_mito in component model (dimensionless).
C STATES(48) is CrPi_cyto in component model (dimensionless).
C ALGBRC(104) is VCK_cyto in component model (dimensionless).
C STATES(49) is ATPi_cyto in component model (dimensionless).
C ALGBRC(105) is VCK_mito in component model (dimensionless).
C STATES(50) is ASP in component model (dimensionless).
C STATES(51) is N0 in component model (dimensionless).
C ALGBRC(124) is reversal_potential in component one_ion (dimensionless).
C ALGBRC(121) is reversal_potential in component one_ion (dimensionless).
C ALGBRC(133) is reversal_potential in component two_ions (dimensionless).
C ALGBRC(131) is reversal_potential in component one_ion (dimensionless).
C RATES(6) is d/dt mNa in component model (dimensionless).
C RATES(7) is d/dt hNa in component model (dimensionless).
C RATES(8) is d/dt jNa in component model (dimensionless).
C RATES(9) is d/dt xKs in component model (dimensionless).
C RATES(31) is d/dt P0 in component model (dimensionless).
C RATES(32) is d/dt P1 in component model (dimensionless).
C RATES(33) is d/dt P2 in component model (dimensionless).
C RATES(34) is d/dt P3 in component model (dimensionless).
C RATES(30) is d/dt N1 in component model (dimensionless).
C RATES(51) is d/dt N0 in component model (dimensionless).
C RATES(35) is d/dt LTRPNCa in component model (dimensionless).
C RATES(36) is d/dt HTRPNCa in component model (dimensionless).
C RATES(14) is d/dt Nai in component model (dimensionless).
C RATES(13) is d/dt Ki in component model (dimensionless).
C RATES(15) is d/dt Cai in component model (dimensionless).
C RATES(29) is d/dt CaSS in component model (dimensionless).
C RATES(28) is d/dt CaJSR in component model (dimensionless).
C RATES(27) is d/dt CaNSR in component model (dimensionless).
C RATES(2) is d/dt V in component model (dimensionless).
C RATES(3) is d/dt C1_RyR in component model (dimensionless).
C RATES(5) is d/dt O2_RyR in component model (dimensionless).
C RATES(4) is d/dt C2_RyR in component model (dimensionless).
C RATES(39) is d/dt C0 in component model (dimensionless).
C RATES(37) is d/dt C1 in component model (dimensionless).
C RATES(40) is d/dt C2 in component model (dimensionless).
C RATES(42) is d/dt C3 in component model (dimensionless).
C RATES(44) is d/dt C4 in component model (dimensionless).
C RATES(10) is d/dt Open in component model (dimensionless).
C RATES(38) is d/dt CCa0 in component model (dimensionless).
C RATES(41) is d/dt CCa1 in component model (dimensionless).
C RATES(43) is d/dt CCa2 in component model (dimensionless).
C RATES(45) is d/dt CCa3 in component model (dimensionless).
C RATES(46) is d/dt CCa4 in component model (dimensionless).
C RATES(12) is d/dt yCa in component model (dimensionless).
C RATES(11) is d/dt OCa in component model (dimensionless).
C RATES(47) is d/dt CrPi_mito in component model (dimensionless).
C RATES(48) is d/dt CrPi_cyto in component model (dimensionless).
C RATES(1) is d/dt ATPi in component model (dimensionless).
C RATES(49) is d/dt ATPi_cyto in component model (dimensionless).
C RATES(26) is d/dt Cam in component model (dimensionless).
C RATES(16) is d/dt ADPm in component model (dimensionless).
C RATES(17) is d/dt Dpsi in component model (dimensionless).
C RATES(18) is d/dt NADH in component model (dimensionless).
C RATES(25) is d/dt ISOC in component model (dimensionless).
C RATES(20) is d/dt AKG in component model (dimensionless).
C RATES(21) is d/dt SCoA in component model (dimensionless).
C RATES(22) is d/dt Succ in component model (dimensionless).
C RATES(23) is d/dt FUM in component model (dimensionless).
C RATES(24) is d/dt MAL in component model (dimensionless).
C RATES(19) is d/dt Oaa in component model (dimensionless).
C RATES(50) is d/dt ASP in component model (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
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)
STATES(1) = CONSTS(197)
STATES(2) = CONSTS(198)
STATES(3) = CONSTS(199)
STATES(4) = CONSTS(200)
STATES(5) = CONSTS(201)
STATES(6) = CONSTS(202)
STATES(7) = CONSTS(203)
STATES(8) = CONSTS(204)
STATES(9) = CONSTS(205)
STATES(10) = CONSTS(206)
STATES(11) = CONSTS(207)
STATES(12) = CONSTS(208)
STATES(13) = CONSTS(209)
STATES(14) = CONSTS(210)
STATES(15) = CONSTS(211)
STATES(16) = CONSTS(212)
STATES(17) = CONSTS(213)
STATES(18) = CONSTS(214)
STATES(19) = CONSTS(215)
STATES(20) = CONSTS(216)
STATES(21) = CONSTS(217)
STATES(22) = CONSTS(218)
STATES(23) = CONSTS(219)
STATES(24) = CONSTS(220)
STATES(25) = CONSTS(221)
STATES(26) = CONSTS(222)
STATES(27) = CONSTS(223)
STATES(28) = CONSTS(224)
STATES(29) = CONSTS(225)
STATES(30) = CONSTS(226)
STATES(31) = CONSTS(227)
STATES(32) = CONSTS(228)
STATES(33) = CONSTS(229)
STATES(34) = CONSTS(230)
STATES(35) = CONSTS(231)
STATES(36) = CONSTS(232)
STATES(37) = CONSTS(233)
STATES(38) = CONSTS(234)
STATES(39) = CONSTS(235)
STATES(40) = CONSTS(236)
STATES(41) = CONSTS(237)
STATES(42) = CONSTS(238)
STATES(43) = CONSTS(239)
STATES(44) = CONSTS(240)
STATES(45) = CONSTS(241)
STATES(46) = CONSTS(242)
STATES(47) = CONSTS(243)
STATES(48) = CONSTS(244)
STATES(49) = CONSTS(245)
STATES(50) = CONSTS(246)
STATES(51) = CONSTS(247)
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