| \n", + " | cellid | \n", + "q | \n", + "boundname | \n", + "laytop | \n", + "laybot | \n", + "layer | \n", + "cell | \n", + "laymidpt | \n", + "rootname | \n", + "
|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "(0, 4736) | \n", + "0.0000 | \n", + "p1_ln__0 | \n", + "104.773590 | \n", + "69.462875 | \n", + "0 | \n", + "4736 | \n", + "17.655357 | \n", + "p1_ln | \n", + "
| 1 | \n", + "(0, 4735) | \n", + "0.0000 | \n", + "p2_ln__0 | \n", + "104.723152 | \n", + "69.486488 | \n", + "0 | \n", + "4735 | \n", + "17.618332 | \n", + "p2_ln | \n", + "
| 2 | \n", + "(0, 4742) | \n", + "0.0000 | \n", + "p3_ln__0 | \n", + "104.679207 | \n", + "69.459198 | \n", + "0 | \n", + "4742 | \n", + "17.610004 | \n", + "p3_ln | \n", + "
| 3 | \n", + "(0, 4739) | \n", + "0.0000 | \n", + "p4_ln__0 | \n", + "104.633415 | \n", + "69.517403 | \n", + "0 | \n", + "4739 | \n", + "17.558006 | \n", + "p4_ln | \n", + "
| 4 | \n", + "(0, 5138) | \n", + "0.0000 | \n", + "p5_ln__0 | \n", + "104.638702 | \n", + "69.421188 | \n", + "0 | \n", + "5138 | \n", + "17.608757 | \n", + "p5_ln | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "
| 178 | \n", + "(2, 7590) | \n", + "-0.0095 | \n", + "19/19__2 | \n", + "51.490337 | \n", + "-8.509663 | \n", + "2 | \n", + "7590 | \n", + "30.000000 | \n", + "19/19 | \n", + "
| 179 | \n", + "(0, 7057) | \n", + "-0.0095 | \n", + "19/20__0 | \n", + "104.024193 | \n", + "67.291756 | \n", + "0 | \n", + "7057 | \n", + "18.366219 | \n", + "19/20 | \n", + "
| 180 | \n", + "(2, 7057) | \n", + "-0.0095 | \n", + "19/20__2 | \n", + "51.473099 | \n", + "-8.526901 | \n", + "2 | \n", + "7057 | \n", + "30.000000 | \n", + "19/20 | \n", + "
| 181 | \n", + "(0, 5637) | \n", + "-0.0095 | \n", + "19/21__0 | \n", + "103.002068 | \n", + "68.493927 | \n", + "0 | \n", + "5637 | \n", + "17.254070 | \n", + "19/21 | \n", + "
| 182 | \n", + "(2, 5637) | \n", + "-0.0095 | \n", + "19/21__2 | \n", + "51.661217 | \n", + "-8.338783 | \n", + "2 | \n", + "5637 | \n", + "30.000000 | \n", + "19/21 | \n", + "
183 rows × 9 columns
\n", + "" + ], + "text/plain": [ + " cellid q boundname laytop laybot layer cell \\\n", + "0 (0, 4736) 0.0000 p1_ln__0 104.773590 69.462875 0 4736 \n", + "1 (0, 4735) 0.0000 p2_ln__0 104.723152 69.486488 0 4735 \n", + "2 (0, 4742) 0.0000 p3_ln__0 104.679207 69.459198 0 4742 \n", + "3 (0, 4739) 0.0000 p4_ln__0 104.633415 69.517403 0 4739 \n", + "4 (0, 5138) 0.0000 p5_ln__0 104.638702 69.421188 0 5138 \n", + ".. ... ... ... ... ... ... ... \n", + "178 (2, 7590) -0.0095 19/19__2 51.490337 -8.509663 2 7590 \n", + "179 (0, 7057) -0.0095 19/20__0 104.024193 67.291756 0 7057 \n", + "180 (2, 7057) -0.0095 19/20__2 51.473099 -8.526901 2 7057 \n", + "181 (0, 5637) -0.0095 19/21__0 103.002068 68.493927 0 5637 \n", + "182 (2, 5637) -0.0095 19/21__2 51.661217 -8.338783 2 5637 \n", + "\n", + " laymidpt rootname \n", + "0 17.655357 p1_ln \n", + "1 17.618332 p2_ln \n", + "2 17.610004 p3_ln \n", + "3 17.558006 p4_ln \n", + "4 17.608757 p5_ln \n", + ".. ... ... \n", + "178 30.000000 19/19 \n", + "179 18.366219 19/20 \n", + "180 30.000000 19/20 \n", + "181 17.254070 19/21 \n", + "182 30.000000 19/21 \n", + "\n", + "[183 rows x 9 columns]" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "spd" ] }, { "cell_type": "code", - "execution_count": null, - "id": "09d82a93", + "execution_count": 14, + "id": "77f90537", "metadata": {}, "outputs": [], "source": [ diff --git a/notebooks/modflow_setup/1_preprocessing.ipynb b/notebooks/modflow_setup/1_preprocessing.ipynb index 3d7339f..1dc5658 100644 --- a/notebooks/modflow_setup/1_preprocessing.ipynb +++ b/notebooks/modflow_setup/1_preprocessing.ipynb @@ -2,8 +2,8 @@ "cells": [ { "cell_type": "code", - "execution_count": null, - "id": "67d0c037", + "execution_count": 1, + "id": "ecbb7115", "metadata": {}, "outputs": [], "source": [ @@ -22,7 +22,7 @@ }, { "cell_type": "markdown", - "id": "02a58ec8", + "id": "070f0869", "metadata": {}, "source": [ "# Let's check out the proposed domain to get some coordiante information" @@ -30,8 +30,8 @@ }, { "cell_type": "code", - "execution_count": null, - "id": "2731b2fd", + "execution_count": 2, + "id": "78dac6b7", "metadata": {}, "outputs": [], "source": [ @@ -40,8 +40,8 @@ }, { "cell_type": "code", - "execution_count": null, - "id": "63ed9177", + "execution_count": 3, + "id": "2fa1b859", "metadata": {}, "outputs": [], "source": [ @@ -50,7 +50,7 @@ }, { "cell_type": "markdown", - "id": "ab335922", + "id": "471b54dc", "metadata": {}, "source": [ "### what's the projection?" @@ -58,17 +58,41 @@ }, { "cell_type": "code", - "execution_count": null, - "id": "5d7fa3f2", - "metadata": {}, - "outputs": [], + "execution_count": 4, + "id": "752d3c21", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "| \n", + " | Id | \n", + "Stage | \n", + "geometry | \n", + "
|---|---|---|---|
| 0 | \n", + "0 | \n", + "103.0 | \n", + "LINESTRING (1519819.286 5034332.958, 1519820.9... | \n", + "
| \n", + " | geometry | \n", + "segname | \n", + "from_id | \n", + "to_id | \n", + "streamwid | \n", + "
|---|---|---|---|---|---|
| 0 | \n", + "LINESTRING (1519819.286 5034482.958, 1519819.2... | \n", + "1001 | \n", + "0 | \n", + "1002 | \n", + "15 | \n", + "
| 1 | \n", + "LINESTRING (1519937.256 5033937.494, 1519964.7... | \n", + "1002 | \n", + "1001 | \n", + "1003 | \n", + "15 | \n", + "
| 2 | \n", + "LINESTRING (1520189.139 5033590.360, 1520233.5... | \n", + "1003 | \n", + "1002 | \n", + "1004 | \n", + "15 | \n", + "
| 3 | \n", + "LINESTRING (1520433.615 5033207.666, 1520420.9... | \n", + "1004 | \n", + "1003 | \n", + "1005 | \n", + "15 | \n", + "
| 4 | \n", + "LINESTRING (1520646.129 5033156.866, 1520692.6... | \n", + "1005 | \n", + "1004 | \n", + "1006 | \n", + "15 | \n", + "
| 5 | \n", + "LINESTRING (1520884.254 5032717.868, 1520880.0... | \n", + "1006 | \n", + "1005 | \n", + "1007 | \n", + "15 | \n", + "
| 6 | \n", + "LINESTRING (1520877.481 5032147.086, 1520871.5... | \n", + "1007 | \n", + "1006 | \n", + "1008 | \n", + "15 | \n", + "
| 7 | \n", + "LINESTRING (1520909.654 5031715.286, 1520910.5... | \n", + "1008 | \n", + "1007 | \n", + "1009 | \n", + "15 | \n", + "
| 8 | \n", + "LINESTRING (1521080.681 5031320.061, 1521118.7... | \n", + "1009 | \n", + "1008 | \n", + "1010 | \n", + "15 | \n", + "
| 9 | \n", + "LINESTRING (1521434.589 5030483.891, 1521474.3... | \n", + "1010 | \n", + "1009 | \n", + "0 | \n", + "15 | \n", + "
| \n", + " | cellid | \n", + "q | \n", + "boundname | \n", + "laytop | \n", + "laybot | \n", + "layer | \n", + "cell | \n", + "laymidpt | \n", + "rootname | \n", + "
|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "(0, 4736) | \n", + "0.0 | \n", + "p1_ln__0 | \n", + "104.773590 | \n", + "69.462875 | \n", + "0 | \n", + "4736 | \n", + "17.655357 | \n", + "p1_ln | \n", + "
| 1 | \n", + "(0, 4735) | \n", + "0.0 | \n", + "p2_ln__0 | \n", + "104.723152 | \n", + "69.486488 | \n", + "0 | \n", + "4735 | \n", + "17.618332 | \n", + "p2_ln | \n", + "
| 2 | \n", + "(0, 4742) | \n", + "0.0 | \n", + "p3_ln__0 | \n", + "104.679207 | \n", + "69.459198 | \n", + "0 | \n", + "4742 | \n", + "17.610004 | \n", + "p3_ln | \n", + "
| 3 | \n", + "(0, 4739) | \n", + "0.0 | \n", + "p4_ln__0 | \n", + "104.633415 | \n", + "69.517403 | \n", + "0 | \n", + "4739 | \n", + "17.558006 | \n", + "p4_ln | \n", + "
| 4 | \n", + "(0, 5138) | \n", + "0.0 | \n", + "p5_ln__0 | \n", + "104.638702 | \n", + "69.421188 | \n", + "0 | \n", + "5138 | \n", + "17.608757 | \n", + "p5_ln | \n", + "
| \n", + " | WellName | \n", + "X | \n", + "Y | \n", + "TopLayer | \n", + "BottomLaye | \n", + "DefaultQ | \n", + "geometry | \n", + "
|---|---|---|---|---|---|---|---|
| 0 | \n", + "p1_ln | \n", + "1.520070e+06 | \n", + "5.032527e+06 | \n", + "1 | \n", + "1 | \n", + "0.0 | \n", + "POINT (1520069.925 5032526.545) | \n", + "
| 1 | \n", + "p2_ln | \n", + "1.520053e+06 | \n", + "5.032516e+06 | \n", + "1 | \n", + "1 | \n", + "0.0 | \n", + "POINT (1520053.252 5032516.080) | \n", + "
| 2 | \n", + "p3_ln | \n", + "1.520072e+06 | \n", + "5.032513e+06 | \n", + "1 | \n", + "1 | \n", + "0.0 | \n", + "POINT (1520071.627 5032513.148) | \n", + "
| 3 | \n", + "p4_ln | \n", + "1.520033e+06 | \n", + "5.032496e+06 | \n", + "1 | \n", + "1 | \n", + "0.0 | \n", + "POINT (1520033.339 5032495.666) | \n", + "
| 4 | \n", + "p5_ln | \n", + "1.520070e+06 | \n", + "5.032489e+06 | \n", + "1 | \n", + "1 | \n", + "0.0 | \n", + "POINT (1520070.059 5032488.559) | \n", + "
| \n", + " | cellid | \n", + "q | \n", + "boundname | \n", + "laytop | \n", + "laybot | \n", + "layer | \n", + "cell | \n", + "laymidpt | \n", + "rootname | \n", + "
|---|---|---|---|---|---|---|---|---|---|
| 53 | \n", + "(0, 2196) | \n", + "-0.005 | \n", + "636__0 | \n", + "110.001999 | \n", + "65.353310 | \n", + "0 | \n", + "2196 | \n", + "22.324345 | \n", + "636 | \n", + "
| 54 | \n", + "(2, 2196) | \n", + "-0.005 | \n", + "636__2 | \n", + "49.060505 | \n", + "-10.939495 | \n", + "2 | \n", + "2196 | \n", + "30.000000 | \n", + "636 | \n", + "
| 55 | \n", + "(0, 1936) | \n", + "-0.005 | \n", + "637__0 | \n", + "110.503960 | \n", + "65.559921 | \n", + "0 | \n", + "1936 | \n", + "22.472019 | \n", + "637 | \n", + "
| 56 | \n", + "(2, 1936) | \n", + "-0.005 | \n", + "637__2 | \n", + "48.660793 | \n", + "-11.339207 | \n", + "2 | \n", + "1936 | \n", + "30.000000 | \n", + "637 | \n", + "
| 57 | \n", + "(0, 2217) | \n", + "-0.001 | \n", + "808__0 | \n", + "107.334236 | \n", + "69.746376 | \n", + "0 | \n", + "2217 | \n", + "18.793930 | \n", + "808 | \n", + "
| \n", + " | cellid | \n", + "q | \n", + "boundname | \n", + "laytop | \n", + "laybot | \n", + "layer | \n", + "cell | \n", + "laymidpt | \n", + "rootname | \n", + "X | \n", + "Y | \n", + "WellName | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "(0, 2196) | \n", + "-0.0050 | \n", + "636__0 | \n", + "110.001999 | \n", + "65.353310 | \n", + "0 | \n", + "2196 | \n", + "22.324345 | \n", + "636 | \n", + "1.518019e+06 | \n", + "5.032596e+06 | \n", + "636 | \n", + "
| 1 | \n", + "(2, 2196) | \n", + "-0.0050 | \n", + "636__2 | \n", + "49.060505 | \n", + "-10.939495 | \n", + "2 | \n", + "2196 | \n", + "30.000000 | \n", + "636 | \n", + "1.518019e+06 | \n", + "5.032596e+06 | \n", + "636 | \n", + "
| 2 | \n", + "(0, 1936) | \n", + "-0.0050 | \n", + "637__0 | \n", + "110.503960 | \n", + "65.559921 | \n", + "0 | \n", + "1936 | \n", + "22.472019 | \n", + "637 | \n", + "1.518020e+06 | \n", + "5.032742e+06 | \n", + "637 | \n", + "
| 3 | \n", + "(2, 1936) | \n", + "-0.0050 | \n", + "637__2 | \n", + "48.660793 | \n", + "-11.339207 | \n", + "2 | \n", + "1936 | \n", + "30.000000 | \n", + "637 | \n", + "1.518020e+06 | \n", + "5.032742e+06 | \n", + "637 | \n", + "
| 4 | \n", + "(0, 2217) | \n", + "-0.0010 | \n", + "808__0 | \n", + "107.334236 | \n", + "69.746376 | \n", + "0 | \n", + "2217 | \n", + "18.793930 | \n", + "808 | \n", + "1.519055e+06 | \n", + "5.032822e+06 | \n", + "808 | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "
| 73 | \n", + "(2, 7590) | \n", + "-0.0095 | \n", + "19/19__2 | \n", + "51.490337 | \n", + "-8.509663 | \n", + "2 | \n", + "7590 | \n", + "30.000000 | \n", + "19/19 | \n", + "1.520824e+06 | \n", + "5.032384e+06 | \n", + "19/19 | \n", + "
| 74 | \n", + "(0, 7057) | \n", + "-0.0095 | \n", + "19/20__0 | \n", + "104.024193 | \n", + "67.291756 | \n", + "0 | \n", + "7057 | \n", + "18.366219 | \n", + "19/20 | \n", + "1.520653e+06 | \n", + "5.032405e+06 | \n", + "19/20 | \n", + "
| 75 | \n", + "(2, 7057) | \n", + "-0.0095 | \n", + "19/20__2 | \n", + "51.473099 | \n", + "-8.526901 | \n", + "2 | \n", + "7057 | \n", + "30.000000 | \n", + "19/20 | \n", + "1.520653e+06 | \n", + "5.032405e+06 | \n", + "19/20 | \n", + "
| 76 | \n", + "(0, 5637) | \n", + "-0.0095 | \n", + "19/21__0 | \n", + "103.002068 | \n", + "68.493927 | \n", + "0 | \n", + "5637 | \n", + "17.254070 | \n", + "19/21 | \n", + "1.520382e+06 | \n", + "5.032503e+06 | \n", + "19/21 | \n", + "
| 77 | \n", + "(2, 5637) | \n", + "-0.0095 | \n", + "19/21__2 | \n", + "51.661217 | \n", + "-8.338783 | \n", + "2 | \n", + "5637 | \n", + "30.000000 | \n", + "19/21 | \n", + "1.520382e+06 | \n", + "5.032503e+06 | \n", + "19/21 | \n", + "
78 rows × 12 columns
\n", + "| \n", + " | cellid | \n", + "q | \n", + "boundname | \n", + "screen_top | \n", + "screen_botm | \n", + "layer | \n", + "cell | \n", + "laymidpt | \n", + "rootname | \n", + "x | \n", + "y | \n", + "WellName | \n", + "datetime | \n", + "enddatetime | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "(0, 2196) | \n", + "-0.0050 | \n", + "636__0 | \n", + "110.001999 | \n", + "65.353310 | \n", + "0 | \n", + "2196 | \n", + "22.324345 | \n", + "636 | \n", + "1.518019e+06 | \n", + "5.032596e+06 | \n", + "636 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 1 | \n", + "(2, 2196) | \n", + "-0.0050 | \n", + "636__2 | \n", + "49.060505 | \n", + "-10.939495 | \n", + "2 | \n", + "2196 | \n", + "30.000000 | \n", + "636 | \n", + "1.518019e+06 | \n", + "5.032596e+06 | \n", + "636 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 2 | \n", + "(0, 1936) | \n", + "-0.0050 | \n", + "637__0 | \n", + "110.503960 | \n", + "65.559921 | \n", + "0 | \n", + "1936 | \n", + "22.472019 | \n", + "637 | \n", + "1.518020e+06 | \n", + "5.032742e+06 | \n", + "637 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 3 | \n", + "(2, 1936) | \n", + "-0.0050 | \n", + "637__2 | \n", + "48.660793 | \n", + "-11.339207 | \n", + "2 | \n", + "1936 | \n", + "30.000000 | \n", + "637 | \n", + "1.518020e+06 | \n", + "5.032742e+06 | \n", + "637 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 4 | \n", + "(0, 2217) | \n", + "-0.0010 | \n", + "808__0 | \n", + "107.334236 | \n", + "69.746376 | \n", + "0 | \n", + "2217 | \n", + "18.793930 | \n", + "808 | \n", + "1.519055e+06 | \n", + "5.032822e+06 | \n", + "808 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "
| 73 | \n", + "(2, 7590) | \n", + "-0.0095 | \n", + "19/19__2 | \n", + "51.490337 | \n", + "-8.509663 | \n", + "2 | \n", + "7590 | \n", + "30.000000 | \n", + "19/19 | \n", + "1.520824e+06 | \n", + "5.032384e+06 | \n", + "19/19 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 74 | \n", + "(0, 7057) | \n", + "-0.0095 | \n", + "19/20__0 | \n", + "104.024193 | \n", + "67.291756 | \n", + "0 | \n", + "7057 | \n", + "18.366219 | \n", + "19/20 | \n", + "1.520653e+06 | \n", + "5.032405e+06 | \n", + "19/20 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 75 | \n", + "(2, 7057) | \n", + "-0.0095 | \n", + "19/20__2 | \n", + "51.473099 | \n", + "-8.526901 | \n", + "2 | \n", + "7057 | \n", + "30.000000 | \n", + "19/20 | \n", + "1.520653e+06 | \n", + "5.032405e+06 | \n", + "19/20 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 76 | \n", + "(0, 5637) | \n", + "-0.0095 | \n", + "19/21__0 | \n", + "103.002068 | \n", + "68.493927 | \n", + "0 | \n", + "5637 | \n", + "17.254070 | \n", + "19/21 | \n", + "1.520382e+06 | \n", + "5.032503e+06 | \n", + "19/21 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
| 77 | \n", + "(2, 5637) | \n", + "-0.0095 | \n", + "19/21__2 | \n", + "51.661217 | \n", + "-8.338783 | \n", + "2 | \n", + "5637 | \n", + "30.000000 | \n", + "19/21 | \n", + "1.520382e+06 | \n", + "5.032503e+06 | \n", + "19/21 | \n", + "2021-01-01 | \n", + "2021-12-31 | \n", + "
78 rows × 14 columns
\n", + "