Added saturation ratio over ice to deposition

open-atmos · Sep 18, 2024 · 4d60bd2 · 4d60bd2
1 parent 0c879cb
commit 4d60bd2
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/PySDM/backends/impl_numba/methods/deposition_methods.py b/PySDM/backends/impl_numba/methods/deposition_methods.py
@@ -57,7 +57,9 @@ def body(
                         temperature, pressure
                     )
 
-                    saturation_ratio_ice = 1.1  # formulae.saturation_vapour_pressure.pvs_ice( temperature )
+                    saturation_ratio_ice = (
+                        ambient_humidity[cid] / ambient_water_activity[cid]
+                    )
 
                     dm_dt = (
                         4
@@ -86,6 +88,8 @@ def deposition(
         water_mass,
         ambient_temperature,
         ambient_total_pressure,
+        ambient_humidity,
+        ambient_water_activity,
         time_step,
         cell_id,
         reynolds_number,
@@ -95,6 +99,8 @@ def deposition(
             water_mass=water_mass.data,
             ambient_temperature=ambient_temperature.data,
             ambient_total_pressure=ambient_total_pressure.data,
+            ambient_humidity=ambient_humidity.data,
+            ambient_water_activity=ambient_water_activity.data,
             time_step=time_step,
             cell_id=cell_id.data,
             reynolds_number=reynolds_number.data,