PySDM-examples 2.65__py3-none-any.whl → 2.67__py3-none-any.whl

PySDM_examples/Abade_and_Albuquerque_2024/__init__.py

@@ -0,0 +1,8 @@
+# pylint: disable=invalid-name
+"""
+mixed-phase example using parcel environment based on
+[Abade & Albuquerque 2024 (QJRMS)](https://doi.org/10.1002/qj.4775)
+"""
+
+from .simulation import Simulation
+from .settings import Settings

PySDM_examples/Abade_and_Albuquerque_2024/settings.py

@@ -0,0 +1,48 @@
+from pystrict import strict
+
+from PySDM import Formulae
+from PySDM.physics import si
+from PySDM.initialisation.spectra import Lognormal
+
+
+@strict
+class Settings:
+    def __init__(self, *, n_sd: int, timestep: float):
+        self.n_sd = n_sd
+        self.timestep = timestep
+
+        self.initial_total_pressure = 1000 * si.hPa  # note: not given in the paper
+
+        # parameters from the paper
+        self.formulae = Formulae(
+            constants={"bulk_phase_partitioning_exponent": 0.1},
+            bulk_phase_partitioning="KaulEtAl2015",
+        )
+        self.initial_water_vapour_mixing_ratio = 1.5 * si.g / si.kg
+        self.parcel_linear_extent = 100 * si.m
+        self.updraft = 0.5 * si.m / si.s
+
+        thd_0 = self.formulae.state_variable_triplet.th_dry(
+            th_std=269 * si.K,
+            water_vapour_mixing_ratio=self.initial_water_vapour_mixing_ratio,
+        )
+        rhod_0 = self.formulae.state_variable_triplet.rho_d(
+            p=self.initial_total_pressure,
+            water_vapour_mixing_ratio=self.initial_water_vapour_mixing_ratio,
+            theta_std=thd_0,
+        )
+
+        self.mass_of_dry_air = rhod_0 * self.formulae.trivia.volume(
+            radius=self.parcel_linear_extent
+        )
+        self.soluble_aerosol = Lognormal(
+            norm_factor=200
+            / si.mg
+            * self.mass_of_dry_air,  # note: assuming per mg of dry air
+            m_mode=75 * si.nm,
+            s_geom=1.6,
+        )
+        self.kappa = 0.6
+        self.initial_temperature = self.formulae.state_variable_triplet.T(
+            rhod=rhod_0, thd=thd_0
+        )

PySDM_examples/Abade_and_Albuquerque_2024/simulation.py

@@ -0,0 +1,49 @@
+import numpy as np
+
+from PySDM_examples.utils import BasicSimulation
+
+from PySDM import Builder
+from PySDM.backends import CPU
+from PySDM.dynamics import Condensation, AmbientThermodynamics
+from PySDM.initialisation.sampling.spectral_sampling import ConstantMultiplicity
+from PySDM.products import AmbientTemperature, ParcelDisplacement, WaterMixingRatio
+from PySDM.environments import Parcel
+
+
+class Simulation(BasicSimulation):
+    def __init__(self, settings):
+        parcel = Parcel(
+            dt=settings.timestep,
+            mass_of_dry_air=settings.mass_of_dry_air,
+            p0=settings.initial_total_pressure,
+            initial_water_vapour_mixing_ratio=settings.initial_water_vapour_mixing_ratio,
+            T0=settings.initial_temperature,
+            w=settings.updraft,
+        )
+        builder = Builder(
+            backend=CPU(settings.formulae, override_jit_flags={"parallel": False}),
+            n_sd=settings.n_sd,
+            environment=parcel,
+        )
+        builder.add_dynamic(AmbientThermodynamics())
+        builder.add_dynamic(Condensation())
+
+        r_dry, n_in_dv = ConstantMultiplicity(settings.soluble_aerosol).sample(
+            n_sd=settings.n_sd
+        )
+        attributes = parcel.init_attributes(
+            n_in_dv=n_in_dv, kappa=settings.kappa, r_dry=r_dry
+        )
+        self.products = (
+            WaterMixingRatio(name="water", radius_range=(0, np.inf)),
+            WaterMixingRatio(name="ice", radius_range=(-np.inf, 0)),
+            WaterMixingRatio(name="total", radius_range=(-np.inf, np.inf)),
+            ParcelDisplacement(name="height"),
+            AmbientTemperature(name="T"),
+        )
+        super().__init__(
+            particulator=builder.build(attributes=attributes, products=self.products)
+        )
+
+    def run(self, *, nt, steps_per_output_interval):
+        return self._run(nt=nt, steps_per_output_interval=steps_per_output_interval)

{PySDM_examples-2.65.dist-info → PySDM_examples-2.67.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM-examples
-Version: 2.65
+Version: 2.67
 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
@@ -238,6 +238,13 @@ Requires-Dist: nbconvert ; extra == 'tests'
     [![launch on mybinder.org](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/open-atmos/PySDM.git/main?urlpath=lab/tree/examples/PySDM_examples/Jaruga_and_Pawlowska_2018/fig_3.ipynb)
     [![launch on Colab](https://colab.research.google.com/assets/colab-badge.svg)](https://colab.research.google.com/github/open-atmos/PySDM/blob/main/examples/PySDM_examples/Jaruga_and_Pawlowska_2018/fig_3.ipynb)  
 
+### 0D parcel-environment condensation/freezing example:
+- [Abade & Albuquerque 2024](https://doi.org/10.1002/qj.4775) (Adiabatic parcel, liquid/solid phase partitioning):
+  - Fig 2: 
+    [![View notebook](https://img.shields.io/static/v1?label=render%20on&logo=github&color=87ce3e&message=GitHub)](https://github.com/open-atmos/PySDM/blob/main/examples/PySDM_examples/Abade_and_Albuquerque_2024/fig_2.ipynb)
+    [![launch on mybinder.org](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/open-atmos/PySDM.git/main?urlpath=lab/tree/examples/PySDM_examples/Abade_and_Albuquerque_2024/fig_2.ipynb)
+    [![launch on Colab](https://colab.research.google.com/assets/colab-badge.svg)](https://colab.research.google.com/github/open-atmos/PySDM/blob/main/examples/PySDM_examples/Abade_and_Albuquerque_2024/fig_2.ipynb)   
+
 ### OD parcel-environment iterative framework mimicking removal of precipitation
 - [Rozanski & Sonntag 1982](https://doi.org/10.1111/j.2153-3490.1982.tb01800.x):
   - Figs 4-6:

{PySDM_examples-2.65.dist-info → PySDM_examples-2.67.dist-info}/RECORD

@@ -1,4 +1,7 @@
 PySDM_examples/__init__.py,sha256=Ht8EcJquTq4qeSaDJ4qFG8SNWA4qMfhfTyGcO415i3A,340
+PySDM_examples/Abade_and_Albuquerque_2024/__init__.py,sha256=agkn4Z6rWIh64M4UweHkacNu-dNOE2dkSfPqvy2zD2k,228
+PySDM_examples/Abade_and_Albuquerque_2024/settings.py,sha256=bKByjymGguWI47OvgVZrTdFoClYb7FmAmOAa6vIhwqg,1613
+PySDM_examples/Abade_and_Albuquerque_2024/simulation.py,sha256=rGIbkOa9XcwLmdjUbyKgFM6n6TLZyZVBD-4l-hWEGUs,1925
 PySDM_examples/Abdul_Razzak_Ghan_2000/__init__.py,sha256=t568yvoJqJgyyhR69eBjs2y8gj4SIgofowulmiWbz9U,165
 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
@@ -142,7 +145,7 @@ PySDM_examples/utils/read_vtk_1d.py,sha256=7_3sagEbX_FxqFuGkFWdwQapd1-9echFghU8T
 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.65.dist-info/METADATA,sha256=6Gjm58LhtEsHGAFAAVXh1OD6LtDKqpHU04LTDPTo5hI,39921
-PySDM_examples-2.65.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-PySDM_examples-2.65.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
-PySDM_examples-2.65.dist-info/RECORD,,
+PySDM_examples-2.67.dist-info/METADATA,sha256=h5s_z6_pGZ14Y3TSgihv2Yjg-YG6T8hi6C0nVI71Eqo,40752
+PySDM_examples-2.67.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+PySDM_examples-2.67.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
+PySDM_examples-2.67.dist-info/RECORD,,