PySDM-examples 2.101__tar.gz → 2.103__tar.gz

{pysdm_examples-2.101 → pysdm_examples-2.103}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: PySDM-examples
-Version: 2.101
+Version: 2.103
 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
@@ -160,3 +160,7 @@ The box environment is also featured in the <a href="https://open-atmos.github.i
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Bolin_1958.html">Bolin_1958</a>
 - below-cloud <mark>kinetic fractionation</mark>:
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Gedzelman_and_Arnold_1994.html">Gedzelman_and_Arnold_1994</a>
+
+## examples depicting extraterrestrial clouds (formulae-only, no simulations yet)
+- Titan (methane clouds):
+  - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Toon_et_al_1980.html">Toon_et_al_1980</a>

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Alpert_and_Knopf_2016/simulation.py

@@ -220,6 +220,7 @@ def simulation(
         environment=Box(dt=time_step, dv=volume),
     )
     builder.add_dynamic(Freezing(singular=False))
+    builder.request_attribute("volume")
 
     if hasattr(spectrum, "s_geom") and spectrum.s_geom == 1:
         _isa, _conc = np.full(n_sd, spectrum.m_mode), np.full(
@@ -230,7 +231,7 @@ def simulation(
     attributes = {
         "multiplicity": discretise_multiplicities(_conc * volume),
         "immersed surface area": _isa,
-        "volume": np.full(n_sd, droplet_volume),
+        "signed water mass": np.full(n_sd, droplet_volume * formulae.constants.rho_w),
     }
     np.testing.assert_almost_equal(attributes["multiplicity"], multiplicity)
     products = (

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Arabas_et_al_2023/make_particulator.py

@@ -27,8 +27,6 @@ def make_particulator(
     volume,
     thaw=False
 ):
-    attributes = {"volume": np.ones(n_sd) * droplet_volume}
-
     formulae_ctor_args = {
         "seed": seed,
         "constants": constants,
@@ -39,6 +37,10 @@ def make_particulator(
     formulae = Formulae(**formulae_ctor_args)
     backend = CPU(formulae, override_jit_flags={"parallel": False})
 
+    attributes = {
+        "signed water mass": np.ones(n_sd) * droplet_volume * formulae.constants.rho_w
+    }
+
     sampling = SpectroGlacialSampling(
         freezing_temperature_spectrum=formulae.freezing_temperature_spectrum,
         insoluble_surface_spectrum=ABIFM_spec,
@@ -65,6 +67,7 @@ def make_particulator(
     env["rhod"] = 1.0
 
     builder.add_dynamic(Freezing(singular=singular, thaw=thaw))
+    builder.request_attribute("volume")
 
     return builder.build(
         attributes=attributes,

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Morrison_and_Grabowski_2007/cumulus.py

@@ -50,7 +50,7 @@ class Cumulus(Common):
         def drhod_dz(z, _):
             lv = self.formulae.latent_heat.lv(T_of_z(z))
             return self.formulae.hydrostatics.drho_dz(
-                self.formulae.constants.g_std, p_of_z(z), T_of_z(z), q_of_z(z), lv
+                p_of_z(z), T_of_z(z), q_of_z(z), lv
             )
 
         theta_std0 = self.formulae.trivia.th_std(

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Morrison_and_Grabowski_2007/strato_cumulus.py

@@ -23,7 +23,6 @@ class StratoCumulus(Common):
             self.formulae.hydrostatics,
             "p_of_z_assuming_const_th_and_initial_water_vapour_mixing_ratio",
         )(
-            self.formulae.constants.g_std,
             self.p0,
             self.th_std0,
             self.initial_water_vapour_mixing_ratio,

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Niedermeier_et_al_2014/simulation.py

@@ -51,15 +51,15 @@ class Simulation(BasicSimulation):
             "multiplicity": ccn_conc_float * air_volume,
             "dry volume": dry_volume,
             "kappa times dry volume": settings.kappa * dry_volume,
-            "volume": None,
             "immersed surface area": immersed_surface_area,
         }
-        attributes["volume"] = settings.formulae.trivia.volume(
+        attributes["signed water mass"] = settings.formulae.trivia.volume(
             radius=equilibrate_wet_radii(
                 r_dry=ccn_diameter / 2,
                 environment=builder.particulator.environment,
                 kappa_times_dry_volume=attributes["kappa times dry volume"],
             )
+            * settings.formulae.constants.rho_w
         )
 
         for attribute, data in attributes.items():

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Shipway_and_Hill_2012/settings.py

@@ -99,7 +99,6 @@ class Settings:
             )
         )
 
-        g = self.formulae.constants.g_std
         self.rhod0 = self.formulae.state_variable_triplet.rho_d(
             p=p0,
             water_vapour_mixing_ratio=self.water_vapour_mixing_ratio(0 * si.m),
@@ -123,7 +122,7 @@ class Settings:
             )
             lv = self.formulae.latent_heat.lv(T)
             return self.formulae.hydrostatics.drho_dz(
-                g, p, T, water_vapour_mixing_ratio, lv
+                p, T, water_vapour_mixing_ratio, lv
             ) / (
                 1 + water_vapour_mixing_ratio
             ) - rhod * d_water_vapour_mixing_ratio__dz / (

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/Szumowski_et_al_1998/simulation.py

@@ -154,6 +154,8 @@ class Simulation:
         )
 
         if self.settings.processes["freezing"]:
+            attributes["signed water mass"] = attributes.pop("water mass")
+
             if self.settings.freezing_inp_spec is None:
                 immersed_surface_area = formulae.trivia.sphere_surface(
                     diameter=2 * formulae.trivia.radius(volume=attributes["dry volume"])

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples/seeding/settings.py

@@ -59,7 +59,7 @@ class Settings:
         v_dry = self.formulae.trivia.volume(radius=r_dry)
         self.seeded_particle_multiplicity = n_in_dv
         self.seeded_particle_extensive_attributes = {
-            "water mass": [0.0001 * si.ng] * n_sd_seeding,
+            "signed water mass": [0.0001 * si.ng] * n_sd_seeding,
             "dry volume": v_dry,
             "kappa times dry volume": 0.8 * v_dry,
         }

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: PySDM-examples
-Version: 2.101
+Version: 2.103
 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
@@ -160,3 +160,7 @@ The box environment is also featured in the <a href="https://open-atmos.github.i
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Bolin_1958.html">Bolin_1958</a>
 - below-cloud <mark>kinetic fractionation</mark>:
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Gedzelman_and_Arnold_1994.html">Gedzelman_and_Arnold_1994</a>
+
+## examples depicting extraterrestrial clouds (formulae-only, no simulations yet)
+- Titan (methane clouds):
+  - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Toon_et_al_1980.html">Toon_et_al_1980</a>

{pysdm_examples-2.101 → pysdm_examples-2.103}/PySDM_examples.egg-info/SOURCES.txt

@@ -131,6 +131,7 @@ PySDM_examples/Szumowski_et_al_1998/plots.py
 PySDM_examples/Szumowski_et_al_1998/simulation.py
 PySDM_examples/Szumowski_et_al_1998/sounding.py
 PySDM_examples/Szumowski_et_al_1998/storage.py
+PySDM_examples/Toon_et_al_1980/__init__.py
 PySDM_examples/Van_Hook_1968/__init__.py
 PySDM_examples/Yang_et_al_2018/__init__.py
 PySDM_examples/Yang_et_al_2018/settings.py

{pysdm_examples-2.101 → pysdm_examples-2.103}/docs/pysdm_examples_landing.md

@@ -124,3 +124,7 @@ The box environment is also featured in the <a href="https://open-atmos.github.i
   - `PySDM_examples.Bolin_1958`
 - below-cloud <mark>kinetic fractionation</mark>:
   - `PySDM_examples.Gedzelman_and_Arnold_1994`
+
+## examples depicting extraterrestrial clouds (formulae-only, no simulations yet)
+- Titan (methane clouds):
+  - `PySDM_examples.Toon_et_al_1980`

{pysdm_examples-2.101 → pysdm_examples-2.103}/pyproject.toml

@@ -4,4 +4,4 @@ local_scheme = "no-local-version"
 version_scheme = "post-release"
 
 [build-system]
-requires = ['setuptools==75.8.0', 'setuptools-scm==8.1.0']
+requires = ['setuptools==75.8.1', 'setuptools-scm==8.2.0']