From 3b1b66bdc61a6f243156bc1bac6774d20d31ca1f Mon Sep 17 00:00:00 2001
From: Ramil Ahtamyanov <ramil.rgk@gmail.com>
Date: Tue, 24 Dec 2024 19:42:41 +0300
Subject: [PATCH] kh fix. add **3
---
obl_legacy.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/obl_legacy.f90 b/obl_legacy.f90
index c062958..284c5f2 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
--
GitLab