test_mf6_load.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "import flopy\n",
    "\n",
    "from IPython.core.interactiveshell import InteractiveShell\n",
    "InteractiveShell.ast_node_interactivity = \"all\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "sim_name = \"mfsim\"\n",
    "sim_path = \"./\"\n",
    "model_name = \"test_mf6_algomes\" #ha!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "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 disu...\n",
      "WARNING: Unable to resolve dimension of ('gwf6', 'disu', 'cell2d', 'cell2d', 'icvert') based on shape \"ncvert\".\n",
      "    loading package ic...\n",
      "    loading package oc...\n",
      "    loading package npf...\n",
      "  loading ims package test_mf6_algomes...\n"
     ]
    }
   ],
   "source": [
    "# load the simulation\n",
    "sim = flopy.mf6.MFSimulation.load(\n",
    "    sim_name, version=\"mf6\", exe_name=\"mf6\", sim_ws=sim_path)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "name = test_mf6_algomes\n",
       "model_type = gwf\n",
       "version = mf6\n",
       "model_relative_path = \n",
       "\n",
       "###################\n",
       "Package disu\n",
       "###################\n",
       "\n",
       "package_name = disu\n",
       "filename = test_mf6_algomesh.disu\n",
       "package_type = disu\n",
       "model_or_simulation_package = model\n",
       "model_name = test_mf6_algomes\n",
       "\n",
       "\n",
       "###################\n",
       "Package ic\n",
       "###################\n",
       "\n",
       "package_name = ic\n",
       "filename = test_mf6_algomesh.ic\n",
       "package_type = ic\n",
       "model_or_simulation_package = model\n",
       "model_name = test_mf6_algomes\n",
       "\n",
       "\n",
       "###################\n",
       "Package oc\n",
       "###################\n",
       "\n",
       "package_name = oc\n",
       "filename = test_mf6_algomesh.oc\n",
       "package_type = oc\n",
       "model_or_simulation_package = model\n",
       "model_name = test_mf6_algomes\n",
       "\n",
       "\n",
       "###################\n",
       "Package npf\n",
       "###################\n",
       "\n",
       "package_name = npf\n",
       "filename = test_mf6_algomesh.npf\n",
       "package_type = npf\n",
       "model_or_simulation_package = model\n",
       "model_name = test_mf6_algomes\n",
       "\n"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "ename": "KeyError",
     "evalue": "None",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyError\u001b[0m                                  Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-23-8196b42f6238>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[0;32m      1\u001b[0m \u001b[0mml\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0msim\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_model\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmodel_name\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      2\u001b[0m \u001b[0mml\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mml\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdisu\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtop\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\mf6\\data\\mfdataarray.py\u001b[0m in \u001b[0;36mplot\u001b[1;34m(self, filename_base, file_extension, mflay, fignum, title, **kwargs)\u001b[0m\n\u001b[0;32m   1414\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1415\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mnum_plottable_layers\u001b[0m \u001b[1;33m==\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1416\u001b[1;33m             axes = PlotUtilities._plot_util2d_helper(\n\u001b[0m\u001b[0;32m   1417\u001b[0m                 \u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1418\u001b[0m                 \u001b[0mtitle\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mtitle\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\plot\\plotutil.py\u001b[0m in \u001b[0;36m_plot_util2d_helper\u001b[1;34m(util2d, title, filename_base, file_extension, fignum, **kwargs)\u001b[0m\n\u001b[0;32m    744\u001b[0m             \u001b[0mfilename\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;34mf\"{filename_base}_{util2d.name}.{fext}\"\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    745\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 746\u001b[1;33m         axes = PlotUtilities._plot_array_helper(\n\u001b[0m\u001b[0;32m    747\u001b[0m             \u001b[0mutil2d\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0marray\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    748\u001b[0m             \u001b[0mutil2d\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmodel\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\plot\\plotutil.py\u001b[0m in \u001b[0;36m_plot_array_helper\u001b[1;34m(plotarray, model, modelgrid, axes, names, filenames, fignum, mflay, **kwargs)\u001b[0m\n\u001b[0;32m   1158\u001b[0m             )\n\u001b[0;32m   1159\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mdefaults\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m\"pcolor\"\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1160\u001b[1;33m                 cm = pmv.plot_array(\n\u001b[0m\u001b[0;32m   1161\u001b[0m                     \u001b[0mplotarray\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1162\u001b[0m                     \u001b[0mmasked_values\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mdefaults\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;34m\"masked_values\"\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\plot\\map.py\u001b[0m in \u001b[0;36mplot_array\u001b[1;34m(self, a, masked_values, **kwargs)\u001b[0m\n\u001b[0;32m    114\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    115\u001b[0m         \u001b[1;31m# use cached patch collection for plotting\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 116\u001b[1;33m         \u001b[0mpolygons\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmg\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmap_polygons\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    117\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mpolygons\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mdict\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    118\u001b[0m             \u001b[0mpolygons\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mpolygons\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlayer\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\discretization\\unstructuredgrid.py\u001b[0m in \u001b[0;36mmap_polygons\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m    479\u001b[0m             \u001b[1;32mor\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_cache_dict\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcache_index\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mout_of_date\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    480\u001b[0m         ):\n\u001b[1;32m--> 481\u001b[1;33m             \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mxyzvertices\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    482\u001b[0m             \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_polygons\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;32mNone\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    483\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\discretization\\unstructuredgrid.py\u001b[0m in \u001b[0;36mxyzvertices\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m    321\u001b[0m             \u001b[1;32mor\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_cache_dict\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcache_index\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mout_of_date\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    322\u001b[0m         ):\n\u001b[1;32m--> 323\u001b[1;33m             \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_build_grid_geometry_info\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    324\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_copy_cache\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    325\u001b[0m             \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_cache_dict\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcache_index\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdata\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mc:\\python39\\lib\\site-packages\\flopy\\discretization\\unstructuredgrid.py\u001b[0m in \u001b[0;36m_build_grid_geometry_info\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m    565\u001b[0m             \u001b[0mycellvert\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    566\u001b[0m             \u001b[1;32mfor\u001b[0m \u001b[0mix\u001b[0m \u001b[1;32min\u001b[0m \u001b[0miverts\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 567\u001b[1;33m                 \u001b[0mxcellvert\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mvertexdict\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mix\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    568\u001b[0m                 \u001b[0mycellvert\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mvertexdict\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mix\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    569\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mKeyError\u001b[0m: None"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQYAAAEICAYAAAC9P1pMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAANtElEQVR4nO3ccaydd13H8feHlYqOMZReEmw71mknVDRh3swpicxsajdC+wcG2zgVU9eAjGggmpEtkxSNQRQTtIolzsESNgp/kGsomYJbliCF3mUw1i4jl9LRWwgrsC0xC9savv5xTvX0/u7tedqee07veL+SJs/znN8953tPbt57znOfu1QVkjToBZMeQNL5xzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMAiDJkSTXTnoOnR8Mg6SGYRBJ7gQuAf49yf8k+fMkW5IcTPJkkvuSvHpg/ZEk705yKMkTSf4tyYsm9x1o1AyDqKrfA74JvLGqXgx8CrgL+FNgCthHLxqrB77sd4HfAn4GuBy4dYwja5kZBi3md4BPV9V/VtVzwN8CPw786sCaf6yqo1X1feCvgO0TmFPLxDBoMT8NPHZyp6p+CBwF1g6sOTqw/Vj/a/Q8YRh00uCf2X4LeOXJnSQB1gPHBtasH9i+pP81ep4wDDrpO8Bl/e29wBuSXJPkhcC7gGeA/x5Y//Yk65L8FHAL8PGxTqtlZRh00l8DtyZ5EngjcAPwD8B3+/tvrKpnB9Z/DPgP4DDwdeAvxzqtllX8H7XoTCU5AvxRVX120rNoeXjGIKkxNAxJbk/yeJKHl3g8ST6YZC7JQ0muGP2YksapyxnDHcDm0zx+HbCx/28n8M/nPpbOZ1V1qR8jnt+GhqGq7ge+f5olW4GPVs9+4KVJXjGqASWN36oRPMdaTr3ZZb5/7NsLFybZSe+sggsvvPCXXvWqV43g5SUt5YEHHvhuVU2d6deNIgydVdUeYA/A9PR0zc7OjvPlpR85SR4bvqo1it9KHOPUu+DWceodcpJWmFGEYQb4/f5vJ64Cnqqq5mOEpJVj6EeJJHcBVwNrkswDfwG8EKCqPkTvT3KvB+aAp4E/XK5hJY3H0DBU1Wn/nLZ6t06+fWQTSZo473yU1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUqNTGJJsTvJokrkkNy/y+CVJ7k3yYJKHklw/+lEljcvQMCS5ANgNXAdsArYn2bRg2a3A3qp6LbAN+KdRDyppfLqcMVwJzFXV4ap6Frgb2LpgTQEv6W9fDHxrdCNKGrcuYVgLHB3Yn+8fG/Qe4IYk88A+4B2LPVGSnUlmk8weP378LMaVNA6juvi4HbijqtYB1wN3Jmmeu6r2VNV0VU1PTU2N6KUljVqXMBwD1g/sr+sfG7QD2AtQVV8AXgSsGcWAksavSxgOABuTbEiymt7FxZkFa74JXAOQ5NX0wuBnBWmFGhqGqjoB3ATcAzxC77cPB5PsSrKlv+xdwI1JvgLcBbylqmq5hpa0vFZ1WVRV++hdVBw8dtvA9iHgdaMdTdKkeOejpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1OgUhiSbkzyaZC7JzUuseXOSQ0kOJvnYaMeUNE6rhi1IcgGwG/gNYB44kGSmqg4NrNkIvBt4XVU9keTlyzWwpOXX5YzhSmCuqg5X1bPA3cDWBWtuBHZX1RMAVfX4aMeUNE5dwrAWODqwP98/Nuhy4PIkn0+yP8nmxZ4oyc4ks0lmjx8/fnYTS1p2o7r4uArYCFwNbAc+nOSlCxdV1Z6qmq6q6ampqRG9tKRR6xKGY8D6gf11/WOD5oGZqnquqr4BfI1eKCStQF3CcADYmGRDktXANmBmwZpP0TtbIMkaeh8tDo9uTEnjNDQMVXUCuAm4B3gE2FtVB5PsSrKlv+we4HtJDgH3An9WVd9brqElLa9U1UReeHp6umZnZyfy2tKPiiQPVNX0mX6ddz5KahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqdEpDEk2J3k0yVySm0+z7k1JKsn06EaUNG5Dw5DkAmA3cB2wCdieZNMi6y4C/gT44qiHlDReXc4YrgTmqupwVT0L3A1sXWTde4H3AT8Y4XySJqBLGNYCRwf25/vH/k+SK4D1VfXp0z1Rkp1JZpPMHj9+/IyHlTQe53zxMckLgA8A7xq2tqr2VNV0VU1PTU2d60tLWiZdwnAMWD+wv65/7KSLgNcA9yU5AlwFzHgBUlq5uoThALAxyYYkq4FtwMzJB6vqqapaU1WXVtWlwH5gS1XNLsvEkpbd0DBU1QngJuAe4BFgb1UdTLIryZblHlDS+K3qsqiq9gH7Fhy7bYm1V5/7WJImyTsfJTUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpEanMCTZnOTRJHNJbl7k8XcmOZTkoSSfS/LK0Y8qaVyGhiHJBcBu4DpgE7A9yaYFyx4EpqvqF4FPAn8z6kEljU+XM4YrgbmqOlxVzwJ3A1sHF1TVvVX1dH93P7ButGNKGqcuYVgLHB3Yn+8fW8oO4DOLPZBkZ5LZJLPHjx/vPqWksRrpxcckNwDTwPsXe7yq9lTVdFVNT01NjfKlJY3Qqg5rjgHrB/bX9Y+dIsm1wC3A66vqmdGMJ2kSupwxHAA2JtmQZDWwDZgZXJDktcC/AFuq6vHRjylpnIaGoapOADcB9wCPAHur6mCSXUm29Je9H3gx8IkkX04ys8TTSVoBunyUoKr2AfsWHLttYPvaEc8laYK881FSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6SGYZDUMAySGoZBUsMwSGoYBkkNwyCpYRgkNQyDpIZhkNQwDJIahkFSwzBIahgGSQ3DIKlhGCQ1DIOkhmGQ1DAMkhqGQVLDMEhqGAZJjU5hSLI5yaNJ5pLcvMjjP5bk4/3Hv5jk0pFPKmlshoYhyQXAbuA6YBOwPcmmBct2AE9U1c8Cfw+8b9SDShqfLmcMVwJzVXW4qp4F7ga2LlizFfhIf/uTwDVJMroxJY3Tqg5r1gJHB/bngV9eak1VnUjyFPAy4LuDi5LsBHb2d59J8vDZDD0ha1jw/ZzHVtKssLLmXUmzAvzc2XxRlzCMTFXtAfYAJJmtqulxvv65WEnzrqRZYWXNu5Jmhd68Z/N1XT5KHAPWD+yv6x9bdE2SVcDFwPfOZiBJk9clDAeAjUk2JFkNbANmFqyZAf6gv/3bwH9VVY1uTEnjNPSjRP+awU3APcAFwO1VdTDJLmC2qmaAfwXuTDIHfJ9ePIbZcw5zT8JKmnclzQora96VNCuc5bzxP+ySFvLOR0kNwyCpsexhWEm3U3eY9Z1JDiV5KMnnkrxyEnMOzHPaeQfWvSlJJZnYr9m6zJrkzf3392CSj417xgWzDPtZuCTJvUke7P88XD+JOfuz3J7k8aXuC0rPB/vfy0NJrhj6pFW1bP/oXaz8OnAZsBr4CrBpwZo/Bj7U394GfHw5ZzrHWX8d+In+9tsmNWvXefvrLgLuB/YD0+frrMBG4EHgJ/v7Lz+f31t6F/Xe1t/eBByZ4Ly/BlwBPLzE49cDnwECXAV8cdhzLvcZw0q6nXrorFV1b1U93d/dT++ejknp8t4CvJfe3678YJzDLdBl1huB3VX1BEBVPT7mGQd1mbeAl/S3Lwa+Ncb5Th2k6n56vw1cylbgo9WzH3hpklec7jmXOwyL3U69dqk1VXUCOHk79bh1mXXQDnoVnpSh8/ZPGddX1afHOdgiury3lwOXJ/l8kv1JNo9tulaXed8D3JBkHtgHvGM8o52VM/3ZHu8t0c8XSW4ApoHXT3qWpSR5AfAB4C0THqWrVfQ+TlxN70zs/iS/UFVPTnKo09gO3FFVf5fkV+jdx/OaqvrhpAcbheU+Y1hJt1N3mZUk1wK3AFuq6pkxzbaYYfNeBLwGuC/JEXqfLWcmdAGyy3s7D8xU1XNV9Q3ga/RCMQld5t0B7AWoqi8AL6L3B1bno04/26dY5osiq4DDwAb+/yLOzy9Y83ZOvfi4d0IXcLrM+lp6F6U2TmLGM513wfr7mNzFxy7v7WbgI/3tNfROfV92Hs/7GeAt/e1X07vGkAn+PFzK0hcf38CpFx+/NPT5xjDw9fTq/3Xglv6xXfT+iwu90n4CmAO+BFw2wTd32KyfBb4DfLn/b2ZSs3aZd8HaiYWh43sbeh99DgFfBbadz+8tvd9EfL4fjS8DvznBWe8Cvg08R+/MawfwVuCtA+/t7v738tUuPwfeEi2p4Z2PkhqGQVLDMEhqGAZJDcMgqWEYJDUMg6TG/wIML6lP9o3LUQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ml = sim.get_model(model_name)\n",
    "ml\n",
    "ml.disu.top.plot()"
   ]
  }
 ],
 "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.9.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}