Add mass in rop matrix for film domain; start should be 0 not 1 becau…

…se it's used as start+i
ReactionMechanismGenerator · Aug 8, 2023 · 12072da · 12072da
1 parent 0941aca
commit 12072da
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/Simulation.jl b/src/Simulation.jl
@@ -271,8 +271,12 @@ function rops(ssys::SystemSimulation, t)
     vdiffs = Array{Any,1}(undef, length(domains))
     vCvave = Array{Any,1}(undef, length(domains))
     vphi = Array{Any,1}(undef, length(domains))
-    ropmat = spzeros(Nrxns, Nspcs)
-    start = 1
+    if any([domain isa FragmentBasedConstantTrhoDomain for domain in domains])
+        ropmat = spzeros(Nrxns, Nspcs + 1)
+    else
+        ropmat = spzeros(Nrxns, Nspcs)
+    end
+    start = 0
     for (k, sim) in enumerate(ssys.sims)
         vns[k], vcs[k], vT[k], vP[k], vV[k], vC[k], vN[k], vmu[k], vkfs[k], vkrevs[k], vHs[k], vUs[k], vGs[k], vdiffs[k], vCvave[k], vphi[k] = calcthermo(sim.domain, ssys.sol(t), t)
         cstot[sim.domain.indexes[1]:sim.domain.indexes[2]] = vcs[k]