diff --git a/src/volume.rs b/src/volume.rs index 1f69f6c..09fc9af 100644 --- a/src/volume.rs +++ b/src/volume.rs @@ -2076,7 +2076,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