diff --git a/src/P3Scheme.jl b/src/P3Scheme.jl index e2ba94f78e..753588eb3e 100644 --- a/src/P3Scheme.jl +++ b/src/P3Scheme.jl @@ -125,7 +125,7 @@ function thresholds(p3::PSP3{FT}, ρ_r::FT, F_r::FT) where {FT} @assert F_r < FT(1) # ... and there must always be some unrimed part if F_r == FT(0) - return (; D_cr = 0, D_gr = 0, ρ_g = 0, ρ_d = 0) + return (; D_cr = FT(0), D_gr = FT(0), ρ_g = FT(0), ρ_d = FT(0)) else @assert ρ_r > FT(0) # rime density must be positive ... @assert ρ_r <= p3.ρ_l # ... and as a bulk ice density can't exceed the density of water