PySDM 2.69__tar.gz → 2.71__tar.gz

{PySDM-2.69 → PySDM-2.71}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.69
+Version: 2.71
 Summary: Pythonic particle-based (super-droplet) warm-rain/aqueous-chemistry cloud microphysics package with box, parcel & 1D/2D prescribed-flow examples in Python, Julia and Matlab
 Home-page: https://github.com/open-atmos/PySDM
 Author: https://github.com/open-atmos/PySDM/graphs/contributors

{PySDM-2.69 → PySDM-2.71}/PySDM/builder.py

@@ -61,9 +61,9 @@ class Builder:
     def _register_product(self, product, buffer):
         if product.name in self.particulator.products:
             raise ValueError(f'product name "{product.name}" already registered')
-        product.set_buffer(buffer)
-        product.register(self)
-        self.particulator.products[product.name] = product
+        self.particulator.products[product.name] = product.instantiate(
+            builder=self, buffer=buffer
+        )
 
     def _resolve_attribute(self, attr_name):
         if attr_name not in self.req_attr:
@@ -121,8 +121,8 @@ class Builder:
             self._resolve_attribute(attr_name)
         self.req_attr_names = None
 
-        for dynamic in self.particulator.dynamics.values():
-            dynamic.register(self)
+        for key, dynamic in self.particulator.dynamics.items():
+            self.particulator.dynamics[key] = dynamic.instantiate(builder=self)
 
         single_buffer_for_all_products = np.empty(self.particulator.mesh.grid)
         for product in products:

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/ambient_thermodynamics.py

@@ -2,7 +2,10 @@
 environment-sync triggering class
 """
 
+from PySDM.dynamics.impl import register_dynamic
 
+
+@register_dynamic()
 class AmbientThermodynamics:
     def __init__(self):
         self.particulator = None

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/aqueous_chemistry.py

@@ -13,12 +13,14 @@ from PySDM.dynamics.impl.chemistry_utils import (
     M,
     SpecificGravities,
 )
+from PySDM.dynamics.impl import register_dynamic
 
 DEFAULTS = namedtuple("_", ("pH_min", "pH_max", "pH_rtol", "ionic_strength_threshold"))(
     pH_min=-1.0, pH_max=14.0, pH_rtol=1e-6, ionic_strength_threshold=0.02 * M
 )
 
 
+@register_dynamic()
 class AqueousChemistry:  # pylint: disable=too-many-instance-attributes
     def __init__(
         self,

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/collisions/collision.py

@@ -23,6 +23,7 @@ from PySDM.dynamics.impl.random_generator_optimizer_nopair import (
     RandomGeneratorOptimizerNoPair,
 )
 from PySDM.physics import si
+from PySDM.dynamics.impl import register_dynamic
 
 # pylint: disable=too-many-lines
 
@@ -36,6 +37,7 @@ DEFAULTS = namedtuple(
 )
 
 
+@register_dynamic()
 class Collision:  # pylint: disable=too-many-instance-attributes
     def __init__(
         self,
@@ -288,6 +290,7 @@ class Collision:  # pylint: disable=too-many-instance-attributes
         )
 
 
+@register_dynamic()
 class Coalescence(Collision):
     def __init__(
         self,
@@ -316,6 +319,7 @@ class Coalescence(Collision):
         )
 
 
+@register_dynamic()
 class Breakup(Collision):
     def __init__(
         self,

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/condensation.py

@@ -7,7 +7,8 @@ from collections import namedtuple
 
 import numpy as np
 
-from ..physics import si
+from PySDM.physics import si
+from PySDM.dynamics.impl import register_dynamic
 
 DEFAULTS = namedtuple("_", ("rtol_x", "rtol_thd", "cond_range", "schedule"))(
     rtol_x=1e-6,
@@ -17,6 +18,7 @@ DEFAULTS = namedtuple("_", ("rtol_x", "rtol_thd", "cond_range", "schedule"))(
 )
 
 
+@register_dynamic()
 class Condensation:  # pylint: disable=too-many-instance-attributes
     def __init__(
         self,

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/displacement.py

@@ -11,9 +11,12 @@ from collections import namedtuple
 
 import numpy as np
 
+from PySDM.dynamics.impl import register_dynamic
+
 DEFAULTS = namedtuple("_", ("rtol", "adaptive"))(rtol=1e-2, adaptive=True)
 
 
+@register_dynamic()
 class Displacement:  # pylint: disable=too-many-instance-attributes
     def __init__(
         self,

PySDM-2.71/PySDM/dynamics/eulerian_advection.py

@@ -0,0 +1,22 @@
+"""
+wrapper class for triggering integration in the Eulerian advection solver
+"""
+
+from PySDM.dynamics.impl import register_dynamic
+
+
+@register_dynamic()
+class EulerianAdvection:
+    def __init__(self, solvers):
+        self.solvers = solvers
+        self.particulator = None
+
+    def register(self, builder):
+        self.particulator = builder.particulator
+
+    def __call__(self):
+        for field in ("water_vapour_mixing_ratio", "thd"):
+            self.particulator.environment.get_predicted(field).download(
+                getattr(self.particulator.environment, f"get_{field}")(), reshape=True
+            )
+        self.solvers(self.particulator.dynamics["Displacement"])

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/freezing.py

@@ -7,8 +7,10 @@ from PySDM.backends.impl_common.freezing_attributes import (
     TimeDependentAttributes,
 )
 from PySDM.physics.heterogeneous_ice_nucleation_rate import Null
+from PySDM.dynamics.impl import register_dynamic
 
 
+@register_dynamic()
 class Freezing:
     def __init__(self, *, singular=True, record_freezing_temperature=False, thaw=False):
         assert not (record_freezing_temperature and singular)

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/impl/__init__.py

@@ -1 +1,3 @@
 """ stuff not intended to be imported from user code """
+
+from .register_dynamic import register_dynamic

PySDM-2.71/PySDM/dynamics/impl/register_dynamic.py

@@ -0,0 +1,21 @@
+""" decorator for dynamics classes
+ensuring that their instances can be re-used with multiple builders """
+
+from copy import deepcopy
+
+
+def _instantiate(self, *, builder):
+    copy = deepcopy(self)
+    copy.register(builder=builder)
+    return copy
+
+
+def register_dynamic():
+    def decorator(cls):
+        if hasattr(cls, "instantiate"):
+            assert cls.instantiate is _instantiate
+        else:
+            setattr(cls, "instantiate", _instantiate)
+        return cls
+
+    return decorator

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/isotopic_fractionation.py

@@ -4,11 +4,13 @@ requires condensation dynamic to be registered (and run beforehand)
 """
 
 from PySDM.dynamics.condensation import Condensation
+from PySDM.dynamics.impl import register_dynamic
 
 LIGHT_ISOTOPES = ("1H", "16O")
 HEAVY_ISOTOPES = ("2H", "3H", "17O", "18O")
 
 
+@register_dynamic()
 class IsotopicFractionation:
     def __init__(self, isotopes: tuple = HEAVY_ISOTOPES):
         self.isotopes = isotopes

{PySDM-2.69 → PySDM-2.71}/PySDM/dynamics/relaxed_velocity.py

@@ -6,8 +6,10 @@ towards the terminal velocity
 
 from PySDM.attributes.impl.attribute import Attribute
 from PySDM.particulator import Particulator
+from PySDM.dynamics.impl import register_dynamic
 
 
+@register_dynamic()
 class RelaxedVelocity:  # pylint: disable=too-many-instance-attributes
     """
     A dynamic which updates the fall momentum according to a relaxation timescale

{PySDM-2.69 → PySDM-2.71}/PySDM/environments/impl/moist.py

@@ -44,21 +44,21 @@ class Moist:
             result[var] = self.particulator.Storage.empty((self.mesh.n_cell,), float)
         return result
 
-    def __getitem__(self, index):
-        """returns a Storage representing the variable (field) at a given index or
-        otherwise a NaN-filled Storage if the index is not found (in order to simplify
+    def __getitem__(self, key: str):
+        """returns a Storage representing the variable (field) at a given key or
+        otherwise a NaN-filled Storage if the key is not found (in order to simplify
         generic code which uses optional variables, e.g. air viscosity, etc.)"""
-        if index in self._values["current"]:
-            return self._values["current"][index]
+        if key in self._values["current"]:
+            return self._values["current"][key]
         return self._nan_field
 
-    def get_predicted(self, index):
+    def get_predicted(self, key: str):
         if self._values["predicted"] is None:
             raise AssertionError(
                 "It seems the AmbientThermodynamics dynamic was not added"
                 " when building particulator"
             )
-        return self._values["predicted"][index]
+        return self._values["predicted"][key]
 
     def _recalculate_temperature_pressure_relative_humidity(self, target):
         self.particulator.backend.temperature_pressure_rh(

{PySDM-2.69 → PySDM-2.71}/PySDM/environments/kinematic_1d.py

@@ -26,7 +26,7 @@ class Kinematic1D(Moist):
         self._tmp["rhod"] = rhod
 
     def get_water_vapour_mixing_ratio(self) -> np.ndarray:
-        return self.particulator.dynamics["EulerianAdvection"].solvers.advectee.get()
+        return self.particulator.dynamics["EulerianAdvection"].solvers.advectee
 
     def get_thd(self) -> np.ndarray:
         return self.thd0
@@ -61,9 +61,7 @@ class Kinematic1D(Moist):
                 v_wet, n_per_kg = spectral_discretisation.sample(
                     backend=self.particulator.backend, n_sd=self.particulator.n_sd
                 )
-                # attributes["dry volume"] = v_wet
                 attributes["volume"] = v_wet
-                # attributes["kappa times dry volume"] = attributes["dry volume"] * kappa
             else:
                 r_dry, n_per_kg = spectral_discretisation.sample(
                     backend=self.particulator.backend, n_sd=self.particulator.n_sd

{PySDM-2.69 → PySDM-2.71}/PySDM/environments/kinematic_2d.py

@@ -85,16 +85,12 @@ class Kinematic2D(Moist):
         return attributes
 
     def get_thd(self):
-        return (
-            self.particulator.dynamics["EulerianAdvection"].solvers["th"].advectee.get()
-        )
+        return self.particulator.dynamics["EulerianAdvection"].solvers["th"]
 
     def get_water_vapour_mixing_ratio(self):
-        return (
-            self.particulator.dynamics["EulerianAdvection"]
-            .solvers["water_vapour_mixing_ratio"]
-            .advectee.get()
-        )
+        return self.particulator.dynamics["EulerianAdvection"].solvers[
+            "water_vapour_mixing_ratio"
+        ]
 
     def sync(self):
         self.particulator.dynamics["EulerianAdvection"].solvers.wait()

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_dry_air_density.py

@@ -2,9 +2,10 @@
 ambient dry-air density
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientDryAirDensity(MoistEnvironmentProduct):
     def __init__(self, name="rhod", unit="kg/m^3", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_dry_air_potential_temperature.py

@@ -2,9 +2,10 @@
 ambient dry-air potential temperature (computed using dry air partial pressure)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientDryAirPotentialTemperature(MoistEnvironmentProduct):
     def __init__(self, unit="K", name=None, var=None):
         super().__init__(unit=unit, name=name, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_pressure.py

@@ -2,9 +2,10 @@
 ambient pressure
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientPressure(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="Pa", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_relative_humidity.py

@@ -2,9 +2,10 @@
 ambient relative humidity (wrt water or ice)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientRelativeHumidity(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="dimensionless", var=None, ice=False):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_temperature.py

@@ -2,9 +2,10 @@
 ambient temperature
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientTemperature(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="K", var=None):
         super().__init__(name=name, unit=unit, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/ambient_thermodynamics/ambient_water_vapour_mixing_ratio.py

@@ -2,9 +2,10 @@
 ambient water vapour mixing ratio (mass of water vapour per mass of dry air)
 """
 
-from PySDM.products.impl.moist_environment_product import MoistEnvironmentProduct
+from PySDM.products.impl import MoistEnvironmentProduct, register_product
 
 
+@register_product()
 class AmbientWaterVapourMixingRatio(MoistEnvironmentProduct):
     def __init__(self, name=None, unit="dimensionless", var=None):
         super().__init__(unit=unit, name=name, var=var)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/aqueous_chemistry/acidity.py

@@ -5,9 +5,10 @@ with weighting either by number or volume
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class Acidity(MomentProduct):
     def __init__(
         self,

{PySDM-2.69 → PySDM-2.71}/PySDM/products/aqueous_chemistry/aqueous_mass_spectrum.py

@@ -8,9 +8,10 @@ from chempy import Substance
 
 from PySDM.dynamics.impl.chemistry_utils import AQUEOUS_COMPOUNDS
 from PySDM.physics.constants import si
-from PySDM.products.impl.spectrum_moment_product import SpectrumMomentProduct
+from PySDM.products.impl import SpectrumMomentProduct, register_product
 
 
+@register_product()
 class AqueousMassSpectrum(SpectrumMomentProduct):
     def __init__(
         self, *, key, dry_radius_bins_edges, specific=False, name=None, unit="kg/m^3"
@@ -58,6 +59,7 @@ class AqueousMassSpectrum(SpectrumMomentProduct):
         return vals
 
 
+@register_product()
 class SpecificAqueousMassSpectrum(AqueousMassSpectrum):
     def __init__(self, key, dry_radius_bins_edges, name=None, unit="dimensionless"):
         super().__init__(

{PySDM-2.69 → PySDM-2.71}/PySDM/products/aqueous_chemistry/aqueous_mole_fraction.py

@@ -2,11 +2,12 @@
 mole fractions of aqueous-chemistry relevant compounds
 """
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 DUMMY_SPECIFIC_GRAVITY = 44
 
 
+@register_product()
 class AqueousMoleFraction(MomentProduct):
     def __init__(self, key, unit="dimensionless", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/aqueous_chemistry/gaseous_mole_fraction.py

@@ -3,9 +3,10 @@ mole fractions of gaseous compounds relevant for aqueous chemistry
 """
 
 from PySDM.dynamics.impl.chemistry_utils import GASEOUS_COMPOUNDS
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class GaseousMoleFraction(Product):
     def __init__(self, key, unit="dimensionless", name=None):
         super().__init__(name=name, unit=unit)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/aqueous_chemistry/total_dry_mass_mixing_ratio.py

@@ -4,9 +4,10 @@ dry aerosol mass summed over all particles in a grid box per mass of dry air
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class TotalDryMassMixingRatio(MomentProduct):
     def __init__(self, density, name=None, unit="kg/kg"):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/collision/collision_rates.py

@@ -4,9 +4,10 @@ Rates of collision events and their deficit wrt expected values in case mutiplic
  or due to unrepresentable breakups from integer overflow
 """
 
-from PySDM.products.impl.rate_product import RateProduct
+from PySDM.products.impl import RateProduct, register_product
 
 
+@register_product()
 class CollisionRateDeficitPerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -14,6 +15,7 @@ class CollisionRateDeficitPerGridbox(RateProduct):
         )
 
 
+@register_product()
 class CollisionRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -21,6 +23,7 @@ class CollisionRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class CoalescenceRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -28,6 +31,7 @@ class CoalescenceRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class BreakupRatePerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(
@@ -35,6 +39,7 @@ class BreakupRatePerGridbox(RateProduct):
         )
 
 
+@register_product()
 class BreakupRateDeficitPerGridbox(RateProduct):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__(

{PySDM-2.69 → PySDM-2.71}/PySDM/products/collision/collision_timestep_mean.py

@@ -7,9 +7,10 @@ import numba
 import numpy as np
 
 from PySDM.backends.impl_numba.conf import JIT_FLAGS
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class CollisionTimestepMean(Product):
     def __init__(self, unit="s", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/collision/collision_timestep_min.py

@@ -3,9 +3,10 @@ Minimal collision timestep used when adaptive timestepping is enabled in the
  `PySDM.dynamics.collisions.collision.Collision` dynamic (fetching a value resets the counter)
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class CollisionTimestepMin(Product):
     def __init__(self, unit="s", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/condensation/activable_fraction.py

@@ -3,9 +3,10 @@ fraction of particles with critical supersaturation lower than a given supersatu
  (passed as keyword argument while calling `get()`)
 """
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class ActivableFraction(MomentProduct):
     def __init__(self, unit="dimensionless", name=None):
         super().__init__(name=name, unit=unit)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/condensation/condensation_timestep.py

@@ -4,9 +4,10 @@ minimum and maximum condensation timestep (fetching a value resets the counter)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class _CondensationTimestep(Product):
     def __init__(self, name, unit, extremum, reset_value):
         super().__init__(
@@ -37,11 +38,13 @@ class _CondensationTimestep(Product):
         return self.buffer
 
 
+@register_product()
 class CondensationTimestepMin(_CondensationTimestep):
     def __init__(self, name=None, unit="s"):
         super().__init__(name=name, unit=unit, extremum=np.minimum, reset_value=np.inf)
 
 
+@register_product()
 class CondensationTimestepMax(_CondensationTimestep):
     def __init__(self, name=None, unit="s"):
         super().__init__(name=name, unit=unit, extremum=np.maximum, reset_value=-np.inf)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/condensation/event_rates.py

@@ -5,9 +5,10 @@ rates of activation, deactivation and ripening events (take into account substep
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class EventRate(Product):
     def __init__(self, what, name=None, unit=None):
         super().__init__(name=name, unit=unit)
@@ -42,16 +43,19 @@ class EventRate(Product):
         return self.buffer
 
 
+@register_product()
 class RipeningRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("ripening", name=name, unit=unit)
 
 
+@register_product()
 class ActivatingRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("activating", name=name, unit=unit)
 
 
+@register_product()
 class DeactivatingRate(EventRate):
     def __init__(self, name=None, unit="s^-1 kg^-1"):
         super().__init__("deactivating", name=name, unit=unit)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/condensation/peak_supersaturation.py

@@ -7,9 +7,10 @@ highest supersaturation encountered while solving for condensation/evaporation (
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class PeakSupersaturation(Product):
     def __init__(self, unit="dimensionless", name=None):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/displacement/averaged_terminal_velocity.py

@@ -4,9 +4,10 @@ average terminal with weighting either by number or volume
 
 import numpy as np
 
-from PySDM.products.impl.moment_product import MomentProduct
+from PySDM.products.impl import MomentProduct, register_product
 
 
+@register_product()
 class AveragedTerminalVelocity(MomentProduct):
     def __init__(
         self, *, radius_range=(0, np.inf), weighting="volume", unit="m/s", name=None

{PySDM-2.69 → PySDM-2.71}/PySDM/products/displacement/flow_velocity_component.py

@@ -5,9 +5,10 @@ absolute values of Courant number on all edges of a cell)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class FlowVelocityComponent(Product):
     def __init__(self, component: int, name=None, unit="m/s"):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/displacement/max_courant_number.py

@@ -5,9 +5,10 @@ absolute values of Courant number on all edges of a cell)
 
 import numpy as np
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class MaxCourantNumber(Product):
     def __init__(self, name=None, unit="dimensionless"):
         super().__init__(unit=unit, name=name)

{PySDM-2.69 → PySDM-2.71}/PySDM/products/displacement/surface_precipitation.py

@@ -2,9 +2,10 @@
 water volume flux derived from sizes of particles crossing bottom domain boundary
 """
 
-from PySDM.products.impl.product import Product
+from PySDM.products.impl import Product, register_product
 
 
+@register_product()
 class SurfacePrecipitation(Product):
     def __init__(self, name=None, unit="m/s"):
         super().__init__(unit=unit, name=name)