update a_w function parameters

CliMA · Oct 5, 2023 · d3c943a · d3c943a
1 parent 17bdda9
commit d3c943a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git 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)