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": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3.6.0 |Anaconda 4.3.0 (x86_64)| (default, Dec 23 2016, 13:19:00) \n",
      "[GCC 4.2.1 Compatible Apple LLVM 6.0 (clang-600.0.57)]\n",
      "numpy version: 1.11.3\n",
      "matplotlib version: 2.0.0\n",
      "flopy version: 3.2.6\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "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",
    "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": {
    "collapsed": false
   },
   "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.exists(modelpth):\n",
    "    os.mkdir(modelpth)\n",
    "\n",
    "# Instantiate MODFLOW object in flopy\n",
    "mf = flopy.modflow.Modflow(modelname=modelname, exe_name=mfexe, model_ws=modelpth, version='mfnwt')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "Set up model discretization"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "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(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": {
    "collapsed": true
   },
   "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": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Newton-Rhapson 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(mf, headtol=headtol, fluxtol=fluxtol, maxiterout=maxiterout,       \n",
    "                               thickfact=thickfact, linmeth=linmeth, iprnwt=iprnwt, 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": {
    "collapsed": true
   },
   "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",
    "def topElev(X, Y):\n",
    "    return 100. - (np.ceil(X)-1) * 0.03\n",
    "\n",
    "grndElev = topElev(X, Y)\n",
    "bedRockElev = grndElev - 3.\n",
    "\n",
    "Steady = [False, False, False]                                                                  \n",
    "nstp = [1, 1, 1]                                                                        \n",
    "tsmult = [1., 1., 1.]\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(mf, nlay, nrow, ncol, nper=3, delr=delr, delc=delc,     \n",
    "                               top=grndElev, botm=bedRockElev, laycbd=0, itmuni=1, lenuni=2,\n",
    "                               steady=Steady, nstp=nstp, tsmult=tsmult, perlen=perlen) \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate upstream weighting (UPW) flow package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "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.\n",
    "ss = 0.000001\n",
    "sy = 0.20\n",
    "hdry = -888\n",
    "\n",
    "upw = flopy.modflow.ModflowUpw(mf, laytyp=laytyp, layavg=layavg, chani=chani, layvka=layvka,\n",
    "                               ipakcb=53, hdry=hdry, iphdry=iphdry, hk=hk, hani=hani,\n",
    "                               vka=vka, ss=ss, sy=sy)\n",
    "                                 \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate basic (BAS or BA6) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "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",
    "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)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate streamflow routing (SFR2) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "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('temp'):\n",
    "    os.mkdir('temp')\n",
    "\n",
    "fpth = os.path.join('temp', 's1.csv')\n",
    "f = open(fpth, 'w')\n",
    "f.write(s1)\n",
    "f.close()\n",
    "\n",
    "dtype = [('k', '<i4'), ('i', '<i4'), ('j', '<i4'), ('iseg', '<i4'), \n",
    "         ('ireach', '<f8'), ('rchlen', '<f8'), ('strtop', '<f8'), \n",
    "         ('slope', '<f8'), ('strthick', '<f8'), ('strhc1', '<f8')]\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('temp', '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,                                                                                                                                                             \n",
    "                1: segment_data,\n",
    "                2: segment_data}\n",
    "                                                                                                                                                                                             \n",
    "# There are 3 stress periods\n",
    "dataset_5 = {0: [nss, 0, 0],\n",
    "             1: [nss, 0, 0],\n",
    "             2: [nss, 0, 0]}\n",
    " \n",
    "sfr = flopy.modflow.ModflowSfr2(mf, nstrm=nstrm, nss=nss, const=const, dleak=dleak, isfropt=isfropt, istcb2=0,\n",
    "                                reachinput=True, reach_data=reach_data, dataset_5=dataset_5,\n",
    "                                segment_data=segment_data, 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": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "gages = [[1,38,61,1],[2,67,62,1], [3,176,63,1], [4,152,64,1], [5,186,65,1], [6,31,66,1]]\n",
    "files = ['CrnkNic.gage','CrnkNic.gag1','CrnkNic.gag2','CrnkNic.gag3','CrnkNic.gag4','CrnkNic.gag5',\n",
    "         'CrnkNic.gag6']\n",
    "gage = flopy.modflow.ModflowGage(mf, numgage=6, gage_data=gages, filenames = files)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate streamflow routing (SFR2) package for MODFLOW-NWT"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "lmt = flopy.modflow.ModflowLmt(mf, output_file_name='CrnkNic.ftl', output_file_header='extended',\n",
    "                               output_file_format='formatted', package_flows = ['sfr'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Write the MODFLOW input files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/Users/JosephHughes/Development/flopy-git/examples/Notebooks\n",
      "FloPy is using the following executable to run the model: /Users/JosephHughes/USGS/Development/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.1.2, 9/15/2016                        \n",
      "                    BASED ON MODFLOW-2005 Version 1.11.0 08/08/2013                       \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): 2017/03/18 19:57:55\n",
      "\n",
      " Solving:  Stress period:     1    Time step:     1    Groundwater-Flow Eqn.\n",
      " 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): 2017/03/18 19:57:55\n",
      " Elapsed run time:  0.134 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": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Instantiate MT3D-USGS object in flopy\n",
    "mt = flopy.mt3d.Mt3dms(modflowmodel=mf, modelname=modelname, model_ws=modelpth,\n",
    "                       version='mt3d-usgs', namefile_ext='mtnam', exe_name=mtexe,\n",
    "                       ftlfilename='CrnkNic.ftl', ftlfree=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate basic transport (BTN) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/JosephHughes/Development/flopy-git/flopy/mt3d/mtbtn.py:492: UserWarning: BTN warning. Laycon has not been set.  A modflow model with a  BCF or LPF package does not exist and laycon was not passed  to the BTN constructor.  Setting laycon to 1 (convertible).\n",
      "  warnings.warn(s)\n"
     ]
    }
   ],
   "source": [
    "btn = flopy.mt3d.Mt3dBtn(mt, sconc=3.7, ncomp=1, prsity=0.2, cinact=-1.0, \n",
    "                         thkmin=0.001, nprs=-1, nprobs=10, chkmas=True, \n",
    "                         nprmas=10, dt0=180, mxstrn=2500)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate advection (ADV) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "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": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "gcg = flopy.mt3d.Mt3dGcg(mt, mxiter=10, iter1=50, isolve=3, ncrs=0,\n",
    "                   accl=1, cclose=1e-6, iprgcg=1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate source-sink mixing (SSM) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "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)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instantiate streamflow transport (SFT) package for MT3D-USGS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "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],\n",
    "                         1: [0, 0, 11.4],\n",
    "                         2: [0, 0, 3.7]}\n",
    "\n",
    "gage_output = [None, None, 'CrnkNic.sftobs']\n",
    "\n",
    "sft = flopy.mt3d.Mt3dSft(mt, nsfinit=650, mxsfbc=650, icbcsf=81, ioutobs=82,\n",
    "                         isfsolv=1, cclosesf=1.0E-6, mxitersf=10, crntsf=1.0, iprtxmd=0,\n",
    "                         coldsf=3.7, dispsf=dispsf, nobssf=5, obs_sf=obs_sf,\n",
    "                         sf_stress_period_data = sf_stress_period_data,\n",
    "                         filenames=gage_output)\n",
    "\n",
    "sft.dispsf.format.fortran = \"(10E15.6)\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Write the MT3D-USGS input files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "FloPy is using the following executable to run the model: /Users/JosephHughes/USGS/Development/bin/mt3dusgs\n",
      "\n",
      " MT3D-USGS - Modular 3D Multi-Species Transport Model [Ver 1.00.00] \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.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:    2   Step Size:   180.0     Total Elapsed Time:   360.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:    3   Step Size:   180.0     Total Elapsed Time:   540.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:    4   Step Size:   180.0     Total Elapsed Time:   720.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:    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",
      " 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",
      " 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",
      " 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",
      " 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.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:    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",
      " 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",
      " Transport Step:    1   Step Size:   180.0     Total Elapsed Time:   41580.    \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:    2   Step Size:   180.0     Total Elapsed Time:   41760.    \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:   41940.    \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:   42120.    \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:   42300.    \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:   42480.    \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:   42660.    \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:   42840.    \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:   43020.    \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:   43200.    \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:   43380.    \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:   43560.    \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:   43740.    \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:   43920.    \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:   44100.    \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:   44280.    \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:   44460.    \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:   44640.    \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:   44820.    \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:   45000.    \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:   45180.    \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:   45360.    \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:   45540.    \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:   45720.    \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:   45900.    \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:   46080.    \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:   46260.    \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:   46440.    \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:   46620.    \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:   46800.    \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:   46980.    \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:   47160.    \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:   47340.    \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:   47520.    \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:   47700.    \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:   47880.    \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:   48060.    \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:   48240.    \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:   48420.    \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:   48600.    \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:   48780.    \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:   48960.    \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:   49140.    \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:   49320.    \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:   49500.    \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:   49680.    \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:   49860.    \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:   50040.    \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:   50220.    \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:   50400.    \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:   50580.    \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:   50760.    \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:   50940.    \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:   51120.    \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:   51300.    \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:   51480.    \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:   51660.    \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:   51840.    \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:   52020.    \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:   52200.    \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:   52380.    \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:   52560.    \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:   63   Step Size:   180.0     Total Elapsed Time:   52740.    \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:   52920.    \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:   53100.    \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:   53280.    \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:   53460.    \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:   53640.    \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:   53820.    \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:   54000.    \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:   54180.    \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:   54360.    \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:   54540.    \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:   54720.    \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:   54900.    \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:   55080.    \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:   55260.    \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:   55440.    \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:   55620.    \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:   55800.    \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:   55980.    \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:   56160.    \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:   56340.    \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:   56520.    \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:   56700.    \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:   56880.    \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:   57060.    \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:   57240.    \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:   57420.    \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:   57600.    \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:   57780.    \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:   57960.    \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:   58140.    \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:   58320.    \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:   58500.    \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:   58680.    \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:   58860.    \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:   59040.    \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:   59220.    \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:   59400.    \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:   59580.    \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:   59760.    \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:   59940.    \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:   60120.    \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:   60300.    \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:   60480.    \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:   60660.    \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:   60840.    \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:   61020.    \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:   61200.    \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:   61380.    \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:   61560.    \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:   61740.    \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:   61920.    \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:   62100.    \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:   62280.    \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:   62460.    \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:   62640.    \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:   62820.    \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:   63000.    \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:   63180.    \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:   63360.    \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:  123   Step Size:   180.0     Total Elapsed Time:   63540.    \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:   63720.    \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:   63900.    \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:   64080.    \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:   64260.    \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:   64440.    \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:  129   Step Size:   180.0     Total Elapsed Time:   64620.    \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:   64800.    \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:   64980.    \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:   65160.    \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:   65340.    \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:   65520.    \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:   65700.    \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:   65880.    \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:   66060.    \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:   66240.    \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:   66420.    \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:   66600.    \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:   66780.    \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:   66960.    \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:   67140.    \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:   67320.    \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:   67500.    \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:   67680.    \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:   67860.    \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:   68040.    \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:   68220.    \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:   68400.    \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:   68580.    \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:   68760.    \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:   68940.    \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:   69120.    \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:   69300.    \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:   69480.    \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:   69660.    \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:   69840.    \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:   70020.    \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:   70200.    \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:   70380.    \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:   70560.    \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:   70740.    \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:   70920.    \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:   71100.    \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:   71280.    \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:   71460.    \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:   71640.    \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:   71820.    \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:   72000.    \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:  171   Step Size:   180.0     Total Elapsed Time:   72180.    \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:  172   Step Size:   180.0     Total Elapsed Time:   72360.    \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:  173   Step Size:   180.0     Total Elapsed Time:   72540.    \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:  174   Step Size:   180.0     Total Elapsed Time:   72720.    \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:  175   Step Size:   180.0     Total Elapsed Time:   72900.    \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:  176   Step Size:   180.0     Total Elapsed Time:   73080.    \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:  177   Step Size:   180.0     Total Elapsed Time:   73260.    \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:  178   Step Size:   180.0     Total Elapsed Time:   73440.    \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:  179   Step Size:   180.0     Total Elapsed Time:   73620.    \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:  180   Step Size:   180.0     Total Elapsed Time:   73800.    \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:  181   Step Size:   180.0     Total Elapsed Time:   73980.    \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:  182   Step Size:   180.0     Total Elapsed Time:   74160.    \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:  183   Step Size:   180.0     Total Elapsed Time:   74340.    \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:  184   Step Size:   180.0     Total Elapsed Time:   74520.    \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:  185   Step Size:   180.0     Total Elapsed Time:   74700.    \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:  186   Step Size:   180.0     Total Elapsed Time:   74880.    \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:  187   Step Size:   180.0     Total Elapsed Time:   75060.    \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:  188   Step Size:   180.0     Total Elapsed Time:   75240.    \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:  189   Step Size:   180.0     Total Elapsed Time:   75420.    \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:  190   Step Size:   180.0     Total Elapsed Time:   75600.    \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:  191   Step Size:   180.0     Total Elapsed Time:   75780.    \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:  192   Step Size:   180.0     Total Elapsed Time:   75960.    \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:  193   Step Size:   180.0     Total Elapsed Time:   76140.    \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:  194   Step Size:   180.0     Total Elapsed Time:   76320.    \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:  195   Step Size:   180.0     Total Elapsed Time:   76500.    \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:  196   Step Size:   180.0     Total Elapsed Time:   76680.    \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:  197   Step Size:   180.0     Total Elapsed Time:   76860.    \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:  198   Step Size:   180.0     Total Elapsed Time:   77040.    \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:  199   Step Size:   180.0     Total Elapsed Time:   77220.    \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:  200   Step Size:   180.0     Total Elapsed Time:   77400.    \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:  201   Step Size:   180.0     Total Elapsed Time:   77580.    \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:  202   Step Size:   180.0     Total Elapsed Time:   77760.    \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:  203   Step Size:   180.0     Total Elapsed Time:   77940.    \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:  204   Step Size:   180.0     Total Elapsed Time:   78120.    \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:  205   Step Size:   180.0     Total Elapsed Time:   78300.    \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:  206   Step Size:   180.0     Total Elapsed Time:   78480.    \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:  207   Step Size:   180.0     Total Elapsed Time:   78660.    \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:  208   Step Size:   180.0     Total Elapsed Time:   78840.    \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:  209   Step Size:   180.0     Total Elapsed Time:   79020.    \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:  210   Step Size:   180.0     Total Elapsed Time:   79200.    \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:  211   Step Size:   180.0     Total Elapsed Time:   79380.    \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:  212   Step Size:   180.0     Total Elapsed Time:   79560.    \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:  213   Step Size:   180.0     Total Elapsed Time:   79740.    \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:  214   Step Size:   180.0     Total Elapsed Time:   79920.    \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:  215   Step Size:   180.0     Total Elapsed Time:   80100.    \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:  216   Step Size:   180.0     Total Elapsed Time:   80280.    \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:  217   Step Size:   180.0     Total Elapsed Time:   80460.    \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:  218   Step Size:   180.0     Total Elapsed Time:   80640.    \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:  219   Step Size:   180.0     Total Elapsed Time:   80820.    \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:  220   Step Size:   180.0     Total Elapsed Time:   81000.    \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:  221   Step Size:   180.0     Total Elapsed Time:   81180.    \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:  222   Step Size:   180.0     Total Elapsed Time:   81360.    \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:  223   Step Size:   180.0     Total Elapsed Time:   81540.    \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:  224   Step Size:   180.0     Total Elapsed Time:   81720.    \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:  225   Step Size:   180.0     Total Elapsed Time:   81900.    \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:  226   Step Size:   180.0     Total Elapsed Time:   82080.    \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:  227   Step Size:   180.0     Total Elapsed Time:   82260.    \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:  228   Step Size:   180.0     Total Elapsed Time:   82440.    \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:  229   Step Size:   180.0     Total Elapsed Time:   82620.    \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:  230   Step Size:   180.0     Total Elapsed Time:   82800.    \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:  231   Step Size:   180.0     Total Elapsed Time:   82980.    \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:  232   Step Size:   180.0     Total Elapsed Time:   83160.    \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:  233   Step Size:   180.0     Total Elapsed Time:   83340.    \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:  234   Step Size:   180.0     Total Elapsed Time:   83520.    \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:  235   Step Size:   180.0     Total Elapsed Time:   83700.    \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:  236   Step Size:   180.0     Total Elapsed Time:   83880.    \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:  237   Step Size:   180.0     Total Elapsed Time:   84060.    \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:  238   Step Size:   180.0     Total Elapsed Time:   84240.    \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:  239   Step Size:   180.0     Total Elapsed Time:   84420.    \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:  240   Step Size:   180.0     Total Elapsed Time:   84600.    \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:  241   Step Size:   180.0     Total Elapsed Time:   84780.    \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:  242   Step Size:   180.0     Total Elapsed Time:   84960.    \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:  243   Step Size:   180.0     Total Elapsed Time:   85140.    \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:  244   Step Size:   180.0     Total Elapsed Time:   85320.    \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:  245   Step Size:   180.0     Total Elapsed Time:   85500.    \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:  246   Step Size:   180.0     Total Elapsed Time:   85680.    \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:  247   Step Size:   180.0     Total Elapsed Time:   85860.    \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:  248   Step Size:   180.0     Total Elapsed Time:   86040.    \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:  249   Step Size:   180.0     Total Elapsed Time:   86220.    \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:  250   Step Size:   180.0     Total Elapsed Time:   86400.    \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",
      " Program completed.   Total CPU time:  000 minutes  1.574 seconds\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(False, [])"
      ]
     },
     "execution_count": 19,
     "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": 20,
   "metadata": {
    "collapsed": true
   },
   "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",
    "# 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\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load output from SFT as well as from the OTIS solution"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "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('..', 'data', 'mt3d_test', 'mfnwt_mt3dusgs', 'sft_crnkNic', 'OTIS_solution.out')\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": 22,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def set_plot_params():\n",
    "    import matplotlib as mpl\n",
    "    from matplotlib.font_manager import FontProperties\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",
    "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",
    "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\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compare output:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEJCAYAAACdePCvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd4VNXWh9/JpPeEQAKElgChSFG4AekWithARPHSFbiK\noALe61U/BC+igoggokgRpF0sV0AUu4YuEGqQJBBCS4f0Ni1zvj8OExKS6ZOQZPb7PPMwmTlz9j5D\nstfZa/3WWgpJkiQEAoFA4HS43O4JCAQCgeD2IAyAQCAQOCnCAAgEAoGTIgyAQCAQOCnCAAgEAoGT\nIgyAQCAQOCnCAAgEAoGTIgyAQCAQOCnCAAgEAoGT4nq7J2CKkJAQWrdufbunIRAIBPWKS5cucf36\ndbPH1WkD0Lp1a2JjY2/3NAQCgaBe0bNnT4uOEy4ggUAgcFKEARAIBAInRRgAgUAgcFKEARAIBAIn\nRRgAgUDQIElNTbX7HCqVipycHAfMpm4iDIBAIKjXaLVaZs6cyZQpU3jsscf47bffAPj3v/8NwIIF\nC8yeY/v27YwbN45JkyaxYcMGMjIyWLNmDV988QXHjx83+dktW7bwxhtv2H8ht4E6LQMVCAQCc6xb\nt46HHnqIoUOHotPpePTRRwkJCeHkyZOcOHGCxMRE9u3bx7p161CpVCxfvpzQ0NBK59i+fTvvv/8+\noaGhbNu2DZVKxeXLl8nMzMTDwwNfX182bdpESUkJzz//fLnM8tChQyQlJaHX62/HpduN2AEIBIJ6\nzalTp+jTpw8Arq6udOzYkUaNGtGtWzfuvPNOANLT03F3d2fSpEkEBgZWOcf777/PypUree655yot\n5n379mXEiBGsXLkSPz8/goOD2bdvX/n7d999NxMnTqzhK6w5xA5AIBDcFhQKhc2frdjKvFOnTvz5\n558MHjwYSZJITk4mLCys0vFdu3bljjvu4JNPPkGtVvPoo49Wen/dunX85z//AWD06NH06tWr0hy1\nWi2vvfYapaWl/PXXXzbPu64hDIBAILgtVFzE7eEf//gHs2bNYsuWLeUuGldXV3Q6HTExMYC8A1i9\nejVKpZIOHTqwfPlynnzyyXJD0aZNG0aOHImPjw+dOnVCqVQC0KpVK1atWsXMmTOZOnUqCoWCuXPn\nOmTedQGF5Kj/hRqgZ8+eohSEQCBwONu3b2fYsGF4eXnd7qnUCJaunWIHIBAInI6RI0fe7inUCUQQ\nWCAQCJwUYQAEAkGDRCSCmUcYAIFAUK+5nYlgW7Zs4fnnn+epp55yiMGpbUQMQCAQ1GtuZyKYn58f\nK1euZPv27Rw8eJDRo0ffjq/AZsQOQCAQ3BYUCtsfFbmdiWCPPPII165dY/v27Tz00EM18C3VLGIH\nIBAIbguOEqDfzkSwpKQkFi1axIcfflgvJaU1ugM4fPgwgwYNAuDkyZP079+fQYMGMXToUDIzM2ty\naIFA4CT84x//YMeOHUyaNIknn3zSaCLYggULyM3NLU8Ey8jIKD+HIRFs3LhxVRLB1q5dW54I9uKL\nL1ZyHz399NPk5uYyY8aM8rHqEzWWCLZ48WI2bdqEj48Pf/75JwMHDmT58uV0796dTz/9lMTERJYu\nXWryHCIRTCAQ1AQiEUymxnYAkZGRfPPNN+U/b9u2je7duwOg0+nw9PSsqaEFAoHAJCNHjmywi781\n1JgBGDVqFG5ubuU/N23aFICDBw/y0UcfMWvWrGo/t3r1anr27EnPnj25du1aTU1PIBA0cEQegHlq\nVQX0xRdf8Oyzz/L999/TuHHjao+ZNm0asbGxxMbGGj1GIBAIDNzOPIDdu3czdepUnn76aUpLSx1z\nQbVIramANm/ezKeffkpMTAzBwcG1NaxAIGjg3M48ABcXF1atWsX7779PYmJiuZu7vlArO4CysjJe\neOEFCgsLeeyxxxg0aBDz5s2rjaEFAkFdxUGJALczD2DYsGHs3r2bnTt30qpVq5r4lmqUGt0BtG7d\nmj///BOgQfvRBAKBDThIgHg78wAOHDjAww8/TLNmzVi/fj2zZ892yDXVFiIRTCAQ1GtuZ0OYS5cu\n8dlnn6HT6eploxjREEYgEDgdIg9ARuwABAKB0yEawsiIYnACgaBBIvIAzCMMgEAgqNfczjwAgPj4\neIYMGWL/hdwGhAtIIBDUa25nHoBarWbNmjX1NmlVGACBQHBbULypMH+QEaR5N7Urp06dYuzYsYD5\nPIAxY8YYzQNYsWIF165do3///uWv9+3bl/DwcFauXEnz5s3x9PRk37595QZg8eLFzJ49u3y3Ud8Q\nBkBQr9Dr9cTExKDT6Th9+jRdunTh1KlTKBQKunbtSlxcHF26dCn/t7r3DK916dKV335LpW3b9vz5\nZya9eoVy6FA6CoWC3r3DOHw4kz59mtGkSS+USheGDwcX4TR1GBUXcXu4XXkAJSUlHDp0iNTUVI4c\nOcLvv//Ovffe65Brqi2EARDUG7KyshgwYACJiYk2fDoKaAOsvfGvC9AOeLr8iA0bbh79+edVX2va\nFD75BB5+WBiCusTtygPw9vZm9+7dAIwbN67eLf4g8gAE9YCEhARmzZrFjz/+aMOnOwIfAvc7bD5t\n2kBCAri7O+yUglpG5AHIiB2AoE5z7tw5+vbtW35HZjkuwP+AEQ6f08WLkJQEnTo5/NSCWkLkAcgI\nAyCos2g0GkaPHk1oaCjx8fF4e3tTUlICwP3338+QIUPo2rVrJT//6dNx7N/fhd2726PTtbFwJD1v\nvvkh7dp14ezZOBQKBbGxJ9BoCtBqlRw8uI+yslQk6aYRevddWL8erLZLglojNTWV5s2b23UOlUpF\nSUlJw61gLNVhevTocbunILiNfP/991LHjh0lFxcXycXFRQKk9u3bSz///LNUVlZW7WfeeUeSAgIk\nSa40VvmRny9JmZmSdP26/MjJkSSVSpI0mqrnWb58uTRixAjJy8tLiouLkzw9/aTHHjtY6XxZWTX8\nBQgsQqPRSDNmzJCeeeYZaeTIkdKvv/4qSZIkjRs3TpIkSfrPf/5j9hzffPONNHbsWGnixInS+vXr\npfT0dGn16tXShg0bpF9++cXo53bs2CE98MAD0sSJE6XTp0875oIcgKVrp9gBCOosO3fupH///iQm\nJuLm5oZWq2XlypXcf3/1/vxLl+Q78/z8yq8rFFBUBN7e4O9v2dgtW7ZEp9MxZswYoqOjeeqpMYwc\n2Yvt2+XlX6mEkBD7rk/gGG5nHsDhw4cJCwvDy8uLyMjI23H5diEMgKBOotfr2blzJ3379gVkGd74\n8eNNKi3+9S949llYtOjma1evQmgoVOhOahEtW7bk6tWr7Ny5kwkTJjBgwACyslzw95dVQJs3w8SJ\nwg1kD4obCh1bkAYNKn9+O/MARo0aRZcuXYiJiWHDhg1Mnz7d5mu6HQgDIKiTxMbG4uHhwQ8//MCJ\nEyf429/+xoQJE3Axor/MyICffoJ+/W6+plRCeLht47dq1YrLly/j4uLCoBuLTVgYeHnJiz/Apk3w\n/vtQT5NAbzsVF3F7uJ39AGJjY+nevTtBQUHodDqHXE9tIgyAoE5y6NAhhg8fzqVLl+jduzdjx44t\nX4irY+NGCAqCG7JsAK5csX384OBgNBoNBQUF+FfwG/XqBd99B2Vl4OoqFv+6wO3sB9C4cePy15cs\nWVLr124vIg9AUCeZPXs2+/bt48yZMwwaNIhdu3bh6lr9/YokQYcOcO7czddcXUGrtW8OHTt25Ouv\nv6Zz587lr82bBzodLFggGwFrXUuCuoHIA5AR+YyCOsmxY8c4ffo0R48eZc+ePezfv9/osRcvyoHf\niq1iCwvtn0PLli25css2IiwMrl+XjY0ReySoB4wcObLBLv7WIAyAoE5SUlLC4MGDiY6OZsyYMQwY\nMMDosfv2ya6ZintZRxmAy5cvV3otOBh++AE6dpSDwGVlyFuCixfhmWfgqafgzjuhQiBRILCGtLS0\nWhtL3MMI6iRXrlxh+/btJCUlMWDAAKPBX4CYGDkIbMBRvvmgoCAKCgoqvebpKSuLQA4Cf7LvDnyu\nxMuDajSVT6BQQEAA5OXZPxknQq/Xs3fvXrP/7wbmz5/PhQsX2LRpEwDR0dEsXLiQzz//nOTkZLy9\nvWnXrh1TpkxhxYoVSJLE4sWLeeWVV9De8BMGBQWxYsWK8qBvTEwMSUlJTJkyBZBr/axdu5bnnnsO\nV1dXNBoNy5Ytw8XFhTlz5lBWVoZareaZZ57hvvvuY8aMGajVakpKSli4cCGtW7euNOdz584xb948\nfH19CQgIYMmSJSxYsIAJEyawZs0a3nrrLaPXGxcXx//93/+xc+dOW77eSggDIKhzlJaWkpeXR1hY\nGElJSWaPj4mRcwAM3HLTbjPu7u5oblnUOzTNR4kv/8cC7uEPfC7dUITcuvgbMPimBg6UJyowiV6v\nZ8qUKWzbto0xY8awdu1ai4xAfn4+KpWKq1ev4unpiUKhYPPmzWzYsIHw8HDuv/9+Dhw4wLJly/j5\n5585dOgQAJ9//jmurq589dVX7Nq1i0ceecToGBcvXiQwMJAPPviAU6dOkZGRwfr163n55Zfp0KED\nAL/++itFRUVcu3aNrVu3kpKSQlpaWhUD8Ntvv/Hkk08yYsQIvv32W0pLS0lMTGT//v3s2bOH7Ozs\n8mBzdHQ0kyZNAiAvL48vv/wST09PG77dqggXkKDOkZKSQvPmzZk2bRrDhw9nypQp6PX6ao9NT5dv\nsA3+eFdXaNbMMfNwd3dHrVYTExMjjx8XR9sn7yKFcObzJgPZW/VD8+bJi/3AgZVf37MH5s93zMQa\nMHv37mXbtm0cOXKEbdu2sXdvNd9xNTzwwAPs3r2b7du389hjj1V7TN++fYmPj2fJkiV079690ntd\nunSpVGVWoajcq0ChUNCxY0cGDRrECy+8wLp16wgJCeHSpUt06NCBhIQEJk2axOrVq/H19eX5559n\nzpw5vPvuu1WSzgCmTp3KlStXmDFjBidPnsT9RmXBvn37MnDgwPJkNH9//3JjBRAYGMiCBQtwc5D6\nQBgAQZ3j8uXL+Pv7W7QQHDwoZ/jqdDBunGN8/wZcXV3ZuXOnbIRGjUJ/3324eLgTRkblAwcOlBf+\nefPkRT4mRn7Mm1fZGr35pjACZhgwYEB59rW52E9F7r33XmJiYsjPz6820Qvg6NGj9OjRgx9++IGl\nS5dWei82NpZ27drx7LPPMnbsWAIDA8nMzATkmkKBgYHs37+f4OBgPvzww3JXTXh4OEeOHKFDhw5s\n2LABpVLJ+fPnuXbtGsuWLWPBggUsqpiZeIM1a9YwceJEPvroI7y9vTl8+DBw0/CUlZUxevRo3nnn\nHZPyZ3sRLiBBnePKlSt07dqVu+66y+xCcPw4GGJmmzfD0qWyn94RpKSkEB8fz4nffiN60CDe6NeP\n1vv23Txg3jz5X2OL+vz58PTTEBUFKlXlKLWgWlxcXFi7dm159rUl7h8ApVKJi4sLXbp0QaVSVXtM\nQUEBTz/9NO7u7kybNo1Vq1YxceJEJEnCz8+Pjz/+mBEj5Oqxer2eFStWMH78eHJzc1m0aBEtWrTg\nueee49NPP6W0tJS3336b0NBQXnzxRZYsWYJOp+Ohhx6iTZs2vPfee2zfvp2SkhJefvll9uzZQ2lp\nKcOGDQOgf//+TJgwAX9/f1xcXMoziIOCgti3bx8bN25kzpw5fP3110ZLnzgCo3kAt1rIW5k9e3aN\nTKgiIg/AOZk/fz56vZ758+ebDQZOnAhbttxMzLJX+1+R5cuX8+mqVVw6f54xHTuyLjgYxY2dSFGP\ngaRsjiEqqrL8tFratoULF+Tnhl2CwKm4ePEi6enp9OnTp1bGs7sfQFxcHF999RVPPPEEt9qIW/1j\nAoEjuXLlCnfffXelMgzGSEiAX3+FJk3kddaReHp68lyLFnQpKWFgYGD54r+h1Ty2BM3n144wfrwF\n9YBmzIBZs8Ti78S0adOGNm0sLU9ei5gqFXrvvfdKmzdvtqssqT2IctDOyYgRI6T//e9/Zo/T6yXJ\nz0+SnnhCLs88frwk6XSOm8dna9dKV4KCJGnMmJs1oAcOlB55RLKuLPTZs5Lk7n7zA/PmOW6SggZH\namqq3eewdO006WBbuXIlf/75Z23ZIoEAgKKiIry9vW+qb4yQkSG3ZfzyS/nnTZsgJ8dx84g8fRqt\nQnEzyAAwaBAtWtzsCWxRzkGHDvIWBarsAiRJ4sekHxm+ZTh+7/jh+ZYnvdb0IqMoo/pzCaowf/58\nxo8fX/5zdHQ0v/zyC+PGjaNPnz7cf//9PPfccxw7doxJkyYxceJEMjMzmTRpEmPHjmXs2LHMmDGj\nkqcjJiaGtWvXlv88btw4VCoVkydPZurUqUycOJHc3Fzy8/OZMmUKkydP5u9//zu//fYbADNmzGDq\n1KmMHTuWSxU1yjc4d+4cTz31FFOnTuXll18GYMGCBVy+fJmPP/7Y6LXGxcUxadIkJk2axM8//2zv\nV2daBdShQwdWrFhh9yACgTUUFRWxYsUKsxLQ+Hi5LWNFCagji7N1+OUXMry9waBAurF4N24Mr7wi\nxx1uNCgzjUIBgwfLPqo335R/nj8flU7F418+zshtI/kl+ReKNcWoy9QcSTtCm2VtOJ152nEXU4+Q\nJNm1Z03M3JAHcP78+Up5ANOmTePf//43n3zyCSqVimXLlvHggw9WygPYsmULAwcOZNeuXSbHMOQB\nrFmzhtmzZ5ORkcHChQt5+eWXWb9+PVu3bkWSpPI8gFWrVvH222+Tnp5e5VyGPIA1a9YwYMCAavMA\npk+fzvTp09mwYUP553Jzc1m2bBkffPCBQxLBrJKBzhf+S0EtkJmZye+//25WApqQIAtsjh+XZaAW\nLcaWcv48/mlppHh7V3mrSRN5p+HiYkUxuLvuAkMvA0lCO/d1Ht32KCcyTqAqU6HT65C4ueKpylT0\n+6yf0+0EysrkwH6lUhsW4Cx5AAMGDMDX15e5c+cyY8YMy74cE1hlAL799lu7BxQIzCFJEg8++KBZ\nCWhiorz4d+0KkyffdMs4hPXrSbv3XrSG3YcklbtuGjeGrCwrz9etG5w6VZ4gNv376eglPWmFN91L\nfu5+SPMk/Nz9ACjUFNJ7bW+nMgLZ2bIrD6xz6TlLHkBpaSnPP/88//jHP+jYsaNlX44JrMoDkISO\nWVALlJSU8OGHHzJjxgyTEtDkZNkAgIObs0gSbNmCW9eujDWUopg/v5IBuHbNynN27SpPVqsFhYLg\nx/zZfKcGnV5uItK3RV/2Py1XPC14tYCOKzuScD2BK/lXmPfHPD59+FMHXFjdp0kT2ZWn01nn0nOW\nPIB33nmH+Ph4Fi1aRN++fXnuuefs+Lat7Acwb9483nzzTbsGtAaRB+Cc+Pj4kJmZia+vr8nj+vSB\nI0dqIAcgLg4efZTLAwbQ6vPPqzijT5yASZPkG3qraNMGLl0iZFEjSnWllGhln1XFxd+AJEm4L3BH\nJ+lo5tuM1DmpdlxQ/UKrlWW1DanfQl3NA7Bq02xY/F999VXbZiUQmEGv11NaWop3Nb73W0lNhbNn\nrQjGWsp330HjxvLiD1W0+/7+cskJq4OV3brx/bAIWvi3KF/8/dz9qiz+ILsCwv3lfpYluhKncgO5\nuVkZX6kHtGnTptYWf2uwyWu6cuVKR89DIABk94+Xl5fZEgB6vVwIrlWrGlgsvvsObvSH7XHXXVUM\ngJ+fXOTT2mDl2W7NmNArnYsZ8YB851/waoHR4/c/vR8PpQcKFCw+sNjmyxEIjGGTARCxAEFNUVRU\nhI+Pj9njsrLkUvvu7tZLBk2SkyO7gG7U/Dl2/Hi1BqCw0Ppg5cKgOPpe1lOsVzMtFvZfMV3jpbl/\nc9yV7uSqctl/xXhHNEHDojYbwthkAEQpCEFNUVxcbNb3D7L7p1kz2ySDJjl4EKKjy5MLOkRFVTEA\nnp6ywbEm/yA5N5kfSk7zUws1OiV8d08zi8pC3Bl2Jx5KD/q17GfDxTgHzpQIduXKFcaNG8ekSZOI\nj4+396szrgL65ptvqn1dkiTKLPxLO3z4MK+88kp5d51JkyahUCi44447WLlypcWV/gTOg2EHEBMT\nY1IBlJoqL7oV78IdogI6cECWoNyIdz13QwpYEYVC3gX89ReEhlpmeJYeWkrb4LYcT5NlS20CLasL\nU6gpRF2mZv3J9fyr778I8w2z/FrqMxVUV5bgLA1hCgoKWLhwIZmZmfz00092S0GNrsArVqyo9vHR\nRx8RHR1t9sSLFy9mypQp5ZKs2bNn89Zbb7Fv3z4kSXJIFpug4VFYWEhmZqbZLODUVIiMrIEs4P37\n4fXXAUi5epX3jLij/P3lCs+WxB+KNcVsjduKhjLKlPC3srBqA7/VcWzaMXzcfOQ2hs4SB5g/v1LG\ntCU4SyLYHXfcgUql4oUXXqBfP/t3hUYNwB9//GHyYY7IyMhKu4hjx44x8EYSzAMPPMCvv/5q9+QF\nDY9Dhw6RnZ1tNgs4NRWaN5fVPw5TAanVssbz999h3rxqW0IaMMQBLOHLv76kR7Me5KnyCFQr6Hck\nw+LFTaFQ4OPmQ746n/Un1zuHGsjwvVRIvjOHsySCxcfH07hxYw4cOMCHH35o0XdjihprCDNq1KhK\nvi9Jksovzs/Pj/z8/Go/t3r1alavXg3ANauzbQT1nYiICMLDw81mAaemQr9+N+++HeJNPHZMXtnf\nfRcAf5XKpAEoMC7gqcSa42to7N2YlIIUyjwkjjTHqqj14sGLeXvf21wrucbiA4tZOtR0r44GgaHZ\njoU4SyKYTqdj+vTpBAQEMGrUKFu+2crYXXfUBBcvXpR69eolSZIkNW/evPz1HTt2SM8//7zZz4ty\n0M7Hxo0bpb///e/SH3/8IZWVlRk9bvBgSdq9W5Li4+Wy0A5h2TJJevZZuWSzJEnFxcWSt7d3tYcO\nHSqPb47knGQpeFGwFPhOoMR8pL5vhEvSwIFWTSvhWoLkscBD/vy6vlZ9VlA3SE5Olg4cOFBr4zmk\nHLQjufPOO4mJiQHghx9+oH///rU1tKAeUVRUhJ+fH4MGDTIpEkhLg48/drAC6NQp6N69vF6Pm5ub\n3S6gr85+RTO/ZigUCmb1nsX+1v+RkxesoF2jdni4ejCo9SC+fuJrqz4rqBvU20SwvXv3Vnrs27eP\nY8eOUWDp/vcG77//PvPmzePuu+9Go9Hw+OOP2zxpQcOluLjY4jyA776TnzusD8CpU3D0KOzZAwoF\nrm+9hU6nqzYQbakL6Mu/vkRbpr2p5Y+IkNtDWqFwcVG4EOQZRGxarPMEgp2Y2swDMBsDmDVrFqdO\nnaJz584olUri4uJo2rQpJSUlrFu3jkcffdToZ1u3bl3eUKZ9+/bs2bPHcTMXNEiKiorM5gFIEuTl\n2VY0zChardxgYO9eWLcOJAkF4P7uu2i1Wjw8PCodbigHYYqLuRe5kn+Fp7o8RVZxlqzlj4yEkydl\nuembb1rcJlIv6SnSFHEk9Yjt19hAmT9/PhcuXGDTDU1wdHQ0Cxcu5PPPPyc5ORlvb2/atWvHlClT\nWLFihayoWryYV155Be2NAlJBQUGsWLGiPE5pkK5PmTIFkPMA1q5dy3PPPYerqysajYZly5bh4uLC\nnDlzKCsrQ61W88wzz3DfffcxY8YM1Go1JSUlLFy4sIoM9Ny5c8ybNw9fX18CAgJYsmQJCxYsKA8u\nv/XWW0avt7S0lB49enD69GlcXe0L45r9dKtWrXj//ffLI9FHjhzhgw8+4L333uORRx4xaQAEAmsp\nLi6msZnVvKhIDv7m5d0sGmY3iYnQsiX4+FQKQBqUQLcaAEM5CENPgupyI3ef382g1oPYeGojeao8\nefG+vykY3EpWBII/fehTFuxdgK+7LxlFGQ06HyC9MJ3JOyezYcQGi6/TWfIAAN59990q57MVsy6g\n5OTkSjKk6Ohozp07R3h4uEMmIBBUxJIdQHY2hIQ4uGjYyZNyzf5bEpDc3NxQq9VV2lP6+sJXX5mO\nQexO2k2xphgFN/z/T++XLVZoKLz4olXT69GsB8fSj3Ek9UiDdwMtPrDY6ut0ljyAnTt3Eh0dTRND\ni1E7MWsA3NzcKvWe/Pnnn3F3d+fatWvl2yeBwFFYYgCuX4dGjRw88KlT8olvSUByc3Nj5syZVRLT\nFAq5EilUH4Mo1ZYScymGA1cPkKvKrey6CQ+H0aOtml4TnyZIklT1XCaQJImuR4+iiIlBERNDz3pS\nWv1o2lGrrhOcJw/gu+++Y9euXRw6dKhcLm8PZl1AK1eu5PHHH0ehUKDX6/H09OTrr7/mvffe49ln\nn7V7AgJBRSwJAmdnQ3CwafeL1Zw9Cy+8ICeBVXDNSJLEjh07OHr0KNHR0UyYMIFBgwbRvLk8rqEm\n0K1eq5hLMQR6BlKqLWVW71mVtfvh4XIig5XlDu5pcw9/v+Pv/PfMf826gZ5JSOCzjMpJY8eKinCN\nieGl8HCWtG1r8bi1jaVZ0hVxljyANWvWADBp0iSmTZtm9fdUBUu0ojqdTjpx4oQUFxcn6XQ6W6Wp\nViPyAJyPwYMHSz/++KPJYzZtkqRWrWS5/vjxkuSQX8nISDmp4BaNfkREhDRq1CjJy8tLmjx5cnlu\nwrffStIDD0hSWZkkaTRVTzfnpzlS8/ebS8xH+tvqv1V+86WXJGnIEPkCQJLmzbNoiv/8+Z9Sn7V9\npKB3g6RZP84yety54mKJP/4w+nD54w+LxhM4jrqaB2B2B5CRkcGGDRvIuWWPu3hxw/ZDCm4PluwA\nrlyBy5fl5w4pAqfRQEoKbNlSLgE1qHM8PDx48803q7Sn9POD4mI5BlFdukLMpRj6tOjDr8m/Vq3k\nGR4uj/Hzz1YFgruGduWzE5+ZdI8U6nREG/pk3mDOjXjd+ykpAOiBl5OS6vQuoKHRpk0b2rSxrABg\nbWLWADzyyCOEh4cTGRlZG/MRODmWxAA0GnnR1esdJAG9cEFWAC1YAG+9VWlRdnNzQ6fTVfLDguk8\ngDxVHonZieW1/Kss1uHhcPiw1eUOuoV2I8Q7hOv/um70mEEnT5Kn09Hey4tzpaXMucXd835KCkFK\nZbkxEEbERyJvAAAgAElEQVSg7pGWlkazZs1qZSyzBkCj0RgtDS0QOBpL+gFcvw5LlshCGodJQNu3\nl5/fsigbKwhnKg9g/5X9dA/tjo+7D+lz0qv66sPD5R3Hl19aNc2okCgu5V1i8KbBbBq5qcp5TxUV\ncbyoCAVUu/gvaduW4rIyVqWn4wJ8kJLSIAyAM+UB5OXlcd9999GlSxcefPBBRlspJrgVsyqgHj16\ncObMGbsGEQgspaioCE9Pzyqyy4pkZ0OTJg6UgCYmytHkaoKyxgyAqR1AzKUYyqQy41JGgwGwEnel\nOz7uPkbPO+rG36kC+Q+7usX9k6goFMhuIIMrqK7xclISHnv2WDU3Qx7A+fPnK+UBTJs2jX//+998\n8sknqFQqli1bxoMPPlgpD2DLli0MHDiQXbt2mRzDkAewZs0aZs+eTUZGBgsXLuTll19m/fr1bN26\nFUmSyvMAVq1axdtvv016enqVcxnyANasWcOAAQOqzQOYPn0606dPZ8OGDeWfO3LkSHmezK1SVlsw\nawD69u1L9+7dadGiBREREeUPgaAmUKvVvPLKKyb7AWRnO1gGmpgI585VW4Pelh3AH5f+IC4rzriv\nvmlTyMyEN96weqoKFBSoC6qc90RhIRdUqvLFfZaJPJ1nb7gXRjRqxIrUVKvnUNOsSE3lp65drZqb\ns+QBREVF8dlnn7Fs2TIWLlxo8fdjDLMGYP78+WzdupUdO3bw1VdflT8EgppApVKxfft2k/0AHG4A\nzp2DOXPk57fUoHdzc6s238XLS64ecetbap2a05mncXNxu5n8dSvu7vLWZcECq5qeAEz/23TmDphb\n5byLrlxBAbgqFEbv/g183L49UV5e7MzOpkyS6twuYGbz5gw9fZqZzZtb/BlnyQNITEykpKQEHx8f\nh7TmNRsDCA4O5oknnrB7IIHAEvR6PU8++aTJfgAONwDnz0O7dtUGZY3tABQKuWpEcTFUXG9OZpw0\nHvytSFQUHD9udTf7jiEd2RK3hWGbh5WXSkhXq9lx/bp89y9JJu/+DayOimLgyZO4KBR1LhawpG1b\nq+fjLHkAbdq04f/+7//w8fFhxowZtny9lVBIkunfwNdffx21Ws2oUaMq1UO566677B7cHD179iS2\nnmQvChyDi4sLarWaAwcOGO0J7OcHV69CRoYDEsFKSuSsspKSavWcjz32GOPGjavWrdC8uSzmqbje\nfnj4Q46mHeVa8TXTtWwefBACAmDrVqumeyL9BIM3DUYv6ZnUfRJLhy6l7/HjHCwo4ImQEHZkZ6O+\nUc7aFJIk4bJnD8obP+tuUTkJHMvFixdJT0+vtZLQlq6dZl1AW7du5X//+x9///vfGTVqFKNGjRKl\nnAU1QtkNSY+bm5vRfgBqNZSWykm7DukFcOmSXJ//P/+p9m1TbSF9feXCdAkJN2/k/0z5k6yiLPO1\nbEJD4b77rJ5uVEgUOaU55TuMMkniYEEBQ4KC+PL6dYvdJgqFgt7+/hgiLHXNDdTQqLf9AC5evFjl\nkZycXBtzEzgZWq0WNzOyntxc2eVyQ/Fnfy+AixflmtJGmpCbMgDe3nLooKIh+jPlT66XXDfvAgoN\nlQPBVuLt5k2rwFacn3me/U/vZ09eHl4uLvycm0sPX1+rXCdr2rdHAn7u1q1OBoMFNU+tdQQTCMxh\niQEoKICgIDkBDByQCHbpEgweLD+vpgm5sSAwgIcH7N4tP9+0CS6k5pJVnEWITwjpc9JN17QJC5N9\nWFYEgA20DmzNE189QUZRBi+cP0+pXs8TISHEFRdbdZ47fH3xVCi479QpuljQhEfQ8BAGQFBnsMQA\n5OfLrvOSEvmOu6TEzkEvXoTWrY1m5ZraAQQHy9WdQTZEV3XHCfAM4GjqUfOljEND4bffjO48TJFX\nmkfC9QTe3v8uf5WU8HBwsFXun4poJYmnmjSx2ngIGgbCAAjqDJbuAPz9HdgL4NIlaNPG6AJsygD4\n+MDGjTcN0YmME5aXbA4NlZsaQLU7D1MUqAso1ZWypSQYgN05OVWyfi3liSZN+G9WVp2UgwpqHmEA\nBHUGS3cA/v4OHNSwA7DBAPj6gkp10xAdTz/OK31fYWjkUPPN28PC5BiAlfWAAD4a/hFDIoeQFzyA\n55s1Q6lQ2Czj3NSxIwCbO3YUcQAnxCYD8OCDDzJy5Mjy5AWBwBFYugMICHDgoJcuwddfG3XFuLm5\nmVUBGTiRcYLYtFjLulkZgsA2xADaBrcl/noi+pKrrExLs8t/r1Qo6OTtzfiEBJtcSIL6jU0G4Nln\nn+Xrr782W7ZXILCGWt8BFBTIt/CG0ubVuGLc3d2NBoENiWAARZoiruRfITk32TIXUFCQ7DcykrRk\nilaBrUhRlYBnGI9rDtjtv3+9VStaeXiwIjVVuIGcDJsMwMMPP4xSqeSOO+5w9HwETkyt7wCuXJHL\nQBvq/1eDJXkAAHGZcbQNaoufh595BRDIYzZpAllZVk/bXekOrSaCwoWv3e62+879nsBALqhU7Lay\n/o6g/mPWAJw9e5bJkyfz6KOP8sgjj5Q/BAJHY6kB8POrnHxlMykp0KKF/NzGGIDBAJzJOoNO0lnX\nzDw01GYpqBQ2nOcDy3B3UdpdxqGphweeCgX3Czmo02HWAEyYMAE/Pz9GjhxZngk8atSo2pibwMnQ\narXlVRGNkZsL27c7KAs4JUWu42BiATZnAAzel7isOEq0JdY1Mw8Lgw8+sFoKmlBcjCs6VuW70qzg\nTzKKMsx/yAwaSeLZZs2EHNTJMFsMTq1W8+GHH9bGXAROjiU7gOvX5fo74IB2kCkpcjew9evlRfhG\nG8iKmAoC+/hU3gF8+tCnDIkcYvn4oaHQuzds22bVdmbauXOUKVwJyYkh//wKFrtnVm46bwPDg4NZ\nlZaGEtEu0pkwuwNo165dtQ0NBAJHY4kBUKsrJ1/ZlQWckgJjx8rPjWjxTQWBK7qATmWe4p1971h3\nN26jFHR/fj5jfVVcC+xn3Y7DBJ916ABgdR1+Qf3G7A5Ar9fTuXNnevTogZeXV/nr3377bY1OTOB8\naDQai2IAv/4KAwY4oB1kSgo89pjJBdiSGMD+Y9cp1hRzKvMUiw8stvxuPDRULkW9YoXFUy6+cdGb\ni7zwUF1CNc/eQIhMY3d3Grm68kBcHC8IOajTYNYAjBw5kpEjR9bGXAROjqUy0EaN5OSraoqFWoch\nBjB8uNFDTBkAT084cwb69wzB7cVIcoPOWnc3HhoK+/ZZNeUD+fkA/HhHR4ad1lCmL0PpojTzKcv4\ne2goJ4uKyncAwg3U8DH7JzRx4sTyjjRarZa+ffsyceLEmp6XwAmxphSEQzAYABOYigHodHJzGgDt\n8r/Imi6Zl39WxOACsoL/u3gRgAfPxOOT9TP3brzXIUFgkOWg+/PzhRvIiTBrAH766Sd69uzJjh07\n+Pbbb/nb3/7Gzp07a2NuAiejVvMACgtBo4Hly00eZmoH0KrVzeeKgBTG/zTMusXYBhlobGEhH0RG\nolQoCEzbxon0E5bLTs0wMDAQV4XC6naMgvqLWQMwd+5c9uzZw/bt2/n222/Zv38/823QLQsE5jBn\nAPR6ed3283PAYKmpcmPf//zHpATTVDloHx+5N8E96+/nkU9esC4HAOQdwOXLFstAc7VaXBUK/pWc\nzMzmzdHpdRRqCh0SBAYIdnMjQKnEMVEFQX3AbAxAo9HQqVOn8p87d+5c3rlJIHAk5gxAcbHchEXp\nCJd3Sgp06wZ//GFSgmnKABjyAM7lJNDUt6n1ipyAgJv9LC2Qgcbk5dHc3Z20GzuSl3q/RHZJNu8N\nec/yMc2Qo9MxvVkzVqSmihiAE2B2B+Dl5VWpt2RsbCze3t41OimBc6LVanF1dSUmJga9Xl/lfYfW\nAUpJkZv6mpFgmjIA7u6gdy0iV5XL4amHkeZZGQMwlIN48UWLDv8tN5eranW5j75tcFvOXDvDsM1W\nup5MMCQ4mI/S0kR5aCfBrAFYvHgxDz30EP3796dfv34MHz6cRYsW1cbcBE6GWq3m4MGDDB8+nClT\nplQxAoYAsEPKQKSlQbNmZt0upgyAQgFe4ecI92vF8C3DbVuEw8LgqacsOvS/WVlIwP2nTjGzeXMi\ngyI5mnrUeteTCZbfuOv/VbSJdArMGoD+/ftz9uxZXnvtNV599VXi4+MZOHBgbcxN4GQkJCSQnJzM\nkSNH2LZtG3v37q30fk6OnAnskDIQ6enQtKnZw0wZAAC3polodVbWAKqIIRBshjS1mhydjl+6dSuv\n/x8RFEGuKtdhyWAA7by88BSBYKfBqAHYvHkzAEuXLmXDhg3Ex8eTmJjI559/ztKl9qWdCwTV0bJl\nS6KiooiOjmbMmDEMGDCg0vupqbIBAAc0gzcYADM7AHd3d9RqtVG3lEvjREq1pbYvwhZKQffl5xPp\n6ckDFRbmAM8AgjyDyJiTYZ3ryQQKhYJmHh5UvVJBQ8SoATh//jwAcXFxVR5nzpyptQkKnIeysjKG\nDx/O7t27Wbt2LS63ZHq5uNyMmdpdBiI9HX780awCR6lUkpKSYtQtpQ9K5Ln271jv/zdgoRR00ZUr\nXFKpmNm8eaXgbGRwJMm5ydaPa4LLKhUjGzUSLiAnwKgBePPNNwEYMWIE69evr/S45557am2CAufB\nUA100KBBVRZ/kMsujB/voGbw6enw+uvycxM9eU+cOEFhYaFRt5Q68AxfJK+0PQgbGgq7dpk1RCeL\niljWtm2VRTkiKIJj6cccGgh+rHFjvraxybygfmFUBrpr1y60Wi3//Oc/0ev1SDeiblqtltdee43x\n48dbPZhWq2XixIlcunQJpVLJmjVr6HCjCJVAYE4GWlQkB4HtLgMhSTddQGZUQH379sXLy6tat5Qk\nSZR6J3C5xNe6GkAVCQuTm9LExhqNbBfpdCiBORcuVFmUIwIj2Bq3lYTrCbbP4RY+i4pi+7VroiSE\nE2DUAJw8eZLff/+drKysSuWgXV1d+ec//2nTYLt370an03Hw4EF++eUXXn/9df73v//ZdC5Bw0Or\n1eLr62v0/eJiWXtvN/n5chd3Hx+zrhcPDw98fX354osvGDBgQKWdSVphGigUFJfl2R6EDQ2Vu4KZ\nMESHCwtp4u7O9WqC0RFBEaw9vtahgWBfV1fKgGWRkcy+cEEYgAaMUQMwd+5c5s6dy8cff8z06dMd\nMlj79u3R6XTo9XoKCgrMpv0LnAtz1UCLiuQ1224sVACBrALS6XTl9bAqkpidSBP13bwVuYd+feQb\neEOMwmKaNDHbHP5Afj6ZGg2/duvG0NOnKy3IEUERRIVEkTnZuppC5ujq48OLSUm8ZKZWkqB+YzYT\neMqUKWzfvp2ioiIkSaKsrIykpCQWLlxo9WC+vr5cunSJDh06cP36db777rsqx6xevZrVq1cDcO3a\nNavHENRfLHEBhYQ4YKCKCiA78gASryfir2/JnBPDKJi1gfEjw1i/3spM5dBQsyqgz9LT0SPr/29d\nkCOCIhweBAZ4uUULFly+LNxADRyzntQnn3yS+fPnM2fOHLZs2cKsWbP466+/bBrsgw8+YOjQoZw7\nd45Tp04xceJEVCpVpWOmTZtGbGwssbGxNLZL5iGob1hiABziAkpPl/WkFtTgMWkAshMpVqZS4HcE\n+i62TZoaGAgqlfyohjJJ4rJazTedO5fr/ysS7h9OVnEWQzYNcVgQGCDa359zpaWiMmgDx6wBOHny\nJMeOHePRRx9l2bJlHDx4kLy8PJsGCwoKIuBGKcfg4GC0Wq2oKyQox5JaQA5zAQ250brRhAIIzBsA\nrUcGeOVC8yO2SVMN5SBefbXat88UFxPk6sqTZ89Wq8pRuijxdvN2aDYwyAlhHiIhrMFj1gA0a9YM\nV1dX2rdvz5kzZ+jUqRMlNmrwZs2axfHjx+nfvz/33nsvb7/9Nj4O+YsWNARqdQdggQIIZNFDWVlZ\nuQquIonXE5nqt4PpWRJla/bbLk3V6WDZsmp3I3OSksjX6aro/yuiUCjIV+c7LAhsOGdTd3eRENbA\nMRsD8PHxYevWrXTr1q1ctplt6IJhJb6+vnz55Zc2fVbQ8LHEAPj4yLWAoqJsCLgaSE+H7t3lpAIz\nKBQKXF1dy3MUDGjKNKQUpLDDbSZd1J/j4hJmuzT1zjvhhx+qlYH+kZfH7BYt+DAlxagBeOqOp+jU\nuBMzomfYOIHquapWM6ZJE1EZtAFj9ld25cqVnDx5ksGDB+Pi4sLAgQNtloEKBKa4dZG9laIiWLLE\nAbWAMjNl/b2FVOcGupx3GS83Ly5rj3LC207XS2goPPJItW/5u7qyPCXFpBumpgLBj4SEsDUrS7iA\nGjBmDcBnn33G4sWLUSgUfPHFF+Tk5PDss8/WxtwEToa5HUB+vnyjDHbWAsrKgm++sfjw6gzAxbyL\nKFBQrM/lmrudrpcmTaB37yov52q1FOl0mNvoRARFcPbaWYdmAwN82r49bgoFK1JTRWnoBopZA1Cd\nVFMgqAnMGQCVSq4BBHbWAkpOhlWrLOrCBdX3BU7OTeaJzk/wx0CJLkftLMRmRAp6vKiIMjCrxIkI\niiA2LdbhgeDG7u5oJIk17dsLJVADxWwMICIigiFDhtCvX79KWZqzZ8+u0YkJnA9LVEBXr8o3zDa7\nf3Q6uRcwWNxUwN3dvcoOIDk3mYigCHw85HnZRWioXAriFmILC7nT19esEqdNYBtyVbnoJb1DA8EA\nUV5ePHPuHC8KN1CDxKwBCA4OBuDixYvlrylsjr4JBMaxJAgcFGRnLaDsbPkkVmS3V+cCSs5NJjIo\nktnHhpGn2wBYHlOoQsVs4Ao7ks/S00lWqXjRhAIIbpaFPvv8WZr4NLF9HtXwj2bN+DwjQySENVDM\n/hmJaqCC2sKUAdBq5abwJmLElpGZKS+4Frh+DBgzAHsu7yEu5wgZbR0QBI6Lq5KYdq60lHUWul8i\ngiI4nHLY4XGAXv7+nC4uFglhDZRarQYqEJjClAEwJIHZvfnMzJQXXCu41QBIksSF3Au4K93J1+Si\nDLbT7RIaetMddePfbK0WN4WCqefOWaTCiQiK4OOjH5fHARxRFRTgTl9flCASwhootVoNVCAwhali\ncA5LAsvKstsA5KpyUaDg0DOH0GoV+PiAfhWcO2djfkJICOTlwdy55S8dKywk1M2NLBPtKCsSERTB\nltNbHFoVFMBLqSTYzY1cnc5h5xTUHYy6gObOncsff/zBokWL+OOPP8ofv/zyCy+88EJtzlHgJJja\nATjMADjABZScm0wL/xY8sOUBcjSyu2X8eDvyE5TKKnGJ2MJC0jUai10vEUER3Btxr+2dyUxwXatl\nWtOmwgXUADEbBJ40aRKff/45OTk5ldLhhQpI4GjMGQCHVA3JzJT97b//Lvvc582zuiJocm4yqjJV\nubvFx2cpW7fK723aBO+/b4NE1SAFvZGgFltYyLDgYItdL5FBkWw6vcnKQS1jcFAQn6SlMUuUhm5w\nmDUAY8aMITU1lS5dupSrf4QKSFATmIsBOMwFNHasbAAslIFWZwAqNoL385Pnp9PZkZ9gaAxzg19z\nc1Hp9bxgRgFkICIogvPZ5xm2eRgbRmwgzNcOVdItfNC2Lb2PHxdKoAaIWQMQHx9PfHw8rq5mDxUI\n7MJUKQiHu4AsKARn4FYDcDH3Iq/3f53n/vYcAB3fk2X8XbrYkZ9QQQqa+dprFJaV8Xu3bgy7pQGM\nMQxloR0dBAaI8vamoKyM7Z078+TZs8IANCDMykBbtGhRG/MQCGrPBRQaal8MIE9OAjPg6wtqtZyb\nYHOTu9BQ2LgR3nyTUA8PhiUkMMwK5Y3SRYmb0s3hQWAAF4WClh4ejDZSklpQfzF7W9+lSxfuuece\nhg0bhpeXV/nrIgYgcDSWyEDtrgTqABXQuexzLNy3kG5h3QjzDcPHx0HZwO7u8Msv9D12jMMFBbxk\nofvHwIBWA3ip10s80O4BOydTlXGhofyelyfcQA0MszuAgoIC2rZtS1JSEnFxccTFxXHmzJnamJvA\nyTBlAAoK4NAhOyuBSpJsAJpYly1b0QDo9DpSClI4k3WmvO6Or6+8Q7ELgwto3jz+LChgbqtWVqtu\nIgJrpiooyAlhhwsKREJYA8PsDmD9+vUA5OXlERgYWOMTEjgvpgxAVhYkJsrPbVba5OWBpye8+67N\nLqCr+Vdxc6nsanGIATCogDZswGvvXt6+csVqd0tEUASnM0/XSCA42s8Pd9EhrMFhdgdw7tw5Onfu\nTOfOnUlLS6Njx44kJCTUxtwEToYpA6BQ3Kz/Y7PSJitL/rAFvYArUrEaaHJuMne3uLuS3t5hLqDM\nTNLUasosVCfdSkRQBL9d/M3hVUEBwjw8cHdxwbaZCeoqZg3AjBkzWLZsGU2aNKFZs2bMnDmTadOm\n1cbcBE6EJEnodDpcXFyIiYlBr6/cjLCoCN57T3b92Nx6MTNT9iHJA1psACpWA03OTSYiMKLS+w7b\nAWRlcaywELUk2eRqiQiKIL0ovUYCwQDFZWX8MzxcuIAaEGYNQHZ2NoMHDy7/efr06RQUFNTopATO\nh06nQ6lUMnXqVIYPH86UKVMqGYGiIvD3t1NpY1AAWSEBhcouoOTcZEK8QyoVXXOIAWjcGLKy8F2w\ngGg/P5tcLRFBEShQoH9D7/BsYIB+AQG8e/WqcAE1IMwaAIVCgUqlKk/+ysjIoMxmsbNAUD1arRal\nUsm2bds4cuQI27ZtY+/eveXvOyQRzKAAssL/D7cYgLxkTmedruRm8fGRDUBCgsW5ZVXx9AQXF77X\naDiRnc3MK1esVtoEeAbgrnTnvo33ObQiqIG32rShiZub6BDWgDBrAKZPn87QoUPJysri1VdfpXfv\n3ky3opa6QGAJWq0WDw8PxowZQ3R0NGPGjGHAgAHl7zskEcyQBGYlFQ3AxdyLpBWmVXKz+PjA9u32\n9yqWWrVixciRbOrWjRUtW9p0DnelO7FpsQ6PAQDc5edHhlbLri5dhBuogWBWBfT000/Ttm1bvv/+\ne7RaLWvWrKnkEhIIHIFGo8Hd3Z21a9cyYcIEBgwYgEuFri8OSQTLzITu3a3+mJubG2q1mpiYGC5k\nXyB+ZnylxitKpVwJFOxQKAEp7drR+9w5xnt42OxmKZPKKNQU1kgMwEeppLGrKw/FxfGCcAM1CMzu\nAAoLCzlw4ACLFi1i5syZfP/99xTbLXkQCCpjKAPh4uLCoEGDKi3+4KAdgI0uIFdXV7766iuGDx9O\n/rf5jP9mfCUXS3j4zcQ0e3oVx3bujG/btqgHDrQ50WrqXVNZcM+CGokBADT18EBvs59LUNcwawAm\nTZpEdnY2AIGBgSgUCqZOnVrjExM4F5b0A3aIC+i776yWgaalpXH69Gk27t5IWWAZh64cquRiCQiA\nfv3sVCgBi3r14icPD7v86xFBNZcMBnC2pIR7AgOFC6iBYNYAnD9/niVLlgAQEBDABx98wF9//VXj\nExM4F5b0A3aIC+jVV+XnVshAIyMj6dSpE2MfGItHlAeFusouFn9/eX52KZSAI0FB/OfcOVZcuWLz\nOSKCIoi/Fu/w1pAGxjZpwm95eUIJ1EAwawC0Wm0l2WdRUVGlvgACgSOwxAA4zAVkpQzUw8ODIUOG\n8PR7T/Ps4GerNF3x94f8fPumJkkSHpLEm5GRzPzqK6t2KBWJCIrgzLUzNZIMBrCuQwdcFQqhBGog\nmA0CT5gwgV69ejF69GgUCgXffPMNkydPro25CZwIS1xAdhWDKymRO8v7+dkkA9XpdCjaKIhsFFnl\n/YAAuVaRPVxWqVAqFJRnPth4k9XCvwXFmmIkpBoJBCsVCnSSxLtt2jD30iVRFK6eY3YH8Oqrr7Jo\n0SLy8/MpKipi8eLFohKowOGY6gWg0YBeD1On2iG1NCSB2VBG1CADjb8Wz4aTG6q4Vgw7AHs2xrGF\nhZQqFPz00UesGD3a5vMoXZR0atyJE/84UWOB4J5+fvz74kXKJEnsAuo5Zg0AwLBhw5g/fz6vv/46\nvXr1Iicnp6bnJXAyDDLQ6jD4/zfd6Hi4aRNY/StoMAA2YDAAJzNOkpSbVMW14ukp2xW12qbTA7IB\n6O3lxdDnn2emjTkABqJCovjz6p81Fgf4d8uWKIBfu3UTweB6jlkDsHz5cgICAggJCaFx48bl/woE\njsSUASgulu+yDU3pbJJa2tAIxoCbmxtqrZp8dT4F6oJqXSv+/va5gf6blcURlYqZO3eyxM7eu1GN\novj89Oc1Fgfo5e8vKoM2EMwagA8//JADBw5QVlZGWVkZer1elIIQOByNRmO2G1hJiR1Sy8xMSE21\nWgIKsgHI1+fT1K9plQCwgYAA2wPBkiRxRa3mi06dWPHoo5Bh3117VKMoLuRcqLGicM09PHBVKERl\n0AaAWQPQtGlT7rrrrtqYi8CJMecC8vWVJZY2Sy0zMmD4cPm5FRJQkKuB5inyKrWBvBVDHMCWekAX\nSkvxUyp56uxZZh46BOnpNu1UDESFRNEyoKVRY+UISvR65ojKoPUeswZg8ODBfPLJJ6SmppKTk1P+\nEAgciakgsMOSwMLCrJaAgrwDKHApIMw3zKhf3c8PXnvNtiD1S0lJFJeVMbN5c5YcPw7Lltm0UzEQ\n1SiKhOsJNRYDABgUGMhiURm03mPWALz77rs8//zztGjRgpCQEBEDENQIlgSB7SIjw64YQKFrIVfz\nrxr1q3t5wS+/yM+tDVL/mJPD5LAw+W66aVMwFMGzcqdiIMgrCL2k53DK4RqJAQAsiYwk0NVV5APU\nc8wagNLSUvR6faWHiAEIHI0lLiC7sFMFVORWRE5pjlG/euPGclE4sD5I3dTDg42ZmfLddNOmkJZm\n006l0pyVbuSp82okBgDQ1deXHJ1OjlsIN1C9xWwimF6vZ8mSJfzwww9otVqGDBnCa6+9hqur2Y8K\nBBZjLghstwHIyJBdQDbg5uZGsUcxXz76Jf1a9qv2mIAAWLIEXnjBOvePXpLI0mgoz04ID4e9e2HD\nBilF+VQAACAASURBVJvmamB0p9H0bNaTZ3s+a9d5jKFUKGjt6cnos2d5UbiB6i0WJYL9/vvvvPji\ni8yePZuDBw/y8ssv18bcBE6EORmoQ+oAffaZTR91c3Oj1LPUbBC4uNj6IPX50lI0FVtAhodDSopd\nQWCQ4wDH047XaBxgStOmdPPxEW6geoxZA/Djjz+ya9cuRowYwWOPPcbOnTv58ccfbR7wnXfe4e67\n76ZHjx6sW7fO5vMIGhamgsB27wCKi6G0FBYtsimwqlVo0Sq1TN452ehiaqsMNLawkHZeXjc19eHh\ncPKkXUFgkJVAvyT/UmO5AAD9AwI4XlRkU/9iQd3AIhdQxa25h4eHyZotpoiJieHgwYMcOHCAkpKS\n8iqjAoG5GICfnx0nz8yUF9bLl22q15CpzURZpuRo6lEWH1jM0qFLqxzj7w8XL1o/taVXr5JcWspL\n4eFyXZ3CQlCp5DftqC0R1SiK9KJ01GXqGosDRN/4T7nn1Cl62O2jE9wOzO4AunfvzqxZs7hw4QIX\nLlxg9uzZdO3a1abBfvrpJ7p06cLIkSN5+OGHeeihh2w6j6DhYSoGYLcLyMZm8AYy1BkgYTKxytaK\noCeKilgUGXnzDtrPT44mv/KKTXM1YHBXqV5X1VgugOeNqPeSyEjiRJOoeolZA7By5Upyc3Pp06cP\nvXv35tq1a6xYscKmwa5fv05sbCxfffUVq1atYuzYsaK0tACwTAZqc9N1QwDYRndKmiqNgEsBJhOr\nbKkIWiZJKIHXkpMr6+nDw2HcOJvmasBN6UarwFYcunqoRuMAPf38eOXCBZEPUE8xawD8/f3ZsGED\nmZmZXL16lU2bNhEcHGzTYI0aNWLo0KG4u7sTFRWFp6cn165dq3TM6tWr6dmzJz179qzynqDhYsoA\nFBbCxo0OqARqI0kFSRSGFppcRCtWBLXUUCWUlNDay6tqC0hDINhOohpF8f6h92s0DjC/dWuaeniI\nQHA9xagB0Gg0TJw4kR07dpS/9vjjjzN58mR0Op1Ng/Xr148ff/wRSZJIS0ujuLiYRo0aVTpm2rRp\nxMbGEhsbKxLOnAhTQeDcXNi3T35uUyVQOySgAEevHUXnozO5iAYGyvOcONFyQxVbWIgS8Nizp/Li\n6SAlUMeQjpzOOl1jNYFADgSnqNXsvOMOEQiuhxg1AG+88QYFBQX06dOn/LVPP/2U3Nxc5tv4i/nQ\nQw9x5513Eh0dzcMPP8zKlStRGrJnBE6NqRiAVmt7khVgVyVQgGxVNpKb6QYrwcFw/bp1JatjCwtJ\nKi2tqqIJD4fNm+1WAt3R5A76tOhTozWB/Fxd8XVx4YG4OLrYrdUV1DZGDcB3333H1q1badKkSflr\nzZs3Z+PGjWzfvt3mARcvXszRo0c5duwYQ4cOtfk8goaFuRjA/v12VgLds8emBVVTpgGg0ceN2Dtp\nLzExMej1+irHNWoEeXnWlaw+UlDA440bVy2rHB4OBpeQjeUgQDYAJ9JP1GgMAKBUr2dM48YiEFwP\nMWoA3N3d8fLyqvK6v78/Hh4eNTopgfNhLhEsIMDOSqAvvCA/t3JBvZR3iWa+zdCqtEyZMoXhw4cz\nZcqUKkbAy0vepWRlWWao1Ho9xwsL2X79ulwE7tYYwNWrdpeD6Ni4I0k5STUaAwAY3aQJX1y7JgLB\n9RCjBkCpVFJYWFjl9cLCQrRabY1OSuB8mAsC250HYKMM9ELOBZr7NadwRCH/3fVfjhw5wrZt29i7\nd2+VY0NCZCWQJYbqVFERZVB9ElXLlnLOgp0xAE9XT9yV7jUaAwDY2KEDbqJRfL3EqAF46qmnmDJl\nCsUVtnXFxcVMmTKFUaNG1crkBM6DqSBwQYGdBsAOGWhSThL5mnykphKREyKJjo5mzJgxDDBU7KxA\no0ZyHMASjhQUcIePT/VdtVq3vpm0ZqcReLD9g6x8YCW+7r415gZyc3FBJ0m82rKlCATXM4wagJde\neomAgADCwsLo3bs30dHRhIWFERQUxBtvvFGbcxQ4AcaCwHq97E6xOdG0qEheSG08wYXcCxRpi8AL\n/Dv6s3v3btauXYuLS9U/nUaNIDvbsvN+kpZGfElJVfcPyEkPfn7w8st2B4K7NOnChlMbatwNNDAw\nkDcvXxaN4usZRg2Ai4sLq1evJi4ujn/961+89tprJCYmsmrVqmp/+QUCezDmAiouvulftwlDIxiF\nwvyx1ZCUk8QHQz8g4IMAvh/1PYMGDTL6+2/YAViSBxBfUsKqdu2M3zG3aQOPPy4/tzMQnJSTVONu\noBXt2gGiUXx9w2wtoNatW9O6detamIrAmTFmAAoK5CQrm6nYCMaGRTQxO5ElB5fgGeJJcXExQUFB\nRo8NCpKbeR09CuPHw/r11RuuXK0WV4WC6efPGw+ctm4tFxeyMxB8Z9ideLt5k/NKzXbx6+TtjSui\nLlB9Q9zKC+oEGo0GV1fXKjJLhwSAc3NtcqWU6cu4mHuRs9fOov6bmqKiIpPHe3vLiz+YzgOILSwk\n1FyUuE0buHTJ7hhA68DWFGoKuffze2tUCqpQKNADLzRrJuSg9QhhAAR1Ao1Gw9KlS6vILB2yAxg0\nSH5upSsltTAVpYuSXFUu2sbaSoKI6mjVSlYAgek8gCOFhaRrNKbLKBt2AHYaAIVCgb+HP0fTjtZo\nDADgoUaN+DAtTcQB6hHCAAjqBNnZ2fz2229VZJaFhXL81uZCcHZKQHs174U0T+Kuk3eZ3QE0agSj\nR5vPA1ifno4euP/UKeMuoDZt4Lff7A4CA+glPUWaohqNAQBs6dgRgB2iLES9QRgAQZ3Aw8ODYcOG\nVZFZ5uVBUpIdheDskIBeyL1AZHAkAL6+vhYZgNxc03kAkiRxQaXivx07olQoqiqADLRufXM7Yacc\ndPH9i3mo3UM1KgUF8HV1pY2nJ4+cOSN2AfUEYQAEdQKtVstbb71VRWaZnn6zMKZNheDsqASalJNE\nE+8mDNs8DNdAV7MuIEtkoJdVKrxdXJiQkGA6c7ZVKzkbeO5cG2ZemTub3sm+K/tqXAoK8ErLligQ\naqD6gjAAgjqBVqvF09OzisxSobDMr26U1FTYvdumOSVmJ/LXtb84knqEC2EXLNoBmDMA08+fR6XX\nV6//r4inp2y4Jk2yfuK3ENUoiiJNUY1LQQEebtQIBbIaSBSHq/sIAyCoExhLBCsogNmz7SgEFx8P\n69bZ5EdPuJ5AWmEauapccr1zze4AQkJkA2CqJ8BPOTk816yZZXfH7dtDYqLdgWCli5IBrQaweeTm\nGncDNfPwQA8sFV3C6gXCAAjqBMbyAAoL5Vr7NhWC02pBI1fztNaPri3TcjH3IgefOYg0T2KCboLZ\nHYC/v2ykxo83HrMIdHVldXq6ZYXToqJg+XKHBIJ7h/dmxZEVteIG6uPvz5wLF0QcoB4gDICgTmDK\nANicB5CebrMCKDk3mTDfMEZsG0FGUYZFQWCFQu4LsGWL/POtMYvrGg1FZWVYnJMcFQURcm9fewPB\nvcN7k5idWCtuoE+jopCAn0UcoM4jDICgTlAjmcApKdC8uU0LZ8L1BFwULuV3zAYDYKwfgIEWLYw3\nr9mfn49Wkkzr/ysSFSW7gOzMBgbZAJTpyxgSOYSvn/ja7vOZorOPD41dXRl86pTYBdRxhAEQ1Am0\nWm21MQC7dgApKXJtfRtIuJ6Aukxdfsfs5eXFDz/8YLQfgIHmzeG///3/9s48qqk7X+CfG0BsWFTE\njU0EpNaKAoKiosa9HautTm0ZW4/WpWpfW+urXc55Z07tGZ2OY51pp/NsX21dOvWNtpVpba1HXtWo\nVbEFZZlxFxWwDLggASOQkPv+uAaCQkhCSJD8Pv9wb3Jz7++XS37f+92b9lkcqqhgZGBg0xVAm8Is\nAJxAT7+eSJJEZnFmm5uAAF6PiABENFB7RwgAgduRZblZJ3BlZSs1gLAwh2sAva15u76dYnFxMQUF\nBVb7AYByudLSpn0WW0tL+amysuUIIDMREUoegxN8AAC+Xr7oanRsytnUpo5ggN/06iWige4DhAAQ\nuJ26ujpUKlWT/aFb1QuguBhOnnRoAc0tzWVD9ob6hXLo0KGEhIRY7QcAigZQXHxvFFCV0UipwcB3\n9mTJqlSKNxmc0htgzcQ1xHaPRUJqcy0g9E400Mq+fUU0UDtGCACB22mpG5jDpSCKi2H+fGXbjgVU\nlmXyS/M5e+Ns/UIZGBhIVFSU1X4AAH36wJdf3hsFlKnT4adSMcXe5ukDB8KMGa1e/AE0kRoull90\niSMYYHK3bqwUPQLaNUIACNyONQFQUaH4QB0qBXHlimKTsdOJerniMpIkcbP6Zv1C6e/vj16vt9oP\nQDkOCgqUbcsooB8rKqg2mdg/ZIh9T8RDhjjNDBTZNZI+AX04OO9gm+cDAHzx8MOAUpDuz+Z0bkG7\nQggAgdtpzgEMigkoPV3ZtrsUhIM+gLzSPEaGj2RK9JT6iBk/P78Ww0ABBg1q2LaMAtp38yaPBwfb\n7gA2Ex+vdAiDVpuBJEliXOQ4Vh5Y6ZJ8gEBvbyTAJMu2h74KXIoQAAK305wGIMtKJI33nbZFdpWC\nMBqVJ+dPP7V7PPml+VTWVDZaJP39/VvMBAbFB6BWK5c3RwFVGo1k6nR8d/267Q5gM0OGQG4uOKkN\n68SoiRwvOU55dTk/Fv7olHNaI61nT8zxUsIM1P4QAkDgdqy1g+zcWVlI7S4FUVSkfHjVKrtNJ/ll\nigCwtJXbkggGSsSSJCm+iwsXFCF24OZNjPbE/1sSEqI0Rl682Cl+gCnRU9Ab9HTr3I3UiNRWn68l\n/nfgQCQgJTCQdcXFQgi0M4QAELgda1nAgYFKOKXdpSAKCiAxUdm203SSV5rHF7O+qA8BBdtNQJKk\naAFz5zb4LX5bcAloof6/tRPGx8NrrznFD9DDrwedvDpRXl3uknBQgEeDgjis0/FUcLDICWhnCAEg\ncDvN5QCUlyt1gBziwgWljIKdDuBqYzUF5QW89n+vNVocfX19qauro6ampsVs4J49YedOZftvf4Oc\nW1V82L+/9fr/1hgypMG54IRw0NdGvkZi70SXhIMC/GPQILyBL65dEzkB7QwhAARux2AwNKkBlJU5\nUP7ZTEEBREc75AAO8A0g65esRoujJEmo1WoWLFjQYjZwTIxlCWsZbwmWnDvn+OKXnAzHjsHYsY59\n/i6m9p/KyasnXaYFdFKpkIF+vr5kV1UJM1A7QggAgdtpzgR09aryNO0QZg3ATgHw85Wf8ZK8moyV\n79SpEzt27GgxG3jAAHjpJcVv8YczxZjA/vBPS0aNgowMOHDAKRnBiX0SCQkMYfag2S7TAhaHhHCx\npoZ+nTsLX0A7QggAgduxJgBapQE40FP3wOUDBKuDKXm1pN7+byYoKIhf/epXLWYDDxgA584pl333\nbAkyDtr/zYSHN3wRTjABSZLEUwOfYt/FfS6LBvrv2FiGBQRwsbpa+ALaEUIACNxOcz6AsjKlyYrd\nWcCyrGgAf/hDw76Ni+begr0U6YqafCr29/fnzTffbDEbeMAApQ/Nb+aY+LfqNuN2DXbc/m9m1CiY\nPt0pkUAATw96Gr3RddFAAN8PHowXwhfQnhACQOB2mtMASkthzx4HsoDLy5VFPyjILiewrkbHzZqb\n6Gp0TZZK6NKlCzqdrsVs4KgoJQl5e0ABeMnseySPBUEOPv2bSU2FwkKnFYYb0msIhjqDS6OBut8R\n8mGdOglfQDtBCACB2zFnAt8dXXPlCvz8s7JtVxbwhQuKA/jtt+1aKLN/ySaxT2KjDGBLwsLCKLah\npIGPj3J5ZlyBV4dAncT6uFY8/QOMHq3UxQBFqDnBDPTOhHeI6RZDVW0Vb+1vfc8BW3gxNJTi2lok\nEOUh2gFCAAjcTnV1Nfn5+fdE1+h0zTdXsUp+vpKKa+fT8qHCQ8iy3GyZhL59+1JYWIjJZGoxFDRy\nsAGu+cKfcwmq8LOvhlFTDBqktLhMTnaaFjBnyBwulF/AaDLy3dnvWjlA23ivf39UgAyYENnB7kYI\nAIHbycvLo6io6J7omqtX4aefHMgCzs1VGvOCXfb/3ed2c/ra6WarZUZERHDp0iUWLlzYYihodfx1\nCK6B5UO40eWWfTWMmkKS4LHH4KmnlH0nOIODHggiWB0MgN6od4kZCOCVO016unt7i4ggNyMEgMDt\nREVF0a9fv3uia8rKlKxau7OA8/KU5Ck77P+3DbfJKsnCx8uH5SnL74kAAkUDyM3NZdu2bS2Ggh5P\nvgBeMqzLRfo6lOBgO8bfHFOnwl//qmw7QQMA2D93P2pvNbIsk/JJikuEwLqYGJ7v04frRqMwBbkZ\nIQAEbsdoNKLRaBpF19TVKb7c7t3tPJksKxpARoZdC+SRoiM84P1AoxLQdxMREUFlZSVpaWlWQ0Hn\nnjpFhZcBn8xgVEjIH8bYX8q6KcaPb3CEOMEPAPBwz4fp7NOZipoKCisKXeYL+J8HHxSmoHaAEAAC\nt2MwGPD19W0UXXPjhlIH6Px5O0NAi4uhthbWrbPrKXnvxb3MT5jfrAMYFAFQVFTEhg0brIaCfl5a\nykC1GsPIa5i+UqJ/7C5l3RRqNUyb1uDgdpIW8MOcH1ChQkZmS+4Wl5mClt8xBXkD64qLScrKcsl1\nBQ0IASBwOzqdDr+74sJLSpQnZrtDQHNzlZh5sMtO/vXpr/l31b+t1sn38/NDrVZz7dq1Zs9z9U6E\ny0m9nqBr/qg+VqJ/vLxwjhlo/nyn9Qcwk9AngWA/ZXBGk9FlpqB3Y2J4JSwMI+CnUonQUDfgFgFQ\nVlZGeHg4p0+fdsflBe2MoqIiIiIiGr124YISBQR2Pj2fOKHY/+2om/Ovsn9xs/omey7sabFdYnh4\nOIsWLWrWCTz9n/+kDngqOBhd8C1mzYJnn1UEmFPMQOPGweXLyraTNACAE4tP4K3ypk6uc6kp6M8x\nMcT7+3PLZMILoQm4GpcLAIPBwOLFi3nggQdcfWlBO6WwsPAeAWA0Kusb2BkCunevkkBgR92cL09+\nSW//3khIzTqAzfj5+bFnz54mncA3DAYydToklGzXuQGh7N8Pn3+uvO8UM5BKBb/9bcO+nbkOzRES\nEMKn0z/FW+XtclPQiaQkXggJoQ7wkSSyq6qEEHARLhcAK1asYMmSJYSEhLj60oJ2SlMC4OJF+I//\nsDMEtLISsrOhb19l3wYTiSzLfJ73OedunLOpWXpiYiLx8fEkJyczduxYUlMbyihMzM0FQEL5YX2S\nFEN4eIMgA+jWzca5WGPp0gYzEDhNCMwZPIeunZX627V1tS7TAkCpFTTEzw+DLOMFwhzkIlwqADZv\n3kyPHj2YMmVKs8d8/PHHJCUlkZSUxNWrV104OoG7aEoAZGcrOU9nzza0hGwRrVYxtL/zjrJvw6K4\n58Ieym+XU22s5vmhz1t9+geIj48nIiICjUbDgQMHeP755zGZTMw9dYoTVVVEd+6MiQYH56xZjZ3Y\njz2maDetQq2G3/8eIiOVfSdFBEmSxOVXLtdrAZtzNjNuyziXaQI5yckk+vtjtpKtKy5GpdUKQdCG\nSLJsV4xFqxgzZgySJCFJEjk5OcTGxrJz50569+7d5PFJSUlkCVWwQ1NRUUFYWBg6nQ7J4lE5OlpJ\nft25U8np2rSpISu4WV58Ec6cgR9+sDl0aPTG0WReycRoMhISEMKV/7RepbK0tJTo6GhMJhOZmZkk\nJyfzzXff8ahFosKrYWH1hd9u3Lg3lPXJJ2HbNhvmYw2jUdF0fvml4TUnCYLhnwxvpAmFBYbx86Kf\n6e3f9O/U2bx6/jx/uis3YKi/P1lJSS65fkfA1rXTpRrAwYMHOXDgAFqtlvj4eD777LNmF3+BZ2B+\n+pdlub68Qnm50s/dsqtWi7bz27dh82Zl8QebFsL0k+kcLT6K+RmoX9d+LX6mV69exMXFkZqaSkJC\nAkajkbS8PED5MamgUdXPoCBFLlny1VdKklur8PZW/B2WOMkUdGzhMYaFDqvfL9YV0+/9fuSV5rX6\n3LawLiaGof7+9fsSiklI0mrxFhqBUxFhoAK3UlhYSHh4eKPyCllZJhISGkw/NjmBt2xRImTAJtu/\n3qBn3jfzqJPrqJPrGBU+qkXzj5knnngCX19fAEzr11ORmIifLDcy/ViyatW94583zwmmoAEDFPXI\n0kam1bbypArHFh5jZNjI+v1qYzUpn6QQtCbIJYIgKykJWaNhqL8/lrpcHYppSNJqhaPYCbhNAGi1\nWgYMGOCuywvaCYWFhfj4+DQqr7B9+0FKSpQF8tlnFd+uVfR6WLtWSQKDFqN/DHUGZn0xi1sGpUNX\nSECIzYs/QFpaGtqsLEzr10NsLOj13AJmmUz8MSrqnuO7dFG0GEsyMpRI1bw8pXeAw4bYzz5TTmbG\nSV3DAA4vOIz8loyfj+Jwvm28TXl1OQkfJZC8Idkl/oGspKRG2oAlZq1A0mqFr8BBXOoDsJfW+ABq\na40sXbqT5OSeHDtWyvDhvTh6tARJkkhJ6W31NXuPd8Y5PPX4Dz74Fm9vb9TqUrKyrjF58sNotavQ\n6xv8AVZ7A9fWwsyZsH9/Q7iQFVv4kaIjPLr1UapqqpAkye6n/+fPnGFDSUnDiq3XK07ZM2dg6VIm\nTZrEhAkTUKlUDB48mPz8fOLi4sjPz+fKlSTee0/T5HkjIy+xYkUxWVmOfZ+9//onnvjXPxqdU+/9\nACenpFFUVMk3CU+36n69dXkRZZ3v5O1YrBi+shqjyYQkSYR0iqC2yptXwldzOsvo9P+fK9N68EM3\n052wKnuWLYmHT97mRa+Advf/39R7WVllrF8/nU6dbI1+uBdb184OKQD0eiN+2w5AP2HhEtiAZZym\nPZh/OncWf9uoBeypbGc7b7GSlbzd7Ptaxt75q0GDtv6vra9p0bD2zVXoO1tks929esigvqVG8qrB\n/2Ywuq5XkSQIKO9BVddrzb5m7b2mXrs+fCbGfk+ieAgsB+PgvWxvFJj4L90lVi1f4NDHPVoAvPHG\nFv74SARKmSkhBARWcHTxR1lqBuTlcWrZMjs+5X3nkzW0xWL1FivRoEXDAaef25LU+XA4opk33bGi\nJH4IAQ9aOeA+Eg5yHUhe/PhEGqMqSh06hUcLAKPRiM+uTRDYyi5Mgo6Ng4u/ZUiiyWRi3759rF27\nlgxLW3yLdEZ5QKmiLTQCsyDQorGqFTgTq0LBnSS0JBzaIbrzGKY+h7fNSTCN8WgBAKDTGSkfMgBd\n4kgCso9QOXQk/llHAKhKsv6avcc74xyefrx+WCoPllykZsQ4fI9qlbVZo8HrkBY0GiW6xfwX7n3N\nSdEvjmLuEmY0GsnJyUGSpHt8AHFxceTm5jZ676GHBpOd/U9kWSI2djDnzuXTv39c/d+TJ/ORJImH\nHhpk9TVr7/X6cC39Lp+/r/4fPPn4wONH6JZ7msBA4QMQiWACgUBgJ+0yEUwgEAgE7Yd2rQEEBwcT\naa53YidXr16lh80lJDsWYu6eOXfw7PmLuTfM/dKlS1b7Vphp1wKgNXiy+UjM3TPnDp49fzF3++cu\nTEACgUDgoQgBIBAIBB6K18qVTuop1w4ZOnSou4fgNsTcPRdPnr+Yu310WB+AQCAQCKwjTEACgUDg\noXQ4AWAymViyZAkjRoxAo9Fw3sNKxCYmJqLRaNBoNDz33HPuHo5LOHbsGBqNBoDz58+TmprK6NGj\nWbp0KSaTyb2Da2Ms537ixAlCQ0Pr7//27dvdO7g2wmAwMGfOHEaPHs2wYcPYuXOnx9z3pubeqvsu\ndzB27Nghz507V5ZlWT569Kg8ffp09w7Ihdy+fVuOj4939zBcypo1a+RBgwbJw4cPl2VZlqdNmybv\n379flmVZXrx4sZyenu7G0bUtd899w4YN8rvvvuvmUbU9GzdulJctWybLsixfv35dDg8P95j73tTc\nW3PfO5wG8OOPP/LII48AkJKS4lFxwbm5uej1eiZPnsz48ePJzMx095DanOjoaNLT0+v3s7OzGTtW\nKXv86KOP8oO5RWQHpKm579q1izFjxrBgwQIqW+ykc38ya9Ysfve73wEgyzLe3t4ec9+bm7uj973D\nCQCdTkeXLl3q9728vDC2uvfe/YFarWbFihXs2bOHjz76iGeeeabDz/3Xv/41PhYN2WVZrm8uHxAQ\nQEVFhbuG1ubcPfdhw4axdu1aDh48SFRUFG+/7ZoqoK7G39+fgIAAKisrefLJJ1m1apXH3Pem5t6a\n+97hBEBgYGAjCWgymRwuqXq/ERsby7PPPoskScTGxtK9e3dKSkrcPSyXolI1/EtXVlbStWtXN47G\ntcyYMaM+FHDGjBmcOHHCzSNqO4qKihg3bhxz5sxh9uzZHnXf7557a+57hxMAo0aN4vvvvwcgMzOT\nuLg4N4/IdWzcuJFXX30VgF9++QWdTkefPn3cPCrXkpCQgPZOaejdu3czevRo9w7IhUyZMoWffvoJ\ngL1793bYmPjS0lImT57MmjVrmD9/PuA5972pubfmvne4PACTycQLL7xAXl4esiyzadMmj2k+X1tb\ny7x58ygsLESSJNasWcPIkSPdPaw259KlS6SlpZGZmcnZs2dZtGgRtbW1PPTQQ2zYsAEvLy93D7HN\nsJz78ePHeemll/Dx8aF37958/PHHBAYGunuITmfZsmVs37690e/6/fff5+WXX+7w972pua9evZrX\nX3/dofve4QSAQCAQCGyjw5mABAKBQGAbQgAIBAKBhyIEgEAgEHgoQgAIBAKBhyIEgEAgEHgoQgAI\nBAKBhyIEgOC+5eWXXyY+Pp74+Hg6derEgw8+WL9/+/Zt4uPjuXnzZpuPQ6PR8NVXX7X5dQQCZ+MZ\nNRIEHZK//OUv9duRkZFs3bqVpKSk+tdycnLcMSyB4L5BaACCDoskSVy7do3Nmzczbdo0Jk6cSExM\nDOPHj2fHjh2MGzeO0NBQ1q1bV/+ZTz/9lKFDh5KQkMDEiRM5ffq0Tdf65ptvGDZsGH379mXBdfcC\ndAAAAtFJREFUggX19ei//vprEhISGDx4MKmpqfUp+ytXruTFF1+s/7zlvkajYebMmQwcOJAPPviA\n9PR0EhMTSUpKYvjw4Rw8eNBZX5HAwxEagMAjOHToEPn5+YSGhhIXF8e2bdvYu3cv+fn5pKSksHz5\ncg4dOsSWLVs4dOgQarWajIwMZs6cycmTJ1s8f2VlJUePHqWmpoaYmBgOHz5Mjx49WLJkCUeOHCEq\nKop9+/bx+OOPc+bMmRbP161bt/rrRkdHs3XrVlJSUsjIyECr1TJmzJhWfycCgRAAAo8gOTmZ8PBw\nAPr168fkyZNRqVRER0dTXV2NXq9n165dnD9/vlH9pBs3bnDjxg2CgoKsnv/pp5/Gy8sLtVpN//79\nKSsrIz8/nwkTJhAVFQXA+PHj6dmzJ9nZ2S2O17KYWVpaGjNmzGDq1KlMmjSJ119/3ZGvQCC4B2EC\nEngEvr6+jfYt6+ibqaurY86cOeTk5JCTk8Px48fJysqiW7duLZ7f8nySJCHLcpNtCU0mEwaDof4Y\nM7W1tY2O8/f3r99evXo1hw8fJikpic2bNzNixIgO2/JQ4FqEABAI7jB58mT+/ve/1/dQ+Oijj5gw\nYYLD5xs/fjwZGRkUFBQAsG/fPoqKihg+fDg9evQgOzsbWZa5desWGRkZTZ7DaDQSGRnJrVu3WLJk\nCevXr+fUqVMYDAaHxyUQmBEmIIHgDlOmTOGNN95g0qRJqFQqAgMDSU9PR5IksrKyWLhwoV2RRQMH\nDmT9+vXMnDkTo9GIWq3m22+/pUuXLjzzzDPs3r2b/v37ExoayogRI2iqMK+3tzfvvfces2fPxsfH\nB5VKxcaNG+/RaAQCRxDloAUCG5k6dSq7du1y9zAEAqchTEACgQ1cuXKFF154wd3DEAicitAABAKB\nwEMRGoBAIBB4KEIACAQCgYciBIBAIBB4KEIACAQCgYciBIBAIBB4KP8P5LU8C8nMWC4AAAAASUVO\nRK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x110fa0eb8>"
      ]
     },
     "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((ts1_mt3d[:,0])/3600, ts1_mt3d[:,1], 'kD', markersize=2.0, mfc='none',mec='k')\n",
    "ax.plot((ts2_mt3d[:,0])/3600, ts2_mt3d[:,1], 'b*', markersize=3.0, mfc='none',mec='b')\n",
    "ax.plot((ts3_mt3d[:,0])/3600, ts3_mt3d[:,1], 'r+', markersize=3.0)\n",
    "ax.plot((ts4_mt3d[:,0])/3600, ts4_mt3d[:,1], 'g^', markersize=2.0, mfc='none',mec='g')\n",
    "ax.plot((ts5_mt3d[:,0])/3600, ts5_mt3d[:,1], 'co', markersize=2.0, mfc='none',mec='c')\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', 'Otis, Site 2', 'Otis, Site 3', 'Otis, Site 4', 'Otis, Site 5',\n",
    "                 'MT3D-USGS, Site 1', 'MT3D-USGS, Site 2', 'MT3D-USGS, Site 3', 'MT3D-USGS, Site 4', 'MT3D-USGS, Site 5',\n",
    "                 ), \n",
    "                loc='upper right', labelspacing=0.25, columnspacing=1,\n",
    "                handletextpad=0.5, handlelength=2.0, numpoints=1)\n",
    "leg._drawFrame = False \n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.6.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}