timflow 0.3.0__tar.gz → 0.4.0__tar.gz

{timflow-0.3.0 → timflow-0.4.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: timflow
-Version: 0.3.0
+Version: 0.4.0
 Summary: Analytic element modeling of groundwater flow
 Author: Mark Bakker, Davíd Brakenhoff
 Maintainer-email: Mark Bakker <markbak@gmail.com>, Davíd Brakenhoff <d.brakenhoff@artesia-water.nl>

{timflow-0.3.0 → timflow-0.4.0}/docs/about/index.rst

@@ -11,8 +11,8 @@ Steady flow
 .. bibliography:: publications.bib
     :list: bullet
 
-    bakker_semi-analytic_2013
-    bakker_analytic_2013
+    bakker2003
+    bakker_analytic_2006
 
 Transient flow
 --------------

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/00userguide/howtos/howto_convert_timml_model.ipynb

@@ -114,7 +114,7 @@
    "id": "c00eb393-d11b-450a-acab-32192e688745",
    "metadata": {},
    "source": [
-    "The same procedure holds for combining other plots. For example combining two contour plots, or a contour plot plus a capture zone. Note that the `plotcapzone` function has moved to the `plots` module. All three can be combined as follows"
+    "The same procedure holds for combining other plots. For example combining two contour plots, or a contour plot plus a capture zone. Note that the `plot_capture_zone` function has moved to the `plots` module. All three can be combined as follows"
    ]
   },
   {
@@ -125,7 +125,7 @@
    "outputs": [],
    "source": [
     "ax0 = ml.plots.contour(win=[-1000, 100, -500, 500], ngr=50, levels=10, labels=False)\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w, nt=10, zstart=[1], hstepmax=20, tmax=10 * 365, color=\"C1\", ax=ax0\n",
     ")\n",
     "ml.plots.contour(win=[-1000, 100, -500, 500], ngr=50, levels=[40], color=\"k\", ax=ax0);"

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/00userguide/howtos/howto_select_a_model.ipynb

@@ -75,7 +75,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "ax = ml.plots.xsection(params=True, labels=True, units={\"k\": \"m/d\", \"c\": \"d\"})\n",
+    "ax = ml.plots.xsection(params=True, labels=True, units={\"kaq\": \"m/d\", \"c\": \"d\"})\n",
     "ax.set_frame_on(False)"
    ]
   },
@@ -134,7 +134,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "ax = ml.plots.xsection(params=True, labels=True, units={\"k\": \"m/d\", \"c\": \"d\"})\n",
+    "ax = ml.plots.xsection(params=True, labels=True, units={\"kaq\": \"m/d\"})\n",
     "ax.set_frame_on(False)"
    ]
   },
@@ -179,7 +179,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "ax = ml.plots.xsection(params=True, labels=True, units={\"k\": \"m/d\", \"c\": \"d\"})\n",
+    "ax = ml.plots.xsection(params=True, labels=True, units={\"kaq\": \"m/d\", \"c\": \"d\"})\n",
     "ax.set_frame_on(False)"
    ]
   },

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/00userguide/tutorials/tutorial0_well_single_layer_aquifer.ipynb

@@ -275,7 +275,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "Alternatively, a capture zone may be computed with the `plotcapzone` function by starting pathlines at the well radius and tracing the pathlines against the flow. The number of pathlines must be specified and the maximum time, which is specified as 20 years here, but since the time units are days, we specify the number of days. "
+    "Alternatively, a capture zone may be computed with the `plot_capture_zone` function by starting pathlines at the well radius and tracing the pathlines against the flow. The number of pathlines must be specified and the maximum time, which is specified as 20 years here, but since the time units are days, we specify the number of days. "
    ]
   },
   {
@@ -291,7 +291,7 @@
     "    labels=False,\n",
     "    figsize=(5, 5),\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w,\n",
     "    nt=20,  # number of pathlines\n",
     "    tmax=20 * 365.25,  # maximum time, d\n",

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/02examples/00_single_layer_model.ipynb

@@ -106,7 +106,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "After the well is added (or any other elements), the model needs to be solved again. A contour plot is created and a 10 strace line are added. The stepsize is given in meters and represents the largest space step that is taken, but it is reduced when certain accuracy constraints are not met. Note that, after running the code cell below, for each trace line it is printed to the screen what the reason was that the traceline was aborted. In this case it was either because the trace line reached a well or reached the maximum number of steps (the default is 100 steps, but this can be changed by specifying the `nstepmax` keyword).  "
+    "After the well is added (or any other elements), the model needs to be solved again. A contour plot is created and a trace line is added. The stepsize is given in meters and represents the largest space step that is taken, but it is reduced when certain accuracy constraints are not met. Note that, after running the code cell below, a message is printed to the screen with the reason why the traceline was aborted. In this case it was because the trace line reached the well. The message an be suppressed by with the keyword `silent=True` or the default `silent='.'`. "
    ]
   },
   {
@@ -120,8 +120,38 @@
     "    win=[-1000, 100, -500, 500], ngr=50, levels=np.arange(39, 42, 0.1), figsize=(6, 6)\n",
     ")\n",
     "ml.plots.tracelines(\n",
-    "    -800 * np.ones(1), -200 * np.ones(1), np.zeros(1), hstepmax=20, color=\"C1\", ax=ax\n",
-    ")"
+    "    -800 * np.ones(1),\n",
+    "    -200 * np.ones(1),\n",
+    "    np.zeros(1),\n",
+    "    hstepmax=20,\n",
+    "    color=\"C1\",\n",
+    "    ax=ax,\n",
+    "    silent=False,\n",
+    ");"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "If you don't want to have contour in the background, simply specify a window for the `tracelines` function and no axes. "
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "ml.plots.tracelines(\n",
+    "    -800 * np.ones(1),\n",
+    "    -200 * np.ones(1),\n",
+    "    np.zeros(1),\n",
+    "    hstepmax=20,\n",
+    "    color=\"C1\",\n",
+    "    silent=False,\n",
+    "    win=[-1000, 100, -500, 500],\n",
+    ");"
    ]
   },
   {
@@ -276,7 +306,7 @@
     "    layers=0,\n",
     "    figsize=(6, 6),\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w, hstepmax=20, nt=20, zstart=0, tmax=5 * 365.25, color=\"C1\", ax=ax\n",
     ")"
    ]
@@ -302,7 +332,7 @@
     "    layers=0,\n",
     "    figsize=(6, 6),\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w, hstepmax=20, nt=20, zstart=0, tmax=20 * 365.25, color=\"C1\", ax=ax\n",
     ")"
    ]

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/02examples/02_wells_rivers_and_recharge.ipynb

@@ -163,7 +163,7 @@
    "outputs": [],
    "source": [
     "axes = ml.plots.topview_and_xsection(win=[0, 20000, 0, 20000], figsize=figsize)\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w1,\n",
     "    hstepmax=50,\n",
     "    nt=10,\n",
@@ -172,13 +172,13 @@
     "    orientation=\"both\",\n",
     "    ax=axes,\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w2, hstepmax=50, nt=10, zstart=30, tmax=250 * 365.25, orientation=\"both\", ax=axes\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w3, hstepmax=50, nt=10, zstart=30, tmax=250 * 365.25, orientation=\"both\", ax=axes\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w3,\n",
     "    hstepmax=50,\n",
     "    nt=10,\n",
@@ -198,16 +198,16 @@
     "axes = ml.plots.topview_and_xsection(\n",
     "    win=[0, 20000, 0, 20000], topfigfrac=0.7, figsize=figsize\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w1, hstepmax=50, nt=10, zstart=150, tmax=50 * 365.25, orientation=\"both\", ax=axes\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w2, hstepmax=50, nt=10, zstart=30, tmax=50 * 365.25, orientation=\"both\", ax=axes\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w3, hstepmax=50, nt=10, zstart=30, tmax=50 * 365.25, orientation=\"both\", ax=axes\n",
     ")\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w3, hstepmax=50, nt=10, zstart=100, tmax=50 * 365.25, orientation=\"both\", ax=axes\n",
     ")"
    ]

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/02examples/03a_modeling_inhomogeneities.ipynb

@@ -202,7 +202,7 @@
    "outputs": [],
    "source": [
     "axes = ml.plots.topview_and_xsection(win=[-200, 800, 0, 800], figsize=figsize)\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w, hstepmax=50, nt=20, zstart=10, tmax=20 * 365.25, orientation=\"both\", ax=axes\n",
     ")"
    ]

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/02examples/03b_modeling_inhomogeneities3D.ipynb

@@ -174,7 +174,7 @@
    "outputs": [],
    "source": [
     "axes = ml.plots.topview_and_xsection(win=[-200, 800, 0, 800])\n",
-    "ml.plots.plotcapzone(\n",
+    "ml.plots.plot_capture_zone(\n",
     "    well=w, hstepmax=50, nt=20, zstart=15, tmax=20 * 365.25, orientation=\"both\", ax=axes\n",
     ")"
    ]

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/02examples/connected_wells.ipynb

@@ -572,7 +572,7 @@
    "source": [
     "ilay = 0\n",
     "levels = [-2.4, -2.2, -2.0, -1.8, -1.6, -1.4, -1.2]\n",
-    "ml2.plots.contour(\n",
+    "ax = ml2.plots.contour(\n",
     "    win=(-50, 50, -50, 50),\n",
     "    ngr=51,\n",
     "    levels=levels,\n",
@@ -581,7 +581,7 @@
     "    linestyles=\"solid\",\n",
     "    color=\"C1\",\n",
     ")\n",
-    "plt.plot(ws.xcp, ws.ycp, \"kx\");"
+    "ax.plot(ws.xcp, ws.ycp, \"kx\");"
    ]
   },
   {
@@ -720,17 +720,16 @@
    "source": [
     "ilay = 1\n",
     "levels = [-2.4, -2.2, -2.0, -1.8, -1.6, -1.4, -1.2]\n",
-    "ml3.plots.contour(\n",
+    "ax = ml3.plots.contour(\n",
     "    win=(-50, 50, -50, 50),\n",
     "    ngr=51,\n",
     "    levels=levels,\n",
     "    decimals=2,\n",
     "    layers=[ilay],\n",
-    "    newfig=False,\n",
     "    linestyles=\"solid\",\n",
     "    color=\"C1\",\n",
     ")\n",
-    "plt.plot(ws.xcp, ws.ycp, \"kx\");"
+    "ax.plot(ws.xcp, ws.ycp, \"kx\");"
    ]
   },
   {

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/03xsections/cross_section_models.ipynb

@@ -536,21 +536,9 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "timflow (3.13.11)",
+   "display_name": "Python 3",
    "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.13.11"
   }
  },
  "nbformat": 4,

{timflow-0.3.0 → timflow-0.4.0}/docs/steady/04benchmarks/test_linesink_discharge.ipynb

@@ -73,7 +73,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "ml1.plots.contour(\n",
+    "ax = ml1.plots.contour(\n",
     "    [-20, 20, -20, 20], 50, [0], np.arange(20, 31, 1), color=\"C0\", labels=False\n",
     ")\n",
     "ml2.plots.contour(\n",
@@ -83,7 +83,8 @@
     "    np.arange(20, 31, 1),\n",
     "    color=\"C1\",\n",
     "    labels=False,\n",
-    "    newfig=False,\n",
+    "    ax=ax,\n",
+    "    linestyles=\"dashed\",\n",
     ")"
    ]
   },

{timflow-0.3.0 → timflow-0.4.0}/docs/transient/00userguide/howtos/howto_select_a_model.ipynb

@@ -45,7 +45,10 @@
     "    tmax=10,\n",
     ")\n",
     "\n",
-    "ax = ml.plots.xsection(params=True)\n",
+    "ax = ml.plots.xsection(\n",
+    "    params=True,\n",
+    "    units={\"kaq\": \"m/d\", \"c\": \"d\", \"Saq\": \"m$^{-1}$\", \"Sll\": \"m$^{-1}$\"},\n",
+    ")\n",
     "ax.set_frame_on(False)"
    ]
   },
@@ -77,7 +80,7 @@
     "    tmax=10,\n",
     ")\n",
     "\n",
-    "ax = ml2.plots.xsection(params=True)\n",
+    "ax = ml2.plots.xsection(params=True, units={\"kaq\": \"m/d\", \"Saq\": \"m$^{-1}$\"})\n",
     "ax.set_frame_on(False)"
    ]
   }

{timflow-0.3.0 → timflow-0.4.0}/docs/transient/03xsections/river_in_cross_section.ipynb

@@ -301,7 +301,7 @@
     "    hstar=10,\n",
     "    sep=\"\\n\",\n",
     "    ax=ax,\n",
-    "    units={\"k\": \"m/d\", \"Ss\": \"m$^{-1}$\", \"c\": \"d\"},\n",
+    "    units={\"kaq\": \"m/d\", \"Saq\": \"m$^{-1}$\", \"c\": \"d\"},\n",
     "    fontsize=8,\n",
     ")\n",
     "ax.set_xlabel(\"x-coordinate [m]\")\n",

timflow-0.4.0/docs/transient/04pumpingtests/confined1_oude_korendijk.ipynb

@@ -0,0 +1,241 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {
+    "editable": true,
+    "slideshow": {
+     "slide_type": ""
+    },
+    "tags": []
+   },
+   "source": [
+    "# 1. Confined Aquifer Test - Oude Korendijk"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Import packages"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "import pandas as pd\n",
+    "\n",
+    "import timflow.transient as tft\n",
+    "\n",
+    "plt.rcParams[\"figure.figsize\"] = (5, 3)  # default figure size"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Introduction and Conceptual Model\n",
+    "\n",
+    "In this example, we will use the pumping test data from Oude Korendijk (Kruseman et al. 1970). \n",
+    "\n",
+    "Oude Korendijk is a polder area South of Rotterdam, the Netherlands. The stratigraphy can be summarized by:\n",
+    "* the first 18 m consist of almost impermeable material,\n",
+    "* the next 7 m is succession of coarse gravel and sands, which are considered as the aquifer layer,\n",
+    "* a layer of fine sands and clayey sediments that are deemed impermeable.\n",
+    "\n",
+    "The well screens the whole thickness of the confined aquifer. Drawdowns were measured in two piezometers located at distances of 30 m and 90 m from the well. The pumping well discharge was constant at 788 m$^3$/d for almost 14 hours. The objective of the pumping test was to estimate the hydraulic conductivity and the specific storage of the aquifer layer. "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "<img src=\"./figs/oude_korendijk.png\" style=\"width:400pt\">"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Load data"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# time and drawdown of piezometer 30m away from pumping well\n",
+    "data1 = np.loadtxt(\"data/piezometer_h30.txt\", skiprows=1)\n",
+    "to1 = data1[:, 0] / 60 / 24  # convert min to days\n",
+    "ho1 = data1[:, 1]\n",
+    "ro1 = 30\n",
+    "# time and drawdown of piezometer 90m away from pumping well\n",
+    "data2 = np.loadtxt(\"data/piezometer_h90.txt\", skiprows=1)\n",
+    "to2 = data2[:, 0] / 60 / 24  # convert min to days\n",
+    "ho2 = data2[:, 1]\n",
+    "ro2 = 90"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Parameters and model"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# known parameters\n",
+    "H = 7  # aquifer thickness in meters, m\n",
+    "zt = -18  # top boundary of aquifer, m\n",
+    "zb = zt - H  # bottom boundary of aquifer, m\n",
+    "Q = 788  # constant discharge, m3/d"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# timflow model\n",
+    "ml = tft.ModelMaq(kaq=60, z=[zt, zb], Saq=1e-4, tmin=1e-5, tmax=1)\n",
+    "w = tft.Well(model=ml, xw=0, yw=0, rw=0.2, tsandQ=[(0, Q)], layers=0)\n",
+    "ml.solve()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Estimate aquifer parameters\n",
+    "The aquifer parameters are estimated using head observations at both piezometers."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# unknown parameters: kaq, Saq\n",
+    "cal = tft.Calibrate(ml)\n",
+    "cal.set_parameter(name=\"kaq\", initial=10, layers=0)\n",
+    "cal.set_parameter(name=\"Saq\", initial=1e-4, layers=0)\n",
+    "cal.series(name=\"obs1\", x=ro1, y=0, t=to1, h=ho1, layer=0)  # Adding well 1\n",
+    "cal.series(name=\"obs2\", x=ro2, y=0, t=to2, h=ho2, layer=0)  # Adding well 2\n",
+    "cal.fit(report=False)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "display(cal.parameters.loc[:, [\"optimal\"]])\n",
+    "print(f\"RMSE: {cal.rmse():.3f} m\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "hm1 = ml.head(ro1, 0, to1)\n",
+    "hm2 = ml.head(ro2, 0, to2)\n",
+    "plt.semilogx(to1, ho1, \"C0.\", label=\"obs at 30m\")\n",
+    "plt.semilogx(to1, hm1[0], \"C0\", label=\"timflow at 30 m\")\n",
+    "plt.semilogx(to2, ho2, \"C1.\", label=\"obs at 90m\")\n",
+    "plt.semilogx(to2, hm2[0], \"C1\", label=\"timflow at 90m\")\n",
+    "plt.title(\"Model Results\")\n",
+    "plt.xlabel(\"time (d)\")\n",
+    "plt.ylabel(\"head change (m)\")\n",
+    "plt.legend()\n",
+    "plt.grid()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Comparison of results"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The performance of `timflow` was compared with MLU (Hemker & Post, 2014), and Kruseman and de Ridder (1970), here abbreviated to K&dR. \n",
+    "\n",
+    "MLU included, compared to `timflow`, three observation wells in the calibration, located at 30, 90, and 215 m from the pumping well. The fit remains poor for the observations at 215 m, which explains the relatively large RMSE. The late-time drawdowns increase more slowly than the computed Theis drawdowns, which is likely caused by leakage (Kruseman and de Ridder, 1970; Hemker & Post, 2014).\n",
+    "\n",
+    "Since `timflow` only includes the first two observation wells (30 and 90 m), the fit to the data is better. However, the estimated parameters are quite similar to those obtained with MLU, while the solution reported by K&dR shows differences."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "t = pd.DataFrame(\n",
+    "    columns=[\"k [m/d]\", \"Ss [1/m]\", \"RMSE [m]\"],\n",
+    "    index=[\"timflow\", \"MLU\", \"K&dR\"],\n",
+    ")\n",
+    "\n",
+    "t.loc[\"timflow\"] = np.append(cal.parameters[\"optimal\"].values, cal.rmse())\n",
+    "t.loc[\"MLU\"] = [63.51428571, 3.58e-5, 0.832886547]\n",
+    "t.loc[\"K&dR\"] = [55.71429, 1.7e-4, 1.292710331]\n",
+    "\n",
+    "t_formatted = t.style.format(\n",
+    "    {\n",
+    "        \"k [m/d]\": \"{:.2f}\",\n",
+    "        \"Ss [1/m]\": \"{:.2e}\",\n",
+    "        \"RMSE [m]\": lambda x: \"-\" if x == \"-\" else f\"{float(x):.3f}\",\n",
+    "    }\n",
+    ")\n",
+    "t_formatted"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {
+    "editable": true,
+    "slideshow": {
+     "slide_type": ""
+    },
+    "tags": []
+   },
+   "source": [
+    "## References\n",
+    "\n",
+    "* Bakker, M. (2013), Semi-analytic modeling of transient multi-layer flow with TTim, Hydrogeol J 21, 935–943, https://doi.org/10.1007/s10040-013-0975-2\n",
+    "* Duffield, G.M. (2007), AQTESOLV for Windows Version 4.5 User's Guide, HydroSOLVE, Inc., Reston, VA.\n",
+    "* Hemker, K. en Post V. (2014) MLU for Windows: well flow modeling in multilayer aquifer systems; [MLU User's guide](https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmicrofem.com%2Fdownload%2Fmlu-user.pdf&data=05%7C02%7CMark.Bakker%40tudelft.nl%7Cad7f16364d2d4fd55dbf08de73832eaa%7C096e524d692940308cd38ab42de0887b%7C0%7C0%7C639075204580287861%7CUnknown%7CTWFpbGZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkFOIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=OBoe8seXZUfoat89Dfr4g6lF%2Bn1FdtXqtp%2F18BMXCn0%3D&reserved=0)\n",
+    "* Kruseman, G.P., De Ridder, N.A. and Verweij, J.M. (1970), Analysis and evaluationof pumping test data, volume 11, International institute for land reclamation and improvement The Netherlands."
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}