From d3c943aae224eea468cfaf1bf1027a318bc86383 Mon Sep 17 00:00:00 2001
From: amylu00 <alu3@caltech.edu>
Date: Thu, 5 Oct 2023 16:09:59 -0700
Subject: [PATCH] update a_w function parameters

---
 parcel/parcel.jl | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/parcel/parcel.jl b/parcel/parcel.jl
index 9cd4323fd4..d1862f7627 100644
--- a/parcel/parcel.jl
+++ b/parcel/parcel.jl
@@ -40,6 +40,7 @@ function parcel_model(dY, Y, p, t)
     grav = CMP.grav(prs)
     ρ_ice = CMP.ρ_cloud_ice(prs)
     ρ_liq = CMP.ρ_cloud_liq(prs)
+    H2SO4_prs = CMP.H2SO4SolutionParameters(FT)
 
     # Get thermodynamic parameters, phase partition and create thermo state.
     thermo_params = CMP.thermodynamics_params(prs)
@@ -87,8 +88,8 @@ function parcel_model(dY, Y, p, t)
     dqi_dt_new_immers = FT(0)
     if "ImmersionFreezing" in ice_nucleation_modes
         Δa_w = T > FT(185) && T < FT(235) ?
-            CMO.a_w_xT(prs, x_sulph, T) - CMO.a_w_ice(prs, T) :
-            CMO.a_w_eT(prs, e, T) - CMO.a_w_ice(prs, T)
+            CMO.a_w_xT(H2SO4_prs, thermo_params, x_sulph, T) - CMO.a_w_ice(prs, T) :
+            CMO.a_w_eT(thermo_params, e, T) - CMO.a_w_ice(prs, 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)