output corespeciesnetconsumptionrates, coreradicalnetterminationrates

hwpang · hwpang · commit 1a761f08644b · 2023-04-18T18:54:17.000-04:00
diff --git a/src/EdgeAnalysis.jl b/src/EdgeAnalysis.jl
@@ -670,7 +670,8 @@ function identifyobjects!(sim,corespcsinds,corerxninds,edgespcsinds,
     (dydt,rts,frts,rrts,cs,corespeciesrates,charrate,edgespeciesrates,
     edgereactionrates,edgerxnradrateratios,corespeciesrateratios,edgespeciesrateratios,
     corereactionrates,corespeciesconcentrations,
-    corespeciesproductionrates,corespeciesconsumptionrates) = processfluxes(sim,corespcsinds,corerxninds,edgespcsinds,edgerxninds)
+    corespeciesproductionrates,corespeciesconsumptionrates,
+    corespeciesnetconsumptionrates, coreradicalnetterminationrates) = processfluxes(sim,corespcsinds,corerxninds,edgespcsinds,edgerxninds)
 
     for i = 1:length(edgespeciesrateratios)
         if @inbounds  edgespeciesrateratios[i] > maxedgespeciesrateratios[i]
