Initialize Kd, Kd_int and Kv_slow using interior values specified by …

…user
NOAA-GFDL · Apr 19, 2018 · 3385857 · 3385857
1 parent 78656bb
commit 3385857
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/parameterizations/vertical/MOM_set_diffusivity.F90 b/src/parameterizations/vertical/MOM_set_diffusivity.F90
@@ -276,6 +276,12 @@ subroutine set_diffusivity(u, v, h, u_h, v_h, tv, fluxes, optics, visc, dt, &
     call MOM_error(FATAL, "set_diffusivity: visc%Kd_extra_T and "//&
          "visc%Kd_extra_S must be associated when DOUBLE_DIFFUSION is true.")
 
+  ! Set Kd, Kd_int and Kv_slow to constant values.
+  ! If nothing else is specified, this will be the value used.
+  Kd(:,:,:) = CS%Kd
+  Kd_int(:,:,:) = CS%Kd
+  visc%Kv_slow(:,:,:) = CS%Kv
+
   ! Set up arrays for diagnostics.
 
   if ((CS%id_N2 > 0) .or. (CS%id_N2_z > 0)) then