diff --git a/obl_inmom.f90 b/obl_inmom.f90 index f70f53c548129e95c2b51f70975c6fd4eb54ccfb..5eb28ddb940ea08e0d61d42432072ea4757edb13 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)