added mix mode #5

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)