{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from myst_nb import glue\n", "%config InlineBackend.rc = {'figure.dpi': 300}\n", "%config InlineBackend.figure_format = 'svg'\n", "import sys\n", "sys.path.append(\"../py_figures\")\n", "import diagrams_2d_chap9 as diagrams_2d\n", "import inter_fig_chap9" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "(content:chap9.5:aquifer_system)=\n", "# Wells in a two-aquifer system" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The solution for steady flow to a well in the bottom aquifer of a two-aquifer system was presented in {numref}`content:chap6.5:well_two_aqu`. In this section, the transient equivalent is presented. The solution is used to investigate the transient head response in a two-aquifer system. When the top aquifer is unconfined, the head response in the top aquifer will be significantly slower than the head response in the bottom aquifer, also referred to as the delayed response of the water table. The solution is quite complicated as it involves solving a system of two coupled differential equations using the Laplace transform method. Although the equations are a bit longer, the Python script of the solution is relatively short.\n", "
      Consider transient flow in a two-aquifer system. The top aquifer has index 0, and the bottom aquifer has index 1. The transmissivity of each aquifer is again approximated as constant in both space and time. The storage coefficients of the top and bottom aquifers are $S_0$ and $S_1$, respectively. The two aquifers are separated by a leaky layer with resistance $c$ ({numref}`Figure %s `). The storage in the leaky layer is neglected. The transient equivalent of the system of two linked differential equations ({eq}`eq_2aqsys6.5`) is\n", "\n", "$$\n", "& \\nabla^2 h_0 = \\frac{S_0}{T_0}\\frac{\\partial h_0}{\\partial t} + \\frac{h_0-h_1}{cT_0}\n", "\n", "& \\nabla^2 h_1 = \\frac{S_1}{T_1}\\frac{\\partial h_1}{\\partial t} + \\frac{h_1-h_0}{cT_1}\n", "$$(eq_2_aqui_sys9.5)\n", "\n", "A well starts pumping with a discharge $Q$ at time $t=0$ in aquifer 1 (the bottom aquifer). The boundary conditions are \n", "\n", "$$\n", "\\lim_{r \\to 0} -2\\pi r Q_{r_0} & = 0 \n", "\n", "\\lim_{r \\to 0} -2\\pi r Q_{r_1} & = Q {\\hskip 2em} t\\ge 0\n", "$$(eq_bctr2aqwell19.5)\n", "\n", "$$\n", "h_0|_{r=\\infty,t} = h_1|_{r=\\infty,t} = 0\n", "$$(eq_bctr2aqwell29.5)\n", "\n", "and the initial condition is\n", "\n", "$$\n", "h_0|_{r,t=0} = h_1|_{r,t=0} = 0\n", "$$(eq_ictr2aqwell9.5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A solution is obtained using Laplace transforms, as explained in {numref}`content:chap5:laplace_transforms`. The solution procedure is a bit more complicated, however, as the system of two linked differential equations needs to be broken down into the solution of two unlinked differential equations. This unlinking is more difficult than for steady flow ({numref}`content:chap6.4:well_sem_con_aqu`) and requires the solution of an eigenvalue problem.\n", "
      The Laplace-transform of the system of differential equations {eq}`eq_2_aqui_sys9.5` is\n", "\n", "$$\n", "\\nabla^2 \\overline{h}_0 = \\frac{S_0 p}{T_0}\\overline{h}_0 + \\frac{\\overline{h}_0-\\overline{h}_1}{cT_0} \\\\\n", "\\nabla^2 \\overline{h}_1 = \\frac{S_1 p}{T_1}\\overline{h}_1 + \\frac{\\overline{h}_1-\\overline{h}_0}{cT_1}\n", "$$(eq_lap_sys9.5)\n", "\n", "which may be written in matrix form as\n", "\n", "$$\n", "\\nabla^2\\left(\n", "\\begin{array}{c}\n", "\\overline{h}_0 \\\\\n", "\\overline{h}_1\n", "\\end{array}\n", "\\right)=\n", "\\left(\n", "\\begin{array}{cc}\n", "\\frac{S_0 c p + 1}{cT_0} & \\frac{-1}{cT_0} \\\\\n", "\\frac{-1}{cT_1} & \\frac{S_1 c p + 1}{cT_1}\n", "\\end{array}\n", "\\right)\n", "\\left(\n", "\\begin{array}{c}\n", "\\overline{h}_0 \\\\\n", "\\overline{h}_1\n", "\\end{array}\n", "\\right)\n", "$$(eq_matdeq9.5)\n", "\n", "The matrix on the right-hand side of the matrix equation is called the system matrix $\\mathbf{A}$ and is written as\n", "\n", "$$\n", "\\mathbf{A}=\n", "\\left(\n", "\\begin{array}{cc}\n", "\\frac{S_0 c p + 1}{cT_0} & \\frac{-1}{cT_0} \\\\\n", "\\frac{-1}{cT_1} & \\frac{S_1 c p + 1}{cT_1}\n", "\\end{array}\n", "\\right) = \\mathbf{V} \\mathbf{W} \\mathbf{V}^{-1}\n", "$$(eq_matrix_A9.5)\n", "\n", "where $\\mathbf{W}$ is a diagonal matrix with the eigenvalues of matrix $\\mathbf{A}$ on the diagonal\n", "\n", "$$\n", "\\mathbf{W} = \n", "\\left(\n", "\\begin{array}{cc}\n", "w_0 & 0 \\\\\n", "0 & w_1\n", "\\end{array}\n", "\\right)\n", "$$(eq_matrix_W9.5)\n", "\n", "and $\\mathbf{V}$ is a matrix with the corresponding eigenvectors \n", "$\\vec{u}=\\left(\\begin{array}{c} u_0 \\\\ u_1 \\end{array}\\right)$ and \n", "$\\vec{v}=\\left(\\begin{array}{c} v_0 \\\\ v_1 \\end{array}\\right)$ as its columns\n", "\n", "$$\n", "\\mathbf{V} = \n", "\\left(\n", "\\begin{array}{cc}\n", "u_0 & v_0 \\\\\n", "u_1 & v_1\n", "\\end{array}\n", "\\right)\n", "$$(eq_col_V9.5)\n", "\n", "The boundary conditions in the Laplace domain are\n", "\n", "$$\n", "\\begin{split}\n", "\\lim_{r \\to 0} -2\\pi r \\overline{Q}_{r_0} &= 0 \\\\\n", "\\lim_{r \\to 0} -2\\pi r \\overline{Q}_{r_1} &= \\frac{Q}{p}\n", "\\end{split}\n", "$$(eq_bctrwell2aqlap19.5)\n", "\n", "$$\n", "\\overline{h}_0|_{r=\\infty} = \\overline{h}_1|_{r=\\infty} = 0\n", "$$(eq_bctrwell2aqlap29.5)\n", "\n", "The solution in the Laplace domain is\n", "\n", "$$\n", "\\left(\n", "\\begin{array}{c}\n", "\\overline{h}_0 \\\\\n", "\\overline{h}_1\n", "\\end{array}\n", "\\right) = -\\frac{Q}{2 \\pi T_1 p}\\left[ b_0\\text{K}_0\\left(r\\sqrt{w_0}\\right) \n", "\\left(\n", "\\begin{array}{c}\n", "u_0 \\\\\n", "u_1\n", "\\end{array}\n", "\\right) + b_1\\text{K}_0\\left(r\\sqrt{w_1}\\right) \\left(\n", "\\begin{array}{c}\n", "v_0 \\\\\n", "v_1\n", "\\end{array}\n", "\\right)\\right]\n", "$$(eq_sol_lap9.5)\n", "\n", "where $b_0$ and $b_1$ are the solution of the following system of two linear equations\n", "\n", "$$\n", "\\left(\n", "\\begin{array}{cc}\n", "u_0 & v_0 \\\\\n", "u_1 & v_1\n", "\\end{array}\n", "\\right)\n", "\\left(\n", "\\begin{array}{c}\n", "b_0 \\\\\n", "b_1\n", "\\end{array}\\right) = \\left(\n", "\\begin{array}{c}\n", "0 \\\\\n", "1\n", "\\end{array}\\right)\n", "$$(eq_b_0_b_19.5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The Python implementation of the solution makes use of the `eig` function of the `lingalg` subpackage of `numpy`. The `eig` function returns the array $[w_0, w_1]$ of eigenvalues and the two-dimensional array $\\mathbf{V}$ with the eigenvectors.\n", "
      In the example below, the top aquifer is unconfined so that the storage coefficient $S_0$ represents the phreatic storage of the top aquifer. The head is computed and plotted in an observation well near the well ({numref}`Figure %s `). The head responds very quickly in the bottom aquifer, but it takes approximately one day before the head in the unconfined top aquifer starts to respond. The head in the bottom aquifer starts to go down again when the head in the top aquifer starts to go down. This delayed response of the water table can be seen especially well in a semi-log plot. For early times (in this case, for the first day), the head response is similar to a Hantush well, as the head in the top aquifer varies very little at early time." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# parameters\n", "T0 = 100 # transmissivity top aquifer, m^2/d\n", "S0 = 0.1 # specific yield top aquifer, -\n", "T1 = 200 # transmissivity bottom aquifer, m^2/d\n", "S1 = 0.001 # storage coefficient bottom aquifer, -\n", "c = 500 # resistance of leaky layer, d\n", "Q = 800 # discharge of well, m^3/d\n", "rw = 0.3 # radius of well, m" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "# Hantush solution\n", "from scipy.integrate import quad\n", "\n", "def integrand(tau, r, T, lab):\n", " return 1 / tau * np.exp(-tau - r ** 2 / (4 * lab ** 2 * tau))\n", "\n", "def hantush(r, t, T, S, c, Q):\n", " lab = np.sqrt(T * c)\n", " u = S * r ** 2 / (4 * T * t)\n", " F = quad(integrand, u, np.inf, args=(r, T, lab))[0]\n", " return -Q / (4 * np.pi * T) * F\n", "hantushvec = np.vectorize(hantush)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "from scipy.special import factorial as fac\n", "def stehfest_coef(M):\n", " assert M % 2 == 0, 'M must be even'\n", " M2 = int(M / 2)\n", " V = np.zeros(M)\n", " for k in range(1, M + 1):\n", " z = np.zeros(min(k, M2) + 1)\n", " for j in range(int((k + 1) / 2), min(k, M2) + 1):\n", " z[j] = j ** M2 * fac(2 * j) / \\\n", " (fac(M2 - j) * fac(j) * fac(j - 1) * fac(k - j) * fac(2 * j - k))\n", " V[k - 1] = (-1) ** (k + M2) * np.sum(z)\n", " return V" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "def stehfest(x, t, func, M=12, **kwargs):\n", " t = np.atleast_1d(t)\n", " x = np.atleast_1d(x)\n", " f = np.zeros((len(x), len(t)))\n", " V = stehfest_coef(M)\n", " for i in range(len(t)):\n", " p = np.arange(1, M + 1) * np.log(2) / t[i]\n", " for j in range(len(x)):\n", " fbar = func(p, x[j], **kwargs)\n", " f[j, i] = np.log(2) / t[i] * np.sum(V * fbar)\n", " return f" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# solution\n", "from scipy.special import k0\n", "def hmaqbar(plist, r, layer):\n", " hbar = np.zeros(len(plist))\n", " for i, p in enumerate(plist):\n", " A = np.array([[(p * c * S0 + 1) / (c * T0), -1 / (c * T0)],\n", " [-1 / (c * T1), (p * c * S1 + 1) / (c * T1)]])\n", " w, v = np.linalg.eig(A)\n", " b = np.linalg.solve(v, [0, 1])\n", " hbar[i] = -Q / (2 * np.pi * T1 * p) * (\n", " b[0] * v[layer, 0] * k0(r * np.sqrt(w[0])) + \n", " b[1] * v[layer, 1] * k0(r * np.sqrt(w[1])))\n", " return hbar\n", "\n", "r = 250\n", "t = np.logspace(-3, 2, 100)\n", "h0 = stehfest(r, t, hmaqbar, layer=0) # stehfest function from Section 5.5\n", "h1 = stehfest(r, t, hmaqbar, layer=1)\n", "h1_hantush = hantushvec(r, t, T1, S1, c, Q) # hantushvec from Section 9.3" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2025-10-08T20:39:48.849407\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic plot semilog\n", "plt.semilogx(t, h0[0])\n", "plt.semilogx(t, h1[0])\n", "plt.semilogx(t, h1_hantush);" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n \n \n \n \n 2025-10-08T20:39:49.491580\n image/svg+xml\n \n \n Matplotlib v3.10.6, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "application/papermill.record/text/plain": "
" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_two_aquifers" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2025-10-08T20:39:49.903836\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.10.6, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "fig = plt.figure(figsize=diagrams_2d.figdef2)\n", "plt.subplot(121)\n", "plt.plot(t, h0[0], label='aquifer 0')\n", "plt.plot(t, h1[0], label='aquifer 1')\n", "plt.plot(t, h1_hantush, '--', label='aquifer 1, Hantush')\n", "plt.xlabel('time (days)')\n", "plt.ylabel('head (m)')\n", "plt.title('(a)', loc='left')\n", "plt.legend()\n", "plt.grid()\n", "plt.subplot(122)\n", "plt.semilogx(t, h0[0], label='aquifer 0')\n", "plt.semilogx(t, h1[0], label='aquifer 1')\n", "plt.semilogx(t, h1_hantush, '--', label='aquifer 1, Hantush')\n", "plt.xlabel('time (days)')\n", "plt.ylabel('head (m)')\n", "plt.title('(b)', loc='left')\n", "plt.legend()\n", "plt.grid()\n", "plt.tight_layout()\n", "glue('fig_two_aquifers', fig, display=False);" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/javascript": "'use strict';\n(function(root) {\n function now() {\n return new Date();\n }\n\n const force = true;\n\n if (typeof root._bokeh_onload_callbacks === \"undefined\" || force === true) {\n root._bokeh_onload_callbacks = [];\n root._bokeh_is_loading = undefined;\n }\n\nconst JS_MIME_TYPE = 'application/javascript';\n const HTML_MIME_TYPE = 'text/html';\n const EXEC_MIME_TYPE = 'application/vnd.bokehjs_exec.v0+json';\n const CLASS_NAME = 'output_bokeh rendered_html';\n\n /**\n * Render data to the DOM node\n */\n function render(props, node) {\n const script = document.createElement(\"script\");\n node.appendChild(script);\n }\n\n /**\n * Handle when an output is cleared or removed\n */\n function handleClearOutput(event, handle) {\n function drop(id) {\n const view = Bokeh.index.get_by_id(id)\n if (view != null) {\n view.model.document.clear()\n Bokeh.index.delete(view)\n }\n }\n\n const cell = handle.cell;\n\n const id = cell.output_area._bokeh_element_id;\n const server_id = cell.output_area._bokeh_server_id;\n\n // Clean up Bokeh references\n if (id != null) {\n drop(id)\n }\n\n if (server_id !== undefined) {\n // Clean up Bokeh references\n const cmd_clean = \"from bokeh.io.state import curstate; print(curstate().uuid_to_server['\" + server_id + \"'].get_sessions()[0].document.roots[0]._id)\";\n cell.notebook.kernel.execute(cmd_clean, {\n iopub: {\n output: function(msg) {\n const id = msg.content.text.trim()\n drop(id)\n }\n }\n });\n // Destroy server and session\n const cmd_destroy = \"import bokeh.io.notebook as ion; ion.destroy_server('\" + server_id + \"')\";\n cell.notebook.kernel.execute(cmd_destroy);\n }\n }\n\n /**\n * Handle when a new output is added\n */\n function handleAddOutput(event, handle) {\n const output_area = handle.output_area;\n const output = handle.output;\n\n // limit handleAddOutput to display_data with EXEC_MIME_TYPE content only\n if ((output.output_type != \"display_data\") || (!Object.prototype.hasOwnProperty.call(output.data, EXEC_MIME_TYPE))) {\n return\n }\n\n const toinsert = output_area.element.find(\".\" + CLASS_NAME.split(' ')[0]);\n\n if (output.metadata[EXEC_MIME_TYPE][\"id\"] !== undefined) {\n toinsert[toinsert.length - 1].firstChild.textContent = output.data[JS_MIME_TYPE];\n // store reference to embed id on output_area\n output_area._bokeh_element_id = output.metadata[EXEC_MIME_TYPE][\"id\"];\n }\n if (output.metadata[EXEC_MIME_TYPE][\"server_id\"] !== undefined) {\n const bk_div = document.createElement(\"div\");\n bk_div.innerHTML = output.data[HTML_MIME_TYPE];\n const script_attrs = bk_div.children[0].attributes;\n for (let i = 0; i < script_attrs.length; i++) {\n toinsert[toinsert.length - 1].firstChild.setAttribute(script_attrs[i].name, script_attrs[i].value);\n toinsert[toinsert.length - 1].firstChild.textContent = bk_div.children[0].textContent\n }\n // store reference to server id on output_area\n output_area._bokeh_server_id = output.metadata[EXEC_MIME_TYPE][\"server_id\"];\n }\n }\n\n function register_renderer(events, OutputArea) {\n\n function append_mime(data, metadata, element) {\n // create a DOM node to render to\n const toinsert = this.create_output_subarea(\n metadata,\n CLASS_NAME,\n EXEC_MIME_TYPE\n );\n this.keyboard_manager.register_events(toinsert);\n // Render to node\n const props = {data: data, metadata: metadata[EXEC_MIME_TYPE]};\n render(props, toinsert[toinsert.length - 1]);\n element.append(toinsert);\n return toinsert\n }\n\n /* Handle when an output is cleared or removed */\n events.on('clear_output.CodeCell', handleClearOutput);\n events.on('delete.Cell', handleClearOutput);\n\n /* Handle when a new output is added */\n events.on('output_added.OutputArea', handleAddOutput);\n\n /**\n * Register the mime type and append_mime function with output_area\n */\n OutputArea.prototype.register_mime_type(EXEC_MIME_TYPE, append_mime, {\n /* Is output safe? */\n safe: true,\n /* Index of renderer in `output_area.display_order` */\n index: 0\n });\n }\n\n // register the mime type if in Jupyter Notebook environment and previously unregistered\n if (root.Jupyter !== undefined) {\n const events = require('base/js/events');\n const OutputArea = require('notebook/js/outputarea').OutputArea;\n\n if (OutputArea.prototype.mime_types().indexOf(EXEC_MIME_TYPE) == -1) {\n register_renderer(events, OutputArea);\n }\n }\n if (typeof (root._bokeh_timeout) === \"undefined\" || force === true) {\n root._bokeh_timeout = Date.now() + 5000;\n root._bokeh_failed_load = false;\n }\n\n const NB_LOAD_WARNING = {'data': {'text/html':\n \"
\\n\"+\n \"

\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"

\\n\"+\n \"
    \\n\"+\n \"
  • re-rerun `output_notebook()` to attempt to load from CDN again, or
    • \\n\"+\n \"
  • use INLINE resources instead, as so:
    • \\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"\\n\"+\n \"
\"}};\n\n function display_loaded(error = null) {\n const el = document.getElementById(null);\n if (el != null) {\n const html = (() => {\n if (typeof root.Bokeh === \"undefined\") {\n if (error == null) {\n return \"BokehJS is loading ...\";\n } else {\n return \"BokehJS failed to load.\";\n }\n } else {\n const prefix = `BokehJS ${root.Bokeh.version}`;\n if (error == null) {\n return `${prefix} successfully loaded.`;\n } else {\n return `${prefix} encountered errors while loading and may not function as expected.`;\n }\n }\n })();\n el.innerHTML = html;\n\n if (error != null) {\n const wrapper = document.createElement(\"div\");\n wrapper.style.overflow = \"auto\";\n wrapper.style.height = \"5em\";\n wrapper.style.resize = \"vertical\";\n const content = document.createElement(\"div\");\n content.style.fontFamily = \"monospace\";\n content.style.whiteSpace = \"pre-wrap\";\n content.style.backgroundColor = \"rgb(255, 221, 221)\";\n content.textContent = error.stack ?? error.toString();\n wrapper.append(content);\n el.append(wrapper);\n }\n } else if (Date.now() < root._bokeh_timeout) {\n setTimeout(() => display_loaded(error), 100);\n }\n }\n\n function run_callbacks() {\n try {\n root._bokeh_onload_callbacks.forEach(function(callback) {\n if (callback != null)\n callback();\n });\n } finally {\n delete root._bokeh_onload_callbacks\n }\n console.debug(\"Bokeh: all callbacks have finished\");\n }\n\n function load_libs(css_urls, js_urls, callback) {\n if (css_urls == null) css_urls = [];\n if (js_urls == null) js_urls = [];\n\n root._bokeh_onload_callbacks.push(callback);\n if (root._bokeh_is_loading > 0) {\n console.debug(\"Bokeh: BokehJS is being loaded, scheduling callback at\", now());\n return null;\n }\n if (js_urls == null || js_urls.length === 0) {\n run_callbacks();\n return null;\n }\n console.debug(\"Bokeh: BokehJS not loaded, scheduling load and callback at\", now());\n root._bokeh_is_loading = css_urls.length + js_urls.length;\n\n function on_load() {\n root._bokeh_is_loading--;\n if (root._bokeh_is_loading === 0) {\n console.debug(\"Bokeh: all BokehJS libraries/stylesheets loaded\");\n run_callbacks()\n }\n }\n\n function on_error(url) {\n console.error(\"failed to load \" + url);\n }\n\n for (let i = 0; i < css_urls.length; i++) {\n const url = css_urls[i];\n const element = document.createElement(\"link\");\n element.onload = on_load;\n element.onerror = on_error.bind(null, url);\n element.rel = \"stylesheet\";\n element.type = \"text/css\";\n element.href = url;\n console.debug(\"Bokeh: injecting link tag for BokehJS stylesheet: \", url);\n document.body.appendChild(element);\n }\n\n for (let i = 0; i < js_urls.length; i++) {\n const url = js_urls[i];\n const element = document.createElement('script');\n element.onload = on_load;\n element.onerror = on_error.bind(null, url);\n element.async = false;\n element.src = url;\n console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n document.head.appendChild(element);\n }\n };\n\n function inject_raw_css(css) {\n const element = document.createElement(\"style\");\n element.appendChild(document.createTextNode(css));\n document.body.appendChild(element);\n }\n\n const js_urls = [\"https://cdn.bokeh.org/bokeh/release/bokeh-3.7.3.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-3.7.3.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-3.7.3.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-3.7.3.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-mathjax-3.7.3.min.js\"];\n const css_urls = [];\n\n const inline_js = [ function(Bokeh) {\n Bokeh.set_log_level(\"info\");\n },\nfunction(Bokeh) {\n }\n ];\n\n function run_inline_js() {\n if (root.Bokeh !== undefined || force === true) {\n try {\n for (let i = 0; i < inline_js.length; i++) {\n inline_js[i].call(root, root.Bokeh);\n }\n\n } catch (error) {throw error;\n }} else if (Date.now() < root._bokeh_timeout) {\n setTimeout(run_inline_js, 100);\n } else if (!root._bokeh_failed_load) {\n console.log(\"Bokeh: BokehJS failed to load within specified timeout.\");\n root._bokeh_failed_load = true;\n } else if (force !== true) {\n const cell = $(document.getElementById(null)).parents('.cell').data().cell;\n cell.output_area.append_execute_result(NB_LOAD_WARNING)\n }\n }\n\n if (root._bokeh_is_loading === 0) {\n console.debug(\"Bokeh: BokehJS loaded, going straight to plotting\");\n run_inline_js();\n } else {\n load_libs(css_urls, js_urls, function() {\n console.debug(\"Bokeh: BokehJS plotting callback run at\", now());\n run_inline_js();\n });\n }\n}(window));", "application/vnd.bokehjs_load.v0+json": "" }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "\n", "
\n" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": "(function(root) {\n function embed_document(root) {\n const docs_json = {\"00010d70-bd92-4dfe-938d-35658ff947ac\":{\"version\":\"3.7.3\",\"title\":\"Bokeh Application\",\"roots\":[{\"type\":\"object\",\"name\":\"Column\",\"id\":\"p1219\",\"attributes\":{\"children\":[{\"type\":\"object\",\"name\":\"Row\",\"id\":\"p1217\",\"attributes\":{\"children\":[{\"type\":\"object\",\"name\":\"Figure\",\"id\":\"p1006\",\"attributes\":{\"width\":380,\"height\":300,\"x_range\":{\"type\":\"object\",\"name\":\"Range1d\",\"id\":\"p1016\",\"attributes\":{\"start\":-5,\"end\":105}},\"y_range\":{\"type\":\"object\",\"name\":\"Range1d\",\"id\":\"p1017\",\"attributes\":{\"start\":-1.5,\"end\":0.05}},\"x_scale\":{\"type\":\"object\",\"name\":\"LinearScale\",\"id\":\"p1018\"},\"y_scale\":{\"type\":\"object\",\"name\":\"LinearScale\",\"id\":\"p1019\"},\"title\":{\"type\":\"object\",\"name\":\"Title\",\"id\":\"p1009\",\"attributes\":{\"text\":\"(a)\"}},\"outline_line_color\":\"Black\",\"renderers\":[{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1049\",\"attributes\":{\"data_source\":{\"type\":\"object\",\"name\":\"ColumnDataSource\",\"id\":\"p1003\",\"attributes\":{\"selected\":{\"type\":\"object\",\"name\":\"Selection\",\"id\":\"p1004\",\"attributes\":{\"indices\":[],\"line_indices\":[]}},\"selection_policy\":{\"type\":\"object\",\"name\":\"UnionRenderers\",\"id\":\"p1005\"},\"data\":{\"type\":\"map\",\"entries\":[[\"t\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"/Knx0k1iUD+sQAH/j2dSPx1M4YKcrFQ/23ENS1E5Vz+VHX+hhBZaP+m9zc4jTl0/6Ggnwap1YD9cbjs8UH1iP1youHILxWQ/kmI1ncNUZz/z1Ut3WTVqPwU77QTGcG0/s7+VkR6JcD9a0iMuKpNyP+6UzEKX3XQ/odY9X1Zwdz/2Xwq9UlR6P0uiWimRk30/yBFIX6mcgD/yoRvzHamCP/1kPRVA9oQ/UYd8twmMhz+aucqdcHOKP4Z3dWyFto0/ubhpRUuwkD/O+aepK7+SP5DAUwwGD5U/FHx0zN2nlz8qxc9Es5KaP7Zp1v6i2Z0/IipGXwTEoD9ACHJwU9WiPyDUgErpJ6U/9j/WxNLDpz+ShY/dGrKqP46WTxHq/K0/kx1JyNTXsD8COEdmleuyP3+AXvLpQLU/axeAx+jftz+8WrOTp9G6P0TO7NRaIL4/qrL+m7zrwD/xWhmq8QHDP/iKryYIWsU/YzZ++x/8xz80PhiTWfHKP+/X83r1Q84/apcT9rv/0D8o1f5aaBjTP5TNXwpEc9U/nvYKiHgY2D9NAM8HMRHbP9O15DS6Z94/oy5V8tIT4T8fyDKY+S7jP8RnhMCdjOU/Eg6PlPI06D8uhRweLjHrP0nqeTSpi+4/uLaxrAEo8T84PhWBpUXzP/7uW2wVpvU/zcWhSI5R+D+XAnoCUVH7P/+8qKvCr/4/UHA4QUg8AUA2Vis1bFwDQOyfTjGrvwVAETEJzEtuCECFPZXhmXELQC2AocwG1A5AbMUZzKZQEUA6bx/UTXMTQH2P7jJf2RVAc2S6RiuLGEBgyFDoCJIbQKbWz8l1+B5AonCnaR1lIUCdVMF9SoojQG3c95Qx8yVAmq3Z4CyoKEB7QcRDnrIrQLD52tUPHS9AXaRUNqx5MUBfagZSYqEzQO7gUHsiDTZA3cq6wlDFOEC0kTwhWtM7QIv/pSPVQT9AjjK2TlOOQUBN2QlxlbhDQJ1kCgoyJ0ZAjyPhFJfiSEA/KzyuPPRLQEYiUObFZk9AXLSCzxKjUUDwuwz7489TQKfOX2VgQVZAAAAAAAAAWUA=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h0_1\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"d/Rr+qLxHD1ZkXj8uQEwPaxpz17dgD89AuuD32mkSj1v4lGPd8FRPTIjgEotjEk9J+ngQHJLVL2idbuxcE18veki0rMIIJO9GYQHrzeQo72TVeLAwVOwvVR5NsYfLra9Dzu0iE+Qtr2YXAGlpcKkvdnwkYMNwbU9cwZbSfVT0z3GLPfHE8XiPYf3Yhy0VOs97yG2lv/h7T38dUzeFmLgPRSfPNzHs/C9yENX+OumFL6mCjyWmFYuvkWKU4ziO0O+sslb/YKcVr7x14MfVhppvht3WltodHq+m0lm+wOJir4d0AAYcmeZvkEblogiSqe+eSUCkEeEtL6dUbSngm3Bvn8dbeqapMy+YyQzmf/W1r4/736Kn7jhvtD12Db80uq+9CvXLObY876wstxPJsX8vp+cUd1ydgS/PbeDmx2cDL/0NrBYnK8Tv3k4kr4Isxq/35a31zTdIb946e3cN5wnvzQtGIQP2S6/VvOt/SPwM78+6idorIQ5v/4l1//gLUC/KGwI66tWRL+V77P7DVxJv4WB7tjtYE+/uu6u6vhFU7/jCQZ214NXv3LTfHNsgVy/hIAF2x0sYb+bDGlfCZJkv8FV8r6YgWi/t8E1ZkALbb/5Ez4HdyBxv/gdvA8dG3S/RmffOFSAd7/gkO9q/Vt7v8goCXZLu3+/0YBh4X9Wgr9ANz8G2CCFvwJXHa0JRoi/VbEloqnQi78lVOCyrsyPvydxje7LI5K//LTic0aolL+aHtcFO3yXv8HLsKQDqZq/g6RqG944nr9WVnxicBuhv7ofGZpyV6O/U1PYH7TWpb+DlYq/r5+ov1ahkHICuau/yhoX0Uwpr79EXH1+iHuxvxP19wJHlLO/zkLo9Mzhtb+5DV2wvWa4v95wN1piJbu/EFDycZgfvr9jmaQUYavAv/s7GI7cZcK/V822A2I/xL+HGmr/xzfGvwoDsdKWTsi/UGxLYwmDyr9/mZE/ENTMvwLsOwZXQM+/Ps+5Aibj0L8IsNrrFDLSv5HHElUBjNO/qtGtWObv1L9UNpiet1zWv+uNbOhn0de/HONBSu9M2b8=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h1_1\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"SOUjbRwgAj7pgbMmp08JPkVbORSZrAE+Y//9HhwIE75qkMefzPU4vqAMBQY8s1C+MyhPUh2PYL5ltD7AimppvkAayxMQwGu+enEjV8hcVb52MSlZEX10PturMDAi1ZA+aGAoOwbTnz7GQuADfsilPu6XnfBflKQ+j5HjPyvoiD7f81lpB7Spvtyks5zQjMO+BsrFVC021L4A8ErOCADkvliK5Pj2TPW+jO5CbrITCL/N++FCuhUbv+k4uEjMAy2/J1DIqIAgPb/n3Jufsl5Lv+Xq3lR6L1i/0aCph4k4ZL/v6HTM4hdwvyUIITEGhXi/47M8fUL3gb/i/yLv8m2Jvwq9E2tbcpG/n4ZjjDVIl785/45E0kyev9RrRv+fRqO/31Ut81UHqL+v1BnOFWetv/nBVRw6sLG//j28RSb1tL/Et4MvDXy4v4MECe4yPby/JVoFl84XwL/vN/A+qCTCvz/+d7rGP8S/BYPCUIZjxr8lXqd8HIrIv0aeP5C0rcq/3NHyHYzIzL+u6WGMD9XOv4a7ZJb7ZtC/BKRsfTJX0b9ZRclwATnSvy2JRX6QCtO/GoxTZ2nK079XWhrvhHfUvwcVFdJVEdW/j5/MJNCX1b/GDwTAawvWvwGTaHggbda/PhoKOlu+1r8/B/HE6wDXv7nlWLjrNte/KkmznqBi179/msW9W4bXv7TmBgdbpNe/t2cqpa6+17/yp/IKJtfXv/HYvSVG79e/K/Zb1kgI2L8TYGk9JCPYv/fl4WeXQNi/4AJDbDlh2L+RhrFMiIXYvyLVlbL1rdi/WFzEgPDa2L+6eHgP6wzZv7aJCo5eRNm/dGRWM8yB2b+1gM4AvcXZvwFGxarAENq/7xNyDGxj2r+7W2pWV77av5GFpAgcItu/TC4Xv1KP27+yYDXFkAbcvympC2tliNy/wHc0NFcV3b+785HO4K3dv1g8dwZuUt6/0BgUzFgD378xNfRZ5sDfv3a+32KiReC/q8uFhESx4L/zSHrSViPhv/VEbwnJm+G/cYx8gH0a4r/Lism3SZ/ivzdgK1X3KeO/CJxjgEW6478=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h0_2\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"L4ZNl5jxUj06NiJTtfNkPUt+evjvnXQ9zAgWKhRvgT16bKW0hjuHPQCWc2N5soA9gOuh78edir1CjPMqI4myvX016zRAC8m9olIjfyqd2b1cWm3ctF/lvRisRMK8B+29RSteU/yE7b29SGyV6hfbvW6uoqBhjuw9acXi5l1VCT659z0r75cYPlUFxmIX5yE+ZAhxY7GSIz5FQqeND3cVPrjuSSm63SW+kEg/kDUKS76uZ3iqLNxjvvm4b2AiLnm+sh+YQpqZjb4xfPijDG6gvpS2P1VSULG++KHne5Jdwb4xfL/I0p/QvlmOVQquet6+E7kLli3Z6r6HXte53s32vknhIVHevAK/P9uJFCbhDb/XkFzLHy4XvyEKwgWFiiG/Y2xoRwH0Kb9ztmN2nc4yv2pE1cVfvzq/toDySrixQr+1fj+/Q7hJv3Q63O6xb1G/Rf9R5jFTV79R+zJn7tBev5KW1jl1H2S/rZ5B/vb/ab9yoiZ3IqFwv853uXQAE3W/jrUVlUF5er8hPVQtxH2Av3tvu562YoS/ZqBS6CwDib+8xvPQm3iOv/smM1H7bpK/CiDl6aknlr/JMuecvHSav42iUNmpZJ+/Xt1cKNiCor85ZPNaxLKlv8xtbH+MSKm/Otq7uxZKrb+GMIDBTt6wv/4FECJLUrO//Jbv184Ctr9KHMlnNPG4v0wOoc1cHry/F+9126qKv79Zi8Oc/5rBv0agkoTaj8O/eo0k1E+jxb/TbAU/hdTHv4nX1+5dIsq/5wzD/n2LzL/6IBZGUA7Pv0tvTj6H1NC/MdJRvOUs0r88XvGAP4/TvxQC9kuJ+tS/M0PtMLdt1r/IMjHJwufXv5XoEsewZ9m/b4ZKp5Xs2r8qsYdOmXXcvzj4g3/5Ad6/STBzJAuR378E5H5DHZHgv5Xbe0uFWuG/0OsBvIkk4r+fH4xEAO/iv4ORpQ7HueO/tjQjjcOE5L9vf/pN4U/lvzGD8OIQG+a/lPAT5kbm5r//Kvgfe7Hnv/TTR9KnfOi/YBUyIslH6b8JT4yh3BLqv7Er7PDg3eq/Oti2dtWo678=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h1_2\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"HsJ4CCogAj7eYosWyU8JPktUDbDmrAE+9fhnlMsHE74Kyr1dp/U4vumk8DMts1C+el2qFRSPYL4EJWL0hGppvoqpefYbwGu+Kyq6bz1dVb7PGGLSxHx0PtzH40f71JA+xvQpycPSnz5lld+aTcilPrbPJvAklKQ+dIn/sE7niD6cXnceILSpvqMS2QPGjMO+Oq+sAxs21L7loFI7+P/jvvfQ58fsTPW+NKG3LLATCL/gHd72vxUbvzHFPmbZAy2/6G63cpQgPb+lrJgvzF5Lv/t/1IiYL1i/wgHu/Ko4ZL/pMV/0BRhwv2Xghb5MhXi/u3eJa4b3gb+Gs9zycG6Jv/+h0GTMcpG/4LCA8flIl7/zjI8rHk6ev/ZP4E6xR6O/Wbzbkw0JqL9sal/UyWmtvxlfKRBQsrG/S/xOUU/4tL/nU5mIwoC4vyp1CtsaRLy/Uq4EcMwcwL8spuykxivCv6slqf7MScS/BRwcQnlxxr/1HFrITZ3IvwrO6+jUx8q/jJHeVMDrzL+Jb+uoCATPvx6cBbIGhtC/tY/elNp/0b8vuA/xxm3Sv7IWS+d1TtO/vM5eEgMh1L+VXwyZB+XUvx4P0s+imtW/MPuaBn5C1r8wYtsyyt3Wv7fPiGM2bte/Iidugd71178pHce4MnfYvwALY/PY9Ni/wTXe2Ypx2b8hFjmc8+/ZvzjoRwiRctq/aF2fF5v72r/6nQEN9Izbv3855MUhKNy/web5Vk/O3L9Ty+jfVIDdvxO4btnCPt6/zISXU+4J379oT7AR/OHfv+wHfL10Y+C/bggIk0nc4L8q01eqXFvhv4tCeZqB4OG/D3naw4Br4r93P0usGPziv8tYmYP/keO/8iKA5eQs5L/T0JzAc8zkv0ddoVVUcOW/uPevJC4Y5r8qQwu2qcPmv2hcsyFycue/KmqJTDYk6L9htHTUqdjovz9eXKuFj+m/45gxbYhI6r8gBsRzdgPrv/na+r4ZwOu/RmOetkF+7L87EBjUwj3tv6OzLT92/u2/3rCYZDnA7r8JpsaM7YLvv/txk7s7I/C/q5fPfF+F8L8=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h0_3\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"xWAnuvhblj2H2wFSTKqoPUiS8Ho9Kbg9N7P7a+Y/xD2CHaUtmmPKPVqFNWuC3cA91UnrEW8v0713jFswUsX3veb42DrwXg++Bj9F0eGzH76O5WBFQyQqvhfc3C+MYTG+IzdJ0e2mML7TjF8Ie0sSvnsr/I2nhzk+6xRvls83Uz6ZdoCfz/VhPurw+GfNpWk+8iRl7aBuaz4tJ6Ql8zBcPvxPA0SEeG++oY4JAnookr4QxJGexYOpvoPJf1YfOr++7z27WG/k0b7mF6YJpIfjvkfofxoRWvS+9QJ0qEg+BL8CboSsBTsTv68dcHeIeiG/GXHMRxh6Lr/PfNNozY85v65EVmPQr0S/nIn59GUzUL+2JEsJUZ9Yv8wgKS2FM2K/Ll/ChIo8ar+uJ19XmXlyv/fWE0USeHm/XcTMtpw2gb/pKbY3StmGvzumND4K042/8OaSbgwqk7+YmmWYbUaYv/XGPrJPVp6/7zA2zQ64or+EO37WI9Smv7eXeXI+iKu/n/0mrbhtsL+3zcJoVWmzv03yuOkvuLa/fzl38iRaur/h/fsJoE2+v3PRXN3MR8G/zxtlddWNw7+VtbJ7G/bFvxUlDuFVfci/oFP13+kfy79SDWCHGdrNv5EG2osZVNC/RZPPB17D0b8SjjKmSTnTv5nrYYiDtNS/08uJ7+4z1r+ZOYGdqbbXv3nfuZsFPNm/EXTGSIDD2r9cN3CTuEzcv2oRnOFl192/VWDik1Bj378VjjU9Jnjgvx5cI40aP+G/mUTBs3UG4r/F0zQ8Ks7ivwTT/GosluO/p+wXBnJe5L90gykw8iblv9IBC0Sl7+W/xbV1r4S45r9j+9PPioHnv52PPc6ySui/i0aRgfgT6b/h/wxSWN3pv1PBOCHPpuq/q2U+OVpw679y07w79znsv3GgBBWkA+2/9ovE8V7N7b+kq4o2Jpfuv+vLR3r4YO+/pWJOPmoV8L9m6+iQXHrwv/pEeLhS3/C/ayTvRkxE8b+WVWnZSKnxv7/KVxpIDvK/xSDpu0lz8r+J0Yl6Tdjyv0Ra9xlTPfO/4kasZVqi878=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"h1_3\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"cDqZRREhAj6Jcio5AFIJPhKrm+XbsQE+yEpguMwCE76s5jhFcfM4vhEz1K1cslC+tDmtYKiOYL56+CLnhGppvjstEOSPwWu+neXuoeRmVb6HuOXPO3d0PhknejBp0pA+aTLggMLOnz5VFawrucWlPmxZJFGgkaQ+QeY2gs7iiD6kuY/a17GpvjoOpcXGisO+hQUzlTE01L6+3zOMnf7jvgv6zMFDTPW+NjmgjK4TCL8Kj8kITRYbv9MNUHzdBC2/98flLPwhPb9NBna7hWBLv/IZiymQMVi/vrtIqMk6ZL8/DJJuMhpwv8zHgrqPiXi/o9X6bIj7gb/2olZ4snWJv4C7NPIjeZG/QO3GXLZTl7+gXCQYxF+ev9cyVSzQVaO/dRLWohkfqL9m/sq3cYutv7LAvJh5y7G//lCgYzkdtb9Dfaai+LW4v/x6PKGRj7y/sOBLY4FRwL/SAOPIXXTCvyhZaM10rMS/S9cLxeH1xr92ewAz6EzJv87CnuIUrsu/hs2YrFsWzr8IpZ59mEHQvwFkfPZOedG/tUbQ4Kax0r9ZZAuYSerTv6UK8S4yI9W/Dav2V6Zc1r/uHl84KpfXv+yaoL9v09i//jazvkIS2r8HDQClc1Tbv1YMVCrDmty/aUElONHl3b9d2Fv9EDbfv62sr57hReC/IucGcnvz4L/Lc4xtx6PhvwRk/bSlVuK/z74xi+oL47+zl6CXY8Pjv1laMFbcfOS/xBI5WiE45b/ayJdeAvXlv+f030xTs+a/V7Frhuxy578zMaq3qjPovzCTXHNu9ei/GorXuxu46b/ONrSUmXvqv5FZrqTRP+u/ABOp5a8E7L8P93RnIsrsvze/QxcZkO2/8uLlk4VW7r/uHMACWx3vvzIVX+iN5O+/gG0OBApW8L9GVwoh8rnwv13uRjr7HfG/J64QsSGC8b+PgTRMYubxv19cki26SvK/WdGcxyav8r/KlrrTpRPzv8HDDUw1ePO/20YXZNPc878xI+l/fkH0v9vuITM1pvS/hzU/N/YK9b+zqTlrwG/1v/7gxcyS1PW/NxpUd2w59r8=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}]]}}},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1050\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1051\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1046\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1047\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1048\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1060\",\"attributes\":{\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1061\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1062\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1057\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1058\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1059\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1070\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1071\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1072\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1067\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1068\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1069\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1079\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1080\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1081\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1076\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1077\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1078\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1088\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1089\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1090\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1085\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1086\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1087\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1097\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1098\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1099\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1094\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1095\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1096\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1106\",\"attributes\":{\"data_source\":{\"type\":\"object\",\"name\":\"ColumnDataSource\",\"id\":\"p1100\",\"attributes\":{\"selected\":{\"type\":\"object\",\"name\":\"Selection\",\"id\":\"p1101\",\"attributes\":{\"indices\":[],\"line_indices\":[]}},\"selection_policy\":{\"type\":\"object\",\"name\":\"UnionRenderers\",\"id\":\"p1102\"},\"data\":{\"type\":\"map\",\"entries\":[[\"x\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"/Knx0k1iUD+sQAH/j2dSPx1M4YKcrFQ/23ENS1E5Vz+VHX+hhBZaP+m9zc4jTl0/6Ggnwap1YD9cbjs8UH1iP1youHILxWQ/kmI1ncNUZz/z1Ut3WTVqPwU77QTGcG0/s7+VkR6JcD9a0iMuKpNyP+6UzEKX3XQ/odY9X1Zwdz/2Xwq9UlR6P0uiWimRk30/yBFIX6mcgD/yoRvzHamCP/1kPRVA9oQ/UYd8twmMhz+aucqdcHOKP4Z3dWyFto0/ubhpRUuwkD/O+aepK7+SP5DAUwwGD5U/FHx0zN2nlz8qxc9Es5KaP7Zp1v6i2Z0/IipGXwTEoD9ACHJwU9WiPyDUgErpJ6U/9j/WxNLDpz+ShY/dGrKqP46WTxHq/K0/kx1JyNTXsD8COEdmleuyP3+AXvLpQLU/axeAx+jftz+8WrOTp9G6P0TO7NRaIL4/qrL+m7zrwD/xWhmq8QHDP/iKryYIWsU/YzZ++x/8xz80PhiTWfHKP+/X83r1Q84/apcT9rv/0D8o1f5aaBjTP5TNXwpEc9U/nvYKiHgY2D9NAM8HMRHbP9O15DS6Z94/oy5V8tIT4T8fyDKY+S7jP8RnhMCdjOU/Eg6PlPI06D8uhRweLjHrP0nqeTSpi+4/uLaxrAEo8T84PhWBpUXzP/7uW2wVpvU/zcWhSI5R+D+XAnoCUVH7P/+8qKvCr/4/UHA4QUg8AUA2Vis1bFwDQOyfTjGrvwVAETEJzEtuCECFPZXhmXELQC2AocwG1A5AbMUZzKZQEUA6bx/UTXMTQH2P7jJf2RVAc2S6RiuLGEBgyFDoCJIbQKbWz8l1+B5AonCnaR1lIUCdVMF9SoojQG3c95Qx8yVAmq3Z4CyoKEB7QcRDnrIrQLD52tUPHS9AXaRUNqx5MUBfagZSYqEzQO7gUHsiDTZA3cq6wlDFOEC0kTwhWtM7QIv/pSPVQT9AjjK2TlOOQUBN2QlxlbhDQJ1kCgoyJ0ZAjyPhFJfiSEA/KzyuPPRLQEYiUObFZk9AXLSCzxKjUUDwuwz7489TQKfOX2VgQVZAAAAAAAAAWUA=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"y\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"c+3rvg0aZLitfbEm2DUtuQ64CRTmjOC5sk1SLMg3gLoE/ek3NCIOuyso8jJczoy78jhJoFB6/rsiCf6FmgpjvKLaqS4Rxb28fk4020qnDr1SMkMk4cJVvW2p5TzzL5a9VtvvWiXY0L0QP7kOHqsDvhEu9Q33LDK+KKZvG/dFW749M43DwP6AvorFXIiK86G+uWOOZJxcwL5g9HNF+ybavj2dFY3DlvK+KJHOFJrNB7+hXkIW3cMbv2VD1b2KzC2/Ndiqt5qwPb+JT5yF7K1Lv3QdWGoLUVi/LFVFe6NBZL+2ob4w+xZwv8rCYE8vfni/JhPSoD/xgb+ZtJfxOWaJv7xhirI4bpG/iV7iPmREl78kJyL4zkmev7Um2tmrRaO/zha1gvkGqL+eXh6yP2etv24Gmuh9sLG/AE5983v1tL+0tJTuUXy4v+kvCgxAPby/JVzjYaQXwL+JHixWNiTCv4OjYR7zPsS/SaKBPTJixr9sS9ItIojIvw9px5Hkqsq/8P/BOKrEzL+Gd79Xz8/OvylUKCN8Y9C/3F3McZdS0b9Drzw/ADPSv65W39zNAtO/UEICt3bA07+5gimE3WrUv1oXdDJcAdW/z51vX8uD1b9zIbUNhfLVv2RGFV1hTta/rnQvVayY1r8aWvhgFNPWv5g8ieqR/9a/fa/Ehkkg179a9f8nazfXvwKO5KIRR9e/ASK9LSZR17/s59gwS1fXvz88iazQWte/xGyXzrJc17/P/9FyoV3Xv6Tif8ANXte/Ls2FWTpe17+er5nRSl7Xv2oZmjNQXte/83lsvFFe17/qYX8cUl7Xv4AfSjBSXte/k2inM1Je17+odx00Ul7Xv+gELDRSXte/DTsrNFJe179vSys0Ul7XvyNMKzRSXte/KEwrNFJe178oTCs0Ul7XvylMKzRSXte/J0wrNFJe1791Qys0Ul7XvylMKzRSXte/KEwrNFJe178nTCs0Ul7XvyhMKzRSXte/KEwrNFJe178oTCs0Ul7XvyhMKzRSXte/KUwrNFJe178oTCs0Ul7XvyhMKzRSXte/KUwrNFJe178=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}]]}}},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1107\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1108\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1103\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_width\":2,\"line_dash\":[6]}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1104\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_alpha\":0.1,\"line_width\":2,\"line_dash\":[6]}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1105\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_alpha\":0.2,\"line_width\":2,\"line_dash\":[6]}}}}],\"toolbar\":{\"type\":\"object\",\"name\":\"Toolbar\",\"id\":\"p1015\",\"attributes\":{\"tools\":[{\"type\":\"object\",\"name\":\"PanTool\",\"id\":\"p1030\"},{\"type\":\"object\",\"name\":\"WheelZoomTool\",\"id\":\"p1031\",\"attributes\":{\"renderers\":\"auto\"}},{\"type\":\"object\",\"name\":\"BoxZoomTool\",\"id\":\"p1032\",\"attributes\":{\"dimensions\":\"both\",\"overlay\":{\"type\":\"object\",\"name\":\"BoxAnnotation\",\"id\":\"p1033\",\"attributes\":{\"syncable\":false,\"line_color\":\"black\",\"line_alpha\":1.0,\"line_width\":2,\"line_dash\":[4,4],\"fill_color\":\"lightgrey\",\"fill_alpha\":0.5,\"level\":\"overlay\",\"visible\":false,\"left\":{\"type\":\"number\",\"value\":\"nan\"},\"right\":{\"type\":\"number\",\"value\":\"nan\"},\"top\":{\"type\":\"number\",\"value\":\"nan\"},\"bottom\":{\"type\":\"number\",\"value\":\"nan\"},\"left_units\":\"canvas\",\"right_units\":\"canvas\",\"top_units\":\"canvas\",\"bottom_units\":\"canvas\",\"handles\":{\"type\":\"object\",\"name\":\"BoxInteractionHandles\",\"id\":\"p1039\",\"attributes\":{\"all\":{\"type\":\"object\",\"name\":\"AreaVisuals\",\"id\":\"p1038\",\"attributes\":{\"fill_color\":\"white\",\"hover_fill_color\":\"lightgray\"}}}}}}}},{\"type\":\"object\",\"name\":\"SaveTool\",\"id\":\"p1040\"},{\"type\":\"object\",\"name\":\"ResetTool\",\"id\":\"p1041\"},{\"type\":\"object\",\"name\":\"HelpTool\",\"id\":\"p1042\"}]}},\"toolbar_location\":null,\"left\":[{\"type\":\"object\",\"name\":\"LinearAxis\",\"id\":\"p1025\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"FixedTicker\",\"id\":\"p1111\",\"attributes\":{\"ticks\":{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"AAAAAAAA+L/NzMzMzMz0v5qZmZmZmfG/zszMzMzM7L9oZmZmZmbmvwIAAAAAAOC/ODMzMzMz07+wmZmZmZm5vw==\"},\"shape\":[8],\"dtype\":\"float64\",\"order\":\"little\"},\"minor_ticks\":[]}},\"formatter\":{\"type\":\"object\",\"name\":\"BasicTickFormatter\",\"id\":\"p1027\"},\"axis_label\":\"head (m)\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1028\"}}}],\"below\":[{\"type\":\"object\",\"name\":\"LinearAxis\",\"id\":\"p1020\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"FixedTicker\",\"id\":\"p1110\",\"attributes\":{\"ticks\":{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"AAAAABQAAAAoAAAAPAAAAFAAAABkAAAA\"},\"shape\":[6],\"dtype\":\"int32\",\"order\":\"little\"},\"minor_ticks\":[]}},\"formatter\":{\"type\":\"object\",\"name\":\"BasicTickFormatter\",\"id\":\"p1022\"},\"axis_label\":\"time (days)\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1023\"}}}],\"center\":[{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1024\",\"attributes\":{\"axis\":{\"id\":\"p1020\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1029\",\"attributes\":{\"dimension\":1,\"axis\":{\"id\":\"p1025\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Legend\",\"id\":\"p1052\",\"attributes\":{\"location\":\"bottom_left\",\"label_text_font_size\":\"9pt\",\"glyph_height\":10,\"glyph_width\":10,\"spacing\":1,\"items\":[{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1053\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 0\"},\"renderers\":[{\"id\":\"p1049\"}]}},{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1063\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 1\"},\"renderers\":[{\"id\":\"p1060\"}]}},{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1109\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 1, Hantush\"},\"renderers\":[{\"id\":\"p1106\"}]}}]}}]}},{\"type\":\"object\",\"name\":\"Figure\",\"id\":\"p1112\",\"attributes\":{\"width\":380,\"height\":300,\"x_range\":{\"type\":\"object\",\"name\":\"DataRange1d\",\"id\":\"p1113\"},\"y_range\":{\"type\":\"object\",\"name\":\"DataRange1d\",\"id\":\"p1114\"},\"x_scale\":{\"type\":\"object\",\"name\":\"LogScale\",\"id\":\"p1122\"},\"y_scale\":{\"type\":\"object\",\"name\":\"LinearScale\",\"id\":\"p1123\"},\"title\":{\"type\":\"object\",\"name\":\"Title\",\"id\":\"p1115\",\"attributes\":{\"text\":\"(b)\"}},\"outline_line_color\":\"Black\",\"renderers\":[{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1153\",\"attributes\":{\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1154\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1155\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1150\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1151\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1152\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_1\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1164\",\"attributes\":{\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1165\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1166\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1161\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1162\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1163\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_1\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1174\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1175\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1176\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1171\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1172\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1173\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_2\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1183\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1184\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1185\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1180\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1181\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1182\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_2\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1192\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1193\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1194\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1189\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1190\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1191\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h0_3\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1201\",\"attributes\":{\"visible\":false,\"data_source\":{\"id\":\"p1003\"},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1202\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1203\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1198\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1199\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1200\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"t\"},\"y\":{\"type\":\"field\",\"field\":\"h1_3\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1210\",\"attributes\":{\"data_source\":{\"type\":\"object\",\"name\":\"ColumnDataSource\",\"id\":\"p1204\",\"attributes\":{\"selected\":{\"type\":\"object\",\"name\":\"Selection\",\"id\":\"p1205\",\"attributes\":{\"indices\":[],\"line_indices\":[]}},\"selection_policy\":{\"type\":\"object\",\"name\":\"UnionRenderers\",\"id\":\"p1206\"},\"data\":{\"type\":\"map\",\"entries\":[[\"x\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"/Knx0k1iUD+sQAH/j2dSPx1M4YKcrFQ/23ENS1E5Vz+VHX+hhBZaP+m9zc4jTl0/6Ggnwap1YD9cbjs8UH1iP1youHILxWQ/kmI1ncNUZz/z1Ut3WTVqPwU77QTGcG0/s7+VkR6JcD9a0iMuKpNyP+6UzEKX3XQ/odY9X1Zwdz/2Xwq9UlR6P0uiWimRk30/yBFIX6mcgD/yoRvzHamCP/1kPRVA9oQ/UYd8twmMhz+aucqdcHOKP4Z3dWyFto0/ubhpRUuwkD/O+aepK7+SP5DAUwwGD5U/FHx0zN2nlz8qxc9Es5KaP7Zp1v6i2Z0/IipGXwTEoD9ACHJwU9WiPyDUgErpJ6U/9j/WxNLDpz+ShY/dGrKqP46WTxHq/K0/kx1JyNTXsD8COEdmleuyP3+AXvLpQLU/axeAx+jftz+8WrOTp9G6P0TO7NRaIL4/qrL+m7zrwD/xWhmq8QHDP/iKryYIWsU/YzZ++x/8xz80PhiTWfHKP+/X83r1Q84/apcT9rv/0D8o1f5aaBjTP5TNXwpEc9U/nvYKiHgY2D9NAM8HMRHbP9O15DS6Z94/oy5V8tIT4T8fyDKY+S7jP8RnhMCdjOU/Eg6PlPI06D8uhRweLjHrP0nqeTSpi+4/uLaxrAEo8T84PhWBpUXzP/7uW2wVpvU/zcWhSI5R+D+XAnoCUVH7P/+8qKvCr/4/UHA4QUg8AUA2Vis1bFwDQOyfTjGrvwVAETEJzEtuCECFPZXhmXELQC2AocwG1A5AbMUZzKZQEUA6bx/UTXMTQH2P7jJf2RVAc2S6RiuLGEBgyFDoCJIbQKbWz8l1+B5AonCnaR1lIUCdVMF9SoojQG3c95Qx8yVAmq3Z4CyoKEB7QcRDnrIrQLD52tUPHS9AXaRUNqx5MUBfagZSYqEzQO7gUHsiDTZA3cq6wlDFOEC0kTwhWtM7QIv/pSPVQT9AjjK2TlOOQUBN2QlxlbhDQJ1kCgoyJ0ZAjyPhFJfiSEA/KzyuPPRLQEYiUObFZk9AXLSCzxKjUUDwuwz7489TQKfOX2VgQVZAAAAAAAAAWUA=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}],[\"y\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"c+3rvg0aZLitfbEm2DUtuQ64CRTmjOC5sk1SLMg3gLoE/ek3NCIOuyso8jJczoy78jhJoFB6/rsiCf6FmgpjvKLaqS4Rxb28fk4020qnDr1SMkMk4cJVvW2p5TzzL5a9VtvvWiXY0L0QP7kOHqsDvhEu9Q33LDK+KKZvG/dFW749M43DwP6AvorFXIiK86G+uWOOZJxcwL5g9HNF+ybavj2dFY3DlvK+KJHOFJrNB7+hXkIW3cMbv2VD1b2KzC2/Ndiqt5qwPb+JT5yF7K1Lv3QdWGoLUVi/LFVFe6NBZL+2ob4w+xZwv8rCYE8vfni/JhPSoD/xgb+ZtJfxOWaJv7xhirI4bpG/iV7iPmREl78kJyL4zkmev7Um2tmrRaO/zha1gvkGqL+eXh6yP2etv24Gmuh9sLG/AE5983v1tL+0tJTuUXy4v+kvCgxAPby/JVzjYaQXwL+JHixWNiTCv4OjYR7zPsS/SaKBPTJixr9sS9ItIojIvw9px5Hkqsq/8P/BOKrEzL+Gd79Xz8/OvylUKCN8Y9C/3F3McZdS0b9Drzw/ADPSv65W39zNAtO/UEICt3bA07+5gimE3WrUv1oXdDJcAdW/z51vX8uD1b9zIbUNhfLVv2RGFV1hTta/rnQvVayY1r8aWvhgFNPWv5g8ieqR/9a/fa/Ehkkg179a9f8nazfXvwKO5KIRR9e/ASK9LSZR17/s59gwS1fXvz88iazQWte/xGyXzrJc17/P/9FyoV3Xv6Tif8ANXte/Ls2FWTpe17+er5nRSl7Xv2oZmjNQXte/83lsvFFe17/qYX8cUl7Xv4AfSjBSXte/k2inM1Je17+odx00Ul7Xv+gELDRSXte/DTsrNFJe179vSys0Ul7XvyNMKzRSXte/KEwrNFJe178oTCs0Ul7XvylMKzRSXte/J0wrNFJe1791Qys0Ul7XvylMKzRSXte/KEwrNFJe178nTCs0Ul7XvyhMKzRSXte/KEwrNFJe178oTCs0Ul7XvyhMKzRSXte/KUwrNFJe178oTCs0Ul7XvyhMKzRSXte/KUwrNFJe178=\"},\"shape\":[100],\"dtype\":\"float64\",\"order\":\"little\"}]]}}},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1211\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1212\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1207\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_width\":2,\"line_dash\":[6]}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1208\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_alpha\":0.1,\"line_width\":2,\"line_dash\":[6]}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1209\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#2ca02c\",\"line_alpha\":0.2,\"line_width\":2,\"line_dash\":[6]}}}}],\"toolbar\":{\"type\":\"object\",\"name\":\"Toolbar\",\"id\":\"p1121\",\"attributes\":{\"tools\":[{\"type\":\"object\",\"name\":\"PanTool\",\"id\":\"p1134\"},{\"type\":\"object\",\"name\":\"WheelZoomTool\",\"id\":\"p1135\",\"attributes\":{\"renderers\":\"auto\"}},{\"type\":\"object\",\"name\":\"BoxZoomTool\",\"id\":\"p1136\",\"attributes\":{\"dimensions\":\"both\",\"overlay\":{\"type\":\"object\",\"name\":\"BoxAnnotation\",\"id\":\"p1137\",\"attributes\":{\"syncable\":false,\"line_color\":\"black\",\"line_alpha\":1.0,\"line_width\":2,\"line_dash\":[4,4],\"fill_color\":\"lightgrey\",\"fill_alpha\":0.5,\"level\":\"overlay\",\"visible\":false,\"left\":{\"type\":\"number\",\"value\":\"nan\"},\"right\":{\"type\":\"number\",\"value\":\"nan\"},\"top\":{\"type\":\"number\",\"value\":\"nan\"},\"bottom\":{\"type\":\"number\",\"value\":\"nan\"},\"left_units\":\"canvas\",\"right_units\":\"canvas\",\"top_units\":\"canvas\",\"bottom_units\":\"canvas\",\"handles\":{\"type\":\"object\",\"name\":\"BoxInteractionHandles\",\"id\":\"p1143\",\"attributes\":{\"all\":{\"type\":\"object\",\"name\":\"AreaVisuals\",\"id\":\"p1142\",\"attributes\":{\"fill_color\":\"white\",\"hover_fill_color\":\"lightgray\"}}}}}}}},{\"type\":\"object\",\"name\":\"SaveTool\",\"id\":\"p1144\"},{\"type\":\"object\",\"name\":\"ResetTool\",\"id\":\"p1145\"},{\"type\":\"object\",\"name\":\"HelpTool\",\"id\":\"p1146\"}]}},\"toolbar_location\":null,\"left\":[{\"type\":\"object\",\"name\":\"LinearAxis\",\"id\":\"p1129\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"FixedTicker\",\"id\":\"p1214\",\"attributes\":{\"ticks\":{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"mpmZmZmZ+b9nZmZmZmb2vzQzMzMzM/O/AQAAAAAA8L+cmZmZmZnpvzYzMzMzM+O/oJmZmZmZ2b+omZmZmZnJvwAAAAAAAMC8\"},\"shape\":[9],\"dtype\":\"float64\",\"order\":\"little\"},\"minor_ticks\":[]}},\"formatter\":{\"type\":\"object\",\"name\":\"BasicTickFormatter\",\"id\":\"p1131\"},\"axis_label\":\"head (m)\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1132\"}}}],\"below\":[{\"type\":\"object\",\"name\":\"LogAxis\",\"id\":\"p1124\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"LogTicker\",\"id\":\"p1125\",\"attributes\":{\"num_minor_ticks\":10,\"mantissas\":[1,5]}},\"formatter\":{\"type\":\"object\",\"name\":\"LogTickFormatter\",\"id\":\"p1126\"},\"axis_label\":\"time (days)\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1127\"}}}],\"center\":[{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1128\",\"attributes\":{\"axis\":{\"id\":\"p1124\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1133\",\"attributes\":{\"dimension\":1,\"axis\":{\"id\":\"p1129\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Legend\",\"id\":\"p1156\",\"attributes\":{\"location\":\"bottom_left\",\"label_text_font_size\":\"9pt\",\"glyph_height\":10,\"glyph_width\":10,\"spacing\":1,\"items\":[{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1157\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 0\"},\"renderers\":[{\"id\":\"p1153\"}]}},{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1167\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 1\"},\"renderers\":[{\"id\":\"p1164\"}]}},{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1213\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Aquifer 1, Hantush\"},\"renderers\":[{\"id\":\"p1210\"}]}}]}}]}}]}},{\"type\":\"object\",\"name\":\"Row\",\"id\":\"p1218\",\"attributes\":{\"align\":\"center\",\"children\":[{\"type\":\"object\",\"name\":\"RadioButtonGroup\",\"id\":\"p1215\",\"attributes\":{\"js_property_callbacks\":{\"type\":\"map\",\"entries\":[[\"change:active\",[{\"type\":\"object\",\"name\":\"CustomJS\",\"id\":\"p1216\",\"attributes\":{\"args\":{\"type\":\"map\",\"entries\":[[\"source\",{\"id\":\"p1003\"}],[\"radio_button_group\",{\"id\":\"p1215\"}],[\"h0_11\",{\"id\":\"p1049\"}],[\"h1_11\",{\"id\":\"p1060\"}],[\"h0_21\",{\"id\":\"p1070\"}],[\"h1_21\",{\"id\":\"p1079\"}],[\"h0_31\",{\"id\":\"p1088\"}],[\"h1_31\",{\"id\":\"p1097\"}],[\"h0_12\",{\"id\":\"p1153\"}],[\"h1_12\",{\"id\":\"p1164\"}],[\"h0_22\",{\"id\":\"p1174\"}],[\"h1_22\",{\"id\":\"p1183\"}],[\"h0_32\",{\"id\":\"p1192\"}],[\"h1_32\",{\"id\":\"p1201\"}]]},\"code\":\"\\n // Hide all lines\\n h0_11.visible = false;\\n h1_11.visible = false;\\n h0_21.visible = false;\\n h1_21.visible = false;\\n h0_31.visible = false;\\n h1_31.visible = false;\\n h0_12.visible = false;\\n h1_12.visible = false;\\n h0_22.visible = false;\\n h1_22.visible = false;\\n h0_32.visible = false;\\n h1_32.visible = false;\\n \\n var active = radio_button_group.active;\\n if (active == 0) {\\n h0_11.visible = true;\\n h1_11.visible = true;\\n h0_12.visible = true;\\n h1_12.visible = true;\\n } else if (active == 1) {\\n h0_21.visible = true;\\n h1_21.visible = true;\\n h0_22.visible = true;\\n h1_22.visible = true;\\n } else if (active == 2) {\\n h0_31.visible = true;\\n h1_31.visible = true;\\n h0_32.visible = true;\\n h1_32.visible = true;\\n }\\n \"}}]]]},\"labels\":[\"S\\u2080 = 0.1\",\"S\\u2080 = 0.01\",\"S\\u2080 = 0.001\"],\"active\":0}}]}}]}}]}};\n const render_items = [{\"docid\":\"00010d70-bd92-4dfe-938d-35658ff947ac\",\"roots\":{\"p1219\":\"cbcfc790-27d9-4028-97f8-5a549c989e40\"},\"root_ids\":[\"p1219\"]}];\n void root.Bokeh.embed.embed_items_notebook(docs_json, render_items);\n }\n if (root.Bokeh !== undefined) {\n embed_document(root);\n } else {\n let attempts = 0;\n const timer = setInterval(function(root) {\n if (root.Bokeh !== undefined) {\n clearInterval(timer);\n embed_document(root);\n } else {\n attempts++;\n if (attempts > 100) {\n clearInterval(timer);\n console.log(\"Bokeh: ERROR: Unable to run BokehJS code because BokehJS library is missing\");\n }\n }\n }, 10, root)\n }\n})(window);", "application/vnd.bokehjs_exec.v0+json": "" }, "metadata": { "application/vnd.bokehjs_exec.v0+json": { "id": "p1219" } }, "output_type": "display_data" } ], "source": [ "inter_fig_chap9.fig_9_10()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{raw} html\n", "\n", "```\n", "\n", "```{figure} ../py_figures/1x1_pixel.png\n", "---\n", "name: fig_two_aquifers9.5\n", "---\n", "Drawdown in an observation well near a pumping well in the bottom aquifer of a two-aquifer system. The top aquifer is unconfined; linear time axis (left) and log-scaled time axis (right).\n", "```" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "## Test using ttim\n", "# from ttim import *\n", "# ml = ModelMaq(kaq=[10, 20], z=[0, -10, -15, -25], c=[c], Saq=[S0 / 10, S1 / 10], tmin=1e-3, tmax=1e3)\n", "# w = Well(ml, 0, 0, rw=0.3, tsandQ=[(0, Q)], layers=1)\n", "# ml.solve()\n", "# httim = ml.head(250, 0, t)\n", "# plt.semilogx(t, httim[0])\n", "# plt.semilogx(t, httim[1])\n", "# plt.semilogx(t, h0, '.', color='C0')\n", "# plt.semilogx(t, h1, '.', color='C1')" ] } ], "metadata": { "celltoolbar": "Tags", "kernelspec": { "display_name": "jupbook", "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.2" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false } }, "nbformat": 4, "nbformat_minor": 4 }