This page was generated from feat_working_stack_examples.py. It's also available as a notebook.

FloPy working stack demo

A short demonstration of core flopy functionality

[1]:

import sys
from pathlib import Path
from pprint import pformat
from tempfile import TemporaryDirectory

[2]:

import matplotlib as mpl
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd

[3]:

from IPython.display import clear_output, display

proj_root = Path.cwd().parent.parent

# run installed version of flopy or add local path
import flopy

print(sys.version)
print(f"numpy version: {np.__version__}")
print(f"matplotlib version: {mpl.__version__}")
print(f"pandas version: {pd.__version__}")
print(f"flopy version: {flopy.__version__}")

3.12.4 | packaged by conda-forge | (main, Jun 17 2024, 10:23:07) [GCC 12.3.0]
numpy version: 2.0.1
matplotlib version: 3.9.1
pandas version: 2.2.2
flopy version: 3.8.0

Model Inputs

[4]:

# first lets load an existing model
model_ws = proj_root / "examples" / "data" / "freyberg_multilayer_transient"
ml = flopy.modflow.Modflow.load(
    "freyberg.nam",
    model_ws=model_ws,
    verbose=False,
    check=False,
    exe_name="mfnwt",
)

[5]:

ml.modelgrid

[5]:

xll:622241.1904510253; yll:3343617.741737109; rotation:15.0; crs:EPSG:32614; units:meters; lenuni:2

Let’s looks at some plots

[6]:

ml.upw.plot()

[6]:

[<Axes: title={'center': 'hk layer 1'}>,
 <Axes: title={'center': 'hk layer 2'}>,
 <Axes: title={'center': 'hk layer 3'}>,
 <Axes: title={'center': 'hani layer 1'}>,
 <Axes: title={'center': 'hani layer 2'}>,
 <Axes: title={'center': 'hani layer 3'}>,
 <Axes: title={'center': 'vka layer 1'}>,
 <Axes: title={'center': 'vka layer 2'}>,
 <Axes: title={'center': 'vka layer 3'}>,
 <Axes: title={'center': 'ss layer 1'}>,
 <Axes: title={'center': 'ss layer 2'}>,
 <Axes: title={'center': 'ss layer 3'}>,
 <Axes: title={'center': 'sy layer 1'}>,
 <Axes: title={'center': 'sy layer 2'}>,
 <Axes: title={'center': 'sy layer 3'}>,
 <Axes: title={'center': 'vkcb layer 1'}>,
 <Axes: title={'center': 'vkcb layer 2'}>,
 <Axes: title={'center': 'vkcb layer 3'}>]

../_images/Notebooks_feat_working_stack_examples_8_1.png

../_images/Notebooks_feat_working_stack_examples_8_2.png

../_images/Notebooks_feat_working_stack_examples_8_3.png

../_images/Notebooks_feat_working_stack_examples_8_4.png

../_images/Notebooks_feat_working_stack_examples_8_5.png

../_images/Notebooks_feat_working_stack_examples_8_6.png

../_images/Notebooks_feat_working_stack_examples_8_7.png

../_images/Notebooks_feat_working_stack_examples_8_8.png

../_images/Notebooks_feat_working_stack_examples_8_9.png

../_images/Notebooks_feat_working_stack_examples_8_10.png

../_images/Notebooks_feat_working_stack_examples_8_11.png

../_images/Notebooks_feat_working_stack_examples_8_12.png

../_images/Notebooks_feat_working_stack_examples_8_13.png

../_images/Notebooks_feat_working_stack_examples_8_14.png

../_images/Notebooks_feat_working_stack_examples_8_15.png

../_images/Notebooks_feat_working_stack_examples_8_16.png

../_images/Notebooks_feat_working_stack_examples_8_17.png

../_images/Notebooks_feat_working_stack_examples_8_18.png

[7]:

ml.dis.plot()

[7]:

[<Axes: title={'center': ' model_top'}>,
 <Axes: title={'center': 'botm layer 1'}>,
 <Axes: title={'center': 'botm layer 2'}>,
 <Axes: title={'center': 'botm layer 3'}>]

../_images/Notebooks_feat_working_stack_examples_9_1.png

../_images/Notebooks_feat_working_stack_examples_9_2.png

../_images/Notebooks_feat_working_stack_examples_9_3.png

../_images/Notebooks_feat_working_stack_examples_9_4.png

[8]:

ml.drn.plot(key="cond")
ml.drn.plot(key="elev")

[8]:

[<Axes: title={'center': ' DRN elev data stress period 1 layer 1'}>,
 <Axes: title={'center': ' DRN elev data stress period 1 layer 2'}>,
 <Axes: title={'center': ' DRN elev data stress period 1 layer 3'}>]

../_images/Notebooks_feat_working_stack_examples_10_1.png

../_images/Notebooks_feat_working_stack_examples_10_2.png

../_images/Notebooks_feat_working_stack_examples_10_3.png

../_images/Notebooks_feat_working_stack_examples_10_4.png

../_images/Notebooks_feat_working_stack_examples_10_5.png

../_images/Notebooks_feat_working_stack_examples_10_6.png

First create a temporary workspace.

[9]:

# create a temporary workspace
temp_dir = TemporaryDirectory()
workspace = Path(temp_dir.name)

Write a shapefile of the DIS package.

[10]:

# write the shapefile
ml.dis.export(workspace / "freyberg_dis.shp")

/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for delc array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for delr array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for laycbd array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for nstp array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for perlen array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for steady array, DIS package
  warn(
/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/shapefile_utils.py:319: UserWarning: Failed to get data for tsmult array, DIS package
  warn(

Write a netCDF file with all model inputs.

[11]:

ml.export(workspace / "freyberg.nc")

/home/runner/micromamba/envs/flopy/lib/python3.12/site-packages/flopy/export/netcdf.py:755: DeprecationWarning: datetime.datetime.utcnow() is deprecated and scheduled for removal in a future version. Use timezone-aware objects to represent datetimes in UTC: datetime.datetime.now(datetime.UTC).
  "date_created", datetime.utcnow().strftime("%Y-%m-%dT%H:%M:%SZ")

initialize_geometry::
model crs: EPSG:32614
initialize_geometry::nc_crs = EPSG:4326
transforming coordinates using = proj=pipeline step inv proj=utm zone=14 ellps=WGS84 step proj=unitconvert xy_in=rad xy_out=deg

[11]:

<flopy.export.netcdf.NetCdf at 0x7fef83a8e600>

Change model directory and external path, modify inputs and write new input files.

[12]:

ml.external_path = workspace / "ref"
ml.model_ws = workspace
ml.write_input()

Util2d:delr: resetting 'how' to external
Util2d:delc: resetting 'how' to external
Util2d:model_top: resetting 'how' to external
Util2d:botm_layer_0: resetting 'how' to external
Util2d:botm_layer_1: resetting 'how' to external
Util2d:botm_layer_2: resetting 'how' to external
Util2d:ibound_layer_0: resetting 'how' to external
Util2d:ibound_layer_1: resetting 'how' to external
Util2d:ibound_layer_2: resetting 'how' to external
Util2d:strt_layer_0: resetting 'how' to external
Util2d:strt_layer_1: resetting 'how' to external
Util2d:strt_layer_2: resetting 'how' to external
Util2d:rech_1: resetting 'how' to external
Util2d:rech_2: resetting 'how' to external
Util2d:rech_3: resetting 'how' to external
Util2d:rech_4: resetting 'how' to external
Util2d:rech_5: resetting 'how' to external
Util2d:rech_6: resetting 'how' to external
Util2d:rech_7: resetting 'how' to external
Util2d:rech_8: resetting 'how' to external
Util2d:rech_9: resetting 'how' to external
Util2d:rech_10: resetting 'how' to external
Util2d:rech_11: resetting 'how' to external
Util2d:rech_12: resetting 'how' to external
Util2d:rech_13: resetting 'how' to external
Util2d:rech_14: resetting 'how' to external
Util2d:rech_15: resetting 'how' to external
Util2d:rech_16: resetting 'how' to external
Util2d:rech_17: resetting 'how' to external
Util2d:rech_18: resetting 'how' to external
Util2d:rech_19: resetting 'how' to external
Util2d:rech_20: resetting 'how' to external
Util2d:rech_21: resetting 'how' to external
Util2d:rech_22: resetting 'how' to external
Util2d:rech_23: resetting 'how' to external
Util2d:rech_24: resetting 'how' to external
Util2d:rech_25: resetting 'how' to external
Util2d:rech_26: resetting 'how' to external
Util2d:rech_27: resetting 'how' to external
Util2d:rech_28: resetting 'how' to external
Util2d:rech_29: resetting 'how' to external
Util2d:rech_30: resetting 'how' to external
Util2d:rech_31: resetting 'how' to external
Util2d:rech_32: resetting 'how' to external
Util2d:rech_33: resetting 'how' to external
Util2d:rech_34: resetting 'how' to external
Util2d:rech_35: resetting 'how' to external
Util2d:rech_36: resetting 'how' to external
Util2d:rech_37: resetting 'how' to external
Util2d:rech_38: resetting 'how' to external
Util2d:rech_39: resetting 'how' to external
Util2d:rech_40: resetting 'how' to external
Util2d:rech_41: resetting 'how' to external
Util2d:rech_42: resetting 'how' to external
Util2d:rech_43: resetting 'how' to external
Util2d:rech_44: resetting 'how' to external
Util2d:rech_45: resetting 'how' to external
Util2d:rech_46: resetting 'how' to external
Util2d:rech_47: resetting 'how' to external
Util2d:rech_48: resetting 'how' to external
Util2d:rech_49: resetting 'how' to external
Util2d:rech_50: resetting 'how' to external
Util2d:rech_51: resetting 'how' to external
Util2d:rech_52: resetting 'how' to external
Util2d:rech_53: resetting 'how' to external
Util2d:rech_54: resetting 'how' to external
Util2d:rech_55: resetting 'how' to external
Util2d:rech_56: resetting 'how' to external
Util2d:rech_57: resetting 'how' to external
Util2d:rech_58: resetting 'how' to external
Util2d:rech_59: resetting 'how' to external
Util2d:rech_60: resetting 'how' to external
Util2d:rech_61: resetting 'how' to external
Util2d:rech_62: resetting 'how' to external
Util2d:rech_63: resetting 'how' to external
Util2d:rech_64: resetting 'how' to external
Util2d:rech_65: resetting 'how' to external
Util2d:rech_66: resetting 'how' to external
Util2d:rech_67: resetting 'how' to external
Util2d:rech_68: resetting 'how' to external
Util2d:rech_69: resetting 'how' to external
Util2d:rech_70: resetting 'how' to external
Util2d:rech_71: resetting 'how' to external
Util2d:rech_72: resetting 'how' to external
Util2d:rech_73: resetting 'how' to external
Util2d:rech_74: resetting 'how' to external
Util2d:rech_75: resetting 'how' to external
Util2d:rech_76: resetting 'how' to external
Util2d:rech_77: resetting 'how' to external
Util2d:rech_78: resetting 'how' to external
Util2d:rech_79: resetting 'how' to external
Util2d:rech_80: resetting 'how' to external
Util2d:rech_81: resetting 'how' to external
Util2d:rech_82: resetting 'how' to external
Util2d:rech_83: resetting 'how' to external
Util2d:rech_84: resetting 'how' to external
Util2d:rech_85: resetting 'how' to external
Util2d:rech_86: resetting 'how' to external
Util2d:rech_87: resetting 'how' to external
Util2d:rech_88: resetting 'how' to external
Util2d:rech_89: resetting 'how' to external
Util2d:rech_90: resetting 'how' to external
Util2d:rech_91: resetting 'how' to external
Util2d:rech_92: resetting 'how' to external
Util2d:rech_93: resetting 'how' to external
Util2d:rech_94: resetting 'how' to external
Util2d:rech_95: resetting 'how' to external
Util2d:rech_96: resetting 'how' to external
Util2d:rech_97: resetting 'how' to external
Util2d:rech_98: resetting 'how' to external
Util2d:rech_99: resetting 'how' to external
Util2d:rech_100: resetting 'how' to external
Util2d:rech_101: resetting 'how' to external
Util2d:rech_102: resetting 'how' to external
Util2d:rech_103: resetting 'how' to external
Util2d:rech_104: resetting 'how' to external
Util2d:rech_105: resetting 'how' to external
Util2d:rech_106: resetting 'how' to external
Util2d:rech_107: resetting 'how' to external
Util2d:rech_108: resetting 'how' to external
Util2d:rech_109: resetting 'how' to external
Util2d:rech_110: resetting 'how' to external
Util2d:rech_111: resetting 'how' to external
Util2d:rech_112: resetting 'how' to external
Util2d:rech_113: resetting 'how' to external
Util2d:rech_114: resetting 'how' to external
Util2d:rech_115: resetting 'how' to external
Util2d:rech_116: resetting 'how' to external
Util2d:rech_117: resetting 'how' to external
Util2d:rech_118: resetting 'how' to external
Util2d:rech_119: resetting 'how' to external
Util2d:rech_120: resetting 'how' to external
Util2d:rech_121: resetting 'how' to external
Util2d:rech_122: resetting 'how' to external
Util2d:rech_123: resetting 'how' to external
Util2d:rech_124: resetting 'how' to external
Util2d:rech_125: resetting 'how' to external
Util2d:rech_126: resetting 'how' to external
Util2d:rech_127: resetting 'how' to external
Util2d:rech_128: resetting 'how' to external
Util2d:rech_129: resetting 'how' to external
Util2d:rech_130: resetting 'how' to external
Util2d:rech_131: resetting 'how' to external
Util2d:rech_132: resetting 'how' to external
Util2d:rech_133: resetting 'how' to external
Util2d:rech_134: resetting 'how' to external
Util2d:rech_135: resetting 'how' to external
Util2d:rech_136: resetting 'how' to external
Util2d:rech_137: resetting 'how' to external
Util2d:rech_138: resetting 'how' to external
Util2d:rech_139: resetting 'how' to external
Util2d:rech_140: resetting 'how' to external
Util2d:rech_141: resetting 'how' to external
Util2d:rech_142: resetting 'how' to external
Util2d:rech_143: resetting 'how' to external
Util2d:rech_144: resetting 'how' to external
Util2d:rech_145: resetting 'how' to external
Util2d:rech_146: resetting 'how' to external
Util2d:rech_147: resetting 'how' to external
Util2d:rech_148: resetting 'how' to external
Util2d:rech_149: resetting 'how' to external
Util2d:rech_150: resetting 'how' to external
Util2d:rech_151: resetting 'how' to external
Util2d:rech_152: resetting 'how' to external
Util2d:rech_153: resetting 'how' to external
Util2d:rech_154: resetting 'how' to external
Util2d:rech_155: resetting 'how' to external
Util2d:rech_156: resetting 'how' to external
Util2d:rech_157: resetting 'how' to external
Util2d:rech_158: resetting 'how' to external
Util2d:rech_159: resetting 'how' to external
Util2d:rech_160: resetting 'how' to external
Util2d:rech_161: resetting 'how' to external
Util2d:rech_162: resetting 'how' to external
Util2d:rech_163: resetting 'how' to external
Util2d:rech_164: resetting 'how' to external
Util2d:rech_165: resetting 'how' to external
Util2d:rech_166: resetting 'how' to external
Util2d:rech_167: resetting 'how' to external
Util2d:rech_168: resetting 'how' to external
Util2d:rech_169: resetting 'how' to external
Util2d:rech_170: resetting 'how' to external
Util2d:rech_171: resetting 'how' to external
Util2d:rech_172: resetting 'how' to external
Util2d:rech_173: resetting 'how' to external
Util2d:rech_174: resetting 'how' to external
Util2d:rech_175: resetting 'how' to external
Util2d:rech_176: resetting 'how' to external
Util2d:rech_177: resetting 'how' to external
Util2d:rech_178: resetting 'how' to external
Util2d:rech_179: resetting 'how' to external
Util2d:rech_180: resetting 'how' to external
Util2d:rech_181: resetting 'how' to external
Util2d:rech_182: resetting 'how' to external
Util2d:rech_183: resetting 'how' to external
Util2d:rech_184: resetting 'how' to external
Util2d:rech_185: resetting 'how' to external
Util2d:rech_186: resetting 'how' to external
Util2d:rech_187: resetting 'how' to external
Util2d:rech_188: resetting 'how' to external
Util2d:rech_189: resetting 'how' to external
Util2d:rech_190: resetting 'how' to external
Util2d:rech_191: resetting 'how' to external
Util2d:rech_192: resetting 'how' to external
Util2d:rech_193: resetting 'how' to external
Util2d:rech_194: resetting 'how' to external
Util2d:rech_195: resetting 'how' to external
Util2d:rech_196: resetting 'how' to external
Util2d:rech_197: resetting 'how' to external
Util2d:rech_198: resetting 'how' to external
Util2d:rech_199: resetting 'how' to external
Util2d:rech_200: resetting 'how' to external
Util2d:rech_201: resetting 'how' to external
Util2d:rech_202: resetting 'how' to external
Util2d:rech_203: resetting 'how' to external
Util2d:rech_204: resetting 'how' to external
Util2d:rech_205: resetting 'how' to external
Util2d:rech_206: resetting 'how' to external
Util2d:rech_207: resetting 'how' to external
Util2d:rech_208: resetting 'how' to external
Util2d:rech_209: resetting 'how' to external
Util2d:rech_210: resetting 'how' to external
Util2d:rech_211: resetting 'how' to external
Util2d:rech_212: resetting 'how' to external
Util2d:rech_213: resetting 'how' to external
Util2d:rech_214: resetting 'how' to external
Util2d:rech_215: resetting 'how' to external
Util2d:rech_216: resetting 'how' to external
Util2d:rech_217: resetting 'how' to external
Util2d:rech_218: resetting 'how' to external
Util2d:rech_219: resetting 'how' to external
Util2d:rech_220: resetting 'how' to external
Util2d:rech_221: resetting 'how' to external
Util2d:rech_222: resetting 'how' to external
Util2d:rech_223: resetting 'how' to external
Util2d:rech_224: resetting 'how' to external
Util2d:rech_225: resetting 'how' to external
Util2d:rech_226: resetting 'how' to external
Util2d:rech_227: resetting 'how' to external
Util2d:rech_228: resetting 'how' to external
Util2d:rech_229: resetting 'how' to external
Util2d:rech_230: resetting 'how' to external
Util2d:rech_231: resetting 'how' to external
Util2d:rech_232: resetting 'how' to external
Util2d:rech_233: resetting 'how' to external
Util2d:rech_234: resetting 'how' to external
Util2d:rech_235: resetting 'how' to external
Util2d:rech_236: resetting 'how' to external
Util2d:rech_237: resetting 'how' to external
Util2d:rech_238: resetting 'how' to external
Util2d:rech_239: resetting 'how' to external
Util2d:rech_240: resetting 'how' to external
Util2d:rech_241: resetting 'how' to external
Util2d:rech_242: resetting 'how' to external
Util2d:rech_243: resetting 'how' to external
Util2d:rech_244: resetting 'how' to external
Util2d:rech_245: resetting 'how' to external
Util2d:rech_246: resetting 'how' to external
Util2d:rech_247: resetting 'how' to external
Util2d:rech_248: resetting 'how' to external
Util2d:rech_249: resetting 'how' to external
Util2d:rech_250: resetting 'how' to external
Util2d:rech_251: resetting 'how' to external
Util2d:rech_252: resetting 'how' to external
Util2d:rech_253: resetting 'how' to external
Util2d:rech_254: resetting 'how' to external
Util2d:rech_255: resetting 'how' to external
Util2d:rech_256: resetting 'how' to external
Util2d:rech_257: resetting 'how' to external
Util2d:rech_258: resetting 'how' to external
Util2d:rech_259: resetting 'how' to external
Util2d:rech_260: resetting 'how' to external
Util2d:rech_261: resetting 'how' to external
Util2d:rech_262: resetting 'how' to external
Util2d:rech_263: resetting 'how' to external
Util2d:rech_264: resetting 'how' to external
Util2d:rech_265: resetting 'how' to external
Util2d:rech_266: resetting 'how' to external
Util2d:rech_267: resetting 'how' to external
Util2d:rech_268: resetting 'how' to external
Util2d:rech_269: resetting 'how' to external
Util2d:rech_270: resetting 'how' to external
Util2d:rech_271: resetting 'how' to external
Util2d:rech_272: resetting 'how' to external
Util2d:rech_273: resetting 'how' to external
Util2d:rech_274: resetting 'how' to external
Util2d:rech_275: resetting 'how' to external
Util2d:rech_276: resetting 'how' to external
Util2d:rech_277: resetting 'how' to external
Util2d:rech_278: resetting 'how' to external
Util2d:rech_279: resetting 'how' to external
Util2d:rech_280: resetting 'how' to external
Util2d:rech_281: resetting 'how' to external
Util2d:rech_282: resetting 'how' to external
Util2d:rech_283: resetting 'how' to external
Util2d:rech_284: resetting 'how' to external
Util2d:rech_285: resetting 'how' to external
Util2d:rech_286: resetting 'how' to external
Util2d:rech_287: resetting 'how' to external
Util2d:rech_288: resetting 'how' to external
Util2d:rech_289: resetting 'how' to external
Util2d:rech_290: resetting 'how' to external
Util2d:rech_291: resetting 'how' to external
Util2d:rech_292: resetting 'how' to external
Util2d:rech_293: resetting 'how' to external
Util2d:rech_294: resetting 'how' to external
Util2d:rech_295: resetting 'how' to external
Util2d:rech_296: resetting 'how' to external
Util2d:rech_297: resetting 'how' to external
Util2d:rech_298: resetting 'how' to external
Util2d:rech_299: resetting 'how' to external
Util2d:rech_300: resetting 'how' to external
Util2d:rech_301: resetting 'how' to external
Util2d:rech_302: resetting 'how' to external
Util2d:rech_303: resetting 'how' to external
Util2d:rech_304: resetting 'how' to external
Util2d:rech_305: resetting 'how' to external
Util2d:rech_306: resetting 'how' to external
Util2d:rech_307: resetting 'how' to external
Util2d:rech_308: resetting 'how' to external
Util2d:rech_309: resetting 'how' to external
Util2d:rech_310: resetting 'how' to external
Util2d:rech_311: resetting 'how' to external
Util2d:rech_312: resetting 'how' to external
Util2d:rech_313: resetting 'how' to external
Util2d:rech_314: resetting 'how' to external
Util2d:rech_315: resetting 'how' to external
Util2d:rech_316: resetting 'how' to external
Util2d:rech_317: resetting 'how' to external
Util2d:rech_318: resetting 'how' to external
Util2d:rech_319: resetting 'how' to external
Util2d:rech_320: resetting 'how' to external
Util2d:rech_321: resetting 'how' to external
Util2d:rech_322: resetting 'how' to external
Util2d:rech_323: resetting 'how' to external
Util2d:rech_324: resetting 'how' to external
Util2d:rech_325: resetting 'how' to external
Util2d:rech_326: resetting 'how' to external
Util2d:rech_327: resetting 'how' to external
Util2d:rech_328: resetting 'how' to external
Util2d:rech_329: resetting 'how' to external
Util2d:rech_330: resetting 'how' to external
Util2d:rech_331: resetting 'how' to external
Util2d:rech_332: resetting 'how' to external
Util2d:rech_333: resetting 'how' to external
Util2d:rech_334: resetting 'how' to external
Util2d:rech_335: resetting 'how' to external
Util2d:rech_336: resetting 'how' to external
Util2d:rech_337: resetting 'how' to external
Util2d:rech_338: resetting 'how' to external
Util2d:rech_339: resetting 'how' to external
Util2d:rech_340: resetting 'how' to external
Util2d:rech_341: resetting 'how' to external
Util2d:rech_342: resetting 'how' to external
Util2d:rech_343: resetting 'how' to external
Util2d:rech_344: resetting 'how' to external
Util2d:rech_345: resetting 'how' to external
Util2d:rech_346: resetting 'how' to external
Util2d:rech_347: resetting 'how' to external
Util2d:rech_348: resetting 'how' to external
Util2d:rech_349: resetting 'how' to external
Util2d:rech_350: resetting 'how' to external
Util2d:rech_351: resetting 'how' to external
Util2d:rech_352: resetting 'how' to external
Util2d:rech_353: resetting 'how' to external
Util2d:rech_354: resetting 'how' to external
Util2d:rech_355: resetting 'how' to external
Util2d:rech_356: resetting 'how' to external
Util2d:rech_357: resetting 'how' to external
Util2d:rech_358: resetting 'how' to external
Util2d:rech_359: resetting 'how' to external
Util2d:rech_360: resetting 'how' to external
Util2d:rech_361: resetting 'how' to external
Util2d:rech_362: resetting 'how' to external
Util2d:rech_363: resetting 'how' to external
Util2d:rech_364: resetting 'how' to external
Util2d:rech_365: resetting 'how' to external
Util2d:rech_366: resetting 'how' to external
Util2d:rech_367: resetting 'how' to external
Util2d:rech_368: resetting 'how' to external
Util2d:rech_369: resetting 'how' to external
Util2d:rech_370: resetting 'how' to external
Util2d:rech_371: resetting 'how' to external
Util2d:rech_372: resetting 'how' to external
Util2d:rech_373: resetting 'how' to external
Util2d:rech_374: resetting 'how' to external
Util2d:rech_375: resetting 'how' to external
Util2d:rech_376: resetting 'how' to external
Util2d:rech_377: resetting 'how' to external
Util2d:rech_378: resetting 'how' to external
Util2d:rech_379: resetting 'how' to external
Util2d:rech_380: resetting 'how' to external
Util2d:rech_381: resetting 'how' to external
Util2d:rech_382: resetting 'how' to external
Util2d:rech_383: resetting 'how' to external
Util2d:rech_384: resetting 'how' to external
Util2d:rech_385: resetting 'how' to external
Util2d:rech_386: resetting 'how' to external
Util2d:rech_387: resetting 'how' to external
Util2d:rech_388: resetting 'how' to external
Util2d:rech_389: resetting 'how' to external
Util2d:rech_390: resetting 'how' to external
Util2d:rech_391: resetting 'how' to external
Util2d:rech_392: resetting 'how' to external
Util2d:rech_393: resetting 'how' to external
Util2d:rech_394: resetting 'how' to external
Util2d:rech_395: resetting 'how' to external
Util2d:rech_396: resetting 'how' to external
Util2d:rech_397: resetting 'how' to external
Util2d:rech_398: resetting 'how' to external
Util2d:rech_399: resetting 'how' to external
Util2d:rech_400: resetting 'how' to external
Util2d:rech_401: resetting 'how' to external
Util2d:rech_402: resetting 'how' to external
Util2d:rech_403: resetting 'how' to external
Util2d:rech_404: resetting 'how' to external
Util2d:rech_405: resetting 'how' to external
Util2d:rech_406: resetting 'how' to external
Util2d:rech_407: resetting 'how' to external
Util2d:rech_408: resetting 'how' to external
Util2d:rech_409: resetting 'how' to external
Util2d:rech_410: resetting 'how' to external
Util2d:rech_411: resetting 'how' to external
Util2d:rech_412: resetting 'how' to external
Util2d:rech_413: resetting 'how' to external
Util2d:rech_414: resetting 'how' to external
Util2d:rech_415: resetting 'how' to external
Util2d:rech_416: resetting 'how' to external
Util2d:rech_417: resetting 'how' to external
Util2d:rech_418: resetting 'how' to external
Util2d:rech_419: resetting 'how' to external
Util2d:rech_420: resetting 'how' to external
Util2d:rech_421: resetting 'how' to external
Util2d:rech_422: resetting 'how' to external
Util2d:rech_423: resetting 'how' to external
Util2d:rech_424: resetting 'how' to external
Util2d:rech_425: resetting 'how' to external
Util2d:rech_426: resetting 'how' to external
Util2d:rech_427: resetting 'how' to external
Util2d:rech_428: resetting 'how' to external
Util2d:rech_429: resetting 'how' to external
Util2d:rech_430: resetting 'how' to external
Util2d:rech_431: resetting 'how' to external
Util2d:rech_432: resetting 'how' to external
Util2d:rech_433: resetting 'how' to external
Util2d:rech_434: resetting 'how' to external
Util2d:rech_435: resetting 'how' to external
Util2d:rech_436: resetting 'how' to external
Util2d:rech_437: resetting 'how' to external
Util2d:rech_438: resetting 'how' to external
Util2d:rech_439: resetting 'how' to external
Util2d:rech_440: resetting 'how' to external
Util2d:rech_441: resetting 'how' to external
Util2d:rech_442: resetting 'how' to external
Util2d:rech_443: resetting 'how' to external
Util2d:rech_444: resetting 'how' to external
Util2d:rech_445: resetting 'how' to external
Util2d:rech_446: resetting 'how' to external
Util2d:rech_447: resetting 'how' to external
Util2d:rech_448: resetting 'how' to external
Util2d:rech_449: resetting 'how' to external
Util2d:rech_450: resetting 'how' to external
Util2d:rech_451: resetting 'how' to external
Util2d:rech_452: resetting 'how' to external
Util2d:rech_453: resetting 'how' to external
Util2d:rech_454: resetting 'how' to external
Util2d:rech_455: resetting 'how' to external
Util2d:rech_456: resetting 'how' to external
Util2d:rech_457: resetting 'how' to external
Util2d:rech_458: resetting 'how' to external
Util2d:rech_459: resetting 'how' to external
Util2d:rech_460: resetting 'how' to external
Util2d:rech_461: resetting 'how' to external
Util2d:rech_462: resetting 'how' to external
Util2d:rech_463: resetting 'how' to external
Util2d:rech_464: resetting 'how' to external
Util2d:rech_465: resetting 'how' to external
Util2d:rech_466: resetting 'how' to external
Util2d:rech_467: resetting 'how' to external
Util2d:rech_468: resetting 'how' to external
Util2d:rech_469: resetting 'how' to external
Util2d:rech_470: resetting 'how' to external
Util2d:rech_471: resetting 'how' to external
Util2d:rech_472: resetting 'how' to external
Util2d:rech_473: resetting 'how' to external
Util2d:rech_474: resetting 'how' to external
Util2d:rech_475: resetting 'how' to external
Util2d:rech_476: resetting 'how' to external
Util2d:rech_477: resetting 'how' to external
Util2d:rech_478: resetting 'how' to external
Util2d:rech_479: resetting 'how' to external
Util2d:rech_480: resetting 'how' to external
Util2d:rech_481: resetting 'how' to external
Util2d:rech_482: resetting 'how' to external
Util2d:rech_483: resetting 'how' to external
Util2d:rech_484: resetting 'how' to external
Util2d:rech_485: resetting 'how' to external
Util2d:rech_486: resetting 'how' to external
Util2d:rech_487: resetting 'how' to external
Util2d:rech_488: resetting 'how' to external
Util2d:rech_489: resetting 'how' to external
Util2d:rech_490: resetting 'how' to external
Util2d:rech_491: resetting 'how' to external
Util2d:rech_492: resetting 'how' to external
Util2d:rech_493: resetting 'how' to external
Util2d:rech_494: resetting 'how' to external
Util2d:rech_495: resetting 'how' to external
Util2d:rech_496: resetting 'how' to external
Util2d:rech_497: resetting 'how' to external
Util2d:rech_498: resetting 'how' to external
Util2d:rech_499: resetting 'how' to external
Util2d:rech_500: resetting 'how' to external
Util2d:rech_501: resetting 'how' to external
Util2d:rech_502: resetting 'how' to external
Util2d:rech_503: resetting 'how' to external
Util2d:rech_504: resetting 'how' to external
Util2d:rech_505: resetting 'how' to external
Util2d:rech_506: resetting 'how' to external
Util2d:rech_507: resetting 'how' to external
Util2d:rech_508: resetting 'how' to external
Util2d:rech_509: resetting 'how' to external
Util2d:rech_510: resetting 'how' to external
Util2d:rech_511: resetting 'how' to external
Util2d:rech_512: resetting 'how' to external
Util2d:rech_513: resetting 'how' to external
Util2d:rech_514: resetting 'how' to external
Util2d:rech_515: resetting 'how' to external
Util2d:rech_516: resetting 'how' to external
Util2d:rech_517: resetting 'how' to external
Util2d:rech_518: resetting 'how' to external
Util2d:rech_519: resetting 'how' to external
Util2d:rech_520: resetting 'how' to external
Util2d:rech_521: resetting 'how' to external
Util2d:rech_522: resetting 'how' to external
Util2d:rech_523: resetting 'how' to external
Util2d:rech_524: resetting 'how' to external
Util2d:rech_525: resetting 'how' to external
Util2d:rech_526: resetting 'how' to external
Util2d:rech_527: resetting 'how' to external
Util2d:rech_528: resetting 'how' to external
Util2d:rech_529: resetting 'how' to external
Util2d:rech_530: resetting 'how' to external
Util2d:rech_531: resetting 'how' to external
Util2d:rech_532: resetting 'how' to external
Util2d:rech_533: resetting 'how' to external
Util2d:rech_534: resetting 'how' to external
Util2d:rech_535: resetting 'how' to external
Util2d:rech_536: resetting 'how' to external
Util2d:rech_537: resetting 'how' to external
Util2d:rech_538: resetting 'how' to external
Util2d:rech_539: resetting 'how' to external
Util2d:rech_540: resetting 'how' to external
Util2d:rech_541: resetting 'how' to external
Util2d:rech_542: resetting 'how' to external
Util2d:rech_543: resetting 'how' to external
Util2d:rech_544: resetting 'how' to external
Util2d:rech_545: resetting 'how' to external
Util2d:rech_546: resetting 'how' to external
Util2d:rech_547: resetting 'how' to external
Util2d:rech_548: resetting 'how' to external
Util2d:rech_549: resetting 'how' to external
Util2d:rech_550: resetting 'how' to external
Util2d:rech_551: resetting 'how' to external
Util2d:rech_552: resetting 'how' to external
Util2d:rech_553: resetting 'how' to external
Util2d:rech_554: resetting 'how' to external
Util2d:rech_555: resetting 'how' to external
Util2d:rech_556: resetting 'how' to external
Util2d:rech_557: resetting 'how' to external
Util2d:rech_558: resetting 'how' to external
Util2d:rech_559: resetting 'how' to external
Util2d:rech_560: resetting 'how' to external
Util2d:rech_561: resetting 'how' to external
Util2d:rech_562: resetting 'how' to external
Util2d:rech_563: resetting 'how' to external
Util2d:rech_564: resetting 'how' to external
Util2d:rech_565: resetting 'how' to external
Util2d:rech_566: resetting 'how' to external
Util2d:rech_567: resetting 'how' to external
Util2d:rech_568: resetting 'how' to external
Util2d:rech_569: resetting 'how' to external
Util2d:rech_570: resetting 'how' to external
Util2d:rech_571: resetting 'how' to external
Util2d:rech_572: resetting 'how' to external
Util2d:rech_573: resetting 'how' to external
Util2d:rech_574: resetting 'how' to external
Util2d:rech_575: resetting 'how' to external
Util2d:rech_576: resetting 'how' to external
Util2d:rech_577: resetting 'how' to external
Util2d:rech_578: resetting 'how' to external
Util2d:rech_579: resetting 'how' to external
Util2d:rech_580: resetting 'how' to external
Util2d:rech_581: resetting 'how' to external
Util2d:rech_582: resetting 'how' to external
Util2d:rech_583: resetting 'how' to external
Util2d:rech_584: resetting 'how' to external
Util2d:rech_585: resetting 'how' to external
Util2d:rech_586: resetting 'how' to external
Util2d:rech_587: resetting 'how' to external
Util2d:rech_588: resetting 'how' to external
Util2d:rech_589: resetting 'how' to external
Util2d:rech_590: resetting 'how' to external
Util2d:rech_591: resetting 'how' to external
Util2d:rech_592: resetting 'how' to external
Util2d:rech_593: resetting 'how' to external
Util2d:rech_594: resetting 'how' to external
Util2d:rech_595: resetting 'how' to external
Util2d:rech_596: resetting 'how' to external
Util2d:rech_597: resetting 'how' to external
Util2d:rech_598: resetting 'how' to external
Util2d:rech_599: resetting 'how' to external
Util2d:rech_600: resetting 'how' to external
Util2d:rech_601: resetting 'how' to external
Util2d:rech_602: resetting 'how' to external
Util2d:rech_603: resetting 'how' to external
Util2d:rech_604: resetting 'how' to external
Util2d:rech_605: resetting 'how' to external
Util2d:rech_606: resetting 'how' to external
Util2d:rech_607: resetting 'how' to external
Util2d:rech_608: resetting 'how' to external
Util2d:rech_609: resetting 'how' to external
Util2d:rech_610: resetting 'how' to external
Util2d:rech_611: resetting 'how' to external
Util2d:rech_612: resetting 'how' to external
Util2d:rech_613: resetting 'how' to external
Util2d:rech_614: resetting 'how' to external
Util2d:rech_615: resetting 'how' to external
Util2d:rech_616: resetting 'how' to external
Util2d:rech_617: resetting 'how' to external
Util2d:rech_618: resetting 'how' to external
Util2d:rech_619: resetting 'how' to external
Util2d:rech_620: resetting 'how' to external
Util2d:rech_621: resetting 'how' to external
Util2d:rech_622: resetting 'how' to external
Util2d:rech_623: resetting 'how' to external
Util2d:rech_624: resetting 'how' to external
Util2d:rech_625: resetting 'how' to external
Util2d:rech_626: resetting 'how' to external
Util2d:rech_627: resetting 'how' to external
Util2d:rech_628: resetting 'how' to external
Util2d:rech_629: resetting 'how' to external
Util2d:rech_630: resetting 'how' to external
Util2d:rech_631: resetting 'how' to external
Util2d:rech_632: resetting 'how' to external
Util2d:rech_633: resetting 'how' to external
Util2d:rech_634: resetting 'how' to external
Util2d:rech_635: resetting 'how' to external
Util2d:rech_636: resetting 'how' to external
Util2d:rech_637: resetting 'how' to external
Util2d:rech_638: resetting 'how' to external
Util2d:rech_639: resetting 'how' to external
Util2d:rech_640: resetting 'how' to external
Util2d:rech_641: resetting 'how' to external
Util2d:rech_642: resetting 'how' to external
Util2d:rech_643: resetting 'how' to external
Util2d:rech_644: resetting 'how' to external
Util2d:rech_645: resetting 'how' to external
Util2d:rech_646: resetting 'how' to external
Util2d:rech_647: resetting 'how' to external
Util2d:rech_648: resetting 'how' to external
Util2d:rech_649: resetting 'how' to external
Util2d:rech_650: resetting 'how' to external
Util2d:rech_651: resetting 'how' to external
Util2d:rech_652: resetting 'how' to external
Util2d:rech_653: resetting 'how' to external
Util2d:rech_654: resetting 'how' to external
Util2d:rech_655: resetting 'how' to external
Util2d:rech_656: resetting 'how' to external
Util2d:rech_657: resetting 'how' to external
Util2d:rech_658: resetting 'how' to external
Util2d:rech_659: resetting 'how' to external
Util2d:rech_660: resetting 'how' to external
Util2d:rech_661: resetting 'how' to external
Util2d:rech_662: resetting 'how' to external
Util2d:rech_663: resetting 'how' to external
Util2d:rech_664: resetting 'how' to external
Util2d:rech_665: resetting 'how' to external
Util2d:rech_666: resetting 'how' to external
Util2d:rech_667: resetting 'how' to external
Util2d:rech_668: resetting 'how' to external
Util2d:rech_669: resetting 'how' to external
Util2d:rech_670: resetting 'how' to external
Util2d:rech_671: resetting 'how' to external
Util2d:rech_672: resetting 'how' to external
Util2d:rech_673: resetting 'how' to external
Util2d:rech_674: resetting 'how' to external
Util2d:rech_675: resetting 'how' to external
Util2d:rech_676: resetting 'how' to external
Util2d:rech_677: resetting 'how' to external
Util2d:rech_678: resetting 'how' to external
Util2d:rech_679: resetting 'how' to external
Util2d:rech_680: resetting 'how' to external
Util2d:rech_681: resetting 'how' to external
Util2d:rech_682: resetting 'how' to external
Util2d:rech_683: resetting 'how' to external
Util2d:rech_684: resetting 'how' to external
Util2d:rech_685: resetting 'how' to external
Util2d:rech_686: resetting 'how' to external
Util2d:rech_687: resetting 'how' to external
Util2d:rech_688: resetting 'how' to external
Util2d:rech_689: resetting 'how' to external
Util2d:rech_690: resetting 'how' to external
Util2d:rech_691: resetting 'how' to external
Util2d:rech_692: resetting 'how' to external
Util2d:rech_693: resetting 'how' to external
Util2d:rech_694: resetting 'how' to external
Util2d:rech_695: resetting 'how' to external
Util2d:rech_696: resetting 'how' to external
Util2d:rech_697: resetting 'how' to external
Util2d:rech_698: resetting 'how' to external
Util2d:rech_699: resetting 'how' to external
Util2d:rech_700: resetting 'how' to external
Util2d:rech_701: resetting 'how' to external
Util2d:rech_702: resetting 'how' to external
Util2d:rech_703: resetting 'how' to external
Util2d:rech_704: resetting 'how' to external
Util2d:rech_705: resetting 'how' to external
Util2d:rech_706: resetting 'how' to external
Util2d:rech_707: resetting 'how' to external
Util2d:rech_708: resetting 'how' to external
Util2d:rech_709: resetting 'how' to external
Util2d:rech_710: resetting 'how' to external
Util2d:rech_711: resetting 'how' to external
Util2d:rech_712: resetting 'how' to external
Util2d:rech_713: resetting 'how' to external
Util2d:rech_714: resetting 'how' to external
Util2d:rech_715: resetting 'how' to external
Util2d:rech_716: resetting 'how' to external
Util2d:rech_717: resetting 'how' to external
Util2d:rech_718: resetting 'how' to external
Util2d:rech_719: resetting 'how' to external
Util2d:rech_720: resetting 'how' to external
Util2d:rech_721: resetting 'how' to external
Util2d:rech_722: resetting 'how' to external
Util2d:rech_723: resetting 'how' to external
Util2d:rech_724: resetting 'how' to external
Util2d:rech_725: resetting 'how' to external
Util2d:rech_726: resetting 'how' to external
Util2d:rech_727: resetting 'how' to external
Util2d:rech_728: resetting 'how' to external
Util2d:rech_729: resetting 'how' to external
Util2d:rech_730: resetting 'how' to external
Util2d:rech_731: resetting 'how' to external
Util2d:rech_732: resetting 'how' to external
Util2d:rech_733: resetting 'how' to external
Util2d:rech_734: resetting 'how' to external
Util2d:rech_735: resetting 'how' to external
Util2d:rech_736: resetting 'how' to external
Util2d:rech_737: resetting 'how' to external
Util2d:rech_738: resetting 'how' to external
Util2d:rech_739: resetting 'how' to external
Util2d:rech_740: resetting 'how' to external
Util2d:rech_741: resetting 'how' to external
Util2d:rech_742: resetting 'how' to external
Util2d:rech_743: resetting 'how' to external
Util2d:rech_744: resetting 'how' to external
Util2d:rech_745: resetting 'how' to external
Util2d:rech_746: resetting 'how' to external
Util2d:rech_747: resetting 'how' to external
Util2d:rech_748: resetting 'how' to external
Util2d:rech_749: resetting 'how' to external
Util2d:rech_750: resetting 'how' to external
Util2d:rech_751: resetting 'how' to external
Util2d:rech_752: resetting 'how' to external
Util2d:rech_753: resetting 'how' to external
Util2d:rech_754: resetting 'how' to external
Util2d:rech_755: resetting 'how' to external
Util2d:rech_756: resetting 'how' to external
Util2d:rech_757: resetting 'how' to external
Util2d:rech_758: resetting 'how' to external
Util2d:rech_759: resetting 'how' to external
Util2d:rech_760: resetting 'how' to external
Util2d:rech_761: resetting 'how' to external
Util2d:rech_762: resetting 'how' to external
Util2d:rech_763: resetting 'how' to external
Util2d:rech_764: resetting 'how' to external
Util2d:rech_765: resetting 'how' to external
Util2d:rech_766: resetting 'how' to external
Util2d:rech_767: resetting 'how' to external
Util2d:rech_768: resetting 'how' to external
Util2d:rech_769: resetting 'how' to external
Util2d:rech_770: resetting 'how' to external
Util2d:rech_771: resetting 'how' to external
Util2d:rech_772: resetting 'how' to external
Util2d:rech_773: resetting 'how' to external
Util2d:rech_774: resetting 'how' to external
Util2d:rech_775: resetting 'how' to external
Util2d:rech_776: resetting 'how' to external
Util2d:rech_777: resetting 'how' to external
Util2d:rech_778: resetting 'how' to external
Util2d:rech_779: resetting 'how' to external
Util2d:rech_780: resetting 'how' to external
Util2d:rech_781: resetting 'how' to external
Util2d:rech_782: resetting 'how' to external
Util2d:rech_783: resetting 'how' to external
Util2d:rech_784: resetting 'how' to external
Util2d:rech_785: resetting 'how' to external
Util2d:rech_786: resetting 'how' to external
Util2d:rech_787: resetting 'how' to external
Util2d:rech_788: resetting 'how' to external
Util2d:rech_789: resetting 'how' to external
Util2d:rech_790: resetting 'how' to external
Util2d:rech_791: resetting 'how' to external
Util2d:rech_792: resetting 'how' to external
Util2d:rech_793: resetting 'how' to external
Util2d:rech_794: resetting 'how' to external
Util2d:rech_795: resetting 'how' to external
Util2d:rech_796: resetting 'how' to external
Util2d:rech_797: resetting 'how' to external
Util2d:rech_798: resetting 'how' to external
Util2d:rech_799: resetting 'how' to external
Util2d:rech_800: resetting 'how' to external
Util2d:rech_801: resetting 'how' to external
Util2d:rech_802: resetting 'how' to external
Util2d:rech_803: resetting 'how' to external
Util2d:rech_804: resetting 'how' to external
Util2d:rech_805: resetting 'how' to external
Util2d:rech_806: resetting 'how' to external
Util2d:rech_807: resetting 'how' to external
Util2d:rech_808: resetting 'how' to external
Util2d:rech_809: resetting 'how' to external
Util2d:rech_810: resetting 'how' to external
Util2d:rech_811: resetting 'how' to external
Util2d:rech_812: resetting 'how' to external
Util2d:rech_813: resetting 'how' to external
Util2d:rech_814: resetting 'how' to external
Util2d:rech_815: resetting 'how' to external
Util2d:rech_816: resetting 'how' to external
Util2d:rech_817: resetting 'how' to external
Util2d:rech_818: resetting 'how' to external
Util2d:rech_819: resetting 'how' to external
Util2d:rech_820: resetting 'how' to external
Util2d:rech_821: resetting 'how' to external
Util2d:rech_822: resetting 'how' to external
Util2d:rech_823: resetting 'how' to external
Util2d:rech_824: resetting 'how' to external
Util2d:rech_825: resetting 'how' to external
Util2d:rech_826: resetting 'how' to external
Util2d:rech_827: resetting 'how' to external
Util2d:rech_828: resetting 'how' to external
Util2d:rech_829: resetting 'how' to external
Util2d:rech_830: resetting 'how' to external
Util2d:rech_831: resetting 'how' to external
Util2d:rech_832: resetting 'how' to external
Util2d:rech_833: resetting 'how' to external
Util2d:rech_834: resetting 'how' to external
Util2d:rech_835: resetting 'how' to external
Util2d:rech_836: resetting 'how' to external
Util2d:rech_837: resetting 'how' to external
Util2d:rech_838: resetting 'how' to external
Util2d:rech_839: resetting 'how' to external
Util2d:rech_840: resetting 'how' to external
Util2d:rech_841: resetting 'how' to external
Util2d:rech_842: resetting 'how' to external
Util2d:rech_843: resetting 'how' to external
Util2d:rech_844: resetting 'how' to external
Util2d:rech_845: resetting 'how' to external
Util2d:rech_846: resetting 'how' to external
Util2d:rech_847: resetting 'how' to external
Util2d:rech_848: resetting 'how' to external
Util2d:rech_849: resetting 'how' to external
Util2d:rech_850: resetting 'how' to external
Util2d:rech_851: resetting 'how' to external
Util2d:rech_852: resetting 'how' to external
Util2d:rech_853: resetting 'how' to external
Util2d:rech_854: resetting 'how' to external
Util2d:rech_855: resetting 'how' to external
Util2d:rech_856: resetting 'how' to external
Util2d:rech_857: resetting 'how' to external
Util2d:rech_858: resetting 'how' to external
Util2d:rech_859: resetting 'how' to external
Util2d:rech_860: resetting 'how' to external
Util2d:rech_861: resetting 'how' to external
Util2d:rech_862: resetting 'how' to external
Util2d:rech_863: resetting 'how' to external
Util2d:rech_864: resetting 'how' to external
Util2d:rech_865: resetting 'how' to external
Util2d:rech_866: resetting 'how' to external
Util2d:rech_867: resetting 'how' to external
Util2d:rech_868: resetting 'how' to external
Util2d:rech_869: resetting 'how' to external
Util2d:rech_870: resetting 'how' to external
Util2d:rech_871: resetting 'how' to external
Util2d:rech_872: resetting 'how' to external
Util2d:rech_873: resetting 'how' to external
Util2d:rech_874: resetting 'how' to external
Util2d:rech_875: resetting 'how' to external
Util2d:rech_876: resetting 'how' to external
Util2d:rech_877: resetting 'how' to external
Util2d:rech_878: resetting 'how' to external
Util2d:rech_879: resetting 'how' to external
Util2d:rech_880: resetting 'how' to external
Util2d:rech_881: resetting 'how' to external
Util2d:rech_882: resetting 'how' to external
Util2d:rech_883: resetting 'how' to external
Util2d:rech_884: resetting 'how' to external
Util2d:rech_885: resetting 'how' to external
Util2d:rech_886: resetting 'how' to external
Util2d:rech_887: resetting 'how' to external
Util2d:rech_888: resetting 'how' to external
Util2d:rech_889: resetting 'how' to external
Util2d:rech_890: resetting 'how' to external
Util2d:rech_891: resetting 'how' to external
Util2d:rech_892: resetting 'how' to external
Util2d:rech_893: resetting 'how' to external
Util2d:rech_894: resetting 'how' to external
Util2d:rech_895: resetting 'how' to external
Util2d:rech_896: resetting 'how' to external
Util2d:rech_897: resetting 'how' to external
Util2d:rech_898: resetting 'how' to external
Util2d:rech_899: resetting 'how' to external
Util2d:rech_900: resetting 'how' to external
Util2d:rech_901: resetting 'how' to external
Util2d:rech_902: resetting 'how' to external
Util2d:rech_903: resetting 'how' to external
Util2d:rech_904: resetting 'how' to external
Util2d:rech_905: resetting 'how' to external
Util2d:rech_906: resetting 'how' to external
Util2d:rech_907: resetting 'how' to external
Util2d:rech_908: resetting 'how' to external
Util2d:rech_909: resetting 'how' to external
Util2d:rech_910: resetting 'how' to external
Util2d:rech_911: resetting 'how' to external
Util2d:rech_912: resetting 'how' to external
Util2d:rech_913: resetting 'how' to external
Util2d:rech_914: resetting 'how' to external
Util2d:rech_915: resetting 'how' to external
Util2d:rech_916: resetting 'how' to external
Util2d:rech_917: resetting 'how' to external
Util2d:rech_918: resetting 'how' to external
Util2d:rech_919: resetting 'how' to external
Util2d:rech_920: resetting 'how' to external
Util2d:rech_921: resetting 'how' to external
Util2d:rech_922: resetting 'how' to external
Util2d:rech_923: resetting 'how' to external
Util2d:rech_924: resetting 'how' to external
Util2d:rech_925: resetting 'how' to external
Util2d:rech_926: resetting 'how' to external
Util2d:rech_927: resetting 'how' to external
Util2d:rech_928: resetting 'how' to external
Util2d:rech_929: resetting 'how' to external
Util2d:rech_930: resetting 'how' to external
Util2d:rech_931: resetting 'how' to external
Util2d:rech_932: resetting 'how' to external
Util2d:rech_933: resetting 'how' to external
Util2d:rech_934: resetting 'how' to external
Util2d:rech_935: resetting 'how' to external
Util2d:rech_936: resetting 'how' to external
Util2d:rech_937: resetting 'how' to external
Util2d:rech_938: resetting 'how' to external
Util2d:rech_939: resetting 'how' to external
Util2d:rech_940: resetting 'how' to external
Util2d:rech_941: resetting 'how' to external
Util2d:rech_942: resetting 'how' to external
Util2d:rech_943: resetting 'how' to external
Util2d:rech_944: resetting 'how' to external
Util2d:rech_945: resetting 'how' to external
Util2d:rech_946: resetting 'how' to external
Util2d:rech_947: resetting 'how' to external
Util2d:rech_948: resetting 'how' to external
Util2d:rech_949: resetting 'how' to external
Util2d:rech_950: resetting 'how' to external
Util2d:rech_951: resetting 'how' to external
Util2d:rech_952: resetting 'how' to external
Util2d:rech_953: resetting 'how' to external
Util2d:rech_954: resetting 'how' to external
Util2d:rech_955: resetting 'how' to external
Util2d:rech_956: resetting 'how' to external
Util2d:rech_957: resetting 'how' to external
Util2d:rech_958: resetting 'how' to external
Util2d:rech_959: resetting 'how' to external
Util2d:rech_960: resetting 'how' to external
Util2d:rech_961: resetting 'how' to external
Util2d:rech_962: resetting 'how' to external
Util2d:rech_963: resetting 'how' to external
Util2d:rech_964: resetting 'how' to external
Util2d:rech_965: resetting 'how' to external
Util2d:rech_966: resetting 'how' to external
Util2d:rech_967: resetting 'how' to external
Util2d:rech_968: resetting 'how' to external
Util2d:rech_969: resetting 'how' to external
Util2d:rech_970: resetting 'how' to external
Util2d:rech_971: resetting 'how' to external
Util2d:rech_972: resetting 'how' to external
Util2d:rech_973: resetting 'how' to external
Util2d:rech_974: resetting 'how' to external
Util2d:rech_975: resetting 'how' to external
Util2d:rech_976: resetting 'how' to external
Util2d:rech_977: resetting 'how' to external
Util2d:rech_978: resetting 'how' to external
Util2d:rech_979: resetting 'how' to external
Util2d:rech_980: resetting 'how' to external
Util2d:rech_981: resetting 'how' to external
Util2d:rech_982: resetting 'how' to external
Util2d:rech_983: resetting 'how' to external
Util2d:rech_984: resetting 'how' to external
Util2d:rech_985: resetting 'how' to external
Util2d:rech_986: resetting 'how' to external
Util2d:rech_987: resetting 'how' to external
Util2d:rech_988: resetting 'how' to external
Util2d:rech_989: resetting 'how' to external
Util2d:rech_990: resetting 'how' to external
Util2d:rech_991: resetting 'how' to external
Util2d:rech_992: resetting 'how' to external
Util2d:rech_993: resetting 'how' to external
Util2d:rech_994: resetting 'how' to external
Util2d:rech_995: resetting 'how' to external
Util2d:rech_996: resetting 'how' to external
Util2d:rech_997: resetting 'how' to external
Util2d:rech_998: resetting 'how' to external
Util2d:rech_999: resetting 'how' to external
Util2d:rech_1000: resetting 'how' to external
Util2d:rech_1001: resetting 'how' to external
Util2d:rech_1002: resetting 'how' to external
Util2d:rech_1003: resetting 'how' to external
Util2d:rech_1004: resetting 'how' to external
Util2d:rech_1005: resetting 'how' to external
Util2d:rech_1006: resetting 'how' to external
Util2d:rech_1007: resetting 'how' to external
Util2d:rech_1008: resetting 'how' to external
Util2d:rech_1009: resetting 'how' to external
Util2d:rech_1010: resetting 'how' to external
Util2d:rech_1011: resetting 'how' to external
Util2d:rech_1012: resetting 'how' to external
Util2d:rech_1013: resetting 'how' to external
Util2d:rech_1014: resetting 'how' to external
Util2d:rech_1015: resetting 'how' to external
Util2d:rech_1016: resetting 'how' to external
Util2d:rech_1017: resetting 'how' to external
Util2d:rech_1018: resetting 'how' to external
Util2d:rech_1019: resetting 'how' to external
Util2d:rech_1020: resetting 'how' to external
Util2d:rech_1021: resetting 'how' to external
Util2d:rech_1022: resetting 'how' to external
Util2d:rech_1023: resetting 'how' to external
Util2d:rech_1024: resetting 'how' to external
Util2d:rech_1025: resetting 'how' to external
Util2d:rech_1026: resetting 'how' to external
Util2d:rech_1027: resetting 'how' to external
Util2d:rech_1028: resetting 'how' to external
Util2d:rech_1029: resetting 'how' to external
Util2d:rech_1030: resetting 'how' to external
Util2d:rech_1031: resetting 'how' to external
Util2d:rech_1032: resetting 'how' to external
Util2d:rech_1033: resetting 'how' to external
Util2d:rech_1034: resetting 'how' to external
Util2d:rech_1035: resetting 'how' to external
Util2d:rech_1036: resetting 'how' to external
Util2d:rech_1037: resetting 'how' to external
Util2d:rech_1038: resetting 'how' to external
Util2d:rech_1039: resetting 'how' to external
Util2d:rech_1040: resetting 'how' to external
Util2d:rech_1041: resetting 'how' to external
Util2d:rech_1042: resetting 'how' to external
Util2d:rech_1043: resetting 'how' to external
Util2d:rech_1044: resetting 'how' to external
Util2d:rech_1045: resetting 'how' to external
Util2d:rech_1046: resetting 'how' to external
Util2d:rech_1047: resetting 'how' to external
Util2d:rech_1048: resetting 'how' to external
Util2d:rech_1049: resetting 'how' to external
Util2d:rech_1050: resetting 'how' to external
Util2d:rech_1051: resetting 'how' to external
Util2d:rech_1052: resetting 'how' to external
Util2d:rech_1053: resetting 'how' to external
Util2d:rech_1054: resetting 'how' to external
Util2d:rech_1055: resetting 'how' to external
Util2d:rech_1056: resetting 'how' to external
Util2d:rech_1057: resetting 'how' to external
Util2d:rech_1058: resetting 'how' to external
Util2d:rech_1059: resetting 'how' to external
Util2d:rech_1060: resetting 'how' to external
Util2d:rech_1061: resetting 'how' to external
Util2d:rech_1062: resetting 'how' to external
Util2d:rech_1063: resetting 'how' to external
Util2d:rech_1064: resetting 'how' to external
Util2d:rech_1065: resetting 'how' to external
Util2d:rech_1066: resetting 'how' to external
Util2d:rech_1067: resetting 'how' to external
Util2d:rech_1068: resetting 'how' to external
Util2d:rech_1069: resetting 'how' to external
Util2d:rech_1070: resetting 'how' to external
Util2d:rech_1071: resetting 'how' to external
Util2d:rech_1072: resetting 'how' to external
Util2d:rech_1073: resetting 'how' to external
Util2d:rech_1074: resetting 'how' to external
Util2d:rech_1075: resetting 'how' to external
Util2d:rech_1076: resetting 'how' to external
Util2d:rech_1077: resetting 'how' to external
Util2d:rech_1078: resetting 'how' to external
Util2d:rech_1079: resetting 'how' to external
Util2d:rech_1080: resetting 'how' to external
Util2d:rech_1081: resetting 'how' to external
Util2d:rech_1082: resetting 'how' to external
Util2d:rech_1083: resetting 'how' to external
Util2d:rech_1084: resetting 'how' to external
Util2d:rech_1085: resetting 'how' to external
Util2d:rech_1086: resetting 'how' to external
Util2d:rech_1087: resetting 'how' to external
Util2d:rech_1088: resetting 'how' to external
Util2d:rech_1089: resetting 'how' to external
Util2d:rech_1090: resetting 'how' to external
Util2d:rech_1091: resetting 'how' to external
Util2d:rech_1092: resetting 'how' to external
Util2d:rech_1093: resetting 'how' to external
Util2d:rech_1094: resetting 'how' to external
Util2d:rech_1095: resetting 'how' to external
Util2d:rech_1096: resetting 'how' to external
Util2d:rech_1097: resetting 'how' to external
Util2d:hk: resetting 'how' to external
Util2d:vk: resetting 'how' to external
Util2d:ss: resetting 'how' to external
Util2d:sy: resetting 'how' to external
Util2d:hk: resetting 'how' to external
Util2d:vk: resetting 'how' to external
Util2d:ss: resetting 'how' to external
Util2d:sy: resetting 'how' to external
Util2d:hk: resetting 'how' to external
Util2d:vk: resetting 'how' to external
Util2d:ss: resetting 'how' to external
Util2d:sy: resetting 'how' to external

Now run the model.

[13]:

success, buff = ml.run_model(silent=True)
assert success, pformat(buff)

Inspecting outputs

First, let’s look at the list file. The list file summarizes the model’s results.

[14]:

mfl = flopy.utils.MfListBudget(model_ws / "freyberg.list")
df_flux, df_vol = mfl.get_dataframes(start_datetime="10-21-2015")
df_flux

[14]:

	STORAGE_IN	CONSTANT_HEAD_IN	WELLS_IN	DRAINS_IN	RECHARGE_IN	TOTAL_IN	STORAGE_OUT	CONSTANT_HEAD_OUT	WELLS_OUT	DRAINS_OUT	RECHARGE_OUT	TOTAL_OUT	IN-OUT	PERCENT_DISCREPANCY
2015-10-22	0.000000	0.000000	0.0	0.0	6276.861816	6276.861816	0.000000	2446.318848	0.000000	3830.650146	0.0	6276.968750	-0.106900	-0.0
2015-10-23	635.447998	0.000000	0.0	0.0	6428.198730	7063.646484	31.594000	2430.337891	1302.403198	3299.415039	0.0	7063.750000	-0.103500	-0.0
2015-10-24	1361.814941	0.000000	0.0	0.0	5397.295898	6759.110840	9.152200	2369.628174	1618.676392	2761.639648	0.0	6759.096191	0.014648	0.0
2015-10-25	677.577209	0.000000	0.0	0.0	5931.377441	6608.954590	180.233307	2353.585449	1498.694702	2576.461670	0.0	6608.975586	-0.020996	-0.0
2015-10-26	697.818298	0.000000	0.0	0.0	8378.572266	9076.390625	1051.897461	2417.248291	3119.501953	2487.737305	0.0	9076.384766	0.005859	0.0
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
2018-10-18	1298.433350	23.293699	0.0	0.0	4240.286133	5562.013184	1606.280396	1909.538574	1930.483154	115.724899	0.0	5562.026855	-0.013672	-0.0
2018-10-19	920.468689	25.997499	0.0	0.0	4082.749512	5029.215820	1659.194702	1892.605835	1279.166382	198.258896	0.0	5029.225586	-0.009766	-0.0
2018-10-20	496.671387	20.047001	0.0	0.0	5053.779297	5570.497559	2428.291016	1933.045898	794.582886	414.594513	0.0	5570.514648	-0.017090	-0.0
2018-10-21	230.320999	9.045700	0.0	0.0	6168.920410	6408.287109	2335.759521	2006.124268	1373.782593	692.638123	0.0	6408.304688	-0.017578	-0.0
2018-10-22	0.000000	317.783386	0.0	0.0	5190.390625	5508.173828	0.000000	745.363403	4762.799805	0.000000	0.0	5508.163086	0.010742	0.0

1097 rows × 14 columns

[15]:

groups = df_flux.groupby(lambda x: x.split("_")[-1], axis=1).groups
df_flux_in = df_flux.loc[:, groups["IN"]]
df_flux_in.columns = df_flux_in.columns.map(lambda x: x.split("_")[0])

df_flux_out = df_flux.loc[:, groups["OUT"]]
df_flux_out.columns = df_flux_out.columns.map(lambda x: x.split("_")[0])


df_flux_delta = df_flux_in - df_flux_out
df_flux_delta.iloc[-1, :].plot(kind="bar", figsize=(10, 10), grid=True)

/tmp/ipykernel_3840/2578515437.py:1: FutureWarning: DataFrame.groupby with axis=1 is deprecated. Do `frame.T.groupby(...)` without axis instead.
  groups = df_flux.groupby(lambda x: x.split("_")[-1], axis=1).groups

[15]:

<Axes: >

../_images/Notebooks_feat_working_stack_examples_23_2.png

Now let’s look at the simulated head.

[16]:

# if you pass the model instance, then the plots will be offset and rotated
h = flopy.utils.HeadFile(model_ws / "freyberg.hds", model=ml)
h.times

[16]:

[np.float32(1.0),
 np.float32(2.0),
 np.float32(3.0),
 np.float32(4.0),
 np.float32(5.0),
 np.float32(6.0),
 np.float32(7.0),
 np.float32(8.0),
 np.float32(9.0),
 np.float32(10.0),
 np.float32(11.0),
 np.float32(12.0),
 np.float32(13.0),
 np.float32(14.0),
 np.float32(15.0),
 np.float32(16.0),
 np.float32(17.0),
 np.float32(18.0),
 np.float32(19.0),
 np.float32(20.0),
 np.float32(21.0),
 np.float32(22.0),
 np.float32(23.0),
 np.float32(24.0),
 np.float32(25.0),
 np.float32(26.0),
 np.float32(27.0),
 np.float32(28.0),
 np.float32(29.0),
 np.float32(30.0),
 np.float32(31.0),
 np.float32(32.0),
 np.float32(33.0),
 np.float32(34.0),
 np.float32(35.0),
 np.float32(36.0),
 np.float32(37.0),
 np.float32(38.0),
 np.float32(39.0),
 np.float32(40.0),
 np.float32(41.0),
 np.float32(42.0),
 np.float32(43.0),
 np.float32(44.0),
 np.float32(45.0),
 np.float32(46.0),
 np.float32(47.0),
 np.float32(48.0),
 np.float32(49.0),
 np.float32(50.0),
 np.float32(51.0),
 np.float32(52.0),
 np.float32(53.0),
 np.float32(54.0),
 np.float32(55.0),
 np.float32(56.0),
 np.float32(57.0),
 np.float32(58.0),
 np.float32(59.0),
 np.float32(60.0),
 np.float32(61.0),
 np.float32(62.0),
 np.float32(63.0),
 np.float32(64.0),
 np.float32(65.0),
 np.float32(66.0),
 np.float32(67.0),
 np.float32(68.0),
 np.float32(69.0),
 np.float32(70.0),
 np.float32(71.0),
 np.float32(72.0),
 np.float32(73.0),
 np.float32(74.0),
 np.float32(75.0),
 np.float32(76.0),
 np.float32(77.0),
 np.float32(78.0),
 np.float32(79.0),
 np.float32(80.0),
 np.float32(81.0),
 np.float32(82.0),
 np.float32(83.0),
 np.float32(84.0),
 np.float32(85.0),
 np.float32(86.0),
 np.float32(87.0),
 np.float32(88.0),
 np.float32(89.0),
 np.float32(90.0),
 np.float32(91.0),
 np.float32(92.0),
 np.float32(93.0),
 np.float32(94.0),
 np.float32(95.0),
 np.float32(96.0),
 np.float32(97.0),
 np.float32(98.0),
 np.float32(99.0),
 np.float32(100.0),
 np.float32(101.0),
 np.float32(102.0),
 np.float32(103.0),
 np.float32(104.0),
 np.float32(105.0),
 np.float32(106.0),
 np.float32(107.0),
 np.float32(108.0),
 np.float32(109.0),
 np.float32(110.0),
 np.float32(111.0),
 np.float32(112.0),
 np.float32(113.0),
 np.float32(114.0),
 np.float32(115.0),
 np.float32(116.0),
 np.float32(117.0),
 np.float32(118.0),
 np.float32(119.0),
 np.float32(120.0),
 np.float32(121.0),
 np.float32(122.0),
 np.float32(123.0),
 np.float32(124.0),
 np.float32(125.0),
 np.float32(126.0),
 np.float32(127.0),
 np.float32(128.0),
 np.float32(129.0),
 np.float32(130.0),
 np.float32(131.0),
 np.float32(132.0),
 np.float32(133.0),
 np.float32(134.0),
 np.float32(135.0),
 np.float32(136.0),
 np.float32(137.0),
 np.float32(138.0),
 np.float32(139.0),
 np.float32(140.0),
 np.float32(141.0),
 np.float32(142.0),
 np.float32(143.0),
 np.float32(144.0),
 np.float32(145.0),
 np.float32(146.0),
 np.float32(147.0),
 np.float32(148.0),
 np.float32(149.0),
 np.float32(150.0),
 np.float32(151.0),
 np.float32(152.0),
 np.float32(153.0),
 np.float32(154.0),
 np.float32(155.0),
 np.float32(156.0),
 np.float32(157.0),
 np.float32(158.0),
 np.float32(159.0),
 np.float32(160.0),
 np.float32(161.0),
 np.float32(162.0),
 np.float32(163.0),
 np.float32(164.0),
 np.float32(165.0),
 np.float32(166.0),
 np.float32(167.0),
 np.float32(168.0),
 np.float32(169.0),
 np.float32(170.0),
 np.float32(171.0),
 np.float32(172.0),
 np.float32(173.0),
 np.float32(174.0),
 np.float32(175.0),
 np.float32(176.0),
 np.float32(177.0),
 np.float32(178.0),
 np.float32(179.0),
 np.float32(180.0),
 np.float32(181.0),
 np.float32(182.0),
 np.float32(183.0),
 np.float32(184.0),
 np.float32(185.0),
 np.float32(186.0),
 np.float32(187.0),
 np.float32(188.0),
 np.float32(189.0),
 np.float32(190.0),
 np.float32(191.0),
 np.float32(192.0),
 np.float32(193.0),
 np.float32(194.0),
 np.float32(195.0),
 np.float32(196.0),
 np.float32(197.0),
 np.float32(198.0),
 np.float32(199.0),
 np.float32(200.0),
 np.float32(201.0),
 np.float32(202.0),
 np.float32(203.0),
 np.float32(204.0),
 np.float32(205.0),
 np.float32(206.0),
 np.float32(207.0),
 np.float32(208.0),
 np.float32(209.0),
 np.float32(210.0),
 np.float32(211.0),
 np.float32(212.0),
 np.float32(213.0),
 np.float32(214.0),
 np.float32(215.0),
 np.float32(216.0),
 np.float32(217.0),
 np.float32(218.0),
 np.float32(219.0),
 np.float32(220.0),
 np.float32(221.0),
 np.float32(222.0),
 np.float32(223.0),
 np.float32(224.0),
 np.float32(225.0),
 np.float32(226.0),
 np.float32(227.0),
 np.float32(228.0),
 np.float32(229.0),
 np.float32(230.0),
 np.float32(231.0),
 np.float32(232.0),
 np.float32(233.0),
 np.float32(234.0),
 np.float32(235.0),
 np.float32(236.0),
 np.float32(237.0),
 np.float32(238.0),
 np.float32(239.0),
 np.float32(240.0),
 np.float32(241.0),
 np.float32(242.0),
 np.float32(243.0),
 np.float32(244.0),
 np.float32(245.0),
 np.float32(246.0),
 np.float32(247.0),
 np.float32(248.0),
 np.float32(249.0),
 np.float32(250.0),
 np.float32(251.0),
 np.float32(252.0),
 np.float32(253.0),
 np.float32(254.0),
 np.float32(255.0),
 np.float32(256.0),
 np.float32(257.0),
 np.float32(258.0),
 np.float32(259.0),
 np.float32(260.0),
 np.float32(261.0),
 np.float32(262.0),
 np.float32(263.0),
 np.float32(264.0),
 np.float32(265.0),
 np.float32(266.0),
 np.float32(267.0),
 np.float32(268.0),
 np.float32(269.0),
 np.float32(270.0),
 np.float32(271.0),
 np.float32(272.0),
 np.float32(273.0),
 np.float32(274.0),
 np.float32(275.0),
 np.float32(276.0),
 np.float32(277.0),
 np.float32(278.0),
 np.float32(279.0),
 np.float32(280.0),
 np.float32(281.0),
 np.float32(282.0),
 np.float32(283.0),
 np.float32(284.0),
 np.float32(285.0),
 np.float32(286.0),
 np.float32(287.0),
 np.float32(288.0),
 np.float32(289.0),
 np.float32(290.0),
 np.float32(291.0),
 np.float32(292.0),
 np.float32(293.0),
 np.float32(294.0),
 np.float32(295.0),
 np.float32(296.0),
 np.float32(297.0),
 np.float32(298.0),
 np.float32(299.0),
 np.float32(300.0),
 np.float32(301.0),
 np.float32(302.0),
 np.float32(303.0),
 np.float32(304.0),
 np.float32(305.0),
 np.float32(306.0),
 np.float32(307.0),
 np.float32(308.0),
 np.float32(309.0),
 np.float32(310.0),
 np.float32(311.0),
 np.float32(312.0),
 np.float32(313.0),
 np.float32(314.0),
 np.float32(315.0),
 np.float32(316.0),
 np.float32(317.0),
 np.float32(318.0),
 np.float32(319.0),
 np.float32(320.0),
 np.float32(321.0),
 np.float32(322.0),
 np.float32(323.0),
 np.float32(324.0),
 np.float32(325.0),
 np.float32(326.0),
 np.float32(327.0),
 np.float32(328.0),
 np.float32(329.0),
 np.float32(330.0),
 np.float32(331.0),
 np.float32(332.0),
 np.float32(333.0),
 np.float32(334.0),
 np.float32(335.0),
 np.float32(336.0),
 np.float32(337.0),
 np.float32(338.0),
 np.float32(339.0),
 np.float32(340.0),
 np.float32(341.0),
 np.float32(342.0),
 np.float32(343.0),
 np.float32(344.0),
 np.float32(345.0),
 np.float32(346.0),
 np.float32(347.0),
 np.float32(348.0),
 np.float32(349.0),
 np.float32(350.0),
 np.float32(351.0),
 np.float32(352.0),
 np.float32(353.0),
 np.float32(354.0),
 np.float32(355.0),
 np.float32(356.0),
 np.float32(357.0),
 np.float32(358.0),
 np.float32(359.0),
 np.float32(360.0),
 np.float32(361.0),
 np.float32(362.0),
 np.float32(363.0),
 np.float32(364.0),
 np.float32(365.0),
 np.float32(366.0),
 np.float32(367.0),
 np.float32(368.0),
 np.float32(369.0),
 np.float32(370.0),
 np.float32(371.0),
 np.float32(372.0),
 np.float32(373.0),
 np.float32(374.0),
 np.float32(375.0),
 np.float32(376.0),
 np.float32(377.0),
 np.float32(378.0),
 np.float32(379.0),
 np.float32(380.0),
 np.float32(381.0),
 np.float32(382.0),
 np.float32(383.0),
 np.float32(384.0),
 np.float32(385.0),
 np.float32(386.0),
 np.float32(387.0),
 np.float32(388.0),
 np.float32(389.0),
 np.float32(390.0),
 np.float32(391.0),
 np.float32(392.0),
 np.float32(393.0),
 np.float32(394.0),
 np.float32(395.0),
 np.float32(396.0),
 np.float32(397.0),
 np.float32(398.0),
 np.float32(399.0),
 np.float32(400.0),
 np.float32(401.0),
 np.float32(402.0),
 np.float32(403.0),
 np.float32(404.0),
 np.float32(405.0),
 np.float32(406.0),
 np.float32(407.0),
 np.float32(408.0),
 np.float32(409.0),
 np.float32(410.0),
 np.float32(411.0),
 np.float32(412.0),
 np.float32(413.0),
 np.float32(414.0),
 np.float32(415.0),
 np.float32(416.0),
 np.float32(417.0),
 np.float32(418.0),
 np.float32(419.0),
 np.float32(420.0),
 np.float32(421.0),
 np.float32(422.0),
 np.float32(423.0),
 np.float32(424.0),
 np.float32(425.0),
 np.float32(426.0),
 np.float32(427.0),
 np.float32(428.0),
 np.float32(429.0),
 np.float32(430.0),
 np.float32(431.0),
 np.float32(432.0),
 np.float32(433.0),
 np.float32(434.0),
 np.float32(435.0),
 np.float32(436.0),
 np.float32(437.0),
 np.float32(438.0),
 np.float32(439.0),
 np.float32(440.0),
 np.float32(441.0),
 np.float32(442.0),
 np.float32(443.0),
 np.float32(444.0),
 np.float32(445.0),
 np.float32(446.0),
 np.float32(447.0),
 np.float32(448.0),
 np.float32(449.0),
 np.float32(450.0),
 np.float32(451.0),
 np.float32(452.0),
 np.float32(453.0),
 np.float32(454.0),
 np.float32(455.0),
 np.float32(456.0),
 np.float32(457.0),
 np.float32(458.0),
 np.float32(459.0),
 np.float32(460.0),
 np.float32(461.0),
 np.float32(462.0),
 np.float32(463.0),
 np.float32(464.0),
 np.float32(465.0),
 np.float32(466.0),
 np.float32(467.0),
 np.float32(468.0),
 np.float32(469.0),
 np.float32(470.0),
 np.float32(471.0),
 np.float32(472.0),
 np.float32(473.0),
 np.float32(474.0),
 np.float32(475.0),
 np.float32(476.0),
 np.float32(477.0),
 np.float32(478.0),
 np.float32(479.0),
 np.float32(480.0),
 np.float32(481.0),
 np.float32(482.0),
 np.float32(483.0),
 np.float32(484.0),
 np.float32(485.0),
 np.float32(486.0),
 np.float32(487.0),
 np.float32(488.0),
 np.float32(489.0),
 np.float32(490.0),
 np.float32(491.0),
 np.float32(492.0),
 np.float32(493.0),
 np.float32(494.0),
 np.float32(495.0),
 np.float32(496.0),
 np.float32(497.0),
 np.float32(498.0),
 np.float32(499.0),
 np.float32(500.0),
 np.float32(501.0),
 np.float32(502.0),
 np.float32(503.0),
 np.float32(504.0),
 np.float32(505.0),
 np.float32(506.0),
 np.float32(507.0),
 np.float32(508.0),
 np.float32(509.0),
 np.float32(510.0),
 np.float32(511.0),
 np.float32(512.0),
 np.float32(513.0),
 np.float32(514.0),
 np.float32(515.0),
 np.float32(516.0),
 np.float32(517.0),
 np.float32(518.0),
 np.float32(519.0),
 np.float32(520.0),
 np.float32(521.0),
 np.float32(522.0),
 np.float32(523.0),
 np.float32(524.0),
 np.float32(525.0),
 np.float32(526.0),
 np.float32(527.0),
 np.float32(528.0),
 np.float32(529.0),
 np.float32(530.0),
 np.float32(531.0),
 np.float32(532.0),
 np.float32(533.0),
 np.float32(534.0),
 np.float32(535.0),
 np.float32(536.0),
 np.float32(537.0),
 np.float32(538.0),
 np.float32(539.0),
 np.float32(540.0),
 np.float32(541.0),
 np.float32(542.0),
 np.float32(543.0),
 np.float32(544.0),
 np.float32(545.0),
 np.float32(546.0),
 np.float32(547.0),
 np.float32(548.0),
 np.float32(549.0),
 np.float32(550.0),
 np.float32(551.0),
 np.float32(552.0),
 np.float32(553.0),
 np.float32(554.0),
 np.float32(555.0),
 np.float32(556.0),
 np.float32(557.0),
 np.float32(558.0),
 np.float32(559.0),
 np.float32(560.0),
 np.float32(561.0),
 np.float32(562.0),
 np.float32(563.0),
 np.float32(564.0),
 np.float32(565.0),
 np.float32(566.0),
 np.float32(567.0),
 np.float32(568.0),
 np.float32(569.0),
 np.float32(570.0),
 np.float32(571.0),
 np.float32(572.0),
 np.float32(573.0),
 np.float32(574.0),
 np.float32(575.0),
 np.float32(576.0),
 np.float32(577.0),
 np.float32(578.0),
 np.float32(579.0),
 np.float32(580.0),
 np.float32(581.0),
 np.float32(582.0),
 np.float32(583.0),
 np.float32(584.0),
 np.float32(585.0),
 np.float32(586.0),
 np.float32(587.0),
 np.float32(588.0),
 np.float32(589.0),
 np.float32(590.0),
 np.float32(591.0),
 np.float32(592.0),
 np.float32(593.0),
 np.float32(594.0),
 np.float32(595.0),
 np.float32(596.0),
 np.float32(597.0),
 np.float32(598.0),
 np.float32(599.0),
 np.float32(600.0),
 np.float32(601.0),
 np.float32(602.0),
 np.float32(603.0),
 np.float32(604.0),
 np.float32(605.0),
 np.float32(606.0),
 np.float32(607.0),
 np.float32(608.0),
 np.float32(609.0),
 np.float32(610.0),
 np.float32(611.0),
 np.float32(612.0),
 np.float32(613.0),
 np.float32(614.0),
 np.float32(615.0),
 np.float32(616.0),
 np.float32(617.0),
 np.float32(618.0),
 np.float32(619.0),
 np.float32(620.0),
 np.float32(621.0),
 np.float32(622.0),
 np.float32(623.0),
 np.float32(624.0),
 np.float32(625.0),
 np.float32(626.0),
 np.float32(627.0),
 np.float32(628.0),
 np.float32(629.0),
 np.float32(630.0),
 np.float32(631.0),
 np.float32(632.0),
 np.float32(633.0),
 np.float32(634.0),
 np.float32(635.0),
 np.float32(636.0),
 np.float32(637.0),
 np.float32(638.0),
 np.float32(639.0),
 np.float32(640.0),
 np.float32(641.0),
 np.float32(642.0),
 np.float32(643.0),
 np.float32(644.0),
 np.float32(645.0),
 np.float32(646.0),
 np.float32(647.0),
 np.float32(648.0),
 np.float32(649.0),
 np.float32(650.0),
 np.float32(651.0),
 np.float32(652.0),
 np.float32(653.0),
 np.float32(654.0),
 np.float32(655.0),
 np.float32(656.0),
 np.float32(657.0),
 np.float32(658.0),
 np.float32(659.0),
 np.float32(660.0),
 np.float32(661.0),
 np.float32(662.0),
 np.float32(663.0),
 np.float32(664.0),
 np.float32(665.0),
 np.float32(666.0),
 np.float32(667.0),
 np.float32(668.0),
 np.float32(669.0),
 np.float32(670.0),
 np.float32(671.0),
 np.float32(672.0),
 np.float32(673.0),
 np.float32(674.0),
 np.float32(675.0),
 np.float32(676.0),
 np.float32(677.0),
 np.float32(678.0),
 np.float32(679.0),
 np.float32(680.0),
 np.float32(681.0),
 np.float32(682.0),
 np.float32(683.0),
 np.float32(684.0),
 np.float32(685.0),
 np.float32(686.0),
 np.float32(687.0),
 np.float32(688.0),
 np.float32(689.0),
 np.float32(690.0),
 np.float32(691.0),
 np.float32(692.0),
 np.float32(693.0),
 np.float32(694.0),
 np.float32(695.0),
 np.float32(696.0),
 np.float32(697.0),
 np.float32(698.0),
 np.float32(699.0),
 np.float32(700.0),
 np.float32(701.0),
 np.float32(702.0),
 np.float32(703.0),
 np.float32(704.0),
 np.float32(705.0),
 np.float32(706.0),
 np.float32(707.0),
 np.float32(708.0),
 np.float32(709.0),
 np.float32(710.0),
 np.float32(711.0),
 np.float32(712.0),
 np.float32(713.0),
 np.float32(714.0),
 np.float32(715.0),
 np.float32(716.0),
 np.float32(717.0),
 np.float32(718.0),
 np.float32(719.0),
 np.float32(720.0),
 np.float32(721.0),
 np.float32(722.0),
 np.float32(723.0),
 np.float32(724.0),
 np.float32(725.0),
 np.float32(726.0),
 np.float32(727.0),
 np.float32(728.0),
 np.float32(729.0),
 np.float32(730.0),
 np.float32(731.0),
 np.float32(732.0),
 np.float32(733.0),
 np.float32(734.0),
 np.float32(735.0),
 np.float32(736.0),
 np.float32(737.0),
 np.float32(738.0),
 np.float32(739.0),
 np.float32(740.0),
 np.float32(741.0),
 np.float32(742.0),
 np.float32(743.0),
 np.float32(744.0),
 np.float32(745.0),
 np.float32(746.0),
 np.float32(747.0),
 np.float32(748.0),
 np.float32(749.0),
 np.float32(750.0),
 np.float32(751.0),
 np.float32(752.0),
 np.float32(753.0),
 np.float32(754.0),
 np.float32(755.0),
 np.float32(756.0),
 np.float32(757.0),
 np.float32(758.0),
 np.float32(759.0),
 np.float32(760.0),
 np.float32(761.0),
 np.float32(762.0),
 np.float32(763.0),
 np.float32(764.0),
 np.float32(765.0),
 np.float32(766.0),
 np.float32(767.0),
 np.float32(768.0),
 np.float32(769.0),
 np.float32(770.0),
 np.float32(771.0),
 np.float32(772.0),
 np.float32(773.0),
 np.float32(774.0),
 np.float32(775.0),
 np.float32(776.0),
 np.float32(777.0),
 np.float32(778.0),
 np.float32(779.0),
 np.float32(780.0),
 np.float32(781.0),
 np.float32(782.0),
 np.float32(783.0),
 np.float32(784.0),
 np.float32(785.0),
 np.float32(786.0),
 np.float32(787.0),
 np.float32(788.0),
 np.float32(789.0),
 np.float32(790.0),
 np.float32(791.0),
 np.float32(792.0),
 np.float32(793.0),
 np.float32(794.0),
 np.float32(795.0),
 np.float32(796.0),
 np.float32(797.0),
 np.float32(798.0),
 np.float32(799.0),
 np.float32(800.0),
 np.float32(801.0),
 np.float32(802.0),
 np.float32(803.0),
 np.float32(804.0),
 np.float32(805.0),
 np.float32(806.0),
 np.float32(807.0),
 np.float32(808.0),
 np.float32(809.0),
 np.float32(810.0),
 np.float32(811.0),
 np.float32(812.0),
 np.float32(813.0),
 np.float32(814.0),
 np.float32(815.0),
 np.float32(816.0),
 np.float32(817.0),
 np.float32(818.0),
 np.float32(819.0),
 np.float32(820.0),
 np.float32(821.0),
 np.float32(822.0),
 np.float32(823.0),
 np.float32(824.0),
 np.float32(825.0),
 np.float32(826.0),
 np.float32(827.0),
 np.float32(828.0),
 np.float32(829.0),
 np.float32(830.0),
 np.float32(831.0),
 np.float32(832.0),
 np.float32(833.0),
 np.float32(834.0),
 np.float32(835.0),
 np.float32(836.0),
 np.float32(837.0),
 np.float32(838.0),
 np.float32(839.0),
 np.float32(840.0),
 np.float32(841.0),
 np.float32(842.0),
 np.float32(843.0),
 np.float32(844.0),
 np.float32(845.0),
 np.float32(846.0),
 np.float32(847.0),
 np.float32(848.0),
 np.float32(849.0),
 np.float32(850.0),
 np.float32(851.0),
 np.float32(852.0),
 np.float32(853.0),
 np.float32(854.0),
 np.float32(855.0),
 np.float32(856.0),
 np.float32(857.0),
 np.float32(858.0),
 np.float32(859.0),
 np.float32(860.0),
 np.float32(861.0),
 np.float32(862.0),
 np.float32(863.0),
 np.float32(864.0),
 np.float32(865.0),
 np.float32(866.0),
 np.float32(867.0),
 np.float32(868.0),
 np.float32(869.0),
 np.float32(870.0),
 np.float32(871.0),
 np.float32(872.0),
 np.float32(873.0),
 np.float32(874.0),
 np.float32(875.0),
 np.float32(876.0),
 np.float32(877.0),
 np.float32(878.0),
 np.float32(879.0),
 np.float32(880.0),
 np.float32(881.0),
 np.float32(882.0),
 np.float32(883.0),
 np.float32(884.0),
 np.float32(885.0),
 np.float32(886.0),
 np.float32(887.0),
 np.float32(888.0),
 np.float32(889.0),
 np.float32(890.0),
 np.float32(891.0),
 np.float32(892.0),
 np.float32(893.0),
 np.float32(894.0),
 np.float32(895.0),
 np.float32(896.0),
 np.float32(897.0),
 np.float32(898.0),
 np.float32(899.0),
 np.float32(900.0),
 np.float32(901.0),
 np.float32(902.0),
 np.float32(903.0),
 np.float32(904.0),
 np.float32(905.0),
 np.float32(906.0),
 np.float32(907.0),
 np.float32(908.0),
 np.float32(909.0),
 np.float32(910.0),
 np.float32(911.0),
 np.float32(912.0),
 np.float32(913.0),
 np.float32(914.0),
 np.float32(915.0),
 np.float32(916.0),
 np.float32(917.0),
 np.float32(918.0),
 np.float32(919.0),
 np.float32(920.0),
 np.float32(921.0),
 np.float32(922.0),
 np.float32(923.0),
 np.float32(924.0),
 np.float32(925.0),
 np.float32(926.0),
 np.float32(927.0),
 np.float32(928.0),
 np.float32(929.0),
 np.float32(930.0),
 np.float32(931.0),
 np.float32(932.0),
 np.float32(933.0),
 np.float32(934.0),
 np.float32(935.0),
 np.float32(936.0),
 np.float32(937.0),
 np.float32(938.0),
 np.float32(939.0),
 np.float32(940.0),
 np.float32(941.0),
 np.float32(942.0),
 np.float32(943.0),
 np.float32(944.0),
 np.float32(945.0),
 np.float32(946.0),
 np.float32(947.0),
 np.float32(948.0),
 np.float32(949.0),
 np.float32(950.0),
 np.float32(951.0),
 np.float32(952.0),
 np.float32(953.0),
 np.float32(954.0),
 np.float32(955.0),
 np.float32(956.0),
 np.float32(957.0),
 np.float32(958.0),
 np.float32(959.0),
 np.float32(960.0),
 np.float32(961.0),
 np.float32(962.0),
 np.float32(963.0),
 np.float32(964.0),
 np.float32(965.0),
 np.float32(966.0),
 np.float32(967.0),
 np.float32(968.0),
 np.float32(969.0),
 np.float32(970.0),
 np.float32(971.0),
 np.float32(972.0),
 np.float32(973.0),
 np.float32(974.0),
 np.float32(975.0),
 np.float32(976.0),
 np.float32(977.0),
 np.float32(978.0),
 np.float32(979.0),
 np.float32(980.0),
 np.float32(981.0),
 np.float32(982.0),
 np.float32(983.0),
 np.float32(984.0),
 np.float32(985.0),
 np.float32(986.0),
 np.float32(987.0),
 np.float32(988.0),
 np.float32(989.0),
 np.float32(990.0),
 np.float32(991.0),
 np.float32(992.0),
 np.float32(993.0),
 np.float32(994.0),
 np.float32(995.0),
 np.float32(996.0),
 np.float32(997.0),
 np.float32(998.0),
 np.float32(999.0),
 np.float32(1000.0),
 ...]

[17]:

h.plot(totim=900, contour=True, grid=True, colorbar=True, figsize=(10, 10))

[17]:

[<Axes: title={'center': 'data Layer 1'}>,
 <Axes: title={'center': 'data Layer 2'}>,
 <Axes: title={'center': 'data Layer 3'}>]

../_images/Notebooks_feat_working_stack_examples_26_1.png

../_images/Notebooks_feat_working_stack_examples_26_2.png

../_images/Notebooks_feat_working_stack_examples_26_3.png

We can write the heads to a shapefile.

[18]:

h.to_shapefile(ml.model_ws / "freyburg_head.shp", verbose=False)

Finally, let’s make an animation of the simulated head over the time domain.

[19]:

f = plt.figure(figsize=(10, 10))
ax = plt.subplot(1, 1, 1, aspect="equal")
for t in h.times[0:-1:10]:
    ax.cla()

    ax.set_title(f"totim: {t:4.0f} days")
    mm = flopy.plot.PlotMapView(model=ml, ax=ax)
    mm.plot_array(h.get_data(totim=t), vmin=0, vmax=20)
    mm.plot_grid(lw=0.5, color="black")

    display(f)
    clear_output(wait=True)
    plt.pause(0.1)

../_images/Notebooks_feat_working_stack_examples_30_0.png

[20]:

try:
    # ignore PermissionError on Windows
    temp_dir.cleanup()
except:
    pass