Apply suggestions

Author: charleskawczynski

src/parameterized_tendencies/les_sgs_models/smagorinsky_lilly.jl

@@ -112,11 +112,10 @@ function horizontal_smagorinsky_lilly_tendency!(Yₜ, Y, p, t, ::SmagorinskyLill
             ᶜρχ = getproperty(Y.c, ρχ_name)
             ᶜρχₜ = getproperty(Yₜ.c, ρχ_name)
             ᶜχ = @. lazy(specific(ᶜρχ, Y.c.ρ))
-            ᶜρχₜ_diffusion =
-                @. p.scratch.ᶜtemp_scalar = wdivₕ(Y.c.ρ * ᶜD_smag * gradₕ(ᶜχ))
+            ᶜρχₜ_diffusion = @. lazy(wdivₕ(Y.c.ρ * ᶜD_smag * gradₕ(ᶜχ)))
             @. ᶜρχₜ += ᶜρχₜ_diffusion
             # Rain and snow does not affect the mass
-            if ρχ_name ∉ (:ρq_rai, :ρq_sno)
+            if ρχ_name == :ρq_tot
                 @. Yₜ.c.ρ += ᶜρχₜ_diffusion
             end
         end
@@ -128,7 +127,7 @@ import UnrolledUtilities as UU
 
 function vertical_smagorinsky_lilly_tendency!(Yₜ, Y, p, t, ::SmagorinskyLilly)
     FT = eltype(Y)
-    (; sfc_temp_C3, ᶠtemp_scalar, ᶜtemp_scalar) = p.scratch
+    (; sfc_temp_C3, ᶠtemp_scalar) = p.scratch
     (; ᶜτ_smag, ᶠτ_smag, ᶠD_smag, ᶜspecific, ᶜh_tot, sfc_conditions) =
         p.precomputed
     (; ρ_flux_uₕ, ρ_flux_h_tot) = sfc_conditions
@@ -172,7 +171,7 @@ function vertical_smagorinsky_lilly_tendency!(Yₜ, Y, p, t, ::SmagorinskyLilly)
                 bottom = Operators.SetValue(C3(FT(0))),
             )
 
-            ᶜ∇ᵥρD∇χₜ = @. ᶜtemp_scalar = ᶜdivᵥ_ρχ(-(ᶠρ * ᶠD_smag * ᶠgradᵥ(ᶜχ)))
+            ᶜ∇ᵥρD∇χₜ = @. lazy(ᶜdivᵥ_ρχ(-(ᶠρ * ᶠD_smag * ᶠgradᵥ(ᶜχ))))
             @. ᶜρχₜ -= ᶜ∇ᵥρD∇χₜ
             # Rain and snow does not affect the mass
             if ρχ_name == :ρq_tot