{ "cells": [ { "attachments": {}, "cell_type": "markdown", "id": "f324071c", "metadata": {}, "source": [ "# Fitting with Varying Fracface Pressure\n", "\n", "More advanced use of `bluebonnet.flow` in a project where pressure variation at the frac-face or bottom-hole is approximately known.\n", "\n", "First, import the necessary libraries" ] }, { "cell_type": "code", "execution_count": 1, "id": "5d796a2f", "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pandas as pd\n", "from scipy import interpolate\n", "\n", "from bluebonnet.flow import FlowProperties, SinglePhaseReservoir\n", "from bluebonnet.fluids import build_pvt_gas\n", "from bluebonnet.forecast.forecast_pressure import (\n", " fit_production_pressure,\n", " plot_production_comparison,\n", ")" ] }, { "cell_type": "code", "execution_count": 2, "id": "6f267dec-0a18-4ba9-9b95-8341cf4668ec", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mf = 0.09720937870967307\n", "CPU times: user 718 ms, sys: 1.47 ms, total: 719 ms\n", "Wall time: 719 ms\n", "0.3642368902150243\n", "pseudopressure_fracface=array(0.09720938) 0.0008833316702989393 Deviation is -1%\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAINCAYAAADP4VPTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABruUlEQVR4nO3deVhUZf/H8fewDTuICIjiFiYuiFsiWmlmovlYLpWpj1tmqWiiuZZCVk8ubdpP09JSK/fcSk0zVHLBfUlTUVFDZXNjl2GZ8/uDmJwAWRxmWL6v65pL5sx9znwPwnw45z7nvlWKoigIIYQQD2Fm6gKEEEKUfxIWQgghiiRhIYQQokgSFkIIIYokYSGEEKJIEhZCCCGKJGEhhBCiSBIWQgghimRh6gIqAq1WS0xMDA4ODqhUKlOXI4QQj0xRFFJSUvD09MTMrOjjBgmLYoiJicHLy8vUZQghhMFdv36d2rVrF9lOwqIYHBwcgNxvqqOjo4mrEUKIR5ecnIyXl5fu860oEhbFkHfqydHRUcJCCFGpFPfUunRwCyGEKJKEhRBCiCJJWAghhCiShIUQQogiSVgIIYQoklwNJYQQlYRWq5CckcXdtEzupWdyLy0LczMVz/i4PfK2JSyEEKKcUhSFVE02d1IzuZOm4VZK7r93UjO5k6rhdlomt1M03EnL5N7fAaH910TZDd3sJSyEEKIiysjK4VaKhvjkDBL+/jc+WcPtVA13UnM//O+kZnI7VYMmW1uibVuSTS31fera3CfVrh61axjm3jAJCyGEMJDsHC0JKRrikjNI+FcQxCdnkJCsIT4lg8T0rGJv04JsPK3SaWCbTj3rdOKcWuPsaE91eyv8k3fx+J0wbLMTsc68g0XGXcwyU3NXzABePw6u3gbZNwkLIYQoppSMLGISM7iZmM7NxAxiEu8/8MggLjmDnH+fByqACi3uFmk0tk/D2zqVhOr+VHNyoIaDmrZ3t+Adux2brHtYau5innEvd6WMvx//fSAAdm+AP/cU8AZmYOMCecFhABIWQgjxt/TMbKLvphN9Jz3337vp3LiXGwY3E++TkpFdxBYUXM3SaGafSprTY1R3dMDdUU2H+3toem83jtm3sc64hcX9W6i02QUHwB4NnDmmv1mVGdhWB7sakKP5Z3mj7uBUK3e5rSvYuYKtC6idoBgjyZaEhIUQospQFIVbqRqu303nrzu5j+t30/nr72C4laIpchvOtpZ4OtnwtHUUbZU/qKm6g2tOAg6aeNRpMaiy70Mm0O8YuDbMXWnPz3Bh97+2pMr9kHdwh+yMfxY3eRFq+OS+lvewcQYz8/zF1GqV+zACCQshRKWjyc7hrzvpXLmVStStNKISUom6ncaVhFRSNIUfHViSTSPrRFo7JtLY+i4NzG/hoY2nWlYcNvdjyRy0HVtPn9zGew/A3m8K3pCtK2iS/3n+eFewrwEONcHBA+w9wN4NzC3zr+veNPdRzkhYCCEqrDRNNpHxKVyMSyHq72C4ciuV6Lvp+S4hzeOkSqWNQyLN7e7RyOoOMQ1ewq2mF3Vd7Gj453ysIz6FZHIf/2KREQf8HRa1n4BWg8GxNjg98HD0BEsb/RVrtc59VGASFkKIci87R8vV22lciEshMi4l99/4ZK7fvV/oOvZqCx6rYcdzdlE8o9lNzewbOKZexSLjTu5posy/G3Z9Hh7zzP36tjdY2EC1ev961M0NguoPXFnk/Wzuo4qQsBBClCtpmmz+jEnmzM0k/oxJ4kJsCpdvpZJZwP0GajJ5wv4WTznG09QylrpKDNUzosno9inVGnfKnavhxFX4abP+ivYe4FI/NwhsnP9Z7vsK+PUHmT45HwkLIYTJPBgMZ28m8ceNRK7cTkP51ykkM7Q4WkJ9DxcaezjQ2eIPnor6BOuUv1Bla+GufnvbtCugeib3iZc/dJwC1RvmXnFU3RvUhcwOZy4fiYWR74wQwihytAoX41M4EX2P43/d4/T1goPBhgza28fR2SmGFuZ/4ZUZhX3KFejxCWateuY2+isVTl79ewWX3A7hGj7g+nhuIHj4/bPBGo/DM+8YZycrMQkLIUSZSErP4sT1e5z86x4nohM5dT2R1HxXIil4ONrQrJYTzzjF8+KVEOxSruQeLdz5V9Pbkf987dEcBm0Ct6a5VxXJaaMyJ2EhhDCI26kaDl+5S8SV2xy+cpdLCfp3D5uhpaVVDP9xuUE7yyjqZZyDpn2w6xaS2yAlHk5dzv3a3h1qtoCaflCzee6Rg3O9fzamtofHOhtlv0QuCQshRKncSdVw+OpdIqLucOjKnXzhAODjYsbbtttopo3ELeVPzLPSIPHBjZz552sHdxi4ATya5d6LIMoVCQshRLFkZOVw+Opdfr94i32XbnExXj8c3LhHn+rXaFTDGps2/+WJetWobmsBs0dAZkpuIyv73PsNvNrm3qfw73sPGnYx0t6IkpKwEEIUSFEUom6lEX7xFuEXb3H4yh294bI9uENvl2t0sbmEj+YP7FKvQRqgbgDNJv6zoY6TQO2YGxA1fAoetkKUexIWQggdTXYOEVF3+PVcPOGRt7iZqH/TW00nazo+XoO346dSI+EApJP7AECV279Q90nIyf7nMtQO44y5C6KMSFgIUcUl3c9ib2SCLiDyrlhSoaWFxXX6V79EgMVFNH2/x7tmtdwb3bb6wK2I3A7ouh2g3lNQp53+DW6iUpGwEKIKupWiYcfZWH49F09E1B2y/x5IyYYMXrG/wKsOZ2iWfgQrzR1I+nulrEhQBeR+3XEKdAkFayfT7IAwOgkLIaqIpPQsdvwZy0+nY4iIuqM30F5DN3vG1TjO89dmY5atgb/n28HSDuo/nTsGUvXH/lnBwd2otQvTk7AQohJL02Tz2/l4fjoVw++XbpGVk5sQtVUJjKh+mmqPt6dZ++dpUMMebjrAkpngXBd8euROrOPVDiysTLwXojyQsBCiklEUhaPX7rHu2HW2n4klPTMHyL16aZjzKV60PIxHytncK5eULKjxSu6Kni1h9KHcK5bkjmjxL4add89IFi5cSL169bC2tsbf358jR448tP28efNo1KgRNjY2eHl5MX78eDIyMh66jhAVTVxSBgv3XOaZT/byylcR/Hj8BpmZGoId97LfbS6HrMfyZsY3uUGhMss9vVT/6X82oFKBW2MJClGgCndksXbtWiZMmMDixYvx9/dn3rx5BAYGEhkZiZubW772q1atYurUqXz77be0b9+eixcvMnToUFQqFZ999pkJ9kAIw9FqFcIv3eK7g9cIv3jr734IBTsrC/7T3JNX2tSi1c8foLpzKXeFOgHQtE/u1J3S7yBKQKUo/x7zsXzz9/fniSeeYMGCBQBotVq8vLwYO3YsU6dOzdd+zJgxnD9/nrCwMN2yt99+m8OHD7N///5ivWdycjJOTk4kJSXh6OhomB0R4hEkpWex/vh1vj/0F3/dyb3RoaHqBm+5HKYjJzAftQ87+7+H4T69BtJuQ9NeuRP4CEHJP9cq1JFFZmYmx48fZ9q0abplZmZmdOnShYiIiALXad++PT/88ANHjhyhbdu2XLlyhe3btzNo0KBC30ej0aDR/DNxe3JyAfMrCmEClxNS+Wb/VTafvMn9rBysyKK/9VGC7PZSO+3vfgiAa79Cs765X/u9arJ6ReVRocLi9u3b5OTk4O6uf/js7u7OhQsXClxnwIAB3L59myeffBJFUcjOzmbkyJG8807h49vPmjWLmTNnGrR2IR7F6euJLNobxc5zcSgK1CCR95x286I2DOuse7khYWYBj3eDFgOh4XOmLllUMhUqLEpj7969fPTRR3z55Zf4+/tz+fJlxo0bxwcffMCMGTMKXGfatGlMmDBB9zw5ORkvLy9jlSyEzsGo2yzcc5kDl/+Z3OG5Ju6M8XXDb8vo3AWOtaDNMGg1JHduByHKQIUKC1dXV8zNzYmPj9dbHh8fj4dHwUMaz5gxg0GDBvH6668D4OvrS1paGm+88QbvvvsuZmb5LwhTq9Wo1WrD74AQxfTHjUTm7ohk/+XbqNASaH6SnrVSadR3Bg3d/+6LuB0MtdvA491lOlBR5irUT5iVlRWtW7cmLCyMXr16Abkd3GFhYYwZM6bAddLT0/MFgrl57qiXFaxvX1QBN+6lM2v7BbadiUWFlj4Wh5hqvxW3jGtw2wLUY4G/w+I5OVUqjKdChQXAhAkTGDJkCG3atKFt27bMmzePtLQ0hg0bBsDgwYOpVasWs2bNAqBnz5589tlntGzZUncaasaMGfTs2VMXGkKYWnaOlqX7rzL/t0vcz8rmOfMTzLTfhKfmCmSQO8R3m2FgaWvqUkUVVeHCol+/fty6dYuQkBDi4uJo0aIFO3bs0HV6R0dH6x1JTJ8+HZVKxfTp07l58yY1atSgZ8+e/O9//zPVLgih51aKhqBVJzhy9S71VLEsdvgGn6xzoAHUTtBhLLR9QwbtEyZV4e6zMAW5z0KUlZuJ9xm27AgX41OxV1vw0XNu9Nz3AqqcLGg3MncuCJtqpi5TVEKV+j4LISoDrVZh/+XbfH/oL8LOx4OixcXOmnVvBuDtZg8ey8CtCTh6mrpUIXQkLIQwkqT0LNYdu84Ph3PvunYmhTfN9/Ca9W5Sn51NfTf73IbeMg+1KH8kLIQoY9fvpvPN/qusO3ad9MwcvFU3+Fj9K73M9mGpaCAHapycC/4vgLmlqcsVokASFkKUkVPXE1ny+xV+ORuLVoGnzU4zzv5XWmefzG2gAB7Nod1oaNZHgkKUaxIWQhiQouT2R/xf2GWOXLurW/5UQ1fmZ+6jWvzJ3OHBfXrkhkSdABkSXFQIEhZCGICiKIRfvMX8sEucjE7EiVTGWoaR0vhV+j3ThsY1HSFqMlz6FfzfhGr1TF2yECUiYSHEI1AUhT2RCcwPu8zp64nUIJEQq2381yIMK20GeNSBmp1zGz/2TO5DiApIwkKIUjpy9S6zfjnPyehEanCP96y2MtA8DEslE7SAuy94+Jq6TCEMQsJCiBK6GJ/C3B0X+O18AgAzrFYz2PzX3CubFKBWG+g0NfcSWOmPEJWEhIUQxRSXlMFnuyL58fgNtAqYm6no94QXA5RqWP6hAS9/6DgFHussISEqHQkLIYqgyc5h6b6rLNh9mZysDIaZ70Lb4Bn++0I3HqthD0lTwK8PNOgkISEqLQkLIR5i94V43v/5HNfupPGC2UHetf0Rd2082NyGGi/lNnKqLXNbi0rPKGFx/fp1VCoVtWvn/kIdOXKEVatW0aRJE9544w1jlCBEiVy7ncb7W8+x+0ICTVXX2GKzAj8lMrfj2qEmNHoeFEWOJESVYZSwGDBgAG+88QaDBg0iLi6O5557jqZNm7Jy5Uri4uIICQkxRhlCFCkrR8uSfVeY99slbLOT+MhyPf3Nw1ApCljawVPjoV0QWMm8EqJqyT+naBk4e/Ysbdu2BWDdunU0a9aMgwcPsnLlSpYvX26MEoQo0pkbSby44ABzd0SSma1lsvsxBpj/hgoFmvWFMUfh6UkSFKJKMsqRRVZWlm5O699++40XXngBAB8fH2JjY41RghCFup+Zw7zfLrJk3xVQtDjbqpnRowl9mj8LG2Kh3Sio96SpyxTCpIwSFk2bNmXx4sX06NGDXbt28cEHHwAQExND9erVjVGCEAU6dT2RCWtPcf12EmPMf6KPw1nsR4Xh6vT3cOGvrjRtgUKUE0YJizlz5tC7d28+/vhjhgwZgp+fHwA//fST7vSUEMaUnaNl4Z4ovth9iebKRZbYLOUx5XrufNfRO8D3JVOXKES5YrRpVXNyckhOTqZatX+miLx27Rq2tra4ubkZo4RSk2lVK5drt9MYv+4UZ6NvM85iA6MtfsYMLdi6wvNzoWkfucpJVHrldlpVc3NzvaAAqFevnrHeXggA1h+7TuhPf1IzK5ot1l/ShKu5LzTvB91mg62LaQsUopwyWlj8+OOPrFu3jujoaDIzM/VeO3HihLHKEFXU/cwcZmw5y4/HbwAw33ElTTKvgk01+M88aNrLpPUJUd4Z5dLZL774gmHDhuHu7s7Jkydp27Yt1atX58qVK3Tv3t0YJYgq7HJCKr0WHuDH4zcwU8HEro/T5I1vwec/MCpCgkKIYjBKn4WPjw+hoaH0798fBwcHTp8+TYMGDQgJCeHu3bssWLCgrEt4JNJnUXFtOXWTaRvP0DjrHB1s/iJgYAgBj8kVeEKU9HPNKEcW0dHRtG/fHgAbGxtSUlIAGDRoEKtXrzZGCaKKyc7R8uHWc4xbc5KBOVtYp/6ACdrlBKjOmLo0ISoko4SFh4cHd+/mzkdcp04dDh06BMDVq1cx0sVYogpJSs/itRXHWLf/DEssP+Vdy1WYo4VmL0Gt1qYuT4gKySgd3J07d+ann36iZcuWDBs2jPHjx/Pjjz9y7Ngx+vTpY4wSRBVxOSGVEd8dw+zORX5Sf0o9VRyYW0H3OdB6mFwSK0QpGaXPQqvVotVqsbDIzaY1a9Zw8OBBGjZsyJtvvomVlVVZl/BIpM+iYtgbmcDYVSdpnXWMBVYLsCcdnLyg3w/g2cLU5QlRrpTL+yzMzMwwM/vnjNerr77Kq6++aoy3FlXE2qPRvLPpLDlahbZumdgnp0Od9vDKd2Bfw9TlCVHhGSUs/vjjjwKXq1QqrK2tqVOnjm6gQSFKQlEU5v12iflhlwDo07IWr/cNhQvNwacnWJTvo1YhKgqjhEWLFi1QPXCuWFEUveeWlpb069ePr776Cmtra2OUJCqBrBwt0zed5ZdjF5htsZJb7d5hTA+/3J+tZn1NXZ4QlYpRrobatGkTDRs25Ouvv+bUqVOcPn2ar7/+mkaNGrFq1Sq++eYbdu/ezfTp041RjqgEMrJyeOO7Y/x+7BTrrd7nVYu9jE2cq/dHiBDCcIxyZPG///2P+fPnExgYqFvm6+tL7dq1mTFjBkeOHMHOzo63336bTz75xBgliQosVZPN6yuOcu/qKTap5+Khugv27vBsqKlLE6LSMsqRxZkzZ6hbt26+5XXr1uXMmdybpFq0aFHsiZAWLlxIvXr1sLa2xt/fnyNHjjy0fWJiIkFBQdSsWRO1Ws3jjz/O9u3bS74jwuSS0rMY9M1hsq8eZL3V+7lB4fo4vP4b1Gxu6vKEqLSMEhY+Pj7Mnj1bbwDBrKwsZs+ejY+PDwA3b97E3d29yG2tXbuWCRMmEBoayokTJ/Dz8yMwMJCEhIQC22dmZvLcc89x7do1fvzxRyIjI1myZAm1atUyzM4Jo7mTqqH/kkPY3tjH9+rZOKrSoU4AvLYTnOuYujwhKjWjnIZauHAhL7zwArVr16Z589y//s6cOUNOTg5bt24F4MqVK4wePbrIbX322WeMGDGCYcOGAbB48WK2bdvGt99+y9SpU/O1//bbb7l79y4HDx7E0tISkKHRK6KElAwGLDnM1YQkdluvwIZM8O4Cr3wvc2ILYQRGm/woJSWFlStXcvHiRQAaNWrEgAEDcHBwKPY2MjMzsbW15ccff6RXr1665UOGDCExMZEtW7bkW+f555/HxcUFW1tbtmzZQo0aNRgwYABTpkzB3Ny8wPfRaDRoNBrd8+TkZLy8vOSmPBO5k6rh1a8PcSkhFQ9Ha9a+7EbdC0uh+1ywkEuuhSiNcnlTHoCDgwMjR458pG3cvn2bnJycfKer3N3duXDhQoHrXLlyhd27dzNw4EC2b9/O5cuXGT16NFlZWYSGFtwhOmvWLGbOnPlItQrDSEzP5L/fHOFuwk3cHd1Y+2Y76la3g4bzTV2aEFWKUfosTEmr1eLm5sbXX39N69at6devH++++y6LFy8udJ1p06aRlJSke1y/ft2IFYs8yRlZDP72CA3if2W/dTCbuyTnBoUQwuiMdmRhCK6urpibmxMfH6+3PD4+Hg8PjwLXqVmzJpaWlnqnnBo3bkxcXByZmZkFjkulVqvljnITS8/MZtiyo3jE/MZ8qwVYoMUmYR8gA08KYQoV6sjCysqK1q1bExYWplum1WoJCwsjICCgwHU6dOjA5cuX0Wq1umUXL16kZs2a5X4Aw6oqK0dL0MoTOFzfzUKrL7BAC81fheflHhwhTKVChQXAhAkTWLJkCStWrOD8+fOMGjWKtLQ03dVRgwcPZtq0abr2o0aN4u7du4wbN46LFy+ybds2PvroI4KCgky1C+IhFEVh6oYzZFzay1eW87AkB5r2hhcXglnBFyQIIcpehToNBdCvXz9u3bpFSEgIcXFxtGjRgh07dug6vaOjo/VGuPXy8mLnzp2MHz+e5s2bU6tWLcaNG8eUKVNMtQviIT7eGcn5k/tZa/UZalUWNOoBfZaAeYX7URWiUimzS2erVatW7HF68mbRK69kPgvjWHHwGqE//clMi2UMsdgFdZ+E/24ASxlcUghDKzeXzs6bN0/39Z07d/jwww8JDAzU9S1ERESwc+dOZsyYUVYliApkb2QCM3/+E4DkTv8D6yehzTAJCiHKCaPclNe3b1+eeeYZxowZo7d8wYIF/Pbbb2zevLmsS3gkcmRRti4npPDywt+5p1F4ubUXc19qLqPHClHGSvq5ZpSwsLe359SpU3h7e+stv3z5Mi1atCA1NbWsS3gkEhZl515aJr0X7mNiyhwc7ezwD/4BtbXcSyFEWSvp55pRroaqXr16gUNxbNmyherVqxujBFEOZeVoGbXyOC8mreQ/5od5KnMf6tsF34kvhDAto1xiMnPmTF5//XX27t2Lv78/AIcPH2bHjh0sWbLEGCWIcuiTnZG4XNvOeKsNAKj+8xnUbm3iqoQQBTFKWAwdOpTGjRvzxRdfsHHjRiD3Lur9+/frwkNULTv/jGP/vt38aPX3sCvtgqDVYNMWJYQolNFGna3IpM/CsP66k8aA/9vBau0U6pjdyh1qfMA6uelOCCMql30WAFFRUUyfPp0BAwboJir65Zdf+PPPP41VgigHMrJyGL3yBNNzFlHH7BaKcx3ou1SCQohyzihhER4ejq+vL4cPH2bDhg26q59Onz5d6DDhonL69NdI/oxJZotld7IdaqF65TuwqWbqsoQQRTBKWEydOpUPP/yQXbt26Q3e17lzZw4dOmSMEkQ5EBF1h6X7rwLQ96X/YjHuFHi2NG1RQohiMUpYnDlzht69e+db7ubmxu3bt41RgjCx5IwsZqw7RE3lNq8+4cVzTdzBQkb9FaKiMEpYODs7Exsbm2/5yZMnqVWrljFKECb20bbzDEv7hl+tpxLqfdnU5QghSsgoYfHqq68yZcoU4uLiUKlUaLVaDhw4wMSJExk8WC6XrOyOXbtL3PGfGWgRhj3p2DjIjZhCVDRGCYuPPvoIHx8fvLy8SE1NpUmTJjz99NO0b9+e6dOnG6MEYSJZOVr+t/EIsy3/vvmy3Who0NG0RQkhSqzM77NQFIXr169To0YNbt++zZkzZ0hNTaVly5Y0bNiwLN/aYOQ+i9L7KjwKi13vMtziF3Kc62EedAgsbUxdlhBVXrkZojyPoih4e3vz559/0rBhQ7y8vMr6LUU5kZieya49v7HWfAcA5v/5TIJCiAqqzE9DmZmZ0bBhQ+7cuVPWbyXKmUV7LzNJ+y3mKgWlSW/wftbUJQkhSskofRazZ89m0qRJnD171hhvJ8qBhOQMVh68zAWtF9kWdqgC/2fqkoQQj8AoAwkOHjyY9PR0/Pz8sLKywsZG/1REeZ9WVZTcD4f+IjXbnC11xjP4tW/A2snUJQkhHoFRwuLBKVZF5afJzmHVkWgAXn+qASoJCiEqPKOExZAhQ4zxNqKcOHD5Ni/e38xlO1+6Nulu6nKEEAZg9FFn+/fvL6POVnKHTp3hXYuVrMiZgkXKDVOXI4QwAKOPOrtx40YZdbaSs/lrN2YqhaTqLcG5jqnLEUIYgIw6KwzqckIKTdOO5D5p2MW0xQghDEZGnRUGc/ZmEp8s+opnVccAcPDtYeKKhBCGIqPOCoM4du0uS7+ez3ztLMxVCpkNnsPMs4WpyxJCGIiMOise2f5Lt3nzm995lyWoVdlkNeqJ1YCVoFKZujQhhIHIqLPikew6F89ry49yJ8uShR7/I7vNG1j2WwEWalOXJoQwoDIfdfZB0dHRnD17VkadrSR2nI3jvVW7idM6EdjUnS/6t0RtYW7qsoQQxVDuRp0F2L9/P08++SR16tShTh25lLIyCDsfzx9r3+NXy80srfcpbw3ojoW50W7bEUIYmVF+uzt37kz9+vV55513OHfunDHeUpShvZEJnFwVymTz1Tiq7jPusXgJCiEqOaP8hsfExPD2228THh5Os2bNaNGiBR9//DE3bsjdvRXN/ku3OfHDdCaarwYgp+M7mD89wcRVCSHKmlHCwtXVlTFjxnDgwAGioqJ4+eWXWbFiBfXq1aNz587GKEEYwKErd9j//XtMMF8LQHbHdzF/ZoqJqxJCGIPRzx3Ur1+fqVOnMnv2bHx9fQkPDy/VdhYuXEi9evWwtrbG39+fI0eOFGu9NWvWoFKp6NWrV6net6o6ezOJbSvmMtXsewCyO07D4pnJJq5KCGEsRg2LAwcOMHr0aGrWrMmAAQNo1qwZ27ZtK/F21q5dy4QJEwgNDeXEiRP4+fkRGBioG6CwMNeuXWPixIk89dRTpd2FKumvO2kM/fYIXbQHAchuNwaLTnJEIURVYpRLZ6dNm8aaNWuIiYnhueeeY+DAgbz44ovY2tqWanv+/v488cQTLFiwAACtVouXlxdjx45l6tSpBa6Tk5PD008/zWuvvca+fftITExk8+bNxXq/qnzpbEJKBi8tiiD6bjrNPWxY6x+FTbvX5YY7ISq4cnnp7O+//86kSZN45ZVXcHV1faRtZWZmcvz4caZNm6ZbZmZmRpcuXYiIiCh0vffffx83NzeGDx/Ovn37HvoeGo0GjUaje56cnPxINVdUyRlZTFy6nei7VtRxsWPp8ABsHKSPSYiqyChhceDAAYNt6/bt2+Tk5ODu7q633N3dnQsXLhS4zv79+/nmm284depUsd5j1qxZzJw581FLrdA02Tm8980mvkgMZpvNMzz52le4OVibuiwhhIkYdfKjsWPH0qVLF7p06cJbb71FVFRUmb9vSkoKgwYNYsmSJcU+qpk2bRpJSUm6x/Xr18u4yvJFURQ+WLOX8fHv4KxKo3eNGOo6GeXvCiFEOWWUT4CdO3fywgsv0KJFCzp06ADkHm00bdqUn3/+meeee67Y23J1dcXc3Jz4+Hi95fHx8Xh4eORrHxUVxbVr1+jZs6dumVarBcDCwoLIyEgee+wxvXXUajVqddUd2+jLX8/w8sWJeJnd4r59HWyHrAdLOaoQoiozSlhMnTqV8ePHM3v27HzLp0yZUqKwsLKyonXr1oSFhekuf9VqtYSFhTFmzJh87X18fDhz5ozesunTp5OSksL8+fPx8vIq+Q5VYj+d+ItG+9/Cz/wKGktnbIZtBrtH62cSQlR8RgmL8+fPs27dunzLX3vtNebNm1fi7U2YMIEhQ4bQpk0b2rZty7x580hLS2PYsGEADB48mFq1ajFr1iysra1p1qyZ3vrOzs4A+ZZXdceu3iFl09u8YH6SLJUa9eD1UP2xolcUQlR6RgmLGjVqcOrUqXyjzJ46dQo3N7cSb69fv37cunWLkJAQ4uLiaNGiBTt27NB1ekdHR2NmJmMVlUT0nXTmf7+W7813oUWF+Utfg1dbU5clhCgnjHKfxfvvv8/nn3/O1KlTad++PZDbZzFnzhwmTJjAjBkzyrqER1LZ77NIzsiiz5cHuZyQyoga55gUYI9V+1GmLksIUYZK+rlmlLBQFIV58+bx6aefEhMTA4CnpyeTJk3irbfeQlXOb/CqzGGh1Sq88f0xfjufgIejNVvGdMDdUTqzhajsymVYPCglJQUABwcHY77tI6nMYbHwlxPUPTiNj5X/8n8je9K8trOpSxJCGEG5vIP76tWrZGdn07BhQ72QuHTpEpaWltSrV88YZYh/+fXsTZoeDKaT+WnaO6XgUuu/pi5JCFFOGaUXeOjQoRw8eDDf8sOHDzN06FBjlCD+5XJCKjfXT6GT+WmyVGpcXl0k4z0JIQpllLA4efKk7ma8B7Vr167YQ3AIw0nOyGLNt58xTPUzAKo+i6Gmn4mrEkKUZ0YJC5VKpeureFBSUhI5OTnGKEH8TatV+OT7LUy4nztib7r/OCx8+5i4KiFEeWeUsHj66aeZNWuWXjDk5OQwa9YsnnzySWOUIP62fM8fDLk+HVuVhhTPDtgGhpq6JCFEBWCUDu45c+bw9NNP06hRI93EQ/v27SM5OZndu3cbowQBHL12l2/C/qCthZo0a3ccBn4HZuamLksIUQEY5ciiSZMm/PHHH7zyyiskJCSQkpLC4MGDuXDhggy5YSR3UjWMXXWSm1oXlvt8he3wrTLmkxCi2Ix+n0VFVNHvs9BqFUZ+G86vl9NoUMOOn8c8iZ1ahhwXoior6eea0QZQ2rdvH//9739p3749N2/eBOD7779n//79xiqhylqx6zDvXx/K21Yb+XJACwkKIUSJGSUsNmzYQGBgIDY2Npw4cUI3ZWlSUhIfffSRMUqosiIu36LhgbfxUN1jsNMpfKpbmbokIUQFZJSw+PDDD1m8eDFLlizB0tJSt7xDhw6cOHHCGCVUSbdSNJxYFcqTZmfRqKxxGrwKrGxNXZYQogIySlhERkby9NNP51vu5OREYmKiMUqocrRahcU/rOLNnNW5C56fCzUambYoIUSFZZSw8PDw4PLly/mW79+/nwYNGhijhCpnZfhphsV9iIVKS3LDXqjbDDZ1SUKICswoYTFixAjGjRvH4cOHUalUxMTEsHLlSiZOnMioUTJvgqH9eTMR9z0Tqa26TYqtF459/0/GfRJCPBKjzcGt1Wp59tlnSU9P5+mnn0atVjNx4kTGjh1rjBKqjPuZOby1+iStc1rQweI89gO/B+uKd7mvEKJ8Mep9FpmZmVy+fJnU1FSaNGmCvb099+/fx8bGxlgllEpFus/inU1nWHU4GndHNb+MbImLS3VTlySEKIfK7X0WAFZWVjRp0oS2bdtiaWnJZ599Rv369Y1ZQqW269QVth4+h0oFn73SQoJCCGEwZRoWGo2GadOm0aZNG9q3b8/mzZsBWLZsGfXr1+fzzz9n/PjxZVlClRGbdJ/EzZPZrn6HmS3T6eAtQ3kIIQynTPssQkJC+Oqrr+jSpQsHDx7k5ZdfZtiwYRw6dIjPPvuMl19+GXNzGcjuUeVoFb5b8RVT2AUq6N+yhqlLEkJUMmUaFuvXr+e7777jhRde4OzZszRv3pzs7GxOnz6NSq7OMZjlvx3jtTufgQqS/Ebg1LCTqUsSQlQyZXoa6saNG7Ru3RqAZs2aoVarGT9+vASFAZ29kUit/e9QQ5VEkn0DnP7zgalLEkJUQmUaFjk5OVhZ/TMWkYWFBfb29mX5llVKRlYO2374nG5mR8jBHMcB34Jl+b6yTAhRMZXpaShFURg6dChqtRqAjIwMRo4ciZ2dnV67jRs3lmUZldaSn8MZdf8rUIHmycnYerY0dUlCiEqqTMNiyJAhes//+9//luXbVSmHrtxh2dF4Glo0I8AtC6dnJpq6JCFEJVamYbFs2bKy3HyVlZKRxdvrTnNXcWSv36d0+099MJc5KgwhJyeHrKwsU5chhEFYWVlhZmaY3gb5hKmA5mw+ys3E+3i52DC9Z1OQyYwemaIoxMXFySjIolIxMzOjfv36en3HpSWfMhXMrjPRvHpuFP6WHni+sAh7CQqDyAsKNzc3bG1t5Yo9UeFptVpiYmKIjY2lTp06j/wzLZ80FcjtVA1/bXyP58yu0cDiHra1ZCIjQ8jJydEFRfXqMkSKqDxq1KhBTEwM2dnZehPPlYZRx4YSpacoCktXrWGYNvfKMcsX54GDh2mLqiTy+ihsbSV8ReWSd/opJyfnkbclYVFBbDp8kVdv/A9zlUKidx8sffuYuqRKR049icrGkD/TFTYsFi5cSL169bC2tsbf358jR44U2nbJkiU89dRTVKtWjWrVqtGlS5eHti9vrt9NJ/uXd6hnFk+K2h3nvp+buiQhRBVTIcNi7dq1TJgwgdDQUE6cOIGfnx+BgYEkJCQU2H7v3r3079+fPXv2EBERgZeXF127duXmzZtGrrzktFqFH75fwiuq3wCwfeVrsHE2bVFCiCK99957tGjRokTrXLhwgXbt2mFtbV3idctahQyLzz77jBEjRjBs2DCaNGnC4sWLsbW15dtvvy2w/cqVKxk9ejQtWrTAx8eHpUuXotVqCQsLM3LlJfftgav8Ea8hTnEhucUIzB/rZOqSRDkydOhQVCpVvkdBc96L8i80NBQ7OzsiIyMJCwtj+fLlODs7m7osoAJeDZWZmcnx48eZNm2abpmZmRldunQhIiKiWNtIT08nKysLFxeXAl/XaDRoNBrd8+Tk5EcrupQuxqcwd2ckmdqm7HvuJ1729zZJHaJ869atW74bYGvUyD9MfWZmpkGuty8tU79/YcpTXVFRUfTo0YO6deuaupR8KtyRxe3bt8nJycHd3V1vubu7O3FxccXaxpQpU/D09KRLly4Fvj5r1iycnJx0Dy8vr0euu6Qys7VMWnOUzGwtzzSqwUsdmskggUakKArpmdkmeZR0pmO1Wo2Hh4few9zcnE6dOjFmzBiCg4NxdXUlMDAQgPDwcNq2bYtaraZmzZpMnTqV7Oxs3fY6derE2LFjCQ4Oplq1ari7u7NkyRLS0tIYNmwYDg4OeHt788svvzy0rnr16vHBBx8wePBgHB0deeONNwDYv38/Tz31FDY2Nnh5efHWW2+RlpamW+/LL7+kYcOGWFtb4+7uzksvvaR7TaPR8NZbb+Hm5oa1tTVPPvkkR48e1b1e0F/imzdv1uvozTs9tHTpUurXr4+1tTUAiYmJvPnmm7i7u2NtbU2zZs3YunWrbr2i6i6OpUuX0rhxY6ytrfHx8eHLL7/UvaZSqTh+/Djvv/8+KpWKTp06MWzYMJKSknRHjO+9916J3s+QKtyRxaOaPXs2a9asYe/evbofkn+bNm0aEyZM0D1PTk42emAs276PRXeH85XNqwT1CZErdYzsflYOTUJ2muS9z70fiK2VYX41V6xYwahRozhw4AAAN2/e5Pnnn2fo0KF89913XLhwgREjRmBtba33QbRixQomT57MkSNHWLt2LaNGjWLTpk307t2bd955h88//5xBgwYRHR390EuOP/nkE0JCQggNDQVy/3Lu1q0bH374Id9++y23bt1izJgxjBkzhmXLlnHs2DHeeustvv/+e9q3b8/du3fZt2+fbnuTJ09mw4YNrFixgrp16zJ37lwCAwO5fPlyoWcKCnL58mU2bNjAxo0bMTc3R6vV0r17d1JSUvjhhx947LHHOHfunG5ytqLqLo6VK1cSEhLCggULaNmyJSdPnmTEiBHY2dkxZMgQYmNj6dKlC926dWPixInY2tqybNkyQkJCiIyMBDDpqN0VLixcXV0xNzcnPj5eb3l8fDweHg+/7+CTTz5h9uzZ/PbbbzRv3rzQdmq1WjdSrimc+OsOvsem4ml2l+BqB6jmUD4OkUX5tHXrVr0Pke7du7N+/XoAGjZsyNy5c3Wvvfvuu3h5ebFgwQJUKhU+Pj7ExMQwZcoUQkJCdOMI+fn5MX36dCD3j6fZs2fj6urKiBEjgNxZMBctWsQff/xBu3btCq2tc+fOvP3227rnr7/+OgMHDiQ4OFhX3xdffEHHjh1ZtGgR0dHR2NnZ8Z///AcHBwfq1q1Ly5a5oymnpaWxaNEili9fTvfu3YHcKx137drFN998w6RJk4r9PcvMzOS7777Tna779ddfOXLkCOfPn+fxxx8HoEGDBrr2s2bNemjdhf3h+aDQ0FA+/fRT+vTJvey9fv36nDt3jq+++oohQ4bg4eGhm8Yh77PMyckJlUpV5GebMVS4sLCysqJ169aEhYXRq1cvAF1n9ZgxYwpdb+7cufzvf/9j586dtGnTxkjVllx6ZjYHV37IGLNzaFTWVBu4DMxk6lljs7E059z7gSZ775J45plnWLRoke75g1MA5E0+luf8+fMEBAToHal26NCB1NRUbty4QZ06dQD0/pgyNzenevXq+Pr66pblnQYu7ArEPP/+XTt9+jR//PEHK1eu1C1TFAWtVsvVq1d57rnnqFu3Lg0aNKBbt25069aN3r17Y2trS1RUFFlZWXTo0EG3rqWlJW3btuX8+fMPrePf6tatq9evc+rUKWrXrq0Lin8rqu7GjRs/9P3S0tKIiopi+PDhusAFyM7OxsnJqUS1m0qFCwuACRMmMGTIENq0aUPbtm2ZN2+e7nwqwODBg6lVqxazZs0CYM6cOYSEhLBq1Srq1aun69uwt7cvd5MxfbNhO29ovgcVaLt+CNUfM3VJVZJKpTLYqaCyZmdnh7d3wRc//HvumOL699AQKpVKb1le2Gi12iJre1Bqaipvvvkmb731Vr62derUwcrKihMnTrB3715+/fVXQkJCeO+99/T6JR7GzMwsX59PQaMI/7suG5uH9wcWVXdRUlNTgdwjIX9/f73X8k51lXcV47fhX/r168etW7cICQkhLi6OFi1asGPHDt1fO9HR0XrD8i5atIjMzEy9jjLIPSw0ZYfRv+07f5Nnzs9AbZbFXc9OuLR73dQliUqmcePGbNiwAUVRdB/4Bw4cwMHBgdq1a5f5+7dq1Ypz584VGm6QO6Nmly5d6NKlC6GhoTg7O7N7924CAwOxsrLiwIEDuquFsrKyOHr0qO70UI0aNUhJSSEtLU0XCKdOnSqyrubNm3Pjxg0uXrxY4NFFcep+GHd3dzw9Pbly5QoDBw4s9npWVlYGGarDECpkWAC6zqWC7N27V+/5tWvXyr6gR5SUnsXl9SEMM7tGmrkTLv2/AunUFgY2evRo5s2bx9ixYxkzZgyRkZGEhoYyYcIEg8178DBTpkyhXbt2jBkzhtdffx07OzvOnTvHrl27WLBgAVu3buXKlSs8/fTTVKtWje3bt6PVamnUqBF2dnaMGjWKSZMm4eLiQp06dZg7dy7p6ekMHz4cAH9/f2xtbXnnnXd46623OHz4MMuXLy+yro4dO/L000/Tt29fPvvsM7y9vblw4QIqlYpu3boVWXdxzJw5k7feegsnJye6deuGRqPh2LFj3Lt3T++CmgfVq1eP1NRUwsLC8PPzw9bW1mRjmFW4S2crq9CfzpKemYkWFRYvzJNBAkWZqFWrFtu3b+fIkSP4+fkxcuRIhg8fruvMLmvNmzcnPDycixcv8tRTT9GyZUtCQkLw9PQEwNnZmY0bN9K5c2caN27M4sWLWb16NU2bNgVyr2bs27cvgwYNolWrVly+fJmdO3dSrVo1AFxcXPjhhx/Yvn07vr6+rF69uthnDzZs2MATTzxB//79adKkCZMnT9b9VV9U3cXx+uuvs3TpUpYtW4avry8dO3Zk+fLl1K9fv9B12rdvz8iRI+nXrx81atTQu1jB2FRKSS/qroKSk5NxcnIiKSkJR0dHg2//lzOxjFp5AjMV/NzfjabNnzD4e4jCZWRkcPXqVb1r7oWoDB72s13Sz7UKexqqskhIyeDdjacBGNXpMZo29zFxRUIIkZ+chjIhRVFY/cNSlua8Sxe3FMY9W/Ble0IIYWoSFia0JeIM/eM+ppXZZebUO46Vhfx3CCHKJ/l0MpHrd9Kw2TkJN1Uid23rU73n+6YuSQghCiVhYQJarcKW7z8nUHWIbMxxGrBMBgkUQpRrEhYmsH53BIPv5V6bneL/Nua1W5q4IiGEeDgJCyO7HJ9M3d8n4qi6zy2n5lTrOsXUJQkhRJEkLIwoO0fLe2v3YUs6GpU1roOWgblcvSyEKP/kk8qIvtwbxf4YOGf9P3a9Wo3qrjLznRCiYpAjCyM5cyOJL8IuARDyYguqN2pv4oqEKFvp6en07dsXR0dHVCoViYmJRn3/r7/+Gi8vL8zMzJg3b55R3/tBmzdvxtvbG3Nzc4KDg8vVvNolIWFhBBlZOZz8bjLBZmvo2aw6L7Yo/ngyQjzM0KFDdVNuPvi4fPmyqUtjxYoV7Nu3j4MHDxIbG2vUeRuSk5MZM2YMU6ZM4ebNm7opXU3hzTff5KWXXuL69et88MEHJqvjUclpKCNYs/FHBmnWYm6hkOLXT6ZIFQbVrVu3fFN7PjixT57MzEysrIw362JUVBSNGzemWbNmhbYpq5qio6PJysqiR48e1KxZ0+DbL67U1FQSEhIIDAws0aCD5ZEcWZSxo5HX6fjnDMxVCrF1XsDBt4epSxIlkZlW+CMrowRt7xevbSmo1Wo8PDz0Hubm5nTq1IkxY8YQHByMq6srgYG5M/+Fh4fTtm1b1Go1NWvWZOrUqWRnZ+u216lTJ8aOHUtwcDDVqlXD3d2dJUuW6CYYc3BwwNvbm19++aXQmjp16sSnn37K77//jkqlolOnTkDukNsffPABgwcPxtHRUfcX/5QpU3j88cextbWlQYMGzJgxI9+kRT///DNPPPEE1tbWuLq60rt37wLfe/ny5bpZ/Ro0aIBKpeLatWtERUXx4osv4u7ujr29PU888QS//fab3roajYYpU6bg5eWFWq3G29ubb775Rvf62bNn6d69O/b29ri7uzNo0CBu375dYB179+7FwcEByJ1eVqVS5Zs+Ic+iRYt47LHHsLKyolGjRnz//fe61yZOnMh//vMf3fN58+ahUqnYsWOHbpm3tzdLly4tcNuGImFRhlI12dxc/zb1zeJJtKhBzf7/Z+qSREl95Fn4Y90g/bYfexfe9gf9ibeY51twOwNbsWKFbsKgxYsXc/PmTZ5//nmeeOIJTp8+zaJFi/jmm2/48MMP863n6urKkSNHGDt2LKNGjeLll1+mffv2nDhxgq5duzJo0CDS09MLfN+NGzcyYsQIAgICiI2NZePGjbrXPvnkE/z8/Dh58iQzZswAwMHBgeXLl3Pu3Dnmz5/PkiVL+Pzzz3XrbNu2jd69e/P8889z8uRJwsLCaNu2bYHv3a9fP10IHDlyhNjYWLy8vEhNTeX5558nLCyMkydP0q1bN3r27El0dLRu3cGDB7N69Wq++OILzp8/z1dffaWbTTMxMZHOnTvTsmVLjh07xo4dO4iPj+eVV14psI727dsTGRkJ5A5/HhsbS/v2+fsqN23axLhx43j77bc5e/Ysb775JsOGDWPPnj1A7lwb+/fv1w2XHh4ejqurqy54bt68SVRUlC6Qy4wiipSUlKQASlJSUonWW7ZssaKEOipKqKOSdv63MqpOPKr79+8r586dU+7fv5//xb///wp8/PCSftsPPQpv++3z+m3n1C+4XQkNGTJEMTc3V+zs7HSPl17Kratjx45Ky5Yt9dq/8847SqNGjRStVqtbtnDhQsXe3l7JycnRrffkk0/qXs/Ozlbs7OyUQYMG6ZbFxsYqgBIREVFobePGjVM6duyot6xu3bpKr169ityvjz/+WGndurXueUBAgDJw4MAi18tz8uRJBVCuXr360HZNmzZV/u///k9RFEWJjIxUAGXXrl0Ftv3ggw+Url276i27fv26AiiRkZEFrnPv3j0FUPbs2aNbtmzZMsXJyUn3vH379sqIESP01nv55ZeV559/XrcNMzMz5ejRo4pWq1VcXFyUWbNmKf7+/oqiKMoPP/yg1KpVq8D3f9jPdkk/16TPoozsP32B569+BCqI9RlKTZ9nTV2SKI13Ygp/TfWvuZMnPaRTWfWvg/jgM6Wv6V+eeeYZFi1apHv+4PzSrVu31mt7/vx5AgIC9PrNOnToQGpqKjdu3NDNJ928eXPd6+bm5lSvXl13agfQTWGckJBQ4nrbtGmTb9natWv54osviIqKIjU1lezsbL05Fk6dOsWIESNK/F4PSk1N5b333mPbtm3ExsaSnZ3N/fv3dUcWp06dwtzcnI4dOxa4/unTp9mzZ4/uSONBUVFRBU7HWhznz5/P1wHfoUMH5s+fD+ROCOXn58fevXuxsrLCysqKN954g9DQUFJTUwkPDy+0ZkOSsCgD99IyWf3TNuZyn1vW9ajZd7apSxKlZWVXdJuyblsEOzu7QueGfjA4SsLS0lLvuUql0luWFzZarbbE2/53TREREQwcOJCZM2cSGBiIk5MTa9as4dNPP9W1sbF59LHTJk6cyK5du/jkk0/w9vbGxsaGl156iczMzGK9R2pqKj179mTOnDn5XivrTvROnTqxd+9e1Go1HTt2xMXFhcaNG7N//37Cw8N5++23y/T9QcKiTMzYcpZtaT4kVZ/PNwOaySCBotxo3LgxGzZsQFEU3Qf+gQMHcHBwoHbt2iap6eDBg9StW5d3331Xt+yvv/7Sa9O8eXPCwsIYNmxYqd/nwIEDDB06VNcxnpqayrVr13Sv+/r6otVqCQ8Pp0uXLvnWb9WqFRs2bKBevXpYWBjuo7Nx48YcOHCAIUOG6NXapEkT3fOOHTvy7bffYmFhQbdu3YDcAFm9ejUXL14s+/4KpIPb4H46HcPWP2IxN1Mx6dXuqGv5Fr2SEEYyevRorl+/ztixY7lw4QJbtmwhNDSUCRMmYGZmmo+Dhg0bEh0dzZo1a4iKiuKLL75g06ZNem1CQ0NZvXo1oaGhnD9/njNnzhT4F35R77Nx40ZOnTrF6dOnGTBggN6RUb169RgyZAivvfYamzdv5urVq+zdu5d169YBEBQUxN27d+nfvz9Hjx4lKiqKnTt3MmzYMF3nc2lMmjSJ5cuXs2jRIi5dusRnn33Gxo0bmThxoq7N008/TUpKClu3btUFQ6dOnVi5ciU1a9Ys9SmwkpCwMKD4pHSUTSNpZ3aOoGe88fNyNnVJQuipVasW27dv58iRI/j5+TFy5EiGDx/O9OnTTVbTCy+8wPjx4xkzZgwtWrTg4MGDuquk8nTq1In169fz008/0aJFCzp37syRI0dK9D6fffYZ1apVo3379vTs2ZPAwEBatWql12bRokW89NJLjB49Gh8fH0aMGEFaWu4lzZ6enhw4cICcnBy6du2Kr68vwcHBODs7P1LQ9urVi/nz5/PJJ5/QtGlTvvrqK5YtW6Z3tFCtWjV8fX2pUaMGPj65Uy8//fTTaLVao/RXAKgURVGM8k4VWHEmNlcUhVVfTGPgvUWkY4PlxD+xtK9u5EpFaTxsUnshKrKH/WwX53PtQXJkYSDb94TT927uTTGpT82QoBBCVCoSFgYQnZBE3fDxWKuyuO7SHrfOo01dkhBCGJSExSPK0SocWTGVZqorpKgc8BzyDcjYT0KISkbC4hH9tHUzvVNXA6Dp9inmThV7sDAhhCiIhMUjOB+bzL2j6zFXKVzz7IGrfz9TlyQegVzrISobQ/5My015paTJzmH82lNcyBpAlkcL3hj0aEMRCNPJuzM5PT3dIHcKC1Fe5N2dbm5uXkTLoklYlNK83y5xIS4FFzs1fQaPQ2WjNnVJopTMzc1xdnbWjXNka2src46ICk+r1XLr1i1sbW0Ncse5hEUpnLpwGa8D7+DMK3zUuyM1HCQoKjoPDw+gdAPjCVFemZmZUadOHYP88SNhUUJJaZkkrRvNAPPDtHFM5PFmr5q6JGEAKpWKmjVr4ubmlm/SHSEqKisrK4MN4yJhUQJarcLOZTPopz1MFhZ4vjzX1CUJAzM3NzfI+V0hKpsKezXUwoULqVevHtbW1vj7+xc5Tsz69evx8fHB2toaX19ftm/fXuL33PLzRvre/RqAWx3ew75e/nH5hRCiMqqQYbF27VomTJhAaGgoJ06cwM/Pj8DAwELPNx88eJD+/fszfPhwTp48Sa9evejVqxdnz54t0fs+9ecMLFRarnn2wLPLGEPsihBCVAgVciBBf39/nnjiCRYsWADk9vp7eXkxduxYpk6dmq99v379SEtLY+vWrbpl7dq1o0WLFixevLjI99MNuDXVgYxqDXEL/t2gk9cIIYSxlXQgwQrXZ5GZmcnx48eZNm2abpmZmRldunQhIiKiwHUiIiKYMGGC3rLAwEA2b95cYHuNRoNGo9E9T0pKAiBO64L7q9+TnJEDGcmPuCdCCGE6ycm5n2HFPV6ocGFx+/ZtcnJydHMA53F3d+fChQsFrhMXF1dg+7i4uALbz5o1i5kzZ+Zb3mjuXzC3aSkrF0KI8iclJQUnJ6ci21W4sDCGadOm6R2JJCYmUrduXaKjo4v1TRWiJJKTk/Hy8uL69evFOh1QWVT0/a4I9T+sRkVRSElJwdOzeOPZVbiwcHV1xdzcnPj4eL3l8fHxuhur/s3Dw6NE7dVqNWp1/hvtnJycyu0Phaj4HB0dq+TPV0Xf74pQf2E1luSP3wp3NZSVlRWtW7cmLCxMt0yr1RIWFkZAQECB6wQEBOi1B9i1a1eh7YUQQuircEcWABMmTGDIkCG0adOGtm3bMm/ePNLS0hg2bBgAgwcPplatWsyaNQuAcePG0bFjRz799FN69OjBmjVrOHbsGF9//bUpd0MIISqMChkW/fr149atW4SEhBAXF0eLFi3YsWOHrhM7Ojpa7xb39u3bs2rVKqZPn84777xDw4YN2bx5M82aNSvW+6nVakJDQws8NSXEo6qqP18Vfb8rQv2GrLFC3mchhBDCuCpcn4UQQgjjk7AQQghRJAkLIYQQRZKwEEIIUSQJi2Io6XDoQhTlvffeQ6VS6T18fHxMXVaZ+P333+nZsyeenp6oVKp8Y7IpikJISAg1a9bExsaGLl26cOnSJdMUW4Ci6h86dGi+/8tu3boZrb5Zs2bxxBNP4ODggJubG7169SIyMlKvTUZGBkFBQVSvXh17e3v69u2b70blokhYFKGkw6ELUVxNmzYlNjZW99i/f7+pSyoTaWlp+Pn5sXDhwgJfnzt3Ll988QWLFy/m8OHD2NnZERgYSEZGhpErLVhR9QN069ZN7/9y9erVRqsvPDycoKAgDh06xK5du8jKyqJr166kpaXp2owfP56ff/6Z9evXEx4eTkxMDH369CnZGyniodq2basEBQXpnufk5Cienp7KrFmzTFiVqOhCQ0MVPz8/U5dhdICyadMm3XOtVqt4eHgoH3/8sW5ZYmKiolarldWrV5ugwof7d/2KoihDhgxRXnzxRZPUU5CEhAQFUMLDwxVFyf1+WlpaKuvXr9e1OX/+vAIoERERxd6uHFk8RN5w6F26dNEtK2o4dCGK69KlS3h6etKgQQMGDhxIdHS0qUsyuqtXrxIXF6f3O+bk5IS/v3+F+h3bu3cvbm5uNGrUiFGjRnHnzh2T1ZI3pYKLiwsAx48fJysrS+977OPjQ506dUr0PZaweIiHDYde2PDmQhSHv78/y5cvZ8eOHSxatIirV6/y1FNPkZKSYurSjCrv96gi/45169aN7777jrCwMObMmUN4eDjdu3cnJyfH6LVotVqCg4Pp0KGDboSKuLg4rKyscHZ21mtb0u9xhRzuQ4iKrnv37rqvmzdvjr+/P3Xr1mXdunUMHz7chJWJknr11Vd1X/v6+tK8eXMee+wx9u7dy7PPPmvUWoKCgjh79myZ9H/JkcVDlGY4dCFKw9nZmccff5zLly+buhSjyvs9qky/Yw0aNMDV1dXo/5djxoxh69at7Nmzh9q1a+uWe3h4kJmZSWJiol77kn6PJSweojTDoQtRGqmpqURFRVGzZk1Tl2JU9evXx8PDQ+93LDk5mcOHD1fY37EbN25w584do/1fKorCmDFj2LRpE7t376Z+/fp6r7du3RpLS0u973FkZCTR0dEl+h7LaagiFDUcuhClMXHiRHr27EndunWJiYkhNDQUc3Nz+vfvb+rSDC41NVXvr+yrV69y6tQpXFxcqFOnDsHBwXz44Yc0bNiQ+vXrM2PGDDw9PenVq5fpin7Aw+p3cXFh5syZ9O3bFw8PD6Kiopg8eTLe3t4EBgYapb6goCBWrVrFli1bcHBw0PVDODk5YWNjg5OTE8OHD2fChAm4uLjg6OjI2LFjCQgIoF27dsV/I0NftlUZ/d///Z9Sp04dxcrKSmnbtq1y6NAhU5ckKrh+/fopNWvWVKysrJRatWop/fr1Uy5fvmzqssrEnj17FCDfY8iQIYqi5F4+O2PGDMXd3V1Rq9XKs88+q0RGRpq26Ac8rP709HSla9euSo0aNRRLS0ulbt26yogRI5S4uDij1VdQbYCybNkyXZv79+8ro0ePVqpVq6bY2toqvXv3VmJjY0v0PjJEuRBCiCJJn4UQQogiSVgIIYQokoSFEEKIIklYCCGEKJKEhRBCiCJJWAghhCiShIUQQogiSVgIYWRDhw4tN3cnC1FcMtyHEAakUqke+npoaCjz589H7oUVFY2EhRAGFBsbq/t67dq1hISE6M2HbG9vj729vSlKE+KRyGkoIQzIw8ND93ByckKlUukts7e3z3caqlOnTowdO5bg4GCqVauGu7s7S5Ys0Q1Y6eDggLe3N7/88ovee509e5bu3btjb2+Pu7s7gwYN4vbt20beY1FVSFgIUQ6sWLECV1dXjhw5wtixYxk1ahQvv/wy7du358SJE3Tt2pVBgwaRnp4OQGJiIp07d6Zly5YcO3aMHTt2EB8fzyuvvGLiPRGVlYSFEOWAn58f06dPp2HDhkybNg1ra2tcXV0ZMWIEDRs2JCQkhDt37vDHH38AsGDBAlq2bMlHH32Ej48PLVu25Ntvv2XPnj1cvHjRxHsjKiPpsxCiHGjevLnua3Nzc6pXr46vr69uWd4c1QkJCQCcPn2aPXv2FNj/ERUVxeOPP17GFYuqRsJCiHLA0tJS77lKpdJblneVlVarBXIn5OnZsydz5szJt62qNtueMA4JCyEqoFatWrFhwwbq1auHhYX8GouyJ30WQlRAQUFB3L17l/79+3P06FGioqLYuXMnw4YNIycnx9TliUpIwkKICsjT05MDBw6Qk5ND165d8fX1JTg4GGdnZ8zM5NdaGJ5MqyqEEKJIcrKzGLRaLTExMTg4OBQ5nIMQQlQEiqKQkpKCp6dnsY5GJSyKISYmBi8vL1OXIYQQBnf9+nVq165dZDsJi2JwcHAAcr+pjo6OJq5GCCEeXXJyMl5eXrrPt6JIWBRD3qknR0dHCQshRKVS3FPrctmEEEKIIklYCCGEKJJJw+K9995DpVLpPXx8fHSvZ2RkEBQURPXq1bG3t6dv377Ex8frbSM6OpoePXpga2uLm5sbkyZNIjs7W6/N3r17adWqFWq1Gm9vb5YvX26M3RNCiErD5EcWTZs2JTY2VvfYv3+/7rXx48fz888/s379esLDw4mJiaFPnz6613NycujRoweZmZkcPHiQFStWsHz5ckJCQnRtrl69So8ePXjmmWc4deoUwcHBvP766+zcudOo+ymEEBWaYkKhoaGKn59fga8lJiYqlpaWyvr163XLzp8/rwBKRESEoiiKsn37dsXMzEyJi4vTtVm0aJHi6OioaDQaRVEUZfLkyUrTpk31tt2vXz8lMDCw2HUmJSUpgJKUlFTsdYQQojwr6eeaya+GunTpEp6enlhbWxMQEMCsWbOoU6cOx48fJysriy5duuja+vj4UKdOHSIiImjXrh0RERH4+vrqhm8GCAwMZNSoUfz555+0bNmSiIgIvW3ktQkODi60Jo1Gg0aj0T1PTk4GYOCSQ1ja2BlozyuGJ+q5MO35xqYuQwhhYiYNC39/f5YvX06jRo2IjY1l5syZPPXUU5w9e5a4uDisrKxwdnbWW8fd3Z24uDgA4uLi9IIi7/W81x7WJjk5mfv372NjY5OvrlmzZjFz5sx8y0/fSMJMnVXq/a2ITkQnMrLjY1SzszJ1KUIIEzJpWHTv3l33dfPmzfH396du3bqsW7euwA9xY5k2bRoTJkzQPc+7eWX+qy2wsy/eDSyVwRvfHwcgWyvDhwlR1Zn8NNSDnJ2defzxx7l8+TLPPfccmZmZJCYm6h1dxMfH4+HhAYCHhwdHjhzR20be1VIPtvn3FVTx8fE4OjoWGkhqtRq1Wp1v+bON3avUTXkqFcgwk0IIKAdXQz0oNTWVqKgoatasSevWrbG0tCQsLEz3emRkJNHR0QQEBAAQEBDAmTNndFNNAuzatQtHR0eaNGmia/PgNvLa5G1DFC7vvk4FSQwhqjqThsXEiRMJDw/n2rVrHDx4kN69e2Nubk7//v1xcnJi+PDhTJgwgT179nD8+HGGDRtGQEAA7dq1A6Br1640adKEQYMGcfr0aXbu3Mn06dMJCgrSHRmMHDmSK1euMHnyZC5cuMCXX37JunXrGD9+vCl3XQghKhSTnoa6ceMG/fv3586dO9SoUYMnn3ySQ4cOUaNGDQA+//xzzMzM6Nu3LxqNhsDAQL788kvd+ubm5mzdupVRo0YREBCAnZ0dQ4YM4f3339e1qV+/Ptu2bWP8+PHMnz+f2rVrs3TpUgIDA42+vxWNKu88lBxYCFHlyeRHxZCcnIyTkxNJSUlVqs/isXe2k6NVOPzOs7g7Wpu6HCGEAZX0c61c9VmI8kWmeRJC5JGwEIXKG7lYjj2FEBIWokhyNZQQQsJCFEolJ6KEEH+TsBBFktNQQggJC1E4ObAQQvxNwkIU6p87uIUQVZ2EhSiS3IojhJCwEIVSyWkoIcTfJCxEofKuhpIDCyGEhIUQQogiSViIQslpKCFEHgkLUSQ5DSWEkLAQhZIDCyFEHgkLUSjV3+ehZGwoIYSEhSiSnIYSQkhYiELJaSghRB4JC1G4vPksTFuFEKIckLAQRZLhPoQQEhaiUHIaSgiRR8JCFEmOK4QQEhaiUCq5hVsI8TcJC1GovKyQLgshhISFKAZJCyGqOgkLUSg5CSWEyCNhIYokp6GEEBIWolD/jA0lhKjqJCxEoeQ0lBAij4SFKJKchhJCSFiIQsltFkKIPBIW4iFkPgshRC4JC1EkOQ0lhJCwEIWS01BCiDwSFqJIcmQhhJCwEIXKO7CQPgshhISFKJSchhJC5JGwEEWS01BCCAkLUSiV3MMthPibhIUQQogiSViIQsnkR0KIPOUmLGbPno1KpSI4OFi3LCMjg6CgIKpXr469vT19+/YlPj5eb73o6Gh69OiBra0tbm5uTJo0iezsbL02e/fupVWrVqjVary9vVm+fLkR9qjik5NQQog85SIsjh49yldffUXz5s31lo8fP56ff/6Z9evXEx4eTkxMDH369NG9npOTQ48ePcjMzOTgwYOsWLGC5cuXExISomtz9epVevTowTPPPMOpU6cIDg7m9ddfZ+fOnUbbv4pOLp0VQqCYWEpKitKwYUNl165dSseOHZVx48YpiqIoiYmJiqWlpbJ+/Xpd2/PnzyuAEhERoSiKomzfvl0xMzNT4uLidG0WLVqkODo6KhqNRlEURZk8ebLStGlTvffs16+fEhgYWOwak5KSFEBJSkoq7W5WSO1nhSl1p2xVTkXfM3UpQggDK+nnmsmPLIKCgujRowddunTRW378+HGysrL0lvv4+FCnTh0iIiIAiIiIwNfXF3d3d12bwMBAkpOT+fPPP3Vt/r3twMBA3TYKotFoSE5O1nsIIURVZmHKN1+zZg0nTpzg6NGj+V6Li4vDysoKZ2dnveXu7u7ExcXp2jwYFHmv5732sDbJycncv38fGxubfO89a9YsZs6cWer9qmzkJJQQwmRHFtevX2fcuHGsXLkSa2trU5VRoGnTppGUlKR7XL9+3dQlmcQ/V0NJXAhR1ZksLI4fP05CQgKtWrXCwsICCwsLwsPD+eKLL7CwsMDd3Z3MzEwSExP11ouPj8fDwwMADw+PfFdH5T0vqo2jo2OBRxUAarUaR0dHvYcQQlRlJguLZ599ljNnznDq1Cndo02bNgwcOFD3taWlJWFhYbp1IiMjiY6OJiAgAICAgADOnDlDQkKCrs2uXbtwdHSkSZMmujYPbiOvTd42ROF0RxamLUMIUQ6YrM/CwcGBZs2a6S2zs7OjevXquuXDhw9nwoQJuLi44OjoyNixYwkICKBdu3YAdO3alSZNmjBo0CDmzp1LXFwc06dPJygoCLVaDcDIkSNZsGABkydP5rXXXmP37t2sW7eObdu2GXeHKyAZ7kMIkcekHdxF+fzzzzEzM6Nv375oNBoCAwP58ssvda+bm5uzdetWRo0aRUBAAHZ2dgwZMoT3339f16Z+/fps27aN8ePHM3/+fGrXrs3SpUsJDAw0xS5VSNJlIYRQKdJ7WaTk5GScnJxISkqqUv0XHT/ew1930tkwKoDWdV1MXY4QwoBK+rlm8vssRPklJ6GEEHkkLESR5NhTCCFhIQql+vtyKMkKIYSEhRBCiCJJWIhC5fVZyGkoIYSEhSic9HALIf4mYSGKJFdXCyEkLEShdKehTFqFEKI8kLAQQghRJAkLUSjdpbNyaCFElSdhIQr1z2koSQshqjoJCyGEEEUqdVhkZ2fz22+/8dVXX5GSkgJATEwMqampBitOmJZKeriFEH8r1RDlf/31F926dSM6OhqNRsNzzz2Hg4MDc+bMQaPRsHjxYkPXKUxA5rMQQuQp1ZHFuHHjaNOmDffu3dObmrR37975ZqUTFZ8cWAghSnVksW/fPg4ePIiVlZXe8nr16nHz5k2DFCZMTzetqqSFEFVeqY4stFotOTk5+ZbfuHEDBweHRy5KCCFE+VKqsOjatSvz5s3TPVepVKSmphIaGsrzzz9vqNpEOSGXzgohSnUa6pNPPqFbt240adKEjIwMBgwYwKVLl3B1dWX16tWGrlGYiNyUJ4TIU6qw8PLy4vTp06xdu5bTp0+TmprK8OHDGThwoF6HtxBCiMqhxGGRlZWFj48PW7duZeDAgQwcOLAs6hLlgNxmIYTIU+I+C0tLSzIyMsqiFlHOqOQ2CyHE30rVwR0UFMScOXPIzs42dD2iHJL5LIQQpeqzOHr0KGFhYfz666/4+vpiZ2en9/rGjRsNUpwwLd19FqYtQwhRDpQqLJydnenbt6+haxFCCFFOlSosli1bZug6RDmkGxtKDi2EqPJkiHJRqH9OQ0laCFHVlerIon79+robtgpy5cqVUhckhBCi/ClVWAQHB+s9z8rK4uTJk+zYsYNJkyYZoi5RDujus5ADCyGqvFKFxbhx4wpcvnDhQo4dO/ZIBQkhhCh/DNpn0b17dzZs2GDITQpTkrGhhBB/M2hY/Pjjj7i4uBhyk8KEZLgPIUSeUp2GatmypV4Ht6IoxMXFcevWLb788kuDFSeEEKJ8KFVY9OrVS++5mZkZNWrUoFOnTvj4+BiiLlEO/DNTnhxbCFHVlSosQkNDDV2HKIfkNJQQIk+p+ixOnDjBmTNndM+3bNlCr169eOedd8jMzDRYcUIIIcqHUoXFm2++ycWLF4HcG/D69euHra0t69evZ/LkyQYtUJiOzJQnhMhTqrC4ePEiLVq0AGD9+vV07NiRVatWsXz5crl0VgghKqFShYWiKGi1WgB+++03nn/+eSB3utXbt28brjphUg9c72bCKoQQ5UGpwqJNmzZ8+OGHfP/994SHh9OjRw8Arl69iru7u0ELFKbzz9VQpq1DCGF6pQqLefPmceLECcaMGcO7776Lt7c3kHtTXvv27Yu9nUWLFtG8eXMcHR1xdHQkICCAX375Rfd6RkYGQUFBVK9eHXt7e/r27Ut8fLzeNqKjo+nRowe2tra4ubkxadKkfDP47d27l1atWqFWq/H29mb58uWl2W0hhKiySnXpbPPmzfWuhsrz8ccfY25uXuzt1K5dm9mzZ9OwYUMURWHFihW8+OKLnDx5kqZNmzJ+/Hi2bdvG+vXrcXJyYsyYMfTp04cDBw4AkJOTQ48ePfDw8ODgwYPExsYyePBgLC0t+eijj4Dco50ePXowcuRIVq5cSVhYGK+//jo1a9YkMDCwNLtfZeTNZyEHFkIIlFKIjo5Wrl+/rnt++PBhZdy4ccpXX31Vms3pqVatmrJ06VIlMTFRsbS0VNavX6977fz58wqgREREKIqiKNu3b1fMzMyUuLg4XZtFixYpjo6OikajURRFUSZPnqw0bdpU7z369eunBAYGFrumpKQkBVCSkpIeZdcqnJcXHVTqTtmqbD0dY+pShBAGVtLPtVKdhhowYAB79uwBIC4ujueee44jR47w7rvv8v7775cqtHJyclizZg1paWkEBARw/PhxsrKy6NKli66Nj48PderUISIiAoCIiAh8fX31+kkCAwNJTk7mzz//1LV5cBt5bfK2URCNRkNycrLeo0oqfMoSIUQVU6qwOHv2LG3btgVg3bp1NGvWjIMHD7Jy5coS9wecOXMGe3t71Go1I0eOZNOmTTRp0oS4uDisrKxwdnbWa+/u7k5cXByQG1T/7lDPe15Um+TkZO7fv19gTbNmzcLJyUn38PLyKtE+VRb/3MEtJ6KEqOpKFRZZWVmo1Wog99LZF154Acj9yz82NrZE22rUqBGnTp3i8OHDjBo1iiFDhnDu3LnSlGUw06ZNIykpSfe4fv26SesRQghTK1VYNG3alMWLF7Nv3z527dpFt27dAIiJiaF69eol2paVlRXe3t60bt2aWbNm4efnx/z58/Hw8CAzM5PExES99vHx8Xh4eADg4eGR7+qovOdFtXF0dMTGxqbAmtRqte4KrbxHVSSXzgoh8pQqLObMmcNXX31Fp06d6N+/P35+fgD89NNPutNTpaXVatFoNLRu3RpLS0vCwsJ0r0VGRhIdHU1AQAAAAQEBnDlzhoSEBF2bXbt24ejoSJMmTXRtHtxGXpu8bYjCydVQQog8pbp0tlOnTty+fZvk5GSqVaumW/7GG29ga2tb7O1MmzaN7t27U6dOHVJSUli1ahV79+5l586dODk5MXz4cCZMmICLiwuOjo6MHTuWgIAA2rVrB0DXrl1p0qQJgwYNYu7cucTFxTF9+nSCgoJ0p8lGjhzJggULmDx5Mq+99hq7d+9m3bp1bNu2rTS7LoQQVVKpwgJyh/w4fvw4UVFRDBgwAAcHB6ysrEoUFgkJCQwePJjY2FicnJxo3rw5O3fu5LnnngPg888/x8zMjL59+6LRaAgMDNSbXMnc3JytW7cyatQoAgICsLOzY8iQIXpXZNWvX59t27Yxfvx45s+fT+3atVm6dKncY1EMMp+FECKPSinFJ8Fff/1Ft27diI6ORqPRcPHiRRo0aMC4cePQaDQsXry4LGo1meTkZJycnEhKSqpS/RcDlhziYNQd5r/aghdb1DJ1OUIIAyrp51qp+izGjRtHmzZtuHfvnl4nce/evfP1D4iKSyX3WQgh/laq01D79u3j4MGDWFlZ6S2vV68eN2/eNEhhwvR0HdxyFkqIKq9URxZarZacnJx8y2/cuIGDg8MjFyWEEKJ8KVVYdO3alXnz5umeq1QqUlNTCQ0N1c1tISo+XQe3XDwrRJVXqtNQn3zyCd26daNJkyZkZGQwYMAALl26hKurK6tXrzZ0jcLE5DSUEKJUYeHl5cXp06dZu3Ytp0+fJjU1leHDhzNw4MBC74oWQghRcZU4LLKysvDx8WHr1q0MHDiQgQMHlkVdohxQqaSDWwiRq8R9FpaWlmRkZJRFLaKckqwQQpSqgzsoKIg5c+bkm75UVC5ym4UQIk+p+iyOHj1KWFgYv/76K76+vtjZ2em9vnHjRoMUJ0xLhvsQQuQpVVg4OzvTt29fQ9cihBCinCpRWGi1Wj7++GMuXrxIZmYmnTt35r333pMroCqpf2bKE0JUdSXqs/jf//7HO++8g729PbVq1eKLL74gKCiorGoT5YWkhRBVXonC4rvvvuPLL79k586dbN68mZ9//pmVK1ei1WrLqj5hQioZSVAI8bcShUV0dLTecB5dunRBpVIRExNj8MKE6f1zGkoOLYSo6koUFtnZ2VhbW+sts7S0JCsry6BFifJFLoYSQpSog1tRFIYOHaqbshQgIyODkSNH6l0+K5fOVg5yFkoIkadEYTFkyJB8y/773/8arBhR3uSmxa5z8cQmyV37xWVlYUbfVrXxcLIuurEQFUSJwmLZsmVlVYcoh2ytzAEIu5BA2IUEE1dTsVy/m87svs1NXYYQBlOqm/JE1TCuS0Nc7dVk5cjVbsUVGZfCkWt3Sc6QfjxRuUhYiEI9VsOekJ5NTF1GhfJ9xDWOXLsrFwWISqdUAwkKIQqWd2+KVtJCVDISFkIY0D+DL5q2DiEMTcJCCANS/X0FmWSFqGwkLIQwIDmyEJWVhIUQBqQbIkXSQlQyEhZCGJCZSk5DicpJwkIIQ5LZBUUlJWEhhAHJhFGispKwEMKA/rnPwsSFCGFgEhZCGJCZnIYSlZSEhRAGJMO6i8pKwkIIA9LdlCcHFqKSkbAQwoDyjixkbChR2UhYCGFAeR3ckhWispGwEMKA/rl0VtJCVC4SFkIYkIwNJSorCQshDEg6uEVlJWEhhAHp7rOQ01CikpGwEMKA5DSUqKxMGhazZs3iiSeewMHBATc3N3r16kVkZKRem4yMDIKCgqhevTr29vb07duX+Ph4vTbR0dH06NEDW1tb3NzcmDRpEtnZ2Xpt9u7dS6tWrVCr1Xh7e7N8+fKy3j1RJcmos6JyMmlYhIeHExQUxKFDh9i1axdZWVl07dqVtLQ0XZvx48fz888/s379esLDw4mJiaFPnz6613NycujRoweZmZkcPHiQFStWsHz5ckJCQnRtrl69So8ePXjmmWc4deoUwcHBvP766+zcudOo+ysqP5UM9yEqK6UcSUhIUAAlPDxcURRFSUxMVCwtLZX169fr2pw/f14BlIiICEVRFGX79u2KmZmZEhcXp2uzaNEixdHRUdFoNIqiKMrkyZOVpk2b6r1Xv379lMDAwGLVlZSUpABKUlLSI+2fqPx+/TNOqTtlq/Ligv2mLkWIhyrp51q56rNISkoCwMXFBYDjx4+TlZVFly5ddG18fHyoU6cOERERAERERODr64u7u7uuTWBgIMnJyfz555+6Ng9uI69N3jb+TaPRkJycrPcQojhkiHJRWZWbsNBqtQQHB9OhQweaNWsGQFxcHFZWVjg7O+u1dXd3Jy4uTtfmwaDIez3vtYe1SU5O5v79+/lqmTVrFk5OTrqHl5eXQfZRVH6qf+ZVNWkdQhhauQmLoKAgzp49y5o1a0xdCtOmTSMpKUn3uH79uqlLEhWErs/CtGUIYXAWpi4AYMyYMWzdupXff/+d2rVr65Z7eHiQmZlJYmKi3tFFfHw8Hh4eujZHjhzR217e1VIPtvn3FVTx8fE4OjpiY2OTrx61Wo1arTbIvomq5Z/JjyQuROVi0iMLRVEYM2YMmzZtYvfu3dSvX1/v9datW2NpaUlYWJhuWWRkJNHR0QQEBAAQEBDAmTNnSEhI0LXZtWsXjo6ONGnSRNfmwW3ktcnbhhCGImehRGVl0iOLoKAgVq1axZYtW3BwcND1MTg5OWFjY4OTkxPDhw9nwoQJuLi44OjoyNixYwkICKBdu3YAdO3alSZNmjBo0CDmzp1LXFwc06dPJygoSHd0MHLkSBYsWMDkyZN57bXX2L17N+vWrWPbtm0m23dROcmos6LSKtuLsx6O3FO7+R7Lli3Ttbl//74yevRopVq1aoqtra3Su3dvJTY2Vm87165dU7p3767Y2Ngorq6uyttvv61kZWXptdmzZ4/SokULxcrKSmnQoIHeexRFLp0VxRUemaDUnbJV6Tbvd1OXIsRDlfRzTaUo8jdQUZKTk3FyciIpKQlHR0dTlyPKsf2XbvPfbw7j4+HAjuCnTV2OEIUq6edaubkaSojKQMaGEpWVhIUQBiSTH4nKSsJCCEOSIwtRSUlYCGFAZnKfhaikJCyEMCAZG0pUVhIWQhiQSsb7EJWUhIUQBiRZISorCQshDChvDm7psxCVjYSFEAYlw32IyknCQggD+uc0lKSFqFwkLIQwIBl1VlRWEhZCGJCZjDorKikJCyEM6J+xoSQtROUiYSGEAanyOrhNXIcQhiZhIYQByaizorKSsBDCgFRyn4WopCQshDAgOQ0lKisJCyEMSE5DicpKwkIIA8oLCzm2EJWNhIUQBvTPfBYmLkQIA5OwEMKA/rmDW9JCVC4SFkIYkAxRLiorCQshDEqG+xCVk4SFEAYk81mIykrCQggDkmlVRWUlYSGEAek6uE1ahRCGJ2EhhAHJqLOispKwEMKA5D4LUVlJWAhRBmRaVVHZSFgIYUAyNpSorCxMXYAQlUne1VCabC1Pz91j4mpEVdCghh3Lh7Ut8/eRsBDCgKrbWeFkY0nS/Syi76abuhxRBdhamRvlfSQshDAga0tz9k7sxNU7aaYuRVQRNpYSFkJUSNXsrKhmZ2XqMoQwKOngFkIIUSQJCyGEEEWSsBBCCFEkCQshhBBFkrAQQghRJAkLIYQQRZJLZ4shbwTR5ORkE1cihBCGkfd5VtwRkiUsiiElJQUALy8vE1cihBCGlZKSgpOTU5HtVIoMvF8krVZLTEwMDg4O/8yEJkQhkpOT8fLy4vr16zg6Opq6HKOrDPtfUfbhUepUFIWUlBQ8PT0xMyu6R0KOLIrBzMyM2rVrm7oMUcE4OjqW6w+aslYZ9r+i7ENp6yzOEUUe6eAWQghRJAkLIYQQRZKwEMLA1Go1oaGhqNVqU5diEpVh/yvKPhizTungFkIIUSQ5shBCCFEkCQshhBBFkrAQQghRJAkLIYQQRZKwEMIA3nvvPVQqld7Dx8fH1GWVqd9//52ePXvi6emJSqVi8+bNeq8rikJISAg1a9bExsaGLl26cOnSJdMUW4ii9mHo0KH5/l+7detm1BpnzZrFE088gYODA25ubvTq1YvIyEi9NhkZGQQFBVG9enXs7e3p27cv8fHxBq1DwkIIA2natCmxsbG6x/79+01dUplKS0vDz8+PhQsXFvj63Llz+eKLL1i8eDGHDx/Gzs6OwMBAMjIyjFxp4YraB4Bu3brp/b+uXr3aiBVCeHg4QUFBHDp0iF27dpGVlUXXrl1JS0vTtRk/fjw///wz69evJzw8nJiYGPr06WPYQhQhxCMLDQ1V/Pz8TF2GyQDKpk2bdM+1Wq3i4eGhfPzxx7pliYmJilqtVlavXm2CCov2731QFEUZMmSI8uKLL5qknsIkJCQogBIeHq4oSu731dLSUlm/fr2uzfnz5xVAiYiIMNj7ypGFEAZy6dIlPD09adCgAQMHDiQ6OtrUJZnM1atXiYuLo0uXLrplTk5O+Pv7ExERYcLKSm7v3r24ubnRqFEjRo0axZ07d0xaT1JSEgAuLi4AHD9+nKysLL3vtY+PD3Xq1DHo91rCQggD8Pf3Z/ny5ezYsYNFixZx9epVnnrqKd3w9lVNXFwcAO7u7nrL3d3dda9VBN26deO7774jLCyMOXPmEB4eTvfu3cnJyTFJPVqtluDgYDp06ECzZs2A3O+1lZUVzs7Oem0N/b2WUWeFMIDu3bvrvm7evDn+/v7UrVuXdevWMXz4cBNWJh7Fq6++qvva19eX5s2b89hjj7F3716effZZo9cTFBTE2bNnTdIfJkcWQpQBZ2dnHn/8cS5fvmzqUkzCw8MDIN8VOfHx8brXKqIGDRrg6upqkv/XMWPGsHXrVvbs2aM3ZYKHhweZmZkkJibqtTf091rCQogykJqaSlRUFDVr1jR1KSZRv359PDw8CAsL0y1LTk7m8OHDBAQEmLCyR3Pjxg3u3Llj1P9XRVEYM2YMmzZtYvfu3dSvX1/v9datW2Npaan3vY6MjCQ6Otqg32s5DSWEAUycOJGePXtSt25dYmJiCA0NxdzcnP79+5u6tDKTmpqq9xf21atXOXXqFC4uLtSpU4fg4GA+/PBDGjZsSP369ZkxYwaenp706tXLdEX/y8P2wcXFhZkzZ9K3b188PDyIiopi8uTJeHt7ExgYaLQag4KCWLVqFVu2bMHBwUHXD+Hk5ISNjQ1OTk4MHz6cCRMm4OLigqOjI2PHjiUgIIB27doZrhCDXVclRBXWr18/pWbNmoqVlZVSq1YtpV+/fsrly5dNXVaZ2rNnjwLkewwZMkRRlNzLZ2fMmKG4u7srarVaefbZZ5XIyEjTFv0vD9uH9PR0pWvXrkqNGjUUS0tLpW7dusqIESOUuLg4o9ZYUH2AsmzZMl2b+/fvK6NHj1aqVaum2NraKr1791ZiY2MNWocMUS6EEKJI0mchhBCiSBIWQgghiiRhIYQQokgSFkIIIYokYSGEEKJIEhZCCCGKJGEhhBCiSBIWQhjZ0KFDy9VdzEIUhwz3IYQBqVSqh74eGhrK/PnzkXthRUUjYSGEAcXGxuq+Xrt2LSEhIXrzJdvb22Nvb2+K0oR4JHIaSggD8vDw0D2cnJxQqVR6y+zt7fOdhurUqRNjx44lODiYatWq4e7uzpIlS0hLS2PYsGE4ODjg7e3NL7/8ovdeZ8+epXv37tjb2+Pu7s6gQYO4ffu2kfdYVBUSFkKUAytWrMDV1ZUjR44wduxYRo0axcsvv0z79u05ceIEXbt2ZdCgQaSnpwOQmJhI586dadmyJceOHWPHjh3Ex8fzyiuvmHhPRGUlYSFEOeDn58f06dNp2LAh06ZNw9raGldXV0aMGEHDhg0JCQnhzp07/PHHHwAsWLCAli1b8tFHH+Hj40PLli359ttv2bNnDxcvXjTx3ojKSPoshCgHmjdvrvva3Nyc6tWr4+vrq1uWN5d1QkICAKdPn2bPnj0F9n9ERUXx+OOPl3HFoqqRsBCiHLC0tNR7rlKp9JblXWWl1WqB3El7evbsyZw5c/Jtq6rOzifKloSFEBVQq1at2LBhA/Xq1cPCQn6NRdmTPgshKqCgoCDu3r1L//79OXr0KFFRUezcuZNhw4aRk5Nj6vJEJSRhIUQF5OnpyYEDB8jJyaFr1674+voSHByMs7MzZmbyay0MT6ZVFUIIUST5E0QIIUSRJCyEEEIUScJCCCFEkSQshBBCFEnCQgghRJEkLIQQQhRJwkIIIUSRJCyEEEIUScJCCCFEkSQshBBCFEnCQgghRJEkLIQQQhTp/wGppS1yHpabsQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "pressure_initial = 12_000.0\n", "pressure_fracface = 5_000.0\n", "\n", "n_times = 1_000\n", "t_end = 20\n", "time = np.linspace(0, np.sqrt(t_end), n_times) ** 2\n", "pressure_v_time = np.full(n_times, pressure_fracface)\n", "pressure_v_time[n_times // 4 : n_times // 2] /= 2.0\n", "pressure_v_time[n_times // 2 :] /= 4.0\n", "\n", "pvt_gas = pd.read_csv(\n", " \"https://raw.githubusercontent.com/frank1010111/bluebonnet/main/tests/data/pvt_gas.csv\"\n", ").rename(\n", " columns={\n", " \"P\": \"pressure\",\n", " \"Z-Factor\": \"z-factor\",\n", " \"Cg\": \"compressibility\",\n", " \"Viscosity\": \"viscosity\",\n", " }\n", ")\n", "flow_properties = FlowProperties(pvt_gas, pressure_initial)\n", "pseudopressure_fracface = flow_properties.m_scaled_func(pressure_fracface)\n", "pseudopressure_initial = flow_properties.m_scaled_func(pressure_initial)\n", "\n", "print(\"mf =\", flow_properties.m_scaled_func(pressure_fracface))\n", "\n", "res = SinglePhaseReservoir(100, pressure_fracface, pressure_initial, flow_properties)\n", "\n", "%time res.simulate(time, pressure_v_time)\n", "\n", "recovery = res.recovery_factor()\n", "\n", "\n", "def resource_left(pseudopressure, pvt):\n", " density = interpolate.interp1d(pvt.pvt_props[\"m-scaled\"], pvt.pvt_props[\"Density\"])\n", " print(max(pvt.pvt_props[\"m-scaled\"]))\n", " p = np.minimum(pseudopressure, max(pvt.pvt_props[\"m-scaled\"]))\n", " return (density(p)).sum(axis=1) / len(pseudopressure[0])\n", "\n", "\n", "density_interp = interpolate.interp1d(\n", " flow_properties.pvt_props[\"m-scaled\"], flow_properties.pvt_props[\"Density\"]\n", ")\n", "remaining_gas = (resource_left(res.pseudopressure, flow_properties)) / (\n", " density_interp(pseudopressure_initial)\n", ")\n", "print(\n", " f\"{pseudopressure_fracface=}\",\n", " res.pseudopressure[-1][0],\n", " f\"Deviation is {(recovery[-1] - 1 + remaining_gas[-1]) / recovery[-1] * 100:3.2g}%\",\n", ")\n", "fig, (ax1, ax2) = plt.subplots(2, 1)\n", "fig.set_size_inches(4, 6)\n", "ax1.plot(time, 1 - remaining_gas, label=\"From resource left\")\n", "ax1.plot(time, recovery, \"--\", label=\"From frac face flow\")\n", "ax1.legend()\n", "ax1.set(\n", " xlabel=\"Time\",\n", " ylabel=\"Recovered gas\",\n", " ylim=(0, None),\n", " xscale=\"squareroot\",\n", " xlim=(0, None),\n", ")\n", "\n", "ax2.plot(time, pressure_v_time, label=\"Pressure\")\n", "ax2.set(xlabel=\"Time\", xscale=\"squareroot\", ylabel=\"Pressure\");" ] }, { "attachments": {}, "cell_type": "markdown", "id": "c26f8388-83d9-4e8b-9d0b-f5ea95e77a60", "metadata": {}, "source": [ "## Build PVT data\n", "\n", "There are two ways to create the necessary pressure-volume-temperature data. \n", "1. Use the well data and an appropriate equation of state\n", "2. Pull data from a pre-made csv file. \n", "\n", "We will show both with the next code cell." ] }, { "cell_type": "code", "execution_count": 3, "id": "fde6e2b6-5fef-4c97-acdc-d59a58ae39b2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " temperature pressure Density z-factor compressibility viscosity \\\n", "0 285.21375 10.0 0.021024 0.999592 0.100043 0.014638 \n", "1 285.21375 20.0 0.042065 0.999181 0.050043 0.014638 \n", "2 285.21375 30.0 0.063123 0.998775 0.033376 0.014639 \n", "3 285.21375 40.0 0.084198 0.998370 0.025043 0.014639 \n", "4 285.21375 50.0 0.105290 0.997967 0.020042 0.014640 \n", "\n", " pseudopressure \n", "0 0.000000 \n", "1 20508.736605 \n", "2 54702.013998 \n", "3 102589.697359 \n", "4 164181.289746 \n" ] } ], "source": [ "# #########################\n", "# Haynesville shale gas play\n", "# Using the well data and bluebonnet.fluids\n", "# #########################\n", "field_values = {\n", " \"N2\": 0.0,\n", " \"H2S\": 0.0,\n", " \"CO2\": 0.0002,\n", " \"Gas Specific Gravity\": 0.58,\n", " \"Reservoir Temperature (deg F)\": 285.21375,\n", "}\n", "gas_dryness = \"wet gas\"\n", "pvt_gas = build_pvt_gas(field_values, gas_dryness)\n", "\n", "\n", "# #########################\n", "# Or, you could using pre-gathered pvt data\n", "# #########################\n", "# pvt_gas = pd.read_csv(\n", "# \"https://raw.githubusercontent.com/frank1010111/bluebonnet/main/tests/data/pvt_gas_HAYNESVILLE%20SHALE_20.csv\",\n", "# index_col=0\n", "# )\n", "\n", "print(pvt_gas.head())" ] }, { "attachments": {}, "cell_type": "markdown", "id": "f7602f6a", "metadata": {}, "source": [ "## Fit wells" ] }, { "cell_type": "code", "execution_count": 4, "id": "141d42a8", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/html": [ "

Fit Result

Fit Statistics
fitting methodNelder-Mead
# function evals41
# data points416
# variables3
chi-square 52370272.2
reduced chi-square 126804.533
Akaike info crit. 4891.15639
Bayesian info crit. 4903.24845
Parameters
namevalueinitial valueminmaxvary
tau 829.929628830 30.0000000 830.000000True
M 52545.088913864.498560000016 13850.5730 100000.000True
p_initial 10403.883010000.0 9539.17738 12000.0000True
" ], "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "well_number = \"20\"\n", "prod_file = f\"https://raw.githubusercontent.com/frank1010111/bluebonnet/main/data/dataset_1_well_{well_number}.csv\"\n", "# This file contains pressure and production data.\n", "# Rows 1 and 2 have information about units, so skip\n", "prod = pd.read_csv(prod_file, skiprows=[1, 2]).rename(\n", " columns={\n", " \"Time (Days)\": \"Days\",\n", " \"Gas Volume (MMscf)\": \"Gas\",\n", " \"Calculated Sandface Pressure (psi(a))\": \"Pressure\",\n", " }\n", ")[[\"Days\", \"Gas\", \"Pressure\"]]\n", "# This file contains an estimate of initial pressure. That's all I need it for here.\n", "well_file = f\"https://media.githubusercontent.com/media/frank1010111/bluebonnet/main/data/WellData_{well_number}.csv\"\n", "pressure_initial = float(\n", " pd.read_csv(well_file).set_index(\"Field\").loc[\"Initial Pressure Estimate (psi)\"].iloc[0]\n", ")\n", "\n", "result = fit_production_pressure(\n", " prod,\n", " pvt_gas,\n", " pressure_initial,\n", " n_iter=40,\n", " pressure_imax=12000,\n", " filter_window_size=30,\n", " inplace_max=100000,\n", ")\n", "result" ] }, { "attachments": {}, "cell_type": "markdown", "id": "bc117055-5b18-42ed-b273-3b15d1f63085", "metadata": {}, "source": [ "## Plot the result" ] }, { "cell_type": "code", "execution_count": 5, "id": "e6f60a36-2958-4f44-b438-d99e022ede5a", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdkAAAINCAYAAAByAAWmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACc7klEQVR4nOzdeVhU5dvA8e8M+77IrgiouO+iiOaOoqVpWZr5ulaWqVlkqZVbVlpZmelPy3KrzLTSzNxxV9x3cUUUF0AF2WWbOe8fo6MTqMCwyv25rrnknPPMc+4Zxrk5z3kWlaIoCkIIIYQocurSDkAIIYR4UkmSFUIIIYqJJFkhhBCimEiSFUIIIYqJJFkhhBCimEiSFUIIIYqJJFkhhBCimEiSFUIIIYqJaWkHUB5otVquX7+OnZ0dKpWqtMMRQghRShRFISUlBS8vL9Tqx1+nSpLNh+vXr+Pt7V3aYQghhCgjrly5QpUqVR5bTpJsPtjZ2QG6N9Xe3r6UoxFCCFFakpOT8fb21ueFx5Ekmw/3mojt7e0lyQohhMj3rUPp+CSEEEIUE0myQgghRDGRJCuEEEIUE7knW0QURSEnJweNRlPaoQghRL6ZmJhgamoqwxOLiSTZIpCVlUVMTAzp6emlHYoQQhSYtbU1np6emJubl3YoTxxJskbSarVERUVhYmKCl5cX5ubm8hehEKJcUBSFrKwsbt68SVRUFP7+/vmaYOFJlK3RolapMFEX7fe3JFkjZWVlodVq8fb2xtraurTDEUKIArGyssLMzIzLly+TlZWFpaVlaYdU4g5dTuCDv07ycmBVBrXyLdK6K+afLMWgov71J4Qo/yry99f+qAR6zw3nbFwKP+2KIkejLdL65UpWCCFEhRXg40RzXyequdgyrlttTE2K9g8OSbJCCCEqjMvxaXwbdp6pPetjY2GKWq3il1cDsTA1KZbzVdw2AlGiBg8eTK9evYr9PL6+vsycObPYzyOEKF+yNVrmbL1Al2928Nfha3wbdl5/rLgSLEiSrdAGDx6MSqVCpVJhbm5OjRo1+Pjjj8nJySnt0B5r0aJFODo65tp/4MABhg0bVuLxtG/fnrfffrvEz5sfqampjBw5kipVqmBlZUXdunWZN2+e/nhCQgKjRo2iVq1aWFlZUbVqVd566y2SkpIM6gkLC6NVq1bY2dnh4eHB2LFjH/tZiYyM5LnnnsPV1RV7e3v69OlDXFycQZlnn32WqlWrYmlpiaenJwMGDOD69esFeo2TJ09GpVLRtWvXXMe+/PJLVCoV7du3L1CdBXXv/9KDj2XLlumP//XXX3Tu3Fn/XgQFBbFhw4Y8X8eDj9q1a+d5PkVR6NatGyqVilWrVhUolormxNUknp29my83nCUzR0sbfxdeblG1RM4tSbaC69q1KzExMZw/f553332XyZMn8+WXX+ZZNisrq4SjKzhXV1fp5f0foaGhrF+/nl9++YXTp0/z9ttvM3LkSFavXg3olnK8fv06M2bM4OTJkyxatIj169fzyiuv6Os4duwYTz/9NF27duXIkSP8/vvvrF69mnHjxj30vGlpaXTp0gWVSsWWLVvYvXs3WVlZ9OjRA632fueSDh06sHz5cs6ePcuff/5JZGQkL7zwQoFfp6enJ1u3buXq1asG+xcsWEDVqiXzhbpw4UJiYmL0jwdbb3bs2EHnzp1Zu3Ythw4dokOHDvTo0YMjR44Y1FGvXj2DOnbt2pXnuWbOnPnI4YKPiqWiyMjWMG3daXrO2cXpmGQcrc34uk8jlgxtga+LTckEoYjHSkpKUgAlKSkp17E7d+4oERERyp07d3IdS8vMfujjTlZOkZctqEGDBik9e/Y02Ne5c2elZcuWBsc/+eQTxdPTU/H19VUURVGOHz+udOjQQbG0tFScnZ2V1157TUlJSdHXkZOTo7zzzjuKg4OD4uzsrLz33nvKwIEDDc7l4+OjfPPNNwbnbtSokTJp0iT99u3bt5Vhw4Ypbm5uioWFhVKvXj3ln3/+UbZu3aoABo97z/tvvZcvX1aeffZZxcbGRrGzs1NefPFFJTY2Vn980qRJSqNGjZQlS5YoPj4+ir29vdK3b18lOTm5QO/jf+OJiopScnJylKFDhyq+vr6KpaWlUrNmTWXmzJkGz23Xrp0yevRog309e/ZUBg0alO/zP069evWUjz/+2GBf06ZNlQ8//PChz1m+fLlibm6uZGfrPlfjx49XAgICDMqsXr1asbS0fOh7tWHDBkWtVhv8v0lMTFRUKpWyadOmh57777//VlQqlZKVlfXY13bPvd9j9+7dlU8++US/f/fu3YqLi4syfPhwpV27dvmurzAAZeXKlQV6Tt26dZUpU6bot++9jsc5cuSIUrlyZSUmJibP8xY0lkd9j5VnE1edUHzGrlF8xq5RRi49rNxMyTC6zkflg7xIx6diVHfihoce61DLlYVDWui3m03dzJ3svKdkDPRz5vfXg/TbT32+lYS03FeVl6Y/Y0S0OlZWVsTHx+u3w8LCsLe3Z9OmTYDu6iQkJISgoCAOHDjAjRs3ePXVVxk5ciSLFi0C4KuvvmLRokUsWLCAOnXq8NVXX7Fy5Uo6duyY7zi0Wi3dunUjJSWFX375herVqxMREYGJiQmtWrVi5syZTJw4kbNnzwJga2ubZx09e/bE1taW7du3k5OTw4gRI+jbty/btm3Tl4uMjGTVqlWsWbOG27dv06dPH6ZPn86nn34K6JqmhwwZgqIoecb67bffcu7cOerXr8/HH38M6K6otVotVapUYcWKFVSqVIk9e/YwbNgwPD096dOnT77fizfeeINffvnlkWVSU1MfeqxVq1asXr2aoUOH4uXlxbZt2zh37hzffPPNQ5+TlJSEvb09pqa6r4jMzMxc4yetrKzIyMjg0KFDeTbFZmZmolKpsLCw0O+ztLRErVaza9cugoODcz0nISGBX3/9lVatWmFmZvbI15yXoUOH8v777/Phhx8CuqvY/v375yr366+/8vrrrz+yrnXr1tGmTZsCnX/EiBG8+uqrVKtWjTfeeIMhQ4Y89GpTq9WSkpKCs7Ozwf7z58/j5eWFpaUlQUFBTJs2zeBKPD09nZdffpk5c+bg4eFRJLE8qUZ0qEH4xXjeC6lN57rujy6s1cKFTeBaG5x8iiwGSbIC0N3fCQsLY8OGDYwaNUq/38bGhh9//FE/3dr8+fPJyMhgyZIl2Njomltmz55Njx49+Pzzz3F3d2fmzJmMHz+e559/HoB58+bluvf0OJs3b2b//v2cPn2amjVrAlCtWjX9cQcHB1Qq1SO/ZMLCwjhx4gRRUVF4e3sDsGTJEurVq8eBAwdo3rw5oPuyW7RokX4R5gEDBhAWFqZPsg4ODtSqVeuh53FwcMDc3Bxra2uDeExMTJgyZYp+28/Pj/DwcJYvX16gJPvxxx8zZsyYfJf/r++++45hw4ZRpUoVTE1NUavVzJ8/n7Zt2+ZZ/tatW0ydOtXg3nZISAgzZ87kt99+o0+fPsTGxur/oIiJicmznpYtW2JjY8PYsWP57LPPUBSFcePGodFocj1n7NixzJ49m/T0dFq2bMmaNWsK9Vq7d+/OG2+8wY4dO2jWrBnLly9n165dLFiwwKDcs88+S2Bg4CPrqly5coHO/fHHH9OxY0esra3ZuHEjb775Jqmpqbz11lt5lp8xYwapqakGn4XAwEAWLVpErVq1iImJYcqUKbRp04aTJ0/qP5/vvPMOrVq1omfPnkUWy5Nic0Qcey/G81H3ugC42Vuy4e22j/7jIj0Bjv0G++fD7Sho+SZ0nVZkMUmSLUYRH4c89Jj6P7/0QxNy/1X/sLK7xnYwLrAHrFmzBltbW7Kzs9Fqtbz88stMnjxZf7xBgwYG85mePn2aRo0a6RMsQOvWrdFqtZw9exZLS0tiYmIMvsBMTU0JCAh46JVgXo4ePUqVKlX0CbYwTp8+jbe3tz7BAtStWxdHR0dOnz6tT7K+vr76LzDQ3du7ceOGfvu5557jueeeK1QMc+bMYcGCBURHR3Pnzh2ysrJo3Lhxgepwc3PDzc2tUOcHXZLdu3cvq1evxsfHhx07djBixAi8vLxyXU0mJyfzzDPPULduXYPPQZcuXfjyyy954403GDBgABYWFkyYMIGdO3c+dCIDV1dXVqxYwfDhw5k1axZqtZp+/frRtGnTXM957733eOWVV7h8+TJTpkxh4MCBrFmzpsBXXmZmZvzf//0fCxcu5OLFi9SsWZOGDRvmKmdnZ2fwOy8KEyZM0P/cpEkT0tLS+PLLL/NMbEuXLmXKlCn8/fffBr/bbt266X9u2LAhgYGB+Pj4sHz5cl555RVWr17Nli1bct3HNSaWJ0F8aiZT/olg9TFdh7k2NV1pV9MVeMji6ooCV/bDwQVwaiVoMnX7LR10jyIkSbYYWZvn/+0trrKP06FDB+bOnYu5uTleXl765sF7HkymRUmtVudKutnZ2fqfraysiuW8eflvs6RKpTLomFNYy5YtY8yYMXz11VcEBQVhZ2fHl19+yb59+/RlHvc+gHHNxXfu3OGDDz5g5cqVPPOM7nZCw4YNOXr0KDNmzDBIsikpKXTt2hU7OztWrlyZ630JDQ3lnXfeISYmBicnJy5dusT48eMNWhj+q0uXLkRGRnLr1i1MTU1xdHTEw8Mj13NcXFxwcXGhZs2a1KlTB29vb/bu3UtQUNBDan64oUOHEhgYyMmTJxk6dGieZYqrufhBgYGBTJ06lczMTIMm82XLlvHqq6+yYsWKPJvMH+To6EjNmjW5cOECAFu2bCEyMjJXz/revXvTpk0bg9sg+YmlvFMUhb+PXmfKP6e4nZ6NWgWvta1GoJ9z3k/ISIYTy+HAArhx6v5+9wbQfCg07AvmRfudJ0m2grOxsaFGjRr5Ll+nTh0WLVpEWlqaPgHv3r0btVpNrVq1cHBwwNPTk3379umbI3Nycjh06BBNmzbV1+Pq6mrQZJicnExUVJR+u2HDhly9epVz587leTVrbm7+2GUF69Spw5UrV7hy5Yr+ajYiIoLExETq1q2b79ecH3nFs3v3blq1asWbb76p3xcZGWlQ5r/vg0aj4eTJk3TocL+1wpjm4uzsbLKzs3NdOZqYmBj8IZGcnExISAgWFhasXr36ofPXqlQqvLy8APjtt9/w9vY2+L0+jIuLC6BLEjdu3ODZZ599aNl7cWVmZj623rzUq1ePevXqcfz4cV5++eU8yxRHc/F/HT16FCcnJ4Ok9ttvvzF06FCWLVum/6PnUVJTU4mMjGTAgAEAjBs3jldffdWgTIMGDfjmm2/o0aNHgWIp764n3uHDlSfYevYmAHU87fmid0MaVMnjSjQuAvZ/D8dXQHaabp+pFdTvDQFDoXJTKKb71ZJkRYH079+fSZMmMWjQICZPnszNmzcZNWoUAwYMwN1d17Fg9OjRTJ8+HX9/f2rXrs3XX39NYmKiQT0dO3Zk0aJF9OjRA0dHRyZOnIiJyf0B4e3ataNt27b07t2br7/+mho1anDmzBn9WEhfX19SU1MJCwujUaNGWFtb5xq6ExwcTIMGDejfvz8zZ84kJyeHN998k3bt2hEQEFCk74uvry/79u3j0qVL2Nra4uzsjL+/P0uWLGHDhg34+fnx888/c+DAAfz8/Azeh9DQUP7991+qV6+e53tlTHOxvb097dq147333sPKygofHx+2b9/OkiVL+PrrrwFdgu3SpQvp6en88ssvJCcnk5ycDOj+CLj3e/nyyy/p2rUrarWav/76i+nTp7N8+XL98WvXrtGpUyeWLFlCixa6Tn0LFy6kTp06uLq6Eh4ezujRo3nnnXf097j37dvHgQMHeOqpp3ByciIyMpIJEyZQvXr1Ql3F3rNlyxays7PzHEsNRd9c/M8//xAXF0fLli2xtLRk06ZNfPbZZwZ/HC1dupRBgwbx7bffEhgYSGxsLKBrtXFw0CWGMWPG0KNHD3x8fLh+/TqTJk3CxMSEfv36AeDh4ZFnP4SqVavqP1f5iaW802oVBvy0j8ibaZibqBkd7M+wttUwe3BKREWBqB2wZxZc2Hx/v0tNXWJt9BJYORV/sEb3Z64ACjuEp6zLawhPfo4/bghPdna2Mnr0aMXe3l5xdHRUQkNDcw3hSUpKUvr27avY29sr3t7eyqJFi3IN4YmPj1eGDBmiVKpUSbG0tFTq16+vrFmzRn/8jTfeUCpVqlQkQ3ge9M033yg+Pj767YULFyqP+69y9uxZpWXLloqVlZV+CE9GRoYyePBgxcHBQXF0dFSGDx+ujBs3zuB8WVlZyvDhwxVnZ2fFzc1NmTZtWpEP4YmJiVEGDx6seHl5KZaWlkqtWrWUr776StFqtYqiKHkOieKBoUj3dOjQQXFwcFAsLS2VwMBAZe3atQbniYqKUgBl69at+n1jx45V3N3dFTMzM8Xf39/gvIpy/7Pk7OysWFhYKL6+vsobb7yhXL161aBuQFm4cOFDX+Pjhr6MHj26WIfwrFu3TmncuLFia2ur2NjYKI0aNVLmzZunaDQafZl27drl+R4/+Lvu27ev4unpqZibmyuVK1dW+vbtq1y4cOGR5+Y/w3XyE8t/lcfvsU2nYpXn5uxSzsf9ZwhZTpaiHF+hKPPaKMoke91jsqOiLPs/RYnaqSgPfP4Ko6BDeFSKUoDeKBVUcnIyDg4O+mEND8rIyCAqKgo/P78KuURURTBp0iS2b9/+0PtdonhFRUVRs2ZNIiIi8Pf3L+1wnkhl/XtMURRWHLyKpbkJzzby0u/XahXU99Z/zUyBwz/D3v9B0hXdPlMraDoAWg4H54f3HSiIR+WDvJTJGZ/mzJmDr68vlpaWBAYGsn///oeWnT9/Pm3atMHJyQknJyeCg4NzlX9w+sB7j7ymXxMiL+vWreOLL74o7TAqrLVr1zJs2DBJsBVUTNIdhiw6wPt/HuejlSe4kZKhP6ZWqyA5BjZPhm/qwYbxugRr7QIdPoTQCHj6yyJLsIVR5u7J/v7774SGhjJv3jwCAwOZOXMmISEhnD17Ns/7Utu2baNfv360atUKS0tLPv/8c7p06cKpU6cMOi507dqVhQsX6refpA4Aong96o88UfxGjBhR2iGIUqAoCisOXWXqmghSMnIwN1XzZocaVLK5+9194zTsmQ3Hfwft3R75lWpA0Ejd/Vazkhuh8Chlrrk4MDCQ5s2bM3v2bEDX09Db25tRo0Y9cp7UezQaDU5OTsyePZuBAwcCuivZxMTEXJNo55c0FwshnmRl7XssNimD8X8d1/ccbuTtyFcvNqSGqy1c3g27v4XzG+8/oWoQtHoLanaFYl6AvqDNxUVyJRsWFkZYWBg3btzINb7wvzOtPEpWVhaHDh1i/Pjx+n1qtZrg4GDCw8PzVUd6ejrZ2dm5pirbtm0bbm5uODk50bFjRz755BMqVaqUZx2ZmZkGwwfu9bQUQghRvBLTswiZuYOkO9mYm6oJ7VyTV1v7Ynp5B6z5AqL33C2pgjo9dMnVu3mpxvwoRifZKVOm8PHHHxMQEICnp6dRc2PeunULjUajHwpyj7u7O2fOnMlXHWPHjs01k03Xrl15/vnn8fPzIzIykg8++IBu3boRHh5uMGzknmnTphlMhyeEEKJkOFqb81yTyhy5ksiM3g3wT9kPi96Aq3dv25iYQ5P/0zULV6peusHmg9FJdt68eSxatEg/WLo0TZ8+nWXLlrFt2zaDJo+XXnpJ/3ODBg1o2LAh1atXZ9u2bXTq1ClXPePHjyc0NFS/nZycbDA1nxBCiKKhKAp/Hr5GgI+Tfvm5cV1rYRa5EZPVPeH6YV1BU0toNgRavwX2Xo+osWwxOslmZWXRqlWroogFFxcXTExMci3qHBcX98iJ4EE32fb06dPZvHlznnOVPqhatWq4uLhw4cKFPJOshYWFdIwSQohidiM5g/F/nSDszA1a+Dmz7LWWqC/vxHLzFLh2UFfIzFo3eUSrt8DuMSvplEFG3yF+9dVXWbp0aVHEgrm5Oc2aNSMsLEy/T6vVEhYW9sjZX7744gumTp3K+vXr8zWTz9WrV4mPj8fT07NI4hZCCFEwa45fp8vMHYSduYG5iZoXPW+i+uU5WNxDl2DNrKH1aBh9HEI+LZcJForgSjYjI4MffvhBfwX530nF703dll+hoaEMGjSIgIAAWrRowcyZM0lLS2PIkCEADBw4kMqVKzNtmm4pos8//5yJEyeydOlSfH199VOV2draYmtrS2pqKlOmTKF37954eHgQGRnJ+++/T40aNQgJefgqOeLJoVKpWLlyJb169eLSpUv4+flx5MiRAq+GI4QwXlJ6NhNXn+Tvo7oVc0Lck5hR6R/sDq/VFVCbQcAQaDOm3CbWBxl9JXv8+HEaN26MWq3m5MmTHDlyRP84evRogevr27cvM2bMYOLEiTRu3JijR4+yfv16fWeo6OhogwnV586dS1ZWFi+88AKenp76x4wZMwDdROjHjx/n2WefpWbNmrzyyis0a9aMnTt3Vvgm4Qcn6TAzM8Pd3Z3OnTuzYMGCAq9Cs2jRoofOE1ueZGdnM3bsWBo0aICNjQ1eXl4MHDiQ69evG5RLSEigf//+2Nvb4+joyCuvvPLIhdOFEHDhRiohM3fw99HrVFHHs6bqb8xLHoHdxbWAChq+BKMO6iaQeAISLBTBlezWrVuLIg4DI0eOZOTIkXke++/UdpcuXXpkXVZWVgVeMLwiuTdJh0ajIS4ujvXr1zN69Gj++OMPVq9enWvpuyddeno6hw8fZsKECTRq1Ijbt28zevRonn32WQ4ePKgv179/f2JiYti0aRPZ2dkMGTKEYcOGFdmtEyGeRN7OVvhYpjNG+Z3nNRtQ38jSHaj1DHT8CNyLdnWsMsGomZL/48qVK8qVK1eKssoyoaItEBAWFqYAyvz58/X7vvrqK6V+/fqKtbW1UqVKFWX48OH6RQHymmD+3oT9S5YsUZo1a6bY2toq7u7uSr9+/ZS4uLiHxvTdd98p9erV02+vXLlSAZS5c+fq93Xq1En58MMP9durVq1SmjRpolhYWCh+fn7K5MmTlezsbP1xHphA/d4k9keOHMn3+7R//34FUC5fvqwoiqJEREQogHLgwAF9mXXr1ikqlUq5du1avusVoqwozu+xE1cTlawcjaLcSVKULZ8qmk8870/cv/AZRYneX+TnLE4FXSDA6OZirVbLxx9/jIODAz4+Pvj4+ODo6MjUqVOLZOHrcklRICutdB5FMIFXx44dadSoEX/99Zd+n1qtZtasWZw6dYrFixezZcsW3n//fQBatWrFzJkzsbe3JyYmhpiYGP2yWtnZ2UydOpVjx46xatUqLl26xODBgx967nbt2hEREcHNm7qZXrZv346Li4u+BSM7O5vw8HDat28PwM6dOxk4cCCjR48mIiKC77//nkWLFvHpp58a/T7ck5SUhEql0jeHh4eH4+joaNDJLjg4GLVabbAguxAVWVaOli83nOHF2VsI/3kKfNsItn+OOjsNPBvD//0Fg/4p0xNJFAWj2wI//PBDfvrpJ6ZPn07r1q0B2LVrF5MnTyYjI6NIv+zKjex0+KyUxnF9cB3MbYyupnbt2hw/fly//fbbb+t/9vX15ZNPPuGNN97gf//7H+bm5jg4OKBSqXINtRo6dKj+52rVqjFr1iyaN29Oamoqtra2uc5bv359nJ2d2b59Oy+88ALbtm3j3Xff5dtvvwV08whnZ2frh41NmTKFcePGMWjQIP05pk6dyvvvv8+kSZOMfh8yMjIYO3Ys/fr100+hFhsbm2sebVNTU5ydnfUd74SoyM7GpvDuskPUu7mGLeZ/4nkpQXegkr+uWbhuz2JbJL2sMTrJLl68mB9//JFnn31Wv69hw4ZUrlyZN998s2Im2SeAoigGs3dt3ryZadOmcebMGZKTk8nJySEjI4P09PRci6U/6NChQ0yePJljx45x+/ZtfetGdHQ0devmvv+iUqlo27Yt27ZtIzg4mIiICN58802++OILzpw5w/bt22nevLn+nMeOHWP37t0GnzONRpOv2B4nOzubPn36oCgKc+fOLXQ9QlQUGq3CTzsjObzpN2aql1LD7G6HQfvK0H4cNHoZTCpWPw+jX21CQgK1a9fOtb927dokJCQYW335ZGatu6IsrXMXgdOnT+Pn5wfoOpd1796d4cOH8+mnn+Ls7MyuXbt45ZVXyMrKemgiS0tLIyQkhJCQEH799VdcXV2Jjo4mJCSErKysh567ffv2/PDDD+zcuZMmTZpgb2+vT7zbt2+nXbt2+rL3hmg9//zzueoxZqLzewn28uXLbNmyxWAicA8PD27cuGFQPicnh4SEhMdOmiLEk+rq7XTm/LKcXjfnMsxUNw2u1tIZdbsxEPAKmJX+wgOlwegk26hRI2bPns2sWbMM9s+ePZtGjRoZW335pFIVSZNtadmyZQsnTpzgnXfeAXRXo1qtlq+++gr13RUuli9fbvAcc3NzNBqNwb4zZ84QHx/P9OnT9dNSPthD92HatWvH22+/zYoVK/T3Xtu3b8/mzZvZvXs37777rr5s06ZNOXv2LDVq1Cj06/2vewn2/PnzbN26NddCEkFBQSQmJnLo0CGaNWsG6N4zrVZLYGBgkcUhRLmRcBGHfycyLf4fUEOO2gKToDdRt3kHLB1KO7pSZXSS/eKLL3jmmWfYvHmzflam8PBwrly5wtq1a40OUBSvzMxMYmNjDYbwTJs2je7du+uXCqxRowbZ2dl899139OjRg927dzNv3jyDenx9fUlNTSUsLIxGjRphbW1N1apVMTc357vvvuONN97g5MmTTJ069bExNWzYECcnJ5YuXcqaNWsAXZIdM2YMKpVKf+8fYOLEiXTv3p2qVavywgsvoFarOXbsGCdPnuSTTz4p8PuRnZ3NCy+8wOHDh1mzZg0ajUZ/n9XZ2Rlzc3Pq1KlD165dee2115g3bx7Z2dmMHDmSl156CS+v8jOnqhDGyk65idmuGXDgJ+y02SioSKvTB9uuE8GhSmmHVzYURZfma9euKR988IHy/PPPK88//7zy4YcfPlFDGZ7kITzcHXJjamqquLq6KsHBwcqCBQsUjUZjUPbrr79WPD09FSsrKyUkJERZsmSJAii3b9/Wl3njjTeUSpUqGQzhWbp0qeLr66tYWFgoQUFByurVq/M1hKZnz56KqampfpiQRqNRnJyclJYtW+Yqu379eqVVq1aKlZWVYm9vr7Ro0UL54Ycf9McpwBCee8fzemzdulVfLj4+XunXr59ia2ur2NvbK0OGDNHHKkR5U+DvsexM5dI/05WUSR73h+P8/LyixJwo3kDLgIIO4Slzi7aXRbJouxDiSVaQ77GcsxtIWvkelTIuA3DJrDq+L30F1TuURKilrkQWbT9+/Dj169dHrVYbDPPIy+NWxBFCCFEO3LpA+j9jsb68mUrATcWerZWH8/TAMWBpXtrRlVmFSrKNGzfWjxVs3LgxKpWKvC6IVSpVrs4wQgghypGMZJTtX6DdOxdrJYdsxYSlqqfx7DmRPk1rlnZ0ZV6hkmxUVBSurq76n4UQQjxhFAWOL4eNH6FKu4EJsFXTiDWeoxjTvzueDlalHWG5UKgk6+Pjo//58uXLtGrVKtdE8jk5OezZs8egrBBCiHIgLgLWjoHLuwFQnKsz1/JVTGuH8OVT1VCrK8ZsTUXB6CE8HTp0ICYmJtc0c0lJSXTo0EGai4UQorzISoPtU1H2zkWlaFBMrVC1ew9V0EiGm5gbzAIn8sfoJKv8Z/q9e+Lj47GxKb8TMhSUdNIWQpRXiqJA9h34tQ/EHUQFbNAEsLvaGD5u0w0ASa+FU+gke28aO5VKxeDBgw0WQNdoNBw/flw/ifuTzMzMDNCtQ2plJfcohBDljCab9BuXIfkaZrdOcRVXPsoazCHz5nzeREaHGKvQSdbBQTdVlqIo2NnZGSQYc3NzWrZsyWuvvWZ8hGWciYkJjo6O+rlsra2tpUlFCFH2KQpKeiLpiXHciE/A4fIGfsrsyDc5vald1YO1LzXB27lo5kKvyAqdZBcuXAjoptN77733jFrtpLy7Nyn8fyeNF0KIMkmTDXduQ04GaLLQRO/j9ZO1OEk13mhfndDONTEzMXq5cQEYPeNTVFQUOTk5+Pv7G+w/f/48ZmZm+Pr6GlN9mZDfGT40Gg3Z2dklGJkQQhSAVgPHfoO9c0GTCWpzshoPJGRPHbIx4Zu+jWjj71raUZZpJTLj04MGDx7M0KFDcyXZffv28eOPP7Jt2zZjT1FumJiYYGJiUtphCCFEbvGRsOpNuLJXt12tPXT/BkvnanxfMxFPBytc7SweWYUoOKPbA44cOWKwKso9LVu25OjRo8ZWL4QQwhhare7KdW5ruLIXjZkNX5q/ycr6s8G5GgANqzhKgi0mRl/JqlQqUlJScu1PSkqSMbJCCFGaEi7C3yP1k0rEVgrkpbj/41JOJfy3XeTZxlUwkYklipXRV7Jt27Zl2rRpBglVo9Ewbdo0nnrqKWOrF0IIUVBaLeyfr7t6vbwbxcyGFR6htLz2FpdyKtGpthsr3giSBFsCjL6S/fzzz2nbti21atWiTZs2AOzcuZPk5GS2bNlidIBCCCEKIOkarHoDonYAkO4VxLDkIey6ZIuJWsV7IbUY1kamRiwpRifZunXrcvz4cWbPns2xY8ewsrJi4MCBjBw5Emdn56KIUQghRH5ErIbVoyAjEcysSWnzES03+5GWpeBub8Hsl5vS3Fe+l0tSkQyE8vLy4rPPPuPff//ljz/+YOLEiUYl2Dlz5uDr64ulpSWBgYHs37//oWXnz59PmzZtcHJywsnJieDg4FzlFUVh4sSJeHp6YmVlRXBwMOfPny90fEIIUaZkpcHqt2D5AF2C9WoCb+zCru0IBgRVo42/C/++1UYSbCkw+kp2x44djzzetm3bAtX3+++/Exoayrx58wgMDGTmzJmEhIRw9uzZXIsQAGzbto1+/frRqlUrLC0t+fzzz+nSpQunTp2icuXKAHzxxRfMmjWLxYsX4+fnx4QJEwgJCSEiIgJLS8sCxSeEEGXK9aPw56sQfx5QcbvpCNJavUeVSo4AjOlSE5VKJfdfS4nRk1Go1bkvhh+cVrCgPYwDAwNp3rw5s2fPBkCr1eLt7c2oUaMYN27cY5+v0WhwcnJi9uzZDBw4EEVR8PLy4t1332XMmDGAruezu7s7ixYt4qWXXnpsnQUdfCyEEMVOUWDv/2DTJNBmg50X+5tMY+h2K6q72bLi9SDMTWXWpqJW0Hxg9G/g9u3bBo8bN26wfv16mjdvzsaNGwtUV1ZWFocOHSI4OPh+gGo1wcHBhIeH56uO9PR0srOz9c3VUVFRxMbGGtTp4OBAYGDgQ+vMzMwkOTnZ4CGEEGXGnUTd1euGD0CbjVK7O9/VWkSfjWakZuZgYaomLTOntKMUFEFz8b2FAh7UuXNnzM3NCQ0N5dChQ/mu69atW2g0Gtzd3Q32u7u7c+bMmXzVMXbsWLy8vPRJNTY2Vl/Hf+u8d+y/pk2bxpQpU/IdtxBClJjrR2D5QEiMBpUJ6R2n8sa5AHYcvQXAK0/5Mb5bbUxl7uEyodh+C+7u7pw9e7a4qs/T9OnTWbZsGStXrjTqXuv48eNJSkrSP65cuVKEUQohRCEoChxcCD910SVYJ1+ieqwgJLwOO87fwtJMzbcvNWZC97qSYMsQo69kjx8/brCtKAoxMTFMnz6dxo0bF6guFxcXTExMiIuLM9gfFxenX+nmYWbMmMH06dPZvHkzDRveXwPx3vPi4uLw9PQ0qPNh8VlYWBisjyuEEKUqKx3+DdVN7g9Q62mUXv/j3YURXElIxNvZiu//L4C6XtJnpKwxOsk2btwYlUrFf/tPtWzZkgULFhSoLnNzc5o1a0ZYWBi9evUCdB2fwsLCGDly5EOf98UXX/Dpp5+yYcMGAgICDI75+fnh4eFBWFiYPqkmJyezb98+hg8fXqD4hBCixCVchGX/BzdOgUoNnSZB69GoVCpm9m3CFxvO8Emv+jham5d2pCIPRifZqKgog221Wo2rq2uhm2tDQ0MZNGgQAQEBtGjRgpkzZ5KWlsaQIUMAGDhwIJUrV2batGmAbsapiRMnsnTpUnx9ffX3WW1tbbG1tUWlUvH222/zySef4O/vrx/C4+XlpU/kQghRJl3aDb//H9xJABs3krr/wLbMmvS8O4KjaiVrZr/ctJSDFI9idJL18fEpijj0+vbty82bN5k4cSKxsbE0btyY9evX6zsuRUdHGwwbmjt3LllZWbzwwgsG9UyaNInJkycD8P7775OWlsawYcNITEzkqaeeYv369TJGVghRdh35Bf55Wzc8x6spFzp9z6AVV7medBR7KzM61Mo9b4Aoewo1TnbWrFn5LvvWW28VtPoyR8bJCiFKjFYDmyfBnu902/WeI6zWFN764zRpWRp8Klnz48AA/N3tSjfOCqqg+aBQSdbPz89g++bNm6Snp+Po6AhAYmIi1tbWuLm5cfHixYJWX+ZIkhVClIjMVPjrNTi7FgCl3Vjm0YcvNp5FUaBV9UrMebkpTjZy/7W0lMhkFFFRUfrHp59+SuPGjTl9+jQJCQkkJCRw+vRpmjZtytSpUwtTvRBCVDyJV2BBiC7BmliQ1Ws+7954ms836BJs/8CqLB7aQhJsOWP0tIrVq1fnjz/+oEmTJgb7Dx06xAsvvJCrY1R5JFeyQohideUALHsZ0m6AjRv0+41/EyozYulhTNQqJveoy4Ag39KOUlDwfGB0x6eYmBhycnJP36XRaHKNdxVCCPEfJ/6AVW+CJhPcG0C/38DRm2eqwIlr1XmqhgtP+buUdpSikIyeFqRTp068/vrrHD58WL/v0KFDDB8+3GC+YCGEEA/QamHrZ/DnK7oEW+tptrVeQpL5/Yl3xnWrLQm2nDM6yS5YsAAPDw8CAgL0MyW1aNECd3d3fvzxx6KIUQghnixZ6fDnUNj+OQBKq7f4znUSg5eeZsTSw+RotKUcoCgqRjcXu7q6snbtWs6dO8fp06dRqVTUrl2bmjVrFkV8QgjxZEmJhd/6wfXDoDYj++mvCT3fgH+ORQJQw822lAMURcnoJHtPzZo18ff3BwzXkxVCCHFXzDFY+hKkXAcrZxJ6LGDwFjOOX72OqVrFxz3r83Jg1dKOUhShIlmqYcmSJTRo0AArKyusrKxo2LAhP//8c1FULYQQT4bT/8CCrroE61KTMz1W0W2VhuNXk3CyNuOXVwMlwT6BjL6S/frrr5kwYQIjR46kdevWAOzatYs33niDW7du8c477xgdpBBClFuKAru+gbC7a1RX74im9wJGzjtOXHIm/m62/DSoOVUrWZdunKJYGD1O1s/PjylTpjBw4ECD/YsXL2by5MkyTlYIUXHlZMI/o+8vUdf8Neg6HUxMORObzHdbLjD9+QbYWZqVbpwi30plnGyrVq1y7W/VqhUxMTHGVi+EEOVT2i3dCjrR4aAyIavzNA669aaVie5rt7aHPXNkBZ0nntH3ZGvUqMHy5ctz7f/999/1HaGEEKJCuXUefgzWJVgLB+Kf+5XnD9Vl0ML9HLiUUNrRiRJk9JXslClT6Nu3Lzt27NDfk929ezdhYWF5Jl8hhHiiXdoFy/pDRiI4+hDR8ScG/ZPEzZRknG3MkbEXFYvRSbZ3797s37+fr7/+mlWrVgFQp04d9u/fn2s+YyGEeKIdWwZ/j9StAVulOesafM3o5dfIytFSy92OHwcF4O0sHZwqEqOSbHZ2Nq+//joTJkzgl19+KaqYhBCifFEU2DYdtk/Xbdbpybd27zJz5RUAguu4MfOlJthaFNnUBKKcMOqerJmZGX/++WdRxSKEEOVPTiasfF2fYGn9NmtqfcbMHVcBeKNddb4fECAJtoIy+rfeq1cvVq1aJeNhhRAVT3qCrgfx5d2gMoHuX0OzwXRXFLaevUXrGi70blaltKMUpcjoJOvv78/HH3/M7t27adasGTY2NgbH33rrLWNPIYQQZU98JCztA/EXwNyOqI5z8WjwNFboppb9um/j0o5QlAFFMhnFQytXqbh48aIx1ZcJMhmFEMJA9F7dIuvp8WBfhc1NZ/Pmpgw613Xnu35NUKulD/GTqsQno3gSZnQSQoh8O/knrBwOmkwUz8bM9viUr9YnAZCt0ZKl0WKpNinlIEVZYVSS3bt3L//88w9ZWVl06tSJrl27FlVcQghRtigK7Poawj4GIMe/G29ljWBtuC7BjuhQnXc715KrWGGg0En2jz/+oG/fvlhZWWFmZsbXX3/N559/zpgxY4oyPiGEKH2abPj3XTi8GICUxq/RN+oZIuKSMTdV80XvhvRqUrmUgxRlUaGH8EybNo3XXnuNpKQkbt++zSeffMJnn31WlLEJIUTpy0iGpX11CValRtv1c3pH9SAiLh0XWwt+H9ZSEqx4qEJ3fLK1teXo0aPUqFEDgKysLGxsbLh27Rpubm5FGmRpk45PQlRQSdd0PYjjToKZNfT+CWo/zcFLCUxbd4bZLzfB08GqtKMUJaig+aDQV7Lp6ekGJzA3N8fS0pLU1NTCVqk3Z84cfH19sbS0JDAwkP379z+07KlTp+jduze+vr6oVCpmzpyZq8zkyZNRqVQGj9q1axsdpxDiCRZ7QjfJf9xJFBs3Irv/DrWfBiDA15k/3giSBCsey6iOTz/++CO2trb67ZycHBYtWoSLi4t+X0HHyf7++++EhoYyb948AgMDmTlzJiEhIZw9ezbPK+T09HSqVavGiy+++MgJMerVq8fmzZv126amMvuKEOIhLmyG5YMgKxVNpZq8Y/ohm/9M5y/3ZGp76C4uVCrp4CQer9DNxfeuHB9ZeSHGyQYGBtK8eXNmz54NgFarxdvbm1GjRjFu3LjHxvT222/z9ttvG+yfPHkyq1at4ujRowWK5R5pLhaiAjm0CNaEgqIh3SuIFxNHcCpBjbW5CXP6N6VDrSfrdpgomBIbJ3vp0qXCPvWhsrKyOHToEOPHj9fvU6vVBAcHEx4eblTd58+fx8vLC0tLS4KCgpg2bRpVq1bNs2xmZiaZmZn67eTkZKPOLYQoB7Ra2DJVN0wHiPHpSbdLfUnMhCpOVvw4KEB/FStEfhm9aHtRunXrFhqNBnd3d4P97u7uxMbGFrrewMBAFi1axPr165k7dy5RUVG0adOGlJSUPMtPmzYNBwcH/cPb27vQ5xZClAM5mfDXq/oEe8BnGK3O9SExEwL9nFk98ilJsKJQKsSNyW7duul/btiwIYGBgfj4+LB8+XJeeeWVXOXHjx9PaGiofjs5OVkSrRBPqvQE3SLr0XtAbcrhxh/z4h5fAF4OrMrkHvUwNy1T1yOiHClTSdbFxQUTExPi4uIM9sfFxeHh4VFk53F0dKRmzZpcuHAhz+MWFhZYWFgU2fmEEGVUwkX49UXdJP8W9tD3Zxr6tKXtrYN0ruPGgCDf0o5QlHNl6s8zc3NzmjVrRlhYmH6fVqslLCyMoKCgIjtPamoqkZGReHp6FlmdQohy5soB3RCd+Atk2niRPWgdVGuPqYmaRYObS4IVRaJMJVmA0NBQ5s+fz+LFizl9+jTDhw8nLS2NIUOGADBw4ECDjlFZWVkcPXqUo0ePkpWVxbVr1zh69KjBVeqYMWPYvn07ly5dYs+ePTz33HOYmJjQr1+/En99QogyIGI1LO4O6fHcsK1Nm4QP+fiB4fgy/7AoKkXSXBwZGcnChQuJjIzk22+/xc3NjXXr1lG1alXq1atXoLr69u3LzZs3mThxIrGxsTRu3Jj169frO0NFR0ejVt//2+D69es0adJEvz1jxgxmzJhBu3bt2LZtGwBXr16lX79+xMfH4+rqylNPPcXevXtxdXU1/sULIcoPRYG9/4MNHwIKx60DeenW66RjSY5Wi1arSIIVRcro9WS3b99Ot27daN26NTt27OD06dNUq1aN6dOnc/DgQf7444+iirXUyDhZIZ4AWg2sHwf7fwDgb7NuhKa8jImJGVN61qNfi7yH9AnxoBKbVvGecePG8cknn7Bp0ybMzc31+zt27MjevXuNrV4IIYyXlabrQbz/BxRUfK4dwOiU/8PV3obfX28pCVYUG6Obi0+cOMHSpUtz7Xdzc+PWrVvGVi+EEMZJidNN8h9zFMXUkndz3uSvrABaVnPmu35NcbWTkQSi+BidZB0dHYmJicHPz89g/5EjR6hcWZZ/EkKUoutHdFewydfAuhKqfssISa5KlevJjO7kj4ncfxXFzOjm4pdeeomxY8cSGxuLSqVCq9Wye/duxowZw8CBA4siRiGEKLjjK2BBV0i+Rrp9NXhlE3i3IKSeB6Gda0qCFSXC6CT72WefUbt2bby9vUlNTaVu3bq0bduWVq1a8dFHHxVFjEIIkX9aDWyaqJsmMSeDME0TQlImcsNMWtZEyTO6d/E90dHRnDx5ktTUVJo0aYK/v39RVFsmSO9iIcqJO4nw56twYRMAc3Ke5WtNH0Z0rCXNw6JIlNgqPPfs2rWLp556iqpVqz50VRshhCh2N8+hLOuHKv4CdxRz3st+nYO2HVjSpxGta7g8/vlCFAOjm4s7duyIn58fH3zwAREREUURkxBCFMy5jSg/dkQVf4FrSiVeyJqEtt5zrH+7jSRYUaqMTrLXr1/n3XffZfv27dSvX5/GjRvz5ZdfcvXq1aKITwghHk5RYNc3sLQPqswULlo34iVlGkNe6MWcl5viaG3++DqEKEZFdk8WICoqiqVLl/Lbb79x5swZ2rZty5YtW4qq+lIj92SFKIOy0slaOQLz03/ptpsNIa3TZ8TfgaqVrEs3NvHEKmg+KNIkC6DRaFi3bh0TJkzg+PHjaDSaoqy+VEiSFaKMSbxC6pKXsE04SQ4mqJ/+EnWL3GtDC1HUSnxaxXt2797Nm2++iaenJy+//DL169fn33//LarqhRACgOyLu0mb3QbbhJPcUuwJtfyYm7X7l3ZYQuTJ6N7F48ePZ9myZVy/fp3OnTvz7bff0rNnT6ytpblGCFG0rofNxXXnR9iQwymtD//UmcG05ztiY1EkC4oJUeSM/mTu2LGD9957jz59+uDiIr34hBBFLycrk4iFI2gYswKAjapWqHvPYVyjaqUcmRCPZnSS3b17d1HEIYQQeUu9iXr5IBrG7EarqFjlPJS2Qz/Dxc6ytCMT4rEKlWRXr15Nt27dMDMzY/Xq1Y8s++yzzxYqMCFExabVKijR4Zj8ORR1SgwaM1v2Nfmc57r1R6WSmZtE+VCo3sVqtZrY2Fjc3NxQqx/ed0qlUknvYiFEgV2JT2Prokn8X+oC1IoGXGtDnyXgWqu0QxMVXIlMq6jVavP8WQghjKEoCn/sjsBx02gGqg4AkFPvBUx7zgJzm1KOToiCM3oIz5IlS8jMzMy1PysriyVLlhhbvRCigohNymDC97/RfONzdFYdIBszEtpPx/SFHyXBinLL6MkoTExMiImJwc3NzWB/fHw8bm5u0lwshHgkRVFYdfQaR/7+jg+Un7BUZZNi6YlN/19Rezcr7fCEMFDiq/AoipJnJ4SrV6/i4OBgbPVCiCfctRvxaFa+ycfq7aCCNJ9O2PX9EaydSzs0IYxW6CTbpEkTVCoVKpWKTp06YWp6vyqNRkNUVBRdu3YtkiCFEE+o+Eiq/DWQF9Qn0aJG6fAhNm1C4REdKoUoTwqdZHv16gXA0aNHCQkJwdbWVn/M3NwcX19fevfubXSAQogny5WEdCb8fZJJNSLx2/U+ZCaDjSvq3j9BtXalHZ4QRarQSXbSpEkA+Pr60rdvXywtZWC4EOLhtFqFn/de5qv1Jxml/RW/y2t1B6oGwQsLwd6zdAMUohgYfU920KBBRRGHEOIJdj4uhQ9WniD6UiQ/mc+iuek53YFWb0GniWBiVroBClFMjL7xodFomDFjBi1atMDDwwNnZ2eDR2HMmTMHX19fLC0tCQwMZP/+/Q8te+rUKXr37o2vry8qlYqZM2caXacQomikZ+Uwfd0Zun27E7Ponay1+IDm6nMoFvbQ91foMlUSrHiiGZ1kp0yZwtdff03fvn1JSkoiNDSU559/HrVazeTJkwtc3++//05oaCiTJk3i8OHDNGrUiJCQEG7cuJFn+fT0dKpVq8b06dPx8PAokjqFEEVj46k4vt9+njdUf/GL+TQqqZLBvQGqYdugTvfSDk+I4qcYqVq1asqaNWsURVEUW1tb5cKFC4qiKMq3336r9OvXr8D1tWjRQhkxYoR+W6PRKF5eXsq0adMe+1wfHx/lm2++KdI6FUVRkpKSFEBJSkrKV3khKrKsHI3+Z01SrHL6846KMsle91g1QlGy0ksxOiGMU9B8YPSVbGxsLA0aNADA1taWpKQkALp3717gRduzsrI4dOgQwcHB+n1qtZrg4GDCw8MLFV9h6szMzCQ5OdngIYR4tMwcDXO2XqDLNztIzcyBi9tQf/8UtdMPgqkV9JwDPWeDmVVphypEiTE6yVapUoWYmBgAqlevzsaNGwE4cOAAFhYWBarr1q1baDQa3N3dDfa7u7sTGxtbqPgKU+e0adNwcHDQP7y9vQt1biEqAkVR2HImjm4zd/LlhrNE30om8vfxsKQXpN0At7owbBs0+b/SDlWIEmd0kn3uuecICwsDYNSoUUyYMAF/f38GDhzI0KFDjQ6wNIwfP56kpCT948qVK6UdkhBl0rm4FAYu2M/QRQe5eCuNujYp7Ks8k0YXfwAUaDoIXg0Dt9qlHaoQpcLoITzTp0/X/9y3b1+qVq1KeHg4/v7+9OjRo0B1ubi4YGJiQlxcnMH+uLi4h3ZqKo46LSwsCnwVLkRFotUqTPnnFL/si0ajVTAzgRm1ztPj+kzU8bfB3BZ6fAsNXijtUIUoVUU+d1lQUBChoaEFTrCgmymqWbNm+itj0C2lFxYWRlBQUKHiKY46hajo1GoViXey0WgV+tRUc8x/IT0vTkKdcRs8G8HrOyTBCkEhr2RXr16d77LPPvtsgeoODQ1l0KBBBAQE0KJFC2bOnElaWhpDhgwBYODAgVSuXJlp06YBuo5NERER+p+vXbvG0aNHsbW1pUaNGvmqUwjxeNvO3qCGmy1VnKwBGNe1FqOc9lHj8GeQmQRqM2g3Fp56W8a+CnFXoZLsvXmLH0elUhV4qbu+ffty8+ZNJk6cSGxsLI0bN2b9+vX6jkvR0dGoH5g8/Pr16zRp0kS/PWPGDGbMmEG7du3Ytm1bvuoUQjzc+bgUPl17mm1nb/JMQ0/mvNwUEq/g+c9bELlFV8irqa73sHvd0g1WiDLG6PVkKwJZT1ZURLFJGczcfI7lB6+gVcDMRMWQVj6Mc92LevNEyEoFEwvo+CG0HAEmRnfxEKLMK/H1ZIUQT5bkjGzmbYtkwe4oMrK1AHSp686E1tZ473wfDuzUFfRuqRv36uJfitEKUbYZnWQ//vjjRx6fOHGisacQQpSgX/Ze5n/bIgEI8HFifLeaNItdAcs+hux0MLOGTpOgxWugNinlaIUo24xOsitXrjTYzs7OJioqClNTU6pXry5JVogyTqtViE/LwtVON2xtcCtfdp67xZDWvnR2S0a1uj9c2asr7NsGnp0FztVKMWIhyg+jk+yRI0dy7UtOTmbw4ME899xzxlYvhChGu87fYvr60ygK/DPyKdRqFdbmpvz2SgDsnQMrP4OcDN24184fQ7MhoC7ykX9CPLGK5Z6svb09U6ZMoUePHgwYMKA4TiGEMMKp60lMX3eGnedvAWBrYcqFm6nUdLeDG6dh1Ztw/bCucPWO0GMWOMr0okIUVLF1fLo3JaEQouy4dCuNbzafY/Wx6yh3ewz3D/RhVMcaVLJSw/YvYfvnoM0GCwfo+hk07g8qVWmHLkS5ZHSSnTVrlsG2oijExMTw888/061bN2OrF0IUkZPXkug5ZzcarW7UXo9GXozpUhOfSjYQcwx+GQGxJ3SFa3aF7t+AvVcpRixE+Wd0kv3mm28MttVqNa6urgwaNIjx48cbW70QwgjZGi1mJrp7qPW87KnraY+LrTnvdqlF/coOkJMJWz6BXd+ANgesnKDbF9DgRbl6FaIIGJ1ko6KiiiIOIUQRSkrP5vsdkaw+dp31b7fF1sIUlUrFsmEtsbG4+9/+4jZY+x7cOqfbrvMsPPMV2LqVWtxCPGlkMgohniCpmTks2h3F9zsukpKRA8DfR6/RP9AHQJdgk67Bxg/h1N3hdzau8PQMqNerlKIW4slldJLNyMjgu+++Y+vWrdy4cQOtVmtw/PDhw8aeQgjxGEl3slm0+xILdkeRdCcbgFrudrzbpSad696dozsnC/bNhW2fQ3YaqNTQ/DXo8AFYOZZe8EI8wYxOsq+88gobN27khRdeoEWLFqjkPo4QJSo5I5u2X2zVJ9dqrjaM7uRPj4ZeqNV3/z9e3H63afisbts7UHf16tmwlKIWomIwOsmuWbOGtWvX0rp166KIRwiRD+lZOVib6/772lua0a6mK2dikxnZ0Z9nGnhici+5Jl+HjR/ByT9129YuukklGvWTSSWEKAFGJ9nKlStjZ2dXFLEIIR7jRnIG83deZNn+K6we9RR+LjYAfPJcfWzNTe9fuWqyYd882DZdt1qOSg0Br+hWzLFyKsVXIETFYnSS/eqrrxg7dizz5s3Dx8enKGISQvxH5M1U5u+4yF+Hr5Gl0fV7WHnkGqGdawK6q1m9qJ2wdgzcPKPbrtJc1zTs1biEoxZCGJ1kAwICyMjIoFq1alhbW2NmZmZwPCEhwdhTCFFhHY6+zffbI9kYEce9lZ+bVnVkVCd/2td0NSycHAObJsCJFbpt60oQPEU3Y5M0DQtRKoxOsv369ePatWt89tlnuLu7S8cnIYpIRraGoYsOkJiu69AUXMeNN9pVJ8DX2bCgJhv2/wBbp0FWCqCCgKHQ8SOwds5dsRCixBidZPfs2UN4eDiNGjUqiniEqLCycrSEnY6ja30PVCoVlmYmvNLaj+iEdIa1rYa/ex59H6J2wLqxcCNCt105AJ6ZAV5NSjZ4IUSejE6ytWvX5s6dO0URixAVUmpmDsv2R/PTrihikjJYOKQ5HWrpZl0a1ck/7yclXISNE+DMGt22lTN0ngKN/0+ahoUoQ4xOstOnT+fdd9/l008/pUGDBrnuydrb2xt7CiGeSNHx6SwOv8TyA1dIydTNzuRqZ0Hq3Zma8pSRDDu+1PUc1mSBykTXNNzhA2kaFqIMUinKve4UhaO++1fzf+/FKoqCSqVCo9EYU32ZkJycjIODA0lJSfJHgzBaSkY27y4/xqbT9zszVXO1YVibavRqUhlLM5PcT9Jq4MjPusn8027q9lXvCCGfgVudkgteiAquoPnA6CvZrVu3GluFEE+8e390gm6B9OiEdBQF2tZ0ZWhrX9r6u94f4/pfUTth/XiIu7sMXaUauuTq30VWyhGijDM6ybZr164o4hDiiXQzJZNf911m9bHr/DPyKWzurobzSa/6OFiZ5d2Z6Z6EKN2QnNP/6LYtHaDdOGj+Kpial8wLEEIYxegku2PHjkceb9u2rbGnEKLcOXYlkSXhl/nn2HX95BF/H73Oy4FVAXIPw3nQnUTY9TXsnXv3vqtad9+1/QdgU6kEohdCFBWjk2z79u1z7Xvw/uyTcE9WiPzIyNaw+th1ftl7meNXk/T7m1R1ZGhrP7rW93h0Bdl3dONdd34NGYm6fdXaQ8g0cK9bbHELIYqP0X39b9++bfC4ceMG69evp3nz5mzcuLFQdc6ZMwdfX18sLS0JDAxk//79jyy/YsUKateujaWlJQ0aNGDt2rUGxwcPHoxKpTJ4dO3atVCxCfEwt1IzGfvncY5fTcLcRM3zTSqz8s1WrHyzNT0aeWFm8pD/bpocOPwzfNcMNk3UJVjX2tBvGQxYJQlWiHLM6CtZBweHXPs6d+6Mubk5oaGhHDp0qED1/f7774SGhjJv3jwCAwOZOXMmISEhnD17Fjc3t1zl9+zZQ79+/Zg2bRrdu3dn6dKl9OrVi8OHD1O/fn19ua5du7Jw4UL9toWFRYHiEuJBORotYWducDommbeDdfMHV3Gy5v8CfajsZMWLzapQyfYxnzFFgTP/QtjH95egs68CHcbfXSUnj17GQohyxeghPA9z5swZAgICSE1NLdDzAgMDad68ObNnzwZAq9Xi7e3NqFGjGDduXK7yffv2JS0tjTVr1uj3tWzZksaNGzNv3jxAdyWbmJjIqlWrCvVaZAiPuOdGSga/77/C0v3RxCRloFLBjvc64O1sXbCKLu2GzZPh6t1WGisnaPOubhF1M8sij1sIUTRKfAjP8ePHDbYVRSEmJobp06fTuHHjAtWVlZXFoUOHGD9+vH6fWq0mODiY8PDwPJ8THh5OaGiowb6QkJBcCXXbtm24ubnh5OREx44d+eSTT6hUSTqRiMfTaBV2nr/J8oNX2Hgqjhyt7u9SZxtz+jb3zntc68PEnYLNU+D8Bt22qRUEvQmt3gIrx6IPXghRqoxOso0bN0alUvHfC+KWLVuyYMGCAtV169YtNBoN7u7uBvvd3d05c+ZMns+JjY3Ns3xsbKx+u2vXrjz//PP4+fkRGRnJBx98QLdu3QgPD8fEJPcXZGZmJpmZmfrt5OTkAr0O8WRZdeQa7644pt9u5uPEgJY+dGvggYVpPhPs7Uu6tV2PLQMU3UxNTQdCu7Fg71kscQshSp/RSTYqKspgW61W4+rqiqVl2Wnyeumll/Q/N2jQgIYNG1K9enW2bdtGp06dcpWfNm0aU6ZMKckQRRmRmaNh46k4LM1M6FxX98dbSH0Pvt50js513ekT4E1drwLcMrh9CXZ+BUeXgvbudIl1e0HHCeBSo8jjF0KULUYn2aJcqN3FxQUTExPi4uIM9sfFxeHhkffwBw8PjwKVB6hWrRouLi5cuHAhzyQ7fvx4gybo5ORkvL29C/JSRDlzJjaZ3w9cYeWRaySmZ1PH057gOm6oVCpsLUzZ+X6Hh8/IlJe8kmv1jrrkWrlpsbwGIUTZU+ghPFu2bKFu3bp5NqUmJSVRr149du7cWaA6zc3NadasGWFhYfp9Wq2WsLAwgoKC8nxOUFCQQXmATZs2PbQ8wNWrV4mPj8fTM+9mOgsLC+zt7Q0e4slzIzmD+Tsu8vS3O+k6cycLd18iMT0bTwdLguu46e+9AvlPsAlRsHqUbjjO4SW6BFu9IwzdCANWSoIVooIp9JXszJkzee211/JMQA4ODrz++ut8/fXXtGnTpkD1hoaGMmjQIAICAmjRogUzZ84kLS2NIUOGADBw4EAqV67MtGnTABg9ejTt2rXjq6++4plnnmHZsmUcPHiQH374AYDU1FSmTJlC79698fDwIDIykvfff58aNWoQEhJS2JcvngBT/ong3xMxAJiZqAiu406f5t609XfFpCBXrQDXDsHuWXB6NSi6GZ6o3lE3DWLVwCKOXAhRXhQ6yR47dozPP//8oce7dOnCjBkzClxv3759uXnzJhMnTiQ2NpbGjRuzfv16feem6Oho/co/AK1atWLp0qV89NFHfPDBB/j7+7Nq1Sr9GFkTExOOHz/O4sWLSUxMxMvLiy5dujB16lQZK1tBaLQKeyJvsfLwNUZ0rEF1V1sAnmtSmdjkDJ5rUplnGnjiZFPA+YAVBS5sht3fwqUHWm2qd4J270PVlkX4KoQQ5VGhx8laWlpy8uRJatTIu/PGhQsXaNCgwROxoLuMky1/FEXh1PVk/j56jb+PXudGiq63+MgONRgTUsu4ynOy4OSfsGcW3IjQ7VObQv0XoNUo8Kj/6OcLIcqtEhsnW7ly5Ucm2ePHjz/0nqcQxSUlI5vvt19kzfHrXIpP1+93tDaje0NPQuo9Zv7gR8lIhsOLIfx/kHJdt8/cFpoNhpbDwaGKccELIZ44hU6yTz/9NBMmTKBr1665huvcuXOHSZMm0b17d6MDFOJxEtOzcLTWNfVamJqwJPwSyRk5WJiq6VjbjV5NKtOhlhvmpoXs55ccA/vmwsGFkHm3o5+tOwS+oVsdRyaREEI8RKGbi+Pi4mjatCkmJiaMHDmSWrV0TXBnzpxhzpw5aDQaDh8+nGuiiPJImovLnos3U1lzPIZ/j8eQnp3Djvc66Fd/WhJ+CQcrMzrVccfWwohRajdOw57ZcPx30Gbr9lXyh9ZvQcO+YCr39IWoaEqsudjd3Z09e/YwfPhwxo8fr5/xSaVSERISwpw5c56IBCvKBq1W4ejVRDZHxLH5dBzn4u7PiW1mouJSfDp+LjYADAzyNeJEGji/UbeWa9T2+/u9W0Lr0VCzK6iNXrxKCFFBGDUZhY+PD2vXruX27dtcuHABRVHw9/fHycmpqOITAoDP15/h+x0X9dumahVP+bvwTANPutT1wMHazLgTZCTDkV9g//e6iSRAt1h67WcgaJQMwxFCFIrRMz4BODk50bx586KoSlRwcckZhJ2+QdjpOF5vV50Wfs4APOXvwtL90bSv5UZwHTfa13QzPrEC3DwHB36Eo79C1t2rY0sH3bzCzV8Dp6Kb0UwIUfEUSZIVorAUReFMbIq+GfjY1ST9MV8XG32SDapWiUMfdS5856UHZSTDqb/gyK/3l5oDcKkFga9Do5fA3Mb48wghKjxJsqLUxCZl0HvuHq4lGo6lbuztSOe67gbDbUxNjEyuWi1c3qVLrBF/Q87dc6pMoGYItHgNqnUAVQFnehJCiEeQJCtKxOX4NLadvUm2RsurbaoB4G5vQY5Wi4Wpmjb+LgTXcadjHTfc7IpwBafEaDj6m645OPHy/f0utaBJf2j4EthJBz0hRPGQJCuKRUa2hr0X49l29ibbz90k6lYaAJVszBna2g+1WoVKpWLRkBb4VrLByrwAC58/TmYKnF0PR3+Bi9uBu6PUzO2gQW9o/H9QJUCuWoUQxU6SrChyH/8TwdL9l8nI1ur3mapVBPg60b6WG1kaLZZqXVKt41lE447vJMLZdboJ+i+EgSbz/jHfNtBkANTpAebWRXM+IYTIB0myolAUReFKwh32RN5i78V4PnmugX7iBwszNRnZWjzsLWlfy5X2tVxpXcMFO8si6A38oLRbcOZfXWK9uP3+hBEATn7QsA80fhmcfIv2vEIIkU+SZEW+xSTdITwynj2R8YRHxht0WOrZuDIdarsB0D+wKs828qK2h51+FqYikxILp//RdV66vPv+snIArrWhzrNQ91lwry/NwUKIUidJVjyUoij6JPnnoau8u+KYwXFTtYomVR0Jqu6Cr8v9IS9VnIq4STY+Es6th4jVcGUf+nusAB4NdUm1Tk9wrVm05xVCCCNJkhV6SenZ7I3SXaWGR8YzqJUvLwdWBaCRtwNqFTSo7EBQdReCqleiua8T1ubF8BHS5OiS6bl1cG4D3DpneLxyANTtqbvH6uxX9OcXQogiIkm2AsvM0eibfvdE3uLU9WQeXC5id+QtfZKt7mrLkYldcLAq4vuq92Qk6RZAP7seLmyCO7fvH1Obgk9rqPW0LrE6VC6eGIQQoohJkq1AMrI13EzJxNtZ15ybmaPllUUH0D6QWKu72hBUvRKtqrsQeHe2JdAt/FCkCVarhdjjcHEbRIbB5T2gzbl/3MoJ/LvoJuSv0Uk31aEQQpQzkmSfYFk5Wo5dTWTPhXjCL97i8OVE6le25683WwNgb2lGSD0P7C3NCKpeiaDqlXC3L8KJIB6kKHA7SpdUL26HqB1wJ8GwjEst3exLtbpBlRZgIh9PIUT5Jt9iT6Cf915mU0QcB6ISuJOtMTgWl5xJtkaL2d1pCuf+X7PiCyTtlm65uIvbdI/EaMPj5nbg+xRUaw/+naFS9eKLRQghSoEk2XJMURTO30jl0OXbvNTcW98TeM+FW+w4dxPQzbDUsnolgqpVolX1Svi52BT9sJp77iTC1QO6hBq1HWJPGB5Xm4F3C11SrdYevJqASTHd4xVCiDJAkmw5c/V2OnsuxLM78hZ7IuO5maKb2ahltUr6Rcv7NPemhZ8zQdUrUdPNDrW6GJKqVgs3z+iS6tX9cOUA3Dqbu5x7A6jWTpdUqwaBhW3RxyKEEGWUJNly4p9j1/lyw1miE9IN9luaqWnu60x61v1OQx1quUGtIg4gPQGuHdIl1Sv7dT9nJucu5+R3vwnYrx3YuhZxIEIIUX5Iki1jUjNzOBCVwO4Lt+jRyItG3o4AWJqZEJ2QjolaRaMqDrSu4UKr6i409XHEwrQIJ9dXFN2sSrHHIeaY7hF7PPf9VAAzG6jcFKo01zUDV2kONi5FF4sQQpRzkmRLmaIoXLyVxtYzN9hy5gYHLiWQrdGNqbE0M9En2aDqlfhpUAAt/JyLZg7g7AxIuAjx5yH+Aty6cP/nB8eoPsi5+v1kWqU5uNWVHsBCCPEI8g1Ziq4n3qHf/L1cjjdsAvZ2tqJ1dRdaVquk32drYUqnOgVc91SrhZTrcOtu8oy/cPfn85B4BYPpCR+kUoNLTfBspJu20LMheDTQjV0VQgiRb5Jki5lGq3ApPo2I68mcjknG0syEtzr5A+Bub0lKRg7mJmoCqznToZYbHWq76TswPVZ2BiRf0zXlJl194HHl/s8PLvn2Xxb2UKkGuPhDJX/dEBoXf90VqywJJ4QQRiuTSXbOnDl8+eWXxMbG0qhRI7777jtatGjx0PIrVqxgwoQJXLp0CX9/fz7//HOefvpp/XFFUZg0aRLz588nMTGR1q1bM3fuXPz9/Ysl/uUHr3D0SiIR15M5G5tiMFa1sqOVPsmaqFUsHtKCaq422Fg88KvIyYLUON290dRY3b8pMbn/fViz7oPUprrOSJVqgEuNu8n0bmK1cZWVaoQQohiVuST7+++/Exoayrx58wgMDGTmzJmEhIRw9uxZ3NzccpXfs2cP/fr1Y9q0aXTv3p2lS5fSq1cvDh8+TP369QH44osvmDVrFosXL8bPz48JEyYQEhJCREQElpYFn+FIURRupGQScT2ZiOuJJKemMr6zny455mSwbfdeLsYmYEEW9cnBziwHPycz/JzM8HUwQdl/EVXGbUi/TYM7t3UzH925revBeycB0uPzH4yZNTh4g0OVu4+7Pzve/de+soxFFUKIUqJSFOUhN+ZKR2BgIM2bN2f27NkAaLVavL29GTVqFOPGjctVvm/fvqSlpbFmzRr9vpYtW9K4cWPmzZuHoih4eXnx7rvvMmbMGACSkpJwd3dn0aJFvPTSS4+NKTk5GQcHBy5NrImdaQ5qTSamSjYWZGOm0jz2+YWiNgM7jwcenrp/bf+zbeUkV6NCCFFC7uWDpKQk7O3tH1u+TF3JZmVlcejQIcaPH6/fp1arCQ4OJjw8PM/nhIeHExoaarAvJCSEVatWARAVFUVsbCzBwcH64w4ODgQGBhIeHp5nks3MzCQz8/69zKSkJABMUq5jaqFLaFrgzt3HfSowMQdTCzCxAFNz3b/39pla6K4qze3AyhEsHcHaCSyddNtWTrqHjZtuW61+9BuWA6SkPLqMEEKIIpOcrJsfIL/Xp2Uqyd66dQuNRoO7u2EvWnd3d86cOZPnc2JjY/MsHxsbqz9+b9/DyvzXtGnTmDJlSq793t+k5u+FCCGEeKKlpKTg4PD41cHKVJItK8aPH29wdZyYmIiPjw/R0dH5elOFKA3Jycl4e3tz5cqVfDVjVRRP4vtSXl9TWY87P/EpikJKSgpeXl75qrNMJVkXFxdMTEyIi4sz2B8XF4eHh0eez/Hw8Hhk+Xv/xsXF4enpaVCmcePGedZpYWGBhYVFrv0ODg5l8oMhxIPs7e3lc5qHJ/F9Ka+vqazH/bj4CnKx9ZibfiXL3NycZs2aERYWpt+n1WoJCwsjKCgoz+cEBQUZlAfYtGmTvryfnx8eHh4GZZKTk9m3b99D6xRCCCGKQpm6kgUIDQ1l0KBBBAQE0KJFC2bOnElaWhpDhgwBYODAgVSuXJlp06YBMHr0aNq1a8dXX33FM888w7Jlyzh48CA//PADACqVirfffptPPvkEf39//RAeLy8vevXqVVovUwghRAVQ5pJs3759uXnzJhMnTiQ2NpbGjRuzfv16fcel6Oho1A/0um3VqhVLly7lo48+4oMPPsDf359Vq1bpx8gCvP/++6SlpTFs2DASExN56qmnWL9+fb7HyFpYWDBp0qQ8m5CFKCvkc5q3J/F9Ka+vqazHXRzxlblxskIIIcSTokzdkxVCCCGeJJJkhRBCiGIiSVYIIYQoJpJkhRBCiGIiSfauOXPm4Ovri6WlJYGBgezfv/+R5VesWEHt2rWxtLSkQYMGrF27toQiFRVVQT6jp06donfv3vj6+qJSqZg5c2bJBVoKCvLezJ8/nzZt2uDk5ISTkxPBwcGP/f9eGgrymv766y8CAgJwdHTExsaGxo0b8/PPP5dgtPcV9Lv0nmXLlqFSqYp1aGVBYlu0aBEqlcrgUZhV21CEsmzZMsXc3FxZsGCBcurUKeW1115THB0dlbi4uDzL7969WzExMVG++OILJSIiQvnoo48UMzMz5cSJEyUcuagoCvoZ3b9/vzJmzBjlt99+Uzw8PJRvvvmmZAMuQQV9b15++WVlzpw5ypEjR5TTp08rgwcPVhwcHJSrV6+WcOQPV9DXtHXrVuWvv/5SIiIilAsXLigzZ85UTExMlPXr15fpuO+JiopSKleurLRp00bp2bNnmYht4cKFir29vRITE6N/xMbGFvi8kmQVRWnRooUyYsQI/bZGo1G8vLyUadOm5Vm+T58+yjPPPGOwLzAwUHn99deLNU5RcRX0M/ogHx+fJzrJGvPeKIqi5OTkKHZ2dsrixYuLK8QCM/Y1KYqiNGnSRPnoo4+KI7yHKkzcOTk5SqtWrZQff/xRGTRoULEl2YLGtnDhQsXBwcHo81b45uJ7y+s9uBRefpbXe7A86JbXe1h5IYxRmM9oRVEU7016ejrZ2dk4OzsXV5gFYuxrUhSFsLAwzp49S9u2bYszVAOFjfvjjz/Gzc2NV155pczFlpqaio+PD97e3vTs2ZNTp04V+NwVPsk+anm9hy2F97jl9YQoSoX5jFYURfHejB07Fi8vr1x/OJeWwr6mpKQkbG1tMTc355lnnuG7776jc+fOxR2uXmHi3rVrFz/99BPz588vc7HVqlWLBQsW8Pfff/PLL7+g1Wpp1aoVV69eLdC5y9y0ikIIUVKmT5/OsmXL2LZtW+E6tZQhdnZ2HD16lNTUVMLCwggNDaVatWq0b9++tEPLU0pKCgMGDGD+/Pm4uLiUdji5BAUFGSwi06pVK+rUqcP333/P1KlT811PhU+yxbG8nhBFqTCf0YrCmPdmxowZTJ8+nc2bN9OwYcPiDLNACvua1Go1NWrUAKBx48acPn2aadOmlViSLWjckZGRXLp0iR49euj3abVaAExNTTl79izVq1cvldjyYmZmRpMmTbhw4UKBzl3hm4uLY3k9IYpSYT6jFUVh35svvviCqVOnsn79egICAkoi1Hwrqt+3VqslMzOzOELMU0Hjrl27NidOnODo0aP6x7PPPkuHDh04evQo3t7epRZbXjQaDSdOnDBYlzxfjO469QRYtmyZYmFhoSxatEiJiIhQhg0bpjg6Ouq7aw8YMEAZN26cvvzu3bsVU1NTZcaMGcrp06eVSZMmyRAeUawK+hnNzMxUjhw5ohw5ckTx9PRUxowZoxw5ckQ5f/58ab2EYlPQ92b69OmKubm58scffxgMz0hJSSmtl5BLQV/TZ599pmzcuFGJjIxUIiIilBkzZiimpqbK/Pnzy3Tc/1WcvYsLGtuUKVOUDRs2KJGRkcqhQ4eUl156SbG0tFROnTpVoPNKkr3ru+++U6pWraqYm5srLVq0UPbu3as/1q5dO2XQoEEG5ZcvX67UrFlTMTc3V+rVq6f8+++/JRyxqGgK8hmNiopSgFyPdu3alXzgJaAg742Pj0+e782kSZNKPvBHKMhr+vDDD5UaNWoolpaWipOTkxIUFKQsW7asFKIu+Hfpg4ozyRY0trfffltf1t3dXXn66aeVw4cPF/icstSdEEIIUUwq/D1ZIYQQorhIkhVCCCGKiSRZIYQQophIkhVCCCGKiSRZIYQQophIkhVCCCGKiSRZIYQQophIkhVC6A0ePJhevXqVdhhCPDEq/AIBQlQUKpXqkccnTZrEt99+i8xPI0TRkSQrRAURExOj//n3339n4sSJnD17Vr/P1tYWW1vb0ghNiCeWNBcLUUF4eHjoHw4ODqhUKoN9tra2uZqL27dvz6hRo3j77bdxcnLC3d2d+fPnk5aWxpAhQ7Czs6NGjRqsW7fO4FwnT56kW7du2Nra4u7uzoABA7h161YJv2IhSp8kWSHEIy1evBgXFxf279/PqFGjGD58OC+++CKtWrXi8OHDdOnShQEDBpCeng5AYmIiHTt2pEmTJhw8eJD169cTFxdHnz59SvmVCFHyJMkKIR6pUaNGfPTRR/j7+zN+/HgsLS1xcXHhtddew9/fn4kTJxIfH8/x48cBmD17Nk2aNOGzzz6jdu3aNGnShAULFrB161bOnTtXyq9GiJIl92SFEI/UsGFD/c8mJiZUqlSJBg0a6Pe5u7sDcOPGDQCOHTvG1q1b87y/GxkZSc2aNYs5YiHKDkmyQohHMjMzM9hWqVQG++71WtZqtQCkpqbSo0cPPv/881x1eXp6FmOkQpQ9kmSFEEWqadOm/Pnnn/j6+mJqKl8xomKTe7JCiCI1YsQIEhIS6NevHwcOHCAyMpINGzYwZMgQNBpNaYcnRImSJCuEKFJeXl7s3r0bjUZDly5daNCgAW+//TaOjo6o1fKVIyoWlSLTuwghhBDFQm6Y5INWq+X69evY2dk9dmo6IYQQTy5FUUhJScHLyytfLTOSZPPh+vXreHt7l3YYQgghyogrV65QpUqVx5aTJJsPdnZ2gO5Ntbe3L+VohBBClJbk5GS8vb31eeFxJMnmw70mYnt7e0myQggh8n3rUJJsASzeE4WVje6+rApQqbj7r0r/MwbHVA+U0W3z4HPu7X+gLICZifruQ4W5iRrTuz+bmagxN1Vjqr7/s42FKdZmJqjVcq9YCCHKGkmyBfDlhnOoLaxLO4xcVCqwMTfF1sIUW0tTbCxMsbPQbdtYmGJnef+Ys7U5zjbmONuaE5+aRcT1ZCJikkjL1OBmb4GHvSWeDpa421vi4WCJh70llWwtMJEkLoQQBSZJtgCeaeCBubUtigIKul5mCoACCopu/4M/o9vGYFvR77+3Dfefp9WCRquQrdWSrdGSnaOQrdGSpdGSo7n/c7ZGS1aOFu3dc6Zm5pCamQPJRf+6TdQq3OwsqOxoRWUnK3ycralX2YGGVRzwsLeUHtdCCPEQMk42H5KTk3FwcCApKalM3ZNVFIXMHC0pGboEm5aZY/hzZg6pGbqfUzNzSM7I5nZaFvFpWcSnZmFnaUpdL3vqetrjZG1OXEoGcUkZxCZnEHv335spmWgf8QlxtbOgYWUHGlZxpGEVXeKtZGtRcm+CEI+g0WjIzs4u7TBEOWJiYoKpqelDLx4Kmg8kyeZDWU2yJSFHo+VWahYxSXe4lniHa7fvEHkzleNXkzh/IxVNHhm4sqMVDas40KCKA42qOFLfywEHa7M8ahei+KSmpnL16lXkK04UlLW1NZ6enpibm+c6Jkm2GFTkJPsod7I0RMQkcfzqvUciF2+lkdcnqpKNOTXcbGnq40RzXyeaVXXGwdqM9KwcYpIy8HG2xtREptwTRUOj0XD+/Hmsra1xdXWVWxoiXxRFISsri5s3b6LRaPD398814URB84HckxWFZmVuQjMfZ5r5OOv3pWRkc/JaMsevJnL8mi7xXkm4o2uijkpgX1QCc++W9XKwJDY5A60CdpamtK7uQpuaLrT1d8Xbuex1MBPlR3Z2Noqi4OrqipWVVWmHI8oRKysrzMzMuHz5MllZWVhaWhpVnyRZUaTsLM0Iql6JoOqV9PtSM3OIupnG6ZhkDl5O4OCl21y8lcb1pAwAzE3UpGTksP5ULOtPxQLg52JDG39dwm1ZvRK2FvJRFQUnV7CiMIpyIQv55hLFztbClAZ379H2aa6bnvJWaiaRN1Lxc7Ghkq0FJ64lsePcTXaev8nh6ESibqURdSuNJeGXMTNREeDjTLtarrSv5Uotd5lDWghRPsg92XyQe7IlKzkjm/DIeHaev8mOc7eITkg3OO5ub0G7mq60r+VG6xouOFhJpyphKCMjg6ioKPz8/Ixu7hMVz6M+P3JPVpR79pZmhNTzIKSeBwCXbqWx/dxNtp29QfjFeOKSM1l+8CrLD17FRK2iaVVHfdKt62kvs18JUY5NmDCBuLg4fvjhhyKpr3379jRu3JiZM2cC0LJlS9577z169+5dJPU/jlzJ5oNcyZYdGdka9kclsP3cTbafu8mFG6kGx11szWnr70q7Wq608XfF2SZ3F3zx5CuvV7KDBw9m8eLFAJiZmVG1alUGDhzIBx98gKnpk39NFBsbS82aNTlx4gQ+Pj5FUmdCQgJmZmb6Cf3XrFnDO++8w9mzZx9671WuZEWFZWlmQtuarrSt6coE4EpCOjvO32T72ZvsvnCLW6lZ/HXkGn8duYZKBY2q6K5y29VypVEVR5keUpR5Xbt2ZeHChWRmZrJ27VpGjBiBmZkZ48ePz1U2Kysrz7GcpcmYmH788UdatWpVZAkWwNnZ2WC7W7duvPrqq6xbt45nnnmmyM7zMDIwUZRr3s7W9A/04YeBARyZ2IXfXmvJ6+2qUdvDDkWBo1cS+TbsPM//bw/NPtnE6GVH+PvoNRLTs0o7dFGCFEUhPSunVB4FbSy0sLDAw8MDHx8fhg8fTnBwMKtXrwZ0V7q9evXi008/xcvLi1q1agG6ZTj79OmDo6Mjzs7O9OzZk0uXLunr3LZtGy1atMDGxgZHR0dat27N5cuXATh27BgdOnTAzs4Oe3t7mjVrxsGDBwGYPHkyjRs3Nohv5syZ+Pr66rcLG1Neli1bRo8ePQz2tW/fnpEjRzJy5EgcHBxwcXFhwoQJBu/r//73P/z9/bG0tMTd3Z0XXnjB4Plvv/22ftvExISnn36aZcuWPTKWoiJXsuKJYW6q1g8fGt+tDrFJGey426y88/xNEtOz+fvodf4+eh0TtYpAP2e61fegSz0P3O3LT5OiKLg72RrqTtxQKueO+DgEa/PCf9VaWVkRHx+v3w4LC8Pe3p5NmzYBujHBISEhBAUFsXPnTkxNTfnkk0/o2rUrx48fR61W06tXL1577TV+++03srKy2L9/v76Hfv/+/WnSpAlz587FxMSEo0ePYmZWsM6EBY0pryvdhIQEIiIiCAgIyHVs8eLFvPLKK+zfv5+DBw8ybNgwqlatymuvvcbBgwd56623+Pnnn2nVqhUJCQns3LnzkfG2aNGC6dOnF+g1FpYkWfHE8nCwpE9zb/o09yZHo+XIlUS2nLnBltM3OBuXwp7IePZExjPh71M0qepIl7oedKnnTnVX29IOXQgURSEsLIwNGzYwatQo/X4bGxt+/PFHfaL65Zdf0Gq1/Pjjj/rEuXDhQhwdHdm2bRsBAQEkJSXRvXt3qlevDkCdOnX09UVHR/Pee+9Ru3ZtAPz9/Qsca0Fj6tKlS646oqOjURQFLy+vXMe8vb355ptvUKlU1KpVixMnTvDNN9/w2muvER0djY2NDd27d8fOzg4fHx+aNGnyyHi9vLy4cuUKWq22SMfE5kWSrKgQTE3UNPd1prmvM2O71iY6Pp31p2JYdzKWI9GJ+sfn689Q3dWGLnd7Nzes7CC9lZ8AVmYmRHwcUmrnLog1a9Zga2tLdnY2Wq2Wl19+mcmTJ+uPN2jQwOBK8NixY1y4cEHfseeejIwMIiMj6dKlC4MHDyYkJITOnTsTHBxMnz598PT0BCA0NJRXX32Vn3/+meDgYF588UV9Ms6vgsaUlzt37gDk2VGtZcuWBmPjg4KC+Oqrr9BoNHTu3BkfHx+qVatG165d6dq1K8899xzW1g+fNc7KygqtVktmZmaxzwgmSVZUSFUrWTOsbXWGta1OXHIGmyLi2BgRR3jkLSJvpjF3WyRzt0Xibm9B57rudK3nSctqzjK/cjmlUqmMarItSR06dGDu3LmYm5vj5eWVq1exjY2NwXZqairNmjXj119/zVWXq6sroLuKfOutt1i/fj2///47H330EZs2baJly5ZMnjyZl19+mX///Zd169YxadIkli1bxnPPPYdarc51TzmvVY0KE9N/ubi4AHD79u2HlsmLnZ0dhw8fZtu2bWzcuJGJEycyefJkDhw4gKOjY57PSUhIwMbGpkSm3CwfnzohipG7vSX/19KH/2vpQ3JGNtvO3mTDqVi2nblBXHImv+yN5pe90ThZ68bvdmvgSavqlTCThCuKgY2NDTVq1Mh3+aZNm/L777/j5ub2yCElTZo0oUmTJowfP56goCCWLl1Ky5YtAahZsyY1a9bknXfeoV+/fixcuJDnnnsOV1dXYmNjURRFfyV59OjRIovpQdWrV8fe3p6IiAhq1qxpcGzfvn0G23v37sXf3x8TE10rgampKcHBwQQHBzNp0iQcHR3ZsmULzz//fJ7nOnny5GOblIuKfEsI8QB7SzOebeTFnJebcnhiZxYOac5Lzb1xtjHndno2yw5cYdCC/QR8spn3Vhxj65kbZOZoSjtsUYH1798fFxcXevbsyc6dO4mKimLbtm289dZbXL16laioKMaPH094eDiXL19m48aNnD9/njp16nDnzh1GjhzJtm3buHz5Mrt37+bAgQP6e7bt27fn5s2bfPHFF0RGRjJnzhzWrVtndEx5UavVBAcHs2vXrlzHoqOjCQ0N5ezZs/z222989913jB49GtA1r8+aNYujR49y+fJllixZglar1fdyzsvOnTvzvC9cHORKVoiHsDA1oUMtNzrUcuOTXlr2RSWw9kQMG07Fcis1ixWHrrLi0FVszHVjd4PruNOxthtOMgGGKEHW1tbs2LGDsWPH8vzzz5OSkkLlypXp1KkT9vb23LlzhzNnzrB48WLi4+Px9PRkxIgRvP766+Tk5BAfH8/AgQOJi4vDxcWF559/nilTpgC6DlL/+9//+Oyzz5g6dSq9e/dmzJgxj52N6XExPcyrr77Ka6+9xhdffGHQIWngwIHcuXOHFi1aYGJiwujRoxk2bBgAjo6O/PXXX0yePJmMjAz8/f357bffqFevXp7nuHbtGnv27OGXX34p6FtdKDLjUz7IjE/iQRqtwv67CXdjRCxxyZn6Y2oVBPg607mOO53ruuPrYvOImkRxKa8zPlV0iqIQGBiob7aG3NMiGmvs2LHcvn37kX8oyIxPQpQiE7VKPx53yrP1OHk9iU0RcWyKiONMbAr7oxLYH5XAp2tPU8PNluA67nSu60ZjbyeZcUqIR1CpVPzwww+cOHGi2M7h5uZGaGhosdX/X5JkhTCCWq2iYRVHGlZx5N0utbiSkE7Y6Tg2nY5j38UELtxI5cKNVOZtj8TF1pyOtd0IruNO25quWBZwaIcQFUHjxo1zzTJVlN59991iqzsv0lycD9JcLAoj6U4228/dZHNEHFvP3iAlI0d/zNrchE513HmmgSfta0nCLWrSXCyMIc3FQpQDDla6nsrPNvIiK0fLgUsJ+mbla4l3+OfYdf45dh2bewm3oSft5ApXiCeKJFkhSoC5qZrWNVxoXcOFST3qcuxqEv8ev87aE7FcS7zD6mPXWX3sOrYWpgTXceOZhl608XeRhGskaagThVGUnxtJskKUMJVKRWNvRxp7O/LB03U4eiWRf4/HsPZEDNeTMlh19DqrjuoSbue6uiblNjVdsDCVhJtf9yYpyMrKKpFZfcSTJT09HaDACyXkRe7J5oPckxUlQatVOHIlkbUndAk3JilDf8zOwpTO9dzp3tCTp2q4Ym4q88g8iqIoREdHk52djZeXV7FPAi+eDIqikJ6ezo0bN3B0dNTP7/ygguYDSbL5IElWlDRdwr3NmrtXuA+OxbWzNKVLXQ+6N/SkdQ0XSbgPkZWVRVRUFFqttrRDEeWMo6MjHh4eBosS3FOukqxGo2Hy5Mn88ssvxMbG4uXlxeDBg/noo4/0L05RFCZNmsT8+fNJTEykdevWzJ0712A5poSEBEaNGsU///yDWq2md+/efPvtt9ja3l+y7Pjx44wYMYIDBw7g6urKqFGjeP/99/MVpyRZUZq0WoVD0bf1Tco3Uu4nXHtLU0LqefBqm2rU8rB7RC0Vk1arJSsrq7TDEOWImZmZ/nZDXgqcD5RS9OmnnyqVKlVS1qxZo0RFRSkrVqxQbG1tlW+//VZfZvr06YqDg4OyatUq5dixY8qzzz6r+Pn5KXfu3NGX6dq1q9KoUSNl7969ys6dO5UaNWoo/fr10x9PSkpS3N3dlf79+ysnT55UfvvtN8XKykr5/vvv8xVnUlKSAihJSUlF9+KFKASNRqvsuxivTFx1Qgn4ZJPiM3aN4jN2jeI3bo3y/opjypmY5NIOUYgnWkHzQaleyXbv3h13d3d++ukn/b7evXtjZWXFL7/8ol/A991332XMmDEAJCUl4e7uzqJFi3jppZc4ffo0devW5cCBAwQEBACwfv16nn76aa5evYqXlxdz587lww8/JDY2Vr/m4bhx41i1ahVnzpx5bJxyJSvKIo1W4eClBBbuvsT6U7H6/bU97HgxwJsXmlbBwdr4jhtCiPsKmg9K9WZOq1atCAsL49y5c4Buod9du3bRrVs3AKKiooiNjSU4OFj/HAcHBwIDAwkPDwcgPDwcR0dHfYIFCA4ORq1W65dHCg8Pp23btgaLCoeEhHD27Flu375d7K9TiOJgolYRWK0S8wY048/hrehc1x0zExVnYlOYuiaCFp9tJnT5UfZHJchQFiFKSakO4Rk3bhzJycnUrl0bExMTNBoNn376Kf379wcgNlb317m7u7vB89zd3fXHYmNjcXNzMzhuamqKs7OzQRk/P79cddw75uTkZHAsMzOTzMz7972Sk5ONfalCFKtmPk7MHxhAUno2q49d49d90ZyJTeGvw9f46/A1qrnaMLxddXo2riwdpYQoQaX6v2358uX8+uuvLF26lMOHD7N48WJmzJjB4sWLSzMspk2bhoODg/7h7e1dqvEIkV8O1mYMCPJl3eg2/Dm8FX0CqmBlZsLFm2m898dxWk4LY/LqU5y4miRXt0KUgFK9kn3vvfcYN24cL730EgANGjTg8uXLTJs2jUGDBuHh4QFAXFycwXiluLg4/QTSHh4e3Lhxw6DenJwcEhIS9M/38PAgLi7OoMy97XtlHjR+/HiDVRqSk5Ml0YpyRaVS0czHiWY+TkzsUY9f917mx11R3EzJZNGeSyzacwl/N1u6N/Ti6QYe+LtLz2QhikOpXsmmp6fnGiRuYmKiH9fm5+eHh4cHYWFh+uPJycns27ePoKAgAIKCgkhMTOTQoUP6Mlu2bEGr1RIYGKgvs2PHDrKzs/VlNm3aRK1atXI1FQNYWFhgb29v8BCivLK1MOX1dtUJH9eRhUOa072hJ+amas7fSOWbzefo/M0OOn+9nXnbI4lLznh8hUKIfCvV3sWDBw9m8+bNfP/999SrV48jR44wbNgwhg4dyueffw7A559/zvTp01m8eDF+fn5MmDCB48ePExERoV8doVu3bsTFxTFv3jyys7MZMmQIAQEBLF26FND1SK5VqxZdunRh7NixnDx5kqFDh/LNN98wbNiwx8YpvYvFkybpTjabIuJYdyKGnedvkaXR/WGrVkG7mq70a1GVjrXdMDWR+7dCPKhcTUaRkpLChAkTWLlyJTdu3MDLy4t+/foxceJEfU9g5e5kFD/88AOJiYk89dRT/O9//6NmzZr6ehISEhg5cqTBZBSzZs166GQULi4ujBo1irFjx+YrTkmy4kmWnJHN2uMx/HHoKgcv3+9t725vQZ8Ab/oEeOPtbF2KEQpRdpSrJFteSJIVFcXFm6n8fuAKKw5dJSFNN1OSSgVt/XVXt53quGEmV7eiApMkWwwkyYqKJitHy6aIOH7bH82uC7f0+13tLOgb4M3LgVXxcpTVbUTFI0m2GEiSFRXZ5fg0lh24woqDV7mVen/8uG8la55vWoVXnvLDxkJWzRQVgyTZYiBJVgjI1mjZHBHHwt2XOHg5Ae3dbw4XWwveDvanT4C3THQhnniSZIuBJFkhDCVnZLPl9A2+2XyOy/G6Ba4rO1rxRvvqvNisCpZmssC8eDJJki0GkmSFyFtWjpal+y4ze2ukvinZ3d6Cfi2q8nKLqrjZW5ZyhEIULUmyxUCSrBCPlpGt4fcDV5i3PZKYJN2EFjbmJozs6M/Qp3yxMJUrW/FkkCRbDCTJCpE/mTka/j0ew5Lwyxy9kgjoOkh99ExdOtVxQ6VSlW6AQhhJkmwxkCQrRMFotQorj1xj+voz3EzRNSO3r+XKlGfr4VPJppSjE6LwJMkWA0myQhROamYOs7dc4KddF8nWKJibqhnerjrD21eXzlGiXJIkWwwkyQphnIs3U5m0+hQ7z+smtqjqbM3E7tKELMqfguYDGdQmhCh21VxtWTK0BXNeboqHvSXRCem8uuQgQxYdIOpWWmmHJ0SxkSvZfJArWSGKTlpmDt890IRsZqJi6FN+vNm+Bg5WZqUdnhCPJM3FxUCSrBBF7+LNVKb8E8H2czcBsLc05Y321RnSyg8rc7lfK8omSbLFQJKsEMVDURTCTt/g8/VnOH8jFdAtQvBWJ3/6yjSNogySJFsMJMkKUbw0WoVVR67xzeZzXL19B9B1jnqnsz/PNqqMiVo6R4myoUSTbGZmJhYWFoV9erkhSVaIkpGVo+W3/dF8t+WCfppG30rWPN3Ak8GtfXGzk2kaRekq1iS7bt06li1bxs6dO7ly5QparRYbGxuaNGlCly5dGDJkCF5eXka9gLJIkqwQJSs9K4eFuy8xb3skKRk5AJibqukfWJU329fA1e7J/+NelE3FkmRXrlzJ2LFjSUlJ4emnn6ZFixZ4eXlhZWVFQkICJ0+eZOfOnYSHhzN48GCmTp2Kq6trkbygskCSrBClIzkjm+1nb7JwdxSHoxMBsDIzYXBrX155yg8XW0m2omQVS5INCgrio48+olu3bqjVD++IcO3aNb777jvc3d155513ChZ5GSZJVojSpSgKuy/E8+XGsxy7Oyeyuama55tU5q1O/ng5WpVugKLCkI5PxUCSrBBlg6IobD59g9lbznPsahKgS7YDW/owokMNnGzMSzlC8aSTJFsMJMkKUfYcvJTAjI1n2XsxAQBHazPeC6lF57ru0kFKFJtiSbKhoaFMnToVGxsbQkNDH1n266+/zn+05YQkWSHKJkVR2HH+FtPWnuZMbAoAJmoVA1r68E5wTRysZQYpUbQKmg9M81PpkSNHyM7O1v/8MDLRtxCiJKlUKtrVdKVV9Uos3nOJPw5d5UxsCov2XOLvo9d4L6Q2/Vp4y3eTKDXSXJwPciUrRPmx6/wtPl5zinNxuhmkejX2Ysqz9eWqVhSJEl+FJzk5mVWrVnHmzBljqxJCCKM95e/C2rfa8MHTtVGrYNXR6wR/s50Np2JLOzRRARU4yfbp04fZs2cDcOfOHQICAujTpw8NGjTgzz//LPIAhRCioExN1AxrW53lrwdRzdWGmymZvP7zId767QhJd7JLOzxRgRQ4ye7YsYM2bdoAukkqFEUhMTGRWbNm8cknnxR5gEIIUVgBvs6sfasNIzpUx0StYvWx63SbuYO9F+NLOzRRQRQ4ySYlJeHs7AzA+vXr6d27N9bW1jzzzDOcP3++yAMUQghjWJqZ8F5Ibf54IwifStZcT8rg5fl7+XbzeTRa6ZIiileBk6y3tzfh4eGkpaWxfv16unTpAsDt27extJSxaUKIsqlJVSfWvtWGF5pVQavAN5vPMeCnfdxIzijt0MQTrMBJ9u2336Z///5UqVIFLy8v2rdvD+iakRs0aFDU8QkhRJGxsTBlxouN+LpPI6zNTdgTGc/Ts3ay8/zN0g5NPKEKNYTn4MGDXLlyhc6dO2NrawvAv//+i6OjI61bty7yIEubDOER4slz4UYqI5ce5kxsCioVvNm+Ou8E18TURBaKFw9XotMq3nvqkz7QW5KsEE+mjGwNH6+JYOm+aABa16jEd/2a4ixzIIuHKJFxsj/99BP169fH0tISS0tL6tevz48//liYqrh27Rr/93//R6VKlbCysqJBgwYcPHhQf1xRFCZOnIinpydWVlYEBwfn6mCVkJBA//79sbe3x9HRkVdeeYXU1FSDMsePH6dNmzZYWlri7e3NF198Uah4hRBPDkszEz57rgHf9WuCtbkJuy/E88ysney5cKu0QxNPiAIn2YkTJzJ69Gh69OjBihUrWLFiBT169OCdd95h4sSJBarr9u3btG7dGjMzM9atW0dERARfffUVTk5O+jJffPEFs2bNYt68eezbtw8bGxtCQkLIyLjfWaF///6cOnWKTZs2sWbNGnbs2MGwYcP0x5OTk+nSpQs+Pj4cOnSIL7/8ksmTJ/PDDz8U9OULIZ5APRp5sfLN1vi52BCTlMHLP+7jkzURZGRrSjs0Ud4pBeTi4qIsXbo01/6lS5cqlSpVKlBdY8eOVZ566qmHHtdqtYqHh4fy5Zdf6vclJiYqFhYWym+//aYoiqJEREQogHLgwAF9mXXr1ikqlUq5du2aoiiK8r///U9xcnJSMjMzDc5dq1atfMWZlJSkAEpSUlKBXp8QonxJzchWxv15XPEZu0bxGbtG6TVnl3IjOaO0wxJlSEHzQYGvZLOzswkICMi1v1mzZuTk5BSortWrVxMQEMCLL76Im5sbTZo0Yf78+frjUVFRxMbGEhwcrN/n4OBAYGAg4eHhAISHh+Po6GgQU3BwMGq1mn379unLtG3bFnPz+/dZQkJCOHv2LLdv3y5QzEKIJ5eNhSnTnm/A/IEB2FuaciQ6kadn7WT9yZjSDk2UUwVOsgMGDGDu3Lm59v/www/079+/QHVdvHiRuXPn4u/vz4YNGxg+fDhvvfUWixcvBiA2VjfXqLu7u8Hz3N3d9cdiY2Nxc3MzOG5qaoqzs7NBmbzqePAcD8rMzCQ5OdngIYSoODrXdWfViNbUcLPlZkomb/xymDd+PkScjKkVBZSvpe7+66effmLjxo20bNkSgH379hEdHc3AgQMN1pt93NqyWq2WgIAAPvvsMwCaNGnCyZMnmTdvHoMGDSpMaEVi2rRpTJkypdTOL4QofdVcbVkz6ilmb7nAvO2RrD8Vy+7IW4zu5E/f5t7YWcqqPuLxCnwle/LkSZo2bYqrqyuRkZFERkbi4uJC06ZNOXnyJEeOHOHIkSMcPXr0sXV5enpSt25dg3116tQhOlrXnd7DwwOAuLg4gzJxcXH6Yx4eHty4ccPgeE5ODgkJCQZl8qrjwXM8aPz48SQlJekfV65ceexrEUI8eSzNTBgTUovVI5+iYRUHUjJy+OTf03SYsY1l+6NlWkbxWAW+kt26dWuRnbx169acPXvWYN+5c+fw8fEBwM/PDw8PD8LCwmjcuDGg6ym8b98+hg8fDkBQUBCJiYkcOnSIZs2aAbBlyxa0Wi2BgYH6Mh9++CHZ2dmYmen++ty0aRO1atUy6Ml8j4WFBRYWFkX2OoUQ5VtdL3tWvtma3/ZH89OuKKJupTHurxMsCb/MhO51CapeqbRDFGVVMXfEeqT9+/crpqamyqeffqqcP39e+fXXXxVra2vll19+0ZeZPn264ujoqPz999/K8ePHlZ49eyp+fn7KnTt39GW6du2qNGnSRNm3b5+ya9cuxd/fX+nXr5/+eGJiouLu7q4MGDBAOXnypLJs2TLF2tpa+f777/MVp/QuFkLck5mtUX7ceVFpMGm9vhfy60sOKlcS0ko7NFECCpoP8pVkX3/9deXKlSv5qnDZsmUGSfJx/vnnH6V+/fqKhYWFUrt2beWHH34wOK7VapUJEyYo7u7uioWFhdKpUyfl7NmzBmXi4+OVfv36Kba2toq9vb0yZMgQJSUlxaDMsWPHlKeeekqxsLBQKleurEyfPj3fMUqSFUL8V3xqpjJh1Qml2vh/FZ+xa5TaH61T5m67oGTlaEo7NFGMCpoP8jWt4oQJE5g1axatW7emR48eBAQE4OXlhaWlJbdv3yYiIoJdu3axbNkyvLy8+OGHH2jYsGFxX4SXGJlWUQjxMOfiUvho1Un2RyUAUNPdlk+fa0BzX+dSjkwUh2KbuzguLo4ff/yRZcuWERERYXDMzs6O4OBgXn31Vbp27Vq4yMswSbJCiEdRFIU/D1/js7WnSUjLAnSzSA1rU40GVRxKOTpRlEpkgYDbt28THR3NnTt3cHFxoXr16k/0IgGSZIUQ+ZGYnsX0dWdYduD+iIQ2/i6M7Vqb+pUl2T4JSnQVnopCkqwQoiBOXktiwa4oVh+7Ts7dYT7PNPRkRPsa1PWS75DyTJJsMZAkK4QojOj4dL7edJa/j13n3jdtu5qujA72p2nV3MMHRdknSbYYSJIVQhjjdEwyc7ZeYO2JGO7NX/Fck8qM7VobDwfL0g1OFIgk2WIgSVYIURQux6fx3ZYL/Hn4KooCVmYmvNm+Oq+2qYaVuUlphyfyQZJsMZAkK4QoSsevJvLxPxEcvKxbBczNzoK3OvnzQrMqWJpJsi3LSiTJ5uTksG3bNiIjI3n55Zexs7Pj+vXr2NvbY2trW6jAyzJJskKIoqYoCv8cj+GL9We4evsOAPaWpjzftAovNKtCPS/7J3rURnlV7En28uXLdO3alejoaDIzMzl37hzVqlVj9OjRZGZmMm/evEIHX1ZJkhVCFJfMHA3L9l9h/s6L+mQLUMvdjpB67rSv7UajKo6YqCXhlgXFnmR79eqFnZ0dP/30E5UqVeLYsWNUq1aNbdu28dprr3H+/PlCB19WSZIVQhQ3rVZh+/mb/HHoKpsi4sjK0eqPOduY07G2G8F13GhZrRKO1ualGGnFVtB8UOBVeHbu3MmePXswNzf8Jfv6+nLt2rWCVieEEAJQq1V0qOVGh1puJN3JZuOpWLadu8mOczdJSMvij0NX+ePQVQBqe9jRws+ZAF9nqjhZ4WZngYuthdzPLYMKnGS1Wi0ajSbX/qtXr2JnZ1ckQQkhREXmYGXGiwHevBjgTbZGy4FLCWyOuMH2czeIvJnGmdgUzsSmsCT8ssHz7CxMqWRrTiVbC1xszXGx1SVfFzsLXB/Ydre3lN7MJaTAzcV9+/bFwcGBH374ATs7O44fP46rqys9e/akatWqLFy4sLhiLTXSXCyEKCtupWZyICqBfVEJHLuayI3kTG6mZJKl0T7+yXeZm6ip5mqDnaUptham2FqaYWthen/77sPGwhRbS1NsLUywtTDDxsIEOwszVGpISM3SHbcwxdJMXWE6aRX7PdkrV67QtWtXFEXh/PnzBAQEcP78eVxcXNixYwdubm6FDr6skiQrhCjLFEUh+U4O8WmZ3ErNIj41k/9v796DojrvN4A/y7K7iLCAUBdJEE0xwQsKhJhiJlFbFA1NgmEa0xiCJtpcViux0cbLyCSZilrbmCYmOqaKtqYoia3jpSJBtKmQ0YBYRIPGatHqgrZlF1SWvby/P5Dzc4UYDnBY2H0+M2fgvOfdc77nuMPjuV9rtOJqYzOuNVpxrcGKq7farjU046at7dHIrlD7qOCvVUPn6wON2gfa1p+3ftdKbaqWcV81NGrV//e//TO39XdtU7WZr0bt47JMX/WtoBdAa7C1JlyTzYHL5psw37ChocmOBqsdDU023Gh24EazHWpVy/zjB4fg8dhB37quPXYLz/bt23HixAk0NjYiISEBM2bMQL9+/eTOqk9gyBKRJ/mmrgFXzE1ovBU2jU12NFpbhoYmGxqtDly/o/26taVv6wVZOl8fWO0d33vuK559KBIr07/9Va2KhqzNZkNMTAz27NmD4cOHd/RjfR5DloioRbPdiSa7A4E6XwgB3LC1BPKNZgea7U7YHE5Y7U7pd+nnrXaXtlv9mh3i1k8HbHaB5lv9m++YT/NtP20u01s+b3M60XrQuvXwdeu4r1qFiKB+GNBf23JY/LZD5H4aNSAErHYnRt8bjNTR3bcnK+vCJ41Gg6amJjkfISIiD9J6GBcAVCpI52+pfT5yP2A0GrFq1SrY7XYl6iEiIvIYsv/7cezYMRQVFeHAgQOIjY1F//79Xabv3Lmz24ojIiLqy2SHbHBwMNLT05WohYiIyKPIDllPvA+WiIhICbLPyRIREVHHyN6THTp06F2f7PHPf/6zSwURERF5Ctkhm5WV5TJus9lw/Phx7N+/HwsXLuyuuoiIiPo82SE7f/78dtvXrVuHr776qssFEREReYpuOyc7depUfPbZZ901OyIioj6v20L2008/xYABA7prdkRERH2e7MPF8fHxLhc+CSFgMplw9epVfPjhh91aHBERUV8mO2TT0tJcxn18fPC9730PEyZMQExMTHfVRURE1Od16lV33oZv4SEiIkB+Hsg+J1teXo7KykppfNeuXUhLS8OSJUvQ3Nwsd3ZEREQeS3bIvvzyyzhz5gyAlgdPTJ8+Hf7+/sjPz8eiRYu6vUAiIqK+SnbInjlzBnFxcQCA/Px8jB8/Hp988glyc3O7dAvPypUroVKpXB520dTUBKPRiNDQUAQEBCA9PR21tbUun6upqUFqair8/f0xcOBALFy4sM1r+A4dOoSEhATodDpER0cjNze303USERF1lOyQFULA6XQCAD7//HM8/vjjAIDIyEhcu3atU0UcO3YMGzZswOjRo13aX3/9dezevRv5+fk4fPgwLl++jKefflqa7nA4kJqaiubmZpSUlGDLli3Izc3F8uXLpT7nz59HamoqJk6ciIqKCmRlZWH27NkoKCjoVK1EREQdJmSaOHGieOGFF8TWrVuFRqMRZ8+eFUIIcejQIREVFSV3dqKhoUEMGzZMFBYWivHjx4v58+cLIYSor68XGo1G5OfnS31Pnz4tAIjS0lIhhBD79u0TPj4+wmQySX0++ugjodfrhdVqFUIIsWjRIjFy5EiXZU6fPl2kpKR0uEaz2SwACLPZLHv9iIjIc8jNA9l7smvXrkV5eTnmzp2LpUuXIjo6GkDLwyjGjRsnO+SNRiNSU1ORnJzs0l5WVgabzebSHhMTg8GDB6O0tBQAUFpaitjYWBgMBqlPSkoKLBYLqqqqpD53zjslJUWaBxERkVJk3yc7evRol6uLW/3617+GWq2WNa+8vDyUl5fj2LFjbaaZTCZotVoEBwe7tBsMBphMJqnP7QHbOr112t36WCwW3Lx5E/369WuzbKvVCqvVKo1bLBZZ60VERAR04pzsxYsXcenSJWn86NGjyMrKwtatW6HRaGTNZ/78+di2bRv8/PzklqGonJwcBAUFSUNkZKS7SyIioj5Idsg+99xzKC4uBtCylzhp0iQcPXoUS5cuxdtvv93h+ZSVlaGurg4JCQnw9fWFr68vDh8+jN/97nfw9fWFwWBAc3Mz6uvrXT5XW1uL8PBwAEB4eHibq41bx7+rj16vb3cvFgAWL14Ms9ksDRcvXuzwehEREbWSHbInT57E2LFjAQA7duzAqFGjUFJSgm3btsm6NeZHP/oRKisrUVFRIQ2JiYmYMWOG9LtGo0FRUZH0merqatTU1CApKQkAkJSUhMrKStTV1Ul9CgsLodfrMWLECKnP7fNo7dM6j/bodDro9XqXgYiISC7Z52RtNht0Oh2Allt4nnzySQAtFyVduXKlw/MJDAzEqFGjXNr69++P0NBQqf2ll17CggULMGDAAOj1esybNw9JSUn4wQ9+AACYPHkyRowYgYyMDKxevRomkwnLli2D0WiUanzllVfwwQcfYNGiRXjxxRdx8OBB7NixA3v37pW76kRERLLI3pMdOXIk1q9fjy+++AKFhYWYMmUKAODy5csIDQ3t1uLeffdd/PjHP0Z6ejoee+wxhIeHY+fOndJ0tVqNPXv2QK1WIykpCc8//zxeeOEFl8PWQ4cOxd69e1FYWIgxY8bgN7/5DT7++GOkpKR0a61ERER3kv2CgEOHDmHatGmwWCzIzMzEpk2bAABLlizB119/7RKCnoIvCCAiIkB+HnTqLTwOhwMWiwUhISFS24ULF6RHG3oahiwREQE98BYeoOXRimVlZdiwYQMaGhoAAFqtFv7+/p2ZHRERkUeSfeHTv/71L0yZMgU1NTWwWq2YNGkSAgMDsWrVKlitVqxfv16JOomIiPoc2Xuy8+fPR2JiIv73v/+53Gc6bdq0NrfKEBEReTPZe7JffPEFSkpKoNVqXdqHDBmCf//7391WGBERUV8ne0/W6XTC4XC0ab906RICAwO7pSgiIiJPIDtkJ0+ejLVr10rjKpUKjY2NyM7Olt4tS0RERJ24hefixYuYMmUKhBA4e/YsEhMTcfbsWYSFheFvf/sbb+EhIiKP1SP3ydrtdmzfvh0nTpxAY2MjEhISMGPGjG994H5fx5AlIiJA4ZC12WyIiYnBnj17MHz48C4V2pcwZImICFD4YRQajQZNTU2dLo6IiMibyL7wyWg0YtWqVbDb7UrUQ0RE5DFk3yd77NgxFBUV4cCBA4iNjUX//v1dpnviCwKIiIg6Q3bIBgcHIz09XYlaiIiIPIrskN28ebMSdRAREXmcDp+TdTqdWLVqFR555BE89NBDePPNN3Hz5k0layMiIurTOhyyv/rVr7BkyRIEBATgnnvuwXvvvQej0ahkbURERH1ah0N269at+PDDD1FQUIC//OUv2L17N7Zt2wan06lkfURERH1Wh0O2pqbG5dnEycnJUKlUuHz5siKFERER9XUdDlm73Q4/Pz+XNo1GA5vN1u1FEREReYIOX10shMDMmTOh0+mktqamJrzyyisu98ryPlkiIqIWHQ7ZzMzMNm3PP/98txZDRETkSTocsrw/loiISB7Zzy4mIiKijmHIEhERKYQhS0REpBCGLBERkUIYskRERAphyBIRESmEIUtERKQQhiwREZFCGLJEREQKcWvI5uTk4KGHHkJgYCAGDhyItLQ0VFdXu/RpamqC0WhEaGgoAgICkJ6ejtraWpc+NTU1SE1Nhb+/PwYOHIiFCxfCbre79Dl06BASEhKg0+kQHR2N3NxcpVePiIi8nFtD9vDhwzAajfjyyy9RWFgIm82GyZMn4/r161Kf119/Hbt370Z+fj4OHz6My5cv4+mnn5amOxwOpKamorm5GSUlJdiyZQtyc3OxfPlyqc/58+eRmpqKiRMnoqKiAllZWZg9ezYKCgp6dH2JiMjLiF6krq5OABCHDx8WQghRX18vNBqNyM/Pl/qcPn1aABClpaVCCCH27dsnfHx8hMlkkvp89NFHQq/XC6vVKoQQYtGiRWLkyJEuy5o+fbpISUnpUF1ms1kAEGazuUvrR0REfZvcPOhV52TNZjMAYMCAAQCAsrIy2Gw2JCcnS31iYmIwePBglJaWAgBKS0sRGxsLg8Eg9UlJSYHFYkFVVZXU5/Z5tPZpnQcREZESOvwWHqU5nU5kZWXhkUcewahRowAAJpMJWq0WwcHBLn0NBgNMJpPU5/aAbZ3eOu1ufSwWC27evIl+/fq5TLNarbBardK4xWLp+goSEZHX6TV7skajESdPnkReXp67S0FOTg6CgoKkITIy0t0lERFRH9QrQnbu3LnYs2cPiouLce+990rt4eHhaG5uRn19vUv/2tpahIeHS33uvNq4dfy7+uj1+jZ7sQCwePFimM1mabh48WKX15GIiLyPW0NWCIG5c+fiz3/+Mw4ePIihQ4e6TH/wwQeh0WhQVFQktVVXV6OmpgZJSUkAgKSkJFRWVqKurk7qU1hYCL1ejxEjRkh9bp9Ha5/WedxJp9NBr9e7DERERHKphBDCXQt/7bXX8Mknn2DXrl144IEHpPagoCBpD/PVV1/Fvn37kJubC71ej3nz5gEASkpKALTcwhMXF4eIiAisXr0aJpMJGRkZmD17NlasWAGg5RaeUaNGwWg04sUXX8TBgwfx85//HHv37kVKSsp31mmxWBAUFASz2czAJSLyYrLzQMlLnb8LgHaHzZs3S31u3rwpXnvtNRESEiL8/f3FtGnTxJUrV1zmc+HCBTF16lTRr18/ERYWJn7xi18Im83m0qe4uFjExcUJrVYr7rvvPpdlfBfewkNERELIzwO37sn2FdyTJSIiQH4e9IoLn4iIiDwRQ5aIiEghDFkiIiKFMGSJiIgUwpAlIiJSCEOWiIhIIQxZIiIihTBkiYiIFMKQJSIiUghDloiISCEMWSIiIoUwZImIiBTCkCUiIlIIQ5aIiEghDFkiIiKFMGSJiIgUwpAlIiJSCEOWiIhIIQxZIiIihTBkiYiIFMKQJSIiUghDloiISCEMWSIiIoUwZImIiBTCkCUiIlIIQ5aIiEghDFkiIiKFMGSJiIgUwpAlIiJSCEOWiIhIIQxZIiIihTBkiYiIFMKQJSIiUohXhey6deswZMgQ+Pn54eGHH8bRo0fdXRIREXkwrwnZ7du3Y8GCBcjOzkZ5eTnGjBmDlJQU1NXVubs0IiLyUF4Tsr/97W8xZ84czJo1CyNGjMD69evh7++PTZs2ubs0IiLyUL7uLqAnNDc3o6ysDIsXL5bafHx8kJycjNLS0jb9rVYrrFarNG42mwEAFotF+WKJiKjXas0BIUSH+ntFyF67dg0OhwMGg8Gl3WAw4Ouvv27TPycnB2+99Vab9sjISMVqJCKivqOhoQFBQUHf2c8rQlauxYsXY8GCBdJ4fX09oqKiUFNT06GNSuQOFosFkZGRuHjxIvR6vbvL6TU8cbv01XXq7XV3pD4hBBoaGhAREdGheXpFyIaFhUGtVqO2ttalvba2FuHh4W3663Q66HS6Nu1BQUG98otBdDu9Xs/vaTs8cbv01XXq7XV/V31ydra84sInrVaLBx98EEVFRVKb0+lEUVERkpKS3FgZERF5Mq/YkwWABQsWIDMzE4mJiRg7dizWrl2L69evY9asWe4ujYiIPJTXhOz06dNx9epVLF++HCaTCXFxcdi/f3+bi6Hao9PpkJ2d3e4hZKLegt/T9nnidumr69Tb61aiPpXo6HXIREREJItXnJMlIiJyB4YsERGRQhiyRERECmHIEhERKYQhe4vc1+Dl5+cjJiYGfn5+iI2Nxb59+3qoUvJWcr6jVVVVSE9Px5AhQ6BSqbB27dqeK9QN5GybjRs34tFHH0VISAhCQkKQnJzcK197KWeddu7cicTERAQHB6N///6Ii4vDH/7whx6s9v919pWieXl5UKlUSEtL6xW15ebmQqVSuQx+fn7yFypI5OXlCa1WKzZt2iSqqqrEnDlzRHBwsKitrW23/5EjR4RarRarV68Wp06dEsuWLRMajUZUVlb2cOXkLeR+R48ePSreeOMN8ac//UmEh4eLd999t2cL7kFyt81zzz0n1q1bJ44fPy5Onz4tZs6cKYKCgsSlS5d6uPJvJ3ediouLxc6dO8WpU6fEN998I9auXSvUarXYv39/r6671fnz58U999wjHn30UfHUU0/1ito2b94s9Hq9uHLlijSYTCbZy2XICiHGjh0rjEajNO5wOERERITIyclpt/8zzzwjUlNTXdoefvhh8fLLLytaJ3kvud/R20VFRXl0yHZl2wghhN1uF4GBgWLLli1KlShbV9dJCCHi4+PFsmXLlCjvW3WmbrvdLsaNGyc+/vhjkZmZqVjIyq1t8+bNIigoqMvL9frDxa2vwUtOTpba7vYaPAAoLS116Q8AKSkp39qfqCs68x31Ft2xbW7cuAGbzYYBAwYoVaYsXV0nIQSKiopQXV2Nxx57TMlSXXS27rfffhsDBw7ESy+91Otqa2xsRFRUFCIjI/HUU0+hqqpK9rK9PmTv9ho8k8nU7mdMJpOs/kRd0ZnvqLfojm3zy1/+EhEREW3+4+wunV0ns9mMgIAAaLVapKam4v3338ekSZOULlfSmbr//ve/4/e//z02btzY62p74IEHsGnTJuzatQt//OMf4XQ6MW7cOFy6dEnWsr3msYpERHdauXIl8vLycOjQoc5d1NKLBAYGoqKiAo2NjSgqKsKCBQtw3333YcKECe4urV0NDQ3IyMjAxo0bERYW5u5y2khKSnJ5gcy4ceMwfPhwbNiwAe+8806H5+P1ISv3NXgAEB4eLqs/UVd05jvqLbqybdasWYOVK1fi888/x+jRo5UsU5bOrpOPjw+io6MBAHFxcTh9+jRycnJ6LGTl1n3u3DlcuHABTzzxhNTmdDoBAL6+vqiursb3v/99t9TWHo1Gg/j4eHzzzTeylu31h4s78xq8pKQkl/4AUFhYyNfmkSL4qsZv19lts3r1arzzzjvYv38/EhMTe6LUDuuuf2+n0wmr1apEie2SW3dMTAwqKytRUVEhDU8++SQmTpyIiooKREZGuq229jgcDlRWVmLQoEHyFt7lS6c8QF5entDpdCI3N1ecOnVK/OxnPxPBwcHS5doZGRnizTfflPofOXJE+Pr6ijVr1ojTp0+L7Oxs3sJDipL7HbVareL48ePi+PHjYtCgQeKNN94Qx48fF2fPnnXXKihG7rZZuXKl0Gq14tNPP3W5PaOhocFdq9CG3HVasWKFOHDggDh37pw4deqUWLNmjfD19RUbN27s1XXfScmri+XW9tZbb4mCggJx7tw5UVZWJp599lnh5+cnqqqqZC2XIXvL+++/LwYPHiy0Wq0YO3as+PLLL6Vp48ePF5mZmS79d+zYIe6//36h1WrFyJEjxd69e3u4YvI2cr6j58+fFwDaDOPHj+/5wnuAnG0TFRXV7rbJzs7u+cLvQs46LV26VERHRws/Pz8REhIikpKSRF5enhuqlv+39HZKhqzc2rKysqS+BoNBPP7446K8vFz2MvmqOyIiIoV4/TlZIiIipTBkiYiIFMKQJSIiUghDloiISCEMWSIiIoUwZImIiBTCkCUiIlIIQ5aIJDNnzkRaWpq7yyDyGF7/ggAib6FSqe46PTs7G++99x74fBqi7sOQJfISV65ckX7fvn07li9fjurqaqktICAAAQEB7iiNyGPxcDGRlwgPD5eGoKAgqFQql7aAgIA2h4snTJiAefPmISsrCyEhITAYDNi4cSOuX7+OWbNmITAwENHR0fjrX//qsqyTJ09i6tSpCAgIgMFgQEZGBq5du9bDa0zkfgxZIrqrLVu2ICwsDEePHsW8efPw6quv4ic/+QnGjRuH8vJyTJ48GRkZGbhx4wYAoL6+Hj/84Q8RHx+Pr776Cvv370dtbS2eeeYZN68JUc9jyBLRXY0ZMwbLli3DsGHDsHjxYvj5+SEsLAxz5szBsGHDsHz5cvznP//BP/7xDwDABx98gPj4eKxYsQIxMTGIj4/Hpk2bUFxcjDNnzrh5bYh6Fs/JEtFdjR49WvpdrVYjNDQUsbGxUpvBYAAA1NXVAQBOnDiB4uLids/vnjt3Dvfff7/CFRP1HgxZIrorjUbjMq5SqVzaWq9adjqdAIDGxkY88cQTWLVqVZt5DRo0SMFKiXofhiwRdauEhAR89tlnGDJkCHx9+SeGvBvPyRJRtzIajfjvf/+Ln/70pzh27BjOnTuHgoICzJo1Cw6Hw93lEfUohiwRdauIiAgcOXIEDocDkydPRmxsLLKyshAcHAwfH/7JIe+iEny8CxERkSL430oiIiKFMGSJiIgUwpAlIiJSCEOWiIhIIQxZIiIihTBkiYiIFMKQJSIiUghDloiISCEMWSIiIoUwZImIiBTCkCUiIlIIQ5aIiEgh/wcvPYgdR0KmogAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, (ax1, ax2) = plot_production_comparison(\n", " prod,\n", " pvt_gas,\n", " result.params,\n", " filter_window_size=30,\n", " well_name=f\"Data well {well_number}\",\n", ")" ] } ], "metadata": { "kernelspec": { "display_name": "bluebonnet", "language": "python", "name": "bluebonnet" }, "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.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }