PySDM-examples 3.0.0rc12__tar.gz → 3.0.0rc16__tar.gz

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pysdm-examples
-Version: 3.0.0rc12
+Version: 3.0.0rc16
 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
@@ -15,7 +15,6 @@ Requires-Dist: PyMPDATA
 Requires-Dist: open-atmos-jupyter-utils
 Requires-Dist: pystrict
 Requires-Dist: matplotlib
-Requires-Dist: Pillow
 Requires-Dist: joblib
 Requires-Dist: ipywidgets
 Requires-Dist: seaborn
@@ -26,11 +25,10 @@ Requires-Dist: pyvinecopulib
 Requires-Dist: networkx
 Provides-Extra: ci-version-pins
 Requires-Dist: PySDM[CI_version_pins]; extra == "ci-version-pins"
-Requires-Dist: PyMPDATA==1.6.1; extra == "ci-version-pins"
+Requires-Dist: PyMPDATA==1.6.3; extra == "ci-version-pins"
 Requires-Dist: open-atmos-jupyter-utils==1.3.0; extra == "ci-version-pins"
 Requires-Dist: pystrict==1.3; extra == "ci-version-pins"
 Requires-Dist: matplotlib!=3.9.1; extra == "ci-version-pins"
-Requires-Dist: Pillow<11.3.0; extra == "ci-version-pins"
 Requires-Dist: joblib==1.5.1; extra == "ci-version-pins"
 Requires-Dist: ipywidgets==8.1.7; extra == "ci-version-pins"
 Requires-Dist: seaborn==0.13.2; extra == "ci-version-pins"
@@ -137,6 +135,7 @@ Example notebooks include:
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Yang_et_al_2018.html">Yang_et_al_2018</a>: polydisperse particle spectrum, activation/deactivation cycles
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Abdul_Razzak_Ghan_2000.html">Abdul_Razzak_Ghan_2000</a>: polydisperse activation, comparison against <mark>GCM parameterisation</mark>
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Pyrcel.html">Pyrcel</a>: polydisperse activation, mimicking example test case from <mark>Pyrcel</mark> documentation
+  - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Strzabala_2025_BEng.html">Strzabala_2025_BEng</a>: ParaView visualisation example
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Lowe_et_al_2019.html">Lowe_et_al_2019</a>: externally mixed polydisperse size spectrum with <mark>surface-active organics</mark> case
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Grabowski_and_Pawlowska_2023.html">Grabowski_and_Pawlowska_2023</a>: polydisperse activation, focus on <mark>ripening</mark>
   - <a href="https://open-atmos.github.io/PySDM/PySDM_examples/Jensen_and_Nugent_2017.html">Jensen_and_Nugent_2017</a>: polydisperse activation featuring <mark>giant CCN</mark>
@@ -148,6 +147,10 @@ Example notebooks include:
 
 The parcel environment is also featured in the <a href="https://open-atmos.github.io/PySDM/PySDM.html#tutorials">PySDM tutorials</a>.
 
+<center>
+  <img src="https://github.com/open-atmos/PySDM/releases/download/tip/parcel_animation_ubuntu-24.04.gif" style="width:80%" />
+</center>
+
 ## 0D box environment
 
 The box environment is void of any spatial or thermodynamic context, it constitutes the most basic framework.

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Abade_and_Albuquerque_2024/simulation.py

@@ -49,9 +49,9 @@ class Simulation(BasicSimulation):
         if settings.enable_vapour_deposition_on_ice:
             builder.add_dynamic(VapourDepositionOnIce(adaptive=True))
 
-        r_dry, n_in_dv = ConstantMultiplicity(settings.soluble_aerosol).sample(
-            n_sd=settings.n_sd
-        )
+        r_dry, n_in_dv = ConstantMultiplicity(
+            settings.soluble_aerosol
+        ).sample_deterministic(n_sd=settings.n_sd)
         attributes = builder.particulator.environment.init_attributes(
             n_in_dv=n_in_dv, kappa=settings.kappa, r_dry=r_dry
         )

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Abdul_Razzak_Ghan_2000/run_ARG_parcel.py

@@ -61,7 +61,9 @@ def run_parcel(
     }
     for i, mode in enumerate(aerosol.modes):
         kappa, spectrum = mode["kappa"]["CompressedFilmOvadnevaite"], mode["spectrum"]
-        r_dry, concentration = ConstantMultiplicity(spectrum).sample(n_sd_per_mode)
+        r_dry, concentration = ConstantMultiplicity(spectrum).sample_deterministic(
+            n_sd_per_mode
+        )
         v_dry = builder.formulae.trivia.volume(radius=r_dry)
         specific_concentration = concentration / builder.formulae.constants.rho_STP
         attributes["multiplicity"] = np.append(

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Alpert_and_Knopf_2016/simulation.py

@@ -227,7 +227,9 @@ def simulation(
             n_sd, multiplicity / volume
         )
     else:
-        _isa, _conc = spectral_sampling.ConstantMultiplicity(spectrum).sample(n_sd)
+        _isa, _conc = spectral_sampling.ConstantMultiplicity(
+            spectrum
+        ).sample_deterministic(n_sd)
     attributes = {
         "multiplicity": discretise_multiplicities(_conc * volume),
         "immersed surface area": _isa,

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Bieli_et_al_2022/simulation.py

@@ -18,7 +18,7 @@ def make_core(settings, coal_eff):
     attributes = {}
     attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
         settings.spectrum
-    ).sample(settings.n_sd)
+    ).sample_deterministic(settings.n_sd)
     collision = Collision(
         collision_kernel=settings.kernel,
         coalescence_efficiency=coal_eff,

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Ervens_and_Feingold_2012/settings.py

@@ -9,4 +9,4 @@ def sampled_ccn_diameter_number_concentration_spectrum(
     return Logarithmic(
         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)
+    ).sample_deterministic(n_sd)

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Grabowski_and_Pawlowska_2023/simulation.py

@@ -17,7 +17,6 @@ class Simulation(BasicSimulation):
         settings,
         products=None,
         scipy_solver=False,
-        sampling_class=ConstantMultiplicity,
     ):
         builder = Builder(
             n_sd=settings.n_sd,
@@ -55,7 +54,9 @@ class Simulation(BasicSimulation):
         kappa = tuple(settings.aerosol_modes_by_kappa.keys())[0]
         spectrum = settings.aerosol_modes_by_kappa[kappa]
 
-        r_dry, n_per_volume = sampling_class(spectrum).sample(settings.n_sd)
+        r_dry, n_per_volume = ConstantMultiplicity(spectrum).sample_deterministic(
+            settings.n_sd
+        )
         v_dry = settings.formulae.trivia.volume(radius=r_dry)
         attributes["multiplicity"] = np.append(
             attributes["multiplicity"], n_per_volume * volume

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Jensen_and_Nugent_2017/simulation.py

@@ -65,7 +65,7 @@ class Simulation(BasicSimulation):
 
         self.r_dry, n_in_unit_volume = Logarithmic(
             spectrum=settings.dry_radii_spectrum,
-        ).sample(builder.particulator.n_sd - n_gccn)
+        ).sample_deterministic(builder.particulator.n_sd - n_gccn)
 
         if gccn:
             nonzero_concentration_mask = np.nonzero(table_3.NA)

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Kreidenweis_et_al_2003/settings.py

@@ -17,7 +17,8 @@ class Settings:
         dt: float,
         n_sd: int,
         n_substep: int,
-        spectral_sampling: spec_sampling.SpectralSampling = spec_sampling.Logarithmic,
+        spectral_sampling_class: spec_sampling.SpectralSampling = spec_sampling.Logarithmic,
+        spectral_sampling_method: str = "sample_deterministic",
     ):
         self.formulae = Formulae(
             saturation_vapour_pressure="AugustRocheMagnus",
@@ -51,13 +52,16 @@ class Settings:
         # note: rho is not specified in the paper
         rho0 = 1
 
-        self.r_dry, self.n_in_dv = spectral_sampling(
-            spectrum=spectra.Lognormal(
-                norm_factor=566 / si.cm**3 / rho0 * self.mass_of_dry_air,
-                m_mode=0.08 * si.um / 2,
-                s_geom=2,
-            )
-        ).sample(n_sd)
+        self.r_dry, self.n_in_dv = getattr(
+            spectral_sampling_class(
+                spectrum=spectra.Lognormal(
+                    norm_factor=566 / si.cm**3 / rho0 * self.mass_of_dry_air,
+                    m_mode=0.08 * si.um / 2,
+                    s_geom=2,
+                )
+            ),
+            spectral_sampling_method,
+        )(n_sd)
 
         self.ENVIRONMENT_MOLE_FRACTIONS = {
             "SO2": 0.2 * PPB,

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Lowe_et_al_2019/settings.py

@@ -7,7 +7,6 @@ from pystrict import strict
 from PySDM import Formulae
 from PySDM.backends import CPU
 from PySDM.initialisation.aerosol_composition import DryAerosolMixture
-from PySDM.initialisation.sampling import spectral_sampling as spec_sampling
 from PySDM.physics import si
 
 
@@ -36,7 +35,6 @@ class Settings:
         n_sd_per_mode: int,
         aerosol: DryAerosolMixture,
         model: str,
-        spectral_sampling: type(spec_sampling.SpectralSampling),
         w: float = 0.32 * si.m / si.s,
     ):
         assert model in ("Constant", "CompressedFilmOvadnevaite")
@@ -46,7 +44,6 @@ class Settings:
         self.formulae = self.backend.formulae
         const = self.formulae.constants
         self.aerosol = aerosol
-        self.spectral_sampling = spectral_sampling
 
         max_altitude = 200 * si.m
         self.w = w

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Lowe_et_al_2019/simulation.py

@@ -7,6 +7,7 @@ from PySDM.dynamics import AmbientThermodynamics, Condensation
 from PySDM.environments import Parcel
 from PySDM.initialisation.hygroscopic_equilibrium import equilibrate_wet_radii
 from PySDM.initialisation.spectra import Sum
+from PySDM.initialisation.sampling.spectral_sampling import ConstantMultiplicity
 
 
 class Simulation(BasicSimulation):
@@ -33,9 +34,9 @@ class Simulation(BasicSimulation):
         }
         initial_volume = settings.mass_of_dry_air / settings.rho0
         for mode in settings.aerosol.modes:
-            r_dry, n_in_dv = settings.spectral_sampling(
+            r_dry, n_in_dv = ConstantMultiplicity(
                 spectrum=mode["spectrum"]
-            ).sample(settings.n_sd_per_mode)
+            ).sample_deterministic(settings.n_sd_per_mode)
             v_dry = settings.formulae.trivia.volume(radius=r_dry)
             attributes["multiplicity"] = np.append(
                 attributes["multiplicity"], n_in_dv * initial_volume

pysdm_examples-3.0.0rc16/PySDM_examples/Matsushima_et_al_2023/__init__.py

@@ -0,0 +1,7 @@
+# pylint: disable=invalid-name
+"""
+Fig. 1 from [Matsushima et al. 2023 (GMD)](https://doi.org/10.5194/gmd-16-6211-2023)
+
+figure_1.ipynb:
+.. include:: ./figure_1.ipynb.badges.md
+"""

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Pyrcel/simulation.py

@@ -13,7 +13,13 @@ from PySDM.physics import si
 
 class Simulation(BasicSimulation):
     def __init__(
-        self, settings, products=None, scipy_solver=False, rtol_thd=1e-10, rtol_x=1e-10
+        self,
+        settings,
+        products=None,
+        scipy_solver=False,
+        rtol_thd=1e-10,
+        rtol_x=1e-10,
+        mass_of_dry_air=44 * si.kg,
     ):
         n_sd = sum(settings.n_sd_per_mode)
         builder = Builder(
@@ -27,7 +33,7 @@ class Simulation(BasicSimulation):
                 initial_water_vapour_mixing_ratio=settings.initial_vapour_mixing_ratio,
                 T0=settings.initial_temperature,
                 w=settings.vertical_velocity,
-                mass_of_dry_air=44 * si.kg,
+                mass_of_dry_air=mass_of_dry_air,
             ),
         )
         builder.add_dynamic(AmbientThermodynamics())
@@ -43,7 +49,9 @@ class Simulation(BasicSimulation):
         }
         for i, (kappa, spectrum) in enumerate(settings.aerosol_modes_by_kappa.items()):
             sampling = ConstantMultiplicity(spectrum)
-            r_dry, n_per_volume = sampling.sample(settings.n_sd_per_mode[i])
+            r_dry, n_per_volume = sampling.sample_deterministic(
+                settings.n_sd_per_mode[i]
+            )
             v_dry = settings.formulae.trivia.volume(radius=r_dry)
             attributes["multiplicity"] = np.append(
                 attributes["multiplicity"], n_per_volume * volume

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Shima_et_al_2009/example.py

@@ -22,7 +22,9 @@ def run(settings, backend=CPU, observers=()):
     )
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
-    attributes["volume"], attributes["multiplicity"] = sampling.sample(settings.n_sd)
+    attributes["volume"], attributes["multiplicity"] = sampling.sample_deterministic(
+        settings.n_sd
+    )
     coalescence = Coalescence(
         collision_kernel=settings.kernel, adaptive=settings.adaptive
     )

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Shima_et_al_2009/example_timing.py

@@ -16,7 +16,9 @@ def run(settings, backend):
     builder = Builder(n_sd=settings.n_sd, backend=backend, environment=env)
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
-    attributes["volume"], attributes["multiplicity"] = sampling.sample(settings.n_sd)
+    attributes["volume"], attributes["multiplicity"] = sampling.sample_deterministic(
+        settings.n_sd
+    )
     builder.add_dynamic(Coalescence(collision_kernel=settings.kernel))
     particles = builder.build(attributes, products=(WallTime(),))
 

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Shima_et_al_2009/tutorial_example.py

@@ -11,7 +11,7 @@ def run(settings, observers=()):
     builder.set_environment(Box(dv=settings.dv, dt=settings.dt))
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
-    attributes["volume"], attributes["n"] = sampling.sample(settings.n_sd)
+    attributes["volume"], attributes["n"] = sampling.sample_deterministic(settings.n_sd)
     coalescence = Coalescence(
         collision_kernel=settings.kernel, adaptive=settings.adaptive
     )

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Spichtinger_et_al_2023/settings.py

@@ -59,7 +59,7 @@ class Settings:
 
         self.r_dry, self.specific_concentration = spectral_sampling.Linear(
             spectrum
-        ).sample(n_sd)
+        ).sample_deterministic(n_sd)
 
         self.t_duration = 7200
         self.dt = dt

pysdm_examples-3.0.0rc16/PySDM_examples/Strzabala_2025_BEng/__init__.py

@@ -0,0 +1,8 @@
+# pylint: disable=invalid-name
+"""
+ParaView visualisation of a parcel-model example based on the test case from
+[Pyrcel package docs](https://pyrcel.readthedocs.io/)
+
+paraview_parcel_model.ipynb:
+.. include:: ./paraview.ipynb.badges.md
+"""

pysdm_examples-3.0.0rc16/PySDM_examples/Strzabala_2025_BEng/paraview_parcel_model.py

@@ -0,0 +1,244 @@
+#!/usr/bin/env pvpython
+"""
+ParaView pvpython script for visualizing sd_products and sd_attributes.
+"""
+
+import argparse
+from collections import namedtuple
+import pathlib
+from paraview import simple as pvs  # pylint: disable=import-error
+
+pvs._DisableFirstRenderCameraReset()  # pylint: disable=protected-access
+
+
+def cli_using_argparse(argparse_parser):
+    """
+    Command line interface using argparse.
+    """
+    argparse_parser.add_argument(
+        "--sd-products-pvd",
+        dest="sd_products_pvd",
+        help="Path to sd_products.pvd",
+    )
+    argparse_parser.add_argument(
+        "--sd-attributes-pvd",
+        dest="sd_attributes_pvd",
+        help="Path to sd_attributes.pvd",
+    )
+    argparse_parser.add_argument(
+        "--output-animation-path",
+        help="Output path for the animation file.",
+    )
+    argparse_parser.add_argument(
+        "--output-screenshot-path",
+        help="Output path for the screenshot file.",
+    )
+
+    argparse_parser.add_argument(
+        "--scalarbar-title-size",
+        type=int,
+        default=30,
+        help="Font size for scalar bar titles.",
+    )
+    argparse_parser.add_argument(
+        "--scalarbar-label-size",
+        type=int,
+        default=30,
+        help="Font size for scalar bar labels.",
+    )
+    argparse_parser.add_argument(
+        "--axes-title-size",
+        type=int,
+        default=30,
+        help="Font size for axes titles.",
+    )
+    argparse_parser.add_argument(
+        "--axes-label-size",
+        type=int,
+        default=30,
+        help="Font size for axes labels.",
+    )
+
+
+parser = argparse.ArgumentParser(
+    description="ParaView pvpython script for visualizing sd_products and sd_attributes PVD files."
+)
+cli_using_argparse(parser)
+args = parser.parse_args()
+
+sd_productspvd = pvs.PVDReader(
+    registrationName="sd_products.pvd", FileName=args.sd_products_pvd
+)
+sd_attributespvd = pvs.PVDReader(
+    registrationName="sd_attributes.pvd", FileName=args.sd_attributes_pvd
+)
+
+setup = {
+    "renderView1": pvs.GetActiveViewOrCreate("RenderView"),
+    "sd_productspvdDisplay": pvs.Show(
+        sd_productspvd,
+        pvs.GetActiveViewOrCreate("RenderView"),
+        "UnstructuredGridRepresentation",
+    ),
+    "sd_attributespvdDisplay": pvs.Show(
+        sd_attributespvd,
+        pvs.GetActiveViewOrCreate("RenderView"),
+        "UnstructuredGridRepresentation",
+    ),
+    "rHlookup_table": pvs.GetColorTransferFunction("RH"),
+    "volumelookup_table": pvs.GetColorTransferFunction("volume"),
+}
+
+animation_setup = namedtuple("setup", setup.keys())(**setup)
+
+
+def configure_color_bar_and_display(
+    display, lookup_table, kind, *, anim_setup=animation_setup
+):
+    """
+    Configure color bar and display settings.
+    """
+    display.Representation = "Surface"
+    display.SetScalarBarVisibility(anim_setup.renderView1, True)
+    color_bar = pvs.GetScalarBar(lookup_table, anim_setup.renderView1)
+    color_bar.LabelColor = [0.0, 0.0, 0.0]
+    color_bar.DrawScalarBarOutline = 1
+    color_bar.ScalarBarOutlineColor = [0.0, 0.0, 0.0]
+    color_bar.TitleColor = [0.0, 0.0, 0.0]
+
+    color_bar.TitleFontSize = args.scalarbar_title_size
+    color_bar.LabelFontSize = args.scalarbar_label_size
+
+    if kind == "prod":
+        display.Opacity = 0.4
+        display.DisableLighting = 1
+        display.Diffuse = 0.76
+        lookup_table.RescaleTransferFunction(90.0, 101.0)
+        lookup_table.ApplyPreset("Black, Blue and White", True)
+        lookup_table.NanColor = [0.67, 1.0, 1.0]
+    else:
+        display.PointSize = 13.0
+        display.RenderPointsAsSpheres = 1
+        display.Interpolation = "PBR"
+        lookup_table.RescaleTransferFunction(1e-18, 1e-13)
+        lookup_table.ApplyPreset("Cold and Hot", True)
+        lookup_table.MapControlPointsToLogSpace()
+        lookup_table.UseLogScale = 1
+        lookup_table.NumberOfTableValues = 16
+        lookup_table.InvertTransferFunction()
+
+
+def configure_data_axes_grid(display, kind):
+    """
+    Configure data axes grid settings.
+    """
+    display.DataAxesGrid.GridAxesVisibility = 1
+    display.DataAxesGrid.XTitle = ""
+    display.DataAxesGrid.YTitle = ""
+    display.DataAxesGrid.XAxisUseCustomLabels = 1
+    display.DataAxesGrid.YAxisUseCustomLabels = 1
+
+    display.DataAxesGrid.XTitleFontSize = args.axes_title_size
+    display.DataAxesGrid.XLabelFontSize = args.axes_label_size
+    display.DataAxesGrid.YTitleFontSize = args.axes_title_size
+    display.DataAxesGrid.YLabelFontSize = args.axes_label_size
+
+    if kind == "prod":
+        display.DataAxesGrid.ZTitle = ""
+        display.DataAxesGrid.GridColor = [0.0, 0.0, 0.0]
+        display.DataAxesGrid.ShowGrid = 1
+        display.DataAxesGrid.LabelUniqueEdgesOnly = 0
+        display.DataAxesGrid.ZAxisUseCustomLabels = 1
+        display.DataAxesGrid.ZTitleFontSize = args.axes_title_size
+        display.DataAxesGrid.ZLabelFontSize = args.axes_label_size
+    else:
+        display.DataAxesGrid.ZTitle = "     Z [m]     "
+        display.DataAxesGrid.ZLabelColor = [0.0, 0.0, 0.0]
+        display.DataAxesGrid.ZTitleColor = [0.0, 0.0, 0.0]
+        display.DataAxesGrid.ZTitleFontSize = args.axes_title_size
+        display.DataAxesGrid.ZLabelFontSize = args.axes_label_size
+
+
+def configure_view_appearance(render_view):
+    """
+    Configure view appearance settings."""
+    render_view.OrientationAxesLabelColor = [0.0, 0.0, 0.0]
+    render_view.OrientationAxesOutlineColor = [0.0, 0.0, 0.0]
+    render_view.OrientationAxesXVisibility = 0
+    render_view.OrientationAxesYVisibility = 0
+    render_view.OrientationAxesZColor = [0.0, 0.0, 0.0]
+    render_view.UseColorPaletteForBackground = 0
+    render_view.Background = [1.0, 1.0, 1.0]
+
+
+def set_camera_view(render_view):
+    """
+    Set camera view settings.
+    """
+    layout1 = pvs.GetLayout()
+    layout1.SetSize(1592, 1128)
+    render_view.CameraPosition = [
+        1548.95,
+        -1349.49,
+        699.27,
+    ]
+    render_view.CameraFocalPoint = [
+        -1.37e-13,
+        3.18e-13,
+        505.00,
+    ]
+    render_view.CameraViewUp = [
+        -0.07,
+        0.06,
+        0.99,
+    ]
+    render_view.CameraParallelScale = 534.08
+
+
+def save_animation_and_screenshot(render_view, animation_path, screenshot_path):
+    """
+    Save animation and screenshot.
+    """
+    animation_scene = pvs.GetAnimationScene()
+    animation_scene.UpdateAnimationUsingDataTimeSteps()
+
+    pvs.SaveAnimation(animation_path, render_view, FrameRate=15)
+
+    if not screenshot_path:
+        return
+
+    time_steps = sd_productspvd.TimestepValues
+    if not time_steps:
+        return
+
+    last_time = time_steps[-1]
+    render_view.ViewTime = last_time
+
+    for reader in (sd_productspvd, sd_attributespvd):
+        reader.UpdatePipeline(last_time)
+
+    pvs.ExportView(
+        filename=str(pathlib.Path(screenshot_path)),
+        view=render_view,
+        Rasterize3Dgeometry=False,
+        GL2PSdepthsortmethod="BSP sorting (slow, best)",
+    )
+
+    pvs.RenderAllViews()
+
+
+configure_color_bar_and_display(
+    animation_setup.sd_productspvdDisplay, animation_setup.rHlookup_table, kind="prod"
+)
+configure_data_axes_grid(animation_setup.sd_productspvdDisplay, kind="prod")
+configure_view_appearance(animation_setup.renderView1)
+configure_color_bar_and_display(
+    animation_setup.sd_attributespvdDisplay,
+    animation_setup.volumelookup_table,
+    kind="attr",
+)
+configure_data_axes_grid(animation_setup.sd_attributespvdDisplay, kind="attr")
+set_camera_view(animation_setup.renderView1)
+save_animation_and_screenshot(
+    animation_setup.renderView1, args.output_animation_path, args.output_screenshot_path
+)

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Yang_et_al_2018/settings.py

@@ -19,9 +19,7 @@ class Settings:
         dt_max: float = 1 * si.second,
     ):
         self.total_time = 3 * si.hours
-        self.mass_of_dry_air = (
-            1000 * si.kilogram
-        )  # TODO #335 doubled with jupyter si unit
+        self.mass_of_dry_air = 1000 * si.kilogram
 
         self.n_steps = int(
             self.total_time / (5 * si.second)
@@ -34,7 +32,7 @@ class Settings:
                 s_geom=1.4,
             ),
             size_range=(10.633 * si.nanometre, 513.06 * si.nanometre),
-        ).sample(n_sd)
+        ).sample_deterministic(n_sd)
         self.dt_max = dt_max
 
         self.dt_output = dt_output

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/Yang_et_al_2018/simulation.py

@@ -84,10 +84,10 @@ class Simulation:
         _sp = self.particulator
         cell_id = 0
         output["r_bins_values"].append(
-            _sp.products["Particles Wet Size Spectrum"].get()
+            (_sp.products["Particles Wet Size Spectrum"].get()).T
         )
         volume = _sp.attributes["volume"].to_ndarray()
-        output["r"].append(self.formulae.trivia.radius(volume=volume))
+        output["r"].append((self.formulae.trivia.radius(volume=volume)).T)
         output["S"].append(_sp.environment["RH"][cell_id] - 1)
         output["t"].append(_sp.products["t"].get())
         for key in ("water_vapour_mixing_ratio", "T", "z"):
@@ -124,4 +124,6 @@ class Simulation:
         for _ in range(self.n_steps):
             self.particulator.run(self.n_substeps)
             self.save(output)
+        for k, v in output.items():
+            output[k] = np.asarray(v)
         return output

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/__init__.py

@@ -11,4 +11,4 @@ except PackageNotFoundError:
     # package is not installed
     pass
 
-assert PySDM.__version__ == __version__
+# assert PySDM.__version__ == __version__

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/deJong_Azimi/simulation_0D.py

@@ -25,7 +25,7 @@ def run_box(settings, backend_class=CPU):
     attributes = {}
     attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
         settings.spectrum
-    ).sample(settings.n_sd)
+    ).sample_deterministic(settings.n_sd)
     builder.add_dynamic(
         Coalescence(
             collision_kernel=settings.kernel,

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/deJong_Mackay_et_al_2023/simulation_0D.py

@@ -35,15 +35,15 @@ def run_box_breakup(
     builder.particulator.environment["rhod"] = 1.0
     attributes = {}
     if sample_in_radius:
-        diams, attributes["multiplicity"] = Logarithmic(settings.spectrum).sample(
-            settings.n_sd
-        )
+        diams, attributes["multiplicity"] = Logarithmic(
+            settings.spectrum
+        ).sample_deterministic(settings.n_sd)
         radii = diams / 2
         attributes["volume"] = Formulae().trivia.volume(radius=radii)
     else:
         attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
             settings.spectrum
-        ).sample(settings.n_sd)
+        ).sample_deterministic(settings.n_sd)
     breakup = Collision(
         collision_kernel=settings.kernel,
         coalescence_efficiency=settings.coal_eff,
@@ -100,7 +100,7 @@ def run_box_NObreakup(settings, steps=None, backend_class=CPU):
     attributes = {}
     attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
         settings.spectrum
-    ).sample(settings.n_sd)
+    ).sample_deterministic(settings.n_sd)
     coal = Coalescence(
         collision_kernel=settings.kernel,
         coalescence_efficiency=settings.coal_eff,

{pysdm_examples-3.0.0rc12 → pysdm_examples-3.0.0rc16}/PySDM_examples/seeding/settings.py

@@ -53,7 +53,7 @@ class Settings:
                 m_mode=1 * si.um,
                 s_geom=1.1,
             )
-        ).sample(
+        ).sample_deterministic(
             n_sd=n_sd_seeding
         )  # TODO #1387: does not to be the same?
         v_dry = self.formulae.trivia.volume(radius=r_dry)