examples/Notebooks/flopy3_MT3D-USGS_example.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# FloPy\n",
    "\n",
    "### MT3D-USGS Example\n",
    "\n",
    "Demonstrates functionality of the flopy MT3D-USGS module using the 'Crank-Nicolson' example distributed with MT3D-USGS.  \n",
    "\n",
    "#### Problem description:\n",
    "\n",
    "* Grid dimensions: 1 Layer, 3 Rows, 650 Columns\n",
    "* Stress periods: 3\n",
    "* Units are in seconds and meters\n",
    "* Flow package: UPW\n",
    "* Stress packages: SFR, GHB\n",
    "* Solvers: NWT, GCG"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3.8.11 (default, Aug  6 2021, 08:56:27) \n",
      "[Clang 10.0.0 ]\n",
      "numpy version: 1.19.2\n",
      "matplotlib version: 3.4.2\n",
      "flopy version: 3.3.5\n"
     ]
    }
   ],
   "source": [
    "import sys\n",
    "import os\n",
    "import platform\n",
    "import string\n",
    "from io import StringIO, BytesIO\n",
    "import math\n",
    "import numpy as np\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# run installed version of flopy or add local path\n",
    "try:\n",
    "    import flopy\n",
    "except:\n",
    "    fpth = os.path.abspath(os.path.join(\"..\", \"..\"))\n",
    "    sys.path.append(fpth)\n",
    "    import flopy\n",
    "\n",
    "print(sys.version)\n",
    "print(\"numpy version: {}\".format(np.__version__))\n",
    "print(\"matplotlib version: {}\".format(mpl.__version__))\n",
    "print(\"flopy version: {}\".format(flopy.__version__))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "modelpth = os.path.join(\"data\")\n",
    "modelname = \"CrnkNic\"\n",
    "mfexe = \"mfnwt\"\n",
    "mtexe = \"mt3dusgs\"\n",
    "if platform.system() == \"Windows\":\n",
    "    mfexe += \".exe\"\n",
    "    mtexe += \".exe\"\n",
    "\n",
    "# Make sure modelpth directory exists\n",
    "if not os.path.isdir(modelpth):\n",
    "    os.makedirs(modelpth, exist_ok=True)\n",
    "\n",
    "# Instantiate MODFLOW object in flopy\n",
    "mf = flopy.modflow.Modflow(\n",
    "    modelname=modelname, exe_name=mfexe, model_ws=modelpth, version=\"mfnwt\"\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Set up model discretization"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "Lx = 650.0\n",
    "Ly = 15\n",
    "nrow = 3\n",
    "ncol = 650\n",
    "nlay = 1\n",
    "\n",
    "delr = Lx / ncol\n",
    "delc = Ly / nrow\n",
    "\n",
    "xmax = ncol * delr\n",
    "ymax = nrow * delc\n",
    "\n",
    "X, Y = np.meshgrid(\n",
    "    np.linspace(delr / 2, xmax - delr / 2, ncol),\n",
    "    np.linspace(ymax - delc / 2, 0 + delc / 2, nrow),\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate output control (oc) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Output Control: Create a flopy output control object\n",
    "oc = flopy.modflow.ModflowOc(mf)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate solver package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Newton-Raphson Solver: Create a flopy nwt package object\n",
    "\n",
    "headtol = 1.0e-4\n",
    "fluxtol = 5\n",
    "maxiterout = 5000\n",
    "thickfact = 1e-06\n",
    "linmeth = 2\n",
    "iprnwt = 1\n",
    "ibotav = 1\n",
    "\n",
    "nwt = flopy.modflow.ModflowNwt(\n",
    "    mf,\n",
    "    headtol=headtol,\n",
    "    fluxtol=fluxtol,\n",
    "    maxiterout=maxiterout,\n",
    "    thickfact=thickfact,\n",
    "    linmeth=linmeth,\n",
    "    iprnwt=iprnwt,\n",
    "    ibotav=ibotav,\n",
    "    options=\"SIMPLE\",\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate discretization (DIS) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# The equations for calculating the ground elevation in the 1 Layer CrnkNic model.\n",
    "# Although Y isn't used, keeping it here for symetry\n",
    "\n",
    "\n",
    "def topElev(X, Y):\n",
    "    return 100.0 - (np.ceil(X) - 1) * 0.03\n",
    "\n",
    "\n",
    "grndElev = topElev(X, Y)\n",
    "bedRockElev = grndElev - 3.0\n",
    "\n",
    "Steady = [False, False, False]\n",
    "nstp = [1, 1, 1]\n",
    "tsmult = [1.0, 1.0, 1.0]\n",
    "\n",
    "# Stress periods extend from (12AM-8:29:59AM); (8:30AM-11:30:59AM); (11:31AM-23:59:59PM)\n",
    "perlen = [30600, 10800, 45000]\n",
    "\n",
    "# Create the discretization object\n",
    "# itmuni = 1 (seconds); lenuni = 2 (meters)\n",
    "dis = flopy.modflow.ModflowDis(\n",
    "    mf,\n",
    "    nlay,\n",
    "    nrow,\n",
    "    ncol,\n",
    "    nper=3,\n",
    "    delr=delr,\n",
    "    delc=delc,\n",
    "    top=grndElev,\n",
    "    botm=bedRockElev,\n",
    "    laycbd=0,\n",
    "    itmuni=1,\n",
    "    lenuni=2,\n",
    "    steady=Steady,\n",
    "    nstp=nstp,\n",
    "    tsmult=tsmult,\n",
    "    perlen=perlen,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate upstream weighting (UPW) flow package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# UPW parameters\n",
    "# UPW must be instantiated after DIS.  Otherwise, during the mf.write_input() procedures,\n",
    "# flopy will crash.\n",
    "\n",
    "\n",
    "laytyp = 1\n",
    "layavg = 2\n",
    "chani = 1.0\n",
    "layvka = 1\n",
    "iphdry = 0\n",
    "hk = 0.1\n",
    "hani = 1\n",
    "vka = 1.0\n",
    "ss = 0.000001\n",
    "sy = 0.20\n",
    "hdry = -888\n",
    "\n",
    "upw = flopy.modflow.ModflowUpw(\n",
    "    mf,\n",
    "    laytyp=laytyp,\n",
    "    layavg=layavg,\n",
    "    chani=chani,\n",
    "    layvka=layvka,\n",
    "    ipakcb=53,\n",
    "    hdry=hdry,\n",
    "    iphdry=iphdry,\n",
    "    hk=hk,\n",
    "    hani=hani,\n",
    "    vka=vka,\n",
    "    ss=ss,\n",
    "    sy=sy,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate basic (BAS or BA6) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a flopy basic package object\n",
    "def calc_strtElev(X, Y):\n",
    "    return 99.5 - (np.ceil(X) - 1) * 0.0001\n",
    "\n",
    "\n",
    "ibound = np.ones((nlay, nrow, ncol))\n",
    "ibound[:, 0, :] *= -1\n",
    "ibound[:, 2, :] *= -1\n",
    "\n",
    "strtElev = calc_strtElev(X, Y)\n",
    "\n",
    "bas = flopy.modflow.ModflowBas(mf, ibound=ibound, hnoflo=hdry, strt=strtElev)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate streamflow routing (SFR2) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Streamflow Routing Package: Try and set up with minimal options in use\n",
    "# 9 11 IFACE # Data Set 1: ISTCB1  ISTCB2\n",
    "\n",
    "nstrm = ncol\n",
    "nss = 6\n",
    "const = 1.0\n",
    "dleak = 0.0001\n",
    "istcb1 = -10\n",
    "istcb2 = 11\n",
    "isfropt = 1\n",
    "\n",
    "segment_data = None\n",
    "channel_geometry_data = None\n",
    "channel_flow_data = None\n",
    "dataset_5 = None\n",
    "reachinput = True\n",
    "\n",
    "\n",
    "# The next couple of lines set up the reach_data for the 30x100 hypothetical model.\n",
    "# Will need to adjust the row based on which grid discretization we're doing.\n",
    "# Ensure that the stream goes down one of the middle rows of the model.\n",
    "\n",
    "strmBed_Elev = 98.75 - (np.ceil(X[1, :]) - 1) * 0.0001\n",
    "\n",
    "s1 = \"k,i,j,iseg,ireach,rchlen,strtop,slope,strthick,strhc1\\n\"\n",
    "iseg = 0\n",
    "irch = 0\n",
    "\n",
    "for y in range(ncol):\n",
    "    if y <= 37:\n",
    "        if iseg == 0:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 1\n",
    "        strhc1 = 1.0e-10\n",
    "    elif y <= 104:\n",
    "        if iseg == 1:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 2\n",
    "        strhc1 = 1.0e-10\n",
    "    elif y <= 280:\n",
    "        if iseg == 2:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 3\n",
    "        strhc1 = 2.946219199e-6\n",
    "    elif y <= 432:\n",
    "        if iseg == 3:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 4\n",
    "        strhc1 = 1.375079882e-6\n",
    "    elif y <= 618:\n",
    "        if iseg == 4:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 5\n",
    "        strhc1 = 1.764700062e-6\n",
    "    else:\n",
    "        if iseg == 5:\n",
    "            irch = 1\n",
    "        else:\n",
    "            irch += 1\n",
    "        iseg = 6\n",
    "        strhc1 = 1e-10\n",
    "\n",
    "    # remember that lay, row, col need to be zero-based and are adjusted accordingly by flopy\n",
    "    #    layer +  row      +      col     +       iseg      +      irch      +     rchlen      +            strtop          +       slope       +     strthick    +     strmbed K\n",
    "    s1 += \"0,{}\".format(1)\n",
    "    s1 += \",{}\".format(y)\n",
    "    s1 += \",{}\".format(iseg)\n",
    "    s1 += \",{}\".format(irch)\n",
    "    s1 += \",{}\".format(delr)\n",
    "    s1 += \",{}\".format(strmBed_Elev[y])\n",
    "    s1 += \",{}\".format(0.0001)\n",
    "    s1 += \",{}\".format(0.50)\n",
    "    s1 += \",{}\\n\".format(strhc1)\n",
    "\n",
    "\n",
    "if not os.path.exists(\"data\"):\n",
    "    os.makedirs(\"data\", exist_ok=True)\n",
    "\n",
    "fpth = os.path.join(\"data\", \"s1.csv\")\n",
    "f = open(fpth, \"w\")\n",
    "f.write(s1)\n",
    "f.close()\n",
    "\n",
    "dtype = [\n",
    "    (\"k\", \"<i4\"),\n",
    "    (\"i\", \"<i4\"),\n",
    "    (\"j\", \"<i4\"),\n",
    "    (\"iseg\", \"<i4\"),\n",
    "    (\"ireach\", \"<f8\"),\n",
    "    (\"rchlen\", \"<f8\"),\n",
    "    (\"strtop\", \"<f8\"),\n",
    "    (\"slope\", \"<f8\"),\n",
    "    (\"strthick\", \"<f8\"),\n",
    "    (\"strhc1\", \"<f8\"),\n",
    "]\n",
    "\n",
    "if sys.version_info > (3, 0):\n",
    "    f = open(fpth, \"rb\")\n",
    "else:\n",
    "    f = open(fpth, \"r\")\n",
    "\n",
    "reach_data = np.genfromtxt(f, delimiter=\",\", names=True, dtype=dtype)\n",
    "f.close()\n",
    "\n",
    "s2 = \"nseg,icalc,outseg,iupseg,nstrpts,   flow,runoff,etsw,pptsw,        roughch,        roughbk,cdpth,fdpth,awdth,bwdth,width1,width2\\n \\\n",
    "      1,    1,     2,     0,      0, 0.0125,   0.0, 0.0,  0.0, 0.082078856000, 0.082078856000,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\\n \\\n",
    "      2,    1,     3,     0,      0,    0.0,   0.0, 0.0,  0.0, 0.143806300000, 0.143806300000,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\\n \\\n",
    "      3,    1,     4,     0,      0,    0.0,   0.0, 0.0,  0.0, 0.104569661821, 0.104569661821,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\\n \\\n",
    "      4,    1,     5,     0,      0,    0.0,   0.0, 0.0,  0.0, 0.126990045841, 0.126990045841,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\\n \\\n",
    "      5,    1,     6,     0,      0,    0.0,   0.0, 0.0,  0.0, 0.183322283828, 0.183322283828,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\\n \\\n",
    "      6,    1,     0,     0,      0,    0.0,   0.0, 0.0,  0.0, 0.183322283828, 0.183322283828,  0.0,  0.0,  0.0,  0.0,   1.5,   1.5\"\n",
    "\n",
    "fpth = os.path.join(\"data\", \"s2.csv\")\n",
    "f = open(fpth, \"w\")\n",
    "f.write(s2)\n",
    "f.close()\n",
    "\n",
    "if sys.version_info > (3, 0):\n",
    "    f = open(fpth, \"rb\")\n",
    "else:\n",
    "    f = open(fpth, \"r\")\n",
    "\n",
    "segment_data = np.genfromtxt(f, delimiter=\",\", names=True)\n",
    "f.close()\n",
    "\n",
    "# Be sure to convert segment_data to a dictionary keyed on stress period.\n",
    "segment_data = np.atleast_1d(segment_data)\n",
    "segment_data = {0: segment_data, 1: segment_data, 2: segment_data}\n",
    "\n",
    "# There are 3 stress periods\n",
    "dataset_5 = {0: [nss, 0, 0], 1: [nss, 0, 0], 2: [nss, 0, 0]}\n",
    "\n",
    "sfr = flopy.modflow.ModflowSfr2(\n",
    "    mf,\n",
    "    nstrm=nstrm,\n",
    "    nss=nss,\n",
    "    const=const,\n",
    "    dleak=dleak,\n",
    "    isfropt=isfropt,\n",
    "    istcb2=0,\n",
    "    reachinput=True,\n",
    "    reach_data=reach_data,\n",
    "    dataset_5=dataset_5,\n",
    "    segment_data=segment_data,\n",
    "    channel_geometry_data=channel_geometry_data,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate gage package for use with MODFLOW-NWT package"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "gages = [\n",
    "    [1, 38, 61, 1],\n",
    "    [2, 67, 62, 1],\n",
    "    [3, 176, 63, 1],\n",
    "    [4, 152, 64, 1],\n",
    "    [5, 186, 65, 1],\n",
    "    [6, 31, 66, 1],\n",
    "]\n",
    "files = [\n",
    "    \"CrnkNic.gage\",\n",
    "    \"CrnkNic.gag1\",\n",
    "    \"CrnkNic.gag2\",\n",
    "    \"CrnkNic.gag3\",\n",
    "    \"CrnkNic.gag4\",\n",
    "    \"CrnkNic.gag5\",\n",
    "    \"CrnkNic.gag6\",\n",
    "]\n",
    "gage = flopy.modflow.ModflowGage(\n",
    "    mf, numgage=6, gage_data=gages, filenames=files\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate linkage with mass transport routing (LMT) package for MODFLOW-NWT (generates linker file)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "lmt = flopy.modflow.ModflowLmt(\n",
    "    mf,\n",
    "    output_file_name=\"CrnkNic.ftl\",\n",
    "    output_file_header=\"extended\",\n",
    "    output_file_format=\"formatted\",\n",
    "    package_flows=[\"sfr\"],\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Write the MODFLOW input files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/Users/jdhughes/Documents/Development/flopy_git/flopy_fork/examples/Notebooks\n",
      "FloPy is using the following executable to run the model: /Users/jdhughes/.local/bin/mfnwt\n",
      "\n",
      "                                  MODFLOW-NWT-SWR1 \n",
      "    U.S. GEOLOGICAL SURVEY MODULAR FINITE-DIFFERENCE GROUNDWATER-FLOW MODEL\n",
      "                             WITH NEWTON FORMULATION\n",
      "                             Version 1.2.0 03/01/2020                        \n",
      "                    BASED ON MODFLOW-2005 Version 1.12.0 02/03/2017                       \n",
      "\n",
      "                    SWR1 Version 1.04.0 09/15/2016                       \n",
      "\n",
      " Using NAME file: CrnkNic.nam \n",
      " Run start date and time (yyyy/mm/dd hh:mm:ss): 2021/10/25 16:36:03\n",
      "\n",
      " Solving:  Stress period:     1    Time step:     1    Groundwater-Flow Eqn.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Solving:  Stress period:     2    Time step:     1    Groundwater-Flow Eqn.\n",
      " Solving:  Stress period:     3    Time step:     1    Groundwater-Flow Eqn.\n",
      " Run end date and time (yyyy/mm/dd hh:mm:ss): 2021/10/25 16:36:03\n",
      " Elapsed run time:  0.119 Seconds\n",
      "\n",
      "  Normal termination of simulation\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(True, [])"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pth = os.getcwd()\n",
    "print(pth)\n",
    "\n",
    "mf.write_input()\n",
    "\n",
    "# run the model\n",
    "mf.run_model()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now draft up MT3D-USGS input files.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate MT3D-USGS object in flopy\n",
    "mt = flopy.mt3d.Mt3dms(\n",
    "    modflowmodel=mf,\n",
    "    modelname=modelname,\n",
    "    model_ws=modelpth,\n",
    "    version=\"mt3d-usgs\",\n",
    "    namefile_ext=\"mtnam\",\n",
    "    exe_name=mtexe,\n",
    "    ftlfilename=\"CrnkNic.ftl\",\n",
    "    ftlfree=True,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate basic transport (BTN) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "btn = flopy.mt3d.Mt3dBtn(\n",
    "    mt,\n",
    "    sconc=3.7,\n",
    "    ncomp=1,\n",
    "    prsity=0.2,\n",
    "    cinact=-1.0,\n",
    "    thkmin=0.001,\n",
    "    nprs=-1,\n",
    "    nprobs=10,\n",
    "    chkmas=True,\n",
    "    nprmas=10,\n",
    "    dt0=180,\n",
    "    mxstrn=2500,\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate advection (ADV) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "adv = flopy.mt3d.Mt3dAdv(mt, mixelm=0, percel=1.00, mxpart=5000, nadvfd=1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instatiate generalized conjugate gradient solver (GCG) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "rct = flopy.mt3d.Mt3dRct(mt, isothm=0, ireact=100, igetsc=0, rc1=0.01)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "gcg = flopy.mt3d.Mt3dGcg(\n",
    "    mt, mxiter=10, iter1=50, isolve=3, ncrs=0, accl=1, cclose=1e-6, iprgcg=1\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate source-sink mixing (SSM) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "# For SSM, need to set the constant head boundary conditions to the ambient concentration\n",
    "# for all 1,300 constant head boundary cells.\n",
    "\n",
    "itype = flopy.mt3d.Mt3dSsm.itype_dict()\n",
    "ssm_data = {}\n",
    "ssm_data[0] = [(0, 0, 0, 3.7, itype[\"CHD\"])]\n",
    "ssm_data[0].append((0, 2, 0, 3.7, itype[\"CHD\"]))\n",
    "for i in [0, 2]:\n",
    "    for j in range(1, ncol):\n",
    "        ssm_data[0].append((0, i, j, 3.7, itype[\"CHD\"]))\n",
    "\n",
    "ssm = flopy.mt3d.Mt3dSsm(mt, stress_period_data=ssm_data)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate streamflow transport (SFT) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "dispsf = []\n",
    "for y in range(ncol):\n",
    "    if y <= 37:\n",
    "        dispsf.append(0.12)\n",
    "    elif y <= 104:\n",
    "        dispsf.append(0.15)\n",
    "    elif y <= 280:\n",
    "        dispsf.append(0.24)\n",
    "    elif y <= 432:\n",
    "        dispsf.append(0.31)\n",
    "    elif y <= 618:\n",
    "        dispsf.append(0.40)\n",
    "    else:\n",
    "        dispsf.append(0.40)\n",
    "\n",
    "# Enter a list of the observation points\n",
    "# Each observation is taken as the last reach within the first 5 segments\n",
    "\n",
    "seg_len = np.unique(reach_data[\"iseg\"], return_counts=True)\n",
    "obs_sf = np.cumsum(seg_len[1])\n",
    "obs_sf = obs_sf.tolist()\n",
    "\n",
    "# The last reach is not an observation point, therefore drop\n",
    "obs_sf.pop(-1)\n",
    "\n",
    "# In the first and last stress periods, concentration at the headwater is 3.7\n",
    "sf_stress_period_data = {0: [0, 0, 3.7], 1: [0, 0, 11.4], 2: [0, 0, 3.7]}\n",
    "\n",
    "gage_output = [None, None, \"CrnkNic.sftobs\"]\n",
    "\n",
    "sft = flopy.mt3d.Mt3dSft(\n",
    "    mt,\n",
    "    nsfinit=650,\n",
    "    mxsfbc=650,\n",
    "    icbcsf=81,\n",
    "    ioutobs=82,\n",
    "    isfsolv=1,\n",
    "    cclosesf=1.0e-6,\n",
    "    mxitersf=10,\n",
    "    crntsf=1.0,\n",
    "    iprtxmd=0,\n",
    "    coldsf=3.7,\n",
    "    dispsf=dispsf,\n",
    "    nobssf=5,\n",
    "    obs_sf=obs_sf,\n",
    "    sf_stress_period_data=sf_stress_period_data,\n",
    "    filenames=gage_output,\n",
    ")\n",
    "\n",
    "sft.dispsf[0].format.fortran = \"(10E15.6)\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Write the MT3D-USGS input files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "FloPy is using the following executable to run the model: /Users/jdhughes/.local/bin/mt3dusgs\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " MT3D-USGS - Modular 3D Multi-Species Transport Model [Ver 1.1.0]   \n",
      " and based on MT3DMS. MT3D-USGS developed in cooperation by \n",
      " S.S. Papadopulos & Associates and the U.S. Geological Survey\n",
      "\n",
      " Using NAME File: CrnkNic.mtnam\n",
      "\n",
      " STRESS PERIOD NO.    1\n",
      "\n",
      " TIME STEP NO.    1\n",
      " FROM TIME =   0.0000     TO    30600.    \n",
      "\n",
      " Transport Step:    1   Step Size:   180.0     Total Elapsed Time:   180.00    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =  0.4865      [K,I,J]    1    2  289\n",
      " Outer Iter.  1  Inner Iter.  2:  Max. DC =  0.5960E-07  [K,I,J]    1    1   36\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    2   Step Size:   180.0     Total Elapsed Time:   360.00    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =  0.9385      [K,I,J]    1    2  394\n",
      " Outer Iter.  1  Inner Iter.  2:  Max. DC =  0.1788E-06  [K,I,J]    1    2  294\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =  0.3934E-05  [K,I,J]    1    2  294\n",
      " Outer Iter.  2  Inner Iter.  2:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  3  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    3   Step Size:   180.0     Total Elapsed Time:   540.00    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =  0.7602      [K,I,J]    1    1    1\n",
      " Outer Iter.  1  Inner Iter.  2:  Max. DC =  0.1788E-06  [K,I,J]    1    2    2\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =  0.8884E-06  [K,I,J]    1    1  642\n",
      " Transport Step:    4   Step Size:   180.0     Total Elapsed Time:   720.00    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =  0.3366E-03  [K,I,J]    1    2    1\n",
      " Outer Iter.  1  Inner Iter.  2:  Max. DC =  0.4657E-09  [K,I,J]    1    2    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    5   Step Size:   180.0     Total Elapsed Time:   900.00    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    6   Step Size:   180.0     Total Elapsed Time:   1080.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    7   Step Size:   180.0     Total Elapsed Time:   1260.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    8   Step Size:   180.0     Total Elapsed Time:   1440.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    9   Step Size:   180.0     Total Elapsed Time:   1620.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   10   Step Size:   180.0     Total Elapsed Time:   1800.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   11   Step Size:   180.0     Total Elapsed Time:   1980.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   12   Step Size:   180.0     Total Elapsed Time:   2160.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   13   Step Size:   180.0     Total Elapsed Time:   2340.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   14   Step Size:   180.0     Total Elapsed Time:   2520.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   15   Step Size:   180.0     Total Elapsed Time:   2700.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   16   Step Size:   180.0     Total Elapsed Time:   2880.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   17   Step Size:   180.0     Total Elapsed Time:   3060.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   18   Step Size:   180.0     Total Elapsed Time:   3240.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   19   Step Size:   180.0     Total Elapsed Time:   3420.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   20   Step Size:   180.0     Total Elapsed Time:   3600.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   21   Step Size:   180.0     Total Elapsed Time:   3780.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   22   Step Size:   180.0     Total Elapsed Time:   3960.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   23   Step Size:   180.0     Total Elapsed Time:   4140.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   24   Step Size:   180.0     Total Elapsed Time:   4320.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   25   Step Size:   180.0     Total Elapsed Time:   4500.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   26   Step Size:   180.0     Total Elapsed Time:   4680.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   27   Step Size:   180.0     Total Elapsed Time:   4860.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   28   Step Size:   180.0     Total Elapsed Time:   5040.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   29   Step Size:   180.0     Total Elapsed Time:   5220.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   30   Step Size:   180.0     Total Elapsed Time:   5400.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   31   Step Size:   180.0     Total Elapsed Time:   5580.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   32   Step Size:   180.0     Total Elapsed Time:   5760.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   33   Step Size:   180.0     Total Elapsed Time:   5940.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   34   Step Size:   180.0     Total Elapsed Time:   6120.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   35   Step Size:   180.0     Total Elapsed Time:   6300.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   36   Step Size:   180.0     Total Elapsed Time:   6480.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   37   Step Size:   180.0     Total Elapsed Time:   6660.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   38   Step Size:   180.0     Total Elapsed Time:   6840.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   39   Step Size:   180.0     Total Elapsed Time:   7020.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   40   Step Size:   180.0     Total Elapsed Time:   7200.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   41   Step Size:   180.0     Total Elapsed Time:   7380.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   42   Step Size:   180.0     Total Elapsed Time:   7560.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   43   Step Size:   180.0     Total Elapsed Time:   7740.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   44   Step Size:   180.0     Total Elapsed Time:   7920.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   45   Step Size:   180.0     Total Elapsed Time:   8100.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   46   Step Size:   180.0     Total Elapsed Time:   8280.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   47   Step Size:   180.0     Total Elapsed Time:   8460.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   48   Step Size:   180.0     Total Elapsed Time:   8640.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   49   Step Size:   180.0     Total Elapsed Time:   8820.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   50   Step Size:   180.0     Total Elapsed Time:   9000.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   51   Step Size:   180.0     Total Elapsed Time:   9180.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   52   Step Size:   180.0     Total Elapsed Time:   9360.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   53   Step Size:   180.0     Total Elapsed Time:   9540.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   54   Step Size:   180.0     Total Elapsed Time:   9720.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   55   Step Size:   180.0     Total Elapsed Time:   9900.0    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   56   Step Size:   180.0     Total Elapsed Time:   10080.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   57   Step Size:   180.0     Total Elapsed Time:   10260.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   58   Step Size:   180.0     Total Elapsed Time:   10440.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   59   Step Size:   180.0     Total Elapsed Time:   10620.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   60   Step Size:   180.0     Total Elapsed Time:   10800.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   61   Step Size:   180.0     Total Elapsed Time:   10980.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   62   Step Size:   180.0     Total Elapsed Time:   11160.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Transport Step:   63   Step Size:   180.0     Total Elapsed Time:   11340.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   64   Step Size:   180.0     Total Elapsed Time:   11520.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   65   Step Size:   180.0     Total Elapsed Time:   11700.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   66   Step Size:   180.0     Total Elapsed Time:   11880.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   67   Step Size:   180.0     Total Elapsed Time:   12060.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   68   Step Size:   180.0     Total Elapsed Time:   12240.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   69   Step Size:   180.0     Total Elapsed Time:   12420.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   70   Step Size:   180.0     Total Elapsed Time:   12600.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   71   Step Size:   180.0     Total Elapsed Time:   12780.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   72   Step Size:   180.0     Total Elapsed Time:   12960.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   73   Step Size:   180.0     Total Elapsed Time:   13140.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   74   Step Size:   180.0     Total Elapsed Time:   13320.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   75   Step Size:   180.0     Total Elapsed Time:   13500.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   76   Step Size:   180.0     Total Elapsed Time:   13680.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   77   Step Size:   180.0     Total Elapsed Time:   13860.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   78   Step Size:   180.0     Total Elapsed Time:   14040.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   79   Step Size:   180.0     Total Elapsed Time:   14220.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   80   Step Size:   180.0     Total Elapsed Time:   14400.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   81   Step Size:   180.0     Total Elapsed Time:   14580.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   82   Step Size:   180.0     Total Elapsed Time:   14760.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   83   Step Size:   180.0     Total Elapsed Time:   14940.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   84   Step Size:   180.0     Total Elapsed Time:   15120.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   85   Step Size:   180.0     Total Elapsed Time:   15300.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   86   Step Size:   180.0     Total Elapsed Time:   15480.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   87   Step Size:   180.0     Total Elapsed Time:   15660.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   88   Step Size:   180.0     Total Elapsed Time:   15840.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   89   Step Size:   180.0     Total Elapsed Time:   16020.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   90   Step Size:   180.0     Total Elapsed Time:   16200.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   91   Step Size:   180.0     Total Elapsed Time:   16380.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   92   Step Size:   180.0     Total Elapsed Time:   16560.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   93   Step Size:   180.0     Total Elapsed Time:   16740.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   94   Step Size:   180.0     Total Elapsed Time:   16920.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   95   Step Size:   180.0     Total Elapsed Time:   17100.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   96   Step Size:   180.0     Total Elapsed Time:   17280.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   97   Step Size:   180.0     Total Elapsed Time:   17460.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   98   Step Size:   180.0     Total Elapsed Time:   17640.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   99   Step Size:   180.0     Total Elapsed Time:   17820.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  100   Step Size:   180.0     Total Elapsed Time:   18000.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  101   Step Size:   180.0     Total Elapsed Time:   18180.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  102   Step Size:   180.0     Total Elapsed Time:   18360.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  103   Step Size:   180.0     Total Elapsed Time:   18540.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  104   Step Size:   180.0     Total Elapsed Time:   18720.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  105   Step Size:   180.0     Total Elapsed Time:   18900.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  106   Step Size:   180.0     Total Elapsed Time:   19080.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  107   Step Size:   180.0     Total Elapsed Time:   19260.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  108   Step Size:   180.0     Total Elapsed Time:   19440.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  109   Step Size:   180.0     Total Elapsed Time:   19620.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  110   Step Size:   180.0     Total Elapsed Time:   19800.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  111   Step Size:   180.0     Total Elapsed Time:   19980.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  112   Step Size:   180.0     Total Elapsed Time:   20160.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  113   Step Size:   180.0     Total Elapsed Time:   20340.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  114   Step Size:   180.0     Total Elapsed Time:   20520.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  115   Step Size:   180.0     Total Elapsed Time:   20700.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  116   Step Size:   180.0     Total Elapsed Time:   20880.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  117   Step Size:   180.0     Total Elapsed Time:   21060.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  118   Step Size:   180.0     Total Elapsed Time:   21240.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  119   Step Size:   180.0     Total Elapsed Time:   21420.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  120   Step Size:   180.0     Total Elapsed Time:   21600.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  121   Step Size:   180.0     Total Elapsed Time:   21780.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  122   Step Size:   180.0     Total Elapsed Time:   21960.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Transport Step:  123   Step Size:   180.0     Total Elapsed Time:   22140.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  124   Step Size:   180.0     Total Elapsed Time:   22320.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  125   Step Size:   180.0     Total Elapsed Time:   22500.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  126   Step Size:   180.0     Total Elapsed Time:   22680.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  127   Step Size:   180.0     Total Elapsed Time:   22860.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  128   Step Size:   180.0     Total Elapsed Time:   23040.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  129   Step Size:   180.0     Total Elapsed Time:   23220.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  130   Step Size:   180.0     Total Elapsed Time:   23400.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  131   Step Size:   180.0     Total Elapsed Time:   23580.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  132   Step Size:   180.0     Total Elapsed Time:   23760.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  133   Step Size:   180.0     Total Elapsed Time:   23940.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  134   Step Size:   180.0     Total Elapsed Time:   24120.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  135   Step Size:   180.0     Total Elapsed Time:   24300.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  136   Step Size:   180.0     Total Elapsed Time:   24480.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  137   Step Size:   180.0     Total Elapsed Time:   24660.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  138   Step Size:   180.0     Total Elapsed Time:   24840.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  139   Step Size:   180.0     Total Elapsed Time:   25020.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  140   Step Size:   180.0     Total Elapsed Time:   25200.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  141   Step Size:   180.0     Total Elapsed Time:   25380.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  142   Step Size:   180.0     Total Elapsed Time:   25560.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  143   Step Size:   180.0     Total Elapsed Time:   25740.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  144   Step Size:   180.0     Total Elapsed Time:   25920.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  145   Step Size:   180.0     Total Elapsed Time:   26100.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  146   Step Size:   180.0     Total Elapsed Time:   26280.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  147   Step Size:   180.0     Total Elapsed Time:   26460.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  148   Step Size:   180.0     Total Elapsed Time:   26640.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  149   Step Size:   180.0     Total Elapsed Time:   26820.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  150   Step Size:   180.0     Total Elapsed Time:   27000.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  151   Step Size:   180.0     Total Elapsed Time:   27180.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  152   Step Size:   180.0     Total Elapsed Time:   27360.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  153   Step Size:   180.0     Total Elapsed Time:   27540.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  154   Step Size:   180.0     Total Elapsed Time:   27720.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  155   Step Size:   180.0     Total Elapsed Time:   27900.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  156   Step Size:   180.0     Total Elapsed Time:   28080.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  157   Step Size:   180.0     Total Elapsed Time:   28260.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  158   Step Size:   180.0     Total Elapsed Time:   28440.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  159   Step Size:   180.0     Total Elapsed Time:   28620.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  160   Step Size:   180.0     Total Elapsed Time:   28800.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  161   Step Size:   180.0     Total Elapsed Time:   28980.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  162   Step Size:   180.0     Total Elapsed Time:   29160.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  163   Step Size:   180.0     Total Elapsed Time:   29340.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  164   Step Size:   180.0     Total Elapsed Time:   29520.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  165   Step Size:   180.0     Total Elapsed Time:   29700.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  166   Step Size:   180.0     Total Elapsed Time:   29880.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  167   Step Size:   180.0     Total Elapsed Time:   30060.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  168   Step Size:   180.0     Total Elapsed Time:   30240.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  169   Step Size:   180.0     Total Elapsed Time:   30420.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:  170   Step Size:   180.0     Total Elapsed Time:   30600.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      "\n",
      " STRESS PERIOD NO.    2\n",
      "\n",
      " TIME STEP NO.    1\n",
      " FROM TIME =   30600.     TO    41400.    \n",
      "\n",
      " Transport Step:    1   Step Size:   180.0     Total Elapsed Time:   30780.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =  0.2062E-04  [K,I,J]    1    1    1\n",
      " Outer Iter.  1  Inner Iter.  2:  Max. DC =  0.6548E-10  [K,I,J]    1    2    3\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    2   Step Size:   180.0     Total Elapsed Time:   30960.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    3   Step Size:   180.0     Total Elapsed Time:   31140.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    4   Step Size:   180.0     Total Elapsed Time:   31320.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    5   Step Size:   180.0     Total Elapsed Time:   31500.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    6   Step Size:   180.0     Total Elapsed Time:   31680.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    7   Step Size:   180.0     Total Elapsed Time:   31860.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    8   Step Size:   180.0     Total Elapsed Time:   32040.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:    9   Step Size:   180.0     Total Elapsed Time:   32220.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   10   Step Size:   180.0     Total Elapsed Time:   32400.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   11   Step Size:   180.0     Total Elapsed Time:   32580.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   12   Step Size:   180.0     Total Elapsed Time:   32760.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   13   Step Size:   180.0     Total Elapsed Time:   32940.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   14   Step Size:   180.0     Total Elapsed Time:   33120.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   15   Step Size:   180.0     Total Elapsed Time:   33300.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   16   Step Size:   180.0     Total Elapsed Time:   33480.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   17   Step Size:   180.0     Total Elapsed Time:   33660.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   18   Step Size:   180.0     Total Elapsed Time:   33840.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   19   Step Size:   180.0     Total Elapsed Time:   34020.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   20   Step Size:   180.0     Total Elapsed Time:   34200.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   21   Step Size:   180.0     Total Elapsed Time:   34380.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   22   Step Size:   180.0     Total Elapsed Time:   34560.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   23   Step Size:   180.0     Total Elapsed Time:   34740.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   24   Step Size:   180.0     Total Elapsed Time:   34920.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   25   Step Size:   180.0     Total Elapsed Time:   35100.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   26   Step Size:   180.0     Total Elapsed Time:   35280.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   27   Step Size:   180.0     Total Elapsed Time:   35460.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   28   Step Size:   180.0     Total Elapsed Time:   35640.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   29   Step Size:   180.0     Total Elapsed Time:   35820.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   30   Step Size:   180.0     Total Elapsed Time:   36000.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   31   Step Size:   180.0     Total Elapsed Time:   36180.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   32   Step Size:   180.0     Total Elapsed Time:   36360.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   33   Step Size:   180.0     Total Elapsed Time:   36540.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   34   Step Size:   180.0     Total Elapsed Time:   36720.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   35   Step Size:   180.0     Total Elapsed Time:   36900.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   36   Step Size:   180.0     Total Elapsed Time:   37080.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   37   Step Size:   180.0     Total Elapsed Time:   37260.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   38   Step Size:   180.0     Total Elapsed Time:   37440.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   39   Step Size:   180.0     Total Elapsed Time:   37620.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   40   Step Size:   180.0     Total Elapsed Time:   37800.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   41   Step Size:   180.0     Total Elapsed Time:   37980.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   42   Step Size:   180.0     Total Elapsed Time:   38160.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   43   Step Size:   180.0     Total Elapsed Time:   38340.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   44   Step Size:   180.0     Total Elapsed Time:   38520.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   45   Step Size:   180.0     Total Elapsed Time:   38700.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   46   Step Size:   180.0     Total Elapsed Time:   38880.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   47   Step Size:   180.0     Total Elapsed Time:   39060.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   48   Step Size:   180.0     Total Elapsed Time:   39240.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   49   Step Size:   180.0     Total Elapsed Time:   39420.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   50   Step Size:   180.0     Total Elapsed Time:   39600.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   51   Step Size:   180.0     Total Elapsed Time:   39780.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   52   Step Size:   180.0     Total Elapsed Time:   39960.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   53   Step Size:   180.0     Total Elapsed Time:   40140.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   54   Step Size:   180.0     Total Elapsed Time:   40320.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   55   Step Size:   180.0     Total Elapsed Time:   40500.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   56   Step Size:   180.0     Total Elapsed Time:   40680.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   57   Step Size:   180.0     Total Elapsed Time:   40860.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   58   Step Size:   180.0     Total Elapsed Time:   41040.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   59   Step Size:   180.0     Total Elapsed Time:   41220.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Transport Step:   60   Step Size:   180.0     Total Elapsed Time:   41400.    \n",
      " Outer Iter.  1  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      " Outer Iter.  2  Inner Iter.  1:  Max. DC =   0.000      [K,I,J]    1    1    1\n",
      "\n",
      " STRESS PERIOD NO.    3\n",
      "\n",
      " TIME STEP NO.    1\n",
      " FROM TIME =   41400.     TO    86400.    \n",
      "\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " Transport Step:    1   Step Size:   180.0     Total Elapsed Time:   41580.    \n",
      "forrtl: error (74): floating underflow\n",
      "Image              PC                Routine            Line        Source             \n",
      "mt3dusgs           000000010B92134A  Unknown               Unknown  Unknown\n",
      "libsystem_platfor  00007FFF6D89DB5D  Unknown               Unknown  Unknown\n",
      "mt3dusgs           000000010B8B331F  _sft1fm_                  523  sft1.f\n",
      "mt3dusgs           000000010B8DAEF5  _gnt1fm_                   23  gnt1.f\n",
      "mt3dusgs           000000010B8E45E8  _MAIN__                   376  mt3d-usgs.f\n",
      "mt3dusgs           000000010B7225EE  Unknown               Unknown  Unknown\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(False, [])"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mt.write_input()\n",
    "\n",
    "# run the model\n",
    "mt.run_model()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Compare mt3d-usgs results to an analytical solution"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define a function to read SFT output file\n",
    "def load_ts_from_SFT_output(fname, nd=1):\n",
    "    f = open(fname, \"r\")\n",
    "    iline = 0\n",
    "    lst = []\n",
    "    for line in f:\n",
    "        if line.strip().split()[0].replace(\".\", \"\", 1).isdigit():\n",
    "            l = line.strip().split()\n",
    "            t = float(l[0])\n",
    "            loc = int(l[1])\n",
    "            conc = float(l[2])\n",
    "            if loc == nd:\n",
    "                lst.append([t, conc])\n",
    "\n",
    "    ts = np.array(lst)\n",
    "    f.close()\n",
    "    return ts\n",
    "\n",
    "\n",
    "# Also define a function to read OTIS output file\n",
    "def load_ts_from_otis(fname, iobs=1):\n",
    "    f = open(fname, \"r\")\n",
    "    iline = 0\n",
    "    lst = []\n",
    "    for line in f:\n",
    "        l = line.strip().split()\n",
    "        t = float(l[0])\n",
    "        val = float(l[iobs])\n",
    "        lst.append([t, val])\n",
    "\n",
    "    ts = np.array(lst)\n",
    "    f.close()\n",
    "    return ts"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load output from SFT as well as from the OTIS solution"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Model output\n",
    "fname_SFTout = os.path.join(\"data\", \"CrnkNic.sftcobs.out\")\n",
    "\n",
    "# Loading MT3D-USGS output\n",
    "ts1_mt3d = load_ts_from_SFT_output(fname_SFTout, nd=38)\n",
    "ts2_mt3d = load_ts_from_SFT_output(fname_SFTout, nd=105)\n",
    "ts3_mt3d = load_ts_from_SFT_output(fname_SFTout, nd=281)\n",
    "ts4_mt3d = load_ts_from_SFT_output(fname_SFTout, nd=433)\n",
    "ts5_mt3d = load_ts_from_SFT_output(fname_SFTout, nd=619)\n",
    "\n",
    "# OTIS results located here\n",
    "fname_OTIS = os.path.join(\n",
    "    \"..\",\n",
    "    \"data\",\n",
    "    \"mt3d_test\",\n",
    "    \"mfnwt_mt3dusgs\",\n",
    "    \"sft_crnkNic\",\n",
    "    \"OTIS_solution.out\",\n",
    ")\n",
    "\n",
    "# Loading OTIS output\n",
    "ts1_Otis = load_ts_from_otis(fname_OTIS, 1)\n",
    "ts2_Otis = load_ts_from_otis(fname_OTIS, 2)\n",
    "ts3_Otis = load_ts_from_otis(fname_OTIS, 3)\n",
    "ts4_Otis = load_ts_from_otis(fname_OTIS, 4)\n",
    "ts5_Otis = load_ts_from_otis(fname_OTIS, 5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Set up some plotting functions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "def set_plot_params():\n",
    "    import matplotlib as mpl\n",
    "    from matplotlib.font_manager import FontProperties\n",
    "\n",
    "    mpl.rcParams[\"font.sans-serif\"] = \"Arial\"\n",
    "    mpl.rcParams[\"font.serif\"] = \"Times\"\n",
    "    mpl.rcParams[\"font.cursive\"] = \"Zapf Chancery\"\n",
    "    mpl.rcParams[\"font.fantasy\"] = \"Comic Sans MS\"\n",
    "    mpl.rcParams[\"font.monospace\"] = \"Courier New\"\n",
    "    mpl.rcParams[\"pdf.compression\"] = 0\n",
    "    mpl.rcParams[\"pdf.fonttype\"] = 42\n",
    "\n",
    "    ticksize = 10\n",
    "    mpl.rcParams[\"legend.fontsize\"] = 7\n",
    "    mpl.rcParams[\"axes.labelsize\"] = 12\n",
    "    mpl.rcParams[\"xtick.labelsize\"] = ticksize\n",
    "    mpl.rcParams[\"ytick.labelsize\"] = ticksize\n",
    "    return\n",
    "\n",
    "\n",
    "def set_sizexaxis(a, fmt, sz):\n",
    "    success = 0\n",
    "    x = a.get_xticks()\n",
    "    #   print x\n",
    "    xc = np.chararray(len(x), itemsize=16)\n",
    "    for i in range(0, len(x)):\n",
    "        text = fmt % (x[i])\n",
    "        xc[i] = string.strip(string.ljust(text, 16))\n",
    "    #   print xc\n",
    "    a.set_xticklabels(xc, size=sz)\n",
    "    success = 1\n",
    "    return success\n",
    "\n",
    "\n",
    "def set_sizeyaxis(a, fmt, sz):\n",
    "    success = 0\n",
    "    y = a.get_yticks()\n",
    "    #   print y\n",
    "    yc = np.chararray(len(y), itemsize=16)\n",
    "    for i in range(0, len(y)):\n",
    "        text = fmt % (y[i])\n",
    "        yc[i] = string.strip(string.ljust(text, 16))\n",
    "    #   print yc\n",
    "    a.set_yticklabels(yc, size=sz)\n",
    "    success = 1\n",
    "    return success"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compare output:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEFCAYAAADkP4z+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABcDUlEQVR4nO2dd3hUVfrHPzOT3kMKoZPQhVAEA0iLhSJiRRQXEVBAUMRFXd3VdcUOKDZsi6DY+GFF0UVRlFCFEGpAWgihJCG9t2n398dlhpTpc9PP53nyZObOveeemUm+973vec/3qCRJkhAIBAJBq0Hd2B0QCAQCQcMihF8gEAhaGUL4BQKBoJUhhF8gEAhaGUL4BQKBoJXh0dgdsEV4eDhdu3Zt7G4IBAJBsyItLY3c3Fyrrzdp4e/atStJSUmN3Q2BQCBoVgwZMsTm6yLVIxAIBK0MIfwCgUDQymjSqR6BQNDyyMrKIj8/v7G70aLw9PQkPDyckJAQh/YXwi8QCBqU/Px8evbsiUajaeyutAgkSaKyspK0tDSHhV+kegQCQYMjRF85VCoVvr6+Th0jhF8gELRI0tPT3W6jsrKyRaalhPALBIJmjU6n4+GHH2b27Nncfvvt/P777wD885//BOCFF16w28b69eu55557mDVrFmvWrOHixYt8+OGHfPnll+zfv9/msV988QX/+c9/3H8jDYjI8QsEgmbN6tWrmTRpEuPHj0ev13PLLbcQHh7OwYMHOXDgACdOnGD79u2sXr2ayspK3nrrLdq2bVujjfXr17N8+XLatm3LunXrqKys5OzZs2RlZeHt7U1AQACfffYZ5eXlPPTQQ+Y6+T///JOUlBSMRmNjvHWXERG/QCBoFFQqlcs/1Tl06BBXX301AB4eHvTp04ewsDAGDBjAoEGDAMjMzMTLy4uZM2daHABdvnw57777LvPnz68h4iNGjODWW2/lnXfeITAwkDZt2rB9+3bz68OHD2fGjBn18OnULyLiFwgEjYJSa0D169eP3bt3M3bsWCRJIjU1laioqBr79O/fn379+vH+++9TVVXFLbfcUuP11atX8/zzzwMwZcoUhg4dCmC+yOj1ep566ikqKio4evSoIv1uTITwCwSCZs38+fN54okn+OKLL8ypGA8PD/R6PQkJCYAc8a9cuRKNRkPv3r156623uOuuu8wXiOjoaG677Tb8/f254oorzFVHXbp04YMPPuDhhx9mzpw56HQ6h8YMmjqqprz04pAhQ4RXj0DQwjh27Bh9+vRp1D5899133HDDDU6XQTZlqn+u9rRTRPwCgaDVcfvttzd2FxoVMbgrEAhaJKKO3zpC+AUCQbOmMev4v/jiC+6//34mT56syIWmoRCpHoFA0KxpzDr+wMBAVq9ezfr169m1axdTpkxpjI/AaUTELxAIGgWVyvWf6jRmHf/NN99MTk4O69evZ9KkScp/SPWEiPgFAkGjoFQ9YWPW8Z85c4ZXXnmFt99+u1lVCImIXyAQNGvmz5/Ppk2bmDlzJnfddZfVOv4XXniBgoICcx3/xYsXzW2Y6vjvueeeOnX8q1atMtfxz58/v0aaaMaMGeTn57NgwQLzuZoD9VrHv2fPHp588kkSEhI4ePAgDz/8MBqNBm9vbz799NM6ebbaiDp+gaDlIer46wdn6vjrLeJftmwZs2fPprKyEoBHHnmEFStWkJCQwO23387SpUvr69QCgUBgk9tvv71Fib6z1Jvwd+vWje+++878fN26dQwcOBCQ82U+Pj4Wj1u5ciVDhgxhyJAh5OTk1Ff3BAJBC0fU8Vun3oR/8uTJeHp6mp+3a9cOgF27dvHOO++waNEii8fNnTuXpKQkkpKSiIiIqK/uCQSCFkJj1vFv3LiRmTNncv/991NRUaHMG2oAGrSq58svv+Sll17if//7nxB1gUCgCI1Zx69Wq1m9ejXLly/nxIkT5qxGU6fBqno+//xz3nnnHRISEoiJiWmo0woEgqaKQoX8jVnHP2HCBDZu3MgPP/xA165dlf+M6okGifgNBgMLFy6kc+fOZnOkMWPG8NxzzzXE6QUCQVOkBfjx79y5k5tuuon27dvz+eefs2DBAkXeU31Tr8LftWtXdu/eDdAiB0gEAkHj05h+/GlpaXz00UeUl5c3K59+4ccvEAgaFFHHXz8IP36BQCCwgfDjFwgEghaIqOO3jhB+gUDQrGnMOn6QUyzjxo1z/400ICLVIxAImjWNWcdfVVXFhx9+2OzmJQnhFwgEjYLqOZX9nawgPXu5JuXQoUNMmzYNsF/HP3XqVKt1/CtWrCAnJ4dRo0aZt48YMYKOHTvyzjvv0LFjR3x8fNi+fbtZ+JcuXcqjjz5qvrtoLgjhFzQrjEYjCQkJ6PV6Dh8+TGxsLIcOHUKlUtG/f3+Sk5OJjY01/7b0mqP7Hz16lIULF+LhIf5N6oPq4u0OjVXHX15ezu7du8nIyCAxMZE//viDa6+9VpH3VN+Iv2hBsyE7O5vRo0dz4sQJF1uYAKwCoi89X33p8apavy+/9tZb75OUdIyICPGv0lRprDp+Pz8/Nm7cCMA999zTbEQfRB2/oBlw/PhxFi1axC+//OLC0b2AHsC7QGcXjj+NRvM/OnVayO+/g3AbcR9Rx18/iDp+QYvh5MmTjBgxwhyBOY4a+AS4x80edOPAgQfZsgVuvhl274aAADebFDQ6rb2OXwi/oMmi1WqZMmUKbdu25dixY/j5+VFeXg7A9ddfz7hx4+jfv3+NvPzhw8ns3BnLli35FBffrUAvJGJjPejXD378EX75Be64Q4FmBYJGRAi/oMmyefNmdDodJ06cQK1WU15eTs+ePXnnnXe47rrrUKvlaSjjx483H3PgwHi2boXi4rrtqVRw8SKYbh4MBvDxgcrKy79Nr73//krOn/+LN95403zsyJGwb58Q/uZCeno6HTp0cKuNyspKysvLadOmjUK9ahqICVyCJssPP/xgLq3z9PRErVbz7rvvMnbsWLPoVyctDZYsgcLCum2VlEBVFURGQliY/BMZCUFBNX+bXuvXL5LMzFT8/C63MWSILPyCpkVjTuD6/vvvmThxIjNnziQ5OVmZN9QAiIhf0CQxGo388MMPjBgxApD/uadPn26zcuKJJ2DePKi+nLNGAxUVUG0xOIfo3Lkz586dq7Ft8GBISpLdhFWul6ALLqG6VHHjClJ8vPlxY07gSkxMJCoqCl9fX7p16+by+2lohPALmiRJSUl4e3vz888/c+DAAa666iruvfdei5E+yCmcTZvkdIwJtdo10Qe5jO/s2bM1tkVFgZ+ffGcRHW35OIHjVBdvd2jMCVyTJ08mNjaWhIQE1qxZw4MPPqjIe6pvhPALmiR//vknEydOJC0tjWHDhjFt2jTibQjFp59CaChcKqsG4Px510QfoE2bNmi1WoqLiwkKCjJvHzJEjvqF8DcdGnMhlqSkJAYOHEhoaCh6vb4+36aiCOEXNEnOnj1LUlISR44cIT4+npUrV1qN9iUJVq+G6gG6hwe0b+/6+VUqFZ07d+b8+fP07dvXvH3AAEhOhilTXG9boCyNuRBLREQEc+bMQaVS8dprrzX4e3cZqQkzePDgxu6CoJEYPXq05OXlJSUnJ0u+vr7Sli1brO57+rQktW0rSSqVJMmXAUmqqHC/D+PGjZM2btxYY9t770nSAw+433Zr5q+//mrsLkjffvutVF5e3tjdUJTqn6s97RQRv6BJUl5eztixY4mLi2Pq1KmMHj3a6r7bt8OwYfDDD5e3lZTIJZru0Llz5zp5/ogIyM6utaNeL+eVMjMhNxfatJFHglvQrNCWhpjAJRA0Qc6dO8f69etJSUlh9OjRVtM8ANu2wfXXw//+J2uwh4cs0O4SGhpKca0JAZGRkJMDZGXBF1/A5s2wcyeEhEC7dnItaE4OnDkDzz0HzWSwrylhNBrZtm2b3e+9pZKRkUF7d/KUDtD6PlVBk6eiooLCwsI6A3TW2LpVrujR6+Gee+RoXwm8vLzQarU1trX1KeL+o49C795w5Ajcfz+cPi0PMOzeLV99EhNhxw5Yvhzef1+ZzrQSjEYjs2fPZuLEicyePRuj0Wj3mMWLFzN9+nTz86uuuorffvuNe+65h6uvvprrr7+e+fPns2/fPv72t78xY8YMsrKymDlzJtOmTWPatGksWLAAqZptWUJCAqtWrTI/v+eee6isrGTWrFnMmTOHGTNmUFBQQFFREXPmzGHWrFn87W9/M88hWLBgAXPmzGHatGmkpaXV6fPJkye5++67mTNnDo899hggzzc4e/Ys7733ns33m5ycXGdw2llExC9ocly4cIEOHTowd+5c1q1bx9SpU1m1apXF6M+UXfnpJ/n555/D66+7n+YBC8KfnEz3u25lb1k8nD8hh//W6NULfvtNLgO65Rb3RppbEdu2bWPdunUkJiYSFxfHvffea7Oay0RRURGVlZVcuHABX19fVCoVn3/+OWvWrKFjx45cf/317Ny5k/fee49ff/2VP//8E4BPPvkEDw8Pvv76a3788Uduvvlmq+c4c+YMISEhvPHGGxw6dIiLFy/y8ccf8+ijj5rN0TZv3kxpaSk5OTmsXbuWCxcukJGRQdeuXWu09fvvv3PXXXdx6623smHDBioqKjhx4gQ7duxg69at5OXl8cwzzwAQFxfHzJkzASgsLOSrr77Cx80/cBHxC5ocZ8+eJSgoyCwA69atY9u2bRb33bMHhg+X0zugXJoHZOHX6XTykyNH4PrrUT3/PLOMq9GF2hB9EzEx8i3Im28q06FWwOjRo5k6dapDYzvVueGGG9i4cSPfffed1fz9iBEj+Ouvv3jttdfM9f0mYmNja9h9q2rN0FOpVPTp04f4+HgWLlzI6tWrCQ8PJy0tjT59+nD8+HFmzpzJypUrCQgI4MEHH+Sxxx5jyZIldSaLAcyZM4dz586xYMECDh48iJeXl7mPY8aMMU8iCwoKMl+kAEJCQnjhhRfwdLVO+RJC+AVNjnPnztG/f3+HBODoUTh3Tvk0D1SL+HNy4IYb4M03UU+fRps2kJfnYCOPPy7XmhYVKdexFoxarWbVqlVs3LjR6l2eJa699lq2bt1KUVGRxQlaAHv37mXIkCH8/PPPrFixosZrSUlJ9OjRg3nz5jFt2jRCQkLIysoCZM+f0NBQduzYQZs2bXj77be59957+fDDD+nYsSN79+6ld+/erFmzBo1Gw6lTp8jNzeXNN9/khRdeYGn1qeSX+PDDD5kxYwbvvPMOfn5+7NmzB7h8wTEYDEyZMoUlS5Y4dMfjLFZTPa+//rrNAx999FHFOyMQgCz8Xbt2ZfHixdx77702B/kOHYK//pIfK5nmAdkfSFtVBfPnw9SpcLfs9mmq7HFoCKJzZ4iLkweBJ09WpmMtHLVa7bTYaTQaVCoVsbGxVFZWWtynuLiY++67D0mSePjhh/nggw+YMWMGkiQRGBjIe++9x6233grIYw2rVq1i+vTpFBQUsHTpUjp16sT8+fNZuXIl5eXlvPzyy7Rt25ZHHnmEiooK9Ho9kyZNIjo6mldffZX169dTXl7O448/ztatW6moqGDChAkAjBo1invvvZegoCDUarV5xm9oaCjbt2/n008/5bHHHmPt2rVMnDjR5c/SGlYXYpk5cybffvstU6ZMofYuKpWKjz76SPHO1EYsxNI6ue+++xg+fDhz5syxu+/QobJxmsEgp3lMmRkl+O9//4v+u+946Px52L/ffEW55hr497/huuscbOjNN+Wr08qVynWuGdMUFmJpaM6cOUNmZiZXX311vZ1DkYVY1qxZw/nz5xk3bhxTp05VvpcCgRUKCgoICwuzu58kwbFjcNtt8M03cNdd8gXA6TVbrODl4cHEvXvhk09q3EaYSzodZfx4eOMN4e7WiomOjia6Cfl82Eygvfvuu+zcubOh+iIQAFBaWoqfnx8JCQk2y/kuXgQvL1n0QS6rz89Xrh/dDh9Gp1LBpEk1tlucxGWL3r3BaAQLawVLksQvKb8w+avJjPp4FM8lPEdBRYGbPRc0ZzIyMur9HDaFv3fv3nUGQQSC+qa0tJQVK1bYreU+dgyuuKJ+KnoAev/2G+t79aoTpUdEOBnxq1Qwdixs2VJjc6W+kunrp7No0yJu7HEjz455lpSCFOI/iaewstDt/jdXJAmOH5d/O0JrquNPTk7m7rvv5u677+bXX3917AOygFN1/IsXL2bx4sUun0wgcISsrCwOHjzI3r17bdZyHz8uB9O//y6ndwwGBTtx6hRBGRn82b17nZciI+VBZae48ko4eND8VGfQccu6Wwj2Dmbf3H34ecorvlwXfR2LNi1i6jdT+Xnaz3XKCls6BgPMmgWffQbTp8PHHzuWumstdfwFBQW8//77SJLEv//9b8aNG+f8h4yT5ZwbNmxw6SQCgTNIksSNN95ot5Tz1Cno0UOeOKtSuW7BbJGPPybj2mspt2C163SqB2Rbz2pXiwf/9yBeGi/WTl5rFn2QCydeG/caZ4vOsun0Jld732zJy5NFH+TfjqbuWksd/+jRowkICOCZZ55hwYIFjn04FnBK+K0UAAkEilJeXs7bb79tt5b7/Hl5AfQ+fWDGDAUjfkmCL74g64YbLk/gqobTqR6A/v3lSWAGA18f/Zpt57bxf5P/Dw913ZtuD7UHL1/7Mk9ufhKjZN+yoCURGela6q611PFXVFSwcOFCHnjgAbcqo5xK9di6DRIIlKK0tJSgoCC7tdxnz8qLooAcHS5frlCO/8gR0GjQ9exZx6sH5PV5LS3mbpPgYIiIIOfoXhb8uoAf7/6RAK8Aq7vf2vtWntnyDDvP7WRUl1FW92uJlJc7n7prLXX8r7zyCkeOHGHp0qWMGDGC+fPnu/Yhu+L7/M9//tOVw5xG+PG3PgwGg6RSqSSDwWB3386dJcnDQ/bf9/BQsBMvvyxJDz8s7dq1Sxo2bFidl1NSJCk62oV2b7lFmv/WOGnhxoUO7b50x1Jp9g+zXThR06Yp+PE3NKmpqdLOnTvr9RzO+PG7ZNnw7rvvunaVEQjsUF5ejq+vr92p+kajbNBWUCBHhuXlCnbip5/gxhstunMCBAa6EPEDfw1ozzf5O3g2/lmH9p8WO41vj31Lha7C+ZMJmhTR0dH1OnnLWVwSfknk+gX1RGlpKf7+/nb3y86Wsyf+/nDy5OW8sNvk58trK44ZY1P4XfEEeik0mccyutLGt41D+3cI6sCV7a7k55SfnT+ZoNnS6HX81mhtJWaChqOsrIyAAOu5bxPp6bLT8YwZCg/u7tole+v4+NR056yGj498x1FV5XizqQWpbKo4wnwnHUgm9pjIr6ddr9duDbSmOv5z585x5513MnPmTI4dO+byZ2Y1Tvruu+8sbpckCYOD/2F79uzhySefJCEhgZSUFGbOnIlKpaJfv368++67rXJ1HYFtTBF/QkKCTXO29HR5ILd66Z8ig7s7d8LIkQBEvveexYhfpboc9Xt7O9bs63++ztwB9xH00n+dsm4Y120cKxJXIElS6wq4Fi+WfxyktdTxFxcX8+qrr5KVlcWmTZtcruyxqrwrVqyw+PPOO+8QFxdnt+Fly5Yxe/Zs8wj7o48+yosvvsj27duRJIkfqi+QKhBcoqSkhKysLLuzdtPToVu3epi1u2MHjBgBixcT+vbbpJ45Y1GAgoIcT/eUactYm7yWB0cukm8XnJgE0DeiL1X6Kk4XnHb4mGbP4sXyspUqlcPi31rq+Pv160dlZSULFy5k5KUAxRWsCv+WLVts/tijW7duNe4a9u3bx5gxYwD5S9q8ebPLnRa0XP7880/y8vLsLsCSng4dOsiDuooN7lZVwYED8srtlwQnMiLCovg4k+f/6uhXjOg8go5BHeWr1WnHRVylUjGu2zg2pbSiyVymz1uSHBb+1lLHf+zYMSIiIti5cycffPCBQ5+NJept6cXJkyfXyG1Vv1UNDAykyMrCFCtXrmTlJfvaHKdnyQiaOzExMXTs2NHurN30dDkjY5qtq0jWcN8+6NlTVvVLgpOdk2Mx7eBMZc+H+z/kXyP/JT/p1g1SU8GJCo9rul7Dzyk/81DcQw4f0+x51rHKJxOtpY5fr9fz4IMP4uPjw2233ebCB3uJ+qoplSRJOnPmjDR06FBJkiSpQ4cO5u3ff/+99NBDD9k9XtTxtz4+/fRT6W9/+5u0ZcsWm7X8Y8dK0saNknTsmCQZjQqd/M03JWnevMvPQfLz87O46/jx8vntkZqfKoUvC5e0eq284emnJWnxYqe6dTznuNTljS5OHdOUEXX89UO91/G7wqBBg0hISADg559/ZtSo1jUbUeAYpaWlBAYGEh8fb3PwPyMD3ntP4YqeQ4dg4ED58eLFGP79b4uDu+B4qufrv77m9t6346m5dGtiividoEdYD4qqirhYetGp4wRNh6ZWx2831VM7x6pSqfDz86NHjx4EBQU5fKLly5czZ84ctFotffr04Y477nC+t4IWT1lZmcN1/EePyo8Vq+g5dAgeeMA8uKgBnka+7a99EXI01fPV0a9Yen21HG9MjLwGrxOoVWqGdhjKngt7uKX3LU4dK2h+ZGRk0L59+3o9h13hX7RoEYcOHaJv375oNBqSk5Np164d5eXlrF69mltusf6H2LVrV3bv3g1Az5492bp1q3I9F7RISktL7dbxSxIUFsqVPHq9QhU9Op1s8N+vH/z8s/lEr3h78y+dDu9adZuOVPWcKTjDuaJzjOk65vJGJwd3TQztMJTdF3YL4bfA4sWLOX36NJ9dqu296qqrePnll/nkk09ITU01B6qzZ89m+fLleHp6smzZMp588knzPI3Q0FBWrFhhHoc0laDPnj0bkOv4V61axfz58/Hw8ECr1fLmm2+iVqt5/PHH0ev1VFVVcf/993PdddexYMECqqqqKC8v56WXXqpTznny5EmeffZZAgICCAoKYvny5bzwwgvmQeMXX3zR6vutqKhg8ODBHD58GA8XZy7aPapLly4sX77cPLKcmJjIG2+8wauvvsrNN99sU/gFAmcpKysjwo6Kl5bKg7qFhQr68J84IS+M/uqrcikhwOLF5tm7tYXfkYh/46mNTOwxsaYDZ7t2svewVisvH+YgwzoOY9muZQ7v35zJLMlk1g+zWHPrGqICHFnRvvXU8QMsWbKkTnvOYjfHn5qaWqOcKC4ujpMnT9KxY0e3TiwQWMKRiD8vD8LDZfFXqxXy4T94UPbMr1VK6OnpaXH2riMR/8YUWfhroNFA27ay0ZATDG4/mAOZB1qFXcqynctITE9k2U7HL3StpY5/w4YNxMXFERkZ6fBnYwm7wu/p6Vljia9ff/0VLy8vcnJyLP5DCATu4Ijw5+aCA2uxO8ehQ7LwA4y5nJpx1a+nQlfBtrPbGBsztu6LHTvChQtOdS/SPxIfDx/OF593+BhJkjhbWcnRsrJmdcHYm7GXgsoCEtMTHT6mtdTx//jjj/z444/8+eef5rJ3V7Cb6nnvvfeYPHkyKpUKo9GIj48P33zzDa+++irz5s1z+cQCgSUcGdzNy6sH4f/rL5g7V474t26VZ40++6zLDp0JaQkMihpEqG9o3Rc7dpQnIjjJgKgBHLp4iM7Bne3ue6q8nOnHjnG2qgpvlQoftZqv+valvwM+SI3Njvt2OH1Ma6nj//DDDwGYOXMmc+fOdfpzMuNIfaher5cOHDggJScnS3q93uU6U2cRdfytj7Fjx0q//PKLzX2++EKS7rxT4Rr+bt3kBiVJNvg3b+4mnTp1qs7uGzZI0o03Wm/usU2PSc8nPG/5xb//XZKWL3e6i//49R/Si1tftLvfybIyqe2OHdJb589LBqNRMhqN0tqLF6WIHTukpOJip8+rNKKOv35wpo7fbsR/8eJF1qxZQ36txS+XLWsdA02ChsWRiD8nR07J9+nj3ILcVtFq5dRLTIz8vNqsUWsOnfZSPQlpCbw54U3LL7qQ6gHo37Y/G07YXve6RK/nxuRkXoiOZk61ksC727ZFo1Jx919/sX/wYAIU87EWOEJ0dDTR0dGN3Q0zdnP8N998M4mJiUiSVONHIKgPHMnxX7gge/CDcwtyW+X0abmix8urjj2Dp6en06mewspCTuSd4Kr2V1newUXhH9B2AIezDtvc5z9paVwdFFRD9E3cGRnJ1UFBPGvBJljQdGgIP367l32tVmvVolkgUBpH/PgrKuRqHqNRoRr+Eydkjx6TK+Rzz8lRf7VyztrYqurZcW4HQzsMxdvDimezi8LfK7wXZ4vOUq4rx8/Tr87rh0pLWZuVxdGrrFxwgCUxMfTdu5dHO3Wig6Oe0k2c1lTHX1hYyHXXXUdsbCw33ngjU6ZMcekzsyv8gwcP5siRI/Tr18+lEwgEzlBaWoqPj49NP/68PPjoIznNo1gNf69el4W/2h2tK4O7CWkJxHeNt34+F4XfS+NFt9BunMg9waB2g+q8/p8zZ3i6SxfCbcwPiPL2Zna7drx09izv9ezpdB8agsdTUliRns7DHTrwWvfuDh3TWur4ExMTzfNcBprsRVzAbqpnxIgRDBw4kE6dOhETE2P+EQjqg6qqKp588kmbfvx5eXIpvGI1/CbhhxqlnGBd+G1F/DvP72RkZxte6e3aQVaWS1etPhF9OJ57vM72AyUl7C0pYU67dnbbeKxTJ9ZmZVHQRMuxV6Sns6l/f1Y4UfnUWur4e/XqxUcffcSbb77JSy+95PDnUxu7wv/qq6+ydu1avv/+e77++mvzj0BQH1RWVrJ+/XqbfvyKl3OePHk51WMq5byU57c2gcvXV3Z5qP1Slb6Kw1mHref3QR5LCA11akEWE33C+3Ast+6Se0vPnePxTp3wdWCUO9LLixvCwvj0Up16U+PhDh0Yf/gwD3fo4PAxraWO/8SJE5SXl+Pv7+/Wimx2Uz0hISHceeedLp9AIHAGo9HIXXfdZdOPX3HhP3UKevSA+HiHUz0qlbzQe1kZVNeZgxcP0qNND/y97BjNdegg1/I7EKFXp094H747XnPMLbOqik0FBfzXdNfiAPPbt+eBkydZ2KFDk1vS8bXu3R1O8ZhoLXX80dHR/Pvf/8bf358FCxa48Olewl5t6FNPPSU99thj0q5du6R9+/aZfxoCUcff+lCpVJJWq7Xpxx8QIEkFBQrV8ZeVSZK3tySZzvXsszVevu2226Rvv/3W4qHt20vS+fM1t721+y3pgR8fsH/eiRMl6aefnO7u/oz9Ur/3+tXY9tyZM9IDx4871Y7RaJS6794tJRYVOd0HdxF1/PWDonX8a9euBeDbb781b1OpVKQ66SkuENjDcCnn7enpWeP2tjpVVXJVz8KFcimn23X8aWnQpcvlJbxqrbRlLeIHCAiQI/7q7L6wm3Hdxtk/b9u2cNF5f/1e4b1IyU9Bb9TjofbAIEl8mJnJT7GxTrWjUqmYGhnJl9nZXOWEvbrANZpdHf+ZM2fq/AjRF9QHOp0OTzujtQUFcmrlUuWe+3X8Z86AqeLCwvqutoTf3192Cq3O7gu7GdphqP3ztm0rD/A6iZ+nH1EBUaQVpgGwtbCQcE9PBrhgxXBXRARf5eRgbIR5OQZFyrEEIHsyVVRUOHWMmL4naDI4IvzFxfK4aFGRQl78aWkQHW2xhh+sD+6CHPFXF/78inxyy3PpFe5Arj0qyiVffoDe4b05lnOM7m2681lWFtMtVI04Qr+AAAI1GnYXF3N1cLBLbbhCmzZtOGmagSdQBE9PT9o5MV4khF/QZHBE+IuKIDgYyssV8uI3Rfz//GedgV2wn+qpLvwHMg8wIGoAapUDK5q2bQu7drnU5Z5tenIq/xQVBgPf5+byshsphNvCw/kxL69Bhb9t27YWSxwFDUeDrbkrENjD0Yg/KEhBL35TxA81PHpM2Ev1VM/xH7h4gCujrnTsvC7m+AG6t+nO6fzT/FpQwKCAANq5MQP3xrAw/peX5/LxguaJEH5Bk8HRiF/RsUhTxF/Lo8eEMxH//sz9XNnOQeGPinIpxw+y8KcUpPBDbi63hIe71IaJuKAgMrVazlopgRS0TFwS/htvvJG7776bxETHF0oQCOzhaMSvaFYiLQ2++UZO81SbuGXCmkkbWEj1XDxg0UrBIi4O7oIs/KfyT/NTXh63uDmhQaNSMbFNG34SUX+rwiXhnzdvHp9//jl+fnWNogQCV2nwiL+4GCorwWQxfmm5xepYs2WGmqmeUm0p54rO0Se8j2PnDg2VBypciLS7hHThPEFEeXnS1dfX6eNrM75NG35z2+JU0JxwSfhvuukmNBqNMG4TKEqDR/znzsl2zJdW27KEo6me5Kxk+oT3wVPj4KCDSgWRkS7ZNnhpvPCPGssoP3cWIbjMNSEhbCsqwiDs1lsNdqt6Tpw4wZIlS8jPz6/hw79hg+0FIQQCZ3FU+MPD4fhx2VfNLbeBCxegUyf5sYX8PtgX/pwc+fGR7CPEtnVuEpV5gLez/aUUa6MLjaMnOU4fZ4l23t609fLiUGkpVwYGKtKmoGljV/hnzJhBXFwcY8aMaXKeHoKWhU6nM7sUWqOgALZvl6sv3Z61e+GCbJEMLg/umlI9ydnJxEY6KfwuDvAeLytDpfFFXZri9LHWuCYkhC2FhUL4Wwl2hb+srIy33367IfoiaOU4EvHn5sIlI0M++wyWL3djApdJ+K1M3gLbg7vVZ+4eyT7CpJ6TnDu/iyWdP+Tl0d+jlNMFygr/Jxcv8pjpDkjQorGb4+/RoweZmZkN0RdBK8cR4a+quhzhuz1rt7rwg9ODu9Vz/Eeyj9Av0skxLxcj/o15eVwX7E9KvnLCHx8Swo6iIvQW1j8QtDzsRvxGo5G+ffsyePBgfKtVEIgcv0BptFqtQzn+zZth9GgFZu1euACmRTtcHNwtK4Os0iwMkoF2Ac5ZLNO2rWwJ7QRlBgP7Skp4rWcXvt6unPBHeHnRyceH/aWlxAnTthaPXeG/7bbbuO222xqiL4JWjqPlnGFh8qxdC6syOkftHL8FHDFpM0X7To+BtW0rD1g4wc6iIgYFBtIvrBtphWkYjAY0auWqe7YUFgrhbwXY/deZMWOG2SJXp9MxYsQIZsyYUd/9ErRCnLFsUITaqR4LODKB60j2EfpFuFDa7EKq5/eCAq4NCcHX05cI/wguFDu/dq81rgkJYUtBgWLtCZoudoV/06ZNDBkyhO+//54NGzZw1VVX8cMPPzRE3wStjAat4y8pAa0W3nrL6qxdcKyO/3jucXqH93a+Dy4M7v5eUMB1oaEAdAvtpmief0xICLuKi9GKPH+Lx67wP/PMM2zdupX169ezYcMGduzYwWIbEZJA4Cr2hN9olPVakYrD9HQ52n/uOfm5hYFdsG3LbJq5ezL/pGNWzLVxMuIv0Ok4WVHBsEu3PN3bdFdU+Nt4etLVx4dDtRcZELQ47Aq/VqvliiuuMD/v27evWERBUC/YE/6yMvDzc6NuvzrV8/tWBnbBMT/+E7kn6BXmgvAHB19eUswBEgoLuTooCK9LgxtKCz/A8KAg/iwuVrRNQdPDrvD7+vqSlJRkfp6UlCQ8egT1gk6nw8PDg4SEBIwW0g2K+vRcuCAveA52c/zWhN/LC4wepRRUFtAp2IX6d5Ntg4NRf/U0D1x26VSSq4OD2VVUpGibgqaHXeFftmwZkyZNYtSoUYwaNYqJEyeydOnShuiboJVRVVXFrl27mDhxIrNnz64j/qaB3ePH66yX4jwZGdC+vd3dbAm/SgW+HU8SHdzdscVXLOFEuuf3wsIawq90jh9ExN9asPvXOmrUKP766y+eeuopnnzySY4dO8aYMWMaom+CVsbx48dJTU0lMTGRdevWsW3bthqv5+fLM3f79IEZM9ys48/MhHbtbEb7YFv4ATzbnaBrgAtpHhMODvBmVFWRrdUysNraujGhMaQWpNbw0HKXHr6+lBkMZFRVKdamoOlhVfg///xzAF5//XXWrFnDsWPHOHnyJJ988gmvv/56g3VQ0Hro3LkzvXr1Ii4ujqlTpzJ69Ogar6eny8IPCiyynpkJO3bYrOgB21U9AOqIE3TydUP4HYz4txcVMSo4GHW1uQLBPsH4eviSXea8w6c1VCoVw4ODRdTfwrE6gevUpRmFycnJdV4TZm2C+sBgMDBx4kQmTpzI6NGjUdeaoaVWyxotSQrYNWRmwiuvwLff2swb2Yv4jaEnaOd1o+v9cDDi31lUxAgLdazd2nQjtSCVtgHKrWE7PCiIXUVFTHbrAxY0ZawK/3OXytxuvfVWbrnllhqvffbZZ/XbK0GrxOTOaZowWJvS0suOnG4XlmVmyjl+GxU9YF/4dcEniFAvcr0fbdvKgxZ22FlUxIoePepsN6V7hnca7nofajE8KIinz5xRrD1B08Oq8P/444/odDr+8Y9/YDQazXlEnU7Hs88+y/Tp050+mU6nY8aMGaSlpaHRaPjwww/p3duFiS+CFom9cs7SUnlw1227BklSJMcvSRIVfidpI7mZ6klIsLlLqV7P8fJyBluYwBATIgu/klwVGMih0lKqjEa83fbFEDRFrAr/wYMH+eOPP8jOzq5hy+zh4cGiRa5FOBs3bkSv17Nr1y5+++03nn76ab799luX2hK0PHQ6HQHVBi9rU1Ym1867TVEReHrKM7DsYEv4M0oy8JQCoNKNqcRt29pdhWtPSQkDAwIsinBMaAzbzznn92OPAA8Pevn5sb+khOGKLnAsaCpYFf5nnnmGZ555hvfee48HH3xQkZP17NkTvV6P0WikuLjY7vR8QevCnjtnaalDWm0fU7QPVhdgMWFL+E/knSBY3wu3Jro6UMdvLb8PsvB/cugTNzpgGVNZpxD+loldd87Zs2ezfv16SktLkSQJg8FASkoKL730ktMnCwgIIC0tjd69e5Obm8tPP/1UZ5+VK1eycuVKAHJylFlaTtA8cCTVEx6uwImqp3msLMBiwqbw554gHDeFv21bh4T/QdNks1qYcvxKMzwoiA15eYq3K2ga2BX+u+66i9TUVDIzMxk0aBB79uyxOvhmjzfeeIPx48fzyiuvcP78ea699lqSk5Px8fEx7zN37lzmzp0LwJAhQ1w6j6B54ojwK5LqqS38Llb1nMg7QVsPN4U/JAQqK+Wfav8HJgySxO7iYj7v08fi4R2DOpJTnkOlvhIfj7rHu0qcGOBt0dgduTl48CD79u3jlltu4c0332Tnzp3ku1hAHRoaSvClW8c2bdqg0+mE74/AjCNePYqnetyo6jmRd4L2Xr3M6+66hMm2wUqe/0hZGVFeXkRYWYtYo9bQObgzaYVpbnSiLj18fSkxGLgoJnK1SOwKf7t27fDw8KBnz54cOXKEvn37UuSil8eiRYvYv38/o0aN4tprr+Xll1/GX5H/ZEFLoMEjfrBb1ePh4YHBYLA4O/ZE7gk6+/d0L+IHm3l+W/l9E/WR7lGpVMQFBpJYUqJou4Kmgd1UT0BAAGvXrmXAgAHm8stSF//SAwIC+Oqrr1w6VtDyaVDhHzjQoV1VKhUeHh7mOQYmtAYt6SXpdI6KZo+7wm8jz7+zqIhrq/nzWKJbaLd6yfPHBQWxp7iYmxUZWBE0JexG/O+++y4HDx5k7NixqNVqxowZwz/+8Y+G6JuglVFbXGtjsmV226QtK0uun3dwXQlL6Z6zhWfpGNSR4ABP91I9YLOkM7GkxOy/b436GuAdGhQkIv4Wil3h/+ijj1i2bBkqlYovv/yS/Px85s2b1xB9E7Qy7EX8JSXw8ssKmLRlZ8N339n16TFhSfjPFJ4hJjTG7MnvFlZSPQU6HZlVVfS2Y4NeX8IfFxjI3uJijAqawAmaBnaF31LJpUBQH9gT/uJi2LBBfuyWSVtW1uVBXSsrb1XH0rq7qQWpRIdEKyP8VlI9+0tLGRQYiMaON1Z9CX+ElxdtPD05WV6ueNuCxsVujj8mJoZx48YxcuTIGrMqH3300XrtmKD1YU/4Kytlcza93g2TNr0eCgrkCQF2KnpMeHl51Yn4UwtSiQmNwd8bZVI91RY7MpFUUsIQB9aZjA6JNtszK22gODQoiD0lJfQWRRgtCrvC36ZNGwDOVKvpFe6cgvrAkcHdggI5z+9ymicvD0JD5SuHGzn+1IJUhrQfQoC6/lI9+0pKuDkszO7hwT7B+Hj4kFOeQ6R/pJudqUlcYCB7iouZERWlaLuCxsWu8At3TkFDYUv4dTp5sXV/fzkt77J3WFaWLLROYE34Y0JjCDDUX6onqaSE57t2daiJmNAYTuefVlz4hwYF8YUTC8ILmgcN6s4pENjClvCbJm+5fbOZlSULrRPUFn5JkjhdcJqY0BhUlQqlemqJa55OR65OR08H17euD3tmgEEBARwrL6fCYMBXkVXuBU2BBnXnFAhsYcukTbEa/uxsWWjtmLNVp7bwF1QWoEJFqE8oOo3cN0ly46IUHg6FhZcHL5DTPFcGBNRYccsW9TXA66vR0NvPjwOlpVwtDNtaDA3qzikQ2MJWxK+Y8GdlQWoq/N//2TRnq05t4TeleVQqFV5esuBrteDt7WKfNBp53CE3V55fgOMDuyZiQmPYeX6nix2wzdCgIBKLi4XwtyDs5vhnzJjBJ598Qn5+fo1p66KqR6A09oRfkcKSrCy45RbYs8fhWWDWhN9EQICc7nFZ+OFyuqea8N/pRNlSt9BufHa4fsbe4gID2VRQUC9tCxoHh4T/zJkzxMbGmqt5RFWPoD6wl+NXLNXTq5fDpZxgWfijQ6LNz/395QvTpQI416g1ezeppIRlMTE2DqhJfaV6QI74nz97tl7aFjQOdoX/0KFDHDt2DA8Pu7sKBG5hy7JB0VRPZCTcd5/Dh9QW/jMFZxgYNdD8XOnZu1laLSUGA918fR0+vGNQR3LKlLdnBujl50e+TkeOVmvVJVTQvLBbFNepU6eG6IdA0HCpHjerelIL66Z6lCzp3FdSwuCAAKfurDVqDZ2CO3G2UPnIXK1SMUQ4dbYo7IbxsbGxXHPNNUyYMAHfahGIyPELlMaRcs7jx+VMjcvZRgWqemrn+P39FSrpvJTqcXZg14Qp3dMr3I3F361gGuC90YEJZYKmj92Iv7i4mO7du5OSkkJycjLJyckcOXKkIfomaGXYEv7iYvjzTzcN2iRJFteVKx02aIOawq836rlQfIEuIV3Mryud6tnnqvCH1G+ef09xcb20LWh47Eb8H3/8MQCFhYWEhITUd38ErRhbwp+dDSdOyI8/+wyWL3fBq6ewUF7e8MUX4aWXXKrqOV90nqiAKLw0l3PdSqd6kkpKeLN7d6ebqM8B3rjAQGaWlNSLH5Cg4bEb8Z88eZIrrriCvn37kpGRQZ8+fTh+/HhD9E3QyrAl/NVtGlw2aDOlecDpqh6TO2ftNA8omOrJyiKjqopKo5GuFtbftUdMaAyphfUj/FHe3gRqNKRUVNRL+4KGxa7wL1iwgLfeeovIyEjat2/Pww8/bF4MXSBQCkmS0Ov1qNVqEhISMBqNNV4vLYVXX5VTPC67BFf36XEwvw813TlTC1KJCakp/IpF/NnZ5jSPK1F1fUb8cHlFLkHzx67w5+XlMXbsWPPzBx98kGLx5QsURq/Xo9FomDNnDhMnTmT27Nk1xL+0FIKC5KjfhoGnbVyo6IGaqZ7UglSiQ6NrvK6I8EdEQHa2ywO7cNmozdL6wEowNDCQPaKyp0VgV/hVKhWVlZXmCOTixYsYXPbEFQgso9Pp0Gg0rFu3jsTERNatW8e2bdvMrysygat6qscJagh/YT2lenx8wM+PpPx8l4W/uj1zfWCq7BE0f+wK/4MPPsj48ePJzs7mX//6F8OGDRPePQLF0el0eHt7M3XqVOLi4pg6dSqjR482v67IBC5TqseJNA/UFP4zBWfqCH9goFx15C5S27YklZW5LPxQv+meKwMDOVJWRlWtNJyg+WFX+O+77z5eeOEFpk2bhk6nY9WqVcyfP78h+iZoRWi1Wry8vFi1ahUbN25k1apVqKuZ7isygSsrC/budaqUE+qmemoLf3CwMsJ/oUcPMBrp6IbpT30Kv79GQw9fXw65ndcSNDZ2hb+kpISdO3eydOlSHn74YX766SfK3L6vFQhqYrJrUKvVxMfH1xB9UCjiz86G+++XHzuw1q4Jk/AXVRZRqa8kwq9mSVFQkDLCn9S3L0OqpVVdob4HeEU9f8vArvDPnDmTvLw8AEJCQlCpVMyZM6feOyZoXdhbdlGRHL8p1eNEKSdcFv4zhWfMdszVCQ6GoiI3+wYkxcQwxE0XTFHZI3AEu8J/6tQpXnvtNQCCg4N54403OHr0aL13TNC6cGS9XUVSPSa7BicwCb+lNA8oF/Hva9uWIenpbrVR7xG/8OxpEdgVfp1OV6N8s7S0tN7KxQStF0eEv7GremrbMZsICnI/4pckiaTAQAafOuVWOzGhMZwuOO1eZ2zQx9+fi1ot+bXWIBY0L+xaNtx7770MHTqUKVOmoFKp+O6775g1a1ZD9E3QinAk1eOWSVt5ubxiuwsVM9WFv094nzqvKzG4e7ayEi+1mvan3RPtTkGdyC3PpUxbhr+XEnamNdGoVAwODCSxuJgJwrCt2WI34v/Xv/7F0qVLKSoqorS0lGXLlglnToHi2PLi12rBaIQ5c9wwaTOleVwYOHUk1VNU5LD1j0WSSkoY4uUFmZmuN4Jsz9wttBun8t27c7CFSPc0f+wKP8CECRNYvHgxTz/9NEOHDiU/P7+++yVoZZjKOS1hyu9/dmllwc8+A6f/BKvP2lU4x+/jI19Pqqqc7FM1kkpKGBIa6rbwA/QK78WJ3BNut2MNMcDb/LEr/B988AHBwcGEh4cTERFh/i0QKIkt4S8rk6Nq0yJwLpm0VR/YdaGOv0pXxbmic3QN6WpxH3cHeJNKShgSGSm/2cpK1xsCeoX14kRe/Qm/qaRTjPU1X+wK/6uvvsrOnTsxGAwYDAaMRqOwbBAojlartbv6Vnm5GyZtpoXMTWLvZB1/kbGIML8wfD0tL4foTkmnJEnsKy1lcFCQfHG6eNG1hi5R38LfwdsbH7WaVDcvUILGw67wR0VFceWVVzZEXwStGHupnoAA2ZzNZZO2ixddsmQG2Z2zUFVoMc1jwp2I/3RFBYEaDW29vKBdO7fTPfWd6gG4OjiYnUpMXhA0CnaFf9y4cbz//vukp6eTn59v/hEIlMTW4K5ik7eiouTHLuT4i9XFdoXfVR2s4ciphPCH9eJk3sl6TcWMCg5muxD+Zovdcs4lS5ZQVVXFQw89ZN6mUqlEukegKI4M7rrFxYsQH+/SoZ6enpR4lFis4TfhTknnvtJSRYU/1DcUHw8fMkszaR/Y3q22rDEyOJj3MjLqpW1B/WM34q+oqMBoNNb4EaIvUBpHUj1u4aIXP8jCX+pZ2nARvwKCWt/pnv4BAWRWVZFzaWUyQfPCrvAbjUaWLVvGNddcw8iRI3n++efR6/UN0TdBK8Le4K7bwn/xYs3BXSfw9PSkzLvMpvC7GvEbJYn9JSUMNgl/x47gpm0D1P8Ar0alYrjI8zdbHJrA9ccff/DII4/w6KOPsmvXLh5//PGG6JugFWGvnFMRn56PPnK6lBNk4a/wqaiXwd1TFRW08fQkzHTR69gRLlxwvqFa9Aqr/wHekSLP32yxm+P/5ZdfSEpKMkdjN954IwMGDKj3jglaF7YGd92O+MvKQK+HV16BpUudnmKrU+nQa/REBURZ3Sc4GC6Z2DpFnaUWlRL+8F5sSdvidju2GBUczONuWkwIGgeHUj3Vb8G9vb1teqrY45VXXmH48OEMHjyY1atXu9yOoGVRrzn+6nYNTpZyAmTpsvAu90atsv7v4mrEX0f4O3WShd/Nipz6TvUAxAUGcrSsjDIx5tfssCv8AwcOZNGiRZw+fZrU1FQeffRR+vfv79LJEhIS2LVrFzt37mTr1q2cP3/epXYELQ9bOX63Uz1u2DUAXKy6iFeJ5YuSCVcHd+sIf2AgaDRu233GhMaQXpxOld4NHwk7+Gg0DAwIYLewb2h22BX+d999l4KCAq6++mqGDh1KTk4OK1ascOlkmzZtIjY2lttuu42bbrqJSZMmudSOoOXhSDnn8eMuBsKmgV0XyajMwKPEdlbUlcFdgyRxoKSEK2vfziiQ7vHUeNIlpAsp+SlutWOPUcHBbC8srNdzCJTHbo4/KCiINWvWAFBZWYmPj4/LJ8vNzeXs2bP89NNPnDlzhptvvpnjx4/XWNFo5cqVrFy5EoCcnByXzyVoXmi1Wvz8/Cy+VlICn34qu3NOnw4ffywHxQ7jRiknQHp5Ouoi2zGSKxH/8fJy2nl7E1r7Tsck/P36OdnTmpjSPX0j+7rVji1Gh4Twqrhzb3ZY/WvWarXMmDGD9evXm7fdcccdzJo1y+VyzrCwMMaPH4+Xlxe9evXCx8enjrjPnTuXpKQkkpKShBlcK8LW4G5BAWzfLj92yZnTjVJOgPNl51EX2Bb+kBBwNvCtk+YxodAAb5/wPhzLOeZ2O7YYFRzM3uJiykWev1lh9a/5P//5D8XFxYwYMcK87b///S8FBQUsdvEfaOTIkfzyyy9IkkRGRgZlZWWEicUcBNjO8et0lyN8l505k5JcKuUEOF96Hinfdo6pTRvnL0hJJSUMtjRqrZDw94vsx5GcI263Y4tADw8GBQayQ5R1NiusCv9PP/3E2rVriYyMNG/r0KEDn376aY27AGeYNGkSgwYNIi4ujptuuol3330XjVP37IKWir0c/44dbjpz3n+//NgJV04ArUFLVnkWxgKjzf3CwuRyTmfGIBKLixkaFFT3BSWFP7t+hR/g+tBQfndzkXhBw2I1x+/l5YWvb10L2qCgILy9vV0+4bJly1w+VtBysTeBKzhYduZUO7R0UC1MzpwulHKmFabRPqA9BVW2hc3XV74rcdRQrspoJLmsjCutpXq++cbpvtamT0QfUvJT0Bq0eGlsVyW5w/WhoSx0c61gQcNi9d9Io9FQYmF5tZKSEnRioWWBwtgS/pISl5bKvUz1RVic5HT+aWJCYxz6mw8Pd3wS16HSUnr4+uJv6Y63c2c4e9bJntbFx8OHLsFdOJl30u22bBEXGEhKRQW5wren2WBV+O+++25mz55NWVmZeVtZWRmzZ89m8uTJDdI5QevB1uBucbGbwu9GOWdKfgrd23RH64CohYVBbq5j7VpN8wB07SoLvwK2yrFtY0nOSna7HVt4qtWMDglhiyjrbDZYFf6///3vBAcHExUVxbBhw4iLiyMqKorQ0FD+85//NGQfBa0Aa4O7RqOc13d55m5pqSygLjZwuuA0PcJ7IEmSXfE35fkdIbGkhDhrVzN/f/lK5+ZKXACxkbENluffLPL8zQarwq9Wq1m5ciVHjhzhiSee4KmnnuLkyZN88MEHqF1KtAoE1rGW6ikru5w/dwnTAiwuVPPA5Yg/ICCgxt2vJZwS/uJi4qxF/ADR0ZCW5nhHrdAvsh+Hsw+73Y49hPA3L+wqeJcuXbj99tu59dZbadeuXUP0SdAKsSb8xcXy5CiXMQ3surDIOsgRf/c23fH391dM+At0OtK1Wq6wMmENkNM9Z8443E9rDIoaxIHMA263Y48r/PwoNxpJraio93MJ3EeE7oImgVarxcPDg4SEBIzGy6WTig/sOlHOaTAaSCtMIyY0hoCAAEpLS23u76jwJ12yafCwdeesUMTfNaQrZboyssuy3W7LFiqViglt2rDRFYtSQYMjhF/QJNBqtbz++utMnDiR2bNnm8VfkYjfNLDrZDlnekk6bXzb4Ofp51Cqx9GqnsSSEttpHlAs4lepVA0W9d8cFsYGIfzNAiH8giZBXl4ev//+O4mJiaxbt45t27YBcsQfEOCGQZsbzpyn80/TLbQbAP7+/g5F/I5U9SQWF1sf2DWhUMQPcrpnf+Z+RdqyxdjQUHYXF1MkVuhr8gjhFzQJvL29mTBhAnFxcUydOpXRo0cDsv9NSgr06QMzZsizd53CjVLO0wWn6dZGFn6lUj2SJLHHVimnCYUifoAr213JgYv1H/EHeHgwKjiYX5w2UxI0NEL4BU0CnU7Hiy++yMaNG1m1apW5ciwz87J7gUsGbW5M3krJTzFH/EpV9ZytrASgk73Z7126wPnzLlzp6jKoXcNE/AA3h4ezwdHJDIJGQwi/oEmg0+nw8fEhPj6+RrmwSnXZpsElg7b0dNi40aWKnhN5J+gd3htwPNVjT/i3FxUxMji4hhW5RXx85AuWAjN4e4X1Irssm7zy+s+/3xQWxs/5+eiMtr2NBI2LEH5Bk8DaBK7iYnj0UTcM2i5ckEUfnDZoO5573Cz8Sg3u7igqYlRIiGMd6NkTTri/fKJGrWFI+yEkpie63ZY92nt7093XV7h1NnGE8AuaBNbq+EtKZK97tRqcXupZp5OVOCrK6YoenUHHmYIzdG/THXAs4g8KgooKsDXBd3tREaOCgx3rRK9eigg/wLCOw9h9YbcibdlDVPc0fYTwC5oEtoTf5Tr+zEw5XaLROJ3jTy1IpWNQR3w85BXnHBncVankqN/awnG5Wi3pVVUMcNQ+olcvOKmMwdqwjsPYnd4wwn9bRATf5ORgVMBrSFA/COEXNAnqZebuhQvQoYNLh1ZP84BjqR6A9u3lYQVL7CgqYnhQEBp7+X0TCkf8ey7swSjVf+69r78/YR4ebBOmbU0WIfyCJoFOp7OY43cr4r9wQfa2d6Gip7bwO5LqAfk6Y034tzuT3wdFhT/SP5IwvzCO5x5XpD173NO2LZ9nZTXIuQTOI4Rf0OiYnC+tCb9bEb9pcNeNih5wPOLv2NGO8Dua3wfZlz8vT3YYVYCRnUey/ex2Rdqyx91t2/Jdbi6VYi3eJokQfkGjYzAYUKvVFpfhdMuL/8IFuOMO+bGTFT3Hco/RK6yX+bm7EX+pXs9fZWX2Z+xWR62G7t0Vy/PHd4kn4WyCIm3Zo4O3N1cGBPCTGORtkgjhFzQ69lbfctmywZTqcbKiR5IkjmYfpV9kP/M2RwZ3wbrw7y4uZmBAAD7O+ktfcQX89Zdzx1ghvms8CWkJSA006CrSPU0XIfyCRseW8BcVybrtkmVDerpLOf6zRWcJ8g4i1DfUvM3RVI814Xeqfr86AwbAoUPOH2eBriFd8dJ41ftSjCZuj4hgS2EheWKp1iaHEH5Bo2NtYBfkVM9338mPnbZsMEX8TnI46zCxbWNrbHMm1WOymKjOH4WFjHEmv29i4EA4eND54yygUqm4pus1bEnbokh79gjy8GBSWBifKrCSmEBZhPALGh1rEb8kybN1PTzk505ZNuj1skFb+/ZOR/zJWcn0j+xfY5uzEX/1bEqJXs+B0lJGuxPxK5SeuT7men49/asibTnCwx06sCI9HYOo6W9SCOEXNDq2ll308ZHF32nLhvPn5clbL7/sdFVPcnZynYjf0Ry/qQKpuPjytq2FhVwVGIifK+tHtm8vLzycmen8sRYY3208v5/5nSp9lSLt2WNYcDARnp5ikLeJIYRf0OjYmrUbFCRbNTht2ZCaCt26ubTy1uGsw8RGupbqUanqlnT+VlDA2NBQ6wfZa3DgQMXy/BH+EfQJ78OOczsUac8RHunYkbcs5b8EjYYQfkGjY62Gv6BA9ulxidOnISZGfuxEVU+lvpIzhWcI9g5mwucTuFgq56e9vb0xGAzoHBiorD3A+6s7wg9yukehPD/AxB4T2Xhqo2Lt2eOOiAhOlJdzWKH5CAL3EcIvaHR0Op3FiD872wUbZhOmiB+cyvEfzjpMz7CevLXnLRLTE1m2cxkgD4w6suA6yPOuTGuonK6ooECn40p3Fg6+6irYs8f142txY48b+eHEDw1W1umlVjO/fXsR9TchhPALGh1rqZ6cHIiMdLHR6hG/E+xN30tc+zj2ZuyloLKghpWxowO8vXtfdlr4MTeXSWFhqB3157HEiBGwc6diA7xXtrsSCalBVuUyMa99e77PzeXcpYVoBI2LEH5Bo2NL+N2O+J2s6NmbsZe4DnF8PeVrxncbzzd3fmN+zdEB3t695QlnAD/m5XFTeLhTfahDp07g56fYDF6VSsWdV9zJV0e/UqQ9Rwj38mJu+/a8rMDCMgL3EcIvaHSs5fizs2WbY6dn7UqSHPF//bXTFT2J6Ylc1eEqlu1cViPVA44P8JqEv1CnY29JCde7k983YYr6FeKufnfx1dGvGizdA/B4p058k5PD6YqKBjunwDJC+AWNjrWIPysLNm1yYdZuQYEs/kuXys8drOgprirmXNE5+kb0tZjqCQ4OptABq+GYGHlw95uLeVwTEoK/K2WctRk5EnYoV4kzoO0AfD192X6uYUzbAMI8PVnUqRP/TE1tsHMKLCOEX9DomGbuJiQkYKy2Vmt6OuzdKz92atbu6dNymkelcqqiZ1/GPgZEDcBT48mO+3YgPSux477LYtuxY0cuODBA6ekpi/8n57K5y+VBilqMGgUJCYrl+VUqFXOvnMvKfSsVac9RHu3YkT3FxWwXXv2NihB+QaNTWVlJcnIyEydOZPbs2WbxLy6WF88CJ2ftJidD375O92P7ue1c3fFqMksya5RymujSpQvnzp1zqK3ogTr2VRVxc1iY0/2wSL9+8lKSCvnzA0wfMJ2fTv7UIIuwm/DVaHi9WzfmnTyJVizI3mgI4Rc0OocPH+b8+fMkJiaybt06tm3bBsiDu4mJLszaPXQI+veX0ztO5Pj/OPMH10ZfazG/D9C5c2fOOjg4aRyRQ3ReGwJMfhPuolLBpEnw00/KtAe08W3DLb1v4b/7/qtYm44wOSKCaB8fljp4ERUojxB+QaMTExNDdHQ0cXFxTJ06ldGjRwPy4G6HDi7M2j18WJ705MSs3QpdBUkZSYzsPNJifh+ci/hP9sgkeFeUE512gBtvVFT4AZ64+gne3vM2FbqGG3BVqVS837MnK9LTSarubSFoMITwCxodvV5PfHw8GzduZNWqVajVagwGeYzW6UyJJMkR/4AB8nMHc/y7zu9iQNQASrWlBHgFkPlYZo38PsgRvyPCf6i0lEo/LSfXtlEqJS9z7bWwfz/k5irWZN/IvgztOJTVB1Yr1qYjdPLxYUWPHvzt2DGK9PoGPbdACL+gCaDT6fD29iY+Ph61Wv6TzM+XfXpSUpwcz7xwAby8ZIM2cLiM8/czv3NtV+tpHrgs/PZKIP+bkcGcDlEE+qnM9fyK4OcHN90Ea9cq2CgsHrOYl7a/RHFVw0bfd0VGcn1oKPceO4ZRuHc2KEL4BY1OcXEx/v7+NbZlZsq5fadLOatH+05M3vr++PdM7DHRapoH5Dp+Pz8/cm1E3DlaLeuys3mgfXtGjYLtSldL3ncfrF6tWHUPwKB2g5jQfQIvb39ZsTYd5c3u3SnU63n41KkGnVPQ2hHCL2h0zp8/T+fOnWtsO336srWxU6WcBw5czu87OLB7NPsoJdoShnYcanHGbnXspXveSU9nSkQE7by9GT0aLo1TK8c118gfzL59ijb70rUv8dGBjzh48aCi7drDS61mQ2wse4qLecpkcCSodxpF+LOzs+nUqRPHFb0PFjRXzp07V0f49XpZs8HJUs7ff4f4eKcGdr/+62vu6HMHapXaZqoHbFf25Ot0vJeRweOdOgFyN/74w8nlIu2hVsOCBbBkiYKNQvvA9rw27jXuXX8vlfqG9dMJ9vBg04AB/Jiby3/OnBGRfwPQ4MKv0+l44IEH8PX1behTC5ooloT/zBl46CEnSzlLSuRIeMwYWewdGNiVJIl1R9Yxpe8UMksyWXNojdVUD8iVPdaE/7m0NKZERNDDzw+A7t3lqqTNmx3sv6PMny/bNyho1Qwwvf90+kb2Ze6PcxtcfMM8Pdk8YACb8vO5+6+/KFf0aimoTYML/+OPP868efNo3769xddXrlzJkCFDGDJkCDk5OQ3cO0FjYEn49+2T3YhPnry89KJdEhIgLg5efVVO8zz3nN1of9PpTXh7eDO843CW7VyGChWLhi2qU9FjYsCAAezevbvO9uTSUtZmZ/Nc1641ts+aBR9/7GD/HcXPD/71L3jsMXl1LoVQqVSsvnk1R3OO8syWZxpc/KO8vdk6cCBeajUjDxzguANOqAIXkRqQjz/+WHrhhRckSZKkMWPGSMeOHbO5/+DBgxuiW4JGpLCwUAoICJCMRmON7TExknTzzZIEkjR9uiTp9Q409tBDkrR0qfzYwT/t6z65Tvr04KeSJEnSiNUjJBYjjVg9wur+Fy9elIKDg6XKykrztgq9XuqXmCitzsios39eniQFB0tSdrZD3XEcnU6S4uIk6d13FW5YkrJKs6R+7/WTHt/0uGQwGhRv3x5Go1F658IFKWz7dunZ1FSpwqEvX1Ade9rZoMI/atQoafTo0dKYMWOk4OBg6aqrrpIyMzOt7i+Ev+Vz+PBh6YorrpAMBoO0ZcsWyWAwSPn5kuTnJ2u36ceucJaXS1JUlCSZgokxY+yeO+FMgtTp9U5Slb5KkiRJyijOkMZ/Nl7KLLH+NylJkjR8+HDpl19+kSRJFqmZx45Jdxw5UufiZeKRRyTp/vvtdsd5jh2TpIgISUpIULzpnLIcafTHo6Ubv7hRulhyUfH2HeF8RYV0a3Ky1PXPP6X3L1wQFwAnaFLCXx0R8QskSZJ++uknafz48dKsWbMkX19fadasWdKvvxqkESMkycNDFn0PDwcaev99SZo0SX787LOXrxjPPmtx9zJtmdTtrW7SD8d/MG/7+89/l0KXhEqLfllk81RLliyR5s+fLxmMRmnRqVPS0KQkqdSGKBUWSlK7dvWiz5K0ebMs/ps3K950lb5K+udv/5QiX42UVuxZYb5ANjQ7CguliYcOSe137pSeOn1aOlZa2ij9aE4I4Rc0ad577z1p0qRJkq+vr5ScnCz5+vpK99+/RYqJkXX7nnskqaLCTiNlZXJuaMeOy9tspHq0eq00+cvJ0rRvp5m3ZRRnSCFLQuymeiRJktLS0qQ2MTHSNTt2SCP27ZNytVq77/OXX2R93rPH7q7O88cf8t3OU09JUkmJ4s0fzDwoTfh8ghT5aqT0xK9PSMdzjit+Dkc4XFIiPXbqlNRu505pQGKi9ERKivRbXp5UpNM1Sn+aMva0UyVJTbd2asiQISQlJbl0rFarZ/78DVx1VSR79mQxdGhb/vwzE5VKxbBhUTa3Obu/Em201v1XrPgRDw8P/PyySErKZdy4viQkvEh5+eWlCm2uvavVwu23y94On3wibzMN6FoY2E3JT2HeT/Pw9fTlmynf4O3hDcCiXxbxyaFPmDlwJq+Pf93q39W5yko+yszk9dRUNJs3s+nOO4kbNMihv8kff5TnX913n1yReanqUxkyMuCJJ+QFDO65ByZMkD38a02Mc4cTuSdYtX8Vnyd/jo+HD9d2vZZhHYfRN7IvfcL7EOITgsqdJSYdxCBJ/FlUxG8FBWwuKOBgaSntvLwYGBDAgIAAYnx96eztTWcfH9p7eeGpbn3TlexpZ4sU/vJyPf7rtkJ0tS9cUsHXHeGD7gr2UFDfqNUGnvz8DSQAyYgRQJKQjEakggL46y+k9u2RJt6ApFYz7oN1TPhwHQA/z7mLnx+YSrm2nOzyXFLyT3G26Dw397qJ8d0noFHJns8FZVn866vRlGuL6dH2Kp665QckQAKK9XoytVrOV1VxoLSUXJ2OOyIieKRDB3auW8fTTz9Nz549GTx4MO3atSMyMhJvb288PDzw8PBAo9Hg4eFhFsS8PG/WrevG1q3tCAuromvXEsLDK2nTpgofHwNeXgY8PY14eBhRqS7PZbisp1KtbZL5sUoFofkXiD34M11Tk4jKPE5hSDuKQttTHNyWSp9AtN5+VHn7o/PyxajSIKnVSCo1xku/JbUGo0p+bA0JiXTPLI74nuSM93nSvbLI8MrCgJEQQxCh+iCCjIH4GL3lH8kLH6M3XpIXGkmNBjWqS781kho1atSSGg0aVKhQSSCh4vIlpPojVZ0tRpWKwoBAckJCyQkJodjfnyJ/f0r8/Cn18cFLr8dHq8VHq8VXq8Vbq8XTYERjNOBhMOBhMKIxGPC49FxjlFBJEiou/Za49PvSNqPpNVBf2l4dFdWe11LXmpdFy/uZ9ukf2Yl/z59l9XuwRasU/ief/IRlEzoDRhyvWK32lRiBb8RFoikQ/t4HaELLLn07l0TO9KJGg8rHGzy9zNsf/XwDj37xI51+/lDeppLwVHvi5+lHG99Qovyj8FBrzPurVCp2732W02nr6R59G8OHPGduX6VSEajR0M7Liw7e3vT39yc2IABNtai2vLyc7du3c+TIEbKyssjKykKr1aLX69Hr9RgMBnQ6XZ33ZTRqKCmJoaysI5WV4VRVtcFg8MZo9MRo9MJo9ESSzD2pduSlT0Kqua3mcxlfqZIYwzk6GzPoYLhIoFRGoFRGgFSOHxVoJCNqjGgwosaARjKiwYAGI3UUywHKvQzkB2jJD6iiyE9PhZeBCi8DlZ5GKi/9NqgkDGoJo1rCoAaDyvT40naVLJymt1O9F7XfoqSSLO6jqrbNqAKDlx9670AMXv7ovQPRe/tj1HghaTwxarwwajyRNF4YPbzk7WqN/JmafqtUYLoYqtS1nqtqXSRt3PFU+7uRqu9X45DLT3qmXGTvYssTCe3RKoVfr9fj+b+PIUgJ4bbyRUqIO4h6RqORZ/A6jMmmAeTJWw549WSWZBL9VjRVhiraB7Yn/dF0V7oqEDQp7GmnQqtENC08PDwoumYWBQN6M+3a0+zsDAx6HwJ7udCaleuiCrjzvPxjhWgdTNt+CJVKRZ8+/Th1KpkePWL566/kOttsvdYa9o+NHcjYsaPQatX4+MjOBN7eTn5VCQlO7Z5Zksnw1cPRG+WrS3RItJMnFAiaJy0y4q/B4sWyIMTHy7+3bq3x8sj7kC8MJmLmQcc7cCxFJIGNXKg1NMDfO3bkte7ibkEx4uMvf7cORvtzf5zLh/vllNCITiOsztYVCJobrTLVY5P4eKsXgdoE/RNKfCy8IAFXunoHYRmNJPH39HReW736cj+rX7AsbbP1Wmvaf80aMPnnjBnjUOQvUjyClkyzFv7w8HC61vI+cZScnBwi7Fk6ZmTIxl4ApaV22zweBmVe9s/tYQSDTwSSd6gDPbWNX2UlfcTapA5THhyMn407KZ1Rx+n805TrypGQQAJ/L396h/duwF7WHw793bdQxHu//N7T0tJsrhvRpIXfHerlbqGZIN5763zv0Lrfv3jvjr/31jezQSAQCFo5QvgFAoGgldFihX/u3LmN3YVGQ7z31ktrfv/ivTtOi83xCwQCgcAyLTbiFwgEAoFlhPALBAJBK6PFCb/RaGTevHkMHz6c+Ph4UlJSGrtLDcqgQYOIj48nPj6eWbNcc/ZrbuzZs4f4+HgAUlJSGDlyJKNGjWL+/PkYFVyTtilS/b3v37+fDh06mL//L7/8snE7V0/odDqmT5/OqFGjiIuLY8OGDa3me7f03l363utlFYBG5Ntvv5VmzJghSZIk/fnnn9LNN9/cuB1qQCoqKqSBAwc2djcalKVLl0r9+vWThg4dKkmSJN10003Sli1bJEmSpAceeED67rvvGrF39Uvt9/7hhx9Kr732WiP3qv756KOPpEceeUSSJEnKzc2VOnXq1Gq+d0vv3ZXvvcVF/Dt27GDChAkADBs2rFVN6Dh06BDl5eWMGzeOa6+9lt27dzd2l+qdbt268d1335mf79u3jzFjxgBwww03sHnz5sbqWr1j6b3/73//Y/To0dx///2UmGaltzCmTJnCCy+8YH7u4eHRar53a+/d2e+9xQl/cXExwcHB5ucajQa9U96+zRc/Pz8ef/xxNm3axAcffMC0adNa/HufPHkynp6e5ueSJJkXPQkMDKSoqKixulbv1H7vcXFxvPrqq2zbto2YmBieM1lUtzACAgIIDAykpKSEO+64gxdffLHVfO+W3rsr33uLE/6goKAaVzyj0YiHR4t0n65Dz549ueeee1CpVPTs2ZOwsDAyMzMbu1sNirraMnslJSWEhIQ0XmcamNtuu43BgwebHx84cKCRe1R/nD9/nmuuuYbp06fzt7/9rVV977Xfuyvfe4sT/hEjRrBx40YAdu/eTWxsbCP3qOH46KOPeOyxxwDIyMiguLiYdu3aNXKvGpZBgwaRcMmd8+eff2bUqFGN26EGZPz48SQmJgLw+++/m8WgpZGVlcW4ceNYunQp9913H9B6vndL792V773FTeAyGo08+OCDHD58GEmS+Pjjj+ndu2U4L9pDq9Uyc+ZMzp07h0qlYunSpVx99dWN3a16Jy0tjalTp7J7925OnjzJnDlz0Gq19OnThw8//BCNRtPYXaw3qr/3/fv3s2DBAry8vIiKimLlypUEBQU1dhcV55FHHuHLL7+s8X/91ltvsXDhwhb/vVt67y+99BL/+Mc/nPreW5zwCwQCgcA2LS7VIxAIBALbCOEXCASCVoYQfoFAIGhlCOEXCASCVoYQfoFAIGhlCOEXNFsWLlzIwIEDGThwIF5eXvTq1cv8/P3332fJkiX13oeEhAT69etX7+cRCJSkdUxpFbRI3n77bfPjrl278sUXXzBkyJBG7JFA0DwQEb+gRbJ48WIWLFgAyBeFp556imHDhtGrVy9Wr17Nfffdx4ABAxgyZAgZGRkApKenm6e/9+/fn5dfftmhc5WWljJ16lQGDhxI79692b59OwBFRUXcc8899OvXj9jYWJ544gmzd5JKpSI3N9fchul5QkICAwYM4Oqrr6Z///4UFRUxZcoUBg4cyJVXXsmcOXNarOWwoOEQwi9oFVRWVrJ7926ef/555s6dyyOPPMKhQ4fo1KkTa9asAWD69Oncd9997Nu3j8TERDZv3sxXX31lt+0LFy6waNEiDh48yAMPPMDixYsBORUVFhZGcnIySUlJHDp0iNdee81ue0eOHOH//u//OHz4MBs2bKCkpISDBw+yd+9eAFJTU13+HAQCEMIvaCVMnjwZkK2Mo6KiGDBggPl5fn4+ZWVlbN26lWeeeYaBAwcybNgwzp07x8GDB+223a1bN4YOHQrAwIEDyc7OBmTPmAULFqBSqfD29mbevHn8/PPPdtvr1KkTXbp0AWDkyJEcPXqU+Ph4lixZwt///ne6d+/uykcgEJgROX5Bq8Db29v8uLqVsQmDwYAkSezatQs/Pz8AcnNz8fHxsdt29fZUKhUmFxSj0Wi2CjY91+l05uem/bRabY32AgICzI+jo6NJSUkhISGBP/74g+uvv56VK1dy00032e2XQGANEfELBMh23sOGDeP1118HoLCwkBEjRvDDDz+43Ob48eN55513kCSJqqoqVq5cydixYwGIiIgwLxK0du1aq228//77zJo1y+zIOH78ePbv3+9ynwQCEMIvEJhZu3at2cp76NCh3H333UybNg2A2bNn88EHHzjV3ttvv012djaxsbHExsbSq1cvnn76afNrDz30EFdeeSXHjh2zap997733YjAYuOKKKxg8eDBFRUUsXLjQvTcqaPUId06BwAF+++03Tp8+zbx58xq7KwKB24iIXyBwgLy8PHP0LxA0d0TELxAIBK0MEfELBAJBK0MIv0AgELQyhPALBAJBK0MIv0AgELQyhPALBAJBK+P/AYgrp4SAM8O8AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# set up figure\n",
    "try:\n",
    "    plt.close(\"all\")\n",
    "except:\n",
    "    pass\n",
    "\n",
    "set_plot_params()\n",
    "\n",
    "fig = plt.figure(figsize=(6, 4), facecolor=\"w\")\n",
    "ax = fig.add_subplot(1, 1, 1)\n",
    "\n",
    "ax.plot(ts1_Otis[:, 0], ts1_Otis[:, 1], \"k-\", linewidth=1.0)\n",
    "ax.plot(ts2_Otis[:, 0], ts2_Otis[:, 1], \"b-\", linewidth=1.0)\n",
    "ax.plot(ts3_Otis[:, 0], ts3_Otis[:, 1], \"r-\", linewidth=1.0)\n",
    "ax.plot(ts4_Otis[:, 0], ts4_Otis[:, 1], \"g-\", linewidth=1.0)\n",
    "ax.plot(ts5_Otis[:, 0], ts5_Otis[:, 1], \"c-\", linewidth=1.0)\n",
    "\n",
    "ax.plot(\n",
    "    (ts1_mt3d[:, 0]) / 3600,\n",
    "    ts1_mt3d[:, 1],\n",
    "    \"kD\",\n",
    "    markersize=2.0,\n",
    "    mfc=\"none\",\n",
    "    mec=\"k\",\n",
    ")\n",
    "ax.plot(\n",
    "    (ts2_mt3d[:, 0]) / 3600,\n",
    "    ts2_mt3d[:, 1],\n",
    "    \"b*\",\n",
    "    markersize=3.0,\n",
    "    mfc=\"none\",\n",
    "    mec=\"b\",\n",
    ")\n",
    "ax.plot((ts3_mt3d[:, 0]) / 3600, ts3_mt3d[:, 1], \"r+\", markersize=3.0)\n",
    "ax.plot(\n",
    "    (ts4_mt3d[:, 0]) / 3600,\n",
    "    ts4_mt3d[:, 1],\n",
    "    \"g^\",\n",
    "    markersize=2.0,\n",
    "    mfc=\"none\",\n",
    "    mec=\"g\",\n",
    ")\n",
    "ax.plot(\n",
    "    (ts5_mt3d[:, 0]) / 3600,\n",
    "    ts5_mt3d[:, 1],\n",
    "    \"co\",\n",
    "    markersize=2.0,\n",
    "    mfc=\"none\",\n",
    "    mec=\"c\",\n",
    ")\n",
    "\n",
    "# customize plot\n",
    "ax.set_xlabel(\"Time, hours\")\n",
    "ax.set_ylabel(\"Concentration, mg L-1\")\n",
    "ax.set_ylim([3.5, 13])\n",
    "ticksize = 10\n",
    "\n",
    "# legend\n",
    "leg = ax.legend(\n",
    "    (\n",
    "        \"Otis, Site 1\",\n",
    "        \"Otis, Site 2\",\n",
    "        \"Otis, Site 3\",\n",
    "        \"Otis, Site 4\",\n",
    "        \"Otis, Site 5\",\n",
    "        \"MT3D-USGS, Site 1\",\n",
    "        \"MT3D-USGS, Site 2\",\n",
    "        \"MT3D-USGS, Site 3\",\n",
    "        \"MT3D-USGS, Site 4\",\n",
    "        \"MT3D-USGS, Site 5\",\n",
    "    ),\n",
    "    loc=\"upper right\",\n",
    "    labelspacing=0.25,\n",
    "    columnspacing=1,\n",
    "    handletextpad=0.5,\n",
    "    handlelength=2.0,\n",
    "    numpoints=1,\n",
    ")\n",
    "leg._drawFrame = False\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "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.8.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}