PySDM-examples 3.0.0rc31__tar.gz → 3.0.0rc33__tar.gz

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pysdm-examples
-Version: 3.0.0rc31
+Version: 3.0.0rc33
 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

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Abade_and_Albuquerque_2024/simulation.py

@@ -22,6 +22,24 @@ from PySDM.environments import Parcel
 
 class Simulation(BasicSimulation):
     def __init__(self, settings):
+
+        dynamics = [
+            AmbientThermodynamics(),
+            Condensation(),
+        ]
+
+        if settings.enable_immersion_freezing:
+            dynamics.append(
+                Freezing(
+                    immersion_freezing=(
+                        "singular" if settings.singular else "time-dependent"
+                    )
+                )
+            )
+
+        if settings.enable_vapour_deposition_on_ice:
+            dynamics.append(VapourDepositionOnIce(adaptive=True))
+
         builder = Builder(
             backend=settings.backend,
             n_sd=settings.n_sd,
@@ -34,20 +52,8 @@ class Simulation(BasicSimulation):
                 w=settings.updraft,
                 mixed_phase=True,
             ),
+            dynamics=dynamics,
         )
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation())
-
-        if settings.enable_immersion_freezing:
-            builder.add_dynamic(
-                Freezing(
-                    immersion_freezing=(
-                        "singular" if settings.singular else "time-dependent"
-                    )
-                )
-            )
-        if settings.enable_vapour_deposition_on_ice:
-            builder.add_dynamic(VapourDepositionOnIce(adaptive=True))
 
         r_dry, n_in_dv = ConstantMultiplicity(
             settings.soluble_aerosol

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Abdul_Razzak_Ghan_2000/run_ARG_parcel.py

@@ -52,9 +52,12 @@ def run_parcel(
     )
     n_sd = n_sd_per_mode * len(aerosol.modes)
 
-    builder = Builder(backend=CPU(formulae), n_sd=n_sd, environment=env)
-    builder.add_dynamic(AmbientThermodynamics())
-    builder.add_dynamic(Condensation())
+    builder = Builder(
+        backend=CPU(formulae),
+        n_sd=n_sd,
+        environment=env,
+        dynamics=[AmbientThermodynamics(), Condensation()],
+    )
     builder.request_attribute("critical saturation")
 
     attributes = {

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Alpert_and_Knopf_2016/simulation.py

@@ -218,8 +218,8 @@ def simulation(
         n_sd=n_sd,
         backend=CPU(formulae=formulae),
         environment=Box(dt=time_step, dv=volume),
+        dynamics=[Freezing(immersion_freezing="time-dependent")],
     )
-    builder.add_dynamic(Freezing(immersion_freezing="time-dependent"))
     builder.request_attribute("volume")
 
     if hasattr(spectrum, "s_geom") and spectrum.s_geom == 1:

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Arabas_and_Pawlowska_2011/simulation.py

@@ -32,11 +32,12 @@ class Simulation(BasicSimulation):
             backend=CPU(settings.formulae),
             n_sd=settings.n_sd,
             environment=environment,
+            dynamics=(
+                AmbientThermodynamics(),
+                Condensation(),
+            ),
         )
 
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation())
-
         builder.request_attribute("radius")
 
         self.builder = builder

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Arabas_and_Shima_2017/simulation.py

@@ -36,16 +36,16 @@ class Simulation:
                 T0=settings.T0,
                 w=settings.w,
             ),
+            dynamics=[
+                AmbientThermodynamics(),
+                Condensation(
+                    rtol_x=settings.rtol_x,
+                    rtol_thd=settings.rtol_thd,
+                    dt_cond_range=settings.dt_cond_range,
+                ),
+            ],
         )
 
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(
-            Condensation(
-                rtol_x=settings.rtol_x,
-                rtol_thd=settings.rtol_thd,
-                dt_cond_range=settings.dt_cond_range,
-            )
-        )
         attributes = {}
         r_dry = np.array([settings.r_dry])
         attributes["dry volume"] = settings.formulae.trivia.volume(radius=r_dry)

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Arabas_et_al_2025/make_particulator.py

@@ -65,16 +65,18 @@ def make_particulator(
         thaw_flag = None
     attributes["multiplicity"] *= total_particle_number
 
-    builder = Builder(n_sd=n_sd, backend=backend, environment=Box(dt, volume))
+    builder = Builder(
+        n_sd=n_sd,
+        backend=backend,
+        environment=Box(dt, volume),
+        dynamics=[Freezing(immersion_freezing=immersion_freezing_flag, thaw=thaw_flag)],
+    )
 
     env = builder.particulator.environment
     env["T"] = initial_temperature
     env["RH"] = A_VALUE_LARGER_THAN_ONE
     env["rhod"] = 1.0
 
-    builder.add_dynamic(
-        Freezing(immersion_freezing=immersion_freezing_flag, thaw=thaw_flag)
-    )
     builder.request_attribute("volume")
 
     return builder.build(

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Bieli_et_al_2022/simulation.py

@@ -8,25 +8,24 @@ from PySDM.initialisation.sampling.spectral_sampling import ConstantMultiplicity
 
 def make_core(settings, coal_eff):
     backend = CPU
-
+    collision = Collision(
+        collision_kernel=settings.kernel,
+        coalescence_efficiency=coal_eff,
+        breakup_efficiency=settings.break_eff,
+        fragmentation_function=settings.fragmentation,
+        adaptive=settings.adaptive,
+    )
     builder = Builder(
         n_sd=settings.n_sd,
         backend=backend(settings.formulae),
         environment=Box(dv=settings.dv, dt=settings.dt),
+        dynamics=[collision],
     )
     builder.particulator.environment["rhod"] = 1.0
     attributes = {}
     attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
         settings.spectrum
     ).sample_deterministic(settings.n_sd)
-    collision = Collision(
-        collision_kernel=settings.kernel,
-        coalescence_efficiency=coal_eff,
-        breakup_efficiency=settings.break_eff,
-        fragmentation_function=settings.fragmentation,
-        adaptive=settings.adaptive,
-    )
-    builder.add_dynamic(collision)
     common_args = {
         "attr": "volume",
         "attr_unit": "m^3",

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Grabowski_and_Pawlowska_2023/simulation.py

@@ -31,10 +31,10 @@ class Simulation(BasicSimulation):
                 w=settings.vertical_velocity,
                 mass_of_dry_air=44 * si.kg,
             ),
-        )
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(
-            Condensation(rtol_thd=settings.rtol_thd, rtol_x=settings.rtol_x)
+            dynamics=(
+                AmbientThermodynamics(),
+                Condensation(rtol_thd=settings.rtol_thd, rtol_x=settings.rtol_x),
+            ),
         )
         for attribute in (
             "critical saturation",

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Jensen_and_Nugent_2017/simulation.py

@@ -46,19 +46,22 @@ class Simulation(BasicSimulation):
                 w=settings.vertical_velocity,
                 z0=settings.z0,
             ),
+            dynamics=[
+                # TODO #1266: order matters here, but error message is not saying it!
+                AmbientThermodynamics(),
+                Condensation(),
+            ]
+            + (
+                []
+                if not gravitational_coalsecence
+                else [Coalescence(collision_kernel=Geometric())]
+            ),
         )
 
         additional_derived_attributes = ("radius", "equilibrium saturation")
         for additional_derived_attribute in additional_derived_attributes:
             builder.request_attribute(additional_derived_attribute)
 
-        builder.add_dynamic(
-            AmbientThermodynamics()
-        )  # TODO #1266: order matters here, but error message is not saying it!
-        builder.add_dynamic(Condensation())
-        if gravitational_coalsecence:
-            builder.add_dynamic(Coalescence(collision_kernel=Geometric()))
-
         self.r_dry, n_in_unit_volume = Logarithmic(
             spectrum=settings.dry_radii_spectrum,
         ).sample_deterministic(builder.particulator.n_sd - n_gccn)

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Kreidenweis_et_al_2003/simulation.py

@@ -25,6 +25,17 @@ class Simulation(BasicSimulation):
                 T0=settings.T0,
                 w=settings.w,
             ),
+            dynamics=(
+                AmbientThermodynamics(),
+                Condensation(),
+                AqueousChemistry(
+                    environment_mole_fractions=settings.ENVIRONMENT_MOLE_FRACTIONS,
+                    system_type=settings.system_type,
+                    n_substep=settings.n_substep,
+                    dry_rho=settings.DRY_RHO,
+                    dry_molar_mass=settings.dry_molar_mass,
+                ),
+            ),
         )
 
         attributes = builder.particulator.environment.init_attributes(
@@ -38,18 +49,6 @@ class Simulation(BasicSimulation):
             **settings.starting_amounts,
         }
 
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation())
-        builder.add_dynamic(
-            AqueousChemistry(
-                environment_mole_fractions=settings.ENVIRONMENT_MOLE_FRACTIONS,
-                system_type=settings.system_type,
-                n_substep=settings.n_substep,
-                dry_rho=settings.DRY_RHO,
-                dry_molar_mass=settings.dry_molar_mass,
-            )
-        )
-
         products = products or (
             PySDM_products.AmbientRelativeHumidity(name="RH", unit="%"),
             PySDM_products.WaterMixingRatio(

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Long_1974/simulation.py

@@ -15,19 +15,20 @@ class Simulation(BasicSimulation):
             n_sd=settings.n_sd,
             backend=CPU(),
             environment=Box(dv=settings.dv, dt=settings.dt),
+            dynamics=(
+                Coalescence(
+                    collision_kernel=settings.kernel,
+                    coalescence_efficiency=settings.coal_eff,
+                    adaptive=settings.adaptive,
+                ),
+            ),
         )
         builder.particulator.environment["rhod"] = settings.rhod
         attributes = {}
         attributes["volume"], attributes["multiplicity"] = ConstantMultiplicity(
             settings.spectrum
         ).sample_deterministic(settings.n_sd)
-        builder.add_dynamic(
-            Coalescence(
-                collision_kernel=settings.kernel,
-                coalescence_efficiency=settings.coal_eff,
-                adaptive=settings.adaptive,
-            )
-        )
+
         products = (
             ParticleVolumeVersusRadiusLogarithmSpectrum(
                 radius_bins_edges=settings.radius_bins_edges, name="dv/dlnr"

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Lowe_et_al_2019/simulation.py

@@ -24,6 +24,7 @@ class Simulation(BasicSimulation):
                 T0=settings.T0,
                 w=settings.w,
             ),
+            dynamics=(AmbientThermodynamics(), Condensation()),
         )
 
         attributes = {
@@ -73,9 +74,6 @@ class Simulation(BasicSimulation):
         if settings.model == "Constant":
             del attributes["dry volume organic"]
 
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation())
-
         products = products or (
             PySDM_products.ParcelDisplacement(name="z"),
             PySDM_products.Time(name="t"),

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Luettmer_homogeneous_freezing/plot.py

@@ -16,10 +16,8 @@ T_frz_bins_kelvin = np.linspace(230, 240, num=100, endpoint=True)
 
 
 def cumulative_histogram(data, bins, reverse=False, density=True):
-    # Compute regular histogram using given bins
     hist, bin_edges = np.histogram(data, bins=bins, density=False)
 
-    # Cumulative sum
     if reverse:
         cum_hist = np.cumsum(hist[::-1])[::-1]
         cum_hist_0 = cum_hist[0]
@@ -27,7 +25,6 @@ def cumulative_histogram(data, bins, reverse=False, density=True):
         cum_hist = np.cumsum(hist)
         cum_hist_0 = cum_hist[-1]
 
-    # Normalize
     if density:
         cum_hist = cum_hist / cum_hist_0
 
@@ -56,8 +53,7 @@ def plot_thermodynamics_and_bulk(
     qv = np.asarray(output["qv"])
     T_frz = np.asarray(output["T_frz"])
     if show_conc:
-        nc = np.asarray(output["ns"])
-        ni = np.asarray(output["ni"])
+        nc, ni = np.asarray(output["ns"]), np.asarray(output["ni"])
 
     if t_lim is None:
         t_lim = np.amax(time)
@@ -67,8 +63,7 @@ def plot_thermodynamics_and_bulk(
     first_T_frz = T[first_ice_idx]
 
     if not show_tf:
-        rc = np.asarray(output["rs"])
-        ri = np.asarray(output["ri"])
+        rc, ri = np.asarray(output["rs"]), np.asarray(output["ri"])
 
     if show_jhom:
         svp = Formulae(
@@ -90,9 +85,7 @@ def plot_thermodynamics_and_bulk(
         radius = np.asarray(output["radius"])
         multiplicity = np.asarray(output["multiplicity"])
 
-    _, axs = pyplot.subplots(
-        1, 4, figsize=(20, 5), sharex=False, constrained_layout=True
-    )
+    _, axs = pyplot.subplots(1, 4, figsize=(20, 5), constrained_layout=True)
 
     # Temperture profile
     iax = 0
@@ -247,30 +240,40 @@ def plot_thermodynamics_and_bulk(
     ax.grid(visible=True)
 
 
-def plot_freezing_temperatures_histogram(ax, simulation):
+def plot_freezing_temperatures_histogram(ax, simulation, plot_rhi=False):
 
     number_of_ensemble_runs = simulation["settings"]["number_of_ensemble_runs"]
 
     for i in range(number_of_ensemble_runs):
         output = simulation["ensemble_member_outputs"][i]
-        T_frz = np.asarray(output["T_frz"])
+        if plot_rhi:
+            var = np.asarray(output["RHi_frz"])
+            bins = np.linspace(1, 1.6, num=60, endpoint=True)
+            cumulative = 1
+        else:
+            var = np.asarray(output["T_frz"])
+            bins = T_frz_bins_kelvin
+            cumulative = -1
         title = "Nucleation rate=" + simulation["settings"]["hom_freezing"]
 
-        # Freezing temperatures
         ax.hist(
-            T_frz,
-            bins=T_frz_bins_kelvin,
+            var,
+            bins=bins,
             density=True,
-            cumulative=-1,
+            cumulative=cumulative,
             alpha=1.0,
             histtype="step",
             linewidth=1.5,
         )
 
-        ax.set_xlim(left=234, right=239)
-        ax.axvline(x=235, color="k", linestyle="--")
+        if plot_rhi:
+            ax.set_xlim(left=1.0, right=1.6)
+            ax.set_xlabel("freezing supersaturation wrt ice", fontsize=ax_lab_fsize)
+        else:
+            ax.set_xlim(left=234, right=239)
+            ax.axvline(x=235, color="k", linestyle="--")
+            ax.set_xlabel("freezing temperature [K]", fontsize=ax_lab_fsize)
         ax.set_title(title, fontsize=ax_lab_fsize)
-        ax.set_xlabel("freezing temperature [K]", fontsize=ax_lab_fsize)
         ax.set_ylabel("frozen fraction", fontsize=ax_lab_fsize)
         ax.tick_params(labelsize=tick_fsize)
     return ax

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Luettmer_homogeneous_freezing/simulation.py

@@ -35,16 +35,23 @@ class Simulation:
             backend=settings.backend,
             n_sd=settings.n_sd,
             environment=env,
-        )
-
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation(adaptive=True))
-        if settings.deposition_enable:
-            builder.add_dynamic(VapourDepositionOnIce(adaptive=True))
-        builder.add_dynamic(
-            Freezing(
-                homogeneous_freezing=settings.hom_freezing_type, immersion_freezing=None
-            )
+            dynamics=(
+                [
+                    AmbientThermodynamics(),
+                    Condensation(adaptive=True),
+                ]
+                + (
+                    [VapourDepositionOnIce(adaptive=True)]
+                    if settings.deposition_enable
+                    else []
+                )
+                + [
+                    Freezing(
+                        homogeneous_freezing=settings.hom_freezing_type,
+                        immersion_freezing=None,
+                    )
+                ]
+            ),
         )
 
         self.n_sd = settings.n_sd
@@ -72,6 +79,7 @@ class Simulation:
             "signed water mass": self.initial_mass,
         }
         builder.request_attribute("temperature of last freezing")
+        builder.request_attribute("supersaturation of last freezing")
         builder.request_attribute("radius")
         builder.request_attribute("wet to critical volume ratio")
 
@@ -133,6 +141,9 @@ class Simulation:
         output["T_frz"] = self.particulator.attributes[
             "temperature of last freezing"
         ].data.tolist()
+        output["RHi_frz"] = self.particulator.attributes[
+            "supersaturation of last freezing"
+        ].data.tolist()
         if not output["radius"]:
             output["radius"] = self.particulator.attributes["radius"].data.tolist()
             output["multiplicity"] = self.particulator.attributes[
@@ -146,6 +157,7 @@ class Simulation:
 
         output = {
             "T_frz": [],
+            "RHi_frz": [],
             "radius": [],
             "multiplicity": [],
         }

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Murray_et_al_2010/simulation.py

@@ -136,10 +136,11 @@ def simulation(
         n_sd=n_sd,
         backend=CPU(formulae=formulae),
         environment=Box(dt=time_step, dv=volume),
+        dynamics=(
+            Freezing(homogeneous_freezing="time-dependent", immersion_freezing=None),
+        ),
     )
-    builder.add_dynamic(
-        Freezing(homogeneous_freezing="time-dependent", immersion_freezing=None)
-    )
+
     builder.request_attribute("temperature of last freezing")
     builder.request_attribute("volume")
     droplet_volume = formulae.trivia.volume(radius=droplet_radius)

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Niedermeier_et_al_2014/simulation.py

@@ -32,12 +32,13 @@ class Simulation(BasicSimulation):
                 override_jit_flags={"parallel": False},
             ),
             environment=env,
+            dynamics=(
+                AmbientThermodynamics(),
+                Condensation(),
+                Freezing(immersion_freezing="time-dependent"),
+            ),
         )
 
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation())
-        builder.add_dynamic(Freezing(immersion_freezing="time-dependent"))
-
         air_volume = settings.mass_of_dry_air / settings.rhod0
         (
             ccn_diameter,

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Pyrcel/simulation.py

@@ -37,12 +37,14 @@ class Simulation(BasicSimulation):
                 w=settings.vertical_velocity,
                 mass_of_dry_air=mass_of_dry_air,
             ),
+            dynamics=(
+                AmbientThermodynamics(),
+                Condensation(rtol_thd=rtol_thd, rtol_x=rtol_x),
+            ),
         )
         if additional_attributes is not None:
             for attribute in additional_attributes:
                 builder.request_attribute(attribute)
-        builder.add_dynamic(AmbientThermodynamics())
-        builder.add_dynamic(Condensation(rtol_thd=rtol_thd, rtol_x=rtol_x))
 
         volume = (
             builder.particulator.environment.mass_of_dry_air

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Shima_et_al_2009/example.py

@@ -19,16 +19,15 @@ def run(settings, backend=CPU, observers=()):
         n_sd=settings.n_sd,
         backend=backend(formulae=settings.formulae),
         environment=env,
+        dynamics=(
+            Coalescence(collision_kernel=settings.kernel, adaptive=settings.adaptive),
+        ),
     )
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
     attributes["volume"], attributes["multiplicity"] = sampling.sample_deterministic(
         settings.n_sd
     )
-    coalescence = Coalescence(
-        collision_kernel=settings.kernel, adaptive=settings.adaptive
-    )
-    builder.add_dynamic(coalescence)
     products = (
         ParticleVolumeVersusRadiusLogarithmSpectrum(
             settings.radius_bins_edges, name="dv/dlnr"

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Shima_et_al_2009/example_timing.py

@@ -13,13 +13,17 @@ from PySDM.products import WallTime
 
 def run(settings, backend):
     env = Box(dv=settings.dv, dt=settings.dt)
-    builder = Builder(n_sd=settings.n_sd, backend=backend, environment=env)
+    builder = Builder(
+        n_sd=settings.n_sd,
+        backend=backend,
+        environment=env,
+        dynamics=(Coalescence(collision_kernel=settings.kernel),),
+    )
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
     attributes["volume"], attributes["multiplicity"] = sampling.sample_deterministic(
         settings.n_sd
     )
-    builder.add_dynamic(Coalescence(collision_kernel=settings.kernel))
     particles = builder.build(attributes, products=(WallTime(),))
 
     states = {}

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Shima_et_al_2009/tutorial_example.py

@@ -11,16 +11,15 @@ def run(settings, observers=()):
         n_sd=settings.n_sd,
         backend=CPU(formulae=settings.formulae),
         environment=Box(dv=settings.dv, dt=settings.dt),
+        dynamics=(
+            Coalescence(collision_kernel=settings.kernel, adaptive=settings.adaptive),
+        ),
     )
     attributes = {}
     sampling = ConstantMultiplicity(settings.spectrum)
     attributes["volume"], attributes["multiplicity"] = sampling.sample_deterministic(
         settings.n_sd
     )
-    coalescence = Coalescence(
-        collision_kernel=settings.kernel, adaptive=settings.adaptive
-    )
-    builder.add_dynamic(coalescence)
     products = (
         ParticleVolumeVersusRadiusLogarithmSpectrum(
             settings.radius_bins_edges, name="dv/dlnr"

{pysdm_examples-3.0.0rc31 → pysdm_examples-3.0.0rc33}/PySDM_examples/Shipway_and_Hill_2012/simulation.py

@@ -1,4 +1,5 @@
 from collections import namedtuple
+from typing import List, Any
 
 import numpy as np
 from PySDM_examples.Shipway_and_Hill_2012.mpdata_1d import MPDATA_1D
@@ -63,15 +64,10 @@ class Simulation:
         )
         self.g_factor_vec = settings.rhod(_z_vec)
 
-        self.builder = Builder(
-            n_sd=settings.n_sd,
-            backend=backend(formulae=settings.formulae),
-            environment=env,
-        )
-        self.builder.add_dynamic(AmbientThermodynamics())
+        dynamics: List[Any] = [AmbientThermodynamics()]
 
         if settings.enable_condensation:
-            self.builder.add_dynamic(
+            dynamics.append(
                 Condensation(
                     adaptive=settings.condensation_adaptive,
                     rtol_thd=settings.condensation_rtol_thd,
@@ -79,19 +75,28 @@ class Simulation:
                     update_thd=settings.condensation_update_thd,
                 )
             )
-        self.builder.add_dynamic(EulerianAdvection(mpdata))
+        dynamics.append(EulerianAdvection(mpdata))
 
         self.products = []
         if settings.precip:
-            self.add_collision_dynamic(self.builder, settings, self.products)
+            self.add_collision_dynamic(dynamics, settings, self.products)
 
-        displacement = Displacement(
-            enable_sedimentation=settings.precip,
-            precipitation_counting_level_index=int(
-                settings.particle_reservoir_depth / settings.dz
-            ),
+        dynamics.append(
+            Displacement(
+                enable_sedimentation=settings.precip,
+                precipitation_counting_level_index=int(
+                    settings.particle_reservoir_depth / settings.dz
+                ),
+            )
         )
-        self.builder.add_dynamic(displacement)
+
+        self.builder = Builder(
+            n_sd=settings.n_sd,
+            backend=backend(formulae=settings.formulae),
+            environment=env,
+            dynamics=dynamics,
+        )
+
         self.attributes = self.builder.particulator.environment.init_attributes(
             spatial_discretisation=spatial_sampling.Pseudorandom(),
             spectral_discretisation=spectral_sampling.ConstantMultiplicity(
@@ -194,8 +199,8 @@ class Simulation:
                 )
 
     @staticmethod
-    def add_collision_dynamic(builder, settings, _):
-        builder.add_dynamic(
+    def add_collision_dynamic(dynamics, settings, _):
+        dynamics.append(
             Coalescence(
                 collision_kernel=settings.collision_kernel,
                 adaptive=settings.coalescence_adaptive,