Avoid NaNs when computing stratification parameter

NOAA-GFDL · May 18, 2018 · d5ce7a4 · d5ce7a4
1 parent 05daede
commit d5ce7a4
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/parameterizations/vertical/MOM_cvmix_ddiff.F90 b/src/parameterizations/vertical/MOM_cvmix_ddiff.F90
@@ -250,6 +250,8 @@ subroutine compute_ddiff_coeffs(h, tv, G, GV, j, Kd_T, Kd_S, CS)
     if (CS%id_R_rho > 0.0)  then
       do k=1,G%ke
         CS%R_rho(i,j,k) = alpha_dT(k)/beta_dS(k)
+        ! avoid NaN's
+        if(CS%R_rho(i,j,k) /= CS%R_rho(i,j,k)) CS%R_rho(i,j,k) = 0.0
       enddo
     endif