Merge pull request #559 from CliMA/al/push_start

push start time for IN05_17

Author: amylu00

papers/ice_nucleation_2024/AIDA_calibrations.jl

@@ -33,7 +33,7 @@ data_file_names = [
 ]
 batch_names = ["HOM", "IN0701", "IN0719", "ACI04_22", "EXP19", "EXP45", "DEP", "IMM"]
 start_time_list = [     # freezing onset
-    [Int32(150), Int32(180)],
+    [Int32(169), Int32(180)],
     [Int32(50)],
     [Int32(35)],
     [Int32(0)],
@@ -190,6 +190,7 @@ for (batch_index, batch_name) in enumerate(batch_names)
             exp_names[exp_index], Nₜ[exp_index],
             UKI_parcel,
             t_profile[exp_index], T_profile[exp_index], P_profile[exp_index], ICNC_profile[exp_index],
+            S_l_profile[exp_index],
         )
 
         ## Comparing AIDA data and calibrated parcel.

papers/ice_nucleation_2024/calibration.jl

@@ -25,7 +25,7 @@ function run_model(p_list, IN_mode, coefficients, FT, IC_list; calibration = fal
     for (exp_index, p) in enumerate(p_list)
         (; prescribed_thermodynamics, t_profile, T_profile, P_profile) = p
         (; const_dt, w, t_max, ips) = p
-        (; aerosol_act, aerosol, r_nuc) = p
+        (; aerosol_act, aerosol, r_nuc, A_aer) = p
         (; deposition_growth, condensation_growth) = p
         (; dep_nucleation, heterogeneous, homogeneous) = p
         (; liq_size_distribution, ice_size_distribution, aero_σ_g) = p
@@ -190,6 +190,7 @@ function run_model(p_list, IN_mode, coefficients, FT, IC_list; calibration = fal
             ice_size_distribution = ice_size_distribution,
             ips = ips,
             r_nuc = r_nuc,
+            A_aer = A_aer,
         )
 
         # solve ODE
@@ -270,8 +271,8 @@ function create_prior(FT, IN_mode, ; perfect_model = false, aerosol_type = nothi
                 error("Aerosol type not supported for ABIFM. Check create_prior function.")
             end
         elseif IN_mode == "ABHOM"
-            m_stats = [FT(260.927125), FT(70), FT(0), Inf]
-            c_stats = [FT(-68.553283), FT(20), -Inf, Inf]
+            m_stats = [FT(255.927125), FT(100), FT(0), Inf]
+            c_stats = [FT(-68.553283), FT(40 ), -Inf, Inf]
         end
     end
 

papers/ice_nucleation_2024/calibration_setup.jl

@@ -200,17 +200,17 @@ function AIDA_IN05_IC(FT, data_file)
     ϵₘ = R_d / R_v
 
     if data_file == "in05_17_aida.edf"
-        Nᵢ = FT(0.05 * 1e6)
-        Nₗ = FT(360 * 1e6) -  Nᵢ
+        Nᵢ = FT(0.647 * 1e6)
+        Nₗ = FT(360 * 1e6) - Nᵢ
         Nₐ = FT(0)
-        r₀ = FT(5.18056 / 2 * 1e-6)
-        p₀ = FT(894.409 * 1e2)
-        T₀ = FT(237.871)
+        r₀ = FT(5.68782 / 2 * 1e-6)
+        p₀ = FT(883.539 * 1e2)
+        T₀ = FT(237.502)
         qₗ = FT(Nₗ * 4 / 3 * FT(π) * r₀^3 * ρₗ / FT(1.2))  # 1.2 should be ρₐ
         qᵢ = FT(Nᵢ * 4 / 3 * FT(π) * r₀^3 * ρₗ / FT(1.2))
         m_l = Nₗ * ρₗ *  4 * π / 3 * r₀^3
         m_i = Nᵢ * ρᵢ *  4 * π / 3 * r₀^3
-        e = FT(30.2314)
+        e = FT(29.811)
         qᵥ = (e / R_v / T₀) / ((p₀ - e) / (R_d * T₀) + e / R_v / T₀ + m_l + m_i)
         q = TD.PhasePartition.(qᵥ + qₗ + qᵢ, qₗ, qᵢ)
         Rₐ = TD.gas_constant_air(tps, q)

papers/ice_nucleation_2024/plots.jl

@@ -109,6 +109,7 @@ function plot_calibrated_parcels(
     exp_name, Nₜ,
     UKI_parcel,
     t_profile, T_profile, P_profile, ICNC_profile,
+    Sₗ_profile,
 )
 
     calibrated_parcel_fig = MK.Figure(size = (1500, 800), fontsize = 20)
@@ -131,7 +132,7 @@ function plot_calibrated_parcels(
         color = :fuchsia,
         linestyle = :dash,
     )
-    # MK.lines!(ax_parcel_1, t_profile, S_l_profile, label = "chamber", color = :blue)
+    MK.lines!(ax_parcel_1, t_profile, Sₗ_profile, label = "chamber", color = :blue)
 
     MK.lines!(ax_parcel_2, UKI_parcel.t, UKI_parcel[5, :], color = :fuchsia)