add comment

CusiniM · CusiniM · commit 9e36eddbec12 · 2025-06-24T10:57:23.000-07:00
diff --git a/src/reactions/geochemistry/Forge.hpp b/src/reactions/geochemistry/Forge.hpp
@@ -37,6 +37,7 @@ constexpr CArrayWrapper< double, 19, 26 > soichMatrix =
   {  0,     0,     0,    0,    0,     0,     0,    0,     0,     0,     0,     0,     0,     0,   0,    0,     -4,   0,    0,     1,    0,    1,     0,     1,     0,   3  }  // Albite: NaAlSi₃O₈(s) + 4H⁺ ⇌ Al³⁺ + Na⁺ + 3SiO₂(aq)
 };
 
+// Must convert these. They should not be the log.
 constexpr CArrayWrapper< double, 19 > equilibriumConstants =
 {
    5.9636,   // CaCO₃(aq) + H⁺ ⇌ Ca²⁺ + HCO₃⁻
@@ -107,9 +108,6 @@ constexpr CArrayWrapper< double, 19 > reverseRateConstant =
    1.0    // Albite: NaAlSi₃O₈(s) + 4H⁺ ⇌ Al³⁺ + Na⁺ + 3SiO₂(aq)
 };
 
-
-
-
 }
 
 using forgeSystemType = reactionsSystems::MixedReactionsParameters< double, int, int, 26, 19, 16 >;

Original file line number	Diff line number	Diff line change
`@@ -37,6 +37,7 @@ constexpr CArrayWrapper< double, 19, 26 > soichMatrix =`
`37`	`37`	`{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -4, 0, 0, 1, 0, 1, 0, 1, 0, 3 } // Albite: NaAlSi₃O₈(s) + 4H⁺ ⇌ Al³⁺ + Na⁺ + 3SiO₂(aq)`
`38`	`38`	`};`
`39`	`39`
	`40`	`+// Must convert these. They should not be the log.`
`40`	`41`	`constexpr CArrayWrapper< double, 19 > equilibriumConstants =`
`41`	`42`	`{`
`42`	`43`	`5.9636, // CaCO₃(aq) + H⁺ ⇌ Ca²⁺ + HCO₃⁻`
`@@ -107,9 +108,6 @@ constexpr CArrayWrapper< double, 19 > reverseRateConstant =`
`107`	`108`	`1.0 // Albite: NaAlSi₃O₈(s) + 4H⁺ ⇌ Al³⁺ + Na⁺ + 3SiO₂(aq)`
`108`	`109`	`};`
`109`	`110`
`110`		`-`
`111`		`-`
`112`		`-`
`113`	`111`	`}`
`114`	`112`
`115`	`113`	`using forgeSystemType = reactionsSystems::MixedReactionsParameters< double, int, int, 26, 19, 16 >;`