diff --git a/scatcluster/processing/scattering.py b/scatcluster/processing/scattering.py
index a5e9750..e9e0436 100644
--- a/scatcluster/processing/scattering.py
+++ b/scatcluster/processing/scattering.py
@@ -545,7 +545,7 @@ def process_vectorized_scattering_coefficients(self) -> None:
             coefficients.order_1.sum(dim=('f1', 'channel')) > 0,
             drop=True,
         )
-        coefficients = self.log(coefficients, waterlevel=1e-20)
+        coefficients = self.log(coefficients, waterlevel=1e-15)
         coefficients = self.nyquist_mask(coefficients)
         coefficients = self.normalize(coefficients)
         print(coefficients)