{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Einstein Toolkit" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import CactusTool\n", "\n", "name = 'BBH006Lev1'\n", "sim = CactusTool.load(name, '/Volumes/simulations/BBH_Catalog')\n", "mp = sim.ThornOutput('multipole')\n", "bbh = sim.ThornOutput('TwoPunctures')\n", "strain = mp.Strain(bbh.ADMMass, bbh.CutoffFrequency, (2,2), -1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## NRSur7dq4" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "setting __package__ to gwsurrogate.new so relative imports work\n", "__name__ = gwsurrogate.new.spline_evaluation\n", "__package__= gwsurrogate.new\n", "setting __package__ to gwsurrogate.new so relative imports work\n", "setting __package__ to gwsurrogate.new so relative imports work\n", "Loaded NRSur7dq4 model\n" ] } ], "source": [ "import warnings\n", "warnings.filterwarnings('ignore')\n", "\n", "import gwsurrogate\n", "import numpy as np\n", "\n", "sur = gwsurrogate.LoadSurrogate('NRSur7dq4')\n", "m1 = bbh.m1\n", "m2 = bbh.m2\n", "q = bbh.mass_ratio\n", "chiA = bbh.chi1\n", "chiB = bbh.chi2\n", "dt = 0.1 # timestep size, Units of M\n", "f_low = bbh.Omega_orb / np.pi # initial frequency, f_low=0 returns the full surrogate\n", "t, h, dyn = sur(q, chiA, chiB, dt=dt, f_low=f_low) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Gravitational Wave" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAClHUlEQVR4nOyddZhc1fnHP2dn1l2zFieeECG4BbfiVqBIsbZQKFJokR+l0BYpTmgpFCtOKQUKLe5OSEICcV9395m5vz/eOXPvHdndhMgGzud59pmdO3dmztj5nvc9ryjLsjAYDAaDYagRt60HYDAYDAZDNIxAGQwGg2FIYgTKYDAYDEMSI1AGg8FgGJIYgTIYDAbDkMS7rQcwFMnLy7NGjRq1rYdhMBgMPwi++uqresuy8sOPG4GKwqhRo5g3b962HobBYDD8IFBKrY923Lj4DAaDwTAkMQJlMBgMhiGJESiDwWAwDEnMHtQg6evro7y8nO7u7m09lO89SUlJlJaWEh8fv62HYjAYtiFGoAZJeXk56enpjBo1CqXUth7O9xbLsmhoaKC8vJzRo0dv6+EYDIZtiHHxDZLu7m5yc3ONOG1hlFLk5uYaS9VgMBiB2hiMOG0dzPtsMBjACJTBYDAYhihGoLYjPB4PM2bMCP3dfPPNHHvsscyYMYMddtiBzMzM0G2ffPIJc+bMCSUcr127lnHjxvH6669v41dhMBiGAg0NoBS88MK2HklsTJDEdkRycjILFy6Mett7773HbbfdxiuvvBJxW3l5OYceeii33347hxxyyBYepcFg2B749lu5vOMOOO64bTuWWBgL6ntOVVUVBx98MH/84x856qijtvVwDAbDEKGkRC4nTdq24+gPY0FtApdccklMS2ZTmTFjBnfddVe/53R1dTFjxozQ9auuuoqTTz653/uceeaZ/OEPf+CEE07YDKM0GAzfF0aOdF8ORYxAbUf05+KLxYEHHsgTTzzBWWedRUpKypYZmMFg2O7w+eD442H69G09ktgYgdoEBrJ0hhJXXnkljz/+OCeeeCIvvfQSXq/5yA0GA8yfD//6F5x//rYeSWzMHtQPgLvuuouMjAzOOeccLMva1sMxGAxDgEBALodyTrwRqO0IvQel/377298O6n5KKR577DGqqqq48sort/AoDQbD9oBeq/7f/23bcfSH8fdsR/j9/pi3zZkzhzlz5riOvffee6H/ExISeOONN7bQyAwGw/aGtqCGslPFWFAGg8HwA0QL1FDGCJTBYDD8ABk3Ti6HctMAI1AGg8HwA6S0FIYNg6KibT2S2BiBMhgMhh8gbW1w+OFwwAHbeiSxMQJlMBgMP0A+/hgeeQSGD9/WI4mNESiDwWD4AaKDJBobt+04+sMI1HaEUorLL788dP22227j+uuvB+D666+npKSEGTNmMHnyZJ5++unQeZ999hm77rorM2bMYNKkSaH7DJa2tjZX/lVeXh6XXHJJ1HPT0tIG/bgXX3zxRp1vMBg2Hzq83DGlDDlMHtR2RGJiIi+88AJXXXUVeXl5Ebdfeuml/PrXv2blypXstNNOnHDCCcTHx3PmmWfy3HPPMX36dPx+P8uXLx/0c1qWRWpqqqsG4E477cRx37E+/7x582hqavpOj2EwGDYdkwdl2Kx4vV7OP/987rzzzn7PGzduHCkpKSEBqK2tpSgYquPxeJg8eTIgVtdtt90Wut/UqVNZt24d69atY8KECZxxxhlMnTqVsrKy0DkrVqygtraWvffeG5BGiLvvvjvTpk3j2muvDZ1nWRa//OUvmTBhAgceeCCHH344zz//PCAJx1dccQW33nrrZnhXDAbDprA95EEZC2oTCSvaAMBJJ8EFF0Bnp0THhHPWWfJXXw/h3S8cRR/65cILL2THHXfst2TR/PnzGTduHAUFBYBYVhMmTGDOnDkceuihnHnmmSQlJfX7PCtXruSxxx5jt912cx1/5plnOPnkk1FKAfCrX/2KX/ziF5xxxhncd999ofP+/e9/s3z5cpYsWUJNTQ2TJ0/m7LPPBmDu3LkcddRRIdE0GAxbn5kz5bK0dNuOoz+2awtKKXWoUmq5UmqVUiqiMJ1S6jKl1BKl1CKl1NtKqSHc+WRwZGRkcMYZZ3DPPfdE3HbnnXcyZcoUdt11V6655prQ8euuu4558+Zx8MEH89RTT3HooYcO+DwjR46MECcQgTrllFNC1z/++OPQ9dNPPz10/IMPPuCUU07B4/FQXFzM/vvvD0BlZSX//Oc/ueiiiwb/og0Gw2ZnxAiYMgWysrb1SGKz3VpQSikPcB9wEFAOfKmUetmyrCWO0xYAsy3L6lRK/QK4Fei/w98g6c/iSUnp//a8vMFbTNG45JJLmDVrFj/96U9dx/Ue1Msvv8w555zD6tWrQ5bS2LFj+cUvfsF5551Hfn4+DQ0NeL1eAg47v9tR1jg1NTXieb/++mt8Ph877bST67i2pgbDggULWLVqFTvssAMAnZ2d7LDDDqxatWrQj2EwGL479fWw994QXDsOSbZnC2oXYJVlWWssy+oFngGOdp5gWda7lmV1Bq9+BgxhY3bw5OTkcNJJJ/HQQw9Fvf2oo45i9uzZPPbYYwC8+uqroTYbK1euxOPxkJWVxahRo5g/fz4gbsG1a9f2+7xPP/20y3oC2HPPPXnmmWcAePLJJ0PH99lnH5599ln8fj9VVVW8++67ABxxxBFUV1eH9rpSUlKMOBkMWwHLgt/+FubNk+vvvAP33y9W1FBlexaoEqDMcb08eCwW5wD/i3WjUup8pdQ8pdS8urq6zTTELcfll19OfX19zNuvu+467rjjDgKBAI8//jgTJkxgxowZnH766Tz55JN4PB6OP/54GhsbmTJlCnPnzmX8+PH9Pudzzz0XIVB333039913H9OmTaOioiJ0/Nhjj2XcuHFMnjyZM844g9133/27vWCDwfCdaGqCW26BI4+U69p5UlYW+z7bGrW9NrBTSp0AHGpZ1rnB66cDu1qW9cso5/4E+CWwr2VZPQM99uzZs615epkRZOnSpUyaNGmzjP2HyFlnncWPfvQjTgiPDomBeb8Nhs1LebldNcKy4Kmn4LTTpFjsmjXbdmxKqa8sy5odfny73YMCKgBnkY7S4DEXSqkDgWsYpDgZDAbD95GuLrnUAcDaNhnKNsr2LFBfAuOUUqMRYfoxcKrzBKXUTOBviKVVu/WHaNA8+uij23oIBsMPmsJC+Ne/YHbQTtEuvo2IcdrqbLd7UJZl+RC33evAUuA5y7K+VUrdoJQ6Knjan4E04J9KqYVKqZe/43N+pzEbBod5nw2GzU96OuTmwvr1cn2ffeQymC45JNmeLSgsy/ov8N+wY9c5/j9wcz1XUlISDQ0N5ObmblRYtWHjsCyLhoaGAROJDQbDxlFfLwUGZs2Cr76CkSNhr70gPn5bjyw227VAbU1KS0spLy9ne4jw295JSkqidCintxsM2yHBjBK+/louy8thxgzYb79tNqQBMQI1SOLj4xk9lHsjGwwGQz/oIAntAHrtNZg7F37zG7luWUNvP2q73YMyGAwGw+DRhWLigrO+DpJYtgwuvVSOD7XtX2NBGQwGww+Avj651CKkBeoXvwBdzKW3FxITt/7YYmEsKIPBYPgB4PfLpe5l6rSWdHcGn2+rDmlAjEAZDAbDD4CDDoI33rBbATn7QenyR0agDAaDwbDVKS6WTguvvCLXtShlZMAee8Ddd8NQy+4wAmUwGAw/AFavhkMPlaaqIP2gDjkEvF647z64666htf8ERqAMBoPhB8Hrr0N7u1Q1B1i5EsaPh9/9DqqqYO1aO9JvqGAEymAwGH4A6CAJzX/+A/feK9Uk3n5bjlVVbf1x9YcRKIPBYPgBoAUqPMz8yy+d5wQYShiBMhgMhu8ha9e6o/LCLSgtUBdeaB+rqqrZ8gPbCIxAGQwGw3ZKVxd8/HFkBYjGRhgzBq64wj6mBerTT+XSfZ9nAKitjd2le1tgBMpgMBi2U+bOlT2ksAbgocoQznbuP/mJiNNOO8n1gMub9xIALS2tW2ysm4IRKIPBYNhOefNNudywwX28vV0uL7rIPlZcLJf33y+XZ5whTQy9XgtYCVyIx2NcfAaDwWDYDOj6eqmp7uMdHZHH58+H44+Hiy+W6yUlkqDr9/uBa4Cf0dNjXHwGg8Fg2Az8/OcSKn7ooe7j3mAZ8FNPtY+9+CJUVoprz7JgwQKxqq6/vg4oBnakqSkskmIbY6qZGwwGw3bKySdHP37YYVJ7r63NPubccwoERLDmzoWf/7wR2BWAioqh1V7XWFAGg8GwnVJdDRMnwk03Rd6WlAQ9PfZ1Z5i5z2cL1nvv2XZKV5fjDkMAI1AGg8GwnXLQQbB8uTtaD+Dxx6VShNOCcgqU328L1C23DA8d7+7u3YKj3XiMQBkMBsN2Smen+1KzfLlchosSwJo1Yl1pgQoELOAvAHQPsWJ8RqAMBoNhO0VH6+lLjdaZE06wj11xBXz7LYwe7W7v7vdbwKvB+/Vt2QFvJEagDAaDYTsllgXV1yd9nm691T5WWCjC9ac/yfmXXAIzZmgLqgGlfk58/JqtNPLBYQTKYDAYtkMsyxammTPdt/X1QXy8u5zRO+/AOefANdfI3lRhIUybpl1/N+LxnENfX8PWGv6gMAJlMBgM2yGWBbfdBh99BH/4g/u24mJoaHA3IHz2WVi4UP73++HDDyE/H84881WgCJ9vZ9rbk7fW8AeFESiDwWDYDomLEzfdnntG3nbttXD99WJJ6WCI8DDzF16Av/8d8vLWAlMBaG4etqWHvVEYgTIYDIbtkL4+WLpUknX33z/ydl1NQrfccCbq6jDz1laYP3+04zFNJQmDwWAwBKmrg9xcsYg2hqoqmDxZ/p8wwX3bVVfBzTfL/319kJAQOw/qtdfsOklDTaCMBWUwGAzbCMuCkSPhyiv7P6e52d18EKC3N/r/AMuW2f/r+/n9UiC2rk5CzbVAWRakpNwRPNd01DUYDAYDUF8vTQdvvx0qKqKfExcH2dmwYoX7uLOSebhA6dt+9Svb1Xf33fDFF5CXBx6Pu/V7evqnwfuFqeA2xgiUwWAwbCPKy+3/162LvN0pPA0N0W9LS3PX3AMRqF13hbvusltu5ObKY1x9NdTUwI03SlFZy4LExER23/0BPJ753/UlbVbMHpTBYDBsI5qb7f/r6iJvb2mJfi7YVtIee0BOTuRtcXFSYSI5Wf7/5z8lau+NN6QNx9Sp4l4MBKCh4Ur8/jiUunFzvKzNhrGgDAaDYRvR6uiwHm4hgd0ZF9xiBSIuDzwgVtLf/+6+bepU+OYbsa7Wr5djTz4p4gQiSq++KtUmZs78NT5fHhUVU+nqyv/Or2lzsl0LlFLqUKXUcqXUKqXUb6Pcvo9Sar5SyqeUOiHaYxgMBsO2YsIEuPxy+d9ZeVzjPBYuUPn5cN55MGJE5P3uuQf++lf5X1tazii+QEDyoJ56CjyepfT0SD/47u4dNvGVbBm2W4FSSnmA+4DDgMnAKUqpyWGnbQDOAp7auqMzGAzfZ+rq4Gc/c+8hbQoTJ0o4eEKCLSROkpNln+iww2D2bPdtLS3w5ZeSkJueHhnlFx/sPeiM4tNYlohUeTmUlx8SOu7zDa0w8+15D2oXYJVlWWsAlFLPAEcDS/QJlmWtC942tGInDUOKZ599lo8//pg77rgDr3d7/kkYthZlZfDEE2LFhJcZcvLccxI5d9tt0W+vr5d6et3doFTk7ePGwX//G/2+n3wChx8Oxx4rrsDeXjti76ij7LJGsSwoHWZeXv6r0PG+vqE1VW63FhRQAjjbdJUHj20SSqnzlVLzlFLz6qLtVhq+l/j9fk499VTuvfde/htrJjAYwpg1S1xrS5b0f57XC3feGd19ByJc48dDW1srP/3pT/nyyy+jnlddHblH5QwzB3fE38qVUF8vceROC2rXXUUMZ8505kHFMXv2E4BUNg8Eho5Ibc8CtVmxLOsBy7JmW5Y1Oz9/aG0UGrYcCxcuDP0g33nnnW08GsP2QmurCNSGDf2fpy2VtWuj397SIoEKJ5+8gUcfTee3v3VvpT/9tFQdLyqS+nqBQIAbb7yRr7/+OiRQaWly6RSovj4oKmoAbmThwtcBePFFCZJITJSoPqcOjRypB6joi+Zr3EZszwJVAQx3XC8NHjNsQXrDMwK3c5YEl8AZGRl8880323g0343W1lZ6whNihjjl5eXcddddW33cneENlDaSnXeWyb4/gVq1Cv78Z/lfR9KF09oqAvXllznAfqwLS4ZqapKcpcRESej94osvuO666zj66KNDghTNgpIisd8A1/Hll/8GRMgqK+Gii6Sr7t/+BmedJftRPl8up532KvDmkPqNb88C9SUwTik1WimVAPwYeHkbj+l7zc0330xaWhpvv/32th7KZmPFihV4PB6OOOIIli5duq2Hs8m0tLQwceJEDjvssG09lI3i4osv5tJLL+Wxxx7bYs9x883iItO8++67ZGVl8fTTT7vOsyzJKXrmGff9e3t72RCmRDU1cun1unsuOVm5UvafwP38TtrbRTj8/iYgnfLycpeLTYeZ5+fLXtUXwQdcv359yIKaPh3OP1/auNtjhvb2NqCAZcvEOrrvPvj972HuXBGqtDRJ3oU4Pv74OL75ZgLQYQRqc2BZlg/4JfA6sBR4zrKsb5VSNyiljgJQSu2slCoHTgT+ppT6dtuNePvG5/Nx++2309fXxz333LOth7PZWLFiBaNHj2bSpElUVlZu8ZX81VfDBRe43SubgzfffJOqqireffddysrKXLdVV8P990eWw9kS+HySEOrM34mFZVm8++67AJu06Bns63nwQbjhBvv6ww8/TF9fX4Qo1tTAp59GJsxeeOGFjB07NrSA6e4W19wf/ygTfbTgBnDvGTU1RT+nu1si9Xy+ViCN3t5eqh1qpt/H3FyxoNY7TLHZs7t46ik47jixhvLy7Mfdf39oa+sDalixQlpoPP44PP+83B4ISIh5YqIFTKOrK42vv94BmGAEanNhWdZ/Lcsab1nWWMuy/hg8dp1lWS8H///SsqxSy7JSLcvKtSxryrYd8fbLggULqK+vJzMzk/fffx8r1rKxH/72N0hJCfDeext/381JUxP85S+yabxmzRrGjBlDSYnE11RWVm6x562pgZtukvyUDz+Mfo5lWbQ6szcHyaeffhr6//PPP3fddsUV8ItfiEhtSbq7u/nnPz/goossfv3rgc+vra2lOVgeYaEOORskxx0HJ500mDGJO6u9XRJXARYvXgzY7l3NqlVyefHFduKrZVk89thj+Hw+Xn31VcC2noYN0DqpsVEu33hDOtlG45e/hF//Gvr6WvB4xFdXW1vrGn9ysrjxurpwWXKJiZWccorcZlluS+7vf++ip+cWAOrqmgkEAvj9duh5IAD//jf8+98WsJaOjuBGFjOMQH2f+fhjmL8ZylndcAMcc4w703xTuO8+WVkFf1ubzNdffw3AeeedR0tLS8QqfSAsC664op2urv/w5pvXbvI4mprggw9iu1UG4qabxIK58EJYvBiqq6spLi6mtLQUkD2RgdhU62fYMHvSeuut6Of86le/Ijc3d6MtipUrVzJ27FgAli9f7rrtoYekOOgrr2z0kKPm5sTi2muv5dRT9yUvb32oYkF/LAuW3J45cybr1q2js3NwOThnnimT66JFA5+rt3Qef1yEwOfzhZ63rKyMjo6O0LmrV9v304VZq6qqQkEDCxYsAOzf5IYN8huNFcnX0CDW1f77a1daJEceCccc46enp5y0NHn99fX1odtnzoTTTpOir7/4BVQ4KsrOn9/AO+9Iwm1cnB1WLs/dAEj4ns9nUV9fHyFQgQAsXRoHXO4YUZwRqO8zf/ub/IA2dQIFmTh/9zt46SVZ6W8qvb0S+dPQAI88sumPEwjAn/60K3FxX7P//scBsGgws4ODJUss2trSgJe49957N3mT+rDD4Jpr7FXsxuDziQ9eT7offhigpqaGoqKikAVVEaukdJD6ehg+HG65ZeOfH6Qq9dSpEJzrXDQ2NvLXv/4Vn8/HXXfdFXH7zTfD//1f9O/WunXrmDJlCqWlpaEJWJOQAD/9afTn7I9bboHMTAiuTQbk2WefBRTd3e+wdu3Abr7VQUVISvozvb11HHTQwGoYCIgLEUR8urv7P3/NGrkcNkzccdqNu88++wBui9lpGen1lzNoYU3wwfLyJDl29Gj5jcYKgOjqEmF64QUR1GgsWQILFzYC53HaaQ8CboE65RRxUZ50Ehx/vHxHJgebQL38ciIHHGB/n52LiV12KQAuDV6Lp6GhAb9fgi3i4+1EXeFPjhEpI1DfZ/bcU1wJA+VH9McLL8jK65BDICtr0x/ngw+kwOQzz9i+503h3Xdh7dppBAI78vnn0hpa/1gHy7//LW6Lffctoa3tTh58cOO3A1etgs8/lw6ihYUbfXfWrZOqz0ccIRPvggW99PX1UVhYGLKgBhKovDz40Y9k8vZvZNL9L38pUV077uju16P58MMP8fl87Ljjjrz//vv4HU+wdq00ofvDHyKFxrIs1q1bx6hRoxg1apTLuv30U7j0UigoEHF1zH39YllS462rCwYTv1BTU0N5eT3Qxbp10j0vvD1EOLLXovjyyzlABosXewZ8nooKGdMuu8gYB/LIaoN4552lwZ8WpNnBsgxOgcrJgd12E+HR75P+PsyePTskVkVFsoDcfXc5J9r+0htvvMGGDT9m+fJm7rxTPBnROPZYuPZaMWumTpXfVrQ8zKoq+c40NjYybtw4ANraRJ1TUuQcne9kWVBVlQAkBu/tpaGhgUAA9t5bFq4HHOAUqDjOP/+D4P+PcNhhQ6fckRGozczBB8tlcO93k/j0U5gxA157DX7+801/HJ3zd8gh/Z83EBI41E16ejXvv59GcnKya7N2MHz0UQPQyVVXHQmcw8svtwx0lwi01+uQQyTxcWNdbVoUJkyASZPg22/lF11YWEhGRgbJyclUVVUN+Dh77SWT0sYG/T33nEzaf/0rhHnhAPjss8+Ij4/nwgsvpK2tzeWqe/11+7w333Tfr6mpiba2NkaNGkVJSYlLZN97T4Rm//3h3HNjWxx+v5/nnnsuNAlXVEhwxdixUq1gIMSiHgkk4vd/BNh7OrGorq4mI2MyPp8HWE1bW3xEvblwtJ7ssotcDvRxnXgizJsnv6f6eigrkwCEnXbaKfh4tkCde6789nRbCrAFas8996S6upru7m46OuT90cIQnoRrWRYXXnghzz77LI888hBpabGtye5usKwu4Eyefvp4lFIuC+r00yVK76qr4JBDLJqamhg9Wlq0t7eLpZOcLOdqC8oueVTGCSd8CSymvr6ehQvh2Wft57Z/P4qpU/WLSKCsLJGhghGozcyoUbLCipEQPiheekmS6kC+dF1dm/Y4e+0lLq2MDNhzz2rOP3/j3HKaK67wodQIDjrofc46SzFixIiNFqi8vP+RkHABBx64E15vFStWbPxXb/FieS1ffCGX0frn9Iee7ydMkBpoq1bJyrWwsBClFHl5eUHffXSqq+Ggg+zosS+/lMno7LPPZty4cazqZ0ZubpbosPHjZeyeKMbCV199xbRp00KTp9NVd+SR8p2YPj3y+6An0ZKSEkpKSqisrAwFsaxdK9bTAQeIqyhoKEZw9913c/LJJ3PAAQfg8/lC79WDD8KBB8Z8WSFkrKOD117n/vsXDLgwqq6uJitrVvCaJEkP9LXS8QN77w1nndXNkiUf9xuwk50NO+1kv+4VK5oBW6CiWcz77CPh5vr2xMTEkHVTW1vLyy/L42mXYLhAVVVVBb8L1/Lgg7mkpUnbi2hIiaNuYDyffJJPdna26zvY0SFCkpgIPT0WPp+P4uJiEhIS6OzsIy5ObgNbmGxXXw0XX9wNLKO+vh6PR1yXZ50l+1X//jeceaZYsZWV+Vx77XzgE3bc8TtufG9GjEBtZpQSF853yflMTJQs9aoqKQL5+OOb9jh77w3XXScbqJ980sKDDw6wpI1BRUUFllXH4Ye3c+aZMHLkyI0WqNbW95kwYT4ej4fs7Ebq6jI2ehxLlsDkyeKC0dcfeeQRLrrookFFvlVWiosuN1dChH//e4kaKCoqAiAvL8+1eg1n8WIJbhgxQvz4y5ZJxNwjjzzCqlWruPHG2L10tKtp5Ej5/xe/iAymWbFiBRMnTgy5cFY4fGQlJXD00TKxXH+9+37aJVRQUEBxcTGdnZ20BE2RykoolkLVBALRFzuWZfHgg7L/sWbNGt5//33231++fzNniktxoC3DsrIyPB7tGlpOevpSMjP7v091dTXJybKfEhcnLqaBvlYeD0yZIu61pUvncP75e/GPf/wj5vmvvCJWw6GHisC3tGwgLi6OHXbYgYSEBJc77fjjJYDmyislVwjku19SUkJBQQEgAqXFpqBAFqTagtHogKKEhGNZv370ICyobqALv1+Rnp7t+i739Mh8kJAA3d0ixDk5OWRnZ9PZ6SMhQcTysstsEbYFKoDXOxHIoKGhgeuugzvuEJdtWZnkcMXF9QFxPPLIND77rBCw8HpNqaPvNXfeKftIm8L8+fJlq6qSzVqlBvblR8Pnk8nMXrktAcZv9ONs2ABnnJECzKKkZARLlkBu7tRBRbtpenth4cIRlJTIarm0tJeenmJXBNVgeOopePhhcc8BvPFGBWeffTZz587l//7v/wa8/+232y6h4mLo7pYExmHBeOGBBErvK+64owSf7LMP/Oc//8Hj8XDiiSfy4osvxiwTo5+3qEg+m/vvd0dddXd3s2HDBsaNG0dGRgaFhYWsXLkydPvDD8fe19STbH5+fkSwR1WVLVB5eZKHFU55eTnLli3jpptuIj4+ntdffx2lZJ/v44+l7txAgRJlZWVkZIwJXqvhtddSeeml/u8jAiWLg9zcVUyZ8nrU1hFODj9cFn+9vav5/PMvgSQeffTRmOf/9a9w662yMDj6aGhoWE9hYSEej4ecnBwadVglsgAJ309qamoiNzc3qkClpYmFeuGF7vvoMPaMjEK6uxtITvYNwsXXFXzM/KgClZhoW+7Z2dlkZWVRVPQ/XnhBFmy33y6eARARnzx5EUo1scce+Xi9Z1BfX8/jj4vLF2SxMncu+P19wDA6Ory89VYxsCfz52fFfD+3NkagtgCTJskPYlP47DMROMsSy2fMGHf462DZsEFWv889p4+sBUYTCGxceOHXX8MHH+QDCRQUjGLKFKiuPoi6urpBF5XcsCFAWdlcLEuqHEyd6gVa+eabKJEC/VBYKO9tdrZMtm+/vYbk5GSOPPJIHn300UHVENPVnquq4OWXp+P1TiY9PR0YWKDWrpWck4ICsUyPOELq9+2yyy6cdNJJtLa2Mm/evKj3DQRkAikttQM8nNUF1qxZg2VZIetphx12CEW5dXRIHs2LL0pl6513dieTaoHKy8ujMPjgOpemp0esL5BxR1tX6DHPmTOHGTNm8OWXX/LAAzKB6VX5QMEIZWVllJbWcvHFFqmpSbz11o4MlM9dXV3NPvvMp6cHiou7GTXqXqZP7/8+GknwrWDEiOf57LPPYiZYNzTId6WzU967tWt9oQVJTk4OTQ5Fqq+Xc594QgImqquhra2N9PT0qAKlSwyFs3btWnJzc/F40oFOTj99VcyQ+EcegUmTFgOyOZiWlhfTgurtVaFxZ2VlEQgs4bDD5LvV3o6rNt+++95PdvbHwXFmhKL4nGHmL78MX3yRDdTS0TE0q/gbgdoCNDSIyGxK5Zy1a6VkiZ7Edthh4M3maOgguzFjQkeAVMrLNy6E1A7WW8XkycPJz4fu7lL8fr9r9dkfixaJT33cONlVvvzyADCJsrIYFTSjUFEhOUx632nECIvVq3s54ogjOOuss2htbY1IUA3n7LPtTeKWFvjww4NISdkXFSwFMJBAlZdLiLlSsvJdscLPggUL2H333dltt90A2UeKxiGHiEtw9Gj5fLOy3AKl96+0QDmDHbQ4DB8u1te8ee79Nz3m3NxccoMJN3of45tvJPVBHlPex3DmzZuHx+Nh+vTp7LLLLsybN4/HH7d4/nnb+hqMQM2YUcHddyuKi4vxeGr7vU9fXx9tbW3k5eWRkACFhcOoqmqJWXFB8+tfw49/LEEZSrVSVDSV7u7uiNwvjRad+npZUGzYMJGcYH90pwVlWbJPmJ0tFkhTk1xvbW0lIyMjJFB1dXV0dMg5CQkSWHHdde7nLC8vp7S0FJ8vEeiktXUNQS9yBKeeChkZq4E6JkywSEvLpM2xqXXMMRLpd/zx8POfi09YW1Dl5QW88YZY1unpuCzWpqYmsrIk+TYpKZWWlpaoeVDLlmUDN2JZuhxGVej9GAoYgdoCdHSImy5WtYD+WLNGJrG44CczdqxYUBv7hXEKlIQrLwE+pqZm4/KP1q8Hj6eXggIPSUlJjB0Lra1S7d2Z8d4fixbJJDB1ajYAo0aNAmSlOVgWLxb3lJ5gTz65np6eh9l///3Zb7/9AAnTjkVvr6xWddzB8GCZ4cTEsaFz8vLyaG5ujmmJFRTY0WP33AMTJnjo6YlnypQpoX2K+YPM0i4sdEegaYHaYQfZx9ECZVlWaKIvLrYtc+deTV1dHTk5OXi9XvKC9W6cQqtL8eTnRw8zX7x4MRMnTiQ5OZlp06bR3t5OWZmPoiLZr/N6+xcov99PRUUFBQVj8flk7D5fNf2tX7SVsHjxnvz+9+JmXbTocX75y9j3AXGLbtggY05N7QZEbGLVUayvl9eg0zXa2rxkBa84BaqjQ35jGRky2cu5tgWVmppKUlIStbW1HH64uNSUkvGEf+RlZWUMHz6czEwFNPDpp31cf33k/l9vL7z/PlRXK9LT/8eyZYr8fL/Lgrr8cinuutNOMHOm+FmzsrLIzMykouJYzj3X9gror+369fD88w/R1yc5i0lJabS2tuL3y+IoJ0fuYztAnInzoqRDpeOGEagtQEmJrFQ2Yv4NsWaN0+qRxmNXX73xddTWrJEVXkpKM9deey3wLrAX+fkbF969fj0kJ9dQUiJL6TFjoL5eAhxqBpktu2aNuF+mTs0CoKcnE6/3Xd55J62fe7nRrintciopeQN4ij322IPs7GzGjBkTyvSPhrZWtLsrNRW83nY8HtsXq62PWJbhX/5i5wTpx4FiJk2ahFKKWbNmxRzDr34lybIabQ1pKioqSElJITs7O/j4JXR1ddHc3OwSKL0Sd1pfdXV16BYx+jXU19ezdCmccIKIO8QWqNWrV4eEcUJwI6O6Oo6iIlkoFRX1L1C1tbX4/X5efPE8DjpIxKa3t4rGxtgLKz0Jr1gxjpdeksVBIFDNQGuehgYRnG+++YbsbOjpSUUpFVWgfD6xlHNzRXQ8HujoSAi9x06B8vkkUnLCBFugWltlnOnp6aEoz/r6evbcUz5PkHPDo/i0QC1ZEgCuZtGieH7/+8j9rfp6mDMHlizZgbRgz4yMjIwIF18gIBGDixenAZ6QYPb2KhISIjvn9vWBz5dCYqLckJSUQktLCwkJEgXa0CD7cc7P5uGH3V+Mjc3x21IMTcfjdo7HI5FemyJQ7e1ugZozR/42lrVrZbV91123c/PNNzsefxBVPB0kJ0N8/NKQ337kSGhoSALUoAWqrCwAtDNunDxGair4fHNYs2bwb1BZmaxYtcvpk0++Ijl5BuPGSfjvzJkz+7VetOWl7w/g8dSjVEHoutP6GDZAoTX7cUqYFIzamDRpEu+//z6BQIC4OPfa74sv3HsWr7/uLjJaVVVFUVFRyN3oDHaorJQJtahIJsS4OHclDadAJSQkkJGRQX19PWvWwL/+JbX4QAIM8vNlYtLPHQgEWLNmDYceeiigBSqdnh5PSAzvuYeYLir9fgH09qaSnS0i2d1dic8nk3dGlIBNHWXY25tMYaG4rQKBampqAvS3bm5pgYkT+6ivr2fmzDiamuIYNWqUK+JR4/HI556cLK83KwuamxNdAqX3oLKyZE8G7ERobUFlBF9AdnY2TU1NrF8v7+GoUbLf43SbdnZ20tjYSGlpKYmJieTk5NDRIaobHhOkc9L6+tqIj9+DffaBMWPGuQRqwgTYd1/Zd7z77pOBX5KamkpaWho+n4hTuAWlhSolxeK+++Cll1axbl1LRGsQZ/PoWbPceQ+xCuBubYwFtYUYNWrTBGrlStm/0gQC8gMI1tQcNBdfLCGlbkH6kNtuy9qox3niCUhLOzc0YZ96KjzySBsQN2gX3+TJbxIXdxR5ebK6T02F+PgO6uoGvz4qLxe3mF4tvv9+Jl1dC1i9Wn5YM2bMYPXq1TEFWFsAtuUDUEcgYJeA1tZHU5SNkMpKCW/W9ez042RkTAq5jMaPH09XV1fUgrPV1e7qF+ETgBYojVOgzjtPXEmZmTLpHnKIu1BpfX19SFz166ivrw+52IJbLhx2mOyXhAtjd3d3qI5fYWEhqamjiYvrC433mGOkE2ss9H5XZ2dSSKC6uuZSVhYINdMLRwtUV1cSWVkE94UaqKvr35fd3AxKickyZ04TF1wAI0aMiGiHAfaCJqhHZGUF8PszQp9XdnY27e3tEaV9CgvhjDMgN1cqjeggmuzsbJqbm7n4YnlPINKC0tGt+fmjOfJISEv7MW1t4ssN/2rquI7e3laSknL58EOIiyugs7MTX1BlnEESAB5PCgkJCUGBiotqQdkCFccFF8DIkY2h97u6WqzqDz6QhO/DDnsfgPfeS+HBB8uAFey884bQY25rjEBtIUaMiB4xNRicSZxVVeLWCu9RMxB77ikledwr+VzWrNm4b55lWdTU1IQEaupUOPXUNDyewbv4uruXUVi43DWW9PR2WlpSB10VvbzcnWRaWSnLXD2E8eMlhH51jJDHnh6xApyWQGLicZx00oOh63riao6yGigrc4d568fJyJgUOqYDHJzh4RodTab573/Ffavzi6qqqkIRePL48gRSbUFcM1pY/vtfXHs1zc3NIasACCUca53VAuXzyQTlrCah3y8tUEopdthBcdhhx3HKKfocu4pHNLQF1dERHxIoy2oiJaWZuBgzjLYSOjsTyMoiOP4mmpv7X7rvthvk5sqEf/zxcVxxhQhUtOLF1dWSqK7jJ+6+uwn4fei9ygwmarW2tvLBByJMn38un+1jj8HEiTKpawsqKyuLpqYmOjpsa3jiRAlk0ujfRHp6Ma+8Aikp42ltlQVLuCtQfw69va2kpsqPPiEhI3iunOwMMwdISclCKUVaWhqW5cXrDZCRIQuPWcGcZy1QycnxLFwIHs9wWlpaOO88KRz8r3/Ze5g+n5hdv/1tAh98kAH4UWqIbEAxSIFSSu2ulLpZKfWGUmpe8PIWpdSeW3qA2yu33y7W0MawcKGUZnF6K4YNE5fOACXiInjtNfkSdrtq21RulNWyejXssouf3t5dQpNnZye88UYc2dnTo9YMi8a8eZPIzJzjOpaX10sgMGzQVtgrr8D//if/y+QraqH3YrQ4xKrmcOqpYgUFg7Hw+Xy0tpaTm5sVOqc/gdJCqIUpPR0KCv7E+PH2yj1agi2I66WtzV3RurIS/vMfu6ROdXW1y4LSLru6ujoefth2P0WjpaUlNNmCHY2oLSgdIPDhhzJ+R2eOUE1FLVCgJ/wNIRfQ/ffLYicWIlCJdHfHhQQKhvH731sxI1n1ij4+XpGbqy2oVzj33A39BgT9978wbZqETxcWjqSyEkpKRlJRUeGqXQjy/b/+ejtNY8yYWmBF6HPWwtPa2hrqXOu0HFpbRSTCLSinQF13nf29BHv/MjlZVgWZmfG0tcmXJzxIQv80e3paSE8XgfJ600NjktvcFlRysoxd9qyu5Lbb2khLEyEOlhckOxsSEx8nP7+DmTNh2bI5dHd38+STVij3LhCAG2+E3l4/Xm8moHj88UxgEl98MWrAaMqtRb8CpZTaXyk1H3gCyANeA+4PXuYC/1BKzVdK7b/FR7qdkZ1t1+oaLMuXS1FXZxCZ1ysru40RqO5ucec88US4u6qS+vrB19mqqIB587xAQsiCamyUx/Z4fhTVFRaOZcHixb+kp+c41/EpU7qA9QMWZ9XEx9sTvPTxkR+9Fg49wfZXbsiJiNAc3nzzuNCKsz+B0lqsrSDLsujo+BPTp9uRDqWlpSQlJUUIVFeXlBrSiZRgv5b6eujq6qKlpcUlUOnp6aFKB7ffDs5iCdddZ0cT+v1+2traXAKlXXyJibLC1xZ5UPNcOVT6/S91mKdKHciyZVeEJqncXPlOxaomIQKl+N3vfMyZowUqi3vuySVWmyc9AX/1VSN//rO2oN7nkEO+GXD/Y/369Xg8Hj74oJiSEkhLm4Tf74+oo6gtFh308OGHPuCYkAWlBaqtrc11rmXJ5Z/+lBo8Jg/gtKBi/ba1uzMxMSd4nwQ6Ot6juTmyZNSECRIa7vfPIz1dFMjrTXO9P+EWVHKyfM4iUEsYO7YlWBzWfr0S6Xkeo0Y14/VCQoK02g0PM3/lFVi3biTJyRbhaWRDJUhiIAvqGuCyYEPAcy3LusOyrL8HL8+1LGsscBlw1ZYf6vbF6tUSIroxSbbaGgiv1F1aunECpY0S6T/kjEirpLk5edAh67ZxUxsSqMJCcTXFx48alEB1doJlJZGb637S3/62GThtUALl84lLS2fBi0C1kpgYCAmUzlWJ5l4DqT59wQX2dXlfJvHxx9NCkW16ku9PoPQkLxNVHl7vjNA5cXFxjB492tWiQcYmJZKcDfa0QDU0EJpYnQKllCI/P5/6+noaGmw3HYhYfP21TKR6IotmQV11lTsXT4urM5KvpqaGjIwMkhz9wv3+afT2/oSOjjbX/WKVKayvryc9PZ7rr/ey555aoBr7vY+2oPS4xYJKYskSX8zak8uXy97uZ59lUlpaSlZWXPAxJGcg3M0XLlDPPZcN3BlaiGjhaW1tdZ0r329obpZZ2unia2tro6PDCllQjz0mhWj1Npb+vSUkyHNkZyfS0lJPWpo/ov5iTo64ebu61pOV5WGnnSAnR8zWjo4OLEvKLu27r+wBzpp1E5mZ8r6lpqYCh/LWWwG6u2WvTbfm6evz0dPTQ2pqKh4PeL1Si8kpUH6/fH8qKsbg893umBNWh24fCvQrUJZlHWBZ1nsDnPOeZVkHbdZRfQ9oapIghW83oquEdjE4thMA2ZDfmP2s2AL1JVOmrBywh0744zgFyuvVJZhKBpWoW13tD47F/escbP8lkEnuvvvs93LlypUkJSUxd64kMGrGjh0bsw3Ixx+7y/WIuMoL1K8zPj6e1NTUqAJVUiLBCXrlLJPhvfzzn0e6zou1YR9ONIEqDFuZ5OXlUVtbR2Oj2z2Yny8TYnu7PdFnOfqy5ObmRt38dz6nxrm/qImPzwf8NDeXBcchx2PlMDc0NJCTU0p5uVj/IjZNEc/lpKWlBa93JD/5SSKffKItqP34zW+OiVl1oaFB3Hb19fWUlpaGogPT0iSkciCBiovrBNIjLCinQOnHTE+HlpZA8H/bxQdwzTXtnH++nNfYKN8rbV02NDTg9XpJTU1hhx1g2LAELCuDX/2q1+VaBVl0vvoqtLYGKCwMMG8e7LefPFBHRwdKwZ/+JB0SSkogI+MNMjOt4GtOA37P3/+eFRHF9/rr3YBFdfVEPB6Ijxfzy+9XJCWJhZWaKlaUz5dAV9d5jlGNDZ4b/TPY2pggiS3EYDPwnVRXyx5J+Mbyz38ujeoGi7YqCgrCXXz/4phjnokobhmL2lpQygIaXJNYcTEEAsMGJVCrV7cE7+MOzpg/vxBYyrffDhz2rgVEWy/r1q1j5MiRnHtunCu6bPjw4TFrBOqETY2M3S1QIBN9NIE6/XTZ19OICNXR0eGudxNNoF57TTbSndZMfr5YA0rpvkhuC0rOyaempiM46dvH9etobIy0RIBQpYRzz/VxrSMHMz5eJibny4smUErlAs2Ul29wPV9/FlRCwqEMHy4h2mJB+UhI6I3ZPqO1tZW0tLH861+K+notsDKwWBGr+nhbWzkFBQUhMYmPzw29FifhAgVtQHqEBdXW1sb48VKhQv82MjKgtdVynafvt+++tey7r5ynFyxaoBobG8nJyWHXXRUrV8Ls2X1AIvfdlxzRx+v992Vvr6srK5QHlRJ8wI6ODgIBsdx7euS1V1RMISFBJhY5PwGl+kICpV3VHR26T1Ri0IJKCo6/l7FjpQrJKae4k3Hfecc9tu1KoJRSJyil7lJKna2U8obd9h2biX8/0UKzMa655GQJZQ7n4IPlxzNYnJ1BGxsbOf7441mxYgUpKSkblQeVnw8jR67B47FDsPXj9vVlD0qg1q+X5I/i4gTXcaU8wETWrRu4l0i4e239+vWMGjWKDRt0ryph+PDhlJWVRY0M1CVvNDJ29z4WxBaocEQI62hpiXe5TEeMGEFtbS1dDj9VdbW4eh1eNAoKJA3hpJOiu/jk9eZTWyuzjlOg9P8NDdEFSn9WH3/siSg0fNNNkpCqqa2tjRCoQCADaAoJ7bRpEggwc2b096K+vp7kZJk4s7NlLHFxcSQkdMcUqJaWltB9srLA4/GEWp7H8hzrx2pp2UB+fr4jvyqDuLjItIef/1wey943bAUSQoEGTgvqmGPg6aftSMn0dGhvV67ztAX13nt9oYVnuEA1NDS4fivyv3wXYgVJQDcpKZnMmgWvvirvSUdHB7W18j15+GGpgLJy5Vy6u6cBWqDisaw+lJJ9RlugZEMpNTWRBx+EY46R3+Ajj/yXK6+0n99ZbT4Y3xMiVp3Brc2AAqWU+iUwF0gBrgA+Vko5fi7svYXGtl2jgxs2xoK67z53YzpNS4uEvw62L9Rhh8mEUlgYoLm5mUmTJjFu3DiSk6dw772/56mnBvc4F14IBxxwE/n5+a4Q8ZtvhuOP/y9tbW0DFmgdO7YCmMKMGe6dbx1NV1ExcIkMLVD6Prp77I03ig9fU1paSnd3d4RwWlYsgWogPt5ybf7HEqiDDpJafhpZrdfR26tc+S0jguW4nZactjycIuOkqqoKj8fjymUCcfE1NS2ivt69QBkzBo47TgSvPwuqqUlFuIsvusid+B3NgkpPTwbKQq8hK0vaVTgtUCdiQYl7Mjtb9uIyMjI49dQbeeCB6PdpaWkhKWlY6PH1fSG2BaUFqqlpLfn5+RQVSaXymTPjgtam24LyeOSx9VfX75cH7u1NDL5O24IK55RTYMKE5a7zxIJK5NxzJ4YqiuiJ3GlB5ebm8uqrsN9+YFkFDCxQPWRmprJgATQ22haUDlxwBkkkJIilpQUK5Pfn9douvvZ2XXg2iZNPht13jw++fy20toqr+qWXpKnqyJF/BqRe47vv+oBKZs+ez/33b9z2xJZiMBbUL4FDLMs6H5gGLATecYjUEMk5HnoUF0u5lO/K229L/sdg224UFcmE0t3dSiAQCK380tJ89PamDlhOxkm0CWzHHWHyZFmuDWRtdHbWA0soLs5yHdfWUE3NwBEbbW0yyeTnyw+3vr6ekSNHUlgo4qVdFcOHR98s7+0Vy9QRSR10fTbQ0eHnPIcLPpZArV7tLjdVW1sbrAXnjorTAuUcQ2OjTJbhFRV++lMJ3qiurmbYsGER1Sfy8/NpbW0mPb3XtaLdcUfJZZk82X7/owlUW5s3QqAqK+30h76+PhobGyM+38cfj2PYsFNCrkeQCS1WJSlphifiqsUmMzOTrq6aqI0ZQayWhAR5Xj3G3Fw5OZYFVVoKhx7ag2U1hlx8V1wh78OwYcMiBOq55yT8WlNU9AH5+QeGLAftVmttbeXkk+3ISBAhnz79Q9d58juSD0J/HkVFIkQ6DFz243LYsEGCenJzcwAfcXGBfi2ojIx0kpIgEJAH6ujoCN3uDDOPj3cLlGXJl/Lmm6UYrtxXW1BJfP451NdnAwncfvte/POf8MYbdh5UT08XEOAPf4D//c8DxNPT4+Waa/pPzt5aDEagiizL+hrAsiyfZVk/Q9pfvqvEWT1E6t4OPT79FP75z8GdGwhIB9wnnoi8TVsOgxWWt9+WqDFtSegJKzPTQinfoB9n//1h3rwTIjbvV62Cr7/eGUgc0M334YeJwM9DY9BogWpsjDGDOTj3XFkd5ucTapQ4atQo8vLkfdN6ogUqfB8qMRG++gp+9jP7WFNTE+np6cTHu/PCYglUXZ09Zrlex7BhS3nuObdlpgXKuQ/V2CjWU3j49DffSFfe8CoSGsmFmsXll/cMOhoO9Oedis8XF2G1/fznkmsHdrFfXanbSWFhoUugTj/dHequ6enpob29nbi4HNLT7fI5WVlZLF48m1tuiT3u5ORESkpsUcvJSWH06NsIVl2K4Ec/gjvuWAN0h/LEVq/WvdMiBeq//5XEVE139wby86tCkWxxcXGkp6fT2tpKc7N779fng4YGHQnnCb2mcIHafXfZv5k4Ua5rF58ua6QXZl6vLyKU2ylQaWlpQYGSwXV2dobOT0qyLaj4eP38qcDhHHigVIK45BLpTwaQnd0E3ElJiZfjjoO//S0TSGDx4rGh6jaBgNQT7O72cNZZ5+DxwK23KiCfxYt3BGJ3Ad6aDEag6pVSo50HLMu6DKk++i6mnl9MvBvxzjQ2SqRZtPlezx+DLNzAH/8oK/NwgUpPTyU+vplB5teyYAG0t1sRK+yPP4a//31XoHhAgfrgg0LgdxEClZoKEyasort7JZ0DtWtFJg+lCIVwjxw5MiK6TOfyRKsqEI6uvnDVVZJUrYkmUN3dEjHnFKK6ujpKSvo48US3ZVRSUoJSyiVQU6dKy4RwcnPtKL5oAiUuv12YOzfdZb35fCKWN93Un0AlU1paTbBwfIjsbFvQ9YQe/vn++Mfg853lyiuSOnaRr0Hn/ey+ezV//rN9XKptT+WRRyLvA2K1TJs2j/JyO4ghJyebxMSHYu51gbs5I0iV71tuiS5QbW3OAAmoqfHQ3X2OKyI2PT09lAflPPfCC+Ghh34Tcu9BdAsqHB0kob/ShYWyH3f55Tdy223uc3/yE7jttgVAL+npYkH19nqIi4uLcPEppV15ToFaRny87CGsWWNXx8/PrwYuY/TohGDVcg+JiXI/bYn5/WJJdXZOJi0tLWr18ijrlq3OYATqbeCs8IOWZV0CvAckhd9mEF55BU4+eXCl62PlQIFdd22wlk9trdxHh3AXB0MKc3Jy6O2t4Ouvq/q7OyDurOZmydEIn8DsL27BgLlQTU0eoM41gYKIzdVXfwI8PWCo+W232W3OnRaUM9kVZJLyer0RAvXhh5Jl7/SpNzc3k5mZyTvvSE0yjRYoZ6CFfvxwCyo3t4i333bXXExMTKSwsNAlUBdeaPdkcpKbK48dXuZII0KjexfZx71eEc26Or2Xk0SiXmIjm/oeTxNnnHGvK/cKRKD0RxZNoCxLukErNc5lQTmFzYkuc7THHn0uCzUzM5NAoKnfIImMMJ9ndnY2dXU5MatPnH02nHuuRBFpgUpPFze6Fijn5xYuOrW1KaxZc5nre6Crh4efm5ICfX0JrjGmpqYSF5cR/F+OrV4tAQb/+Y8kXHd1dZGbm0tnpwhLfHwcWVlZtLZG/k5GjIBx48qCryOdffaBceMUqampdHR0UFwsYeaTJkFqagdwKNOny283Pj4ej+fnrFghq6a99rL7UrW0dAIpJCdLHpTfDxkZ2cH7yTm6H1Rv7668//6ZoUCfuDip0LJo0eAXxFuSwe5BRTXULcu6GBi1OQf0fWLNGvGDx3LPONFfhmgClZEhK5/BfmFqakRE9MSvc47EBfY0X30Vw+/iQFtZfn9lhEDZE3X+gBZUS0siXm9zxP6Kc1wDCdQrr9hhsOvWrSMhIYHCwkJmzZK9mGAZPjweD8XFxREuvrIycfE5Ldrm5maysrLIyXF/PllZWfj9flc7+rg4We06Iyxra2vJySniwAOl+oeT4cOHDzoXqqHBoq6uLqoFJSv2HJKSfDj0ByA07vAyRyBJvuHtzO3XJxO632+7+Jyfb2enuFPz8jxUV1eHuiZnZUXfG7Lr8I1wWSaZmZn4fA1RBcqyLFpbW1m06BjOPdf9ehsb7+DSS6PvGtTWQk+PRPppgZJwcHkNXV1drijVcNHp7KwJHddoF1/4uVJxP5G0NPugUoqMjHoOOeTBUFkhj0dc3vX1bo/FsGH2nlZWVhYffjiHv/7V/Xq++ALefTcrNI5nn5U9NS1QJSVw1VWyd+r3dwCvu4od+/238e23Up7E67Wj+N5/vwTooLExIyRQ6enyPAkJ8j3OzSXYXbuYxYtnhx4zEJgcfOyoH8FWZ0CBsiyr17KsmD4Yy7IG/iX+QNG/+8EIi16sRuvyoBQ89ZTsAwyE+M7lccrLy/F6vaE9hiuvvBK4hZycJwd8nGhJuhqnBTWQQHV0JJGUFBklBTB37mzgnajVv50493/WrVvHiBEjiIuLo6BAotmcrreioiLXyh9sC8h5nhYo7WbTRCt3VFwMjz8uBXhBygs1NDRQXJxBQkJkAmt4PtakSXb/ICcTJ8L48X0EAoF+BCqXlJTIzGo97mgCBZCY+COeffaKiEomOiChpSW6BaVFaNiweHw+X+jzjeXi0wJ1ww3TXK8xMzOT3t56urrcpbtA9lf8fj91daNdzf6ys7OxrGYaG6MLVHMzeDyycNARj06BAncTze5ut+i0tspCyClQGRkZtLW1cdppks6hkfBxD2lp7tDFnBzIzX0vVJPRGWau3Z25ublceqlUDNfvxbp1M101+wD+/nd45JGdADsQA2yB6ugQ67y3F1pa2oFjaW52fk/i8fvluxEfbwtUV1dP8LWlOAQqjeTkKgoLZe/z7LPB77ddO+FJxDNnuoNGthUmUXcLsjEClZws/vRoFhRIxYSpUwd+nLo6cdNogZL227LJO3LkSK6//nrq6zsiqgyEk5AAe+5ZD6zux4IaWKC6u9NJSYm+vklKSgZGD2hBOQVK50CBuCjeeMOdAFtYWBixF1FfL1aQo9jCRglUuIu2sbERy7IoKMgPtRN3Ulpa6hKo8nKiRrNdcAHcf/83oXGHIwKVSVJS5PuXmyv7ldpVGU5i4liamsZE7IMedBA8+qhsvNfU1JCSkuKaHLVAFRfLzKv3oW69VfKEwtGTcnt7vCtiMDMzk+7uGhISrIhIVr1v1teX4vpM5PU209QUffne0gJKtZKZmUlCcDNFC5RehDk/+wUL7CClQCBAW1tkVXFtQf3xj3DaafZx7cJLSXGH/qekjGPlytGhPSanQIXv+TrfC8vqihrF5/H0hcZx8smyCNUC9frrklKwZAm0t3cAL/D111I6XVxyCfj9MhBnmHlXly/4mMncdx/8+teQl2cxY8bxrtdYUNCHTo4uKnL/PkACeLY1RqC2IBuzd3TccTBvnjt5zsm33/bf8kCTlyf+4xNOEIFyFgEFWLToMKCT9ev7V80pU+BXv3oXWBwhUMnJUt4lPf3JAfegpkw5mmnTnot6myTv5vcrUH6/CIDTgtICpZQkneqcFJCVdDQLKifHLRJaoEpK5Dbt0ogmUPfdJxOR1mLnRn1eXqQLt7S0lLa2NlpaWkIliWLlQMWqIgGyuo+LO5HTT/9LxG1HHikNCFtaWlxljjTx8SJ44c87eTKceaa8nmgpBH6/uEzHjElxjW/CBLEEw9EWVEtLXIRAWda9NDR0RORP2YVQk133kdfR0m8eVCDQ5Io6vPRS6TgdzYIC+zOX55TnjRSoNrq73R1md9sNsrLuJiPDXXbF79+fL7/8Q8jjoStPhFtQF18sbmH9ugKBzqgCFRcnC8XUVEn/WL9eqkmEB0l0d7cjrTDE12t3z5UHdbr4urpEqRIS4jjgAGl2qPfaAgHYYw9J/n344TXApYDs8b74IkA3w4fb+QTvv79t278bgdqCDBsmE+sAuayD4s9/drcMj0V8vGT+6z2oEneHPgoL5Qu+bNnAG2OxorxAcnFycuIGFKjW1g0UFiZEvU1EJ5X162OHFba1SR2ykhLZhK6pqWGklGtGKTvQQFNYWEhdXV2o4RtI7TFnJelAIBCa2K+8UtwoeiLTk73zddXXy2SiFw/O8Ozw5wd3uLsWtWhJrl98AeedtzuwU1SBUkqRnZ1FW1vk+3PxxXDttbFdfF5vPuCLaBjY2QmffCJjrqmpiQgxnzFDirIeeqiYENqCWrRIipGG703U19eTkVFAd7eKECiwQtaSE32suzsxqgXV1hZ9WjrhBEhO/iy0/wSS73fUUfaelP5sLEvcWLpNiSw4+rjppue56CL7MUWg4khOls7Bml13Ba/3D+Tmur+7ugistrA8HmleuMMOtgWVm5vL0qWyB63fC7+/PaIGphaotLQ04uLiSEqSskapqamuMPPERIJ7or0oJTFpdvdcsaCuv95Oo+julu++1ysRt198AfHxI1i16n7eekvceWVluplpD8nJPh59VLskk+jstAND9ttv23bX/U4CpZR6Ryl1o1Iqe+Czf3jk5Ij1dMYZA597+un9C1BBgTzWQJXIv/5aVvzt7VZUC2r0aPllrVzZ3O/jXHMNXHvtaSgVF1HhAODf/wbL+km/AlVRARUVv0SpiVFv1/PMhg2xS2RkZcmP6ec/t3OLRjlip8NdbIWFhViW5epV9ZvfuN1TbW1tWJYV1fLQx5wTa2OjLscj150W1K23uic2sMPdnQIVzYJSimA798KoLj4An++OkFsnHIm0jC5QSuWiVFPE5LJ6teylvftudAtKo8ejLah33pFoxHB3XX19PVlZkoUSaQ2N5cILU1m82H0f/d6Wlva6mv2JQD3JpZfOj/o9v+MO8HqfcglURYVMuM7+WSAT/SOPyH4L2AuOCRM8rrSA9PT0UCUQp5jLvk9KqCySJj5erjvDzP/9b3EPagtKh5lr919WVhY+X31El1oRrJ6QizUpSY5pF59ToERMerGshNC5c+acQ0bGs4Dkth10kH5Ny8jMvIeUFLEwr78ekpLy6OnZK+TNCQTgqqtGAIW88opsQOkowIYGO6PdsiL3p7Ym39WCeh8oRsLNDd+Br7+OngOlGTZMfnQDVaZ4801pTdHQ0ExnZ2eEQI0bJ5PZ2rX9Z+FVVoLPZ5Gfnxfaw3Ly6KPQ2HhWvwK1fLmfnp5LUWp41NunTIExYz6kpmZw9aCcOVCacBebnlj76/brrAC+eLGswnXfomguvoYGtwXkFKjZs+1OphpnPlZKiiQaT4yi0Vq0UlKGu8LENdJO41QqKyOtq7lzZeJqaVFRBaqoqAvLejOiFJWznFA0gXrySWnvEBeXRlpaWsiC0loe7n6rr68nPz+JJ56QnlcaGVMGL72UFRGooV18//jHBq5yNOoRgVrEpElLI4TVsggWT61zCdT994vgJiQkkpmZGbKgwgvF6s/zo48m869/2Y+bnp6OzyfvvVN03n67j76+9bS0uD84ae4XvR9UY2MjSUlJpKSkuARKLKgf8/HHbl/ZX/4CO+30l1CuVX8CJRZUDyACJZVV2ujpkQXE8uVSr0/G9jWlpQ+QlGTX6EtN1dEiAeLi5L1cuFAaFKak2MqsVGRAzuaohrOpfCeBsizr95ZlnWNZ1vTNNaDvG1dcgatAYyxqamIHSIAdOTdQkm1NjZ64ZJM+XKDGj5cZqry8J+K+TmprwettirnCzs8Hny+7X4HasEHcD0VF8VFv3203OPHEV6mtnR8KZw7n3XelgkBZmTsHShNuQenxOvehZs2SxGWNnqyysrLo7ZX6hzoqPFpPqPB+TFqgcnNzWbFCIvycrq/i4mKUUpSXlzNqFDz4oLRsD0c/Znr6qMgbkUnWsrz4fJFiq62Azs7EqAJ11FHLgdMiPh8tUA0NAerr6yM+3+XLJW8sKcldTSJWnbyGhgaGDUvltNPscH/Q76MsBMK9fHqBEC0PCvL49NOMiEmxqkom29raI10ClZEh4tXRIQsG/dmEC5R+H15+eYSrWouIgy4fZB9XSgcfuN/buLh0oJNAwHYh77WXLEJ0mSMQV6oWvKysrFBovZMddoC4uOUhgdplF3GpaYHad1+xzqVwbTvwIy6+WPasWlth9eqTaGqSxd+559r9zpqb/cTHD8eyCEXxpabKe93X1x0SKPnJ7cHVV9uVYpOTpZ7a66+LuxfEAyT7U1uf7X4PSil1qFJquVJqlVLqt1FuT1RKPRu8/XOl1KitOb5vvpFJtj98PhGewQjUQAEXNTVibZWVyYwbLlDDh3tJSbkXpZb1+zjSaiNyj8I5nt7eTJqammM+hnbdlZZGWgea4uIS+vr8IfdIOCtWSM+cuDixoLxebyjxGMQt4Zxwwl1TliUBJk7/v1OgnJXBIXpPqCOPlOKhGsmByiE+Pp7XXpMfsHPijo+Pp7CwkLKyMvr6Ym8yi674SU4uiXq7tqh7eiITq23BzIkqUHqiDI+yTEmRfcrKyi4CgUCEQDU12QVWi4qKIiyo8PVIfX09qalj+eADd8fd/gRKJuqRHHVUCW+8EX6f3bj//mNYvjz8PnIZCLS6BEoLkI7kG8iCyshQriAJca9FClQgICd5PI44dWC//VYCR7lcwO3t8vt1VjKfPVv2ae3XdRrnnOMOqXzmGSgrGx9y8V16qViEWqCmT5eagLYF9SVTp4pZ1tgI8+efQHu75EE5w8xXrTqChQv/RyAQKVA9PZ3stpvUNZQ8qFLefdd+4Z2dMmifz957qqsT8TvkEIko3JoMSqCGarsNJT0b7gMOAyYDpyilJoeddg7QZFnWDsCdxEg63lIMGzZwmHl9vR0aHotddhGhGyjUXAvU0mDs9cQw31JaGowb9xA+32cDPo7PF5mkq8nPh0DAS1NT7BCfigqx0kaMiF4XprUVLr/8AuBXMSP5tHDk5opADR8+3OVynDpVwvM14RZUZ6fsJzhddM4Cq9F6HYWXO/rlL6XWmcbpZorVzE/nQt19t0weUQpmExcHyclvM2xYdGtWa0tXV+R74xSoaHtpN9xwGHBLhEApJWJTXS2KHU2gtLXktKD6c/G1te3Gvvu6K/c7BSp6mHkeixZ5XfXpvF4vycm9UZ/Hfow216JJG2Gtrbo9iQhUb6+MWWu3tqCysrwRUXxQzjnn1DLaUdDNsmRjSiwmG2lL8bbr+5GSYoeZ64XBk08Scl/KezGN//zH7Re87jooKzvYVU5JHk+i+CoqrNAemgjU0XzySXLo9YEUhAZ3FF9PTyA4djnu90NWVjKwEK+3jQ8/hF/8Qlv9Yvp//rk7yvWII9zFgauqJKVjzz1lT6uubvMEfw3EgAI1xNtt7AKssixrjSVlfZ8Bjg4752hAByI/Dxyg1NaLS9EC1V9wg88nkUjhYbw33HADV199NXV1dWRnS5uE8IrY4WiB+uqrrygpKYnIyQDIzx/Lhg39t9U94QTo7X2tX4ECyXPqCa+CGaSurg9oYdiwrKi3p6eDZUmBylgCVV8vrpKkJHcOlGbNGtkM16v31NRU0tPTQ3tQ0VpdOC0oXeDUOY+HC1RLi/vzcwpUeLkljc6FamyUiSI8mg6kooJlHc3ee0e3Zjs6ICGhm46OyB5XtuDmRrWgKirSgOSoeWqPPAIHHiglzWNZUOAWqIkTxZp1FnLt6uqio6MDj0fei8goPh8ZGW0R+0lSnqko4j5yPyt4jvu4LVCtES4+kAVAQUFByMW3887yWnT0ZnOzVDPJyvJEEahVnH32Khxbm/T16QG4F1cVFWOA/V0WlBao8F5QGllAdNHXF+dyBXd3g9/fHhKoP/1J8pFSU1Px+/3ceaefnXeWc6Ug743cd58Kjk+O+/1d+Hw+l0D19vpRyodSEv17110waZIHmMnkybY1nptbB0ippbQ0Kb0UjnPxJ++jVIgvKJD0DsuSIJEtxWBLHQ3Vdhsl6HdYKA8ei3qOZVk+ZFkX8S1SSp2vlJqnlJpXN9hqqoNg2DAdERT7nNJSaWew337u44sWLeKWW25h9uzZvP/++zz7rLtBXzTefdfi/PO/4NVXX+VAZ2y1gwULbmfZsv43xm68sYPe3odiCtSxx8Kttz4CrIq5D3XccZ8AhVF/tCCr+ZwcC8iPWU3C2cfJmQOl+eQTCSd2ltlx5kI5LTCNU6CUEuvUOcc7Bcrnkwn7hhvs2wdjQZWWllJWVhbav4q2JGppaaG7uztqiDnA3nvDH/5wL4HA+xGNJgsL4aSTyoAVEQLl90N7uwdojOo6PeIISE2VyIXwz3f0aEJlfAoLC2lpaaGrq4vERLEenMEB+rH1VOA05HQV8F/+8mauucb9/K2traSkFEfcByA7W6akcAvKFhW3QM2eLcm4Y8bYAhVtP7OpqYmsrKwIF5+IQyoVFd0uV6xSzcAV7Lyzz/U4L788GbipXwsqEJD38b775Hb5fMTd7XQ1d3dDX58tUL29UlEmOVlEsb3dFyru2tHRQVycP2Rx2tZLHx0dHa5E3d5ei7g4eTGzZonI6P2+lpYWdt9dhOukk+aSkHArAJdfbo83Pd3OfNe5XBqnYB1xhGxhxKpYvzkYTL1tV7sN4GdKqTuQdhv7M4h2G0qpPSzL+uS7DXXLYlnWA8ADALNnzx7wNQ2W0aOlbEh7e+SPcSCef/555s2bx8knn8ycOXPweJqZNWspl122jqysLHp6emhqaqKyspLKykrWr1/PvHnzqK6uprCwkKucIVIOMjN7aGjIJBAIRK2R5/dDeXnsHCgQq2b48GTAoqmpKWqYtKzeu6NacZphw+Koq8ujomJh1Ntzc2Wztqenh6qqKlcEH9gCER7JpwUqOVmswTFj7NvDeyh9/LH7ObOyskL319rrFLja2lr22muvmM8P4uJra2ujpqaP3NzoQSKyv/MADz98JJddFvUUR/PBJpcrKDMTTj11Ps8991WEQInFp4AmGhsjTe6FC+Gjj2TmC99jnDvX/l8LZ01NDaNGjeKuu8Slqtc9WqACgUwyMtwuIqUkujBWHlRCgrgLwn8TuidU+N1Gj4YDDljE229XuQSqsFA+XxAXn9/vp6mpia++yuXBB2VPJzfXTsy++253Ww15T8/lpJP2d1mP3d3NwG3stNNZrnFIs8MOl0AdeCBUVVncfntDqFDsunW2Va8tKJCmhTp4orsbfL52V5g5QGKifGYdHXYNxvb2drxef8i1ZxeC6aWjo4Nf/zoz9Hx9fYGQQH30kcw93d0FwOfMm9fCkiUSoNTd3U5SUoCUFHjtNbv/U1eXvcf7f/9nv/Z//UsKCrzzjohhUZH8fdb/bsF3YjAW1OZot/HKJoxtMFQAzhjm0uCxqOcE988ygUGUb908HHsszJ8vVlIsbrtN6r2FJ/IBzJ49m6+//pq77rqL+PgmvvxyA6eccgqHHXYYxxxzDD/96U+55ppreOqpp1i9uo7c3Af4v/97nhUrVjBhwoSoz5ebK1ZLrMCEL76AiRPHAAfGFKiODnjxxZ2BvWJaUM8/Px04J+oeiaagQBEfXxzTxXfnnRJBpFu5h1tQ0VxsznJHEyfKCtsZRdfc3ExaWhreGP1QnBZUuIswEAjQ0NAQmiSLi6UkTHg7DR2cUlXV028nXUikoSF6+ZCnn4ZHHjkYUFHf45qaTiArIhpOn6pUc1QX3y23wLPPziEhIaHfz0YvOnSgxA03iKWv0VUkfL60qK8xMzOT998/kKuvdh+XSuZ97L13pEDl5SUxcuQvQqKjmTED5sx5GXCHmXd3S2rFhg222NbV1bFkiRTx1ZarbUG53a3OKD5nmLl02R1DZ6f7s+nriwc6XcJ78cVwzTUd9PX1hQQK3GHm0EhWVrtr36ary6Kvr9UVZg6QkJAevN0WKHGl2gI1cybMnftP4G06OjrYZx9xv1qWRV/fi8yZ8xYgc8tvfgOWlQHsQkNDL3Fxsgh98cUL8XonhZoX6uaOPp/9mtvb7T5T2ojff38JmNgaDEagNke7jS3lBvwSGKeUGq2USgB+DLwcds7LwJnB/08A3rHCHfrbmIoKcWEkxXgn09LS+NWvfsXOO49ir72O49tvv+WTTz7hq6++YtWqVXR2dtLY2Mjzz3/Gt98eyZQpx0dsvDopLvYABax19olwYEcKNsSM4vN64dlnxwL7xBSoBQumER+/f9Q8Ks2PfwzDhn0wYD2+aDlQENuCcvYyCkevpjXXXiu1DjXRBEoLYWNjI4FAIDRJxseLmyl8L0UL1K67ro9Z5FestEba2qJX2vjiC5g/vwiwogrN1VcfDtwTYUF5vZK4mZ5eHbOieWdnEgUFBYRvx06bZltR4RGR4S03tCv8Zz/r4O9/jxx/ZmYmNTUjIkp0tba2MmrUAj74IHJvLjs7E5/vP4R9zHR0QGVlM+np6a6csZYWKfL6n//YAlVbWxuaTJ1h5tnZ2Xz0kbiz9GJQC5TH43cl0opAfcEDDzh6rAA9PV7CLSiITNKFcIF6kssvvwunN3fevA7gzgiB8nrlemdnIOTia29vJz4+EHLxxcdDcbEX6KGjo4Nly8Td3dPTg2W9w5w5iwA7ii8pSdS3q6sdj0ciS+vrxxIXN9Ll2kxJiXS15+TIdzFaT7MtzWD3oL5ru40tIghBl+MvgdeBpcBzlmV9q5S6QSl1VPC0h4BcpdQq4DIgIhR9S9LSInscjz8e+5yBcqA0BQXQ0OBl8uTJ7L777syaNYuxY8eSHCwIpqMFB2o0NmFCDpDCl19Gjxntr5K5JjER0tL89NcTqqsrNWahWM1558HMmR/E3IPaay9JaNQ5UOECpYXDuW04bNgwmpub6enp4aab5BznyjVcoKqr3W4KZ0+o8FJFOkrMuYp/5hnp3upElzuaPv0zV58kJyKiDXR0eKJGREkFC5k9or3HiYmdRAszHzlS2rwMG7YmqkBlZ0NvbwoFBe7PtqtL9hT05B4uUOEtN7QFtdde6a4kXY3k/zRHzYOKFtghY8umvn4mH37oPv7738MDD/wpYsHkDJJwljtqa9P9mOT2pqYmcnJy+PprqUihx6QFKiHBHegjofAd9Pa63bNdXXGEC9TvfgcTJ8qCJJoFlZCQQHJycoS7s6CgDagPCdTEiZKyoOv/HX54WagJZEdHB1On/iVUEWXlSnjuuenAcNrb27nlFkmFkGi/YfT0iJtOC1RiYnJw/O2uPKienhkcdpg9prQ0mUQefth+b198UTwQ/ax5txhbpN2GUmqtUmqN/gMygv+vDV7fbFiW9V/LssZbljXWsqw/Bo9dZ1nWy8H/uy3LOtGyrB0sy9rFsqzN+vwDkZYmvYhWrox9TnX14AWqvzwoLVD9hasDnHhiNgkJv+YbHcMa43GgLqYFBdp6yY86eXZ2gt+fSFpa/9GCgQBkZ4+jvDxSoHp7ZX+ooQHWrl2Lx+MJTfya9HRxoTr7CjmrSdTWyuM4V8fhAqXr6Wm72tkTauxYCQfWuqitBuf7cvPNkQ0Ji4qKUEqxdGkTsQrHV1ZWEh8vahCtioi0ipdBRXuP4+PbUCo/VNk7nFg9obKzwbI85OWNch3XT6Gtwfz8fJRSLoFyGg5aoL74Ipevv458fmla2Bg1D+rbb38e1U2UnZ1NT88N3HprIOw+0mrDuTAAQtUSnBXN6+rqIvo7aQtKH9OBEomJiSiVTny8W6Da2tpQqovubvcU+corirS0e1wC5fFI63TwkJOTQ1KS1OdzrqVSU/fimWfOCOURdXbC7bd7gemhPag5cyQyrqRELMSxYys57jg5v729nby8ttBe6sqV8MwzY4BCV5CECNQfuPfeEwE7zDwQEEu5o6ONOXOk+K9lKfz+HNfirLZ2JiD3cea1hZdp2lr0K1BKqX8rpfqtEqGUmq6UeiHs8Bxgv+Df/kBH8H99/AeDxyMTeX+5UNXVA4sKSF5Ff3WxBitQs2fHsdNOn7B06byot9fWQkJCJ9nZqTEnP4DCQkUsC0rvCWVm9p8scf/98I9/3E5DQ1xEuLp2r+Xl2TlQ4ftGSok/3ulicwpUeJkiiKwAnpdHqOo4uMsdTZ4sq3c9L4a3HIfIgrUgq+aCghHcfvuvI1p9ayorKykoqI2IlHK+/rw8cY9Ge4+lEWRkhOT998v7kZY2Ouo+o36vMjNHuY6HC5TuJdafQOXk5HDeeR5XcIUmKyuLvr76iIg8iQwsjNrIU6pJtFBf746ekzDztgiBUkpW+m1thKJFa2trg/l+co5lWTEFCiA19b/MnPl62PO14vF04+hbCYglkZPT6LKG7L2rFHJzcxkzRkKvd9vNvl9yciHl5dNCv9HmZrj99nxg1wh3fGrwAZcs8bryoNraduORR+S6veiRKL74eBEoifb0hgJWdKmjzExITPyCvr46nntO9ydTofO++MId5HLeeXK/K66A22/fdgVjBwpweBx4XinVBrwDLENq1mcAExHxSQN+47yTZVnrndeVUv7wYz8kBkrWPfJI+8fUH2HGQwT19bJiCt8PCaevD4YN24/PPnsp6u377w9vv/0ygUD/SldQEEdcXE7UybO1FeLiWsjJid7bR+PszltVVeUKgnAK1Nq1ayMCJDQvvCCv6eST5bozWTeaQInw2Pnczn2s9HS3QCUmlmJZzlJTkQKVl0dUC6KwcDI1NbFbbVRWVjJ2bGz3b1YWlJZ6+OgjT9T3WKlmLCvyi9PQIBNgQUESa9ZEWlCHH27h8ezPyJF7uY5rYytc7PV+3t//7l5J19XVkZeXx/r10b9zmZmZ9PWtYdQomexkld8XbFiYFjWqVVc0b252W1BS9qk5QqDAbvseHx9PTk4OtbW1oZBpkHyt3t7emAKVk/MZI0cmIVvY+vY2vN5el0D5/fIeJCbu7LKgbIFKjRmxmpERHxwLrkvoCgnUm29KWa+nn5Yv7N//PpWXX5aWFx0dHZSV7c3vfidFpcPDzHUlCbGgvKE+YNdfL/t3O+4I48efi8djF4JNSFhIRkYibW0ysYwZE+npOf102ZfcVvRrQVmW9QIwHrgWSAfOA/4QvEwH/g+YYFnWQKla2zJXaptTUNC/QN1yi9tFFYsNG8R/Hiue4IYbZBUcJXLcRWUlvPjiH6mu3pXOzkjv7THHQHb23JgVtjXPPQfFxcdEnTynToWsrDHsuGP0QAyNPd/kRQRKaKskN1cEarQz1d/BX/8q0X4apwUlbjL3+eEuvjFjJActWk+o665zV+/Qe1DOCu/RmhbK8fGh8UejqqoqVLYpWtjOf/8LDzyggq3Qo7nq3qOk5JGI442NMmnm52dEvV9SUhN+/3sUF7vfmNRUiQRzRpyG1+NzBjXU19eTnV1MT090gZJUiNv49ttAaMJsCypDb29KPwLVQnOze8pobbXw+RqjCtTjj0ukGriTdTX6++kUKKfwJCaOp77e/aNpa2ujuPhJLr/cPtbRIVX1fb79IvKgBBGol1+W37yzkWZmpngi9M/NjtjtDglUfLxYRkrJflFPD64w88REFTPM3O3i84YWEmPH2iWXMjIyaGlpYc89dWv5w5kxQ7J/zj/frlfp8dippeElp7Y2g9mDsoL7PD+zLGvX4H7PrsHrrw4yIu6pzTDW7ZY99rC/JOH4/YMvGbJ+vUQgxaqHpVT0igXhODvirlkTuSVXXQ3l5XURvaTCSUyEnJzoBWMDgQBNTU0xk3QjxxJZTSIxUZJV8/IkByqWQOXnuwVC70VUV1fzox/B0Y7aIpZlRXSh3Xdfye0YG1xcOgUqWqHYrKws4h2mRF6eiEJ4r6TMTNkwiLaotiyLyspKUlMnkZrqbroYTk5OdCs1Pv59xo//b8RxXa4oJyeHlpYWV28sgDVr6oHTgLGu47Nmwf/+565o4hSojz6S7qx6cqyvrw+5CWOFmVuWFRIlsAvF9vQkRxUoee+baW93T00nn9yFZT0UVaD22UcaMYJd7uhnP7PbRzgFauedxdJwBgZUVNzHF1+c73rM1tZWhg9fyI9+ZB/TopaWplwCNWUKzJjxIcnJARITE2lpkaAdpzc6M9NtQdkC1RWRB6UFqrdXkZAgv6XOzk6Sk22BCregzjlHwuptgRKB//hjEfC33oKFC5+mpiaHigqoq5PitdnZXnJyJDfu22/lEf1+21UTHqyytdkqxWIty7pwazzPUOWGGyI30TXz5kl79ddfj367k4EKxv7ud/DQQwM/TkoKJCf7gXxWh/dCAKZNsygru2xAgfrwQ6ipuYW6ushSR08+2YllPU1aWuSE4sQpUOGRfLvvLk3UEhPXAcR08eXnu6P4EhMTyc7Oprq6mmuukVp6mvb2dgKBQL/5P06BamyMbLURHjhy8cVi1YZbrmlpI4LjcVeBAFmhd3R0MHJkVrBMTvjtIpwvvRQ72KGpqRulJkZ0atVWo3Y3hYdEr1rVCDxBZWX0PlNOtEBZlsXXX8t+hNbK+vp60tJkMovl4oOdOfjgxNDkpyt6T59ezy67RN5HLKg7ueqqN13HDzqoAngiatDOp5/aUZTagvrgA9uCcQpUXJx7r0VIwbLcBQPFxTeGr76yj2mBysiIc+1B7bQT7Ljj3ykoELekXXbLvm9ubiIez8qQteW0oHQemxYoy0pEKUVvryIxkZCXIzExLhRmfvrp0NrqB9bT3t7O5Mli/YpAPcAFF8gX6skn4bLLJGqxo2M4bW1Szby314/Pt5DW1pGhfa6bbnK/K6mpcM454e/V1uU7CZRSyqOUGquUKh74bEM0dLpOlJ6AEQwkUI88IqvcwSCPVcCqVatcx30+mSwDgaqISujhbNgANTWHUl8fGeLz+ee9wAnk5/dfPDAnB373O4uEhK9j5kLpHKj+LKjWVlyFR2VirYlIfnaWOdK0tUmriAcewHVbLAsqfBWfmysZ9eEbydOmeYGriIsrJxwtxmPG5OHxREbxNTSIOIvYRBeourodeeutuYQHYx5wgOzHxapo3t2tc8TcqnLTTTBqlLv6elFREX19fTQ2NjpywQg1hRw1yuKTTyLLdIF+H5P44ouk0PdcT+zXXLOS88+PvI8I1HoyM92u4a+/bgHSogrUnXcScsVpC8oZxecUqO5uKZT6v//Z9w8EUvD73QLV2trK+vWnuiwtW6C8LtG3LKipaSYnp8B1nlOgiooU8fE7hvLtdtkFrrvuL8CHIYHS7eN7ehSpqan09SlHJXNITvaELCiPB9LTPSQlJdLR0cHKlbKYaW1tB97k+OP9ofP8ftu6b29vIS5OkoBhCn5/lsu9XFxszwft7bE9P1uLTRYopdSeSO27FUCZUqpZKfWuUup2pdRpSqlJW7Mo61DmjTekZXk015xuWxSjHJuLrCzxU0cTKPmRDC4aEKCw0IPXWxJhQUm4tQJqBrSg9Dzd1BQZa1NV1Qc0kJ/fv4vP64Xrr1eUllZFCNT118sPWScUx7KgotXDGzZsGOXlnSQnw7332sftqtZZoWOpqdJptizoenf2hAoPsqitrY0QqLIyifQLN0Z32y0LuJm2tsj4IB14UFJSTE5OpEA5O/HGEqjOThG+cOvroovg6qvtqLbwSD7pFRYgEHAvHiorZaXttASduVDO1iRtbW309fVRVJTJ7rs7LWGbaC03tED1lwcFY3j99dGh6D/LglNOmQVcE1Wg0tPtoIeCggIaGhpoa7NCAqXfu+zsbLxeiXJ01rT0+5Px+Zpdj9nW1kZKioWzBKIWnqyseFpaWkI1/xYsgNdffwnLOhwgIg9K7pNFd3c33cEVU3w8WFY1SvWGovZyc6Wtxdixkpy/++5PcumltkAdddSKUBv5d9+VSLyUlDw6Ojr45z9l77itrQsYR1WVfLbRBErC4sVHWFVVzF6OWJnUVHmuG2+MeJu3Cd/FgrobieD7E/A74DWk1NClSPTfN+hv5w+c5GT58UczEKqrZeUd7QcejlKxc6Gam2VvYDD5VCDuwDFj/hlhQTmTdAcSKD1XtLZG9nuqrbWA+n7r8DmfMzd3JuXlbktj9Wpx3a1du5aEhARXHygnp50m4uwUeYk+E5PK6c3Tk7Vzbywuzh0q7uwJ9fvfuwtmRrOg6upETMMtmfj4kcCIiNcFtgVVXCwCFS4yzhJL0QRK9iVEUcO1S1uSsSyouroaoJnubnf7B2erDY1ToJytSXQOVE/PDjz2GBFuRtCLgGgCNZYDDpjFC+HJKUh4fkLCbvznP4eHFgw9PeDzxRFeKFaTkWFXOy8oKMCyRFiiWVBer/wetaBZFvT1JdLX1+x6zLa2NtLS5HXpyX3WLOlaO3VqS/A5RL20ECUnyxs0ZYq44JwRj1lZ2cB73HuvCNSiRfDGG3uSkjIqVM0jL08Ktu68s4SaZ2fPY+edCT3PsGHxoUjeefOkmWFKSlooSELG3Q3cx2WXyQJAVznXr6Gvr4sDDvCzww6twetJOAvKrFwpWUVh25bbjO8iUJOAP1uW9X+WZf3BsqwfWxLzmg0cgISeb6kafNsVel6NViyhulq+mINNhPvyS6LmnGhLbLACddhhMGtWU4QFZUcbDt6C6u3NpDcsG7WhQQH1AwZJgBSgLCu7ISJgQ1cyX7duHSNHjoxa2BZkIioocK/8hw0bRn29+C6cQ4gmUPocpwWmq0n89Kd2cdRAIHoX2lgtNx54oBj4rF+BKioq4swz4aCD3Lc7K1hEC3aQwIMG17mazEzJmYslUNXV1Xg8rbS0uB0cgxWolhZboNasGc9ZZxE1GVmspObQfUC7WLNpa/MQK8UuPd0fPFeux2q1YZ8v7qhAQIf/JzBhQncoUbapqSlUvFafrwUqEIDDDnuRvr5/hVqa9PX10d3dHRHxl5QkCa7DhqU6XovtytMCdeyx8I9/uMeYk5MN7MLSpfIZLloEn39+CGlp9m9MxFL+0tLSWLNmHMuW2RZURUUxN96oq6ATfO6EUJg5aAvKS0KC/Bi0BVVcDJMmrQc6+P3vGznySFmYJiTIhtx777nHq+vybWu+i0BVAxF9KSzLarEs613Lsm6zLOvU7/D43xv0yj6aQB10EDErWcd6LO2rdtLSIivJwbr4KipAqUNYt249fY4wwrFj4YADXkep1QOGmefnQ1JSL5AQEWUmFRJWDsqCkuaHuVRUVNDlWIrX1sptq1atYoyzHHkYzc0SsfXll/axwsJCOjvljRqMQOXlRTYtbGjoZMECezLTdfjC3UyxWm40N3vwelsp06aAg/Xr1wf7UaVz1VWSGOkkOVmSQvPyogc7SLBBU3Bc9v26usTiyMzsX6DGjbs8ok1CU1NkNJ5ToEaMkMc+4ww7H8zvz8Tjid6nTCyoNoYPrw0Jn1hQWcHbI+8DkT2h9MtOTu5z1eHTSF8xERL5bHqZO/fT0Hva1NREZmZmaIHjFCiPB/baaxmBwOehRHEddZiRIedrN9/ixZLmER+fFxyXDCwhwRe8lBcZLa5ZPosuWltFyfXeaFqaHbHR0yMBU7fdJhbUF1/8hr//3RaotWuHcd11Mh79k01LS3RZUO3tPcTFJeL1yuLjssskR2///eGqqz4AqmhpaQm+xtcZOdIfGrMjNXDI8F0E6lmkMoRhAFJS5McYTaCOPx5+uxHVAV95xd2bSLPbbvKDnjNncI/z3HPw9NNnEQikh+rcgeQEjRjxLEVFcTGrfWsSEuChh54H5kZEih133BPAucE9hf7Jz4eeHpnhnAVsZU/NYvny5RGdgZ34fOIzd1bZkIlVhMQp2rEE6qCD3Jn/WVlZVFRkMGsWoWKnOgcqXKBSUmR1He6ma2yEpKSOqBbU+vXrQ3UFLSuy4+5RR0nob0FBdKGRid7P+efPc4VCO/euMjMzUUpFFajRo7sIN5D3318KrzrJyMggOTmZqqoq4uIIWT069LyvLyNmvyuxWAL87GcPctZZ9rjj4+X9iyVQ4T2h9NonKyt6Rstpp4nLKznZTqB25kI5O92CfN/0eDs7oalpLJAREiYdaThrViP//KfdK+zTTyUYIz7evWDo6ZH3NyFBXtBPfhKZ3Cq/gy5aW0XMtEBlZNhmpNbe7m6xoAKBeBITbRdfWpr8Hnt75U8CJVKDhWTlvm1tXXg8CSHBGjZMrD55LvmNnXzyMO64YwpwKMcdJ6J8yilwYSjWuocxY4ZGPe1BC5RSapFS6h9KqV8rpQ4GHgEmKKWMlTQITj89ekRMRYU7+mwg3n9far/Fyj4bbFiKPccWuNx8GzbAkiVdMQMSwtHBBuEWVENDA5mZmQOKHIiV0N6eBCjXWPbfHyZPbqKjoyNm6xCQyTguLrJgLHzDqadWRghUSkoKSWGl46+7TjqaOl9XU5P86u3oyegCpVT0ZN3GRkhL64tqQTmbL/761/0HycQWKDj++EZmzXI/J4irzuPxkJWVFVWgAoH9IlIS/vAHSeB0vzblyoW65hqppqCDPLq7U2JWykhISCApKcm1eGlubiYpqTg0xmiE94QaPhzGjbuX4uIo2dCI+2qnnWS/RT6biVxxxZxQRGv4vuGnn0r4Ncie0m23nQzsFxIofTl2rPSa0i487eorLEwLjk8G2N7eANzE9OldofPCvdFaoNrbYwuUUiJS3d2QkpIOeF1RfGlpomA9PfL7T04WS6ujo4OjjpK5weerIy4uPiRQn38uqQH/+Aecd94hQD6trc6gkiTy82WLwP7pJRJlTbVN2BgLqhc4EbgV+B9S9mg88LhS6l9KqZOUUiP7e4AfMvfcE+nG8fulEdv11w/+cQoKxI0TXiPswQcJrVIH+zjB/1yBEr/7Hcybd1u/LjUnL7wwDbjFJVCNjfDMMxeSnHz8oB4jPx/8fgVkugTqH/+AWbMWAPQrUDrIwSlQYkF9yQknfO5yicZqyx1OVlZWKPhjIIECCZB48EH3sbo6yMkJRFhQlmWF9tVAJuqODvdC5be/lags6F+gWloKWLTIvl94Pb3c3FxXFF8gEKCmpoaamn259lrnmGIvepwC9e9/S85eZWUlOTk5NDd7YgoUyPv473+fEMpFk0rmdZxxRuwSUAUFXoqK9ueMM+R6cTEkJT3I8OHRB1hZKXmGVVXyXik1jLKyYaH3s7a2NmbRY3t/qyVCoDyeHN56y1542FXeZXPKbslSD1zNPvtIVF9HhzvEHLRAzSM1VT4LLVBZWe6FUlKS3JacLGab24ISMevtlcVUW5stUMXFkrDc3d1CUdFfQwuNN9+UBVBTE9TVyXNZli/YcXkdK1emh4rF3nGHPY7e3qERgD1ogbIsazYStTcd6Q91N/AZEqZzLPAMsEYpVa+UelMpdYtS6qTNP+Ttl/BKAzU14kseMWLwjxErF+rjj4nouzOYx0lIKHWJQlWVH7+/irG6rMIArFuXBRzqEqjKSqirm0BKygB9P4IcdBA88IBFampCaCx6slwerLXSn4sPxIKJtKDyWbfOveqOJVD33y+Tip6wsrKyaG+X8KzBCFRmZmQC6G23wT77rKW1tTXkNgIRmvb29pAFpYfjNEIXLxZrFqILlH68u+7awZVPVFQEV15p13YMjwBsbGzE5/ORmxvner7OTnHf3XNPxEtz1ePTEYdVVVUUFRXx+OP9t5LJzMykuTkzlGLR3NxMcfFSHnss+l6qvB9ZdHR8FQpSqKqCysoccnOjf59WrZISREuXitWYni4Wmr5/uED95S+2O8vOt22NcPE1NORz0EF2SHpHh7xH+flZodci5zUAmSHXZWenO8QctKfhFPbc81VAgliGD59MZqZbybRAJSXZAqUtqIwMWTA5A1K0QK1fD088AS0tFkVFi0Kufv2dtO/jJxDQjQ9H0tOT4FqYTJgwtEqmbtQelGVZfsuyFluW9bhlWZdZlrW/ZVk5wGjgOKRO38fABOAK4OnNPuLtlBtvlAnQmQSpJ6DNIVCD7SkV/jh5eVNCIgCwbl0fUDFogSos9BBe0VxHFObmDq6G09SpcN55ih12KAoJ1LvvyqT/3nvdpKWlhdqPxyI/3x0sIC6dJ7njjkNd5zU0NEQN3EhKwlXRISsri87ONBISLNdEp5SKev9nnnG3xwY49VTYZx9ZiW7QHzaRva2c+UWa2lp77yyaQOn3Oz/f43rd48ZJbUcdwRYuUNoSys9PoKfHDg9vapK9vGhNM0tLS0MdjXNz3QJVXGyXiIpGVlYWSrW4wswzMnJiWmv6Pq2tJ/Hii7Kie+ihAA0N75GXF/0Lrj8fvQbIyCgJHdcJxU6BWrCAUIh7NAtKC1R+viiotpza26WUmDNPDnRE43zuvlveiGgWlMfjITMzM/S5eTzQ3l4d0Q35V7+SbrVZWfEkJBzOscfaAnXccQk0Nkoww8MPy36YFqjPP5dthKamNHy+nUJVNLSrz7bO/VhWH3198gF8+61ybT3ExweAWs48M0rewDZgs5Q6sixrvWVZL1qW9TvLso62LGsEkAccPNB9fyhkZcmXxLnK/y4CFR5aXFk5uGRfTX6+ZJ7vvnsdCxYscDxOHFAxaBdfaWkikEdDQ6RA6dIvA9HTIxF4o0fvxrfBmjg1NTJ5rF+/iMmTJ0d0fg3ntdeknp4mPj4ej6cEj8f9RklFhEgLSkfi6c9HC9zcue2hfb3a2lry8vKidgj+8ENZmWtaW6XDaUmJFLZbtmxZ6LbwxGNnhQaNM+la7/OFW0IAhYUJrvs1NTmtgtgCVVycHDofIlvbOxk9ejRtbW3B2ooyzsrKSoqLi/nzn+V1xkLq8TW7BGrZsuuYOTP2fUSQL+PRR2WBU1nZDXRSUhK93IqzaSFAaqoIWXq6iIjP53MJlDOKz36vbIHS7lMdTq4F6uabxbKNj48nLS0t5DqtqakBOvD7Rd1PPhlX4IomELibl146BYCnn7Zobr4sotXGVVdJ2kVmZgq9vf+jtNRPe3s7CQkJpKXFB/cWpcrI88/bAqWDJDo6elmy5LehRNtoFtTo0UvJz5ekve5ud9+nb74ZDRTQ1radBUlsLJZlNVqWtRFOp+83ukqPMyluUwRq+nT5Uh1+uPv4hg0b9zher0SK7b33GKqqqqisrKSri2D78fIBXWoasaC8VFfbNYW0QBUXD86PXVdHsC7bMWzYsIGmpqZQPtbSpe8xyxkFEIPk5MgAEaUKAHcZ+VguvvBcNZnQFrHXXvb+UX97GXl5MtlrN+5XX8Gee0JHx0SUUixxlBHRYjV+/PjgpVhfegzhVUGiBTs0NjaSmppKfr6HpibbMr/mGrdFE0ughg+XyVcLlBbmaGkKWkjXrl1LQQF4PBbV1dUUFJRy5ZXuhUE40vyxIRSR19zcTCCQFWFhOJH9wxbq6rRAdQFNMfPywi2o3FyL5OR5ZGREd8ump4uVEwiItXLXXfVAc4RAFRVJMIQWqLQ0+zMqKioKvZdVVVV4PL3BpoVETRsAiIsbTmOjiOdrr/mxrDMjBKq1VaIXPZ5s4GRWrOiio6OD1NRU1q2Tqu2rVsnWQHy8RPtJDqIEX3R19QF2kIQWqLFj4cQTISsrjcmTX6Ck5GnX7a++6h7rCy+E+Si3EVtMoAxutEHizEXdbz/Zp4hR9SUqXq8djqrp6ZFSSv3EEUTl88+hu1uyUL/66iuUgj32mEth4ReDCiQAac3g9a6kvt52CaSnB4DPQtFOA1FYKIEOyclSvHTRokVUV0N8vEVb23pm9rfcDvL++9InR7usAgHw+bLx+ezyHYFAIKYFFS5QYkHtzbx5dq2bgQTKsuwJX7tgR4xIYvTo0S6BWrJkCSNGjAhVsR4xQlIHtLD09Ehy8PTp9uOHC40WWskhsy2g2lpnAAzBQIZm/EHl1AnCJ52UTnm5tBl3jjfay9MCtW7dOm6+Gb74oo6+vj5SU+VL3Z9rWbrqLmTPPeW6JByn99uzTNy5zTQ0BIJj66U/gdIWlBaoWbPWk5h4kEugnFF8zgTccePgJz9RgM8lUBK9KJ+PtrYeflj+5DXb+3JVVVUkJvpClkhnp9uVr0lKCtDXJ8rR3t6Hs9WG5uCD4cc/ht7eYcAzfPRRH+3t7aSlpVFVBbfeKnOIFihdJsnv7w4+dx/OflBnnw3l5eL+e+45KCjIoK6ujubmOvLzP2CHYL3gQABmzNCjuISPPhoae1FGoLYS2oJyFm7YaSdc/WYGy003SeioJjFRosguumjjHufOO+HBBycRFxfHV199RVIStLTcz+zZMXavo3DCCTBx4nH09a1xHGsCdic3d+AkXRDRLSkBn09UYuHChVRVQUaG/OIHY0GtWwePPmqXk5IJ20NXl/1Dq6+vJxAIRFSCAHF5/vSndnCBCNHjPPSQPbENJFDyHHLp7G48efJkl0AtXbqUSY6eFtpi0vqTlCQr2lMdCRzRgh1ycnI44giJrNMWSbhA5ebmhlqMgOx/ZWdnU1ycTkmJvYIeM0aKqEZzE+sivbpor75MS5M3ayCB6u29l5degt7e3mDzwLR+I/9EoOx9KxGq2AKVmCgBEr/4hVzPz8+nubmZ3t7ekJXj/Nzy82VR0NUlLrslS0ThnAKVkZFBUlIcr74qOUIgxZifeMIeo1OgUlICocja7GyphRhOcrLC50sIPpcP6I7Yz9RBEvHxIo5+f3cw8jEzlIPW0yMClZAQKVAdHb2AJyRQ6eny29Jh7/n5+bz55q9ZvPj3zJkzN1QM98gjRRiFJeTnb0TuyxbECNRWIjlZzPOdd7aPff65e09qsLz9NvzrX999TMOHQ0VFHBMmTOTLL79k0aJWlizJY8aMKH0Q+iE7OzsszFwm0sFaYXos9fXJlJSU8PHHH7P33jB69Lukp6cz3WlKxEAXXtcR3QkJsN9+L9LR8WqohI2erKIFXHg8sjrWJYdycvKBYpKS3O6xaOIGIlAJCfaeRm2tPGZOjgjU8uXL6e3txe/3s2zZMldHX5AJ8+abY7++aAKVm5vL2LESjq6jxpzBFfp++nyQYI0RI0bQ2yt5T7rEzW67yR5arIoQmZmZrF27lmXL4MILi4HpJCZKJEZ/AqWLpPb09IS+I93dyYOyoCQUGnbZ5Q2UuilmZROlxBLUib/vv38M8Do1NTWh4JQRDv/3mWdKb7WCAvjjH+Gcc+JJSEgIhXO3traGAiEOP9z2frS32wuBcBff5MlfcPHFEmjS2xsZJAFyzO9PCD7W4AQqEOiKEKjeXplPcnJsgRo3ronPP++jt/dTnC3fFyyQNJarr5b3Jy8vn95eP3198vvMyLA/P7tWaMBVXWZbYgRqK3LzzdKzBeSLvM8+4uLbWEaMkB+Y5p//lMcKr2QwEMOHy49h992P5J133uHGG1diWe+x114HDPoxAgFYuvReVq48JHTslFMKgD8OqsyRcyxlZYoDDzyQt99+m7PPDtDcfCn77ruvqzlgLMIFKjMTfvSjNfT0zAtNjHrFG2uisyx7v8Hvz0N8+XKf1lYJQ47VguSAA+S93HVXuV5ZadcH3Hnnnent7WXBggV8++23dHV1udyWSskqV4/9pZfk/Vixwn78nJwcVz6Tjkbs65MAER2IWVMTaUGBXTtPV7CQKvJSaR/EjdVfgdDRo0ezZs0afD6YN68UmIBShcH3M/b95Pn3YdQoLx98IG/uIYesDtU3jEZqaippaTdy2mmy05+Q8BGFhQv7Tfp+5BF48UX5v7OzCN3rbP369aSnp8fs/9XaKt+V9PT0CAsK5P3RQSAdHXZD0MLCQtra2mhvb6eqqopdd63grLOit9rQlJS0AB9hWRbd3f0LlNcrKw6/v5OWlhaysrJCrv3eXnjqKYl01QKlVDtjx7YCHZxxxpsha3LhQqmrV1Ehr7WgIJ/e3m58viIeffR2Pv9ccqXAWfU/ENHkclthBGorEgjIBieIS6q3d+P3jUAEqqrKjsz55hvJg4q2+u0PXRl5n31Opbu7m+efn49S9ey//87939FBXBx0dIygpUXcQJYFixenAfEx3WHRuOwySbY89NBDaWxs4vrr72TVqlX8KFo4VBS090evAsvLweudCqiQS0qveGMJ1EknSZNEgJoaLYplwceVB44lUHFx7iCNSy6xm0fuGdyA+eijj/j4449dxzQi0PL/hg0yfuecqvscaZzlew4/HJ5+Wr5fN94oUWD248qHrEPbtQUVFyeWlnZF/uQn4nKOxeTJk/n2229D73N6+iR+/vNEqqr6jx6V70A31dUe1qwRl+3FF1dz5JGx7wNQXJxKY+NafD5YuLCA/Pz+GxPdeafdldjnSweaWbVqFRs2bGDkyJGuKNCyMlkovvOOuFWzsyMFSltQl1wijw12mDnYVvjixYvp6ekhK2sMK1faFnS0feX99luBZR1JV1cXl176GrBPTIHyeER4+vraQx2gnS4+jd7HLC/v4a9/DQAj2WWX5lDouLakenrk//z8fPr6uoF4urtT6OtzJ2hPm1YPvGsE6ofI3XfLHkddnb3i3VSBsix7xb1ypeS9DLYiukbnyqSl7cjpp58O7MCoUb5BWSxOMjNb6O4WMaqpgd5eD7A+pjssGrvsIqWNjj32WPLzd+TGGy8nOflifmw7xvslLU1W8jpIYu5cuPxyWaZrgRrIgioocAscQG/v6uB1OdBfhfdLLpF2CSBlrbR/v6ioiPHjx/Pf//6X1157jdLS0ojmi8OH28+5YYPsqzgLdxcVFdHa2kpnZyeWZYVcfPHx9rjj4qR7sLN5oH6eNWvW0NLSQktLSyj/qrDQjrisq+u/5cuOO+4YdJc1ExfXTUbGJLxeeYwoUfch5DtQG3xdPUAC8fHD+s2DAsjK2pWPPz6Gr76C+fNvwuvtvyqJM3S8uzsRpVpYuXKlq+ahxuORahgrV8qeYX5+bIFKS4vMgwL7e/Duu+8C8PXXBzJpki1Q0RaLTmtWu1zDBeqUU+DSS2GvvbzALmRllYXGM3KkLErPPFNcdpJbKUK2fr2f//u/fGAGK1eOR2c1OMPMPR79eQQAsUZXrnSXYLMskQQjUD9AdJTMggUwf76suKdO3fjHGTtWJiW9Ib98uYQqbyxTpkiRzcMPh8cee4zi4n3Ze++NyPYNkpvbgd9fQl9fn8P1uG6jLKjWVtnsr6tL5PrrHwXgqqvOiNnYLhqVlXbZqLIyKC62ACuUd1RdXU1aWlpo1RnO8OEShdfeLi7TqVMvp7Pza2BgCwrE5fL667J4eOYZd0rBaaedxjvvvMPLL7/MiSeeGJHXVVoqIhMI2CkDzlO0qFZVVdHW1obP5wtNbsXFct+6OmmK6dw+SE1NpaCggLVr17Jy5UqAUI6b04KqqCCieKyTHYOz2NdfL8SyKkhNHcftt9tWSyzkOyACVVnpA/Zin32m8MEH/d8vL28KlZWn8s47/uCY+w97zsx0FpdVZGT4Wbp0KcuXL2ecjnwJ4myPUlcn+4dOgXLuQelWHiDn6nqNOkXglVeko9DIkZn4/bIPee210X/XS5fOBtazfHkdzz47Gzgv4vt99NHws5/ByJFZwJd0dtaEBEopWYQqBW+9Ja5HLVB9fTqZycutt+4W+lyciboeD8EE/P8AsokdXt3mm29ygB2MQP0Q0QL11VcSIDFx4saFmGv22Ucmll12kclw+fJNs8QSE8Wtk5wMXV2Kysq4UNjpxjBsWA8wgsbGZoLGCunpjVFbI8SiulpcU2+9BUlJMwA49dR+fE5RcE7o69bB6NGSva8tqLKysn4tID2PrVghk9b06TVUVMikri2oWE0TQdqlr1snE9kpp0jlec0FF1zAmDFjKC4u5rIo/VWOPVasPp9PBCps0R9yKVVVVYVcldpCHTNGXMf/+pcsOsKrjIwZM4a1a9eGkqCnTJkCiPVTWyvPWVHRfx7dTkH/39/+9jcsaz6FhWn85S+y/9UfIlCdxMf3Bcc18L6VjFF+GG+9JbHjEydmDfA8dsDR/vvD+PFN/Oc//6Grq4sZdvw0IN/79HQRqBdfFEEYjAWVmGiXZxo+fDipqal8+umneDweRo+WqI+UFLFsHEGaIZKTs4ERrF9fz6JFk0lI2Dci6bupScLIq6qygXNYv74On89HZmYmPT1SoumNNyKj+Hp6dHHOeAIBFZEHNWuWWF5S0/J+4HbX7c8+63o3TZDED5HsbFlZvf66uPu0b3tjcU7E7e2w777uVhEbwzvvyDgSEqTKszO0ebBMmdINvM+6dQ3k5cGwYQspLOwe8H5Oxo6VCeDbb8UKSEqSCX9jePFFCVbo65Ow40mTJIdHW1Br1qzpt4STFvnlyyVLPy5uH8rLy/H7/ZSVlZGXlxdRBd3JqFFiNWn3rfOp8vLyWLZsGRs2bIhqhe26K5x/vnwO++8fWYnAKVDaVamPTZ4s6QuLFskEGj7577DDDixZsoRvv/2WhISE0Htw770y3qoqWUnrPcloFBQUMGvWLJ555hngJB5+OIX16/svcwTS4iEhIYHJk+eTkrKGxERRwYG8vzvtJK/t009l42W33fpvOZ2fL2JrWVJk+OyzbQtgjz32iDhfV5/fbz95/6R6fROWZbmCJLRANTeL+1TX5VNKhazKmTNnkp8vilBeLuOIlgdVWCiPWVnZRG9vPMnJ/ohzbrlFxvPRR/HA31m1Sj5rHeTxl7+I1yM8D6q3VwtUgkugjjpK3I7XXy/3HT16NKmpmSQmKs46y16UBQIwe7YehQmS+MFy/PGSe1FUJFnsm8of/yh9cNLTceVqbCyvvy6Vs3t7ReQGWYLPxQkn+IDDaG0t44ADYPz4iyku3riIDY9HLMpvvpG/iRP739uIRnOzCO5nn8lKdOJE+UGuXbsWy7JYs2ZNvyWcdthBKjFMniz5acuX74fP56OyspKVK1eywwDm5eTJMplpyynczSPll6K/qEBABGbVKvlsf/Ur9+1ajKqrq13deAHOOUes8iVL5DnDn2LnnXemqqqKF198kSlTpoSi4VJT5dykJHnOKPO4i18GS5IfccQRNDSU4vfTb8kikIl82LBhzJjxV3JyXiI5eTQJCQN7DmbOnARU0NWVCrSy775T+j3/mmvs/TSAE044gdLSUg4//PAIFx9IErTHI5ZDQ4MED9TV1dHR0UFvb29ov+h3v5PFSnW17C868xhPOklqYR999NGhsPm//EXEN7iGcDFsmC1QPl9iqCGik6QkccfpaMCyMgnlzMzMDO0x9/bazQ21kHZ06JpN4rXQ3wE5x77u8XjYa69yRo36lkceIZRAfcopdoSxEagfMBdfLH/h1Y43lrY2yQyP9kPYGObMkS/86afLHsqmoN1mFRWV9PZKbbKNCZDQzJ4tVtyZZ25cl2GNTi1at072s448EiZNmsSKFSuoqqqitbW1X4FKSpLcoOHDxc02ebIVfLx1rFixot+WH3r8Y8aI6Gdk9G+RRGP33eHPf47ePj0vLw+v1xvVgho5Ulx7CxdGF4zdg6GJK1euZD9HBEV1tbi3Vq2STfeB9kPPOussFixYwKOPPhtqbBgWjBiVgoICampqWLeugsTEsZSWDty3bNy4cSQmfgssZ4cdriUtrZ/aSEjEY2amdI/Nzob58/NYvXp1aI8onH//W9yqP/6xuNTy8/NpaWkJib+uPDF+vIiZ3t9yRlZedNFFvPPOO/zmN79h8mSJQtUJ29EEOD1dFgZVVR1YVjxZWQkR52gDXT/fihWyB1pQUEBcnFhNvb0igoWFhHqb+XwbuPzyvwEvAvbe07Jl0m7jiCPsaikZGWl4PAnBMdl7j3Y/MyNQP1hycmRV9l05+GDZOygullYDm8o++4gb44UXZKW4Kci+zOv86U97kZFhsX79cQO2i4/G4YdLsMTkySKYG8v06WIVfP65JK+OGSNVKHw+H8899xxgb27HoqdHRArgwANlw+Gbb76hoqJiwPvutJO9wt5778E3jwSJwJs2DR54QD6P8ATuuLg4hg0bRmVlZbC0TqIrt2fZMimxEy0qf/bs2UyfPh2Px8OZZ54ZOh4fL8/3j3/YwRL9oZRixowZ5OWlcsMNYi0M5mMeMWIE8+YdzJdfvsnkyUu46qrBPdcvfvEuMJFLLhl4g3XtWmkz8vHHMrmnpUnDxP6KDOv90tGj7UoTuuKHFqjly+1eU+AWKI/Hw3777Ud8fDxFReKiTU2VzzJaHlRpKWRmvklDQxVQR15eZF6XFqimJlCqj7a2ZsAOkklIkO/oRx/Z/ZuknFUdXu9avN5W3n3XCnlUNmyQqjP6tYKMb8kS+fw/+cSu7H7//foMn9mDMnw39t7bjtzbZ59Nf5zUVGm0d8wx4ibZFJKTk0lIqGDlyrH09Ch6ej6LCO0dDEccIRbUIApHRCUxUWrYvfmm/YPUZZIeeOABwN7sj8Wbb8qeXHo6HHdcMSkpKTwWDIkKr/4Qi7fftkVuY9DiMm5c9JDvUaNGsWbNmlCwR/jke9990QVKKcU777zDN998E9o3AYlmGz1aJqaB8pLCufhiu7TQQOywww7U1y/FspKZNq2Pc88d3P1uv/2PfPPNN1xwwQUDnltfL9an9gIMZMC/8YakBYC8D1qQtEDlBU2hTz+VBaAO245VASMQkMjc+fPdLeWdTJsGBxxwP6tWPQkUcMQRayPOcQqU12uLhPZIpKdHPnZubi61tS28884ssrL2Z84cFQqy0ZaU5FbJ/7rskTaSnCH/U6Z0AguNBWX4bsTHy37LO+9s+v6T5sc/FpdHPwFqA1Ja+hYQoLS0B/hwkwQqMdGuxLCp3HqrREtqd8Xo0aMpKSlh6dKljB49ekDL7ogjpAHfBx9AcrKHmTNn8vnnnwOw2yAjUXJznYU3B88vfylVsINaGsHEiRNZvnw5K1eujNhXmTZN7hvLYMjJyYlaof6II+RyYwVqY5Cxyl5Kefl+DHZxHhcXx5QpUwZstQK2IC1c6L4eC72deMgh8p5pgVoabKSkr+s0pfJymeBjFKTAssTF+7//RS+4q5k4cWKwPQdRA3b23lsWGjfdBD/96d8AXD3Iqqqk+sxhh0ntSdBdk1v58suT8Hh+xJNP2pa83rfSYebgbkdfVmYHWHk8Jg/KsBkpKZEopI1xJW0pxo1rYMqUU/nDH94C/KEK2Fub8eNl41tHJCml+NnPfgbAuYNYuislVRW0wBwRnMGnT5++SW7LjSErS8Qp1r7OxIkTqa2tZf78+QO6GwfLH/8oFS90i/AtgdRSlPC3F16YiaP92GajpETcX6tXi6jESHULMWaMROQ9+aRc1y6+b76RPklaoHTO1NFHS+RcLOHzeMTdmZsb+72sr4d77rkOuAX4F5YVWQR5yhS44AKxbGfPlhcRHx9PnENVLEvC+3WeXW5uLo2NtcFzh/OTn9g1FsPzoEAi+047LXJ8fj8sWZIEDB8yAhW7uJXBsBGMGzeOTz55jIYGKTS7KRbUluKaa67hmGOOYeomZEVfdNFFNDU1ccp3NVM3A87xh+f2bCoZGdKSYUuy0047MX78cFpbb+bII3/Dzjtv/hWVJKFKekFwPTIgzsLN2oW2cOFCsrKyQtFxWqAaGgZeCI4cKS7zWPunErKeCEwBjiA1NbKcRlubiOy330J1tVRC2dXhVvjNb2z3r04zzM3NpalJCgbGx2cGL+U2LVAHHGAHwZx0kliQTz5pi9bDDzu/BxlGoL4rSqkc4FlgFLAOOMmyrKYo570G7AZ8ZFnW4Aq7GTaaKVOm0NbWxmuvvUZ+fv5GVTLf0sTFxTFt2rRNum9aWhq33nrrZh7RprHvvvuG/t9///234Ug2Do/Hw4IFC1BKkZy85cz9iRMlfUNXe9gYsrKyQlX5nWWotIvvmmsksEB3qo1GUZG4yjs7o0fpJiaKcIwatS8rV0JaWuR78fHH4r7LyYHRo0fx3HPPueo2vvKK7YrX+1V5eXk0NNQCAbxe6S+lBWr2bLGMnG69jg4RposustNKAgFJ/Jc8L1PNfHPwW+Bty7LGAW8Hr0fjz8AmxIQZNgZdneDNN99kxx13HNS+gWHjSE5O5q233uK5557bZi7UTSUlJYVkXYZhC/HQQ1KJZFNQSoXcpk6Bys2V/LSCgoHTMHRe46efxnoOCT9vaBDXXbRQdC06jY2QkqI48cQTXdVLEhPtmoPagho5cmSwIaUPj0csPy1QSok4dXfLH0hwyOGHwz332AFJ557rzIMzYeabg6MBXQnsMeCYaCcF2863baUx/WBxNhXce++9t+FIvt8ccMABnHjiidt6GEOS7Ozvth+rIzx3dvj+PB4JQGlq6r9qO0iY+bJl4k6LRWam3ZgyWkSgszpYtIKzCQkiNNOm2cEYtqCO4NRTpSSXFqiaGnF5jhljByDFxUkEX3e3VB7R3njdkdgI1OZhmGVZOk21Gtj4zFAHSqnzlVLzlFLz6jali+APnNTUVC677DKGDRvGGWecsa2HYzBsNFdeeSWXXnppKKhG8+KLUp9xoLgUpQauiXnKKeJ2GzEiukA5K2nFEqjsbLHq9DrFrnBSw8EH5zBvnpQ/A7G2HnhAov+cUXy1tSJOX35pF/zV7eyhzwjUYFBKvaWU+ibK39HO8yxpmTpAAf/+sSzrAcuyZluWNTu/v74DhpjcfvvtVFVVRbSSMBi2B0aOHMkdd9xBdphyzJsnl7pX2HfhxhtFFNavtwvPOnEeiyZQ0hrEfWzUqFGMHTuW5OTfUVU1mZ12ssPhnT0eo4WZx8XZeVAFBTB+fABYN2QEakgHSViWFbPvplKqRilVZFlWlVKqCF3T37BNMXtPhu8b114rFkl/XYA3hq6u6OIEUm3ikUck4i5aVa5//UuCNXbbTSpJ7LGH/OY+/PBDpk/P58EHPdTUSFh8YaG7R1w0gaqvt0POPR4IBOT3a4IkvjsvA7puy5nAS9twLAaD4XtKUhIcdNDmyTf87W8lwu/SS6PfnpYGZ50lbsCwXoYhmpqknJcuKAtSlzEx0cvChVL5Ys0aOR7NgjrySKkcA26xqqqCVasUUDBkLKjtWaBuBg5SSq0EDgxeRyk1Wyn1d32SUupD4J/AAUqpcqXUd6ghbjAYDJuODmzQbsNovPWWlKxatCjytr/+VSq/gDugAsRa6gz2LdTt4aUlhwRI/PSncuzgg+1q+VrA7r3XfhylkoeMQA1pF19/WJbVAETEy1iWNQ8413HdhJQZDIYhgQ606C+P/aCD5PLss93t2EHce8HemRHVMrxecR+C7drLybEbLmqamiTv6eqr7Z5rgYCI2Oefy+MYgTIYDIYfGIccItXFg62k+iVabUynKIULVHy8lDTS/ztpaLBrCd58s+RAdXWJWIFYVOecYwuU2YMyGAyGHxjx8dLrLEpT5RB33CGRetGKn+gIvp12ikz0/fJLSVbWzwNSReK006RP1cknyzGduFtfL6KlgzG0peX1WkPGgjICZTAYDEOISy6RZN5oPaW01fTGG5GFa1NSJDdq+XLbhRgXB089Jf+HR/Hl50vtwrlz5fqzz8plfLxx8RkMBoMhCkq5o++cDBsm+1jalefknnvkfs72WUqJMPn9tlXljNxTyu4LNW2auAL9/h4jUAaDwWDYOE4+WZpxHn20Luxq8/zz0g+rsxMuvNDOtfJ6YwtUV5e03wCdBwVer3fICJRx8RkMBsN2xMqVUFEReTw5WUobXXGF28LS1pgWqMMPt5OOtdsPRNyqq8Hj8ZogCYPBYDBsHGVlEghRWRl5m7OOnzOKr6REIgJPPVWu77qr7QbU4nXzzXYfq4SEeGNBGQwGg2Hj0PX5orWVd5ZPcv6/fLlYXEceKdfr6sTl98c/2qHsgYCUWOrqGlouPrMHZTAYDNsJmZnw+99DtI42WpQSEtz7TCCljzIyJNz8oYfgqqtEjPR5V18tAjZ16tASKGNBGQwGw3bEddfBfvtFHn/4YUm4dbr6QBJwx46F66+X61qUVq0Sy2ncOLnu98ul12v2oAwGg8GwGVEK/vAHaZro5J135FLvS+nAiGnTpO3HTTfZ9wdjQRkMBoNhM/PSS3DRRZGdf7Uw6YCI8H5QutyRFqj4eBMkYTAYDIbNyOLF8OijcPfd7uPhYebhAnXeefK/saAMBoPBsEXQe0/33OM+Hi5QBx1k72E5xUrX9jMCZTAYDIbNio7i0y3cNaNHy+WPfiSXkyfbOU9xcWI5XXstDB8ux4ZSkIQJMzcYDIbvAVqgwlttvBTWa7ymRqL3br9d+kXFxblFzexBGQwGg2Gzol104WHmPh989ZUk6III1rnnSmfe9HQ7aXflSrnduPgMBoPBsFk57jhx0+2yi/v4z38Os2fDM8/Idb3vtHAh9PWJyw9sK2ooCZRx8RkMBsP3AKWkYkRHh/v4/PlymZIil1qgjjhCisNedZXsSTmj+IbKHpSxoAwGg+F7gtcb2WlXC5LeowoPMy8rcx8fShaUESiDwWD4HqMtIy1QzhYbcXFSh895ngmSMBgMBsNWIdyC2mcfu9is05pKT5dLY0EZDAaDYatw8MFyOX26XI4cCcceK//HxYkwXXIJ5OfLMSNQBoPBYNgq3HijROjpGn3V1XI5d65YVUq586BMkITBYDAYtgkffgiXXQb77iu9o1pbpX6fbiNv9qAMBoPBsE3Q+04ffigVJbTrTwdPGBefwWAwGLYJWoguuEASdc8/X66bMHODwWAwbFPC86CWLnUfN3tQBoPBYNgmhAvU3Lnu416vF8uyCAQCW39wYRiBMhgMhh8Qu+4Ke+wh/zvFSheZjQ+WQx8Kbj4jUAaDwfADIj8fDjxQ/lcKCgulq66u1ecNdjg0AmUwGAyGrUptrQRK3H+/XO/rg7Y2dzVzMAJlMBgMhq3M4sXwu9/ZbTYaGqQVR3e3XNcCNRQCJbZbgVJK5Sil3lRKrQxeZkc5Z4ZS6lOl1LdKqUVKqZO3xVgNBoNhqKD3nV57TS533NF9XO9BGYH6bvwWeNuyrHHA28Hr4XQCZ1iWNQU4FLhLKZW19YZoMBgMQwudB/WnP8nlSSfJpRaoxMREAHp6erbyyCLZngXqaOCx4P+PAceEn2BZ1grLslYG/68EaoH8rTVAg8FgGGrEhc36uqGhPp4UDOfr1j6/bcj2LFDDLMuqCv5fDQzr72Sl1C5AArA6xu3nK6XmKaXm1dXVbd6RGgwGwxAhXKBeeMF9XAuUsaAGQCn1llLqmyh/RzvPsyzLAqwYD4NSqgh4HPipZVlRs88sy3rAsqzZlmXNzs83RpbBYPh+Mm0azJ5tNyjU6OtDyYLybusB9IdlWQfGuk0pVaOUKrIsqyooQLUxzssAXgWusSzrsy00VIPBYNguSE+HnXeG9evl+tixkryrGUoCNaQtqAF4GTgz+P+ZwEvhJyilEoB/A/+wLOv5rTg2g8FgGJI0NIhI3XabXG9vh6oq+3YdJGEE6rtxM3CQUmolcGDwOkqp2UqpvwfPOQnYBzhLKbUw+Ddjm4zWYDAYhgDl5XDrrXaL95oaePdd+/ahtAc1pF18/WFZVgNwQJTj84Bzg/8/ATyxlYdmMBgMQ5aEBLl8+WVp/T51KnzzjX27cfEZDAaDYZugBerRR+XykEMgNdW+3QiUwWAwGLYJWqA0X3wBHR32dbMHZTAYDIZtQrhAffih+7rZgzIYDAbDNiE3F0aPBr8/+u1DycVnBMpgMBh+QHi9MGmStN0A6Q3V2WnfbgTKYDAYDNsEy4JRo+CYY+R6XR1UVtq3e71ePB7PkHDxmT0og8Fg+IHx179CWZn8v2KFiJSTpKQkurq6tv7AwjACZTAYDD8glBIr6n//k+vr1sGGDe5zMjMzaW5u3tpDi8C4+AwGg+EHyLx5cllQEHlbTk4OjY2NW3dAUTAWlMFgMPwAOfHE2LcNFYEyFpTBYDD8wOjttTvrRqOxsZFvvvmGlpYWMjMzt97AwjAWlMFgMPzAiI+PbFzo5IADpMzptg41NwJlMBgMBhd//vOfqaysZNiwfhuVb3GMQBkMBoPBRXx8PEVFRdt6GEagDAaDwTA0MQJlMBgMhiGJESiDwWAwDEmMQBkMBoNhSGIEymAwGAxDEiNQBoPBYBiSGIEyGAwGw5DECJTBYDAYhiTKsqxtPYYhh1KqDlgfdjgPqN8GwxkMZmybhhnbpmHGtmmYscVmpGVZ+eEHjUANEqXUPMuyZm/rcUTDjG3TMGPbNMzYNg0zto3HuPgMBoPBMCQxAmUwGAyGIYkRqMHzwLYeQD+YsW0aZmybhhnbpmHGtpGYPSiDwWAwDEmMBWUwGAyGIYkRKIPBYDAMSYxAAUqpPyullimlFiml/q2UynLcdpVSapVSarlS6hDH8UODx1YppX7rOD5aKfV58PizSqmE7zi2E5VS3yqlAkqp2Y7jpymlFjr+AkqpGcHb3guOTd9WEDyeGBzTquAYR22hsY1SSnU5nv9+x207KaUWB8dwj1JKBY/nKKXeVEqtDF5mb6GxHaT+v717D7GjPOM4/n2a1VSKWpN6iUYxC2qrFbSKFUGQeJdqjEqMFNPGP6RK/2spgW1roCAo+IcgGBDiJdHEew0twdY7tKaVxNQLmpqNUrPGWEPrFVOiT/94n01eT87ZPXN25szs5veB4cx537k8+8yceee8Z3bGbH3EsN7M5mZ1teYt6mrd39rE+mCWj3fMbGOUF97GZTOzpWY2ksVwSVZXKI8VxNb2mNKEvLWJtS856Ym77/MDcAEwEOO3ALfE+InAP4DpwBxgGJgWwzAwCOwf05wY8zwELIzxZcANE4zte8AJwHPA6R2mORkYzt63nRa4EVgW4wuBB6uIDTgWeK3DPH8HzgQMWAtcHOW3AktifMnoNqggtlOBI2P8+8BIg/JW+/42Tty3Ab/tdRtXEM9S4JdtygvnsYLYOh1Tas9byzr7lpNeBn2DAtz9T+6+K96uA2bH+DxgtbvvdPe3gc3AGTFsdvct7v4/YDUwL8545gKPxPz3ApdPMLY33H3TOJNdEzGMZ17EBCnGcydyltZlbLuZ2SzgIHdf5+nTcR978pPHVlne3P1ld38v3r4OHGBm08dZXL/yVvv+1kmsawGwapzpxtrG/VIoj1UEMMYxpa0a89a3nPRCDdTeriOdvQAcBbyb1W2Nsk7lM4H/ZjvmaHnVrmbvA8fd0YXwm+xgujvuiPGjiLkKc8zsZTN73szOzta/NZsmz8/h7r4txt8HDq8ortyVwAZ335mV1Zm3Ju9vZwPb3f2trKzoNq7Cz6MbbXnWLVw0j1XLjynQjLyNqisnXRmoO4B+MbOngCPaVA25+xMxzRCwC7i/abGNMe8Pgc/d/bWs+MfuPmJmBwKPAteSzsj6Fds24Bh332FmpwG/N7OTul2nu7uZjfv/DxPM20mkrpcLsuK681aLLmO9hq+fBE1oG5cRG3An8DvA4/U2UmPQFz0eU/qSt6lin2mg3P28serN7KfAj4Bz4ys2wAhwdDbZ7CijQ/kO4NtmNhBntfn0Pcc2joW0fHty95F4/cTMHiB9jb+PPX/PVjMbAA6OmEuNLb6R7Izx9WY2DBwf68+7OvL8bDezWe6+Lbo7PuhiPT3lzcxmA48Di9x9OFterXmjT/tb0Vjjb74COC2bp5dtXFi3eTSzu4A/xNuieawktnbHlH7lrYCxclU7dfGRrmIBfgVc5u6fZ1VrgIWWruKaAxxH+iHzJeA4S1dQ7U9qJNbETvgscFXM/xOgsrNlM/sG6XeB1VnZgJl9J8b3I31ARr9drYmYiBifyRrjMuM61MymxfggKW9bogvvYzM7M7rPFrEnP3lsleUtrqb6I+mCjL9k5bXnjebub+cBb7r77i6oHrdxqeJEZtR8vr69us5jRbG1PaY0IW8t+paTnvTzioymDqQfUd8FNsawLKsbIl3lsonsqhrgEuCfUTeUlQ+SPgybgYeB6ROMbT6pX3gnsB14Mqs7B1jXMv23gPXAK6SLAG4HpkXdNyOmzRHjYBWxkX7beT1yuQG4NJvndNKBZBi4gz13M5kJPA28BTwFzKgotl8Dn2XbeiNwWBPy1oT9rUO89wA/aykrvI0riGsF8GpsszXArF7zWEFsbY8pTchbm1j7kpNeBt3qSEREGkldfCIi0khqoEREpJHUQImISCOpgRIRkUZSAyUiIo2kBkpkErF0x3WPYdf4cxRa9hfZsp8rc9kivVADJVITMzvF0iMjji0464ek2zAtalne0qyBuardjGY2P5tmaUv14ljuhwXjEamEGiiR+pwC3ER6BEMRn7n7Snd/oEP9F6TGpp3ron4v7r7K3VeS/pFZpHZqoESmnseBC1tuBYSZHQFcBDxWS1QiBamBEqlBdK/dHW+fzbrd7ilh8SuBr2jpAoz3HvUijbfP3M1cpGEeA2YB1wM3A29E+XDHObr3AemGuItJjxQZtZh0x+9/l7AOkcrpG5RIDdz9FeDFePvn+E1ppbu/ONZ8BSwHTjCzswDi9btRLjIpqIESmZrWkp5MPHqxxGLSw/LWdpxDpGHUQIlMQZ4eYLgCWGBmhwJXAyvc/ct6IxPpnhookalrOXAQ6XHjB6LuPZlk1ECJ1KfSh7G5+5uk37nOB/7q7puqXJ9I2XQVn0h9Po3XGRWuYwkwF3imwnWIVEINlEh9XiL9v9KQmR1CuoPD2+7+t7JW4O4vAC+UtTyRflIXn0hN3P1fpFsPHQDcCawCbqg1KJEGMfdKu8FFpERxl/FB4AeAu/uOEpc9EzBgA7DF3c8pa9kivVAXn8jkczTpbhBfUu5neASYHuNbSlyuSE/UQIlMLr8ADonxr0pe9sXAtBj/T8nLFilMXXwiItJIukhCREQaSQ2UiIg0khooERFpJDVQIiLSSGqgRESkkdRAiYhII/0fixW5pxSyZZQAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "plt.plot(strain.t-1700, strain.y.real, 'k', label='ETK')\n", "plt.plot(t, h[(2,2)].real, 'b--', label='NRSur7dq4')\n", "plt.ylabel('$h_{+}^{(2,2)}$', fontsize=18)\n", "plt.xlabel('t [M]', fontsize=18)\n", "# plt.xlim(-100,100)\n", "plt.legend()\n", "plt.tight_layout()\n", "plt.show()\n", "# plt.savefig('/Users/liuyu/Downloads/GW_BBH002.png')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }