From c443e9c535dc47936e8d64902c2524201d208997 Mon Sep 17 00:00:00 2001 From: Joseph D Hughes Date: Tue, 17 Aug 2021 13:49:04 -0500 Subject: [PATCH 1/2] refactor(plot): remove deprecated dis and disu thickness property --- .docs/main.rst | 2 +- .docs/tutorials.rst | 12 + autotest/t007_test.py | 26 +- examples/Notebooks/flopy3_export.ipynb | 102 +-- .../Notebooks/flopy3_shapefile_export.ipynb | 200 +++-- .../Notebooks/flopy3_working_stack_demo.ipynb | 707 +++--------------- flopy/modflow/mfdis.py | 34 - flopy/modflow/mfdisu.py | 39 - 8 files changed, 242 insertions(+), 880 deletions(-) diff --git a/.docs/main.rst b/.docs/main.rst index 95aa7948bc..0c21558f5a 100644 --- a/.docs/main.rst +++ b/.docs/main.rst @@ -6,7 +6,7 @@ FloPy Documentation =================== -**Documentation for version 3.3.4 --- release candidate** +**Documentation for version 3.3.5 --- release candidate** Return to the Github `FloPy `_ website where the code resides. diff --git a/.docs/tutorials.rst b/.docs/tutorials.rst index 0e82787392..c113a69c06 100644 --- a/.docs/tutorials.rst +++ b/.docs/tutorials.rst @@ -53,6 +53,18 @@ Contents: _notebooks/tutorial02_mf +.IPYNB_CHECKPOINTS Tutorials +---------------------------- + +Contents: + +.. toctree:: + :maxdepth: 2 + + _notebooks/tutorial01_mf-checkpoint + _notebooks/tutorial02_mf-checkpoint + + SEAWAT Tutorials ---------------- diff --git a/autotest/t007_test.py b/autotest/t007_test.py index 865bee4d50..db2016c686 100644 --- a/autotest/t007_test.py +++ b/autotest/t007_test.py @@ -735,7 +735,7 @@ def test_mg(): botm=botm, ) bas = flopy.modflow.ModflowBas(ms, ifrefm=True) - t = dis.thickness + t = ms.modelgrid.thick # test instantiation of an empty basic Structured Grid mg = flopy.discretization.StructuredGrid(dis.delc.array, dis.delr.array) @@ -1503,20 +1503,28 @@ def test_model_dot_plot(): import matplotlib.pyplot as plt loadpth = os.path.join("..", "examples", "data", "mf2005_test") - ml = flopy.modflow.Modflow.load("ibs2k.nam", "mf2k", model_ws=loadpth) + ml = flopy.modflow.Modflow.load( + "ibs2k.nam", "mf2k", model_ws=loadpth, check=False + ) ax = ml.plot() - assert isinstance(ax, list) - assert len(ax) == 20 + assert isinstance(ax, list), "ml.plot() ax is is not a list" + assert len(ax) == 18, "number of axes ({}) is " "not equal to 18".format( + len(ax) + ) plt.close("all") # plot specific dataset ax = ml.bcf6.hy.plot() - assert isinstance(ax, list) - assert len(ax) == 2 + assert isinstance(ax, list), "ml.bcf6.hy.plot() ax is is not a list" + assert len(ax) == 2, "number of hy axes ({}) " "is not equal to 2".format( + len(ax) + ) # special case where nlay != plottable ax = ml.bcf6.vcont.plot() - assert isinstance(ax, plt.Axes) + assert isinstance( + ax, plt.Axes + ), "ml.bcf6.vcont.plot() ax is is not of type plt.Axes" plt.close("all") @@ -1783,12 +1791,12 @@ def main(): # test_netcdf_classmethods() # build_netcdf() # build_sfr_netcdf() - test_twri_mg() + # test_twri_mg() # test_mg() # test_mbase_modelgrid() # test_mt_modelgrid() # test_rotation() - # test_model_dot_plot() + test_model_dot_plot() # test_get_lrc_get_node() # test_vertex_model_dot_plot() # test_sr_with_Map() diff --git a/examples/Notebooks/flopy3_export.ipynb b/examples/Notebooks/flopy3_export.ipynb index c31812584c..4a99a27c36 100644 --- a/examples/Notebooks/flopy3_export.ipynb +++ b/examples/Notebooks/flopy3_export.ipynb @@ -18,9 +18,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "3.8.10 (default, May 19 2021, 11:01:55) \n", + "3.8.11 (default, Aug 6 2021, 08:56:27) \n", "[Clang 10.0.0 ]\n", - "flopy version: 3.3.4\n" + "flopy version: 3.3.5\n" ] } ], @@ -57,7 +57,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/mbase.py:354: DeprecationWarning: xul/yul have been deprecated. Use xll/yll instead.\n", + "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/mbase.py:352: DeprecationWarning: xul/yul have been deprecated. Use xll/yll instead.\n", " warnings.warn(\n" ] } @@ -195,8 +195,6 @@ "text": [ "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/reference.py:2074: DeprecationWarning: crs has been deprecated and will be removed in version 3.3.5. Use CRS in shapefile_utils instead.\n", " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n" ] @@ -222,7 +220,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 8, @@ -357,7 +355,9 @@ "output_type": "stream", "text": [ "initialize_geometry::proj4_str = +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs\n", - "initialize_geometry::self.grid_crs = +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +type=crs\n" + "initialize_geometry::self.grid_crs = +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +type=crs\n", + "initialize_geometry::nc_crs = epsg:4326\n", + "transforming coordinates using = proj=noop ellps=GRS80\n" ] }, { @@ -368,18 +368,10 @@ " warnings.warn(\n" ] }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "initialize_geometry::nc_crs = epsg:4326\n", - "transforming coordinates using = proj=noop ellps=GRS80\n" - ] - }, { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 12, @@ -434,8 +426,8 @@ "text/plain": [ "\n", "root group (NETCDF4 data model, file format HDF5):\n", - " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.4\n", - " date_created: 2021-08-06T21:19:00Z\n", + " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.5\n", + " date_created: 2021-08-17T18:19:00Z\n", " geospatial_vertical_positive: up\n", " geospatial_vertical_min: -25.0\n", " geospatial_vertical_max: 4.832500457763672\n", @@ -500,8 +492,6 @@ "text": [ "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/reference.py:2074: DeprecationWarning: crs has been deprecated and will be removed in version 3.3.5. Use CRS in shapefile_utils instead.\n", " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n" ] @@ -511,8 +501,8 @@ "text/plain": [ "\n", "root group (NETCDF4 data model, file format HDF5):\n", - " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.4\n", - " date_created: 2021-08-06T21:19:00Z\n", + " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.5\n", + " date_created: 2021-08-17T18:19:00Z\n", " geospatial_vertical_positive: up\n", " geospatial_vertical_min: -25.0\n", " geospatial_vertical_max: 4.832500457763672\n", @@ -585,13 +575,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n", "error getting data for cell_by_cell_flowstorage at time 1097.0:list index out of range\n" ] }, @@ -600,8 +584,8 @@ "text/plain": [ "\n", "root group (NETCDF4 data model, file format HDF5):\n", - " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.4\n", - " date_created: 2021-08-06T21:19:00Z\n", + " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.5\n", + " date_created: 2021-08-17T18:19:00Z\n", " geospatial_vertical_positive: up\n", " geospatial_vertical_min: -25.0\n", " geospatial_vertical_max: 4.832500457763672\n", @@ -738,13 +722,7 @@ "name": "stdout", "output_type": "stream", "text": [ - "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n", "error getting data for cell_by_cell_flowstorage at time 1097.0:list index out of range\n" ] }, @@ -755,21 +733,9 @@ "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/zonbud.py:3224: PendingDeprecationWarning: ZoneBudgetOutput will be deprecated in version 3.3.5\n", " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/zonbud.py:3198: PendingDeprecationWarning: ZoneBudgetOutput will be deprecated in version 3.3.5\n", - " warnings.warn(\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ + " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/zonbud.py:3292: PendingDeprecationWarning: ZoneBudgetOutput will be deprecated in version 3.3.5\n", - " warnings.warn(\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ + " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/zonbud.py:3174: PendingDeprecationWarning: ZoneBudgetOutput will be deprecated in version 3.3.5\n", " warnings.warn(\n" ] @@ -779,8 +745,8 @@ "text/plain": [ "\n", "root group (NETCDF4 data model, file format HDF5):\n", - " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.4\n", - " date_created: 2021-08-06T21:19:00Z\n", + " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.5\n", + " date_created: 2021-08-17T18:19:00Z\n", " geospatial_vertical_positive: up\n", " geospatial_vertical_min: -25.0\n", " geospatial_vertical_max: 4.832500457763672\n", @@ -1108,9 +1074,7 @@ { "cell_type": "code", "execution_count": 19, - "metadata": { - "scrolled": false - }, + "metadata": {}, "outputs": [ { "name": "stderr", @@ -1446,20 +1410,8 @@ "output_type": "stream", "text": [ "initialize_geometry::nc_crs = epsg:4326\n", - "transforming coordinates using = proj=noop ellps=GRS80\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "transforming coordinates using = proj=noop ellps=GRS80\n", + "error getting data for cell_by_cell_flowstorage at time 1.0:list index out of range\n", "error getting data for cell_by_cell_flowstorage at time 1097.0:list index out of range\n" ] }, @@ -1468,8 +1420,8 @@ "text/plain": [ "\n", "root group (NETCDF4 data model, file format HDF5):\n", - " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.4\n", - " date_created: 2021-08-06T21:19:00Z\n", + " Conventions: CF-1.6, ACDD-1.3, flopy 3.3.5\n", + " date_created: 2021-08-17T18:19:00Z\n", " geospatial_vertical_positive: up\n", " geospatial_vertical_min: -25.0\n", " geospatial_vertical_max: 4.832500457763672\n", @@ -1538,9 +1490,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.10" + "version": "3.8.11" } }, "nbformat": 4, - "nbformat_minor": 1 + "nbformat_minor": 4 } diff --git a/examples/Notebooks/flopy3_shapefile_export.ipynb b/examples/Notebooks/flopy3_shapefile_export.ipynb index 59ba549d6f..122fbd7ad0 100644 --- a/examples/Notebooks/flopy3_shapefile_export.ipynb +++ b/examples/Notebooks/flopy3_shapefile_export.ipynb @@ -21,11 +21,11 @@ "name": "stdout", "output_type": "stream", "text": [ - "3.8.10 (default, May 19 2021, 11:01:55) \n", + "3.8.11 (default, Aug 6 2021, 08:56:27) \n", "[Clang 10.0.0 ]\n", "numpy version: 1.19.2\n", "matplotlib version: 3.4.2\n", - "flopy version: 3.3.4\n" + "flopy version: 3.3.5\n" ] } ], @@ -145,10 +145,6 @@ "name": "stderr", "output_type": "stream", "text": [ - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", @@ -210,13 +206,6 @@ "execution_count": 8, "metadata": {}, "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "wrote temp/shapefile_export/wel.shp\n" - ] - }, { "name": "stderr", "output_type": "stream", @@ -226,6 +215,13 @@ "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n" ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "wrote temp/shapefile_export/wel.shp\n" + ] } ], "source": [ @@ -255,7 +251,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 9, @@ -321,7 +317,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 12, @@ -551,52 +547,52 @@ { "data": { "text/plain": [ - "[,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ]" + "[,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ]" ] }, "execution_count": 18, @@ -712,49 +708,49 @@ " \n", " 0\n", " 0\n", - " -1172.274885\n", - " 277780.407448\n", - " 5.093921e+06\n", + " -923.150882\n", + " 274071.891678\n", + " 5.091823e+06\n", " \n", " \n", " 1\n", " 1\n", - " -1013.143739\n", - " 276089.394787\n", - " 5.095465e+06\n", + " -962.082825\n", + " 277320.002989\n", + " 5.089461e+06\n", " \n", " \n", " 2\n", " 2\n", - " -800.181374\n", - " 274300.382227\n", - " 5.092448e+06\n", + " -981.479964\n", + " 276032.262286\n", + " 5.098156e+06\n", " \n", " \n", " 3\n", " 3\n", - " -1012.078046\n", - " 277121.455475\n", - " 5.089068e+06\n", + " -843.991158\n", + " 277428.015559\n", + " 5.088979e+06\n", " \n", " \n", " 4\n", " 4\n", - " -912.590841\n", - " 275054.084064\n", - " 5.095245e+06\n", + " -898.895335\n", + " 277774.557126\n", + " 5.095573e+06\n", " \n", " \n", "\n", "" ], "text/plain": [ - " wellID q x_utm y_utm\n", - "0 0 -1172.274885 277780.407448 5.093921e+06\n", - "1 1 -1013.143739 276089.394787 5.095465e+06\n", - "2 2 -800.181374 274300.382227 5.092448e+06\n", - "3 3 -1012.078046 277121.455475 5.089068e+06\n", - "4 4 -912.590841 275054.084064 5.095245e+06" + " wellID q x_utm y_utm\n", + "0 0 -923.150882 274071.891678 5.091823e+06\n", + "1 1 -962.082825 277320.002989 5.089461e+06\n", + "2 2 -981.479964 276032.262286 5.098156e+06\n", + "3 3 -843.991158 277428.015559 5.088979e+06\n", + "4 4 -898.895335 277774.557126 5.095573e+06" ] }, "execution_count": 21, @@ -807,7 +803,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOkAAAEICAYAAAC6WNNwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAeRklEQVR4nO3deZwV1Zn/8c+3F2g2aWRTaQHFPUYdpgET4xZ3JgnqjBmzjEuMDom/RJIxGs1MJnmNictMfkmc/H4qCTExLnELiRmVSIxiYoKkUREQFQSUll1AwGbp5Zk/zmkp2ts0S/et0/Tzfr3uq6vq1Kl6qm49dU7VvbdaZoZzLl0leQfgnNsxT1LnEudJ6lziPEmdS5wnqXOJ8yR1LnF7TZJKMkmH5LDewyW9IGmDpC9Lul3Sv8WyUyTVFjumzkzSJZL+tBv1npb0+Y6IaRdi+Jaku9t7uTuVpJIWSzq9vVe+syRNlHRFXutvwzXA02bWx8xuNbPxZvYfeQfV3naUPJJek3RYsWPaXXkfz7uqs7SkZwOP5R1EK4YBc/MOoiNJKttB2QigxMxeK2JIXUqbSSrpF8BQ4LeSNkq6RtLxkv4saZ2kWZJOycz/tKQbYvlGSb+V1F/SPZLWS/qrpOGZ+S12ExdKWi3pPyWVZMqPAdaZWa2kQyRNk/ROnPf+FuGeLmm+pLWS/p8kxWWMkPQHSW/HevdIqsysY7Gk6yS9HOveKakiU/4xSS/G7f1zjAlJfwBOBX4Ut/UwST+TdEMr+/IASQ9LWiVpkaQv78T+L5H0dUmvx/gfkLRvLLtN0kOZeW+W9KSCUyTVSro+bvNiSZ/JzNtX0l0xljck/Wvzfo+t5rOSvi9pDXA/cDvwobid6zIh/h3wmKSD4v5pXsZPJK3MrO9uSRMy654kaZmkt+LxUtrWvmixX86Q9Eo8Fn4EKFPW6vtd6HiO0x+UtDwu7xlJH9jBuq+NcW+Q9Kqk0zLF3eJ+3SBprqTqTL0dHmetMrM2X8Bi4PQ4PAR4GxhLSPIz4vjAWP40sAAYAfQFXgZeA04HyoC7gDszyzbgKWDfuPNeAz6fKf86cGMcvg/4RlxvBfCRFsv5H6AyLmcVcHYsOyTG2R0YCDwD/KDF9s0BDoxxPAvcEMtGAiuBMUApcHGcv3tme7Px/ixT9xSgNg6XADOBbwLdgIOBhcBZbez7CcB0oCrGfwdwXyzrGffXJcCJwGqgKrPuBuD/xnonA+8Ch8fyu4DfAH2A4XE5l8WyS2LdL8X3rEec9qcC8U1p3gbgTeBv4/CrcfuOzJT9TRz+ddyOXsAgYAbwz5l1v289LdY5AFgP/ANQDnwlxvv5XXi/T2+xzM/FfdEd+AHwYivrPhxYAhwQx4cDI+Lwt4DNhNwoBW4Epu/McbbD7d2NJL0W+EWL8t8BF2cO2m9kyr4HPJ4Z/3h2BxCS6+zM+BeBJzPjfwROzBxYE4kHYosYjO2T9gHg661sz7nAC8BPCQm4FRifKR8LvB6HbwP+Iw5/knDS2QL8fheTdAzwZos4riNzwmol1nnAaZnx/YF6oCyOjwbWAG8An8rMdwrhwO3VYp/8WzyAtgBHZcr+mXBtDSFRWsZ6CS2Sh3CSeBuoiOO/AL4K7EdI0luA8cBBwDrCiWpwXHePzHI+BTy1C0l6Edsf/AJqs+9Dofd7R0naYv7KeDz1LVB2SDxmTgfKW5R9q/m4iONHAZtarLfgcbaj1+5ckw4DLohdm3Wx6/MRwsHTbEVmeFOB8d4tlrkkM/wGcABA7KIcAfw5ll1DeENmxK7E51osZ3lmuK55PZIGSfpl7KKsB+4mnI1/RrjebTWGuL3/EuvdS+hJNAIPsmuGAQe02G/XEw7atupNztSZF9c/GMDMZhBaLBGSMGutmb1bYLsGEFrzN1qUDcmMZ/dHa04D/mxmm+P4NMLJ4SRC6/U0oQU/GfijmTXF7SkHlmW26Q5Ci7qzDsjGZ+GIf298B+93QZJKJd0ULynWE5KJQnXMbAGhd/MtYGVczwGZWVoegxXa/pq+teOsVTubpNmfyiwhtKSVmVcvM7tpJ5dVyIGZ4aHA0jh8FqFVbQQws+VmdrmZHUA48/9/7dzHLjfGbTjGzPYBPgvIzJ4htELvxaBwI+TnQImkPwIbgO8QrsnGm1lfM+tpZnfs4jYuARa12G99zGzsTtQ7p0W9CjN7K8Z7JaGLtpRwEsvqJ6lXZrx5364mtMbDWpS9lRlv+fOoQj+XGgs8mhmfRuh2nxKH/wScQEjSaZnt2QIMyGzPPmbW6jVgAcvIHDOSxPbHUMH3ewfb8mlgHKF17EvowtKizrbKZvea2UcI+8+Am3ch9taO9VbtbJKuIFxDQTgrfVzSWfEMVBFvUlTtQqAtfU1SP0kHAlcRblRAvCnRPJOkCzLrWUvYQY07sfw+wEZgnaQhwNcKzHNlXPZP43K/B1xN2O7xhO7qYZL+Em9YnLeL2zgDWB9vOvSI++5oSaPaqHc78B1JwwAkDZQ0Lg4fBtxAOAj/CbhG0nEt6n9bUjdJJwIfAx6MJ70H4nL7xGV/lfDetmYFUCWpW2baOWTeHzObT+gpfRZ4xszWx3p/T0xSM1sGPAF8T9I+CjfGRkg6uY39kPUo8AFJ58dW6suELnaztt7v7PHcPP8WQte9J/Dd1las8Ln4RyV1J1x/bmLnjsFmV0qqUrj5dz3bjvVW7WyS3gj8a+ya/CPhrHM94ebMEsJO2JOPc35DuKnyIuENmBTPjmcQbkw0GwU8J2kj8AhwlZkt2onlf5twA+iduPxfFZjnXmAqoas2gnCCuIPQZb4cOJbQzRkOvE64o1u5sxsYE+PjwHHAIkJr9hPCmXtHfkjY1ickbSDcRBoTD867gZvNbFZMkOuBX8QDCELXay3hbH0PoSfwSiz7EuFG0kJCi3cv4QTVmj8QPmpaHu+YHg1sNLM3W8w3DXg7M30aoUV6ITPPRYTu9ssxvofY/nJph8xsNXABcBMhsQ4l3IRp1tb7/d7xLOlqwr2ONwg9iZcJ+/g9CnfIH4+j3eN6VxP27yDCft9Z9xJOUgvjq+AnAdtp66K1o1+EVuuQAtNHAzOKsP7hhBtHpwP7AMtame924JLM+JPAqLz33w626xTiTasOWv41wC15b2dnetHGDavWXql/meHfi7kyC92zRZIugHCtI+nYWPxrwmeiSBoAHEY4E3ZVi4E78w6iK0g2Sc1shpk93vacu0/SfcBfCHcb75N0GfAZ4DJJswjdu3Fx9t8Bb0t6mfC57tfM7O12iuNxhQ/WW752pRtVVGb2gJnN66jlSzqxlX2ysaPWmSrFZtg5l6hkW1LnXNDqF6f3dgMGDLDhw4fnHYbbBTNnzlxtZgPzjqPYumySDh8+nJqamrzDcLtA0httz7X38e6uc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS5wnqXOJK3qSxscazlZ4ROb7vvITfx52q6QFkl6SNDJTdpWkOQrPN5qQmX6cpOnNy5Q0ukib41yHy+trgada+HV9IecQfml/KOGRJbcRnkRwNOEJCaMJP9KeIulRC08kuAX4tpk9LmlsHD+lg7fBuaJIsbs7DrjLgulApaT9gSMJj3GsM7MGwmM5mp8zZISnKkB4HEmbD3dyrrPIoyU1wvN6DLjDzCa2KB/C9o89rI3T5hAenNWf8PCnsUBzd3kC8DtJ/0U48Xy448J3rrjySNITzGyppEHAVEmvWHi0ZrNCj1E0M5sn6WbCw8I2ArMID38G+ALwFTN7WNIngUmEZxZtR+GfPl0BMHTo0PbbIuc6UNG7u2a2NP5dCUwmXGNm1bL9s0mriN1XM5tkZiPN7CTC83Lnx3kuZtsT4R4ssMzmdU80s2ozqx44sMv9LNF1UkVNUkm9JPVpHgbOJHRjsx4BLop3eY8H3rHwrFZi64ukocD5hP8NAyGJm5/b+lG2Ja9znV6xu7uDCf8yoXnd95rZFEnjAczsdsLDlscS/ulTHXBppv7D8Zq0HrjSzNbG6ZcDP4zPot1M7NI6tzfosg8iq66uNn8yQ+ciaaaZVbc9594lxY9gnHMZnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS5wnqXOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziWuLO8AHGyuf5e7Z1zLvBXP0qd7fy4acwsHDxiZd1guEUVvSSUtljRb0ouSagqUS9KtkhZIeknSyEzZVZLmSJoraUKLel+S9Gosu6UIm9Ju7vzLV3ih9gk2blnDsvXz+cFTn2Ft3bK8w3KJyKslPdXMVrdSdg5waHyNAW4Dxkg6GrgcGA1sBaZIetTM5ks6FRgHHGNmWyQN6vhNaD9zlz9NQ9OW98ZFCQtW1TBq2MdzjMqlIsVr0nHAXRZMByol7Q8cCUw3szozawCmAefFOl8AbjKzLQBmtjKPwHdXRVnv7caNJnp1q8wnGJecPJLUgCckzZR0RYHyIcCSzHhtnDYHOElSf0k9gbHAgXGew4ATJT0naZqkUYVWLOkKSTWSalatWtVuG7SnPjv6RspLKygvraB7WU8O7j+SI/Y7Ie+wXCLy6O6eYGZLY5d0qqRXzOyZTLkK1DEzmyfpZmAqsBGYBTTE8jKgH3A8MAp4QNLBZmYtFjIRmAhQXV29XVmejqs6i+vO/A3zV81gn4oBHDfkLEqUYifH5aHoSWpmS+PflZImE64xs0lay7YWEqAKaK4zCZgEIOm7cd7mOr+KSTlDUhMwAEinuWzDkMojGFJ5RN5huAQV9XQtqZekPs3DwJmEbmzWI8BF8S7v8cA7ZrYs1hkU/w4Fzgfui3V+DXw0lh0GdANauzHlXKdS7JZ0MDBZUvO67zWzKZLGA5jZ7cBjhOvNBUAdcGmm/sOS+gP1wJVmtjZO/ynwU0lzCHd+L27Z1XWus1JXPZarq6utpuZ9H9O6hEmaaWbVecdRbH53wrnEeZI6lzhPUucS51+w70LWbanjq398kL+ufIMDevXl+yd+kqP23T/vsFwbvCXtQi578i6efus11m6pY+6aZfzD43ewdktd3mG5NniSdhH1TY38dcUbbG1qfG+aGcxe/VaOUbmd4UnaRZSphG6l21/dNFoTld175BSR21mepIlqencrTWvqaK/PsSXxnePH0aO0nIrSMnqWdeO0qiP4YP8h7bJ813H8xlFizIwN332SuvtngUT50YPpN/ECSnp12+Nl/+Nh1Ry57368sGoJQ3pXclrVEcRvf7mEeZImZvOUV6n71RxoaAKgfu4KNtz4B/recHa7LP+YAVUcM6CqXZblisO7u4mpn7scNtVvm7C1ka0v+aNUujJP0sSUHdwfemQ6OGUllB86IL+AXO48SRPTY9wH6H7CQdC9DPUsp3RIX/b5xml5h+Vy5NekiVFpCZU/HEfjm+uwLQ2UHbQvKi/NOyyXI0/SBEmibFi/vMNwifDurnOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5x/rXAxDQ0NDBt2jSWL1/OiBEjGDNmjP8wu4vzJE1IU1MTt956KwsXLqS+vp7p06ezePFiLrzwwrxDczny7m5Cli5dyqJFi6ivDz/63rp1K9OmTXtv3HVNnqQJaWhoKNi1bWxsLDC36yo8SRNSVVVF3759KS0Nvx8tLy/nqKOOoqKiIufIXJ78mjQhZWVlXHvttTz44IMsX76cQw45hHPPPTfvsFzOPEkT07t3by699NK2Z3Rdhnd3nUucJ6lziSt6kkpaLGm2pBcl1RQol6RbJS2Q9JKkkZmyqyTNkTRX0oQCda+WZJL8GZhur5HXNempZra6lbJzgEPjawxwGzBG0tHA5cBoYCswRdKjZjYfQNKBwBnAmx0dvHPFlGJ3dxxwlwXTgUpJ+wNHAtPNrM7MGoBpwHmZet8HrgHa5z8cuaJ5Ze1yfr9kHrUb1+YdSpLyaEkNeEKSAXeY2cQW5UOAJZnx2jhtDvAdSf2BTcBYoAZA0ieAt8xsln/PtXP53gtTuW32M5SXlNDQ1MR/n3whZw/7QN5hJSWPlvQEMxtJ6NZeKemkFuWFsszMbB5wMzAVmALMAhok9QS+AXyzrRVLukJSjaSaVatW7dFGuD23aP1qbps9jc2N9Wyo38Kmxnq+9Mwv2+3fPe4tip6kZrY0/l0JTCZcY2bVAgdmxquA5jqTzGykmZ0ErAHmAyOAg4BZkhbH+Z+XtF+BdU80s2ozqx44cGD7bpjbZSvq1lNesn1nbmtjI5sa/LvKWUVNUkm9JPVpHgbOJHRjsx4BLop3eY8H3jGzZbHOoPh3KHA+cJ+ZzTazQWY23MyGE5J8pJktL85Wud11eOVgmjK3EEoQQ/v0o2f5nv8v1r1JsVvSwcCfJM0CZgCPmtkUSeMljY/zPAYsBBYAPwa+mKn/sKSXgd8CV5qZ32noxPpV9OIXZ1zKwB69EeKQykHcc+ZleYeVHHXV/n91dbXV1LzvY1qXkyZrokQ7bjMkzTSz6iKFlIwUP4JxXVBbCdqV+Z5xLnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS5wnqXOJy+ufCO/VGjesZvMbz1PaZyDdhx6H/ztGtyc8SdvZ5kU11N5yBiCssZ4+oy9g8GWTPFHdbvPubjtbdsc/0bRpPU2b3sG21rHhrw+xad5TeYflOjFP0nbWsG7p+6bVv/1mDpG4vYUnaTurOHgMlJZvm2BNVBw8Kr+AXKfnSdrO9h9/DxXDRoJKUPdeDL70x3Qf8oG8w3KdmN84amdl+wxk6Df/jDXUQ2mZ3zBye8yTtIOorLztmZzbCd7ddS5xnqTOJa7oSSppsaTZkl6UVFOgXJJulbRA0kuSRmbKrpI0R9JcSRMy0/9T0itx/smSKouzNc51vLxa0lPN7Dgzqy5Qdg5waHxdAdwGIOlo4HJgNHAs8DFJh8Y6U4GjzewY4DXgug6O37miSbG7Ow64y4LpQKWk/YEjgelmVmdmDcA04DwAM3siTgOYDlTlEbhzHSGPJDXgCUkzJV1RoHwIsCQzXhunzQFOktRfUk9gLHBggfqfAx4vtGJJV0iqkVSzatWqPdoI54olj49gTjCzpZIGAVMlvWJmz2TKC32waGY2T9LNhK7tRmAW0JCdSdI34rR7Cq3YzCYCEwGqq6ttzzfFuY5X9JbUzJbGvyuByYRrzKxatm8hq4DmOpPMbKSZnQSsAeY3zyTpYuBjwGfMzBPQ7TWKmqSSeknq0zwMnEnoxmY9AlwU7/IeD7xjZstinUHx71DgfOC+OH42cC3wCTOrK8rGOFckxe7uDgYmx6/KlQH3mtkUSeMBzOx24DHC9eYCoA64NFP/YUn9gXrgSjNbG6f/COhO6D5DuME0vgjb41yHK2qSmtlCwscnLaffnhk24MpW6p/YyvRD2itG51KT4kcwzrkMT1LnEudJ6lziPEmdS5wnqXOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS5/9VLQdbVz/L5mX/Q0n3gfQ86POUlO+Td0guYZ6kRbZpyUOsq7kEGjdBSXfqXr+NAWe8QElZ77xDc4ny7m6RbZh9TUhQgKYtNG1ewea3JucblEuat6RFZo2btx+naVvSdjGNDVt49U+TeGfFfAYOH8WIUReiEm83WvIkLbIewy/h3QX/DY3hGd5SGd33OzvnqIrPmpp46sef4e3al2hq2MyS2Y+xpvYlRv/9TXmHlhw/bRVZn6NvoPcR11G2zwfpNvCj9D/lGUp7Ds07rKJbt+JV1rw1m6aG0LNorN/Ewpr7aajf3EbNrsdb0iKTSuhz5PX0OfL6vEPJlTU1ILVsI4Q1NeYST8q8JXW5qNzvCHpWHkBJaTkApeUVHHD4qZR375VzZOnxltTloqS0nDO+8DCzptwcbhwdNIqjT5+Qd1hJ8iR1uenWs5JR59+YdxjJ8+6uc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnElf0JJW0WNJsSS9KqilQLkm3Slog6SVJIzNlV0maI2mupAmZ6ftKmippfvzbr0ib41yHy6slPdXMjjOz6gJl5wCHxtcVwG0Ako4GLgdGA8cCH5N0aKzzdeBJMzsUeDKOO7dXSLG7Ow64y4LpQKWk/YEjgelmVmdmDcA04LxMnZ/H4Z8D5xY5Zuc6TB5JasATkmZKuqJA+RBgSWa8Nk6bA5wkqb+knsBY4MA4z2AzWwYQ/w7qsOidK7I8vhZ4gpktlTQImCrpFTN7JlOuAnXMzOZJuhmYCmwEZgENu7LieFK4AmDo0K738zDXORW9JTWzpfHvSmAy4Rozq5ZtLSRAFdBcZ5KZjTSzk4A1wPw4z4rYJSb+XdnKuieaWbWZVQ8cOLC9Nsm5DlXUJJXUS1Kf5mHgTEI3NusR4KJ4l/d44J3mrmxsfZE0FDgfuC9T5+I4fDHwmw7dEOeKqNjd3cHAZEnN677XzKZIGg9gZrcDjxGuNxcAdcClmfoPS+oP1ANXmtnaOP0m4AFJlwFvAhcUY2OcKwaZWd4x5KK6utpqat73Ma1LmKSZrXxst1dL8SMY51yGJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnD8cOwfPLf4Nzy78JT3KevOJY65mSOXheYfkEuZJWmTPvn4/v5z572xt3ASIeSue5d/OfpyBfYblHZpLlHd3i+z3r/4kJiiAUd+4mZo3f5trTC5t3pIWWam23+WihJKSzv82zHt5BU88/ioNDcZHThrOmA95z6C9eEtaZB//4FcpL60AQJRSUd6L44ef10attC1euIa7fzaTt2rXs2L5Bh6Z/DI1M5a0XdHtlM5/Cu9kjq06gy+dfCfPLf41Pcr7cNrhl9G3x+C8w9ojLzz/FvX1Te+N19c3MmP6EqpHH7iDWm5neZLm4PDBH+bwwR/OO4x2072ijJISaNqWp1RU+KHVXry76/bYCSceREVFOaWlQiXQrVspZ431j5Xai5/u3B7r27eCf7n2ZJ6f+RYNDU0cc+z+DBzUO++w9hqepK5d9NmngpNPHZF3GHsl7+46lzhPUucS50nqXOI8SZ1LnCepc4nrsg/HlrQKeKMdFjUAWN0OyymWzhRvy1iHmVmX+yc+XTZJ24ukms70VPXOFG9nirUjeXfXucR5kjqXOE/SPTcx7wB2UWeKtzPF2mH8mtS5xHlL6lziPEmdS1yXTVJJB0p6StI8SXMlXRWn3y/pxfhaLOnFFvWGStoo6erMtL+VNFvSAkm3Kv4rc0nd4/IWSHpO0vBMnYslzY+vi9s7VknHSPpLnH+2pIpixLo78Uoql/TzGNc8SdcVa992CmbWJV/A/sDIONwHeA04qsU83wO+2WLaw8CDwNWZaTOADwECHgfOidO/CNwehy8E7o/D+wIL499+cbhfe8VK+AniS8Cxcbw/UFqMWHcz3k8Dv4zDPYHFwPBixZv6q8u2pGa2zMyej8MbgHnAkObyeMb+JHBfZtq5hDd9bmba/sA+ZvYXC0fJXcC5sXgc8PM4/BBwWlzuWcBUM1tjZmuBqcDZ7RjrmcBLZjYr1nnbzBqLEetuxmtAL0llQA9gK7C+WPGmrssmaVbsKv0N8Fxm8onACjObH+fpBVwLfLtF9SFAbWa8lm0H5BBgCYCZNQDvEFq196YXqLPHsQKHASbpd5Kel3RNHrHuQrwPAe8Cy4A3gf8yszV5xJuiLv9kBkm9CV3YCWa2PlP0KTKtKCE5v29mG+Nl0XuLKLBYa6NsR3XaI9Yy4CPAKKAOeFLSTCBbp0Nj3cV4RwONwAGELuofJf2+jXW3e7yp6tItqaRywkF0j5n9KjO9DDgfuD8z+xjgFkmLgQnA9ZL+D+FMXZWZrwpYGodrgQMzy+wLrMlOL1CnPWKtBaaZ2WozqwMeA0YWK9bdiPfTwBQzqzezlcCzQHUx401a3hfFeb0IZ9y7gB8UKDubcJC3VvdbbH/j6K/A8Wy7uTE2Tr+S7W9uPGDbbm4sIrQa/eLwvu0Va1zm84SbMGXA74G/K0asuxnvtcCdsV4v4GXgmGLFm/or9wBy2/DQHTTCXdAX46v5APgZMH4HdVsmaTUwB3gd+BHbvslVQbgTvIBwl/LgTJ3PxekLgEvbO1bgs4QbXHOAW4oV6+7EC/SO654bE/RrxYw39Zd/LdC5xHXpa1LnOgNPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXufwEJPh7r/KZQhwAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOkAAAEICAYAAAC6WNNwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAeO0lEQVR4nO3dfZxWdZ3/8dd77hiYuBMYUmAgFU0xNRqV1lRMvKOMYtd2a/OukmV1S9vSzGqzuy1t28xfm8QvrTAxNaU0b5IsIUtiweROTFBRJlBAUW5lbvjsH9/vwGG6Bhjkus53Zj7Px+N6zDnne77nfM65zud8v+dc13VGZoZzLl1leQfgnNs9T1LnEudJ6lziPEmdS5wnqXOJ8yR1LnFdJkklmaRDc1jv4ZL+LGmjpE9KmiLpi7FsrKSGUsfUmUm6UNKj+1DvEUkfL0ZMHYjhGkk/3d/L3asklbRC0rj9vfK9JWmqpEl5rX8PrgQeMbPeZnaDmU02s6/mHdT+trvkkfS0pMNKHdO+yvt47qjO0pKeBdyfdxDtGA4syTuIYpJUsZuyQ4AyM3u6hCF1K3tMUkm3AHXAvZI2SbpS0hhJf5T0qqQFksZm5n9E0tdi+SZJ90oaIOlWSRsk/a+kEZn5LXYTn5W0TtK3JJVlyo8GXjWzBkmHSpol6bU47+1twh0naZmk9ZL+R5LiMg6R9FtJL8d6t0rql1nHCkmfk/RkrPsjSdWZ8vdKeiJu7x9jTEj6LXAq8L24rYdJ+rGkr7WzLw+SdJektZKek/TJvdj/ZZKukvRMjP8OSQfEshsl/Twz77WSHlYwVlKDpKvjNq+Q9M+ZeftKmhZjeV7SF1r3e2w1/yDpO5JeAW4HpgDvjNv5aibE9wD3S3pL3D+ty/ihpDWZ9f1U0uWZdd8kabWkv8bjpXxP+6LNfjld0lPxWPgeoExZu+93oeM5Tr9T0otxebMljdrNuj8b494o6S+STssUV8X9ulHSEkn1mXq7Pc7aZWZ7fAErgHFxeAjwMjCekOSnx/FBsfwRYDlwCNAXeBJ4GhgHVADTgB9llm3A74AD4s57Gvh4pvwq4Btx+Dbg83G91cC72iznV0C/uJy1wFmx7NAYZw9gEDAbuL7N9i0GhsU4/gB8LZaNBtYAJwDlwAVx/h6Z7c3G++NM3bFAQxwuA+YD/wFUAQcDzwJn7mHfXw7MAYbG+H8A3BbLesX9dSFwErAOGJpZdzPw37HeKcBm4PBYPg34JdAbGBGX87FYdmGs+4n4nvWM0x4tEN+DrdsAvAC8Iw7/JW7fEZmyt8fhX8TtqAFqgbnAv2TW/TfrabPOgcAG4B+ASuBTMd6Pd+D9HtdmmR+N+6IHcD3wRDvrPhxYCRwUx0cAh8Tha4DXCblRDnwDmLM3x9lut3cfkvSzwC1tyn8NXJA5aD+fKfs28EBm/JzsDiAk11mZ8UuAhzPjvwdOyhxYU4kHYpsYjF2T9g7gqna25/3An4GbCQnYCEzOlI8HnonDNwJfjcMfJJx0tgG/6WCSngC80CaOz5E5YbUT61LgtMz4gUATUBHHjwdeAZ4HPpSZbyzhwK1ps0++GA+gbcCRmbJ/IVxbQ0iUtrFeSJvkIZwkXgaq4/gtwL8DbyYk6XXAZOAtwKuEE9XguO6emeV8CPhdB5L0fHY9+AU0ZN+HQu/37pK0zfz94vHUt0DZofGYGQdUtim7pvW4iONHAlvbrLfgcba7175ckw4Hzo1dm1dj1+ddhIOn1UuZ4a0Fxt/UZpkrM8PPAwcBxC7KW4E/xrIrCW/I3NiV+Gib5byYGd7Suh5JtZJ+FrsoG4CfEs7GPyZc77YbQ9zeT8d60wk9iRbgTjpmOHBQm/12NeGg3VO9GZk6S+P6BwOY2VxCiyVCEmatN7PNBbZrIKE1f75N2ZDMeHZ/tOc04I9m9nocn0U4OZxMaL0eIbTgpwC/N7PtcXsqgdWZbfoBoUXdWwdl47NwxO8Y3837XZCkcknfjJcUGwjJRKE6Zrac0Lu5BlgT13NQZpa2x2C1dr2mb+84a9feJmn2pzIrCS1pv8yrxsy+uZfLKmRYZrgOWBWHzyS0qi0AZvaimV1sZgcRzvzf19597PKNuA1Hm1kf4COAzGw2oRXaEYPCjZCfAGWSfg9sBL5OuCabbGZ9zayXmf2gg9u4EniuzX7rbWbj96Le2W3qVZvZX2O8lxK6aKsIJ7Gs/pJqMuOt+3YdoTUe3qbsr5nxtj+PKvRzqfHAfZnxWYRu99g4/ChwIiFJZ2W2ZxswMLM9fcys3WvAAlaTOWYkiV2PoYLv92625cPABELr2JfQhaVNnZ2Vzaab2bsI+8+AazsQe3vHerv2NklfIlxDQTgrnSPpzHgGqo43KYZ2INC2rpDUX9Iw4DLCjQqINyVaZ5J0bmY96wk7qGUvlt8b2AS8KmkIcEWBeS6Ny745LvfbwGcI2z2Z0F09TNJj8YbFBzq4jXOBDfGmQ8+4746SdNwe6k0Bvi5pOICkQZImxOHDgK8RDsLzgCslHdum/pclVUk6CXgvcGc86d0Rl9s7LvvfCe9te14Chkqqykw7m8z7Y2bLCD2ljwCzzWxDrPf3xCQ1s9XAQ8C3JfVRuDF2iKRT9rAfsu4DRkmaGFupTxK62K329H5nj+fW+bcRuu69gP9sb8UKn4u/W1IPwvXnVvbuGGx1qaShCjf/rmbnsd6uvU3SbwBfiF2TfyScda4m3JxZSdgJb+TjnF8Sbqo8QXgDbopnx9MJNyZaHQf8SdIm4B7gMjN7bi+W/2XCDaDX4vLvLjDPdGAmoat2COEE8QNCl/li4BhCN2cE8Azhjm6/vd3AmBjnAMcCzxFasx8Szty7813Ctj4kaSPhJtIJ8eD8KXCtmS2ICXI1cEs8gCB0vdYTzta3EnoCT8WyTxBuJD1LaPGmE05Q7fkt4aOmF+Md06OATWb2Qpv5ZgEvZ6bPIrRIf87Mcz6hu/1kjO/n7Hq5tFtmtg44F/gmIbFGEm7CtNrT+73jeJb0GcK9jucJPYknCft4B4U75A/E0R5xvesI+7eWsN/31nTCSerZ+Cr4ScAu9nTRWuwXodU6tMD044G5JVj/CMKNo3FAH2B1O/NNAS7MjD8MHJf3/tvNdo0l3rQq0vKvBK7Lezs704s93LBq75X6lxm+VMqVWeiePSfpXAjXOpKOicW/IHwmiqSBwGGEM2F3tQL4Ud5BdAfJJqmZzTWzB/Y8576TdBvwGOFu422SPgb8M/AxSQsI3bsJcfZfAy9LepLwue4VZvbyforjAYUP1tu+OtKNKikzu8PMlhZr+ZJOamefbCrWOlOl2Aw75xKVbEvqnAva/eJ0Vzdw4EAbMWJE3mG4Dpg/f/46MxuUdxyl1m2TdMSIEcybNy/vMFwHSHp+z3N1Pd7ddS5xnqTOJc6T1LnEddtrUpe+PzxzB3cv+AZNLds4bvj78g4nN56kLklPr5nDz+Z/kcaW8Cu4P62YkXNE+fHurkvS0hcfpbFl247xppbXdzN31+YtaTdhZmx58mGa162gx4h6qocfm3dIu9WnehCV5T12JKcK/7SzW/Ak7SZeunkSG+feDvFroLXn/T/6nnRhvkHtxokH/yOPPvMz1m4KH42Wqft2+rrtd3fr6+utu3yZYVvDYl74yjuxxi07pqmymkOnbiL8bDdNzS2NLHlxNo3NWzm8dgx9e9XON7P6PdfsWkrekkpaQXgkSQvQ3Hanxx97f5fwaI4thN9wPh7LLiP8AFvA/zez6+P0Ywm/96wmPHzrEgvP/nFAy8Z1qLxyl2eGWEsT1rQNVe35iZJ5qSiv4pghneYZ1kWTV3f3VAu/ri/kbMIv7UcSHllyI+FJBEcREvR4wo+0H5R0n4UnElwHfNnMHpA0Po6PLfI2dBo96o6FbItZXkmPumMoSzhB3U4pdvQnANMsmAP0k3QgcAThMY5bzKyZ8FiO1ucMGeGpChAeR7LHhzt1J+U1/Rh21W+pGnoUZT370uvwUxjyqV/lHZbbS3m0pEZ4Xo8BPzCzqW3Kh7DrYw8b4rTFhAdnDSA8/Gk80HpReTnwa0n/RTjx/F3xwu+cetQdw4ivLcg7DLcP8kjSE81slaRaYKakpyw8WrNVoTsZZmZLJV1LeFjYJmAB4foT4F+BT5nZXZI+CNxEeGbRLhT+6dMkgLq6uv23Rc4VUcm7u2a2Kv5dA8wgXGNmNbDrs0mHEruvZnaTmY02s5MJz8tdFue5gJ1PhLuzwDJb1z3VzOrNrH7QoG73s0TXSZU0SSXVSOrdOgycQejGZt0DnB8fAjYGeM3Cs1qJrS+S6oCJhP8NAyGJW5/b+m52Jq9znV6pu7uDCf8yoXXd083sQUmTAcxsCuFhy+MJ//RpC3BRpv5d8Zq0CbjUzNbH6RcD343Pon2d2KV1rivwLzO4TkNSt/wyQ4ofwTjnMjxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS5wnqXOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4kreZJKWiFpkaQnJM0rUC5JN0haLmmhpNGZssskLZa0RNLlbep9QtJfYtl1JdgU50qiIqf1nmpm69opOxsYGV8nADcCJ0g6CrgYOB5oBB6UdJ+ZLZN0KjABONrMtkmqLf4mOFcaKXZ3JwDTLJgD9JN0IHAEMMfMtphZMzAL+ECs86/AN81sG4CZrckjcOeKIY8kNeAhSfMlTSpQPgRYmRlviNMWAydLGiCpFzAeGBbnOQw4SdKfJM2SdFyhFUuaJGmepHlr167dbxvkXDHl0d090cxWxS7pTElPmdnsTLkK1DEzWyrpWmAmsAlYADTH8gqgPzAGOA64Q9LBZmZtFjIVmApQX1+/S5lzqSp5S2pmq+LfNcAMwjVmVgM7W0iAoUBrnZvMbLSZnQy8AizL1Lk7dpHnAtuBgcXbCudKp6RJKqlGUu/WYeAMQjc26x7g/HiXdwzwmpmtjnVq4986YCJwW6zzC+DdsewwoApo78aUc51Kqbu7g4EZklrXPd3MHpQ0GcDMpgD3E643lwNbgIsy9e+SNABoAi41s/Vx+s3AzZIWE+78XtC2q+tcZ6XueizX19fbvHl/8zGtS5ik+WZWn3ccpZbiRzCum2le8QqN8xvYvnFb3qEkKa8vMzgHwIavP8yWny9ElWVQJg748T9R+Vb/LkqWt6QuN43zVrLl7kWwrRnb1Iht2MarV/wq77CS40nqctPS8Bpq86n49tUb8gkmYZ6kLjcVR9Ri2zM3LstFxRGD8wsoUZ6kLjeVh9fS5wvjoKocysuoOHgA/f77nLzDSo7fOHK56jXxbfScMAp7vZmymqq8w0mSJ6nLncrLkCdou7y761ziPEmdS5wnqXOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEucPx3Yd1rK9mRfWLwYz6g44ivKyyrxD6tI8SV2HvN60mW/95h9Yu2kFIAbUDOXK0++iZ2XvvEPrskre3ZW0QtIiSU9ImlegXJJukLRc0kJJozNll0laLGmJpMsL1P2MJJM0sMib0W098OT/8OKGZ9jWvIVtzZtZs/E57l10fd5hdWl5taSnmtm6dsrOBkbG1wnAjcAJko4CLgaOBxqBByXdZ2bLACQNA04HXih28N3ZixuW07x9247x5u2NvLThmRwj6vpSvHE0AZhmwRygn6QDgSOAOWa2xcyagVnABzL1vgNcCdjfLNHtN4fXvpOq8p47xqvKe3JY7ZgcI+r68khSAx6SNF/SpALlQ4CVmfGGOG0xcLKkAZJ6AeOBYQCS3gf81cwWFDd0N/awCzh++ATKVE6Zyhk97GxOf+vFeYfVpeXR3T3RzFZJqgVmSnrKzGZnylWgjpnZUknXAjOBTcACoDkm7OeBM/a04nhSmARQV1f3RrejWypTGeedcC0fqv8qhlFZ3iPvkLq8krekZrYq/l0DzCBcY2Y1EFvIaCjQWucmMxttZicDrwDLgEOAtwALJK2I8z8u6c0F1j3VzOrNrH7QoEH7d8O6mYryKk/QEilpkkqqkdS7dZjQ+i1uM9s9wPnxLu8Y4DUzWx3r1Ma/dcBE4DYzW2RmtWY2wsxGEJJ8tJm9WJqtcq64St3dHQzMkNS67ulm9qCkyQBmNgW4n3C9uRzYAlyUqX+XpAFAE3Cpma0vZfDO5aGkSWpmzwLHFJg+JTNswKXt1D9pL9Yx4g2E6FxyUvwIxjmX4UnqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS5wnqXOJ8yR1LnGepM4lzpPUucR5kjqXOE9S5xLnSepc4jxJnUucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEudJ6lziPEmdS1zJk1TSCkmLJD0haV6Bckm6QdJySQsljc6UXSZpsaQlki7PTP+WpKfi/DMk9SvN1jhXfHm1pKea2bFmVl+g7GxgZHxNAm4EkHQUcDFwPHAM8F5JI2OdmcBRZnY08DTwuSLH71zJpNjdnQBMs2AO0E/SgcARwBwz22JmzcAs4AMAZvZQnAYwBxiaR+DOFUMeSWrAQ5LmS5pUoHwIsDIz3hCnLQZOljRAUi9gPDCsQP2PAg8UWrGkSZLmSZq3du3aN7QRzpVKRQ7rPNHMVkmqBWZKesrMZmfKVaCOmdlSSdcSurabgAVAc3YmSZ+P024ttGIzmwpMBaivr7c3uiFmxvSn53LvcwsZ1LM3V4w+g7reB7zRxTq3i5InqZmtin/XSJpBuMbMJmkDu7aQQ4HWOjcBNwFI+s84L3H8AuC9wGlm9oYTcG98b+Ej3LDwt2xtbqJc4ncNf+GRiZ9mYM83lWL1rpsoaXdXUo2k3q3DwBmEbmzWPcD58S7vGOA1M1sd69TGv3XAROC2OH4W8FngfWa2pSQbA9y89A9sbW4CoMWMbdub+c3KpaVavesmSt2SDgZmSGpd93Qze1DSZAAzmwLcT7jeXA5sAS7K1L9L0gCgCbjUzNbH6d8DehC6zxBuME0u9saUaddznBDlSvFenOvMSpqkZvYs4eOTttOnZIYNuLSd+ie1M/3Q/RVjR1x+zLv5ytz72NrSRIXK6F3ZgzPqjswjFNeF5XHjqMs4761jqO3Zm/ufX8zA6jdxydGn0LdHz7zDcl2MJ+kbdObwUZw5fFTeYbguzC+gnEucJ6lzifMkdS5xnqTOJc6T1LnEeZI6lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpc4T1LnEue/Jy0Ba95C4ytzUVkVlQccj8p8t7u950dLkbW8/iLrHn4n1vQqsJ2K3kcyYOzvUHl13qG5TsK7u0W2YcEVbH99Nda8AWveRNNrC9m8/Pt5h+U6EU/SImvZtAysaeeE7a/Tsnl5fgG5TseTtMiqak+D8p0PJ1N5DVWDTskxItfZeJIWWe9RX6L6oPeBykGV9Br5SaqHfjDvsFwn4jeOikxlVfQ/YTp23DRQGfKHZ7sO8iQtEf/Yxe0rP607lzhPUucS50nqXOI8SZ1LnCepc4nzJHUucZ6kziXOk9S5xHmSOpe4kieppBWSFkl6QtK8AuWSdIOk5ZIWShqdKbtM0mJJSyRdnpl+gKSZkpbFv/1LtDnOFV1eLempZnasmdUXKDsbGBlfk4AbASQdBVwMHA8cA7xX0shY5yrgYTMbCTwcx53rElLs7k4AplkwB+gn6UDgCGCOmW0xs2ZgFvCBTJ2fxOGfAO8vcczOFU0eSWrAQ5LmS5pUoHwIsDIz3hCnLQZOljRAUi9gPDAszjPYzFYDxL+1RYveuRLL46cZJ5rZKkm1wExJT5nZ7Ey5CtQxM1sq6VpgJrAJWAA0d2TF8aQwCaCurm7foneuxErekprZqvh3DTCDcI2Z1cDOFhJgKNBa5yYzG21mJwOvAMviPC/FLjHx75p21j3VzOrNrH7QoEH7a5OcK6qSJqmkGkm9W4eBMwjd2Kx7gPPjXd4xwGutXdnY+iKpDpgI3Japc0EcvgD4ZVE3xLkSKnV3dzAwQ1Lruqeb2YOSJgOY2RTgfsL15nJgC3BRpv5dkgYATcClZrY+Tv8mcIekjwEvAOeWYmOcKwWZWd4x5KK+vt7mzfubj2ldwiTNb+djuy4txY9gnHMZnqTOJc6TtAS2bV7P1o1r6K6XFu6N8UfYFdH27S3Muf1TrFx0H0gMGHYsYz96CxVVPfdcuQv489qVNGxaz6gDDuLgvgPzDqfT8iQtomfm3kbDkl+zvSX8m4mXVy5gwYPX8o73XZNvYCXwlbn3cctTcygvK6Nl+3auP+mDvOctb8s7rE7Ju7tF9ErDQlqatu4Y3968jZdXPpFfQCWy/NU1THtqDltbmtjUtI2tLU1c/ugd3t3fR56kRdTvwCMor9z5Lw7LKqrof+CROUZUGi9t3UhlWfku07a1NLOtpUPf4nSRJ2kRjRxzHoMPPZHyih6UV/ak7+DDOHb85/IOq+iO7P9mjJ2tZrnEyL61VFdU5hhV5+XXpEVUVl7ByRfczOb1DWxvaaL3gBGorOufF/tX13DrGR/jkkems2brRo484EB++O7z8g6r0/JvHLlOo7t+48hb0mjRokUsWLCAPn36cNppp1FTU5N3SM4BnqQAzJ49mzvvvJPGxkbKy8t57LHH+NKXvkR1dfWeKztXZF3/Amkv3HvvvTQ2NgLQ0tLC5s2bWbhwYc5RORd4kgLbt2/fZdzMaGlpySka53bl3V3glFNOYebMmTQ2NiKJiooKRo0alXdYXUJzcwv3/uJJnlzyEjU1VUw8923UDfcnrnaEJylwzjnn0KtXLx5//HH69u3LxIkT6dOnT95hdQl337mYBX/+K01N23nt1deZ+v05fPqqsfTv3z2+v7w/eJICkhg3bhzjxo3LO5QuZ/Gi1TQ17Xo58cyyddQfP6ydGq4tvyZ1RdWjqk07IKiu9rahIzxJXVG9/++PorKyjPJyUVVVzuDBvTli1OC8w+pU/JTmimrU297MJZ88kWeWv0xNTRXHvP0gysu9begIT1JXdEOG9mXI0L55h9Fp+SnNucR5kjqXOE9S5xLXbX+qJmkt8Px+WNRAYN1+WE6pdKZ428Y63My63T/x6bZJur9ImteZfuPYmeLtTLEWk3d3nUucJ6lzifMkfeOm5h1AB3WmeDtTrEXj16TOJc5bUucS50nqXOK6bZJKGibpd5KWSloi6bI4/XZJT8TXCklPtKlXJ2mTpM9kpr1D0iJJyyXdoPivzCX1iMtbLulPkkZk6lwgaVl8XbC/Y5V0tKTH4vyLJFWXItZ9iVdSpaSfxLiWSvpcZllFjzd5ZtYtX8CBwOg43Bt4GjiyzTzfBv6jzbS7gDuBz2SmzQXeCQh4ADg7Tr8EmBKH/wm4PQ4fADwb//aPw/33V6yEH04sBI6J4wOA8lLEuo/xfhj4WRzuBawARpQq3tRf3bYlNbPVZvZ4HN4ILAWGtJbHM/YHgdsy095PeNOXZKYdCPQxs8csHCXTgPfH4gnAT+Lwz4HT4nLPBGaa2Stmth6YCZy1H2M9A1hoZgtinZfNrKUUse5jvAbUSKoAegKNwIZSxZu6bpukWbGr9HbgT5nJJwEvmdmyOE8N8Fngy22qDwEaMuMN7DwghwArAcysGXiN0KrtmF6gzhuOFTgMMEm/lvS4pCvziLUD8f4c2AysBl4A/svMXskj3hR1+9+TSnoToQt7uZltyBR9iEwrSkjO75jZpnhZtGMRBRZreyjbXZ39EWsF8C7gOGAL8LCk+UC2TlFj7WC8xwMtwEGELurvJf1mD+ve7/Gmqlu3pJIqCQfRrWZ2d2Z6BTARuD0z+wnAdZJWAJcDV0v6N8KZemhmvqHAqjjcAAzLLLMv8Ep2eoE6+yPWBmCWma0zsy3A/cDoUsW6D/F+GHjQzJrMbA3wB6C+lPEmLe+L4rxehDPuNOD6AmVnEQ7y9upew643jv4XGMPOmxvj4/RL2fXmxh228+bGc4RWo38cPmB/xRqX+TjhJkwF8BvgPaWIdR/j/Szwo1ivBngSOLpU8ab+yj2A3DY8dAeNcBf0ifhqPQB+DEzeTd22SVoPLAaeAb7Hzm9yVRPuBC8n3KU8OFPno3H6cuCi/R0r8BHCDa7FwHWlinVf4gXeFNe9JCboFaWMN/WXfy3QucR162tS5zoDT1LnEudJ6lziPEmdS5wnqXOJ8yR1LnGepM4l7v8A2nE+naJ8wKEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -1058,14 +1054,14 @@ { "data": { "text/plain": [ - "array([ 0.75141507, 1.14833122, 0.96556452, 0.68202242, 0.47354219,\n", - " 2.53208661, 2.88199067, 0.94960797, -0.38487508, 0.28049792,\n", - " 0.53176173, 1.43783564, -0.86276298, 1.06649661, 0.17319538,\n", - " 1.41823675, 0.71377984, 0.48296745, 1.27987357, 2.82573732,\n", - " 2.43129513, 1.94906869, 0.97500437, 1.86407191, 0.58033449,\n", - " 1.11092039, 0.67836559, 1.86015599, 0.95727419, 3.29341463,\n", - " -1.29286744, 0.76073929, 1.39946125, 1.36837727, -0.42501107,\n", - " 3.39964208, 1.04062549, -0.84626529, 1.53551352, 2.33455969])" + "array([ 3.05828465, -0.28932373, 0.4753333 , -0.64006398, 1.33554241,\n", + " 0.67785004, 1.64601032, -2.06629308, 2.33285487, 2.17162016,\n", + " 0.72198926, 0.23620621, -0.94067248, 0.94042514, -0.17516254,\n", + " -0.07060064, -1.84321572, 0.77205522, 2.04546285, 0.32892863,\n", + " 1.63786462, 0.37460045, 0.72719062, 1.971905 , 1.88292559,\n", + " 2.34865765, 1.99634041, 0.46903137, 1.75419429, 1.52706647,\n", + " 1.20447388, 0.87579679, 0.42315797, 1.48726631, 0.49144473,\n", + " 0.04416208, 2.91365807, 2.27138945, 2.75618211, 3.451797 ])" ] }, "execution_count": 28, @@ -1157,27 +1153,23 @@ "execution_count": 32, "metadata": {}, "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "wrote temp/shapefile_export/freyberg.shp\n" - ] - }, { "name": "stderr", "output_type": "stream", "text": [ - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n" ] }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "wrote temp/shapefile_export/freyberg.shp\n" + ] + }, { "data": { "text/plain": [ @@ -1278,9 +1270,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.10" + "version": "3.8.11" } }, "nbformat": 4, - "nbformat_minor": 2 + "nbformat_minor": 4 } diff --git a/examples/Notebooks/flopy3_working_stack_demo.ipynb b/examples/Notebooks/flopy3_working_stack_demo.ipynb index fbbe22e029..38ef50a088 100644 --- a/examples/Notebooks/flopy3_working_stack_demo.ipynb +++ b/examples/Notebooks/flopy3_working_stack_demo.ipynb @@ -20,12 +20,12 @@ "name": "stdout", "output_type": "stream", "text": [ - "3.8.10 (default, May 19 2021, 11:01:55) \n", + "3.8.11 (default, Aug 6 2021, 08:56:27) \n", "[Clang 10.0.0 ]\n", "numpy version: 1.19.2\n", "matplotlib version: 3.4.2\n", "pandas version: 1.0.5\n", - "flopy version: 3.3.4\n" + "flopy version: 3.3.5\n" ] } ], @@ -71,7 +71,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/mbase.py:354: DeprecationWarning: xul/yul have been deprecated. Use xll/yll instead.\n", + "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/mbase.py:352: DeprecationWarning: xul/yul have been deprecated. Use xll/yll instead.\n", " warnings.warn(\n" ] } @@ -388,42 +388,6 @@ "needs_background": "light" }, "output_type": "display_data" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAO8AAAEICAYAAAC3RqM3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAA+c0lEQVR4nO29eZxcZZn3/b2qqru6053uTtKdkIUkJGHHECAiAsoOYXlAXp4RfVV8XYZhZvSRUXRk8FFEHUBH8JmPC6KjMOMyMuMgDAqIKAKPLCaQBJSwh4RAtu50eu+u5Xr/OKe665w61XXq3Kdr6T7fz+f61Fnuc9fdXfWre7vu6xZVJSIiov6IVbsAERERwYjEGxFRp0TijYioUyLxRkTUKZF4IyLqlEi8ERF1SiTeIojIchFREUkUuf8PIvJ9H/ncJiJfDr+E5th/36pqlyMiGJF4bURkq4ic6Te9qv6jqn50Kss0UxCRd4vIH0RkSEQeqnZ56gXPWiUiYqoQkbiqZlyXe4BvAIcBp1e8UHVKVPMCIvJvwFLgv0VkQEQ+k3f7fSKyTUT2isg1ec9cKyI/yjs/2a49ekVku4j8fx7vM1tEfici/ywWt4nIt0TklyLSLyJPiMjKvPSHicgDItIjIs+LyLvz7p0nIn+2n9shIlfZ1ztF5B67HD0i8oiIlPycReR8EXlaRPrs8l+bd++XIvJxV/rNIvIuH+W8TUS+IyK/EpFB4DT3e6vqb1T1DuCNUuWMyENVI7NcRLcCZ+adLwcU+B7QDBwNjAKH2/evBX5kHy8F+oH3Ag3APGCNfe824Mv2tSeBL+e9x21Ytc7xWK2gHwP/bt9rAbYDH7LvHQvsBY60778JvMM+ngMcax9fD9xil6MBeAcgRf5mBVbZx6cCb8H6QV8N7ALeZd97N/BE3nNHA91Ao49y3gbsB06y826a5DP4KPBQtb8L9WJ1V/OKyA9EZLeIPOsz/bvtGupPIvKTAG/5RVUdVtVNwCasL66b9wG/UdWfqmpKVbtVdWPe/UXA74H/UNXPuZ79L1V9UlXTWOJdY1+/ANiqqj9U1bSqPgX8HPif9v0UcISItKnqPvt+7vpCYJldlkfUVsZkqOpDqvqMqmZVdTPwU+AU+/ZdwMEicrB9/gHgZ6o65qOcAHep6v+18x4pVZYIf9SdeLF+ydf5SWh/2a4GTlLVI4ErA7zfzrzjIaDVI82BwMuT5HE+Vu19Sxn5LwPeZjd/e0WkF+tH4gD7/iXAecBrIvJ7EXm7ff1rwEvAr0XkFRH57CTlGkdE3mY36feIyH7gCqATQFVHgTuA99tN8PcC/+aznGDVzBEhU3fiVdWHsZqa44jIShG5T0Q22H28w+xbfwl8S1X32c/unixrg2JtB1ZOcv97wH3Ar0SkpYw8f6+qHXnWqqp/DaCqf1TVi4D5wC+wxIWq9qvqp1R1BfA/gE+KyBk+3u8nwN3AgarajvVDI3n3b8cS5RnAkKo+5qecNtHStSmg7sRbhFuBj6vqccBVwLft64cAh4jI/xWRx0Vkshp7F7Ai4Pv/GDjTbqInRGSeiKxxpfkY8Dxwj4g0+8jzHqyyf0BEGmx7q4gcLiKNIvI+EWlX1RTQB2QAROQCEVklIpJ33T2668VsoEdVR0TkeOD/zb9pizULfJ2JWnfScvp4T+wyx0WkCavPHBORJhFp8Pv8TKXuxSsircCJwH+IyEbgu1h9PrC+DAdjDca8F/i+iHQUyep64HN20++qcsqgqtuwmrCfwmoVbMTVN7b7nZdj1VR32V/WyfLsB84G3oM1CrsTuBFI2kk+AGwVkT6sJu777esHA78BBoDHgG+r6kM+/oy/Aa4TkX7g89g1uYt/xRrUGh9l91FOP3wAGAa+gzXANozVWomYBPExllFziMhy4B5VPUpE2oDnVXWhR7pbgMdV9Tb7/EHgs6r6x0qWd7ogIpcBl6vqydUuS8Q0qHlVtQ94VUT+AsCeP83Ver/AnlcUkU6sZvQr1ShnvSMis7Bq51urXZYIi7oTr4j8FKs5eKiIvC4iH8EaSPmIiGwC/gRcZCe/H+gWkT8DvwM+rard1Sh3PSMi5wB7sMYFgky3RUwBddlsjoiIqMOaNyIiwqKuFiZ0dnbq8uXLq12MGcWft+1ynB+xdEGVShIuGzZs2KuqXdUuhwl1Jd7ly5ezfv36ahdjRrDmYzcDcMgJzuvrv/l3VShN+IjIa9Uugyl1Jd6IqeHUdTcWXlzVWPmCRJRFJN4Zxonv/nrBtUim9Ukk3mnO0Vfe7Dj361gdUftE4p2GuAUblMRI4TTi2o/exPrvfzKU/CPMiMQ7jQgq2rG2BI196fHz0bkJmnuypGbJJE9FVJtIvHXKsX/lIVQfa5UGF8VoeSPreW90bvR1qCeiT6sOOPoThUKNG+Y5Mi/mOm8kud9b1BG1SUnx2kvXHsZa4pUA/lNVv+BKcxHwJaz1nmngSlV91L63FSu+UwZIq+pa17NXYUV/6FLVvaZ/0HRh9ScnBGvaeG3sc/ZdU63+cmwY0qjpXMP4qXlHgdNVdcBeIP2oiNyrqo/npXkQuFtVVURWY60FPSzv/mlewhSRA4GzgG3B/4TpRb5oTeh4eaIPO9RlVk+7BXz0/7qZTf88PZw16pmS4rUXkQ/Yp7mIhOpKM5B32uK+Pwk3A5/BCnA2I1n70ZucF9rMarp80ZogdjaNfVHtW6v46vOKSBzYAKzCign1hEeai7GiUczHCriWQ7GCoSnwXVW91U5/IbBDVTdZEVuKvvflWBEoWLp0qZ/i1ixeg0x+VoZoHMQjkE1iVF3nAQuWh/jUflT7Vh9f4lUrwv0aO4TMnSJylKo+60pzp33vnVj939zWISep6hsiMh94QES2AOuBa7DCp5R671uxF4CvXbu2btYvvsWj+esnKFNjnzJWpPadtds5oJRJmtWITT3RAFU9U9Zos6r2irWXzDrAM26yqj5sR3PsVNW9qvqGfX23iNyJFWB8H3AQkKt1lwBPicjxqrrTK9964ODrJwQ7aYAqH8RDjm7c1FNYdWeS0YrQesbPaHMXkLKF24xVo97oSrMKeNkesDoWy1222w5zGlPVfvv4bOA6VX0Gq3mde34rsLZeR5vzRWtCY7+Ovw7PCy6sfQcXfqxNPR4JyyDm0Zw+5m9u5ulvR03nauGn5l0I3G73e2PAHap6j4hcAaCqt2AFAL9MRFJYkf8utYW8AKspnXuvn6jqfVPxh1SSZd//muO80XC6PCfacomPqqPpPDxXGJ6bMJ5aGpnr/OGIj4FGlXTN4We0eTNwjMf1W/KOb8RVG9vXX8F7exB3uuWl0lSLIz7nUasuD55fc3dh81XjZnIbnms+GixR97fuiDysPLj0sSvyzg4NnE/DsOtCCLWXe9Q5MaTgQ7z7VyRof6Ww7RuJtn6JxJuHU7T+GetK07hn4l+ZWjRGahHMfdxspaxbqKNtMZr2lVabUtwrq+dw50fe8VKk3nplxor3/If/V8G1VsMNNlKLxswyAKOBKjeLHigc/+s5vDNwfsm+wr75Ce+/icd/FC0RrAYzQrzXPXuhx9XlRnnGhp0iS7ea1WCmoj3gicFAz/WuinnWvrP2FPbN083RqFUtMS3Fu/V1984nbzXKb/CgVKDnek4Y82w6j7aHI4LW7UP2K2jMcNArxBo/ojJMG/G++Pqi8WOT1u9Ayim2g85+lWdfXGKQIwwuyj+OM3urmaNYTrSmHPDYxIja8Dw/GxdG1BJ1L9580Qald7R5/LUj6R4iLo+mrU7xDx1QJKEL95xtjvl/KOy3ZlvN/LfyRRuU0TwXzshZozrUlXh7Rp/jJy+9zXHtrT6+x5d1/JF/7S1sOudEa4p2BGtWe9G1vj+0vAA06WyHLHlgP+nW0qPgQ11xz37vqOGqp4jwqCvxmvL7PatCzS+oaPuXi6PpPOe5cJrB8X6XQ/Qkq7X8MjLHX184qn0rz7QV7/6s80+7qO1pI/E29BT+q0xXzoYiWnujOMmYL7jq3DRQcG1kThQstlaZVuLdMjYxh7k0sa9k+qZEipF04fBWa4NzvvaEI15hw6OHGJUt6VoYUGzZn5tsc5LYsMdCXcPdHRP9zjy7/jiKNpb+OsTHlExj1HSuBaaFePNFG4RZeWIdy5qGdoPY0fsLL/6+3ShPTUx8VDKSRpPBy5kYMHcm8Yrp/NYP38QffxA5bFSKuhPvebPecF2JsyVAPn0jEyNdswzmlhq3JR3n2UQILsyuDLItSWTUI5RGhckfwBqeZ/4jF2FGXc3Mz00ebpxH30iTQ7h+Oe7kFxzn2cUjxOIh+AWL07TB30dSTMyZtmaHZZsMfT6xROs18hxRXepKvOUSE3XYN1bd4eu5nkmmkLKLR8guNgtz0bYtQ9u2DMMLTeNtYI0o25ZtSZZOXyq7VMZh8zb2mZcxYkqou2ZzMY5LFjozdGfNfpsOaXVuLH3Iul38+JngrpZt28KpvUz6u/mMdc4quNbYbe7AEVEZ6q7m7Vi03XH+h9F24iEsSh1ONzhsx2Cbr+fSB3l/2Yd2zHZY3yrzMmoyHki47qZzZlYDmRbzjT1TLVJgaz4eTkigiNLUbc37h1Gz0VuAo7smBr+29M6fJGVp3nbmnxznj+ww659nG10ibYwTGytdc0smi8a9f5MzBiNzjb2pgvO9qyN/6GpSl+LtWLQdXj0q8PP5ojUh8Up5X97XTxeW/LZwikUTscJzw3ncrGvgK9uQRDLl1f4jC6xmtaTrJuLujKIuxevFxrG5rGksDJG4IuGswR5a82M+seM0o/cqV7Rutp+XL6JGDv5XM1+tTMA+cHw4Taa58CuQE21EbTNtxJsjv77pimUx3U/v8e7ljvOVi/bw2iuld254x9rneGT9RNP5uGNetg6Oge3PLfd+KAiGNTRA495wfKsjKktdi3dezBk9Yns6yTKDPT+OmbvDcd6bMq+BxkVrgr3AIJtMEBs1q6XdPtCJAX+LK1KtCRoGCt/ba9eFE973dR7/8aeCFTDCN3U32pzjvIM8N2wIzGEdu43z+FP3AQ6be1Bp/2qAFy8r8huam8M1JD6YIj6YIjZiPlWlMmHznommlapJXde8fujPppgdmxhlvWvwAE7veI7f9pqNBp9z3nrTohWg+ToVQQxbxPHBYEsWm3YOMnLAxGqi5h0Tq41SHSE4lkSEwrQXb467Bn2GtJiEJS291us5vcZ5Ne52Ttu8tq6BpfcHr8kSHkJVw1o7X7QRtce0FO/2jMsBwWdrcTjbQHOsUAQ50ZoylrEGz1Yfsp0tu1eUTK+CZ+071lHoBukl3nLIuiJuNO7zt/KoZccog4vN3TIjyqfuxRt0gGqypvOzvea1dD450QYmb0S5wHkjBLJNDf63Qy9C24vOwcOzj7+OXz/5ebNMIyalrsX7tmWvsnOHWQC61wbmhFQab17tncvi2R7re31SrPb1S2ZW4UccS5m5arp/QJr3mMYUiQhC3Y42B6E/0+Swt85+1ddzW/vnel7PIA474YCtIZa2PEbnOQeSMs0Jsk3mtbSKOCyx339L5+zjrzN+/4ji1HXN64c4Vi1zWOOb/HH4IKO8Hn75YMf5OYf82Sg/d426/exmDvx16UGrxFCatEeNCnh6TJWLJqyBrnRbkni/edSNiKmhZM0rIk0i8qSIbBKRP4nIFz3SXCQim0Vko4isF5GT8+5tFZFncvfyrn9NRLbYz90pIh2h/VVYos0J14T5zYPMbx5kKG02KLOgud9hp5y1ybhsjT2jDjMdXdaEjAvXhExrw7idcfr1xvlFeOPnZ3oUOF1VB0SkAXhURO5V1cfz0jwI3G1vqL0auAM4LO/+aR673j8AXK2qaRG5Ebga+Psgf8S/9B5bcO3c2ZuDZDXO/OZge//kc+bi58ePe1JmURjdqx69vJ28n1PPrVDig4U1arrN7AcqY7pTW0RZ+NlcW4HchF+DbepKkz8h2OK+XyTfX+edPg78z1LPeHHA4jeg94Igj47z4vYFBdfmryot3vtfOMKz6ZwvWhPym8AN/eEFds82xAw9vqGhN7gbakQ4+OogiUgc2ACsAr6lqk94pLkYuB6YD5yfd0uBX4uIAt9V1Vs93uLDwM+KvPflwOUAS5eWXhBQimcHnfsOHda2mxcpFG857E9NrDI6fN6uSVL6w7TfmthfGKYnNSf4SihtcEo9Em5t4OtboqoZYI3dL71TRI5S1Wddae60770T+BJwpn3rJFV9Q0TmAw+IyBZVfTj3nIhcgxW//MdF3vtW4FaAtWvXBpo0+fjz7x0/Pnbe60GyGOedK190nA9n/DUV5zYMejadOzc528PDC5I09pXvg5xbE2y69tY9jaRNCQghoHtE+JT1E6+qvSLyELAO8FwZoKoPi8hKEelU1b2q+oZ9fbeI3AkcDzwMICIfBC4AzrCb54G45sh7+MqfCpvO+aINwqLmXsf5hYf10psqXYPtGOlgcVOv572Hf320UZncuBfy+yXV0URD70QNrYkY8aFUQS1bLpIuHCQ8851f4TcPX2OUb0QhJcUrIl1AyhZuM1aNeqMrzSrgZXvA6ligEegWkRYgpqr99vHZwHX2M+uwBqhOUdVQFpRmdOKLfE/fmkB5LFzQC8C2/o4C8ZbLg2+Y7bLgJjZcOMgUVLxhPQ9WxEmAxL4MmdnmsbEi/OGn5l0I3G73e2PAHap6j4hcAaCqtwCXAJeJSAoYBi61hbwAqymde6+fqOp9dr7fBJJYTWmAx1X1iqB/yDVH3sN1z15YMt1T3Us8m8450ZrykxO+7zg/9lf/ECifhv3hz6/Gxpy1YmZWI7ExH6PWcfFsOudEO1M557RZute1nvmpzaP3q+q6Sry/n9HmzcAxHtdvyTu+EVdtbF9/BfBsJ6pquFv2eXDOwue4/01v/+WtgxNeU29ZvoO9w2ZTOW7R+mXv0TFHv3dkToyROTHmbTYTb3aWswaMj2Z8zQNLKmPcdJ4p7OnJ8Mh9Tj/41kXbzPbeKYNp72EFsLzNGduqxzBCxm3H32b0vBd+t9L0Q6YliRiGx8k0FX41gq4Pnq4oSkpD2DUjINNWvJ/r3DL++oHXTgmczzeO8ZzBMmJsnrO5+cbpMPfp0s9lkw3ERqdOQOVujaINce+ms0cNf+Yp/8hvfh+sC1GrKJAKwYsvKNNOvDnRlktn86Cj6XzmwnAcLZ467x89ry+75WuTPte9uoV5m70dRTKtTk8oMV0lZLifUarTOQLv7ltPV7IoI1WseafVqqLPH3W3cR5nLnw+NOGGgmsjMlOhxUZTDgsafVIyGSSTKesbdOYp3j9k9YoqpFxWSaZdzVsOL3fPq3YRHDQOOD/9/hWzmP1q6Vk0bYh51r5ePs2mSGZmjzDnowgpDf9/7JdpVfNOxvaBDoftH/MXSO1zm981tQWLKcSUnuPCXdBuusIILG+rfFOfUTyyjTPja6XAiMYdVkmm5X/5xdSAwz6/6JfVLlIhtmhNSbUkSLUkGOswc45I7O5zWOOOXuOyxQfHCuyctdca51srZBHGiDuskkw78cYOeKF0ohJs3r/YYRc+8rchlMzJa5d/xuj5nGhNiO0bJLZvkMQb/uJLT0Z8OOMwd3D36YgijGQTDqskM7rPm88FB1m7/I1W+APwQ2yocHrIl3AVa6DLnd8+87XKMDF4Fh9KhxIcvt6wpoqq59BSe9/UCpMTba0gWWeNNbCsmbbnzOZ2E7uCB8DzwnTEe7qgKqQq3M/NZ9o1mwEOPdB7C8+xdMJhbc21sS41NhQnNhRn91vN80q80e20Hd3+Hkx7jyJnknGHpTqCrQvOtDSO23SZMsoijGijwyrJtK55D25oHT9+6Mi7OHHTJVUsTSGxIbNf7VlbAu6vFBPIFumT5ok48cY+MgeZuepmWqbvKiNFGItq3vApVvv6fn72Lod9emOgKD2T8uqV5e+kFxsYJjYwzKwX3SHBApCIO800u+4Bh8UGa6NlM1UokNKEwyrJtK55/fLwzsIFTofONg9nE4RdbxMWPFFYK8YGQtqRz+FkEVywXgHsZhpZFUay1ev/T9ua14s/HP1zx/mhc3aza7StSqWxKFb7Zhpk3Patbjd/o0zGJVx/JF911vDpljij86K9iSDnYRV3WCWZUeLNceic3Rw6x3w/3rBo6BOH9Rxu+LFksk4bNa8l0y1x0i1mX85sIjZup59xg3GZqk0k3inmgeGEw754yF3VLlIBL3z+k8EfjonTEiH0hFx5Jl/zN2KtSe8mZE6w0w1FGNEGh1WS6fcfzcOvt1WxWnj78ByHveexvwyzeGXRfaJr58JMSEvRBoecNupvkGlo6eyi99JzW8bNKyDddMEasKpezRsNWOWxsmOif7djpKN6BZmMsEQLkDJf2J91R9zwiMAxXVEVUlX0yJvWNe9kvNY/x2FNCX9f5Jqqff0Qj1lmGBYHrJjQ+Ta8qLX0Q0B8ZHpuAZpFGM0mHFZJZpR4r9x4KVduvJQnd5jvvBA26VnqsO7VPh8s5myRE60po2MwOkbL02bB6sEKgpdvZ538FfPyVZEgA1YicqCI/E5EnrM37vuEff1aEdlhb8i3UUTOK5XXjBBvTrTlEnPtwfnCvi5e2NcVOJzrZLxyVfkOG0VJhuDVZIvWGBEQIT46/RbxW5E0yu7zpoFPqerhwAnA34rIEfa9m1V1jW2/KpXRtBfvOQeZ7aELE6KtGYaHnTboM2Z9saZzMum0YrV5OdiiLUU9174aoNmsqm+q6lP2cT/wHLA4yPtPe/GWS8/oLIf1jgXfoCtMUq1CqlXYeW6gz9mJxCYsae5wkZnV4LBscmbEfVaEdDbuMKDT3qM6Z5cXe15ElmPFRM9t3Pcxe7/qH4jInFLvP3OGBj14Y197wfnK+SH4DIdIqjXgOtlEAtJ5A0WxcASV7CncgXCmouq5/nuvqq4t9ayItAI/B65U1T4R+Q7WBn1qv34da/fMosy4mveQrj0c0rWH9mQNfgnVad1vCSHPWDw04c56ZodxHvG+4QJbt/pzIZSu8ihCWmMO84O9Sf3PgR+r6n8BqOouVc2oahb4HtaGfJMyY8SbE60p3dvm0L1tTsm4y0F45dOGg1a5fmbOGgJ6/AwPO2vtgMT7RwpsOqHg1WyeFLE25voX4DlVvSnv+sK8ZBdTZBfOfGZEs/m5i7/ARY9+zCiP7m0luyCVxS3MsZT5KPOw+cql+MDM2RJFVRjzIVgXJwEfAJ4RkY32tX8A3isia7B+E7YCf1Uqoxkh3nJ54c35VXvvVz79KVZ87euOa3O2OF0WR484kOSLO83eKATvqnxmkmhzWDVveY1XVX0Uz8hilJwacjNjxHvXyd/0rH0PnOeMnDiWqZ2R0kyzsu8IZU6wHVwmyAacY02nPRc6xHf3Oc+BTGd1l1ZWg1yft1rMGPHmWNDUP3G8vJ+X+2tr1wSwRGvMSF5tbepplfd8fOc+4wGwkUWFixpOXXcjD93390b5VhrV6v7Yl/xURaRJRJ4UkU22O9cXPdJcZM9PbbTntk7Ou7dVRJ7J3cu7PldEHhCRF+3XKe9U3nXyN42eT/TGHbbyhptKPxSA1z52VfCHR0adwg1KU9Iyw/A4MjTqtOHpE4FDETLZmMMqiZ93GwVOV9WjgTXAOhE5wZXmQeBoVV2DNTfl3mn6NNvlK3/+67PAg6p6sP38ZwOUf0pYPt+5fjXT1xjK7gZBeeX/8TlqHIY7Y060ZRLf62xKSyqDxqb3ZIYqgaaKwqLku6nFgH3aYJu60gyojvvetbjvF+Ei4Hb7+HbgXX4KHDYrZxdfaJ7pa7SEWyPsX5kct93rloX/Bn4X8k/Sh5ZUxnvP3mmJkM7EHVZJfP1UiEjcHtbeDTygqk94pLlYRLYAv8TpGaLAr0Vkg8tVbIGqvgmWvyfgOcQrIpfnXM327DGfp/WLu/atNM07Yw5780SftWExf+JU2mnZENYFZzMOi+8uf9uUxp7RcTv7+OvMy1RBVCGTFYdVEl/itT0/1gBLgONF5CiPNHeq6mFYNeiX8m6dpKrHAudiraB4ZzkFVNVbVXWtqq7t6gpnccCGPUsc1jviz385Pbdy61K3fNkgNA44Q7qGsRVJzlOrTI+tbKv3/zYn2HpGgYzGHFZJyno3Ve0FHgLWTZLmYWCliHTa52/Yr7uBO5lw+9qV8yqxXysSEe7WtbeXTlQGUzVo5YeCprNqKIvux/MJKz8g3dU2bg19xeeE66v2FdKZmMMqiZ/R5i4R6bCPm4EzgS2uNKtsty9E5FigEegWkRYRmW1fbwHOZsLt627gg/bxB4HaiwwHkBWHpTtqsD8XhshCCrwOoMm4wzLz/EXcqDesZnP1Rpv9jFAsBG4XkTiW2O9Q1XtE5AoAVb0FuAS4TERSwDBwqaqqiCwA7rR1nQB+oqr32fneANwhIh8BtgF/EeYfZsIR7VbA9SOO28V///GYKpfGiXu9967zl7Hgnq2lH2xssFwo3YQk2FzTPNbbT2ZBRzh51jiKkKlwbZtPSfGq6masNYfu67fkHd8I3OiR5hXg6CL5dgNnlFPYqSYn2nJZecNNvPxZwz6qB0MHpVznMO+PhuFFw56+Cdif1kTMEVlS8gLCn3PcF7h/Q4E7Qe2h1q4J1WLGeVjlWNVR3dFkP2z7yN+z9F8KfhP941XTmhDUzXISJMAuDrWEVniEOZ/pPYtehA3nmodeSc/OOGzZt8JfIlgWyUanxX02h4ts7akts5zWHELEjZaEw9LttTOHHgRVyGZiDqskM1K85ZLoizmsmt5W+UjGst3nLg854/BrEwljXrkG0azTKsmMbTYX41e/Pc5xXjtrjCaQsFqaOZFmsuGEibWJ7+qdIYNWEjWbq8HP3n6L43xbfwfr9x5YpdIUJ96fcFjvYZWL7OhGW2Y5L8TiVhA7Q7JNDQV21klfNs53ylHQbMxhlWTGijfHtv4OtvV3mGWiAios++Y/hVKmfIJswO0gFnOaT2RokqgaIcTFatwzROOeISRTG12QwGRdVkFmtHjdtW/Z2KKtGdybhoXQz5SxlMMCx8WyyYl2WqDWaHO+VZIZLV4v5q/29tJseV0c1v585XrDk9W+DYNZGgaz7HtbgH2M3Lj29ZX+QX/PFWk6x3f3OSz5cm2F1Q0DyYjDSqYvvt1J2evbI/EWQRsmLFuDMxo50RozOmrZZM1kv6Qz4xbfUfvz6MaoFLjP+qDYdidlr2+PxJtHa+MorY2jrFi7vdpFKSDblB237uNCGG7OiTYouY3Mcma4qbdksgV29gl1sEihzD7vJNudlL2+PRIvE6INxOy0w5bdfkO4hbN57YpPl0yz7/QV/jILuoghF2LHMPKkDAw7LLFrv1F+VUPLbzbn49ruxNf69nxmvHgfO9uf2PpXOL/wDYNK59PVG6za+c5JBDi7dcKam8zfLGhcLLfIsxmkt987bb3i2uUCn3sVubc7CfLWkZNGmTQM1tbURvNeZ4CAkdVLaXq1J3iG7lo5jGB2U+ATXStIYT+35F5FXtudYK9vV9U3/a5vn/E172Rk56Qctv84fwHepqrpnKNxT5zGPXF6Djf8+NrbnBYg8FwB7oX8Y/7+Z/G9dVgj64SLas5KUWy7EwKsb4/Ea5NVcdhBq/3tBN/5sOESvTLIibZs3E3nhgZoaTEvUG4D7pz59dqaZB8kjcccVvOeVuU7aeS2OzndDoe8UUTOw1rffpaIvAicZZ9PSiRe4IFTbzbOQ8diDlv6va+GUDInL/5vwzXDDQ3GThbjrpUhhJnNdM52WL0hAQasVPVRVRVVXW2HQ16jqr9S1W5VPUNVD7ZfS/Z9IvEa0Pp6itbXUyz/SbVL4kFQb6tiA1wB/aGLUsYuhLVc+0rWaZUkGrAKQOvrtbWpVv+BDa7zBXTd+2rpB9OZ4mFwWl2LEAYDOHHkFtpnMuY7GNYiWnnB5hOJtxQDzi/33mOztPrrDk8Jba86R4P7l4TQ5x4M2dc4pOgYsZGJH8laDY0TibdG6Eg6N38+btVrbNjo0/GhXggacK6l2bv29apRB3z6RBchX7Q1TZVr3qjPa/Pgaebxl+P9cYcddNPXSz/kk0N//kUO/XnAmie3HDDkPXlDwTVaK6k6i7hRxSWBUc0bkN3HOWuclldhpDNcB47AYs0RVqTIvoGJY7+xsbzIOv8/8Td6yBwwN3h+VUYUYlX0P4nEWwaxtDXauu28BE3BosSWxLmgv7D5OXD+AK2/LAxirq6R4N3nrWD+fVvNCpMv2oBoZ8f4sew3a07XJFGzuTb4+Ynf9rweS8u4cGuSMCr8ZNJpoz6b2EXme7WzwyHcckm3NRXYGaddHzi/qSKaKqox/rxrwfhx88JBRrdXZ7uOkTdbaFo4UVtlbT/avnMHabvXwENq2DkwF0bzWvf1GueRHwpW0rXlQ+6JhhgMMABRzeuiWO3rZmSB989s435x2KFfMB8Iy2ZlXLhl415369PXeNLydPc47Y2dvp7Tdu8fnHR7Y93GcK5mzRuJ1wfJA4v3/eY+lxm3RY+E8OmJOmxk56zSzwDiGgxK9mXZf+JS4+Lo0LDT+n32g9+cZC/lRGLc4nsDrYarCVwfFVLhxkIk3jLItGQcNrgi/KmX1/629KL7yUj2ZUn2mf2IZHfttqzXfJG89PQ5LCwHjlqh3FVFYRKJtwTHL97G8Yu38Y4T/lztohQw1iYOS/ntmhcJsJ4TrSk6OGzZS9sMM9ICO+PUfzQuX2hoNGBVkxy/ONgXr3Nzhr2rJ+ZCNQ6HXHeT54jwC1/wv0po/k/dO8xnGWsrPee6/21LaH/C258zG8IgUz4axP8ZSL68l9GVnRPnz1t96LEVJSPBVJ1owKrGeO7iLxjnofHCvXTLJqGQUMZ6zEPZ6P5+pxm6MEo87rAwdiRMPr9zXLiTUTO1b5Vr3pLiFZEmEXlSRDbZcWYL3H5E5CIR2WwvLF4vIie77sdF5GkRuSfv2hoReTzvmePD+ZMqT6ol5rD2l013qbdEa0r7U7tof2qX/zW8xZb8ZbMOE8NIkV7UY0xnofabzaPA6ao6YMfeeVRE7lXVx/PSPAjcraoqIquBO4DD8u5/AivEZVveta8CX1TVe+1IAl8FTjX4WypCYp/zX/bmidC5yceDQjjOFD5of8rM/cshzjB29/NaI5wt/c9ofGV3bTedlapu11JSvKqqQG5+oME2daXJnz9oyb8vIkuA84GvAPmdPGVCzO3AG2WWvWI8tn05AMmOUTL7zGqd0WUT86zLbrsBSRSKQ8to0cZca9qH5sdoD1Auya+dDQQrDc7/j/YPIIYRLBMDhXPT56y9lvvXX2uUbxjU/JJAEYkDG4BVwLdU9QmPNBcD12PFmz0/79Y3gM8A7jgnVwL3i8g/YTXfTyzy3pcDlwMsXWo+b1kOOdEGJT0r2K+ytKTRwcKPJtXi7OU0dyuj7WZum2IaFgfQnNjD2CbUlUfja3vJzmsrkrj61PySQFXNqOoaYAlwvIgc5ZHmTlU9DCvS+5cAROQCYLeqbvDI9q+Bv1PVA4G/w4qo5/Xet6rqWlVd29XV5ae4oXD78Z7FKclwlzDcJST3mZeh67H4uI22hSCMTMZhYugWqdnshHDNMrIsxD2CK0FuVVG++XpO5AcisltEns27dq2I7HAFpZuUstqAqtorIg8B64Bni6R5WERWikgnVqS8C+2CNAFtIvIjVX0/VnjLT9iP/Qfw/XLKUmsMdwWrATUdczSdD/jlhJtg2iASa7bd5ZnVPgt54bXgGULwnRaK5ldna3c9cHu2+eQ24JvAv7qu36yqvveJ9TPa3CUiHfZxM3AmsMWVZpUdjxYRORZoBLpV9WpVXaKqy4H3AL+1hQtWH/cU+/h04EW/ha42sZQ4rOcI8xVHB/yy0SHcIMx7Zph5zwwztrjDuDxeDhK+HitWE4+lnFZGALqaJUDcZrAqOMAgMr6Fn5p3IXC73e+NAXeo6j0icoVdkFuAS4DLRCQFDAOX2gNdk/GXwP8RkQQwgt2vrTUO6nTtdtfZzSt/WGaWqUvrmvHXXMzGvZtm854JYYc/KHugSmKxomLNv647dxObW3LHyklJzy5shpxx2vU8+LurjfI1wnu0uVNE1ued36qqt/rM8WMichmwHmsnwUk7X35GmzdjbYbkvn5L3vGNwI0l8nkIeCjv/FHguFLvX01euOR/c87vrzTKo/3wwh/Ynr1mMYqbd6fyjo2ysjDtt2acz2vGsP/q2u831jMAHuKtBTwGrEpud1KE72CNFan9+nXgw5M9ELlHTgEtO7P2K4webpZXtkFc52b5qXthQCYTyOlCUyE1e9MeTYmG+hi4EtWgfd4CVHV8cl5EvgfcM0lyIHKPDMSKE70Hflp2ZseFG4TG3ozDZu3yJ5BMk7cfZqal0WHpYw4OXLYcmkoHEm62x9UCjEn4A2BVIMhosxf25mI5LqbIgLDjvYO/3czg/lO+UfSexiZsaEHRZJMy5/FG5jzeyEiH+UeRGBgbt5olJpYZIFkdtzNPqaKfswIZdZoPROSnwGPAoSLyuoh8BPiqiDwjIpuB07CmTyclajb7pLXBKYjVp7zIpkeC12RzHg8nckQ2GRt/NRGteoz+ikmkSCDW7FoJNTxixX8OSFhN1DAJUiZVfa/H5bIdC6Ka1weT1b757Dna+9+ZenCew2ZvMhduNhkbF245ZNYe5jiPzZqFNJqXR5qbnOY3z2IDW7ka2rbG7cYzK6FjOWmowypJVPNOAW+eltfv3dFe4BdaLi2uvq9f0cYHx8i0eIsoNstfeJ3JcPsxG5Ev4mwGYvXx1azphQkR/ki3WoJ98ySzfAYXOj+SML4cDbtdG1cvXQx7g/tvhiVa7emdyHNuR8n0DTv2kVpsNl8cKqpVFW/UbDYk3ZodF25Zz7m7fgJ9B4W/UXd80F8/uGhN3NjgMDHos+bQnl6HcMslsXfQYeuOvMa4TIGwnTTyrZJE4i2DOY1DDjv1DD8LeaH/sEmmVoQCjyu/NPamHFasiVwOuakgTaXJjoya55dOO21/f+mHoKj7ZMOOEFZ8hElWnVZBIvH6xG8858noX5F12L4jy6+xx2bHx82UbH+/w3xHdmz0biG488v09hqX0S9VqX0VJJ11WCWJ+rxTRLzfKa6hxf5+lTNNQnykMG0YYrXewP6CzWmHfp+1oB/CWGjg3jpldAxazAfWpgpBkTCWRAYkEm8IaHP4IQRTLRNt6f0tCZp7DJtkmWBfsmzfALG2iZiy+REijQauXC6Z2jeAeO31W8vknDSqRNRsDsgL+7t4YX8XSw/zt9XHZKSTTjONjuFeypedFUJf2I7FbEQuEEAIZDqaC+ysk74cSt7+USSTcVgliWreMnlhf0jRPOLWL/a+t2SY/ULwJnFi/0jBNY2biT/mjjnl05c51txMdthD4GF9qQeHAJDBIegwH/U2Jqp564fvr73NPJO4jgu3HIbnFgpy4EDzeM6x9jaHxecECV/nRJLJcdNUCFvCDA6NC7fWkGzWYZUkqnkrSHxUyAQMShcq9mqe2IGLoc9w0Mrl/xxrbTUTrNufOoTpqilD1XtJY4WIat4yefiMr3leb+oYcVhiufmu8vGU08Zm+2sOF3UWKCOcjWe+Lv9l96BTMYpG2xgecVi2t852DFSsgcB8qyCReAOSjKfH7ZAVbxrnFx9zmt+Qor1HdfhKl50dsIkdj0M8bhx7GXB6axkuCwTLe8xt69Z83rycvlHLDzvfKkgk3gAUq319k8w6rP8t5v3Cfas7HNZ7pGHf1RatKeN7GRXbSsUvrhH0xMs1EKNfsZrN+VZBIvFOIe6mc6Y1AzHzPm9+EACNWcL1Q7bNe4Q2OzDoMK+1veUyLtoyyXa73B8TcajV+V/VgljYlSQasAqJpXP3sa3He8VLptX8Q0325kVjNG1y5nsyJZNguGOgO4CdX9FKS3PxueNESB5lU0p1B6wi8YbMaJ8ryuHcFIwFb+Dki9YIu9mabZ9FbHfwiBteXlVqOiLsCr6e3dtN7IAa3mAsh3oE9KsgUbPZgJWz9zrstGU+48anvGvOWMpphKFbkUD9TR11CjI3ImyKjow6zHeI2Cbv0K/S219g5y4vGf4pHAI2m4tsdzJXRB4QkRft15ILlyPxBuSBU282zuOA38THrWmfuVJjKXVY3yH+NuiKzZtb9J6RaHODXjkLY8lca8u4yUBIweYDowVLHn1yG9aWQfl8FnhQVQ/G2jL3s6UyiZrNFaJ5u9e/Orhgh+YX1qatO8zEIckQApuHGc7VLfYQYm2FihJokMrez2u56/JFTOxPfTvWBgV/P1k+kXinEOkP5987sHAKGkitzqV2sdZZ6M49JR9zrzIaJyTRih1xUvf3I7M93qeGUFWvPm/Q7U4WqOqbdr5vikjJTn8k3ikgqGgbB7KMtU4IdSLwnJl4062uGqt1Po27zTzAsn3O503jWok7TKxfEgnnWmLbmeTcw6/m3ueuNypTKfrZd/8D6Z91ui7vVVV3k3hKiMRryJvDTmeIo1dtZ/PTBxnl6Y4WGYT25/Yb55HD7V1lutVJwb7AYTSHw/AAK5OQRbpLRBbate5CoOQuVNGAlQG3HPdvxnnsPDvtsJ4TfU7jFOkutz+3Pxzh5jbgLhLyJihhxIgmFnNakTzPPbyKOwiWz91Ye1Zjv95V6oGo5q0gMVeF1b/U329nxytpeld4f1TNeyf6mpmkub9wUC8hTaUdTWfNid90EMyrRq3zPY7s7U5Oxeofvw58AbgBuMPe+mQb8Bel8ikpXhFpAh4Gknb6/1TVL7jSXIS1LWEWSANX2lt45u7HsfYc3aGqF+Rd/zjwMfuZX6rqZ0qVp96QrCWo0U6leaeZuBoMHaEKcK2Cka556E6zPUNNnBYKdiscHqlKc3iqKbLdCcAZ5eTjp+YdBU5X1QERaQAeFZF7VfXxvDQPAnerqorIauAOIH9fjU8AzwHjE48ichrW8PhqVR31M7pWT+REa8pou+Um2NytpJtK59mzpoO5G3sLrjfsc86JakMcMV3C5nLT9C1c9yBTSE3z+7bcUHBN5AZK7/Nen/jZXFuB3NBig23qSpM/9NiSf19ElgDnA18BPpmX7q+BG1R11M4jjG2iq8LmrYudF+aMQbdZczEn2qCk2ydqLJmC7640NxV4YXmho6PFm85hiDbnW1xFH+Nq4avTJSJxEdmINQL2gKo+4ZHmYhHZAvwS547e3wA+Q+EQyyHAO0TkCRH5vYi8NUD5q85dJ3/TOI8D74w5LCixtGV+1/gWQ0dHnTYUQgga15I+4+WGVViCV2v4+qaoakZV1wBLgONF5CiPNHeq6mHAu7D6v4jIBcBuVd3gkW0CmAOcAHwaq7Ne0C4UkctFZL2IrN+zp7QTQa2g87xrpYGDUw7TjHnzOidaI3r3Q+9+pIgPcVm4/H39il+1SDN+dNRpYzW8/3AFKetnXlV7sdy2is5vqerDwEoR6QROAi4Uka3AvwOni8iP7KSvA/+lFk9i1czuCW9U9VZVXauqa7u6QorcWEE0puM2tMjcf1kTTsuP7zwZ6WLRFm3RmpJbv2u6py9Yfedx29NtnN90paR4RaRLRDrs42bgTGCLK82qXK0pIscCjUC3ql6tqktUdTnwHuC3qvp++7FfAKfbzxxiP7M3hL+p4rz6fud84uz2YWa3D9O6otco37E2cVg2hCnSTGty3NKd5u6HYQk2FwMqjEAAbjwadNMCP6PNC4Hb7emeGHCHqt4jIlcAqOotwCXAZSKSAoaBS7X0EN8PgB/Yy6LGgA/6eKammd1uvspleF54i9BjY84+od/g67H2NrL7C4PBqbu5GoYowgjaZseOWrfsSu577Rvm+dUJUk96Wbt2ra5fv750wiqx+r8Lg5/1by2MJZVtcYpq6d3iKdr0rEJxeAWmS+6f+AzbXsob+PcQl1cUjsRLOwrLuL8PaW1xPuux6L7A1REPkWP5Lns+79XH9ii3uBdDxOKeAd+Kidf9PReRDaq61jNxnRB5WFUJyYtltf1dSucjZjWuQ7RBcLkYxro6Ua/dD3xSMD0URkDyBtfUUoWjNdYakW9zBWjojzksudPf/GZiyLtV5A5ANzInYPNVLUuvXFwy6aTkonXYFsa6YPd0VdbnwNW6ZVcav3e9ENW8IfPMhV/yvH7wV26a9Lm970jR+Yi3qGOuyLBq+pMbsKckTcmCpq9ms+YDVkFdKhsawOfuDCIy7TytIvGGyOb/cR2Ct3j94t4kLD5qJtbutxSOKM99Nmwn6fLJuVLq4BAxw7nlmTRIlU8k3hpi7ztSzPtD8Pmg/hUtpROViysUjTQ2eg5I+SW0aIsB9kOabrVvJN46JJaCbF4Lu2VnlsEDQxCuK7i5JBvRnt6Sj2kmU7TpLLOaHccZd1D1cglj18FpQiTeCvHiNZ/0vL7sh4UrYfzSsjOEEVx76ii9bD6Jnb3m+YWIe4mg7t6LhLAF6XQhEm+dkNwf8g50hrsuZD38leOzZ5d8Lj5vjmftW6ng5dOp6RyJN2RU1dgdb8wj3HLSh/uxuzk98axT+IMHNtOyo/xYzH6WAPolP7xOsS1AIyYnEm+Vee1Dztjar9mvR31q8qmlycgmJn48EiMhbGwWkmjV3t0+1pQ0EuxMHV12Ezlp1DuudbKxMX9iHVxcJLxMVh0mHeZ9TB0cGhduULJDQ2SHhsjsMN8LeboQ1bxTQBhN58lo7JvoH6abzX5/4/snRKXNjchgeE1jAB3y6WIp4hlYzqtvbcp0WWUUibeOyBdtOWSSQnx0Qhit28IRRGyWc9eFyaaM/BLKliszhEi8U4Rp7fvs172nlt553o2B8wxDtDlxxZNJI2cNCB7D+azEpTyQ/pnRe08HIvFOU4Y73T8cQuu20s9pS9Kz6RxWjegQvA/xxmbPJtvfH8p7Tzci8U4DwlzA70VswXzUMFSOaS0d1bSFROKtM4YWhP+RZVtctWpLkvjuwkgavgnJCeL+fd8PJZ/pSiTeOmP9Dz/F2g99PfDz8b0eTVC3eD2QjnbP2rcgkkYIwg1btF4eVdNhxDkS7xQy1VNGfmnbao1SDy1sZraXeMugIIyr6cZhU+gWOV3cIIsROWnUIet/+Clf6dq2ZsaFa0Q6Del0QUyrcpHmZodN1YZhqjrthQtRzTttUNcnGQ/D1yKsMKzxePDNs0swE0RajEi8dYpX7bvmb4L7Q0va6WucndtKbHdvyed0bKz4fG0Y8Zzz32sGC9WLSLwzEG12RYocGkMbzb4KmiqspUMJxk4k2mJEfd4pppJfvI3fdnplDS50ime0I8bedxxg/D6STDoslODreZzT9qHxfmsk3OJENe8MYLTD4DfatWg/dsB8dJ/53kYTGca4L5rPDUQk3mlG815X37UhQK1o13bZrnZi3QbOGkWIxBoOkXinGY/dcRVvf/c/TZqme+1c5q3vKbwRchNVR0a5f/jfQs0zYoJIvBWgVpw13Gj+gFJGQx0BiUQ79UTinSE09TidNQZXtDPrNcP9jfKQxgbu2/+D0PKLKE0k3mnM4AETVWmTRys5KDKnnXu33hxehhGBiMRbISrZdM4XbVlkKWg6Z+e1cf+GLxqXKSJ8Sn7KItIkIk+KyCYR+ZOIFHySInKRiGwWkY0isl5ETnbdj4vI0yJyj8ezV4mIikin2Z8SkWPzP3tH4SiX+zd8MRJuDeOn5h0FTlfVARFpAB4VkXtV9fG8NA8Cd6uqishq4A7gsLz7nwCeAxwRiUXkQOAswEeMh/qnVgeuHGTh/qcjwdYDJcWrlotLbmSjwTZ1pckf+WjJvy8iS4Dzga8A7irhZuAzwF3lFjzCnKFlrTzyi09XuxgRAfHV5xWROLABWAV8S1Wf8EhzMXA9MB9LrDm+gSXQ2a70FwI7VHXTZLWRiFwOXA6wdOlSP8WN8KBvWYJN/+fvql2MiDDJ9yEtZUAH8DvgqEnSvBP4jX18AfBt+/hU4B77eBbwBNBun28FOku9/3HHHaf1DuP70U+tRUwOsF7L+O7XopU1LKmqvcBDwLpJ0jwMrLQHoE4CLhSRrcC/A6eLyI+AlcBBwCb73hLgKREx95qf4eQ+2Ijpj5/R5i4R6bCPm4EzgS2uNKvEbvuKyLFAI9Ctqler6hJVXQ68B/itqr5fVZ9R1fmquty+9zpwrKruDPFvq0mmSliRaGcefvq8C4Hb7X5vDLhDVe8RkSsAVPUW4BLgMhFJAcPApRp9k6aM6F8bASD19EVYu3atrl+/vtrFMKac6aJ6+nzqCRHZoKprq10OEyIPqxolEm1EKSLx1hCRYCPKoa6azSKyh4n9p910AnsrWJxyqeXy1XLZYGrKt0xVu0LOs6LUlXgnQ0TW13IfppbLV8tlg9ovX7WIAtBFRNQpkXgjIuqU6STeW6tdgBLUcvlquWxQ++WrCtOmzxsRMdOYTjVvRMSMIhJvRESdUnXxikiHiPyniGwRkedE5O0i8jX7fLOI3JlbGGGnv1pEXhKR50XknLzrx4nIM/a9f85bKJEUkZ/Z158QkeV5z3xQRF607YOm5RORs0Rkg12ODSJyei2VL++ZpSIyICJX1Vr5RGS1iDwmVsilZ0SkaSrLV9dUe00icDvwUfu4EWvN8NlAwr52I3CjfXwEsAlIYi0pfBmI2/eeBN4OCHAvcK59/W+AW+zj9wA/s4/nAq/Yr3Ps4zmG5TsGWGQfH4UVbIBaKV/eMz8H/gO4qpbKh+Xxtxk42j6fN9Wfbz1btYXbBryKPXBWJM3FwI/t46uBq/Pu3W9/oAuBLXnX3wt8Nz9N3pdjr/0FGE9j3/su8F6T8rmuC9CN9UNTM+UD3gV8DbgWW7y1Uj7gPOBHHmmmpHz1btVuNq8A9gA/FCu65PdFxL39+oexfmkBFgPb8+69bl9bbB+7rzueUdU0sB/rF71YXibly+cS4GlVHa2V8tn3/h5wR5irifIBhwAqIveLyFMi8pkpLl9dU23xJoBjge+o6jHAIPDZ3E0RuQZIAz/OXfLIQye5HvSZoOXLXT8Sqzn4VyXKUOnyfRG4WZ0BA2upfAngZOB99uvFInLGFJavrqm2eF8HXteJgHb/ifVhYw8wXAC8T+12j53+wLznlwBv2NeXeFx3PCMiCaAd6JkkL5Py5aJl3glcpqov5+VTC+V7G/BVsUIPXQn8g4h8rIbK9zrwe1Xdq6pDwK/s9FNVvvqm2u124BHgUPv4Wqz+2Drgz0CXK+2ROAesXmFiQOOPwAlMDGicZ1//W5wDGnfYx3Ox+mNzbHsVmGtYvg67fJd45FP18rmeuxbngFXVy2fn8xRWgMIE8Bvg/KksXz1b9QsAa4D1WKOMv7D/0S9h9Vc22nZLXvprsEaZn8cecbSvrwWete99kwnvsSaskdWXsEYsV+Q982H7+kvAh0zLB3wOq2m4Mc/m10r5Soi3JsoHvB/4k12Wr051+erZIvfIiIg6pdp93oiIiIBE4o2IqFMi8UZE1CmReCMi6pRIvBERdUok3oiIOiUSb0REnfL/A9ZL6Bb3qWBxAAAAAElFTkSuQmCC\n", - "text/plain": [ - "
" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPgAAAEICAYAAAByNDmmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfQElEQVR4nO2de7hcVXn/P99cEIwgaAJGQgwQEJFCgBi1oNUIikBBShUoCC1WfnjFSkpFfSrUKwqPSkVjRAst3vASpeEmotxaCSSYBJBQYkCIqLkgShRCcs7398deA3uGmTP7zLnM3nPez/OsJ3uvy97v5Mx33rXetfdask0QBL3JuG4bEATByBECD4IeJgQeBD1MCDwIepgQeBD0MCHwIOhhQuAtkDRDkiVNaFH+QUkXF7jOJZI+NvwWDp30+WZ2245g5AiBJyQ9IOmQovVtf8L2P46kTWMFSedLuk/SY5JWSjq52zb1Ck29UxCMFJLG2+5ryP4T8NfA/wEvA66RtMr2/466gT1GeHBA0n8B04H/lrRR0lm54hMlPShpvaQP5dqcI+my3PnBkv5X0qOSHpL0903us62kn0q6UBmXSLpI0pXJey2WtHuu/l6SrpP0iKR7Jb0lV3a4pF+kdr+WNC/lT5a0KNnxiKSbJbX9O0s6QtLPJf0x2X9OruxKSe9pqL9C0psK2HmJpC9JukrSn4DXNt7b9kdsr7Tdb3sxcDPwynY2BwWwHSl7XPcB4JDc+QzAwFeAbYD9gE3AS1L5OcBl6Xg68BhwAjAReD4wK5VdAnws5d0GfCx3j0uAR4A5ZL2prwPfSmWTgIeAf0hlBwDrgZem8t8Ar0rHOwAHpONPAvOTHROBVwFq8ZkNzEzHrwH+guxHf1/gd8CbUtlbgMW5dvsBG4CtCth5CfAH4KB07a3b/B22SZ/tsG5/J3ohVc6DS/qapLWS7ipY/y3J090t6Rsd3PJc24/bXg4sJ/tyN3Ii8GPb37S92fYG28ty5S8EbgS+Y/vDDW2/b/s221vIBD4r5R8JPGD7P2xvsX0H8D3gb1P5ZmBvSdvZ/n0qr+VPBV6UbLnZSTkDYfsG23c686IrgG8Cf5WKfwjsIWmPdP5W4Nu2nyxgJ8APbf9PuvYTbUyZT/b/fG07m4P2VE7gZB7hsCIV0xfybOAg2y8F3tfB/X6bO/4z8JwmdXYBfjnANY4g80zzB3H9FwEvT13tRyU9SvZD8oJUfixwOPArSTdKqnVpPwOsAn4kabWkDwxg11NIenkaPqyT9AfgdGAygO1NwOXASam7fwLwXwXthMzDF7HhM8A+wFuK/CgF7amcwG3fRNatfQpJu0u6RtLSNObcKxW9HbjI9u9T27UDXXoIZj0E7D5A+VeAa4CrJE0axDVvtL19Lj3H9jsAbN9u+2hgR+AHZALE9mO2z7S9G1ng6v2SXlfgft8ArgB2sf1csh8j5covJRPu64A/2/5ZETsTbf9vJZ0LvBF4ve0/FrA3KEDlBN6CBcB7bB8IzAO+mPL3BPaU9D+SbpU0kOf/HbBbh/f/OnBIGg5MkPR8SbMa6rwbuBdYJGmbAtdcRGb7WyVNTOllkl4iaStJJ0p6ru3NwB+BPgBJR0qaKUm5/MaodTO2BR6x/YSkOcDf5QuToPuBC3jaew9oZ4F7kmw+O93vUNsbirYL2lN5gUt6DvCXwHckLQO+TDYGhSzoswdZAOkE4GJJ27e41CeBD6du5rzB2GD7QbLu8plkvYtlNIzVU5fzNDKP90NJW7e55mPA64HjgYfJuvLnAc9KVd4KPCDpj2Td6ZNS/h7Aj4GNwM+AL9q+ocDHeCfwb5IeA/6V1CNo4D/JAnFPzR4UsLMInyALVN6nbBZjo6QPDqJ90AJVcagjaQawyPY+krYD7rU9tUm9+cCtti9J59cDH7B9+2ja2ysoewDlNNsHd9uWoBiV9+BpvHa/pDcDpPnlmvf8AWneVdJksi776m7YWXUkPZvMyy/oti1BcSoncEnfJOt6vljSGklvIwv+vE3ScuBu4OhU/Vpgg6RfAD8F/jnGeINH0huAdWRxik6mGoMuUckuehAExaicBw+CoDiVetlk8uTJnjFjRrfNGFPcufZ3ded/seNOXbJkeFm6dOl621O6bcdIUymBz5gxgyVLlnTbjDHBbhdeAMDODflL3nvm6BszAkj6VbdtGA0qJfBgZKiJOeg9QuBjjBDz2CIE3uOEoMc2IfAeZCRFvduFF7C6R8bhY4EQeA8R3jpoJAReUULMQRFC4BUgxBx0SluBp9cabyJ7/W8C8F3bH2moczTwUbL3hbcA77N9Syp7gGy9sj5gi+3ZDW3nka1CMsX2+qF+oF4hRB0MB0U8+CZgru2NkiYCt0i62vatuTrXA1fYtqR9yd4l3itX/tpm4pW0C3Ao8GDnH6G3qIKwI9BWHdoKPC1UsDGd1lbqdEOdjbnTSY3lA/BZ4CyyRf3GJFUQdFBdCo3BJY0HlgIzydY4W9ykzjFkq6LsSLbIYA2TLQBo4Mu2F6T6RwG/tr08W12o5b1PI1sJhenTpxcxt7T0kpjDi1eDQgJ3thPFrLTc0UJJ+9i+q6HOwlT2arLxeG0boINsPyxpR+A6SSuBJcCHyJb6aXfvBaRFBmbPnl2Zd1t7ScxBdRlUFN32o5JuIFu2uOm65LZvSqucTra93vbDKX+tpIVki/z/HtgVqHnvacAdkubY/m2z61aBEHVQNopsaTOltlBhWg30EGBlQ53aKp5IOoBsx4sNkiZJ2jblTyLz2HelBfZ3tD3D9gxgDdnOHJUU924XXjAmxT0WP3PVKOLBpwKXpnH4OOBy24sknQ5gez7ZIvwnS9oMPA4clyLqO5F122v3+obta0big4wm8cUOqkKRKPoKYP8m+fNzx+eRLZXbWGc1zbf6aaw3o12dbhFiDqpMPMnWhBB10CuEwHOEsINeY8wKPMQ8PMR8eLkZEwIPMQdjlZ4UeAg6CDJ6RuAh6iB4JpXf+GCsPmRSJuL/v7xUyoPfufZ38WUKgkFQeQ8elIP44S0nIfAg6GFC4EHQw4TAg2EjuunlIwQeBD1MpQTeK1vXBsFoUSmBB0EwOELgQdDDVE7g8eZSuYlAW7monMCDIChOJQUeXjwIilFJgQdBUIwQeBD0MCHwYNiJQFt5qKzAYxweBO2prMCDIGhPCDwIepgQeBD0MCHwYESIQFs5qLTAI9AWBANTaYEH5Sa8ePcJgQdBD9NW4JK2lnSbpOWS7pZ0bpM6R0taIWmZpCWSDs6VPSDpzlpZLv8zklamdgslbT9snyooDeHFu0sRD74JmGt7P2AWcJikVzTUuR7Yz/Ys4FTg4oby19qeZXt2Lu86YB/b+wL/B5zdgf1BEAxAW4E7Y2M6nZiSG+pstF3Lm9RY3uK6P7K9JZ3eCkwrbHWOCLQFQWsKjcEljZe0DFgLXGd7cZM6x0haCVxJ5sVrGPiRpKWSTmtxi1OBq1vc+7TU7V+ybt26IuYGQZAoJHDbfan7PQ2YI2mfJnUW2t4LeBPw0VzRQbYPAN4IvEvSq/PtJH0I2AJ8vcW9F9iebXv2lClTipgbBKVG0tckrZV0Vy7vzSnG1S9p9gBtD5N0r6RVkj7Q7l6DiqLbfhS4AThsgDo3AbtLmpzOH07/rgUWAnNyxp4CHAmcmOviD5roppebCLQ9g0t4pobuAv4GuKlVI0njgYvInOXewAmS9h7oRkWi6FNqEW5J2wCHACsb6syUpHR8ALAVsEHSJEnbpvxJwOvTB0HSYcC/AEfZ/nM7O4KgV0hO8JGGvHts39um6Rxgle3Vtp8EvgUcPVCDIruLTgUuTb8e44DLbS+SdHoybD5wLHCypM3A48Bxti1pJ2Bh0v4E4Bu2r0nX/QLwLOC6VH6r7dML2NOU1e89sxSe4v4z5tWd7/r587tkSTDSvOG1z/b6R/rr8u5Yselu4Ilc1gLbC4bpljsDD+XO1wAvH6hBW4HbXgHs3yR/fu74POC8JnVWA/u1uO7MdveuEo3CDnqfdY/0cfM1L6jLe84LH3yiYTp4OFGTvAGHtpXaH7wshJgDAGM2u799xeFjDbBL7nwa8PBADULgbQgxDw+7XXhBzwVDDWxmVAV+O7CHpF2BXwPHA383UIMQ+AjT6gcixubVpx/zRAceXNI3gdcAkyWtAT5CFnT7d2AKcKWkZbbfIOmFwMW2D7e9RdK7gWuB8cDXbN890L166mWTXvMQvUYZgqDDiQ2bG1Kxdj7B9lTbE21Ps/3V9BzJNNvPsr2T7Tekug/bPjzX9irbe9re3fbH290rPHgQdIgRm90s7lUeesqDjwTRlQ5aYeAJj69LZSMEHowqvdRN70c8yfi6VDZ6TuAjMQ4fCS8e0fnqY8QT/RPqUtkon0VBUBGyabLyee08IfAg6BBbbC7huDtPz3XRIabLyk6vjMP7EU94q7pUNsKDB0GHGPFkePDuMNxePAJtQSMGNntCXSob5bMoCCpCv8UT/RO7bcaA9KwHrwrhxatL9iTb+LpUNkLgQVfohUBbCDwoRHjxamLEE55Yl8pGTwu8CoG2oLpkQbZye/AIsgVBh9hicwkfT83T0x58JAgvHtToR2zqn1CXykYIPOgaVQ+0RZCtR4mHXgKoregSAu8q8Vx6uamyF3eHXfQWWxc9T9J1ku5L/+7Qom3T7bhb0fMCD4KRwogt/ePrUkEu4ZlbF30AuN72HmTbcQ+071iz7bibEgIPgg6x6ciDN9u6iGwLokvT8aVkm3gOmRB40HWq2k03YovH1SWypZCX5FKrLbMb2cn2bwDSvzu2vG377bifIgTeIRFoCwzNuujra9tdpzRc+5LVGHA77kbGhMAj0BaMBLZ4sn98XRoCv5M0FSD9u7b5PVtvx92MMSHwKhFevDpkHnxcXRoCVwCnpONTgB82VhhoO+5WjBmBhxcPhpsWY/C2pK2Lfga8WNIaSW8DPgUcKuk+4NB0jqQXSroqNd0JuEXScuA24MrcdtxNKd+zdcGYpIqbE9rwZN/gu+W2T2hR9LomdR8GDk/HLbfjbkXbnxxJW0u6TdJySXdLOrdJnaMlrahNvks6OFfWdGK+6MT+cFKVt8uim14NjOjrH1eXykYRizYBc23vB8wCDpP0ioY61wP72Z4FnApc3FDebGJ+MBP7QVA6bDrqoo8mbS1yxsZ0OjElN9TZaLuWN6mxvAUjMrEfBKOH2NI3vi6VjUI/OZLGS1pGFrq/zvbiJnWOkbQSuJLMi9doNTFfaGJf0mm1hwbWrVtX6EMF1aRqD7zY0NevulQ2Cgncdl/qfk8D5kjap0mdhbb3IvPEH80VDWpivsl1F9QeGpgyZcpgmlaaGIeXHwN9HleXysagLLL9KHADz3xQPl/nJmB3SZPTeauJ+UIT+8NNVQJtY5VqeXGxpW9cXSobRaLoUyRtn463AQ4BVjbUmSlJ6fgAYCtgQ5uJ+bYT+0FQZrIuermj6EXmwacCl0oaT/aDcLntRZJOB7A9HzgWOFnSZuBx4DjblrQTsDBpfwLwjdzE/KeAy9Mk/4PAm4fzgwXBSGNEXwm9dp62Are9Ati/Sf783PF5wHlN6rScmLe9gSYT+1Vk18+fH2PmYaQyD704292kzJT752eMEz8a5cf9qktlY0wKfCS8QwTbxh429PeNq0tlI55FD4Ih4P5uWzAwIfAg6JhydsvzlK9PMUpUIogzhqnEfLjB/ePqUtkon0UVJpZxGoP0N6SSMaYFHl48GBKOKHowDIQXLy/qU10qGxFkC4JOsaCEXjtPePCgtFQi0NbBGHyIWxcdJuleSasktV0kJQQ+zMQyTmMId9xFv4QOti5K74NcRPbq9d7ACZL2HuhGY17gEWgLhoQbUpEmnW9dNAdYZXu17SeBb6V2LRnzAg+CoaB+1SVGduuinYGHcudrUl5LQuAjQDyXPkYwqK8+MbJbFzUbAwzYbwiBB6Wm9IG24XvQpcgKR2uAXXLn04CHB7poCJzqvF0WgbZyoc6DbM0ossLR7cAeknaVtBVwfGrXkhB4UHrK7MXVX58Ktelw6yLbW4B3A9cC95CtrnT3QPeKB12CoFNcXNR1zTrcuiidXwVc1VivFeHBg0pQVi/eiQcfTcKDB0GndOjBR5Pw4ImyB9p2/fz5Mf1WRkr+umh48BITgi43Mozr67YVAxMCLxkh6opRQq+dJ7roOeK59HJTxkBbBNmCoFfxU4+nlpbw4A3E5oTBYAgPHgQ9ipylMhMevGTE8+bVosnbZKUiBB5UilIF2hxd9IDW4/Dw1tWnjF47T3jwJsR0WbkpjRevgAdvK3BJW0u6TdJySXdLOrdJnaMlrZC0LC1Tc3BD+XhJP5e0KJc3S9KtuTZzhucjBcHoIMov8CJd9E3AXNsbJU0EbpF0te1bc3WuB66wbUn7ApcDe+XKzyB7f3W7XN6ngXNtXy3p8HT+miF8liAYXQzqK3cYva3AbRvYmE4npuSGOhtzp5Py5ZKmAUcAHwfen2/G04J/Lm2WnulFYo68c3a78IJSDKXK6LXzFAqypfWYlwIzgYtsL25S5xjgk2SrQR6RK/occBawbUOT9wHXSjqfbKjwly3ufRpwGsD06dOLmBsEo0bZBV4oyGa7z/YsskXe5kjap0mdhbb3IlvP+aMAko4E1tpe2uSy7wD+yfYuwD8BX21x7wW1FSqnTJlSxNxhoQzeISg3tbfJ8qlsDCqKbvtR4AaeuStDvs5NwO6SJgMHAUdJeoBskfa5ki5LVU8Bvp+Ov0O2qHsQVAr1uy4VaiOdIemuFLR+X5Py10j6QwpAL5P0r53aVySKPkXS9ul4G+AQYGVDnZmSlI4PALYCNtg+2/Y02zPIVoD8ie2TUrOHgb9Kx3OB+zr9EEHQFZqviz4gqff7djKHth9wpKQ9mlS92faslP6tUxOLePCpwE8lrSBbtvU624sknS7p9FTnWOAuScvI9k46LgXnBuLtwAWSlgOfII2zg6AoXZ8PT1H0fCrAS4Bbbf85rZJ6I3DMSJlYJIq+Ati/Sf783PF5wHltrnMDWfe+dn4LcGBxU0ef1e89s/tfoqDUNAmyTZa0JHe+oGF3k7uAj0t6PvA42Yqp+fo1Xpmc38PAvHbLI7ciHlUNgg6Rm46719ue3aqN7XsknQdcRzb9vBzY0lDtDuBF6dmTw4EfAM268W2JR1WDYAh0EkW3/VXbB9h+Ndkuo/c1lP+x9mxJWgd9YgpaD96+ThqNJWK6rNx0dQhloM/1qQCSdkz/Tgf+BvhmQ/kLckHrOWQ63dCJidFFD4IhUHRqrIHvpTH4ZuBdtn9fC1in2NbfAu+QtIVsnH58gaB1U0LgBYhgW9CM7EGXwevO9qua5OWD1l8AvjA06zJC4EEwBCr/skkQBC1w4bnvrhFBtqDydG341NmDLqNKePAgGAqdBdlGjfDgBYnpsnLTFS9u0Jb+ulQ2woMHQYcIo/7yiTpPCDwIOqX2oEuJiS560DOMfjfdqK+vLpWN8OBB0CnhwXuLCLQFjai/vy6VjfDgQdApNmwpX7c8T3jwQRJePHgKA3399alkhMCDnmJ0A22G/r76VDJC4B0QXjwAMg++pa8+lYwYgwdBp9hQwqmxPCHwIOiY8gfZQuBB0CkGl9yDxxg86DlGLdBW66LnU8kIgXdIBNoCMN6ypS4VocDWRZJ0oaRVklak3YI6IgQeBJ1iBu3BC25d9EayddD3INvx50udmhgCD4IOsY37+upSAYpsXXQ08J/OuBXYXtLUTmyMIFvQk+x24QUjPox6jN9fe92WbzduSLD1MGxdtDPwUO58Tcr7zWBtDIEHQYfYbrmN9gBtimxdpGZNB29hdNGHRATayk1Z17Jvt3URmcfeJXc+jWwTwkETAg+CUabd1kXAFcDJKZr+CuAPtgfdPYcCApe0taTbJC1PYf1zm9Q5OoXzl0laIunghvLxkn4uaVFD/nsk3Zuu++lOPkAQVJDvSfoF8N/kti6qbV8EXAWsBlYBXwHe2emNiozBNwFz01amE4FbJF2dons1rgeusG1J+wKXA3vlys8A7gG2q2VIei1ZtHBf25tqv2pB0Cn3nzHvGXk6Yx4dbus1YhTYusjAu4bjXm09eArVb0ynE1NyQ52Nuc3RJuXLJU0DjgAubrj0O4BP2d6UrrG2o08QBEFLCo3BUxd7GbAWuM724iZ1jpG0ErgSODVX9DngLKDxbfg9gVdJWizpRkkv68D+rhOBtqDMFBK47T7bs8iieXPS0ziNdRba3gt4E/BRAElHAmttL21y2QnADsArgH8GLq/tiZxH0mlpXL9k3bp1xT5VEATAIKPoth8FbgBazv/ZvgnYXdJk4CDgKEkPAN8C5kq6LFVdA3w/DQFuI/PwjQ8NYHuB7dm2Z0+ZMmUw5gbBmKdIFH2KpO3T8TbAIcDKhjoza943PRi/FbDB9tm2p9meARwP/MT2SanZD4C5qc2eqc36YfhMo05008tNk47hmKGIB58K/FTSCuB2sjH4ooaw/rHAXWmcfhFwXC7o1oqvAbtJuovMu59SoE0QtGTXz5/fbRNKh6qkqdmzZ3vJksbHdstDWZ+cGks0myoDnjFVJmmp7dmjYVM3iSfZgqCHCYEHPUV00+uJt8mGmVZdxPjidRdJpXuibTQIgQ8jq997Jmoh8GDkaPWjGkQXPRhDjMXpshB4EPQw0UUfJWJsHnSD8ODBmGKsddND4MPMWIzUBuUluuhdJrrugyei5sUJDx4EPUwIfASIbnq5GUvj8BB4EPQw8TbZCDISniLG5hnDMA6Pt8mCIKg2EUUPKkFEzjsjBB6UlhD10AmBV4z7z5jX8+Pw4RZ2szjTWImkxxh8BKlSALMXsT3m/wYh8AoSXdeBCWE/TXTRg0oTQh6YEHhFaebFe31sHmIePCHwoPSEsDsnxuAjzGh+OXtlbL7r589/ahwd4h4a4cGDrhNbP40cIfBg1AlBjx4h8B6jbA/ChJi7Swh8FLA9Zp6cqhHCLgch8GDYCFGXjxB40BEh5moQAh8lqtxNDzFXl7bz4JK2lnSbpOWS7pZ0bpM6R0taIWmZpCWSDm4oHy/p55IWNWk7T5IlTR7aRwlqDNd8+Or3nhnirjhFPPgmYK7tjZImArdIutr2rbk61wNX2LakfYHLgb1y5WcA9wDb5S8saRfgUODBoXyIqlAVLx6i7h3aCtzZo0Qb0+nElNxQZ2PudFK+XNI04Ajg48D7Gy7/WeAs4IeDNTwYHkLMvU2hMbik8cBSYCZwke3FTeocA3wS2JFM0DU+RybibRvqHwX82vbygbyapNOA0wCmT59exNygBSHmMUj+md92Cdge+CmwzwB1Xg38OB0fCXwxHb8GWJSOnw0sBp6bzh8AJre7/4EHHuiqQ9a7GfEUDAywxIP47lc1DeplE9uPAjcAhw1Q5yZg9xQ0Owg4StIDwLeAuZIuA3YHdgWWp7JpwB2SXjAYe4JnUvvDBgEUi6JPkbR9Ot4GOARY2VBnplI/W9IBwFbABttn255mewZwPPAT2yfZvtP2jrZnpLI1wAG2fzuMn62UjJT4QthBM4qMwacCl6Zx+DjgctuLJJ0OYHs+cCxwsqTNwOPAcY5v24gR/7VBUWJnky4wmKmyKv19qoSkMbGzSTzJVlJC2MFwEAIvESHqYLipVBdd0jrgVy2KJwPrR9GcwVJm+8psG4yMfS+yPWWYr1k6KiXwgZC0pMxjqjLbV2bboPz2lZlYdDEIepgQeBD0ML0k8AXdNqANZbavzLZB+e0rLT0zBg+C4Jn0kgcPgqCBEHgQ9DBdF7ik7SV9V9JKSfdIeqWkz6TzFZIW1l52SfXPlrRK0r2S3pDLP1DSnanswtzLL8+S9O2Uv1jSjFybUyTdl9IpQ7VP0qGSliY7lkqaWyb7cm2mS9ooaV4urxT2SdpX0s+ULQ92p6StR9K+nqfb76sClwL/mI63Invn/PXAhJR3HnBeOt4bWA48i+x1018C41PZbcArAQFXA29M+e8E5qfj44Fvp+PnAavTvzuk4x2GaN/+wAvT8T5kC1pQFvtybb4HfAeYVyb7yJ6sXAHsl86fP9J/315P3Rb3dsD9pGBfizrHAF9Px2cDZ+fKrk1/9KnAylz+CcCX83VyX6D16UvyVJ1U9mXghKHY15AvYAPZj1Fp7APeBHwGOIck8LLYBxwOXNakzojYNxZSt7vouwHrgP9QturqxZImNdQ5lewXG2Bn4KFc2ZqUt3M6bsyva2N7C/AHMs/Q6lpDsS/PscDPbW8qi32p7F+AxpVxS2EfsCdgSddKukPSWSNsX8/TbYFPAA4AvmR7f+BPwAdqhZI+BGwBvl7LanIND5DfaZtO7avlv5Ss6/n/2tgw2vadC3zW9Ytklsm+CcDBwInp32MkvW4E7et5ui3wNcAaP72I43fJvhCkoMiRwIlOfaxUf5dc+2nAwyl/WpP8ujaSJgDPBR4Z4FpDsa+2iuxC4GTbv8xdpwz2vRz4tLJlst4HfFDSu0tk3xrgRtvrbf8ZuCrVHyn7ep9ujxGAm4EXp+NzyMaHhwG/AKY01H0p9UG21TwdhLkdeAVPB2EOT/nvoj4Ic3k6fh7Z+HCHlO4HnjdE+7ZP9h3b5Dpdt6+h3TnUB9m6bl+6zh1ki3JOAH4MHDGS9vV66r4BMAtYQhY9/UH6Y6wiGz8tS2l+rv6HyKLn95IiqSl/NnBXKvsCTz+ltzVZxHgVWSR2t1ybU1P+KuAfhmof8GGybuiyXNqxLPY1tDuHeoGXwj7gJODuZMunR9q+Xk/xqGoQ9DDdHoMHQTCChMCDoIcJgQdBDxMCD4IeJgQeBD1MCDwIepgQeBD0MP8frwEkefY0+AYAAAAASUVORK5CYII=\n", - "text/plain": [ - "
" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPgAAAEICAYAAAByNDmmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfP0lEQVR4nO2dfbhcVXn2f3c+EIwgaA4YCTGEgIgUAqRRC1qNqAgURKqAILzFSvETKymF2quFqlUULpWKxogWWvzCjygNICIKSF8JJJAEkPCKASGi5gNRohCSc+73j70G9gwz5+wz58yZvec8v+taV/ZeH3s/kzP3PGs9a++1ZJsgCHqTCd02IAiCzhECD4IeJgQeBD1MCDwIepgQeBD0MCHwIOhhQuAtkDRTkiVNalH+T5IuKXCdSyV9ZPQtHDnp883uth1B5wiBJyQ9IOnQovVt/7vtv+2kTeMFSZ+Q9JCkP0j6paQPddumXiEEHowpkiY2yf4SsLftHYC/AN4m6c1ja1lvEgIHJP03MAP4H0mbJJ2VKz5R0oOSNuQ9i6RzJV2eOz9E0v+V9GjyRv+nyX22l/RjSRcp41JJF0u6StJjkpZK2iNXf29J10l6RNK9kt6aKztc0s9Su19JWpDyp0pakux4RNJPJA35d5Z0hKQ7khd9SNK5ubKrJL2vof4qSW8qYOelkj4v6WpJfwRe03hv2/fa/mMuawCIocNoYDtS9rjuA8ChufOZgIEvAtsB+wObgZek8nOBy9PxDOAx4ARgMvB8YE4quxT4SMq7FfhI7h6XAo8A84BJwFeAr6eyKcBDwN+ksgOBDcBLU/mvgVem452AA9Pxx4CFyY7JwCsBtfjMBman41cDf0b2o78f8FvgTansrcDSXLv9gY3ANgXsvBT4PXBwuva2LWw5G9iUbFoDTO/2d6IXUuU8uKQvS1on6a6C9d+aPN3dkr7axi3Ps/247ZXASrIvdyMnAj+0/TXbW2xvtL0iV/5C4Ebgm7b/uaHtd2zfansrmcDnpPwjgQds/6ftrbZvB74N/HUq3wLsI2kH279L5bX8acCLki0/cVLQYNi+wfadtgdsrwK+BvxlKv4esKekPdP524Fv2H6ygJ0A37P9v+naT7S4/8eB7cl+IP6b7EchGCGVEziZRzisSMX0hTwHONj2S4EPtHG/3+SO/wQ8p0md3YBfDHKNI8h6AQuHcf0XAS9LXe1HJT1K9kPyglR+LHA48EtJN0p6Rcr/JHAf8ANJaySdPYhdTyHpZWn4sF7S74HTgakAtjcDVwAnpe7+CWQiLGInZB5+SJxxB/A4cF6RNsHgVE7gtm8i69Y+haQ9JH1f0vI05tw7Fb0TuNj271LbdYNdegRmPQTsMUj5F4HvA1dLmjKMa95oe8dceo7tdwHYvs320cDOwHfJBIjtx2yfaXsW8FfAByW9tsD9vgpcCexm+7lkP0bKlV9GJtzXAn+y/dMidiaG+387icH/P4OCVE7gLVgEvM/2QcAC4HMpfy9gL0n/K+kWSYN5/t8Cs9q8/1eAQ9NwYJKk50ua01DnvcC9wBJJ2xW45hIy298uaXJKfy7pJZK2kXSipOfa3gL8AegHkHSkpNmSlMvvL3C/7YFHbD8haR7wtnxhEvQAcCFPe+9B7SxwTyRNkPR3knZKgcd5wHuA64u0Dwan8gKX9ByyqZVvSloBfIFsDAqZJ9iTLIB0AnCJpB1bXOpjwD+nbuaC4dhg+0Gy7vKZZL2LFTSM1dM4+DQyj/c9SdsOcc3HgNcDxwMPk3Xlzweelaq8HXhA0h/IutMnpfw9gR+SBax+CnzO9g0FPsa7gX+T9BjwL6QeQQP/RRaIe2r2oICdRTiGbIjzWLr2f6QUjBAViL+UDkkzgSW295W0A3Cv7WlN6i0EbrF9aTq/Hjjb9m1jaW+vIOlk4DTbh3TblqAYlffgtv8A3C/pLQCpm1fznt8lzbtKmkrWZV/TDTurjqRnk3n5Rd22JShO5QQu6WtkXc8XS1or6R1kwZ93SFoJ3A0cnapfC2yU9DPgx8A/2N7YDburjKQ3AOvJ4hTtTDUGXaKSXfQgCIpROQ8eBEFxmr4KWVamTp3qmTNndtuMccWd635bd/5nO+/SJUtGl+XLl2+w3ddtOzpNpQQ+c+ZMli1b1m0zxgWzLroQgF0b8pe9/8yxN6YDSPplt20YCyol8KAz1MQc9B4h8HFGiHl8EQLvcULQ45sQeA/SSVHPuuhC1vTIOHw8EALvIcJbB42EwCtKiDkoQgi8AoSYg3YZUuDptcabyF7/mwR8y/a/NtQ5Gvgw2fvCW4EP2L45lT1A9hpgP7DV9tyGtgvIViHps71hpB+oVwhRB6NBEQ++GZhve5OkycDNkq6xfUuuzvXAlbYtaT+yd4n3zpW/ppl4Je0GvA54sP2P0FtUQdgRaKsOQwo8LVSwKZ3WVup0Q51NudMpjeWD8CngLLJF/cYlVRB0UF0KjcGVLVa/nGyt6ottL21S5xiyVVF2JltksIbJFgA08AXbi1L9o4Bf2V6ZrS7U8t6nka2EwowZM4qYW1p6SczhxatBIYHb7gfmpOWOFkva1/ZdDXUWp7JXkY3Ha9sAHWz7YUk7A9dJWg0sAz5EttTPUPdeRFpkYO7cuZV5t7WXxBxUl2FF0W0/KukGsmWLm65LbvumtMrpVNsbbD+c8tdJWky2yP/vgN2BmveeDtwuaZ7t3zS7bhUIUQdlo8iWNn21hQrTaqCHAqsb6tRW8UTSgWQ7XmyUNEXS9il/CpnHvistsL+z7Zm2ZwJryXbmqKS4Z1104bgU93j8zFWjiAefBlyWxuETgCtsL5F0OoDthWSL8J8saQvZovXHpYj6LmTd9tq9vmr7+534IGNJfLGDqlAkir4KOKBJ/sLc8flkS+U21llD861+GuvNHKpOtwgxB1UmnmRrQog66BVC4DlC2EGvMW4FHmIeHWI+vNyMC4GHmIPxSk8KPAQdBBk9I/AQdRA8k8pvfDBeHzIpE/H/X14q5cHvXPfb+DIFwTCovAcPykH88JaTEHgQ9DAh8CDoYULgwagR3fTyEQIPgh6mUgLvla1rg2CsqJTAgyAYHiHwIOhhKifweHOp3ESgrVxUTuBBEBSnkgIPLx4ExaikwIMgKEYIPAh6mBB4MOpEoK08VFbgMQ4PgqGprMCDIBiaEHgQ9DAh8CDoYULgQUeIQFs5qLTAI9AWBINTaYEH5Sa8ePcJgQdBDzOkwCVtK+lWSSsl3S3pvCZ1jpa0StIKScskHZIre0DSnbWyXP4nJa1O7RZL2nHUPlVQGsKLd5ciHnwzMN/2/sAc4DBJL2+ocz2wv+05wKnAJQ3lr7E9x/bcXN51wL629wP+H3BOG/YHQTAIQwrcGZvS6eSU3FBnk+1a3pTG8hbX/YHtren0FmB6YatzRKAtCFpTaAwuaaKkFcA64DrbS5vUOUbSauAqMi9ew8APJC2XdFqLW5wKXNPi3qelbv+y9evXFzE3CIJEIYHb7k/d7+nAPEn7Nqmz2PbewJuAD+eKDrZ9IPBG4D2SXpVvJ+lDwFbgKy3uvcj2XNtz+/r6ipgbBKVG0pclrZN0Vy7vLSnGNSBp7iBtD5N0r6T7JJ091L2GFUW3/ShwA3DYIHVuAvaQNDWdP5z+XQcsBubljD0FOBI4MdfFHzbRTS83EWh7BpfyTA3dBbwZuKlVI0kTgYvJnOU+wAmS9hnsRkWi6H21CLek7YBDgdUNdWZLUjo+ENgG2ChpiqTtU/4U4PXpgyDpMOAfgaNs/2koO4KgV0hO8JGGvHts3ztE03nAfbbX2H4S+Dpw9GANiuwuOg24LP16TACusL1E0unJsIXAscDJkrYAjwPH2bakXYDFSfuTgK/a/n667meBZwHXpfJbbJ9ewJ6mrHn/maXwFPefsaDufPfPXNAlS4JO84bXPNsbHhmoy7t91ea7gSdyWYtsLxqlW+4KPJQ7Xwu8bLAGQwrc9irggCb5C3PH5wPnN6mzBti/xXVnD3XvKtEo7KD3Wf9IPz/5/gvq8p7zwgefaJgOHk3UJG/QoW2l9gcvCyHmAMCYLR4YuuLosRbYLXc+HXh4sAYh8CEIMY8Osy66sOeCoQa2MKYCvw3YU9LuwK+A44G3DdYgBN5hWv1AxNi8+gxgnmjDg0v6GvBqYKqktcC/kgXd/gPoA66StML2GyS9ELjE9uG2t0p6L3AtMBH4su27B7tXT71s0mseotcoQxB0NLFhS0Mq1s4n2J5me7Lt6ba/lJ4jmW77WbZ3sf2GVPdh24fn2l5tey/be9j+6FD3Cg8eBG1ixBY3i3uVh57y4J0gutJBKww84Yl1qWyEwIMxpZe66QOIJ5lYl8pGzwm8E+PwTnjxiM5XHyOeGJhUl8pG+SwKgoqQTZOVz2vnCYEHQZvYYksJx915eq6LDjFdVnZ6ZRw+gHjC29SlshEePAjaxIgnw4N3h9H24hFoCxoxsMWT6lLZKJ9FQVARBiyeGJjcbTMGpWc9eFUIL15dsifZJtalshECD7pCLwTaQuBBIcKLVxMjnvDkulQ2elrgVQi0BdUlC7KV24NHkC0I2sQWW0r4eGqenvbgnSC8eFBjALF5YFJdKhsh8KBrVD3QFkG2HiUeegmgtqJLCLyrxHPp5abKXtxtdtFbbF30PEnXSfp5+nenFm2bbsfdip4XeBB0CiO2DkysSwW5lGduXXQ2cL3tPcm24x5s37Fm23E3JQQeBG1i05YHb7Z1EdkWRJel48vINvEcMSHwoOtUtZtuxFZPqEtkSyEvy6VWW2Y3sovtXwOkf3dueduht+N+ihB4m0SgLTA066JvqG13ndJo7UtWY9DtuBsZFwKPQFvQCWzx5MDEujQCfitpGkD6d13ze7bejrsZ40LgVSK8eHXIPPiEujQCrgROScenAN9rrDDYdtytGDcCDy8ejDYtxuBDkrYu+inwYklrJb0D+DjwOkk/B16XzpH0QklXp6a7ADdLWgncClyV2467KeV7ti4Yl1Rxc0Ibnuwffrfc9gktil7bpO7DwOHpuOV23K0Y8idH0raSbpW0UtLdks5rUudoSatqk++SDsmVNZ2YLzqxP5pU5e2y6KZXAyP6BybUpbJRxKLNwHzb+wNzgMMkvbyhzvXA/rbnAKcClzSUN5uYH87EfhCUDpu2uuhjyZAWOWNTOp2ckhvqbLJdy5vSWN6CjkzsB8HYIbb2T6xLZaPQT46kiZJWkIXur7O9tEmdYyStBq4i8+I1Wk3MF5rYl3Ra7aGB9evXF/pQQTWp2gMvNvQPqC6VjUICt92fut/TgXmS9m1SZ7Htvck88YdzRcOamG9y3UW1hwb6+vqG07TSxDi8/Bjo94S6VDaGZZHtR4EbeOaD8vk6NwF7SJqazltNzBea2B9tqhJoG69Uy4uLrf0T6lLZKBJF75O0YzreDjgUWN1QZ7YkpeMDgW2AjUNMzA85sR8EZSbropc7il5kHnwacJmkiWQ/CFfYXiLpdADbC4FjgZMlbQEeB46zbUm7AIuT9icBX81NzH8cuCJN8j8IvGU0P1gQdBoj+kvotfMMKXDbq4ADmuQvzB2fD5zfpE7LiXnbG2kysV9Fdv/MBTFmHkUq89CLs91Nyky5f37GOfGjUX48oLpUNsalwDvhHSLYNv6wYaB/Ql0qG/EsehCMAA9024LBCYEHQduUs1uep3x9ijGiEkGccUwl5sMNHphQl8pG+SyqMLGM0zhkoCGVjHEt8PDiwYhwRNGDUSC8eHlRv+pS2YggWxC0iwUl9Np5woMHpaUSgbY2xuAj3LroMEn3SrpP0pCLpITAR5lYxmkc4ba76JfSxtZF6X2Qi8levd4HOEHSPoPdaNwLPAJtwYhwQyrSpP2ti+YB99leY/tJ4OupXUvGvcCDYCRoQHWJzm5dtCvwUO58bcprSQi8A8Rz6eMEg/rrE53duqjZGGDQfkMIPCg1pQ+0jd6DLkVWOFoL7JY7nw48PNhFQ+BU5+2yCLSVC7UfZGtGkRWObgP2lLS7pG2A41O7loTAg9JTZi+ugfpUqE2bWxfZ3gq8F7gWuIdsdaW7B7tXPOgSBO3i4qKua9bm1kXp/Grg6sZ6rQgPHlSCsnrxdjz4WBIePAjapU0PPpaEB0+UPdC2+2cuiOm3MlLy10XDg5eYEHS5kWFCf7etGJwQeMkIUVeMEnrtPNFFzxHPpZebMgbaIsgWBL2Kn3o8tbSEB28gNicMhkN48CDoUeQslZnw4CUjnjevFk3eJisVIfCgUpQq0Obooge0HoeHt64+ZfTaecKDNyGmy8pNabx4BTz4kAKXtK2kWyWtlHS3pPOa1Dla0ipJK9IyNYc0lE+UdIekJbm8OZJuybWZNzofKQjGBlF+gRfpom8G5tveJGkycLOka2zfkqtzPXClbUvaD7gC2DtXfgbZ+6s75PI+AZxn+xpJh6fzV4/gswTB2GJQf7nD6EMK3LaBTel0ckpuqLMpdzolXy5pOnAE8FHgg/lmPC345zLE0jO9SMyRt8+siy4sxVCqjF47T6EgW1qPeTkwG7jY9tImdY4BPka2GuQRuaJPA2cB2zc0+QBwraQLyIYKf9Hi3qcBpwHMmDGjiLlBMGaUXeCFgmy2+23PIVvkbZ6kfZvUWWx7b7L1nD8MIOlIYJ3t5U0u+y7g723vBvw98KUW915UW6Gyr6+viLmjQhm8Q1Buam+T5VPZGFYU3fajwA08c1eGfJ2bgD0kTQUOBo6S9ADZIu3zJV2eqp4CfCcdf5NsUfcgqBQacF0q1EY6Q9JdKWj9gSblr5b0+xSAXiHpX9q1r0gUvU/Sjul4O+BQYHVDndmSlI4PBLYBNto+x/Z02zPJVoD8ke2TUrOHgb9Mx/OBn7f7IYKgKzRfF31QUu/3nWQObX/gSEl7Nqn6E9tzUvq3dk0s4sGnAT+WtIps2dbrbC+RdLqk01OdY4G7JK0g2zvpuBScG4x3AhdKWgn8O2mcHQRF6fp8eIqi51MBXgLcYvtPaZXUG4FjOmVikSj6KuCAJvkLc8fnA+cPcZ0byLr3tfObgYOKmzr2rHn/md3/EgWlpkmQbaqkZbnzRQ27m9wFfFTS84HHyVZMzdev8Yrk/B4GFgy1PHIr4lHVIGgTuem4e4Ptua3a2L5H0vnAdWTTzyuBrQ3VbgdelJ49ORz4LtCsGz8k8ahqEIyAdqLotr9k+0DbryLbZfTnDeV/qD1bktZBn5yC1sO3r51G44mYLis3XR1CGeh3fSqApJ3TvzOANwNfayh/QS5oPY9MpxvbMTG66EEwAopOjTXw7TQG3wK8x/bvagHrFNv6a+BdkraSjdOPLxC0bkoIvAARbAuakT3oMnzd2X5lk7x80PqzwGdHZl1GCDwIRkDlXzYJgqAFLjz33TUiyBZUnq4Nn9p70GVMCQ8eBCOhvSDbmBEevCAxXVZuuuLFDdo6UJfKRnjwIGgTYTRQPlHnCYEHQbvUHnQpMdFFD3qGse+mG/X316WyER48CNolPHhvEYG2oBENDNSlshEePAjaxYat5euW5wkPPkzCiwdPYaB/oD6VjBB40FOMbaDNMNBfn0pGCLwNwosHQObBt/bXp5IRY/AgaBcbSjg1licEHgRtU/4gWwg8CNrF4JJ78BiDBz3HmAXaal30fCoZIfA2iUBbAMZbt9alIhTYukiSLpJ0n6RVabegtgiBB0G7mGF78IJbF72RbB30Pcl2/Pl8uyaGwIOgTWzj/v66VIAiWxcdDfyXM24BdpQ0rR0bI8gW9CSzLrqw48Oox/jdtddt/UbjhgTbjsLWRbsCD+XO16a8Xw/XxhB4ELSJ7ZbbaA/SpsjWRWrWdPgWRhd9RESgrdyUdS37obYuIvPYu+XOp5NtQjhsQuBBMMYMtXURcCVwcoqmvxz4ve1hd8+hgMAlbSvpVkkrU1j/vCZ1jk7h/BWSlkk6pKF8oqQ7JC1pyH+fpHvTdT/RzgcIggrybUk/A/6H3NZFte2LgKuBNcB9wBeBd7d7oyJj8M3A/LSV6WTgZknXpOhejeuBK21b0n7AFcDeufIzgHuAHWoZkl5DFi3cz/bm2q9aELTL/WcseEaezlhAm9t6dYwCWxcZeM9o3GtID55C9ZvS6eSU3FBnU25ztCn5cknTgSOASxou/S7g47Y3p2usa+sTBEHQkkJj8NTFXgGsA66zvbRJnWMkrQauAk7NFX0aOAtofBt+L+CVkpZKulHSn7dhf9eJQFtQZgoJ3Ha/7Tlk0bx56WmcxjqLbe8NvAn4MICkI4F1tpc3uewkYCfg5cA/AFfU9kTOI+m0NK5ftn79+mKfKggCYJhRdNuPAjcALef/bN8E7CFpKnAwcJSkB4CvA/MlXZ6qrgW+k4YAt5J5+MaHBrC9yPZc23P7+vqGY24QjHuKRNH7JO2YjrcDDgVWN9SZXfO+6cH4bYCNts+xPd32TOB44Ee2T0rNvgvMT232Sm02jMJnGnOim15umnQMxw1FPPg04MeSVgG3kY3BlzSE9Y8F7krj9IuB43JBt1Z8GZgl6S4y735KgTZB0JLdP3NBt00oHaqSpubOnetlyxof2y0PZX1yajzRbKoMeMZUmaTltueOhU3dJJ5kC4IeJgQe9BTRTa8n3iYbZVp1EeOL110kle6JtrEgBD6KrHn/maiFwIPO0epHNYguejCOGI/TZSHwIOhhoos+RsTYPOgG4cGDccV466aHwEeZ8RipDcpLdNG7THTdh09EzYsTHjwIepgQeAeIbnq5GU/j8BB4EPQw8TZZB+mEp4ixecYojMPjbbIgCKpNRNGDShCR8/YIgQelJUQ9ckLgFeP+Mxb0/Dh8tIXdLM40XiLpMQbvIFUKYPYitsf93yAEXkGi6zo4IeyniS56UGlCyIMTAq8ozbx4r4/NQ8zDJwQelJ4QdvvEGLzDjOWXs1fG5rt/5oKnxtEh7pERHjzoOrH1U+cIgQdjTgh67AiB9xhlexAmxNxdQuBjgO1x8+RUjRB2OQiBB6NGiLp8hMCDtggxV4MQ+BhR5W56iLm6DDkPLmlbSbdKWinpbknnNalztKRVklZIWibpkIbyiZLukLSkSdsFkixp6sg+SlBjtObD17z/zBB3xSniwTcD821vkjQZuFnSNbZvydW5HrjStiXtB1wB7J0rPwO4B9ghf2FJuwGvAx4cyYeoClXx4iHq3mFIgTt7lGhTOp2ckhvqbMqdTsmXS5oOHAF8FPhgw+U/BZwFfG+4hgejQ4i5tyk0Bpc0EVgOzAYutr20SZ1jgI8BO5MJusanyUS8fUP9o4Bf2V45mFeTdBpwGsCMGTOKmBu0IMQ8Dsk/8ztUAnYEfgzsO0idVwE/TMdHAp9Lx68GlqTjZwNLgeem8weAqUPd/6CDDnLVIevddDwFgwMs8zC++1VNw3rZxPajwA3AYYPUuQnYIwXNDgaOkvQA8HVgvqTLgT2A3YGVqWw6cLukFwzHnuCZ1P6wQQDFouh9knZMx9sBhwKrG+rMVupnSzoQ2AbYaPsc29NtzwSOB35k+yTbd9re2fbMVLYWOND2b0bxs5WSTokvhB00o8gYfBpwWRqHTwCusL1E0ukAthcCxwInS9oCPA4c5/i2dYz4rw2KEjubdIHhTJVV6e9TJSSNi51N4km2khLCDkaDEHiJCFEHo02luuiS1gO/bFE8FdgwhuYMlzLbV2bboDP2vch23yhfs3RUSuCDIWlZmcdUZbavzLZB+e0rM7HoYhD0MCHwIOhhekngi7ptwBCU2b4y2wblt6+09MwYPAiCZ9JLHjwIggZC4EHQw3Rd4JJ2lPQtSasl3SPpFZI+mc5XSVpce9kl1T9H0n2S7pX0hlz+QZLuTGUX5V5+eZakb6T8pZJm5tqcIunnKZ0yUvskvU7S8mTHcknzy2Rfrs0MSZskLcjllcI+SftJ+qmy5cHulLRtJ+3rebr9vipwGfC36XgbsnfOXw9MSnnnA+en432AlcCzyF43/QUwMZXdCrwCEHAN8MaU/25gYTo+HvhGOn4esCb9u1M63mmE9h0AvDAd70u2oAVlsS/X5tvAN4EFZbKP7MnKVcD+6fz5nf779nrqtrh3AO4nBfta1DkG+Eo6Pgc4J1d2bfqjTwNW5/JPAL6Qr5P7Am1IX5Kn6qSyLwAnjMS+hnwBG8l+jEpjH/Am4JPAuSSBl8U+4HDg8iZ1OmLfeEjd7qLPAtYD/6ls1dVLJE1pqHMq2S82wK7AQ7mytSlv13TcmF/XxvZW4PdknqHVtUZiX55jgTtsby6LfansH4HGlXFLYR+wF2BJ10q6XdJZHbav5+m2wCcBBwKft30A8Efg7FqhpA8BW4Gv1LKaXMOD5Lfbpl37avkvJet6/t0QNoy1fecBn3L9Ipllsm8ScAhwYvr3GEmv7aB9PU+3Bb4WWOunF3H8FtkXghQUORI40amPlervlms/HXg45U9vkl/XRtIk4LnAI4NcayT21VaRXQycbPsXueuUwb6XAZ9QtkzWB4B/kvTeEtm3FrjR9gbbfwKuTvU7ZV/v0+0xAvAT4MXp+Fyy8eFhwM+Avoa6L6U+yLaGp4MwtwEv5+kgzOEp/z3UB2GuSMfPIxsf7pTS/cDzRmjfjsm+Y5tcp+v2NbQ7l/ogW9ftS9e5nWxRzknAD4EjOmlfr6fuGwBzgGVk0dPvpj/GfWTjpxUpLczV/xBZ9PxeUiQ15c8F7kpln+Xpp/S2JYsY30cWiZ2Va3Nqyr8P+JuR2gf8M1k3dEUu7VwW+xranUu9wEthH3AScHey5ROdtq/XUzyqGgQ9TLfH4EEQdJAQeBD0MCHwIOhhQuBB0MOEwIOghwmBB0EPEwIPgh7m/wPdIhI49ELX8AAAAABJRU5ErkJggg==\n", - "text/plain": [ - "
" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" } ], "source": [ @@ -531,16 +495,6 @@ "text": [ "wrote data/freyberg_dis.shp\n" ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n" - ] } ], "source": [ @@ -575,8 +529,6 @@ "text": [ "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/utils/reference.py:2074: DeprecationWarning: crs has been deprecated and will be removed in version 3.3.5. Use CRS in shapefile_utils instead.\n", " warnings.warn(\n", - "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfdis.py:659: PendingDeprecationWarning: ModflowDis.thickness will be deprecated and removed in version 3.3.5. Use grid.thick().\n", - " warnings.warn(\n", "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/flopy/modflow/mfwel.py:280: DeprecationWarning: phiramp_unit will be replaced with iunitramp for consistency\n", " warnings.warn(\n" ] @@ -584,7 +536,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 8, @@ -723,24 +675,11 @@ "Util2d:rech_70: resetting 'how' to external\n", "Util2d:rech_71: resetting 'how' to external\n", "Util2d:rech_72: resetting 'how' to external\n", - "Util2d:rech_73: resetting 'how' to external" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "\n", + "Util2d:rech_73: resetting 'how' to external\n", "Util2d:rech_74: resetting 'how' to external\n", "Util2d:rech_75: resetting 'how' to external\n", "Util2d:rech_76: resetting 'how' to external\n", - "Util2d:rech_77: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_77: resetting 'how' to external\n", "Util2d:rech_78: resetting 'how' to external\n", "Util2d:rech_79: resetting 'how' to external\n", "Util2d:rech_80: resetting 'how' to external\n", @@ -868,25 +807,13 @@ "Util2d:rech_202: resetting 'how' to external\n", "Util2d:rech_203: resetting 'how' to external\n", "Util2d:rech_204: resetting 'how' to external\n", - "Util2d:rech_205: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_205: resetting 'how' to external\n", "Util2d:rech_206: resetting 'how' to external\n", "Util2d:rech_207: resetting 'how' to external\n", "Util2d:rech_208: resetting 'how' to external\n", "Util2d:rech_209: resetting 'how' to external\n", "Util2d:rech_210: resetting 'how' to external\n", - "Util2d:rech_211: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_211: resetting 'how' to external\n", "Util2d:rech_212: resetting 'how' to external\n", "Util2d:rech_213: resetting 'how' to external\n", "Util2d:rech_214: resetting 'how' to external\n", @@ -1020,25 +947,13 @@ "Util2d:rech_342: resetting 'how' to external\n", "Util2d:rech_343: resetting 'how' to external\n", "Util2d:rech_344: resetting 'how' to external\n", - "Util2d:rech_345: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_345: resetting 'how' to external\n", "Util2d:rech_346: resetting 'how' to external\n", "Util2d:rech_347: resetting 'how' to external\n", "Util2d:rech_348: resetting 'how' to external\n", "Util2d:rech_349: resetting 'how' to external\n", "Util2d:rech_350: resetting 'how' to external\n", - "Util2d:rech_351: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_351: resetting 'how' to external\n", "Util2d:rech_352: resetting 'how' to external\n", "Util2d:rech_353: resetting 'how' to external\n", "Util2d:rech_354: resetting 'how' to external\n", @@ -1168,24 +1083,12 @@ "Util2d:rech_478: resetting 'how' to external\n", "Util2d:rech_479: resetting 'how' to external\n", "Util2d:rech_480: resetting 'how' to external\n", - "Util2d:rech_481: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_481: resetting 'how' to external\n", "Util2d:rech_482: resetting 'how' to external\n", "Util2d:rech_483: resetting 'how' to external\n", "Util2d:rech_484: resetting 'how' to external\n", "Util2d:rech_485: resetting 'how' to external\n", - "Util2d:rech_486: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_486: resetting 'how' to external\n", "Util2d:rech_487: resetting 'how' to external\n", "Util2d:rech_488: resetting 'how' to external\n", "Util2d:rech_489: resetting 'how' to external\n", @@ -1329,24 +1232,12 @@ "Util2d:rech_627: resetting 'how' to external\n", "Util2d:rech_628: resetting 'how' to external\n", "Util2d:rech_629: resetting 'how' to external\n", - "Util2d:rech_630: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_630: resetting 'how' to external\n", "Util2d:rech_631: resetting 'how' to external\n", "Util2d:rech_632: resetting 'how' to external\n", "Util2d:rech_633: resetting 'how' to external\n", "Util2d:rech_634: resetting 'how' to external\n", - "Util2d:rech_635: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_635: resetting 'how' to external\n", "Util2d:rech_636: resetting 'how' to external\n", "Util2d:rech_637: resetting 'how' to external\n", "Util2d:rech_638: resetting 'how' to external\n", @@ -1494,25 +1385,13 @@ "Util2d:rech_780: resetting 'how' to external\n", "Util2d:rech_781: resetting 'how' to external\n", "Util2d:rech_782: resetting 'how' to external\n", - "Util2d:rech_783: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_783: resetting 'how' to external\n", "Util2d:rech_784: resetting 'how' to external\n", "Util2d:rech_785: resetting 'how' to external\n", "Util2d:rech_786: resetting 'how' to external\n", "Util2d:rech_787: resetting 'how' to external\n", "Util2d:rech_788: resetting 'how' to external\n", - "Util2d:rech_789: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_789: resetting 'how' to external\n", "Util2d:rech_790: resetting 'how' to external\n", "Util2d:rech_791: resetting 'how' to external\n", "Util2d:rech_792: resetting 'how' to external\n", @@ -1653,23 +1532,11 @@ "Util2d:rech_927: resetting 'how' to external\n", "Util2d:rech_928: resetting 'how' to external\n", "Util2d:rech_929: resetting 'how' to external\n", - "Util2d:rech_930: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_930: resetting 'how' to external\n", "Util2d:rech_931: resetting 'how' to external\n", "Util2d:rech_932: resetting 'how' to external\n", "Util2d:rech_933: resetting 'how' to external\n", - "Util2d:rech_934: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_934: resetting 'how' to external\n", "Util2d:rech_935: resetting 'how' to external\n", "Util2d:rech_936: resetting 'how' to external\n", "Util2d:rech_937: resetting 'how' to external\n", @@ -1810,24 +1677,12 @@ "Util2d:rech_1072: resetting 'how' to external\n", "Util2d:rech_1073: resetting 'how' to external\n", "Util2d:rech_1074: resetting 'how' to external\n", - "Util2d:rech_1075: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_1075: resetting 'how' to external\n", "Util2d:rech_1076: resetting 'how' to external\n", "Util2d:rech_1077: resetting 'how' to external\n", "Util2d:rech_1078: resetting 'how' to external\n", "Util2d:rech_1079: resetting 'how' to external\n", - "Util2d:rech_1080: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_1080: resetting 'how' to external\n", "Util2d:rech_1081: resetting 'how' to external\n", "Util2d:rech_1082: resetting 'how' to external\n", "Util2d:rech_1083: resetting 'how' to external\n", @@ -1844,13 +1699,7 @@ "Util2d:rech_1094: resetting 'how' to external\n", "Util2d:rech_1095: resetting 'how' to external\n", "Util2d:rech_1096: resetting 'how' to external\n", - "Util2d:rech_1097: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:rech_1097: resetting 'how' to external\n", "Util2d:hk: resetting 'how' to external\n", "Util2d:vk: resetting 'how' to external\n", "Util2d:ss: resetting 'how' to external\n", @@ -1862,13 +1711,7 @@ "Util2d:hk: resetting 'how' to external\n", "Util2d:vk: resetting 'how' to external\n", "Util2d:ss: resetting 'how' to external\n", - "Util2d:sy: resetting 'how' to external\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + "Util2d:sy: resetting 'how' to external\n", "FloPy is using the following executable to run the model: /Users/jdhughes/.local/bin/mfnwt\n", "\n", " MODFLOW-NWT-SWR1 \n", @@ -1880,19 +1723,13 @@ " SWR1 Version 1.04.0 09/15/2016 \n", "\n", " Using NAME file: freyberg.nam \n", - " Run start date and time (yyyy/mm/dd hh:mm:ss): 2021/08/06 16:45:54\n", + " Run start date and time (yyyy/mm/dd hh:mm:ss): 2021/08/17 13:28:12\n", "\n", " Solving: Stress period: 1 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 2 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 3 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 4 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 5 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 5 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 6 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 7 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 8 Time step: 1 Groundwater-Flow Eqn.\n", @@ -1908,13 +1745,7 @@ " Solving: Stress period: 18 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 19 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 20 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 21 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 21 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 22 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 23 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 24 Time step: 1 Groundwater-Flow Eqn.\n", @@ -1934,13 +1765,7 @@ " Solving: Stress period: 38 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 39 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 40 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 41 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 41 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 42 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 43 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 44 Time step: 1 Groundwater-Flow Eqn.\n", @@ -1960,13 +1785,7 @@ " Solving: Stress period: 58 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 59 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 60 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 61 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 61 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 62 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 63 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 64 Time step: 1 Groundwater-Flow Eqn.\n", @@ -1986,13 +1805,7 @@ " Solving: Stress period: 78 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 79 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 80 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 81 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 81 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 82 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 83 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 84 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2011,13 +1824,7 @@ " Solving: Stress period: 97 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 98 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 99 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 100 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 100 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 101 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 102 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 103 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2035,13 +1842,7 @@ " Solving: Stress period: 115 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 116 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 117 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 118 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 118 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 119 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 120 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 121 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2057,13 +1858,7 @@ " Solving: Stress period: 131 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 132 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 133 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 134 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 134 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 135 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 136 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 137 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2080,13 +1875,7 @@ " Solving: Stress period: 148 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 149 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 150 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 151 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 151 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 152 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 153 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 154 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2101,13 +1890,7 @@ " Solving: Stress period: 163 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 164 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 165 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 166 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 166 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 167 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 168 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 169 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2125,13 +1908,7 @@ " Solving: Stress period: 181 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 182 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 183 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 184 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 184 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 185 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 186 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 187 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2152,13 +1929,7 @@ " Solving: Stress period: 202 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 203 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 204 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 205 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 205 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 206 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 207 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 208 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2178,13 +1949,7 @@ " Solving: Stress period: 222 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 223 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 224 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 225 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 225 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 226 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 227 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 228 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2204,13 +1969,7 @@ " Solving: Stress period: 242 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 243 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 244 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 245 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 245 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 246 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 247 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 248 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2230,13 +1989,7 @@ " Solving: Stress period: 262 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 263 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 264 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 265 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 265 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 266 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 267 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 268 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2256,13 +2009,7 @@ " Solving: Stress period: 282 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 283 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 284 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 285 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 285 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 286 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 287 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 288 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2282,13 +2029,7 @@ " Solving: Stress period: 302 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 303 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 304 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 305 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 305 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 306 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 307 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 308 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2308,13 +2049,7 @@ " Solving: Stress period: 322 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 323 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 324 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 325 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 325 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 326 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 327 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 328 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2333,13 +2068,7 @@ " Solving: Stress period: 341 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 342 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 343 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 344 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 344 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 345 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 346 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 347 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2349,13 +2078,7 @@ " Solving: Stress period: 351 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 352 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 353 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 354 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 354 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 355 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 356 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 357 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2374,13 +2097,7 @@ " Solving: Stress period: 370 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 371 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 372 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 373 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 373 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 374 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 375 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 376 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2398,13 +2115,7 @@ " Solving: Stress period: 388 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 389 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 390 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 391 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 391 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 392 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 393 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 394 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2423,13 +2134,7 @@ " Solving: Stress period: 407 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 408 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 409 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 410 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 410 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 411 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 412 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 413 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2448,13 +2153,7 @@ " Solving: Stress period: 426 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 427 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 428 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 429 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 429 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 430 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 431 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 432 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2468,13 +2167,7 @@ " Solving: Stress period: 440 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 441 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 442 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 443 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 443 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 444 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 445 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 446 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2494,13 +2187,7 @@ " Solving: Stress period: 460 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 461 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 462 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 463 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 463 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 464 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 465 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 466 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2519,13 +2206,7 @@ " Solving: Stress period: 479 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 480 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 481 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 482 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 482 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 483 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 484 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 485 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2540,13 +2221,7 @@ " Solving: Stress period: 494 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 495 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 496 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 497 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 497 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 498 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 499 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 500 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2566,13 +2241,7 @@ " Solving: Stress period: 514 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 515 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 516 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 517 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 517 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 518 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 519 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 520 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2592,13 +2261,7 @@ " Solving: Stress period: 534 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 535 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 536 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 537 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 537 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 538 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 539 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 540 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2611,13 +2274,7 @@ " Solving: Stress period: 547 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 548 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 549 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 550 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 550 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 551 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 552 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 553 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2630,13 +2287,7 @@ " Solving: Stress period: 560 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 561 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 562 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 563 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 563 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 564 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 565 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 566 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2655,13 +2306,7 @@ " Solving: Stress period: 579 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 580 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 581 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 582 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 582 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 583 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 584 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 585 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2680,13 +2325,7 @@ " Solving: Stress period: 598 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 599 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 600 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 601 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 601 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 602 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 603 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 604 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2705,13 +2344,7 @@ " Solving: Stress period: 617 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 618 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 619 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 620 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 620 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 621 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 622 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 623 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2731,13 +2364,7 @@ " Solving: Stress period: 637 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 638 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 639 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 640 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 640 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 641 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 642 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 643 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2756,13 +2383,7 @@ " Solving: Stress period: 656 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 657 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 658 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 659 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 659 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 660 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 661 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 662 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2781,13 +2402,7 @@ " Solving: Stress period: 675 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 676 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 677 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 678 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 678 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 679 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 680 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 681 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2806,13 +2421,7 @@ " Solving: Stress period: 694 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 695 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 696 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 697 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 697 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 698 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 699 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 700 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2830,13 +2439,7 @@ " Solving: Stress period: 712 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 713 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 714 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 715 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 715 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 716 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 717 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 718 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2851,13 +2454,7 @@ " Solving: Stress period: 727 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 728 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 729 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 730 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 730 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 731 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 732 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 733 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2866,13 +2463,7 @@ " Solving: Stress period: 736 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 737 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 738 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 739 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 739 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 740 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 741 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 742 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2889,13 +2480,7 @@ " Solving: Stress period: 753 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 754 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 755 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 756 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 756 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 757 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 758 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 759 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2912,13 +2497,7 @@ " Solving: Stress period: 770 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 771 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 772 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 773 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 773 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 774 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 775 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 776 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2936,13 +2515,7 @@ " Solving: Stress period: 788 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 789 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 790 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 791 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 791 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 792 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 793 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 794 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2961,13 +2534,7 @@ " Solving: Stress period: 807 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 808 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 809 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 810 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 810 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 811 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 812 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 813 Time step: 1 Groundwater-Flow Eqn.\n", @@ -2985,13 +2552,7 @@ " Solving: Stress period: 825 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 826 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 827 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 828 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 828 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 829 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 830 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 831 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3011,13 +2572,7 @@ " Solving: Stress period: 845 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 846 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 847 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 848 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 848 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 849 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 850 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 851 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3035,13 +2590,7 @@ " Solving: Stress period: 863 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 864 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 865 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 866 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 866 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 867 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 868 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 869 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3059,13 +2608,7 @@ " Solving: Stress period: 881 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 882 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 883 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 884 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 884 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 885 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 886 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 887 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3083,13 +2626,7 @@ " Solving: Stress period: 899 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 900 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 901 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 902 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 902 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 903 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 904 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 905 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3107,13 +2644,7 @@ " Solving: Stress period: 917 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 918 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 919 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 920 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 920 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 921 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 922 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 923 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3131,13 +2662,7 @@ " Solving: Stress period: 935 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 936 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 937 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 938 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 938 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 939 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 940 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 941 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3154,13 +2679,7 @@ " Solving: Stress period: 952 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 953 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 954 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 955 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 955 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 956 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 957 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 958 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3178,13 +2697,7 @@ " Solving: Stress period: 970 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 971 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 972 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 973 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 973 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 974 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 975 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 976 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3201,13 +2714,7 @@ " Solving: Stress period: 987 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 988 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 989 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 990 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 990 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 991 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 992 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 993 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3225,13 +2732,7 @@ " Solving: Stress period: 1005 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1006 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1007 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1008 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 1008 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1009 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1010 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1011 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3250,13 +2751,7 @@ " Solving: Stress period: 1024 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1025 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1026 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1027 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 1027 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1028 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1029 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1030 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3269,13 +2764,7 @@ " Solving: Stress period: 1037 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1038 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1039 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1040 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 1040 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1041 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1042 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1043 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3295,13 +2784,7 @@ " Solving: Stress period: 1057 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1058 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1059 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1060 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 1060 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1061 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1062 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1063 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3319,13 +2802,7 @@ " Solving: Stress period: 1075 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1076 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1077 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1078 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ + " Solving: Stress period: 1078 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1079 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1080 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1081 Time step: 1 Groundwater-Flow Eqn.\n", @@ -3344,15 +2821,9 @@ " Solving: Stress period: 1094 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1095 Time step: 1 Groundwater-Flow Eqn.\n", " Solving: Stress period: 1096 Time step: 1 Groundwater-Flow Eqn.\n", - " Solving: Stress period: 1097 Time step: 1 Groundwater-Flow Eqn.\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - " Run end date and time (yyyy/mm/dd hh:mm:ss): 2021/08/06 16:46:07\n", - " Elapsed run time: 12.900 Seconds\n", + " Solving: Stress period: 1097 Time step: 1 Groundwater-Flow Eqn.\n", + " Run end date and time (yyyy/mm/dd hh:mm:ss): 2021/08/17 13:28:26\n", + " Elapsed run time: 14.047 Seconds\n", "\n", " Normal termination of simulation\n" ] @@ -4890,9 +4361,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.10" + "version": "3.8.11" } }, "nbformat": 4, - "nbformat_minor": 1 + "nbformat_minor": 4 } diff --git a/flopy/modflow/mfdis.py b/flopy/modflow/mfdis.py index 9a1135dbdc..805d5e72f6 100644 --- a/flopy/modflow/mfdis.py +++ b/flopy/modflow/mfdis.py @@ -21,8 +21,6 @@ ITMUNI = {"u": 0, "s": 1, "m": 2, "h": 3, "d": 4, "y": 5} LENUNI = {"u": 0, "f": 1, "m": 2, "c": 3} -warnings.simplefilter("always", PendingDeprecationWarning) - class ModflowDis(Package): """ @@ -287,8 +285,6 @@ def __init__( ) self.start_datetime = start_datetime - # calculate layer thicknesses - self.__calculate_thickness() self._totim = None @property @@ -633,36 +629,6 @@ def getbotm(self, k=None): else: return self.botm.array[k, :, :] - def __calculate_thickness(self): - # thk = [] - # thk.append(self.top - self.botm[0]) - # for k in range(1, self.nlay + sum(self.laycbd)): - # thk.append(self.botm[k - 1] - self.botm[k]) - self.__thickness = Util3d( - self.parent, - (self.nlay + sum(self.laycbd), self.nrow, self.ncol), - np.float32, - self.parent.modelgrid.thick, - name="thickness", - ) - - @property - def thickness(self): - """ - Return cell thicknesses. - - Returns - ------- - thickness : array of floats (nlay, nrow, ncol) - - """ - warnings.warn( - "ModflowDis.thickness will be deprecated and removed " - "in version 3.3.5. Use grid.thick().", - PendingDeprecationWarning, - ) - return self.parent.modelgrid.thick - def write_file(self, check=True): """ Write the package file. diff --git a/flopy/modflow/mfdisu.py b/flopy/modflow/mfdisu.py index 961b44e5fa..1eeaac4889 100644 --- a/flopy/modflow/mfdisu.py +++ b/flopy/modflow/mfdisu.py @@ -15,8 +15,6 @@ ITMUNI = {"u": 0, "s": 1, "m": 2, "h": 3, "d": 4, "y": 5} LENUNI = {"u": 0, "f": 1, "m": 2, "c": 3} -warnings.simplefilter("always", PendingDeprecationWarning) - class ModflowDisU(Package): """ @@ -475,9 +473,6 @@ def __init__( self.start_datetime = start_datetime - # calculate layer thicknesses - self.__calculate_thickness() - # get neighboring nodes self._get_neighboring_nodes() @@ -485,40 +480,6 @@ def __init__( self.parent.add_package(self) return - def __calculate_thickness(self): - # set ncol and nrow for array readers - nrow = None - ncol = self.nodelay.array - nlay = self.nlay - thk = [] - for k in range(self.nlay): - thk.append(self.top[k] - self.bot[k]) - self.__thickness = Util3d( - self.parent, - (nlay, nrow, ncol), - np.float32, - thk, - name="thickness", - ) - return - - @property - def thickness(self): - """ - Return cell thicknesses. - - Returns - ------- - thickness : array of floats (nodes,) - - """ - warnings.warn( - "ModflowDisU.thickness will be deprecated and removed " - "in version 3.3.5. Use grid.thick().", - PendingDeprecationWarning, - ) - return self.__thickness.array - def checklayerthickness(self): """ Check layer thickness. From b79b58b095a8ca548b724dc5535843a9076876ae Mon Sep 17 00:00:00 2001 From: Joseph D Hughes Date: Tue, 17 Aug 2021 14:58:14 -0500 Subject: [PATCH 2/2] refactor(thickness): remove deprecated mfdis and mfdisu thickness property --- autotest/t050_test.py | 36 ++++++++++++++++++------------------ 1 file changed, 18 insertions(+), 18 deletions(-) diff --git a/autotest/t050_test.py b/autotest/t050_test.py index 9516b87307..d844c95b0b 100644 --- a/autotest/t050_test.py +++ b/autotest/t050_test.py @@ -175,7 +175,7 @@ def test_vtk_export_packages(): # totalbytes = os.path.getsize(filetocheck) # assert(totalbytes==1019857) nlines = count_lines_in_file(filetocheck) - assert nlines == 8496 + assert nlines == 8491, "nlines ({}) not equal to 8491".format(nlines) # upw with point scalar output output_dir = os.path.join(cpth, "UPW") @@ -184,7 +184,7 @@ def test_vtk_export_packages(): # totalbytes1 = os.path.getsize(filetocheck) # assert(totalbytes1==2559173) nlines1 = count_lines_in_file(filetocheck) - assert nlines1 == 21215 + assert nlines1 == 21215, "nlines ({}) not equal to 21215".format(nlines) # bas with smoothing on output_dir = os.path.join(cpth, "BAS") @@ -528,7 +528,7 @@ def test_vtk_vti(): # totalbytes = os.path.getsize(filetocheck) # assert(totalbytes==1075) nlines = count_lines_in_file(filetocheck) - assert nlines == 21 + assert nlines == 17, "nlines ({}) not equal to 17".format(nlines) # with point scalar dis.export(output_dir + "_points", fmt="vtk", point_scalars=True) @@ -536,7 +536,7 @@ def test_vtk_vti(): # totalbytes1 = os.path.getsize(filetocheck) # assert(totalbytes1==2474) nlines1 = count_lines_in_file(filetocheck) - assert nlines1 == 38 + assert nlines1 == 29, "nlines1 ({}) not equal to 29".format(nlines1) # with binary dis.export(output_dir + "_bin", fmt="vtk", binary=True) @@ -554,7 +554,7 @@ def test_vtk_vti(): # totalbytes3 = os.path.getsize(filetocheck) # assert(totalbytes3==1606) nlines3 = count_lines_in_file(filetocheck) - assert nlines3 == 41 + assert nlines3 == 37, "nlines3 ({}) not equal to 37".format(nlines3) # force .vtu filenametocheck = "DIS.vtu" @@ -563,7 +563,7 @@ def test_vtk_vti(): # totalbytes4 = os.path.getsize(filetocheck) # assert(totalbytes4==5723) nlines4 = count_lines_in_file(filetocheck) - assert nlines4 == 129 + assert nlines4 == 125, "nlines4 ({}) not equal to 125".format(nlines4) # vector filenametocheck = "vect.vti" @@ -954,17 +954,17 @@ def test_vtk_export_true2d_nonregyz(): if __name__ == "__main__": - test_vtk_export_array2d() - test_vtk_export_array3d() - test_vtk_transient_array_2d() + # test_vtk_export_array2d() + # test_vtk_export_array3d() + # test_vtk_transient_array_2d() test_vtk_export_packages() - test_vtk_mf6() - test_vtk_binary_head_export() - test_vtk_cbc() - test_vtk_vector() + # test_vtk_mf6() + # test_vtk_binary_head_export() + # test_vtk_cbc() + # test_vtk_vector() test_vtk_vti() - test_vtk_vtr() - test_vtk_export_true2d_regular() - test_vtk_export_true2d_nonregxy() - test_vtk_export_true2d_nonregxz() - test_vtk_export_true2d_nonregyz() + # test_vtk_vtr() + # test_vtk_export_true2d_regular() + # test_vtk_export_true2d_nonregxy() + # test_vtk_export_true2d_nonregxz() + # test_vtk_export_true2d_nonregyz()