diff --git a/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py b/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
index 1490cfcb..bd451d28 100644
--- a/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
+++ b/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
@@ -54,7 +54,7 @@ def drhod_dz(z, rhod):
             if not apprx_drhod_dz:  # to resolve issue #335
                 qv = self.qv(z)
                 dqv_dz = Derivative(self.qv)(z)
-                drhod_dz = drhod_dz / (1 + qv) - rhod * dqv_dz / (1 + qv)
+                drhod_dz = drhod_dz / (1 + qv) - rhod * dqv_dz / (1 + qv)**2
             return drhod_dz
 
         z_points = np.arange(0, self.z_max + self.dz / 2, self.dz / 2)