PySDM 2.53__tar.gz → 2.54__tar.gz

{PySDM-2.53 → PySDM-2.54}/.github/workflows/tests+artifacts+pypi.yml

@@ -132,7 +132,7 @@ jobs:
             python-version: "3.8"
       fail-fast: false
     runs-on: ${{ matrix.platform }}
-    timeout-minutes: 40
+    timeout-minutes: 45
     steps:
       - uses: actions/checkout@v4.1.1
         with:

{PySDM-2.53 → PySDM-2.54}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.53
+Version: 2.54
 Summary: Pythonic particle-based (super-droplet) warm-rain/aqueous-chemistry cloud microphysics package with box, parcel & 1D/2D prescribed-flow examples in Python, Julia and Matlab
 Home-page: https://github.com/open-atmos/PySDM
 Author: https://github.com/open-atmos/PySDM/graphs/contributors
@@ -842,8 +842,10 @@ Description: # PySDM
           https://github.com/darothen/superdroplet
         - NTLP (FORTRAN)   
           https://github.com/Folca/NTLP/blob/SuperDroplet/les.F
-        - CLEO (C++)
+        - CLEO (C++)    
           https://yoctoyotta1024.github.io/CLEO/
+        - droplets.jl (Julia)   
+          https://github.com/emmacware/droplets.jl
         
         ### non-SDM probabilistic particle-based coagulation solvers
         
@@ -856,6 +858,10 @@ Description: # PySDM
         - PyBox: https://github.com/loftytopping/PyBox
         - py-cloud-parcel-model: https://github.com/emmasimp/py-cloud-parcel-model
         
+        ### non-Python cloud microphysics open-source software
+        
+        - CloudMicrophysics.jl: https://github.com/CliMA/CloudMicrophysics.jl
+        
 Keywords: physics-simulation,monte-carlo-simulation,gpu-computing,atmospheric-modelling,particle-system,numba,thrust,nvrtc,pint,atmospheric-physics
 Platform: UNKNOWN
 Classifier: Development Status :: 4 - Beta

{PySDM-2.53 → PySDM-2.54}/PySDM/backends/impl_numba/methods/terminal_velocity_methods.py

@@ -30,22 +30,17 @@ class TerminalVelocityMethods(BackendMethods):
 
     @cached_property
     def terminal_velocity_body(self):
+        v_term = self.formulae.terminal_velocity.v_term
+
         @numba.njit(**{**conf.JIT_FLAGS, "fastmath": self.formulae.fastmath})
-        def terminal_velocity_body(*, values, radius, k1, k2, k3, r1, r2):
+        def terminal_velocity_body(*, values, radius):
             for i in numba.prange(len(values)):  # pylint: disable=not-an-iterable
-                if radius[i] < r1:
-                    values[i] = k1 * radius[i] ** 2
-                elif radius[i] < r2:
-                    values[i] = k2 * radius[i]
-                else:
-                    values[i] = k3 * radius[i] ** (1 / 2)
+                values[i] = v_term(radius[i])
 
         return terminal_velocity_body
 
-    def terminal_velocity(self, *, values, radius, k1, k2, k3, r1, r2):
-        self.terminal_velocity_body(
-            values=values, radius=radius, k1=k1, k2=k2, k3=k3, r1=r1, r2=r2
-        )
+    def terminal_velocity(self, *, values, radius):
+        self.terminal_velocity_body(values=values, radius=radius)
 
     @cached_property
     def power_series_body(self):

{PySDM-2.53 → PySDM-2.54}/PySDM/backends/impl_thrust_rtc/methods/physics_methods.py

@@ -68,28 +68,6 @@ class PhysicsMethods(ThrustRTCBackendMethods):
             ),
         )
 
-    @cached_property
-    def __terminal_velocity_body(self):
-        return trtc.For(
-            ("values", "radius", "k1", "k2", "k3", "r1", "r2"),
-            "i",
-            """
-            if (radius[i] < r1) {
-                values[i] = k1 * radius[i] * radius[i];
-            }
-            else {
-                if (radius[i] < r2) {
-                    values[i] = k2 * radius[i];
-                }
-                else {
-                    values[i] = k3 * pow(radius[i], (real_type)(.5));
-                }
-            }
-            """.replace(
-                "real_type", self._get_c_type()
-            ),
-        )
-
     @cached_property
     def __volume_of_mass_body(self):
         return trtc.For(
@@ -145,17 +123,6 @@ class PhysicsMethods(ThrustRTCBackendMethods):
             ),
         )
 
-    @nice_thrust(**NICE_THRUST_FLAGS)
-    def terminal_velocity(self, *, values, radius, k1, k2, k3, r1, r2):
-        k1 = self._get_floating_point(k1)
-        k2 = self._get_floating_point(k2)
-        k3 = self._get_floating_point(k3)
-        r1 = self._get_floating_point(r1)
-        r2 = self._get_floating_point(r2)
-        self.__terminal_velocity_body.launch_n(
-            values.size(), [values, radius, k1, k2, k3, r1, r2]
-        )
-
     @nice_thrust(**NICE_THRUST_FLAGS)
     def explicit_euler(self, y, dt, dy_dt):
         dt = self._get_floating_point(dt)

{PySDM-2.53 → PySDM-2.54}/PySDM/backends/impl_thrust_rtc/methods/terminal_velocity_methods.py

@@ -155,3 +155,21 @@ class TerminalVelocityMethods(ThrustRTCBackendMethods):
         self.__power_series_body.launch_n(
             values.size(), (values, radius, num_terms, prefactors, powers)
         )
+
+    @cached_property
+    def __terminal_velocity_body(self):
+        return trtc.For(
+            param_names=("values", "radius"),
+            name_iter="i",
+            body=f"""
+            values[i] = {self.formulae.terminal_velocity.v_term.c_inline(
+                radius="radius[i]"
+            )};
+            """.replace(
+                "real_type", self._get_c_type()
+            ),
+        )
+
+    @nice_thrust(**NICE_THRUST_FLAGS)
+    def terminal_velocity(self, *, values, radius):
+        self.__terminal_velocity_body.launch_n(n=values.size(), args=[values, radius])

PySDM-2.54/PySDM/dynamics/terminal_velocity/rogers_and_yau.py

@@ -0,0 +1,14 @@
+"""
+Rogers & Yau, equations: (8.5), (8.6), (8.8)
+"""
+
+
+class RogersYau:  # pylint: disable=too-few-public-methods
+    def __init__(self, particulator):
+        self.particulator = particulator
+
+    def __call__(self, output, radius):
+        self.particulator.backend.terminal_velocity(
+            values=output.data,
+            radius=radius.data,
+        )

{PySDM-2.53 → PySDM-2.54}/PySDM/formulae.py

@@ -48,6 +48,8 @@ class Formulae:  # pylint: disable=too-few-public-methods,too-many-instance-attr
         isotope_equilibrium_fractionation_factors: str = "Null",
         isotope_meteoric_water_line_excess: str = "Null",
         isotope_ratio_evolution: str = "Null",
+        isotope_diffusivity_ratios: str = "Null",
+        isotope_relaxation_timescale: str = "Null",
         optical_albedo: str = "Null",
         optical_depth: str = "Null",
         particle_shape_and_density: str = "LiquidSpheres",
@@ -79,8 +81,11 @@ class Formulae:  # pylint: disable=too-few-public-methods,too-many-instance-attr
         )
         self.isotope_meteoric_water_line_excess = isotope_meteoric_water_line_excess
         self.isotope_ratio_evolution = isotope_ratio_evolution
+        self.isotope_diffusivity_ratios = isotope_diffusivity_ratios
+        self.isotope_relaxation_timescale = isotope_relaxation_timescale
         self.particle_shape_and_density = particle_shape_and_density
         self.air_dynamic_viscosity = air_dynamic_viscosity
+        self.terminal_velocity = terminal_velocity
 
         self._components = tuple(
             i
@@ -298,6 +303,8 @@ def _pick(value: str, choices: dict, constants: namedtuple):
         for name, cls in choices.items():
             if name == value:
                 obj = cls(constants)
+                break
+            name = value
     else:
         parent_classes = []
         for name in value.split("+"):
@@ -317,7 +324,7 @@ def _pick(value: str, choices: dict, constants: namedtuple):
 
     if obj is None:
         raise ValueError(
-            f"Unknown setting: '{name}'; choices are: {tuple(choices.keys())}"
+            f"Unknown setting: {name}; choices are: {', '.join(choices.keys())}"
         )
 
     obj.__name__ = value  # pylint: disable=attribute-defined-outside-init

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/__init__.py

@@ -31,6 +31,8 @@ from . import (
     isotope_equilibrium_fractionation_factors,
     isotope_meteoric_water_line_excess,
     isotope_ratio_evolution,
+    isotope_diffusivity_ratios,
+    isotope_relaxation_timescale,
     latent_heat,
     optical_albedo,
     optical_depth,
@@ -42,5 +44,6 @@ from . import (
     trivia,
     ventilation,
     air_dynamic_viscosity,
+    terminal_velocity,
 )
 from .constants import convert_to, in_unit, si

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/constants.py

@@ -43,6 +43,7 @@ ONE_THIRD = 1 / 3
 ONE_HALF = 1 / 2
 TWO_THIRDS = 2 / 3
 ONE_AND_A_HALF = 3 / 2
+TWO_AND_A_HALF = 5 / 2
 NaN = np.nan
 
 default_random_seed = (

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/constants_defaults.py

@@ -310,15 +310,15 @@ CRAIG_1961_INTERCEPT_COEFF = 10 * PER_MILLE
 """ [Bohren 1987](https://doi.org/10.1119/1.15109) """
 asymmetry_g = 0.85  # forward scattering from cloud droplets
 
-""" TODO #1266 """
-diffussion_thermics_D_G11_A = 1e-5 * si.m**2 / si.s
-diffussion_thermics_D_G11_B = 0.015 / si.K
-diffussion_thermics_D_G11_C = -1.9
+""" [Grabowski et al. 2011](https://doi.org/10.1016/j.atmosres.2010.10.020) """
+diffusion_thermics_D_G11_A = 1e-5 * si.m**2 / si.s
+diffusion_thermics_D_G11_B = 0.015 / si.K
+diffusion_thermics_D_G11_C = -1.9
 
-diffussion_thermics_K_G11_A = 1.5e-11 * si.W / si.m / si.K**4
-diffussion_thermics_K_G11_B = -4.8e-8 * si.W / si.m / si.K**3
-diffussion_thermics_K_G11_C = 1e-4 * si.W / si.m / si.K**2
-diffussion_thermics_K_G11_D = -3.9e-4 * si.W / si.m / si.K
+diffusion_thermics_K_G11_A = 1.5e-11 * si.W / si.m / si.K**4
+diffusion_thermics_K_G11_B = -4.8e-8 * si.W / si.m / si.K**3
+diffusion_thermics_K_G11_C = 1e-4 * si.W / si.m / si.K**2
+diffusion_thermics_K_G11_D = -3.9e-4 * si.W / si.m / si.K
 
 """
 [Pruppacher & Rasmussen (1979))](https://doi.org/10.1175/1520-0469(1979)036<1255:AWTIOT>2.0.CO;2)
@@ -343,11 +343,32 @@ FROESSLING_1938_A = 1
 FROESSLING_1938_B = 0.276
 
 
+""" fit coefficients from [Hellmann & Harvey 2020](https://doi.org/10.1029/2020GL089999) """
+HELLMANN_HARVEY_T_UNIT = 100 * si.K
+HELLMANN_HARVEY_EQ6_COEFF0 = 0.98258
+HELLMANN_HARVEY_EQ6_COEFF1 = -0.02546
+HELLMANN_HARVEY_EQ6_COEFF2 = 0.02421
+HELLMANN_HARVEY_EQ7_COEFF0 = 0.98284
+HELLMANN_HARVEY_EQ7_COEFF1 = 0.003517
+HELLMANN_HARVEY_EQ7_COEFF2 = -0.001996
+HELLMANN_HARVEY_EQ8_COEFF0 = 0.96671
+HELLMANN_HARVEY_EQ8_COEFF1 = 0.007406
+HELLMANN_HARVEY_EQ8_COEFF2 = -0.004861
+
+
+""" terminal velocity formulation from Rogers & Yau (equations: 8.5, 8.6, 8.8) """
+ROGERS_YAU_TERM_VEL_SMALL_K = 1.19e6 / si.cm / si.s
+ROGERS_YAU_TERM_VEL_MEDIUM_K = 8e3 / si.s
+ROGERS_YAU_TERM_VEL_LARGE_K = 2.01e3 * si.cm**0.5 / si.s
+ROGERS_YAU_TERM_VEL_SMALL_R_LIMIT = 35 * si.um
+ROGERS_YAU_TERM_VEL_MEDIUM_R_LIMIT = 600 * si.um
+
+
 def compute_derived_values(c: dict):
     """
     computes derived quantities such as molar mass ratios, etc.
 
-    water molar mass is computed from from molecular masses and VSMOW isotope abundances
+    water molar mass is computed from molecular masses and VSMOW isotope abundances
     (and neglecting molecular binding energies)
     for discussion, see:
     - caption of Table 2.1 in [Gat 2010](https://doi.org/10.1142/p027)

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/diffusion_thermics/grabowski_et_al_2011.py

@@ -10,16 +10,16 @@ class GrabowskiEtAl2011:
     @staticmethod
     def D(const, T, p):  # pylint: disable=unused-argument
         """eq (10)"""
-        return const.diffussion_thermics_D_G11_A * (
-            const.diffussion_thermics_D_G11_B * T + const.diffussion_thermics_D_G11_C
+        return const.diffusion_thermics_D_G11_A * (
+            const.diffusion_thermics_D_G11_B * T + const.diffusion_thermics_D_G11_C
         )
 
     @staticmethod
     def K(const, T, p):  # pylint: disable=unused-argument
         """eq (12)"""
         return (
-            const.diffussion_thermics_K_G11_A * T**3
-            + const.diffussion_thermics_K_G11_B * T**2
-            + const.diffussion_thermics_K_G11_C * T
-            + const.diffussion_thermics_K_G11_D
+            const.diffusion_thermics_K_G11_A * T**3
+            + const.diffusion_thermics_K_G11_B * T**2
+            + const.diffusion_thermics_K_G11_C * T
+            + const.diffusion_thermics_K_G11_D
         )

PySDM-2.54/PySDM/physics/isotope_diffusivity_ratios/__init__.py

@@ -0,0 +1,6 @@
+""" heavy-to-light diffusivity ratios for water molecules """
+
+from PySDM.impl.null_physics_class import Null
+from .grahams_law import GrahamsLaw
+from .stewart_1975 import Stewart1975
+from .hellmann_and_harvey_2020 import HellmannAndHarvey2020

PySDM-2.54/PySDM/physics/isotope_diffusivity_ratios/grahams_law.py

@@ -0,0 +1,14 @@
+""" [Graham's law](https://en.wikipedia.org/wiki/Graham%27s_law)
+    see also eq. (21) in [Horita et al. 2008](https://doi.org/10.1080/10256010801887174)
+"""
+
+
+class GrahamsLaw:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    def ratio_2H(const, temperature):  # pylint: disable=unused-argument
+        return (
+            (2 * const.M_1H + const.M_16O) / (const.M_2H + const.M_1H + const.M_16O)
+        ) ** const.ONE_HALF

PySDM-2.54/PySDM/physics/isotope_diffusivity_ratios/hellmann_and_harvey_2020.py

@@ -0,0 +1,39 @@
+"""
+Ratios of diffusivity in air of heavy vs. light isotope using fits provided in
+[Hellmann & Harvey 2020](https://doi.org/10.1029/2020GL089999)
+"""
+
+
+class HellmannAndHarvey2020:
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    def ratio_2H(const, temperature):
+        return (
+            const.HELLMANN_HARVEY_EQ6_COEFF0
+            + const.HELLMANN_HARVEY_EQ6_COEFF1
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT)
+            + const.HELLMANN_HARVEY_EQ6_COEFF2
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT) ** (const.TWO_AND_A_HALF)
+        )
+
+    @staticmethod
+    def ratio_17O(const, temperature):
+        return (
+            const.HELLMANN_HARVEY_EQ7_COEFF0
+            + const.HELLMANN_HARVEY_EQ7_COEFF1
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT) ** (const.ONE_HALF)
+            + const.HELLMANN_HARVEY_EQ7_COEFF2
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT) ** (const.TWO_AND_A_HALF)
+        )
+
+    @staticmethod
+    def ratio_18O(const, temperature):
+        return (
+            const.HELLMANN_HARVEY_EQ8_COEFF0
+            + const.HELLMANN_HARVEY_EQ8_COEFF1
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT) ** (const.ONE_HALF)
+            + const.HELLMANN_HARVEY_EQ8_COEFF2
+            / (temperature / const.HELLMANN_HARVEY_T_UNIT) ** (const.THREE)
+        )

PySDM-2.54/PySDM/physics/isotope_diffusivity_ratios/stewart_1975.py

@@ -0,0 +1,40 @@
+"""
+Temperature-independent ratio of vapour diffusivity in air for heavy vs. light isotopes
+assuming same collision diameters for different isotopic water molecules, see:
+eq. (8) in [Stewart 1975](https://doi.org/10.1029/JC080i009p01133),
+eq. (1) in [Merlivat 1978](https://doi.org/10.1063/1.436884),
+eq. (6) in [Cappa et al. 2003](https://doi.org/10.1029/2003JD003597),
+eq. (22) in [Horita et al. 2008](https://doi.org/10.1080/10256010801887174),
+eq. (3) in [Hellmann and Harvey 2020](https://doi.org/10.1029/2020GL089999).
+
+All functions return constants, so there is a potential overhead in computing them on each call,
+but this variant is provided for historical reference only, hence leaving like that.
+"""
+
+
+class Stewart1975:
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    def ratio_2H(const, temperature):  # pylint: disable=unused-argument
+        return (
+            (
+                (2 * const.M_1H + const.M_16O)
+                * (const.Md + const.M_2H + const.M_1H + const.M_16O)
+            )
+            / (
+                (const.M_2H + const.M_1H + const.M_16O)
+                * (const.Md + (2 * const.M_1H + const.M_16O))
+            )
+        ) ** const.ONE_HALF
+
+    @staticmethod
+    def ratio_18O(const, temperature):  # pylint: disable=unused-argument
+        return (
+            ((2 * const.M_1H + const.M_16O) * (const.Md + 2 * const.M_1H + const.M_18O))
+            / (
+                (2 * const.M_1H + const.M_18O)
+                * (const.Md + (2 * const.M_1H + const.M_16O))
+            )
+        ) ** const.ONE_HALF

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/isotope_equilibrium_fractionation_factors/__init__.py

@@ -1,12 +1,12 @@
 """ factors describing partitioning of water isotopologues under thermodynamic
 phase equilibrium, defined as ratios alpha of isotopic ratios R """
 
+from PySDM.impl.null_physics_class import Null
 from .barkan_and_luz_2005 import BarkanAndLuz2005
 from .horita_and_wesolowski_1994 import HoritaAndWesolowski1994
 from .majoube_1970 import Majoube1970
 from .majoube_1971 import Majoube1971
 from .merlivat_and_nief_1967 import MerlivatAndNief1967
-from .null import Null
 from .van_hook_1968 import VanHook1968
 from .lamb_et_al_2017 import LambEtAl2017
 from .ellehoj_et_al_2013 import EllehojEtAl2013

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/isotope_meteoric_water_line_excess/__init__.py

@@ -1,5 +1,5 @@
 """ definitions of meteoric water line excess parameters """
 
+from PySDM.impl.null_physics_class import Null
 from .barkan_and_luz_2007 import BarkanAndLuz2007
 from .dansgaard_1964 import Dansgaard1964
-from .null import Null

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/isotope_meteoric_water_line_excess/barkan_and_luz_2007.py

@@ -15,3 +15,10 @@ class BarkanAndLuz2007:  # pylint: disable=too-few-public-methods
         return np.log(
             delta_17O + 1
         ) - const.BARKAN_AND_LUZ_2007_EXCESS_18O_COEFF * np.log(delta_18O + 1)
+
+    @staticmethod
+    def d17O_of_d18O(const, delta_18O):
+        return (
+            np.exp(const.BARKAN_AND_LUZ_2007_EXCESS_18O_COEFF * np.log(delta_18O + 1))
+            - 1
+        )

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/isotope_ratio_evolution/__init__.py

@@ -1,5 +1,5 @@
 """ isotope-ratio evolution formulae (incl. Rayleigh distillation) """
 
+from PySDM.impl.null_physics_class import Null
 from .merlivat_and_jouzel_1979 import MerlivatAndJouzel1979
-from .null import Null
 from .rayleigh_distillation import RayleighDistillation

PySDM-2.54/PySDM/physics/isotope_relaxation_timescale/__init__.py

@@ -0,0 +1,6 @@
+"""
+isotope relaxation timescale formulae
+"""
+
+from PySDM.impl.null_physics_class import Null
+from .miyake_et_al_1968 import MiyakeEtAl1968

PySDM-2.54/PySDM/physics/isotope_relaxation_timescale/miyake_et_al_1968.py

@@ -0,0 +1,13 @@
+""" eq. 28 in [Miyake et al. 1968](https://doi.org/10.2467/mripapers1950.19.2_243) """
+
+
+class MiyakeEtAl1968:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    # pylint: disable=too-many-arguments
+    def tau(const, e_s, D, M, vent_coeff, radius, alpha, temperature):
+        return (radius**2 * alpha * const.rho_w * const.R_str * temperature) / (
+            3 * e_s * D * M * vent_coeff
+        )

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/optical_albedo/__init__.py

@@ -2,5 +2,5 @@
 alternative formulations of cloud albedo
 """
 
+from PySDM.impl.null_physics_class import Null
 from .bohren1987 import Bohren1987
-from .null import Null

{PySDM-2.53 → PySDM-2.54}/PySDM/physics/optical_depth/__init__.py

@@ -2,5 +2,5 @@
 alternative formulations of optical depth
 """
 
+from PySDM.impl.null_physics_class import Null
 from .stephens_1978 import Stephens1978
-from .null import Null

PySDM-2.54/PySDM/physics/terminal_velocity/__init__.py

@@ -0,0 +1,5 @@
+""" terminal velocity formulae """
+
+from .rogers_yau import RogersYau
+from .gunn_kinzer_1949 import GunnKinzer1949
+from .power_series import PowerSeries

PySDM-2.54/PySDM/physics/terminal_velocity/gunn_kinzer_1949.py

@@ -0,0 +1,6 @@
+""" [Gunn & Kinzer 1949](https://doi.org/10.1175/1520-0469(1949)006%3C0243:TTVOFF%3E2.0.CO;2) """
+
+
+class GunnKinzer1949:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass

PySDM-2.54/PySDM/physics/terminal_velocity/power_series.py

@@ -0,0 +1,6 @@
+""" power-series formulation """
+
+
+class PowerSeries:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass

PySDM-2.54/PySDM/physics/terminal_velocity/rogers_yau.py

@@ -0,0 +1,22 @@
+"""
+Rogers & Yau "A short course in cloud physics" (equations: 8.5, 8.6, 8.8)
+"""
+
+import numpy as np
+
+
+class RogersYau:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    def v_term(const, radius):
+        return np.where(
+            radius < const.ROGERS_YAU_TERM_VEL_SMALL_R_LIMIT,
+            const.ROGERS_YAU_TERM_VEL_SMALL_K * radius**const.TWO,
+            np.where(
+                radius < const.ROGERS_YAU_TERM_VEL_MEDIUM_R_LIMIT,
+                const.ROGERS_YAU_TERM_VEL_MEDIUM_K * radius,
+                const.ROGERS_YAU_TERM_VEL_LARGE_K * radius**const.ONE_HALF,
+            ),
+        )

{PySDM-2.53 → PySDM-2.54}/PySDM.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.53
+Version: 2.54
 Summary: Pythonic particle-based (super-droplet) warm-rain/aqueous-chemistry cloud microphysics package with box, parcel & 1D/2D prescribed-flow examples in Python, Julia and Matlab
 Home-page: https://github.com/open-atmos/PySDM
 Author: https://github.com/open-atmos/PySDM/graphs/contributors
@@ -842,8 +842,10 @@ Description: # PySDM
           https://github.com/darothen/superdroplet
         - NTLP (FORTRAN)   
           https://github.com/Folca/NTLP/blob/SuperDroplet/les.F
-        - CLEO (C++)
+        - CLEO (C++)    
           https://yoctoyotta1024.github.io/CLEO/
+        - droplets.jl (Julia)   
+          https://github.com/emmacware/droplets.jl
         
         ### non-SDM probabilistic particle-based coagulation solvers
         
@@ -856,6 +858,10 @@ Description: # PySDM
         - PyBox: https://github.com/loftytopping/PyBox
         - py-cloud-parcel-model: https://github.com/emmasimp/py-cloud-parcel-model
         
+        ### non-Python cloud microphysics open-source software
+        
+        - CloudMicrophysics.jl: https://github.com/CliMA/CloudMicrophysics.jl
+        
 Keywords: physics-simulation,monte-carlo-simulation,gpu-computing,atmospheric-modelling,particle-system,numba,thrust,nvrtc,pint,atmospheric-physics
 Platform: UNKNOWN
 Classifier: Development Status :: 4 - Beta