diff --git a/src/P3Scheme.jl b/src/P3Scheme.jl index 4758c4c2c5..d6dacb9ab9 100644 --- a/src/P3Scheme.jl +++ b/src/P3Scheme.jl @@ -124,13 +124,13 @@ 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 == 0 - - return() + + return (; D_cr = 0, D_gr = 0, ρ_g = 0, ρ_d = 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 - + P3_problem(ρ_d) = ρ_d - ρ_d_helper( p3, @@ -206,20 +206,4 @@ function p3_mass( end end -#= """ - p3_mass(p3, D) - - - p3 - a struct with P3 scheme parameters - - D - maximum particle dimension - - Returns mass(D) refime specifically for no Rime (q_rim, F_r = 0) -""" -function p3_mass(p3::PSP3, D:: FT) where {FT <: Real} - if D_th_helper(p3) > D - return mass_s(D, p3.ρ_i) # small spherical ice - else - return mass_nl(p3, D) # large nonspherical unrimed ice - end -end =# - end