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