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