{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Multivairate Model" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {"nbsphinx": "hidden"}, "outputs": [], "source": [ "import warnings\n", "warnings.filterwarnings(action='ignore')\n", "\n", "# for python 2 and 3 compatibility\n", "# from builtins import str\n", "# try:\n", "# from pathlib import Path\n", "# except:\n", "# from pathlib2 import Path\n", "#--------------------------------------------\n", "import sys\n", "ppath = \"../..\"\n", "\n", "if ppath not in sys.path:\n", " sys.path.append(ppath)\n", "#--------------------------------------------" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from __future__ import print_function, division, unicode_literals\n", "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "\n", "plt.style.use(['seaborn-paper', 'seaborn-whitegrid'])\n", "plt.rcParams['font.sans-serif']=['SimHei']\n", "plt.rcParams['axes.unicode_minus']=False\n", "\n", "import numpy as np\n", "\n", "import pygeopressure as ppp" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1. Calculate optimized multivariate model coefficients" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Create survey with the example survey `CUG`:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# set to the directory on your computer\n", "SURVEY_FOLDER = \"M:/CUG_depth\"\n", "\n", "survey = ppp.Survey(SURVEY_FOLDER)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Retrieve well `CUG1`:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "well_cug1 = survey.wells['CUG1']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Get Velocity, Shale volume and Porosity logs:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "vel_log = well_cug1.get_log(\"Velocity\")\n", "por_log = well_cug1.get_log(\"Porosity\")\n", "vsh_log = well_cug1.get_log(\"Shale_Volume\")\n", "\n", "obp_log = well_cug1.get_log(\"Overburden_Pressure\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Calculate optimized multivariate model parameters:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "a0, a1, a2, a3 = ppp.optimize_multivaraite(\n", " well=well_cug1, \n", " obp_log=obp_log, \n", " vel_log=vel_log,\n", " por_log=por_log, \n", " vsh_log=vsh_log, \n", " B=well_cug1.params['bowers']['B'],\n", " upper=1500, lower=3500)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "View velocity, porosity, shale volume and effecive pressure used for optimization, and Velocity predicted by the optimized model (blue line):" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEPCAYAAACgFqixAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXl8E3Xexz+Ts2nSNkkL9G5paQdo5QY55T5EBBUvvPVZvB9XUYTddWFdHnxWxWN3H13xRsUDcVFRRLkXSrmPQluGQo/Q+76b5prnjzBDJleTNkeTzPv1yquZmd/MfDNNft/f7/v7HgRN0+Dh4eHhCW0E/haAh4eHh8f/8MqAh4eHh4dXBjw8PDw8vDLg4eHh4QGvDHh4eHh4wCsDHh4eHh4AIm9dmCTJ+wGoABwFsBZAKYCPAWgAbALQCmA5gDDLbYqiWr0lEw8PDw+PfQhvxBmQJDkCwDMA8gB0AzhEUVT+1WOvw9z5SwFMApBiuU1R1P95XCAeHh4eHqd4xUxEUVQegM8AEABGAVhBkuQ2kiQTAMRRFHUewDkAgwHEWm3z8PDw8PgYr5mJYFYENICVFEW1kyQ5AcBDALquHhcCUADQWmzLHV3s5MmTfKh0D4wdO5bozXn8s3WN3jxf/tm6Bv/d9R6uPltvKgPArBCWAfgAQDQAI4DzJEmOhdkspAHQbrF9xdnFxo4d6/BYYWEhhg0b5vR4TEwM8vLycOjQIZvjAoEATz/9NFQqVa+v78/jJ0+edHjMFZw9W3dl6e/te3NOX56vO8/Wkt58rr6c5697evu7a0+2n376CeXl5Xj88ccdtnHlOv29jTvP1tvKgAagJUnyewAGAI8B0AN4B+ZZwaMwrylYbnuUzs5OvP76607bvPjii7h8+TKUSqWnb8/Dw9MPUSqVOH/+PGiaBkH0alISdHhNGVAUdQDAgaubn1sdvr+HbY8RHh6OtWvXoqCgAJmZmVi/fr1Nm9deew0A8N1339kcW7NmDf9l4eEJMpRKJbq7u6HVaiGTyfwtTr/A2zODfgNBEBCJRPjv//5vNDY2wmAwwGg0sq/y8nIMGDAARqMRu3fvBgBERUVhz549AICGhgZUVFQ4vH5DQwPkcjmSk5N98nl4eHh6T1RUFACgubmZVwZXCRllwKBWq6FWqwEAJpMJ+/btg9FoRHt7O6RSKQAgKysL+fn5aGlpQU5Ojt3rjBw5EjExMWBcc8ViMUSikHucPDwBCWMSbm5uRlxcnJ+l6R+EbO9F0zTWrVvX6/PPnj0LALjxxhsxYcIEFBYWIj4+3lPi8fDweBGFQgGhUIiWlhZ/i9JvCAll0NXVhS1btgAAJBIJlEolBAIBEhISoFAoAACtra2IjIzEjBkzEBsba3ONvnhL8PDw9C8IgoBSqURzc7O/Rek3hIQyCAsLY9/rdDrU1tZyjj/yyCNob2/nO3senhCCVwZcQiZRnUgkgkKhgEBg+5E//vhjbNmyBZcvX/aDZDw8PP4gKiqKVwYWBP3MoKOjAxs2bAAAaLVamEwm9tjvfvc7iMViEAQBjUaD9PR0f4nJw8PjY5RKJfLz8/0tRr8h6JWBXC7H2rVrWZt/bW0t8vLykJOTgw8//JDT9ueff8bEiRNBEAQIgsDkyZMRHh7uJ8l5eHi8iWWsAU8ImIlomkZVVRWam5tRUVEBrVaLtLQ0h21zc3Nx+PBhSKVSPtiMhyeIsXQv5QmBmUFHRwfef/99REVF4dy5c2hvb0dbWxvS0tIgFAohFAohEolQUVGBpqYmyOVyqFQqlJSUoLS0FIDZ82DQoEH8AjMPTxBhGXjGEwLKQKFQcMxEDC0tLTAajRAKhaBpGn//+98BmJVHR0cHYmJiIJebk6jSNI2amhrs3r0bEydOZN1ReXh4ApeIiAgIBAI0NzeziiGUCXplYI1er8drr70Gg8Fg93hqaioIgkBCQgK7j6ZpNDQ0QCqVQigU+kpUHh4eL2IZa8ArgxBUBmKxGH/6059gMpnYxSOdToeioiIcO3aMNQ2VlJSw50RERGD+/PnIysryk9Q8PDzeQKlU8lHIVwk5ZcAgEAggk8mwa9cunDt3zuFMQSwWY8WKFSgsLPSxhDw8PN4mKioKVVVV/hajXxCyygAAamtr0dHRwSqC6dOnQ6fTITc3F9dffz3vWsrDE+QolUp+oHeVkFYGOp0OFy9eZLcPHDgAiUSCxMREVFZW4ttvv4XJZML8+fP9KCUPD4+3UCqVrKk41AlJZVBQUIBvv/3W7jGdTofy8nLMmjWLXUweNGgQLl265GMpeXh4vA1T5rajo8PPkvifoA86s8fw4cOxdu1ajBo1ymGbvXv34uOPP0ZYWBjEYrEPpePh4fEVvDK4RkgqA4bJkyez76Ojo22O//nPf0ZMTIwvReLh4fEhcrkcIpEI7e3t/hbF74SkmYhhwIABWLNmDc6dO4dt27bZHF+3bh3kcjmEQiGMRiP27NmDyZMnY8yYMX6QloeHx9MQBAGVSsXPDBDiygAwl77Mzc11eNzyS9LR0YHt27dDpVJh8ODBvhCPh4fHy6hUKrS2tvpbDL8T0sqgtrYW//rXv1xqS5IkWzt5wIAB3hSLh4fHhyiVSj7WACG+ZjBw4ECsWbPGpc6doijk5uYiNzfXboEcHh6ewIQxE1nWOglFQnpmAAAGA4HOzidRXl6P9PQYCIUUCgouQKGIRHJyKghCCECIxsZmREcPgFgsxTffhMNoBEaOBMaP9/cn4AlUWlpacOLECdA0bfMCgIaGBhQXF8NkMsFkMsFoNHJe2dnZUKlUUKlUiIiI8POnCVxUKhVMJhPa2tpCOkdRyCuDrVu1WLkyDADjNURefVmTYLPn1luBf//bi8L5kP3797vctq6uDjU1NV5tv2XLFpfb94akpCQUFhZCJBJh4cKFEIl891Po7OzE8ePH0djYiLy8vF5fx53Yl8zMTCxbtqzX9wpmGPfSUE9YF9LKoLq6GgbDLqxdG4NXXpkLvd7x49ix4xJmzBgCgQDsK5gSmM6YMcPlttbpwF1pr9FocOTIkV5I5h06OzthMplAEAROnToFgiDYUbnRaMTZs2chl8sxdepUh8WQXOX06dM4e/YsaJrG8OHD0dnZiWPHjkEgEEAul7OV9dra2np1/XvvvRdDhgxht939/4Q6TJGbpqYmpKSk+Fka/xHSykCpVGLy5MmIisrH9On7sHv3XACARELjttsIvP02MGiQuW1hoR4ymR+FDXDmz5/vclqPnjozrVaLgoIC1pxSWVmJ2NhYdtvS1MJgua3T6TgzofLycof3SkpK6rMyuHjxIsrKygAAGo3GbhuhUMg+H4IgIBAIUFNTg7i4OCQlJWHgwIF9koHHMRKJBGFhYWhqavK3KH4lpJVBWFgY0tLS8MUXXyA393l2v9EIfP454EPLAY8b6HQ6XLlyhd1uaWmx6fwB84i8Lzz88MOsB1lfmD9/PrKzs7F161aHbYxGI3799Vd2Ozs7Gx0dHTAajYiMjOSVgZeRy+W8MvC3AP6Gpmk899xz2L6dwsmTYwEARiOBZcu24r33ZtqNTObxL5GRkViyZAm7bT2TqK6uxsaNG/t0j7i4OEilUo9UtVMqlaiurkZsbCyqq6t7bP/nP/8ZAoGAN/f4ELlcHvLlL0NeGRw7dgy//vorZs+WscoAALZuvR3MQC4j4yLuvXcLRCIRjEYj5s+fj+uvv95PEnuHdevWgSAIAOD8tX4B5lKidXV1Nu2sz2VG6zU1NWhsbAQAjinH3l/g2oIzTdPs9Zxhb4F6+vTpnO0DBw70eJ1FixZBIpEgPDwcMpnMowOB1NRUREZGoru7G93d3WhqasJvv/3GaRMREYExY8bgP//5DwDu52KeBfPX3jNnqKqqcmkBXq1WY8KECUH3Xe4NcrmcM9sMRUJKGeTk5KC6uhoSiYR112O8OQYPjsSaNTnYvn0ITp8exDmvvd08OszKysLChQshFArZxUdXOqtAwB0f666uLuzbt88rcsTExECv16O+vh4EQWDRokWIioridHwSiYTj/ePKCNp6gdzXo+6wsDDEx8dz9k2aNMluW6PRCJPJhAsXLiAzM5PjclpZWYnNmzd7RKbGxkbs3LmTjZ0hCAJ6vR6JiYkh56qqUCjQ3t4OvV4fsokpQ0oZMCO9sLAwdpEuIyMDgHkRcdIkBf7971oAXGVw223mvEVnz57F+fPnQdM023kOHz4cd9xxh+8+hJd46qmnWOX20UcfobOz0y9y1NfXAwBbivDDDz/s8RySJNkCRQRBQC6XB3S6kB9//LFPLqeuMHv2bIwYMYL9LjPK5tKlS5DL5V69d3+E+cxNTU0huz4TUspg6NChDo8JhUIMGzYMY8fqcP4899iKFU+huzu47beWHcDTTz/NOcaYcD777DO34gWcMXfuXBtTB/O+tLQU7e3tDs0glhAEgc7OTrZIEU3TUKvVAa0MZsyYAaPRiMbGRoSHh7PPxmAwID09HWKxmO3ArU1N9lCr1XjwwQdtZrKWJqfw8HCIRCLU1taGZIQ98/1vbm7mlQGPmU8/leD77wHLGtmbNwO33+4/mXzBP/7xD2i1Wp/db9euXb06TyKRYOXKlW6biQIJlUqF22+/nf1cGzduZBeeS0pK3L5eY2Mj3nrrLbvHFAoFCILAkiVLkJ6e3ie5AxmZTAaBQBDSHkUhpwza2trQ3d2NXbuK8cknLSAIGgKBCQKBCVptFfbunWVzzi+/5GLp0uCOTPSlIrAkKysL46/m9GBmIGVlZUhOTmbbWMcNlJaWco5VVVVBJBIhLi7OI94//Y3HHnsMV65cgVar5cyWDAYDvv766z5de9CgQSAIAkeOHMHRo0fR3t6O06dPc2YQcrkcN998c9Csj9lDIBBAqVTyyiCUePPNNwEA3313K86dm+DSObm5k/DZZ6fxt795U7LQJD8/H/n5+Tb7mXUdwDxSlkqlnM7I8n1HRweqqqowefLkoFQGAHDu3Dk0NTVxvK86OzvZWhu95fLlyy6162vMhiMkEgnmzZvnlWu7i0ql4pVBKLF27Vp89NFHmDdvF8LDu1BZGQejUYjKStvcQ5a8+upo7NpVhP/5n0u48cYbfSRt4BEbG8uxOdM0Da1Wi7CwMHYb4K4RVFVV2QSNySzCvadNm8amW7BWCMHuj8/MiubPn8++r6iogEajweXLl/ukCDzJzTff3OuiTydPnvSwNL1DpVI5jBAPBUJOGQDA4sWLcfnyhzh61D3/aq02FqNGBefI01O4ElTlCpYjtG+++QYAbJRMcnIyHnroIY/cr7+g0+nw6quvclx9LWM2+ivbt2/H9u3be2ynVqvx9NNP90uTk1KpZHNI9Uf5vE1IKoMBAwbglVf+gFtvBc6dM++rrq6ERiPHxo3ctYH33weWL7dcpAxO/+u1a9e63NaVkbhGo8FXX30FnU4HwH4AG+PX7urodtq0aTb79u/f73JWVK1Wi6NHj3L2TZ8+HaNGjYJEIrGJX/AmBoMBRqMRZWVl+Oqrr1hZgGufs66ujlNrw17gXHR0NLKysiAUCiESiSAUClFXV4f09HQMHTrU7U6tt7OsYJidqVQq6PV6dHR0BK250RkhqQwYxo+/Vo+gsLAFBw4cBXArAOCHH5qweLHKf8IFOMnJyVi1ahWA3mU59UZ7g8HAug1WVVUhKSkJ9fX1+Pvf/862ueeeezBgwAA2k6Wnqa2tRU1NDf5tkft87ty5iI2NtUmIZ/25XM0sGwwdsz9gUlk3NTXxyiCUKSgoQHn5BTDKYMkS8xfjo4/24sEHpzs5k8dfNDc3Y9euXWhtbUV+fn6P6TEsaWlpYWckkZGRbA3cL7/8ktPu/vvv73PW0mPHjuGXX35xeFyj0aCsrIxN8U0QBJYuXdqne/K4j2Vdg6SkJD9L43t4ZXCV4cOHY+jQkVi/nrtfp0sJySCc/oa9amBCoRBDhgxBRUUFYmNj2RQjllG1ltG1JpMJZWVlThcJVSoVnnjiCY+mJJgwYQImTDB7rul0Ovzv//4v5zhFUQCAsWPHYsCAARAIBOjq6kJXVxc6OjrYhXJGyVm+D6aUKP4mLCwspFNZ88oAwI4dO3D8+HEYDEIAL3GOPfFEOp54Arj77lioVOb6Bn/8IxBs6Utefvllv907OjqaM4Jn6g1YdnRMBygQCDjvmQjktrY29phQKOS8t3zFx8cjLCwMnZ2dUKvVEAqFIAgCERERMJlMuO6667yam0YikeCRRx4BRVEQiURsmmqj0YiIiAicP38eYWFhqKioQFNTE+sxxOQrsn5vqQBpmvZoR2YvwC/YCWX30tD5Lzth4cKFGDx4MLq6uvA//2O/zddfX1s/GD0aWLzYR8L5CGYB2TLAy957mqZBURSbQM2Vc4qKijBkyBB2/969e3HeIudHQ0NDn+Wvqqpy+xzLojZ//vOffTbKTkpKQmJiIsrKyrBp0yanbSdNmgSRSISUlBQ21sLeDMHy/eXLl5GZmel0JsHsA7iL+6G+3sArAx4AcOon/Yc/VGPx4ljI5cB11/lQKB/hyZlBTIy5njTT2eh0Ohw/fpzT2cbFxXHaWB7r6upic/IA5s7T0ejUldTUzmDusW7dOgDAyJEjccstt/Tpmj3xwQcfoLKyEgRB9Bg0lpubC8CcUkIqlaK2thaAfaXLzBC0Wi1+/vlnr34GwDfPyteoVCpUVFT4Wwy/wCsDKyIiAHulaOfMacPEibG+F8hHTJgwAceOHfPItZjMo33BcrbQ10CguLg4CK0KVnd2dkImk4Gmadx000026aW9yfLly9n3+fn5yMnJsZmRVFZWcrbtRWl7ktmzZ2Ps2LEoLi7GsGHDXEoSGIwolUrWucD6OxPs8MrgKvX19Xj55ZfR1mbf3372bHOq66QkIyhKGHT1kCmKQmRkJLttnQ+oLwXb/U1PJqQPPviAs/3ggw8iNTXVK7L88ssvHlO6vYGZtTEw/+MzZ84gOTmZXWsJVSw9ikKtyiGvDK7CjMTGjz+O48fHO2x35YoQ/SQDgEd59tlnXW7ryK7c1taGkpISGI1GdHR0YN++fZg+fToaGhrcSgtcW1trt731KJXZrqmpYROuWR/fuXOny/dl2LRpE1566SWPjwzb29vtKoLFixfbmMEIgkBFRQUSEhLcHp2Xl5cjISEBsbGxNp1/TxQWFrrVPtiwjDXglUGIIpFIAADV1YNsjsXEtOCHHyowefJwX4sVUERERGDEiBEAzCPOwYMH4+TJk24VaklMTERnZyfHTMR40CQnJ9sNBtLr9eju7rbZzxS8cReVSuWV0bFEIrGbWuLHH3+0aTtlyhS0tbWhsbERAoEANE1Do9GwI3dr7yrLbZqmkZ2d7XH5QwGmql4oLiJ7XBmQJDkBwGtXNx8AsBFAK4DlAMIAbHK0TVFUq6flcQRN03jvvffYBTmGO+74Fm+++TxnX319FKZNU2Dhwp/w8ssJGD16VEjYUtvb29kfheXnbWho4HjiMJ3b4cOHceHChT7d0/K61tTV1eH2229HVlYWZ78zD5g5c+bY3e8PrxmJRII1a9a43L6goAAkScJkMqGzsxMURcFoNLqU/8k6+duDDz5o40HEQBAEZDIZ1Gq1y7IFK0KhEFFRUbwy8BAZABYAmAfgDIBpMHf6DwBIAbASgNTB9v95QR67EASBJ554AgCwZ88eHDp0CAAQEdFut73JJMRPPy3CxIlnMXq0r6T0Hf6MM+iJhQsXIirKnDOKIAhcunSJ06nV1tZCJpPZdZW0tw8w24Rramo46yBSqdQnkac0TaOkpIRNRkcQBBQKBSdeoKurC99++y3nvJiYGIwePRrXXXcdm720uLjY7qzIGkcurFKpFNHR0UhMTOSz8V5FrVbzysATUBS1GQBIkhwH4CRFUfkkSUoA3ANgIEVR551s+4XZs2cjPj4ew4YNQ319G5z1iy+9NBJDh34XdOkCXE1UR9M0NmzY4NMayTt27PDIdRhzikKhAE3TOHTokM3nWL16NaRSqUfu54jy8nJ8/vnnbp9XX1/vsELchAkTMG/ePDZWYPhws0kzFGawnkapVDqdoQYrXlkzIEkyA0AigNKru4QAFAC0TrZ7rMLtbHFLq9V65LjZzGxrPhg0SIsbb+zA7NltGD58uM21PHX/3h7vK+5ce86cOWx9Ams6OjpQW1vLsYsbDAZ2gfT06dN+y8HPROs2Nzez+wYOHAiCIBAfH4+MjAwUFxd7/L72nu2dd96Js2fPoqamhiNPbzl27BiOHTsGtVqNqVOn9tpc19vvmbe/n87o6b6uyGbZxmAwoLGxEQUFBRxl6u51+kMbd/DGmoEM5jWDhwE8RJLkWJjNQBoA7U62r/R0bWc23p5swK4eNxpNdo9v3FiJJUvSANj3MPDU/Xt7vK8FQjyVJbSiosImaKe5uZnNAnqdCxF7lu2dwSic5uZmREVF2Q3C6unHUltbi5kzZyIrK8up90hfnq+jZxUVFYW8vDybtNpCoRBZWVloaWmBSqVyKzJ6+PDh0Ol0vV4P8UcKa29/d12RzbINTdPIy8uzcVhw9zr9oY07z9YbM4M/ABgAs0I4COBZmGcBjwLoBvCOk22/IxRyvUiGDdNi+HADPv1Ujd9+A155BYgK7nLIfSIhIQEJCdyqcf5MYd3W1saWOnXEvn37sG/fPgBAWloa7rvvPq+ZV/R6PV555RWnbWbPno1Jkyb1qWPm6T2hmsraG2sGawBYukxYG0fv72Hbb9A0fdUd8VqissLCMFj+tgYPBl54wfeyBSq7d+9GTk6Ov8VwmeLiYvz1r3/FsmXLoFKpOMVl+gIzY3KWyprht99+w2+//Wazf/HixRgdjN4L/QxLZRBKqaz5OAOYc+ds3brVotQgdzH1hRcAmawaCxbEYtIk38sXyMyZMwcJCQmcKfi7777b55QV1iN3Ji7AnheRK9421nz11VcYO3Ysbrrppj7PEkwmE37++We3kumJRCKOdxFgjkewF5MAmGM8nnvuOX7B2AOEhYVBJpOhsbHR36L4FF4ZwOz/vWTJEnz//fd2a81u2AAAsXjzzXasXPkGlEolnnnmGf6H1wsIgsBTTz3ltI21eYSmaeTm5nI8aaz/T95YlM7Ly8PkyZP77H9fVFRkE89iSWJiIpYtW2aTUZSiKDZPkHVcAI93UavVHlnYDyR4ZXAVsViMNWvW4Oeff8bcub9h1655Nm3GjGkCSa7BHXcQ4H+PvoMgCEyePBmTJ0+2e7ywsBBDhw4FwFUSlgvKlvuMRiO+++47XLp0yel9V69e7ZFIZJIk8eyzz+KNN96wew978RDMe8sZD4/vUKlU/Mwg1ImMjERS0kW7xw4eTMLBg8DmzT/hq69mhNTikjWHDx926PMe6CQmJuKWW27xaEoKhUKBhx9+GGfPnkVnZyfkcjnUajVycnJw8OBBj92HYfr06S7XTOaxRaVSobS01N9i+BReGVgxdepUyGQybNpkgtFovzMYPfoU3niD67K1ONiq3Vixb98+/Oc///HKtVNSUjhZQouKimxSONuDcT/V6/UQi8WsO6mj7KpJSUkYP348rrvuOr+ko0hOTkZycjK7bRn1PXDgQDavEJOLSafTsVXXYmNjMWiQbd6s2tpau95DBw4c4NR6uPnmm9m0E8wsiscxKpUK7e3t7HcrFAhZZfDFF1/g8uXLnH0dHTK8/vqLAMY5PG/ChEYsXboEALcwi16v95qsvsCyY4qOjuZUxwLMmUG9RVlZGcrKynpsp1Ao2FoAls++qKiIrewFXAsuk0gkDoPj+gPWUd9GoxEGg4FTUS4jI8Om9vOWLVtcUpaWbN++nX2fnJxst6a0yWSCVqvF3r177daRdvQXMD9zgiCQlpaGe+7xWzIBj2HpUeROxt1AJmSVwX333ce+r6kBVq2qwaZNtiMvS+688xsMG3YBxcUjbOzRTOWuQOXee+9FW1sbxGIxBg0aZNNRANeCuUpKSpCSkmK3JKa99hqNBomJiQ7bW+4DzG6YzPO0fs5FRUWc7SNHjvSqbOawYcNQWVnJjsQLCwvR0dGBqVOnYvr06e4/QBehaRrff/89CgoKIJfLERERwaY+UCgUaG9vR0xMDMaMGYPa2lrodDqbXEvjx4+32Wf5qqioQFJSktOyl45epaWlSE9Pt6kz7egv87pw4QKnKE6gwzgN8MogiGE6Ha1Wi61bt8JoNOGRRx4C4FwRAMCWLXdh+fIPQBDXUjKvWbOGzQcTyGzevNnfInicpKQkpKWlsduWs4m6ujo2bblQKIRQKITBYHAp8rmvKBQKGAwGREVFITMzE5mZmQC4ik+n00Gv17O5k/Ly8oIqedqNN96IMWPG9FvlERERAaFQGFKLyCGnDNatW2eVN0cAgjCBpp0vFiYnlyE5WYMBA7guguvWrcPq1au9IitP37hy5QquXLmCp556yqbIi/WagSNPJU9DEATmzp2LuXPn9tiWGWBs2bLF22J5FalUajMLZALrpFIpZs2a5WcJbSEIAiqVKqgUcE+EnDJg8sl3d3ezQWb33EPhxx9/xIkTadi69Q675yUmVmD27L0cl1KpVIqVK1eGXK1UX7BgwQJcf/31LrX1x2Kwrxg2bJhLGWU1Gg0++eQTH0hkZtq0aZxBVX19PZvbybLTHzNmTI/V1vqam8hb8MogRPj66685RUKMRiPGjSvGhAmvw2QSYPVqboGbw4cn4/Bh8+jx7rvP46uvgquSFNPh6PV6fPTRR2yVMHt25e7ubrZ+APNiUnkwHQFBEKiuroZUKoVSqURUVBSbs4imaezfv9+pPDt37uxVyUpr7r//fhv7NmAulCOTydgRa1JSks+9Rjo7O1FVVQWDwQCDwYCTJ0+ipKSk19e77rrroFarUVdXZzeNBkEQaG5uxsWLFyESiTilNpln0dXVBYlEwllUBsz/M5FIhOXLl9uU6ASCUyGrVCobJ5NgJmSVwYMPPsjZZgKXDh06hL179zo8TybT48knPZOvpj8iFovx+OOPO23jzg//0qVL2Lx5M2pqanDxov34DW/ibt2AgQMH4vHHH/drorrecu7cOafHCYLAqlWrsGTJEodtgrFT7y0qlQrNzc0wmUxLrSCDAAAgAElEQVReKYPa3whZZWBNZWWlS7bZhx+ug1QausFm7jJkyBDceeedbmchzc3NxZUrPWY19zhLly6F0Wi0O/r1BCaTCXPmzMHu3bs9el2JRIKVK1dCJBLxHbqHUKlUMBqNaGtrYyvtBTMhrwzq6+uxf/9+tLe3Y/To0RCJRDh+/LjD9u++G4933wXspDDi8SCPPPIIALOS/uCDD3x233/961+45557kJGR4ZXr79+/H0eOHOmx3VNPPYW6ujq+U/cjlu6lvDIIAZqamnDuXD5OnRqLlhbzP5ym7Xs3jBpVijvukGPp0uA1EzHU1dXh3XffdescS9fNaxlgvYdAIMC0adNQV1eHgoICj133yy+/BGCOHXn00b6X2TAYDNiwYYNb2VPfeeedPt/XGcw6T2946KGHkJKS4mGJ+h+Mm3FjYyMnQj5YCSllYDQa0dXVhbq6Onz22Wfsfo0mGT/9tKjH88+cScX583pIpR9i927zD5tZaNPpdPjtt98wduxYTJ061WufwVu87Kzws4vYSwznTUwmEyflgjsw9nNv1zsGzOmoe3I/rqqqwvvvv+91WRj68j/69NNPPScIgEWLev7t+QOxWIyIiIiQ8SgKKWXQ0NCATZs22SwOpqZqoFQ2oblZ1eM1DAYxBg2KR1SUkfVOGT16NMrLy5GRkQGZTOYt8UMGlUqF8PBwCIVCCAQC9q+1VxCzqNfa2ooBAwaw6SemTZvW4wJwYWGhTxSBq8TFxeHFF19EV1cXK/ulS5dQXV2NU6dO2bTvy8je3zz55JMcb6f+6loKhJZ7aUgpg4EDB2LlypV2j91xRyHS0lQwmYDCwiJ89dWPeOON523aKZVNuHz5OCZMmADAPMKKiopi69UGKoxrqfXo3jJYqLGxEQcPHkRLSwsiIiI4x+2lmmhtbXWax98RTU1Nbv8ANRoNxGIxJBIJjEajQ2VA07TDhHuTJk3CvHm2qcu9wddffw2Konp9fm8UQVZWFhYsWGDjKgxccyGmKApDhw7l7Lc8brnPkmBdtFapVKirq/O3GD4hpJSBMwgCYHKaFRefgULRbrddc7MKf/mLuePUaoF+NLj0CM5+8AMHDsTSpUt9VtP4hx9+wJkzZ5zKabnNmAEty2xa+sr3RG5uLnJzcwEATz/9NBtE5SkoisLXX39t99jgwYMhl8uxZMkS1pPJH8XpRSJRyGTpdAWVSuUXl2h/wCsDO9xxxx2Ijz8FazO6XN4JpTIMBEHgrrsI8L8Z77JkyRLWJ95e0RoGJkakra2NDZZj2r333ntu3XPRokUYO3ZsH6R2DEmSWLt2Lfbv32+z1sEEm50/fx5DhgzBvffe6xUZeNxDrVajq6sLWq3W36J4HV4ZWGEymfDGG2+go6MTwBjOsY6OcIwefRBz5uzFo48+CoEgsDOVegOTyYSPP/4Yzc3N0Ol07D5mJC8QCNj91kilUrfrFQ8ePBgymQx1dXW9XkwGgPT0dGRlZXl8NmCP+Ph4REREoK2tDRKJBDqdDlKpFBMnToTJZMLBgwd7XNBnRu9M5Lcz/vSnPzmspsbjHMtU1sEOrwys+Oijj9DZ2Qmj0X6+oUOHpmH69G5ER/OKwJq8vDxs27at1+e7owhUKhWkUini4uJgNBoRHR2N2267jT1uucAaHR2N+Ph4zvn+tHFnZmZixYoVdo/RNA2ZTIa8vDxotVqHdXjdqZ+xfv16zvbtt9+OrKws1wUOYRhl0NjYGPRRyLwysGL58uVXZwffOGyzfv0crF8P/OUv5tFbTEwMdDodBAIBSJL0lah+w2AwYP369XZHpCkpKRAIBGytAEezAGuGDBkCwNwZVlRU9DgtZ0Zqlvml3GHhwoW9Os8b7NmzB4cOHfL4dZVKJYxGI5YtW4bw8HDWA0sgEECr1XK2GU8tHi7h4eGQSCRoamryyazRn/DKwAq9Xo/bb9+BH39c5rCNUGjATTftYLdJkoRCoQjoQBxPxBkAcKlimT1EIhGb11+tViMuLg4EQUCr1bJlLA8fPgzA7IY5bty1anRVVVWckb+9Cl5xcXGckpP9qf7E7NmzMXv2bM6+/Px8lJSUsM8BsE0aaD1SNRgMOHHiBKqrq2E0GtlZhS/jFxjCwsLwwgsvBHxGX8tU1rwyCCFoGujuFuP55yfh/PkuFBfbxgzce+8XyMi4DLFYjAcf/B2bibOwsLBfl1jsCes0yc5cTJkFW+s2Bw4cYL1x3OXChQu4cOGCS22rqqo4ZRzdRSaTISsrC1qtFrGxsX6pUldaWopjx47Z7GdG6OfPn/favZOSknDXXXfZzUhLURRbscza/dQZwepaCpgHJ/yaQYhx992AOVed4zJ35eVJeOaZTDbOIFjw1MzAlyiVSmRnZ6OhocEmZz7TgR09epRdi2AKztM0jfr6ehgMBnbdwdeo1WpkZ2fDYDCw9n+mdnNLSwuio6MRHR0NlUrFqTlsLwDNXa5cuYINGzb0+Tq9ZeLEiZg/f77f7u8uSqWyX80kvQWvDCz44x+B6mpz7EBYGFBeXo7TpxM5be69l8SECbF+ktB7uFJAxWQyoaKiAqWlpRyTmLXbp/WMoaysjDXRiMVit2omO6uzzLy6u7uhUqk4+yoqKmxiFJh4A61Wi9LSUgDgtElNTbVJbe4NdDod3nrrLc6+cePGISYmBgKBAGq1GpMmTUJ1dTUGDRrEGaGnpqb2WNO4vLwcCQkJDoPLnP0tLS1FamoqEhMT3VpDCPaZQUtLC4xGo79F8SohrQxomsYPP/yAs2fPAgDS0tIwa1Yx9u2bjl27ZgBItDnn0UdjMWUKMHy4b2X1Nrt27cKgQYPYTI3WvvyWHa3RaITBYMDhw4e9WvzD1cL0HR0daGlp4eyLioqyOf/SpUtoa2tDa2urzTVEIpFPFAFgTje9evVq/O1vf2P3EQSBrq4uTrvOzk60tLSApmmXO2Zr99rIyEhMnz7doeKwvD8AdHV1obW1FQUFBU7vaX2svLzcbvv4+HhERka6JHt/hRloMPWog5WQVgYEQbCKAACKi4sBALW1js1EAPDqq59h06YHvCqbr2EWZ+Pj4+1G91q+Ojs7odFoQBAE0tLSHHYszN/29nZERESgqqrKbkfsCFfjBuRyOetV5MzGTdM0pk2bhnHjxvl9JCuVSu3Oxmiaxs6dO+2uJ7hKSkoKlEol+xzKy8sdzqqYezK0trayCsgR9o61traioaGB04YgCIjF4oBXBszCMePIEKyEtDIAgGeffRZ6vR6XL19GWloaaJrGQw91obHxNPLzI3D//UNszvnsswfw2WfA998XYOrUQUHhZeCKmYjBE+ko7M08AODixYs4c+aMTd4eS3u/NR0dHejo6HBJlr179yI3NxeNjY3svlGjRmH+/PkQi8Ved7FkSn4ePnwYBoMBWVlZiImJcZgzacqUKRzXT0sT0dixY+06LfRF0fkjBUZ/JyoqCkKhkFcGwQ5TtMKybmxRURH27NlztYNx3Enu2HEIZ85UYejQoRgxYoQvxA0a7M0+ALBeSnfffTcAc0bL8vJyTluTycSJX2hra2MT51liOfKtqqpCW1sburq6bMwxZ86c4awdkCTJ3t+T0DSNv/71r5x9+fn5rGlOrVazayB33XUXmpqabOJWBAIBJBKJx2XjcQxBEFCr1Whvt5+vLFgIeWVgj+LiYnakKRbroNdzf3wvvPA6FIpr9sMLFy6goaEhaEdGfaGsrMylcqL+JCYmBomJiZg0aRIGDnRuIuwLRqMRK1aswMmTJzkmMMtZCsMnn3zSp3up1Wo8/fTTfCCZh4iOjg5691JeGdhh/vz5UCgU2L17N+6++xt8/vn97LGhQ4uRkTETYjEwceI4CIWAUAiUlFxCVxfAlzPgkpKS0qsayNbtT506hfb2drZzYzKUAuY6FVFRUWw1KsYN09Ffg8GAvXv3steur69HfX29jfdRRESEw7QR7mLtupueno7hw4dzIoCrq6uh1+thNBoRFRXFma0mJye7HNQYzCYbf6FWq6HRaPwthlfhlYEDpkyZgjFjxqClpQOff35t/4ULaXjssTQ7ZzDpFHwjX39n586dOHr0qE/vac8Hf/z48WxxHJFIxJbklEgk0Ov1dtcgHnjgAajVasjlco/Jtnr1aly6dAlbt24FAFy+fNktT6ybb745oCPcAx21Wo3Ozk4YDAY2xXiwEZyfykPIZDI0NdkO9adM6cDAgSbQtAAmkwAikQRtbR147rnA9prwJAsWLMCCBQsAeGbBGQB+++03tyOcjx8/7lb7lStXIjw83K1zXEEqlSIrKwtZWVkoKSnhlF11he3bt9tEXc+cORM33HCDJ8XkcQDjJNLY2OhVU6I/6VEZkCQZBmAYABpAAUVRrmUeC1BaWlrw9ttvAzCP8l9+2XYBOSdHjlGjNABMAEyIjx8AiYTG5s3A5s3XZgcrVgBeSo3vM8rKytDU1MSOoE0mE6qqqtjFNMbsYjQaERcXB6FQaOO2WFlZaeMN1Nraip07d/rnQzmhsLDQa/UMGIqKiuzuT09PR2LitdgWSzMRAOTk5HCSA+7btw/79u2zuc6oUaNQVVWF5uZmjhnKunyodZK6U6dOQaFQsKU3CYKATCYLieSLPRHSyoAkSSGAlVfbnAZAAFhMkmQbRVFvOTov0LHnlWKPM2eSLd7bb3PoENDLvG0+J1DSUdx8882s3zrTYWk0GqSmptrN2e8oCrekpATp6ekQCAQQi8UQiUTsy5vQNI3p06dj6tSp7BqG0WjE3//+d7dNR45wVB3OFZqbm1nvrZEjRyIiIgKZmZkhvxCtUCggEok4sRTBhrNv/l0AtlIUdcli308kSWaRJLmEoqgfvCybXxAIBFi7di3rBvj88xvwxhsvcNr8+qsOCQnE1dGWeVRVVHQZgwenAzCX0AQAiySZ/R5vxRm0trbi4sWLbIpqa+wFQlnGfDhLW8G8dDodNm3a5LL8PbFixQqXBwXuUlpa6tBEtGTJEs6ovbKyEklJSTYxBj39LSkpwZAhQzijfsvOXCqVOnRP5Ref7UMQBBQKRWgqA4qivmTekyT5NoARFsdmeVkuv3Lw4EF2dKRQ2AYz/ec/r0MmM3u00DSN5557DgaDDvxvyJbIyEjI5XIolUqXz6mvr3drKt7R0YHnnnuO0+ExiqO0tNTtgjtvvfUWVq1aBamHC1xXVVU5VATLly+3KcAjEAgwdOhQt+/T0NAQFIGQ/Q2FQmHXDThYcHVOTAB4AkDQFgLV6/VYt26dTfH0o0evt2m7fv2fAAD79uVDqTRBo9G4XHS9v7J582aYTCa2pkB/tOcD5tQTlrUMALNtvba21m57V7N8pqam4oEHHvCqOSQuLg5r167FF198wTEHpaamgqIotvC6wWBATk6O1+Swx8yZM5GWloaGhgZUVlaivr7eJjEbk14iPT2dVbbM83Kn8lqgEhERYRMAGUy4qgzGA1hlsf2IF2TxK2KxGC+99BKOHqXx2mvAtm09l7ibOdNcOnD16r9h7tzATml96ZLZGsjkZ/IGM2bMQHZ2NgDbCGRrWdLT02320zQNuVwOmVUwhzPTxowZM+zu96c55L777uuxzeTJk0FRFKucrWFMZG+++aZHZLJcjN6zZ0+vrqFQKDjb0dHRbIR/MKBQKNDW1gadTheUUeCuKgMdgFoAXT01DBS+++47uwVEPv74IWg07vlzDxs2EoA5r45KpeJ4gPBcY//+/Rg6dKhNCmaAu+ArFAohkUhsEt4BZu+lzs5OCASCgC4mZI3RaGTrLlhniLUMmrN+Pf7445yKbjRNo6SkBBKJxOezu88tA3JchCAI3HTTTV6QxvMw60iNjY2IjQ2+NPbOvIniKYqqvLppubpIkyQZR1FUlXdF8y6M1wjTATFmngkTjrutDDQaLbq6NKiurkZmZmZAK4MRI0ZAoVCwWUw9zXvvveeV61rzxBNPBJQLYHFxMbZt2waaplllaTQasX//fp8lSMvOzkZraysUCgUkEglbDN4VrN1gnWE9AAgUmJlPQ0NDaCkDACRJkk8BeJWiqANXd6gA/BHAFgABrQwWL16MxYsXs6Opbdu2obm5GU1Nrv8AGC5c6IDRWI3q6mrcdtttXpDWd+Tl5dndL5PJOCNQ6+ye9hLPjRgxAtOmTUNRURFIkrTbCdh7z5hwXO0ogsEDJiMjAy+++CJnn7ufi/n/FBYWIj8/3+3qXO6W2rSMx+js7ERra6vdiG7rID7rNoEyw5NKpZBKpUG7iOzMm2gfSZLnATx+NfAMANoBvEJRVMBnbGI6GqZgd3Z2NoYNG4bz5wuwcWMXtFrXkwxt23YrXnzRXEbw1VdftTl+6623QiwWe0Bq/2Gd6RMAsrKyWHOOZZUtprbAiRMnUFhYCL1ej/3799utrmXZ4c+ePRsqlQqtra2oq6uzayZiYDoUazt1IFJTUwOdTsfOVpnP1tjYiIqKCs4+AKioqMDJkyc5rrYGg8GmwI+3OXnyZK/PTUtLw/3339/n6/gSgiAQHR0desoAACiKqgPwmo9k8Ssmkwnff/89PvpICK3W9TJmTz75Dm67bQxaW0dg0KBBSE5O5kSRMgRjDdX8/Hynxy07KFeqRH3zzTcAzO6oTAfhrMgKQRCYPn06q9ADlby8PGg0Gk50MUEQ0Gq12L17tx8l8x7FxcV4+eWXIZFIMG/ePH+L4zJqtTpoYw1COjcRRVH4+uuvOftiYtxbL4iObsChQ78hMzMTbW1t0GqD1vu2TyQmJnJMS8xfmqahUqnYEpUCgQClpaUYOnSoyx4bga5o586da3d/YWEhMjMzObUbrKFpmp21HTp0CE1NTejs7IRcLmfrPPdnnH22/oharfaqx50/CUll8NFHHzn0F46JqXPrWuvWrUF8fB1EImD+/GZs3Gg/0jbUYZ63VCq1Ge1XVVWhoKCA3TaZTNi1axd+//vfB22GSFcRCoUcV1om4ppZH2hubsYHH3xgc15dnXvfY39h7Sbc34mOjkZnZye0Wm3ArHW4Skj+0v7rv/4LJpOJ9SCiaRoXLlzAzp07YTJ1YvLkHBw+PAWJiVcQEdEGo1GIixcdJ+uqrDR7UXzwwQAMGrQOIpEJt99+O7KysnzyeTzBhAncOAl75hmTyYSTJ09CIpGwphl7awnOYNwnGQQCAaKjo2EwGDjFQ9rb27F+/XqH15k9ezamTp3q1r0DhY0bN7LrLr5i5Eize3RLSwub+8nyO5Camsp+n+05C1y4cIGNlnZUxc7R9unTpz3wCXwDE9nd0NCAhIQEP0vjWTyuDEiSnIBr6wxrAKwGUArgYwAaAJsAtAJYDiDMcpuiKNerpfcRgUCA/fv34+DBg1b7gXnzdiMtrQRffNFzcBDDqlV/g1TaDYUiHL///e9BEATHBtzfI5TdKcDuyam9yWRCXV0dFi1axK61FBcXY/DgwZzOiHkvEAiC0q3PksceewyArTdRXl6e26k17JGZmcn6zNM0DY1Gg7Nnzzo959y5czYptAHzTI9R8JazPstYCJPJBJIkcdddd/VZdn/DKwP3yACwAMA8AKMArKQoKh8ASJJ8HeZMqFIADwBIsdr+Py/I45BZs2Zh1ixzmiXmh2cymfD00z/giy9udfk6cnk7ZLJuCAQC6PV6bNiwgfNjIAgC2dnZATVT8DU//fSTw2Pz5s3DpEmTfChN/8EypbqnYNJeeAJGEYwdOxZTpkzhBBRap8rW6/VsUGEgxRdYEhYWBrlcHpSLyIQzbw0GkiSHAJgIQAAAFEX1WJmDJMm/wpzGohKAGsDTMMcs3EeSpATA/wIYSFHU/cw2RVHPO7reyZMnaWdFR3qy4fV0XKPR4MiRIwCAv/zF9QyeKlUjfv/7f9rsj4+PR0xMDLutVCqdjmj7Kn9nZyfGjh3bq1/YyZMnaWedsSvExsZyUktbYjQakZyczB5nsPzuiUQi1vTkrj22N/Zbd8/p7fPt6XvrDFdkrKqqsnF17OrqgkDgPJ0Kk36kP7Bo0aI+fXd7er6uPEd32uzduxcymczuAMXT9+prG3e+t67ODD4E8BVcTFRHkmQGgEQAd1AU1X7VdPQQrqWzEAJQWFxPCKDHGoPOAnB6CtBxJYBn+PDh+Pe//92TGByamtRoaYlAVBQ3SrSyshKVlZWcfX/84x8dxhv0VX5/+2o3NzejubmZ3bbs6E0mE65cucJu2xsVLlmyhLU5e6oymifP6cvz7W1AnCsyWh+naRpvv/02mzjOerBH07TNuk2g09MzcuU5utOmqKgIVVVVdtt7+l59bePO99ZVZZAD4AOKono0fJMkKYN5zeBhAPcAeB9ANAAjgPMkSY6F2SykAdBusX3F/hW9w86dO1FYWMgpkm4ymaDX6zFu3AmcODGu54tc5a23nsNLL/0PRCITRowYgdTUVIwePZrTprCwMOADz5zBuNTGxsayNm+GYIgQDhQIgsCCBQvsKgnAHPvxz3/+02cpLoKR6OhonD9/njUBBwtOlQFJkvtgLnc5EMD9JEleAnqsZ/AHAANgVggHSZL8HoABwGMA9ADegXlW8CiAbqttn3HDDTdg0qRJrH2zqKgImZmZV5N/AePGmVBe3nPmUjME6uvjMXRoOyZOnIi4uDivyt6fqa6u9kjVtBtvvJE1G1mnrADMMSJFRUVsmuWFCxcCMKc+GDZsWI9mkv6A0WjEqVOnoNVqUVBQ4HMPor7AZJ+VyWQQCARobGyEWq122jlaHsvOzmbrN/h7Vusu0dHR0Ov1aGtrszF9BjI9RSDPdPeCFEWtgdmLiME6leH9PWz7BGs7o0wmY1MbREQAjFVj7969rMfR668/j44ObvqDRx/diNRUE0aMSMagQYMC1ve4t5XOLL1HmL80TePAgQPo7u4GTdNobGxEZGQkOwszmUzQ6XQOawEDwC+//OKW/Dt27ABg/qEOGTLE44VpvEFZWRkrN2BO0RAeHg6xWAyBQNDrTjIjIwMxMTEIDw/HlClT3Bq99nYWF0qzP2YtsKGhIXSUAQNJkndRFPXN1ffvUhT1pHfF6j8wHkeFhYVYtUqBW24x4ddfr40633//mkkkO/sctmypdivbY6BhNBqxbds2bNmyxa9yJCQksF4qcXFxmDlzJseTJRCm8CkpKRg0aBBqamoAOK4lkZKSgtTUVABAU1OTw2SCDEVFRayitaxNMHjwYDzwwAMekDy0UalUIAgC9fX1GDx4sL/F8Rg9KgOSJKcDuJkkyaqr7dO8LlU/5cEHW/Drr46LdQwbdgHl5cqgHiEJhULceuutHlvgPXLkCNra2tjssSaTCY2NjVAqlZx9TOF45j1zzGg0or6+HsePH7fJ+W9JRUWFXQ+amJgYzjVNJhPa29sBAM8884zHFPsnn3wCjUZj99ioUaMgFAohFAohl8sxatQo0DSN0tJSNhjMkltvde72HEqjdH8gFAqhVCqDzr3UlZnBDJjjAWbCvAj8kjcF6s8UF9sWSf/qq59x9903Xf0B3uEHqQKbiRMn2uxzpzNjUjYPGzaMVQDV1dV2UzTYo76+nn0fHx+Pu+++G4D5B99bl1B7LFu2DN3d3XZjBm644Qa2Ipjl2kgwVtMKFmJiYkJPGVAU9TJJkhRFUV/31DbYOXSIhvWSwLJlN+GJJ9rwu9+p8Prr/pHL1+zbt8/vZiJPMHy4OTvtlStXkJGRgRkzZrCRuZ4mLCwMYWFhWLx4MX788UfOsX/84x9euacjZs6ciRtuuMGn9ww21Gq1R4P3+gOuupaeJElyAwATgPcoigrOtH1WlJSUYOfOndi7dwD27p0JgmgCEGPTrrk5Ahs2RAS0MvCEB5C3GD9+PKKjozl1EyxfVVVVrGcKcG0xu7a2FlVV3BpMAoEAiYmJ6O7uxo033ujTzwEAo0ePxujRo1FZWYmWlhZotVoYDAbU19fj2LFjbPZWwCy/NzKyWtY7doWJEyeyaab7+zqMr4iJicGxY8dgMBiCJpmiq59iHYBXYY4H+AeARV6TqB/BFLPYseNGdHY6j4l74ol38fLL5kyREokEf/jDH3whYkBBEARWrlwJsVjssNIZ4L7NOywsrFdBZ/4kPj6eo8AA2CinwsJC3HnnnQC4uX7seXBZtsnLy3PbG8sZR44cYaPze2LkyJG45ZZbPHbv/kp0dDRomkZTU1NAl7m1xFVlUEVR1GmSJMUA7Od+DkI0Gg0KCwvx4ouF2L//Buzf79jT9tNPH8GXX+6BSGT+QW7fvh3p6emsKSKYEIvFGDdunMORujWW9XGPHj1qc9zyXIFAgFOnTuG7776z6fxSUlIwZ84czjmA2czDRNWqVCokJycHxAjWaDQiLy8PV65cgVwuR3x8PGeRXKPR2JjjMjIybJSINXV1dRCJRJxEid6GJElMmDABaWmh4V9imbAu1JTBTJIkT8KcNsJIkuTeHgLPgoIbbriBta3eeWchpFLg22+B1att20okWpw6dRJCIY2bbroJANhFwUCgt3EGfW1vPeKlaRqtra3s7MHyJZfL2ZTZlqPh48ePsx1feHg4xo8fjwMHDrgs37333oshQ3xTh6KjowMbNmyw2R8TE4Ps7GzWXVYgEMBoNGLWrFnsNkEQbAyCMwoLC3HHHb1zZuhLnEGoKAIAiIiIgFgs5jggBDquKoPxMLuUdgKopihK7z2R+i9pacD337cAuNbJm9NQGNltkUiEkSNHBnXqCU9ibzaRmprqVod022232bQfPHgwtmzZAqPR2GMuns2bNwMAhg4dioEDB7JZNWmaxqhRo/ocWFRZWYlLly7Z2OmTk5Px0EMPOZzFhIeHB0WN52CEMSEHk0eRq8rg06ttPwGwHuZ00yHHO++8gyNHnuLs6+qSITVVitGjRyMlJQVGo5FXBF6CpmmcOHGCU1CHIAjU1dWhvr7eZu1h8uTJnNgDo9HI5qBiokd1Oh0uX77MFji6cOEC55779gOG8DkAACAASURBVO3DmjVr0BcYN9fhw4ezPurh4eEgCAJHjx7lzHIs5VWr1X26L493CVVl0AyghaKonSRJ3uNNgfozTz31FJYtoxEdfW0k98Ybz2PlytfR0HCtcPn8+fPt+s/z9B2j0chm5ATMHaher0d3dzdaWlpw/vz5Xl87KysLS5cu9fh6w9q1a3Hq1Cm7xWEA4PrrrwcAm1oArqSX5/Ef0dHRKCkp8bcYHsNVZdAA4EaSJAcCCIziqh4gJwc4dAgQCoG6OjXM+ecIXH89jaNHr3UYFRVLcN99lzFjxgw2kjQQ0iH0B+x5xDAjeHttxowZw9nev38/mpubodPp3C5U/vjjj2PgwIE++T+NGTMGY8aMwRtvvMFGODMwi+rjx4/HvHnzWHkuXLjAVsizlJH/XvUPmHrIXV1dAVfL2R4uKQOKotaQJPklAJqiKMrLMnmNhoYG/PLLLxAIBE6TpAEATQMvv2y5qDrI4j33x/jll5kYOHAHWzpSIBDAZDJh+fLlPXp+9Bd8FWcwefJku2Y0SxNPXV0dWw/ZXido+be+vh46nQ5isRgpKSkQCoUQiUTsy/ocoVAIlUoFtVqNQYMs/6e+4fnnn0dOTg7q6uqg1WrR3d2N7u5uGAwGaDQa5ObmssqOMX8xOJspWCrI+vp6VFZW2ricWrZz1VXUFyxaFJie6pYJ65iSrYFMTymsRwL4GwDm19tBkuQqiqIuODmt3xIREYGJEyeivLwcly9f5nQUfR3Jv/32s/jLX17GqFGjsGTJEk+JHHQcPnzYI9d57LHHEB0djXHjxuHChQu47rrr3Drfn3EGU6ZMAWDumJkZELNWUFlZic8/t070y9MfsXQvDXplAOCvAB6kKKoWAEiSHADgXQABmYRHIpFgyJAhiI2NxcCBAzkFuysrK5GQkMAZfZ07V49vvrGNOLbH4sXmFANnzpzBmTNnMGzYMCQnJ7P3qKmpQUNDA2eUNmnSpH6z2PzAAw9AJBI5HE1abpeVlSEpKcnuMXvvKyoqEB8fz0k6xxy397eqqgqxsbEOg6p27dqFjRs3cuS3rlAXFxeHtrY2G5PMtGnT2LrX/uIf//gHO/OxZsqUKViwYAEIgkBNTQ3i4uLsBugx6wuW+5lXRUUFkpKSHMaAOIsPKS0tZTOkMm0Bx/8ry+8E833vqb3le4IgbP5HgYJUKoVCoQga99KelIERQBhJkskW+/p/1ZAeUCgUNsFgYrHYxj3x66/NL8DW/3r9+iq89NK1IjY//rgYLS1RmDlzP9u+oaEBSqUSBEGgo6MDCoWC8+PrTwuETCpeg8HQY7BSe3s7qwwA8496x44dfVq89RSJiYmQSCQQCARQKBTo7OyETCZjlb9Op/Prj5emaURHRztUBjk5OX2+x7Rp05CSkmJjIrOUgUEoFHIGJNb/W1fp6urqdTrnQCtuY0kweRT1pAxaAFgbk1u9JItfaW9vx6ZNmwBY2lVTsXHjjKstrP3ebauZHTgwHenpl5CaWonIyEjo9XrU1ZnX2wUCAWbPnt1vK3D159xEPSGRSDimF+u/TNAWEz/gz2ygjFxTp06FTqdj15k8ycGDB9mCTK4QGRnJztYMBgN+/vlnm+do+d7Roj9zDXvIZDJERUXZzFBEIpHbJr7+RHR0NMrLgyMpQ0+Vzh72lSD+5tixYzYjxsJC90c6H3/8X7j11m149FEZ5s6dy5ZuLCws7LeKoLdkZGRg9uzZAMzBdowN1RpvF7jvTdRsRUWFW+09hUAgwL333stuO0uW5+rnomka//znPx3ONnqitdX747uuri5OfIglgawMYmJikJeXFxTeg8GRbs8DTJ8+HSRJArgW/PPQQzT++78NOHFCgM5OA1paXBtRnjgxFWvWVARd52+NZUUtwJyR054duqmpCaWlpU6vFRkZyY7ueRzD1G+4fPkyTp065W9xbHjyySfZgkD2khBavwcC30xkMBjQ0tICpVLpb3H6BK8MrmJtOwWAlBSAST1//PhZCIVCUJQU99yT5fA6d965BU89NRSjRo3yprgeRyqVQiwWQywWQygUcoKemEW/nmyjp0+fdut+lp1CeHg4G0xWV1fHloLsiXPnzqGxsREJCQkYOnQohEIhdDodmpubAYA1DwkEAsyaNcuvJqKWlhYcOnQIJ06cAGD+zhmN11KZZGdnc2ZXls+htbXVredrjVKptJsHyvplNBrZ9azo6GgMGDDApvMmCAJhYWHsQjcQutXVLD2KeGUQIigUCgwbNgxjxgA5OZfwzjvcxGY7dlzGtGnpUCju9JOEfWO1vex7DuhPZp8ZM2agsLAQJEnCaDRCo9Hgiy++sNvWXsZUeyxduhTZ2dkuy+sq+fn5rCIAgDlz5mD8+PEczyBLrJ/D4sWLXbpPXzrmUO3Ue4tKpYJAIEB9fT3S09P9LU6f4JWBCxQVFWHLli0wmYDXXnsRWq1thsuFC9OxcSPw6KN+ENDDMDOCn376yaOmiKysLDZQxxnuzAys2+v1emRlmWduBEHAYDBw8g1NmjQJ8+bN80unJ5fLERERgba2NiQmJkKr1XIWeq3tzo6ew5kzZ9DS0uJ1edPS0nDXXXcBsB/wx/wNdnOoMwQCAVQqVVB4FPHKoAe+//57nD17FgBw5UoStFrHYeePPQZkZOzHxInXB1x4uj1vIuaH7sxLxB3y8/M9ch3rBcf29nZ2oR4Ap3MSi8X9JlfUyJEjMWLECPz6669oa2tDXV0dCgoKfCoDU9fZkZeQyWTi1InYsGGDjfeQpWnLmlWrViHMujZskBMs9ZB5ZdADO3fego0bXa/c9MMPZ/Cf/xyATCbDihUrAqYkXnx8PLq6ujgeKZYBQv2Jc+fOOTw2btw4myRvRqMRV65cAU3TEAqFSEhI8IWYdiEIAhkZGdBqtaBp2ufKoLOz0+1zIiIiOGst9kxaOp0OUqkUH3/8sdN21vsEAgHGjh3rtkz9iejoaI8NdPxJYPRUfmTuXMAq2NUhQ4YUQaVqQWRkJJ555hnOaLW/s3DhQnYmoNFocObMGYcmgK6uLoSFhYEgiH7nY21pk7dHWFgYXnjhBR9Jw4XJXJqQkMDmUMrIyIBIJIJAIABFUR6vTiaRSBATE2N3Edje+87OTnb2wDBr1iwoFAqHAwOaplFcXMwpbmOvrfU+ZruystLFT9M/iY6ORktLCyebbiDCK4MeWLrUnLTuzJkCvPnmOXz++V0O2166lAHA7PnxxhtvsNNunU6HrKysfr0w9+GHH/pbhF6TmZmJrKws6HQ61qTFpL3IycmBXq/HgAEDoFKp2I7YH2RkZOC+++6DVCqFwWBgs7NWV1ejvr4eJEmyOYpMJhNKSkr63MHodDq3O1trkwcTjMnw4osvQiqVsttMwkB7nkeuEAzKALB9boEGrwwc8Omnn6KsrAwAUFERjw8+WA6g53rG//73rfjmGwkb7RoWFoaEhASboinBgL1cN76su8tw8eJFdo1GpVJh+vTp7LHJkyf7XB5HREREICIigt0+efIkSktLfZLGQyAQQCKROIwsNplMTtcCLHnttdfcuvc999yDjIwMt84JJBiniPr6+oCyBljDKwMHPPTQQwCAPXuAqzXYXYKihiEn5zN2hKpWq7F06VIvSOh/7HUsvmbBggVoaGhgE+EFUpnI2NhYEAThE2VgMpmg1Wq9fh97fPnll06PB2oKawa5XA6ZTIb6+nq/pEX3FLwy6IGZM4H/+z/gzJl6GAxGNuI2J2eq3fbPP/8ayssNiIiIwIoVK3wpakAQGRkJmqbR3d0NnU7H2U8QBGbMmAGpVOr1dBT9gYSEBCQkJGDMmDF2j7v6uXQ6Hd555x3o9Xp0dXVBIpGwAZQdHR0ekVUkEuGWW25h3XZ7K6s9AjkCmSEmJoZXBsGOQAA89RSQn1+DrVu3gsnuW1MziF0jsOSVV/6EhQt345VXBvhW0AChvb2dnUlY2peZ/Dg//PCDz2QhCAKrVq3i2L+9SUFBAb799ttendtTWgfL9zqdjqNoewvjQUQQBB544AHExsb2+ZrBSkxMTMCvffDKwEUEAgHWrl0LmgZMJmD1aqC7G1i40Fwe05IdO+bgT386jMbGRgD910XTEnfiCWJjY5GWlmbzmXJzc3s811MxC5YkJCTgd7/7nVvn+GM2MXz4cKxduxZVVVV4//333Tq3J5PciBEjIJfLodFoXErCRxAE1Go1x1ssMTERGRkZATnL8jcxMTE4d+6cV77fvoJXBi5gHnS5/gOJi6vEnj2HcPz4cXbEZrmo2R+x/hJHR0c7LIKi1WpRVlZmk5DOE4hEIk7xG1eUaEVFhU3Q3MCBA/H444/3u0yS27ZtQ15ensPjTIyEu4OH/Px8pKenIzw8HHFxcZDL5QB6Li5j+be+vv7/2zvv+KiqtPF/Z5KZ9AKZhJbQhBxAKSFR3gCCiNIWC7A2dl2778+yiFIWK+vqrr91XwtredFd+4qr/ER0RUARexQBBUTghiKEEtJIQiB1yu+Pm3uZcqelzUxyv59PPpl7zzn3nntm7nnOec55noeioiIKCgoCClLjfL2GhgY2btzoM4/iLwqgR48eqmvxnJycoJ41HElPT8dqtbbIjiNc0IVBAPz3fwduawBQXNybmhqYO3em6q8klGEWA+H2229XR4neopwpf7/88gv9+vXzOO+c32az+V04BNma2NmIrDW7kXJzc8N+MXLWrFnk5uayZcsWTeO5QEaWv/71rzXdQiifjxw5okY685bH2/9Dhw6pQWr85XW+3i+//MLAgQN95lFiHrj/bo4fP+73mcMdZUdRTU1NiGvScnRhEADLl8Ndd2mrFWTLS1dPmPff/yiPPvpAR1WvTUhPD3yNo7q6OqCYr0uXLgX8q2Rmz56tfv7444/59ttvXWYdgW553Lp1KydPnmTu3LkB5Q8V+/fv92lFDbJHU+cO87LLLgvYE67dbic7O7tFdautrW1RpLMTJ06QkZHRont2BmGQkpJCdHR0h8SGaC90YdBKDh06xNVXb+Hf/75GPffllxM4erSMPn30ReRgufjii1V31CDPNJ5//vmAA7fs3buXL774IqzVcl9++aXfPDk5OfTr149hw4a1uSpOp+0xGo2kpaXpwqArM3jwYBwO111FX301gUWL3kKIQm6//fagRt1dHcUFhrPxzvXXX89TTz0V8DU+//xzPv/8c5dziYmJ3H333SH3sKnlEHDs2LH07duX7Oxsl06/M0bH68xYLJagvO2GG7owaANGjz7N22/LC3ajRv3I+PFfc/XVoxk37ho/JXUCITk5uXknlzxL8BfQvlu3bsybN6+DahccS5cupaamhieffBKz2UxjYyMFBQUUFBTQvXt38vPzycvLC3U1dVqAxWJh3759oa5Gi9GFQQuxWq3s2rWLMWPOob4+QT2/bVsOy5bV6CEckVU8K1eupLy8nE2bNvndIaOk19XV8d133wGoUbeUtNjYWDIzM306yKusrHQZgU+ZMoX8/PzWPk6boSwyutsCnDhxgo8//pjCwkJAtsn48ccfPVRESltYLBamTJnSATXWCQSLxUJDQ4Oms79IQBcGAdLYCDU1stM6mw1OnWri2Wc3UV8/wiPv++9vZtSoUSQnJ4egpi1DS30RKFo+ipTjjvbk6ByhTOk0GxoaOrQO/ujduzf33nsva9euZdu2bS5pTU1NLnGlfbF3717VetfXzq6WkJCQwD333KOrqYLA2UdR3759Q1yb4NGFQQB8+OE6LrlkmtvZOOAWzfxPPrmApqa/k5ZWSV5eHr/61a/avY7twahRo1w6d6VjqaqqIiUlRT2nbBdMSEggOjrao0OqqKhQg6R767Ccj6uqqlRB6m/XjTve/Px88cUX6ufZs2eHPM6E0WjUjGVw9tlnq/vvq6ur6datm4eANRgM1NXV0djYqGkHcvLkSVJTU73agfiyGaiull2w19fX88gjj7RnE6iE+3bgQFG8l5aVlenCoLMyeLBnmEt/PPPMPJYufZgffvghYoWB+6jVH7/61a8YNGiQR+e1d+9edXHU255352NJktStqJdffrmmKwZnWmJNHGq7j+joaO69916fedyfSxGadrud9evXc+DAAa8O6A4fPtzmdW4rzGYz559/fkRY5geDyWQiISHB75pWuNKlhYHzy2W1WlV/+M6j3a+++gpJkvjd7/ry+uuz/V/UicOHsxg1qm0chbU39913HwClpaUtjm2wZs2atqxSUCxevDjiQo26s3nzZj766KNQV4OkpCRuvvlmnwLb/bMkSQwZMkQzn/N/LTqDozqF5ORkXRhEEuvXr1cXKN1RfsjulqCbNw8J+j4vv3wjt9zyDx5++GEyMjLC2ufLX/7yl1BXISCys7Pp3r07AHFxcZx//vns2bMn4gRBWVkZzz//vIeKx2QyqZ9DtdZRU1MT1Fbe1tJZ1EQgC9LS0tJQV6NFdElhMHXqVKZOnQrIqpCPPvrIJRC4lkuEq66SvU2uXTuNTZvGBHQfs7kBi6UM8AziHm4o1q3KtP3QoUMuvmRaghACSZJaXH7o0KFceeWVrapDuJKens5DDz2k6e5DmZkWFha6WBI7q1SeeeaZFscnyMrKUl2Gg7xoLUkS+fn5VFVVufgNUuw+nD9rrVMcO3aMPn36aM4kDAYDWVlZqr8kdzrbzKCwsJCmpibVjXik0KWEQXl5OZWVlQH5zPHG9OnrAhIGI0du4/LL38dgkF+McN9ZFOz6gDPeVAkHDhxoVZ12796tuctJ2W6qkJycrKooIgnntvrxxx/Zvn272vGCvLV0x44dVFVVqR5w24KSkhJKSkpchE9UVJS6/bct11Puuecel+hunR3lWSsqKiLO5XeXEgY//fRTQK4AtNi1ayjvvBP4KLWpyYTSN/Xo0SPsRwlXXnkljY2NpKamYjQaKS8v54MPPgiobHtHPHMenRqNRubNm0dcXJzake7evTviBIE7AwcOJDU1VZ0VgLwI3LdvX3bs2NGmwiDQWAcmk4lLLrkE8BzpK3+nTp2irq6OsrIyevXqpYbXdHZgp6C1jtDZtq4qg77y8nJdGIQz48aNUwOlAxQUFPgtc91119GtWzduvjlw/e2KFU3MmZON2bxUPRfq3Sv+eOeddzr8nqNHj2b8eDlinPMoed++fc1uPrTj9TY0NFBfX09ycnLIt4i2FSkpKep2XZvNhsPhoK6ujqysLDIzM5k+fbpfW4KKigpef/31NqtTU1MTq1atCji/L9fcvuhMawYxMTHEx8dH5CJy53iTAsRsNjN58mRAHnVt2rSJ9PR0kpKSqKmpIT4+nhMnTqi68rS0NF577TUAzjorFvhDQPeZO9fE7bc/R0aG/IMwmUxceumlbf9AbciUKVOwWq1s3brVo/N1HnU7HA6sVivx8fF069bNr598u93udZvjDz/8wA8//NCqeoer24mW4nA4eOGFF6ioqNDU0bdFHON+/fp5Ve2dPn2axMRE9VxlZaWmPYNS1mQykZOTQ0lJCb169fJ5X2fB7fwbC9QrbaSghMCMNLqUMHCmpqaG4cOHq7OExsZGkpOTiY+Px26306dPH3XXCsijrk8+2aAef/31Nzz88FKP6yq88cYtVFebXFQZ4Ux+fj4NDQ0e/tid9fPKs5w4ccKlbU6ePMmePXs6rK7OHdJzzz2Hw+Fg4cKFEbejCGT9/fLly72mt1VHGRUVRWpqKmazGYvFomlMqJxTDAQBdbbiDbPZTG5urottSLB0pgVkkIVBINHmwo02FwZCiBRgGWAGlgP3AieRzXVjgde8HUuS1GH+X4cNG8awYcPUY3+GS+7pEydOZPHiWsrLy7n33q9YseI3LvlraswYjbBz537i442UlZWRkJBAWlqa110VoWTdunXAGT/64H1vuNFopLKyUo3q5EsH3aNHDxoaGtQ4w6Wlpa1eW3A3VjKZTDz77LMuO3FMJhPz588PWzXSihUrVLcTsbGxbTLid0cIoSkgfQkYq9WquhBxbmOz2cy0adM6nY6/PbBYLGzfvh273R5R7dUeb8q5wINACrADGI7c6f8O6AcsAmK8HD/bDvVpF0wmEyaTib59+3L//b/B2wall19+n0GD4rDZbOzbt4/8/PywtDfwNwJUUHTZ27dv95vXbDbTq1cvqqqqOHjwYCtr6Mnvf/97TCYTe/fuRQjhosaIiooKW0EAMHfuXNXYcefOnfznP/9p9TVnzpyJwWCguLjYp8rG12K7yWRyKavkjY6OjvhF+o7CYrFgs9lUdyKRgqE9doAIIfoBfwbiJEmaI4QwA48BGZIkXevtWJKkBd6uuXXrVocvT4D19fXExsaGLL2urp7cXM9Yrnfe+SwWS4XLuSFDhjBihKuDu9bev7a2ltzc3Ba9rf7aNti6aOU/ceIEp06dwuFwsH//fk6dOuW3XO/evYmKisJut6shE1NSUigrK/NqA5GXl4fZfCbyXGpqKomJia1+hpa2r7e2tVqtlJSUsGnTJqxWK1lZWQG5kEhLS6N///4YDAZSUlJUfzjOBPtsbVG2Nfds799uIHVryzw2m401a9Ywfvx4evfuHdL6BNO27TV0qgB+Bm5sPo4CEoF6H8d+dSfBqHE6Ov2FF9Z6nPvTn/7M1VdfyeDBg9v9/q3Vu7pfe/Xq1QGN/rWIjY11WWS02WzqKF05ryxSOp8DXCJFHTt2zOPa/qw7t2zZon6Oj49n5syZmu0WrD+j1rSv1n0OHTrkslPHmyBQjMOcA88r6zPKubq6uhbXbeHChS5qy5b4eWpNOWj73647gdStLfMMGTKE9evXExcX1+LfXlvlCaZt22PNYA7wOfAW8EchxGhkNVERcEoIkYusFnI/Dl/PWs3s2weDXYOa0aMHpKdDXd0kl/NGow0wUVRUxGD3QmGIov9XOpjJkyczefJkza2M+/btIysri8bGRq/bP53LFBUVqTGTve0+Uj5/8803VFdXq5HOvFnoKh2jN0aOHElWVlZYquTKysrYs2cP0dHRHtbu8fHxqp7Z4XBw4YUXqm3hywIYUMsVFRWpXjPNZrNaxh2DwUBNTQ01NTVqenV1NWVlZWq6+39vPodqa2vVzQfBlOuMqieDwUBaWlrE7Shqj5nBHuB1oA7IAxYjzwJuBRqA53wchzWPPeZ5rqRE/pPl3Rns9igeeWQ+77wTesdjgfC3v/0t1FVoU7Zv38727dsZOXJk2K0dbN++3atvLEUoK6xevbrd6pGcnExcXJyLUG5oaHDZXqz1X+tcfX09H374YYvq0ZnsDBTS09N1YSBJ0s+As8/ma92y+DsOW156Ce65B+rqwGAAk0kOdtPY6GDjxq955ZXzkKQYNX9Tk5ni4uIQ1jhwHnzwQV599dWQuz42mUykpKSou2CUDsdsNjNu3DgPmwc4MxKuqKigurrapUw47uYYMmQIKSkpREVFqQFRQH6e4uJiNc5Be30X3bp1Izk5mcmTJ5OVleWSFolqonAkLS2N/fv3e9jphDPhNWSKAM4+2/OczWbno482ctNN37F48SKXtNtvv41Tp95l2LBYTp8+zfDhw9VtluGE0WjkxhtvZN++fTQ1NWEwGNi0aVO77ALyRVNTk+aIKjMzk7S0NFV/npSUpHb0DQ0NDBw4kIEDB3ZoXVtKZmamqjZzprGxkVdffbXN7zd16lRSU1MBWeAMHjw47GZLnY309HTq6uqora0Ny63kWui/iFby0092rruukhMnrsNu1x6Fnj6dSnn5Eex2O0OGDAlLYdCasJcdwZEjR3jppZdwOBwYjUauu+46F8O3zkCg4S6DZf369YBsW3PFFVe0yz10XHEOgakLg06Ms9Xo229fwe7dwwCL1/w227ccPGgjJiZGc5tjOLB0qac1tcPh4MCBA1itVurr69Xp7tGjR/n+++87uorU1NSQl5fHli1beOaZZwCYPn06BoNBXXQ2Go0u/50XL8MJh8Oheg7dtm1bh7Tnrl27+Oc//8nNN9/c7vfq6ihbfsvLy+nXr1+IaxMYujBoAT169FA7z8WL4ZZbGtmw4QQlJa5eCjdvtpOXZwQeaJVONVQYDAbOOussj/Mmk4np06cHfJ1gn91XfofDQWlpKVVVVZw8eZK1a+UtvcFECBs4cCB5eXkB528r3n//fU1X4bNnz2bw4MHExMRgMplITEwkMTExJPp7nbYhOjqabt26qTuzIgFdGLSSXbscvPmmGfB0V3vuubLa6Ior3mH06AEMGhR5AS+CJVxCN/oiVPryXr16aQqDYDyDtgVz5szhnHPO6dB7dkUsFgsVFRX+M4YJujBoBbL2wb8KYuXKK1m5EnbteoOzzvIM+GIwGLj88svbvoJB8Pjjj2vuYXe3HwB5kfejjz5Sj4cOHeqiF504caLLtcvKykhPTw+4Ll988UXQ9b/88svJyspy2UnkK6hKKBwHjh49mszMTGw2G++8805AVtjjx49X7QwCxd+i/+bNm9m8ebNmWlFRUVD38kZUVBR/+MMfOv3gxxcWi0XdGRYJ6MLAD9XV1Tz99NPq8bFjPfnXv35LbW1wi0L/9V/byc8/DWQAkJGRwZw5c9T0UHs1bakVa25uLmPHjvXpgyVYtUWPHj3aTK0UTkRHR6vuCRYs8Op5RcXXc/3000988803mp5Hjx8/3kY1bjk2m42dO3e6hNcEvH7ujFgsFqqrqyMmBKYuDLzQ2NiI3W4nNjaWJUuWsGfPHrKzs5k2zRa0IIiPP83YsRsoLT3FyJEjyc/Pj0h3y84onf/+/fvZt28f4N2q2Gq1qo7YvOVx/qy4FVfOOxMbG8vChQuDHi1HAp999lmLI/GFI4FGyuuMRmfguqPIX6yHcEAXBl7YsGEDO3bswG63YzAYVDfNBsNE4IKgrjVnziqSk2WVgGIZe//997dxjdueK6+8kvj4eA83CAcPHmTAgAHqyM7d5YG7+4G9e/eSnZ3tcyD//wAAIABJREFUcd7b5z179jB06FBNNwadmUmTJjFp0hm3Jt5mBlrOJR0OB01NTTgcDgoLCxFC8Omnn6pGj+5lItHffqShC4NOwowZM5gxY4Z6rLyY8+ZBsNvb33jjWuLjT5OQEM2yZXsZNswWESPbtgqFOXDgQAYMGKAeK7MAxRNpfHy8i34/KioqItqnpdjtds21GPdZUlNTkxrnwDnt5ZdfDsuFyYyMDPWzu/D2FRujsxIfHx9RITB1YRAkLbUfqa1NoLYW/vGPMiZO/NLF58wFF1zg8iKFmoyMDK+qGofDQWNjI2az2WO0qdXBARw/fpxXXnnFpUNQgttHRUUxYsQIxo0b1x6PEpasW7dO9a7qy6mbYmCnlScuLs7rzMpqtao6an8O4pSybdFh+fMo605cXFzAcTQilUjyUaQLgyAxmWDqVAfr1weutrBYyrjzzudZsmQJMTGTgEkeeUK9gOxMaWkpt912m1cB1ZZ2A10R91mnN1prZ1BUVMQrr7zSkiq2ObNmzVI/u68VOQ88OhtpaWkcOXIk1NUICF0YBInBAOvWyYKgogI2bwZ/9lfl5em89NJirrnGyKBB9rC0ih09erSLqmbXrl1et8WVlZVRIrtqDYiW5D9+/LjaUfhyY+1wOKisrKSwsBCHw+ERg2H27NkMHz484Hu3F4rB2ahRowIeDftrN29O0JzLuW/z1bqG88J9eXk53bt3Z/v27R4eVP2RmprK7bffHvTOmc7oqE4hPT09YkJg6sKgFaSlwbRpsufS3bt38+233/Hpp2ZOnUrigw8udcl7+HAcQsANN7xCv35FzJkzh169emlGqgoFl1xyScB522pmUFxcTF1dnaoyMhgMlJaWtsjOwBurVq3i8OHDAY3E2wvF5QSg/r/yyivbJHhJW5ZzLjtlypQOu2dnRgmBWVVVFfa+tHRh0IbceOMN3HijbJQ1btxmNm8+1yPPK6/cwD33PMm7774LyEFYBg0axJEjR7BYLEEZZ0Uq+/bto6Kigk8//VQNvh4MJpOJW2+9Vd3ltH//frKzs9WFZ2f/ROEwAzMYDMyfP9/FXkVrcX7+/PmdXofe1XDeUaQLgy7IgQMHmDp1naYwAGhqMmE2m4mKiqKqqoo9e/ZQU1PTqvCFkYLD4eDNN99s1TWampp47rnnXM6tWbOGxYsXtzgOb3vjHM4T5AX0ESNGYLfb1b9169aparGjR49SW1vLokWLCCY+tU54kZKSQnR0NOXl5WRnZ4e6Oj7RhUEbY7fbWbHi3zz6qKcXUIXExAbVbuHQoUMYDAYuu+wyNVRhZ8ZgMDBlyhQ2btzoEfLRHWd9t7vqaOzYsZhMJtVxXUpKCp988gk//vijS74HH3wwLHS1WVlZxMfHq3p4u91OZmYm0dHRqprM2ZZj8ODBamhMncjFYDBgsVgiYkeR/ktrY4xGI7feeit/+pN2+vz5TxETc9rlnMlkighz9bYiPz+f/Px8l3NKvAJnfK0dFBQUuBwbjUamT5+u2ifExMSQmZkZFoJAIScnh2+++UY93rhxI3fddRcmk0lTnZWYmIjZbO7IKuq0A7ow6ML06dOLqqpqUlM99b8Oh+tLn5qaysSJE8NCtx0q3n77bfbs2dOqa9xwww3ExsbSv39/ABISEsLO5Yfil0ihtraWx5oDa6empjJv3rwu/TvorFgsFvbt2xf2ITB1YdBOrFu3HrjS4/yyZfO5774/YzbLKpKqqiref/99lzzz5s3z6fitM/HJJ5+0WhAALrOKtLQ0zj33XMaMGdPq67YlK1euBOC8885j7NixLqohfV2g82KxWKivrw/7EJi6MGgHHA4HQ4cOYezYGgoKXN0oDx++QxUECkuWLOHAgQOddmveZ5991mauLbRISUlh0qRJCCGIjo4OSz27sy+g77//PqDIZr179+60v4muRKSEwAy/tybC2LoV7r0XTp7sy/HjcOgQgIG0tMFUVLiqKeLiapk9+z1A9r550UUX4XA4+OmnnyguLqampgaQdcXDhg3r4CfxjbvFqM1m00xz/q98Hj9+vLqToqCggAMH5JgONptN091yQkKCujVUuZ/y19jY6HJvkN2Mr169GpPJxKJFi9rwqdsOrW2FiYmJ5ObmeqyNKLPCkydPquE9k5OTufRS2XZFy3mf87m6ujo1VkKgzgGVz4q/KG/30QmeSAmBqQuDVjJ/Pnz9NYCrxHcXBAB/+MPfGDp0KDExMRgMBo4ePaq+aNXV1epiZ0v23reWhx9+OOC8/joj9zRnHzsNDQ1+r3/69Gm/eR566CH1HpFg8BQXF6eGSn3ttdc4ePAgp06d4ttvv/XIW1lZ6XHuxIkTvP76614FrvN/q9XKxo0b/boK1/rscDhYtWqVpmfUlhIdHc3ixYu71CYJZ5QQmOG+iKwLg1by+edQXQ1790oMHiyw2cBqhbKyE7z88ussWzZfzbt27RJmzaomNbWe2NhYF98/oe7QlI4qENrSN1FTUxN/+ctfvJbNysoiOjraw5Ppxx9/zNSpUwOuQziRl5fHwYMHmThxIvX19WzatEkz34UXXsj5558f9PXbwgK5I+/ZFYiEHUW6MGglUVGyS+vkZLvq2rq+vp4XX3yGbt3AZGqkqUneHrhpUwx33ikxefJGILgOuLNiMpl46KGHWLZsGdXV1R7phw8f9jh31VVXReyCq8PhUDcMKOoho9GI2WwmPj6eO+64Q50hhpPzQp3WYbFYwv771IVBOxAbG8tDDz3Eu+9WM3ZsFV98cWYG8NVX5zNx4hdER9s4evQoffr0CWFNO54dO3bw3nvvBV1OCEF+fn5Y61wDwWAwcNNNN7F8+XL13HnnnceECRNU9xo6nQ+LxUJVVVVIVMCBoguDdmD9+vX85z97ee65OzXT58yZTd++PcPKD00wawahYP/+/Zw+fZrExERsNhs2m40ZM2aEjaO/QNi9ezeffvopaWlp5OTkcPz4cYqLi/nuu+/47rvvQl09FxYsWEBiYmKoq9FpUHYUhWNQIgVdGLQDU6dOZerUqVx6KezYYWXRItdm/uyzT7jiijlh5biqPdcMdu3axdChQ9UFypKSEr788kvVm6OymyU1NZVevXpRWlpKXFwcP/zwg3oNq9XKkSNHSExMZObMmURFRbm43I4EkpKSqKioUDuEvLw8srOziY6O5uDBg9TV1REdHc2cOXNISkpSw38GSyjWDHR847y9NFyj+OnCoI34618zeO01rRTPJj55ssqvX56OJhxmBqWlpRQWFtK7d28MBgPnnHMOUVFRGAwGMjMzyc7OVncmReKoNTMzk6VLl7qE/Ny1a5eH0eFTTz0VkvolJiaSk5OjC4N2wDkEZo8ePUJdHU10YdBGFBcHtm3urbe2cPXV4bdw7G9moASRAVllM2jQIDVNax+687bFffv2MXDgQGw2Gx999BGNjY0uWxnLyspcyh47doxjx465nFPiACQlJWGxWLj22msjdv+7wWBQDeNGjRpFXV0dH3/8sdf83bp1cwkTqsyk3G0FlP91dXXEx8d7tSUAPLaO5uTk0K9fPw7JhjI67YCyo0gXBp2cp58+ytChyQB89x24+WFTeeuto+za9RFRUQ7y8/ODDiISKgwGg6rWKikpCcpdRlJSEsnJyTz++OMeBmPBIITg6quvbnH5cOPQoUO89tprfvf0a9kd+GPy5MmMHz8+6HLhqsLoDFgsFg4fPszZZ58d6qpooguDduC88+CZZ2D1avj0U9e0Dz64jKYmM2PGfE99fX1oKhgCzGYzDzzwACC7b37iiSeCDqsoSRIPP/ww8fHxOBwOUlJSmDFjRkBGauFGVVUVBQUFXgVBz549ycvLIzc3N+hr63r/8CQjI4Nt27aFbcxnXRi0gro6WLMGGhuhuDiZHTvAaJT/MjNhxw7PMnPmvMu4cUe5++7wUxV1FEajkUWLFrF8+fKgYiMrKEJEUa9YLBby8vLauprthtVqZcWKFR7qMWeOHz/Ohx9+yIcffuiRNnHiRC644IJ2rKFOe5Ceno7dblfdzoQbujBoBf37Q2mpchSYvcCBAwPJydnDE088ofrdMRgM2Gw2PvzwQ2bPns1ZZ53VXlX2SqgWkJX99SB38uXl5fzyyy8A9OnTx8XB25IlS4iJifG4Rrgb87gTHR3Nbbfdxt///neqqqrU87GxsWrAHrvdTs+ePcnKytK8xueffw74jvkQKcycOTPUVegQFI8D7lHvwgVdGLSCu+6C++8PrkxpaTpWq1V1JKZw6aWXkp2dHbZhG30xfvx4n7rmsrIyNbazs58i50VRo9FIU1MTcXFx9O3bl9jYWAoLC0lISGD06NEYDAaOHz8e8UZnCgaDgbvuuguQ1w7eeust6uvrqa+vZ8yYMcTHx1N6ZqTBoUOHOHjwIL169XIJnzhhwgSXxXqQty+mpaXR2NjI5s2bO+6hdHyixNjQsrQPB3Rh0Aruu0/+A3l0OmTIUA4cgKIiaGoCsxni48HZrf7Ro5mA7J540qRJREdHk5ycTElJSUjd23akbyJnT6Q2m401a9awc+dOzbKnT5/mq6++AuDLL7/0SO/fvz/nnXdewHUJNVu2bGHNmjWkpKRgMpk0/dUIIejXr59qZxDsrinn9p4xY0aLy3ZEOYCtW7e2qFykYTAYyMjI0GcGnZU33oDf/Q4g8BdBUXVs2LABm82G0WhskUOy9qK8vJzi4mLVOZwyklcoLi728EDp7P3SYDCoMwWHw8Hx48dV9YdyTvnc2NjoVRAEwsGDBzl48KB6fMkllzB69OgWX6+9UUaFWqPDgQMHkpeXR0NDA4WFhRw7dkxtR19eYd3PlZSUqDPMYD3MVlRUuKjmvOG+8F1RUcGRI0fUNPfv2mQydTnXK1qkp6dTWFgY6mpooguDFmK323nssb/zwAPz/Wd24/HHb+Xmm1/iN7+ZQt++fTGZTJoO2TqSUBqd/fa3vyUhIUHtlH755Rd13UTZH68IpqioKNWLqSKkImn3zOTJk8nPz6e6ulq1NfDmVvrEiRMkJSVRVlaG3W7ngw8+CFm924uYmBguvvjiUFejw0hPT2fr1q3YbLaw28arC4MWUlVVxbnnDqNv32MUFfX2X8CJysru1NZ2UxcBe/bsyYgRI9qhluFD3759vaYpi6DOC8nOoTCV83V1deruoylTppDvzZgjzLFarbz44ovq4rDWqB3kdti9e7c64o50brvtNs3zoR4IdSQZGRk4HA4qKipcXNiHA7owaCFRUVEkJCSQlRVLUVHw5WtqzFRVydPxqqoqqqqqImZ02xKKfDTS9OnTVT9DMTExHD58WF0odlZHWK1WCgsLMRqNEeWgzp2YmBhmzZpFbGwstbW1HoFqDAYDDoeDdevWhbWXy2D53//9XwCP7y/cYlW3J8pGitLSUl0YdBaSk5PZtm0bw4bt4ptvbgm6fJ8+R0lOTmb+/PmqqiOU9O/f30X33p6kpKSo6hC73c4XX3yhHs+cOZOMjAz69++vWXbw4MEdUsf2JCYmhhEjRlBQUMCWLVuAM0LPeWagxDhwTteKcOa+My1cmTp1KgaDgQEDBqjrZgaDgb1794a4Zh1HQkICMTExLjvFwgVdGLQQg8HAHXfcwfDhLSt/7FgvYmIOhU3ErtYKgtzcXFUHmpKS4rJFtri4mJ49e3qU0bK+tVgsnDhxolV1iRTGjh3L0KFDXVRicKZdSkpKPEaPWkLBPa20tJSff/65ParcKtavX695vqvYGSgkJyf7NDgMFbowaCXvvQctGazm5+dz663XhY2ztfT0dJcfaFycHMO5rq4uoPLK9sCUlBQGDBjgklZdXe3hk8jbc6empgZc585AU1OTurPIvU3q6+tVVxtK2qlTp9i9e7d6bLfbO8wDbkZGBgaDgaysLCZNmuSyG6mwsBAhhIdjPC1Hee7P2VW2liqkpKTowqAzYbfbKSws5KmnGoCRQZd/++0TVFb+lREjRgS9F7w9cP9xBioEFEaOHInJZKJ3797k5OS4pLXELqGrkJGRwbRp0zTTtNpt//797NDyc9IBKKqNkpISVb3VFnTFmcH+/fuxWq3qjrJwIHxqEuYoztUU/ayyyLdp02Utut499+Rw0UU5Ybe9rHv37i4Wwsr/yspKTQExZswYF5VQdXW1uktKoaysLCgfROXl5ZSWlqrXycjI4JprriElJSVsZlJthdVqZdWqVezevVuNhTxy5EhiY2Nd2k1RAQXjlC8xMVEzAJC7qsl5O2tDQ4NqQ6KV7v5Z+W+1WomKitJM70oOGQNBWTMrLy/XVJ+GCl0YBIjRaGTBggUe5x98EJ59Fo4eLWHChB7YbGC1wvff/8Devb9gtxux24188MEZofHZZ3DBBeHldqIjLZC9oXQeTz/9tIvOu7S0lGXLlnktN2TIkIjbifX00097GJ717t2b3r17M378eBITE722m7+RdCginekWyIGTnCy7ui8tLdWFQWfCaIR582D37hMMHXomaMXw4UmsWCFb1kZHR7sIgzDcSBAWGAwGGhoaGDRoEPHx8Xz99dd+y3Tr1i3ibDQ+/fRTTQvkI0eOEBMTE5FR3HQCR/mOw23dQBcGQWK1Wl32fivT4vr6ekpLT3H99bFs2mSkoSELu30RDocBu93oco2rrpJjHnjZPdnlOHz4MC+//HKLylZWVvLOO+9w2223qY7vFPWWolKKj48PL92sj7rs37+fLVu2cM4559DY2KiqWNzVY76OrVarx6KyVn7nHUnuLkd02peMjIzOLwyEECnAMsAM/B14CDgIvAwUAa8BJ4FbgFjnY0mSwtODkxMFBQUUFBTQ0NDgkVZR0Z21a3/ffORbDTRgANTWQvOmnZDTUncUQggyMjIwGo2a6Ypu1N+awcSJE9XPzl5OnbHb7ZSUlFBUVERSUhJmsxmbzUZdXR0rV650cYBnt9tVA65Zs2Z57HAKJe6d7pQpUzAajWrnbLVa2bZtG8ePH1djN3z77bdt7gffWWDefffd+oykAwlHH0XtMVw6F3gQSAEuBRZJkvQzgBDib8AiIAb4HdDP7fjZdqhPmzJhwgQmTJjA3r17WbFihUvaU0/dzrJlDo+X3WqFH3+UOHRIcMUVZ86bAgub3CG4rxm88sorPq2GFSRJQpIkl+AyI0aMcPHD315rDC3NH2qU35A/nJ9Lcb1RUlKixkCw2Wxs2LDBZ1jMBQsW6J18GJKRkcGmTZtobGzEbDaHujpAOwgDSZI2CCH6AUuAdOAsIUR34E6glyRJO4UQZmAukOF2HDEMGjSIiy66iA0bNqjnHn30UZc8druBdeumIUnZVFcLl7TU1HpMplhiYqzcccc2hg//2WWkMHDgQIa31KKtBbR0ZjBhwgQPD6ZRUVGqKk0J3KPYGfhTd+j4pkePHmpAdZvNxunTp3E4HKxdu1Yz/xNPPBHU9bOzsyNKsEYqijFhWVlZ2HhzNfgLxt0ShBCJyJ3/UUmS3hBCnAdcDPSXJOkWIUQc8DRgdDp+SpKk/+Ptmlu3bm37inYycnNzW9Sz6m0bGC1pX71tA0P/7bYfgbZtmwsDIcQc4HMgCfhBkqTuQojpyJZZdcDXyGqhScAp52NJkv7cppXR0dHR0QmI9lgz2AO8jtzxPyWEWA1Ygf8GmoDngETgVqDB7VhHR0dHJwS0i5pIR0dHRyey0N4PqKOjo6PTpdCFgY6Ojo6OLgx0dHR0dMLUHYUQ4logVZKkZ9zOZyBbLHcD3pUk6W9u6WOAJ4H+wHHgckmSDruVXw0MAiRgrlv6GOARoC+y1fRdkiRJbuXfbC6/GXjQLX0MsLS5bDZwoyRJRW7lFYtrc/P1i9zKPwlkAt8Df5UkaYtb+dXNz7cGeMEt3fn+Lytpbvd1sfR2T0M2nX4LGAYUADf4yT8U+DMwuLlNf+1uSe78fTqVbwDigUr3OmmUGePrHs1blx9vPry2+dq+nsE9f28/13e2qn+4uc4e7RlMOweR5lJX59+rv7JOeVYCC3z8FrXa/xzgt5IkLdG4nq/6jgHuAhKAOyRJ8gjgHMD77dcjQQD1vxZIBd7WyufWrr8DXtDI4/y9Lwfu9XG/lcDfkH8f7tdxfi9f1crTnO8c4LfIfYC39r0GmIA8kM8CDF7up3wHv9d6Ni3CThgIIUYAFwDbNZIXAf9A/vLO00gfBbyP7AZjCjCtOb9z+feBL4HrNdJHIf94Pka2oJ4NPOZW/m3gGDBOI31Uc55ebtd1Lr8IOB9Z6GjV/xNgpWK1rVH+E2An0MNZEDjfX6Oscl8tS2/3tH7IbkNeAHIDyN/UnH++1vU1vk+l/GLgZ2BVAGVG+boHcic+DdmWZXrzsa9ncM9v8HN9Z6v6OciDEa32DKadA01T6qr1e/ZXFiHERUAenvgsB/wRuFGjnL+yvwHuAEYAl2vUx9/7HahHAq95m+8xEdjhI59zu25Dfidj3fI4f+87gOEaeZQ2PhdY6OVe6nup4YXB+Rn/CNwEPOAtjyRJbwFvCSEuBf4HmKVRJ+fv4Edk4eFRb3fCTk0kSdIO5K2pWvSUJGkV8AygFRE9BxDIo8IL8PzB9ZQk6a/AVuSGck/PAcYiv3APA+9plP8nUIj8Jbin5wALkP0wrdKqP/LW2/Obr6FV/6uANUKItUKI3hrlewAzgLs00nOAe4QQ77ml9ZQkaSfwE+DupMc9rSfyIGE18ojIX/5RyD+2R4By9/wa32ev5vJRyD9QjzpplMnxc483JUmqR34ht/l7Bo38/p5hA/K7sgT5e/XWnsG0c0BpTnXNQ7sD9VpWCBGNPOjRep98lfs18qzqlebRasBlgfXIHc7vkWevLgTwfnu7bsB5ne5h8JbPrV23Ng+g3PM4f+/vaeVxauPXgAYvdcoB7m7eZj9KK49zmwMDA2iHy4HNWnXizHcwD1jnJY8HYScM/KBEyTAiv+TuLJQk6SbgX8BESZK+91I+G+imka6UfwzZ9sE9hJlSvgp52uWevhB5BLFUI00pfyfylFPLM9FC5A7qSuQp5Q0a5RcCtzffxz3duf43uJUDuQNO8PJMUcjPXN/8pxz7y78I+Kz5ntdr5HdHiZDT0JxXq07uLPR3DyHEYCCz+Tv19wzu+QN5hgrkGdkovLdnoO0cTJp7Xd3xVfZO5IGTPchyc5FnTLcC9wVZ9r+QZ8/vAuM1yvrCZzsEmdfgL5/SrsjvmrdrVSDPYEd4yXMnsibCjqwR0MqzUJKkm5F/Xxd4yaO0+S3IM0+vzyaEEMA+oNZLPuU7+H9Avq9rORN2aiI/7BRC5CJP8ao00ucKId5AHp3/x0v5schfykov5SuQG+5HYDRyozqXX4QsZb/SSJ+L3FkkN9fxTmR1iFoeuAgYgzzCd09X/DMdRp4iugYOlssvAT4CajTS5wIvIs+anNOUdotpvjZe0oqQrcIzgauRf+D+8l/TXC/lmdzzu6OUr0YW6DkBlJnr6x7N7kwe54wA3OnrGTTy+7u+YlX/FvL021t7BtrOAadp1NUdX9cdhbz+lY/cESwOsFwlcLQ5rTHIe56DPMNKQlZ3vOGl3sE+S0vzauZza9frveRx/t7/KIQYjTybdb6fcxvnecmjvJfdgYe81Nu5zXO9XEdhBvLof7yXaynfQTLwvJ9rqYT1zEAIcZEQwtm94yvIet2rgW800uuRO/FsYIAQ4tca5ZcjN1Y/jfR65FH3MuQp2z6N8ucjLwIJjfR65EUrO7ACKNQoX4X8xWil1yNbav8TsAD7NcpPQpb43TXS65unojfhqltW2m2RlzopaS82H2chq8kyA8jfgKx7X4i83uCe3x2l/CBkR4ZadXKn3s897m2+1uNCiMuR1XS+nsE9f52f6ytW9U8gt6t7G2i1jb9202oTrXLudXXHa1lJkq6XJOkeYC2eumJf91yO3Im/jfwuBHxP5Hfjc+S1OZ+CwMf7rdVGvurgLa8D+beglU9tV+TRv1Ye5Xv/H2R10t3uedzaeLRWHlzfy395uZdzm8/0ch2FMZIkbfXRBq8ifwerkT1H+7qWim6BrKOjo6MT3jMDHR0dHZ2OQRcGOjo6Ojq6MNDR0dHR0YWBjo6Ojg66MAgYIcQcIYTvKPeeZaKFEN00zs8VQkS1Xe10dHR0WkeXEQZCiK+b9xYjhPhSCOHPqMW57FCgryRJ9UKIz4K47W+9nFcMnXR0dHTCgi4jDIANwPlCiGSgTpKk00GUnYscsxnkvct+EUJYALskSZXuaZIk7QOig51p6OgEg5DZ0+zqACFErBBiZbMvHZ12QggxSwjxgBBiihAiYiI4RpoFcmtYh+xYLgb4RAgxEdk74CfIHhT/j4+yTZIkqUJACNEf+Auy6fguZP82xcBZkiRd3JztNzQb+gghfo9sadsNuFmSpApk47ihzf913BBCTJAk6ctQ1yOSkSRJEkJs4oxPo0zgIUmSdoewWl2BJODfQCmyQ8uIoCvNDDYhm45fiGwtCPBVs4veBiHE2T7KGpw+90B2B/FI8+yiGNkPyzZkNxaKK9rdkiQpLiFGInf6jyC7kVCuqVv8ecfgP4tOABxEtsY3AEOdBYEQIk4IsSBkNeu8COSB9nxkq/WIoMsIg+aRfQlwrpOLZ7PTf18ds1EIobRVHLJn0T80H9e5/QeYJknSx07HzyGvEyxBdiIFsun6rmCfo7PTPMWeDQwTQvx3qOvTCTgAnIXsqtn5N4kkSXVoO7HTaR1WSZL2IM/ILgx1ZQKlK6mJQHbu5OwLJVcI8SKAJEm+OuY3kJ18/V/goCRJPwkhioQQF+ImRIQQ0/F03Xsx8gsZBxxpnoXUSJKk5QSsyyKEuAx5gFLCGd82L4S0UpHPL8jO9TZJktQg5IAt5yP746kA4oQQ90qS9Jivi+hoI4S4B9nP1kXIqqF3gb7NfcPZnFlrDHu6rG+i5jWDiZIk/SnA/Bcjq5Xq/eRLkiSpxk+eyyRJej/w2nYad+UCAAABCklEQVQNhBD/gyxwz0FW6z0vSZI3b506ASCE6AM8IEnSbc3H8cjrXSuRO7H/IHvPfcR5XUzHP0KIWUCMJEn/FkK8IUnStaGuU2vossJAJ/wQQgxB9oEfj7y28kHzYrtOG9HcgR1F9gQ7VpKk5UKIB4A3JUn6JbS1iyyEEE8jR0IzAA9LknR3iKvUKnRhoBNWCCHMwEhJkjaHui46Or4QQkxCjh2SihwJLaIHLl1tzUAn/BmKHFlKRyeskSQpGAPUsEefGejo6OjodJ2tpTo6Ojo63tGFgY6Ojo6OLgx0dHR0dHRhoKOjo6ODLgx0dHR0dNCFgY6Ojo4OujDQ0dHR0UEXBjo6Ojo6wP8HSutZrXoOf3EAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, axes = plt.subplots(ncols=4, nrows=1, sharey=True)\n", "axes[0].invert_yaxis()\n", "\n", "ppp.plot_multivariate(\n", " axes,\n", " well_cug1,\n", " vel_log, por_log, vsh_log, obp_log, 1500, 3500, a0, a1, a2, a3, \n", " well_cug1.params['bowers']['B'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2. Pressure Prediction with multivaraite model" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Multivaraite pressure prediction:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "pres_log = well_cug1.multivariate(vel_log, por_log, vsh_log)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Post-process predicted pressure:" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "pres_log_filter = ppp.upscale_log(pres_log, freq=20)\n", "\n", "pres_log_filter_smooth = ppp.smooth_log(pres_log_filter, window=1501)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "View predicted pressure:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAHyCAYAAADFr1rYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXlcFEfeh5+e4b5vBAUE0QbviKDG+4zGJBo17ibmjtmcJrrZ3Nn1zW7y7q45NrsxeXOZjckm2VzmMLf3EY0IURHFNiogIoKC3PfMvH8MMzIwwAAzzAzU48cP09XVVd+e49vV1VW/knQ6HQKBQGAtVPYWIBAIehfCVAQCgVURpiIQCKyKMBWBQGBVhKkIBAKrIkxFIBBYFWEqToAsyzGyLM8zk37YTNo6K9Z7iyzLk7uizdbIsvxPWZYXWZh3taV5bYW5z8qCY8bIsrxbluVdsizf3pQWL8vyflmWt8iyHGd9pd1HmIpzMBC40kx6q0FGiqLcYa1KFUVZryjKrg6yDcS8Npshy/IfgUk9WacV6MqAsDuA+cB04A9Naf8E7gJuBv5iHWnWxcXeAqyNLMurgcnozy0L+BvwKpALuCmKslyW5QTg/wBP4E+Kovwoy/IU4O/ojfZpRVG+lWV5KfD7pqLvAy4AaxVFuVqW5anAYkVRHmi6Cv0LuFdRlMtkWfYC3gNCgR8VRXnGjM4Y4EegDFCAZOAa4DNFUUY07V8LfNlUd4gsy8OA3yqKUtTO+R9WFGVE0+t/A39WFCVbluWvgRuAhKZydcAriqK8K8vyLUAiMBL4SlGU15q9l4cVRdnQtP0/wFygEfgNejO5t7k24DywDogFjimKck8bOsOA9UAA8K2iKH+RZXlb03taYjgPWZavBx4GioFK4DrgbUBq6z1oKt8V+BTwRf+ZHpZl2QN4FwgH0oDHgE3oP/c7gM8Bn6bz2N70eexSFGVly++CoijpsizfDdzS9H7cpChKjizLbwBDgVrgGkVRqs1o6/D7AaAoyn1N+f2BiqbkGEVRfmlKD27n/FcCSwEtsAz4H8AfiALOAYWKoixv6/ju0FtbKqmKokxD/+XxB6YBHzV7E/+G3vlnoX+zAZYATwMz0H9JQP9h3ApcC7g2pZm74ngAIYqiXNa0/Ttgj6IoU4CZsiwHtaGzCliD/gv4BRDfonydoihvASuBzxVFmdGeoZjRtwGY0/RjalQUpRz9F/k6YDlwf7O8S4DbDYbSBoXABGAv+h/Mm2a0LQCqFEWZDvjKsjyyjbIeR2/QEwAfWZZbfhcN57EU/edwEvg/RVEaFUXJpwNTARYBGYqizABqmtJ+B/ykKMpUwA+YCDSgN4Gz6A33QFPeEmAc+s8e4EaafRdkWQ4FbmzS/wdgrCzLvsAuYCr6C9DUNrRZ+v0w8BTwj2a6DLTXKLgZmIPeLH3Qv5/voL+QrkJ/AbEJvdVUMpr+ZqFvJu9TFGVHs/3xwHPARsCrKe1l4HbgEy6ZyrPoDeg1oOUVx73Z60bg+WbbQ4Drm668QUBkGzrPoL+SGP42/6G4mz2ic2xC32qbDHzXrNzXgAfQm6GBdxVFOddWQbIsS+i/iF8Aw9F/Oc0xBJguy/JW9K2fts5dBn5pev2YoijaFvsN538AfasyFP0P1lIGApnNypCa6kxtSktr0pcLjEBvAqMBQ9/HBkVRGtG3JAGewfS7EAscAVAUZR/wGfrPcC76FlJ/2n+PLPl+IMvyNGCgoij/bUpqftHwa+s49Eb0X/S3SJVNaflNGvPo2JS7TG81FUOLYTiwDX2LoDnHgRuarqZvN6XNBG4D7gRWN6VNURRlIfrWxEPor2qGq8q0ZuU1KorS0GxbAZ5rKv9Z9E13c5j7YH2a/k5vllYDeFt4vDFNUZRaoB6Yh/42CvQts6XACy2Oq6R9xgBhiqIsAA61o+048F5TC+Fx9D9acxxrKhNgiyzL4U1aA5tuT8Ob9o1VFGWaoiiLFUWpMVdQG+ShNw3Qm6EO/UVmXFNaMnpTOIjegLKB6GZ1GP4a3s+W34VTwDBZlqWmH/5a9IZyTlGUa9H/gJvT/LOy6PvR1BH7v+hbGwYqZVn2b2p9upo7ronhiqJchf5zt1o/myX0VlMZKcvyduA0+qtKy1uWJ4H3ZVn+mUsfTB6wBfgKfTMRoLQpz7+ADxVFOQuUybL8CvrbKgMty38LuE6W5V3o78vbumXRmfm7pekJzsBm+Q4AobIsb5dl+bp26jWXtgUY2awVsgF9f8H/oL89ae870Lys40CELMu7gWFcurq21PYlIMuyvAN9ayivjbL/Djwgy/JO9P0KheivrC8DN6G/HQE43/S04ztZlu9q5zxb8hn6W5JtQGBT2lvAxCZtZYqi7GzSnwP8SlPLo0XZhtctvwsXgA+An9C3CtagbwXNavrcW7ZAmpdp6ffjf4EQ4EtZlrfKsuwOvIm+7+dT4KN2zt9DluWf0N/qfGFmv81mEku9bZayLMt/AjINnYuOQNNV+L/oP0ip6e83iqK0bC30Spp+2M2/aOmKojxswXESepN3Q39x+FVRlEfM5PsH+lsXw/vbqCjKbGto7wk6+/2QZfkyILTpAcMt6Pt6mh+7QlGUI+aO7Ql6nakIBAL7YvdHys0eLZYDdzY9oRAIBE6K3U0F/RiEh9H39t+MvsNLIBA4KY7QUdtPUZRM9I/yYu0tRiAQdA9HaKnUNv1VY/6xKQDp6emi80cgsBNJSUkWj2txBFPJlGU5Cf3tT1uPHwFISkrqGUVWICsri8TExI4zOhDOptnZ9IJzak5PT+9UfkcwlX8Dr6Af9PU7O2sRCATdxO6m0vS05yZ76xAIBNbBETpqBQJBL0KYikAgsCpOZSpaXcuJrAKBwNFwKlPJLMrsOJNAYAcef/xxKivbn+idmprKxo0be0RPRkYGX3xhbh6h7XEqU9mavdXeEgSCbiFJXQtj0pFptdw/cuRIFi5c2KW6uovdn/50hvNV5+0tQSBok3Xr1vHrr79y5swZXn/9dc6ePcuhQ4dISEjg3XffRaVSERMTQ2pqKhs2bKC2tpY//elPFBUV8eqrr6JWq3nssceorq5mzZo1uLm5sWLFCg4fPkxmZiYvvfQSS5YsISQkhKeeego3NzfuuOMOTp06ZbI/ISGB1NRUsrKyuOWWW6isrGT16tVoNBrmzZvHFVdcYdP3walMxdvNm4s1Fwn0DOw4s0DQgu0529mes71V+rSB05g2cFqH+dvKZ2DmzJlMnz6dXbt28d1333H+/HluvfVWnn32WV566SU2b97MkSP6iASVlZWsXauf5vbss8+yZs0aCgoKeO+995g8eTI+Pj4sXbqU/v37Ex8fT2pqKitXrsTHx4eKigqWLFlCamoqe/bs4Z577jHZb8DQKvrqq69YtGgREydO7PRAtq7gVKayUF7Ia2mv8fjkx+0tReCEdGQK3c0/cOBAzpw5g06n4/Dhw3h4eBAaGookSahUKtzc3Ix5U1JSjK8lSUKr1WIIQxITE8Ptt9/OV199RV5eHosWLUKtVqPRaADYsmULubm5TJ06lQMH9CF1m+9vSfPwJllZWTYfme5UplLTWMP+s/up19Tjpnbr+ACBwE5ERUURHR0NwHXXXccf//hHqqqqiIqKatWvcvvtt/PYY48hSRKPPfYYFy9eZN26ddTW1jJz5kwApkyZwlNPPcW0adOIiIhg165dnD17loCAgFb7Fy9ebFL+Nddcw+rVq/n4449tfusD6F3MGf6npaXpGjWNut999Tvduwff1Tk6R48etbeETuNsmh1Vb2pqqu7+++/X1dTUtNrnqJrbIy0tTafrxG/VqZ7+qFVqQr1DOXjuoEmTTiBwJJKTk3n55Zfx8PDoOHMvxKlMBSDAI4BxA8ax+dRme0sRCARmcDpTmRQ9CU8XT344+YO9pQgEAjM4nakkRSSh1WkZGjpUjLAVCBwQpzMVV7UrCxIWcMOIG/jg8Af2liMQGPn888/ZsmULAJs2bWLv3r1m861du5Zjx47ZVEthYSGvv/66TetoC6czFQMeLh4EeARwrrLNlToFArsxe/ZsJkyYYNy2ZG5Qe3RkRC33h4eHc9ddd7WZ35Y41TiVltx+2e28kf4GT0x+wt5SBAITPv/8c/z8/Jg5cyaff/65cRi9YeDZJ598QmVlJRMnTiQ5OZm//e1vSJLE3XffTXh4OE888QReXl4sWbKEmpoadu/eTX5+Ptdeey2jRo3iD3/4A66urlx55ZWoVCqT/SkpKeTn57N+/XqeeOIJNBoNf/zjH2loaGDkyJHcdJNtY6I5tamEeIVQ3VBNVX0V3m5txswWCAD44tgXHDx3sMvHj+43moUJnZ+kd+211xqH0efl5XHy5ElSUlJITk7mueee4/jx46xcuZJ+/frx1FNP8eCDD6LValm8eDHDhg3Dz8+Po0ePMmvWLBISEqirq2Px4sVkZWWxfft2nnnmGZP9BgyD7Hbv3s3o0aNZunQpaWlpXT5/S3FqUwG4edTNvJfxHnePvdveUgQOzsKEhV0yBWvQchh9TEwMrq6X1lfX6XRotVokScLX15eVK1eyd+9e0tPTWbFiBSqVCq1WH0/o0KFDbN26lWXLlvHOO+8AmOxvj6ysLMaOHWvdk2uB0/ap7MrdRXVDNUOCh3Cy5CQarfl5DwKBI2AYRm/oyG3OsmXLeOmll3jiiSe44447qK+v56233iIzM5P4+HgAJkyYwL/+9S/eeustAgICKC0t5YMPPqC6urrV/pZMnDiRAwcO8Ic//MEi4+kuTrOWcnp6uq75RKht2dvQ6rTMjJvJrtxdXKi+wLWJ19pRoSnOuBSDs2l2Nr3gnJrT09M7te6P07ZUxg0Yx89nfgb0A+J2nd4lhu4LBA6A05qKl6sX1Q36pp8kSUyNmcrO3J12ViUQCJzWVADcXdypbdSvmnq1fDUbj/dM/E+BQNA2Tm0q4/qPIzU/FQCVpGJ42HAyCjPsrEog6Ns4talMjplMbECscVsM3Rc4A82H80P7Q/q7gj0j6YOTm4qXqxdR/lHGbTe1G5G+keSU5thPlEDQgvvuu6/d/S2H9HeXlpH0e2KuUXOcfvBbS24bfRvP73mep6c/bW8pgj7CddddR3BwMOHh4eTk5PDKK69w5swZNm/ejI+PDydOnGDNmjXcfPPNZo9vPqT/vffeY9SoUfTr149///vf3HHHHcbI+cuXLyc4OJi//vWvuLu7c/fddxMYGGiyf+TIkSaR9Ldu3WoyhD8hIYE1a9ZQU1PDVVddxfTp063+fvQ6U/F198XdxZ0L1RcI8QqxtxyBA3Hs2DHOnev6BNR+/fqZDIM34Ofnx9y5c3F1daW4uNhk36233kpaWhqPPPIIABcvXmy3jjlz5vD+++8TERHBFVdcgbu7uzFy/k8//cQ111zDuXPn+OSTT5AkySSy/k8//cTIkSOBS0P0Z8yYYTKE//3336e2tpZ+/fqRmZkpTMVS7rjsDtb9so5HJz1qbykCByIhIcGsKXQXQ9hId3d3ysvLqa+v58yZM10qKzw8nAsXLlBcXMz111/PF1980SpyflJSktE0zEXWb0nLIfxXX301l19+OZs2beqSxo5w6j6Vtgj3CaeyvpKq+ip7SxH0AZpHx09LS+OVV17h559/NqalpKTw0EMPsXWrZStsDhkyxLh+T2hoKKdPn+arr76ioqKiVd6O9oPpEP6rr76aDRs28NBDD+Hn59eZ07SczkTJtuf/pojeZnlhzwu6Bk2DSdqJ4hO6tfvWWhYu3AY4Y9R0Z9PsbHp1OufU3Kuj6bdFbEBsqyntg4IGcbrsNA2aBjupEgj6Jr3CVCbHTGZX7q5W6b8Z/hs+OvKRHRQJBH2XXmEqIV4hlNSUtEofEzGGAwUHxERDgaAH6RWmAk3r0epax4qYN3ge3534zg6KBIK+Sa8xlRFhI1AuKK3SZ8bOFAuPCZwKS0fA5ufn87//+789oKhz9BpTWTJ0CYmhrYPfSJLExKiJ7D692w6qBALL6Gq0/ZaLvTsCvWbwW3tv7rWJ1/Lwjw8zKXpSDyoS9CWys7NZs2YNrq6uVFRU4O/vj7u7O5WVlTz11FPGaPlXXHEFAQEBJtHzs7KyjNH2lyxZAphG2x8zZgyvvvoqlZWV3HnnnXh4ePD888/j7+9PYGCgnc+8Nb3GVNpDJakY3W80vxT8wpiIMfaWI7ATX3wBB7seTJ/Ro2FhG3Gzi4qK8PHxYenSpaxdu5brr7+eQ4cOkZGRwfvvv2+Mlv/AAw9w5MgRk+j5L7zwgjHavo+PD5s3bzaJtn/s2DHc3d3x8vLi8OHD5Obm8uSTT1JaWso333zT9ROyEX3CVAB+O/y3PLn1SWEqfZiFC9s2he4SHR3N7bffzldffcWRI0dQqVS4u7sbnzwaBoYZ0DWLng/6ofRtRduXJIkbb7yRwMBAMjMzyc3NxcXFxSQavyPRZ0zFVe1KfFA8WeezzPa9CATdobS0lHXr1lFbW8vw4cON6ZIksWzZMv76178iSRILFy7ksssuM7n9AZg6dSpPPfUU06ZNa1X29ddfz9///ne0Wi0rVqxg0aJFvPjii7i7uxuH8zsUnRl+a8//7Q3Tb07OxZw299U01Oge2/SYReV0F2ccju1smp1Nr07nnJr75DD95rz5y5ttDnbzcPEgwjeC7IvZPaxKIOg79DpTGRw0mBMlJ9rcf/tlt/P2gbd7UJFA0LfodaYydeBUduTuaHO/j5sPvu6+FFQU9KAqgaDv0OtMJcY/psMYtXeOuZM3f3mzZwQJBH2MXmcqkiS1enzXkkDPQFSSiuLq4jbzCASCrtHrTAVg2sBpxtUL20K0VgQC29ArTWX2oNl4u3m3myfcJ5zqhmrK68p7SJVA0DfolaZiKcvHLGfdL+vsLUMg6FX0aVOJ9o/mfPV5ahpq7C1FIOg19GlTAbh19K28c/Ade8sQCHoNfd5UhgQPIbs0m3pNvb2lCAS9gl5tKl8f/9qifDeOvJH/ZPzHxmoEgr5BrzaVfWf2WZRvZPhIjhQdEct5CARWoFebSqh3KIWVhRblXTZyGR8c/sDGigSC3k+vNpUpMVPYmbvTorxjIsaQUZhBo7bRxqoEgt5NrzaVEWEjOFx02OL8vx3+Wz7KFIuPCQTdoVebilqlNrsWUFsk90/ml4Jf0Gg1HWcWCARm6dWmAvr4KZ1hydAlfHr0UxupEQh6P73eVOIC4zqVf0LUBPbl7+tUC0cgEFyi15tKV1iUuIjPsz63twyBwCkRpmKGSdGT2H16t2itCARdQJhKGyxIWMBGZaO9ZQgETocwlTaYGjOV7Tnb240gJxAIWtMnTEWj1fD33X/v1DGSJDF/yHy+/fVbG6kSCHonfcJU1Co1FfUVnT5uZuxMNp/aLForAkEn6BOmAvqlOSrqOmcskiRxRfwV/HjyRxupEgh6H33GVC6Pupy9Z/Z2+rgrBl3B9ye+F60VgcBC+oyppPRPsTgUQnMkSWJm3Ey2Zm+1gSqBoPfRZ0zFw8UDV7Vrl46dP3i+xQGfBIK+Tp8xFYDHJj3WpeMkSWLqQP0jZoFA0D59ylS6wwJ5AV8e+1L0rQgEHSBMxUIkSWLawGlsy9lmbykCgUMjTKUTXCNfw1fKV6K1IhC0gzCVTiBJErPiZrEle4u9pQgEDkufM5Wq+irOVZ7r8vHzB8/nm+PfiNaKQNAGfc5UtDot6w+u7/LxkiQxe9BsNp3aZEVVAkHvoc+Ziq+7b5fmATVnXvw8vv31W9FaEQjM0OdMBcBV5dqthcMkSWJe/Dy+P/G9FVUJBL2DPmkqYyLG8EvBL90qY86gOfx48kfRWhEIWuBiq4JlWb4JCAA+AtYD5cCdgEd724qilNtKk4FJ0ZNYd2Ad4waM63IZkiRx5eAr+fbXb5k/ZL4V1QkEzo1NTEWW5ZHAVCADeLjpvztwMxDTwfZaW2hqTqBnIIkhid0uZ1bcLFb9sIorB1+JJElWUCYQOD82uf1RFCUDeBeQgH6KomQCh4FYC7Z7BGu0LiRJ4hr5GjYeF7FsBQIDNrv9QW8oALVNf9WATwfb3u0VmJWVZWWJ3aefrh9rDq0hXhNv0lqpra11SL3t4WyanU0vOKfmzmJLUzGQKctyEvrbm9NAZTvbee0VlJjY/VsWW3Cr160crz3OwoSFxrSsrCyH1dsWzqbZ2fSCc2pOT0/vVH5bm4oO+DfwKvpWye+AOuCVdradjmkDp7Hq+1VcI1+DSuqTD9QEAiM2MxVFUXYAO5o2b2qxu6Ntp+PaxGv5POtzFg9dbG8pAoFd6dOX1cyiTI4XH7dKWVNipohVDQUC+riphHiFWDVS/pKhS/js6GdWK08gcEb6tKn08+nXrRnLLZkYPZG9Z/aK1oqgT9OnTcWANYfaXzf0Oj458onVyhMInI0+byoJIQkcu3DMauVNiJpAan4qjdpGq5UpEDgTfd5UpsZMZUfujo4zdoIbRtzAN6e/sWqZAoGz0OdNJco/ivmDrTshMCkyCaVUoV5Tb9VyBQJnoM+bCuiNxdpcG3st7xx8x+rlCgSOjjAVGzHYfzAnS05S01BjbykCQY8iTMWG3Jl0J+sOrLO3DIGgRxGmYkPig+IpqCigsr7S3lIEgh5DmEoTtgoLedfYu3g97XWblC0QOCLCVJpYf2g9p8tOW73caP9oyurKKK0ttXrZAoEjIkyliQkDJrAjx7rjVQzcM/Ye/m///9mkbIHA0RCm0sSQ4CEoxYpNyo7wjaBR28j5qvM2KV8gcCSEqTRh68DV9yTfw2tpr9m0DoHAERCm0owInwjOVpy1SdkhXiG4ql1tVr5A4CgIU2nGtIHTOHTukM3Kv2esaK0Iej/CVJoxLGwY8wbPs1n5/h7+BHgEkFOaY7M6BAJ7I0ylh7kr6S7eSH/D3jIEApshTKWH8XbzJtI30mqxcQUCR0OYih1YPmY5b/3ylr1lCAQ2QZiKHfBw8WBw0GAOFx62txSBwOoIUzGDNcNLtsUto29h/aH1Nq9HIOhphKmY4dtfv6WkpsSmdbip3RgVPor9+fttWo9A0NMIUzHDpOhJ7MrdZfN6bhhxAx8c/sDm9QgEPYkwFTOMiRjDgXMHbF6PWqVmSswUtmZvtXldAkFPIUzFDC4qlx5bYmNhwkK+PPalzeK5CAQ9jTCVNvB186Wstszm9UiSxIKEBXypfGnzugSCnkCYShtcI19DWZ3tTQVgRuwMtudsR6PV9Eh9AoEtEabSBomhiUT7R/dYfctGLOP9w+/3WH0Cga0QpuIgJPdPJqMwg7rGOntLEQi6hTAVB2L5mOViSQ+B0yNMxYFICEngTPkZKuoq7C1FIOgywlQcjLvH3i0COQmcGmEqHfDxkY97tL5o/2iqG6opri7u0XoFAmshTKUDjhcf7/H1kO9NvpdX9r/So3UKBNZCmEoHjOs/jtT81B6tM9Q7FHe1O2fKz/RovQKBNRCm0gGXR13Onrw9PV7vvcn3igXIBE6JMJUO8Hbzpqqhqsfr9XX3JcI3QoSdFDgdwlQswFXlSr2mvsfrFWEnBc6IMBULuC/lPiRsu4KhOTxcPEgMSeSXgl96vG6BoKsIU7EAw+qC9uCmUTfx3qH37FK3QNAVhKk4OC4qFyZFT2J7znZ7SxEILEKYihOwKHERn2d9LgI5CZwCYSpOgCRJXC1fzcbjG+0tRSDoEGEqnUCr09qt7llxs9iavVUEchI4PMJUOsHT25+2a/03jryR/2T8x64aBIKOEKbSSezZUhgbOZbMoswen4skEHQGYSqdYGT4SDIKM+yq4a6xd/FG+ht21SAQtIcwlU4wOWYyO3N32lVDfFA8xTXFXKy5aFcdAkFbCFPpBGHeYRRVFdlbBven3C9CIwgcFmEqnSTMO8zeEgjzDsNV5UpeWZ69pQgErRCm0kkeHP+gvSUATaER0kRoBIHjIUzFSfF19yXaP5ojRUfsLUUgMEGYihNzx2V38PaBt+0tQyAwQZiKE+OqdmXcgHHsyt1lbykCgRFhKk7OkqFL+PTop2KyocBhEKbSBTRaDUfPH7W3DABUkoqFCQv5UvnS3lIEAkCYSpdQq9R8ePhDe8swMj12Ojtzd9KobbS3FIFAmEpX0aFzqFuOW0bdwvqD6+0tQyAQptJV4gLjyC7NtrcMI6P6jeLXkl+pbqi2txRBH0eYSheZEjOFHTk77C3DBLEOs8ARcLG3AHuj0+k4cuQIOp0OtVqNJF2Kmp+ens7JkydbHaPRSBSXBHE49GOiS6Lx8PBolaeoqIiSkpJWdVmLyMhI4uLiTNIGBgykoq6C4upigr2CrVaXQNAZ+pypaLVaXnjhBaqru3ObIJF1dBhxU+NIS0tjxowZJoah0+lQqVSoVCp0Oh3+/v4MGTKk++KboVarzabfl3Ifa1PXsnraaqvWJxBYSp8zFZVKxcMPP9ypYy5evMiJEydQqVR8/fXXqNVaGhpciSWW2tpavv32W7PHjRo1ChcXF9RqNa6uPbPER4hXCD5uPuSU5jAwYGCP1CkQNKfPmYolVFVVUVBQgE6nQ6vVkp2dzb59+0zy+PhUUFwcSHCw+bgmAwcO5KqrrsLFpeff4nuS7+HPO/7M32b9rcfrFgj6vKnk5+fT0NBg3JYkiczMTNLS0to9LikpnU2b5jB37neomnV3S5KETqcjJyeHHTt2MHPmTFtJbxMvVy8GBw3m0LlDjOo3qsfrF/Rt+rSp6HQ69u7dS2VlpXFbp9ORl9dxnBJXVw1jxqSze/dkJk/ehVqtYvDgwQwcOBCVSkVhYSH+/v5Gc2reAdzytSRJqFQq41/D/+bbLfN4enoSHNx2Z+wto2/hkU2P8GK/F7v69ggEXaLXmsqxY8f46KOPbFpHv35F1NR4smfP5UycuAdFUVAUpVW+MWPGALTqzG3+WqvVotVqTV633Da81ul0hIaGsmDBgja1uahcmBIzha3ZW5kRO8OKZy0QtE+vNZWEhARWr+7+ExCdTkddXR1lZWVkZGRQX19vcmsUGnuEM5oi9u9PIjk5HYDAwEA0Gg0rVqywS5+KgQXyAlb9sIppA6ehksSQJEHP0GtNxVpIkoSHhwdVVVXs2bOn1X4vvGgfJIgJAAAgAElEQVSI/wbdziQSEycyalQUQUFBnDp1itLSUuOtjlarRaPRmNz6GDC0Wpr/9fPzw8fHp9vafzv8t/w387/cMOKGbpUlEFiKMJUmzp07R2lpqbEPQ61Wm/RnFBUVER4eTmFhIX5+fsTExKAoCvX19ejQERFxjh9/rCUr6yeraZIkiUcffRR3d/culzF+wHg+O/oZixIX4eHSepCeQGBthKk0cfr0aU6dOmVsUWi1WnJycszmLS8v5/Dhw8ZtL7yoV1XhprXuWBSdTsfatWtZtWoVKlXXb18MawU9MO4BK6oTCMwjTKWJlJQUUlJSTNLKy8uprq42uS1paGjgnXfeMckXQwyKdI7Y0CTgdLd0DB48mAEDBjBlypRuldOc+KB4SmpKKKkpIcgzyGrlCgTmEKbSDn5+fvj5+ZmkffPNN63y9aMfP+T15/6r/VmyRN85nJWVRWJiYo/otIT7U+7n5X0vi+H7ApsjTMUCdDodRUVFNDQ0MHr0aIYNG8b69Zdil+ScGsRQf3eOHNnCkSNbjOk+Pj4sXbqUqKgoe8g2IcQrBD93P05dPEVcYFzHBwgEXUSYigWUlJTw2mvmQwo0NKg5fnwwc+f+aExbsGABjY2NDBkyBG9v756S2SF3j72b1dtXs2b2GntLEfRi+rSp5OXlcfToUeNTHgPN+1B0Oh1lZWVtllFWFkBExDmTtC+/1MeLNXer1JL4+HiWLVvWFfmdxtPVk+Fhw9mfv5/k/sk9Uqeg79GnTcXNzQ1fX1/jiFWA+vp6CgoKWuUdNGgQjY2N5ObmmqQHBZVw4MBok7QZM2YQEhLiUH0qBpaNWMZDPz7E2MixZsfMCATdpU+bSnh4OOHh4SZphYWF5OfnA5fm6BgmCZoLsqRS6fD0rKWqyhNv7xoAtm7dyrRp0xzSVNQqNfMHz+fr419ztXy1veUIeiF92lTMER4ezk033cS5c+eorKzkiy++ADDOuTHHyJEZZGaOIHnczwCoULF9+3a2b98OwLRp00zytxcBrvk+b29vkpOTrd6imD1oNqu+X8W8wfNwUYmvgMC6WP0bJctyCmDoCbwZeB0oB+4EPID1bW0rilJubT1doaqqijfffNPi/H5+FVRWepNHHvXUM5jBJvszMjIYNmyYcbulSbQ1g9mW3HbZbbxz8B2Wj1neI/UJ+g62uEwNBuYCc4CDwGT05nEzEAM8DLi3sb3WBno6jYeHB7///e9N0pr/2IuKinjvvfdM9qtUOiI0UexV72plKiUlJezaZbo0aXBwsEXD7zUaDRMmTOjsKXTIyPCRfHj4QyrrK/Fx694cI4GgOVY3FUVR3geQZXkskK4oyhFZlt2AG4AwRVEy29l2CN5++23OnTvXccZmhIRcoKwkDE2oxqL8xcXFrdJuvPHGVkG0fX19O6WjM9yXch+vpL7Co5MetVkdgr6HTW6oZVkeDAwAcpqS1IAPUNvOtl0GdGg0Gp5//nlqa2s7ztwOYWGFFBaG4xLqQj31uOFmsn/u3Lm4uLggSRIhISGEhoa2KsPDw6NHn8gM8BuARqehoKKACN+IHqtX0LuRrL3KnizLnsAHwG3ArcAu9Lc304FKYHdb24qiPNtWuenp6TovLy+rau0OpaWlbNu2zRiKUqNRsWvXJGKmrUeHjjj0o1aTkpKIjY3t1oRAW1JRX8Hbyts8OOJBamtrzS434qg4m15wTs3V1dUkJSVZfLWzRUvlcSAUfWftLmAl+lbJ74A64JV2ttvFkR7RZmdnm8S2ra11x8OjjgEMoCauhqVjlzqU3vZIq0+DEPC44OE0msHx5ldZgjNqTk9P71R+W/Sp/An4U7Ok91pkuamDbYcjNze31czklhQVhREWVoQrrriecuXjUx+b7L/++usZNGhQm+v12JPlY5bz+ObHWR4tngQJuo8YpGABMTExJqEpT5w4wffff2/S2VpcHEJ8/K9tlvHhhx+2SpsxYwaTJ0+2rtgu4KZ2Y/yA8aSdT3O6q6jA8ejTppKfn89bb71llbLi4k5x/LhMSsp+k3QfHx9UKhVardYYCT84OJgbb7zRoYbJLxm6hJs/vJkbdY6lS+B89GlT6d+/f6vg2GVlZbz77rut1kHuiKCgi5SV+dHQ4IKra6Mx3bD8R2BgICtWrHDYH6wkSVwx4Ao+OfoJS4cttbccgRPTZ03lv//9r9nlNLqDv38Z1dVe+PvrBwZPnToVb29vIiMjcXFxoaioCJ1OhyRJhIWFOZzBJIUmsS5vHQvkBbi7dD0urqBv02dMpaysjPT0dOMPuV+/fvTr189k0mBHfyVJ4uTJk5w8edJsHaWlgUZDKaSQLTu24GLjt3jhwoWMGmW9VQiXj1nOW7+8xX0p91mtTEHfos+YCoBarW61HIYh5IFOp+PixYvk5eWZGEl9fT2NjY3GbcA4ziAgIMA48raqyhNf3wpjXaWUUk89UZiP+hYeHm6M1G/oa1Gr1cb/7u7ueHp6IkkSI0aMoH///tZ+O8ySEJLA+xnvU1pbSoBHQI/UKehd9GpT+eyzz8jMzLRZ+c2H8nt51XDxYiAajQq1Wks00fzCL22aSmFhocX1GBaHHzBggDFNkiRcXV25/vrrrb5g2f0p9/NK6is8OeVJq5Yr6Bv0alNZvHgxixcvRqfT8d1331FefmkSdHV1tUVrJluKJEFKSip7945n0qQ9eOJJLV0b+h8REcGcOXOMrSlXV1ezS6ZKkmSTcS/hPuG4ql05XXaaaP9oq5cv6N30alMxIEkSV155pUmaTqczWUHQcBvSMs9LL71kcT0hIcWo1VqKi4MIDi5BLalB0sdXaSvIkzkKCgpMAmu35JFHHsHT09NiXV3h3uR7eXbns/x11l9tWo+g99EnTMUckiQRGBjYYb7mj5wNxqDVannhhRfMTkJMTt7P5s0zmTv3R8J14eTr8ulP6/6QpKQkk07jpKSkbpyN9fFx8yE+KJ5D5w4xqp/1OoIFvZ8+ayqdpaKigo0bNwJ6czH0b5SXl1NUVGTM5+KiIT7+JIcOjWDIqBrKMY07FRkZyciRI03S/Pz8TG5pHIVbRt/Cwz8+zD/m/sPeUgROhDAVC/H09GTs2LHGbUmS+P77780OkouPP0la2hgKTiUSF5djsu/s2bOcPXu2U/U+9NBDdpkz5KJyYWbcTH48+SNzBs3p8foFzkmvN5UNGzZw7NixNsefGGgrpGPz1zqdjurqaovqTUr6he+/v4KYmFzU6q6Hl6ipqeGZZ54xu0+SJJYuXWpyW2Z4PXz4cKuEW5g/eD6rfljFzNiZqFWONxlS4Hj0elO55ppruOqqq1qt5WN43RKdTkdtbS1r13YvsqUkgbt7HVqtGrW6seMDuoBOp+Ojjz4iNDSUkJAQ43gXSZJISEjAzc2t40I6QJIklo1YxgeHP+CmUQ4/oVzgAPR6U+nKGA5vb29jB21xcTG7d+822xrQ6XQUFBRQWlraqoxjx2QCA0tN5gEZGDZsGHFxcYwZM6bzJ2QHkvsn8+nRT6lpqMHT1bZPnQTOT683le4SHBzMggULWqXX1dVRVFTEvn37WpmKogyhosKX5OS0VscFBAQQGxtr0ZMnR+KusXfxevrrrBy/0t5SBA6OMJV2KCgo4I033ujUMRUVPhQURDBt2g5j2kEOMopRSEiUlpby9ddf4+vry+233w5c6rdxdXXFkUJmNicuMI6y2jIuVF8gxCvE3nIEDowwlXaIiIjgqaeeYuPGjZSXl3PmzBmTEJLmaGhwwcuryiRNhYrznCeMMGNaRUUF69atM+lMjYuLM9sqchTuT7mftalrWT1tdceZBX0WYSrAxx9/3GbfSGcJCiqlrs6Ds2cjiIzUr8kcQwzHOW5iKgAJCQkkJibi6elJRITjR7MP9gomwCOAEyUniA+Kt7ccgYMiTAWYN2+eSQvEsHbyyy+/3KXyJk3azZYtM3FxaSAs7AL++LcaBAeQlpZGcXExwcHBzJ8/v8v6e5K7x97Nk1uf5Pk5z9tbisBBEaZC2wt2tYwK15Ly8nL+8Y/Wo00lCWbM2MqmTbNJTk4lKEjfAtKhQ8J0fEx2dnarEbaOjLuLOxMGTGBb9jamx063txyBAyJMpRv4+fm1Mp4LFy5QU1NDdnY2v/1tGLfdNp65c78nWAqmhBKCCWbmzJlERESQm5uLVqvlwoULbNmyhdjYWOLi4ux0NpazKHERq35YxZSYKWJAnKAVwlSsyNmzZ9mxYwfHjx83pg0aFE9OTgzDYhuMrZQtW7aYPX737t2dqm/evHmkpKR0XXAXkSSJm0fdzLuH3uW2y27r8foFjo0wFStQWFjIa6+9ZnZfRMQ5jh1LIDY216KyvLy8jPN8Oppc2Dw+TE8zJmIMnx79lIq6Cnzdbbfes8D5EKbSDtu2bWPnzp3dKqOwMJzQ0KI297cMcVlbW8sDDzyAv79/t+rtCe5PuZ+XU1/miclP2FuKwIFwzAV+7czp06d5+umnzRpKZ0MTFBaG0a/fuTb3azQatFqtcfi/Vqtl/fr1Fgd0sieRvpG4qd3IKc2xtxSBAyFaKmaIjo5u88lPXl4eJ0+eZMeOHWb3t8TPr4LS0kD69TMfk3b06NGkpKQ4xTgVc9yXfB+rt69mzew19pYicBCEqXQSnU5HQ0MD48aNo7GxkYaGBoqKiqioqMDPz4+CggKT/AkJWWzZMpNHHy1CkqBeV0/KZSnExMTY6Qysi6erJ2MixvDT6Z+YGD3R3nIEDkC7piLLciKwDKgDJEANvKsoivmFb/oALi4u7NmzB1dXV9RqNSqVyhhjpaqqqlV+V1cNgwef4L33vBg5MpMtbGHmoZkW1TVv3jyioqKM8XMlSUKj0VBUVGRyGxYVFYW7+6XFvyyJDWMunozBJNu7xVOr1a32/2bYb1j5/UomRE1AJYk76r6O1Na9uyzLS4Ao4J+Komib0lyAh4F0RVF+7DGVQHp6us5R4rg2Njby3HPPUV9fb/Exe/ZMICrqNOeivmAIQ/DH8o7Y8PBwY6iF0tJSs/OPDIGwzUXdb+91821L+nEeeOAB/Pz8WqX/fOZnTpSc4MaRN3ZYhrXIyspyugXlnVFzeno6SUlJFncmttdS+UZRlBoAWZbHA0Oa0vOB7j0ScVL+85//tFqdsHkkOcPCZOaYMGEvW7bMJMY7idygTEZiOop20aJFDB06tN3o/ramO1/48QPG83nW51TVV+Ht5m1lZQJnok1TMRhKE48BG5rt69qCNk5A85CRLa/cCxcuNMln+Gt4anP06FHOnz+PTqcjIyPD5FhJgunTt7F58xW4XLGtVb0bNmxgw4YNrdKbm4tWq2X+/PkmsXIdiftS7uOV/a/wyMRH7C1FYEcs7agNBu4FagAd8K7NFPUwp06dMq4WqNPpqKur6/bYlLbQaNS4uGhapfv4+HD99dcbt9vrB5EkybjQu4Hmr93d3e0WACraPxqtTsuZ8jMM8BvQ8QGCXomlppIFvABdXHLPgamqqqKsrAy49AMeP368yfbevXutUldqagpjx6aRSyj11OOGPoZsZWUlb775ZofHu7q6EhBwaX1jc0G8IyMjueaaa6yityvcn3I/z+x8hr/N+pvdNAjsi6Wm4go8iv4JkA643WaKepgRI0YwYsSIdvPMmaNfnqK2tpa0tLQ25+50RG2tB76+5SwbsQydTsfJkyepqanp+MAmGhoaOH/+fLt5CgsLOXDgQKt0FxcXHnnkEVxdXTutuzP4uPkwNHQoqfmppPTv+XlJAvvTpqnIsnwTcEJRlL2KotzWLH2mLMsLFUX5okcU9jA6nY7XXnvNZIEwazFsWCb79o1HpfrZmKZSqfD09ESlUqFWq03W9zG3rMjChQsdfqDcjSNvZNX3q0iOTHaoxdEEPUN7HbXvybJ8myzLV6F/4iMBkcBhRVH+21MCbYlOp+P111839qnYmoiIQo4eHUpj46W+Fa1Wa3Z8izmuuuoqQkIcPz6sSlKxZOgSPjn6CUuHLbW3HEEP0+7tj6Io/waQZTmoabv1cnxOyPHjx/nwww97vN7ycl/q691Qq1t31lrC119/zddff92pYzoKNGUrJsdM5uEfH+bqIVeLZT36GBb1qfQWMzEwZMgQsz+2iooKPv/881bp5kafVldXtxqS3x4lJYHs2zeOWbM205N3BE8//TSBgYGsWLGix29F7k2+l1f3v8pDlz/Uo/UK7ItFpiLLcjwwnqZZzYqiONUj5YaGBj744AMaG9tfKdDSH50kSZ0yFNA/+Zkz50fU6ksD5M5whgFY99GrSqUyWYUR4OLFi/z5z39ulffaa6+1aSjL2MBY6jR15Jfn09+vv83qETgWlj79eQv4ECd9pOzi4sLMmTPNDkM/f/48GzdutLkGd/daE0MBUFC6ZCpJSUnMnj3buN18blDL//bmgXEP8Jcdf+Hvs/9ubymCHsJSU/kJeNMwB8henD171vi65cLpzWm5/cEHH1i8sLqt0GpbT7TzwIMaavCkdZ9DbGxsm4vKl5WV8dlnn9GvXz9mzJhhQ9Xdx8fNh5HhI9mTt4fLoy63txxBD9DRLOVt6MelhAE3ybJ8AkBRFLt8kw2dlOZaHOb6PQwEBAQYB421d/XOz8+3hsxWFBWFEhzculsqmmjyyGOIcVrVJbKzs42v1Wo1M2demtlsOAcfHx9KSkrMnpOlac3Tq6urqa2txcPDo4Mz6hw3jLiBB79/kPEDxotZzH2Ajp7+ONQaDL/73e/sLYEzZ86wbt26Th2TmxvDoEEnWqVHEskudpk1leYEBASQlnZpXeb2Zhx3Jq1lemNjI+fOnbP6GkSGQNnrD64XgbL7AJZ21P5GUZSPml6/qijKvbaV5bgMGDCA1atXk5ubyzvvvGPRMSNGZLBlyyxSUvYRGlpsTFejRkPHj5eLi/XH3HPPPYSFhXWQu+vYclr+2MixbMjaQFltGf4ejh9/V9B1OjQVWZanAlfLslzQlN/xF6axIeXl5ezatQtJkkhOTgZg//797R7j4VHP3Lnf8fPP48nJaSApKQ3DEsrNWykzZ85k0qRJNtNubx4c9yD/3PdP/jT1T/aWIrAhlrRUpgExwHRAAzxlS0GOjouLC6GhocZtnU7H9OnTTQI2SZJEYWEhv/76qzFNrdYxceJe8vMj2bJlFrNnbwbg1cdfxc3NredOwI6E+4QT6BHIsQvHSAhJsLccgY3o0FQURXlalmWltwzN7y5eXl4WLeB15MgRE1Mx0L//WY4dk43bf/3rX032h4SEcO+99zrE42BbcNfYu3h006O8eMWLvfYc+zqWPlJOl2X5eUALvKYoyikbanJqjh8/TkNDAwUFBUyYMKHTYRMuXLhgdqBaSwYPHmyccTxw4EDjrZij46Z2Y86gOXzz6zdcNeQqe8sR2ABLTeUvwN8Bd+BfgPg2AO+++67Jo19LkaTur+kzatQoQH/7Za+gTF1l3uB5rPp+FXMGzcFN3Tdu/foSlppKgaIoB2RZdgXO2FKQI3P69GnKysr46aefujWzWaczP1Zj4sSJREZGmqSpVK3zDhw40OpjSXqau8bexetpr7Ni3Ap7SxFYGUtNZbosy+nol+jQyLK81V4D4OzJmTNn2LRpU7fLcXWtp7bWDQ+PepYsWUJk/0g2nNrArDGzrKDSOUgISeC/mf/lfNV5Qr1DOz5A4DRYairJ6B8lVwPnFEVpvUZEH+Dyyy8nIiKCPXv20NjYSFVVVat5Noa/KpWqzRG6CQkKx44lMHp0Bp9++ikAW9jCmY1nGDZsGEuWLOmZE7IzD457kBf3vshfZvzF3lIEVsRSU3mnKe+/gWeBm20lyNHx9vYmPDycrKwsSkq6FhEiLOw8mZnDjNvjxo2jvL6cx658DHcX93aO7F0EegYywG8AGYUZjAy33WxpQc9iqamUAmWKonwvy/INthTkSNgytGTzztp9+/ZRSCH3HLiHGLq/HGpkZKRx0ffmC7/PnTuX2NhYYz7DEH1Jkmweu7Yt7hhzBw/98BAvzX1JPGLuJVhqKsXAXFmWwwDr/8IcFI1Gg6enJ/3762OBGH6E58+fN7tKYGfQ6Ux/QFFEkUaaVUyl+Wzu5rQV7a5///4sX7682/V2BReVC1cNuYovlS9ZmLCw4wMEDk9Hs5RHAX8D3IByIBp4vgd0WRWdTscPP/xAXV2dxccYrppBQUGt9nl5eXHiROsJgp3BxaWRykpvfHz08Wk98KAOU32DBw82q6ktOlq2tK39zWdA24PZg2bz+x9+z7z4eX3q9q+30lFL5c/ALYqiFAHIshwC/B9wna2FWZsBAwa0al1YsnawIV92djZHjhyxmp6UlFR27pzCnDmbjWmTmaz/O3kyiYmJaDQaPDw8kCQJFxcX/P1770S8e8bew2tpr/Hg+AftLUXQTToyFQ3gIctydLM0pwuIIUkSw4cPN0lraGhg/fr1xhCTzZcx7WhtHWvg4VFP//75nDo1kLi4HADj4mK7du1i165dFpe1fPlyk6U9WtJRC8fX1xcvLy+L67MFg4MHU5ZZRmFlIeE+4XbVIugeHZlKGfB0i7RyG2npUVxcXIxxQ9qKsGYu7fz583z88cdW0TB0aBabNs0ymkpXaS8cpiWtsZSUFJKSkrqlwRo8MO4BnvvpOZ6d+ay9pQi6QUdBmnptRB1Jkrq0KFdQUBArV65s03h27NhBZmYmWq3W2HpoaxVCSQJf30oqK73w8bEs3GVUVBQuLpc+Ni8vLxYvXtwrnpwEeAQQGxjLwXMHGd1vtL3lCLqIpU9/BE2oVKp2+zbmz5/P/PnzOwx4tH79enJycoiLO0V2diwjRrTfX/PII48AeuNy9iH67XHb6Nv4/Q+/F4+YnRhhKjYiNzfXotuksLCipoFwelPRNP0z9K8YWLNmTZtleHt7ExUVZdwODAw0rv/sbKhVahYkLGBD1gYWD11sbzmCLiBMpQMyMjJaTR5sL95rRUVFp54SSRKoVJfKq6WWDDKYwIRWeQMDA/Hx8TGu69O/f3+io6NbTTp09pbMjNgZPPTDQ1w5+EqxuqETIkylA9zc3PD01H+xO4pGr9Pp8PLyQpIkysrKyMvLs6iO5gPhvPGmBvN9MBcvXuTixYvG7fz8fOLj41uNZ+kN3J9yP6/sf4U/XP4He0sRdBJhKh2QkKAPe7h//36+/fZbq5ff2KhCpWq9nJIWLap2nt5PmzYN0I+etXS1xJEjRxqXKnF0YgNj0Wg15JbmEhPQ/VHGgp5DmIqFJCcntxldLTU11dgqqa2t7dRo27NnI4mMNB1WH0YYRRTRj35tHrd9+3aL6zBw4sQJbr/99k4fZy9WjFvB6m2reW7Oc/aWIugEwlSsQExMDN99912Xji0pCSIqyvQ2KYYYjnLUrKlceeWV6HQ6hg8fbvcBa7bGy9WLcQPGsTV7KzNi+1z4HqdFmIoVCA8PZ/Xq1cbtmpoa0tPTiY2NpbGxkXfeeQe1Wm2cLdwcb+9qqqu9CQ6+1Ffii2+rpz8GDLdgnTGxESNGsGjRos6cksOwOHExK79fyZSYKbioxNfVGRCfkhXQ6XQUFxcbR7bqdDqqq6tRFMX4lMbwxKa0tNRkMJy3dyVlZa37OS7jMkAfa2Xu3Lk9cyIOiCRJ3Jl0J2+kv8G9yX12DTunQpiKhbz11ls2WWu5psYTT0/zo2mXLFnCsGHDzO7rSwwPG84nRz4RoSedBGEqFtLZeCPNR9SWlZWRnZ2NJEmUlJSwc+dOY77ycj8GDsw1OXbWrFlcfvnl3Rfdi1g5fiXP73lezAtyAoSp2JiGhgbUajXx8fGUlpaSlZXVYr8rLi6mIRk2b97M5s2bTdLuu+8+QkJCbK7XUQn0DGRQ0CDSzqYxNnKsveUI2kGYio05dOgQ27Zto76+3hhmoTlRUXmcOjWIK68sZu7cuYSGhuLuLgIVmeOWUbew6odVjIkYg0pyuggcfQZhKjZm7NixjB07Fq1Wa+y4/eSTT4z7IyIKycpKJD29kPz8dSbH3vTATXi4eBDpG9my2D6JWqXm+uHX837G+9w06iZ7yxG0gTCVHkKlUpGRkcGxY8da7Zs2bTt7915Obm4MSUlpuLjoHzv/81//JIccxjCmzXJXrVqFn5+fzXQ7GhOiJvCl8iXldb0irE+vRJhKB+h0OtauXdvl5TgsQaWCiRP3cP58MNu3T8PNrZ4hQ44TFq7jonSg3WP/8Y9/mE2Pi4vjppt659X8wXEP8s+f/8mSsL6xPpKzIUylAyRJ4v777+fIkSNoNBrjeBNz/5tPODx37hw+Pj4cPXqUCxcuWFRXaGgxs2ZtpbbWnV9/jefIkeHkSLFkhPgSHZ1LQED7V+e7776b0NBQk8XNeiMRvhH4e/hzqvwUibQds0ZgH4SpWIC5GLftsW7dOs6c6fqS0x4edU1Bm46wX5eGz4XLOXFiOOXl/uh0EBhYSmxsNoGBpSbHvfbaaxaV7+LiYhK+oXl83pZpAPHx8SxbtqzL52ML7h57N8s/Ws6VKVf2agN1RoSpWJnS0lKrPr2JlQaSH7qfsaH1xrSLFwPIzo7ll1/0I3EjIwsYPPi4sS+mI8w9hWqPEydO8PTTLUMVw4QJE+wWDMpN7cbkiMlsPL6Ra+Rr7KJBYB5hKlbGzc2NIUOGoFKpCAoKYt++fd0qL5hg6qk3SQsMLCUwUN/XotPBmTP92bp1BtHReSQkKN2qrzPs3buXvXv3dvq45tMO/P39jeElOsukfpN4K+ct5gyag4eLcwem6k0IU7EyXl5epKSk4OvrS2JiIufOnTNZNlWn06HRaCxe4VBCIpK2HylLEkRF5TNgQD7ffTe3R02lI1JSUox9Tc1vp0pKSszecnWF+1Lu45XUV3jo8oe6VY7AeghTsQJ1dXV8+eWXJmlFRUXGGcXNo99XVFR0u77m5UmSRFFRCPv3j+ayy47h5ubWKtJ/bW0tbm5uBAQEtBtZr8wAACAASURBVPlj1ul01NfXW0WfAS8vL6ZOnWq18swRFxhHg7aBvLI8ovyjOj5AYHOsbiqyLPsD/0S/VOprwOPo1wq6E/AA1re1rSiKUw4+cHV1bdWRu3PnToqLi21SX0NDIxcuhHDyZBwVFb4EB5cwY8aPuLhoqK83f0x9fb1NFppvD0vDaXaXFSkreHrH06yZ3XZwcEHPYYuWSjLwR8AfyABGoDePm4EY4GHAvY3ttTbQY3NUKhVDhw41SZMkicTERE6dOkVxcTEHDhywOOyjOUpL/cjJiaW4WL+28+zZUcyYUUBgYB0ajRcwniFDhhAdHd1+Qe3Q0bIijoq3mzdjI8eyI2cHUwfatmUk6Birm4qiKJtlWY4BHgM+VxTliCzLbsANQJiiKJntbPcaTp8+zZkzZzh48CDV1ZYtFAb6gE8jR47kwgVXcnJiKSrSTyJMTIR774X4eH0/ip5Y6wt3Uq4beh0rv1/JxOiJIpiTnbHVu1+MfiEbQ0BUNeAD1Laz7d1RoS1n+DoyWq2WqqoqBg0aZEwz9GPU1dUREBCARqPhwoULnD59uukY+Oab/nz9tR8hIbUMH56Kh8cBVCqJwiANhZXTaTwWbDPNtbW1TvUet9Q7PWA6f/nmL/w2/rd2VNU+zvYedwVb9KksBrYDHwL/I8vyGPS3P6eBSlmWk9Df7rTc7vAG3Nma5i31Gkbe7tu3jx9//LFV/iNHhhMUVMLUqeeNT4mqqhrRaGDUgFGccz/HpMRJNtPrbLc/LfUmksiBbQcIjQklxMsxw0Q423sMkJ6e3qn8tmipHAPeBWqAscAj6FslvwPqgFfa2XZa3n33XbKzs7tVRnR0Lnv3TiAmZhNqtY6EhAT8/PxITk5mYOxAnt7xNEuGivku7fHg+Ad5ce+LPDPjGXtL6bPYok/lCDC/WVLLWW0dbTslN998s8l2e1ekmpoaUlNTqaioMLkK+PtXMHZsOps3z2LWrM0cO3YMWZZJT09nw4YN7KjewTu57wCXbqUkSSIlJaVVR3FfJcgziIEBA/ml4BfGRLQ9u1tgO0SPVg+yfft2duzY0W6ekJBiLrvsIKmpKUyYsA9FuTSYzQMPsnKz8MJ0aY7c3FxGjhzJkCFDRExb9Iu8r/phFZf1u0zMC7IDwlRszP79+9m2bZtx5rCPjw8qlQqVSmVMaxlWISzsPBkZI1qVFUMMpzlNAqbD2iMiIoiNjSUoKMim5+IsqFVqfjPsN3yY+SE3jOhVDxWdAmEqVqSqqsp4W1JTU0NFRQXx8fFER0cb1/zRarXG1xqNxvjX8H/btm0m6yU3J4IIswuMTZ482ek6/2zNxOiJbNy8kauHXI2vu6+95fQphKl0kcbGRhoaGowtjsbGRl599VXqm4a0toyv0pmZwVlZMkFBrY2lrbWVP/74Y+PrwMBAVqxYIZr9wAPjHuBf+/7Fk1OetLeUPoUwlS6Snp7Otm3bjI+JDahU+h++VqtFrVYb0w2tlPbQaCR2755EcHAxY8a0H/GtLS5evMif//xn3NzcCAwMNKa3N3HP29u7V0aJi/SNxMfNh+PFxxkSPMTecvoMwlQ6iU6n4/Tp05SVlVFXV9du3s60Turq3Ni8eSbjxqUSEtK9OUOenp6MHj2auLg4k3RzxqLT6VCr1ZSVlVFVVUVpaakxrTfEvr0n+R4e3fQoL17xomi99RDCVFqwZ88eNm3a1OP1Hjw4mokTf+owZKQl1NTUdDnWSXOio6O57bbbuq3Hnrip3bgi/gq+Pv41V8tX21tOn0CYSgsuv/xyi1cHtOQRsaWEhRVSUBBhkalUU40XXmYfH5u7GltyhVar1QwaNKhTYTOdhbnxc3noh4eYPWi2CObUAwhT6SR5eXm8/fbbVi83NjaXH36YQ3z8CVxdNe3mTSON8YznyJEjrfYtXry4zf4TczFUDIYjSZKxP6g3IoI59RzCVDpJVFQUq1evNruvvLycqqoq3njjjS6VnZycyv79KfzmN2eMI2SbtzIMncKNRxupbahlbPhYk0fSDQ0NuLu7M3jw4C7V35snusUFxtGobeR02Wmi/bseHkLQMcJUrIifnx9+fn6sXr26zWH6jY2NaLVaJEkiIyODkpIS9uzZA0BQUCkajZrMzFLy8n5osx4dOtJJx6fQxyR92bJlREREWPekehErxq1g9bbVPDfnOXtL6dUIU7EhpaWl/Oc//zG2MAwD3wytj7KyslbHpKTsY+fOqcyevbnVPgOuuBITF8PvF/4eX18xsMtSvFy9GD9gPFtObWFm3Ex7y+m19N6baDuTn5/P5s2bcXNzw8PDAy8vL/z8/AgICMDHx8doKM3HkgC4uzfg41NJZaWXuWKJiIggPDycorNFpB1M63bg6L7GosRFbDy+kQaNZYHHBZ1HmIqN8PX1ZciQIYwfPx4/Pz/Onj1LXl4eeXl5JmElzQ3J79evgPPnw8yWW1BQQGFhIZf5X8ah7EOdXsOnryNJEncl3cXr6a/bW0qvRZiKjfDz82PEiBFs376dY8eOmUSx72hkbXBwCUVFoe3mWZC0gMXTF+Pq6mo1zX2FxNBESmpKOFd5zt5SeiWiT8WK5OTksH79+m6X4+dXQWWlLyUl/9/emcdFdd77/z0zwLAMqyibbCIeBUQRFVzivmVvmusva1ObtWkSl2hM25u2N8mtadObNkuTxmymSdq0SdvsRqPGLZGowRWFAyqLgCuCCMKwzPz+GGZkmWEGODAzzPPOKy/mbM/5zjjnM8/yXUIIC+tY2jQ9PR2VSsWZM2cIDg7u8708leXZy3lm5zM8M+8ZZ5sy6BCioiAJCQmW5WZrqz+5ubl8/vnnDrU1a9Y2Nm2az6RJexky5EpqhEOHDlleHzlyhFWrVnWIMRI4RpA2iNRhqew6uYupsY45OwocQ4iKg3zxxRd8//33A3Y/jcbA/Pmb+Prr2aSkHCUmpmt5j8bGRtasWcO1117LhAkiy1lPuWPsHSzbsIysmCw0aiHMSqFyl9WD3NxcY2ZmprPNcJjeJDg+e/YstbW1lohm0/9GXnpJR0PDdpKS7OfAffDBBxk2zPokb3/Y7EyUsHf/qf3knsrl3gn3KmRV97jbZwymHnZmZqbD0Ziip+IiGAwG8vLyKC0ttZTsMBMdDbt2TcHLq4X4+I5FB0ooIYEEy3ZjYyMCx8mIyuA/+f+h6nIVQ/z7r/yJJyFExQU4cuQI//rXv7o9Z8qUHDZvntdFVNQj1MxNmovWW4tarebcuXOcO3fOctzX11fkrbXD8uzlvLD7BZ6a/ZSzTRkUCFFxAVJTU0lNTcVoNJKTk4Ner7cE+JWXl1NUVERrqwajsWsP9PKJy2yo28DkmMmAyQ+j/ZDW39+flJQUkUukG4b4DyE2KJb9p/aTEZXhbHPcHiEqLoBer2f79u3o9Xr27dvX5XhTkzdbtswhO3t3h/3BwcGkk86euj1EnYwiICAAb29vi4DMnTuXiIiIAXkP7s7dGXezYuMKxkeOFwLcR4SoDCClpaXU1dVZts3Z9Jubmy3zKNHR0R1ihYqKGvjmm6uYMWM7gYH1Hdozu/pXUcX5y+ctdZjB5GQ3GDK3DRQatYZb027l74f/zh3pdzjbHLdGiMoAUlxcbJnvOHr0qN3zKyqiOHIkmwULNnbIsZKSksLixYst261bW3li5hNiWbSPTI2dyqfyp9TqawnSCkHuLcJNfwCZNWsWixcvZvHixXbd6w0GOHQonfnzN3VJ2nT06FHefvtty/ad6XfS1NrUHyZ7HMuzl/PCdy842wy3xqN7KkeOHOGLL76wbHee5OxM+7G2eehi7TWY8qaYc922P15V5VhSa4NBg59fI7aG96WlpTz55JOWbVuJowQ9I1IXSYhvCAXnCxgdPtr+BYIueLSoSJJEQkJCFzGxNVHXPijQ/NfW62PHjpGUlERzczOvvvpqj23z8mrFYFBRVxeATldv87xp06Yxb968HrcvsM1PJ/6UxzY9xp8W/klM2vYCjxYVLy8vvLz65yPQ6XSWMqSdexF6vZ7f//73dnOhTJv2LVu2zGXRog1YSx/7q1/9alDnlXUW3hpvbpBu4D/5/+HmlJudbY7b4dGiojRlZWWsX78eo9GIXq9n27ZtXXoy5m1fX18aGhq6bU+rbWbChH3s2DGDmTN3dBkKPf3005bXMTExNtuZPn060dHRgKkX5uXlhZ+fX2/eoscwJ3EOqzetZkHSAlE2tYcIUXGQ48eP250PaWpqIjo6mtbWVmpqajqkJuiuG200Gjl8+LDVY5GRZ9HrtWzZMpeZM7fj7W09KVNFRYXN9v/5z392a7cZtVpNcnKyxdbY2FiHy5UMRlZkr+D5757nVzN/5WxT3AohKg5SU1PTwf3dFhqNhv37e1ey1Bbx8ScJCalh06Z5pKYe6eKqb2Y3u8kiy+F2x4wZQ1hYGP7+/uh0Os6cOcPw4cMtxz09X0tUYBRD/IeQdzaPtGGDrx5SfyFExUF6EiF97bXXcvToUbZu3cr58+cVuX9w8CWuvnoDhw6N5dixZKZO3YWfX8fgwTrqMGJEhWOTixUVFZw7d46xY8eSnp7ulhG0/c39mfezcuNKnl/0vJi0dRAhKv2ESqXioYceAqCoqIg9e/ZY5lTOnj3bwbPW8TZh3LjD1Nf78e230wgLqyI9/TBeXiY/liEMoYoqwgkHTN659913n3JvygPxUntx29jbePfQu9w17i5nm+MWCFEZAJKTkzsU+Dp69Cgffvhhr9sLCGhg3rwtnD49jO3bZ6BSQXx8KcPj6inxzicjNsPiG/POO++gVqvJysoiLCzMst/8q9t+Arm+3vbStSeTPTybzws/50LDBcL8wpxtjssjRMUJpKSkOOSsZn7YW1tbMRqNrFu3jtOnryRrjow8S2TkWVpb1ZSVxXE0ZxFFLanUa+KIiakgLq4MHx9TKYrjx493e6/Q0FC8vLzIzMwU3XwrrMhewZ9y/sTTc562f7KHI0TFhTH3KMy+KA888IDlWPvscJs2bcLX9yCJiSW08DXTWuZTXh7D7t1ZNDV54+3dQkJCCbGxJ7n99tsYNWqU1fvl5+cLQbHBEP8hJIUlsadijyXNhMA6QlTcFLVazfnz5/nLX/7SYf8sZqH2aiUhoYyEBFPkc1OTF8XFI/jqq/kcO7aXRx6pZvLkyUJAeshd4+5ixYYVZEZliuDNbhCi4iI4kv3NEdRWYkR9fFqQpEIkqZCdO6fxj3/ksGHDBofbnDlzJgDnz5+nqKiow7yMSqWipaWFwMBAvL29OwzPwOSUZx7GDR06lJtuuqlvb9CJqFVq7plwD2/se4MHJj5g/wIPRYiKi2DO/qY07UMBjEYj33xTxaFDE7jxxnp8fHzw8fGxPPTHjh1jxIgRHSZvS0pKyM3NBeh2edzsGLhw4UIiIyMtwzZzb0itVqPVahV/fwNNekQ6Hxd8zOm600TqIp1tjksiRKUf2bt3L+vXr3e2GR04dSqCs2cjqKo6ZP/kTqSmpqJSqaiurqalpYWsrCwCAgLQaDTEx8f3WxyVq7E8ezlrdq7hd/N+52xTXBLP+BY4iUmTJjFp0iQAGhoaOpQ77dyDANPka25uLq2trZb9BoOBhoYGamtrLdsGg4GgoCAiIyMxGAwYjUZ27NjRrS2XLuk4fDgNf/9IfvlL0GqTOiwtA1RWVhITE9NlrqVzygeAUaNGeWwwY5A2iAlRE9hWso1ZCbOcbY7LIURlgHAkgM9oNHLw4EFqa2vtnltRUUF+fr7VYwYMNNCAsW4oBQVjqKkJISCgjrS0PIKCvsNe56m9V6018TMajSQmJg6K4UxvWZyymBUbVzA1dio+Gh9nm+NSCFFxIZqbm7njDuv5US9evEhpaSl5eXmWHob5IVer1ZZew8WLWvbvT+T7c2cZpQ1jzJgCQkNrurQ3fPhwy/XmthoaGvDz83NI1N59990u+9rn1rUlRt3ts3eut7c3DzzwgEsMs1QqFT+b9DNe2fsKy7OXO9scl8L5/zoCC2VlZWzcuLHLfkfih/R6H3bunI6/fwNjxuRxJu3fTGWuzfPLy8sB0zxJVpYpCLG0tJT4+Hig+6jq7jBPznYeWnXOjGdrX3fnmtM2uAqjhozio/yPKK0pJT4k3tnmuAyu8y/kYWzevJm8vDxLRvy+0NKiYfNmU2oEne4yAIEEcolLBNJ9LpDW1laam5vb2mmxvHaE9r2JqKgo/P39e2G9e7M0aylPfP0Ezy18ztmmuAxCVDA9TK2trTZ/Oe1hDhJ88803+8M8u+zYMYOpU3dZBAUgnnhKKSWN7kP2CwoKKCgo6PacxMREm8fMQ6dZs2YRFxfXM8MHAX7efswdMZfPCz/nulHXOdscl8BjRWXjxo2WeBhH8qS4KgcPphMXV0poaMcez1CGcoQjvWozMDCQZcuWdQgRENjmmuRrWLlxJXMS5+Dv7Xm9tc54rKhMmTKFjAxTiUtruWJtTSR2Pm4wGJzWQwGorIzi6qu7zsOoUTOMYQ61MXv2bFQqFefOnWPYsGGoVCq+++67LnMj5tdjxowhMFCkWGzPsuxlvLj7RX4+/efONsXpeKyoBAUFKVLBr7a2lqysLFpaOqZ5rK6uprS01OJz0l+MHi2zadNcpk7dRUBAx5y39oY+ZrZu3dqje3755ZcOnZeWlsbNN3tG4ui44Dj8vPyQz8tI4ZKzzXEqHisq7fn3v/9NXl6es83oFYmJJUREnGb37iz8/BqYODHXkrTJmajVao4ePcqRI12HYPaqCPzkJz9xy/mZByc9yOpNqz2+tIcQFeDmm2/u8y9qS0sLzz77bI9WT5TC37+R2bO3c/58GFu2zEWSZBISSgfcDoAbb7yRyMhIhg4datln7QFzdJ874aPx4abRN/HBkQ+4Je0WZ5vjNISoKIRGoyE4OFixnLS9ITz8AgsXfsXWrTMJDq7pMnnbW8zzLJ0xGo34+/szd+7cDs54zc3NVFZWWi2+1vlvREQEAQEBitjpCsxMmMnjmx7n6uSrPbYesxAVhWifkxZQLIl0aWlph7rJjjB58l4OH04jO3tPn+8PcPbs2W6Pt5+oHjJkCL6+vl2W5239nT59OklJSYrY6SqsmGLKEvebWZ5ZilaISh/pbnXIlgv69u3b7QYA9oWAgMtcvmxa2rzMZWqoIZroPrUZGBhISEiIQ+e279WMHTuWiRMn9une7kakLpIIXQQHTh9gfOR4Z5sz4AhR6cShQ4f46KOPnG1GnzE/11q0yMh9FpVLly4RExPD2LFjbdzvipC0jykKDg7m8uXLVs+zt6+7OCAwFW9zVe6bcB8rNq7g+UXPo1Z5lq+PEJVOpKenEx0dTU5OjkPndw7KM1NTU2P1l72xsdFmdLGSmE3RoMGAofuT26FWqy1zHJ17WuXl5Zw82bWQWXd+Pn3d113MUEBAAOnp6S45watRa7hr3F2s27+Oeybc42xzBhQhKlYIDw/n+uuv71MbtuZUmpubSUxMtETyVldXs2fPHoKCgkhLS2PXrl19uq8ZrbaJxkYffH2b8MEHPXq0aImKirKc0/mBbWhosMTvmCOfNRoNarWaGTNm2KzX7KyH2tUTdU+Mnsin8qecrT/LsADHHBEHA0JUBhhvb29L4iYzV199Nd99953VCOXeEhV1imPHkklLO0IccZRRRjLJnDp1qtvrLly4YHV/UVFRn2165JFHCAvzrLo5K6es5Hff/I5n5j3jbFMGDCEqTmDt2rVdEkQrzYgRxRw8mM6WLXNIn7CH4tDPSSbZ/oU9xMvLi9WrV+Pt7a1424OBYN9gxkWO86gscUJU7FBTU8MLL7zgbDN6xbhxh9DrvTlwIIOLl9KomXiUkBD7CZh6QktLC2vWrFG0TVfA29ub1atXK5K/5ZbUW1i+YTlTY6cqYJnro7LnMu0q5ObmGntSJN3Z2JpTaWxstNRRfvnllwfMnqYmb3bunE5KylGios70673uu+8+i69KaGhov93HnQrKF1UV8VnhZ1wdcrXb2GwmNzeXzMxMhyevRE9lgCkoKGDDhg3o9foBva+PTzMREWdobLSfK9cWSv56exrJQ5JpMbRwsu4kY3AvUekp4tsxgJSVlVFUVMSIESMs+xxxZVfib3W1H01NSSxatB8YbtU+c4pJWzQ3N/Pb3/7Wsr148WKb78MW7Zfe+7Jyo1KpKC8vx2Aw9LkdR4iPj+9zZrulWUu5/4P7mT9xvkuvWvUVISr9SElJCVu2bLFs19bWOpRUuj84f34Ily9X2RWOnvDxxx8DtnPMgvXE1ubjnY9ptdpuo5M7n3/p0iVqarom9e7u2s7+RI4SFhbWZ1Hx9fJlZtRM/pP/H25OGbwpITxeVC5cuMD69eu7/aLZS9Rkjfr6eku94/bnqdVqy4NlzZW/vwgPryIi4iw7d05n2rRvUCKhW2xsbI9+ce2dGxQUxHXXXedwm+40p2JmauRU1pWvY0HSAgK1gzPRlceLik6nIzs7G+j+S2/rmK39ZWVlxMfHd/hl7nxNcXEx27dv743ZvcIveRs+3iPJy0sjPb3v+WNOnDihgFVXeOihhwb1sMDMo1Me5Y85fxy0AYceLyo+Pj6MHDlS8XYbGxtJSEjo9pz4+HhmzZpFa2sre/bswWAwcPLkSWRZVtwegIDLcRw+KrF4WrHdc8PDwwHTPEpPMv7HxsZaemMAc+fOdcuES/2JOeBw/6n9ZERlONscxfF4UXEFjEYjdXV1GI1GwsLCmDJlSpdzbM0FtC/e1X6y1LzvwIEDAJw5M4x9+zKInf0UaaP+C39/k7t+YWGh1XZ7khcmLCyMhQsXAiaRbm9z+2RNgivcN+E+lm9YzguRLwy6gEMhKi6Al5cX8+fP75e2b7zxRgAefRR27YKnd2Rx51V3Wkp1vvfee5aqAp1x17SO7oBGrWHJ+CW8tf8t7p1wr7PNURQhKi5AWVkZH330EbW1tR1KmvY1aXZgYCCjRo1Cp9NRWTmCbdtO4FPtw+ufv05qcKrd+Zx169ZZXgt3fOXJjM7kU/lTztSdIUIX4WxzFEOIygBx6tQpXnvttQG956VLl8jNzQUgIOA4r7wSTPLYWraznfP0LO2lo+74ERER/PSnP+2VvZ7Io1MeZc3ONfx+/u+dbYpiDBpR+dvf/saxY8ecbYbLEhdXzsaNY0hJOUq2JrvL8fDwcCIiIuwmrLZF+3kZaz0ga/M2ISEhJCdbD3JUqVQeUUY12DeYzOhMvi7+mjmJc5xtjiIMGlG54447nG1CBxzxocjPz6e6upqWlpYe197pDRqNKVmTiq5icf78eack7fbzM4UNdPa0VavVPPLIIx0mfgcri1MWs2LjCqbFTkPrpXW2OX1m0IiKO9JedCZPnozRaKSoqKhf0lkOH55BSkoUc+bMsHnOuXPnOvRUWltb+eabbzqcs2jRIuLj44mMjFTcRk9FpVLx8OSHeWnPS6yausrZ5vQZISouQnl5OXv2mLLfJycnK5IUqT15eccpKAhh27bt9MW/bMOGDcoZNcDodDqbx2644QabQ7GBYGTYSIxGI8cuHGNkmPJ+UwOJEBUXYdiwYZas859++qni7YeE1CJJMhs3LiQ+vhRJKlDEVd8VsJa/tvNrnU7H3XffbXOeyDwMcyZLs5by+ObH3b7CoceJSmtrKx999BFNTU2W5dvGxkZKSx2v6Ddx4kTUVp7I9pOR1dXVFBfb91y1RkpKCgB79+7t1fW2iImpJCamkpKSeL74ahYj4k6RmtI/3ruOMG3aNMLDwy1Z8f38/DpMztp7sMrKypgyZQparfvPQwBovbTcIN3Av47+i8Wpi+1f4KJ4nKio1WoyMjIsBdWNRiMXL17k3LlzHbxSDQaD5Xj7/Wq1mtDQUDQaTbf3aWlpYciQIQ7ZVF9fz86dO3v7lnpMQkIp9QnruVQ4g/37x5GRcZChQ4fi7+9PUFCQzZIY1tIqLFiwgODgYKsi29/o9fpBIyhm5iTOYfWm1SwcudBtKxx6nKioVCqrFfGysrIUvU9PImgbGhoGVFQA4ojjyKgvqP76F4BpkrY3FBQU8KMf/cgiSvbEVmCflVNW8tyu53hy9pPONqVXeJyouCJ+fn785jcDF7FaUVHBpUuXeP7Q84Q2jiImJo7s7DrKy8uJioqyxA599tlndtvS6XSW1aof/vCHJCYm9rf5g54IXQTDg4az79Q+JkRNcLY5PUaIigdirt8TXhnOU7+N4IMPWnj22YtMm1ZLfLyPZa7p1ltvZcOGDd0mQoqMjLSkl9y7dy979+61XA90+9pM+/32Mshdc801HhEqcM+Ee1j25TLGRYxDo3av3l+/iYokSR8CfwCeBGqB+wBf4K+2tmVZdk5aNBfGaDTy5ptvcunSJcs2YNnuC/nk88jOhxnGMBIS1HzxxXjefdfA1Kk5BARctnpNSEgIGRkZREREdImM7vzXkded91kTnfbbnjK8UqvU3J1xN2/se4MHJj7gbHN6RL+IiiRJ84BJwCrgMUAL3AXE29n+c3/Y46oUFhby/vvvO+3+iSRSg6kXotEYyMzch17vw+bNc7nmmi8t/iztH/CamhoqKirQ6/UYDAYMBgOtra0YDAby8/MJDg5Gq9ViNBqpqKjAYDAQHBzMgw8+OOgmVfubjKgMPpE/4XTdaSJ17uNsqLioSJLkBSzB1AMZIctyniRJPsDtwDA72x7FqFGjejWXsnv3bkWc0ALb/muPVttEYmIJFRXRDB9eCcCMGVe8cHfu3ElhYSGFhYVW22xsbOTWW2/Fz8/PUo4kKSnJI9zt+4OVU1by9I6neXb+s842xWH6o6fyMPAicDVgLtyrAXRAYzfbAfYaHojC5kqhZCH2Xbt2KZqw2h7JyUVs2LCQ4OCLBAbWs337doKDg5kzZw433nij3XkSg8FAfX09YJqErqyspLKyUnE7B6rYvZL0xuZoYzRvbnuTqRHuUAKPEwAAF1BJREFUUYysP0RlPBAHTAEmSpI0AdPcSRlQJ0lSJqbhTuftk/Yadqckx0omZbbWjtFo5MiRI90mz26//+LFiw4HLXp7tzBv3ma+/34ier0viYnFJCQU90tMUm+YPXs2sbGxnDlzhpCQkA7vMzIy0qWjm3vzvRg9ejQrNq7gjuQ78PXy7SfLbGNOn+EoiouKLMtLACRJ+jVwP7AaU6/kfkAPvNzNtqATOTk5fPXVVwN+Xz8/PVdd9S0Gg4ri4gS2bZsFwIgRJ4iPL0Otdk5ly4SEBIqLiykuLuby5csWT2hzL2nmzJnEx8c7xbb+QqVSsTRrKS/ufpHV01Y72xy79NvqjyzLT7W9/FGnQ/a2Be2YMmWK1Zy19igvL7d4DZeUlHDq1Cmg43DFz8+PgADTqLNzNLIZtdpIUlIxSUnFtLaqOXEika1bZ6NSGUhKOk5c3MleBSj++te/7nN8izuW6OgtI0JHoFFpKKoqInmI8wIfHUH4qQwSjh8/znvvvdera4soYiQjreZZaY9GYyA5+TjJycdpadFw7NhINm2az/z5m3osLE899ZT9k+zg6+uLJElOCRFwBg9PfpjVm1bz/KLnXTrgUIjKICEpKanXXrlv7HuDuYlzaTzdSHx8PM3NzVbPszZ389RTWu6/fwxa7ZWl5bVr1zp8b5VKxcKFC1Gr1RZxaD/xm5KSYnPlKD8/32MEBUwBhz8c80P+kfcPbht7m7PNsYkQFQ/HaDQyI24G20q2MdlnMlqtFh8fH4ur/s6dO2lqauoQWGk6ZuS113zw96/n9deP2L3PypUru81n0htc+de6v5iZMJNfbP4FC0cuJMwvzNnmWEWIyiBg/fr1fUqTYMTIVrZSRplD59fVBbBz53TS0w8TE+PYUvFzzz3Xa/vMZGRkcMMNN/S5HXdn1dRVPLfrOX4797fONsUqQlTckP3793PyZMcV+IyMDKvDE3Mxse6wN5fSnvx8icrKGObM+Rqt1vowqb/Yv38/+/fv77Dv9ttvd2rGNmcwxH8Io8NH823Zt0yLm+Zsc7ogRMVFqaqq4p///KfV+JgLFy4ofr9AArnIRYIJtnq8sVHLzp3TiY8vY+7crxW/v5mEhATLipStvC5jxoxh7Nix5Ofne5ygmLkz/U6WfrmUSTGTLIXhXAUhKi5IQ0MDX331lSVRlHkZ2DwpGR4ernjm+xRSUGN90rO4OJ7CQonp07+xGWhoJjExkdjY2C5Bg0OGDOmwnG2mJ9vm16KU6pVk2S/veZkVU1Y425wOCFFxQby9vRk3blyHfbYmJcvLy9m1a1ef7xlgI0qisDCZS5d0LFjwlUPLxmbHtM5MnTrV0gMxvxedTmfxM7GVW9a87YmTsvaQwiU+kT/hRPUJRoSOcLY5FoSouCBeXl6WPLX2GDNmDNOnT2fLli0cPnwYMPUObC0L9xS9XotOV9enDPyAIsLXmf6qP+1OLM1ayuObHncp3xUhKoMAPz8/srOzCQw0RRwbDAZ27NihSNtjx+Zx9OgYNm5cQFraYWJiTinSbk+56667umSVc7dgwv7A18uXH4z+AR8c+YBb0m5xtjmAEBW3p6KigjfeeKNf75GSko8kFZCXl0Ze3th+EZfx48czZcoUhg0bpmi7nsDsxNn8fPPPWThyISG+Ic42R4iKuxMdHc2ECRPYt29fv95HozEybtxhxo7N4/vvM6mv1zFqVO8Kno0dOxYfHx90Oh1XXXUVYKpS4Crdd3dk5ZSV/N+u/+N/5/yvs00RouLuqFQq5s2bZ5nY7Typ6UheWPP2iRMn+OLiFyydsNTuA/7rX/vz//5fGSdPlpGTk2PXzuuvv96SWzY1NdWj3OsHgqEBQxkZNpKckzlMie15AKqSCFEZBPj5+REXF9fndqqqqggzhqEKUBHuH97tuXPmwOHDQWi1DaSnp3Po0KFuzzdn5vf19bUsL5sxvw4KCnLpXCiuzo/H/ZilXy5lYvREvDXOSw4uRGUQU1JSwl//+lebx631RiqpZOmmpaSoU2htbbV5rcEAGzYsYtGigz3KrdLY2Mjrr79u9diNN97I+PHjHW5L0BGVSsXPJv2Ml/e+zPLs5U6zQ4jKIOHy5cv8+c9/pqGhweFrrLn1D2MYBRQgtUrdXqtWmyodVlTEEBvreKpLlUpFdHS01WPff/89BQUF3HLLLWJ+pZeMGTqGT+RPKK4uJjHUOTWYhKgMEvz8/HjgAVMph6amJgoLC6mqquoSK2OLlJQU9Ho9x48fx4DB7vlGI5SWxrFgwaYux8zzNT/4wQ8IDQ3t0fvw8vISgtJHlmcv5+ebf+60Qu9CVAYJKpWK4OArcTtmV3ZHo3qrq6v59ttvWbBgATlf5dBEEz5ciSkJCwuzFCErLi5mz57hSJKMRnNFgJKSkrjllls8otiXK+Pr5ct1o65zWqF3ISqDiObmZtauXUtVVVWf2plC19WDCxcudAhkNBg0GAwdV3COHz/OmjVrulx79913Exsb2yebBD1j3oh5rN60mgVJCwj2tR4k2l8IURlENDc34+PjY4mxAetVANvvb4/RaLRaAdBaFzo19Qj79mWwbdsMsrJ24+ent2mXqPnjHFZNXcVzOc/x1Oy+p+7sCUJUBhH+/v7cf3/vixI4kki6uLiYd955B5UKMjP3c+lSAN99l42PTzOZmd/j69tkOTcmJoZ777231/YI+sawgGEkhCSwu3w3WcOzBuy+QlQGOQcPHuTjjz/ut/YDA+uZPXs7tbWB5ORMRavVk5mZi1bbREVFBU8++WSv2l26dGmPJ3kFXVkyfgnLvlxGZnQmXuqBedyFqLgZBoOBF154gbq6Oss+W2kDzK/NE6dKRS5bIyjoErNnb+PixSC+/XYavr4NZGbuQ6ttsn+xFV588cU+23TvvfdaJpc9FbVKzYOTHuQve//CI1mPDMg9hai4GWq1mqVLl1q2bc2T2Js/MRgMvP7669TV1WE0Gi01gtpTTTWh9Ky3EBxcy5w5W6mpCSYnJxuDQU1q6hEiIs71qB1HGDp0aJdk2ub3qlKp2Lx5s0PtzJkzZ1BPJKcMTeHjgo8pu1hGXHDfPa/tIUTFDbE2mdob7r77bgwGAy0tLbz00ktdju9nP3OY06u2Q0IuMmvWDpqbvcjLS+XAgfHExFQyenQ+Xl72/WAc4dy5c5w7171YLViwgMjIyG7PCQ/vPiRhMLA8ezlPfP0Ef1z4x36/lxAVD8NgMKDXm1ZqNBoNGo0Gb29vVq9eTWFhIRcuXGDv3r0YjUZ8m31RaVT4qHwsiZ9s1W22hbd3CxkZBwGoqIhm27ZZaLVNZGTsQ6frPjWloyQnJ3P77bcr0tZgxd/bn3kj5vGp/Ck3SP1bkUCIiodx6tQpm5UMW1tbLb0glUpFvCaeClUFIzUjUalUlqVh81yNeejkKDExlcTEVFJf78e+fRNobvZh/Pj9hIXV9Oq9zJo1i5kzZ/bqWk/kmuRrWPXVKuYmziXAx3r6UCVQ9fSXx1nk5uYaMzMznW2Gw7hjnd/ONl9uvsyfcv7Ef8/4726vW7t2LadPn+7x/fR6bw4cyKC2NoiUlKNER1f2OW2lu5KSksLixf3v/VpeW87bB97miRlPOHxNbm4umZmZDv/LiJ6KwCrV1dVs2bKF3NO5fHD6A6t5WVQqlSUvbm/QapvJytpDS4saWR5Nfv5o1GojcXFlJCSU4OVlO0raFmq1miVLlrjsxKuzf2yGBw0n1DeUvLN5pA1L65d7CFERWMXPz4+EhAS0R7TkV/VvLlgvLwOpqUdJTT1Ka6uKsrI4du6cTmurFxERp0lOPoavr22PXTPJycmo1WrKy8sJCTGlVbQVUOdooJ2tnry9Hr75eEBAgGIT60rxwMQHeHTjozy/6HnUKuWTZQlREVjF19eXiRMn8v7E922eYzQa+eSTTzqkWzA/rK2trRiNRsrKynrkH6PRGElMLCUxsRSjEc6ciWDPnkno9VqSk4uIjy+zOkTS6XScOmXKm1tRUcGuXbt6LQi26KlAqVQqbr75ZuLj43t1v/7CS+3FHWPv4J2D77Bk/BLl21e8RYHHYE5v0B16vd4S4FhTU4Msy10KjXVuc+jQoRZvWnPMUnV1LWvX+rJ58zx0ukukph4hKOiKA6Cvr6+lwJpOp7Ncl5WVZSl3YjQaaW1tpaWlBZVK5dGOcVnDs/is8DPOXz5vN8tfTxGiIuhXtFqtJSlTdHS0w/WMrGFe6CktbeL992dQUHCJuLg6kpMP4+vbzMWLF2lubu7gbbx161a2bt3aoZ3o6Gh0Oh233HKLR+fKXTllJX/Y9QfWzO0aWd4XhKgIFOftt9+mtLS03++TmAhnzw7lm29G09TkQ1JSAyNGdK2O2JnKykoAnn76aYfv5ePjw2OPPYaX1+B5ZEL9QkkdmsrO0p1cFX+VYu0Onk9I4DIsWbIEMA036uvruxw3z0HU19fzySefUFVVZXHI6ynDhp1j2LBzGAxQUDCaffsymDDBsWx3AJIkdZkTMRqNGAwGjEYj4eHhLFy4sFe2uQO3j72dZRuWkTU8S7FC70JUBHYxGA2sL1rPdaOu69F1er2eV199FbgyN9KTHLo9QaWC2tpgwsOtF65PS0vrsixuMBi6LImrVCrS0tKIi4tDo9EQFBSEXq9HrVZb/h9M6S5VKhUPTXqIV/a+oliybCEqAruoVWpyTub0WFR8fX1ZtWpVn+69bt06ysrKuj2noUHL9u2zSEvLY/jwCqvn5OXlOXQ/o9HI4cOHbfrf+Pr68vjjjzvUlrsghUt8VPARJTUlJIQk9Lk9ISoCh9CoNbQaWtGoe+dz8e2333LmzJkeX+fv749Op+sw+dqeiopoDh1KZ+bMbfj7N/bKthkzZjB79uxeXTtYWJa1jCe+foLnFj7X57aEqAgcYnzkeA6cPkBmtOOhEpWVlRw+fBij0UhJSQm1tbUd0lp2fg1YTcFgjdZWdVvGOT0LF27AvIgTGhrKTTfdREhIiKVgvcA+ft5+zE6czReFX3DtqGv71JYQFQGAZX7B1lCjpbmF18+/zvmo81b9S3bv3t3fJlqorIzkwIEMJk/eTXj4hQ7HqqureeuttyzbDz/8MKGhoR69dOwo1426jkc3PsqcxDn4efv1uh0hKgLANGHn6+tr9dddpVIRRBD6aj1BQUFWr1+wYAFAl55H578NDQ0O1V62RW1tIEVFySxa9CWO6MRbb73FbbfdxvDhw3t9T09iadZSXtrzEqunre51G0JUBIBJOJKTk7sNdosaE6XIRJ5ZgBylfTY3cydJpZrj9OC8wUhCSAL+3v7U6msJ0lr/AbGHEBWBwyghKL2hY95dp5jgUTw8+eE+XS8GmgKBQFGEqAgEAkURoiIQCBRFiIpAIFAUISoCgUBRhKgIBAJFEaIiEAgURYiKQCBQFCEqAoFAUYSoCAQCRRGiIhAIFEWIikAgUBQhKgKBQFGEqAgEAkURoiIQCBRFiIpAIFAUISoCgUBRhKgIBAJFEaIiEAgURYiKQCBQFCEqAoFAUYSoCAQCRRGiIhAIFEWIikAgUBQhKgKBQFGEqAgEAkURoiIQCBRFiIpAIFAUISoCgUBRhKgIBAJFEaIiEAgURYiKQCBQFCEqAoFAUYSoCAQCRRGiIhAIFEWIikAgUBQhKgKBQFGEqAgEAkURoiIQCBTFS+kGJUnKAn4DlABvA08CtcB9gC/wV1vbsizXKm2PQCAYWBQXFWA88Jgsy0ckSfoD8BigBe4C4u1s/7kf7BEIBANIfwx/MoAVkiR9DIyXZTkPOAwkApF2tgUCgZvTHz2VVbIs17UNg75p26cBdEBjN9sB9hrOzc1V2NT+xd3sBfez2d3sBfe0uSf0h6jcDrwGhAG/liQpE9Pwpgyo62b7ZHeNZmZmqvrBVoFAoDD9ISqNbUOfFmAFsAZTr+R+QA+83M22QCBwc1RGo9HZNggEgkGE8FMRCASKIkRFIBAoSn/MqSiKJEnDcAMHOUmSJgPPtm3eBazFxW02I0nSh8AfaOeo6Mo2S5KUBtwJ/BEX/260rYIuw7S6+Qgu/r2QJOlHQAjwT3rpqOryooLJOc4dHOSSgUXAAuAAcBWmfwhXthlJkuYBk4BVuMfnDPA/wD3AE7i+zXcADwHpwH5gBi76vZAkKR2YCRyi63PnsKOqOwx/3MJBTpblv8my3AhMBHJlWT6Ci9ssSZIXsATTL5DeHT5nSZL+C0gB1gEj3MDmjZgewKXABlf+XsiyfAh4B1DRB0dVdxCVHjnIORNJkpKB4ZjinsD1bX4YeBEwAFFt+1zd5tuBqzF1yUPb9rmyzdmYhhL/Aqa07XNle83+YL12VHWH4U+eow5yzkSSJD9Mcyo/AZa4g82Y4rTiMH3ZJ0qSNAFT19yVba4GKjB9tpluYHMa8DQQBLziBvaayeuto6o7iMo63MNB7hfAUEzCshNYjovbLMvyEgBJkn6Nyc7VuLjNwKvAu0AgcB0mB0tXtvltYBum3uANuL69AEZMz90r9MJRVTi/CQQCRXGHORWBQOBGCFERCASKIkRFIBAoihAVgUCgKEJUBAKBoghREXgskiTd2cvr/No8ewVWEEvKLoAkSb8BpgEXgSOyLP+Pcy0yIUnSTOATWZZDJEn6DNgry/JTNs7dKsvy7G7aygKyZFl+sY82/V2W5dv70kZbO/MwZSf0A1bKspwuSdIh4DlMoQsngSHAWlmWP7Vy/QpMbvf5fbVlsCF6Kq7D/8qyvBiYIElSoLONaYe3JEkJmAImu6PDr5MkSVvbb8uyvLuvgtLWTreCIknSjyVJ+rEDTc2VZfmDttdxkiTpMAXNGQGjLMt3YXJWW2bj+hcwBQsKOuEOHrWehhZ4XZIkGUiSZflOSZLGYvLEDAD+BqwH/g6cB+pkWV4tSdIjmCoZhAL3Yorr2SrL8g5zL6LtQf8SiJJleYUkSbMxPRhhwO9lWd5txZ6DmCJrzwNIkrQO+I0sy2XWeieSJIUDK4FESZLWAM/JslzV1uuZJcvyk5IkJQLLZVleJknSi5hSGEwAfoQp9uRxWZZlSZK+A/YBB2RZfq2tfcs9JUmaismTuRl4HlOvYzFglCTJT5blVyVJuhWYBUQCj8qyfKLN1IZ2Zh8FfgjkcSX2hbbPu0GSJO+2z7sOOCfL8mpZlg2SJOlt/Bt6NKKn4jr8QpKk94A3MH3hy2VZNo/5/xuowRSDMRVT7EgMJoFY23bOOEyh9U9j+vK3p30vYr8syyvaXv8GOAeUcyXYrTOHMD3shzA9cO3b6jJ2lmX5vCzLvwBOyLL8S1mWqzqfL8tyMRDd9rBGy7JcAhRiipBtAOa1nZ+MScBes3HPSuDNtr/Xy7K8EfgQ+FCW5VfbzvklcAE4iynFgzUOtr3Hg23bKkmS3sLkpv4LTM/Jm8Be4Pp214lk7FYQouI6PCPL8p2yLH/Ytr2x3TEjpl/zX2EKPQdTT0SFKQ4GTHEZe4CfA1mYYk28JUkKACLMDcmyvKldu62YcpOsAY7ZsOsUpqA48wNnbjeJjg9V5wfM3gO3BVNSqM1t289iep9fYoqEBTgky/K5btr9byAc+KTdNa107IFfkmX5l5g+v4p2+7XtXudh6o2ZP1ujLMt3t/17HMYkRvcAOZgE2Jw2QoOgC2L44xp0/sXvvP0M8FLb69fajq8CLmF6IADmA0mYhgDlwA5MwpPddp61dn8LvI/p4bA6AdvGHkw9iQhMYvBLoIjuey2FkiT9FXivk5CZ+QA4ASS0bR/HNEyJ44qAWVtFaL/vGDAZmN5u/w7gjbbh1n3Ay5IkfYBJRJa2u/ZLSZKWYBLJOmATUI91MbyA6XO9G4iVJEkLPIqpZyXohFj9EViQJOl2TFGo5mHOt21JhQYlkiTdYG1lx4Hr/ICrZFn+qh/McnuEqAgEAkURcyoCgUBRhKgIBAJFEaIiEAgURYiKQCBQFCEqAoFAUYSoCAQCRfn/Ytb7zBdq6bEAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig_pres, ax_pres = plt.subplots()\n", "ax_pres.invert_yaxis()\n", "\n", "well_cug1.hydro_log().plot(ax_pres, color='green', linestyle='--', \n", " zorder=2, label='hydrostatic')\n", "\n", "well_cug1.get_log(\"Overburden_Pressure\").plot(ax_pres, color='g', \n", " label='lithostatic')\n", "\n", "pres_log.plot(ax_pres, label='mutlivaraite', zorder=1)\n", "\n", "pres_log_filter_smooth.plot(ax_pres, label='smoothed', zorder=5, color='b')\n", "\n", "ax_pres.set(xlim=[0,100], ylim=[5000,0], aspect=(100/5000)*2)\n", "ax_pres.legend()\n", "fig_pres.set(figheight=8)\n", "fig_pres.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.15" } }, "nbformat": 4, "nbformat_minor": 2 }