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 179 entries in the algebraic variable array.
C There are a total of 17 entries in each of the rate and state variable arrays.
C There are a total of 230 entries in the constant variable array.
C
C
C VOI is time in component environment (day).
C ALGBRC(11) is U_Pd_PdPs in component Q_Pd (mole_per_day).
C STATES(1) is Q_Pd in component Q_Pd (mole).
C ALGBRC(1) is C_Pd in component Q_Pd (molar).
C ALGBRC(9) is U_Pd_PdEx in component Q_Pd (mole_per_day).
C CONSTS(156) is P_Pd_InPd in component Q_Pd (mole_per_day).
C ALGBRC(10) is v_PdPs in component Q_Pd (mole_per_day).
C CONSTS(1) is v_PdPs_star in component Q_Pd (mole_per_gram_day).
C CONSTS(175) is M_Pd_PdPs in component Q_Pd (molar).
C CONSTS(2) is M_Pd_PdPs_star in component Q_Pd (molar).
C CONSTS(3) is Y_Pd_InPd in component Q_Pd (molar).
C CONSTS(4) is D_Pd in component Q_Pd (litre_per_day).
C CONSTS(155) is k_PdEx in component Q_Pd (first_order_rate_constant).
C CONSTS(5) is T_Pd in component Q_Pd (day).
C CONSTS(6) is T_Pd_star in component Q_Pd (day).
C STATES(2) is Q_Ma in component Q_Ma (gram).
C STATES(3) is Q_Mc in component Q_Mc (gram).
C CONSTS(7) is V_Ru in component model_parameters (litre).
C CONSTS(8) is k_SoEx in component model_parameters (first_order_rate_constant).
C ALGBRC(12) is C_Ps in component Q_Ps (molar).
C ALGBRC(49) is U_McPs_PsAm in component Q_Ps (mole_per_day).
C ALGBRC(37) is U_Ps_PsMa in component Q_Ps (mole_per_day).
C ALGBRC(50) is U_Ps_PsMc in component Q_Ps (mole_per_day).
C ALGBRC(36) is U_MaPs_PsAm in component Q_Ps (mole_per_day).
C CONSTS(9) is M_Ha_McMa in component Q_Ps (molar).
C STATES(4) is Q_Ps in component Q_Ps (mole).
C CONSTS(177) is P_Ps_InPs in component Q_Ps (mole_per_day).
C CONSTS(180) is P_Ps_SaPs in component Q_Ps (mole_per_day).
C ALGBRC(13) is P_Ps_PdPs in component Q_Ps (mole_per_day).
C ALGBRC(63) is P_Ps_MaMd in component Q_Ps (mole_per_day).
C ALGBRC(79) is P_Ps_McPs in component Q_Ps (mole_per_day).
C ALGBRC(14) is U_Ps_PsEx in component Q_Ps (mole_per_day).
C ALGBRC(77) is U_Mc_McPs in component Q_Ps (mole_per_day).
C ALGBRC(15) is v_Ma_PsAm in component Q_Ps (mole_per_day).
C ALGBRC(16) is v_Mc_PsAm in component Q_Ps (mole_per_day).
C ALGBRC(17) is v_PsMa in component Q_Ps (mole_per_day).
C ALGBRC(18) is v_PsMc in component Q_Ps (mole_per_day).
C CONSTS(10) is v_PsAm_star in component Q_Ps (mole_per_gram_day).
C CONSTS(11) is v_PsMa_star in component Q_Ps (mole_per_gram_day).
C CONSTS(12) is v_PsMc_star in component Q_Ps (mole_per_gram_day).
C CONSTS(157) is k_PsEx in component Q_Ps (first_order_rate_constant).
C CONSTS(13) is Y_Ps_InPs in component Q_Ps (molar).
C CONSTS(14) is Y_Ps_SaPs in component Q_Ps (molar).
C CONSTS(15) is Y_Ps_PdPs in component Q_Ps (dimensionless).
C CONSTS(16) is Y_Ps_MaMd in component Q_Ps (mole_per_gram).
C CONSTS(17) is Y_Ps_McPs in component Q_Ps (dimensionless).
C CONSTS(18) is D_Ps in component Q_Ps (litre_per_day).
C CONSTS(19) is M_Ps_PsAm in component Q_Ps (molar).
C CONSTS(20) is M_Ps_PsMa in component Q_Ps (molar).
C CONSTS(21) is M_Ps_PsMc in component Q_Ps (molar).
C CONSTS(22) is M_Ha_PsMa in component Q_Ps (molar).
C CONSTS(23) is M_Hc_PsMc in component Q_Ps (molar).
C CONSTS(24) is J_Ha_PsAm in component Q_Ps (molar).
C CONSTS(25) is J_Hc_PsAm in component Q_Ps (molar).
C ALGBRC(76) is U_Mc_McEg in component Q_Mc (gram_per_day).
C ALGBRC(62) is U_Ma_MaMd in component Q_Ma (gram_per_day).
C ALGBRC(55) is U_Hc_PsMc in component Q_Hc (mole_per_day).
C ALGBRC(83) is U_Ha_McMa in component Q_Ha (mole_per_day).
C ALGBRC(35) is C_Ha in component Q_Ha (molar).
C ALGBRC(48) is C_Hc in component Q_Hc (molar).
C CONSTS(26) is k_FlEx in component model_parameters (first_order_rate_constant).
C CONSTS(27) is J_Ha_McAm in component model_parameters (molar).
C STATES(5) is Q_Pu in component Q_Pu (mole).
C ALGBRC(2) is C_Pu in component Q_Pu (molar).
C CONSTS(158) is P_Pu_InPu in component Q_Pu (mole_per_day).
C ALGBRC(3) is U_Pu_PuEx in component Q_Pu (mole_per_day).
C CONSTS(28) is Y_Pu_InPu in component Q_Pu (molar).
C CONSTS(29) is D_Pu in component Q_Pu (litre_per_day).
C CONSTS(159) is k_PuEx in component Q_Pu (first_order_rate_constant).
C ALGBRC(19) is C_Am in component Q_Am (molar).
C ALGBRC(52) is U_Am_AmMc in component Q_Am (mole_per_day).
C ALGBRC(39) is U_Am_AmMa in component Q_Am (mole_per_day).
C ALGBRC(78) is U_Mc_McAm in component Q_Am (mole_per_day).
C STATES(6) is Q_Am in component Q_Am (mole).
C CONSTS(181) is P_Am_InAm in component Q_Am (mole_per_day).
C ALGBRC(23) is P_Am_UeAm in component Q_Am (mole_per_day).
C ALGBRC(38) is P_MaAm_PsAm in component Q_Am (mole_per_day).
C ALGBRC(51) is P_McAm_PsAm in component Q_Am (mole_per_day).
C ALGBRC(80) is P_Am_McAm in component Q_Am (mole_per_day).
C ALGBRC(20) is U_Am_AmAb in component Q_Am (mole_per_day).
C ALGBRC(21) is U_Am_AmEx in component Q_Am (mole_per_day).
C ALGBRC(22) is v_UeAm in component Q_Am (mole_per_day).
C CONSTS(182) is v_AmAb in component Q_Am (mole_per_day).
C ALGBRC(24) is v_AmMa in component Q_Am (mole_per_day).
C ALGBRC(25) is v_AmMc in component Q_Am (mole_per_day).
C CONSTS(30) is v_UeAm_star in component Q_Am (mole_per_litre_day).
C CONSTS(31) is v_AmAb_star in component Q_Am (mole_per_litre_day).
C CONSTS(32) is v_AmMa_star in component Q_Am (mole_per_gram_day).
C CONSTS(33) is v_AmMc_star in component Q_Am (mole_per_gram_day).
C CONSTS(160) is k_AmEx in component Q_Am (first_order_rate_constant).
C CONSTS(34) is Y_Am_InAm in component Q_Am (molar).
C CONSTS(35) is Y_Am_UeAm in component Q_Am (dimensionless).
C CONSTS(36) is Y_Am_PsAm in component Q_Am (dimensionless).
C CONSTS(37) is Y_Am_McAm in component Q_Am (dimensionless).
C CONSTS(38) is M_Am_AmMa in component Q_Am (molar).
C CONSTS(39) is M_Am_AmMc in component Q_Am (molar).
C CONSTS(40) is M_Am_AmAb in component Q_Am (molar).
C CONSTS(41) is M_pH_AmAb in component Q_Am (dimensionless).
C CONSTS(42) is phi_pH_AmAb in component Q_Am (dimensionless).
C CONSTS(43) is M_Ha_AmMa in component Q_Am (molar).
C CONSTS(44) is M_Hc_AmMc in component Q_Am (molar).
C CONSTS(45) is J_Am_UeAm in component Q_Am (molar).
C CONSTS(46) is D_Am in component Q_Am (litre_per_day).
C CONSTS(47) is pH in component model_parameters (dimensionless).
C CONSTS(48) is f_Po_Ma in component model_parameters (dimensionless).
C CONSTS(49) is f_Ni_Fe in component model_parameters (dimensionless).
C STATES(7) is Q_Li in component Q_Li (mole).
C ALGBRC(4) is C_Li in component Q_Li (molar).
C CONSTS(183) is P_Li_InLi in component Q_Li (mole_per_day).
C ALGBRC(64) is P_Li_MaMd in component Q_Li (mole_per_day).
C ALGBRC(26) is U_Li_LiEx in component Q_Li (mole_per_day).
C CONSTS(50) is Y_Li_InLi in component Q_Li (molar).
C CONSTS(51) is Y_Li_MaMd in component Q_Li (mole_per_gram).
C CONSTS(52) is D_Li in component Q_Li (litre_per_day).
C CONSTS(161) is k_LiEx in component Q_Li (first_order_rate_constant).
C ALGBRC(30) is U_Sd_SdHa in component Q_Sd (mole_per_day).
C STATES(8) is Q_Sd in component Q_Sd (mole).
C ALGBRC(27) is C_Sd in component Q_Sd (molar).
C CONSTS(184) is P_Sd_InSd in component Q_Sd (mole_per_day).
C ALGBRC(28) is U_Sd_SdEx in component Q_Sd (mole_per_day).
C CONSTS(53) is Y_Sd_InSd in component Q_Sd (molar).
C CONSTS(54) is D_Sd in component Q_Sd (litre_per_day).
C CONSTS(162) is k_SdEx in component Q_Sd (first_order_rate_constant).
C ALGBRC(29) is v_SdHa in component Q_Sd (mole_per_day).
C CONSTS(55) is v_SdHa_star in component Q_Sd (mole_per_gram_day).
C CONSTS(185) is M_Sd_SdHa in component Q_Sd (molar).
C CONSTS(56) is M_Sd_SdHa_star in component Q_Sd (molar).
C CONSTS(57) is T_Sd in component Q_Sd (day).
C CONSTS(58) is T_Sd_star in component Q_Sd (day).
C ALGBRC(34) is U_Fd_FdHc in component Q_Fd (mole_per_day).
C STATES(9) is Q_Fd in component Q_Fd (mole).
C ALGBRC(31) is C_Fd in component Q_Fd (molar).
C CONSTS(186) is P_Fd_InFd in component Q_Fd (mole_per_day).
C ALGBRC(32) is U_Fd_FdEx in component Q_Fd (mole_per_day).
C CONSTS(59) is Y_Fd_InFd in component Q_Fd (molar).
C CONSTS(60) is D_Fd in component Q_Fd (litre_per_day).
C CONSTS(163) is k_FdEx in component Q_Fd (first_order_rate_constant).
C ALGBRC(33) is v_FdHc in component Q_Fd (mole_per_day).
C CONSTS(61) is v_FdHc_star in component Q_Fd (mole_per_gram_day).
C CONSTS(187) is M_Fd_FdHc in component Q_Fd (molar).
C CONSTS(62) is M_Fd_FdHc_star in component Q_Fd (molar).
C CONSTS(63) is T_Fd in component Q_Fd (day).
C CONSTS(64) is T_Fd_star in component Q_Fd (day).
C CONSTS(65) is phi_pH_FdHc in component Q_Fd (dimensionless).
C CONSTS(66) is M_pH_FdHc in component Q_Fd (dimensionless).
C CONSTS(67) is pm in component model_parameters (dimensionless).
C CONSTS(68) is f in component model_parameters (dimensionless).
C STATES(10) is Q_Fu in component Q_Fu (mole).
C ALGBRC(5) is C_Fu in component Q_Fu (molar).
C CONSTS(164) is P_Fu_InFu in component Q_Fu (mole_per_day).
C ALGBRC(6) is U_Fu_FuEx in component Q_Fu (mole_per_day).
C CONSTS(69) is Y_Fu_InFu in component Q_Fu (molar).
C CONSTS(70) is D_Fu in component Q_Fu (litre_per_day).
C CONSTS(165) is k_FuEx in component Q_Fu (first_order_rate_constant).
C ALGBRC(45) is U_Ha_HaAs in component Q_Ha (mole_per_day).
C ALGBRC(41) is U_Ha_AmMa in component Q_Ha (mole_per_day).
C ALGBRC(42) is U_Ha_PsMa in component Q_Ha (mole_per_day).
C ALGBRC(47) is U_Ha_HaVa in component Q_Ha (mole_per_day).
C CONSTS(188) is P_Ha_WrHa in component Q_Ha (mole_per_day).
C CONSTS(190) is P_Ha_LaHa in component Q_Ha (mole_per_day).
C CONSTS(189) is P_Ha_SrHa in component Q_Ha (mole_per_day).
C ALGBRC(40) is P_Ha_SdHa in component Q_Ha (mole_per_day).
C ALGBRC(68) is P_Ha_MaMd in component Q_Ha (mole_per_day).
C STATES(11) is Q_Ha in component Q_Ha (mole).
C ALGBRC(43) is U_Ha_HaEx in component Q_Ha (mole_per_day).
C CONSTS(71) is R_Ha_AmMa in component Q_Ha (dimensionless).
C CONSTS(72) is R_Ha_PsMa in component Q_Ha (dimensionless).
C CONSTS(73) is R_Ha_McMa in component Q_Ha (mole_per_gram).
C CONSTS(74) is M_Ha_HaAs in component Q_Ha (molar).
C CONSTS(75) is M_Ha_HaVa in component Q_Ha (molar).
C ALGBRC(44) is v_HaAs in component Q_Ha (mole_per_day).
C ALGBRC(46) is v_HaVa in component Q_Ha (mole_per_day).
C CONSTS(76) is v_HaAs_star in component Q_Ha (mole_per_gram_day).
C CONSTS(77) is v_HaVa_star in component Q_Ha (mole_per_gram_day).
C CONSTS(166) is k_HaEx in component Q_Ha (first_order_rate_constant).
C CONSTS(78) is Y_Ha_WrHa in component Q_Ha (molar).
C CONSTS(79) is Y_Ha_SrHa in component Q_Ha (molar).
C CONSTS(80) is Y_Ha_LaHa in component Q_Ha (molar).
C CONSTS(81) is Y_Ha_SdHa in component Q_Ha (dimensionless).
C CONSTS(82) is Y_Ha_MaMd in component Q_Ha (mole_per_gram).
C CONSTS(83) is D_Wr in component Q_Ha (litre_per_day).
C CONSTS(84) is D_Sr in component Q_Ha (litre_per_day).
C CONSTS(85) is D_La in component Q_Ha (litre_per_day).
C CONSTS(86) is J_Am_HaVa in component Q_Ha (molar).
C CONSTS(87) is J_Ps_HaVa in component Q_Ha (molar).
C ALGBRC(67) is U_As_MaMd in component Q_As (gram_per_day).
C ALGBRC(81) is U_Mc_McMa in component Q_Ma (gram_per_day).
C ALGBRC(54) is U_Hc_AmMc in component Q_Hc (mole_per_day).
C ALGBRC(58) is U_Hc_HcVa in component Q_Hc (mole_per_day).
C STATES(12) is Q_Hc in component Q_Hc (mole).
C ALGBRC(53) is P_Hc_FdHc in component Q_Hc (mole_per_day).
C ALGBRC(56) is U_Hc_HcEx in component Q_Hc (mole_per_day).
C CONSTS(88) is R_Hc_AmMc in component Q_Hc (dimensionless).
C CONSTS(89) is R_Hc_PsMc in component Q_Hc (dimensionless).
C ALGBRC(57) is v_HcVa in component Q_Hc (mole_per_day).
C CONSTS(90) is v_HcVa_star in component Q_Hc (mole_per_gram_day).
C CONSTS(167) is k_HcEx in component Q_Hc (first_order_rate_constant).
C CONSTS(91) is Y_Hc_FdHc in component Q_Hc (dimensionless).
C CONSTS(92) is M_Hc_HcVa in component Q_Hc (molar).
C CONSTS(93) is J_Am_HcVa in component Q_Hc (molar).
C CONSTS(94) is J_Ps_HcVa in component Q_Hc (molar).
C ALGBRC(59) is C_Ma in component Q_Ma (gram_per_litre).
C ALGBRC(60) is P_Ma_AmMa in component Q_Ma (gram_per_day).
C ALGBRC(61) is P_Ma_PsMa in component Q_Ma (gram_per_day).
C ALGBRC(82) is P_Ma_McMa in component Q_Ma (gram_per_day).
C ALGBRC(65) is U_Ma_PoEx in component Q_Ma (gram_per_day).
C ALGBRC(7) is U_Ma_MaEx in component Q_Ma (gram_per_day).
C CONSTS(178) is k_MaMd in component Q_Ma (first_order_rate_constant).
C CONSTS(95) is Y_Ma_AmMa in component Q_Ma (gram_per_mole).
C CONSTS(96) is Y_Ma_PsMa in component Q_Ma (gram_per_mole).
C CONSTS(97) is Y_Ma_McMa in component Q_Ma (dimensionless).
C CONSTS(168) is k_MaEx in component model_parameters (first_order_rate_constant).
C CONSTS(176) is k_PoEx in component model_parameters (first_order_rate_constant).
C STATES(13) is Q_As in component Q_As (gram).
C ALGBRC(8) is C_As in component Q_As (gram_per_litre).
C ALGBRC(66) is P_As_HaAs in component Q_As (gram_per_day).
C ALGBRC(71) is U_As_PoEx in component Q_As (gram_per_day).
C ALGBRC(69) is U_As_AsEx in component Q_As (gram_per_day).
C CONSTS(169) is k_AsEx in component Q_As (first_order_rate_constant).
C CONSTS(179) is k_AsMd in component Q_As (first_order_rate_constant).
C CONSTS(98) is Y_As_HaAs in component Q_As (gram_per_mole).
C ALGBRC(70) is C_Mc in component Q_Mc (gram_per_litre).
C ALGBRC(72) is P_Mc_AmMc in component Q_Mc (gram_per_day).
C ALGBRC(73) is P_Mc_PsMc in component Q_Mc (gram_per_day).
C ALGBRC(74) is U_Mc_McEx in component Q_Mc (gram_per_day).
C CONSTS(170) is k_McEx in component Q_Mc (first_order_rate_constant).
C CONSTS(99) is Y_Mc_AmMc in component Q_Mc (gram_per_mole).
C CONSTS(100) is Y_Mc_PsMc in component Q_Mc (gram_per_mole).
C ALGBRC(75) is v_McEg in component Q_Mc (gram_per_day).
C CONSTS(101) is v_McEg_star in component Q_Mc (first_order_rate_constant).
C CONSTS(102) is M_Mc_McEg in component Q_Mc (gram_per_litre).
C STATES(14) is Q_Ac in component Q_Ac (mole).
C ALGBRC(84) is C_Ac in component Q_Ac (molar).
C CONSTS(191) is P_Ac_InAc in component Q_Ac (mole_per_day).
C ALGBRC(160) is P_Ac_AmMa in component Q_Ac (mole_per_day).
C ALGBRC(120) is P_Ac_AmMc in component Q_Ac (mole_per_day).
C ALGBRC(168) is P_Ac_PsMa in component Q_Ac (mole_per_day).
C ALGBRC(132) is P_Ac_PsMc in component Q_Ac (mole_per_day).
C ALGBRC(176) is P_Ac_HaAs in component Q_Ac (mole_per_day).
C ALGBRC(152) is P_Ac_HaVa in component Q_Ac (mole_per_day).
C ALGBRC(88) is P_Ac_HcVa in component Q_Ac (mole_per_day).
C ALGBRC(92) is P_MaAc_PsAm in component Q_Ac (mole_per_day).
C ALGBRC(96) is P_McAc_PsAm in component Q_Ac (mole_per_day).
C ALGBRC(144) is P_Ac_McMa in component Q_Ac (mole_per_day).
C ALGBRC(100) is P_Ac_McAm in component Q_Ac (mole_per_day).
C ALGBRC(104) is U_Ac_AcAb in component Q_Ac (mole_per_day).
C ALGBRC(108) is U_Ac_AcEx in component Q_Ac (mole_per_day).
C ALGBRC(112) is U_Hf_AmMa in component Q_Ac (mole_per_day).
C ALGBRC(116) is U_Hf_AmMc in component Q_Ac (mole_per_day).
C ALGBRC(124) is U_Hf_PsMa in component Q_Ac (mole_per_day).
C ALGBRC(128) is U_Hf_PsMc in component Q_Ac (mole_per_day).
C ALGBRC(136) is U_Hf_HaAs in component Q_Ac (mole_per_day).
C ALGBRC(140) is U_Hf_McMa in component Q_Ac (mole_per_day).
C CONSTS(103) is M_Ac_AcAb in component Q_Ac (molar).
C CONSTS(200) is v_AcAb in component Q_Ac (mole_per_day).
C CONSTS(104) is v_AcAb_star in component Q_Ac (mole_per_litre_day).
C CONSTS(192) is Y_Ac_LaAc in component Q_Ac (dimensionless).
C ALGBRC(148) is Y_Ac_HaVa in component Q_Ac (dimensionless).
C CONSTS(194) is Y_Ac_PsVa in component Q_Ac (dimensionless).
C ALGBRC(156) is Y_Ac_AmMa in component Q_Ac (dimensionless).
C CONSTS(195) is Y_Ac_AmMc in component Q_Ac (dimensionless).
C ALGBRC(164) is Y_Ac_PsMa in component Q_Ac (dimensionless).
C CONSTS(196) is Y_Ac_PsMc in component Q_Ac (dimensionless).
C ALGBRC(172) is Y_Ac_HaAs in component Q_Ac (dimensionless).
C CONSTS(197) is Y_Ac_McMa in component Q_Ac (dimensionless).
C CONSTS(198) is Y_Ac_PsAm in component Q_Ac (dimensionless).
C CONSTS(199) is Y_Ac_McAm in component Q_Ac (dimensionless).
C CONSTS(105) is Y_Ac_WrAc in component Q_Ac (dimensionless).
C CONSTS(106) is Y_Ac_InAc in component Q_Ac (molar).
C CONSTS(107) is Y_Ac_StAc in component Q_Ac (dimensionless).
C CONSTS(108) is Y_Ac_CeAc in component Q_Ac (dimensionless).
C CONSTS(109) is Y_Ac_HeAc in component Q_Ac (dimensionless).
C CONSTS(110) is Y_Ac_PsAc in component Q_Ac (dimensionless).
C CONSTS(193) is Y_Ac_HcVa in component Q_Ac (dimensionless).
C CONSTS(171) is k_AcEx in component Q_Ac (first_order_rate_constant).
C CONSTS(111) is D_Ac in component Q_Ac (litre_per_day).
C CONSTS(112) is J_pH_AcAb in component Q_Ac (dimensionless).
C CONSTS(113) is phi_pH_AcAb in component Q_Ac (dimensionless).
C CONSTS(114) is f_Hf_AmMa in component model_parameters (dimensionless).
C CONSTS(115) is f_Hf_AmMc in component model_parameters (dimensionless).
C CONSTS(116) is f_Hf_PsMa in component model_parameters (dimensionless).
C CONSTS(117) is f_Hf_PsMc in component model_parameters (dimensionless).
C CONSTS(118) is f_Hf_HaAs in component model_parameters (dimensionless).
C CONSTS(119) is f_Hf_McMa in component model_parameters (dimensionless).
C CONSTS(120) is f_Lc_Le in component model_parameters (dimensionless).
C CONSTS(121) is f_Ce_Fd in component model_parameters (dimensionless).
C STATES(15) is Q_Pr in component Q_Pr (mole).
C ALGBRC(85) is C_Pr in component Q_Pr (molar).
C CONSTS(201) is P_Pr_InPr in component Q_Pr (mole_per_day).
C ALGBRC(161) is P_Pr_AmMa in component Q_Pr (mole_per_day).
C ALGBRC(121) is P_Pr_AmMc in component Q_Pr (mole_per_day).
C ALGBRC(169) is P_Pr_PsMa in component Q_Pr (mole_per_day).
C ALGBRC(133) is P_Pr_PsMc in component Q_Pr (mole_per_day).
C ALGBRC(177) is P_Pr_HaAs in component Q_Pr (mole_per_day).
C ALGBRC(153) is P_Pr_HaVa in component Q_Pr (mole_per_day).
C ALGBRC(89) is P_Pr_HcVa in component Q_Pr (mole_per_day).
C ALGBRC(93) is P_MaPr_PsAm in component Q_Pr (mole_per_day).
C ALGBRC(97) is P_McPr_PsAm in component Q_Pr (mole_per_day).
C ALGBRC(145) is P_Pr_McMa in component Q_Pr (mole_per_day).
C ALGBRC(101) is P_Pr_McAm in component Q_Pr (mole_per_day).
C ALGBRC(105) is U_Pr_PrAb in component Q_Pr (mole_per_day).
C ALGBRC(109) is U_Pr_PrEx in component Q_Pr (mole_per_day).
C ALGBRC(113) is U_Hf_AmMa in component Q_Pr (mole_per_day).
C ALGBRC(117) is U_Hf_AmMc in component Q_Pr (mole_per_day).
C ALGBRC(125) is U_Hf_PsMa in component Q_Pr (mole_per_day).
C ALGBRC(129) is U_Hf_PsMc in component Q_Pr (mole_per_day).
C ALGBRC(137) is U_Hf_HaAs in component Q_Pr (mole_per_day).
C ALGBRC(141) is U_Hf_McMa in component Q_Pr (mole_per_day).
C CONSTS(122) is M_Pr_PrAb in component Q_Pr (molar).
C CONSTS(210) is v_PrAb in component Q_Pr (mole_per_day).
C CONSTS(123) is v_PrAb_star in component Q_Pr (mole_per_litre_day).
C CONSTS(202) is Y_Pr_LaPr in component Q_Pr (dimensionless).
C ALGBRC(149) is Y_Pr_HaVa in component Q_Pr (dimensionless).
C CONSTS(204) is Y_Pr_PsVa in component Q_Pr (dimensionless).
C ALGBRC(157) is Y_Pr_AmMa in component Q_Pr (dimensionless).
C CONSTS(205) is Y_Pr_AmMc in component Q_Pr (dimensionless).
C ALGBRC(165) is Y_Pr_PsMa in component Q_Pr (dimensionless).
C CONSTS(206) is Y_Pr_PsMc in component Q_Pr (dimensionless).
C ALGBRC(173) is Y_Pr_HaAs in component Q_Pr (dimensionless).
C CONSTS(207) is Y_Pr_McMa in component Q_Pr (dimensionless).
C CONSTS(208) is Y_Pr_PsAm in component Q_Pr (dimensionless).
C CONSTS(209) is Y_Pr_McAm in component Q_Pr (dimensionless).
C CONSTS(124) is Y_Pr_WrPr in component Q_Pr (dimensionless).
C CONSTS(125) is Y_Pr_InPr in component Q_Pr (molar).
C CONSTS(203) is Y_Pr_HcVa in component Q_Pr (dimensionless).
C CONSTS(126) is Y_Pr_StPr in component Q_Pr (dimensionless).
C CONSTS(127) is Y_Pr_CePr in component Q_Pr (dimensionless).
C CONSTS(128) is Y_Pr_HePr in component Q_Pr (dimensionless).
C CONSTS(129) is Y_Pr_PsPr in component Q_Pr (dimensionless).
C CONSTS(172) is k_PrEx in component Q_Pr (first_order_rate_constant).
C CONSTS(130) is D_Pr in component Q_Pr (litre_per_day).
C CONSTS(131) is phi_pH_PrAb in component Q_Pr (dimensionless).
C CONSTS(132) is J_pH_PrAb in component Q_Pr (dimensionless).
C STATES(16) is Q_Bu in component Q_Bu (mole).
C ALGBRC(86) is C_Bu in component Q_Bu (molar).
C CONSTS(211) is P_Bu_InBu in component Q_Bu (mole_per_day).
C ALGBRC(162) is P_Bu_AmMa in component Q_Bu (mole_per_day).
C ALGBRC(122) is P_Bu_AmMc in component Q_Bu (mole_per_day).
C ALGBRC(170) is P_Bu_PsMa in component Q_Bu (mole_per_day).
C ALGBRC(134) is P_Bu_PsMc in component Q_Bu (mole_per_day).
C ALGBRC(178) is P_Bu_HaAs in component Q_Bu (mole_per_day).
C ALGBRC(154) is P_Bu_HaVa in component Q_Bu (mole_per_day).
C ALGBRC(90) is P_Bu_HcVa in component Q_Bu (mole_per_day).
C ALGBRC(94) is P_MaBu_PsAm in component Q_Bu (mole_per_day).
C ALGBRC(98) is P_McBu_PsAm in component Q_Bu (mole_per_day).
C ALGBRC(146) is P_Bu_McMa in component Q_Bu (mole_per_day).
C ALGBRC(102) is P_Bu_McAm in component Q_Bu (mole_per_day).
C ALGBRC(106) is U_Bu_BuAb in component Q_Bu (mole_per_day).
C ALGBRC(110) is U_Bu_BuEx in component Q_Bu (mole_per_day).
C ALGBRC(114) is U_Hf_AmMa in component Q_Bu (mole_per_day).
C ALGBRC(118) is U_Hf_AmMc in component Q_Bu (mole_per_day).
C ALGBRC(126) is U_Hf_PsMa in component Q_Bu (mole_per_day).
C ALGBRC(130) is U_Hf_PsMc in component Q_Bu (mole_per_day).
C ALGBRC(138) is U_Hf_HaAs in component Q_Bu (mole_per_day).
C ALGBRC(142) is U_Hf_McMa in component Q_Bu (mole_per_day).
C CONSTS(133) is M_Bu_BuAb in component Q_Bu (molar).
C CONSTS(220) is v_BuAb in component Q_Bu (mole_per_day).
C CONSTS(134) is v_BuAb_star in component Q_Bu (mole_per_litre_day).
C CONSTS(212) is Y_Bu_LaBu in component Q_Bu (dimensionless).
C ALGBRC(150) is Y_Bu_HaVa in component Q_Bu (dimensionless).
C CONSTS(214) is Y_Bu_PsVa in component Q_Bu (dimensionless).
C ALGBRC(158) is Y_Bu_AmMa in component Q_Bu (dimensionless).
C CONSTS(215) is Y_Bu_AmMc in component Q_Bu (dimensionless).
C ALGBRC(166) is Y_Bu_PsMa in component Q_Bu (dimensionless).
C CONSTS(216) is Y_Bu_PsMc in component Q_Bu (dimensionless).
C ALGBRC(174) is Y_Bu_HaAs in component Q_Bu (dimensionless).
C CONSTS(217) is Y_Bu_McMa in component Q_Bu (dimensionless).
C CONSTS(218) is Y_Bu_PsAm in component Q_Bu (dimensionless).
C CONSTS(219) is Y_Bu_McAm in component Q_Bu (dimensionless).
C CONSTS(135) is Y_Bu_WrBu in component Q_Bu (dimensionless).
C CONSTS(136) is Y_Bu_InBu in component Q_Bu (molar).
C CONSTS(213) is Y_Bu_HcVa in component Q_Bu (dimensionless).
C CONSTS(137) is Y_Bu_StBu in component Q_Bu (dimensionless).
C CONSTS(138) is Y_Bu_CeBu in component Q_Bu (dimensionless).
C CONSTS(139) is Y_Bu_HeBu in component Q_Bu (dimensionless).
C CONSTS(140) is Y_Bu_PsBu in component Q_Bu (dimensionless).
C CONSTS(173) is k_BuEx in component Q_Bu (first_order_rate_constant).
C CONSTS(141) is D_Bu in component Q_Bu (litre_per_day).
C CONSTS(142) is phi_pH_BuAb in component Q_Bu (dimensionless).
C CONSTS(143) is J_pH_BuAb in component Q_Bu (dimensionless).
C STATES(17) is Q_Vl in component Q_Vl (mole).
C ALGBRC(87) is C_Vl in component Q_Vl (molar).
C CONSTS(221) is P_Vl_InVl in component Q_Vl (mole_per_day).
C ALGBRC(163) is P_Vl_AmMa in component Q_Vl (mole_per_day).
C ALGBRC(123) is P_Vl_AmMc in component Q_Vl (mole_per_day).
C ALGBRC(171) is P_Vl_PsMa in component Q_Vl (mole_per_day).
C ALGBRC(135) is P_Vl_PsMc in component Q_Vl (mole_per_day).
C ALGBRC(179) is P_Vl_HaAs in component Q_Vl (mole_per_day).
C ALGBRC(155) is P_Vl_HaVa in component Q_Vl (mole_per_day).
C ALGBRC(91) is P_Vl_HcVa in component Q_Vl (mole_per_day).
C ALGBRC(95) is P_MaVl_PsAm in component Q_Vl (mole_per_day).
C ALGBRC(99) is P_McVl_PsAm in component Q_Vl (mole_per_day).
C ALGBRC(147) is P_Vl_McMa in component Q_Vl (mole_per_day).
C ALGBRC(103) is P_Vl_McAm in component Q_Vl (mole_per_day).
C ALGBRC(107) is U_Vl_VlAb in component Q_Vl (mole_per_day).
C ALGBRC(111) is U_Vl_VlEx in component Q_Vl (mole_per_day).
C ALGBRC(115) is U_Hf_AmMa in component Q_Vl (mole_per_day).
C ALGBRC(119) is U_Hf_AmMc in component Q_Vl (mole_per_day).
C ALGBRC(127) is U_Hf_PsMa in component Q_Vl (mole_per_day).
C ALGBRC(131) is U_Hf_PsMc in component Q_Vl (mole_per_day).
C ALGBRC(139) is U_Hf_HaAs in component Q_Vl (mole_per_day).
C ALGBRC(143) is U_Hf_McMa in component Q_Vl (mole_per_day).
C CONSTS(144) is M_Vl_VlAb in component Q_Vl (molar).
C CONSTS(230) is v_VlAb in component Q_Vl (mole_per_day).
C CONSTS(145) is v_VlAb_star in component Q_Vl (mole_per_litre_day).
C CONSTS(222) is Y_Vl_LaVl in component Q_Vl (dimensionless).
C ALGBRC(151) is Y_Vl_HaVa in component Q_Vl (dimensionless).
C CONSTS(224) is Y_Vl_PsVa in component Q_Vl (dimensionless).
C ALGBRC(159) is Y_Vl_AmMa in component Q_Vl (dimensionless).
C CONSTS(225) is Y_Vl_AmMc in component Q_Vl (dimensionless).
C ALGBRC(167) is Y_Vl_PsMa in component Q_Vl (dimensionless).
C CONSTS(226) is Y_Vl_PsMc in component Q_Vl (dimensionless).
C ALGBRC(175) is Y_Vl_HaAs in component Q_Vl (dimensionless).
C CONSTS(227) is Y_Vl_McMa in component Q_Vl (dimensionless).
C CONSTS(228) is Y_Vl_PsAm in component Q_Vl (dimensionless).
C CONSTS(229) is Y_Vl_McAm in component Q_Vl (dimensionless).
C CONSTS(146) is Y_Vl_WrVl in component Q_Vl (dimensionless).
C CONSTS(147) is Y_Vl_InVl in component Q_Vl (molar).
C CONSTS(223) is Y_Vl_HcVa in component Q_Vl (dimensionless).
C CONSTS(148) is Y_Vl_StVl in component Q_Vl (dimensionless).
C CONSTS(149) is Y_Vl_CeVl in component Q_Vl (dimensionless).
C CONSTS(150) is Y_Vl_HeVl in component Q_Vl (dimensionless).
C CONSTS(151) is Y_Vl_PsVl in component Q_Vl (dimensionless).
C CONSTS(174) is k_VlEx in component Q_Vl (first_order_rate_constant).
C CONSTS(152) is D_Vl in component Q_Vl (litre_per_day).
C CONSTS(153) is phi_pH_VlAb in component Q_Vl (dimensionless).
C CONSTS(154) is J_pH_VlAb in component Q_Vl (dimensionless).
C RATES(1) is d/dt Q_Pd in component Q_Pd (mole).
C RATES(4) is d/dt Q_Ps in component Q_Ps (mole).
C RATES(5) is d/dt Q_Pu in component Q_Pu (mole).
C RATES(6) is d/dt Q_Am in component Q_Am (mole).
C RATES(7) is d/dt Q_Li in component Q_Li (mole).
C RATES(8) is d/dt Q_Sd in component Q_Sd (mole).
C RATES(9) is d/dt Q_Fd in component Q_Fd (mole).
C RATES(10) is d/dt Q_Fu in component Q_Fu (mole).
C RATES(11) is d/dt Q_Ha in component Q_Ha (mole).
C RATES(12) is d/dt Q_Hc in component Q_Hc (mole).
C RATES(2) is d/dt Q_Ma in component Q_Ma (gram).
C RATES(13) is d/dt Q_As in component Q_As (gram).
C RATES(3) is d/dt Q_Mc in component Q_Mc (gram).
C RATES(14) is d/dt Q_Ac in component Q_Ac (mole).
C RATES(15) is d/dt Q_Pr in component Q_Pr (mole).
C RATES(16) is d/dt Q_Bu in component Q_Bu (mole).
C RATES(17) is d/dt Q_Vl in component Q_Vl (mole).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 1.000
CONSTS(1) = 0.0576
CONSTS(2) = 0.264
CONSTS(3) = 0.0091
CONSTS(4) = 1.000
CONSTS(5) = 1.000
CONSTS(6) = 0.66
STATES(2) = 1.000
STATES(3) = 1.000
CONSTS(7) = 1.000
CONSTS(8) = 1.000
CONSTS(9) = 0.0248
STATES(4) = 1.000
CONSTS(10) = 0.0144
CONSTS(11) = 0.0576
CONSTS(12) = 0.0576
CONSTS(13) = 0.0091
CONSTS(14) = 0.0010
CONSTS(15) = 1.0
CONSTS(16) = 0.0067
CONSTS(17) = 0.0067
CONSTS(18) = 1.000
CONSTS(19) = 0.0289
CONSTS(20) = 0.0224
CONSTS(21) = 0.0224
CONSTS(22) = 0.0248
CONSTS(23) = 0.0248
CONSTS(24) = 0.0165
CONSTS(25) = 0.0165
CONSTS(26) = 1.000
CONSTS(27) = 0.0165
STATES(5) = 1.000
CONSTS(28) = 0.0091
CONSTS(29) = 1.000
STATES(6) = 1.000
CONSTS(30) = 0.00165
CONSTS(31) = 1.10
CONSTS(32) = 0.0528
CONSTS(33) = 0.0528
CONSTS(34) = 0.0588
CONSTS(35) = 2.0
CONSTS(36) = 1.257
CONSTS(37) = 0.0084
CONSTS(38) = 0.00135
CONSTS(39) = 0.00135
CONSTS(40) = 0.0132
CONSTS(41) = 7.5
CONSTS(42) = 7.85
CONSTS(43) = 0.0159
CONSTS(44) = 0.0159
CONSTS(45) = 0.00621
CONSTS(46) = 1.000
CONSTS(47) = 7.000
CONSTS(48) = 0.35
CONSTS(49) = 1.000
STATES(7) = 1.000
CONSTS(50) = 0.0015
CONSTS(51) = 0.00021
CONSTS(52) = 1.000
STATES(8) = 1.000
CONSTS(53) = 0.0062
CONSTS(54) = 1.000
CONSTS(55) = 0.2179
CONSTS(56) = 0.416
CONSTS(57) = 1.000
CONSTS(58) = 1.04
STATES(9) = 1.000
CONSTS(59) = 0.0062
CONSTS(60) = 1.000
CONSTS(61) = 0.1646
CONSTS(62) = 0.332
CONSTS(63) = 1.000
CONSTS(64) = 0.83
CONSTS(65) = 22.9
CONSTS(66) = 5.97
CONSTS(67) = 1.000
CONSTS(68) = 1.000
STATES(10) = 1.000
CONSTS(69) = 0.0062
CONSTS(70) = 1.000
STATES(11) = 1.000
CONSTS(71) = 1.793
CONSTS(72) = 1.291
CONSTS(73) = 0.0086
CONSTS(74) = 0.0268
CONSTS(75) = 0.055
CONSTS(76) = 0.053
CONSTS(77) = 0.1646
CONSTS(78) = 0.0062
CONSTS(79) = 0.0062
CONSTS(80) = 0.0025
CONSTS(81) = 1.0
CONSTS(82) = 0.0062
CONSTS(83) = 1.000
CONSTS(84) = 1.000
CONSTS(85) = 1.000
CONSTS(86) = 0.00861
CONSTS(87) = 0.01465
STATES(12) = 1.000
CONSTS(88) = 1.793
CONSTS(89) = 1.291
CONSTS(90) = 0.1646
CONSTS(91) = 1.0
CONSTS(92) = 0.055
CONSTS(93) = 0.00861
CONSTS(94) = 0.01465
CONSTS(95) = 118.91
CONSTS(96) = 149.48
CONSTS(97) = 149.48
STATES(13) = 1.000
CONSTS(98) = 112.5
CONSTS(99) = 118.91
CONSTS(100) = 149.48
CONSTS(101) = 15.439
CONSTS(102) = 34.694
STATES(14) = 1.000
CONSTS(103) = 0.338
CONSTS(104) = 7.86
CONSTS(105) = 1.000
CONSTS(106) = 0.0167
CONSTS(107) = 1.000
CONSTS(108) = 1.000
CONSTS(109) = 1.000
CONSTS(110) = 1.000
CONSTS(111) = 1.000
CONSTS(112) = 6.45
CONSTS(113) = 6.48
CONSTS(114) = 0.526
CONSTS(115) = 0.526
CONSTS(116) = 0.711
CONSTS(117) = 0.711
CONSTS(118) = 0.306
CONSTS(119) = 0.711
CONSTS(120) = 2.250
CONSTS(121) = 1.000
STATES(15) = 1.000
CONSTS(122) = 0.338
CONSTS(123) = 7.86
CONSTS(124) = 1.000
CONSTS(125) = 0.0135
CONSTS(126) = 1.000
CONSTS(127) = 1.000
CONSTS(128) = 1.000
CONSTS(129) = 1.000
CONSTS(130) = 1.000
CONSTS(131) = 6.48
CONSTS(132) = 6.45
STATES(16) = 1.000
CONSTS(133) = 0.338
CONSTS(134) = 7.86
CONSTS(135) = 1.000
CONSTS(136) = 0.0114
CONSTS(137) = 1.000
CONSTS(138) = 1.000
CONSTS(139) = 1.000
CONSTS(140) = 1.000
CONSTS(141) = 1.000
CONSTS(142) = 6.48
CONSTS(143) = 6.45
STATES(17) = 1.000
CONSTS(144) = 0.338
CONSTS(145) = 7.86
CONSTS(146) = 1.000
CONSTS(147) = 0.0098
CONSTS(148) = 1.000
CONSTS(149) = 1.000
CONSTS(150) = 1.000
CONSTS(151) = 1.000
CONSTS(152) = 1.000
CONSTS(153) = 6.48
CONSTS(154) = 6.45
CONSTS(155) = CONSTS(8)
CONSTS(156) = CONSTS(3)*CONSTS(4)
CONSTS(157) = CONSTS(26)
CONSTS(158) = CONSTS(28)*CONSTS(29)
CONSTS(159) = CONSTS(8)
CONSTS(160) = CONSTS(26)
CONSTS(161) = CONSTS(26)
CONSTS(162) = CONSTS(8)
CONSTS(163) = CONSTS(8)
CONSTS(164) = CONSTS(69)*CONSTS(70)
CONSTS(165) = CONSTS(8)
CONSTS(166) = CONSTS(26)
CONSTS(167) = CONSTS(26)
CONSTS(168) = CONSTS(26)
CONSTS(169) = CONSTS(26)
CONSTS(170) = CONSTS(8)
CONSTS(171) = CONSTS(26)
CONSTS(172) = CONSTS(26)
CONSTS(173) = CONSTS(26)
CONSTS(174) = CONSTS(26)
CONSTS(175) = CONSTS(2)*(CONSTS(5)/CONSTS(6))
CONSTS(176) = CONSTS(8)/2.00000
CONSTS(177) = CONSTS(13)*CONSTS(18)
CONSTS(178) = CONSTS(168) - CONSTS(176)
CONSTS(179) = CONSTS(168) - CONSTS(176)
CONSTS(180) = CONSTS(14)*CONSTS(18)
CONSTS(181) = CONSTS(34)*CONSTS(46)
CONSTS(182) = CONSTS(31)*1.00000 ** 0.250000*CONSTS(7) ** 0.750000*(1.00000+CONSTS(41)/CONSTS(47) ** CONSTS(42))
CONSTS(183) = CONSTS(50)*CONSTS(52)
CONSTS(184) = CONSTS(53)*CONSTS(54)
CONSTS(185) = CONSTS(56)*(CONSTS(57)/CONSTS(58))
CONSTS(186) = CONSTS(59)*CONSTS(60)
CONSTS(187) = CONSTS(62)*(CONSTS(63)/CONSTS(64))
CONSTS(188) = CONSTS(78)*CONSTS(83)
CONSTS(189) = CONSTS(79)*CONSTS(84)
CONSTS(190) = CONSTS(80)*CONSTS(85)
CONSTS(191) = CONSTS(106)*CONSTS(111)
CONSTS(192) = CONSTS(120)*CONSTS(105)
CONSTS(193) = CONSTS(121)*CONSTS(108)+ (1.00000 - CONSTS(121))*CONSTS(109)
CONSTS(194) = CONSTS(110)
CONSTS(195) = CONSTS(193)
CONSTS(196) = CONSTS(193)
CONSTS(197) = CONSTS(193)
CONSTS(198) = CONSTS(194)
CONSTS(199) = CONSTS(194)
CONSTS(200) = CONSTS(104)*1.00000 ** 0.250000*(CONSTS(7) ** 0.750000/(1.00000+CONSTS(47)/CONSTS(112) ** CONSTS(113)))
CONSTS(201) = CONSTS(125)*CONSTS(130)
CONSTS(202) = CONSTS(120)*CONSTS(124)
CONSTS(203) = CONSTS(121)*CONSTS(127)+ (1.00000 - CONSTS(121))*CONSTS(128)
CONSTS(204) = CONSTS(129)
CONSTS(205) = CONSTS(203)
CONSTS(206) = CONSTS(203)
CONSTS(207) = CONSTS(203)
CONSTS(208) = CONSTS(204)
CONSTS(209) = CONSTS(204)
CONSTS(210) = CONSTS(123)*1.00000 ** 0.250000*(CONSTS(7) ** 0.750000/(1.00000+CONSTS(47)/CONSTS(132) ** CONSTS(131)))
CONSTS(211) = CONSTS(136)*CONSTS(141)
CONSTS(212) = CONSTS(120)*CONSTS(135)
CONSTS(213) = CONSTS(121)*CONSTS(138)+ (1.00000 - CONSTS(121))*CONSTS(139)
CONSTS(214) = CONSTS(140)
CONSTS(215) = CONSTS(213)
CONSTS(216) = CONSTS(213)
CONSTS(217) = CONSTS(213)
CONSTS(218) = CONSTS(214)
CONSTS(219) = CONSTS(214)
CONSTS(220) = CONSTS(134)*1.00000 ** 0.250000*(CONSTS(7) ** 0.750000/(1.00000+CONSTS(47)/CONSTS(143) ** CONSTS(142)))
CONSTS(221) = CONSTS(147)*CONSTS(152)
CONSTS(222) = CONSTS(120)*CONSTS(146)
CONSTS(223) = CONSTS(121)*CONSTS(149)+ (1.00000 - CONSTS(121))*CONSTS(150)
CONSTS(224) = CONSTS(151)
CONSTS(225) = CONSTS(223)
CONSTS(226) = CONSTS(223)
CONSTS(227) = CONSTS(223)
CONSTS(228) = CONSTS(224)
CONSTS(229) = CONSTS(224)
CONSTS(230) = CONSTS(145)*1.00000 ** 0.250000*(CONSTS(7) ** 0.750000/(1.00000+CONSTS(47)/CONSTS(154) ** CONSTS(153)))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = CONSTS(159)*STATES(5)
RATES(5) = CONSTS(158) - ALGBRC(3)
ALGBRC(6) = CONSTS(165)*STATES(10)
RATES(10) = CONSTS(164) - ALGBRC(6)
ALGBRC(1) = STATES(1)/CONSTS(7)
ALGBRC(10) = CONSTS(1)*(STATES(2)+STATES(3))
ALGBRC(11) = ALGBRC(10)/(1.00000+CONSTS(175)/ALGBRC(1))
ALGBRC(9) = CONSTS(155)*STATES(1)
RATES(1) = CONSTS(156) - (ALGBRC(11)+ALGBRC(9))
ALGBRC(27) = STATES(8)/CONSTS(7)
ALGBRC(29) = CONSTS(55)*STATES(2)
ALGBRC(30) = ALGBRC(29)/(1.00000+CONSTS(185)/ALGBRC(27))
ALGBRC(28) = CONSTS(162)*STATES(8)
RATES(8) = CONSTS(184) - (ALGBRC(30)+ALGBRC(28))
ALGBRC(31) = STATES(9)/CONSTS(7)
ALGBRC(33) = CONSTS(61)*STATES(3)*((1.00000 - ( VOI*CONSTS(68))/1.00000)+(( VOI*CONSTS(68))/1.00000)/(1.00000+CONSTS(66)/CONSTS(67) ** CONSTS(65)))
ALGBRC(34) = ALGBRC(33)/(1.00000+CONSTS(187)/ALGBRC(31))
ALGBRC(32) = CONSTS(163)*STATES(9)
RATES(9) = CONSTS(186) - (ALGBRC(34)+ALGBRC(32))
ALGBRC(12) = STATES(4)/CONSTS(7)
ALGBRC(18) = CONSTS(12)*STATES(3)
ALGBRC(48) = STATES(12)/CONSTS(7)
ALGBRC(50) = ALGBRC(18)/(1.00000+CONSTS(21)/ALGBRC(12)+CONSTS(23)/ALGBRC(48))
ALGBRC(55) = CONSTS(89)*ALGBRC(50)
ALGBRC(19) = STATES(6)/CONSTS(7)
ALGBRC(25) = CONSTS(33)*STATES(3)
ALGBRC(52) = ALGBRC(25)/(1.00000+CONSTS(39)/ALGBRC(19)+CONSTS(44)/ALGBRC(48))
ALGBRC(54) = CONSTS(88)*ALGBRC(52)
ALGBRC(57) = CONSTS(90)*STATES(3)
ALGBRC(58) = ALGBRC(57)/( (1.00000+CONSTS(92)/ALGBRC(48))*(1.00000+ALGBRC(19)/CONSTS(93)+ALGBRC(12)/CONSTS(94)))
ALGBRC(53) = CONSTS(91)*ALGBRC(34)
ALGBRC(56) = CONSTS(167)*STATES(12)
RATES(12) = ALGBRC(53) - (ALGBRC(54)+ALGBRC(55)+ALGBRC(58)+ALGBRC(56))
ALGBRC(62) = CONSTS(178)*STATES(2)*CONSTS(48)
ALGBRC(64) = CONSTS(51)*ALGBRC(62)
ALGBRC(26) = CONSTS(161)*STATES(7)
RATES(7) = (CONSTS(183)+ALGBRC(64)) - ALGBRC(26)
ALGBRC(67) = CONSTS(179)*STATES(13)*CONSTS(48)
ALGBRC(35) = STATES(11)/CONSTS(7)
ALGBRC(44) = CONSTS(76)*STATES(2)
ALGBRC(45) = ALGBRC(44)/(1.00000+CONSTS(74)/ALGBRC(35))
ALGBRC(66) = CONSTS(98)*ALGBRC(45)
ALGBRC(71) = CONSTS(176)*STATES(13)*CONSTS(48)
ALGBRC(69) = CONSTS(169)*STATES(13)*(1.00000 - CONSTS(48))
RATES(13) = ALGBRC(66) - (ALGBRC(67)+ALGBRC(69)+ALGBRC(71))
ALGBRC(59) = STATES(2)/CONSTS(7)
ALGBRC(70) = STATES(3)/CONSTS(7)
ALGBRC(75) = ( CONSTS(101)*STATES(2)*CONSTS(48)*ALGBRC(70))/(ALGBRC(70)+ ALGBRC(59)*(1.00000 - CONSTS(48)))
ALGBRC(76) = ALGBRC(75)/(1.00000+CONSTS(102)/(ALGBRC(70)+ ALGBRC(59)*(1.00000 - CONSTS(48))))
ALGBRC(72) = CONSTS(99)*ALGBRC(52)
ALGBRC(73) = CONSTS(100)*ALGBRC(50)
ALGBRC(74) = CONSTS(170)*STATES(3)
RATES(3) = (ALGBRC(72)+ALGBRC(73)) - (ALGBRC(76)+ALGBRC(74))
ALGBRC(16) = CONSTS(10)*STATES(3)
ALGBRC(49) = ALGBRC(16)/(1.00000+CONSTS(19)/ALGBRC(12)+ALGBRC(48)/CONSTS(25))
ALGBRC(17) = CONSTS(11)*STATES(2)
ALGBRC(37) = ALGBRC(17)/(1.00000+CONSTS(20)/ALGBRC(12)+CONSTS(22)/ALGBRC(35))
ALGBRC(15) = CONSTS(10)*STATES(2)
ALGBRC(36) = ALGBRC(15)/(1.00000+CONSTS(19)/ALGBRC(12)+ALGBRC(35)/CONSTS(24))
ALGBRC(13) = CONSTS(15)*ALGBRC(11)
ALGBRC(63) = CONSTS(16)*ALGBRC(62)
ALGBRC(77) = ALGBRC(76)*(1.00000 - (1.00000/(1.00000+ALGBRC(35)/CONSTS(27))+1.00000/(1.00000+CONSTS(9)/ALGBRC(35))))
ALGBRC(79) = CONSTS(17)*ALGBRC(77)
ALGBRC(14) = CONSTS(157)*STATES(4)
RATES(4) = (CONSTS(177)+ALGBRC(13)+ALGBRC(63)+CONSTS(180)+ALGBRC(79)) - (ALGBRC(36)+ALGBRC(49)+ALGBRC(37)+ALGBRC(50)+ALGBRC(14))
ALGBRC(24) = CONSTS(32)*STATES(2)*(1.00000 - CONSTS(48))
ALGBRC(39) = ALGBRC(24)/(1.00000+CONSTS(38)/ALGBRC(19)+CONSTS(43)/ALGBRC(35))
ALGBRC(22) = CONSTS(30)*CONSTS(7)*(CONSTS(49)/(1.00000+ALGBRC(19)/CONSTS(45)))
ALGBRC(23) = CONSTS(35)*ALGBRC(22)
ALGBRC(38) = CONSTS(36)*ALGBRC(36)
ALGBRC(51) = CONSTS(36)*ALGBRC(49)
ALGBRC(78) = ( 1.00000*ALGBRC(76))/(1.00000+ALGBRC(35)/CONSTS(27))
ALGBRC(80) = CONSTS(37)*ALGBRC(78)
ALGBRC(20) = CONSTS(182)/(1.00000+CONSTS(40)/ALGBRC(19))
ALGBRC(21) = CONSTS(160)*STATES(6)
RATES(6) = (CONSTS(181)+ALGBRC(38)+ALGBRC(51)+ALGBRC(23)+ALGBRC(80)) - (ALGBRC(20)+ALGBRC(39)+ALGBRC(52)+ALGBRC(21))
ALGBRC(81) = ALGBRC(76)/(1.00000+CONSTS(9)/ALGBRC(35))
ALGBRC(83) = CONSTS(73)*ALGBRC(81)
ALGBRC(41) = CONSTS(71)*ALGBRC(39)
ALGBRC(42) = CONSTS(72)*ALGBRC(37)
ALGBRC(46) = CONSTS(77)*STATES(2)
ALGBRC(47) = ALGBRC(46)*((1.00000 - CONSTS(48))/( (1.00000+CONSTS(75)/ALGBRC(35))*(1.00000+ALGBRC(19)/CONSTS(86))*(ALGBRC(12)/CONSTS(87))))+( ALGBRC(46)*CONSTS(48))/(1.00000+CONSTS(75)/ALGBRC(35))
ALGBRC(40) = CONSTS(81)*ALGBRC(30)
ALGBRC(68) = CONSTS(82)*ALGBRC(67)
ALGBRC(43) = CONSTS(166)*STATES(11)
RATES(11) = (CONSTS(188)+CONSTS(189)+ALGBRC(40)+ALGBRC(68)+CONSTS(190)) - (ALGBRC(41)+ALGBRC(42)+ALGBRC(83)+ALGBRC(45)+ALGBRC(47)+ALGBRC(43))
ALGBRC(60) = CONSTS(95)*ALGBRC(39)
ALGBRC(61) = CONSTS(96)*ALGBRC(37)
ALGBRC(82) = CONSTS(97)*ALGBRC(81)
ALGBRC(65) = CONSTS(176)*STATES(2)*CONSTS(48)
RATES(2) = (ALGBRC(60)+ALGBRC(61)+ALGBRC(82)) - (ALGBRC(62)+ALGBRC(81)+ALGBRC(65))
ALGBRC(112) = CONSTS(114)*ALGBRC(41)
ALGBRC(148) = ( CONSTS(105)*CONSTS(188)+ CONSTS(192)*CONSTS(190)+ CONSTS(107)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(156) = ALGBRC(148)
ALGBRC(160) = ALGBRC(156)*ALGBRC(112)
ALGBRC(116) = CONSTS(115)*ALGBRC(54)
ALGBRC(120) = CONSTS(195)*ALGBRC(116)
ALGBRC(124) = CONSTS(116)*ALGBRC(42)
ALGBRC(164) = ALGBRC(148)
ALGBRC(168) = ALGBRC(164)*ALGBRC(124)
ALGBRC(128) = CONSTS(117)*ALGBRC(55)
ALGBRC(132) = CONSTS(196)*ALGBRC(128)
ALGBRC(136) = CONSTS(118)*ALGBRC(45)
ALGBRC(172) = ALGBRC(148)
ALGBRC(176) = ALGBRC(172)*ALGBRC(136)
ALGBRC(152) = ALGBRC(148)*ALGBRC(47)
ALGBRC(88) = CONSTS(193)*ALGBRC(58)
ALGBRC(92) = CONSTS(198)*ALGBRC(36)
ALGBRC(96) = CONSTS(198)*ALGBRC(49)
ALGBRC(140) = CONSTS(119)*ALGBRC(83)
ALGBRC(144) = CONSTS(197)*ALGBRC(140)
ALGBRC(100) = CONSTS(199)*ALGBRC(78)
ALGBRC(84) = STATES(14)/CONSTS(7)
ALGBRC(104) = CONSTS(200)/(1.00000+CONSTS(103)/ALGBRC(84))
ALGBRC(108) = CONSTS(171)*STATES(14)
RATES(14) = (CONSTS(191)+ALGBRC(160)+ALGBRC(120)+ALGBRC(168)+ALGBRC(132)+ALGBRC(176)+ALGBRC(152)+ALGBRC(88)+ALGBRC(92)+ALGBRC(96)+ALGBRC(144)+ALGBRC(100)) - (ALGBRC(104)+ALGBRC(108))
ALGBRC(113) = CONSTS(114)*ALGBRC(41)
ALGBRC(149) = ( CONSTS(124)*CONSTS(188)+ CONSTS(202)*CONSTS(190)+ CONSTS(126)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(157) = ALGBRC(149)
ALGBRC(161) = ALGBRC(157)*ALGBRC(113)
ALGBRC(117) = CONSTS(115)*ALGBRC(54)
ALGBRC(121) = CONSTS(205)*ALGBRC(117)
ALGBRC(125) = CONSTS(116)*ALGBRC(42)
ALGBRC(165) = ALGBRC(149)
ALGBRC(169) = ALGBRC(165)*ALGBRC(125)
ALGBRC(129) = CONSTS(117)*ALGBRC(55)
ALGBRC(133) = CONSTS(206)*ALGBRC(129)
ALGBRC(137) = CONSTS(118)*ALGBRC(45)
ALGBRC(173) = ALGBRC(149)
ALGBRC(177) = ALGBRC(173)*ALGBRC(137)
ALGBRC(153) = ALGBRC(149)*ALGBRC(47)
ALGBRC(89) = CONSTS(203)*ALGBRC(58)
ALGBRC(93) = CONSTS(208)*ALGBRC(36)
ALGBRC(97) = CONSTS(208)*ALGBRC(49)
ALGBRC(141) = CONSTS(119)*ALGBRC(83)
ALGBRC(145) = CONSTS(207)*ALGBRC(141)
ALGBRC(101) = CONSTS(209)*ALGBRC(78)
ALGBRC(85) = STATES(15)/CONSTS(7)
ALGBRC(105) = CONSTS(210)/(1.00000+CONSTS(122)/ALGBRC(85))
ALGBRC(109) = CONSTS(172)*STATES(15)
RATES(15) = (CONSTS(201)+ALGBRC(161)+ALGBRC(121)+ALGBRC(169)+ALGBRC(133)+ALGBRC(177)+ALGBRC(153)+ALGBRC(89)+ALGBRC(93)+ALGBRC(97)+ALGBRC(145)+ALGBRC(101)) - (ALGBRC(105)+ALGBRC(109))
ALGBRC(114) = CONSTS(114)*ALGBRC(41)
ALGBRC(150) = ( CONSTS(135)*CONSTS(188)+ CONSTS(212)*CONSTS(190)+ CONSTS(137)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(158) = ALGBRC(150)
ALGBRC(162) = ALGBRC(158)*ALGBRC(114)
ALGBRC(118) = CONSTS(115)*ALGBRC(54)
ALGBRC(122) = CONSTS(215)*ALGBRC(118)
ALGBRC(126) = CONSTS(116)*ALGBRC(42)
ALGBRC(166) = ALGBRC(150)
ALGBRC(170) = ALGBRC(166)*ALGBRC(126)
ALGBRC(130) = CONSTS(117)*ALGBRC(55)
ALGBRC(134) = CONSTS(216)*ALGBRC(130)
ALGBRC(138) = CONSTS(118)*ALGBRC(45)
ALGBRC(174) = ALGBRC(150)
ALGBRC(178) = ALGBRC(174)*ALGBRC(138)
ALGBRC(154) = ALGBRC(150)*ALGBRC(47)
ALGBRC(90) = CONSTS(213)*ALGBRC(58)
ALGBRC(94) = CONSTS(218)*ALGBRC(36)
ALGBRC(98) = CONSTS(218)*ALGBRC(49)
ALGBRC(142) = CONSTS(119)*ALGBRC(83)
ALGBRC(146) = CONSTS(217)*ALGBRC(142)
ALGBRC(102) = CONSTS(219)*ALGBRC(78)
ALGBRC(86) = STATES(16)/CONSTS(7)
ALGBRC(106) = CONSTS(220)/(1.00000+CONSTS(133)/ALGBRC(86))
ALGBRC(110) = CONSTS(173)*STATES(16)
RATES(16) = (CONSTS(211)+ALGBRC(162)+ALGBRC(122)+ALGBRC(170)+ALGBRC(134)+ALGBRC(178)+ALGBRC(154)+ALGBRC(90)+ALGBRC(94)+ALGBRC(98)+ALGBRC(146)+ALGBRC(102)) - (ALGBRC(106)+ALGBRC(110))
ALGBRC(115) = CONSTS(114)*ALGBRC(41)
ALGBRC(151) = ( CONSTS(146)*CONSTS(188)+ CONSTS(222)*CONSTS(190)+ CONSTS(148)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(159) = ALGBRC(151)
ALGBRC(163) = ALGBRC(159)*ALGBRC(115)
ALGBRC(119) = CONSTS(115)*ALGBRC(54)
ALGBRC(123) = CONSTS(225)*ALGBRC(119)
ALGBRC(127) = CONSTS(116)*ALGBRC(42)
ALGBRC(167) = ALGBRC(151)
ALGBRC(171) = ALGBRC(167)*ALGBRC(127)
ALGBRC(131) = CONSTS(117)*ALGBRC(55)
ALGBRC(135) = CONSTS(226)*ALGBRC(131)
ALGBRC(139) = CONSTS(118)*ALGBRC(45)
ALGBRC(175) = ALGBRC(151)
ALGBRC(179) = ALGBRC(175)*ALGBRC(139)
ALGBRC(155) = ALGBRC(151)*ALGBRC(47)
ALGBRC(91) = CONSTS(223)*ALGBRC(58)
ALGBRC(95) = CONSTS(228)*ALGBRC(36)
ALGBRC(99) = CONSTS(228)*ALGBRC(49)
ALGBRC(143) = CONSTS(119)*ALGBRC(83)
ALGBRC(147) = CONSTS(227)*ALGBRC(143)
ALGBRC(103) = CONSTS(229)*ALGBRC(78)
ALGBRC(87) = STATES(17)/CONSTS(7)
ALGBRC(107) = CONSTS(230)/(1.00000+CONSTS(144)/ALGBRC(87))
ALGBRC(111) = CONSTS(174)*STATES(17)
RATES(17) = (CONSTS(221)+ALGBRC(163)+ALGBRC(123)+ALGBRC(171)+ALGBRC(135)+ALGBRC(179)+ALGBRC(155)+ALGBRC(91)+ALGBRC(95)+ALGBRC(99)+ALGBRC(147)+ALGBRC(103)) - (ALGBRC(107)+ALGBRC(111))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = CONSTS(159)*STATES(5)
ALGBRC(6) = CONSTS(165)*STATES(10)
ALGBRC(1) = STATES(1)/CONSTS(7)
ALGBRC(10) = CONSTS(1)*(STATES(2)+STATES(3))
ALGBRC(11) = ALGBRC(10)/(1.00000+CONSTS(175)/ALGBRC(1))
ALGBRC(9) = CONSTS(155)*STATES(1)
ALGBRC(27) = STATES(8)/CONSTS(7)
ALGBRC(29) = CONSTS(55)*STATES(2)
ALGBRC(30) = ALGBRC(29)/(1.00000+CONSTS(185)/ALGBRC(27))
ALGBRC(28) = CONSTS(162)*STATES(8)
ALGBRC(31) = STATES(9)/CONSTS(7)
ALGBRC(33) = CONSTS(61)*STATES(3)*((1.00000 - ( VOI*CONSTS(68))/1.00000)+(( VOI*CONSTS(68))/1.00000)/(1.00000+CONSTS(66)/CONSTS(67) ** CONSTS(65)))
ALGBRC(34) = ALGBRC(33)/(1.00000+CONSTS(187)/ALGBRC(31))
ALGBRC(32) = CONSTS(163)*STATES(9)
ALGBRC(12) = STATES(4)/CONSTS(7)
ALGBRC(18) = CONSTS(12)*STATES(3)
ALGBRC(48) = STATES(12)/CONSTS(7)
ALGBRC(50) = ALGBRC(18)/(1.00000+CONSTS(21)/ALGBRC(12)+CONSTS(23)/ALGBRC(48))
ALGBRC(55) = CONSTS(89)*ALGBRC(50)
ALGBRC(19) = STATES(6)/CONSTS(7)
ALGBRC(25) = CONSTS(33)*STATES(3)
ALGBRC(52) = ALGBRC(25)/(1.00000+CONSTS(39)/ALGBRC(19)+CONSTS(44)/ALGBRC(48))
ALGBRC(54) = CONSTS(88)*ALGBRC(52)
ALGBRC(57) = CONSTS(90)*STATES(3)
ALGBRC(58) = ALGBRC(57)/( (1.00000+CONSTS(92)/ALGBRC(48))*(1.00000+ALGBRC(19)/CONSTS(93)+ALGBRC(12)/CONSTS(94)))
ALGBRC(53) = CONSTS(91)*ALGBRC(34)
ALGBRC(56) = CONSTS(167)*STATES(12)
ALGBRC(62) = CONSTS(178)*STATES(2)*CONSTS(48)
ALGBRC(64) = CONSTS(51)*ALGBRC(62)
ALGBRC(26) = CONSTS(161)*STATES(7)
ALGBRC(67) = CONSTS(179)*STATES(13)*CONSTS(48)
ALGBRC(35) = STATES(11)/CONSTS(7)
ALGBRC(44) = CONSTS(76)*STATES(2)
ALGBRC(45) = ALGBRC(44)/(1.00000+CONSTS(74)/ALGBRC(35))
ALGBRC(66) = CONSTS(98)*ALGBRC(45)
ALGBRC(71) = CONSTS(176)*STATES(13)*CONSTS(48)
ALGBRC(69) = CONSTS(169)*STATES(13)*(1.00000 - CONSTS(48))
ALGBRC(59) = STATES(2)/CONSTS(7)
ALGBRC(70) = STATES(3)/CONSTS(7)
ALGBRC(75) = ( CONSTS(101)*STATES(2)*CONSTS(48)*ALGBRC(70))/(ALGBRC(70)+ ALGBRC(59)*(1.00000 - CONSTS(48)))
ALGBRC(76) = ALGBRC(75)/(1.00000+CONSTS(102)/(ALGBRC(70)+ ALGBRC(59)*(1.00000 - CONSTS(48))))
ALGBRC(72) = CONSTS(99)*ALGBRC(52)
ALGBRC(73) = CONSTS(100)*ALGBRC(50)
ALGBRC(74) = CONSTS(170)*STATES(3)
ALGBRC(16) = CONSTS(10)*STATES(3)
ALGBRC(49) = ALGBRC(16)/(1.00000+CONSTS(19)/ALGBRC(12)+ALGBRC(48)/CONSTS(25))
ALGBRC(17) = CONSTS(11)*STATES(2)
ALGBRC(37) = ALGBRC(17)/(1.00000+CONSTS(20)/ALGBRC(12)+CONSTS(22)/ALGBRC(35))
ALGBRC(15) = CONSTS(10)*STATES(2)
ALGBRC(36) = ALGBRC(15)/(1.00000+CONSTS(19)/ALGBRC(12)+ALGBRC(35)/CONSTS(24))
ALGBRC(13) = CONSTS(15)*ALGBRC(11)
ALGBRC(63) = CONSTS(16)*ALGBRC(62)
ALGBRC(77) = ALGBRC(76)*(1.00000 - (1.00000/(1.00000+ALGBRC(35)/CONSTS(27))+1.00000/(1.00000+CONSTS(9)/ALGBRC(35))))
ALGBRC(79) = CONSTS(17)*ALGBRC(77)
ALGBRC(14) = CONSTS(157)*STATES(4)
ALGBRC(24) = CONSTS(32)*STATES(2)*(1.00000 - CONSTS(48))
ALGBRC(39) = ALGBRC(24)/(1.00000+CONSTS(38)/ALGBRC(19)+CONSTS(43)/ALGBRC(35))
ALGBRC(22) = CONSTS(30)*CONSTS(7)*(CONSTS(49)/(1.00000+ALGBRC(19)/CONSTS(45)))
ALGBRC(23) = CONSTS(35)*ALGBRC(22)
ALGBRC(38) = CONSTS(36)*ALGBRC(36)
ALGBRC(51) = CONSTS(36)*ALGBRC(49)
ALGBRC(78) = ( 1.00000*ALGBRC(76))/(1.00000+ALGBRC(35)/CONSTS(27))
ALGBRC(80) = CONSTS(37)*ALGBRC(78)
ALGBRC(20) = CONSTS(182)/(1.00000+CONSTS(40)/ALGBRC(19))
ALGBRC(21) = CONSTS(160)*STATES(6)
ALGBRC(81) = ALGBRC(76)/(1.00000+CONSTS(9)/ALGBRC(35))
ALGBRC(83) = CONSTS(73)*ALGBRC(81)
ALGBRC(41) = CONSTS(71)*ALGBRC(39)
ALGBRC(42) = CONSTS(72)*ALGBRC(37)
ALGBRC(46) = CONSTS(77)*STATES(2)
ALGBRC(47) = ALGBRC(46)*((1.00000 - CONSTS(48))/( (1.00000+CONSTS(75)/ALGBRC(35))*(1.00000+ALGBRC(19)/CONSTS(86))*(ALGBRC(12)/CONSTS(87))))+( ALGBRC(46)*CONSTS(48))/(1.00000+CONSTS(75)/ALGBRC(35))
ALGBRC(40) = CONSTS(81)*ALGBRC(30)
ALGBRC(68) = CONSTS(82)*ALGBRC(67)
ALGBRC(43) = CONSTS(166)*STATES(11)
ALGBRC(60) = CONSTS(95)*ALGBRC(39)
ALGBRC(61) = CONSTS(96)*ALGBRC(37)
ALGBRC(82) = CONSTS(97)*ALGBRC(81)
ALGBRC(65) = CONSTS(176)*STATES(2)*CONSTS(48)
ALGBRC(112) = CONSTS(114)*ALGBRC(41)
ALGBRC(148) = ( CONSTS(105)*CONSTS(188)+ CONSTS(192)*CONSTS(190)+ CONSTS(107)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(156) = ALGBRC(148)
ALGBRC(160) = ALGBRC(156)*ALGBRC(112)
ALGBRC(116) = CONSTS(115)*ALGBRC(54)
ALGBRC(120) = CONSTS(195)*ALGBRC(116)
ALGBRC(124) = CONSTS(116)*ALGBRC(42)
ALGBRC(164) = ALGBRC(148)
ALGBRC(168) = ALGBRC(164)*ALGBRC(124)
ALGBRC(128) = CONSTS(117)*ALGBRC(55)
ALGBRC(132) = CONSTS(196)*ALGBRC(128)
ALGBRC(136) = CONSTS(118)*ALGBRC(45)
ALGBRC(172) = ALGBRC(148)
ALGBRC(176) = ALGBRC(172)*ALGBRC(136)
ALGBRC(152) = ALGBRC(148)*ALGBRC(47)
ALGBRC(88) = CONSTS(193)*ALGBRC(58)
ALGBRC(92) = CONSTS(198)*ALGBRC(36)
ALGBRC(96) = CONSTS(198)*ALGBRC(49)
ALGBRC(140) = CONSTS(119)*ALGBRC(83)
ALGBRC(144) = CONSTS(197)*ALGBRC(140)
ALGBRC(100) = CONSTS(199)*ALGBRC(78)
ALGBRC(84) = STATES(14)/CONSTS(7)
ALGBRC(104) = CONSTS(200)/(1.00000+CONSTS(103)/ALGBRC(84))
ALGBRC(108) = CONSTS(171)*STATES(14)
ALGBRC(113) = CONSTS(114)*ALGBRC(41)
ALGBRC(149) = ( CONSTS(124)*CONSTS(188)+ CONSTS(202)*CONSTS(190)+ CONSTS(126)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(157) = ALGBRC(149)
ALGBRC(161) = ALGBRC(157)*ALGBRC(113)
ALGBRC(117) = CONSTS(115)*ALGBRC(54)
ALGBRC(121) = CONSTS(205)*ALGBRC(117)
ALGBRC(125) = CONSTS(116)*ALGBRC(42)
ALGBRC(165) = ALGBRC(149)
ALGBRC(169) = ALGBRC(165)*ALGBRC(125)
ALGBRC(129) = CONSTS(117)*ALGBRC(55)
ALGBRC(133) = CONSTS(206)*ALGBRC(129)
ALGBRC(137) = CONSTS(118)*ALGBRC(45)
ALGBRC(173) = ALGBRC(149)
ALGBRC(177) = ALGBRC(173)*ALGBRC(137)
ALGBRC(153) = ALGBRC(149)*ALGBRC(47)
ALGBRC(89) = CONSTS(203)*ALGBRC(58)
ALGBRC(93) = CONSTS(208)*ALGBRC(36)
ALGBRC(97) = CONSTS(208)*ALGBRC(49)
ALGBRC(141) = CONSTS(119)*ALGBRC(83)
ALGBRC(145) = CONSTS(207)*ALGBRC(141)
ALGBRC(101) = CONSTS(209)*ALGBRC(78)
ALGBRC(85) = STATES(15)/CONSTS(7)
ALGBRC(105) = CONSTS(210)/(1.00000+CONSTS(122)/ALGBRC(85))
ALGBRC(109) = CONSTS(172)*STATES(15)
ALGBRC(114) = CONSTS(114)*ALGBRC(41)
ALGBRC(150) = ( CONSTS(135)*CONSTS(188)+ CONSTS(212)*CONSTS(190)+ CONSTS(137)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(158) = ALGBRC(150)
ALGBRC(162) = ALGBRC(158)*ALGBRC(114)
ALGBRC(118) = CONSTS(115)*ALGBRC(54)
ALGBRC(122) = CONSTS(215)*ALGBRC(118)
ALGBRC(126) = CONSTS(116)*ALGBRC(42)
ALGBRC(166) = ALGBRC(150)
ALGBRC(170) = ALGBRC(166)*ALGBRC(126)
ALGBRC(130) = CONSTS(117)*ALGBRC(55)
ALGBRC(134) = CONSTS(216)*ALGBRC(130)
ALGBRC(138) = CONSTS(118)*ALGBRC(45)
ALGBRC(174) = ALGBRC(150)
ALGBRC(178) = ALGBRC(174)*ALGBRC(138)
ALGBRC(154) = ALGBRC(150)*ALGBRC(47)
ALGBRC(90) = CONSTS(213)*ALGBRC(58)
ALGBRC(94) = CONSTS(218)*ALGBRC(36)
ALGBRC(98) = CONSTS(218)*ALGBRC(49)
ALGBRC(142) = CONSTS(119)*ALGBRC(83)
ALGBRC(146) = CONSTS(217)*ALGBRC(142)
ALGBRC(102) = CONSTS(219)*ALGBRC(78)
ALGBRC(86) = STATES(16)/CONSTS(7)
ALGBRC(106) = CONSTS(220)/(1.00000+CONSTS(133)/ALGBRC(86))
ALGBRC(110) = CONSTS(173)*STATES(16)
ALGBRC(115) = CONSTS(114)*ALGBRC(41)
ALGBRC(151) = ( CONSTS(146)*CONSTS(188)+ CONSTS(222)*CONSTS(190)+ CONSTS(148)*(CONSTS(189)+ALGBRC(40)+ALGBRC(68)))/(CONSTS(188)+CONSTS(190)+CONSTS(189)+ALGBRC(40)+ALGBRC(68))
ALGBRC(159) = ALGBRC(151)
ALGBRC(163) = ALGBRC(159)*ALGBRC(115)
ALGBRC(119) = CONSTS(115)*ALGBRC(54)
ALGBRC(123) = CONSTS(225)*ALGBRC(119)
ALGBRC(127) = CONSTS(116)*ALGBRC(42)
ALGBRC(167) = ALGBRC(151)
ALGBRC(171) = ALGBRC(167)*ALGBRC(127)
ALGBRC(131) = CONSTS(117)*ALGBRC(55)
ALGBRC(135) = CONSTS(226)*ALGBRC(131)
ALGBRC(139) = CONSTS(118)*ALGBRC(45)
ALGBRC(175) = ALGBRC(151)
ALGBRC(179) = ALGBRC(175)*ALGBRC(139)
ALGBRC(155) = ALGBRC(151)*ALGBRC(47)
ALGBRC(91) = CONSTS(223)*ALGBRC(58)
ALGBRC(95) = CONSTS(228)*ALGBRC(36)
ALGBRC(99) = CONSTS(228)*ALGBRC(49)
ALGBRC(143) = CONSTS(119)*ALGBRC(83)
ALGBRC(147) = CONSTS(227)*ALGBRC(143)
ALGBRC(103) = CONSTS(229)*ALGBRC(78)
ALGBRC(87) = STATES(17)/CONSTS(7)
ALGBRC(107) = CONSTS(230)/(1.00000+CONSTS(144)/ALGBRC(87))
ALGBRC(111) = CONSTS(174)*STATES(17)
ALGBRC(2) = STATES(5)/CONSTS(7)
ALGBRC(4) = STATES(7)/CONSTS(7)
ALGBRC(5) = STATES(10)/CONSTS(7)
ALGBRC(7) = CONSTS(168)*STATES(2)*(1.00000 - CONSTS(48))
ALGBRC(8) = STATES(13)/CONSTS(7)
RETURN
END