{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "\n", "\n", "# Standalone Fishbone-Moncrief C Code\n", "\n", "We start with the NRPy+ expressions generated in the [Tutorial-FishboneMoncriefID](Tutorial-FishboneMoncriefID.ipynb), and output them to the C file \"FishboneMoncriefID/FMstandalone.h\".\n", "\n", "Further, $\\Gamma = \\alpha u^0$ is given by (as shown [here](Tutorial-u0_smallb_Poynting-Cartesian.ipynb)):\n", "$$\n", "\\Gamma = \\alpha u^0 = \\sqrt{\\frac{1}{1 - \\gamma_{ij}v^i_{(n)}v^j_{(n)}}}.\n", "$$" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2020-12-17T02:14:18.841310Z", "iopub.status.busy": "2020-12-17T02:14:18.837158Z", "iopub.status.idle": "2020-12-17T02:14:24.365260Z", "shell.execute_reply": "2020-12-17T02:14:24.364668Z" }, "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Wrote to file \"FishboneMoncriefID/FM_standalone.h\"\n" ] } ], "source": [ "import sympy as sp # SymPy: The Python computer algebra package upon which NRPy+ depends\n", "from outputC import lhrh # NRPy+: Core C code output module\n", "import indexedexp as ixp # NRPy+: Symbolic indexed expression (e.g., tensors, vectors, etc.) support\n", "import finite_difference as fin # NRPy+: Finite difference C code generation module\n", "import FishboneMoncriefID.FishboneMoncriefID as fmid\n", "\n", "# Step 1: Set up the Fishbone-Moncrief initial data. This sets all the ID gridfunctions.\n", "fmid.FishboneMoncriefID(\"Spherical\")\n", "\n", "gammaDD = ixp.zerorank2()\n", "\n", "DIM = 3\n", "for i in range(DIM):\n", " for j in range(DIM):\n", " if i<=j:\n", " gammaDD[i][j] = fmid.IDgammaDD[i][j]\n", " else:\n", " gammaDD[i][j] = fmid.IDgammaDD[j][i]\n", "\n", "# gamma_{ij} v^i_{(n)} v^j_{(n)}\n", "Gammacontraction = sp.sympify(0)\n", "for i in range(DIM):\n", " for j in range(DIM):\n", " Gammacontraction += gammaDD[i][j] * fmid.IDValencia3velocityU[i] * fmid.IDValencia3velocityU[j]\n", "\n", "Gammafactor = sp.sqrt(1 / (1 - Gammacontraction))\n", "\n", "# -={ F-M quantities: Generate C code from expressions and output to file }=-\n", "FishboneMoncrief_to_print = [\\\n", " lhrh(lhs=\"alpha\",rhs=fmid.IDalpha),\\\n", " lhrh(lhs=\"betaU0\",rhs=fmid.IDbetaU[0]),\\\n", " lhrh(lhs=\"betaU1\",rhs=fmid.IDbetaU[1]),\\\n", " lhrh(lhs=\"betaU2\",rhs=fmid.IDbetaU[2]),\\\n", " lhrh(lhs=\"Gammafactor\",rhs=Gammafactor),\\\n", " lhrh(lhs=\"Gamma_times_ValenciavU0\",rhs=Gammafactor*fmid.IDValencia3velocityU[0]),\\\n", " lhrh(lhs=\"Gamma_times_ValenciavU1\",rhs=Gammafactor*fmid.IDValencia3velocityU[1]),\\\n", " lhrh(lhs=\"Gamma_times_ValenciavU2\",rhs=Gammafactor*fmid.IDValencia3velocityU[2]),\\\n", " lhrh(lhs=\"uKS4U1\",rhs=fmid.uKS4U[1]),\\\n", " lhrh(lhs=\"uKS4U2\",rhs=fmid.uKS4U[2]),\\\n", " lhrh(lhs=\"uKS4U3\",rhs=fmid.uKS4U[3]),\\\n", " lhrh(lhs=\"uBL4U1\",rhs=fmid.uBL4U[1]),\\\n", " lhrh(lhs=\"uBL4U2\",rhs=fmid.uBL4U[2]),\\\n", " lhrh(lhs=\"uBL4U3\",rhs=fmid.uBL4U[3])\n", " ]\n", "# print(fmid.uKS4U[3])\n", "fin.FD_outputC(\"FishboneMoncriefID/FM_standalone.h\",FishboneMoncrief_to_print,params=\"outCverbose=False,CSE_enable=False\")" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2020-12-17T02:14:24.369070Z", "iopub.status.busy": "2020-12-17T02:14:24.368440Z", "iopub.status.idle": "2020-12-17T02:14:24.370477Z", "shell.execute_reply": "2020-12-17T02:14:24.370973Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Writing FishboneMoncriefID/FM_standalone.c\n" ] } ], "source": [ "%%writefile FishboneMoncriefID/FM_standalone.c\n", "\n", "#include \"stdio.h\"\n", "#include \"stdlib.h\"\n", "#include \"math.h\"\n", "\n", "const double a = 0.9375;\n", "const double M = 1.0;\n", "const double r_at_max_density = 12.0;\n", "const double r_in = 6.0;\n", "\n", "int main(int argc, const char *argv[]) {\n", "\n", " // Step 0a: Read command-line input, error out if nonconformant\n", " double xx0,xx1,xx2;\n", "/*\n", " if(argc != 4) {\n", " printf(\"Error: Expected three command-line arguments: ./FM_standalone r theta phi\\n\");\n", " exit(1);\n", " }\n", " xx0 = strtod(argv[1],NULL);\n", " xx1 = strtod(argv[2],NULL);\n", " xx2 = strtod(argv[3],NULL);\n", "*/\n", "\n", "// printf(\"# Output: r,th,ph, alpha, betaU0, betaU1, betaU2, Gamma, Gamma*vValenciaU0, Gamma*vValenciaU1, Gamma*vValenciaU2\\n\");\n", " for(double xx0=1.6;xx0<50.0;xx0+=0.2) {\n", " xx1 = 1.56463634120e0; //M_PI/2.0;\n", " xx2 = 0.0;\n", " double alpha,betaU0,betaU1,betaU2,Gammafactor,Gamma_times_ValenciavU0,Gamma_times_ValenciavU1,Gamma_times_ValenciavU2;\n", " double uKS4U1,uKS4U2,uKS4U3,uBL4U1,uBL4U2,uBL4U3;\n", "#include \"FM_standalone.h\"\n", " if(xx0 < r_in) {\n", " Gammafactor = 1.0;\n", " Gamma_times_ValenciavU0 = Gamma_times_ValenciavU1 = Gamma_times_ValenciavU2 = 0.0;\n", " uKS4U1 = uKS4U2 = uKS4U3 = 0.0;\n", " uBL4U1 = uBL4U2 = uBL4U3 = 0.0;\n", " }\n", " printf(\"%e %e %e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e %.15e\\n\",\n", " xx0,xx1,xx2,\n", " alpha,betaU0,betaU1,betaU2,\n", " Gammafactor,\n", " Gamma_times_ValenciavU0, // util1(1) in FMtorus.f90; util(1,i,j,k) near the write statement\n", " Gamma_times_ValenciavU1, // util1(3) in FMtorus.f90.\n", " Gamma_times_ValenciavU2, // util1(2) in FMtorus.f90.\n", " uKS4U1,uKS4U2,uKS4U3,\n", " uBL4U1,uBL4U2,uBL4U3);\n", " }\n", " return 0;\n", "}" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2020-12-17T02:14:24.376108Z", "iopub.status.busy": "2020-12-17T02:14:24.375824Z", "iopub.status.idle": "2020-12-17T02:14:24.733826Z", "shell.execute_reply": "2020-12-17T02:14:24.734257Z" } }, "outputs": [], "source": [ "!gcc -O2 FishboneMoncriefID/FM_standalone.c -o FM_standalone -lm" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2020-12-17T02:14:24.739653Z", "iopub.status.busy": "2020-12-17T02:14:24.737603Z", "iopub.status.idle": "2020-12-17T02:14:24.849288Z", "shell.execute_reply": "2020-12-17T02:14:24.849011Z" } }, "outputs": [], "source": [ "!./FM_standalone > out.txt" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2020-12-17T02:14:24.864880Z", "iopub.status.busy": "2020-12-17T02:14:24.861157Z", "iopub.status.idle": "2020-12-17T02:14:25.608280Z", "shell.execute_reply": "2020-12-17T02:14:25.607793Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "--2020-12-16 21:14:25-- http://astro.phys.wvu.edu/zetienne/torus_cuts.csv\n", "Resolving astro.phys.wvu.edu (astro.phys.wvu.edu)... 157.182.3.45\n", "Connecting to astro.phys.wvu.edu (astro.phys.wvu.edu)|157.182.3.45|:80... connected.\n", "HTTP request sent, awaiting response... 200 OK\n", "Length: 16445 (16K) [text/csv]\n", "Saving to: ‘torus_cuts.csv’\n", "\n", "torus_cuts.csv 100%[===================>] 16.06K --.-KB/s in 0.01s \n", "\n", "2020-12-16 21:14:25 (1.40 MB/s) - ‘torus_cuts.csv’ saved [16445/16445]\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAEWCAYAAAC9qEq5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABPVklEQVR4nO3dd3xV9fnA8c+TnZCEQAhhB2QP2SBWUMEFto66N/ZnReuqda+2bm3raJ1VW0WsojhQFBBRZKgoMkQQlI2AbAgJZOc+vz/OueEmZFzgruQ+79frvpJ7vmc85+bmPOf7Pd/zPaKqGGOMMaEQE+4AjDHGRA9LOsYYY0LGko4xxpiQsaRjjDEmZCzpGGOMCRlLOsYYY0LGkk6UEZF7ReR/Ydz+DyJyfC3lM0Xk96GLKHKI4xUR2S0i88IdjwkvEbleRH4RkcUBXOdiEdksIn8K1DoPVr1MOiKyTkS2iUgjn2m/F5GZPu9VRPaJyF4R2SQiT4hIrE/5TBEpcst3iMh7ItLSj23fKyL3HkLMvvHsFZFcd/rxbtnEKvP3cafPrG599ZWq9lTVmXD4CdD97Dw+n+leEfnQZ90qIn+ssswf3en3+qxj46HvUUANBU4C2qjqYKj4Xq9y9+1jEWnlndlNUn8TkZ3u628iIv5sKFL22/1fPvEQlxURuU5EvheRAhHZ4v5fXxDoOMPkXuAaVe3jnVD1f0ZEWovIjyLylPt59BSRT0Rkl4jkisgCETnVO7+7riuBvx5OYCLSWUTeFJHtIpInIitF5GkRaVPXsvUy6bhigT/WMU8fVU0FjgPOB/6vSvl1bnkXIAN4MtBBVheP+8rwmb4dOFpEMn2mjQZWBDmehuAXn880VVVP8ylbAVxWZf6Afq4iEheodQE5wDpV3eeu+3jgYeAMoCmwFhjvM/8Y4EygD9AbOA24KoDx1Mn3RC4MngJuBG4GMoHWwD3AyDDGFEhNgaU1FYpIDjAbmKSqN6hzp/+HwHSgBdAcuAHIq7LoUqDxof7tRKQT8A3wC9BPVdOBY4DVOCdOtarPSecfwC0iklHXjKq6CvgS6FtD+S7gXaCXiAwSka1VakVnVVfFFZFmIvKRe0axS0TmiMihfKYlwPvABe56Y3GS5Os1LSAiSSLyP/cMN1dEvhWRbLeslYhMcmNaJSJX1rCOqSJyXZVpi0XkLPf3biIy3V3PTyJyXg3rGS4iS3zeTxeRb33ezxGRM93f14nIiSIyErgLON89i/f9fHNE5EsRyXfP2prV+MnV7lsgRUR6utvuCSS50xGnpjwVaOVTU2olIoki8k9xmjZ+cX9PdJc5XkQ2isjtIrIFeEVELheRL6p8Jur+cyIip4rIMnd/NonILdV8hlcA/8E5+dgrIvcBvwHeVtUfVLUEeAA4VkQ6uouNBh5X1Y2qugl4HLj8ED8r31i6uzWGXHGaQ0/3KRsrIs+LyBQR2QcMdz+zd92z3rUicoPP/PeKyAQRGefu/w8iMtAtew1oB3zo7vNtIvKMVK65lkk1LQsi0gW4BrhAVaeraqGqlqvqF6p6uc98vxOR5e6214jIVT5l3r/lbeK0nGwWkTPdv9cK93t/V5V9edv9v8sXkSUi0kVE7nSX3yAiJ/uzbT/+Bt7jj6eG8o44Ced1Vb3NndYM6AC8pKol7utLVf2iyuLedR7qCdO9wJeqepOqbgRQ1W2q+k9VfbOuhetz0pkPzAQO+AeuSkS6AcOAVTWUNwPOBhap6rfATuBkn1kuBcYBqOq9qnqvO/1mYCOQBWTjHEQPdVyhcew/Kz8F52zkl1rmHw00BtrinOVdDRS6ZW+6cbUCzgEeFpER1axjPHCh942I9MA5257sHpCnA2/gnDFdADznzlPV10BncZJwPM5ZdysRSRORZGAgMMd3AVX9GOcs/i23htLHp/gi4HfudhPw429ci9fY/7mOdt97Y9gHjKJybekX4G5gCM5JSh9gMM4ZtFcLnLPQHJzaRl3+C1ylqmlAL2BG1RlU9b84f8O5bhze5g/f5jLv773cnz0B32S92J3mzOw0O13kR3z7N+D8/T4EPsH5/K8HXheRrj6zXQQ8BKQBX7nzL8apaZwA3Cgip/jMfzrOdzIDmAQ84+7zpcDPwGnuPv9dVa/z/i1wzpp3Ax9UE+oIYIOqzq9jl7bhJO90nO/UkyLS36e8Bc6JSGvgL8BLwCXAAJxjxp9FpIPP/KfhfIeaAIuAaTjH0dbA/cALB7Ht2pwEFANbqyk7AifhvKCqf/GZvhPnGPc/N3lm17Du7UCRu40KInKRe6JR06udO+uJOCfph6Q+Jx1wviTXi0hWDeUL3bOx5TgJ6rkq5U+Jc21lMbAZuMmd/irOFw8RaYqTBN6oZv2lQEsgR1VLVXWO1j6Y3UKfP+BTvgWq+hXQ1P3nvgw3ydWiFCfZdHLP8Baoap6ItMWp6t6uqkWq+h3OGXTVZiaAiUBfcarpABcD76lqMc4/yzpVfUVVy1R1Ec4X7dyqK1HVQpzaw7E4/6yLcWqWx+AcvFeq6s469sfXK6q6wl3vBGqoobpaVfnHqFob+x9woXswvcB9X5eLgfvds7ftwH04Jx5eHuCvqlrsxliXUqCHiKSr6m5VXejHMgAfA+eJSG83ef8F56QmxS1PBfb4zL8HSBVxruuoam9Vre57W5sh7nofdc+UZwAf4XNyAnzgnkF7gCOBLFW9351/Dc6B2/e6yheqOkVVy3EO2L4nGNVy/6ffB653v3tVNQO2VFlmo/sdKPJ+p1V1sqquVscsnGQ6zGexUuAhVS3FSYzNgH+par6q/gAsqxLvHFWdpqplwNs4J5yP+izfXtzWFz+2XdO+z8epgd/mbWqtohfQCHjLd6J77BkOrMOp9W4Wkdki0rnKfIXAbcAkEfnOZ/obqppRy+tnd9ZKn70419Vy3ZrpS3XtX71OOqq6FOcf4o4aZumP8w90PnAUzh/K1w3uh9laVS92DzDgHJhOc8/2z8P5om2uZv3/wDmz+MStPtcUR0U8Pn/AG6opfw24DueLM7Ga8qrzTgPeFKcJ6O/ugbUVsEtV833mXY9zJlaJO89k9h8gLmR/k14OcJTvAR3nYNyihnhmAcfjJJ5ZOEn+OPc1q459qcr3YFKA8zesyS9V/jEm+Ba6/yircGpVK1V1gx/bb4XzmXmtd6d5bVfVIj/W43U2cCqwXkRmicjR/iykqp/iXPB9F+dAsg7Ix6nFAuzFOYv2Sgf21nHiU5dWODUI32adqt8f388whyqJH6fG73uWXfXvmSS1XAtzv8fvAG/U0lyzE+eEr4KqtsE5ICbi1gpFZJSIfO02leXi/B18m2t3uskQ9rcU+NYuCqn8/atatqOa5VP93HZNBuH8T97rfhZVTQJeBmb4nDACoE5T63Wq2hHnb7OPKiew7mfvPZHq50c8VVX67FX1GXWuUf8TqC7eSup10nH9Fac3xgEHVXCyv3sgmotzplgnddrH5wJn4fxhXqthvnxVvVlVj8BpQrhJRE44+F2o8BpOO/UUVS2oI8ZSVb1PVXsAv8KpmVyG0yTXVETSfGZvB2yqYVXjcWoCR+M0M3zuTt8AzKpyQE9V1T/UsJ6qSWcWdSedUA1xPg6nKbS62mN1MfyC8w/r1Y7KTZ1Vl9nH/toHIlIpMavqt6p6Bk5z1fs4tTe/qOqzqtpZVbNxkk8c+y8u/0Dls/A+7rTD8QvQVipfm6z6/fHd/w3A2irfkzRVPRX/VPf5P41z8fueasq8ZgBtxL0+VB1xrsO9CzwGZLsHxilUbrIMisPZtnvS8D5OE161PWpV9SacE+4ZIlLTsW8D8Cz7m2O9st11v+97giIiF0vl62lVX97mtc9wjo2HpN4nHXU6CbyF00ujNo8CV1Y9INRiHE4V9EjgvepmEJHfiEgntzljD1BODRf+/KGqa3EO0nfXNa84F++PdC845uE0E3jcL9pXwCPidDboDVxBzc1KU3AOsPfjXF/xxv8R0EVELhWRePc1SES617Cer4CuONc/5rlNEzk4NczZNSyzFac5Itjfw7dwrtFVd7DfCmSKSGOfaeOBe0Qky73e9xdqb5ZbDPQUkb4ikoRzoRUAEUlw/5kbu00wefj5HXH/fr3E0Q54EafpZ7c7yzicE53W4nSlvhkY68+6q2yj4gXMw6mN3Ob+zY/HuY5RU41jHpAvTseKZBGJdWMe5GcIW3GuUXjjuQrnf+DiKrWtSlT1J5zrJ2+KyEnebeOcgHkl4NR6tgNlIjKKytdqg+mwtq1OE7d3PTW5Duck8TMRyRaRJiJyn3tMinG/u/+Hc83Vl7c2Uuw7UVVf18o9Qau+vM1r9wLDxLkNpTVUXBev6dhQSb1POq77ObDprBJVXYJz8LvVz3VOxDloTqyl1tEZ+BSnmWMu8Jyqfl7DvH5Rp/dNbR0IvFrgNEHk4VyzmsX+GtmFQHucs9aJONcfPq1he8U4SfVEfK5buU1vJ+NU83/BaSL5G84/UnXr2QcsBLw9rcD5TNar6rYa9uFt9+dOEfH3OsdBU6dn06dazfUXVf0RJ8mscZuHWgEP4nRU+R5YgrNfD9ay/hU438FPgZVA1d5ClwLrRCQPp7PAxQAi0q7KGWRVSTh/k704B/e5wJ99yl/AuYi/BKf2MxmfC9ni9BS7uKa4cVoHCqu82uIkmVHADpzroJe5n1N1+16OU8vui9OlewfONcTG1c1fjUdwEnyuOL36LsRJQr/4nGHfVcOy1+J0m34C2IXT7PgATnP6z+53+Aack43dOB0gJvkZ12EJ0LaVWo7Rbi1lDM5341Oc5tX27u95ON+JYg7s0Vhrz7g6g3K+70cBbYDFIpKPcw33Fyp/P6slh9f827CJyGqcXkfVHrCNMSZYROQXnM4EAR1BRJzONs+oavNArtdfDaWmE3AicjbOmcYB3VuNMSYE7gD+KiILArVCd10PA7cHap0HHYPVdA4kztAzPYBLVXVamMMxxjQgbnNhdU2Gc1R1VKjjCTVLOsYYY0LGmteMMcaETCAHK4xozZo10/bt24c7DGOMqVcWLFiwQ1VrGvXloEVN0mnfvj3z59c1TJMxxhhfIrK+7rn8Z81rxhhjQsaSjjHGmJCxpGOMMSZkLOkYY4wJGUs6xhhjQsaSjjHGmJCJmi7TxjQEHo+HzRvXs2P7Vso0DiToj4YxDZUqKfHQsXtvEhJqe4JCYFnSiTATJ05kyJAhtGxZ7bObTJRbvfx7ZOO3dFv8OAmFW5GaHzljTK08EsfmjuezaPfptGzTjnadeoRku9a8FkHKy8s566yzGDTI3+dfmWiTV1jKEd/cRWLBZks45rDEaBktV79FbGom7//3SX5etSw02w3JVsxB2bSppidLm6gnMcR4Suqezxg/xGgZIjEkpaYyZ8rbdS8QiG2GZCvGL7GxsXTp0oXMzMxwh2KMiSLJKWns2VnTA34Dy5JOhCkrK2PXrl3YIyeMMSEjgsdTHpJNWdKJIEVFRaxZswZVJTc3N9zhGGPC6PL3Czlx3L5whxFw1nstgpSX7z/TWLduHU2aNAljNMYE3uXvF/Lq4lJu/VUCfz8pqWL6xjwPbZ/cy+ejUzi+fRxyX15FWXIctM+I4ff947np6MSK6etyPXT4196K9+mJ0K1ZDHcNTeSMbvF+x+S7La/EWCi6Jx2A9v/MZ/0e5V8jE7nhqMRK8/3p4yL++U0JJ3SI5dPLGgHgUeWJuSWM/a6Utbke4mMgJyOG07rE8eCIpAO2FW0s6UQQj2d/b6SZM2fSr1+/MEZjTHAkxcFT35Rw7aAEcjJqbmx5ZlQSZ/eIo7AUPlldxnVTi0hNEMYMqHxPyQcXJDO4dSy7C5W/fVnC2RMK+eL/hCFt/D+8ebflVfXup3aNhf8sLK2UdIrKlHHfl5LTuPLc988q5l/flPD0qCSObhNHUZmydJuHrzeGpvkq0lnzWgTxvY4zderUMEZiTPD8qm0sfVrEcNeMolrna5wELVJj6NAkhqsGJtA7O4Zpq8sOmK9pstAiNYbuWbG8dFoSCbHwwY9lzFxXRuz9eWzYU7lr+bjFJTR+NI99Jfv/37zb8r6yUysfGi/oGc+a3R6+2bh/++8sK6VJEhzXvnJye//HMq7ol8AlvRPo2DSGns1jOb9XPE+OrFzLefW7Eno8u5eEB/Jo80Q+98woosxT/bXc6audfdmYV3lf3lpaSspDeeQVO8tt3evh8vcLyfpHPmmP5HHMy/uYvf7AzyycrKYTQXxrOuvWrQtfIKbeOX7sgW3/5/WM55pBCRSUKqe+XnBA+eV947m8bwI7CjycM6HwgPI/DEzg/F7xbNjj4dKJB5bPvLzRIcUqwGMnJXHc2AL+NKScga1ia51fVZmxtpzl2z10yaz9PDkuBuJjodQDx7ePo3PTGF5eVMpfj99fQ3lpYSkX9YqnUYL/ozmkJQoX9IrnpYWlHOXWoF5cUMrv+yfw447KiaBlmjBrfRmb8jy0Tq8+3skrSvm/SUU8ODyRs3vEsWizh6snFyLAA9U0wZ1wRCwtU4XXvy/l9qH79+XVxaWc2S2O9EShsFQZ/moB3bNimHpxChlJwltLSznptQK+u6oR3bNq/5xDxWo6ESQuLo7s7GwAtm0LTfdFY8JhWE4cZ3SL45ZPaq7t/H5SEakP55HwYD4nvlaACNx4VM3DtRSVKffNKiavGE48wkkMYwbE8/J3JXjcVoQfd5Tzxc/lBzTRebflfT0wq/iA9Y8ZkMCbS0vJL1Z+3FHO1xvL+V3fA68dPXlKEoWl0PbJvXR9Zi+j3y/k9e9LK9ViHv2yhLO7x3HnsES6ZDo1oXuPS+SxuSWUlB9Y24kR4ZLe8bz2fWnFtK17PXyyuozRfZx9eeuHUvKKlbfOSWZgq1g6NY3h7mMTOaZtLC8sKD1gneFiNZ0IkpqayoABA5gyZQp79+6tewFjXLXVOlLipdbyZikxtZa3bVx7+aH624mJ9HxuH5N+KqV/ywPPwh8a4XQI2LLXwx2fFnNuj7iKWoavk18rIEagsAyaJAlPnpLIyE7OfKP7xHP3jGKmrSpjVOd4/rOwlAEtY+hXZXvebXk1TT6wFjS4tXMgH7+0lB93eDita9wBzXAA3ZrFsuQPjfhui4cvfi7nq41l/P7DQp78OoY5v2tEcrzww7Zyzu9ZuUZzXPs4isqKWb3LU22tZHSfeP72ZQkLN5fTv2Usry8ppXkj4cQjnHm/3VTOlr1KxqP5lZYrLofk+MgZo8+SToTxXtcpLy9nz549NG7cOMwRGRMcXTJjuWpAPLd/WszUi1MOKM9OFTo1jaFT0xjevyCGLk/vpV/LWI7NqXzYeuWMZAa0iiUjyUmgvjJTYjinh9MsdsIRcYxbXMqDIyr3QPPdVl3GDEjg+fklbNijvH5Wco3ziQj9WsbSr2Us1x+VwBc/lzHslQIm/FDK6L6HNrhm96xYBraKYdxiJ0mPW1zKJb3jiY1xEopHoXtWDBPPP/CzTPG/M1/QBbV5TUReFpFtIrK0hvJuIjJXRIpF5Baf6V1F5DufV56I3OiW3Ssim3zKTg3mPoTS9u3b+fTTTyveb9myJYzRGBN8fz0ukV/yPby4oPahfZqlxHDtoARumFp0wI3TrdOdhFE14XhdNSCeD1eU8cL8UgrLlAt7HfoR+JLe8azc6SEtEU7q6P81ku7NnNi27XNi79k89oAL/LPWlZEcBx1rSX6j+yQwfmkpCzeXs3irh8v67N+Xga1iWbPbQ3oiFcna+2qVFjlXUoIdyVhgZC3lu4AbgMd8J6rqT6raV1X7AgOAAmCizyxPestVdUpgQw6f8vJySktLEXe4+rVr14Y5ImOCK6tRDHcck8g/v657PLnrBjsX7d9cenC9sYa2i6NrZgy3TC/igp7xpCUeelNTeqKw6aY0vr86lZgaHitx9oQCHv+qmLkbylif6+GrDWVcOrGI+Bj4dRenlnbn0ATeXV7Go18Us2JnORN+KOXeWcXcfHQCCbE1x3dhrzh2FypXTCqkf8sYejXfn/gu7h1Ph4wYfv1GAZ+sLmNdrtPb7pE5xbz/Y+Rc0wlq0lHV2TiJpabybar6LVDbJ3ICsFpV1wc6vkjj7b2WkZEBwJQpDSafGlOjPx2dQLOUuhNBdmoMl/WJ5y8zi2vsWlyTK/vHU1LOAR0IDkXjJKk1cY3sGMfHq8s4a0IhXZ7Zy7lvF5IQC7MuT6GHe63m1M7xvHx6Eq8uLqXXc/v407QirhmYUKmXXXUyU2L4dZc4vtvi4bLelWtsSXHCrMtTGNgylt99UEiXp/dy1oRC5v1STk7jyKnpSLDH+BKR9sBHqtqrlnnuBfaq6mPVlL0MLFTVZ3zmvRzIA+YDN6vq7hrWOwYYA9CuXbsB69dHdt7atGkTbdq0oUOHDqxdu5a+ffuyaNGicIdlIsiCBQsY8OGIcIdR79w2vYjpa8pYdFVquEOJOAtOm8GyOZPYs3Mr1z3w7wPKRWSBqg4M1PYiJ/1VQ0QSgNMB3zG3nwc6An2BzcDjNS2vqi+q6kBVHZiVlRXMUAPCW9MREZKSkvjll1/CHJEx9dueIuXbTeW8uKCEPw0J3dMxTc0ivffaKJxazlbvBN/fReQl4KNwBBYMKSkptGzZkuTkZDIzM+1eHWMO0xlvFvDNpnIu6BXPJb0jqAtXFIv0pHMhMN53goi0VNXN7tvfAtX2jKuPMjMz6dmzJwUFBaSnp7Np0yaKiopISrJBAo05FMG4v8gcnmB3mR4PzAW6ishGEblCRK4Wkavd8hYishG4CbjHnSfdLWsEnAS8V2W1fxeRJSLyPTAc+FMw9yHUVBURoXv37gB89dVXYY7IGGMCJ6g1HVW9sI7yLUCbGsr2AQc8QlNVLw1MdJFn9erVzJgxg86dO/PQQw/x8ssvs29fw3uehjEmekV0R4Jo4/F4Km58O/LIIwH46aefwhmSMcYElCWdCOLbe61p06akpKTw7rvvhjkqY4wJHEs6EcRby/GOSBAXF8eyZcvCGZIxxgSUJZ0I4q3pxMQ4f5ZWrVqRn59/wFhTxhhTX1nSiSBNmjShVatWJCc7o9d27twZVWX16tVhjswYczju/LSI7MfykfvyGPtd3ePMNWSRfp9OVGnZsiWdOnWqqOn069ePDz/8kM8++4xOnTqFOTpjDt/l7xfy6uIDh1ocf3YyP+4o575ZJfRrEcPCKsPVLN5STt8XnJ6cG/6USpv0mIrHBaz9YyrtMyL3/PmbjWU8+mUJ75+fzFFtYmmcKJR5lCfmlvDfRaWsz/XQtnEMNx6VwLWDK4+asGJnOddPLWLO+nJS4oVzesTx+MlJB/XU00hjSSeCeDwePB4PsbHOoIDHHnssYD3YTMMyrF0sE86t/CyajCTh0S/KyUoRftzhqXhQmdcLC0rIaSys3xOcpuaScq11dOfDsXKXhxih0kPi7v6siBcXlvLib5Lo0yKWuRvKGfORMzDole6gpHtLlBPGFdA7O5avrmjErkLl/z4oJLeokDfPOfCZOfVF5J4eRKHFixfzxRdfsGPHDgCGDRtGXFycjUhgGpSEWGiRGlPplRTnHPDTE3EeuubzfJ2CUuX1JaVc0W9/LWBdrodhrxQA0OFfe5H78jh+rFMTUlUe+6qYI/6VT8IDeXR8Kp9/fl358dPt/5nPPTOKuGZyIZl/z2fYK86ycl8e//u+cvPXieP2cfn7hRXvP/ixlH4v7CXloTwyHs1j8Et7WbS5vNp9vfz9Qi6dWIRHnXXLfXkAvLq4lJuPTuC33eM5okkMF/eO5/f9Enhozv4431hSyo4C5Y2zkunbIpYRHeJ49tQk3vqhjLW7PQf3oUcQq+lEkKq91xISEujatStLlzaYkX5MENz4cRHfban+oBdMfVvE8s+RgT8hGjMgnl+/UcDjpySREi+8ubSUVmkxDMvZX/Npmy58cEEyZ7xZyLzfN6JtY6moqTz3bSl//ryYf41MYnj7WD5bW86NHxeRliBc0X9/4nrqmxJuOjqBuVekUObnMXzLXg/nvl3IgyMSObdHPEVlyqItHuJqOH3/18gk+rWI4eZPitl40/4mw6IySKpy9E2Oh/V7lPW5HnIyYvhyQzlHt4mlcdL+GtjJHeOIEfhyQxkdmtTPAUwt6UQQ3/t0vBISEio9TdSY+m7munJSH86reN86PYafrtt/QB7aLo426TG8tbSU3/VL4MUFpVzZv/JgnbExQtNk5/8kq5HQInX/Uf/RL4u5fnBCxbNzOmfG8tMODw/NKa6UdAa1juXe4w8uaW7OV0o9cF7P+IrrSN2zan6CaOMkqUgavjGO6hzHU9+UcEKHOHo1j2HepnJeXuRc6/ol30k6m/M9tEit3OQXH+vs9+b8+tuj1ZJOBKla0wGnR1thYSHbt2+nPjyewYReMGobwXRUm1hePXP/NZ3qaglX9o/npYWlDGgVy3dbyvnoomSWbqu7OpJXrGzMU47NqZwIjmsfy7++KaGgVEmJd/6/Brfy/3HTXr2zYzilYyy9ntvLSR3jOD4nlrO6x9P2IB+S9q+RSVz9USF9X9iHAK3ShCv6xfPolyXE1N8+An6xazoRpLqaTv/+/QH45JNPwhKTMYGWHAedmsZUvKrreXZZnwQWbSnnpmlF/LZ7HM1SAn+oqq4HmABVb4sr9cl1sTHC1ItTmDG6EYNaxfLu8jK6PLOXj1Yc3OOgmyYLE85NofDuNNbdmMq6G1MrEtcRTZyfLdNi2LK3cjCl5cquQqVlWv3NTJZ0IkjLli1p1aoVKSn7e6Ycd9xxAMyZMydcYRkTck2ThXN6xPPZ2nKu7F/9tQvvNZxyn6SQnii0SRdmr698jWvWunI6NJGKWk5NmjcSfvFpuiouU5Ztr1zDEhEGt47lrmGJzP5dI47LieWV7w4u6fjuQ5v0GGJEGL+0lGNzYslq5ByWj2kby9yN5eQV749n+poyPArHtK2/jVSWdCJIu3btaNu2LY0a7X8GyIgRzqOJFy9eHK6wjAmLl05LYvutqYzoUP0BNqexECMwZWUZ2/Z52FPkHJzvHJrI0/NKeGlBCSt3lvPC/BKen1/CXUMT69zmiUfE8e8FJczdUMbSbeVc/kEhJeX7D/pfbSjjgVnFfLOxjJ/3ePhsTRnfb/XQo9n+Q2m3Z/byzLzabwD9dlM5b/9QyupdHuZuKOOcCQV8t6Wcp3yaSi86Mp5mKcJF7xayeEs5n68t49opRZzfM44OTervobv+pssGqLS0lLKyskrTUlJSaNy4MXv27AlTVMaER1KcVHSlrk52agyPnJDIo18Wc+O0Ioa1i2Xm5Y34w8B49pUoD39RzDVTlLbpwqMnJlbqRFCTx05O5MoPlVP+V0DjJOGuoYls37c/6TROFOZuLOfZb0vYXaS0SBUuPjKePx+3P6H9tNPDjoLarz8Vlyv3zSpm9W4PCbFwbE4cX/1fI47M3n+dKTVB+PTSFK6fWsTR/91HcrxwTvc4njilfl3Dq0qiZVyvgQMH6vz588MdRq1mz57Ncccdx6BBg5g3b17F9EsuuYTZs2fz888/hzE6EwkWLFjAgA9HhDsM04AsOG0Gy+ZMYs/OrVz3wL8PKBeRBao6MFDbq791tAaout5rAD179mTDhg3k5uaGISpjjAmcYD+u+mUR2SYi1d7dKCLdRGSuiBSLyC1Vyta5j6X+TkTm+0xvKiLTRWSl+7NJMPchlKrrvQb7R52eMGFCyGMyxphAqjXpiMgVInKrz/tNIpInIvkicrUf6x8LjKylfBdwA/BYDeXDVbVvlardHcBnqtoZ+Mx93yDUVNM55phjAKf5zRhj6rO6ajpXAy/7vN+mqulAFnBhXStX1dk4iaWm8m2q+i1wMP0NzwBedX9/FTjzIJaNaDXVdI4++mgAvv/++5DHZIwxgVRX0hFV3enz/m0AVS0CkqtfJGAU+EREFojIGJ/p2aq62f19C5Bd0wpEZIyIzBeR+du3bw9mrAHRsWPHA+7TAYiNjSU1NdU6EhhQxSPW6dQEhkoMaGgHD60r6WT4vlHVhwFEJAZoFqSYvIaqan9gFHCtiBxbdQZ12qNq7H6nqi+q6kBVHVgfhpDp0KED2dnZByQdcO7h2bNnD+XloR/Y0USOBCmjIKNruMMwDURJcgu0pCCk26wr6XwiIg9WM/1+IKjjsqjqJvfnNmAiMNgt2ioiLQHcn9uCGUcoFRQUUFxcXNHM5mv48OEArFmzJtRhmQjSun1nVh31CHub9LQajzksnpgE1vW8gd3bNqEeD7Fxofk+1bWVW4H/iMgqwHtLfB9gPvD7YAUlIo2AGFXNd38/GSfRAUwCRgOPuj8/CFYcofb555+zbNkyMjMzDyi7+OKLefbZZ1m+fDmdO3cOQ3QmEjTNzCQvtyU/9L2XuEZNQOyuB3OIPOUU5W5hx9rv2btnF1kt24Zks7UmHVXdB1woIkcAPd3Jy1R1tT8rF5HxwPFAMxHZCPwViHfX/W8RaYGTwNIBj4jcCPTAabqb6F5QjwPeUNWP3dU+CkwQkSuA9cB5/u1q5PP2XvN2kfbVu3dvRISZM2dy+umnhzo0E0FyjujMlrXLmfXW66SkZRAbF3dA5xNj/KGqFO3bS3xiIiN+e2lItulXfUpV14jIDlXNE5F0f1euqrX2cFPVLUCbaorycGpU1S2zEzjB3xjqk5p6rwE0atSI+Ph4xo8fzxNPPBHq0EwEERGOOuF0GqU1Yd2K7ynclx/ukEw9FRMTQ0bX3vQbehKZ2a1Css2DacSbCfT3+WkCrLakA5Cdnc0vv/wSypBMhBIRjjzqWI486oD+NcZEtENpELZ6fJDUdHOoV48ePSgvL2ft2rWhDMsYYwLGrkJGkN69e9OyZUuSk6u/Ber4448H4J133glhVMYYEziWdCJIx44dadq0KUlJ1Q9dfvbZZwNOLzdjjKmPDiXpRMezEMJgz549FBUVVTSzVdW5c2eaNQv2PbnGGBM8B5N0pMpPE2BTp05l9erV7N27t8Z5hg0bxqpVq0IYlTHGBM7BJJ3zq/w0AVbbfTpe3bp1Y+XKlTYOmzGmXvI76ajqCt+fJvDq6jINkJ7u3CY1fvz4kMRkjDGBdNDXdETk8iDEYdifdGqr6VxwwQUAfPrppyGJyRhjAumgko6IdASeFpEaHydgDl1d9+kAtG/fnoSEBJYsWRKqsIwxJmDqenLon0TkZxHZIyL5wAyc8dN+cJ8gullE/hSSSKPAr371K7Kzs0lMTKx1vjZt2rB9+/ZqR6M2xphIVldN5ypgkKo2VtU0Vc1R1SdUtZn7BNE+OCM9mwDo1KkT6enpJCQk1DrfgAED8Hg8fPPNNyGKzBhjAqOupHOfqm6tqdB91s3fAhtS9Nq6dSuFhYV11mCuuuoqANavXx+KsIwxJmBqTTqqWmcXKX/mMf6ZNGkSGzdupKSkpNb5jj32WJKTk62mY4ypd+q6pjPIfeaN9/1lIvKBiDwlIk2DH1508acjAUB8fDzt27e3MdiMMfVOXc1rLwAlACJyLM4D1MYBe4AXgxta9PHnPh2vlJQUNm7cSF5eXrDDMsaYgKkr6cSq6i739/OBF1X1XVX9M9ApuKFFn4NJOiNGjABgwoQJQY3JGGMCqc6kIyLeB72dgNNl2qvOB8CJyMsisk1EltZQ3k1E5opIsYjc4jO9rYh8LiLLROQHEfmjT9m9IrJJRL5zX6fWFUd94c8wOF4XX3wx4FwHMsaY+qKuxDEemCUiO4BCYA6AiHTCaWKry1jgGZwmuersAm4AzqwyvQy4WVUXikgasEBEpqvqMrf8SVV9zI/t1yunnHIKWVlZdXaZBujTpw9xcXHMnz8/BJEZY0xg1NV77SHgZpzkMVT3j7kfA9xS03I+y8/GSSw1lW9T1W+B0irTN6vqQvf3fGA50Lqu7dV3nTp1Ijk5mdjYWL/mb9++PTt27KjxUQjGGBNp6uq99hdV/VpVJ6rqPp+ircAjwQ2tIob2QD/At3/wdSLyvdt816SWZceIyHwRmb99+/Zgh3rY1q9fX+vzdKq69dZbKS0t5aeffgpyZMYYExh1XTwYKiIP+U5wx12bReXrO0EhIqnAu8CNqurtpvU80BHoC2wGHq9peVV9UVUHqurArKysYId72N577z22bdvm9/A2w4cPB2DWrFnBDMsYYwKmrqRzOtBHRJ4AEJHOwJfAv1X1/mAGJiLxOAnndVV9zztdVbeqarmqeoCXgMHBjCOUDqb3Guxvjnv88RrzrjHGRJS6rukUAb8F2ovIeOBT4FZV/XcwgxLnqPtfYLmqPlGlrKXP298C1faMq4/8ebSBLxEhMzOTNWvW2HUdY0y9UNc1nZuA63Gup5wMLAI6iMhNblmt3EQ1F+gqIhtF5AoRuVpErnbLW4jIRuAm4B53nnTgGOBSYEQ1XaP/LiJLROR7YDjQYEa5PpTEMWTIEMrLy/niiy+CEJExxgRWXV2m03x+f6qaabVS1QvrKN8CtKmm6Aug2jYmVb3U3+3XNwdb0wHnfp133nmHV155hWHDhgUrNGOMCYhak46q3heqQAycd955/OMf//C7yzTA6aefjogwc+bM4AVmjDEBUueAn3WtwJ95jH+OOOIIEhISDqqmExMTQ7du3cjNzbXrOsaYiFfX0e2V2gpFJBanB5kJgOXLlx/UfTpe1113Hbt372bVqlVBiswYYwKjrqTzlTvO2ffuxfyPRGSUiMxyL+SvA74KfpjR4e233yY3N9fvLtNeJ598MgBvvfVWMMIyxpiAqeuazhi3N1kGEIszMsA7OF2VfwIKVDXyb/WvJw72Ph2vTp06ER8fzzPPPMM999wTjNCMMSYg6hwp2h0JwDsawFoRGayqnwQ3rOh0MKNMV9W5c2eWLVvGvn37aNSoUaBDM8aYgDjoo5uq3lF1moj0DEw40e1QazoAp57q3Mb06quvBjQmY4wJpIM/pa7eawFaT1Tzd8y16lxzzTWAXdcxxkS2QCWdgz81NwcYM2YMjRs3PqSaTocOHUhNTWXBggVBiMwYYwIjUEnHbhAJgJycHGJjYw8p6QCce+65FBQUsGPHjgBHZowxgRGopGMCYP78+RQXFx9y0hkzZgyqyiefWD8PY0xkClTSKQnQeqLaW2+9xb59++qesQaDBg0iIyODJ598MoBRGWNM4NTZZdrLfUJnZyDJO819HDWqOiTwoUWfw+m9BhAbG0vz5s2ZP38+e/fuJTU1NZDhGWPMYfOrpiMivwdmA9OA+9yf9wYvrOjkvU/nUJMOwDnnnAPA888/H5CYjDEmkPxtXvsjMAhYr6rDcUYmyA1WUNHqcGs6AH/84x8BeOONNwISkzHGBJK/SafIfYooIpKoqj8CXYMXVnQKRNJp3rw5WVlZLF261EadNsZEHH+TzkYRyQDeB6aLyAfA+mAFFa3uuOMOUlJSDns9J510EmVlZcyZMycAURljTOD4lXRU9beqmquq9wJ/Bv4LnFnXciLysohsE5GlNZR3E5G5IlIsIrdUKRspIj+JyCoRucNnegcR+cad/paIJPizD/VBq1atiImJOayaDsDDDz8MwOzZswMRljHGBIzfXaZFpImI9AbygY1ALz8WGwuMrKV8F3AD8FiVbcUCzwKjgB7AhSLSwy3+G/CkqnYCdgNX+LsPkW7mzJmUlpYedtLJyclhyJAhvPfeewGKzBhjAsPf3msPAN8DTwOPu6/Hal2Iii7Vu2op36aq3wKlVYoGA6tUdY2qlgBvAmeIczQegfN4BYBX8aPGVV+8+eabh3VzqK/+/fuzaNEiPv744wBEZowxgeFvTec8oKOqHqeqw93XiCDG1RrY4PN+ozstE8hV1bIq06slImNEZL6IzN++PfIf+xOIjgRel19+OQCPPVbnuYExxoSMv0lnKc6D3OoVVX1RVQeq6sCsrKxwh1OnQPY2GzRoEGlpaXz55ZcBW6cxxhwuf5POI8AiEZkmIpO8ryDGtQlo6/O+jTttJ5AhInFVpjcIgazpAJx44okUFRUxderUgKzPGGMOl79J51WcC/iPsv+azuPBCgr4Fujs9lRLAC4AJqlTFfgcOMedbzTwQRDjCKlAJ5277roLgH/84x8BWZ8xxhwuf8deK1DVpw525SIyHjgeaCYiG4G/AvEAqvpvEWkBzAfSAY+I3Aj0UNU8EbkOZ7idWOBlVf3BXe3twJsi8iCwCKf7doPw97//nTfeeCNgSWfgwIG0b9+etWvXoqoBW68xxhwqf5POHBF5BJgEFHsnqurC2hZS1QvrKN+C00RWXdkUYEo109fg9G5rcLKysgI+isDtt9/OH/7wBxYtWkT//v0Dum5jjDlY/jav9QOGAA9zEF2mzcGZNGkS5eXlAa2RnHfeecTHx/PnP/85YOs0xphD5VdNxx3k0wTZW2+9hcfjCWjSadq0KS1atGDq1Kn2uANjTNj5e3NohojcICJPiMhT3lewg4s2ge5I4PX73/8eVeWhhx4K6HqNMeZg+du8NgVoDywBFvi8TAAFK+ncdtttxMTE8OqrrwZ0vcYYc7D87UiQpKo3BTUSE7RHESQlJdGvXz8WLFjA8uXL6d69e1C2Y4wxdfG3pvOaiFwpIi1FpKn3FdTIolCwajpARUeCJ598MuDrNsYYf/mbdEqAfwBz2d+0Nj9YQUWrl19+GQhO0jnjjDM46qijmDlzZkVyM8aYUPM36dwMdFLV9qrawX0dEczAolF6ejoQnKQDcN1117Fy5UomT54clPUbY0xd/E06q4CCYAZiYNy4cUFd/znnnEN8fDxXX311ULdjjDE18bcjwT7gOxH5nMojEtwQlKii1Ntvvw0Er6bj7VAwb948vv32WwYNGhSU7RhjTE38rem8DzwEfIV1mQ6aYHYk8Hr00UcBuPXWW4O2DWOMqYm/IxLYDR4hEIqkM3z4cLKyspgzZw579uyhcePGQduWMcZU5e+IBJ1F5B0RWSYia7yvYAcXbUKRdABuueUWPB4PN91kt14ZY0LL3+a1V4DngTJgODAO+F+wgopWoerKfMstt5CTk8PcuXODdkOqMcZUx9+kk6yqnwGiqutV9V7g18ELKzp98IHzPLpg13RiYmJ44IEHWL58OR9//HFQt2WMMb78TTrFIhIDrBSR60Tkt4ANVxxgcXHOJbZQPGzt/PPPJyMjgyuuuCLo2zLGGC9/k84fgRTgBmAAcCnOo6JNAD3zzDNAaJJOQkICRx99NJs3bw76/UHGGOPlV9JR1W9Vda+qblTV36nqWar6dV3LicjLIrJNRJbWUC7uYxJWicj3ItLfnT5cRL7zeRWJyJlu2VgRWetT1tf/3Y1s77//fki398ILLyAi3HbbbSHdrjEmetWadERkqIhc5vP+HRGZ4b5G+LH+scDIWspHAZ3d1xiczgqo6ueq2ldV+wIjcEZD+MRnuVu95ar6nR9x1Auh6r3m1bZtW4YPH87WrVt5/fXXQ7JNY0x0q6umcx+VB/bsCtwK3AvUeXqsqrOBXbXMcgYwTh1fAxki0rLKPOcAU1W1wQ/DE+qkA/DKK68gItx8880h26YxJnrVlXTSVXWZz/uVqrrATSZpAdh+a2CDz/uN7jRfFwDjq0x7yG2Oe1JEEmtauYiMEZH5IjJ/+/btAQg3uLzdl0OZdNq1a8dvf/tbtm7dyvz5NnC4MSa46ko6Gb5vVPUsn7fZAY+mCrfWcyQwzWfynUA3YBDQFLi9puVV9UVVHaiqA7OysoIaayDExDh/jlAmHXBqO02aNOG+++4L6XaNMdGnrqTzo4gccD+OiPwG+CkA298EtPV538ad5nUeMFFVS70TVHWz2xxXjHPT6uAAxBERpk2bVvdMQZCens7111/PRx99xAsvvBCWGIwx0aGupPMn4AkReUVErndfY4En3LLDNQm4zO3FNgTYo6qbfcovpErTmveajzjVgTOBanvG1UfhaF7zuvbaa4mJiakYIscYY4Kh1qSjqquA3sAcoL37mg30VtUVda1cRMbjPG20q4hsFJErRORqEfE+0GUKsAbneT0vAdf4LNsepxY0q8pqXxeRJcASoBnwYF1x1BePPPIIEJ6k07x5cy655BL27t1rXaiNMUEj0TL21sCBAzXSL5T37duXxYsX8/jjj4dlMM7i4mIaN26Mx+Nh165dpKbaoBPGRDsRWaCqAwO1Pn9HJDAhEM7mNYDExETuvvtuSktLueiii8ISgzGmYbOkE0Ei4VrKPffcQ6dOnZg5cyb1oZu5MaZ+8fd5OpeKSFqVab8JTkjRK9w1He+2P/jgAwoLC+3ajjEm4Pyt6TwNzBGR7j7T7g9CPFHNew0lnEkHoEePHlxzzTWMHTuWZ599NqyxGGMaFn+Tzlrg/4B3RORcd1p4j4wN0JQpU4DwJx1wmtliY2O5+eab2bt3b7jDMcY0EP4mHVXVhcBxwBgReQyIDV5Y0SmSehJmZWVx1113UVxczJlnnhnucIwxDYS/SWczgKruAE4BFOgZrKCi1R133AFERk0H4P7776dt27Z89tlnIX/sgjGmYfI36Yz1/qKqHlW9FWe0ABNAM2bMACIn6YAzNI+IcPnll1NaWlr3AsYYUwt/k86d1Uy7I5CBmMjovVZV9+7deeCBB9izZw+PPvpouMMxxtRzcbUVisgo4FSgtYg85VOUDpQFM7BoFI7n6fjj7rvv5ocffuD+++9n4MCBjBo1KtwhGWPqqbpqOr/gPMStCFjg85qEc23HBFAkdSSo6umnnyYuLo4zzzyTrVu3hjscY0w9VdeAn4tV9VWgk6q+6vN6T1V3hyjGqJGd7TyiKNJqOgCZmZk8+OCDlJSUMGzYsHCHY4ypp/y9pjNYRKaLyAoRWSMia0VkTVAji0KTJk0CIjPpANx8880cc8wxrFy5kjFjxoQ7HGNMPeRv0vkvzjN0huI8sXOg+9MEUCR2JKjqs88+IyMjg5deeomJEyeGOxxjTD3jb9LZo6pTVXWbqu70voIaWRT6wx/+EO4Q6pSYmMjMmTOJj4/nxhtvZPdua2U1xvjP36TzuYj8Q0SOFpH+3ldQI4tCX331FRDZNR2APn36MGPGDDZv3swFF1xASUlJuEMyxtQTtXaZ9nGU+9P3QT4KjAhsONGtPjSveQ0dOpRnnnmGq666il/96ldE+gPyjDGRwa+ajqoOr+ZVZ8IRkZdFZJuILK2hXETkKRFZJSLf+9aeRKRcRL5zX5N8pncQkW/cZd4SkQR/9qE+iNT7dGoyZswYevfuzYIFCxg9enS4wzHG1AP+Pk8nW0T+KyJT3fc9ROQKPxYdC4yspXwU0Nl9jQGe9ykrVNW+7ut0n+l/A55U1U7AbsCfOOqF+pZ0AL7++msyMzMZN24cf/7zn8MdjjEmwh3M2GvTgFbu+xXAjXUtpKqzgV21zHIGME4dXwMZItKyppnFORqPAN5xJ70KnFlXHPVFp06dwh3CQUtOTmbJkiU0atSIBx98kOeeey7cIRljIpi/SaeZqk4APACqWgaUB2D7rYENPu83utMAkkRkvoh8LSJnutMygVx3+1Xnr/feffddoH7VdABatmzJ/PnzSUxM5Pbbb2fp0mpbU40xxu+ks09EMnE6DyAiQ4A9QYvKkaOqA4GLgH+KSMeDXYGIjHET1/zt27cHPsIAq08dCarq1q0by5cvJz09nVNOOcUSjzGmWv4mnZtwxlvrKCJfAuOA6wOw/U1AW5/3bdxpqKr35xpgJtAP2InTBBdXdf7qqOqLqjpQVQdmZWUFINzguuiii4D6mXQAOnTowLRp08jLy6Nv375Mnz493CEZYyJMnUlHRGJxnhh6HPAr4Cqgp6p+H4DtTwIuc3uxDcG5CXWziDQRkUR3+82AY4Bl6lQFPgfOcZcfDXwQgDgiwqJFi8IdwmHr1asXY8eOxePxMHLkSKZOnRrukIwxEaTOpKOq5cCFqlqmqj+o6lJV9etpXiIyHpgLdBWRjSJyhYhcLSJXu7NMAdYAq4CXgGvc6d2B+SKyGCfJPKqqy9yy24GbRGQVzjWe//q3q5GvPjev+Tr77LN57733UFV+85vfVFyrMsYYf28O/VJEngHeAvZ5J6rqwtoWUtVany7q1lyurWb6V8CRNSyzBhjsR8z1Tn3sMl2TM888k0mTJnHGGWdw7rnnMnnyZHsOjzHG72s6fYGewP3A4+7rsSDFFLUaSk3H6ze/+Q0zZsygWbNmnH322UyePDncIRljwiyoIxKYg9OnTx+g4SQdgOOOO46lS5fSvXv3ilpPJD+szhgTXP6OSNBYRJ7wdj8WkcdFpHGwg4s2b7zxRrhDCIrmzZszc+ZMWrduzTvvvEPv3r0pKCgId1jGmDDwt3ntZSAfOM995QGvBCuoaNXQmtd8paWlsWrVKgYNGsTSpUtp3bo1P/74Y7jDMsaEmL9Jp6Oq/lVV17iv+4AjghlYNDrttNOAhpl0AOLj45k3bx5jxowhNzeXXr162XUeY6KMv0mnUESGet+IyDFAYXBCik6qyooVK4CGm3S8XnjhBcaNG0dSUhJnnXUWTz/9tF3nMSZK+Jt0rgaeFZF1IrIOeAbnJlETIL4H3YaedAAuvfRS1q9fz8knn8wNN9xATk4Oq1atCndYxpgg87f32mJV7QP0Bnqraj/sAW4B5b1HJ5pkZmYyadIkRo8ezYYNG+jatSsPP/xwuMMyxgSRvzUdAFQ1T1Xz3Lc3BSGeqOWbdKKhpuMlIowdO5bXX3+d+Ph47r77bnr27Mkvv/wS7tCMMUFwUEmniug5MoaAiHD00UdX/B5tLrroIjZv3syAAQNYtmwZ3bt3Z8qUKeEOyxgTYIeTdOzKbwDFx8fz0ksvAdGZdACaNGnC/PnzefHFF2nRogW//vWvOfHEE/n222/DHZoxJkBqTToiki8iedW88tn/FFETINaDy3HllVeyZMkSHn74YT7//HMGDx7MqFGjyM/PD3doxpjDVGvSUdU0VU2v5pWmqv4OFmr8sG/fPk4++WQgems6vhISErjzzjv57LPPyM7O5uOPPyYzM5N77rknKjtdGNNQHE7zmgmg8vJyNm/eDFjS8XX88cezZcsW7rvvPkSEhx56iO7du/PFF1+EOzRjzCGwpBMhorX3mr/+8pe/kJuby6WXXkp+fj7Dhg1jyJAhTJgwIdyhGWMOgiWdCBFtN4ceiuTkZMaNG8eqVat4+OGH+fbbbzn//PNp06YN48ePD3d4xhg/WNKJEHadwn8pKSnceeedrF69mmHDhrFp0yYuuugisrKy+M9//hPu8IwxtbCkEyHi4+MZPnw4YDUdf7Vv357Zs2ezbt06Tj75ZHbu3MmVV17JUUcdxdixYykstOEBjYk0QU06IvKyiGwTkaU1lIuIPCUiq0TkexHp707vKyJzReQHd/r5PsuMFZG1IvKd++obzH0IlfT0dB5//HHAks7BysnJYdq0aWzfvp2nn36a3Nxcfve735GamspJJ53E0qXVfv2MMWEQ7JrOWGBkLeWjgM7uawzwvDu9ALhMVXu6y/9TRDJ8lrtVVfu6r+8CHXS4NOTn6YRCZmYm1113HcuXL+e+++6jSZMmfPrppxx55JG0a9eOxx9/3JoxjQmzoCYdVZ0N7KplljOAcer4GsgQkZaqukJVV7rr+AXYBmQFM9Zw27JlS0Xzmjk8MTEx/OUvf2HHjh1MnjyZfv36sXHjRm655RY6dOjAvffeyzfffBPuMI2JSuG+ptMa2ODzfqM7rYKIDAYSgNU+kx9ym92eFJHEmlYuImO8j9jevn17IOMOuPLycvLynLFUraYTOKeeeioLFy5k586dPP/883Tr1o3777+fIUOGkJGRwWWXXcbKlSvDHaYxUSPcSadWItISeA34nap620XuBLoBg4CmwO01La+qL6rqQFUdmJUV2RUlu08nuJo0acLVV1/NtGnTmD9/PiNGjKCgoIDXXnuNLl260KxZMx5++GF27twZ7lCNadDCnXQ2AW193rdxpyEi6cBk4G636Q0AVd3sNscVA68Ag0MYb9DYfTqh079/fz777DOKior43//+x5AhQ8jLy+Puu+8mOzubfv36cckll7Bw4cJwh2pMgxPupDMJuMztxTYE2KOqm0UkAZiIc73nHd8F3NoP4hyZzwQaRNckq+mEXkxMDBdffDFz586luLiYhQsXcuutt7Jy5Upef/11BgwYQEpKCscccwzPP/88paWl4Q7ZmHov2F2mxwNzga4islFErhCRq0XkaneWKcAaYBXwEnCNO/084Fjg8mq6Rr8uIkuAJUAz4MFg7kOoNGrUiFGjRoU7jKglIvTr149HHnmEvLw83n33XUaOHElCQgJfffUV11xzDc2aNePss8/m+uuvt8ctGHOIJFqG0x84cKDOnz8/3GHU6uuvv+boo49mypQploAiyPr16/noo49YvHgxkydPrniqaXx8PB07duSEE07gyiuvpE+fPmGO1JjAE5EFqjowUOuzxxNEIGteiyw5OTlce+21gNPLcOLEiYwbN4558+bx448/8uOPP/Lss8/SsWNHBg8eTEpKCueccw4nnXQSsbGxYY7emMhiNZ0IsWLFCnr27ElZWRkff/wxp5xySrhDMn4oLi5mwoQJrF27loULF/Lpp5+yb98+wDl5aNasGb169eLKK6/klFNOoWnTpmGO2JiDYzWdBsrj8VBWVhbuMMxBSkxM5NJLL614X1ZWxuTJk3nvvff4+uuvWb9+PZ9//jmff/45ANnZ2SQlJdGrVy+OPfZYzjzzTLp06RKu8I0JOavpRIhly5bRs2dPAKZNm1bxFFFT//3888+sXr2ab775hjfeeINly5ZRXl5eUR4bG8uJJ55I3759ycjIoEePHpx00kkkJyeHMWpjHFbTaaDsPp2Gq127drRr147hw4dzxx13AE5z6vvvv8/s2bNZu3YtW7Zs4YknnqjULTspKYns7Gx69uzJZZddRpcuXejUqRNpaWnh2hVjDpvVdCLEkiVL6N27NwDTp0/nxBNPDHNEJtRKSkr48MMPmT59OgsXLmTt2rXs3r27Uq0InJpRRkYGrVq1okuXLgwdOpTTTjuNnJwc4uLsPNIEVqBrOpZ0IsSmTZu4/vrrmThxIp9++iknnHBCuEMyESI/P5+1a9eyYsUKnn/+eVavXs327dspKCioNF9CQgIxMTFkZGTQokULcnJy6N69O8OGDWPo0KGkp6eHaQ9MfWZJ5xBFetIBmD17Nscdd5wlHeMXj8fDsmXLWLFiBXv27GHp0qWMHz+eXbt2UVxcfMD8GRkZFBUV0bhxY5o3b067du3o2LEjxx57LAMHDqRVq1bEx8eHYU9MJLNrOg2Ux+OhpKQEsGs6xj8xMTH06tWLXr16VUzzPgjQ4/GwatUq5s2bR15eHvv27eOHH35g6tSp5ObmsnXrVpYsWQLAU089BTjfOxEhKSmJ9PR0mjZtSosWLTjhhBMqOjmkp6fTtWtXS07mkFnSiRALFy7kpJNOAizpmMMXExNDly5dau2O/fPPPzN//nxKS0vJz89n8eLFfPLJJ+zatYvdu3ezdetWli1bxowZM6pdf2JiIo0aNWLo0KF07doVVWXXrl20adOG9u3b06lTJ7p27UpmZqZ9p00FSzoRwnqvmVDz9qqrTV5eHlu3bmXHjh189913zJgxg02bNrFz505yc3PZu3cvX331FZMnT65xQNTExESaN29OXFwc+fn5pKamkpGRQdOmTWnevDm//vWvadGiBeXl5SQkJNC+fXvatGljtakGypJOhLBRpk0kSk9PJz09nc6dO3P00Ufzhz/8odr5VJVt27axbNkyVq9ezfr169m4cSObN2+ma9eu7Nmzh0WLFrF161Zyc3NZt25dxbJvvvlmtesUEWJjY+nZsydNmjSpaCZMS0ujcePGNGnShObNm3PqqafSpEkTCgsLSU1NpWXLlrRs2dKSVoSypBMhoqVDh2mYRITs7Gyys7P9fux6fn4+a9asQUTIzc1l0aJFLF26lB07dlTUpPbt20e7du3Izc3l559/ZteuXZVO0ACee+65GrfhHZQ1LS2N7du3U1JSQnJyMikpKTRq1IjmzZtzwgknkJ6ezrp164iLi6uogfXs2ZPMzEyaN29+WJ+NqcySToSwmo6JNmlpaZVG5j722GP9XjY/P58NGzawefNm0tLS2L17N3PnzmXNmjXk5uayZ88e8vPz8Xg8dOnShfz8fH7++Wdyc3MpKyvD4/FUnOhNmjSpxu0cd9xxzJw5E4B+/fqxadMm0tLSSE1NJTU1leHDh/Pgg87TVW666aZKSS05OZk+ffpUjBj/4YcfEh8fX1GWkpJCVlZWRVIrLS2NitqZJZ0I0bp1ay644ALefPNNSzrG1CEtLY0ePXrQo0ePimmHMkjuvn37KCwsJD8/n6VLl7J582Z27tzJY489Rrdu3bj99tsr5j377LPZtGkTe/fuJT8/n/z8/EotFNOnT2fz5s0UFBRQWFgIwOjRoyuSztlnn33Ada9rrrmGZ599ltLSUhISEoiNja2UlK677jpuvvlm8vLyOOOMM0hKSqr0Ovfcczn11FPJzc3lySefPKD89NNPJzMz86A/l2CypBMhcnJyuOKKKyzpGBNCjRo1olGjRjRr1owOHTpUTJ88eTIJCQmVnmt1zz331Loubxd0cJrLi4qKKrVgzJs3j8LCwoqkVFBQwBFHHFEx/4MPPliprKCggLZt2wJU1M527dpFUVERRUVFFBcXM2jQIAB27NjB/ffff0BM33//vSUdU73S0lL27NkT7jCMMUCLFi1YtmzZIS8vIgcM2Nq3b98a509ISODuu++usbxp06bMmjWrxvJOnTpRXl5OSUlJRUIqKiqiZcuWBx17sAX1cdUAIvKyiGwTkaU1lIuIPCUiq0TkexHp71M2WkRWuq/RPtMHiMgSd5mnpAFUDb744gvOOeccwK7pGBNuLVq0YMuWLeEO46DExMSQlJRERkYG2dnZ5OTkkJCQEO6wDhD0pAOMBUbWUj4K6Oy+xgDPA4hIU+CvwFHAYOCvItLEXeZ54Eqf5Wpbf71g9+kYEzlatGjB7t27qx1OyByeoDevqepsEWlfyyxnAOPUOep+LSIZItISOB6Yrqq7AERkOjBSRGYC6ar6tTt9HHAmMDUY8d944428+OKLB3TTTEtLq+h1snr16gOWa9y4Mc2aNUNVWbNmzQHlGRkZZGZm4vF4WLt2rfVeMyaCeJulevfu3SBG7n733Xfp1q1buMMAIuOaTmtgg8/7je602qZvrGb6AURkDE7tqc47r2vSpk0bmjdvfsBTPZs3b07Hjh0BKrpm+mrRogUdOnTA4/GQl5d3wHpbt25Nu3btKCsrIz8/H3DuKRg2bBhHHnnkIcVqjAmMkSNHcskll1BUVBTuUAIiMTEx3CFUiISkEzSq+iLwIjijTB/KOm655RZuueWWgMZljIlsrVu35rXXXgt3GA1SKK7p1GUT0NbnfRt3Wm3T21Qz3RhjTISLhKQzCbjM7cU2BNijqpuBacDJItLE7UBwMjDNLcsTkSFur7XLgA/CFr0xxhi/Bb15TUTG43QKaCYiG3F6pMUDqOq/gSnAqcAqoAD4nVu2S0QeAL51V3W/t1MBcA1Or7hknA4EQelEYIwxJrDsyaHGGGNqFOgnh0ZC85oxxpgoYUnHGGNMyFjSMcYYEzKWdIwxxoRM1HQkEJHtwPoqk5sBO8IQTiSwfY9Otu/R53D3O0dVswIVTNQkneqIyPxA9sqoT2zfbd+jTbTue6TttzWvGWOMCRlLOsYYY0Im2pPOi+EOIIxs36OT7Xv0iaj9juprOsYYY0Ir2ms6xhhjQsiSjjHGmJCJ2qQjIiNF5CcRWSUid4Q7nmASkZdFZJuILPWZ1lREpovISvdnk3DGGAwi0lZEPheRZSLyg4j80Z0eDfueJCLzRGSxu+/3udM7iMg37vf+LRFJCHeswSIisSKySEQ+ct9Hxb6LyDoRWSIi34nIfHdaxHznozLpiEgs8CwwCugBXCgiPcIbVVCNBUZWmXYH8JmqdgY+c983NGXAzaraAxgCXOv+naNh34uBEaraB+gLjHSfV/U34ElV7QTsBq4IX4hB90dguc/7aNr34ara1+f+nIj5zkdl0gEGA6tUdY2qlgBvAmeEOaagUdXZwK4qk88AXnV/fxU4M5QxhYKqblbVhe7v+TgHoNZEx76rqu5138a7LwVGAO+40xvkvgOISBvg18B/3PdClOx7DSLmOx+tSac1sMHn/UZ3WjTJdp/CCrAFyA5nMMEmIu2BfsA3RMm+u81L3wHbgOnAaiBXVcvcWRry9/6fwG2Ax32fSfTsuwKfiMgCERnjTouY73zQnxxqIp+qqog02L7zIpIKvAvcqKp5zkmvoyHvu6qWA31FJAOYCHQLb0ShISK/Abap6gIROT7M4YTDUFXdJCLNgeki8qNvYbi/89Fa09kEtPV538adFk22ikhLAPfntjDHExQiEo+TcF5X1ffcyVGx716qmgt8DhwNZIiI92SzoX7vjwFOF5F1OE3nI4B/ER37jqpucn9uwznZGEwEfeejNel8C3R2e7MkABcAk8IcU6hNAka7v48GPghjLEHhtuP/F1iuqk/4FEXDvme5NRxEJBk4Ceea1ufAOe5sDXLfVfVOVW2jqu1x/rdnqOrFRMG+i0gjEUnz/g6cDCwlgr7zUTsigYicitPuGwu8rKoPhTei4BGR8cDxOEOcbwX+CrwPTADa4Tzy4TxVrdrZoF4TkaHAHGAJ+9v278K5rtPQ9703zgXjWJyTywmqer+IHIFz9t8UWARcoqrF4Ys0uNzmtVtU9TfRsO/uPk5038YBb6jqQyKSSYR856M26RhjjAm9aG1eM8YYEwaWdIwxxoSMJR1jjDEhY0nHGGNMyFjSMcYYEzKWdIwJExG5Q0QuFpF7RURFpJNP2Y3utIG1rcOY+saSjjEhJo4Y4BTgE3fyEpwbGb3OBX4IdWzGBJslHWNCQETau89vGodzh3hbIEFVt7uzvI870rmIdAT2ADvCEasxwWRJx5jQ6Qw8p6o9gQE4zzXxygM2iEgvnBrPW2GIz5igs6RjTOisV9Wv3d9HAlOrlL+Jk3DOZP9QJsY0KJZ0jAmdfT6/DwbmVSn/CLgU+FlV80IWlTEhZM/TMSbERKQn8KP7vJsKqlogIrcDK8ITmTHBZ0nHmNAbBXxcXYGqvhniWIwJKRtl2pgQE5HpwGU+jw82JmpY0jHGGBMy1pHAGGNMyFjSMcYYEzKWdIwxxoSMJR1jjDEhY0nHGGNMyFjSMcYYEzL/D29d44NuHfb/AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEWCAYAAACT7WsrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABLfUlEQVR4nO3dd3hUZfbA8e+ZSaPX0CFEijSpAVQUEVHBArqKgqiAfV1dy8+1rK6FFcvuumJdYdeGDbGDoigoCNK7NOktgHQCIT3n98e9EyZ9EjKZlPN5nnkyc9ucO5m5577ve+/7iqpijDHG5MUT6gCMMcaUXZYkjDHG5MuShDHGmHxZkjDGGJMvSxLGGGPyZUnCGGNMvixJmFMiIk+KyPshfP81ItKvgPmzROSW0ouo7BDH2yJyWEQWhTqeUyUiHhG5V0RGlOA2zxOR9SJyTESalNR2S4KIjBKRuQXMXykie0TkvmDGYUkCEJFtIrJPRKr5TbtFRGb5vVYRSRSR4yISLyL/FhGv3/xZIpLszj8gIp+LSOMA3vtJEXmyGDH7x3NcRI640/u5877IsXwXd/qsvLZXXqlqR1WdBaeesNzPLtPvMz0uIlP9tq0ick+Ode5xpz/pt41dxd+jEnUOcCHQTFV7Qdb3epO7b9/5HxjdpPK8iBx0H8+LiAT6ZiJSw/1dbHO/mztE5FMR6Z1juSdFpGUx9udl4DZgXAkm/j8D36pqDVXd7caX7X8oIhHu7/kXEakpIrVF5C0R2esmlw0i8nDODYtIUxFJF5FWecz7QkT+dSqBq2oX4FbgiVPZTmEsSZzkBe4pZJkuqlodOA+4Frgpx/y73PltgdrAiyUdZF7xuI/aftP3A2eJSD2/aSOBDUGOpyLY7feZVlfVy/3mbQBuzLF8iX6uIhJWUtsCYoBtqprobrsf8AwwBKgLbAU+8lv+NuAKoAvQGbgcuD2QNxKRSOBH4AzgMqAm0B6YBAxyl/mriJzrrhImIo+JyJkBbv9x4GygL9APeFJEhhSwfKCfY11gTQHbiQQ+x/k9X6SqCTi/6+o4+1cLGAxsyrmuqsYDM4EbcmyzLnAJ8G6AMRZkNVDL/4S1pFmSOOmfwAMiUruwBVV1E/AL0DWf+YeAz4BOItJTRH7PUer4g4iszLmeiNQXka9F5IiIHBKROSJSnP9RKvAlMMzdrhcnqX2Q3woiEiUi77tnkEdEZLGINHTnNRGRKW5Mm0Tk1ny28a2I3JVj2koR+YP7vJ2I/OBu5zcRuSaf7ZwvIr/6vf5BRBb7vZ4jIle4z7eJyAARGQj8FbjWPUv2/3xj3LPAYyLyvYjUz/eTK9hioKqIdHTfuyMQ5U5HnJLot0ATv5JIExGJFJFxIrLbfYxzDz5ZZ60i8pCI7AXeljyqGdzSSmv3+SUistbdn3gReSCPz/Bm4H84JwvHReQpnIP3J6q6RlVTgb8Dff3OdEcCL6jqLvcA9wIwKsDP5gagGXCFqq5W1QxVTVTVT1X1SXeZl4CBON/LN4A1qroggO/N7TiJpr+qHlDVNcD5wD9F5Jwcn9GfRGQjsDHAuMOAzLxmiEhVYKq7zKW+ZAv0BD5U1cOqmqmq61X103y2/y45koS7/2tV9VcReVhENrv/y7UicmWAcfv4Yi/Jk4tsLEmctASYBeT6weUkIu2Ac8nj7MGdXx+4CliuqouBg8BFfovcAEwEUNUn/X5E/wfsAqKBhjgHveL2mzKRk2e9F+OccewuYPmROGdFzYF6wB1AkjtvkhtXE+Bq4BkR6Z/HNj4ChvteiEgHnLPZb9wD6A/Ah0ADnB/K6+4yOS0A2rhJMxznrLaJONUZVYA4YI7/Cqr6Hc5Z8sduCaCL3+zrgNHu+0YQwP+4AO9x8nMd6b72xZCIczDzL43sBh4FzsQ5qegC9AIe89tmI5wz2hics/nCvAncrqo1gE44Z/DZqOqbOP/D+W4cvioJ/+oj3/NO7t+OgH9yXelOcxYWWSUi1+UT0wBgut+BND/q9zfDfZ7v98bdl/GqepaqHvHbv42q2lZVc9bZXwH0BvL6XmUjIs3c5XbkMTsSJ+EnA0NUNclv3gJgrIiMFpE2hbzNF0B9/2SG8/v3lSI24xxLagFPAe9LANXUfva7MV5YhHWKxJJEdo8Dd4tIdD7zl4lIIrAOJ6G8nmP+y+K0DawE9gD3u9PfBa6HrKLmxTgHy5zSgMZAjKqmqeocLbhzrWXuWf8REXnZf4aqzgPqisjpOAe1iQVsx/fe9YDW7lngUlVNEJHmQB/gIVVNVtUVOGeoOatdwPlBdBWRGPf1COBzVU3BOYvdpqpvq2q6qi7HKW0NzbkR9we5GKdqoQfO5/mLG8eZwEZVPVjI/vh7W1U3uNudTD4lQFcTv8/0SB6lnfeB4W7yGua+LswIYIyq7lPV/TgHA/+zy0zgCVVNyXEwyk8a0EFEarpns8sCWAfgO+AaEensJtvHcQ7WVd351YGjfssfBaqLOO0SqtpZVfP63gLUB/b6XohIV/fzSxCR39zJ9wDf45x0/BHoIk51U0Hfm6J6VlUPFfY5itMesBOYpaoz8likBnAW8G4ecdyNUyq/C1jrlq4H5fU+bhyf4P5e3KTSA/f3r6qfqOput0TyMU4JqFeA++rb/oPAFBFZEeh6RWFJwo+qrga+BnI1Qrm64/yQrsU5W6mWY/6fVbW2qjZV1RHuAQGcA8nl7tn0NcAcVd2Tx/b/iVM6+V5EtkgejWE543Hfr7aq/jmP+e/hfJHPx/khFuQ9YDowya0S+Yd7IGwCHFLVY37Lbgea5tyAu8w3uNVcOGeHviquGKC3/wEY52DQKJ94ZuPUPfd1n8/CaQs6z31dFHv9np/A+R/mZ7ffZ1pbVSf7z1TVHTj/o2dwktXOAN6/Cc5n5rPdneazX1WTA9iOz1U4ddrbRWS2iJwVyEruwfAJnOS8zX0cwyklAhzHaUvwqQkcL+RExecgzgmO771WuO1kf8A5K0dVn1HVn91F0lX176q6oJDvTVEF8v9AVR/AaTu8SETi8ljkgBvPuyJycY51k9x96YFzYjUZ+MQ9AczLu8BQEYnCOTmYrqr7AETkRhFZ4feb6ISTcAMiTtuL76SjW6DrFYUlidyewLliINdBEEAdk4H5OGdihXLrd+fj/GBuwK+KIsdyx1T1/1T1NJzGsPtF5IKi70KW94A7gWmqeqKQGNNU9SlV7YDTQHgZztnPbpwSSQ2/xVsA8fls6iOcM+2zcOrrf3Kn7wRm5zgAV1fVP+aznZxJYjaFJ4nS6tJ4Ik7VYF6ls7xi2I2TJH1akL3qL+c6iZw8u0dEsiVSVV2sqkNwqs++xDlIBURVX1PVNqraECdZhOFURYLTgOtfTdeFAhp1c5iJc8DNeeKUVwxPquq2HJPz+94UVcDfAVXdiLPveVZNqernOMeCT0Xk/HyWScA5YagGxObzVnOBQzgXDFyPW9Xklpz+i3MiV89NqqvJXiVYmIZAHeDLAJN5kVmSyEGdRumPcS6NK8hzwK05f8AFmIhTLDwD52qJXETkMhFp7Rbvj+LU2ebZqBYIVd2Kc1B9tLBlxWksPkOcRu4EnCqNTPdMeR7wrDiN252Bm8m/mmUazgFxDE77gC/+r4G2InKDiIS7j54i0j6f7cwDTscpei9Sp7EyBqcE93M+6/wOtJTiNfYXxcc4bUx5HZx/B+qJSC2/aR8Bj4lItNte9TgFV1OtBDq6VTZRwJO+GeJcjjlCRGqpahrO/yqg74j7/+skjhbABOAlVT3sLjIR58SkqTiXxv4f8E4g23bX3QN84b6H1409r7P0vOT3vQm2FJx2qjyp6kc4B/GvRKQPgIj8zf3uRrj7eA9wBPjNnT9L/C5rdw/eE4Hnca6SmurOqoaT1Pa7643mZPtQoML99iMoLEnkbQy5q5KyUdVfcQ5Wfwlwm1/g/Ai+KOCsvg0wA6fYPx94XVWLe0bli3Ouutd/F6IR8CnOQWcdztm6r8QzHGiJc/b7BU79eV71uLj1t5/jNGR+6Df9GM6BdZi7nb04P5rIfLaTCCzDuQIm1Z08H9juK6rn4RP370ERCbSevsjc6oYZedV7q+p6nKSwxa1CaAI8jXNhxCrgV5z9erqA7W/A+Q7OwKmjztk4ewOwTUQScBqnRwCISAtxrmRqkc+mo3D+J8eBRTif59/85o/HOYD9inNG+407DXf7aySfG9nc6rLzgbXuegk4B82eOFWsBcrve1MKMinkOKiq7+IkzG9EpBfOgf1tnCqp3TiNxpeq6nF3leY4bWj+JuKUID/2tXGo6lqcK8jm45xcnJHHelny+fx9V00GLalKkEooJg8ishnnqpQ8D7DGmNIlIr6keUdJlF7cK6Ymq+rZpxxcYO93DfCqqjYI1nsE7dpak52IXIVzBpLrckVjTMj8G+cqxf0ickaApe58qeounDa9oBORpTiXzj4U1PexkkTwidMVRgfgBlWdHuJwjKnwxLmz+9u85qnTK4IJkCUJY4wx+bKGa2OMMfmqUG0S9evX15YtW4Y6DGOMKVeWLl16QFXz7GmiQiWJli1bsmTJklCHYYwx5YqIbM9vnlU3GWOMyVfQk4SIDBSnW+hNefVFJCJ3iMivbv8lc8XtFVREWopIkjt9hYi8EexYjTHGZBfU6ia3i4fXcO5I3AUsFpEp7p2GPh+q6hvu8oNxrlse6M7brKpdgxmjMcaY/AW7JNEL2KSqW9yuFSbhdHKVxe0gy8fXl4kxxpgyINhJoinZu+7dRR69q4ozmtRm4B9k71gvVkSWu90hn5tzPXfd20RkiYgs2b9/f16LGGOMKaYy0XDtdl/cCuf2ct+IXXuAFqraDWfwng9FpGYe605Q1ThVjYuOzm+sIGOMMcUR7Etg43F6RPRpRv7jEIBTHfUfyOoV0tdb4lK3pNEWpzdNYyq01NRUNq9bxYk0QIoyvIAxflQJk3TqRzekcbMYPJ6ilwuCnSQW44xVHIuTHIbhjDecRUTauIN/AFyKO4C5OEOIHlLVDBE5Dacb7S1BjrdU7Nq1i6VLlzJkyJDCFzaV0uZ1q6i95h1O3/gRHk0PdTimnFLxkFqlIVs73s3uba3peub5hEfkO3xGnoJa3aSq6TgDdkzHGaNgsqquEZEx7pVMAHe5/aSvwKlWGulO7wuscqd/itOV76Fgxltahg4dyhVXXMHRo0cLX9hUSifSoKElCHOKRDOJPLGHtkufRKrUYdpHb5CellakbQT9jmtVnYYz6pT/tMf9nt+Tz3qf4QyvWOE0beq03c+bN49Bg/IcP91UdiKWIEyJ8WSmIh4vW9YuZ/f2jbRoneeIrXmvG8S4TD5uv/12ACZOzGuIZGOMCQ4RDyeOFa0Gw5JECNSvXx+AlStXhjgSY0xlk5lZtAH4LEmEwNtvvw1AtWoFDqNtjKngRn2ZxICJiaEOo0AVqhfY8uLIkSMA/Pbbb2RmZhbrsjRjyqpRXybx7so0/nJ2BP+4MCpr+q6ETJq/eJyfRlalX8sw5KmTnS1UCYOWtT3c0j2c+8+KzJq+7UgmsS8dz3pdMxLa1ffw13MiGdIuPOCY/N/LJ9ILyY85t161HHeM7UeVlwZG8ufekdmWu++7ZMYtTOWCWC8zbnRO7DJV+ff8VN5ZkcbWI5mEeyCmtofL24bxdP+oXO9VnlmSCAFfce/YsWPMnz+fPn36hDgiY0pWVBi8vDCVP/WMIKZ2/idBrw6K4qoOYSSlwfeb07nr22SqRwi39ch+meZXw6rQq6mXw0nK87+kctXkJObeJJzZLPBDmO+9fHLefdKilvC/ZWnZkkRyujJxVRoxtbIvPWZ2Ci8tTOWVQVGc1SyM5HRl9b5MFuzKCDie8sJOYUPAf8jYzz//PISRGBMcZzf30qWRh7/+mFzgcrWioFF1D7F1PNweF0Hnhh6mb859VVfdKkKj6h7aR3v57+VRRHjhq/XpzNqWjndMAjuPZq9nn7gylVrPJZCYevK35nsv36Nh9eyHv2Edw9lyOJOFu06+/6dr06gTBee1zJ6Mvlyfzs3dIri+cwSt6nro2MDLtZ3CeXFg9lLEuytS6fDacSL+nkCzfx/jsR+TSc/Mu3u6HzY7+7IrIfu+fLw6japjE0hIcdb7/Xgmo75MIvqfx6jxbAJ93krk5+3BuxLOShIh4N9wtHDhwhBGYsqTfu/krru+pmM4d/aM4ESacskHJ3LNH9U1nFFdIzhwIpOrJyflmv/HuAiu7RTOzqOZ3PBF7vmzRhWv3UyAf10YxXnvnOC+MzOIa+ItcHlV5cetGazbn0nbegWfu4Z5INwLaZnQr2UYbep6eGt5Gk/0O1kC+O+yNK7rFE61iMDvVq8RKQzrFM5/l6XR2y2hTFiaxi3dI1h/IPuBu3ENYfb2dOITMmlaM+94v9mQxk1Tknn6/Eiu6hDG8j2Z3PFNEgL8PY8qqQtO89K4uvDBqjQeOufkvry7Mo0r2oVRM1JISlPOf/cE7aM9fDuiKrWjhI9Xp3HheydYcXs12kcX/DkXh5UkQqB169ZZzzdt2hTCSIwJnnNjwhjSLowHvs+/NHHLlGSqP5NAxNPHGPDeCUTg3t753xGcnK48NTuFhBQYcJpzIL+tRzhvrUgl0y2hrz+QwdwdGbmqrHzv5Xv8fXZKru3f1iOCSavTOJairD+QwYJdGYzumrvt48WLo0hKg+YvHuf0V48z8sskPliVlq2U8NwvqVzVPoxHzo2kbT2npPHkeZH8a34qqRm5SxMeEa7vHM57q07e7Pb78Uy+35zOyC7Ovny8Jo2EFOXjq6sQ18RL67oeHu0bSZ/mXsYvLdpNcoGykkQIDB48mKeffprw8HAOHDgQ6nBMOVHQWX3VcClwfv2qngLnN69V8Pzien5AJB1fT2TKb2l0b5z7LHdsf6cBeu/xTB6ekcLQDmFZZ/H+LnrvBB6BpHSoEyW8eHEkA1s7y43sEs6jP6YwfVM6g9qE879lafRo7KFbjvfzvZdP3Sq5Sxm9mjoH3o9Wp7H+QCaXnx6Wq1oKoF19L7/+sRor9mYyd0cG83alc8vUJF5c4GHO6GpUCRfW7Mvg2o7ZSwzntQwjOT2FzYcy8zzrH9klnOd/SWXZngy6N/bywa9pNKgmDDjNWXZxfAZ7jyu1nzuWbb2UDKgSHpw+vixJhICvTaJZs2Zs3bqVbdu20bJly9AGZUwQtK3n5fYe4Tw0I4VvR1TNNb9hdaF1XQ+t63r4cpiHtq8cp1tjL31jsh+a3h5ShR5NvNSOchKev3pVPVzdwakmuuC0MCauTOPp/tmvUPJ/r8Lc1iOC/yxJZedR5YM/VMl3ORGhW2Mv3Rp7ubt3BHN3pHPu2yeYvCaNkV2L1j+ST/toL3FNPExc6STViSvTuL5zOF6PkwAyFdpHe/ji2tyfZdXAL/YqEqtuCoGXX34ZgIEDnQH4tm7dGspwjAmqJ86LZPexTCYsTS1wufpVPfypZwR//jY528UdAE1rOgf4nAnC5/Ye4UzdkM74JWkkpSvDOxX/iHl953A2HsykRiRc2CrwOv729Z3Y9iU6sXds4M3VoDx7WzpVwqBVAclqZJcIPlqdxrI9Gaz8PZMbu5zcl7gmXrYczqRmJFnJ1fdoUiM4h3NLEiGQmOg0QPoufV28eHEowzEmqKKreXi4TyTjFhScJADu6uU0Ek9aXbSrdc5pEcbp9Tw88EMywzqGUyOy+FUvNSOF+PtrsOqO6njy6ab9qskneGFeCvN3prP9SCbzdqZzwxfJhHvg0rZOKeiRcyL4bF06z81NYcPBDCavSePJ2Sn831kRRHjzj294pzAOJyk3T0mie2MPnRqcTFQjOocTW9vDpR+e4PvN6Ww74lyN9eycFL5cH5w2CUsSIeC7uqlevXrUq1eP9957L8QRGRNc950VQf2qhR+4G1b3cGOXcB6flZLvpaL5ubV7OKkZ5GqwLo5aUVJgohnYKozvNqfzh8lJtH31OEM/SSLCC7NHVaWD29ZwSZtw3hocxbsr0+j0eiL3TU/mzriIbFdh5aVeVQ+Xtg1jxd5MbuycvUQUFSbMHlWVuMZeRn+VRNtXjvOHyUks2p1BTK3gHM6tTSIEfEVpEaFKlSqsX78+xBEZU3LeuSJ3PX5UmLDjvhrZpukTuQaaBGDC5SfXb1lb8l0up/hjStdGHno2zV1FVNg2tt1bo8D5Offp1h4R3BpAMhrZNaLA9om8PisgzzYHn3pVPfznsirO6GylwEoSIeArSYgInTt3Jj093RKFMcV0NFlZHJ/BhKWp3HfmqZciTHaWJEKgU6dOgJMkLrroIgAmT54cypCMKbeGTDpB33cSubJ9ONd3DtIlPpWYJYkQuOyyywAnSVxzzTUA/PTTT6EMyZhya9aoaiQ9WpO3h1TJt6HZFJ8liRBr3LgxVapUYdeuXaEOxRhjcrEkEQLPPPMM4JQkwBnz+tixY7muDTfGmFCzJBECae5A5L4k0bNnT37//XcrTRhjypygJwkRGSgiv4nIJhF5OI/5d4jIryKyQkTmikgHv3mPuOv9JiIXBzvW0uJ/dRNATEwMAOPGjQtVSMYYk6egJgkR8QKvAYOADsBw/yTg+lBVz1DVrsA/gH+763YAhgEdgYHA6+72yr2c1Ur9+vUDYPbs2SGIxhhj8hfskkQvYJOqblHVVGASMMR/AVX1H1ewGuA7gg4BJqlqiqpuBTa52yv3cpYkatSoQY0aNdiwYUMowzLGmFyCnSSaAjv9Xu9yp2UjIn8Skc04JYk/F3Hd20RkiYgs2b9/f4kFHkw9e/YETiYJgLZt23Ls2DESEnKPxWuMKfsemZFMw38dQ55K4J0VhfdTVV6UiYZrVX1NVVsBDwGPFXHdCaoap6px0dHRwQmwhA0aNAjIniQGDBgAwIcffhiSmIwpKaO+TEKeSsj1mLQ6jSdnJSNPJdB9/PFc663cm5G1rG8Iz7k70pGnEth2JDPX8mXJwl3pPPdLKhMui2LP/1Xn2o7hpGcq//glhdNfPU7U0wm0eeU4ry3KnTw2HMzg4vcTqTo2gfr/OMYdXydlG3Y11ILdd1M80NzvdTN3Wn4mQVaXJEVdt9zIeXUTwMiRI3n++edZvXp1qMIypsSc28LL5KHZ+yWqHSU8NzeD6KrC+gOZWQPr+IxfmkpMLWH70eAcIFMztMDeV0/FxkOZeIRsgxo9OjOZCcvSmHBZFF0aeZm/M4PbvnY6AvT1+3Q8Vblg4gk6N/Qy7+ZqHEpSbvoqiSPJSUy6Ov/+m0pTsEsSi4E2IhIrIhE4DdFT/BcQkTZ+Ly8FNrrPpwDDRCRSRGKBNsCiIMdbKp566qlc09q3b0/79u1tbAlTIUR4oVF1T7ZHVJhzgK4ZiTNIkN/4EifSlA9+TePmbif7Xtp2JJNz33bG7Y596TjyVELWON+qyr/mpXDaS8eI+HsCrV4+xrgF2YcjbTnuGI/9mMyd3yRR7x/HOPdtZ115KoH3V2U/ox8wMZFRX54c4/ur9Wl0G3+cqmMTqP1cAr3+e5zlezLy3NdRXyZxwxfJZCpZJSFwxqb+v7MiuLJ9OKfV8TCiczi3dItg7JyTcX74axoHTigf/qEKXRt56R8bxmuXRPHxmnS2Hi4bpaegliRUNV1E7gKmA17gLVVdIyJjgCWqOgW4S0QGAGnAYWCku+4aEZkMrAXSgT+pat7/pXLGvxdYf+eeey6TJk0iLS2N8HDrg8acdO93yazYW/pf/66NvIwbGFX4gkV0W49wLv3wBC9cHEXVcGHS6jSa1PBwbszJkkXzmsJXw6owZFISi26pRvNaklUSeH1xGn/7KYWXBkZxfksvM7dmcO93ydSIEG7ufjLRvLwwlfvPimD+zVVJD/CYu/d4JkM/SeLp/pEM7RBOcrqyfG8mYfmcUr80MIpujTz83/cp7Lq/etb05HSIynGErRIO248q249kElPbwy87MzirmZdaUSePBRe1CsMj8MvOdGLrhL7DwqB3Fa6q04BpOaY97vf8ngLWHQuMDV50oZHz6iafiIgIEhIS+PLLLxk6dGgoQjOmRMzalkH1Z05ehNG0poff7jp5AD2nRRjNanr4eHUao7tFMGFpGrd2z35i5PVI1jjU0dWERn5jTT/3Swp394rIGjuiTT0vvx3IZOyclGxJomdTL0/2K1qS23NMScuEazqG07K28555jUftUytKsg7y/jEOahPGywtTuSA2jE4NPCyKz+Ct5U5V8+5jTpLYcyyTRtWzHwfCvc5+7zlWNtolbDyJEMivJHHttdfy6quv8sknn1iSMNkE42w+mHo38/Ku31gJeZ2F39rdGZe6RxMvK/Zm8PV1VVi9r/DT/YQUZVeC0jcm+4H7vJZeXlqYyok0pWq489vq1aTot1Z1bujh4lZeOr1+nAtbhdEvxssf2ofTvIiD+rw0MIo7vk6i6/hEBGhSQ7i5WzjP/ZKKpxz1Q2hJIgR8JYmc+vTpg9frZcGCBaUckTElq0qYMwZzQW7sEsEjM1O4f3oyV7YPc8evLtl6+GoRuY/GAuTsJi3N7229HuHbEVVZvDuTGVvS+WxdOg/PTOGToVW4rG3g1cB1qwiTh1YlNUPZl6g0qSG8scQpSZxWx/lsGtfwsPNo9n1Oy1AOJSmNa5SNTFImLoGtbPr27QvkLkmICM2aNSM+Pj7fRGJMRVG3inB1h3Bmbs3g1u5517372iAy/H4ONSOFZjWFn7dnb6OZvS2D2DqSVYrIT4Nqwm6/qpyUdGXt/uy/NxGhV1Mvfz03kp9HV+O8GC9vryjeGNIRXqFZTQ8eET5anUbfGC/R1ZxDb5/mXubvyiAh5WQ8P2xJJ1OhT/OycQ5vSSIEfAMN5UwS4DReZ2Zm8v3335d2WMaUuv9eHsX+v1Snf2zeB8SYWoJHYNrGdPYlZnI02TmYPnJOJK8sSuW/S1PZeDCD8UtS+c+SVP56TsHjRwMMOC2MN5amMn9nOqv3ZTDqqyRSM04epOftTOfvs1NYuCudHUczmbklnVW/Z9Kh/snDZbtXj/NqHvc8+Fscn8Ena9LYfCiT+TvTuXryCVbszeBlv6rD684Ip35V4brPkli5N4Oftqbzp2nJXNsxjNg6ZePwXDZSVSWTmOheipdHkrj99tt5//33Wbt2LQMHDizt0IwpVVFhknVpbF4aVvfw7AWRPPdLCvdOT+bcFl5mjarGH+PCSUxVnpmbwp3TlOY1hecGRGZrtM7Pvy6K5NapysXvn6BWlPDXcyLZn3gySdSKFObvyuC1xakcTlYaVRdGnBHO3847mYB+O5jJgRMFl/ZTMpSnZqew+XAmEV7oGxPGvJuqcUbDk+0k1SOEGTdU5e5vkznrzUSqhAtXtw/j3xeXnTYoqUhjGMTFxemSJUtCHUahunTpwqpVq1i4cCG9emXvjkpViY2NpUePHnz22WchitCE2tKlS+kxtX+owzAVyNLLf2TOJ69z0dCb6NCjT7Z5IrJUVePyWq9slGcqmYISs4hw5pln8u2332bdmW2MMaFiSSIE8rtPwqdevXokJSUxadKk0gzLGGNysSQRAvndJ+Hzxz/+EYAPPvig1GIyxpi8WJIIgcKSRKdOnYiKimLhwoWlGZYxxuRiSSIEfJfAFqRTp04cOXKE8jJGhilhqqjYz9OUjEwJAy3evVf2LQyB/v2dq1byK0kAXHnllQBMnDixVGIyZUuYpJNapWGowzAVxInap5OZXLwBzSxJhMDBgweBgpPEnXfeSVRUFDt37sx3GVNx1Y9uyNaOd5PpCX0voKb8ypQwjtfpyIaeT7N3xxYEweMtWn9WdjNdCIwd63RsW1CSqF27Nueffz5ff/01L774YoHLmoqncbMYdm2OZdmgbxBP0TupMwYAzSQzOYG92zZzdN8uMjWDGrXqFmkTliRCoLCGa5/OnTvz7bff8uOPP3LBBReURmimjPB4PLTr0pNPxz9P4rGjVKlWA/FYwd8UT0Z6GokJR+jRdxCNY1oXaV1LEiEQ6F3ul1xyCc8//zyvvPKKJYlKqFbdaK6+/SGW/jydQ/t2k5lZIcbcMiEQVaUazVt3oMtZ/fEU8WTDkkQIFHYznU/fvn2JjIxk9uzZpRGWKYNq1Y2m/xXXhzoMU4lZ+TWEAmln6N69O0eOHGHbtm3BD8gYY3KwJBECQ4YMAQJLEiNHjgTg3//+d1BjMsaYvFiSCIFzzjkn4GVHjx5NeHg4v/76axAjMsaYvAU9SYjIQBH5TUQ2icjDecy/X0TWisgqEZkpIjF+8zJEZIX7mBLsWEvLnj17gMBKEhEREYwYMYLly5eTkpIS7NCMMSaboCYJEfECrwGDgA7AcBHpkGOx5UCcqnYGPgX+4TcvSVW7uo/BwYy1ND333HNAYEkC4Oqrr+bo0aN88sknwQzLGGNyCXZJohewSVW3qGoqMAkY4r+Aqv6kqifclwuAZkGOqcwINEn069cPgKeffjqI0RhjTG7BThJNAf9+JXa50/JzM/Ct3+soEVkiIgtE5Iq8VhCR29xllpSXzvACvQTWp1q1asTExLBhwwaSkpKCGZoxxmRTZhquReR6IA74p9/kGHdIveuAcSLSKud6qjpBVeNUNS46OrqUoj01xRkydsSIEagq48aNK/mAjDEmH8FOEvFAc7/Xzdxp2YjIAOBRYLCqZrXOqmq8+3cLMAvoFsxgS0ug3XL4e/DBBwF45513ghGSMcbkKdhJYjHQRkRiRSQCGAZku0pJRLoB43ESxD6/6XVEJNJ9Xh/oA6wNcrylYvjw4UDRkkStWrWIjY1l48aNJCYmBis0Y4zJJqhJQlXTgbuA6cA6YLKqrhGRMSLiu1rpn0B14JMcl7q2B5aIyErgJ+A5Va0QSaJXr15A0ZIEwAMPPICqWjcdxphSE/S+m1R1GjAtx7TH/Z4PyGe9ecAZwY0uNLZs2VKs9W655RYef/xx3nvvPS655JISjsoYY3IrMw3XlYmvi42iliQiIiLo378/n3zySdYNecYYE0yWJEKgOA3XPv379ycjI4O//e1vJR2WMcbkUqwkISIVom0gVE4lSdx2222Eh4fz2WeflXRYxhiTS6FtEiLSALgUqA14gc5AVRH5D7ARSAS+9l2uagp3KknC4/Fw7rnn8uOPPzJ37twidRZojDFFFUhJ4muc7jVqAlVw7lfoinO1Ug2cpPF+cMKrmIpzM52/MWPGAFiVkzEm6AK5uilZVf+Yx/SXfU9ExK7JLIKbb76Zl156qVglCYA+ffrQuHFjVq1aRUZGBl6vt4QjNMYYR6ElCVXtG8Ay55VMOJVDly5dgOJVN/mMGzeOQ4cO8cMPP5RUWMYYk4td3RQCa9c67f6nkiSuuOIKoqOjefbZZ0sqLGOMyeWUkoSIPFlCcVQqr7766ilvIyIigtatW/Pzzz+zbNmyEojKGGNyO9WSxNISiaKSOZWrm/z5GrAfeOCBU47JGGPyckpJQlWnllQglUlJJYkBAwZQr149Zs+ebZ3+GWOCotAkISKviMjL+T1KI8iKpqSSBMCdd95JZmYmjz766ClvyxhjcpLCrtkXkZEFzVfVd0s0olMQFxenS5YsCXUYhQoPDyc9PZ34+HiaNGlySttKS0ujatWq1K5dm/IyMp8xpmwRkaXuAG+5BHIJ7LsFPUo+3IrvzjvvLLFthYeH89hjj3HgwAEWLlxYYts1xhgoQpuEiESLyL9EZJqI/Oh7BDO4iqpDhw5AyVQ3Adx///3UrFmTF198sUS2Z4wxPkVpuP4ApyuOWOApYBvOyHOmiHyXrJZUkqhRowYDBgzg448/Zu7cuSWyTWOMgaIliXqq+iaQpqqzVfUmoH+Q4qrQJkyYAJRckgB46KGHAPjTn/5UYts0xpiiJIk09+8eEbnUHZu6bhBiMsXQq1cv2rRpw6pVq7Lu6DbGmFMVyCWw4e7Tp0WkFvB/wAPA/4D7ghhbheR/NVlJliTA6c8J4Pbbby/R7RpjKq9AShLxIvI/IAlIUNXVqnq+qvZQ1SlBjq/CyczMzHpe0knikksuoWnTpsydO5f4eBvewxhz6gJJEu1xGqgfA3aKyEsicmagbyAiA0XkNxHZJCIP5zH/fhFZKyKrRGSmiMT4zRspIhvdR4H3a5QXwUwSAK+88goAb7/9dolv2xhT+QRyn8RBVR2vqufjDD60BXhRRDaLyNiC1hURL/AaMAjoAAwXkQ45FlsOxKlqZ+BT4B/uunWBJ4De7vs+ISJ1irR3ZZDX6+Xee+8FgpMkrrzySgYPHsy//vUvDh06VOLbN8ZULkXqu0lVdwNvAv8BjgG3FLJKL2CTqm5R1VRgEjAkxzZ/UtUT7ssFQDP3+cXAD6p6SFUPAz8AA4sSb1nk8XiIjY0N6nuMGTOGo0ePMmTIkMIXNsaYAgSUJEQkSkSGisjnwCacS18fBgrrU6IpsNPv9S53Wn5uBr4tyroicpuILBGRJeWhW4r09HQWLVoEBKckAc6gRi1atGDu3LmsXr06KO9hjKkcArm66UNgB3ANzg11LVV1lKp+p6oZJRWIiFwPxAH/LMp6qjpBVeNUNS46OrqkwgmalJQUPvjgAyB4SQJg/PjxAIwYMSJo72GMqfgCKUl8B7RS1aGq+pmqJhdh+/FAc7/Xzdxp2YjIAOBRYLCqphRl3fIm2A3XPgMHDqR9+/asWrWKn376KWjvY4yp2AJpuJ6oqsdyTheRISLSu5DVFwNtRCRWRCKAYUC2y2bdm/LG4ySIfX6zpgMXiUgdt8H6IndauVZYr7slyVdiGTVqVKm9pzGmYjmVQYd6A4+JyLf5LaCq6cBdOAf3dcBkVV0jImNEZLC72D+B6sAnIrJCRKa46x4C/o6TaBYDY9xp5VpplSQAunXrxjXXXMOOHTusTydjTLEUOp5EeVIexpM4dOgQ9erVA+Do0aPUrFkzqO+XmJhIu3btiI6OZtGiRYSFhQX1/Ywx5c8pjSchIj1FpJHf6xtF5Ct3ZDrru6mIatasyf333w8EvyQBUK1aNZ544gmWL1/OrbfeGvT3M8ZULIFUN40HUgFEpC/wHDAROApMCF5oFVNYWBiNGzcGSidJAIwePZoaNWrw7rvvsn379lJ5T2NMxRBIkvD6tQVcC0xwr3L6G9A6eKFVTCdOnCj19gGv18v48eNRVQYPHlz4CsYY4wooSYiIryL7AsB/NDqr4C6ihIQEvvrqK6D0ShIAw4cP54wzzmDVqlW8+66NOmuMCUwgSeIjYLaIfIXTE+wcABFpjVPlZIqgNK9uyunrr7/G4/Fw9913k5aWVvgKxphKL5D7JMbijCHxDnCOnrwcyoMzroQpgmCOJ1GYFi1a8MILL3Ds2DH++c8i3dhujKmkArm66XFVXaCqX6hqot+s34FngxdaxeRfkgiFe++9l6FDhzJmzBjmzJkT0liMMWVfINVN5+TsElxEGgKzyd4+YQIQyuomn5deeomMjAwuueQSkpOL0suKMaayCSRJDAa6iMi/AUSkDfAL8IaqjglmcBVR06ZNeeABp5YuVEmicePG3HHHHRw/fpzLLrssJDEYY8qHQNokkoErgZYi8hEwA/iLqr4R7OAqorCwMOrUccZOClWSAHj55Zdp0aIFM2fOZOLEiSGLwxhTtgXSJnE/cDewEKeTveVArDvs6P1Bjq/COXz4MLNmzQJCmyREhNmzZ+P1ern11ls5cOBAyGIxxpRdgVQ31XAfUcDLOEmiht/DFMGBAwf44YcfQh0GAC1btuS5554jNTWVP//5z6EOxxhTBhV6M5yqPlUagVQWZaHh2t8DDzzAoUOHePbZZxk4cCA33nhjqEMyxpQhAXXwVxLLGEdZSxLgjIl9zjnncNNNN/HZZ5+FOhxjTBkSSHXT2wXNFBEv8N+SCafiC+XNdPkJCwvj9ddfR1UZPnw48fHlfgBAY0wJCSRJzBOReBFZ5Q4K9LWIDBKR2SKyCtgGzAtumBVHqG+my88ZZ5zBuHHjSEtLo1evXmRklNjw5caYciyQS2BvA9oDl+FcCvsW8CkwFrgc6K6qdwYzyIqkXbt2PPjgg0DZKUn43H333QwePJjdu3fTr1+/UIdjjCkDAhq+VFUTVHWHqm5V1c+BV1T1e1Xdrqr7gxxjhRIWFkaVKlWAspckAL744gvatm3L3LlzmTDBhgsxprIr1hjXqvpwzmki0vHUw6n4du/ezYwZM0IdRr48Hg9Lly6ld+/e/PnPf2bhwoWhDskYE0LFShL5eC+viSIyUER+E5FNIpJXcukrIstEJF1Ers4xL8NtB1khIlNKMNaQ2bt3L7/88kuZLEX4VK9ena+//pomTZowYMAAZs+eHeqQjDEhUpJJItdRz73y6TVgENABGC4iHXIstgMYBXyYxzaTVLWr+6gQQ6qV1YbrnOrXr8///vc/EhMTufDCC1m9enWoQzLGhEBJJgnNY1ovYJOqblHVVGASMCTbSqrbVHUVUD6OnqfIlyTKcknCp3///rz++utZVzzt2LEj1CEZY0pZSSaJvDQFdvq93uVOC1SUiCwRkQUickWJRhYi/vdJlAd33HEHTz31FElJSXTu3JlDhw4VvpIxpsIIKEmIo3khi6WWQDw5xahqHHAdME5EWuUR221uIlmyf3/Zv9CqPJUkfB5//HHuuusujh49ymWXXUZKSkqoQzLGlJJAL4FVYFohy5yZx+R4wD+5NHOnBURV492/W4BZQLc8lpmgqnGqGhcdHR3opkOmd+/e/OUvfylXSQLglVde4dlnn2X+/PkMHTrUBisyppIoSnXTsmL00bQYaCMisSISAQwDArpKSUTqiEik+7w+0AdYW8T3L3M8Hg9hYWHlLkkAPPzww7z22mtMnTqV2NhYEhISQh2SMSbIipIkegPzRWSz20XHr263HPlS1XTgLmA6sA6YrKprRGSMiAwGp3NAEdkFDAXGi8gad/X2wBIRWQn8BDynquU+SWzatInp06eHOoxiu/POOxk2bBh79+6ldevWHD58ONQhGWOCqNCuwv1cXJw3UNVp5KiqUtXH/Z4vxqmGyrnePOCM4rxnWRYfH8+yZcsIDw8PdSjF9tFHHxEWFsb7779Pq1atWLt2LY0aNQp1WMaYIAg4Sajq9mAGUlmUx4brvLz33ntUq1aN8ePHZyWKmJiYUIdljClhAVc3iciTQYyj0igvN9MF4o033uBvf/sbqampXHTRRWzbti3UIRljSlgggw55RORNILIU4qnwfPdJlPeShM+YMWOYNWsW+/bto2fPnrzzzjuhDskYU4ICKUl8DRxS1UeCHUxlISIVJkkA9OnTh7lz53Ls2DFGjx7NfffdF+qQjDElJJAk0QP4PNiBVBYDBgzg/vvvx+v1hjqUEtWxY0eWLFlCjRo1GDduHP369bOBi4ypAAJJEucDE0Skd7CDqSzKW9ccgerUqRM7d+6kVatWzJ49m5iYGPbt2xfqsIwxpyCQkenW4lz++s/gh1PxrVy5ku+++y7UYQRNrVq12LBhA5dddhnx8fH069ePtWvL/e0txlRagXbLsRu4NMixVAo7duyo8AdNj8fD1KlTmTFjBgcPHiQuLo6777471GEZY4qhKHdctwlaFJVIRa1qyssFF1zA8uXLqVOnDq+++iodOnSwO7SNKWeKkiReEJF1IvJ3EekUtIgquIpyM12gmjRpwubNm+nVqxfr1q2jcePGfPTRR6EOyxgToICThKqej9OIvR+nj6VfReSxoEVWQVWmkoRPVFQUCxcuZOzYsaSlpXHdddcxatQou/rJmHKgSIMOqepeVX0ZuANYATxe8Bomp7CwMCIjI/F4gj3eU9nz17/+lfXr19OqVSveffdd+vfvz8aNG0MdljGmAEXplqO9iDwpIr8CrwDzyKNjPlOwyy+/nNtuu63C3ScRqDZt2rBx40beeecdVqxYQbt27bjiiitsICNjyqiinM6+BRwGLlbVfqr6H1W1i+CLQVUrTZtEXkSEkSNH8ssvv1C/fn2++uor6taty6RJk0IdmjEmh0CHL+2sqmcBP7qXw5pimjt3LtOmTatQHf0VV6dOndizZw/33HMPycnJDB8+nF69enHgwIFQh2aMcQVakrhJRNoANwczmMpgx44dbNmyJdRhlBkej4dx48axadMmTj/9dBYvXky3bt34+OOPK2UjvzFlTSC9wD7hLrcA8IiINVafgorWC2xJiY2NZf369cycOZPo6GiGDRtGdHS0VUEZE2KBdMvxFDAD+Bj4QVXHBD2qCqyy3SdRVP3792fx4sU8+OCDHDp0iOHDh9OuXTtWrFgR6tCMqZQCrW7qpap3Ar2CGUxlYCWJwnm9Xp5//nm2bNlC9+7d+e233+jWrRt9+/blyJEjoQ7PmEolkOqmIcBeAFX9m4gsFJEt7uPqoEdYwVStWpXq1atbkghAy5YtWbp0KTNmzKBZs2bMmTOHVq1a8cILL5CYmBjq8IypFAIpSTwIfOX3OhLoCfQD/ljYyiIyUER+E5FNIvJwHvP7isgyEUnPmXREZKSIbHQfIwOItcy7+uqrGTFiRKW9T6I4LrjgAnbu3MnChQvp2bMnDzzwALVq1eKaa64hISEh1OEZU6EFkiQiVHWn3+u5qnpQVXcA1QpaUUS8wGvAIKADMFxEOuRYbAcwCvgwx7p1gSeA3jjVXE+ISJ0A4i3zKvt9EsXVq1cvvvvuO958802qVavGJ598Qp06dbjyyivtslljgiSQJJHtwKyqd/m9jC5k3V7AJlXdoqqpwCRgSI7tbVPVVUDOGwcuxmkoP6Sqh4EfgIEBxFumTZs2jW+++cYu7zwFN910E4cPH+b555+nRo0afPnllzRs2JCxY8daycKYEhZIklgoIrfmnCgitwOLClm3KeBfCtnlTgtEQOuKyG0iskREluzfvz/ATYfOzp07iY+PD3UY5Z7H4+HBBx/kyJEjvPbaa7Rr147HHnuMZs2accEFF/DLL7+EOkRjKoRAksR9wGgR+UlEXnAfs3CqiO4NYmwBUdUJqhqnqnHR0YUVbELPrm4qeXfeeSdr1qxhyZIlnHfeefz444+cc845NG/enPHjx1upzZhTEMh9EvtU9Wzg78A29zFGVc9S1d8LWT0eaO73upk7LRCnsm6ZZfdJBE+PHj2YOnUq8+bNo1evXsTHx3PHHXdQo0YNXnjhBetE0JhiKMp4Ej+q6ivu48cAV1sMtBGRWBGJAIYBUwJcdzpwkYjUcRusL3KnlWuWJILvrLPOYuHChezZs4drr72WjIwMHnjgAZo2bcrw4cOZMWNGqEM0ptwI6qAGqpoO3IVzcF8HTFbVNSIyRkQGA4hITxHZBQzFGcxojbvuIZzSy2L3McadVq7VrVuX2rVrW5IoBQ0bNmTSpEkkJibyww8/0K9fPyZNmsSFF15I3bp1ueeee6yh25hCSEWqr42Li9MlS5aEOoxCjR49mpkzZ7Jjx45Qh1LpLF++nIceeoiffvqJ9PR0RIS+ffvyj3/8g549e1ryNpWSiCxV1bi85lW+4dHKCDsYhUa3bt34/vvvOXHiBM888wzNmzdn3rx59O7dmzZt2nDxxRczb968UIdpTJlhSaKUffTRR3zzzTehDqPSCw8P55FHHmH79u3s37+f//73v4SFhfH999/Tp08fateuzYgRI1i/fn2oQzUmpCxJlLL4+Hj2799vJYkypFatWtxyyy2sX7+e2bNn079/f5KSkvjwww9p37495513HuPHj2fPnj2hDtWYUmdJopRVpDagiqhv377MnDmTpKQkJk2axCWXXMLevXu54447aNKkCXXq1GHYsGEsXbo01KEaUyosSZQyuwS2fPB4PFx77bV88803rF+/nhUrVtC3b19OnDjBxx9/TFxcHNWrV+fqq69m+fLllvxNhWVJopRZkih/RIQuXbowe/ZskpOTmTJlChdddBEAn3/+Od27d6dZs2Z07NiRJ598kkOHyv2V2sZksSRRypo2bUr9+vUtSZRTIsLll1/O9OnTOX78OHv27OGtt96idevWrF27lqeeeop69eoRHR3NH/7wB37++WcrZZhyze6TCIERI0awaNEiNm7cGOpQTAk6evQor776Kp9++ilr1qwhLS0NgMaNG9O1a1datGjBHXfcQdeuXUMbqDE5FHSfhCWJELjuuutYvHixJYkK7pdffmHdunXMmDGDL7/8MqvvqMjISNq3b8+gQYO47777KA8dU5qKzZJEGfKf//yHRx55hAYNGrBhw4ZQh2NKSVpaGpMnT+aDDz5g4cKF2dotOnbsSHR0NF27dmX06NF07tw5hJGaysjuuC5D9u7dy9GjR61NopIJDw9nxIgRTJs2jYMHD3L06FE+/PBDxo4dS9OmTZk9ezbjxo2jS5cuREREcPrpp3Pvvfeyfv36rIsdjAmFsFAHUNnYeBIGoGbNmgwfPhyAv/71r1mX1n7++edZVZEbNmzgpZdeonbt2oSHh9OxY0cGDBjAddddR2xsbIj3wFQWVpIoZZmZmZYgTC5Vq1Zl9OjRTJ06lb1795KWlsbSpUt56623GDhwIAkJCcyaNYvHHnuM0047jcjISC688ELeeOMN5s2bZ73ZmqCxkkQps/skTCC8Xi/du3ene/fujB49GoBdu3YxadIkvv/+e1atWsXixYuzjY1RpUoVmjdvTrdu3RgwYABXX301tWvXDtEemIrCGq5L2VtvvcWjjz5KnTp1WLt2bajDMeWYqrJ9+3a++OILJk+ezIYNGzh8+HBWlabX66VDhw40adKEqKgo+vbty+DBg2ndunWIIzdljV3dVMYMHTqUtWvXsmbNmlCHYiqYzMxMFixYwNSpU/F4PCxfvpzZs2dz4sSJrGU8Hg9169blhhtuoGvXrjRp0oS4uDgrdVRiBSUJq24KAVW16iYTFB6Ph7PPPpuzzz47a5qqsmbNGqZOncrcuXNZv349+/fv5/XXX8827ndERAQNGjSgTZs29O7dmxtuuIHWrVsTERERil0xZYSVJErZ2LFjef7552nRogWrV68OdTimEktPT2fTpk1MnDiRuXPnsmnTJvbv3096enrWMl6vlypVqlC9enViYmLo1KkTZ555JhdeeCExMTEhjN6UJCtJlCEHDx4kKSnJShIm5MLCwmjXrh3PPPNMtun79+9n0aJFHDlyhHXr1jFx4kR+//139u7dy8KFC3nzzTcBp7uRdu3acfDgQVq0aEHXrl05++yzOeecc6hRo0YodskEgSWJUmZXN5myLjo6mksvvTTr9dNPPw1AcnIyP//8Mz///DMHDx4kOTmZNWvW8Ouvv7Jq1Sq+/vrrrHVq1apFjx49aNmyJXv37qVTp0707NmTc889l4YNG5b6PpniC3qSEJGBwEuAF/ifqj6XY34kMBHoARwErlXVbSLSElgH/OYuukBV7wh2vMHma4+wJGHKm6ioKC666KKsbtL9bd68mVmzZrFo0SLWr19PeHg4iYmJfPHFFxw+fJhp06ZlLevxeGjbti29e/cmOjqao0eP0rlzZ+Li4ujWrRuRkZGluVumEEFNEiLiBV4DLgR2AYtFZIqq+l/7eTNwWFVbi8gw4HngWnfeZlXtGswYS5t1sWAqolatWtGqVStuvvnmXPO2bdvG3LlzWbJkCWvXrmXbtm1EREQwc+ZM4uPjc3WlHhYWRqdOnejSpQs1a9YkMTGRjh070qNHD+Li4qhWrVpp7ZYh+CWJXsAmVd0CICKTgCGAf5IYAjzpPv8UeFUq8Gl2ly5daNiwoZUkTKXRsmVLWrZsyfXXX59r3vHjx5k3b15WAtm6dSu7d+/G4/EwY8YMdu/enSuJeDweOnfuTNu2bQkLCyMlJYVWrVrRrl07unbtSqdOnQgPDy+t3avwgp0kmgI7/V7vAnrnt4yqpovIUaCeOy9WRJYDCcBjqjon5xuIyG3AbQAtWrQo2eiD4JZbbmHq1Kns2LEj1KEYE3LVq1fPtwoLnCSyePFili1bxpo1a9i0aRN79uyhevXqLF++nK1bt2a7GsunadOmWceD5ORkmjZtymmnncbpp5/OGWecQe/eve3S3gCV5YbrPUALVT0oIj2AL0Wko6pm66RGVScAE8C5BDYEcRaZ3SdhTGCqV6/O+eefz/nnn5/n/IyMDH777TdWrlzJmjVr2Lx5Mzt37qRVq1bs2rWL5cuXc/jwYZYvX55r3QYNGlCzZk1at25NkyZNaNy4MY0aNaJFixYMHjwYgNTU1EqfTIKdJOKB5n6vm7nT8lpml4iEAbWAg+qUMVMAVHWpiGwG2gJl+0aIQtx///389NNPtG3bNtShGFPu+boe6dChQ77L+LovWb58OWvWrGHXrl00adKEGTNmMGfOHKpWrcqqVav4/fffycjIoFWrVllJ4pJLLmHBggU0atQo69GtWzceffRRAObPn094eDiNGjWiQYMGFTKhBDtJLAbaiEgsTjIYBlyXY5kpwEhgPnA18KOqqohEA4dUNUNETgPaAFuCHG/QHTt2LGtYS2NM8IlIVrvIlVdemTU9OjqaOXPmMH36dBo1akRmZiYHDx7M1qPu9ddfT5cuXdi7dy979+5l3bp12aq3brrpJtavX5/1unbt2gwZMoR33nkHgEceeQQRITo6OuvRqlWrctV/VlCThNvGcBcwHecS2LdUdY2IjAGWqOoU4E3gPRHZBBzCSSQAfYExIpIGZAJ3qOqh3O9Svth9EsaUDb6zft9Jm8fjyTqQ+4waNarAbXz44Yfs2rWLvXv3smfPHvbv35+tluCzzz7L1W4ycuRI3nnnHVSVBg0aUKNGjWxJ5LLLLuOqq64iIyODb7/9lnr16mU9ateujdfrLcFPoXBBb5NQ1WnAtBzTHvd7ngwMzWO9z4DPgh1fabP7JIwpG3xJIjU1tdjb6NatG926dct3/oYNG1BVjhw5wv79+9m/f39WR4rp6emMGDGC/fv3c+DAAXbv3s3KlStp2bIlV111FUeOHOHyyy/Ptj0R4bnnnuPBBx9k37593HTTTZx33nn85S9/KfY+FKYsN1xXSFaSMKZs8F0meypJIhAiQp06dahTp062UkZ4eDjjxo3Ld72aNWuycOFCDh48mO3Ru7dzgWhiYiK7d+9m586d+W6jJFiSKGV9+vRh9uzZoQ7DmEovZ3VTWRMeHk6vXr3ynR8bG8uyZcuCHocNX1rKbr31Vk4//XQrSRgTYiVR3VQZWJIIEUsSxoRWaVU3lXeWJErZjTfeyNy5cy1JGBNiZb26qaywJFHKUlNT7Y5rY8oAq24KjCWJUma9wBpTNlh1U2AsSZQyu0/CmLLBqpsCY0milNl9EsaUDVbdFBhLEqVs0KBBNp6EMWWAVTcFxpJEKbvllluIiYkJdRjGVHpW3RQYSxKlLC0tjczMTCtJGBNiVt0UGEsSpezKK69k6dKlliSMCTGrbgqMJYlS5huv15KEMaFl1U2BsSRRynxVTZYkjAktq24KjCWJUuYrSRhjQsuqmwJjSaKU2X0SxpQNXq8Xj8dj1U2FsCRRyoYNG0aDBg0sSRhTBkRERFhJohCWJErZTTfdRKNGjSxJGFMGhIeHW5IohCWJUnb06NFsg6IbY0LHShKFC3qSEJGBIvKbiGwSkYfzmB8pIh+78xeKSEu/eY+4038TkYuDHWtpGDx4MOvWrbOShDFlQEREhLVJFCKoSUJEvMBrwCCgAzBcRDrkWOxm4LCqtgZeBJ531+0ADAM6AgOB193tlWt2n4QxZYdVNxUuLMjb7wVsUtUtACIyCRgCrPVbZgjwpPv8U+BVcY6gQ4BJqpoCbBWRTe725pd0kBs2bKB///4cOHAg17yYmBjCwsI4fPgwhw4dyjU/NjYWj8fDwYMHOXLkSK75p512GiLCgQMHOHr0KKmpqXafhDFlREREBF988QWLFy8OdSinrH///rzyyislvt1gJ4mmwE6/17uA3vkto6rpInIUqOdOX5Bj3aY530BEbgNuA2jRokWxgoyMjCQ2NjbPeaeffjqRkZHEx8fj8eQueLVr146wsDB27NhBWFjuj7N9+/Z4PB62bt2adV128+bNuf3224sVqzGm5DzwwAPMmDEj1GGUiObNmwdlu8FOEkGnqhOACQBxcXHFulMtJiaGOXPmlGhcxpiy7/bbb7cTtkIEu+E6HvBPb83caXkuIyJhQC3gYIDrGmOMCaJgJ4nFQBsRiRWRCJyG6Ck5lpkCjHSfXw38qE7r7hRgmHv1UyzQBlgU5HiNMcb4CWp1k9vGcBcwHfACb6nqGhEZAyxR1SnAm8B7bsP0IZxEgrvcZJxG7nTgT6qaEcx4jTHGZCcVqcO5uLg4XbJkSajDMMaYckVElqpqXF7z7I5rY4wx+bIkYYwxJl+WJIwxxuTLkoQxxph8VaiGaxHZD2zPY1Z9IHefGxVfZd1vsH23fa98TmXfY1Q1Oq8ZFSpJ5EdEluTXcl+RVdb9Btt32/fKJ1j7btVNxhhj8mVJwhhjTL4qS5KYEOoAQqSy7jfYvldWtu8lrFK0SRhjjCmeylKSMMYYUwyWJIwxxuSrQicJERkoIr+JyCYReTjU8QSTiLwlIvtEZLXftLoi8oOIbHT/1glljMEiIs1F5CcRWSsia0TkHnd6hd9/EYkSkUUistLd96fc6bEistD97n/sdtVf4YiIV0SWi8jX7uvKst/bRORXEVkhIkvcaUH5vlfYJCEiXuA1YBDQARguIh1CG1VQvQMMzDHtYWCmqrYBZrqvK6J04P9UtQNwJvAn939dGfY/Beivql2ArsBAETkTeB54UVVbA4eBm0MXYlDdA6zze11Z9hvgfFXt6ndvRFC+7xU2SQC9gE2qukVVU4FJwJAQxxQ0qvozzngc/oYA77rP3wWuKM2YSouq7lHVZe7zYzgHjaZUgv1Xx3H3Zbj7UKA/8Kk7vULuu4g0Ay4F/ue+FirBfhcgKN/3ipwkmgI7/V7vcqdVJg1VdY/7fC/QMJTBlAYRaQl0AxZSSfbfrXJZAewDfgA2A0dUNd1dpKJ+98cBDwKZ7ut6VI79BudE4HsRWSoit7nTgvJ9D+rIdKbsUFUVkQp9vbOIVAc+A+5V1QTnxNJRkfffHbGxq4jUBr4A2oU2ouATkcuAfaq6VET6hTicUDhHVeNFpAHwg4is959Zkt/3ilySiAea+71u5k6rTH4XkcYA7t99IY4naEQkHCdBfKCqn7uTK83+A6jqEeAn4Cygtoj4TgIr4ne/DzBYRLbhVCX3B16i4u83AKoa7/7dh3Ni0Isgfd8rcpJYDLRxr3aIwBk7e0qIYyptU4CR7vORwFchjCVo3LroN4F1qvpvv1kVfv9FJNotQSAiVYALcdpkfgKudhercPuuqo+oajNVbYnz2/5RVUdQwfcbQESqiUgN33PgImA1Qfq+V+g7rkXkEpx6Sy/wlqqODW1EwSMiHwH9cLoL/h14AvgSmAy0wOlC/RpVzdm4Xe6JyDnAHOBXTtZP/xWnXaJC77+IdMZppPTinPRNVtUxInIazhl2XWA5cL2qpoQu0uBxq5seUNXLKsN+u/v4hfsyDPhQVceKSD2C8H2v0EnCGGPMqanI1U3GGGNOkSUJY4wx+bIkYYwxJl+WJIwxxuTLkoQxxph8WZIwJshE5GERGSEiT4qIikhrv3n3utNKfAB7Y0qCJQljgkQcHuBi4Ht38q84N3/5DAXWlHZsxgTKkoQxJUhEWrpjmEzEuQu2ORChqvvdRb7E7Y1YRFoBR4EDoYjVmEBYkjCm5LUBXlfVjkAPnL79fRKAnSLSCadE8XEI4jMmYJYkjCl521V1gft8IPBtjvmTcBLEFZzsXsGYMsmShDElL9HveS9gUY75XwM3ADtUNaHUojKmGGw8CWOCREQ6Auvd8R6yqOoJEXkI2BCayIwJnCUJY4JnEPBdXjNUdVIpx2JMsVgvsMYEiYj8ANzoN6SkMeWOJQljjDH5soZrY4wx+bIkYYwxJl+WJIwxxuTLkoQxxph8WZIwxhiTL0sSxhhj8vX/Lec7J3641V8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZwAAAEWCAYAAABSaiGHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA110lEQVR4nO3dd5gV9dn/8fdnKyAdEaQINiKgiLpifOwdE40+iUajMZhoNMVETfJLTHkimmaqJSZPYmKPsSR5VGLHGruCgkpsiBIQkSKytO3374+ZhbPLLnuW3XMWDp/Xde21Z2a+Z+aes7Nzz7ecGUUEZmZmuVbU1QGYmdmWwQnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHkDRZ0l+6cPuzJB28geWPSjozfxFtOpS4VtIySc91dTwdJalI0nmSTm2j3EGSXpO0QtKQdqw/r8eypJmS3pN0/ka8V5Kul/SBpCtzEV9HSHpH0uGtLDtJ0nJJT7Xn71MwCSf9cBZJ2ipj3pmSHs2YDkmrJK2U9K6k30gqzlj+qKSqdPkSSf8nadsstj1Z0uSNiDkznpWSPkznH5wuu71Z+d3T+Y+2tL7NVUSMjYhHoeMnjPSza8j4TFdK+mfGukPSuc3ec246f3LGOuZv/B51qv2BI4BhETEB1h7Xs9N9uy/zHz49if1c0tL05+eSlO3GJPVK/y/eSY/N/0j6u6R9mpWbLGnkRuzPFcBZwGVtXER8Hbg3InpFxAJJ+0l6slkMnfp32phjLyJ2B74IXLgRm9wVOAkYHRHnZMTR5AIr3c9lkk5Op4+TNENSZXqeeljS9i3szx8k3dDC/N0lVUvqvxExAxARtwJbp5OfzvZ9BZNwUsXAuW2U2T0iegIHkfyxv9Bs+Tnp8lFAX+DSzg6ypXjSn74Z8xcD+0oakDFvEvBGjuMpBAsyPtOeEXFsxrI3gM81K9+pn6ukks5aFzACeCciVqXrPhj4KXAc0B94G7g5o/xZwPHA7sA44Fjg7Gw2JKkceBjYDTgG6A2MBm4Bjk7LfE/SAelbSiT9QNJHs1z/D4H/Ag4EDgYmSzquleL9gVkZ088AO3fkJJlDrwB9Mi9es9QfWBIR77dWQNKRwB3A5yPiFkk7ATcA3wT6ANsDvwPqW3j79cAnMy/CU6cBd0XEB+2Mt4mIqCX5vxnQVtlGhZZwfgl8S1LftgpGxGzgSWB8K8s/AP4B7Cppb0nvN6sNfVLSzObvk7S1pLskfZhWlR+XtDGfcw3JgdZ4VVNMkiBvau0NkrpJ+kt6ZfuhpOclDUqXDZE0JY1ptqQvtrKOeyWd02zeTEmfTF/vImlqup7XJbV4dSPpEEkvZ0xPlfR8xvTjko5PX78j6XBJE4HvASelV++Zn+8ISU8qaWJ5QNLWbJzngR6SxqbbHgt0S+eT/nPeCwzJqCENkVQu6TJJC9Kfy9IT9NorbUnfkbQQuFbS6ZKeaPaZRHrCQNLHJP073Z93JX2rhc/wDODPJBceKyVdRJII/hYRsyKiBvgRcKCkHdO3TQJ+HRHzI+Jd4NfA6Vl+NqcBw4DjI+KViKiPiFUR8feImJyWuRyYSHJc/gGYFRHPZHHcnE2StA6NiCURMQs4BPilpP1biKUEaGiciIh64CHgqHR9Lf6d0uJlkm5IP9tZkioyYhoi6R+SFkt6W9LX0/ktHnuSPi/p1XRdc9L9aK4xzvZeaDTZx+YkHQPcBpwSEXeks8cDb0fEQ5FYERH/iIj/NH9/RDwNvAt8KmOdxcApwA2SdlRSO1qqpKZ0UzbnzmYaaMd+F1rCmQY8Cqz3z9ucpF2AA4DZrSzfmuQP9WJEPA8sBY7MKHIayZUGETE54x/ym8B8YCAwiOQg3tj7B93Auqvxo0iupBZsoPwkkque4SRXHV8C1qTLbknjGgKcAPxU0qEtrONm4DONE5LGkFxl353+k08F/gpsQ3LS+X1aprnGK9KtJZWSXG0PUdJk0x2oAB7PfENE3Edy9X5rWjPZPWPxKcDn0+2WkcXfeANuZN3nOimdboxhFcmJMbOWtAD4PvBRkn/43YEJwA8y1jmY5Ip1BEktoy1XA2dHRC+SppWHmxeIiKtJ/oZPp3E0NttkNpE1vt41/T0WyEzUM9N5SWHpJUmntBLT4cD9jbWpDYiM341X1q0eN+m+/DEi9o2IDzP2782IGBURzRPzMGAM0PwkejfwsfS9rf2dAD5Bcrz3BaYAV6brLQL+SfKZDAUOA86TdNQGjr1FrKvtfR64VNKezeJaDFSRNH1mJY3lsBb2sdGxJMflCRFxT8b8F4BdJF2aXtT1bGNTmecQSP7GpcA9JMfOz0jOCaNJzhuTs92H1Dzgv7R+LapFhZZwAH4IfE3SwFaWvyBpFfAqSXL6fbPlVyjpS5kJvAd8I51/PfBZACXV+qNITrzN1QLbAiMiojYiHo8N37DuhbQ28qGkKzIXRMRTQH9JHyE5aNZrj21h2wOAndKr0+kRUSlpOLAf8J2IqIqIGSRXzs2blgBuB8ZLGpFOnwr8X0RUk/zjvRMR10ZEXUS8SFILPLH5SiJiDUmt4UBgL5LP88k0jo8Cb0bE0jb2J9O1EfFGut7baKVmmhqS8Zl+2EIt7C/AZ9JEeHI63ZZTgYsjYlFELAYuIrnoaNQAXBgR1WmMbakFxkjqHRHLIuKFLN4DcB/waUnj0sT9Q5ITf490eU9geUb55UBPKenHiYhxEdHScQtJm/zCxglJ49PPr1LS6+nsc4EHSE7oXwZ2V9KktqHjJmuSfkVyEns0Ih5sYd+PVNstBk9ExD1prehGkgsEgL2BgRFxcUTURMQc4E+krQgtiYi7I+KttDbxGMm+H9CszBrg28AUSTOy2Mf+JBeCXwP+XyvFDgHeJPmfydzWHJLmyKEk/wdLJF23gcRzI3BQmsQh+Z//a3pumh0RU9NjdjHwG5Kuhva4jOR8V9lYm92Qgks4EfEKcBdwQStF9iT5pzwJ2Adonpm/HhF9I2JoRJya/iEgOSkdm2byTwOPR8R7Laz/lyS1pgfSKnhrcayNJ91e34j4egvLbwTOITkAb29hefOy9wO3pM0+v0hPqkOADyJiRUbZuSQHbRNpmbtZ90/4GdY1440A9sk8mZOcWAa3Es9jJP8cB6avHyU5oA9Kp9tjYcbr1SR/w9YsyPhM+0bEbZkL0+aH2SRXtG9GxLwstj+E5DNrNDed12hxRFRlsZ5GnyK5Wp8r6TFJ+2bzpvQkfCFJon8n/VlBUnsFWElyNd6oN7CyjYueRktJTh6N25oRSb/iJ4HydN5PI+JfaZG6iPhRRDzTxnGTtYj4Fkn/6ZGZTWHpsiUkfVYT2lhN82Olm5J+tRE0uxghaYEY1NqKJB0t6RklTcgfkvzNtm5WpoR1FyB7ZLGPH5Ccd64nqTm35H+AauAOpU23Ge9/JiI+HREDSZLfga2tJz3W/wV8Nk1Kx5NeuEoaJOmWtEm3kuQc196m6s8DlUD/iPi/tgoXXMJJXUgycmS9EypAerVyG/A0yRVimyJpD3+a5J/vNDKaYZqVWxER34yIHUiq9t+QdFj7d2GtG4GvAPdExOo2YqyNiIsiYgxJ5+wxJFc0C0hqSr0yim9H0r7bkptJagD7kvRvPJLOnwc81uxk3jMivtzKeponnMdoO+Hk6/bljR2vLdUaW4phAckJq9F2NG3ebP6eVayrdSCpSVKOiOcj4jiSJsI7SK5WsxIRv4uInSNiEEniKSFpboWkoz2zKXJ3mna+b8hDJCf6NptHImlGfqfZ7NaOm3aJiDdJ9qelptq1zWq0/1iZR9L/kXn89oqIFteXnuj/AfwKGJQm38amqEyDgH7AHVkmdiKijqR5r6V9hOT4+RhJE/nf0gvHltbzPPB/rGtSbcn1JOesT5Hs//R0/k9J9nm3iOhN0oKT9YjG1GjgkYhY3mZJCjThRDIg4FaSoZUbcgnwxeYngw24gaTqvBvJH3k9ko6RtFPahLGcpI271Y7BtkTE2yQn6NauhDK3fYik3dKOwUqSZpuG9Ar+KeBnSgYWjAPOoPWmpHtITq4Xk7RpN8Z/FzBK0mmSStOfvSWNbmU9TwEfIbkifS6SjuIRJDXLf7XynveBkVk0m3TUrSR9ci2d6N8HBkjqkzHvZuAHkgam/Xs/ZMNNcTOBsWmzVDcy2sYllUk6VVKfSEb6VJLlMZL+/XZVYjvgKuDyiFiWFrmB5CJnqJJO9G8C12Wz7vS97wG3p9soTmOvaON9jVo7bjZGNUlfXUvb+Hj6uqW/04Y8B6xQMrije7p/u0raO2N9mcdeGUnNbjFQJ+lomvbjNmpMBu1qPqT1fQTWtjZMJLlw/msa7/6SvihpG1jbF/0Jkj7TtV+paLaqf5BcIF1Eknwa9SKpES+XNJTWm/c2pJR27HdBJpzUxazfXNZERLxMcuLL9oO+neQf6vYN1DZ2Bh4k+UM+Dfw+IjbqSi8jzidiXYfohgwG/k5yAnuVpBbRWBP7DDCS5Kr8dpL+huZt5I3bqyZJqIeT0U+V/gMcSdJssoCk6eLnpM0tLaxnFUknZ+OIKkg+k7kRsaiVffhb+nuppGz7NdotItZExIPRQn9LRLxGkmDmpE0vQ4AfkwxKeQl4mWS/fryB9b9Bcgw+SNIW/0SzIqcB76RNGV8iaZpE0nZKRklt18qqu5H8TVaSnECfJml+afRHkivnl0lqCXen80jXP0utfOkybRI8BPh3+r5K4HWSvo82v2vR2nGzkRpo+fz0ArCtpMGt/J02FF89Sa1/PEnT3BKSvszGhNXk2EuP96+TXJQsIxm4MqWFVTeOXm1vgm1tHzNj/pBkMMIokguCSpIE87KklST9WrcDv0jfMpzkQi9zHatIks4wmjZzXkTSxbCc5O/dapNYeoHUUk25mHbst7KsAVpK0lsko4taPFmbWcdJakyqX2peU5J0LfCviLi2S4JrRsmglCsjYpt2vm8USXPnqLQlozNi+TPJsPn7O2N9bWyrO8mF1NUR0XzwVYsKuYbT6SR9iqTNc70hrGbWqX5DUhNZ3ELN5UqSZq4uJ2k6SV/Id9r73rQW/DvgcTUbobqxIuLMPCWbT5MMnHmfdvQ/uoaTJSW3kxkDnJaPP6iZbb4k3Uuz4dOpn0bET/Mdz6bCCcfMzPLCTWpmZpYXnXmTwc3C1ltvHSNHjuzqMMzMNhvTp09fkn7RtEO2uIQzcuRIpk2b1tVhmJltNiTNbbtU29ykZmZmeeGEY2ZmeeGEY2ZmeeGEY2ZmeeGEY2ZmeeGEY2ZmebHFDYs2KwQNDQ28N38uSxa/T12UgNr7GBMzIIIS1bH1wEFsO2wERUW5rYM44Zhtht569SU0/3l2mflryta8jzr06BnbUoWKqOk+iLfHfo15s0cyZvw+9O7f3od+Zs9Namaboco1tezw7PcoX/2ek41tNEUD5avfY9T0yZT0Gsjf//QLKpctzdn2nHDMNkcqoqihpu1yZlkoaqhBRcVULlvCrGmP5247OVuzmZltVrr16MmyxQtztn4nHDMzS4mG+vqcrd0Jx8xsE3X6HWs4/IZVXR1Gp/EoNTPLq9PvWMP1M2v5f/9Vxi+O6LZ2/vzKBoZfupJHJvXg4JEl6KLKtcu6l8DIvkWcuWcp39i3fO38dz5sYPvLV66d7l0Ou2xdxPf2L+e4XUqzjilzW43Ki6HqB70BGHnZCuYuDy6fWM7X9ylvUu78+6q47NkaDtu+mAc/txUADRH85ukarptRy9sfNlBaBCP6FnHsqBJ+fGi39ba1pXDCMbO861YCVzxbw1f3LmNE39YbWq48uhufGlPCmlp44K06zrm3ip5l4qy9ypqUu/Pk7kwYWsyyNcHPn6zhU7et4YkviI8Oy/4U17itRs2/2bRdH/HnF2qbJJyquuCGl2oZ0adp6Ysfq+byZ2v47dHd2HdYCVV1wSuLGnhmfu6aqzYHblIzs7z7r+HF7D64iO89XLXBcn26weCeRWzfr4izK8oYN6iI+9+qW69c/+5icM8iRg8s5k/HdqOsGO58rY5H36mj+OJK5i1vOnT8hpk19LmkklU1sd62Gn8G9Wx6ejx5bClzljXw7Px12//7v2vp1w0OGtk0sd3xWh1n7FHGZ8eVsWP/IsZuU8xJu5Zy6cSmtZvrZ9Qw5ncrKftRJcN+s4IfPFxFXUPQkqlvJfsyv7Lpvtz6Si09flJJZXXyvvdXNnD6HWsY+MsV9PpZJftds4p/zV3/M+sKruGYFZCDr1u/vf/TY0v5yt5lrK4NPnbT6vWWnz6+lNPHl7FkdQMn3LZmveVfrijjpF1Lmbe8gdNub7r80dO32qg4BfzqiG4cdN1qzv9oPRVDijdYPiJ4+O16Xl3cwKgBG75OLimC0mKobYCDR5awc/8irnmxlgsPXlcz+dMLtZyyaylblWV/h4Ze5eLkXUv50wu17JPWnK6aXsuZe5bx2pKmSWDbXuKxuXW8W9nA0N4tx3v3G7V8YUoVPz6knE+NKeHF9xr40t1rEPCjFprdDtuhmG17iptequU7+6/bl+tn1nL8LiX0LhdraoNDrl/N6IFF3HtqD/p2E7e+UssRN65mxtlbMXrghj/nXHMNx8y6xAEjSjhulxK+9UDrtZwzp1TR86eVlP14BYffuBoJztunrNXyVXXBRY9VU1kNh++QJIWz9irlmhk1NERSA3htST1P/Kd+vWa5xm01/vzoser11n/WXmXc8kotK6qD15bU88z8ej4/fv2+okuP6saaWhh+6Uo+cuVKJt2xhpteqm1Se7nkyRo+NbqE7x5QzqgBSQ1o8kHl/OrpGmrq16/lFEl8dlwpN75Uu3be+ysbeOCtOibtnuzLrbNqqawObj2hOxVDitmpfxHfP7Cc/YYX88fpteutM99cwzErIBuqcfQo1QaXb92jaIPLh/fZ8PKN8fPDyxn7+1VMeb2WPbdd/+r7J4cmnf8LVzZwwYPVnDimZG3tItORN66mSLCmDvp1E5ceVc7EnZJyk3Yv5fsPV3P/7DqO3rmUP79Qy17bFrFHs+01bqtR/+7r134mDE1O4je/UstrSxo49iMl6zW9AeyydTEvf3krZixs4In/1PPU/DrO/OcaLn2miMc/vxXdS8WsRfWcNLZpTeagkSVU1VXz1gcNLdZGJu1eys+frOGF9+rZc9tibnq5lm22EofvkJR9/t16Fq4M+l6yosn7quuhe2nX32/PCcfMusyoAcWcvVcp33mwmntP7bHe8kE9xU79i9ipfxF3nFzEqN+uZI9tizlwRNNT17XHdWevIcX07ZYkzkwDehRxwpikKeywHUq4YWYtPz606UizzG215ay9yvjfaTXMWx7c9MnurZaTxB7bFrPHtsV8bZ8ynvhPHQdcu5rbZtUyaXzrtbQNGT2wmIohRdwwM0nQN8ys5bPjSikuSpJJQ8DogUXcftL6n2WP7Aft5Yyb1MysS114UDkLVjRw1fQN36pn6x5FfHXvMr5+bxURTZuchvZOkkXzZNPo7L1K+ecbdfxxWi1r6oLP7LrxZ9/PjivlzaUN9CqHI3bMvk9k9NZJbItWJbGP3aZ4vc78x96po3sJ7LiBxDdp9zJufqWWF96rZ+b7DXxu93X7UjGkmDnLGuhdztpE3fgzpFfXn+67PgIz26IN3KqIC/Yr57Jn2r433DkTkg76W15p36ir/bcr4SMDivjW1CpOHltKr/KNb17qXS7e/UYvXvpST4paeSzEp25bza+fqubpeXXM/bCBp+bVcdrtVZQWwcdHJbWz7+5fxj9ereOSJ6p5Y2k9t82qZfJj1Xxz3zLKiluP7zO7lrBsTXDGlDXsuW0Ru26zLumdOq6U7fsW8fG/ruaBt+p458NkVN3PHq/mjte6vg/HCcfMutz5+5axdY+2k8CgnkV8bvdSfvhodavDh1vzxT1LqalnvcECG6NPN20waU3csYT73qrjk7etYdSVKznxb2soK4bHTu/BmLRv5mM7l3LNJ7px/cxadv39Ks6/v4qvVJQ1GU3XkgE9ivj4qBJmLGzgc+Oa1tS6lYjHTu9BxbbFfP7ONYz67Uo+edsanltQz4g+XX+6V/OqaaGrqKiIadOmdXUYZh0yffp09vrnoV0dxmbl21OrmDqnjhfP7tnVoWySph/7MM/ffQPDth/FMad9tckySdMjoqKj2/CgATMraMurgjeWJn1EVxy95d5WZlPQ5XUsSRMlvS5ptqQLWlheLunWdPmzkkY2W76dpJWSvpW3oM1ss3HcLas58LpV/PfoUj47bhMYqrUF69IajqRi4HfAEcB84HlJUyLi3xnFzgCWRcROkk4Gfg6clLH8N8C9+YrZzDYvnf3dIdt4XV3DmQDMjog5EVED3AIc16zMccD16eu/A4dJydAQSccDbwOz8hOumZltrK5OOEOBeRnT89N5LZaJiDpgOTBAUk/gO8BFbW1E0lmSpkmatnjx4k4J3MzM2qerE05HTAYujYiVbRWMiKsioiIiKgYOHJj7yMzMbD1dPUrtXWB4xvSwdF5LZeZLKgH6AEuBfYATJP0C6As0SKqKiCtzHrWZmbVbVyec54GdJW1PklhOBk5pVmYKMAl4GjgBeDiSLw8d0FhA0mRgpZONmdmmq0sTTkTUSToHuB8oBq6JiFmSLgamRcQU4GrgRkmzgQ9IkpKZmW1murqGQ0TcA9zTbN4PM15XASe2sY7JOQnOzKyTfPfBKq6ZUcuiVcG1x3Xj9I28Y/TmrMsTjpltWU6/Yw3Xz1z/RpI3f6o7ry2p56LHathjcBEvNLsFzcyF9Yz/Y/JE03nn92RY76K1t/x/+9yejOy76Y6BenZ+HZc8WcMdJ3Vnn2HF9CkXdQ3Bb56u4eoXa5n7YQPD+xRx3j5lfHVC00T0xtJ6vnZvFY/PradHqThhTAm/PrJbu55WuqlwwjGzvDtgu2JuO7Hps2T6dhOXPFHPwB7itSUNax8y1uiP02sY0UfMXZ6b+z/W1McG79LcEW9+0ECRaPKAt+8/VMVVL9Ry1THd2H1wMU/Pq+esu5KbfH4xvcHoyprgsBtWM25QMU+dsRUfrAm+cOcaPqxawy0nrP/Mm03dpntJYGYFq6wYBvcsavLTrSQ52fcuJ3lgWsbzcVbXBje9XMsZe6y7+n/nwwYOuHY1ANtfvhJdVMnB1yU1oIjgV09Vs8PlKyj7USU7XrGCy55p+sjokZet4AcPV/GVu9cw4BcrOODa5L26qJK/vNT0UQmH37CK0+9Ys3b6ztdq2eOPK+nxk0r6XlLJhD+t5MX36lvc19PvWMNpt1fREMm6dVElANfPrOWb+5bx36NL2aFfEaeOK+XMPcr4yePr4vzry7UsWR389ZPdGT+4mEO3L+F3H+vGrbPqeHtZQ/s+9E2AazhmBeK8+6qYsbDlk16ujB9czGUTO/+GmGftVcrH/7qaXx/VjR6l4pZXahnSq4gDRqyr8QzvLe48uTvH3bKG587ciuF9tLaG8vvna/mfR6q5fGI3DhlZzENv13PefVX0KhNn7LkuaV3xbA3f2LeMp8/oQV2W5++FKxs48W9r+PGh5Zw4ppSquuDFhQ2UtHL5fvnEbuwxuIhvPlDN/G+sayasqoNuzc7A3Uth7vJg7ocNjOhbxJPz6tl3WDF9uq2reR25YwlFgifn1bF9v82rH8gJx8zy7tF36un508q100N7F/H6OetOxvtvV8Kw3kXc+kotn9+jjKum1/LFPZveeLO4SPTvnpyIB24lBvdcd8a/5MlqvjahbO2zb3YeUMzrSxr4yePVTRLO3kOLmXxw+xLmeyuC2gb49NjStf1Gowe2/uTPPt20NmFkxnj0ziVc8WwNh21fwq7bFPHcu/Vc82LSt7VgRZJw3lvRwOCeTZv5SouT/X5vxeb3aBknHLMCkYuaRq7sM6yY649f14fTUu3gi3uW8qcXatlrSDEzFtZz1yndeWVR29WQyupgfmVw4IimSeCgkcVc/mwNq2uDHqXJSXzCkOwfEd1o3KAijtqxmF1/v5Ijdizh4BHFfHJ0KcPb+YCzyyd240t3rWH8H1chYEgvccYepVzyZA1Fm994gKw44ZhZ3nUvgZ36b/gE/bndy/juQ9V84/4q/nt0CVv3KAI6t9+ipZFeApo/l7I2Y7PFReLeU3vw/IIGHpxTxz9ereOCh6r524ndOWZU9o8/6N9d3HZiD2rqg0WrgiG9xB+mJTWcHfoln822vYqYt7zpPtfWBx+sCbbttfllJQ8aMLNNUv/u4oQxpTz0dj1f3LPlvorGPpv6jHNy73IxrLf419ym/VmPvVPP9v20tnbTmm22Egsymquq64J/L2560pfEhKHFfO+Acv71+a04aEQx185Yf6h3NsqKxbDeRRRJ3PxKLQeOKGbgVsmpeb/hxTw9v57K6nXxTJ1TR0PAfsM3v/qCE46ZbbL+dGw3Fv+/nhy6fcsn1xF9RJHgnjfrWLSqgeVVyYn5u/uX89vnavjT9BreXFrPH6fV8L/Tavje/uVtbvPwHUr4w/Qanp5XxyuL6jn9zjXU1K874T81r44fPVbNs/Pr+M/yBh6aU8dL7zcwZut1p9NdrlzJlc/VtLT6tZ5/t56/zarlrQ8aeHpeHSfctpoZC+u5IqNp9JTdStm6hzjlH2uYubCeR96u46v3VHHS2BK277f5nb43vxRpZluMbiVaO1y6JYN6FvGzw8q55Mlqzru/igO2K+bR07fiyxWlrKoJfvpENV+5JxjeW1xyeHmTAQOt+dWR5Xzxn8FRf1lNn27ie/uXs3jVuoTTp1w8Pb+e3z1fw7KqYHBPcepupfzPQeuS2etLG1iyesPNf9X1wUWPVfPWsgbKiuHAESU89YWt2G3Qun6lnmXiwdN68LV7q9j36lV0LxUnjC7hN0dtPv11mRTNGysLXEVFRUybNq2rwzDrkOnTp7PXPw/t6jCsgEw/9mGev/sGhm0/imNO+2qTZZKmR0RFR7ex+dXJzMxss+SEY2ZmeeGEY2ZmeeGEY2ZmeeGEY7Y5iqBBHmRqnaNBJRC5vxnoRiUcSad3chxm1g5lqmN13490dRhWIFb3/QgNVZVtF+ygdiccSTsCv5U0KAfxmFkWho7cmTcn/JSV/ca6pmMbrUElrOw3ljf2/jEL/zOHiAaKS7O/PU97tXmkSjofOB/oQ5KgPgAuBGZJKgNWAb+IiEtzFqWZNdF/wADmlZQya/xkSrbqB3LruG2EaKChqpKF77xF5ZIFrFm5gq0HD8vZ5rK5NDob2Dsi3m82/zcAkrYBHgCccMzyaMzuFdx90//y9r9nslWfvhQVt//Ox2YADQ31rKpczrAdPsK4fQ7O2XaySTgXtZBs1oqIRZJ+3okxmVkWSsvK+fipX2bao/eycN4camur236TWQtKSkoZs+d+7H3wxyjvnrtHV7eZcCLi5s4oY2adr7SsnH2PPL6rwzDLSocafiUN7qxAzMyssHW0p/HqTonCzMwKXocSTkR8vLMCMTOzwpbNsOj+G1oeER90XjhmZlaoshmlNh0Ikkd9NxfADp0akZmZFaRsRqltn49AzMyssLXrnhiS+gE7A2ufbxoR/+rsoMzMrPBknXAknQmcCwwDZgAfBZ4G/JxbMzNrU3tGqZ0L7A3MjYhDgD2ADzsagKSJkl6XNFvSBS0sL5d0a7r8WUkj0/lHSJou6eX0txOfmdkmrD0JpyoiqiBJAhHxGtCh+6NLKgZ+BxwNjAE+I2lMs2JnAMsiYieS+7U13kZnCXBsROwGTAJu7EgsZmaWW+1JOPMl9QXuAKZKuhOY28HtTwBmR8SciKgBbgGOa1bmOOD69PXfgcMkKSJejIgF6fxZQHdJ5R2Mx8zMciSb7+H8P+DmiPjvdNZkSY+QPK7gvg5ufygwL2N6PrBPa2Uiok7ScmAASQ2n0aeAFyLCdy80M9tEZTNoYAjwtKR3gJuBv0XEYzmNqh0kjSVpZjtyA2XOAs4C2G677fIUmZmZZWqzSS0izge2A34A7Aa8JOk+SZMk9erg9t8FhmdMD0vntVhGUglJzWppOj0MuB34XES8tYF9uCoiKiKiYuDAgR0M2czMNkZWfTiReCwivkySFC4FzgNafU5Olp4Hdpa0ffr00JOBKc3KTCEZFABwAvBwRETan3Q3cEFEPNnBOMzMLMfadfNOSbsBF5OMLKsGvtuRjUdEHXAOcD/wKnBbRMySdLGkT6TFrgYGSJoNfANoHDp9DrAT8ENJM9KfbToSj5mZ5Y4iYsMFpJ1Jah4nA/UkI8luiYg5uQ+v81VUVMS0adO6Ogwzs82GpOkRUdHR9WQzaOA+ksECJ0XEKx3doJmZbZmyuXnnjvkIxMzMCttGP4BN0oOS7pV0TGcGZGZmhaldd4tu5nPAtiQ38TQzM9ugdiccSaXArsC7ETGd5AFtZmZmG9Rmk5qkP6Tf5kdSH2AmcAPwoqTP5Dg+MzMrENn04RwQEbPS158H3kjv0LwX8O2cRWZmZgUlm4RTk/H6CJK7RRMRC3MRkJmZFaZsEs6Hko6RtAewH+kdotP7mnXPZXBmZlY4shk0cDZwBTAYOC+jZnMYyb3MzMzM2pRNwlkTERObz4yI+9MRa2ZmZm3KpkltqqSRzWdK+gJweadHZGZmBSmbhPMN4IH0Jp4ASPoucD5wUK4CMzOzwpLNvdTukVQN3CvpeOBMYAJwYEQsy3F8ZmZWILJ9ANtDJN/BeRTYATjUycbMzNqjzRqOpBVAAALKSUanLZIkkoeB9s5tiGZmVgiyaVLrlY9AzMyssG304wnMzMzaI5ubd17dGWXMzGzLls0XPydI+hKwCGgA5kfENEkHA/2BbUgeV2BmZtaqbBLOJOCzJEOhi4GxklaTDCKYDawGvp6zCM3MrCBkM2jgBeCFxun0pp3/AYZGROQwNjMzKyDtHjQQEXXABCcbMzNrj40apRYR85vPk/TbjodjZmaFqjOHRe/XiesyM7MC4+/hmJlZXjjhmJlZXnRmwlEnrsvMzApMZyYcP4zNzMxatVEJR9IRzedFxHUdjsbMzApWuxOOpL7ArZK6d344ZmZWqLK5eedJkv4laaakWcAbwF+AdyS9JOlpSZ/Z2AAkTZT0uqTZki5oYXm5pFvT5c9KGpmx7Lvp/NclHbWxMZiZWe5lcy+1i0jup7aQ5OadSyJijaQfAn2AQcAfgJvbu3FJxcDvgCOA+cDzkqZExL8zip0BLIuInSSdDPwcOEnSGOBkYCwwBHhQ0qiIqG9vHGZmlnvZJJw/RcSzzWdGxIfAh8BcSTdt5PYnALMjYg6ApFuA44DMhHMcMDl9/XfgyvRpo8cBt0RENfC2pNnp+p7eyFg26LzzzuO6665bb/4222zDkCFDaGho4KWXXlpv+eDBgxk8eDC1tbXMmjVrveVDhgxhm222obq6mldffXW95cOHD2fAgAGsXr2aN954Y73lI0aMoF+/fqxcuZLZs2evt3yHHXagd+/eVFZWMmfOnPWW77TTTvTs2ZNly5Yxd+7c9ZaPGjWKHj16sHTpUubNm7fe8tGjR1NeXs6iRYtYsGDBesvHjh1LaWkpCxcuZOHChestHzduHEVFRSxYsIBFixatt3z8+PEAzJs3j6VLlzZZVlRUxLhx4wCYO3cuy5Y1fep5aWkpY8eOBWDOnDlUVlY2WV5eXs7o0aMBmD17NitXrmyyvEePHowaNQqAN954g9WrVzdZ3rNnT3baaScAXn31Vaqrq5ss7927NzvssAMAs2bNora2tsnyfv36MWLECABeeuklGhoamiwfMGAAw4cPB2DGjBk052PPxx5kd+yNHTuW8ePHc9lll623n/nUZpNaRPy6M8q0YiiQeSTNT+e1WCa9j9tyYECW7wVA0lmSpkmatnjx4o0M1czMOiKbGs5akvoBOwPdGudFxL86O6jOFhFXAVcBVFRUbNRNRy+77LIuvzowM9ucZZ1wJJ0JnAsMA2YAHyVpvjq0A9t/FxieMT0snddSmfnpoxH6AEuzfK+ZmW0i2jMs+lxgb2BuRBwC7EHSh9MRzwM7S9peUhnJIIApzcpMIRm0AHAC8HD6aIQpwMnpKLbtSWpez3UwHjMzy5H2NKlVRUSVJCSVR8Rrkj7SkY1HRJ2kc4D7SZ4mek1EzJJ0MTAtIqYAVwM3poMCPiBJSqTlbiMZYFAHfNUj1MzMNl3tSTjz0y993gFMlbQMWH9oSTtFxD3APc3m/TDjdRVwYivv/Qnwk47GYGZmuZd1womI/05fTpb0CElfyr05icrMzApO1n04km5sfB0Rj6XNXdfkJCozMys47Rk0MDZzIr1LwF6dG46ZmRWqbO6l9l1JK4BxkiolrUinFwF35jxCMzMrCNncaeBnEdEL+GVE9I6IXunPgIj4bh5iNDOzAtCeJrXvS/qspP8BkDRc0oQcxWVmZgWmPQnnd8C+wCnp9Mp0npmZWZva8z2cfSJiT0kvAkTEsvTuAGZmZm1qTw2nNh2ZFgCSBpI8H8fMzKxN7Uk4VwC3A9tI+gnwBPDTnERlZmYFpz13GrhJ0nTgMEDA8RGx/lObzMzMWtCu5+EAbwKVje+TtF1E/KfTozIzs4LTnufhfA24EHgfqCep5QQwLjehmZlZIWlPDedc4CMRsbTNkmZmZs20Z9DAPGB5rgIxM7PC1mYNR9I30pdzgEcl3Q1UNy6PiN/kKDYzMysg2TSp9Up//yf9KUt/IP1OjpmZWVvaTDgRcRGApBMj4m+ZyyS1+CROMzOz5trTh9PSnaF9t2gzM8tKNn04RwMfA4ZKuiJjUW+gLleBmZlZYcmmD2cBMA34BDA9Y/4K4PxcBGVmZoUnmz6cmcBMSX+NiNrG+ZJ6AYdLmhgRZ+cySDMz2/y1515qtZJ2ByYCR6XvfQi4JkexmZlZAcmmD+ckkiSzBzATuA840XccMDOz9simhrMTcGVETG+zpJmZWSuyGRZ9XVvJRtLQTorHzMwKVDYJ554sytzV0UDMzKywZdOkNk/SM8AikkdKzwf+TDIkuh8wEJidswjNzKwgZJNwjgV2A/oCxcB44FHg28DrwGpgRi6CMzOzwpHN93ACeClj1iOSKiLiqtyFZWZmhaY991JbKyJO7eiGJfWXNFXSm+nvfq2Um5SWeVPSpHReD0l3S3pN0ixJl3Q0HjMzy62NSjid5ALgoYjYmeQLpBc0LyCpP8ljrfcBJgAXZiSmX0XELiTfD9ovveebmZltorJOOJImd/K2jwOuT19fDxzfQpmjgKkR8UFELAOmAhMjYnVEPAIQETXAC8CwTo7PzMw6UZsJR1KRpKuB8k7e9qCIeC99vRAY1EKZoSSPtm40P52XGV9fkoEND3VyfGZm1omyGaV2FzArItr97BtJDwKDW1j0/cyJiAhJ7X56qKQS4GbgioiYs4FyZwFnAWy33Xbt3YyZmXWCbBLOXsCPNmblEXF4a8skvS9p24h4T9K2JN/zae5d4OCM6WEkQ7IbXQW8GRGXtRHHVWlZKioq/FhsM7MukE0fziHAVZL26eRtTwEmpa8nAXe2UOZ+4EhJ/dLBAkem85D0Y6APcF4nx2VmZjnQZsKJiH+TdN7/spO3fQlwhKQ3gcPTaSRVSPpzuu0PSGpXz6c/F0fEB5KGkTTLjQFekDRD0pmdHJ+ZmXUiJd/rzKKg1CsiVuQ4npyrqKiIadOmdXUYZmabDUnTI6Kio+tpz/dwdu7oxszMbMvVnoTza0mvSvqRpF1zFpGZmRWkrBNORBxCMoBgMfBHSS9L+kHOIjMzs4LSrlvbRMTCiLgC+BLJHaJ/mIugzMys8LTn1jajJU2W9DLwW+ApfDsZMzPLUjZf/Gx0DXALcFRELMhRPGZmVqCySjiSxkXEvpJ2c7IxM7ONkW2T2hck7QyckctgzMyscGVzt+gL03LPAEWSPFDAzMzaLZtHTF8k6RNp2akR8c/ch2VmZoUm2ya1CRHxFZKnbpqZmbVbNk1qx5E8II2I+B9Jz0qak/6ckPMIzcysIGRTw/k2TR8dUA7sTfKcmi/nICYzMytA2QyLLouIzMc8PxERS4GlkrbKUVxmZlZgsqnh9MuciIhzMiYHdm44ZmZWqLJJOM9K+mLzmZLOBp7r/JDMzKwQZdOkdj5wh6RTgBfSeXuR9OUcn6O4zMyswGTzPZxFwH9JOhQYm86+OyIezmlkZmZWULK+eWeaYJxkzMxso7TreThmZmYbywnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzywgnHzMzyossSjqT+kqZKejP93a+VcpPSMm9KmtTC8imSXsl9xGZm1hFdWcO5AHgoInYGHkqnm5DUH7gQ2AeYAFyYmZgkfRJYmZ9wzcysI7oy4RwHXJ++vp6WH+Z2FDA1Ij6IiGXAVGAigKSewDeAH+c+VDMz66iuTDiDIuK99PVCYFALZYYC8zKm56fzAH4E/BpY3daGJJ0laZqkaYsXL+5AyGZmtrGyfgDbxpD0IDC4hUXfz5yIiJAU7VjveGDHiDhf0si2ykfEVcBVABUVFVlvx8zMOk9OE05EHN7aMknvS9o2It6TtC2wqIVi7wIHZ0wPAx4F9gUqJL1Dsg/bSHo0Ig7GzMw2SV3ZpDYFaBx1Ngm4s4Uy9wNHSuqXDhY4Erg/Iv43IoZExEhgf+ANJxszs01bVyacS4AjJL0JHJ5OI6lC0p8BIuIDkr6a59Ofi9N5Zma2mVHEltWlUVFREdOmTevqMMzMNhuSpkdERUfX4zsNmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXjjhmJlZXigiujqGvJK0GJjbbPbWwJIuCGdT4H3f8myp+w3e943d9xERMbCjAWxxCaclkqZFREVXx9EVvO9b3r5vqfsN3veu3nc3qZmZWV444ZiZWV444SSu6uoAupD3fcuzpe43eN+7lPtwzMwsL1zDMTOzvHDCMTOzvNiiE46kiZJelzRb0gVdHU8uSbpG0iJJr2TM6y9pqqQ309/9ujLGXJE0XNIjkv4taZakc9P5Bb//krpJek7SzHTfL0rnby/p2fTYv1VSWVfHmguSiiW9KOmudHqL2G8ASe9IelnSDEnT0nldesxvsQlHUjHwO+BoYAzwGUljujaqnLoOmNhs3gXAQxGxM/BQOl2I6oBvRsQY4KPAV9O/9Zaw/9XAoRGxOzAemCjpo8DPgUsjYidgGXBG14WYU+cCr2ZMbyn73eiQiBif8f2bLj3mt9iEA0wAZkfEnIioAW4BjuvimHImIv4FfNBs9nHA9enr64Hj8xlTvkTEexHxQvp6BckJaChbwP5HYmU6WZr+BHAo8Pd0fkHuu6RhwMeBP6fTYgvY7zZ06TG/JSecocC8jOn56bwtyaCIeC99vRAY1JXB5IOkkcAewLNsIfufNivNABYBU4G3gA8joi4tUqjH/mXAt4GGdHoAW8Z+NwrgAUnTJZ2VzuvSY74knxuzTVdEhKSCHiMvqSfwD+C8iKhMLngThbz/EVEPjJfUF7gd2KVrI8o9SccAiyJiuqSDuzicrrJ/RLwraRtgqqTXMhd2xTG/Jddw3gWGZ0wPS+dtSd6XtC1A+ntRF8eTM5JKSZLNTRHxf+nsLWb/ASLiQ+ARYF+gr6TGC85CPPb3Az4h6R2S5vJDgcsp/P1eKyLeTX8vIrnQmEAXH/NbcsJ5Htg5HbVSBpwMTOnimPJtCjApfT0JuLMLY8mZtO3+auDViPhNxqKC339JA9OaDZK6A0eQ9GE9ApyQFiu4fY+I70bEsIgYSfK//XBEnEqB73cjSVtJ6tX4GjgSeIUuPua36DsNSPoYSTtvMXBNRPykayPKHUk3AweT3KL8feBC4A7gNmA7kkc2fDoimg8s2OxJ2h94HHiZde353yPpxyno/Zc0jqRzuJjkAvO2iLhY0g4kV/79gReBz0ZEdddFmjtpk9q3IuKYLWW/0/28PZ0sAf4aET+RNIAuPOa36IRjZmb5syU3qZmZWR454ZiZWV444ZiZWV444ZiZWV444ZiZWV444Zh1MUkXSDpV0mRJIWmnjGXnpfMqNrQOs82BE45ZF1GiCDgKeCCd/TLJFxUbnQjMyndsZrnghGOWR5JGps9guoHkm9/DgbKIWJwWuYP0ruWSdgSWA0u6IlazzuaEY5Z/OwO/j4ixwF4kzyVpVAnMk7QrSU3n1i6IzywnnHDM8m9uRDyTvp4I3Nts+S0kyeZ41t2exGyz54Rjln+rMl5PAJ5rtvwu4DTgPxFRmbeozHLMz8Mx6yKSxgKvpc+rWSsiVkv6DvBG10RmlhtOOGZd52jgvpYWRMQteY7FLOd8t2izLiJpKvC5jEf+mhU0JxwzM8sLDxowM7O8cMIxM7O8cMIxM7O8cMIxM7O8cMIxM7O8cMIxM7O8+P9TZo3Pe0+w8gAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEWCAYAAACaBstRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABMWUlEQVR4nO3dd3wUdfrA8c+z6RBC6J0QBJGqdDlBiiJY4awoKnpYTz099e7n2fXuLGfDesp5FjyxKyKKICKoIF2qSJFeQgglIaQnz++PmYRNSNmE7G7YPO/Xa1/ZnZmdeWZ3M898v9/5fkdUFWOMMcYXnmAHYIwx5vhhScMYY4zPLGkYY4zxmSUNY4wxPrOkYYwxxmeWNIwxxvjMkkYtJiIPi8j/grj9NSIypJz5c0TkusBFVHOI400ROSAii4Idz7ESEY+I3CEiYytYbrCI/Coih0SkZRW39ZaI/KOc+eki0t59PkBE0kRkuYh0r8K2uorIz+46BlQlXn8RkSEisqOc+V+IyF4Reaoy6z3uk4aIbBGRZBGp6zXtOhGZ4/VaReSw+2PZKSLPikiY1/w5IpLlzk8RkU9FpIUP235YRB6uQsze8aSLyEF3+hB33mcllj/ZnT6ntPUdr1S1q6rOgWNPYO5nV+D1maaLyBde61YRub3Ee253pz/stY4y/8kCbCAwHGitqv2g6He90d23r70Pqm6SeVJE9rmPJ0VEfN2YiNRz/y+2uL/NbSLysYj0L7HcwyLSrgr78wJwAzChghOBPwHTVbWequ4SkdNEZF4VtlcmVY1V1U3u85+AeGA9ML4Kq/sDsAmId9eFiLRzf1fh7msRkRfdZNhKRCJF5BkR2eF+l1tEZELJFYtItIgcFJFhpcx7TkQ+rkK8RVT1fGAocLeIxPv6vuM+abjCgNsrWOZkVY0FBgOX4XzZ3m5155+I8yN6rrqDLC0e9xHvNX0vMEBEGnlNG4fzozbl2+X1mca6/xSF1gNXl1i+Wj/XwoNENUkAtqjqYXfdQ4DHgFFAQ2Az8J7X8jcAo4GTgR7A+cCNvmxIRKKA2UB34DwgDugMvA+c7S5zr4gMct8SLiL3i8ipPq7/QeB3wOnAEOBhERlVxuINgTVerxcAHUWkoS/bqgpVLQB+ARpVtGwpGgJr3XUcRUQ8wGs4+z1YVXcCfwP6AP2Aeu68ZaXElQV8QInfrXvCeznwdhXiLbmN1e5Tn/c9VJLGU/iYLVV1IzAPOKWM+fuBT4BuItJXRPaUKJVcKCIrSr5PRBqLyDT3zGC/iPzg/mAqKweYAoxx1xuGk+TeLesN7hnJ/9wzzIMislhEmrnzWorIVDemjSJyfRnrmC4it5aYtkJELnSfnyQi37jrWScil5axnqEissrr9Tcistjr9Q8iMtp9vkVEzhSRkcC9wGXumZf355sgIvPEqa6YKSKNy/zkyrcYqCMiXd1tdwWi3emIU1KdDrT0Kqm0FJEoEZkgIrvcxwT3IFtUMhGR/xORJOBNEblGRH4s8ZmoiHRwn58jIr+4+7NTRO4u5TMcD7yOc/KQLiKP4BzMP1LVNaqaA/wdOF1ETnDfNg54RlV3uAemZ4BrfPxsrgJaA6NVdbWq5qvqYVX9WFUfdpd5HhiJ87t8FVijqgt8+N3ciJN4hqlqiqquwTm7fUpEBpYSSzhQdABW1XzgW2CE1/oLP/d7xakZ2CJHV3s1EJEv3c95odfnVOz78FLgbruyisVbQhjwJk6CGKKqe9zpfYHPVHWXOrao6qQy1vE2cJGI1PGaNgLn2D1dRK4VkbXufm5yP+/KUiqx76GSNJYAc4Cj/gFLEpGTgEHAxjLmNwYuAn5W1cXAPuAsr0WuAiYBqOrDXv9UdwE7gCZAM5yDYFXHaJnEkbOLEcBqYFc5y48D6gNtcM4YbgIy3Xnvu3G1BC4GHiutuItz1np54QsR6YJztvule0D9BpgMNMU5cLziLlNS4ZlhYxGJwDnrbSlO9UcMzj/QD95vUNWvcc6iP3BLCCd7zb4CuNbdbiQ+fMfleIcjn+s493VhDIdxDm7epZVdwH3AqTgnGSfjnB3e77XO5jhnmwk4Z/sV+S9wo6rWA7rhnOEXo6r/xfkOf3LjeMid5V3dVPi8m/u3K+CdbFe405yFRVaKyBVlxHQmMKOwVFMO9fqb7z4v83fj7strqjpAVQ967d8GVT1RVUsm19ZAF2Bbie1+CZxTYlpzoDHQCue7nCginbzmjwEeARrg/K//s4J92w70luIl/HK5pZ8+pcRb6F2gE07C3Oc1fQFwp4j8UUS6i5Rdjaiq84HdwIVek68CJqtqHpDMkdLhtcBzItLL131wbQfOLC8Ob6GSNAAeBG4TkSZlzF8mIoeBtTgJ5pUS818Qp21hBc6XdKc7/W3gSij6kYzAOXiWlAu0ABJUNVdVf9DyB/Za5pYKDorIC94z3B9KQ/ef4GrcJFWOXJxk0cE9S1yqqmki0gY4Dfg/Vc1S1eU4Z7Alq2kAPgNOEZEE9/VY4FNVzcb5UW5R1TdVNU9Vf8YpjV1SciWqmolz9n460Bvn85znxnEqsKHEP1BF3lTV9e56P6SMEqKrpddnerCU0tD/gMvdZDbGfV2RscCjqpqsqntxDkRXec0vAB5S1Ww3xorkAl1EJE5VD6jqUdUSZfgauFREerjJ90Gcg3fhGWgskOq1fCoQW3ggUNUeqlra7xacg29S4QsROcX9/NJEZJ07+XZgJs5JyM3AyeJUT5X3u/GZiDyNc/Cao6qzStn3s+TokvsD7uc+FyexeH/fn6nqIvfA+i7l/27A+S2kACkicocP8d6Gc0KZStnVRGfhlA4Plpj+OPAkzme1BNgpIuPK2VzRSaSIxOFUUb4NoKpfqupvbollLs53NKjMNZXuDpzq+AO+LBwyScOtm5sG3FPGIr1w/rEuA/oDdUvM/5OqxqtqK1Ud6x4gwPkxne+ebV8K/KCqu0tZ/1M4ZzQz3WJiWXEUxeNuL15V/1TK/HeAW3GK8p+VMr/ksjOA990qlH+5B8aWwH5VPeS17Facs7Ni3GW+xK0Wwzl7LKwSSwD6ex+QcX7wzcuIZy5OPe3p7vM5OG1Jg93XlZHk9TwD5zssyy6vzzReVT/0nqmq23C+o8dwktd2H7bfEuczK7TVnVZor1v37KuLcM6at4rIXPHxihv3QPoQTrLe4j4O4ZQiAdJxzjYLxQHpFZy4FNqHc8JTuK3l6rSzXQhEudMeU9Xv3UXyVPXvqrqggt+Nz1T1bpz2xLNEpE+JeSk4bTj9vCYfKFEyKvm9VOZ3A04bUBugpapO8CHeF3E+s+Y4B/HSnAc8JCLF2k/dE7uXVfU0nPbTfwJviEjnMtbzDjBUnAsfLgZ+c0/cEJGzRWSBONXGB3F+W5Wtwn0EeBinxFyhkEkaroeA6ynloAjgZuMPgZ9wztQqpE798E84/0BX4VWlUWK5Q6p6l6q2By7AKX6eUfldKPIO8EfgK1XNqCDGXFV9RFW74DQ4nodzZrILp8RSz2vxtsDOMlb1Hs6Z+ACc+v7v3OnbgbklDsixqnpzGespmTTmUnHSCNRwy5NwqhJLK72VFsMunKRZqC3FqwpLvucwR87+EZFiiVVVF6vqKJzqtik4pSefuAeajqraDCd5hONUXYLTeOxdrXcyxRuUy/MtzsG65IlUaTE8rKpbSkwu63dTKaq6AWd/Sqv2LFlF1aBEvCW/l8rqDCwo44SwVKqahHNsKC1egPk4yej5sqoGVTVTVV/GOcsvdT2quhWnSvdKnGPQ21B0AcMnwNNAMzfRf0XxakxfdAY+1zIa80sKqaShTiP3BziX7ZXnCeD6kv/Q5ZgE/BXn6pJPS1tARM4TkQ5udUAqTp2vT19CaVR1M85B9r6KlhWn8bm7OI3maThVIAXumfR84HFxGst74FxWWFa1zFc4B8hHcdoXCuOfBpwoIleJSIT76FvOmdF8nLrcfsAidRo/E3BKeN+X8Z49QLtSqiCq2wc41QalHaz3AI1EpL7XtPeA+0Wkidve9SDlV2utALq6VTzROGdwAIhzqeVYEamvqrk435VPvxH3++smjrbAROB5VS2sUpiEc6LSyj0jvQt4y5d1u+/dDXzmbiPMjb1PBe8rVNbvpiqycdquStvGuSWmPeJ+poNwLxQ4hu1GuNuurLLiBcCtMroQp83lIgBx+qsMEZEYEQl3q6bqAYWlh7dE5K0Sq3obp+bhNI6U5CJxSoJ7gTwROZvi7a++CqcS+x5SScP1KEdXPRWjqqtwDl5/8XGdn+H8U3xWzll/R2AWTjXBT8ArqlqlMy6vOH9UpzG2Is2Bj3EOQmtxzuYLS0SXA+1wzsI+w6l/L1lnXLi9bJykeCZe7TZuFcRZOFUQu3CK/k/iVl2Usp7DOJcQFl7pA85nslVVk8vYh8J/+H0i4ms9f6W5Z3azSmt/UNVfcZLEJrcariXwD5x655XAKpz9KrPjmKqux/kNzgI2AD+WWOQqYIuIpOE0do8FEJG24lwp1baMVUfjfCfpwCKcz/MBr/mvAV+4Ma7GOTN/rXCmOB0pS+1Y51avDcW57PRLnN/ROpyrfEq9Sq7E+0v93VRRAaUfl5YBLbxO9JJwzs534RxEb3K/v6oKo2oneWXFW0RVv8GpFn9bRM7HqS57BmcfUoBbgIvU7TuCU01Wsm/KJzjVR98Wlobc/8s/4ZwAHcC5aGRqWXGIc6XbvSWmFV4Z6vO+i29VnkZEfsO56qXUA64x5tiJSGFivKlkiUVE3sQ52dsM/E9VW1fTNsNxkt6vqvrXSr73MaAncIFbejzWWCJxSqs9qmN9PmyvH84JSHyJts8yhWJJo9q5xUqllMsjjTHV6lmcK532ytHDiLyEUxVTbdwrwJJwGqQnVmEVrwMxwC7xsbNjeVQ1R1U7ByhhTMGpoXjA14QBVtKokDhDd3QBrlLVGUEOx5haT5ze8dVW0ihjG21xqutK08W9Eq9W8ntJQ0RGitODeGNpl6GKyOkiskxE8kTk4hLzxonIBvdR3nXMfqOqQ1S1qSUMY2oGVZ3jz4ThbmObFh+SxvtRaxMG+Lmk4TayrMcZeG0HTqevy1X1F69l2uFcU343MFVVP3anN8RpgOyDUzW0FOjtdbWIMcaYAKvOAdZK0w/YWHhVgIi8j9MRpihpFF7zLSIlW+9HAN+oMxYUIvINztg371GGxo0ba7t27aoxfGOMCX1Lly5NUdWyRtMoxt9JoxVOx7BCO3Cu1a/qe4/qtCciN+CO+dO2bVuWLFlStUiNMaaWEpGtFS/lOO6vnlLViaraR1X7NGniU6I0xhhTRf5OGjtxOqoUak3ZQ1hU53uNMcb4gb+TxmKcYbIT3U4rYyinx2IJM3DGw2kgIg1weiTbFUzGGBNEfk0a6gxLfCvOwX4t8KGqrhGRR0XkAgB3DKMdOMNsvyYia9z37se50cxi9/FoYaO4McaY4Aipzn19+vRRawg3xpjKEZGlqurTAJXHfUO4McaYwPH3JbfGmEoqKChg946tpOzdQ56Gg2934TTmaKqESx6NmzSjResEPJ5jLydY0giCX3/9lV27djFsWGm36ja13W9rVyI7FnPSimeIzNyDHNPtKUxtpuIhJ6YZm7veRtK2TpwyYChhYWEVv7EcVj0VBJ07d+aMM47lpn4mlKVl5tJ+4b1EZey2hGGOiWgBURm7OXHpw2hUHHOnTqag4Nh+U5Y0gqBOHeduoNnZVblRmAl54sFTkFPxcsb4yFOQg3jCWLnwOw7sTar4DeWtq5piMpUwYcIEALZtq9WDZRpjAszjCSMjPfXY1lFNsZhKSElJAWDGDOuraIwJLLXqqeNPYbKIjo4OciTGmGC6ZkomZ046HOwwKsWungqCnBynvnr/fuvgbkLPNVMyeXtFLn/5XST/Gn7kxGhHWgFtnkvnu3F1GNIuHHkkrWheTDi0i/dwXa8I7hwQVTR9y8ECEp9PL3odFwUnNfZw78AoRp0U4XNM3tsqFBUGWffHAdBuwiG2pirPj4ziT/2jii3356+zmLAwhzMSw5h1dV0AClR59qcc3lqey+aDBUR4ICHew/knhvOPYaF9MmhJIwgKr16YP39+kCMxxj+iw+GFhTnc0jeShPiyKzReOjuai7qEk5kLM3/L49bpWcRGCjf0jiy23OdjYujXKowDmcqT83K46MNMfvyDcGpr3w9hhdsqVLL3S9v6wuvLcosljaw8ZdLKXBLqF1/60bnZPL8whxfPjmZA63Cy8pTVyQUs2JHvczzHK6ueCoLCoVvmzZsX5EiM8Y/ftQnj5OYe7p2dVe5y9aOheayHxAYebuwTSY9mHmb8lnfUcg1jhOaxHjo3CeM/50cTGQaf/5rHnC15hD2axvbU4vX0k1bkUP+JNA7nHBkmqXBbhY9mscUPf2O6RrDpQAELdxzZ/se/5NIgGga3K56cpvyax/iekVzZI5ITGnro2jSMy7pF8NzI4qWMt5fn0OXldCL/nkbrZw9x/+ws8gpKH7rpm9+cfdmRVnxfPlidS51/ppGW7bxvT3oB10zJpMlTh6j3eBqnvXGY77ce/Zn5i5U0gqAwaRw6dCjIkZjjyZC3jq77vrRrBH/sG0lGrnLOuxlHzb/mlAiuOSWSlIwCLv4w86j5N/eJ5LJuEWxPLeCqz46eP+eaulWKVYCnh0cz+K0M/nxqPn1alt+hTFWZvTmftXsLOLFR+eey4R6ICIPcAhjSLpyODT288XMuDw05UkL4z7JcrugWQd1I33vT14sSxnSL4D/LcunvlmAmLs3lul6R/JpS/EDeop4wd2seO9MKaBVXerxfrs/lD1Oz+MfQKC7qEs7Puwu46ctMBPh7KVVYZ7QPo0Ws8O7KXP5v4JF9eXtFLqNPCicuSsjMVYa+nUHnJh6mj61DfLTwwepchr+TwfIb69K5ybF13POFlTSCoGPHjoDTTyOUBow0xtughHBGnRTO3TPLLm1cNzWL2MfSiPzHIc58JwMRuKN/ZJnLZ+Upj8zNJi0bzmzvHNhv6B3BG8tzKHD/l35NyefHbflHVXEVbqvw8fe5R/eTuqF3JO+vzuVQtvJrSj4LduRz7SlHt508NyKazFxo81w6nV5KZ9yUTN5dmVusFPHEvBwu6hzO3wZFcWIjpyTy8OAonv4ph5z8o//vPSJc2SOCd1bmFk3bk17AzN/yGHeysy8frMklLVv54OIY+rQMo0NDD/edHsVpbcJ4bWnuUev0BytpBMHo0aOZPHkyAFu2bCExMTHIEZnjQXln/XUipNz5jet4yp3fpn7586vqyTOj6PrKYaauy6VXi6PPgv85zGnQTkov4J5Z2VzSJbzoLN/bWe9k4BHIzIMG0cJzI6IY2cFZbtzJEdw3O5sZG/M4u2MEry/LpXcLDz1LbK9wW4UaxhxdCunXyjkQv7c6l19TCji/U/hR1VgAJzUOY9XNdVmeVMCP2/KZvyOP677I5LkFHn64ti4xEcKa5Hwu61q8RDG4XThZedn8tr+g1FLBuJMjeHJeDst259OrRRjvrsqlaV3hzPbOsot35pOUrsQ/UbyWIjsfYiICM0aZJY0gW7JkiSUNE7JObBTGjb0j+L9Z2UwfW+eo+c1ihQ4NPXRo6GHKGA8nvphOzxZhnJ5Q/ND05qgYercMIz7aSYDeGtXxcHEXp1rpjPbhTFqRyz+GFb8CyntbFbmhdyT/XpLD9lTl3QtjylxOROjZIoyeLcK4rX8kP27LY9CbGXy4Jpdxp5RdWipP5yZh9GnpYdIKJ8lOWpHLlT0iCPM4CaFAoXMTD59ddvRnWcf3i8mOiVVPBcFzzz1X9Dwq6ugftzGh5KHBUew6VMDEpeUPjdK4jodb+kbyp+lZR1XbtopzDvglE0ahG3tH8MX6PF5bkktmnnJ5t6ofQa/sEcGGfQXUi4LhJ/jeRtC5sRNb8mEn9q5Nw45qoJ67JY+YcDihnOQ17uRI3ludy7Ld+azYU8DVJx/Zlz4tw9h0oIC4KIqSbeGjZb3AHM4taQTB4cNHGjR3794dxEiM8b8mdT3cc1oUExZUPJ7Wrf2cRuf3V1fuaqCBbcPp1MjD3d9kMaZrBPWiql5VExcl7LyzHitvisVTxrD0F32YwTPzs/lpex5bDxYwf3seV32WRYQHzj3RKSX9bWAkn6zN44kfs1m/L58P1+Ty8Nxs7hoQSWRY2fFd3i2cA5nK+KmZ9GrhoVvTI4lrbI8IEuM9nDs5g5m/5bHloHO11+M/ZDPl18C0aVjSCALvs6jZs2cHMRJjAuPPAyJpXKfiA3mzWA9XnxzBg3Oyy7w0tSzX94ogJ5+jGsCron60lJt4Rp4Qzte/5XHhh5mc+FI6l3yUSWQYzL2mDl3ctopzOkbwxgXRvL0il26vHObPM7L4Y5/IYld5laZRHQ/nnhjO8qQCru5RvMQUHS7MvaYOfVqEce3nmZz4YjoXfpjJol35JNQPzOHcbvcaBN27d2f16tUAtG3blq1btwY5IlOTLF26lN5f2L1WKuuv32TxzaY8fr4xNtih1EhLz5/NvE9eZfS1d9C2Y9di8ypzu1drCA+Cwh7hkZGRHDhwIMjRGHN8S81S1u9z2kxeODu0h/CoCax6Kgi6d+8OQFxcHBkZR3fIMsb4btT7GZz+1mF+3zmCK3sE6BKiWsySRhCcf/75ADRt2pT8/PxiDePGmMqZc01dMu+L481RMWU2XJvqY0kjCArbkRISEgCnDtsYY44HljSC4LHHHgPg7LPPBiAvL3CDjRljzLGwpBEEubnO9dRdunQBYPPmzcEMxxhjfGZJIwgKq6fatGmDx+Phyy+/DHJExhjjG0saQVCYNMLDw/F4PNamYYw5bljSCILCpCEixMXFsW/fviBHZIwxvrGkEQR9+jgdL0WEZs2aWV8NY0LA32Zl0ezpQ8gjaby1vOJxto5X1iM8CM4++2w++ugjRITExETWrl3L1q1biy7BNeZ4ds2UTN5ecfTgee9dFMOvKfk8MjeHns09LCsx3MeKpHxOec3ps7T9z7G0jvMUDTe++fZY2pVzr/FgW7gjjyfm5TDlshj6tw6jfpSQV6A8+1MO//05l60HC2hT38Md/SO5pV/xsbHW78vntulZ/LA1nzoRwsVdwnnmrOhK3XUwkCxpBIH3JbZdunThq6++YsmSJZY0TMgY1DaMDy8pfi+K+GjhiR/zaVJH+DWloOhGQ4VeW5pDQn1ha6p/xsPLyddyR5c9Fhv2F+ARit3k6b5vs5i4LJeJ50VzcvMwftqezw3TnIENr3cHVUzPUc6YlEGPZmHMH1+X/ZnKHz7P5GBWJu9ffPQ9M2qCmpu6Q9ijjz4KONVTl1xyCQAej30VJnREhkHzWE+xR3S4c8COi8K5aZLX/TUycpV3V+UyvueRs/AtBwsY9KZTdZv4fDrySFrRfdJVlafnZ9P++UNE/j2NE144xIQFxW/f2m7CIe6fncUfv8yk0b8OMehN573ySBr/W1m8+ujMSYe5ZsqRe6R//msuPV9Lp84/04h/Io1+/0nn5935pe7rNVMyueqzLArUWbc8kgY49/a+a0Akv+8cQfsGHsb2iOC6npH884cjcU5elUtKhjL5whhOaR7GsMRwXj4nmg/W5LH5QEGp2ws2v5c0RGQk8DwQBryuqk+UmB8FTAJ6A/uAy1R1i4hEAK8Dvdw4J6nq4/6ONxC8G8I7dOgAWF8NU747vs5ieVLpBy1/OqV5GBNGVv8ggDf0juDcyRk8MyKaOhHC+6tzaVnPw6CEIyWPNnHC52NiGPV+Jouuq0ub+lJUUnhlcS4PfJfN8yOjGdoujG8353PH11nUixTG9zqSeF5YmMOdAyL5aXwd8nw8BielF3DJR5n8Y1gUl3SJICtP+TmpgPAyzuueHxlNz+Ye7pqZzY47j1S5ZeVBdIkjbEwEbE1Vth4sICHew7zt+QxoHUb96CMloLNOCMcjMG97HokNjn2Y9+rm16QhImHAy8BwYAewWESmquovXouNBw6oagcRGQM8CVwGXAJEqWp3EakD/CIi76nqFn/GHAjeSaNBgwZERUXx0UcfceeddwY5MmOqx5wt+cQ+llb0ulWch3W3HjmgDmwbTus4Dx+szuXanpFMXJrL9b2KDzYY5pGi+3g3qSs097pX9xPzsrmtX2TRvTM6NgpjXUoB//whu1jS6NsqjIeHVC7p7T6k5BbApV0jitpRSrufd6H60VJ00PeO8eyO4bywMIczEsPp1tTDop35vPGz09az65CTNHYfKqB5bPEqs4gwZ793H6qZt63wd0mjH7BRVTcBiMj7wCjAO2mMAh52n38MvCQiAihQV0TCgRggB0gjBHjfw0REiIiIYN26dUGMyNR0/jjb96f+rcN4e/SRNo3SztKv7+Xc17t3yzCWJ+Uz7YoYVidXXBxIy1Z2pCmnJxQ/kA9uF8bzC3PIyFXqRDgH4n4tfb9da6EezTyMOCGMbq+kM/yEcIYkhHFh5wjaVPImR8+PjOamaZmc8tphBGhZTxjfM4In5uXgqZlt3D7xd9JoBWz3er0D6F/WMqqaJyKpQCOcBDIK2A3UAf6sqvtLbkBEbgBuAOeGRscD75IGQOPGjdm+fXt5bzHmuBIT7tzDujxXnxzJ377N5s4ZWfy+c7h7/+/qrccv7QokAUreey7Xa7NhHmH62Dos3lXArE15fLI2j3u+zeajS2I470Tfh15vGCN8eEkdcvKV5MNKy3rCq0uckkb7Bs5n06Keh+2pxfc5N1/Zn6m0qFczM0tNbn3tB+QDLYFE4C4RaV9yIVWdqKp9VLVPkyZNAh1jlZx22mnAkaTRrl078vPz7YZMplZpGCNc3CWCbzfnc32v0uvuC9sw8r2Oq3FRQus44futxdt45m7JJ7GBFJUyytK0rrDLq+onO0/5ZW/xA7eI0K9VGPcOiuL7a+syOCGMN5dX7R7ckWFC6zgPHhHeW53L6QlhNKnrHHpPaxPGTzvyScs+Es83m/IoUDitTc28uNXfSWMn0MbrdWt3WqnLuFVR9XEaxK8AvlbVXFVNBuYBPt2OsKYbPnw4cCRpFA5c+MMPPwQtJmOC4T/nR7P3L7EMSyz9AJlQX/AIfLUhj+TDBaRmOQfXvw2M4sVFOfxnaQ4b9uXz2pIc/r0kh3sHln//bYAz24fz6tIcftqex+rkfK75PJOc/CMH7fnb8/j73GwW7shjW2oB327KY+WeAro0PnK4POmldF5aVH4HvsU78/loTS6/7S/gp+15XPxhBsuT8nnBq6rxiu4RNK4jXPFJJiuS8vlucx63fJXFZV3DSWxQM8/p/Z3KFgMdRSQRJzmMwUkG3qYC44CfgIuB2aqqIrINGAa8IyJ1gVOBCX6ONyAKb7pUmDT69+/PK6+8Yu0aptaJDpeiS3FL0yzWw+NnRPHEvGzumJHFoLZhzLmmLjf3ieBwjvLYj9n88SulTZzwxJlRxRrBy/L0WVFc/4Uy4n8Z1I8W7h0Yxd7DR5JG/Sjhpx35vLw4hwNZSvNYYWz3CB4YfCQhrdtXQEpG+VVp2fnKI3Oz+e1AAZFhcHpCOPP/UJfuzY60s8RGCrOuqsNt07MY8N/DxEQIF3cO59kRNbcNS7Rk5V51b0DkHJyDfRjwhqr+U0QeBZao6lQRiQbeAXoC+4ExqrpJRGKBN4EuONWQb6rqU+Vtq0+fPrpkyRI/7k31aNy4Mfv27WP37t00b96c/fv306hRI55++mnuuuuuYIdngmzp0qX0/mJYsMMwIWbp+bOZ98mrjL72Dtp27FpsnogsVVWfanL8Xmmmql8BX5WY9qDX8yycy2tLvi+9tOmhoGRDeMOGDWnUqBHr168PZljGGFOhmllpFuJKJg1whkn//PPPgxWSMcb4xJJGEJSWNOrUqUNKSkqwQjLGGJ9Y0giC0pJG+/btyc/PZ8+ePcEKyxhjKmRJIwjOPPPMo6b17NkTgFmzZgU6HFPTqFIgNfMafXN8KpBw0OrpOGlJIwhOP/10oHhJY+DAgQDMnz8/KDGZmiNS8siI7xTsMEwIyYjvREFW9YzCZEkjCAp7fnsnjaFDhxabZ2qvVu068tuAJ0lv0NVKHOaYFEg46Q26sr7vP0jatgkAj6fy43F5s19kEDz1lNPdxDtpxMXFkZCQgL/7zZiar2GjRkAX1mb/HYmJQ8TO7UwVaQEFWWkkbfmNA0nb0IICYuMbHtMqLWkEQVmJoWvXrqxZsybA0ZiaqGGjRuS178DHE/9FQX4+UXXqFjvJMMZXqkp+Xi4Zh1IZcsFY6jc8tjH6LGkEQWlXTwFkZmayatUqDh8+TN26dYMRmqlBmrZK4JIb/4/l87/l4L5ktJoaMk3tU6duHCd07clJPQcc88lHlZKGiNyvqv84pi3XYmUljU6dOvHdd98xe/Zszj///GCEZmqYJi3bMvzia4MdhjFFfKosFZEwEWkgIo1F5AzgXhHpJyIN3bGjTCWUlTQGDRoEwJw5cwIdkjHG+KTCpCEizwApwFJgIfAgcDXwqjtth4iUO5CgKe68884Djk4aZ511FgA///xzwGMyxhhf+FLSOBtooqrtVfUEVR2sqh+rai9VTcS5SdJZ/g0ztJx66qmlTm/cuDHh4eFs3LgxwBEZY4xvfEkaf1TVvLJmqmoOcHv1hRT6kpKSgKNLGgAdO3YMdDjGGOOzCpOGqs6pjmXMERMmTABKTxrnnnsuycnJ5OfnHzXPGGOC7Zh6DYlIr+oKpDYqLWl06dKF7OxsVq1aFYSIjDGmfMfa1fTmaomilinr6ikAj8f5St57772AxmSMMb44pqShqtdXVyC1SXlDhZx99tkALFq0KFDhGGOMzyrs3FdRFZSqLqu+cGqH8koaTZs2JSIignXr1gU6LGOMqZAvPcKfKWeeAsOqKZZa46KLLuKTTz4pszt/8+bN2b17d4CjMsaYilWYNFR1aCACqU169epVbtLo0qUL27dvZ+vWrSQkJAQ4OmOMKVul2jREpJuIXCoiVxc+/BVYqFJVtm7dWu4yo0ePBmD9+vUBiMgYY3znc9IQkYeAF93HUOBfwAV+iitk5eXlMXHiRKD0Ng2gaLBCa9cwxtQ0lSlpXAycASSp6rXAyUB9v0QVwryvnCorabRs2ZL4+HimT58eqLCMMcYnlUkameoM6J8nInFAMtDGP2GFroKCI/dEKCtpiAhhYWF89913gQrLGGN8UpmksURE4oH/4Ixuuwz4yR9BhTJfkgZAhw4dyMzMJD09PRBhGWOMT3wZGv1lETlNVf+oqgdV9VVgODDOraYyleDrPcD79esHwLRp0/wZjjHGVIovJY31wNMiskVE/iUiPVV1i6qu9HdwoSgqKopLLrkEKL+kUXjPjRkzZgQkLmOM8YUvo9w+r6oDgMHAPuANEflVRB4SkRP9HmGICQ8Pp0uXLkD5SWPoUKd7zJIlSwISlzHG+MLnNg1V3aqqT6pqT+ByYDSw1l+Bharc3FyfbrIUERFB7969iYqKCkBUxhjjm8r00wgXkfNF5F1gOrAOuNBvkYWo9PR03n333XJLGYWGDh3K6tWrycnJCUBkxhhTMV8awoeLyBvADuB64EvgBFUdo6qf+zvAUON99VRF2rdvT3Z2Nl9++aUfIzLGGN/5UtL4GzAf6KyqF6jqZFU97OsGRGSkiKwTkY0ick8p86NE5AN3/kIRaec1r4eI/CQia0RklYhE+7rdmqoySaNr164AfPTRR/4KxxhjKsWXAQurPIqtiIQBL+NcorsDWCwiU1X1F6/FxgMHVLWDiIwBngQuE5Fw4H/AVaq6QkQaAblVjaWmKG9Y9JIGDhyIiLB48WJ/h2WMMT6p8k2YRGSt+7i1nMX6ARtVdZOq5gDvA6NKLDMKeNt9/jFwhjhH1LOAlaq6AkBV96nqcX/j7MqUNDweD40bN2b79u1+jMgYY3xX5aShqp2BgcDmchZrBXgf8Xa400pdRlXzgFSgEXAioCIyQ0SWichfqxprTdKgQQMuu+yyotu6VqR79+5kZ2dXODKuMcYEQqWThog0EpHfi0gv9+zfX6204ThJaaz79/cickYp8dwgIktEZMnevXv9FEr1iYqKIjEx0afqKYARI0YAMHXqVH+GZYwxPvHl6qlpItLNfd4CWA38AfifiNxRwdt3UnxQw9butFKXcdsx6uN0ItwBfK+qKaqaAXwFHHXrWVWdqKp9VLVPkyZNKtqdoMvIyODXX3/1OWlcd911iAgHDhzwc2TGGFMxX0oaiaq62n1+LfCNqp4P9MdJHuVZDHQUkUQRiQTGACVPmacC49znFwOz1WktngF0F5E6bjIZDPzCcS4lJYUpU6b43LbRsGFDunfvzrx58/wcmTHGVMyXpOF9xdIZOGf8qOohoNwjn9tGcStOAlgLfKiqa0TkUREpvIHTf4FGIrIRuBO4x33vAeBZnMSzHFjmx6qwgClMFr6WNACaNm3Kt99+S27ucX/xmDHmOFfhJbfAdhG5Dae6qBfwNYCIxAARFb1ZVb/CTTRe0x70ep4FXFLGe/+Hc9ltyPB1lFtv7du3Z9asWXzxxRdceKF1wjfGBI8vJY3xQFfgGuAyVT3oTj8VeNM/YYWuqpQ0Lr30UgA+/fRTv8RkjDG+8qWk0UJVbyo5UVW/E5GT/BBTSKtMP41CQ4YMwePxMH/+fD9EZIwxvvOlpPGZiPQuOVFEHsEZi8pUQuvWrRkzZgzh4b7ka0dYWBgtWrRg27ZtfozMGGMq5kvSuAT4SEQGAIjjVeB0YIgfYwtJMTExtGzZ0ufOfYWGDBlCfn4+GzZs8FNkxhhTMV9uwrQU594Z/xORkThDfTQBRqpqmn/DCz0HDx5kzZo1lW4Q//Of/wzYTZmMMcHlS+e+hjhXTo3DuZIpF7gRqOvOM5Wwc+dOZsyYUem2jVNOOYX4+HimT5/up8iMMaZivtSRLAWWAJOAQzid+hZ5TTeVUJWGcHDaNeLj45k8eXKVLts1xpjq4Ev1VKKqtncfiV6vE1W1fSCCDCWVGRq9pL59+5Kfn8/cuXOrOyxjjPFJlUe5NVVT1ZIGwFVXXQXAm29a9xhjTHD40qZRYSW6L8sYR1U69xU699xz8Xg8VtIwxgSNL50FmonIE0AyzlhTO3DuE34Z0BBoCsT6LcIQ07lzZ8aMGcPXX39d6fd6PB7atGnD9u3byc/PJywszA8RGmNM2XypnhoN7MJJDA2AK3EGEBwO1MMZxvwav0QXgmJiYmjcuHGl+2kUGjNmDAUFBaxcubKaIzPGmIr5co/wbcALha/dW7FuVtWx/gwsVCUlJbFixYoqXwF1yy238OSTTzJnzhx69uxZzdEZY0z5Kn26697rwq6aqqItW7bwww8/kJ9ftdudt2nThnbt2jFp0qRqjswYYypWpToSVT3qEiAR+eTYwwl9x9IQXqhRo0YsX76cPXv2VFdYxhjjk+q85NZKHz6ojo55Y8aMAeCFF16oYEljjKle1Zk0rJuyD6qjpHHTTc5I9Z999lm1xGSMMb6yzn0BVh1JIzY2lpYtW7J+/fpj6ixojDGVVZ1Jo+pHwVqkf//+XH755UREVHin3HINHz6c/Px8G8DQGBNQ1Zk0/q8a1xWyoqOjiY2NrXI/jUIPPujcZn358uXVEJUxxvimSkcuEbmm5DRVnXnM0dQCmzdvZvny5cdcrdS+fXv69u3LtGnTqikyY4ypWKWThoicALwoIs38EE/I27hxI4sXL66WtoiePXuyYMECli1bVg2RGWNMxXwZsPDPIrJNRFJF5BAwG3gIWCMiaSKyW0T+7PdIQ0R1Nlyfd955ADzxxBPVtk5jjCmPLwMW3gj0VdWSPcmeBRCRpsBM4Llqji0kVcfVU4XOO+88IiMjmTnTagaNMYHhS/XUI6UkjCKqmgw8WX0hhbZjuQlTSSJC//79SU1NZfXq1ce8PmOMqYgvd+57rzqWMY7qLGkA3H777QD885//rJb1GWNMeSqsnhKRYao6W0QuLG2+qn5a/WGFrrPOOosxY8Ywf/78alnf73//e+Li4lizZk21rM8YY8rjS/XUYPfv+aU8zvNTXCErMjKSqKioalufx+Ph9ttvZ82aNSQnJ1fbeo0xpjS+VE895P69tpTHH/wfYmhZs2YNS5curdarqC655BIKCgpsAENjjN/53E9DRKJE5AoRuVdEHix8+DO4ULRu3TpWr15drUmjW7duxMTEMGHChGpbpzHGlKYynfs+B0YBecBhr4ephOpuCC9c1/Dhwzl8+DBfffVVta3XGGNK8qWfRqHWqjrSb5HUEv5IGgB///vfmTp1Kv/4xz8455xzqnXdxhhTqDIljfki0t1vkdQS1dlPw1uPHj1o0qQJixYtIjc3t1rXbYwxhXwZRmSViKwEBgLLRGSdiKz0ml7R+0e679koIveUMj9KRD5w5y8UkXYl5rcVkXQRubsS+1Vj+aukAXD55ZeTn5/Pq6++Wu3rNsYY8K2kcR7O5bVnAx2Aszhyue355b1RRMKAl933dgEuF5EuJRYbDxxQ1Q44Q5GU7F3+LBAyN4249NJLueyyy4iMjKz2dT/88MPExcWxYMGCal+3McaAb5fcbi18AI1wGsMvABq508rTD9ioqptUNQd4332/t1HA2+7zj4EzxD0NF5HRwGYgZHquhYWF4fF4/FLSaNCgAVdeeSWffPIJ+/btq/b1G2NMZS65fRDn4N4IaAy8KSL3V/C2VsB2r9c73GmlLqOqeUAq0EhEYnFu7PRIBXHdICJLRGTJ3r17fd2doFm0aBFLliwhPz/fL+sfP3482dnZXHPNNX5ZvzGmdqtMQ/hYnNFuH3I7/J0KXOWfsAB4GHhOVdPLW0hVJ6pqH1Xt06RJEz+GUz3Wrl3Lhg0bihrEq1uvXr2oX78+X3/9NXl5eX7ZhjGm9qpM0tgFRHu9jgJ2VvCenUAbr9etS3lP0TIiEg7UB/YB/YF/icgW4A7gXhG5tRLx1kj+bAgvdPXVV5OXl8fjjz/ut20YY2qnyiSNVJwbL70lIm8Cq4GDIvKCiJQ1fsVioKOIJIpIJDAGmFpimanAOPf5xcBsdQxS1Xaq2g6YADymqi9VIt4aKRBJ47HHHsPj8fDyyy/7bRvGmNqpMp37PnMfheZU9AZVzXNLBzOAMOANVV0jIo8CS1R1KvBf4B0R2Qjsx0ksIauwWsrjqdLt2X0SGxvLgAEDmDdvHnPnzmXw4MEVv8kYY3zgc9JQ1be9X4tIPWA4MEJVbyznfV8BX5WY9qDX8yzgkgq2/bCvcR4P/FnKKPT888/Tr18//ve//1nSMMZUm0qd7orIySLyfyIyG/gS6A684ZfIQtR1113HhRdeSEREhF+307t3b6666iomT55sl98aY6qNLz3CLxWRN0XkZ+BOYBtwiaqerqqPqOpCv0cZggJR2rj77rvJyMjgxhvLLAgaY0yl+FI91RF4SVWX+juY2uDbb79l8eLFxMbG+n1b3bp1o2nTpnz66afs37+fhg0b+n2bxpjQ5kv11FsVJQwRKdlhz5Thl19+Ydu2bX7rp1HSfffdh6ryxz/+MSDbM8aENl+Shi83aJh2rIHUFv4a5bYst956K7GxsXz88cekpaUFZJvGmNDlS9LYLiILRGSqiEwRkZdE5BQRedud9hOw0d+BhopA9NPw5vF4uP/++8nPz+e6664LyDaNMaHLlzaN83GukorH6WtxCk4fjb8C64AMYLk/ggtFgS5pAPzlL3/hmWeeYdasWWRmZhITExOwbRtjQosvo9yqqq5U1e9V9TtVfQ740h3zaa6qLlZVu+uPjyIjI4mIiPBr576SPB4PH374IQcOHOC1114L2HaNMaFHAtUgGwh9+vTRJUuWBDuMCl1wwQVs376dn3/+OaDbHTp0KEuXLmXdunW0aNEioNs2xtRcIrJUVfv4smzgTndNEVUNaPVUoauvvppDhw5x6aWXBnzbxpjQUJn7aTzsxzhqjU8//ZRglYauvfZaEhIS+PHHH1m40PpkGmMqz5ce4R4R+S/OUOjmGP36668kJSUFpaQBMHnyZMC5n7gxxlSWLyWNacB+Vf2bv4OpDYJx9ZS33/3ud/zud79j8+bNvPnmm0GJwRhz/PIlafQGPvV3ILVFYT+NYPrwww+JjIzkqaee8tttZ40xocmXpDEUmCgi/f0dTG0QiPtpVKRVq1ZMmjSJtWvX8u9//ztocRhjjj++9NP4BRgBPOX/cEJfXFwc0dHRQaueKnTppZcybNgw7rrrLpYutbEojTG+8el0V1V3Aef6OZZa4Y477mDw4MFBTxoiwn333UdOTg4XXHBBwAZQNMYc3ypTR9LRb1HUMjXlAD1s2DBGjBjBrl27uOOOO4IdjjHmOFCZpPGMiKwVkb+LSDe/RRTi3nzzTZYtWxb0kkahzz77jLp16/Liiy8GvIe6Meb443PSUNWhOI3ie4HXRGSViNzvt8hC1Pr169m3b1+NSRoxMTF8/PHHqCpnn312jSkFGWNqpkpdwqOqSar6AnATzsi2D/ojqFAW7H4apRk5ciTjx49nz549TJw4MdjhGGNqsMoMI9JZRB4WkVXAi8B8oLXfIgtRBQUFiEiNShoAEydOZPjw4dxxxx189913wQ7HGFNDVaak8QZwABihqkNU9d+qmuynuEJWTejcVxqPx8M777yDx+NhxIgR7Ny5M9ghGWNqIJ+Shoj0UNUBwGz38ltTRc2bN6du3bo1rqQB0KxZM5544glyc3MZMGBAjU1wxpjg8bWk8QcR6QiM92cwtcHdd99Nnz59amTSALjtttsYOXIk27dv58ILLwx2OMaYGsaXUW4fcpdbAHhExBq/j1Gw7qfhq2nTptGyZUs+//xzHn300WCHY4ypQXwZRuQRYBbwAfCNqtpR5Bg8++yzrFixokYnjbCwMJYtW0ZsbCxPPvkkq1atCnZIxpgawtfqqX6q+kegnz+DqQ02bdrEoUOHgh1GhZo1a8batWuJj4/nvPPOY+vWrcEOyRhTA/hSPTUKSAJQ1QdEZKGIbHIfF/s9whBTE/tplKV169ZMnTqV3bt307lzZ7Zt2xbskIwxQeZLSeOvwOder6OAvsAQ4GY/xBTSamo/jbL07t2bBx54gMzMTHr06MH+/fuDHZIxJoh8SRqRqrrd6/WPqrpPVbcBdf0UV8gqvIz1eEkaAA888AC33HILqampdOrUibS0tGCHZIwJEl+SRgPvF6p6q9fLJtUbTuhLTEwkNjY22GFU2ksvvcS4ceNISUmhY8eOZGZmBjskY0wQ+JI0ForI9SUnisiNwKKK3iwiI0VknYhsFJF7SpkfJSIfuPMXikg7d/pwEVnqDoy4VESG+RBrjXfPPffQuXPn46qkUeitt97iiiuuIDk5mYsuusgShzG1ULgPy/wZmCIiVwDL3Gm9cdo2Rpf3RhEJA14GhgM7gMUiMtW9G2Ch8cABVe0gImOAJ4HLgBTgfFXd5Q7FPgNo5fOe1XDHY9IAePfddxkyZAg33ngjAwcOZNq0abRo0SLYYRljAsSXfhrJqvo74O/AFvfxqKoOUNU9Fby9H7BRVTepag7wPjCqxDKjgLfd5x8DZ4iIqOrPXkOWrAFiRCTKl52qye6//37WrFlz3CYNgOuvv55///vfLFu2jA4dOrBy5cpgh2SMCZDK3E9jtqq+6D5m+/i2VoB3I/oOji4tFC2jqnlAKtCoxDIXActUNbvkBkTkBhFZIiJL9u7d62NYwbN9+3YyMzOP66QBcOONN/Lggw+SkZFBr169mDZtWrBDMsYEQKXupxEMItIVp8rqxtLmq+pEVe2jqn2aNKn57fI1fQiRynjkkUd4++23UVUuuOACJkyYEOyQjDF+5u+ksRNo4/W6tTut1GVEJByoD+xzX7cGPgOuVtXf/BxrQByPl9yW5+qrr2bOnDlERERw11138fLLL9vd/4wJYf5OGouBjiKSKCKRwBhgaollpgLj3OcX4wy/riISD3wJ3KOq8/wcZ8CEWtIAGDRoEOvXr2fEiBHceuutjB49muRku9WKMaHIr0nDbaO4FefKp7XAh6q6RkQeFZEL3MX+CzQSkY3AnUDhZbm3Ah2AB0Vkufto6s94A6F79+7Uq1cvpJIGQEJCAtOmTePhhx9m6tSptG3blunTpwc7LGNMNZNQqkro06ePLlmyJNhhVKhv3740bdqUL7/8Mtih+MXjjz/Offfdh6pyxRVXFN0R0BhTM4nIUlXt48uy9p8cBKGUqEvzt7/9jUWLFtGwYUMmT55MixYt2LBhQ7DDMsZUA0saAXbzzTezfv36kKueKqlPnz7s2bOHiy66iJSUFH73u9/x3nvvhXzCNCbUWdIIsD179pCbmxvySQMgPDycjz/+mFWrVpGYmMgVV1xBixYtmDNnTrBDM8ZUkSWNAAvFq6cq0qVLF3766Sduvvlm9uzZw9ChQxkxYgSpqanBDs0YU0mWNAIslDr3VUZYWBivvPIKixYtomXLlsycOZPGjRtz1113WZWVMccRSxoBVhtLGt769u3Lzp07efLJJwkPD+fZZ59l4MCBLFy4MNihGWN8YEkjwAYMGEBsbGytTRqF/vrXv5KamsqLL77Ib7/9xqmnnkrz5s159913gx2aMaYcljQC7N5776VFixa1PmkAREZGcuutt7JhwwbGjBlDcnIyV155JS1atGDy5MnBDs8YUwpLGkFQW9s1ylKvXj3ee+89tmzZwtChQ0lKSmLs2LE0a9aMGTNmWJuHMTWIJY0Au/jii9m8eXOww6iR2rZty+zZs9myZQsjR44kKyuLkSNH0qNHD26++Wa7N7kxNYAljQBLTU0lPz/fShrlSEhIYPr06SQnJ/PWW2+RkZHBq6++Snx8PP3792fWrFnBDtGYWsuSRoAVFBQgIpY0fBAVFcW4cePYsGEDjz/+OM2bN2fRokUMHz6chg0b8tprr5GdfdR9uYwxfmRJI8AK6+ctafjO4/Fwzz33sGvXLpYtW8aZZ55JRkYGN910Ey1atGD48OE8//zz5OXlBTtUY0KeJY0Aq+39NI5Vz549+eabb8jIyGDGjBmce+65fPvtt9xxxx1ER0fTr18/3n333aLP2RhTvSxpBNiIESOoW7dusMM47nk8Hs466yzeeecdkpKSuO2222jcuDGLFy/myiuvJD4+nj/96U/MmjWLjIyMYIdrTMiw+2kEQefOnenRowcffPBBsEMJORs3buTxxx9nw4YNLF68mKysLESEhIQELrjgAm699VY6duwY7DCNqVHsfho1nPXT8J8OHTrw3//+l++//56UlBRefPFFTjjhBLZu3coLL7zAiSeeSN26dbnpppuYM2eONaQbU0mWNAJs8ODB7Nixw5JGANStW7eox3lWVhYTJ05k6NChhIeH8/rrrzN06FDq1atHkyZNGDlyJK+88oqNvGtMBSxpBFh2drb1cA6CyMhIrr/+embPnk1qair79u1jypQpnHrqqaSlpTFjxgxuueUW4uPjadiwIffeey/Tp08nKSkp2KEbU6NY0ggwu3qqZqhfvz6jRo3i+++/Jzs7m19++YV77rmH3r17ExkZyVNPPcU555xDixYtiIyMLGoTefHFF9mzZ0+wwzcmaKwhPMD69u3LL7/8wujRo21E1xrs8OHDLFy4kKeffpoVK1awZ88e8vPzASfhd+3alV69epGamkr//v0555xz6N69Ox6PnYeZ409lGsLD/R2MKc5KGseHunXrMmzYMIYNG1Y0bfXq1Xz00Ufs2LGDpKQkZs6cSVJSEp9//jn33nsvIkL9+vUZNGgQo0eP5qSTTiIhIYFWrVoFcU+MqV6WNALs0ksv5emnn7akcRzq1q0b3bp1KzZt3bp1fPnll8ybN4/Vq1ezc+dOZs6cyRdffFG0jMfjoV69erRo0YITTzyRUaNGMWjQIBITEwkPt39Bc3yx6qkgaN++PaeddhrvvPNOsEMxflBQUMCmTZuYPXs2n376Kb/99ht79uwhPT292EUQ4eHhiAjx8fE0bdqUdu3a0alTJy644AK6d+9OgwYN7OTCBERlqqcsaQRYTk4OnTp1YtCgQUyaNCnY4ZgA27BhAzt37mTLli3MmjWLOXPmsH//frKyso66qi46OhqA+Ph4mjdvTps2bUhMTGTUqFG0b9+eFi1aEBUVFYzdMCHGkkYN1q1bNzZv3szFF1/M22+/HexwTA2ydetWFixYgMfjYfv27cyaNYulS5eSmppaZidEj8dDVFQUcXFxNG7cmJYtWzJ27Niiu0M2a9aMTp06WXIx5bKG8BrMRrk1ZUlISCAhIaHo9Z133ln0vKCggG3btrF8+XJiYmLYsWMHM2bMYOXKlezfv5/U1FSSk5NZs2YN33zzzVHrFhEiIyOpX78+Z5xxBs2aNWPv3r1ERkbStm1b2rdvX/Sw2xGb8ljSCDAbfdVUhcfjoV27drRr165o2vjx449a7uDBg6SkpJCUlMTs2bP55Zdf2L17N3v37uXgwYPk5eWxcOFC9uzZw+HDh0vdVnh4OA0bNiQzMxMRoU6dOtSrV4/4+Hjat2/PmWeeSaNGjUhOTqZx48a0bduWhIQEGjdubJcc1wKWNALMShrGn+Lj44mPj6dDhw4MHDiw3GX37t3L2rVr2bBhA5s2bWL37t3k5ubSunVr9u3bx6xZs9i/fz8HDhwgOTmZgoICFi9eXO5Am5GRkbRq1Yr69euzbds2IiMjqVOnDrGxscTFxXHSSSdx6qmnEhcXx2+//UazZs1o3rw5zZs3p2nTpjRt2pSIiIjq/lhMNbKkEWDWT8PUFE2aNKFJkyacfvrpPi2vqhw8eJBDhw6xb98+vvvuO7Zv387evXvZt28fqampxMbG0rRpU1JTU9m0aRMZGRns3bu3qGPkjz/+yOuvv17udmJiYoiNjeXAgQOEh4cTERFBZGQk0dHRdOrUic6dOxMeHs7q1aupV68ecXFxxMfH06BBA7p27UpiYiJhYWGkp6fToEEDGjRoQKNGjYouLDDHxpJGgN1www08/vjjljTMcUdEig7Cbdu2pWfPnpV6f15eHgcPHiQjI4N9+/bx008/kZycTEpKCgcOHCAtLY0mTZrQqFEjUlJSmDNnDllZWeTk5JCVlUV6ejpZWVmsWLGCQ4cOkZOTU+l9iI2NJT4+Ho/Hw549e4iLi6NTp07UqVOHOnXqcP/999O7d29Wr17NpEmTiImJoU6dOkV/zz//fJo1a8bOnTtZs2YNUVFRREdHExUVRVRUFO3btycqKors7Gzy8vKIiooKub44obU3x4G7776bCRMmWNIwtU54eDiNGzcGqFLSKSk9PZ29e/eSnJzM3r17SUlJoV69ekRHR7N7924WLFhAeno6hw8fJiMjg4yMDBITE4mJiWHXrl0cPHiQvXv30qVLF9LS0khKSiq6Sm3jxo28+OKLZGVlFdvmvHnzaNasGTNnzuQPf/jDUTGtXLmS7t27M3HiRP70pz8BR65wi4qKYvny5SQkJPCf//yHl156qWh6YeKZPHkycXFxfPrpp8ycOZPIyMhij/vvv5/w8HDmzp3L2rVriYyMpHv37vTt2/eYPsvKsKQRYPv377fGcGOqQWxsLLGxsSQmJpY6/7rrriv3/a+//jrXX389b731VrELDABGjx5NZmYmBQUFZGVlkZmZSUZGBk2aNAHgvPPOY968eWRlZZGdnV30t23btgCcdtpp/Otf/yo2Lzs7m7i4OAAaNGhA+/bti+ZlZGSwf//+ogsJfvnlF6ZMmUJOTk7RIzc3lwceeACA999/n1dffRVwrrILZNLwez8NERkJPA+EAa+r6hMl5kcBk4DewD7gMlXd4s77GzAeyAf+pKozytvW8dBPIyEhgb179zJ27Fj+85//BDscY2qtb775hrPOOos5c+YwePDgYIdToYKCgqKkkp6eTnp6Ojk5OdStW5dGjRod07prTD8NEQkDXgaGAzuAxSIyVVV/8VpsPHBAVTuIyBjgSeAyEekCjAG6Ai2BWSJyoqrm+zNmf7Orp4ypGQr7xGzbti3IkfjG+3LmwlJWMPi7eqofsFFVNwGIyPvAKMA7aYwCHnaffwy8JM4RdRTwvqpmA5tFZKO7vp+qO8jvvvuOMWPGlHrXthNOOAFwLk9MS0srNk9EaN++PUDR2ELewsLCioq9u3fvJiMjg+zsbETEkoYxQdamTRsA/vKXv/DEE09UsHTNd9VVV3HPPff4fTv+ThqtgO1er3cA/ctaRlXzRCQVaOROX1DivUeNMS0iNwA3AEX1iZVVr1492rZty+7du4+a16VLFwA2bdp01M13wsLCiuZHRESQkpJSbH5kZGTRfI/Hw4EDBxARTjjhBMaOHVulWI0x1SMmJoZHHnmEVatWBTuUatGsWbOAbOe4bwhX1YnARHDaNKqyjj59+rB48eJqjcsYU/M9+OCDwQ7huOPvPv87gTZer1u700pdRkTCgfo4DeK+vNcYY0wA+TtpLAY6ikiiiETiNGxPLbHMVGCc+/xiYLY6rcVTgTEiEiUiiUBHYJGf4zXGGFMOv1ZPuW0UtwIzcC65fUNV14jIo8ASVZ0K/Bd4x23o3o+TWHCX+xCn0TwPuOV4v3LKGGOOd3Y/DWOMqeUq00/DxjE2xhjjM0saxhhjfGZJwxhjjM8saRhjjPFZSDWEi8heYGuJyY2BlFIWrw1s32sn2/fa6Vj2PUFVm/iyYEgljdKIyBJfrwoINbbvtu+1je27//fdqqeMMcb4zJKGMcYYn9WGpDEx2AEEke177WT7XjsFZN9Dvk3DGGNM9akNJQ1jjDHVxJKGMcYYn4V00hCRkSKyTkQ2ioj/74MYRCLyhogki8hqr2kNReQbEdng/m0QzBj9RUTaiMh3IvKLiKwRkdvd6SG//yISLSKLRGSFu++PuNMTRWSh+9v/wL01QcgRkTAR+VlEprmva8t+bxGRVSKyXESWuNMC8nsP2aQhImHAy8DZQBfgchHpEtyo/OotYGSJafcA36pqR+Bb93UoygPuUtUuwKnALe53XRv2PxsYpqonA6cAI0XkVOBJ4DlV7QAcAMYHL0S/uh1Y6/W6tuw3wFBVPcWrb0ZAfu8hmzSAfsBGVd2kqjnA+8CoIMfkN6r6Pc79SLyNAt52n78NjA5kTIGiqrtVdZn7/BDOQaQVtWD/1ZHuvoxwHwoMAz52p4fkvotIa+Bc4HX3tVAL9rscAfm9h3LSaAVs93q9w51WmzRT1d3u8yQgMHeeDyIRaQf0BBZSS/bfraJZDiQD3wC/AQdVNc9dJFR/+xOAvwIF7utG1I79BufEYKaILBWRG9xpAfm9+/XOfabmUFUVkZC+vlpEYoFPgDtUNc058XSE8v67d7Q8RUTigc+Ak4Ibkf+JyHlAsqouFZEhQQ4nGAaq6k4RaQp8IyK/es/05+89lEsaO4E2Xq9bu9Nqkz0i0gLA/Zsc5Hj8RkQicBLGu6r6qTu51uw/gKoeBL4DBgDxIlJ4UhiKv/3TgAtEZAtO1fMw4HlCf78BUNWd7t9knBOFfgTo9x7KSWMx0NG9miIS597jU4McU6BNBca5z8cBnwcxFr9x67L/C6xV1We9ZoX8/otIE7eEgYjEAMNx2nS+Ay52Fwu5fVfVv6lqa1Vth/O/PVtVxxLi+w0gInVFpF7hc+AsYDUB+r2HdI9wETkHp94zDHhDVf8Z3Ij8R0TeA4bgDI+8B3gImAJ8CLTFGTL+UlUt2Vh+3BORgcAPwCqO1G/fi9OuEdL7LyI9cBo9w3BOAj9U1UdFpD3OGXhD4GfgSlXNDl6k/uNWT92tqufVhv129/Ez92U4MFlV/ykijQjA7z2kk4YxxpjqFcrVU8YYY6qZJQ1jjDE+s6RhjDHGZ5Y0jDHG+MyShjHGGJ9Z0jDGj0TkHhEZKyIPi4iKSAeveXe40/qUtw5jahJLGsb4gTg8wAhgpjt5FU5HtEKXAGsCHZsxx8KShjHVRETaufdvmYTTQ7cNEKmqe91FpuCOtCwiJwCpQEowYjWmqixpGFO9OgKvqGpXoDfOfQ0KpQHbRaQbTonjgyDEZ8wxsaRhTPXaqqoL3Ocjgekl5r+PkzBGc2QoCGOOG5Y0jKleh72e9wMWlZg/DbgK2KaqaQGLyphqYvfTMMYPRKQr8Kt7r4siqpohIv8HrA9OZMYcG0saxvjH2cDXpc1Q1fcDHIsx1cZGuTXGD0TkG+Bqr9tvGhMSLGkYY4zxmTWEG2OM8ZklDWOMMT6zpGGMMcZnljSMMcb4zJKGMcYYn1nSMMYY47P/B+d0qwTXXumjAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "\n", "import matplotlib.pyplot as plt\n", "import mpmath as mp\n", "import csv\n", "\n", "# Download torus_cuts.csv:\n", "URL = \"http://astro.phys.wvu.edu/zetienne/torus_cuts.csv\"\n", "outfile = \"torus_cuts.csv\"\n", "try:\n", " with open(outfile,\"w\") as file:\n", " file.write(urllib.request.urlopen(URL).read().decode(\"utf-8\"))\n", "except:\n", " try:\n", " with open(outfile,\"w\") as file:\n", " file.write(urllib.urlopen(URL).read().decode(\"utf-8\"))\n", " except:\n", " # If all else fails, hope wget does the job\n", " !wget -O $outfile $URL\n", "\n", "def file_reader(filename,list_of_cols,delim=\" \"):\n", " with open(filename) as file:\n", " reader = csv.reader(file, delimiter=delim)\n", " data = list(zip(*reader))\n", "# print(data)\n", " # data is a tuple of strings. Tuples are immutable, and we need to perform math on\n", " # the data, so here we convert tuple to lists of floats:\n", "# data_output = [[sp.sympify(0) for i in range(len(list_of_cols))] for j in range(len(data[0]))]\n", " data_output = [[sp.sympify(0) for i in range(len(data[0]))] for j in range(len(list_of_cols))]\n", " for i in range(len(data[0])):\n", " for j in range(len(list_of_cols)):\n", "# print(i,j,data[list_of_cols[j]][i])\n", " data_output[j][i] = float(data[list_of_cols[j]][i])\n", " return data_output\n", "\n", "NRPy_data_output = file_reader('out.txt', [0,7,8,9,10])\n", "std_data_output = file_reader('torus_cuts.csv',[0,4,1,3,2])\n", "\n", "\n", "ylabels = ['Lorentz Gamma_{KS}=G','G*v^r_{KS,Val.}','G*v^{\\\\theta}_{KS,Val.}','G*v^{\\phi}_{KS,Val.}']\n", "\n", "for i in range(len(ylabels)):\n", " # https://matplotlib.org/gallery/text_labels_and_annotations/legend.html#sphx-glr-gallery-text-labels-and-annotations-legend-py\n", " fig, ax = plt.subplots()\n", " plt.title(\"NRPy's FM solve with FMtorus.f90: \"+ylabels[i])\n", " plt.xlabel(\"r/M\")\n", " plt.ylabel(ylabels[i])\n", " ax.plot(NRPy_data_output[0], NRPy_data_output[i+1], 'k--', label='NRPyFMSolve')\n", " ax.plot(std_data_output[0], std_data_output[i+1], 'k-', label='FMtorus.f90')\n", " legend = ax.legend(loc='upper right', shadow=True, fontsize='x-large')\n", " legend.get_frame().set_facecolor('C1')\n", " plt.show()" ] } ], "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.9.1" } }, "nbformat": 4, "nbformat_minor": 2 }