|
14 | 14 | "metadata": {}, |
15 | 15 | "outputs": [], |
16 | 16 | "source": [ |
17 | | - "import numpy as np\n", |
18 | 17 | "import matplotlib.pyplot as plt\n", |
| 18 | + "import numpy as np\n", |
| 19 | + "\n", |
19 | 20 | "import timml as tml" |
20 | 21 | ] |
21 | 22 | }, |
|
59 | 60 | "metadata": {}, |
60 | 61 | "outputs": [], |
61 | 62 | "source": [ |
62 | | - "l = 40.0 # length building pit in m\n", |
63 | | - "b = 30.0 # width building pit in m\n", |
| 63 | + "length = 40.0 # length building pit in m\n", |
| 64 | + "width = 30.0 # width building pit in m\n", |
64 | 65 | "\n", |
65 | 66 | "h_bem = -6.21 # m\n", |
66 | 67 | "offset = 5.0 # distance groundwater extraction element from sheetpiles in m" |
|
73 | 74 | "outputs": [], |
74 | 75 | "source": [ |
75 | 76 | "xy = [\n", |
76 | | - " (-l / 2, -b / 2),\n", |
77 | | - " (l / 2, -b / 2),\n", |
78 | | - " (l / 2, b / 2),\n", |
79 | | - " (-l / 2, b / 2),\n", |
80 | | - " (-l / 2, -b / 2),\n", |
| 77 | + " (-length / 2, -width / 2),\n", |
| 78 | + " (length / 2, -width / 2),\n", |
| 79 | + " (length / 2, width / 2),\n", |
| 80 | + " (-length / 2, width / 2),\n", |
| 81 | + " (-length / 2, -width / 2),\n", |
81 | 82 | "]" |
82 | 83 | ] |
83 | 84 | }, |
|
211 | 212 | "\n", |
212 | 213 | "tml.HeadLineSink(\n", |
213 | 214 | " ml,\n", |
214 | | - " x1=-l / 2 + offset,\n", |
215 | | - " y1=b / 2 - offset,\n", |
216 | | - " x2=l / 2 - offset,\n", |
217 | | - " y2=b / 2 - offset,\n", |
| 215 | + " x1=-length / 2 + offset,\n", |
| 216 | + " y1=width / 2 - offset,\n", |
| 217 | + " x2=length / 2 - offset,\n", |
| 218 | + " y2=width / 2 - offset,\n", |
218 | 219 | " hls=h_bem,\n", |
219 | 220 | " layers=np.arange(last_lay_dw + 1),\n", |
220 | 221 | ")\n", |
221 | 222 | "tml.HeadLineSink(\n", |
222 | 223 | " ml,\n", |
223 | | - " x1=-l / 2 + offset,\n", |
| 224 | + " x1=-length / 2 + offset,\n", |
224 | 225 | " y1=0,\n", |
225 | | - " x2=l / 2 - offset,\n", |
| 226 | + " x2=length / 2 - offset,\n", |
226 | 227 | " y2=0,\n", |
227 | 228 | " hls=h_bem,\n", |
228 | 229 | " layers=np.arange(last_lay_dw + 1),\n", |
229 | 230 | ")\n", |
230 | 231 | "tml.HeadLineSink(\n", |
231 | 232 | " ml,\n", |
232 | | - " x1=-l / 2 + offset,\n", |
233 | | - " y1=-b / 2 + offset,\n", |
234 | | - " x2=l / 2 - offset,\n", |
235 | | - " y2=-b / 2 + offset,\n", |
| 233 | + " x1=-length / 2 + offset,\n", |
| 234 | + " y1=-width / 2 + offset,\n", |
| 235 | + " x2=length / 2 - offset,\n", |
| 236 | + " y2=-width / 2 + offset,\n", |
236 | 237 | " hls=h_bem,\n", |
237 | 238 | " layers=np.arange(last_lay_dw + 1),\n", |
238 | 239 | ")\n", |
|
248 | 249 | "\n", |
249 | 250 | "print(\"\\nDischarge =\", np.round(Qtot.sum(), 2), \"m^3/dag\")\n", |
250 | 251 | "\n", |
251 | | - "y = np.linspace(-b / 2 - 25, b / 2 + 1100, 201)\n", |
| 252 | + "y = np.linspace(-width / 2 - 25, width / 2 + 1100, 201)\n", |
252 | 253 | "hl = ml.headalongline(np.zeros(201), y, layers=[0])\n", |
253 | 254 | "y_5cm = np.interp(\n", |
254 | 255 | " -0.05, ml.headalongline(np.zeros(201), y, layers=0).squeeze(), y, right=np.nan\n", |
|
298 | 299 | "metadata": {}, |
299 | 300 | "outputs": [], |
300 | 301 | "source": [ |
301 | | - "x = np.linspace(0.0, l / 2 + 1100, 201)\n", |
| 302 | + "x = np.linspace(0.0, length / 2 + 1100, 201)\n", |
302 | 303 | "hl = ml.headalongline(x, np.zeros(201), layers=[last_lay_dw, last_lay_dw + 1])" |
303 | 304 | ] |
304 | 305 | }, |
|
337 | 338 | "metadata": {}, |
338 | 339 | "outputs": [], |
339 | 340 | "source": [ |
340 | | - "x = np.linspace(-l / 2 - 25, 0.0, 201)\n", |
| 341 | + "x = np.linspace(-length / 2 - 25, 0.0, 201)\n", |
341 | 342 | "hl = ml.headalongline(x, np.zeros(201), layers=[last_lay_dw, last_lay_dw + 1])" |
342 | 343 | ] |
343 | 344 | }, |
|
378 | 379 | "source": [ |
379 | 380 | "xoffset = 50\n", |
380 | 381 | "zoffset = 15\n", |
381 | | - "x1, x2, y1, y2 = [-l / 2 - xoffset, 0.0, 0, 0]\n", |
| 382 | + "x1, x2, y1, y2 = [-length / 2 - xoffset, 0.0, 0, 0]\n", |
382 | 383 | "nudge = 1e-6\n", |
383 | 384 | "n = 301" |
384 | 385 | ] |
|
581 | 582 | "\n", |
582 | 583 | "tml.HeadLineSink(\n", |
583 | 584 | " ml,\n", |
584 | | - " x1=-l / 2 + offset,\n", |
585 | | - " y1=b / 2 - offset,\n", |
586 | | - " x2=l / 2 - offset,\n", |
587 | | - " y2=b / 2 - offset,\n", |
| 585 | + " x1=-length / 2 + offset,\n", |
| 586 | + " y1=width / 2 - offset,\n", |
| 587 | + " x2=length / 2 - offset,\n", |
| 588 | + " y2=width / 2 - offset,\n", |
588 | 589 | " hls=h_bem,\n", |
589 | 590 | " layers=wlayers,\n", |
590 | 591 | ")\n", |
591 | 592 | "tml.HeadLineSink(\n", |
592 | 593 | " ml,\n", |
593 | | - " x1=-l / 2 + offset,\n", |
| 594 | + " x1=-length / 2 + offset,\n", |
594 | 595 | " y1=0,\n", |
595 | | - " x2=l / 2 - offset,\n", |
| 596 | + " x2=length / 2 - offset,\n", |
596 | 597 | " y2=0,\n", |
597 | 598 | " hls=h_bem,\n", |
598 | 599 | " layers=wlayers,\n", |
599 | 600 | " order=5,\n", |
600 | 601 | ")\n", |
601 | 602 | "tml.HeadLineSink(\n", |
602 | 603 | " ml,\n", |
603 | | - " x1=-l / 2 + offset,\n", |
604 | | - " y1=-b / 2 + offset,\n", |
605 | | - " x2=l / 2 - offset,\n", |
606 | | - " y2=-b / 2 + offset,\n", |
| 604 | + " x1=-length / 2 + offset,\n", |
| 605 | + " y1=-width / 2 + offset,\n", |
| 606 | + " x2=length / 2 - offset,\n", |
| 607 | + " y2=-width / 2 + offset,\n", |
607 | 608 | " hls=h_bem,\n", |
608 | 609 | " layers=wlayers,\n", |
609 | 610 | ")\n", |
|
619 | 620 | "\n", |
620 | 621 | "print(\"\\nDischarge =\", np.round(Qtot.sum(), 2), \"m^3/dag\")\n", |
621 | 622 | "\n", |
622 | | - "y = np.linspace(-b / 2 - 25, b / 2 + 1100, 201)\n", |
| 623 | + "y = np.linspace(-width / 2 - 25, width / 2 + 1100, 201)\n", |
623 | 624 | "hl = ml.headalongline(np.zeros(201), y, layers=[0])\n", |
624 | 625 | "y_5cm = np.interp(\n", |
625 | 626 | " -0.05, ml.headalongline(np.zeros(201), y, layers=0).squeeze(), y, right=np.nan\n", |
|
642 | 643 | "metadata": {}, |
643 | 644 | "outputs": [], |
644 | 645 | "source": [ |
645 | | - "x = np.linspace(0.0, l / 2 + 1100, 201)\n", |
| 646 | + "x = np.linspace(0.0, length / 2 + 1100, 201)\n", |
646 | 647 | "hl = ml.headalongline(x, np.zeros(201), layers=[0, last_lay_dw, last_lay_dw + 1])" |
647 | 648 | ] |
648 | 649 | }, |
|
682 | 683 | "metadata": {}, |
683 | 684 | "outputs": [], |
684 | 685 | "source": [ |
685 | | - "x = np.linspace(-l / 2 - 25, 0.0, 201)\n", |
| 686 | + "x = np.linspace(-length / 2 - 25, 0.0, 201)\n", |
686 | 687 | "hl = ml.headalongline(x, np.zeros(201), layers=[0, last_lay_dw, last_lay_dw + 1])" |
687 | 688 | ] |
688 | 689 | }, |
|
724 | 725 | "source": [ |
725 | 726 | "xoffset = 50\n", |
726 | 727 | "zoffset = 15\n", |
727 | | - "x1, x2, y1, y2 = [-l / 2 - xoffset, 0.0, 0, 0]\n", |
| 728 | + "x1, x2, y1, y2 = [-length / 2 - xoffset, 0.0, 0, 0]\n", |
728 | 729 | "nudge = 1e-6\n", |
729 | 730 | "n = 301" |
730 | 731 | ] |
|
0 commit comments