C
C There are a total of 132 entries in the algebraic variable array.
C There are a total of 94 entries in each of the rate and state variable arrays.
C There are a total of 82 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is c1 in component c1 (molecules).
C ALGBRC(1) is v1 in component v1 (flux).
C STATES(2) is c2 in component c2 (molecules).
C ALGBRC(11) is v6 in component v6 (flux).
C CONSTS(81) is v13 in component v13 (flux).
C STATES(3) is c3 in component c3 (molecules).
C ALGBRC(8) is v2 in component v2 (flux).
C STATES(4) is c4 in component c4 (molecules).
C ALGBRC(9) is v3 in component v3 (flux).
C STATES(5) is c5 in component c5 (molecules).
C ALGBRC(12) is v7 in component v7 (flux).
C ALGBRC(13) is v8 in component v8 (flux).
C STATES(6) is c6 in component c6 (molecules).
C ALGBRC(15) is v10 in component v10 (flux).
C ALGBRC(16) is v60 in component v60 (flux).
C STATES(7) is c7 in component c7 (molecules).
C ALGBRC(10) is v4 in component v4 (flux).
C ALGBRC(125) is v5 in component v5 (flux).
C STATES(8) is c8 in component c8 (molecules).
C ALGBRC(19) is v12 in component v12 (flux).
C ALGBRC(20) is v14 in component v14 (flux).
C ALGBRC(21) is v62 in component v62 (flux).
C STATES(9) is c9 in component c9 (molecules).
C ALGBRC(22) is v15 in component v15 (flux).
C ALGBRC(126) is v107 in component v107 (flux).
C ALGBRC(127) is v110 in component v110 (flux).
C ALGBRC(128) is v113 in component v113 (flux).
C ALGBRC(129) is v116 in component v116 (flux).
C ALGBRC(130) is v119 in component v119 (flux).
C ALGBRC(131) is v122 in component v122 (flux).
C ALGBRC(132) is v125 in component v125 (flux).
C STATES(10) is c10 in component c10 (molecules).
C ALGBRC(17) is v11 in component v11 (flux).
C STATES(11) is c11 in component c11 (molecules).
C STATES(12) is c12 in component c12 (molecules).
C ALGBRC(103) is v106 in component v106 (flux).
C ALGBRC(108) is v109 in component v109 (flux).
C ALGBRC(112) is v112 in component v112 (flux).
C ALGBRC(116) is v115 in component v115 (flux).
C ALGBRC(118) is v118 in component v118 (flux).
C ALGBRC(120) is v121 in component v121 (flux).
C ALGBRC(122) is v124 in component v124 (flux).
C STATES(13) is c13 in component c13 (molecules).
C ALGBRC(18) is v61 in component v61 (flux).
C STATES(14) is c14 in component c14 (molecules).
C STATES(15) is c15 in component c15 (molecules).
C ALGBRC(23) is v16 in component v16 (flux).
C ALGBRC(29) is v22 in component v22 (flux).
C ALGBRC(39) is v32 in component v32 (flux).
C ALGBRC(41) is v34 in component v34 (flux).
C ALGBRC(44) is v37 in component v37 (flux).
C ALGBRC(46) is v39 in component v39 (flux).
C ALGBRC(96) is v102 in component v102 (flux).
C STATES(16) is c16 in component c16 (molecules).
C STATES(17) is c17 in component c17 (molecules).
C ALGBRC(67) is v63 in component v63 (flux).
C ALGBRC(73) is v69 in component v69 (flux).
C ALGBRC(87) is v79 in component v79 (flux).
C ALGBRC(89) is v80 in component v80 (flux).
C ALGBRC(91) is v81 in component v81 (flux).
C ALGBRC(93) is v82 in component v82 (flux).
C STATES(18) is c18 in component c18 (molecules).
C ALGBRC(14) is v9 in component v9 (flux).
C ALGBRC(68) is v64 in component v64 (flux).
C STATES(19) is c19 in component c19 (molecules).
C ALGBRC(69) is v65 in component v65 (flux).
C ALGBRC(70) is v66 in component v66 (flux).
C ALGBRC(71) is v67 in component v67 (flux).
C ALGBRC(72) is v68 in component v68 (flux).
C ALGBRC(99) is v105 in component v105 (flux).
C STATES(20) is c20 in component c20 (molecules).
C ALGBRC(106) is v108 in component v108 (flux).
C STATES(21) is c21 in component c21 (molecules).
C ALGBRC(110) is v111 in component v111 (flux).
C STATES(22) is c22 in component c22 (molecules).
C ALGBRC(31) is v24 in component v24 (flux).
C ALGBRC(42) is v35 in component v35 (flux).
C ALGBRC(45) is v38 in component v38 (flux).
C ALGBRC(75) is v71 in component v71 (flux).
C STATES(23) is c23 in component c23 (molecules).
C ALGBRC(24) is v17 in component v17 (flux).
C STATES(24) is c24 in component c24 (molecules).
C ALGBRC(77) is v72 in component v72 (flux).
C ALGBRC(32) is v25 in component v25 (flux).
C ALGBRC(47) is v40 in component v40 (flux).
C STATES(25) is c25 in component c25 (molecules).
C ALGBRC(25) is v18 in component v18 (flux).
C ALGBRC(26) is v19 in component v19 (flux).
C ALGBRC(27) is v20 in component v20 (flux).
C ALGBRC(28) is v21 in component v21 (flux).
C STATES(26) is c26 in component c26 (molecules).
C ALGBRC(33) is v26 in component v26 (flux).
C ALGBRC(38) is v31 in component v31 (flux).
C ALGBRC(84) is v78 in component v78 (flux).
C ALGBRC(78) is v73 in component v73 (flux).
C STATES(27) is c27 in component c27 (molecules).
C STATES(28) is c28 in component c28 (molecules).
C ALGBRC(34) is v27 in component v27 (flux).
C ALGBRC(35) is v28 in component v28 (flux).
C STATES(29) is c29 in component c29 (molecules).
C STATES(30) is c30 in component c30 (molecules).
C ALGBRC(40) is v33 in component v33 (flux).
C ALGBRC(48) is v41 in component v41 (flux).
C ALGBRC(95) is v83 in component v83 (flux).
C STATES(31) is c31 in component c31 (molecules).
C ALGBRC(43) is v36 in component v36 (flux).
C STATES(32) is c32 in component c32 (molecules).
C ALGBRC(30) is v23 in component v23 (flux).
C ALGBRC(76) is v103 in component v103 (flux).
C STATES(33) is c33 in component c33 (molecules).
C ALGBRC(100) is v104 in component v104 (flux).
C STATES(34) is c34 in component c34 (molecules).
C ALGBRC(114) is v114 in component v114 (flux).
C STATES(35) is c35 in component c35 (molecules).
C ALGBRC(37) is v30 in component v30 (flux).
C ALGBRC(117) is v117 in component v117 (flux).
C STATES(36) is c36 in component c36 (molecules).
C ALGBRC(119) is v120 in component v120 (flux).
C STATES(37) is c37 in component c37 (molecules).
C ALGBRC(121) is v123 in component v123 (flux).
C STATES(38) is c38 in component c38 (molecules).
C STATES(39) is c39 in component c39 (molecules).
C STATES(40) is c40 in component c40 (molecules).
C STATES(41) is c41 in component c41 (molecules).
C ALGBRC(50) is v43 in component v43 (flux).
C ALGBRC(80) is v75 in component v75 (flux).
C ALGBRC(85) is v85 in component v85 (flux).
C STATES(42) is c42 in component c42 (molecules).
C ALGBRC(36) is v29 in component v29 (flux).
C STATES(43) is c43 in component c43 (molecules).
C STATES(44) is c44 in component c44 (molecules).
C ALGBRC(49) is v42 in component v42 (flux).
C ALGBRC(83) is v84 in component v84 (flux).
C STATES(45) is c45 in component c45 (molecules).
C ALGBRC(51) is v44 in component v44 (flux).
C ALGBRC(52) is v45 in component v45 (flux).
C ALGBRC(53) is v46 in component v46 (flux).
C ALGBRC(54) is v47 in component v47 (flux).
C STATES(46) is c46 in component c46 (molecules).
C STATES(47) is c47 in component c47 (molecules).
C ALGBRC(58) is v51 in component v51 (flux).
C ALGBRC(86) is v86 in component v86 (flux).
C ALGBRC(104) is v93 in component v93 (flux).
C STATES(48) is c48 in component c48 (molecules).
C STATES(49) is c49 in component c49 (molecules).
C ALGBRC(56) is v49 in component v49 (flux).
C ALGBRC(57) is v50 in component v50 (flux).
C STATES(50) is c50 in component c50 (molecules).
C STATES(51) is c51 in component c51 (molecules).
C ALGBRC(55) is v48 in component v48 (flux).
C ALGBRC(59) is v52 in component v52 (flux).
C ALGBRC(60) is v53 in component v53 (flux).
C ALGBRC(61) is v54 in component v54 (flux).
C ALGBRC(62) is v55 in component v55 (flux).
C STATES(52) is c52 in component c52 (molecules).
C STATES(53) is c53 in component c53 (molecules).
C ALGBRC(94) is v90 in component v90 (flux).
C ALGBRC(97) is v91 in component v91 (flux).
C ALGBRC(101) is v92 in component v92 (flux).
C STATES(54) is c54 in component c54 (molecules).
C STATES(55) is c55 in component c55 (molecules).
C ALGBRC(66) is v59 in component v59 (flux).
C ALGBRC(98) is v94 in component v94 (flux).
C ALGBRC(115) is v101 in component v101 (flux).
C STATES(56) is c56 in component c56 (molecules).
C STATES(57) is c57 in component c57 (molecules).
C ALGBRC(64) is v57 in component v57 (flux).
C ALGBRC(65) is v58 in component v58 (flux).
C STATES(58) is c58 in component c58 (molecules).
C STATES(59) is c59 in component c59 (molecules).
C ALGBRC(63) is v56 in component v56 (flux).
C STATES(60) is c60 in component c60 (molecules).
C ALGBRC(109) is v98 in component v98 (flux).
C ALGBRC(111) is v99 in component v99 (flux).
C ALGBRC(113) is v100 in component v100 (flux).
C STATES(61) is c61 in component c61 (molecules).
C STATES(62) is c62 in component c62 (molecules).
C STATES(63) is c63 in component c63 (molecules).
C ALGBRC(74) is v70 in component v70 (flux).
C STATES(64) is c64 in component c64 (molecules).
C STATES(65) is c65 in component c65 (molecules).
C STATES(66) is c66 in component c66 (molecules).
C ALGBRC(79) is v74 in component v74 (flux).
C ALGBRC(82) is v77 in component v77 (flux).
C STATES(67) is c67 in component c67 (molecules).
C STATES(68) is c68 in component c68 (molecules).
C STATES(69) is c69 in component c69 (molecules).
C STATES(70) is c70 in component c70 (molecules).
C ALGBRC(81) is v76 in component v76 (flux).
C STATES(71) is c71 in component c71 (molecules).
C STATES(72) is c72 in component c72 (molecules).
C ALGBRC(88) is v87 in component v87 (flux).
C ALGBRC(90) is v88 in component v88 (flux).
C ALGBRC(92) is v89 in component v89 (flux).
C STATES(73) is c73 in component c73 (molecules).
C STATES(74) is c74 in component c74 (molecules).
C STATES(75) is c75 in component c75 (molecules).
C STATES(76) is c76 in component c76 (molecules).
C STATES(77) is c77 in component c77 (molecules).
C ALGBRC(102) is v95 in component v95 (flux).
C ALGBRC(105) is v96 in component v96 (flux).
C ALGBRC(107) is v97 in component v97 (flux).
C STATES(78) is c78 in component c78 (molecules).
C STATES(79) is c79 in component c79 (molecules).
C STATES(80) is c80 in component c80 (molecules).
C STATES(81) is c81 in component c81 (molecules).
C STATES(82) is c82 in component c82 (molecules).
C STATES(83) is c83 in component c83 (molecules).
C STATES(84) is c84 in component c84 (molecules).
C STATES(85) is c85 in component c85 (molecules).
C STATES(86) is c86 in component c86 (molecules).
C STATES(87) is c87 in component c87 (molecules).
C STATES(88) is c88 in component c88 (molecules).
C STATES(89) is c89 in component c89 (molecules).
C STATES(90) is c90 in component c90 (molecules).
C STATES(91) is c91 in component c91 (molecules).
C STATES(92) is c92 in component c92 (molecules).
C STATES(93) is c93 in component c93 (molecules).
C STATES(94) is c94 in component c94 (molecules).
C ALGBRC(2) is Raf_activation in component Raf_activation (molecules).
C ALGBRC(3) is MEKPP_total in component MEKPP_total (molecules).
C ALGBRC(4) is ERKPP_total in component ERKPP_total (molecules).
C ALGBRC(5) is Ras_GTP_total in component Ras_GTP_total (molecules).
C ALGBRC(6) is ShcP_total in component ShcP_total (molecules).
C ALGBRC(7) is EGF_EGFR_totalx2 in component EGF_EGFR_totalx2 (molecules).
C CONSTS(1) is k1 in component reaction_rates (second_order_rate_constant).
C CONSTS(2) is kd1 in component reaction_rates (first_order_rate_constant).
C CONSTS(3) is k2 in component reaction_rates (second_order_rate_constant).
C CONSTS(4) is kd2 in component reaction_rates (first_order_rate_constant).
C CONSTS(5) is k3 in component reaction_rates (first_order_rate_constant).
C CONSTS(6) is kd3 in component reaction_rates (first_order_rate_constant).
C CONSTS(7) is k4 in component reaction_rates (second_order_rate_constant).
C CONSTS(8) is kd4 in component reaction_rates (first_order_rate_constant).
C ALGBRC(124) is k5 in component reaction_rates (first_order_rate_constant).
C CONSTS(9) is k6 in component reaction_rates (first_order_rate_constant).
C CONSTS(10) is kd6 in component reaction_rates (first_order_rate_constant).
C CONSTS(11) is k8 in component reaction_rates (second_order_rate_constant).
C CONSTS(12) is kd8 in component reaction_rates (first_order_rate_constant).
C CONSTS(13) is k10b in component reaction_rates (second_order_rate_constant).
C CONSTS(14) is kd10 in component reaction_rates (first_order_rate_constant).
C CONSTS(15) is k13 in component reaction_rates (first_order_rate_constant).
C CONSTS(16) is k15 in component reaction_rates (first_order_rate_constant).
C CONSTS(17) is k16 in component reaction_rates (second_order_rate_constant).
C CONSTS(18) is kd16 in component reaction_rates (first_order_rate_constant).
C CONSTS(19) is k17 in component reaction_rates (second_order_rate_constant).
C CONSTS(20) is kd17 in component reaction_rates (first_order_rate_constant).
C CONSTS(21) is k18 in component reaction_rates (second_order_rate_constant).
C CONSTS(22) is kd18 in component reaction_rates (first_order_rate_constant).
C CONSTS(23) is k19 in component reaction_rates (first_order_rate_constant).
C CONSTS(24) is kd19 in component reaction_rates (second_order_rate_constant).
C CONSTS(25) is k20 in component reaction_rates (second_order_rate_constant).
C CONSTS(26) is kd20 in component reaction_rates (first_order_rate_constant).
C CONSTS(27) is k21 in component reaction_rates (first_order_rate_constant).
C CONSTS(28) is kd21 in component reaction_rates (second_order_rate_constant).
C CONSTS(29) is k22 in component reaction_rates (second_order_rate_constant).
C CONSTS(30) is kd22 in component reaction_rates (first_order_rate_constant).
C CONSTS(31) is k23 in component reaction_rates (first_order_rate_constant).
C CONSTS(32) is kd23 in component reaction_rates (first_order_rate_constant).
C CONSTS(33) is kd24 in component reaction_rates (first_order_rate_constant).
C CONSTS(34) is k25 in component reaction_rates (second_order_rate_constant).
C CONSTS(35) is kd25 in component reaction_rates (first_order_rate_constant).
C CONSTS(36) is k28 in component reaction_rates (second_order_rate_constant).
C CONSTS(37) is kd28 in component reaction_rates (first_order_rate_constant).
C CONSTS(38) is k29 in component reaction_rates (first_order_rate_constant).
C CONSTS(39) is kd29 in component reaction_rates (second_order_rate_constant).
C CONSTS(40) is k32 in component reaction_rates (first_order_rate_constant).
C CONSTS(41) is kd32 in component reaction_rates (second_order_rate_constant).
C CONSTS(42) is k33 in component reaction_rates (first_order_rate_constant).
C CONSTS(43) is kd33 in component reaction_rates (second_order_rate_constant).
C CONSTS(44) is k34 in component reaction_rates (first_order_rate_constant).
C CONSTS(45) is kd34 in component reaction_rates (second_order_rate_constant).
C CONSTS(46) is k35 in component reaction_rates (first_order_rate_constant).
C CONSTS(47) is kd35 in component reaction_rates (second_order_rate_constant).
C CONSTS(48) is Km36 in component reaction_rates (molecules).
C CONSTS(49) is Vm36 in component reaction_rates (flux).
C CONSTS(50) is k37 in component reaction_rates (first_order_rate_constant).
C CONSTS(51) is kd37 in component reaction_rates (second_order_rate_constant).
C CONSTS(52) is k40 in component reaction_rates (second_order_rate_constant).
C CONSTS(53) is kd40 in component reaction_rates (first_order_rate_constant).
C CONSTS(54) is k41 in component reaction_rates (second_order_rate_constant).
C CONSTS(55) is kd41 in component reaction_rates (first_order_rate_constant).
C CONSTS(56) is k42 in component reaction_rates (second_order_rate_constant).
C CONSTS(57) is kd42 in component reaction_rates (first_order_rate_constant).
C CONSTS(58) is k43 in component reaction_rates (first_order_rate_constant).
C CONSTS(59) is k44 in component reaction_rates (second_order_rate_constant).
C CONSTS(60) is kd44 in component reaction_rates (first_order_rate_constant).
C CONSTS(61) is k45 in component reaction_rates (first_order_rate_constant).
C CONSTS(62) is k47 in component reaction_rates (first_order_rate_constant).
C CONSTS(63) is k48 in component reaction_rates (second_order_rate_constant).
C CONSTS(64) is kd48 in component reaction_rates (first_order_rate_constant).
C CONSTS(65) is k49 in component reaction_rates (first_order_rate_constant).
C CONSTS(66) is k50 in component reaction_rates (second_order_rate_constant).
C CONSTS(67) is kd50 in component reaction_rates (first_order_rate_constant).
C CONSTS(68) is k52 in component reaction_rates (second_order_rate_constant).
C CONSTS(69) is kd52 in component reaction_rates (first_order_rate_constant).
C CONSTS(70) is k53 in component reaction_rates (first_order_rate_constant).
C CONSTS(71) is k55 in component reaction_rates (first_order_rate_constant).
C CONSTS(72) is k56 in component reaction_rates (second_order_rate_constant).
C CONSTS(73) is kd56 in component reaction_rates (first_order_rate_constant).
C CONSTS(74) is k57 in component reaction_rates (first_order_rate_constant).
C CONSTS(75) is k58 in component reaction_rates (second_order_rate_constant).
C CONSTS(76) is kd58 in component reaction_rates (first_order_rate_constant).
C CONSTS(77) is k59 in component reaction_rates (first_order_rate_constant).
C CONSTS(78) is k60 in component reaction_rates (first_order_rate_constant).
C CONSTS(79) is k61 in component reaction_rates (first_order_rate_constant).
C ALGBRC(123) is C in component C (flux).
C CONSTS(80) is RT in component C (molecules).
C RATES(1) is d/dt c1 in component c1 (molecules).
C RATES(2) is d/dt c2 in component c2 (molecules).
C RATES(3) is d/dt c3 in component c3 (molecules).
C RATES(4) is d/dt c4 in component c4 (molecules).
C RATES(5) is d/dt c5 in component c5 (molecules).
C RATES(6) is d/dt c6 in component c6 (molecules).
C RATES(7) is d/dt c7 in component c7 (molecules).
C RATES(8) is d/dt c8 in component c8 (molecules).
C RATES(9) is d/dt c9 in component c9 (molecules).
C RATES(10) is d/dt c10 in component c10 (molecules).
C RATES(11) is d/dt c11 in component c11 (molecules).
C RATES(12) is d/dt c12 in component c12 (molecules).
C RATES(13) is d/dt c13 in component c13 (molecules).
C RATES(14) is d/dt c14 in component c14 (molecules).
C RATES(15) is d/dt c15 in component c15 (molecules).
C RATES(16) is d/dt c16 in component c16 (molecules).
C RATES(17) is d/dt c17 in component c17 (molecules).
C RATES(18) is d/dt c18 in component c18 (molecules).
C RATES(19) is d/dt c19 in component c19 (molecules).
C RATES(20) is d/dt c20 in component c20 (molecules).
C RATES(21) is d/dt c21 in component c21 (molecules).
C RATES(22) is d/dt c22 in component c22 (molecules).
C RATES(23) is d/dt c23 in component c23 (molecules).
C RATES(24) is d/dt c24 in component c24 (molecules).
C RATES(25) is d/dt c25 in component c25 (molecules).
C RATES(26) is d/dt c26 in component c26 (molecules).
C RATES(27) is d/dt c27 in component c27 (molecules).
C RATES(28) is d/dt c28 in component c28 (molecules).
C RATES(29) is d/dt c29 in component c29 (molecules).
C RATES(30) is d/dt c30 in component c30 (molecules).
C RATES(31) is d/dt c31 in component c31 (molecules).
C RATES(32) is d/dt c32 in component c32 (molecules).
C RATES(33) is d/dt c33 in component c33 (molecules).
C RATES(34) is d/dt c34 in component c34 (molecules).
C RATES(35) is d/dt c35 in component c35 (molecules).
C RATES(36) is d/dt c36 in component c36 (molecules).
C RATES(37) is d/dt c37 in component c37 (molecules).
C RATES(38) is d/dt c38 in component c38 (molecules).
C RATES(39) is d/dt c39 in component c39 (molecules).
C RATES(40) is d/dt c40 in component c40 (molecules).
C RATES(41) is d/dt c41 in component c41 (molecules).
C RATES(42) is d/dt c42 in component c42 (molecules).
C RATES(43) is d/dt c43 in component c43 (molecules).
C RATES(44) is d/dt c44 in component c44 (molecules).
C RATES(45) is d/dt c45 in component c45 (molecules).
C RATES(46) is d/dt c46 in component c46 (molecules).
C RATES(47) is d/dt c47 in component c47 (molecules).
C RATES(48) is d/dt c48 in component c48 (molecules).
C RATES(49) is d/dt c49 in component c49 (molecules).
C RATES(50) is d/dt c50 in component c50 (molecules).
C RATES(51) is d/dt c51 in component c51 (molecules).
C RATES(52) is d/dt c52 in component c52 (molecules).
C RATES(53) is d/dt c53 in component c53 (molecules).
C RATES(54) is d/dt c54 in component c54 (molecules).
C RATES(55) is d/dt c55 in component c55 (molecules).
C RATES(56) is d/dt c56 in component c56 (molecules).
C RATES(57) is d/dt c57 in component c57 (molecules).
C RATES(58) is d/dt c58 in component c58 (molecules).
C RATES(59) is d/dt c59 in component c59 (molecules).
C RATES(60) is d/dt c60 in component c60 (molecules).
C RATES(61) is d/dt c61 in component c61 (molecules).
C RATES(62) is d/dt c62 in component c62 (molecules).
C RATES(63) is d/dt c63 in component c63 (molecules).
C RATES(64) is d/dt c64 in component c64 (molecules).
C RATES(65) is d/dt c65 in component c65 (molecules).
C RATES(66) is d/dt c66 in component c66 (molecules).
C RATES(67) is d/dt c67 in component c67 (molecules).
C RATES(68) is d/dt c68 in component c68 (molecules).
C RATES(69) is d/dt c69 in component c69 (molecules).
C RATES(70) is d/dt c70 in component c70 (molecules).
C RATES(71) is d/dt c71 in component c71 (molecules).
C RATES(72) is d/dt c72 in component c72 (molecules).
C RATES(73) is d/dt c73 in component c73 (molecules).
C RATES(74) is d/dt c74 in component c74 (molecules).
C RATES(75) is d/dt c75 in component c75 (molecules).
C RATES(76) is d/dt c76 in component c76 (molecules).
C RATES(77) is d/dt c77 in component c77 (molecules).
C RATES(78) is d/dt c78 in component c78 (molecules).
C RATES(79) is d/dt c79 in component c79 (molecules).
C RATES(80) is d/dt c80 in component c80 (molecules).
C RATES(81) is d/dt c81 in component c81 (molecules).
C RATES(82) is d/dt c82 in component c82 (molecules).
C RATES(83) is d/dt c83 in component c83 (molecules).
C RATES(84) is d/dt c84 in component c84 (molecules).
C RATES(85) is d/dt c85 in component c85 (molecules).
C RATES(86) is d/dt c86 in component c86 (molecules).
C RATES(87) is d/dt c87 in component c87 (molecules).
C RATES(88) is d/dt c88 in component c88 (molecules).
C RATES(89) is d/dt c89 in component c89 (molecules).
C RATES(90) is d/dt c90 in component c90 (molecules).
C RATES(91) is d/dt c91 in component c91 (molecules).
C RATES(92) is d/dt c92 in component c92 (molecules).
C RATES(93) is d/dt c93 in component c93 (molecules).
C RATES(94) is d/dt c94 in component c94 (molecules).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 4962
      STATES(2) = 5e4
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 0
      STATES(11) = 0
      STATES(12) = 8.1e4
      STATES(13) = 0
      STATES(14) = 1.2e4
      STATES(15) = 0
      STATES(16) = 0
      STATES(17) = 0
      STATES(18) = 0
      STATES(19) = 0
      STATES(20) = 0
      STATES(21) = 0
      STATES(22) = 1.1e4
      STATES(23) = 0
      STATES(24) = 2.63e4
      STATES(25) = 0
      STATES(26) = 7.2e4
      STATES(27) = 0
      STATES(28) = 0
      STATES(29) = 0
      STATES(30) = 4e4
      STATES(31) = 1.01e5
      STATES(32) = 0
      STATES(33) = 0
      STATES(34) = 0
      STATES(35) = 0
      STATES(36) = 0
      STATES(37) = 0
      STATES(38) = 0
      STATES(39) = 0
      STATES(40) = 0
      STATES(41) = 4e4
      STATES(42) = 0
      STATES(43) = 0
      STATES(44) = 4e4
      STATES(45) = 0
      STATES(46) = 0
      STATES(47) = 2.2e7
      STATES(48) = 0
      STATES(49) = 0
      STATES(50) = 0
      STATES(51) = 0
      STATES(52) = 0
      STATES(53) = 4e4
      STATES(54) = 0
      STATES(55) = 2.1e7
      STATES(56) = 0
      STATES(57) = 0
      STATES(58) = 0
      STATES(59) = 0
      STATES(60) = 1e7
      STATES(61) = 0
      STATES(62) = 0
      STATES(63) = 0
      STATES(64) = 0
      STATES(65) = 0
      STATES(66) = 0
      STATES(67) = 0
      STATES(68) = 0
      STATES(69) = 0
      STATES(70) = 0
      STATES(71) = 0
      STATES(72) = 0
      STATES(73) = 0
      STATES(74) = 0
      STATES(75) = 0
      STATES(76) = 0
      STATES(77) = 0
      STATES(78) = 0
      STATES(79) = 0
      STATES(80) = 0
      STATES(81) = 0
      STATES(82) = 0
      STATES(83) = 0
      STATES(84) = 0
      STATES(85) = 0
      STATES(86) = 0
      STATES(87) = 0
      STATES(88) = 0
      STATES(89) = 0
      STATES(90) = 0
      STATES(91) = 0
      STATES(92) = 0
      STATES(93) = 0
      STATES(94) = 0
      CONSTS(1) = 3e-3
      CONSTS(2) = 0.228
      CONSTS(3) = 1e-3
      CONSTS(4) = 6
      CONSTS(5) = 60
      CONSTS(6) = 0.6
      CONSTS(7) = 1.038e-5
      CONSTS(8) = 0.0996
      CONSTS(9) = 3e-3
      CONSTS(10) = 0.3
      CONSTS(11) = 1e-4
      CONSTS(12) = 12
      CONSTS(13) = 3.25581
      CONSTS(14) = 0.66
      CONSTS(15) = 130.2
      CONSTS(16) = 600000
      CONSTS(17) = 0.001
      CONSTS(18) = 16.5
      CONSTS(19) = 0.001
      CONSTS(20) = 3.6
      CONSTS(21) = 0.0015
      CONSTS(22) = 78
      CONSTS(23) = 30
      CONSTS(24) = 1e-5
      CONSTS(25) = 2.1e-4
      CONSTS(26) = 24
      CONSTS(27) = 1.38
      CONSTS(28) = 2.2e-5
      CONSTS(29) = 0.0021
      CONSTS(30) = 6
      CONSTS(31) = 360
      CONSTS(32) = 36
      CONSTS(33) = 33
      CONSTS(34) = 0.001
      CONSTS(35) = 1.284
      CONSTS(36) = 1e-4
      CONSTS(37) = 0.318
      CONSTS(38) = 60
      CONSTS(39) = 7e-5
      CONSTS(40) = 6
      CONSTS(41) = 2.4e-5
      CONSTS(42) = 12
      CONSTS(43) = 0.0021
      CONSTS(44) = 1.8
      CONSTS(45) = 4.5e-4
      CONSTS(46) = 0.09
      CONSTS(47) = 4.5e-4
      CONSTS(48) = 2e14
      CONSTS(49) = 61200
      CONSTS(50) = 18
      CONSTS(51) = 9e-5
      CONSTS(52) = 0.003
      CONSTS(53) = 3.84
      CONSTS(54) = 0.003
      CONSTS(55) = 2.574
      CONSTS(56) = 0.007
      CONSTS(57) = 12
      CONSTS(58) = 60
      CONSTS(59) = 0.00111
      CONSTS(60) = 1.0998
      CONSTS(61) = 210
      CONSTS(62) = 174
      CONSTS(63) = 0.00143
      CONSTS(64) = 48
      CONSTS(65) = 3.48
      CONSTS(66) = 2.5e-5
      CONSTS(67) = 30
      CONSTS(68) = 0.00534
      CONSTS(69) = 1.98
      CONSTS(70) = 960
      CONSTS(71) = 342
      CONSTS(72) = 0.00145
      CONSTS(73) = 36
      CONSTS(74) = 16.2
      CONSTS(75) = 5e-4
      CONSTS(76) = 30
      CONSTS(77) = 18
      CONSTS(78) = 0.04002
      CONSTS(79) = 0.01002
      CONSTS(80) = 50000
      CONSTS(81) =  CONSTS(15)*1.00000
      CONSTS(81) = 0.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = CONSTS(81)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(2) -  CONSTS(2)*STATES(3)
      ALGBRC(8) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      RATES(3) = ALGBRC(1) -  2.00000*ALGBRC(8)
      ALGBRC(9) =  CONSTS(5)*STATES(4) -  CONSTS(6)*STATES(5)
      RATES(4) = ALGBRC(8) - ALGBRC(9)
      ALGBRC(11) =  CONSTS(9)*STATES(2) -  CONSTS(10)*STATES(6)
      RATES(2) = CONSTS(81) - (ALGBRC(1)+ALGBRC(11))
      ALGBRC(12) =  CONSTS(9)*STATES(5)
      ALGBRC(13) =  CONSTS(11)*STATES(5)*STATES(14) -  CONSTS(12)*STATES(15)
      RATES(5) = ALGBRC(9) - (ALGBRC(12)+ALGBRC(13))
      ALGBRC(15) =  CONSTS(13)*STATES(6)*STATES(16) -  CONSTS(14)*STATES(10)
      ALGBRC(16) =  CONSTS(78)*STATES(6)
      RATES(6) = ALGBRC(11) - (ALGBRC(15)+ALGBRC(16))
      ALGBRC(17) =  CONSTS(3)*STATES(10)*STATES(10) -  CONSTS(4)*STATES(11)
      RATES(10) = ALGBRC(15) -  2.00000*ALGBRC(17)
      ALGBRC(18) =  CONSTS(79)*STATES(16)
      RATES(13) = ALGBRC(18)
      RATES(16) = - (ALGBRC(15)+ALGBRC(18))
      RATES(86) = ALGBRC(16)
      ALGBRC(19) =  CONSTS(5)*STATES(11) -  CONSTS(6)*STATES(8)
      RATES(11) = ALGBRC(17) - ALGBRC(19)
      ALGBRC(20) =  CONSTS(11)*STATES(8)*STATES(14) -  CONSTS(12)*STATES(17)
      RATES(14) = - (ALGBRC(13)+ALGBRC(20))
      ALGBRC(21) =  CONSTS(78)*STATES(8)
      RATES(8) = (ALGBRC(12)+ALGBRC(19)) - (ALGBRC(20)+ALGBRC(21))
      RATES(87) = ALGBRC(21)
      ALGBRC(10) =  CONSTS(7)*STATES(23)*STATES(12) -  CONSTS(8)*STATES(7)
      ALGBRC(23) =  CONSTS(17)*STATES(22)*STATES(15) -  CONSTS(18)*STATES(23)
      ALGBRC(14) =  CONSTS(9)*STATES(23)
      ALGBRC(24) =  CONSTS(19)*STATES(24)*STATES(23) -  CONSTS(20)*STATES(25)
      RATES(23) = ALGBRC(23) - (ALGBRC(10)+ALGBRC(14)+ALGBRC(24))
      ALGBRC(26) =  CONSTS(23)*STATES(27) -  CONSTS(24)*STATES(28)*STATES(25)
      ALGBRC(34) =  CONSTS(23)*STATES(36) -  CONSTS(24)*STATES(35)*STATES(28)
      ALGBRC(35) =  CONSTS(36)*STATES(28)*STATES(41) -  CONSTS(37)*STATES(42)
      RATES(28) = (ALGBRC(26)+ALGBRC(34)) - ALGBRC(35)
      ALGBRC(36) =  CONSTS(38)*STATES(42) -  CONSTS(39)*STATES(43)*STATES(45)
      RATES(42) = ALGBRC(35) - ALGBRC(36)
      ALGBRC(27) =  CONSTS(25)*STATES(25)*STATES(43) -  CONSTS(26)*STATES(29)
      ALGBRC(37) =  CONSTS(25)*STATES(35)*STATES(43) -  CONSTS(26)*STATES(37)
      RATES(43) = ALGBRC(36) - (ALGBRC(27)+ALGBRC(37))
      ALGBRC(50) =  CONSTS(58)*STATES(46)
      ALGBRC(49) =  CONSTS(56)*STATES(44)*STATES(45) -  CONSTS(57)*STATES(46)
      RATES(46) = ALGBRC(49) - ALGBRC(50)
      ALGBRC(51) =  CONSTS(59)*STATES(47)*STATES(45) -  CONSTS(60)*STATES(48)
      ALGBRC(52) =  CONSTS(61)*STATES(48)
      RATES(48) = ALGBRC(51) - ALGBRC(52)
      ALGBRC(53) =  CONSTS(59)*STATES(49)*STATES(45) -  CONSTS(60)*STATES(50)
      ALGBRC(54) =  CONSTS(62)*STATES(50)
      RATES(45) = (ALGBRC(36)+ALGBRC(52)+ALGBRC(54)) - (ALGBRC(49)+ALGBRC(51)+ALGBRC(53))
      RATES(50) = ALGBRC(53) - ALGBRC(54)
      ALGBRC(56) =  CONSTS(65)*STATES(52)
      ALGBRC(55) =  CONSTS(63)*STATES(51)*STATES(53) -  CONSTS(64)*STATES(52)
      RATES(52) = ALGBRC(55) - ALGBRC(56)
      ALGBRC(57) =  CONSTS(66)*STATES(53)*STATES(49) -  CONSTS(67)*STATES(54)
      RATES(49) = (ALGBRC(52)+ALGBRC(56)) - (ALGBRC(53)+ALGBRC(57))
      ALGBRC(58) =  CONSTS(65)*STATES(54)
      RATES(54) = ALGBRC(57) - ALGBRC(58)
      ALGBRC(59) =  CONSTS(68)*STATES(55)*STATES(51) -  CONSTS(69)*STATES(56)
      ALGBRC(60) =  CONSTS(70)*STATES(56)
      RATES(56) = ALGBRC(59) - ALGBRC(60)
      ALGBRC(61) =  CONSTS(68)*STATES(51)*STATES(57) -  CONSTS(69)*STATES(58)
      ALGBRC(62) =  CONSTS(71)*STATES(58)
      RATES(51) = (ALGBRC(54)+ALGBRC(60)+ALGBRC(62)) - (ALGBRC(55)+ALGBRC(59)+ALGBRC(61))
      RATES(58) = ALGBRC(61) - ALGBRC(62)
      ALGBRC(63) =  CONSTS(72)*STATES(59)*STATES(60) -  CONSTS(73)*STATES(61)
      RATES(59) = ALGBRC(62) - ALGBRC(63)
      ALGBRC(64) =  CONSTS(74)*STATES(61)
      RATES(61) = ALGBRC(63) - ALGBRC(64)
      ALGBRC(65) =  CONSTS(75)*STATES(60)*STATES(57) -  CONSTS(76)*STATES(62)
      RATES(57) = (ALGBRC(60)+ALGBRC(64)) - (ALGBRC(61)+ALGBRC(65))
      ALGBRC(66) =  CONSTS(77)*STATES(62)
      RATES(62) = ALGBRC(65) - ALGBRC(66)
      ALGBRC(29) =  CONSTS(29)*STATES(31)*STATES(15) -  CONSTS(30)*STATES(32)
      ALGBRC(73) =  CONSTS(29)*STATES(31)*STATES(17) -  CONSTS(30)*STATES(63)
      ALGBRC(43) = ( CONSTS(49)*STATES(40))/(CONSTS(48)+STATES(40))
      RATES(31) = ALGBRC(43) - (ALGBRC(29)+ALGBRC(73))
      ALGBRC(67) =  CONSTS(17)*STATES(17)*STATES(22) -  CONSTS(18)*STATES(18)
      ALGBRC(31) =  CONSTS(17)*STATES(22)*STATES(33) -  CONSTS(33)*STATES(34)
      ALGBRC(42) =  CONSTS(46)*STATES(30) -  CONSTS(47)*STATES(24)*STATES(22)
      ALGBRC(45) =  CONSTS(17)*STATES(22)*STATES(40) -  CONSTS(33)*STATES(39)
      ALGBRC(75) =  CONSTS(17)*STATES(22)*STATES(64) -  CONSTS(33)*STATES(65)
      RATES(22) = ALGBRC(42) - (ALGBRC(23)+ALGBRC(31)+ALGBRC(45)+ALGBRC(67)+ALGBRC(75))
      ALGBRC(30) =  CONSTS(31)*STATES(32) -  CONSTS(32)*STATES(33)
      ALGBRC(76) =  CONSTS(9)*STATES(32) -  CONSTS(10)*STATES(63)
      RATES(32) = ALGBRC(29) - (ALGBRC(30)+ALGBRC(76))
      ALGBRC(74) =  CONSTS(31)*STATES(63) -  CONSTS(32)*STATES(64)
      RATES(63) = (ALGBRC(73)+ALGBRC(76)) - ALGBRC(74)
      ALGBRC(68) =  CONSTS(19)*STATES(24)*STATES(18) -  CONSTS(20)*STATES(19)
      ALGBRC(77) =  CONSTS(34)*STATES(24)*STATES(65) -  CONSTS(35)*STATES(66)
      ALGBRC(32) =  CONSTS(34)*STATES(24)*STATES(34) -  CONSTS(35)*STATES(35)
      ALGBRC(47) =  CONSTS(52)*STATES(24)*STATES(39) -  CONSTS(53)*STATES(38)
      RATES(24) = ALGBRC(42) - (ALGBRC(24)+ALGBRC(32)+ALGBRC(47)+ALGBRC(68)+ALGBRC(77))
      ALGBRC(70) =  CONSTS(23)*STATES(20) -  CONSTS(24)*STATES(69)*STATES(19)
      ALGBRC(80) =  CONSTS(36)*STATES(69)*STATES(41) -  CONSTS(37)*STATES(70)
      ALGBRC(79) =  CONSTS(23)*STATES(67) -  CONSTS(24)*STATES(66)*STATES(69)
      RATES(69) = (ALGBRC(70)+ALGBRC(79)) - ALGBRC(80)
      ALGBRC(81) =  CONSTS(38)*STATES(70) -  CONSTS(39)*STATES(71)*STATES(72)
      RATES(70) = ALGBRC(80) - ALGBRC(81)
      ALGBRC(71) =  CONSTS(25)*STATES(71)*STATES(19) -  CONSTS(26)*STATES(21)
      ALGBRC(82) =  CONSTS(25)*STATES(71)*STATES(66) -  CONSTS(26)*STATES(68)
      RATES(71) = ALGBRC(81) - (ALGBRC(71)+ALGBRC(82))
      ALGBRC(69) =  CONSTS(21)*STATES(26)*STATES(19) -  CONSTS(22)*STATES(20)
      ALGBRC(72) =  CONSTS(27)*STATES(21) -  CONSTS(28)*STATES(19)*STATES(26)
      ALGBRC(25) =  CONSTS(21)*STATES(26)*STATES(25) -  CONSTS(22)*STATES(27)
      ALGBRC(28) =  CONSTS(27)*STATES(29) -  CONSTS(28)*STATES(25)*STATES(26)
      ALGBRC(33) =  CONSTS(21)*STATES(26)*STATES(35) -  CONSTS(22)*STATES(36)
      ALGBRC(38) =  CONSTS(27)*STATES(37) -  CONSTS(28)*STATES(35)*STATES(26)
      ALGBRC(84) =  CONSTS(27)*STATES(68) -  STATES(26)*CONSTS(28)*STATES(66)
      ALGBRC(78) =  CONSTS(21)*STATES(26)*STATES(66) -  CONSTS(22)*STATES(67)
      RATES(26) = (ALGBRC(28)+ALGBRC(38)+ALGBRC(72)+ALGBRC(84)) - (ALGBRC(25)+ALGBRC(33)+ALGBRC(69)+ALGBRC(78))
      ALGBRC(85) =  CONSTS(58)*STATES(73)
      RATES(41) = (ALGBRC(50)+ALGBRC(85)) - (ALGBRC(35)+ALGBRC(80))
      ALGBRC(83) =  CONSTS(56)*STATES(44)*STATES(72) -  CONSTS(57)*STATES(73)
      RATES(44) = (ALGBRC(50)+ALGBRC(85)) - (ALGBRC(49)+ALGBRC(83))
      RATES(73) = ALGBRC(83) - ALGBRC(85)
      ALGBRC(39) =  CONSTS(40)*STATES(35) -  CONSTS(41)*STATES(38)*STATES(15)
      ALGBRC(87) =  CONSTS(40)*STATES(66) -  CONSTS(41)*STATES(17)*STATES(38)
      ALGBRC(40) =  CONSTS(42)*STATES(38) -  CONSTS(43)*STATES(40)*STATES(30)
      RATES(38) = (ALGBRC(39)+ALGBRC(47)+ALGBRC(87)) - ALGBRC(40)
      ALGBRC(86) =  CONSTS(59)*STATES(47)*STATES(72) -  CONSTS(60)*STATES(74)
      ALGBRC(88) =  CONSTS(61)*STATES(74)
      RATES(74) = ALGBRC(86) - ALGBRC(88)
      ALGBRC(44) =  CONSTS(50)*STATES(33) -  CONSTS(51)*STATES(15)*STATES(40)
      ALGBRC(91) =  CONSTS(50)*STATES(64) -  CONSTS(51)*STATES(17)*STATES(40)
      RATES(40) = (ALGBRC(40)+ALGBRC(44)+ALGBRC(91)) - (ALGBRC(43)+ALGBRC(45))
      ALGBRC(90) =  CONSTS(59)*STATES(72)*STATES(75) -  CONSTS(60)*STATES(76)
      ALGBRC(92) =  CONSTS(62)*STATES(76)
      RATES(72) = (ALGBRC(81)+ALGBRC(88)+ALGBRC(92)) - (ALGBRC(83)+ALGBRC(86)+ALGBRC(90))
      RATES(76) = ALGBRC(90) - ALGBRC(92)
      ALGBRC(46) =  CONSTS(50)*STATES(34) -  CONSTS(51)*STATES(15)*STATES(39)
      ALGBRC(93) =  CONSTS(50)*STATES(65) -  CONSTS(51)*STATES(17)*STATES(39)
      RATES(39) = (ALGBRC(45)+ALGBRC(46)+ALGBRC(93)) - ALGBRC(47)
      ALGBRC(41) =  CONSTS(44)*STATES(25) -  CONSTS(45)*STATES(15)*STATES(30)
      ALGBRC(96) =  CONSTS(9)*STATES(15) -  CONSTS(10)*STATES(17)
      RATES(15) = (ALGBRC(13)+ALGBRC(39)+ALGBRC(41)+ALGBRC(44)+ALGBRC(46)) - (ALGBRC(23)+ALGBRC(29)+ALGBRC(96))
      ALGBRC(89) =  CONSTS(44)*STATES(19) -  CONSTS(45)*STATES(17)*STATES(30)
      RATES(17) = (ALGBRC(20)+ALGBRC(87)+ALGBRC(89)+ALGBRC(91)+ALGBRC(93)+ALGBRC(96)) - (ALGBRC(67)+ALGBRC(73))
      ALGBRC(48) =  CONSTS(54)*STATES(30)*STATES(33) -  CONSTS(55)*STATES(35)
      ALGBRC(95) =  CONSTS(54)*STATES(30)*STATES(64) -  CONSTS(55)*STATES(66)
      RATES(30) = (ALGBRC(40)+ALGBRC(41)+ALGBRC(89)) - (ALGBRC(42)+ALGBRC(48)+ALGBRC(95))
      ALGBRC(94) =  CONSTS(63)*STATES(77)*STATES(53) -  CONSTS(64)*STATES(78)
      ALGBRC(97) =  CONSTS(65)*STATES(78)
      RATES(78) = ALGBRC(94) - ALGBRC(97)
      ALGBRC(100) =  CONSTS(9)*STATES(33) -  CONSTS(10)*STATES(64)
      RATES(33) = ALGBRC(30) - (ALGBRC(31)+ALGBRC(44)+ALGBRC(48)+ALGBRC(100))
      RATES(64) = (ALGBRC(74)+ALGBRC(100)) - (ALGBRC(75)+ALGBRC(91)+ALGBRC(95))
      ALGBRC(101) =  CONSTS(66)*STATES(53)*STATES(75) -  CONSTS(67)*STATES(79)
      RATES(75) = (ALGBRC(88)+ALGBRC(97)) - (ALGBRC(101)+ALGBRC(90))
      ALGBRC(98) =  CONSTS(68)*STATES(55)*STATES(77) -  CONSTS(69)*STATES(80)
      ALGBRC(102) =  CONSTS(70)*STATES(80)
      RATES(80) = ALGBRC(98) - ALGBRC(102)
      ALGBRC(103) =  CONSTS(7)*STATES(25)*STATES(12) -  CONSTS(8)*STATES(88)
      ALGBRC(99) =  CONSTS(9)*STATES(25) -  CONSTS(10)*STATES(19)
      RATES(25) = (ALGBRC(24)+ALGBRC(26)+ALGBRC(28)) - (ALGBRC(25)+ALGBRC(27)+ALGBRC(41)+ALGBRC(99)+ALGBRC(103))
      ALGBRC(104) =  CONSTS(65)*STATES(79)
      RATES(47) = (ALGBRC(58)+ALGBRC(104)) - (ALGBRC(51)+ALGBRC(86))
      RATES(53) = (ALGBRC(56)+ALGBRC(58)+ALGBRC(97)+ALGBRC(104)) - (ALGBRC(55)+ALGBRC(57)+ALGBRC(94)+ALGBRC(101))
      RATES(79) = ALGBRC(101) - ALGBRC(104)
      ALGBRC(105) =  CONSTS(68)*STATES(77)*STATES(81) -  CONSTS(69)*STATES(82)
      ALGBRC(107) =  CONSTS(71)*STATES(82)
      RATES(77) = (ALGBRC(92)+ALGBRC(102)+ALGBRC(107)) - (ALGBRC(94)+ALGBRC(98)+ALGBRC(105))
      RATES(82) = ALGBRC(105) - ALGBRC(107)
      ALGBRC(108) =  CONSTS(7)*STATES(27)*STATES(12) -  CONSTS(8)*STATES(89)
      ALGBRC(106) =  CONSTS(9)*STATES(27) -  CONSTS(10)*STATES(20)
      RATES(27) = ALGBRC(25) - (ALGBRC(26)+ALGBRC(106)+ALGBRC(108))
      ALGBRC(109) =  CONSTS(72)*STATES(83)*STATES(60) -  CONSTS(73)*STATES(84)
      RATES(83) = ALGBRC(107) - ALGBRC(109)
      ALGBRC(111) =  CONSTS(74)*STATES(84)
      RATES(84) = ALGBRC(109) - ALGBRC(111)
      ALGBRC(112) =  CONSTS(7)*STATES(29)*STATES(12) -  CONSTS(8)*STATES(90)
      ALGBRC(110) =  CONSTS(9)*STATES(29) -  CONSTS(10)*STATES(21)
      RATES(29) = ALGBRC(27) - (ALGBRC(28)+ALGBRC(110)+ALGBRC(112))
      ALGBRC(113) =  CONSTS(75)*STATES(60)*STATES(81) -  CONSTS(76)*STATES(85)
      RATES(81) = (ALGBRC(102)+ALGBRC(111)) - (ALGBRC(105)+ALGBRC(113))
      ALGBRC(115) =  CONSTS(77)*STATES(85)
      RATES(55) = (ALGBRC(66)+ALGBRC(115)) - (ALGBRC(59)+ALGBRC(98))
      RATES(60) = (ALGBRC(64)+ALGBRC(66)+ALGBRC(111)+ALGBRC(115)) - (ALGBRC(63)+ALGBRC(65)+ALGBRC(109)+ALGBRC(113))
      RATES(85) = ALGBRC(113) - ALGBRC(115)
      ALGBRC(116) =  CONSTS(7)*STATES(34)*STATES(12) -  CONSTS(8)*STATES(91)
      ALGBRC(114) =  CONSTS(9)*STATES(34) -  CONSTS(10)*STATES(65)
      RATES(34) = ALGBRC(31) - (ALGBRC(32)+ALGBRC(46)+ALGBRC(114)+ALGBRC(116))
      ALGBRC(118) =  CONSTS(7)*STATES(35)*STATES(12) -  CONSTS(8)*STATES(92)
      ALGBRC(117) =  CONSTS(9)*STATES(35) -  CONSTS(10)*STATES(66)
      RATES(35) = (ALGBRC(32)+ALGBRC(34)+ALGBRC(38)+ALGBRC(48)) - (ALGBRC(117)+ALGBRC(118)+ALGBRC(33)+ALGBRC(37)+ALGBRC(39))
      ALGBRC(120) =  CONSTS(7)*STATES(36)*STATES(12) -  CONSTS(8)*STATES(93)
      ALGBRC(119) =  CONSTS(9)*STATES(36) -  CONSTS(10)*STATES(67)
      RATES(36) = ALGBRC(33) - (ALGBRC(34)+ALGBRC(119)+ALGBRC(120))
      ALGBRC(22) =  CONSTS(16)*STATES(9)
      ALGBRC(122) =  CONSTS(7)*STATES(37)*STATES(12) -  CONSTS(8)*STATES(94)
      RATES(12) = ALGBRC(22) - (ALGBRC(10)+ALGBRC(103)+ALGBRC(108)+ALGBRC(112)+ALGBRC(116)+ALGBRC(118)+ALGBRC(120)+ALGBRC(122))
      ALGBRC(121) =  CONSTS(9)*STATES(37) -  CONSTS(10)*STATES(68)
      RATES(37) = ALGBRC(37) - (ALGBRC(38)+ALGBRC(121)+ALGBRC(122))
      ALGBRC(123) = ( CONSTS(80)*1.00000)/(CONSTS(2)/( CONSTS(1)*STATES(1))+1.00000)
      ALGBRC(124) = TERNRY(ALGBRC(123).LT.3100.00, 1.55000, TERNRY(ALGBRC(123).GT.100000., 0.200000, ( ALGBRC(123)*- 1.35000e-05+1.55000)/1.00000)
      ALGBRC(125) =  ALGBRC(124)*STATES(7)
      RATES(7) = ALGBRC(10) - ALGBRC(125)
      RATES(18) = (ALGBRC(125)+ALGBRC(14)+ALGBRC(67)) - ALGBRC(68)
      ALGBRC(126) =  ALGBRC(124)*STATES(88)
      RATES(19) = (ALGBRC(70)+ALGBRC(68)+ALGBRC(72)+ALGBRC(99)+ALGBRC(126)) - (ALGBRC(69)+ALGBRC(71)+ALGBRC(89))
      RATES(88) = ALGBRC(103) - ALGBRC(126)
      ALGBRC(127) =  ALGBRC(124)*STATES(89)
      RATES(20) = (ALGBRC(69)+ALGBRC(106)+ALGBRC(127)) - ALGBRC(70)
      RATES(89) = ALGBRC(108) - ALGBRC(127)
      ALGBRC(128) =  ALGBRC(124)*STATES(90)
      RATES(21) = (ALGBRC(71)+ALGBRC(110)+ALGBRC(128)) - ALGBRC(72)
      RATES(90) = ALGBRC(112) - ALGBRC(128)
      ALGBRC(129) =  ALGBRC(124)*STATES(91)
      RATES(65) = (ALGBRC(75)+ALGBRC(114)+ALGBRC(129)) - (ALGBRC(77)+ALGBRC(93))
      RATES(91) = ALGBRC(116) - ALGBRC(129)
      ALGBRC(130) =  ALGBRC(124)*STATES(92)
      RATES(66) = (ALGBRC(77)+ALGBRC(79)+ALGBRC(84)+ALGBRC(95)+ALGBRC(117)+ALGBRC(130)) - (ALGBRC(78)+ALGBRC(82)+ALGBRC(87))
      RATES(92) = ALGBRC(118) - ALGBRC(130)
      ALGBRC(131) =  ALGBRC(124)*STATES(93)
      RATES(67) = (ALGBRC(78)+ALGBRC(119)+ALGBRC(131)) - ALGBRC(79)
      RATES(93) = ALGBRC(120) - ALGBRC(131)
      ALGBRC(132) =  ALGBRC(124)*STATES(94)
      RATES(9) = (ALGBRC(125)+ALGBRC(126)+ALGBRC(127)+ALGBRC(128)+ALGBRC(129)+ALGBRC(130)+ALGBRC(131)+ALGBRC(132)) - ALGBRC(22)
      RATES(68) = (ALGBRC(82)+ALGBRC(121)+ALGBRC(132)) - ALGBRC(84)
      RATES(94) = ALGBRC(122) - ALGBRC(132)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(2) -  CONSTS(2)*STATES(3)
      ALGBRC(8) =  CONSTS(3)*STATES(3)*STATES(3) -  CONSTS(4)*STATES(4)
      ALGBRC(9) =  CONSTS(5)*STATES(4) -  CONSTS(6)*STATES(5)
      ALGBRC(11) =  CONSTS(9)*STATES(2) -  CONSTS(10)*STATES(6)
      ALGBRC(12) =  CONSTS(9)*STATES(5)
      ALGBRC(13) =  CONSTS(11)*STATES(5)*STATES(14) -  CONSTS(12)*STATES(15)
      ALGBRC(15) =  CONSTS(13)*STATES(6)*STATES(16) -  CONSTS(14)*STATES(10)
      ALGBRC(16) =  CONSTS(78)*STATES(6)
      ALGBRC(17) =  CONSTS(3)*STATES(10)*STATES(10) -  CONSTS(4)*STATES(11)
      ALGBRC(18) =  CONSTS(79)*STATES(16)
      ALGBRC(19) =  CONSTS(5)*STATES(11) -  CONSTS(6)*STATES(8)
      ALGBRC(20) =  CONSTS(11)*STATES(8)*STATES(14) -  CONSTS(12)*STATES(17)
      ALGBRC(21) =  CONSTS(78)*STATES(8)
      ALGBRC(10) =  CONSTS(7)*STATES(23)*STATES(12) -  CONSTS(8)*STATES(7)
      ALGBRC(23) =  CONSTS(17)*STATES(22)*STATES(15) -  CONSTS(18)*STATES(23)
      ALGBRC(14) =  CONSTS(9)*STATES(23)
      ALGBRC(24) =  CONSTS(19)*STATES(24)*STATES(23) -  CONSTS(20)*STATES(25)
      ALGBRC(26) =  CONSTS(23)*STATES(27) -  CONSTS(24)*STATES(28)*STATES(25)
      ALGBRC(34) =  CONSTS(23)*STATES(36) -  CONSTS(24)*STATES(35)*STATES(28)
      ALGBRC(35) =  CONSTS(36)*STATES(28)*STATES(41) -  CONSTS(37)*STATES(42)
      ALGBRC(36) =  CONSTS(38)*STATES(42) -  CONSTS(39)*STATES(43)*STATES(45)
      ALGBRC(27) =  CONSTS(25)*STATES(25)*STATES(43) -  CONSTS(26)*STATES(29)
      ALGBRC(37) =  CONSTS(25)*STATES(35)*STATES(43) -  CONSTS(26)*STATES(37)
      ALGBRC(50) =  CONSTS(58)*STATES(46)
      ALGBRC(49) =  CONSTS(56)*STATES(44)*STATES(45) -  CONSTS(57)*STATES(46)
      ALGBRC(51) =  CONSTS(59)*STATES(47)*STATES(45) -  CONSTS(60)*STATES(48)
      ALGBRC(52) =  CONSTS(61)*STATES(48)
      ALGBRC(53) =  CONSTS(59)*STATES(49)*STATES(45) -  CONSTS(60)*STATES(50)
      ALGBRC(54) =  CONSTS(62)*STATES(50)
      ALGBRC(56) =  CONSTS(65)*STATES(52)
      ALGBRC(55) =  CONSTS(63)*STATES(51)*STATES(53) -  CONSTS(64)*STATES(52)
      ALGBRC(57) =  CONSTS(66)*STATES(53)*STATES(49) -  CONSTS(67)*STATES(54)
      ALGBRC(58) =  CONSTS(65)*STATES(54)
      ALGBRC(59) =  CONSTS(68)*STATES(55)*STATES(51) -  CONSTS(69)*STATES(56)
      ALGBRC(60) =  CONSTS(70)*STATES(56)
      ALGBRC(61) =  CONSTS(68)*STATES(51)*STATES(57) -  CONSTS(69)*STATES(58)
      ALGBRC(62) =  CONSTS(71)*STATES(58)
      ALGBRC(63) =  CONSTS(72)*STATES(59)*STATES(60) -  CONSTS(73)*STATES(61)
      ALGBRC(64) =  CONSTS(74)*STATES(61)
      ALGBRC(65) =  CONSTS(75)*STATES(60)*STATES(57) -  CONSTS(76)*STATES(62)
      ALGBRC(66) =  CONSTS(77)*STATES(62)
      ALGBRC(29) =  CONSTS(29)*STATES(31)*STATES(15) -  CONSTS(30)*STATES(32)
      ALGBRC(73) =  CONSTS(29)*STATES(31)*STATES(17) -  CONSTS(30)*STATES(63)
      ALGBRC(43) = ( CONSTS(49)*STATES(40))/(CONSTS(48)+STATES(40))
      ALGBRC(67) =  CONSTS(17)*STATES(17)*STATES(22) -  CONSTS(18)*STATES(18)
      ALGBRC(31) =  CONSTS(17)*STATES(22)*STATES(33) -  CONSTS(33)*STATES(34)
      ALGBRC(42) =  CONSTS(46)*STATES(30) -  CONSTS(47)*STATES(24)*STATES(22)
      ALGBRC(45) =  CONSTS(17)*STATES(22)*STATES(40) -  CONSTS(33)*STATES(39)
      ALGBRC(75) =  CONSTS(17)*STATES(22)*STATES(64) -  CONSTS(33)*STATES(65)
      ALGBRC(30) =  CONSTS(31)*STATES(32) -  CONSTS(32)*STATES(33)
      ALGBRC(76) =  CONSTS(9)*STATES(32) -  CONSTS(10)*STATES(63)
      ALGBRC(74) =  CONSTS(31)*STATES(63) -  CONSTS(32)*STATES(64)
      ALGBRC(68) =  CONSTS(19)*STATES(24)*STATES(18) -  CONSTS(20)*STATES(19)
      ALGBRC(77) =  CONSTS(34)*STATES(24)*STATES(65) -  CONSTS(35)*STATES(66)
      ALGBRC(32) =  CONSTS(34)*STATES(24)*STATES(34) -  CONSTS(35)*STATES(35)
      ALGBRC(47) =  CONSTS(52)*STATES(24)*STATES(39) -  CONSTS(53)*STATES(38)
      ALGBRC(70) =  CONSTS(23)*STATES(20) -  CONSTS(24)*STATES(69)*STATES(19)
      ALGBRC(80) =  CONSTS(36)*STATES(69)*STATES(41) -  CONSTS(37)*STATES(70)
      ALGBRC(79) =  CONSTS(23)*STATES(67) -  CONSTS(24)*STATES(66)*STATES(69)
      ALGBRC(81) =  CONSTS(38)*STATES(70) -  CONSTS(39)*STATES(71)*STATES(72)
      ALGBRC(71) =  CONSTS(25)*STATES(71)*STATES(19) -  CONSTS(26)*STATES(21)
      ALGBRC(82) =  CONSTS(25)*STATES(71)*STATES(66) -  CONSTS(26)*STATES(68)
      ALGBRC(69) =  CONSTS(21)*STATES(26)*STATES(19) -  CONSTS(22)*STATES(20)
      ALGBRC(72) =  CONSTS(27)*STATES(21) -  CONSTS(28)*STATES(19)*STATES(26)
      ALGBRC(25) =  CONSTS(21)*STATES(26)*STATES(25) -  CONSTS(22)*STATES(27)
      ALGBRC(28) =  CONSTS(27)*STATES(29) -  CONSTS(28)*STATES(25)*STATES(26)
      ALGBRC(33) =  CONSTS(21)*STATES(26)*STATES(35) -  CONSTS(22)*STATES(36)
      ALGBRC(38) =  CONSTS(27)*STATES(37) -  CONSTS(28)*STATES(35)*STATES(26)
      ALGBRC(84) =  CONSTS(27)*STATES(68) -  STATES(26)*CONSTS(28)*STATES(66)
      ALGBRC(78) =  CONSTS(21)*STATES(26)*STATES(66) -  CONSTS(22)*STATES(67)
      ALGBRC(85) =  CONSTS(58)*STATES(73)
      ALGBRC(83) =  CONSTS(56)*STATES(44)*STATES(72) -  CONSTS(57)*STATES(73)
      ALGBRC(39) =  CONSTS(40)*STATES(35) -  CONSTS(41)*STATES(38)*STATES(15)
      ALGBRC(87) =  CONSTS(40)*STATES(66) -  CONSTS(41)*STATES(17)*STATES(38)
      ALGBRC(40) =  CONSTS(42)*STATES(38) -  CONSTS(43)*STATES(40)*STATES(30)
      ALGBRC(86) =  CONSTS(59)*STATES(47)*STATES(72) -  CONSTS(60)*STATES(74)
      ALGBRC(88) =  CONSTS(61)*STATES(74)
      ALGBRC(44) =  CONSTS(50)*STATES(33) -  CONSTS(51)*STATES(15)*STATES(40)
      ALGBRC(91) =  CONSTS(50)*STATES(64) -  CONSTS(51)*STATES(17)*STATES(40)
      ALGBRC(90) =  CONSTS(59)*STATES(72)*STATES(75) -  CONSTS(60)*STATES(76)
      ALGBRC(92) =  CONSTS(62)*STATES(76)
      ALGBRC(46) =  CONSTS(50)*STATES(34) -  CONSTS(51)*STATES(15)*STATES(39)
      ALGBRC(93) =  CONSTS(50)*STATES(65) -  CONSTS(51)*STATES(17)*STATES(39)
      ALGBRC(41) =  CONSTS(44)*STATES(25) -  CONSTS(45)*STATES(15)*STATES(30)
      ALGBRC(96) =  CONSTS(9)*STATES(15) -  CONSTS(10)*STATES(17)
      ALGBRC(89) =  CONSTS(44)*STATES(19) -  CONSTS(45)*STATES(17)*STATES(30)
      ALGBRC(48) =  CONSTS(54)*STATES(30)*STATES(33) -  CONSTS(55)*STATES(35)
      ALGBRC(95) =  CONSTS(54)*STATES(30)*STATES(64) -  CONSTS(55)*STATES(66)
      ALGBRC(94) =  CONSTS(63)*STATES(77)*STATES(53) -  CONSTS(64)*STATES(78)
      ALGBRC(97) =  CONSTS(65)*STATES(78)
      ALGBRC(100) =  CONSTS(9)*STATES(33) -  CONSTS(10)*STATES(64)
      ALGBRC(101) =  CONSTS(66)*STATES(53)*STATES(75) -  CONSTS(67)*STATES(79)
      ALGBRC(98) =  CONSTS(68)*STATES(55)*STATES(77) -  CONSTS(69)*STATES(80)
      ALGBRC(102) =  CONSTS(70)*STATES(80)
      ALGBRC(103) =  CONSTS(7)*STATES(25)*STATES(12) -  CONSTS(8)*STATES(88)
      ALGBRC(99) =  CONSTS(9)*STATES(25) -  CONSTS(10)*STATES(19)
      ALGBRC(104) =  CONSTS(65)*STATES(79)
      ALGBRC(105) =  CONSTS(68)*STATES(77)*STATES(81) -  CONSTS(69)*STATES(82)
      ALGBRC(107) =  CONSTS(71)*STATES(82)
      ALGBRC(108) =  CONSTS(7)*STATES(27)*STATES(12) -  CONSTS(8)*STATES(89)
      ALGBRC(106) =  CONSTS(9)*STATES(27) -  CONSTS(10)*STATES(20)
      ALGBRC(109) =  CONSTS(72)*STATES(83)*STATES(60) -  CONSTS(73)*STATES(84)
      ALGBRC(111) =  CONSTS(74)*STATES(84)
      ALGBRC(112) =  CONSTS(7)*STATES(29)*STATES(12) -  CONSTS(8)*STATES(90)
      ALGBRC(110) =  CONSTS(9)*STATES(29) -  CONSTS(10)*STATES(21)
      ALGBRC(113) =  CONSTS(75)*STATES(60)*STATES(81) -  CONSTS(76)*STATES(85)
      ALGBRC(115) =  CONSTS(77)*STATES(85)
      ALGBRC(116) =  CONSTS(7)*STATES(34)*STATES(12) -  CONSTS(8)*STATES(91)
      ALGBRC(114) =  CONSTS(9)*STATES(34) -  CONSTS(10)*STATES(65)
      ALGBRC(118) =  CONSTS(7)*STATES(35)*STATES(12) -  CONSTS(8)*STATES(92)
      ALGBRC(117) =  CONSTS(9)*STATES(35) -  CONSTS(10)*STATES(66)
      ALGBRC(120) =  CONSTS(7)*STATES(36)*STATES(12) -  CONSTS(8)*STATES(93)
      ALGBRC(119) =  CONSTS(9)*STATES(36) -  CONSTS(10)*STATES(67)
      ALGBRC(22) =  CONSTS(16)*STATES(9)
      ALGBRC(122) =  CONSTS(7)*STATES(37)*STATES(12) -  CONSTS(8)*STATES(94)
      ALGBRC(121) =  CONSTS(9)*STATES(37) -  CONSTS(10)*STATES(68)
      ALGBRC(123) = ( CONSTS(80)*1.00000)/(CONSTS(2)/( CONSTS(1)*STATES(1))+1.00000)
      ALGBRC(124) = TERNRY(ALGBRC(123).LT.3100.00, 1.55000, TERNRY(ALGBRC(123).GT.100000., 0.200000, ( ALGBRC(123)*- 1.35000e-05+1.55000)/1.00000)
      ALGBRC(125) =  ALGBRC(124)*STATES(7)
      ALGBRC(126) =  ALGBRC(124)*STATES(88)
      ALGBRC(127) =  ALGBRC(124)*STATES(89)
      ALGBRC(128) =  ALGBRC(124)*STATES(90)
      ALGBRC(129) =  ALGBRC(124)*STATES(91)
      ALGBRC(130) =  ALGBRC(124)*STATES(92)
      ALGBRC(131) =  ALGBRC(124)*STATES(93)
      ALGBRC(132) =  ALGBRC(124)*STATES(94)
      ALGBRC(2) = STATES(45)+STATES(46)+STATES(48)+STATES(50)+STATES(72)+STATES(73)+STATES(74)+STATES(76)
      ALGBRC(3) = STATES(51)+STATES(77)
      ALGBRC(4) = STATES(59)+STATES(83)
      ALGBRC(5) = STATES(42)+STATES(28)+STATES(70)+STATES(69)
      ALGBRC(6) = STATES(33)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(38)+STATES(39)+STATES(40)+STATES(64)+STATES(65)+STATES(66)+STATES(67)+STATES(68)+STATES(91)+STATES(92)+STATES(93)+STATES(94)
      ALGBRC(7) =  2.00000*(STATES(5)+STATES(7)+STATES(8)+STATES(11)+STATES(15)+STATES(17)+STATES(18)+STATES(19)+STATES(20)+STATES(21)+STATES(23)+STATES(25)+STATES(27)+STATES(29)+STATES(32)+STATES(33)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(63)+STATES(64)+STATES(65)+STATES(66)+STATES(67)+STATES(68)+STATES(88)+STATES(89)+STATES(90)+STATES(91)+STATES(92)+STATES(93)+STATES(94))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END