From 58125633d7253aca8e22fe83e56eaa6722f7eee1 Mon Sep 17 00:00:00 2001 From: Irena Vankova Date: Tue, 8 Oct 2024 08:31:51 -0600 Subject: [PATCH] replaceconfig_flux_attenuation_coefficient_runoff with config_flux_attenuation_coefficient_subglacial_runoff in subglacial kpp loop --- components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F b/components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F index 1777e1f6d1e9..0334667f2f10 100644 --- a/components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F +++ b/components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F @@ -3574,8 +3574,8 @@ subroutine ocn_compute_KPP_input_fields(statePool, forcingPool, & do k = maxLevelCell(iCell), minLevelCell(iCell), -1 zBot = zTop - layerThickness(k,iCell) if (k == minLevelCell(iCell)) then - transmissionCoeffTop = exp( max(zTop / config_flux_attenuation_coefficient_runoff, -100.0_RKIND) ) - transmissionCoeffBot = exp( max(zBot / config_flux_attenuation_coefficient_runoff, -100.0_RKIND) ) + transmissionCoeffTop = exp( max(zTop / config_flux_attenuation_coefficient_subglacial_runoff, -100.0_RKIND) ) + transmissionCoeffBot = exp( max(zBot / config_flux_attenuation_coefficient_subglacial_runoff, -100.0_RKIND) ) fracAbsorbedSubglacialRunoff = transmissionCoeffTop - transmissionCoeffBot end if zTop = zBot