PySDM 2.85__tar.gz → 2.87__tar.gz

{pysdm-2.85 → pysdm-2.87}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.85
+Version: 2.87
 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
@@ -75,7 +75,7 @@ The package features a Pythonic high-performance implementation of the
   Super-Droplet Method (SDM) Monte-Carlo algorithm for representing collisional growth 
   ([Shima et al. 2009](https://rmets.onlinelibrary.wiley.com/doi/abs/10.1002/qj.441)), hence the name. 
 
-There is a growing set of example Jupyter notebooks exemplifying how to perform 
+There is a growing set of [example Jupyter notebooks](https://open-atmos.github.io/PySDM/PySDM_examples.html) exemplifying how to perform 
   various types of calculations and simulations using PySDM.
 Most of the example notebooks reproduce results and plot from literature, see below for 
   a list of examples and links to the notebooks (which can be either executed or viewed 
@@ -134,7 +134,7 @@ Conda users may use ``pip`` as well, see the [Installing non-conda packages](htt
 For development purposes, we suggest cloning the repository and installing it using ``pip -e``.
 Test-time dependencies can be installed with ``pip -e .[tests]``.
 
-PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM-examples) package.
+PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM/blob/main/examples) package.
 The examples have additional dependencies listed in [``PySDM_examples`` package ``setup.py``](https://github.com/open-atmos/PySDM/blob/main/examples/setup.py) file.
 Running the example Jupyter notebooks requires the ``PySDM_examples`` package to be installed.
 The suggested install and launch steps are:
@@ -189,7 +189,8 @@ To run all pre-commit hooks, run `pre-commit run --all-files`.
 The `.pre-commit-config.yaml` file can be modified in case new hooks are to be added or
   existing ones need to be altered.  
 
-Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki).
+Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki)
+  and in [PySDM HOWTOs](https://github.com/open-atmos/PySDM/tree/main/examples/PySDM_examples/_HOWTOS).
 
 Issues regarding any incorrect, unintuitive or undocumented bahaviour of
 PySDM are best to be reported on the [GitHub issue tracker](https://github.com/open-atmos/PySDM/issues/new).

{pysdm-2.85 → pysdm-2.87}/PySDM/physics/constants_defaults.py

@@ -428,3 +428,12 @@ one_kelvin = 1 * si.K
 bulk_phase_partitioning_T_cold = 235 * si.K
 bulk_phase_partitioning_T_warm = 273 * si.K
 bulk_phase_partitioning_exponent = np.nan
+
+
+BOLIN_ISOTOPE_TIMESCALE_COEFF_C1 = np.nan * si.dimensionless
+"""
+Coeffitient c1 used in [Bolin 1958](https://https://digitallibrary.un.org/record/3892725)
+for the falling drop evaporation timescale of equilibration with ambient air void of a given
+isotopologue; in the paper timescale is calculated for tritium with assumption of no tritium
+in the environment around the drop (Table 1).
+"""

{pysdm-2.85 → pysdm-2.87}/PySDM/physics/dimensional_analysis.py

@@ -1,12 +1,12 @@
 """
-A context manager (for use with the `with` statement)
-for use in unit tests which disables Numba and enables Pint
+A context manager (for use with the `with` statement) which
+enables Pint physical-units checks and disables Numba in `PySDM.formulae.Formulae`
 """
 
 from importlib import reload
 
 from PySDM import formulae
-
+from PySDM import physics
 from . import constants, constants_defaults
 from .impl import flag
 
@@ -17,9 +17,11 @@ class DimensionalAnalysis:
         reload(constants)
         reload(constants_defaults)
         reload(formulae)
+        reload(physics)
 
     def __exit__(*_):  # pylint: disable=no-method-argument,no-self-argument
         flag.DIMENSIONAL_ANALYSIS = False
         reload(constants)
         reload(constants_defaults)
         reload(formulae)
+        reload(physics)

{pysdm-2.85 → pysdm-2.87}/PySDM/physics/isotope_relaxation_timescale/__init__.py

@@ -4,3 +4,4 @@ isotope relaxation timescale formulae
 
 from PySDM.impl.null_physics_class import Null
 from .miyake_et_al_1968 import MiyakeEtAl1968
+from .bolin_1958 import Bolin1958

pysdm-2.87/PySDM/physics/isotope_relaxation_timescale/bolin_1958.py

@@ -0,0 +1,17 @@
+""" Timescale derived for tritium with assumption of zero ambient concentration - see text above
+    Table 1 [Bolin 1958](https://digitallibrary.un.org/record/3892725) """
+
+import numpy as np
+
+
+class Bolin1958:  # pylint: disable=too-few-public-methods
+    """Implementation of timescale used by Bolin 1958"""
+
+    def __init__(self, const):
+        assert np.isfinite(const.BOLIN_ISOTOPE_TIMESCALE_COEFF_C1)
+
+    @staticmethod
+    # pylint: disable=too-many-arguments
+    def tau(const, radius, r_dr_dt):
+        """timescale for evaporation of a falling drop with tritium"""
+        return (-3 / radius**2 * r_dr_dt * const.BOLIN_ISOTOPE_TIMESCALE_COEFF_C1) ** -1

{pysdm-2.85 → pysdm-2.87}/PySDM.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PySDM
-Version: 2.85
+Version: 2.87
 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
@@ -75,7 +75,7 @@ The package features a Pythonic high-performance implementation of the
   Super-Droplet Method (SDM) Monte-Carlo algorithm for representing collisional growth 
   ([Shima et al. 2009](https://rmets.onlinelibrary.wiley.com/doi/abs/10.1002/qj.441)), hence the name. 
 
-There is a growing set of example Jupyter notebooks exemplifying how to perform 
+There is a growing set of [example Jupyter notebooks](https://open-atmos.github.io/PySDM/PySDM_examples.html) exemplifying how to perform 
   various types of calculations and simulations using PySDM.
 Most of the example notebooks reproduce results and plot from literature, see below for 
   a list of examples and links to the notebooks (which can be either executed or viewed 
@@ -134,7 +134,7 @@ Conda users may use ``pip`` as well, see the [Installing non-conda packages](htt
 For development purposes, we suggest cloning the repository and installing it using ``pip -e``.
 Test-time dependencies can be installed with ``pip -e .[tests]``.
 
-PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM-examples) package.
+PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM/blob/main/examples) package.
 The examples have additional dependencies listed in [``PySDM_examples`` package ``setup.py``](https://github.com/open-atmos/PySDM/blob/main/examples/setup.py) file.
 Running the example Jupyter notebooks requires the ``PySDM_examples`` package to be installed.
 The suggested install and launch steps are:
@@ -189,7 +189,8 @@ To run all pre-commit hooks, run `pre-commit run --all-files`.
 The `.pre-commit-config.yaml` file can be modified in case new hooks are to be added or
   existing ones need to be altered.  
 
-Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki).
+Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki)
+  and in [PySDM HOWTOs](https://github.com/open-atmos/PySDM/tree/main/examples/PySDM_examples/_HOWTOS).
 
 Issues regarding any incorrect, unintuitive or undocumented bahaviour of
 PySDM are best to be reported on the [GitHub issue tracker](https://github.com/open-atmos/PySDM/issues/new).

{pysdm-2.85 → pysdm-2.87}/PySDM.egg-info/SOURCES.txt

@@ -299,6 +299,7 @@ PySDM/physics/isotope_ratio_evolution/gedzelman_and_arnold_1994.py
 PySDM/physics/isotope_ratio_evolution/merlivat_and_jouzel_1979.py
 PySDM/physics/isotope_ratio_evolution/rayleigh_distillation.py
 PySDM/physics/isotope_relaxation_timescale/__init__.py
+PySDM/physics/isotope_relaxation_timescale/bolin_1958.py
 PySDM/physics/isotope_relaxation_timescale/miyake_et_al_1968.py
 PySDM/physics/latent_heat/__init__.py
 PySDM/physics/latent_heat/constant.py
@@ -411,8 +412,11 @@ PySDM/products/size_spectral/total_particle_concentration.py
 PySDM/products/size_spectral/total_particle_specific_concentration.py
 PySDM/products/size_spectral/water_mixing_ratio.py
 docs/markdown/pysdm_landing.md
+docs/templates/README.md
 docs/templates/custom.css
 docs/templates/index.html.jinja2
+docs/templates/syntax-highlighting.css
+docs/templates/theme.css
 examples/MANIFEST.in
 examples/README.md
 examples/pyproject.toml
@@ -483,6 +487,8 @@ examples/PySDM_examples/Bieli_et_al_2022/__init__.py
 examples/PySDM_examples/Bieli_et_al_2022/make_fig_3.ipynb
 examples/PySDM_examples/Bieli_et_al_2022/settings.py
 examples/PySDM_examples/Bieli_et_al_2022/simulation.py
+examples/PySDM_examples/Bolin_1958/__init__.py
+examples/PySDM_examples/Bolin_1958/table_1.ipynb
 examples/PySDM_examples/Bolot_et_al_2013/__init__.py
 examples/PySDM_examples/Bolot_et_al_2013/fig_1.ipynb
 examples/PySDM_examples/Bulenok_2023_MasterThesis/__init__.py
@@ -606,6 +612,7 @@ examples/PySDM_examples/Yang_et_al_2018/__init__.py
 examples/PySDM_examples/Yang_et_al_2018/fig_2.ipynb
 examples/PySDM_examples/Yang_et_al_2018/settings.py
 examples/PySDM_examples/Yang_et_al_2018/simulation.py
+examples/PySDM_examples/_HOWTOs/dimensional_analysis.ipynb
 examples/PySDM_examples/deJong_Azimi/__init__.py
 examples/PySDM_examples/deJong_Azimi/box.ipynb
 examples/PySDM_examples/deJong_Azimi/cloudy_data.py
@@ -692,6 +699,7 @@ tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_freezing.py
 tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_gui_settings.py
 tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_initialisation.py
 tests/smoke_tests/kinematic_2d/arabas_et_al_2015/test_spin_up.py
+tests/smoke_tests/no_env/bolin_1958/test_table_1.py
 tests/smoke_tests/no_env/gedzelman_and_arnold_1994/__init__.py
 tests/smoke_tests/no_env/gedzelman_and_arnold_1994/test_fig_2.py
 tests/smoke_tests/no_env/lamb_et_al_2017/__init__.py

{pysdm-2.85 → pysdm-2.87}/README.md

@@ -35,7 +35,7 @@ The package features a Pythonic high-performance implementation of the
   Super-Droplet Method (SDM) Monte-Carlo algorithm for representing collisional growth 
   ([Shima et al. 2009](https://rmets.onlinelibrary.wiley.com/doi/abs/10.1002/qj.441)), hence the name. 
 
-There is a growing set of example Jupyter notebooks exemplifying how to perform 
+There is a growing set of [example Jupyter notebooks](https://open-atmos.github.io/PySDM/PySDM_examples.html) exemplifying how to perform 
   various types of calculations and simulations using PySDM.
 Most of the example notebooks reproduce results and plot from literature, see below for 
   a list of examples and links to the notebooks (which can be either executed or viewed 
@@ -94,7 +94,7 @@ Conda users may use ``pip`` as well, see the [Installing non-conda packages](htt
 For development purposes, we suggest cloning the repository and installing it using ``pip -e``.
 Test-time dependencies can be installed with ``pip -e .[tests]``.
 
-PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM-examples) package.
+PySDM examples constitute the [``PySDM-examples``](https://github.com/open-atmos/PySDM/blob/main/examples) package.
 The examples have additional dependencies listed in [``PySDM_examples`` package ``setup.py``](https://github.com/open-atmos/PySDM/blob/main/examples/setup.py) file.
 Running the example Jupyter notebooks requires the ``PySDM_examples`` package to be installed.
 The suggested install and launch steps are:
@@ -149,7 +149,8 @@ To run all pre-commit hooks, run `pre-commit run --all-files`.
 The `.pre-commit-config.yaml` file can be modified in case new hooks are to be added or
   existing ones need to be altered.  
 
-Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki).
+Further hints addressed at PySDM developers are maintained in the [open-atmos/python-dev-hints Wiki](https://github.com/open-atmos/python-dev-hints/wiki)
+  and in [PySDM HOWTOs](https://github.com/open-atmos/PySDM/tree/main/examples/PySDM_examples/_HOWTOS).
 
 Issues regarding any incorrect, unintuitive or undocumented bahaviour of
 PySDM are best to be reported on the [GitHub issue tracker](https://github.com/open-atmos/PySDM/issues/new).

pysdm-2.87/docs/templates/README.md

@@ -0,0 +1 @@
+dark mode .css from https://github.com/mitmproxy/pdoc/tree/main/examples/dark-mode

pysdm-2.87/docs/templates/custom.css

@@ -0,0 +1,9 @@
+img {
+    max-width: 100%;
+}
+
+mark {
+    background-color: #859bed;
+    color: black;
+}
+

{pysdm-2.85 → pysdm-2.87}/docs/templates/index.html.jinja2

@@ -91,7 +91,8 @@
         <pre><code>git clone https://github.com/open-atmos/PySDM.git
 pip install -e PySDM[tests] -e PySDM/examples[tests]
 pytest PySDM</code></pre>
-        <p>(the above should be a viable way to set up development environment for PySDM, see also our <a href="https://github.com/open-atmos/python-dev-hints/wiki">Python dev hints Wiki</a> for further information)</p>
+        <p>(the above should be a viable way to set up development environment for PySDM, see also our <a href="https://github.com/open-atmos/python-dev-hints/wiki">Python dev hints Wiki</a> 
+          and <a href="https://github.com/open-atmos/PySDM/tree/main/examples/PySDM_examples/_HOWTOS">PySDM HOWTOs</a> for further information)</p>
         <p>
             PySDM-examples is also available on PyPI and can be installed using pip:
         </p>

pysdm-2.87/docs/templates/syntax-highlighting.css

@@ -0,0 +1,80 @@
+/* monokai color scheme, see pdoc/template/README.md */
+pre { line-height: 125%; }
+span.linenos { color: inherit; background-color: transparent; padding-left: 5px; padding-right: 20px; }
+.pdoc-code .hll { background-color: #49483e }
+.pdoc-code { background: #272822; color: #f8f8f2 }
+.pdoc-code .c { color: #75715e } /* Comment */
+.pdoc-code .err { color: #960050; background-color: #1e0010 } /* Error */
+.pdoc-code .esc { color: #f8f8f2 } /* Escape */
+.pdoc-code .g { color: #f8f8f2 } /* Generic */
+.pdoc-code .k { color: #66d9ef } /* Keyword */
+.pdoc-code .l { color: #ae81ff } /* Literal */
+.pdoc-code .n { color: #f8f8f2 } /* Name */
+.pdoc-code .o { color: #f92672 } /* Operator */
+.pdoc-code .x { color: #f8f8f2 } /* Other */
+.pdoc-code .p { color: #f8f8f2 } /* Punctuation */
+.pdoc-code .ch { color: #75715e } /* Comment.Hashbang */
+.pdoc-code .cm { color: #75715e } /* Comment.Multiline */
+.pdoc-code .cp { color: #75715e } /* Comment.Preproc */
+.pdoc-code .cpf { color: #75715e } /* Comment.PreprocFile */
+.pdoc-code .c1 { color: #75715e } /* Comment.Single */
+.pdoc-code .cs { color: #75715e } /* Comment.Special */
+.pdoc-code .gd { color: #f92672 } /* Generic.Deleted */
+.pdoc-code .ge { color: #f8f8f2; font-style: italic } /* Generic.Emph */
+.pdoc-code .gr { color: #f8f8f2 } /* Generic.Error */
+.pdoc-code .gh { color: #f8f8f2 } /* Generic.Heading */
+.pdoc-code .gi { color: #a6e22e } /* Generic.Inserted */
+.pdoc-code .go { color: #66d9ef } /* Generic.Output */
+.pdoc-code .gp { color: #f92672; font-weight: bold } /* Generic.Prompt */
+.pdoc-code .gs { color: #f8f8f2; font-weight: bold } /* Generic.Strong */
+.pdoc-code .gu { color: #75715e } /* Generic.Subheading */
+.pdoc-code .gt { color: #f8f8f2 } /* Generic.Traceback */
+.pdoc-code .kc { color: #66d9ef } /* Keyword.Constant */
+.pdoc-code .kd { color: #66d9ef } /* Keyword.Declaration */
+.pdoc-code .kn { color: #f92672 } /* Keyword.Namespace */
+.pdoc-code .kp { color: #66d9ef } /* Keyword.Pseudo */
+.pdoc-code .kr { color: #66d9ef } /* Keyword.Reserved */
+.pdoc-code .kt { color: #66d9ef } /* Keyword.Type */
+.pdoc-code .ld { color: #e6db74 } /* Literal.Date */
+.pdoc-code .m { color: #ae81ff } /* Literal.Number */
+.pdoc-code .s { color: #e6db74 } /* Literal.String */
+.pdoc-code .na { color: #a6e22e } /* Name.Attribute */
+.pdoc-code .nb { color: #f8f8f2 } /* Name.Builtin */
+.pdoc-code .nc { color: #a6e22e } /* Name.Class */
+.pdoc-code .no { color: #66d9ef } /* Name.Constant */
+.pdoc-code .nd { color: #a6e22e } /* Name.Decorator */
+.pdoc-code .ni { color: #f8f8f2 } /* Name.Entity */
+.pdoc-code .ne { color: #a6e22e } /* Name.Exception */
+.pdoc-code .nf { color: #a6e22e } /* Name.Function */
+.pdoc-code .nl { color: #f8f8f2 } /* Name.Label */
+.pdoc-code .nn { color: #f8f8f2 } /* Name.Namespace */
+.pdoc-code .nx { color: #a6e22e } /* Name.Other */
+.pdoc-code .py { color: #f8f8f2 } /* Name.Property */
+.pdoc-code .nt { color: #f92672 } /* Name.Tag */
+.pdoc-code .nv { color: #f8f8f2 } /* Name.Variable */
+.pdoc-code .ow { color: #f92672 } /* Operator.Word */
+.pdoc-code .w { color: #f8f8f2 } /* Text.Whitespace */
+.pdoc-code .mb { color: #ae81ff } /* Literal.Number.Bin */
+.pdoc-code .mf { color: #ae81ff } /* Literal.Number.Float */
+.pdoc-code .mh { color: #ae81ff } /* Literal.Number.Hex */
+.pdoc-code .mi { color: #ae81ff } /* Literal.Number.Integer */
+.pdoc-code .mo { color: #ae81ff } /* Literal.Number.Oct */
+.pdoc-code .sa { color: #e6db74 } /* Literal.String.Affix */
+.pdoc-code .sb { color: #e6db74 } /* Literal.String.Backtick */
+.pdoc-code .sc { color: #e6db74 } /* Literal.String.Char */
+.pdoc-code .dl { color: #e6db74 } /* Literal.String.Delimiter */
+.pdoc-code .sd { color: #e6db74 } /* Literal.String.Doc */
+.pdoc-code .s2 { color: #e6db74 } /* Literal.String.Double */
+.pdoc-code .se { color: #ae81ff } /* Literal.String.Escape */
+.pdoc-code .sh { color: #e6db74 } /* Literal.String.Heredoc */
+.pdoc-code .si { color: #e6db74 } /* Literal.String.Interpol */
+.pdoc-code .sx { color: #e6db74 } /* Literal.String.Other */
+.pdoc-code .sr { color: #e6db74 } /* Literal.String.Regex */
+.pdoc-code .s1 { color: #e6db74 } /* Literal.String.Single */
+.pdoc-code .ss { color: #e6db74 } /* Literal.String.Symbol */
+.pdoc-code .bp { color: #f8f8f2 } /* Name.Builtin.Pseudo */
+.pdoc-code .fm { color: #a6e22e } /* Name.Function.Magic */
+.pdoc-code .vc { color: #f8f8f2 } /* Name.Variable.Class */
+.pdoc-code .vg { color: #f8f8f2 } /* Name.Variable.Global */
+.pdoc-code .vi { color: #f8f8f2 } /* Name.Variable.Instance */
+.pdoc-code .vm { color: #f8f8f2 } /* Name.Variable.Magic */

pysdm-2.87/docs/templates/theme.css

@@ -0,0 +1,20 @@
+:root {
+    --pdoc-background: #212529;
+}
+
+.pdoc {
+    --text: #f7f7f7;
+    --muted: #9d9d9d;
+    --link: #58a6ff;
+    --link-hover: #3989ff;
+    --code: #333;
+    --active: #555;
+
+    --accent: #343434;
+    --accent2: #555;
+
+    --nav-hover: rgba(0, 0, 0, 0.1);
+    --name: #77C1FF;
+    --def: #0cdd0c;
+    --annotation: #00c037;
+}

{pysdm-2.85 → pysdm-2.87}/examples/PySDM_examples/Arabas_et_al_2015/settings.py

@@ -1,12 +1,8 @@
 from typing import Iterable
 
-import numpy as np
 from PySDM_examples.Morrison_and_Grabowski_2007.strato_cumulus import StratoCumulus
 
 from PySDM import Formulae
-from PySDM.dynamics.collisions.breakup_efficiencies import ConstEb
-from PySDM.dynamics.collisions.breakup_fragmentations import Gaussian
-from PySDM.dynamics.collisions.coalescence_efficiencies import ConstEc
 from PySDM.physics import si
 
 
@@ -39,12 +35,3 @@ class Settings(StratoCumulus):
         self.simulation_time = 90 * si.minute
         self.dt = 5 * si.second
         self.spin_up_time = 1 * si.hour
-
-        # additional breakup dynamics
-        mu_r = 10 * si.um
-        mu = 4 / 3 * np.pi * mu_r**3
-        sigma = mu / 2.5
-        vmin = mu / 1000
-        self.coalescence_efficiency = ConstEc(Ec=0.95)
-        self.breakup_efficiency = ConstEb(Eb=1.0)
-        self.breakup_fragmentation = Gaussian(mu=mu, sigma=sigma, vmin=vmin, nfmax=10)

pysdm-2.87/examples/PySDM_examples/Bolin_1958/table_1.ipynb

@@ -0,0 +1,267 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "e3f1edc815e974f4",
+   "metadata": {},
+   "source": [
+    "[![preview notebook](https://img.shields.io/static/v1?label=render%20on&logo=github&color=87ce3e&message=GitHub)](https://github.com/open-atmos/devops_tests/blob/main/test_files/template.ipynb)\n",
+    "[![launch on mybinder.org](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/open-atmos/devops_tests.git/main?urlpath=lab/tree/test_files/template.ipynb)\n",
+    "[![launch on Colab](https://colab.research.google.com/assets/colab-badge.svg)](https://colab.research.google.com/github/open-atmos/devops_tests/blob/main/test_files/template.ipynb)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "bb84c52c57a729ec",
+   "metadata": {},
+   "source": "### based on Table 1 from B.Bolin 1958 \"On the use of tritium as a tracer for water in nature\" (https://digitallibrary.un.org/record/3892725)"
+  },
+  {
+   "cell_type": "code",
+   "id": "5e4e58050cc01f64",
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2024-12-05T21:30:20.214060Z",
+     "start_time": "2024-12-05T21:30:20.209298Z"
+    }
+   },
+   "source": [
+    "import sys\n",
+    "if 'google.colab' in sys.modules:\n",
+    "    !pip --quiet install open-atmos-jupyter-utils\n",
+    "    from open_atmos_jupyter_utils import pip_install_on_colab\n",
+    "    pip_install_on_colab('PySDM-examples')"
+   ],
+   "outputs": [],
+   "execution_count": 1
+  },
+  {
+   "cell_type": "code",
+   "id": "95f360dc62f373f9",
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2024-12-05T21:30:21.622672Z",
+     "start_time": "2024-12-05T21:30:20.218383Z"
+    }
+   },
+   "source": [
+    "import numpy as np\n",
+    "import pandas\n",
+    "from PySDM.physics import in_unit, si\n",
+    "from PySDM import Formulae"
+   ],
+   "outputs": [],
+   "execution_count": 2
+  },
+  {
+   "cell_type": "code",
+   "id": "e29a92d6680c4e51",
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2024-12-05T21:30:21.689116Z",
+     "start_time": "2024-12-05T21:30:21.687254Z"
+    }
+   },
+   "source": "radii = np.asarray([0.005, 0.01, 0.025, 0.05, 0.075, 0.1, 0.15, 0.20]) * si.cm",
+   "outputs": [],
+   "execution_count": 3
+  },
+  {
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2024-12-05T21:30:24.420354Z",
+     "start_time": "2024-12-05T21:30:21.698737Z"
+    }
+   },
+   "cell_type": "code",
+   "source": [
+    "formulae = Formulae(\n",
+    "    terminal_velocity='RogersYau',\n",
+    "    drop_growth='Mason1951',\n",
+    "    diffusion_thermics='Neglect',\n",
+    "    saturation_vapour_pressure='AugustRocheMagnus',\n",
+    "    ventilation='Froessling1938',\n",
+    "    particle_shape_and_density='LiquidSpheres',\n",
+    "    air_dynamic_viscosity='ZografosEtAl1987',\n",
+    "    constants={'BOLIN_ISOTOPE_TIMESCALE_COEFF_C1': 1.63},\n",
+    "    isotope_relaxation_timescale='Bolin1958',\n",
+    ")\n",
+    "const = formulae.constants\n",
+    "any_non_zero_value = 44.\n",
+    "temperature = const.T0 + 10 * si.K\n",
+    "pressure = const.p_STP\n",
+    "pvs = formulae.saturation_vapour_pressure.pvs_water(temperature)\n",
+    "v_term = formulae.terminal_velocity.v_term(radii)\n",
+    "eta_air=formulae.air_dynamic_viscosity.eta_air(temperature)\n",
+    "D=formulae.diffusion_thermics.D(T=temperature, p=pressure)\n",
+    "\n",
+    "air_density = pressure/const.Rd/temperature\n",
+    "assert abs(air_density - 1)/air_density <.3\n",
+    "Re = formulae.particle_shape_and_density.reynolds_number(\n",
+    "     radius=radii,\n",
+    "     velocity_wrt_air=v_term,\n",
+    "     dynamic_viscosity=eta_air,\n",
+    "     density=air_density,\n",
+    ")\n",
+    "Sc = formulae.trivia.air_schmidt_number(\n",
+    "    dynamic_viscosity=eta_air, \n",
+    "    diffusivity=D, \n",
+    "    density=air_density,\n",
+    ")\n",
+    "F = formulae.ventilation.ventilation_coefficient(sqrt_re_times_cbrt_sc=Re**(1/2) * Sc**(1/3))\n",
+    "\n",
+    "r_dr_dt = formulae.drop_growth.r_dr_dt(\n",
+    "    RH_eq=1,\n",
+    "    T=temperature,\n",
+    "    RH=0,\n",
+    "    lv=0,\n",
+    "    pvs=pvs,\n",
+    "    D=D,\n",
+    "    K=any_non_zero_value,\n",
+    "    ventilation_factor=F\n",
+    ")\n",
+    "adjustment_time = formulae.isotope_relaxation_timescale.tau(radius = radii, r_dr_dt = r_dr_dt)\n",
+    "\n",
+    "\n",
+    "pandas.options.display.float_format = '{:>,.2g}'.format\n",
+    "data = pandas.DataFrame({\n",
+    "    'radius [cm]': in_unit(radii, si.cm),\n",
+    "    'adjustment time [s]': adjustment_time,\n",
+    "    'terminal velocity [m/s]': v_term,\n",
+    "    'distance [m]': v_term * adjustment_time,\n",
+    "})\n",
+    "\n",
+    "data # pylint: disable=pointless-statement"
+   ],
+   "id": "314f42c310883cfa",
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "   radius [cm]  adjustment time [s]  terminal velocity [m/s]  distance [m]\n",
+       "0        0.005                  1.7                      0.4          0.69\n",
+       "1         0.01                  5.4                      0.8           4.3\n",
+       "2        0.025                   20                        2            40\n",
+       "3         0.05                   48                        4       1.9e+02\n",
+       "4        0.075                   81                      5.5       4.4e+02\n",
+       "5          0.1              1.2e+02                      6.4       7.6e+02\n",
+       "6         0.15                2e+02                      7.8       1.6e+03\n",
+       "7          0.2                3e+02                        9       2.7e+03"
+      ],
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>radius [cm]</th>\n",
+       "      <th>adjustment time [s]</th>\n",
+       "      <th>terminal velocity [m/s]</th>\n",
+       "      <th>distance [m]</th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.005</td>\n",
+       "      <td>1.7</td>\n",
+       "      <td>0.4</td>\n",
+       "      <td>0.69</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>1</th>\n",
+       "      <td>0.01</td>\n",
+       "      <td>5.4</td>\n",
+       "      <td>0.8</td>\n",
+       "      <td>4.3</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>2</th>\n",
+       "      <td>0.025</td>\n",
+       "      <td>20</td>\n",
+       "      <td>2</td>\n",
+       "      <td>40</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>3</th>\n",
+       "      <td>0.05</td>\n",
+       "      <td>48</td>\n",
+       "      <td>4</td>\n",
+       "      <td>1.9e+02</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>4</th>\n",
+       "      <td>0.075</td>\n",
+       "      <td>81</td>\n",
+       "      <td>5.5</td>\n",
+       "      <td>4.4e+02</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>5</th>\n",
+       "      <td>0.1</td>\n",
+       "      <td>1.2e+02</td>\n",
+       "      <td>6.4</td>\n",
+       "      <td>7.6e+02</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>6</th>\n",
+       "      <td>0.15</td>\n",
+       "      <td>2e+02</td>\n",
+       "      <td>7.8</td>\n",
+       "      <td>1.6e+03</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>7</th>\n",
+       "      <td>0.2</td>\n",
+       "      <td>3e+02</td>\n",
+       "      <td>9</td>\n",
+       "      <td>2.7e+03</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "</div>"
+      ]
+     },
+     "execution_count": 4,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "execution_count": 4
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.9.6"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}