diff --git a/parcel/parcel.jl b/parcel/parcel.jl
index 7e12f6aedb..c3f0c67882 100644
--- a/parcel/parcel.jl
+++ b/parcel/parcel.jl
@@ -93,16 +93,16 @@ function parcel_model(dY, Y, p, t)
         J_immersion = CMI_het.ABIFM_J(prs, aerosol_type, Δa_w)
         if "Monodisperse" in droplet_size_distribution && "ImmersionFreezing" in ice_nucleation_modes
             r_l = cbrt(q_liq / N_liq / (4 / 3 * π) / ρ_liq * ρ_air)
-            A_aer = 4 * π * r_l^2 * N_liq
-            dN_act_dt_immersion = max(FT(0), J_immersion * A_aer)
+            A_aer = 4 * π * r_l^2
+            dN_act_dt_immersion = max(FT(0), J_immersion * N_liq * A_aer)
             dqi_dt_new_immers = dN_act_dt_immersion * 4 / 3 * π * r_l^3 * ρ_ice / ρ_air
         end
         if "Gamma" in droplet_size_distribution && "ImmersionFreezing" in ice_nucleation_modes
             λ = cbrt(32 * π * N_liq / q_liq * ρ_liq / ρ_air)
             #A = N_liq* λ^2
             r_l = 2 / λ
-            A_aer = 4 * π * r_l^2 * N_liq
-            dN_act_dt_immersion = max(FT(0), J_immersion * A_aer)
+            A_aer = 4 * π * r_l^2
+            dN_act_dt_immersion = max(FT(0), J_immersion * N_liq * A_aer)
             dqi_dt_new_immers = dN_act_dt_immersion * 4 / 3 * π * r_l^3 * ρ_ice / ρ_air
         end
     end