Mercurial > bmsapp / file diff

--- a/src/Utils.cpp	Fri Nov 18 12:17:32 2022 +0100
+++ b/src/Utils.cpp	Fri Nov 18 16:57:02 2022 +0100
@@ -771,6 +771,22 @@
 }
 
 
+double Utils::ResCO2(double T)
+{
+    double F = T * 1.8 + 32;
+    return round((3.0378 - 0.050062 * F + 0.00026555 * F * F) * 1000000.0) / 1000000.0;
+}
+
+
+double Utils::CarbCO2toS(double CO2, double T, double SFactor)
+{
+    double sugar = round((SFactor * (CO2 - ResCO2(T)) * 4.014094) * 1000000.0) / 1000000.0;
+    if (sugar < 0)
+        sugar = 0;
+    return sugar;
+}
+
+
 double Utils::GetPressureBar(double gl, double T)
 {
     if (gl <= 0)