Revert to returning zero when dividing by zero

In get_field(::ClimaAtmosSimulation, ::Val{:diffuse_fraction}
CliMA · Feb 27, 2025 · f32cf97 · f32cf97
1 parent 67d8d6a
commit f32cf97
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/experiments/ClimaEarth/components/atmosphere/climaatmos.jl b/experiments/ClimaEarth/components/atmosphere/climaatmos.jl
@@ -188,7 +188,11 @@ function Interfacer.get_field(sim::ClimaAtmosSimulation, ::Val{:diffuse_fraction
         direct_flux_dn = radiation_model.face_sw_direct_flux_dn[1, :]
         FT = eltype(total_flux_dn)
         diffuse_fraction =
-            clamp.((total_flux_dn .- direct_flux_dn) ./ (max.(total_flux_dn, eps(FT))), zero(FT), one(FT))
+            clamp.(
+                ((x, y) -> y > zero(y) ? x / y : zero(y)).(total_flux_dn .- direct_flux_dn, total_flux_dn),
+                zero(FT),
+                one(FT),
+            )
     end
     return CC.Fields.array2field(diffuse_fraction, lowest_face_space)
 end