{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"1D plotting utilities\n",
"===============\n",
"\n",
"**Suhas Somnath**\n",
"\n",
"8/12/2017\n",
"\n",
"### Table of contents:\n",
"* [use_scientific_ticks](#use_scientific_ticks)\n",
"* [plot_curves](#plot_curves)\n",
"* [plot_line_family](#plot_line_family)\n",
"* [plot_complex_spectra](#plot_complex_spectra)\n",
"* [rainbow_plot](#rainbow_plot)\n",
"* [plot_scree](#plot_scree)\n",
"\n",
"#### Import necessary packages:"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from __future__ import division, print_function, absolute_import, unicode_literals\n",
"import numpy as np\n",
"from warnings import warn\n",
"import matplotlib.pyplot as plt\n",
"import subprocess\n",
"import sys\n",
"\n",
"\n",
"def install(package):\n",
" subprocess.call([sys.executable, \"-m\", \"pip\", \"install\", package])\n",
"# Package for downloading online files:\n",
"try:\n",
" import sidpy\n",
"except ImportError:\n",
" warn('sidpy not found. Will install with pip.')\n",
" import pip\n",
" install('sidpy')\n",
" import sidpy"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"use_scientific_ticks() \n",
"------------------------\n",
"Often scientific plots look ugly because of the way tick marks are formatted for small or large values.\n",
"use_scientific_ticks() is a handy function that can convert tick marks on 1D plots to scientific form"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAEYCAYAAABRMYxdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd4XOWZ+P3vrWpLLrKqu+UmdzBYNqa5ggubYJJ1EpNGEjYsm7Dv7pLNhpQfISTZDdlClg0pJCQhhUBCKA6huRvHuMjG3ZYsF9ylUbFsy+q63z/mjBnEjOqMzpT7c11zMXPmPGV0cR7f52lHVBVjjDHGmFiS4HYFjDHGGGNCzQIcY4wxxsQcC3CMMcYYE3MswDHGGGNMzLEAxxhjjDExxwIcY4wxxsQcC3BMzBCRfxCRMhG5JCJZPczruIjcEqq6GWPcJyL7RWRemPJ+VUTu8vv8HRGpEJFzIjLSaZcSw1G2CUxsHxwTKUTkOJAHNAMtwAHg18ATqtraQdpk4AIwW1V3h6guf6eqq0XkIWCcqn6yp/kaY0JDRG4Cvg9MwdteHAT+WVW390LZD9FOmyAiI4FiYJSqloe7PiYw68ExkeaDqtofGAV8D/gK8GQn0uUBfYD9YaybMSYCiMgA4GXg/4BMYBjwLaDBzXr5GQlUWnDjLgtwTERS1RpVXQl8DLhLRKaKSKqI/JeInHCGon4iIn1FpADv3RLAeRFZCyAi/ysiJ0XkgojsEJGbffmLyK9E5Dt+n+eJyKm29RCRJcDXgI85Xcw97h0yxvRYAYCq/l5VW1S1TlXfUNU9ACLyeRE5KCIXReSAiFzrHL8y9CwiCSLygIgcEZFKEfmDiGQ63+WLiIrIXU57UyEiX3e+C9gmiMh6Efk7J/9VwFDn+1/55ZfknJspIr8UkTMiUi0iL/buny8+WIBjIpqqbgNOATfj7dEpAKYD4/DetT2oqiV4u6kBMlR1gfN+u3NuJvA08EcR6dPF8l8D/h14VlX7qerVPfxJxpieKwFaROQpEVkqIoN8X4jIR4CHgE8DA4DbgcoAefwjcAcwFxgKVAOPtznnJmACsBB4UEQmddQmqOpqYClwxvn+MwHK/g2QhrfdygUe7cJvN51kAY6JBmfwBin3AP+iqlWqehFvI7MiWCJV/a2qVqpqs6r+N5CKt7EyxkQxVb2AN/hQ4GeAR0RWikge8HfA91V1u3qVquo7AbK5F/i6qp5S1Qa8QdFyXy+L41tO79BuYDfQ4xscERmCNwC6V1WrVbVJVTf0NF/zfkkdn2KM64bh/X81DdghIr7jAgRdlSAi/wrcjffuTPHezWWHtabGmF6hqgeBzwCIyETgt8APgBHAkU5kMQp4QUT8FzC04J3P53PO7/1loF8PquwzAqhS1eoQ5GXaYT04JqKJyEy8Ac6LQB0wRVUznNdAVQ3Y4Djzbf4N+CgwSFUzgBq8QRFALd6AyWdwO9WwpYbGRDBVPQT8CpgKnATGdiLZSWCpX3uSoap9VPV0Z4rsfm05CWSKSEYP8jCdYAGOiUgiMkBEPgA8A/zW6SL+GfCoiOQ65wwTkcVBsuiPd7m5B0gSkQfx9uD47AJucyb7DQb+uZ3qlAH5ImLXizERQEQmisiXRGS483kEcCewBfg58K8iMkO8xonIqADZ/AT4ru87EckRkWWdrEK32wRVPQu8CvxIRAaJSLKIzOlqPqZj1mCbSPNnEbmI9y7n68D/AJ91vvsKUApsEZELwGqCz6l5HXgN72TEd4B6J0+f3+AdUz8OvAE8206d/uj8t1JEdnbx9xhjQu8icB2wVURq8QY2+4Avqeofge/iXVhwEW/vb2aAPP4XWAm84bQ5W5w8O6OnbcKngCbgEFBO+zdYpptsoz9jjDHGxBzrwTHGGGNMzLEAxxhjjDExxwIcY4wxxsQcC3CMMcYYE3Nso78Qyc7O1vz8fLerYUzE27FjR4Wq5rhdj0hjbYgxndPZNsQCnBDJz8+nqKjI7WoYE/FEJNC2+XHP2hBjOqezbYgNURljjDEm5liAY4wxxpiYYwGOMcYYY2KOBTjGGGOMiTkW4BhjjDEm5liAY4wxxpiYE5IAR0SWiEixiJSKyAMBvk8VkWed77eKSL7fd191jheLyOKO8hSR0U4epU6eKe2VISJZIrJORC6JyA/b1GuGiOx10jwmIuIczxSRVSJy2PnvoFD8nYwxoScivxCRchHZF+R7ca7vUhHZIyLX+n13l3OdHxaRu7pRdsD0IpIiIk+ISImIHBKRv+3erzPGdFePAxwRSQQeB5YCk4E7RWRym9PuBqpVdRzwKPCIk3YysAKYAiwBfiQiiR3k+QjwqJNXtZN30DKAeuD/Af8aoPo/Bj4PjHdeS5zjDwBrVHU8sMb5bIyJTL/i3Ws3kKW8e43fg/e6R0QygW8C1wGzgG925Wamg/RfB8pVtQBvG7ahC7/HGBMCodjobxZQqqpHAUTkGWAZcMDvnGXAQ87754AfOr0ly4BnVLUBOCYipU5+BMpTRA4CC4CPO+c85eT742BlqGotsElExvlXWkSGAANUdYvz+dfAHcCrTl7z/MpYD3yl63+a6HXw7AXWF3vYd6aGC3VNpKUkMjq7H3MLcpiZP4ikRBvdNJFBVTf69woHsAz4taoqsEVEMpzrfx6wSlWrAERkFd5A6fcisgj4FpAKHAE+q6qX2uS7OFh64HPARKd+rUBFCH5qVDlcdpG1h8rZe7qGmrom+iYnMjo7nTkFOcwanUmytSEmzEIR4AwDTvp9PoX3jibgOaraLCI1QJZzfEubtMOc94HyzALOq2pzgPODlRGsYRnmpA9Udp6qnnXenwPyAmUgIvfgvSNk5MiRQYqJLluPVvLIa4fYeeI8ACMz08hMT+FcTT3rDnn4yYYjjMxM474F41h+7XASEsTlGhvToUBt1LBgx0UkG/gGcIuq1orIV4D7gYc7k6+IZDifvy0i8/AGSPepalnbisViG7LzRDXfe/UQ245VATAisy+Z6amUXahnfYmHn248yrCMvnxh/lhWzBxJorUhJkzsUQ0dUFUVEQ3y3RPAEwCFhYUBz4kWdY0tfOvP+3lm+0mGDuzD//vAZJZNH0p2v9Qr59Q2NLOuuJwnNh7l357bw592nOK/PnI1IzLTXKy5MSE3G++w0l+daXkpwFtdSJ8EDAc2q+r9InI/8F/Ap9qeGEttSH1TC9/9y0F+s+Udcvun8rXbJnLHNcPI7d/nyjl1jS2sLy7niTeP8vUX9vH8ztP890euJj873cWam1gVigDnNDDC7/Nw51igc06JSBIwEKjsIG2g45VAhogkOb04/ucHK6O9eg8PUnaZiAxR1bNOV3Z5O/lEvXM19Xz2V9s5ePYC984dyz8tHE/flMT3nZeemsQHrhrK30wbwh+LTvHtlw+w7PG/8sSnZlCYn+lCzY3plGDtzGneHYr2HV8PCN6hpzv9MxGR64CfOh8fbCd9JXAZeN45/kfenSsYkzwXG7j7qe3sOVXD524czZcWFZCe+v5/XvqmJLJ02hCWTB3MS7vO8M2V+1n2+F/58Sev5Yax2S7U3MSyUAyCbgfGO6ubUvBOGl7Z5pyVgG+FwXJgrTMevhJY4ayAGo13EuC2YHk6adY5eeDk+VIHZQTkDEFdEJHZznygTwfJy7+MmHOi8jLLf7KZk1WX+eVnZ/LA0okBgxt/IsJHZ47gpftuZGDfZD7x861sOhx3UwxM9FgJfNpZTTUbqHGu/9eBRSIyyJkcvMg5tgW40TdvT0TSRaRAVbeq6nTntTJYeqfd+TPvBj8Lee+cxJhy5nwdH/3pWxwuu8QTn5rBgx+cHDC48Sci3HHNMP58303k9k/lrl9sY83B943gGdMzqtrjF3AbUIJ3rPnrzrGHgdud933w3sWU4g1gxvil/bqTrhhY2l6ezvExTh6lTp6pnSjjOFAFXMI7Tj7ZOV4I7HPK+CEgzvEsvKunDgOrgcyO/gYzZszQaFN2oU5vfmStTv/W67r7ZHW38qi81KCLH92gE77xim49WhniGppYBBRpCNodfff6/j1wFmhyru+7gXuBe53vBe+qzCPAXqDQL+3nnDajFO9EYt/xBXhvtPY4r9uDlB0s/Shgo5N2DTCyo98RjW1I5aUGnf9f63TqN1/TouPdu/7PX27UD/7fmzr+a6/oxpLyENfQxKLOtiG+f9BNDxUWFmpRUZHb1ei0usYWlv9kM8cqann687OZPiKj40RBVF5q4CM/fYuq2kZe/MKNNp5u2iUiO1S10O16RJpoa0MamltY8cQWDpy5wG//7jpm9mCYuuZyEx974i1OV9fx/BduYHxe/xDW1MSazrYhtk4vDqkqX31+DwfOXuCHH7+mR8ENQFa/VH75mZkI8He/LuJyY3OHaYwx0e2hlft5+8R5fvCx6T0KbgAGpiXz5GdmkpqcyN1PFXGhvilEtTTxzAKcOPT0thO8uOsM999SwIKJAVfAd9morHQe//i1HPFc4sGX9ockT2NMZHp+5yl+v+0kX5g3lqXThoQkz2EZffnJJ6/l9Pk6vvr8Xmx0wfSUBThx5qjnEt95+SA3j8/mi/PHdZygC24Yl80/LhjPcztO8efdZ0KatzEmMpysusyDL+1nZv4gvrRoQkjzLszP5EuLCvjLnrP8ccepjhMY0w4LcOJIS6vypT/uJiUpgf9cfnVYNun7p4XjmT4igwdf2kfFpYaQ52+McY+q8m/P7QHgfz46PSyb9N07ZyzXjc7k2y8f4FxNfcjzN/HDApw48vTWd3j7xHkeun0ygwf26ThBNyQmCP+5/CpqG1p4aKUNVRkTS57bcYq3jlbytdsmhW2Dz4QE4ZG/vYqmlla+8eLesJRh4oMFOHGi/EI933+tmJvGZXPH9GEdJ+iB8Xn9uW/BOF7ec9b2xzEmRlTVNvLdVw5SOGoQK2aO6DhBD+Rnp3P/rQWsPlhu++OYbrMAJ058//ViGppb+fYdU3G2nw+re+aMYWRmGt9cuY/G5tawl2eMCa//WVXMxfpm/v3D03rlGXSfvXE043L78a0/H6C+qSXs5ZnYYwFOHNh7qobndpziszflM7qX9qjpk5zIgx+YzBFPLc9uP9ErZRpjwqP43EWe3nqCT143koJe2qMmOTGBb90+hRNVl/nNW+/0SpkmtliAEwe++8oBsvulcF+IV011ZOGkXGblZ/LY2lLqGu0OzJho9R+vHqRfahL/fEtBr5Z747hs5hTk8KP1pVy0vXFMF1mAE+M2H6lgy9Eqvjh/HP37JPdq2SLCvy6egOdiA0+9dbxXyzbGhMbOE9WsL/Zw77yxDEpP6fXyv7xoAtWXm3hy07FeL9tENwtwYpiq8oNVh8kbkMqds0a6UodZozOZNyGHH68/YruTGhOFfrD6MJnpKdx1fb4r5U8bPpClUwfz8zePUVXb6EodTHSyACeGbT5SybbjVXxh3jj6JLf/hPBw+tdFE6ipa+JnG4+6VgdjTNfteKeajSUe7pkzpsMnhIfT/bcWcLmxmR+vL3WtDib6WIATo1SVR1eVMHhAHz4W5iWdHZk6bCB/M20IT26yOzBjoskPVpeQlZ7Cp68f5Wo9xuf150PXDOept96h7IJt/mc6xwKcGLWptIKid6r5wvyxrvbe+PzTLeO53NhiqyGMiRJFx6t483AFfz93DGkp7vXe+Px/C8fR3NLKrzYfd7sqJkpYgBOjHl9XGhG9Nz4Fef2ZPyGHX7913Pa0MCYKPL6ulOx+KXxytru9Nz6jstJZMnUwv9vyDpcamt2ujokCFuDEoP1nathytIrP3phPapL7vTc+98wZS2VtI8/vPO12VYwx7Sgtv8S6Yg+fmp0fEb03Pp+/eQwX6pv5w/aTblfFRAELcGLQLzYdJy0lkRUz3Vk5FczsMZlMGzaQn795lNZWdbs6xpggfrX5GClJCXxidmS1IdeMHMSs/Eye3HSM5hbbId20LyQBjogsEZFiESkVkQcCfJ8qIs86328VkXy/777qHC8WkcUd5Skio508Sp08U7pThohMEJFdfq8LIvLPzncPichpv+9uC8XfqTeUX6znz7vPsHzGcAam9e6+Nx0RET4/ZwxHK2pZc6jc7eqYGNHetex3zjwRqfE750G/79ptvzpRfrC26k2/8s6IyIs9+6W94/zlRv604zR3TB9Kdr9Ut6vzPp+fM4bT5+t4dd85t6tiIlyPAxwRSQQeB5YCk4E7RWRym9PuBqpVdRzwKPCIk3YysAKYAiwBfiQiiR3k+QjwqJNXtZN3l8tQ1WJVna6q04EZwGXgBb86P+r7XlVf6enfqbf8dssJGlta+eyNo92uSkC3TR3MsIy+PLHxiNtVMTGiE9eyz5t+1/TD0On2K6j20qvqzX71egt4vvu/svc8ve0EdU0tfO6myGxDFk7MZUx2Ok9sPIqq9QSb4ELRgzMLKFXVo6raCDwDLGtzzjLgKef9c8BC8T7xcRnwjKo2qOoxoNTJL2CeTpoFTh44ed7RzTL8LQSOqGpUL/Gpb2rhd1veYeHE3F575lRXJSUm8LmbRrP9eDX7Tte4XR0Te7p6LQdtv0RkhohsEJEdIvK6iAzpSnofERmAt92K+B6cppZWfr35HW4cl8XEwQPcrk5ACQnC3TePZu/pGnaeqHa7OiaChSLAGQb4z/g65RwLeI6qNgM1QFY7aYMdzwLOO3m0LaurZfhbAfy+zbH7RGSPiPxCRAYF+uGRZuXuM1TWNnJ3hN55+Sy/djipSQk8vc0ewmlCLtC17HO9iOwWkVdFZIpzLGD7ICLJwP8By1V1BvAL4LsB8uxM+3IHsEZVL3Ttp/S+V/ed49yFej4XoT3APndMH0a/1CR+t9XaEBNc3E8ydubw3A780e/wj4GxwHTgLPDfQdLeIyJFIlLk8XjCXteO/G7LOxTk9eP6sVluV6VdA9OS+eDVQ3np7dO23NOETJBr2WcnMEpVr8YbuHTUmzIBmAqsEpFdwDeA4d2s2p0ECboisQ0ZmZnG/Am5blelXempSdxxzVD+sucs5y/b5qEmsFAEOKcB/81WhjvHAp4jIknAQKCynbTBjlcCGU4ebcvqahk+S4GdqlrmO6CqZaraoqqtwM94/5CW77wnVLVQVQtzcnICndJrDpy5wO5TNayYORLvyFxk+/h1I6ltbGHlrjNuV8XEjvddyz6qekFVLznvXwGSRSSb4O2DAPv95uxMU9VFIjLCb+Lwve2kB8ApYxbwl0AVjqQ25KjnEluPVfGxmSNISIiCNmTWKBqaW/mTbTthgghFgLMdGO+sbkrB20W8ss05K4G7nPfLgbXqnR22EljhrIAaDYwHtgXL00mzzskDJ8+XulmGz/vurtqMtX8I2Nelv4gLnt1+gpSkBD58bdve8ch0zYgMJg7uz++2vmMTBU2otNdTMtiZk4eIzMLb9lUSvP0qBnJE5HonTbKITFHVk35Bz0/aSe+zHHhZVSP++QLPbj9JYoLwkRnd7ajqXZOHDuCakRk8bW2ICaLHAY4z3+U+4HXgIPAHVd0vIg+LyO3OaU8CWSJSCtwPPOCk3Q/8ATgAvAZ80ek5CZink9dXgPudvLKcvLtcBoCIpAO38v7VDd8Xkb0isgeYD/xLT/9O4VTf1MILb59m6dTBZKSluF2dThERPjF7FPvPXGDPKZtsbHom0LUsIvc6vSzgDTT2ichu4DFghXoFbGucCcPLgUecNLuAG9qW20FbBe3PCYoYjc2t/GnnKRZOzCV3QB+3q9NpH581kiOeWrYeq3K7KiYCiUW+oVFYWKhFRUWulP38zlPc/4fdPP3567hhbLYrdeiOi/VNXPfva/jgVUN5ZPlVblfH9BIR2aGqhW7XI9K42Ya8uvcs//C7nfzyMzOZPzGy59/4q2ts4bp/X828Cbk8duc1blfH9JLOtiFxP8k4Fjyz7ST5WWlcPyayJxe31b9PMrdfPZSVu8/YZGNjXPT77ScZMrAPcwrcnQfUVX1TEvnwtcN5bd85qmttsrF5Lwtwolxp+SW2Ha/iY1EyubitjxQOp66phddsV1JjXHGy6jJvHvbwkcIRJEbB5OK2PlI4nMaWVl7ee9btqpgIYwFOlHtuxykSE4S/nREdk4vbunbkIEZlpfH8zlNuV8WYuOR7+O1HC6NjcnFbk4cMYOLg/taGmPexACeKtbYqK3edZs74bHL7R8/EQH8iwoevGc5bRys5c77O7eoYE1dUlRd3nWb26CyGD0pzuzrdIiJ8+NphvH3iPEc9l9yujokgFuBEsW3HqzhTU88d10Rn743Ph64Zhiq8uMv2szCmN+0+VcOxilo+FOVtyLLpw0gQePFta0PMuyzAiWIvvn2a9JREFk0e7HZVemRkVhqz8jN5fudp28/CmF704tunSUlKYMm06G5D8gb04abxOTz/9mlaW60NMV4W4ESp+qYW/rL3LIunDqZvSqLb1emxD107jNLyS+y1B3Aa0yuaWlr58+4z3DopjwF9kt2uTo99+JphnKquY/tx2xPHeFmAE6XWHSrnYn1z1Hct+9w2bQgpSQlXJjwaY8Jr0+EKKmsbo36I22fRlDzSUxKtDTFXWIATpV54+zS5/VOjamO/9gzsm8ytk/NYufsMzS2tblfHmJj3wtunyUhLZm6U7X0TTFpKEkunDeGVvWepb2pxuzomAliAE4XOX25kXXE5y6YPjcp9K4L54FVDqapt5K2jlW5XxZiYdqmhmTcOnOMDV3l7TmPFB68eysWGZjaWuP9kduO+2Pk/O468uu8cTS3Ksumx0bXsM29CDv1Sk3h5t23YZUw4rTpwjvqmVu6IsTbkhrFZDEpL5uU91oYYC3Ci0l/2nGV0djpThg5wuyoh1Sc5kVsn5/Ha/nM0NtswlTHh8pc9ZxkysA/XjhzkdlVCKjkxgSVTB7P6YBl1jTZMFe8swIkylZcaeOtoJbdNGxyVj2boyAeuGkJNXRN/La1wuyrGxKQL9U1sLKngtmlDSIihIW6fD1w1lMuNLawrLne7KsZlFuBEmdf3l9HSqvzNtKFuVyUsbh6fw4A+Sfx5zxm3q2JMTFp9oIzGllZumzbE7aqExXWjM8nul8LL1obEPQtwoswre73DU5OG9He7KmGRkpTA4imDWbW/zFZCGBMGr+w9y9CBfbhmRIbbVQmLpMQElk4dwtpD5dQ2NLtdHeMiC3CiSOWlBjYfqYjZ4SmfD9hKCGPCwjc8tTRGh6d8PnDVEOqbWllzyIap4pkFOFHk9f1ltCoxOzzlYyshjAkP3/DU31wVm8NTPjPzM8kbkMrLu22YKp5ZgBNF/rL3TEwPT/kkJ3qHqdYeKqeh2YapjAmVWB+e8klIEJZOHcKGEg+XG22YKl6FJMARkSUiUiwipSLyQIDvU0XkWef7rSKS7/fdV53jxSKyuKM8RWS0k0epk2dKD8o4LiJ7RWSXiBT5Hc8UkVUictj5r+trKSsvNfDWkdhdPdXW4qmDudTQzOZS2/TPdCzYtez3vYjIY047sEdErvX77i7nWj8sInd1o+z3pReR/k5dfK8KEflBz35lz/gPT8VFGzJlMA3NrWwotqHueNXjAEdEEoHHgaXAZOBOEZnc5rS7gWpVHQc8CjzipJ0MrACmAEuAH4lIYgd5PgI86uRV7eTd5TL86jZfVaeraqHfsQeANao6HljjfHbV6oPe4amlU2O7a9nnhrFZ9EtN4rV959yuiokega5ln6XAeOd1D/Bj8N7MAN8ErgNmAd/syg1NsPSqetGpy3RVnQ68Azzfg9/WY+sOlTurp6L7yeGdNTN/EIPSknltv7Uh8SoUPTizgFJVPaqqjcAzwLI25ywDnnLePwcsFO8txDLgGVVtUNVjQKmTX8A8nTQLnDxw8ryjm2W0xz8v/zJc88b+MoZl9I25zf2CSU1KZMHEXFYd9C6LN6aHlgG/Vq8tQIaIDAEWA6tUtUpVq4FVeG+EEJFFIvKWiOwUkT+KSL8A+QZN7yMiBUAu8Gb4fl7H3jhQRna/VK4Z4XqHdK9ISkzg1sl5rD1YbhuHxqlQBDjDgJN+n085xwKeo6rNQA2Q1U7aYMezgPNOHm3L6moZAAq8ISI7ROQev3PyVNU3w/UckBfoh4vIPSJSJCJFHk/4ukFrG5p5s7SCWyfnxUXXss/iKYOpqm1k+/Eqt6tiIl+wa9mnS22NiGQD3wBuUdVrgSLg/i7k628F8Kyqvi9S7602pKG5hQ3FHm6dnBvTq6faWjxlMBcbmtl8xDYOjUfxPsn4JqfxWgp8UUTmtD3BaZQCdiGo6hOqWqiqhTk54Xsi75uHPTQ2t7JoSsA4K2bNm5BDSlICr1sXs+lYh9dyF83GOzz+VxHZBdwFjOpmXiuA3wf6orfakLeOVHKpoZlFk+NjeMrnxnHZpKck8vr+MrerYlwQigDnNDDC7/Nw51jAc0QkCRgIVLaTNtjxSrxdy0kByupqGaiq77/lwAu8O3RV5nRf4/zX1c0U3jhQxsC+yczKz3SzGr0uPTWJOeOzeWN/GQFufo25op1r2aerbY3gHXryzaOZrKp3i8h1fhOHb28nPQAicjWQpKo7QvJDu2nVgTLSUhK5fmyWm9XodX2SE5k3MZdVB87ZUHccCkWAsx0Y76xuSsF7t7KyzTkr8d4BASwH1jo9IyuBFc4KqNF4JwBuC5ank2adkwdOni91pwwRSReR/gAikg4sAvYFyMu/jF7X3NLK2kPlLJyYS1Ji/HW4LZ4ymNPn69h3+oLbVTERqoNr2Wcl8GlnNdVsoMYZhn4dWCQig5zJxYucY1uAG0VknF8ZBaq61S/oWdlOep87CdJ701taW5VVB8qYW5BDn+TEjhPEmCVTBlNxqZGdJ6rdrorpZUkdn9I+VW0WkfvwXtSJwC9Udb+IPAwUOY3Ak8BvRKQUqMIbsOCc9wfgANAMfFFVWwAC5ekU+RXgGRH5DvC2kzddLUNE8oAXnDktScDTqvqak9f3gD+IyN14Vz98tKd/p+7afrya85ebuHVyfA1P+dwyKY/EBOG1/WeZNnyg29UxkSngtSwi9wKo6k+AV4Db8C4yuAx81vmuSkS+jfemCuBhVa0CEJHPAL8XkVTnu28AJf4Ft5dUsdAbAAAgAElEQVTe8VGnXNfsOV1D+cWGuBvi9pk/MZeUxARe23eOmXHWCx7vxLr+Q6OwsFCLit63/UaPPfznA/x26zu8/f9uJT21x/FoVLrziS1U1Tby+r/0dFqFiQQisiPIUu64Fq425PuvHeKnG4+y4xu3kJGWEvL8o8Fnf7mNI55aNnx5Xlwt1IhVnW1D4m/MI4qoKm8cOMdN47LjNrgBWDgpl+Kyi5ysuux2VYyJOqsOlHHd6My4DW4AFk7K40TVZY54LrldFdOLLMCJYIfOXeRUdR2L4nR4ymfhJO/vX3PQVkIY0xXHKmo5XH4pboe4fRZMzAVg9UF7+GY8sQAngq06UIbIu//Ax6vR2emMyUm3JwMb00WrDni3WIj3AGdoRl8mDxlgN0lxxgKcCPbGgXNcO3IQOf1TOz45xt0yKY8tR717eRhjOueN/WVMHjKA4YPS3K6K626ZlMuOd6qprm10uyqml1iAE6HOOEuj4/3Oy2fBxFyaWpQ3S+zBecZ0RsWlBnacqLY2xLFwUh6tCuuKrSc4XliAE6FWHfB2pcb7/BufwlGDGNg32cbQjemkNQfLUCVul4e3NW3YQHL6p9pQdxyxACdCrT5YxpicdMbkBHq+X/xJSkxg3oQc1hWX246kxnTC6oPlDHPmnhhISBAWTMhlY7HHHr4ZJyzAiUCXG5vZerSKhc7Mf+O1cFIeVbWN7Dp53u2qGBPRGppb+GtpBfMn5ti+L34WTsrlYkOzPcA3TliAE4H+WlpJY0sr8ydYgONv7vgcEhPEVkIY04Ftx6q43NhibUgbN43PJiUpgdXWhsQFC3Ai0NpD5fRLTaLQthV/j4FpyczMH8Qam4djTLvWHionNSmBG8Zmu12ViJKWksSNY7NYc7DcHuAbByzAiTCqyvricm4a573TMO91y6Q829XYmA6sL/Zw/dgs+qbE38M1O7LA2dW4tNx2NY519i9ohDl07iJna+qv7Lxp3sv3d1lrKyGMCehYRS3HKmpteCoI39xGW00V+yzAiTC+f7jnTchxuSaRaUxOP8Zkp9sYujFB+NoQu0kKzHY1jh8W4ESYdYfKmTpsALkD+rhdlYi1cFIuW45WUmu7GhvzPuuLyxmX248RmbZ7cTC2q3F8sAAnglTXNrLzRDULrGu5XfMneHc13nyk0u2qGBNRahu8W0zMtx7gds2bmEurwpulFW5XxYSRBTgRZONhD63qvfhMcDPyB5GWksh623LdmPfYVFrh3WLC2pB2XT08g4y0ZGtDYpwFOBFk3aFyMtNTuHp4httViWipSYncMDabDSUeW+ppjJ/1xd4tJmbaFhPtSkwQbh6fw8aSClptZ/SYFZIAR0SWiEixiJSKyAMBvk8VkWed77eKSL7fd191jheLyOKO8hSR0U4epU6eKd0pQ0RGiMg6ETkgIvtF5J/8zn9IRE6LyC7ndVso/k7taWlVNpR4mFfg3czOtG/ehBxOVddxxFPrdlWMi9q7jv3OmSciNX7X84N+37XbdnWi/GDt1Jt+5Z0RkRe7/ys7R1VZd8jDzeOzSU60e9eOzCvIoeJSAwfOXnC7KiZMenwViEgi8DiwFJgM3Ckik9ucdjdQrarjgEeBR5y0k4EVwBRgCfAjEUnsIM9HgEedvKqdvLtcBtAMfElVJwOzgS+2qfejqjrdeb3S079TR3adPE/15SYbnuqkuQXeOQbWxRz3OrqOfd70u54fhk63XUG1l15Vb/aVB7wFPN/9n9g5B85e4NyFehue6qQ51obEvFCE+bOAUlU9qqqNwDPAsjbnLAOect4/BywU7wNSlgHPqGqDqh4DSp38AubppFng5IGT5x3dKUNVz6rqTgBVvQgcBIaF4O/RLesOlZMg3scRmI6NyExjbE46G0o8blfFuKiH13HQtktEZojIBhHZISKvi8iQrqT3EZEBeNussPfgrC/2Xgu2xUTn5PRPZeqwAdaGxLBQBDjDgJN+n0/x/gbmyjmq2gzUAFntpA12PAs47+TRtqyulnGFM5x1DbDV7/B9IrJHRH4hIoOC/fhQWXuonBmjBjEwLTncRcWMeRNy2Xq0isuNtlzcBL2Ofa4Xkd0i8qqITHGOBWwbRCQZ+D9guarOAH4BfDdAnp1p++4A1qhq2MdB1h4qZ9qwgeT2ty0mOmteQS47T5ynpq7J7aqYMIj7gVoR6Qf8Cfhnv0box8BYYDpwFvjvIGnvEZEiESnyeLp/F1B+oZ4DZy8wz5aHd8m8CTk0trSy5agtF493Qa5jn53AKFW9Gm/g0lFvygRgKrBKRHYB3wCGd7NqdwK/b6feIWlDai438faJalse3kXzJuTQ0qpsOmzLxWNRKAKc08AIv8/DnWMBzxGRJGAgUNlO2mDHK4EMJ4+2ZXW1DJw7tT8Bv1PVK2Pkqlqmqi2q2gr8DG9X9Puo6hOqWqiqhTk53W9YNjoXl29eiemcmfmZ9E1OvNI1b+JTsOvYR1UvqOol5/0rQLKIZBO8bRBgv9+cnWmqusiZ0OybOHxvO+l99crG23b8JVjdQ9WGbCqtoFVhrgU4XTJ9RAYD+iSxocTm4cSiUAQ424HxzuqmFLwTele2OWclcJfzfjmwVr3re1cCK5wVUKOB8cC2YHk6adY5eeDk+VJ3ynDm5zwJHFTV//GvbJvx9g8B+7r1l+mkDSUesvulMnnIgHAWE3P6JCdy/dgsG0OPY+1dx37nDHbOQ0Rm4W33KgnedhUDOSJyvZMmWUSmqOpJv6DnJ+2k91kOvKyq9eH47f42lJQzoE+SbTHRRUmJCdw8Pse2nIhRSR2f0j5VbRaR+4DXgUTgF6q6X0QeBopUdSXeBug3IlIKVOFtCHDO+wNwAO9qiC+qagtAoDydIr8CPCMi3wHedvKmq2WIyE3Ap4C9Tjc0wNecO7zvi8h0QIHjwN/39O8UjLd71MP8Cbkk2PLwLps3IYe1h8o5VlHL6Ox0t6tjet+NBLiOgZEATiCyHPgHEWkG6oAVzs1PwLYLQESWA4+JyEC87eQPgP34Cdb2+Z2yAvheGH7ze6gqG0squGl8Nkm2PLzL5k7I4S97z3Lw7EUmD7WbzFjS4wAHrnT7vtLm2IN+7+uBjwRJ+10CTOALlKdz/CgBhoy6WoaqbsLbFR3o/E8FOh4Oe0/XUH25ybqWu8l/ufjo7NEu18b0tvauY79zfgj8MMh3wdqZXcCcTpQfML3z3byO0odCSdklzl2otyHubvL93TaUeCzAiTEW7rtsY4kHEbhpXLbbVYlKo7LSGZ1ty8VN/Nro/L8/xwKcbskb0IdJQwbYfjgxyAIcl20o8TBt2ECy+qW6XZWoNbcgh7eOVFLf1OJ2VYzpdRtKPBTk9WPIwL5uVyVqzZuQw453qrlYb8vFY4kFOC6qqWti18nzzLHN/Xpk3oQcGpptubiJP5cbm9l2rMrakB6aV5BDc6vy11JrQ2KJBTgu2lxaQUur2vybHpo9JovUpARbLm7iztajVTS2tFob0kPXjhpE/1RbLh5rLMBx0YYSD/1Tk5g+wpZ29kSf5ERmj8m6MhfBmHixocRDn+QEe3p4DyUnJnDjuGw2FNty8VhiAY5LvEs7Pdw4zp78GwrzJuRwtKKWE5WX3a6KMb1mY4mH2WOy6JOc6HZVot68CTmcqanncPklt6tiQsT+ZXVJafklztTU28qHELmyXNy6mE2cOFl1maMVtTb/JkR8w3y2mip2WIDjkg1Xlnba8vBQGJ2dzsjMNJuHY+KGrw2x+TehMWRgXybk9bc2JIZYgOOSDSUexuakM3xQmttViQkiwrwJ3uXiDc22XNzEvo0lHoZl9GWM7eAdMvMm5FB0vJrahma3q2JCwAIcF9Q3tbDtWBVzC+zp4aE0tyCHuqYWio5Xu10VY8KqqaWVzUcqmTshB+cxWyYE5hbk0NhiW07ECgtwXLDlaCUNza02PBVis8dkkZwotquxiXk736nmUkOzzb8JsRn5g+ibnGhtSIywAMcFG0sqSE1KYPaYLLerElPSU5OYmZ9py8VNzNtQ4iEpQbhhnLUhoZSalMgNY23LiVhhAY4LNpSUM2t0pi3tDIO5BTkcOneRczX1blfFmLDZeNjDtSMHMaBPsttViTlzJ+RwvPIyxytq3a6K6SELcHrZqerLHPHU2pN/w8S37N7uwEys8lxsYN/pC7Z6Kkx8w34bD1sbEu0swOllG0sqACzACZOJg/uT2z+VDdY4mRi1qdTZYsLm34RFfnY6o7LS7CYpBliA08s2lngYOrAP43L7uV2VmCQizC3IYdPhCppbWt2ujjEht6HYQ1Z6ClOGDnC7KjFrbkEOm23LiahnAU4vampp5a+lFcwpsKWd4TR3Qg41dU3sPlXjdlWMCanWVuXNwxXcPD6bhARrQ8JlbkEOlxtb2GFbTkQ1C3B60a6T57nY0GzDU2F207hsEgRb6mlizv4zF6isbbT5N2FmW07EhpAEOCKyRESKRaRURB4I8H2qiDzrfL9VRPL9vvuqc7xYRBZ3lKeIjHbyKHXyTOmtMnpqQ7GHxAThhnG2/004ZaSlcPWIDBtDj1HhaG96UJdeKwu8KzABbrb5N2Hl23LCApzo1uMAR0QSgceBpcBk4E4RmdzmtLuBalUdBzwKPOKknQysAKYAS4AfiUhiB3k+Ajzq5FXt5N1bZfTIxsMerhmRwcC+trQz3OYW5LD71HmqaxvdrooJoXC0N50sN19E1gf4KuRltWdjSQVThw0gu19qT7MyHfBtOVF2wbaciFah6MGZBZSq6lFVbQSeAZa1OWcZ8JTz/jlgoXgnoSwDnlHVBlU9BpQ6+QXM00mzwMkDJ887erGMbqu81MDe0zX29PBeMqcgB1V4s7TC7arEja+/sJc/bD8Z7mLC0d4gIp8UkW0isktEftqFYKTLZXXXhfomdpyottVTvcTXVlsvTu95aOV+frPlnZDlF4oAZxjg36qdco4FPEdVm4EaIKudtMGOZwHnnTzaltUbZbyHiNwjIkUiUuTxtH8RnKquY/igvjb/ppdcPdzbU2bDVL2jpq6JZ7af5GT15XAXFfL2RkQmAR8DblTV6UAL8Imu1qeT7c57dKUNOV1dxwhrQ3qNb8sJa0N6R11jC09vO8E7IdxgMSlkOcUhVX0CeAKgsLBQ2zv36hEZbPzy/F6pl4HEBOHm8dlsKPGgqrZqLcw2l1bQ0qrR2kO5EJgBbHf+P+kLlAOIyAvAaCAFGCkiu5w0/6uqv+xpwV1pQyYNGcD6L89Htd3TTIj4tpx440AZLa1Koq1aC6stxyppbG4NaRsSih6c08AIv8/DnWMBzxGRJGAgUNlO2mDHK4EMJ4+2ZfVGGT0iIvYPbS+aW5CD52IDB89edLsqMW9DiYf+qUlMH5ER7qLC0d4I8JSqTndeE1T1IQBV/ZDTq3MbUOR3zi+7WVaPWRvSe97dcuK821WJeRuKPaQmJTBrdGbI8gxFgLMdGO+sPErBO7FuZZtzVgJ3Oe+XA2vVexuyEljhrEQYDYwHtgXL00mzzskDJ8+XerEME0VsDL13qCobSzzcMC6L5MSw7zwRjvZmDbBcRHIBRCRTREZ1sj5dLctEkStbThRbGxJuGw97mD0mK6TPaOxxa+SMO98HvA4cBP6gqvtF5GERud057UkgS0RKgfuBB5y0+4E/AAeA14AvqmpLsDydvL4C3O/kleXk3VtlmCiSN6APEwf3tzH0MDviucSZmnrmFuSGvawwtTcHgG8Ab4jIHmAVMKSTVepSWT379aa3Xdlywh79Elanqi9z1FMb8iFusfHc0CgsLNSioiK3q2Ha+I9XD/KLTcd4+8FF9Eu1KWfh8OSmY3z75QNs+sp8hg9K6/B8EdmhqoW9ULWoYm1IZPrB6hIeW3OYHd+4lUHpIdkSzbTx9NYTfO2Fvay+fw7jcvt3eH5n2xDbydjEtLnjc2hqUd46Uul2VWLWhhIPY3LSOxXcGBNt5hTk0KqwybacCJsNJeUMHdiHsTmhfUajBTgmps3IH0RaSqINU4VJfVMLW49W2t4sJmbZlhPh1dTSyubSSuZOCP0zGi3AMTEtNSmRG8Zmsb6k3JbXhsG2Y1U0NLfa3iwmZrXdcsKElu8ZjeG4SbIAx8S8uQU5nKyq43hl2DehizsbSzykJCVw3ZjQLe00JtLMLcih/GIDh87ZlhOhtrEkfM9otADHxLwry8WLy12uSezZUOJhVn4maSk2gdvELttyInw2lHiYHqZnNFqAY2LeqKx08rPS2HjYJgmG0pnzdRwuv8ScgtDfeRkTSWzLifCoqm1k7+masA1xW4Bj4sLcghzeOlJJfZNtRRIqbzp7g0Tp4xmM6ZK5E3LYfryK2obmjk82nfLmYQ+q4WtDLMAxcWFOQQ51TS0UHa92uyoxY2NJBXkDUpmQ1/G+FcZEO9tyIvQ2llSQkZbMtGEDw5K/BTgmLswek0VKYoLtSBoizS2tbCqtYM740C/tNCYSXdlywtqQkFBVNh72cNO47LA9yNQCHBMX0lOTmDl6kD1TJkR2n6qhpq7JhqdM3PBtOWETjUPj4NmLeC42hHWLCQtwTNyYW5BDcdlFztbUuV2VqLexxIOI92GExsSLuQU5vFN5meMVtW5XJept7IU5fBbgmLjhu5BsJUTPbTzs4arhGfZsHhNXbLl46Gws8TBxcH/yBvQJWxkW4Ji4MSGvP3kDUtlYYsvFe+L85UZ2nzzP3PHWe2Piy5UtJyzA6ZHahmaKjleHfYjbAhwTN0SEuQU5vHnYQ3NLq9vViVqbSitoVe+yWWPizdyCHDYfqaSh2bac6K4tRytpbAn/I14swDFxZU5BDhfqm9l96rzbVYlaG0s89O+TxNXDM9yuijG9zrac6LmNJR76JidSmD8orOVYgGPiyk3jskkQ2GDDVN2iqmwo8XDz+GySEq35MPHnypYTNkzVbetLPFw/NovUpMSwlmMtlIkrGWkpTB+RYZMEu+ng2YuUXWhg3oRct6tijCuubDlhbUi3HKuo5Z3Ky8zrhSHuHgU4IpIpIqtE5LDz34D9TSJyl3POYRG5y+/4DBHZKyKlIvKYODuGBctXvB5zzt8jItf2oIz/FJFDTj4viEiGczxfROpEZJfz+klP/kYm8swtyGXPqfNU1Ta6XZWos855YOm8CNn/RkQ+4VzDe0Vks4hcHeS8X4nIMb/rerpzPGib0snyA6YXkfl+Ze0SkXoRuaPnv9hEgrkFORw6d5FzNfVuVyXqrDvka0PCf5PU0x6cB4A1qjoeWON8fg8RyQS+CVwHzAK+6RcI/Rj4PDDeeS3pIN+lfufe46TvbhmrgKmqehVQAnzVr9pHVHW687q3G38XE8HmFGSj+u6zlEznbSj2MGXoAHLDuLSzi44Bc1V1GvBt4Il2zv2y33W9yzkWsE3pgoDpVXWdryxgAXAZeKOLeZsIZVtOdN/6Eg9jc9IZmZUW9rJ6GuAsA55y3j8FBLpDWQysUtUqVa3GG1gsEZEhwABV3aKqCvzaL32wfJcBv1avLUCGk0+Xy1DVN1TV99S0LcDwHv4tTJS4angGGWnJtly8i2rqmthxorpXupY7S1U3O9c8dO86DtamICJfFpHtTs/Mt7qa3s9y4FVVvdzFupkI5dtyYoPdJHVJXWMLW45W9toQd08DnDxVPeu8PwfkBThnGHDS7/Mp59gw533b4+3l215eXS3D3+eAV/0+jxaRt0Vkg4jcHOB8E8USE4Sbx+ewocRDa6u6XZ2oselwBS2tyvzInX9zN++9jtv6rhOsPCoiqc6xgG2HiCzC2yszC5gOzBCROQHyDNb2+FsB/L7zP8NEOt+WE5sOV9iWE13w1tEKGptbe60N6TDAEZHVIrIvwGuZ/3lOD0nI/7UIV74+IvJ1oBn4nXPoLDBSVa8B7geeFpEBQdLeIyJFIlLk8VgkH03mjM+m4lIDB89dcLsqUWNdcTkD+yYzfUTkLQ8Xkfl4A5yvBDnlq8BEYCaQ2c55Pouc19vATift+G7UawgwDXg9yPfWhkSpOQU51NQ1sftUjdtViRrrDnlIS0lk5ujwLg/36TDAUdVbVHVqgNdLQJlfd+4QoDxAFqeBEX6fhzvHTvPe7mTfcdrJt728uloGIvIZ4APAJ5xAClVtUNVK5/0O4AhQEORv84SqFqpqYU5O5HTbm47NvTKGbsNUndHaGjnLw0Xki36Td4eKyFXAz4Flvmu3LVU96wwjNQC/xNszA8HbDgH+w2/OzjhVfbJt2e2k9/ko8IKqNgWpl7UhUcq35YTNw+kcVWVdcTk3jM0O+/Jwn562VCsB34qlu4CXApzzOrBIRAY5E38XAa87Q1AXRGS2s7Lp037pg+W7Evi0s3JhNlDj5NPlMkRkCfBvwO3+Y+MikiMiic77MXjv2o725I9kIk/ugD5MGjKADSWBYnLT1oGzF/BcbIiI4SlVfdxvAm8S8DzwKVUtCZbG74ZJ8M7D2+d81V6b8jkR6eekGyYiuf5lq+qZdtL73IkNT8Uk23Kia454ajlVXderc/h6GuB8D7hVRA4DtzifEZFCEfk5gKpW4V3dsN15PewcA/gC3juvUrw9Ja+2ly/wCt5goxT4mZO+u2X8EOgPrJL3LgefA+wRkV3Ac8C9fnmZGDK3IIei49Vcamju+OQ4t95ZHh7uZ8d0w4NAFvAj5zou8n0hIq84vSwAvxORvcBeIBv4jnM8WJvyBvA08JaT7jm87UVbAdM75efj7d3ZEIofaiLP3IJcdp86T7VtOdEhXxvSmwGOOCMzpocKCwu1qKio4xNNxNh8pIKP/2wrT3xqBoumDHa7OhHtb3+8maaWVlbed1OP8xKRHapaGIJqxRRrQ6LP2yeq+dCPNvPYnddw+9VDO04Qxz75862UXahn1f1ze5xXZ9sQ28nYxK3CUZmkpSSy0ZZ6tuv85UbePlFtuxcb08a7W05YG9Ke2oZmth2rYv7E3m1DLMAxcSslKYEbxmazvtiD9WQGt/Gw9+nhkbT/jTGRwLac6JzNR7xPD+/tHdAtwDFxbcHEXE5V13G4/JLbVYlY6w+VMygt2Z4ebkwA8yfk4LnYwL4ztlw8mHXF5aSnJFKYn9mr5VqAY+LaAqfLdPXBMpdrEpl8y8PnFOSQmCBuV8eYiDN/Qi4JAqsP2orMQFSVDcUebhyXTUpS74YcFuCYuDZ4YB+mDhvAGmucAtp96jyVtY0RsTzcmEg0KD2FGaMGscZukgIqLrvI6fN1vT7/BizAMYaFE/PYeaKayksNblcl4qw+WEZigliAY0w7Fk7KY/+ZC5ytqXO7KhFn9QFv4LfQAhxjet8tk/JQhfXFthKirTUHy5mZP4iBacluV8WYiOX7x9t6gt9v9cFyrh6RQe6APr1etgU4Ju5NHTaAvAGprDlkXcz+TlZd5tC5i9wyKdAzdI0xPuNy+zEyM421hyzA8Vd+sZ5dJ89ziwu9N2ABjjGICAsm5rKxxPukW+Plm3i90AIcY9olIiyclMtfSyuoa2xxuzoRY63To3XLZHfaEAtwjME7D+eSsxmV8VpzsJxxuf0YnZ3udlWMiXi3TMqjobmVTaX2AF+f1QfLGZbRl4mDAz3lJPwswDEGuHFcNqlJCbZc3HGhvoktRytZOMkmFxvTGTPzM+mfmmSrqRx1jS1sKvVwy6RcvM+47X0W4BgD9E1J5KZx2aw5VGa7GgMbij00tyq32vCUMZ2SkpTAnIIc1hwqt12Ngb+WVlDf1Ora8BRYgGPMFQsn5XGyynY1BlhzsIzM9BSuGTnI7aoYEzUWTsq1XY0daw6V0S81ietGZ7lWBwtwjHEssKWeADS1tLL2UDnzJ+Ta7sXGdIHtauzV2qqsPljO3IKcXt+92J8FOMY43t3VOL7H0IuOV3OhvplbbP6NMV1iuxp77Tldg+dig+tz+CzAMcaPb1fjqtpGt6vimjUHy0hJTODmXn7yrzGxYMFE29V4zcEyEgTXd0C3AMcYP7dOzqNVids7MFVl9cEyrh+bRb/UJLerY0zUuXWy8wDfA/HZhgCsOlBGYX4mg9JTXK2HBTjG+JkydADDMvry+v5zblfFFUc8tRyvvGzDU8Z007jc/ozJSee1OG1DTlX7dkB3vw3pUYAjIpkiskpEDjv/DbjkQkTucs45LCJ3+R2fISJ7RaRURB4TZ7F8sHzF6zHn/D0icm0PynhIRE6LyC7ndZtfmq865xeLyOKe/I1MdBERFk8ZzMbDFdQ2NLtdnV636kB07V4sIvNEpMbvOn4wyHmjRWSrc10/KyIpzvFU53Op831+F8sPmF5EPuFXp10i0ioi03v2a020WDJlMFuOVnH+cvwNdUdSG9LTHpwHgDWqOh5Y43x+DxHJBL4JXAfMAr7pFwj9GPg8MN55Lekg36V+597jpO9uGQCPqup05/WKk9dkYAUwxTn3RyKS2K2/jolKi6fk0djcGpcP33xt31muGj6QoRl93a5KV7zpdx0/HOScR/Be7+OAauBu5/jdQLVz/FHnvK4ImF5Vf+erE/Ap4Jiq7upi3iZKLZ4ymBZnJVG8eW3fOcbn9mNsTj+3q9LjAGcZ8JTz/ingjgDnLAZWqWqVqlYDq4AlIjIEGKCqW9S7s9qv/dIHy3cZ8Gv12gJkOPl0p4z2ftMzqtqgqseAUrxBk4kThfmZZKWnxN0w1enzdew+VcOSqYPdrkpIOb22C4DnnENt2xRfW/McsNDpKU4Ukf8Uke1Ob/HfB8k+YPo259wJPBOK32Kiw1XDBzJkYJ+4a0MqLjWw/XgVSyOkDelpgJOnqmed9+eAQH1Sw4CTfp9POceGOe/bHm8v3/by6moZAPc5jdcv/Hp8guX1PiJyj4gUiUiRxxN/d/uxKjFBuHVyHmsPldPQHD8Pznt9n7cxXjp1iMs16bLrRWS3iLwqIlMCfJ8FnFdV35ij/zV95Xp3vq9xzr8bqFHVmcBM4PMiMiiox44AABt5SURBVDpA3sHS+/sY8PtAFbc2JDZdGeou8cTVUPcb+8toVVgSIW1IhwGOiKwWkX0BXsv8z3N6SEK+P3W48sU7dDUWmA6cBf67qxmo6hOqWqiqhTk5tqQ2liyeOphLDc1sLq10uyq95rV955g4uH+0PVxzJzBKVa8G/g94MUT5LgI+LSK7gK14g5bxXc1ERK4DLqvqvkDfWxsSuxZPGUxDcysbSuIncH1131lGZaUxaYg7D9dsq8MAR1VvUdWpAV4vAWXOMBDOfwMNOJ4GRvh9Hu4cO+28b3ucdvJtL68ulaGqZaraoqqtwM94dxgqWF4mjtwwNov+qUlx08VcfrGe7e9UsXhKZHQtt0dEvuibvAv0U9VLAM48umQRyW6TpBLvcLZv3bv/NX3lene+H+icL8A/+s3tGa2qb4jId/3Kbi+9zwqC9N6Y2DYzfxCZcTTUXXO5ibeOVLJk6mDXHq7ZVk+HqFYCvhVLdwEvBTjndWCRiAxyhoEWAa87Q1AXRGS2M2b9ab/0wfJdifeuSkRkNt4u5LPdKcMXQDk+BOzzK2OFszpiNN67tm3d+/OYaJWalMj8ibmsOlBGSxw8OG/VgTJUYem0yA9wVPVxvwm8rX4rI2fhbdMq25yvwDpguXOobZvia2uWA2ud818H/kFEkp28C0QkXVW/7ld2e+kRkQTgo9j8m7iUlJjALZNyWXuwnMbmVrerE3arD5bR3KoRNcTd0wDne8CtInIYuMX5jIgUisjPAVS1Cvg2sN15PewcA/gC8HO8E3mPAK+2ly/wCnDUOf9nTvrulvF98S4f3wPMB/7FyWs/8AfgAPAa8EVVjZ+JGOaKJVMHU1nbyNZjsT9M9dq+c4zOTmdCXmR0LXfBcmCfiOwGHgNW+AUYr4jIUOe8rwD3i0gp3uGmJ53jTwJZzvH7eXfF5s/xtgE7RWQf8FMg0M6HwdIDzPn/27vz8KrqM4Hj3zcLCWEJBEIIWdgEwcgeFndGERBEcBlF6zKdVtvaaevYdqx15tFxmWo7o91sq4U6WkVxrbgii1KKgCasYQ9rFnZCQghkfeePe8Lc4s1G7s259+T9PM99vPfknN/vPdfnvrznd87vHKBAVXcFZ1dNpJl6YW9OVNbwt3zvn6b6KO8AqYnxjEhPdDuUM8TJBaaVsrOzNScnx+0wTBCdqqplzOOLmDUqjf+6fpjb4YTM8Yoqsh9fzN2XD+CBqUNC3p+I5Kpqdsg7ijCWQ7ynqqaO7McXMWloCk/f4t3bIJVX1jD6sUXcNi6TR64LdJ1/cDU3h9idjI1pQMcO0Vw1NIWP8w5QU+vdIeZPNvuGlqdGwPU3xkSSDjFRTMnqzSebD3K62rsnApZu9Z2GC5fp4fWswDGmEdcOT+XYySo+3+nd01TvrS8mMymB4WE0tGyMV1w7og/llTWenk313vpiUrrGMbZfktuh/B0rcIxpxBWDk+kSF8P7G4rdDiUkjpRXsiL/CDNGpIbNzAdjvOTigT3onhDL+xv2N71yBCo9Vc2ybYe5dngfoqLCK4dYgWNMI+Jjo7n6At9pKi/OhPhw437qFK4bEfBelsaYVoqNjmLqhaks2XKQU1XeO021cNMBqmrruG5En6ZXbmNW4BjThGtHpFJ22pszId5bX8zglM6c3zviZk8ZEzFmDE+loqqWpVu992yqcD7FbQWOMU249LxkusbH8P56bw0xFx0/xZd7SsLyyMsYLxk/oAc9O8d57lT3kfJKPt95NGxPcVuBY0wTvDoT4gMn2c6wAseYkIqOEqYN683SrYco99CzqT7auJ/aOg3bU9xW4BjTDLNGpVFeWcPiLQfdDiVoFqwvZkR6In17RNSzp4yJSDNHplFZU8dHG70zErwgzE9xW4FjTDNMGNCD1MR43l7jjceS7TpcTl5RmY3eGNNGRmd2o1+PBM/kkOIIOMVtBY4xzRAdJcwalcay7Yc5fKLS7XBa7e01RUSJnZ4ypq2ICDeMTmflrqMUHT/ldjit9s5aX6EWrqenwAocY5rthlFp1NYp762P7AsFa+uUt9YUcvngZFK6xrsdjjHtxvWjfMXAX9ZG9iiOqvJmbiHj+yeR2SPB7XAaZAWOMc00KKULw9MTeXttoduhtMrKnUfZX3qam8akux2KMe1KRlIC4/on8faaQiL5OZBr9pWw+8jJsM8hVuAY0wI3jEojr6iMbQdOuB3KOXsjt4Cu8TFMGpridijGtDs3jk5j5+GTbCgsdTuUc/ZGTiEJHaKZNizV7VAaZQWOMS0wY0QfYqIkYkdxyk5X83HeAWaOTCM+NtrtcIxpd64ZlkpcTBRvr4nMHHKqqpb3N+xn2rBUOsXFuB1Oo6zAMaYFenSOY+L5vXhnTVFEPmH8gw37qaypC/uhZWO8qmt8LJOzevPu+mIqayLvvloLNx2gvLImInKIFTjGtNDN2ekcOlHJkgi87fqbuYUM6tU5LG+rbkx7cXN2OscrfKOpkebN3EIykjoyLsyeHB6IFTjGtNCVQ3rRu2s881bvczuUFsk/VE7u3hJuHJMelrdVN6a9uGRgTzKTEnglwnJIwbEKVuw8wg2j0sPuyeGBtKrAEZEkEVkkIjuc/3ZvYL27nHV2iMhdfsvHiMhGEckXkV+Lk3Ubald8fu2sv0FERreij/kiss557RGRdc7yfiJyyu9vf2jNd2S8JyY6ilvGZvDXHYcpOFbhdjjN9vKqvcRGS0QMLTeXiPzY77eaJyK1IvKVQ0sR+V8R2e237khneYM5pZn9B9xeRP7Br691InJaRGYFZ69NpIuKEm4dl8kXu4+RfyhyJizM+2IfAtwyNsPtUJqltSM4PwGWqOogYInz+e84yeZhYDwwDnjYrxD6PXA3MMh5TW2i3Wv81r3H2f6c+lDVW1R1pKqOBN4C3vYLe2f931T12+f43RgPmz0uAwFe/SIyjsAqqmp4K7eQacNS6dk5zu1wgkZVf+H3O34QWKaqxxpY/cd+v+t1zrKAOaUFAm6vqp/6xXUlUAF80sK2jYf9Y3Y6sdHCvNUFbofSLJU1tcz/soBJQ1Po062j2+E0S2sLnJnAi877F4FARyhTgEWqekxVS4BFwFQRSQW6quoq9d0Q4CW/7RtqdybwkvqsAro57ZxLH4DvCAy4GXi1dV+FaU9SEzty5ZAUXs8ppKom/C82fnddMScqa7hjQl+3QwmlW2n577ihnFI/OvSlMzLzny3d3s9NwEeqGjnDfSbkenaOY0pWb97MLYiIh/h+tPEAx05WccdFkZNDWlvgpKhq/ZPDDgCBbqyRBviXqIXOsjTn/dnLG2u3sbZa2ke9y4CDqrrDb1l/EVkrIstE5LIA+2QMXxufyZHyShZtDu8HcKoqL63cy9DUrozpG/AscsQTkQR8o7NvNbLaE06x8oyI1A9jBcwdIjIZ36jMOGAkMEZELg/QZkO5x99s7ADKBHDb+EzKTtfwwYbwfwDnSyv3MKBnJy4Z2NPtUJqtyQJHRBY757bPfs30X88ZIQn6rRlD1a6fs4/69gOZqjoKuB+YJyJdA20oIveISI6I5Bw+fDiEIZpwdPngZNK6deTlVXvdDqVRa/aVsGV/GXdM6Ovli4tnACsaOT31IDAEGAskAQ800d5k57UWWONsO6ilQTmjOcOAhQ383XJIO3bRgB4M6NmJl1eHdw7JKyplzb7jfG1C34i4uLhekwWOqk5S1QsDvN4FDvoN56YCgebNFgH+VySlO8uKnPdnL6eRdhtrq6V9ICIxwA3AfL/9rVTVo877XGAnMLiB7+Z5Vc1W1ezk5ORAqxgPi44Sbp/Ql5W7jrKpOHzvSvrnlXvpEhfDzJHeeLCmiHzX7+Ld+p1qdJREVfc7p5EqgRfwjcxAw7lDgJ/5XbNznqrODdB3Q9vXuxl4R1WrG4jLckg7JiLccVFf1u47Tu7eErfDadArq/cSHxvFTaMja4JCa09RLQDqZyzdBbwbYJ2FwGQR6e5c+DsZWOicgioTkQnOdTB3+m3fULsLgDudmQsTgFKnnXPpA2ASsFVVz5zGEpFkEYl23g/Ad9S269y+HuN1t43PpFOHaOYs3+12KAEdKjvNhxsPcOOY9LC/62hzqeqzfoVHsYgkAlcQOP8AZw6U6q+5mwXkOX9qLKf8s4h0drZLE5FeZ/fdyPb1zuW6INOO3JydQdf4GOYsD89/ZkpOVvGXtcXMHJFGYkKs2+G0SGsLnCeBq0VkB75i4UkAEckWkTkAzpDxY8CXzutRv2Hke4E5QD6+kZKPGmsX+BBfsZEP/NHZ/lz7gMBHfZcDG8Q3bfxN4NuNDHubdi6xYyy3jM3kvfXFFB8/5XY4XzF3xW5q6ur4+iX93A4llK4HPlHVk/4LReRDvxGeV0RkI7AR6Ak87ixvKKd8AswDVjrbvQl0CdB3wO2d/vvhG91Z1uo9NJ7VKS6G2yf05eNNB9h79GTTG7SxF1fu4VR1Ld+4rL/bobSYRPITTcNJdna25uTkuB2GcUFhSQVX/OIzvnFpf346bajb4ZxReqqaS55cyj8M6cVvbh3ldjhniEiuqma7HUe4sRzSfh0qO82lT33K7HEZPDrzQrfDOaOiqoaLn1xKdt8k5twVPj/Z5uYQu5OxMa2U3j2B6cNSmbd6H2WnA15q4YqXV+2lvLKGb18xwO1QjDGN6NU1nlmj+vB6TgHHTla5Hc4Zr31RwPGKar4zcaDboZwTK3CMCYK7LxtAeWUNr4bJrddPV9fywordXDE4maw+9twpY8LdNy8bwOnqOv68MjxmVFXV1DFn+S7G9U+K2NtLWIFjTBAMS0/k4oE9+OPy3VRU1bgdDm/kFnKkvCpij7yMaW8Gp3ThyiG9eOHz3ZwIg5HgBeuLKS49HdE5xAocY4Lk/qsHc6S8khc/d/cIrLKmlueW7WRUZjfG9w//J/4aY3z+ddJgjldUuz4rs6a2jt99ls+Q3l2YODhyb19gBY4xQZLdL4krh/TiD8t2UnrKvSOwV1bto7DkFD+4apCXb+xnjOcMS09k2rDezFm+y9Vrcd7ILWTX4ZPcNymyc4gVOMYE0Q8nD6b0VLVr97QoO13Nb5bu4NLzenJFBB95GdNe3X/1YE5V1/L7z/Jd6b+iqoZnFm1nTN/uTMnq7UoMwWIFjjFBlNUnkWuHpzL3b7s5Ul7Z5v3/4bOdlFRU85NrhkT0kZcx7dV5vbpww+h0Xly5l/2lbX9vrbnLd3PoRCU/nRb5OcQKHGOC7P6rB1NZU8dvl7btEdj+0lPM/dtuZo3sw4VpNnPKmEj1g6sGoar8avGOplcOoiPllfxh2U6mZKUwpm/kX79nBY4xQTYguTOzx2bw51V72XqgrM36ffqT7ajCDyef32Z9GmOCLyMpgTsv6sf8nAI2FB5vs35/tXgHp2vq+LepQ9qsz1CyAseYEPjR5PPpEh/Dw+9uoi3uFp679xhv5BbyT5f0IyMpIeT9GWNC675Jg+jZOY7/eHcTdXWhzyEbC0t5ZfVebhuXycDkziHvry1YgWNMCHTv1IEHpg5h9e5jvJ5TENK+qmvr+OnbeaR168gPrhoU0r6MMW2jS3wsD00byvqC47y0ck9I+6qprePBdzbQo3McP5rinRFgK3CMCZFbsjOYMCCJx9/fEtKLBX+zZAfbDp7g0ZlZnnliuDEGZo7sw8Tzk3nq420UHKsIWT/P/XUXeUVlPDzjAhI7RtYTwxtjBY4xIRIVJTx143Bq6pQfvr6e2hAMM+fuLeG3n+Zz05h0rhqaEvT2jTHuERH+6/phREcJ/zp/HTW1dUHvI6+olGcWbWf68FSmD0sNevtusgLHmBDq26MTj1x3AZ/vPBr0+1qUnKzi+6+uJTWxIw/PuCCobRtjwkOfbh154voLydlbwjOLtwe17bLT1Xzv1bX06NyBJ2ZdGPHTws9mBY4xIXZzdgYzR/bh6UXb+XTroaC0WVunfP+1tRw+UcmzXxtNl3jvDCsbY/7ezJFpzB6bwbOf7uTjvP1BabOuTrl//noKjlXwm1tH0y2hQ1DaDSdW4BgTYiLCz24YxtDUrnzv1bVsLm7d1HFV5d//ksfyHUf4z5lZjMzoFqRIjTHh6pHrshid2Y375q9jfUHrpo6rKo99sJnFWw7y0PShjPPoM+uswDGmDSR0iGHuXWPpHBfDHXNXs/3giXNqR1V58uOtvPrFPu6dOJBbx2UGOVJjTDiKj43muTuySe4Sx51/+oK8otJzakdV+eXiHbywYg9fv6Qf/3Rxv+AGGkaswDGmjfROjGfe3eOJjhJmP7+KL/cca9H21bV1PPSXPJ5btovbJ2TyYw9N5zTGNC25SxzzvjmBznEx3PbHVXyef6RF29fWKY++v5lfLdnBTWPS+Y/pF3juuht/rSpwRCRJRBaJyA7nv90bWO8uZ50dInKX3/IxIrJRRPJF5NfifNMNtSs+v3bW3yAio5vRxxMiUiAi5WfFFCci8522VotIP7+/Pegs3yYiU1rzHRnjb0ByZ+Z/6yISO8bytT+uZs7yXc2aXVVYUsGtz69i3up9fGfiQB6b6b0LAs+FiAwRkZUiUikiP2pkvf7O7zzf+d13cJY3mAea2X9jeWS4E9smJ8/Fn9teGvP/MpISeO2eCaR0jefOP33B7z7Lb9bsqgOlp7l9zuozIzc/v3E4UVHeziGtHcH5CbBEVQcBS5zPf0dEkoCHgfHAOOBhv0Lo98DdwCDnNbWJdq/xW/ceZ/um+njPWXa2bwAlqnoe8AzwlNPWBcBsIMuJ53ciEt2ib8WYRvTv2Yl37r2Yywf35PEPtjDr2RUs2nwwYJI6WHaa/164jUlPL2PL/jJ+NXskD0yN/IfgBdEx4PvAfzex3lPAM87vvQTf7x8ayAMt0FAeiQFeBr6tqlnARKC6hW0bE1BGUgJv3Xsxk7NS+PnH25jx2xV8nHeA6gA55NCJ0/xy8Xau/J/PWFtQwi9uGs7DM7I8X9wASGtuIy8i24CJqrpfRFKBz1T1/LPWudVZ51vO5+eAz5zXp6o65Oz1Gmq3fltVfdW///rX2X3Ur+csK1fVzn6fFwKPqOpKJxkdAJJxiilV/dnZ6zX2XWRnZ2tOTk4Lv0HTnqkq764r5n8WbaPg2Cm6J8Qypm8SvRPjqKyuY/vBE+QVl1Fbp1w7PJUHpw0lrVtHt8NuNRHJVdXsILf5CFCuql8pdJyR4cNAb1WtEZGL8P2mpzSSB6KAJ/HlljjgWVV9LkDbDW1/DXCbqt7e3H2wHGLOxQcb9vOLhVvZc7SCxI6xZPftTmq3eKpq6th+sJy8olJq6pQpWSk8NO0CMntE/qNcmptDWnvb0xRVrZ+zdgAIdKexNMD/XvWFzrI05/3Zyxtrt7G2Ai1vzJltnKRXCvRwlq9qTlsicg++kSQyM+1iT9MyIsKsUWlMH57K4s0H+WTzQTYXl5Gz9xhxMVH07dGJ704cyI1j0unbo5Pb4UayHsBxVa1xPvv/phvKAzcApao6VkTigBUi8omq7j6r7Ya2HwyoUwAlA6+p6s/PDsxyiGmt6cNTmZKVwtKth1i46SCbikvJ3VdCXEwUmUkJfOuKAdw4Op0BHnm+VEs0WeCIyGKgd4A/PeT/QVVVRIJ+q9ZQtRsMqvo88Dz4jr5cDsdEqNjoKK4Zlso1HruLaISbDAwXkZucz4n4To2fXeA0JAa4FBgLVABLnKPOJf4rWQ4xwRATHcXkrN5Mzgr0T3X71WSBo6qTGvqbiBwUkVS/U0mB7mJWhG+Yt146vtNTRc57/+VFzvuG2i0CMgJs01Afjalvq9AZWk4EjjbShzEmDIjId/FduwcwTVWLm9jkKNBNRGKcURz/33RDeUCA76nqwrP6fgKYDqCqIxvZvhD4q6oecbb7EBiN75pCY0wbaO1FxguA+hlLdwHvBlhnITBZRLo7F/5OBhY6p6DKRGSCc478Tr/tG2p3AXCnM5tqAr4h5P0N9dGC2G8ClqrvgqQFwGxndkR/fEdtXzTr2zDGhJyqPquqI51XU8UNzu/6U3y/c/hqTgmUBxYC3xGRWAARGSwinVT1ofq+m7H9MBFJcAqfK4DNrdlvY0zLtPYanCeB10XkG8Be4GYAEcnGN3vgm6p6TEQeA750tnlUVetvAHIv8L9AR+Aj59Vgu8CHwDQgH9+w79cBGutDRH4O3AYkiEghMEdVHwHmAn8WkXx8MzFmO21tEpHX8SWjGuC7qlrbyu/JGBMCItIbyAG6AnUich9wgaqWOaMm33SKoAeA10TkcWAtvt8/NJAHgDlAP2CN30XKswKE0FAeKRGRp/HlJAU+VNUPgrv3xpjGtGoWlfl/NgPCmOYJxSwqL7AcYkzzNDeH2J2MjTHGGOM5VuAYY4wxxnOswDHGGGOM51iBY4wxxhjPsQLHGGOMMZ5js6iCREQO45vS3pSeQMuecR+Z2st+gu1rS/VV1eRgBOMllkO+or3sJ9i+tlSzcogVOG1MRHLawxTZ9rKfYPtq2lZ7+X/QXvYTbF9DxU5RGWOMMcZzrMAxxhhjjOdYgdP2nnc7gDbSXvYTbF9N22ov/w/ay36C7WtI2DU4xhhjjPEcG8ExxhhjjOdYgWOMMcYYz7ECp42IyFQR2SYi+SLyE7fjCRURyRCRT0Vks4hsEpEfuB1TKIlItIisFZH33Y4llESkm4i8KSJbRWSLiFzkdkztjeUQb7IcEsI+7Rqc0BORaGA7cDVQCHwJ3Kqqm10NLAREJBVIVdU1ItIFyAVmeXFfAUTkfiAb6Kqq17odT6iIyIvAclWdIyIdgARVPe52XO2F5RDLIZHOjRxiIzhtYxyQr6q7VLUKeA2Y6XJMIaGq+1V1jfP+BLAFSHM3qtAQkXRgOjDH7VhCSUQSgcuBuQCqWmXFTZuzHOJBlkNCywqctpEGFPh9LsSjP1h/ItIPGAWsdjeSkPkl8G9AnduBhFh/4DDwgjOUPkdEOrkdVDtjOcSbLIeEkBU4JiREpDPwFnCfqpa5HU+wici1wCFVzXU7ljYQA4wGfq+qo4CTgGevATHhwXKIp7iSQ6zAaRtFQIbf53RnmSeJSCy+xPSKqr7tdjwhcglwnYjswXe64EoRedndkEKmEChU1fqj6DfxJSvTdiyHeI/lkBCzAqdtfAkMEpH+zsVVs4EFLscUEiIi+M6zblHVp92OJ1RU9UFVTVfVfvj+fy5V1dtdDiskVPUAUCAi5zuLrgI8ecFnGLMc4jGWQ0KfQ2JC3YEBVa0RkX8BFgLRwJ9UdZPLYYXKJcAdwEYRWecs+6mqfuhiTKb1vge84vzjugv4usvxtCuWQyyHeECb5xCbJm6MMcYYz7FTVMYYY4zxHCtwjDHGGOM5VuAYY4wxxnOswDHGGGOM51iBY4wxxhjPsQLHGGOMMZ5jBY4xxhhjPOf/AN8TIvA1jyAIAAAAAElFTkSuQmCC",
"text/plain": [
"