Use thermovariabledict to access volume

ReactionMechanismGenerator · Jan 23, 2024 · 2624531 · 2624531
1 parent 4997f01
commit 2624531
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Interface.jl b/src/Interface.jl
@@ -534,7 +534,7 @@ function evaluate(ri::FragmentBasedReactiveFilmGrowthInterfaceConstantT, dydt, d
     addreactionratecontributions!(dydt, ri.fragmentbasedrxnarray, ri.rxnarray, cstot, kfs, krevs, Vfilm, domainfilm.indexes[3], ri.Mws, domainfilm.indexes[1]:domainfilm.indexes[2])
     if hasproperty(ri.domain2, :epsilon)
         epsilon = ri.domain2.epsilon
-        dydt[ri.domain2.indexes[3]] = dydt[domainfilm.indexes[3]] / domainfilm.rho / (1 - epsilon) * epsilon
+        dydt[ri.domain2.thermovariabledict["V"]] = dydt[domainfilm.indexes[3]] / domainfilm.rho / (1 - epsilon) * epsilon
     end
 end