revisited Fig 2 from Gedzelman & Arnold 1994; plot for AMS Madison Summit 2026

PySDM/attributes/isotopes/bolin_numbers.py

@@ -25,7 +25,7 @@ def __init__(self, builder, *, heavy_isotope: str):
             Heavy isotopologue identifier (entry of ``HEAVY_ISOTOPES``).
         """
         self.moles_heavy = builder.get_attribute(f"moles_{heavy_isotope}")
-        self.moles_light = builder.get_attribute("moles light water")
+        self.moles_light = builder.get_attribute("moles light water")  # TODO #1787
         self.cell_id = builder.get_attribute("cell id")
         self.molality_in_dry_air = builder.particulator.environment[
             f"molality {heavy_isotope} in dry air"

PySDM/backends/impl_numba/methods/isotope_methods.py

@@ -64,19 +64,16 @@ def body(
             moles_heavy_molecule,
             bolin_number,
             molality_in_dry_air,
-        ):
+        ):  # pylint: disable=too-many-locals
             for sd_id in range(multiplicity.shape[0]):
+                Bo = bolin_number[sd_id]
                 mass_ratio_heavy_to_total = (
                     moles_heavy_molecule[sd_id] * molar_mass_heavy_molecule
                 ) / signed_water_mass[sd_id]
-                if bolin_number[sd_id] == 0:
+                if Bo == 0:
                     dm_heavy = 0
                 else:
-                    dm_heavy = (
-                        dm_total[sd_id]
-                        / bolin_number[sd_id]
-                        * mass_ratio_heavy_to_total
-                    )
+                    dm_heavy = dm_total[sd_id] / Bo * mass_ratio_heavy_to_total
                 dn_heavy_molecule = dm_heavy / molar_mass_heavy_molecule
                 moles_heavy_molecule[sd_id] += dn_heavy_molecule
                 mass_of_dry_air = (

PySDM/physics/constants_defaults.py

@@ -103,7 +103,7 @@
 [Kaercher & Lohmann 2002](https://doi.org/10.1029/2001JD000470) """
 HAC_ice = 1.0
 """ thermal accommodation coefficient for vapour deposition as recommended in
-[Pruppacher & Klett 2010](https://doi.org/10.1007/978-0-306-48100-0) """
+[Pruppacher & Klett](https://doi.org/10.1007/978-0-306-48100-0) """
 
 C_cunn = 0.7
 """ Cunningham correction factor as used in
@@ -222,12 +222,12 @@
 """ 〃 """
 
 dv_pk05 = 0.0 * si.metres
-""" Delta v for diffusivity in [Pruppacher & Klett 2010](https://doi.org/10.1007/978-0-306-48100-0)
+""" Delta v for diffusivity in [Pruppacher & Klett](https://doi.org/10.1007/978-0-306-48100-0)
 eq. 13-14 """
 
 lmbd_w_0 = 6.6e-8 * si.metre
 """ Mean free path of water molecules as in table 13.1 in
-[Pruppacher & Klett 2010](https://doi.org/10.1007/978-0-306-48100-0) """
+[Pruppacher & Klett](https://doi.org/10.1007/978-0-306-48100-0) """
 
 d_l19_a = 0.211e-4 * si.metre**2 / si.second
 """ [Seinfeld & Pandis](https://archive.org/details/0237-pdf-atmospheric-chemistry-and-physics-2nd-ed-j.-seinfeld-s.-pandis-wiley-2006-ww)

PySDM/physics/diffusion_kinetics/grabowski_et_al_2011.py

@@ -2,7 +2,7 @@
 as in [Grabowski et al. (2011)](https://doi.org/10.1016/j.atmosres.2010.10.020)
 """
 
-from .pruppacher_and_klett_2010 import PruppacherKlett
+from .pruppacher_and_klett_2005 import PruppacherKlett
 
 
 class GrabowskiEtAl2011(PruppacherKlett):

PySDM/physics/diffusion_kinetics/lowe_et_al_2019.py

@@ -1,10 +1,10 @@
 """
 as in [Lowe et al. 2019](https://doi.org/10.1038/s41467-019-12982-0)
-uses eq. 13-14 in [Pruppacher & Klett 2010](https://doi.org/10.1007/978-0-306-48100-0)
-with Delta v = 0 and no corrections for thermal conductivity
+uses eq. 13-14 in [Pruppacher & Klett](https://doi.org/10.1007/978-0-306-48100-0) with Delta v = 0
+and no corrections for thermal conductivity
 """
 
-from PySDM.physics.diffusion_kinetics.pruppacher_and_klett_2010 import PruppacherKlett
+from PySDM.physics.diffusion_kinetics.pruppacher_and_klett_2005 import PruppacherKlett
 
 
 class LoweEtAl2019(PruppacherKlett):

PySDM/physics/diffusion_kinetics/pruppacher_and_klett_2010.py → PySDM/physics/diffusion_kinetics/pruppacher_and_klett_2005.py

@@ -1,5 +1,5 @@
 """
-as in [Pruppacher & Klett 2010](https://doi.org/10.1007/978-0-306-48100-0) (eq. 13-14)
+as in Pruppacher and Klett 2005 (eq. 13-14)
 with reference to [Okuyama and Zung 1967](https://doi.org/10.1063/1.1840906)
 """
 

PySDM/physics/isotope_relaxation_timescale/zaba_et_al.py

@@ -22,6 +22,7 @@ def tau(
 
     @staticmethod
     def bolin_number(
+        const,
         D_ratio_heavy_to_light,
         alpha,
         D_light,
@@ -42,8 +43,10 @@ def bolin_number(
             * (1 - relative_humidity)
             / D_ratio_heavy_to_light
             / (
-                (1 + rho_v * D_light * Fk)
+                (1 + rho_v * D_light * Fk / const.rho_w)
                 * relative_humidity
                 * (1 - alpha * R_vap / R_liq)
+                - relative_humidity
+                + 1
             )
         )

docs/bibliography.json

@@ -251,7 +251,7 @@
   },
   "https://doi.org/10.1063/1.1840906": {
     "usages": [
-      "PySDM/physics/diffusion_kinetics/pruppacher_and_klett_2010.py"
+      "PySDM/physics/diffusion_kinetics/pruppacher_and_klett_2005.py"
     ],
     "label": "Okuyama & Zung 1967 (J. Chem. Phys. 1967)",
     "title": "Evaporation—Condensation Coefficient for Small Droplets"
@@ -789,8 +789,7 @@
   "https://doi.org/10.1007/978-0-306-48100-0": {
     "usages": [
       "PySDM/physics/constants_defaults.py",
-      "PySDM/physics/diffusion_kinetics/lowe_et_al_2019.py",
-      "PySDM/physics/diffusion_kinetics/pruppacher_and_klett_2010.py"
+      "PySDM/physics/diffusion_kinetics/lowe_et_al_2019.py"
     ],
     "title": "Microphysics of Clouds and Precipitation",
     "label": "Pruppacher & Klett 2010 (Springer Dordrecht)"