Checked values for howell factor and added it to deposition growth term

open-atmos · Dec 11, 2024 · 0e839c4 · 0e839c4
1 parent 211f1f0
commit 0e839c4
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/PySDM/backends/impl_numba/methods/deposition_methods.py b/PySDM/backends/impl_numba/methods/deposition_methods.py
@@ -68,8 +68,7 @@ def body(
                     thermal_conductivity = formulae.diffusion_ice_kinetics__K(Ka_const,  radius, lambdaK, temperature, rho)  
 
 
-                    howell_factor = 1. / ((latent_heat_sub / Rv / temperature - 1.) * latent_heat_sub * diffusion_coefficient / temperature / thermal_conductivity +  Rv * temperature  /  pvs_ice)
-                    # print( f" {howell_factor=}, {latent_heat_sub= }, {Rv=}, {Dv_const=}, {diffusion_coefficient=}, {Ka_const=}, {thermal_conductivity=}, {pvs_ice=}")
+                    howell_factor = 1. / ( (latent_heat_sub / Rv / temperature - 1.) * latent_heat_sub * diffusion_coefficient / temperature / thermal_conductivity +  Rv * temperature  /  pvs_ice)
 
                     saturation_ratio_ice = (
                         ambient_humidity[cid] / ambient_water_activity[cid]
@@ -86,6 +85,7 @@ def body(
                         * np.pi
                         * capacity
                         * diffusion_coefficient
+                        * howell_factor
                         * (saturation_ratio_ice - 1)
                     ) * rho_vs_ice