Ude condensate in eaporation/sublimation

Author: trontrytel

src/MicrophysicsNonEq.jl

@@ -193,9 +193,11 @@ function conv_q_vap_to_q_lcl_icl_MM2015(
     dqsl_dT = dqcld_dT(qᵥ_sat_liq, Lᵥ, Rᵥ, T)
     Γₗ = gamma_helper(Lᵥ, cₚ_air, dqsl_dT)
 
-    # compute the tendency
-    tendency = (qᵥ - qᵥ_sat_liq) / (τ_relax * Γₗ)
+    sat_excess = qᵥ - qᵥ_sat_liq
+    timescale = τ_relax * Γₗ
 
+    # compute the tendency
+    tendency = ifelse(sat_excess < 0, q_lcl / timescale, sat_excess / timescale)
     return ifelse(limit_MM2015_sinks(tendency, q_lcl), zero(tendency), tendency)
 end
 function conv_q_vap_to_q_lcl_icl_MM2015(
@@ -219,8 +221,11 @@ function conv_q_vap_to_q_lcl_icl_MM2015(
     dqsi_dT = dqcld_dT(qᵥ_sat_ice, Lₛ, Rᵥ, T)
     Γᵢ = gamma_helper(Lₛ, cₚ_air, dqsi_dT)
 
+    sat_excess = qᵥ - qᵥ_sat_ice
+    timescale = τ_relax * Γᵢ
+
     # compute the tendency
-    tendency = (qᵥ - qᵥ_sat_ice) / (τ_relax * Γᵢ)
+    tendency = ifelse(sat_excess < 0, q_icl / timescale, sat_excess / timescale)
 
     limiter = limit_MM2015_sinks(tendency, q_icl) || INP_limiter(tendency, tps, T)
     return ifelse(limiter, zero(tendency), tendency)