{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Earth Density Models\n", "\n", "Test models of Earth density and electron fraction built into ASTERIA." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from astropy import units as u\n", "\n", "import numpy as np\n", "import matplotlib as mpl\n", "import matplotlib.pyplot as plt\n", "\n", "from asteria.matter import PREM, SimpleEarth" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Density vs. Radial Distance\n", "\n", "Compute radial distance and calculate the Earth's density using the PREM model and a layered approximation to PREM in which the density is treated as constant in each layer." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "r = np.linspace(0, 6370.99, 6372) * u.km\n", "prem = PREM()\n", "earth = SimpleEarth()\n", "\n", "rho = prem.density(r)\n", "rho_c = earth.density(r)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot Density Layers" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAAFgCAYAAACmDI9oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsnXmcFMX5/981M3uxLMsNy7XLKZewCKKgIIhHVMQjXmg8vkaMUWPwTDSHftFvjMErxiOan8YrgsYzJtFIUDxJ8AJBQRRBRRQQuWGX3Z3+/dHTszU1Vd09s7PsLtT79Wqmu+qpo3tm2c88+9RTwnEcLBaLxWKxWCwWi0ukqSdgsVgsFovFYrE0J6xAtlgsFovFYrFYJKxAtlgsFovFYrFYJKxAtlgsFovFYrFYJKxAtlgsFovFYrFYJKxAtlgsFovFYrFYJKxAtlgslr0UIcQDQoh1QoglmrorhBCOEKJjU8zNYrFYmhIrkC0Wi2Xv5UHge2qhEKIncDjwxe6ekMVisTQHrEC2WCyWvRTHcV4DvtNU3QZcBdidpCwWy15JrKknYLFYLJbmgxBiCvCV4ziLhBBBtucD5wMUFRWN7NmzZ0ZjxeNxIpGW4aexc809LWWeYOfaWDTFXJcvX/6t4zidguysQLZYLBYLAEKIVsAvgCPC2DuOcx9wH8CoUaOcd955J6Px5s2bx4QJEzKcZdNg55p7Wso8wc61sWiKuQohPg9jZwWyZY9HCFEShR6tYVAB7FsAg/Ogq4ACAXkCYsprxIE6B2odqAFq47DLcV+rauGralhcBR9uhY8c+NpxnJ1NfZ8WSw7oC/QGPO9xD+A9IcRox3G+adKZWSwWy27ECmRLiyUhfHuWwJBCGJbvCuD++dC1EEryoSAPIofhBttHAJE4jH1K537Bl07iiAN1wCFCOLsgvguqqmFzNXxdDZ8khPTirbDUcf9sbYW0pdniOM5ioLN3LYRYBYxyHOfbJpuUxWKxNAFWIFtaBEKIju3hmDZwagmMbgPtDoNIlPqVpkJ5bfCYIeqiQB5Q6BZFgWKg2IFuwEgHTpOF9Bgh6rbCui3w+mZ4bAu87DjO1hxN2WLJCCHELGAC0FEIsRq41nGc+5t2VhaLxdL0WIFsaXYoYviANtD2MIjkEewBbi6oYl0S0tH2UAacEodTakgRzW9sgVmb4d9WNFt2B47jTA2or9hNU7FYmj01NTWsXr2aqqqqpp6KL6WlpSxdurSppxGKxpxrYWEhPXr0IC8vL6v2ViBbmhQhRJv2cEJLFsOZIovmKCmi+eQ4nGwQzS84jtO8/1e2WCyWPZjVq1dTUlJCRUUFQRlempKtW7dSUlLS1NMIRWPN1XEcNmzYwOrVq+ndu3dWfViBbNntFAnRvyP8qhOceCgUF9A8xLBu/N2VBDZINO8Ehgnx3bfw4Hq4qcZx1u2mqVksFosFqKqqavbi2OIihKBDhw6sX78+6z6sQLY0OkIIUQKHdYJrO8EB4yAWo/F2qRGGc911Q1DFs+NTlw2yaG4NFEP7HnBZDVw2Soiqb2HOerhuu+O8l4PhLBaLxRKAFccth4a+V1YgWxoFIURRB/ifTnD5wdC7GESU3AhU3WK8pvgvK6z4VoVztuLZ87IXAF2hsAscWwvHjhGi7ltY/C383yZ41nGc2iyHsFgsFovFghXIlhySJ0TXznBNRzhzArQtouGhE6oYDtuXzi7XItokdNVyk5B3pNdsRLPAXfjXwQ3HqOwHf90BzhAhvtkA96yF3zuOsyWLri0Wi8Vi2atpGXsRWpotQohoFyGuGCHE1gnw9XD4SQ9oW0x93uHQfSUOL2dxFHP+YqGxj2jshaZt0Phh2+jsdfMx9Snbe7HHcptM8HLMlYCogLL9YMZE2LyvEGvbC3Faht1ZLBaLpRkSjUaprKxk6NChnHzyyezYsSOt/JRTTmHTpk0ArFq1iqKiIiorK5PHww8/DEBFRQXjxo1L6d/rQ8fXX3/N5MmTAdiwYQMTJ06kdevWXHzxxSl2jz/+OMOGDWPIkCFcddVVyfIvvviCiRMnMmLECIYNG8Y///nPwL4OO+wwNm7cmO3jahBWIFuyolCIHn2EeGks7BoJM7tB6wKyE8XZiGE/EeuVR3H/RJIH5CtHgeZQbcLY5SXGiCXG083JJKAzEc3ZPNNWQC/oPBpmHShETS8hHhJCtMmgK4vFYrE0I4qKili4cCFLliwhPz+fP/7xj2nl7dq146677kq26du3LwsXLkweZ511VrJu69atfPnllwCB6dZuvfVWpk2bBrgp1K6//npuvvnmFJsNGzZw5ZVXMnfuXD788EPWrl3L3LlzAbjhhhs45ZRTeP/995k9ezYXXnihb18AZ555JnfffXemjykn2BALS0Z0EOLUbnDHwdDZyz6RCXKoRFAYRNBGHTpBGhSK0ZiL9NRyXexxnVSum7duoZ9Q6jMJyRC4P+QdIdYezuoDZw0U4tM1cPYWx3krZDcWi8VikfnTn+Czz3LbZ58+kBCgYRg3bhwffPBBWvno0aNZvnx5qD5OOeUUHn/8ca644gpmzZrF1KlTeeSRR7S2Tz31FDfccAMAxcXFHHzwwXz66acpNp999hkDBgygU6dOgOsBfuqpp5g0aRJCCLZscaP+Nm/eTLdu3Xz7ApgyZQrjxo3jF7/4Raj7ySXWg2wJRAjRpqcQDx4gRM3+MLsXdC4ks3hgXdiBVxfGqyp7gj0vrue1lQ+/EAtTebaH3ziqB9jzZuu8z37ect2z0HncwxDBzYbRF/qNgTf3E2JHmRD/K4TID9mFxWKxWJoBtbW1vPDCC+y7774p5XV1dbz66qtMmTIlWbZixYqUEIvXX389WXfSSSfx9NNPA/D8889z7LHHasdbuXIl7dq1o6CgwHde/fr1Y9myZaxatYra2lqeffbZpIf6uuuu49FHH6VHjx4cffTR/OEPfwi8z3bt2lFdXc2GDRsCbXON9SBbjLQRYmw3ePAQ6N+KzL5NBYVA4FMne4ZNHmJdf0HzycTeD9l7Kwj25so28theWZT0BXuep1l9VrKdXJfJgj+BK87LoKgr/Hor/KqvEPO/hrN3OE76V3iLxWKxpJKBpzeX7Ny5k8rKSsD1IP/whz9MKV+1ahWVlZUcfvjhyTZeiIWO9u3b065dO2bPns2gQYNo1aqV1u7rr79OeoX9aNeuHffccw+nnnoqkUiEsWPH8lnC0z5r1izOOeccLr/8cubPn8+ZZ57J/PnzA/vs3Lkza9asoUOHDoG2ucQKZEsKQohIF7i6G/xiDBSpHk7ftmQuij1B7InviFJn6kd3HSasIlchFn45kMOUywSJ5jhQq9iqolgnpIPEsvdelYJoA2P7wCfDhNi4Bn76rePo/8ZmsVgslibDizU2lW/evJmjjjqKu+66i0suuSRUn6eeeioXXXQRDz74oO+4YbfYPvbYY5Oe6Pvuu49oNArA/fffz4svvgjAmDFjqKqqYsOGDZSWlvr2V1VVRVFRUaixc4kNsbAk6SzEhftD1X5wQxkUeds9B6GGSMjlulABz4PphRp4m4boPMe6kANd1opoyPJchldkMn6m85DbxEgNK9GFp5hCMcKGYIjEGD2h3Sh4uFKILe2FOCJEU4vFYrE0E0pLS/nd737HzTffTE1NTag2J5xwAldddRVHHnmk0WbAgAGsWrUqVH/r1rkbvW7cuJG7776b8847D4BevXolF+wtXbqUqqoqOnbs6NuX4zh88803VFRUhBo7l1gPsoUOQkzpBbNGQqtMBJXOW2zyIMvZGFSxjM952LIwc5VFeDap1MD1yMaVw/PUBnmwdQvwwu6+53mMvXuQQzHipD4b1bOcqVc5BnSHks7wr6FCrP0CjtriOO/7NLNYLBZLM2H48OEMHz6c2bNnM27cuGQMsse5556b4l0uKSnhZz/7mW+fxcXF9O3bl08//ZR+/foBUFFRwZYtW9i1axfPPvssL730EoMHD+anP/0pixYtAuDXv/41AwYMAOCWW25h2rRp3HbbbQghePDBB5O73Zn6evfddznwwAOJxXa/XLUCeS+mjRAHlMM/RkGHsKEUfnHBant5O+mIxsavj7Bi2AvRMAlwU7uGENWUqYJXFdK1hjn5iWbdhiOyEHcSc4knrusShyqW1WwYYYVyPlAOXcrgvYFCfPIFHL7DcT73aWaxhOfdd+n23HOwfDlEo+lHQQEUFaUerVq5R2kpFBaC3fbXspexbdu2UOXPP/988nznzp3aNjqPcEVFBUuWLNHaX3zxxTz44IPJTBYmj/KsWbO05YMHD+bNN99MKdu6datvX4888kgyHdzuxgrkvZDWQvTvCXMOhPIwYRR+3mKUci8kIKyXWA0X0NnJZV42CFPOYV2bTMrDECQsPVQhLXt9awkvmk2CWR5H9i6T6NvzLJu8ynKMs9/9FAB9oX83WNVXiP98CcfscpzvfJpZLMGsXk3p4sXw9ddQV1d/xP0+kRL5+a5Qlo/27aFzZ+jUyT06d3aFtcViaTAnnHDCbs8mMXToUCZNmrRbx/SwAnkvolCITj3hxQNgvzBp2jIJo/C8xWE8xTpRrBsjjBg2ifbGJMw8QC9ovfCOPMkmSDTrMlQEieU86r3Ycp+yB1kobfy8ygJ345GBcGAP+LZCiL99DlMdx9G7JiyWII47jnlVVQwbNiyluCA/nz7l5VBdzaqPP2bnpk1QXQ07d0J1NUXxOBVFRbB5M5+tWEH1xo2wciVs2QJbttBaCHomRPGn27dTU1hYL5w7d6akTx96VFZC9+4s/+or6urqUsYvLS1N5mddtmwZjuP+VHz++ecsXbqUdu3a0bVrVxzHYdmyZSltW7duTc+ePRvriVksTY4XT7y7mNZE2ULACuS9AiFEcTk8dSAcWUx2wtjPW+ydm+zCimJv4V4YMZytEM6lgA6TTk13rQudUEVzHKhJHKb4Yj+xLMcsRxLnteg9yGpZkEe5BMRgOK4HbO8hxD1fwU8dx6n1aWaxaPnHP/7Bo48+mlJWXl7OnXfeCbEYdz72GB9//HFK/cCBA5k5cyYAt158MZ9/8UV9ZVERIwYOZMZ558H69dx43XWsW70aVqyAqiqormZsu3ZcnYihvG7JErbGYvWhG61aMenoo5l++eUAXH311dTWuh/tDRs28OSTTzJ58mR+9KMfUVtbm7KNrsd9991HWVlZzp6RxWJpGqxA3oMRQojucPfB8KMSEEEpS0zC2OQtVu117YK8xHJmhqD45KC5727CjqkTxOq5auPFVRdQL253Ke3CimU5BMPzVptilYVk7yeUI0BbEMPgwh7woy5CXLHWcW73aWKxpHHkkUcyePDglLLCwsLk+bRp09ixY0dKvZyn9aKLLkpLPVVSUgL9+sHAgVz6+9+nruSvraW0uhpiMVi9mp+//TZ1a9bAmjVQUwM1NbR/6SX44gvo3ZtrDz0Up6ICysp45/33GTVqVHLVfSwWY8aMGcmuv/vuO5YuXZpMaWWxWFo2ViDvoZQIMWgYLOgOrcPk8jOlaJOvPSFr8hbrPMWqiIziiuIoeqGYKzEcNgyioegEqZ8gDmMjlHrvi4QXNlGHK5bjSls1TEIVuLJX2Vvcp25IApkL5Q4QbQu3DRLi0pWwX5Xj7P4tjywtku7duzNixAhj/T777OPbftCgQb71Q4cO9a0fdtJJ7kldnSuSV650j88+g3feoXLzZre+oID2BQWU5+fDwIHQpg2itDRt7k0VK2mxWHKPFch7GEII0QPuHwX/04pgURhGGOehF7w6MawTunmJQzeWapt2PyHrmsKDbBpXLdNlqggSzeqz8USrJ3DlGONd1McZq+Ja7ddDzs/s4IZy6ObjfSEKik+OAX2gVwdY31WIy79xnNsM5hZL8yMahZ493WP8eLfMcWD9eli2DJYtIzJnDjzzjCumwbUdNsw9hg6FNm2S8crCZtewNBLffPMN06dP5+2336agoICKigpuv/32ZCq1hvLss88yYMCAtL/s7I1YgbwHkYnXOEismoSxGoOMxsarKwwxjkqY8kx/9TT2r6pMYpF1IRUm0WwSzHJ5FCii3rNcpdQ5pHuW1bk5uB5qL+7ZJPrDLOQrBVEJtw4SYrr1JltaNEIkF/YxfjwrBgyg59ixbjzzhx/C4sUwdy784x8AvNuqFdd9/DE3/+Uv7DNkSBNP3rIn4jgOJ5xwAmeffTazZ88GYOHChaxduzanAnny5MlNJpBra2ubJOexDruT3h6AEEL0FOKBUfBRjwBx7HkETd5fTxjnY971TbdjnNdfPlCcONSd8XQ7vqlzUK91c9Ddk9+RKyZdfXXK9W8S+RuDxg+as3qPOhtdW/W5xIDWuM8+T7LV7fane0+98Jc8ZRyh6cvvi4znTR4F67sIcanB1GJpeeTnw6BBcNJJ8L//C7NmwcyZcOaZiLw8+PprnDVrmnqWlj2UV155hby8PC644IJkWWVlJePGjcNxHK688kqGDh3KgQceyOOPPw7AvHnzmDBhAieddBIDBw7kjDPOSP6l4+c//zmDBw9m2LBhXHHFFbz11lv87W9/48orr6SyspIVK1YY57Jq1SrGjRvHfvvtx3777cdbb72VHG/8+PGccMIJDB48mAsuuIB4InVj69atufzyy9lvv/2YNGkS69evB2DChAlcc801HHLIIfz+97/n888/Z9KkSQwbNoxJkybxRWIh7nHHHcfDDz8MwL333ssZZ5yR4yecSvOQ6ZasCes1NoVIqKEQOgGtttf1U0hqXHGmnuIge1O7MORKJE+65hpevvFGbd+ZZrTQtVGfgZplQtdO51n2vPcFuKEX1ZK99xkxeYM9r7MgNUWc7EFWr02e6VIQIxLe5BUwwuZOtuxxxGJuTPLAgYhdu2DOHJxam9Blb+FqxWkCMG7cOI4++miqq6u57rrr0uoPO+wwJk2axJYtW7hR+X2iXqssWbKEkSNHauuefvppFi5cyKJFi1i1ahUTJ05kfCJc6P333+fDDz+kW7duHHTQQbz55psMHjyYZ555hmXLliGEYNOmTbRt25YpU6YwefJkTvLi8w107tyZOXPmUFhYyCeffMLUqVN55513AFiwYAEfffQR5eXlfO973+Ppp5/mpJNOYvv27ey3337ccsstzJgxg//93/9N3vOmTZt49dVXATj22GM566yzOPvss3nggQe45JJLePbZZ7nvvvs46KCD6N27N7fccgv/+c9/fOfYUKwHuYWSiddY9fqpIjhG/cI51eur8xh77fKo9xbHFHs/b6gqusN4lv28wapN+/JyfrZ0Kaf86U9cuXgxpz/6KP0nTeKiN97gZ8uX02v//em1//5c/OabXPree1z85pt0HjAAAex/9tmc/dRTTHvhBX62fDmTb7oJARx9443kFRVx2fvvc7qSlgpg4hVXMH3BAi5ftIgjr7suq3mH8RSbnqWpP7/3SH1P1ffc+1zEpP51nyPTZ0/2Jo+Gb6032bInIyLuT4ITdqMTiyWHvPHGG0ydOpVoNErnzp055JBDePvttwEYPXo0PXr0IBKJUFlZyapVq2jTpg2FhYWcd955PP300ynZYcJQU1PDtGnT2HfffTn55JP56KOPknWjR4+mT58+RKNRpk6dyhtvvAFAJBLh1FNPBeAHP/hBshxIlgPMnz+f008/HYAzzzwzadelSxdmzJjBxIkTueWWW2jfvn0WTyo81oPcAikRYtBwWNAtQ6+xKqjkzT3CeowFrmcyZmijm4PuOqxnOZN6ue8O/frxyMkn8/T55/OTt99mv9NP556DD2bwlCkces01PH7WWfxx/HjidXX0mzSJo37zGx5JfGPuVlnJ70eMoLa6mis//pg3//AHXrz6ag66+GJuT6xal7fQ7n/44XTs3587Ro9GCME5f/sbfcaNY+Xrr4delKe7vyBPsdq3ei2UMi9WuQbXq6zzBKtzkkW0nB5O9Sj7LeSz3mTL3kAkkd4trmw8Ytlz8fP4FhQU+Na3adMm0GOsMmTIEJ588kltnRc2YZqLRzQaTcb5LliwgLlz5zJ79mzuvPNOXn755dBzue222+jSpQuLFi0iHo+npGdUF6maFq3K5cXFxcaxZLvFixfToUMH1uyGUCbrQW5h9BDiTyPho6CQCp3w9c69WFNdPmM/72IR9fGtft5Pnag2xdnq7HXlfmPovKAbV65k7ZIlOI7D2g8/5NO5cwH4ZvFi2ldUUFRayg/++lcuW7yYY2+7jS5DhiTbrpg7l+otW6irrmbdRx/Rvrzc+CUDYMARRzDgiCO49P33mf7ee3QeOJBO/fsb5xbm/sI8Fz/Psmyrjp+PG6tcqNgEeZTl2HT1WWTqTe4qxE8NphZLi6RT+/Yc37UrHUpLm3oqlj2UQw89lOrqav70pz8ly95++21effVVxo8fz+OPP05dXR3ffvstr732GqNHjzb2tW3bNjZv3szRRx/N7bffzsKFCwE3j/jWxPoaPzZv3kxZWRmRSIRHHnkkZUfKBQsWsHLlSuLxOI8//jgHH3wwAPF4PCnwH3vssWS5ytixY5OLEP/yl78k7RYsWMALL7zA+++/z80338zKlSsD59kQrAe5hSCEKBgAH+4Lff3S0Ad5jVWBI7+ahFwh6R7jIG9xJl7loDpdv35tBVBbXV0/z3icuupq9/7icSKxGEdefz2fvfIKj554Iu3Kyzl/3rykvdcWXG9QNBYzessjQEQI5t14I/+97760+eg8yKrHF4OdzoNsspc9uroxdF7lPNz3VfYoq2Opfyz2hLKc8SIbb/JwuL1CiCNXOc7RGjOLpcVR1qULP+zVCzp0aOqpWPZQhBA888wzTJ8+nd/+9rcUFhYm07yNHz+e+fPnM3z4cBzH4Xe/+x1du3ZN2w7dY+vWrRx33HFUVVXhOA633eZm5jzttNOYNm0ad9xxB08++SRz5swBSFkYCHDhhRfy/e9/n7/+9a9MnDgxxQM8ZswYfv7zn7N48eLkgj1wvcQffvghI0eOpLS0NLmQUOWOO+7g3HPPZebMmXTq1Ik///nPVFdXM23aNP785z/TrVs3brnlFs4991xefvll7r33Xu0cG4oVyC2AYiF6VsLHZVAU5DXWeSihPr5YF04hhwvI7Quo9xar/clj6K6D7LLpw9Qmm/aFpaVs+eorAEaec46PZT3xmhoisRhxZRHO8n/9i8Ovv56Ff/kLu7Zvp023btTV1LB9/frQYli181uYp7aV+5MFsNyHKrZlOy9zRTX1+ZDVz5kslOUvUbWJOnkcuX9dNKZIjDcIjhoixFcfQX/HcXZoTC2WFkMcqK6rI7+mBruXnqWx6NatG0888YS2bubMmcycOZOtW7e6O0riZoiYMGFC0ubOO+9Mni9YsCCtj4MOOiglnrhv377asfr3788HH3yQvJbDRVq1amUUv9dffz3XX3998nrr1q3MmzcvxaaiokIb7rFo0aLk+ZQpU5gyZQqQe2HsYUMsmjkdhDiqEj7v5iOO5T+j6649AaR6iP0W37VOtFP/nC+P6fenfJOdXKbzdAct8AuT9k0dSzeP1373O753441c8MYbydhB1Vad34L77mP6Bx9w2qOPppR/OmcOix57jAvnz2f6Bx9wxpNPUpD4zynTezA9B92c1DamPtT2sp08TiH+i/nUOXjeZDWlnOlzqBIFKqDbKNhYIkR/g5nF0iJYtmoVp7z7LouXL2/qqVgslhwg/AK7LU1LNyH+byBcU4RZZGTjNZYFDlJdHq7XOKzHOBNvsV9ZNh7goL5zYRtEJj85YWzVxXdh6sJsdW1qrzuXt6neSb0HWF7EF9e08VLCqfYoNjocYAs4y+C09Y6jd41Ymj2jRo1yvDRPYfFytLYEgua69O9/56of/pAZt9/OiKlTd9/ENLSU59pS5gnuXLt06RK4vXlzQPYgN3cae65Lly5Ne8+EEO86jjMqqK0NsWiGCCEiveHVYXCw3xtkWoQH5sVUOs9lDFcY6+yRyvCpV21MZUFtw3qEw9TlUgyHGQ/MAjaMrfoMdKESwlBuGstvlz6hOfeI4HqTZaEsjw2pgjdsbLIp5KINiOHweA8hDlntOBdpzCyWZo1I/CXKpnnbs3Ecx24l3kJoqAPYhlg0M2JClAyBdQN9xLHpz+dQn6HCLzuFXNYKNztF0J/8MdSrNn7tTN5uU3YLdUxdW9N9NdV/X2GeucnW1Jcp3MJUruvDb0y/+gj1nxFIHVf9QhUh9fOnjq/7DMjjFgJD4cL+QiwUQtgv75YWRTIPsk3ztsdSWFjIhg0bGiy8LI2P4zhs2LAhJf1cpthfQs2IUiGGj4R3OkLMJPBy5TWWwylM3lu1XBVSKmGFs6l9mH4zFb67Uyj7/Zepuwc1BEHn5VXbmDzIqrdW7UMdU7Xxq4f6bayrqN9dLxNvsjyOSNTpnlcM6A/DW8GGQiH6VjnOtxozi6XZYTcK2fPp0aMHq1evTm6R3FypqqpqkDDcnTTmXAsLC+nRo0fW7a1AbiZ0FuKHw+FPJWD8441JHHu7noWJNRa43sCoph+kMt2rWq+79hPT2YjiMAK3KURz0CYfQbaqvSp81Xo/sazLYhEkjE1t/OrB9fLW4YZdeJ8th/rwCd3nr0YqcyR7U1xyBOgBbVrBN+2EGL/Rcd7SmFkszYoO7dsztXt3yqJRcByf/8ktLZW8vDx69+7d1NMIZN68eYxIbGrV3GnOc7UCuRnQS4hHh8MZXtYIHSZxHCN9e2jZPhOvse5P7RjqTW3DCuNMPcumdmH7zTW6fhsSfyzbZSqWdV5lk6e5IULZIxtvsrwLnzwnv7jk9hAdDm92FeKybxznNo2ZxdJs6FBezuk9e8ITT8ALL8CQITB0KOy7L1RUQMRGNFosLQkrkJuYfkL8dwiMNuXNVEWuKnijmrpMvcYNFca6+lyJ4jB1zcVPY7q/oBzIap1cr/YRFEKhE7lhhG8YoZytN1mOWdaFXPiJ5NbAMLi1pxD7fOk4jZPsci9GCPEAMBlY5zjO0ETZTOBYYBewAvgfx3E2Nd0sWwa1bduy5eabab1iBfkffwxLlsB//uNWFhe7QnnkSPfo1KlpJ2uxWAKxArmJEEKIfrCwPwzzE8ec+P8fAAAgAElEQVSmxU5yXmPZNluvcUOFcVhvcbaiuCFieHcI6IbGH6t2fhku/EIovDI/r3K2QjlTb7JuTgJXdelEspoaziMfGAw/6iVEqy8c5yyNiSV7HgTuBB6WyuYAVzuOUyuEuAm4GvhZE8ytRfHFF1/w08su45JLLmHkmWcCULxzJwXLl1O7cCFb3n4bXnvNNe7eHYYPp/XYseRXVlID2u19S0pKyMvLY9euXWzbto28vLwWk77LYmnpWIHcBAghIgNgaV8YYPqjmymkQhXHqpdOPrL1GudCGGfqLc7Gu+zXvrHahIknDtNGfYZhQyq8Op0AVm3VOem8z2GEss42rDfZizOWhXA+6XHJ3rlp8V4eMBjOLBei+HPH+b7GxJIFjuO8JoSoUMpeki7/A5y0O+fUUsnPzwfcbXI9pk+fzqRJk/ika1euevllqKuDzZth5Ur45z+5um9fxnbrxuKyMq595x1o2xai9S6TGTNmMGLECN5++21++9vfAu6OZUOHDt29N2ex7IXYjUJ2M0KI6D6wog+UZyqOvYV4OnGs7oSXC6+xnxDORBjnWhRnU9+YXmTdT1DQT1XYelNcs84+yFa1UcfIpL3aj1y+g3qhK4tklGvTxiKmxXsk2nwML650nKMMJpYMSQjkv3shFkrd88DjjuM8amh7PnA+QJcuXUbOnj07o7G3bdtG69atM51ykxA0V8dxWLJkCTt37kyW9e7dm06dOrF161aWLl2aYi9qahgWi9F7zRpqFi1iydq1EI2ys6yMHeXl7OjVi/7Dh1NaWsp3333HwoUL+fe//82pp57K8OHDGzTX5kJLmSfYuTYWTTHXiRMnhtooxArk3YgQIm8wrCqHbpmKY91W0d4htykkdavghnqNsxHGfkI3G1Ectq4xRXC26ARpkJ2pLoxYNglY1SZboWxqK9c5QDX1XmJP8OpEsieU1b79RHIdsAxeW+k4hxhMLBlgEshCiF8Ao4ATnRC/KPb2nfQaRDwOy5bB/Pnw1luwbp27qK+yEiZMgAMP5OtNmzj//PO57LLLmDhxYtPNNYe0lHmCnWtj0RRztTvpNTOEELEh8GUv6JKpOM431KniuBXBIRVhvcYmz3Mm4Rdhxwiy9+sjExpDQAepBt2cGxp/7Nn4hVV4tqYwCl29qV81DEJtq7uHQtzPYhX1n1N1rvI26Jks3osCA2F8XyHeXOE4B2lMLA1ECHE27uK9SWHEsaWBRCIweLB7nHuuG4Lx5pswbx7ceivk5xMZMgQ2bSJeUxPYncViaThWIO8GhBDRQbAyU3EsCLcYT+CKY92OdNl6jXMhjBsSbqGbUxC5EMB+QjXTMcPEK6uiN2gO6ng6gSvbZiqUg8RvmNhk79zLsrIdveD17L1MLKbFeyaRPADG9hZizkrHOVxjYskSIcT3cBflHeI4zo6mns9ehxDQp497/OAHrmd53jzazZvHtUDFvfe6AvqII6B376aercWyx2IFciOTWJD3SQX0CCOO5awVYRbjxXC9dUKxQzrP1mucC2HsJ3D9yoLEbraCORMRHWQbFHusa69rY/IQq30EeYrlctM8MvEoZ+tNlssjuFkutifaqIv3kMrlxXvy+AGe5MMqhHh+leMcqzGxBCCEmAVMADoKIVYD1+JmrSgA5gh3s4v/ODbFXtMgBAwaBIMGkT9tGqPefx9eeQX+9S/4+9+hf39XKB9yCBQVBfdnsVhCYwVyIyKEEP3ho77QOxNxrO6MJ9fL9vmJIxNxHMZrrBOpYURvtsI4jCgOW5cLL3JYgu4tKJRCZ6c+C11Iha6NKq5NHl8UG93YDfUm68Yspj4VnFyn5kvOIyuRPLlciNmfO85pGhOLD47jTNUU37/bJ2IJpBZ4p66O8jPPpOzHP3bDL/71L7jrLrj/fhg3zhXL++zT1FO1WPYIrEBuRPrBe/1gnyBxrIpgVRx7hyyOi6h/87IRx9l4jcMKY5NI1bXLVBTnQgznQkRnGn6h8+xmE38cxlNs8gbrxsqlN1m2M8Ul1+Au4AN9XHI2IjkGDIJTewqx40vHOVdjYrG0eKqqqvi///s/pk2bxpQpU+DYY2HyZFi+3BXKr78Oc+ZA37607dkTxo6FROo5i8WSOVYgNxJ9hXhzAFQ2hjguxryBiFqGcm6y8/MaN0QYZ+ot9usjrLDdHV7koDEaK/7Y5ClW7RsihFGuw3qTw4jkfFyv7w7MojdbkTwY/qeHEFtWO850jYnF0qKJJvIjf/nllyxcuJBu3brRuXNnqisqWDp+PIweDW+/Da+/zsZXXuHb+fPpePzx7Jw4kY/Xr0/rr1evXrRv357t27fzySefABCLxRg4cCCxmJUGFov9KWgEeghx21AYa9ohTyeOo6Su6teJ4wjuYjy1LdJ5Y3iNg4SxWma6bgxR3BRe5GwX7+3u+GN1DqawCj/xm6k3WRXxqj24n/NiXJEMuRPJecBA+GknId5a7zhPaEwslhZLXl4e+fn5vPjii7z44oucffbZnHTSSWzcuJFf/epXKbabgJ8XFnLME0/wzUMP8auvvoKuXd0trxNceumlHHrooXz++ecp7S+55BIOP9yue7VYrEDOMR2FOKESppserJqRwivTieOIcp4rcZwLr3FYr7LOJky/JpqDFzkX3mOdnfoM/GKK5TYmr7Jf3LPaV5D4lev8BLVfyIU8jwiuSJYX7+VCJBcAA2F2iRDvb3WcTzQmFkuLJBaLceedd7Jx40YAOnfuDED79u256aabUmwXLFjAmClToLqasmef5abnnoOqKigrc+OUBw2iW/fuAJSXl3PTTTdRVVXFtddey/bt27FYLFYg55TWQvSthKcKDPVhF+Sp4jiKG3Oca3Gcjde4MYSxn+AMI3QbKprDtA/KWBHUn5/32GRj8gDr7HXlfkI5bNhFWG9ytnHJuRbJAmgDYiAsEUK0s2nKLHsSZWVllJWVpZTl5+czePDglLJ169bRvn17AAp//GMGn3OOG6f8zDPw0ENuCrmTToKxYykuLmbw4MFUVVUBYNNeWywuViDnCCFE0Uj4sBSETiBlK44zTePmJ5KDxHQmfZjq1bJM7fzKg+pM4+SCMGIf9IvwwtiqNtmEVHhlOruwYRVqv0HeZPm8OYnkzpA/GJYDPbBY9naKiuD44+GYY9zsF089Bb/7nRt2ceKJcNhhFBQU8MADD1AshWFYLHszpjVklgwZCB93gYIgcSy/ZiuOdfZy3+o4aqiGn8gO6gPFXr32KxMGO1P7oDr1Oehiu/36ashhIpM5ZXvPavq/IDs/G9N9BX22TO+pbjz1XDe+t/BUvT/5WeYZxtH9JyaAXtC9Qoh/aqotlr2TvDw4/HC45x645hpo0wbuvhvOPx/xwgt0atuWVq1aNfUsLZZmgRXIOaBCiL9VQM8w4tg71E1AMhHHKGUYzv1EVJA4ChJW8nkY8azb5c/Uvzo/tSwT4RlRjlwJZLXfTOfj936YnoWuLIwI9rMxfb50Y+uEqUkco5ybPgNyW3k3yExFsu65RYF94KiuQlylqbZY9l6EgDFj4Oab4YYboFMnnLvv5onDD+fDe++F2trgPiyWPRwrkBtIVyEuHQDH6jJW6H6Jm8SxfOwOcazOSVfnJ3x19bq5mGyChKRaFiSwwwjnxiaTufjdVybPxyvPpVDW1WXiTdbZ+N2DKpJ1z0YVyep4us8ZCdtBcFM7IQ7WVFssezdCwPDhcNNNMGMGj3z5JR/cey/86EduzLIVypa9GCuQG0BbIQ4cCLfqUrF7v8y9c++Xt0kce2W7WxwToo3azk/0hrVRy8KK4mzFsE4A5uoIO26Q9zgbsayWZRJWAf5jBb3HpvKGiuRi6dz0s5KvtNPN16MVMBDmFQnR2WBisezdCIHYbz8YPJj4iSdC27Zw553w4x/Da6+BXbhn2QuxAjlLCoXoOBDeMC1n0AmGfPzFcZTsxLGurzD2uvIgr3GuhbF6namAVMlGxDaUTMfM9l5zJZR1ny9Tvd8cciWSdZ9jnUhGKZN3kpTvW0UA7SE6AD4VQtiFyRaLgUgkQry83A29+PWvobAQZs6Eyy6DDz5o6ulZLLsV+8siC4QQsX1hRQeI6oSQ7s/R6oI8r06OudSlcgsjatX6sMI4F3WmetVGdx1mXD+yFb65EMy52ixEvVddlgq1rVCug8q8ctVGV59JnV8qOJTzoAwXujmr2S2gPmOFwP1C6QB1Ur/ez5Ius0V3KNkB7wHDsFgsaUQiEeLxOAgB++8PI0e6WS8efRR+8Qv3+pxzoKKiiWdqsTQ+1oOcBf1gQQ9oE1Yce7vkCcVOFsy6TUCaShz79adrp6v3E9A6zyBKmUnEmsb0s9N5KxtKtv0H2emegd/zCvvsdZ8JXX3QZ8E0z6DPT9jx1bJiqUz+mfHKopo+TM83AvSDfXsIca+m2mLZ64lEIql5kCMROPRQ+OMf4dxz4eOP4ZJL4Lbb4Ntvm26iFstuwHqQM6SrEFdWwgjTn3J1wiCqlKu/6IvR/5JX+1TLTWPqbBsqmiG9D796U5swZUF9ZmPjZxumfWNvFOLnWVa9yoJUz6yQztX2YTb60M1FaOrC5EX2a6O289tSWy7zPMnyFwXPQ+wt2tuljOnZqM81kdni/LZCPLzJcd7EYrEkeeihh3Ach++++y6trs2xxxI7/HCqZs1ix3PPwSuvwJQpMHky5OdTWlpKNBqlqqqKHTvS9+dp164dQrtLgMXSPLECOQMKheg4Cm4yZawIsyhP9Q7ubs9xULmuzq+d2jab/vza6sikPlf/HYeZp2mTkCBbtT6T7aTRlAUJZfU6rEhWBfbuEskR3Nj8KvSfG2/RniqS5f5lCoG+MFcI0dpxHLtM32JJ0Lp1a+bOncvtt9+eVveHP/yBiooK5nTtyn27dsGXX7phFzNmQM+e3P/003Tu0oXnn3+ehx9+OK39Y489RklJye64DYslJ1iBnAG94T3TTnk6AWlalOfVF5L6J2K1v+YgjsMKXZN9mD5MfYWpy1YMN0Q4h/Ec68Rq0NbQQXXy81Q9wbK9SShn6k1uiMc4lyLZwf2iWUf9bnrqfXue5BqlzhSP3AUKyuEZ4FgsFkuSgQMHctFFF6WVd+jQAYBhw4Zx0WWXuYUrV8ILL8DatZTcfDNcdBEjR47UCuGCggJuvfVWBg0axFFHHdWo92Cx5AIrkEOSCK3w3QxEFn9BuY7zSM3ninSeC3GshoBkKprD1umucyGMw4jlsIK6IWI4qG8I7zmWxWCmYtnPq+wnPjO9zqXHOFuRrGvvAAW4IjlOagiF16e3aK9W6UcnkhPxyJPbCnGQDbWwWOrp3r073bt3N9aXl5dTXl5eX3D++fDSS/DII3DJJfQ56ij6nHkmtG6d1va9996jqKioMaZtseQcu0gvBIVCdOxnCK0wZaxQPcNCso3g/rKX20B6P7qyTMWxLLh1i7xM/TREHJvuW21j6ttUFmYXuqbYLEQ3ftDGJpnuqOdXrn4GTe+z7lrXXjeurl1Dyk2fXUhvL9urG4mon/WoVOb3GQb3ZzARamEdBRZLtkSjcNRRcO+9bjzyiy/CBRe4McpK/mQhhJslw2JpAViBHILehtAKnRhQf2nLZbpf8joxgKEsSBx7okHX3iSSTELFJGL8RE+QUDPNJeh+goRmWDEsGuEIM15YwRz0vP3KdZ+nMLmQdddBdXJ50LyzEcmmnwv550ceX73PmKZ/03OUQi0sFktDKClxvcm33QZdusCtt8IvfwlffZU0SaaRs1haAFYgB9BViCvLDaEVul/oar5j9Rd42EV5fsLSJEiCxEaYfoL6C7I1iccwos4ry2QHPR3ZiFm1fSa2YcfyE8O6vnTlfvPUtQ/7PnrnYbeUNo0bZG/63Kvz8fv5iJAqktVDtx21V64ih1poqi0WS6b06eNuLnLhhbBiBVx8MfzlL7BrF0KI1DRyFkszxv5p0YdCITrun0FohSnu2KOIdHEUJI797HT1YW2DBK9arrv2E0Wm9n5lfnMyEUbMyjb5BVDSBtp3go6doV0n9zoadVN+RqIghPuXwXgdxONQWwtbNsGGde7x3bewbSvU1tT3G5S1QrVR701N5ea9OoZ2YTfb8OJw5b50fZjGEPjHEofJQhEmVlmtD8rKEcUNkahG/9nxdtqTU1TIz0FGCrWwWS0sllwQibhhFwceCPffD7Nnw6uvUhaNUlpa2tSzs1hCYQWyD73hvTZkFlqh894J3F/C6ta4fv2hKTOJ1zBhFbkUx0HzCmpv6i+o3NSPWicEdC6D/oOhoh+UtIVozC0XYTqANBXVtYdi5yQEZ0JAb/oWViyHT5fCxm+1XWhFZMq8pXLdYjm13E8468bzE8Z+wlu1U8V2Q0SyOj/dczEJaC+zRS36vMdqiI/fGDarhcXSCLRrB1dcAYcdBnffzY3r1sGWLbB9OxQXB7e3WJoQK5AN+GWt0AnDGKm/jGXBHMOcsUItM43RGOJYPfer87MLK4R181HLMxHFXlmnrrDP0IQYbgexmPKtJqmM1Z5E+tcf79pJ/lOPk+rPFTiICOTHoHN3h87dYczEetG8cb0rmj9eAls2hstcIT8bnVdZbuMnik2iVSdwVftcpHsLa6v27SeIdeMUYt6OWreJiFdus1pYLLuRykq480547DF4+ml4913QpJKzWJoTViBraGhohepJLpSuTV5YXZnOXj7PhTgOErgmMRxGGJv6yvRLgloWicCAwTDyIDdEIpLyIDwxnJih8ESwNJIngoXf6CpOqmh24kqdkxTQIuKQF00VzXW1sHYNLHgNvvgs2SptdFXsyqIwbDo3r6wh3mQ/z3Jji2T1HkwC2qMVqTvtqWLYhlpYLM2A/Hxu/e47uh96KKd+8gnMmEGPnj1h5Eh3gZ/F0sywAllD7waGVsj1RaT/4g4SwH51QWW5KA9jFySO/e5H7VtF17agEEYcCEMqoVWJrG0lAZwihr1rKelXijA2zFh+1x1ZmpJ67siviUMSyep1NM+hWy+H43/gxjVv3QQL/wMfLnS9zWFjir06WeSFCZNQ+9GJbT+RLI+VS5GslunaBIlkdSc9VfzaUAuLpXmwfPly6vr2hdtvhyeeoO0997iL+S68EMaMaerpWSwpWIGs0FmIc/fLYWhFVNNW7UMnPv3sTPKuscVxYwpj0/Ped6TrgS0okrRrxBPAkYQAVsSwiKSK4OS1OpIkmpPX0qUaZqGEWCRFMNR7lJ24RiTHU84jEYfS9g6HHA3jj4Jtm+GVF2DVJ+GFskm8Il3jc+31ESSS1bGC7HIhkrPxMudTv8teQ0Mt+sLkEiEGbXWcpVgslpyRzIOclwdnnMGnQtDxv/+F3/wGxo+HH/9Yu8GIxdIUWIEsIYQQo+CexgitSBlH0xeasjDiWBXUQbZ+5brzbMVxGBFv6qOwCCZNht4D3KwSgF4URyJKGT7iWBopxausm4UOSQyDIcQirr924vUiOZ54FW6ZwKGkrcOUqW5WjA/fhzfnuud+XlWdKJbL/MRtkL2fwM5GJAcJb1N7dR5B/eUq1KIQ6AkvuS8WiyVXqBuFVHXrBrfcAk8+6Wa6+PBDuPRSGD68CWdpsbhYgSzRBX7dEfLDCNqGhlaoBAlhXVmQyN5d4jjMPMMI44r+MPEoaO3tyuKFTHgC2CiORfo5grT447RYZN1M/FBEcsq1n0hOF8Yp9YmyWL7D8NEOw/aHDWthznOwfq2/Z1gnlMG/TZA32SSwG0Mkq2OglJvaNnaoRQ/o0VaIgzc5zhtYLJacEI1G0/Mgx2Jw2mluLPItt7ibixx3HJx1FuTnN81ELRasQE4ihIgdAL/SbSagE4oNCa3Q9Sv3n6kQ9ptnmPIgmyAx7Gfr92XAOx9SCYd8D2Le/4WetziS2Dg4EpFeJcFsEsqgvwa0oRZC965L6EIrkuWkeorlV1kERySRnOJRrkt7FTh07Oow9UdQtQP+8SR8tco/zEAWyrpruY3fdZAwzkYkm2zl+fuFU5gEcWOGWuQBveBvQHssFktO6NmzJx06dNBX9u8Pv/89PPQQPPccvPeemyKuT5/dO0mLJYEVyAl6wN3tIernDfVecx1a4VcXVgi3BHGs1lf0hSO/7y7AA1yxq3qLtUei97TQigzCLEBf5ovqQZbKAsMrIMV7nPQoRxNe5Eii74hbjiuYC1vFOfFM2L4FnnsMNqz3F4lymZ/YVa9VQaoK7GxFsm5uQcI5jBBWX+Uxsgm1UBFAN2jXQYiTNzjOXzUmFoslQ6666iqqqqp47LHHAPjoo49Ys2YNACNGjGDQoEFsnTqV59evhxdfhFNOgYMOgv33Z/8DDqB///589913vPjii2l9jx07loqKit15O5Y9HCuQASFEq4PhvCDvsffa2KEVurogIay28bMNM0aQuMVgG6Ztp65w3FQ3GwWQKowjUddK+yp7kRMCWSuUFXEsv3qzyCjVm+oxlj3KTv21NtxC8Rx7dRHVexx1XyOJVxFJCOc6hBOndanD6RfAd+vg2cdg+1azx1cVtyaxi3Kts0XTLoxINvUP6eOr7U1lpr7VdmFCLdSxdV7kKFAODwFWIFssOaKqqopZs2YBsGHDBhYtWgRAYWEhgwYNYtu2bcz6z3/cxXrffuvGJv/jH7S94Qb69+/Ppk2bku1lysrKrEC25BQrkIFyeLxEo5B0wjSPVMNchVaYxjMJTr/+/No3pThu0xaOPx3adkjo0jRh7O31nIn3WKSKZq8ONGJZeVKZ5kGGdBGcrNZkr0gRw5rwihRJ56S/xmvBSQhl6hDE6dDF4dzpsOZzeP5xqK6uv4Ow3mRdnZ9olc/DimS/dkh1urn69ZnLUIsapZ3uOXSBos5CTF/nOLdjsVgaTNu2bXn++ecBmDdvHhMmTEipLysrS9bjOPDyy/DHP8Lzz0PfvvTZf//6eoWrr76a4uJifvnLXzbmLVj2EvZ6gZwvRIcxMDmM91j2HKsH5Da0Qq2Tr8OIXFP7oHYme137TET1pGNg8H6eMBYkxbBOGKveY5NYBoNoRiqTnppfHuS0rNcSprAKSA+tUOs8MSw0HuWU8Ardq3BfZY9yQih3r3A4/0p4Yw4s/G94b7Ja3xgi2dSHLoQiSDg3VqiFem3yUkeACpgphLjDcRzVyWyxWBoTIWDSJBg4EG66CWbMgOOPh7PPdhf3KezatYuCgoImmKhlT2SvF8g94Z+6rIs6wakuzJNtGiu0ImxZGHHt166xxHG7DnDquW4e42RYRIoYjqH1IoNPeIXBeyykrzBpGSwU4ZwyYdM7JXuIMxDJKYv1JOEcr0sM479ID1C8x4kAgBShXEtEOIw/0qFyf5j1/6Cqqv4uwohd1VZXp2vnJ5JNYjWsSA4q04lYnQ2keoqF0taLRa5R2umeXUeIlcHvgCuwWCy7n+7d4eab4YEH4Nln3XRwV14JZWVppmlZMiyWLAlYvr9nUyxE33IYHUbAJiRaigiWj8YIrTC1M7UNsvUbM4w4FqTng/YTzocdA2demBDHkQhEY64gjsQgmgeRPOk85l57NtH8xHV+ankkr94+2S5RLven618uS/aVL7XNq28bzYOIVJe085tHnjT3WHqbZLncTyx1/BR7dd55ylhRiERo0x6mXQEjDgj/3qjvo+799vuchPl86fow9Wtq7/ezEeZLWoFiIx+6Z6D7vyDhRb5UCLHHuaaEEA8IIdYJIZZIZe2FEHOEEJ8kXts15RwtFsBN+XbBBXDNNbBmDfz0p/Daaykmwu8vgRZLhuzVArknvFSkKdf94lQ3D5F/KReRmQjIpC6TslyL4yA7U7t2HeBHV8KQkSAiIiHkVGGcpxGKssg0iMi09rKgDdMupohg+VqZSzSmt0sTyRpRa7ov9TmoolkV2yl16vw8mwiRqGD8kXDOxVBYqH9fgkQzSl3Q50D3qyis8NXNx68sTJ2pf1Uky/V5IcYAaAuRnu6CvT2NB4HvKWU/B+Y6jtMfmJu4tliaB2PGwB13QEUFzJwJd90Fu3Y19awseyB7bYhFGyFGHQB9TL/k5dcYqd5jeWGed65rF6bvsCIirBA2/YLPRhyrczD1L9eNOwwqx1C/CE8Om4jElGs1hEItk14DF+xpziH9GvBduKfFSfzd3WeBnvyaEmrhLcxzUs+9+GQ1J7J3LZTwCgTJTBck7t+pq39HHLdtm/Zxpl3hbjLy8WL/kAf1WrbDp52ujZ9tQ8MqdGXqXP1s8oDq+ieV9tmWn3BALPIpMSHOr3WcLewhOI7zmhCiQik+DpiQOH8ImAf8bLdNymIJonNnuPFGePRRdxe+Tz+Fn/+cAw880MYgW3KG2FvjdYYKsbYcOquySA2jAHdFvEkgFyt1aNrrXv1sclGmG9e7Pz9b3XkYQR0RcPL/QNcekBJrHImCSBwp14bMFSlliX7k+OMUwSyLYnnRnmmhnlemzN5XJHtxxk7qtVpmXKAnn0vCOWCzkDQBHa+VzutIxiTHazXXcRwHPnof5v49Rdan3ZljONc8gRQ7k73OFvzHybTMVKfayNc1QFXiPJ44HOna8z/FpVf1OcSBZTDvM8eZyB5EQiD/3XGcoYnrTY7jtJXqNzqOow2zEEKcD5wP0KVLl5GzZ8/OaOxt27bRurVuFUjzw8419+RiniUffUTPJ5/EEYLVJ5/M1oEDczS7VFrKMwU71yAmTpz4ruM4o4Ls9kqBXCLE/mNggW4TSzmWWOB6j2PStSyOo7gr5XViVfUqy32qNro6XVk24lhnk2txXFAAZ18MRcW4SjlFEKteY3lRnnSeUub1oclakeZZ9kS0cg4+AlkVxzphLGMQyWl5kX0Esk4cqx5mWeiqQlmuS3iLXUHsXddJ165IxnH4di3M/hPUSfkX/ERsGJEs+ckDRXLcUO9kWaaOHUZYe2J3O3qB7OAKZLVMl7JiG/AqtHEcZ6umukXSEIEsM2rUKOedd97JaGxdmvHCK+QAACAASURBVK/mip1r7snZPL/+Gn77W5wVK+CiixBHHdXwPhVayjMFO9cghBChBPJeGYPcDR7M05TrvKtyPLJKoaFcZ2/yJvuV6aSbSc6FFcdhbHX2JrsuZTDtck8cR5RYXV1srhSPrI2nNdSrC/FMsb3Cr70uttcUM+wXF63GExva6eKijQsGlfhsP3uhG0eNjXZjkzt2cRfwFRdn9j7rPncN+Uw1VlnQF0W5rMjHJibVmWzA/ULcHW7VVO1JrBVClAEkXtc18XwsFn/KymDmTK768kt+feONTT0byx7CXheDLIRofQgMDhKZUL+ltFcne4/VuGR8XtX+/dro2gWVNYU4FsCwkXDIUQmnrSeOPZGrepFlD3HSCyyfh90oJGL2GofaKMTgSQ6FxpOspnpTwypE4jy5lbQX5QrJdG/qIUTCVtTbO3HXS+zdlxd/7OVKlt8ZJ3EggFoKCuOcOx2eeRRWf17//snxvZAyM2Msblh7kUF9JmVh5ifPU7aTv+x6P9fy5iFR3A2+dfNAsusG5wDTfIZu6fwNOBv4beL1uaadjsUSgvx8yNO5viyW7NjrBHI3uLWVplznPRaGA+pXxquEEbImD5epLEhYB3nagrxifuOa+jzieNhnX9wsFcl4Y5NAluOOPbEclcRwVC+MjQv1hOFaJ5b9BLL61IJwEspJI5DlEApIDauQwyicOO5GKejFsRNPCFw5UMATyJH087hwhbZRIAPUEiHOiWfBG/+C9xfohag0Wsq5KlbD2qOxzUWZXKeKZ3U+aj9FwA5NX6AXyLr7aQ+xdkKcuNFxnjZ01WIQQszCXZDXUQixGrgWVxg/IYT4IfAFcHLTzdBisViahr1OIHeHc3RxJapMkr3HXr1Q6kyCVddfkIhVxwkjonXjqP0Gid1M7QUw+TToMwD08caxdEFs9Bor52rssV+8sXGxniKMA7eaDiuQg+KQE95inSjWxR/jJIRsRCOQvQwW3jhCaofZm+zdW5pABqhDOHEOPtIhrxAWvJZ7kawTxg31GOv68xvXTzR7XmQva4XOixwDagPmEQW6wV1AixfIjuNMNVRN2q0TsVgslmbGXiWQ2wlxwgH1qU+T+IlF9YDMvcc6Ie03dpg+M2mbS3F8wg+gZx/qxXEy7jdAIAd5jXUL9vy2mFa9ykC6Bxnp3BPFcoYL05Mz4SeKIySFsPDqpLAKrUA2pHxL2giSCcjSBLLsTVbfsYjrWVbfRacWQZwDDnGIxeCtlxtHJPvZNjSsQiWMoJavwfUib9f0Bal/STJ5kQXQFboWCtGtynHWGLqyWCy7G7tRiCWH7FUCuRvcpW74AemC04svlutlcWwSvaZXdZwgG10/YcpMc8iVOP7+WdC9AlccywviVIGc3E7aW3SXqENovMa5EMiJWWcSg4xInGYokD0BjANESctk4RtiIXmLtQK5LtUmmQtZJ5Dlc+W+nLrEB1h5J+NApA5BHSMPcohG4PV/B4tkP7Ern5vaSTNI1mdTZhLNOgHrF3ohe4pVL7LwaS9TCHSF+4HcL5e3WCxZMaFnT2JFuu2/LJbM2WsEcr4QZeOgLIwg9QvB0O28Jb/q2uVKMJva6voIahMklNXy4083iOOkEFbEcKAHWbNgL2URns9iPUi1xWujEcgpr8pTzsiLrAmxkMVyUggnxjIt0ksRt3WJMk8YC+W6LiF+vf7qUvtO9qEKZKHxKnv9eFd1VI5x2FUD/301czGs4idcTWI62zLdmGG9yB4FuAJZdx951OdF9sp0dt3hCCFE1HG8N8ZisTQlx/TpA610q4wslszZawRyN7hfl5ZNJ3Tlw4tZhNR8yGofQYI5rBjWnfu1awyhrPZxzMnQqx/14jjNe6ymZ5PDKDQhFzqvsdzGKJA98StSxbIco5wikEVqebJMuctMNgoRktwybhSihFwkxbJDMu44GYMsvcYT/ae8ahbmyQI5LoC6+oV/YBDICaTgW0Edo8c71O6Cd+f7i1TpSaQ8PZ3XWUdjCGJTma5cvRf1Z9bv58Qkvksh0gl+Atyuv2uLxbI7qaqthdpaYwpWiyUT9gqBLISIjoEjg8QiuN4jkz9RjT0O8jv6Ceewgtkk4IP6MPgPA/tX+5jwPeg7iFTPcZr3WMk9LHuLTYI5zWusE8qKRxnS7eWy5LksjJHOTYJY+L+ZSU8xpAhmkbh2nHrhrC7QE6ANj/BCJ5KiWRLbspzzsl6YvMmQaC9fx+ttU0iUeeHNEQdBnLGHOWzbCsuWaK0zikfOdahFY3mRBW6YxE7qfw48GzkEI6iP7jADK5AtlmbBr998k/yCAm5o6olY9gj2CoHcES4q1URO+InJoMPUDz71av8mG1M/QWVBnjCTrandgEGw7/7UC1R5Ew5hEMfGBXq6a40YTstqoXiMg7zJpiwWRi+y6amq1Yo8MmWvEN6rvEBPKF5jKc5YFsp442iEcZo3Wbof9TqeEMtpccgKcUDUIkQdRxwP676G7zakC1JpVinnmYpkk12QZzpTL3LYV3XXTHnOYRfrdYKS1kLsu81xFmOxWJoWu0jPkkP2CoHcHa73iyv28FucV6ix1/WhtlWFqJ+dqV9d+4YIZT9x7J23bQdHft/TlRHFG6wLq5DrVC+zLJil6xShHCVVIKu5jzUCOVkvnYMinpWnJ6R3WBduoUUTViELYTke2XFwNwCRvMdy/DEiUa4syEsKYunViacL44iAeMK36YVVqNfelOVcZqb7SVwL4px+Ptx3M+yq8RfJJvHqZyvX+wlbdXZ+glhno7s7Uxtw/2JUQ6oAFpCSDs5PzOcBZfAgMFIzrMVi2c2YwrwslkzZ4wVyayGGjoU2OgGqCk+dOPbqoopt0KuJMIJZVxdG9OrsMhHSSdGQB2dckNBl3g55QvYaG8RxSgiGvIOeYUc937AKH2Esi98UL7KoL097JbUMyNjboIZayHHHwitLqNIU77HnWRakxx+rQlleoCct+AsMs1Cu4yTENKnCOfWGpKOWaF6csy6C+28P5yWWCbINaqMSVhD7zSlM23xcgawjhrtYTxXgav/dYIQQoshxnJ2GriwWy24gw//RLRZf9niB3AXu120+qRPHCamVJo7zlGu1D13fYYVzkG02oljXNqhvuf7MCyGWR704TmankI8Q4ljrPZYEs3wue41TvMyC+oV7ariFkK6lcyAlu0XyOnGH6gYhQUI56TVOFpAed0yqt1gNtVC9xwhSFuR54lleoEfcfW6yNxlh8CLXpV6DaxOJJV51CQ5VailuE+eEM+CZv/iHMoC/SNXZEmCXKy+yTpCb+lU9xerPhNqfZytTDKLMjUW+EovF0nTYEAtLDtkbBPLIMCLT70HkG9qofZnGkK+DBHM27cIKZVM/ctlxU6GkFEkcBx0+4lhN96b1HhsW6iWFshdfrAjjUF5k6utAI5INT1SIetGbYiIJYlkci0Rd0rvreY+lnfSE0HuPvXAKTzCD4jHWeJOT6d4kjzK4554o9vDEcTLcwvMYyzgpR48+DvsfDAveUB6LZO0nVlVbAuxy7UX261NHIe7207Ig9vqLYfYwy3PpBOdhBbLF0qQcUVGByNO5xCyWzNmjBXKxEKMOro+OSCIvwPFeTYfsWZbRiVFVcgWJYZ1g1QlXU1s/+2yE9PD95XRuft5ixQOcqTg2eY+NIRaCzLzIqa8Pvvg+R+zfn24d25DqYZDO1XJtOjdwxa5UnhTHkrdY5z12BIHeY1UY67zJsudY9ijj9S29enjuUXA/zN4cI1FSBHOiXETqOHCiwxefwTdrzEI1SPDKghVSxaufsPUT10Hjm/rze5X/P1B/NuSfGUeqV+fWHtraMAuLpWmZVF4O+fnBhhZLCLRLePYUOsF1Yb4BePmNPeRfkpnmU/QTxUGCWVenE9FBfZp+wfvVFxXBuCNBeFtIq+ncjII4hDhOSwMnL+aTPNXJxYDquNJr2pw081B2+HvwhfdYs2GbuQ95Tol519aRfp86b7m60DBtvhHqHJH+HJJfHDRjCGUsda4pzzmS2pcui4gc7oLy/kr3LJ8L4e6c6Glv02cqm8+f7jMc9GVT14fuy2WYL6W6sQuUa48I5vUHMkVAR/ihYViLxbIb2FxdzeaqqqaehmUPYY/2IHeEw8KISmE4QO9tDvNqGlO+9mubiSjWCQxTvWmM085zoypSU7rphFSGnmO/OGS/MAsfL/KtTy7ggRcWAnDeMSOZfvJYVn2zmclXP8KSh6aDENw863W27dzF0D5deefjrzhjxuMUFeQx/0+X8NHKtVx2x3Ns27mLjm2LefCXUynrWMqEC+9k7L69efODlUwZN4TLT5+YfFbbdlTzk1uf5p2lXyCE4Npzj+D7E/Zl1r/f4zcPv4zjOBwzZiA3/fh7IKD1Eddy2SkH8a8Fn3DLj4+kqCDGZXe96I5ZWsSDV06mrENxqvfYe0dSsllIad68OmMWi1hq5gpt7LFT/xp3pPfCUUIwHGL5cY45Gf7+hDmeN1uPsikUwq+/TLzIKkFe5Dyg2tA2Csg+edPPWke4DLgzxHQsFksj8Nv//hcRjfKbpp6IZY9gjxXIeUJ0nZC+t0eguJSvCzT1aK51fYcVzn79+nnPsikzCenKA6BNO0hblJciiHXeXVMsskkce7mN5dAK1QsqJFGshFVEory7/Gv+/OIi/nvPj3CAA358H4eM6Ee7Nq2SNiDwtp4+adII7nx6Pjf/ZAqjBvWipraOn9z6DM/NnEandiU8Puc9fnHvCzzwyzNACDZtr+LVe6fXP5xEeMX1f55DaesiFj92NTgOG7dsZ82Gbfzsnn/y7gOX0a6kkCOm38uzbyzj+HFD2J4Q5zPOO4KamloO+cl9PPebM+hU2orH5y7iF39+jQeunJx4RxLqVwgl7EIoC/OgfuGeGn9clxDHsfo+VMHsKyXlOOT66977OHTpDt98ZRammQpjnZ0sunV1DQ2rMM3Nu4b6zUEiAWOY+ugIFUII4TiO34O2WCwWSwtgjxXIneEqv+wVHlFNmYf3cDL1FvuNaWrbmN5jvzGKiuDgw6kXoqY/5ZsW56WFXviJYyUcIS3FmyBFkGu8yG8s+YoTxg2muFUhiAgnHjKE1xd/zpSDhyRuTNOeeq/4x1+uY8ln33D4JXcDUBePU9ahNCmsTz18FCm5kgFw+Pfby5l9wzmJPqFd2za89uoiJuzXn07t2wAOZxw5itcWruL4Q4YTjUb4/oThICJ8/OValqxcx+GXPwQ41NXFKetQ4o4pL8JTF+2B4jGWhHHSo5x4TYkzTnxy47XufePFHMdcwZ+MPdaRiJkWUSCOEHV8/0y456aEwxm98IVwgtRPuPrZq3WZCusw5FO/e57c1guzqAvouxhECRwOvJTBsBaLxWJphuyxArkDnBNGOOoW6HiHaXFeUJ9BAjpbwSxfZ+M91pWlhVYgLZ5LW6gnxd2mhEUo4RJBu+mpIlnNYGESyAgczzvseYpx62J5ecQdkuVVuxJyJhJLPDD3C4BDhCF9ypj/56s0TwuKi4rqRaX0sBxARBMe8kSJkxSxibJkKIigMD9GNC8fnDiOiDCkdxfm//FCt6fkwr049d5jaTEeifnGAeo0wliQkm9M1vPxxLxTMlkkBDIJcex5mL2Fel56upT3M54U0mFCLXTkwovs15eJIG+y6VUNp5KJkBpmoSOKu+4BK5AtFoulxbNHLtITQhS2h3Zhxa3OG5vpOlg/oRsknP3EsMlO7TtoLjq7EX6hFXiiVBbAUtgFkphSF5+liGVNBgtVHKcsgFPGkQV4JMb4yj48+/pSdlTXsb26jmde+5Bxw/vSpUMp6zZuY8OWKqpr4/z9zY+Sor+kVSFbd9ZAJMo+FWWs37SN+YtXgYhSU/f/2Xv3eDmqMt/7u3rv7NzvgUi4C5IASQgYUAElAuqMiONx8C4D4hmOHh3Hl5lz5owzjjPveR3njJ7zUfHo6CB4FxQcQQcdR5mIF0QC4X5VuSYBQkhCEnLZu3u9f1St6qeeWququnfv7OzO+n0+napa96runf7Wr5+1ynLPw0+1wzsaDXyu+atfejyf+faN2Rg3b9vFS5YdxU/X/pZntu6k2YJv/vttnHHSi9pwmp7b4sNfwMYtO7jp3sfbfT76jApNafiP9bXXD1+RTnzuepvieeiQmdykvcD72RiAhuHIxbDw4PBnyfc5r8qv81kvK1NWtpObVDm+AXUsX7pdXxsHwMk1hhwVFRUVtY+rLx3k+XDhNE+6/pLTq1dI1Q2v6ORLX+eXfXl3AtyhtsrKDQ7CaWeDP7RCwJfP0Q0t/5bFD+vwC88EvbLQipL1j09achgXnrOSU/44mQv1n1//Uk489nDA8Dfv/j1e8u7/zZGLFrDkiBdk9S489zTe87FvMnXKJG664i+5+h/fwwc+fiVbt+9kpNnig28/m+OPPlT0JSe2JV7lX//n1/G+f/g6S9/ytww0Gnzk4nN545kn8bH3/yde+d5PYq3ltact5Q9WnZjVSdxZy9DkBlf//bv4wP/5Dlt37GRkpMUH33w6xx+5kATJnBWstz73GPLLvRFwkt0b7kIq2ueSlw3kFUMt3vAO+Pw/hkMcytzhqrROyvvc6E5d5VDbQ8BO2n9nsv0yF9xpFgxONeaYndY+WDKcqKioMdA5Rx+d/iQaFTV69SUgHwB/XgaMTvrPqI4L5lM3DlZV/aq+Q3BcF55f8wcuGiGF08rQCl8MsgZoTyhGJRxr+BYT+QLrH1/ytjO55O1nibTknfvAW8/iA287Oz3pdt4fnr2SPzz75KzciiVHcuPlHypck9WX/Y/glZ8xYzpf/v8uTpJsGyrf/vsv4+2//1KRnoDl9hs/lR63AMuKxYdx4z99AG+IhXuISO7BICYfcgEKhgOQLNc3dqEUwRhk9aoItZg8pcWyF8OdtxavUN1QClR+HdjtpkwVzIbqD3jKuvxB8g8N8bU9CCyAvwHeGWgqKipqjHT6IYdEQI7qmfoOkI0x5nR4YSdAqX9GDT1auuOxePbLxlWWVifPV843hsmT4ahjof0zvIDWYGiFjgn2/FxfCJ/wAbSnzdKJeh5Yzpzesn3axyDS3JUwKk2XSaUXJLAp0Dr8Mq6MfLl2HPzawL57WZII1wGxhdwEPTduPWnP7bsHm2RAPNCG5WxljNQRdiAs32vc9XWQPAANBc+0eMWr4a5bw8DaSxe519Ac6kce662Get2uVgM4AN7gP7OoqKix1NM7dkCjwYHjPZCovlDfAfJMeOV0z3eX/vJzWKDzIQHksrqdbn19dFPG98VdNs5Q+69/W8pFxpBBUm5VikBohfeloDeDVA3CdeC4qq86kIw7OQXKASgOXimXrDBKhl9kYOweFZ2m5WBUw3GzDalZv6EQC7WkW+YeDxSdZHcjkJuE524ObDsuOhl4uvybTSbwZSBM/lwzMHbvVYOBSS1e9kr45X8Ur16nsBwCWl9eJ2WqVFVviGRNZEP+XPRx6BznwnRjzBxr7ZYuhhcVFdWlPnnrrVjgY+M9kKi+UN/9FnEA/F3oZ1IpXxkJz/K4TJ2U6QaYq/rqFJRnzoKDDiVxImXoAwJEK0MrPGEU3nAKGTYhnWcJx3qiWeCVc5gVWGfxzwOiL6PyNUwPtMvrkBFvmIma5OZguzCZTt8gGHGehvyNQCDu2he64iuj36/C+yDTTLFu6bUWn4nczZPhxafCQCP8Gav6nJb9HZR9lqv+1rq9YdXbkGtgKP6/4RvTZGAh/GnpYKOiosZGcRnyqB6p7wB5AbzE90UYcmGNyh/A/2Vc98u5Tvk6HmYdmPblV7XxB29POdELPxKMNMx5oE7CZnAVCwG1uTwfEHpc56D77ANjB8AhUJbnaYrnEHxJIBaAHgRjHwS7cSnYzy2dN1C8Ll5I9gBzAdQlAKt0d43cuRdCLXR9eRNlaAzAWa+r9/mt+1kta6uqfp2/m7q/G/j+X5Bl9M1z6Pzmw8WBLqOiosZQEY+jeqW+AmRjzLQZxQiJdr7Y6i9Vd1y1vFuVK1VVr6ytOvVGA8oLFsLcA0jdY+UcFyCtKuTB4zhqGPVBtdf9le5waNvwlBX7GRjrYwHKBfdYucvBlwZij4ucA2PjuYa+/QE1Js/WB8nBukaMU4Cwz1mW59NQ10ODN/IaJNd38XIYmlTvM1n22e0GdqvSRitDfhUbHyiXjcEAs2HhGAwtKiqqRGPx/0HU/qu+AuQZcLoPcPUfzYAnTebJOnVBuM4XftmYqvro1JHzfWmfcx4e91jCngBmn4voC62QjqMOZwiCoW9fwpwPmiVca/CTbrI6lqCsYTjkHudAuMJFDoGxKyvDK3KhFtoxDt0clECyXvM4eCMiX6Y8X7vI+to5F7kBZ7++/PNa58uq08/8aP8O625DN8oNwv9HSE2DAWPMnIrhRkVFRUXto+qrSXpz4G0+4tdfZPJn0tDLV7+q3bLynQBzN+5aVVuTp8CseYTdY1+ohYNCDO2wCwlQrk7ApZXQmi3fVgXHGuY8TrJ0bn37UDzO9l1+BWoV8k37Ylo5oc09+W6gvXyb8RwzkN9Pl3/Lv7t6JQu5ooXrKp1w1xpJ9wfyS7y5rWnQXt5tgOyHRzmBr7CsWyNfVi79Jj8rtMBajlqSPEtQT1qD8OQ7X5rO95UPSfffTZnQ+NxWfxIaVD9VbxIwO3ns9LcrikZFRfVIbzjmmPEeQlQfqa8AeRacVbesDyI7fXpeN33ptNG6bLJMWbkzz3HM53EvMwD2ucbyZ3dfaIJyiwuQLUFZQ7MHfkvh2LnClO+DAuMGfiCuAc0Sht2qFbptHxg7eM6Om55996hpHxTLbXpL5wCYZhiMse1rq5d3c2CcA2EBwY2B9r7xfCZwazenLvKAZekKuOt2z6WlN2AcKtMrwPbJvbuDwIhKk/v6pkCqAcyGtxEBOSpqr+nkRYviJL2onqmvAHkmHBRyc0OOkNwf7dPz6kJvGTCXpdepJ89Pbo9aAu2f4j1xpTpswv3UruE5A0oJtCpUIJfnK6PdYu0gh45l2506yCINAsDsu8AShtPdHDS32uPTy75hwg5y5syWOcfCQQYByPLY5sFYp+dgWVxzK8Yo03LvV3oucqymAY0G2MTSPvUsuPv26uXayuC3DtCOBnp9LnfZ1skBsvHk1dFMOK3LIUdFRXWhJ557Dttqceh4DySqL9Q3gGyMmXam53x8P4/WgdGyNJkeAnJfmU4guE67oXyt41akBqBe91iGUHjdY/kY6YDDXBZa4dKkUxtaBSNzkEvgOKsHtdxkd2WCDrLL17dNWsI9BnLQ7Opl8CtxqsRBppEH5ipItgKAsen7KAEZPyhnoRY1wim0o+xzkXPrIhsmT7PMmQebn+0cjMvSOqlfF3jrytUbKMmv6ssAs2B+F91HRUV1qc+tXctIs8n/Gu+BRPWFGtVFJoZCE/S09AlrNKrpK9ZWpzDbSVkfTPvA+jQXeOLcvxycNvzusZyM5nOPc2BdFVohXWPdpi+MIhRm4cat2zX5dHeuOVgnnx6auFd3FQvpRPsm57l+smO1L938snCSYFl5/RoEr3PQtZc3R+l7XvhlQHxecjdXrn4DY5LJep3e0JV9XqvKjKbdTtoqK6vh2Vc2nag3O9BMVFRUVNQ+rL5xkOtO0PN4f0D5yhZV7fi2Ve10UqebMm5/9jyYMo1kcl4KNV4ozq1uICCyECbhe0iIcGE1fHtBWcKbDwoD0IgnD/L7XsfYl0a77pQ5MGUuDAwJcLRgmzCyG3Zugj3b2uEUrp4LtbDphDsZapG5w6mDq93kXJgFlMciC0cY2mk5xzgQauFCKfSNiZ64p8Mp3GcF9bmw6jNkW7zgEEujAc1WduWyz2Dd2OCyPF3GV7ZO/U5U9Xdd5+8xnaj3amIcclTUXpONMchRPVLfAPIsOLNuWR9I6gvRCbyW9dFJ/TquVydOF8jJeUX3r+gySmcxbVlDcCgGOVdOgXcutEI7oDUdVG94hQJfHV6RG0OqSdNg5sEwNBMGBov5Ps04MIVeC81h2L0VnnscWsPtfh0UZ5P4BBi761UIqWjUgOQ0zwoQppEeO1c3FGoh3t+WiEk2DQXOvvfbAXCzCMsZKCd1Gw3LytPg5p/5L19dMA4BdTdArOtWbUMaIE7Ui4qaKOr2+zoqyqe+AeSZsKjK5dGusdzXE/TK1I2r1Emdbsfgc70WHUYCQDkwlj+dOyBWP7/7HhgRDKfQwFzyk790MwsA3gEcy7SCQ+xxjKe/AGYdUgy5qLzyIu7YGGhMhkkHwvQDoLkHtjwCuza3x2tT51neSEg3uSNI9sHxgIDj9Id+Hxhbd/08LrJzj61wjGWavEkKwbJMo8WylfDrn3XnGIfAOFTel14XeEMK1Y8T9aKiJpiigxzVI/UFIFdN0HNb6Y2GylaldTSuLvN0mU7KSs2dDwMD5AFWxiDnJmApEM4ASbjAGscL8FwGyoH9XFueLaKcF5glKKtjMwBzj4aps/Nlsj7FOZW5yKGl3gYnw/xjkrTnNyawnMGvg830OmducgUkW9uGade3hGOa4lyaKQg32mOSYOzCLUwjuamRy77J98m3ooVc1k3fzOQc5AY0DNNm+MMsxNWuvV5xHdUpGwLvTurGiXpRURNHbz7uOOzw8HgPI6pPpKlhQmpfnaCnVQbh3TrLZWM+5RWiVC7sIeQiK0AKLQlXB4q1O4wDNw3PA8X+C1sCcDzQ3pftA8x4ASxaCdPmesrI83CgXCI9WS+bqJfmNRowYyEsOhkmzxbjFGXkOTc8E/Yk8GdbecMwUGwneH3V9cvdQBhV1/cSn0r5edDvbaP9vhkDxxzvuXTlV7bjMqO96ex0HGU31P0+Uc8Y8/8YY+4xxtxtjPmmMWbKeI8pKqpMyxcu5IQDDxzvYUT1ifoCkCfKBL06qgPKvjxfvSOPoQ1wBeByxz4IEsCq4drrHFeBsoM7D7hp1xjyRuksWQAAIABJREFU28o4ZMTY07ILT4DZhyvQk5MF1ZVzMBt8uXN3VdI23KoWLq/RgAVLYMFx4tqKfAnNXkgWW9/5Z9dftOt15QMg7Q2ZCbwKD5Ihn+/G0mi/Jy8+tfg5DH0+dVrob2gsgbhK8v8J/X9G3f7FE/V6Ny5j5tV4jeox18aYg4EPACuttUtJ/pt8ay/GHxU1Vvrd5s08tHnzeA8jqk/UFyEW+/oEvU6goZs83zkNTYZJQy5Bw0wAdLNGKtxFL2BWgLIGulzdACTWhWPpGhfAWDqi4vzQ51Em076wcjIetMcoY4+nzErc600PJhP63FPzZMhFFlIhV7tohbcylCILtXBtyfCKRrucEf24/Rxka+gX758sJ0MqpHsswyyMYe4B7QADQzHswJdWcsW7npgXGkPVNqTRPlFvTjJR7+qSLjrV+vRV9gEeAA4bZT+DwFRjzDAwLe0zKmqf1RV33MGuPXv4+HgPJKov1BeAPBUWVrk72gGS+2M9Qa+TMp3UKcs78SWOfwTQaMDJOaMB17HMYazlHiunWrumEsTKtpmLGYDjOUfD9Pn+PHl+0AEYeySBGNJJaq5NEXvcGEjc5K2PwvYn230XJuQ10jbE9XJxx3rrVrCglRKZnLDXIBd3rMFYh1boyXktS/FhIiFQpp0m3tNGw3LQIbD+CXXJ6G1scVWdbtqQ9XT90U7UmwondjGcMt1nrS1t0xizdjQdWGvXGWM+ATwG7AR+ZK39kaefi4GLARYuXMjq1as76mf79u0d1xkvxbH2Xr0e55bnnmPPnj1jcu4T5ZpCHGuv1BeAPCX5JdMrh0LqR3Jvmaq0TlQHbHs1Hh+sLz3JHQgQzsBYAmdDpGto1g1rd1G17Ys99oVzuLZK42/VNudwyvMC5i1JJ+IZNRY37tDKFRLA8ZfxTdCTZc2gcofTl3OTZx+ehGI890QbdpHgC/nJeukkPL2Vy7u5c0Lk+VxkCcbZhEzPUm/6VbhJ8rzn+jPVSq7HKa+A737Dc6kph8tegHGdMiEILqtbZ6JeWX9Tej9R72U9KhOUMWYu8AfAkcAW4NvGmHdaa78my1lrvwB8AWDlypV21apVHfWzevVqOq0zXopj7b16Pc6fzJ7Nzl27xuTcJ8o1hTjWXqlRXWTf1yTPefigUR+bkvxeqxPA7baOTJs6A/FwEAWwPrfY6ySrlw9+QeVL91fCsHKPZXplaIVuT8Lx4hI4Fu3kLlSD2hP1JORnE/Qa/vYkaMtznHkwzFiUHDdkOeloi358WzxltENfuObqJka+P0b3J987PJ8Tcc31ZwKT3XgtOqz42RzN59iX3s3fUl11Op46Y5qchCf0RMaY9wCXGmPeaoz5vjHmvb5y1tpdo+zqbOBha+1Ga+0w8B3g1FG2GRU1phrr7/Go/UsTHpCNMfPq2OAhGNagHKpXdzsa1WmjTplZc9KICg1VOcDS4CSdZA1O6mf3QnsBUNYAp91jeUbBEAvlJss6c46EqXPyZbOxiEl5WbV0DWftMOfa1S9dxMHyoDiHNL0x4D9HY2D2oTDtwPy55Mq489N9y2srrqduJ+f8qmuvPwM+p9gLxGIMof5zoTswKH7L6QQ2y/LGEog7Uafn4zQpbEJ3ozNJQhreb619HXBCD9uWegx4qTFmmjHGAGcB941RX1FRUVH7nCZ8iMUQHFrn2ycEsr3+gq0Lzr360veVPXqJLKAhUxy7/YYCq9zofVdMwZTXOa4DycI9zsYX2ErwNSZZTm36gWE4zoGxHI86lxxkBqQfKZ1VN2SPjM5CL9K+soeFDLT35x4Ju7ZAa08+3Zr0MssHhFh1nE7Uc+dhDdnaxzL+OAuxsPkQCyuuvy/EIheW4fnFQF4v36S9dN8Yy4IDYOPGwpXuOE7YV6aqXqifqm1Ze50e6/YmgTHGDFjrAtZHpU3WWmuM+V/p8e4etFmQtfZmY8zVwG0kIdhrSUMpoqL2Vb1z2TKaO3eO9zCi+kQT3kGeBseXnYQJ7DuVxRf2Ut201y1EH3WsKyShUTqAtMGpEFes3MKCe6xDEnxOpYZkipAlz6IsxEC7yC5t/uLiuYEHjgfy+Vn50LJvvoudlvGFWMj0XPkB//6BS/PnLK+N10UnX67gFpfdjOibB3WehbAKKF4nj3Oce9/I9wUcfVz55axSL/5WOm2jqrzn94va/aQuxMLORhTUpwCstd9Lj7/To3YLstZ+xFq7xFq71Fp7vrV2TGA8KqpXOmbBAo6dH5/NE9UbTXhAngLLfelVX1wuX1+Asf5yNhVlfOmdgLIB5h1AOSBpqMoqa9hS0IpsR8FVLmQDUU//VO+B8axveUYarEXagmNdDEkeiL1wrHDGB8ydKDsvtXp2LmxEHcv9gUnJI6/luVW56BqYfXHJuXzf+67fz0BYhf5lILtuvvfbd/NCsv526PKFs3p+E9lLhf6fqNN/A5gGL+zFOKy19wMYYxakxz/tRbtRUf2g+595hrufeWa8hxHVJ5rwgDwZjhuNW9TJBejkS7GbsdRJD+XJtML6x9krTVOOn9cpzMWrinoFgNIOpAQr8mV8oRa6bMhFlaEVk2eKcWtoJF9eHuuY5NEoc46Vm6wh2eeSzzwYGkNFiNblkoT2+eXO06jz91230I2Mwft+618VQr8U5H4NkJCcHM+ZX/yMjubz7Evv5F3spmyd8etPeahcA5gCnucMjkqX97i9qKgJr2/efTdfuvvu8R5GVJ+oHwC5ljPTyy/cTlWn7V6V8U/QUyAswVf/vO/tWLahYNgLa+pnfh13LPsrwBrFrQTG+ccU4TEHbfjhuFdgrJW5ye5YQ3KjuG8MHHh8e7+QH3KPNXDLevo9KgHcUgzUdWU9KLynBRgfm4l6naiqjW768N3OdaIpvZ9MN5b/ZUVFTUyN1f/zUful+gGQaz0kJKSJ/OfkG3v5BD2KwKTjS6WLWMhT8ctZO8p1zMr6PDaxXxaLLMfqyk1fKMYwUMzX6b7jsZA3vMIXPiHSB4ZgYEqaHsKv0HX3QLOG2dxNhHofXVuyTe81Uu9rwZVWYRrpvjGw4ICKa+Y5216V7WWbrsyAJ62TNibD4rrjqqlunoUSFdX3in8YUb3ShAfkIZgZylOIUFqmKm1vabR9H3yEa8gTfhCcoCfqSBfR91O/18XUkKYgORheQTGt4B4LzTrE34eG9azpvQDHWV8KkgvXSV1HY2DeUcU07ar73gPt6JeFZ+jrrX9RkGMthFnINuTnQo9ZtQscfLjnEhWTKtUJxI6lRtvHZPBckVFpPP+bioraJxX/KKJ6qQkPyGVP0XMaKwjW6NhpvdH269PsObqkBpqyCXp4gEeCU4MCLOUAVMMa+Xxf2TruqTGJ29oYzI9V19PhDHv757ZQOIkcm0wfmkHhvSi46sVJcPn31LePp10N2qH3UaZVTdTTNzft/fkvKLtQ/aXAJzd3PBkW9Ljbv+xxe1FRE18xxCKqh5rwgDzY5VP05H7dP6kqIO4WmMva6rTOlGnpTsDZy9Jy4KTAtxCrKkfkAWMffPn6yf08HxhTyDWdd5SnTmiSXh5OTz11bz4ATN8QeK6FTPOtaCG3hXZN/pxzfemtbktff73vca19bnLuVwb/uc+dR1B1PttlZcb6K9DXftnffJ2/+6EePk0PwFobZyJFRSlduGIF71vuXdgqKqpjTfgHhVTax/gxT++X1ev2C3k8vsgnTXY5AbgNglegxSoHuADEkIdB1a4vXMBbX/UxNCNN8jmpuo/8PdMvf/nL4nn1UCMjIwwOOnfbkD0ExI1F7xvTftDHjBfAc4+nadkJkD3Gwu1mWa69NCNry4g21DXN0siXR6S3O0ib9twEyfI+x1mkz5hVHHoddVOnqg1xNXvSdp00rclj8H+tMWYl8Fck4RuD6VCstTYSQtR+qSPmzGkvARoVNUpN+E/SvvQUvar+x6quLDOQW563CnJUDGqhYY+LWBqfSnu/LEZWj9wbOiDSh2b4yxRA2nPOwIwZCVyvXr2aVatWcd5557FkyRLe8Y53YNOn5B1xxBF85CMf4aSTTmLZsmXcf//9AOzYsYOLLrqIk08+mRNPPJFrr70WgC996Uu86U1v4txzz+XVr361/7rpc86ly9AJ32RFSs/J68LrsmU3I6Xy3KDovnVxdc5Tptbopseq+/cd+LSVthu4Favd76Sx+e/m68AVwB8C5wKvS7dRUful7nzqKW596qnxHkZUn2jCO8hlKvtSg/DdwViBc50v4lH3YQIHZQ6gr06ZC5zvsSRdu4yqH68T7Gl7qngykje8ILDv0dq1a7nnnntYtGgRp512Gr/4xS84/fTTAViwYAG33XYbn/3sZ/nEJz7BZZddxkc/+lHOPPNMLr/8crZs2cIpp5zC2WefDcBNN93EnXfeybx5nngC00gf5ez2m8V0dz6TZ8Ouzcm+dV6nOH+ZntmgHoi1rmya5nWORb9lYTjtgmKr6wfWsia/1JvMGYsZ5r1qt6odma/3qzRGTsRGa+11Y9N0VNTE07/cdx9bt23jxeM9kKi+0IR3kEejuo5TJ230ok5ZflnelKku3wenPghycOxzLn1gFHBDQ05jrqtQ26Ey8sTczMNA+8EwkKJOOeUUDjnkEBqNBitWrOCRRx7J8t74xjcC8OIXvzhL/9GPfsQ//MM/sGLFClatWsWuXbt47LHHAHjVq17lh+Oysfqu39R5+bw6dbx91bg++rjOtQs50jovc8XVPUtg1LkuapTpRZ3RaG/3V0MfMcZcZox5mzHmje413oOKioqK6gdNeAd5X/vS6mQ8vS47c5YoWAhtkK04YBZpGnhrDaomUGl49gJbCXwNDBXrVcUhBzR58uRsf2BggJGRkUKeTLfWcs0117B4cX4Z25tvvpnp06eHO5Jxxrn44lyhJH9olqcunnomX4+yMrp8FVyn+6GblDJI9rRpjG0b34Feq1zfsXKcx0Nj9P/Uu4AlJFMx3M8SFvjO2HQXFbWPy5gsbC4qarTqSwd5X4PmvaXpwRWhwQ+oPg9d/4xuivveK+xzqEvGoNPKYD4Uo5wrD1XhFd3oNa95DZdeemn2n+7atWs7qO2BVR/ku3gE701K4LpWXmtPucJ7GNivG8OcHRdvcAww1L4XiRobnWCtXWmtvcBa+670ddF4Dyoqary0v373R42NJjwgj+YPot/+mAYLvweY4m7lz/SeclWxzD6wKw25CLiUvnEMTvGU68Lt7lIf/vCHGR4eZvny5SxdupQPf/jD9SsHwxdKgD94bUra8jr2nnq1VedmqqIeYzeZfF/5u+1kHGM05l8ZY44bm6ajoqKi9m9N+BCLva19+ct5oM6SHmOqGsDbUV3aDwep02QA/rdv3w7AqlWrWLVqVZb+mc98JtuXscgrV65k9erVAEydOpXPf/7zhTYvvPBCLrzwwsBAvIOryFaT93L1Qj8ZVgUhjNWnNRR+kd8fmPC333nVvLXcmzoduMAY8zCwOx1KXOYtar/VH69cSXPz5vEeRlSfaL8E5E48sb2hXo2jI0Cu5TLuPac22Hdj3Kl/7ygEyKMJxK1Vd+ze4zoRL/tinPG+OKaAfm+8BxAVtS/poJkzYXh4vIcR1SfaLwE5al9SB3HLdeNuo+prDC9pt2/Xvgqo+1o4l7X20TFoNipqwmrNunXsevZZTh/vgUT1hfrsR9B62te+fHs1nmYnDeVm+oYq7o0rVdF3a0Qk7WvvXA/lDa9gdG9Bnbq1PgfdqRU4pVz3NdP2psbir2FM1n825svGmDnieK4x5vIx6CoqakLo+oce4urf/Ga8hxHVJ9ovAdlpvL+IRyPf2Jtj+stSLdrqsl5JuVazZn32cYCuGFsIkEvrVZ3vWF0P1a677ja/3+zgrZsIGs/byICWW2u3ZOOwdjNw4vgNJyoqKqp/NOEBeW+7OvsyghVZ0hZ3gxCpoUfuW8++bR/7XEhri216y1SMo7m7WG4MXU+f5FrJncl33WS63or9wrUpact2cD1rqeT9rluPeg5yN9pX/gb3gXE0jDFz3YExZh4xbC5qP1dcBzmqV5rwgOzT/vrn8fzzZbl1oEdCrVWOoK2Gpcr/mEoArgB57riZH5Nup0NY/spXvsLy5cs54YQTOP/883n00Uc566yzWL58OWeddVb2lLwLL7yQSy65hFe+8pX8xV/8BTt27OCiiy7i5JNP5sQTT+Taa6+t7Ct846DymyOe8hXnVHatc32odO97GNj33QxV/krQ3h+Jc2UyjdH/Sf8b+KUx5n8aY/5f4JfAP45NV1FR+75MnKcS1UNFt2EcZak/eadO2ee2kHeKQ9BZ6VjWdA2tFYOyxfTcVpxBrl5VOyaJHRmcnC/n8nz9BnTPPffw0Y9+lF/84hcsWLCAZ599lgsuuIA/+qM/4oILLuDyyy/nAx/4AN/97ncBePDBB/nxj3/MwMAAH/rQhzjzzDO5/PLL2bJlC6eccgpnn312+Gl6EiarXPvh7f684Pun+/G0me3Ll69/W9wv/CqAZyzkz1H127LlIRad+9ETW2NxLtbarxhj1gBnknzy32itvXcMuoqKiora7zThAXm8v0Q1uNYB2aoynYCz1PM7XF0f+Gjo8biJOQjygJOGqCDkepSD5YxyyVVyebkywO5tCSCXuakOlnVdoRtuuIHzzjuPBQsWADBv3jxuuukmvvOd5Mm8559/Pv/9v//3rPyb3vQmBtK18370ox9x3XXX8YlPfAKAXbt28dhjj3HssceWnLR3sO0xOz3/bDhPn6csF2q7KtkbMhM4lmlljnjuhixsYofUzd/x3v7bH+//a3xKgThCcVQU8L6XvITmxo3jPYyoPtGEB+QySfzyQacOkXRlugXUTsbTTX6tNio41Q/DngYKP6eXQVnIma4AuWywobNK03dugukLRB0DtgVmIF8OVLo+NVv5E5zMl+6wtZZrrrmGxYsXl9ZvV2gF9j3Xe/cWf152nVrtY9+vAYVfC3S5gCOsPwOlIR7qffaGa7TLND1h22MFmL1qt6qd0F9Dnf7HKBw7KipKaN60aTB16ngPI6pPNOFjkOt88YR8t731pTWaL/BO3bdWSyTYVn33z9uwx2nUP79nxwFAszYFPJ9brfv3lQF2bw2DoT5PHwymOuuss/jWt77Fpk2bAHj22Wc59dRTufLKKwH4+te/zumn+1fQfM1rXsOll16aTQBZu3att5z/3Hz7AnqzpexUXpXTG4onLg2tqfFe566xbt/XbrH+7p3+bsZSgduLyvJ1ynUKxFojY3d/EBUVleqmxx7jJ48/Pt7DiOoTTXgHeRgYqijTrfPjfMmQvzne8o1rZBgGJ7lcWRIFki69AhMciMmYYlc+5wCjynnac4AYCs0o9ONgrwUju2HSlOQcske0Bd6ZgIt8/PHH81d/9VecccYZDAwMcOKJJ/LpT3+aiy66iI9//OMccMABXHHFFd7L8OEPf5gPfvCDLF++HGstRxxxBN///vcD10zeetliugRlF15R6d5WpWlg9kCyvpnxD16UF8ehcBx9Pjapt2N7d0TYC4oMtVEXiDtt2/dLlNYe6Pmid8aYSzzJW4FbrbW397q/qKh9XT/53e94euNGzhrvgUT1hfoBkFsoJ1xjU5mz5PbrAHAVMPv63ZtgbYFdz8OUqVbQgMdR9cFNBjnypVunnZdBsIJjDWdGtyVgXUKsF77F1d78Ozjg2HyMsoRl2wQGKmORL7jgAi644IJc2g033FAo96UvfSl3PHXqVD7/+c8XyhUk3fLMPUeBqcjf+oioJ8qFXGDtHIfCJIIOvRXjEttcns2PI+RUB9Ng8ybfxSnW6KZMp4A71tBdB7x3w1h46ivT1/fS43OAW4D3GGO+ba2NK1pERUVFdakJH2KxByoXqS3z5kaj0BfjWH4hV+Vt3ypL+aDMA2zWBl4CngrAFACwMmDz9q38twIQpvl7nhP1QzG9EvybxRuDsZY8R6DSSR7ZnYRXhEIv5Nbn/hZCV9S18L3Huk2Zp2+SdPu5z4VMaxX62PR08fJMFPUawC2wGzZ3OZwyzQdOstb+mbX2z0hg+QDgFcCFY9BfVFRU1H6jCQ/Iu2FHKC+EDr4yVWm9VC/dMa0N61xDElZbnuMWeeDyuIm5dgJgXIBjD8DJeqF4ZQ2/Pgje8ZTqizwI6xCSQkjJGMraZCzyOOgkp9vNv3OFPXVa7XqF665iy7VrXPqe6L8KBcQ5WBblrCgj2y64y8n2ySc8l6iY1BPtjXe4kz58ZffAWARGHpY0nWkYONxauxPYPQb9RUVFRe03mvCAvAs26i+kquN+1kPZgk8BGNW3DdqZDDnIgM8pzNK1w5xzH0Xd3H4IwuSYRfnnniimFfalcywBcwzlg+PccauYbluJKy7zZf3csb7unhuWstCM0I2NfhVPLFDG8/kQx9bCk+s8zXWhsXzn6rTtygR+C6g9vt3wYM2inegbwK+MMR8xxnwE+AXwTWPMdOLSb1FRUVGjUj/EID8KHFNVLvRFVghgHgNZRr98m6+Mr84zTyWsYrwQpVxjrzusOxVwpCfqOSc6F/frRiZBLY039u1j0rQG7bjk9MzksSu/9XGYfWg6FheD7PoWS73REOOSxz1WLTj23GBsvE9cH1S+2IbeH+vZz7m86ubH50QXTyb/0u3JMQXDLJIHhLTah95eqvLqpndTthvo1rd0ncgCu+CuLqqWt2vt/zTGXA+cTvIX8R5r7Zo0+x297i8qal/XB089FfvUU9UFo6JqaMID8m64z8KrOkEf+SUnv/h6AbGdKtRmWV++PJk2MgKTGmlqAYp9gOMJhSiAlIYnBcYZQLfAmhR4rZqkJ8ErMElPgrN8sIiD5e3rYfpCmDQ5HbdJ8/UkPQ3JzfS4R7dDhZADArDsibveuSV5el7OUdZQqt8D9f7p91a/h2U3QlnoRkvAr8xTr6x/CeKIOnmg377Fc7nKLmXNtKp2Omm/qqyu41sSUv91hPpKAfnuDoZRW9baW4Fbx6LtqKiJphlTpsBQ1bpWUVH11A8hFnd0U899kXWzFrL+Eqw69uV1AgCdwIUFtm4iDzc+SAY/RBVeIl7ZFz6hwxly6QoQfesVF2KO1Rn5HpKx8e78uWVlK8IrbAtao5y8l12DZg04bhb3W0141v3abj119XsUiEWuDK8IvFe+1Svk+yxvgOSntfDeK4AW+4/9rjvA7QXIjpX0/xOd9N8EtvcwxMIYc1svykRF9ZtufPhh/vXhh8d7GFF9ognvIO+Ee31hEpbikmxy38n3xTcWax93015ZnbK8hx+CBS+gCKAFaNY/u0sgVq9cXqPdlnHpRiy5JtJJl2WTjnEulMK1Q7schgQrBpSLnKa1hvOhFq0mNFzb2kmWx7TTLOl4Tb3QC59jnOW1CMKydpE33tc+Z31ToEFYg3zO/ZcgbYvjy8Fw4KYmdxNUBN32e67PXcKxfsGD91Rcywp1A72jBeWq+rbGfkjpVdvQ2YhKdawx5s6SfAPM7mF/UVETQj9/7DE2bNjAOeM9kKi+UD8A8sN1XGDp+Egc6sZBrupntGshjxbOH7oXTn552pJ2ifXP9raVBIxmQBv4+ncQm7uSApJdGQe51tAOk2i11y9G7WcxyCL22G0dSEvIlaEWU+fD5Ont8/CuiZweO4CXMCziZr2gnAsv8EiDZ5am4Ti9Ljue8YdWZOArJvC5beZWq/cxB7U+SLb+MRQcZPleKgc5d2Mkb6J02faxtf4JeqMF2Lrt6Hes7rasPfUX07H2gLXW7qkuWVtLapTp+YNJoqKiovYnTXhABp4eYe89Ta9qq+tVtVuVVlXPV6cwUa/VAiNgp9VMwVPFIWewLJ1iDUqN/NYIN5lWHoY1HLt45JyjrNIcMIdikDMoN7DxLli4ov2EPVCQLGOOHWh6QDnLr4k/ITAuc3B3boUtvxXvmiujQyvk9ZBQrCA6l6b60oAr20WV1a60L6xGvkct6ZbLYwu2txP06kBpr8C7jnx9hUIvpIZ7DKvW2kd72V5Ixpg5wGXAUpJTu8hae9Pe6DsqqlvZ0YTQRUUJTfgY5NSaKRjBoS9beWw9+2Olbr7kO60j0/bsAloKirw/zavwAK97GDjWccrabfRCWwjQnHvabB9r9zRzTEUc8VN3wMie9jnJGGNvzHHaZqupILDsYouxtkbSsSo4ljHJ+hx3b4Nn72+PqSy0onD+PvdWXhMZ2uEB29x7o1+h91LCuQZj0W5Lu8iWZ58ufjZH8zmu206v1Ol46oxpN+zqfkTjqk8BP7TWLgFOAO4b5/FERZVqbz65Nqr/NeEBGeo9TU8HD0h18qVYV3WgoFfj8QFJth6ydPwc1HjdQgnCriENgeqWovRn+JAjqX/uF2egj3V67qEgAhKfvA327GzXKUzWEzAsr5JtpXnNBHxbTc9rpF0mdyORXpNWmpe7Ts32OT6/GZ65V+QpiA6FVuitL/Y4g2AH0+p98aapGyX5fubeWwfGTVHeFttQrvKvf0ZQnYJyVdk6fyeh4zr9uTKFu+8O29s1Nk/RA8AY0zDG9PyXQGPMLJIn8n0RwFq7x1rrWZ8kKioqqj/VDyEWPA9bLRygQw/c1oVAhL7E0h/js7Jlku11Gmtcp0wndcrybvkZLD0pLZWDIuH+NcRV8brEAYC2LuRBhFsYq/oKxBxXxSL7tnLN41w8sth/+g6YtximzmnnyXCKDDZDIRbuila9IRLa9Q2EBFgL2zbAc4+p6+vKVoRW+EIttMMroV9Cufd99N3Q+F4CeOUNS+7zIPtvZuOzLXjYs1ZDN/DbibqpU9aObm9EpOvbiqqxWGAXPDT60RVljHk/8BFgjzHmGeBSa+1lPWr+hcBG4ApjzAkkS8n9qbU29+RSY8zFwMUACxcuZPXq1R11sn379o7rjJfiWHuvXo/zPx18MJOHhsbk3CfKNYU41l6pLwB5O9wEvL6TOu7LzpB8AU5SeaP/nY56AAAgAElEQVRZyaIbeK5qo25Zt7/tORhpwmBDwEwWi5yWzMUiC1DLYKnZhuHcS65iQbt9CcZGOqMN/HDsxjJA+5EtrfzEPLeaRQbLCozl/rMPwNBMWLAkXdXCinbIr2xhabfnrp4PmnNOtOc2yxd73GrC03dBM33abwiOcwAbCC2RoRbaPc7eIx8Aa1DWoF0B0PpmSucVfpmwbH+ufbnK4LFboK2qp9vvdBuS/nkqBMqhc34Orq7oolv9GbDMWvukMeYg4O+NMYdYa/+2B20PAicBf2KtvdkY8yngfwAfloWstV8AvgCwcuVKu2rVqo46Wb16NZ3WGS/FsfZePR/nr34FU6Zw5Bic+0S5phDH2iv1RYjFVvhmnS9kq15OVTNouv2irfuF3kmdUH++vtc/Cpkj6P2JHA/8akiCwpXLldWxvDaflsGjhjblSur23diyMQpYBM9+WmbPNli/JnkYRy6/bpiFeuXCK0RdGV4h4Xf7k7BhTQLHuRAF2ueQu+4agkvOOwfA0jGW74F6D7N2Q45xiaucizum3Z7v8wTcuaazv4Gqz28n9cdC+v8L/f9Gy1dJaQTYDP/S+9EBsB14GsBauwF4N/DGHrX9BPCEtfbm9PhqEmCOitpn9ZPf/Iarf/vb6oJRUTXUF4C8Bf5tuEY535dw6MtvLNTJl30dCK5q85YbRSnt+NlWAm85V9EDR7KsF6YCYOydBNZU+yXAXIDrKkgWbioWNt0PG+9NJvDpMp1M0MsuoRhjFnds23nDO+HJ22Hro6K8KhOCY5emzzf33ihI1tfOdyNSAGXf++J5f7Prr95rGY+ciz1Olne7yz3k2Hf5yi5tjXJj/bcZ6qOsX31j7Sv7PNimtU92P6pSfQ74tjHm6PT4sKTL0csmY37cGLM4TToLuLekSlTUuOuWJ57ghnWedSajorpQX4RYWGs3n2pMa7ICfks+TMH92F730c6+tI7GVVLf5Y22TFlf6x5PQ41lfLB8qAapw4mhvSSau1WQsNpoL7mWhVZI0EpDMULhFdnxQH7fG2rhC62gHe9Mo52XhViIB37I4z3bkgl8A5NhzgthyiwRZiFwxkIhzMJm/4itvODpddqxCbY+Is5Xwq3vWAMuHrBV21xdDyRXucfBXwfkbaH4jJTdGEnX3R23LMO7Yc+edku5S1tyXDdPl6kDs6HjTvoJucT6yoW0HZ6r0V1XstZ+1hjzJHCZMWY5MAv4ujHmTcDt1trRxj7/SdreEPA74F2jbC8qamxV56FPUVE11ReADLAdnp0LC/bViXqdwG4d+drzpW15BuYdKKHUQbEE4rR25iobMQEv9HJAnAK0BOMW0DDknrLnJt3J/cIkPQ3HVZDcaLcF5B4m4o4hCXfYdF9SduYimL4QBgZpxyC7q0cenAsXPYXO5h547gnY+Uw+L3NxZXkPOHcch6yvvQRVAax1QVmu0lH2q0EublnVyX1e4KF7ygG127yytLL0ThWC3dFO0HsOxvSRz9ba7wDfSVeyOA44ETgVeC9w5ijbvh1YOepBRkXtTdX9VTAqqkJ9A8jPwS/ZByfqVdUvA9wyCK47hhuuh/MuhCweVoJxBjjaPSYAZvLVoOAih8DYTdKDdL+p9jUMN8S2DJI1GPvcZEQasG1d8oJkMuCU2ckT+YZmpNdAPFjE3Uzsfg52bkocaQnA8jpJPPK6yE3Pfh041qEoEoaFkxsC5dxSdR73OBfSolziHEgLhFSgbC387McEVffrqhPHt6psp9uQejBB79sVXfRE1toR4M709eW90WdUVFRUP6tvADmdqPf6UJiE3Pq+zOpM1CuD126AuCytTr3Q2KTWPwbDe2DSUKu9rJuDnww2JSyLcAsvGLuVLQwFFzkExg4I5aoTWbrHKc7KBiA5O1sPGEtozrUD2RJvciw7n01enUjGNGswli5yZXiFcJO94RUKklsKft2+jIn2OchyLEH3uFnM8/26kItbb0HLsvVZ2LM7K9m+HHUuZRdl94Zs4OUUCr2QGgGehe/2fnRRUVE+GcDEMIuoHqlvANlN1JtcUc7nPOkv6V7+eZUB8ViAsQ/Yb/8VnPwK/A5yAZbJ5/mcYxmLnFsdIQDGWboOr5Awq9JaJLHSXkg27bQMen1uMuQcZRcq4C6UXubNKQv9UFfbakxyWSEXuZvwCnkjoiFZ5ynHOAfK4r2tdI9b5B+lI/IKYRXaQYaffN//d6X3u4XgOs5yWf91+q7Tr1TdCXqtZHWJqKiovaC/WLUK7rkHHngAjj4aBgYq60RFhdQ3gLy3J+qFILcKVrvt19dunXFC8nSzlaeD0ZP1srAK6SS6H5VDccjNvH3WGPQAcgCMrRhhti9c3xwU62MRd5x7YAhiP+277AEhkI879sWr1SK5EhfZ5xrL/VIHWZWRcIwnrQyUc3343OO0nH5yYKG+gm0H2K0Ww3tg3aPhS1gFnqFynbjFvXahXZmQS+y7qfZpLCfoRUVFeXTkkfDzn8Of/zlMnQrHHw/LlyevI4+ERl8s3BW1l9Q3gAy9n6hXF2zrltf1ytryle3UcXb7zSY8uQ4OOtQKMHbQ0wIjfo7P4K7EQXagLX9qlzHImYMsHWWRT0vsCy/OAa10s1sIkHehGTZNkw6yBGPIHgZiRPxxmXtc9rOcjL/1ucjBCXoKUL37oYl4PliWDrLYL7jBCopz5WW+Ox/hLhduhkZUPQHawJ2/LsJiXfe4E1gtc4ur6nS6ddrXJ+hFRUXl9cNZs9j4utdx/vHHw513Jq816fqTM2bA0qUJLJ9wAhx6aPn/+1H7vfoKkEc7UW8P+Yl6vrJ1QyN89bptq6yMrx0fZP/4OnjnfwWTQU4Kxtlybk1oNRJozcUiu3hj6SbLMiaFYhmPLCfhpa50Y5Bc2AWQj0GmPS8PVMgH5Jd3E3V9+7UcZOUe5wjH6ATPhfc4yDqtEJMccpM9E/IKQOycYe0SWwqhFXpinncVC/lAlEBcsgzLyGC7fW62Bb/6acklGkVaVZudAHOncv8n6LXVQ6AcOqe9NUEvKioq0R133MEjjzzC+f/lv8DppyeJmzbBXXe1gflXv0rSZ8+GZcsSWF62DBYtisAclVNfAXK3E/XcsfyxPOTilrVbVb5OW2VpVW5yGSxv3gS7n4cp0yUYtwToQm5yXuYiK3fSucgogMvcY+kcaxgGf5iFgmRj05dahs6V87nGuX0HdXTpIPsc4jS9joNcNlkv6CBXwHEhBtmB8ghFaNYu8gilK1n43GPvpD0ZamF5al3y64S+YnWd4TpObqdusa9OGciG0mUd/Y6HQi+k4gS9qKjx0ebNm/nMZz7DjBkzuPDCC2H+fL63bRuPNhqwYgVs2wbr1jF382bece+98POfc82GDWwYGICDD05A+eCDOfCoo3jzm98MwJVXXslNN93E3XffnfVzyCGH8IY3vAGAr371q2zdupWhoSHe+ta3MmvWrPE49ageq68Aue5EPedb+tQiWy14VJP1fFBeBdFlZQnkher52vjJv8I5b0aBsYRKBUgY1YrY17HHWbyyhGUNw8YTZuEKCEjOSS8F51arSIE5A1+ZDoVVLawGYgg7yBUKAbN3eTd5rEImclDqg+UyOA7tS7dX5Mu+NPz63GP9Wchik5M8a+EH3wnDqNwvu7Sd5nUKzN1u9fJusv04QS8qat/Ucccdx7333sstt9zCvHnzsvQHH3yQO++8M1f2oEMPhY99DNav574PfYiH7r4bbr8dbrkFgCMPOACefBKWL+eeNWu4//772bJlS1Z/x44d2f4dd9zBhg0beO655zjuuOM43bnXURNafQXI6US95mTFvxog06//rJAVeXuAqSV16247GndJnW5dZJ9+ez/s3AFTp7cAHYsswyrU8m+Zw5iCqHaOgUKohZvs52DYCjjNOcsuXYVb+Jxj3ANJDEU32Q3EB8aNABAraM5OxeTLugpB95giGJc6yM2StAAcS4AthELIvEBoRQbBnjCKXDseYJafkdYItCxPPgHbtubfKf2uadVxeX2q22YvtUe0H3KTQ+OywDbYOkZDi4qKCujcc8/l3HPPLaT/2Z/9WbjSwQfz11/+crJvLTz2WDsc45e/hH//d/4nsH7yZBa99KVJOMayZUmIRqpPfOITrFu3jve85z00m/oWOmqiqq8AGWAT3D8XjtdQGXJZtUZEei8cZN9xHVe5UzD29Y/qxwLfvwrOexeYHBjLB4Uo51DGJLv9LJRChVrkoJh2WmMQf2gF7ZCK7PHXvifpOYcYVUa4xhqUfRP2CmEVeEAY/wcjywu4xy4v5CCHgDhzdrVTHIBj78uKuGP1fujQCg3DBSfa5x7LEIsk9vi6b5ZDbh2o9dxulOZ10nadsqHyoRe0b8OqxvAMXBXoMioqal+VMXD44cnr3HOTuToPPwx33snu666DG26A669Pyh5xRHuFjKVLGUiXlBsZCf3+FDXR1HeAvBE+9kL4WtXqhxq5QhAt8zuJQ65Tpg401xmDrk/J+Dc8Adu2wKy5lrCLLF1VvVX7WaiFC0aVUOzSRNhFtpqFGmAureHfFlavcG6yB5QLYOx5kLjLczHKufNT4NuuUEwvuMcirdMVLBzEeqFYOrvaNfaEROj6eim33H5gsh6qrrX87n7Ytat4RcbSPS4rF3J068Kx3ob8H+vJC91kPwN/H2gmKipqoqjRgKOOgqOO4tG5czny9NPht79tO8w//CFcdx0Yw7wjjuDvTzmFQxoN2LkzWWYuakKr7wB5M3zrefjqTIU8GkCbJCfvA04Xx6zrdLvVY/CNq04ZX7tV0Iyn/e9+Hc5/n3aRHSiKWGIaxS2WYKhFwyT7bj3iLE2EXTQGKYXkbIKe7+UgV66FLMMuJCiDf4KeO09xBY0ppzQg5zS70IkCRNt2XsEdFmBccI0VJPvCJ3LALNM0MFeVq3j5oFm4x60W/NAz9azXoFzHCa56y0J9V9XbLcrp/urA9xYY3m3tYx0OLyoqal/X4CAsXpy83vQmGB5OHkpy110M3Xkny267DX796+QBJS96UdthPvZYGBoa79FHdai+A2Rr7fBSYzbOgAM1GBbKkv+ic7A5DAwxuhAL12YV7JbV67R+XRd5y7PwzAY4YJEl8btEeIJtki3pJkMsgqEWqZsLKh453WZhF2lWEJJ9+IE6TuHYWjAtMPLx0wKUC+EWtGHZB8TZfug3Al9yTfdYwrAG5lzYg4TlkvjjEPT6YLkAxk3x8kF1M71mqu8MlC13rQH3C2I3UFxxVYNtdQPMdR1keWzJroD3VdWuBTbCjYEuo6Ki+kmTJiVrKy9dyvB553HTT3/KUa0WBz/1FNxxB1x9NXzrWwlYL1nSXlJu8eIkLWqfVl++Q8/AFw6Dv66CS8lnvi/BOu6ur11fWlldKsr4+vb1FaoHxbrXfhPefQkY96CIQogFqqYv1EKUySbjiYl7WSyyDLGgCMnaKXZ52WhdiIUMqwhM0pNLvdl0bM5tdq5xlq7PqQa2+eKOoQjKhWNfzLGGZOXklrrJCmB9YRQaoH1PzCsNsUiB2bag1aLVhBv/LXhlcvt1bkjLypUd+1Sn3Tr1WyVlQqEXusxG+NsOu4+Kiprg2r17Nx//1Ke4+OKLOfj88+H885NQi3vuaYdkfOMbyf/9kycnrrID5vhY7H1SfQnIT8H/2QN/PUWla5B0U8B8yCcCCmqFTdRxhvUYqoDbV6YKqH0usq+953fAg3fB4uVWgLGrURFqYU07ryVGkE3GkyEWg+07EXcMKSTL/xAE0rg8F3PsC7HI8gQMm0CIRQaqDorrvFMB5YBYjLvWWsginMI5zB2tf1wHjtUqFoWHgoiwidLQCgHMaXjKf1wPLZs769z+aNzjOmCrb16r2vS1WwXSvtUrgGzlm6qxboPWdmt/HhhOVFRUn6qRPsZ67dq17N69m8MPP5yTTz4ZVq7kXx5/nObpp8PKlckqGY8+ytEPPsiK22+nZS3fefbZ5Ml+boLggQey5NhjWbp0KXv27OG6664r9Ld06VKWLFnCjh07+MEPfgDQ7jOqJ+pLQLbWbj7JmOcXwbRcutg6RNJf9C6v7lP1qtLLXOVOytRtO9SGr+6/fReOXAxDk1skE/YMycoUDoQDoRYZJBsyFzfnHNN2jJ1bnIPkdKVp6R5rWG6RjxPOtjLmOAVLB8pZ/LMnxCKDftFujnCqbmvkoRiPXgcZ8gAsHWPtHuecYwnL6dY30S4IxyPF8jkgboo+PaEW3tAKB8uWzRvhnrXlV6YKlDtxectguG5/VX34xjfiSfOdo69tCzwDD9bsMioqqo80efJk5s+fzy233MItt9zCmWeemcHqV77ylcLqFueccw4r3vY2WmvX8uU/+ZPkCX+rVyeZg4Oc97KXsfTCC9n9ohfxZbcMndD555+fAbLLnz59OldeeeWYnuf+pL4EZIBN8P2D4M1VYRYhSBY/kNcOs3BlfX1UtVEXfkN9VbUbAuhvfwnefjEY7SYWQi200vxcHLJ0jtN9CECyak7DsnSPC26ywQ/Kabp0kbP4agmybvx4wix8cm15cK0qxKJOHHLlyxOLHITjZhGeJQi7dO+KFjq0opmFVlz1RT/g1oHUTtKqYLhOmSoQD0F3GZTr8Apf/RawEf5XYHhRUVF9rIGBAb74xS9m6yA7RxngqquKqz42Gg0YHGTgjDO45tZbk0T3WOy776Zx113w+c8zw1quWbAATjwxeS1bBtOnZ+0fcMABXHPNNXz1q1/lercEXVRP1LeA/DT83TC8uWreqAuzCOUNEgbeEABrheC0qo26IE1JuRBAu/1NT9UJtXChFBXA3DB55xiKk/McCINyjwWeaPdYH+sl3nKT89w+nhCLdJxyLWQr2y5Rbr1kAb4yry4k156gZ9PraIuhFbXgWIdjSCj2hF2UhFbs3iMuoe/yEM73gXQVoOq0bsqE2g+1tQe/LP7H12g9D2yGbwSaiYqK6nMNDAxk6yFLDZWsYGGMaecfdFDyevWrk+MNGzBr1zJ0223w85/Dj3+cLD23ZAmcdBKcdBLm6KMZGhpiaGgoPqSkx+pbQN5p7b0vNaa5YBRP1dsFTE/T67rIZW5tFYaVQXOnYBxq13cu1aEWMpTCCPgy+RhkyK9/7CDZQbA3BtniDbVwsFtwj2XcsQRi4RRnIOxA3MG+c7ldmrqYpb8NyCSBSRn81oDkWrHHgZfPEdYwXDjW4OuZkJdrUwN0eWhFHdilZjmfa1vVrq/tsvbK2rAkq9fIY/mq08azsNFaG+LsqKioqM7kgPm1r02WD7r/fli7Fm67Db72teQ1cyaceCLLZs2i8fu/P94j7iv1LSADPAO3zodTNBD6QFR/ocsv2jrucRnc+lQGv75+y/KqYFnmh+C9GGohJ+CpUAqbOsXZQzugMP2/zDn2gvGgSmuIY+kei5jjHCjLfecQO1A1bYiGtsMsr1jORZZXjTwQ+9IKbrIDXwHN2WQ9BcaVLrJzjD1xybXhWMYcizAK/cq51PVDK3oBynXK++qXtVWnDbeVDrH+/6AVqKP7eQYurxhOVFRUVHcaHMyWlOP882HrVrj99gSW165lxebNrAC4+25eMHUqzJkDxx0X118ehfoakJ+Gvzsa/rXqJEdoR65q7QY6eR5OCGx9kFpWRren9335ddN8+5ueggfuhCUnWBWDrEMr0uNsAh7FsAqnHCQPtPchAMYy5rhVdJMlFGew7MJA5L4D4nQQWciFCKvIxSAbP/Xkrh7loJwLp9DHEoK7cZF1aIUnzrgSjpVT7JuMJ/Ks7T60oi4U+yDb10ZIZX1XjVen7QqMpUV+4l6ojT3A0/Dx4GCjoqKieqnZs+GMM5KXtex+4AGev/lm5jzwAPN/9jO4774Ejpctgxe/OFlB46CDxnvUE0p9Dcjb4Ic7wM5WvOhzYOUXo3OMGiRfjlXgquVrvywkw1cnlFbWhx5jWR++8/nRtbDoMJg114cFIUlYTsMrpJucucFqXwKzBmNIwLUF7QeCOCBukJusVwizsO36LszCwbEDYx2DnJUPKBd/nFaqXL2CNpTqPA3ABRdZxRmHJubVhmPPZDxv6IV7WR75Te9CK+R+p3XLIDpUv5Nt1QtVNjT2Z2Fn09pNFUOMioqK6r2M4fv33MOXrr6aq6++mvt++lMOmDu3HY7xhS8kr4MPTkB55Uo4/vjkQSdRQfU1IFtrW8ca88gsOLIqzKJJcWEFJ/dkvRD4VoGw7LcT+JXH1Cjva1/n+9qR+V/9HPyXP4fBIQnJITdZqDGYd4ezdAHCMv7YlTMN5SKL8I3MNXbbOu5xi1xMcnacnocOtdAXuUyVE/WsZ9sq7tdykRUga/fY+9CQKjj2Tc7TrxbbtsD3vlkOqXVAuSzfV06mlbXbTb0y7Q6kW9pxyboPqRawCb7XYbdRUVFRPZNb1aLZbNIaGoKTT05eABs2wK23wpo1cP31cO21MGVK8qASB8wLFozj6PdN9TUgA2yA9x0O1+sT9QHyAP7Jertpr4lcx0UuA9I68BsC3U76CdWtguiREfj6P8H574cGctKefimnGNqQLGEZ63GOUxjOOcYKiGWYhY47LnOPczDsnGvnJis4ztFO1a2NPBQoFpqk142L7H16Xh33WEy4qwXH2kluZnA8Mgxf/Ww1bPry64Kyr05ZOz5HN1S/CsZ97bpfiVrpy1Lss6z9XcB6eH9J11FRUVFjKrd6xuc+9zmazSarVq3KjoeH01v9efPg9NM5fvJkzmo04JZb+LRbN3nePDjsMDjsME589at5+apV7Nmzh3/6p38q9HXKKafw0pe+lB07dvDFL34xS3/Na17D4sWLx/Q896b6HpC3WPuDk40ZXqie+yEBWG71F7XLa5HAs69OJ25yCGZ95TpxkX3lQu35wFqmb9kCP/g2vPbNJJP2Wk2y2F43QQ/a8cdSbrUKIxYPMdZznF5d0xCOsc9NTkdVAGaTtOlzj1HHGRQLYM7eaI+bXLioCo969aAQCcZyW3CNU/DNQbRvYl6zJhyPgHUxx23n2LaSm6PhkdzZZvt1nN1Q2ao6Oq0OMJe1Fepfb8sCiXyT83z9PAnrhq3dWNJUVFRU1Jjq6KOPZuHChdx9993MmTMnS7/jjjvYvTv/O9msM86ACy+E97yHtW95CzzzTPJ65BGwlgNXr4Y3vIHW8uWsvfnmwkS/Qw89FIDh4WHWrk1i8TZtSiLMIiBPMK2DfzoA/iQUQuGkJ+tJQJZLvtWRD0I7gWYf0Ibar+sih8pCsf3fPgB33AwrXmopPlKaNih7ITnduucSZ2EW8thBsKUYf6zDKwy58ItcGY97XAi1SKE5F26hzrp2DLItpuccYrcV0KxdZA3GhVUsPE5xaBWLOg8MyT1Rb0TlJcBtLfzwatiyuXCWBeCtgt5OQFpvq0C3LM9X39evLFtncl7ZGJrABnhPoNuoqKiovaLjjjuOyy67DIDV7ol84HWAMxnDFd/6Vvt4x45kZYw1a2DNGqb87GdcYQy86EVJGMbJJ8NRR2Xfl3PmzOGKK64A4N3vfjetloaBia39ApA3wP94Ebx/puJAnyusvyjdsfv5Va92UQa3nQCxzq8KxwiV97VbB5J9+Tf+CA46BBYeYklQQPUsn4qX+7uw+bAKCBxLMLZFN9kXVmFa5Bxlr3ucnp1e/1imZeNw8EyFbBGUa4daSAgOTcwTeUHXWKV5YVjHJI/k0z0rVmAtd90CD93vPetKOB4tKJflhcC2qkydrXSIfS9dxzfmTTC82drve08kKioqaiJp+nQ47bTkZS389rcZLPPNb8I3vpGEYpx8MpxyShLDPHkyAMcffzyHHXbYOJ9Ab7VfALK19vljjLljBqyogtERkgl5Pg0Dkz11gv2KclTU8cFvCFzLnOWy9kLjK+vrqsvhgvfD7HnSUxOST8XT6TKsIrfMW5rnJt1lYCzAWr5yccgpKDdlXLIbg1vSTbnHofWPfeshl14tqBVuUbYWss9NDr6Ea5xBsQLg0nAL5RxLOBYPA3n0IVj9w3L4DQGi72ZS1ytrq6z9OjBeNT5fHaeQe2ypt4ZLC1gPn69RNCoqKmpiyRg4+ujk9da3Jusur1kDv/41/PSn8G//loRenHACnHIKl1xwAcyfP96j7qn2C0AGWA8XHgG3a/jV4KmdYivy9tCG505dZA2zZXV846qbFoLxOjBMIP0r/xfe9acwY5YHkjMXWU7MQ4RRCIQJLfmm3WQNyG40Xte4IdqR+849lq61dJRpg3OW55Fvibcs3VIOyCkk52BZ7FcCcsA1DsFyyFGWUKzgeN3DcK1nxYq6+zLN10ZZ/brlyupVlQnVCU3Iky9Zz/fD4Q6w6+EvPFlRUVFR/aXZs+Gss5LX8DDcc08CyzffDLfckpR597vhDW8Y33H2UPsNIG+39o4Tjdm+CGZoUPVBssMyXWaEZLZfCHy1qkBYlqsbVuFL023oPgns67J4yrQsfOnTcNEHYdoMDyRnE+4GRU2busQ2X64As66cB4Rl+EUWQ6xXsmipNgUYu3woOsn66vVkHeQUvQqATBt25X7o5QA4N1lPTM7Tad6l3WSazz22PPkEXPPVMER2As1V9erW0eXqwrCvjbI6Zc+DrhN7bIENcJe19vmSpqKioqL6T5MmwYoVyeuP/xgef5wPv/a1LPr613lvBOSJqXXwtwfBJ0Jg6eQm6/mWfNtF+6J14gj7+gnBb7cusk8h6NVt4dmXZZotuPyTiZM8faaD5LS0SWGwIXq0FhqtPDTnwircPgKcPWcnV7IoOMvkXeUcLEu32LUlnqqn45B9FJS7kArVrNuKSNbgWsjptjApLwDGeqtdY1qe/VD8sYo5TuF4w2PJ48XrwHEdUNZlfGVD5crSfG2UpfnSdbst2usb6/6ds1zWPmn99XCBJysqKipq/5ExcNhhPNtqMe35/vIL9itA3gif3Ar/OFc9E0R/Eab7FbIAACAASURBVBp17PYdtvXSRR4NEJel1YFeXa+qfLMFV3wqiUmeOUf86CyXbXO1s+XexBUIhlikwJsBrgJkM5AHYAnIrqwxJI+Vbok2lWtsm6KsPEM6cJDFDUBuH3IOsdv6Vq+os9Rbbk3kgGusHWTvA0O0s2x54nfwna+NHRxTUrZOubK6obbqtgntB4NY9QI/HPv63gjbt1t7e6CLqKioqP1KjUaDVrM53sPoqfYrQLbWNo805sdz4NVVLrJ7ep6ERFemVy5yWZmytsvSdLt49vU4QuP3lW+24EuXwvnvhTkLUrjLAbJozdgUWNM8r4Ms9nXMMaSALKBbxx0XQiw84RauHZyr7LlqXYVYBPZHC8iFGOSQa6wcZPmUvQIgJ30+/AB876qxgeNuANgHn55PUmEb6isEzq6O/O1DxiEj8uqMex38nScrKioqar9Uo9HAxmXeJrY2wLuPgsenqHQNq6Evape/B/+KFj7wdHUJlKkC607TdH6oX33uvrFAsc2Wha98Fs55M7xwscU0miWALANVWu3VLHxucgGQLdCAhsQcGWqhwyoUIEu32Pc0PaduQiyyNAHIOQdZwrKlEJdcFnNcutybdo1FGW8McvoQEAu//incfGNncBy8HIHyIbCtA7llwBwag+5T52mVPVZausdlbW+F1kb4ZKCpqKioqP1OJjrIE1+7rH1iqTFPHw4HakjV0DlMEkrhi0WWK1qg6tXZ6jpSewuSy/bx1NVt/uu3YMUp8PLXWIxpkscJB8jSQW6k8cQ27yCTLuXmBeQ0NEKmS0D2hVXI2GMdgwyMbpKezadVxSB3MkmvAMq+h4No17gKkFu0mnDtN+Cxhwtnke2HILVqvwqOdVmtqrQ6YF1WVtYJuccyr87YNsBPrLV1VoKb8DLGDABrgHXW2teN93iioqL2Ta048ECmzZ073sPoqfY7QAZYD39yCFylT14DpUvzHRsSN0o70bKdUNuh0AZUnobTumlV/YfKdgPJa38NG56AN70LGjkESbHDOca2lUBxQ+e51lzohAij0ECcpRvhRFct8+YLs1D7taXgOINhl1YTkLPJeuoBIRqQJQy7PA3Fpesgt9izC77yOdixrXAW2f7egGN9FX1wG+q/Cpir+nTt7vSMxcnnHvt+LNwNrIeLSprqN/0pcB8wa7wHEhUVte/qwqVLYVZ//TdR9fTlvtSz1n7radgV+iKVGib/RencJivyfKBQ5nL5YED3XwcOdNpoy/r269R5cj1c9n9g1/NAKwW35rD4iX9YbUdIlh6T5WRIgM7Tr+F8W4U8X7qv7ZplZbtNz7nI8wnVdfWae8Ljb8o+VJnQtdDp2XGLzRvhC5/oPRyHPtt1ylal+erqtn15vjJyKx3jMve46tyegId3WftEoPu+kjHmEOAc4LLxHktUVNS+LdtocNcDD7Dmqqu4dfVqdmzfPt5DGrWMLUxA2j+0wJi3rYRv+Cx0F1LhPMbB9OWOG7QfJDIITBV5iK28+5B5umwor9M0X3tV+b4yvv06dRoG3nIRHHiwS2gkDxBpDJCFUTQG206xDKuQZQqPmU5dZJcPAefYJH3KkIpSB1mfQZVsSk4+B1mtYIFNbhaQDjFFtzjbIo5lzHEo3CLsGlsL998J/35t9Y2a3tdl6wB16GbMl65XiZBly8r5nN068O1AeAd+OHY3uk2V5nOPnwduhqN3WPtbT3bfyRhzNfAxYCbw574QC2PMxcDFAAsXLnzxlVde2VEf27dvZ8aMGT0Y7dgrjrX3mijjhDjWKv3bhz7Ebfffnx2fefjhrLrkEnYeemhpvfEY6ytf+cpbrbUrq8rtt4AMsNyYLYfCbI1IDrEkDA6pdAnJ01UenvohwK2C3zLwDcGqD8zL2hotJOu2l50EZ/y+Y16TwLAD4ka6H4LiXJoEYzfpzrMWcnD1ijTEoieQ7EDYA8duX8cjFybrIYC3DJQ9gFy1/jEtaFlGhuG734D1j1aDa+g4tF8GotTMD/XfCTCHxqzz3P4ISXhFyD12Dw1pia3vhuIBWPMba09mP5Ax5nXAa621/9UYs4oAIEutXLnSrlmzpqN+Vq9ezapVq7oe595UHGvvNVHGCXGsVdqxbRvr7rwTnnqKj/3f/8txW7fy3xYtgte/Ht75TpjiC0gdn7EaY2oB8n4Zg+z0GPzBC2C1fvx0i8RFdl+SaYRsMB5lJzCNImbJ+lbl6y9gX56sJ+vLNALpOi1UL9SGb1/WcXm+tu+8DR64B97ybpi7IEWUhsKShoo1LsQZt5IWGwPJ1oZcZUMOikOucS7m2B1DPUi25JxjuQZyBsoCmPXKFWVOctVKFhKGg7HHiWv8+O/gum8mS/GJkef2xxqOKcn3QWfdtLp5vj53Beo7eK7qAxIH+gl4rSerX3Ua8HpjzGtJplrMMsZ8zVr7znEeV1RU1D6o6TNncsxppwGw7JFHOHD2bLZs2cKtX/wiXHNNAsrHHJPkL1vGgQceOJ7DraX9GpC3WPvT44xZfyQs8sGtTHNP13OQ7BjEoVyLHH5Vgm0VEPvGUgeSdXv6XLqB5LJyoXq7d8NXPyvcZPf8MgnALtQCS8FJduX0w0EKMOwJsyi4xgqanTqeqCcgOEvyhFZ0sxZy7il7EpBDK1iUu8by/XH7ZdCrzjJYbrRwXKes3vrSQmV0O3LKqO8VCteQagGPwvW7rN3oye5LWWv/EvhLAOEgRziOioqq1CWXXALAfffdxye//3247z647TaYPx8OO4y//Ju/iYA8EfQYvHoR3K3Nf+0YG5I4xTIXeboqHwJiyH8R18kLpflgt64zXBeSfXXK6sk8v5sswLdh89BciEPWYRYy3EK7yKEVLDzQnJ1EXUj2wLGGYeka+9zjXLrnyXmgjj0hFRks13ON3XEIJnsNx1V963rdAHNIeqyWcvd4OFBPaxvYR+DNJV1HRUVFRSkdddRRfPGqq2B4GK67Dr73PZgyheb69fzqpptYtnz5eA+xVPs9IG+39p7Fxtx7NBxX5dw2SUIvfC4ylC/7hmpL5yHK+OqF0nwQPl6QHKrr3ORjl8ErXweDk1rk10huqJALFV6hwyxslYvs9j2hFtm+HC2ePMjFGOeugsqToRV6uTfv5DwdcqFAudRFTtrb9XzyRLwNj9cH1jI4rgvVdfqqU69qbKF2Qnm+tnZ6yuiXbkurBTwCl1trd3iy9wtZa1cDq8d5GFFRURNMQ0NDbaf4ve+Fc86BSy/lpk9+kr/fuJFPX3/9+A6wQvs9IAM8Cq86CNbpeZTuS1QC6Aj5B4TIL1r3eGoHzL0ItahK0+MdC0guq6/rhepa4N674L674OTT4SVnWBoDKUQaQ8491uEVOswC8iEVuQl8KrRCbjuOPdZScBxa+9iVlRP1nEPsC7nwuslWuchJneE9cMP3Ele+LuCGjuvU7RaOKUmr6q8sLVRG99mkvTKFC4GSZeq6x5uh+QS815MVFRUVFdWJDjsM/vEfGfjwh+Gf/5nm+vXjPaJSRUAGdlm7/oXG/HQJnKFDKDTMui/bhqcMJEtB9TLUoo47XAXDo4XkUB08ZUOgLPd//XNY8wt45e/DcSdZGg0Nyto99oVXCFguOMoKknUYhfdYj1qdjX6anu+4chULW+0e661N6rdG4Jc3wG2/CkOh268LxqHjXsJxWdt16vvyfO3KtixF91jut6jXfhN4BD5qrR32ZEdFRUVFdSpjGHjRiwBo7t4Ng/suhu67I9vLegxefwhsmaU41X2pShB1j6DeG6EWVZBcF4a7gWSXXsc9lmVDkCyPWxZ+cj3c+CP4vT+EI19kMRKUCy5yCsQt7RjXWf/YHcsRIJxlFCQrFRzjLKOdrifrla6H7EkruMlJ/VYTbv81/OLHyTXLDUsP0zO6unBcVrcbOO4WhEN1Qnm+fstCK9KporX6eQaGn4S/JSoqKiqqZxoYSn6Hb0VAnhgasfa5w4z59lJ4c5mL7ORzkV1+t6EWIZjFk7Y3ILmsXzzHobpl7Q2PJLG0Q0Nwxmtg8TLLwGCKMg2D1z22HjCWLrILpfCFW2DSJOUaVwKyQjHpIIdWsPDCMZS7x0m7e/bA2pvglp/lJ+Dp6+qOy2A3VMdXvgyqu4VjOkzT8FsG27qtTkIrXB11ebN2HoX32f15ofioqKioMVAjBeSRNWsYWrp0nEcTVgRkocfhjw6B8+apxSrcl6kEUr3sGyIPugu1CEFyFbzqtF5AstuvC81U1JXlfW3s2QP//j348ffgRcfCK34Pps2wGJNOWjMSjCF52IjxOMol4RWlkFxXHcBxZWiFDMGwWAtbnkmuwYYn6kNuGRzXaWO84LhqLL7x6vKynqXz0AofHFvgSdj+tLX/XNJ9VFRUVFQXOuplL+Njr3oVh//qV/zzV77C3ZdeCkceCUccAQsXcsKJJ/KSl7xkvIcZAVnKWrt7kTGfnAOXaBfZOcZjGWqhYbIOJOv2KCnbCST7xunrt8pNlvUlOJeN/cH74KH7YPpMOPtcOPSFaZwyLeEMNz2OMtWQ7Hrr9EEhQGmoRWHFCoFtXlBug3ZzBB64G278YXKj0AswLqvnqxOC27J+OkmvKlsFxz4Y1u31IrQCkpvfR+EtnqyoqKioqFFq+qJFLP3a12itX8+//97vMWfjRvjNb5LvxEmTmP7yl/OSWbPg2GPHdZwRkJU2wH9bCO9ZBNPKXF6n0YRa6DTZhi/PVqTRg3Rf/3pf1wsdh8rXBeXt25IHYDQMHLUYTn45zD8wXf0CyJxlUCDsQizUKhZVcchVsvoHewToShgOwbHNqo40YcNj8OsbYd1jflALpYXAMVSvDIz18VjAcajvEByXnWNofL0KrbDAo/DQZmv37fWHoqKioia4GosW8Ref/nTyqOkdO5KHidx8M9x8M/bSSzGf/Wy7cKsFO3cmr+efh0mT4KCDxnR8EZCVrLWtWcacNhfWTtN5FCG1V6EWOm28ITlUzge9ZXXdvobiKlCW7bcsPHR/8jLAghfAS18Bhx5lGRy0Kfc6d5l2qx3FIUvpM5VyAOwOdUiFhGOystbCnl3wm3uTVTy2ba0PxS69CozrHteF21C90UJzWVpo/GVgbeldaMUWaP0GVnqyo6KioqLGStOnw8tfDi9/OV97xzu46ppraPzgBzAyAs1kqdNvv/jFDDUafPaRR/jl5s189sYbmfXCF47ZkCIge/SctbcfZsxVS+EtvQy16AR0fYgWCl/QdUebHiongTY0hqrjMlCWZX31XdrGJ+F730rSp06DFS9JHOaZcy2Dg2CMGJ3xwfBYhFiIfBJuHhmBzRvhwbvhrrXJGsZlABxKq4LcUJmyOiHg9u3XyauTrh/rjCrnA2BfOVneB8fy1UloRRN4CN43Yu1znuyoqKioqL2gs9/5TszOnbQGB2Hy5GQW/9AQjbPPhpkzGbniCrbecAPPrlsXAXk89Di8fT6c26tQi8H0JWEQUaasjzKQ1WndQrIPXOuGX+jxV51jqD+pKlfZpT//PPzyP5LX/9/emQfbddR3/vOTrM2WbEmWLMmrZFsYDCFAHJuEzeyYSYAKTCBhCGEGGCAZwjAkAyELTKqmAhUyITUYhy1ABjAzrI6pCWO2ckiFxTY2XiTZki3ZsiVLsixZ23tP773f/HHOsY76dfc557573yJ9P1Wn7jndv+7+nbvU/d7f/XW3Uejds1bDxZfC2ovh9OWlaMbjGjgbSa7jSRXrXojhvbvh3k2weQM8+kj8vmL3FKtLCeYmIRwrm8ni2DNtwnFjvhylf6kVW+GeXe7XRKqFEEJMEauvuorXX3VVsv7yBx7ghu99j/HR0YH6IYGcoF+pFpT1RyhSLbrkIzeJ5HCMsE3bchraxOzqAjcl1tv4Wh+j3lcooMM+wvEe789h547i+OF3j9XPXwBLlsDylXDmWcXjGcvg1MUwf74//sJUAeeqr0oAjwzDoQOwby/s3QWP7IK9e4o86eozmhO9TfUpgdhUnivrIozD65xo7SqO29qm2qT8HKX4dyYmsivBHPat1AohhJjdzJk7F4AxCeTpox+pFnNq9U35yGFZRUrE9iKGc0K7SSSH1zG/Yv02laVEr0dsuwjmev3wcHHs2QNsSNvnaBLAXe2boqZtRHDb9jlx3IsAztXFyicrjmPi90jQxoP6jqkV71BqhRBCzHwqgTzoCHKHKfwnJw/Ab++AwzlhEn4x123Ha2VOIZJjES9IiwQi9rn6tuWhr7k2OYETRurCfmN95PrNtR+PtIvZxYRZW/vU0bWfnF3sPtrcf5uy3OuRE7lhfe4175c4jo2f6rPe3+FaWfgZC8Vx1b4hteLvItVCCCFmGBddcAHvu/hizl6xYqDjSCA34O7jW+BZR2J1TPwSH2Wi+KmLlHGO/1uYoP1UieTJCOjYdSiGehF6YX+h8AntxxNtQ7qI4F6PprG73EvM71ifYVlMCKdev9Tz3UYcp3zpVRx3ee9X4jh8XusotUIIIU5cli1dyq8uX86SRYsGOo4Ecgsec7/1Hrg29iUb+zI+ykTxVj+OUgjpehtIi4rwsR8iuU2bVF2T6MkJr1Q/ufK2IjIUoTnhOghCP1LR7tz9kChLleeEcZvXN6xven1j44Ztcq9Tqs/ce75+PszxPzbDPseY+NlKPZ9KrRBCiNnHwaEhbt63jwMHDgx0HAnklmyH17dNtQAYIS2SAYY4Nvu+i0AK61MCNSci2wjosN+u0eSYfe5eu5bnhHBMNMfEc1shnWqT6jvXvkkUdylveo5zr2msjybh3FQO+dUqUm2b3ifV+dHyyL0eYWpFVR7iKLVCCCFmI/fv3MkH7r6bLfffP9BxJJBb0pRqURcG9S/0MDeyXn84UpYTErn6mMjoIoZpaJMbJ3cdi/LFIsqx9mF5GyHcJHrbit1eRXUXn3J95cpjr1Euz5jEde790DXdIvwMNI3T1HdYNkbxo7JeFj7XIw0+1vtUaoUQQsxO5pxSrC8xdvRog+Ukxxlo7ycYuVSLmBCovthDAVoXH02T9iYjknspT9U1ieguQjks60U4tqlvErm90CSu++V7rLxJGEO6bapNV+HcVTT3Io5jn6PDQVnbSXmx51mpFUIIMXuZO28eAD+74YaBjiOB3JHt8PoH4NHYF29MeIxxvEiu7Opf8EMcLwYgLyxy9ZMVyTFxEraLtW3qr7rOpQV0EY2p9l3sclHjLtHkfvkVq4v9w5B7HlPXYVkX4dxUN2hxfCjwI/RlLLCP3V/FOLAFblFqhRBCzE7O++VfBuDwwYMDHUfrIHfE3cfnm110KuxaAaeE6+iOM3F95NHycW5pU9VVttWGBwtq/XTdSITIecoGJgqU2C56dZs2G4YQsUld1+8jtkFI3T7sN1aXsqvTZc3jtqTG6tU+JiDblOf6z/3QCev7UddPcXykdh6uWFGJ47FI+9S/PNth/2a4PFIthBBiFrBw6VLWLF2Keddv4G5IIPfAiPujS82euRBuWsxE4VUJ34pqExFjYsi+enmrv4inUiTn6mL9xfqszlNCuY1wjgn0sH3TBiExm5DBfpR6HzMlQKuytiK4qW1O4DbVt4kax+omI47rE1lTueyxFStS4vgRGLsbLnL3sYiJEEKIWcInnvc8WL9+oGMoxaJH9rnfvBF+L5UiHhMJbZZ/GwnKUn21rY9dN03e61IXXjf99Z/rq+s6wSGx53Qq6TJ+7r5yr1OTyK3K2kzaS9XH2qfaxdqE572K49FEX1VZ20l5AAeBjfDsIfdHEiZCCCFmC6ecAmODjXUogjwJHna/+gKzF1wKr54b1FVf2mEENrUddUX1pb+g1q7LttL1sjpdUi5ibep+VH3kUihSfdTLcm3r9xCrC/tJRY+nWiTHiP3giNnk6tqWxYRwzqZJOLepS4nt1DhdxHHs6DopbwTYBP95n/uPItVCCCFmGV/YvJnF+/fzygGOoQjyJNnm/pr74L62m4g4EycVhRGyEfLRsZzISQmVJptUXZPIaRJKqTFSZbH+UusMh8/hdG4S0tWnmE2qn1h5k23sOhdVjtWnxp1KcRw+f5U4Du859ZqPAXfDN3e6/02kWgghxCzkpj17uO6nPx3oGBLIfWATXLoThmNf0G1WtgiFABST9qZKJDetLNFVBMdsehGBqfJeNgmJ3VNXER1rn1vlIudXL/cdK0+tcpFr2yR+2wrn6RDHo4k+UnnH98H9W91fFakWQggxSzl85Ai79u/nse3bBzaGBHIfcPehu+GJ+yA6pzImnEYpRHJ9c4Xx4HG4PGLRMoIyEudtxFOsLtZHG6Ectk+JrSbbsLyNEM4J35RAHvQSb22iyE0+x8raLvPWq/gN7cPzVDQ69x7sIo5jr9Moxz4v4ecmxIGdMLwJnhipFkIIMYu56lVF3GPk0KEGy96RQO4TB923boRXDCXqY4IlJ5KrY4RmkZyL8MX6zU3S6xJNhok+5SLFqT5jtm3q2ojmfkWPcz70EkXOCeJU21S0ODcOEZvJCOfU650T07H7ahLH4fu1qzjeD74JLnX32OaXQgghZjGnnXEGAD4e+xboDxLIfeQR9+s3wodHI3UpoVCJ5JzoOEpaJPfjvKkO4v7F7q+NTThOboWKJkHbRvimBHKX6PFkosht7VJ1XVazaLt6RUz8hm3avkdCm5R9zIcjTEyrCJ/n8EckNdsYQ8BG+I2D7vcmTIQQQsxi1q5ezavXrGHRwoUDG0OrWPSZB93/60Vmz74EfjW15nG4ikQlqOdS/GKpvvzr7auJSfW3QpuVKZrWS67KY3axulzbnA0Zu/p1yr5uQ6I+ZdemTS/kxum1TUq4NpWlRGq9LGeTq28Syk1tUuL4MPml4ZxjPyDDtilxPApsgv+xx/0bCRMhhBCznPXnnsv6886D004b2BgSyAPgXnjOIth5AayMibxK/MZEMqRF8iiFqDi1VtarSK6L25RYzW0uUq9vI5RpsIuNnxPLYduwjxS9iNpeaTtWV1FclcdEbpNdqk2vEeXUeU4cj1O8j3PCeJyJu+Q1ieMxYDP8eLv7uxMmQgghTgBGx8cZGh3l1LGxgaVCKMViALj7+J1w/nZ4LCVa6n8ZVzb1nOSUeBgHDtXKU0vJtRUv9fMuuclt2qd86bKKQ/2+c6kYsT7CYxD0Ml7TPaX6CN8TYXloO5l0i1z71Psr9n4Ny8co3r9N4rhr5Hgc2AIbtsCvJEyEEEKcINz485/zW7fcwq7duwc2hiLIA8Ldh8xsjcGuc+C0rpHk+sYj4xyLvFa/aA4BiyhewHr7ULjkosS9RJPrtvRgE/oZsw3tY/6EfobnMQYlktuM2SSYu7ZPtcn9OGpr00YMx85TQrs6qlz6qjwnjsMfkLkfHaU43nIPPNk9upCMEEKIEwkrvvEHOUlPAnmAuPvhU8zWzIFda2DhZERyRWVvFBOc5pdHRdeUCyJ2MZE7GaFMxg6OFz9t8o9TwjsmmmNt+k1b8dq1PiVu29hPVhiH17l7bCOOq5UqqvIqel6367paRWWzFR64Gy6ROBZCiJMDs0F+qxdIIA+YUfcD883OngM7V8H8lEgOI6GVUDiF4/Ng6iLQKJaBG6fd5D3Ii+ZUNDm0JWEf8zEVCYZ8/ylxHbZr6rdtm7Z0VWBt7buK4qq8jTifjDAO7bvYVceh4LoueOvXR2t9tBXH22DnBrjQ3ccSZkIIIU5QtMzbLGfE/dE74bzdcDQnYMK/kivREMtDrl+PMlGENPVNwi513hRdTIm10CaWp5zqI3W/MWL3mbNvatOvvtrYT1cecjhG6rpuHxsvPK/nxx+M+Bx7P48E1wTnIePAA/DIXXCBu8dWVhSTwMzOM7Pvm9kGM7vTzP5gun0SQogKm1PIV6VYnAAccd+1yOz8p8C2sxKR5PC6shmhSKOoUisq4RDmJR+kWOFibq39IKPJOb+b8o/r99cmnSIm1NrmH+eE6lSQEpkpm1x9275ydqm+cmK4rV11PUqRAlSVpybvjVH8CAzbN4nj+2HPnXCuu48kzMTkGAX+i7vfYmZLgJvN7AZ3v2u6HRNCiHVr1vD6c87h9CVLBjaGBPIUcsR950Kzs58MD66GBTFh6BwvhOFYKsW8Wl2del7yYYp0i/oLGxOZsbo2aRThNQ1taLALx247WS+si4nm2PUgSQnRtvaxuphNrqxJQKeiyKk2uahxrF+nmIhXF72xfONKHI8G5fU2Mcq0ih1l5PhowkxMEnffAewozw+Y2QbgHEACWQgx7VywZg0XnHPOQNdBNs1rmXoWmJ3xJHjoHDg1JeDqQrguRquc5DlB+ZygzVyKVS5ifcBEQZqLxoZ1uetUm1i/qfHatGlq30/7GF0/NW3su4riqryNMO+XMM7Vhesbx1ImquvYShXQLI7vhW2b4CLlHE8dZrYWuBF4irs/FtS9FXgrwKpVq37p2muv7dT3wYMHWbx4cX8cHTDytf/MFj9Bvg6KXn1ddPPNLL32Wna9+92MrVnTqe3zn//8m939siY7RZCngWH3/aeYrXJ46FxYEhNs9ahwPTo6WpbXJ+8Zx0RFJVjHKFIuFnH8ahhto8mp614jyiTsQl+6TNaLta/T1b4fdOm7TbQ3167XdIu2wji0z0WNwyXc6mI4FL9Hg7K6wE7dR7mU2z13wxPdfXBJZ+I4zGwx8FXgXaE4BnD3TwCfALjsssv8yiuv7NT/D37wA7q2mS7ka/+ZLX6CfB0Uvfr6z3fcwXvuv5+rL7yQ8664ov+OIYE8bYy6HzSzVePw0HmwNDZbssozDpeBq0RGlXIR2ldtjCIPtIomw0SRWl3nRGpKKBNcx8RsWzsi9uFYMR/7nXvctBpGr7QVs03t20aLY2P2KoxzfdWjxlVdbrOQkUgfoW1IuUPeHZvhqVrKbeows3kU4vgL7v616fZHCCEqNEnvBMfdj5jZ6hHYcCGsi617XAmLWK5wbPIetcfJRJPD+jbCsSlanIoqNCrn7gAAFnxJREFU18dqE12ul6fqe02j6If6aiOG24yVEsW5si5CN3c92ahxrKzrZDzKNpvg21vdX5YxE33GikVGPw1scPe/nm5/hBAihpZ5O4Fx9+FN7hfeBV+plrqaYEN6ua4RjuVzhpG98aDdkfKI/f0NE/vPRSBj9rH+2i5hFvaRE4CpiGMoyMKjqX0TsfZN47T1PbTr9XlrKgvHSvkXq68v03aIibvixe6/yjfuKo69HOM2+GOJ42nhWcAbgBeY2a3l8fLpdkoIIYDHNwoZ5J+KiiDPELa5/9tVZm9/Ily9mHgENJeXHNtUhMCuSzQ5lV8cRnzDCG69Tb2MwC70L+yjqmuKCqc+Gm1tB/HR6tJnl+huWNcU/U31kfphkfOpsmkTNaZWlso3bhLHe2F8Izz/UfcbE2ZigLj7D5naRWCEEKI1VYrFIJFAnkE87P7xpWY/uQR+tAJO6SKSq/zOcCm4qq6X3OSqfRuhHGsT9hsTvk0pFTGx1iaVootIbVIB/RLRTcK1aayYAK6u24jrLsI4rO+aa1xFjmPjNE3GexAO3A3rj7g/nDATQghxErNuzRrefP75LDv99IGNIYE8w9jnfvMCs5VPgHvPg2W5yXthjrBRROzCVS6o2VGzraLJCyiEddhX2CYnlMmU5aLKbcVyqv8uojnFoKPIbURvri4ndrtGi6uytsLYmbiu8WSixjExXzEGbIE774Gna41jIYQQKdasWMErV68GbRRycjHsvs/MVhyBf70YLo+9SJXQ6Ec0ebg86huMdBXKoU2sLCWKQ7HXVjDXbZqEYkw09/L/8WTSNNrYdRXFXXzqKoxHOH7VCegWNQ4fcykVI8BG+PQD7m9OmAkhhBAAHBkZYc+RI6waHmb+gMbQJL0ZiruPb3a/4g74SDWxLka4S1n97+uj5RGboBabxDdEMTFqNLALRVXqb/aYTVgWRhNz7VMT3lL95mxjNrEJdk1Hr2O18T98LkI/m6KxubLchL1w7KMU/yyES7Kl/KiixmFKRdgudf8HwG+DN0gcCyGEaMPPt2zhHbffzgM7dgxsDEWQZzjb3d+z0uzGS+AbZ4ClJp+NMblocj0afKQ8X8TEX1BNayLHbOplsfJcBDkW8exnDnKXKHLbPtu0zwnmVH1Tm7CsKdIcCt0xitc+rI8J/clGjav63XB0I/ziY+4bMqZCCCHERLSKxcnNbvfrFptdeAlsXA0LUqIu3FjEaZ+bHBPKh0vbrkK53h+BHS3Kw35S9eH4Mf+aGMRHq0mYpur6IYq79FmP8B6JlMdEcXU+mVzjqv022LkJ1o+6H8yYCiGEEFG0zJvgoPtWM1v6BLhtHTxhLumlzNpEkyshHbYNhXK15u0pFJP55kRsq+uK1MQ8Erax8no/YX2TaI5dh+261LUVuE10iSB3qesijOvllTCObewRCt7qsR9R42FgC/zTVni5dsYTQgjRFe2kJ47D3YeAS84ye9t6uHppIuUCmqPJxrFocmwyXH0iHxSCaJRjQjlcQQMmitU2q1WkxHJYF+s/HCtmm+uvbV0XukSB24zbNRUjdR0K2CG6C2Pn2GoWRGzbRI13w9EtcNVe9+9mTIUQQogkWgdZRNnlfs0Cs2vXwU/XwsVdo8kVlVCuIso5gRoKZaMQyrE3UC6qXJV3FcthfW6sJtsmYs9Vil4jyL3YtBXFVVn4Q6K+yUe9XWziXfWYEsYxu5zfZdT4O1uLqLGWcBNCCNEz684+m3euW8dZy5cPbAwJ5FnKsPs+YH0v0eSKSgiOkBbKofCp/2YbKh/nAfNJp19UZeG49fJUFLpO16hwm2XiYkwmmty1bS+CuKksjOoOcywlol5eF8ixdpUwjvUb6yeGosZCCCH6zYply3jxypWwePHAxpBAnuVMJppcF7KVUE5N5AuFcj1PuVpObg7FWsqhUIbmSHFTdDlsF7aPMdOSW7tEkLukasSE6yiFMM6J2lxdThhX501R4xFgs6LGQggh+syhoSG2HzzI+UNDj+8I3G+0DvIJwLD7vo3u62+Bd+wjv+hJuPaxB9fVRL426yeH6+JW2xFXa+iG7WLRxlgUMhwrF6kM/WsT1RwUOV+afG9a6zj1vIXP/xDF8z+U6D+3nnO1MkX12oU+VNdj5MXxOLALjt4ML7rP/cUSx0IIIfrJpm3beM9dd7HtwQcHNoYiyCcQu9w/vsDsS+vgprVwUS6aXI8exyLKY+VRRZNjv6Qq0VTV1furduebQ5GrPDfSrimtou2kvMmuhTwVxO6hyS5WHvshMUYhascjdWEaRayvcY693mH7+nXTXOFa1Pi7W+EqCWMhhBCDRKtYiNaUuckXn2X29vXwsVxuchuhXOWuzi2P2GS+6u0ZW0u5Wkqsqp9P8aYLUy6aUjEI6lOCM3avbSfnTYYuKRFd62OiuEpriT0PuTSKsH6UuLCO2eZQrrEQQoipQqtYiJ6pRZN/fB48YT75SGvbiHK1PFxsHWWCNjBRUFeRZTg2uS8llmGizznRnKpvatOVXiLTbdqkRH8VvR3h+Ml2dbsu0eJKGMeEcHjdRsQfBrbBd+5TrrEQQogpQOsgi0lRRpMvOcPs8vPh+rNhZRi9rdNGKFfRS2iXfpESy/XJfVVf9ZUwYv0QOa/b5Oh1RYs2TFYwp0T9KBNTJ+r2sXzlWN/1VIx6GkXYLrRv8n8EeADuvh9ecth9W0MTIYQQoi+cftppXL50KUtOPXVgY0ggnwTsd/8JcNaZZr9+PnxpFZyWyk+GYwKpnltcF8pVWRXNrFa+gN7FcrW+cnU9t+wz5meXPOSwzXTQRgxXArb+PIRtYo+pFAvIR4tjfrUVxkeBh+DhbfDSA+63NTQRQggh+sqFV1zBn37oQ/C0pw1sDAnkk4hH3P8RWHyW2dsugL9dCfNyWTyx3OJQQBnH/v6vhG0qqgwTxW0qFSMUim1Fc90+Rj/zkXvNO86J4bCfVJQ4NlZd5Na3gw7rU22a/B4DdsKBbfDqR91vaGgihBBCDIaVK+FVrxroEBLIJyG73K8xs79bDX+xFt67DObmhHIq9YKgvIpWQiGQ50baWNAuFV2uP1bkRHO1vnPYf+xeBkksyluJ1SYxnGqfy02u21QiNpeDXD9vk7lV+b8HRrbCW3e7f65FMyGEEGJWI4F8kuLuDvyJmX3wXPjEWnjjErA2QhnSecpVXbXebnVdRZZjq2DU+4e4QG4rmutUY1ViPRy/l2hyKHxjawu3aR/2VS+LRXljbSofxjJtekmjqBgH9sLYVvjzh+G/l+8ZIYQQ4oRHAvkkp1x14E1m9vtr4asXwEtPo11Ob1NUmVp5mF9cCfHcShh1NdZFNNfHTaUbDJJc6kWTGI61DyO+scl2qTHq43QRxo+Bb4VrHoR3uvtUPn1CCCHEtCOBLABw90PAyxaYnXUuXHcOXL64XEM5J5brwiuc1FfVp8QyHFtfuarPRZfbiGYi9TmblH193By5aG9Y1jZnuU3qRFPfXUVx9UNiP4w/CN+4H/6dux9paieEEEKciEggi+MYdt8FPNPMFq6BvzgH/tOZsCC36kVFOKmvOk+JZTg+GlpFly3SR4yYaA7pKpK70EYcN9WnJuOFW4LH7Nv01YRTrEqxCx57CP50N3zM3WMBaiGEEOKkQQJZRHH3IeAPgT9cavaiNfDJ1bB2EekVKh5vy8RIbywNo6qDeHS5Ku8immO+1CfuDYo2fXcVw7F2sbougrje5iD4Q3DLTvjdg+53dOhCCCGEOKGRQBaN7HP/DrBurtmyc+Gas+HVS8uVL9oI1dTkPjhe/MWivW1EM5HHcPx+R4/D/nPlXcRw2F8ugtxV8FepGo/A8IPwsR3w/vKHkBBCCCFqSCCL1oy5Pwq81szsTHjj2fCRVbA8t411SEwsV+cxYdhWNNfr66tVhH5NZgWLul+x81AIhzZd+ovZ9CKKKX05AuyEbTvgLfu0hrEQQgiRRQJZdKZc7uuzwGcXmV14NnzubHjWkhaT+o7rh3gqBsTTInKiuV7flECby0uO9Re77kWotm0/WUFctR0D9sPYg/D1B+E/jrnv7bE7IYQQ4qRCAllMiiPu9wLPMbN5Z8IbzoQ/WgnrF8OcNhP76sQEYVvRTGCTK5us0A1pivy2aTsZMVzv6yhwAMZ2wy274b89Bt/S+sVCCCFENySQRV8o11P+THmwxOw5K+DPV8Jzl8K8U8gv45bsl3aiuW4f6yNkMjnJk4nqVo/9EuZOsc33XjjyCFy/Gz5w2P2uPnQvhBBCnLRIIIuBcMD9n4EXASwwu2Al/MkK+M3lcPpCehPLFSmB2TZ1ouqjX/Q7FaNprHHgMLAXdu+BT+6Cjyh9QgghhOgfEshi4Ay7bwPeArzFzBafBW87E965As49FaxrKkaKtsK032MNknpe9QEY3wMbd8NfPgpfdveRKXJDCCGEOKmQQBZTirsfBP4K+Cszm3M6vHw5vH0JXHE6LFsEc+YxuQhzow8D6neyVH6NU+QSH4KxA7DrMbhxL3z0oPu/TqN7QgghxEmDBLKYNtx9HLi+PAAwszOXw8tPh9dNpWieanJieD986TH4nrsfmEYXhRBCiJMWCWQxo3D3R4B/KA8gLZrDiX8zUTyHE/MiYviLpRg+OH1eCiGEEKKOBLKY8aRE81J43kJ42ny4dCGsnwerFsLp82HBPLD5YHM4fuMQ6G8OciV8xynyhI+Cj8D4CAwNw/5h2DEMdw/D7UPws8fghxLDQgghxMxGAlnMSkrR/LXyiGJmi+fCuYvhSQvgFxbAUxbC+vmweh6cajDHis1N5gBz5hSCek41xHh5eLlD9Hi5/8YIHBiBh4bhnmH4+RDceQDuctihrZvFiYKZvQz4KMXO7p9y97+cZpeEEGLKkEAWJyxlpHZjeXx9mt0RYtZgZnOBjwEvBrYDPzWz61xrbAshThLmNJsIIYQ4ybgc2Ozu95bLCV4LvHKafRJCiClDEWQhhBAh5wAP1K63A1eERmb2VuCt5eVBM9vUcZwVwJ6ePJx65Gv/mS1+gnwdFNPh6wVtjCSQhRBChMTmsk5YQtzdPwF8oudBzG5y98t6bT+VyNf+M1v8BPk6KGayr0qxEEIIEbIdOK92fS7w0DT5IoQQU44EshBCiJCfAuvNbJ2ZzQdeB1w3zT4JIcSUoRQLIYQQx+Huo2b2+8C3KZZ5+4y73zmAoXpOz5gG5Gv/mS1+gnwdFDPWV3OfkFYmhBBCCCHESYtSLIQQQgghhKghgSyEEEIIIUQNCWQhhBBTjpm9zMw2mdlmM3vvNPnwGTPbZWZ31MqWm9kNZnZP+bisLDcz+9vS35+b2TNqbd5Y2t9jZm8cgJ/nmdn3zWyDmd1pZn8wg31daGY/MbPbSl8/WJavM7Mfl+N+uZz8iZktKK83l/Vra329ryzfZGYv7bev5RhzzexnZnb9DPdzq5ndbma3mtlNZdmMe/3LMZaa2VfMbGP5nv2VmeprFnfXoUOHDh06puygmPi3BbgQmA/cBlw6DX48F3gGcEet7MPAe8vz9wIfKs9fDvxfijWinwn8uCxfDtxbPi4rz5f12c81wDPK8yXA3cClM9RXAxaX5/OAH5c+/G/gdWX5NcDby/N3ANeU568DvlyeX1q+LxYA68r3y9wBvAfeDXwRuL68nql+bgVWBGUz7vUvx/kc8ObyfD6wdKb6mjsUQRZCCDHVzIitrN39RmBvUPxKii94ysdX1co/7wU/Apaa2RrgpcAN7r7X3R8FbgBe1mc/d7j7LeX5AWADxW6HM9FXd/eD5eW88nDgBcBXEr5W9/AV4IVmZmX5te4+7O73AZsp3jd9w8zOBf4N8Kny2mainxlm3OtvZqdT/PD8NIC7j7j7vpnoaxMSyEIIIaaa2FbW50yTLyGr3H0HFMIUOKssT/k8pfdS/rX/dIrI7Iz0tUxbuBXYRSFstgD73H00Mu7jPpX1+4Ezp8jXvwH+CBgvr8+coX5C8SPj/5nZzVZs8Q4z8/W/ENgN/H2ZuvIpMztthvqaRQJZCCHEVNNqK+sZRsrnKbsXM1sMfBV4l7s/ljNN+DQlvrr7mLs/jWIHxsuBJ2XGnRZfzezXgF3ufnO9ODPmdL/+z3L3ZwBXAb9nZs/N2E6nr6dQpC193N2fDhyiSKlIMd3PaxIJZCGEEFPNTN7K+uHyL17Kx11lecrnKbkXM5tHIY6/4O5fm8m+VpR/rf+AIrd0qZlVm5PVx33cp7L+DIq0l0H7+izgFWa2lSLF5wUUEeWZ5icA7v5Q+bgL+DrFD4+Z+PpvB7a7+4/L669QCOaZ6GsWCWQhhBBTzUzeyvo6oJox/0bgm7Xy3yln3T8T2F/+Vfxt4CVmtqycmf+SsqxvlLmunwY2uPtfz3BfV5rZ0vJ8EfAiipzp7wOvSfha3cNrgO95MUvrOuB15eoR64D1wE/65ae7v8/dz3X3tRTvv++5++tnmp8AZnaamS2pziletzuYga+/u+8EHjCzS8qiFwJ3zURfG5nKGYE6dOjQoUOH++Oz1++myE99/zT58CVgB3CUImL1HyjySr8L3FM+Li9tDfhY6e/twGW1fv49xeSszcCbBuDnsyn+Xv45cGt5vHyG+vpU4Gelr3cAf1aWX0ghHDcD/wdYUJYvLK83l/UX1vp6f3kPm4CrBvg+uJJjq1jMOD9Ln24rjzurz8tMfP3LMZ4G3FS+B75BsQrFjPQ1d2iraSGEEEIIIWooxUIIIYQQQogaEshCCCGEEELUkEAWQgghhBCihgSyEEIIIYQQNSSQhRBCCCGEqCGBLIQQQgghRA0JZCGEEEIMDDO70syuL89fYWa5rYfb9rnVzG43s8tq1yt67Ou1Zra58lEIKPbMFkIIIYRoTbm7n7n7eJd27n4d/ds18fnuvmeynbj7l83sYeA9ffBJnCAogiyEEEKIRsxsrZltMLOrgVuA88zs42Z2k5ndaWYfrNm+zMw2mtkPgd+olf+umf3P8vyzZvaaWt3B8nGNmd1oZrea2R1m9pwOPi4ys38ys7eU/m40s0+V/XzBzF5kZv9iZveY2eV9eFrECYoEshBCCCHacgnweXd/urtvo9j2+DKKLaafZ2ZPNbOFwCeBXweeA6zuOMZvA99296cBv0ixtXYbFgP/CHzR3T9Zll0MfLT074ll38+miBb/cUe/xEmEBLIQQggh2rLN3X9Uu/5NM7sF+BnwZOBSCiF6n7vf4+4O/K+OY/wUeJOZfQD4BXc/0LLdN4G/d/fP18ruc/fby1SQO4Hvlj7dDqzt6Jc4iZBAFkIIIURbDlUnZraOIhL7Qnd/KvAtYGFZ7S36GqXUIWVO83wAd78ReC7wIPAPZvY7LX37F+Cqsq+K4dr5eO16HM3DEhkkkIUQQgjRC6dTCOb9ZrYKuKos3wisM7OLyuvfSrTfCvxSef5KYB6AmV0A7CrTJD4NPKOlP38GPAJc3eEehIgigSyEEEKIzrj7bRSpFXcCn6GI4OLuQ8BbgW+Vk/S2Jbr4JEXe8k+AKzgWnb4SuNXMfga8miKHuC3vAhaa2Ye73Y0Qx2NFKo4QQgghxOzAzLYCl/VjmbeyvyuB97j7r/WjPzH7UQRZCCGEELON3cB3q41CJoOZvZYiLePRSXslThgUQRZCCCGEEKKGIshCCCGEEELUkEAWQgghhBCihgSyEEIIIYQQNSSQhRBCCCGEqPH/Aca6NSVKmiKEAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, [ax0, ax1] = plt.subplots(1,2, figsize=(10,5))\n", "patches = []\n", "colors = []\n", "\n", "for r1 in (np.linspace(6371, 3, 1000) * u.km):\n", " c = mpl.patches.Circle((0,0), r1.value)\n", " patches.append(c)\n", " colors.append(prem.density(r1).value)\n", " \n", "p = mpl.collections.PatchCollection(patches)\n", "p.set_array(np.asarray(colors))\n", "p.set_clim(1, 13)\n", "p.set_cmap('gist_heat')\n", "ax0.add_collection(p)\n", "ax0.set(aspect='equal',\n", " xlim=(-6371., 6371.),\n", " ylim=(-6371., 6371.))\n", "ax0.axis('off')\n", "ax0.text(0,0, 'inner\\ncore', horizontalalignment='center', verticalalignment='center')\n", "ax0.text(0,1800, 'outer core', horizontalalignment='center')\n", "ax0.text(0,5000, 'mantle', horizontalalignment='center', color='w')\n", "\n", "ax1.plot(r, rho, 'r-', label='PREM (1981)', alpha=0.7)\n", "ax1.plot(r, rho_c, 'k--', label='Const. approx.', alpha=0.7)\n", "ax1.set(xlabel='radius [km]',\n", " ylabel=r'$\\rho$ [g cm$^{-3}$]',\n", " ylim=(0,14))\n", "ax1.legend()\n", "ax1.grid()\n", "\n", "fig.tight_layout()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Electron Fraction vs. Radial Distance\n", "\n", "Plot electron fraction $Y_{e}$ as a function of radial distance from the core." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "y_e = prem.y_e(r)\n", "y_e_c = earth.y_e(r)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAAEYCAYAAABBfQDEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3XuYJlV57/3vjxkOHkBABkVAByOKeAKZIO5Eg25BcOvAjm4FkUMSQzSySbYahSSiosmOvnnF+AYTMaJgUFASZUQMHtHoFpxGRmBAZETcjJAwCCIqAoP3+0etninaPs0w/fTTPd/PdT1XV61atequnp7Vd69aVZWqQpIkSVJni9kOQJIkSRomJsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyNMOS/O8kfzoEcbwnyWtmOw5JkoadCbLmjSQ3Jrk7yc+S/GeSDyd5eNt2SZJftm23JfnXJLv09n1bkvva9tHPT3rbq7W5sFe2MMmtSSZ8mHiSRcAxwAfa+lZJzm+xVpIDx9TfPslZrd1bk7xtzPZ9kvx7kjuTrE5ySm/bpG0D/w/wF0m2mvY3VdJAJXllkpHWB92S5HNJfnsGj3dgktUbuM+NSV4wUzE9GO18KsmbZjsWzW0myJpvXlJVDweeCfwm8Je9bSe0bU8AHg787Zh9z6uqh/c+24/Z/hPg0N76i4A7pojnOOCiqrq7V/Z14FXAf4xT/zTgocBiYH/g6CS/19v+MeBrwI7A7wCvTbJ0Om1X1S3Ad4GlY7dJmn1JXg+8F/hr4FHAY4H3A4fNZlzDpj9QMY5jgdvb19k4vuYJE2TNS1X1I+BzwFPH2fYT4NPAPhvY7EfpRoNHHQOcPcU+hwJf7R373qp6b1V9Hbh/nPovAd5dVb+oqhuBDwG/39u+GDinqu6vqu/TJcRPmWbbAJcA/22KmCUNWJJHAKcCr6uqf62qn1fVfVX1mar6s1Zn6yTvTXJz+7w3ydZt24HtqtIb2tWnW/p/XCd5UZJrktyV5EdJ3pjkYXT95GN6V84e8yDOYYckFyZZk+SOtrxb2/Y/klw+pv4bkny6d25/m+T/tqt1/5jkIWPO7c1J/gP48ATHfyjwMuB1wJ5JlvS2LW4jy8e3790tSd7Q2/62dgXuvPY9+naSZ/S239iOfyXw83YF8cnt6uRPkqwcHaxoV/NWJPmfbX1Bkm/0r/hp+Jkga15KsjvdCO8V42x7JPC7wKoNbPbTwHPbNIjtgecAF0yxz9OA6zbwOBmz3E/y3wsck2TLJE8Cng18cQPavhZ4xpS1JA3as4FtgE9NUucvgAPo/rh/Bt1Vpv5VskcDjwB2Bf4AOD3JDm3bh4A/qqpt6fqUL1fVz+n+iL+5d+Xs5gdxDlvQJa+Poxv9vhv4+7ZtGbBHkif36r+KbuAB4F3AE9u5PaGdQz+hfDTdlbPHAcdPcPyXAj8DPglczAMHNEY9D9gTOBg4acxUkcPavjvSXa37dJIte9uPpBtg2J6ub/4M8HlgZ+B/AuckeVJV3dvO7dR2vicBC4C/miBuDSETZM03n043d/jrdCO3f93b9r4kdwK3ATvRdWh9L28jAaOfr4zZ/ku6DvEVwBF0Hf4vp4hne+CuDYj/3+g67W2TPIFu9Pihve0X0o2Q3E03XeJDVbV8A9q/q8Ukabg8EritqtZOUuco4NSqurWq1gBvB47ubb+vbb+vqi6iSxaf1Nu2d5LtquqOqvr2pj6BqvpxVf1LuwJ2F11C+Dtt2z3AeXSJI0meQndF7MIkAf4Q+F9VdXvb96/p+tlRvwLeWlX3jJmy1ncs3VS5++kS3CPHJLgAb2+j81fRJfNH9rZdXlXnV9V9wHvo/mA5oLf9fVV1Uzv+AXRT9f6mXb37Ml3/fGQ736uBd9L9wfNG4OgWl+YIE2TNN4dX1fZV9biq+uMxHemJVfUI4OnADsBuY/b9RNt39PO8cdo/m25UYjrTK6Cbo7ztBsR/Il3yez3d6PTHgdUASXakS6BPpeu4dwdemOSPN6D9benmUksaLj8GdppifutjgB/21n/Yyta1MSbB/gVdEgfd6OqLgB8m+WqSZ2+CmB8gyUOTfCDJD5P8lO5+ie2TLGhVzgJe2RLio+n63HuARXQDAZePDlDQ9XWLes2vqaoJByTaVcPnAee0ogvo+smxU8pu6i2P/f6t21ZVv6Lre8fd3spvavX67e3aWz+L7o+Ai6rq+oli13AyQdZmp40cvJPu8mOmqj/GvwO70N1A8/Vp1L+S7rLhdGO7vaqOqqpHV9VT6P6Pfqttfjxwf1WdXVVrq2o1cC7dL73pejLwnQ2oL2kwvkl3RerwSercTDfFYNRjW9mUqmp5VR1GNx3g08AnRjdteKgTegPdiPWzqmo74LmtPC2GS4F76aanvZL10ytuoxsYeEpvgOIR7abqdacwxbGPpusvP9PmKd9AlyCPnWaxe2957Pdv3bYkW9ANovS392O4Gdi91eu396Pe+vvpRpVfmBl8EolmhgmyNldn0f2i2KAnOlRV0d1It7QtT+Ui2iXGUe1mlG3a6lZJthlN1JP8RpJHtps6DqWba/fOVvd7XZW8MskWSR5NN93jO9Npu/kduptyJA2RqrqTbs7t6UkOb6OxWyY5NMm7W7WPA3+ZZFGSnVr9f56q7XbT2FFJHtGmD/yU9Tfy/ifwyHQ3CW6ILVv/MvpZSHeF6m7gJ+2K11vH2e9sunnJa9sNxaOjtR8ETkuyc4t51yQv3IB4jqGbcrJP7/NS4L+1+05GvaV9b58C/B7dtI9R+yX53XYufwrcA1w6wfEuA34OvKn9Ox1I97vh3Bb/0cB+dE8yOhE4K+2xo5obTJC1WWo3UbwPeEuv+BV54HOQfzbaWY/Zd2VVrZzmoc4GXpR2N3ZzHd0vkV3pbiS5m/WjQvsBV9HNFf7fwFGjx6qqn9LdXPi/6KZurACu5oE3fkzYdrrnPu9NN3okachU1XuA19PdeLeG7pL+Caz/P/tOYITuytRVwLdZ/wf0VI4GbmxTH15DmwtcVd+lS7xvaNMbHtOS6an6uIvo+pfRz9vobiJ+CN2I8KV00yTG+ijdTYIfHVP+Zrobpy9tMX6R9fOnJ5XkALqpDKdX1X/0Pstam/15xl9tZV8C/raqPt/bdgHdoMMddN+v321/UPya9jtkKd1NjrfRjRYfU1XfTfLY9r04pqp+VlUfo/t3O63FO53vr2ZZpjcIJmljJflr4Naqeu8sx/H/At+vqvfPZhySNl9tsOBW4JmDnJebZDHwA2DL8W6ETPdSpidU1asGFZOGmw+7lmZYVf35bMcAUFVvmLqWJM2o1wLLvWlNw26gUyySHJLkuiSrkpw0QZ2Xp3uY+cokH+uVH5vk+vY5tle+X5KrWpvv24ibriRJ0gxLciPwJ3Q380lDbWBTLNpjXr4HHET36JTlwJFVdU2vzp50d9Y+v6ruSLJzVd3aJvuPAEvo7iK9HNiv1fkW3X+4S+nmRL2vqrwJSZIkSRtlkCPI+wOrquqGNrn9XH79/fJ/SDfJ/g6Aqrq1lb8Q+EJ7BNYdwBeAQ9pNR9tV1TfbEwXOZvJH5EiSJEmTGuQc5F154EO2VwPPGlPniQBJvkH3Wsa3VdW/TbDvru2zepzyX5PkeNrrKR/2sIftt9dee230iUjSMLr88stvq6pFU9f8dfaRkjYH0+0nB5kgjzc3eOz8joV070g/kO4B3f+e5KmT7DudNrvCqjOAMwCWLFlSIyMj04takuaIJD+cutb47CMlbQ6m208OcorFah74Bpuxb6gZrXNBe4/8D+ie6brnJPuu5oGvCx6vTUmSJGnaBpkgLwf2TLJHkq2AI4BlY+p8mu5d6rS3BD2R7nWRFwMHJ9khyQ7AwcDFVXULcFeSA9rTK46he9C3JEmStFEGNsWiqtYmOYEu2V0AnFlVK5OcCoy0N96MJsLX0L0G88+q6scASd5Bl2QDnFpVt7fl1wIfoXt7z+fwNbqSJEl6EAb6opCquojuUWz9slN6y0X3ms3Xj7PvmcCZ45SP0L22UpIkSXrQBvqiEEmSJGnYmSBLkiRJPSbIkiRJUo8JsiRJktRjgixJkiT1mCBLkiRJPSbIkiRJUo8JsiRJktRjgixJkiT1mCBLkiRJPSbIkiRJUo8JsiRJktRjgixJkiT1mCBLkiRJPSbIkiRJUo8JsiRJktRjgixJkiT1mCBLkiRJPSbIkiRJUs9AE+QkhyS5LsmqJCeNs/24JGuSrGifV7fy5/XKViT5ZZLD27aPJPlBb9s+gzwnSZIkzS8LB3WgJAuA04GDgNXA8iTLquqaMVXPq6oT+gVV9RVgn9bOjsAq4PO9Kn9WVefPWPCSJEnabAxyBHl/YFVV3VBV9wLnAodtRDsvAz5XVb/YpNFJkiRJDDZB3hW4qbe+upWN9dIkVyY5P8nu42w/Avj4mLK/avuclmTr8Q6e5PgkI0lG1qxZs1EnIEnzlX2kJK03yAQ545TVmPXPAIur6unAF4GzHtBAsgvwNODiXvHJwF7AbwI7Am8e7+BVdUZVLamqJYsWLdq4M5Ckeco+UpLWG2SCvBrojwjvBtzcr1BVP66qe9rqB4H9xrTxcuBTVXVfb59bqnMP8GG6qRySJEnSRhlkgrwc2DPJHkm2opsqsaxfoY0Qj1oKXDumjSMZM71idJ8kAQ4Hrt7EcUuSJGkzMrCnWFTV2iQn0E2PWACcWVUrk5wKjFTVMuDEJEuBtcDtwHGj+ydZTDcC/dUxTZ+TZBHdFI4VwGtm+FQkSZI0jw0sQQaoqouAi8aUndJbPpluTvF4+97IODf1VdXzN22UkiRJ2pz5Jj1JkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6hlogpzkkCTXJVmV5KRxth+XZE2SFe3z6t62+3vly3rleyS5LMn1Sc5LstWgzkeSJEnzz8AS5CQLgNOBQ4G9gSOT7D1O1fOqap/2+ade+d298qW98ncBp1XVnsAdwB/M1DlIkiRp/hvkCPL+wKqquqGq7gXOBQ57MA0mCfB84PxWdBZw+IOKUpIkSZu1QSbIuwI39dZXt7KxXprkyiTnJ9m9V75NkpEklyYZTYIfCfykqtZO0SZJjm/7j6xZs+ZBnookzS/2kZK03iAT5IxTVmPWPwMsrqqnA1+kGxEe9diqWgK8Enhvkt+YZptdYdUZVbWkqpYsWrRow6OXpHnMPlKS1htkgrwa6I8I7wbc3K9QVT+uqnva6geB/Xrbbm5fbwAuAfYFbgO2T7JwojYlSZKkDTHIBHk5sGd76sRWwBHAsn6FJLv0VpcC17byHZJs3ZZ3An4LuKaqCvgK8LK2z7HABTN6FpIkSZrXFk5dZdOoqrVJTgAuBhYAZ1bVyiSnAiNVtQw4MclSYC1wO3Bc2/3JwAeS/Iouqf+bqrqmbXszcG6SdwJXAB8a1DlJkiRp/hlYggxQVRcBF40pO6W3fDJw8jj7/R/gaRO0eQPdEzIkSZKkB8036UmSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSz0AT5CSHJLkuyaokJ42z/bgka5KsaJ9Xt/J9knwzycokVyZ5RW+fjyT5QW+ffQZ5TpIkSZpfFg7qQEkWAKcDBwGrgeVJllXVNWOqnldVJ4wp+wVwTFVdn+QxwOVJLq6qn7Ttf1ZV58/oCUiSJGmzMMgR5P2BVVV1Q1XdC5wLHDadHavqe1V1fVu+GbgVWDRjkUqSJGmzNcgEeVfgpt766lY21kvbNIrzk+w+dmOS/YGtgO/3iv+q7XNakq03adSSJEnarAwyQc44ZTVm/TPA4qp6OvBF4KwHNJDsAnwU+L2q+lUrPhnYC/hNYEfgzeMePDk+yUiSkTVr1mz8WUjSPGQfKUnrDTJBXg30R4R3A27uV6iqH1fVPW31g8B+o9uSbAd8FvjLqrq0t88t1bkH+DDdVI5fU1VnVNWSqlqyaJGzMySpzz5SktYbZIK8HNgzyR5JtgKOAJb1K7QR4lFLgWtb+VbAp4Czq+qT4+2TJMDhwNUzdgaSJEma9wb2FIuqWpvkBOBiYAFwZlWtTHIqMFJVy4ATkywF1gK3A8e13V8OPBd4ZJLRsuOqagVwTpJFdFM4VgCvGdQ5SZIkaf4ZWIIMUFUXAReNKTult3wy3Zzisfv9M/DPE7T5/E0cpiRJkjZjvklPkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqQeE2RJkiSpxwRZkiRJ6jFBliRJknpMkCVJkqSejUqQk5ye5CNt+eAN2O+QJNclWZXkpHG2H5dkTZIV7fPq3rZjk1zfPsf2yvdLclVr831JsjHnJEmSJMHGjyDfC9zQlp8/nR2SLABOBw4F9gaOTLL3OFXPq6p92uef2r47Am8FngXsD7w1yQ6t/j8AxwN7ts8hG3dKkiRJ0sYnyL8AHpFkS+Cx09xnf2BVVd1QVfcC5wKHTXPfFwJfqKrbq+oO4AvAIUl2Abarqm9WVQFnA4dv0JlIkiRJPVMmyEleNU7xW4Hv040If2yax9oVuKm3vrqVjfXSJFcmOT/J7lPsu2tbnqpNkhyfZCTJyJo1a6YZsiRtHuwjJWm96YwgH53k79oUCQCqam1Vvb+qjq+qC6d5rPHmBteY9c8Ai6vq6cAXgbOm2Hc6bY7GfEZVLamqJYsWLZpmyJK0ebCPlKT1ppMgHwLcDXw5yc4P4lirgd1767sBN/crVNWPq+qetvpBYL8p9l3dlidsU5IkSdoQUybI1TkJ+Dvga+0y3P5JHrqBx1oO7JlkjyRbAUcAy/oV2pziUUuBa9vyxcDBSXZoN+cdDFxcVbcAdyU5oD294hjggg2MS5IkSVpn4XQqJXkx8Gq6p1c8E3gV8JQkd1TVE6bTRlWtTXICXbK7ADizqlYmORUYqaplwIlJlgJrgduB49q+tyd5B12SDXBqVd3ell8LfAR4CPC59pHmrCPPuJSVN98522Foll325y/gIVstmLqitJm55uafctyHv8Uv77t/tkPRLHrBkx/Fe16xz4y1P2WCnOQGupHc06rqC2O27Tb+XuOrqouAi8aUndJbPhk4eYJ9zwTOHKd8BHjqhsQhDbORH97OXo/ejv0et8PUlTVvLdjCR7pL4/n+mp9x61338N/33ZVHPGTL2Q5Hs2TvXbab0fanM4L8oqr67ngbqmr1eOWSNl4VPGfPnXjTIXvNdiiSNHRG78R/3fN+gyfsvO2sxqL5azpzkMdNjiXNjAJ8H6Qkja977QGM/yAradPY2BeFSJohVUXs+CVpUg4kaCaZIEtDxhFkSZrY6ACy3aRmkgmyNGSq7PglaSLVZiHHkQTNIBNkaRjZ8UvSuBxB1iCYIEtDZPTmEzt+SRrfugTZjlIzyARZGiJ2/JI0ufXPsLCj1MwxQZaGiB2/JE1u3ZU2u0nNIBNkaYjY8UvS5GrqKtKDZoIsDREffy9JU3AqmgbABFkaIs5BlqTJ+Zg3DYIJsjRE7PglaXI+5k2DYIIsDZFycp0kTWrdVDQzZM0gE2RpCNnxS9L41o8g21Fq5pggS0PEjl+SJrd+KtosB6J5zQRZGiJ2/JI0OecgaxBMkKUhYscvSZNbd6uGHaVm0EAT5CSHJLkuyaokJ01S72VJKsmStn5UkhW9z6+S7NO2XdLaHN2286DOR9rUvPlEkqYw+kIlM2TNoIWDOlCSBcDpwEHAamB5kmVVdc2YetsCJwKXjZZV1TnAOW3704ALqmpFb7ejqmpkhk9BmnFlxy9Jk3IgQYMwyBHk/YFVVXVDVd0LnAscNk69dwDvBn45QTtHAh+fmRCl2WXHL0mTcyqaBmGQCfKuwE299dWtbJ0k+wK7V9WFk7TzCn49Qf5wm17xlviGBc1hPgdZkia37kqbv+41gwaZII/3k7x+rn2yBXAa8IYJG0ieBfyiqq7uFR9VVU8DntM+R0+w7/FJRpKMrFmzZmPil2beuldN2/FrsOwjNVesu9I2q1Fovhtkgrwa2L23vhtwc299W+CpwCVJbgQOAJaN3qjXHMGY0eOq+lH7ehfwMbqpHL+mqs6oqiVVtWTRokUP8lSkmbHuMW+zHIc2P/aRmivWTbGwo9QMGmSCvBzYM8keSbaiS3aXjW6sqjuraqeqWlxVi4FLgaWjN9+1Eeb/QTd3mVa2MMlObXlL4MVAf3RZmlNGO/4t7PglaVzr79Wwo9TMGdhTLKpqbZITgIuBBcCZVbUyyanASFUtm7wFngusrqobemVbAxe35HgB8EXggzMQvjQQdvySNLn1c5BnORDNawNLkAGq6iLgojFlp0xQ98Ax65fQTbvol/0c2G+TBinNIjt+SZoeu0nNJN+kJw0Rbz6RpMmVNzNrAEyQpSFSPghZkiblzcwaBBNkaYjY8UvS5HyKhQbBBFkaJnb8kjSp9VPR7Cg1cwZ6k95c9tvv+jJ333v/bIehee7+NjSyhRmy5pivfPdW3vjJ78x2GNoM/KL9Lrab1EwyQZ6mg/d+NPfeb4Ksmbflgi14/l47z3YY0gZ51HbbcOjTHj3bYWgz8bgdH8Y2Wy6Y7TA0j5kgT9MpL9l7tkOQpKG192O2452HP222w5CkTcI5yJIkSVKPCbIkSZLUY4IsSZIk9ZggS5IkST0myJIkSVKPCbIkSZLUY4IsSZIk9ZggS5IkST0myJIkSVKPCbIkSZLUY4IsSZIk9ZggS5IkST0DTZCTHJLkuiSrkpw0Sb2XJakkS9r64iR3J1nRPv/Yq7tfkqtam+9LkkGciyRJkuanhYM6UJIFwOnAQcBqYHmSZVV1zZh62wInApeNaeL7VbXPOE3/A3A8cClwEXAI8LlNHL4kSZI2E4McQd4fWFVVN1TVvcC5wGHj1HsH8G7gl1M1mGQXYLuq+mZVFXA2cPgmjFmSJEmbmUEmyLsCN/XWV7eydZLsC+xeVReOs/8eSa5I8tUkz+m1uXqyNnttH59kJMnImjVrNvokJGk+so+UpPUGmSCPNze41m1MtgBOA94wTr1bgMdW1b7A64GPJdluqjYfUFh1RlUtqaolixYt2uDgJWk+s4+UpPUGNgeZbnR39976bsDNvfVtgacCl7T77B4NLEuytKpGgHsAquryJN8Hntja3G2SNiVJkqQNMsgR5OXAnkn2SLIVcASwbHRjVd1ZVTtV1eKqWkx3093SqhpJsqjd5EeSxwN7AjdU1S3AXUkOaE+vOAa4YIDnJEmSpHlmYCPIVbU2yQnAxcAC4MyqWpnkVGCkqpZNsvtzgVOTrAXuB15TVbe3ba8FPgI8hO7pFT7BQpIkSRttkFMsqKqL6B7F1i87ZYK6B/aW/wX4lwnqjdBNzZAkSZIeNN+kJ0mSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9JsiSJElSjwmyJEmS1GOCLEmSJPWYIEuSJEk9A02QkxyS5Lokq5KcNEm9lyWpJEva+kFJLk9yVfv6/F7dS1qbK9pn50GciyRJkuanhYM6UJIFwOnAQcBqYHmSZVV1zZh62wInApf1im8DXlJVNyd5KnAxsGtv+1FVNTKjJyBJkqTNwiBHkPcHVlXVDVV1L3AucNg49d4BvBv45WhBVV1RVTe31ZXANkm2numAJUmStPkZZIK8K3BTb301DxwFJsm+wO5VdeEk7bwUuKKq7umVfbhNr3hLkoy3U5Ljk4wkGVmzZs1GnoIkzU/2kZK03iAT5PES11q3MdkCOA14w4QNJE8B3gX8Ua/4qKp6GvCc9jl6vH2r6oyqWlJVSxYtWrQR4UvS/GUfKUnrDTJBXg3s3lvfDbi5t74t8FTgkiQ3AgcAy3o36u0GfAo4pqq+P7pTVf2ofb0L+BjdVA5JkiRpowwyQV4O7JlkjyRbAUcAy0Y3VtWdVbVTVS2uqsXApcDSqhpJsj3wWeDkqvrG6D5JFibZqS1vCbwYuHpwpyRJkqT5ZmAJclWtBU6gewLFtcAnqmplklOTLJ1i9xOAJwBvGfM4t62Bi5NcCawAfgR8cObOQpIkSfPdwB7zBlBVFwEXjSk7ZYK6B/aW3wm8c4Jm99tU8UmSJEm+SU+SJEnqMUGWJEmSekyQJUmSpB4TZEmSJKnHBFmSJEnqMUGWJEmSekyQJUmSpB4TZEmSJKnHBFmSJEnqMUGWJEmSekyQJUmSpB4TZEmSJKnHBFmSJEnqMUGWJEmSekyQJUmSpB4TZEmSJKnHBFmSJEnqMUGWJEmSekyQJUmSpJ6BJshJDklyXZJVSU6apN7LklSSJb2yk9t+1yV54Ya2KUmSJE3HwkEdKMkC4HTgIGA1sDzJsqq6Zky9bYETgct6ZXsDRwBPAR4DfDHJE9vmKduUJEmSpmuQI8j7A6uq6oaquhc4FzhsnHrvAN4N/LJXdhhwblXdU1U/AFa19qbbpiRJkjQtAxtBBnYFbuqtrwae1a+QZF9g96q6MMkbx+x76Zh9d23Lk7bZa/t44Pi2+rMk123wGcBOwG0bsd9smWvxwtyLea7FC3Mv5rkWL8xezI/b2B030z4S5l7Mcy1emHsxz7V4Ye7FPJvxTqufHGSCnHHKat3GZAvgNOC4Ddh3vBHwGqeMqjoDOGPKKCeRZKSqlkxdczjMtXhh7sU81+KFuRfzXIsX5mbMm2MfCXMv5rkWL8y9mOdavDD3Yp4L8Q4yQV4N7N5b3w27eEHxAAAJoElEQVS4ube+LfBU4JIkAI8GliVZOsW+k7UpSZIkbZBBzkFeDuyZZI8kW9HddLdsdGNV3VlVO1XV4qpaTDelYmlVjbR6RyTZOskewJ7At6ZqU5IkSdpQAxtBrqq1SU4ALgYWAGdW1cokpwIjVTVhYtvqfQK4BlgLvK6q7gcYr80ZPI0HdflxFsy1eGHuxTzX4oW5F/NcixfmZsybwlw877kW81yLF+ZezHMtXph7MQ99vKkad8quJEmStFnyTXqSJElSjwmyJEmS1GOCPA3D9DrrJGcmuTXJ1b2yHZN8Icn17esOrTxJ3tfivjLJM3v7HNvqX5/k2BmMd/ckX0lybZKVSf5kDsS8TZJvJflOi/ntrXyPJJe145/Xbgyl3Tx6Xov5siSLe22N+4r0GYp7QZIrklw4R+K9MclVSVYkGWllw/xzsX2S85N8t/08P3uY4x20Yekn51of2Y41p/pJ+0j7yElinj/9ZFX5meRDd/Pf94HHA1sB3wH2nsV4ngs8E7i6V/Zu4KS2fBLwrrb8IuBzdM+RPgC4rJXvCNzQvu7QlneYoXh3AZ7ZlrcFvgfsPeQxB3h4W96S7rXnBwCfAI5o5f8IvLYt/zHwj235COC8trx3+3nZGtij/RwtmMGfjdcDHwMubOvDHu+NwE5jyob55+Is4NVteStg+2GOd5AfhqifZI71ke14c6qfxD7SPnLimOdNPznQg83FD/Bs4OLe+snAybMc02Ie2PlfB+zSlncBrmvLHwCOHFsPOBL4QK/8AfVmOPYLgIPmSszAQ4Fv072h8TZg4difC7qnqDy7LS9s9TL2Z6Vfbwbi3A34EvB84MJ2/KGNt7V/I7/e+Q/lzwWwHfAD2o3Nwx7voD8MWT/JHO4j2/HmTD+JfaR95Pq251U/6RSLqY33iuxdJ6g7Wx5VVbcAtK87t/KJYp+Vc2qXqfalG20Y6pjbpbgVwK3AF+hGCn5SVWvHOf662Nr2O4FHDjjm9wJvAn7V1h855PFC99bLzye5PN1rjmF4fy4eD6wBPtwu0f5TkocNcbyDNuznNWf+neZKP2kfOePxwtzqI2Ge9ZMmyFOb9BXZQ26i2Ad+TkkeDvwL8KdV9dPJqo5TNvCYq+r+qtqHbtRhf+DJkxx/VmNO8mLg1qq6vF88ybGH4nsM/FZVPRM4FHhdkudOUne2Y15Id9n+H6pqX+DndJcKJzLb8Q7aXD2vofp3mkv9pH3kr+0zE+ZSHwnzrJ80QZ7aVK/IHgb/mWQXgPb11lY+UewDPackW9J1+udU1b/OhZhHVdVPgEvo5kdtn2T05Tr946+LrW1/BHD7AGP+LWBpkhuBc+kuIb53iOMFoKpubl9vBT5F90t2WH8uVgOrq+qytn4+3S+CYY130Ib9vIb+32mu9pP2kfaRPfOqnzRBntpceJ31MuDYtnws3fy10fJj2p2iBwB3tssbFwMHJ9mh3U16cCvb5JIE+BBwbVW9Z47EvCjJ9m35IcALgGuBrwAvmyDm0XN5GfDl6iZOTfSK9E2qqk6uqt2qe0X7Ee34Rw1rvABJHpZk29Flun/PqxnSn4uq+g/gpiRPakX/le7NnkMZ7ywY9n5yqP+d5lo/aR9pHzmeeddPDnrS81z80N1p+T26OVZ/McuxfBy4BbiP7q+sP6CbG/Ul4Pr2dcdWN8DpLe6rgCW9dn4fWNU+vzeD8f423aWRK4EV7fOiIY/56cAVLeargVNa+ePpOsNVwCeBrVv5Nm19Vdv++F5bf9HO5Trg0AH8fBzI+ju0hzbeFtt32mfl6P+rIf+52AcYaT8Xn6a7u3po4x30hyHpJ5ljfWQ71pzqJ7GPtI+cOO5500/6qmlJkiSpxykWkiRJUo8JsiRJktRjgixJkiT1mCBLkiRJPSbIkiRJUo8JsiRJktRjgqx5JcmBSS5sy0uTTPaaS5Icl+Tvp9HujUl2asv/Z4q6f74hMW9KLc6rkizpre+0kW29Ismq0e+nJEmbCxNkDb32lp0N/lmtqmVV9TebOp6q+i9TVJm1BLl5XlWNPNhGquo84NWbIB5JkuYUE2QNpSSLk1yb5P3At4Hdk/xDkpEkK5O8vVf3kCTfTfJ14Hd75etGh5O8JMllSa5I8sUkj5ri+I9M8vlW/wN0b/wZ3faz9nWXJF9LsiLJ1Umek+RvgIe0snNavU8nubzFfXy/nSR/leQ7SS4djSnJo5J8qpV/J8l/aeWvSvKt1vYHkizYgO/nQ5L8W5I/bN/b7yb5pxb3OUlekOQbSa5Psv9025UkaT4yQdYwexJwdlXtW1U/pHvV5hK615z+TpKnJ9kG+CDwEuA5wKMnaOvrwAFVtS9wLvCmKY79VuDrrf4y4LHj1HklcHFV7QM8A1hRVScBd1fVPlV1VKv3+1W1H7AEODHJI1v5w4BLq+oZwNeAP2zl7wO+2sqfCaxM8mTgFcBvtePdD4y2P5WHA58BPlZVH2xlTwD+ju57uVc7l98G3sjsj4BLkjSrFs52ANIkflhVl/bWX95GYBcCuwB70/2R94Oquh4gyT8Dx/9aS7AbcF6SXYCtgB9Mcezn0kajq+qzSe4Yp85y4MwkWwKfrqoVE7R1YpL/3pZ3B/YEfgzcC4zO770cOKgtPx84ph37fuDOJEcD+wHLkwA8BLh1inMYdQHw7qo6p1f2g6q6CiDJSuBLVVVJrgIWT7NdSZLmJUeQNcx+PrqQZA+60c3/WlVPBz4LbNM21zTa+v+Av6+qpwF/1Nt3MpO2W1Vfo0ukfwR8NMkxY+skORB4AfDsNiJ8Re/Y91XV6DHuZ/I/WAOc1Uam96mqJ1XV26ZxDgDfAA5Ny6ybe3rLv+qt/2qKOCRJmvdMkDVXbEeXMN/Z5uoe2sq/C+yR5Dfa+pET7P8IukQW4NhpHO9rtCkMSQ4FdhhbIcnjgFvbtIUP0U2HALivjSqPHveOqvpFkr2AA6Zx7C8Br23HWJBku1b2siQ7t/Id2/Gn4xS6Eev3T7O+JEmbNRNkzQlV9R260deVwJl0o6JU1S/pplR8tt2k98MJmngb8Mkk/w7cNo1Dvh14bpJvAwcD/3ecOgcCK5JcAbyUbk4vwBnAle0mvX8DFia5EngHcOk47Yz1J8Dz2nSHy4GnVNU1wF8Cn29tfYFumsl0/SmwTZJ3b8A+kiRtlrL+Cq+kuS7JjcCSqprOHwHTae9A4I1V9eJN0Z4kSXOBI8jS/LIG+NLoi0IejCSvoJuWMd4NipIkzVuOIEuSJEk9jiBLkiRJPSbIkiRJUo8JsiRJktRjgixJkiT1/P9QcjLXDJEhOAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, [ax0, ax1] = plt.subplots(1,2, figsize=(10,4), sharex=True, sharey=True)\n", "\n", "ax0.plot(r, y_e)\n", "ax0.set(xlabel='radial distance [km]',\n", " ylabel='$Y_e$',\n", " ylim=(0.4,0.6),\n", " title='PREM (1981)')\n", "\n", "ax1.plot(r, y_e_c)\n", "ax1.set(title='Const. Layer Approx.')\n", "\n", "fig.tight_layout()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }