Gradient model is updated

mom-ocean · Oct 18, 2024 · b90c189 · b90c189
1 parent d8e2f3c
commit b90c189
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/parameterizations/lateral/MOM_thickness_diffuse.F90 b/parameterizations/lateral/MOM_thickness_diffuse.F90
@@ -318,7 +318,8 @@ subroutine thickness_diffuse(h, uhtr, vhtr, tv, dt, G, GV, US, MEKE, VarMix, CDp
     if (CS%Grad_L_Scale > 0.0) then
       !$OMP do
       do k=1,nz ; do j=js,je ; do I=is-1,ie
-        KH_u(I,j,k) = CS%Grad_L_Scale*VarMix%L2grad_u(I,j)*VarMix%UH_grad(I,j,k)
+        KH_u(I,j,k) = 1.0*CS%Grad_L_Scale*VarMix%L2grad_u(I,j)*VarMix%UH_grad(I,j,k)
+!!        print*, "KH_u=", KH_u(I,j,k)
       enddo ; enddo ; enddo
     endif
   endif
@@ -424,7 +425,8 @@ subroutine thickness_diffuse(h, uhtr, vhtr, tv, dt, G, GV, US, MEKE, VarMix, CDp
     if (CS%Grad_L_Scale > 0.0) then  !< Gradient model
       !$OMP do
       do k=1,nz ; do J=js-1,je ; do i=is,ie
-        KH_v(i,J,k) = CS%Grad_L_Scale*VarMix%L2grad_v(i,J)*VarMix%VH_grad(i,J,k)
+        KH_v(i,J,k) = 1.0*CS%Grad_L_Scale*VarMix%L2grad_v(i,J)*VarMix%VH_grad(i,J,k)
+!!        print*, "KH_v=", KH_v(i,J,k)
       enddo ; enddo ; enddo
     endif
   endif