From a95bdd5f98df4de842231be96cb21d93c3da641d Mon Sep 17 00:00:00 2001
From: Ramil Ahtamyanov <ramil.rgk@gmail.com>
Date: Sat, 21 Dec 2024 03:11:13 +0300
Subject: [PATCH] added mix mode #5

---
 obl_inmom.f90 | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

diff --git a/obl_inmom.f90 b/obl_inmom.f90
index f70f53c..5eb28dd 100644
--- a/obl_inmom.f90
+++ b/obl_inmom.f90
@@ -164,6 +164,23 @@ module obl_inmom
       end do
       kh = kh * 10000.0
       km = km * 10000.0
+    ! obl_pph mixing mode (dasha constants)
+    else if (kh_km_mode == 5) then
+      pphParams%Km_0 = 7.0 * 0.01
+      pphParams%Kh_0 = 5.0 * 0.01
+      pphParams%alpha = 25.0 / 7.0 !Nuzhno tak!
+      pphParams%Kh_unstable = 0.05
+      pphParams%Km_unstable = 0.05
+      do j = 1, ny
+        do i = 1, nx
+          if (lu(i, j) > lu_min) then
+            call pph_kh(kh(i,j,:), rit(i,j,:), pphParams, nz)
+            call pph_km(km(i,j,:), rit(i,j,:), pphParams, nz)
+          end if
+        end do
+      end do
+      kh = kh * 10000.0
+      km = km * 10000.0
     end if
     
     ! print *, "Kh first:", kh(3,3,1:4)
-- 
GitLab