diff --git a/obl_legacy.f90 b/obl_legacy.f90
index c062958d18ea11febfaa19d89bf39c1908d8eb56..284c5f260c64529f52039da22d9a15d851b70b81 100644
--- a/obl_legacy.f90
+++ b/obl_legacy.f90
@@ -301,7 +301,7 @@ module obl_legacy
 
     do k = 2, nz
       if (rit(k) >= 0.0) then
-        kh(k) = kh_0 / (1.0 + 5.0 * rit(k)) + kh_b(k)
+        kh(k) = kh_0 / (1.0 + 5.0 * rit(k))**3 + kh_b(k)
       else
         kh(k) = unstable
       end if