add isotopes to simulation kinematic 2d

Author: AgnieszkaZaba

examples/PySDM_examples/Strzabala_2026/__init__.py

@@ -0,0 +1,56 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "initial_id",
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": [
+    "from PySDM_examples.Arabas_et_al_2015 import Settings, SpinUp\n",
+    "from PySDM_examples.utils.kinematic_2d import Simulation"
+   ]
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
+   "source": [
+    "settings = Settings()\n",
+    "storage = Storage()\n",
+    "simulation = Simulation(settings, storage, SpinUp=SpinUp)\n",
+    "advectees_init_profiles = {\n",
+    "    \"HDO vapour mixing ratio\": np.full(shape=settings.grid[1], fill_value=1e-6),\n",
+    "    \"H18O vapour mixing ratio\": np.full(shape=settings.grid[1], fill_value=1e-6),\n",
+    "    \"H17O vapour mixing ratio\": np.full(shape=settings.grid[1], fill_value=1e-6),\n",
+    "}\n",
+    "simulation.reinit(additional_advectees_initial_profiles=advectees_init_profiles)"
+   ],
+   "id": "c340b54ecdfefc85"
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 2
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython2",
+   "version": "2.7.6"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

examples/PySDM_examples/Strzabala_2026/kinematic_2d_isotopes.ipynb

@@ -1,4 +1,6 @@
 import numpy as np
+from numpy import ndarray
+
 from PySDM_examples.utils.kinematic_2d.make_default_product_collection import (
     make_default_product_collection,
 )
@@ -32,7 +34,11 @@ def __init__(self, settings, storage, SpinUp, backend_class=CPU):
     def products(self):
         return self.particulator.products
 
-    def reinit(self, products=None):
+    def reinit(
+        self,
+        products=None,
+        additional_advectees_initial_profiles: Dict[str, ndarray] = None,
+    ):
         formulae = self.settings.formulae
         backend = self.backend_class(formulae=formulae)
         environment = Kinematic2D(
@@ -70,6 +76,7 @@ def reinit(self, products=None):
             initial_profiles = {
                 "th": self.settings.initial_dry_potential_temperature_profile,
                 "water_vapour_mixing_ratio": self.settings.initial_vapour_mixing_ratio_profile,
+                **(additional_advectees_initial_profiles or {}),
             }
             advectees = dict(
                 (