diff --git a/examples/PyMPDATA_examples/Jaruga_et_al_2015/__init__.py b/examples/PyMPDATA_examples/Jaruga_et_al_2015/__init__.py
new file mode 100644
index 00000000..e69de29b
diff --git a/examples/PyMPDATA_examples/Jaruga_et_al_2015/fig19.ipynb b/examples/PyMPDATA_examples/Jaruga_et_al_2015/fig19.ipynb
new file mode 100644
index 00000000..f126a912
--- /dev/null
+++ b/examples/PyMPDATA_examples/Jaruga_et_al_2015/fig19.ipynb
@@ -0,0 +1,637 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "id": "8519709f-5065-4493-9af6-b5857dfc41b2",
+ "metadata": {},
+ "source": []
+ },
+ {
+ "cell_type": "markdown",
+ "id": "76ffcb39-d58b-4c7b-b795-a4206a63716d",
+ "metadata": {},
+ "source": []
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "id": "9ef2a5c3-3ac9-418b-92dd-e422ba66af66",
+ "metadata": {},
+ "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('PyMPDATA-examples')"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "id": "69583de6-8b30-4a74-8cd0-d3ade00ec579",
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "from open_atmos_jupyter_utils import show_plot"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 16,
+ "id": "28a63fae-2714-4fd3-af80-acb56dd0ce76",
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "image/svg+xml": [
+ "\n",
+ "\n",
+ "\n"
+ ],
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "application/vnd.jupyter.widget-view+json": {
+ "model_id": "495b5735ebc841fc815c45fb93eb955d",
+ "version_major": 2,
+ "version_minor": 0
+ },
+ "text/plain": [
+ "HBox(children=(HTML(value=\".\\\\fig_19.pdf
\"), HTML(value=\"