ISWI2016_Py1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline #this is a formatting statement for jupyter notebooks "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "# Using Python for your science\n",
    "\n",
    "\n",
    "This is an introduction to moving your science from an IDL environment to Python. If you havne't done so yet, please download and install the the latest versions of anaconda (www.continuum.io/downloads) and sunpy (docs.sunpy.org/en/stable/guide/installation/).\n",
    "\n",
    "Let's make sure that your instillation of anaconda is up to date. From the terminal, type:\n"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "conda update conda"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "  "
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "conda update anaconda"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, lets launch an ipython (interactive python) session"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "ipython"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Getting started with Python"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The first thing to note is that by itself, Python is just a platform. You can think of Python as an engine - it is very powerful but unless you connect it to a transmission, axle, and wheels, it isn't very effective at getting you down the road. \n",
    "\n",
    "Fortunatly in most cases you don't need to reinvent the wheel and you can just find some that someone else created that will fit your needs. \n",
    "\n",
    "The Anaconda package comes with literally thousands of tools to help you work with data. \n",
    "\n",
    "Let's explore some basic tools. NumPy is the fundamental package for scientific computing with Python. Let's load it in to our current session. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Notice the syntax we are using: import \"the library\" as \"what we want to call it.\" \n",
    "Now whenever we want to call a tool in numpy, we can use \"np.\"\n",
    "\n",
    "Next we will want to plot our data, so lets grab some plotting tools. We will use a package called matplotlib. Matplotlib is a 2D plotting library, which produces all kinds figures. We don't need everything that matplotlib can do, so let's import just the tool that produces line plots. This tool is called pyplot. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Again notice the syntax we are using: import \"library.tool\" as \"what we want to call it.\"\n",
    "\n",
    "Okay, we have the tools 'np' and 'plt' defined. Loading in the tools and libraries you want to use is something you have to do every time you write a program or load a python session. \n",
    "\n",
    "Now lets put our tools to use:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "dt = 0.1\n",
    "Fs = 1/dt\n",
    "t=np.arange(0,Fs,dt)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can dynamically define variables (i.e. we don't have to specify the type) and use conventional assignments to make new variables. \n",
    "\n",
    "We also use a tool within np to create t: np.arange. Np.arange is similar to IDL's FINDGEN( ). It creates an evenly spaced array with the syntax of output=np.arange(start,stop,step). \n",
    "\n",
    "So what does t look like? We just type the variable or use a print statement:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 0. ,  0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.9,  1. ,\n",
       "        1.1,  1.2,  1.3,  1.4,  1.5,  1.6,  1.7,  1.8,  1.9,  2. ,  2.1,\n",
       "        2.2,  2.3,  2.4,  2.5,  2.6,  2.7,  2.8,  2.9,  3. ,  3.1,  3.2,\n",
       "        3.3,  3.4,  3.5,  3.6,  3.7,  3.8,  3.9,  4. ,  4.1,  4.2,  4.3,\n",
       "        4.4,  4.5,  4.6,  4.7,  4.8,  4.9,  5. ,  5.1,  5.2,  5.3,  5.4,\n",
       "        5.5,  5.6,  5.7,  5.8,  5.9,  6. ,  6.1,  6.2,  6.3,  6.4,  6.5,\n",
       "        6.6,  6.7,  6.8,  6.9,  7. ,  7.1,  7.2,  7.3,  7.4,  7.5,  7.6,\n",
       "        7.7,  7.8,  7.9,  8. ,  8.1,  8.2,  8.3,  8.4,  8.5,  8.6,  8.7,\n",
       "        8.8,  8.9,  9. ,  9.1,  9.2,  9.3,  9.4,  9.5,  9.6,  9.7,  9.8,\n",
       "        9.9])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 0.   0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.   1.1  1.2  1.3  1.4\n",
      "  1.5  1.6  1.7  1.8  1.9  2.   2.1  2.2  2.3  2.4  2.5  2.6  2.7  2.8  2.9\n",
      "  3.   3.1  3.2  3.3  3.4  3.5  3.6  3.7  3.8  3.9  4.   4.1  4.2  4.3  4.4\n",
      "  4.5  4.6  4.7  4.8  4.9  5.   5.1  5.2  5.3  5.4  5.5  5.6  5.7  5.8  5.9\n",
      "  6.   6.1  6.2  6.3  6.4  6.5  6.6  6.7  6.8  6.9  7.   7.1  7.2  7.3  7.4\n",
      "  7.5  7.6  7.7  7.8  7.9  8.   8.1  8.2  8.3  8.4  8.5  8.6  8.7  8.8  8.9\n",
      "  9.   9.1  9.2  9.3  9.4  9.5  9.6  9.7  9.8  9.9]\n"
     ]
    }
   ],
   "source": [
    "print(t)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "But what are the dimensions of t? How many elements are in it? \n",
    "\n",
    "In IDL we would have to call a function to return this information, such as MAX( ). In Python, since it is an object oriented programming language, the variable t is actually more than just a simple array. It is self-aware! So let's ask t about itself. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<function min>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.min"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Um... what just happened here? Python just said \"yep, that is a valid function of t\" but didn't actually return the value. To get the value out, we need to use a different syntax:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.min()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Also we can ask:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "9.9000000000000004"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.max()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4.9500000000000002"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "8.3325000000000014"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.var()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are all \"t.something( )\" - \"something( )\" is a function of the object \"t\".\n",
    "To find out how many elements are the array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "100"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.size"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To see the shape of the array:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(100,)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "t.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you want to see all of the functions already available in t, type \"t.\" and then 'tab'. This will show you dynamically all of the functions that exist. Go ahead, try it!\n",
    "\n",
    "You'll note that .size and .shape don't have brackets after them, but .min and .var do.  This is because the size and the shape of the array are attributes of the array.  The reason .min (and other methods) have brackets is that these functions have different options that you may want to use.\n",
    "\n",
    "Let's create a 2d random array\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "rand_nn = np.random.randn(t.size, t.size)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can find the minimum values along one dimension of the array..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-2.89518052, -2.38851405, -2.65438088, -2.42616205, -2.98637711,\n",
       "       -2.61917924, -3.04400683, -1.99551291, -2.54481928, -2.75898872,\n",
       "       -2.00846906, -2.82244842, -1.55895129, -3.62905027, -2.15917153,\n",
       "       -2.27459189, -2.18018518, -3.02437466, -2.42973237, -2.77289929,\n",
       "       -2.61672753, -2.00925514, -2.61912512, -2.58307325, -2.45159215,\n",
       "       -2.53986963, -2.56444877, -1.82938401, -2.29653528, -2.55571806,\n",
       "       -2.94060788, -2.46487206, -2.05620166, -2.37845073, -2.43397154,\n",
       "       -2.14122798, -3.25831335, -2.2645945 , -2.41222808, -2.84049392,\n",
       "       -4.08426662, -2.80068056, -3.2978516 , -2.49876356, -2.60339569,\n",
       "       -2.02333597, -2.80719912, -2.09827113, -2.52161858, -2.48798412,\n",
       "       -2.11527522, -2.04119045, -2.11583545, -2.58030957, -2.03842926,\n",
       "       -2.63240217, -2.38451646, -2.13098785, -2.81934984, -2.13577515,\n",
       "       -2.48779106, -2.38749245, -2.46409163, -3.27733676, -2.60068506,\n",
       "       -2.67887403, -2.04875584, -2.54092582, -2.37890464, -1.89240065,\n",
       "       -2.19670845, -2.55542198, -2.5086707 , -2.12344642, -2.06953826,\n",
       "       -2.00131686, -2.17518508, -2.82904241, -2.58758977, -2.43779972,\n",
       "       -2.51828533, -2.18251377, -2.35536258, -2.49192285, -2.78535396,\n",
       "       -2.20373031, -2.50341687, -2.18704082, -2.2692945 , -1.94611425,\n",
       "       -2.24584646, -2.44274415, -2.44158504, -2.35887649, -2.79038533,\n",
       "       -2.11925456, -3.33111796, -1.94091676, -2.15708445, -2.88000045])"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rand_nn.min(axis=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "or the other one..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-2.06035846, -2.42968466, -2.14733083, -2.41222808, -2.94060788,\n",
       "       -2.09774047, -3.27733676, -2.46409163, -2.81169946, -2.53165595,\n",
       "       -3.33111796, -2.00255171, -2.53416261, -2.00131686, -1.86039766,\n",
       "       -1.79499492, -2.6139758 , -2.80068056, -2.77289929, -2.45927262,\n",
       "       -2.43397154, -2.02778759, -2.56444877, -2.4963192 , -2.42973237,\n",
       "       -2.05636977, -2.61672753, -1.80562577, -2.80719912, -2.82904241,\n",
       "       -2.50341687, -1.87812105, -2.74139925, -2.44158504, -3.02437466,\n",
       "       -3.62905027, -2.63240217, -2.43325322, -2.42289984, -2.61917924,\n",
       "       -2.61912512, -2.79038533, -2.2645945 , -2.54092582, -2.88000045,\n",
       "       -3.2978516 , -1.88906553, -2.78535396, -2.38709958, -2.35887649,\n",
       "       -2.58758977, -2.29667604, -2.38851405, -2.96679138, -2.38125347,\n",
       "       -2.60386995, -2.45159215, -2.55542198, -2.68562536, -1.90493965,\n",
       "       -2.11583545, -2.18248954, -2.60339569, -2.29653528, -2.74144036,\n",
       "       -2.82244842, -2.13098785, -2.38451646, -3.04400683, -1.93908917,\n",
       "       -2.89518052, -2.67887403, -2.77061606, -2.15917153, -2.43857984,\n",
       "       -2.05672415, -2.31948536, -2.32891487, -2.43779972, -2.54481928,\n",
       "       -2.98637711, -2.37890464, -2.34942342, -2.58030957, -2.33437098,\n",
       "       -2.30671029, -2.55571806, -2.60193192, -2.24663844, -3.25831335,\n",
       "       -2.65438088, -2.18251377, -2.40651613, -2.50890559, -2.30407762,\n",
       "       -2.29585467, -2.1856842 , -2.84049392, -4.08426662, -2.41878524])"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rand_nn.min(axis=1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Or the whole array."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-4.0842666160402041, -4.0842666160402041)"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rand_nn.min(), rand_nn.min(axis=None)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The min( ) function with the axis option finds the minimum value along different dimensions of the array rand_nn."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We used \"axis=None\".  There is a special value in Python called **None**, which indicates **no value**.  In this case, this means that axis has no assigned value.  The axis function then resorts to a default behavior, which is to give you the minimum value over the entire array."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "Now let's create an array of random numbers with the same number of elements as t. We will use the Numpy tool 'random', but there are lots of ways to generate random numbers. In IDL we might use RANDOMN( ) or RANDOMU( ) depending on the distribution we want to sample from. Numpy.random had many ways of calculating random numbers too. \n",
    "\n",
    "Try typing \"np.random.\" and then 'tab'. \n",
    "\n",
    "Let's sample from a normal distribution:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "rand_n = np.random.randn(t.size)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Notice that we can use the t.size attribute of 't' without having to define another variable or call a separate function.\n",
    "\n",
    "Let's plot our variables. To do this we will use the matplotlib's plot function. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEACAYAAACnJV25AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmUHUd59p/S7KtGtixZshZLlmQZ432XMR5sSByTYAi7\nSSAEThLORxxIDuHD5ICUhBMwCQnHhARjICbHmMWGBPxh43W8gAGDjbHkRd40I1nLSLMvmpFmpr4/\n3infmr7V1dXd1X2393fOHM1c3enb08vTTz31VpWQUoJhGIapHhaVegcYhmEYv7CwMwzDVBks7AzD\nMFUGCzvDMEyVwcLOMAxTZbCwMwzDVBn1aTcghGgC8CCAxvnt3Sql3JZ2uwzDMEwyhI86diFEq5Ry\nUghRB+CnAK6RUv4y9YYZhmGY2HiJYqSUk/PfNoFcO496YhiGKRFehF0IsUgI8TiA/QDullI+6mO7\nDMMwTHx8OfY5KeVZAFYBuEAI8Sof22UYhmHik7rzVEdKOSqEuB/AFQCe0v9PCMHxDMMwTAKklCLO\n+1M7diHEUiHE4vnvWwC8AcAzITvHX1Li05/+dMn3oVy++FjwseBjYf9Kgg/HvgLATUKIRaAHxXek\nlD/2sF2GYRgmAamFXUr5JICzPewLwzAM4wEeeVoCuru7S70LZQMfiwJ8LArwsUiHlwFKTh8khMzr\nsxiGYaoFIQRk3p2nDMMwTHnBws4wDFNlsLAzDMNUGSzsDMMwVQYLO8MwTJXBws4wDFNlsLAzTJkx\nNgY8/XSp94KpZFjYGabMuO024BOfKPVeMJUMCzvDlBl9fcDkZPT7GCYMFnaGKTN6e4HDh0u9F0wl\nw8LOMGUGCzuTFq8LbTAMk56+PqChodR7wVQy7NgZpoyYm+OMnUkPCzvDlBEHDwJHj3IUw6SDhZ1h\nyojeXmDdOhZ2Jh0s7AxTRvT1AZs3s7Az6WBhj8HzzwPf+U6p94KpZnp7gQ0bgNlZYGam1HvDVCos\n7DH45S+BW27JbvtTU9ltm6kM+vqAtWuBlhZ27UxyWNhjMD2drfiedhrwm99kt32m/OntBdasYWFn\n0sHCHoPpafrKioEB4IYbsts+U/4ox97aysLOJIeFPQZTU9k69ulp4FvfAiYmsvsMprxhx874gIU9\nBllHMVNTwEUXAd/9bnafwZQv4+M0MOm440jYeZASkxQW9hhkKeyqAuJDH+I4plbp6yO3LgQ7diYd\nLOwxmJrKLmOfngaamoArr6Qb/Mkns/kcpnzp7aV8HeCMnUkHC3sMsnTs09NAczNQXw984APAV7+a\nzecw5Yty7AA7diYdLOwxyLLzdGqKHDtAwn7zzXxj1xq6Y+eMnUlDamEXQqwSQtwnhNghhHhSCHGN\njx0rR7J27ErY164Fzj8f+MEPsvkspjxhx874wsd87DMA/lpK+RshRDuAXwsh7pJSPuNh22WFqmOX\nkjq4fG9bCTsAnHUWsGuX389gyhvO2BlfpHbsUsr9UsrfzH8/DuBpACek3W45MjVF82VnMYeHytgV\nra3cFK812LEzvvCasQshTgRwJoBf+NxuuaAqYrKIY/SMHeCMtdaYmQH27QNWraKfWdiZNHgT9vkY\n5lYAfzXv3KuOLIU9GMVwU7y22LuXBiY1NtLP/GBn0uBlzVMhRD1I1P9bSvm/Ye/bunXrK993d3ej\nu7vbx8fnhhL0LGrZg8LON3ZtoefrAJ3/wcHS7Q9TOnp6etDT05NqG74Ws/46gKeklF+0vUkX9kok\na8fOGXvlMzMDXH898NGPxvu9oLBzi612CZrebdu2xd6Gj3LHiwG8B8BlQojHhRCPCSGuSLvdciTP\njJ1v7MpkcBD45Cfj/57ecQpwxs6kI7Vjl1L+FECdh30pe6am6IbLI2PnKKYymZ4mQY5bEtvbC5xx\nRuFnPv9MGnjkaQymp4HFi/PJ2DmKqUzUQz/uw3/7dmDTpsLP7NiZNLCwx2BqioQ9r4ydb+zKQz30\n4zyUJyaAJ56gKZsVfP6ZNLCwx0A5dq5jZ8JQ10acc/fIIxTDtLUVXmPHzqSBhT0G09NAVxfXsTPh\nKMce59z19ADByl9+sDNpYGGPgYpiOGNnwkji2B94wCzs/GBnksLC7sjsLM0T09aWT8bOjq0yiZux\nT04Cjz8ObNmy8HVusTFpYGF3RAlvVuWOwYy9uRk4coQeJkzlENexP/IIcPrpC/N1gB07kw4WdkeU\n8DY355OxC0GfxTd3ZaGuDdfzZsrXARZ2Jh0s7I4o4W1qyidjBzhnr0TiRjGmfB0oRHFSets1poZg\nYXdERTFZOnY9Ywc4Z61E4kQxk5PAY48V5+sArX27aBFw9Kjf/WNqAxZ2R7KOYoIZO8AdqJVInHLH\nn/8cOO00oL3d/P/8YGeSwsLuSB6OnaOYyieOYw/L1xWcszNJYWF3pFQZO9/YlUWcjD0sX1dwi41J\nCgu7I3lUxQQzdr6xK4+pKYpWoh7IY2Ph+bqCHTuTFBZ2R7KOYkwZOzv2ymN6GjjmmOgH8q23Apdf\nDnR0hL+Hzz+TFBZ2R/KuYwfYsVciU1PAkiXR5+2mm4D3vc/+HnbsTFJY2B3hOnbGBRdhf+klYMcO\n4I1vtG+LH+xMUljYHeE6dsaF6WkSdtt5++Y3gXe9C2hstG+LHTuTFBZ2R7iOnXEhyrFLScIeFcMA\n/GBnksPC7oiKSriOnbER1Xn68MN0DZ1zTvS22LEzSWFhd0RFJVzHzthQjj3svKlOU5eFrlnYmaTU\nl3oHKoVS1bHv3ev/s5jssDn2yUng+9+nhatd4CiOSQo7dkeyjGLm5oCZGaChYeHrHMVUHraM/Sc/\noQhm5Uq3bbFjz5evfhX48Y9LvRd+YMfuyPQ00NlJwu47ipmepgqJYPOco5jKQ1XFmIR9715gwwb3\nbbW2AqOj/vaNsfPQQ9S5feWVfra3axdtb906P9uLAzt2R1QU09Tk37Gb8nWAm+KViC1jHxsjc+AK\nO/Z8GRmhidl8zYF/443ADTf42VZcWNgd0TtPp6b8LoAQJux5OfYHH8ym36AW0aOY4DUyOmqfQiAI\nP9jzZXgY2LOHBpD5YHQUmJjws624sLA7osRXLYAwM+Nv21NTxR2nQD439k9/Clx2GfCzn2X7ObXC\n9DRNAiZE8SIZ7NjLm5ERYPNmcu0+GBsr3YOZhd0RfQCR75JHm2PP8sIYHASuvhpYurR0zqLaUA9p\nU2srrmPnPpZ8GR4GrrqKplP2ATv2CkAvRzRVxnz848APfpB823lHMVICH/wg8Ja3AJdeysLuC3Wd\nmB7K7NjLm5ERFvYFCCG+JoQ4IIT4rY/tlSO6YzcJ+7PPAnffnWzbpeg8/fKXgd5e4HOfA9raWNh9\nIGWhwsl07pJk7CzsBWZmgCuuyGaB77k5Oj/nnUf39q5d6bdZ8cIO4BsAftfTtsoSXXxNwj46Cvzi\nF8m2HZaxZxXFDA4Cf/d3wLe/TX9TWxswPu7/c2qNI0dI1BctMre2kjh27jwtcOAAjQXIoqN/fJzO\nWX09tWB9uPaKF3Yp5cMAhnxsq1zRoxhTxj42BjzxRLIbMe8oZt8+GiSzcSP93N7Ojt0HeqvO9FDm\njD0dahT2yIj/bQ8PA11d9H13t58O1IoX9logKooZGyOBfOyx+NuOimJ8Nz2Hh4HFiws/cxTjB/3h\nb3LbnLGnY98++jeLQVsjI4V7wpdjHxsr3X2V68jTrVu3vvJ9d3c3um0r+ZYZLlHM615HccxrXpN8\n2zoNDYWyuai5u+OguxOAhP3QIX/br1V8O/a8hf3gQep7+fSn8/vMOChhz8qxK2E/5RSKZvr6gDVr\nkm1vbo6EvbU1/u/29PSgJ2WToWTCXmlEVcWMjQFveEOyJlxYxg4UmuM+hX1kpFjY2bGnRz+PwRhF\nSjrG5TxA6cUXabRkuQt7Vo5d3RNCFFz7H/9xsu2Nj1NfS5L7Kmh6t23bFnsbPqMYMf9Vldjq2Ofm\n6Aa87LJkHahhjh3I5ubmKCYb9PMYPG8TEyT6dXXu28vbsU9NAfv30/VcjmTp2PUoBqCcPU0cMzoK\nLFuWTZTqgq9yx28B+BmATUKIPiHE+31st5ywRTGqR/3kk8m579+ffNtBsuhAM0UxLOzpCTp2XdhH\nR+Pl6wAJu+/pK2xMTVFJ4cBAPp8Xl337SHyzcOzBe+Lss4Enn0y+vdFRmlqioaE003X4qoq5Wkq5\nUkrZJKVcI6X8ho/tlhO2KEZ1igkBnH9+fNce17E/8ggwOxvvM3Q4iskG/RoJPpDHxuLFMAA15Rsb\ns1nYxYS6ppUzLjf27SPzlIdj7+ykc5YU1Z9SqnuLq2IcsVXF6J1iF1wQX9ijMvagsL/zncBTT8X7\nDB127Nlg6zxN4tiBfHP2ShD2zZvzcezt7enGdiizx8LuiXvvzWawje6qgxm77saSCHvcKObQIboQ\nk2LK2HmAUnps5Y5JHLvaTl45ezkL+9wc0N8PbNpkduwTE8BzzyXfftCxd3Skd+xK2EsxyKzqhP2a\na5KPAA1jZoZilvr5GqKwKAagKOZXv4oXlcSJYiYm6EYfSjEczORO2LGnJwvHnucgpXIW9kOH6Pgt\nXWp27LffDnz0o8m3b7onxsaS92/ows6OPSVS0hwPaZ60JvQbFrBHMUuXAscdBzzzjPv2oxy7LhCq\n3jyNY+eMPRt8Z+xA/o590aLyFPZ9+4AVK8hVmxz7wEC6Tt+gY29spAqmpP0bSthbW1nYU9PfTyLo\nW9iDC03bHDsAXHhhvFaDSx27Qgl7WsfO5Y7EzAxw//1+tqUbgGBLK03Gnqewr1pV3sLe2Wl27IOD\n9JWU4D0BpItj2LF7RM3IZjoZk5PAO96RbLtBxx7M2IMjCi+4APj5z923HyeK8eHYg81O1SooRb1t\nqXn2WeB97/OzLVu5YxrHnmfn6bp15S3sYY59aCidsAdbsUB6YeeqmHkOHkw3tF0taWU6GQcP0nzp\nUeI1Nwf8/d8vfC0ovCbHrt+0p54K7Nzpvt9JohifGXtdHX1+Lc5LMjmZThB0bFFMpWTs5S7sUY49\n6eAqk2NPUxnDVTEaX/wicN11yX9/1y7q5DQJ++goNbuj3M/oKA2p1pc1ixvFHHccPUhciVMVc+gQ\nPUSSOnY14CUY/dRqHHP4MP3dPmrFbVFMpWTsStjLrfUW5diVqCcthczCsbOwzzM8nK5k6aWXgA0b\nwoVdfYYNdRL0iydO5ylAwh6n5RF8cOiYopgNG5I7dtVJJAKTP9SqsKtj62O0ZdCxV2LGftxx1IGa\nRa14GqIcu7ofkrS+jhyhr+CEXb6EvebLHcfG0gv76aebT4Z6LWrUmhI3/QEQdNSmOnb9pj322HjN\nwuCDQ8fk2DdudBP2668vvtCDMYyi2oR95046D1deCXz+8zRXvgl1bH3EMcFyR19VMXlm7M3NJKDl\nFseoNQQ6OigeCd5bg4N07JOcR+XWg2YnTRTDVTEao6PACy8kz8l27QoX9riOPSjsUVGMftM2NNDP\nrq46bufpxo1uUcznPlcsaKYmJ1B9g5T27gXWr6d1XXt7gUsuMS93xo69QLkL+4oV1B/U2lp8rQ4O\nAiedlOw8BksdFRzFeGJsjC6ul1+O/7tzczR/8qtfbRf2JI49bhQDxMvZ43SeHjzo5thnZuhmOHBg\n4eumTiKg+gYpjY0By5cDf/iHwJe+RPGVyc0p0TQJwshIvEgti4w9787TchR2KQvCDphz9qGh8HMc\nRVgrlqtiPDE6Shdykjhm716aTW3ZsnSOXTkBWxQT1XkKxBf2OHXsGzZE/x1799LDziTstRDFBIU0\nrEmsxNckCP/+78BnPuP+mbXi2GdmgIcfzmefABLxxsZCBt7ZuVDYDx8m8T/hhORRTJjZ4aoYD4yN\nAWeeCTz/fPzf3bWLevQ7OsydKz6jmKg6diCesNsydlMUs24dXcx65U6Q3bvp31oW9vb2ws9hC4Pb\nHPv+/fGa9rWSsd9/P/DmN+ezT8BCtw4UT907OAgccwz1qSSJYrJy7Czs84yO0jzISRz7Sy8BJ54Y\nfjJGRyn7ziKK8eHYXaIYKenCXbo0vOxLESbstoy92oQ96NhNAjk5SeJpEoT+/nhlpfoApaYmqrRQ\ncwZVk2O/9146Xnktp7h3b7Gw69f+0BC11o85xq9jTyrsyvQ1NbGwAygIe1rHHlYVc8IJ+XSeAn6F\nXd3YIyP0c2MjXci2nH33bsriXTP2ahf2sLKzw4dpGH2YsMcpK9WvEyEKojw3R8dWb0G4Uo4Z+333\n0fuefTaf/Qo69mDJo+7Yk2bsPqMY/SFe8+WOanm5M89M7thtwj46SgvTujj2xkb3jF1KOvnBm9ZX\nxq43xQ8dIrcOkOu2PaT6+oBzz63dKGZ83N2xr1plFoS4wh5s2SlRnpig8xhnWTxFuTn24WHg6aeB\nq64qnbAHHbsSdhfHvmMHcO21C18La8Umdey6sNd8uaNaXm7jRlpUN27Jo4pimppIbIMjCUdHgdWr\n3TpPV660RzF6xh520y5d6idj1wVJF3YXx17Lwh4nilm92uzYDxyIF8UEH9DqoZw0X1fbCBP2uTm/\nJaouwt7TA1x0EXDGGfFmME2Di2NXUUxUxr5jB3DrrQtfC3PsaYRdne+aj2JUTt3eTsITt+RRRTFC\nmE+IcuwuUUwwsrFFMaaOUyCbKCaOY1fC3t+/cHh4WJ5YbXXsrsIeFsWotT/TOvbJyeT5OmDvPL3v\nPuBtb0u2XRNK2I85ho6L6YFy773A5ZfTEnXl4tiHhtwd++AgmcCZmcJrvh273udW88KuC+SGDfHi\nmJkZ6mBZvZp+tgm7SxSzapV752mYG3MVdimpk82lKiauY9+4kbar/81xHPvhw8C73x39N5Qjccod\nTVHMwAAJhRpu7kLQAChhT+PYbRn7vn3UuvWFEnYhgOOPN7t2JeybN5eXY3fN2AcGSC96ewuvZZ2x\n17Sw60+5jRvjdaDu3k2DURob6ecwYXeJYsIcezCK0YXd5MZchf3oUVqZaVHImdBv7IMHabsACXvY\n33L4MIn58uX0pccxYcJuGqC0bx/w7W8vdDeVQhLHrrds+vvp2EW1jHTCMva0jj1M2AcH6dr3NWGX\nXtVjimP27aMS0LPOolGevb3uD700uGbsyuzYYlwl/Lq+ZJmx17yw645948Z4jn3XLsrXFWkc+/h4\ndBRTX0+uZmbGHsUcOhR909nydaDQOpibK45iwhz7nj3UT7BoUbGwxyl3VPFEXmVtPnGtipmcJEFo\nbFzozpSwR7WMdPLO2AcH6drwMc+N6pdS16JJ2O+7D+juLkzzvHq13xZDGFGOXZU71tfTebZNYDY4\nSMdU15csMnYl7M3NZN5sS2XOztrHpCShbIRdd74bNsRz7KoiRtHZufCESEk/ByMWEybHbhJfJbhh\njr25mcQiapY8W74OkDirFkIwign7W3bvpocYYHbsruWOStjjTEFcLsQZoNTSUpzP9vfTKGbbAzRI\nWMYedo24YMvY1f6qMQtpOHKExnmolqNJ2FUMo8gjjpmYINHTr9kwxw5E5+yDg8A557g5dh9RjBDh\n157iP/4D+NSn4n+OjbIR9rSOXRf24OjTw4fpou3qIiG1PR0nJkgMp6YK7zOJry7sYW7MJY6JEnag\n0KR3dey7dxf6G3RhV/PRm+qpbcLe32/fv3IkTlVMa2vxqMUDB0jY4zj24BKHehSThWNX+7VnT7Jt\n6wT3PSjsUhYLex4dqMqt6zMvhmXsQHTOPjhIK5zp+hJWUNDURC3luHFT8HxHlTzu3+/fPJWNsOuu\n5qST4pU8qlJHRbAJpZ6gQkSP2JyYoN/X32eqM1cu2nbT+hJ25driOHaTsKsL2JTnV5uwu9axK8ce\nFHbl2G3HOUgWUYyt83RwkFqXPhx7UNhXrlwo7C+8QMbg5JMLr+Xh2IMxDBA+8hSILnkcHATOP7/g\n2KUMF/awCrsogi20qJx9eNj/Os1lI+x68yVuyWMwigkTdiC6M2xigk6E/r6wKGZ62t7MdhF220LW\nCiVKaR17WJMTMF98yvlECbtyc+XC3ByJYVtb4bUox542ilHOTnXg65+ZZefp6afn49jvuAO47LKF\nzvnkk7MX9v37i4Xd5thdophzz6WO35kZuuabmqhFbyJJHBM831HCPjJS5cKuu5o4cUxU52lQ2G2O\nfXy8WNhtUUwejt0UxdicZF+fWdjD8nUg3LEvXRot7AMDwOtfXz5rpqrBbnrLJKw57OLYXYRdlazq\nwqfOWxrHruacMbVeBwdpoFAWjl0X9qeeAv7hH4C//MuFv7N5M0UxWS6jNzxccOMK3bHPztLxVde1\nTdilpP9bsYLui95e+z0BJHPscYW9qh170Pm6dqBOT5N4rlpVeM0m7IsXx3fspigmqvMU8BvFjI7S\nxawu8ijHbuo8DSt1BMwDlAYGgFNOif4blLD4EBgfmITU5NilLEwCZhJ2VRXjEsWYWnU+HLstPszD\nsQ8O0vQBn/88RRg6S5fS/mXZua6Mlk5LC7ntI0fouHR0FEZ+2zL2yUl6X0tLoaTa1ooF8hH2kRH/\ngwO9CLsQ4gohxDNCiJ1CiI8n2UZSx75rF4m6PqQ/bRSjoiBbFKMy9rw6T19+mfZJ/Z1q/0xuKWkU\nMzm5cHsDA+TKohy7EnR90EcpMZ0TU7njkSN0POvrwx27axRjevj7yNgBuo6C52Bujs7/6aebH6g3\n3UTT67oSFPbjjqO/+x3vAN70JuB97yv+HSGy70A1TZ4mRCGO0WMYwJ6xDw7SeQYKgyCjHHvNRjFC\niEUAvgTgdwGcCuDdQojNcbeT1LHv3Als2rTwtaCw69u2dZ6q3u/GRjfHPj2dPopxzdj7+goxjNrH\nxsbiC2Z0lNyMcvbLl1NOKaX9Ila1yXqcElfY+/rs78sLV8d++HBh8YZgEz5uVUxWjh2g/Qieg9FR\nEowTTyTHHnzA33IL8JOfuH9G8Dqsq6Prt66OllgMI+sOVNWCDqLmZDcJe5hjV6OJgewdu8sYCkW5\nRjHnA3hOStkrpTwK4NsAroq7keDB2LCBeuKjeO65aGHXt21z7HqzzzVjzyuK6esrjDpVmERHuXWV\n9ba1UdY8NmaPYtR79QeFimJchL2lJR/Hfvgw8Itf2N/jKuwqhgEWOvaJCXLE7e3uUYzp4e8jYwfM\nwq4Era3NPJ/89u3xxoKYDMaXv0wjj+vrw38v6w5U08ypQGEVJb0iBrALu/4QcHXsPqpiosody9Kx\nAzgBgN4Y3DP/WiyCByOqk1NhcuzBAUqunae6O3Cpismz83T37oWOPbiPCr3jVKHimLjCPjjonrFf\ncEE+wv7ww8A119jfExycBBRiEd3Z6o5dF3YVwwiRzrGrz8zCsesitWrVwpx9eJiiuzhjQUzC/uY3\nF3dcBlEdqFkR17HbVlHS3+vq2LOOYubm6HqdnraPTo2L5Vnsn61bt77yfXd3N7q7u1/5OW7zRbFz\nJ/DWty58Larz9KWXzNsKCru6YG117DbH7jJ1bxzHft55C18Pc+yq41ShhH1kZGH1UBD9Ajx6lI7/\nmjUkgLbIaPdu4NJL81kHc2go2t2YHHJ9PX3p51KVOgILnZ7qOAXcM3bT8fExCRgQLeyrV9M5OPNM\n+nnHDhLcF16gB5leqRNn/10otWOPE8Xo712/nozIwIBfx25aVMUm7GNjBc0ZH6d96enpQU9Pj/uH\nGvAh7C8D0KVk1fxrRejCHiRuUb9i5056+uoER56OjtJADsAexegnJE4de14Z++/93sLXTX+L3nGq\n0B277SLWj7ma43rRIhKWgweLt6t/5sUXAzffbP87fOCSRwYHJymU0KrjrUodAbNjB/xEMT4ce1A8\ng8KuO/bt24EtW4DbbzfXgZtIKuwnnUSf7WJQkmBz7CqKSSLszc10X/z2tzSpWRhxhV3FuXqprU3L\n1OAo5dwXLy42vdu2bXPfgXl8RDGPAtgghFgrhGgE8C4AP4y7kaBjb2kpTH4VxsQE3YxBwYly7HGj\nmKR17O3thZWhwnCNYoaGiqMYW8aukySK0TuaTFUZirk5mjL5oouo+R/VnHzssXR1z0kdO1DcCtQd\n+5IldF3Mzi4UdtXkD16HN964cA6esChmfDz5sngKlyhGr4zZvh049dR4018nFfaGBmDtWrf+sCTY\nHLuKYvS4yDbDY9Ddb9wI/OpXfqtiTA9xm7Ars9Xe7jdnTy3sUspZAB8GcBeAHQC+LaV8Os42ZmeL\nRwqqNSNtovj88+QYgqsXJS13tHWemqpiDh8Od4fqb4hy7a5RDOCWsfsUdlUaZhIWxYEDdGEuXkw3\nzf794dt/5hmagOmLXwx/TxRK2G0PhzBhD3ag6o69ro5+Z2SkUBGjXm9vLzYDX/gC8OijhZ/DHPvB\ng8mXxVO4RDG6Y9+xA3j1q+NNppdU2AG6LqNWLkpKlGMPinVDQ/gMj8H3bthAsazPKMZk9Gyxssr4\nOzr81rJ7qWOXUt4ppTxZSrlRSvnZuL9var4A0XGMqeMUKMwtrp7aeswTp/N0ZKR4OlNFUxOJTGOj\nvWrAh7ArV+ni2G2dp2FzYij0QUpBYQ/7G/QHydq19g7Ur38duPpqKp9LOgXB8HB0K8hV2HXHDhTi\nGN2xA8VxjJR0nHWxDSt3PHAgXb4OmIVdjyBsjj0PYTfN5e8L0wAlILyOHQiPY/Q6dqAQ4fosdzT1\nublEMUmnCA6jLEaehnVAJhV2NbpM/a7ryFNTxn70KG0v+NBpbiaxi7ppo4TdNWMHzMIeFJw9e/w4\ndv0msDl2XdjXrAmvZT96lAbNfOpTVGP9nvckm8tbPchsN0ESxw7YhV1/gA4N0XHSo5gwx55myl6F\nKQoLc+z9/XSsV66MNy1HGmHPcjGJsBhLz9iDlTs2YQ86drWtMHxEMbZyRxXFVLywHzoE/Nu/LXwt\nrHPJRdiDHacK/UC51rGbMnbTDQsUhD3qps06itEF59Ahem/wRnAVdt15uUYxro799tupguLkk2mx\nhk9+ksrp4jY/fQq7ybEPDi6sigGKj7N6eEU5dnXe0jr2Y44hEdNXsjKVO0pZiGGEyM+xZ7lebhLH\nHlbyaMp3AAemAAAfxElEQVTYAb+OPW7GrqKYssvY47Jjh1nYXTq7gpgGJymCwq4OtrogTJ0rurC3\ntxcydJPwNjfTjZ3WsaeNYvSH1Pbt5uOhRp9GucckGburY7/xRuADHyj8/OEPkwB95jPh+2NC/b1J\nhV2/wfQBSkBhOHpUFGMS9jDHDqR37HV1tG/6Sla6SLW20tfAQCGGAQrC7tJZXa5RTJRjjxPF6AUB\nAJU8qrl4wshD2JVjL7uMPQ4jI1RFoQur7ygGCBf2+nq6CUwHUXcH6oT395sv+KYmutHycuwNDcWf\nFXSS995LU6sGWb6cxKilxd4fkKQqxsWxv/wy8MgjwNveVnhNCOCf/gm44YZ4S+8NDdFCy7aVqUwD\nlABzFOOasQcd+4YNxVUxprligPSOHSh+uAarQVTOrhw7QNdHc/PC/QyjHB37kSP0UNKnQlaEjTwF\n3KOY5mbqBD/++PB9yKMqRnWeVrRjHxmhDFC/mcMcuy2bGhig7eg3oI4afXr0KImnfgOHlTwG3UFX\nF90UYY49z4xdzaSnE3SSwRVuFJ2dJOi2Jidgd+xpOk//679oMqlgk3rtWhL7L3zBvl86Q0P0eUkc\nu63cEaC/9+DBwnTFClMUc+650VFMXR2JUlrHDpiFXRcplbNv314QdsA9jilHx24rE128mFqhqnpO\nxyTshw/TQyL43o98JHwudsBfVUzVd54qp6UvopHEsasYJmxUnTpQatv6+8Jy9mBpVVcXXTxhwj4w\nkE8U09lpfoDpgjM6Sjf1xRcXv08Icu1phD1pFDM3B3zta8AHP2j+/WuvBb7yFbdyOTWRWZSwRw1Q\nUgQ7T485hq6rzs6FN3vwAbp7d7Gwh/XFtLb6d+xqXnGTY9ejGCAfYc+q8zQsXwdIDPv6imMYwJyx\nqwehyyhcHV9VMbZyx7KsY4+Lcsq6sNsy9rALxhbDAIXRp6YDHUfYDxwI7zwF8olizjqLOh+D6ILz\nwAM0X0vYzbl8uT1LBMKrYlQUE8xqZ2bo9ZUr6eeuroL4Kh54gM7FOeeYPzOOax8fp7/v2GOz6zx9\n+unih6gpijn7bBIPNSDL5NjVZ/py7Oo60ucVV6xeTZOjNTUtnCzOtTKmHKMYm2Pv7KSoxjSXjcmx\nB0sdXWlpoc/RO65tJKmKKds69jiYhD2pYw+riAEKT1rTgQ6LYoIOwRbFqNfycOxCLFxIRKE6eI8e\nBe65xxzDKFwdu6mOva2NhCR44e3dS3+fcrdCFLv2736XatdtTunaa4H//M9o167fBFEZe9Jyx507\nF1bEAOYoZv16uo6UgIQ59pYW/47d1GG4ahVw550L3TpQ2VGMzbGre9rk2MOE3fTeKISIl7PbMnZT\nJ3ZVRTFLl7o79rAmjItjV8Ie3LavKAZI79jT3FBCFAZSheXrijhRjJTFFQSmOMY0ylXP2WdngR/8\noHiStiDKtf/rv9rfpzrKbDfBzExxn4rCVBWjv++YY8idmRy7ul6OHqWH/cqVdEzVMcnasesd2CaR\nWr2a/l/P1wH3aQXKMYqxOfaGBnpohgl7WBSThDiiaxL2hgYaB6PWe9Cpmjr2kRHgVa9K79jjCLur\nYw/rPLVFMVFurKuLBGR62vz/aSdP6uqiofovvxwedwDxhH1ioripb+pANQn7mjUFYf/pT6niQA0E\nsfH2t9P7bbgIu5pbxNRCcHHsgD2K2buXjmVDA71PVZzkmbGHCTtQLOxqetqoksdyjGJsjh0oTGMR\nxLQ8Xhphtzn2wUHg4x8HHnyQ+pPijsnJqo4912l7gYKw61Pn2qpiTFOmSukWxfT1mQ90XMe+fn3x\ne12FXQi6oFSZXpC0wr5kCXDbbTToxzYfyVvfGu2qVJNaj2EUppLHMMeuopjbbot264rFi6MvbBXF\nBOfb17HNtumSsQPFwq5HMX19hWmRdbHNI2O3CbuavTQYxSxZQg+hgwfDK8iA8oxioiZP6+yMl7Fn\n4djvvBP44Q+BH/+Yrs+wjnsl7MF9qJo69tHR9I5971464S4DC9IKe1TG7nLT2ka7+nDst91mj2EA\n6uy75BL7e9TxNgl73Chmbg74/vfdhV0NHLPh4thtwh6M9oKOvaODykJtjj24nqw6JmGO/eqr6din\nJUrYW1uBK64ATjut+Hdd4pgsHPv4OPDnf55sm+r3kzh20wyPwWgxDrbr7aGHgD/7M+DJJ4E77gA+\n8QmaCz+ISctUKXZbW5VEMaecUpyxxxH2KLcOFJc76sTpPD10KF0Uo7YTJuxpbiiALuLdu6OF3QV1\nvE0VBK7CrjpPH32UHr6vepXbZ0d1iAILhT3svWGDk4Box65aVybHrhYODzp2FcWEnccPftC+uIkr\n+vEPzkGuuOMO832k4hgA+PWvgQsvLG4Jh7U4XAhz7Hv3At/6lvl3Dh2K7tR1ceym49DQQK1jfS6i\nrKKYBx8EXvta+v7Vrwb+9m+La+UBc3+hcutCVImwr19PJ00tnBy33DEqXweSO/Zgxi5lus5T2+cB\nfhz7ihVmlxCXKMfukrErxx4nhgHs8YrCZZSebRplk7AHb8Jjjy0W9uZmirkmJ+NHMb7o6CisahVX\npFRlzC23kKtXfTI6WXSejo7S+TCVCt58M4mgjSjHvnx5+CIiW7bQaGdFFlHMwAANCjvjjOhtmEoe\n9dlWK76OfXSUbs6VKwsXV9woZt++QqYYhhIK185TKc1RDGAXdlfHHra8mo+M/fLL4w+8MKFchanZ\n6urYV6yg3//Od+IJe1sbiYutXjhtFBM1pQBAg6kuuKD4d1XzXhd2lyjGF0IUHq5JhP0rX6Gy0nvv\npRw+aDSyiGJUq8rUuhoeptaDjSjHfuONwJveZP6/iy4qFvYkdexA+PX28MPU+rFN06EwaZk+KV9b\nG12PvtY9zVXYpSwI7QknFIQ9brnj6Gj0YBvVXHd17NPThXUx9fcB4XPFqM+JIkvH/v730w3rg7o6\nasbu2RMdxUxPk9AFa77r6uihXVdXWIPTBZd64eFhEtg0nae2ckeABME0xFyVPNqimCwdu/q8/v74\nwr5lCz38H30UOP108/WYVtiDC4UDBUE3xZ7qWNpKgaMce0tL8XTaii1bgJ/9rPBzFlHMgw9G91sp\nTMKuO/ZFi/yWjeYq7GrkYH09DaiIcuxhI7ZGRqIjkKg69uDFZlqpxcWxu0QxtnUz02bsmzdTn4Uv\n2trohouqitmzpyDgQdasIbcetxUR1YE6NBQ9QCmuYzfloSZMjl1/2GXt2IHCmIi4IrV+PcUwav6b\noLBLme7BVF9PD8OpqYWvq3NkuvbVazbXnmZJwbPOor44ZQCyiGIeeqiQr0cR5th1g+ozZ89V2HWn\nrRz7zAxdEKYBJbbszsWx26KY4MVmWlvR5tjr6ymGKLVj902YsAcduymGUXzkI8CHPhT/s6M6UH1W\nxUgZT9i7uqjvYGamUGLnMkDJJ0kde5Dg9XjkCAlzmuX7TK42StjXr48Wdptjt9HYSOL+y1/Sz76F\nfXwceOop4Pzz3bYR5tj1sSU+c/ZchV1veihhV51dJndnK+qP49hdohjTRdTeTk2ksBu2r88+M5zt\n84BCnuyS0eVFWxuJdlDY1bqWqoTsrrvM9f0A8Ja3hP+fjagOVBXFNDdTFmkayefq2Ken6eYPa8oH\nWbKEStpWry5cq21tdDwmJvJx7FkJe9pWI2C+V6OimMsvtwt72ELWrqicfXqarpWkDwnTQ+uRRyhq\ndD1uUVEMUMGOXRdkJexh+TpgF/Yox65u/oGBYmFvbqYbUm86moRdDdlP68TChL3c3DpAF7EpY29s\npP8bHgauuw743vfiL5IRhWsUYysPcxV2U75uY8kS4IknCjEMUOjQ7O/P37GbBua4kpWwx3Xsr389\n8KtfhW8zjWMHCjm7Kg9NWmBgutbixDCAub/QFMX4GqRU8ijGtqpPmihG3fwvv1y8fX2OFUXYRaQW\nK0hDmLD7uKF809ZGLQmTI1y2DPjYx6gaoafHPDFZGlyjGCDc3duEXbWujh6NF8MAdA6Dwg4UxDYv\nx757N31WGiebhbCbatlHR6k1Gibs55xD7wnrQI3qPI1COfZDh9K1cMKE3bXjFAgvd9SjmIp27EFh\ntzl2tXLK0aPF23HptOzooINpem+w5DHsIqo1x66Ogak0bNkyKvHq6YkuN02CLYqZnqYHjnLZYQ8B\nm7ADhRssiWPfv79Y2NV6snk8pJctA559Np37BPKNYlatCo9iliyhUblhcUyazlOABiktWUKuPWmp\nI1AcxUxPU4XRli3u23CJYio6Y1ciu3Il1aNHibTpSefi2IHCDW660YMXd9hF5EvYTXXs5SrsQpib\n+tddR/Orq/nXfWOLYlTNrxK0MHcTlcuqOCauY1fHw+TY9++nh45Lf0sali2jwXlp3CeQnWM3RTFr\n1hSbGn2yrHPOCY9j0jp2gMT39tv9OvZf/5pG87pokCKqjt30OWkoWRTT1ETfv/CC3WEFD8jsLN2Y\nLie8o4NuZFPn5JIlCycKCotizj2XRlOmodIce1eXuULiwgvt60OmxXZhK4cX9V4Xxz45Gd+xqxvQ\nJOy7d9N59DFIzMayZfRAKkdhD3Psq1cXO/axMXp/fT3dX1k5doDimHvv9Svs99wTL4YB3DtPKzJj\nD/4hJ5xAK9bYHHvwgKgb16WaobMz/CZft27hXBJhwv65z8VrcpnQ5xrRKdeMPU2zNQ02x646TvX3\nJhF21YmVJIoBzFFMX18+D2i1MlK5CrurY9ed6jnnhAu7L8c+OZnumOmtkV27gOuvB/7iL+Jto6rr\n2IOxSxJhd83XATpQYe/dtImatYq0PfA21Fwjam4cRbk69rTCkRRb52lwNfo0GXuaKCbYYbxsGQl7\nHg/o5mb629Ken+CAuSw7T1evtgv7SSfR+4LTVahy1jjnyMRpp6W/ppXgzs3RSO+Pfcx9cjtFzdSx\nAyTszzwTL4pxKXVUxBH2tDWzUfgexp0VpXbs5RrFLF9Oc64EH8RK2PN6QC9bll7YVeGAakFmGcWs\nWVMcxejCLoS5A3Vigs5P2nirvp4GEaU5ZqqVd/31ZMb+5m+Sb0OnaurYg52eJ5xAU3tGOXb9gLh2\nnALl49gBs7AfPFgY5l0utLdXRhSTVtjjOvbly8mEmF7fsye/B7QPYW9ooAeREuJSRjGAOWf3ka8r\nPvvZ8MnCXFi0iK6bbduAm25KNkI3WASiFn0vy4xdCPE2IcR2IcSsECJyOQFTFAO4laeFbcOGTdjX\nraMbUo1eLIWwHziQbWdkErq7gXe+szSfHTeKCQr70aNUnWITqaTljoD5hl62jK6hSnLswMLrMYso\nZnqa4pTjj7c7dsCcs/vI1xXnn5++PLerC9i6NXodiDBWraJOdqU3hw9Ta0K/bsrJsT8J4C0AHnB5\nsymKAeJl7L4ce2Mj5X+qAzUPYQ+WPO7fXzw7Yqk5+2zg93+/NJ8dJ4oxvVe5dVvzPaljD+PYY+nz\n8nLsf/RH8UY8huFb2E1FDp2dhXmZ9MIBF2H36dh98MADwIc/nPz3u7rooaDmrjFFymWTsUspn5VS\nPgfAKQkzRTFA/Izd1bFv2GBfgEKPY9ixl564UUzwvbZFNhRJM/Yw6uspTstL2N/2NvPyd3HJwrHr\nMYKqU29spC/9Hg4K+7p1NKZFH4jo07H74KST3OcVCuOyy4D77qPvgzEMUF6OPRZhUUxWjv3tb7fP\nVb5pU2EtyKw7T01T95ajYy8ltgvbJYqxLYunUH02vhw7QPFIuVU3RZG1Y9cn3wtO3xEU9ro6KuVU\nc9sD5efYfXDZZcD999P3wYoYwG/GHjmvoBDibgC6/AgAEsAnpZQ/ivNhAwNb8aUvUedNd3c3Lr20\nG01N2VXFRLFxI83/AZTOsbOwF4gaeeoi7HEcu6+yzmXLyq+6KYoshF0XJX0OKPVZysgNDxe3Olas\noEIKVU5abo7dB695DZnNw4ftjr2npwc9PT2pPitS2KWUb0j1CRpCbMU//uPCDPS888LXLQTo5OpP\n+5ERf/HFpk00SyGQj7AHJzvav5+jGB19Ol41T5DCpSrGVdgHB/069uXLi+czKneCwu6yroCNYOep\n7tiD8zIFHTtQmGJEUY2OvaOD1kf92c/Mjl1l7N3d3eju7n7l9W3btsX+LJ9RTGTOrlbk1nnoIfvc\nI8FVb+JEMVGUMmOXkh17ENt0vHE6T22kqYoJgx17dBSjX/smYVeOXVGNjh0AXvc6ytnDOk8nJwtr\nHqQhbbnjm4UQuwFcCOB2IcQdtve7dnrqpOk8jWLVKnKC4+P5D1AaGSFX6ktcqoWwOMZl5GncKMan\nsNd6xh7WeQoUr1hWq44dKHSgmqIYVS/vY93TtFUx/yOlXC2lbJFSrpBS/p7t/UmcdprO0ygWLaLK\nmeeey9+xs1s3YxL2ubniefvTZuw+o5hLLgEuvdTPtvKinDpPgdpx7BddRCtx7dlTfAwAf5UxuVbF\n+BB2n44doDjm2Wf9OjgTwTp2ztfNmC7s0VFyb/oAIVOz1UXYk04CZuO1rwXe+14/28qLcotiasWx\nt7TQgKnbbzfroa9a9lyF3UcU49OxAyTsTzxBTek0i/lGwY7dDZNjD8YwALW2TJUYpXDslYhvYW9u\nLoz8BcI7T/W52HVWrFgo7NXq2AHK2V96yaxjNe3YfQv7449n7w6Cdew8OMmM6cIOVsTo79UfAn19\n0UPHs8jYKxFd2H0s6yfEwns1zLHrc7HrrFy5MIqpVscOUM4OhF/TPmrZy17Y9aoYKbOJYh57LHt3\nEBxazYOTzJgce7AiRhF8COzYAZx6qn377NgJ344dWNiCCus8NcUwAHVADwwUHH/WfV6l5Lzz6G+r\nGseeNoqZnqYmuM8KhE2bqL4864uosXHhjHrs2M24RjHAwptgdpb6Sk45xb79LModK5EshF2vZQ/r\nPA0T9vp6mndHjT6t5iimsRH45382Tw1RkRl72ijGt1sH6GJasiSfi0i/mdixmzE5ljAx0N/74ot0\nPKOa7+zYCX1Odp+OPSqKCTuXwMIO1GqOYgBagck08rkiHXtaYffdcQpQNrhpU/7Czo7dTBzHrg9S\ncolhAM7YFfqc7D4de1gUE+XYgYUlj9Xs2G1UZMaexG2rG1Hl676FHSBhz8MdsGOPJkzYozpPn3rK\nTdhVDnzkSOWNFvWNuh7zduymhzRQW449jJpx7HV15CwOHzaXSfkgT8c+NEQPqf5+FnYTpgt7cDC6\n2erq2BsbqeSuqSn9smuVTl7C7tJ5CrBjB2ooYwcKHV5ZOfYrrwT+4A/8bzeIupGGhijfrXXHaMLk\n2Ht7aYm1IEmEXQi6nmo5hlH4FnYVxczOUitbCXNbG7WQjhzhjD0KX449cnZHnyR122q0YBadpwCt\nGnR25MJ+6VG17Jyvh2MS9l27gBNPNL93bIyEZOdO+6IqOq2txbNH1iJZOXblttXCFEIUKmOGh80P\naaAwSEnKwmLWtUZFZuxJ3ba6YLLoPM0TdSNxvh5O0LFIGS7sKmN/4QV6ULo23dmxEyoanJryU0Ks\nHLspMlUdqFGOfe9eil2bm7MdCV6u1EzGDhSEPasoJi+UsLNjDyfo2Pv76fybmuXqJnCNYRStrbVd\n6qhQawTU1/sRUd2ABYVdXftRGfu+fbWbrwMVmrEnnczfdsFUErqws2M3E5xn/aWXaE1ME0mFva2N\nHTtA1+P+/f76emz3qepAtQn78uX0oFGTvtUiFenYk7qCanPsPLNjOCpeUVMvZCHs7NgJ38Jui2L0\njD1M2BsaqPrpxRdr17FXZMaeFL0qptId+9AQO3YbjY1kAKam6OewfB0ouPskws6OPV/H7hLFAJSz\nP/ccO/a0VISwq6qYauo8Zccejn5xRzn2wUESAteKGIAduyJPx+7SeQpQzv7cc7Xr2FesoDmP0lIx\nwl4Njl0vd2THHo7egfrSS+GOvaMDeP55mqo3jgNnx07k7dgHB4tXwgpS64590aL0C4sDFSbs7Nhr\nA13Yd+2yO/a5uXgxDMDCrujqAg4dyq/zdPdueo+tr23FCnpY16pj90VFCXuld56q5ujBgzT/NGNG\nRTGzsyQGa9eGvw9IJuwcxZCwS5lf52lvrz2GAcixv/hi7Tp2X1SUsFd6uWN9PYlKRwePfLShHPve\nvVQlESY8ao77uMLe3s6OHSiIbF5RzK5d0cK+YgUttsGOPR25TimQlNZWcnDj437yp1LS1VX5f0PW\nqGoXWwyj6OiIL+wf+lChnLKWyVPYFy8G9uwJb30pVqygf9mxp6MihL2tjToc9fknKpWuLuC440q9\nF+WNqmW3dZwqbr7ZvBKNDSUetU59PQloXlHM7KxbFAOwY09LxQj7vn2Vna8rurq4IiYKFcX090c7\n9t/5nXz2qVrp6vLr2MMm61P3bpSwq3uDhT0dFeF/q03YuSLGjopiXBw7kw6fwl5XR/0eBw6YHbv+\nbxiNjdSi5SgmHRUl7JXccapYsoQdexR6FBPl2Jl0+BR2IPxe1VdTimLFCnbsaamYKGZmpjoc+zvf\nSQNqmHBUFOPSecqkIwthHxgoFnaV57sI+xln8D2SlooQdlWaVg2O/Y1vLPUelD+dnTRKcd8+YPXq\nUu9NdeNb2FWEYqr86upyE/ZvftPf/tQqqaIYIcR1QoinhRC/EULcJoTIRHpVs6waHDsTTUcHLU59\n/PE04x+THVk49pYW83lbvNhN2Jn0pM3Y7wJwqpTyTADPAfhE+l0qRgl7NTh2JprOTloViTtOs+fd\n76Y1f33R3h5+nx5/PBcO5EWqKEZKeY/2488BvDXd7phhx15bdHTQACLO17Nnyxa/22trCxf2H/2I\nF3DPC59VMX8K4A6P23uF5mZaEJeFvTZQwsDCXnnYhL2lhe5jJnsiHbsQ4m4AeoGeACABfFJK+aP5\n93wSwFEp5bey2Ekh7BcMU12o88xRTOVhi2KY/IgUdinlG2z/L4T4EwBXArgsaltbt2595fvu7m50\nd3dH/cortLayY68VVGUFO/bKgw1Yenp6etDT05NqG0KmmA1JCHEFgH8B8Fop5UDEe2Waz1q/Hvjq\nV4HLL0+8CaaCaG+nypg1a0q9J0wcrr2WJvvikkV/CCEgpYwVYqXN2K8H0A7gbiHEY0KIL6fcXiiL\nF9MUrkxtcP/9XMNeiXAUUx6kcuyxPiilY+/tJffGnS8MU748/TTN8HjeeaXek+ohiWOvGGFnGIap\nRUoRxTAMwzBlBgs7wzBMlcHCzjAMU2WwsDMMw1QZLOwMwzBVBgs7wzBMlcHCzjAMU2WwsDMMw1QZ\nLOwMwzBVBgs7wzBMlcHCzjAMU2WwsDMMw1QZLOwMwzBVBgs7wzBMlcHCzjAMU2WwsDMMw1QZLOwM\nwzBVBgs7wzBMlcHCzjAMU2WwsDMMw1QZLOwMwzBVBgs7wzBMlcHCzjAMU2WwsDMMw1QZLOwMwzBV\nBgs7wzBMlcHCzjAMU2WkEnYhxN8LIZ4QQjwuhLhTCHG8rx1jGIZhkpHWsV8npTxDSnkWgP8H4NMe\n9qnq6enpKfUulA18LArwsSjAxyIdqYRdSjmu/dgGYC7d7tQGfNEW4GNRgI9FAT4W6ahPuwEhxD8C\neC+AYQCvS71HDMMwTCoiHbsQ4m4hxG+1ryfn//0DAJBS/p2Ucg2AmwH8ZdY7zDAMw9gRUko/GxJi\nNYAfSylPC/l/Px/EMAxTY0gpRZz3p4pihBAbpJTPz//4ZgBP+9oxhmEYJhmpHLsQ4lYAm0Cdpr0A\n/kJKuc/TvjEMwzAJ8BbFMAzDMOVB5iNPhRBXCCGeEULsFEJ8POvPK1eEEKuEEPcJIXbMd0BfU+p9\nKjVCiEVCiMeEED8s9b6UEiHEYiHE94QQT89fHxeUep9KhRDio0KI7fMFGjcLIRpLvU95IoT4mhDi\ngBDit9prS4QQdwkhnhVC/EQIsThqO5kKuxBiEYAvAfhdAKcCeLcQYnOWn1nGzAD4aynlqQAuAvB/\navhYKP4KwFOl3oky4IugwoNTAJwBS19VNSOEWAmqrDtbSnk6qA/wXaXdq9z5Bkgvdf4vgHuklCcD\nuA/AJ6I2krVjPx/Ac1LKXinlUQDfBnBVxp9Zlkgp90spfzP//Tjo5j2htHtVOoQQqwBcCeDGUu9L\nKRFCdAK4REr5DQCQUs5IKUdLvFulpA5AmxCiHkArgL0l3p9ckVI+DGAo8PJVAG6a//4mUKGKlayF\n/QQAu7Wf96CGxUwhhDgRwJkAflHaPSkp/wrgYwBqvZNnHYBDQohvzMdSNwghWkq9U6VASrkXwL8A\n6APwMoBhKeU9pd2rsmCZlPIAQAYRwLKoX+DZHXNGCNEO4FYAfxWYkqFmEEK8EcCB+RaMmP+qVeoB\nnA3g36WUZwOYBDW9aw4hRBfIna4FsBJAuxDi6tLuVVkSaYayFvaXAazRfl41/1pNMt+8vBXAf0sp\n/7fU+1NCLgbwJiHEiwBuAfA6IcQ3S7xPpWIPgN1Syl/N/3wrSOhrkdcDeFFKOSilnAXwfQBbSrxP\n5cABIcRyAJifQbc/6heyFvZHAWwQQqyd791+F4BaroD4OoCnpJRfLPWOlBIp5bVSyjVSyvWga+I+\nKeV7S71fpWC+ib1bCLFp/qXLUbsdyn0ALhRCNAshBOhY1GJHcrAV+0MAfzL//fsARJrC1JOA2ZBS\nzgohPgzgLtBD5GtSylo8URBCXAzgPQCeFEI8DmpOXSulvLO0e8aUAdcAuFkI0QDgRQDvL/H+lAQp\n5S/nBz0+DuDo/L83lHav8kUI8S0A3QCOFUL0gaZC/yyA7wkh/hQ0EPQdkdvhAUoMwzDVBXeeMgzD\nVBks7AzDMFUGCzvDMEyVwcLOMAxTZbCwMwzDVBks7AzDMFUGCzvDMEyVwcLOMAxTZfx/nlGUoaPh\n0rYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10e6a3d10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(t,rand_n)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Matplotlib generates a plot and then in a separate command you have to display it to the screen. Yeah, it's a pain but you get use to it. \n",
    "\n",
    "Let's now generate some more complex data:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "199"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = np.exp(-t/0.05)\n",
    "r_dat = np.convolve(rand_n, r)*dt\n",
    "r_dat.size"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The convolution returned an array that is twice the size of our other working data sets. We should trim that down. The syntax for segmenting an array uses square brackets:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "r_dat=r_dat[:t.size]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This means we want r_dat from the beginning element to t.size (which is 100). Let's create one more data set to work with using some more of NumPy's functions (sin and pi). "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "ss = 0.1*np.sin(2*np.pi*t) + r_dat"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What does all of this generated data look like? Why don't we make five different plots to highlight some different plotting capabilities.\n",
    "\n",
    "First, we will define a plotting space, but instead of showing the plot this time, we are going to plot some more values. First, we will plot the magnitude spectrum of the generated data with a linear scale and a dB (logarithmic) scale. This can all be done within matplotlib. Next, We can easily show the wrapped and unwrapped phase spectrum. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The 'subplot' command means that we want a 3 by 2 grid of plots and we want to operate on the first sextant. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEPCAYAAACgFqixAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXeYVNXZwH8vTUE6CEgHWRUUFFBBQV0VjdiwxAoau8YW\nNTEqfhFMNBaMEWM0agRNFDR2sAAqrFgRKdKb9LIs0nvb9/vjvZeZnZ1yZ3fKLnt+z3OfnbntnJm9\nc95z3iqqisPhcDgqNpWy3QGHw+FwZB8nDBwOh8PhhIHD4XA4nDBwOBwOB04YOBwOhwMnDBwOh8NB\nBoSBiJwlInNEZJ6I3Bfl+Pki8pOITBGRH0SkR7r75HA4HI6iSDrjDESkEjAPOB1YCUwELlfVOWHn\n1FDVbd7rjsD/VLV92jrlcDgcjmKke2VwPDBfVZeo6m7gTaBP+Am+IPCoCRSmuU8Oh8PhiCDdwqAZ\nsCzs/XJvXxFE5AIRmQ2MBK5Lc58cDofDEUGZMCCr6geeaugC4JFs98fhcDgqGlXSfP8VQMuw9829\nfVFR1a9FpK2I1FfVdeHHRMQlUXKkFVWVbPfB4cgW6V4ZTATaiUgrEakGXA6MCD9BRA4Ne90FqBYp\nCHxUNSvbgAEDym27+fnKHXeUn8/cv/8A6tdXtm3LbLsOR0UnrcJAVfcCtwNjgJnAm6o6W0RuFpGb\nvNMuFpEZIjIZ+AdwaTr7VNGYMQPefz/bvQjOjh2wbh1MnJjtnjgcFYt0q4lQ1VHA4RH7Xgx7/STw\nZLr7UVFZvhw2bMh2L4KzY4f9/eorOPnk7PbF4ahIlAkDclknNze33La7bBls2QK7d2e+7ZLQoYO1\n+/XXWWne4aiwlIUI5Cu9COSfRORrL/AsJjt2wE8/pa+/0SjPwmD5cvu7cWPm2y4Jbdvm0qULfPst\n7N2blS44HBWStAoDLwL5OeBXwJHAFSJyRMRpC4GTVfVozK305Xj3HDsWbr45Hb3dP/GFwfr12e1H\nUDZuhHbt4JBDYPr0bPfG4ag4lIUI5O9V1Z+3fk+UoLRwCgpg5cq09HW/ZPlyqFat/NgNNm6EOnXg\npJOcqsjhyCRlIgI5jBuAT+PdcM0ayM+HwgqYtOKuu5Kf4S9bBu3bZ3ZlsGVLyQfyDRugbl3o1Alm\nzrR9L70Ev/1t6vrncDiKk3ZvoqCIyKnAtUDPWOcMHDiQzz4zY+iIEblccEFuxvpXFnj1VbjoouBe\nNtu2wdatpnbJ5Mrgyy/hL3+B779P/lp/ZZCTAyO8iJRvvoH//hfOPBMuvDA1fczLyyMvLy81N3M4\n9gPKRASyiHQCXgLOUtWYc9iBAweyeLEZF9u2TXVXyzbbt9tAuWhRcGGwYgU0awb16ye3MigogEaN\nQu+HDoU337QBOXx/LNasgbVrg7cXzsaN0LSpCbD5823fvHnwxBO2Ojj7bDjggJLdO5zc3NwiRvKH\nH3649Dd1OMoxZSECuSXwLnCVqv4c72Zbt9pAU6kSrFqVtj6XSVavtr8LFwa/ZvlyaNHC1C5BVwa7\nd0PLlkXb+d//TC13/PG22khEQQH88kvwfobjq4lat7b/8c6dMHcuXHMNHHigsxc5HOmiLEQg/wmo\nDzzvF7iJdb/Vq00Y5OSUn0Fh504YNqz098nPt7+LFgW/ZtkyaN4c6tULvjIoKLA+Dx1q77dvN/3/\n229D1aoh76R4rFljg/qePcH76uOriapUMaH0g/c0NGwITZpUvEmAw5Ep0h5noKqjVPVwVc1R1ce9\nfS+q6kve6xtVtYGqdlHVzqp6fKx7rV5tg9XRR5efQeGLL+C665IP+opk1SqbMScjDJYvN2FQt25w\nYZCfb+e/+qr5+X/1lX3fdesGVzetWWN/S2K09oUBmND/6CM4/HAQMXfT8vJ/dzjKG2Uh6OxwEflW\nRHaIyD3x7pWfbwNNp07RVwbJBCnt3Alt2oTSH6SLMWOsrdmzS3ef/Hzo3j15NZG/MgiqJvLbadIE\nPv0URo+GX/3KjtWvb3mDEuELg5LYDXw1EZjd4OOPTRhA6oTB2LElt2k4HPsrZSHobC1wBzAo0f0W\nLTLd9WGHFR8U5s+HE04I3rc5c2Dx4pCRMgjnnw/jxwc/H0wYHH44TJmS3HWR5OdD166mi9+xA7p1\nM8NqPJYtC9kMgs7SV682QfDgg/Cb38ArryQvDAoKLLahJANu5Mpg5szUCoMFC+Ccc0zIORyOEGUh\n6OwXVZ0EJNQwz5gBBx9s3iaRg8KSJcnZEWbMsL+zZgU7f+VKGDUKbrwx+Gpi2TIbGK+9FiZPhl27\n4J57zBNqyJDgfQUTBs2bmx599GjTpc+dG/+axYvNEOuvDDZsgDvuSNxO48ZwwQUmMAcNMiEEya0M\ncnJKZkQOFwbt2tnfcGHg205KgircdBPUqBHsczgcFYmyFnQWl+nTTRgcckjxgX/NGti8Of71W7aY\n/7uqCYPq1YOrbz7/3FYGHTvCY48Fu2bMGDjjDBtMp0wxv/nx481XPtkUzfn5NmNv0waeecb2LV0a\n+3xVEwatWoVWBhMnwssv27FE7YB91zfeCJUr2/t69YILgyOOSH5lsHev/Y9q17b3OTn21xcGpTUg\n//ijqc5uvtmpiRyOSMpM0FkQpk4dSKtWprpYsSIX1VzEq021Zo0NJKrs2xfJn/5kA+nFF5swOPvs\n5ITBGWeYPv2SSyCIW/q4cXD66dC5M0ydagPxXXeZZ8xTTwVr1ydcGLz0EpxySnxhsHateeTUrWsr\nkg0bTJju3GnfVax4gdWr4cQTox+rX98ETDy2bjVVXqtWyQ+4mzdDzZrmOgx2j3bt4FCv/FFp1UQr\nV5qQatQIvvoqj4ED80p+M4djPyPdK4Okyl4mYs+egZx44kAefXQgNWvmFpmlrlljg1AsP/gffoDh\nw+HXvzY1y/TpcOmlRdVEq1dHN9CqmjDo1ctWBn5KjETMnm3G7gYNbFb9ww8miHJykrNVQFFhUKsW\nXH21qaFisXixnQuhOINp0+x9vOt8NVE0gqiJ1qyxFUXDhskLg3AVEZgr6/z5Fl8AyQuDV14JGbMh\nFExXvz7UqJHLwIED920OR0Un60FnESSsQXvwwfY3cmAoKLC/sVRFL78M990HV15pQVRr1kDv3mZQ\n9P3hn3zSdOSRzJxpA9Khh9qstUePxLl3VM3Ae9hh9r5zZ+jb11RTrVqZ4Alqe1ANDdLHHWc2iHbt\n4q8MfHsBmDG3WjX47jsTKPGuC1cTRZKMMGjQIHlhEO5JFI1GjeyeQeMXnn7aUln4rF5t32FQ24fD\nUZHIetCZiDQWkWXA3cCDIrJURGrGuqev3mjWrGgAlD8D3LQp+nXz5pm//Gmnme64fXubYTdpEvLd\nHzOmuDApLIQHHjAh4nPSSUW9imbMKD6wr1plhkp/cBs0CAYMsNdVqphA+DluvHWIDRtMiFSvbmqn\nwYPNkBxvUF+0KCQMwPoxfz6cdVb8lYHvTRSNoMKgUSMTBskakCNXBpFUqWL39QV/NL74wlRhYJ8l\nPC4jfGXghIHDUZSyEHS2WlVbqGpdVa2vqi1VdUu0e1WvHloZtG1bVKXjC4NYK4P58009U6eO6cSP\nOsr2t29v6pyVK21QjxQmjz9uA8dDD4X2nXSSBWP53Hqr2QfCCV8VgLXt991/H1RVFG223qyZ7Y81\nSw5fGYAJg7ZtTWceS4hs325brNl5EANy5Mpg3ToLYAvCxo3xVwaQ2Ih83XVmKN+929oPt3H4wqBB\nAycMHI5Ish505p3zrIjMF5GpInJMrHs1bhwaUNu1MxWPjz8IRRMGW7bY7LqZ58f029+a6yRAhw7m\n9vnZZ8WvHzUKnnvO1ErVqoX2H3usDeR+9bC1a4sHdc2dG/KCiUZphUHVqjawxRoYFy0K2QzABvKO\nHW1FEWtl4KtRYhngg0QgFxQUtRl8/DEEVclv2BB/ZQDx3Ut377bV4sqVoclB+Mpg9erQysB5Ezkc\nRUmrN1FY0NnpwEpgooh8qKpzws7pDRyqqjki0g34F9A92v169AjNttu1K6q3LyiwwTeaMFiwIKTv\nB7j88tCxa64x1VHLliYgJk+2/YsWWdDVu++GhIhPtWo2sE6bZquEdeuKl5WMXBlEkpNjHkZBiGXU\nbdHCZvktWhQ/Fm1l0LFj6JpY7cRSEUFoZRDPYytcTbR2rf2Pli2zgbpq1dj3hsRqIohvRF62zNR6\nK1fawF+1avGVQePGobiLwsLQM+FwVHSyHnTmvf8PgKpOAOqISFR/ltdfD822w1cGu3ebeqdly+g2\nA19FFI2jjrKkbNOmmaePL0w+/RTOPRd6xqiu0Lix6cRVbYCMtjJIJAyCrgzmzAkFYIUTy27gxxiE\nC4NevSySONHKIJ4wqFbNDOlboirxjHA10bp1ZlsRiW+n8ElkQIb4wsAf+H1hcMwxJtT9uApfTVSl\nChx0UGz7ksNRESkLQWeR56yIck4xfJtBYaHNQOvXt1lltJXB/PnRB1Ofc86xWfFRR4Wu37Ahfu7+\nhg1NGGzbZn780VYGqVITTZli3kiRxBrYCwrMeF2rVmjfXXfZyqppUzseLXFeopUBRDe+hgex+cKg\nWjWrO7BihcVm+IPynDlE5Z134IMPTIjEI1rAoc/ixdbmypX2GQ87zAb+devMtrJhQ+j+zojscBSl\n3C6SDzrIlvsrVoRUE7VrxxYGsVYGPvXr2+DpzxYTzVJ9YeAPKOHCYNcum7HHK8DTvHnw8p1Tp0YX\nBrFUPn7kcTSqVLFVTbQBdcGCxMIg0og8Y4Z9Fj9eY+XK0D0aNrR8UTk5Jgy++caiuCNZvhyuv94i\ng2++OX77zZrFFwbHHRdaGTRubKujRYvsf1W/fiia2gkDh6MoZaHS2QqgRYJzAIoEB+Xm5tKuXS4/\n/2xpDA4+2AbzWMLgN79J3NmaNW2mX1howsCPfI1Gw4amrogmDBYvthl4vIpcvqpi40YbYL/+OrpK\nau1aOyfcGOzTsqUFw0XiF7WJhS9EwgXGkiWmLktUqjLciLxwoUVlH3SQ2VqOOKLoiqhBA/tMqjYg\nb9kS3QC9apWt3G68MX7bYN/rihhhi4sWmafYhx+GhEGbNvb/qFYtcqWXx+DBeXH/xw5HRSLdwmBf\n0BmwCgs6uyLinBHAbcBbItId2KCqq6PdLDJS1Lcb1KwZEgbRZspBVgZgxsSDDgp5HyVaGUyfHl0Y\nrF1b1I00Fv7AWr26GaK3bLH2w5kyxXTf0Qyd4bPkXbtscD7iiFDq6li0aGEz+ZNOCu274w5TJcVT\np/l99j/zqFGW0qNVK7vf0qX2nfm5hXJzze4yY4bZYJYsse8p0gAdRD0V/pljCYPFiy3C/PnnTRh0\n7BhaGdStW1QY5OTkcvbZuVzhPY2u7KWjopP1oDNV/QRYJCILgBeBW4Pe/9BDTRj4eupoK4ONG22Q\nbdo02D39eySjJmrYsKgBOYhXDIQGVt/NcXUUEegLg2iEz5I/+wyuuspeJxIGt95qeZpef93ejx1r\nsRb33hu8z2DC5/DDzT139mzb2rcPnfvUU6beatPGBuQffzRbxfbtRe+5apXZAoLQpIl979FsHosX\nW2BhYaE9F/7KYNGikFtptM/hcDgykKhOVUcBh0fsezHi/e0luXe7dlaOsVq12MJg3jxbFcRyhYwk\nGWGwZo0NKG3aFF0ZJCsM/BiG1auL2xmmTjVVTDQaNw6lZ1iyJBTRvHy5DYqxOOUUC5I74wz73v7y\nF4uODlJoPlIYdOtmAmDWrOLCwKdNGzu+d29IcNaoETqezMqgShXrc35+UVXYrl1mNG7e3FYP06bZ\n4L9nj6mNcnKKuuc6YeBwFKXcGpABunSBL7+0WXEsA/LcuaY6CYpvRE5mZRApDIK4SEJoQPLTNsRa\nGUQzHoMNjA0b2sC4bJmpnNavT2wzADjySHjzTUvWt3Yt+9QlQfsMJgzatrWBdskS+Omn6MLgkEMs\nRcTRRxdfRUFywgCiq4qWLrWVUpUq9nfbNhv8TznF1HlffulWBg5HPNImDESknoiMEZG5IjJaRKLO\nlUXkFRFZLSLTkm0jJ8cykU6fboNJuDeQT6JI4Eh8gZKMMGjbtmQrA98zJ5aaaN06G9g7dIh9D99u\n4NtKFi5MrCbyOflk06+/8ELIyyYRfm4g3yjctq2tbFq3tmjjaMKgUiWzK3Ttat9LpBtuMmoiiC4M\n5swJxVX4KsFGjeyZuPtuWx04YeBwxCadK4P7gc9V9XBgLPBAjPOGYmUxS0SvXjbgn3dedDVRssLA\nv0eiAb1OHdN9r1plA92WLaEazMmoidavj70y+Oorc82MF7nr2w2WLjUBtWCBCYfIqOlY9O1rht6g\ndOpknkPr1tkgX6+e7W/f3oRaLMGVkwPHHx9Kpx1OSVYG4e6lgwaZt9i119r7pk3tu/W/t9tuMyHm\n1EQOR2zSKQz6AK95r18DLoh2kqp+DQSs0BudZs3sh58qYVBQYAOdn0c/GiI2+M6bZ39r1gy1nazN\n4JdfbEUSKQzGjzc1Rzz8gXHZMjv3++/tXvH6XhqOPtpWH1OnFrVvdOhgnyeWF9Wrr1oakFQJg/CV\nwd//bt/V1VeHjocP/LVqmY2kV6/QPicMHI6ipFMYNPJdRFU1H4gTz5safBXPqlXmNrlnj7mVxksL\nEUmtWjawBtH5+8LAj3721R8l8Sbq0KG4MPjyS1PlxKNpU1sVrFxp5375ZTAVUUmpVs1sGG++WTT2\noUMHWx3EMtQ3bBiqvBYuDFTt/1UaYbB+fdG+NG1aPHq8Y0dz4fVp3Di2i6rDUREplTeRiHwGhOcR\nEkCB/4tyepzKu8GIDDrLjdBv+DaDH36wIK5hw0yNEZ6WIRG1a9vgGmQwb9jQfOhLKwx27zaDrl/g\nvqDAPHvmzLGI2ng0bWqfs2FDG4ynTrX0GunkhBOsWFB4kNhFFyXuKxS3GWzcaAImMr4iHuHCYMcO\ncyUNH+h7944fMAiwbFke+fl53Htvcm07HPsrpRIGqhrD6RE8o3BjVV0tIk2AOCVJgpGoPOEBB9hM\n84cfbHAcODA5FRGY4Jg4MfjKAEouDOrVs1ntrl1w5pmm6ti50wyhrVrZ4JrI3bNZM/u87dub2kY1\nvSsDMGEwaFBRNVH16sFWYJErg2RVRFA0vmLDBvsew1cktWtbmvF4nHaaTSZ69oQ+fVzQmcORTjXR\nCOAa7/VvgA/jnCsEKHmZCBEbzL/+Gh580NwdSyIMklETgQmD8EGuJDaDDh1sYJw504TB3XfD73+f\n+B5Nm5pqrGVL2ypXzowwgPi5l2IRTRgk40kEoZWBqglT34idLCeeCN9+W7JrHY79jXQKgyeAM0Rk\nLlbP4HEAETlERD7yTxKRYcC3wGFeyctrS9No7do2Uz7rLJttx3PLjEaywsBXcUSuDJKJM1i71tQa\nO3daMrfjjoObbrJUDonwvYZatDAjeqtW6RcGTZqYQEj2u4WQMNi1y4zKydoLIJTuYtOm0gmDE05w\nwsDh8ElbBLKqrgN6Rdm/Cjg37P2VkeeUBt8bKCfHdOnhuuQg1K5tg3JQYVC/vq1ISmMz8COoGzWy\nHD6/SsLRtl49UyW19NIBnn66GUvTTUkHUf97+ukncwU94YRgtoZwREIVz9avD/a/ika3bhbUt2tX\nya53OPYnshp0JiLNRWSsiMwUkekicmdp261Vy+oSVK5sA2WyLpa+sTkZYQChQU7VZqz+7DUe1avb\nwLZ9u53fuLG5QHbtGry/IrY68IXBSy9ZZHZZxV8Z/PyzGc1/+CF5NRHYd7V6delWBrVr24osaMU5\nh2N/JttBZ3uAe1T1SOAE4DYRSSJ5RHFq1YqflyfI9RBMGDRpEvKrr1PHBrmtW22mnqjEo0+9ehYQ\nJWID3M6dsRPTxeLXvy7bAiCccGFw7rnw4otF/f+DkgphAJaGI17lNoejopDORHV9AD9k6jUgDxMQ\n+/DiD/K911tEZDZW5SxGPazE1KmT/GAaTjLCIDfXZrf++YsWBVcR+fgrC7AB7vDDLYAtGZ54Irnz\ns4kvNBcssMpr119fsvs0amTqwNIKg/vvT3yOw1ERSKcwKBJ0JiJxg85EpDVwDDChNI0+8USwWgKx\n8NU7QYRB5coh46evJiqJMPDdIhs3Lj8z/JJSt659Rz//HIoYLgn+ymDz5sRJ+RwOR2LKRNCZiNQE\n3gF+p6oxF+2Jgs4gekWwZEhmZRCOLww2bEheGPhJ4q64woKo9mdq1LAgu9mzEweGxaNxY9P179xp\n+ZKSJS8vj7y8vJJ3wOHYz8h60JmIVMEEwX9VNV4sQsKgs1RQUmFQt655BpVkZeDbF0oyqJU3ROy7\n2ry5dC6w/spAtWRqosjJhAs6c1R0ykLQ2RBglqoOTmNfAlO1qhmAkxUGHTpYAZe1a5MTBr4BuSJR\np46t4KKV8gxKqgzIDofDSKfN4AngfyJyHbAEuBQs6Ax4WVXPFZEeQF9guohMwVRJ/b3qaFmjRYvk\n7Q5161rk8PjxyQmSG24IXktgf6Fu3ZK5k4bjG5APOsgJA4cjFWQ16ExVvwHK3FA4e7Zl2EyWE0+E\nTz6xVM1BiVYMZn+nbt3S2QsgtDKoV88JA4cjFWQ76OwAEZkgIlO8oLMB6epPMpREEIBF0y5fnpya\nqCJSt67Vry4NtWpZMaH8fCcMHI5UkNWgM1XdCZyqqp0xt9LeInJ8GvuUVk480f46YRCf++6zQLnS\n4AfpqboU1A5HKigLlc62eS8PwNRWpa57kC1ycswY7IRBfI47LvnkdNFo1Kh4+mqHw1Eysl7pTEQq\necbjfOAzVZ2Yxj6lFRG44IKSpXZ2JE/jxk5F5HCkiqwHnalqIdBZRGoDH4hIB1WdVZp+ZZN//zvb\nPag4NG5sHkUOh6P0ZD3oLOxem0RkHHAWEFUYBIlAdlQcGjc2g31JcBHIDkdRRDU9KnoReQJYp6pP\niMh9QD1VvT/inIbAblXdKCLVgdHA46r6SZT7abr66iifPPssfPcdDB9e+nuJCKrqrA+OCks6hUF9\n4H9AC7ygM1XdEBF01hEzLlfytrdU9dEY93PCwFGEpUut/KVfhrM0OGHgqOikTRikGicMHOnECQNH\nRSed3kQOh8PhKCdkNQI57NxKIjJZREakqz+lIVuGxmwaOCviZ3Y4KjLZLnvp8ztieBCVBSriwFgR\nP7PDUZEJJAy8Wf6RItJWRIIKkEARyCLSHDgbcB76DofDkSVixhl4ap3bgCuAasAa4ECgsYh8Dzyv\nquPi3Dto2cu/A/cCLomDw+FwZImY3kRedPF/gJGquiHiWFfgKuAMYG/4IUIRyK+qav2wa9aqaoOI\n+5wD9FbV20UkF/i9qp4Xoz/OlciRVpw3kaMik844g9lAblgE8jhVbR9xzl+BfsAeoDpQC3hPVUtR\nKt3hcDgcyZJQ/y8i74nIOUnYCnwSlr1U1f6q2lJV2wKXA2OdIHA4HI7ME2SAfx64EpgvIo+LyOEB\n7/0EcIaIzAVOBx4HK3spIh+VqLcOh8PhSAuB1USeQfkK4EFgGfAy8Lqq7k5f9xwOh8ORCYK6ljbA\nVD43AFOAwUAX4LO09axo+2eJyBwRmeclvctEm81FZKyIzPRKct6ZiXbD2s9KIJ6I1BGRt0VktvfZ\nu2Wo3btFZIaITBORN0SkWhrbesXLqjstbF/gIEmHY38kiM3gfeAroAZwnqqer6pvqeodQM10d9Cz\nVTwH/Ao4ErhCRI5Id7uYUfseVT0SOAG4LUPt+mQrEG8w8Iln7D8amJ3uBkWkKXAH0EVVO2Euz5en\nscmh2PMUTjJBkg7HfkeQlcGzqtpBVR9T1VXhB1T12DT1K5zjgfmqusRTSb2JBbSlFVXNV9Wp3ust\n2KDYLN3tQvYC8bwCQyep6lAAVd2jqpsy1Hxl4CARqYJNPFamqyFV/RpYH7E7UJCkw7G/EqS4TT0R\nuShi30Zguqpmos5UM8xG4bMcExAZQ0RaA8cAEzLUZLYC8doAv4jIUGxV8CPwO1Xdns5GVXWliPwN\nWApsA8ao6ufpbDMKQYMkHY79kiArg+uxGWpfb3sZuA/4RkSuinehiBwgIhNEZIqndx8Q5ZxTRGSD\npx+fLCLRSmZmDRGpCbyDDYpbMtDeOcBqb1Ui3pYpqmC2oH+qahdsYL4//iWlR0TqYjPzVkBToKaI\nXJniNhqJyIUicpuIXIcJu3gEDsBxNgjH/kAQYVAVaK+qF6vqxUAH7IfSDRMKMVHVncCpqtoZm1n3\nFpFos/rxqtrF2x6JOLYCaBn2vrm3L+14Kot3gP+qarE4iTTRAzhfRBYCw4FTReQ/GWp7ObBMVX/0\n3r+DCYd00wtYqKrrVHUv8B5wYipuLCKnisho4GOgN3AI9gzfARwmIg976rHVItLYuyZhmdYInA3C\nUe4JoiZq7i+fPQqAFqq6TkQSupWq6jbv5QFee9FmXPFmvxOBdiLSCliFGRavCNDvVDAEmKWqgzPU\nHqraH+gPtmrCUnRkJBDPixZfJiKHqeo8LD4kE0bspUB3ETkQ2Om1OzFF9z4buFFVl4bv9FR/IzHv\nuDMIBUk+QYwgyVio6tfe8xlOH+AU7/VrQB4ZWGU5HCUlyMogT0Q+EpHfiIj/I8kTkYOADQmu9V0k\npwD5wGeqGu1HfoKITBWRj0WkQ/gBb6Z4OzAGmAm8qaqZ8HDpganFTvPUXJNF5Kx0t1sGuBN4Q0Sm\nYqqUv6a7QVX9AVuFTAF+wiYHL6Xo3vdGEQTDgG+Bw4BngdpYUGSxIMlSUMQGATgbhKNMkzDoTEQE\nuAjo6e36Bng32RqU3lL8A+B2VZ0Vtr8mUKiq20SkNzBYVQ+Lcr1LVOdIK6VJVOetDEZ6rrGIyLpE\niRrDjrln25FWgjzbcVcGIlIZyxf0rqre7W3vlKQYseeiOA44K2L/Fl+VpKqfAlVFpH6UW6CqWdkG\nDBhQodrdnz7zrbfeSs+ePTnuuOPo27cvffr04YUXXqBfv35ceeWV+85LA0nZIDL9Pb/1lnLggQPY\nudM9X/t720GJazNQ1b0iUigidVR1Y+C7eohIQ2C3qm4UkeqYbvbxiHMaq7ec9ozLoqrrkm3L4YjG\nuHHjmDVvjEVhAAAgAElEQVRrFjt27KBZs2YUFBRQuXJlbr75Zjp16pTKpiI9v0psg8gEs2fDjh0w\nZgyce262e+MoCwQxIG8BpovVN9jq71TVIOkZDgFe86KIKwFvqeonInKz3UJfAn4tIr8FdgPbgcuS\n/RAORywOPPDAfX9btWpF5cqVARARqlatmpI2PBtELtBARJYCA7BJz9ueG+sS4NKUNJYiZs2CFi1g\n+HAnDBxGEGHwnrcljapOJ4proqq+GPb6n8A/S3L/TJGbm1uh2s1m26lut6CggKeffhpV3fcaTC2z\nZs2alLShqrFiInqlpIE0MGsW3HFHLo8+Ctu2QY0amW1/f3m+ykvbQQiUtdRT8bRU1blJ3VzkAGA8\nVjazCvCOqj4c5bxnMR/wrcA16qWBiDhHk9F/ORwADz9c7HErwoABFgcpImiWKp1l+tneswdq1YJ1\n6+CCC+D66+HSMrVucaSSoM92EG+i84CngGqq2kZEjgH+rKrnB+xIDTVPocqYJ9Kdaq6E/vHemIfR\nOWIZMgeravco93HCwJE2KpIwmDsXzjkHFiyAoUNh5Eh4r0Rrf0d5IOizHURNNBDLBZQHoKpTRaRt\n0I5o4qCzPlitZVR1glgK5X1GZYejNNx5Z3zT1rPPPpvW9kVkMZbLqxBzpshoXq1ozJoF7b0CtBde\nCHfdBRs3Qh2XMKNCEyTobHcUT6LCoA0ECDqLTES3ggxlBw3C1VdDilTLjizQtWtXunbtyo4dO5g8\neTI5OTnk5OQwdepUdu3alYkuFGK1wDuXBUEAJgw6eKGddevCqafCBx9kt0+O7BNkZTDTSxpWWURy\nsAjVb4M2oKqFQGc/6ExEOmhY0FlZ58MPTSD0KrOmQEc8fvOb3wDwwgsv8PXXX1Olij3yt9xyCyed\ndFImuiAELCKVKWbNgjPPDL2/4gp47TXwvipHBSWIMLgDK3W5E0ucNhr4S7INqeomEfGDzsKFwQqg\nRdj7mInoBg4cuO91bm5u2q3zO3bApk2mW3XCoHyzfv16Nm3aRP36Fs84atQoFixYUOSZShMKfCYi\ne4GXVPXldDeYiFmzTDXk0707/P732euPo2yQUBh4Ov8HvS0pggSdYcE5twFviUh3YEMse0EGfrhF\n8NVD8+dntFlHGrj//vvp3Lkzp556KqrK+PHj+dvf/rZv5ZDI66gU9FDVVSJyMCYUZqsV1ylCpiY6\ne/eaAfmIsJp9zZvDL7/Y5McLy3CUY/Ly8sjLy0v6uiDeRIcBfwBaEyY8VPW0hDcX6YhlbAwPOns0\nIugMEXkOWzFsBa5V1clR7pVxb6LJk6FrVzjvPBiR0UrEjnSQn5/PhAlWn6hbt240adJk37FMeBOJ\n1fPYrKpPR+zP2LO9cKHZCJYsKbq/XTv4+GM4/PCMdMORQVLpTfQ28C+swM3eZDoRJOjMe397MvfN\nFAUF0KSJqYkc5ZPFixfTunVrAJo0aUKfPkUrpqoqK1akpzyGiNQAKqnqFi/L75lA2pYgQQg3HofT\npg0sWuSEQUUmiDDYo6ovpL0nZZCCAjjxRPjkE1tee5kMHOWIe++9l8LCQvr06UPXrl05+OCD2bFj\nBwsWLGDcuHF88cUX6VQRNQbe97KSVgHeUNUx6WosCLGEQdu2JgwcFZcgwmCkiNwKvI8ZkQEIkkxO\nrLD7f7AfRSHwsqo+G3HOKVgSr4Xerve0eLWzrFBQAK1aQYMGsHy5vXaUL95++21mzZrFG2+8wZAh\nQ1i1ahU1atSgffv2nH322Tz44IP78helGlVdhFX4KzPMmgXRnKj8lYGj4hJEGPgOZ/eG7VMgSODZ\nHuAeL1CtJjBJRMao6pyI88YHjWjOJGvWQKNGpk+dP98Jg/JKhw4dePTRR7PdjTLBrFlw883F97dp\nAz/+WHy/o+KQ0P9ZVdtE2QJFIKtqvp9nSK2Y/GyiB5RlJQ1AIgoKTBjk5DiPIkf5R9VSV/vRx+G4\nlYEjoTAQkRoi8n8i4nv+5IhI0klvvZqzxwATohyOWfYym4QLA2dEdpQEETlLROaIyDwRuS+bfVm2\nDGrXtqjjSJzNwBFETTQUmASc6L1fgXkYfRS0EU9F9A7wO2+FEM4kLCOqX/byA6w2bTEyHXRWUAAH\nH2xqoq+LeYY7yjMl9cVOBq+Ox3NYTeWVwEQR+TCKmjQjjB4NRx0V/ViDBrB7N2zYEF1YOPZ/gsQZ\n/Kiqx4rIFFXt7O37SVWPDtSASBVMcHyqqoMDnL8I6BppoM5GnEHr1jBuHGzZApddZvpWR/lEVXnj\njTdYuHAhDz30EEuXLiU/P5/jj7d0QemIM/CCKAeoam/v/f3WFX0i4ry0P9tTplgKis8/h6Nj/HI7\ndbK0FJ07p7UrjgwT9NkOkjNllxc9rN6NDyXMqygAQ4BZsQSBXyfWe11myl6qhtREhx5qS+i9SUVZ\nOMoSt956K9999x3Dhw8HoFatWtx2223pbjYyCeNyspCEcc0ay076z3/GFgTgVEUVnSBqogHAKKCF\niLwB9MBquyZERHoAfbGymVMwgdIfaEUZL3u5dSuIwEEH2fsGDUzn6sUvOcoZEyZMYPLkyXT2pr31\n6tXLVNbSQKRLBbp7t61qr7gicQGb8mxE/vprqFTJ4oIqCiNHWqGiCy6wscqnpCrQILmJPhORyUB3\nzOvnd6r6S5Cbq+o3QNxQrbJa9tK3F/j4RmQnDMonVatWZe/evYj3q1mzZg2VKqU9megKoGXY+0BJ\nGFPJ669DYSE8EiByp00by1tUHnn+ebN1VBRhoAp/+INpKwYNsq1HDzsWOZkIGlQZ89fgef94Deta\nVf1YVT/yBYEYzePdXESai8hYEZkpItNFJGqlERF5VkTmex5FZSJIx48x8PFjDRzlkzvvvJMLL7yQ\ngoICHnzwQXr27En//v3T3exEoJ2ItBKRasDlWGLGhLz1Fjz1VOlraYwaZampg0TPl2c10fffw88/\nZ7sXmWP2bKtdPXcu/Pa3tvp7++3S3TPeymCQ5w3xIebxswY4EGgHnIp5SAzA9KCxSBh05nkQHaqq\nOV7Zy39hq5Cs4tsLfJx7afmmb9++dO3alS+++AJV5YMPPqB9NIf7FKKqe0XkdmAMNvF6RVVnJ7pu\n3Tq4/XZLm/7II2b4ffxxG6xj8eabMH06hMfW7d1rBuO//S1Yf8urmqigwBLvSZmMVkoP779v6qHK\nleGqq2DpUpg0CS65pOT3jLkyUNVLgD8Bh2NqnK8wwXADMBc4TVU/i3fzgEFnRcpeAnXCjcrZIpow\ncCuD8se6dev2bY0aNeKKK67gyiuvpHHjxqxbl34/BVUdpaqHq2qOqkamb4/Kk0/CRRfB8OGweLGp\nP556Kv41330HL71kOmSfSZPgkEMsRXUQWre29goD1zEsG0yYACefbCljdu/Odm8ywwcfmFOAT6tW\nxTPRJktcm4FXkSzpOgbRiBN0FqvsZVZrIEfaDJyaqHzStWtX37WOpUuXUq9ePVSVDRs20LJlSxaV\nsanwypXw8sswbZq9r1sXbr01sfF37lwrxDR2bKiK2ZgxRSuaJeKggywoLT8fmjYtWf/TwWOPwSmn\nxLYHfP+95VtassS2du3S2x9V2LULDjggve3EYtkyW8GdfHJoXyqEQUbK8SUIOiuTRNoMnHtp+WTR\nokUsXLiQXr16MXLkSH755RfWrl3LRx99xJnJjJQZ4i9/geuug2Zh6+dOnSwYLN6Pfd48u+6tt0L7\nRo9OThhAUbvBrl2mh542LXvP/dKlMGAAvPNO7HMmTLBqbe3apU+Vu2IFXHyxBe3VrGkqtWytoD74\nAM49F6qETeXTvjJIBV7Q2TvAf1X1wyinlMmylwUFRYNvatSAhg2de2l55fvvv+fll0MVJ6tXr84H\nH3zAIYcckpb2vEI2NwIF3q7+qjoq3jULFtjgG+nRU6kSnH46fPGFDfiR7NxpK4o//hGOPRZeeMGq\nlk2dWnT2GATfbtCjBzz4IHz6qale8vOtKM6rr2Y2QvnJJ+0zTYiWxAYTUhMnQrduNmFLhzAoLDQj\nfMeO8NBD9h0df7x9v12KVWtJP++/D7/7XdF9TZtatbpdu6BatZLdN+3CgARBZ5TRspeRNgOwWdPi\nxU4YlEeaNm3KI488Qr9+/QD45ptvOOaYY/Y9U2mqafB0ZFWzWCxZYobAu+6ymJZIevUyY3A0YfDz\nz/ZMtmljtQrGjLFBsnt3m8QkQ5s2Vg3t009tlTF5sk2C1q6FgQOhd2+7f61ayd23JKxcCcOG2WDf\nqVP0gW72bPudNmgQfWXwzjs2aAdxrY3FM8/A9u3mvunPxs84Az77LPPCYO1aswVFrviqVDH70LJl\nJhRLQpBEdSIi/UTkIe99Sy9SOCFhQWenicgUEZnsJe66WURuAlDVT4BFIrIAeBG4tWQfJbVE2gwA\nWrYs/VLMkR2GDx/OmjVruPDCC/e5mPrRyGkksH/Lscead8gDD0Q/3quXrQyiqSbmzYPDvGxel11m\ng3hJVERgE55vv4Vrr7UYhYYNbX+DBvDsszY7Pv98GxzTzaBBNiM/9FDbfDtKOL6KCEwYRLqXjhhh\ng/mmTSXrw7RpZrN4/fWiahlfGGSakSPtWahevfixUquKVDXuBryAeRPN9t7XAyYmui7Vm3U1cxxy\niOqyZUX39e+v+uc/Z7QbjgzhPV+pfF4HAIuAqVjJ2DpxztWffkrcx3btNOp5Tzyh+vvf2+tVq1Tr\n1FFt1Sr6uYn44gtVUH344ejH9+xR7dtX9eyzVQsLk79/UFavVq1XT3X5cnt/ww2qzz1X/Lzw/dOn\nqx5xRNHjhx2meuSRqs8/H7uttWtVp00rvn/HDtWjjlJ99dXixzZuVK1ZU3XbtmCfpzSsW6f6xhuq\n11yj2rCh6ttvRz/vqqtUhwwpvj/osx1ETdRNVbt46SRQ1fVeAE1CROQV4Fxgtap2inK8TFY5UzX9\nW7SVwcSJ2emTo3Sceuqp+6KPwxk7dmyJ7ykin2FV/PbtwlKuPAg8D/xZVVVEHgGeBq6Pda/33hvI\ne+/Z61j2MF9V1CnilzR3runMwWp2H3sszJhhs/hk6dLFbA8PxvAhrFzZ7AZHH21JHE87Lfk2gvD3\nv8Pll4cM6d27Q14eRKaT+v77ULEe3/jtl6hdtw5WrTI7zL33wi23FI9FWLIEfvUrO3/mzKLHPvkE\n6teHq68u3r/ate07+Oqrkq3AkuGSS8xu1KePrRwPi5rTObQyKHFG3kTSAnMFrQxM9t4fDEwJImmA\nnpg76bQYx08BRgS8V3LitBSsW2ezq0g+/VT1jDMy1g1HCvnxxx/3bV9//bXefffdeu+99+47TopX\nBlr02W0V6zegSTzb776r2rt38f09e6qOGxd6//bbqn/4Q6BblpgXX1Q9//z03HvPHtUmTVTnzAnt\nmz5dNSen6HmbNqkedJDqzp2hfc2aqS5ebK8//VQ1N1d1715bVX37bdHrp01Tbd5c9emnbca9dGnR\n4zffrPrUU7H7OXBg+r/nvXtVa9VS/eWXxOe+/LKtHiIJ+mwHeZD7Ykbe5cCjWMDZJUFurgl+CJ4w\nGBnwPkG/v1IzZ449PJHMnGnLTsf+wXHHHbfvdaqFAdAk7PXdwLA45wbq77p1NjCED36qqo0aqa5Y\nkeynLx1bt9oAumBB6u/9xReqnTsX3bdnj2rt2kUHxS++UO3Ro+h5J59s+1VtsL7vPns9aJDq1VeH\nzvv0U/vehg+395dfrvrKK6HjhYWqrVubEIrFN9+oHn10cp8tWebONZVfEMaMUT311OL7gz7bQcpe\nvgH8EXgMWAVcoKqlzIJRhDJX5SyaJxGYmmjpUlMjOcoX4ZHIv/zyC6NHj2bjxo3pbPJJEZkmIlOx\nSc/dpb1hvXpwxBGmGvHZsMFy1KTJQzYmNWrADTfAP/6R+nu/+aZlWQ2ncmVTf/3wQ2jf11+H1GM+\n4R5FEyaEjl9zDXz4oamN7rkHbrzR1EeXX27HzzzTjO4+CxaY99KRR8bu5/HHm3fh6jSGx06ZEry+\nRGkNyDFtBiJSP+xtATA8/JimpuZA4CpnmSQy4MynZk37EcQ67ii7hEciV6lShTZt2vDKK6+krT1V\njaJpLj1nnGG6bD9+YP580yFnIy/PrbfCMcfAn/9sOvRUsGsXvPeeubRG0r27DfC9e9sA/NxzZkMJ\nxxcGqiY4/v1v29+woQVqdehg8RJTpxZ14T3jjFAW0MqVQ9Hb8b7XKlUgN9e8vK68stQfPSqTJwd3\nX23Z0lJyFBaajSFZ4hmQJ2HGMMHS8K73XtcFlgJtkm+uKBoWjayqn4rI8/EETaaCzmKtDMCk79Kl\nThiUN2bPns2BBx64731eXh6jRo3i88jRpIxz+eU2qP31r/aDnzs3tkEx3bRoYUbtoUOLB0GVlM8+\ns9VPy5bFj3XrZqmqweIxrruuuDH90EPNtfbnn839MjytxkMPmSC58srig3zz5ra6mjTJZvyjRwcb\n4M84wwRHKoTBuHEm8MLdRqdMCf7dHnigrR5XrSoawR6YRHok4GXg7LD3vYEXg+igvPNbA9NjHGsc\n9vp4YHGc+wRTnKWAhx9WffDB6Mf69FF9552MdcWRIjpHKqEj9pFGA3KiLdlnu1Mn1S+/tNd/+pPq\nQw8ldXlK+fZb1bZtVbdvT839+vZV/cc/oh/Lzzd3008+sTa3bi1+zqRJ9v28/rrqxRcn1/bdd6v+\n5S9mk6ldW3XNmsTXzJtntof164sf891ig7Bzp9mD/vvf0L7CQtUGDZKzBx1/vNkywgn6bAdZTHRX\nCwzzhcenQKASEiIyDPgWOExElorIteEBZ1iVsxme2+ozlJEqZ0FWBo7yQX5+PpMmTWL79u1MmTKF\nyZMnM3nyZPLy8ti2bVup7y8i/jO8V0S6RBx7wKvTMVtEUuaA2K+fBUFB0YCzbNC9u83Y+/Urff6i\nbdvgo49ip2Fu3Bjq1LFI7RdeiB5d7aek+P774vaERPzqV7Yi+O47+079gLt45OTYau2yy4pmjH3u\nOVvdrF8frO2vv7bqiiPCql0sWxaKLA5K69YltxsEiTNYKSL/B3iPH32BlUFurqpxF09aRqucrVkD\nPXtGP+YbkR3lg9GjR/Pqq6+yfPly7rnnnn37a9WqxV//+tdUNDEduBCLnt+HiLQHLgXaY/m2PheR\nHG+mViquuMJ83J99NvvCQMTURGedZYbZZ54puf3i449NRdM4TgJ7P3NpLN/+OnVMSIwYERKYQTnp\nJLMlvP22CYag/O1vcPbZZnN45hnLOjtokMU9zJkDJ5yQ+B6ffGJxEMOGhdJuTJli9oJkvs9SGZET\nLR2A+sBgYIq3DQbqB1l2pHIjg2qi3NyQe1okb7+teuGFGeuKI0W8k0C3RynVRMA4oEvY+/uB+8Le\nf4oFcKbk2T7tNHsWDzpIdcOGpC9POevXW7TuoEElv8eFFxZ174zVzo4d8c/p3l21cuXoaqREnHGG\natWqquPHJ3fdunXmdt63r8U6zJ9vrqz//new6484QvWHH1RPOMFcRFVN/de/f3L9eO451VtuKbov\n6LMdpAbyOqBE5qFEEcjeOc9idoitwDXqFcPJJtHyEvm4lUH54vXXX6dfv34sXryYp58unjMufLWQ\nYpoB34W99+t0pIR+/Wz2WbOmzYazTd26ltyuY0fo2zd5V9eCAqvF8OqridtJRLt2lrU12SR9YCuO\n778P5TsKSr16ljfo2mvNw6ldO2jf3hLpJWLhQouW7trV8j6NGGGG6cmTzSU2GVq1shVWSUgoDERk\nHOZVVARVDRKIPhT4B14lsyj3LnMlLxcsMPesWNb4VOQNd2SOrVu3ArBlS8nLaMRLO6GqI0vTP59k\nPeUuushcO487LhWtp4bmzS319bffWu7/ZHj9dUvUlwoX1cMPL3lW1UsvNS+tqlWTv/aww+Cbb0Lv\n27e36nOJ+OQT83KqVMmEQe/epgKcMsX+JkOrVjBrVh4DB+YldyEEUhN1Ddt6YDlWngyy7PCujxeB\n/C/gsrD3swnzMIo4N7n1UgnYulW1Y8foCbF89u5VPfDAki1BHWUX0q8mGkUK1USqqpdeqnr99SX7\nvOni0UdDSfOCUlio2qFDyEOqtGzdaqkqss28eapt2iQ+r3dv1f/9z14XFqoeeqipiurWTT4Z4MaN\npjoMvy7osx1ETTQpYtc3IvJD1JOTp8yUvFSFm26yIJpb4yTRrlTJ/KuXLjV/aEf5YM2aNbz88sss\nXryYPWFuH0OGDEllM+GmvhHAGyLyd+yZbgek6ncDWI7+XbtSecfSc8IJsZPcxWLCBCugc9JJqelD\nSdRD6aBNG/P537Ytdp+2bTNPomHD7L2IrQ4eftgij5M1xteubauadeui18WIRxA1UXgkciVshZAV\nLWU6g87++U/L9Pjtt4n/Ab7doLwJg61bzdPhrruy3ZPMsWIFDBkCI0b04bTTTqJXr15UrlyZGTNm\nMHPmzFIXTBKRCzBVaEPgIxGZqqq9VXWWiPwPmAXsBm71ZmkpIycnlXdLDccdBz/9ZNXXgtYIHjLE\nAsiyEUWdTqpUMdvB3LmxU0qMG2ceQ+G2kPPOs6ytv/99ydr1VdnJCoMgy99FWIrpRcB8YAzQM8iy\nQ5NXE80hC2qiNWtsSRY06dZ116m+9FLaupM2hg1TrVSpbHifZIp77rEgpMqVj9Y//MHy5EeDchR0\nVtbp3Ll4htBYbNligWQrV6a3T9ni0kutFkEsbr3V6lGEs2uXjUfxrovH+eervvde6H3QZztI0Fl7\nVW2rqm1UNUdVzwSSyeovxK74NAK4GiBRyctofPedhcM//njw2qfRkkq9/rqF+ActF1dePYree8/U\nXOFGrv2ZrVvhtdcsQdmtt57LjBmfpD33vMNiAb77LvF5YD79PXpkPtFepmjfHmbNin5s717zHDr3\n3KL7q1a1FfxZZ5WszTZtYPr05K8LIgy+jbIv0L86UQSylqLk5caNlg/kV7+ygblnT/tS4y3Ef/7Z\nvITC/zmqlszqhhuCtlw+y19u22Y5VG65Bb78Mtu9yQzDhtlA07o1DB06mNGjz+Wnn6pTq1ZtatWq\nRe1UZVdzFOGEE0zdGoShQ6PXdd5fiOdeOnasubB3iJKr+de/tsI6JeGGGywCel2yqURjLRmAJph9\nYDbQGejibbnAnCDLjlRuRCyl+/UrGlyxe7cFfcTzSBg82IrWXHJJaN9331ntgmSs9l98oXrKKcHP\nzwSLF6u++Wbs4++/b4FKY8eqduuWmT7t3ZuZdqJRWGjqIT+Ax+fcc0OeG+FQQjUR8GtgBrCXot5E\nrYBtwGRvez7OPVL98bPKwoWqTZsm/k0tXWq5dyLrM+xP/PSTavv20Y9deaXqs8+mp92bb7ZcS6rB\nn+14D/lvMHe5zd5ffxsBXBTk5t59zsJsAfMIc7ULO34KsCHsR/N/Me6z74O+8Ybq4YcXd+/85z/j\nRwefeabqf/5jVZSmTLF911+v+vjjyXzNFl3YurW93rvX2n3sMatENWNG5gfB7dtVu3RRrVHDIiGj\n0a+f9XPbNnM927w59f3YvNkqLR15pP3Iq1TJXlK/8ePtGfH/F5MmTdJJkybpXXdN0ksusdcLFizQ\n3bt3q2qphMHhQA4wNoowiFndTGM82/sDhYWqjRuHKo7F4sknrYbx/sz27eaKvmtX0f3r19vENEgF\ns5KQn2+/wfnzUyAMNPSgXhzkRjGurQQs8H4YVbHi4EdEnBOo9CWge/aofvihVViaNKn4F7B5s30B\nP/8c/VjNmuaHO3iw6nnnmS9y3bpWRDwZduxQrVbN7nnJJTbT/sMfzHDTurVVX5o9O7l7loYbbzRD\n1RVXqP7978WP79xpRjo/+2HPnqqjR4eO79plaQQiH9hkWL1a9dhjTRj89JO9/+QTC7Pfs6fk9y0p\nl15adNbVrVs3rVq1qh5xRBc98MAuWrVqVe3cubO2bdtWR48enY44g1bEyNYb5do0fxuZ58ILQ1XE\nYtG5c+y0L/sT7dqpzppVdN+//qX661+nt92//lX1ootSIAyAft7f3wP3RG6Bbm7RxJ+GvS8SiKMh\nYZCw9CWgLVvawDtiROwv4I9/DC2PwvngA1OTqJq0btHCBtE+fZL8hj2aNLFAmSuuKJq+d+9eC1pr\n0ED1z39O/xJ46FAbcDdtUv3qK6sTG7kyGTVK9cQTQ+8ffLBozpOhQ83LKNr3Fo43iS7GvHkWKPPQ\nQ0VVA4WFlmtl2LBkPlHJ2bXL+vLuuyb8Nm4MHbvwwgt1xowZumePzcjGj5+pF198sf7888969NFH\np0sYbPZWu+OI44G3PwqDJ59UveOO2MfnzFE95JDsTBQyzXnn2TMZTrduqh99lN52t21TbdkyuDCI\nF2dwkPe3ZpxzEhEZVLYcq1sQyQleecAVwL2qGtX+/t57lr8jHrfdZj69AwcWDW3/+GM45xx7feCB\n8Kc/WZDZyBImEzjxRAtQ+7//K+ofXamS9eH88+H66y0L4ltvmc9xqpk6Fe69F/LyLPy+Rw/7bGPH\nmpeVz7vvWvoCn1NOsaAWMI+Gxx6D//3Psi6eeKIZr8KZO9fy4AwbZvlnTjkldGzjRjj9dPsebrqp\n6HUi1s6dd1qYf+XKKf34RXj7bbj6amjSxOI/nnmm6P9/3rx5HOnVMOzZE/LzOzBnzhzatm2b8N4l\nTEexEqvit95Lbf2BiHTQsIJO4WSqcFOmOPFEi2dRtfidb76xbKt+HqXhw9P/TJQVIo3Is2eb00sy\nmVGTIS8vj7y8PMCM+YE9H4NIjJJuwMXAS2Hv+wHPRpxTE6jhve4NzItxr8AS8dJLVZ95JvS+sNAy\nCc6ZE9q3a5cVsYk1200FO3ZYFsTrr08+rDwRBQWmkoo0Gr/wgi0NffbsseIbCxeG9m3ZYnaDrVvN\nmNq9u/Xvxx9VDz7YinBv2WLH+vSxfQMH2gqiXTubcfjceKMZq2JRWGirkpL6TAdh504L+8/Li33O\npSEqRsMAACAASURBVJdeqrfccovm5eXpLbfk6VFH/VYvueQS3bFjhx577LEpXxkkczyZZ7u8sH27\navXq9oy2bm2ZgE85xfYXFpqzx4QJ2e5lZnj1Vctm6vPHP9qWKYI+20Ee8oOB/sBLwBB/C3RzUxON\nCntfTE0U5ZpFREmRDeiAAQP2bePGjYv54b/7zgb/JUvs/ZQppsZI9YAchM2bbUmYyn/+rl3247r/\n/uLHNm0yFcmyZSYAzjxTtVev4ud17676+eeqRx+tOnJkaP+LL5rxr3Ztu/all0ww+Fx2meq999rr\nzz83dVu4OiYan31mxtxUqAQ2bzZBGM4LL5jQjce2bdv0qaee0gsuuEBzcy/QevVu0f79++tDDz2k\nDzzwQKqEQdew9w2BSt7rttgKuW6Ma0v/xZRBPv5Yddo0+93t2WP2tYsuslTN2fo9ZoMJE8w+sm2b\n6t/+Zr/PTNoUUykMvgWewAp1XOxvgW4OlQkZkKthBuT2EecEKn2Z7A/m6aetNN6SJaqPPKJ6553J\nf4mp4pdfzMOmTx8z3JbW2+iOOyy5VazB9bbbVE8+2ewWjz0W3TB8330mEDp1Kq7nHz8+tpfD6tW2\n0hg3zmbjH3+cuL+FhTYr7NVL9dNP7fPv3m22jKuvLiqM4uEL1pYtQ54q27aZ4E9mlrl7t9kNwqOR\nSyoMgAu8gX47sArPRgZchLmcTgZ+JKx0bJR7BO98OWbHDtXTT7fJRqyysvsjGzeqHnCAavPmNgZM\nm5bZ9lMpDKYGuVGc688C5mKpLO739t0M3OS9vs370UzxBE/KMjv+/e82YB15ZFHvmWywaZPNYI85\nxvoUzwgei127VB94wJbY0Wqu+sybZ7OwefNin/PJJ/bfjxebEIvXX7cCIFddFfya7dtVhwyxlUhO\njhngjzvOavg2aVJ8th/Jjh0mTK6/3v6v7dpZCoOnnw7mBDBv3jy9+OKLtX379tqmTRutXr2NNmoU\nSilZ2pVBabaKIgxU7Xdw0UXm8liRePhh01hkg1QKg0fizWoytZX0BzN4sC3LElVHyhSFhTarbtQo\nOW+C+fOt2PVZZyXvChuNzZttBVES1U1hoQ3Ia9eW7Nrvviu6TP7jH4sGAkayZ4+54V18cai/jz5q\n3lyNGwebafXo0UM///xz7dixoy5evFjPOGOAHnvsn/Ydd8LAsb+SSmGwGSj0lsGbvPebgtw8lVtp\nfjDhBs+ywvffm2E2MkI2km3bVJ96ymIrBg/eP/Ws27ebe2y0yODx403l1atXcYE+cKDqb38brI0u\nXbqoqupRRx2lqqpTp6pWrtxFH3hAdeZMJwwc+y9Bn+0g9QxKWDPIEJGzgGewALRXVPWJKOektfRl\n9eqluz4vLy/lrn7duoVcPq+6ylL+btkCBx0ERx4JRx0F776bx7vv5tKtm+UTipbDJF2k4zPH4sAD\nLUfNhRfCihV5tGuXy5o15sr688/w0ENW5jHSPXfAgOBtHHDAARQWFpKTk8Nzzz1Hs2bNaN58C7t3\nW4nBkiIiTwLnATuBn4FrVXWTd+wB4DpgD/A7VR1T8pZSTyb/x2Wl7Yr4mYOSMFGdiHSJsh0qIkFq\nIVQCngN+BRwJXCEiR0Scs6/0JWZL+FeJPkka8X12U81JJ8FHH1myqvbt4bTTrGTfTz9B//4wZkwe\nI0dafEUmBQGk7zPHont380sfNCiPF16w2InLLrMYh2uuKX2cxuDBg9m2bRvPPvsskyZN4r///S9v\nvvkagwaVOgPtGOBIVT0Gs4s9ACAiHTCni/bYROd5kbKVsT/T/+Oy0HZF/MxBCfITex5LUOcnRe2I\nGXzriMhvE8x2jgfmq+oSABF5E+iD5Sry6YNXI1lVJ4hIHRFprEmksi7PdOtmWzQGDrTCFxWF++6D\n7dvtc6ea47xiwTVr1mTo0KFFjpUm8ElVPw97+z3mbQdwPvCmqu4BFovIfOz3MKHkrTkc6SOIMFgJ\nXK+qM2HfjOfPwB+B97CZUSyCRCCXmdKXjv2P888/P+7xESNGpLK564Dh3utmFE317j/XDkfZJJFR\nAZgRax8J3E4JFoE8Ejgx7P3nRInUxML/3ea2tG1xnuPPgGlh23Tv73lh5zwIvBv2/h/AlWHv/02M\nbL/Z/txu2/+3lBiQgZki8gLwpvf+MmCWiByA1XaNxwqgZdj75t6+yHNaJDgHVS1T+lZH+UBEKgNn\nAFcAnYCPgeH+SjcIqhrXxCwi1wBnA6eF7Q70XHv3d8+2I+sEqXR2DRZFfJe3LfT27QZOTXDtRKCd\niLQSkWrA5Vg9hHBKVfrS4YiHqu5V1VGq+hssPcoCIE9Ebk/F/T1vuXuB81V1Z9ihEcDlIlJNRNoA\n7YAfUtGmw5EOxFumpq8B+7EMJuRa+riI3IwtXV7yznkOi1TeirnmTU5rpxwVCm8Vew62OmiNDdRD\nVDXqTD3Je8/HUq2s9XZ9r6q3esceAK7HJk5lzrXU4QgnoTAQkRzgMaADcKC/X1UT5/51OLKMiPwH\nOAr4BPPumZHlLjkcZZIgaqKhwAtY4MypmBvo6+nsVCQicpaIzBGReSJyXwbbfUVEVovItEy16bXb\nXETGishMEZkuIndmqN0DRGSCiEzx2h2QiXbD2q8kIpNFJJUuPv2wspS/A74VkU3etllE/OCwxSLy\nk/e5M6bKcc91Zp5rr+398dkO0m7wZzuAN9Ek7+/0yH2Z2AhQOjONbfcEjiFgLdsUttsEOMZ7XRNL\n9Jepz+zXlqiM+c0fn8HPfTc20UhYBjXF7S4E6mW4TfdcZ/C59tp0z3acLcjKYKcXSTxfRG4XkQsp\nXfWzZNkXuKaquzGvpj6ZaFhVvwbWZ6KtiHbz1UvJoVYZazYZ8lFX1W3eywOwOJT0GpU8RKQ55pHz\n70y0F9k8wVbJqcQ91xl8rr023bMdhyAn/Q6oAdwJdAWuAn5T4q4lT7TAtQoTvCMirbFZXEYiV73l\n7BQgH/hMVSdmol3g75hXTkZ+oBEo8JmITBSRGzPUpnuuM/hce226ZzsOCYWBqk5U1S2qulxVr1XV\ni1T1+5R11SNbOrWyjIjUBN7BPFGi1s5NNapaqKqdMb/4bl7EeVoRkXOA1d6sUbwtHe3UE5EjRaSt\nt9r16aGqXbDZ220i0jPJ+xbTwYvIABFZ7j3Tkz2vOgfZea5h/3624xD42Y4pDERkhLd9JSKLRGSj\niOwUkW3eg3+biNRJYad/B8yKsj9I4Np+h1giwHeA/6rqh5luXy3z5jjM5Tfd9ADOF5GFWDqHUz0v\noFLj5brqLyLTMT3xi8D/gCUi8raInKqqqwBUdQ3wPsVTpiRiKJaMMZKnVbWLt42KOOae6yw817D/\nPNtBSObZjrcyOAE4BdOvvQRci2VfvAgLvT8Q+FBE4id/CUACnVqQwLV0kg1pDlZrepaqDs5UgyLS\n0BfwIlIdi9ydE/+q0qOq/VW1pZq78uXAWFW9OkW3fwdTx5ykqoerak9VPVZVWwCPAxeLiB8XcBBw\nJpaIMZn+x9LBx3tu3HOdQfbTZzsuIlLDW4UFerbjpaNoguUWOgdLxRstjP9vItKw1L0O6dSKrTRU\nda9YtOgYQoFrs1PQZkJEZBiQCzQQkaXAAFUdmoF2ewB9gemejlOB/lFml6nmEOA1T4VSCXhLVT9J\nc5tpReOkklDVSSKyDnjf06dWAd7Q1AWH3S4iV2E1kH+vqhvD2nbPdeaea9gPn+0ANMaebSXAsx0o\nAllE6mGFv5/EUlpPwgp/J8pNFOTe5wC9VfV2EcnFfjTnRTkvG8YXRwVCS5EjSERaASNVtZP3/mDg\nF1VVEXkEOERVr49xrXu2HWklyLMd14DsBWpchPmq3o7pRq/DPIpeTUEfIQmdWib9c8O3AQMGVKh2\n98fP3LFjRwoLC5k6dSrHHHMMzz33HCeffPK+46lGVddo6MYvA8clOH+/+a7LctsV8TMHJZ4B+T9Y\nPvYuwCpV7YrpXv+hqpdglctKjWZRp+aoOFSpUgUR4cMPP+T222/ntttuY/PmzalsoogOXkSahB27\niCTtEA5HpolnM+iHJY77HVBDRLZixuSdIvJ/wJIM9M/hSAm1atXiscce4/XXX2f8+PEUFhaye3ep\ntZxAdB08tsI9BigEFmMlXdPGjz9Cfj6ce246W3Hsz8QUBqq6b9UgIicDfwC+UdUnRKQtls46pajq\nl8CXqb5vaclWEetsFs/e3z7zW2+9xbBhw3jllVdo0qQJS5cu5d57703JvVX1yii7026Q9dm6FS69\nFNavh4kToV27YNe556vitB2EtKewDtQJSzE8HksFXAV4R1UfjjhHy0JfHfsnIoJmqchMaZ/tu+6C\ntWvh2GPhrbdg/HioEqRslaNCEPTZDupNdBi2MmhN2GpCVU+LdU2yiEgNVd0mVpnqG+BOVf0h7LgT\nBo4S895773HfffdRUFCwz7AmImzatAkov8Lg22/h17+G6dOhXj0480w47TTo3z/FnXSUW1ItDH4C\n/oW5lO7196vqpNJ0MkZbNbBVwm81LHeIEwaO0tCuXTtGjhxJ+/btix3bvh1q1Ci5MBCRV4BzsbQD\nvmtpPeAtLCvpYuBSDYsziLi+RM/2jh1wzDHwyCMmEACWLYOuXWH0aOjcuSSfxrG/EVQYBM3UuEdV\nX1DVH1R1kr+Vso9FkOwlkXJUABo3bhxVEACsXFnq20dLR3E/8LmqHg6MBR4odSthTJwIZ58NRx0V\nEgQALVrAoEFw001QWJjKFh37O0GFwUgRuVVEDhGR+v6Wyo5oFpJIOSoOxx57LJdddhnDhw/nvffe\n27cBeJqiEqPR01H0AV7zXr+GBW2WmkWL4KKL4IIL4JJLYNiw4udcdRWIwPDhqWjRUVEIambyU1aH\nu18okPLSl6q6SUT8JFJFEtcNHDhw3+vc3Nwyb513lB02bdpEjRo1GDPGovHz8/PJz89n2rRpLEmP\nk3QjVV0NlsdfRBqV9oa7d5sQ6NMH3ngDqlePfl6lSvD009C3rwmOWOc5HOGUFW+ihsBuVd3oJZEa\nDTyuYblDnM3AWLcOli41XbEjNXz0EZx3XukMyFHSUaxT1fphx9eqaoMY1+qAAaEqjLEmOo8+Ct98\nAx9/bDP/RFx8sXkYPZBSBZWjrJOXl0deXt6+9w8//HDqDMgAInIU0AHLVgqAqqYqzXBHbCkdnkTq\n0YhznDAAhgyBV16xQSEIe/fCjBlw9NHp69O8eebbXinTtcKSYMeOHbzyyivMnDmTHTt27Ns/ZMgQ\nhg2Dvn1TLgxmA7mqutqLRh6nqlGNFkGe7Vmz4JRT+P/2zjw8iirrw79DIGELEAiQQAiRzQCCgEoQ\nVFBBcWNRUBmXoKDDiMrIIriz+Imi4oCOC4g648bIiICirOKooIBEZFEQZEtCSFhC2CHpnO+P02VX\nV1dVV3VXdxKo93n6Sae7+t7b3bfvuWe9WLcOSE01vfRPtm8HunSR1zYIWy9xqag46kAmOTz6Ve/t\nSkjBurBLVysw80aWmu8dmLm9VhBEgz17otfXokXAmhCPXS8oEOfhyZPWrv/yS6B//9D6ssoNNwCr\no3ZeVWjcdddd2LdvHxYvXozu3bsjJycH8fHxAACHqlJoS0IvADDYez8TQMi1+z0eYOhQYMIE64IA\nEAF9113AsGESeeTiYobVvdwAAFcD2MfM9wC4EDrlpqMBM3D6tPNtpqcD27Y5264RH38soX+hUFAg\ntuOffrJ2/WefAXl58h7VFBZK5Mn/wsz39niA3buBHTvCaydcNm8G3jXJ+d2+fTsmTZqEGjVqIDMz\nEwsXLsRqrwQL14HsLUexCkArItpDRPdAzkroRURbIb+d50Nt/803JYls2DD7r332WaBKFeDyy6O7\n4XGpeFgVBieZuRRACRHVAlAAoIlTgyCiFCL6mog2E9FGInrY6NpvvwVuCihwHR5FRbLTXrjQ2XaN\nyM2VxTgU8vNF5f/+++DXlpQACxbIgl2kiXDfvVt2iwMGiPYQKnl5IpzKWhh88w3w+uvGz1epUgUA\nUKdOHWzatAlFRUUoKCgAEL5mwMx/YeZGzBznLbr4LjMXMnNPlgN1rmHmw6G0feYM8PzzwCuvhGaG\nq1EDmD0buP12oHNn+ZxcXPSwOr1+IqI6kFK86wBkQSqaOkUJgJHM3BZywtpwIkrXu3DzZmDXLgd7\nhuy2AXEkRoO9e8URHAoFBRJRYkUYfPcdcN55csvLCxzDxRcDn38O3HOPlDEIhZ075W9ZC4PsbPGN\neDz6z99///0oLCzEpEmT0KdPH7Rp0waPPvoogPA1g0jyySdAy5aSSBYqRMCoUcB770mEkVWT0ZEj\nkrPguurODSwJA2Z+gJkPM/ObkOPiMr3mIkdg5n0sB0aD5YDs3wA01rv2jz+kOqOTFBQA7duL3dvJ\nhWHxYt9iqSYczaCgQMIFV60KnlQ0d674C5KSAoVBXh7QqJE4GJcuBQYPDi1Jadcuab+shUFOjixy\n27frPz906FAkJCSge/fu2LFjBwoKCjDMa3cpr8KAWRbj0aOdaa93b4lC+9e/gl8LAG+9BTz6aHia\no0vFwTTPgIjSmXkLEXXSea4TM2c5PSAiSgPQAYCuS/KPP0StP35cVGCFoiKxGf89hFqqBQVAs2ay\nOC5Z4p/RGQ6vvgpcey3w0EO+x44dk8UnHM2gXTugfn3ZCbdvr39daan4C5YtAzZs0NcMGjWS++3b\nSyz6wYPSrh127ZJaON99Z/utOEp2ttTm2bgROP983+MffPAB7rzzTkydOlX3dSNHjnTKgawLEe0C\nUAQpZV3MzIYHkmtZtkxMfddd59x4nngCuPNOYMgQ82J2Z84A06YBzzwjoam9ewMxMaH3+/33QNu2\n8h25lE+CaQajvH9f1rm95PRgvIc3/xfACK+G4Mf48eOxatV4AOMxf/43fs+tXy921VDYv1/s8Dfe\n6KypKD8fAQlNSumDUDSD0lLgwAFZsC+7zNxU9NNPQHy8OMaTk82FAaB/jRV27ZKx5OeH79g/eFC+\ni1DIyRHBu2GD/+PHjx8HABw9etTvtnHjRixYsADjx49HVtb48AZuTikkxLSjHUEAAC+9JFqBlZwC\nq3TtKhFJs2ebX/fxxzJ3nnkGqFlTP9PZKuvXS8RZy5bAk0/K9+xSDnHqaLVwbxAtZRFEEOg9z6Wl\nzNWrM7dty/z99+zHhx8yV6vGITFxIvOTTzLv3Mlcvz6zxxNaO1qaNGEeOND/sRUrmJs1Y05Ott/e\n/v3MdevK/VmzmAcNMr523Djmxx+X+y+8wDxypP/zffowz5vn+//qq5kXLbI/pquuYl6yRN7T1q32\nX6/m4YflOxwzhrmgwPrrPB7muDjmf/2LuW9f+/126cIsP4WIzOudAOoFuSZgTOvXMzdqxHzqlP33\nE4wlS5jbtDGe56WlzBdc4JsP//sfc1paaGPxeJgzMphnzmTesYP5vvuYExKY+/dnnj6defPm0N+H\nizWszu1gZqKbgwiSuWFLIx/vAPiVmacZXbBvn+xSWrYM9Bvk5kpE0IkTQPXq9jouKJA209KAhg0l\njj8jI/C6wkLrai6zvmaQmyvFxbxVEWyPU0keuvxyQFWdI6DvuXOlZAEgu/716/2v2btXHlcIRzNI\nSxMz286dQKtW9ttQ2L8fePppMfmkp4tt28rJXQcOiBaUkSGx+GoeftgwMA0AMH369IiaiSBlW5YS\nkQfADGaeaeVFL74o5sW4OOcH1LOn/Ebmz9fPQVm8WLSRa66R/6+4AmjTRkJcR4yw19fMmWJeuvde\niYaaMQOYOBFYvhxYsUJCX6dNk2gnl7IlWG0iJYizAYCukOqLgCSerQLgiDAgom4A7gCw0Vu5lAE8\nzsyL1Nf98QfQvLk4LLXCICdH/h44YC8xB5BFtls3ua+YivSEQffuwAcfGNvp1RQVid1Vz0zUooUk\nnp08aa9ujBJWCkgbp09L7Lj2/W7bJkJRiUBJTg78vPTMRHYd8x6PfO6pqRKxFK4TubBQ/CHjxkn5\n5dmzrQmD7GwgJUU+k7w88Sl588lwkfdDWLlyJX799VfcdtttAIA5c+agTRuphRhhB3I3Zs4jovoQ\nofAbS2E7P9R1t1q27IFFi3qYhsqGA5H4Dp57TiLTtGYoxWmtfnzyZBEO997r+2yDUVAAPPWU+D7U\nYbFJSRLVdMcdkucyeLAIpUgIvnMRbTkKy1hRHwAsAZCs+j8ZwGIrr3XqBoDfe4/5jjuYJ0wQs46a\nW25hBpjXrbOvRvXowbx8udz//nvmDh30r4uPZ1640FqbW7eK6SQ2lvnkSd/jI0Ywv/wyc1ISc26u\nvXHOnu1vdrr5ZjGPafniC+brrvP9v2kTc3q67/+SEuYqVZjPnPE9NnWqmGnssHs3c+PGcn/yZObR\no+29XkvnzsyrVsn99euZW7e29rp585hvuknuX3SRrw01GRkZXFxc/Of/Z86c4YyMDGZmrl07cmYi\n9p/Dz0BCqE3NRPfey/z009bee6h4PMwtWjB/953/42vWMKekMJ8+HfiaO+4Qk6oVSkvleitz4oYb\nmF95xVq7LvaxOret5hk0YWa1ESEfgM39tzFENIuI8olog9l1as1Aa9LIzZXoogMH7PevNr906SKm\nD60j8+hRueXnW2szP192240by85VQdmRJyTYjyhSjxMQx61eFE8w53BBAVC3rmSmKuh9psFQTESA\nz0wUDmozXJs2olUdCwgjCETRDADR2rROZGm78M9TzQDg2LFjKCwsBLO1PkKBiKp7gyJARDUAXANg\nk9lrdu0C5s2zb46xS6VKEnmnDbp48kk5JS02NvA148eLSceKA3jKFDFNqurvGfL886J5HA4pLc/F\nKawKg+VEtJiIBhPRYAALASxzcBx6h4MEYGYmys0VE0O4wiAmRkweRlFAdoRBgwZA06b+beXmioCo\nW9d+RJFWGHToIEXItGiFQUKCmKSUekba54HQfAZaYRCumejQIflcABFUbdsCv/wS/HU5OVJaAxBh\nsHFj4DXjxo1Dx44dMXjwYGRmZqJTp054/PHHceKE/sLnEA0BfO81ff4IKWRn6i2aMkUOpqnr6Gkh\n+mRmiplG+d6+/VZMjEOG6F/fooVUQp0yxbzdd98V/8LixeLjC8YFF4g58IUX7I3fxVmsJp09COAt\nSE2iCyGOsIfMX2Ud1j8cJAAjYVBaKv+3b29fGJSUyI6knqq4cJMm/rt5QBZxwJ4waNgwUBiEoxmo\nfQbKOBVfiRrtYk/kv/OPhDA47zz5fjjEbNXSUvke1A76Tp2ALAuZLGrNoF07fc3gnnvuwerVq9G/\nf3/cfPPN+OGHH5CZmYkjR4BatUIbczCYeSdL8cWOzNyOmU3rE+Xmip/kkUciMx4tNWvKwj99unxv\nTzwhu38z4fjUU+IUNjod7vPPRbNYtEg2PVaZMEGcy+Fql06iBIGcK1iudsLMc5n5Ee/ts0gOyghF\nGGidnQUFQO3aMvnsCoODB2UBUifU6AmDvXslSUdvcgweHFg4Ti0MlAJhzL6FOFTNoGFD3/+NG8sC\nol2Agy32Rs/bdSCrhYGykw01s/roUXGmq01XF10kJZuDodUMNmzQF0pxcXFITk5GQkICfv/9d3z7\n7bc4ejRywsAuL7wgcyma5aYffBD497+l7MXBg+LUNSMlRcqX/J9OXeFDh2T88+f7J/5ZISVFoox6\n9gzUysuKWbMkqu1cyYuwdNIZEXWBlK9uDSAWQAyA48wc1Z9RUdF4vPGGRLHs29cDpaU9UKmSLIgp\nKUBiomTl2kFregGkLe2Oe+9eMVvoCYMff5RQz4sv9j2Wny9mnNhYn13/4EHxa1Sr5ozPoFo12d0p\niWjqsdoVBvHxsoCqI3GCsXOnb/Eg8vkN9Ewcw4dLxEjPnvptFRYGvq5TJ8niDoZaM6hfXz4XtYAA\ngLfffhvTpk1DTk4OOnTogFWrVqFRo0bo1Ssz0qGllvj1V0n02rw5uv02aSLZxXffLaHIVrKMx42T\nRXL0aNEIFWbMEHNPZ1updT6GD5eih1deKWGnTZuG1o4TnDghWlKnTqK1TJ9edmOJFlY1g9cADAKw\nDUA1AEMB/DNSgzIiPX08JkwYj2efHY/atXv8KbFzcmSXnJhoXzMoKAgswWBkJurUSV8Y5OYGFs/T\nMxMp/gLAGZ8BIIugnhZjVxhoTUlWUGsGgHl46bp1UgPJiEOHAnM4LrhAag2Znd1QWurbDCjomYqm\nTZuGtWvXomnTplixYgV++eUXdOzYEbffPh7p6eONOwgTIupNRFuI6HciGmt03YgR4rwti0Noxo6V\nRfxm06wiH/XrAw88AEya5HvszBkR3OGauP7+d/ksevTwmWbLgunTgUsvFbPdxx8DW7aU3ViihR0z\n0XYAMczsYeZ3IWcUO4n2cJAAmjf33Vf7DZRFNpgwKC0F3njD/zGlFIUaI82gY8dAYXDkiESjmAkD\nxUykXoQTEuwLA6VNNVq/QUmJaCDa96Q2A+kJA+Uaq8KgpETaUe++zZzIeXnmB+ConccKcXGyA9Xz\nASjs3y9mnqpVfY/pRRRVrVoV+/ZVRU4OcPr0aaSnp2Pr1q22NCG7EFElyEbqWgBtAQwyqsa7d68s\nsGXBhRcCn35qr0T2yJFSHv333+X/OXPENOTEcawjRkh9MLWwiSYHD0opkP/7PxF8Y8dKwb6zHatf\n/wkiigWwnoimENEjNl4bFIPDQQJQCwP14mZVGOzaJT84dQib3m7bSDNo315eW1Li/7jSthpl4VYW\na2UHq9YM7JiJTpwQFVq7cGkFV36+fA7aImTqhV6pWKrFjjDIzZXPTe1sNAovZZbvat06/89OjVF2\ndzAncna2v0AC9IVBSkoKVq48jJKSfujVqxf69u2Lpk2bRtSBDKAzgG3MvJuZiwHMBtBX78Jp0/z9\nJeWdhATZxU+YIN/vK6+IgHCKRx8VP0YoWfHhMnkyMHCgL5v+oYfEfLd8efTHEk2sLuh3ea99r2xV\nWQAAH+lJREFUEMBxyME2tzg1CNY5HETvOq1moEwUq8Jg61b5+9tvvsf0hEHjxtK2uqTz3r2SaVu3\nrn8OQm6ujMtIGFStKj+cvLzwNANFg9Fmi2qFgdGuP1g0EWDPiaw1EQHGZqLCQrHjp6QY28T1NANA\nhIGZEzknx99EBEhpEe04PvvsMxQW1kFx8XhMmjQJQ4YMwbx58yLtQG4MQL2tyIFBaXYjX0p5ZsQI\nMf299Zb4mq6/3rm269eX6qqhFp8MxsqVvnItarKzJTT26ad9j8XFSTjtyJHG52WcDQR1IBNRDIDn\nmPkOAKcATAjyEtsQUW8A/4AInFnMrBtxHMxMVK+eCANm/UqPamFw6aVyv6Ag8OCQuDigTh1f4lhp\nqSykycmywCuPA7IYde4savaZM7JTPn5cxqDEWCumor17faUs7GoGeiYiQBbCZaqMj2AmoOJi6Vev\nVLUdzUBPGBiZifbtk+8rI0Oc7RdeGHiNngMZkO/m7beNx6GnGaSm+h/x6PF40LZtW/Trt8VbpqP7\nn9/NkSORMxPZQV2OokePHujRo0eZjcUq8fHAmDGibb/2WmgnsZkxerSYZh97zNnS13PmiLO6pESq\nqdap43vu5Zel5Ia6bhcg/pSpU0WA3H23c2OJBKGWowgqDJjZQ0RNiSiWmc+EMjgzVHbVqwHsBbCW\niOYzc4DLRisMlFhnxYEYFyc7caPd3tatIjSCaQaAz1SUnCw2xPh4aVsRBgq5ubL4KJnGzZv7Fm5F\nIClO5NxcX216u5qB2TitaAbKQq/kKuhFjSQn+382Zuzc6R9JAsj7zM6W3ZO6fUWQZmSI3+Cvfw1s\nz0gzaN9enHenT+vXrtHTDJKSpD3lNTExMTj//POxZcseAKkoKICfMIigZpAL/0z9FO9jAaiFQUVi\n+HBJ8svMdL7t1FQ54vaf/xTnuhkvvihzLpip6pVXZMFfskT+vvmmREcB8jv/97/1kxaJRDv4y1+A\nW2/191GVN7SbiQna6o0GWJXlOwCsJKKniGikcrM9Sn0s21XVBdn0fAaAualo61agb19rwkBtfsnN\n9S2wesKgcWPZJSumIu0uPjVVhIF6obarGVgZJ2AsDBo0kP727NF/HrAXTaSnGVStKhqH1vmuaAZd\nuhg7kY00g2rVRMAahQzraQYxMfIe1eMoLCzEF1+0BXA1MjP7oE8fuUXYTLQWQAtlMwXgdgALItZb\nGVC9uiyg6oOmnGTsWIlSOnHC+JqiIl9JC+/xFbq8+66Ev65aJY7uMWMkakg5h+ONN6Rwn1GyXLdu\noqm89lro76c8Y1UY/AHgC+/18aqbE1i2q+rV0jl6VEwfiqoXTBj062dPMwBkgVUmSCjCQDETqYVW\nnTrijFb7JZRT0PTQJpwpNG4si56SZGUkDCpXFjPazz8HqsAKemYixbSmRU8YAPqmIkUzaNdOhGJR\nUeDr9EJLFS66yNiJrKcZAIGmokmTJqF+/S+Qmvo0evcehVGj5BZJMxEzeyB+tiUANgOYzcwWdS8X\nAGjdWhbhWbOMr3n9dfFXdOtmfKRnSYmUy54507epbN9ebh9+KOHLr70mAsKMyZMlOTDU5MryjKWk\nM2Z23E8QLorPQFlgFZNMYqL+aVlHj8oX2L27LJhK+WgnNIOUlEBhoG6zaVNg4UL/kM8qVWRXdfSo\nZE8DouquW6d/2lp+fuAOGPAlsR08KO9dPVYtycnSvtnzagdyUZEk2u3dG/gZ7d6tnxSk/hwUFM2g\ncmVxCK9dG+gwNTITAbKLM6pRpKcZAD5hwMwgIlx+eXccOiRRIg0ayDwAgH/8A4iPj9yJ7yxl2G3m\n47qoGT0auOsu8U1ozZsnT0ok1vLl8vsePFjMkNrrPvlE5v1ll/k/PmaMRAudPClmzNatzcfSurUk\nT06eHLxGU0XDVDMgoplE1M7guRpEdC8RBUlgD4otu6py2779Gz9hoGCkGfz+u0SZxMaK2WHrVpkA\np0/rmwmsagZKwlswzWDNGlmE1CGfWr/Bli0iNPTi6o2EFuAvuIw0A0AW+6ws4+cTE0UzOeP1DH3z\nje/MAjWlpf6fiZrzzgsUBopmAPicyFrMDg5KT9fXUMzGoQiDK6+8Eq+++iqysvagVi35LvLzgaVL\nlyIzMxPLlnXA1Kk6J7y4lBu6dpW5P29e4HPvvCNzqm1b0QwSEyX/QQ2zmJEefzzw9VddJZupRx+1\nnkswfrwENdgt31LeCWYm+ieAp4joNyKaQ0SvE9E7RPQdJC8gHnJmcThYtquqhUGfPj1w4oTUK1Kb\nCYyEwdatvnoprVuLqWj/frFx60UeGS2wamFQXCw7ciW5zMxnUFgYuAhr/QbbtkmyzfM65cyMzESA\nvxPZTBgkJUlop9HzlSrJj06Z5EpctTYTVC/RSyEtLTDXQNEMAGO/gZlmkJ6unwGq1KTSG4ciDBYt\nWoSYmBgMGTIIhw83wqxZbfDii+dh6NChiI2NRWrqu5g6VWeVcSlXjB4tiWBqiotFm37sMfmfSP+6\nhQtFU+itkyZLJMX3Lr/cd8BVMBo1EkdypMJeywpTYcDM65n5VgCXQATDd5CFeigzX8jM05g5rGPQ\nQ7WrEsnimJVlTTPYsiVQGJjtttWagZGZKC/Pt9s30wzq1BG7tHYHq9YMmKX0ghLp8Mcf/tdqTU9q\nFMF1+rTs7BMT9a9LThbbqZEwUK5R/AbLlolZRysMtOUf1FjRDFavDiwkZ+RABmRhP3gQATWEjPwF\nymv27JHM4wceeADPPLMSvXvvxtNPL8dVV/2M3bt3Y+bMmfB4OkbEgUxEzxBRDhFleW9OZ+yfU/Tr\nJ7/XVat8j739tsy3Ll18j/XvL/NNiaxklhPdxo3T3/QpbX/1lb3xPPqo9G+3vlh5xmoJ62PM/A0z\nf8zM85jZwK0YGsy8iJnPZ+aWwcr8qklKkmqhVoSBkWZgtMAqiWcej7GZSL0oNm4sk/XMmUBhQCSa\ng3YRVherKygQE1ZqKjBsWODuJpiZKDtbxpuUZBzvrSzIZsJAcczn5sr7uPHGQDORYhrTI5hmkJIi\n/hL1NcXFYrIzcuRWqiTZoErpAwUjfwEQ6EDOyQGaNq2C9PRkHD7sCyyPcDTRVGbu5L0tCn65ixEx\nMZLx/PLL8hsbNUrs9trdeUyMlJEYMEB+czfeKL/zAQPM2zcSFEakpkpk4tlUwM7hNJHokpQkMcF2\nhUGbNsE1g9hY2anm5/trBvXryy5VsaUrfVeuLNdkZ+sniDVtGriAqovVbdsmh4cAktn5n//4zDWl\npT4HsR6KZmBmIgKsCQPFifz111I9MjU1UDMw25GnpPiEIgCcOiVRUupdf0aGOJEVFH+B2Q9Sz1S0\nZYt/7okanwNZ/lcqmzZo4O/ziXDSmc0lxsWMe+6Rw3g6dxbN+eef9WshDRokAmDZMuC224D337dW\njdUu48ZJDkR5qHrrBGUuDIhoABFtIiIPEXWy89rkZFl0ggmD0lJZbBVh0KqVTCa9KBk1TZpImGRh\noX8UUO3asjhrndeKqUhPGNx4ozjC1Kg1g23bxMENiMC54w45bIRZrqld27h2jeIzsCIMqlTxP8hH\n75q8PPkh9ezpOzNBjfZ9q1GEorIrVz4LtbaSnu6/yzcLK1U4//xAYbB2LXDJJfrXx8dLwplS2VbR\nIho2FGEFyLw4ccLaaVwh8iARrSeit4modsR6OUeoWVP8afffD3z2mfk8JpLf0913+5uRnKRVK+Dq\nqwOLX1ZULIWWKhBRdWY2Sf8IiY0A+kNOUrOFYnoIJgxycsQUoJgDqlWTBevHH82dRikpsuBoM3YV\nU5GeMNiyRXbD6hR3QEw/WtSawfbtPmEASEx09+6i8vbvby60rGoGzZqJD8CsbIAScbR8uTjWTp7U\nNxOZVUtQhGKLFv7+AoXmzf3PbjbzFyikpwNz5/o/tnateXifoh0kJoowqFJlN7KytqGoqCeOHj2J\nI0dKUL16fMhlFIhoKeRoyz8fAsAAngDwOoCJzMxE9CyAqQAMDpSsmOUoyoKhQ8t6BP48/jjQq5eE\nLGsz8suKiJWjAAAi6grgbQA1AaQS0YUA/srMYRfdVfwPRHatdj77uCIUAH1hoDYRKbRuLWe+9jeJ\nKmzSRJyd2l2wWhio6+ykpflCSK28m4QEn6N42zb/sdSuLUcHdusmJi0zYaCUwjDLMQDkc9IL61ST\nnCzOt5gYWbQPH7anGQDyo1B8Amp/gUKLFpINqmBFM9Caifbtk2xTIzMR4BMGnToBW7bMxLPPzsDx\n44dQr94f2LAhB2PHDkOtWqGXomTmXhYvnQngc7MLKmo5inOd9u2lqF2vXrLBMUrojCaRLkfxCqQm\n+0EAYOZfAFxha4QRIClJFma1+aRuXVnA1NUFjYTB4cP6BdsUFGGgXWAbNoRujkNamlxvFAKqResz\nUGsGyvtbskR26WZtxseLj8MsbNQqSUkylp49RaDVqSMOXrVd1MxnAPhHVu3bp68ZqKOlzMJKFVq1\nEu1J+V4VE5GZ0FWEgccDHDr0T6xcuRK1atVCgwZAzZotUVBQEDHnMRGpReDNAGyewedSURg+XJLd\nrr22YmcmWzYTMXO2ZvNuuZirmTrNzKY7JjVaVbplyx5op0mJq1xZzEHqQ+6NhAEQ3Pyya5dUNlSj\naAbaqJq0NOnLailfxWeghJUqDmQ1zZvLEYDqMxiMxrpmDfDww9b6NkJZuK++Wv4S+fwG6eky1mDC\n4LzzfKF6SoSTmkaNJLv52DGxA1sxE9WoId/V7t1i7jLzFygowmDfPqBy5TjEx8vhCw0bAl99tRwH\nDhSgWrXxiNCmfAoRdQBQCmAXAJ3yfC5nC088Ib/lm24SzVp7nkhFwOqQs72mIiaiKgBGALBcY8WG\nOm2Kniq9eHHgdYqpSC0MtAknbdrI32AOZEBfM8jPD8x+VWr12NUM8vMlcUrrZ1AIliKvjNUpzSA2\nVjIzFVJSfMLgyBEREGY7aq1m0LGj//OVKvnOPmjf3pqZCPCZihRhoOeHUZOaKqHH2dlAvXrd8dxz\nz+HkyZPweJbiv/99Hb17D0VBgQgDq6q0VZi5nBc6dnESIgkH79VLHMoPPVTWI7KPVTPRMADDIQXk\ncgF08P7vNI6E4mn9BkaaAZG5mUjZ/WoX2KQksePHxflXa2zcWGztVs+xVTQDPRORXYzGapfYWFmk\n1bt5dURRMK0A8BcGepoBIFrQ9u1y34pmAPiEAbM9zSAnB7j00udRv359tGvXDtnZb+G8867HLbc8\nG8kcA5dzjEqVpMLqxIm+iLWKhNWkswPMfAczN2TmBsx8JzMfdGIARNSPiLIBdAHwBRHZzAUMRC0M\njh2Tnbe2wmbt2hI1U726cTuNGvnMJGqUzGftoli5si980QqKZuCUMIiLc+YQEO37VSqjAsGdx4B8\nbgcPSiSSns8A8PcbWNUMlPDSXbvkvQYTfIowyM4GUlMr4b777sOcOXMwaNAM1KyZgaNHKewcA7PQ\naCJ6jIi2ecu5XBNeTy4VgTZtJJxVKZFRkTA1ExHRqxDbvi7MHKaFGmDmeQAcLQ6jFgYffiiqm54N\nL9jh3bGxspBpF/2GDcVEdMEFga9JS7MuDOLjJc79t9/0/QV2SEnxCS+nSUnxRfJY0QxiYkQo7tlj\nrBk0b+47RMSOZjB7tjWtAJDv7sABETpffdUDkyYtQElJCd544yIADbBzZ1e0axd2gRnd0Ggiag3g\nVgCtIcUXlxFRS2ZtIQ6Xs41nnpG5unq1JFhWFIJpBj8BWGdyCxsimuLdOa0nok+JKGzFXREGpaVy\nVN2oUaG3tWKFfLFqlMVeb4f84ovWHciVKomfYO3a8DWD5s31zxdwAruaAeDzCeTnBzcT2fUZWBUG\nyiE3P/4IeDxFqFWrFubOnYteve7GJZesxpYty8M2EzHzVmbehkATZ19Ina0SZt4FYBvkICeXs5xa\nteTMg+HDfZn40aSgQDK1s7P9z0sJRrBCdf8yu4U7aC9LALRl5g6QH0zYCpYiDBYulGgVpXZ9KLRq\nFbjbVnwCeovixRcbO4L1SEgQJ2e4wqB7d2D+/PDaMEJxIAPWNANABNO6dfL56x1XqTUTWdEMkpMl\noW/xYmvCABBT0fr1QExMCfLy8vDJJ5/g2mtvREGB/FAiWIpCe2hTLgwObXI5+7jzTlkfhg4NLMoY\nSZiBPn2kjlNGhr0T6Cz5DIjocyJaoLm9T0QjiCis00CZeRkzK/LrR4hKHRbKATcvvyxagdOmk9hY\nWcStLIrBqFtXkqfCFQZEkVvY7DqQAREGP/xgnITTtKmY2s6csW4mIhLtYONGEbpWSE2VPIPRo5/G\ntddeixYtWuDKKy9Bbu4O1KzZ0pJmQERLiWiD6rbR+/cma6NwOdcgAj7+WPyBeuco2GHzZvF9WuGz\nz6R68bp18vsyOvVRD6uhpTsA1Afwsff/2wAcBdAKkl15l/UuTbkXcgZyWCQm+g5mGTgw/EHp0bCh\nNXNJMBISRNMoz1EtDRvKpCoutmcmmjrV2C9TpYoIlZ07rZuJABEGVoUHIEKHCBg6dCD+9jeZDKdO\nAYcONcOAAZ9a+txDDI3OBaCuqWp4aBPglqM4G6leHfj8c6kikJQEPPig/YJ5J05Iwchu3aQci9nG\ntqREHNdDhnyDiRO/sT1eq8KgKzOrFfPPiWgtM19CRJuDvdhK0hkRPQGgmJk/MmrH6g8mMVGSk6ZM\nMS7uFi4PPSTVE8Olbt3wnceRpnJlEQh5efY0g0OHzNPzW7SQ4yxjY/VNSXq0bi2T3iqpqTIGj+cU\nZsyYhc2bN+PUqVPwePIwf/4+eDx9//RdOID6p7oAwIdE9ArEPNQCwBqjF7rlKM5OEhPFrNm7t5yB\n0KiRWAFmzLDm43vnHTGJ7t0b3P/5zjvy2xwzpgeIevz5uOUcGmYOeoMkmKWq/k8F8Jv3/s9W2gjS\n/mAAKwHEmVzDVtmxgzk+nrmw0PJLyowHHmDOzCzrUQSnc2fmr79mjo1l9niCX793LzPAPGqU8TV/\n+xvz6NHMjRtbH0dhobRtlUWLmC+9lHnAgAH85JNPcrNmzfi9997jatV6cfXqD/PXX8t13vkVytzt\nB/ENnASQB+Ar1XOPAdju/f1cY9KG9TfkUmE5dYp52zbmiROZL7iAuajI/PrTp5lTU5lXr2betYu5\nQQPm777Tv/b4ceZGjZjXrAl8zurctjrhrwewB8AKAN8A2A3gBgA1APzdShsmbfeGnHBWL8h15p+c\nhv37bV1eZrzxBvPMmWU9iuD078/8/PPMaWnWrvd4mOPimF96yfial15i7t6duV07R4aoS0kJ8+7d\nzB06dGBm5nbezi699AwDGbx2rVwXqjBw4uYKg3OL0lLmYcOYr7uOubjY+Lp33mHu2dP3/8KFsnHK\nzw+8dvJk5oED9duxOretJp19CaAlgL9DSlGcz8wLmfk4M//Dmg5iyKuQaqhLvccDvh5mewCMD4Ip\nbwwbVv7K8uqRkqJfwdWISpVEDdYLK1Vo0UIcXU4kyhkREyOmoipee2GdOnWwadMmJCQUAYhcoToX\nFyOI5IS04mJjs4/HI2c3qJ3P118P3HKLZDirOX5cTnwLt6KKnXJKFwFI877mQiICM/87vO4BZg4z\njsYlGjRuDHz6qRwcbpU+ffxLfGtp3jzwFLRIcf/996OwsBCTJk1Cnz59kJd3DMBEVxi4lAlVqgBz\n5kj45yefALfe6v/83LlSW03rFn3iCQmiGDvWVzttxgzgiius1TAzxYr6AOB9AKsgB3a86r1Nt/Ja\nC21PBPALgJ8BLAKQZHCdsT7lEnHef1+MiiNHOtfm8ePS5r33OtemVZ5+Wvo+dkz+R+g+gwGQ8tQe\nAJ1UjzcFcAJAlvf2ukkb0XrbLuWMH39kbtjQ3/Szfz9zixZiFtJj9GjxtzEznzwpvoKsLOM+rM5t\nq4XqLgbQjZkfYOaHvLewS1F4mcLMFzJzRwALATzjULuOEcqpQRW5X72+lQgiJ3IrFKpXl+gKtZko\nUu/59OnT+Oijj/Dcc89h4sSJWLduIogmmtamsohSjuJ/Os9tZ+ZO3lvYB0E5TXmaX2d7v0Z9Z2QA\nmZmSrQxIKGmfPmIOMqpkMGaMlGXZswd47z0J39ZWBg4Fq8JgEwAT62/oMLM6naIGpP57ucKduD5f\ngZPCABC/gdpMFKn33LdvX8yfPx+VK1dGjRo1UK9eDcTF1Qg7IZGNy1HA4LFyQ3maX2d7v2Z9T5gA\nbNokSWqDBon59LnnjNtp0AC47z7xHbzwgpiOnMCqzyARwK9EtAbAaeVBZu7jxCC8Z8TeDeAwgCud\naNPFWRRh4ESinZp27ZwXMHrk5ORg0aJFf/7/7bdSdyrCpBFRFoAiAE8x8/cR79GlwlG1quzwL79c\nysrMmWN+Tjkg2kFamuQgdO3qzDisCoPx4XQSLOmMmZ8E8CQRjQXwULj9uThP9eqSeNa0qbPtTp8e\nmUqrWrp27YqNGzeinfdovPPPl2MKrRDiSX17Ibk5hd7S1vOIqI1GE3ZxASDmoi+/lETW2Njg1ycm\nAm++CbRt69wYSPwLNl9EdBmAQczs6AE3RNQEwJfM3E7nObf0r0tEYeaQxRIRrQAwipmz7D7vzm2X\nSGNlblsOLSWijgD+AmAggJ0APg19aH7ttmBmpSBAPxgcpxnOD9Xl3IWITHUZZt7tZHeqfhMBHGLm\nUiJqBilHscNgDO7cdilzgh1u0wrAIO/tAID/QLQJJ+36z3v7KYVkNgc52dbFxRb5kDnVAhL5M4uZ\nbVQ3MoeI+kFCrRMhJ/WtZ+brAFwBYCIRnYHM7b8y82Gn+nVxcRpTMxERlQL4DsAQZfdORDuYuVmU\nxufiEhZE9B8AxZB5fB2A3cw8omxH5eJS/ggWWnozpPjWCiKaSURXowzC5YioNxFtIaLfvU7maPU7\ni4jyiWhDtPr09ptCRF8T0WZv7XyncjqC9RtHRKuJ6Gdvv1HN+SCiSt6SJAscbLYNy5ndb0ESxAJy\nqIloFxH94n3fhpVFncad19GZ196+z8a5baVf63PbSmYaJP7/LwA+B3AcwBswqcLo5A0isLZDMjqr\nAFgPID1KfV8GoAOADdHoT9VvEoAO3vs1AWyN4nuu7v0bAzlsqHMU3/cjAD4AsMDBNrPM/vc+tgNA\nQpS/Y3deR3Fee/s8q+a2xX4tz22rheqOM/NHzHwT5JCOnwFEayfTGcA2Zt7NzMWQw2/6RqNjlrjw\nwmj0pel3HzOv994/BnGqR+XIRGY+4b0bB/EpRSXShYhSINVx33a46QuJ6Ij3dhRAe+U+ER1Ruof1\nBEyncOd1FOe1t8+zbW5b6h4W57btHwAzFzLzDGa+2vawQkN7lmwOzqGzZIkoDbKLWx2l/ioR0c8A\n9gFYysxro9EvgFcAjIHDP1BmjmHmWt5bPDNXVt1XytQxpGruWiK6z8n+TXDndRTntbfPs2puW8Ty\n3I72bsjFBkRUE8B/AYzgKCUrMXMpS52oFAAZRNQm0n0S0Q0A8r27RkL0/VLdmLkTZPc23JtH4xIh\nymJeA+7cRpC5XRGEQS7kZDUF07NkzxaIqDLkB/M+M8+Pdv/MfARymFHvKHTXDUAfItoBOWf7SiIK\nuzy6VZg5z/t3P4DPICacSOPO6zKY14A7t42urQjCYC2AFkTUlIhiAdwOOV82WpSFNAeAdwD8yszT\notUhESUSUW3v/WoAegHYEul+mflxZk5lCVm+HcDXzHx3pPsFACKq7t2pgohqALgGUpgx0rjzOoq4\nczv43C73woCZPQAeBLAEcjzmbGbWzVJ2GiL6CHKOQysi2kNE90Sp324A7gBwlTckLIuIorGLSYaE\nEa+H2HIXs5xydzbTEMD3XlvyjwA+Z+Ylke7UnddRndeAO7eDzu2QahO5uLi4uJxdlHvNwMXFxcUl\n8rjCwMXFxcXFFQYuLi4uLq4wcHFxcXGBKwxcXFxcXOAKAxcXFxcX2DjpzMU5iMgD4Bf4ztLtx8x7\nynZULi7h487tioubZ1AGENERVZE0vedjvElJLi4VCnduV1xcM1HZEFAGgIgyiWg+ES0HsMz72Ggi\nWkNE69WHcRDRE0S0lYi+JaKPiGik9/EVRNTJe78eEe303q9ERFO8h3usV6oXElF372vmENFvRPS+\nqo9LiGil9/ofiagmEf2PiNqrrvmOiNpF6DNyqZi4c7uC4pqJyoZqRJQF+eHsYOZbvI93BNCOmYuI\nqBeAlszcmYgIwAJvxcETAG4F0B5ALIAsAD8Z9KOofUMAHGbmDG8dnJVEpKSldwDQBlLWdyURdYXU\nzZkNYCAzZ3nrm5yE1GO/B8AjRNQSQBwzb3TmI3E5S3DndgXFFQZlwwlvWVktS5m5yHv/GgC9VD+s\nGgBaAqgF4DNmPg3gNFk7Ru8aAO2IaKD3/1retooBrFEqG3rrtqQBOAJgLzNnAX8eRAIi+i+Ap4ho\nNIB7Abxn6127nAu4c7uC4gqD8sVx1X0CMJmZZ6ovICKzw9xL4DP9VdW09RAzL9W01R3AadVDHvjm\nRIC6z8wniWgpgH4ABgK4yGQsLi5q3LldznF9BmWDldLBiwHcS1J6FkTUiIjqA/gWQD+SA77jAdyk\nes0uABd77w/UtPUASS15EFFLIqpu0vdWAElEdJH3+ppEpMyVWQCmQ3ZdRUYNuJyzuHO7guJqBmVD\n0BAuZl5KROkAfhCzKo4CuJOZfyaiTwBsAJAPYI3qZS8B+MTrRFuoevxtiIqc5bXRFkB2QLrjYuZi\nIroNwGsktd9PAOgJMQFkkZwd/K6dN+xyzuDO7QqKG1pawfFGYhxl5qlR6q8R5ICO9Gj053Lu4s7t\n6OKaiVwsQ0R3AfgBwONlPRYXFydx57arGbi4uLi4wNUMXFxcXFzgCgMXFxcXF7jCwMXFxcUFrjBw\ncXFxcYErDFxcXFxc4AoDFxcXFxcA/w/eYVKh9Ail2QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10e7312d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.subplot(3, 2, 1)\n",
    "plt.plot(t, ss)\n",
    "\n",
    "plt.subplot(3, 2, 3)\n",
    "plt.magnitude_spectrum(ss, Fs=Fs)\n",
    "\n",
    "plt.subplot(3, 2, 4)\n",
    "plt.magnitude_spectrum(ss, Fs=Fs, scale='dB')\n",
    "\n",
    "plt.subplot(3, 2, 5)\n",
    "plt.angle_spectrum(ss, Fs=Fs)\n",
    "\n",
    "plt.subplot(3, 2, 6)\n",
    "plt.phase_spectrum(ss, Fs=Fs)\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using matplotlib, we can showcase our generated data set in several different ways. Of course this is just an example and there are many more types of plots that matplotlib can generate. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# Basic Programming Techniques in Python\n",
    "\n",
    "To showcase some other basic programming techniques using the skills we just learned, lets analytically calculate $\\pi$. \n",
    "\n",
    "The easiest way of doing this is using the Monty Carlo Method through finding random points inside a circle inscribed in a square. \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np #always good to be explicit about our libraries used"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we will initialize our counting variables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "total = 0\n",
    "inside = 0"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we will set up a for-loop to calculate our points. Loops in Python are a bit different than in IDL. First, Notice that in Python loops are tab delimited. Unlike IDL there is no need for 'begin' and 'end' statements since the tabbing indicates where the loop starts and ends. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "for ii in range(10000):\n",
    "    \n",
    "    # Generate two coordinates in the interval 0-1:\n",
    "    x_coord = np.random.uniform()\n",
    "    y_coord = np.random.uniform()\n",
    "    \n",
    "    # Calculate the distance to the origin\n",
    "    r = np.sqrt(x_coord**2 + y_coord**2)\n",
    "    \n",
    "    # Count this if it is inside the circle.\n",
    "    if r< 1:\n",
    "        inside += 1\n",
    "    \n",
    "    total += 1\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You will see somethig else too: the variable we are iterating over, 'ii', is nowhere to be found in the loop. This is because the iteration of 'ii' over the list that the range( ) function produces is implicit. We can make the interation explicit using a slightly different syntax. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Pi=', 3.1336)\n"
     ]
    }
   ],
   "source": [
    "print('Pi=',4.0*inside/(1.0*total))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this case, 'iteration' now contains the iteration number (0,1,...,9999) while ii contains what ever is in the list at that iteration. This could be any data type. To drive this point home, let's show a trivial example:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "J\n",
      "a\n",
      "c\n",
      "k\n",
      "\n",
      "(0, 'J')\n",
      "(1, 'a')\n",
      "(2, 'c')\n",
      "(3, 'k')\n"
     ]
    }
   ],
   "source": [
    "for value in 'Jack':\n",
    "    print(value)\n",
    "print \n",
    "for iteration, value in enumerate('Jack'):\n",
    "    print(iteration, value)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A few other things to notice: we are using np.random again, as well as np.sqrt. To define an exponent python uses a double astrix. We also encounter an if statement. For conditionals, python drops the 'then' from the statement and uses tabbing to denote the begining and ending.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Writing A Function\n",
    "\n",
    "Lets take this same python script and make it into a simple function we can iterate over. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np #always good to be explicit about our libraries used"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def pitest():\n",
    "    # Create the coordinate\n",
    "    xtemp = np.random.uniform()\n",
    "    ytemp = np.random.uniform()\n",
    "    \n",
    "    # Calculate the distance to the origin\n",
    "    r = np.sqrt(xtemp**2 + ytemp**2)\n",
    "\n",
    "    # Count this if it is inside the circle.\n",
    "    if r< 1:\n",
    "        return True\n",
    "  \n",
    "    return False"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The Python syntax for functions use the same tabbing to start and end the function that loops and conditional statemetns do. This particular function pitest( ) does not have any variable input, but if it did, it would go in the parenthesis. Now let's put the new function in a loop and see how it works:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Pi=', 3.1228)\n"
     ]
    }
   ],
   "source": [
    "count = 0\n",
    "inside = 0\n",
    "\n",
    "for i in range(10000):\n",
    "    is_inside = pitest()\n",
    "    \n",
    "    inside += is_inside\n",
    "    count += 1.\n",
    "    \n",
    "print('Pi=',((inside/count)*4.))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Python has boolean data types, **True** and **False**.  These are also understood to have the integer values 1 and 0 respectively. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Running Your Own Function\n",
    "So now I have writen a script, how do I run it from the command line? Let's say that the function we wrote, __'pitest,'__ is saved with the name __'Pifunction.py'__. First, change to the directory where your program is located. "
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "cd /Users/mkirk/code/python/ "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Then import the name of your script. "
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "import Pifunction as Pifunc"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "And then run your function like we have all along:"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "Pifunc.pitest()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "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.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}