From ae2c21e83c8d64e3874efc72549d298c39158541 Mon Sep 17 00:00:00 2001
From: Ramil Ahtamyanov <ramil.rgk@gmail.com>
Date: Sat, 21 Dec 2024 02:25:08 +0300
Subject: [PATCH] add pph mixing modes: lake&inmom constants
---
obl_inmom.f90 | 18 ++++++++++++++++--
1 file changed, 16 insertions(+), 2 deletions(-)
diff --git a/obl_inmom.f90 b/obl_inmom.f90
index b51725f..d51587a 100644
--- a/obl_inmom.f90
+++ b/obl_inmom.f90
@@ -116,8 +116,20 @@ module obl_inmom
end if
end do
end do
- ! obl_pph mixing mode
+ ! obl_pph mixing mode (lake constants)
else if (kh_km_mode == 2) then
+ 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
+ ! obl_pph mixing mode (inmom constants)
+ else if (kh_km_mode == 3) then
pphParams%Km_0 = 7.0 * 0.01
pphParams%Kh_0 = 5.0 * 0.01
pphParams%alpha = 5.0
@@ -135,7 +147,7 @@ module obl_inmom
kh = kh * 10000.0
km = km * 10000.0
! obl_pph_dyn mixing mode
- else if (kh_km_mode == 3) then
+ else if (kh_km_mode == 4) then
do j = 1, ny
do i = 1, nx
if (lu(i, j) > lu_min) then
@@ -144,6 +156,8 @@ module obl_inmom
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