From 5f9a61f82c7cf0644549552f7a163adbe3c961ca Mon Sep 17 00:00:00 2001 From: Nicholas Hannah Date: Fri, 17 Jul 2015 11:06:42 -0700 Subject: [PATCH] Do full denominator check to avoid possible divide-by-zero. #194 --- src/parameterizations/vertical/MOM_set_diffusivity.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/parameterizations/vertical/MOM_set_diffusivity.F90 b/src/parameterizations/vertical/MOM_set_diffusivity.F90 index b1c6e9a227..a79d797594 100644 --- a/src/parameterizations/vertical/MOM_set_diffusivity.F90 +++ b/src/parameterizations/vertical/MOM_set_diffusivity.F90 @@ -1653,10 +1653,10 @@ subroutine add_LOTW_BBL_diffusivity(h, u, v, tv, fluxes, visc, j, N2_int, G, CS, ! Diffusivity using law of the wall, limited by rotation, at height z, in m2/s. ! This calculation is at the upper interface of the layer - if (D_minus_z /= 0.) then - Kd_wall = ( ( von_karm * ustar2 ) * ( z * D_minus_z ) )/( ustar_D + absf * ( z * D_minus_z ) ) - else + if ( ustar_D + absf * ( z * D_minus_z ) == 0.) then Kd_wall = 0. + else + Kd_wall = ( ( von_karm * ustar2 ) * ( z * D_minus_z ) )/( ustar_D + absf * ( z * D_minus_z ) ) endif ! TKE associated with Kd_wall, in m3 s-2.