{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:08:46.621921Z", "iopub.status.busy": "2024-09-05T16:08:46.621758Z", "iopub.status.idle": "2024-09-05T16:08:47.720487Z", "shell.execute_reply": "2024-09-05T16:08:47.719986Z" } }, "outputs": [], "source": [ "import flopy.modflow as fm\n", "import mfexport" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## export a model\n", "\n", "#### load the model using flopy" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:08:47.723094Z", "iopub.status.busy": "2024-09-05T16:08:47.722592Z", "iopub.status.idle": "2024-09-05T16:08:47.907659Z", "shell.execute_reply": "2024-09-05T16:08:47.907163Z" } }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/mbase.py:97: UserWarning: The program mf2005 does not exist or is not executable.\n", " warn(\n", "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/pyproj/crs/crs.py:141: FutureWarning: '+init=:' syntax is deprecated. ':' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6\n", " in_crs_string = _prepare_from_proj_string(in_crs_string)\n" ] }, { "data": { "text/plain": [ "MODFLOW 3 layer(s) 137 row(s) 171 column(s) 12 stress period(s)" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "m = fm.Modflow.load('lpr_inset.nam',\n", " model_ws='data/lpr')\n", "m" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### define the location of the model grid\n", "* spacing in meters" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:08:47.936918Z", "iopub.status.busy": "2024-09-05T16:08:47.936348Z", "iopub.status.idle": "2024-09-05T16:08:47.947729Z", "shell.execute_reply": "2024-09-05T16:08:47.947218Z" } }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/utils/crs.py:114: PendingDeprecationWarning: the prj argument will be deprecated and will be removed in version 3.6. Use prjfile instead.\n", " warnings.warn(\n", "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/utils/crs.py:121: PendingDeprecationWarning: the epsg argument will be deprecated and will be removed in version 3.6. Use crs instead.\n", " warnings.warn(\n", "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/utils/crs.py:140: PendingDeprecationWarning: the proj4 argument will be deprecated and will be removed in version 3.6. Use crs instead.\n", " warnings.warn(\n", "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/utils/crs.py:121: PendingDeprecationWarning: the epsg argument will be deprecated and will be removed in version 3.6. Use crs instead.\n", " warnings.warn(\n", "/home/runner/micromamba/envs/mfexport_ci/lib/python3.12/site-packages/flopy/utils/crs.py:140: PendingDeprecationWarning: the proj4 argument will be deprecated and will be removed in version 3.6. Use crs instead.\n", " warnings.warn(\n" ] } ], "source": [ "grid = mfexport.MFexportGrid(delr=m.dis.delr.array * .3048, # grid spacing in meters\n", " delc=m.dis.delc.array * .3048, \n", " xul=557571, yul=446166, # upper left corner in CRS\n", " epsg=3070 # epsg reference for projected CRS\n", " )" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### export all data arrays and stress period data (`MfLists`) to shapefiles, rasters and PDFs\n", "* array data are written to compressed GeoTIFFs\n", "* transient lists are written to shapefiles\n", " * only model cells in the list are included\n", " * stress period data are included as attributes; only periods where stresses change are included\n", "* all data are displayed in PDFs using `matplotlib.pyplot.imshow`" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:08:47.949807Z", "iopub.status.busy": "2024-09-05T16:08:47.949426Z", "iopub.status.idle": "2024-09-05T16:09:04.424263Z", "shell.execute_reply": "2024-09-05T16:09:04.423713Z" } }, "outputs": [], "source": [ "%%capture\n", "mfexport.export(m, grid, output_path='postproc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### export a package" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:04.426734Z", "iopub.status.busy": "2024-09-05T16:09:04.426546Z", "iopub.status.idle": "2024-09-05T16:09:07.334326Z", "shell.execute_reply": "2024-09-05T16:09:07.333814Z" } }, "outputs": [], "source": [ "%%capture\n", "mfexport.export(m, grid, 'wel', output_path='postproc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### summarize model input" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:07.336800Z", "iopub.status.busy": "2024-09-05T16:09:07.336444Z", "iopub.status.idle": "2024-09-05T16:09:07.739160Z", "shell.execute_reply": "2024-09-05T16:09:07.738652Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "summarizing lpr_inset input...\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
packagevariablelayerminmeanmaxperiod
0DISbotm0.01033.5000001049.3065271066.699951NaN
1DISbotm1.0991.2399901004.0924571063.699951NaN
2DISbotm2.0988.2399901000.7403421042.599976NaN
3DISmodel_topNaN1071.0799561105.8371961190.459961NaN
4BAS6ibound0.01.0000001.0000001.000000NaN
5BAS6ibound1.01.0000001.0000001.000000NaN
6BAS6ibound2.01.0000001.0000001.000000NaN
7BAS6strt0.01069.8000491087.2566701107.099976NaN
8BAS6strt1.01069.9000241087.2477061107.099976NaN
9BAS6strt2.01070.0000001087.2474501107.099976NaN
10UPWhani0.00.0000000.0000000.000000NaN
11UPWhani1.00.0000000.0000000.000000NaN
12UPWhani2.00.0000000.0000000.000000NaN
13UPWhk0.01.005800241.21944810000.000000NaN
14UPWhk1.019.238001249.383617354.209991NaN
15UPWhk2.019.238001243.610578354.209991NaN
16UPWss0.00.0010000.0041131.000000NaN
17UPWss1.00.0010000.0010000.001000NaN
18UPWss2.00.0010000.0010000.001000NaN
19UPWsy0.00.1200000.1227421.000000NaN
\n", "
" ], "text/plain": [ " package variable layer min mean max period\n", "0 DIS botm 0.0 1033.500000 1049.306527 1066.699951 NaN\n", "1 DIS botm 1.0 991.239990 1004.092457 1063.699951 NaN\n", "2 DIS botm 2.0 988.239990 1000.740342 1042.599976 NaN\n", "3 DIS model_top NaN 1071.079956 1105.837196 1190.459961 NaN\n", "4 BAS6 ibound 0.0 1.000000 1.000000 1.000000 NaN\n", "5 BAS6 ibound 1.0 1.000000 1.000000 1.000000 NaN\n", "6 BAS6 ibound 2.0 1.000000 1.000000 1.000000 NaN\n", "7 BAS6 strt 0.0 1069.800049 1087.256670 1107.099976 NaN\n", "8 BAS6 strt 1.0 1069.900024 1087.247706 1107.099976 NaN\n", "9 BAS6 strt 2.0 1070.000000 1087.247450 1107.099976 NaN\n", "10 UPW hani 0.0 0.000000 0.000000 0.000000 NaN\n", "11 UPW hani 1.0 0.000000 0.000000 0.000000 NaN\n", "12 UPW hani 2.0 0.000000 0.000000 0.000000 NaN\n", "13 UPW hk 0.0 1.005800 241.219448 10000.000000 NaN\n", "14 UPW hk 1.0 19.238001 249.383617 354.209991 NaN\n", "15 UPW hk 2.0 19.238001 243.610578 354.209991 NaN\n", "16 UPW ss 0.0 0.001000 0.004113 1.000000 NaN\n", "17 UPW ss 1.0 0.001000 0.001000 0.001000 NaN\n", "18 UPW ss 2.0 0.001000 0.001000 0.001000 NaN\n", "19 UPW sy 0.0 0.120000 0.122742 1.000000 NaN" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df = mfexport.summarize(m)\n", "df.head(20)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Export model results\n", "\n", "#### cell by cell flows\n", "* written to compressed GeoTIFFs" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:07.741290Z", "iopub.status.busy": "2024-09-05T16:09:07.741112Z", "iopub.status.idle": "2024-09-05T16:09:07.802897Z", "shell.execute_reply": "2024-09-05T16:09:07.802378Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Exporting cell budget info...\n", "file: data/lpr/lpr_inset.cbc\n", "binary grid file: None\n", "stress period 0, timestep 4\n", "wrote postproc/rasters/STORAGE_per0_stp4.tif\n", "no data found at (4, 0) for b' CONSTANT HEAD'\n", "wrote postproc/rasters/CONSTANT HEAD_per0_stp4.tif\n", "no data found at (4, 0) for b'FLOW RIGHT FACE '\n", "wrote postproc/rasters/FLOW RIGHT FACE_per0_stp4.tif\n", "no data found at (4, 0) for b'FLOW FRONT FACE '\n", "wrote postproc/rasters/FLOW FRONT FACE_per0_stp4.tif\n", "no data found at (4, 0) for b'FLOW LOWER FACE '\n", "wrote postproc/rasters/FLOW LOWER FACE_per0_stp4.tif\n", "no data found at (4, 0) for b' WELLS'\n", "wrote postproc/rasters/WELLS_per0_stp4.tif\n", "no data found at (4, 0) for b' RECHARGE'\n", "wrote postproc/rasters/RECHARGE_per0_stp4.tif\n", "no data found at (4, 0) for b' STREAM LEAKAGE'\n", "wrote postproc/rasters/STREAM LEAKAGE_per0_stp4.tif\n" ] }, { "data": { "text/plain": [ "['postproc/rasters/STORAGE_per0_stp4.tif',\n", " 'postproc/rasters/CONSTANT HEAD_per0_stp4.tif',\n", " 'postproc/rasters/FLOW RIGHT FACE_per0_stp4.tif',\n", " 'postproc/rasters/FLOW FRONT FACE_per0_stp4.tif',\n", " 'postproc/rasters/FLOW LOWER FACE_per0_stp4.tif',\n", " 'postproc/rasters/WELLS_per0_stp4.tif',\n", " 'postproc/rasters/RECHARGE_per0_stp4.tif',\n", " 'postproc/rasters/STREAM LEAKAGE_per0_stp4.tif']" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "mfexport.export_cell_budget('data/lpr/lpr_inset.cbc', \n", " grid,\n", " kstpkper=(4, 0), output_path='postproc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### heads\n", "* compressed GeoTIFFs\n", "* shapefiles of equipotential contours" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:07.805053Z", "iopub.status.busy": "2024-09-05T16:09:07.804709Z", "iopub.status.idle": "2024-09-05T16:09:07.845807Z", "shell.execute_reply": "2024-09-05T16:09:07.845261Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Exporting heads...\n", "file: data/lpr/lpr_inset.hds\n", "stress period 0, timestep 4\n", "wrote postproc/rasters/wt_per0_stp4.tif\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/runner/work/modflow-export/modflow-export/mfexport/array_export.py:165: MatplotlibDeprecationWarning: The collections attribute was deprecated in Matplotlib 3.8 and will be removed in 3.10.\n", " for i, c in enumerate(ctr.collections):\n" ] }, { "data": { "text/plain": [ "['postproc/rasters/wt_per0_stp4.tif', 'postproc/shps/wt_ctr_per0_stp4.shp']" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "headsfile = 'data/lpr/lpr_inset.hds'\n", "\n", "mfexport.export_heads(headsfile, grid, hdry=m.upw.hdry, hnflo=m.bas6.hnoflo,\n", " kstpkper=(4, 0),\n", " output_path='postproc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### drawdown\n", "* compressed GeoTIFFs\n", "* shapefiles of equipotential contours" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:07.848222Z", "iopub.status.busy": "2024-09-05T16:09:07.847610Z", "iopub.status.idle": "2024-09-05T16:09:07.882042Z", "shell.execute_reply": "2024-09-05T16:09:07.881482Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Exporting drawdown...\n", "file: data/lpr/lpr_inset.hds\n", "from stress period 0, timestep 4\n", "to stress period 8, timestep 4\n", "\n", "\n", "wrote postproc/rasters/wt-ddn_per8_stp4.tif\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/runner/work/modflow-export/modflow-export/mfexport/array_export.py:165: MatplotlibDeprecationWarning: The collections attribute was deprecated in Matplotlib 3.8 and will be removed in 3.10.\n", " for i, c in enumerate(ctr.collections):\n" ] }, { "data": { "text/plain": [ "['postproc/rasters/wt-ddn_per8_stp4.tif',\n", " 'postproc/shps/wt-ddn_ctr_per8_stp4.shp']" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "mfexport.export_drawdown(headsfile, grid, hdry=m.upw.hdry, hnflo=m.bas6.hnoflo,\n", " kstpkper0=(4, 0),\n", " kstpkper1=(4, 8),\n", " output_path='postproc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### SFR results\n", "* shapefile of SFR water balance results\n", "* PDFs of simulated stream flows and stream/aquifer interactions" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "execution": { "iopub.execute_input": "2024-09-05T16:09:07.884080Z", "iopub.status.busy": "2024-09-05T16:09:07.883745Z", "iopub.status.idle": "2024-09-05T16:09:08.226470Z", "shell.execute_reply": "2024-09-05T16:09:08.225956Z" } }, "outputs": [], "source": [ "%%capture\n", "mf2005_sfr_outputfile = 'data/lpr/lpr_inset.sfr.out'\n", "outfiles = mfexport.export_sfr_results(mf2005_sfr_outputfile=mf2005_sfr_outputfile,\n", " model=m,\n", " grid=grid,\n", " kstpkper=(4, 0),\n", " output_length_units='feet',\n", " output_time_units='seconds',\n", " output_path='postproc'\n", " )" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.5" } }, "nbformat": 4, "nbformat_minor": 2 }