add pph mixing modes: lake&inmom constants

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)