02-Code-Build-and-Test.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2>Funwave-TVD</h2>\n",
    "Funwave-TVD is a code to simulate the shallow water and Boussinesq equations written by Dr. Fengyan Shi. It is capable of running a variety of complex ocean simulations. For more information, see Dr. Fengyan Shi's website: https://fengyanshi.github.io/shi_website/_build/html/funwave.html\n",
    "\n",
    "We will use it as our example scientific code for this tutorial because it is reasonably complex, real scientific code but can be compiled and run quickly."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/home/jovyan/agave\n",
      "Collecting setvar\n",
      "  Downloading https://files.pythonhosted.org/packages/69/fb/e93e07768694deb9b1cc7fdf7217775753b6b42a3e10a9fa7c2c2f24ceab/setvar-1.1.2.tar.gz\n",
      "Building wheels for collected packages: setvar\n",
      "  Running setup.py bdist_wheel for setvar ... \u001b[?25ldone\n",
      "\u001b[?25h  Stored in directory: /home/jovyan/.cache/pip/wheels/c7/cc/6f/fa8b8dcd762a31d58ba0c3dd0992200b5a32eddecd18fc7d83\n",
      "Successfully built setvar\n",
      "Installing collected packages: setvar\n",
      "Successfully installed setvar-1.1.2\n",
      "\u001b[33mYou are using pip version 9.0.1, however version 10.0.0 is available.\n",
      "You should consider upgrading via the 'pip install --upgrade pip' command.\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "!mkdir -p ~/agave\n",
    "\n",
    "%cd ~/agave\n",
    "\n",
    "!pip install --upgrade setvar\n",
    "\n",
    "import re\n",
    "import os\n",
    "import sys\n",
    "from setvar import *\n",
    "from time import sleep\n",
    "\n",
    "# This cell enables inline plotting in the notebook\n",
    "%matplotlib inline\n",
    "\n",
    "import matplotlib\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The parameter file for Funwave-TVD is called \"input.txt\". Here we create it using writefile. You don't need to understand the details of how this file for this tutorial."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Writing file `input.txt'\n"
     ]
    }
   ],
   "source": [
    "writefile(\"input.txt\",\"\"\"\n",
    "!INPUT FILE FOR FUNWAVE_TVD\n",
    "  ! NOTE: all input parameter are capital sensitive\n",
    "  ! --------------------TITLE-------------------------------------\n",
    "  ! title only for log file\n",
    "TITLE = VESSEL\n",
    "  ! -------------------HOT START---------------------------------\n",
    "HOT_START = F\n",
    "FileNumber_HOTSTART = 1\n",
    "  ! -------------------PARALLEL INFO-----------------------------\n",
    "  ! \n",
    "  !    PX,PY - processor numbers in X and Y\n",
    "  !    NOTE: make sure consistency with mpirun -np n (px*py)\n",
    "  !    \n",
    "PX = 2\n",
    "PY = 1\n",
    "  ! --------------------DEPTH-------------------------------------\n",
    "  ! Depth types, DEPTH_TYPE=DATA: from depth file\n",
    "  !              DEPTH_TYPE=FLAT: idealized flat, need depth_flat\n",
    "  !              DEPTH_TYPE=SLOPE: idealized slope, \n",
    "  !                                 need slope,SLP starting point, Xslp\n",
    "  !                                 and depth_flat\n",
    "DEPTH_TYPE = FLAT\n",
    "DEPTH_FLAT = 10.0\n",
    "  ! -------------------PRINT---------------------------------\n",
    "  ! PRINT*,\n",
    "  ! result folder\n",
    "RESULT_FOLDER = output/\n",
    "\n",
    "  ! ------------------DIMENSION-----------------------------\n",
    "  ! global grid dimension\n",
    "Mglob = 500\n",
    "Nglob = 100\n",
    "\n",
    "  ! ----------------- TIME----------------------------------\n",
    "  ! time: total computational time/ plot time / screen interval \n",
    "  ! all in seconds\n",
    "TOTAL_TIME = 3.0\n",
    "PLOT_INTV = 1.0\n",
    "PLOT_INTV_STATION = 50000.0\n",
    "SCREEN_INTV = 1.0\n",
    "HOTSTART_INTV = 360000000000.0\n",
    "\n",
    "WAVEMAKER = INI_GAU\n",
    "AMP = 3.0\n",
    "Xc = 250.0\n",
    "Yc = 50.0\n",
    "WID = 20.0\n",
    "\n",
    "  ! -----------------GRID----------------------------------\n",
    "  ! if use spherical grid, in decimal degrees\n",
    "  ! cartesian grid sizes\n",
    "DX = 1.0\n",
    "DY = 1.0\n",
    "  ! ----------------SHIP WAKES ----------------------------\n",
    "VESSEL_FOLDER = ./\n",
    "NumVessel = 2\n",
    "  ! -----------------OUTPUT-----------------------------\n",
    "ETA = T\n",
    "U = T\n",
    "V = T\n",
    "\"\"\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here is a simple bash script to run the code...."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Writing file `run.sh'\n"
     ]
    }
   ],
   "source": [
    "writefile(\"run.sh\",\"\"\"\n",
    "#!/bin/bash\n",
    "export LD_LIBRARY_PATH=/usr/local/lib\n",
    "mkdir -p rundir\n",
    "cd ./rundir\n",
    "cp ../input.txt .\n",
    "mpirun -np 2 ~/FUNWAVE-TVD/src/funwave_vessel\n",
    "\"\"\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## FOLLOWING ALONG AT HOME  \n",
    "\n",
    "If you are following along at home using the docker-compose stack, you will need to run the following cell to get the hostname and port of your tcp tunnel so Agave can contact your system without a public IP address."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.tcp.ngrok.io:16068\n",
      "Reading file `ngrok_port.txt'\n",
      "Reading file `ngrok_host.txt'\n",
      "Reading file `ngrok_host.txt'\n",
      "Reading file `ngrok_port.txt'\n",
      "VM_PORT=16068\n",
      "VM_MACHINE=0.tcp.ngrok.io\n",
      "VM_IPADDRESS=52.15.72.79\n"
     ]
    }
   ],
   "source": [
    "if os.environ.get('USE_TUNNEL') == 'True': \n",
    "    # fetch the hostname and port of the reverse tunnel running in the sandbox \n",
    "    # so Agave can connect to our local sandbox\n",
    "    !echo $(ssh -q -o StrictHostKeyChecking=no -o UserKnownHostsFile=/dev/null sandbox 'curl -s  http://localhost:4040/api/tunnels | jq -r '.tunnels[0].public_url'') > ngrok_url.txt  \n",
    "    !cat ngrok_url.txt | sed 's|^tcp://||'\n",
    "    !cat ngrok_url.txt | sed 's|^tcp://||' | sed -r 's#(.*):(.*)#\\1#' > ngrok_host.txt\n",
    "    !cat ngrok_url.txt | sed 's|^tcp://||' | sed -r 's#(.*):(.*)#\\2#' > ngrok_port.txt\n",
    "\n",
    "    # set the environment variables otherwise set when running in a training cluster\n",
    "    os.environ['VM_PORT'] = readfile('ngrok_port.txt').strip()\n",
    "    os.environ['VM_MACHINE'] = readfile('ngrok_host.txt').strip()\n",
    "    os.environ['AGAVE_SYSTEM_HOST'] = readfile('ngrok_host.txt').strip()\n",
    "    os.environ['AGAVE_SYSTEM_PORT'] = readfile('ngrok_port.txt').strip()\n",
    "    !echo \"VM_PORT=$VM_PORT\"\n",
    "    !echo \"VM_MACHINE=$VM_MACHINE\"\n",
    "    setvar(\"VM_IPADDRESS=$(getent hosts ${VM_MACHINE}|cut -d' ' -f1)\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here's how the average developer might run the code when doing their research."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Warning: Permanently added '[52.15.72.79]:16068' (ECDSA) to the list of known hosts.\n",
      "input.txt                                     100% 1753     1.7KB/s   00:00    \n",
      "run.sh                                        100%  141     0.1KB/s   00:00    \n",
      "Warning: Permanently added '[52.15.72.79]:16068' (ECDSA) to the list of known hosts.\n",
      " BATHY_CORRECTION DOES NOT EXIST. USE DEFAULT VALUE\n",
      " INI_UVZ DOES NOT EXIST. USE DEFAULT VALUE\n",
      " WaterLevel DOES NOT EXIST. USE DEFAULT VALUE\n",
      " ETA_LIMITER DOES NOT EXIST. USE DEFAULT VALUE\n",
      " EqualEnergy DOES NOT EXIST. USE DEFAULT VALUE\n",
      " PERIODIC DOES NOT EXIST. USE DEFAULT VALUE\n",
      " DIFFUSION_SPONGE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " DIRECT_SPONGE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " FRICTION_SPONGE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " OBSTACLE_FILE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " BREAKWATER_FILE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " DISPERSION DOES NOT EXIST. USE DEFAULT VALUE\n",
      "         DISPERSION Default:  DISPERSION\n",
      " Gamma1 DOES NOT EXIST. USE DEFAULT VALUE\n",
      "        Gamma1 Default:  1.0: DISPERSION\n",
      " Gamma2 DOES NOT EXIST. USE DEFAULT VALUE\n",
      "              Gamma2 Default:  1.0: Full nonlinear\n",
      " Beta_ref DOES NOT EXIST. USE DEFAULT VALUE\n",
      "               Beta_ref Default:  -0.531\n",
      " Gamma3 DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                      Gamma3 Default:  1.0: NOT fully linear\n",
      " VISCOSITY_BREAKING DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                   VISCOSITY_BREAKING Default:  SWE Breaking\n",
      " SWE_ETA_DEP DOES NOT EXIST. USE DEFAULT VALUE\n",
      "               SWE_ETA_DEP Default:  0.8\n",
      " Friction_Matrix DOES NOT EXIST. USE DEFAULT VALUE\n",
      "             Friction_Matrix Default:  constant Cd\n",
      " Cd DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                  Cd_fixed Default:  0.0\n",
      " Time_Scheme DOES NOT EXIST. USE DEFAULT VALUE\n",
      "       Time_Scheme Default:  Runge_Kutta\n",
      " CONSTRUCTION DOES NOT EXIST. USE DEFAULT VALUE\n",
      " HIGH_ORDER DOES NOT EXIST. USE DEFAULT VALUE\n",
      "  HIGH_ORDER                     NOT DEFINED, USE FOURTH-ORDER\n",
      " CFL DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                       CFL Default:  0.5\n",
      " DT_fixed DOES NOT EXIST. USE DEFAULT VALUE\n",
      " FroudeCap DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                 FroudeCap Default:  3.0\n",
      " MinDepth DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                MinDepth Default:  0.1 m\n",
      " MinDepthFrc DOES NOT EXIST. USE DEFAULT VALUE\n",
      "             MinDepthFrc Default:  0.1 m\n",
      " SHOW_BREAKING DOES NOT EXIST. USE DEFAULT VALUE\n",
      "            SHOW_BREAKING Default:  TRUE\n",
      " Cbrk1 DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                    Cbrk1 Default:  0.65\n",
      " Cbrk2 DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                    Cbrk2 Default:  0.35\n",
      " WAVEMAKER_Cbrk DOES NOT EXIST. USE DEFAULT VALUE\n",
      "            WAVEMAKER_Cbrk Default:  1.0\n",
      " WAVEMAKER_VIS DOES NOT EXIST. USE DEFAULT VALUE\n",
      "           WAVEMAKER_VIS Default:  FALSE\n",
      " T_INTV_mean DOES NOT EXIST. USE DEFAULT VALUE\n",
      "             T_INTV_mean Default:  LARGE\n",
      " STEADY_TIME DOES NOT EXIST. USE DEFAULT VALUE\n",
      "             STEADY_TIME Default:  LARGE\n",
      " C_smg DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                     C_smg Default:  0.0\n",
      " nu_bkg DOES NOT EXIST. USE DEFAULT VALUE\n",
      "                    nu_bkg Default:  0.0\n",
      " FIELD_IO_TYPE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " NumberStations DOES NOT EXIST. USE DEFAULT VALUE\n",
      " OUTPUT_RES DOES NOT EXIST. USE DEFAULT VALUE\n",
      "         OUTPUT_RES NOT FOUND, OUTPUT_RES=1: full resolution\n",
      " DEPTH_OUT DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Hmax DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Hmin DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Umax DOES NOT EXIST. USE DEFAULT VALUE\n",
      " MFmax DOES NOT EXIST. USE DEFAULT VALUE\n",
      " VORmax DOES NOT EXIST. USE DEFAULT VALUE\n",
      " MASK DOES NOT EXIST. USE DEFAULT VALUE\n",
      " MASK9 DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Umean DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Vmean DOES NOT EXIST. USE DEFAULT VALUE\n",
      " ETAmean DOES NOT EXIST. USE DEFAULT VALUE\n",
      " WaveHeight DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SXL DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SXR DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SYL DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SYR DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SourceX DOES NOT EXIST. USE DEFAULT VALUE\n",
      " SourceY DOES NOT EXIST. USE DEFAULT VALUE\n",
      " P DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Q DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Fx DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Fy DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Gx DOES NOT EXIST. USE DEFAULT VALUE\n",
      " Gy DOES NOT EXIST. USE DEFAULT VALUE\n",
      " AGE DOES NOT EXIST. USE DEFAULT VALUE\n",
      " OUT_NU DOES NOT EXIST. USE DEFAULT VALUE\n",
      " TMP DOES NOT EXIST. USE DEFAULT VALUE\n",
      "      EtaBlowVal Default:  100xmax_depth\n",
      " ----------------- STATISTICS ----------------\n",
      "  TIME        DT\n",
      "  0.4428E-01  0.4428E-01\n",
      "  MassVolume  Energy      MaxEta      MinEta      Max U       Max V \n",
      "  0.1937E+04  0.1247E+08  0.2996E+01 -0.4596E-06  0.4050E-01  0.4050E-01\n",
      "  MaxTotalU   PhaseS      Froude \n",
      "  0.4050E-01  0.1066E+02  0.3800E-02\n",
      "   PRINTING FILE NO.     1  TIME/TOTAL:        0.044 /       3.000\n",
      " ----------------- STATISTICS ----------------\n",
      "  TIME        DT\n",
      "  0.1017E+01  0.4414E-01\n",
      "  MassVolume  Energy      MaxEta      MinEta      Max U       Max V \n",
      "  0.1927E+04  0.1246E+08  0.1900E+01 -0.1727E-05  0.7329E+00  0.7325E+00\n",
      "  MaxTotalU   PhaseS      Froude \n",
      "  0.7338E+00  0.1050E+02  0.6987E-01\n",
      "   PRINTING FILE NO.     2  TIME/TOTAL:        1.017 /       3.000\n",
      " ----------------- STATISTICS ----------------\n",
      "  TIME        DT\n",
      "  0.2037E+01  0.4473E-01\n",
      "  MassVolume  Energy      MaxEta      MinEta      Max U       Max V \n",
      "  0.1906E+04  0.1246E+08  0.7393E+00 -0.2109E-05  0.9186E+00  0.9148E+00\n",
      "  MaxTotalU   PhaseS      Froude \n",
      "  0.9215E+00  0.1024E+02  0.8998E-01\n",
      "   PRINTING FILE NO.     3  TIME/TOTAL:        2.037 /       3.000\n",
      " ----------------- STATISTICS ----------------\n",
      "  TIME        DT\n",
      "  0.3032E+01  0.4566E-01\n",
      "  MassVolume  Energy      MaxEta      MinEta      Max U       Max V \n",
      "  0.1890E+04  0.1246E+08  0.7188E+00 -0.9413E+00  0.7352E+00  0.6932E+00\n",
      "  MaxTotalU   PhaseS      Froude \n",
      "  0.7378E+00  0.1019E+02  0.7239E-01\n",
      "   PRINTING FILE NO.     4  TIME/TOTAL:        3.032 /       3.000\n",
      " Simulation takes   13.767519200002425      seconds\n",
      " Normal Termination!\n"
     ]
    }
   ],
   "source": [
    "!scp -o \"StrictHostKeyChecking=no\" -P $VM_PORT input.txt run.sh $VM_IPADDRESS:.\n",
    "!ssh -p $VM_PORT $VM_IPADDRESS bash run.sh"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here's how they might retrieve the data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Warning: Permanently added '[52.15.72.79]:16068' (ECDSA) to the list of known hosts.\r",
      "\r\n"
     ]
    }
   ],
   "source": [
    "!ssh -p $VM_PORT $VM_IPADDRESS tar -C rundir -czf output.tgz output\n",
    "!rm -fr output.tgz output\n",
    "!scp -q -r -P $VM_PORT $VM_IPADDRESS:output.tgz .\n",
    "!tar xzf output.tgz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "eta_00001  eta_00003  eta_00005  u_00002  u_00004  v_00001  v_00003  v_00005\r\n",
      "eta_00002  eta_00004  u_00001\t u_00003  u_00005  v_00002  v_00004\r\n"
     ]
    }
   ],
   "source": [
    "!ls output"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now let's use Python to graph the output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsMAAACoCAYAAAAFIOPyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnV+MLdlZ3de365y+F9sQxgKNJp5xxg9WJAcpEFkExVJk\n4aA4AWEeopGJQANxNC9AcESEx7ygvM1DhOAheRgZkolAwRNAsoVQCDJYKFJkGAwRsR0nFth4nLEH\nEhPA4t4+p/aXh/3v27t2nVPdfU5333vWT+rp2lW7qnZVH91ZvXrtb4uqghBCCCGEkFPE3fQACCGE\nEEIIuSkohgkhhBBCyMlCMUwIIYQQQk4WimFCCCGEEHKyUAwTQgghhJCThWKYEEIIIYScLFcSwyLy\nThH5tIh8RkSePdSgCCGEEEIIuQ7ksnWGRWQA8D8BfBuAlwH8NoDvVtVPHm54hBBCCCGEHI+rOMPf\nDOAzqvoHqnoO4OcBvOswwyKEEEIIIeT4rK5w7hsAfN60Xwbwt3edMLzutbp6/euvcEtCCCGEEEL2\nc/75l/9EVb9+X7+riOFFiMgzAJ4BgOGRR/BXf+S9x74lIeQUkYt11339ZUeEbNe5c8cuPL75+8ve\nwbcXu8Sxnefsvv+uV3fhexFCyCX57Hv/xeeW9LtKTOILAJ4w7cfjvgpVfV5V36qqbx1e99or3I4Q\nQgghhJDDchVn+LcBvFlE3oQggt8N4B8fZFSEENLjAoboTvOyZ122/fe0J87tvvMvMpY96C4rdYlr\n3J7ePkp7jUl/3Xk8nz4zlskj7xoyXWNCyJG5tBhW1a2I/CCAXwUwAPgZVf3EwUZGCCELBeWs/ruo\nYDXtSuxKv0/3PnPX3NVnbt9StNOw++wLmhPCpk8ltrX5jkYsK+qx2+uLdu4nsz+v7u8FbV+KY0LI\ngblSZlhVfwXArxxoLIQQElhibu4Tpb3r7BC1Krpb9O4Sx+2x6ng4Jj1R3LmPXNAp1r0urqAqoZn6\ntyJ37rvEc6z+1VosV+JYOmPYNUArwuPmzlcwJ7wJIeSScAU6QgghhBByshy9mgQhhCzism4wsNwR\n3hWD2OcEW2e34wxPnN/kcjb9pXGSJ7e9ojOcz477VTXcQ03/JkKhih0OsYSMcDWs2gnOTnF+9Jkf\nVC9SYZ83jdlefolLTIeYEHIFKIYJITfPMWIRne3FAjgda9umjzTiWExfK3gFSRAbAWz722sBlxbD\nVXIhCdy8Q0o/aBHA6cQkeHNfe675isdhxa+NUMSIhM0cd4Wx1v1LZ538oFUWzDGkKCaEXAGKYULI\nzbFHBF/YCZ4Rw9q4tV0R3HN/TXsifp0RxaKlaxS+qb/EbckiWjtieN4pds2z+hkn2DrEVgxrFMJ1\nG8Y5liyQQ1sBb9xjRdW/colb1zgJYyNO9wrj9mGu6hRTEBNCLggzw4QQQggh5GShM0wIuRku4wpf\nNhbRdY613t9pZzfY6SQGYZ1eiMIZpzi4v6XtnM+XdvGYdYNd4xS3bnAP3zjBfhKTkNxPAXjv4rH6\nuCrgvdROsUlC5BiFL82JSzzXbn6GCg2VJ2yswbq5h4hN0B0mhFwQimFCyPVzLCHcRiMuEotw5R45\n+uCaHDCK0C0xiSCEk5hN2+nSg/Ohjzl/cD4PxcW+7gpiOAleu3/M4jcIZY33HL0LejWJZS9wToIg\njvvVW7EclXEak4lQZGHsrZoFSgeZF8Spu2JeHAOzgjgd6sIMMSHkAlAME0Kuj+vMCEt73FwjCd/W\nCa7Er0JikKxyeqPYTU5wEr9JwCbxO6TjUfxaZ9h1BLGDcZqxWxBPXGFItW/0Dj5ePwnlJI69eoxG\n7I7i4DUIYiCIYy9GDHupnGRAs2ssCqgHstJvnWFvRHF6zyg54uwSt85wK2aluQ5Kky4xIeSqMDNM\nCCGEEEJOFjrDhJDrYY8r3D/nChnhXj5Y+m1xOnWGnakQ4fzUCW6c38HEJgbnMcSxrKJTvIpOrYPW\nbeMKh7a/UEzCa3R8UTLCW+9KDAKxHe+39Q6jk+wUD04x+tIWEYhxhr3ECEU8jgE5P9zWJFbfWLou\nRSr6sYlDZIj3QneYELIHimFCyK1gr8a5TEZ4T93gIn6DEE6xCMQogxgB7KKYTBPikth1UQin2MMQ\nhW5qryRsJ4G7Eg8nfiKGXY5HlOvMMVZ54RBzaMVwEslbdVhX4tdh9A7beP1RBSIuxzO8dxhF8oS7\nEBFx8HFc6iXr4vBD0xCVACDOCmKgzjzE9pwgBpZliBsW1SEmhJAdUAwTQo7PISbM7bhmNyPc9m3E\ncnKDARQhbPK64mwu2FdO8GoYsxOcXOAkbgfnsRKP9TCW/jLm46vo/K5dPB6FsBXBboGVmcTvqJIF\ncdq/8UNub9Vh6x22bgAArFWwGQcMLijaJIxHlybeKZwXbNN9VBCs4ORAAynu7AGIK2pV/R5B3D7W\nvgxxD06oI4QcGGaGCSGEEELIyUJnmBBy+1hSPQImIzx3jYUZYXHTihHD4Ctn2MYi1oOvYhHrYcRK\nijN85kasovO7Sk5xaseYxFpSTMJjLaOJSUTHGaXaRMtofAyvoZLERofYdti4Oiax8QO2Gu6/9Q4r\n8Tj3of8oDoNz2IxDvr8Tc/0Ymyil5STf3/lotrpivSZ3GOhkiCXavjOZGBUTmQAuXF2CEEIuA8Uw\nIeR4XLaUWu/8ue3U3jlhLn3NZ4Rt7eAkhJ3NAQ8mBuFK5nc9jBjE42xIYnfEmRHHazfizG2z+F25\nEQPM+RLOT9EIu50YxGPU+g95JSbh4CH5+EaHMEkuCtatH7BxDuc+/HO/kQFb5+HGOKHODZBxyCJ8\nMw7VqxwFwOhKNgIOaQbdCAd4o1FdnSEOF2rELYDZOsS9csW7JtQ1sNQaIeQyUAwTQm4PF5gJtXtR\nDTTiF5NqEfb4NCMchHDJAYecMFCc4LPo9CYn+GwIKduVeNwZtlkM33FbrNyItYy5PUQ3GEB0hT0G\nlAl0A0qGuIdXhzFPmAvb3ojhjRuyOL4vK6x0yGL8vqywVZcF91aD+E5Ocah2MZRXmdzy6PJuYy8A\nUNVJhjgswRfPdUkY21BvNXuu/DzyrsYd3sVlqksQQkgDM8OEEEIIIeRkoTNMCDkOx4hI7D1H6237\nF/omMyyulDNLrnCKRaSMcKoYsRrGnBlOrnCqFnHmtlUs4o7b4s6wxZkLHupaRtxx28oJXjftAcUp\nDjGJeVc44aOXMarDRocci9joEL5iLGItIzY64L6E9sqNOPer7Azf96uQE46xiXNp/rcwDsAwIhQY\nDm5wsoLDSnRltM6XcsRAyAyLIK9YF96/cXP3ObvHKLXGqAQhpGGvGBaRJwD8ewCPIvwT8ryq/pSI\nvB7ABwE8CeCzAJ5S1S8fb6iEEBKYXW55UmfYHjMT5NJf8nNbq/JpKSPcTpIDkIVwErtncfuOKzGI\nO26DO/F4EsJ33QZAFMNixfA2XF9C/wGKwUQkWmFsRTAAjBBsdJVrD290FcRwFLUbHXDPr/P97vsV\nBpRcshPFfT9gJ+OQF+0I7RiTcArAQ+NYghCV6r1r9XOSaenhShynPabUGiGEHJklzvAWwI+o6sdF\n5KsB/I6I/BqA7wPwEVV9TkSeBfAsgPcdb6iEkIeaC9QVXnQtkwnOrrARadUqck7DYhrVQhr1JLmU\nEU5C+Mw4w1811OL3jttm8XvHbYIYliKG77pNFr9nMsLB40xS3WFfieGhsTFTVnjMK885nOuQRfK5\nDtjoCvf8GkAQw058rjYxxDrHdpEP17rBe/BpgTkVqEh+j+GYlhHH3zrEilwx1uwSsXuJusOEEHIR\n9maGVfUVVf143P5zAJ8C8AYA7wLwQuz2AoDvOtYgCSGEEEIIOQYXsgNE5EkA3wTgYwAeVdVX4qEv\nIsQoeuc8A+AZABgeeeSy4ySEPEQszgtfqK5w07+JUIQ6wqarM+0YkxhM9QjrDLfLK5+ZmMQdN2Y3\nOLS3eI077zjD5wCAu26DsxiVAIA1xuAOpzrFWFhNIg5+hMNdDe4wAKxlwEbL9e/5NZxqcYbRX+HO\nt6XbkvPsBCsINP7QFMir26l6KJBjEiIS3utoX7SWn3e7wtyeleEW1x2enBe7MxtMCFnAYjEsIq8D\n8IsA3quqfyZS/oFSVRXp/7Ojqs8DeB4A7rzxCf7TRMgpcOy/WvfywpPjRey2pdQki+MSk7AT5oY8\ngS5lhn2OSZzFsml1RriI4SSE70QxfFc2uOs2OSZx151nAQwg54dTHGJtyqwBmAhXjyAoU1xiE4Xx\nOordja5wjhFrKeJ3UM2xibm/B1ox7NXlUnGhLfBOzHaaQCdQ1ZgdDpPrREo5Ng0vu0xcFBuVgBHF\ne+oOH5o9YpoQclosEsMiskYQwj+nqr8Ud39JRB5T1VdE5DEArx5rkISQh5xD5IXbusLmuskVLvO4\natE2EccxU7uS4gyvzApyd9w2ryiXJsslJzgJ4SR+X+vuV2J4LVvclU2e0HaGsBpdEr1n4iu9OjTP\nPsZHTL7xOcKiGxukahAj1jrgngbxOzjFoL6saOdRCeK0eEdyez0kr2oHAN6NWKnLzvAorlThaN5b\n2C7lHPJkuiSOBbXgTXWH94nSfeKVK9IRQq7A3sywBAv4pwF8SlV/whz6MICn4/bTAD50+OERQggh\nhBByPJY4w28D8L0Afl9Efi/u+zEAzwF4UUTeA+BzAJ46zhAJISfJRfPC7a6qegSmpdVMNQlBaNtS\nakOTEy754RErZzPD21hLuGSEUzQCQHaFa2e4LM+8hseZeKzjWB1CRd9BZlxOAUZVjLG5Fo+NljiF\ng1Yxi7wdrY+0hHJph6WcUwxihMOoglWuViFYeZ+PD+qNa+wxesmv2TmNUYkS2s1RiTj2qu6weSYA\nXed3NjdMCCEHYq8YVtX/gvk/Ur7jsMMhhJw0F51Y18sK7zy//+f9IHxtubEQkUjtlSuLaqzFx6+Y\n+XWpVFpZRKONRaQvALgrW9yV0YhhxVqAsygYgxgWuDkxDMBDMUZFGHSt5nrEQyOGgSh448u5Kw5w\nwBhXx0jjTXGJtY7YyoBtHJ+XEBHZRnG8Eo9RXH5Pg1Ooxgyxd9V7zbGI2Z8HamGbyqQtzQ1TGBNC\nDgBXoCOEXBtHjXROVp/rLLKRDku9T2JG2LbFZIZDhrh2hq34dUYcn03Eb3KGUzsI4btpRTsAaxEM\nSYyKgzMJttYhHjXUUPZR/G7Uw0lxgMP3+UoUoziMEJzFCXYbGbCRoXqeTcwJh+u7Kj+9FVfea3w3\ntp3ebX7X+T8olSWq6hLHqxO8dzU6QgjBgswwIYQQQgghDyt0hgkhDz4zdYmrLlLXFe4dtzGJyg2W\nUGoNiHlilAiFg8YllouzOkBz6bQhH/elPzQ7EckVXkuKLQyVM+yah0r39Xa/+vzso4brpxjFWjxG\njNho+Of+TMK2He9aRmzi/w5cfD77vCvxOTbhmvc09y7D92DN6pzzu7BuMCGEHBOKYULI7eVQfz3f\nK34Rv7dieFrv1yHsT0smD7Eu8ICyaMZatjnDu5ZtleM9E4/BZISTEE4COGwJBun/4W4AMGqZAedi\n5hgxt3smIVN8Fu/vVaIgDzGNcx2wlq1ZhCOO3zyPE62et/du6n3l/Xmx6rZ9/1f8gTIjTAg5AoxJ\nEEIIIYSQk4XOMCHkoSOXVWsNyh2zqXousD1WXNDGORUPJz47q844rIl2CWTrQqSqEXai3JwrbI97\nHfN5Xst1oFpfvy21Jh7QepzpGVJ/GxFpHfL2/bTHLCIKtT+E+DMxQyWEkBuHYpgQ8uByyb+62xXT\n7HdLTwQupRXC1XWb74dm7vo9QXyh6zbRCMv8+7zgD4gxCELIDUAxTAh5cLlCFHXXxK+rMqorE9pi\njd+87HHsk7+rYhAJJdMQ8rejGjHduMT2WGhrvk73+vH+o3lJox5Gis+/w0v+lkIhTAi5AZgZJoQQ\nQgghJwudYULIyaEq+avdb/ELF4Pw6rLb6qMrnNppkYvkzI4QbLTkeF2sBJGWY95oyAGn6hIpG1zd\nD744wvDYqM+l1EZVbBTV/cJXHI8dZ2z7BU5x71303l/vvRJCyG2GYpgQcm2kv6ov1kqXjEGoxntp\nu79f0FZVKrHn0bTNca8utFHaoxixCYeNrnJu+ExHeDhszApxg4zYmMF6UxrNicKrTCbxWTyK+PWq\n2ECzON4gRCQ2ZgU5bwTvBgM2uqrGO0Ly8fTsud2+G/PsqT0nfif748/k0hPnLngeV58jhCyBYpgQ\n8uCTYqpz4rmzLq9WAhfxu+TvlQCuMr9h/8aX5YzXOua6vRsdQu3euMjFvUbBJUd4jKneEfUiGQME\nDsh1i3uM0Hx0VMWI4CgDwAaCjTojhgfc0zXu6Tq2V9joUI3XqyttP9RiPz57+27q91XeX1cYz/32\no813Qgi5AZgZJoQQQgghJwudYULI7eFAC5V1L62S/z4f3MtdTnHY3ka3dKU+O8Er8VjJAJdiCTrg\nnl/nurz3/BrOebgUe+g4vCNGk+n12ER3GEg1iZHHN0zODRRnOLjBYxz7Bi46w+HM5Aqfx/a5DjiP\nY7bjt07xRofifEenOL2L1hnuOcFpn2onKnFI6CgTQg4AxTAh5PBctLpWJ8bQXkNUoFb9CMqf30Xr\na6jU2VTVsNiDaaqaP/N7h1GAwfl4XDCqwzb+/X8rLtfX3arDfb/K7SB+FfejuIQD4AGfJsCJg3eb\nakLbGiPO4sS4tYwYVLFOi1yoYlgQdh1NhGOjLl9/owPOMVRieKOrLH7v6Rr3fIlN3I9C2Irj+36F\nbYxZbNVh68NXuK8rYhfA6B28L/lite899knt8AuJ1LGJ3raJT0gvd9xyUcFNEU0IMVAME0KunaRd\nF7NPXHcyw3kSXTpuRJhqEMASBWhwO4Owi6dDzCpzWzfARWF8Hys4KM6NYJ3ke6MgBlI1CYe1pIzx\nCndlg3uxw5mMcOKzOHbiMSyYYTaaCXznGPKEt3Md4OGy2D3XoRLDKUOcxPs9v8Z9v8J9H/53cO5X\n2PoB98fQ3nqHrQ753VTC2Dt480uFqsD7IpbzOzc/h+rRlmaGOXGOEHJEmBkmhBBCCCEnC51hQsjt\nY0nMou2jKDs0HlTTVM0OplOt/pzvo6M5ptSFKJx32Kblhcc6udtbkjiVKvMIEYuNC+fclS02OmAt\nwfm96za4h3Vur2WLQXypOwxftXuM8R7hfq5qp2oRaTwpD7xJ1S10hY1f4X5qR2f4PDrD98dVbMeM\nsV/hfBxybGL0Lkc0UkTCVxnh+r1WueGJUyxT13eJq0vnlxByQBaLYREZALwE4Auq+h0i8noAHwTw\nJIDPAnhKVb98jEESQh5ADpUbRrmO5KyqNkI4ZiRsGeEyX67UHE7newcPD8lLHXuMIhCfLtBOW5ti\ny66NZlLZ1gchvI5icyNbrN02i997GoRwilYEMay57TpCeBA/WULZxiRGFIGaaghPJ8RFsasrbPyQ\nYxEpI2zF8Lkfcvt8DJPpNmOaUOewjdujj8/uS/bae4HGGEUWu/aXEqAWs+bnUu0D88KEkOvhIs7w\nDwP4FICvie1nAXxEVZ8TkWdj+30HHh8h5CFlb27YiulelQkraiZraSQlXCZkqaIUxHUCUYHPUV8H\nEWCbL+gxJ4hVBX6oaxBv/IBtdIJXLtQcXvsoft0aaz8aJ3jE2m3NCnS+EsfJdd5dZ7gIY69SiV8f\nJ9NtjNi1dYW33mUBHNr1hLnzMQjhJH636rAZB2xSZngcMPoi/L2dQOeD85sn2HkzaQ4w4riZKJdf\n7oJ9e8Qs88KEkIuyKDMsIo8D+HYAHzC73wXghbj9AoDvOuzQCCGEEEIIOS5LneGfBPCjAL7a7HtU\nVV+J218E8GjvRBF5BsAzADA88sglh0kIeSDZUzd4sjyzLZW243oCgU4iFU3tYC+Ai3WFvYRiE+nX\n/5gfjpXUkidblwQzDVuDeHQeHpIrKmx1xEpGbIfQXnmPc7fCKlaqWLtwPDnDThRrGbMDvJYRg/hq\n+eW0lPMubGzCV5nhUBd4Y+oCb3TIdYI3PuR/Ux3hrbqYCY4Z4XEImeexOMlb70o1idFkhkeXoxHp\nPaeccHpv8AL1xgn2dXWJ3p8HZM45bmnO3esK0zUmhHTYK4ZF5DsAvKqqvyMib+/1UVUV6f8zpKrP\nA3geAO688Qn+U0TIKXLZ/LAV001sQrCj7nDan0SY0yyI81EFvAt7RBWqHmoEproy4W5wIRoBhPq+\nW++wikp6VIdBBpzHzMVKPFZuzMdX4mMUosQg1m7M4teJh6syw/U/k23bNwIwxSTKQhguRzfS8Y26\nXHqt1A0eqnYW0+MQ22WS3HYsYtibjHAQwyUjrColGhHfv1rxm7btBDqgEsfSTqrrlmJjTpgQcjiW\nOMNvA/CdIvIPAdwF8DUi8rMAviQij6nqKyLyGIBXjzlQQsgDzg5B3M0PV4toYFJHOAnicFjL/nRu\nxyWuHEqH3EdFIHB5hpeqwDmFc6UOcRKDg/MYnGIbj228CxPgkvh1oRpEErErF1zflSvOsDM1jEPb\nVwttuAXqzZuXGSbwGbGKepW4rR9y/22s/pDEb3J9rfgNX6ZiRGeSHGCcYPNetXKCpbjB6XW3brAR\nx7NCuNrXc5J3vCgKYULIHvZmhlX1/ar6uKo+CeDdAH5dVb8HwIcBPB27PQ3gQ0cbJSGEEEIIIUfg\nKnWGnwPwooi8B8DnADx1mCERQh5admSIZ93hfBBVTMI6xZPIRLpJvl88IVmIIjEqUWIUNjYBVajX\n3HZOkZJgo3NwzmOIZdlEFIMrzrBIWE45V4VwPrvBQHB9xbbNdjoe9s/nhn2TFw77yndvKjokl9hr\n7fQmJ1ij652cYNVULq1UiFDjDCc3GPH9TWIREydY8mp83SoSKRph98/GJOgKE0IOz4XEsKp+FMBH\n4/b/AfCOww+JEPLQMxOZsMKm0j1thriD2EU2sug153sFRMp+K75zjCKcJw6x7nA6XMSwOIWIg4uT\n80S0EsuCIIAli11U4jdtp1tbUVw9zw6Vp40otDlir0HoV0skW3GsRQADUa9asRu3S3m0sJ3bJtag\nqZTaXAxiSSwCpT27PDMnyhFCjghXoCOE3AxWsCxxim2jJ4yNwE7CuBLFKvV9xNxXAPXIYlmdQsSe\nL0WcShTE6VRRQLTo7LgtRvy2bfvIrtl/GbQjhu0xbSo8TNvI71fj+VUOuFk1LhfayBPi9ohfsy0z\nwnjSF811InSCCSGHZlGdYUIIIYQQQh5G6AwTQm6emSxxt5Sw3bDVJtrriXEhc8R1oVM8SriF6Z/v\nnVzgJopROcf2eLq0fRh7rH3mCzrEbWRCjYNr+1RXte6uNtfQ9rjUDm6v4sNFnOBOn/oBGIkghFwv\nFMOEkNtDZw5c1WzFcdvBTpjrxDBkojzNqYKpOLY78s2l1mtR+Gp1XnutNg9cD+MgawjPieLe8Z4Y\n3TmJbeY9x+2u4N11r7nrguKXEHL9UAwTQm4ve4RPJ2ocz5NdzXhyI1Ah8xecCOP2dp2B7rjWteq5\n3s16w90leOeuM3v9/sMv0vwUu4SQa4aZYUIIIYQQcrLQGSaEPLgsdBH7juSsr3xh5IDXIoQQcr1Q\nDBNCDs81/6n70n9+X/rn/AWRgcmVbroE2C59vicR0T2/c70l8ZO9Y9l3vWPC32EIIaAYJoRclmsU\nvJeaVGWU1ay4tQUe2uP94g/zx+L9pHPtyb0W5ngXs0+o7hC2VdUMe17zvXe9NBGwGrrsylan/f2H\n3adNDy6W5945RTIhJwUzw4QQQggh5GShM0wIWcZtcIJn9zeubNu/Lglc9qX9nZXQxNTPrY7rjv6d\nMbZ9J+fOnbOHiUs6026d3t55KtOay9V5TZ3l9ppir21L1LUV7tIxWwN58kOZf67ZutOHZk+JP0LI\nwwXFMCFkntsggIHlCzO0ArU9vxW/VVuq64ivz5N2e+e1mjF1xPKSOMUiuuK205Z+uxWwOhG/YTun\nH0Sgrr2W+W2juVZbU1m1HrM9PKkbPVMv2h7OYz4mO8ZACHnwoRgmhPS5BiF8IQEMzOeAew7wDkc2\nCN094tebc725X9zOx5tzRVMfLcd3iOd2bN3na+i7u3V7t7iVHJLLx63AdaZ/3FYbqjN9RQB10UFv\nxbFr7q21OJ4zYKuVBu2L0KZjatou1yGMKYgJeahgZpgQQgghhJwsdIYJITXXFI24cixiLhKRzu9F\nFXydAd7l7ooH4Js2wj7x7bma28UZrq9dxywUE6fYPlvvHVQvwDxqLxOc3F8gO7N1tEErZ1hdcXez\nM5z6u45TrPW9KsfWAbbMhDo1EYtOhlinj1o5xPZnn3MR0/dgu1xbbIIOMSEPBRTDhJDCTQvhHSL4\nIrGI2UxwI36r9g7BW/X1jfjt9E0COR+vhLd2YxM5VtF7Dz0qQVwmwBUxnG4glZjNwrdqB9Fat2MH\nh0ocq0MW0pL2W3GopW8QwSXqkDPEzVizdtZyDfuIO2MTM4K4Ou9YMDJByEMBxTAhJHDbMsI9ITzn\nBs85wTkjLLud4J74HacCGADciCyWU1839sRxEcMTceyNHZozxeXhLlxNQmRa/SGKWXU6I37rth+k\naiOKYz8AcAKx/Ye4nUS9OWaz1mXSXhLWOpmL1jrF1QS7uF2J24WCGKBLTAhZBjPDhBBCCCHkZFnk\nDIvI1wL4AIBvQPg9+J8A+DSADwJ4EsBnATylql8+yigJIcfjpqMRvTH07LyeK9xmbRu3NWWEsxPs\n675ujP1bZ3hM7rBxc2Pf4AQbpzf1Nc6xeDX9tXKWReN1c0xC65iER+USz5Ld1jYmEbY1Ob3JJc7u\nrUAHzc6xH8K+9Kwa2z66v84j9o/3HUzwd8AkjuKH8h41DivHJrz03eGm2kR+/GQdtx+HdMJtyRCn\nMdAdJuSBZGlM4qcA/CdV/UcicgbgNQB+DMBHVPU5EXkWwLMA3nekcRJCTpDZRTRa2gxxzAjXpdTq\nGEUljpMYTgI2ClpXtcu2G2ux60YjnON2EctaxybSMVt6bfRNJlr720ARwaZGmQqgQ1KciDnhKIaH\nsJ1jEanKYjpDAAAXo0lEQVQdxa4MQQAn8as+iOWkIHUIQ0j9PebVpcZhqRmKmnhKyBcLJP5gNHdC\nOWFyQfPo1yVsCSEnxV4xLCJ/BcDfBfB9AKCq5wDOReRdAN4eu70A4KOgGCbkweKmXeElOeFOv9kJ\nc+Y6E6d4SUZ4InBjOwrgsl3ErxsVsq3FrtsaMZy283FfieHgEhtxbGsU70CNMyxSZqHpII0YdoAz\n7VXY9qvY9gIdAYn/N/ADoFrsWdWSJ463C4I4jH46MDuxL2aIs6scfyHJi2vERTVsRBiALUbRHGgE\nMSfUEUIOwJLM8JsA/DGAfysivysiHxCR1wJ4VFVfiX2+CODR3ski8oyIvCQiL41/8ZXDjJoQQggh\nhJADsCQmsQLwtwD8kKp+TER+CiESkVFVFel7P6r6PIDnAeDOG5+4Jh+KEPLQoTtMt+gC5+oRTbt3\n3Nb+bStAJFc4OcM2FpFcYbeNfbfhmGx96esVsjVO8KiA97EdtpOzDO+Da5xjE2YbZUz5UfNqbyYu\nIQKk/c5BB4G44HXo4IG4DwDUO+hKcp7aDwJduXLLEHGGi287bCs8anc4HIshX1vJwqOyejW9//je\n1f4cFbXp3/yAU8yCbish5JgsEcMvA3hZVT8W27+AIIa/JCKPqeorIvIYgFePNUhCyBG4hl9ND1pK\nrT1XO9e3fwdvYxT2eikyUWWGS3m03oS6kh8OQthlsRuEcG5vfewTT/Y+COU8Gc8Dozfi2McsQupf\ni+GW/IQ5JiGAcyYm4YIQThli54DBAyuXx++9CxPhADh18FHyAkH4plcU2hL0rRG8anR4+M2iTNbL\ngjgN1pZGSz+TmZiD7BHHLLVGCDkGe8Wwqn5RRD4vIn9dVT8N4B0APhm/ngbwXPz+oaOOlBByuuwT\n7q0wTm2t23PO8LQuce0U166x1vniKISTM5yFcG6PUQxHhTiOQQAn8duIYbWieA5nfnFwrhLDMrjo\nBEfxOwyAL8vG6ZByv8XfdUiCGABcNJqjKJbgQmeHWqQIY2cEbnyvrTNcub8pM9xmgFHaXSgwCSFH\nZGk1iR8C8HOxksQfAPh+hH9PXxSR9wD4HICnjjNEQgghhBBCjsMiMayqvwfgrZ1D7zjscAghpE/X\nFGzM053xisYpbtttHeJ6yWWtIhO2dFpwirVUl4iusGzG2H8ERg/ZRivZe2A7AmNoa8oM+3RcAbXL\nuM0Ql4QLNYSHkhkeHTAMkFwLTaE6VHWMFah8YYhms9m7mI9OTrMDRLTUbBaFpMoUHlUptd577dcJ\nLu+5ii10+trKa4QQcgy4HDMhp8Y1KYuDLrLR7tuXF57tU4+tfNcqRtGfcJciE1oWzkjt0YfoAxCE\n8DiGbDBQhHAUv0kIa+rv47F4fU2RiV3l1STmhBEn0rkRGIL4VReFbzp/NdSlfEWAsSljJijXG6PY\ntRERJ9Pcb3p/iiK0NVQQrvraYadhVQL4crnheh8X4SCEXB6KYULIrWGngF7CLv3Yc5FbJ9mK3+QW\np/2Vi1yc4TQZr4jlWCkiit0shE1mWEcPbLelrVrEsvdREO/ODKdqEuqiE5z6DwrFajLRrsxZi+em\nh/MSBGrlgpcV6vIvBum464jd5pcUO/KJuLVcUTByEQ5CyCFYUmeYEEIIIYSQhxI6w4SQB56u6wvM\nxiK66Mw2Oi6xL3ZniE2YWELaTjGJ5PKmjPDoSzQCgI4joBq+AzkzrDucYVGBSlpfOaznptYBdmJi\nEGOMVSQnOy4JZ8YrxolWjRGJaLmmqhvVaHa8q/a9VbEJsx/opBvo8hJCbgCKYULIw8kFMstdkazY\nndu1eC21dT2KIIbZzm1fohHxuNrM8DjWE+gmwduwfHIatGCAYgwl1BAyxzKOKGsau/547Hi9AsOe\nZ9R5xdqK3b3vnqKXEHKLoBgmhJBDY13dKD6t+A19fP29rSAxJ8RVi9DN5w3lOiJhwY1UV1gb5zqN\nb5/4JYSQE4GZYUIIIYQQcrLQGSaEPJzs+lN8W4VNOt0FxoG9YJkLF8qXheu0hXPTMsrRi0hZ4WpA\ne+6nOv986bpz903juyjp/F4VvHbfBd49IYTcNBTDhJCHjly61lQHE5gJW72TZGY7nqdWY7pSG1cl\n1OHNf2dzCELYik8xpdCAuGBGjDU4FwRxmhCHmPfdV2c4b7tynXYfYl3hZix5nOm7kyJo8/PFZ3d7\nxO4OcZvOq66NbuyYEEJuDIphQsgDTyt+e0audX81LjQhjUizYi195d2mTzgWqy2IlNq9YUcQpmLE\nrvhys2EIFRySOI7nZrGtAvXODLbJEkehm1aByzWHzXjgJC/CgXTMOsbOGac3jB9mPOn9tO/Evpve\ne0Pz3ux5XSiOCSG3AGaGCSGEEELIyUJnmBByOnTiD2G/QKXkcPMSxTYq4CWfE45H+znFDLJTG5Zo\nkyF4Dep92Nbg1Ioq1A3AkEqjpVq8MXYBQGytt9Y2zbGI6AQPQ4hhJCd4GAA3FKd6GIDBlfENLri/\n0SlWF53kZI0IiluM6BQ71C7wxEEv76YeKwgh5NZDMUwIuTW0izBMsH9+b/vZtXmjYKuu1+Zce8Kt\nEnkCdWmJZZMRHkJyQYcoANOxKH7FK+BNyHgYYimzNNgVBFto+ufXCcRJWIwDMeZgF+3o/QHPxhyc\nC2I7iWE3hPYqXn9woX8+7oBBwhcAHVzMRKd2EL+aSq8lcZzW+JAmX7zjPeY2ynudHq9/2L2IxhyM\nVxBCDgHFMCHkdpIs00h3rYYmIzwR01ZU2Ql12d2Nx2NZ3rzKnCQ3OPZ3RdtCpWqrE2DlILFzauch\nqkJ1aMa+CivDIYlpgYipN+xNXeLeq0m5YCA4xUYMSxK+QxHjuiptXbnwZTLHtq3R5S7tehJd6xKr\nQ+Mq784UW5Y4yZO3QAFMCDkwzAwTQgghhJCThc4wIeTmsRGHS52P2kIUU2EiYv/8Lsm9zG4vYhTA\nrCusJjqgWpxU1RyVAAC/EjgFdGWtU5TljldaD89JcIVT9QcfKk1oilGMY6g2kfr7xht15iEAYIj5\n4CoG4YIbnNqrAbpKGWEXohGpvQousF8ZJ3ione/gDtdOcd42sYnsIPcqT6Ds60YrrsJsroYQQvZD\nMUwIOQqtGK1oxes+FvQPpdPUNupYhM1ZtNnVNmOc2kn0qY1JxFytqXimCviorF1TCQ0OsHWGsR1D\nOTVnlmN2LohiIIhXXRCTsIt3RAGc2ppywkAWwln8RiGssb8fXBDEOUMsOTcMmJhEE4VI77X3HufK\nsvXywot+B7qoWN7RnzljQkgLYxKEEEIIIeRkWeQMi8g/B/BPEbyZ3wfw/QBeA+CDAJ4E8FkAT6nq\nl48ySkLISbHTVQaK42ibTSSiKi6RnF57fjOBLkQjSiwiRSVCW+BjjMGpQFVzrMAp4AeTuKjyF3bM\nZTDiPdQ4w+I88opzPkygqwxM1XJ+oirlJtnpzTGJ7Ay76PaaCXTRDQ5tgR/KhDm/qmMRfihRCaCN\nSTSTCaWORrTO8WRFuvxu+v3noLtLCDkke8WwiLwBwD8D8BZV/UsReRHAuwG8BcBHVPU5EXkWwLMA\n3nfU0RJCHl56ueHU1Orb7HLK+Xhef9kcswI76VV73JXjmoRw1U4xAoWHwMWDQcMKSjrCQUTL6sip\n8kPKI4uE0msulWLz0JUCoy83Uw0l2mZQJ5W4xlBKrWksnTYvhiVHIwAEITwINP7fQIe0r9xrkhM2\nSzlPahDbNjDNCO861j5nu6PXl3lhQsgVWZoZXgH4KhHZIDjC/xvA+wG8PR5/AcBHQTFMCFnKIXLD\nxiFuaw4rJOvFOD+udpQbA1eqXLAUpxZhLpzL20EIe6PMbKUzn3RqWjZZFOJKRlgGAUZX6gin7aEs\nyrEvMwyTGVaRuPBHygRLzA0ncRyEsM0E+yYjHNpx/FEI+3w8usNZHNeT6apstWsmzTVZ7Lzdqy2c\n+k+edf419N/NBfsTQk6evZlhVf0CgH8F4I8AvALg/6nqfwbwqKq+Ert9EcCjvfNF5BkReUlEXhr/\n4isHGjYhhBBCCCFXZ0lM4hEA7wLwJgB/CuA/isj32D6qqiL9v1Wp6vMAngeAO298gn/PIuSmaaIH\nx2TvinJ7zq1iDe0CHK1TLM3xzqpzkxXpbFbVWgPZOtbc3xq1NibROsUunupT6TTxECc5NqGjQpyG\nVe0QYhfiXb6B+uIOz2KWS07VKqpFNOzy0EOMRaT7r1wsnWYywhMnWConGDYm0as00TjB+TW2GeA2\nM9yLSUjno9le85Jca9aYDjUhDwxLYhJ/D8AfquofA4CI/BKAvwPgSyLymKq+IiKPAXj1iOMkhJwC\nObQ7oyTmhHwUUPmsqIbTNLQUR86xCYcq+yBAFbWo4hIAPHTyZ7QifhVeBNkPEIF4E+l1DjJqzvDK\nAIhXyFgCyeI1P5OoQk1euP1FIr+aFFUQiVGFkiFWMyEu1A2WOjaxaibIGYGbhHARv2USXb6eKbNm\nV6DLMQl7fCKOzQPl2AT67BOUzAsTQg7AEjH8RwC+RUReA+AvAbwDwEsAvgLgaQDPxe8fOtYgCSEP\nKXtc6omT23OHe+dncWwmrdn7mT5pW4yIm8aTBT7ucXFcOUMsAhm11BEW1GLYI2SGYzUK8a0Yjm3j\nDIeNPZnhRFr4wpX7T8SwaybEDY0YjoI59fcD8gPmahKNU5y3jeDdX3dY94rfK7nCdGMJIZdgrxhW\n1Y+JyC8A+DiALYDfRYg9vA7AiyLyHgCfA/DUMQdKCCGEEELIoVlUTUJVfxzAjze77yO4xISQB5Gp\n/Xk09tYNbtm3PPNMXrgtvVY5vzDOrw/bkqqZ9coCN7dIbqyH5upoABBWVha46PR6AcRLKiYRjpt7\nBWcYwFAiHFAt78fH7T11lm1cIritKTaBmPFtnOEm82tzv97mgJNTPNMfzbnd2ISNo+xwjmfrCh84\nHsGsMCFkF1yOmRBy8+wR5pPllPfUHU6rMdvCwwoppydxbNITc4LYXibvFa2jAKOZMJdiEEn8OikC\nOB6HSq4jnNspFqG9rEAHc/+UGwaK+G3bUwFbZ4Dr0mmdjHB30Y1G/Aoq8dstpdbkiC1zEQnGIwgh\nx4RimBByLRzTHa7ywyj7qzrEMBPqXMjmJvdWfUcQG9FmaxQjurDlWHCKa+dXTFuz4M3HtfSHJnc4\nOcWN0p979nR6mkBnBagr45vWAu4vomGP73R/TfWI1vnNx8yKdDkn3Iy7vMvWgd/x3MDtdoUJIQ8k\ne+sME0IIIYQQ8rBCZ5iQU+Yac8N76Y3FlFrr1h1ucwzWMa2bVXUJhWR3GIhOsXFr1YUTqss3f/rP\nTrSEgdk8so1FhEiEcY41lFFLdYZTPrgYnoKdlSRyt/J0tkJDLm2W203mV6ZRB1uarRuLkE7bXKut\n11yON9UjBJN2/UzT/ZWz23OFb5Pze5vGQghZDMUwIeTa2LsIRz+k2z/fiueoH7MWMfGI3mVTZCLH\nJKIQribYAeVvZ4JQlxhGCJuogI1FhLbJDGuMYaSbJ3Gc2yYakdv7VVUbPVBT2q1b69dEF2wmuBXH\nc+3JpDdMIxhtRnjJhLnqF41m+6pZYUYkCCFLoBgm5NTZI0CPwVXzw90JdcZQFWAiuHoucc4Qo7jB\nVSedXkB8426mexshDSfFdc5OsHkUc22bF85tYPfPoycYZwRnzynu5n53VIRonedWDNfHmoywYNLO\nY2+eZXbC3IOQE6bwJuSBhZlhQgghhBBystAZJoQEepndI7LTHe6NZc4dtv2Ts5giEzMuZHFnTeUG\nlRJRtn3SjhSDGKautK1ZjDZ/HK9rr5Xc4jQW+6yHcoZtbGJabaLub9s5JmHv0R5P29k13hGLQH1u\n9WjtcWCZK7zDhb12V5iOMCEPPBTDhJDbyZwgBrLi2Zchzqc1gmVvbMJNowxZZMX91b0VJZPb5o87\nsQhtYxFtPKN65HqHtg/TE8NtuxNzmOZ8m+MLog6LYhF2LO2Ym+2rCmFCCLkMFMOEkMINuMP51r37\ntorVdl6QIZ6tQ5wOLXWKJ2LWjCeJYzs0h8m51nWu2jNOcOk/VX9T0VhvV5lmNO2e+EXT7mWEGwE8\ne8yOpxXBnbH2n6d9GdjJjU2UozAn5KGAmWFCCCGEEHKy0BkmhNRYd/Ua2VthonWtbWeV2mXuXb+N\nTTTOMaRxa02uNw8wde44w5ULrc397KV02q6+dx5v8ix7nOGqT899lWl7iRNcHd9zHeB6YhGsHEEI\nuSoUw4SQPtccmQCaHG6POaG+LzbR3idFHqxgbHLBE8HaLpPciuOqb92eCOM977W69wxLa/DORhd2\n7e+J356Ibb9jJhLRRlSWiODOeS0UwoSQQ0AxTAiZZy6ze2QuJYqbzm3WtusW2w0xfaQRrJ3Mbzpb\ntbm1SnaHpzcyQrx9Nu1uLmOX2OwJ1uq41vvm+s+5uzvGsTgTvOsaDRTAhJBDw8wwIYQQQgg5WegM\nE0KWcYOxCWBBtQlgcabY5oPzddpYhL237Mn4VmPT6bU6wyrHp7bjkgXXZh3S1iGfyxcvcZPNdndI\nS2MQM2PrndvjRpdVpitMyEMPxTAhZDk3FJsApoLoQuK403lSsswK5FbcdYTzTpHUG1srgvPQemO7\nPF3huGDf5J67Juntvd/lYxCz17xObvr+hJBrhWKYEHI5doZwj88+wbRrAl3pkC6255o7brZ3stue\nd6IL+ixiqYC7bCZ3iV29cAw3LnYtt2kshJAbgZlhQgghhBBystAZJoQcjou4bEd2kS/kPu6LYOxw\nRWePXHOM5EIcyA29VQ6v5baOixByKxHdV/DykDcT+WMAXwHwJ9d2U/Kg8HXg54LU8DNBevBzQXrw\nc0F6/DVV/fp9na5VDAOAiLykqm+91puSWw8/F6SFnwnSg58L0oOfC3IVmBkmhBBCCCEnC8UwIYQQ\nQgg5WW5CDD9/A/cktx9+LkgLPxOkBz8XpAc/F+TSXHtmmBBCCCGEkNsCYxKEEEIIIeRkoRgmhBBC\nCCEny7WJYRF5p4h8WkQ+IyLPXtd9yc0jIj8jIq+KyH83+14vIr8mIv8rfn/EHHt//Jx8WkT+/s2M\nmhwbEXlCRH5DRD4pIp8QkR+O+/nZOGFE5K6I/JaI/Lf4ufiXcT8/FyeOiAwi8rsi8suxzc8EOQjX\nIoZFZADwrwH8AwBvAfDdIvKW67g3uRX8OwDvbPY9C+AjqvpmAB+JbcTPxbsB/I14zr+Jnx/y8LEF\n8COq+hYA3wLgB+LPn5+N0+Y+gG9V1b8J4BsBvFNEvgX8XBDghwF8yrT5mSAH4bqc4W8G8BlV/QNV\nPQfw8wDedU33JjeMqv4mgP/b7H4XgBfi9gsAvsvs/3lVva+qfwjgMwifH/KQoaqvqOrH4/afI/xP\n7g3gZ+Ok0cBfxOY6fin4uThpRORxAN8O4ANmNz8T5CBclxh+A4DPm/bLcR85XR5V1Vfi9hcBPBq3\n+Vk5QUTkSQDfBOBj4Gfj5Il/Dv89AK8C+DVV5eeC/CSAHwXgzT5+JshB4AQ6cuNoqO/HGn8nioi8\nDsAvAnivqv6ZPcbPxmmiqqOqfiOAxwF8s4h8Q3Ocn4sTQkS+A8Crqvo7c334mSBX4brE8BcAPGHa\nj8d95HT5kog8BgDx+6txPz8rJ4SIrBGE8M+p6i/F3fxsEACAqv4pgN9AyH3yc3G6vA3Ad4rIZxFi\nlt8qIj8LfibIgbguMfzbAN4sIm8SkTOEYPuHr+ne5HbyYQBPx+2nAXzI7H+3iNwRkTcBeDOA37qB\n8ZEjIyIC4KcBfEpVf8Ic4mfjhBGRrxeRr43bXwXg2wD8D/BzcbKo6vtV9XFVfRJBP/y6qn4P+Jkg\nB2J1HTdR1a2I/CCAXwUwAPgZVf3Eddyb3Dwi8h8AvB3A14nIywB+HMBzAF4UkfcA+ByApwBAVT8h\nIi8C+CRCtYEfUNXxRgZOjs3bAHwvgN+P+VAA+DHws3HqPAbghTj73wF4UVV/WUT+K/i5IDX8t4Ic\nBC7HTAghhBBCThZOoCOEEEIIIScLxTAhhBBCCDlZKIYJIYQQQsjJQjFMCCGEEEJOFophQgghhBBy\nslAME0IIIYSQk4VimBBCCCGEnCz/H6VNP7bVtmE1AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f20185c9588>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "data = np.genfromtxt(\"output/v_00003\")\n",
    "fig = plt.figure(figsize=(12,12))\n",
    "pltres = plt.imshow(data[::-1,:])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}