diff --git a/parameterizations/lateral/MOM_thickness_diffuse.F90 b/parameterizations/lateral/MOM_thickness_diffuse.F90 index c3e251a7e9..bcba8f2021 100644 --- 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