PySDM-examples 2.112__py3-none-any.whl → 2.113__py3-none-any.whl

PySDM_examples/Abade_and_Albuquerque_2024/settings.py

@@ -14,6 +14,7 @@ class Settings:
         timestep: float,
         enable_immersion_freezing: bool = True,
         enable_vapour_deposition_on_ice: bool = True,
+        inp_frac: float = 0.1,
     ):
         self.n_sd = n_sd
         self.timestep = timestep
@@ -24,14 +25,21 @@ class Settings:
 
         # parameters from the paper
         self.formulae = Formulae(
-            constants={"bulk_phase_partitioning_exponent": 0.1},
+            constants={
+                "bulk_phase_partitioning_exponent": 0.1,
+                "NIEMAND_A": -0.517,
+                "NIEMAND_B": 8.934,
+            },
             bulk_phase_partitioning="KaulEtAl2015",
             particle_shape_and_density="MixedPhaseSpheres",
             diffusion_coordinate="WaterMassLogarithm",
+            freezing_temperature_spectrum="Niemand_et_al_2012",
         )
         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
+        self.freezing_inp_frac = inp_frac
+        self.freezing_inp_dry_radius = 0.5 * si.um
 
         thd_0 = self.formulae.state_variable_triplet.th_dry(
             th_std=269 * si.K,

PySDM_examples/Abade_and_Albuquerque_2024/simulation.py

@@ -57,10 +57,24 @@ class Simulation(BasicSimulation):
         )
         del attributes["volume"]
 
-        # TODO #1524
         if settings.enable_immersion_freezing:
-            attributes["freezing temperature"] = np.full(
-                shape=(settings.n_sd,), fill_value=250
+            trivia = builder.particulator.formulae.trivia
+            n_inp = int(settings.n_sd * settings.freezing_inp_frac)
+
+            rng = np.random.default_rng(seed=builder.particulator.formulae.seed)
+
+            attributes["freezing temperature"] = rng.permutation(
+                np.pad(
+                    builder.particulator.formulae.freezing_temperature_spectrum.invcdf(
+                        cdf=rng.uniform(low=0, high=1, size=n_inp),
+                        A_insol=trivia.sphere_surface(
+                            diameter=2 * settings.freezing_inp_dry_radius
+                        ),
+                    ),
+                    (0, settings.n_sd - n_inp),
+                    mode="constant",
+                    constant_values=(trivia.C2K(-38)),
+                )
             )
 
         self.products = (

{pysdm_examples-2.112.dist-info → pysdm_examples-2.113.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: pysdm-examples
-Version: 2.112
+Version: 2.113
 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-2.112.dist-info → pysdm_examples-2.113.dist-info}/RECORD

@@ -1,7 +1,7 @@
 PySDM_examples/__init__.py,sha256=TnieTsOWPG-6J5Z9S_LgBRil-kbBWmsl3q2dlA4equk,228
 PySDM_examples/Abade_and_Albuquerque_2024/__init__.py,sha256=JlXdzgcWzQEj87r0IatZAHx4j3hgDsyxP8LisaJvP_U,279
-PySDM_examples/Abade_and_Albuquerque_2024/settings.py,sha256=KVll9MY-EEYbxCrDs8PDasHg1YLMhLVIO634WXeb9sQ,2015
-PySDM_examples/Abade_and_Albuquerque_2024/simulation.py,sha256=mbUH6mr76TnXDsjhc20-8dBpD1uGuW0IuD7WSW7OUHM,2771
+PySDM_examples/Abade_and_Albuquerque_2024/settings.py,sha256=GtHwC6ZnHy9yEl4vNS2PdBQK6jxm_boAlNnych93OOQ,2307
+PySDM_examples/Abade_and_Albuquerque_2024/simulation.py,sha256=9lBU1M4nKIn1zJFhXsJbdms2Nzasb-98F_a9HCvwccE,3429
 PySDM_examples/Abdul_Razzak_Ghan_2000/__init__.py,sha256=2M2S0Qff2PJhxzEiV7eIGhJHCVkUywM0moIBT_BXJIg,311
 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
@@ -154,7 +154,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.112.dist-info/METADATA,sha256=Dc2BNvc42BXVtcaIB7ZNKD1TJMRRNh-qAG1_cdwfC-A,12686
-pysdm_examples-2.112.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
-pysdm_examples-2.112.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
-pysdm_examples-2.112.dist-info/RECORD,,
+pysdm_examples-2.113.dist-info/METADATA,sha256=D9MY9yOy4EwDTXJhse7wKgP3xTmTdMdqPDtrIUWSlZQ,12686
+pysdm_examples-2.113.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
+pysdm_examples-2.113.dist-info/top_level.txt,sha256=axw2inuuYKPWkQhzBYtkZJTPDlo5LrpGLOA0iSAf9v0,15
+pysdm_examples-2.113.dist-info/RECORD,,