Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 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