diff --git a/src/parameterized_tendencies/microphysics/precipitation.jl b/src/parameterized_tendencies/microphysics/precipitation.jl index ac8fd950f38..a68e2716119 100644 --- a/src/parameterized_tendencies/microphysics/precipitation.jl +++ b/src/parameterized_tendencies/microphysics/precipitation.jl @@ -93,11 +93,10 @@ function precipitation_tendency!( # update precip in cache for coupler's use # 3d rain and snow + T_f = FT(273.15) @. ᶜT[colidx] = TD.air_temperature(thermo_params, ᶜts[colidx]) - @. ᶜ3d_rain[colidx] = - ifelse(ᶜT[colidx] >= FT(273.15), ᶜS_ρq_tot[colidx], FT(0)) - @. ᶜ3d_snow[colidx] = - ifelse(ᶜT[colidx] < FT(273.15), ᶜS_ρq_tot[colidx], FT(0)) + @. ᶜ3d_rain[colidx] = ifelse(ᶜT[colidx] >= T_f, ᶜS_ρq_tot[colidx], FT(0)) + @. ᶜ3d_snow[colidx] = ifelse(ᶜT[colidx] < T_f, ᶜS_ρq_tot[colidx], FT(0)) Operators.column_integral_definite!( col_integrated_rain[colidx], ᶜ3d_rain[colidx],