{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "tags": [ "remove-input" ] }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\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\n", "import matplotlib as mpl\n", "from matplotlib.lines import Line2D\n", "mpl.rcParams['contour.negative_linestyle'] = 'solid'\n", "from myst_nb import glue\n", "\n", "import inter_fig_chap10" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "(content:chap10:interface_flow_strip_island)=\n", "# Interface flow below a strip island" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "Two-dimensional steady interface flow in the vertical plane is explored further by considering interface flow in a cross-section of a strip island; the saltwater is at rest. The distance between the two shores (i.e., the width of the island) is $2L$, and the infiltration rate is equal to $N$ ({numref}`Figure %s `). Flow is unconfined, but the small variation in the phreatic surface is neglected and the top of the aquifer is approximated as horizontal at the same elevation as the sea bottom. The origin of the coordinate system is chosen at the top of the aquifer in the center of the island. The sea bottom is horizontal and the head along the sea bottom is set equal to zero. \n", "\n", "$$\n", "h\\vert_{x>L,z=0} = h\\vert_{x<-L,z=0} = 0\n", "$$(eq_hbalance10.6)\n", "\n", "The stream function along the interface is set equal to zero. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n", "application/papermill.record/text/plain": "

" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_flowcrosection" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = diagrams_2d.interface_island()\n", "glue('fig_flowcrosection', fig, display=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{glue:figure} fig_flowcrosection\n", "---\n", "name: fig_flowcrosection10.6\n", "---\n", "Two-dimensional interface flow in a vertical cross-section of a strip island.\n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " The complex specific discharge potential solution to the stated problem is obtained from {cite:author}`bruggeman1999analytical` ({cite:year}`bruggeman1999analytical`, Eq. 910.03), modified for the coordinate system chosen here. The solution has the following form\n", "\n", "$$\n", "\\omega^2 +2i N\\zeta\\omega + kN \\zeta^2/\\alpha - N(N + k/\\alpha)L^2 = 0\n", "$$(eq_omegavdveer10.6)\n", "\n", "Note that Eq. ({eq}`eq_omegavdveer10.6`) is a second-order polynomial in $\\omega$. A solution for $\\omega$ for a given value of $\\zeta$ is obtained using the standard formula for the root of a parabola, which gives\n", "\n", "$$\n", "\\omega = -iN\\zeta + \\sqrt{N^2(L^2 - \\zeta^2) + Nk(L^2 - \\zeta^2)/\\alpha}\n", "$$(eq_omega-iN10.6)\n", "\n", "The specific discharge potential along the interface equals $\\phi_i=-kz_i\\alpha$, so that an equation for the elevation of the interface as a function of $x$ may be obtained by substituting $\\omega=-kz_i\\alpha$ for $\\omega$ and $x+iz_i$ for $\\zeta$ in Eq. ({eq}`eq_omegavdveer10.6`), which gives (after considerable algebra)\n", "\n", "$$\n", "z_i = -\\sqrt{ \\frac{N\\alpha}{k}L^2 - \\frac{N\\alpha}{N\\alpha+k} x^2}\n", "$$(eq_vdveerinterface10.6)\n", "\n", "The shape of the interface is an ellipse and the locations of the shorelines are the foci of the ellipse. \n", "The location $x_\\text{tip}$ of the tip of the interface is obtained by setting $z_i=0$ in Eq. ({eq}`eq_vdveerinterface10.6`), which gives\n", "\n", "$$\n", "x_\\text{tip} = \\pm L\\sqrt{1 + \\frac{N\\alpha}{k}}\n", "$$(eq_vdveertip10.6)\n", "\n", "The elevation of the interface at the center of the island is obtained by setting $x=0$ in Eq. ({eq}`eq_vdveerinterface10.6`), which gives\n", "\n", "$$\n", "z_i|_{x=0} = -L\\sqrt{\\frac{N\\alpha}{k}}\n", "$$(eq_ziatzero10.6)\n", "\n", "A flow net is constructed as follows. First, the location of the tip is computed with Eq. ({eq}`eq_vdveertip10.6`). Second, an array of $x$ values is created from $-x_\\text{tip} + \\varepsilon$ to $-x_\\text{tip} - \\varepsilon$, where $\\varepsilon$ is a small number ($\\varepsilon$ is used to make sure that the first point is numerically not to the left of $-x_\\text{tip}$ and the last point is numerically not to the right of $+x_\\text{tip}$). Third, a grid of $z$ values is computed between the interface and the top of the aquifer for every $x$-value. Fourth and final, $\\omega$ is computed on the grid and a flow net is contoured. \n", "
\n", "In the example below, the width of the island equals 2000 m; the flow net is shown in {numref}`Figure %s `. The freshwater zone is a bit more than 63 m thick at the center of the island, while the outflow zone is only 2 m on each side. Note that, contrary to other flow nets in this chapter, the vertical scale is exaggerated in {numref}`Figure %s ` (otherwise the figure would be unreadable), which makes the flow net resemble flow in an anisotropic aquifer while this is not the case. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# parameters\n", "k = 10 # hydraulic conductivity, m/d\n", "L = 1000 # half width of island, m\n", "N = 0.001 # areal recharge, m/d\n", "rhof = 1000 # density freshwater, kg/m^3\n", "rhos = 1025 # density of saltwater, kg/m^3\n", "alpha = rhof / (rhos - rhof)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "length outflow zone: 2.00 m\n", "elevation of interface at center of island: -63.25 m\n" ] } ], "source": [ "# solution\n", "def omega(zeta, k=k, alpha=alpha, L=L, N=N):\n", " om = -1j * N * zeta + np.sqrt(N ** 2 * (L ** 2 - zeta ** 2) + \n", " N * k * (L ** 2 - zeta ** 2) / alpha)\n", " return om\n", "\n", "# compute interface\n", "ng = 201\n", "xtip = L * np.sqrt(1 + N * alpha / k)\n", "print(f'length outflow zone: {xtip - L:0.2f} m')\n", "xi = np.linspace(-xtip + 1e-12, xtip - 1e-12, ng)\n", "zi = -np.sqrt(N * alpha / k * L ** 2 - N * alpha / (N * alpha + k) * xi ** 2)\n", "zcenter = -L * np.sqrt(N * alpha / k)\n", "print(f'elevation of interface at center of island: {zcenter:0.2f} m')\n", "# compute grid for contouring\n", "xg = xi * np.ones((50, ng))\n", "zg = np.zeros_like(xg)\n", "for i in range(ng):\n", " zg[:, i] = np.linspace(zi[i], -1e-12, 50)\n", "om = omega(xg + 1j * zg)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic plot\n", "plt.subplot(111, aspect=5, xlim=(-xtip, xtip))\n", "plt.contour(xg, zg, om.imag, 20)\n", "plt.contour(xg, zg, om.real, 20)\n", "plt.plot(xi, zi, 'k');" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "editable": true, "scrolled": true, "slideshow": { "slide_type": "" }, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n", "application/papermill.record/text/plain": "

" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_notevertexagg" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "fig = plt.figure(figsize=diagrams_2d.figdef2)\n", "plt.subplot(111, aspect=5)\n", "plt.contour(xg, zg, om.imag, 20, colors='C0')\n", "plt.contour(xg, zg, om.real, 20, colors='C1')\n", "plt.plot(xi, zi, color=diagrams_2d.C1dark)\n", "plt.xlabel('x (m)')\n", "plt.ylabel('z (m) - VE:5x')\n", "plt.tight_layout();\n", "glue('fig_notevertexagg', fig, display=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{glue:figure} fig_notevertexagg\n", "---\n", "name: fig_notevertexagg10.6\n", "---\n", "Flow net for two-dimensional interface flow in a cross-section of a strip island with isotropic hydraulic conductivity (note vertical exaggeration VE=5$\\times$). \n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "      The solution for an anisotropic hydraulic conductivity is obtained using the approach of {numref}`content:chap10_vertical_anisotropy`. In addition, the infiltration rate $\\tilde{N}$ in the transformed domain is\n", "\n", "$$\n", "\\tilde{N}=N\\sqrt{k_x/k_z}\n", "$$(eq_ntransdom10.6)\n", "\n", "The resulting equation for the elevation of the interface is \n", "\n", "$$\n", "z_i = -\\sqrt{ \\frac{N\\alpha}{k_x}L^2 - \\frac{N\\alpha}{N\\alpha+k_z} \\frac{k_z}{k_x}x^2}\n", "$$(eq_zelevation10.6)\n", "\n", "Hence, the elevation of the interface at the center of the island is\n", "\n", "$$\n", "z_i|_{x=0} = -\\sqrt{ \\frac{N\\alpha}{k_x}L^2}\n", "$$(eq_zi_center_2d10.6)\n", "\n", "This result may be surprising, as it means that the depth of the interface at the center of the island only depends on the horizontal component $k_x$ of the hydraulic conductivity and not on the vertical component. Conversely, the equation for the tip of the interface only depends on $k_z$ and not on $k_x$\n", "\n", "$$\n", "x_\\text{tip} = L\\sqrt{1 + \\frac{N\\alpha}{k_z}}\n", "$$(eq_xtipsq10.6)\n", "\n", "This means that the length of the outflow zone is larger when the vertical hydraulic conductivity is smaller. \n", "
     \n", "The solution for two-dimensional flow in the vertical plane is compared to the Dupuit solution. The Dupuit solution for the interface is obtained from {numref}`content:chap4:unconfined_flow` as\n", "\n", "$$\n", "z_{i} = -\\sqrt{\\frac{N\\alpha}{k_x}(L^2-x^2)}\n", "$$(eq_zfracalpha10.6)\n", "\n", "so that the Dupuit elevation of the interface at the center of the island is\n", "\n", "$$\n", "z_{i}|_{x=0} = -\\sqrt{ \\frac{N\\alpha}{k_x}L^2}\n", "$$(eq_dupelzi10.6)\n", "\n", "which is exactly the same as the two-dimensional solution (Eq. {eq}`eq_zi_center_2d10.6`) irrespective of the vertical hydraulic conductivity! Everywhere else, the Dupuit solution is conservative: the thickness of the freshwater zone is smaller than in the two-dimensional solution. \n", "
     \n", "In the example below, the interface obtained with the Dupuit solution is compared to the two-dimensional solution for the case that the hydraulic conductivity is isotropic and for a case that the hydraulic conductivity is anisotropic. The Dupuit solution is again a very good and conservative approximation ({numref}`Figure %s `)." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# additional parameters for anisotropic case\n", "kx = k\n", "kz = kx / 20" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# solution\n", "xtipaniso = L * np.sqrt(1 + N * alpha / kz)\n", "xaniso = np.linspace(-xtipaniso + 1e-12, xtipaniso - 1e-12, 400)\n", "zianiso = -np.sqrt(N * alpha / kx * L ** 2 - \n", " N * alpha / (N * alpha + kz) * kz / kx * xaniso ** 2)\n", "xdupuit = np.linspace(-L, L, 200)\n", "zidupuit = -np.sqrt(N * alpha / kx * (L ** 2 - xdupuit ** 2))" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "tags": [ "remove-output" ] }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# basic plot\n", "plt.subplot(111, aspect=5)\n", "plt.plot(xi, zi, label='2D, $k_z=k_x$')\n", "plt.plot(xdupuit, zidupuit, '--', label='Dupuit')\n", "plt.plot(xaniso, zianiso, label='2D, $k_z=k_x/20$')\n", "plt.legend();" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "scrolled": true, "tags": [ "remove-input", "remove-output" ] }, "outputs": [ { "data": { "application/papermill.record/image/svg+xml": "\n\n\n", "application/papermill.record/text/plain": "

" }, "metadata": { "scrapbook": { "mime_prefix": "application/papermill.record/", "name": "fig_compinterface" } }, "output_type": "display_data" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=diagrams_2d.figdef2)\n", "plt.subplot(111, aspect=5)\n", "plt.plot(xi, zi, color=diagrams_2d.C1dark, label='2D, $k_z=k_x$')\n", "plt.plot(xaniso, zianiso, 'C2', label='2D, $k_z=k_x/20$')\n", "plt.plot(xdupuit, zidupuit, 'C1--', label='Dupuit')\n", "plt.legend()\n", "plt.grid()\n", "plt.xlabel('x (m)')\n", "plt.ylabel('z (m) - VE:5x')\n", "plt.tight_layout();\n", "glue(\"fig_compinterface\", fig, display=False)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "editable": true, "scrolled": true, "slideshow": { "slide_type": "" }, "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", "\n", "const 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", " \"

re-rerun `output_notebook()` to attempt to load from CDN again, or
use INLINE resources instead, as so:

\\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", " },\n", "function(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 = {\"1ede088a-211a-4a24-96cc-cbb881efb049\":{\"version\":\"3.7.3\",\"title\":\"Bokeh Application\",\"roots\":[{\"type\":\"object\",\"name\":\"Column\",\"id\":\"p1066\",\"attributes\":{\"children\":[{\"type\":\"object\",\"name\":\"Div\",\"id\":\"p1065\",\"attributes\":{\"styles\":{\"type\":\"map\",\"entries\":[[\"margin\",\"0 auto\"],[\"width\",\"750px\"]]}}},{\"type\":\"object\",\"name\":\"Figure\",\"id\":\"p1006\",\"attributes\":{\"width\":750,\"height\":150,\"x_range\":{\"type\":\"object\",\"name\":\"DataRange1d\",\"id\":\"p1007\",\"attributes\":{\"start\":-1200,\"end\":1200}},\"y_range\":{\"type\":\"object\",\"name\":\"DataRange1d\",\"id\":\"p1008\",\"attributes\":{\"start\":-70,\"end\":5}},\"x_scale\":{\"type\":\"object\",\"name\":\"LinearScale\",\"id\":\"p1015\"},\"y_scale\":{\"type\":\"object\",\"name\":\"LinearScale\",\"id\":\"p1016\"},\"title\":{\"type\":\"object\",\"name\":\"Title\",\"id\":\"p1013\"},\"renderers\":[{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1046\",\"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\":[[\"x_values\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"ihOE6ftPj8D8gzj+0v+OwG707BKqr47A32ShJ4FfjsBR1VU8WA+OwMNFClEvv43ANba+ZQZvjcCnJnN63R6NwBiXJ4+0zozAigfco4t+jMD8d5C4Yi6MwG7oRM053ovA4Fj54RCOi8BSya325z2LwMQ5Ygu/7YrANaoWIJadisCnGss0bU2KwBmLf0lE/YnAivszXhuticD8a+hy8lyJwG7cnIfJDInA4ExRnKC8iMBSvQWxd2yIwMQtusVOHIjANp5u2iXMh8CnDiPv/HuHwBl/1wPUK4fAi++LGKvbhsD9X0AtgouGwG7Q9EFZO4bA4ECpVjDrhcBSsV1rB5uFwMQhEoDeSoXANpLGlLX6hMCoAnupjKqEwBpzL75jWoTAi+Pj0joKhMD9U5jnEbqDwG/ETPzoaYPA4DQBEcAZg8BSpbUll8mCwMQVajpueYLANoYeT0UpgsCo9tJjHNmBwBpnh3jziIHAjNc7jco4gcD+R/ChoeiAwG+4pLZ4mIDA4ShZy09IgMCmMhvATfB/wIkThOn7T3/AbfTsEqqvfsBQ1VU8WA9+wDS2vmUGb33AGJcnj7TOfMD8d5C4Yi58wOBY+eEQjnvAwjliC7/tesCmGss0bU16wIr7M14brXnAbtych8kMecBSvQWxd2x4wDSebtolzHfAGH/XA9Qrd8D8X0Atgot2wOBAqVYw63XAxCESgN5KdcCmAnupjKp0wIrj49I6CnTAbsRM/Ohpc8BSpbUll8lywDaGHk9FKXLAGGeHePOIccD8R/ChoehwwOAoWctPSHDAiBOE6ftPb8BQ1VU8WA9uwBiXJ4+0zmzA3Fj54RCOa8CkGss0bU1qwGzcnIfJDGnANJ5u2iXMZ8D8X0AtgotmwMAhEoDeSmXAiOPj0joKZMBQpbUll8liwBhnh3jziGHA4ChZy09IYMBQ1VU8WA9ewNhY+eEQjlvAaNych8kMWcD4X0AtgotWwIjj49I6ClTAGGeHePOIUcBA1VU8WA9OwGDcnIfJDEnAgOPj0joKRMBA1VU8WA8+wIDj49I6CjTAQOPj0joKJMAAAAAAAABAPcDj49I6CiRAoOPj0joKNECA1VU8WA8+QKDj49I6CkRAgNych8kMSUBg1VU8WA9OQCBnh3jziFFAkOPj0joKVEAAYEAtgotWQHDcnIfJDFlA4Fj54RCOW0BQ1VU8WA9eQOgoWctPSGBAIGeHePOIYUBYpbUll8liQJDj49I6CmRAyCESgN5KZUAAYEAtgotmQDiebtolzGdAcNych8kMaUCoGss0bU1qQOBY+eEQjmtAIJcnj7TObEBY1VU8WA9uQJAThOn7T29A5ChZy09IcEAASPChoehwQBxnh3jziHFAOIYeT0UpckBUpbUll8lyQHDETPzoaXNAjOPj0joKdECoAnupjKp0QMghEoDeSnVA5ECpVjDrdUAAYEAtgot2QBx/1wPUK3dAOJ5u2iXMd0BUvQWxd2x4QHDcnIfJDHlAjPszXhuteUCoGss0bU16QMQ5Ygu/7XpA5Fj54RCOe0AAeJC4Yi58QByXJ4+0znxAOLa+ZQZvfUBU1VU8WA9+QHD07BKqr35AjBOE6ftPf0CoMhvATfB/QOIoWctPSIBAcLiktniYgED+R/ChoeiAQI7XO43KOIFAHGeHePOIgUCq9tJjHNmBQDiGHk9FKYJAxhVqOm55gkBUpbUll8mCQOI0ARHAGYNAcMRM/Ohpg0D+U5jnEbqDQIzj49I6CoRAGnMvvmNahECqAnupjKqEQDiSxpS1+oRAxiESgN5KhUBUsV1rB5uFQOJAqVYw64VAcND0QVk7hkD+X0AtgouGQIzvixir24ZAGn/XA9Qrh0CoDiPv/HuHQDaebtolzIdAxi26xU4ciEBUvQWxd2yIQOJMUZygvIhAcNych8kMiUD+a+hy8lyJQIz7M14brYlAGot/SUT9iUCoGss0bU2KQDaqFiCWnYpAxDliC7/tikBUya325z2LQOJY+eEQjotAcOhEzTnei0D+d5C4Yi6MQIwH3KOLfoxAGpcnj7TOjECoJnN63R6NQDa2vmUGb41AxEUKUS+/jUBS1VU8WA+OQOBkoSeBX45AcPTsEqqvjkD+gzj+0v+OQIoThOn7T49A\"},\"shape\":[201],\"dtype\":\"float64\",\"order\":\"little\"}],[\"z_values\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"WQb3JFyoyL6xKnUTAdghwB+DH67hKynA9Yjo9SfALsAtQP/1cLUxwK5stUyYvzPA+8WYfuiTNcBLhYn1Gj83wHRyP/R+yTjAfaFxYt04OsARZYnWbpE7wK0cVc1k1jzAzNXXsTsKPsAPYTWk7i4/wCbpKM8MI0DAUjWQZYioQMCHfW98eChBwDPber5Vo0HA9n35focZQsD8B5IJZ4tCwInEIyxC+ULAPpFyMF1jQ8D8dvRq9MlDwP/2J3k9LUTAE893QmiNRMDa25fJn+pEwH6MsdgKRUXA/1InkMycRcCTxdvdBPJFwDd3jOHQREbAmkLNQUuVRsCiu210jONGwN08fvyqL0fApTu4n7t5R8A2aLWU0cFHwHpCHKv+B0jAMJKyblNMSMD8KxtG345IwAwO4o2wz0jA8zpesNQOScBYW9o6WExJwFn6YvBGiEnAXEmK2qvCScAoFWVYkftJwGxy+isBM0rA8jpXhgRpSsDWbG8SpJ1KwECY8f7n0ErAHpIrB9gCS8B8cRt7ezNLwPZOxEbZYkvAHTfb+PeQS8BIK97I3b1LwCHVpJyQ6UvAMqp4DRYUTMBemMBscz1MwAXpSsitZUzAYc0+7smMTMDg9L1wzLJMwFSfPam510zADsydu5X7TMBfcASZZB5NwAj/gQIqQE3AI/2Di+lgTcAr5BmcpoBNwLQiEHNkn03ApKvkJya9TcA1LZis7tlNwN67Xs/A9U3AlHgyPJ8QTsCHfEp+jCpOwJQdeAGLQ07Al21sE51bTsCDqufkxHJOwIgt1IoEiU7A2UNO/12eTsAoO5oi07JOwP/OCbxlxk7AfAnSehfZTsAKktL26epOwNpOT7He+07AByqdFfcLT8AgucJ5NBtPwER1DR+YKU/A0iKcMiM3T8CZ+d7N1kNPwLwRDvezT0/AhY2WobtaT8CB7X6u7mRPwOnxwuxNbk/ADWOnGdp2T8B1EQbhk35PwDhWkt17hU/AJlQWmZKLT8AhM6mM2JBPwOeH3iBOlU/A8hPvrfOYT8DAAtt7yZtPwJvDhcLPnU/AKZnLqQafT8An8pBJbp9PwCmZy6kGn0/Am8OFws+dT8DAAtt7yZtPwPIT763zmE/A54feIE6VT8AhM6mM2JBPwCZUFpmSi0/AOFaS3XuFT8B1EQbhk35PwA1jpxnadk/A6PHC7E1uT8CB7X6u7mRPwIWNlqG7Wk/AuxEO97NPT8CZ+d7N1kNPwNEinDIjN0/AQ3UNH5gpT8AfucJ5NBtPwAcqnRX3C0/A2k5Psd77TsAKktL26epOwHwJ0noX2U7A/s4JvGXGTsAnO5oi07JOwNlDTv9dnk7AiC3UigSJTsCCqufkxHJOwJdtbBOdW07AlB14AYtDTsCGfEp+jCpOwJR4MjyfEE7A3rtez8D1TcA1LZis7tlNwKOr5CcmvU3AtCIQc2SfTcAq5BmcpoBNwCL9g4vpYE3AB/+BAipATcBecASZZB5NwA3MnbuV+0zAU589qbnXTMDg9L1wzLJMwGDNPu7JjEzABOlKyK1lTMBdmMBscz1MwDGqeA0WFEzAINWknJDpS8BHK97I3b1LwBw32/j3kEvA9U7ERtliS8B8cRt7ezNLwB2SKwfYAkvAQJjx/ufQSsDWbG8SpJ1KwPE6V4YEaUrAa3L6KwEzSsAnFWVYkftJwFtJitqrwknAWPpi8EaIScBWW9o6WExJwPE6XrDUDknADA7ijbDPSMD8KxtG345IwC+Ssm5TTEjAekIcq/4HSMA0aLWU0cFHwKI7uJ+7eUfA3Dx+/KovR8Cgu210jONGwJlCzUFLlUbANXeM4dBERsCSxdvdBPJFwP5SJ5DMnEXAfYyx2ApFRcDZ25fJn+pEwBPPd0JojUTA/PYneT0tRMD6dvRq9MlDwDuRcjBdY0PAhsQjLEL5QsD6B5IJZ4tCwPN9+X6HGULAMdt6vlWjQcCFfW98eChBwFA1kGWIqEDAJukozwwjQMAIYTWk7i4/wMXV17E7Cj7AphxVzWTWPMAIZYnWbpE7wHOhcWLdODrAanI/9H7JOMBFhYn1Gj83wPXFmH7okzXAq2y1TJi/M8AmQP/1cLUxwOWI6PUnwC7AAYMfruErKcB4KnUTAdghwFkG9yRcqMi+\"},\"shape\":[201],\"dtype\":\"float64\",\"order\":\"little\"}]]}}},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1047\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1048\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1043\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x_values\"},\"y\":{\"type\":\"field\",\"field\":\"z_values\"},\"line_color\":\"#1f77b4\",\"line_width\":2}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1044\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x_values\"},\"y\":{\"type\":\"field\",\"field\":\"z_values\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.1,\"line_width\":2}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1045\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x_values\"},\"y\":{\"type\":\"field\",\"field\":\"z_values\"},\"line_color\":\"#1f77b4\",\"line_alpha\":0.2,\"line_width\":2}}}},{\"type\":\"object\",\"name\":\"GlyphRenderer\",\"id\":\"p1055\",\"attributes\":{\"data_source\":{\"type\":\"object\",\"name\":\"ColumnDataSource\",\"id\":\"p1049\",\"attributes\":{\"selected\":{\"type\":\"object\",\"name\":\"Selection\",\"id\":\"p1050\",\"attributes\":{\"indices\":[],\"line_indices\":[]}},\"selection_policy\":{\"type\":\"object\",\"name\":\"UnionRenderers\",\"id\":\"p1051\"},\"data\":{\"type\":\"map\",\"entries\":[[\"x\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"AAAAAABAj8CCjd4Vme+OwAQbvSsyn47AhqibQctOjsAINnpXZP6NwIrDWG39rY3ADFE3g5ZdjcCO3hWZLw2NwA9s9K7IvIzAkfnSxGFsjMATh7Ha+huMwJUUkPCTy4vAF6JuBi17i8CZL00cxiqLwBu9KzJf2orAnUoKSPiJisAf2OhdkTmKwKFlx3Mq6YnAIvOlicOYicCkgISfXEiJwCYOY7X194jAqJtBy46niMAqKSDhJ1eIwKy2/vbABojALkTdDFq2h8Cw0bsi82WHwDJfmjiMFYfAtOx4TiXFhsA2eldkvnSGwLgHNnpXJIbAOpUUkPDThcC8IvOliYOFwD6w0bsiM4XAwD2w0bvihMBCy47nVJKEwMRYbf3tQYTAReZLE4fxg8DHcyopIKGDwEkBCT+5UIPAy47nVFIAg8BNHMZq66+CwM+ppICEX4LAUDeDlh0PgsDSxGGstr6BwFRSQMJPboHA1t8e2OgdgcBYbf3tgc2AwNr62wMbfYDAXIi6GbQsgMC8KzJfmrh/wMBG74rMF3/AxGGstv52fsDIfGniMNZ9wMyXJg5jNX3A0LLjOZWUfMDUzaBlx/N7wNjoXZH5UnvA3AMbvSuyesDeHtjoXRF6wOI5lRSQcHnA5lRSQMLPeMDqbw9s9C54wO6KzJcmjnfA8qWJw1jtdsD2wEbvikx2wPrbAxu9q3XA/vbARu8KdcACEn5yIWp0wAYtO55TyXPACkj4yYUoc8AOY7X1t4dywBB+ciHq5nHAFJkvTRxGccAYtOx4TqVwwBzPqaSABHDAQNTNoGXHbsBICkj4yYVtwFBAwk8uRGzAWHY8p5ICa8BgrLb+9sBpwGjiMFZbf2jAcBirrb89Z8B4TiUFJPxlwICEn1yIumTAhLoZtOx4Y8CM8JMLUTdiwJQmDmO19WDAOLkQdTNoX8BIJQUk/ORcwFiR+dLEYVrAaP3tgY3eV8B4aeIwVltVwIjV1t8e2FLAmEHLjudUUMBQW397YKNLwHAzaNnxnEbAgAtRN4OWQcBAx3MqKSA5wADvisyXJi7AAJ9ciLoZFMAAoFyIuhkUQIDvisyXJi5AgMdzKikgOUCgC1E3g5ZBQIAzaNnxnEZAYFt/e2CjS0CgQcuO51RQQJDV1t8e2FJAgGniMFZbVUBw/e2Bjd5XQGCR+dLEYVpAUCUFJPzkXEBAuRB1M2hfQJgmDmO19WBAkPCTC1E3YkCIuhm07HhjQIiEn1yIumRAgE4lBST8ZUB4GKutvz1nQHDiMFZbf2hAaKy2/vbAaUBgdjynkgJrQFhAwk8uRGxAUApI+MmFbUBI1M2gZcduQCDPqaSABHBAHLTseE6lcEAYmS9NHEZxQBR+ciHq5nFAEGO19beHckAMSPjJhShzQAgtO55TyXNABBJ+ciFqdEAA98BG7wp1QPzbAxu9q3VA+MBG74pMdkD0pYnDWO12QPCKzJcmjndA7G8PbPQueEDoVFJAws94QOQ5lRSQcHlA5B7Y6F0RekDgAxu9K7J6QNzoXZH5UntA2M2gZcfze0DUsuM5lZR8QNCXJg5jNX1AzHxp4jDWfUDIYay2/nZ+QMRG74rMF39AwCsyX5q4f0BeiLoZtCyAQNz62wMbfYBAWm397YHNgEDY3x7Y6B2BQFZSQMJPboFA1MRhrLa+gUBSN4OWHQ+CQNCppICEX4JAThzGauuvgkDMjudUUgCDQEoBCT+5UINAyHMqKSChg0BG5ksTh/GDQMRYbf3tQYRAQsuO51SShEDAPbDRu+KEQECw0bsiM4VAviLzpYmDhUA8lRSQ8NOFQLoHNnpXJIZAOHpXZL50hkC27HhOJcWGQDRfmjiMFYdAstG7IvNlh0AwRN0MWraHQK62/vbABohALCkg4SdXiECqm0HLjqeIQCgOY7X194hApoCEn1xIiUAk86WJw5iJQKJlx3Mq6YlAINjoXZE5ikCeSgpI+ImKQBy9KzJf2opAmi9NHMYqi0AYom4GLXuLQJYUkPCTy4tAFIex2vobjECS+dLEYWyMQBBs9K7IvIxAkN4VmS8NjUAOUTeDll2NQIzDWG39rY1ACjZ6V2T+jUCIqJtBy06OQAYbvSsyn45AhI3eFZnvjkAAAAAAAECPQA==\"},\"shape\":[200],\"dtype\":\"float64\",\"order\":\"little\"}],[\"y\",{\"type\":\"ndarray\",\"array\":{\"type\":\"bytes\",\"data\":\"AAAAAAAAAIDimEXOaOMhwDC3clLjOynALZADz5vTLsAvs9+glcAxwIVV33j0yzPAzcee8FahNcCF8gXSfk03wEkefN3A2DjAbIvp0OlIOsA0hN8RNaI7wIev9z7W5zzAwllbbUscPsBZ46gRkUE/wNnrgh2iLEDA9wHn/FyyQMAX0VTuhzJBwML6JO6brUHAINeblQAkQsDcguhsD5ZCwOEsrnYWBEPA1tEAK1puQ8D77TwHF9VDwG6OMs2COETABRC1hM2YRMDjY4RNIvZEwCUj7QuoUEXA6QLs94GoRcCLSsUU0P1FwD5Zo5WvUEbA5lHLMjuhRsAffDNyi+9GwGlIsuW2O0fARM+KYNKFR8A0Q8Ml8c1HwCvDbBAlFEjACJ3Ntn5YSMC77TOJDZtIwBwXFO3f20jA4sz6VAMbScBRDcRVhFhJwCgZdrlulEnAX40QkM3OScAqdpM+qwdKwPYHeIwRP0rA/zjMrwl1SsD7aBtYnKlKwN5bSLjR3ErAS9N3j7EOS8Bu3iYxQz9LwLRnhIyNbkvA+H4iM5ecS8CSRxFfZslLwEM5c/gA9UvAZnqYmmwfTMB7dK2YrkhMwO9XBwLMcEzA/wcYpsmXTMAL1REYrL1MwA18QrJ34kzA3BAsmTAGTcBFwWG+2ihNwDm+LeN5Sk3AdAoFmxFrTcAPcc1NpYpNwNx3+Tk4qU3AkL99ds3GTcDS7KL0Z+NNwO7ntoEK/03AHPyfyLcZTsC2IVRTcjNOwLKHNow8TE7AfT9dvxhkTsDmwcAbCXtOwBDcVrQPkU7AhG8agS6mTsDzTwFgZ7pOwBxs4RW8zU7AkFRGTy7gTsAbGzihv/FOwJ5v9IlxAk/AKMuacUUST8D2Z8yqPCFPwEq0QHNYL0/AtN5O9Jk8T8BrDG1DAklPwGW9pWKSVE/Az9UDQUtfT8COufW6LWlPwCPcp5o6ck/Aph1XmHJ6T8C5RJta1oFPwMfcqXZmiE/AFriRcCOOT8CqTW+7DZNPwN8jmrkll0/AznHLvGuaT8AlGz4G4JxPwLgjyMaCnk/ABbLtHlSfT8AFsu0eVJ9PwLgjyMaCnk/AJRs+BuCcT8DOccu8a5pPwN8jmrkll0/Aqk1vuw2TT8AWuJFwI45PwMfcqXZmiE/AuUSbWtaBT8CmHVeYcnpPwCPcp5o6ck/Ajrn1ui1pT8DP1QNBS19PwGW9pWKSVE/AawxtQwJJT8Cz3k70mTxPwEq0QHNYL0/A9WfMqjwhT8Any5pxRRJPwJ1v9IlxAk/AGxs4ob/xTsCPVEZPLuBOwBxs4RW8zU7A808BYGe6TsCDbxqBLqZOwA/cVrQPkU7A5cHAGwl7TsB9P12/GGROwLGHNow8TE7AtSFUU3IzTsAc/J/ItxlOwO7ntoEK/03A0uyi9GfjTcCQv312zcZNwNt3+Tk4qU3AD3HNTaWKTcB0CgWbEWtNwDi+LeN5Sk3ARcFhvtooTcDbECyZMAZNwAx8QrJ34kzACtURGKy9TMD/BximyZdMwO5XBwLMcEzAenStmK5ITMBlepiabB9MwEM5c/gA9UvAkUcRX2bJS8D3fiIzl5xLwLNnhIyNbkvAbN4mMUM/S8BK03ePsQ5LwN1bSLjR3ErA+mgbWJypSsD+OMyvCXVKwPUHeIwRP0rAKHaTPqsHSsBejRCQzc5JwCcZdrlulEnAUA3EVYRYScDizPpUAxtJwBwXFO3f20jAuu0ziQ2bSMAGnc22flhIwCvDbBAlFEjANEPDJfHNR8BEz4pg0oVHwGZIsuW2O0fAHXwzcovvRsDlUcsyO6FGwDxZo5WvUEbAiUrFFND9RcDnAuz3gahFwCMj7QuoUEXA4GOETSL2RMADELWEzZhEwGuOMs2COETA+e08BxfVQ8DU0QArWm5DwN8srnYWBEPA2oLobA+WQsAc15uVACRCwMH6JO6brUHAFdFU7ocyQcD2Aef8XLJAwNjrgh2iLEDAVeOoEZFBP8C9WVttSxw+wIWv9z7W5zzAMoTfETWiO8Bni+nQ6Ug6wEQefN3A2DjAfPIF0n5NN8DBx57wVqE1wHtV33j0yzPAILPfoJXAMcAMkAPPm9MuwAa3clLjOynAtphFzmjjIcAAAAAAAAAAgA==\"},\"shape\":[200],\"dtype\":\"float64\",\"order\":\"little\"}]]}}},\"view\":{\"type\":\"object\",\"name\":\"CDSView\",\"id\":\"p1056\",\"attributes\":{\"filter\":{\"type\":\"object\",\"name\":\"AllIndices\",\"id\":\"p1057\"}}},\"glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1052\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#ff7f0e\",\"line_width\":2,\"line_dash\":[6]}},\"nonselection_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1053\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.1,\"line_width\":2,\"line_dash\":[6]}},\"muted_glyph\":{\"type\":\"object\",\"name\":\"Line\",\"id\":\"p1054\",\"attributes\":{\"x\":{\"type\":\"field\",\"field\":\"x\"},\"y\":{\"type\":\"field\",\"field\":\"y\"},\"line_color\":\"#ff7f0e\",\"line_alpha\":0.2,\"line_width\":2,\"line_dash\":[6]}}}}],\"toolbar\":{\"type\":\"object\",\"name\":\"Toolbar\",\"id\":\"p1014\",\"attributes\":{\"tools\":[{\"type\":\"object\",\"name\":\"PanTool\",\"id\":\"p1027\"},{\"type\":\"object\",\"name\":\"WheelZoomTool\",\"id\":\"p1028\",\"attributes\":{\"renderers\":\"auto\"}},{\"type\":\"object\",\"name\":\"BoxZoomTool\",\"id\":\"p1029\",\"attributes\":{\"dimensions\":\"both\",\"overlay\":{\"type\":\"object\",\"name\":\"BoxAnnotation\",\"id\":\"p1030\",\"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\":\"p1036\",\"attributes\":{\"all\":{\"type\":\"object\",\"name\":\"AreaVisuals\",\"id\":\"p1035\",\"attributes\":{\"fill_color\":\"white\",\"hover_fill_color\":\"lightgray\"}}}}}}}},{\"type\":\"object\",\"name\":\"SaveTool\",\"id\":\"p1037\"},{\"type\":\"object\",\"name\":\"ResetTool\",\"id\":\"p1038\"},{\"type\":\"object\",\"name\":\"HelpTool\",\"id\":\"p1039\"}]}},\"toolbar_location\":null,\"left\":[{\"type\":\"object\",\"name\":\"LinearAxis\",\"id\":\"p1022\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"FixedTicker\",\"id\":\"p1061\",\"attributes\":{\"ticks\":[-50,-25,0],\"minor_ticks\":[]}},\"formatter\":{\"type\":\"object\",\"name\":\"BasicTickFormatter\",\"id\":\"p1024\"},\"axis_label\":\"z (m) - VE:5x\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1025\"}}}],\"below\":[{\"type\":\"object\",\"name\":\"LinearAxis\",\"id\":\"p1017\",\"attributes\":{\"ticker\":{\"type\":\"object\",\"name\":\"FixedTicker\",\"id\":\"p1060\",\"attributes\":{\"ticks\":[-1000,-500,0,500,1000],\"minor_ticks\":[]}},\"formatter\":{\"type\":\"object\",\"name\":\"BasicTickFormatter\",\"id\":\"p1019\"},\"axis_label\":\"x (m)\",\"axis_label_text_font_style\":\"normal\",\"major_label_policy\":{\"type\":\"object\",\"name\":\"AllLabels\",\"id\":\"p1020\"}}}],\"center\":[{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1021\",\"attributes\":{\"axis\":{\"id\":\"p1017\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Grid\",\"id\":\"p1026\",\"attributes\":{\"dimension\":1,\"axis\":{\"id\":\"p1022\"},\"grid_line_color\":\"#b0b0b0\"}},{\"type\":\"object\",\"name\":\"Legend\",\"id\":\"p1058\",\"attributes\":{\"location\":\"bottom_right\",\"items\":[{\"type\":\"object\",\"name\":\"LegendItem\",\"id\":\"p1059\",\"attributes\":{\"label\":{\"type\":\"value\",\"value\":\"Dupuit\"},\"renderers\":[{\"id\":\"p1055\"}]}}]}}]}},{\"type\":\"object\",\"name\":\"Slider\",\"id\":\"p1062\",\"attributes\":{\"js_property_callbacks\":{\"type\":\"map\",\"entries\":[[\"change:value\",[{\"type\":\"object\",\"name\":\"CustomJS\",\"id\":\"p1063\",\"attributes\":{\"args\":{\"type\":\"map\",\"entries\":[[\"source\",{\"id\":\"p1003\"}],[\"slider\",{\"id\":\"p1062\"}]]},\"code\":\"\\n const data = source.data;\\n const kz_kx = slider.value;\\n \\n const kx = 10;\\n const kz = kx * kz_kx;\\n const L = 1000;\\n const N = 0.001;\\n const rhof = 1000;\\n const rhos = 1025;\\n const alpha = rhof / (rhos - rhof);\\n const ng = 201;\\n \\n function linspace(start, stop, num) {\\n const arr = [];\\n const step = (stop - start) / (num - 1);\\n for (let i = 0; i < num; i++) {\\n arr.push(start + step * i);\\n }\\n return arr;\\n }\\n \\n const xtip = L * Math.sqrt(1 + N * alpha / kz);\\n const xi = linspace(-xtip + 1e-12, xtip - 1e-12, ng);\\n const zi = xi.map(x => {\\n const term = N * alpha / kx * L * L - N * alpha / (N * alpha + kz) * kz / kx * x * x;\\n return -Math.sqrt(term);\\n });\\n \\n data['x_values'] = xi;\\n data['z_values'] = zi;\\n source.change.emit();\\n \"}}]]]},\"width\":570,\"align\":\"center\",\"title\":\"kz/kx\",\"start\":0.01,\"end\":1,\"value\":1,\"step\":0.01}},{\"type\":\"object\",\"name\":\"Div\",\"id\":\"p1064\",\"attributes\":{\"align\":\"center\",\"text\":\"

^0.1⁄₁₀ | ²⁄₁₀ | ⁴⁄₁₀ | ⁶⁄₁₀ | ⁸⁄₁₀ | ¹⁰⁄₁₀

\"}}]}}]}};\n", " const render_items = [{\"docid\":\"1ede088a-211a-4a24-96cc-cbb881efb049\",\"roots\":{\"p1066\":\"c169144d-cbf2-4cf0-a211-7c7640cbc4f5\"},\"root_ids\":[\"p1066\"]}];\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": "p1066" } }, "output_type": "display_data" } ], "source": [ "inter_fig_chap10.fig_10_22()" ] }, { "cell_type": "markdown", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "```{raw} html\n", "\n", "```\n", "\n", "```{figure} ../py_figures/1x1_pixel.png\n", "---\n", "name: fig_compinterface10.6\n", "---\n", "Comparison of the interface position below a strip island for the anisotropic, two-dimensional solution and the Dupuit solution. \n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```{exercise}\n", ":label: ex_10.8\n", "Compute the ratio of the depth of the interface and the head at the interface at $x=0$ for the anisotropic case. Next, compute ratio of the depth of the interface and the head at the top of the aquifer at $x=0$ for the anisotropic case. \n", "\n", "

\n", "

Solution

\n", "\n", "Factor based on the head at the interface at $x=0$ m: $40$\\\n", "Factor based on the head at the top of the aquifer at $x=0$ m: $38.5$\n", "\n", "

\n", "```" ] } ], "metadata": { "celltoolbar": "Tags", "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.7" }, "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 }