{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example 3 - Calculating the grain boundary resistivity and activation energy\n", "\n", "**Running the code in this notebook (under Mott-Schottky conditions for the 6 defined temperatures) takes approximately 3 minutes (iMac with 4 Ghz i7 processor).**\n", "\n", "The resistivity ratio has been calculated for the system where defects are specific to their sites, there is a site with two positively charged defects, followed by two sites which each have one negatively charged defect. The central positively charged defects have a segregation energy of -1.0 eV and the negatively charged defects are fixed to their bulk mole fractions under the Mott-Schottky approximation. The figure below shows the grain boundary resistivity for the system at a range of temperatures between 773.15 K and 1273.15 K. The resistivity decreases as temperature increases as expected. \n", "\n", "The function to calculate the resistivity ratio takes two arguments. The first argument is whether the space charge potential is positive or negative, and the second argument is a limit for the potential to cut off the sites included in the calculation when the potential is less than or greater than the limit.\n", "```\n", "c_o.calculate_resistivity_ratio_new('positive', -4e-2)\n", "```" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from pyscses.defect_species import DefectSpecies\n", "from pyscses.set_of_sites import SetOfSites\n", "from pyscses.constants import boltzmann_eV\n", "from pyscses.calculation import Calculation, calculate_activation_energies\n", "from pyscses.set_up_calculation import calculate_grid_offsets\n", "from pyscses.grid import Grid\n", "\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "boundary_conditions = 'periodic'\n", "site_charges = False\n", "systems = 'mott-schottky'\n", "core_models = False\n", "site_models = 'site_explicit'" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "alpha = 0.0005\n", "\n", "conv = 1e-8\n", "grid_x_min = -2.0e-8\n", "grid_x_max = +2.0e-8\n", "bulk_x_min = -2.0e-8\n", "bulk_x_max = -1.0e-8\n", "\n", "dielectric = 1\n", "\n", "index = 111\n", "\n", "b = 5e-9\n", "c = 5e-9\n", "\n", "temp = [ 773.15, 873.15, 973.15, 1073.15, 1173.15, 1273.15 ]" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "valence = [ +1.0, -1.0 ]\n", "site_labels = ['site_1', 'site_2']\n", "defect_labels = ['defect_1', 'defect_2']\n", "mole_fractions = [ [ 0.2, 0.2 ] ]\n", "initial_guess = [ [ 0.2, 0.2 ] ]\n", "mobilities = [ 1.0, 1.0 ]" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "data = '../input_data/example_data_2_one_seg_energies.txt'" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "--- 165.04243683815002 seconds ---\n" ] } ], "source": [ "par_resistivity_ratio_list = []\n", "per_resistivity_ratio_list = []\n", "\n", "limits, laplacian_limits = calculate_grid_offsets( data, grid_x_min, grid_x_max, 'single' )\n", "\n", "for m in mole_fractions:\n", " for t in temp:\n", " \n", " defect_species = { l : DefectSpecies( l, v, m, mob ) for l, v, m, mob in zip( defect_labels, valence, m, mobilities ) }\n", "\n", " all_sites = SetOfSites.set_of_sites_from_input_data( data, [grid_x_min, grid_x_max], defect_species, site_charges, core_models, t )\n", " if site_models == 'continuum':\n", " all_sites, limits = SetOfSites.form_continuum_sites( all_sites, grid_x_min, grid_x_max, 1000, b, c, defect_species, laplacian_limits, site_labels, defect_labels )\n", " if systems == 'mott-schottky':\n", " for site in all_sites.subset( 'site_2' ):\n", " site.defect_with_label('defect_2').fixed = True\n", " if systems == 'gouy-chapman':\n", " for site in all_sites.subset( 'site_2' ):\n", " site.defect_with_label('defect_2').fixed = False\n", " grid = Grid.grid_from_set_of_sites( all_sites, limits, laplacian_limits, b, c )\n", " \n", " c_o = Calculation( grid, bulk_x_min, bulk_x_max, alpha, conv, dielectric, t, boundary_conditions )\n", " c_o.form_subgrids( site_labels )\n", " if systems == 'gouy-chapman':\n", " c_o.mole_fraction_correction( m, systems, initial_guess )\n", " c_o.solve(systems)\n", " c_o.mole_fractions()\n", " c_o.calculate_resistivity_ratio( 'positive', 2e-2 )\n", " \n", " per_resistivity_ratio_list.append(c_o.perpendicular_resistivity_ratio)\n", " par_resistivity_ratio_list.append(c_o.parallel_resistivity_ratio)\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VVXWx/HvCiH03kQ6SJMqBEFKcAZpomAXuwzoCzaEKTrj6IxlHEdnQBxQbKPYHQQUUQRsCSAtinQQkC5SlSp9vX/cQ+YaUyE3N+X3eZ775N5zdvZdO4S7svc5Zx1zd0RERABioh2AiIjkHUoKIiKSQklBRERSKCmIiEgKJQUREUmhpCAiIimUFKRQMLO/mtlrwfPaZrbfzIqcZp8vm9kjORNhzjmd8ZnZWDO7P5LvIXmbkoJki5mtN7Ofgg+EbcEHY+lox5Ud7r7R3Uu7+/FoxxIJWR2fmd1sZrNSfe9gd384u+9hZp+b2aDTi1zyAiUFORUXu3tpoA0QD/w5ux2YWWyOR5VPWEiG//cK889HoktJQU6Zu28BpgLNAcysnJm9aGZbzWyLmT1ycnkh+Kt0tpmNNLNdwF/Dto02sz1mttLMup3sPwv9zTKzf5rZD2a2zsx6h31vPTNLNLN9ZjYDqBy2r66Z+ckPXjOraGYvmdl3QV/vhr9H+JiD7zsr9c/CzCqY2RQz2xH0McXMaobt/9zM/mZms4GDQP00+lhvZveY2WLggJnFmtmZZjYh6Hedmd0V1v5cM0s2s73BrG1EOuO72cy+DX4W68zsOjNrCowFzgtmfT8GbVOWxMxshZldFPZ+sUEcbcLfw8z+BnQBRgd9jTazMWb2r1Tjm2xmw37xiyR5ipKCnDIzqwVcCCwMNr0MHAPOAs4BegDhSwrtgW+BasDfwratJfSh/RdgoplVzEZ/q4LvfRx40cws2PcG8GWw72HgpgyG8ipQEmgGVAVGZj76X4gBXgLqALWBn4DRqdrcANwKlAE2pNPPNUAfoDxwAngfWATUALoBd5tZz6DtKGCUu5cFGgD/Td2ZmZUCngJ6u3sZoCPwtbuvAAYDc4JloPJpxPJmEM9JPYGd7v5VeCN3vw+YCdwR9HUHMA645uSMyMwqAxcQ+neRPExJQU7Fu8FflrOAROBRM6tGKEHc7e4H3H07oQ/X/mHf9527/9vdj7n7T8G27cCT7n7U3d8m9CHfJ4v9bXD354N17XFAdaCamdUG2gH3u/thd08i9OH6C2ZWHegNDHb3H4I4ErP7A3H3Xe4+wd0Puvs+Qkmva6pmL7v7smD8R9Pp6il33xT8fNoBVdz9IXc/4u7fAs+H/QyOAmeZWWV33+/uc9Pp8wTQ3MxKuPtWd1+WxWG9AfQ1s5LB62sJJYpMuft8YA+hREYQ8+fuvi2L7y1RoqQgp+ISdy/v7nXc/bbgA6wOUBTYamY/BknjWUJ/eZ+0KY2+tvjPqzJuAM7MYn/fn3zi7geDp6WD7//B3Q+k6jcttYDd7v5DJmPOkJmVNLNnzWyDme0FkoDy9vOzc9Iaf2rhbeoAZ54cf/Az+BOhmRbAQKARsNLMFoQv9ZwU/AyuJjQr2GpmH5hZk6yMyd3XACuAi4PE0Jfs/aU/Drg+eH49oRmZ5HE6mCU5ZRNwGKjs7sfSaZNWSd4aZmZhiaE2MDmL/aVnK1DBzEqFJYba6bz/JqCimZV39x9T7TtAaFkJADM7I4P3/C3QGGjv7t+bWWtCy2oW1iYrJYnD22wC1rl7wzQbuq/mf0s0lwHvmFmlNNpNA6aZWQngEUKzjS5ZjOfkElIMsDxIFJnFfdJrwFIzawU0Bd7NwvtJlGmmIDnC3bcC04F/mVlZM4sxswZmlnoJJbWqwF1mVtTMriT04fHhafSHu28AkoEHzSzOzDoDF2cQ91Tg6eBgcVEzSwh2LwKamVlrMysO/DWDty1D6DjCj8Exkb9kFmcWzAf2BQefS5hZETNrbmbtAMzsejOr4u4ngJMJ7UR4B2ZWzcz6BccWDgP7w9psA2qaWVwGMbxF6FjOEDKeJWwj1cFzd98MLCA0Q5gQtmQoeZiSguSkG4E4YDnwA/AOoXX+jMwDGgI7Ca3DX+Huu06jv5OuJXQgejehD+hXMmh7A6H1+ZWEjnHcDeDu3wAPAR8DqwkdQ0nPk0CJYBxzgY+yGGe6gmMlFwGtgXVB3y8A5YImvYBlZraf0EHn/ml88MYAw4HvCP0suhL6gAf4FFgGfG9mO9OJYSswh9AB6rczCHcUcEVw5tVTYdvHAS3Q0lG+YbrJjkSLmd0MDHL3ztGORSIjmHW9BtRxfdjkC5opiEhEmFlRYCjwghJC/qGkICI5Lrg47kdCy31PRjkcyYaILR+ZWWN+vgZZH3jA3X/xCxIcOJtDaE30nYgEJCIimcqVYwrBudpbCJ2utyGNfTOAQ8B/lBRERKInt65T6AasTZ0QAncCEwhdvZmpypUre926dXMwNBGRgu/LL7/c6e5VMmuXW0mhP2lcHm9mNYBLgV+RQVIws1sJ1Yyhdu3aJCcnRyhMEZGCyczSu6r/ZyJ+oDm4MKYvMD6N3U8C9wQX36TL3Z9z93h3j69SJdNEJyIipyg3Zgq9ga/SKYQVD7wVFLasDFxoZsfcXZfDi4hEQW4khWtIp7Kiu9c7+dzMXgamKCGIiERPRJePgnor3YGJYdsGm9ngSL6viIicmojOFIIKlZVSbRubTtubIxmLiIhkTlc0i4hICiUFERFJUWiSwq79h3no/eXsPZTeXRBFRKTQJIXZa3fx8hfr6D4ikY+X6zaxIiJpKTRJoW+rM3n39k5UKBnHoFeSueONr9i5/3C0wxIRyVMKTVIAaFmzPJPv6Mxvuzdi+rJtXDAikUkLN6NS7yIiIYUqKQDExcZwZ7eGfDi0M/Url2LY24sY8PICtvyo28eKiBS6pHDSWVXLMH5wR/568dnMX7ebHiMSeWXOek6c0KxBRAqvQpsUAIrEGDd3qsf0YQm0rVuRB95bxtXPzWHtjv3RDk1EJCoKdVI4qWaFkowb0I5/XdmKb7btp/eomYz5bA1Hj2dYvFVEpMBRUgiYGZe3rcnHw7vSvWk1npi2ir6jZ7N0y55ohyYikmuUFFKpUqYYY65rw7M3tGXX/sP0GzObv09dwaGjx6MdmohIxCkppKNnszOYMbwrV7atybOJ39J71Ezmfrsr2mGJiESUkkIGypUoymOXt+T1Qe05fsLp/9xc7pu0hH0qlSEiBZSSQhZ0Oqsy0+5O4JYu9Xhz/kZ6jEzikxUqlSEiBY+SQhaViCvCfX3OZuJtnShbvCgDxyVz15sL2aVSGSJSgCgpZFPrWuV5/87ODO/eiKlLt3LBiETeXbhFpTJEpEBQUjgFcbEx3NWtIR/c1YW6lUtx99tf85uXF/CdSmWISD6npHAaGlUrwzuDO/LARWcz99vd9BiZxKtzN6hUhojkW0oKp6lIjPGbzqFSGefULs/97y6l/3Nz+ValMkQkH1JSyCG1Kpbkld+cyxNXtGTl93vpNWomT3+uUhkikr8oKeQgM+PK+Fp8/NuudGtSlcc/WsUlY1QqQ0TyDyWFCKhapjjPXN+Wsde3Yfu+UKmMf3y0UqUyRCTPU1KIoF7Nq/PxsK5c3qYGz3y+lgtHzWT+ut3RDktEJF1KChFWrmRRHr+iFa8NbM/REye46tk53P/uUpXKEJE8SUkhl3RuGCqVMbBzPV6bt4EeI5P4bOX2aIclIvIzSgq5qGRcLPdfdDYThnSkTPFYBry8gLvfWsjuA0eiHZqICKCkEBVtaldgyp1duPuChnywJFQq472vVSpDRKJPSSFK4mJjuPuCRky5swu1KpZk6FtfM2hcMlv3qFSGiESPkkKUNT6jDBOHdOTPfZryxdpddB+RxGsqlSEiUaKkkAcUiTEGdanPtLsTaFWrHH9+dyn9n5/Lup0Hoh2aiBQySgp5SO1KJXltYHsev7wlK7bupdeTSYxNXMsxlcoQkVyipJDHmBlXtavFJ8O7cn7jKjw2dSWXPD2bZd+pVIaIRJ6SQh5VtWxxnr0hnmeua8P3ew7Td/RsnpimUhkiEllKCnlc7xbV+Xh4ApeeU4Mxn63lwqdmsmC9SmWISGRELCmYWWMz+zrssdfM7k7V5jozW2xmS8zsCzNrFal48rPyJeP455WteOU353Lk2AmuHDuHB95byv7Dx6IdmogUMBFLCu6+yt1bu3troC1wEJiUqtk6oKu7twAeBp6LVDwFQUKjKky7O4EBnery6twN9BiRyGerVCpDRHJObi0fdQPWuvuG8I3u/oW7/xC8nAvUzKV48q1SxWL5y8XNeGdwR0oWi2XASwsY9vbXKpUhIjkit5JCf+DNTNoMBKamtcPMbjWzZDNL3rFjR44Hlx+1rVOBD+7qzF3dGvL+ou/oPiKR9xd9p1IZInJaLNIfImYWB3wHNHP3bem0+RXwNNDZ3Xdl1F98fLwnJyfnfKD52Iqte7l3wmIWbd7DBU2r8cglzTmjXPFohyUieYiZfenu8Zm1y42ZQm/gqwwSQkvgBaBfZglB0ta0elkm3taJP/dpyqw1O+g+IpE35m1UqQwRybbcSArXkM7SkZnVBiYCN7j7N7kQS4EVXiqjeY1y/GnSEq59YS7rVSpDRLIhoknBzEoB3Ql98J/cNtjMBgcvHwAqAU8Hp61qXeg01alUijduac9jl7Vg2Xd76flkEs8nfctxzRpEJAsifkwhp+mYQtZt23uI+yYt5eMV22hdqzxPXNGShtXKRDssEYmCvHRMQaKkWtniPH9jW0b1b82GXQfo89Qsxny2hqMqsCci6VBSKODMjH6tazBjeFe6n12NJ6at4pIxs1n+3d5ohyYieZCSQiFRuXQxxlzXhmeua8O2vYfpO3oWI2Z8w5FjmjWIyP8oKRQyvVtUZ8awBPq2OpOnPlnNxf+exaJNP0Y7LBHJI5QUCqEKpeIYcXVr/nNzPHt+OsqlT8/msakqyy0iSgqF2q+bVGP68ASuiq/F2MRQWe4vN6gst0hhpqRQyJUtXpTHLm/JqwPP5fDRE1wxdg4Pvr+Mg0dUllukMFJSEAC6NKzC9GEJ3NChDi/NXk+vJ2fyxdqd0Q5LRHKZkoKkKFUslof6NeftWzsQY3Dt8/O4b9IS9h06Gu3QRCSXKCnIL7SvX4mpQxMY1Lkeb8zfSM+RSXyum/mIFApKCpKmEnFF+PNFZzNhSOhmPje/tIDfjV/EnoOaNYgUZEoKkqE2tSsw5c7O3P6rBkxauIXuIxOZsTzNKugiUgAoKUimihctwu97NuG92ztRsVQct7ySzF1vLtQtQEUKICUFybLmNcox+Y7ODLugEVOXbqX7iEQ+WLxVtwAVKUCUFCRb4mJjGHpBQ96/szM1KpTg9je+YshrX7F936FohyYiOUBJQU5JkzPKMnFIR+7t3YRPV22nx8gkJi3crFmDSD6npCCnLLZIDIO7NuDDu7pQv3Iphr29iIHjktm656dohyYip0hJQU7bWVVLM35wRx646Gy+WLuTHiOSeGv+Rs0aRPIhJQXJEUVijN90rse0uxNoVqMs905cwg0vzmfT7oPRDk1EskFJQXJUnUqleGNQBx65pDkLN/5AzyeTeGXOek6c0KxBJD9QUpAcFxNjXN+hDtOHdyW+bkUeeG8Z/Z+by7qdB6IdmohkQklBIqZG+RKMG9COJ65oyYrv99LrySSeT/qW45o1iORZSgoSUWbGlfG1+Hh4V7o0rMzfPlzB5c98wept+6IdmoikQUlBckW1ssV5/sZ4RvVvzYZdB+jz1CzGfLaGo8dPRDs0EQmjpCC5xszo17oGM4Z3pfvZ1Xhi2ioufXo2y7/bG+3QRCSQaVIws4vNTMlDckzl0sUYc10bnrmuDd/vOUTf0bMYMeMbjhzTrEEk2rLyYX81sNrMHjezJpEOSAqP3i2qM2NYV/q2OpOnPlnNxf+exaJNP0Y7LJFCLdOk4O7XA+cAa4GXzWyOmd1qZmUiHp0UeBVKxTHi6ta8eFM8P/50hEufns1jU1dy6OjxaIcmUihlaVnI3fcC7wBvAdWBS4GvzOzOCMYmhUi3ptWYPqwrV8XXYmziWi58aiZfbtgd7bBECp2sHFPoZ2aTgM+BosC57t4baAX8NrLhSWFSrkRRHru8Ja8OPJfDR09wxdg5PPj+Mg4eORbt0EQKjazMFC4DRrp7C3d/wt23A7j7QWBgRKOTQqlLwypMG5bADR3q8NLs9fR6ciZfrN0Z7bBECoWsJIXv3T0pfIOZ/QPA3T+JSFRS6JUuFstD/Zrz9q0diDG49vl53DdpCfsOHY12aCIFWlaSQvc0tvXO6UBE0tK+fiWmDk1gUOd6vDF/Iz1HJvH5qu3RDkukwEo3KZjZEDNbAjQxs8Vhj3XA4twLUQq7EnFF+PNFZzNhSEdKFovl5pcW8Lvxi9hzULMGkZxm6d0IxczKARWAvwP3hu3a5+5ROy0kPj7ek5OTo/X2EmWHjh7n35+uZmzit1QqFcffLm1B97OrRTsskTzPzL509/jM2mW0fOTuvh64HdgX9sDMKmYhgMZm9nXYY6+Z3Z2qjZnZU2a2JpiFtMmsXyncihctwu97NuG92ztRsVQct7ySzF1vLmT3gSPRDk2kQIjNYN8bwEXAl4ADFrbPgfoZdezuq4DWAGZWBNgCTErVrDfQMHi0B54JvopkqHmNcky+ozPPfL6W0Z+tZvaanTzUrzkXtjgDM8u8AxFJU7ozBXe/KPhaz93rB19PPjJMCGnoBqx19w2ptvcDXvGQuUB5M6uezb6lkIqLjWHoBQ15/87O1KhQgtvf+Iohr33F9n2Hoh2aSL6VlYvXJpvZNWZW8jTepz/wZhrbawCbwl5vDraljuFWM0s2s+QdO3acRhhSEDU5oywTh3Tknl5N+HTVdrqPSGLSws2kd7xMRNKXlVNS/wV0AVaY2TtmdoWZFc/qG5hZHNAXGH+KMeLuz7l7vLvHV6lS5VS7kQIstkgMQ85vwId3daFBlVIMe3sRg8Yl8/0ezRpEsiMrBfES3f02QscQngWuArJzonhv4Ct335bGvi1ArbDXNYNtIqfkrKqlGT+4I3/u05TZa3fSfWQi/03epFmDSBZlqSCemZUALgcGA+2Acdl4j2tIe+kIYDJwY3AWUgdgj7tvzUbfIr9QJMYY1KU+Hw1NoGn1svzhncXc9NICtvz4U7RDE8nzsnJM4b/ACuDXwGiggbtnqTqqmZUidEX0xLBtg81scPDyQ+BbYA3wPHBbtqIXyUDdyqV465YOPNi3Gcnrd9NzZBKvz9ugWYNIBtK9eC2lgVlP4GN3zxMF7nXxmpyKTbsPcs+ExXyxdhcdG1TiH5e3pFbF0zl3QiR/yerFaxld0fxrd//UzC5La7+7T0xre6QpKcipcnfemL+Rv3+4khPu3NOrCTd0qENMjK5rkIIvq0kho4vXugKfAhensc8JWxISyQ/MjOva1+H8xlW5d8Ji/jJ5GR8s2crjl7ekbuVS0Q5PJE/IyvJRPXdfl9m23KKZguQEd2d88mYe/mA5R4+f4Hc9GjOgUz2KaNYgBVRO1D46aUIa297JfkgieYeZcVW7WswY1pWODSrzyAcruOrZOazdsT/aoYlEVUals5uY2eVAOTO7LOxxM5Dli9dE8rIzyhXnxZviGXl1K9Zs30/vUTMZm7iWY8dPRDs0kajI6JhCY0IF8crz8+MK+4BbIhmUSG4yMy49pyadzqrM/e8u5bGpK5m6ZCtPXNmKRtXKRDs8kVyVlWMK57n7nFyKJ1M6piCR5O5MWbyVv0xexv5Dx7ir21n8X9cGFC2Spes8RfKsnDymcKmZlTWzomb2iZntMLPrcyBGkTzHzLi41ZnMGJZA92bV+Of0b7hkzGyWf7c32qGJ5IqsJIUe7r6X0FLSeuAs4PeRDEok2iqVLsaYa9sw9vo2bNt7iL6jZzFyxjccOaZjDVKwZSUpFA2+9gHGu/ueCMYjkqf0al6dGcO6clHL6oz6ZDV9R89iyWb9F5CCKytJ4X0zWwm0BT4xsyqA6hFLoVGhVBxP9j+H52+MZ/eBI1zy9Gwe/2glh4/licovIjkq0wPNkHJP5j3ufjy42U5Zd/8+4tGlQQeaJZr2HDzKwx8s550vN9OwammeuLIVrWuVj3ZYIpk67QPNZvbr4OtlwPlAv+B5L6BjDsUpkq+UK1mUf17ZipcGtGP/4WNc9vRs/v7hCg4d1axBCgbVPhI5Bb9qXJVpwxL4+4creDbpW2Ys38YTV7akbZ2K0Q5N5LRkafkoL9HykeQ1s1bv5J4Ji/luz08M6FiP3/VsRMm4jP7eEsl9OXadgpkNDa5TMDN7wcy+MrMeOROmSP7XuWFlpg1L4IYOdfjP7HX0HjWTud/uinZYIqckK2cf/Sa4TqEHUAm4AXgsolGJ5DOli8XyUL/mvHVrBwD6PzeX+99dyoHDx6IcmUj2ZCUpnKwlfCHwirsvC9smImE61K/E1KFd+E2nerw2bwM9RiYxa/XOaIclkmVZSQpfmtl0QklhmpmVAXRZp0g6SsbF8sDFZzP+/86jWGwM1784jz9OXMzeQ0ejHZpIprKSFAYC9wLt3P0gEAcMiGhUIgVAfN2KfDi0C7cm1OftBZvoOTKJz1dtj3ZYIhnKSlJw4GzgruB1KXQ/BZEsKV60CH+6sCkThnSkdLFYbn5pAb8bv4g9BzVrkLwpK0nhaeA84Jrg9T5gTMQiEimAzqldgSl3deb2XzVg0sItdB+ZyMfLt0U7LJFfyEpSaO/utxPUO3L3HwgtIYlINhSLLcLvezbh3ds6UbFUHINeSebutxbyw4Ej0Q5NJEVWksJRMytCaBmJoCCeDjSLnKIWNcsx+Y7ODO3WkCmLt9J9ZBIfLd0a7bBEgKwlhaeASUBVM/sbMAt4NKJRiRRwcbExDOveiMl3dKZa2WIMfu0rbn/jK3btPxzt0KSQy2qV1CZAN0LXJ3zi7isiHVh6VOZCCpqjx0/wbOJanvpkDaWLx/Jg32Zc1LI6ZrocSHJOVstcZJgUgmWjZe7eJCeDOx1KClJQfbNtH78fv4hFm/fQs1k1Hr6kOVXL6EQ/yRk5UvvI3Y8Dq8ysdo5FJiJpalStDBOGdOTe3k34bNUOuo9IYuJXm8lvRSslf8vKMYUKwDIz+8TMJp98RDowkcIotkgMg7s24MO7utCgSimG/3cRg8Yl8/0e3exQckemxxTMrGta2909MSIRZULLR1JYHD/hvDR7Hf+cvoqiRWK4v8/ZXBlfU8ca5JTkyDGFvEhJQQqb9TsP8IcJi5m/bjcJjarw98taUKN8iWiHJflMjt1PQUSiq27lUrx1Swce7NuM5PW76TkyidfnbdCxBokIJQWRfCAmxripY12m3Z1Ay5rluG/SUq57YR6bdh+MdmhSwCgpiOQjtSqW5PVB7Xn00hYs3ryHnk8mMe6L9Zw4oVmD5Ix0k4KZNTSzl81shJnVNLOpZnbAzBaZWbvcDFJE/sfMuLZ9baYNSyC+bkX+MnkZ/Z+by7qdB6IdmhQAGc0UXgK+AL4D5gH/IXQ7zt8Bo7PSuZmVN7N3zGylma0ws/NS7S9nZu8HiWaZmek+DSJZVKN8CcYNaMfjV7Rk5fd76fVkEs8mruXYcZUmk1OX7tlHZva1u7cOnq9x97PS2pdh52bjgJnu/oKZxQEl3f3HsP1/Asq5+z1Bob1VwBnunm7ZSJ19JPJL2/Ye4v53lzJ9+TZa1izHPy5vSdPqZaMdluQhOXH2UfifG3sz2JdeAOWABOBFAHc/Ep4QAg6UsdCJ16WB3YDudC6STdXKFufZG9oy+tpz2PLDT1z871mMmPENh48dj3Zoks9kNFM4CKwhVASvQfCc4HV9dy+VYcdmrYHngOVAK+BLYKi7HwhrUwaYDDQBygBXu/sHafR1K3ArQO3atdtu2LAhG0MUKVx2HzjCw1OWM2nhFhpVK83jV7Sida3y0Q5Louy0L14zszoZfaO7Z/jJbGbxwFygk7vPM7NRwF53vz+szRVAJ2A4ocQzA2jl7qlnJim0fCSSNZ+u3MZ9k5aybe8hBnaux/DujSkRVyTaYUmUnPbykbtvyOiRhRg2A5vdfV7w+h2gTao2A4CJHrIGWEdo1iAip+nXTaoxfVgC/c+tzfMz19FrVBJz1u6KdliSx2V0SupAM/t92OstZrbXzPaZ2eDMOnb374FNZtY42NSN0FJSuI3BdsysGtAY+DabYxCRdJQpXpRHL23Bm7d0AOCa5+fyp0lL2HfoaJQjk7wqowPNgwmdhnrSdncvC1QBrsli/3cCr5vZYqA18KiZDQ5LKg8DHc1sCfAJcI+778zWCEQkU+c1qMRHQxMY1Lkeb83fSI+RSXy2cnu0w5I8KKNjCsnh609m9id3fzR4Pt/dz82lGH9GxxRETs/CjT9wz4TFfLNtP5e0PpMHLm5GxVJx0Q5LIiwnTkn92ekKYQkhBqh8euGJSLScU7sC79/Zmbu6NWTK4q10H5HIlMXfqcCeABknhelm9kga2x8CpkcoHhHJBcViizC8eyPev7MzZ5YvwR1vLOT/Xv2SbXt1M5/CLqOk8HuggZmtMbMJwWMNcBahUhciks81rV6WSbd15I+9m5D4zQ4uGJHIfxds0qyhEMvKndfqA82Cl8vdfW3Eo8qAjimIRMa3O/Zz74QlzF+/my4NK/PopS2oVbFktMOSHKI7r4lItp044bw+bwOPTV2JA3/o2Zgbz6tLTIxuAZrf6c5rIpJtMTHGDefVTSnL/df3l3PVs3NYu2N/tEOTXKKkICK/ULNCScYNaMe/rmzF6u376T1qJmM+W8NRleUu8JQURCRNZsblbWsyY3gC3ZpU5Ylpq7hkzGyWfbcn2qFJBCkpiEiGqpYpzjPXt+WZ69qwbe9h+o2ezT+nreLQUZXlLoiUFEQkS3q3qM7HwxPo17oGoz9bQ5+nZvLlhh+iHZbkMCUFEcmy8iXj+NdVrXh5QDt+OnKcK8Z+wYPvL+PgEd0bq6BQUhCRbDu/cVWmD+/K9e3r8NLs9fR8MonZa1TLsiBQUhCRU1K6WCwPX9Kct2/tQGxMDNe9MI97Jyy/UuM5AAAQUElEQVRmz08qy52fKSmIyGlpX78SU4d24f+61ue/yZvoMTKRGcu3RTssOUVKCiJy2ooXLcIfezfl3ds7UaFkHLe8ksydby5k1/7D0Q5NsklJQURyTMua5Zl8R2eGXdCIj5ZupfvIJN77eosK7OUjSgoikqPiYmMYekFDptzZhVoVSzL0ra8ZNC6Z7/eoLHd+oKQgIhHR+IwyTBzSkfsubMrstTvpPiKRN+dv1Kwhj1NSEJGIKRJj3JJQn4+GJtCsRln+OHEJ170wj427DkY7NEmHkoKIRFzdyqV4Y1AH/nZpcxZv3kPPJ5N4cdY6jp/QrCGvUVIQkVwRE2Nc174O04cl0KF+RR6espwrxn7B6m37oh2ahFFSEJFcdWb5Evzn5nY8eXVr1u88QJ+nZvHvT1arLHceoaQgIrnOzLjknBrMGN6V7s2q8a8Z39B39GyWblFZ7mhTUhCRqKlcuhhjrm3Dsze0Zef+w/QbM5vHpq5UWe4oUlIQkajr2ewMPh7WlcvOqcHYxLVcOGomC9bvjnZYhZKSgojkCeVKFuWJK1vx6sBzOXzsBFc9O4e/vLeUA4dVljs3KSmISJ7SpWEVpg9L4Kbz6vLK3A30GJlE0jc7oh1WoaGkICJ5Tqlisfy1bzPG/995FCsaw43/mc/vxi9iz0GV5Y40JQURybPi61bkw7u6cNv5DZi0cAsXjEzko6XfRzusAk1JQUTytOJFi/CHXk147/ZOVC5djMGvfcntr3/Fjn0qyx0JSgoiki80r1GOyXd04nc9GjFj+TYuGJHI2wtUYC+nKSmISL5RtEgMd/y6IR8O7UzjamW4Z8IS+j83l7U79kc7tAJDSUFE8p2zqpbhrVs78PfLWrBi6156PzmTUR+v5vAxXfR2upQURCRfiokxrjm3Nh//tis9mlVj5Mff0OepWbro7TQpKYhIvla1THFGX9uGl25ux09HjnPl2Dn8ceISnb56iiKaFMysvJm9Y2YrzWyFmZ2XRpvzzexrM1tmZomRjEdECq5fNanKjOEJ3NKlHm8v2Ei3EYlMWfydDkRnU6RnCqOAj9y9CdAKWBG+08zKA08Dfd29GXBlhOMRkQKsZFws9/U5m8l3dKZ6ueLc8cZCBo5LZvMPutNbVkUsKZhZOSABeBHA3Y+4+4+pml0LTHT3jUGb7ZGKR0QKj+Y1yjHpto78uU9T5qzdRfcRSbww81uO6Z4NmYrkTKEesAN4ycwWmtkLZlYqVZtGQAUz+9zMvjSzG9PqyMxuNbNkM0vesUM1UEQkc7FFYhjUpT4zhidwXoNKPPLBCi55WvdsyEwkk0Is0AZ4xt3PAQ4A96bRpi3QB+gJ3G9mjVJ35O7PuXu8u8dXqVIlgiGLSEFTs0JJXrwpntHXnsP3ew7Td/QsHpmyXNVX0xHJpLAZ2Ozu84LX7xBKEqnbTHP3A+6+E0gidOxBRCTHmBkXtTyTT4Z35ep2tXlh1jp6jEzis5VasU4tYknB3b8HNplZ42BTN2B5qmbvAZ3NLNbMSgLtSXUwWkQkp5QrWZS/X9aC8YPPo0RcEQa8vIA73viK7fsORTu0PCPSZx/dCbxuZouB1sCjZjbYzAYDuPsK4CNgMTAfeMHdl0Y4JhEp5NrVrcgHd3VmePdGTF+2jQv+lcib8zdy4oROX7X8dg5vfHy8JycnRzsMESkg1u7Yz58mLmHeut2cW7cij17WnLOqlol2WDnOzL509/jM2umKZhEp1BpUKc1bt3bg8ctbsmrbPnqPmsnIGd8U2jpKSgoiUuiZGVe1q8Unv+3KhS2qM+qT1fQeNZN53+6Kdmi5TklBRCRQuXQxRvU/h3G/OZcjx05w9XNzuXfC4kJVR0lJQUQkla6NqjB9WAL/l1Cf8V9uptuIz5m8qHDUUVJSEBFJQ8m4WP54YVMm39GJM8uX4K43F3LzSwvYtLtg11FSUhARyUCzM8sx6bZO/OXis0lev5seI5N4Lmltga2jpKQgIpKJIjHGgE71mDG8K53OqsSjH66k35jZLN6cusZn/qekICKSRWeWL8HzN8bzzHVt2LHvMJeMmc1D7xesOkpKCiIi2WBm9G5RnY9/25Vr29fmpS9CdZQ+WbEt2qHlCCUFEZFTULZ4UR65pAXvDD6PUsWKMHBcMre//hXb9+bvOkpKCiIip6FtnYpMubMLv+vRiBkrttFtRCKvz9uQb+soKSmIiJymuNgY7vh1Q6bdnUDzM8tx36SlXPXsHL7Zti/aoWWbkoKISA6pV7kUb9zSnn9e2Yo1O/bT56mZjJi+ikNH808dJSUFEZEcZGZc0bYmnwzvysUtz+SpT9dw4aiZzFmbP+ooKSmIiERApdLFGHF1a14deC7HTjjXPD+X349fxA8HjkQ7tAwpKYiIRFCXhlWYdncCQ85vwMSFW7hgRCLvfb0lz9ZRUlIQEYmwEnFFuKdXE6bc2ZlaFUsy9K2vufE/89m4K+/VUVJSEBHJJU2rl2XCkI482LcZCzf+SI8nExmbuJajeaiOkpKCiEguKhJj3NSxLjOGJ5DQsAqPTV1J39GzWbQpb9RRUlIQEYmC6uVK8NyN8Yy9vi27Dxzmkqdn89fJy9gf5TpKSgoiIlHUq/kZfDy8Kzd2qMO4OevpPiKRGcujV0dJSUFEJMrKFC/Kg/2aM2FIR8oWL8otryQz5LUv2RaFOkpKCiIieUSb2hWYcldn/tCrMZ+u3M4F/0rk1bm5W0dJSUFEJA8pWiSG284/i+nDEmhVqzz3v7uUK8Z+warvc6eOkpKCiEgeVKdSKV4deC4jrmrF+l0H6fPUTF6ctS7i7xsb8XcQEZFTYmZc1qYm5zeuyt8+WEGdiiUj/p5KCiIieVzFUnH866pWufJeWj4SEZEUSgoiIpJCSUFERFIoKYiISAolBRERSaGkICIiKZQUREQkhZKCiIiksLx6n9D0mNkOYEMUQ6gM7Izi+0eDxlzwFbbxQuEbcx13r5JZo3yXFKLNzJLdPT7aceQmjbngK2zjhcI55qzQ8pGIiKRQUhARkRRKCtn3XLQDiAKNueArbOOFwjnmTOmYgoiIpNBMQUREUigpiIhICiWFVMxsmJktM7OlZvammRU3s3pmNs/M1pjZ22YWF7QtFrxeE+yvG93oT42ZDQ3Gu8zM7g62VTSzGWa2OvhaIdhuZvZUMObFZtYmutFnjZn9x8y2m9nSsG3ZHqOZ3RS0X21mN0VjLFmVzpivDP6dT5hZfKr2fwzGvMrMeoZt7xVsW2Nm9+bmGLIrnTE/YWYrg3/LSWZWPmxfvh9zjnN3PYIHUANYB5QIXv8XuDn42j/YNhYYEjy/DRgbPO8PvB3tMZzCmJsDS4GShO7E9zFwFvA4cG/Q5l7gH8HzC4GpgAEdgHnRHkMWx5kAtAGWhm3L1hiBisC3wdcKwfMK0R5bNsfcFGgMfA7Eh20/G1gEFAPqAWuBIsFjLVAfiAvanB3tsWVzzD2A2OD5P8L+nQvEmHP6oZnCL8UCJcwsltAH5Vbg18A7wf5xwCXB837Ba4L93czMcjHWnNCU0IfeQXc/BiQCl/HzsaUe8yseMhcob2bVczvo7HL3JGB3qs3ZHWNPYIa773b3H4AZQK/IR39q0hqzu69w91VpNO8HvOXuh919HbAGODd4rHH3b939CPBW0DZPSmfM04PfbYC5QM3geYEYc05TUgjj7luAfwIbCSWDPcCXwI9hv1SbCc0oCL5uCr73WNC+Um7GnAOWAl3MrJKZlST0V3ItoJq7bw3afA9UC56njDkQ/vPIb7I7xoI09tQKy5h/Q2gWCIVnzNmipBAmWFPuR2gqeSZQijz8l2BOcPcVhKbU04GPgK+B46naOFCgz10uDGMs7MzsPuAY8Hq0Y8nLlBR+7gJgnbvvcPejwESgE6Hlg9igTU1gS/B8C6G/qgn2lwN25W7Ip8/dX3T3tu6eAPwAfANsO7ksFHzdHjRPGXMg/OeR32R3jAVp7KkV6DGb2c3ARcB1wR8AUMDHfKqUFH5uI9DBzEoGxwa6AcuBz4ArgjY3Ae8FzycHrwn2fxr2C5dvmFnV4GttQscT3uDnY0s95huDM3Q6AHvClmDym+yOcRrQw8wqBLPKHsG2gmAy0D84o64e0BCYDywAGgZn4MUROqFichTjzDYz6wX8Aejr7gfDdhXYMZ+WaB/pzmsP4EFgJaG19lcJnZlQn9AvyxpgPFAsaFs8eL0m2F8/2vGf4phnEkp+i4BuwbZKwCfAakJnJFUMthswhtDZGUsIO4MlLz+ANwkdJzpKaI144KmMkdCa9JrgMSDa4zqFMV8aPD8MbAOmhbW/LxjzKqB32PYLCc0e1wL3RXtcpzDmNYSOEXwdPMYWpDHn9ENlLkREJIWWj0REJIWSgoiIpFBSEBGRFEoKIiKSQklBRERSKClIvhJUo/3azDaa2Y7g+dd5uUKtmV1mZk0i2H+MmX1mZqXNLNbMfgzbd3FQ7bOWmd1tZjdGKg4pGGIzbyKSd7h7e0i5QjXe3e+IbkQhZhbr/6uPldplwAlC17/kRH+pXQwku/v+sCvvCUpBjwS6u/smM3uBUMHDV7IahxQ+milIgWFmvc1sjpl9ZaH7XJQKtm82s0fNbJGZLTCzNmY23czWmtktQZsLgr+2pwZ/WY85WfE2k34fM7OFwKVmNjjof5GZjTezEmbWhdCFUCNPzmjMbJaZtQ76OMPM1gTPB5nZu2b2GcGV0mZ2r5nND+4F8EA6Q7+O/12NffJn8SvgaUIXZK0DcPf9wHeWT+6BIdGhpCAFQlCq415CV2S3ARYDQ8OarHP3VoRKJ79I6MrejsDDYW3aA0MI1dlvCvTLQr/b3f0cdx8PjHf3dsH7rAVudveZwIfAMHdv7e7rMxnKOcBl7t7NzC4EagdxtQY6mlnHNL6nI/BV2OuSwASgn7uvTtU2GeiSSQxSiGn5SAqKjoQ+zL8I/sCPA2aF7T9Zu2YJoRuuHAAOWOgOZKWDfXNPfmib2VtA52B7Rv2+Hfa8pZk9BJQHygBTTmEc0z10rwYI1VbqDSwMXpcGGgFfpPqesv7zmj6HgHnAAOC3qdpuB+qeQlxSSCgpSEFhwEfufkM6+w8HX0+EPT/5+uT/g9Q1XzwL/R4Ie/4KoeWapWY2iNBd29JyjP/N0otn0J8Bj7j7i+n0c9KJNF5fAXxqZn9w98fD9hUHfsqkPynEtHwkBcUXQFczqw9gZqXMrGE2++hgZrXNrAhwFaEZQXb6LQV8b2ZFgWvDtu8jNHM4aT3QNnh+BembBgwMO4ZR08wqp9FuTeqzr4KZUB9ggP38XtKNCBV7FEmTkoIUCO6+jVBFzLfNbBGhD/NG2exmPqF7cC8nVDVzcjb7fYBQ2eXZQR8nvQn8KezU2SeAoWb2FaF7Pac3pg8J3eZ1rpktIXSv8NJpNP0AOD+N799J6CZRD5pZn2DzeYQqw4qkSVVSRQidfQTc4e6XZNo4jzGzmsAL7p7hXQLNrB1wm7sPyJ3IJD/STEEkn3P3zcDLYQfM01MR+EsuhCT5mGYKIiKSQjMFERFJoaQgIiIplBRERCSFkoKIiKRQUhARkRT/Dw4kL4Q+oliHAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAasAAAEWCAYAAADYRbjGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecVdW5//HPdwaGMnQGsdCLQSxBGaoajcYajd1gVFAxaqJXTWISvUmu93rTNMX8cqMxxgaWWKMSY++oCAyKHWQEpKj0jvTn98fe6DhOOcDMnDMz3/frdV6zzzp7r/PsEXlYe6/1bEUEZmZmuSwv2wGYmZlVx8nKzMxynpOVmZnlPCcrMzPLeU5WZmaW85yszMws5zlZmWWRpOclnZtunyXppQyPu03SL2s3OpB0g6RfbMdx3SStlpRfW99hjUuTbAdglmskzQY6A5uBNcBjwEURsTqbcWVDRFyQyX7p7+zciHg6PW4O0Gpbv0PSwcAdEdFlm4O1Bs0jK7OKHRsRrYD9gGLg59vagaSc+MdgrsRhtiOcrMyqEBHzSUZWewFIOlvSe5JWSZop6fyt+0o6WNI8ST+V9Alwq6T2kh6RtEjSsnQ7o1GDpH6SnpK0VNJ0SadmeNxZkl6WdK2kJcB/p+3npLEvk/SEpO5pu9J9F0paKektSVvP97PLjZKK0viXpzGNl5Qn6XagG/Cv9NLfTyT1kBSSmkj6tqSScjH+QNK4st8hqTD9Xe+a9rNa0q6S1krqWObY/dLfZ9NMfh/WMDhZmVVBUlfgaOD1tGkhcAzQBjgbuFbSfmUO2RnoAHQHziP5f+zW9H034FPgLxl8byHwFHAXsBMwArheUv8MQx8CzCS5nPkrSccB/wmcCHQCxgP/SPc9HPgasDvQFjgVWFJBnz8C5qXHd077i4g4E5hDOhqNiGvKHfcv4CuS+pZp+056bp+JiDXAUcBHaT+tIuIj4Pk0pq3OBO6OiI0Z/i6sAXCyMqvYQ5KWAy8BLwC/BoiIf0fEB5F4AXgSOLDMcVuAKyNifUR8GhFLIuKBiFgbEauAXwEHZfD9xwCzI+LWiNgUEa8DDwCnZBj/RxHxf+mxnwIXAL+JiPciYlN6PgPS0dVGoDXQD1C6z8cV9LkR2AXoHhEbI2J8ZFBcNCLWAg8DpwGkSasfMC7DcxkDnJEem5/2c3uGx1oD4WRlVrHjI6JdRHSPiO+nf+Ej6ShJr6aXwZaTjLqKyhy3KCLWbX0jqaWkv0n6UNJK4EWgXQaz5LoDQ9JLbsvT7zqdZOSWibkV9Pf/yvS1FBCwW0Q8SzLauw5YKOlGSW0q6PN3QCnwZHoJ9PIMY4FkFHVauv0d4KE0iWXiYaC/pJ7AYcCKiJi0Dd9tDYCTlVmGJDUjGd38HugcEe2AR0n+0t+q/EjjR8BXgCER0YbkchvljqnIXOCFNGFufbWKiO9lGG75OOYC55frr0VEvAIQEX+OiIFAf5LLgT/+UocRqyLiRxHRC/gW8ENJh1byfeU9BXSSNIAkad1VyX5f6idN/veSjK7OxKOqRsnJyixzBUAzYBGwSdJRJPd7qtKa5D7VckkdgCsz/K5HgN0lnSmpafoaJGmP7Yz9BuAKSXsCSGor6ZR0e5CkIemEhTXAOpLLmV8g6RhJfSQJWEEytX/rfguAXpV9eXp/6T6S0VkHkuRVkQVAR0lty7WPBc4iSZJOVo2Qk5VZhtJ7TheT/Ct/GcnlrOruu/wJaAEsBl4FHt+G7zqcZGLFR8AnwNUkyXJ7Yn8wPf7u9HLk2ySTGSCZLPJ3knP6kGRyxe8q6KYv8DSwGpgAXB8Rz6Wf/Qb4eXqZ8bJKwrgL+AZwX3rfrKI4p5FM/JiZ9rVr2v4ySWJ8LSI+zPzMraGQH75oZvWBpGeBuyLipmzHYnXPycrMcp6kQSSXDrumo05rZHwZ0MxymqQxJJcfL3Wiarw8sjIzs5znkZWZmeU8F7isIUVFRdGjR49sh2FmVq9MmTJlcUR0qm4/J6sa0qNHD0pKSqrf0czMPiMpo6UIvgxoZmY5z8nKzMxynpOVmZnlPCcrMzPLeU5WZmaW85yszMws5zlZmZlZznOyyrLX5yzjhhc+yHYYZmY5zYuCs+zB1+czdsKHtGiaz6jhPbIdjplZTnKyyrL/OqY/H69Yx3//6x3aFxbwra/umu2QzMxyji8DZlmT/Dz+77R9GdSjAz+6dyovvr8o2yGZmeUcJ6sc0LxpPjeNKqbPTq254I4pTJ27PNshmZnlFCerHNGmeVPGnDOIolbNOPvWSZQu9DPmzMy2crLKITu1bs7toweTn5fHyJsn8dHyT7MdkplZTnCyyjHdOxYy5pxBrFq3iZG3TGLZmg3ZDsnMLOucrHLQnru25e+jipmzdC1n3zaZtRs2ZTskM7OscrLKUUN7deT/TtuXN+ct54I7XmPDpi3ZDsnMLGucrHLYEXvuzG9O3JsX31/EZfe9wZYtke2QzMyywouCc9y3B3Vj6ZqNXP34NDoUFnDlsf2RlO2wzMzqlJNVPXDBQb1Ysno9N700i6JWBVx0SN9sh2RmVqecrOoBSfzn0XuwdM0Gfv/k+3QobMZ3hnTLdlhmZnXGyaqeyMsTV5+8D8vWbuDnD71F+5ZNOWrvXbIdlplZnfAEi3qkaX4e158+kH27teeSu6fySunibIdkZlYnnKzqmRYF+dw8qpgeRS357tgS3pq3ItshmZnVOiereqhdywLGnjOEdi0LOOvWScxavCbbIZmZ1Sonq3pq57ZJHcEAzrx5IgtWrst2SGZmtcbJqh7r1akVt509iGVrNjDy5kmsWLsx2yGZmdUKJ6t6bp8u7bhxZDGzFq9h9JjJfLphc7ZDMjOrcU5WDcD+fYr404gBTJmzjAvveo2Nm11H0MwaFierBuLovXfhf4/bi2enLeSnD7zpOoJm1qB4UXADcsbQ7ixds4E/PvU+Ra2a8Z9H75HtkMzMaoSTVQPzH4f0Ycnq9dz44kw6FhZw/kG9sx2SmdkOc7JqYCRx5bF7snTtRn7z2DTaFxZwanHXbIdlZrZDnKwaoLw88YdTvsrytRu4/IE3ad+ygMP6d852WGZm280TLBqogiZ53HDGQPbu0o6L7nqNiTOXZDskM7Pt5mTVgBU2a8KtZw2iS/sWnDu2hHc/WpntkMzMtouTVQPXobCAsaOH0KpZE0bdOok5S9ZmOyQzs21Wq8lK0pGSpksqlXR5BZ83k3RP+vlEST3KfHZF2j5d0hHV9SnpzrT9bUm3SGqatveTNEHSekmXlfv+2ZLekjRVUkmZ9v+WND9tnyrp6Jr9zdSt3dq1YOw5g9m4eQtn3jKRhatcR9DM6pdaS1aS8oHrgKOA/sBpkvqX2200sCwi+gDXAlenx/YHRgB7AkcC10vKr6bPO4F+wN5AC+DctH0pcDHw+0pC/XpEDIiI4nLt16btAyLi0W3/DeSWvp1bc8tZg1i4cj1n3TKZletcR9DM6o/aHFkNBkojYmZEbADuBo4rt89xwJh0+37gUElK2++OiPURMQsoTfurtM+IeDRSwCSgS9q+MCImA43+b+f9urXnr2fsx/sLVvHdMSWs2+g6gmZWP9RmstoNmFvm/by0rcJ9ImITsALoWMWx1faZXv47E3g8gxgDeFLSFEnnlfvsIklvppcU21d0sKTzJJVIKlm0aFEGX5d9B39lJ/5w6leZOGspF//jdTa5jqCZ1QMNcYLF9cCLETE+g30PiIj9SC4rXijpa2n7X4HewADgY+APFR0cETdGRHFEFHfq1KkGQq8bxw3YjSuP7c+T7y7g5w+9TTIYNTPLXbW5KHg+ULZ0Qpe0raJ95klqArQFllRzbKV9SroS6AScn0mAETE//blQ0oMklxlfjIgFZfr8O/BIJv3VJ2fv35Olazbwf8+W0qGwgJ8c2S/bIZmZVao2R1aTgb6SekoqIJkwMa7cPuOAUen2ycCz6T2nccCIdLZgT6AvyX2oSvuUdC5wBHBaRFR7bUtSoaTWW7eBw4G30/e7lNn1hK3tDc0PD9ud0wZ34/rnP+Cm8TOzHY6ZWaVqbWQVEZskXQQ8AeQDt0TEO5KuAkoiYhxwM3C7pFKSWXsj0mPfkXQv8C6wCbgwIjYDVNRn+pU3AB8CE5I5GvwzIq6StDNQArQBtki6lGQmYRHwYLpvE+CuiNh6n+saSQNI7mnNJsORWn0jiV8evxfL1mzgl/9+j46tCjhh3y7ZDsvM7Evk+xU1o7i4OEpKSqrfMQet27iZs2+dzOTZS/n7yGK+3m+nbIdkZo2EpCkVLB36koY4wcK2UfOm+dw4ciD9dmnN9+6cwpQPl2U7JDOzL3CyMgBaN2/KbWcPZuc2zTnntsm8v2BVtkMyM/uMk5V9pqhVM24fPYRmTfI48+aJzFvmOoJmlhucrOwLunZoydjRg/l0w2ZG3jyJJavXZzskMzMnK/uyfju34eazBjF/+aecfdtkVq/flO2QzKyRc7KyCg3q0YHrT9+Pdz5ayfm3l7B+k+sImln2OFlZpQ7dozPXnLQPL5cu4Yf3vMHmLV7mYGbZUZvllqwBOGlgF5au2cCvHn2Pdi2b8svj9yJdSG1mVmecrKxa3/1aLxavWc/fXphJUatm/OCw3bMdkpk1Mk5WlpHLj+zH0tUb+H/PzKBjqwJGDuuR7ZDMrBFxsrKMSOI3J+7NsrUbuXLcO7RvWcCxX90122GZWSPhCRaWsSb5efzlO/syqHsHfnjvVF58v348cNLM6j8nK9smzZvm8/dRxfTu1IoL7pjC1LnLsx2SmTUCTla2zdq2aMrYcwbTsVUBZ986idKFq7Mdkpk1cE5Wtl12atOc288ZQn6eGHnzRD5e8Wm2QzKzBszJyrZbj6JCbjt7MKvWbeLMmyexbM2GbIdkZg2Uk5XtkL12a8uNI4uZs3Qt54yZzNoNriNoZjXPycp22LDeHfnziH15Y+5yLrjjNTZs2pLtkMysgXGyshpx5F478+sT9ubF9xfx4/vfYIvrCJpZDfKiYKsxIwZ3Y8maDfzuiem0b1nAlcf2dx1BM6sRTlZWo75/cG+WrN7ALS/PolPrZlz49T7ZDsnMGgAnK6tRkvj5N/dg2drPR1jfGdIt22GZWT3nZGU1Li9PXHPyPixbu4GfP/QW7Vs25ai9d8l2WGZWj3mChdWKpvl5XH/6fgzo2o5L7p7KKx8sznZIZlaPOVlZrWlZ0IRbzhpE944tOW/sFN6evyLbIZlZPeVkZbWqXcsCxo4eTNsWTRl1yyRmLV6T7ZDMrB5ysrJat0vbFowdPZgAzrx5IgtWrst2SGZWzzhZWZ3o3akVt509iGVrNjDy5kmsWLsx2yGZWT1SbbKSdKwkJzXbYft0acffzixm5uLVjLxlIotWrc92SGZWT2SShL4NzJB0jaR+tR2QNWwH9C3i+tMHMn3BKo6/7mWmfbIy2yGZWT1QbbKKiDOAfYEPgNskTZB0nqTWtR6dNUiH9e/MvecPY+PmLZz81wk8N21htkMysxyX0eW9iFgJ3A/cDewCnAC8Juk/ajE2a8D26dKOhy/an24dWjJ6zGRufXkWES5+a2YVy+Se1XGSHgSeB5oCgyPiKOCrwI9qNzxryHZp24L7LhjGIf068z//epf/evgdNm3240XM7MsyGVmdCFwbEXtHxO8iYiFARKwFRld1oKQjJU2XVCrp8go+bybpnvTziZJ6lPnsirR9uqQjqutT0p1p+9uSbpHUNG3vl166XC/psnLfP1vSW5KmSiop095B0lOSZqQ/22fwe7LtUNisCX87cyDnfa0Xt7/6IWffNpmV6zxT0My+KJNk9UlEvFi2QdLVABHxTGUHScoHrgOOAvoDp0nqX2630cCyiOgDXAtcnR7bHxgB7AkcCVwvKb+aPu8E+gF7Ay2Ac9P2pcDFwO8rCfXrETEgIorLtF0OPBMRfYFn0vdWS/LzxH8evQe/PXFvJnywhJOuf4W5S9dmOywzyyGZJKvDKmg7KoPjBgOlETEzIjaQ3O86rtw+xwFj0u37gUOVPADpOODuiFgfEbOA0rS/SvuMiEcjBUwCuqTtCyNiMrAt/1wvG9cY4PhtONa204jB3Rh7zmAWrFzH8de9zJQPl2Y7JDPLEZUmK0nfk/QW0E/Sm2Ves4A3M+h7N2Bumffz0rYK94mITcAKoGMVx1bbZ3r570zg8QxiDOBJSVMknVemvXNEfJxufwJ0rujgdFZkiaSSRYsWZfB1Vp3hfYp48ML9ad28CafdOJGHXp+f7ZDMLAdUNbK6CzgWeDj9ufU1MJ3OnquuB16MiPEZ7HtAROxHMlK8UNLXyu+QjtQqnKYWETdGRHFEFHfq1GmHgrbP9e7Uige/vz8DurXj0num8sen3vdMQbNGrqpkFRExG7gQWFXmhaQOGfQ9H+ha5n2XtK3CfSQ1AdoCS6o4tso+JV0JdAJ+mEF8RMT89OdC4EGSy4wACyTtkva5C+CFQHWsfWEBd4wewskDu/DnZ2bwH/94nXUbN2c7LDPLkupGVgBTgJL055Qy76szGegrqaekApIJE+PK7TMOGJVunww8m45kxgEj0tmCPYG+JPehKu1T0rnAEcBpEVHt/GdJhVsXNksqBA4H3q4grlEko0urYwVN8vjdyfvw0yP78cibHzPixlddosmskVJtXl6RdDTwJyAfuCUifiXpKqAkIsZJag7cTlIhYykwIiJmpsf+DDgH2ARcGhGPVdZn2r4J+JB09Af8MyKukrQzSXJtA2wBVpPMJCwiGU1B8sTku8r01RG4F+iW9nlqRFR5t7+4uDhKSjLJ4bY9Hn/7Yy69ZyodC5tx81nF9Nu5TbZDMrMaIGlKudnYFe9XXbKSNA74B/BwurbKKuBkVfvemreC0WMms2b9Jv7ynf34er+dsh2Sme2gTJNVJlPX/wAcCLwn6X5JJ6cjIrM6tXeXtjx80f5071joEk1mjUwmhWxfiIjvA72AvwGn4gkHliUu0WTWOGVUyFZSC+Ak4AJgEJ8vmDWrc1tLNJ3vEk1mjUYmhWzvBd4DDgH+AvSOCFdbt6zKzxNXlCvRNGeJb6maNVSZjKxuJklQF0TEc5lMCzerK18o0XT9y5TMdokms4aoqnJLh6SbhcBxkk4s+6qb8MyqN7xPEQ9duD9tmjfhO393iSazhqhJFZ8dBDxLUmKpvAD+WSsRmW2HXmmJpgvumMKl90xl5qLV/OCw3UnqIptZfVdpsoqIK9PNq9LK559Jq0qY5ZT2hQXcPnoIP3vwLf78bCkzF6/h96d8leZN87MdmpntoEzuWT1QQdv9NR2IWU0oaJLHNeVKNC1ctS7bYZnZDqp0ZCWpH8nDD9uWu0fVBvCiYMtZkvjewb3pWdSSS++ZygnXveISTWb1XFUjq68AxwDt+OIjQvYDvlv7oZntmCP32oX7zh/Opi1bOOn6V3humteym9VXmdQGHBYRE+oonnrLtQFz1ycr1jF6zGTe+3glvzimP2cN7+GJF2Y5oiZrA54gqY2kppKekbRIUi4/fNHsC3Zu25x7zx/GoXu4RJNZfZVJsjo8IlaSXBKcDfQBflybQZnVtMJmTfjbGQM5/6DPSzSt+NQlmszqi0ySVdP05zeB+yJiRS3GY1Zr8vLEFUftwdUnpSWa/uoSTWb1RSbJ6l+SpgEDgWckdQI8F9jqrW8P6sbY0YNZtGq9SzSZ1ROZPCLkcmA4UBwRG4E1wHG1HZhZbRreu4gHvz/cJZrM6olqawOma6wOJq0PCBxJkrzM6rWtJZr27daOS++Zyh+fnM6WLX6Yo1kucm1Aa9S2lmj6+UNJiaYPFq/hDy7RZJZzqq0NGBFn1104ZnWvoEkeV5+0D706teLqx6cxb9mn/H3kQHZq7UItZrkik4cvXpKus5KkmyS9JunwugjOrK5I4oKDevPX0wfy/ierOOG6V3jv45XZDsvMUpnMBjwnXWd1ONAROBP4ba1GZZYlR+61M/ddMIxNW7Zw8l9dosksV2SSrLbWpTkaGBsR75RpM2tw9tqtLQ9feAA9igoZPWYyt748i+rKkplZ7cokWU2R9CRJsnpCUmvAtWqsQdu5bXPuu2AY30hLNP3i4bddosksizJJVqOBy4FBEbEWKAA86cIavJYFTbghLdF0x6tzXKLJLIsySVYB9AcuTt8X4udZWSPhEk1muSGTZHU9MAw4LX2/Criu1iIyy0HfHtSN20cPcYkmsyzJJFkNiYgLSesBRsQykkuBZo3KsN4defD7w2nboinf+ftEHnx9XrZDMms0MklWGyXlk1wOJC1k6zvN1iglJZqGs1/3dvzgnjf4g0s0mdWJTJLVn4EHgZ0k/Qp4Cfh1rUZllsPatSxg7DlDOLW4C//3bCn/cffrrNu4OdthmTVoVdUGBCAi7pQ0BTiUZH3V8RHxXq1HZpbDtpZo6t2pFb91iSazWlflyEpSvqRpETEtIq6LiL84UZklJHH+Qb254QyXaDKrbVUmq4jYDEyX1K2O4jGrd47Y84slmp6dtiDbIZk1OJncs2oPvCPpGUnjtr4y6VzSkZKmSyqVdHkFnzeTdE/6+URJPcp8dkXaPl3SEdX1KenOtP1tSbdIapq295M0QdJ6SZdVEEO+pNclPVKm7TZJsyRNTV8DMjlfa7y2lmjq2amQc8eUcMtLLtFkVpOqvWcF/GJ7Ok5nEF4HHAbMAyZLGhcR75bZbTSwLCL6SBoBXA18W1J/YASwJ7Ar8LSk3dNjKuvzTuCMdJ+7gHOBvwJLSRY0H19JqJcA7wFtyrX/OCLu355zt8Zp57bNuff8Yfzgnqlc9ci7zFy8miuP3ZOm+Zn8m9DMqpLJY+1fqOiVQd+DgdKImBkRG4C7gePK7XMcMCbdvh84VJLS9rsjYn1EzAJK0/4q7TMiHo0UMAnokrYvjIjJwJfq5EjqAnwTuCmD8zGrVsuCJvz19IFccFBv7nh1Due4RJNZjajNf/LtBswt835e2lbhPhGxCVhB8hiSyo6tts/08t+ZwOMZxPgn4CdUvG7sV5LelHStpGYZ9GUGJCWaLj+qH9ectI9LNJnVkIZ4feJ64MWIGF/VTpKOARZGxJQKPr4C6AcMAjoAP62kj/MklUgqWbRo0Q6GbQ3NqYO6cvvoISxenZRoeuWDxdkOyazeqs1kNR/oWuZ9l7Stwn0kNQHaAkuqOLbKPiVdCXQCfphBfPsD35I0m+Ry4iGS7gCIiI/TK4rrgVtJLj9+SUTcGBHFEVHcqVOnDL7SGpukRNP+tEtLNF123xssWb0+22GZ1TuVJitJfdNZcX+U1EXSY5LWSHpD0qAM+p4M9JXUU1IByYSJ8rMIxwGj0u2TgWfTe07jgBHpbMGeQF+S+1CV9inpXOAI4LSIqLYcVERcERFdIqJH2s+zEXFG2tcu6U+RTMx4O4PzNatQz6JC/n3xgXzv4N489Pp8Dv3jC9w9aY7LNJltg6pGVrcCrwAfAROBW0juJ10G/KW6jtN7UBcBT5DMtrs3It6RdJWkb6W73Qx0lFRKMhq6PD32HeBe4F2Se08XRsTmyvpM+7oB6AxMSKeb/xeApJ0lzUv7/7mkeZLKz/wr705JbwFvAUXAL6s7X7OqtCjI56dH9uPRSw5k986tufyfb3HK3yYw7RMvIjbLhCpbCyJpakQMSLdLI6JPRZ9Zori4OEpKSrIdhtUDEcEDr83n14++x4pPNzL6gJ5ccmhfCptlspLErGGRNCUiiqvbr6qRVdlLaeX/+eeq62bbSRInD+zCMz88iFOLu3DjizM57I8v8MQ7n3ghsVklqhpZrSVZ3ySgd7pN+r5XRBTWSYT1hEdWtr1KZi/l5w+9zbRPVnFov53472/tSdcOLbMdllmdyHRkVVWy6l7VgRHx4XbG1iA5WdmO2Lh5C7e9PJtrn36fLRFccujunHtgT1e/sAZvh5OVbRsnK6sJ85d/yv+Me4cn313A7p1b8cvj92Zwzw7ZDsus1uzwPStJoyX9uMz7+ZJWSlol6YKaCtTMPrdbuxbcOLKYm0YWs2b9Zk792wR+fN8bLF2zIduhmWVVVdcYLiCZrr7VwohoQ7Lo9rRajcqskftG/8489cOv8b2De/Pg6/M55A/Pc89kr82yxquqZKWIWFLm/X0AEbEOaFGrUZkZLQuafL42a6fW/PQBr82yxquqZNWu7JuI+DWApDyShbJmVgd279yae84fyu9O3odZi9fwzT+/xK8ffY816zdlOzSzOlNVsnpSUkWVG64CnqyleMysApI4pbgrz/zwIE4Z+PnarCff+STboZnViaqS1Y+B3ukTeR9IX6VAH5KSS2ZWx9oXFvDbk/bh/guG0aZFU867fQrnjpnMvGV+BIk1bNVOXZfUi+SJvQDvRsQHtR5VPeSp61bXNm7ewq0vz+Lap2YAcMk3+jL6AK/NsvrF66zqmJOVZYvXZll9VhO1Ac2sHvDaLGsMnKzMGoita7MuOMhrs6zhcbIya0BaFjTh8qO+uDbrVK/NsgbAycqsASq7NuuDRas55s8v8ZtH32PtBq/NsvrJycqsgdq6NuvZHx3MyQO78LcXZ3LYH1/02iyrl5yszBq4smuzWjVrkq7NKvHaLKtXnKzMGoniHh145OIDuOKofrxcupjD/vgiN7zwARs3+8HflvucrMwakab5eZx/UG+e/tFBHNC3iN8+No1j/vwSk2cvzXZoZlVysjJrhHZr14K/jyzm7yOLWb1+E6fcMIGf3O+1WZa7nKzMGrHDyqzN+udr8zn0D89z7+S5XptlOcfJyqyR27o2698XH0ifnVrxkwfe5Ns3TmD6J6uyHZrZZ5yszAyAr+zcmnvOG8Y1J+9D6cLVfPPP4/nNY16bZbnBycrMPpOXJ05N12adtF8X/vZCsjbrqXcXZDs0a+ScrMzsS9oXFnD1yZ+vzfru2BKvzbKscrIys0p5bZblCicrM6vS1rVZT/3wa16bZVnjZGVmGenSvuWX1mb99P43Wea1WVYHnKzMbJtsXZt1/kG9eOC1eRzyh+e5t8Rrs6x2OVmZ2TZrWdCEK47a4/O1Wfd7bZbVLicrM9tun63NOslrs6x2OVmZ2Q7JyxOnDurKMz86mBP3281rs6z563hbAAATr0lEQVRWOFmZWY3oUFjANSd/lfsuGEZhs3y+O7aEU254hUfe/MhT3W2H1WqyknSkpOmSSiVdXsHnzSTdk34+UVKPMp9dkbZPl3REdX1KujNtf1vSLZKapu39JE2QtF7SZRXEkC/pdUmPlGnrmcZTmsZXUHO/FbOGbVCPDvz74gO58tj+LFi5novuep39f/ss/+/pGSxctS7b4Vk9VWvJSlI+cB1wFNAfOE1S/3K7jQaWRUQf4Frg6vTY/sAIYE/gSOD6NKlU1eedQD9gb6AFcG7avhS4GPh9JaFeArxXru1q4No0rmVpnGaWoab5eZy9f0+ev+xgbj1rEP13bcO1T7/P/r99lov/8TpTPlxKhGcPWuZqc2Q1GCiNiJkRsQG4Gziu3D7HAWPS7fuBQyUpbb87ItZHxCygNO2v0j4j4tFIAZOALmn7woiYDGwsH6CkLsA3gZvKtAk4JI2HNL7jd+xXYdY45eWJr/fbidvOHsxzlx3MyGE9eG76Qk766wSO/ctL3Dt5Lus2bs52mFYP1Gay2g2YW+b9vLStwn0iYhOwAuhYxbHV9ple/jsTeDyDGP8E/AQoe0G9I7A8jaeyuLd+13mSSiSVLFq0KIOvM2u8ehYV8otj+vPqFYfyqxP2YuOm4CcPvMnQ3zzDbx57j7lLXXfQKtcQJ1hcD7wYEeOr2knSMcDCiJiyvV8UETdGRHFEFHfq1Gl7uzFrVAqbNeH0Id15/NIDufu8oQzv3ZGbxs/ia797jnPHlDB+xiJfIrQvaVKLfc8HupZ53yVtq2ifeZKaAG2BJdUcW2mfkq4EOgHnZxDf/sC3JB0NNAfaSLqDZFTWTlKTdHRVUdxmtoMkMbRXR4b26sjHKz7lrolz+MekOTx98wJ6dSpk5NDunDSwC62bN812qJYDanNkNRnom86sKyCZMDGu3D7jgFHp9snAs+k9p3HAiHS2YE+gL8l9qEr7lHQucARwWkRUO082Iq6IiC4R0SPt59mIOCP9/ufSeEjje3j7fgVmlold2rbgR4d/hZcvP4Q/fXsAbVs05b//9S5Df/0Mv3jobWYscGWMxk61OdxORy1/AvKBWyLiV5KuAkoiYpyk5sDtwL4ks/ZGRMTM9NifAecAm4BLI+KxyvpM2zcBHwJb/1T/MyKukrQzUAK0Ibk3tRroHxEry8R5MHBZRByTvu9FMnmjA/A6cEZErK/qXIuLi6OkpGT7f1lm9gVvzlvOmFc+5F9vfsSGTVsY3rsjo4b34NB+O9EkvyHewWicJE2JiOJq9/O14ZrhZGVWO5asXs89JXO589U5zF/+Kbu1a8HpQ7sxYlA3OhR6CWR952RVx5yszGrXps1beGbaQsZOmM3LpUsoaJLHsfvsyqjh3dmnS7tsh2fbycmqjjlZmdWdGQtWcfurH/LAlHms2bCZAV3bMWp4d47eexeaNcnPdni2DZys6piTlVndW7VuI/98bT5jJsxm5qI1dCws4LTB3Th9aDd2adsi2+FZBpys6piTlVn2RAQvly7htldm88y0BeRJHLFnZ0YO68GQnh1ICtNYLso0WdXmOiszszohiQP6FnFA3yLmLl3LHRM/5J7Jc3n0rU/4SufWjBzeneMH7EZhM/+VV195ZFVDPLIyyy3rNm5m3BsfMeaV2bzz0UpaN2/CKQO7cuaw7vQsKsx2eJbyZcA65mRllpsigtfmLGfshNk8+tbHbNwcHLR7J0YN787Bu+9EXp4vEWaTk1Udc7Iyy30LV63j7klzuXPihyxYuZ5uHVpy5tDunFrclbYtXdYpG5ys6piTlVn9sXHzFp58ZwFjJsxm0qylNG+axwn77saZQ3vQf9c22Q6vUXGyqmNOVmb107sfreT2V2fz4OvzWbdxC4N7dGDk8O4csefONHVZp1rnZFXHnKzM6rcVazdy35S5jJ3wIXOWrmWn1s04fUh3ThvSlZ1aN892eA2Wk1Udc7Iyaxi2bAleeH8RYybM5vnpi2iaL47aaxdGDe/Oft3ae81WDfM6KzOz7ZCXJ77ebye+3m8nZi1ewx2vfsi9JXMZ98ZH7LVbG0YO7cG3BuxK86Yu61SXPLKqIR5ZmTVcazds4qHXkzVb0xesol3Lpnx7UFfOGNKdrh1aZju8es2XAeuYk5VZwxcRTJy1lLETZvPEOwvYEsGh/Tozanh3DuhT5EuE28GXAc3MapgkhvbqyNBeHfl4xafcNXEO/5g0h6dvXkCvToWMHJpUft+pjSdk1DSPrGqIR1ZmjdP6TZt57K1PGDNhNq/PWQ7AVzq3/qxW4ZCeHWhZ4HFBZXwZsI45WZnZtE9W8vz0Rbw0YzGTZi9lw6YtFOTnMbB7ew7oW8SBfYvYc9e25LvE02ecrOqYk5WZlbVu42Ymz17K+BmLGT9jMe99vBKAdi2bsn/vJHEd0LeILu0b9wQN37MyM8ui5k3zObBvJw7s2wmARavW88oHi3nx/cW8VLqIf7/1MQA9iwo5oE+SuIb17kib5q5RWBGPrGqIR1ZmlqmIoHThasbPWMxLpYt5deYS1m7YTH6eGNC1HQf0SUZeA7q2o0kDL/nky4B1zMnKzLbXhk1beG3OMl6asZjxpYt5a95ytgS0btaEob07JpcM+xTRs6iwwU2Pd7KqY05WZlZTlq/dwCsfLEnvdy1i3rJPAditXYvP7nXt37uI9oUFWY50xzlZ1TEnKzOrDRHBh0vWMr50MS/NWMQrHyxh1bpNSLDXrm0/S14Du7enWZP6VwLKyaqOOVmZWV3YtHkLb85fwfh0osbrc5azaUvQomk+g3t24MC+RRzYtxO7d25VLy4ZOlnVMScrM8uG1es38eoHSxg/YxHjSxczc9EaAHZq3eyzWYYH9CnK2aoaTlZ1zMnKzHLB/OWf8vKMxbw4YxEvly5m2dqNAPTbufVnyWtIz460KMiNS4ZOVnXMycrMcs2WLcG7H69Mp8gvYvKsZWzY/HlVjQN3L+LAPp3Yc9c25GWpqoaTVR1zsjKzXPfphs1Mmr2Ul2YsYvyMxUz7ZBUA7Vs2ZXifIg7sU/dVNVzBwszMvqBFQT4H7d6Jg3ZPqmosXLWOl0uTclAvzVjMv99Mqmr0Kir87F7XsN4daZ0DVTU8sqohHlmZWX0WEcxYuJoX31/ES6WLmThzKZ9uTKpq7Nu13WeFeL/apWaravgyYB1zsjKzhmT9ps289uFyxs9Iktdb81cQaVWNYVuravTtRI+OLXdoiryTVR1zsjKzhmzZmq1VNZL7XfOXf15V43cn78PwPkXb1W+myapWKyRKOlLSdEmlki6v4PNmku5JP58oqUeZz65I26dLOqK6PiXdmba/LekWSU3T9n6SJkhaL+myMvs3lzRJ0huS3pH0P2U+u03SLElT09eAmv/tmJnVH+0LC/jmPrvw25P24aWffp3nLjuY/z1uT/bctQ27tGtR699faxMsJOUD1wGHAfOAyZLGRcS7ZXYbDSyLiD6SRgBXA9+W1B8YAewJ7Ao8LWn39JjK+rwTOCPd5y7gXOCvwFLgYuD4ciGuBw6JiNVpYntJ0mMR8Wr6+Y8j4v6a+W2YmTUckuhZVEjPokLOHNajTr6zNkdWg4HSiJgZERuAu4Hjyu1zHDAm3b4fOFTJxc/jgLsjYn1EzAJK0/4q7TMiHo0UMAnokrYvjIjJwMayX5zuujp92zR9+ZqomVkOqs1ktRswt8z7eWlbhftExCZgBdCximOr7TMdJZ0JPF5dgJLyJU0FFgJPRcTEMh//StKbkq6V1KyS48+TVCKpZNGiRdV9nZmZbaeG+FSv64EXI2J8dTtGxOaIGEAyChssaa/0oyuAfsAgoAPw00qOvzEiiiOiuFOnTjUTvZmZfUltJqv5QNcy77ukbRXuI6kJ0BZYUsWxVfYp6UqgE/DDbQk0IpYDzwFHpu8/Ti8TrgduJbn8aGZmWVKbyWoy0FdST0kFJBMmxpXbZxwwKt0+GXg2vec0DhiRzhbsCfQluQ9VaZ+SzgWOAE6LiC3VBSepk6R26XYLkkkb09L3u6Q/RTIx4+3t/B2YmVkNqLXZgBGxSdJFwBNAPnBLRLwj6SqgJCLGATcDt0sqJZm1NyI99h1J9wLvApuACyNiM0BFfaZfeQPwITAhXaD2z4i4StLOQAnQBtgi6VKgP7ALMCadtZgH3BsRj6R93SmpEyBgKnBBLf2azMwsA14UXEO8KNjMbNvlxKJgMzOzmuCRVQ2RtIjkMmS2FAGLs/j9da2xnS/4nBuLxnbO3SOi2unUTlYNhKSSTIbSDUVjO1/wOTcWjfGcM+HLgGZmlvOcrMzMLOc5WTUcN2Y7gDrW2M4XfM6NRWM852r5npWZmeU8j6zMzCznOVmZmVnOc7KqJyT9IH2i8duS/pE+6bhn+oTl0vSJywXpvpU+gbk+kXRJer7vpGWykNRB0lOSZqQ/26ftkvTn9JzflLRfdqPPTPpU64WS3i7Tts3nKGlUuv8MSaMq+q5cUck5n5L+d94iqbjc/tv01PBcVMk5/07StPS/5YNba5Wmn9X7c65xEeFXjr9Intk1C2iRvr8XOCv9OSJtuwH4Xrr9feCGdHsEcE+2z2E7znkvkgLCLUlqWD4N9AGuAS5P97kcuDrdPhp4jKSe41BgYrbPIcPz/BqwH/B2mbZtOkeSx9jMTH+2T7fbZ/vctvGc9wC+AjwPFJdp7w+8ATQDegIfkNQFzU+3ewEF6T79s31u23jOhwNN0u2ry/x3bhDnXNMvj6zqjyZAi/RRKi2Bj4FDSJ6wDMkTl49Ptyt7AnN9sgfJX8ZrI3kw5wvAiXzx3Mqf89hIvAq021o9P5dFxIskRZzL2tZzPILk4aFLI2IZ8BTp425yUUXnHBHvRcT0Cnbf5qeG56JKzvnJ9M82wKukTzengZxzTXOyqgciYj7we2AOSZJaAUwBlpf5w172qcmVPYG5PnkbOFBSR0ktSUYVXYHOEfFxus8nQOd0O5MnU9cX23qODency2ss53wOyagZGs85bxMnq3ogvWdxHMklgV2BQnL4X841ISLeI7k08iTwOMmjWjaX2yeABr32ojGcY2Mn6Wckj0K6M9ux5DInq/rhG8CsiFgUERuBfwL7k1wG2vpMsrJPTa7sCcz1SkTcHBEDI+JrwDLgfWBBmYdj7gIsTHfP5MnU9cW2nmNDOvfyGvQ5SzoLOAY4Pf2HCTTwc95eTlb1wxxgqKSW6b2nQ0keTPkcyROWIXni8sPpdmVPYK5XJO2U/uxGcr/qLr54buXPeWQ6Y24osKLMpbT6ZlvP8QngcEnt01H44WlbQ7DNTw2vLyQdCfwE+FZErC3zUYM95x2S7RkefmX2Av4HmEZyL+d2kplCvUj+EJcC9wHN0n2bp+9L0897ZTv+7Tzn8SRJ+Q3g0LStI/AMMINkhmCHtF3AdSSzpd6izIyyXH4B/yC5D7mR5B7E6O05R5J7HqXp6+xsn9d2nPMJ6fZ6YAHwRJn9f5ae83TgqDLtR5OMtj8Afpbt89qOcy4luQc1NX3d0JDOuaZfLrdkZmY5z5cBzcws5zlZmZlZznOyMjOznOdkZWZmOc/JyszMcp6TlVkNSKvbT5U0R9KidHtqLle8l3SipH612H+epOcktZLURNLyMp8dm1YP7yrpUkkjaysOaxiaVL+LmVUnIobAZxUJiiPiouxGlJDUJD6vH1neicAWkvV7NdFfeccCJRGxukylFdJHXlwLHBYRcyXdRFKoeGymcVjj45GVWS2TdJSkCZJeU/KcscK0fZ6kX0t6Q9JkSftJelLSB5K+m+7zjXR08lg6ErluawX9avr9raTXgRMkXZD2/4ak+yS1kHQgyQLTa7eOACW9JGlA2sfOkkrT7XMlPSTpOdLKGJIulzQpfRbTf1Vy6qfzefWNrb+LrwPXkyx0nQUQEauBj1RPnkFm2eFkZVaL0pJRl5NU4NgPeBO4pMwusyLiqySPiLiZpJLDcOB/y+wzBPgeyXOO9gCOy6DfhRGxb0TcB9wXEYPS7/kAOCsixgOPAj+IiAERMbuaU9kXODEiDpV0NNAtjWsAMFzS8AqOGQ68VuZ9S+AB4LiImFFu3xLgwGpisEbMlwHNatdwkiTzSjogKgBeKvP51tpub5E8iG8NsEbJE3NbpZ+9ujWZSLobOCBtr6rfe8ps7yPpKqAd0Bp4ZDvO48lInpUFSe3Bo4DX0/etgN2BV8od0ya+WPNuHTAROBv4Ubl9FwI9tiMuayScrMxql4DHI+LMSj5fn/7cUmZ76/ut/3+Wr4kWGfS7psz2WJLLbm9LOpfkKcMV2cTnV1uaV9GfgF9GxM2V9LPVlgrenww8K+knEXFNmc+aA59W0581Yr4MaFa7XgEOktQLQFKhpL7b2MdQSd0k5QOnkoygtqXfQuATSU2B75RpX0Uy0tpqNjAw3T6Zyj0BjC5zj6yLpKIK9istPxsyHTl+Ezhb0qgyH+1OUqTZrEJOVma1KCIWkFTYvkfSGyRJZvdt7GYScANJBfrpwLht7Pe/SB4v8XLax1b/AP6zzBT73wGXSHoNaF/FOT0K3A+8Kukt4F6SS4Hl/Rs4uILjF5M8PPR/JH0zbR5GUmnerEKuum6WwyR9A7goIo7PdizbSlIX4KaIqPKp1pIGAd+PiLPrJjKrjzyyMrNaERHzgNvKTBSpTAfgyjoIyeoxj6zMzCzneWRlZmY5z8nKzMxynpOVmZnlPCcrMzPLeU5WZmaW8/4//iNXpfiHDGYAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot( temp, per_resistivity_ratio_list)\n", "plt.title('Perpendicular resistivity')\n", "plt.xlabel( ' Temperature (K)' )\n", "plt.ylabel( ' GB resistivity' )\n", "plt.show()\n", "\n", "plt.plot( temp, par_resistivity_ratio_list)\n", "plt.title('Parallel resistivity')\n", "plt.xlabel( ' Temperature (K)' )\n", "plt.ylabel( ' GB resistivity' )\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The activation energy can then be calculated from the Arrhenius relationship as described above using the ```calculate_activation_energy``` function in the source code. The function takes the resistivity ratio and temperature as input, calculates $x$ as $\\frac{1}{T}$ and $y$ as $ \\mathrm{ln} \\left( \\frac{1}{r_{GB}}\\right) $. The slope at each point is calculated using a central difference method and the activation energy is calculated as the slope * Boltzmann constant. " ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEWCAYAAABMoxE0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VNX5x/HPl7DvsipLWIPsKERA61YVxQVR0Z9oVbRWXEoX64ZVUVFbtVqr1WqxLtXWqqBicNeKWjckIEvCGvZN9n0JJHl+f9wbHWOWiclkJsnzfr3mlTv3nnvnOZlknjn33HuOzAznnHPux6oR7wCcc85Vbp5InHPOlYknEuecc2XiicQ551yZeCJxzjlXJp5InHPOlYknEldpSbpD0r/C5WRJuyQllfGYz0q6u3wiLD/lVb8ijv0zSe+V93Fd9eGJxCFpuaS94QfV+vDDtGG84yoNM1tpZg3NLDfesZSH8D05Kf95edVPUkdJJqlmxLH/bWYnl+W4rnrzROLyDTOzhkB/IBW4tbQHiPxwqm4U8P+nOKjOf3eJwv/w3feY2RrgbaA3gKQmkp6StE7SGkl3559ekXSppM8kPSRpM3BHxLpHJW2XtEDSifnHj+J4n0p6QNJWScsknRqxbydJH0vaKel9oEXEtu9905bUTNIzktaGx5oc+RqRdQ7361rwdyHpIElvSNoYHuMNSe0itn8k6R5JnwF7gM6FHGOspCVhzPMknV1g+xWS5kds7y/peSAZmBK2Em+MrJ+k8yWlFzjOtZLSwuXTJX0taYekVZLuiCj6SfhzW3jsIwv+TiQdJWl6+P5Nl3RUgTrfFb7HOyW9J6kFRZB0hqRZkrZJ+lxS34htyyVdL2lO+FovSapbin1vkjQH2B3+XvqH9d4paWJ4vLvD8hmShkXsX0vSJkmHFxW7KwUz80c1fwDLgZPC5fZAJnBX+Pw14O9AA6AV8BVwZbjtUiAH+BVQE6gXse5aoBZwPrAdaBbl8Q4AVwBJwNXAWkDh9i+APwN1gGOBncC/wm0dAQNqhs/fBF4CDgrjOC7iNT4tUH8DuobLzwJ3h8vNgRFAfaARMBGYHLHfR8BKoFdY/1qF/G7PA9oQfGk7H9gNHBKxbQ1wBCCgK9Ch4HtSsH5hPDuBlIjt04GR4fLxQJ/wNfsC64GzCvs9FfydAM2ArcDF4WtdED5vHlHnJUC38P3+CLi3iL+rw4ENwKDw/RwV1qtORB2/Cn8/zYD5wFWl2HcWwd9rPaA2sAL4Tfh+nwPsj3gvbwReiohtODA33v97VeUR9wD8Ef9H+E+5C9gW/jP+LfznbA1kA/Uiyl4ATA2XLwVWFjjWpUR8+Ifrvgo/mKI5XlbEtvrhh97BBN/Qc4AGEdtfoJBEAhwC5AEHFVLXbz80I9YVmkgK2fcwYGvE84+A8aX8Xc8ChofL7wK/KeY9KTSRhM//BYwLl1MIEkv9Io71F+Chwo5T8HcSvk9fFdj/C+DSiDrfGrHtGuCdIl73ccIvJBHrFvJdUl8OXBSx7X7giVLs+/OIbccSJOXIv7tP+S6RtAl/R43D55OAG+P9v1dVHn5u0eU7y8w+iFwhqQ/Bt7t1kvJX1wBWRRSLXM63xsL/1tAKgn/kDlEc75v8BTPbE5ZrSHAaa6uZ7S5w3PaFvH57YIuZbS1kW9Qk1QceAoYStGwAGklKsu86vQurf+QxLgF+R/ABDt/VJT/OJT8yvBeAB4HxwIUELaU94WsOAu4lOD1Zm6AFNzHK47Yh+L1GWgG0jXj+TcTyHoI6FaYDMErSryLW1Q5fo6hj5W+LZt/I330bfvh39+12M1sbnoIcIek14FSC1osrB95H4oqziqAF0cLMmoaPxmbWK6JMYcNHt1VEpiBoTayN8nhFWQccJKlBgeMWFXczSU0L2baboKUDgKSDi3nN64BDgUFm1pjgWy8Ep6HyFTl8tqQOwJPAGIJTQ02BjIj9VwFditi9pGG53wdaSjqMoFX3QsS2F4A0oL2ZNQGeiHjNko67luBDPFIywbf90loF3BPxXjc1s/pm9p9y2jeyLuv44d9dwS8Z/wQuIjil+IUF/YGuHHgicUUys3XAe8CDkhpLqiGpi6TjSti1FfDrsEPzPKAH8FYZjoeZrQDSgTsl1ZZ0NDCsiLLrCC4Y+FvYYV5LUn4SmA30knRY2LF7RzEv2wjYS9Ax3Qy4vaQ4C2hA8GG3EUDSZYQXMYT+AVwvaYACXcPkA0G/xg867yPqeICglfEngv6F9wvEvcXM9kkaSNBiybeR4LRfUcd+C+gm6cL8jn2gJ/BGVDX+vieBqyQNCuvXILwQoFEM9v0CyAXGhHEPBwYWKDOZ4KrE3wDP/Yj6uCJ4InEluYTglMI8gk7XSQR9EMWZRnDefhNwD3CumW0uw/HyXUjQ+bqF4EO9uA+Diwk67hcQdNr+FsDMFhGcDvoAWExwHr0ofyHoK9oEfAm8E2WchK81j+D00xcEiaEP8FnE9okEv58XCM7fTyZICgB/BG4Nr1i6voiXeAE4CZhoZjkR668BxkvaCYwDXo54zT3ha34WHntwgZg3A2cQtMY2E3RSn2Fmm0pT9/BY6QQXTjxK8F5nEfTHlPu+ZrafoIP9coK+vosIkl92RJm9wCtAJ+DVUlbHFSP/ahjnyoWkS4FfmNnR8Y7FVW+SphF03j8TsW4c0M3MLopfZFWPt0icc1WCpOMkHRye2hpFcOnzOxHbmxG0WCbEK8aqyhOJc66qOJSgD2wbwam5c8P+MiRdQdCB/7aZfVL0IdyP4ae2nHPOlYm3SJxzzpVJtbghsUWLFtaxY8d4h+Gcc5XKjBkzNplZy5LKVYtE0rFjR9LT00su6Jxz7luSCo5yUCg/teWcc65MPJE455wrE08kzjnnysQTiXPOuTKJaSKRNFTSQklZksYWsr1OOItZlqRpkjqG6weGM6PNkjRb4axyktpLmqpgJrlMST4MtHPOxVnMEomC6VMfIxj3vydwgaSeBYpdTjDHRFeCeR/uC9dnAKlmdhjBXBB/VzCFag5wnZn1BAYDvyzkmM455ypQLFskAwlmu1sajsz5IsH0lpGGE8wRAMEosCdKkpntiRjNtC7hvANmts7MZobLOwmm5myLc865uIllImnL92cwW80PP/S/LRMmju0E82QTzkOQCcwlmMc5cphswtNghxMMWf4DkkZLSpeUvnHjxjJXxjnnKpNF63dy/zsLqIhhsBK2s93MpoUz5x0B3BxOQgSApIYE8wr81sx2FLH/BDNLNbPUli1LvDHTOeeqhNw844mPl3DGI5/y4vRVrNu+L+avGcs729fw/aku2/HD6Trzy6wO+0CaEEym8y0zmy9pF8HMcumSahEkkX+bmU9O45xzoaUbd3H9xNnMXLmNob0O5u6ze9OiYZ2Yv24sE8l0IEVSJ4KEMZLvT/kJwbzSowhmkDsX+NDMLNxnlZnlhFOPdgeWh/MxPwXMN7M/xzB255yrNPLyjGc/X8797y6gTs0kHh55GGf2a8P3p7CPnZglkjAJjAHeBZKAp80sU9J4IN3M0giSwvOSsgimTx0Z7n40MFbSAYL5pa8xs03hPN0XA3MlzQrL/t7M3opVPZxzLpGt3LyHGybNZtqyLZzQvRV/PKcPrRvXLXnHclQt5iNJTU01H7TROVeVmBn/nraSP7w1nySJ24b15LwB7cq1FSJphpmlllSuWoz+65xzVcnabXu56ZU5/G/xJo5JacG9I/rStmm9uMXjicQ55yoJM2Ni+mruemMeuWbcfVZvfjYoucL6QoriicQ55yqB9Tv2cfOrc/lwwQYGdWrGn87tR3Lz+vEOC/BE4pxzCc3MeH3WWm5PyyQ7J5dxZ/Tk0qM6UqNGfFshkTyROOdcgtq0K5tbXpvLu5nr6Z/clAfO60fnlg3jHdYPeCJxzrkE9Nbcddw6OYNd+3K4+dTu/OKYziQlUCskkicS55xLIFt372dcWiZTZq+lb7smPHheP1JaN4p3WMXyROKccwnig3nrufm1uWzbs5/rhnTjquO7UCspYYdE/JYnEueci7Ptew8wfso8Xpm5mu4HN+Kflw2kZ5vG8Q4rap5InHMujj5etJGxr8xhw85sfnVCV351Qgq1ayZ+KySSJxLnnIuDXdk53PPmfP7z1Uq6tmrIqxcNoF/7pvEO60fxROKccxXs8yWbuHHSHNZs28uVx3bm2iHdqFsrKd5h/WieSJxzroLs2Z/D/e8s5NnPl9OxeX0mXXUkAzo0i3dYZeaJxDnnKkD68i1cP3E2yzfv4dKjOnLT0O7Uq115WyGRPJE451wM7TuQy5/fX8ST/1tK26b1+M8VgzmyS/N4h1WuPJE451yMzFq1jetensWSjbu5cFAyvz+tBw3rVL2P3apXI+eci7PsnFwe+e9invh4Ka0a1eG5nw/k2G4t4x1WzMT0YmVJQyUtlJQlaWwh2+tIeincPk1Sx3D9QEmzwsdsSWdHe0znnIunzLXbGf7oZzw2dQlnH96Wd357bJVOIhDDFomkJOAxYAiwGpguKc3M5kUUuxzYamZdJY0E7gPOBzKA1HDe90OA2ZKmABbFMZ1zrsIdyM3jb1OX8NcPF3NQg9o8NSqVE3u0jndYFSKWp7YGAllmthRA0ovAcCDyQ384cEe4PAl4VJLMbE9EmboECSTaYzrnXIVatH4n1708m7lrtjP8sDbcMawXBzWoHe+wKkwsE0lbYFXE89XAoKLKhK2P7UBzYJOkQcDTQAfg4nB7NMcEQNJoYDRAcnJy2WvjnHMF5OYZEz5ZykPvL6JR3Zo8/rP+nNrnkHiHVeEStrPdzKYBvST1AP4p6e1S7j8BmACQmppqJRR3zrlSWbpxF9dNnM3XK7dxau+Dueus3rRoWCfeYcVFLBPJGqB9xPN24brCyqyWVBNoAmyOLGBm8yXtAnpHeUznnIuZvDzjmc+Xc/87C6hbK4mHRx7Gmf3aICXmpFMVIZaJZDqQIqkTwYf9SODCAmXSgFHAF8C5wIdmZuE+q8LTWR2A7sByYFsUx3TOuZhYuXkP10+azVfLtnBC91b88Zw+tG5cN95hxV3MEkmYBMYA7wJJwNNmlilpPJBuZmnAU8DzkrKALQSJAeBoYKykA0AecI2ZbQIo7JixqoNzzgGYGf+atpI/vjWfJIn7z+3LeQPaVetWSCSZVf3ug9TUVEtPT493GM65SmjNtr3cNGkOn2Zt4piUFtw3oi9tmtaLd1gVQtIMM0stqVzCdrY751w8mRkT01dz1xvzyDXjnrN7c+HAZG+FFMITiXPOFbB+xz7GvjKHqQs3MqhTMx44rx/tm9WPd1gJyxOJc86FzIzXZ63l9rRMsnNyuX1YT0Yd2ZEaNbwVUhxPJM45B2zcmc2tk+fybuZ6+ic35YHz+tG5ZcN4h1UpeCJxzlV7b85Zx22vZ7ArO4ebT+3OL47pTJK3QqLmicQ5V21t3b2f217P4I056+jbrgkPntePlNaN4h1WpeOJxDlXLb0/bz03vzqX7Xv3c92Qblx9fBdqJsV0Zo0qyxOJc65a2b73AHdOyeTVmWvocUhjnvv5QHq2aRzvsCo1TyTOuWrj40UbuWnSHDbuyuZXJ3TlVyekULumt0LKyhOJc67K25Wdwz1vzuM/X62ia6uG/P3iAfRr3zTeYVUZnkicc1Xa50s2ccPEOazdvpcrj+3MtUO6UbdWUrzDqlI8kTjnqqQ9+3O47+0F/POLFXRq0YBJVx3JgA7N4h1WleSJxDlX5aQv38L1E2ezfPMeLj2qIzcN7U692t4KiZViE4mkdgRDux8DtAH2AhnAm8DbZpYX8widcy5K+w7k8uB7C/nHp8tod1A9Xhw9mMGdm8c7rCqvyEQi6RmCOdXfAO4DNgB1gW7AUOAWSWPN7JOKCNQ554oza9U2rnt5Fks27uZng5L5/Wk9aFDHT7pUhOJ+yw+aWUYh6zOAVyXVBpJjE5ZzzkUnOyeXR/67mMc/WkLrxnV5/vKBHJPSMt5hVSvFJZJTJW0zs9WFbTSz/UBWbMJyzrmSZazZzvUTZ7Pgm52cN6Adtw3rSeO6teIdVrVT3J04bYAvJP1P0jWSSp3iJQ2VtFBSlqSxhWyvI+mlcPs0SR3D9UMkzZA0N/x5QsQ+F4Tr50h6R1KL0sblnKvcDuTm8ZcPFnHWY5+xefd+nhqVyp/O6+dJJE6KTCRmdi3BqatbgT5A/gf3KEkljmomKQl4DDgV6AlcIKlngWKXA1vNrCvwEEFfDMAmYJiZ9QFGAc+Hx6wJPAz81Mz6AnOAMdFW1jlX+S38Zidn/+0z/vLBYk7vewjvX3ssJ/ZoHe+wqrVixwawwMdmdjXQjuDD/rfA+iiOPRDIMrOl4WmwF4HhBcoMB/4ZLk8CTpQkM/vazNaG6zOBepLqAAofDRTMd9kYWItzrsrLyc3jbx9lMeyvn7Ju2z6euKg/D488nKb1a8c7tGovqksaJPUhuAz4fILWws1R7NYWWBXxfDUwqKgyZpYjaTvQPHyNfCOAmWaWHcZyNTAX2A0sBn5ZRMyjgdEAycl+TYBzldmSjbu47uXZzFq1jVN7H8zdZ/WmecM68Q7LhYq7/DeFIHmMBHIJWhQnm9nSCooNSb0ITnedHD6vBVwNHA4sBf5KkNTuLrivmU0AJgCkpqZaBYXsnCtHeXnG058t40/vLqRurSQeHnkYZ/ZrQ3BCwiWK4lok7wD/Ac4v4jLgkqwB2kc8bxeuK6zM6rD/owmwGb69GfI14BIzWxKWPwwg/7mkl4EfdOI75yq/FZt3c8PEOXy1fAsndm/FH8/pQ6vGdeMdlitEkYnEzLrkL0vqAKSY2QeS6gE1zWxnCceeDqRI6kSQMEYCFxYok0bQmf4FcC7woZmZpKYEd8+PNbPPIsqvAXpKamlmG4EhwPxoKuqcqxzy8ox/T1vBH99eQJLEn87ty7kD2nkrJIGV2Eci6QqCvoZmQBeClsUTwInF7Rf2eYwB3gWSgKfNLFPSeCDdzNKAp4DnJWUBWwiSDQRXYnUFxkkaF6472czWSroT+ETSAWAFcGlpKuycS1xrtu3lxkmz+SxrM8ektOC+EX1p07RevMNyJZBZ8d0HkmYRXIE1zcwOD9fNDS/NrRRSU1MtPT093mE454pgZrycvoq73phPnhm3nN6DCwcmeyskziTNMLPUkspFc9VWtpntz39Dw74M77x2zpWL9Tv2MfaVOUxduJHBnZvxp3P70b5Z/XiH5UohmkTysaTfE9zLMQS4BpgS27Ccc1WdmTF51hpufz2T/bl53DGsJ5cc2ZEaNbwVUtlEk0jGEtyBPhe4EngL+Ecsg3LOVW0bd2Zzy2tzeW/eegZ0OIgHzutHpxYN4h2W+5FKTCThnCNPhg/nnCuTN+es49bJc9m9P5ffn9ady4/uTJK3Qiq14m5InEJwQ987ZnagwLbOBFdLLTezp2MaoXOuStiyez/jXs/gjTnr6NeuCQ+c14+U1iUO2+cqgeJaJFcAvwP+ImkLsJFgYqtOBMPHP2pmr8c+ROdcZfde5jf8/rUMtu/dz/Und+Oq47pQM6nYof5cJVLcDYnfADcCN4bDux9CMNXuIjPbUyHROecqte17DnDnlExe/XoNPQ5pzHM/H0jPNo3jHZYrZ1EN2mhmy4HlMY3EOVelfLRwAze9ModNu/bz6xO6MuaEFGrX9FZIVeQTGjvnytXOfQe45835vDh9FSmtGvLkJan0bdc03mG5GPJE4pwrN59nbeKGSXNYt30vVx7XmWtP6kbdWknxDsvFWDRjbQ0D3gwvA3bOuR/Ysz+He99ewHNfrKBTiwZMvOooBnQ4KN5huQoSTYvkfIIrt14hGHhxQYxjcs5VItOWbubGV+awYvMeLvtJR248pTv1ansrpDqJ5obEiyQ1Bi4AnpVkwDPAf6IYSt45V0Xtzs7h/ncW8M8vVpDcrD4vjh7M4M7N4x2Wi4Nor9raIWkSUI9gzvazgRskPWJmf41lgM65xPP5kk3c9MocVm/dy6VHdeTGoYdSv7Z3uVZX0fSRnAlcRjA/yHPAQDPbIKk+MI9gulvnXDWwKzuHe9+ez7++XEnH5vV5afSRDOzULN5huTiL5ivECOAhM/skcqWZ7ZF0eWzCcs4lmk8XB62Qtdv3cvnRnbj+5EO9L8QB0fWRjCpm23/LNxznXKLZue8Af3hrAf/5aiWdWzRg0lVHMqCDt0Lcd0q8zVTSTkk7CjxWSXotHLyxuH2HSlooKUvS2EK215H0Urh9WjgUC5KGSJohaW7484SIfWpLmiBpkaQFkkaUvtrOuWh8smgjpzz0CS9NX8mVx3bmrd8c40nE/UA0p7b+AqwGXgBEMK96F2Am8DRwfGE7SUoCHgOGhPtPl5RmZvMiil0ObDWzrpJGAvcRXG68CRgWztHem2De97bhPrcAG8ysm6QaBHPJO+fK0Y59B7jnjfm8lL6KLi0bMOnqo+if7PeFuMJFk0jONLN+Ec8nSJplZjeFMycWZSCQZWZLASS9CAwn6KDPNxy4I1yeBDwqSWb2dUSZTILZGeuYWTbwc6A7fDtXyqYo6uCci9LUBRu4+dW5bNi5j6uP78JvTkzxu9NdsaIZQW2PpP+TVCN8/B+wL9xW3NztbYFVEc9X812r4gdlzCwH2A4UvBB9BDDTzLIl5Q/Yc5ekmZImSmpd2ItLGi0pXVL6xo0bS6ykc9Xd9j0HuH7ibC57djqN69XktWt+wk1Du3sScSWKJpH8DLgY2ACsD5cvklQPGBPD2JDUi+B015XhqppAO+BzM+sPfAE8UNi+ZjbBzFLNLLVly5axDNO5Su+DeesZ8tDHvPb1Gsb8tCtTfnU0/dr7QIsuOsWe2gr7OYab2bAiinxazO5rgPYRz9uF6wors1pSTaAJsDl87XbAa8AlZrYkLL8Z2AO8Gj6fSNDP4pz7Ebbt2c/4KfN49es1dD+4EU+NOoI+7ZrEOyxXyRTbIjGzXIKhUX6M6UCKpE6SahN00qcVKJMG5F9efC7woZlZeArrTWCsmX0WEY8BU/iug/9Evt/n4pyL0nuZ3zDkoU9Im72WX5+YQtqYoz2JuB8lms72zyQ9CrwE7M5faWYzi9vJzHIkjSG44iqJYMDHTEnjgXQzSwOeAp6XlAVsIUg2EJwy6wqMkzQuXHeymW0Abgr3+QvB9L+XRVlX5xzB3Ol3pGWSNnstPQ5pzDOXHkHvtp5A3I+n4Et+MQWkqYWsNjM7oZD1CSk1NdXS09PjHYZzcfdOxjpunZzBtj0H+NUJKVx9fBeftdAVSdIMM0stqVw0d7b/tHxCcs7Fy+Zd2dyelskbc9bRu21jnr98ED0O8bnTXfmIZtDG1sAfgDZmdqqknsCRZvZUzKNzzpXZm3PWMe71DHbsO8D1J3fjyuO6UCvJWyGu/ETTR/Iswfwjt4TPFxH0l3gicS6BbdqVzbjXM3hr7jf0aduEF84bzKEHN4p3WK4KiiaRtDCzlyXdDN92oufGOC7n3I9kZkyZs47bX89gd3YuNw49lNHHdKamt0JcjESTSHZLak54F7ukwQR3oDvnEsyGnfu4bXIG72aup1/7pjxwbl9SWnsrxMVWNInkdwT3e3SR9BnQkuCeD+dcgjAzXp+1ljumZLJnfy43n9qdy4/u5K0QVyGiuWprpqTjgEMJRv9daGYHYh6Zcy4q63fs45bX5vLB/A30T27K/ef2o2urhvEOy1Uj0U6yPBDoGJbvLwkzey5mUTnnSmRmvDpzDXdOySQ7J49bT+/BZT/pRFINxTs0V81Ec/nv8wTzj8wC8jvZjWD+dudcHHyzfR83vzqHqQs3ktrhIO4/ty+dW3orxMVHNC2SVKCnlXQLvHMu5syMiTNWc9cb8ziQm8e4M3oy6qiO3gpxcRVNIskADgbWxTgW51wx1m7by9hX5/LJoo0M7NSM+0f0pWOLBvEOy7no7iMB5kn6CsjOX2lmZ8YsKufct8yMF6ev4p4355ObZ9x5Zi8uHtyBGt4KcQkimkRyR6yDcM4VbvXWPdz86lz+t3gTR3Zuzn0j+pLcvH68w3Lue6K5/PdjSR2AFDP7QFJ9gmHhnXMxYma88NVK/vDmfADuPqs3Fw5M9laIS0jRXLV1BTAaaEZw9VZb4AmCSaWcc+Vs1ZY93PTKHD5fspmfdG3Ovef0pX0zb4W4xBXNqa1fEtxHMg3AzBZLahXTqJyrhvLyjH9NW8G9by+ghsQfzu7DBQPbI3krxCW2aBJJtpntz/9jDudW90uBnStHKzbv5sZJc5i2bAvHpLTg3hF9adu0XrzDci4q0QzE87Gk3wP1JA0BJhLMm14iSUMlLZSUJWlsIdvrSHop3D5NUsdw/RBJMyTNDX/+YDZGSWmSMqKJw7lElZdnPPvZMob+5X/MW7uD+0f05bmfD/Qk4iqVaFokY4HLgbnAlcBbwD9K2klSEvAYMARYDUyXlGZm8yKKXQ5sNbOukkYC9wHnA5uAYWa2VlJvgnnf20Yc+xxgVxSxO5ewlm3azU2T5vDV8i0cf2hL/nhOHw5p4gnEVT7RXLWVBzwZPkpjIJBlZksBJL0IDAciE8lwvru8eBLwqCSZ2dcRZTIJWkN1zCxbUkOCEYlHAy+XMibn4i43z3jms2U88N5CaiXV4IHz+jGif1vvC3GVVrSDNv4YbYFVEc9XA4OKKhNOmLUdaE7QIsk3AphpZvk3Q94FPAjsiUXQzsXSko27uHHSHGas2MqJ3Vtxz9l9OLhJ3XiH5VyZxDKRlJmkXgSnu04Onx8GdDGza/P7U4rZdzRBq4Xk5OTYBupcCXLzjKc+XcqD7y2ibq0kHjq/H2cd5q0QVzXEMpGsAdpHPG8XriuszOrwarAmwGYASe2A14BLzGxJWP5IIFXScoLYW0n6yMyOL/jiZjYBmACQmprqV5m5uMnasJMbJs3h65XbGNKzNfec1ZtWjb0V4qqOaG5I7AbcAHSILG9mP7iSqoDpQIqkTgQJYyRwYYEyacAo4AuCWRc/NDOT1BR4ExhrZp9FvObjwONhXB2BNwpLIs4lgpzcPJ783zIe+mAR9Wsn8fDIwzizXxtvhbgqJ5oWyUSCO9mf5Lv5SEoU9nmMIbj96F/bAAAa8ElEQVTiKgl42swyJY0H0s0sDXgKeF5SFrCFINkAjAG6AuMkjQvXnWxmG6J9fefiadH6ndwwcTazV29naK+Dueus3rRsVCfeYTkXEyppmhFJM8xsQAXFExOpqamWnp4e7zBcNZCTm8ffP1nKwx8spmHdmowf3ovT+xzirRBXKYWf/6kllYumRTJF0jUE/RWRw8hvKUN8zlU5C77ZwQ0T5zB3zXZO73sI48/sRfOG3gpxVV80iWRU+POGiHUGdC7/cJyrfA7k5vH4R0v464eLaVy3Fn/7WX9O63NIvMNyrsJEc0Nip4oIxLnKKHPtdm6YOId563YwrF8b7jyzF80a1I53WM5VqGiu2qoFXA0cG676CPi7mR2IYVzOJbT9OXk8NjWLx6Zm0bR+bZ64aABDex8c77Cci4toTm09DtQC/hY+vzhc94tYBeVcIstYs53rJ85mwTc7Ofvwtow7oycHeSvEVWPRJJIjzKxfxPMPJc2OVUDOJarsnFwe/TCLv320hOYNavPkJakM6dk63mE5F3fRJJJcSV3y7y6X1JlS3E/iXFUwZ/U2bpg4h4XrdzKifzvGndGTJvVrxTss5xJCNInkBmCqpKWACO5wvyymUTmXILJzcnn4g8X8/ZOltGhYm6cvTeWE7t4KcS5SNFdt/VdSCnBouGphxEi8zlVZs1Zt44aJs1m8YRf/l9qOW07vSZN63gpxrqAiE4mkE8zsw3ASqUhdJWFmr8Y4NufiYt+BXB76YBFPfrKU1o3r8uxlR3D8oa3iHZZzCau4FslxwIfAsEK2GeCJxFU5M1Zs5cZJs1mycTcXDGzPzaf1oHFdb4U4V5wiE4mZ3R4ujjezZZHbwhF9nasy9h3I5cH3FvKPT5fRpkk9nr98IMektIx3WM5VCtF0tr8C9C+wbhJQqQdydC7f9OVbuHHSHJZt2s3PBiVz82k9aFgnoed8cy6hFNdH0h3oBTQp0E/SGPBZeVylt3d/Ln96dyHPfL6Mtk3r8cIvBnFU1xbxDsu5Sqe4r12HAmcATfl+P8lO4IpYBuVcrE1bupkbX5nDis17uOTIDtw0tDsNvBXi3I9SXB/J68Drko40sy8qMCbnYmZ3dg73v7OAf36xguRm9fnPFYM5skvzeIflXKUWzVewryX9kuA017entMzs5zGLyrkY+HzJJm56ZQ6rtuzl0qM6cuPQQ6lf21shzpVVNP9FzwMLgFOA8cDPgPmxDMq58rQrO4d7357Pv75cScfm9Xn5yiMZ2KlZvMNyrsqoEUWZrmZ2G7DbzP4JnA4MiubgkoZKWigpS9LYQrbXkfRSuH2apI7h+iGSZkiaG/48IVxfX9KbkhZIypR0b7QVddXTZ1mbOOWhT/j3tJVcfnQn3v7NsZ5EnCtn0bRI8ucd2SapN/ANUOJtvpKSgMeAIcBqYLqkNDObF1HscmCrmXWVNBK4Dzgf2AQMM7O14Wu+C7QN93nAzKZKqg38V9KpZvZ2FPVw1cjOfQf449sLeGHaSjq3aMCkq45kQAdPIM7FQjSJZIKkg4DbgDSgYbhckoFAlpktBZD0IjAciEwkw4E7wuVJwKOSZGZfR5TJBOpJqmNme4CpAGa2X9JMoF0Usbhq5JNFGxn7yhy+2bGP0cd25ndDulG3VlK8w3KuyoomkTxjZrnAx5Runva2wKqI56v54Smxb8uYWY6k7UBzghZJvhHAzIIDRUrKvyz54cJeXNJoYDRAcnJyKcJ2ldWOfQe45435vJS+ii4tGzDp6qPon3xQvMNyrsqLJpEsk/QO8BLwoZlZjGP6lqReBKe7Ti6wvibwH+CR/BZPQWY2AZgAkJqaWmExu/iYunADv391Lut37OOq47rw25NSvBXiXAWJJpF0J7gx8ZfA05KmAC+a2acl7LcGaB/xvF24rrAyq8Pk0ATYDCCpHfAacEn+pFoRJgCLzewvUcTvqrDtew9w9xvzmDhjNSmtGvL4NT/hsPZN4x2Wc9VKNPOR7AFeBl4O+0oeJjjNVdLXvelASjjA4xpgJHBhgTJpwCjgC+BcwhZPeNrqTWCsmX0WuYOkuwkSjs8ZX839d/56fv/aXDbt2s+Yn3blVyd2pU5Nb4U4V9GiuhtL0nEEV1MNBdKB/ytpn7DPYwzBFVdJwNNmlilpPJBuZmnAU8DzkrKALQTJBmAM0BUYJ2lcuO5koDZwC8F9LTMlATxqZv+Iph6uati2Zz/jp8zj1a/XcGjrRvzjkiPo065JvMNyrtpSSV0ekpYDXxO0StLMbHcFxFWuUlNTLT09Pd5huDIyM96a+w13TMlk6+79XHN8F8ackELtmtHcDuWcKy1JM8wstaRy0bRI+prZjnKIybkfbfmm3YxLy+STRRvpeUhjnrn0CHq39VaIc4mguGHkbzSz+4F7JP2g2WJmv45pZM4RTDj1+EdLePzjJdRJqsEdw3py0eAO1EzyVohziaK4Fkn+eFp+TsjFxUcLN3B7WiYrNu9h+GFtuOW0HrRq7FPhOJdoihtGfkq4uMfMJkZuk3ReTKNy1dq67XsZP2Ueb2d8Q+eWDfj3LwbxE59wyrmEFU0fyc3AxCjWOVcmB3LzePaz5Tz0wSJy84wbTjmUXxzTyS/pdS7BFddHcipwGtBW0iMRmxoDObEOzFUv05dv4dbXMli4ficndm/FHWf2on2z+vEOyzkXheJaJGsJ+kfOBGZErN8JXBvLoFz1sXlXNve+vYCJM1bTtmk9Jlw8gCE9WxPeI+ScqwSK6yOZDcyW9BrBXCS58O3w8HUqKD5XReXlGS9OX8V97yxgd3YOVx/fhV+d0NVnLHSuEormv/Y94CRgV/i8XrjuqFgF5aq2jDXbuWVyBrNXbWNw52bcNbw3Ka0bxTss59yPFE0iqWtm+UkEM9slyU9eu1Lbse8Af35vEc99sZxmDWrz0Pn9OOuwtn4ay7lKLppEsltSfzObCSBpALA3tmG5qsTMSJu9lrvfnM+mXdlcPLgD1518KE3q1Yp3aM65chBNIvktMFHSWkDAwQQDODpXoqwNuxj3egafL9lM33ZNeHqUD7DoXFUTzTDy0yV1Bw4NVy00swPF7ePc3v25PDp1MRM+WUq9WkncfVZvLhiYTFINP43lXFUT7SUyhwI9gbpAf0mY2XOxC8tVZv+dv57b0zJZvXUv5/Rvy82n9qBlI7/Qz7mqqsREIul24HiCRPIWcCrwKeCJxH3P6q17uHPKPN6ft56UVg15cfRgBnduHu+wnHMxFk2L5FygH/C1mV0mqTXwr9iG5SqT/Tl5/OPTpTzy38UIcfOp3fn50Z2o5SP0OlctRJNI9ppZnqQcSY2BDXx/LnZXjX2xZDO3vZ5B1oZdnNKrNeOG9aJt03rxDss5V4Gi+cqYHs6h/iTBUCkzCeZYL5GkoZIWSsqSNLaQ7XUkvRRunyapY7h+iKQZkuaGP0+I2GdAuD5L0iPymxDiYsPOfVz70iwuePJLsnNyefrSVP5+caonEeeqoWiu2romXHxC0jtAYzObU9J+4VAqjwFDgNXAdElpZjYvotjlwFYz6yppJHAfwaXFm4BhZrZWUm+Ced/bhvs8DlwBTCPosxkKvF1yVV15yM0z/j1tBX96dyHZB/L49QldueanXalby0foda66KtXARma2vBTFBwJZZrYUQNKLwHAgMpEMB+4IlycBj0qSmX0dUSYTqCepDtCMIJF9GR7zOeAsPJFUiNmrtnHr5AzmrtnO0V1bcOfwXnRp2TDeYTnn4iyWI+S1BVZFPF8NDCqqjJnlSNoONCdokeQbAcw0s2xJbcPjRB6zLYWQNBoYDZCcnFyGarjtew7wp/cW8O9pK2nZsA5/veBwzuh7iA9t4pwDYptIykxSL4LTXSeXdl8zmwBMAEhNTf3BnPOuZGbGqzPX8Ie35rN1z34uO6oT1w5JoVFdH9rEOfedWCaSNXz/6q524brCyqyWVBNoAmwGkNQOeA24xMyWRJRvV8IxXTlYtH4nt07O4KtlW+if3JTnLh9IrzY+tIlz7odimUimAymSOhF82I8ELixQJg0YRXAV2LnAh2Zm4VVibwJjzeyz/MJmtk7SDkmDCTrbLwH+GsM6VDu7s3N45L+LeerTZTSsW5P7RvThvAHtqeFDmzjnihCzRBL2eYwhuOIqCXjazDIljQfSzSwNeAp4XlIWsIUg2QCMAboC4ySNC9edbGYbgGuAZwnmRXkb72gvF2bGu5nrGT8lk7Xb93F+antuOrU7zRrUjndozrkEJ7Oq332Qmppq6enp8Q4jYa3cvIfb0zKYunAj3Q9uxD1n92ZAh2bxDss5F2eSZphZaknlErqz3cVWdk4uEz5eyqNTs6hZQ9x6eg8uPaojNX1oE+dcKXgiqab+t3gj417PZNmm3Zze9xBuO70nBzepG++wnHOVkCeSamb9jn3c9cY83pizjo7N6/PczwdybLeW8Q7LOVeJeSKpJnJy83juixX8+f1F7M/N49qTunHlcZ19aBPnXJl5IqkGZqzYyq2TM5i/bgfHdWvJ+OG96NC8QbzDcs5VEZ5IqrCtu/dz3zsLeHH6Kg5pUpcnLurPKb0O9qFNnHPlyhNJFZSXZ0ycsYp7317Azn05XHlsZ359YgoN6vjb7Zwrf/7JUsXMX7eDWydnMGPFVo7oeBB3n9WHQw9uFO+wnHNVmCeSKmJXdg4Pvb+IZz9fTpN6tXjgvH6M6N/WT2M552LOE0klZ2a8OXcdd70xjw07s7lwYDI3nHIoTev70CbOuYrhiaQSW7ZpN+Nez+B/izfRq01jnrhoAIcnHxTvsJxz1Ywnkkpo34Fc/vbREp74aAl1atbgzjN7cdHgDiT5CL3OuTjwRFLJTF24gdtfz2Tllj2cdVgbfn96D1o18qFNnHPx44mkkli7bS93vTGPtzO+oUvLBrxwxSCO6tIi3mE555wnkkR3IDePZz5bxl8+WEyeGTeccihXHNOZ2jV9hF7nXGLwRJLAvlq2hVsnz2XR+l2c1KM1tw/rSftm9eMdlnPOfY8nkgS0eVc2f3hrAa/MXE3bpvV48pJUhvRsHe+wnHOuUDE9PyJpqKSFkrIkjS1kex1JL4Xbp0nqGK5vLmmqpF2SHi2wzwWS5kqaI+kdSVWmoyAvz/j3tBWc8ODHpM1ewzXHd+H93x3rScQ5l9Bi1iKRlAQ8BgwBVgPTJaWZ2byIYpcDW82sq6SRwH3A+cA+4Dagd/jIP2ZN4GGgp5ltknQ/wfzud8SqHhUlY812bpmcwexV2ziyc3PuOqsXXVv50CbOucQXy1NbA4EsM1sKIOlFYDgQmUiG810SmAQ8Kklmthv4VFLXAsdU+GggaTPQGMiKXRVib8e+Azz47kKe/3IFzRrU4eGRh3FmvzY+tIlzrtKIZSJpC6yKeL4aGFRUGTPLkbQdaA5sKuyAZnZA0tXAXGA3sBj4ZWFlJY0GRgMkJyf/+FrEiJmRNnstd70xny27s7nkyI5cO6QbTerVindozjlXKpWqs11SLeBq4HBgKfBX4Gbg7oJlzWwCMAEgNTXVKjDMEmVt2MVtkzP4Yulm+rVrwjOXHkGfdk3iHZZzzv0osUwka4D2Ec/bhesKK7M67P9oAmwu5piHAZjZEgBJLwM/6MRPVHv35/Lo1MVM+GQp9Wolcc/ZvRl5RLIPbeKcq9RimUimAymSOhEkjJHAhQXKpAGjgC+Ac4EPzay41sMaoKeklma2kaAjf365Rx4DH8xbz+1pmazZtpcR/dtx82ndadGwTrzDcs65MotZIgn7PMYA7wJJwNNmlilpPJBuZmnAU8DzkrKALQTJBgBJywk602tLOgs42czmSboT+ETSAWAFcGms6lAeVm3Zw51T5vHB/PV0a92Ql0YPZlDn5vEOyznnyo2KbwBUDampqZaenl6hr7k/J48n/7eUv364mBoSvz0phct+0olaST60iXOucpA0w8xSSypXqTrbK4vPl2zitskZLNm4m6G9DmbcsJ60aVov3mE551xMeCIpRxt27uMPb85n8qy1JDerzzOXHsFPu7eKd1jOORdTnkjKQW6e8a8vV/DAuwvJzsnj1yemcM3xXahbKyneoTnnXMx5IimjWau2cevkuWSs2cExKS0YP7w3nVo0iHdYzjlXYTyR/Ejb9xzg/ncX8MJXK2nVqA6PXdif0/oc7EObOOeqHU8kpWRmvDJzDX98az7b9h7g5z/pxG9PSqFRXR/axDlXPXkiKYWF3+zktskZfLV8C/2Tm/L8WX3o2aZxvMNyzrm48kQShd3ZOTzy38U89ekyGtWtyf0j+nLugHbU8KFNnHPOE0lxzIx3M7/hzinzWLd9HyOPaM9NQ7tzUIPa8Q7NOecShieSIuTk5jH6+Rl8uGADPQ5pzKMX9mdAh4PiHZZzziUcTyRFqJlUg47NG3DbGT0ZdWQHavrQJs45VyhPJMUYN6xnvENwzrmE51+znXPOlYknEuecc2XiicQ551yZeCJxzjlXJp5InHPOlYknEuecc2XiicQ551yZeCJxzjlXJjKzeMcQc5I2AiviHUcUWgCb4h1EOfG6JJ6qUg/wulSUDmbWsqRC1SKRVBaS0s0sNd5xlAevS+KpKvUAr0ui8VNbzjnnysQTiXPOuTLxRJJYJsQ7gHLkdUk8VaUe4HVJKN5H4pxzrky8ReKcc65MPJE455wrE08kFUjSbyRlSMqU9NtwXTNJ70taHP48KFwvSY9IypI0R1L/+Eb/fUXU5Q5JayTNCh+nRZS/OazLQkmnxC9ykPS0pA2SMiLWlfp9kDQqLL9Y0qhKUJfjJW2PeH/GRewzNHxvsiSNTaC6nBf+jeVJSi1QvtC/qXjXpTT1kNRR0t6I9+SJiG0DJM0N6/GIJFV0XaJmZv6ogAfQG8gA6hPMTPkB0BW4HxgblhkL3Bcunwa8DQgYDEyLdx2iqMsdwPWFlO8JzAbqAJ2AJUBSHOM/FugPZESsK9X7ADQDloY/DwqXD0rwuhwPvFHIMZLC96QzUDt8r3omSF16AIcCHwGpJf1NJUJdSlmPjpHlChznq/BvTuHf4KkV/Z5E+/AWScXpQfAhtMfMcoCPgXOA4cA/wzL/BM4Kl4cDz1ngS6CppEMqOugiFFWXogwHXjSzbDNbBmQBAysgzkKZ2SfAlgKrS/s+nAK8b2ZbzGwr8D4wNPbRf18p61KUgUCWmS01s/3Ai+ExKlRhdTGz+Wa2sJDiRf1Nxb0upaxHocK/scZm9qUFWeU5Sn4f48YTScXJAI6R1FxSfYJvuu2B1ma2LizzDdA6XG4LrIrYf3W4LhEUVReAMeEpoKfzT6mQ2HXJV9r3IZHrVFRdAI6UNFvS25J6hesSuS5FqYzvS1E6Sfpa0seSjgnXtSWIPV9C18MTSQUxs/nAfcB7wDvALCC3QBkDEv567GLq8jjQBTgMWAc8GK8Yy6KyvA/RKFCXmQRjJ/UD/gpMjltgLt86INnMDgd+B7wgqXGcYyo1TyQVyMyeMrMBZnYssBVYBKzPP2UV/twQFl/Dd9/yAdqF6xJCYXUxs/VmlmtmecCTfHf6KqHrEirt+5DIdSq0Lma2w8x2hctvAbUktSCx61KUyvi+/EB4am5zuDyDoH+nG0HM7SKKJnQ9PJFUIEmtwp/JBH0KLwBpQP4VP6OA18PlNOCS8KqhwcD2iNMVcVdYXQr04ZxNcAoMgrqMlFRHUicghaAjMZGU9n14FzhZ0kHhKbyTw3WJoNC6SDo4/8ofSQMJ/v83A9OBFEmdJNUGRobHSGRF/U1VqrpIaikpKVzuTFCPpeHf2A5Jg8P37BK++5tMPPHu7a9OD+B/wDyCK0lODNc1B/4LLCa4+qlZuF7AYwTfUOYScaVHIjyKqMvzYaxzCP55D4kof0tYl4XE+eoT4D8EpxQOEJx7vvzHvA/Azwk6ebOAyypBXcYAmeF79iVwVMRxTiNoIS8BbkmgupwdLmcD64F3S/qbinddSlMPYET4nswiOPU4LOI4qQRfxpYAjxKORJKIDx8ixTnnXJn4qS3nnHNl4onEOedcmXgicc45VyaeSJxzzpWJJxLnnHNl4onEVWmSpoWjqq6UtDFilNWO8Y6tKJLOkdQ9hsevIWmqpIaSakraFrFtWDhybntJv5V0SazicFVHzXgH4FwsmdkgAEmXEtwDMia+EQUk1bRgwMvCnAPkAQvK6XgFDQPSzWyXpG8/A8Kh2B8ChpjZKkn/IBiQ87lo43DVk7dIXLUl6VRJX0iaKeklSQ3C9asl/SEc3HC6pP6S3pO0RNIVYZmTwm/1b4ff4B+LuGu8uOPeK+lr4GxJV4XHny1poqR64aB9pwEP5becJH0q6bDwGAdLygqXfyFpsqSphHfVSxor6atw4MxxBesc+hkF7pKW9FPgbwQ39i0DsGA4lbVKsLlwXOLxROKqpXCIl7EEd+X3J7gb/zcRRZZZMLjhl8BTBHcmHwXcFVFmEHA1wdwYPYDhURx3g5kdbmYTgYlmdkT4OkuAS83sf8BbwLVmdpiZLS+hKocD55jZiQomEksO4zoMOErSUYXscxTBXdT56gOvAMPNbHGBsunAMThXDD+15aqrowgSwOdhQ6I28GnE9vzxmeYCNc1sN7BbwQx3DcNtX+Z/0Et6ETg6XF/ccV+KWO4raTzQFGgEvPEj6vGeBfOhQDDe16nA1+HzhgQDAH5eYJ/GZrYn4vk+YBpwGXBdgbIbCCZfcq5InkhcdSXgHTO7uIjt2eHPvIjl/Of5/zcFxxeyKI67O2L5OYJTSRmSfkEwG15hcvju7EHdYo4n4G4ze6qI4+TLK+T5ucCHkm40s/sjttUF9pZwPFfN+aktV119DhwXjriKpAaSUkp5jMGSksPRW/+PoOVRmuM2AL6RVAu4MGL9ToIWSr7lwIBw+dxi4nkXuDyiT6adgmHiC8oqeNVa2OI6HbhM359/vhvfjeLsXKE8kbhqyczWE4zK+pKk2QQJoFspD/MV8ATBKMgLgbRSHnccwbDnn4XHyPcf4PcRlyn/CfiNpJkE88MXVae3gEnAl5LmAi8TnN4q6E2C+dsL7r+JYLrgOyWdHq4+kmAkYeeK5KP/OvcjSDoJGGNmCTuPdlEktQP+YWbFzjEv6QjgGjO7rGIic5WVt0icq2bMbDXwbMRFA0VpBtxeASG5Ss5bJM4558rEWyTOOefKxBOJc865MvFE4pxzrkw8kTjnnCsTTyTOOefK5P8Bie3Wx8q2Sl8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEWCAYAAABbgYH9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecFfXVx/HPd3eBBekdFpBqWVSKiIIlJhZARYoNTBQVu0RJTPJgNGp8kjxRY7CioqhYIiAKYkWNBVSkdxBZAYWlFyki/Tx/zKxe1y2zuHfv3d3zfr3ui7kzv5l7xrvu2Zn5nd9PZoZzzjlX3FISHYBzzrmyyROMc865uPAE45xzLi48wTjnnIsLTzDOOefiwhOMc865uPAE48ocSR9KujJcvkzSxxH3e0bS3+IbHUh6TNJf4nTshZJOjcexnSuqtEQH4MoHSSuABsB+4FvgLWCQme1IZFzxJuky4EozOylnnZldW0zHfgZYZWa3xRy7bXEc27ni4FcwriT1NLOqQEegE3BbIe1/QpL/UVSG+fdbtniCcSXOzLIJrmCOApB0uaTFkrZLWibpmpy2kk6VtErS/0haCzwtqZak1yVtkLQlXG4S5bMlHSHpXUmbJS2RdGHE/VpJel/SJkkbJb0gqWbM9qaSXglj2iTpYUlHAo8BXSTtkPRN2Pb7W3HheZ8Tc5y08Bgdw/cvSVoraaukSZLahuuvBn4N/Ck89mvh+hWSTg+XK0m6X9Lq8HW/pEq5/rveLGm9pDWSLi/g/GtIGhG2y5b0N0mp4bbLJH0s6V/h97FcUo8i7PuJpKGSNgF3SkqVdF/433m5pEGSLPxvc4Gkmbli+72kV6N8j65keYJxJU5SU+AsYHa4aj1wDlAduBwYmvMLNtQQqA0cClxN8HP7dPi+GfAd8HCEzz0EeBf4D1Af6AcMk5QZJWzg/4DGwJFAU+DO8LipwOvAV0BzIAMYZWaLgWuBKWZW1cxq/vSwvAj0j3nfDdhoZrPC928BbcJ4ZwEvAJjZ8HD5nvDYPfM49q3ACUB7oB3QmR9fNTYEaoTxDgQekVQrn/N/BtgHtAY6AGcCV8ZsPx5YAtQF7gFGSFIR9l1GcAv178BVQI8w7o5A75i2E4AWYfLOcQnwbD5xu0QyM3/5K+4vYAWwA/iG4BfxMKByPm3HAzeFy6cCe4D0Ao7dHtgS8/5DguceAJcBH4fLFwGTc+37OHBHuPwM8LeI59MbmB0udwE2AGl5tPv+82PWff85BL90twNVwvcvALfn85k1AQNq5Bdv+N/59HD5S+CsmG3dgBUx/12/i42ZINGfkMfnNgB2x35fBEnxg5hzzIrZViWMs2HEfb/O9XnvA9fEvD89PF5a+P5R4O/hcltgC1Ap0T/j/vrpy+93upLU28zey70yvJ1yB3AYwdVJFWB+TJMNZrYrpn0VYCjQHcj5i7uapFQz21/A5x8KHJ9zqyqUBjxXWOCSGgAPACcD1cI4t4SbmwJfmdm+wo6Tm5llSVoM9Axvc51L8Fd+zpXR34ELgHrAgXC3usDWCIdvTJDMc3wVrsuxKVfMO4GqeRznUKACsOaHixJSgJUxbdbGnNPOsF1VgivPwvaNXc6Ju6DtI4EXJd1GcPUyxsx25xG3SzBPMC6hwmcCLwOXAq+a2V5J4wluSeXIPeT3zcDhwPFmtlZSe4LbbaJgK4GPzOyMgwj1H2EcR5vZZkm9+eG23EqgmaS0PJJMlOHKc26TpQCLzCwrXH8x0IvgL/gVBLeztvDDeRZ27NUEyWFh+L5ZuK6oVhJchdQ9iCQaZd/c57EGiH2m1vRHjc0+k7SHINlfHL5cEvJnMC7RKgKVCG4x7QuvZs4sZJ9qBLd3vpFUm+DqJ4rXgcMkXSKpQvg6Ltf9/II+cwewVVIG8MeYbdMIfin+U9IhktIlnRhuWwc0kVSxgGOPIjjn6wieD8V+5m5gE8FV3T9y7bcOaFnAcV8EbpNUT1Jd4Hbg+QLa58nM1gDvAPdJqi4pJez08Is47TsGuElSRtiR4n/yaPMsQYLfa2aR6pxcyfME4xLKzLYDNxL8UtlC8NfohEJ2ux+oDGwEPgPeLsJnnUnwcH81wW2duwkSXGH+SvDAeSvwBvBKzHH3Az0Jnqd8DawieN4DwfOEhcBaSRvziWsNMAXoCoyO2fQswW2tbGARwbnGGgFkSvomvOrL7W/ADGAewS3HWeG6g3EpwR8Diwi+p7FAozjt+wRBUppHcGX6JkEngdjbn88R9EIscsJ0JUfhgzLnnEtK4VXtY2Z2aMy6ygSdEjqa2dKEBecK5FcwzrmkIqmypLPCupcMglug43I1uw6Y7sklufkVjHMuqYS9BD8CjiB41vYGQbf1beH2FQQdHXqb2ez8juMSzxOMc865uPBbZM455+KiXNfB1K1b15o3b57oMJxzrlSZOXPmRjOrV1i7cp1gmjdvzowZMxIdhnPOlSqSviq8ld8ic845FyeeYJxzzsWFJxjnnHNxEdcEI6m7gkmdsiQNyWN7JUmjw+1TJTWP2XZLuH6JpG7huqaSPpC0SMHc4zfFtL8gXHdAUqd4npdzzrnCxS3BhEONP0IwcVAm0D+PiZ0GEszj0Zpg+PW7w30zCcaLakswJPuw8Hj7gJvNLJNgIqUbYo65AOgLTIrXOTnnnIsunlcwnQkmIVpmZnsIRoztlatNL4K5HSAYAO+0cBa8XgQzAu42s+VAFtDZzNZYONNfOHDhYoLZ+DCzxWa2JI7n45xzrgjimWAy+PFEQavCdXm2CeeK2ArUibJveDutAzC1KEFJulrSDEkzNmzYUJRdnXPOFUGpfMgvqSrBJFWDc8YnisrMhptZJzPrVK9eoXVCefpwyXpGfLycDdt9Ej3nnMtPPAsts/nxTHRNwnV5tVklKY1gxr5NBe0rqQJBcnnBzF4hAd7/fD3PTvmKf7y5mJNa16VPhwzObNuAKhXLdd2qc879SDx/I04H2khqQZAc+vHTqU0nAAMIJls6H3jfzEzSBOA/kv5NMD93G2Ba+HxmBLDYzP4dx9gLdFevo7jkhEMZPyeb8bNXM3j0HKpUTKV724b07pDBia3rkppS2Oy9zjlXtsV1NGVJZxHMPpgKPGVmf5d0FzDDzCZISieYma4DsBnoZ2bLwn1vBa4g6Dk22MzeknQSMJlgdr4D4cf82czelNQHeAioB3wDzDGzbgXF16lTJ/u5Q8UcOGBMX7GZ8XOyeX3eGrbv2ke9apXo1a4xvTtk0LZxdYK86JxzZYOkmWZWaDlIuR6uvzgSTKxde/fzwefrGTc7mw+WrGfvfqNN/ar07pBB7w4ZZNSsXGyf5ZxzieIJJoLiTjCxvtm5hzfmr2H87Gymr9gCwPEtatOnQwY9jm5EjcoV4vK5zjkXb55gIohngon19aadvDonm3Gzs1m28VsqpqVw2hH16dMhg1MPr0/FtFLZmc85V055gomgpBJMDjNj3qqtjJudzWtzV7Pp2z3UrFKBs49uRJ8OGRx7aC1/XuOcS3qeYCIo6QQTa+/+A3yctZFxs7J5Z9Fadu09QNPalenTPoNeHTJoVa9qQuJyzrnCeIKJIJEJJtaO3fuYuGAt4+dk80nWRg4YtGtSg94dMujZrjF1q1ZKdIjOOfc9TzARJEuCibVu2y4mzFnNuNnZLFqzjdQUcUqbuvTukMGZmQ2pXDE10SE658o5TzARJGOCibVk7XbGz8nm1dnZrN66i0MqptL9qOB5TZdWdbyY0zmXEJ5gIkj2BJPjwAFj6vLNjJ+dzZvz17B99z7qV6tEr/aN6dOhCUc2quadA5xzJcYTTASlJcHE2rV3P/9dHBRzfrhkPfsOGIc3qEbvDhn0at+Yxl7M6ZyLM08wEZTGBBNr87dBMee4WauY9fU3SHBCizr06ZBB96MbUj3dizmdc8XPE0wEpT3BxPpq07eMn72acbNXsWLTTiqmpXDGkQ3o0yGDUw6r58Wczrli4wkmgrKUYHKYGXNWfsP42dm8Nm8Nm7/dQ60qFTjnmGDwzY7NavrzGufcz+IJJoKymGBi7d1/gMlLN/DKrGzeXbSO3fsOcGidKvRuHwy+2aLuIYkO0TlXCnmCiaCsJ5hY23ft5e2wmPPTLzdhBu2b1qRPhwzOOaYRdbyY0zkXkSeYCMpTgom1Zut33xdzfr52O2kp4heH1aN3hwzOyGxAegUv5nTO5c8TTATlNcHEWrxmW1jMuZq123ZRtVIaPY5qSJ8OGRzf0os5nXM/5QkmAk8wP9h/wJi6bBPjZmfz1oK17Ni9j4bV04Nizo4ZHNGweqJDdM4liaRIMJK6Aw8QTJn8pJn9M9f2SsCzwLHAJuAiM1sRbrsFGAjsB240s4mSmobtGwAGDDezB8L2tYHRQHNgBXChmW0pKD5PMHn7bs9+3lu8jvGzs/noiw3sO2Ac0bAafTpkcG77xjSq4cWczpVnCU8wklKBL4AzgFXAdKC/mS2KaXM9cIyZXSupH9DHzC6SlAm8CHQGGgPvAYcB9YFGZjZLUjVgJtDbzBZJugfYbGb/lDQEqGVm/1NQjJ5gCrdpx27emL+GV2ZlM2dlUMzZpWVYzHlUQ6p5Madz5U4yJJguwJ1m1i18fwuAmf1fTJuJYZspktKAtUA9YEhs29h2uT7jVeBhM3tX0hLgVDNbI6kR8KGZHV5QjJ5gimb5xm8ZPzub8XOy+WrTTiqlpXBGZgP6dszg5Db1qJDqxZzOlQdRE0xaHGPIAFbGvF8FHJ9fGzPbJ2krUCdc/1mufTNid5TUHOgATA1XNTCzNeHyWoLbaD8h6WrgaoBmzZoV5XzKvRZ1D+F3ZxzG4NPbMOvroJjz9XmreX3eGmofUpGexzSid4cM2jf1Yk7nXHwTTNxIqgq8DAw2s225t5uZScrz0szMhgPDIbiCiWugZZQkjj20FsceWou/nJPJpC82MG52Ni9OX8nIKV/Rou4hYTFnYw6t48WczpVX8Uww2UDTmPdNwnV5tVkV3iKrQfCwP999JVUgSC4vmNkrMW3WSWoUc4tsfXGejMtbxbQUTs9swOmZDdi2ay9vz1/LuNnZ3P/fLxj63hd0bJZTzNmYWodUTHS4zrkSFM9nMGkED/lPI0gO04GLzWxhTJsbgKNjHvL3NbMLJbUF/sMPD/n/C7QBDgAjCR7mD871efcCm2Ie8tc2sz8VFKM/g4mf1d98x6tzgsE3v1i3g7QUcerh9enTIYPTjqzvxZzOlWIJf8gfBnEWcD9BN+WnzOzvku4CZpjZBEnpwHMEz1I2A/3MbFm4763AFcA+glthb0k6CZgMzCdINgB/NrM3JdUBxgDNgK8IuilvLig+TzDxZ2YsXrOdcbNX8eqc1azfvptqldI46+jgec3xLWqT4sWczpUqSZFgkp0nmJK1/4Ax5cugmPPtBWv4ds9+GtdI59z2GfTtmMFhDaolOkTnXASeYCLwBJM4O/fs491FQTHnpKUb2X/AyGxU/ftizgbV0xMdonMuH55gIvAEkxw27tjNa3NXM352NnNXbSVFcGLruvRun0G3oxpStVKp7OzoXJnlCSYCTzDJ58sNO3h1djbj5mSzcvN3pFdI4czMhvTpmMHJreuS5sWcziWcJ5gIPMEkLzNj5ldbGDc7m9fnrWHrd3upW7Ui5xzTmD4dMjimSQ0v5nQuQTzBROAJpnTYvW8/Hy7ZwPjZ2fx38Xr27D9Ay3qH0CecmbNp7SqJDtG5csUTTASeYEqfrd/t5a35a3hldjbTlge90DsdWove4cycNat4Madz8eYJJgJPMKXbqi07w2LObLLW76BiagrXndqKQb9q7QNvOhdHnmAi8ARTNpgZC1dv44nJy3h1zmraNa3J0Avb0bJe1USH5lyZFDXB+J95rtSTxFEZNXigXwceubgjKzZ+y9kPfszzn31Fef4DyrlE8wTjypSzj2nExMGn0Kl5LW4bv4CBI2ewYfvuRIflXLnkCcaVOQ1rpDPy8s7c2TOTT7I20u3+SbyzcG2iw3Ku3PEE48qklBRx2YkteP23J9GoRjpXPzeT/xk7jx279yU6NOfKDU8wrkxr06Aa464/ketPbcVLM1dy1gOTmfnVlkSH5Vy54AnGlXkV01L4U/cjGH1NFw6YccFjn3LfO0vYu/9A4Ts75w5agQlGUhNJf5D0qqTpkiZJGibpbEmenFypclzz2rx108n07diEh97Pou+wT8lavyPRYTlXZuWbJCQ9DTwF7AHuBvoD1wPvAd2BjyWdUhJBOldcqqVX4F8XtOOx33Rk1ZadnPPQZJ6dssK7MzsXB/kWWko6yswW5LujVBFoZmZZ8Qou3rzQsnxbv20Xfxw7j4++2MAvDqvHvecfQ32fh8a5QhVHoWUPSU3y22hmewpLLpK6S1oiKUvSkDy2V5I0Otw+VVLzmG23hOuXSOoWs/4pSeslLch1rHaSpkiaL+k1SdULis25+tXTeeby4/jfXm2ZunwT3e6fxNsL1iQ6LOfKjIISTGNgiqTJkq6XVK8oB5aUCjwC9AAygf6SMnM1GwhsMbPWwFCCW3GE7foBbQluxw0LjwfwTLgutyeBIWZ2NDAO+GNR4nXlkyQu6dKc1397Mk1rV+Ha52fxh5fmsn3X3kSH5lypl2+CMbPfAc2A24CjgXmS3pY0QFKUydM7A1lmtszM9gCjgF652vQCRobLY4HTFEzy0QsYZWa7zWw5kBUeDzObBGzO4/MOAyaFy+8C50WI0TkAWtevysvXdeW3v2rNK7NW0eOByd+P1uycOzgF9gSzwEdmdh3QhOAqYzCwLsKxM4CVMe9XhevybGNm+4CtQJ2I++a2kB8S2AVA07waSbpa0gxJMzZs2BDhNFx5USE1hZvPPJyXru1CisRFw6dw99ufs2efd2d27mBE6mos6WjgLoJbXruBW+IZ1EG6Arhe0kygGkHvt58ws+Fm1snMOtWrV6S7fq6cOPbQ2rx508lc1Kkpj374JX2GfcLSddsTHZZzpU5B3ZTbSPqLpIXAC8C3wJlmdoKZPRDh2Nn8+CqiSbguzzaS0oAawKaI+/6ImX1uZmea2bHAi8CXEWJ0Lk9VK6Xxz/OOYfglx7Jm6y7Oeehjnv5kOQcOeHdm56Iq6ArmbaAScJGZHWNm/zCzZUU49nSgjaQWYZfmfsCEXG0mAAPC5fOB9y3oNz0B6Bf2MmsBtAGmFfRhkuqH/6YQPDd6rAixOpenM9s2ZOLgUzixdV3++toiBjw9jbVbdyU6LOdKhYIe8rcys9vMbIGkQyWdDiCpcpSH/OEzlUHARGAxMMbMFkq6S9K5YbMRQB1JWcDvgSHhvguBMcAigkR3g5ntDz//RWAKcLikVZIGhsfqL+kL4HNgNfB00f5TOJe3etUqMWJAJ/7e5yhmrNhCt/sn8cY8787sXGEKndFS0lXA1UBtM2slqQ3wmJmdVhIBxpMXWrqiWrZhB78bM5e5K7+hT4cM/tqrLdXTKyQ6LOdKVHHOaHkDcCKwDcDMlgL1f154zpVOLetVZey1XRh8ehsmzF1Nj/snM3XZpkSH5VxSipJgdod1LMD3D+P9SacrtyqkpjD49MMYe20XKqSKfk98xv+9uZjd+/YnOjTnkkqUBPORpD8DlSWdAbwEvBbfsJxLfh2a1eLNm06mf+dmPD5pGb0f+ZQla707s3M5oiSYIcAGYD5wDfAmQS8t58q9KhXT+EefoxkxoBMbtu+i58Mf8+TkZd6d2TkiPOQvy/whvytOG3fsZsjL83lv8TpObF2Hf13QjkY1Kic6LOeK3c9+yB+OSNxT0k+6yEhqGXY3vuLnBupcWVG3aiWeuPRY/tn3aGZ//Q3dhk5iwtzViQ7LuYQp6BbZVcDJwOfhbJZvSnpf0nLgcWCmmT1VIlE6V0pIol/nZrx548m0ql+VG1+czU2jZrN1p4/O7MqfSLfIwnlaGgHfAV+Y2c74hlUy/BaZi6d9+w/w6Idfcv9/l1K/WiXuu6AdXVvXTXRYzv1sxVkHg5mtMLMpZjanrCQX5+ItLTWF357Whleu60rlCqlc/ORU/vb6Inbt9e7MrnyIlGCccwevXdOavHHjyVxywqE8+fFyej/yCYvXbEt0WM7FnScY50pA5Yqp/G/vo3j68uPY9O0eej38CcMnfendmV2ZVmiCCXuSeSJyrhj88vD6TBx8Cr88oh7/ePNzLn7yM7K/+S7RYTkXF1ESx0XAUkn3SDoi3gE5V9bVPqQij/3mWO49/xjmr9pK9/snMX52NuW5Js2VTYUmGDP7DdCBYAKvZyRNCacdLnTIfudc3iRxQaemvHXTKRzeoBqDR89h0Iuz+WZnnhOxOlcqRe1Ftg0YC4wi6K7cB5gl6bdxjM25Mq9ZnSqMvqYLf+x2OBMXrKX7/ZP5eOnGRIflXLGI8gzmXEnjgA+BCkBnM+sBtANujm94zpV9qSnihl+2ZvwNJ1I1PY3fjJjKX19b6N2ZXakX5QrmPGComR1tZvea2XqAsB5mYMG7OueiOiqjBq//9iQu69qcpz9ZQc+HPmbh6q2JDsu5gxblGcwAM5uUz7b/FrSvpO6SlkjKkjQkj+2VJI0Ot08NRwzI2XZLuH6JpG4x65+StF7SglzHai/pM0lzJM2Q1Lmwc3Mu2aRXSOXOc9sy8orObP1uL70f+YRHP/yS/d6d2ZVCUW6RbZe0LddrpaRxkloWsF8q8AjQA8gE+kvKzNVsILDFzFoDQ4G7w30zgX5AW6A7MCw8HsAz4brc7gH+ambtgdvD986VSr84rB4TB5/CGZkNuPvtz+k//DNWbvZBNFzpEuUW2f3AH4EMoAnwB+A/BA/8CxrssjOQZWbLwhkxRwG9crXpBYwMl8cCp0lSuH6Ume02s+VAVng8wqupzXl8ngHVw+UagA9j60q1WodU5JGLO3LfBe1YtGYbPR6YzNiZq7w7sys1oiSYc83scTPbbmbbzGw40M3MRgO1CtgvA1gZ835VuC7PNma2D9gK1Im4b26DgXslrQT+BdxSSHvnkp4kzju2CW/ddDKZjarzh5fmcv0Ls9jyrXdndskvSoLZKelCSSnh60JgV7gtmf6Uug74nZk1BX4HjMirUVjDM0PSjA0bNpRogM4drKa1q/Di1ScwpMcRvLd4Hd3un8RHX/jPr0tuURLMr4FLgPXAunD5N5IqA4MK2C8baBrzvkm4Ls82ktIIbm1tirhvbgOAV8LllwhvqeVmZsPNrJOZdapXr14hh3QueaSmiGt/0YrxN5xIzSoVGPDUNO54dQHf7fHuzC45FZhgwgfrvcysp5nVNbN64XKWmX1nZh8XsPt0oI2kFpIqEjy0n5CrzQSCxABwPvC+BTeYJwD9wl5mLYA2wLRCzmU18Itw+VfA0kLaO1cqtW1cgwmDTuKKE1swcspXnPPQZOav8u7MLvkUmGDMbD/Q/2AOHD5TGQRMBBYDY8xsYTjV8rlhsxFAHUlZwO+BIeG+C4ExwCLgbeCGMBYkvQhMAQ6XtEpSTi3OVcB9kuYC/wCuPpi4nSsN0iukcnvPTJ4feDzf7t5Pn2Gf8PD7S707s0sqhc5oKWkoQQX/aODbnPVmNiu+ocWfz2jpyoJvdu7htvELeH3eGo49tBZDL2xPszpVEh2WK8OizmgZJcF8kMdqM7NfHWxwycITjCtLXp2TzW3jF3DggHFHz7Zc0KkJQa9/54pX1ASTVlgDM/tl8YTknIunXu0z6NS8NjePmcOfXp7He4vX8X99j6ZO1UqJDs2VU1Eq+RtIGiHprfB9ZsxzD+dcEsmoWZn/XHkCt551JB8u2UC3+yfzwefrEx2WK6eidFN+huBBfePw/RcERY3OuSSUkiKuOqUlrw46kbpVK3L5M9O5ddx8du7Zl+jQXDkTJcHUNbMxwAH4vneYd7x3Lskd2ag64284katObsF/pn3NOQ9+zNyV3yQ6LFeOREkw30qqQ1i1L+kEgiFdnHNJLr1CKreenckLVx7Prr376fvopzzw3lL27T+Q6NBcORAlwfyeoPCxlaRPgGcBn8nSuVKka6u6vDX4FHoe04ih733BBY9PYcXGbwvf0bmfodBuyvD9MC6HAwKWmNneeAdWErybsiuPXpu7mlvHzWffAeMv52TS77im3p3ZFUnUbspRrmAgGNerHdCRYF6XS39OcM65xOnZrjETf3cKHZrV5JZX5nPVszPZuGN3osNyZVCUbsrPEQx/fxJwXPgqNHM555JXoxqVee6K4/nLOZlMWrqBbkMn8d6idYkOy5UxhRZaEiSTTPNZjpwrU1JSxMCTWnBS67oMHj2HK5+dQf/Ozbjt7CM5pFKUXw3OFSzKLbIFQMN4B+KcS4zDG1Zj/A1dufYXrRg1/WvOfnAys77ekuiwXBkQqQ4GWCRpoqQJOa94B+acKzmV0lIZ0uMIRl11Anv3Gxc8NoV/v/sFe707s/sZolwH3xnvIJxzyeH4lnV4a/DJ3DlhIQ/+dykfLVnP0Iva07Je1USH5kqhQq9gzOwjYAVQIVyeDpT6ofqdc3mrnl6Bf1/YnmG/7shXm3dy9oMf8/xnX+GPYV1RRelFdhUwFng8XJUBjI9nUM65xDvr6EZMHHwKnZrX4rbxCxg4cgbrt+9KdFiuFInyDOYG4ERgG4CZLQXqxzMo51xyaFA9nZGXd+bOnpl8krWR7vdPZuLCtYkOy5USURLMbjPbk/MmrOr3a2XnyomUFHHZiS1448aTaFQjnWuem8mfxs5lx24fndkVLEqC+UjSn4HKks4AXgJei3JwSd0lLZGUJWlIHtsrSRodbp8qqXnMtlvC9UskdYtZ/5Sk9ZIW5DrWaElzwtcKSXOixOici6Z1/WqMu/5Erj+1FWNnruKsByYz86vNiQ7LJbEoCWYIsAGYD1wDvAncVthOklKBR4AeQCbBEDOZuZoNBLaYWWtgKHB3uG8m0A9oC3QHhoXHg2B+mu65P8/MLjKz9mbWHngZeCXCuTnniqBiWgp/6n4Eo6/pghF0Z/7XxCXendnlKUovsgNm9oSZXWBm54fLUW6RdQayzGxZeIttFNArV5tewMhweSxwmoJR93oBo8xst5ktB7LC42Fmk4B8/2wK978QeDFCjM65g3Bc89q8eeNZZ213AAAbbElEQVTJnNexCQ9/kEXfYZ+StX5HosNySSbqYJcHIwNYGfN+VbguzzbhRGZbgToR983PycC6sDPCT0i6WtIMSTM2bNgQ8ZDOudyqpVfg3gva8dhvOrJqy07OeWgyz05Z4d2Z3ffimWASpT8FXL2Y2XAz62RmnerVq1eCYTlXNnU/KujOfELLOtz+6kIGPD2dddu8O7OLb4LJBprGvG8SrsuzTdg7rQawKeK+PxEeoy8w+qCjds4VWf3q6Tx92XH8b6+2TFu+iW73T+Kt+WsSHZZLsCiFlodJekLSO5Lez3lFOPZ0oI2kFpIqEjy0zz2G2QRgQLh8PvB++HxnAtAv7GXWAmgDTIvwmacDn5vZqghtnXPFSBKXdGnOGzeeTLPaVbjuhVncPGYu23eVifkJ3UGIMhbZS8BjwBPA/qgHNrN9kgYBE4FU4CkzWyjpLmCGmU0ARgDPScoieHDfL9x3oaQxwCJgH3CDme0HkPQicCpQV9Iq4A4zGxF+bD/84b5zCdWqXlVevq4rD/13KQ9/kMWsr7cw8vLONKtTJdGhuRJW6JTJ4dSYx5ZQPCXKp0x2Lr6mLd/M1c/NIC1FPH1ZZ45uUiPRIbliUJxTJr8m6XpJjSTVznkVQ4zOuTKuc4vavHxdV9IrpHLR8Cl8uGR9okNyJShKghkA/BH4FJgZvvzPfudcJK3qVeWV67vSou4hDBw5g5dmrCx8J1cmFPoMxsxalEQgzrmyq361dEZf04Xrnp/JH8fOY+3WXQz6VWuCumhXVkXpRVZB0o2SxoavQZIqlERwzrmyo2qlNEYMOI6+HTO4790vuHX8Avb5EDNlWpReZI8CFYBh4ftLwnVXxiso51zZVDEthfsuaEejGuk88sGXrN+2m4f6d6ByxdTCd3alTpQEc5yZtYt5/76kufEKyDlXtknij92OoGH1dG6fsJCLn/yMEQOOo/YhFRMdmitmUR7y75fUKueNpJYUoR7GOefyckmX5jz662NZtHob5z36KV9v2pnokFwxi5Jg/gh8IOlDSR8B7wM3xzcs51x50P2ohvznquPZsnMPfR/9hPmrtiY6JFeMogzX/1+CoVpuBH4LHG5mH8Q7MOdc+XDsobUZe21XKqV5rUxZk2+CkfSr8N++wNlA6/B1drjOOeeKRev6VRl3fVea1zmEK71Wpswo6CH/Lwhuh/XMY5vhM0Y654pR/erpjL7mBK5/YRZ/HDuPddt2ccMvvVamNMs3wZjZHeHiXeGskt8LRzh2zrliVS29AiMGHMeQl+fxr3e+YM3WXdzV6yhSUzzJlEZRHvK/nMe6scUdiHPOQVgrc2E7rj+1FS9M/Zprn5/Jd3u842pplO8VjKQjgLZAjVzPXKoD6fEOzDlXfkniT92PoGGNdO7wWplSq6ArmMOBc4CaBM9hcl4dgaviH5pzrry7NKZW5vxHP2XlZq+VKU2izAfTxcymlFA8Jcrng3GudJixYjMDR86gQmoKz1x+HEdl+LwyiVSc88HMlnSDpGGSnsp5FUOMzjkXSafmtXn5ui5USkvhosen8NEXGxIdkosgSoJ5DmgIdAM+ApoA26McXFJ3SUskZUkaksf2SpJGh9unSmoes+2WcP0SSd1i1j8lab2kBXkc77eSPpe0UNI9UWJ0zpUOretXY9z1XTm0ziEMfGY6Y2euSnRIrhBREkxrM/sL8K2ZjSQoujy+sJ0kpQKPAD2ATKC/pMxczQYCW8ysNTAUuDvcNxPoR9DJoDswLDwewDPhutyf90ugF9DOzNoC/4pwbs65UiSnVuaElnX4w0tzeeSDLAq7ze8SJ0qC2Rv++42ko4AaQP0I+3UGssxsmZntAUYRJIBYvYCR4fJY4DQFVVW9gFFmtjuswckKj4eZTQI25/F51wH/NLPdYTsfb8K5MqhaegWeuuw4+nTI4N6JS/jLqwvYf8CTTDKKkmCGS6oF/AWYACwivNIoRAYQO97DqnBdnm3MbB+wFagTcd/cDgNODm+1fSTpuLwaSbpa0gxJMzZs8Pu4zpVGFdNS+PeF7bju1FY8/9nXXPf8THbt9VqZZBMlwTxtZlvM7CMza2lm9c3s8bhHVnRpQG3gBIIRoMcojzEmzGy4mXUys0716tUr6Ridc8VEEv/T/Qj+em5b3l28jouf+Iwt3+5JdFguRpQEs1zScEmn5fULuwDZQNOY903CdXm2kZRGcPttU8R9c1sFvGKBacABoG4R4nXOlUIDujbn0V93ZMHqbZz3mNfKJJMoCeYI4D3gBmCFpIclnRRhv+lAG0ktJFUkeGg/IVebCcCAcPl84H0LnthNAPqFvcxaEEwXMK2QzxsP/BJA0mFARWBjhDidc6Vc96Ma8cKVx7Npxx76DPuUBdk+r0wyiDIfzE4zG2NmfYH2BEPFfBRhv33AIGAisBgYY2YLJd0l6dyw2QigjqQs4PfAkHDfhcAYguc9bwM3mNl+AEkvAlOAwyWtkjQwPNZTQMuw+/IoYIB59xLnyo3jctXKTPJamYQrtJIfQNIvgIsIugfPAEabWV6DYJYqXsnvXNmzbtsuLnt6OkvXbefu847hvGObJDqkMqfYKvklrQAGA5OBo83swrKQXJxzZVOD6umMueYEjm9Zm5u9ViahCppwLMcxZrYt7pE451wxqZZegacv68yfxs7l3olLWLP1O/56rs8rU9IKGq7/T2Z2D/B3ST9J/2Z2Y1wjc865nyGolWlPwxqVeeyjL1m/bTcP9u9AeoXUwnd2xaKgK5jF4b/+kMI5VyqlpIghPY6gYfVK/PX1Rfz6yak8eWknavm8MiWioCmTXwsXd5rZS7HbJF0Q16icc64YXXZiCxpUT+em0XM477FPGXl5Z5rWrpLosMq8KHUwt0Rc55xzSavH0T/UyvR91GtlSkK+CUZSD0kPARmSHox5PQPsK7EInXOumOTUylRMDWplJi/1Wpl4KugKZjXB85ddwMyY1wSCuWGcc67UaV2/Gq9c35Wmtatw+dPTeWWWzysTLwU9g5kLzJU0jmAumJxK+lSgUgnF55xzxa5B9XTGXNuF656fye/HzGXN1l1cf2orijbcoitMlGcw7wCVY95XJhibzDnnSq3qYa1Mr/aNuXfiEm5/daHPK1PMohRappvZjpw3ZrZDkne/cM6VehXTUhh6YXsaVk/n8UnLWLdtl9fKFKMoVzDfSuqY80bSscB38QvJOedKTkqKuOWsI7mjZybvLl7Hr5+c6vPKFJMoCWYw8JKkyZI+BkYTjJLsnHNlxuUntmDYxR2Zn73V55UpJlGG659OMCfMdcC1wJFmNjPegTnnXEnrcXQjnh94PBu376bvo5+ycLXXyvwcUa5gAA4HMoGOQH9Jl8YvJOecS5zOLWrz8nVdqZAiLnr8M6+V+RmiDNd/B/BQ+PolcA9wboE7OedcKdamQTVeuf5EmtSq7LUyP0OUK5jzgdOAtWZ2OdAOqBHXqJxzLsEa1ghqZTq3qM3vx8xl2Ic+r0xRRUkw35nZAWCfpOrAeqBplINL6i5piaQsSUPy2F5J0uhw+1RJzWO23RKuXyKpW8z6pyStD6dGjj3WnZKyJc0JX2dFidE55/JTPb0Cz1zemXPbNeaet5dwxwSvlSmKKHUwMyTVBJ4gGCpmBzClsJ3Civ9HgDOAVcB0SRPMbFFMs4HAFjNrLakfcDdwkaRMoB/QFmgMvCfpsHA0gWeAh4Fn8/jYoWb2rwjn5JxzkVRMS+H+i9rTqMYPtTIP9PNamSii9CK73sy+MbPHCJLFgPBWWWE6A1lmtszM9gCjgF652vQCRobLY4HTFIzV0AsYZWa7zWw5kBUeDzObBGyO8PnOOVcscmplbj8nk3cWreM3T07lm51eK1OYqL3IADCzFWY2L2LzDGBlzPtV4bo825jZPmArUCfivnkZJGleeButVl4NJF0taYakGRs2eO8Q51x0V5zUgkcu7si87K2c9+inrNritTIFKVKCSXKPAq2A9sAa4L68GpnZcDPrZGad6tWrV5LxOefKgLOObsRzV3Rmw/bd9B3mtTIFiWeCyebHnQGahOvybCMpjaB32qaI+/6Ima0zs/1hh4QnCG+pOedccTu+ZR3GXteVtLBW5uOlGxMdUlKKZ4KZDrSR1EJSRYKH9hNytZkADAiXzwfet6Af4ASgX9jLrAXQBphW0IdJahTztg+wIL+2zjn3cx0WUytz2dPTGDfba2Vyi1uCCZ+pDAImAouBMWa2UNJdknIKNUcAdSRlAb8HhoT7LgTGAIuAt4EbYuajeZGgF9vhklZJGhge6x5J8yXNIygI/V28zs055+CHWpnjmtfmd6Pn8uiHX3qtTAyV5/8YnTp1shkzZiQ6DOdcKbd7337+8NI8Xpu7mgFdDuX2nm1JTSm7k5dJmmlmnQprF6UOxjnnXAEqpaXyQFgrM3zSMtZt2839/dqX+1qZstSLzDnnEiYlRfw5rJWZuGgtl4zwWhlPMM45V4yuOKkFD/fvyNyVWzn/sSnlulbGE4xzzhWzs49pxLMDO7N+2y76DvuURau3JTqkhPAE45xzcXBCWCuTmiIufHwKn2SVv1oZTzDOORcnQa1M1+9rZcbPLrBevMzxBOOcc3HUqEZlxlzbhWMPrcXg0XN47KPyUyvjCcY55+KsenoFRl7RmXOOacQ/3/qcv762qFzMK+N1MM45VwIqpaXyYL8ONKqRzhOTl7N2664yXyvjVzDOOVdCUlLErWdn8pdyUivjCcY550rYwJNa8FD/Dt/XymR/812iQ4oLTzDOOZcA5xzTmGcHdmbdtl30HfZJmayV8QTjnHMJckLLOoy9tispKpu1Mp5gnHMugQ5vGNTKZNQMamVenVN2amU8wTjnXILl1Mp0bFaLm0bN4fEyUivjCcY555JAjcoVeHZgZ84+phH/V0ZqZbwOxjnnkkSltFQe6teBRtXTefLj5azbtouhF5XeWpm4XsFI6i5piaQsSUPy2F5J0uhw+1RJzWO23RKuXyKpW8z6pyStl7Qgn8+8WZJJqhuPc3LOuXhKSRG3nZPJbWcfyVsL1nLpiGls3bk30WEdlLglGEmpwCNADyAT6C8pM1ezgcAWM2sNDAXuDvfNBPoBbYHuwLDweADPhOvy+symwJnA18V6Ms45V8KuPLklD/XvwJyV33D+Y5+WylqZeF7BdAayzGyZme0BRgG9crXpBYwMl8cCp0lSuH6Ume02s+VAVng8zGwSsDmfzxwK/Ako3TcunXMO6NmuMSOv6MzasFZm8ZrSVSsTzwSTAayMeb8qXJdnGzPbB2wF6kTc90ck9QKyzWxuIe2uljRD0owNGzZEOQ/nnEuYLq2CWhkhLnxsCp+WolqZMtGLTFIV4M/A7YW1NbPhZtbJzDrVq1cv/sE559zPlFMr06hmOgNKUa1MPBNMNtA05n2TcF2ebSSlATWATRH3jdUKaAHMlbQibD9LUsOfEb9zziWNxjUr89K1Xb+vlRk+KflrZeKZYKYDbSS1kFSR4KH9hFxtJgADwuXzgfct+C82AegX9jJrAbQBpuX3QWY238zqm1lzM2tOcEuto5mtLd5Tcs65xImtlfnHm8lfKxO3Ohgz2ydpEDARSAWeMrOFku4CZpjZBGAE8JykLIIH9/3CfRdKGgMsAvYBN5jZfgBJLwKnAnUlrQLuMLMR8ToP55xLJjm1Mg2rpzPi4+Ws376Lf1+YnLUySvZLrHjq1KmTzZgxI9FhOOfcQXly8jL+9sZiOjevzROXdqJGlQol8rmSZppZp8LalYmH/M45Vx5deXJLHoyplVmdZLUynmCcc64UOzemVqZPktXKeIJxzrlSrkurOrx0bZcfamW+TI5aGU8wzjlXBhzRsPoPtTJPTWPC3NWJDskTjHPOlRWNa1bmpWu60qFZLW58cTZPTFqW0FoZTzDOOVeG1KhSgWev6MzZRzfi728u5q7XF3EgQbUyPh+Mc86VMekVUnmofwcaVE/nqU+Ws37bbu67sF2J18p4gnHOuTIoJUXc3jOTxjXT+dsbi9mwYzdPXFJytTLgt8icc65My6mVmf31Fi54vGRrZTzBOOdcGZdTK7Pmm130HfYpn68tmVoZTzDOOVcOdG1Vl5eu64JhXPDoFGasyG/exuLjCcY558qJoFbmRNo3q0nDGulx/zx/yO+cc+VIRs3KPDfw+BL5LL+Ccc45FxeeYJxzzsWFJxjnnHNx4QnGOedcXHiCcc45FxdxTTCSuktaIilL0pA8tleSNDrcPlVS85htt4Trl0jqFrP+KUnrJS3Idaz/lTRP0hxJ70hqHM9zc845V7C4JRhJqcAjQA8gE+gvKTNXs4HAFjNrDQwF7g73zQT6AW2B7sCw8HgAz4TrcrvXzI4xs/bA68DtxXtGzjnniiKeVzCdgSwzW2Zme4BRQK9cbXoBI8PlscBpkhSuH2Vmu81sOZAVHg8zmwT8pATVzGLHPjgESNwkCM455+JaaJkBrIx5vwrIXd3zfRsz2ydpK1AnXP9Zrn0zCvtASX8HLgW2Ar/Mp83VwNXh2x2SlhR6JolXF0iOOVB/Pj+X5FRWzqWsnAck97kcGqVRmarkN7NbgVsl3QIMAu7Io81wYHhJx/ZzSJphZp0SHUdx8HNJTmXlXMrKeUDZOJd43iLLBprGvG8SrsuzjaQ0oAawKeK+BXkBOK+I8TrnnCtG8Uww04E2klpIqkjw0H5CrjYTgAHh8vnA+xZMID0B6Bf2MmsBtAGmFfRhktrEvO0FfF4M5+Ccc+4gxe0WWfhMZRAwEUgFnjKzhZLuAmaY2QRgBPCcpCyCB/f9wn0XShoDLAL2ATeY2X4ASS8CpwJ1Ja0C7jCzEcA/JR0OHAC+Aq6N17klQKm6pVcIP5fkVFbOpaycB5SBc1FwweCcc84VL6/kd845FxeeYJxzzsWFJ5gkIekmSQskLZQ0OFxXW9K7kpaG/9YK10vSg+FQOvMkdUxs9D/I5zzulJQdDuMzR9JZMe3zHBIoEfIahuhgvgNJA8L2SyUNyOuzkuxcTpW0Neb7uT1mnwKHe0rguVwQ/owdkNQpV/v8hplK6LkU5TwkNZf0Xcx38ljMtmMlzQ/P48GwOD05mZm/EvwCjgIWAFUIOl68B7QG7gGGhG2GAHeHy2cBbwECTgCmJvocCjmPO4E/5NE+E5gLVAJaAF8CqQmM/xSgI7AgZl2RvgOgNrAs/LdWuFwryc/lVOD1PI6RGn4nLYGK4XeVmSTnciRwOPAh0Kmwn6lkOJcinkfz2Ha5jjMt/JlT+DPYo6S/k6gvv4JJDkcS/ILaaWb7gI+Avvx4KJ2RQO9wuRfwrAU+A2pKalTSQechv/PIT75DAiWC5T0MUVG/g27Au2a22cy2AO+S99h5cVXEc8lPlOGe4i6vczGzxWaW1ygc+f1MJfxcingeeQp/xqqb2WcWZJtnKfx7TBhPMMlhAXCypDqSqhD8ddwUaGBma8I2a4EG4XJew/AUOpROCcjvPAAGhbeSnsq5NUPynkeson4HyXxO+Z0LQBdJcyW9JaltuC6ZzyU/pfF7yU8LSbMlfSTp5HBdBkHsOZL6PDzBJAEzW0wwkvQ7wNvAHGB/rjZGkg/gWcB5PAq0AtoDa4D7EhXjz1EavoOocp3LLOBQM2sHPASMT1hgLscaoJmZdQB+D/xHUvUEx1RknmCShJmNMLNjzewUYAvwBbAu59ZX+O/6sPnPHUonbvI6DzNbZ2b7zewA8AQ/3AZL2vOIUdTvIJnPKc9zMbNtZrYjXH4TqCCpLsl9Lvkpjd/LT4S3+DaFyzMJnh8dRhBzk5imSX0enmCShKT64b/NCJ5b/IcfD6UzAHg1XJ4AXBr2ZDoB2Bpz6yOh8jqPXM+H+hDcSoODGBIoAYr6HUwEzpRUK7wVeGa4LhnkeS6SGub0RJLUmeD3wiaiDfeUbPL7mSpV5yKpnsI5sCS1JDiPZeHP2DZJJ4Tf2aX88DOZfBLdy8BfwQuYTDA0zlzgtHBdHeC/wFKCHlm1w/UimMztS2A+Mb1PEv3K5zyeC+OcR/A/daOY9reG57GEBPeGAV4kuDWxl+De9sCD+Q6AKwgeLmcBl5eCcxkELAy/s8+ArjHHOYvgavpL4NYkOpc+4fJuYB0wsbCfqUSfS1HOg2Cw3oUEt5lnAT1jjtOJ4I+0L4GHCUdkScaXDxXjnHMuLvwWmXPOubjwBOOccy4uPME455yLC08wzjnn4sITjHPOubjwBOPKLUlTw5Fqv5a0IWbk2uaJji0/kvpKOiKOx0+R9IGkqpLSJH0Ts61nOBpxU0mDJV0arzhc2RC3KZOdS3ZmdjyApMsI6lgGJTaigKQ0CwYLzUtfgmnBPy+m4+XWk2BK8x2Svv/9EA57PxQ4w8xWSnqSYDDTZ6PG4cofv4JxLg+SekiaImmWpNGSDgnXr5L0j3BgyOmSOkp6R9KXkq4K25weXgW8Ff7F/0hMpXxBx/2npNlAH0nXhsefK+klSZXDAQ/PAobmXGlJ+lhS+/AYDSVlhctXShov6QPCkQQkDZE0LRx09Pbc5xz6NbkqwyX9EhhGULS4HMCCoWVWK4nmInLJxxOMc7mEw90MIRiJoCPBCAQ3xTRZbsHAkJ8BIwiqsbsC/xvT5njgOoL5SY4EekU47noz62BmLwEvmdlx4ed8CVxmZpOBN4HfmVl7M1tRyKl0APqa2WkKJnlrFsbVHugqqWse+3QlqBzPUQV4GehlZktztZ0BnIxz+fBbZM79VFeCxPBpeOFREfg4ZnvOGFbzgTQz+xb4VsGshFXDbZ/lJABJo4CTwvUFHXd0zPIxku4CagLVgNcP4jzesWBOGgjGROsBzA7fVyUYPPHTXPtUN7OdMe93AVOBy4Gbc7VdTzAxlnN58gTj3E8JeNvMLsln++7w3wMxyznvc/6fyj0Gk0U47rcxy88S3JJaIOlKghkM87KPH+5EpBdwPAF/M7MR+Rwnx4E83p8PvC/pT2Z2T8y2dOC7Qo7nyjG/RebcT30K/CIcxRZJh0hqU8RjnCCpWTgi7oUEVypFOe4hwFpJFYCLY9ZvJ7iiybECODZcPr+AeCYCA2Oe+TRRMCR/blm5e9GFV2hnA5dLGhCz6TB+GBnbuZ/wBONcLma2jmCk29GS5hIkhsOKeJhpwGMEI0svASYU8bi3Ewwx/0l4jBwvAn+O6U59L3CTpFlArdwHiTmnN4GxwGeS5gNjCG6T5fYGcGoe+28kmPr5r5LODld3IRid2bk8+WjKzhUzSacDg8wsaedKz4+kJsCTZta9kHbHAdeb2eUlE5krjfwKxjn3PTNbBTwT01khP7WBO0ogJFeK+RWMc865uPArGOecc3HhCcY551xceIJxzjkXF55gnHPOxYUnGOecc3Hx/+w4aNwI0MDEAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "per_ea = calculate_activation_energies( per_resistivity_ratio_list, temp ) \n", "par_ea = calculate_activation_energies( par_resistivity_ratio_list, temp ) \n", "\n", "plt.plot( temp[1:-1], per_ea[1:-1] )\n", "plt.title( 'Perpendicular activation energy')\n", "plt.xlabel( ' Temperature (K)' )\n", "plt.ylabel( ' activation energy (eV)' )\n", "plt.show()\n", "\n", "plt.plot( temp[1:-1], par_ea[1:-1] )\n", "plt.title( 'Parallel activation energy')\n", "plt.xlabel( ' Temperature (K)' )\n", "plt.ylabel( ' activation energy (eV)' )\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.1" } }, "nbformat": 4, "nbformat_minor": 2 }