PyPI - PySDM-examples - Versions diffs - 2.39__py3-none-any.whl → 2.41__py3-none-any.whl - Mend

PySDM-examples 2.39py3-none-any.whl → 2.41py3-none-any.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (37) hide show

PySDM_examples/Abdul_Razzak_Ghan_2000/aerosol.py CHANGED Viewed

@@ -16,6 +16,7 @@ CONSTANTS_ARG = {
 class AerosolARG(DryAerosolMixture):
     def __init__(
         self,
+        water_molar_volume: float,
         M2_sol: float = 0,
         M2_N: float = 100 / si.cm**3,
         M2_rad: float = 50 * si.nm,
@@ -35,13 +36,19 @@ class AerosolARG(DryAerosolMixture):
         )
         self.modes = (
             {
-                "kappa": self.kappa({"(NH4)2SO4": 1.0, "insoluble": 0.0}),
+                "kappa": self.kappa(
+                    mass_fractions={"(NH4)2SO4": 1.0, "insoluble": 0.0},
+                    water_molar_volume=water_molar_volume,
+                ),
                 "spectrum": spectra.Lognormal(
                     norm_factor=100.0 / si.cm**3, m_mode=50.0 * si.nm, s_geom=2.0
                 ),
             },
             {
-                "kappa": self.kappa({"(NH4)2SO4": M2_sol, "insoluble": (1 - M2_sol)}),
+                "kappa": self.kappa(
+                    mass_fractions={"(NH4)2SO4": M2_sol, "insoluble": (1 - M2_sol)},
+                    water_molar_volume=water_molar_volume,
+                ),
                 "spectrum": spectra.Lognormal(
                     norm_factor=M2_N, m_mode=M2_rad, s_geom=2.0
                 ),
@@ -51,7 +58,7 @@ class AerosolARG(DryAerosolMixture):
 @strict
 class AerosolWhitby(DryAerosolMixture):
-    def __init__(self):
+    def __init__(self, water_molar_volume: float):
         nuclei = {"(NH4)2SO4": 1.0}
         accum = {"(NH4)2SO4": 1.0}
         coarse = {"(NH4)2SO4": 1.0}
@@ -69,19 +76,25 @@ class AerosolWhitby(DryAerosolMixture):
         )
         self.modes = (
             {
-                "kappa": self.kappa(nuclei),
+                "kappa": self.kappa(
+                    mass_fractions=nuclei, water_molar_volume=water_molar_volume
+                ),
                 "spectrum": spectra.Lognormal(
                     norm_factor=1000.0 / si.cm**3, m_mode=0.008 * si.um, s_geom=1.6
                 ),
             },
             {
-                "kappa": self.kappa(accum),
+                "kappa": self.kappa(
+                    mass_fractions=accum, water_molar_volume=water_molar_volume
+                ),
                 "spectrum": spectra.Lognormal(
                     norm_factor=800 / si.cm**3, m_mode=0.034 * si.um, s_geom=2.1
                 ),
             },
             {
-                "kappa": self.kappa(coarse),
+                "kappa": self.kappa(
+                    mass_fractions=coarse, water_molar_volume=water_molar_volume
+                ),
                 "spectrum": spectra.Lognormal(
                     norm_factor=0.72 / si.cm**3, m_mode=0.46 * si.um, s_geom=2.2
                 ),

PySDM_examples/Abdul_Razzak_Ghan_2000/run_ARG_parcel.py CHANGED Viewed

@@ -6,7 +6,6 @@ from PySDM_examples.Abdul_Razzak_Ghan_2000.aerosol import CONSTANTS_ARG, Aerosol
 from PySDM import Builder, Formulae
 from PySDM import products as PySDM_products
 from PySDM.backends import CPU
-from PySDM.backends.impl_numba.test_helpers import scipy_ode_condensation_solver
 from PySDM.dynamics import AmbientThermodynamics, Condensation
 from PySDM.environments import Parcel
 from PySDM.initialisation import equilibrate_wet_radii
@@ -47,10 +46,12 @@ def run_parcel(
         T0=T0,
     )
-    aerosol = AerosolARG(M2_sol=sol2, M2_N=N2, M2_rad=rad2)
+    aerosol = AerosolARG(
+        M2_sol=sol2, M2_N=N2, M2_rad=rad2, water_molar_volume=const.Mv / const.rho_w
+    )
     n_sd = n_sd_per_mode * len(aerosol.modes)
-    builder = Builder(backend=CPU(), n_sd=n_sd, environment=env)
+    builder = Builder(backend=CPU(formulae), n_sd=n_sd, environment=env)
     builder.add_dynamic(AmbientThermodynamics())
     builder.add_dynamic(Condensation())
     builder.request_attribute("critical supersaturation")
@@ -79,7 +80,6 @@ def run_parcel(
     attributes["volume"] = builder.formulae.trivia.volume(radius=r_wet)
     particulator = builder.build(attributes, products=products)
-    scipy_ode_condensation_solver.patch_particulator(particulator)
     output = {product.name: [] for product in particulator.products.values()}
     output_attributes = {

PySDM_examples/Alpert_and_Knopf_2016/__init__.py CHANGED Viewed

@@ -2,6 +2,7 @@
 box-environment example based on
 [Alpert & Knopf 2016 (Atmos. Chem. Phys. 16)](https://doi.org/10.5194/acp-16-2083-2016)
 """
 # pylint: disable=invalid-name
 from .simulation import Simulation, simulation
 from .table_1 import Table1

PySDM_examples/Arabas_and_Shima_2017/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 condensation-evaportaion parcel example based on
 [Arabas and Shima 2017 (Nonlin. Processes Geophys. 24)](https://doi.org/10.5194/npg-24-535-2017)
 """
 # pylint: disable=invalid-name

PySDM_examples/Arabas_et_al_2015/settings.py CHANGED Viewed

@@ -28,10 +28,7 @@ class Settings(StratoCumulus):
     def __init__(
         self,
         formulae=None,
-        rhod_w_max: float = 0.6
-        * si.metres
-        / si.seconds
-        * (si.kilogram / si.metre**3),
+        rhod_w_max: float = 0.6 * si.metres / si.seconds * (si.kilogram / si.metre**3),
     ):
         super().__init__(formulae or Formulae(), rhod_w_max=rhod_w_max)

PySDM_examples/Arabas_et_al_2023/__init__.py CHANGED Viewed

@@ -2,6 +2,7 @@
 box-model and 2D prescribed-flow immersion-freezing examples based on
 [Arabas et al. 2023](https://arxiv.org/abs/2308.05015)
 """
 from .make_particulator import make_particulator
 from .plots import make_freezing_spec_plot, make_pdf_plot, make_temperature_plot
 from .run_simulation import run_simulation

PySDM_examples/Arabas_et_al_2023/plots.py CHANGED Viewed

@@ -28,9 +28,11 @@ def make_temperature_plot(data):
         xy2.plot(
             datum["t"],
             np.asarray(datum["qi"]) / qi_unit,  # marker='.',
-            label=f"Monte-Carlo ({labels[v['singular']]})"
-            if not plotted[v["singular"]]
-            else "",
+            label=(
+                f"Monte-Carlo ({labels[v['singular']]})"
+                if not plotted[v["singular"]]
+                else ""
+            ),
             color=colors[v["singular"]],
         )
         plotted[v["singular"]] = True

PySDM_examples/Bartman_et_al_2021/label.py CHANGED Viewed

@@ -1,11 +1,11 @@
 def label(settings):
     lbl = str(
         {
-            k.replace("condensation_", ""): f"{v:.1e}"
-            if isinstance(v, float)
-            else str(v).zfill(2)
-            if isinstance(v, int)
-            else v
+            k.replace("condensation_", ""): (
+                f"{v:.1e}"
+                if isinstance(v, float)
+                else str(v).zfill(2) if isinstance(v, int) else v
+            )
             for k, v in settings.items()
         }
     )

PySDM_examples/Berry_1967/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 box-model coalescence-only example based on
 [Berry 1967 (J. Atmos. Sci. 24)](https://doi.org/10.1175/1520-0469(1967)024%3C0688:CDGBC%3E2.0.CO;2)
 """
 # pylint: disable=invalid-name

PySDM_examples/Bieli_et_al_2022/__init__.py CHANGED Viewed

@@ -1,4 +1,5 @@
 """
 collision-only box-model example from [Bieli et al. 2022](https://doi.org/10.1029/2022MS002994)
 """
 # pylint: disable=invalid-name

PySDM_examples/Bolot_et_al_2013/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 figures from Bolot et al. 2013 (Atmos. Chem. Phys. 13)
 https://doi.org/10.5194/acp-13-7903-2013
 """
 # pylint: disable=invalid-name

PySDM_examples/Ervens_and_Feingold_2012/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 aerosol distribution from
 [Ervens and Feingold 2012 (Atmos. Chem. Phys. 12)](https://doi.org/10.5194/acp-12-5807-2012)
 """
 # pylint: disable=invalid-name

PySDM_examples/Ervens_and_Feingold_2012/settings.py CHANGED Viewed

@@ -7,8 +7,6 @@ def sampled_ccn_diameter_number_concentration_spectrum(
     n_sd: int = 11, size_range: tuple = (0.02 * si.um, 2 * si.um)
 ):
     return Logarithmic(
-        spectrum=Lognormal(
-            s_geom=1.4, m_mode=0.04 * si.um, norm_factor=100 / si.cm**3
-        ),
+        spectrum=Lognormal(s_geom=1.4, m_mode=0.04 * si.um, norm_factor=100 / si.cm**3),
         size_range=size_range,
     ).sample(n_sd)

PySDM_examples/Kreidenweis_et_al_2003/settings.py CHANGED Viewed

@@ -70,11 +70,13 @@ class Settings:
         self.starting_amounts = {
             "moles_"
-            + k: self.formulae.trivia.volume(self.r_dry)
-            * self.DRY_RHO
-            / self.dry_molar_mass
-            if k in ("N_mIII", "S_VI")
-            else np.zeros(self.n_sd)
+            + k: (
+                self.formulae.trivia.volume(self.r_dry)
+                * self.DRY_RHO
+                / self.dry_molar_mass
+                if k in ("N_mIII", "S_VI")
+                else np.zeros(self.n_sd)
+            )
             for k in AQUEOUS_COMPOUNDS
         }

PySDM_examples/Lowe_et_al_2019/aerosol.py CHANGED Viewed

@@ -8,7 +8,9 @@ from PySDM.physics import si
 @strict
 class AerosolMarine(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.2, Acc_N2: float = 134):
+    def __init__(
+        self, water_molar_volume: float, Forg: float = 0.2, Acc_N2: float = 134
+    ):
         Aitken = {
             "palmitic": Forg,
             "(NH4)2SO4": (1 - Forg),
@@ -48,7 +50,7 @@ class AerosolMarine(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Aitken),
-                "kappa": self.kappa(Aitken),
+                "kappa": self.kappa(Aitken, water_molar_volume=water_molar_volume),
                 "nu_org": self.nu_org(Aitken),
                 "spectrum": spectra.Lognormal(
                     norm_factor=226 / si.cm**3, m_mode=19.6 * si.nm, s_geom=1.71
@@ -56,7 +58,9 @@ class AerosolMarine(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3, m_mode=69.5 * si.nm, s_geom=1.7
@@ -69,7 +73,9 @@ class AerosolMarine(DryAerosolMixture):
 @strict
 class AerosolBoreal(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.668, Acc_N2: float = 540):
+    def __init__(
+        self, water_molar_volume: float, Forg: float = 0.668, Acc_N2: float = 540
+    ):
         # note: SOA1 or SOA2 unclear from the paper
         Aitken = {
             "SOA1": Forg,
@@ -116,7 +122,7 @@ class AerosolBoreal(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Aitken),
-                "kappa": self.kappa(Aitken),
+                "kappa": self.kappa(Aitken, water_molar_volume=water_molar_volume),
                 "nu_org": self.nu_org(Aitken),
                 "spectrum": spectra.Lognormal(
                     norm_factor=1110 / si.cm**3, m_mode=22.7 * si.nm, s_geom=1.75
@@ -124,7 +130,9 @@ class AerosolBoreal(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3,
@@ -139,7 +147,9 @@ class AerosolBoreal(DryAerosolMixture):
 @strict
 class AerosolNascent(DryAerosolMixture):
-    def __init__(self, Acc_Forg: float = 0.3, Acc_N2: float = 30):
+    def __init__(
+        self, water_molar_volume: float, Acc_Forg: float = 0.3, Acc_N2: float = 30
+    ):
         Ultrafine = {
             "SOA1": 0.52,
             "SOA2": 0,
@@ -178,7 +188,7 @@ class AerosolNascent(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Ultrafine),
-                "kappa": self.kappa(Ultrafine),
+                "kappa": self.kappa(Ultrafine, water_molar_volume=water_molar_volume),
                 "nu_org": self.nu_org(Ultrafine),
                 "spectrum": spectra.Lognormal(
                     norm_factor=2000 / si.cm**3, m_mode=11.5 * si.nm, s_geom=1.71
@@ -186,7 +196,9 @@ class AerosolNascent(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3, m_mode=100 * si.nm, s_geom=1.70

PySDM_examples/Lowe_et_al_2019/aerosol_code.py CHANGED Viewed

@@ -19,7 +19,9 @@ class AerosolMarine(DryAerosolMixture):
     # NMODE                   = [1 1];
     # DENSITY                 = [1.841, 1.78, 1.77, 0.852, 1.5, 2., 2.65, 2.165]; %[gcm-3]
-    def __init__(self, Forg: float = 0.2, Acc_N2: float = 137):
+    def __init__(
+        self, water_molar_volume: float, Forg: float = 0.2, Acc_N2: float = 137
+    ):
         Aitken = {
             "palmitic": Forg,
             "(NH4)2SO4": (1 - Forg),
@@ -59,7 +61,9 @@ class AerosolMarine(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Aitken),
-                "kappa": self.kappa(Aitken),
+                "kappa": self.kappa(
+                    mass_fractions=Aitken, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Aitken),
                 "spectrum": spectra.Lognormal(
                     norm_factor=223 / si.cm**3, m_mode=19.6 * si.nm, s_geom=1.68
@@ -67,7 +71,9 @@ class AerosolMarine(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    mass_fractions=Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3, m_mode=69.5 * si.nm, s_geom=1.68
@@ -100,9 +106,11 @@ class AerosolBoreal(DryAerosolMixture):
     # DENSITY                     = [1.841, 1.78, 1.77, 1.5, 1.5, 2., 2.65, 2.165]; %[gcm-3]
     # DENSITY                     = [1.72, 1.78, 1.77, 1.2, 1.4, 2., 2.65, 2.165]; %[gcm-3]
-    def __init__(self, Forg: float = 0.668, Acc_N2: float = 540):
-        # TODO #604: SOA1 or SOA2 unclear from the paper
-        # TODO #604: CAN'T FIND WHERE NH4NO3 PROPERTIES ARE DEFINED IN ICPM
+    def __init__(
+        self, water_molar_volume: float, Forg: float = 0.668, Acc_N2: float = 540
+    ):
+        # TODO #1247: SOA1 or SOA2 unclear from the paper
+        # TODO #1247: CAN'T FIND WHERE NH4NO3 PROPERTIES ARE DEFINED IN ICPM
         INORG_MASS_RATIO = 0.1515 / 0.1559
         Aitken = {
             "SOA1": Forg,
@@ -129,8 +137,8 @@ class AerosolBoreal(DryAerosolMixture):
                 * si.gram
                 / si.mole,
                 "NH4NO3": Substance.from_formula("NH4NO3").mass * si.gram / si.mole,
-                "SOA1": 190 * si.g / si.mole,  # TODO #604: 190 OR 200?
-                "SOA2": 368.4 * si.g / si.mole,  # TODO #604: 368.4 OR 200?
+                "SOA1": 190 * si.g / si.mole,  # TODO #1247: 190 OR 200?
+                "SOA2": 368.4 * si.g / si.mole,  # TODO #1247: 368.4 OR 200?
             },
             densities={
                 "SOA1": 1.2 * si.g / si.cm**3,
@@ -149,7 +157,9 @@ class AerosolBoreal(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Aitken),
-                "kappa": self.kappa(Aitken),
+                "kappa": self.kappa(
+                    mass_fractions=Aitken, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Aitken),
                 "spectrum": spectra.Lognormal(
                     norm_factor=1110 / si.cm**3, m_mode=22.65 * si.nm, s_geom=1.75
@@ -157,7 +167,9 @@ class AerosolBoreal(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    mass_fractions=Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3,
@@ -184,8 +196,10 @@ class AerosolNascent(DryAerosolMixture):
     # DENSITY                 = [1.841, 1.78, 1.77, 1.2, 1.24, 2., 2.65, 2.165]; %[gcm-3]
-    def __init__(self, Acc_Forg: float = 0.3, Acc_N2: float = 30):
-        # TODO #604: CAN'T FIND WHEN PHI IS MULTIPLIED FOR KÖHLER B IN ICPM CODE
+    def __init__(
+        self, water_molar_volume: float, Acc_Forg: float = 0.3, Acc_N2: float = 30
+    ):
+        # TODO #1247: CAN'T FIND WHEN PHI IS MULTIPLIED FOR KÖHLER B IN ICPM CODE
         Ultrafine = {
             "SOA1": 0.52,
             "SOA2": 0,
@@ -203,8 +217,8 @@ class AerosolNascent(DryAerosolMixture):
                 "(NH4)2SO4": 3,
             },
             molar_masses={
-                "SOA1": 190 * si.g / si.mole,  # TODO #604: 190 OR 200?
-                "SOA2": 368.4 * si.g / si.mole,  # TODO #604: 368.4 OR 200?
+                "SOA1": 190 * si.g / si.mole,  # TODO #1247: 190 OR 200?
+                "SOA2": 368.4 * si.g / si.mole,  # TODO #1247: 368.4 OR 200?
                 "(NH4)2SO4": Substance.from_formula("(NH4)2SO4").mass
                 * si.gram
                 / si.mole,
@@ -224,7 +238,9 @@ class AerosolNascent(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(Ultrafine),
-                "kappa": self.kappa(Ultrafine),
+                "kappa": self.kappa(
+                    mass_fractions=Ultrafine, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Ultrafine),
                 "spectrum": spectra.Lognormal(
                     norm_factor=2000 / si.cm**3, m_mode=11.5 * si.nm, s_geom=1.71
@@ -232,7 +248,9 @@ class AerosolNascent(DryAerosolMixture):
             },
             {
                 "f_org": 1 - self.f_soluble_volume(Accumulation),
-                "kappa": self.kappa(Accumulation),
+                "kappa": self.kappa(
+                    mass_fractions=Accumulation, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(Accumulation),
                 "spectrum": spectra.Lognormal(
                     norm_factor=Acc_N2 / si.cm**3, m_mode=100 * si.nm, s_geom=1.703

PySDM_examples/Merlivat_and_Nief_1967/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 figure from Merlivat and Nief 1967 (Tellus)
 https://doi.org/10.3402/tellusa.v19i1.9756
 """
 # pylint: disable=invalid-name

PySDM_examples/Niedermeier_et_al_2014/__init__.py CHANGED Viewed

@@ -1,3 +1,4 @@
 """ based on [Niedermeier et al. 2014](https://doi.org/10.1002/2013GL058684) """
 from .settings import Settings
 from .simulation import Simulation

PySDM_examples/Niedermeier_et_al_2014/simulation.py CHANGED Viewed

@@ -68,9 +68,11 @@ class Simulation(BasicSimulation):
                     data[:-1],
                     np.full(
                         settings.in_sampling_n,
-                        data[-1]
-                        if attribute != "multiplicity"
-                        else data[-1] / settings.in_sampling_n,
+                        (
+                            data[-1]
+                            if attribute != "multiplicity"
+                            else data[-1] / settings.in_sampling_n
+                        ),
                     ),
                 )
             )

PySDM_examples/Pierchala_et_al_2022/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 based on Pierchala et al. 2022 (Geochim. Cosmochim. Acta)
 https://doi.org/10.1016/j.gca.2022.01.020
 """
 # pylint: disable=invalid-name

PySDM_examples/Pyrcel/profile_plotter.py CHANGED Viewed

@@ -57,11 +57,11 @@ class ProfilePlotter:
                 settings.formulae.trivia.radius(volume=np.asarray(volume)) / si.um,
                 output["products"]["z"],
                 color="magenta" if drop_id < settings.n_sd_per_mode[0] else "blue",
-                label="mode 1"
-                if drop_id == 0
-                else "mode 2"
-                if drop_id == settings.n_sd_per_mode[0]
-                else "",
+                label=(
+                    "mode 1"
+                    if drop_id == 0
+                    else "mode 2" if drop_id == settings.n_sd_per_mode[0] else ""
+                ),
             )
         axR.legend(loc="upper right")
         axR.set_xlabel("Droplet radius [μm]")

PySDM_examples/Shipway_and_Hill_2012/settings.py CHANGED Viewed

@@ -1,9 +1,9 @@
 from typing import Iterable
 import numpy as np
+from numdifftools import Derivative
 from scipy.integrate import solve_ivp
 from scipy.interpolate import interp1d
-from scipy.misc import derivative
 from PySDM import Formulae
 from PySDM.dynamics import condensation
@@ -60,8 +60,8 @@ class Settings:
         self.t_max = t_max
         t_1 = 600 * si.s
-        self.rho_times_w = (
-            lambda t: rho_times_w_1 * np.sin(np.pi * t / t_1) if t < t_1 else 0
+        self.rho_times_w = lambda t: (
+            rho_times_w_1 * np.sin(np.pi * t / t_1) if t < t_1 else 0
         )
         apprx_w1 = rho_times_w_1 / self.formulae.constants.rho_STP
         self.particle_reservoir_depth = (
@@ -106,9 +106,9 @@ class Settings:
             water_vapour_mixing_ratio = self.water_vapour_mixing_ratio(
                 z_above_reservoir
             )
-            d_water_vapour_mixing_ratio__dz = derivative(
-                self.water_vapour_mixing_ratio, z_above_reservoir
-            )
+            d_water_vapour_mixing_ratio__dz = Derivative(
+                self.water_vapour_mixing_ratio
+            )(z_above_reservoir)
             T = self.formulae.state_variable_triplet.T(
                 rhod[0], self.thd(z_above_reservoir)
             )

PySDM_examples/Singer_Ward/aerosol.py CHANGED Viewed

@@ -8,7 +8,7 @@ from PySDM.physics import si
 @strict
 class AerosolBetaCaryophylleneDark(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.8, N: float = 400):
+    def __init__(self, water_molar_volume: float, Forg: float = 0.8, N: float = 400):
         mode = {
             "(NH4)2SO4": (1 - Forg),
             "bcary_dark": Forg,
@@ -38,7 +38,9 @@ class AerosolBetaCaryophylleneDark(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(mode),
-                "kappa": self.kappa(mode),
+                "kappa": self.kappa(
+                    mass_fractions=mode, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(mode),
                 "spectrum": spectra.Lognormal(
                     norm_factor=N / si.cm**3, m_mode=50.0 * si.nm, s_geom=1.75
@@ -50,7 +52,7 @@ class AerosolBetaCaryophylleneDark(DryAerosolMixture):
 class AerosolBetaCaryophylleneLight(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.8, N: float = 400):
+    def __init__(self, water_molar_volume: float, Forg: float = 0.8, N: float = 400):
         mode = {
             "(NH4)2SO4": (1 - Forg),
             "bcary_light": Forg,
@@ -80,7 +82,9 @@ class AerosolBetaCaryophylleneLight(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(mode),
-                "kappa": self.kappa(mode),
+                "kappa": self.kappa(
+                    mass_fractions=mode, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(mode),
                 "spectrum": spectra.Lognormal(
                     norm_factor=N / si.cm**3, m_mode=50.0 * si.nm, s_geom=1.75
@@ -93,7 +97,7 @@ class AerosolBetaCaryophylleneLight(DryAerosolMixture):
 @strict
 class AerosolAlphaPineneDark(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.8, N: float = 400):
+    def __init__(self, water_molar_volume: float, Forg: float = 0.8, N: float = 400):
         mode = {
             "(NH4)2SO4": (1 - Forg),
             "apinene_dark": Forg,
@@ -123,7 +127,9 @@ class AerosolAlphaPineneDark(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(mode),
-                "kappa": self.kappa(mode),
+                "kappa": self.kappa(
+                    mass_fractions=mode, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(mode),
                 "spectrum": spectra.Lognormal(
                     norm_factor=N / si.cm**3, m_mode=50.0 * si.nm, s_geom=1.75
@@ -136,7 +142,7 @@ class AerosolAlphaPineneDark(DryAerosolMixture):
 @strict
 class AerosolAlphaPineneLight(DryAerosolMixture):
-    def __init__(self, Forg: float = 0.8, N: float = 400):
+    def __init__(self, water_molar_volume: float, Forg: float = 0.8, N: float = 400):
         mode = {
             "(NH4)2SO4": (1 - Forg),
             "apinene_light": Forg,
@@ -166,7 +172,9 @@ class AerosolAlphaPineneLight(DryAerosolMixture):
         self.modes = (
             {
                 "f_org": 1 - self.f_soluble_volume(mode),
-                "kappa": self.kappa(mode),
+                "kappa": self.kappa(
+                    mass_fractions=mode, water_molar_volume=water_molar_volume
+                ),
                 "nu_org": self.nu_org(mode),
                 "spectrum": spectra.Lognormal(
                     norm_factor=N / si.cm**3, m_mode=50.0 * si.nm, s_geom=1.75

PySDM_examples/Srivastava_1982/__init__.py CHANGED Viewed

@@ -2,6 +2,7 @@
 box-model examples feat. analytic solution for breakup from
 [Srivastava 1982 (JAS)](https://doi.org/10.1175/1520-0469(1982)039%3C1317:ASMOPC%3E2.0.CO;2)
 """
 from .equations import Equations, EquationsHelpers
 from .example import (
     add_to_plot_simulation_results,

PySDM_examples/Srivastava_1982/example.py CHANGED Viewed

@@ -137,10 +137,10 @@ def get_breakup_coalescence_analytic_results(equations, settings, m0, x, x_log):
 def get_analytic_results(equations, settings, mean_mass, mean_mass_ratio):
     res = {}
-    res[
-        SimProducts.Computed.mean_drop_volume_total_volume_ratio.name
-    ] = compute_drop_volume_total_volume_ratio(
-        mean_volume=mean_mass / settings.rho, total_volume=settings.total_volume
+    res[SimProducts.Computed.mean_drop_volume_total_volume_ratio.name] = (
+        compute_drop_volume_total_volume_ratio(
+            mean_volume=mean_mass / settings.rho, total_volume=settings.total_volume
+        )
     )
     res[SimProducts.PySDM.total_numer.name] = equations.M / mean_mass_ratio
     return res

PySDM_examples/Srivastava_1982/simulation.py CHANGED Viewed

@@ -70,9 +70,9 @@ class Simulation:
                     if step != 0:
                         particulator.run(steps=1)
                     for prod in self.settings.prods:
-                        self.simulation_res[n_sd][prod][seed][
-                            step
-                        ] = particulator.products[prod].get()
+                        self.simulation_res[n_sd][prod][seed][step] = (
+                            particulator.products[prod].get()
+                        )
                 np.testing.assert_allclose(
                     actual=self.simulation_res[n_sd][

PySDM_examples/Szumowski_et_al_1998/gui_viewer.py CHANGED Viewed

@@ -93,13 +93,15 @@ class GUIViewer:
                 self.spectrumPlots[key] = (
                     _SpectrumPlot(r_bins, self.settings.spectrum_per_mass_of_dry_air)
                     if key == "size"
-                    else _TemperaturePlot(
-                        self.settings.T_bins_edges, self.settings.formulae
-                    )
-                    if key == "temperature"
-                    else _TerminalVelocityPlot(
-                        self.settings.terminal_velocity_radius_bin_edges,
-                        self.settings.formulae,
+                    else (
+                        _TemperaturePlot(
+                            self.settings.T_bins_edges, self.settings.formulae
+                        )
+                        if key == "temperature"
+                        else _TerminalVelocityPlot(
+                            self.settings.terminal_velocity_radius_bin_edges,
+                            self.settings.formulae,
+                        )
                     )
                 )
                 clear_output()

PySDM_examples/Szumowski_et_al_1998/simulation.py CHANGED Viewed

@@ -162,11 +162,11 @@ class Simulation:
                 )
             if self.settings.freezing_singular:
-                attributes[
-                    "freezing temperature"
-                ] = formulae.freezing_temperature_spectrum.invcdf(
-                    np.random.random(immersed_surface_area.size),  # TODO #599: seed
-                    immersed_surface_area,
+                attributes["freezing temperature"] = (
+                    formulae.freezing_temperature_spectrum.invcdf(
+                        np.random.random(immersed_surface_area.size),  # TODO #599: seed
+                        immersed_surface_area,
+                    )
                 )
             else:
                 attributes["immersed surface area"] = immersed_surface_area

PySDM_examples/Van_Hook_1968/__init__.py CHANGED Viewed

@@ -2,4 +2,5 @@
 figure from Van Hook 1968 (J. Phys. Chem.)
 https://doi.org/10.1021/j100850a028
 """
 # pylint: disable=invalid-name

PySDM_examples/__init__.py CHANGED Viewed

@@ -2,6 +2,7 @@
 PySDM_examples package includes common Python modules used in PySDM smoke tests
 and in example notebooks (but the package wheels do not include the notebooks)
 """
 from importlib.metadata import PackageNotFoundError, version
 try:

PySDM_examples/deJong_Mackay_et_al_2023/__init__.py CHANGED Viewed

@@ -2,6 +2,7 @@
 box- and single-column breakup-focused examples from
 [de Jong et al. 2023](https://doi.org/10.5194/gmd-16-4193-2023)
 """
 # pylint: disable=invalid-name
 from .plot_rates import plot_ax, plot_zeros_ax
 from .settings1D import Settings1D

PySDM_examples/utils/__init__.py CHANGED Viewed

@@ -1,6 +1,7 @@
 """
 reusable commons for examples
 """
 from .basic_simulation import BasicSimulation
 from .dummy_controller import DummyController
 from .notebook_vars import notebook_vars

PySDM_examples/utils/notebook_vars.py CHANGED Viewed

@@ -1,4 +1,5 @@
 """ helper routines for use in smoke tests """
 from pathlib import Path
 import nbformat

{PySDM_examples-2.39.dist-info → PySDM_examples-2.41.dist-info}/METADATA RENAMED Viewed

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM-examples
-Version: 2.39
+Version: 2.41
 Summary: PySDM usage examples reproducing results from literature and depicting how to use PySDM from Python Jupyter notebooks
 Home-page: https://github.com/open-atmos/PySDM
 Author: https://github.com/open-atmos/PySDM/graphs/contributors
@@ -12,12 +12,14 @@ Project-URL: Source, https://github.com/open-atmos/PySDM/tree/main/examples
 Platform: UNKNOWN
 Description-Content-Type: text/markdown
 Requires-Dist: PySDM
+Requires-Dist: PyMPDATA
 Requires-Dist: open-atmos-jupyter-utils
 Requires-Dist: pystrict
+Requires-Dist: matplotlib
 Requires-Dist: joblib
 Requires-Dist: ipywidgets
 Requires-Dist: seaborn
-Requires-Dist: ghapi
+Requires-Dist: numdifftools
 Requires-Dist: pyvinecopulib
 Requires-Dist: vtk
 Provides-Extra: tests

{PySDM_examples-2.39.dist-info → PySDM_examples-2.41.dist-info}/RECORD RENAMED Viewed

@@ -1,75 +1,75 @@
-PySDM_examples/__init__.py,sha256=MQ7MlkaF_WPwhzoaaZrTB4OLY8iesN2iWD-1RV-T2j4,339
+PySDM_examples/__init__.py,sha256=Ht8EcJquTq4qeSaDJ4qFG8SNWA4qMfhfTyGcO415i3A,340
 PySDM_examples/Abdul_Razzak_Ghan_2000/__init__.py,sha256=t568yvoJqJgyyhR69eBjs2y8gj4SIgofowulmiWbz9U,165
-PySDM_examples/Abdul_Razzak_Ghan_2000/aerosol.py,sha256=n6mffeceJ0ls4aW16YMlTX_HCwy5W8GAROTuF_IxvwY,2841
+PySDM_examples/Abdul_Razzak_Ghan_2000/aerosol.py,sha256=BHZL7DILM2WGK5v251N0pHd_bdk8V1PnhSpItkPh1HI,3405
 PySDM_examples/Abdul_Razzak_Ghan_2000/data_from_ARG2000_paper.py,sha256=MRyurAnXq5ESE8gA5rzMsVx1ZhwcOz76o2xPRAcTPMQ,11585
 PySDM_examples/Abdul_Razzak_Ghan_2000/data_from_CloudMicrophysics_ARG.py,sha256=6vuCL55D2KnUutnqZvdv1Wnpgt0rz4BTFELKzOQP1uo,11547
-PySDM_examples/Abdul_Razzak_Ghan_2000/run_ARG_parcel.py,sha256=BKBMy_QxLOg0sEt05vz13BaYVKxXuV4jG-2NaIqy3Hc,5485
-PySDM_examples/Alpert_and_Knopf_2016/__init__.py,sha256=7LPh9TVrxIsUu_n_xAdJHqkXstHyNCeBK0Lqkf6wdoM,263
+PySDM_examples/Abdul_Razzak_Ghan_2000/run_ARG_parcel.py,sha256=kd1IH1XhhS6x-6wkdAHK4R6aIIgI_UXcYo1dDtivbOU,5402
+PySDM_examples/Alpert_and_Knopf_2016/__init__.py,sha256=NIqxh_uW0b83mrhYojbd1sMh6-lpz8wf8uuhtNtr0mU,264
 PySDM_examples/Alpert_and_Knopf_2016/simulation.py,sha256=KfXJ3Q2D1EG5EgvNfMc5HymzmpGNJ8f74hqkrR7KYAQ,10113
 PySDM_examples/Alpert_and_Knopf_2016/table.py,sha256=z3d4Pms0nwe4Jb10FilOtNcE4cpUx3zGtDSeStH1l6o,228
 PySDM_examples/Alpert_and_Knopf_2016/table_1.py,sha256=edJfFNqOYwZvOL4nKM5Ir5kFppDqbfz4hPUESqCnNaI,3792
 PySDM_examples/Alpert_and_Knopf_2016/table_2.py,sha256=oxeB_UNLWzrSHKoWOot-kAe5X_i9n18SKDOh7fF9KkQ,1857
-PySDM_examples/Arabas_and_Shima_2017/__init__.py,sha256=zap8wLAO5OwLNGGZF2E1FWd4X4viYCHDB29Inh7i2-Q,185
+PySDM_examples/Arabas_and_Shima_2017/__init__.py,sha256=hR7OQjmlyxwzrO61gQsfJaatHYI8BJBdwGSmgBZn8HI,186
 PySDM_examples/Arabas_and_Shima_2017/example.py,sha256=dWzGbRCkdwjWethxtdEjn6IjEssFrmW37sEK8y36Q3I,253
 PySDM_examples/Arabas_and_Shima_2017/settings.py,sha256=RxPoFAZ7K9EL9B8dVwj41wngTnCVPaFu5tOgFm_K8gU,2383
 PySDM_examples/Arabas_and_Shima_2017/simulation.py,sha256=p8zb3-dC0er8d5o7CKSlaM9lAgIzFwO7LqovvuQUrdk,3965
 PySDM_examples/Arabas_et_al_2015/__init__.py,sha256=uPSOIkvpr_UnLacbdQ5DYGW6k59JqvMVTKSuz6ff7bQ,244
 PySDM_examples/Arabas_et_al_2015/example.py,sha256=YB-ABtUByHQnv0jTGJvROhxNX_aQnM_PgoilG7mSEmQ,828
 PySDM_examples/Arabas_et_al_2015/example_benchmark.py,sha256=iUn5qGO0dB_s0DcYlOtllFURbRa6X9H9HW9AkaiHFoc,2343
-PySDM_examples/Arabas_et_al_2015/settings.py,sha256=OFEanJVLKj7EE3mZGCKtFRr7M84zkYt0jthMOwVXgfo,1623
+PySDM_examples/Arabas_et_al_2015/settings.py,sha256=HNMezPT8SK2bZlcLdoxoFAzczkz836fX0vRIb5XjntM,1599
 PySDM_examples/Arabas_et_al_2015/spin_up.py,sha256=OxISF9coJpQ3x-lUnyIkPbHsQ46U37ouQM_ijibZ1Zo,837
-PySDM_examples/Arabas_et_al_2023/__init__.py,sha256=v9P3lzVSv4KauX_uIGNoqaKNroj8joXCGzITXBakUaE,306
+PySDM_examples/Arabas_et_al_2023/__init__.py,sha256=zZzcM2x30oOKy8BXBtsDFHnrbemEj9lfBdproC_ccpM,307
 PySDM_examples/Arabas_et_al_2023/commons.py,sha256=BEg1kqBX0626TYTWAaJdMrOAspHdgkPkPoDm2cR8ZXk,574
 PySDM_examples/Arabas_et_al_2023/curved_text.py,sha256=DoUvIuHkZ8xrgyXwdwEjUZLMeFseyOG7HfPS78CRNnY,5338
 PySDM_examples/Arabas_et_al_2023/frozen_fraction.py,sha256=8ok3PbeXMwn87hCQt5tsNcHMCbYZHdW7dnCL-Fcx3LQ,783
 PySDM_examples/Arabas_et_al_2023/make_particulator.py,sha256=sPjM5wkljDajIh0XAWH8N22OBublwvnzknAv3rfIydI,2025
-PySDM_examples/Arabas_et_al_2023/plots.py,sha256=UEDAcWTEQ0U5yxBd4-sJNgkp46WzEUZNZ0SlThlCvS4,5007
+PySDM_examples/Arabas_et_al_2023/plots.py,sha256=fw7E3bkFiXLkB6QR8nK4qCMVDxW46B77LE-KpdepK_Q,5047
 PySDM_examples/Arabas_et_al_2023/run_simulation.py,sha256=LsZt0l4GMg9td2NwlyIo5agGiuiunkFRJ6rRVPuqfSs,1353
 PySDM_examples/Bartman_2020_MasterThesis/__init__.py,sha256=la2o79L6WLPueZ3EI5ThGtNRuSkXvKQwJJmvrNhhHa8,136
 PySDM_examples/Bartman_2020_MasterThesis/fig_4_adaptive_sdm.py,sha256=bbXDRKH7ekiwEd8dgNRcavlWzoQGfZvXKL_qOT6TtWU,2337
 PySDM_examples/Bartman_2020_MasterThesis/fig_5_SCIPY_VS_ADAPTIVE.py,sha256=CJ4oXcQe8LGxBoD352dyYQICTWkxPvfUVAhBVf5Uxls,4192
 PySDM_examples/Bartman_et_al_2021/__init__.py,sha256=G-GE8actgWv5-0x49DPJmAPR88e0ddlDGqIzzo5RQHU,88
-PySDM_examples/Bartman_et_al_2021/label.py,sha256=RGsSa-8AEzPbCU-YW_cZsOAkKnuqT67wuO7xYG3779c,507
-PySDM_examples/Berry_1967/__init__.py,sha256=t8YiMrYS6v_T0dGXmEhxauCQ8kvBABjgQ5U3wFTc2NY,184
+PySDM_examples/Bartman_et_al_2021/label.py,sha256=sHaZcODS8lrKxWha9Zd-VALHHqFTr_PsUrMPRunTxzA,523
+PySDM_examples/Berry_1967/__init__.py,sha256=WK67_TDhxgbKgiZWsdthnKvFcchvpJ8WuNcIP2b1cXY,185
 PySDM_examples/Berry_1967/example.py,sha256=r0WlZ0R4Y4RiDIsLlq3RG_R9YSvCSrizWGjg3mh0WgQ,2506
 PySDM_examples/Berry_1967/example_fig_6.py,sha256=wD11VB2Ocn-83jkLp55H5xquojdIw_-KdYtii5ZtkJY,3104
 PySDM_examples/Berry_1967/settings.py,sha256=9VmVS4zBTGQff-xBdxb9AjUB9iNNm-qMznURH0ESfwE,1617
 PySDM_examples/Berry_1967/spectrum_plotter.py,sha256=WEEaryYBDydStKyx6H915LVTxeFFKj0pmpdVotktjI4,1023
-PySDM_examples/Bieli_et_al_2022/__init__.py,sha256=74DQZDe3vAi5xmgsPBDS7SS8xRI_98_znPzwXitCWCs,139
+PySDM_examples/Bieli_et_al_2022/__init__.py,sha256=YdI2papu13HSfH973XxfrwGZkP_nbiXO4lgGHime_7M,141
 PySDM_examples/Bieli_et_al_2022/settings.py,sha256=qJxoGChzqta-c_O2yv4E3wEM4yoECa2dyhKmE8AGrtI,1679
 PySDM_examples/Bieli_et_al_2022/simulation.py,sha256=i-EEZvt37VuUXOEzbFBy7YrJez7_YHiI-gCtsPj_EbM,1304
-PySDM_examples/Bolot_et_al_2013/__init__.py,sha256=67DJlNlh37-EAj38-LyIBe9tNsgPGI0ks8ZSkt4c5iU,135
+PySDM_examples/Bolot_et_al_2013/__init__.py,sha256=6yZ-PdTBwwN9f4HqOe-E0m-CyM47mgg2E4ezVKVELF4,136
 PySDM_examples/Bulenok_2023_MasterThesis/__init__.py,sha256=GPHWvxMGulWIaSH7CQlESd2jNTRJeMcXRqjzzJ-AAf8,164
 PySDM_examples/Bulenok_2023_MasterThesis/setups.py,sha256=UVTp0oBMblcR46J29aylfZ6g6-8jlhXCb90ODdcfaSE,3424
 PySDM_examples/Bulenok_2023_MasterThesis/utils.py,sha256=Y0OXKd8D2vJ0vAstl_f68St39MRY89sMqaK1ASo-h1w,8019
-PySDM_examples/Ervens_and_Feingold_2012/__init__.py,sha256=cpE8Gj0BAeNWjKH6pGcVAgQ_3bFGHTG-pcm4gwmF2gc,158
-PySDM_examples/Ervens_and_Feingold_2012/settings.py,sha256=s5eHSjK4Y0QIB4hwm8ODkElHjLorK0F6GUy5KZOPmIc,464
+PySDM_examples/Ervens_and_Feingold_2012/__init__.py,sha256=F9Q-77npT6TyIK1vdNS4RPeBb8g4mFBUYz5qg4pj2wo,159
+PySDM_examples/Ervens_and_Feingold_2012/settings.py,sha256=DtkcDqkcO2WAkH63P1RQ19cejn0oAamdsZefN3y0Gy8,442
 PySDM_examples/Grabowski_and_Pawlowska_2023/__init__.py,sha256=_UWULcA_8_vA7zYlW_QuyS4SYe5827jyZBum-ZrGwGk,225
 PySDM_examples/Grabowski_and_Pawlowska_2023/settings.py,sha256=nErcGLTGJTai32jpRhrhpxoJ4Acwh1ccrwGBpQQjGcY,3273
 PySDM_examples/Grabowski_and_Pawlowska_2023/simulation.py,sha256=uXH0zfr4M4XZ81u5GNVFLlbQTSHoA15ScxdwR43X6w8,4515
 PySDM_examples/Graf_et_al_2019/__init__.py,sha256=nT6uOv3PODePRUzFoyOfH_eOCjJkg5z9qS3tXy3USPc,140
 PySDM_examples/Jaruga_and_Pawlowska_2018/__init__.py,sha256=VJz_xFOhVvnqYHKbAPRnZQ3Vky6JEOZxBx8kV0kHIcg,161
 PySDM_examples/Kreidenweis_et_al_2003/__init__.py,sha256=9QbrxhPe0QJxboGKeqHhq9rrvQ5VUcjTpqWgWsxnM1E,225
-PySDM_examples/Kreidenweis_et_al_2003/settings.py,sha256=70cH8qdhxGR2WrkptbSQJIzszJIxNJuAIAewJzp7g5o,2788
+PySDM_examples/Kreidenweis_et_al_2003/settings.py,sha256=iS7rZINn4EyxSsy_0CvQeGAz6ZyB-don1OCY7q_M21c,2836
 PySDM_examples/Kreidenweis_et_al_2003/simulation.py,sha256=UeU6SLhNAYNXVkw14PvPaM1efKmXYUzqkvpYI6L4ODs,4624
 PySDM_examples/Lowe_et_al_2019/__init__.py,sha256=JpFMpZWcBOXcQnYSXNoe4cc0giQJwJAm27dOewsR3wE,236
-PySDM_examples/Lowe_et_al_2019/aerosol.py,sha256=mtATQA5paEMf_Wo_RJf2D52Q7Mo8T2SjIqG9FGSF6pQ,6356
-PySDM_examples/Lowe_et_al_2019/aerosol_code.py,sha256=_HRzqJLst8bo9PBl_uNlb60LrSk0WB_0TYiI-JJmemw,8985
+PySDM_examples/Lowe_et_al_2019/aerosol.py,sha256=8h0e1QAo4CQyyKokYScfuiWbKVJuLHfJir_r3_KAnNk,6827
+PySDM_examples/Lowe_et_al_2019/aerosol_code.py,sha256=HqjEfWwsdKFnrCvRNpyUrfX4JGJLdS53AbHTJ0NgyzQ,9667
 PySDM_examples/Lowe_et_al_2019/plot_helper.py,sha256=Uc35rjIMvigpadiWhKvyNj6b9xauVKss17TkaO58fRQ,6450
 PySDM_examples/Lowe_et_al_2019/settings.py,sha256=J-tFTXjCMer0UJZT3o8xkyEcMWSYzZ6jLtahvaqiHYE,3096
 PySDM_examples/Lowe_et_al_2019/simulation.py,sha256=zy1RITIcYUhkgI8Cqcu13P6WTQ4NMjtNmTjKW_Xa7jQ,5143
-PySDM_examples/Merlivat_and_Nief_1967/__init__.py,sha256=Zjh0KimgWCpzQa6hrYhI_ByEJAFpCat9Dd4x9syHpsc,126
+PySDM_examples/Merlivat_and_Nief_1967/__init__.py,sha256=bbjcqXvQxmMwZIRTAOGWwFGkKGnQeGsi7x9kI-VdQpQ,127
 PySDM_examples/Morrison_and_Grabowski_2007/__init__.py,sha256=sae2OEASd-sh-f0U_8ObIGitWGiXAPXrr7aTn950xoM,181
 PySDM_examples/Morrison_and_Grabowski_2007/common.py,sha256=aOSCsXSCRPLEgV2LThpmAdF_gGANPAypUDxwGe5C2Do,5221
 PySDM_examples/Morrison_and_Grabowski_2007/cumulus.py,sha256=kyypwfFdYqY_QyoPr88ScvNb53rMTPvseRL7KVad9Fo,4490
 PySDM_examples/Morrison_and_Grabowski_2007/strato_cumulus.py,sha256=NsYP-AOHcY-HDvhg1dWAHWVFZzlaxpDoQJNbrSphdUg,1139
-PySDM_examples/Niedermeier_et_al_2014/__init__.py,sha256=Buya98WtrA4u9YWx4okx5018wSKmFbLTTyql63xQCVs,147
+PySDM_examples/Niedermeier_et_al_2014/__init__.py,sha256=Oy4G3rNeJtEvLy6F7vznxKL3vLYhPdhp2SoP4YyoPUQ,148
 PySDM_examples/Niedermeier_et_al_2014/settings.py,sha256=fNNORJlWJjKIJDou-UMj3D9EQvxhwa3at5wvemoQhi4,1614
-PySDM_examples/Niedermeier_et_al_2014/simulation.py,sha256=mG3OOF3b1M3OfmARL5i7vvAQk4ENXo099DBo-Iqgm8M,3344
-PySDM_examples/Pierchala_et_al_2022/__init__.py,sha256=KWD7hYVy3zIFrdxmPk7FhPurgbEbAxus-C7vLnwhFro,139
+PySDM_examples/Niedermeier_et_al_2014/simulation.py,sha256=F1DKol-ujtJ5PaNKyKATquuTgHwces2it4jXDu-4dho,3408
+PySDM_examples/Pierchala_et_al_2022/__init__.py,sha256=vPjlbuWAJ7karsAxnZi2QoJTvduugqCcvDgPe1aPxMA,140
 PySDM_examples/Pierchala_et_al_2022/commons.py,sha256=qqdQkEoURw7ewSCmsDz62lXDk20BQ_9xDmvxwOvznSg,1016
 PySDM_examples/Pyrcel/__init__.py,sha256=Hzj0Dq23fGz_QPHd9Fse76MG1LPJqk_Lpd2qzXCN8RY,208
-PySDM_examples/Pyrcel/profile_plotter.py,sha256=ijrzwnwpLXEYF_94xImyliQoV_mF56FrlKaw4XD9iXg,2180
+PySDM_examples/Pyrcel/profile_plotter.py,sha256=Yquz_ylj8W6KGecFY162PB4z7few2wldxJs5x7YyuW8,2196
 PySDM_examples/Pyrcel/settings.py,sha256=Kg-7_pN1Ah7P_Zc6Vg1mbk-oWZTDcDYM-qGec1BPsEI,2147
 PySDM_examples/Pyrcel/simulation.py,sha256=73Ao6u7F6e-pIrUhcYy7ZI3fVhDYg1TVSKbogv2uAVc,3774
 PySDM_examples/Pyrcel/tutorial_settings.py,sha256=C6Vxo-HAaxnRYTLv4yxjZC4qwTU7yuWvQR4ey6SDZ_s,2096
@@ -86,48 +86,48 @@ PySDM_examples/Shima_et_al_2009/tutorial_settings.py,sha256=nTe_rKAIPpXwnh52mDBr
 PySDM_examples/Shipway_and_Hill_2012/__init__.py,sha256=fVgTFoG-sHDYawCJhZdUrqqziLBiJKExw3_tAQeXQ5o,247
 PySDM_examples/Shipway_and_Hill_2012/mpdata_1d.py,sha256=GfbGXQFAVJ2rYnkhwKKESRLPDfJJeL0Qs1h1pZUjlWo,1907
 PySDM_examples/Shipway_and_Hill_2012/plot.py,sha256=Q6lnGWE3oSDKV_F1qKUIihhbYABQ3S1A1qg7Q2rXwn4,3606
-PySDM_examples/Shipway_and_Hill_2012/settings.py,sha256=j3YQ6Pba4lE6TQDJkWW0rId-zQgkn-fbXrTIgX9whRY,6097
+PySDM_examples/Shipway_and_Hill_2012/settings.py,sha256=350S9UiiRRDE240R6Q-XGKbF6n3-ImAzroeTtvrXvaI,6099
 PySDM_examples/Shipway_and_Hill_2012/simulation.py,sha256=jXEuytNmngdFaGIJSm0diauN11X1-EbW5Pbec7zptFY,9727
 PySDM_examples/Singer_Ward/__init__.py,sha256=wLus4IEzXkeOKCR3Dn2tM9PqNW9X7tjBiMymLvBow-c,31
-PySDM_examples/Singer_Ward/aerosol.py,sha256=Cu6o2bnitj_PDaLGoEXCWcml-2QrY8d3FLKN37WQbVk,5329
-PySDM_examples/Srivastava_1982/__init__.py,sha256=2m9a7s2gzNh4-dI3mwUfF9jEL4QK09m5vp3M45yFyug,508
+PySDM_examples/Singer_Ward/aerosol.py,sha256=H7TcvI94Wc-IvhTMFshFaJRJY3eUjihAlSdNxVKZzL0,5805
+PySDM_examples/Srivastava_1982/__init__.py,sha256=u7RxFGW08DIfwKU2yHpEQiiNlfz7F3GyYR4TV9LA7Ww,509
 PySDM_examples/Srivastava_1982/equations.py,sha256=RQHvXbznreupzhOdfSQpc859DUNhThjPEKt7sXTeEAw,3262
-PySDM_examples/Srivastava_1982/example.py,sha256=yIgZ60BeYhI7QRMmoe9l_Utc98wnWQ-4CCAL5hSIioU,8566
+PySDM_examples/Srivastava_1982/example.py,sha256=fGaLDvqz4BBLLG8NqHB2D2aAYFPOJhEg1r5DBPUntB8,8576
 PySDM_examples/Srivastava_1982/settings.py,sha256=f5R-MyEKEw4Lbfix_4eyjiP9o_SqabM4IOH1fvQ3jjA,2275
-PySDM_examples/Srivastava_1982/simulation.py,sha256=SbQSclTWZ_8eX8Ct2OoNx6GRWCQVloHv-9iKAj0jbjk,3034
+PySDM_examples/Srivastava_1982/simulation.py,sha256=V1ofPcEAI2wrxYptZxo8ibe32nekGeZvhL6jrcg2izc,3036
 PySDM_examples/Szumowski_et_al_1998/__init__.py,sha256=L9UucV9lyb08Hkzs1zYI9cxmN7mcogIFNsYgmnUpfY8,315
 PySDM_examples/Szumowski_et_al_1998/fields.py,sha256=tcR2MAB7azOLF955-ACHJU6v7aJxxpDXVdSIMQTbHEg,1584
 PySDM_examples/Szumowski_et_al_1998/gui.py,sha256=Vjt_ahdSl74NEjbZ-WxhK9eU3H2bqifjibPp6hy2DME,1309
 PySDM_examples/Szumowski_et_al_1998/gui_controller.py,sha256=AsP78e9bcGSP81kfDet_5fxpthlHsjDyJxHknvN2EgI,4283
 PySDM_examples/Szumowski_et_al_1998/gui_settings.py,sha256=3-Ox7GBHbEe30_znfry_UqAVU6ZTvJsCSRwMIw6_cJI,14400
-PySDM_examples/Szumowski_et_al_1998/gui_viewer.py,sha256=qnqE-Seh6rXieV_YGpBWXRY2O-KZMR_TTsZO3Xq5d_E,13088
+PySDM_examples/Szumowski_et_al_1998/gui_viewer.py,sha256=JZFAfD9FXHS5msVYAzccoF5rhHlb6hdwY28Q6lAdJD4,13164
 PySDM_examples/Szumowski_et_al_1998/make_default_product_collection.py,sha256=L7aUMbg0CVGGZXti8btqIh6UkKKkD-WblwN0EqURcLI,4666
 PySDM_examples/Szumowski_et_al_1998/mpdata_2d.py,sha256=41F3D4Kx_Um01wqOlkDSu8gzbwy47OewK535PK-Ii54,4137
 PySDM_examples/Szumowski_et_al_1998/plots.py,sha256=fzeluxlLF9VbSMmcu2wvd9WPl3ApIuA0jouwOh6rJHA,7001
-PySDM_examples/Szumowski_et_al_1998/simulation.py,sha256=HepfGl5HrBCk0nRQxdwL1EkAAoPcywuQ9UR5qCwfxcM,9927
+PySDM_examples/Szumowski_et_al_1998/simulation.py,sha256=bEkYxTiP8S7MYzz2SlPLqeOGLiGibs4uaE7lTOBD8AA,9941
 PySDM_examples/Szumowski_et_al_1998/sounding.py,sha256=vk3KHRPrBbZpz6QCqj2RtWyEpl_N7jSMR-zoDWM_ZVQ,1440
 PySDM_examples/Szumowski_et_al_1998/storage.py,sha256=NaGoujCRLubSEbhLA0QZDaSkATpJuF09HCnGdWQzcxQ,2579
-PySDM_examples/Van_Hook_1968/__init__.py,sha256=KUZAyoCgL4Mg2dKI9u-Aw58PXoUrWQsN7J4Sz8L-72M,118
+PySDM_examples/Van_Hook_1968/__init__.py,sha256=nTGyMcf6D4l6Ib0hMwyx4es1nX5dpuPbM9AYEleaIPE,119
 PySDM_examples/Yang_et_al_2018/__init__.py,sha256=zY7OnJ0d5eyWh3j0VvuHJMYOum9jY596-Ek2po0efmE,227
 PySDM_examples/Yang_et_al_2018/settings.py,sha256=J62kB5pEThmMcWSLsLKgcv0Pii6VLPDbfz6HPgs7O9o,2185
 PySDM_examples/Yang_et_al_2018/simulation.py,sha256=sG2kO_l66ec5XRYIDygB6_493jqKh6fF2hc_TyAh8z8,3928
-PySDM_examples/deJong_Mackay_et_al_2023/__init__.py,sha256=3LiuVu6tWH-xu3v59GqlxD7Map3Eq4ul1SlEW6zw3aI,500
+PySDM_examples/deJong_Mackay_et_al_2023/__init__.py,sha256=dxlC4ZYZDCfV5Zo2ichuMqE_ZrkrdyIsru8wijKJXyU,502
 PySDM_examples/deJong_Mackay_et_al_2023/plot_rates.py,sha256=6BJ9V4uaJu-6nPfkLqwevDEIXvo_stRok8unhBPbCHY,2421
 PySDM_examples/deJong_Mackay_et_al_2023/settings1D.py,sha256=2I0tK_GgWkQn2ZQjZdX4iUUPpOtL1dV8dOlEyaZ3rKQ,2951
 PySDM_examples/deJong_Mackay_et_al_2023/settings_0D.py,sha256=ksclyxZ0sZNeLS-QZC7Dop6feEmQA1KBIZNuhvxcMxs,1918
 PySDM_examples/deJong_Mackay_et_al_2023/simulation1D.py,sha256=xDA7I3EBoeAw_sHfp7MuBzoDrCbkF4wOKDR2BRePZxk,2430
 PySDM_examples/deJong_Mackay_et_al_2023/simulation_0D.py,sha256=EceMIsvm41kG_UocPCAW4lped1J3sBNqf4CLDf-f9CU,4522
 PySDM_examples/deJong_Mackay_et_al_2023/simulation_ss.py,sha256=fYE5LRPbt5qIYC3NaXhC-qsq9GDETBxmXf2Ih3ezFyw,5864
-PySDM_examples/utils/__init__.py,sha256=l5g2Dp3DYUFqm-UMtUaKQkxQnV-_vIv7qz2iqh14uyc,257
+PySDM_examples/utils/__init__.py,sha256=MayUnfhdnjWwc5Ept-yYfCS_tuSR8yqe-RdBl8FO3Xk,258
 PySDM_examples/utils/basic_simulation.py,sha256=4EdgD0p8EKNNcqYYVSlGm1MG2ac7_6heQdwS-kp-Tmw,698
 PySDM_examples/utils/dummy_controller.py,sha256=f1AjVNXPnDhhD15GkBFtJthDf7jZivDrppj79guYc_w,741
-PySDM_examples/utils/notebook_vars.py,sha256=iT8xwGpnkivmFm3AOEDqomcwIf9rz14UHn_A-cBQHyc,1276
+PySDM_examples/utils/notebook_vars.py,sha256=isvt3euX_v1pxOlWPQ720_1O3YrJoVPEmifIvucQjCo,1277
 PySDM_examples/utils/progbar_controller.py,sha256=IXk7rkXUOBKP7EQe6zLHUpnCPwtE2LuEeho5lxdsewA,470
 PySDM_examples/utils/read_vtk_1d.py,sha256=7_3sagEbX_FxqFuGkFWdwQapd1-9echFghU8TSNvGuw,966
 PySDM_examples/utils/widgets/__init__.py,sha256=U67WZ2YDcXSocMraI7pzcMxCfcxzXNqehDzoG6ogsh4,519
 PySDM_examples/utils/widgets/freezer.py,sha256=tFlZWShXUzA2X0cxgeE44N1g9kD33q3KQWZBqkgBP0A,407
 PySDM_examples/utils/widgets/progbar_updater.py,sha256=3uhG32QkAiSM3dTV7DoROJl7fAUNA6ZsJEIKSg3tdmk,267
-PySDM_examples-2.39.dist-info/METADATA,sha256=aiRIxT1XgspsGfO6yrahNW5We1hKyMzzmKktshX9R1g,30431
-PySDM_examples-2.39.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-PySDM_examples-2.39.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
-PySDM_examples-2.39.dist-info/RECORD,,
+PySDM_examples-2.41.dist-info/METADATA,sha256=JEge5iOnR-32BkTOp6EmumVHYLQdQtvv7LUO9y9dJyg,30488
+PySDM_examples-2.41.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
+PySDM_examples-2.41.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
+PySDM_examples-2.41.dist-info/RECORD,,

{PySDM_examples-2.39.dist-info → PySDM_examples-2.41.dist-info}/WHEEL RENAMED Viewed

File without changes

{PySDM_examples-2.39.dist-info → PySDM_examples-2.41.dist-info}/top_level.txt RENAMED Viewed

File without changes

PySDM-examples 2.39__py3-none-any.whl → 2.41__py3-none-any.whl

PySDM-examples 2.39py3-none-any.whl → 2.41py3-none-any.whl