{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Compare NIST energy levels data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Uncomment the next line if you have ipympl installed and want interactive plots\n", "# %matplotlib widget\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "from ryd_numerov.elements import Strontium88Singlet" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "element = Strontium88Singlet(use_nist_data=True, nist_n_max=60)\n", "element_without_nist = Strontium88Singlet(use_nist_data=False)\n", "\n", "energies_with_nist = []\n", "energies_without_nist = []\n", "labels = []\n", "n_list = []\n", "\n", "l_int2str = {0: \"s\", 1: \"p\", 2: \"d\", 3: \"f\", 4: \"g\", 5: \"h\", 6: \"i\", 7: \"j\"}\n", "for n in range(5, 25):\n", " for l in range(n):\n", " if not element.is_allowed_shell(n, l):\n", " continue\n", " for _j in np.arange(abs(l - element.s), l + element.s + 1):\n", " j = float(_j)\n", " if (n, l, j) not in element._nist_energy_levels: # noqa: SLF001\n", " continue\n", "\n", " labels.append(f\"{n}{l_int2str.get(l, ',' + str(l))}_{j:.1f}\")\n", " energies_with_nist.append(element.calc_energy(n, l, j, unit=\"hartree\"))\n", " energies_without_nist.append(element_without_nist.calc_energy(n, l, j, unit=\"hartree\"))\n", " n_list.append(n)\n", "\n", "energies_with_nist = np.array(energies_with_nist) + element.get_ionization_energy(unit=\"hartree\")\n", "energies_without_nist = np.array(energies_without_nist) + element_without_nist.get_ionization_energy(unit=\"hartree\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAGdCAYAAAD60sxaAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAiYdJREFUeJzt3Xd4U2X7wPFvdtO9Wzqg7L2XLEFFQHDgQMSFuF73wIkDcIIKbl7xRVH8uVHEASJbZcuUPQsFSksH3SNNcn5/pKRpyWg6KJD7c11cNCd3njzn9PTkznmWSlEUBSGEEEIIH6Ju6AoIIYQQQpxtkgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jrahK3AuslqtpKamEhQUhEqlaujqCCGEEKIaFEUhPz+fuLg41Gr393gkAXIiNTWVxMTEhq6GEEIIIWrg6NGjJCQkuI2RBMiJoKAgwHYAg4ODG7g2QgghhKiOvLw8EhMT7Z/j7kgC5MTpZq/g4GBJgIQQQojzTHW6r0gnaCGEEEL4HEmAhBBCCOFzJAESQgghhM+RBOgct7tNW/KXLm3oagghhBAXFOkEfRZlfPAhmTNmVNqmb9qU5r8vdBpfsmcPAMceehiAwMsuI3HGh05jTceOk/nRfylatx5zZiba6GhCrrqKyPv+g0qvd1kna2kpJ994g7wFC7GWlRHYrx+xkyaijYysyS4KIYQ4z+06kcuNH62lwGQB4PJ2MdzZryljZq1j26QhhBh1DVzDuqFSFEVp6Eqca/Ly8ggJCSE3N7fao8AyP/4f+UuWYDp0CJWfH6jVWE6dAovtBNI1aULARReRt2gRKApKaSmKSgUlJfYy9M2b0WTOHHvyceDyIZQdPWp/vtWG9Whc1Kfg7785+tDDUFpaabs6PJwWi/9AExjo9HXHJ04k78d59nraKquj1epVLt9LCCHEuenS6Ss5lFFof2zUafjo1m4Mah0NwJGsQm6fvZ4jWcX2mAC9hh/v70ubRrZr/sVvriAlu8j+/LZJQ+j80mL7z1UToEFvreRIViGnkwkVEBGoRwGKSi0kRQbwn4ubMbJrvP01VquVGz5aw+ajufZtvZuGo9eq+edwNgF6Ldd3T+Dpoa3RaqrfWOXN57c0gdWRzI8/pmTHDqxFRViys7FkZlZKKspOnCDnu++w5uZizctDKS2tlPwAmA4eYv/AQezp1JmDI66slPxoEhLcJiQqo/8ZyQ+ANTublHvucfoaS34+ed/PrZz8AJSVcfDKq6qz20IIIc4hjskPQHGZhTs//4fMghIURWHgWysrJT8AhSYLw9//2/74dPKjAhpHGMkrLnP7no7JD4ACZBaYyCowUVxmITmzgPHfb2XZ7nR7zG2fbqiU/ACsT87m7/2ZlJRZySo0MXtVMm8s2lP9nfeSJEB1RDGZ3Ad4ev40iwV1cBCmgwcrbz52jJx5P7l8Wc7337t8rmTbv063F23f7roaJ09iOnbMQ2WFEEKc66wKfLshhYMn893GfLMhhWbPLrBvU4CUrGJG/2+t2/IDDBoC9BqXz5eUWfHXa1m0I832XlYrqw9mAbYkyxWzVeGz1Ycxma1u37+mJAGqI5qI8Dory5KZ5XR78PArXL7GfOqU6wKtVsrST56x+dT8n93WI/O//3X7vBBCiHOLzsWnekp2Mb/9m+b2tV+sPYyzVCM1p8TJ1gqNQoz4ad2nEwWlZkL9bU1nv21PtW/31AfHbFVYd8j5Z2JtSQJUR7RR0VDPC6eq/fxcPmctLnb5HEDxv9vO2Fa4eLHb1xRt2VqtegkhhDg3BBqcj226rE0MOo37z6jjp9x/jgCYLc7vxuSWuG8mAxjVw7bG5rcbjnqIrGxDcrZX8dUlCVAd8WvVCnVQEGjraWCdh3LVLjo5n2bJzDxzY4n7rF4pLHT7vBBCiHOfn07F0A6x9GkW5jauzOq5qWnZnvQztgX6aalOK1WrGNv6XCnZnhMtR9VJrmpChsHXEUOzphSFhmItKKiX8jUehqUrRUVunzfn5Lp93mmZVovnoGqwuqubRoPaYKherFpd6S6YV7HFxeBqwKNKhdporFlsSQm4uWio/f1rFltaembn9BrGqoxG+7o4VpMJzOa6ifXzQ6W2fYdSTCaUuoo1GFBpNN7HlpWhlLm+UKr0elTlXyS8ijWb3fbxU+l0qHQ672MtFttgCFexWq19CguvYq1WFDdfbryJRatFfTpWUVDc3Wn2Jtabv3sfukbUBYOTpqjv7u0DwMl894mETq2m2GkjWIWqnawBik2ePyvUDjefTGbvPltc3dWqLUmA6khpylHKjh93e6LXhiUtDUtensuRYKWHD7t9vSYslLL0dFLuGEfcG1Mxdurk8T1Vxtr9YSqKQnGZhSPduruM8evXj5iPZtofH+nTD0qdX5D13brR6PM59scpl16GkpPjNFbXti1x31V0DD96xXCsac7bv7VNmxL/8y/2x8evux5zcrLTWHVsLImLl9gfp465mbLdu53GqkJDafxXxciKE3fehWnzZqexGPxo8s8/9ofpDzxIyerVzmOBJv9WdGA/+fh4ipcvdxmbsHYdmoAAADImPEfRggUuY+NXrEQbEQFA5suvUPjDDy5j4xYsRJdou6Wd/eY08r/8P5exsT/8iKFVKwBOfTCDvFn/cxkb8+VX+JWfnzmfzCb3/fdcxkbP+gRj794A5H31FafeeMNlbOR77xNwySUA5M+bR/bkyS5jI954g8ArhgNQ8PtCsp55xmVs+OTJBF13PQCFK1aQ+egjLmPDnnmG4FtuBaB4/XpO3nO3y9iQRx4l9G7b8yX//kv6rbe4jA2+517CHrbNF1a6bx9pN1zvMjbo1tsIf/ppAMqOHiV1xHCXsQE33EDkxEkAmLOyOH7JIJex/iNGEDVlKgCWwkKO9bnIZazx0kuJfrfi9yrXCBujTlOtRTzdaRYdQFp+5SS8c6Ltzs/Gw+770sSHGslLc91R2pUbuifw4bL95JS4/rKidtgvg9a7fWzXqH6mZJEEqI6U7NplS350OnD8ZqnRuP127g11UJDL5zz1AQq65FKUsjJMyckcHn0TLVau8Ph+2qgor+vo6FR+Ed1eX8nvbmJOrt/IJRP/sD9eWFriclRAztbtXOYQuyAnx2UbbsGe/bRziP0tLR1XYxSKk49Uiv0l+QiupvkqS0uvFDt/z34MLmItOTmVYn/cuh1XKaVSWlIp9rv1G3H3J+8Y+9Xfa3HXBb/fcz+RF2L7XX6+9C9i3MQOf2Eex6IaAzBr8TIS3MTeMflHtsS3BeDDpYto7ib20Vd/4K8k24fcOyt+o42b2BemfMXClicAeOOv+bhL1ae+9RVz2+UAMGnNPFx/5MJH737D7BW2D4Zn1//AQDexn8/4nv+utZ0xD2z+HneTQvz08Q9M3Wprgr5z24+MchP7+2fzeGmv7XcxatcC7nQT+/eX83kmxZZgDt+/nIfdxG74/jcez7AlmBcf3sQEN7Fb5y/ioYLOAHQ9vpvX3cTuXbiMi822o5qQkcIsN7HJS/9ioM52XgbnZvCdm9jjf69lkMM5LNcIm83PDSI8OMBFdPVsP5Z3xrY+U5Zy74BmZBe6H418e98kJsxzPToY4PPVyWw/nsuTQ1rTtXEY6w5m8dOW4xSVuf+cM1sV2k9axLLxgwgPMHD0lOs7mlX1alp3g4wcSR+gOlK6Y4ftlmjV2+p1lPwQHOz+m4GbZgIAXUw0RRvL7zAoCgV/r/L4lrqYaG9qeIYVTtqKqwoyV/+PwGitfGzdfYfQK5WPh9rNWANtlVu+VR+7K6fq+7irX9X6u+PNcQkrc5/8ts2omFIhutT9t7vLjlbchYovcjOyELjm4Er7z80KMtzG3rh/mf3n1rmpbiLhxv0V34g7ZB9xG3vDwTX2n3ud3O829srDFfvW78QOt7FDjm5x+rMzjmU5voczjnV0rLszjvvueEyccTymjsfaGcfflePv0BnHc8Dx3HDG8dxyPOec8XTOOvKla0R1rpme5Jee+ZlzIreUl37b7eQdKxvTq7HH8ufc1ZuEUH9u/3QDWQWlLNx+gl2peZgsnudULiy1MO6zDexI9e4uU0yw6wFAtSEJ0PkiLw9L3pmZvZ2HDs0AecsqLozuRpSdpo2s3R2g3Sfc1FecNU1yT1Q7tuWp6o/OaJJX/XJjC6s/jDXUVP2LY4Cl+h+Ojh9Eni58WsXq9GdnHMty92FXlTd19+aYeHOsvfkdenVueHHOiQr1fc38c7/7LyrV0btpBC9c2Zb8UjN70vIJ9NN4NQB6d1o+Fuu5sQCFNIHVEf8+F1G0dl29la+OiWHfwEFO7wL5d+9ua3pz1flSa7uxazmVU16YmpArR5D6zDNu71AZu3SuXaXrqBO1qB2zF3chzV5MOGbxoiOjlerHenNpPDcuozVTX/vpzbH25nfozbnhzTknHNTBNdNfp6aozPnvylIHv5eOkxYREWjAX6+hbaNgMvJLaBIRQEpWIZ5uAsUE6UnPN6HXqKp1xwjAw/RCtSIJUF2pr+Hv5VQAxcVOL4SFf/8NajdnicWKYrHYR7WgKOQtWeJx3qKAAQNqWl2g/m5bCu8U6KvfpyDPUP3YQm31f79FLntMnMmbAa/eXM7rZ3jCme/hej7cyixU/wJcBlT3aBehI6iaR7FQ6wfVbGLKMwRANVuuvDnnRIW6uGZq1a6v6+76Hp9ONFS4T7jzSy3kl9pG13V7ZQlqla33R3XSmfTyztnhAXrS8qp33rnbn9qSJrA6Ut+rp1vSPbQNu0uAFAXTseME9Oxpf2zNL7DdNXKj8G/3/Q48iQmRBOhckGsMqXZslr/7eUIqx4ZWvw561x34qypUGz0HlSvx4hLmmBJ4SobMLn52xrEsb5I3b+ruzTHx5lh78zv05tzw5pwTFerimml207ykcpNMmK3w6apk9F6O0LJWM/mBioQj1IvV5NWSAJ37QoYNa9D393S3xtCkMdoY2/gfbWIiodddCwHuv6UVbz1z9mhvhAdV/6It6k9OYPU/uNKCK/p9ebqonQiuGE/mKdaxXE/JR0ZgxYgPT3d4sv1C7T+Xeejgma+tON9LPdynKVIbnP7sjGNZju/hjGMdHevujOO+Ox4TZxyPqeOxdsbxd+X4O/QU66lcR96cc6JCXVwzIwP1rp8L8pxgWeuxf87pRd33pNvmy6tOrqV19+W+liQBqiNWN5OUnZX3dzMMXlU+S7S5fDZoXWwsAFoPCZA5s3Yd5kKMrv8QTzuf+3DUJ+/6fLhXoK8YfO8pocg0htp/9nTnIy2g4q6nycOl5ERAxYdniYfkI9Wh3EIPyUdaQEVikKd1/+Fx0uEORo7B/czpWcZgpz8741jWSQ93SRzr6Fh3Zxz33fGYOON4TB2PtTOOv6s0D+U6ngOO54YzjueW4znnjDfNkb50jajONdOTjALXn0Wxwa7/Rk7nIp6Wy/DE3Q2b012TruwU6zL2dM+M0//7u1lktbYkAaojZanORz2oDAYin3zS/YuNRnAzxw/gvokLKNnjfKItAH2zZgD49+gBYJ9R1q99e/fvWUvV6TNZv6unnb+8OS6e/oj1Dh0rPcX6Wyo60nu67ASYK5Jud8OCAYLMFTPy6jzGVpRr8DCqKshcMfrRz8M0A46xRov7+VACzCanPzvjWJbjezjjWEdPsY777nhMnHE8po7H2hnH31WAh3IdzwF/D8es0mg4D515vfng8aVrRF0sen5RswiXz5nMrtPJ08+EODRP1WQGZr2bXstqle35do1sTaSnq+OYCJ2eYPv0pqh67EsqCVAdCRkxnLi3pxNUZVZVpbSUzGnT3L+4uLjSMHaV0UmW7qHDcsTd97h87vSU9NYC21DashO2ZC3gInfTxtV+GLyfTvrYnwv0Dh+6nj5M/B2GZnuKDShzOGe9iNV4+E4f6BjrYQi6v0OszkOs0aGzr97qPrHyc/iw9/Pwwe9YltFDh2LHOjrW3RnHfQ/0FOtwTAM8xDr+rryJ9fcwbN8xVu/FnFeiQl1cM9cecD0NwtpqrKqeEFbRMlBQWv1pHU6zuBndFeqnw2S2Mn3xXgDCA8qXhXHyktMzR3eIq59ZoOEcSYBmzJhBUlISfn5+9O7dmw0bNriMnTVrFgMGDCAsLIywsDAGDx58RryiKEycOJFGjRphNBoZPHgw+/e7nySttnIXLCR1/BPkL1hYabs6KIjmS9yvug4QcVfFnLBO19Bx12FZpaLQzbIJ1vLkylpgW8NFE2LLvh3X13GmtsPgQ/0lAToXmDQVvwdPzQmO/V08xRY7NOd4bIZziLV4SJccR5dZVO4vUcUOsWUq9/esirUV+2ZSuT83SzR6pz8741iW43s441jHYg+j6Bz33dOIO8djWuChKdDxd1XsIdbxHPDUF8ox1vGcE9VXF9dMjZsmrMgg9+fy+kNZ3D/INqe7iprdie7pZtbm7OIyrugQy7s3dQEgq8DWK65ltK0ZOSJAb78DZVUUDFo1zaPcN1fXRoMnQN999x3jx49n0qRJbN68mc6dOzN06FBOnjzpNH7lypWMGTOGFStWsHbtWhITExkyZAjHjx+3x7z55pu8//77zJw5k/Xr1xMQEMDQoUMpqcZkgTUVccdY2u7ZTetNGwGIfta2dlDQ5YOx5OejDnPoG6BW25bIOM3Pj9xff3NbvsrNgnnqkBAC3a2706kTisVCyJUjADClpJC3ZAlK1Tkh1GrQV/yBBA0e7LZOnlRnngdfat/3Rl32ATKpK36nnmKL9A7Jh4fLX67BMalxn3zkGyrO3zIPyYfjUPxStfvRIrkOsZ4SlXxdRR2KdO4TigKH5ws8xDqW5fgezjjWMdfDlAOO++5pegLHY+p4rJ1x/F05/g6dx1acA47nhjOO55bjOecp1hNfukZUd24cdzrGhzrdrqLyoqWqKv+fNuvvQwDotSo0agiq0gymd5E1nP6d/nM42239luxK59FvtwK2TtEKkJxViFoFsSF+FJbfdWoUYqTUbKVbk/rrUN/gCdDbb7/NPffcw7hx42jXrh0zZ87E39+f2bNnO43/6quveOCBB+jSpQtt2rThk08+wWq1sqx8lmNFUXj33Xd54YUXuOaaa+jUqRNffPEFqampzJ8/v972w1pYSMnu3eQutN0BOjnVtihj7ryfOHzd9ehiHEZbWK2VJyAsKcHskMA5a+5Ssl2fVNacHLfrduUvWEDZiTTKTg+lN5mw5hdU3Gk6PT+Q1WqfTFGl13u8QyTOE96MoqgU6+H7n0PDveKhiVZxuJvhKdbqcJfEU6zj34riYbisY1mey1U7/7mW5Vaqo8djVvG81cPdrcp1qH593fZYBSqdAzU+j8TZFB7gPPlUqSo3NSlV/gdYujud9YdsnzWlZgWzFeLDjEQEVCTjJg/Z6zVd4gBbZ+ogvzO/7IQYdVzTJQ5/vcbe5+nRy1rip9OwMzWPJhH+dE4IodBkplvjULo1vkATIJPJxKZNmxjscKdBrVYzePBg1q5dW60yioqKKCsrIzzcdtstOTmZtLS0SmWGhITQu3dvl2WWlpaSl5dX6Z+3infsJPna60h7cWKl7f69e9F2z26a/jQPdXkd0enwv+giIu77zxnlqPz9Ky4e5f8bu3e3dZR2Q6V1/U3Zv3dv9AnxKGW2zLrxnDm2YfDliY+2Uaz9YqwODvY4P1B1+ek899731MThyOrFDVnlLHWd9OZ9vKm/N8fFrHZ/R8WkrbggmjzcJSl0GL1TqnOfAOcaKtrmiz3E5vhXzAtTpHd/Lmf7VwwIKPBw5yPLYb6ZfA8ju0451CHX6H7QQbbDyK9sD6PAHMtyfA9nHOuY5WGuHMd9dzwmzjge0xwPdXD8XTn+Dp1xPAcKPYzscjy3HM85Zzyds4586RpRnWumJydynbd0WBXnzWN6jYqEMNv5E+Kn49VrOwAQaNAQ5q/jwMkCFFTEBBswaNVOh647bnvkUtuivGUWhfwqq8MH+2nZ9OLljL+8NUUOd6Nu7JHI5hcv57aLmpBZYGLH8TwMWjUzb+vu1b57q0EToMzMTCwWCzExleeiiImJIS0trVplPPPMM8TFxdkTntOv86bMKVOmEBISYv+XmJjo7a4QUJ7oOP4DCL/tNgBUKhUJ77xtCy4ro2j9erJmflxRQHnSEXjpJbRav474t6fbO0OX7NxJ+C0322Kc/NNER6ONsg1nDRk1yt7cpjrd16d8hJk+IZ62e3YT0LsXAIZWLQEwHz2GJiIC1enEp6ys4udaqM76MNsjmlZ6XOrmwrg3rPLvJd9N/4Xk4NhKjzP8XF/oj1cZClz1sbtyqr6Pu/pVrb+jqvtd9bi4szmqpdvn8xw+dNfHtHUbe8Jh39fEdnAbezAs3v7z3/Hu+4ttj2hm/3llQhe3sZujK9aKX5nY1W3s2riKOi5P7OY29s+4inXlVyS4j10Z38Xpz844luX4Hs441tGx7k7r4LDvjsfEaazDMXU81s44/q4cf4fOOJ4DJzwMmXc8t/I8JKOezllHvnSN8GZNLVf6t7DVrerNPYNWTXCVOzJTr+/IG9d3JjXH1hpQbLYwoKWtNeHBS1owtH0sRr2GxDAj00d1YcnjA0mKdPKlxOG9QoxaooIMBBq0aNUqQv11dEsMBSCvxMxFry9lxd6TXN4uxp6AZBSU4qfT8MrIDmybNITuSWFc0aER0dWYt6g2zuv7lFOnTuXbb7/lp59+wq8ai3u6MmHCBHJzc+3/jh6t/qJ/1aVYLBSsXVvp1rAmIgL/Pn1sD8qbxDSBgWgCAzn1zbcohbZOy0pJCbkLFtpWmnfyz3LyJKkvTgSVity5c7Hmly+cWP7XpG9ecUE8eMVw2zIYgLFdO/t2a2EhiqJUzCd0lpq/rFX/4p0NBzgdW+Wxu9byqt+6FMX1leWMWDff2KqW4za2ymO3d46r7PcZx8UNj7EOT3uOdWh28dSk5NCXxOJFE5jVwzd6xzpaPFyiHL/xeyrX8W/PYx0cYz0051Qqy4tYT3crHPfd0+/NWqmJ0X0dHH9Xioe+W1Yvmuwq1dHDaVar8/sCvkbUhRKzBZXKdsfntECDBrPFSkKYf6XaPPvjdp78Yas9tqTMyuXv/AnAe0v38/PWVIa2j2XbsVxO5BZz2dsrOZBReMZ7nm7K0qhgzcEs7uzXlIJSM4nh/uQUlbH5aA5gS8rS8kqZ+PNOVh/IJNioRa2CNQ4j1/JLyth6NKde+/6c1qAJUGRkJBqNhvQqyzykp6cTG+s6awaYNm0aU6dOZfHixXTqVPGt6/TrvCnTYDAQHBxc6V9Nne4LVLLbdgfIdOwYJbt3c3L62+R++x3BI2wdkVX+/lhOnaLodLNc+YWz4O+/Of7cc5hSUioKVanw79rF1kHZyT//Xr2Inzat4gJltt12PN3HR+VwUTYlJ1O0bn15sQ4XQvvkC7ZtSh1M7Ohs4daqWlRZNdqguJ4/JCm/csd4d3OjJBRWnsQx0s1q2jHF2W4fuyun6vu4q1/V+juqut9Vj4s7rXKOu33e4DBMu62HFb2Dyirq3CY7xU0khJdWNBW3zz7iNjamqOKYtvMQm1BQcTFse8p9bFJ+xd+5p/o2czimrT0cB8dj6un4OpbVzMPvzbGOjnV3xnHfHY+JM47H1PFYO+P4u3L8HTrjWF/Hc8MZx3PL4GGaAU/H1JEvXSOqc8305Mt1hyvliEadmoJSCxYFDDoN0UG2L7cqbOtsqVUq/LRqNCpoHO7PDd0SAHhueFveuL4jqw5kolLBsz/+i6Io+Os19iTKMT/Wa9Q0jQrk1QW7adfI1uqQkmVLlk53pLYqtrj4UD8aBRvIKTYT5q/ng+X7WbIrnT1peYz/fhsxwQaGtHM/S3ldaNAESK/X0717d3sHZsDeobnP6TsjTrz55pu88sorLFq0iB7lk/ud1rRpU2JjYyuVmZeXx/r1692WWVdO9wVKvvY6wNYZOvna68hbsABD+/YEljfVKYWFtk7HKhXx776DLs7WccyceoK8eT9hdkjgAocMIfb552n+6y9O/8W/8zbGtm1AUeyTHgJgsWDo1JGSHTsr1dFS3sfJfHp1eECt1aJSqdCGhtrqVwcj5ozVaM8+GBJX6XGpxnXT25Eqt5IL3Iy4ORZYuVN4lp/rPhTp/uFuH7srp+r7uKtf1fo7qrrfVY+LO/tDE9w+7ziaaE9YY7ex+bqKW/J7wpq4jc32q2jm2B2e5DY23WGG5F3h7ss9HljRvLDHQ6zjMd0T7n7fDgY3sv/srjkSYH9IvNOfnXEsy/E9nHGso7vzwRZbse+Ox8QZx2Oa7mE2asfflePv0GkdHM4Bx3PDeWzFvnkavefpnHXkS9eI6lwzPSkuU6o8rrgvVlhq5s7+SYDt7pPFqmC2KFgUKxbF1nzWPj6Edo2CmbpoD8//tIPmUQH8+lB//ntLN8xWKDJZCPXX8cSQloT62/p6tYkNZMYtXflgTBcubRPNw99sASAh3J+YIANj+9reUwWYLFYahweQVWibKyou1MAdfZOYMG87V3+4miKTmTnjetVJfyhPGnyyhvHjxzN27Fh69OhBr169ePfddyksLGTcuHEA3H777cTHxzNlyhQA3njjDSZOnMjXX39NUlKSvV9PYGAggYGBqFQqHnvsMV599VVatmxJ06ZNefHFF4mLi2PkyJH1vj+n+wJVlTnzY7LnzOHklCkV3fFVKmImTUIxmSg7dswWqCigVqPS6VCsVigrI+zGUWgjIiDC9QyfZeknbYueHjpUsdFspvTf7VjiKi7Kurg49I1tFyqlzHYCGtq1w5JzCktGKdrYWKwmE9bc3Fofi+p8m/Gl29vSBHa6Dg3TBKachSYwRZrAymPdhkoTmIv61cUdoNP0GjUKCoEGLaeKbNf6rAITBofEolVsEEcyC1GrVZRZbE1nRp2Gw+V3biyKwpaUHG78eC3W8uM+5bqOzFmTzPTFFXPr7Ukr4J4vNvHWDZ14ZWQHHrq0Bb1fX0bvpuH8tOU4TSNtiV5kkJ6cojI2lA+VVwGKomb8kNaMH9K6zva9uho8ARo9ejQZGRlMnDiRtLQ0unTpwqJFi+ydmFNSUlA7XFQ++ugjTCYTN9xwQ6VyJk2axOTJkwF4+umnKSws5N577yUnJ4f+/fuzaNGiWvUTqq2Ie+/BWlhA1ief2pq7LBZQFNJfeQW/du0IHjGCvF9/BbUaQ4sWmDMy0ISFUpZ+kqzZs7FkZ3Ni0mSnZeviGtH4U9u0AU2++RrFVEbK2LG02rCeA4MHV2oMbrF8mcMrbdv1SUkkvP2jfeu+ARfX+f67csbF2s0FwLvYKs95Eet2KHOtYt18MHkTW4XHWC+GSFdKgLwa/l39WM/lVj/58C6xqn6s4kVC4U253iQqjvvuTR3qdkoCx+kAql+uV7Ee+NI1oi4YdWpMZismiy0dPFVkm2xQrVahUsEN3ROZ/Msu29tjuwtUUt6J5+v1KczbfAwFiAnyY86dvezlZhaUcsPMtSSEGfHTadFpVAxpF8uC7SdoExvEF3f1IiLA1rwWVn5n6PSQ/MLyEV+nCstoHxdM/5aRzPorGX8/jadZJuqVSlHcpNM+Ki8vj5CQEHJzc2vVH8hR7oIFnHxrGtFPPYmhRUtOTJxIyfbtRD/zDBFjb+fY44+T//si1EFBtFq7hqING0ifPp3SPbYpw6OfeYaMt992WrYuJoZmv/7Cnq7dSHjvXdRBwfYEaP+gS9BFR9N80e9nvK4s/SQHBg7E0Lo1zX6eD4BitbKnk22USNsd22u1zwUlJjpMXkJwieu2dbNaQ5HD8Fp3sVa1ptIii97EBpUWoHJxqisqVaXhyd7EBpqKULtZ9yjP4Xa4N7H+piK0dRWrD6iYUsFUgs7NMgXexBbojVjLR+T4mUvQm+smtkhvtA+T9iZWbzbh52YpiiKtAXP58GxvYrVmE/5uYku0Bvuwb69irWb8Ta77qJi0OkrKZ4D2JlZtNRNYR7Flah3FpydAtFoJNp3ZAbYmsd783fvSNWLH5MsJ9KvdgqgtnlvgdE0xjRoubhnFZ+N60fTZBRWLjgIWBfy0KnQaDb8+3B8ArUZFQljF8flh41Ge/OFf1FU6WDv6++lLSAy3vSbp2QX0ax5BcmYh13WL58MVByvFGnUaru7ciD3pBfz8YL9a7bMjbz6/G/wOkK84+dY0Iu65m5DyTtCxLzzP4dE3kTVzJv7du1GwYiVgu5uj0moJ6NsX/5UrseTkYj5+HNOhQy5v5yooKGYz2vBwcn/9laAhQwA48dLLKEVFBF5ccUfnyB3jCBo8mPBbb0EbFgpA6b595Pw0H2OnjmTP+QIsFrQxte+AdvpbxXeLXnIZsyGmDZP63G1/PGfJ6/hZnH/g/RvRjGcGPGB//L/l0whxcZHdF5rAo4Mesz/+YOV7xBSfchp7JCiG+y57yv74rVUf0cRFB9V0Yxh3DH3e/vi1Nf+jVc4xp7G5+gBuGl6x7y+u/5xOWYecxpZodFx71RT742c2fU2v9D1OYwGuGFmxvtxjW39gQOq/LmNHXvkapeXLGNy/fT6XH93oMvamKyaTW37xHrd7IVclr3EZO/by5zhZvqL5LXuWcMOBP13G/ufSJ0kp759x/f4/uXXvEpexjw58hH3l/UlGJK/l7p0LXMY+3e8+tke1AGBIyj88+O9PLmMnXnQn/8TaRj5efHwbT2z5zmXsaz1vY1X5cPGL0nfz/D//5zJ2etfRLG3SE4CumQd4eZ3zSVwBZnS6lt+a2S72bbMO8+bqmS5jP2k/gh9bXgJAs9xU3vvzfZexX7a+nK/aDgUgoSCTj5e7Xn/whxYD+bTDVQBEFucxZ8nrLmN/bdqX/3a29WcMKSviWzd/y0sSe/B295sAMFjL3P7d/x3Xidd73W5/LNcImxKzldou/BBs1HOq0MR9g5ozuG00y3af5L8rD2KxwuXtYjmZX4KCrSUxyE/L+ucuo93EP+jWJJyU7CK2Hs3huZ8qvvwqQJnZikZlS6KiAgwkhPuz9egpruoUx09bU2kfF0SIUc/yPSft/X0+GNOVJ+Zuo21sEN9vtO2/XqPik7E9uf/LTRj1GsqsClGBDTfhriRAZ4lSXEzGe++R/sqrlbZbTp3i8A2jUBkMtnl7tFqyv/g//Hv2wJKZibl84dKcn3+mefldmqpUWi1oNFhKSsj/YzH5f9jWHive+A/qgAAi7q1YKLV42zYUk4nwW29Bpdejb94M08FDnJgwoVKZp5Oo2vBzsyqwve5VF7B0c0PSu9gqz3kR6+qbnbNyvIt10wvIm9gqPMY6lO1NrNpDrMqLclVelVvxvNrqRayHctVexDqW62nfvClX7U251prVwe05SdV9cx9baX+8KNerWA986RpRnWumJ0WlZhTgo5UH+WhlxV0XvVbNqB4J3Pn5P/ZtIUYdbyyyfdHKKy6jRXQgg9vF0KV83h6ANxbtYfmek5SWf6FNyy8lwE/Lu6O78tMW22i+rAITOUVmsgtNjP9uKzHBBnQaNXqNmm3HbP1J1SrbSvWbjpyi1Gyl0GRh5d4M7h7QtNb7XFOSAJ0lAf37kff7Iozdu6Nv2ZL833+v1NFYKS2F0lJKc3PJzs0j4/33sRYVgdVK6OjRRD3ysK0jtBsqgwGlqMg+p5AlL5/IB+63JUinY7RaFHPFEFV1sPNZY4Ovvqo2u1tt6qoXBzdt4lUvDe5az1VV7pepVG4ublVj3XSdrFqO29gqj91e2qrs9xnHxQ2PsQ5Pe451SFRc3ecup6Litr7Giw88z0lCRVkaD6tGqfEisbJ6UQdr9ZOwSsmHF7FqD7PAOO67p9+bN4mV4+/K8XfotFzHc8BjHZyss1CdWA986RpRF5pGBbA3LZ8wfz25xWWYrQoaNdzcuzEalYqNhyvucp3ILeGbDbbpC3adyOOzYb1Yuiu90h0gs0Wxd46OCjJg1KnZf7KAx77bSqMQW5Onv17LBzd3JT7UyH/+bxP70vPZnZbPqO4JdE4MZcI8W3l/7c9kb3oBL13dntcW7qao1Myo7t5PPFxXJAE6S7SRkWgiIyk7kUrJjh1oo6NRGQxYsrNpsXwZKpWK/QMuRpeYiPnkSfs8POqwMBq9NLla72Hs0IHClSvtj5XiYjKmv40mLJywG64HQBMUROCAAbbnFYWSbdtApUITGmpbtNVoxJqfT/Ynn+L/7ju12ucSZw3RVTTPTa302ODi1jZAk7zKM3kHlhW5jE0oqDz3RoSbvgBV501xN49K1XKqvo+7+lWtv6Oq+131uLjT0sXtdXvZ1jJKsV2o2pxyP1dOUFkxueVLP7TxMAdPeEkBaYHRALTNPuw2NqboFIdDbcPJPc0DFF+Qye7y2X/beIhtkpfG1ujW5bHu96153gnWxdvmDPM0D1DL3OP8SXf7z+60PnWURU372N/DHcc6ujsfbLEV+x5fkOk2tvI8QM6bcU5z/F2FlxS4r4PDOeBpHiDHc8vgpu8YeD5nHfnSNaIumsAGtYomLdc2s3Jeie1Y3X5REhOuaEvP15ZSXFaR9MYG+5FRYPusGdU9kYGtoigoNVe6A+RIq1GRX2Lmivf+xqooHD1lOycOZRYy4v1VvHVDJ7699yJ6vb6MQIOGmGA//PUadBoVKpWKqzvHsfZgFi/9ugutWsWwjrFEBUkT2AWvYNVqQq4YRllaOkX//IMqIADLiRNoQkPRRUfb48qOHiXw8ssp3rQJS3Y21lOn2N2mLQH9+1O0ebPL8tts3kTjmR8BkD5tOtmffELjr74k5ZZbCejV0x7nOAqs9NAhsFoJHDSIxPLXWvLz2dezF0UbNtR6n6szp8X+0MpzrJSqtfi5mESt6pTy+Tp/Qlxc4I4GRld6nGkIJqYkx2lsmn/EGY8bFziftDCzytpJRwOjXX5AVl0ZPDk4lo4uEoWq0/vvD42n58l9TmOr2huWSNQJ19MWOM6bsiu8CXFFrifVy3foFLorIomm+a4/pLP8Ku4e7ohoSms3H2ppARGVYl31hQI4GlQxb8rO8Kb0dtMXKtlh3p1dEUkMTN3mMvaAw7m2O7wJw1Jcn+OOc/t4mjNot8McPAeqnM9V7YpIsv+c7GHOoJ3hFU0DjsfEmR0Oy0U4HmtPsY6/Q2cc65vvYS0wx7mI3M3VA56PqSNfukbUxTxAeq2akjILOo0KjUpFeLCeuy9uhl6r5o5+ScxYcQAVKu4b1Ixx/Zqy+fAp7vj8H37aepxf/z3zi9eul4fZfy4ymZm96jCJ4Ua+uvsibvt0PRc1jeCNGyomI07JKiIjv5QFj/SnfVwIP289jgoVXRJDCTRoWf3spQBc8+EqYoIbbmQ2yCgwp+pjFNjpkVVht9+ONSeH3IULbc1VgNrfH3VYGObjx1EHBtB640bM2dmkPvMshevWoQ0PI37GDFJuvc1l+c0XLkAXF0f2V1+R8c67WAsK0ESEg0ZLy6VLUOnPHFmQ/d33pE+aRPDIkZRs20bZ8ePomzShNDkZlVpNm39df5hUh6IoFJdZbE15Lqg0GlQOy264i0WtRu0wlYFXscXFrm/hq1SoHRab9Sq2pMQ2oaWravj71yhWKS1FsbhunvAmVmU02ucXUUymSk2gtYr187PPMq6Uldnnlap1rMGASqOp+1i93t4c7FWs2YxiMrmO1ensa+d5FWuxuJ1xvcaxVqvbiUy9itVq7dcORVHss8vXOtaLv3tfukYYdZpazwW0bHc6czceY+ORbDILTGjVKox6DWuevZS5G4/x1h97UangiSGtuat/UzLyS+j52jKeuLwlM1YePKM8tUqF2aKgVtuWyogI1KNWqcjMLyXUX8fssT3p6rBsxaYj2Vz/0Vo2PHcZ0cF+rDmQyc2frMdfr6GkzELL6CCevaINL8zfwbh+Sdw9wP3add6SUWDnIEVRMLZvjyYwkNwffkDfpAllJ06g0mqxFhQQOGAAOd9+i7WkBEtBIdrwcFR6PZrgYCyncjAdOOCyJVkXHY02Opq8hQtJf32K/Q/TWlKKSmXCkp/vtP/Q6Rmh8375hegnxhM4aBC5v/1G6UczUbR1d2o4/oE3WKzR/Sy2NY71Ym4pb2JVBkO1+wd4FavXO02Gax3r8KF6wcVqtZX60dVZrEaDqprnsFexanX9xKpU9RILco2oS3vS8lmfnMWAllGsOpDJM8Na88yP23n+p+3c0D2RkjJLpVmWd52wNdk1jw5i0aMXs2RXOtOX7LU/b1WU8rl7FcID9GQXmFCrVPRvGUlcqJGxn21g+ZODiHQxmstc3o8szF9H8+gwOieEcs8XGzFblbOy3pc7kgCdJdqoSPQtmlO8ZQv+F11E/uLFhI4aRcGKFQQOGIC1fOFTFEh7aTLF//5L2dFj9mQmaPDl+Hft6rTs0xfdrM/nEHLddaiAnLlzMXbpQtH69Zz67nuiHrj/jNfpG9luwesSE4m46y4Aoh991LZKfR3cGCwus9Bu4h+1LkcIIXzBrpeH4q+v3cfypiOnGNw2hjUHs7ixRyKjezbmtYW72X0in/4tIvE3aFCh4tO/DzH19932S71GrSIpMoAxvRtzuYt1uO6es5Eys5VCk5nBbaMZ2iGW5btP8tnqZJ4a2sYWVD4T9qoDmVzXLYHvNh4lLtSPjHwTHeI1DGkXy6erkgnVa+jWuGEToPN6NfjziX/XbpiSD2Ps2pXC1atRh4Sg0uvQhIVRtHkzgRcPsMfm/baAqIcfwf+ii0CjQR0UhCYwAH2TJk7/6eLjUUwmSnbuJGjgxahDbO368W9MBUWhYNkyp3Xy69QRoGIFeMCcl2dfjkMIIcT5pXuTMFbsPcnxnGJu7JHApiPZ5Beb6RAfglqtonlkIAWlZtLzSokN9mPiVe0qvT7QoCUpMuCMf3GhRg5lFpBfasaqwKRfdnHR68vJKChlhsMkh9HBtrvGm4/YOuNvOXKKm3omYrZa+Xt/Jjd+vJZGIX72dcQaktwBOkvC7xjL4TE349+vn609Pz+fU//3JahURD32GCFXXUXqU0/bhrCr1WDQU7JnD5jNRD/xhMfyzadywGKhZN8+TIdsJ2PesuVgtWLOrhixsKdLV4KvvJK4V19Bn2jriGg5eZK0114n4KLepL36GkC1mz7c8dOq2fjCZViOuxlFY9Cjiazo4Ok2VqdFEx1Ts9i0E2Bx0Q6vUaOJbVSz2JPpUOa6n4wmPr5msZkZUOq6L4k3scTGoinvJ2PJzgY3fTO8io2KQlN+nlhyc6HAzYgib2LDw9GUNzF4FVuQD7luVjcPCUETGFiD2AJwtzZeSDCawCDvY4uLIdv1aCICA9GUf5nxKtZkggzXI4+8ijUa0YTbJru0WCyQ5mbkmjex3vzd+9A1oi7mATpVWErf5pH8ui2Vy9/5C4tVwajT8MKItrYAlW0V+K2ThhBYvkr77X2SnJaVklXEr/+mcnHLKCyKglWBHk3COJBRwNLxA4kMNDBl4W4+X3OYRTvSGNYhlsTwAJ4Z1oaPVh5gYOtoTuaXsnhXOvFhRpY8PhA/nYb/W3uY95btd/qeZ5MkQGeJsWNHEj54nxPPPY9SVIQmPBxdYiIBF/Um65NPKu7ClH/4pD78iO1heDgh1ZiTpyzV1ns/8/0P7NvSJ00CqNQpUykpoXSvbVSNSqWyJVuKwqkvv+TU//0faDSojEbUhtoPTSwxW+nx6jJ+n/+k63qrNFx9zRv2xwvnP+myP4sFFVeOfKtasVZUjKh2LIxwmFl5wfwnXd4aVYDhlWKfcjmfS9XY3+Y/haaasb/8/Aw6xXXHZseZoOf/8iwGF6NiAG4fPIGMQFsfsB9/fR5/i+vOtI/3f5A9kbZRQt8tmEiwm2HEr/a4ldUJXQD46vfJhJe6TlRmdLiG31rY7nJ+/serxBTnuIz9stVgvmpnG3ny8dI3aOxmGPEvSX34qIttiof3V7zjdsj6ivguvNnzVgDe/GsGHbOTXcaui2nLS31szcKT1n7KRelnLnB82vbwpjx98YMAPP3Pl1xyfKvL2P0h8TxyyeMA3L/1R64+vNZlbEpgFP8Z/AwAt+xaxK37lrqMTTeGcsfQFwC48sDfPLjjZ5ex2YZAbrliMgD9jm3lhY1fuozN0/kzesTLALTJTOadVTNcxhZpDFx/le0LVFRBFl8sneIytlStZeTVU+2P5RphY2sCq10StCklh21Hc1CrVIQadSSGGzmUUciKvRnc0D2BiAA9/notF7+5goGtohjYKoqh7WMx6s8cgWbQqfnncDafrU4mp8j2OWLUa/jx/r6V+vyUmq3kl1QMKrhvYDOKTWYmzNuO2apgMlvP2grv3pAE6CwKuuQS0gwGYl58gfBbbuH4+PHkzvsJa2EhWbNm2WIGX0b89OkUbdjAickvYcnM5NDIkTSdOxeNmx7t2ijbNyR1UBDWkhK0ERH49+jhdLRC2JibK14XE42hWTNK9u7DmpeHsXs39E2bUbxpUz0cASGEEPUpLbeEyVe3r3RX54Nl+/nvygPc0D2Bz8b1wmyxsu5QNn/vz+DtJft4d+k+fn6oP1+uO8KMFQecltuneSRrDmZx60VNaB5VMVtRRkEpg9vGMKpHxdQGKpXKvsJ73ynLGN2zMc0qvcbkstP02STD4J2oj2Hwp+3rfRFRjz1K2Jgx9m2ZH/+P3HnzMB05Qsj11xH32mv256xFRezt1Zv4t6cT7GF5iuQbR2Ps2JHYF23fBBWrlQOXXErYLbcQ6bAchqNjjz+OUlxinwcI4PBNYzC0bl3tCRhdsVqtZBeZMB88c2ilndEPbVzFLWC3sXod2sTGNYs9nOz2lrU2qWnNYo+mgMn1cGpt8+Y1i009DsWuhyd7E0uTJmjLRyaZ09LgdIf72sbGx6MtH+Fizsx03/TjTWxMDNry5ievYnNOQZabZqKICLShoTWIzYEs13MnERGONjTM+9iCAkh3vp4UACEhaCMjvY8tKQF3TT/exAYEoI21za1jNpvhiJuJKb2J9ebv3oeuEeH+etS17H/Z5eXFPDGkNbddVDEv04wVB/hh0zFWPDnojPgik5lOkxfz4c1duahZBDlFzuvqp9Pwny830SUhhJeu6QCA1arQd+pybu/bhAcGtXD6uge/3kyJycKnd1TMR3fdf1fTplEwr1/bsRZ76pwMgz+HBV5yCZkzP0bbqBGGFi0p2b2L7M8/J/T668j65FPKTpywrwWmCQ6m4K+/wGrl1Hffc/yRR52WqfI30mbzZiLuGEvqsxPw69DBvrCptbiY0OuutcemPvMM2ugYop8YD0D4bbdz5PbbyZr9GYGDBpK3YCHFO3cS+7LrxQmrqzpNYFUXOvzp1wnVXujw24WTqr3Q4ed/vFbthQ5nLnur2gsdvrfy3WovdPjG3/+t9mKoL639pNqLoT634QvPi6Fqbd+2xm/6ttqLoT6wbV61F0O9a8ev1V4M9Zbdf1R7MdTr96+o9mKoVx5aXe3FUAcf+afai6H2P76t2ouh9kzbVe3FUDtmHKj2YqitTqVUezHUxnlp1V4MNbowu/qLoZYW8O3vk13GVloM1VzK/N+edxlbdTFUuUbY1EUT2GVtYpix/ADxoX60jA5iZ2oen65KZlSPBMA2T1BKdhG9moaz/lA2U37fjdmqcN+Xm9FpVOg0tvePDzWyZPzASmXf3b8pT8zdRseEULokhvDpqsMUmSovZzH+u63EhPjxzDDbqLA7+yUx+uN1zPrrEJe0iebXbalsP57LlOs60dAkATrLYl54gYz33yPt5ZexZGWjjY4mdPSNRD3wANFPPknRpk1kvPseGTNmoJSWom/ShPjp0yg7eZKif/5xWqYmzPaNMnj4cMzZp8j44H0sGZkY2ral8az/2b/tAZSlngBVxR+Yf7euxE97y/ae77yDPqkJiR9+gF+rVvV7IIQQQtS5l65pzyu/7uLBr7bYl70ID9AzpHxoe7BRx6Idaby7dD+lZguNw/25tXcTZv51kBu6J9iTGa3mzB5RV3WO4/ipIib9vINCkwUVEBGgZ8aKA4wf0opgPx3Hc4orTebYvUk4793UlWmL9zLl9922OaKAl3/byasjO9I0MqDej4kr0gTmRH02gdWUpaAQS5bz9YBUWi06h5EE5wqZCVpmgq5RrMwEXbtYmQm6ZrHnwDWiLmaCzi0qY/j7f9OneQS3XtSEiAA9yZmFNInwp0mE62Sj39Tl3Nm/KXf1d786+960fN5Zso8buifQMiaQ46eKeX7+DtrEBvHRrd1dvu6jlQf578oDTB/VmcRwf6Yv3sfe9Dz7yLC6Ik1gFyBNYACawIbLlGtCpVLZJvXSe5FE1ltsUD3FerF0oVexXvxpnjOx1Zwd93yM9a/mbL7exKIFY3U7gnoTC/h5MY2FN7GG6s2e7XWsXCPqzEd/HiQu1I9pozrbtyWGV54RO7OglGd++JdVBzKJCjLw5JDW1S6/dWwQM2+rSHSaRATw5JDWPP7dVswWK1rNmU14iqIwe3UyD1/agiHtbU3hb4/uTI9Xl7J4VzpXd47zdjfrhCRA54my1FQOXul6OHzz335FF9cwJ5EQQohzw9Ld6VzcMooHvtrE+kPZxAT7cVufJozpVdHh+8m520jPK+Wbey9Cp1Yz+dedZBaUMPX33UxfvNdpuUvGDyQ+1PmXhvySMgL9tE6TH4Cj2cVk5JfSr0VFd4xgPx1dEkPZfOSUJEDCPW10NM1+muf2eSGEEL4tJbuIL9cf4e7+TXlgUAv+PZbL5F92otOouaF7AocyCli5N4OfH+xH58RQAN64vhOD3/6TBy9pwajuCU7LjQlyfgcyu9DEB8sPMKZXotPnATIKbM2sUVWGvkcFGsgocN20W98kATpPqLRa9E2aeA4UQgjhsxRFoWN8CE+Xj8LqEB/CvvR8vlp/hBu6J3DgZAFatYqO8SH217SIDiTYT0t4gJ4kLzol55eUMe7zf2gRHchjg8+/gTOSAAkhhBAXCKsVNqfkkPRs5ekjjLq6Wd/xw+X7Wb7nJLtO5FFmtnJR8wg+vq27ffi8M5EBtjs/Q9/9iyKThR5JYbw6siMZBaW0a9RwA41kxUshhBDiAjG4XTRdEkLY8PxlbHj+Mr68qzcACWG2/jvNowMxWxW2H6+YYPRgRgF5Ja5HfDoyWRQuaxtDkJ8OVCo+ub2nx1FcC3ecQAVc1jaG+Q/2w6jTcusn69mScopuTRpuRXi5AyTqzelh8EIIITyri2HwDwxqwfUfrWHuxmOM6NiIWX8fRFW+HeCy6X/SMjqQWz5Zh8msEBagI9Cgxa+ad4juGdCU2z7dgEalwqhTk19aRn6pbUqJiAADGrWt/pdOX8nTQ9swtH0Ms1cfZmDrKBbvTGNo+1juH9SMUTPXEuavs89P1BAkARL1prjMQruJfzR0NYQQ4rxgmwm6dh/LnRND+fi27ry5aC/vLt2HxapwebsYru1W0bn5ZH4pCWH+HMwooKjUQkZ+abXX5tpxPI+tR3Psj3u9tsz+899PX2Ifcn8oo5D8kjL7CLCnhrTmj51pTJi3nbySMgINWga2im7QBVIlARJCCCEuEP2mLud4TuVJKBfvSufF+Tt4ZaRtDa8rOzXiNYd1uEbOWE2H+GC3kyAqisI7S/bxzT9HMWjVJIYbOZFbws6XhjmNPzx1BACvLthle4//rqZdXAifjO1Bl8RQHvxqM6Wu1lM7SyQBEvXGqNOw6+WhDV0NIYQ4Lxjr4G7ILw/1w1I+S/XD32zBVGZly9EchndsZI/p1rhyv5tujcPYdcLNosPAzD8P8dmaw/aZnJ/4fhtFJgslZRaXd3F+3ZbK56sPA/B/d/bi522p3P7pepY7WZS1IUgC1AAOXHoZZampZ2wPu3kMsRMn1rhcRVHI/OADTs2dizUvH2O3rjSaNAl9UpLb12V/9RXZn87GnJmJoU0bYl94HmOn2i9UZ58JWgghxFkRUd6UdexUEf8kZzOodTRNIvy5qFk4p1e+euW3XUyYt51Qfx3DOzZC7dDvaP6W4zz30/ZKZdr6c1qJDNTbZ3Ie0yuRF3/e6XYm509WJXNV50b8tCWVIKOO10Z2ZPmek3y/8WiDjwADSYAaRNIPc8Fh7abS/ftJufMugoY6v5VYXVmffEL2/31J3NQp6BISyHjvfVLuvodmC35DbXDevpu3cCEnp75B7OTJGDt3InvOF6TcfQ/Nf1+INiKiVvURQgjRMOZuPEZEgIHNR7K55+LmqFQqFm4/AUDXJqG2Yej5pew+kcf3G4/SPs6WjAxuF0OX8gkST0vNKebmT9bz1g0VX4z9dBo0apXLmZxNZis7judy/8BmrDqQxZoDWbSPC6Ffi0g2HMpm69Ecbr2oYee2kwSoAWjDwys9zpw1C13jxvj36mm7i/PhDHLmzcOSmYkmNJSgoUOJfeF5t2UqikL2F18Qed99BF12GQBxb0xlf7/+5C9dSsiIEU5fl/X5HEJHjSL0+usAiH1pMgV//knOj/OIvPeeWu2njAITQojqq4tRYABWq8IPm47RpXEIy/dkcEP57M6p5X2DtqbksHp/Jj2SwjiRW8K2ozm8eb0tuQk0aAk0VE4NsgptszW3jwvheE4xOUUmUnNsszsfyChgZ2ouSREBBJS/7tLpK/nPxc2wWBWigvy4s19TPli+n6TIANQqFZuPniIm2NCgI8BAEqAGp5hM5P3yK+F33IFKpSJv0R9kz5lD/NvTMbRogTkzk5I9ezyWU3bsGJaMTAL69rFv0wQFYezUieKt25wmQIrJRMnOnZUSHZVaTUCfPhRv3VrrfZNRYEIIUX11MQoMYNWBTI7nFBMb4segVlHEBNsW6B3esRGvLtiNxaowfclesgpMxIb48f6YrrSMqd4Cr28v3sePm49VvNf+TEbsX8U391xEn+a2VoNDGYUUlFbMK3TfwGYUm8xMmLedU4UmjHoNc8b1atARYCAJUIPLX7YMS34+IddeC0DZiRNoIyMJ6NMHlU6HLi6uWv1xzBmZAGiqNFtpIiMxZ2Y4f82pHLBYnLwmgtLk5BrsjRBCiIZ2casoVj1zCRe/uYKZt1as3B5XvpjpbRc1odBkYcH2E0QHGRha3q/HlahAWwKVUVDK9Bs7M/1G20rzN368lnaNgpl8dftK8YenjsBktvL6wj1kFpSiUqkYP6S17d/3W8krNtMsKrAud7lGJAFqYDk//EjggAHoYmyLmQYPG0r2F19w4PIhBA7oT8DFFxN0ySWotOffr8pPq2bjC5c1dDWEEOK84Ketu8UZ5m48RkSggUvbnLlQdqfEUIa2j+W2Pk24bPqf7E3Lp4PD2mBVJYYbiQoy2PvxgG0dMHf9ePRaNR3iQ1hzINOeYFmtCmsOZHF733NjXcvz71P1AlJ2/DiFa9eS8MH79m26Ro1o/vtCCtespXDNGtJefpnsT2fT5P++QKXTuSxLGxUJgCUrC53DyvCWzEwMbds6f01YKGg0WLKyKm23ZGahjYysxZ7ZlJit9Hh1medAIYQQ5U1gtU+CTvcBur5bAlqHNbrmbjwKQOqpYlKyipi/5Th+OjXx5XeGZqw4wIwVB5yWGRPsZ+/HkxhuZPrifWf047l51jqGto9lbN8kAO7u35Qn5m6jY0IoXRJD+HTVYYpMZkZ1d71y/NkkCVADypn3E5qIcAIHDqy0Xe3nR9CllxB06SWE3Xwzh4YPp2TfPozt27soCXQJCWiiIilcuw6/8oTHUlBA8b//EjrmJqevUen1+LVvT+HadQQNHgyAYrVSuG4dYbfcUkd7KYQQ4mw63Qfoxh4JlbYHG3V0SQzl7SX7eGvxXlrHBvHp2J6EBegBuKV3Y0Y4zBfkyKBV882GFPtMzj2Tws7ox3Mkq4jsQpP98VWd48guNPHOkn1k5JfSNi6YOXf2IiqoerNO1zeVcnpiAGGXl5dHSEgIubm5BAfXzzwFitXKgcGDCRkxgugnnrBvz5n3E1gt+HXqhNpoJGfePLI/+5wWK5ajDXO/aFzmrFlkzfrENgw+PoGM99+ndO/eSsPgj9wxjqDBgwm/1Zbg5C1cSOqzE4h96SWMnTqSPecL8hYtovnCBbW+C3R6FJi1qMhljEqjQeUwRN9dLGo1aj+/msUWF4OrU12lQm001iy2pMS2/LKravj71yhWKS1FsbgeQedNrMpotI8sUUwmFLPrRQ+9ivXzQ6W2fbtUyspQysrqJtZgQKXR1H2sXm9vSvYq1mxGMZlcx+p09ruzXsVaLCilpXUfa7WilJTUTaxWi0pv+3BUFAWluLhuYr34u/ela0RdjQLzZd58fssdoAZSuGYt5tQThFx3XaXtmuAgMmfNwjT1DRSrFb+WLUn86L8ekx+AiLvvRiku5sTESVjz8jB270birP9VmgOoLCUFy6lT9sfBw4djzj5FxgfvY8mwNZc1nvW/OmkCOz0R4u5OvV3GBAy8mMYff2x/vKf3IJcXTv+ePWnyf1/YH+8bOKzSvjjy69CBpj/MtT8+MGyk08knAfQtmtP8t9/sjw9eNwbTgYNOY3VxcbRYXtGsl3zzOEp27HAaqwkLo9XaNfbHR+56gKJ//nEaqzIaabNls/1xysMPUvjnX05jAdru2W3/+dhTT5L/h+vRdq03b7JfZFMnvkju/PkuY1uuWW2fpiFt6uuc+vobl7HNly5FnxAPQPq775A9e7bL2Ga//oKhZUsAMj6eSeaMGS5jk+Z+j7GjbZr+rP+bw8m3prmMbTxnDv69ewGQPfc70l951WVswsyPCBo0CICc337lxHPPuYyNf/cdgofZ5uXKW76U44897jK20euvE3qdbRBD/ppVHLvvfpexMS++QHj53dXC9Zs5Onasy9jop54k4q67ACjevpvDo250GRv54INEPfwQYJtX7NBVV7uMDb/zTmKefgoA07HjHCy/++uM4+Ss5uxs9vft5zI2ZORI4qZOAWyJx96LXP/dBw0dSsJ779ofyzVCNARJgBpIYP9+lT7ETgsaPNjeHOUtlUpF1COPEPXIIy5jHP8wTwu/9Rb7HSEhhBC+ZUNyNnd8tsHl87tert0kvecqaQJz4mw0gfkSt7ehNZpKd6gu5NvbXsWWllaaLbw2sY7NWlaTCarZBOYx1rFZy5vmMk+xjs1a3sReaE1gjk1K3sR606zlIRatFnU1m7W8ivXm796HrhENpaTMQlqu6/MgKTLgLNamdrz5/JYEyIlzMQEq2riRlHv/4/L5Nps3ncXaCCGEEOce6QN0AfLr0IFmP81r6GoIIYQQFwRJgM4Taj8/9E3OjcmjhBBCiPNd3U07KYQQQghxnpAESAghhBA+RxIgIYQQQvgcSYCEEEII4XMkARJCCCGEz5EESAghhBA+RxIgIYQQQvgcSYCEEEII4XMkARJCCCGEz5EESAghhPBhSc8u4I+daQ1djbNOlsIQQgghLhDvLNnHe8v2V9rWLCqA5U8MqnXZR7OL+GD5ftYczCIjv5SYYD9Gdo3noUtaoNe6vp9SUmbhtQW7+fXfVExmKxe3jOKVkR2ICjLUuk61IQmQEEIIcQFpFRPIl3f3tj/WquumsedgRgFWBV6/tiNJEQHsTc9nwrx/KTaZeX5EO5eve+W3XazYc5L/3tyNID8dE3/ZwX1fbuLH+/vWSb1qShIgIYQQ4gKiUauJDvJz+lxyZiHP/PAvW4/l0Djcn0lXuU5cqhrUOppBraPtjxtH+HMooxlfrj/iMgHKKynj+41Hee+mrvRtEQnAWzd0ZvDbf7I55RTdGod5sWd1SxIgIYQQ4gJyOLOQXq8txaBT061xGE8Pa0N8qBGrVeG+/9tEZJCe+Q/0I7+kjJd/21Wr98ovMRNq1Lt8fsexXMosCv3Kkx+AFtGBxIca2XykYRMg6QQthBBCXCC6NA5l2qjOzLmzF6+O7MjR7CJunLmWglIzqw5kcjCjgLdv7EK7uGB6N4vgqaGta/xehzMLmbPmMDf3buwyJqOgFL1GTYhRV2l7ZKCejILSGr93XZA7QEIIIcQF4pLW0aTlljD1992s3JdBsclCmcXKx38eIMzfQKNQP2KCK5rHujXx/g7M1+tTmLvpKFuP5qAoMLxjo2q97ou1h/n4z0NkFJSiValoFhXg9XvXJbkDJIQQQlwgcovKuP6jNWg1aj4f14ul4wfSJMKfrMKyOnuPk/klHM4spF1ssMfYqEADJouV7/45yqu/7ebRwS1Z8HB/UMHv29PIbMC7QJIACSGEEBeIj/48SFyoH9NGdaZLYijhAXqyC8toGR1Ii+hATuSU0GfKMt5ftp+Hv9lCj1eXArBi78lqlZ+WW8LPW1Pp3zKK50a09RjfISEEnUbFRysPcFOvRG7skYharaLIZMHfoOX7jUdrtb+1IU1got4oikJxmaWhqyGEEOcFo06DSqWqVRnfbkihT/MI7vhsA5uPnMJiVVCAqzvHEeavp2lkACnZRXy08iDXdo0nKdyffScL+P6fo4zo2IgBLaNclp2WW8JN/1tLfJiR54e3ZcvRUwBkFpTY+/ik5ZZw8yfrePvGLnRJDCXYT8f13RL49p+jXBuoZ/uxXCb9soNujUNJigxg85GcWu1vbUgCJOpNcZmFdhP/aOhqCCHEeWHXy0Px19fuYzmvpIzfd6ShVkGYv54mEUb2puezYm8GN3RP4OPbujPs3b8pMVtYdyiLSVe3Z+zsDXRrHManq5LdJkB/78/gcFYRh7OKuGjKMvv2y6b/xeGpIwAos1g5lFFIsaniy+/9g5rz7T9H+d+fh5i58hAXt4rklZEd+PTvZA5mFNZqf2tDEiAhhBDiAqFRq+gSH8K8B/rZt03+ZSdfrT/CDd0TaBYVSFSQgRt7JPLo4JYAHJ46gtmrkpm9Otlt2aN6JDKqR6L98dqDWYyZtY5tk4bYtyWG+9uTodP8dBoAvrirN91r0Om6vkgCJOqNUadh18tDG7oaQghxXjCWJwq1ER3kR8vooErbmkcH8vuOE7Uuu6bC/PVo1KozOjxnFJQSFdhwy2FIAiSEEEJcILo3CeNQZkGlbckZhcSHGittO91/p+JxDjqNmnYTFzktt2dSOHPu7FWjOum1ajrEh7DmQCZD28cCYLUqrDmQxe19m9SozLogCZCoN9IHSAghqq8u+gDd1b8p13+0hhkrDjCiYyO2Hcvhmw0pTLmuY6W4TYdPMfPPgwxpF8OqA5ks3H6CD27qQru4EKfl+jncnTqZX0JGfilHsmz9d/am5RNg0BAfaiTU3zYr9M2z1jG0fSxj+yYBcHf/pjwxdxsdE0LpkhjCp6sOU2QyM6p74hnvdbZIAiSEEEJcIDonhvLxbd15c9Fe3lu2n8QwIxOvasfIrvGV4u4e0Iztx3J5b+l+Av20vDCiLcM7xVXrPb5al1JpxfkbP14LwFs3dLL3ETqSVUR2ockec1XnOLILTbyzZB8Z+aW0jQtmzp29GnRFeJWiKEqDvfs5Ki8vj5CQEHJzcwkO9jzRk3BOhsELIUT11cUw+OroN3U5d/Zvyl39m9b7e51t3nx+yx0gUW9UKlWtb+cKIYQQ9UE+nYQQQggBwPwtx3nup+1On4sPNbJk/MCzXKP6I01gTkgTmBBCCF9UUGomM9/5+lxajYqEMP+zXCPvePP53eBrgc2YMYOkpCT8/Pzo3bs3GzZscBm7c+dOrr/+epKSklCpVLz77rtnxEyePBmVSlXpX5s2bepxD4QQQogLQ6BBS1JkgNN/53ry460GTYC+++47xo8fz6RJk9i8eTOdO3dm6NChnDzpfFG2oqIimjVrxtSpU4mNjXVZbvv27Tlx4oT936pVq+prF2qkLD2d4089zb7eF7GncxcOXXU1xdt31EnZ2V99xYFLL2NPp84k3zia4n//9fiavEWLOHjFcPZ06syhq66m4M8/66QuQgghxLmqQROgt99+m3vuuYdx48bRrl07Zs6cib+/P7Nnz3Ya37NnT9566y1uuukmDAbXQ+e0Wi2xsbH2f5GRkfW1C16z5OZyZMzNqLRaEmf9j2YLfiP6mWfQhNS+qS1v4UJOTn2DyAcfpOm8H/Fr3ZqUu+/BnJXl8jVFm7dw/IknCb3hepr+NI/AwZdx9KGHKdm3r9b1EUIIIc5VDZYAmUwmNm3axODBgysqo1YzePBg1q5dW6uy9+/fT1xcHM2aNeOWW24hJSXFbXxpaSl5eXmV/tWXrE8+QduoEXFTXsfYqRP6hAQC+/dD37ixPSb76685MHQoezp1Zl+//hx75NHqlf35HEJHjSL0+uswtGhB7EuTUfv5kfPjPJevyf6/Lwjs35+Iu+7C0Lw50Y8+il+7tpz66uta76sQQghxrmqwBCgzMxOLxUJMTEyl7TExMaSlpdW43N69e/P555+zaNEiPvroI5KTkxkwYAD5+fkuXzNlyhRCQkLs/xIT629myvzlKzB2aM+xRx9jX99+HLr2Ok59/739+eLtO0h/7XWiHn6E5r8vpPGs/+Hfo4fHchWTiZKdOwno28e+TaVWE9CnD8Vbt7p8XfHWbZVeAxDYr7/b1wghhBDnuwtuGPwVV1xh/7lTp0707t2bJk2a8P3333PXXXc5fc2ECRMYP368/XFeXl69JUFlR49y6ptvCb/jDiL/c6894VHp9IReO5KyE6mojUYCBw1CExiALj4ev3btPJZrPpUDFguaiIhK2zWREZQmu17h15yZiSaichOhJjICc2ZmjfZPCCGEOB80WAIUGRmJRqMhPT290vb09HS3HZy9FRoaSqtWrThw4IDLGIPB4LZPUV1SFAVj+/ZEj38cAL927Sjdv5+cb78l9NqRBPTthy4ujoOXX07AgP4EDhhA0ODBqI1GDyULIYQQoroarAlMr9fTvXt3li1bZt9mtVpZtmwZffr0cfNK7xQUFHDw4EEaNWpUZ2XWhjYqEn2L5pW2GZo3o+zECQA0gQE0nfcj8dOnoY2KIuP9Dzg0ciQWD/2StGGhoNFgqdLh2ZKZhdZNJ3BtZCSWrMp3ezy9RgghhDjfNegosPHjxzNr1izmzJnD7t27uf/++yksLGTcuHEA3H777UyYMMEebzKZ2Lp1K1u3bsVkMnH8+HG2bt1a6e7Ok08+yZ9//snhw4dZs2YN1157LRqNhjFjxpz1/XPGv2s3TMmHK20zHT6MLq5iETqVVktA377EPPUUzX6eT9nxVArXrXNbrkqvx699ewrXVsQpViuF69Zh7NLF5euMXTpXeg1A4Zo1bl8jhBBCnO8atA/Q6NGjycjIYOLEiaSlpdGlSxcWLVpk7xidkpKCWl2Ro6WmptK1a1f742nTpjFt2jQGDhzIypUrATh27BhjxowhKyuLqKgo+vfvz7p164iKijqr++ZK+B1jOTzmZjJnfkzwFcMo/nc7p76fS6OXXwIgf8UKyo4ew79nDzTBwRT89RdYrRiael60LuKOsaQ+OwG/Dh0wdupI9pwvsBYXE3rdtfaY1GeeQRsdQ/QTtj5P4bfdzpHbbydr9mcEDhpI3oKFFO/cSWx5feqCtajI9ZMaDWqH5ke3sWo1aj+/msUWF4OrSc9VqkpNjF7FlpSA1eq6Gv7+NYstLQWL64VkvYlVGY32BRatJhOYzXUT6+eHqvzvUzGZUOoq1mBApdF4H1tWhlJW5jpWr0el1XofazajmEyuY3U6VDqd97EWC0qp8xl3wfZFSKXXex9rtaKUlNRJLFot6tOxioJSXFw3sd783fvQNUKcXbIUhhP1vRRG/ooVZLz9DqYjR9AlJBB+x1jCbrwRgKJNm8h49z1K9u1DKS1F36QJkff9h2CHzt3uZH/5FVmzP8WSkYmhbVtin38OY+fO9ueP3HY7uvh44qZOsW/LW7SIjHffo+z4cfRJTYh+8kkCB9bdei+727R1+VzAwItp/PHH9sd7unZzeeH079mTJv/3hf3xvj59sZw65TTWr0MHmv4w1/74wKWXUZaa6jRW36I5zX/7zf744JVXYjpw0GmsLi6OFssrmm2TbxhFyQ7nk1hqwsJotXaN/fGR226n6J9/nMaqjEbabNlsf5zyn/9Q+OdfTmMB2u7Zbf/52KOPkf/HHy5jW2/eZL/Ipj47gdz5813GtlyzGm14OABpL7/Mqa+/cRnbfOlS9AnxAKS/+RbZLubvAmj26y8YWrYEIOODD8mcMcNlbNLc7zF27AhA1qefcvKtaS5jG8+ZQ0DvXoBtEtD0V151GZsw8yOCBg0CIGfeT5x47jmXsfHvvkPwsGGA7e/j+GOPu4xt9Prr9i8Z+StXcuy++13Gxrz4AuG33AJA4foNpIwd6zI2+qkniSgfuFG8fTuHR93oMjbywQeJevghAEr37+fQVVe7jA2/805inn4KANOx4xx0mIqkqrCbxxA7cSIA5uxs9vft5zI2ZORI+3XFWlTE3m7dXcYGDR1Kwnvv2h/LNULUFVkN/hwXdMklBF1yidPn/Lt3r/QH7K3wW28h/NZbXD7vrOzgYcPsF3shhBDnt7TcEqb+vpuV+zIoNllIigjgrVGd6JQQWuuyv1h7mI//PERGQSltGwXz0tXt6ZLovtwF/55g+pK9HDtVTNOIAJ69og2XtImudV1qS+4AOSGLodYtub1dg1hpAvM+VprAbLHSBFaz2HPgGlEXcovKGP7+3/RpHsGtFzUhIkBPcmYhTSL8aRIRUKuyf92WyhPfb+PVazvQNTGU2auTWfDvCZY/OYjIQOcjqTcdyebGj9fx9NDWXNY2mp+3pjLzz4P89vAAWscG1ao+znjz+S0JkBPnYgJUlprKwSuvcvl8899+rdSRWgghhO+Z+vseNh3JZu59fV3G/N/aw3y6KpnU3BKC/bT0TArno1tdN1meds2M1XROCOHlazoAYLUq9Jm6jLF9k3hgUAunr3nw680UmyzMvqOnfdvIGatpFxfM69d29HLvPJMmsAuQNjqaZj+5XtJCG93wtxOFEEI0rKW707m4ZRQPfLWJ9YeyiQn247Y+TRjTy7bc0r/Hcpj86y7evrEz3ZuEkVNUxj+Hsz2WazJb2XE8lwcGVUzjolar6Nciks1Hcly+bsuRU9w1oFmlbRe3imLxzpqv+FBXJAE6T6i0WvRNmjR0NYQQQpzDUrKL+HL9Ee7u35QHBrXg32O5TP5lJzqNmhu6J5CaU4y/TsNlbWMINGhJCIMO8SEeyz1VZMJiVc5o6ooKNHAwo9Dl6zIKSokM1Fd5jZ7MAtdNumeLJEBCCCHEBUJRFDrGh/D0sDaALbnZl57PV+uPcEP3BPq3jCI+zMjFb65gYKsoBraKYmj7WIx6TQPX/Oxr0IkQhRBCCFF3ooP8aBlduXNx8+hAUnNsndIDDVp+e7g/79/UleggA28v2ccV7/1FbrHrQQEAYf56NGrVGXduMgpKiXLRARpsd4gyCyoPDMgoMLnsNH02SQIkhBBCXCC6NwnjUGZBpW3JGYXEh1aMTNNq1PRvGcmE4W1Z9NgAjp0qZu1B9wtg67VqOsSHsOZARZzVqrDmQBbdmoS6fF3XJmGVXgOwan8G3ZqEebFX9UOawES9URSF4jLXw7OFEEJUMOo09mkoauqu/k25/qM1zFhxgBEdG7HtWA7fbEhhynW2EVfLdqeTkl1Er6bhhBh1rNibgVVRaBYV6LHsu/s35Ym52+iYEEqXxBA+XXWYIpOZUd0T7THjv9tKTIgfz5Q3wd3ZL4nRH69j1l+HuKRNNL9uS2X78VymXNepVvtZFyQBEvWmuMxCu4muZygWQghRYdfLQ/HX1+5juXNiKB/f1p03F+3lvWX7SQwzMvGqdozsapu1PdioY9GONN5dup9Ss22SxPfHdKVVjOc5ea7qHEd2oYl3luwjI7+UtnHBzLmzF1FBFc1Zx3OKKyVx3ZuE895NXZm+eC9v/bGXpEh//ndbj3qZA8hbMg+QE+fiPEDnoyKTWRIgIYSoprpIgHydzAMkzglGnYZdLw9t6GoIIcR5wajzvZFYDUkSIFFvVCqVfJsRQojzwPGcYi5/+0+Xzy8ZP7BSR+oLgXw6CSGEED4uJsjAwkcGuH3+QiMJkBBCCOHjtBo1SZG1Wyz1fCPzAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifIwmQEEIIIXyOJEBCCCGE8DmSAAkhhBDC50gCJIQQQgifU6MEqLCwsK7rIYQQQghx1tQoAYqJieHOO+9k1apVdV0fIYQQQoh6V6ME6MsvvyQ7O5tLL72UVq1aMXXqVFJTU+u6bkIIIYQQ9aJGCdDIkSOZP38+x48f57777uPrr7+mSZMmXHnllcybNw+z2VzX9RRCCCGEqDMqRVGUuijogw8+4KmnnsJkMhEZGcl9993Hs88+i7+/f10Uf1bl5eUREhJCbm4uwcHBDV0dIYQQQlSDN5/f2tq8UXp6OnPmzOHzzz/nyJEj3HDDDdx1110cO3aMN954g3Xr1rF48eLavIUQQgghRJ2rUQI0b948PvvsM/744w/atWvHAw88wK233kpoaKg9pm/fvrRt27au6imEEEIIUWdqlACNGzeOm266idWrV9OzZ0+nMXFxcTz//PO1qpwQQgghRH2oUR+goqKi87JvT3VJHyAhhBDi/FPvfYDMZjN5eXlnbFepVBgMBvR6fU2KFUIIIYQ4K2qUAIWGhqJSqVw+n5CQwB133MGkSZNQq2W1DSGEEEKcW2qUAH3++ec8//zz3HHHHfTq1QuADRs2MGfOHF544QUyMjKYNm0aBoOB5557rk4rLIQQQghRWzVKgObMmcP06dO58cYb7duuuuoqOnbsyMcff8yyZcto3Lgxr732miRAQgghhDjn1Kh9as2aNXTt2vWM7V27dmXt2rUA9O/fn5SUlNrVTgghhBCiHtQoAUpMTOTTTz89Y/unn35KYmIiAFlZWYSFhdWudkIIIYQQ9aBGTWDTpk1j1KhR/P777/Z5gDZu3MiePXv44YcfAPjnn38YPXp03dVUCCGEEKKO1HgtsMOHD/Pxxx+zd+9eAFq3bs1//vMfkpKS6rJ+DULmARJCCCHOP/U6D1BZWRnDhg1j5syZTJkypcaVFEIIIYRoKF73AdLpdPz777/1URchhBBCiLOiRp2gb731VqedoIUQQgghzgc1Xgpj9uzZLF26lO7duxMQEFDp+bfffrtOKnchyvjgQzJnzKi0Td+0Kc1/X1irchVFIfODDzg1dy7WvHyM3brSaNIk9B76ZGV/9RXZn87GnJmJoU0bYl94HmOnTrWqixBCCHGuq1ECtGPHDrp16wbAvn37Kj3nbokMYWNo2YLGs2dXbNDW6NdQSdYnn5D9f18SN3UKuoQEMt57n5S776HZgt9QGwxOX5O3cCEnp75B7OTJGDt3InvOF6TcfQ/Nf1+INiKi1nUSQghxdr2zZB/vLdtfaVuzqACWPzGoTsr/cPl+lu85ya4Teeg0arZPHurxNYqi8M6SfXzzz1HyisvokRTGqyM70jQywONr61ONPnlXrFhR1/XwLRot2qioMzYrikLmhzPImTcPS2YmmtBQgoYOJfaF590WpygK2V98QeR99xF02WUAxL0xlf39+pO/dCkhI0Y4fV3W53MIHTWK0OuvAyD2pckU/PknOT/OI/Lee2q5k0IIIRpCq5hAvry7t/2xtg7X5DRZFIZ3bES3xmF8t/FotV4z889DfLbmMNNHdSYx3J/pi/dx++z1LHl8IH46TZ3VzVu1uvVw4MABDh48yMUXX4zRaERRFLkDVA2mI0fYP+BiVAYDxi5diB7/OLq4OPL/WEz2nDnEvz0dQ4sWmDMzKdmzx2N5ZceOYcnIJKBvH/s2TVAQxk6dKN66zWkCpJhMlOzcWSnRUanVBPTpQ/HWrXWyn0IIIc4+jVpNdJCf0+eSnl3AKyM7sHRXOusOZREdbGDCFW0Z3rFRtcoef3krAOZWM/lRFIXZq5N5+NIWDGkfC8DbozvT49WlLN6VztWd46pVTn2oUVqYlZXFZZddRqtWrRg+fDgnTpwA4K677uKJJ57wqqwZM2aQlJSEn58fvXv3ZsOGDS5jd+7cyfXXX09SUhIqlYp333231mWebcbOnYib8jqJn8widtIkyo4d4/Ctt2IpKKTsxAm0kZEE9OmDLi4OY6dOhDmst+aKOSMTAE2VZitNZCTmzAznrzmVAxaLk9dEYM7MrNnOCSGEaHCHMwvp9dpSBry5nEe/3cLxnOJKz7+9eC9XdIjl90cHMLJLPA9/s4UDJ/PrpS5Hs4vJyC+lX4tI+7ZgPx1dEkPZfORUvbxnddUoAXr88cfR6XSkpKTg7+9v3z569GgWLVpU7XK+++47xo8fz6RJk9i8eTOdO3dm6NChnDx50ml8UVERzZo1Y+rUqcTGxtZJmWdb4MUXEzxsGH6tWxM4oD+J//sYa14++Yt+J3jYUKylpRy4fAgnXnyRvCVLUMzmhq6yEEKI80SXxqFMG9WZOXf24tWRHTmaXcSNM9dSUFrxWTK8YyNu6tWYZlGBPDGkNR3jQ/h8zeF6qU9GQQkAUYGV+6JGBRrIKCitl/esrholQIsXL+aNN94gISGh0vaWLVty5MiRapfz9ttvc8899zBu3DjatWvHzJkz8ff3Z7ZjB2EHPXv25K233uKmm27C4KJjr7dlNjRNcDD6pCRMR1LQNWpE898XEjtxIiqDH2kvv8yRW29DKStzW4Y2ypZZW7KyKm23ZGaijTyzrxGANiwUNBonr8lCGxnp9DVCCCHObZe0jmZEp0a0bRTMwFZRfDauF3nFZSz4N9Ue061x5XU6uzUO48DJgrNd1QZXoz5AhYWFle78nJadne0yManKZDKxadMmJkyYYN+mVqsZPHiwfUV5b9W0zNLSUkpLKzLRvLy8Gr1/TVgLCzEdPUrI1VcDoPbzI+jSSwi69BLCbr6ZQ8OHU7JvH8b27V2WoUtIQBMVSeHadfi1bQuApaCA4n//JXTMTU5fo9Lr8WvfnsK16wgaPBgAxWqlcN06wm65pW73sajI9ZMaTaVRam5j1WrUfn41iy0uBlervqhUqI3GmsWWlIDV6roaDn8nXsWWloLFUiexKqPR3jfPajKBm7uKXsX6+aEq71ypmExu71Z6FWswoNJovI8tK3P7ZUGl16MqH3HpVazZjGIyuY7V6VDpdN7HWiwopa6/Aau0WlR6vfexVitKSUmdxKLVoj4dqygoxcV1E+vN370PXSPqQ4hRR9OoAA5nuTkW9Sgq0HaMMwpKiQ6uON4ZBaW0a9SwS03VKAEaMGAAX3zxBa+88gpgG/putVp58803ueSSS6pVRmZmJhaLhZiYmErbY2Ji2FONjr91WeaUKVN46aWXavSe3kp/400CLxmELi4e88mTZH74ASq1muArR5Az7yewWvDr1Am10Ujur7+g8vNDF+e+k5hKpSL89tvJnDkTfVITdPEJZLz/PtroaHtyA3DkjnEEDR5M+K22BCfijrGkPjsBvw4dMHbqSPacL7AWFxN63bV1sq+KolBcZuFIt+4uY4wDBhA947/2xyl9+6OUOL9wGnr0IHb2Z/bHRy8bjPWU8zZkffv2NPrmW/vjYyOuxJKa6jRW17w5cT/Ntz9OvWEUZQcPOo3VxMWRsOgP++MTt9yKaedOp7HqsDAS//zL/jjt7nso3bjRaazKz0hjh75qJx96mOK//3YaC9Dk3+32nzOefIqiJUtcxiauW2+/yGa+8CKFv/ziMjZh5Z9owsMByHrtdQq++85lbPzvi9DGxwNwavo75M353GVso3k/oW/RAoCc/84kd+ZHLmNjv/4GQ4cOAOR+Noecd1zPKxbz6Wz8yhdkzv/2W7Jff91lbNSHM/C/+GIACn7+mawXX3QZGzltGgFDbMN7Cxf/QeaTT7qMjXjlFQKvGQlA0V9/kfHQgy5jw597jqCbxgBQ8s8/pN91p8vY0MfHEzJuHAClO3aQdvMYl7Eh991P6AMPAGA6cIATbv6Gg8feQVh5X03z8eMcv2KYy9jA0aOJeP4FACzZ2RwbNNBlbMDVVxP56muALfE4elFvl7H+l19O1PSK36tcI2yMOk2dDyQqLDVzJKuIa7tWJJFbjp7i+u4JlR63j6ufZCQx3EhUkIE1B7JoHxcCQH5JGVuP5nDrRU3q5T2rq0YJ0Jtvvslll13Gxo0bMZlMPP300+zcuZPs7GxWr15d13WsdxMmTGD8+PH2x3l5eSQmJtbLe5nT00h94kksOTlowsPx796NpO++RRsejiY4iMxZszBNfQPFasWvZUsSP/ov2rAwj+VG3H03SnExJyZOwpqXh7F7NxJn/a/SN6eylBQsDheD4OHDMWefIuOD97FkZGJo25bGs/5XZ01gxWUW2k38g9/dxPy5L4NJEysuFj+VmXE+dgH+Sc7mGYfYbwtNhLiI3XE8l8scYj8/VUyMi9gDJwsY7BA782QBrv4sU08VM8Qh9r3jubRyEXuq0MRQh9g3krNxNcVkcZmZdg6xL+3LoJeLWKBS7HM70xngJrb7q0sp1drOg/FbUrncTWz/N1aQawgE4IFtR7nKTezgt//iZIAtWbprRzI3uIm9+sPVpATbPjBu2X2AW93Ejv54LfvCjgNw/f693O0mduzsDWz/NRuAKw/txnXqAfd/uYl/lto+OAcf2YG74RqPf7eNVatsP/c/vg13E1E8N28HSzfZvvH3TNvFy25iX/ltN7/tsv3uOmYc4E03sdMW7+XHZFtsq1MpvOcmdsaKA3yVZottnJfGx25iZ69O5tNcW2x0YTZz3MR+s+Eo/y0/10JKC/jWTez8Lam8XR5rMJcy303sHzvTed3hHJZrhM2ul4fir6/dvHCvLdjFZW1jiA81cjK/hHeW7EejVlUabbVwexqd4kPpkRTG/K2pbDuaw5vXV28C3OM5xeQUmUjNKcFqVdiZmgtAUkQAAQZb3S+dvpKnh7ZhWIdYVCoVd/ZrygfL95MUGUBiuJHpi/cRE2xgSDtXR/zsqPFq8Lm5uXz44Yds27aNgoICunXrxoMPPkijRtUbSmcymfD39+eHH35g5MiR9u1jx44lJyeHn3/+2e3rk5KSeOyxx3jsscfqrMzTZDX4ulFksn2oG8yub91bVWrKNDr7Y3exikqNqcaxJsDVqa6iVKuvUazeUoZKcX17+3Ti4W2szlKGuq5iNXoo/1aps5hRK66by7yJNWl0KCpbs5bWakZjrZvYMo0Oaw1iNVYLWqvr5rIytRarWuN1rNpqQecm1qzWYqlJrGJFZ3HdDGdRazCrtV7HqhQr+jqKtao0lGnKP5AVBYPFdfOed7HV/7v3pWtEXSRAD329mQ3J2eQUlREeoKdHUhhPDW1NkwjbpINJzy7glWvas3hXOuuTs4kOMvDsFW24slP1hqM/8f02ftx87Izt39xzEX2aR9jf460bOjGqh+1GwumJEL/ecJS8kjJ6JoXxyjUdaBYVWKt9dcabz+8aJ0B1oXfv3vTq1YsPPvgAAKvVSuPGjXnooYd49tln3b7WWQJU2zJPkwSobpxuAhNCCOFZfTSBVZX07AI+vq07Q9s7H0l9vvPm87vGqWZOTg4bNmzg5MmTWKt08Lr99turVcb48eMZO3YsPXr0oFevXrz77rsUFhYyrrzd+/bbbyc+Pp4pU6YAtjs8u3btsv98/Phxtm7dSmBgIC3K+xd4KvN8VLRxIyn3/sfl8202bzqLtak+lUpV628zQgghRH2o0afTr7/+yi233EJBQQHBwcGVMlaVSlXtBGj06NFkZGQwceJE0tLS6NKlC4sWLbJ3Yk5JSUHtMIV3amoqXbt2tT+eNm0a06ZNY+DAgaxcubJaZZ6P/Dp0oNlP8xq6GkIIIS5gM1YcYMaKA06f65kUzpw73fVMPP/UqAns9AzQr7/+utPh8Oc7aQITQgjha3KKTOQUOe8T5qfTEBviqvv5uaPem8COHz/OI488ckEmP0IIIYQvCvXXE+qv9xx4gajRTNBDhw5lo4v5TIQQQgghznU1ugM0YsQInnrqKXbt2kXHjh3R6XSVnr+6fFZjIYQQQohzUY36ADl2TD6jQJUKi5tp+c8H0gdICCGEOP/Uex+gqsPehRBCCCHOJ171ARo+fDi5ubn2x1OnTiUnJ8f+OCsri3bt2tVZ5YQQQggh6oNXCdAff/xRadX0119/nezsbPtjs9nM3r176652QgghhBD1wKsEqGp3oQZcRUMIIYQQosZqNAxeCCGEEOJ85lUCpFKpzliorb4XbhNCCCGEqGtejQJTFIU77rgDg8EAQElJCffddx8BAQEAlfoHCSGEEEKcq7xKgMaOHVvp8a233npGTHUXQhVCCCGEaCheJUCfffZZfdVDCCGEEOKskU7QQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA550QCNGPGDJKSkvDz86N3795s2LDBbfzcuXNp06YNfn5+dOzYkYULF1Z6/o477kClUlX6N2zYsPrcBSGEEEKcRxo8Afruu+8YP348kyZNYvPmzXTu3JmhQ4dy8uRJp/Fr1qxhzJgx3HXXXWzZsoWRI0cycuRIduzYUSlu2LBhnDhxwv7vm2++ORu7I4QQQojzgEpRFKUhK9C7d2969uzJhx9+CIDVaiUxMZGHH36YZ5999oz40aNHU1hYyG+//WbfdtFFF9GlSxdmzpwJ2O4A5eTkMH/+/BrVKS8vj5CQEHJzcwkODq5RGUIIIYQ4u7z5/G7QO0Amk4lNmzYxePBg+za1Ws3gwYNZu3at09esXbu2UjzA0KFDz4hfuXIl0dHRtG7dmvvvv5+srCyX9SgtLSUvL6/SPyGEEEJcuBo0AcrMzMRisRATE1Npe0xMDGlpaU5fk5aW5jF+2LBhfPHFFyxbtow33niDP//8kyuuuAKLxeK0zClTphASEmL/l5iYWMs9E0IIIcS5TNvQFagPN910k/3njh070qlTJ5o3b87KlSu57LLLzoifMGEC48ePtz/Oy8uTJEgIIYS4gDXoHaDIyEg0Gg3p6emVtqenpxMbG+v0NbGxsV7FAzRr1ozIyEgOHDjg9HmDwUBwcHClf0IIIYS4cDVoAqTX6+nevTvLli2zb7NarSxbtow+ffo4fU2fPn0qxQMsWbLEZTzAsWPHyMrKolGjRnVTcSGEEEKc1xp8GPz48eOZNWsWc+bMYffu3dx///0UFhYybtw4AG6//XYmTJhgj3/00UdZtGgR06dPZ8+ePUyePJmNGzfy0EMPAVBQUMBTTz3FunXrOHz4MMuWLeOaa66hRYsWDB06tEH2UQghhBDnlgbvAzR69GgyMjKYOHEiaWlpdOnShUWLFtk7OqekpKBWV+Rpffv25euvv+aFF17gueeeo2XLlsyfP58OHToAoNFo+Pfff5kzZw45OTnExcUxZMgQXnnlFQwGQ4PsoxBCCCHOLQ0+D9C5SOYBEkIIIc4/5808QEIIIYQQDUESICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM+RBEgIIYQQPkcSICGEEEL4HEmAhBBCCOFzJAESQgghhM/RNnQFhBBCCFE33lmyj/eW7a+0rVlUAMufGFSrchVF4Z0l+/jmn6PkFZfRIymMV0d2pGlkgNvXfbH2MB//eYiMglLaNgrmpavb0yUxtFZ1qSuSAAkhhBAXkFYxgXx5d2/7Y6269o09M/88xGdrDjN9VGcSw/2Zvngft89ez5LHB+Kn0zh9za/bUnn1t928em0HuiaGMnt1Mrd/up7lTw4iMtBQ6zrVljSBCSGEEBcQjVpNdJCf/V94gB6ouIvTd8oyWj3/O71eW8rkX3Z6LE9RFGavTubhS1swpH0sbRsF8/bozqTnlbJ4V7rL132yKpmbeiVyY49EWsYE8drIjhj1Gr7feLTO9rU25A6QEEIIcQE5nFlIr9eWYtCp6dY4jKeHtSE+1MjvO9KYvSqZ92/uSquYIDLyS9l9Is9jeUezi8nIL6Vfi0j7tmA/HV0SQ9l85BRXd4474zUms5Udx3N5YFBz+za1WkW/FpFsPpJTJ/tZW5IACSGEEBeILo1DmTaqM82iAjiZX8p7S/dx48y1/PH4xaTmFBMVZKB/i0h0GjXxocZq9cfJKCgBIKpKs1VUoIGMglKnrzlVZMJiVc5o6ooKNHAwo7BmO1fHzokmsBkzZpCUlISfnx+9e/dmw4YNbuPnzp1LmzZt8PPzo2PHjixcuLDS84qiMHHiRBo1aoTRaGTw4MHs37/fRWlCCCHEheGS1tGM6NSIto2CGdgqis/G9SKvuIwF/6YyvGMjSsosXPzmCp798V8W7UjDbLE2dJUbTIMnQN999x3jx49n0qRJbN68mc6dOzN06FBOnjzpNH7NmjWMGTOGu+66iy1btjBy5EhGjhzJjh077DFvvvkm77//PjNnzmT9+vUEBAQwdOhQSkpKztZuCSGEEA0uxKijaVQAh7OKiAs1svzJQbxyTQf8dBpe/HkHN368ljIPSVBUoB/AGXd7MgpKz7grdFqYvx6NWkWmF68521SKoigNWYHevXvTs2dPPvzwQwCsViuJiYk8/PDDPPvss2fEjx49msLCQn777Tf7tosuuoguXbowc+ZMFEUhLi6OJ554gieffBKA3NxcYmJi+Pzzz7nppps81ikvL4+QkBByc3MJDg6uoz31PYqiUFxmaehqCCHEecGo06BSqeq0zMJSM32nLuexwS0Z169ppecOZhRw2fQ/+e3h/nSID3FZhqIo9Hp9GfcOaMY9FzcDIL+kjO6vLmXaqM5O+wABXDNjNV0SQnjpmg4AWK0Kfacu5/a+TXhgUIs62sPKvPn8btA+QCaTiU2bNjFhwgT7NrVazeDBg1m7dq3T16xdu5bx48dX2jZ06FDmz58PQHJyMmlpaQwePNj+fEhICL1792bt2rVOE6DS0lJKSyuy1Lw8z53ChGfFZRbaTfyjoashhBDnhV0vD8VfX7uP5dcW7OKytjHEhxo5mV/CO0v2o1GruLpzHHM3HsWqKHRJDMOo0zB/y3H8dLa+QO6oVCru7NeUD5bvJykygMRwI9MX7yMm2MCQdjH2uJtnrWNo+1jG9k0C4O7+TXli7jY6JoTSJTGET1cdpshkZlT3xFrtY11p0AQoMzMTi8VCTExMpe0xMTHs2bPH6WvS0tKcxqelpdmfP73NVUxVU6ZM4aWXXqrRPgghhBDnihO5JTzyzRZyisoID9DTIymMnx7oS0SggWCjjo9WHuTV33ZjURRaxwbx6diehJUPk3fnvoHNKDaZmTBvO3klZfRMCmPOuF6V5gA6klVEdqHJ/viqznFkF5p4Z8k+MvJLaRsXzJw7exEVdG40gckoMGDChAmV7irl5eWRmHhuZKjnM6NOw66XhzZ0NYQQ4rxgdDGhoDc+vLmby+eGto9laPvYGpWrUqkYP6Q144e0dhmz+tlLz9g2tm+S/Y7QuaZBE6DIyEg0Gg3p6ZUnUkpPTyc21vkvKTY21m386f/T09Np1KhRpZguXbo4LdNgMGAwnBsZ6YVEpVLV+nauEEIIUR8a9NNJr9fTvXt3li1bxsiRIwFbJ+hly5bx0EMPOX1Nnz59WLZsGY899ph925IlS+jTpw8ATZs2JTY2lmXLltkTnry8PNavX8/9999fn7sjhBBCnHc2JGdzx2eup5/Z9fKws1ibs6fBv56PHz+esWPH0qNHD3r16sW7775LYWEh48aNA+D2228nPj6eKVOmAPDoo48ycOBApk+fzogRI/j222/ZuHEj//vf/wDbXYfHHnuMV199lZYtW9K0aVNefPFF4uLi7EmWEEIIIWw6JYSw8JEBDV2Ns67BE6DRo0eTkZHBxIkTSUtLo0uXLixatMjeiTklJQW1w0Juffv25euvv+aFF17gueeeo2XLlsyfP58OHTrYY55++mkKCwu59957ycnJoX///ixatAg/P7+zvn9CCCHEucxPpyHJw6ruF6IGnwfoXCTzAAkhhBDnH28+vxt8JmghhBBCiLNNEiAhhBBC+BxJgIQQQgjhcyQBEkIIIYTPkQRICCGEED5HEiAhhBBC+BxJgIQQQgjhcyQBEkIIIYTPkQRICCGEED6nwZfCOBednhw7Ly+vgWsihBBCiOo6/bldnUUuJAFyIj8/H4DExMQGrokQQgghvJWfn09ISIjbGFkLzAmr1UpqaipBQUGoVKoalZGXl0diYiJHjx716fXE5DjYyHGwkeNgI8fBRo6DjRwHm7o4DoqikJ+fT1xcXKWF1J2RO0BOqNVqEhIS6qSs4OBgnz6hT5PjYCPHwUaOg40cBxs5DjZyHGxqexw83fk5TTpBCyGEEMLnSAIkhBBCCJ8jCVA9MRgMTJo0CYPB0NBVaVByHGzkONjIcbCR42Ajx8FGjoPN2T4O0glaCCGEED5H7gAJIYQQwudIAiSEEEIInyMJkBBCCCF8jiRAQgghhPA5kgDVkezsbG655RaCg4MJDQ3lrrvuoqCgwG38ww8/TOvWrTEajTRu3JhHHnmE3Nzcs1jrujFjxgySkpLw8/Ojd+/ebNiwwW383LlzadOmDX5+fnTs2JGFCxeepZrWL2+Ow6xZsxgwYABhYWGEhYUxePBgj8ftfOHt+XDat99+i0qlYuTIkfVbwbPE2+OQk5PDgw8+SKNGjTAYDLRq1eqC+Nvw9ji8++679utiYmIijz/+OCUlJWeptvXjr7/+4qqrriIuLg6VSsX8+fM9vmblypV069YNg8FAixYt+Pzzz+u9nvXN2+Mwb948Lr/8cqKioggODqZPnz788ccfdVchRdSJYcOGKZ07d1bWrVun/P3330qLFi2UMWPGuIzfvn27ct111ym//PKLcuDAAWXZsmVKy5Ytleuvv/4s1rr2vv32W0Wv1yuzZ89Wdu7cqdxzzz1KaGiokp6e7jR+9erVikajUd58801l165dygsvvKDodDpl+/btZ7nmdcvb43DzzTcrM2bMULZs2aLs3r1bueOOO5SQkBDl2LFjZ7nmdcvb43BacnKyEh8frwwYMEC55pprzk5l65G3x6G0tFTp0aOHMnz4cGXVqlVKcnKysnLlSmXr1q1nueZ1y9vj8NVXXykGg0H56quvlOTkZOWPP/5QGjVqpDz++ONnueZ1a+HChcrzzz+vzJs3TwGUn376yW38oUOHFH9/f2X8+PHKrl27lA8++EDRaDTKokWLzk6F64m3x+HRRx9V3njjDWXDhg3Kvn37lAkTJig6nU7ZvHlzndRHEqA6sGvXLgVQ/vnnH/u233//XVGpVMrx48erXc7333+v6PV6paysrD6qWS969eqlPPjgg/bHFotFiYuLU6ZMmeI0/sYbb1RGjBhRaVvv3r2V//znP/Vaz/rm7XGoymw2K0FBQcqcOXPqq4pnRU2Og9lsVvr27at88sknytixYy+IBMjb4/DRRx8pzZo1U0wm09mq4lnh7XF48MEHlUsvvbTStvHjxyv9+vWr13qeTdX54H/66aeV9u3bV9o2evRoZejQofVYs7OrOsfBmXbt2ikvvfRSndRBmsDqwNq1awkNDaVHjx72bYMHD0atVrN+/fpql5Obm0twcDBa7fmxRJvJZGLTpk0MHjzYvk2tVjN48GDWrl3r9DVr166tFA8wdOhQl/Hng5och6qKioooKysjPDy8vqpZ72p6HF5++WWio6O56667zkY1611NjsMvv/xCnz59ePDBB4mJiaFDhw68/vrrWCyWs1XtOleT49C3b182bdpkbyY7dOgQCxcuZPjw4WelzueKC/E6WResViv5+fl1dp08Pz5pz3FpaWlER0dX2qbVagkPDyctLa1aZWRmZvLKK69w77331kcV60VmZiYWi4WYmJhK22NiYtizZ4/T16SlpTmNr+5xOhfV5DhU9cwzzxAXF3fGRe98UpPjsGrVKj799FO2bt16Fmp4dtTkOBw6dIjly5dzyy23sHDhQg4cOMADDzxAWVkZkyZNOhvVrnM1OQ4333wzmZmZ9O/fH0VRMJvN3HfffTz33HNno8rnDFfXyby8PIqLizEajQ1Us4Y1bdo0CgoKuPHGG+ukPLkD5Mazzz6LSqVy+6+6H3Du5OXlMWLECNq1a8fkyZNrX3FxXpk6dSrffvstP/30E35+fg1dnbMmPz+f2267jVmzZhEZGdnQ1WlQVquV6Oho/ve//9G9e3dGjx7N888/z8yZMxu6amfVypUref311/nvf//L5s2bmTdvHgsWLOCVV15p6KqJBvb111/z0ksv8f33359xw6Gm5A6QG0888QR33HGH25hmzZoRGxvLyZMnK203m81kZ2cTGxvr9vX5+fkMGzaMoKAgfvrpJ3Q6XW2rfdZERkai0WhIT0+vtD09Pd3lfsfGxnoVfz6oyXE4bdq0aUydOpWlS5fSqVOn+qxmvfP2OBw8eJDDhw9z1VVX2bdZrVbAdgd17969NG/evH4rXQ9qcj40atQInU6HRqOxb2vbti1paWmYTCb0en291rk+1OQ4vPjii9x2223cfffdAHTs2JHCwkLuvfdenn/+edRq3/jO7uo6GRwc7JN3f7799lvuvvtu5s6dW6d3yX3jbKqhqKgo2rRp4/afXq+nT58+5OTksGnTJvtrly9fjtVqpXfv3i7Lz8vLY8iQIej1en755Zfz7tu/Xq+ne/fuLFu2zL7NarWybNky+vTp4/Q1ffr0qRQPsGTJEpfx54OaHAeAN998k1deeYVFixZV6j92vvL2OLRp04bt27ezdetW+7+rr76aSy65hK1bt5KYmHg2q19nanI+9OvXjwMHDtgTQIB9+/bRqFGj8zL5gZodh6KiojOSnNNJoeJDy1ZeiNfJmvrmm28YN24c33zzDSNGjKjbwuukK7VQhg0bpnTt2lVZv369smrVKqVly5aVhsEfO3ZMad26tbJ+/XpFURQlNzdX6d27t9KxY0flwIEDyokTJ+z/zGZzQ+2G17799lvFYDAon3/+ubJr1y7l3nvvVUJDQ5W0tDRFURTltttuU5599ll7/OrVqxWtVqtMmzZN2b17tzJp0qQLZhi8N8dh6tSpil6vV3744YdKv/v8/PyG2oU64e1xqOpCGQXm7XFISUlRgoKClIceekjZu3ev8ttvvynR0dHKq6++2lC7UCe8PQ6TJk1SgoKClG+++UY5dOiQsnjxYqV58+bKjTfe2FC7UCfy8/OVLVu2KFu2bFEA5e2331a2bNmiHDlyRFEURXn22WeV2267zR5/ehj8U089pezevVuZMWPGBTEM3tvj8NVXXylarVaZMWNGpetkTk5OndRHEqA6kpWVpYwZM0YJDAxUgoODlXHjxlX6MEtOTlYAZcWKFYqiKMqKFSsUwOm/5OTkhtmJGvrggw+Uxo0bK3q9XunVq5eybt06+3MDBw5Uxo4dWyn++++/V1q1aqXo9Xqlffv2yoIFC85yjeuHN8ehSZMmTn/3kyZNOvsVr2Peng+OLpQESFG8Pw5r1qxRevfurRgMBqVZs2bKa6+9dl59GXLFm+NQVlamTJ48WWnevLni5+enJCYmKg888IBy6tSps1/xOuTqen9638eOHasMHDjwjNd06dJF0ev1SrNmzZTPPvvsrNe7rnl7HAYOHOg2vrZUiuJD9xWFEEIIIZA+QEIIIYTwQZIACSGEEMLnSAIkhBBCCJ8jCZAQQgghfI4kQEIIIYTwOZIACSGEEMLnSAIkhBBCCJ8jCZAQQgghfI4kQEIIIYTwOZIACSGEEMLnSAIkhBBCCJ8jCZAQQgghfM7/A2dbOyxrpvAAAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots()\n", "\n", "ax.set_ylabel(\"Energy\")\n", "if False: # Change to True if you want to scale energies by n^2\n", " energies_with_nist *= np.array(n_list) ** 2\n", " energies_without_nist *= np.array(n_list) ** 2\n", " ax.set_ylabel(\"Energy * n^2\")\n", "\n", "for i, label in enumerate(labels):\n", " if energies_with_nist[i] < 0 or energies_without_nist[i] < 0:\n", " print(\n", " f\"Skipping negative energy for {label}: \"\n", " f\"{energies_with_nist[i]} (with NIST), \"\n", " f\"{energies_without_nist[i]} (without NIST)\"\n", " )\n", " continue\n", " ax.plot([0, 1], [energies_with_nist[i]] * 2, \"C0\", zorder=-10)\n", " ax.plot([0, 1], [energies_without_nist[i]] * 2, \"C3--\", zorder=10)\n", " ax.text(1.05, energies_with_nist[i], label, color=\"C0\")\n", " ax.text(-0.05, energies_without_nist[i], label, color=\"C3\", ha=\"right\")\n", "\n", "ax.set_xlim(-0.25, 1.25)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAdW9JREFUeJzt3XlYVOX7P/D3GZYZQFZBFgFxN9JAUcklt0iCcs8tzR1b8JtGWS6lZaW2uJWUWZpaVlYfs9JyQw33NdFyX3AHdwhU1vv3hxfzE5kDMzDDAL5f18VVc55zznM/Z44z99xnU0REQERERESVnsbaARARERGReTCxIyIiIqoimNgRERERVRFM7IiIiIiqCCZ2RERERFUEEzsiIiKiKoKJHREREVEVwcSOiIiIqIqwtXYAFV1+fj4uXrwIZ2dnKIpi7XCIiIjoASMi+O+//+Dn5weNpviaHBO7Ely8eBEBAQHWDoOIiIgecOfOnYO/v3+x8zCxK4GzszOAuxvTxcXFytEQERHRgyY9PR0BAQH6nKQ4TOxUxMfHIz4+Hnl5eQAAFxcXJnZERERkNcacEqaIiJRDLJVWeno6XF1dkZaWxsSOiIiIyp0puQiviiUiIiKqIpjYEREREVURTOyIiIiIqggmdkRERERVBBM7IiIioiqCiR0RERFRFcHEjoiIiKiKYGJHREREVEU8EIldjx494O7ujmeeecbaoRARERFZzAOR2I0ePRpLliyxdhhEREREFvVAJHYdOnQw6sG5RERERJWZ1RO7xMREdOnSBX5+flAUBStWrCgyT3x8PIKCgqDT6RAeHo5du3aVf6BEREREFZzVE7vMzEyEhIQgPj7eYPuyZcsQFxeHyZMnY9++fQgJCUFkZCQuX76snyc0NBSNGzcu8nfx4sXyGgYRERGR1dlaO4CoqChERUWpts+cORMxMTEYOnQoAGDevHlYtWoVFi5ciHHjxgEA9u/fb7Z4srKykJWVpX+dnp5utnWXVtC4VaptydOfKsdIiIiIqCKzemJXnOzsbOzduxfjx4/XT9NoNIiIiMD27dst0ue0adPwzjvvWGTdxVFL3pi4ERERkbGsfii2OFevXkVeXh68vb0LTff29kZKSorR64mIiEDv3r3xxx9/wN/fv9ikcPz48UhLS9P/nTt3rtTxExEREZWnCl2xM5f169cbPa9Wq4VWq7VgNERERESWUaErdp6enrCxsUFqamqh6ampqfDx8bFo3/Hx8QgODkaLFi0s2g8RERGRuVToxM7e3h5hYWFISEjQT8vPz0dCQgJatWpl0b5jY2Nx6NAh7N6926L9EBEREZmL1Q/FZmRk4MSJE/rXp0+fxv79++Hh4YHAwEDExcVh8ODBaN68OVq2bInZs2cjMzNTf5UsEREREd1l9cRuz5496Nixo/51XFwcAGDw4MFYtGgR+vbtiytXrmDSpElISUlBaGgoVq9eXeSCCnOLj49HfHw88vLyLNoPERERkbkoIiLWDqIiS09Ph6urK9LS0uDi4mKxfoq73QnvY0dERPTgMiUXqdDn2BERERGR8ZjYqeBVsURERFTZMLFTwatiiYiIqLJhYkdERERURTCxIyIiIqoimNip4Dl2REREVNkwsVPBc+yIiIiosmFiR0RERFRFMLEjIiIiqiKs/kixiqoyPVKsuKdWEBER0YODFTsVPMeOiIiIKhsmdkRERERVBBM7IiIioiqCiR0RERFRFcHETgVvUExERESVDRM7Fbx4goiIiCobJnZEREREVQQTOyIiIqIqgokdERERURXBJ09QifhkCyIiosqBFTsiIiKiKoKJnQre7oSIiIgqGyZ2Kni7EyIiIqpsmNgRERERVRFM7IiIiIiqCCZ2RERERFUEb3dCZaJ2KxSAt0MhIiIqb6zYEREREVURTOyIiIiIqggmdkRERERVBBM7IiIioiqCiZ0KPnmCiIiIKhsmdir45AkiIiKqbJjYEREREVURTOyIiIiIqggmdkRERERVBBM7IiIioiqCiR0RERFRFcHEjoiIiKiKYGJHREREVEUwsSMiIiKqIpjYEREREVURVT6xO3fuHDp06IDg4GA88sgj+Omnn6wdEhEREZFF2Fo7AEuztbXF7NmzERoaipSUFISFhSE6OhpOTk7WDo2IiIjIrKp8Yufr6wtfX18AgI+PDzw9PXH9+nUmdpVE0LhVBqcnT3/KqHYiIqIHidUPxSYmJqJLly7w8/ODoihYsWJFkXni4+MRFBQEnU6H8PBw7Nq1q1R97d27F3l5eQgICChj1EREREQVj9UTu8zMTISEhCA+Pt5g+7JlyxAXF4fJkydj3759CAkJQWRkJC5fvqyfJzQ0FI0bNy7yd/HiRf08169fx6BBgzB//nyLj4mIiIjIGqx+KDYqKgpRUVGq7TNnzkRMTAyGDh0KAJg3bx5WrVqFhQsXYty4cQCA/fv3F9tHVlYWunfvjnHjxqF169YlzpuVlaV/nZ6ebuRIiIiIiKzL6hW74mRnZ2Pv3r2IiIjQT9NoNIiIiMD27duNWoeIYMiQIejUqROee+65EuefNm0aXF1d9X88bEtERESVRYVO7K5evYq8vDx4e3sXmu7t7Y2UlBSj1rF161YsW7YMK1asQGhoKEJDQ3Hw4EHV+cePH4+0tDT937lz58o0BiIiIqLyYvVDsZbWtm1b5OfnGz2/VquFVqtFfHw84uPjkZeXZ8HoiIiIiMynQlfsPD09YWNjg9TU1ELTU1NT4ePjY9G+Y2NjcejQIezevdui/RARERGZS4VO7Ozt7REWFoaEhAT9tPz8fCQkJKBVq1ZWjIyIiIio4rH6odiMjAycOHFC//r06dPYv38/PDw8EBgYiLi4OAwePBjNmzdHy5YtMXv2bGRmZuqvkrUUHoolIiKiysbqid2ePXvQsWNH/eu4uDgAwODBg7Fo0SL07dsXV65cwaRJk5CSkoLQ0FCsXr26yAUV5hYbG4vY2Fikp6fD1dXVon0RERERmYPVE7sOHTpARIqdZ9SoURg1alQ5RURERERUOVXoc+ysKT4+HsHBwWjRooW1QyEiIiIyChM7FbwqloiIiCobJnZEREREVYTVz7Ej6wsat8rg9OTpT5VzJERERFQWrNip4Dl2REREVNkwsVPBc+yIiIiosmFiR0RERFRFMLEjIiIiqiKY2BERERFVEbwqVgWfFUsl4dXERERU0bBip4IXTxAREVFlw8SOiIiIqIpgYkdERERURTCxIyIiIqoimNip4JMniIiIqLJhYqeCF08QERFRZcPEjoiIiKiKMOk+docPH8YPP/yAzZs348yZM7h16xa8vLzQtGlTREZGolevXtBqtZaKlYiIiIiKYVTFbt++fYiIiEDTpk2xZcsWhIeHY8yYMXj33XcxcOBAiAgmTpwIPz8/fPDBB8jKyrJ03ERERER0H6Mqdr169cLYsWPx888/w83NTXW+7du3Y86cOZgxYwYmTJhgrhiJiIiIyAhGJXbHjh2DnZ1difO1atUKrVq1Qk5OTpkDIyIiIiLTGHUoVi2pu3PnjknzVya83QkRERFVNiZfFZufn493330XNWvWRLVq1XDq1CkAwFtvvYUFCxaYPUBr4e1OiIiIqLIxObF77733sGjRInz44Yewt7fXT2/cuDG++uorswZHRERERMYzObFbsmQJ5s+fjwEDBsDGxkY/PSQkBEeOHDFrcERERERkPJMTuwsXLqBevXpFpufn5/OiCSIiIiIrMjmxCw4OxubNm4tM//nnn9G0aVOzBEVEREREpjPpyRMAMGnSJAwePBgXLlxAfn4+li9fjqNHj2LJkiVYuXKlJWIkIiIiIiOYXLHr1q0bfv/9d6xfvx5OTk6YNGkSDh8+jN9//x1PPPGEJWIkIiIiIiOYXLEDgMceewzr1q0zdyxEZKSgcasMTk+e/lQ5R0JERBWJyRU7ALh58ya++uorTJgwAdevXwdw93myFy5cMGtwRERERGQ8kyt2Bw4cQEREBFxdXZGcnIwRI0bAw8MDy5cvx9mzZ7FkyRJLxFnu4uPjER8fj7y8PGuHQkRERGQUkyt2cXFxGDJkCI4fPw6dTqefHh0djcTERLMGZ0188gQRERFVNiYndrt378bzzz9fZHrNmjWRkpJilqCIiIiIyHQmJ3ZarRbp6elFph87dgxeXl5mCYqIiIiITGdyYte1a1dMmTJF/5QJRVFw9uxZvPHGG+jVq5fZAyQiIiIi45ic2M2YMQMZGRmoUaMGbt++jfbt26NevXpwdnbG+++/b4kYiYiIiMgIJl8V6+rqinXr1mHr1q1ISkpCRkYGmjVrhoiICEvER0RERERGMimxy8nJgYODA/bv3482bdqgTZs2loqLiIiIiExk0qFYOzs7BAYG8t5uRERERBWQyYdiJ06ciAkTJuCbb76Bh4eHJWKiBwgfjUVERGQ+Jid2c+fOxYkTJ+Dn54datWrBycmpUPu+ffvMFpw53Lx5ExEREcjNzUVubi5Gjx6NmJgYa4dFREREZHYmJ3bdunWDoiiWiMUinJ2dkZiYCEdHR2RmZqJx48bo2bMnqlevbu3QiIiIiMzK5MTu7bfftkAYlmNjYwNHR0cAQFZWFkQEImLlqIiIiIjMz+T72NWpUwfXrl0rMv3mzZuoU6eOyQEkJiaiS5cu8PPzg6IoWLFiRZF54uPjERQUBJ1Oh/DwcOzatcukPm7evImQkBD4+/tj7Nix8PT0NDlOIiIioorO5MQuOTnZ4FWxWVlZOH/+vMkBZGZmIiQkBPHx8Qbbly1bhri4OEyePBn79u1DSEgIIiMjcfnyZf08oaGhaNy4cZG/ixcvAgDc3NyQlJSE06dP47vvvkNqaqrJcRIRERFVdEYfiv3tt9/0/79mzRq4urrqX+fl5SEhIQG1a9c2OYCoqChERUWpts+cORMxMTEYOnQoAGDevHlYtWoVFi5ciHHjxgEA9u/fb1Rf3t7eCAkJwebNm/HMM88YnCcrKwtZWVn614aei0tERERUERmd2HXv3l3//4MHDy7UZmdnh6CgIMyYMcNsgQFAdnY29u7di/Hjx+unaTQaREREYPv27UatIzU1FY6OjnB2dkZaWhoSExPx4osvqs4/bdo0vPPOO2WOnYiIiKi8GZ3Y5efnAwBq166NPXv2lMtVpVevXkVeXh68vb0LTff29saRI0eMWseZM2cwcuRI/UUT//d//4cmTZqozj9+/HjExcXpX6enpyMgIKB0AyAiIiIqRyY/UqxOnTq4fv16pbldSMuWLY0+VAsAWq0WWq3WcgERERERWYhJiZ2dnR0OHDhgqViK8PT0hI2NTZGLHVJTU+Hj42PRvuPj4xEfH8/HpxHdh08LISKquEy+KnbgwIFYsGCBJWIpwt7eHmFhYUhISNBPy8/PR0JCAlq1amXRvmNjY3Ho0CHs3r3bov0QERERmYvJNyjOzc3FwoULsX79eoSFhRV5pNjMmTNNWl9GRgZOnDihf3369Gns378fHh4eCAwMRFxcHAYPHozmzZujZcuWmD17NjIzM/VXyRIRERHRXSYndv/88w+aNWsGADh27FihttI8amzPnj3o2LGj/nXBhQuDBw/GokWL0LdvX1y5cgWTJk1CSkoKQkNDsXr16iIXVJgbD8USERFRZWNyYrdx40azBtChQ4cSH/E1atQojBo1yqz9liQ2NhaxsbFIT08vdM8+IiIioorK5HPsiIiIiKhiMrliB9w9fPrjjz/i7NmzyM7OLtS2fPlyswRmbTwUS0RERJWNyRW7H374Aa1bt8bhw4fxyy+/ICcnB//++y82bNhQpQ5Z8qpYIiIiqmxMTuymTp2KWbNm4ffff4e9vT3mzJmDI0eOoE+fPggMDLREjERERERkBJMTu5MnT+Kpp+7eiNTe3h6ZmZlQFAWvvPIK5s+fb/YAiYiIiMg4Jid27u7u+O+//wAANWvWxD///AMAuHnzJm7dumXe6KwoPj4ewcHBaNGihbVDISIiIjKKyYldu3btsG7dOgBA7969MXr0aMTExKB///54/PHHzR6gtfAcOyIiIqpsTL4qdu7cubhz5w4AYOLEibCzs8O2bdvQq1cvvPnmm2YPkIiIiIiMY3Ji5+Hhof9/jUaDcePGmTUgIiIiIiodoxO79PR0o+ZzcXEpdTAVCe9jR0RERJWN0Ymdm5tbsc+CFREoilJlEiE+UoyIiIgqG6MTu3ufESsiiI6OxldffYWaNWtaJDAiIiIiMo3RiV379u0LvbaxscGjjz6KOnXqmD0oIiIiIjKdybc7ISIiIqKKyeSrYh8UvHiCqOIJGrdKtS15+lPlGAkRUcVUpopdcRdTVHa8QTERERFVNkZX7Hr27Fno9Z07d/DCCy/Aycmp0PTly5ebJzIiIiIiMonRid39t/wYOHCg2YMhIiIiotIzOrH7+uuvLRkHEREREZURr4olIiIiqiKMSuxeeOEFnD9/3qgVLlu2DEuXLi1TUERERERkOqMOxXp5eeHhhx9GmzZt0KVLFzRv3hx+fn7Q6XS4ceMGDh06hC1btuCHH36An58f5s+fb+m4LY63OyEiIqLKxqjE7t1338WoUaPw1Vdf4bPPPsOhQ4cKtTs7OyMiIgLz58/Hk08+aZFAyxufFUtERESVjdEXT3h7e2PixImYOHEibty4gbNnz+L27dvw9PRE3bp1q/Q97YiIiIgqg1I9ecLd3R3u7u7mjoWIiIiIyoBXxRIRERFVEUzsiIiIiKoIJnZEREREVQQTOyIiIqIqolSJXW5uLtavX48vvvgC//33HwDg4sWLyMjIMGtwRERERGQ8k6+KPXPmDJ588kmcPXsWWVlZeOKJJ+Ds7IwPPvgAWVlZmDdvniXiJCIiIqISmFyxGz16NJo3b44bN27AwcFBP71Hjx5ISEgwa3DWFB8fj+DgYLRo0cLaoRAREREZxeSK3ebNm7Ft2zbY29sXmh4UFIQLFy6YLTBr45MniIiIqLIxObHLz883+PzU8+fPw9nZ2SxBkXkFjVtlcHry9KfKORIiIiKyJJMTu86dO2P27NmYP38+AEBRFGRkZGDy5MmIjo42e4BERBUFfyQRUUVncmI3Y8YMREZGIjg4GHfu3MGzzz6L48ePw9PTE99//70lYiQiIiIiI5ic2Pn7+yMpKQnLli1DUlISMjIyMHz4cAwYMKDQxRREREREVL5MTuwAwNbWFgMGDMCAAQPMHQ8RERERlZLJtzuZNm0aFi5cWGT6woUL8cEHH5glKCIiIiIyncmJ3RdffIFGjRoVmf7www/z5sREREREVmRyYpeSkgJfX98i0728vHDp0iWzBEVEREREpjM5sQsICMDWrVuLTN+6dSv8/PzMEhQRERERmc7kiydiYmIwZswY5OTkoFOnTgCAhIQEvP7663j11VfNHqC53Lp1Cw899BB69+6Njz/+2NrhEBEREZmdyYnd2LFjce3aNbz00kvIzs4GAOh0OrzxxhsYP3682QM0l/fffx+PPvqotcMgIiIishiTEztFUfDBBx/grbfewuHDh+Hg4ID69etDq9VaIj6zOH78OI4cOYIuXbrgn3/+sXY4RERERBZh8jl2BapVq4YWLVqgcePGZUrqEhMT0aVLF/j5+UFRFKxYsaLIPPHx8QgKCoJOp0N4eDh27dplUh+vvfYapk2bVuoYiYiIiCoDkyt2mZmZmD59OhISEnD58mXk5+cXaj916pTJ6wsJCcGwYcPQs2fPIu3Lli1DXFwc5s2bh/DwcMyePRuRkZE4evQoatSoAQAIDQ1Fbm5ukWXXrl2L3bt3o0GDBmjQoAG2bdtmUmxERERElYnJid2IESPw119/4bnnnoOvry8URSlTAFFRUYiKilJtnzlzJmJiYjB06FAAwLx587Bq1SosXLgQ48aNAwDs379fdfkdO3bghx9+wE8//YSMjAzk5OTAxcUFkyZNMjh/VlYWsrKy9K/T09NLMSoiIiKi8mdyYvfnn39i1apVaNOmjSXiKSQ7Oxt79+4tdFGGRqNBREQEtm/fbtQ6pk2bpj8Mu2jRIvzzzz+qSV3B/O+8807ZAiciIiKyApMTO3d3d3h4eFgiliKuXr2KvLw8eHt7F5ru7e2NI0eOWKTP8ePHIy4uTv86PT0dAQEBFumLiMhYQeNWGZyePP2pco6EiCoykxO7d999F5MmTcLixYvh6OhoiZgsZsiQISXOo9VqodVqER8fj/j4eOTl5Vk+MCIiIiIzMDmxmzFjBk6ePAlvb28EBQXBzs6uUPu+ffvMFpynpydsbGyQmppaaHpqaip8fHzM1o8hsbGxiI2NRXp6OlxdXS3aFxEREZE5mJzYde/e3QJhGGZvb4+wsDAkJCTo+83Pz0dCQgJGjRpVbnEQERERVQYmJ3aTJ082awAZGRk4ceKE/vXp06exf/9+eHh4IDAwEHFxcRg8eDCaN2+Oli1bYvbs2cjMzNRfJWspPBRLRERElY3JiR0A3Lx5Ez///DNOnjyJsWPHwsPDA/v27YO3tzdq1qxp0rr27NmDjh076l8XXLgwePBgLFq0CH379sWVK1cwadIkpKSkIDQ0FKtXry5yQYW58VAsERERVTYmJ3YHDhxAREQEXF1dkZycjJiYGHh4eGD58uU4e/YslixZYtL6OnToABEpdp5Ro0bx0CsRERFRCUx+pFhcXByGDBmC48ePQ6fT6adHR0cjMTHRrMFZU3x8PIKDg9GiRQtrh0JERERkFJMrdrt378YXX3xRZHrNmjWRkpJilqAqAh6KJd43jIiIKhuTK3ZardbgY7aOHTsGLy8vswRFRERERKYzObHr2rUrpkyZgpycHACAoig4e/Ys3njjDfTq1cvsARIRERGRcUxO7GbMmIGMjAzUqFEDt2/fRvv27VGvXj04Ozvj/ffft0SMVsFz7IiIiKiyMfkcO1dXV6xbtw5bt25FUlISMjIy0KxZM0RERFgiPqvhOXZERERU2ZiU2OXk5MDBwQH79+9HmzZt0KZNG0vFRUREFQQvJCKqPEw6FGtnZ4fAwEA+jYGIiIioAjL5HLuJEydiwoQJuH79uiXiISIiIqJSMvkcu7lz5+LEiRPw8/NDrVq14OTkVKh93759ZgvOmvisWCIiIqpsTE7sunfvboEwKh5ePEFERESVjcmJ3eTJky0RBxERERGVkcmJHQDcvHkTP//8M06ePImxY8fCw8MD+/btg7e3N2rWrGnuGImI6AHFK3KJTGNyYnfgwAFERETA1dUVycnJiImJgYeHB5YvX46zZ89iyZIlloiTiIiIiEpg8lWxcXFxGDJkCI4fPw6dTqefHh0djcTERLMGZ0188gQRERFVNiYndrt378bzzz9fZHrNmjWRkpJilqAqgtjYWBw6dAi7d++2dihERERERjE5sdNqtUhPTy8y/dixY/Dy8jJLUERERERkOpMTu65du2LKlCnIyckBACiKgrNnz+KNN95Ar169zB4gERERERnH5MRuxowZyMjIQI0aNXD79m20b98e9erVg7OzM95//31LxEhERERERjD5qlhXV1esW7cOW7duRVJSEjIyMtCsWTNERERYIj4iIiIiMpJRiZ2HhweOHTsGT09PDBs2DHPmzEGbNm3Qpk0bS8dHREREREYyKrHLzs5Geno6PD09sXjxYnzwwQdwdna2dGxWxWfFmgdvLkpERFR+jErsWrVqhe7duyMsLAwigpdffhkODg4G5124cKFZA7QWPiuWiIiIKhujErtvv/0Ws2bNwsmTJwEAaWlpuHPnjkUDIyIiIiLTGJXYeXt7Y/r06QCA2rVr45tvvkH16tUtGhgRERERmcao2514eHjg6tWrAICOHTvC3t7eokERERERkemMSuwKLp4AgMWLF/MwLBEREVEFxIsniIiIiKoIky+eUBSFF08QERERVUC8eIKIiIioijD5kWKnT5+2RBxEREREVEZGJXaffPIJRo4cCZ1Oh08++aTYeV9++WWzBGZtfPIEERERVTZGJXazZs3CgAEDoNPpMGvWLNX5FEWpMokdnzxBRERElY1Rid29h195KJaIiIioYjLqPnZEREREVPEZVbGLi4szeoUzZ84sdTBEREREVHpGJXZ///13odf79u1Dbm4uGjZsCAA4duwYbGxsEBYWZv4IiYiIiMgoRiV2Gzdu1P//zJkz4ezsjMWLF8Pd3R0AcOPGDQwdOhSPPfaYZaIkIiIiohKZfI7djBkzMG3aNH1SBwDu7u547733MGPGDLMGR0RERETGMzmxS09Px5UrV4pMv3LlCv777z+zBEVEREREpjM5sevRoweGDh2K5cuX4/z58zh//jz+97//Yfjw4ejZs6clYiQiIiIiI5j8SLF58+bhtddew7PPPoucnJy7K7G1xfDhw/HRRx+ZPUBzCAoKgouLCzQaDdzd3QudM0hERERUVZic2Dk6OuKzzz7DRx99hJMnTwIA6tatCycnJ7MHZ07btm1DtWrVrB0GERERkcWYnNgVcHJywiOPPGLOWIiIiIioDKz+5InExER06dIFfn5+UBQFK1asKDJPfHw8goKCoNPpEB4ejl27dpnUh6IoaN++PVq0aIGlS5eaKXIiIiKiiqXUFTtzyczMREhICIYNG2bw4otly5YhLi4O8+bNQ3h4OGbPno3IyEgcPXoUNWrUAACEhoYiNze3yLJr166Fn58ftmzZgpo1a+LSpUuIiIhAkyZNWG0kIiKiKsfqiV1UVBSioqJU22fOnImYmBgMHToUwN2LN1atWoWFCxdi3LhxAID9+/cX20fNmjUBAL6+voiOjsa+fftUE7usrCxkZWXpX6enp5syHCIiIiKrMflQbGZmpiXiMCg7Oxt79+5FRESEfppGo0FERAS2b99u1DoyMzP199fLyMjAhg0b8PDDD6vOP23aNLi6uur/AgICyjYIIiIionJicmLn7e2NYcOGYcuWLZaIp5CrV68iLy8P3t7eRWJISUkxah2pqalo27YtQkJC8Oijj2LQoEFo0aKF6vzjx49HWlqa/u/cuXNlGgMRERFReTH5UOy3336LRYsWoVOnTggKCsKwYcMwaNAg+Pn5WSK+MqtTpw6SkpKMnl+r1UKr1VowIiIiIiLLMLli1717d6xYsQIXLlzACy+8gO+++w61atXC008/jeXLlxu8iKG0PD09YWNjg9TU1ELTU1NT4ePjY7Z+DImPj0dwcHCx1T0iIiKiiqTUtzvx8vJCXFwcDhw4gJkzZ2L9+vV45pln4Ofnh0mTJuHWrVtlDs7e3h5hYWFISEjQT8vPz0dCQgJatWpV5vUXJzY2FocOHcLu3bst2g8RERGRuZT6qtjU1FQsXrwYixYtwpkzZ/DMM89g+PDhOH/+PD744APs2LEDa9euLXE9GRkZOHHihP716dOnsX//fnh4eCAwMBBxcXEYPHgwmjdvjpYtW2L27NnIzMzUXyVLRERERHeZnNgtX74cX3/9NdasWYPg4GC89NJLGDhwINzc3PTztG7dGg899JBR69uzZw86duyofx0XFwcAGDx4MBYtWoS+ffviypUrmDRpElJSUhAaGorVq1cXuaDC3OLj4xEfH4+8vDyL9kNERERkLiYndkOHDkW/fv2wdetW1fPP/Pz8MHHiRKPW16FDB4hIsfOMGjUKo0aNMjXUMomNjUVsbCzS09Ph6uparn0TERERlYbJid2lS5fg6OhY7DwODg6YPHlyqYMiMoegcatU25KnP1WOkRAREZUPkxO73Nxcg09jUBQFWq0W9vb2ZgnM2ngoloiIiCobkxM7Nzc3KIqi2u7v748hQ4Zg8uTJ0GhKfdGt1fFQLD2oWOkkIqq8TE7sFi1ahIkTJ2LIkCFo2bIlAGDXrl1YvHgx3nzzTVy5cgUff/wxtFotJkyYYPaAiYiIiMgwkxO7xYsXY8aMGejTp49+WpcuXdCkSRN88cUXSEhIQGBgIN5//30mdkRERETlyORjpdu2bUPTpk2LTG/atCm2b98OAGjbti3Onj1b9uisiE+eICIiosrG5IpdQEAAFixYgOnTpxeavmDBAgQEBAAArl27Bnd3d/NEaCU8x46I6MHG802pMjI5sfv444/Ru3dv/Pnnn/pq1p49e3DkyBH8/PPPAIDdu3ejb9++5o2UqAop6QtDrZ1fJkREVByTE7uuXbvi6NGj+OKLL3D06FEAQFRUFFasWIGgoCAAwIsvvmjWIImIiIioZCYldjk5OXjyyScxb948TJs2zVIxVQi8jx0RERFVNiZdPGFnZ4cDBw5YKpYKJTY2FocOHcLu3butHQoRERGRUUy+KnbgwIFYsGCBJWIhIiIiojIo1SPFFi5ciPXr1yMsLAxOTk6F2mfOnGm24IiIiIjIeCYndv/88w+aNWsGADh27FihtuIeNUZERERElmVyYrdx40ZLxEFU7nhLESIiqmpMPseuwIkTJ7BmzRrcvn0bACAiZguqIuCTJ4iIiKiyMTmxu3btGh5//HE0aNAA0dHRuHTpEgBg+PDhePXVV80eoLXwqlgiIiKqbExO7F555RXY2dnh7NmzcHR01E/v27cvVq9ebdbgiIiIiMh4Jp9jt3btWqxZswb+/v6FptevXx9nzpwxW2BEREREZBqTK3aZmZmFKnUFrl+/Dq1Wa5agiIiIiMh0Jid2jz32GJYsWaJ/rSgK8vPz8eGHH6Jjx45mDY6IiIiIjGfyodgPP/wQjz/+OPbs2YPs7Gy8/vrr+Pfff3H9+nVs3brVEjFaBZ8VS0RERJWNyRW7xo0b49ixY2jbti26deuGzMxM9OzZE3///Tfq1q1riRitglfFEhERUWVjcsUOAFxdXTFx4kRzx0JEREREZVCqxO7mzZvYtWsXLl++jPz8/EJtgwYNMktgRJUdn2xROtxuRESlZ3Ji9/vvv2PAgAHIyMiAi4tLoefDKorCxI6IiIjISkw+x+7VV1/FsGHDkJGRgZs3b+LGjRv6v+vXr1siRiIiIiIygsmJ3YULF/Dyyy8bvJcdEREREVmPyYldZGQk9uzZY4lYiIiIiKgMTD7H7qmnnsLYsWNx6NAhNGnSBHZ2doXau3btarbgiIiIHlS8kIhKw+TELiYmBgAwZcqUIm2KovCGvkRERERWYvKh2Pz8fNW/qpTUxcfHIzg4GC1atLB2KERERERGMTmxe1DwyRNERERU2Rid2EVHRyMtLU3/evr06bh586b+9bVr1xAcHGzW4IiIiIjIeEafY7dmzRpkZWXpX0+dOhV9+vSBm5sbACA3NxdHjx41e4BEZH48KZuIqGoyumInIsW+JiIiIiLr4jl2RERERFWE0YmdoiiFngtbMI2IiIiIKgajz7ETEQwZMgRarRYAcOfOHbzwwgtwcnICgELn3xGRdfEcOiKiB5PRid3gwYMLvR44cGCReQYNGlT2iIiIiIioVIxO7L7++mtLxkFEZHWsdBJRZceLJ4iIiIiqCJOfFVsZnT59GsOGDUNqaipsbGywY8cO/bmBREREVBir15XXA5HYDRkyBO+99x4ee+wxXL9+XX8BCBEREVFVUuUTu3///Rd2dnZ47LHHAAAeHh5WjoiIiIjIMqx+jl1iYiK6dOkCPz8/KIqCFStWFJknPj4eQUFB0Ol0CA8Px65du4xe//Hjx1GtWjV06dIFzZo1w9SpU80YPREREVHFYfWKXWZmJkJCQjBs2DD07NmzSPuyZcsQFxeHefPmITw8HLNnz0ZkZCSOHj2KGjVqAABCQ0ORm5tbZNm1a9ciNzcXmzdvxv79+1GjRg08+eSTaNGiBZ544gmLj42IyhfPCyIyDv+tVF1WT+yioqIQFRWl2j5z5kzExMRg6NChAIB58+Zh1apVWLhwIcaNGwcA2L9/v+ryNWvWRPPmzREQEAAAiI6Oxv79+1UTu6ysrEI3W05PTzd1SEQPtIr8hWHJ2NTWba71Ez0oKvJnSGVg9cSuONnZ2di7dy/Gjx+vn6bRaBAREYHt27cbtY4WLVrg8uXLuHHjBlxdXZGYmIjnn39edf5p06bhnXfeKXPsRET3svSXVXHrt3TSWdLYrPlFzb4rTt/l0T9/YFXwxO7q1avIy8uDt7d3oene3t44cuSIUeuwtbXF1KlT0a5dO4gIOnfujKefflp1/vHjxyMuLk7/Oj09XV/tIyIiosqrpMSvKlQLK3RiZy4lHe69l1arhVarRXx8POLj45GXl2fh6IiI6EFUFZIIqngqdGLn6ekJGxsbpKamFpqempoKHx8fi/YdGxuL2NhYpKenw9XV1aJ9ERERPSiY0FpWhU7s7O3tERYWhoSEBHTv3h0AkJ+fj4SEBIwaNcq6wRER0QPtQT2PjSo2qyd2GRkZOHHihP716dOnsX//fnh4eCAwMBBxcXEYPHgwmjdvjpYtW2L27NnIzMzUXyVrKTwUS0RU8bH6Q1SY1RO7PXv2oGPHjvrXBRcuDB48GIsWLULfvn1x5coVTJo0CSkpKQgNDcXq1auLXFBhbjwUS0RERJWN1RO7Dh06QESKnWfUqFE89EpERGbHil/pcLtVXFZP7IiIiIiMVRHPbaxICa3VnxVbUcXHxyM4OBgtWrSwdihERERERmFipyI2NhaHDh3C7t27rR0KERERkVF4KJaIiKyGt+0gMi9W7FTwUCwRERFVNkzsVPBQLBEREVU2TOyIiIiIqggmdkRERERVBBM7IiIioiqCiZ0KXjxBRERElQ0TOxW8eIKIiIgqG97HjoioAqgMjyoiooqPFTsiIiKiKoKJHREREVEVwcROBS+eICIiosqGiZ0KXjxBRERElQ0TOyIiIqIqgokdERERURXBxI6IiIioimBiR0RERFRFMLEjIiIiqiKY2Kng7U6IiIiosmFip4K3OyEiIqLKhokdERERURXBxI6IiIioimBiR0RERFRFMLEjIiIiqiKY2BERERFVEUzsiIiIiKoIJnZEREREVQQTOyIiIqIqgomdCj55goiIiCobW2sHUFHFxsYiNjYW6enpcHV1tXY4RFVG0LhVBqcnT3+qnCMhIqp6WLEjIiIiqiKY2BERERFVEUzsiIiIiKoIJnZEREREVQQTOyIiIqIqgokdERERURXBxI6IiIioimBiR0RERFRFVPnE7ujRowgNDdX/OTg4YMWKFdYOi4iIiMjsqvyTJxo2bIj9+/cDADIyMhAUFIQnnnjCukERERERWUCVr9jd67fffsPjjz8OJycna4dCREREZHZWT+wSExPRpUsX+Pn5QVEUg4dJ4+PjERQUBJ1Oh/DwcOzatatUff3444/o27dvGSMmIiIiqpisnthlZmYiJCQE8fHxBtuXLVuGuLg4TJ48Gfv27UNISAgiIyNx+fJl/TyhoaFo3Lhxkb+LFy/q50lPT8e2bdsQHR1t8TERERERWYPVz7GLiopCVFSUavvMmTMRExODoUOHAgDmzZuHVatWYeHChRg3bhwA6M+hK86vv/6Kzp07Q6fTFTtfVlYWsrKy9K/T09ONGAURERGR9Vm9Ylec7Oxs7N27FxEREfppGo0GERER2L59u0nrMvYw7LRp0+Dq6qr/CwgIMDluIiIiImuo0Ind1atXkZeXB29v70LTvb29kZKSYvR60tLSsGvXLkRGRpY47/jx45GWlqb/O3funMlxExEREVmD1Q/FlgdXV1ekpqYaNa9Wq4VWq7VwRERERETmV6Erdp6enrCxsSmSlKWmpsLHx8eifcfHxyM4OBgtWrSwaD9ERERE5lKhEzt7e3uEhYUhISFBPy0/Px8JCQlo1aqVRfuOjY3FoUOHsHv3bov2Q0RERGQuVj8Um5GRgRMnTuhfnz59Gvv374eHhwcCAwMRFxeHwYMHo3nz5mjZsiVmz56NzMxM/VWyRERERHSX1RO7PXv2oGPHjvrXcXFxAIDBgwdj0aJF6Nu3L65cuYJJkyYhJSUFoaGhWL16dZELKswtPj4e8fHxyMvLs2g/REREROZi9cSuQ4cOEJFi5xk1ahRGjRpVThHdFRsbi9jYWKSnp8PV1bVc+yYiIiIqjQp9jh0RERERGY+JnQpeFUtERESVDRM7FbwqloiIiCobJnZEREREVQQTOyIiIqIqgomdCp5jR0RERJUNEzsVPMeOiIiIKhsmdkRERERVhNVvUExERNYVNG6VwenJ058q50iIqKxYsVPBc+yIiIiosmFip4Ln2BEREVFlw8SOiIiIqIpgYkdERERURTCxIyIiIqoimNgRERERVRFM7FTwqlgiIiKqbJjYqeBVsURERFTZMLEjIiIiqiKY2BERERFVEUzsiIiIiKoIJnZEREREVQQTOxW8KpaIiIgqGyZ2KnhVLBEREVU2TOyIiIiIqggmdkRERERVhK21A6joRAQAkJ6ebtF+8rNuGZyenp6u2lZSe0HMlmxnbIyNsTE2xsbYKkpspe3bXLFZSsH6C3KS4ihizFwPsPPnzyMgIMDaYRAREdED7ty5c/D39y92HiZ2JcjPz8fFixfh7OwMRVEs3l96ejoCAgJw7tw5uLi4mNRelmUt3c7YGNuDEhu3C2NjbNbv29qxmZuI4L///oOfnx80muLPouOh2BJoNJoSs2NLcHFxKXZnKa69LMtaup2xMbYHJTZuF8bG2Kzft7VjMydXV1ej5uPFE0RERERVBBM7IiIioiqCiV0Fo9VqMXnyZGi1WpPby7KspdsZG2N7UGLjdmFsjM36fVs7NmvixRNEREREVQQrdkRERERVBBM7IiIioiqCiR0RERFRFcHEjoiIiKiKYGJHREREVEXwyRMPsEOHDmHu3LnYvn07UlJSAAA+Pj5o1aoVRo0ahRo1amDhwoVF2lu3bo0hQ4bA1dUVK1asMNjepk0bfPvtt6rLenl5WWfQsOy4u3XrhvT0dNXlW7Vqhe+++84qfVt7uxMRWdquXbsMfr62bNkSKSkp2LlzZ6G28PBw+Pj4AAByc3Px77//FmoPDg6GnZ0dAJS4fEXB251Y2dWrV8v0Rb5//37VnTg7O1t1Wa1Wi969e6NZs2aIjIyEt7c3ACA1NRXr1q3D7t27odVq4erqioiIiELtCQkJ+O+//1CtWjVcu3YN4eHhhdq3b9+O3NxceHl54cknnyyy7K1bt7Bs2TLs27ev1OO2t7cv9h+w2tgdHBwwY8YMi4x7586d8PT01M9z//IrV67EtWvX0KRJE/Ts2bNc+y7Y7itXrsS5c+dKvV2Bkj/c1JavU6dOmfb10r7nBcsDKPW/pbKM25jYyrrdrdl3cWMPDQ21Wt8P8v5mre8VY35clmW7F7c/tWrVCv3798fWrVsRGBhY6PPvzJkz8PLywvXr16EoCjw8PAAA169fh4igX79+qFmzJubPn4+0tDTcy9XVFSNHjsT58+exbNkyg8v3798fX3zxBRwdHVERMLGzot27dyMyMhKOjo4mf5Hv2LEDiqIgKyuryE589uxZNGvWDFevXkVqaqrBJCA/Px8xMTGYO3euwdj8/f1x+/ZtXL16FYqiFGoTEQQEBODWrVtITk4u8py8Fi1a4MqVK2jQoAHWrl1bZNmePXvi999/h4+PT6kSGB8fH3h6emLv3r0mjz0xMRFubm7YtWsX6tWrZ9Zxp6enIygoCE5OTjh79myR5UNCQmBrawt7e3ts3769XPsu+PBZsWIFFEUp1XZ99NFH4efnp1/H/R9uPXr0wMWLF7F9+3aDH6y2trbw8vJC586dy/U937lzJ2rUqAERwZUrV0z+t1SWcRsTW1m2e35+PqpXr46rV6+We9/GvOc6nQ4AuL+V4/72yCOP4MyZM3Bycir375WSflyWZbuXtD/Z2NigQYMG+P7779GwYUPcq3fv3li5ciWaN2+OTZs2wcbGBgCQl5eHhIQE9OvXD7dv38bs2bOL/OBfu3YtRo8eDQcHB3z//feIiIgosvz//d//oV27dvjyyy9RIQhZTXh4uIwcOVLy8/OLtOXn50vNmjXF3d1d0tLSirR37dpV3N3dpU2bNkXajhw5Iq6uruLn52dw2bS0NNFoNAaXLaDVasXe3l61XafTqbbrdDr59ddfxcHBwWD7I488IjY2NqUad1pamvj6+oqrq6scOXKkSHtJY9fpdNKpUyfp3LmzwdjKMu6C5bVareqyf/zxh+h0unLvW0SkVatWotFoSr1da9SoIdWqVZPVq1dLbm6uvi03N1fWrFkj1apVkxo1ahhcPiQkRLy9vaVXr15F2iz9nqelpYmXl5d4eXmV6t9SWcZtTGxl2e5t27YVnU4nzzzzTLn3XdLYW7duLR4eHtKtW7dy7/tB3t+qVasmtWvXtsr3iru7u/j7+1vks72k/cnGxkYeffTRIm0iIm5ubrJgwQKpVq2awXYPDw9xcnIy2CYiUq1aNXF3d1dt37Jli7i5uam2lzcmdlak0+nk8OHDxbarfZFXq1ZNli1bppo8abVa1TYRkdq1a4utra1qu4eHh/j4+Ki2u7m5iZeXl8G2oKAgGTNmjPj6+hpst7OzEz8/P9V1l5TAODo6FpvAFDf2Ro0ayWuvvabaXpZxi4jUqFFDXF1dVfvu37+/1KpVq9z7FhGxt7dXfU9ESt6uzs7Oxe5TDg4Oqh+cOp1Ofv7552LbLfWeF6xfbfmS/i2VZdzGxFaW7V6tWjVZuHCh6peKJfsWKX7sDg4O8sMPP6i2c3+z3P7m6OhYbGyW+l4p6cdlWbZ7SfuTm5ub6rpdXFxk3rx5Ur16ddW4Skrs1H6Qi4js2rVLXFxcVNvLGy+esKKCc2MaNWpksF2n0+lP2ryfVqvFyZMn4ebmZrC9WrVqyM7OVu27d+/e+PDDD9G1a1eDJfO0tDRoNBqMHj0ajz/+eJH2jIwM2NvbY9asWUXamzVrhtmzZ6Nly5b47bffiiybm5uLiIgI1diKGzcA2NjY6Evypo59ypQp6N+/P2xtbfHJJ5+YddwF7VlZWQaXr127Nr7//ns0adKk3PtOSEhAdna2/tyf0mzX/Px82Nvbq7bb29sjJyfHYJuPjw/27Nmj+lxFS77nAODo6AhROeukpH9LZRm3MbGVZbtrtVpkZ2cjPz+/3PsGih+7m5sbTp06pfqec39zM9he1v3N3d0dGRkZqu2W/F5xdXVFVlZWqfoGit/uJe1PrVu3xpo1a/DLL7/g8ccf15+ukp6ejpCQELz88svo0aNHkeX+/vtv2Nvbw8nJCVevXoWnp2eh9qtXr8LZ2RmZmZn4+++/0bRp0yLLv/jii+jSpYvquMobEzsreu211zBy5Ejs3bvX5C/yRo0a4c0338QzzzyD9PT0QjtxQfKUnZ2tmgR89dVXGDZsGDIzMzFjxowiJ6lu3rwZZ8+exaxZs/DZZ58hLy8PwN1/eGFhYVi6dClOnz6Njz76CK+++qr+nC4RgY+PD5599lmcPHkSvXr1KrLs0KFDsXTpUri5uZUqgcnOzoZGozH4D7iksZ8/fx6Ojo4IDAy0yLgnT56MWrVqqS4/ZcoU/Pvvv1bpu1evXli2bBkaNGhQqu0K3P1SUftws7e3R25ursHln3jiCUyfPh2PPPKIwWTfku95wYUjAEr1b6ks4zYmtrJs906dOuHll1/Go48+WqrPAUu+523btsWbb76J1q1b48CBA9zfyml/y87Oxu3bt0v1A7Gs3ysl/bgsy3YvaX/avn07QkND0a9fP+Tm5uqT4+zsbNjY2MDHxwc//vgj1q1bhxo1agAALl++jJs3b6Jdu3ZITU2Fr68vmjRpUmjdBw8eRIMGDdCwYUOEhYXB3d29yPKRkZGq56tbAy+esLJly5Zh1qxZ2Lt3b5Ev4ri4OJw+fRpz5sxBSkpKoS9yb29vBAUFYd++fUV2YltbWwwfPhw1a9ZEfHx8kWV9fHwwZswYvP7660bFmJOTg6tXrwIAPD09i/ziOn36dKEkpXbt2iUuW9px+/j4YNSoUTh37hwWLlxo0bGXZdzGLG+Nvj/44INSb9eBAwfi/PnzWLduncEPt4iICNSsWRNLly41uHz79u2RlpaGffv2WeU9F5FS/Vsq67hLiq0s2/3GjRvw9/dHampqufdtzNibNWuGM2fOcH8r5/2tdevWmDt3rlW+Vwp+XFris724/angcz09PR179+4t9PkYFhYGFxcXHDlyxODVto0aNUJ+fj7WrFmDHTt2FGnv3LkzNBpNsctXJEzsKojSfpEXtxOXtGxFUJYEprKP3ZrKsl1L+nAraXlrv+el/bdU1nFbcrtbs29jluf+Zp39zZrfK5bc7vxcLx4TOzJowoQJSElJwcKFCw22f/bZZ7h69SomTZpksP3XX39FWloaBg0aZPKy1mTJcZe0vDX7JiKqyvbs2YNbt26hXbt2RdouXbqEnJwcBAYGlmrdZV3e7MrhAg0qpfj4eHnnnXdU21esWCGLFy822LZ7927566+/SrWsiMigQYOkY8eOqu2dOnWS2rVrq7Y3bNhQNBpNqZYty7hFyjZ2S467pOWt2bdI2bfrxYsX5cyZM6Va3prveUntlhy3MbGVpX9r9l3S8tzfrLO/WfN7xZLbvaRlGzVqpPr5WFybOdrLGxO7CqwsX+Ql7WglJQHWVNYEpjKP3ZrKul3L0m7t97ws/5bK2m7J7W7Nvktq5/5mnf3Nmt8rltzuJS174cIFSU5ONti2a9cu2bRpk+qyv/zyiyxatEi1vaTlyxsPxVZRFy9eRE5ODmrVqmXtUMrdgzx2Syppu+7evRu3bt1C+/btS7W8JWOz5LotPe6y9G/Nvsvav7W3e1lis+S6q+q4K0P/lQUTuwdYSc8TLO6B8Tdu3MDvv/+OgQMHQqPRFGm/cuUKEhIS0LlzZ3h4eODq1atYsGABsrKy0Lt3bzz00EMWG1dJLDnu/Px8HDx4ENevX0dISEiRsUdERGDLli1W6dva252IyNJKesauOaWmpuofv1aRMLGzMmMe0K0mMzMTe/fuNXgyKHD3MvDk5GQEBATA1tYW2dnZ+OWXX5CVlQVfX1/07dtX9Tm1t27dwpo1a9C8eXOD6966dSvatm0LnU4HFxcXPP/885g8eTJsbGywa9cuPPHEE0hPT4e7uzvWrVuH3r17w9bWFvn5+bh48SI2btyIM2fOWGTcxY398OHDmDNnjupzFMsybgD4888/ER0dDUVR4ObmVmjst27dwoULF1CjRg1ERUWVa98F233Lli3Izc0t9uHjanJzc3Hx4sUSP8Ty8vL0MQHAzp07kZWVhbCwMKxatarc3/OsrCxER0fD09OzxAfOW2LcrVq1gp2dXan6NrZ/a/ZdUv9///0397dy3t9ExCrfK9HR0XBxcbFK3+3bt8fEiRPxww8/GHzGbv/+/fHFF1/A0dGxyHqTkpLQtGlTDBs2DB4eHhg2bFihW5icPXsWLVu2hFarRYcOHfDll1/ilVdeweeffw5FUdC2bVv8/vvvRZ7hbTXWOP5Ldx0/flzq1KkjOp1O2rdvL3369JE+ffpI+/btRafTSb169eT48eMGl83OzpbBgwcLAGnRooUsWLCgUPuWLVsEgGg0GqlXr56cOnVKwsLCxMnJSRwdHcXW1lb69u2r+ky/IUOGSPPmzSUtLc3g3zPPPCMA5KeffpIvv/xSatWqJU899ZRkZWVJRESEDBgwQADIRx99JP7+/jJixAj9+nv16iWOjo6lGrfI3RNoAUjdunVNHrtGoxGdTidHjx41+7hFRNq1aycAJD09vcjYw8PDpUGDBgafdWjpvkVE+vXrJ9WrVxdFUaRWrVrSsmVLadmypdSqVUsURZG2bdtKamqq6nbfv3+/KIoijz/+uPTu3VvWr19fqP3gwYOi1WrFxsZG2rVrJ9evX5ennnpKFEURAGJnZ2eV99zR0VE8PDwkLCysVGMvy7gVRZE6depIy5YtLbLdL168KC1bthQA5d63Me+5g4MD97dy3t+CgoKkVq1aVvlecXd3l4CAAIts9yNHjoi/v78AMNi3VquVoKAg1WfsNmjQoNDn4b2mTp0qAOSpp57SP3/522+/1bcPGzZMAMgnn3wiHTp0kG7duknjxo1ly5Yt8tdff0lwcLBMmDBBdVzljYmdFUVEREi3bt1UH4jcrVs31YfVT548Wf8lPXHiRHF1dZWRI0fq25988kkBIAcOHJAxY8bIQw89JN26dZPs7Gy5c+eOaDQa6dKli2psBR+OGo3G4B8AURRFP/+VK1ekZcuW0rlzZ3F3d5fExETRaDSSnZ0tGo1Gdu7cqZ83PDxcdDpdqcYtIvL888/rk0ZTx67T6aRDhw4ycOBAs4/7zp074ubmpm+/f+wFz6+sWbNmufctIvL444+Lvb296sPDW7dubfBh8gVef/11ASCxsbEycOBAsbe3l6lTp+rbCxLP3377Tfr27SutW7eWxx57TM6fPy9t27YVDw8PiYmJKbJeS7/nd+7cEV9fX/H09CzV2Msy7jNnzkj16tXFx8fHItv9ueeekxYtWoiiKOXed0ljj4qKEmdnZ3n22WfLve8HeX9zdXWV2rVrW+V7xdPTU/z9/S3y2d6tWzeJjIwUAAb7trOzk8jISIPr7dGjh7Rr107s7OykR48eRf6qVatW6LN12bJl4uTkJF999ZWIiNSsWVPffuHCBVEURX7//Xf9/CtXrpSGDRuqjqu8MbGzIgcHBzl48KBqu4uLiwAQd3f3In8ajUYcHBz0VwEdP35c6tWrJ0OGDJH8/Hz9P04RkYyMDFEURTZv3qxft6+vr3h4eBQbm5ubm2zatMngn729faF/CCIi6enp0qpVK7GxsdEndiJ3H6B88uRJ/XwlPQi6uHEXjP3evk0Ze1BQkLz55psSGBho9nF36tRJHB0dC7XfO/agoCCZOXOm6sOkLdm3iIiTk1Ox273gqrKmTZsa/Lu//61bt4qXl5e89dZbIiLi7e2tb7927ZooiqKvNjg4OMhXX30lderUMdi3Jd/zgv59fHwM9t20adNix16WcRf0rZbMl3W7+/r6yqpVq0Sj0ZR73yWNvVq1ajJv3jzV95z7m2X2N61WK/7+/gb7LmnsZf1e0Wq1quMu63b38vKS9evXi0ajMdi3k5OTeHt7G+zX1tZWWrduLXZ2djJkyJAifzY2NkU+Wzds2CDVqlWTzz//XLRabaF2R0fHQkd9kpOTxdHRUXXc5Y3PirUiNzc3JCcno3Hjxgbbs7Ky4OTkhFmzZhVpGzlyJJ5//nl88sknAIB69eph06ZN6NSpE5577jlkZmbq53VycoKTkxN8fX3101566SW89dZbqs/0y8rKQmRkpOqVV35+fkhOTi40zdnZGWvXrkWNGjUwZMgQ/fQffvihUN/VqlVTfUB2SeMGgJiYGOTm5upfmzL21157Da+++ipyc3MNPkOyLOPu3Lmz/hE6hsb+2muvIS4uDvb29uXeNwDY2toaPL+kwKlTp2BnZ4du3boZbD948GCh161bt8aGDRsQERGBnJwcpKWl6ds8PDzg6Oiov3rNzc0NGo0Gly5dMrhuS77nwN2Hj1+7ds3gug8dOoQOHTrg7NmzBsdelnEX9H3lyhWDfQNl2+7Xrl3TnxRe3n2XNHatVgt3d3fV95z7m2X2N1dXV1y+fNlg3yWNvazfK66urvqnTZjaN1D8dr9+/br+HDZDfT/xxBP49ddfDT5jNygoCBcuXECfPn3w9ddfF+l35cqVRd6vjh07YuXKlXj66adhZ2eH7OxsfVu3bt3g5uamf52RkQGtVqs67nJn7czyQfbWW2+Ju7u7zJw5U5KSkiQlJUVSUlIkKSlJZs6cKba2tqql5dq1a8sXX3xR5L49Fy5ckAYNGujPbSnw2WefSXp6uv713r17xdXVVcLDw8XW1lZ/foatra2Eh4fLyJEjZc6cOaqxDx8+XIKDgw22jR8/XurVq6d6T6E2bdqInZ1dqcYtUrgsXpqxT5s2Tezs7Mw+7vT0dNXYCnTv3l3c3d2t0nezZs3E0dFRli9fXuhQSVpamixfvlzs7e2lffv2qsv7+PgYXP+///4r3t7e4uTkVKj9jTfekGvXronI3X3dxcVFnJycrPKe9+7dWzQajcGx161bVzw8PGTUqFFmH7eISJ8+fVT7Lut212g0Eh0drf+3Vp59lzT2l156SXx9fcXZ2Zn72z1jt/T+NnLkSFEUxSrfKyNGjCh13yLFb3c7Oztp3ry5Prb7+964caO+2unh4SGNGjWSRo0aiYeHhwCQwMBAuXHjhsF+O3XqJK6urgbbNm7caLCid6+vv/5aWrdurdpe3pjYWdn06dPF19dXFEXRn0ulKIr4+vpKZGSkvP322waXGz58uPTt21eGDBlSpO38+fPi4uJS7I44bdo0iY6OFpG752NdvHhRLl68KNnZ2QbnP3funOTl5elfX79+Xf755x/V9aenp6vesDEzM1Pee++9Uo1b5O5FAPXr1zfYZsrYy3vcInfHfufOHav0ff36dYmJiRF7e3v9RSQ6nU40Go3Y29tLkyZNVL9sRES6dOmieojnn3/+MXio+F7dunXTz1Pe7/m7774rgYGBBseu0WikSZMmcufOHYuMe9asWVKzZk2LbPeOHTvqLwoq775LGvudO3ekbdu2+veb+1v57G9z586VunXrWu17pWHDhhb5bB84cKB4eXmp7usFn+uHDx+WhQsXytSpU2Xq1KmycOFCSUpKkszMTNV+N23aVOg8xvutWLFC+vfvr9r+xx9/yMaNG1XbyxsTuwri1KlTsm3bNtm2bZucOnWqxPmTk5Nl9erVqu0XLlwo9k7Zp06dkosXLxodn7Ozc6Hzte7XuHFjOXv2rMntp06dEkdHR/npp5+MGrdI+Y7dUuM2Zv2W6jstLU02bNggOp1OZs2aJRs2bDB4svP9kpKSZOHChartBw8eLPZDe+fOnXLw4EGT93UR873nBWP/7rvv5LvvvjNq7OYad2n6Lmv/1uy7oP/t27dbre8HeX8riMNa3yvlvd1N/U6Ljo5WnX/atGmq1T1j2l988UW5cuWK0bGYGxO7SqSkL/ridtTi2kRKThLuPxHfnO0lLVuWcRvTXtzYLTnuktqt2bdIyftESR9eJbUX975Y8z0vqd2S4zYmtrL0b82+S+qf+5t19jdrfq9YcruX5TutpLjK2m5pRW9fTxWWlHAv6cTERNy+fdvkNgBITk5GTk5OmeKzlLKM25j2ijx2ayppu3z77bdIT08vdXtx74u13/Pi2i05bmNiK0v/1uy7pP65v1lnf7Pm94olt3tZPtdLiqus7ZbGxI6ISqWif7hZirXHbc3+H9S+relBHTeVHhM7IiIioiqCiR0ZRVGUSrnusrJ0bMWt35p9ExFR5cTErhKx5hexseX+4s71KO26K8O4AfOP3Zp9WxuTTipPD/L+Zs2xP8jb3ZKY2FUi5fFFnJ6ejhUrVuDw4cOFph86dKjQ3c3v98UXX8Db2xvu7u76u5536tQJN2/eLNRuyJ9//omaNWuqrrs8xv3ff/+VadwADI69uHEDd8fu7Oxstb6L2+7WVJGTTqp6HuT9zZpjt2bfEyZMgIeHh8G2xx57DA4ODgCAnj176n84L1myBFlZWYXaDRk4cKD+KRlWYfHrbslkubm58vfff8v169cLTd+8ebPqTS1FRKZOnSo3btwQd3d3/eXvQ4cOlfT0dH3b/Xr37i2ffvqpLF26VK5cuSL169fXP5Xh559/FhGR9evXy1NPPSV16tSROnXqSHBwsLz00ksyZ86cIn86nU7Gjx8vc+bMEUVR5PLly+U2bhExeuwF4xa5e9fwevXqlWncpoz93r5v3bpV5m1e1u0uIgb3jaVLl0pGRobqMi+88IJ+Wzdt2lT/vr3zzjuSmZlZqN2QqVOnytWrV8vtPVcbm6ljN8e4C/osy3Y31Pf98Vmqb5GyjZ37m3r7/cy1zcvze+X+cVmi73vb77Vo0SJZuXKlvu+xY8eKq6urtGrVSpKTkyU3N1d++uknmTJlikyZMkUWL14s165dk7S0NIN/dnZ2cvToUUlLSxONRiOpqamq8VYUisgD/FOlghgzZgyaNGmC4cOHIy8vD+3bt8e2bdvg6OiIlStXokOHDjh69Cg+/fRTHD58GFeuXEFAQACefvppgxWXvn37Ys6cOfDx8UGPHj2QkpICLy8vg337+PhgzZo1CAkJwXfffYfJkycjKSkJixcvxvz58xETE4PRo0fjmWeeQatWrQDc/aWTmZmJ6tWrw9nZGVeuXMGtW7fg5uaGjIwM5OTkQKPRID8/H23atIG9vb3Bvh955BGTxg3cfQaj2rhNGbu5xw0AN27cgKIo0Gq1yMrKQuvWrQ2OfevWrdi1a5dV+gaAyMhIBAUFoW/fvgCAPn364H//+x98fHzwxx9/ICQkBDdu3MCCBQv0293V1RU9evSAq6urwXW2bNkSv//+O7y9vdG0aVNcunQJNWrUKDKfqfu6Od9zAPjggw+KHXtgYKBFxm1M36Zud2v2bUr/1uz7Qd7frPm9Yu7tbkrfDRs2xOeff45OnTph+/btiIiIwKxZs7By5UrcunULp0+fRkpKCho2bAgA+PvvvwEAGo3hA5j5+fkA/v9h408++US1Gjdo0CCD08sbE7sKwN/fHytWrEDz5s2xYsUKxMbGYuPGjfjmm2+wYcMGxMXFoV+/fmjevDlatWpV6AHKiqJAUZRCJe17/19RFPTt21e1bPz999/j2LFjCAgIwKBBg+Dn54fp06fj7NmzCA4OhpubG8aNG4dRo0YVWi4+Ph5Tp07FRx99hM8++wwLFixAw4YNcfv2bXz44Yf48ssvcfHiRcTExKg+eP6nn34yadwA9GMv+EdW2rEvWbIEp0+fNtu4ASApKQnPPPMMGjRogD///FN17J9++qnV+gaAFStWYOnSpWjdujXWrVuHPn36YNmyZfjxxx9x9uxZvPnmm+jatStcXFzQvHlzAMAvv/yi36aGzou5/8PvtddeQ7Vq1YrMN2PGDCQkJFjlPQeAjRs3qo49KSkJx48ft8i4AeDrr78263a/t28RwdixY8utb1PGPmvWLKxatYr7Wznvb/Pnz7fa98ratWvN+tl+f5rSr18/1b5/+OEHHDlyBIGBgXjjjTdw6dIlLFmyBP/++y9CQ0MRFRWFxYsXw93dHQCwcuVKTJs2DTdv3sTrr7+OcePGYciQIfq4fvrpJyxbtgw6nQ6ZmZlwdnY2+N4oioLr168bjKnclXuNkIrQarVy7tw5ERGJiYmR0aNHi8jdR6Q4OztLnTp15K233iqy3KRJk8THx0eaNWsmq1ev1peOv//+e/Hx8ZHHHntMNBqNREVFSffu3Q3+1a9fX5YtWyYZGRni5eUlCQkJIiKyf/9+qV69ujg5Ocnx48eL9H3s2DFxcnKSOnXqyL59+4q079mzR3Q6XbFl+rKMu06dOrJu3bpSj93Jycli4w4KCpIOHTqojt2S27ykvkVEdDqd/o7sL7/8sowcOVJERI4ePSpubm7SuHFjiYmJkdzcXP0yJ0+elP79+0uDBg3kl19+kbp168q8efMkKSlJkpKS5O2335Zq1apJvXr1RKPRSOPGjSU0NLTIn6IoVnvPu3fvXuzYNRqNxcbdtGlTs2/3X3/9VZo3by5arbbc+zb1Pef+Vv77mzW/V8z92X78+HHp06ePuLu7l9i3l5eX/vMxNDRUlixZIiIiJ06cEAAGn7d98OBB0el00qlTJ/nuu++KtC9dulTat28viqJUikOxTOwqgMDAQFmzZo3k5uZKQECArFy5UkTuPvDZzc1NHBwcVL/oFUWRzZs3F2lLTEyURo0aSVBQkFy9elW17/j4eLG1tRU3NzcJCQnRP3T+k08+kQ4dOkj//v3lww8/LLLcRx99JH379hUHBwfZtWtXkfadO3eKg4ODxcbt4OAgDz/8cKnHbs1x39v3I488Yta+dTqdvPnmm9KxY0dp1KiRBAcHy9NPPy1fffWV/gvE19dXtm7dKiIiDRo0kB9//FFERI4cOSLOzs6i0+nkyJEjRdZ/5MgR0el00qJFC1m1alWR9lWrVkmzZs2K/fCz5nte0tgBWGzcJfVd1u0OwGp9lzR27m/W2d+s+b1iye1uZ2dXbN/PPvusNGvWTIYPHy6Ojo76eX/99VfRarX6H9L3SkhIkMaNG4uDg4McO3asSPvRo0fFwcFBkpOTJT8/X7XvisLWuvVCAoChQ4eiT58+8PX1haIoiIiIAADs3LkTjRo1gru7OzZv3ox69eoVWm7Lli1QFEV/ntW9XF1dkZycrH/cSmZmJn788UecOHECfn5+6NevH6pXr46XXnoJ4eHhOHv2LJ544gn9eQZ16tTBe++9h40bN+L999/Hpk2b9KXpHTt2YOvWrXj11Vfxzz//oGfPnujXrx8++ugjAMDSpUsxbdo0/Ti++eYbzJs3D2fPnkWtWrUwatQo9OvXr0zjfuyxx5CYmGhw7J999hlOnTqFrKws1W3+0ksvoWXLljh37pzZxr137168+OKLqFu3LgYNGoTo6Gj069cP33zzDaZNm4b8/Hz07NkTU6ZM0W/zzp07m63vQYMGITc3F3/++SccHBxw/PhxPPvss8jOzsZrr72GhQsXYvXq1ejZsyeeffZZ1K9fH9euXUNUVBSAu+ea1KtXDw4ODjh8+LD+MG+Bw4cPIyQkBElJSahdu3aRbVq7dm0cOnRIf7jIEEu95/fv72qKG7uTk5PFxl1S32Xd7jqdTvVcK0v3XdLYub9ZZ38rj+8VS/Rd0na3sbFB9erVVfuOj4/Hm2++iXPnzuF///ufft69e/eiT58+ePnll/H222/j0UcfBXD3s3XKlCn44IMPMHr0aMydOxfvvvtuofPovvrqKwQEBCA1NRUrVqxASkoKgLvnardq1QotW7YsdnuUN55jV0H8/PPPOHfuHHr37g1/f38AwOLFi+Hm5oZLly5h0qRJ6NOnT6Gd8aeffoKLiwuqVauGV155BW5ubujatStSU1NRt25dPPLII9i2bRvOnTuHxx57DDdv3kSDBg1w8uRJ2NraYseOHQY/OO5VUnteXh6uXr2K27dv60/Wz87ORosWLbBy5Ur89ttvePnllxETE4OHHnoIR48exVdffYU5c+Zg2LBhpR73O++8g08//RT29vZ45ZVX9CetpqamwsfHBwBQv359DB8+HIMHD9ZPM1Zpxp2bm4u6devi0qVLiIyMxNatWzFmzBh89NFHeOWVV6DRaPDxxx/j4Ycfhq2tLS5dugSNRoM6deqge/fuGDJkCGxsbErVd05ODtzc3DBixAh8+OGHAO4+Q3Lu3LnYsWMHbty4gU6dOqFdu3b4+OOPMWfOHJw7dw5DhgxB06ZNAdw9x8XZ2RnOzs54/fXX8X//93+Ftnt8fDymT5+Ot99+G/Xq1cPbb7+NsLAw/Xs+YsQIHDhwABMmTMD27dsLffi1bt0a3bp1g729vUXe80GDBuHOnTv466+/DG6z1NRUfPHFFxg/frzq2I8cOYLVq1eXatz//PMPfvvtN7i5uRU55yknJwfbt29Hq1atLLLdn3vuOezbtw87d+6Eh4cHrl69igULFiArKwu9e/fGQw89hJycHIu95//88w/27dtXaMx16tTBmjVrUL9+/XLrW0SwadMmnDhxAr6+voiMjISdnZ1F9rc+ffogJycH27ZtAwBs3ry50I/X2NhYtGrVqtixl3Z/a926Ne7cuYN//vnH4L5+L3N/rwwaNAgXLlxAz549ERkZiTZt2mDDhg34+OOP9T9cR44cWaa+i9vuAwYMwLlz59CmTZsin5+PP/54idvj3gskCs6TK0iBCs7pK3g9fPhwAMCuXbtw7Ngx1KtXD4cOHUJgYKD+NlKpqak4e/Ys2rRpg//973/F/rgqT0zsKpB7q2q+vr7o378/qlevrnq1DlD0pNK6devi3LlzyM7Oxo4dOxAeHo6BAwfi9OnT+OOPP+Dq6oqMjAw8/vjj8Pb2xm+//QZAvapmrOPHj+uvburbty+OHDmCWrVqoVmzZnjxxRcRExOjnzcyMhJHjx5FcnJysessbtz3j71u3boAoB/74sWLsW/fPixduhRpaWmIiopCTEwMoqOj8dlnn2HXrl3FVtRsbY0rZt877kaNGiE6OhoffvghevbsiaSkJISFhWHx4sUYMGAA9uzZg/bt2yMvLw8tWrTA9u3b9RW1NWvWIDg4GKtXr4azs3Op+g4NDcU///yDOnXqALh7orVOp8O5c+fg7e2NdevWYciQIbhw4UKx6y1puwP/f9sXfJgeOHAAeXl5cHBwwLVr1xAeHl7ow2/nzp3w9/fHn3/+WeRXuil9q73n9evXx4oVK1TXnZSUhGbNmiEvL6/Ufd/b/73jzs/Ph7e3Nw4fPgyNRoNnn30Wn332mT7BS01NhZ+fX7F9l7b/vXv3Ii0tDQDg5uaGdevWoXfv3rC1tUV+fj4uXryILVu2oFmzZmbv+8CBA7h9+zZGjhxZ5H6LcXFxeP311/U/qF5++WWz952WloY///wTnTp1wvXr1xEdHY1du3bB09MT165dQ4MGDZCYmKh69aQxfavtbzY2Npg9ezZGjhyJX3/9FT179sTTTz+Nhx56CMeOHcPKlSuxfPlyPP3002Yf95UrV2BjY4OOHTtixIgR6NGjh+oV8MUpzfdK9erVceXKFYSEhOD48eOIj4/HSy+9hL59+8LGxgZLlizBtGnTMHr06FL3fX//Bdv97NmzAAAXFxc4Ojri/PnziI6OxtWrV7Fnzx707NkT3333Hfbt21fkh2VBVU3th9+9Ll++jF9//RW3bt0CADz00EP4+++/cfPmTXz99ddFKqxHjx7FsGHD4Ofnh59++qnE9ZcHJnZWFBwcjC1btsDDwwPnzp1Du3btcOPGDX1Vzc7ODtu3by+xgiMiWLduHY4cOQLg7o4YGRmJlJQU1KhRA3Xr1sW8efPwxBNP6JepX78+0tLScPnyZXz11VfFVtVM5enpiTVr1iAsLAze3t5Yu3YtQkJC9O0ajQYiUqaKmrFjz8nJwS+//IKFCxdi/fr1cHR0RFZWFjp27IikpKQiFbVZs2bhxRdfxDvvvFOqeBwdHfVXZAGAvb09/v77bzz88MNo27YtWrRogfnz5yMzM9NgRS0sLAx+fn7YsmWLwYre5cuX8fnnnxtsf/fdd7F06VK0adMGAHDp0iXUrFkTmZmZcHBwQHJyMho1aoQlS5YUW1E7c+ZMieO8desWVqxYgStXrui3+/fffw8XFxcsWbKkyO0A0tPTMWjQINy+fRtr1qwpsr6CitqkSZNK7NvQe16jRo1i72J/5MgR9O/fH2fOnFGtqv3yyy8IDw83edybNm3CyZMnMXfuXNy8eRPjxo2DoihYu3Yt3N3d9VXkDRs2ICQkpEhVrU+fPmjUqFGptvvq1asRHh6OTz75BF988QXmzJmDJ598El9++SUAYNiwYbhx44b+assC91bUAJT6PX/++edRs2bNIj+Ezpw5Az8/P9jZ2UFRFJw6dUq1olaWvgv+nb/00kv466+/sHLlStSuXRu7d+9GTEwMWrVqhc8//1y1omYMQ/tbjx49cPDgQdSuXRuPPvooevTogTfeeEO/TI8ePXD8+PFiq2plGffChQuxYsUK/PHHH3BxccHAgQMxYsQING7cGMDdKz137dpVYlXN1HGPGTMGY8aMQUxMDDZu3Ijo6GjMmDEDL730EgBg0aJFmD59Ol5//XWDn1HGVNXU+v/222/RtGlTzJs3D4qi4IMPPsBff/2FP/74A8ePH8fjjz8ORVFw7tw5s1fVnJ2dkZiYqK+43m/v3r3o0KED/vvvP5PXbRHlcSIfGXbvya8DBgyQ1q1by82bN0VE5L///pOIiAjp37+/fv6MjAxZuHChTJgwQebOnSvXrl0z2Pbpp58Wukmtn5+fHDx4sFDfOp1OtFqtiNy96eX8+fMLtb/33ntSr149/eslS5ZI69atxd/fX9q0aSNTp06VU6dOGWz38vKSjh07isjdm/G++eabRcZdu3ZtGT16tHh6eoqdnZ107dpVfv/9d8nLy5NRo0ZJYmKi6nYrqV3tpOIzZ86Iu7u7eHl5iUajkf3794uNjY18++23+nmWL18unp6e8txzz8n333+vH9tDDz0kDRs2lPHjx8vs2bNV293c3PQnCh87dkw0Go3+hGkHBwf56quvJCgoSERE8vLyxM7OTlJSUkREZO7cuaIoioSFhUnbtm3FxsZGnnvuOenbt6/+QgtXV1fVdl9fXwkODpY///xTNmzYIB07dpQOHTrox7Zw4UKxtbUVnU4n7du3lz59+kifPn2kffv2otPppG7dugZPaDaWg4NDkf3sXgcOHFC9sGT//v2i0WhK3beiKKLRaERRlCJ/Go1GAAgA0Wg0+u3233//6ZdPSUkpdf9+fn6yc+dO/es7d+5Ily5dJDQ0VK5duyZ//PGHABBFUcTd3V327NkjtWvXlvr160vdunXFwcFB9uzZU6q+3d3d5dChQyIikp2dLRqNplAsr732mri6uha5obWNjY3+htZz5swpVd8iIs8//7yEhobqYyhga2srjz32mP7z7Nq1a9KyZUtRFEX/769Ro0Ym30j7Xvf+O2/YsKH8+uuv+raWLVvKe++9J7Vr15YVK1aIRqORrl27yhtvvCE9evQQOzs7+f3330vdt6urqyQlJYmISI0aNfT/f29sACQiIkJ++OEHycrKKnVf97t33KmpqfLBBx9Io0aNRKPRSIsWLWTAgAFia2srYWFh4uLiIt988404OzvLiBEj5PnnnxcHBweZPXt2qfp2cHCQM2fO6F/b2dkV+je/ceNGURRFatSoIQEBAaIoijz11FMSHh4uNjY20rt3b8nJySlV346OjoUubMjKyip0IUWrVq1Eq9WqXpDSunVr6dSpk8yePVvGjRsn48aNk9mzZxe6EG3nzp0G26tXry6bNm1SjW3jxo1SvXr1Uo3LEpjYWdG9/0Dr1Kkja9euLdQeFBQkfn5+IiJy9uxZqVWrlri6ukqLFi1Eo9GIp6ennDp1Ss6ePStBQUH6Ng8PDwEgDRs2lKZNm0q1atX0TzQo4OrqKl5eXiJy94Np//79hdobNWok9vb2IiLy5ZdfioODg7z88svy+eefy5gxY0Sj0cirr75qsH3EiBGiKIo0aNBA4uLixMHBQdq2bSsxMTHSrl07ASBLly4VkbtfRsuWLZPIyEixsbERPz8//Rdw/fr1Zfr06XLp0qUi262kdrWrxQqubCrY1nZ2doUuf3/11VcFgPTq1Ut8fHxk+vTpUr16dXnvvfdk6tSp4ujoKPb29sW2Ozo6yogRI6R27doybtw4CQwMlM8//1w8PDykRo0a8sorr4iIyMWLF0VRFLl165aIiDRv3lxsbW31sXzzzTcSHh4uIiLXr18XJycn/WtD7U2aNJH69euLra2tKIoirVu3LpR8N23aVJo3by5paWlFtktaWpp069ZNIiIiZNmyZTJmzBjp16+f9OvXT8aMGSM//vijZGVlSVZWlmq7r6+v6pdlUlKSzJkzRzw9PfW3brj3b9myZcUmVikpKfLOO++otnt4eEi3bt0kOTnZ4N/jjz8uAGT37t2ybt06CQsLk+bNm+vvhJ+SkiIACiV7BbKzs+Wvv/6Sc+fOGWx3cnLS788FcnJypHv37vLII4/Io48+KgAkPT1dPvroI/H395cRI0bo5x06dKhER0fLhg0b9D/Wrly5ItOnT5d33nlHDh06JFevXjXYbm9vL+vWrdOvq1q1anLy5En964IEIygoqNCfoihSs2ZNCQoKktq1axcZU+3atQ1eHWioffny5RIQEKB/moqI6PfBgn+HL774ogQHB+v3x3PnzklYWJi88MIL+mXy8/Nlw4YNMn/+fPn9998lOzu7UJ/3tq9cubLQj9caNWoU+nfs5OQkmzdvFq1WK+Hh4TJ9+vRC65oyZYo0adJE/zoxMVGeffZZadu2rQwYMEBWrFhR6GkO97e3bdtWxo0bJyIikZGRRZJjRVHE29tbunXrJnZ2dlK9enUZPXq0HDx4UD7++GNJTk5W3bYltat9viUmJsrgwYNFURT9Z/eGDRtEp9NJfHy8fr6vv/5a/P395a233pItW7aIyN2rQqOioiQyMlJiY2NV29zd3fU/qi9cuCCKohS6cjc8PFycnJz0V45Onz5doqKiROTuD92goCCZMGGCLFiwQIYOHSpPPvmkREdHy6hRo2T9+vUicvdpPIbaq1evLnv37tX3dePGDVEURf8UCicnJ7GzszO4zVJTUyU0NFQASK1ataRly5bSsmVLqVWrliiKIi1btpTw8HBRFMVgu6+vrwQEBMjy5csLfX6mpaXJ8uXLJSgoSEaNGqX6npU3JnZWVFJVTVEUfVXt/oqeoijSrl076d+/v8FqX506daRx48by9ttvy9tvvy2rV68utO7g4GCpU6eOiBiuqtnZ2UnDhg1FxHBFz97eXl/RM9Q+f/58qV69ugQHB4tOpxN7e3upVauWPPvss8VW1CZPnqyvbqhV9BRFkfXr16u2F3cpfu3ateXPP/8UkaIVNZG7tycoSHgNVfR8fHzEx8dHtf3nn38WDw8Pefrpp2Xq1KmSn58v33//vQQEBIhOpxM3Nzf55ZdfDFbUtFqtBAYG6l/fX9HTarVSo0YN1fa1a9eKn5+f3L5922ACUlJFreDLUq2iFxgYqB+HoXZ3d3dxdXWVmTNnSlJSkqSkpEhKSookJSXpE4yC9/b+ilrBf9WUVNFr1aqVKIqi2u7l5SX3HqC4t6r277//6j/0DVXzkpKSiq32PfTQQwb7LkjuCsYnYriqtnjxYv22MFTR02q1Uq1aNYPtdnZ2Ym9vr//CW7lypf6Hgojoq1OGKmr//vuv6iPqCip6PXv2lJ49e5ZY8Tt//rx06tRJnnzySbl06VKRxO7+ipqI6L84Re5W9Aq+WAsqek5OTnLixAnV9oKKWI8ePcTd3b3QjwpXV1f55ptvxNvb22BFLSQkRP/ZaqiipyiK/j5rhtptbW3F2dlZBg0aJO+++65Uq1ZNBg4cKO+//74MGjRIAOiTvfuragX7UseOHQ1W8xRFERsbG9VqX0mPtXJwcJBp06bpX99fVXv//fcFgMGK3mOPPSYAJCAgwGC1z8bGRry8vOS9996Tli1byuDBg6VRo0by559/yurVq0VRFOnVq5e+r/urap9//rnY2NioVvSioqKkVq1aBtsL3vuDBw/KqVOnpG/fvtK0aVN9Xy4uLvrP7vv16tVLHn74YXFzcyvSduTIEfHw8BAPDw/Vat+jjz4qderUEXt7e9FoNKLT6USn04lGoxF7e3t58cUXi300WnljYmdFiqJIkyZNVKtqiqLok4j7K3qKosjKlSslICDAYLVv69atEhAQoNr3hQsXJCgoSNq1a6daVSso1xuq6Lm7u+s/GA21nzhxQvWwW0n3X1IURZYtWyYi6hW9HTt2qLZPmDBB9ZDim2++KV5eXgYravPmzRNFUWT48OH6+e+v6BX8g1ZrT05OFkdHR4N9//fff9KnTx/Vipq3t7dMmTJF//r+il7NmjX1v8QNtZ8+fbpQbPcrrqIm8v+/7NQqel5eXuLl5VVsxa9+/fri6+urT9TuPTz6zDPPqFbUCg5DG6rmJSUlyYcfflhs+6uvvlpsYufk5FSkalOQeLm7u0uTJk1EURSD1bzevXsXW+2LjY0VtbNacnJyxMbGplD7/VW1gsPqahU9Pz8/qVWrlsH2t99+W9q3by/du3c32P+ECRMkPDzcYEXt33//FUVRxN/fX7WiV5DsGlPxy8/Pl6lTp4qPj4/Y2NgUW1ETkUKVJUMVPQAyePBg1fbq1atLw4YNZciQITJkyBD9Z4aISNeuXSU8PFwiIyMNVtS0Wq3+lAhDFT17e3sJDg5Wbf/0008lODhY+vXrJ87Ozvp93M7OTlq3bl1sVQ2A2Nvbi42NTZFqXsF2+frrrw1W+wrai/v89Pf3L7aqFhQUpE9w7q/oBQcHy5AhQ+Shhx4yWO2bN2+e/qbSI0eOlKysLPnoo4/E3t5e/37eW0G+v6rWvn170Wg0qhU9BwcHCQsLM9i+fft20Wq1+sS4Vq1ahW7UHhkZKR4eHgarajqdTnx9fVWrao6OjsXef3TPnj1SrVo1SUtLkw0bNsh3330n3333nWzYsMHg56G1MbGzooJqmlpVDYD+A/v+ip6iKLJ3717R6XQGq33JycnFfsmL3P1H98YbbxisqkVFRekTHEMVvdDQUHF3d1dtnzp1aqFDHfcq6eaWxlb01Npr1aqlWt3Jy8uT999/32BFrXr16lKtWjX55ZdfRMRwRa+gAqDWvmrVKv0Xhhq1itro0aOlcePGqufIde/eXezt7VXbV69eLXXr1lXt96233hJ3d3eDFbWZM2cKAHnxxRdVl7/3vExD7j2H7tSpU7Jt2zbZtm2bnDp1Sjp37izvvvuu6rIFFb2SzpFTay+p4tekSZMiP5xE7iZeOp1OvL299cvff46cj49Pof3t/vbz588Xm1Q2atSo0N3s76+qubi4iLe3t4gYrui5uLjo9zlD7Xv37pWaNWsa7DszM1Pu3LljsKL277//FnuOnDHthuzZs0dmz54tiqJIdHS0wYqaiOgrMCKGK3qKougr2Iba169fb/AwsojIoUOHxMPDQwYMGKBaUSv4EWWooldw42S19hMnTuh/wOXn50tKSopcvHhRf/i4uKqaoihy4sQJmT9/vsFz5BRF0Sf+htrnz5+vT5QMiY2Nlfr16xdbVevTp49+/nsreg4ODrJlyxb92O6v9p0+fdrgD9fbt29Lenq6DB48WNq3by+HDx82WFXT6XT6YoVI0YqeVqsttC/f375ixQqpWbOmHDx4sMi5enfu3JEXXnjBYFUNgHTr1k21qubi4iIuLi6q27SinUNXEiZ2FVhxFT1FUaROnTpiZ2dnsNr3119/qX7YG6Okip6dnZ14e3urttvb2xu8a7qx4za2omdIfn5+kQqmsUqq6Lm4uIiDg4Nqe0BAgP4cOlOVVNFbsWKFtGrVSrV9zZo1hZJMQ6ZPn26woubr6yvOzs7FVvQ8PDz0ybwhv/32m/j6+hpsW758ucyePVuGDh2quu5hw4apVvRcXFxEURTV9lWrVhWb2L3wwgv681Xv5+TkJJ06dSqUnN17jpyDg0ORxO3e9gMHDhTb9yuvvCLt27dXbbezs5PIyEj96/sregXndaq1nzlzRvVH3NmzZ/Xb/P6KWkFipnaOnLHtaoYMGSK9e/eWevXqFamoidz94VpwkVVJFT1D7cnJycX+0EhMTJTatWsbrKi1aNGi2HPkQkJC9Pu6ofYvv/xS6tevr9q3oijSr18/1bbizpEDYDB5Kmh3cnISJycn1b6PHDkiDRo0UK2qabVafZJ8f0XP399fPvnkE/H39zdY7du0aZP4+/ur9r1v3z7x8vLSf77cX1Vzd3eXsWPH6l/fX9Hz9vYudJ7c/e2nTp1Sfc8L9nVDVbURI0ZIrVq1VM+Rc3Z2Fmdn52LPoXvhhRdk8+bNBvf727dvy+LFi1W3S3ljYleBFVfRe/vtt6V169b68+jur/a99tprqh8sxiquord79+4S20urpIpeSe1lUVJFb/DgwTJ58mTV9iFDhkhGRkaZYlCr6Bnbboz7K2oiJVf0Cn4Bq7V7eHjI5MmTVfss7jy5kip6rVq1Uj3cWbDu4qpme/bsUW1v0qSJLFu2rMgJ6wXJm52dXbHn0AUGBpbp/MCGDRvqz/sUKVrRq1Wrlnh6eqq279ixQ/XL1lDfBRW1gkPJIqJa0TO2vTRjVxRFnnzyyWIrevb29qrtO3bs0Fc6i+vbUEXt0KFDUr16ddVz5Ap+MKu1a7Va+frrr0s17pLOkdNoNPLxxx+rtqelpRU5p9mYvguqasVV9Lp06SL29vbSrFkzg9W+Jk2ayLBhw0rs+9ixYwaraiVV9CIjI0Wr1aq2b9q0SfUUo+K2eXHVPHt7exk5cqTExMSotvfv318CAwP1CWu7du3kwoUL+vWX5ap6S+B97IgIwN0br/bo0QMXL15ESkpKoTuz+/j4YMyYMRARzJkzx2B7586d0bNnT9X1nzp1Cq+++qrBG/V++eWX+Pbbb1VvIPrNN99gy5Yt+OKLLwy2f/fdd9i1axc6depkct9vvPEGduzYgbp162LhwoWF2nJzc/HQQw/hxIkTRW7aWtDepk0b7Nq1C7/++qvJfQN3b+a7b98+bNq0yWB7u3btcOfOHezatatI22+//YZvvvkG58+fx/jx403u+9y5c5g8eTIWLlwIEcH06dPxySef4MqVKzhw4ACCg4P18xpqP3HihMH1GtN/nz598Pfff6Nt27YAgKioKPTp00ff3qRJE1y/fh2dO3c22N6zZ08kJyfj7bffNrnvc+fO4ZVXXoGdnR1WrVqFjIwMAICtrS1atGiBsWPHokmTJnjzzTcNtrdv317/xART+9ZoNNi7dy8+/fTTIvtbQXvB/fkMKbipfGn6Bu7eULdr166wt7dH69at8emnn+KTTz7BxIkTkZ2dDV9fX7i6uqJdu3aF2nJycvDwww9j7NixBh/3ZUzff//9NyIjI3H16lUoioKAgAD88ssv+vvDLViwAO+99x7OnDlTpP23337D1q1bcePGDYM3fj516hTi4uKQmJgIDw+PQvsuANy5cweLFy9GgwYNCt3DMywsTH/fzfT0dOzdu7dI++DBg5GTk4NFixbh5s2bGDNmDA4dOoRNmzYhMDDQ6JuQlxvr5ZREVJHc+4vXUEXvXoba7z3fTe1P7VetJe9jV1LfOTk5smXLFtX227dvF3s1cVn6Fil57AXnyZVX34YqemrtlnzPMzIy5Pbt26rt5urbUEXvXobaLf2eF8dSfRdU9AwpaDNX32oVvQKG2kvqG/ech2vuqlqNGjXkwIED+tf5+fnywgsvSGBgoJw8ebLCVeyMe24SEVV6xvzSL1C7du0iTzy5t7pjqL1GjRp4+OGHkZCQYLDvgl/UhuK4t29D7u3bkOL6BoDZs2er9l1S/1euXMHMmTOt0jcAXLt2TXXsvr6+mDJlCrZu3WqwvTR9h4WF6Z9Nami739tuyff8+vXrpX7PTelbURT9UwoK3DtuQ+2WfM+tta/rdDr9Iwjv77+gzVx9Fzzx5H7Fjb2kfb1Tp07YuHEjLl++rK+qtW3bVl9Vu337NkQEhw4dMljN+/bbb9GoUSOD1b7bt29j1apVaNKkCYC7+8znn3+OUaNGoX379vjuu+8MjsdqrJ1ZElH5sHSVoX379qrnwRWsGypXtlqy74L+i+u7LJUlS/ZdUv9dunSRkSNHFnsem6X6Fqm47zn3N8u85xV5Xy+4MX+Be6tq69evF39/f/2tUu6v5m3durXYal9oaKgoiuFzdGNjY8XNzY0VOyIqf76+vvjss8/QrVs3g+0l/dr+7bffICKq7U2aNEFiYqJq3zNnzoSPjw/at29frn0Dd8+ViY2NxcSJEw2uu7jqjiX7Bso29latWiE5OdkqfQMV9z3n/la693znzp1ITU1V7b8i7+s5OTlYvHix/vW9VbWnn34aTZs2xYULF3D06NEi1bz33nsPAFSrfdHR0di/f7/BfufOnYv8/HzMmzdPNbZyZ62MkojKV5cuXfR30zekpF/bKOZeciX9Grdm3yX1X1J1x5J9W3rs7Lv8+y6p/4q8vxX0XVz/FXVfb9GihSxZssRgm4ODgzg7Oxc6r/Lec+Q8PT1FUf5/Ra6in0NXEiZ2RA+IxMTEQrfVuJ+vr6+89957qu0F96dS8/fff6t++Fmz75L69/Pzkx9++EH1Id+W7Fuk6m73B7XvkvqvyPubn5+fzJw5U3X9FXlfnzp1qv4pFfdzdnaWfv36FVl3bGys+Pv7i6Ojo8F+C9oTExMrVWLHQ7FED4jHHnus2PbmzZsjKytLtT04OFj1diTA3UMfonL3JGv2XVL/YWFh+Pfff9G3b99y7xuoutv9Qe27pP4r8v4WFhaGI0eOqK6/Iu/r48ePN3jLHwBo1KgRoqOj8f333xeaPnfuXABQPYxa0N61a1fVmCoiJnZEBAAYO3YsMjMzVdsnTJiA6Oho1fZ69eph48aN7LsS9c++H6y+S+p/7NixuHbtGp599tly7xuw3Nh79OiB77//Hs8991yRtrlz52LPnj3YuXOnwWUr5Dl0JeANiomIiIiqCI21AyAiIiIi82BiR0RERFRFMLEjIiIiqiKY2BERERFVEUzsiIiIiKoIJnZEREREVQQTOyIiIqIqgokdERERURXx/wD0OWgEgns3sgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "energies_diff = np.abs(np.array(energies_with_nist) - np.array(energies_without_nist))\n", "\n", "fig, ax = plt.subplots()\n", "ax.bar(np.arange(len(energies_diff)), energies_diff, tick_label=labels, color=\"C0\")\n", "ax.set_ylabel(\"Energy difference (Hartree)\")\n", "ax.set_xticklabels(labels, rotation=90)\n", "\n", "ax.set_yscale(\"log\")\n", "\n", "fig.tight_layout()\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.13.1" } }, "nbformat": 4, "nbformat_minor": 2 }