{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# [Prody](http://prody.csb.pitt.edu) is an open-source Python package for protein structural dynamics analysis with an expansive and well-written API." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "scrolled": true }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "@> PDB file is found in working directory (3ub5.pdb.gz).\n", "@> 4055 atoms and 1 coordinate set(s) were parsed in 0.06s.\n" ] } ], "source": [ "import prody\n", "prot = prody.parsePDB('3UB5')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The default visualization in ProDy uses matplotlib." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcUAAAE1CAYAAACWU/udAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvWlwI+d17v80NoIgMNxnOMOZ4c4hOftwVtsla7FTtpMopdiRlDjXdq5TtrNUKfEXqyp2Konvtcex44rL/mBXRbF0XUl8Fd+ULTlO/WNJlmzLkmbRjKTReEhi4QaCBEnsaAANoPv/gfO2uhvdQHdj4fb+qqZUAtnoBgi8T5/znnMeRhAEUCgUCoVCASybfQEUCoVCoWwVqChSKBQKhXIXKooUCoVCodyFiiKFQqFQKHehokihUCgUyl2oKFIoFAqFchcqihQKhUKh3IWKIoVCoVAod6GiSKFQKBTKXWwGf5+Ov6FQKBTKdoTR80s0UqRQKBQK5S5UFCkUCoVCuQsVRQqFQqFQ7kJFkUKhUCiUu1BRpFAoFArlLlQUKRQKhUK5CxVFCoVCoVDuQkWRQqFQKJS7UFGkUCgUCuUuVBQpFAqFQrkLFUUKhUKhUO5CRZFCoVAolLtQUaRQKBQK5S5UFCkUCoVCuQsVRQqFQqFQ7kJFkUKhUCiUu1BRpFAoFArlLrbNvgAKpZEIgoBisQiGYWCxWMAwusy4KRTKLoGKImVXQMSwUCggl8tBEAQAgMVigcVigc1mg9VqFf+fiiWFsjthyOKgE0O/TKFsNlIxFAQBDMMgn89D+rkXBAHK74HFYoHVahX/UbGkULY9ur68VBQpOxI1MSSCxnGc+JjWseS/0u9HNptFPp9HR0cHFUsKZfuh60tK06eUHYUgCCgUCigWi6LwWSzG6smIwCmFLpvNIhqNYs+ePcjn87Kf0ciSQtkZUFGk7AiIGBYKBQDQFMNMJoO5uTlYrVa43W60tLSgqalJl3iR31E+L4kmC4UCFUsKZZtDRZGyrVETQzXBYVkWfr8fyWQSBw8eBM/ziEajWFxcRC6Xg9VqRUtLi+yfw+GQPRfDMCV7j+Rx6X+l1wZsiCXHcbKfU7GkULYmdE+Rsi0hYuj3+3H48GFNMUyn0/D7/WBZFgMDA+ju7kahUADP87LfLxQKYFkWqVQK6XQa6XQaHMfBZrOJIsnzPFiWxdjYWNXXLv3HMAw4jkM0GsX+/fupWFIo9YHuKVJ2HjzPiwU0ALC0tIT+/v6S30ulUvD5fMhmsxgaGkJnZ2dZYbHZbNizZw/27Nkje7xQKCCdTiOVSiEejyORSODKlSuw2Wxi+pX8s9vtul6DmoCTyLWnpwf5fF41srTZbGKEScWSQqkPVBQp2wKlGGpFhslkEj6fD/l8HoODg+jo6KhKOGw2G1pbW9Ha2gq3243l5WUcOXIE+XxejCjD4TBSqRQKhQLsdnuJWNps+r9mavugJKJUq5pVRpVWq1XzvaFQKJWhokjZ0vA8L1aTAtpiGI/H4fP5wPM8hoaG0N7eXtfrstvtaGtrQ1tbm+xxjuNEsVxeXkYqlUKxWITD4UBLS4tMMK1Wq65zab1mQRDEm4VyYkkEk4olhVIZKoqULYleMSwUCrh+/ToYhsHQ0BBaW1vrdk1ahTZSHA4HHA6HTJRJlEfEMhgMgmVZFItFNDU1idEkx3EoFot1EcvZ2Vn09/dTsaRQKkBFkbJlIGnCfD4PnucBaC/8kUgEPp8PHMfh5MmTdRVDgh5R1DquqakJTU1N6OjoEB8XBAG5XA7pdBrRaBSZTAY3btwAz/NwOp1iROl2u+FyuXT3W6q9Z5FIBIODg6JYKqFiSaFsQEWRsunoFUNBELC+vg6/34+mpiaMjY3h1q1b8Hg8m3HZVcMwDJxOJ5xOJ5qbm5HNZnHs2DEIgoBsNitGlpFIBCzLgud5NDc3y8Syubm5KrEE5JGlEiKSpMiHiiVlp0NFkbJpkMW4UCjg9u3bGBsb0xTDtbU1+P1+NDc3Y2JiAm63G4D56M0MjToXwzBobm5Gc3Mzurq6xMcFQUAmkxHFcm1tDSzLAoBMLFtaWmoulhzHyX6fFPVQsaTsNKgoUhqOVAxJv2AsFlOdFBMOhxEIBOB2u3H8+HG4XC7Z75gRKrMLdyMFWOv8LpcLLpcL3d3d4uM8z8vEMhwOI5PJAABcLhdaWlrEPszm5mbdr1+vWLIsi3g8jgMHDlCxpGx7qChSGoaaGGpFhisrKwgEAmhtbcXJkyfR3Nys+pwWi0VMue5WLBaLGCFKIcMG0uk0isUifD4fMpkMLBaLKJbkn9PpNC2WxWIRLMuCYZiykaWaPRcVS8pWg4oipe5oOVaoiWEoFMLs7Cza29tx+vRpOJ3Oss+9E9OntcJiscDtdsPtdmN+fh7Hjx8H8I6IpdNpxONxLC0tIZvNysSV/NMzF7bc35T8nOd55HI52eMMw6hO7zE6wJ1CqSVUFCl1Q00M1RY8nufBcRxeeeUVdHZ2YnJyEk1NTbrOYTZSLGcdpcV2E0UtrFYrPB5PSYFSsVgUU7BG5sJWei+1boDIOclABunvaw0loFDqDRVFSs0xIobBYBDz8/PgeR5nz56Fw+EwdK6dIlQANv11WK1WzVF3ZC7s+vo65ufnZXNhye9wHKf771dpiLpULKPRKACIPpbS/UoqlpRaQ0WRUjP0imGxWMTi4iIWFxexb98+nD9/HteuXdM9O1RKI/cU6ynAW3lh15oLm8/nwbIslpeXwbIs3n77beTz+arnwkr/C2z4WJJoUxlZCoIgK+6hc2Ep1UJFkVI1xLEiHo9jdXUVAwMDmmK4sLCAxcVF7N+/HxcuXBDngpIiDb3TXAhGhYrjOAQCAeTzeXG/zeVy6ZpPupOi0lpgt9vR2toqRogDAwMANsSSuI2srKwgnU5XNRdWEATNma7KyFKayiVpWCqWFCNQUaSYRullyPM80ul0yaJTKBQwPz+PpaUl9Pb24uLFiyWLIYn4jIqi3kiR4zjMzs5ibW0NBw8exJ49e8CyLJaWlsTqzGqmyOxmlHuKdrsd7e3tJfNnpaPuQqGQ+L47HI4SsZR+Dnie1/w76E3DSq+RellSykFFkWIYLWNfq9UqE6h8Po/5+XksLy+jt7cXly5d0hQ9s2nQStEbiQzX19fR19eH4eFhcXpOZ2en7DVJp8isr6+DZVkIgiC2LzidzpIFlqK/aMnsXNhMJiOmTo3MhZX+V3o+YONGLZ/PA9hwVuE4Dt3d3VQsKVQUKfqp5HJPhI3jOMzNzSEcDuPQoUO4ePFixcWs1qKoFMORkREx2lAbZ6Y1RUbaGJ9IJMCyLK5evSr2+pEIx+12yyoydxPV3CTomQubTCaxvr6OcDhck7mw0v8CG3uWpF1EKpYEGlnuLqgoUipSSQwJhUIBiUQCV69eRV9fHy5duqR7sbJYLKb265RiWk4MzSDt3Wtvb0cikcCZM2fEXr9UKlXSviBNBbrdblMFRNuJekTO0rmw0WgUHR0d6OjoUI3oM5mMbC4sef/1jrojaftye5ZULHcPVBQpmug19s1ms5idncX6+jqsVqshMSSQQhujkEix1mKodS6CVq9foVBAOp1GKpXC6uoqAoFASZGJkeKe7UC908nK4hk9c2FXV1dV58KSIerS661mz5KK5c5jZ3wrKTVFr5dhJpNBIBBAPB5Hf38/hoeHcf36dVNiVE0T/sLCAjKZDPr7++sihgQ91ac2mw2tra0lVlYcx4kVmcFgEOl0WpYKdDgc4vi77VbcU29R1POeVDMXNpPJoKWlxdDrqCSW+XweHMeJn5lQKISDBw9SsdwGUFGkiBAxvHPnDgYHB2Gz2VS/sCzLwu/3I5lMYmBgAOPj4+KXX22/Tg9GRZFEhqFQCD09PTh58uSWFhOHwyGmAAnSVGAsFgPLsrh+/ToAqKYCt+riSXoF6/n8Zl+7nrmwkUgE6XQay8vLNZkLC0C2h14oFBAOh9Hb2ysTS+k1UrHcOlBR3OWoeRkmEgnVhSidTsPv9yOdTmNwcBBHjx6V/U41X2K9opjL5cRUbX9/Pw4fPoyWlpaGCGKt+xSlqUC32w2WZXHixAlZdJNMJrG8vCwb5L3Vinu2QqRoFOlc2FQqhdbWVnR1demeC2vkvSd7lmqvQfr9UxNLaY+l1r4npbZQUdyllDP2tVqtsogvlUrB5/Mhm81iaGgInZ2dNf9iViq0UYrh6OgoGIZBIBDYERNtpEgX4L1794qPb9XinkbuKdYDaX9srefCkmO1qq/LDVEne+XK1681F5aKZW2gorjL0ONyb7FYUCwWkUwm4fP5kM/nMTg4iI6Ojrp98bQKbXK5HAKBACKRiEwMKx1X6VzbkWqKe4rFoqE+PyNsx0jR6POXmwsrrYRVzoU1m8Wo5DhSLBZl73sqlYLdbofH45EJJhVL41BR3CXo9TIENu5sb9++DYvFgqGhoZLJJPVAmT5ViuGRI0dUr9VsK4cZturiolbcQ258SHFPPp/HzZs3ZZN7iGBWO7mnEZHiZouiFlqFVWQubCqVQiQSQTwex5UrV6qaCwtoi2UkEoHL5UJTU1PZyJKKZWWoKO5wjIhhLBaDz+dDKpXC4OAgDh061LDrJKKoVwwJdB6pOgzDyIp7lpeXMTk5KSvuSaVSWFtbK2ldMFrc04hIcbtFomQubGtrqyhWIyMjNZ8LSygWi+L+oxRpZKmEiKTUdYSKJRXFHYteY19g4y7T5/PBZrNhZGQES0tLcLlcDb1eYiMVCAQwMDBQUQwJFoulpE+Moo2eyT1Gi3t2Q/q0GorFovj8anNhlVG9dC4sGXWnNRdWeg41EdWThuU4Tvaz3S6WVBR3GEQMQ6GQWFquVfVGxNDhcGBsbEzcq1pZWTHdWkGeW+8XiESGy8vL6OrqwuTkpKEv306JFDd7wammuIf0+dWLnSC65fZylVE9QToXNpVKlcyFlUaWhULB0H6xUbFkGAbZbFaMeneyWFJR3CEoI8NYLAaGYeB2u0t+b21tDX6/H83NzZiYmCj5nWo8CvW6XUjTpAMDA2JkavQL1kg/xXqzFcVdT3EPKTSZm5uDw+GQVcFqRTZG2AmRopn3QM9cWNJnub6+jmQyKd4Im90vLieWPp8P/f394g3zwsICvvWtb+E73/mO4de2laGiuM3RMvZVtlUIgoDV1VX4/X643W4cP35cM0WqPNYIlUQxl8vB7/cjGo3K0qTBYNBUGtRopJhKpeD1epFKpcT+wO3QIL/VkBaYsCyLrq4utLW1ydKAapN7qlms68VWFUUtpHNhidPLjRs3cPToUbF1RLpfLAiC6bmw0nOSCJW8llgsJk4H2klQUdymkCHd0tJs6YecCJsgCFhZWUEgEMCePXtw8uRJNDc3l31u0pJhBqV9FEEphmNjYyVtIGYiPr3Vp+l0Wuy1HBwchNvtRjabRSqVKtlDky4eZA+Nog15/8ulAcsV90jTgJtxY9KISLTec27JnqLD4Sg7FzaVSslG3RkprioUCrLXEY/HS1pUdgJUFLcZao4Val9oi8WCSCSCubk5tLe34/Tp03A6nbrOYbVaTRevKAU1m80iEAhoiqH0OLMDwcsdR0bSpdNpcfAAKWwgszKVe2hk8VhfX8fc3Bzy+bxYHZjP55FIJGqSFtxJlBMVM8U90sISnufrvq9Yz+euV3+olEpDzSvNhU2lUlheXkY2mwUA2ag7t9sNp9NZ0hqTSCTQ1tZW19e1GVBR3CbotW/ieR6hUEjcM5ycnERTU5Ohc1mtVvHLYRQiblIxHBwc1BRD6XFm9tS0RDGTyYjzWYeGhtDV1SWev9x5tJq0yUDvcDgsSwsq77RdLteuS8GaFSy9xT3ZbBZXr16VNcST93s72HI1QhSrff+lSOfCSm9WWJbF7du3IQgCbt++jZWVlaoixampKTzyyCPi//v9fvzd3/0dPvaxj+GRRx7B7Ows+vv78fTTTzekV5pARXGLY0QMg8Eg5ufn0d3djeHhYWSzWcOCCFRXvEI25Ml81EpiSDAzmQYoFdNsNgu/3494PI6hoSFMTEzURKRIWtDhcGB8fByAelqKZdktO6O0XtS6j1BZ3BONRnHu3Dnk83mZ20UqlUKhUBCLe6Rp2K0UxespPKuWWhZpSefCSrly5QoOHTqEubk5/PKXv8TNmzextraGH/7whxgfH8fp06fx53/+57rPc+TIEdy8eRPAxo1Db28vHnroIVy+fBkPPPAAHn/8cVy+fBmXL1/GV77ylZq9vkpQUdyiCIKAZDJZcbZhsVhEMBjEwsIC9u7di3PnzsHhcGB1dRXpdNrUuc0U2pDIMBKJoK+vDydOnDC0UFaTPiXVeGTPcnBwUHTu0DrGLNLFRystJY10IpEIFhYWxNFfUqE02qBdi2uu1/M3QvDtdjva2tpkKTtp20Kti3tqhbRPsR40omqZ3Ph4PB4cO3YM3/jGN3D58mWcOnUKH/rQh3Dnzh0sLCyYfv7nn38eQ0ND6Ovrw49+9CO8+OKLAICPf/zjuPfee6ko7makxr4zMzPo7e1VTR0Ui0UsLCwgGAyip6cH58+fl6WSbDZbQ2yciBjGYjEMDAxAEAS0t7c3rLWiWCwiGo3i+vXrZfcsG4lWGwOpzEylUrIGbafTiebmZrHM3mhl4Gazma0k5doW9BT3uN3uul9/vdOn9R6DB6gPB4jH42hra0NLSwsmJycxOTlp+vm///3v4/d///cBbPRJ79+/HwDQ09ODlZUV8xduAiqKWwQ1Y1+bzVYiFIVCAfPz81haWkJvby8uXLigGm1UkwLVEylK05RSMYrH46bToEb9FGdnZ7GysgK73Y4LFy40RAyrOYfWNJNsNotYLIa1tTUEAgGxMlC6f1ZNCrYRkdxm34goqVTcQ6qOQ6EQWJbFtWvXTFtDVaLeomi0cd/sOZTrTCKRKJn5agaO4/DMM8/gy1/+csnPNmM4ABXFTaacy71UnPL5PObn57G8vIze3l5cunSp7BfBarWK+5BGKSdQSjFUpim1WjL0nFPPHXs+n8fs7CzC4TD6+vpw6tQp+Hy+Lbco64Us3larFSsrKzh27BiAjc+FmlURKTaR9lc2IgVbjnqKbq2jOGVxCc/zeP3113H69GlZcY805V1tcU8j+iDr/RlQE0USKVbLf/3Xf+HMmTPYt28fAGDfvn0IhULYv38/QqGQrAirEVBR3CR4ni9r3wRspECz2SxmZmYQDodx6NAhXLx4UdddoVlxIscqI8VKYkgw2+NYqdCmUChgbm4Oy8vLOHz4MC5dugSLxSI2J+80LBaLZgpWWkKfSqVKxn653W7Z/tl29jts1Ii3cilvcnMiHeBtpLin3pFiI6pbtSLFWojiv/3bv4mpUwB48MEH8dRTT+Hxxx/HU089hd/5nd+p+hxGoKLYQPR4GRI4jsPa2prYX0dEQC+1mEoD6BdDtWPNnlMKSReHQiEcOnSo5H1o9OzTzRZgrWKTXC4n7leur6+L+2ekECifz4vVyLUWmXruaW32XFIzxT1SoXS5XHUX9kaJovIctUifptNp/PSnP5WNinv88cfx8MMP44knnkBfXx+efvrpqs5hFCqKDcCIGEpngra2tqKrq8uUhVM1okia92/fvo14PF6xmlNKrUSxWCzK9k61IuSdNPvULNKxX8r9M5ZlkUwmUSwWMTU1VTLMuxb9fvW8UWiEoBgVXaPFPel0Grdv3y5phq/V61KL4mqN2jlIxFwNLS0tWF9flz3W2dmJ559/vqrnrQYqinVE6WUIaIshichIFeeRI0cQDoeRTCZNndusWGSzWZmnol4xlJ7XjBiT6yVVtYuLi2XFkNDoSHE77V2SfjOn04lQKISTJ08CkA/zlqYEiU2RNNLRIxj1FK7tNAxcq7jnypUr6O/v1xwpKH3PzRT3NCJSVO5bbnbGpJ5QUawDUjG8desWjh49qimGZPJKIpEoSU9WE+0Z/WJlMhkEAgExMozFYuLGtxGq8TfMZDJ49dVXsX//fly8eFHX3S+NFI2j5hZPUoIkBRuJRMQ+V2ULgzLK2Ql7ivWEYRjVyTFkpCBxupifn5cV90hvTspF8o1Kn6oZCGynm0S9UFGsIWqOFfF4XPWDk06nxZmcg4ODqpNXqhFFvSjFkIiyz+cz9XxGC3ykk3h4ntdsMdFiM/wU67VQ1+t16LleaUqQOC8A8haGRCKBUCiEbDYLq9UqLtgcx4nu8bVmO0WKapCtEzW0RgqWK+6RCiUp7tmMloxsNqt7lvJ2g4piDSjncq9cjFKpFPx+PzKZDAYHB2UzOZXYbDbTbRWVkEaoZtKkWuhNn/I8j6WlJczNzWHv3r04f/68ON/S6PmMRorFYhGRSAQul6umezvVsBWuQQ1pik+aOZCmYMn+szIFS6pgq1mwt3ukaKYIqVJxDzEcJsU9PM+LAul2u+sy/EEpirFYrCY9ilsRKopVoOVlKIW0GhDrIo7jMDQ0hI6Ojopf9npEikoxrNVsUEIlkRIEAaFQCLOzs+jq6hLH0pnFyLXzPI+FhQUsLCygra0NwWBQjHqk6UG3260pziQy3aoipkWtr1eagg2FQjh9+jQYhhFTsOl0GgsLC2BZFjzPl5jf6rWI2u6RYq1Sm+WKe6ampuBwOJDJZGSTe5TveTU3gEpRrFXj/laEiqIJ9Iih9Hdv3LgBQRAwODgo+0BXotpIkQiyxWKRiWEtB2Ur0RJFQRCwvLyMQCCAzs5OnD17tmE+hdKotKenBxcuXJANsVYWnvj9fhQKBdnsTHIHvhnp2mpp1OxTrRSsIAiqrgt6vCvrPcJsu4suWXs6OjpkkaUy7b20tIRsNmu6uGe3eCkCVBQNYUQMY7GYaGo7NDRkqmil2kjRarUinU5jfn7elBiaiYiUoig1OW5razNlZWUWaVTa3d0tzoclqSiCVuGJtPdvdXVVtM+ZmpqCx+PZVibEmznmTVpoYsS7sqWlpeyeXC2o97DuzWqs17KFkhb3rK+v6y7uUe5b1mqazVaEiqIOKrncS4lGo/D5fLBarRgeHsbc3Jxq1ZYeqqmsJAv4W2+9heHhYcORITm30S80KbQRBAHhcBh+vx+tra2GTI6rRXru9vZ2U1GpVu/f9evXceDAAWQyGdlCLi2CIP/dKkO9t2pkW8m7kizcqVQKV65ckXlXSiP3aqi3rVOj2iX0nsNscQ8ZVk/2LneqwTBARbEsgiCIZenEvkltoRMEAZFIBD6fDw6HA0eOHBHHRTWyrQKQ7xk6nU4cO3asxBdND+S6jX6hGYYBy7J47bXX4PF4cOrUKTQ3Nxs+vxkEQcD6+jq8Xi88Hk9dhJjcgSsXBGk7A9lLEwRBtq/jdrvrMlFmp0G8Kzs6OuByuRCLxTA4OCjzrlxZWRFTsMS70kyv32ZPzKkFtRDeSsU9y8vLCAaD8Pl8+MIXvoD29nZ0dnbi6aefxrFjxzAyMmKo+jgWi+GP//iPcevWLTAMg3/+53/GkSNHNtVcmEBFUQWpsa/X60V3d7csWpD+3tramuhyPzExUSJA9awglZLJZGRN9xMTE7h165bpSNNolEoEaWZmBrlcDufPnzcdIZshEonA6/XC6XTixIkTFc9da2GSLuQE6b5OPB5HMBgUh3orC3vqzXYVYiJaerwrlelAqVBqDU7neb6u02DqnZ4F6ueSId0jJuba4+Pj+OAHP4gvfelL4DgOXq8XP/zhD5FIJPDjH/9Y93M/9thj+MAHPoAf/OAH4DgOLMviS1/60qaaCxOoKEpQc7l3OBwloiYIAlZXV+H3+9HS0oJjx46V5O4J9RZFpRiSQQFA9aPe9B5LorPm5maMj49jZmamYYIYj8fFfVO1m5JaY6TQRqudQTrUm/gqFgoFcByHQCBguEKzEls1faqHSvvaWoO8pe0LS0tLqrNJ3W73jthTbES0q+aEc//99+PBBx80/HzxeBw///nP8eSTTwLYuKF0OBybbi5MoKIIubEvIB/FJhU1adHInj17dEUk1Vg4EdQWBpZl4ff7VcVQeu56RorRaBRerxcOhwNHjx6F2+2WjbQzipF2h2QyCa/XKy50p06dMnVOo9Si+lQtVZXP53Hz5k243W7RAYOkB6URj9vtNtwkvx1bSAhmF3yy0Kp5V0qLp2KxGKxWK9bW1mriXamkEaII1DcToGZNVU1LRiAQQHd3N/7oj/4Ib7zxBiYnJ/GNb3xj082FCbtaFMuJIcFms4HjOLGK0WjRSLWRIonYyIdSjxgSzM4hlZ5XjVgshpmZGdjtdoyPj8uis2qKg0gLSblFJJ1Ow+v1guM4DA8Po729Hb/61a8Mn2urCQXDMLBareju7i5JD5KIZ3V1FYFAQGyS17KK2knU8u8knU1K3mOv14u2tjY0NTWpeikqi6eMplobNUaunmh5KZoVxUKhgNdffx3f/OY3ceHCBTz22GO4fPmy7Hc2w1yYsCtFsZyxr/L34vG4eAdz5swZw+0ENpsNuVzO9LUSceI4ThTDoaGhsmIoPXc19lHKY+PxOLxeLywWC8bGxkpSVuQ4s5FUuYpXkiZOp9MYHh6W9cE1kkb3KapVCyrnlEqtoqTRznZpFylHI1KDNptN00uRVMGSNHexWJT1rxJ7KK1rbIQBcL0/j1q2UWaLYA4ePIiDBw/iwoULAICPfOQjuHz58qabCxN2lSgaEUMyj5P0Vo2NjZk6Zy32FH/961+L/Y56xJBQbaRIIr5EIgGv1wtBEDA8PFy3SRZqUWYulxPdQ4aGhtDd3a36+s1EFGaO2QrN++XmlJKoUjpgmgxpD4VCNRm91kg2s3nfbrejvb1dNQWrtIcCIFbBEsFsamoSDaA34/prRa0Nhnt6enDo0CFMTU3hyJEjeP755zExMYGJiYlNNRcm7BpRJM3YQHkxXFxcxMLCAvbu3Ytz584hk8lgYWHB9HnNiiJJk5Km+8OHDxtewKsttEmlUlhcXESxWMTw8HDd+5KkokiKTtbX1zE4OIixsbGyaWIjpe/k729W3DZbFLWwWCyqEU88HofP50OhUMDi4iLS6TQEQRDdL2oxBqxebLXqUC17KOJdqaw0zufzcLvd4DiuJt6Vate/GcMB0ulM2/XyAAAgAElEQVR0VUVt3/zmN/HRj34UHMdhcHAQ3/3ud8Hz/KaaCxN2jShqDegGIHr4BYNB9PT0iJNPAMiqUc1gVJike4ZDQ0NgGAZ79uwxtVhZrVbZ5Ba9kEIPQRBw9OjRhvUKkYgmGAxiZWUF/f39GBkZqbhoNTJ622qioQebzYampiaZWbV09Jp0DBiZAysVy3qn/8qxXcawSQuipExNTcHlcoFhGE3vymr2hDfDIYNMGarmfTt16hSuXbtW8vhmmgsTdo0oAqX7XYVCQRTDAwcOqNoW2Ww20/6A5Hg9osqyrLhnJk2TRqNR06JsVJDJ0PJsNouOjg643e6GCSIpKLl58yb6+/tx6dIl3V+6RnoqboX0aS3QGr1WKBTEvcrl5WWkUilxH00qlrVqF6lEI1wy6ikqgiCgtbW1ZE+YTIhR7gkbjd43w2CYsB1vEPWwq0SRUCgUMDc3h1AohIMHD+LSpUuaH6xq9wQrHU/EkGVZVSupRvQaSq9heHgYHR0dCAaDDRk6QJwrFhcXwTAMjh8/bliId4pQ1fM16F3AbDab6mQTaStDOBxGJpMRhZXjOESjUVPtIpXYLpGiFmrpWekIQS3vyng8LovelYPTyfvciEIepcFwI/YxN5NdJYqFQgF+vx/hcLiiGBKqqaYEtIWpkhhWOr6acxOkjf9DQ0OyazCbeiVUusNXc66YmpoydfdpJFIk+8bz8/MlqUK3213x81AvAa7nXXe116vWygC8E91Ho1Gsra1hdnZWNge22tQgufbt7KdoJBLV410pbctxOBzidzyZTNatgEqZPk0mk6qV5zuFXSWKPp8Pzc3NhlJz1aIUVb1iSKgmUtUSxWw2K/NU1Gr8r6adQ2sxk1pIdXV1yfZvzaZB9QgVccwIBALYt28fzpw5I862JZNlUqkUeJ6XpbDcbrcshbVTotJaQG4qHA4HRkZGAJRvF5HOKNXbIL8TIsVqhUrLxYXjOCwuLiKZTGJhYUEsoDLrXanFbrKNAnaZKI6NjTVs70lJOp2G3+8Hy7IYGhpCZ2enrg9qNRGbUtik7Q2Dg4MYHx+vWNFpBtIKIl1spM4VWhZSZqNy0vSvBhnJ5/P50NbWJpoacxwHQRBUU4UkhZVMJhEKhZDNZsVGblJhuN3aGuoVbSmfu1y7CHnvotEoFhcXZXNgpfto0ve1EUMW6j0Nph6iS95nl8sFh8MhFlGRFCzxriSfX6l3pXRwuh6UxTyJRIKK4k7B7IdfatZrlHQ6jUwmg1u3bhkSQ0It0qe5XA6BQACRSKRie0MtzisVVKlzhdvtLjsNqJy4VTqfmphGIhHMzMygpaVFt1uHdPC0tACFNHInEgmsrq4iGAyKjvLSRX23uWDoFS2t6kzyvqZSKQSDQXFGKYnWWZYVb2C24/vaaGsqaQpWy7tydXVVTHUT70rp4HTl9SojxVgstmNto4BdJopmISlMI9NBSCVnJpOB3W7HuXPnTIlqNelTstdw/fp19Pf348iRI4Ya/6udmxqNRjEzM6PbuaJW6dNEIoHp6WlYrVZxJmulYypBGrk9Hg96enrQ2toqtjVouWBoRT+bQaMiRaNoNciTaD0cDmNxcRGBQEBWcEL+bWa7iF7qKeaFQkHXzV4l70qtm5KWlhbkcrmSSLFeAzy2Alv/E1VDzH44jYiiVAxJZHjlyhXTkaaZiC2fz2N2dhbhcBgMw+DixYuGz11NpFgsFvHmm2/C4XAYcq4wK4rkuHQ6jZmZGRQKBYyMjNT9iytta1C6YKgtNNKo0uPx1GzodCXquQdajwhOGq2Hw2EMDg7C5XLJCk5WVlbg8/nEiTHKdpGdXB0ppdo9SzXLM+lNSTqdRjabxeuvv45EIoHvfve78Hg84t/GzCi2/v5+eDweWK1W2Gw2XLt2DZFIZEt4KQK7TBTNoidaUxNDpdOGmbtaI+KUz+cxNzeHlZUV9PX14dKlS3j11VdNLRBmBIo4V6RSKYyNjYkT7+t5TmBjYfD7/cjn83Wfi6onwlSLfsp5K5JFplgs1q3vbKtGikaeX6vgJJfLyZwvSLuIWmHPTqMeLRnSmxIAWF1dxfnz55HNZuF0OvH9738fgUAAf/iHf4hwOIyRkRH8+7//u6Fz/OxnP5NNBLp8+fKW8FIEdpkoVhspqiEVQ9Ljp+a0YTYFqudY0ne5vLyMQ4cO1aS61ogYS5v+R0ZGsLS0pCulo8RooQ0Zkr62tobDhw+LE4DqidnqU61ye5K+SiaTyOVyuHHjhlhBaLRSU4vtFilKqZRhkfb8SRdZpfnw3Nyc6h7adq8kbpQ1FQA4nU685z3vwcsvv4wPfehDeOSRRwBsVKNWy1bxUgR2mSiaRU2Y9IhhueP1Uk6cCoUC5ufnZUMIapU20hO1SfscSYTGMAyWl5dN7w3qOU56E9Df348DBw6gra2tYYUYtVhIBUFANBpFMplEoVBAPp8XU4GFQgHJZBKRSAT5fF5mbUZSTna7HXa7HYIgoKmpCT09PZviKtAIUTQ74lDLfJhElQsLC2BZFlevXi1pw6lFwVQjBLfeY97U/r7KPUWj2xQMw+A3fuM3wDAMPv3pT+NTn/rUlvFSBKgo6sJut4sLkzQqGhoaKiuGhFo34Etntfb29uLixYs1/2KUu2alc4Wyz7GavcFyI/V4nsf8/DwWFxdlEfHMzExDx7zVAq/XC6/XW/L44uJi2eMKhYKqFVkwGMTY2Bh6e3tL0oSNbMmox/PXcn9QuodG3stTp05pprar8VNsxOSXek+00bKNqma//pe//CV6e3sRDofx/ve/v8SBaDO9FIFdJorVpE9TqRTefPNNQ2IoPd5spChNKRaLRSwuLmJxcVFzVqsatbJIkjpXDAwMaLZ2VFswo0QQBCwtLWF2dhY9PT24ePGi7HU3eiB4teciAxxqzdraGmKxmJgm9Hg8YoqwXu/PVo0U9T63xWLRTG3n83nZcAfip9jc3CyrgtVqjm/UXNJ6nqPWtlEA0NvbCwDYu3cvHnroIVy5cmXLeCkCu0wUzUAqCLPZLI4dO2ZIDAnVzk8VBAHz8/NYWFgQR6LpvTskvohGvzjS11goFDA7OysW8FRyrqhGFKWLN2n49/l86OzsFBvv1a51O+0NzczM1Px6bTYbjh07Ju7lchyHZDIpeiuSHstamxBvt0hRSqVIzm63V5wDu7KygkwmU9Ic73a7GyaK9YxGay2KpBLb4/EgnU7jv//7v/HXf/3XePDBB7eElyKwy0TRyJc3lUrB5/Mhl8uhu7sbHMeZrmo0K4pkPmg6nQbHcbKRaHoh02XMfDkFQUAgEMDS0pKhAh6pQbHRayXHra+vY2ZmBh6PB2fOnNFs+Fcep5dqBjlUI2jJZBJLS0uyxw4cOIDm5mZxj9Rms4l7hzabTbxZ4DgOV69eVX3eixcvwul0Ip1Oo7m5GQ6HA52dnejs7MSePXuwtraGoaEhmTMDKT6pZlZpI6tPa42Z9GalObDS+aQcx6FQKMDn84mCWc0c2HLXVC/URDEej5sWxZWVFTz00EPic//BH/wBPvCBD+DcuXNbwksR2GWiqAepGJLWilgshmAwaPo5jY5q43keoVAIs7Oz2Lt3L1paWjA8PGz63Eb3M8nQ7HQ6DQCG9yyJEBuFYRhkMhlcvXoVDocDx48fR0tLi67jtkv6dGZmRvb/brcbJ06cgCAISCQSOHz4sOaxd+7cKXls3759OH78OCwWC1555RXE43E0NTXhXe96l+xGgmEYVRPicrNKpZGPx+NRvSHbrpNmgNru+ak1x8fjcSwuLqKtrU3zvSX/bVTPqlHU9iw5jit7k1qOwcFBvPHGGyWPd3Z2bgkvRWAXiqLWoqYmhoR620cRyNDq2dlZWbpwbW2tIc3/UueKffv2oaWlBQMDA4bPaUYUyfufzWZx+vRpQ7MVG+mnWA2xWKykqm50dFSX0LIsi9nZWdlj/f394r5uIBAQS+NzuRy8Xi+OHTsGoHwVpNasUmVLw+zsrGiOSxZzj8dT1z2/elPvEWw8z8sidunj0jmwCwsL4DhObBeRCqYe15Z6omYwvNPZdaKopJwYEuotilLniI6OjpJh2UTY6iWK0vNLxTgcDpuKBCpVkUrJZDLwer1gWRb79+9HMpk0PGzYaPSWy+XEaSgej8fQKDYzkSIZBzc1NSV7vLW1VXdBwfT0tOy8TqdTFFSgNAJdWFhQdT/Ri1pLA4kqyV7l2tqa2FLC87wsBVtrX8V6UO/9OK3n15oDSyL2dDqtOgmJiCVxbWnEjaBa+nSzq0Prza4TRbKo6RFDgt1u173Iq6ElikrnCK29M3K8mYWmnCgqHSSUYkwiMKN303oiN47j4PP5xLaO7u5uJJNJJBIJQ+cCNv6meiLTQqGAQCCA1dVV9PX1wWq1yhYgZdM8GcWmpJwokkKMeDwu+6f295eKWrnnjUajCIVCssdcLhf8fj+y2SxWV1dVXz/Z+6lVilMaVZJG+Wg0inA4jP3795f4/SnHr7lcri21mG4lL0VAe+Qay7IlrhdWq1UcfRePxw23i+hFaTBMItqdzK4TxVQqhZmZGXAcJ7ZWVKLa9JxSmKRi1NraWtY5Qu34as5NIM4V5RwkzIpiuUIbIkzhcLikraNWVatKSH9jMBjE4cOHcfHiRfA8j0KhIPv7K9Na8/PzskIUj8eDfD5fYnkFbESf09PTWF1dVe0jVNLR0SG7EdMSC0EQSiJMYMMBJBKJlD1HI9wMBEFQ3U9TG7/GsqxqleZmLbLbwUtROl9XmlUoFAqIRCJIJBJYXl5GKpVCsViE0+mUTeyp9kZE2ae4070UgV0oitFoFIcOHdIlhoRq725JpCcIAtbW1uDz+eB2u3XbGdVSFKPRKLxeLxwOB44dO1a2kIUca6biVSluxWIR8/PzZStZa20dJTUWVrayqJ1HLa0lLURJJpOIRqOidZR0YQ8EAlheXtZ9vXrsuwAgHA4jGo3qel4lykKbeqAVhZYbv6bmIi9dzMv1/tWS7SCKWthsNvHm4siRIwBK20XC4XDJjYi0sEcPyvTpTnfIAHahKB4+fNi0wJjFarUim83iypUrcLlcumyUpNRiTFwikcDMzIy4ICvHX6lRTeRG3mNp8c7+/fvLVrJWYx0lPY7cfHi9XpmxsBmUhSikVWLv3r2yXrX19fWyz2Oz2bBnzx60trait7dX1/sPwPS4q9HRUbERfSvNPtWKKrV6/3K5HBYXF+tiFdWI9Gk9p80oRVdPu4jRVhylKMbjcSqKlHcwszcTiUTg9XrBcRzOnDmjq8VASTWRYj6fRzAYRFNTk2E7JbPnJaK4vLwMv9+Prq4uXT2WtfBTjMfjmJ6eRlNTU9lIvJoohKQMpY4Nra2tuHHjRsnvOp1OtLS0oK2tTVbUoxczg9X379+PoaEhXb/77LNv4Otffw6hUBz797fis599H377t0/qPlctBFdrMc/n87h+/ToYhpFZRdUqqtzOkSKgXgSjhtaNiLSwZ35+XmzBku6rcxxXkj6lorjDqKYaz8icQZKmtNvtmJiYwJtvvmlKEIGNKMOoOLEsC6/Xi3g8jo6ODhw9etTwec2IlCAIiMfjCIfDsFqtJcU7lc5n1oEil8vh5s2bKBaLuiNhM2h9fnp6ejA+Po5f//rXsscLhYLo3kGGUJOiHqfTCavVCkEQUCgUkEqlRHsewsDAAARBwMrKClKplK5rVL6Hajdzzz77Bv7X//oJYrGM+NjSUhyf//wzAGBIGOuV5iRRORkLBmhHlVIDYjLertJ3tdGR3FZ6fq1WHOW+eiqVwhtvvIFCoYCnnnpK/F6xLGso26W87rNnz6K3txc//vGPEQgE8Oijj2J9fR2Tk5P43ve+t6k2X7tOFM1CKlArfYlisRi8Xi+sVmvNFmer1ao7fZrNZuHz+ZBMJjE8PCym+cye14gYx2IxzMzMiFHUxMSEofOZEeFsNou5uTkkEgmcPHnS0F6xGcq1ZPT396NYLGJ6elp8rFAo4MaNGzh27BgEQQDP8ygWi0gmk6op11/+8pcYGBgQo8rm5maMjo5idHQUgiBgZmam7NxUt9stG/Tw7LNv4Ktf/W+srCRhtTIoFgW0tTUjmcyhWCx9r7PZPP7hH57DyZPtsFqt6OnpKbvwNnrijFZUKTUglhaeKN0vSDsDsLE413PxrXckWg+HDOW+eiQSwblz58CyLBKJBJ555hl4vV488MADyGazOHXqFL773e8aOsc3vvENjI+Pi5Xmn/vc5/CXf/mXePTRR/GZz3wGTzzxBP7kT/6kpq/LCLtOFOvhqQhspBW8Xi8YhsHo6KhqhZbZBcRqtVZsCZE6VwwODmJiYgIMw2Btbc106lVvE34ymRT75MbGxmC323H79m3D5zNSaCNtr9i3bx+cTmfdBVEPQ0ND4Hle5oCRz+dVU6tqkL2edDoti4LIQrVv3z6Ew2Ekk0nZcWfOnClpe3j22Tfw+c8/g2x247NTLG6IuTQ6VCMUiuPtt98GsJHxIEMA1Kjnwm/k+6JlQEzcL5TtDG63G9lsVmznqUdEt5UjRSOQCtgPfvCDmJmZwfve9z584hOfQLFYLBlZWInFxUX853/+J/7qr/4KX//61yEIAl544QX867/+K4ANL8W/+Zu/oaK4HdASRSIIgiBgeHhYM99uNP2qPDabzar+LJ/PIxAIYG1tTTbhRHleM1SaYUpStMRcmDjNcxxnei+yUvpUrb0ikUhUNYbPCHqEe3h4GMViEYFAwNQ5SLsQgaRWiVuD2kK4vr6OvXv3yv72X//6c6IgGqG9/Z3oSavtg/ydttpsUilSB3llO0MqlUIgEEAsFsPq6ip4ni8bVZqhEaJodtyaWRKJBMbHxwFsrA+HDh0ydPxf/MVf4O///u/Fm7r19XW0tbWJ6+LBgwcb9l3WgoqiTpSimEwm4fP5UCgUMDw8XLEfrBpRVBNkqdFuX18fLl68qLqA1KPHkUyEicfjGB4eRldXV038FMuhbK+QVrFutdmnDMPgyJEj4Hkec3Nzmr9HpsZYLBaZ+KysrMhE0Wazlbg13Lp1CwsLC+L/z83NIRqNwuVyianXUMi4I7rFAvzmb74zfzWdTiObzcLhcCAWi2FtbQ2rq6tIJBJwuVw4dOhQ3foM6xWFkveT9P6RAQflokrpXqVeoav3GLl6i66WwbDZ3tcf//jH2Lt3LyYnJ/Hiiy/W4Arrw64TxWrTp2QSjpHmf+nxeotOpEjFifT7BYNBXc4VZodzk2Ol4iaNSgcHBzE+Pl5TP0U1pO0V7e3tqu0VW3H2KcMwGB8fh8PhwOLiomxsGllgSdVkNpvFz372M/HYeDyOlZWVkshPEATEYjFYrVaMjo4iHA7LBgUUi0UcOXJEXNy7ulxYXWUNXXdzsw1nz3bLHnvppZdgsVhKbsxYlsXS0hL6+/sNnUMvjfBqJN8dragyn8+XHb1G/jU1NZVcayPGyNW7ulX5/NU4ZLz88st45pln8JOf/ATZbBaJRAKPPfYYYrGYWEm7uLgoK6zaDHadKJqF53nMzs7CarVWHAunRjW9huRY4ql44MABXLp0SdcXwmazmRYMIsbFYhFzc3MIhUJiyrLcl71WkRsp3NHTXrGVIkXp7w4PD2s6nESjUSwuLqKlpQWtra3iQG8AeP3119Hd3Y0zZ84A2HgvXnvtNfHnLS0tGBwclFW7ktL6sbExdHZ24lOfuoSvfvVn4Dj9f/90uvQzyvO85meIbB+EQiFMTEyYrrBWoxF9hJWe3263o729XdwaAN4ZvZZKpRCPxxEMBpHL5WCz2WRCWY9CGCl6WzJq+fzVRIpf/vKX8eUvfxkA8OKLL+JrX/sa/uVf/gW/93u/hx/84Ad49NFHN91LEdiFomj0zpNlWfj9fkSjUbS1teHYsWOmi2XMiCIZCRcOh+FyuQwZDAPVRYoMwyASiWBxcREHDhzQbSFV7d19Op3G9PQ0eJ7XVcHbSD9FoDa9eRzH4cqVK+J1q4nJ6uoqXnnlFbAsW/LZSafTWFtbw8GDB7G4uCg+HggEsHfvXiwtLaG7m8XDDw/iP/9zHtEoB4sF4PmNPUOWLSCXK33PpPuJeslms8hms7h9+zbOnTsn+1kikcAbb7yBXC6Hrq4uDA4O6h4T1shI0QjS0WtkQALwTlRJjMnT6TSuXbumK6o0QyP2LGvppajFV77yFTz66KP4/Oc/j9OnT+OTn/xkTZ/fKLtOFAF9d/uZTAZ+vx/JZBJDQ0PYu3cvYrFYVelXI+JEetPIsPA9e/aY8lQ00+NI9u+8Xi+cTqdhITZLNptFJpPBrVu3MDIyojs1vVUjxXLE43GZkJPGaSXlBqTbbDaMjY1hbW1NVogljSjPnu0uSYcCwLVrq/i//9ePfP6da7DbLbL9RPJZJ6/XarWira1Nc3qP8nOWy+Vw7do1McUbCoUQCoXgcrkwPDxcMU22FSJFIyijylQqhcnJSTGdrRVVejweuFwuwwK3GcMBUqlUTWaf3nvvvbj33nsBbHgsXrlypernrBW7UhTLkc1m4ff7EY/HMTQ0JLY2RKPRhnkqSoeFT05Owmaz4dq1a6bOa6TNQTkebWRkBMlksu6CKN2rtNlsOH/+vKGbD6ORIjFRTiQS4j6fXkf0au/ws9ksotFoidWTUaxWK0ZGRmC323H8+HFcvXrV0PFEKH/yk0VEo1m0tTnwm795WCaggiDg4MGDYr8kx3Elno5SpH2DgiCIEaKUzulpdL/9NoJHj+LN0VFMTk6iu7tb9X3dqpGiEcjcUWVUKTV2JgMdAJREleXMhzcjPUumOe1kqCjehfT5RaNR1SKSRhgNazlXCIJQVQpUD2Shbm5uFs8diURk+1y1Rq294tVXXzW8EOqN3sgNh9frRWdnJzo6OsCyLObm5sCyrJgWI0KpNhXFTKS4traGYDCIaDSKTKZ8j6AeXC4X3vve94r/b3ZAxNmz3fgf/+Oekik8UkhqtqOjA8lksmy/7MrKCmKxGFpaWsCybElE2Tk9jfd+8YuwFArgbTa89IUv4Do2Bg4MDg5i//79MpHabpGiEdRsoiqZD0vdL8w62BhhNxoMA7tUFKULG8dxCAQCWF9fL7EyklKtKFqtVk1LoUrOFfW8W5YOCp+YmJA5RFTTzlEOQRCwtLSE2dnZkiHh5G9T60gxHo9jamoKzc3NOHPmDBwOBziOU3VwSCaTsqko5O6dWEcZWRyWl5d1N+7rhbgiABtp15///Oemn6ucIEqpZFMFoOINVPfbb8NSKMDC80ChgO6338b66ChSqRTefPNNeL1enD59WkzP1TtSrHd1qFEqmQ8ro8pMJoP5+XnZXmUt0RLFreSJWQ92pSgCEFNBpOldafiqpBZGw8p9IyJIDMPUdV6nGul0WhxUPjo6qjp0oNpWB+Wipqe9gqR7jdwBl4veWJbF9PQ0isUixsfHxfdY7XVpDU6W3r1HIhHkcjkx9apsr1CytrZW8fpJSnH16FGsj45WfK1kwlEymZTtH5ZDT6QnpdqbQDVWjx4Fb7MBdyPFVcU8XpZlcfv2bVy8eBFA/SM5QRDq9vw8z9dMPLSiytdeew1NTU0y709pVGlkW0ANpcFwOp2uaXVxPfj2t7+Nb3/72wA2btL6+/tl7U562JWiuLCwgNnZ2bJN70qqjZqki0wqlYLX60WhUDDsXFEt0tmoIyMjZVtLqnnNyogvFothenpalp5Vw0xaSE0UOY6Dz+dDLBbD6Oio4RYa6XNL94TW19dFT85kMin61knHsUndMLq6umRN9krUUorlhFEQBFP7y729vThw4ADC4TAWFxexurpa9veHh4fR1NSE+fl5mZ8jwzAYGRnBoUOH8Pzzzxu6Bv7CBdz42tfQ9Ktfad4AJJNJ3L5921RUbpR6R6H1TG1aLBbYbDb09PTIHuc4Tvxckm0BACXGznpmvm5Hg+HPfOYz+MxnPoN8Po/7778fn/3sZw0/x64Uxa6urpL9i0rUwmg4m83izTffRDabxfDwsOFZnWbTSQzDIJfLYXZ2Fuvr67IConJUEymSEXEsy4pj8KSRWi3PKRVFaU9luXQ4Oc4MUocBafqVDHdIJpOyRm+n06k5pk+ZUuwLBCpGi2axWCzo6elBT08PIpFI2SjT5/PhnnvuwYEDB5BMJrG0tIRisYhDhw6Zymjcc889aGlpwcrAAF6XFJwoGRsbQ1tbmzgwnexNSu2iPB5P1SPY6k0jolw1HA4HOjs7S5wvzPgpKlsy6tGOUS8ee+wx3H///fjt3/5tw8fuSlF0uVwNNRrOZrMIBAKIRCI4ceJEyVg0PRCRMXr3WSgUwHEcrl69qitNrDyn2fdJEATcvn0bmUwGo6OjsubncpixjyKiGAwGMTs7a6in0ijlUrVq49h4nsfKygpu3rypeowypej8wAewx+Mp24phhnQ6jVwuJ/6rNJc1n8/jzp07OHHiBDwej2wfE9goTGMO3YLlxIuAKwGwe8C/eS+weFz2/thsNly6dElMu5X7PHV3d+PgwYNgGEaW5u7t7dW0iyJROYnMt0pl5FYaBm6xWMTiMYLUT5GIpTKqJBZnhEQiseUjRQB48sknMTc3h29961umjt+Volht47be4zmOg9/vRyQSweHDh8FxnKxs3Qik+d/I3MWFhQUsLCzAYrHgzJkzhv3PzIgiaa9IJpM4cOAAjh8/Xtf2CmBj345lWSSTSdV9ylpitPqUYZiyTgLro6N46QtfeGdPURDQVYfr9/v98Pv9mj8/duwYOI6T2V4Fg0H09vbKoo54PI65uTksWf4/WM79BIzt7r5jSwKWcz8BDwALG64aDMPg9OnTYuFIOp3GrVu3VM9vs9lw9OjRkj1ohmHK2kUpI3NBEGRtDR6PZ1O8+baSKKqh5adYLBbFPXSy1ZLJZPDFL34RPT09KBaLuH37NkZHRw21amWzWdxzzz3I5XIoFAr4yNbDnAsAACAASURBVEc+gr/927+tuZfi9evX8bWvfQ2/+MUvTEfqu1IUzaJ3v4sIw+rqKgYGBsTh0NLJI0bRK1DSys59+/bh4sWLuHXrlqm9GSMCQGayLi0toa+vT2x5MNNeoVcUE4kEpqam4HA40NzcjLGxMUPnqnd1I7DRsB4OhzV/3traivXRUVnKdG1traHzXDs6OnDw4EHxsyP137x16xbe/e53Y21tDbOzs+L+ovU3X3xHEO/C2AqwXHgWPABh4RiOHj0qppeLxaJoAq3G2NhYyT5zJRNgrchcWhhFClBIqpBElmYNcvXSCC/FevQPS2f0Li0t4dSpU2AYBt/5znfw5JNP4tatW/jiF7+I6elpdHR04Kc//amu521qasILL7wAt9uNfD6P97znPfjgBz+Ir3/96zX1UvzWt76FSCSC++67DwBw9uxZ/NM//ZOh59iVomh2IbTb7WWjNalzxeHDh2XDuqtd5CqJoiAICIfD8Pl86OjokEVMZtOget4nrfaKSCRi6vXqeZ8ymQxmZmaQy+XEytlf/epXhs5DIhCjGLlRiEajeOutt8r+Tjwex4EDB0qiSa33oKPQgd5CLxyCAxzDIWgLImKr3C6hhdVqlUVoLS0tMlFkWRbPPfdc6Wt2qad3GYsAy8VnwFx8BrPCPthyn8L+wm9gampKlhKWVtwy73oXDh48WPJcZm5a1NoapKnCZDKJ1dVVsCwLlmVx584d2b5arYRmq0eKeiDCzjAM+vv70dvbi6GhIfzpn/4pABiqxmcYRvyb5PN55PN5MAxTcy9Fo4bHauxKUTSLzWZDPp8v6QcqFotYWFjA4uIiDh48qLqfVYtCHa0SeWnT/5kzZ0o81urRbyidvNPR0YHz58/LLITM3gSUOy6fz4vpaDXLqkZQSRRJoc/S0lKJEbAWS0tLaG5urtjY31HoQF++D1ZsfLaahCb05fsAwJQwNjU14ezZs+JiNTMzg5WVlZLfU3297B6gRUMY7/5JsswKbjv/HvGVOObm3hEpZcVt4j/+Q/XvWKtoSy1VWCwWcePGDbGQaGVlBT6fD8ViUfRVJFGlmVmlO0EUAfm6FY/HZdWuRi3DisUiJicn4fV68Wd/9mcYGhracl6KABVFQyiFiaREFxYWxCipXiPR1IQtHo9jZmYGNptNtem/3LHVoDb9Ru2cZkVRuQhLJ9/09fUZKhZqNG+99RZCoZDh4/RMuukt9IqCSLDCit5CrylRzOVyougEg0H4fD7Zz8v1T/Jv3ivfU9SAZ3KY9zwJ4M/Fx6QVt0yxiJarV8Hdc0/JsfU2MNbqSyWzSol5NZlVqizqKSfYm2HrVG8SiURV7WNWqxU3b95ELBbDQw89hDt37tTw6mrHrhTFaj0VpSnDvXv3lkRJ5TD7RZcO9jba52hWoJSkUind7RVm3TmkkaIgCFheXobf7y8xFt4sykWKxWIRy8vLqj8z0qAvRRpBOgT1AgStxythsVhEE2FlAUyl/klh4Rh4AJYLz4KxVEgnK1KtpOKWKRYBhwOF97xH9bB67stp1QaU81Uk/X/SqTLS/j+PxyOuA/XeUzRrWK4Xtc94rVoy2tracN999+GVV17Zcl6KwC4VRbNYrVasra2Jm8xGKx2rmVdotVrBsixu3bqFdDptyEXCrG0VIZPJwOfzIZ1O626vMJs+JYU26+vrmJmZQWtra90rSo1QThSVJewEow36UqQRJMdwaBJKR3lxDKfz6t+hubkZvb29iEQimJ+fl/2tGIbB+MqK6kg2KaIwVooYWXkZf/P99yN99iycr76KwnveA/7CBdXD6hkpGh3xZrfbVafKSPv/ZmdnRSNxYkTMsqzmtKNqr7/R6dlqIsXV1VXY7Xa0tbUhk8ngpz/9KT73uc/hvvvu21JeisAuFUWjH1AyniwYDIqzM5X7dnogaUyjH2aO47C6uop0Oo2JiQlNV4Fy5zUjUPl8HrlcDjdu3MDw8LCh85oVxXw+j5mZGbhcLhw/fnxLjpXSEkVpgQqwsWfHcZzmzE+jBG1B2Z4iABRRRNBmfB9maGgIhw4dwp07d0r+TgzD4Nf79mFvmZFswN2/MRFG0q+Id/YUAUAo2Db6F+9y/vx5cV+Pe/e7y15jvSPFap9bq/8vl8vB5/Mhn8+LLQ1k2hGJKKvtqSwWi3W9UdQyGNbbb6wkFArh4x//OIrFIniex8MPP4zf+q3fwsTExJbyUgR2qSgC+qsISRGLy+VCX99GUYMZQQTeSb/q/TAXCgXMzs5iZWUF7e3t6OzslKV09GKxWAxFitL2CpvNhrNnzxr+AhoVxWw2C6/Xi/X1dfT29pryjtQbWZD09/LysswVQ09BRbmfK0XxwIEDGBgYQIRlIfzHf4AvIzB6IPuGtag+JQMd1Gaz8jxf2j+pIuI2mw0cx0FYOIYi6U1UaegX7v7svvvuM/Tdqfds0no8N8MwcDqdaG5uxp49e8SWFNJTmUqlxLYXEk1Ki3rKWUVJUc4lrTVqolhN+vTEiROqg/G3mpcisItFsRKxWAwzMzOw2+04evQo3G43wuFwVVZKeocsS6tZSWvH2tqa6XOTxasSZCrM3NycOBXmxo0bpnoc9UanhUJB7OkcHh5GS0uLqWn/pECn0oKyvr6O6elptLW1YWBgAJlMRhxjlsvlYLfbZYO+XS6X7DnL3UwpRbGlpQXNN2+iLxAA99WvIvjGG5gbGNCMEqX7jgAwmUwiPDEB/759YjtDxBapqgWDQCoty6Hsn1Si9pkSFo6h7XkHut9+GzmPBy3ZWeC9feh+8EHDN5PSRnIyYq5WNGLPTxoJavVUkqIepVWUtKhHbaj3ZhgM53I5zZnFOwkqigoSiQS8Xi8AlDhX1NtTked5LC0tYW5uTiwsIR/MaipIKxW9lGuvMHveSpEimbhDhJ8MZlfub+mF7EVqLXSpVArT09NgGAYnT54UTXPdbrdsSgoZqJxMJsVJOdJB3zabTfP6lKLYMTUF10c/CnAcBLsdc3/1V2UFUdx3tFjAAGB4HnscDvQ++yyyDzyA1157reQcZqnXmEPxdeTzYAQBYBgITz+NoMeD4j33wOl0ip9xYGNQ+aFDh1AsFsX3nfxLJBLi7zU1NW1rUVSjnAExKeoh49fIYHqphdlmeClu1arvWrJrRVF5x08qOvP5PEZGRlTTBLUQRbXFSBAErKyswO/3o7OzU7WatZpzlxM20l7hcrlw+vTpkrv5avoN1Zp7pa917969uHDhguzLZ2SijZ7r5DgOXq8XyWRSViRkZKCydJwYccm4evVqiUUPmR1J2PP66wDHgSkWIQhC2b1E2b6jIIARBDAAeI7D2v/7f2g/fx7nzp3Da6+9VnKerYT4OgQBArAhjIUCCs89h7dU5mneuXNHV2k+mdlaK8/ArWxgrPYZJF6fxJVlfX0diUQCLpdLFlXWalC6WqTYiAlQW4FdK4oElmXh8/nAsmzFis5aGA0rjyf+gh6Pp2wBTzWRotqxJHICUGIuXIvzqolUNBrF9PQ03G43JicnVRc4LTGthPImR+qWMTg4iPHx8apacUjqq6urC16vF0ePHpUNqE4kErLzW61W5N/9bjQ7HOBzuYp7ibLB4JJIkbfZ4O3tBf/aazh9+jTOnTuHl156ydTrqBXl2kvE13E3UuQZBrzVCv/hw1VHuclkctuIYq3Tm8qeykwmgyNHjkAQBPGGLRQKIZvNwmazlRT1GH2tyj7IzeiL3Cx2rSjmcjnMzMwgHo/rno5Sy/Qp2bN0OBw4ceJExU3zWokiKWjR215RC1GU9jeS/Vk9xxmBiKIgCAiFQggEAjhw4IBs1F4tsVgsskUqn8/jueeeE38uCALecLng/PKX4Xz11Yr9idLClpzHg/03b8IZiSBw//0bx0WjeOGFF9DS0oIDBw5gdXVV9eZhbGwMPM8jHo8jEolUZYxNkEbBldpLlK+jKZk03JspxWKxiBNmgNpFK9t94gzpU7TZbCWD0vP5vHjDtri4KBuULo0qyxXPFQoF2f7hdnHIqAW7VhRXV1fR0dFhKIKohSgmEgncuHEDPM+X7FlWOrYaUczn85iamsL6+rqh9opq06dvv/02UqmU7r5KM9ZR5LhoNIq5ubm69jbqrVom/aj8hQu4s39/2d/t7e1Ff38/Xr77//f+7d/Ccvdz1jY3h8Thw6KopNNpsXFc+nmURm+p48fFSuXW27dhf/llLI+NYWVw0NRrlqZr9bSXVCrQ0aKlpQVtbW3weDzYs2cPAoEAjhw5ItpGLS0twe/3i5GQ1MzZ6I3PVk6f6qGc6NrtdrS3t8tueMmgdLIFQDwVm5qaSjwVGYYpSZ9uJy/Fatm1okg2941gdsEGNtIdi4uLYFkWJ0+eNNzvY7YBv1gsIhgMIhKJYN++fYZHpJmJFAuFApaXlxEOhzExMaHL0JhgRoTT6bQ4GaPevY1GBoLrneMoCALu3LkDp9Mpig55tyzFouZepFQQldHb0ugoOqenceHu40MGhwYQlKlSpf+j2fYSNfr7+3H48GHx/+fm5tDc3AyPx6MaCSWTyZLpMtLK4XITXyo5cFRLI2aTGrVk0xqUTop6yKB0aQuXzWaD1WpFPB433bi/sLCAj33sY1hZWQHDMPjUpz6Fxx57DJFIBI888ghmZ2fR39+Pp59+2nQfZC3ZtaLYqA3jXC4Hv9+PWCyGvXv3wuPxmPrDGxVkaXtFT08PPB6PqhOBnvPqFSme5xEMBjE/P4/Ozk50d3fLBgjrwUihDcdx8Pl8iMfjcLvdGBsb27Rmf7M3SwAQDofR2tqKbDYrig6JFHmrtaLwaEVv1Q4N0EqVVupflB5vZLTd+vq6TBS1oi21SEhaiLK6ugq/3192uHe9IzmgvmtMNZ83gnRQOumnBDZutt566y1YLBYsLS3h05/+NOLxOOx2O/73//7fOHXqFE6ePIne3l5dr9Fms+Ef/uEfcObMGSSTSUxOTuL9738/nnzySTzwwAN4/PHHcfnyZVy+fBlf+cpXqn5d1bJrRbHe5PN5zM7Oip6KY2NjSCQSWFhYqOt5SXuF1+sVK1ltNpuq+4Ee9ESK0nN2dXXh/PnzyGazFd3d1dAjwjzPi04U5L016xlpFLVIMR6P49q1a6afs1AoYH19HcBG6vH1P/ojDL7wAop2OxI6bmS0ordqozotUdVKjyr7LI2OtotEIrI9QyP7h5WGe8fjcXG4t91uB8/z2LNnjyxlSNnAZrPBYrGgt7cXTU1N+PnPf44f/vCHeOWVVzAyMoJf/epXePLJJ/H000/rer79+/dj/90tBI/Hg/HxcQSDQfzoRz/Ciy++CGDDNuree++loriZmP0SVLrLlE6DkfbfAdXvSVaCVHdqWUiZoZJIxeNxTE1NlYy/4ziu5tZR0gHhUu9GcpxRUayVn6LX69U1HKESndPT6HvpJQy8+CKYYhGMIKBragr9L76I2Xvvxdx736sqLFrRm5GoTo2cxwOBYTYqSCuIqjKqnL3nHsNRKvE8lO6zVyNWWsO9OY7D9PS0ODiCpAylEaXb7d6y1ZaNaI1Q7immUikMDQ3h4YcfxsMPP2z6eWdnZ3Hjxg1cuHABKysrolj29PSYvnGvNbtWFM2iNapNmjok02CUX6p6+BoC8sb0StWdRrFarcjlciWPsyyL6elpFItFVceMagp01MSNCL7H41EtojHb32gG5fXVYm9K2fTOABAAWAQBTKGAoeeew+DPfobr//N/InH4sKoAagmmmaKXzulpnH7qKTA8D1gsuPHxj5d9HmVUCcBUlMqyrO7iM7M4HA4xZUhSsIVCAel0WmxtIGPYlC4YW2EwfSP2K5UD08n2TzWkUil8+MMfxj/+4z+WVLKaNf2uB7tWFKu1jyJfDhK9BAIBdHd3l7WRqjZSVEapmUwGXq8XmUwGIyMjddmkVoob2ceLxWIYHR2VNRiXO87s+ViWxdTUVMV2jkaJotrn5siRI6L/HkkhKtsRhoaGSvwKpSib3ons8sBGzyIAFIuYfOIJ8BYLLHd7GNXSkmZtqrSuhxcENFUwTJb2JwoMg+jAAObe+17N61C7RtIP2giU2R6bzYbW1lZZMQmp2EylUqKbiLRik0SVSheMeqfxN8NgOJlMmppHTMjn8/jwhz+Mj370o/jd3/1dAMC+ffsQCoWwf/9+hEKhqkW3VuxaUTSL3W5HPp8XnTO8Xi9aW1s1m9GlVBspkuNJ2icSiWBoaMiwa4aZc0qb4ck+XrlzmnXmIOKmV3wJZtKnZtJQaulTp9OJd73rXVh95hn0SaI9gWHA2+146QtfgJYcSkVUsNk2hoZbLJi9915EBwbQHghg8IUXAJ7fEEaeh0UQYLk7KUaZlqzGpkqK2n4kwzBobW1FV1cXurq60NraiitXrqCnpwdr3d248YlPYPKJJ8DwPE4/9RRe+sIXcOehh0per5givivs1y5fhuXd78a+ffvEYQgsyyISieAXv/gFHA4HxsbGqjK4VaLHOkpasUkKxogLBql+XVlZEV0wiEjWe4+yUaIoJR6Pm77pFgQBn/zkJzE+Po7Pfvaz4uMPPvggnnrqKTz++ONbxjYKoKJoGJvNhmg0iqmpKTidTpw8eVL3tPpqvygWi0Ucnt3f32+4vcKsCMTjcbz66quaaWGtazV7A5BKpXD16lX09/dXFF/pdTai0AZQjwQYhsFBrxfWu3uBJPVZbj9NKWA3PvEJeHI5pJxO2ONxJA4fRuB970N0YABn/vmfwfA8BKt1I4q8KyjKtGStbKqU+5G4dAln+vtLbsCsViv6+/vR19eH5f/zfwA9gi1JEaNQgPvaNdw5dAizs7Ml10Gm4Ny4cQP33nuv4dehhVlfU+KC4XQ6ZRWbag3zV69elbWJuN1u3Wbk5aj3dBm1daIaL8WXX34Z3/ve93D8+HGcOnUKAPClL30Jjz/+OB5++GE88cQT6Ovr01248/+z9+XhbdVnuq92L5JleY3j3ZZkObuTOE5Yy7RAB5j2FjoZ2mee0qedaWcG5rZPp0x5ysC0U2amM3Bpy4XShc5tb3tbynSgUOjQFgqkQBKThSSOF8mWd0u2rH1fzjn3D+f340g6ko7WBOL3efJAHEvnHB2d3/f7vu/93rfcuGyDYiEByu/3w+FwQC6XY8eOHWXvfRAQoXCPxwONRlOQAz3J+PLpfzkcDkxOTkIqlWYtCwsh3yBFNFGnp6fBsiyuvPLKvK4xn3ItwzCwWq2w2WyQy+VJ9lG5CBbZvjeJq66CUqkEF41uBK0cBJXUAKb0+bC0fXtapjf7gQ/A19WF7gvybr/aVYWfN0zDoYyhXv4MbgjcgCH1EABxjFOx5dWkfqTTSRmywIZ7el9fH1iWBcuymJycRKC/H31Zjp1WIr7w+awODua8d6UmqJV6JIM/JtLU1ETL/aRPKWZMRCzyfY7zRSaD4UIzxauuuirjWvDyyy8X9J7lxGUbFAHxC3coFILFYkEsFkNTUxPq6uoqEhA5jsPa2hpmZmbQ1NSElpYWbNmypaBdYj7lTJ/Ph6mpKahUKphMJthstrx3uPk85B6PB1NTU1Cr1dizZw8mJibyvkYx95L4KM7NzaG9vR0HDhxIUvogBAuO49IGwSUSCdbW1mC32+HxeDAzM4Pm5mb6b5FIBAs6HQL33w/dmTOiJM6EAlj3a69BFovRLIqfbfUcOYLnd6jwRGMjIsqNBd3DePCM+xkAwJB6KD3DA2B65hl6Ho1mM973wAOQxuPg5HKs3nQTxvftyzub9Hg8OHXqFGQyGV5++eWNoJWD7Zqq75qNUZsKk8mU1/nlQiUMjDOZEGcaE8llF0VwMWyjihnef7fhsg6KuRCJRGC1WuHz+aDX69HY2Ijl5eWSkmUyQWi8YmJioqSi4KkIh8OwWCyIRqMwGo3QarUIBoNlI7DwGayEREOcufNFrkyRlLz5EnCxWAwSiSRtvo1lWbrDX1pagtPpRCQSSQq6ZrMZZrMZSqUSMpkM4XB44x/6+tLk1Orq6mip22KxwOPxABAoUQLouTC3xQHgpFL6c5Jl/a8PttKASBDn4ng58DKu3nL1ht3QhQxPqL/Yu7CwoZjDskAshi3PPovmX/+64N5j6ndKiO2q0+nAsiw8g4MFjYjs2bOHUvdLhXIGxWxBK9uYCFGWIZZl/J4mv4pR7vKpUFD0+/2bQfFyRjweh9VqhdPpTHNYkMvlgiMKYpFppIPA7/fDbDZDKpWmsS3LZR9FrtflcqWJo5djjCQej2NmZgZutxsGgyGpN1MoizRTpkgCL8uy2LFjR85xlWg0ioWFBbhcLgSDwZz3Wsx8os/nw8mTJ9Hb24sDBw5gdHQ0KTCS4GB65hlILxBqWIkEs+97H9TXXw+dQgHPnj1gn34atnrhR9YZc6a5UAj1F2e3b0e3XJ5kUSVjGLROThbMVM0GImAeDAY39HsPHIDj0KGNe+x2Z32tTCbDlVdeWRaVoosVFDMhk10U6VPyx0SIEo3T6SzLmIhQUGQYpiT90HcDLuugmM1uqLu7GwaDIe3BIezTQpEpKPLHK4xGoyA1vdRGw3xlmO7ubkHiTqGjFULgGwt3d3djYGAg7XiFEmZS2aeJRAIzMzNwuVyi2KvAxj148803CxrET+3TpX5uHMfBarUi8uqr2LO8jCW9HjMtLUnnnFpOnb/2WrRpNBvmugYDQnv2oOX032FVHk47fr0s/fsSv/JK4Je/BBuNgpNIENVo4DQa8eo//ANMo6PY+pvfgEskAIkE/QcOQDU0hNOnT+d97dnQ3d2dJt1GsvD6+nrKNE2FUqlEXV1d2WT7yjkAX6qAK5PJBMdESG/S7XZjYWEBsVgMVVVVWcdE8kFqJkrcZy4XXNZBkYBlWSwtLWFxcRHt7e1ZiSyl8FTkB6dYLAar1Qq3253TwqpU9lF8e6VUZZhSHpOA3xsVMhbmo1ilIY7jsLS0hIWFBXR1dYlm6Pr9foyOjhYcEFPLlFtvuw2zs7NJ2Vuj2YyDF36vUS6H9gc/wHRzM/x+P1iWFVSgcY6NQaFQYMuWLZBccQU+13wf/umVryDCvXOeSk6GG7Q3JAXiRrMZW+x2xP/qr6D43/8b0gtjEsRx4w2jEb1tbXSEouqee7B12zaETSZRpr9isba2hu7ubvp3fp9t69atMBgM+O1vfyv4mTafPw/70hJkV11VUgNdgnIFxXL2/KRSKQ2WRCRdzJhIPr6KqZkiCYiXynB9uXHZB8WVlRXMzs6itbU162JNUCpPRX5W2tPTI5g1Cb220CyVBDen0wmLxSLaXqnYTNHj8cBsNqOmpqZk0nNCkEgk8Pl8OHbsGBobG0XdSwKfz4fR0dGCP1uhMuXa1VdDp9MhHA7TTUXz+fOQxuMbA/HxOIIvvADvhTm+9vZ2rK2twSfQa3v77bexf/9+NDU14SbDTZD/+gV8M/wabPUytHkYHHb3obV7CFVVVQiFQklBWiKVbsw4chykKcQdlc+3UUJlWXCxGJjf/x7+m28uyfA/wfr6elYvvunp6bSfNZrNuPKf/xmSWAzcU09h5rvfhaW/H5FIJImQotFoLknd0koQYfjvL2ZMJNVNhJ9VppZFU4NiJBJJ8lZ8r+OyDopWqxXhcDgv771SZIp2ux3j4+M5s1Kh10YikYKOm0gkYLFYUFNTk5e9UqELTjgcpiVhIRm4UiIQCGB+fh4AMDQ0JGpulOO4DfLHBRZltoCo0+kQCoUE+4t6vV6QReoU0HGMajR0hlHCcYjyPhObzUZLhlKpFHa7Pelc5372M+i8Xiiuvx4fPPgJ3PqhXwCxGBiZDK/d95dw4h3fQ36Q5gBwMtlG8FMoksYkkrRNZTK8VVMDvPKK4PB/aqDMp41w9OhR3HDDDWnfJa/XmyYar1AoMOT1QhKPb5CB4nF0Wa3YcmHzIERIkclkabql5XbAyIZKeCmK2fBl8lXkj4nMzs4ikUigqqqKfn6RSCSp9345MU+Byzwo6vX6vANcoUGRlBBXVlZQX1+f99wfUFgpMxKJYHp6Gi6XC21tbTAYDHm9Pl/wSTsqlQp79uwp20wV3zqqtbUVUqk0Z0AkpSCGYeD1enHy5Mmc91OpVMIvIHMmk8kwNzeHhEjhbZXfD04i2cgUJRLoZmeh4o1LKBQK7N69G36/HzKZjPoxNprNuPpCoMIjjyD47LMIPfccZH/4A07V1cHZ2Zl0HBKkJQwDKJWY//zn0chxCOzfD+eF708mbVPTM8+kZb1AuuMFc+AAJQsRZHo2WJbFiy++CL1ej46ODlRXV4NlWZw7dy6pV6VSqXD11VdDpdEAjz4KluMgUSqRuOqqpHuRSkhJJBK0dLi8vExL1qnzp+Wc7eODYZiyaqQWk4lmGhOJRCLw+/00WK6trcFms+EXv/gF9aYUIuDkwqc+9Sk8//zzaGlpwdjYGABcsj6KBJd1UCwEhciXuVwuWCwWqNVqdHZ2oqqqqiAmVz4BOZFIwGq1Yn19HXq9PmP5qlTg92VJL+/EiRNlGedgWRYLCwtYXl6mknNra2uCgYsPkh2SDPHEiRNpn2d3dzcikUiSYn8m9X7+BkWM8LZj+3awCgWQSICTStHz6qtpGqYnTpzA7t27sWvXLmg0GkxOTiZlfmwshrWnnsL8xz8O9c03IxGJALyskpzLkfvvx6FoFOw118BRVwdVR8dGNn3BuiyTtik/64VSidgVV6D5zTfTAuVkyrWSjd7MzExGjdfp6WlMT0+jubmZbjT4GWjPxz62Yes0MgLXf/4nAi+8gKbbbgM7MpL1cyWaqXxyGj8jWl1dxczMDBiGQU1NDdRqNRKJBGKxWFmCV7nLp6V+f4lEgurqalRXV6OlpQXxeBxbtmxBdXU1lpaW8PLLL2N6ehojIyNQKBTYtWsXvv3tb4sKkJ/85Cdx11134ROf+AT92de//vVL0keRYDMolhFC4xXLy8tF9wWzgc/w5FtXEk/EdAAAIABJREFUraysFDVKkgl8L8Xm5uakXl6h+qdijtXS0pJUes7GWuUHQwBwTD2F0NSj2AE34qjHiuRGuGV70d/fD5PJhPX19bLY2PCJNDXr6+j7/e/T5NjcbjdeffVVNDY2wmQyoaOjI60869q9G9FoFB6PJ+Mmad1gwK8AwOWC3OeDx+NJyuwyKd/wz1F5ww0w3HorHPX1YJ9+OqtKTiQSgUwmg9FohN1up/0rITgcDgDJBCVOoUDk4EGwF2b3Yvv2wd7cjIZt2wr6rDNlRESoIZFI4Pz581Tgm59RFkvouRjD9eV4f6VSiT/5kz+BSqWCTqfDQw89hFAohImJCdHHv+aaa9Lk+y5VH0WCyzoolqtBT4bgI5FI2niFXC5/Z9A7T2QLikQmzWq1CjI8yzFv6PV6YTabUVVVJUiiKUb/NBV+vx+Tk5MZjyU035gaDCUSCYKLzyM29XUosbExUcKDTu5ptDZvQb/pZkgkEjQ1NaGqqiqtf1uKz5A/WN9z5EjGQON0OvHGG29s/CW1PNvaCuTIivlIJBJppc5sXov8rPfsb36DnoEBzN9/P5rHxjKWh6PRKB1z2Lt3L/7whz/kPK+k3mciAfnrryN2ISsk83jS48chf/11JK66KmfGmAsSiQS1tbWoqqrC8vIyhoaGKHOTlA5tNhsl9JA+Zb6Enkr0FCsZdH0+H13DampqsG/fvqLe/1L1USS4rINiMRCacxIzXiGXywteWDOVT4n6jVqtzujWUcyCLmRZRZRvBgYGMpZmC2WukqxPIpEgGo3CYrEgFAphYGAgY8OfP6fI7xuS9yH3wT/xLUiRnKnLEEet5ylIJJ+hx29pacHCwkLS7w0MDGBmZqYkGXe+BsDZyrOFskXFei3Ozc0BBgOcWfrRpC9VXV0tSOJqbW1FOByGz+ejP0vtffJ7hxzHQTM2hpq77gIulDlDzz1XdGAEkoMWn7lJRhyA4gg9lXCxKCfjNpXd6vF4yka0uZR8FAku66BYrKci6QsyDIO5uTnY7Xb09vZmHa8gMk2FIDWwEXNhADnNhYsd/CcZF+lTEiWabJ9hMUExHo9jaWkJNpsN/f39aG1tzXksMmRMzlUikaQtWmzYLvh6JmxLWszq6+vTguLMzExJh5gLNQDmo1RWUaUAf56tvr4+KTvt6upCU1MTvF4vva9OoxEz3/0uumdnwVx9dVLAY1kWdadOAbEYJAwDLhZLyiSLgZhMrhhCz8WwdiolUs/f7/ejt7e3ZO9/qfooElzWQbFQkKBIGIILCwtob2/HoUOHcj5spZBqi0ajmJ6eRiAQgNFoFMXcKjYoLiwswGazobOzk/Ypxbwu36DIcRzi8ThGR0dFf6YE0WgU0WgUMpks4w5UVr0FTNiW9vM46nHixAns378f8XgcbgEJstQMsZTzfPlAJpOhpqZmg1zy/PMlsYoqFlqtNmmWbd++fThx4gS8Xi/a2trQ0NBAf0+r1WL79u00ixfqsHMX2LJQKsHFYmmZZDEotLwpltATCATAMAy0Wi0NlqUk9FRCXYb/7PDLp6XApeqjSLAZFAuATCbD6uoqVlZW0NzcnNd4RTFBkQTEEydOoL+/H9u2bROd7RYSFImRstfrRW1tbV4D8UD+QZG4ZTAMg/379+fUKSXnyLIsqqqqIJVKcerUKQCAWq2mbiYajYbufDXbPg/v2/8IjnmnX8hAsUG2cTjw+uuvQyqVoq2tDX19fbBarYLHLWeG1t/fD41GQzc+qbj66quhVCpx9uxZrBiNMKQQZqRSKerr6xGLxQRfXygaGhrAcRzkcjmVEXO73ejt7U3zWVQqlbjiiityimNnAsuyiAwNIfTccyXrKfLfu1Q9PyFCz6lTp9DT04NwOEyl2EpF6CmnPF0mFDOn+LGPfQyvvvoq1tfX0dHRga9+9auXrI8iwWUdFAv5crlcLrhcLjAMU5BCSyE9RZZlaUYqkUjyyp4I8g2KfPuohoYGdHV15c14E8s+DYfDMJvNiMfj2LZtGywWS85jpfYNFQoFtl8gq7AsS4kTKysr8Pv9Gz0qjQYazS5UG7+ExNz3wIbtiEt1WOauh1u2FwDoWIfX6816/FKY+VZXV6e5bwAbZVq1Wk2ZsKnsvWPHjkGhUGycq1BvkmXhcrnQ2NiIK664AkeOHMlL9OGqC7JqR48eTfocenp60NraSpVSXC4XAGB5eRk2my2JmEIkxQotI5LFnx0ZKUnJlI9yE2FYloVGo0kKJKUi9JT73DMZDBeaKf7sZz8T/Pml6KNIcFkHRUC8ADUZr5DJZGhpaUFzc3NBkmX59BT5IwhNTU0YGRnB6OhoQQ+F2AAViUSSmLNarRZjY2MFlV5zsU8TiQRmZ2fhcDhgMBgo0SFXhpmrbyiVSgWFlAOBAHw+HxzMDgR0/wJWu/Ee7guLez4QY+abC9lYyIFAACdOnEB9fT3a29vpID+wcY/4QS5Tb9LpdNIyntigqFQq6cxp6r1zvfAC2p1OSK+5BoqREahUKoTDYezYsSOp30YkxQjbU6yBMx/l9justHN9NkIP+dxSCT18Ug/5LMo9jiGU2Xu93pKWTy91XPZBMReEPAaJNFIhEJudEs3Q6urqtIy0kBJKrmDMD1B6vT6pHFYoYSbT6ziOw/LyMubn5wV7lJnso4RGLMR+DlKpNMk30e124/Tp04IODWKQL3u0EJCepWvXLqC/P+/Xk15wPiXUWCyGN998U/Bc9hJVnQcfRPhXvwK3Y0eSpVpqv41hGNpvEzJwzqY0QzY85QDDMGWXgRN77kqlEg0NDbTnCmQn9KhUKjAMU7bgKPS+Pp/vklKcKTcu+6CYKVMkEmIejydtvKJY/dNs4BvvCmmGkkBTKmd6vqtEJhJNoSQdoaDocrlgNpuzSt0J2S4VGgxTEYvFcOLECVr6KwalYI8CyYQdYKM0G9VoMPSjH230LJ9+uqCeJcuysFgsRZ8fOSe+qo7sD38At3171vsgk8myGjjzlWZIoCTZEcdxZc0UL6Y2ai5kI/Q4HA7EYjGcOXMmSaGnVIQeIQPjUCgkSk/4vYLLPiimIpFIYH5+no5XmEwmwVnDQgfwM4Gv42kwGDL6/5EAVUhQTMX6+josFgsaGhqykoVKkSmGQiFMTU0BQE5Bcv54BQmIqfOGhcDhcOQU/yaQyWS44YYbsLa2hpMnTxZ8zFxIUnWRSsEBG8PsEgkkLAspx11UVilBarnYvHUr6o8cwdbXXoP01ltFk2ByKc04nU7Mzc0hHA5DoVAgkUiUnMF5qQdFIZDPLZFIIJFIwGg0Jn1upSL0ZLKNerd9XsVgMyheQKqnYjYySykyRbLI8y2kMgVhoWMXs0D4/X5MTU1BoVBgz549OW1hCs0UZTIZwuEwJicn4Xa7RZv9kl5ktr4hcIG2HwhgaWkJS0tLUCgUaGtrS/NQDAaDMJvNeVlvMQyDSCSCtrY2XHfddXjjjTdEeS3m2kD09PRgbm6OZoc16+vvqLpccNAgguGQSsFyXME9y2Kxw++H6uhRzPX0pJWL4XBg2//8nxvl1P/zfxB89lng0KGCjkN6j7W1tdiyZQsAYGFhARzHoaqqKmnB55vpajQaqFSqvDdK78agSMB3yOB/bgTFEnoylWUvtQH7cmIzKAKw2+2wWq1p2p2ZUApPxXg8jrW1NczPz2Pr1q2iGaXFzBuyLIuxsbGc6jCpKCRTZC8wINfW1mA0GkX5RQKgZbOFhQU0NTWlzb8Rt5HJyck0AXCiftPQ0IDm5mbq2OF2uyGRSKjmplicPn0aV1xxBR1HOXr0aNp9VygU2Lt3LzQaDRQKBbxeL86ePYtgMChYrm5tbUWTxYIuMs4hlYKTSsECNFPEBZHw03fcAZXfX7KeZT5zlY1mMwYeeADSRALdMhkt35LX8d002FgMC//3/yLS2Ij+/v6SkFg4joNKpUJraytaW1vpzyKRCO23ET1fpVJJF3uNRpMzMyo3iaecASRXlahYQk9qUCz39VyKuOyD4uzsLILBYEZ5NCHk4yWXCo7jwDAM3nrrLTQ1NeVtIVXISEcikcDc3BxCoRAMBgNaWlry+qLnG4gdDgf1bmxvb0dHR0fO1/D7hl1dXfB6vUku4sRr0OFw5OwHHj9+nP5/IRsYfvCYamzE4OAgtFot9u/fj9HR0aQNwq5du+jiE41Gsbi4CKlUisHBQcjlcpw9e5b+rkQiwfT0NDpfeumdcQ4A1j/6I4SampJ6imICV1tbG1wulyjZuXznKpvPn6eehlKOSyvfppZT7SYTnDMziEQi2LVrV87zyQWhniLfzYG/4PMzI7vdjkgkArlcnhQo+ZlRuYPipah7KpbQE4/HoVQqqadnc3NzybxQX3zxRXzuc58DwzD4i7/4C9xzzz0led9S47IPin19fXkHmUIzRTL7F4vFsH379iSXbLHId6SDsDzb29uhVqvThqzFQCqVitoEBAIBTE1NQS6XY2hoCMFgEE6nM+c5pvYNlUolmpub0xY+p9NJ+5JiUUhATAoeANabm9HU1ISmpibs27cPZ86cQSwWQ39/P7Zs2ULJSufOnaMBc3x8XPBanU4n0N+PPl5Amb/22jRBbjFYX18XvTlrt1jymqvMNnIil8szsm/tdjtMJhOkUikNDoUEiXyCi0qlgkqlSnqeMmVGGo0GsVgMarW6LAGsEmLdhdjOCUGI0ENmYmUyGX7+85/j9ddfh91ux+HDhzE0NIShoSG8//3vz/scGIbBnXfeid/97nfo6OjA8PAwPvShD2FbgS4o5cRlHxQLKQ3kGxRTBbQXFxcLplOLzdqcTifMZjN0Oh3NRtfW1gqio8tksqxzbrFYDNPT0/D7/RgYGKAPWTgcLmrekA+VSkU93kpNcuJ/pkJD+ceOHcP1119Py3nvf//7qbIL2QgkEgnRJeZCxzlSy59iAmJfXx9YloXdZII+j7nKbOdIvvtC7Fv+d5P8P/kv/x7n+g4Wq9ySKTPy+/2Yn5+Hy+WCw+GARCJJYr1mEvkWi0p4KebiABT7/vX19WhsbMRDDz2Es2fP4rHHHsP999+P06dP4/e//z2uu+66vN93dHQUer0efX19AIDbb78dzz777GZQfK9AbGCKx+OYmZlJc80opieZq3xKFmmZTIbdu3cnUanJeee7y8vUU2RZFvPz81hZWUFfXx8GBweTFrJsc4qFjFhIpVKMjIzAarWmiXUXA0Jump6ezpghzczM0AdYJpMhEAhQx4zOzs68Gar5jnN0Li1hWET5UyqVoqGhAQqFAi6X6x2Zuv7+vANxoSMn/IWbVADI/RYKlOQPPxiVI4uTy+XQ6XRwuVx04WcYRnAmkE9KyUd04N1mMJyK1J6ix+OBTqeDyWSCyWTCxz72sYLed3l5GZ2dnfTvHR0dSW2OSwmXfVAsZDea6zX8YNHd3Z1GMilWFFzotXyRcH62lvraQlmk/NcRssvMzAxaW1uTzH5TX1fqeUO1Wo1du3Zh586dsNlsmJmZQX19PWpraxGPxxEMBrGyspLXe66uruLgwYPo7OzEol6Po0oltKdPJwUPInfGcRzOnTuXFJTX1tbyOl4hqH3rLVHlT5Zlsb6+Lvge/CCXmnXqdDrodDrEYjHY7fa8vp+p3pPhcJjq1pLAxv9+8AMl+S+QHChL7f3JBz/gymSyNPWjTKIDNTU1SX3KTKIDl2JPUSyyeSleLrjsg2IpwXEcbDYbZmdn0dbWljVYFPrQy2SyJGIF37Yql0h4KYbwfT4fJicnqdloNnISGa0o9bwhMTfWaDQ4ePAgzXxdLleSHJpYrK+v4/e//z10Oh20Wi06/vRP4bvxRjgv9C8lEgm6u7upYHmuPmk5UApZOdKvVZ87J0i68Xg86OrqEnVvlEoldDod+vv7MTc3l7QRySUplytQulwueDwetLe30xJxpoyyEOQKXJlEB8hMoMPhgNVqpcPz/EBZ6aBV7vf3+Xwl8VJsb2/H4uIi/fvS0hLa29uLft9y4LIPisUszvy+B+nh1dfXY3h4OOscYbHl01AoBI7jsLKygrm5OdEWS8VkirFYDOfOnUMkEoHJZMpoLMwH34dRTN8wF0g2HI1GMTg4SLMRhmFgNpsxMzNT8HuHQiGEQiEEAgGa4Wq1WkrimJ+fRzQaLanrRDakKhAV2oeUSqWoqamBTCaDVCqFWq1Gl90umHVyHIf5+fm097jmmmtQU1OTcTG225M9KvMRH+efZywWw9TUFFiWxb59+6BUKnNmlIUEykKyOfLZqdVq6hovJDoQiUSoEDp/lrJUqETQ5b9/MQ4ZfAwPD8NisWB2dhbt7e148skn8dOf/rTo9y0HLvugWChIgCEOD0I9vEyQy+UFu7fLZDL4/X4cP34cWq02ZwAWOud8wDAMlpeXsb6+jp07d4oe5yCUer/fj7fffht1dXXQarWoq6vLu6dJytGrq6vo6+tLYtByHIfTp0+nLcyFIpc7RqlgMBgwMDCAN954Q9C7keO4NFJRvj2+4eFhNDU10UWOLOKx970PeOIJsPF4zqyzoaEh5wYoVRg/36BIWNKLi4vo7+9PM53NlFHyS/GkIiGVSnMGylKVOIVEB2w2G4LBIGpra6mhciwWS1KZKVR0ACh/UEyVkPT5fOjp6Sn6feVyOR599FHceOONYBgGn/rUp6irzaWGzaBYICQSCcbGxhCLxWA0GvOquxeaKQYCAUxPTyMSiWD//v1ZpdKEkE9Q5JeCW1paoNPp6BB1rteRxUoqleKKK65AJBKBz+eD0+mkYuo1NTW0RFVXVyeYhRCXEKvVii1btuDAgQNpi1k8Hi9ZQLwYyJZ59vb20vnSfLF37960+0UX8Ztugue//gvh//5vWLu64Myy6CmVSsRisawbr2L6f4FAAJOTk9BoNBgeHs5ZGiy0R8kPlOXs+xEVnpaWFhrcs6nM8AMl8ajMBiFt0lKDfw6lyhQB4KabbsJNN91UkvcqJy77oJjvbi2RSMBqtcLv90Ov14vuwfCRb8bGH3no7OzE+vp63gExn+O63W7asxseHgYAnDlzJutrsvUNa2pqUFNTQ3fTJGPxer2UsMMwDDUGrqurg0QiwczMDFQqFYaGhjKWoGQyWVE92ouBpqYm9Pb2IhqNZh2rEJp1FIPOzk5a4uMjHA5jZWUFy8vL8AWDwDXXiHq/c+fOIRaLobq6mho319XV0Xvi8XiSfl9MaZ1hGMzOzsLtdmNgYEDUazIhn0AZj8ephyXDMCXpUfLBl2EjyKYyQwLl2toawuFwVtGB1OutBPx+/ybR5nKEGE9FlmWxuLiIxcVFdHd3o7W1FVqttqxzjnxdVDLyEIlECs6McgWPcDiMqakN5/vt27fTnl2uGbx85g2B5LLT1q1bAbzjAuB2uzExMYFwOEz7eWtra6irqxOkxstkMhw8eBBvvPGGqM+gpaUljS26d+9eNDc3w2w2Y3FxsWwOKAS1tbUIBAJlOU5dXR06OjqSymBk41EIQchkMkGtVoPjOITDYfj9fng8HiwsLCAWi0Eqlaa9by6bIafTCYvFgq1bt2L//v1lkRETCpSEKNXY2Eh9I4HkTJeUXwsNPAzDiO4hKpVKNDY2JukBx+NxGihT5dgIkadcma7QGljKTPHdgs2gmAMcx2F1dRVWqxUtLS04ePAgJbsUKvWWKyjyS5dbt25NYrEWy1wVOmeS/a6vr8NoNKYp7WQ6ZiktnYCNRWt5eRk9PT1UKYYYAy8vL1OtU5KpkECp0+lw7bXX4vjx4xn7WTqdDvv27UNVVRXOnTuXRChxPPccuLExNPzRH6Hr+utx7Nixgnu+YuD1eqFQKGCz2Ur+3j6fD0ePHoVEIoFWq0VtbW1ORq5UKkVraytYlsXq6ir9uVwupxUJiURCM/7W1lZwHIe5uTmMj4+nLabT09O0f8zPKKPRKMxmM1iWxZ49ewoy6S4EDMPQSsuOHTuSqiypGSX5/0QikeQnKjZQFhuwFApFVjm2WCyGkydPQiKRJM1S1tbWFl1WFepXXm5eisBmUASQOVMkZUS1Wp02fkBsbQpBtqBI/Abr6uoESTTFzjjyg0Yus1+C1EBX6mDodDoxPT1NtWDJgymRSNKo8WTY2ufzYXFxEX6/n9rq9Pf3w2q1CiretLa20kW4u7ubBsVGs/mdofj/9//w2n33IVpmiyaPxwO/31/Wki/HcfB4PGmlTT6amprQ3t6OLVu2QKFQIBqN4qWXXqLPAtksdXZ2Jn0PE4kEzp49KzgP2t3dDb1eD5/PB4/Hg8XFRbrBiEaj6OjoQEdHR8UC4vr6Oqanp9HR0ZHmngJkLr0C75B3yPed3C+O4yCTyZJ6lQTlIMIQOTatVgu73Y7h4eE00QEiQC/GwDkTMhkMb5ZPN4FAIACz2QwASWVEPooNTqkLIrE3ArL7DYop9WY7LnngyQhJLi9FPko9b0gMlaVSKXbt2iVKvirTsLXf74fP50NjYyNir72GhrNnk0YX+PequrqaDpwLybqVy7eQ9D9jsVjJA6JYBwylUon+/n60t7enBSYiY8cvz09MTGBychLNzc1UP/f06dOCBKG+vj6qe1pdXY3W1lYEAgGMj4/TjD4YDGJ8fBzRaBRVVVVpPcpSlVL54x3ZetJCEJKiI88N2Qzy+5TAO4GS9CnLgVyiA7kMnEmgzPSsC5F4iED45YTNoMgDXxXGeEHlIxOKGasgJrpAsrmw0WhMKpsIoZgHTiqVIhKJ4NSpU5BIJKJHSAhKNW+YSCQoycJgMBRdnpHJZFTYuMdmg+prXwMXjSYNpq+urtKsMhgMoqmpCUtLSyUZis+GqqoqtLa2ora2Fmtra6iurk4aYi4G9fX1YBgGypMnRTtgdHR0oL+/P+N7dnZ2pvWsiYJRJuUeuVyO3bt3J5F7GIaB1WqFx+PB4OBgmtMCsYEim5mlpSUaKPnl8XwDJWk9zM/PC453FAoxgdLr9cLj8aCjowPxeDwpo0x9bSHIlYUKGTinig4Q9nd1dXUSoUepVKaRhArdfL/bsRkUsfFlm5mZwdraGvr6+rKqwhDI5XIEg8GCj8lxHGZnZ7GysiLKXLhYxONxLC0twel0Ys+ePTmDL/88WZZFbW0tjh07htraWjpzqNFo8ioV8Reszs5O6PX6kl+z9MgRIBaDlGUhYRh0TE9D+8EPYuvWrZiamoJSqUR7ezv8fj+qq6vhGRzE61/5CprPn8f6jh0IDA6C5BRiNj01NTXQarVp/cHOzk709vZCqVRienoaXq8Xu3btglKphM1my1pl0Gq16Ovrw9jYWMa+dVdXF6RSKebn5zGQR7abaxPU0tKyIXcnMnBrNBrs27cvqZpCSpbt7e0ZiTR8G6jU0QWfz5cUKFUqVdL4TqZAGQqFMDk5ierqauzfv79kbhKZQIIceZYDgQDdaApllPxZykIcRAoZx8gkOkBIU3wDZ9I7XV9fh0QioWXTYp/R//zP/8RXvvIVTExMYHR0FPv376f/9q//+q/4wQ9+AJlMhkceeQQ33nhjUccqBTaDIkBd2zP11IRQaPmU4zjY7XbaA8gkBVcqENbs0tISHYkQExBT+4bbt28Hx3EIBoPw+Xyw2Wwwm83gOI6OUmi12owuAx6PB2azmfoSlmvBYq+5BlAqwcVigFKJto99DD6Og9VqFSyFx2Ix+IaH4fV6Iff5oI5E6KD10tJSzo0PUcJJRWtrKzweDxVL5wsOdHR0UIseAtLbO3PmDAKBAJaWltDb20tL6qnga6+KzXYJ+SVbxkEqCP39/XR8I9Nn0N7ejp07d9LsIhqNUmuvQog0/NGFTIFyeXkZkQv3iE+2Wl1dxdraWkbd33KBcABSe5apz0BqkCQBkvy+GAcRoXGPQpBKmgJA7c98Ph+8Xi8eeughHD16FOFwGPfeey/27t2LvXv3oqenJ+8guWPHDjz99NP47Gc/m/Tz8fFxPPnkkzh//jxWVlbwgQ98gAqhXExsBkVsDEnnG+AKMRrmz//V1NSgt7e34F1YLmsdjuOwvr4Oi8WC5uZmjIyM0D5LrvfN1DckNjsajYbqFrIsS0tgi4uLCAQC9PeITNry8jId8yhkvjIfsCMjiL7wAiSvvQbbwACsCgX6LngzCn1eSqWSeiUSRCIROByOoioBp06dgtFoxPDwcNJDHo1GBfuJOp0OFouFfg8dDgccDoeoY2WTgFMoFGhvb0dDQwMYhoHdbofFYgEAei9JX4+/GKvVahiNRhgMBhqMVlZWEIlEIJPJsG3bNjqjSxbU5eVl9Pf3J83iFYtcgXJ9fR3nz5+nGdH6+jqi0Sjq6upQVVVVtupLPB6ndnBiNgD87JAglfmay0GknGo25Bj19fVob2/H448/Drvdjk9/+tM4dOgQTp06hR//+Mf4zne+QzfXYjE4OCj482effRa33347VCoVent7odfrMTo6ikOHDpXikgrGZlAsEPlkioRQwrIszVaOHz9e8M6P2Edleq3f76elwr1799IHlmGYrAQP8mCSYCgma5ZKpYLEF4/Hg/n5eTp+oFKpsLS0RDNKMeodhcKh12MaG+zK4Z6evBcSsgBnsr4SA5Zl0dfXRz9DlmUxMTGB+fl5wfesqqoqasETkoBrbW3F3r17Bd+Xz+JdWlqixJnUcRf+/R0cHEQoFIJSqaSZvt/vx+TkJNX8rcQuXyKRQKFQwOPxIBQKYXh4GGq1OqlHSQI4yfrJNZUiUDocDkxPT9PRoULfL191nkgkQh1ESi06AGyUZ/nB3e/3o6WlBbfccgtuueWWkh4L2LCTOnjwIP17R0dHQYL+pcZmUCwQYoIiIdF4PB4YjcakIV3y+kKCIrGPSn0tnyhkMpnShm4rMW9IMlSr1YqtW7diz549kEqliMfjtAS2trZGF1e+Jmqx7EOiQwtANJs1E1QqFXbu3JlTyQcQZn6SAXqycL399ttZLa2i0Sh6e3sFdVALQVVVFfbt25fTVYXRAAAgAElEQVRx4czG4vX7/UlZP19piARKPpHGZDKlEWnKCTL8397eDoPBQL8zQqox/NIrCZTke5dvoIzFYpicnAQAKlheaggFSoZhsLKygqWlJRiNxpIJo6dCyCFDrNLQBz7wAUFRkX/+53/Ghz/84YLP6WJgMyiisEZytqDIsiwWFhboILoQiaaUIx0sy2Jubg42my2rfZSQL2Ip5w39fj/MZjOqq6vTFg2FQpGm3iHUKyI0fRIoxSw8xD5rfX0der0+6RiFQHr8OKRHjqBGoOxDeoWEhNNoNqcxP6uuuw5VVVU4ffo0GIaBVCrN6bl47tw5Krguxp+xpqYGVVVVcLlcgv8+ODiY9wLJZ/ES8MddFhYWEAgEwDAM4vE4GhsbMTAwUPaSOEEsFoPZbAbDMKJ7liqVCs0XyucE/O+dzWZDOBymgZJklfxKBuEBzM3NlZTRKgaRSATj4+Oorq5O0obNJowOFBYohQyGxfZnX3rppXwuC8Claye1GRQLhFDw4KvfZDPfBYpTpiEBlRxvZmYmq39j6jmXet6QaLOGw2EYjUbRWUPqgkVo+j6fD263G/Pz84jFYpThSRYtUrrjC4a3tbVheHi46JKS9PhxqG6+GYjFsFUuR+M//APN/vr7+2l/5NVXX0UgEEibcxzyelE1NARgY0Gbnp5OI9Vkwrlz59DX1wen05n1u6FSqVBdXY2Ojg50dXXh7bffTvp3vV5fssWFHygjkQglV7W1tSESidBAKZVKk4JKbW1tycp7pR6zyBYo/X4/7HY7DZQ1NTXwer3UP7RSM3vEGm5xcVFwVCtX6bUQB5FUdmu5B/c/9KEP4eMf/zi+8IUvYGVlBRaLBQcOHCjb8cRiMyiieMox8I6uopD6jRCKzRR9Ph+mpqZQU1OD/fv35zWcnG/fMBMIs5UwLMXaSmUCn6bPZ8WFQiH4fD44HA46kKxSqRAKhVBTU4Pdu3cXVSrlg4x0SBgGEo5LGm9wOBxQnDiB6uPHoerrQ8BopMxPCcMACgUk111H+6nLy8t59yStViudGcuE+vp6NFksSPzwh7APDABlVuDhE2n0en2aDCCwsaCSjHJubg7BYJB6CpIyZW1tbd7fj3A4jImJibKPWQht0Obn57G0tETnQU+ePAmFQkGvJ5Ngd7Eg10yebbEtlmIdRFJ5Cl6vtyihdoJnnnkGf/u3fwuHw4Gbb74Ze/bswW9+8xts374dhw8fxrZt2yCXy/HYY49ddOYpAEjyHNB8T05zchyHWCyW9+vefPNN7Nmzh5Z08smSpqenodFoRNkx8RGJRHDixAlIpVLs3LlT9PHIA3H27FmEQiHU1tbSnpKQ0HYukFm0lpYWdHd3V+zLnEgkMD09DY/Hg6amJiqgTEZDSEaZaTQkF/iZIiOTJQ3CC5VKnUajaDWZUiHTeRDs27dP0CWjEPCJNH19fXnd50QiQbMvn8+XV6AkLQi73Y6BgYGK6m+GQiFMTExArVZDr9cnXXMsFku6plAoRAMlua5CAyWRXSS9Q7GzxPlCKFDGYjGcPXsWQ0ND1JD60UcfRUdHBz75yU+W5TwuAkTdlM1MsUDE43FEo1G8/fbbgiLauUAYpGJBrHaIY0Rzc7OogJjaN9y5cyedN/R6vVRom4xRkH5epsWKyNHJ5fKKijrzS2hdXV0YGBhI05wU0kTlEyrEZCpkpCP+u99htKYmKdhkkoTL1/y3WGQ6D41GA71eX5KASAQtvF5vwUQauVyeJm5NNjE+nw9Wq5W6QPDvUyKRoG4WQh6a5QLJDldXVzPOOwqN8BALKD6JjHglkmvKFSj52WG5Wbz8jJLfLx0YGIBSqaTP2k9+8hP85V/+ZdnO41LFZlBEfuVT/jC8QqFIGnnIB4RBmguktzA3N4eOjg4cPHgQi4uLOQOq2HlDAkKo8Hq9mJ2dRTAYhFwuT3qoV1ZW4PP58jZVLhZer5eKpGcqoQmxKfklPbIA868plVBBwI6MQDYygq1PPomWX/4Sa9u2wckrlZZLEk4sUs/Dt3cvurq60NzcDK1Wm3OGNef7XyhTd3R0JLE7SwEhFwgSKD0eD6xWK6LRKGpra5FIJLC6ulpU9iUWgUAAExMT0Ol0efemhSygxAZKADQ7rHRGHIvFMDExAblcTp8rjuPw9NNP48EHH8SDDz5YllGMSx2b5dMLiMViWbX+iP7jzMwMWlpa0NPTg7GxMRgMhoLYd4T11tfXl/F33G43pqamoNVq0d/fT5v8S0tLYBgG3d3dGc+V3zcsdDGJx+Pwer1YWlqCy+WCXC5PKlHybYHKgVgsBovFgkgkgoGBAUFh9nzBHw3xer3UtzGVoh986SU0HD4MaTwOKJVY//nPEdi5E5JjxxD/3e+w0NeXlB2q1WqqM1kJNJrN2O5wQPMnf4LE8DC9JlLSK2TsIBKJYGpqClKpFEajsaz3NhX8MYuOjg5aeiV/iAFvrg1NvmBZFrOzs3A6nYIaraUE/7vn9/sRCAQQjUZRXV2Nzs7OrBWaUoOsZXq9nvZRnU4n/u7v/g4ymQyPPvpo0SzuSxCiPtjNoHgB2YKi1+ulpBa9Xk8zw7GxMfplzhcOhwNutxtGgbIbf9hfiPJOZOJShZ1LPWLhdrthsVig0+nQ29sLuVxO2aFerxc+n4+yQ/ljFMVKUbEsS4kdpSDw5AJhHhJBZ9noKIw/+xmaTp+GhOPAyWSI33cfEnffDWDjc37llVeSAqDBYMDAwABYlsWvf/3rspxnKjlrYGAABoNB8HdJ74sfVDJpiHIcRwlTmYg05QIZs0gkEjCZTFmrLvzsi1xTaj8vn0Dp9XoxOTmJ1tZWqiVbCaQSlyQSSVKPkpST+UzeUn3/4/E4dQ8xmUy0XPriiy/iq1/9Ku69914cPny4IoH5ImCzp5gPhCyZyDB4PB4X3EUWwyAVem0ikcDMzAxcLlfasD8f5Z43DIfDsFgsYFkWO3bsSBKRFpLcIuxQsvtkGCZJDzVVQiwbXC4XLBZLmrdiOaFSqdDY2IhwOIyaM2ew9777IInFAI4DK5GAlclwrKoK0lOn0NDQgLq6OnR3d2NiYoK+x+zsLDweT0aT41yQSCSora3F7t27cfbsWWqozEfq92VpaSljUMwkX5cqtk02OnV1dSVl8eYCv5fV19cninCWqUxJrok/SpEtSyb9Up/Pl2Y6XG6Ew2FqpcXvHfKvKVPftVgmL5F97O3tpVJtPp8P99xzD9xuN37729/mLeH2XsRmUBQAPzjxywupKET/lIBPtCE7x4WFBXR1dQmaofJBgmKp5w0LGYIni3ltbS0leLAsS0kvYok8JBBzHFe0Gk2+IGLlDQ0N2Ol0QhKPQ8KyYCUSrO7cifE//VM4e3shtduRSCQQDAap2gvfkFesVmkqVCoVTCYTOjo6YLVaBQOiEEKhUF79Q/6GhnzHPR4Purq6EIvFMD4+npT5kz+lHoMgpJKqqqqixywyEV9Sh/OJ3JtUKsXq6mpZ+qXZwM8OTSZT1p58tr5rtkBZU1OTcf7QbDYjFoth7969UKlU4DgOR44cwT333IPPf/7zuOOOOyqWKV/q2AyKPJCy3eLiIrq6umAwGLJ+UcSSZbK9luzeGhsbMTIyIqr0SLLMUs0b8nftHR0dRQ/B81mfBJmIPGq1GrFYjHpYlrJ0R5Rp2GuuATsykvbzeF0dXBYL3Lt3Y8dtt6Gmpga+vXvRIJNBwnFg5fKNgJgydrG4bRv6+/tRXV2d5FaRL5RKJbZs2ULntIgIgljkO85DQLQ7Ozs70zZgxFbI5/PB6XRidnYW8XicWoalCijkA0JSs9lsZSWVCAVKYhweDodRVVWF5eVlrK+vi7KkKhZkxEOj0RTMLM0VKMkzxQ+UGo2G9uW7urrQ1tYGiUSCUCiE+++/HxaLBc8991xGbsLlis2e4gXY7XZMTEygubmZ9s9yYXl5GfF4HD09PXkfz+1249SpU2hsbITRaBRt9ksWrbfeeov28Orr6wuaNQTeYXZqNBr09fVVVLHDZrPBarVSZiERcC4FkYc/bwilEtEXXgA7MvLOzyMRgOMAqRRQqRD+1a9wXCqFw+FImzvMNRcICOufZoJWq4VUKsWuXbuSSvJjY2NJ6jcymQw6nQ4ymQxyuRxyuZz+v1qtRmtra173vFAiDb9ETv4QR3d+UMl2Ln6/HxMTE2hsbERvb29Fs5JMAt58AXGfz5dmSVVsoOT3anNlh6UCISh5vV6srKxQJq/ZbIbb7UZjYyMee+wxfPrTn8add955uWWHmz3FfCCTyfIer5DL5QiHw3kdh4iEe71eKJVK7NmzR9Tr+H1DhUKBQ4cOIRQKCc4akqCSre9AxMOj0SgGBwdLwuwUi2AwiKmpKahUKgwPD9PFmW8J5PF4sLCwkCbzJpbIw1em4WKxjYxxZASx3/4Wqmh0g0ADQMKy4GIxTD/xBBwf+QiAdMeJTHOBBGKCJj0vqRQcx0GlUsFmsyEYDFKPw/n5+aTfNRgM0Ov1+X68aeAvzgaDIW9WYaYSOfHWJHZULMsmiYeTgE+Ew7dt21bR7xmxSuM4TlCiLVVAPJPJMdHjzWejFgqFMD4+Dq1WWzH3EGBjTSIl4q6uLnR0dIBhGCQSCTz22GMYGxtDVVUVnnzySUxOTuLOO+/Etm3bKnJu7xZsBsULaGhoyLsUmg/Rhi8S3tvbC5PJhKNHj+Z8Xba+odCsIXmgrVYrgsFgmhOFUqmkSiGp5rflRiKRoAuk0KxjJu+8TEQevh5qmldditlw5OBBTI6NoaqnBzvJz1kWnFQKRibLOnOYaz4xV9BMOi+WRXNzM12gPR4PbDYb1tbWkoheVVVV6O3tFf/hZoDP58Pk5CQaGhpKujhLpVJBb03SS15ZWYHb7UYkEoFGo0FbWxvd1JU7OylUwDuXd6PX68Xi4mLWQEk2IDabTdCpppxgWZZuuHft2kWrT+Pj4/jyl7+Mj370o/jFL34BuVyOQCCAt99+u6Ln927BZvn0AshuKh+Qh2THjh0Zf4eIVhNJtN7eXrowvfnmm7jiiiuyvrbYecNYLAav1wuv14v19XUEg0FUV1djy5YtNKiUS0+SIFWNZuvWrUUFYv7i6/V6k4yN+UQe2ejoO2bDra10EyAbHYX0yBEwOh1mRkexajLlLHlmK4/yM0VOJsPqTTchevgw3q6uTiNiiS2z6nQ6SKVSVFdXJy2+YsvbhEjj9/thMpkqnqGRMQuj0Zg2n0cCKp9JWapAGYlEMDk5CYVCAaPRWJbvNl+4nlxTNBqFUqlEKBRCXV0djEZjRcliPp8PExMT2LJlCzV/jsfj+MY3voEXX3wR3/ve97Br166Knc8lis05xXzAsmzeTNJgMAiLxZKxBEp0I6uqqmAwGNJKs5mCYqlHLAjJQKVSob+/HyzLJs0aFjNCkQtEuFyj0aC/v79sAZifJXu9XjpDGI1GodVqodfr08rJpZgpVKvVqK2tRZPFgs5XXoHmF78AGAacQoFjDzyApc5O+rtiy6z19fW48sorAWws8uQ++Xw+JBKJNHZoajmZZNSl2IDkg9Qxi0zzpXw7Kp/Ph0AgUPTIAV83tJAScTEg8nA2mw1btmyhfb1YLFbwpkYsiPiAy+XCtm3b6HjJ5OQk7rrrLrz//e/HfffdVzGuwCWOzZ5iuZFpJCMajcJisSAUCmFgYCBriYJPqS91MIzH4zRbMBqNSedRU1NDZ5JIfyi1P0keZK1Wm7fMFt9OqhI9S0JK0el09POPRqPo7Oyk/VNCzS+VyMDQ0BBaW1tx8uRJnK+rA8Oy2HGhjMpxHLqsVqz191OxebFl1u3bt9PPmriGkHtFdGuFysk1NTVwOp3UYLiSC2E4HMbk5CRUKlXOMQsh30a+gk0mTdRM30G+gHcl+3fAxsZ4fHwcOp0OIyMjSZtJPpM31QqN33ct9D4FAgGMj4+jubmZGkozDIPvfOc7ePLJJ/H4449fElZM7zZsBsUikNpTJGQJYvbb2tqaNZBIpVLaYynlvCHLsnTX3NPTkyaeLXQemfqTXq8XMzMzgv1JIVJSpdVoMh27v78/rV/K7w+tr6/DarUWfKxdu3ahra0NJ06coPOJqb3Hqba2JPcV/r9zCgWihw5BpVJRw2Jgg1yTbVRBIpFArVZDrVZj69atAN4ZRl9aWoJarUYoFMLp06fTysnl6OVxHIeFhQXYbLainB0yiYeTQEn0Q/mD+RqNBmtra7Db7RVjdxLwxcMHBwcFLZYkEglqamqSNqBCIy+JRCIto8y2qSDHXltbSxIVmZ+fx5133ok9e/bg9ddfr2j59r2EzfLpBRRjH3Xo0CHY7XZYrVZs3boV3d3dohagEydOYMeOHZDL5SUJhsCGfuH09DSlvpdy10z6kyRYRqPRJIk3hmFgtVrLcuxc4A/gizm2zWbDyZMnCzqWQqHADTfcgDNnzmBpaSnp3/g9QwBp/UOiWVr/P/4HnZ1kGIZKluU7gsIn0vCvO1eJspDsPxWkPUBkACtxv8lg/vr6Omw2G6RSqaAebzk3YvzssK+vr+jNRurIi9/vTyuTk9lQwmolVl5kY/2jH/0I3//+9/Gtb30L1157bYmu9D2HzZ5iPig0KP7hD3+AUqmEWq2GwWAQXQrhOA5nzpxBVVUVlQ4rpt8WCoVgsVggkUhgMBgqsksku9719XUsLi4iFovROa9szNBSgpRK4/E4BgYGRM97RqNRvPzyy3mbAAMbw+FEfYYPvsJNtv7hzp07ix6YLoRIkyqyzXcN4Wf/uQIK2fx4PJ6Kj/MICXjzSS8+n6/gMYpcEJMdlgpCs6HhcBgMw6CtrQ3BYBA9PT2IxWK466670N3djQcffLCsgubvAWwGxXzBL2PlAtFFXVtbw8GDB/M2+yXGni6Xi37p+Ua5YkteiUQCs7OzcLvd0Ov1ZTMmFQLDMFhYWMDq6ioVkk5lhvJ9DcnCWwoboFylUjGwWq0YHx8v6jwIqqqqcOWVV0IqlSIQCED18MNofuSRDbk4qRRjhw9j8iMfQV1dHQ4ePFhUv6+URJpswuHkfvEDisvlgtlsxtatW9HZ2VlR4WixAt6p7NBU4fpCennEWopk5JUceufrpba1tSEQCOC73/0unn/+edjtduzevRu33HILhoeHMTIyUjGP03chNoNivhATFEkQcjgcMBgMmJ6eTmuwC0GMTinLslQKzev1JnkakkDJdzYgYw6dnZ1ob2+vKMtwfX0dMzMzlAKe7fr5GQphhiqVSrrokusSi3xLpdmu46WXXsprMyQEmUyGpqYmJBIJSuRpsliw5c//nCrqhJ57DtELupOF3qdwOIypqSnI5XIYjcayEWkyZV7RaBRSqbTsFkup4At4Dw4OFiTgze/lkT9i5OtYlsX8/DwcDgdMJlNZs0Ohc15ZWcHi4mKSLJ7D4cAXvvAF1NTU4OGHH8ba2hpOnDiBkydP4u///u/R0dFRsXN8l2EzKOaLbPZRhPJNglBHRwekUineeust7N69O+sCVcy8IfE0JL28SCSS5GxQ6XmoYDAIs9kMhUIBg8FQcFmKPxBNFt5cFlSFlkpzXc8rr7xS8OtlMhkOHTqE+vp6SuQh18S9+SbUJ04gPDIC2VVX5aXIwwdfM7QYMkshIGMWVqsVTU1NkEgkaaMhZIynWDavEEhmSnwWS7nxyyVfp1AosLS0dFGk6SKRCBVNNxgMlHfwwgsv4Gtf+xr+8R//Ebfddtt71eKpXNgMivkiU1AkD2Z9fX3arN3p06czLtClHrGIRCJ01KCpqYkGFoZhkggUarW65A8wv0wrpEZTLDItUIQVG4lEqHtHqVV4Zmdncf78+bxfJ5FIcODAgYwuKsA718WfNSRyaGL6rsTLs7GxET09PRUlL/HHLAwGQ9L3nj8akum6CtXjBd5xdohGozCZTBXb+JHyP7ECUyqVlPHLzyjLdR/4s558+ziPx4MvfelLCAaDePzxxwsWg8+GSCSCa665BtFoFIlEAh/96Efx1a9+FbOzs7j99tvhdDqxb98+/PjHP363zj1uBsV8EY/Hk4gXoVAIU1NTAACj0ShYtjl37hy6u7uTyiqlDoZk1GNtbY32z/ggDzLJKAnTkF92FUOgEAL/Ia10mZZlWSoaLpfLqRtIqfuTHMfh6NGjcLlcol+j0Wiwffv2glw9siny8Pt4ZBTGZDJV1POPr5WaT2bKvy7CouQHFLF9cqIA1d3dTZ0dKgUy+9fU1ISenh7K7ky9LgBpG4BiN6KxWAwTExO0PK5QKKih9Ze//GV88YtfxJ//+Z+XLWMlGx21Wo14PI6rrroK3/rWt/Dwww/j1ltvxe23346/+qu/wu7du/HXf/3XZTmHMmMzKOYLEhTJ0DvJirKpY0xMTKC1tRUNDQ0l9zfkOA5ra2uYnZ1FW1sbOjs7RT8QZMaLZCjEMocESTFsV5/PB7PZDLVaXVY1GiEQy5tYLJaUiQv1J7MRQ8QiGAziyJEjSebNBFVVVaivr4dWq6X/LfVOmT8Xura2Br/fT8Wq82GGFotSj1lkGw1J3djwBbwHBgZKwhgVC5ZlqZeomJ4pwzBJgTJVajBf+TpCnuL7twYCAdx3332Ym5vDE088gU6eOlK5EQqFcNVVV+Hxxx/HzTffDLvdDrlcjqNHj+IrX/kKfvOb31TsXEqIzaCYL2KxGObn57G4uIju7m5RWZHFYoFWq0Vzc3PROqV8+P1+mM1mVFdXQ6/XF70I8xl5JKPMJO9GnDxCoRCMRmNFSRWEVbqysiJasJzfx/N6vQU7a3g8HoyPj8Pv96O+vh5dXV3Q6XQVY/ORcqVSqYTBYACApA1AJBJJGjUoZXBmGIaWx00mU1nveepoSDAYpCNRhNVaiQ0AAbG1am5uFj1jLIRC5Ovi8TimpqbAsixMJhO9n2+++SbuvvtufOYzn8FnP/vZivUzGYbBvn37MD09jTvvvBN33303Dh48SH0+FxcX8cd//McYGxuryPmUGJtBMV9MTEyAYRjRforARj9KJpNRS51iAyIJSMFgEAMDA2VdnPjybmR8IpFIIJFIUNp7KcqTYlFKVinpT5JgybJsUnaSWu6KRqMwm81gGAYDAwMVJS8RBxW73Z61XJlK5CGjBqm+hvkSXi7mmAUR8JZKpWhpaaH911J6a2YCPzssl62V0AaAMMolEgkcDgf6+vro+hGJRPDAAw/g1KlTeOKJJ0piHVYIPB4PPvKRj+BrX/saPvnJT15WQXFT5o0Hg8EgWD4TAimTarVaWCwWLCwsoLa2lpYn823G8zMkYi1V7sWJL+9WW1sLi8WClpYWNDQ0wO/3Y3p6mpYn+eMTpS4d8kul27dvL7p/lsn/j+ziFxYWknbxDMPA4/HAYDCIthkqFfhEmgMHDmTNCPjWRoRowSfyrK2tYXp6WjSRJx6Pw2w2IxaLYffu3RXdCOQS8E61bOJ7a4qVQ8sGfna4f//+smViQvJ14XAYExMTiEQiqKmpwZkzZ3D48GH09/fj1KlTOHz4MF566aWysHnFor6+Htdddx2OHj0Kj8eDRCIBuVyOpaUlahX2XsVmpsiDGPsofs+QXyolTWp+1kVEtXORQtbX16m1VHd3d0UZhoTRyjAMjEajIIuW79Tg9XqRSCToBqAYNl4hpdJSwul0YmpqCgqFAgqFIm1wvRwbAALSty4HkYZPvErtd5ENWyAQwNzcHHp7e3Nq9JYaRMC7trYWer1e9OLPnzXkVzZSM+Vs30W+q0Sl1XiAd7Lyrq4uSiKKx+N44IEH8NZbb2Hnzp1YWFjA3NwcRkZG8L3vfa9i5+ZwOKBQKFBfX49wOIwbbrgBX/rSl/CjH/0It912GyXa7Nq1C3/zN39TsfMqITbLp/kiV1DMd94wEymELLhyuRyzs7OQy+WC1lLlBBlK5qvRiEWmDQBZdMXoapaqVFoI4vE4zYIHBgaSFkZ+3zW1PFmKeTxCnrJarRVlVyYSCfj9fjidTiwvLydlk5Ui8vDFw0sl4J1pNITfxyOZMskOW1pacgpOlBoMw2B6ehrBYBDbtm2jz/r4+DjuuusufPCDH8S9996blPn6fL6KigWcPXsWd9xxBxiGAcuyOHz4MO6//35YrVbcfvvtcLlcGBoawk9+8pOKkqBKiM2gmC8yeSqWcsQiEonA5XJhcXERwWAwKUiWa8YwFQ6HQ7QajVjw2ZOkd0JUa/hZF79UmmnMpVzgj5f09PRgy5YtOe+j0JxhIXJ8QDqRppKzXkJjFnyJt3ITeYhMWiXEw4VGXqLRKDiOQ2dnJ5qbm8vmGiIEj8eDyclJdHR0UPIewzB49NFH8V//9V/47ne/i3379lXkXC5zbAbFfJEaFEs9b0j6KITd2tbWRne6Ho+HzkCVasYwFWTuUqFQQK/Xlz0z5ZNCPB4PQqEQGIZBc3Mz2tvbc5a6SolAIICpqSnU1NRAr9cXNV7C708KzYXW1dWhurqa3jOxRJpygYxZEGeFTJ95OYg8+Y46lBrEkb65uRk6nY7KKAaDQdF+jYWCZVnMzMzA6/Vi27ZttDVhtVpx5513YmRkBP/0T/9UtudwcXERn/jEJ7C6ugqJRILPfOYz+NznPgeXy4U/+7M/o5vDp556Kqtd2XsIm0ExX5CgmKlvWAzcbjcsFgvdKWdbWPgzhmQHX11dnVTqymdhKrcaTS6QUqlOp0Nrayvtd/n9fnAcl1R2zddxPReIo4Pb7c5p+FwM+N5/Xq+X9idVKhXcbjdaW1vR399f8ZJdsWMWmRR5hMqTqSAC3qRXXslrZ1mW3ne+Iz0f8XicBkkyy6tQKJICZaEbUhKMSTVGIpGAZVn84Ac/wA9/+EM88sgjuPrqq0txqRlhs9lgs9mwd+9e+P1+7Nu3D7/85S/xwx/+EA0NDbjnnnvw9a9/HW63G//2b/9W1nO5RLAZFPMF2SmXct4wHA7DYrGAZdmMRBYx5xUOh5PILhSUsTcAACAASURBVPwRA1J2TT1Xfrmwo6Oj5NqRuSCmVEpmu/gi6AqFIq3XVQhImbgcupm5QObP/H4/tFotwuEwzbr4YwblypTLOWaRSuQhTijk+6hWq2Gz2eD3+wsW8C4GJBjzA5JY8D1DidZwPiVlPpGHH4yXl5dx5513Qq/X48EHH6z4ZwIAH/7wh3HXXXfhrrvuwquvvoq2tjbYbDa8733vo8pd73FsBsV88fOf/xy//OUvsX//fhw4cAA7d+4suKfCMAwtG+n1+qyqOIUgk6MGCSYymQyzs7Oora1Ff39/RftXxbJKycKUKhYudhifuEnIZDIYjcaKkgI4jsPq6ipmZ2fTiDR8UkhqppyPDFo2xONxqo9bSc1QQuSx2+2w2+2QyWRpAu/lJvKQ7JD4PJYi8KSOhgjZUJHRECIRxxcBYFkWTz75JB555BE89NBDuP766y+KiPfc3ByuueYajI2NoaurCx6Ph16fTqejf3+PYzMo5ot4PI4zZ87g2LFjOH78OMbGxlBbW4v9+/djeHgYBw4cyOlfx18USYZSqbIR8Wecm5tDOBym5sf8rKvcPbxysEqFhvGFggmAi9q7K4RIQ+TC+OMThfSU+axWsSSiUiKRSMBisSAcDmNwcBDV1dWUyMN3eCkXkaeY7DBfpFZuSOmVZVls3boV4XAYvb29iEQi+NznPof6+np885vfvGh9u0AggGuvvRb33nsvbr31VtTX1ycFQZ1OB7fbfVHOrcLYDIrFguM4uFwuHD9+HEePHsXx48dhs9mg1+sxPDyM4eFhDA0NUVLF+fPnEQ6HoVar0dfXV3F2ISHxkEURAO0HkQcYQNKCWypyQaVZpamZMskoNRoNOjo6oNPpivIuzPdciGC70WgsevHLpFvLz7r43y2iCkPsvCrtYCBWwDvV/Nfr9aZ5GubbLyc9Y6/Xe1FKtcFgEBMTE9BqtWhtbYXf78e3v/1t/Pd//zfW1tYwNDSE2267DcPDwzkt5sqBeDyOW265BTfeeCO+8IUvAAAGBgY2y6fZfmkzKOYHhmEwNTWFo0ePYnR0FKdOnUI8HodMJgPDMHj88cexc+fOipIK+ESWXCQe/uiEkKB2vrt3juOwtLSE5eXlizKAH4vFYDabEY/H0d/fn9QT4hOUyPWVWiXE4/FgamoKzc3N1FWh1EhlhfKDCenvDQwMZLWwKgfIZ88wDEwmU0Fl6lIQeYhYfiW/d/wRl8HBQUrgcrvd+OIXv4hEIoFvfOMbsNvteOutt3DixAncfffdMJlMFT3HO+64Aw0NDfjmN79Jf3733XejsbGREm1cLhf+/d//vWLndRGxGRTLjXg8jocffhg//elP8dGPfhRSqRSjo6OwWq1ob2+n2eT+/fuh1WpL/tDyTXeLyc7Igkv+8BVr+ELhqbiYA/j8zJjYaQkRjVLLXKXyniQCAOFwGCaTqSSGx/nA7/fj/PnzUCqVqKqqQiAQSGLyEn3XcgQKfougr6+v5N5+mYg8JEiq1WrY7Xb4fL6kUYdKIRwOY3x8HBqNBv39/ZDJZOA4Di+99BLuu+8+fOlLX8LHP/7xi24A/Prrr+Pqq69O2qT/y7/8C0ZGRnD48GEsLCygu7sbTz31VMVbDRcJm0Gx3GAYBv/xH/+BT3ziE0m7ZDKXdvToURw7dgxvvfUWwuEwdu7cSQPltm3bCs5a+HNvQv6KxYLjuDR/Rr5kXU1NDRYXF4sOxoXC5/NhamoKWq0WfX19ec/NFeM9yQ8IF6N3R8YsiEwZf8yC79LAJ1+VkuzCL9USz79KgBB5VldXYbPZKJGHT76qhCLPysoKFhcXMTAwQMvkfr8f9957L1ZWVvD973+/bNqgn/rUp/D888+jpaWFCnJfxjOHhWAzKF5KiEajOH36NI4dO4Zjx45hYmIC9fX1NEgODw+L0qAkOqmtra0Vnf1KJBLwer1YXFyEy+WCQqFII/GUe4FMJBKYnp5GIBCAyWQqmW6lWO/JUCiEyclJVFVVpTnRVwJutxtTU1P/v70zj2rqTMP4k7DKIiLIIosoURZFbSCgrXasOONSBzvao7hOq4y0lUqlbj2OFcYB0RlaRRytC3WZ0apUh+NSxLqOIybiMkoJiCzKJgiyCAgJyZ0/8N7eIGtIbli+3zmeI9GTfIFw3/tuzwN7e/sOD5O0tPNKvzc6WHbkfbADQksC3tpGoVAgOzubWfMwMTFpc5Cnpd5rV6ivr4dUKmV+9vr6+qAoCjdu3MDatWuxYsUKBAUFafX38fr16zAzM8OSJUuYoLh27dq+unOoDiQodmfoacFbt24x/cmysjKMGDGCmXQdM2YMs6OXk5OD6upq6Onpca6TCrxZKuXz+YxQOH1R0lRpsjltrTlog5a8J+npQnt7e9jZ2bVaUtYGmlyzaP7eqqur2xV4V1fAW1PQMmnt7ZtqY5CHvevLNhx/9eoVIiIikJaWhgMHDmDo0KEaea/tkZeXh5kzZzJBsQ8PzagDCYo9jcbGRqSnpzNl1wcPHoDP58PY2BjPnj3Dvn37tGpz0xKdmSptrTTJ1j/tbDCvra1FZmamzrIz9iCNlZUVM/HK7nPR740t7aYJuFqzaEngHQDMzc2ZsqWnpyfnZTlaRLumpgaenp5q3Qx0RjC8OTKZDFKpFPr6+iql4tTUVKxatQqLFi3CypUrOe2lNw+K7PWKPrZzqA4kKPZkKIrC8ePHERkZiUmTJmHgwIFITU3FkydP4OLiwmSTQqFQKwMV9FRpQUFBq4MsHUEul6sM8TRfxKeFBprD7p3pQpqOzs7q6+tbHaShXVDojEuT3pN0746+IHM9yl9VVYX09HQYGhoypWNN9yfbgi4Va0ONiK1d29IgT//+/VFbW4ucnBwIBAKmZy+TyRAdHY0bN25g37598PDw0NiZOkpbQRHoUzuH6kBMhns6paWluHLlioqtEy0ynJKSgjNnziAiIgJyuRxjxoxhepNubm5duntll0p9fX279FwGBgawtrZm3gN7EZ82xW2uf/rq1StkZ2dj8ODBnGfGnXHSaMlAli4pV1ZW4smTJ532nmSvuOiid8cW8B45cqSKdRHbVaOoqEilP6mpvjI7O9SW8TGfz2fOTEPf4FRUVCA7O5v5uV28eBEymQyDBw9GVFQUAgICcOXKFc4rFq1ha2uL4uJipnzKtUl2b4Rkir2Auro63Llzh1HiefToEWxsbBi5Oh8fH1hZWbV7t02XShsaGuDm5sbZVCk9NUl7/SkUCpiZmcHS0pLJTLiQatPGIE1nvCdpe6X23Cy0BS1i3VEBb3X6k21BZ4dsiyUuKSsrQ1ZWFrNmIpPJcP36dezfvx/37t2DsbEx3Nzc4Ovri7lz53K6c0jTPFPswzuH6kDKp30VelKQ7k2KxWJmp4suu44aNYopyTU2NiInJwfl5eVdKpWqC3vFhM6O2HqTVVVVKmLa9O6kpoIGW5HGzc1N66XalrwnaWNXepCIy3KpJlVhlErlG/qu7JuA/v37v+GEolAokJWVhbq6OkYijksaGxvx6NEjyGQyeHh4MDdgWVlZCAkJwYQJExAeHg5DQ0M8ffoUt2/fhpubG7y8vDg95/z583H16lWUlZXB1tYWERER+OCDD/rqzqE6kKBI+BW5XI4HDx4wcnVpaWkwMTGBs7Mz7ty5gz/96U8IDg7WyWRhZmYmrK2t4eLi0qbXX0sZF7vHpY5kXUVFBZNZc21vRL8+PchjZmbGDPKwfQy1qVtLv7423DRomjuh0P1JCwsL8Pl8lJSUwMnJiXMnE+BXNxFnZ2dmqlmhUGD//v04cuQI4uLi8Pbbb3N6JoLWIEGR0DolJSUICwtDZmYm3nnnHUil0jZ1XTUNu1SrriIM3QeisxL2oAv9p7USaEcGabQJ+/Vbyo7actTQhPdkSwLeXEIbXtfW1sLY2BhyuVxFkk/be6/s7NTT05OZis7Pz8dnn30GT09PbN26ldPPRVJSEkJDQ6FQKBAUFIT169dz9tp9BBIUCa1z7Ngx9OvXD7NmzVJxiKd1XcViMe7evQsAEAqFTH+yq0a5dGn36dOnGDZsGGxsbDQadNm7k1VVVUx/ki19Vlpairy8PAwdOrRDggmapqSkRK01C015T9K9My52PluCzs6cnJwY1xm6P9l879XMzIzJljW1G0rvPbJ7l0qlEv/85z+xe/duxMTEwN/fn3OlohEjRuDixYtwdHSESCTCsWPH4OnpydkZ+gAkKBK6Bi33lpqayggM0FOhvr6+ndZ1ffnyJTIyMtC/f3+4urpyUqqle1xVVVUoLy9HeXk59PT0YGNjA0tLyw5bM2mC+vp6FZ9HTfQNO+M9qQkB765AZ6d0dtxeAGf/7OjVifb6k+09X3Z2NqqqqlQ0U589e4aVK1fCxsYG3377rcpUKlekpKQgPDwcFy5cAABs2bIFAPDVV19xfpZeDAmKmiA8PBz79u1jdpWioqIwY8YMAE0f3AMHDkBPTw+xsbGYOnWqLo/KCW3putLZZHNdV5lMhtzcXFRXV8Pd3V1Fr5OrM+fl5eH58+fMVG1z6TO6dEdfbDUZsLlcs2jNe1JfXx+1tbVwdnaGs7Mz573Tlnp36kAPKbFL5gYGBm/sTzaHnqxl+y1SFIVTp05h27ZtiIqKwsyZM3Um4p2QkICkpCTs378fAHDkyBGIxWLExcXp5Dy9FLKnqClWrVqF1atXqzyWnp6OH374Ab/88guKioowZcoUPHr0iPMxeq7h8/lwcXGBi4sL5s+fD6BJ1/X+/ftISUnB9u3bmbUCHx8fAE29khMnTsDHx4fziw49SGJnZweRSMQEAysrKyY4sd00aG9AWvWEDpTqCiTU1NQw2bFIJNL654PH48HU1BSmpqawt7dHQ0MDpFIplEolnJyc8PLlS0gkEpVFfHrlRRs/G3Z2OHbs2C7LE+rp6cHS0lJFXYdtF1ZQUICGhgbmJsfc3BwVFRWorKzEqFGjmMna8vJyfPnll9DT03tjF5jQtyFBUU0SExMRGBgIIyMjDB06FAKBABKJBOPHj9f10TjHyMgIfn5+8PPzA9AUZMRiMVasWAGKomBvb4/58+e3quuqDdjydO0tgfN4PJiYmMDExAT29vYAVE2M8/LyGMcJ9hBPW+VHpVKJ3NxclJeXw93dXWUJngvYvdvhw4e/cdFnKw3Ri/ia9p4sLy9HVlYWnJ2d4e7urrUbIkNDQwwaNIip5tA3OaWlpfjll1/A5/NhYGCAf/zjH2hsbISZmRkOHz6MjRs3Yt68eTq3eAIABwcH5OfnM18XFBRozW2D0DYkKHaAuLg4HD58GD4+PoiJiYGlpSUKCwsxbtw45v84OjqisLBQh6fsPlAUhZiYGMTGxuKdd94BoKrrevjwYTx48AAGBgbw9vZmyq6aWImgKArFxcV48uRJlwZpWlI9Yffv6IyE7t+x1yboNRNbW1vOFXmAJrFqqVQKExMTiESiFoNbS0pDdLZcVlaGnJwctQXe6b2/hoYGjWSH6lBaWorS0lIIhUKYm5tDqVSisbERsbGxePz4MYyMjLBjxw7cuHEDwcHBnO8cNkckEiErKwu5ublwcHDADz/8gKNHj+r0TH0V0lMEMGXKFDx79uyNxyMjIzFu3DhYW1uDx+Nh48aNKC4uRnx8PEJCQjBu3DgsWrQIALBs2TJMnz4dH374IdfH75FQFIWqqircvn2bmXbtqq5rbW0tMjIyYGJiAoFAoHUpLrZrPLs/yefz4ezsjEGDBqm1O9mV89Bu8Gy/P3VRx3uSzg6HDBnCudck0PQZSE9Ph6WlJYYNGwY+nw+KovCf//wH69atQ2hoKD766CPw+XzU1NTg7t27cHZ2houLC6fnbInz58/jiy++gEKhwNKlS7FhwwZdH6m3QQZtNA1bYqn5dNjUqVMRHh6udvk0JiYGq1evxvPnz2FtbQ2KohAaGorz58/DxMQEBw8ehFAo1Nh76Y7Q04G05+SdO3cgk8kwduzYNnVd6VJlWVkZ3N3ddTI9WFpaiuzsbDg5OcHU1FRFrcbIyEglkGhDrYYribjWvCdp4QEAGDlyJOeTrewbAg8PD+YzUFdXh02bNiEzMxMHDhzAkCFDOD0XoVtBgqImoMV2AeDbb7+FWCxmBmwWLFgAiUSCoqIi+Pv7IysrS62LUX5+PoKCgpCRkYE7d+7A2toa58+fx86dO3H+/HmIxWKEhoZCLBZr+u11e1rTdfX29oavry9qamqQmJiIv/71rzqZquzImkVDQ4NKNimXy1UMmruyf8cW8NZV77K4uBjZ2dkwNTWFQqF4YzdU296Tr169Qnp6OszNzeHq6sr8DkokEoSFheHjjz/GihUrOPtsnDx5EuHh4ZBKpZBIJMzAGdA3J9a7EWT6VBOsXbsW9+/fB4/Hg4uLC7777jsATXfDc+fOZdYPdu3apfbd+apVq7Bt2zbMmjWLeSwxMRFLliwBj8fDuHHjUFlZqRKg+womJiaYOHEiJk6cCODXAZLk5GRs2rQJ1dXVsLe3x+bNm1vUddUWFEWhsLAQBQUFEAgEbU4vGhkZwcbGhnEwoPc/6d5kTU2NimRdR70Zq6urkZGRgUGDBumkdymXy/Ho0SM0NjbC19eXyQ7Z+4UFBQVa856kPwv5+fkq5eKGhgZERUVBLBbj2LFjcHNz6/J77QyjRo3CqVOnEBwcrPJ4X51Y72mQoNgOR44cafXfNmzY0OW6f2JiIhwcHDBmzBiVxwsLC+Hk5MR8TQ/y9LWg2BwejwczMzPExsYiPDwcc+bMQWNjIx48eIBbt27hu+++Y3Rd6QEeX19fRjlFE9TW1kIqlaq9ZkEvoJubm8PR0RGAqjdjSUkJU5ZsyZaJFvCurKzEyJEjOXMzYUOvrrQ0zMTn89t8f6WlpW9I8vXv379TNzL19fWQSqXo16+fys/gwYMH+PzzzzF79mxcvnyZcy1fAK36LJKJ9Z4BCYoc0NYgT1RUFJKTk3Vwqp6LhYUFxGIxcxGlp1i9vb2ZNZAXL15ALBYz065FRUUQCATw8fGBSCSCUCjsdLaizTWL5t6MbNkz9jSooaEhampqYGdnB6FQyHmWwc4OhUJhh3uHbXlPVlRUIC8vr0O2U2y/yxEjRjC7pnK5HNu3b8dPP/2EAwcOYPTo0Zp70xqCTKz3DEhQ5ICff/65xccfPnyI3NxcJkssKCiAUCiERCIhe0vt0FZWwePxYGVlhRkzZjDqQ2xd1x9//JHJ8GldV5FIBIFA0GoJkus1Cx6Ph379+qFfv36ws7NjluBfvnwJR0dH1NXV4fbt29DT01PZndTm+kNb2aE6GBsbw9jYGLa2tgBUnVCKiopUnFBoF5S8vDwYGBjAx8eHyZwzMzMREhKC9957D9evX+fEdqutG112G4TQ8yBBUYd4eXmhtLSU+drFxQWpqamwtrZGQEAA4uLiEBgYCLFYDAsLi06XTjdu3IjExETw+XzY2Njg4MGDGDx4cJ+cbOXz+fDw8ICHhweWLl3KXIBpXdfw8HA8fvwYDg4OKrquCoUCe/fuxdSpU+Hl5cW5mwagKuDdfAmevYRfWFjI7E6yZc+6WkKUy+XIzMyEQqGAt7e31oIOXRo3MzNjbgBpWbfCwkJkZmZCX18fJiYmCA8Px8iRI1FQUIAzZ85g9+7d8PX11cq5WqK1G9226OqN7tdff42BAwfiiy++ANDUvrGxsUFoaGinz0JoHTJ92o1gB0WKohASEoKkpCSYmJjg+++/V5li6wjV1dVMiS82Nhbp6enYs2cPmWxthea6rsnJySgrK8P48ePh7+8PPz+/N3RdtQk7GHVUwJu9hM/WPlXXcorODmk3eq6hvwdKpRLu7u4wNDREfX094uPjcfr0aRQWFqJ///7w8PCAr68vgoODYWZmxvk5W2LSpEn4+9//zvzednViPS8vD7Nnz8bdu3ehVCoxfPhwSCQSrWrp9jLI9GlPIy8vj/k7j8fDrl27uvR87J5XbW0tcyEkk60tQ+u6Ojg44NSpU/D29kZ0dDQKCwvf0HWl9yZFIpFW7Kdoe6nOBqOWJOvYllM5OTmora2FoaGhiqRb84DLDkbazA7bgs6Q2d8DpVKJY8eO4ciRI9i+fTvee+89pjQukUi0LtjQEU6fPo3PP/8cz58/x/vvv4+xY8fiwoULXZ5Yd3FxgZWVFe7du4eSkhK89dZbJCBqAZIp9nI2bNiAw4cPw8LCAleuXMGgQYMwc+ZMrF+/HhMmTAAA+Pv7Y+vWrZ3ORHsz9+/fx9ixY994nKIolJaWMkM8YrGYcd/QhK5rQ0MDMjIyoKenBzc3N61d5NmSdVVVVZDJZMyQi1KpRHFxMVxdXXWSHdIycTKZDB4eHkzALi4uRkhICJycnBATE8O520p34Pjx47h58yaePXuGP/7xj0zPnNAhyPJ+X6CjDf8tW7agvr4eERERJChqGLauq1gsxv/+9z8YGBhAKBQygbI9Xdf2BLy1DUVRqKysxKNHjyCXy6Gvr6+yW0gPumhbtq0liymKonDy5El88803iI6OxvTp0zmTj1uzZg3OnDkDQ0NDuLq64vvvv8eAAQMA6GYRXyaTwcvLC3K5XG2xkD4MCYqEX3n69ClmzJiBtLQ0BAcHY9KkSYz1k5ubG65evdqp8ml3u1h0J9rTdRWJRPD29mZ0XR8/foyKigqYm5tDIBDoZLeOlqlzdXVlRAboIRc6m2y+W2hhYaGxTFahUCArKwt1dXXw9PRkMu3nz58jLCwM/fr1Q2xsrMpKBxckJydj8uTJ0NfXx7p16wAAW7duRXp6OubPn8/0B7lcxP/kk08wYMAAREdHa/21ehkkKPZ1srKyMHz4cADAzp07ce3aNSQkJODcuXOIi4tjBm1WrlwJiUTSqefujheL7oxSqUROTg4zxHPnzh00NDTAwsICT58+RUxMDKZMmcL590kmkyEjIwM8Hg9ubm7t9g7p3UJ6iIe2YuqKpFtlZSUyMjLg6OgIBwcHJjs8d+4cNm/ejE2bNmHOnDk6t3g6ffo0EhIS8K9//Uvj2scdRalUQigU4uTJk8zvNqHDkEGbvs769euRmZkJPp+PIUOGYM+ePQCAGTNm4Pz58xAIBMxka2f53e9+x/x93LhxSEhIAEBUO1qDz+dDIBBAIBBg8eLFkEqlWL58ORwdHTFt2jQcPXoUmzZtYiTb6LUQKysrrQUDepiHnR22R/PdwpYk3dpz0qChBeCrqqowevRoZt2lsrIS69atw8uXL/Hzzz/rpK/ZEvHx8Zg3bx4A3Szip6enY+bMmfjDH/5AAqIWIUGxF/Pjjz+2+LgmJlvZ6Ppi0RP573//ix07dqjsh9J9xZSUFKSkpGDHjh2oqqqCp6cnU3b18vLq8iQonR3y+XyVJXh1aEnSje2kUVxcrGJgTGeUdXV1kEqlsLOzg7e3N5MdXr16FV999RW+/PJLLF68mBM914705SMjI6Gvr4+FCxdq/Tyt4enpiZycHJ29fl+BBEVCq/SUi0VPJCgo6I3HeDweHBwc8OGHHzK+nHK5HA8fPkRKSgr27t2Lhw8fwtTUVG1dVzo7FAgEjFO9pjEwMICVlRWzLkBL1lVWVqK0tBRpaWlobGyElZUV8vPzUVJSAjc3N0RERCA3Nxfnzp1T0f3VNu0t4h88eBBnz57FpUuXmO8zUZzqvZCgSGgVXV4siP1OE/QUq1AoVNF1lUgkKrqurq6uTDbZkq6rTCaDVCqFnp5el7PDzkJL1ikUCuTn58PJyQlOTk6ora3FxYsXER8fz+x/fvDBB0hNTYWhoWG3KJsmJSVh27ZtuHbtmoqaUUBAABYsWICwsDAUFRUhKyuLU0UdgvYggzYEtUhKSkJYWBiuXbumknFoymdSKpWCz+cjODhYRRWEDPK8CVvXVSwW4969e6AoCkKhEN7e3nj27Blu376NuLg4rWWHbUFRFJ48eYLS0lJ4eHgw+4X19fWIjIxEamoq9u/fD0tLS0gkEty6dQtDhw7Fxx9/zPlZmyMQCNDQ0MBkvePGjWN685GRkYiPj4e+vj62b9+O6dOn6/KohPYh06cE7cHVxaK5VJaupv56ErSua3JyMr7++msoFAoYGxvDxsZGRdfVwsJC6xOdtbW1SE9Ph6WlJYYNG8b0CO/du4eVK1ciMDAQYWFhffqmhsAZZPqUoD0eP37c6r9pwmeyNcggT/vweDzw+XxER0dj27ZtmDFjBqPreuvWLVy+fBlbt25FXV0dRo0axQgMaFLXlaIo5Ofno7i4GO7u7rCwsADQ1CP929/+hkuXLuHQoUMYNWqURl6vIxCBfEJHIEGRoDOI/Y72MDExwc2bN5kgR+u6uri4IDAwEECTpNz9+/cZXdeMjAxYWFgwmaSvr69auq6vXr1Ceno6zM3N4ePjw2SB6enpCAkJwbRp03D9+nXOdUrXrFmDzZs3A2gSyP/LX/6CPXv24KeffkJWVhaysrIgFovx6aefEoH8PgwJigSdoQv7HTZJSUkIDQ2FQqFAUFAQ1q9fr9bzdFfay/qMjIzg5+cHPz8/AG/quu7du7dTuq70Skl+fj7c3NxgaWkJoEmtZteuXUhISMCePXt0JidIBPIJHYEERUKPQlNTfwqFAitWrMDFixfh6OgIkUiEgIAAeHp6auHUPQMejwdbW1sEBAQgICAAwK+6rrdu3cKRI0ewevXqFnVdnzx5gqtXr2L8+PEQiURMdpibm4vPPvsMvr6+uHHjhlZNkDtCc4F8oKkkz14BoUvyJCj2TUhQJHRLtGW/QyORSCAQCDBs2DAAQGBgIBITE/t0UGwJfX19jB49GqNHj8by5ctBURSqq6uZlZCEhAQ8ePAASqUSAQEBGDp0KBwcHGBqaor4+HjEx8cjNjYW7777Lifnba8kHxkZicjISGzZsgVxcXGIiIjg5FyEngOZPiX0SRISEpCUlIT9+/cDAI4cOQKxWIy4uDgd3hm0VgAABa1JREFUn6zn8OLFCyxfvhzm5uZYuXIl0tLSGF3X7Oxs/P73v0dsbGy3Mf1lo2mBfEKPgEyfEggE7WFkZISgoCBMmzYNAPDWW29h8eLFAJoMggcOHMiJTFtHYQvkJyYmwt3dHUBTST4uLg6BgYEQi8WwsLAgAbEPQ4IioU9CZLq6jqmpKRMQm8O1H2RH0KZAPqH3QMqnBK1y+/ZtLFu2DBKJBAqFAr6+vjh+/Din+2kt0djYiBEjRuDSpUtwcHCASCTC0aNHMXLkyC4979KlS3H27FnY2NggLS0NQFOZcd68ecjLy4OLiwtOnDjBTGYSCATO6FD5tPvUNgi9Enqq889//jPWrl2LRYsW6TwgAk0DJHFxcZg6dSo8PDwwd+7cLgdEAPjoo4+QlJSk8lh0dDQjd+fv70/MYTVITEwMeDweysrKADSthaxcuRICgQCjR4/G3bt3dXxCQk+DZIoErSOTySASiWBsbIybN2/2ekmvvLw8zJw5k8kU2YMbxcXFmDRpEjIzM3V8yp5Pfn4+goKCkJGRgTt37sDa2hrnz5/Hzp07GQPt0NBQsohPoCGZIqF7UF5ejpqaGrx8+RL19fW6Pg7nlJSUMIMbdnZ2KCkp0fGJegerVq3Ctm3bVBR3WlvEJxA6CgmKBK0THByMzZs3Y+HChVi3bp2uj6NTeDye1kW4+wKJiYlwcHDAmDFjVB5vbRGfQOgoZPqUoFUOHz4MAwMDLFiwAAqFAm+//TYuX76MyZMn6/ponGFra8vIhhUXF8PGxkbt58rPz8eSJUtQUlICHo+H5cuXIzQ0tFcO87S1iB8VFYXk5GQdnIrQ2yE9RQJBwzTvKa5ZswZWVlZYv349oqOj8eLFC2zbtk2t5y4uLkZxcTGEQiFevnwJb29v/Pvf/8bBgwcxcOBA5jUqKiqwdetWTb6tbsPDhw/h7+/PmP4WFBRg8ODBkEgk2LRpE1nEJ7RGx0o0FEV15g+BQGiDwMBAys7OjtLX16ccHByo/fv3U2VlZdTkyZMpgUBA+fv7U+Xl5Rp7vYCAACo5OZkaMWIEVVRURFEURRUVFVEjRozQ2Gt0d4YMGUI9f/6coiiKOnv2LDVt2jRKqVRSKSkplEgk0vHpCN2IDsU5kikSCD2UvLw8vPvuu0hLS4OzszMqKysBNN3oWlpaMl/3dlxcXJCamgpra2tQFIWQkBAkJSUxi/i6cuUgdDs6lCmSoEgg9EBqamrwm9/8Bhs2bMDs2bMxYMAAlSBoaWmJiooKHZ6QQOh2kJUMAqE3IpfLMWfOHCxcuBCzZ88G8OswD4AuD/MQCH0ZEhQJhB4ERVFYtmwZPDw8EBYWxjweEBCAQ4cOAQAOHTqEWbNmqf0a9fX1jKHwyJEjsWnTJgBN3oh+fn4QCASYN28eZDJZ194MgdANIeVTAqEHcePGDUycOBFeXl6MA0VUVBT8/Pwwd+5cPH36FEOGDMGJEycwcOBAtV6DoijU1tbCzMwMcrkcEyZMwI4dO/DNN99g9uzZCAwMxCeffIIxY8bg008/1eTbIxC0CekpEgiErlFXV4cJEyZg9+7deP/99/Hs2TPo6+sjJSUF4eHhuHDhgq6PSCB0FNJTJBAI6qFQKDB27FjY2Njgt7/9LVxdXTFgwADo6zfpfRClGEJvpbOZIoFA6EPweLwBAE4D2AjgIEVRgtePOwH4iaIo3VueEAgahGSKBAKhVSiKqgRwBcB4AAN4PB4tDekIgKSKhF4HCYoEAkEFHo836HWGCB6P1w/AbwFI0RQcP3z93/4IIFE3JyQQtAcpnxIIBBV4PN5oAIcA6KHpxvkERVF/4fF4wwD8AGAggHsAFlEU1aC7kxIImocERQKBQCAQXkPKpwQCgUAgvIYERQKBQCAQXkOCIoFAIBAIryFBkUAgEAiE15CgSCAQCATCa0hQJBAIBALhNf8H4sDkJ0h2tm8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "prody.showProtein(prot)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "With the latest version of ProDy, `py3Dmol` will be used if it has been imported." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "@> ProDy v1.10.11 is available, you are using 1.10.10.\n" ] } ], "source": [ "prody.checkUpdates()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import py3Dmol\n", "prody.showProtein(prot)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Configuration\n", "\n", "A number of keyword arguments are supported." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The background color is configurable." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,backgroundColor=\"black\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The height and width of the viewer can be changed." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": true }, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot, height=400, width=200,backgroundcolor='gray')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The viewer can be zoomed and centered on a given [AtomSelectionSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomSelectionSpec)." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot, zoomto={'resn':'ATP'})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Styles\n", "\n", "Custom styles can be applied by providing [AtomSelectionSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomSelectionSpec) and [AtomStyleSpec](http://3dmol.csb.pitt.edu/doc/types.html#AtomStyleSpec) objects.\n", "\n", "A single style may be applied." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,style={'cartoon':{'colorscheme':'chain'},'stick':{'colorscheme':'chain'}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For more complex styling, generate and manipulate the py3Dmol view object directly." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "view = prody.view3D(prot)\n", "view.setStyle({'cartoon':{'colorscheme':'chain'}})\n", "view.setStyle({'hetflag':True},{'stick':{'colorscheme':'magentaCarbon'}})\n", "view.setStyle({'bonds': 0},{'sphere':{'radius':0.5}})\n", "view.addLabel('Chain A',{'inFront':True,'showBackground':False,'fontColor':'blue'},{'chain':'A'})\n", "view.addLabel('Chain P',{'inFront':True,'showBackground':False,'fontColor':'green'},{'chain':'P'})\n", "\n", "view.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# GNM/ANM\n", "\n", "Squared fluctuations (or any other alpha-carbon length numerical vector) can be visualized in a RWB colorscheme. Red values are larger." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "@> Kirchhoff was built in 0.04s.\n", "@> 20 modes were calculated in 0.23s.\n" ] } ], "source": [ "gnm, sel = prody.calcGNM(prot,cutoff=7.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Individual modes may be visualized." ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,flucts=prody.calcSqFlucts(gnm[0]))" ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "Styles can still be applied." ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,flucts=gnm[1].getArray(),style={'cartoon':{}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Motions\n", "\n", "It is also possible to provide an array of vectors (one per a calpha) that can be used to animate motions. The animate kwarg takes the 3Dmol.js animate options." ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "@> Hessian was built in 1.27s.\n", "@> 20 modes were calculated in 0.60s.\n" ] } ], "source": [ "anm, sel = prody.calcANM(prot)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "scrolled": false }, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,vecs=anm[0].getArray(),\n", " flucts=prody.calcSqFlucts(anm[0]),\n", " style={'cartoon':{}})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Multiple Structures" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "application/3dmoljs_load.v0": "
\n

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n jupyter labextension install jupyterlab_3dmol

\n
\n", "text/html": [ "
\n", "

You appear to be running in JupyterLab (or JavaScript failed to load for some other reason). You need to install the 3dmol extension:
\n", " jupyter labextension install jupyterlab_3dmol

\n", "
\n", "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.showProtein(prot,prot,flucts=[gnm[1].getArray(),prody.calcSqFlucts(anm[0])],vecs=[anm[0].getArray(),[]])" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'1.10.10'" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prody.__version__" ] }, { "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.9" } }, "nbformat": 4, "nbformat_minor": 1 }