diff --git a/PySDM/backends/impl_thrust_rtc/test_helpers/cpp2python.py b/PySDM/backends/impl_thrust_rtc/test_helpers/cpp2python.py
index 8d3e7623bc..00b37421fa 100644
--- a/PySDM/backends/impl_thrust_rtc/test_helpers/cpp2python.py
+++ b/PySDM/backends/impl_thrust_rtc/test_helpers/cpp2python.py
@@ -197,7 +197,7 @@ def make(self):
         power, sqrt,
         arctanh as atanh,
         arcsinh as asinh,
-        sinh,
+        sinh, tanh,
         maximum, minimum,
         where,  # TODO #1295
     )
diff --git a/PySDM/formulae.py b/PySDM/formulae.py
index 4b796f5784..9667c243de 100644
--- a/PySDM/formulae.py
+++ b/PySDM/formulae.py
@@ -31,7 +31,7 @@ def __init__(  # pylint: disable=too-many-locals
         seed: int = None,
         fastmath: bool = True,
         diffusion_coordinate: str = "WaterMassLogarithm",
-        saturation_vapour_pressure: str = "FlatauWalkoCotton",
+        saturation_vapour_pressure: str = "MurphyKoop2005",
         latent_heat_vapourisation: str = "Kirchhoff",
         latent_heat_sublimation: str = "MurphyKoop2005",
         hygroscopicity: str = "KappaKoehlerLeadingTerms",
diff --git a/tests/unit_tests/dynamics/condensation/test_diagnostics.py b/tests/unit_tests/dynamics/condensation/test_diagnostics.py
index 4c8f01cb9b..c034e4e954 100644
--- a/tests/unit_tests/dynamics/condensation/test_diagnostics.py
+++ b/tests/unit_tests/dynamics/condensation/test_diagnostics.py
@@ -5,7 +5,7 @@
 
 import numpy as np
 
-from PySDM import Builder
+from PySDM import Builder, Formulae
 from PySDM.backends import CPU
 from PySDM.dynamics.condensation import Condensation
 from PySDM.impl.mesh import Mesh
@@ -79,7 +79,8 @@ def __init__(  # pylint: disable=too-many-locals
             rho=rho,
             eta=eta,
         )
-        builder = Builder(n_sd=n_sd, backend=backend(), environment=env)
+        formulae = Formulae(saturation_vapour_pressure="FlatauWalkoCotton")
+        builder = Builder(n_sd=n_sd, backend=backend(formulae), environment=env)
 
         builder.add_dynamic(Condensation(max_iters=max_iters))
         self.particulator = builder.build(
diff --git a/tests/unit_tests/dynamics/test_vapour_deposition_on_ice.py b/tests/unit_tests/dynamics/test_vapour_deposition_on_ice.py
index 819e0fee0a..7e9644cef3 100644
--- a/tests/unit_tests/dynamics/test_vapour_deposition_on_ice.py
+++ b/tests/unit_tests/dynamics/test_vapour_deposition_on_ice.py
@@ -234,9 +234,9 @@ def test_growth_rates_against_spichtinger_and_gierens_2009_fig_5(
             assert (np.diff(mass_rate) > 0).all()
 
         assert 1.8e-15 * si.kg / si.s < dm_dt[230 * si.K][0] < 4.0e-15 * si.kg / si.s
-        assert 7.0e-17 * si.kg / si.s < dm_dt[200 * si.K][0] < 1.0e-16 * si.kg / si.s
+        assert 6.0e-17 * si.kg / si.s < dm_dt[200 * si.K][0] < 1.0e-16 * si.kg / si.s
         assert 1.3e-12 * si.kg / si.s < dm_dt[230 * si.K][-1] < 1.5e-12 * si.kg / si.s
-        assert 6.0e-14 * si.kg / si.s < dm_dt[200 * si.K][-1] < 1.2e-13 * si.kg / si.s
+        assert 5.0e-14 * si.kg / si.s < dm_dt[200 * si.K][-1] < 1.2e-13 * si.kg / si.s
 
     @staticmethod
     @pytest.mark.parametrize("diffusion_coordinate", DIFFUSION_COORDINATES)