brb gonna fix mohler param

CliMA · Feb 1, 2024 · 89a4e61 · 89a4e61
1 parent aa4c9ea
commit 89a4e61
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 17 deletions.
diff --git a/parcel/Deposition_Nucleation.jl b/parcel/Deposition_Nucleation.jl
@@ -41,7 +41,7 @@ Sₗ = FT(e / eₛ)
 IC = [Sₗ, p₀, T₀, qᵥ, qₗ, qᵢ, Nₐ, Nₗ, Nᵢ, x_sulph]
 
 # Simulation parameters passed into ODE solver
-r_nuc = FT(5e-11)             # assumed size of nucleated particles
+r_nuc = FT(1.25e-7)             # assumed size of nucleated particles
 w = FT(3.5 * 1e-2)                                # updraft speed
 α_m = FT(0.5)                              # accomodation coefficient
 const_dt = FT(0.1)                           # model timestep

diff --git a/parcel/parcel.jl b/parcel/parcel.jl
@@ -82,23 +82,11 @@ function parcel_model(dY, Y, p, t)
             S_i,
             T,
         )
-        dN_act_dt_depo = max(FT(0), AF * N_aer - N_ice) / const_dt
+        #dN_act_dt_depo = max(FT(0), AF * N_aer - N_ice) / const_dt
+        # println("dS_l_dt = ", dY[1])
+        dN_act_dt_depo = max(FT(0), N_aer * aerosol.a_warm * (dY[1] * ξ))
         dqi_dt_new_depo = dN_act_dt_depo * 4 / 3 * FT(π) * r_nuc^3 * ρ_ice / ρ_air
-    end
-    if "DepositionNucleation" in ice_nucleation_modes
-        Δa_w = CMO.a_w_eT(tps, e, T) - CMO.a_w_ice(tps, T)
-        J_deposition = CMI_het.deposition_J(aerosol, Δa_w)
-        if "Monodisperse" in droplet_size_distribution
-            A_aer = 4 * FT(π) * r_nuc^2
-            dN_act_dt_depo = max(FT(0), J_deposition * N_aer * A_aer)
-            dqi_dt_new_depo = dN_act_dt_depo * 4 / 3 * FT(π) * r_nuc^3 * ρ_ice / ρ_air
-        end
-        if "Gamma" in droplet_size_distribution
-            A_aer = 4 * FT(π) * r_nuc^2
-            dN_act_dt_depo = max(FT(0), J_deposition * N_aer * A_aer)
-            dqi_dt_new_depo = dN_act_dt_depo * 4 / 3 * FT(π) * r_nuc^3 * ρ_ice / ρ_air
-        end
-        println("J_deposition = ", J_deposition)
+        # println("dN_act_dt_depo = ", dN_act_dt_depo)
     end
     if "ActivityBasedDeposition" in ice_nucleation_modes
         Δa_w = CMO.a_w_eT(tps, e, T) - CMO.a_w_ice(tps, T)
@@ -111,6 +99,7 @@ function parcel_model(dY, Y, p, t)
         println("N_aer = ", N_aer)
         println("A_aer = ", A_aer)
         println("J_deposition = ", J_deposition)
+        println("dN_act_dt_depo = ", dN_act_dt_depo)
     end
 
     dN_act_dt_immersion = FT(0)