Merge pull request #2905 from CliMA/zs/hotfix

fix dividing by 0 in percentile_bounds_mean_norm

Author: szy21

src/prognostic_equations/edmfx_boundary_condition.jl

@@ -56,8 +56,9 @@ function percentile_bounds_mean_norm(
     gauss_int(x) = -exp(-x * x / 2) / sqrt(2 * pi)
     xp_low = Distributions.quantile(Distributions.Normal(), low_percentile)
     xp_high = Distributions.quantile(Distributions.Normal(), high_percentile)
-    return (gauss_int(xp_high) - gauss_int(xp_low)) / (
+    return (gauss_int(xp_high) - gauss_int(xp_low)) / max(
         Distributions.cdf(Distributions.Normal(), xp_high) -
-        Distributions.cdf(Distributions.Normal(), xp_low)
+        Distributions.cdf(Distributions.Normal(), xp_low),
+        eps(FT),
     )
 end