notebooks/learning/gwf_fhb.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pathlib import Path\n",
    "\n",
    "import flopy as fp\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 1. Load and explore model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "loading simulation...\n",
      "  loading simulation name file...\n",
      "  loading tdis package...\n",
      "  loading model gwf6...\n",
      "    loading package dis...\n",
      "    loading package npf...\n",
      "    loading package ic...\n",
      "    loading package sto...\n",
      "    loading package chd...\n",
      "    loading package wel...\n",
      "    loading package oc...\n",
      "    loading package obs...\n",
      "  loading solution package ex-gwf-fhb...\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "sim_name = ex-gwf-fhb\n",
       "sim_path = C:\\Users\\buckl\\Documents\\Software\\mf6.4.1\\examples\\ex-gwf-fhb\n",
       "exe_name = C:\\Users\\buckl\\Documents\\Software\\mf6.4.1\\bin\\mf6.exe\n",
       "\n",
       "###################\n",
       "Package mfsim.nam\n",
       "###################\n",
       "\n",
       "package_name = mfsim.nam\n",
       "filename = mfsim.nam\n",
       "package_type = nam\n",
       "model_or_simulation_package = simulation\n",
       "simulation_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package ex-gwf-fhb.tdis\n",
       "###################\n",
       "\n",
       "package_name = ex-gwf-fhb.tdis\n",
       "filename = ex-gwf-fhb.tdis\n",
       "package_type = tdis\n",
       "model_or_simulation_package = simulation\n",
       "simulation_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package ex-gwf-fhb\n",
       "###################\n",
       "\n",
       "package_name = ex-gwf-fhb\n",
       "filename = ex-gwf-fhb.ims\n",
       "package_type = ims\n",
       "model_or_simulation_package = simulation\n",
       "simulation_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "@@@@@@@@@@@@@@@@@@@@\n",
       "Model ex-gwf-fhb\n",
       "@@@@@@@@@@@@@@@@@@@@\n",
       "\n",
       "name = ex-gwf-fhb\n",
       "model_type = gwf6\n",
       "version = mf6\n",
       "model_relative_path = \n",
       "\n",
       "###################\n",
       "Package dis\n",
       "###################\n",
       "\n",
       "package_name = dis\n",
       "filename = ex-gwf-fhb.dis\n",
       "package_type = dis\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package npf\n",
       "###################\n",
       "\n",
       "package_name = npf\n",
       "filename = ex-gwf-fhb.npf\n",
       "package_type = npf\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package ic\n",
       "###################\n",
       "\n",
       "package_name = ic\n",
       "filename = ex-gwf-fhb.ic\n",
       "package_type = ic\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package sto\n",
       "###################\n",
       "\n",
       "package_name = sto\n",
       "filename = ex-gwf-fhb.sto\n",
       "package_type = sto\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package chd_ts\n",
       "###################\n",
       "\n",
       "package_name = chd_ts\n",
       "filename = ex-gwf-fhb.chd.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = chd\n",
       "\n",
       "\n",
       "###################\n",
       "Package chd\n",
       "###################\n",
       "\n",
       "package_name = chd\n",
       "filename = ex-gwf-fhb.chd\n",
       "package_type = chd\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package wel_ts\n",
       "###################\n",
       "\n",
       "package_name = wel_ts\n",
       "filename = ex-gwf-fhb.wel.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = wel\n",
       "\n",
       "\n",
       "###################\n",
       "Package wel\n",
       "###################\n",
       "\n",
       "package_name = wel\n",
       "filename = ex-gwf-fhb.wel\n",
       "package_type = wel\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package oc\n",
       "###################\n",
       "\n",
       "package_name = oc\n",
       "filename = ex-gwf-fhb.oc\n",
       "package_type = oc\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package obs_0\n",
       "###################\n",
       "\n",
       "package_name = obs_0\n",
       "filename = ex-gwf-fhb.obs\n",
       "package_type = obs\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "workspace = Path(\"C:\\\\Users\\\\buckl\\\\Documents\\\\Software\\\\mf6.4.1\\\\examples\\\\ex-gwf-fhb\")\n",
    "sim_name = \"ex-gwf-fhb\"\n",
    "sim = fp.mf6.MFSimulation.load(\n",
    "    sim_name=sim_name,\n",
    "    version=\"mf6\",\n",
    "    exe_name=\"C:\\\\Users\\\\buckl\\\\Documents\\\\Software\\\\mf6.4.1\\\\bin\\\\mf6.exe\",\n",
    "    sim_ws=workspace\n",
    ")\n",
    "sim"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['ex-gwf-fhb']"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim.model_names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "name = ex-gwf-fhb\n",
       "model_type = gwf6\n",
       "version = mf6\n",
       "model_relative_path = \n",
       "\n",
       "###################\n",
       "Package dis\n",
       "###################\n",
       "\n",
       "package_name = dis\n",
       "filename = ex-gwf-fhb.dis\n",
       "package_type = dis\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package npf\n",
       "###################\n",
       "\n",
       "package_name = npf\n",
       "filename = ex-gwf-fhb.npf\n",
       "package_type = npf\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package ic\n",
       "###################\n",
       "\n",
       "package_name = ic\n",
       "filename = ex-gwf-fhb.ic\n",
       "package_type = ic\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package sto\n",
       "###################\n",
       "\n",
       "package_name = sto\n",
       "filename = ex-gwf-fhb.sto\n",
       "package_type = sto\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package chd_ts\n",
       "###################\n",
       "\n",
       "package_name = chd_ts\n",
       "filename = ex-gwf-fhb.chd.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = chd\n",
       "\n",
       "\n",
       "###################\n",
       "Package chd\n",
       "###################\n",
       "\n",
       "package_name = chd\n",
       "filename = ex-gwf-fhb.chd\n",
       "package_type = chd\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package wel_ts\n",
       "###################\n",
       "\n",
       "package_name = wel_ts\n",
       "filename = ex-gwf-fhb.wel.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = wel\n",
       "\n",
       "\n",
       "###################\n",
       "Package wel\n",
       "###################\n",
       "\n",
       "package_name = wel\n",
       "filename = ex-gwf-fhb.wel\n",
       "package_type = wel\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package oc\n",
       "###################\n",
       "\n",
       "package_name = oc\n",
       "filename = ex-gwf-fhb.oc\n",
       "package_type = oc\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "\n",
       "###################\n",
       "Package obs_0\n",
       "###################\n",
       "\n",
       "package_name = obs_0\n",
       "filename = ex-gwf-fhb.obs\n",
       "package_type = obs\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model = sim.get_model('ex-gwf-fhb')\n",
    "model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['dis', 'npf', 'ic', 'sto', 'chd_ts', 'chd', 'wel_ts', 'wel', 'oc', 'obs_0']"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.package_names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdoAAAC1CAYAAAD4MT6UAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAlrUlEQVR4nO3deVRTd9oH8G+AJEAhIFCIVqBYHfcdl7Su1REtWludUTtWaXXs6OCCeNzeWm2d6cCxi6J19J22YnveutGjoqJYXtywRawoFZdaq7b41gY7RUCpEJbn/cPhTiKLhnJJdL6fc3LM/f2e3Dy/X0Ie7829NxoRERAREZEqXBydABER0cOMhZaIiEhFLLREREQqYqElIiJSEQstERGRilhoiYiIVMRCS0REpCIWWiIiIhWx0BIREamIhZaoiQwaNAgajQYajQbfffedo9P5VTZt2oSuXbvC09MTGo0Gvr6+jk7pvjz++OPKa+DIddB/FhZacrjXX39d+eDSaDQYNmxYjZjs7GybGI1Gg9LSUgdka5vv5MmTa43Zv3+/EtOmTZsmzlBdmZmZePHFF3H69Gncvn37vh5z92v3xRdf1IgZNWqUTcyiRYsaO3Uih2ChJaeTnp6O77//3qbt/fffd1A2Nb3wwgvK/V27dqGsrKxGTFJSknJ/woQJAIA1a9YgIyMDGRkZaN68ufqJqiQlJQXVl0j/05/+hCNHjiA1NdWudXzwwQc2yz/88AP27dvXaDkSORMWWnI6VVVV+PDDD5XlkpISbNq0yYEZ2Wrbti26desGACgqKsL+/ftt+isqKrBjxw5lubrQdu7cGf369UO/fv2g1+ubLN/Gdu3aNeX++PHj0b9/f/Tt29eudWzbtg03b95Uljds2IDKyspGy5HImbDQklPx9vYGACQmJqKqqgoAsHXrVty8eVPpq828efPw5JNPonnz5tDr9fDy8kKPHj3w9ttvo6Kiwia2etfk448/josXL2LkyJHw8vJCQEAAoqOjUVJScs88rbdqt23bZtOXnp6OgoICAHeKa8eOHQHU/x1tcnIyhg4dimbNmkGv16Nt27Z44403bHbNrlmzRnm89RZhv379lPFU27dvnxK7cOHCe47HbDZj9uzZeOKJJ6DX6+Hr64tBgwbZbJl/99130Gg0SExMVNqefvppaDQaDBo06J7PUc3b2xslJSXYvHkzANv/WNX3Gt9PjtV++eUXzJ49G48++ii8vLzw7LPP1vu9uIggMTERTz31FAwGAzw8PNC1a1ckJCQo70OiBhMiB1u2bJkAEADy0ksviVarFQCSkpIiIiJ9+vQRAPLKK68ocQDk9u3byjr0er1Nn/Xt5Zdftnm+6nZfX18JCgqqET98+PB75vz999+LRqMRAOLt7S2lpaVK35QpU5R1vfnmm0r7wIEDlfYrV64o7a+99lqduffv31/KyspEROTUqVNK+9SpU0VExGKxiLu7u9J+7dq1Guvcs2dPvWO5fPmyGI3GOnNYuHChiIhcuXKlzpiBAwfW+xzWsdOmTRMA0qtXLxER2bdvnwAQV1dXmTp1ao3ntSfHapGRkTViWrZsKX5+fsqytcmTJ9e57vHjx9vEhoaG1roOorrwnUIOZ11oFy5cKM8//7wAkOeff15Onz6t9B0/frzOQrt8+XLZvHmzpKamyqFDh2T79u1KgdZoNHL16lUl1nodPXv2lJ07d8qaNWvE09NTad+1a9c9837yySeV+B07doiISHl5uc2H+aVLl5T42gqt9ZiaN28uH374oaSmptoUivj4eBERqaysFB8fHwEgHTp0EBGRrKwsm/F8+umnIiIyZMgQASAuLi5y48aNescxfPhw5fGDBg2SXbt2ybvvvmtTwI8dOyalpaWSkZEhI0aMUNpXr14tGRkZcvr06XqfwzpH65xPnz4tY8aMEQDyzDPP1Hgv2JujiEhqaqrS5uHhIatWrZKdO3dKeHi4TR7VkpKSlLa2bdvK5s2bZffu3dK3b1+lfcuWLUo8Cy3Zi+8Ucri7P1xTUlIEgGi1Whk3bpwAkC5duoiI7Qe2daE9evSojB49WoxGo7i5udXYKklOTlZirdsvXryotL/66qtK+5QpU+6Z95o1a5T4F154QUT+vXUGQHr37m0TX1uhnTNnjtL2X//1X5KRkSEZGRmye/dupb1Tp07KOp555hnlPw83btyQVatWCQDp2LGjAJB58+ZJZWWleHt728xbXX7++Wdly1yv18s///lPpW/evHlKDnPmzFHao6KilPaDBw/ec55Ear5uvXr1EgAybtw4ZQ/G9u3bay209uY4Y8YMpW3+/PlK7DfffFNroR09enSN/zhkZGTI+++/r7SPHDlSiWehJXvxO1pyOsOHD0dwcDDKy8uV7z+nTZtWZ/zx48cxePBgJCcnw2w21/hOFgAKCwtrtPn5+aF169bKcu/evZX7ly9fvmeev//97+Hq6goA2L17N0pLS2s92rg+33zzjXL/b3/7G/r374/+/ftj1KhRSvvXX3+t3B8wYAAAQESQlZWFzMxMAEBMTAwA4IsvvkBubq5yoFH//v3rff6LFy8qRxA/8cQT8Pf3V/qs58M6z8ZQ/Xpu27YN5eXlMBqNNmP+NTlav3a9evVS7rdp0wbNmjWrsX7rsc2ePVt5Dazfc+fPn7drfETWWGjJ6bi4uODll19Wlt3d3fHiiy/WGb9+/XqUl5cDAEaOHIm9e/ciIyPD5hzX+zmgxd4LEAQFBWHw4MEAgFu3biE5ORk7d+5UxjBu3Di71leXiooK5RSi6kIL3DmfNTMzE66urpgwYQICAwNx8uRJHDp0SIm5V6Gtj5oXZJgwYQIeeeQRZTkqKgpubm52r8feHBs6pvs5QI6oLiy05JSmTJkCF5c7b8+xY8fWe+WhH374QbkfFxeHESNGoF+/fsjPz6/3OQoKCvDtt98qy1lZWcr9Vq1a3Vee1kcfz58/XznauH///njsscfu+fjf/OY3yv3ExETIna9zbG4lJSXK6UDh4eHw8PAAAGzfvh15eXno3LkzvLy8YDKZUFZWhnXr1inrtC7MtWndurVSfC5duoSff/5Z6bOeD+s8G4O3tzfGjx+vLE+dOrXRcrR+7U6cOKHc//bbb5XXx5r12A4ePFjra3Dp0iV7hkdkw/7/QhI1gdDQUKxduxZmsxm/+93v7hlbLS4uDlFRUdi3b1+N81tr84c//AFLlizB//3f/2HVqlVK++jRo+8rzzFjxmDGjBmwWCy4evWq0n4/u42rnz8hIQEAMHfuXBQUFKBLly4oLCzEpUuX8NlnnyE0NBQbNmwAAGi1WvTt2xcHDx5Ebm4uAMBkMin/Jicn48KFCwDu7Ga914Ux/P39ERERgdTUVJSVlWHcuHGYO3cuLl26hL///e9KnPV/KBrL/PnzERwcjICAgHqvnmVvjs8++6zyn4333nsPLVu2RGhoKN58881a1z9x4kQkJycDACZNmoRXX30Vbdq0wU8//YSLFy8iJSUFI0aMwLJlyxpr6PSfxiHfDBNZqetI09oANQ+GysrKUg6Wqb5pNBoxmUzKcmJiYo11+Pn5ScuWLWscOPXb3/5Wqqqq7jv/Z5991ubxbm5u8tNPP9WIa8jpPQAkKirKZj1Lly616f/4449FROTw4cM27S+99NJ95X/p0iW7Tp1pjIOh6lLXe8HeHK2PjK6+Pfroo8pR23d/9NV3eg8AWbZsmRLLg6HIXtx1TA+83r17Y8eOHejcuTPc3d3RsWNHJCUl1XrNZGve3t7IyMjAqFGj8Mgjj8DPzw/Tp0/H9u3b7fou7+6tvSFDhiAgIOC+H798+XLs2bMHw4cPh7+/P7RaLR577DH069cP8fHxeOONN2zi794dXL1FGx4ebvM95/1+P9uqVSucPHkSM2fORFhYGLRaLQwGAwYMGICtW7ciPj7+vseiFntzTEpKQnR0NPz9/eHp6YmIiAgcOXKkzq8gPvroI3z88ccYOHAgfHx8oNPpEBISgiFDhmD16tX485//3ASjpIeVRuRfh/MR/YeoLqKhoaEP/K/oEJHz4xYtERGRilhoiYiIVMRCS0REpCK7Cu26devQpUsXGAwGGAwGmEwmm9+QLC0tVQ5A8PLywtixY2ucy5iXl4fIyEh4enoiMDAQ8+fPr3Eln0OHDqFHjx7Q6/Vo3bo1Nm7c2PAREt1F/nVuJL+fJaKmYFehbdmyJeLj45GdnY0TJ07g6aefxujRo3H27FkAd84D3L17N5KSknD48GFcu3YNY8aMUR5fWVmJyMhIWCwWfPHFF/joo4+wceNGLF26VIm5cuUKIiMjMXjwYOTk5CAmJgZ//OMf7+ucSCIiImfzq4869vPzw1tvvYXf/e53ePTRR7Fp0yblAgNff/012rdvj8zMTPTt2xf79u3DyJEjce3aNQQFBQG4c/m8hQsX4qeffoJOp8PChQuRkpKCM2fOKM8xYcIEFBYWIjU19dekSkRE1OQafGWoyspKJCUloaSkBCaTCdnZ2SgvL8fQoUOVmHbt2iEkJEQptJmZmejcubNSZAEgIiICM2bMwNmzZ9G9e3dkZmbarKM6pvqi6XUpKytTrgcL3Lm2bUFBAfz9/VW9ZisRETkvEcHNmzfRokUL5bKuTc3uQpubmwuTyYTS0lJ4eXlhx44d6NChA3JycqDT6WqcEB4UFASz2QwAMJvNNkW2ur+6r76Y4uJi3L59W7nO693i4uJqnNhPREQEAFevXkXLli0d8tx2F9q2bdsiJycHRUVF+PTTTxEVFYXDhw+rkZtdFi9ejNjYWGW5qKgIISEh+Prrr21+VsvRLBYL1qxZAwCYNWsWdDqdgzO6w1nzAphbQzlrbs6aF8DcGso6t4QE4F8/puVwWi0wfXoZVq5cCW9vb4flYXeh1el0ym949uzZE19++SUSEhIwfvx4WCwWFBYW2mzV5ufnw2g0AgCMRiOOHz9us77qo5KtY+4+Ujk/Px8Gg6HOrVkA0Ov1yi+cWPP397frcnhqs1gscHd3BwAEBAQ4zR+Ls+YFMLeGctbcnDUvgLk1lHVuFRV3bs5AowGqy4Ijv0L81Tusq6qqUFZWhp49e0Kr1SI9PV3pu3DhAvLy8mx+XSQ3NxfXr19XYtLS0mAwGNChQwclxnod1THV6yAiInqQ2LVFu3jxYowYMQIhISG4efMmNm3ahEOHDmH//v3w8fHB1KlTERsbCz8/PxgMBsyaNQsmkwl9+/YFAAwbNgwdOnTApEmTsGLFCpjNZixZsgTR0dHK1uj06dPx3nvvYcGCBZgyZQoOHDiAbdu2ISUlpfFHT0REpDK7Cu3169cxefJk/Pjjj/Dx8UGXLl2wf/9+/Pa3vwUArFy5Ei4uLhg7dizKysoQERFh83uRrq6u2LNnD2bMmAGTyYRHHnkEUVFRWL58uRITFhaGlJQUzJ07FwkJCWjZsiU++OADRERENNKQiYiImo5dhfbDDz+st9/d3R1r167F2rVr64wJDQ3F3r17613PoEGDcOrUKXtSIyIickq81jEREZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIrsKbVxcHHr16gVvb28EBgbiueeew4ULF2xiSktLER0dDX9/f3h5eWHs2LHIz8+3icnLy0NkZCQ8PT0RGBiI+fPno6Kiwibm0KFD6NGjB/R6PVq3bo2NGzc2bIREREQOZFehPXz4MKKjo3Hs2DGkpaWhvLwcw4YNQ0lJiRIzd+5c7N69G0lJSTh8+DCuXbuGMWPGKP2VlZWIjIyExWLBF198gY8++ggbN27E0qVLlZgrV64gMjISgwcPRk5ODmJiYvDHP/4R+/fvb4QhExERNR03e4JTU1Ntljdu3IjAwEBkZ2djwIABKCoqwocffohNmzbh6aefBgAkJiaiffv2OHbsGPr27YvPPvsM586dw//+7/8iKCgI3bp1w1/+8hcsXLgQr7/+OnQ6HdavX4+wsDC88847AID27dvj6NGjWLlyJSIiIuwaoMVigcVisesxarLOhXndH+bWMM6am7PmBTC3hrLOR6t1YCJ3cZZcNCIiDX3wt99+izZt2iA3NxedOnXCgQMHMGTIENy4cQO+vr5KXGhoKGJiYjB37lwsXboUu3btQk5OjtJ/5coVtGrVCidPnkT37t0xYMAA9OjRA6tWrVJiEhMTERMTg6KiolpzKSsrQ1lZmbJcXFyM4OBgLFq0CO7u7g0dIhERPcBKS0sRHx+PoqIiGAwGh+TQ4IOhqqqqEBMTg6eeegqdOnUCAJjNZuh0OpsiCwBBQUEwm81KTFBQUI3+6r76YoqLi3H79u1a84mLi4OPj49yCw4ObujQiIiIGo1du46tRUdH48yZMzh69Ghj5tNgixcvRmxsrLJcvUU7KyEBAXcdaOVIFq0W7yxYAACYN28edDqdgzO6w2KxKLvqnSkvgLk1lLPm5qx5AcytoaxzW7ECKC93cEL/otUCs2c7OosGFtqZM2diz549OHLkCFq2bKm0G41GWCwWFBYW2mzV5ufnw2g0KjHHjx+3WV/1UcnWMXcfqZyfnw+DwQAPD49ac9Lr9dDr9TXadeXl0DlRobWm0+mc6o+lmrPmBTC3hnLW3Jw1L4C5NVR5ufMUWmdh165jEcHMmTOxY8cOHDhwAGFhYTb9PXv2hFarRXp6utJ24cIF5OXlwWQyAQBMJhNyc3Nx/fp1JSYtLQ0GgwEdOnRQYqzXUR1TvQ4iIqIHhV1btNHR0di0aROSk5Ph7e2tfKfq4+MDDw8P+Pj4YOrUqYiNjYWfnx8MBgNmzZoFk8mEvn37AgCGDRuGDh06YNKkSVixYgXMZjOWLFmC6OhoZYt0+vTpeO+997BgwQJMmTIFBw4cwLZt25CSktLIwyciIlKXXVu069atQ1FREQYNGoTmzZsrt61btyoxK1euxMiRIzF27FgMGDAARqMR27dvV/pdXV2xZ88euLq6wmQy4cUXX8TkyZOxfPlyJSYsLAwpKSlIS0tD165d8c477+CDDz6w+9QeIiIiR7Nri/Z+zgRyd3fH2rVrsXbt2jpjQkNDsXfv3nrXM2jQIJw6dcqe9IiIiJwOr3VMRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEUstERERCpioSUiIlIRCy0REZGKWGiJiIhUxEJLRESkIhZaIiIiFbHQEhERqYiFloiISEV2F9ojR45g1KhRaNGiBTQaDXbu3GnTLyJYunQpmjdvDg8PDwwdOhQXL160iSkoKMDEiRNhMBjg6+uLqVOn4tatWzYxp0+fRv/+/eHu7o7g4GCsWLHC/tERERE5mN2FtqSkBF27dsXatWtr7V+xYgVWr16N9evXIysrC4888ggiIiJQWlqqxEycOBFnz55FWloa9uzZgyNHjuCVV15R+ouLizFs2DCEhoYiOzsbb731Fl5//XX84x//aMAQiYiIHMfN3geMGDECI0aMqLVPRLBq1SosWbIEo0ePBgB8/PHHCAoKws6dOzFhwgScP38eqamp+PLLLxEeHg4AWLNmDZ555hm8/fbbaNGiBT755BNYLBZs2LABOp0OHTt2RE5ODt59912bgnw/LFotLBqNvcNUjUWr/fd9i8WBmdiyzsWZ8gKYW0M5a27OmhfA3BrKOh+rjziHc5ZcNCIiDX6wRoMdO3bgueeeAwBcvnwZTzzxBE6dOoVu3bopcQMHDkS3bt2QkJCADRs2YN68ebhx44bSX1FRAXd3dyQlJeH555/H5MmTUVxcbLNb+uDBg3j66adRUFCAZs2a1cilrKwMZWVlynJxcTGCg4OxaNEiuLu7N3SIRET0ACstLUV8fDyKiopgMBgckkOjHgxlNpsBAEFBQTbtQUFBSp/ZbEZgYKBNv5ubG/z8/GxialuH9XPcLS4uDj4+PsotODj41w+IiIjoV7J717GzWrx4MWJjY5Xl6i3aWbNmISAgwIGZ2bJYLHjnnXcAAPPmzYNOp3NwRnc4a14Ac2soZ83NWfMCmFtDWee2YgVQXu7ghP5FqwVmz3Z0Fo1caI1GIwAgPz8fzZs3V9rz8/OVXclGoxHXr1+3eVxFRQUKCgqUxxuNRuTn59vEVC9Xx9xNr9dDr9fXaNfpdE71hrTmrLk5a14Ac2soZ83NWfMCmFtDlZc7T6F1Fo266zgsLAxGoxHp6elKW3FxMbKysmAymQAAJpMJhYWFyM7OVmIOHDiAqqoq9OnTR4k5cuQIyq1erbS0NLRt27bW72eJiIicld2F9tatW8jJyUFOTg4A4MqVK8jJyUFeXh40Gg1iYmLw17/+Fbt27UJubi4mT56MFi1aKAdMtW/fHsOHD8e0adNw/PhxfP7555g5cyYmTJiAFi1aAAD+8Ic/QKfTYerUqTh79iy2bt2KhIQEm13DREREDwK7dx2fOHECgwcPVpari19UVBQ2btyIBQsWoKSkBK+88goKCwvRr18/pKam2hz5+8knn2DmzJkYMmQIXFxcMHbsWKxevVrp9/HxwWeffYbo6Gj07NkTAQEBWLp0qd2n9hARETma3YV20KBBqO+MII1Gg+XLl2P58uV1xvj5+WHTpk31Pk+XLl2QkZFhb3pEREROhdc6JiIiUhELLRERkYpYaImIiFTEQktERKQiFloiIiIVsdASERGpiIWWiIhIRSy0REREKmKhJSIiUhELLRERkYpYaImIiFTEQktERKQiFloiIiIVsdASERGpiIWWiIhIRSy0REREKmKhJSIiUhELLRERkYpYaImIiFTEQktERKQiFloiIiIVsdASERGpiIWWiIhIRSy0REREKmKhJSIiUhELLRERkYpYaImIiFTk1IV27dq1ePzxx+Hu7o4+ffrg+PHjjk6JiIjILk5baLdu3YrY2FgsW7YMJ0+eRNeuXREREYHr1687OjUiIqL75uboBOry7rvvYtq0aXj55ZcBAOvXr0dKSgo2bNiARYsW3fd6LBYLLBaLWmnazToX5nV/mFvDOGtuzpoXwNwayjofrdaBidzFWXLRiIg4Oom7WSwWeHp64tNPP8Vzzz2ntEdFRaGwsBDJyck1HlNWVoaysjJluaioCCEhIZg7dy70en1TpE1ERE6mrKwMK1euRGFhIXx8fBySg1Nu0f7zn/9EZWUlgoKCbNqDgoLw9ddf1/qYuLg4vPHGGzXaV65cqUqORET04Pj5559ZaH+txYsXIzY2VlkuLCxEaGgo8vLyHDa5D5ri4mIEBwfj6tWrMBgMjk7ngcF5sx/nrGE4b/ar3rvp5+fnsBycstAGBATA1dUV+fn5Nu35+fkwGo21Pkav19e6i9jHx4dvSDsZDAbOWQNw3uzHOWsYzpv9XFwcd+yvUx51rNPp0LNnT6SnpyttVVVVSE9Ph8lkcmBmRERE9nHKLVoAiI2NRVRUFMLDw9G7d2+sWrUKJSUlylHIREREDwKnLbTjx4/HTz/9hKVLl8JsNqNbt25ITU2tcYBUXfR6PZYtW8Yjju3AOWsYzpv9OGcNw3mznzPMmVOe3kNERPSwcMrvaImIiB4WLLREREQqYqElIiJSEQstERGRilhoiYiIVPRQFtr/5N+xjYuLQ69eveDt7Y3AwEA899xzuHDhgk1MaWkpoqOj4e/vDy8vL4wdO7bGVbjy8vIQGRkJT09PBAYGYv78+aioqLCJOXToEHr06AG9Xo/WrVtj48aNag+vScTHx0Oj0SAmJkZp45zV7ocffsCLL74If39/eHh4oHPnzjhx4oTSLyJYunQpmjdvDg8PDwwdOhQXL160WUdBQQEmTpwIg8EAX19fTJ06Fbdu3bKJOX36NPr37w93d3cEBwdjxYoVTTK+xlZZWYnXXnsNYWFh8PDwwBNPPIG//OUvsD75g3MGHDlyBKNGjUKLFi2g0Wiwc+dOm/6mnKOkpCS0a9cO7u7u6Ny5M/bu3Wv/gOQhs2XLFtHpdLJhwwY5e/asTJs2TXx9fSU/P9/RqTWJiIgISUxMlDNnzkhOTo4888wzEhISIrdu3VJipk+fLsHBwZKeni4nTpyQvn37ypNPPqn0V1RUSKdOnWTo0KFy6tQp2bt3rwQEBMjixYuVmMuXL4unp6fExsbKuXPnZM2aNeLq6iqpqalNOt7Gdvz4cXn88celS5cuMmfOHKWdc1ZTQUGBhIaGyksvvSRZWVly+fJl2b9/v3z77bdKTHx8vPj4+MjOnTvlq6++kmeffVbCwsLk9u3bSszw4cOla9eucuzYMcnIyJDWrVvLCy+8oPQXFRVJUFCQTJw4Uc6cOSObN28WDw8P+e///u8mHW9jePPNN8Xf31/27NkjV65ckaSkJPHy8pKEhAQlhnMmsnfvXnn11Vdl+/btAkB27Nhh099Uc/T555+Lq6urrFixQs6dOydLliwRrVYrubm5do3noSu0vXv3lujoaGW5srJSWrRoIXFxcQ7MynGuX78uAOTw4cMiIlJYWCharVaSkpKUmPPnzwsAyczMFJE7b3IXFxcxm81KzLp168RgMEhZWZmIiCxYsEA6duxo81zjx4+XiIgItYekmps3b0qbNm0kLS1NBg4cqBRazlntFi5cKP369auzv6qqSoxGo7z11ltKW2Fhoej1etm8ebOIiJw7d04AyJdffqnE7Nu3TzQajfzwww8iIvL3v/9dmjVrpsxj9XO3bdu2sYekusjISJkyZYpN25gxY2TixIkiwjmrzd2FtinnaNy4cRIZGWmTT58+feRPf/qTXWN4qHYdWywWZGdnY+jQoUqbi4sLhg4diszMTAdm5jhFRUUAoPxyRXZ2NsrLy23mqF27dggJCVHmKDMzE507d7a5CldERASKi4tx9uxZJcZ6HdUxD/I8R0dHIzIyssa4OGe127VrF8LDw/H73/8egYGB6N69O95//32l/8qVKzCbzTZj9vHxQZ8+fWzmzdfXF+Hh4UrM0KFD4eLigqysLCVmwIAB0Ol0SkxERAQuXLiAGzduqD3MRvXkk08iPT0d33zzDQDgq6++wtGjRzFixAgAnLP70ZRz1Fh/sw9Voa3vd2zNZrODsnKcqqoqxMTE4KmnnkKnTp0AAGazGTqdDr6+vjax1nNkNptrncPqvvpiiouLcfv2bTWGo6otW7bg5MmTiIuLq9HHOavd5cuXsW7dOrRp0wb79+/HjBkzMHv2bHz00UcA/j3u+v4ezWYzAgMDbfrd3Nzg5+dn19w+KBYtWoQJEyagXbt20Gq16N69O2JiYjBx4kQAnLP70ZRzVFeMvXPotNc6pl8vOjoaZ86cwdGjRx2dilO7evUq5syZg7S0NLi7uzs6nQdGVVUVwsPD8be//Q0A0L17d5w5cwbr169HVFSUg7NzTtu2bcMnn3yCTZs2oWPHjsjJyUFMTAxatGjBOXuIPVRbtA35HduH1cyZM7Fnzx4cPHgQLVu2VNqNRiMsFgsKCwtt4q3nyGg01jqH1X31xRgMBnh4eDT2cFSVnZ2N69evo0ePHnBzc4ObmxsOHz6M1atXw83NDUFBQZyzWjRv3hwdOnSwaWvfvj3y8vIA/Hvc9f09Go1GXL9+3aa/oqICBQUFds3tg2L+/PnKVm3nzp0xadIkzJ07V9mTwjm7t6aco7pi7J3Dh6rQ8nds7xz2PnPmTOzYsQMHDhxAWFiYTX/Pnj2h1Wpt5ujChQvIy8tT5shkMiE3N9fmjZqWlgaDwaB8sJpMJpt1VMc8iPM8ZMgQ5ObmIicnR7mFh4dj4sSJyn3OWU1PPfVUjVPHvvnmG4SGhgIAwsLCYDQabcZcXFyMrKwsm3krLCxEdna2EnPgwAFUVVWhT58+SsyRI0dQXl6uxKSlpaFt27Zo1qyZauNTwy+//FLjB8hdXV1RVVUFgHN2P5pyjhrtb9auQ6ceAFu2bBG9Xi8bN26Uc+fOySuvvCK+vr42R4M+zGbMmCE+Pj5y6NAh+fHHH5XbL7/8osRMnz5dQkJC5MCBA3LixAkxmUxiMpmU/upTVYYNGyY5OTmSmpoqjz76aK2nqsyfP1/Onz8va9eufaBPVbmb9VHHIpyz2hw/flzc3NzkzTfflIsXL8onn3winp6e8j//8z9KTHx8vPj6+kpycrKcPn1aRo8eXetpGN27d5esrCw5evSotGnTxuY0jMLCQgkKCpJJkybJmTNnZMuWLeLp6fnAnKpiLSoqSh577DHl9J7t27dLQECALFiwQInhnN05A+DUqVNy6tQpASDvvvuunDp1Sr7//nsRabo5+vzzz8XNzU3efvttOX/+vCxbtoyn91Rbs2aNhISEiE6nk969e8uxY8ccnVKTAVDrLTExUYm5ffu2/PnPf5ZmzZqJp6enPP/88/Ljjz/arOe7776TESNGiIeHhwQEBMi8efOkvLzcJubgwYPSrVs30el00qpVK5vneNDdXWg5Z7XbvXu3dOrUSfR6vbRr107+8Y9/2PRXVVXJa6+9JkFBQaLX62XIkCFy4cIFm5iff/5ZXnjhBfHy8hKDwSAvv/yy3Lx50ybmq6++kn79+oler5fHHntM4uPjVR+bGoqLi2XOnDkSEhIi7u7u0qpVK3n11VdtTjHhnN35O6ntcywqKkpEmnaOtm3bJr/5zW9Ep9NJx44dJSUlxe7x8PdoiYiIVPRQfUdLRETkbFhoiYiIVMRCS0REpCIWWiIiIhWx0BIREamIhZaIiEhFLLREREQqYqElIiJSEQstERGRilhoiYiIVMRCS0REpKL/B00fZn32Fc0RAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 500x150 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(5, 1.5))\n",
    "ax.set_title(\"Map View of Model\", fontweight=\"bold\")\n",
    "mapview = fp.plot.PlotMapView(model=model, ax=ax)\n",
    "mapview.plot_ibound()\n",
    "mapview.plot_grid()\n",
    "\n",
    "chd_bc = mapview.plot_bc('CHD')\n",
    "wel_bc = mapview.plot_bc('WEL')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 2. Breakdown: simulation level packages"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2(a). [Temporal Discretization (TDIS) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mftdis.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = ex-gwf-fhb.tdis\n",
       "filename = ex-gwf-fhb.tdis\n",
       "package_type = tdis\n",
       "model_or_simulation_package = simulation\n",
       "simulation_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "time_units\n",
       "{internal}\n",
       "('days')\n",
       "\n",
       "\n",
       "Block dimensions\n",
       "--------------------\n",
       "nper\n",
       "{internal}\n",
       "(3)\n",
       "\n",
       "\n",
       "Block perioddata\n",
       "--------------------\n",
       "perioddata\n",
       "{internal}\n",
       "(rec.array([(400., 10, 1.), (200.,  4, 1.), (400.,  6, 1.)],\n",
       "          dtype=[('perlen', '<f8'), ('nstp', '<i4'), ('tsmult', '<f8')]))\n",
       "\n"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim.get_package('tdis')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "time_units = \"days\"\n",
    "nper = 3\n",
    "perioddata = [\n",
    "    # first stress period, transient\n",
    "    (400.0, 10, 1.0),\n",
    "    # second stress period, transient\n",
    "    (200.0, 4, 1.0),\n",
    "    # third stress period, transient\n",
    "    (400.0, 6, 1.0)\n",
    "]\n",
    "filename = \"ex-gwf-fhb.tdis\"\n",
    "pname = \"ex-gwf-fhb.tdis\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2(b). [Iterative Model Solution (IMS)](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfims.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = ex-gwf-fhb\n",
       "filename = ex-gwf-fhb.ims\n",
       "package_type = ims\n",
       "model_or_simulation_package = simulation\n",
       "simulation_name = ex-gwf-fhb\n",
       "\n",
       "Block nonlinear\n",
       "--------------------\n",
       "outer_dvclose\n",
       "{internal}\n",
       "(1e-09)\n",
       "\n",
       "outer_maximum\n",
       "{internal}\n",
       "(50)\n",
       "\n",
       "\n",
       "Block linear\n",
       "--------------------\n",
       "inner_maximum\n",
       "{internal}\n",
       "(100)\n",
       "\n",
       "inner_dvclose\n",
       "{internal}\n",
       "(1e-09)\n",
       "\n",
       "rcloserecord\n",
       "{internal}\n",
       "(rec.array([('inner_rclose', 1.e-06, 'strict')],\n",
       "          dtype=[('inner_rclose_label', 'O'), ('inner_rclose', '<f8'), ('rclose_option', 'O')]))\n",
       "\n"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sim.get_package('ims')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "outer_dvclose = 1e-9\n",
    "outer_maximum = 50\n",
    "inner_maximum = 100\n",
    "inner_dvclose = 1e-9\n",
    "rcloserecord = [1e-6, \"strict\"]\n",
    "filename = \"ex-gwf-fhb.ims\"\n",
    "pname = \"ex-gwf-fhb\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 3. Breakdown: model level packages"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['dis', 'npf', 'ic', 'sto', 'chd_ts', 'chd', 'wel_ts', 'wel', 'oc', 'obs_0']"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.package_names"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(a). [Structured Discretization (DIS)](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfdis.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = dis\n",
       "filename = ex-gwf-fhb.dis\n",
       "package_type = dis\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "length_units\n",
       "{internal}\n",
       "('meters')\n",
       "\n",
       "\n",
       "Block dimensions\n",
       "--------------------\n",
       "nlay\n",
       "{internal}\n",
       "(1)\n",
       "\n",
       "nrow\n",
       "{internal}\n",
       "(3)\n",
       "\n",
       "ncol\n",
       "{internal}\n",
       "(10)\n",
       "\n",
       "\n",
       "Block griddata\n",
       "--------------------\n",
       "delr\n",
       "{constant 1000.0}\n",
       "\n",
       "delc\n",
       "{constant 1000.0}\n",
       "\n",
       "top\n",
       "{constant 50.0}\n",
       "\n",
       "botm\n",
       "{constant -200.0}\n",
       "\n"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('dis')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "length_units = \"meters\"\n",
    "nlay = 1\n",
    "nrow = 3\n",
    "ncol = 10\n",
    "delr = 1000.0\n",
    "delc = 1000.0\n",
    "top = 50.0\n",
    "botm = -200.0\n",
    "filename = \"ex-gwf-fhb.dis\"\n",
    "pname = \"dis\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(b). [Node Property Flow (NPF) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfnpf.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = npf\n",
       "filename = ex-gwf-fhb.npf\n",
       "package_type = npf\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "save_specific_discharge\n",
       "{internal}\n",
       "(True)\n",
       "\n",
       "\n",
       "Block griddata\n",
       "--------------------\n",
       "icelltype\n",
       "{constant 0}\n",
       "\n",
       "k\n",
       "{constant 20.0}\n",
       "\n"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('npf')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "save_specific_discharge = True\n",
    "# layer is confined\n",
    "icelltype = 0\n",
    "k = 20.0\n",
    "filename = \"ex-gwf-fhb.npf\"\n",
    "pname = \"npf\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(c). [Initial Conditions (IC) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfic.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = ic\n",
       "filename = ex-gwf-fhb.ic\n",
       "package_type = ic\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block griddata\n",
       "--------------------\n",
       "strt\n",
       "{constant 0.0}\n",
       "\n"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('ic')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "# starting head\n",
    "strt = 0.0\n",
    "filename = \"ex-gwf-fhb.ic\"\n",
    "pname = \"ic\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(d). [Storage (STO) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfsto.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = sto\n",
       "filename = ex-gwf-fhb.sto\n",
       "package_type = sto\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "storagecoefficient\n",
       "{internal}\n",
       "(True)\n",
       "\n",
       "\n",
       "Block griddata\n",
       "--------------------\n",
       "iconvert\n",
       "{constant 0}\n",
       "\n",
       "ss\n",
       "{constant 1e-06}\n",
       "\n",
       "\n",
       "Block period\n",
       "--------------------\n",
       "steady-state\n",
       "None\n",
       "transient\n",
       "{internal}\n",
       "(True)\n",
       "\n"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('sto')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "# SS array is read as storage coefficient rather than specific storage\n",
    "storagecoefficient = True\n",
    "# confined storage is used\n",
    "iconvert = 0\n",
    "# storage coefficient (since flag is set)\n",
    "ss = 1e-6\n",
    "transient = {0: True}\n",
    "filename = \"ex-gwf-fhb.sto\"\n",
    "pname = \"sto\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(e). [Constant-Head (CHD) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfchd.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = chd\n",
       "filename = ex-gwf-fhb.chd\n",
       "package_type = chd\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "ts_filerecord\n",
       "{internal}\n",
       "(rec.array([('ex-gwf-fhb.chd.ts',)],\n",
       "          dtype=[('ts6_filename', 'O')]))\n",
       "\n",
       "\n",
       "Block dimensions\n",
       "--------------------\n",
       "maxbound\n",
       "{internal}\n",
       "(3)\n",
       "\n",
       "\n",
       "Block period\n",
       "--------------------\n",
       "stress_period_data\n",
       "{0: <flopy.mf6.data.mfdataplist.PandasListStorage object at 0x000001BACF825F90>}\n"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('chd')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "rec.array([((0, 0, 9), 'chdhead'), ((0, 1, 9), 'chdhead'),\n",
       "           ((0, 2, 9), 'chdhead')],\n",
       "          dtype=[('cellid', 'O'), ('head', 'O')])"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('chd').stress_period_data.get_data(key=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "For the constant-head package...\n",
    "\n",
    "```python\n",
    "maxbound = 3\n",
    "# list of ((layer, row, column), head)\n",
    "stress_period_data = [\n",
    "    ((0, 0, 9), \"chdhead\"),\n",
    "    ((0, 1, 9), \"chdhead\"),\n",
    "    ((0, 2, 9), \"chdhead\")\n",
    "]\n",
    "filename = \"ex-gwf-fhb.chd\"\n",
    "pname = \"chd\"\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = chd_ts\n",
       "filename = ex-gwf-fhb.chd.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = chd\n",
       "\n",
       "Block attributes\n",
       "--------------------\n",
       "time_series_namerecord\n",
       "{internal}\n",
       "(rec.array([('chdhead',)],\n",
       "          dtype=[('time_series_names', 'O')]))\n",
       "\n",
       "interpolation_methodrecord\n",
       "{internal}\n",
       "(rec.array([('linearend',)],\n",
       "          dtype=[('interpolation_method', 'O')]))\n",
       "\n",
       "\n",
       "Block timeseries\n",
       "--------------------\n",
       "timeseries\n",
       "{internal}\n",
       "(rec.array([(   0., 0.), ( 307., 1.), ( 791., 5.), (1000., 2.)],\n",
       "          dtype=[('ts_time', '<f8'), ('ts_array', '<f8')]))\n",
       "\n"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('chd_ts')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "For the [time series](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfutlts.html) attached to the constant-head package...\n",
    "\n",
    "```python\n",
    "time_series_namerecord = [\"chdhead\"]\n",
    "interpolation_methodrecord = [\"linearend\"]\n",
    "# list of (ts_time, ts_array)\n",
    "timeseries = [\n",
    "    (0.0, 0.0),\n",
    "    (307.0, 1.0),\n",
    "    (791.0, 5.0),\n",
    "    (1000.0, 2.0)\n",
    "]\n",
    "filename = \"ex-gwf-fhb.chd.ts\"\n",
    "pname = \"chd_ts\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(f). [Well (WEL) Package](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfwel.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = wel\n",
       "filename = ex-gwf-fhb.wel\n",
       "package_type = wel\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "ts_filerecord\n",
       "{internal}\n",
       "(rec.array([('ex-gwf-fhb.wel.ts',)],\n",
       "          dtype=[('ts6_filename', 'O')]))\n",
       "\n",
       "\n",
       "Block dimensions\n",
       "--------------------\n",
       "maxbound\n",
       "{internal}\n",
       "(1)\n",
       "\n",
       "\n",
       "Block period\n",
       "--------------------\n",
       "stress_period_data\n",
       "{0: <flopy.mf6.data.mfdataplist.PandasListStorage object at 0x000001BACF834B10>}\n"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('wel')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "rec.array([((0, 1, 0), 'flowrate')],\n",
       "          dtype=[('cellid', 'O'), ('q', 'O')])"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('wel').stress_period_data.get_data(key=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "For the well package...\n",
    "\n",
    "```python\n",
    "maxbound = 1\n",
    "# list of ((layer, row, column), q)\n",
    "stress_period_data = [\n",
    "    ((0, 1, 0), \"flowrate\")\n",
    "]\n",
    "filename = \"ex-gwf-fhb.wel\"\n",
    "pname = \"wel\"\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = wel_ts\n",
       "filename = ex-gwf-fhb.wel.ts\n",
       "package_type = ts\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "parent_file = wel\n",
       "\n",
       "Block attributes\n",
       "--------------------\n",
       "time_series_namerecord\n",
       "{internal}\n",
       "(rec.array([('flowrate',)],\n",
       "          dtype=[('time_series_names', 'O')]))\n",
       "\n",
       "interpolation_methodrecord\n",
       "{internal}\n",
       "(rec.array([('linearend',)],\n",
       "          dtype=[('interpolation_method', 'O')]))\n",
       "\n",
       "\n",
       "Block timeseries\n",
       "--------------------\n",
       "timeseries\n",
       "{internal}\n",
       "(rec.array([(   0., 2000.), ( 307., 6000.), ( 791., 5000.), (1000., 9000.)],\n",
       "          dtype=[('ts_time', '<f8'), ('ts_array', '<f8')]))\n",
       "\n"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('wel_ts')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "For the [time series](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfutlts.html) attached to the well package...\n",
    "\n",
    "```python\n",
    "time_series_namerecord = [\"flowrate\"]\n",
    "interpolation_methodrecord = [\"linearend\"]\n",
    "# list of (ts_time, ts_array)\n",
    "timeseries = [\n",
    "    (0.0, 2000.0),\n",
    "    (307.0, 6000.0),\n",
    "    (791.0, 5000.0),\n",
    "    (1000.0, 9000.0)\n",
    "]\n",
    "filename = \"ex-gwf-fhb.wel.ts\"\n",
    "pname = \"wel_ts\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(g). [Output Control (OC)](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfgwfoc.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = oc\n",
       "filename = ex-gwf-fhb.oc\n",
       "package_type = oc\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block options\n",
       "--------------------\n",
       "budget_filerecord\n",
       "{internal}\n",
       "(rec.array([('ex-gwf-fhb.cbc',)],\n",
       "          dtype=[('budgetfile', 'O')]))\n",
       "\n",
       "head_filerecord\n",
       "{internal}\n",
       "(rec.array([('ex-gwf-fhb.hds',)],\n",
       "          dtype=[('headfile', 'O')]))\n",
       "\n",
       "\n",
       "Block period\n",
       "--------------------\n",
       "saverecord\n",
       "{internal}\n",
       "(rec.array([('head', 'all', None), ('budget', 'all', None)],\n",
       "          dtype=[('rtype', 'O'), ('ocsetting', 'O'), ('ocsetting_data', 'O')]))\n",
       "\n",
       "printrecord\n",
       "None\n"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('oc')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "```python\n",
    "budget_filerecord = [\"ex-gwf-fhb.cbc\"]\n",
    "head_filerecord = [\"ex-gwf-fhb.hds\"]\n",
    "saverecord = [\n",
    "    (\"head\", \"all\"),\n",
    "    (\"budget\", \"all\")\n",
    "]\n",
    "printrecord = None # default value\n",
    "filename = \"ex-gwf-fhb.oc\"\n",
    "pname = \"oc\"\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3(h). [Observation (OBS)](https://flopy.readthedocs.io/en/latest/source/flopy.mf6.modflow.mfutlobs.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "package_name = obs_0\n",
       "filename = ex-gwf-fhb.obs\n",
       "package_type = obs\n",
       "model_or_simulation_package = model\n",
       "model_name = ex-gwf-fhb\n",
       "\n",
       "Block continuous\n",
       "--------------------\n",
       "continuous\n",
       "{internal}\n",
       "(rec.array([('h1_2_1', 'head', (0, 1, 0), None),\n",
       "           ('h1_2_10', 'head', (0, 1, 9), None)],\n",
       "          dtype=[('obsname', 'O'), ('obstype', 'O'), ('id', 'O'), ('id2', 'O')]))\n",
       "\n"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.get_package('obs_0')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Reverse-engineer\n",
    "\n",
    "There is also a flow observation in the `ex-gwf-fhb.obs` file.\n",
    "\n",
    "```python\n",
    "continuous = {\n",
    "    \"ex-gwf-fhb.obs.flow.csv\": [\n",
    "        (\"icf1\", \"flow-ja-face\", (0, 1, 1), (0, 1, 0))\n",
    "    ],\n",
    "    \"ex-gwf-fhb.obs.head.csv\": [\n",
    "        (\"h1_2_1\", \"head\", (0, 1, 0), None),\n",
    "        (\"h1_2_10\", \"head\", (0, 1, 9), None)\n",
    "    ]\n",
    "}\n",
    "filename = \"ex-gwf-fhb.obs\"\n",
    "pname = \"obs_0\"\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.11.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}