add P3 ice nucleation modes

Author: haakon-e

src/BulkMicrophysicsTendencies.jl

@@ -32,6 +32,7 @@ import ..Microphysics1M as CM1
 import ..Microphysics2M as CM2
 import ..MicrophysicsNonEq as CMNonEq
 import ..P3Scheme as CMP3
+import ..HetIceNucleation as CM_HetIce
 import ...ThermodynamicsInterface as TDI
 
 export MicrophysicsScheme,
@@ -665,6 +666,9 @@ to be non-Nothing, eliminating runtime type checks and dynamic dispatch.
     vel = mp.ice.terminal_velocity
     pdf_c = mp.ice.cloud_pdf
     pdf_r = mp.ice.rain_pdf
+    heterogeneous = mp.ice.heterogeneous
+    deposition_condfreeze = mp.ice.deposition_condfreeze
+
 
     # Only compute ice processes if there is ice mass/number present
     if (q_ice > zero(q_ice) || n_ice > zero(n_ice))
@@ -684,23 +688,8 @@ to be non-Nothing, eliminating runtime type checks and dynamic dispatch.
         # Only compute if there is ice present
         if L_ice > zero(L_ice) && N_ice > zero(N_ice)
             coll = CMP3.bulk_liquid_ice_collision_sources(
-                p3,
-                logλ,
-                L_ice,
-                N_ice,
-                F_rim,
-                ρ_rim,
-                pdf_c,
-                pdf_r,
-                L_lcl,
-                N_lcl,
-                L_rai,
-                N_rai,
-                aps,
-                tps,
-                vel,
-                ρ,
-                T,
+                p3, logλ, L_ice, N_ice, F_rim, ρ_rim, pdf_c, pdf_r,
+                L_lcl, N_lcl, L_rai, N_rai, aps, tps, vel, ρ, T,
             )
 
             # Add collision tendencies
@@ -730,9 +719,66 @@ to be non-Nothing, eliminating runtime type checks and dynamic dispatch.
         end
     end
 
-    # TODO: When ice number concentration is tracked, add dn_ice_dt to return tuple:
-    # return (; dq_lcl_dt, dn_lcl_dt, dq_rai_dt, dn_rai_dt, dq_ice_dt, dn_ice_dt, dq_rim_dt, db_rim_dt)
-    return (; dq_lcl_dt, dn_lcl_dt, dq_rai_dt, dn_rai_dt, dq_ice_dt, dq_rim_dt, db_rim_dt)
+    # --- -------------------- ---
+    # --- Ice Nucleation Modes ---
+    # --- -------------------- ---
+
+    # --- Ice Nucleation (Deposition) ---
+    # Assume nucleated particles have diameter 1μm --> nucleated mass (per particle)
+    D_nuc_ice = 1e-6 # 1μm
+    m_nuc = p3.ρ_i * CO.volume_sphere_D(D_nuc_ice)  # [kg]
+
+    # TODO: The parameterixation should return a rate, `∂N/∂t`, not number changes `ΔN`
+    n_dep = CM_HetIce.P3_deposition_N_i(deposition_condfreeze, T) / ρ  # [particles / kg air]
+    m_dep = n_dep * m_nuc  # [kg ice / kg air]
+    dn_ice_dt += n_dep
+    dq_ice_dt += m_dep
+    dq_rim_dt += m_dep
+    db_rim_dt += m_dep / 900
+
+    # --- Heterogeneous Ice Nucleation (Immersion freezing) ---
+    # Assume mass loss is mean condensate mass
+    m_lcl = q_lcl / n_lcl  # mean liquid mass
+
+    # TODO: The parameterixation should return a rate, `∂N/∂t`, not number changes `ΔN`
+    inpc = CM_HetIce.INP_concentration_mean(heterogeneous, T) / ρ  # [particles / kg air]
+    n_het = max(FT(0), inpc - n_ice)  # [particles / kg air]
+    q_het = n_het * m_lcl  # [kg ice / kg air]
+
+    dn_ice_dt += n_het
+    dq_ice_dt += q_het
+    dq_rim_dt += q_het
+    db_rim_dt += q_het / 900  # = ρ^* TODO: make this a parameter?
+    dn_lcl_dt -= n_het
+    dq_lcl_dt -= q_het
+
+    # --- Cloud Droplet Condensation Freezing ---
+    # ref: `homogeneous_freezing` in `parcel/ParcelTendencies.jl`
+    # get mean diameter of cloud droplets, then convert to volume
+
+    # --- Ice Sublimation / Deposition ---
+    # Deposition/sublimation of cloud ice
+    ∂ₜq_ice_dep = CMNonEq.conv_q_vap_to_q_lcl_icl_MM2015(icl, tps, q_tot, q_lcl, q_ice, q_rai, zero(q_ice), ρ, T)
+    # No ice deposition above freezing (lack of INPs)
+    ∂ₜq_ice_dep = ifelse(T > tps.T_freeze, min(∂ₜq_ice_dep, zero(T)), ∂ₜq_ice_dep)
+    ∂ₜn_ice_dep = ifelse(∂ₜq_ice_dep < 0, (N_ice / q_ice) * ∂ₜq_ice_dep, zero(∂ₜq_ice_dep))
+    dq_ice_dt += ∂ₜq_ice_dep
+    dn_ice_dt += ∂ₜn_ice_dep
+
+    # --- Ice Self-collection (Aggregation) ---
+    # TODO: Implement P3 ice self-collection (aggregation)
+    # This process is currently missing in P3_processes.jl
+    # S_ice_agg = CMP3.ice_self_collection(p3, q_ice, n_ice, ...)
+    # dn_ice_dt -= S_ice_agg
+
+    # --- Rain Heterogeneous Freezing ---
+    # TODO: Implement heterogeneous freezing of rain
+    # This process is currently missing in P3_processes.jl
+    # S_rai_frz = ...
+    # dq_rai_dt -= S_rai_frz
+    # dq_ice_dt += S_rai_frz
+
+    return (; dq_lcl_dt, dn_lcl_dt, dq_rai_dt, dn_rai_dt, dq_ice_dt, dn_ice_dt, dq_rim_dt, db_rim_dt)
 end
 
 end # module BulkMicrophysicsTendencies

src/parameters/MicrophysicsP3.jl

@@ -338,7 +338,7 @@ function ParametersP3(toml_dict::CP.ParamDict; slope_law = :powerlaw)
     vent = VentilationFactor(toml_dict)
     ρ_rim_local = LocalRimeDensity(toml_dict)
     name_map = (;
-        :density_ice_water => :ρ_i,
+        :density_ice_water => :ρ_i,  # TODO: Check if this should be 900 or 916.7
         :density_liquid_water => :ρ_l,
         :temperature_water_freeze => :T_freeze,
         :P3_wet_growth_timescale => :τ_wet,