From 3385857b6fe19e7c7e35a0370914cdc339c131fa Mon Sep 17 00:00:00 2001 From: Gustavo Marques Date: Thu, 19 Apr 2018 08:55:22 -0600 Subject: [PATCH] Initialize Kd, Kd_int and Kv_slow using interior values specified by user --- src/parameterizations/vertical/MOM_set_diffusivity.F90 | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/src/parameterizations/vertical/MOM_set_diffusivity.F90 b/src/parameterizations/vertical/MOM_set_diffusivity.F90 index 7f192d65d9..b9905977d5 100644 --- 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