PySDM 2.120__tar.gz → 2.121__tar.gz

{pysdm-2.120 → pysdm-2.121}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pysdm
-Version: 2.120
+Version: 2.121
 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
 Author-email: "https://github.com/open-atmos/PySDM/graphs/contributors" <sylwester.arabas@agh.edu.pl>
 License-Expression: GPL-3.0-only

{pysdm-2.120 → pysdm-2.121}/PySDM/backends/impl_numba/methods/condensation_methods.py

@@ -388,7 +388,8 @@ class CondensationMethods(BackendMethods):
             pvs,
             D,
             K,
-            ventilation_factor,
+            mass_ventilation_factor,
+            heat_ventilation_factor,
         ):
             if x_new > formulae.diffusion_coordinate__x_max():
                 return x_old - x_new
@@ -410,9 +411,8 @@ class CondensationMethods(BackendMethods):
                 RH,
                 lv,
                 pvs,
-                D,
-                K,
-                ventilation_factor,
+                D * mass_ventilation_factor,
+                K * heat_ventilation_factor,
             )
             dm_dt = formulae.particle_shape_and_density__dm_dt(r=r_new, r_dr_dt=r_dr_dt)
             return (
@@ -469,12 +469,13 @@ class CondensationMethods(BackendMethods):
                 ):
                     Dr = formulae.diffusion_kinetics__D(DTp, r_old, lambdaD)
                     Kr = formulae.diffusion_kinetics__K(KTp, r_old, lambdaK)
-                    ventilation_factor = formulae.ventilation__ventilation_coefficient(
+                    mass_ventilation_factor = formulae.ventilation__ventilation_coefficient(
                         sqrt_re_times_cbrt_sc=formulae.trivia__sqrt_re_times_cbrt_sc(
                             Re=attributes.reynolds_number[drop],
                             Sc=Sc,
                         )
                     )
+                    heat_ventilation_factor = mass_ventilation_factor  # TODO #1588
                     args = (
                         x_old,
                         timestep,
@@ -487,7 +488,8 @@ class CondensationMethods(BackendMethods):
                         pvs,
                         Dr,
                         Kr,
-                        ventilation_factor,
+                        mass_ventilation_factor,
+                        heat_ventilation_factor,
                     )
                     r_dr_dt_old = formulae.drop_growth__r_dr_dt(
                         RH_eq,
@@ -495,9 +497,8 @@ class CondensationMethods(BackendMethods):
                         RH,
                         lv,
                         pvs,
-                        Dr,
-                        Kr,
-                        ventilation_factor,
+                        mass_ventilation_factor * Dr,
+                        heat_ventilation_factor * Kr,
                     )
                     mass_old = formulae.diffusion_coordinate__mass(x_old)
                     dm_dt_old = formulae.particle_shape_and_density__dm_dt(

{pysdm-2.120 → pysdm-2.121}/PySDM/backends/impl_numba/methods/deposition_methods.py

@@ -57,12 +57,13 @@ class DepositionMethods(BackendMethods):  # pylint:disable=too-few-public-method
 
                     capacity = formulae.diffusion_ice_capacity__capacity(diameter)
 
-                    ventilation_factor = formulae.ventilation__ventilation_coefficient(
+                    mass_ventilation_factor = formulae.ventilation__ventilation_coefficient(
                         sqrt_re_times_cbrt_sc=formulae.trivia__sqrt_re_times_cbrt_sc(
                             Re=reynolds_number[i],
                             Sc=schmidt_number[cid],
                         )
                     )
+                    heat_ventilation_factor = mass_ventilation_factor  # TODO #1588
 
                     Dv_const = formulae.diffusion_thermics__D(temperature, pressure)
                     lambdaD = formulae.diffusion_ice_kinetics__lambdaD(
@@ -93,9 +94,8 @@ class DepositionMethods(BackendMethods):  # pylint:disable=too-few-public-method
                             RH=saturation_ratio_ice,
                             lv=latent_heat_sub,
                             pvs=pvs_ice,
-                            D=diffusion_coefficient,
-                            K=thermal_conductivity,
-                            ventilation_factor=ventilation_factor,
+                            D=mass_ventilation_factor * diffusion_coefficient,
+                            K=heat_ventilation_factor * thermal_conductivity,
                         )
                         * formulae.constants.rho_w
                     )

{pysdm-2.120 → pysdm-2.121}/PySDM/backends/impl_numba/test_helpers/scipy_ode_condensation_solver.py

@@ -132,12 +132,13 @@ def _make_solve(formulae):  # pylint: disable=too-many-statements,too-many-local
             r = jit_formulae.trivia__radius(v)
             Dr = jit_formulae.diffusion_kinetics__D(DTp, r, lambdaD)
             Kr = jit_formulae.diffusion_kinetics__K(KTp, r, lambdaK)
-            ventilation_factor = jit_formulae.ventilation__ventilation_coefficient(
+            mass_ventilation_factor = jit_formulae.ventilation__ventilation_coefficient(
                 sqrt_re_times_cbrt_sc=jit_formulae.trivia__sqrt_re_times_cbrt_sc(
                     Re=reynolds_number[i],
                     Sc=schmidt_number,
                 )
             )
+            heat_ventilation_factor = mass_ventilation_factor  # TODO #1588
             sgm = jit_formulae.surface_tension__sigma(T, v, dry_volume[i], f_org[i])
             r_dr_dt = jit_formulae.drop_growth__r_dr_dt(
                 jit_formulae.hygroscopicity__RH_eq(
@@ -147,9 +148,8 @@ def _make_solve(formulae):  # pylint: disable=too-many-statements,too-many-local
                 RH,
                 lv,
                 pvs,
-                Dr,
-                Kr,
-                ventilation_factor,
+                mass_ventilation_factor * Dr,
+                heat_ventilation_factor * Kr,
             )
             dm_dt = jit_formulae.particle_shape_and_density__dm_dt(r, r_dr_dt)
             dy_dt[idx_x + i] = jit_formulae.diffusion_coordinate__dx_dt(m, dm_dt)

{pysdm-2.120 → pysdm-2.121}/PySDM/backends/impl_thrust_rtc/methods/condensation_methods.py

@@ -25,7 +25,7 @@ ARGS_VARS = (
     "_pvs",
     "Dr",
     "Kr",
-    "ventilation_factor",
+    "mass_ventilation_factor",
 )
 
 
@@ -125,10 +125,9 @@ class CondensationMethods(
                         RH=args("_RH"),
                         lv=args("_lv"),
                         pvs=args("_pvs"),
-                        D=args("Dr"),
-                        K=args("Kr"),
-                        ventilation_factor=args("ventilation_factor"),
-                    )};
+                        D=f'{args("mass_ventilation_factor")}*{args("Dr")}',
+                        K=f'{args("mass_ventilation_factor")}*{args("Kr")}',
+                    )}; // TODO #1588
                     auto dm_dt = {phys.particle_shape_and_density.dm_dt.c_inline(
                         r="r_new", r_dr_dt="r_dr_dt"
                     )};
@@ -167,7 +166,7 @@ class CondensationMethods(
             real_type Dr=0;
             real_type Kr=0;
             real_type qrt_re_times_cbrt_sc=0;
-            real_type ventilation_factor=0;
+            real_type mass_ventilation_factor=0;
             real_type r_dr_dt_old=0;
             real_type dm_dt_old=0;
             real_type dx_old=0;
@@ -186,12 +185,13 @@ class CondensationMethods(
                     Re="reynolds_number[i]",
                     Sc="_schmidt_number",
                 )};
-                ventilation_factor = {phys.ventilation.ventilation_coefficient.c_inline(
+                mass_ventilation_factor = {phys.ventilation.ventilation_coefficient.c_inline(
                     sqrt_re_times_cbrt_sc="qrt_re_times_cbrt_sc"
                 )};
+                auto heat_ventilation_factor = mass_ventilation_factor; // TODO #1588
                 r_dr_dt_old = {phys.drop_growth.r_dr_dt.c_inline(
-                    RH_eq="RH_eq", T="_T", RH="_RH", lv="_lv", pvs="_pvs", D="Dr", K="Kr",
-                    ventilation_factor="ventilation_factor",
+                    RH_eq="RH_eq", T="_T", RH="_RH", lv="_lv", pvs="_pvs",
+                    D="mass_ventilation_factor*Dr", K="heat_ventilation_factor*Kr",
                 )};
                 dm_dt_old = {phys.particle_shape_and_density.dm_dt.c_inline(r="r_old", r_dr_dt="r_dr_dt_old")};
                 dx_old = dt * {phys.diffusion_coordinate.dx_dt.c_inline(

{pysdm-2.120 → pysdm-2.121}/PySDM/physics/constants_defaults.py

@@ -681,6 +681,12 @@ def compute_derived_values(c: dict):
     - [IAPWS Guidelines](http://www.iapws.org/relguide/fundam.pdf)
     """
 
+    c["M_1H2_16O"] = c["M_1H"] * 2 + c["M_16O"]
+    c["M_2H_1H_16O"] = c["M_2H"] + c["M_1H"] + c["M_16O"]
+    c["M_3H_1H_16O"] = c["M_3H"] + c["M_1H"] + c["M_16O"]
+    c["M_1H2_17O"] = c["M_1H"] * 2 + c["M_17O"]
+    c["M_1H2_18O"] = c["M_1H"] * 2 + c["M_18O"]
+
     c["Mv"] = (
         (
             1
@@ -689,19 +695,15 @@ def compute_derived_values(c: dict):
             - 1 * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_17O"])
             - 1 * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_18O"])
         )
-        * (c["M_1H"] * 2 + c["M_16O"])
+        * c["M_1H2_16O"]
         + 2
         * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_2H"])
-        * (c["M_2H"] + c["M_1H"] + c["M_16O"])
+        * c["M_2H_1H_16O"]
         + 2
         * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_3H"])
-        * (c["M_3H"] + c["M_1H"] + c["M_16O"])
-        + 1
-        * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_17O"])
-        * (c["M_1H"] * 2 + c["M_17O"])
-        + 1
-        * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_18O"])
-        * (c["M_1H"] * 2 + c["M_18O"])
+        * c["M_3H_1H_16O"]
+        + 1 * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_17O"]) * c["M_1H2_17O"]
+        + 1 * Trivia.mixing_ratio_to_specific_content(c["VSMOW_R_18O"]) * c["M_1H2_18O"]
     )
 
     c["eps"] = c["Mv"] / c["Md"]

{pysdm-2.120 → pysdm-2.121}/PySDM/physics/drop_growth/__init__.py

@@ -2,5 +2,6 @@
 Formulation of the coupled heat-moisture diffusion problem
 """
 
-from .mason_1951 import Mason1951
+from .fick import Fick
+from .howell_1949 import Howell1949
 from .mason_1971 import Mason1971

pysdm-2.121/PySDM/physics/drop_growth/fick.py

@@ -0,0 +1,13 @@
+"""
+Fickian diffusion only drop growth
+"""
+
+
+class Fick:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    # pylint: disable=too-many-arguments
+    @staticmethod
+    def r_dr_dt(const, RH_eq, T, RH, lv, pvs, D, K):  # pylint: disable=unused-argument
+        return (RH - RH_eq) / const.rho_w / (const.Rv * T / D / pvs)

pysdm-2.121/PySDM/physics/drop_growth/howell_1949.py

@@ -0,0 +1,19 @@
+"""
+single-equation approximation of the vapour and heat diffusion problem
+as proposed in [Howell 1949](https://doi.org/10.1175/1520-0469(1949)006%3C0134:TGOCDI%3E2.0.CO;2)
+same as in [Mason 1951](https://doi.org/10.1088/0370-1301/64/9/307)
+"""
+
+
+class Howell1949:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    # pylint: disable=too-many-arguments
+    @staticmethod
+    def r_dr_dt(const, RH_eq, T, RH, lv, pvs, D, K):
+        return (
+            (RH - RH_eq)
+            / const.rho_w
+            / (const.Rv * T / D / pvs + lv**2 / K / T**2 / const.Rv)
+        )

{pysdm-2.120 → pysdm-2.121}/PySDM/physics/drop_growth/mason_1971.py

@@ -11,10 +11,9 @@ class Mason1971:  # pylint: disable=too-few-public-methods
 
     # pylint: disable=too-many-arguments
     @staticmethod
-    def r_dr_dt(const, RH_eq, T, RH, lv, pvs, D, K, ventilation_factor):
+    def r_dr_dt(const, RH_eq, T, RH, lv, pvs, D, K):
         return (
-            ventilation_factor
-            * (RH - RH_eq)
+            (RH - RH_eq)
             / const.rho_w
             / (const.Rv * T / D / pvs + lv / K / T * (lv / T / const.Rv - 1))
         )

{pysdm-2.120 → pysdm-2.121}/PySDM/physics/isotope_kinetic_fractionation_factors/jouzel_and_merlivat_1984.py

@@ -10,12 +10,12 @@ class JouzelAndMerlivat1984:  # pylint: disable=too-few-public-methods
 
     @staticmethod
     def alpha_kinetic(
-        alpha_equilibrium, relative_humidity, heavy_to_light_diffusivity_ratio
+        alpha_equilibrium, saturation_over_ice, heavy_to_light_diffusivity_ratio
     ):
-        """eq. (11) or eq. (14)"""
-        return relative_humidity / (
+        """eq. (11)"""
+        return saturation_over_ice / (
             alpha_equilibrium
             / heavy_to_light_diffusivity_ratio
-            * (relative_humidity - 1)
+            * (saturation_over_ice - 1)
             + 1
         )

{pysdm-2.120 → pysdm-2.121}/docs/bibliography.json

@@ -558,7 +558,7 @@
   },
   "https://doi.org/10.1088/0370-1301/64/9/307": {
     "usages": [
-      "PySDM/physics/drop_growth/mason_1951.py"
+      "PySDM/physics/drop_growth/howell_1949.py"
     ],
     "label": "Mason 1951 (Proc. Phys. Soc. B 64)",
     "title": "Spontaneous Condensation of Water Vapour in Expansion Chamber Experiments"
@@ -836,7 +836,8 @@
   },
   "https://doi.org/10.1175/1520-0469(1949)006%3C0134:TGOCDI%3E2.0.CO;2": {
     "usages": [
-      "PySDM/backends/impl_numba/methods/deposition_methods.py"
+      "PySDM/backends/impl_numba/methods/deposition_methods.py",
+      "PySDM/physics/drop_growth/howell_1949.py"
     ],
     "title": "The Growth of Cloud Drops in Uniformly Cooled Air",
     "label": "Howell 1949 (J. Atmos. Sci.)"
@@ -844,7 +845,11 @@
   "https://doi.org/10.1029/JD089iD07p11749": {
     "usages": [
       "PySDM/physics/isotope_kinetic_fractionation_factors/jouzel_and_merlivat_1984.py",
-      "tests/unit_tests/physics/test_isotope_kinetic_fractionation_factors.py"
+      "examples/PySDM_examples/Jouzel_and_Merlivat_1984/__init__.py",
+      "examples/PySDM_examples/Jouzel_and_Merlivat_1984/fig_8_9.ipynb",
+      "examples/PySDM_examples/Jouzel_and_Merlivat_1984/thermodynamic_profiles.py",
+      "tests/unit_tests/physics/test_isotope_kinetic_fractionation_factors.py",
+      "examples/PySDM_examples/Fisher_1991/__init__.py"
     ],
     "title": "Deuterium and oxygen 18 in precipitation: Modeling of the isotopic effects during snow formation",
     "label": "Jouzel & Merlivat 1984 (J. Geophys. Res. Atmos. 89)"
@@ -862,5 +867,13 @@
     ],
     "title": "IAEA Reference Sheet for International Measurement Standards: VSMOW & SLAP",
     "label": "IAEA 2006"
+  },
+  "https://doi.org/10.3402/tellusb.v43i5.15414": {
+    "usages": [
+      "examples/PySDM_examples/Fisher_1991/fig_2.ipynb",
+      "examples/PySDM_examples/Fisher_1991/__init__.py"
+    ],
+    "title": "Remarks on the deuterium excess in precipitation in cold regions",
+    "label": "Fisher 1991 (Tellus B)"
   }
 }

{pysdm-2.120 → pysdm-2.121}/examples/PySDM_examples/Bolin_1958/table_1.ipynb

@@ -20,70 +20,140 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
    "id": "5e4e58050cc01f64",
    "metadata": {
     "ExecuteTime": {
-     "end_time": "2024-12-05T21:30:20.214060Z",
-     "start_time": "2024-12-05T21:30:20.209298Z"
+     "end_time": "2025-05-15T21:53:10.897690Z",
+     "start_time": "2025-05-15T21:53:10.890113Z"
     }
    },
-   "outputs": [],
    "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": 9
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
    "id": "95f360dc62f373f9",
    "metadata": {
     "ExecuteTime": {
-     "end_time": "2024-12-05T21:30:21.622672Z",
-     "start_time": "2024-12-05T21:30:20.218383Z"
+     "end_time": "2025-05-15T21:53:10.907788Z",
+     "start_time": "2025-05-15T21:53:10.904663Z"
     }
    },
-   "outputs": [],
    "source": [
     "import numpy as np\n",
     "import pandas\n",
     "from PySDM.physics import in_unit, si\n",
     "from PySDM import Formulae\n",
     "from IPython.display import display"
-   ]
+   ],
+   "outputs": [],
+   "execution_count": 10
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
    "id": "e29a92d6680c4e51",
    "metadata": {
     "ExecuteTime": {
-     "end_time": "2024-12-05T21:30:21.689116Z",
-     "start_time": "2024-12-05T21:30:21.687254Z"
+     "end_time": "2025-05-15T21:53:10.922280Z",
+     "start_time": "2025-05-15T21:53:10.920362Z"
     }
    },
-   "outputs": [],
    "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": 11
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
    "id": "314f42c310883cfa",
    "metadata": {
     "ExecuteTime": {
-     "end_time": "2024-12-05T21:30:24.420354Z",
-     "start_time": "2024-12-05T21:30:21.698737Z"
+     "end_time": "2025-05-15T21:53:11.229114Z",
+     "start_time": "2025-05-15T21:53:10.930920Z"
     }
    },
+   "source": [
+    "formulae = Formulae(\n",
+    "    terminal_velocity='RogersYau',\n",
+    "    drop_growth='Howell1949',\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",
+    "mass_ventilation = formulae.ventilation.ventilation_coefficient(\n",
+    "    sqrt_re_times_cbrt_sc=Re**(1/2) * Sc**(1/3)\n",
+    ")\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  *mass_ventilation,\n",
+    "    K=any_non_zero_value\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",
+    "display(data)"
+   ],
    "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",
@@ -169,90 +239,26 @@
        "  </tbody>\n",
        "</table>\n",
        "</div>"
-      ],
-      "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"
       ]
      },
      "metadata": {},
      "output_type": "display_data"
     }
    ],
-   "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",
-    "display(data)"
-   ]
+   "execution_count": 12
   },
   {
    "cell_type": "code",
-   "execution_count": null,
    "id": "0bef57bc-7e19-4a3f-8f77-6ec762e8f236",
-   "metadata": {},
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2025-05-15T21:53:11.238952Z",
+     "start_time": "2025-05-15T21:53:11.237500Z"
+    }
+   },
+   "source": [],
    "outputs": [],
-   "source": []
+   "execution_count": null
   }
  ],
  "metadata": {

pysdm-2.121/examples/PySDM_examples/Fisher_1991/__init__.py

@@ -0,0 +1,11 @@
+# pylint: disable=invalid-name
+"""
+based on the [Jouzel and Merlivat 1984](https://doi.org/10.1029/JD089iD07p11749) calculations
+of delta_{18O} and delta_D used in d-excess formula and solved
+by integration of the eq. (1) in
+[Fisher 1991](https://doi.org/10.3402/tellusb.v43i5.15414) and eq. (3) in Jouzel and Merlivat 1984
+
+
+fig_2.ipynb:
+.. include:: ./fig_2.ipynb.badges.md
+"""