From 06d950c07c2e0c25a3d31c130126d1189b3ec59a Mon Sep 17 00:00:00 2001
From: dcherian <deepak@cherian.net>
Date: Wed, 22 Mar 2023 09:54:34 -0600
Subject: [PATCH] Fix combo of constant and lat dependent background
 diffusivity.

Add constant background viscosity KD to latitudinal dependent
formulations, if any.

Closes #238
---
 src/parameterizations/vertical/MOM_set_diffusivity.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/parameterizations/vertical/MOM_set_diffusivity.F90 b/src/parameterizations/vertical/MOM_set_diffusivity.F90
index 6d35616b3a..d0e934c840 100644
--- a/src/parameterizations/vertical/MOM_set_diffusivity.F90
+++ b/src/parameterizations/vertical/MOM_set_diffusivity.F90
@@ -538,7 +538,7 @@ subroutine set_diffusivity(u, v, h, u_h, v_h, tv, fluxes, optics, visc, dt, Kd_i
 
     ! Copy the 2-d slices into the 3-d array that is exported.
     do K=1,nz+1 ; do i=is,ie
-      Kd_int(i,j,K) = Kd_int_2d(i,K)
+      Kd_int(i,j,K) = Kd_int(i,j,K) + Kd_int_2d(i,K)
     enddo ; enddo
 
     if (CS%limit_dissipation) then