Fix error in enthalpy_second_derivatives h_sa_sa

castelao · Aug 4, 2024 · 11a18fd · 11a18fd
1 parent 37f157c
commit 11a18fd
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/volume.rs b/src/volume.rs
@@ -2098,7 +2098,7 @@ pub fn enthalpy_second_derivatives(sa: f64, ct: f64, p: f64) -> Result<(f64, f64
                         + 3.0 * H303 * s2
                         + tau * (-H113 + tau * (-H123 - H133 * tau))
                         + pi * (-H104 - H114 * tau - H105 * pi)))));
-    let h_sa_sa = 1e8 * 0.25 * GSW_SFAC * dynamic_h_sa_sa_part / (s * s * s);
+    let h_sa_sa = 1e8 * 0.25 * GSW_SFAC * GSW_SFAC * dynamic_h_sa_sa_part / (s * s * s);
 
     let dynamic_h_sa_ct_part = pi
         * (H111