Leapfrog_load.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b9afb838-081f-4893-a80d-ed78d1e3fb94",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "gridgen executable was found at: C:\\mffiles\\gridgen.EXE\n",
      "Pandas version = 2.1.1\n",
      "Numpy version = 1.26.0\n",
      "Flopy version = 3.4.3\n",
      "Matplotlib version = 3.8.0\n"
     ]
    }
   ],
   "source": [
    "import os\n",
    "import sys\n",
    "import shutil\n",
    "import platform\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "import flopy\n",
    "from flopy.discretization import VertexGrid\n",
    "from flopy.utils import Raster\n",
    "from flopy.utils import GridIntersect\n",
    "from flopy.utils.gridgen import Gridgen\n",
    "gridgen_exe = \"gridgen\"\n",
    "if platform.system() in \"Windows\":\n",
    "    gridgen_exe += \".exe\"\n",
    "gridgen_exe = flopy.which(\"gridgen\")\n",
    "if gridgen_exe is None:\n",
    "    msg = (\n",
    "        \"Warning, gridgen is not in your path. \"\n",
    "        \"When you create the griden object you will need to \"\n",
    "        \"provide a full path to the gridgen binary executable.\"\n",
    "    )\n",
    "    print(msg)\n",
    "else:\n",
    "    print(\"gridgen executable was found at: {}\".format(gridgen_exe))\n",
    "\n",
    "print(f\"Pandas version = {pd.__version__}\")\n",
    "print(f\"Numpy version = {np.__version__}\")\n",
    "print(f\"Flopy version = {flopy.__version__}\")\n",
    "print(f\"Matplotlib version = {matplotlib.__version__}\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "58bc833a-910b-4ea5-b4c9-d4de7738dd42",
   "metadata": {},
   "outputs": [],
   "source": [
    "ws = os.path.join('Leapfrog_load') # here we are making a path not creating the folder\n",
    "gis_f = os.path.join(ws,'GIS') # creating a sub-directory path for our gis input/output\n",
    "model_f = os.path.join(ws,'model') # creating a sub-directory path for our model input/output\n",
    "plots_f = os.path.join(ws,'plots') # creating a sub-directory path for our plots\n",
    "for path in [ws,gis_f,model_f,plots_f]:\n",
    "    if os.path.exists(path): # here we are asking if the path exists on the computer. \n",
    "        shutil.rmtree(path)# if it does exist, delete it and all the files in it\n",
    "        os.mkdir(path) # then remake it\n",
    "    else:\n",
    "        os.mkdir(path) # if it doesn't exist then make the folder"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "c78f8903-ea94-4faa-be8a-431fcf0c118b",
   "metadata": {},
   "outputs": [],
   "source": [
    "path_to_files = os.path.join('files','Leapfrog_Export')\n",
    "flist = [x for x in os.listdir(path_to_files)]\n",
    "_ = [shutil.copyfile(os.path.join(path_to_files,file),os.path.join(model_f,file)) for file in flist]    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "0d7da0d2-8c09-4f6a-b046-bfa86047131d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['Leapfrog_Export.bas6',\n",
       " 'Leapfrog_Export.dis',\n",
       " 'Leapfrog_Export.list',\n",
       " 'Leapfrog_Export.lpf',\n",
       " 'Leapfrog_Export.nam',\n",
       " 'Leapfrog_Export.zon',\n",
       " 'Leapfrog_Export_conc.txt',\n",
       " 'Leapfrog_Export_conc.ucn']"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "flist"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "390ee38b-e1ee-43bb-bf6a-81e3d62b8e0b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LIST 10 WB19BS1E0003_PT.list\n",
      "DIS 11 WB19BS1E0003_PT.dis\n",
      "BAS6 12 WB19BS1E0003_PT.bas6\n",
      "LPF 13 WB19BS1E0003_PT.lpf\n",
      "ZONE 14 WB19BS1E0003_PT.zon\n"
     ]
    }
   ],
   "source": [
    "# Now look at the flow.dis input file\n",
    "nam_file = 'Leapfrog_Export.nam'\n",
    "_ = [print(line.rstrip()) for line in open(os.path.join(model_f,nam_file))]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1823be1e-e842-4c0b-8bd5-394c78bc96bc",
   "metadata": {},
   "source": [
    "The filenames in the .nam file do not match the rest of the names of the files in the directory so this load will fail. You will have to rename the files in the namfile to match those in the rest of the directory in order to get it to load. Unsure why Leapfrog did/does this. Or you could rename the files in the directory to match those listed in the namfile."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "b8ae96f1-7c41-476c-a58c-513bccf71562",
   "metadata": {},
   "outputs": [],
   "source": [
    "new_nam = 'WB19BS1E0003_PT'\n",
    "_ = [os.rename(os.path.join(model_f,file),os.path.join(model_f,new_nam+\".\"+file.split(\".\")[-1])) for file in flist]    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "bd994062-6dd8-4c72-9026-420d45ac4859",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Creating new model with name: WB19BS1E0003_PT\n",
      "--------------------------------------------------\n",
      "\n",
      "Parsing the namefile --> c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.nam\n",
      "\n",
      "--------------------------------------------------\n",
      "External unit dictionary:\n",
      "{10: filename:c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.list, filetype:LIST, 11: filename:c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.dis, filetype:DIS, 12: filename:c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.bas6, filetype:BAS6, 13: filename:c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.lpf, filetype:LPF, 14: filename:c:\\Working\\NCGRT_IntroFlopy\\Leapfrog_load\\model\\WB19BS1E0003_PT.zon, filetype:ZONE}\n",
      "--------------------------------------------------\n",
      "\n",
      "ModflowBas6 free format:True\n",
      "\n",
      "loading dis package file...\n",
      "   Loading dis package with:\n",
      "      17 layers, 50 rows, 75 columns, and 1 stress periods\n",
      "   loading laycbd...\n",
      "   loading delr...\n",
      "   loading delc...\n",
      "   loading top...\n",
      "   loading botm...\n",
      "      for 17 layers and 0 confining beds\n",
      "   loading stress period data...\n",
      "       for 1 stress periods\n",
      "adding Package:  DIS\n",
      "   DIS  package load...success\n",
      "loading zone package file...\n",
      "   reading data for \"litho1    \" zone\n",
      "   reading data for \"litho2    \" zone\n",
      "   reading data for \"litho3    \" zone\n",
      "   reading data for \"litho4    \" zone\n",
      "   reading data for \"litho5    \" zone\n",
      "   reading data for \"litho6    \" zone\n",
      "   reading data for \"litho7    \" zone\n",
      "   reading data for \"litho8    \" zone\n",
      "   reading data for \"litho9    \" zone\n",
      "   reading data for \"litho10   \" zone\n",
      "   reading data for \"litho11   \" zone\n",
      "   reading data for \"litho12   \" zone\n",
      "   reading data for \"litho13   \" zone\n",
      "   reading data for \"litho14   \" zone\n",
      "   reading data for \"litho15   \" zone\n",
      "   reading data for \"litho16   \" zone\n",
      "   reading data for \"litho17   \" zone\n",
      "adding Package:  ZONE\n",
      "   ZONE package load...success\n",
      "removing Package:  ['ZONE']\n",
      "   LIST package load...skipped\n",
      "loading bas6 package file...\n",
      "adding Package:  BAS6\n",
      "   BAS6 package load...success\n",
      "loading lpf package file...\n",
      "   loading IBCFCB, HDRY, NPLPF...\n",
      "   loading LAYTYP...\n",
      "   loading LAYAVG...\n",
      "   loading CHANI...\n",
      "   loading LAYVKA...\n",
      "   loading LAYWET...\n",
      "   loading parameter \"kx1\"...\n",
      "   loading parameter \"kx2\"...\n",
      "   loading parameter \"kx3\"...\n",
      "   loading parameter \"kx4\"...\n",
      "   loading parameter \"kx5\"...\n",
      "   loading parameter \"kx6\"...\n",
      "   loading parameter \"kx7\"...\n",
      "   loading parameter \"kx8\"...\n",
      "   loading parameter \"kx9\"...\n",
      "   loading parameter \"kx10\"...\n",
      "   loading parameter \"kx11\"...\n",
      "   loading parameter \"kx12\"...\n",
      "   loading parameter \"kx13\"...\n",
      "   loading parameter \"kx14\"...\n",
      "   loading parameter \"kx15\"...\n",
      "   loading parameter \"kx16\"...\n",
      "   loading parameter \"kx17\"...\n",
      "   loading parameter \"kz1\"...\n",
      "   loading parameter \"kz2\"...\n",
      "   loading parameter \"kz3\"...\n",
      "   loading parameter \"kz4\"...\n",
      "   loading parameter \"kz5\"...\n",
      "   loading parameter \"kz6\"...\n",
      "   loading parameter \"kz7\"...\n",
      "   loading parameter \"kz8\"...\n",
      "   loading parameter \"kz9\"...\n",
      "   loading parameter \"kz10\"...\n",
      "   loading parameter \"kz11\"...\n",
      "   loading parameter \"kz12\"...\n",
      "   loading parameter \"kz13\"...\n",
      "   loading parameter \"kz14\"...\n",
      "   loading parameter \"kz15\"...\n",
      "   loading parameter \"kz16\"...\n",
      "   loading parameter \"kz17\"...\n",
      "   loading parameter \"ss1\"...\n",
      "   loading parameter \"ss2\"...\n",
      "   loading parameter \"ss3\"...\n",
      "   loading parameter \"ss4\"...\n",
      "   loading parameter \"ss5\"...\n",
      "   loading parameter \"ss6\"...\n",
      "   loading parameter \"ss7\"...\n",
      "   loading parameter \"ss8\"...\n",
      "   loading parameter \"ss9\"...\n",
      "   loading parameter \"ss10\"...\n",
      "   loading parameter \"ss11\"...\n",
      "   loading parameter \"ss12\"...\n",
      "   loading parameter \"ss13\"...\n",
      "   loading parameter \"ss14\"...\n",
      "   loading parameter \"ss15\"...\n",
      "   loading parameter \"ss16\"...\n",
      "   loading parameter \"ss17\"...\n",
      "   loading parameter \"sy1\"...\n",
      "   loading parameter \"sy2\"...\n",
      "   loading parameter \"sy3\"...\n",
      "   loading parameter \"sy4\"...\n",
      "   loading parameter \"sy5\"...\n",
      "   loading parameter \"sy6\"...\n",
      "   loading parameter \"sy7\"...\n",
      "   loading parameter \"sy8\"...\n",
      "   loading parameter \"sy9\"...\n",
      "   loading parameter \"sy10\"...\n",
      "   loading parameter \"sy11\"...\n",
      "   loading parameter \"sy12\"...\n",
      "   loading parameter \"sy13\"...\n",
      "   loading parameter \"sy14\"...\n",
      "   loading parameter \"sy15\"...\n",
      "   loading parameter \"sy16\"...\n",
      "   loading parameter \"sy17\"...\n",
      "   loading parameter \"an1\"...\n",
      "   loading parameter \"an2\"...\n",
      "   loading parameter \"an3\"...\n",
      "   loading parameter \"an4\"...\n",
      "   loading parameter \"an5\"...\n",
      "   loading parameter \"an6\"...\n",
      "   loading parameter \"an7\"...\n",
      "   loading parameter \"an8\"...\n",
      "   loading parameter \"an9\"...\n",
      "   loading parameter \"an10\"...\n",
      "   loading parameter \"an11\"...\n",
      "   loading parameter \"an12\"...\n",
      "   loading parameter \"an13\"...\n",
      "   loading parameter \"an14\"...\n",
      "   loading parameter \"an15\"...\n",
      "   loading parameter \"an16\"...\n",
      "   loading parameter \"an17\"...\n",
      "   loading hk layer   1...\n",
      "   loading hani layer   1...\n",
      "   loading vka layer   1...\n",
      "   loading ss layer   1...\n",
      "   loading sy layer   1...\n",
      "   loading hk layer   2...\n",
      "   loading hani layer   2...\n",
      "   loading vka layer   2...\n",
      "   loading ss layer   2...\n",
      "   loading sy layer   2...\n",
      "   loading hk layer   3...\n",
      "   loading hani layer   3...\n",
      "   loading vka layer   3...\n",
      "   loading ss layer   3...\n",
      "   loading sy layer   3...\n",
      "   loading hk layer   4...\n",
      "   loading hani layer   4...\n",
      "   loading vka layer   4...\n",
      "   loading ss layer   4...\n",
      "   loading sy layer   4...\n",
      "   loading hk layer   5...\n",
      "   loading hani layer   5...\n",
      "   loading vka layer   5...\n",
      "   loading ss layer   5...\n",
      "   loading sy layer   5...\n",
      "   loading hk layer   6...\n",
      "   loading hani layer   6...\n",
      "   loading vka layer   6...\n",
      "   loading ss layer   6...\n",
      "   loading sy layer   6...\n",
      "   loading hk layer   7...\n",
      "   loading hani layer   7...\n",
      "   loading vka layer   7...\n",
      "   loading ss layer   7...\n",
      "   loading sy layer   7...\n",
      "   loading hk layer   8...\n",
      "   loading hani layer   8...\n",
      "   loading vka layer   8...\n",
      "   loading ss layer   8...\n",
      "   loading sy layer   8...\n",
      "   loading hk layer   9...\n",
      "   loading hani layer   9...\n",
      "   loading vka layer   9...\n",
      "   loading ss layer   9...\n",
      "   loading sy layer   9...\n",
      "   loading hk layer  10...\n",
      "   loading hani layer  10...\n",
      "   loading vka layer  10...\n",
      "   loading ss layer  10...\n",
      "   loading sy layer  10...\n",
      "   loading hk layer  11...\n",
      "   loading hani layer  11...\n",
      "   loading vka layer  11...\n",
      "   loading ss layer  11...\n",
      "   loading sy layer  11...\n",
      "   loading hk layer  12...\n",
      "   loading hani layer  12...\n",
      "   loading vka layer  12...\n",
      "   loading ss layer  12...\n",
      "   loading sy layer  12...\n",
      "   loading hk layer  13...\n",
      "   loading hani layer  13...\n",
      "   loading vka layer  13...\n",
      "   loading ss layer  13...\n",
      "   loading sy layer  13...\n",
      "   loading hk layer  14...\n",
      "   loading hani layer  14...\n",
      "   loading vka layer  14...\n",
      "   loading ss layer  14...\n",
      "   loading sy layer  14...\n",
      "   loading hk layer  15...\n",
      "   loading hani layer  15...\n",
      "   loading vka layer  15...\n",
      "   loading ss layer  15...\n",
      "   loading sy layer  15...\n",
      "   loading hk layer  16...\n",
      "   loading hani layer  16...\n",
      "   loading vka layer  16...\n",
      "   loading ss layer  16...\n",
      "   loading sy layer  16...\n",
      "   loading hk layer  17...\n",
      "   loading hani layer  17...\n",
      "   loading vka layer  17...\n",
      "   loading ss layer  17...\n",
      "   loading sy layer  17...\n",
      "adding Package:  LPF\n",
      "   LPF  package load...success\n",
      "\n",
      "   The following 3 packages were successfully loaded.\n",
      "      WB19BS1E0003_PT.dis\n",
      "      WB19BS1E0003_PT.bas6\n",
      "      WB19BS1E0003_PT.lpf\n",
      "   The following 1 packages were not loaded.\n",
      "      WB19BS1E0003_PT.list\n"
     ]
    }
   ],
   "source": [
    "# Load an MF2005 model\n",
    "ml = flopy.modflow.Modflow.load(\n",
    "    f\"{new_nam}.nam\",\n",
    "    model_ws=model_f,\n",
    "    verbose=True,\n",
    "    version=\"mf2005\",\n",
    "    check=False,\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "54a2d342-6b88-436a-a3b1-f8932cbac163",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "17"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dis = ml.get_package(\"dis\")\n",
    "dis.nlay"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1a8e1170-d32c-4223-9234-dd36a6430f9a",
   "metadata": {},
   "source": [
    "# Do a vtk export for viewing with paraview"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "61e1bf77-aa3e-4ee7-b4f5-02d424fa6ecc",
   "metadata": {},
   "outputs": [],
   "source": [
    "ml.export(plots_f, fmt=\"vtk\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e7099664-713a-413a-abf2-657016aff6d8",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}