finesse-simulation-O4/generate_from_finesse-virgo.ipynb

179 lines
99 KiB
Text

{
"cells": [
{
"cell_type": "code",
"execution_count": 8,
"id": "8382eebe-ea67-402e-b4da-64d9f616f0c1",
"metadata": {},
"outputs": [],
"source": [
"from finesse.virgo import Virgo\n",
"from rich.console import Console\n",
"from matplotlib.pyplot import figure, show\n",
"from numpy import save, array, geomspace\n",
"from finesse.analysis.actions import FrequencyResponse"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "77a9edfb-d910-4264-a29c-1674a4841c8d",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Adjusting recycling cavity lengths...\n",
"Pretuning...\n",
"Optimizing demodulation phases...\n",
"Optimizing lock gains...\n",
"Running RF locks...\n",
"Switching to DARM DC lock...\n",
"Done.\n"
]
}
],
"source": [
"console = Console()\n",
"\n",
"virgo = Virgo()\n",
"virgo.make()"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "2d780140-9125-4a81-927e-daf298504add",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGhCAYAAACZCkVQAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAATGRJREFUeJzt3Xl8VNXdP/DPnSXLZCUJ2UhCkE22JJJNVBSVxbqLChhlkaqte0trW/v0V7VP+7TP8/RpaTXWilaRGsUNl4rKoiyiEAIEkCVhCUnIvpBMkskymbm/PwJUTcxsd+Zun/frxUvJOeee7wxfki93zj1HEEVRBBEREZHKGeQOgIiIiEgKLGqIiIhIE1jUEBERkSawqCEiIiJNYFFDREREmsCihoiIiDSBRQ0RERFpgknuAALF6XSitrYWEREREARB7nCIiIjIDaIooqOjA8nJyTAYhr8Xo5uipra2FqmpqXKHQURERF6orq5GSkrKsH10U9REREQAGHhTIiMjB7Xb7XZs2LABc+fOhdls9ujano51t7+rfr60+/J65RLImKWaS4155aqP1vIKCFzcUs7j7bWUmleu2tWYW3rKK2/Gutu/tbUVY8aMOf9zfDi6KWrOfeQUGRn5nUWNxWJBZGSkV3+Qnox1t7+rfr60+/J65RLImKWaS4155aqP1vIKCFzcUs7j7bWUmleu2tWYW3rKK2/GepJbANxaOsKFwkRERKQJLGqIiIhIE1jUEBERkSawqCEiIiJNYFFDREREmsCihoiIiDRBN490n2O3288/Hvbtr3/9v55e05Ox7vZ31c+Xdl9er1wCGbNUc6kxr1z10VpeAYGLW8p5vL2WUvPKVbsac0tPeeXNWE9zyx2CKIqi271VqLCwEIWFhXA4HCgvL0dRUREsFovcYREREZEbbDYbCgoK0N7ePuQ+c1+n+aLmHKvViqioKDQ3N3/n5nsbN27EnDlzvNpwyJOx7vZ31c+Xdl9er1wCGbNUc6kxr1z10VpeAYGLW8p5vL2WUvPKVbsac0tPeeXNWHf7t7S0ICkpya2iRncfP5nN5mHfPFftvlzb2/6+xjxcuy+vVy6BjFmqudSYV676aC2vgMDFLeU83l5LqXnlql2NuaWnvPJmrDv54C4uFCYiIiJNYFFDREREmsCiRgIVzV2o6ZI7CiIiIn1jUeOj14qrMPcvO/BeJd9KIiIiOfEnsY8uGxcHAChrN6D6jE3maIiIiPSLRY2PUmMsuHRsLADgrT21MkdDRESkXyxqJLAgexQA4O29Neh3OGWOhoiISJ9Y1Ejg6knxCDOJaOjoxZayJrnDISIi0iUWNRIINhmQN3JgY+bXd1fLHA0REZE+saiRyIyEgY+dPitrRIO1R+ZoiIiI9IdFjUQSQoGc0dFwOEW8WcK7NURERIHGokZCC7JTAABrS6rhdOrinFAiIiLFYFEjoWumJCAixITq1m58caJF7nCIiIh0hUWNhEKDjLg5a+Dx7td2V8kcDRERkb6wqJHYorxUAMCGQ/Vo6eyVORoiIiL9YFEjsSnJUchIiYLdIWLdvhq5wyEiItINFjV+sCg3DcDAYZeiyAXDREREgcCixg9uyExCqNmIE01dKKk8I3c4REREusCixg8iQsy4ITMJAPB6MfesISIiCgQWNX6yKG/gI6gPD9aivdsuczRERETax6LGTy5KjcaEhHD02J14v5QLhomIiPyNRY2fCILwtQXD1VwwTERE5Gcsavxo/vRRCDIZcLjOiq9qrHKHQ0REpGksavwo2hKEa6YkAuAOw0RERP7GosbPzu0w/H5pLbp6+2WOhoiISLtY1PjZjAtikR5rQWdvPz48WCd3OERERJrFosbPBEHAwrMLhl8v5kdQRERE/sKiJgBuzR4Fk0HA3qo2lDd0yB0OERGRJrGoCYD4iBBcPSkewMB5UERERCQ9k9wBBJrdbofdPniH33NfG6rNnWu6Gnv79GR8cqgB6/bW4JEr0tyay9V1fWn35fXKJZAxSzWXv/PKm/7u9PM2d9SYV0Dg4pZyHm+vpdS8ctWuxtzSU155M9bT3HKHIGp8V7jCwkIUFhbC4XCgvLwcRUVFsFgsAY/DKQJP7TWirU/AkvEOZMdp+m0nIiKShM1mQ0FBAdrb2xEZGTlsX80XNedYrVZERUWhubl5yDfFbrdj48aNmDNnDsxms0fXdnfsXz89jqc/O4n89GgUJDW77O/qur60+/J65RLImKWaKxB55Wl/d/p5mztqzCsgcHFLOY+311JqXrlqV2Nu6SmvvBnrbv+WlhYkJSW5VdTo7uMns9k87Jvnqt2Xay/KT8czW05i16k2zIl2fy5fYx6u3ZfXK5dAxizVXP7MK2/7u9PP29xRY14BgYtbynm8vZZS88pVuxpzS0955c1Yd/LBXVwoHECjokNx+fiRAICdjXzriYiIpMSfrAF2x9kdhnc1CbA7nDJHQ0REpB0sagLs6kkJiAsPQoddwJayZrnDISIi0gwWNQFmNhow/6JkAMDaPadljoaIiEg7WNTI4PbsUQCA7ceaUdvWLXM0RERE2sCiRgbpsWEYF+mEUwTeKKmWOxwiIiJNYFEjkxnxA9sDvVlyGg6nLrYKIiIi8isWNTLJjBURFWpCTVs3th9rkjscIiIi1WNRIxOzAbgpc2DB8OvF/AiKiIjIVyxqZLTg7ILhTUca0NTRK3M0RERE6saiRkYTEyOQlRqNfqeIt/fy8W4iIiJfsKiR2bkdhtfuroZOzhYlIiLyCxY1Mrs+IxlhQUZUNHdh58lWucMhIiJSLRY1MgsLNuHGrIG1NWt3V8kcDRERkXqxqFGARbkDH0Gt/6oebbY+maMhIiJSJxY1CpCREoVJSZHo63di3b4aucMhIiJSJRY1CiAIwvkFw68Xc8EwERGRN1jUKMRNmaMQbDKgrKEDpdVtcodDRESkOixqFCLKYsZ105IAcIdhIiIib7CoUZBFeWkAgA8O1KKzt1/maIiIiNSFRY2C5KaPwAUjw2Drc+D90lq5wyEiIlIVFjUKIggC7sgduFvzWjH3rCEiIvIEixqFuTU7BUFGAw7WtONQrVXucIiIiFSDRY3CxIQFYd7URADA2hIecklEROQuFjUKdG7PmvcP1KHXIXMwREREKsGiRoFmXBCL9FgLunod2NssyB0OERGRKrCoUSBBEM4/3v1lI/+IiIiI3MGfmAp1W3YKzEYBlZ0CjtR1yB0OERGR4rGoUai48GDMvjAeAPDGHi4YJiIicoVFjYItyEkBALy3vw7dfVwxTERENBwWNQp2yQUxiA0W0dHTjw8P1skdDhERkaKxqFEwg0HAjAQnAO4wTERE5AqLGoXLHynCZBCwp/IMyuq5YJiIiOi7sKhRuMgg4MqJIwHwbg0REdFwWNSowKLcgQXD6/bVoMfOBcNERERDYVGjApeOjcWo6FC0d9vx0VdcMExERDQUFjUqYDQIWJQ7cB7Ua7uqZY6GiIhImVjUqMTtOakwCEDxqVYcb+SCYSIiom9jUaMSiVEhuOrCBADA68W8W0NERPRtLGpUpCB/4COot/ee5oJhIiKib2FRoyJXTIhHUlQIztjs+ORQvdzhEBERKQqLGhUxGgQsyDm7YJh71hAREX0DixqVWZA7sGB458lWVDR3yR0OERGRYrCoUZlR0aGYNTEeAPDGnhqZoyEiIlIOFjUqdEdeGgDg7b016HfKHAwREZFCsKhRoSsnjkRCZDDO2Ow42CrIHQ4REZEisKhRIZPRcH7B8BeNLGqIiIgAFjWqtSAnFYIAlLcbUNlqkzscIiIi2SmuqKmursasWbMwefJkZGRk4M0333SrTW9SYyyYOS4WAPBmCRcMExERKa6oMZlMWLlyJQ4fPowNGzbgRz/6Ebq6uly26dGC7BQAwFt7a9DHFcNERKRziitqkpKSkJWVBQBITExEXFwcWltbXbbp0VUXjkSEWURLVx82H2mQOxwiIiJZeVzUbNu2DTfccAOSk5MhCALefffdQX0KCwuRnp6OkJAQ5Ofno7i42Kvg9uzZA4fDgdTUVI/a9MJsNCA/XgQAvLabh1wSEZG+mTwd0NXVhczMTCxfvhzz588f1L527VqsWLECzz33HPLz87Fy5UrMmzcPZWVliI8f2DQuKysL/f39g8Zu2LABycnJAIDW1lYsWbIEq1atGtRvuLZzent70dvbe/73VqsVAGC322G32wf1P/e1odpc8XSsu/1d9bPb7ZgR78SmGgO2H2vCycZ2pI6wuDXel9crl0DGLNVcaswrV320lldA4OKWch5vr6XUvHLVrsbc0lNeeTPW09xyhyCKouh2728PFgSsW7cON9988/mv5efnIzc3F8888wwAwOl0IjU1FQ8//DB+8YtfuHXd3t5ezJkzB/feey8WL17sdtvXPfnkk3jqqacGfb2oqAgWi2WIEepVeNiA8nYD5o5y4ro0rq0hIiLtsNlsKCgoQHt7OyIjI4ftK2lR09fXB4vFgrfeeusbhc7SpUvR1taG9957z+U1RVFEQUEBJk6ciCeffNLttm8b6k5Namoqmpubh3xT7HY7Nm7ciDlz5sBsNruM05ex7vZ31e9cu2NUJla8dQjxEcHY+pOZMBkNLsf78nrlEsiYpZpLjXnlqo/W8goIXNxSzuPttZSaV67a1Zhbesorb8a627+lpQVJSUluFTUef/w0nObmZjgcDiQkJHzj6wkJCTh69Khb19ixYwfWrl2LjIyM8+t11qxZg2nTpg3b9m3BwcEIDg4e9HWz2Tzsm+eqfTiejnW3v6t+86YkIfajY2js6MX2E2cwd0qi2+N9eb1yCWTMUs2lxrxy1UdreQUELm4p5/H2WkrNK1ftaswtPeWVN2PdyQd3SVrUSOGyyy6D0zn0RyjDtelZkMmA23JS8PetJ/FacdWgooaIiEgPJH2kOy4uDkajEQ0N33y8uKGhAYmJ/EHrT4tyBw653FLehJq2bpmjISIiCjxJi5qgoCBkZ2dj8+bN57/mdDqxefNmzJgxQ8qp6FvGxIVhxgWxEEVgLR/vJiIiHfK4qOns7ERpaSlKS0sBABUVFSgtLUVVVRUAYMWKFVi1ahVWr16NI0eO4P7770dXVxfuvvtuSQOnwe7IH7hb82ZJNfod/JiOiIj0xeM1NSUlJbjyyivP/37FihUABp5wevnll7Fw4UI0NTXh17/+Nerr65GVlYWPP/540OJhkt68KQkYYTGjrr0HW8ubcPm4GLlDIiIiChiPi5pZs2bB1VPgDz30EB566CGvg/InrW2+9/X/GgDckpWMf3xRiVd3VmLG6IjvHM+NrAIzlxrzylUfrW2QBuhrkzSl5pWrdjXmlp7yypuxitt8Tw0KCwtRWFgIh8OB8vJyTW6+93UN3cB/lZogQMST0x2IHvxUOxERkWoEbPM9NbFarYiKitLk5nvfbi94cTd2nzqDh2eNwbjeY9zISsa51JhXrvpobYM0QF+bpCk1r1y1qzG39JRX3oxV/OZ7aqDFzfe+3X5n/mjsPnUGb++rw2OTuJGVEuZSY1656qO1vAL0tUmaUvPKVbsac0tPeeXNWCk335P0kW5ShmumJiIq1Iza9h4cbRPkDoeIiCggWNRoUIjZiPnTRwEAvmhgUUNERPrAokaj7sgb2LPm0BkBjR29LnoTERGpH4sajZqQEIHstGg4IeDtvTVyh0NEROR3LGo0bEHOwEdQb5SchtOpi4fciIhIx3T39JOWN9/7ttkTYxFqFHG6rQdbyuoxc1yc1zErATffk6Y/N98bTE+bpCk1r1y1qzG39JRX3ozl5nte0Nvme9/2VoUB2+sNyIxxYvlEngdFRETqws33hqCnzfe+3r763Y347wMmmAwCtj92OeLCg72KWQm4+Z40/bn53mB62iRNqXnlql2NuaWnvPJmLDffk4AeNt/7uuQwIDMlCvtPt+Pd/Q24f9ZYn2JWAm6+J01/br43mJ42SVNqXrlqV2Nu6SmvvBnLzffIIwtzUgAAr++u4oJhIiLSLBY1OnDdtASEB5tQ2WLDlydb5A6HiIjIL1jU6IAlyISbspIBAK8VV8kcDRERkX+wqNGJczsMf3KoHi2d3GGYiIi0h0WNTkwdFYWMlCjYHSLe3nta7nCIiIgkx6JGR87drXmtuBo6eZKfiIh0hEWNjtyQmYywICMqmrtQfOqM3OEQERFJSnf71OjpmIRvtwcbgOszkrC25DSKdlVhXiS3HPf3XGrMK1d9tLaVPaCv7eyVmleu2tWYW3rKK2/G8pgEL+j9mIRvq+4E/njQBKMg4jfZDoSraw8rIiLSGR6TMAS9HpMwVPv853biYI0VN6Y58Idls1WzOyePSZCmP49JGExP29krNa9ctasxt/SUV96M5TEJEtDbMQlDtRfkj8bj7xzEF40GGI0m1XyDOIfHJEjTn8ckDKan7eyVmleu2tWYW3rKK2/G8pgE8smNmckICzaiuUfAzopWucMhIiKSBIsaHQoLNuHGjCQAwOu7uWcNERFpA4sanVqUO3DI5cYjjWjq4A7DRESkfixqdGpyUiRGh4vod4p4aw/v1hARkfqxqNGxSxKcAAYOuXQ6dfEQHBERaRiLGh27KFZEeLAJVa027DjRLHc4REREPmFRo2PBRuDmrIEFw0W7qmSOhoiIyDcsanRuYc7ZBcOHG9DY0SNzNERERN5jUaNzFyZG4KK0aPQ7RbxZwgXDRESkXrrbUVjPB1p+V9vC7FHYV9WG14qrcM8laTAYBDdeReDxQEtp+vNAy8H0dPCgUvPKVbsac0tPeeXNWB5o6QUeaOlanwP49R4juh0CfjjJgUnRmk4JIiJSER5oOQQeaDn84XC/+fAo1uyswpxJ8Xi2IMv1i5YBD7SUpj8PtBxMTwcPKjWvXLWrMbf0lFfejOWBlhLggZZDty2ekY41O6vwaVkTWrsdSIgMcf0CZMIDLaXpzwMtB9PTwYNKzStX7WrMLT3llTdjeaAlSW5CQgRyRo+Awynijd3VcodDRETkMRY1dF5BfhoA4PXd1XBwh2EiIlIZFjV03rXTkhAVakZNWze2HWuSOxwiIiKPsKih80LMRsyfPgoAdxgmIiL1YVFD33Dn2Y+gPj3aiPp27jBMRETqwaKGvmFcfATy0mPgcIpYywXDRESkIixqaJBzC4bX7q7igmEiIlINFjU0yDVTEzHCYkZtew+2lDXKHQ4REZFbWNTQICFmI26dPnB6NxcMExGRWrCooSHdcfYjqM/KGlHb1i1zNERERK6xqKEhjR0ZjosviIFTBBcMExGRKuju7Ce73S75sfb+Pm79u/r50u5ODAuyR2HnyVa8vrsKP5w5GiajvDWwL39Gcs2lxrxy1cfXvFKiQMUt5TzeXkupeeWqXY25pae8AoCunl6PxnqaW+7Q/CndhYWFKCwshMPhQHl5OYqKimCxWOQOSxX6ncCv9xjR1S/gnokOTIvRdKoQEZGX2nqBvx4y4vo0J6bHSfuzwmazoaCgwK1TujVf1JxjtVoRFRWF5ubmId8UJR637qqfL+3uxvCHj8vw4o5KXDEhDi8snu7ytfmTL39Gcs2lxrxy1UeKvFKaQMUt5TzeXkupeeWqXY25pZe86rU7UPDibhyosSLZIuKjH8+CJSRYsrlaWlqQlJTkVlGju4+ffDny3tdre9vf15iHa3c19q4ZY/DijkpsO9aMhk47UkbIf5fLlz8jueZSY1656uNLXilVoOKWch5vr6XUvHLVrsbc0nJeiaKIn687hAM1VkSHmnHPxG5YQoIlzS1PrsWFwjSsMXFhuGRsLEQuGCYiom95accpvLO3BkaDgL8szEBsiLzxsKghl/69w3A1+h1OmaMhIiIl2HG8Gb9bfwQA8MtrJ+GSsbEyR8Sihtwwd3Ii4sKD0NjRi81HucMwEZHeVbfa8GDRXjicIuZfNArLL02XOyQALGrIDUEmA27LTgXAHYaJiPTO1tePe18pQZvNjoyUKPzX/GkQBEHusACwqCE3LcodKGq2HWtCdatN5miIiEgOoijisTcP4Gh9B+LCg/D3xdkIMRvlDus8FjXklvS4MFw2Lg6iCLy+m3driIj06NktJ/DhwTqYjQL+dlc2kqJC5Q7pG1jUkNvOLRh+o+Q07FwwTESkK58ebcAfN5QBAJ68cQpy02NkjmgwFjXktjmTExAXHoymjl5sOtwgdzhERBQgJ5o68ehrpRDFgX/g3pk/Wu6QhsSihtxmNhqwICcFAFBUzI+giIj0oKPHjvteKUFHbz9yRo/AkzdMkTuk78SihjxyR14aBAHYfqwZVS1cMExEpGVOp4gfry3FiaYuJEaG4Nm7piPIpNzSQbmRkSKlxlgwc/xIAMBrXDBMRKRpf95Ujk1HGhFkMuDvi7MRHyHzlsEusKghjxXkDSwYfrOkGn39XDBMRKRFHx2sw9OfHgcA/P6WachMjZY3IDewqCGPXT0pHiMjgtHc2YeNXDBMRKQ5ZfUd+Mmb+wEAyy8dg1uzU2SOyD0sashjZqMBC3PO7jBcXClzNEREJKUuO3B/USlsfQ5cOi4Wv7z2QrlDcptJ7gACzW63w263D/n1r//X02t6Mtbd/q76+dLuy+sFgNumJ6Fwy3HsON6C4/XtGB1r8eo6nvA1ZjnmUmNeuerjz7ySS6DilnIeb6+l1Lxy1a7G3FJjXnX39uLlYwZUt3cjJToEf759GkSnA3anQ/I4PM0tdwiiKIpu91ahwsJCFBYWwuFwoLy8HEVFRbBY/P8DWA+eO2LAkTYDrk524sbRXFtDRKR2604ZsKXOgCCDiB9NdWBUmNwRATabDQUFBWhvb0dkZOSwfTVf1JxjtVoRFRWF5ubmId8Uu92OjRs3Ys6cOTCbzR5d29Ox7vZ31c+Xdl9e7zmbjjTi/qJSxISZsf2nV/j9MT8pYg70XGrMK1d9/J1XcghU3FLO4+21lJpXrtrVmFtqy6t39tXg5+8cAgCsvH0KrssY5dc43O3f0tKCpKQkt4oa3X38ZDabh33zXLX7cm1v+/sa83DtvrzeOVOSkBh5FPXWHnxa3oIbMpO9uo6nfIlZrrnUmFeu+vgrr+QUqLilnMfbayk1r1y1qzG31JBX+6rO4P+9fwQAMG+UE9dljFLM9yxPrsWFwuQ1k9GABWdP7y7axT1riIjUqMHagx+s2YO+fidmXzgS16SqdzkBixryyaLcVBgE4MuTLTje2Cl3OERE5IEeuwM/WLMHjR29mJAQjv+9bRoMgtxReY9FDfkkOToUV12YAIB3a4iI1EQURfzHuq9QWt2GqFAzVi3JQXiwulelsKghn9158cAOw2/tqUaP3fVjf0REJL+XdpzC23tPwyAAzxRchNGxCnjUyUcsashnV4wfiZQRobD29OOD/bVyh0NERC58fqwZv1s/sDD4l9dOOn+mn9qxqCGfGQwCCvIH7ta8yo+giIgUrbKlCw8W7YXDKeLW6Sn4/mVj5A5JMixqSBILclJhNgoorW7DVzXtcodDRERD6Oztx72vlKC9247M1Gj87papEAQVrwz+FhY1JIm48GBcMzUJAO/WEBEpkdMpYsXaUpQ3dCI+IhjPL85GiNkod1iSYlFDkrnr7EdQ75XWoKNHPeezEBHpwV82H8OGww0IMhrw3OJsJESGyB2S5FjUkGTyxsRgfHw4bH0OvLuvRu5wiIjorI+/qsNfNh8DAPzulqmYnjZC5oj8g0UNSUYQBNx59m7NP3dWQSfHihERKdrReitWvLEfALD80jG4PSdV5oj8h0UNSeqW6SkINRtR1tCBksozcodDRKRrrV19uPeVEtj6HLh0XCx+ee2FcofkVyxqSFJRoWbcePZgy1d3VsocDRGRftkdTjz46l5Ut3YjLcaCZ+6YDpNR2z/2tf3qSBbndhhef7AeLZ29MkdDRKRPv/vwCL482YKwICNWLcnBiLAguUPyOxY1JLmMlGhkpEShz+HEW3tOyx0OEZHuvLG7Gi9/cQoA8KeFWZiYGCFvQAHCoob84tyC4aLiKjidXDBMRBQoeyrP4FfvfgUA+PHsCZg3JVHmiAKHRQ35xQ2ZyYgIMaGyxYbPjzfLHQ4RkS7Ut/fgh//cgz6HE9dMScTDV42TO6SAYlFDfmEJMuHW6SkAgH9ywTARkd/12B34wZoSNHX0YmJCBP5vQSYMBu0cgeAOFjXkN+c+gtp8tBF17d0yR0NEpF2iCPzqvcPYf7od0RYzVi3JQViwSe6wAk53r9hut8NuH7yF/7mvDdXmzjU9Getuf1f9fGn35fW6Kz0mBLnpI7D71BkU7TyFR3y8DRqImKWeS4155aqP3HnlD4GKW8p5vL2WUvPKVbsacyuQebWlTsB7lXUwGgT8dWEGkiLNAfm+48tYT3PLHYKo8W1fCwsLUVhYCIfDgfLychQVFcFiscgdlm7sbRaw+pgRUWYRT2Q7YNTXnVAiIr872ibguSMGiBBwa7oDlydp68e6zWZDQUEB2tvbERkZOWxfzRc151itVkRFRaG5uXnIN8Vut2Pjxo2YM2cOzGazR9f2dKy7/V3186Xdl9frid5+Jy7/41a0dtlReEcm5k5O8PpagYpZyrnUmFeu+ighr6QWqLilnMfbayk1r1y1qzG3AhFzZYsN85/bCWtPP+ZnJeEP86dCELz/16MSv2e1tLQgKSnJraJGdx8/mc3mYd88V+2+XNvb/r7GPFy7L6/XHWYzsDA3DX/bcgKvl9TguswUCa7p35j9MZca88pVHznzyl8CFbeU83h7LaXmlat2NeaWv2Lu6LHjh0WlsPb0Iz1cxG9umoKgIGk22FPS9yxPrsWFwuR3BXlpEARg+7FmnGrukjscIiLVczpF/Hjtfhxv7ERCZDCWT3Qg2MQf6XwHyO9SYyy4YsJIAMBrxVUyR0NEpH5/3FCGTUcaEGQy4Nk7shCl/RMQ3MKihgLirvzRAIA3SqrRY3fIHA0RkXq9V1qDZ7ecAAD8z60ZyEiJkjki5WBRQwFx5YXxSI4KwRmbHR9/VS93OEREqrS/ug0/e+sAAOCHV4zFzReNkjkiZWFRQwFhNAi4I29gMz7uMExE5LkGaw/ufaUEvf1OXH1hPB6bN1HukBSHRQ0FzMLcVJgMAkoqz+BovVXucIiIVKPH7sB9r5SgsaMXExLCsXJRFow6OwLBHSxqKGDiI0Mwd8rAPjWv7uSCYSIid4iiiJ+/feD8EQgvLMlFRIi6HmsPFBY1FFB3nl0wvG5fDbp6+2WOhohI+f629QTeK62FySDg2TunIy2Wu+J/FxY1FFCXjI3FBXFh6Oztx3ultXKHQ0SkaJsON+B/PykDADxx4xRcMjZO5oiUjUUNBZQgCCjI//eCYZ2c0kFE5LGy+g48+vo+iCJw18VpWHzxaLlDUjwWNRRwt2WnIMhkwOE6K0qr2+QOh4hIcVq7+nDPK7vR1efAjAti8cQNU+QOSRVY1FDARVuCcH1GEgDgn1wwTET0DXaHEw+8ugfVrd1Ii7Hg2Tunw2zkj2t38F0iWdx19jbqvw7U4kxXn8zREBEpx5PvH8LOk60ICzLihaU5GBHGMxDcxaKGZHFRajSmJEeit9+JN0qq5Q6HiEgR1nx5Cq/uqoIgAH9ZdBEmJETIHZKqsKghWQiCgCUzBu7W/HNXJRxOLhgmIn374ngznvzgMADgsXkTMXtygswRqQ+LGpLNjZmjEBVqRnVrN7aWN8odDhGRbCpbuvBA0V44nCJuzkrG/VeMlTskVWJRQ7IJDTLi9uwUAMArX/I8KCLSp44eO76/ugRtNjsyU6Lwh1szIAg8AsEbLGpIVucWDG8tb0JlS5fM0RARBZbDKeLR10txvLETCZHBeH5JDkLMRrnDUi0WNSSr9LgwXDFhJESRp3cTkf78zydH8enRRgSbDHh+cQ4SIkPkDknVWNSQ7M4tGH6j5DR67A6ZoyEiCox39p7G37eeBAD8z20ZyEyNljcgDWBRQ7KbNTEeKSNC0d5tx/v7eR4UEWnfvqoz+MU7BwEAD8wai5uyRskckTawqCHZGQ3C+bU1a77keVBEpG117d24b80e9PU7MXtSAn46d6LcIWkGixpShAU5qQgyGXCwpp3nQRGRZtn6+nHP6hI0dfRiYkIEVi7KgsHAJ52kwqKGFCEm7N/nQa3h491EpEFOp4gfvV6KQ7VWxIYF4YWlOQgPNskdlqawqCHFWDIjHQDwrwN1aOnslTcYIiKJ/XFDGTYcbkCQ0YC/L85GaoxF7pA0h0UNKUZWajQyUqLQ53BiLc+DIiINeXvPaTy75QQA4A+3TkNOeozMEWkTixpSlMVnFwy/urOK50ERkSaUVJ7B41970mn+9BSZI9IuFjWkKDdkJiPaYkZNWzc+PcrzoIhI3Vp6gAeKStHncOKaKYl80snPWNSQooSYjViYkwoAWMMdholIxTp6+vH8USPO2OyYOioSf1qYySed/Ex3y67tdjvsdvuQX//6fz29pidj3e3vqp8v7b68Xn9bmJOM57efxLbyJhyrb0N6bBiAwMYs1VxqzCtXfdSaV8MJVNxSzuPttZSaV67a1ZZbDqeIH60tRX23gJHhQfhbQRbMguiX+JWQV96M9TS33CGIGt/prLCwEIWFhXA4HCgvL0dRUREsFq44V7q/HzHgcJsBs5KcuCXdKXc4REQeeeeUAVvrDDALIh6Z6kBauNwRqZfNZkNBQQHa29sRGRk5bF/NFzXnWK1WREVFobm5ecg3xW63Y+PGjZgzZw7MZrNH1/Z0rLv9XfXzpd2X1xsIW8qbcO+afYgMMeHzx65AaJAxoDFLNZca88pVHzXn1XcJVNxSzuPttZSaV67a1ZRbr+8+jf/3/mEAwLIJDvxs0WzN55U3Y93t39LSgqSkJLeKGt19/GQ2m4d981y1+3Jtb/v7GvNw7b68Xn+6elIS0mLKUNVqw/pDjViUl3a+LZAxSzWXGvPKVR815pUrgYpbynm8vZZS88pVu9Jz64sTzXjqX0cAAI9eNRYXdJfpKq+8GetOPriLC4VJkQwGAXddPFDIvMLzoIhIBU42deL+f+5Fv1PETVnJeHDWBXKHpDssakixFuSkIthkwOE6K/ZWnZE7HCKi79Rus+Oe1SVo77bjorRo/PetGRAEPukUaCxqSLGiLUG4MTMZwMDdGiIiJbI7nHigaA9ONnchOSoEzy/OQYjZKHdYusSihhRt6SXpAID1B+vQ2MHzoIhIWURRxJPvH8KO4y2wBBnx4rJcjIwIljss3WJRQ4o2dVQUckaPgN0h4rVingdFRMqy+otTeHVXFQQB+OuiizApafinc8i/WNSQ4i27NB0A8Nru0+jnljVEpBCfHW3Eb/418Oj249+7ELMnJ8gcEbGoIcWbNyURiZEhaOnqw74WLrwjIvkdrrXioaK9cIrAwpxU3DuTTzopAYsaUjyz0YDFMwZO795aZ+Dj3UQkq/r2Hix/eTe6+hy4dFwsfnvLVD7ppBAsakgV7shLQ5DJgOouAaXV7XKHQ0Q61dXbj++v3o16aw/GxYfj2TuzYTbyR6lS8E+CVCEmLAg3ZiQBAFbvrJI5GiLSI4dTxKOv78OhWitiw4Lw0rJcRIUqd3djPWJRQ6qx+OJUAMAnhxpQ394jczREpDe//fAwNh1pRLDJgFVLc5Aaw8ORlYZFDanG5KRIjI0Q0e8U8eoubsZHRIGz+otTeGnHKQDAnxZkYXraCHkDoiGxqCFVuTxp4Jnuol1V6LE7ZI6GiPTg06MNeOqDQwCAn10zEded/SiclIdFDanKtBgRSVEDj3d/sL9W7nCISOMO1bbjoaJ95x/dvv+KsXKHRMNgUUOqYhSAO/MG1ta8/MUpPt5NRH5T396D779cAhsf3VYNFjWkOgtzUhBiNuBQrRUllTy9m4ikx0e31Yl/QqQ60RYzbrloFADg5bML94iIpOJwinjkNT66rUYsakiVzp3e/fGhetS2dcsbDBFpyn/+6zA2H+Wj22rEooZU6cLESMy4IBYOp4h/7uTj3UQkjZd3VODlL04BAP68kI9uqw2LGlKt86d3F/PxbiLy3SeH6vHU2VO3f3bNRFw7jY9uqw2LGlKt2ZMSMCo6FGdsdry7r0bucIhIxfZWncEjr+2DKAJ35PHRbbViUUOqZTQIuPvs3ZoXP6/g491E5JVTzV24Z3UJevuduHLiSPznTXx0W61Y1JCqLchNRXiwCccaO7G1vEnucIhIZVo6e7H0pWK0dvVh2qgoPFMwHSY+uq1a/JMjVYsMMWNh7sBmfC9+XiFzNESkJt19DixfXYLKFhtSRoTixWU5CAs2yR0W+YBFDaneskvSYRCA7ceacbTeKnc4RKQCDqeIR17fh/3VbYgKNePlu/MQHxEid1jkIxY1pHqpMRZ87+xTCi9u590aIhqeKIp46oND2Hi4AUEmA15YmoNx8eFyh0USYFFDmnDPZWMAAO+V1qKxo0fmaIhIyZ7fdhKvfFkJQQBWLsxCbnqM3CGRRFjUkCZclDYC2aNHoM/hxJovuRkfEQ3t/f21+P1HRwEA/3HtJO5FozEsakgzzt2t+efOSnT3cTM+IvqmnSdb8NM39gMA7r40HffMvEDmiEhqLGpIM+ZOSURqzMBmfO/sOy13OESkIMcaOnDfKyXoczjxvamJ+NV1k+UOifyARQ1phtEg4O5LBu7WvPh5BZxObsZHREBdezeW/qMY1p5+ZI8egT8vzILRwM31tIhFDWnKgtxURISYcLKpC1vKG+UOh4hk1mbrw5IXi1Hb3oMLRobhhSU5CDEb5Q6L/IRFDWlKeLAJBXlpAIAX+Hg3ka519zlwz+oSHGvsRGJkCNZ8Px8jwoLkDov8iEUNac7SS9JhNAj44kQLDtW2yx0OEcmg3+HEw6/tRUnlGUSGmLB6eR5GRYfKHRb5GYsa0pzk6FBcd/YxTd6tIdIfURTxy3UHselII4JNBry4LBcTEyPkDosCgEUNadK9Zx/VfH9/LU6fsckcDREF0h83lOGNktMwCMAzBdO5uZ6OsKghTZqWEoVLx8XC4RR5t4ZIR17aUYHCz04AAH4/fxrmTE6QOSIKJBY1pFk/vGIsAGDt7mqc6eqTORoi8rf399fiN/86DAB4bN5ELMxNkzkiCjQWNaRZl42Lw5TkSHTbHVj95Sm5wyEiP/r8WDN+8kYpRBFYdkk6Hpg1Vu6QSAYsakizBEE4f7dm9RenYOvrlzkiIvKHg6fb8YM1JbA7RFyXkYRfXz8ZgsDN9fSIRQ1p2vemJiItxoIzNjve2F0tdzhEJLHjjZ1Y+lIxuvocuHRcLP60IBMG7hasWyxqSNNMRgPuvXzgSahV2ytgdzhljoiIpFLdasNdL+xCa1cfpo2KwnN3ZSPYxN2C9YxFDWne7dkpiAsPQk1bNz48UCd3OEQkgUZrD+56cRfqrT0YHx+O1cvzEBFiljsskhmLGtK8ELMRyy5JBwA8t/UERJEHXRKpWZutD4tfLEZliw2pMaFY8/18xPD4AwKLGtKJxRenIyzIiKP1HdhS3iR3OETkpc7efix7aTfKGjoQHxGMV79/MRKjQuQOixSCRQ3pQpTFjDvOHnT53JYTMkdDRN6wO4EHikpRWt2GaIsZa76fj7RYi9xhkYKwqCHd+P7MMTAbBeyqaMWeyla5wyEiD9gdTqwuN+DLk60ICzLi5bvzeJ4TDcKihnQjKSoUt05PAQD8dfNxmaMhInc5nSIeX3cIB88YEGQy4IWluchKjZY7LFIgFjWkKw/MGgejQcDW8ibsr26TOxwickEURfz6/a/w3v46GAQRTy/KxIyxsXKHRQqluKKmuroas2bNwuTJk5GRkYE333zzfFtbWxtycnKQlZWFqVOnYtWqVTJGSmqUFmvBTVnJAICnP+XdGiIlE0UR//mvI/jnzioIAnDnWCeumjhS7rBIwUxyB/BtJpMJK1euRFZWFurr65GdnY1rr70WYWFhiIiIwLZt22CxWNDV1YWpU6di/vz5iI1l1U7ue/DKcVi3rwabjjTgcK0Vk5Mj5Q6JiL5FFEX898dl+MeOCgDAb2+cjPDGAzJHRUqnuDs1SUlJyMrKAgAkJiYiLi4Ora0DizqNRiMsloGV7r29vRBFkXuOkMfGjgzH9RkDd2ue+eyYzNEQ0VD+vOkYnts68KTif948FQtyUmSOiNTA46Jm27ZtuOGGG5CcnAxBEPDuu+8O6lNYWIj09HSEhIQgPz8fxcXFXgW3Z88eOBwOpKamnv9aW1sbMjMzkZKSgsceewxxcXFeXZv07aErxwEAPvqqHuUNHTJHQ0Rf9+yWk/jr5oF/cPz6+slYfPFomSMitfD446euri5kZmZi+fLlmD9//qD2tWvXYsWKFXjuueeQn5+PlStXYt68eSgrK0N8fDwAICsrC/39g09M3rBhA5KTB/4F3draiiVLlgxaNxMdHY39+/ejoaEB8+fPx2233YaEhIRB1+rt7UVvb+/531utVgCA3W6H3W4f1P/c14Zqc8XTse72d9XPl3ZfXq9cpIz5gtgQzJ0cjw2HG/HXTeX484IMv8ylxrxy1UdreQUELm4p5/H2WkrNq3Nf31wj4P2qgfVuP5s3HovzU77xfVtNuaWnvPJmrKe55Q5B9OHzG0EQsG7dOtx8883nv5afn4/c3Fw888wzAACn04nU1FQ8/PDD+MUvfuHWdXt7ezFnzhzce++9WLx48Xf2e+CBB3DVVVfhtttuG9T25JNP4qmnnhr09aKiovMfYZG+ne4C/veACQJE/DLLgfhQuSMi0rctdQLWnRo4kPK6VAfmpnB5AQE2mw0FBQVob29HZOTwayAlXSjc19eHPXv24PHHHz//NYPBgNmzZ+PLL7906xqiKGLZsmW46qqrBhU0DQ0NsFgsiIiIQHt7O7Zt24b7779/yOs8/vjjWLFixfnfW61WpKamYu7cuUO+KXa7HRs3bsScOXNgNnt2KJqnY93t76qfL+2+vF65+CPmPb378GlZEw4LaVh27VTJ51JjXrnqo7W8AgIXt5TzeHstpeZVUXE11n15BABw/8zRWDF3ok9xK4Ge8sqbse72b2lpcTsGSYua5uZmOByOQR8HJSQk4OjRo25dY8eOHVi7di0yMjLOr9dZs2YNpk2bhsrKStx3333nFwg//PDDmDZt2pDXCQ4ORnBw8KCvm83mYd88V+3D8XSsu/19jXm4dl9er1ykjPnR2RPwaVkT3t9fh0evnoD0uDC/zKXGvHLVR2t5BQQubinn8fZaSsqrNTsr8cQHAwXN1clO/HjOBE3llp7yypux7vwMc5fiHum+7LLL4HQ6h2zLy8tDaWlpYAMiTctMjcaVE0fis7ImrNxUjpWLLpI7JCJdeXlHBZ784DAAYPklo5HhPAFBEGSOitRK0ke64+LiYDQa0dDQ8I2vNzQ0IDExUcqpiCTzk7O3ud/bX4uyej4JRRQoL2w/eb6g+eEVY/GLayaA9Qz5QtKiJigoCNnZ2di8efP5rzmdTmzevBkzZsyQcioiyUwdFYXvTU2EKAJ/2lgmdzhEuvD3rSfw2w8HPnJ66Mpx+Pk1E3mHhnzm8cdPnZ2dOH7839vLV1RUoLS0FDExMUhLS8OKFSuwdOlS5OTkIC8vDytXrkRXVxfuvvtuSQMnktKKORPw8aF6fHKoAQdOt2FSQpjrQUTklee2nsT/bRr4OfLo1ePxo9njWdCQJDwuakpKSnDllVee//25J4yWLl2Kl19+GQsXLkRTUxN+/etfo76+HllZWfj444+H3EuGSCnGJ0Tg5qxRWLevBv+3oRwvLObaGiJ/+OS0gPXVAwXNT+ZMwMNXj5c5ItISj4uaWbNmuTya4KGHHsJDDz3kdVD+xM33PItBSfwd84OzxuCD/bXYWt6EnSeaJJlLjXnlqo/W8grQ1yZpcuWVKIr488ZyrK8e2Ifmp3PG4weXp3+jv9ZyS0955c1YxW2+pwaFhYUoLCyEw+FAeXk5N9+jYa09YcAXjQaMjRDx8BQHFy0SScApAu+eMmBr/cAyzptGO3BVsqZ/9JCEPNl8T/NFzTlWqxVRUVFobm7m5ntexKwEgYi5rr0HV/95O+wOEfdPcuDh22dz8z1uvqfIedSy+V6/w4n/eO8w3tlXCwC4bYwDv1k89N8rreWWnvLKm7GebL6XlJQU+B2F1YCb73kXg5L4M+a0ODPuung0XtpxCh9UGfAjo4mb73HzPUXPo+TN93r7HVjx1gF89FU9jAYBf7hlCoJqS/k9SwXzKOl7lifXkvSRbiIteOjKcQgPNuF0l4D3D9TJHQ6RKtn6+nHP6hJ89FU9gowGPHvndNyclSx3WKRxLGqIviU2PBg/vHwMAOBPm46jx+6QOSIidWnvtmPxi8XYfqwZliAjXro7F/OmcANW8j8WNURDWDojDdFBIurae/CPHRVyh0OkGk0dvbjj+Z3YU3kGkSEmrPl+Pi4dFyd3WKQTLGqIhhBiNuL6tIEzyJ797ARaOntljohI+Rq7gQXP78LhOiviwoOx9gczkD16hNxhkY6wqCH6DtlxIqYkR6Cztx9/2XxM7nCIFG1fdRtWfmXE6bYepMda8Pb9MzApafgnVYikprunn7j5nmcxKEkgY7bb7TAIwE+uHovla0pRtKsKd+am4IKRnh2foMa8ctVHa3kF6GuTNH/k1eajjfjR2gPo6RcwLTkCqxZPR2x40KAx/J6l/HmU/D3LHZrfp4ab75Gvnj9qwKEzBkwZ4cR9FzrlDodIUb5oEPDGSQNECJgU7cTdE5wINsodFWkJN98bAjff40ZW3s5V1daH65/5Av1OEasWX4RZE0Z6dR215JWrPlrLK0Bfm6RJlVeiKOLpz07g6c9OAgDmZyXhsuBqXDPPu7xy1a7G3NJTXnkzlpvvSYCb73kXg5IEMmaz2YwLky1YftkYPL/tJH63vgyXT0xAsMmzf4qqMa9c9dFaXgH62iTNl7zq7Xfgl+u+wtt7TwMAHrlqHB6aNQYffVTtc165aldjbukpr7wZy833iALs4avGIT4iGKdabHhhOx/xJv1q7erD4heK8fbe0zAaBPzulqlYMXciBB6URgrAoobIDREhZvzy2kkAgGc+PY7atm6ZIyIKvOONHbi5cAeKT7UiItiEfyzLxZ35o+UOi+g8FjVEbropKxm56SPQbXfgd+uPyB0OUUBtP96MW579AlWtNqTGhOKdBy7BFR6sLyMKBBY1RG4SBAFP3jgFBgH48EAdth9rkjskooDYXi/g3jX70NHTj9z0EXj3gUsxPiFC7rCIBmFRQ+SBKclRWDIjHQDwy3UHYevrlzcgIj/q7XfgV+8dxlsVRjicIuZPH4V/3pOP2PBguUMjGhKLGiIP/XTeRCRHhaC6tRsrN3GnYdKm+vYeLHp+J9aWnIYAET+dMx7/d3umx0/+EQUSixoiD4UHm/DbW6YCAF7YfhIHT7fLHBGRtIorWnH9059jX1UbIkNMuO9CJ35w+Rg+4USKp7t9anhMgmcxKEmgj0kYbq6ZY2Nw3bREfHiwHj97az/e/mE+zMbB/0ZQY1656qO1vAL0tZ39cONEUcSaXdX4/Udl6HeKmJgQjr/cPhVlez73e165aldjbukpr7wZy2MSvMBjEshfOuzAf5UaYesXcH2aA3NGafqvEmlcrwNYe9KAPc0Dxfn0WCcWjeWRByQ/HpMwBB6TwC3H/THXO/tq8PN3DsFsFPDu/RdjwreeCFFjXrnqo7W8AvS1nf1Q447UdeDRtftR0WKD0SDgZ3PH4+5LRkMQhIDllat2NeaWnvLKm7E8JkECPCbBuxiUJNDHJAw314Lc0fjkcBM+PdqIn759CO8+eMmQCynVmFeu+mgtrwB9bWdvNpthMplQVFyFpz44jL5+JxIjQ/B0wUXITY/xeh4ekzCYnvLKm7E8JoFIIQRBwB9unYYRFjOO1Fn5NBSpRkdPPx5+bR/+Y91X6Ot34qoL47H+0ZlDFjREasGihshH8REh+P38aQCA57aewO5TrTJHRDS8Ux3AzX/7Ev86UAeTQcAvr70QLyzJQUxYkNyhEfmERQ2RBK6ZmoTbslMgisCP15aivVs9T2iQftgdTvx503Gs/MqIqtZujIoOxRs/nIH7Lh8Lg4GPa5P6saghksgTN0xGakwoTp/pxs/e2g+drMEnlTje2IFbnt2BZ7eehAgBN2UmYf2jMzE9bYTcoRFJhkUNkUQiQswoLJiOIKMBnxxqwD92nJI7JCI4nSJe2lGB6/76Ob6qsSIq1IRl4x34423TEBWqrgW3RK6wqCGSUEZKNH51/SQAwO/XH8G+6jZ5AyJdO9nUiTtW7cRTHxxGb78Tl08YiQ8fugQXxfEuImmT7h7pJvK3xRePxq6TrfjwYB0eXXsAD42XOyLSG7vDiVXbT2LlpmPo63ci1GzE49deiMUXj0Z/Pw9hJe1iUUMksXOPeR+us6KiuQsvlhlxU78TKttag1Tq4Ol2/PztAzhcZwUAzBwfh/+6ZRpSY7iTOmkfP34i8oOIEDNWLclGeLAJJzsEPPHBYS4cJr/q7O3H7z48jJsKP8fhOiuiLWb83+2ZeGV5Hgsa0g3d3anhgZaexaAkSjrQ0h2jR4TgT7dOxg+K9uPtvbUYNzIc91yW7rcYeKCl99R88GBfXx/e21+H//mkHI0dvQCA66cl4lfXTkRsePCgj5uUmleu2tWYW2rOKyX+LHSH5s9+4oGWJLetdQLeOWWEABF3T3AiM1bTf+UogGq6gLcqjDjZMbDHTFyIiPnpTkwZwRwj7eCBlkPggZY8HE6Ouex2OzZs2Igv7WlYu6cWZqOAFxZPxyVjYyWPgQdaek9tBw+2dPbiL5uP4/WS0xAhINRswANXXIC7L01HsGn4VQVKzStX7WrMLbXlla/X4oGWMuCBlt7FoCRKOtDSHYIAPHXjFLT3OPHxoXo8UFSKonsvRmZqtF9i4IGW3lP6wYO2vn68sL0Cf996Al19DgACrp2agF9dPwXJ0aF+jYEHWnpP6Xkl9bV4oCWRxhkNAv5yRxYuHReLrj4Hlr1UjPKGDrnDIpXodzhRtKsKV/zvFvxpYzm6+hzIGBWJhyf34y8LMz0uaIi0ikUNUYAEm4z4++IcZKZE4YzNjkXP78ThWqvcYZGC9TucWLfvNOb+eRt+ue4gmjp6kRZjwTMFF+GtH+RjXJTcERIpi+4+fiKSU3iwCauX52Hxi8U4WNOOO1btxCvL89z+KIr0od/hxPv7a/H0p8dR0dwFABhhMeORq8fjzvzRCDIZVPUUEFGgsKghCrBoSxBevTcfy/5RjL1VbbjrhV14cVku8sbEyB0ayay334H3S2vx7JYT3yhm7r38AiyZkY7wYH7LJhoO/4YQySAyxIxXvp+P5S/vRnFFK+56YRf+9/YM3JQ1Su7QSAbtNjteLa7EyztOnd9rhsUMkef4N4VIJuHBJqy+Ow8/WrsPnxxqwKOvl+JUsw2PXD0OgiDIHR4FQHWrDS9+XoE3Sqph63MAABIig3H3pWNw18WjWcwQeYh/Y4hkFBpkxN/uzMZ/f3wUf992En/eVI7jTZ34/fxpCOYyfk1yisCnZU14o6QGn5U1wnl2p7ALEyNw78wLcENmMoJc7DVDRENjUUMkM4NBwOPXTkJ6XBj+37tf4YP9tThU046VCzLkDo0k1GjtQdGuSqzea8SZnfvOf33m+DjcO/MCzBwfxzt0RD5iUUOkEHfkpWFCQjgeLtqHk81duO35XbgxVcA1Tl1s+q1JPXYHNh9pxLp9NdhS1oh+pwhAQHSoGbdlp+CO/DSMHRkud5hEmsGihkhBskfH4MNHZuInb+7Hp0cb8VaFEdUvl+C/b81EelyY3OGRG5xOEbsqWrFu32l8dLAeHb3/PlAyOy0ak4Ja8IuCqxFuCZExSiJtYlFDpDAjwoLwwpIc/OPzE/ifj49iV8UZzFu5DT+eMwHLLx3D9RYK1O9woriiFZ8cqscnhxpQb+053zYqOhQ3ZSXjlotGIT0mBOvXr0ew2ShjtETapbuixm63S36svb+PW/+ufr60S3lUfaAEMmap5vLlOgU5yRDqDmFj20jsrDiDP3x0FEW7KvHY3AmYNzl+0PoLqfLKVR+t5RXgXdzdfQ58cbIFG4804tOjTThj+/fYiBATvjclATdlJSEnbQQMBsHreaSM2ZtxgcorV+1qzK1AxayEvPJmrKe55Q7Nn9JdWFiIwsJCOBwOlJeXo6ioCBaLRe6wiNwmikBxk4B/VRlgtQ/8cBwTIeKaFCcmRong2tLAEEWgrhs42ibgSJuAk1YB/eK/3/wwk4ipI0RkxoqYECXCzBtqRJKw2WwoKChw65RuzRc151itVkRFRaG5uXnIN0WJx6276udLu5RH1QdKIGOWai4p86qrtx+rPj+FF3ecQo/dCQDISInE/ZdfgCsnjoTT0S9JXrnqo7W8AoaOWxRFVJ/pRknlGRSfOoPPj7egwdr7jXHJUSG4elI85k2OR3ZaNEzG4SsZKd8fb68l1/crd/poLbcCFbMS8sqbse72b2lpQVJSkltFje4+fvLlyHtfr+1tf19jHq5dyqPqAyWQMUs1lxR5FW0247FrJmHxjDF4busJvL67CgdOW3F/USlGRYdiYc4ojOiTLq9c9dFaXjlF4GRLD/adbsSuilbsPtU6qIgJMRtw8QWxuHz8SFw+YSTGjgzz6jFsKd8fb68l1/crd/poLbcCFbMS8sqbse7kg7t0V9QQqV1iVAievHEKHrpqHP7xeQVe3VWFmrZu/GnTcRgEIzZ17MUNmaMwZ3ICokLV9c0/UJxOERUtXTh4uh0Ha9qxv/oMDlQb0bfzy2/0MxsFTBsVhdwxMbhsXBxy02MQwkW+RIrFooZIpeLCg/Gzay7EI1ePx4cH6vDqrkrsrWrD1vJmbC1vhtko4NJxcbh8/EjMHB+HcfHhutvcTRRFNHX24lhDJ8obOnCssRPHGjpwpK4DnV971HqAAEuQEdPTRiA3PQZ5Y2KQlRqN0CAWMURqwaKGSOVCzEbcmp2CGzMS8NLb69EVMxEfH2pEWUMHtpQ1YUtZEwAgMTIEuWNiMC05AraOgY3h1HYbfyiiKOKMzY6qVtvAr5YuVLXacKrZhvLGDrTZhn5yIthkwJTkSGSkRGNSYhjOnNiPZfPnICQ4KMCvgIikwqKGSEMSQoFrrxyLH8+9EMcbO/Dp0UZsP9aM4opW1Ft78MH+WnywHwBMeObwpxg7Mhxj48MwdmQ4xsWHY0xcGJKjQxFhVsYdHadTxBlbHxqsvWjo6EGjtWfg/8/+t6atG9WttiHuuvybIACjYywYnxCB8fHhmJAQgYmJA/9/bmGv3W7H+rr9MBqU8bqJyDssaog0alx8BMbFR+C+y8eix+7Ansoz2Fc18Kv4RCM67EBZQwfKGjoGjTUbBUSajFhTW4yEqFBEhZoRGWJGVKgZYUECTjQLCD7SiNCQIAQZDQgyCQgyGiGIDtR0AYfrrDAYTHCKIhyiCFEU0dvvRHefA912B7r7HOixO2DrG/hl7bGjzWZHm60Pbd12tNvsA//ttsPh5jERiZEhSIuxIDXGgrQYC0bHWjAufqBY4zoYIn1gUUOkAyFmIy4dF4dLx8XBbrfjww/X46JLr8TJ1h6caOzEiaZOnGjsQkVLF5o7e2F3iGhxCGipbAPQNsQVjVh9rPQ7ZjMBB3ZKGn9ceBDiI0KQEBmMhMgQxEcO/H9SVAjSYsKQMiKUhQsRsagh0iNBAJKjQzF6ZCSunBj/jTa7w4na1k6s++QzjJk8Ha3d/bB296O92z5wR6WrF6dqGxEWGQ27Q0Sfwwm7w4m+/oFf9r5ehIaEwGAQYBAECAJgEAQEmwwIDTIixGyEJciIUPPAr5AgI6JCzYgONWOEJQhRloH/j7YEIdpiRkxYEMwu9n8hIgJY1BDRt5iNBiRHh2JMBHDttMQhN0Fbv349rr02f5i2KzSxCJmI1IX//CEiIiJNYFFDREREmsCihoiIiDSBRQ0RERFpAosaIiIi0gQWNURERKQJunuk2263w24ffBbMua8N1ebONT0Z625/V/18affl9colkDFLNZca88pVH63lFRC4uKWcx9trKTWvXLWrMbf0lFfejPU0t9whiKLo3h7kKlVYWIjCwkI4HA6Ul5ejqKgIFotF7rCIiIjIDTabDQUFBWhvb0dkZOSwfTVf1JxjtVoRFRWF5ubmId8Uu92OjRs3Ys6cOR5vGubpWHf7u+rnS7svr1cugYxZqrnUmFeu+mgtr4DAxS3lPN5eS6l55apdjbmlp7zyZqy7/VtaWpCUlORWUaO7j5/MZvOwb56rdl+u7W1/X2Mert2X1yuXQMYs1VxqzCtXfbSWV0Dg4pZyHm+vpdS8ctWuxtzSU155M9adfHAXFwoTERGRJrCoISIiIk1gUUNERESawKKGiIiINEE3C4XPPeRltVqHbLfb7bDZbLBarV6t+PZkrLv9XfXzpd2X1yuXQMYs1VxqzCtXfbSWV0Dg4pZyHm+vpdS8ctWuxtzSU155M9bd/h0dHQD+/XN8OLopas69KampqTJHQkRERJ7q6OhAVFTUsH10s0+N0+lEbW0tIiIiIAjCkH1yc3Oxe/dur67v6Vh3+7vq52271WpFamoqqqurXT73ryS+/BnJNZca88pVH63lFRC43JJyHm+vpdS8Gq5drbmlp7zyZqw7/UVRRHZ2NsrLy2EwDL9qRjd3agwGA1JSUobtYzQavf7L4ulYd/u76udre2RkpKq+QfjyZyTXXGrMK1d9tJZXQOByS8p5vL2WUvPKnXa15Zae8sqbse72DwoKclnQAFwo/A0PPvhgwMa6299VP1/b1SaQr0equdSYV676aC2vgMC9Jinn8fZaSs0rT+ZSCz3llTdjpcwtQEcfP9E3nTs2wp1tp4ncxbwif2FukTt4p0angoOD8cQTTyA4OFjuUEhDmFfkL8wtcgfv1BAREZEm8E4NERERaQKLGiIiItIEFjVERESkCSxqiIiISBNY1BAREZEmsKihQW655RaMGDECt912m9yhkIZUV1dj1qxZmDx5MjIyMvDmm2/KHRJpQFtbG3JycpCVlYWpU6di1apVcodEMuIj3TTIli1b0NHRgdWrV+Ott96SOxzSiLq6OjQ0NCArKwv19fXnz3IJCwuTOzRSMYfDgd7eXlgsFnR1dWHq1KkoKSlBbGys3KGRDHinhgaZNWsWIiIi5A6DNCYpKQlZWVkAgMTERMTFxaG1tVXeoEj1jEYjLBYLAKC3txeiKIL/VtcvFjUas23bNtxwww1ITk6GIAh49913B/UpLCxEeno6QkJCkJ+fj+Li4sAHSqojZW7t2bMHDocDqampfo6alE6KvGpra0NmZiZSUlLw2GOPIS4uLkDRk9KwqNGYrq4uZGZmorCwcMj2tWvXYsWKFXjiiSewd+9eZGZmYt68eWhsbAxwpKQ2UuVWa2srlixZgueffz4QYZPCSZFX0dHR2L9/PyoqKlBUVISGhoZAhU9KI5JmARDXrVv3ja/l5eWJDz744PnfOxwOMTk5Wfz973//jX6fffaZeOuttwYiTFIhb3Orp6dHnDlzpvjKK68EKlRSEV++Z51z//33i2+++aY/wyQF450aHenr68OePXswe/bs818zGAyYPXs2vvzySxkjI7VzJ7dEUcSyZctw1VVXYfHixXKFSiriTl41NDSgo6MDANDe3o5t27Zh4sSJssRL8mNRoyPNzc1wOBxISEj4xtcTEhJQX19//vezZ8/G7bffjvXr1yMlJYUFD7nkTm7t2LEDa9euxbvvvousrCxkZWXh4MGDcoRLKuFOXlVWVmLmzJnIzMzEzJkz8fDDD2PatGlyhEsKYJI7AFKeTZs2yR0CadBll10Gp9MpdxikMXl5eSgtLZU7DFII3qnRkbi4OBiNxkGL6BoaGpCYmChTVKQFzC3yB+YVeYpFjY4EBQUhOzsbmzdvPv81p9OJzZs3Y8aMGTJGRmrH3CJ/YF6Rp/jxk8Z0dnbi+PHj539fUVGB0tJSxMTEIC0tDStWrMDSpUuRk5ODvLw8rFy5El1dXbj77rtljJrUgLlF/sC8IknJ/fgVSeuzzz4TAQz6tXTp0vN9nn76aTEtLU0MCgoS8/LyxJ07d8oXMKkGc4v8gXlFUuLZT0RERKQJXFNDREREmsCihoiIiDSBRQ0RERFpAosaIiIi0gQWNURERKQJLGqIiIhIE1jUEBERkSawqCEiIiJNYFFDREREmsCihoiIiDSBRQ0RERFpAosaIiIi0oT/Dy8ENUljyQhrAAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"result = virgo.get_QNLS([5, 5000, 1000])\n",
"\n",
"Figure = figure()\n",
"_ = Figure.gca().loglog(result.x0, abs(result[\"NSR_with_RP\"]))\n",
"Figure.gca().grid(True, \"both\", \"both\")\n",
"show()"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "8261c88c-8104-4ef9-9a0f-c186891c8b27",
"metadata": {},
"outputs": [],
"source": [
"save(\n",
" \"sensitivities/finesse-virgo.npy\",\n",
" array([result.x0, abs(result[\"NSR_with_RP\"])]),\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "e6519eb5-8dbc-4656-8036-3a4522787cff",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/demagny/codes/python/finesse/finesse-simulation-04/.venv/lib/python3.13/site-packages/finesse/analysis/actions/base.py:157: UserWarning: Signal frequency (fsig) was set to None but simulation needs it. Setting default value of 1 Hz\n",
" self.__sim.__enter__()\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAi4AAAGhCAYAAABGRD9PAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAWZFJREFUeJzt3XtclFX+B/DPM8MAIhdFBARBvN8FUxi1LGnxVtmqWW4XQ00zQ7uw5epumba7tbtdllLK0l3tHmortFnmJRULksQw74jiFUHRcASUy8z8/nDll3GZZ555Zp5n5vm8X6/f6/eCOec537Gzx+NzvuccwWq1WkFERETkBnRKB0BEREQkFicuRERE5DY4cSEiIiK3wYkLERERuQ1OXIiIiMhtcOJCREREboMTFyIiInIbXkoHIDeLxYKSkhIEBARAEASlwyEiIiIRrFYrLl++jIiICOh0zb9X8biJS0lJCaKiopQOg4iIiCQ4deoUOnbs2OznHjdxCQgIAHDtiwcGBjb6vK6uDhs3bsSoUaNgMBhEP9feemLLiylnq0xLn0v9vkpyVcxytuPIs5TqW+xX6m/LHfuVrTKe1q8AbY1ZzuxXmZmZmDFjRsPf483xuInL9eWhwMDAZicufn5+CAwMtPs/lj31xJYXU85WmZY+l/p9leSqmOVsx5FnKdW32K/U35Y79itbZTytXwHaGrOc3a8A2EzzYHIuERERuQ1OXIiIiMhtcOJCREREboMTFyIiInIbHpece11dXR3q6uqa/P0v/789z7OnntjyYsrZKtPS51K/r5JcFbOc7TjyLKX6FvuV+ttyx35lq4yn9StAW2OWs/uVGILVarWKLq1i6enpSE9Ph9lsRmFhIT7++OOGDGUiIiJSt+rqajzwwAO4dOlSk7uCr/OYict1JpMJQUFBKC8vb3Y79KZNmzBy5Ei7t4DZU09seTHlbJVp6XOp31dJropZznYceZZSfYv9Sv1tuWO/slXG0/oVoK0xy5n9KisrS9TExWOXigwGQ4t/SLY+l6ue2PJiyjnynaR+XyW5KmY523HkWUr1LfYr9bfljv3KVhlP61eAtsYsZ/UrMZicS0RERG6DExciIiJyG5y4EBERkdvgxIWIiIjcBicuIl2pNeOniwI8bBMWERGRW+HERaQ3vinCvw7rkfLJHpwzXVU6HCIiIk1S5cSluLgYiYmJ6NOnD/r374+qqiqlQ4K/jxd0ghWbDp5D0uvbsXrXKb59ISIicjFVTlymTp2KF198EQcOHMD27dvh4+OjdEiYk9gVz/Y3o39kIExX6zFv7U94+N95OHWxWunQiIiINEN1E5f9+/fDYDBg+PDhAIDg4GB4eanjnLyI1sDqmQn44x294OOlw44j5Ridlo2V3xXDbOHbFyIiImeTfeKSnZ2NcePGISIiAoIgIDMzs1GZ9PR0xMTEwNfXF0ajEXl5eQ2fHTlyBP7+/hg3bhxuuukmvPTSS3KH6BAvvQ6P3toVG566FQmdg1Fda8bi/x7Afe/koujcZaXDIyIi8miyv8qoqqpCbGwspk+fjokTJzb6PCMjA6mpqVi2bBmMRiPS0tIwevRoHD58GKGhoaivr8eOHTtQUFCA0NBQjBkzBvHx8Rg5cmST7dXU1KCmpqbhZ5PJBMD5t0N3DPLGB1MH4dNdp/GPjYXIP/Ezxr6xA3NGdMXM4TGAxSyqHd4O3ZiWblqVUpe3Q0vD26HlKc/boRvT0pjl8bdDC4KAdevWYfz48Q2/MxqNiI+Px9KlSwEAFosFUVFRmDt3LubPn4/c3FwsWrQIX3/9NQDglVdeAQA8++yzTbaxaNEiLF68uNHvXXk79M81QMYxHQ5WXHuBFelnxf1dzYjyd0nzREREbk/s7dAuTR6pra1Ffn4+FixY0PA7nU6HpKQk5ObmAgDi4+Nx7tw5/PzzzwgKCkJ2djZmzZrV7DMXLFiA1NTUhp9NJhOioqKQmJjY5Bevr6/H1q1bkZiYaFfujK1691qt+O/eMrz89RGcqa7HP/d5IbGDBS8/NBytfb0lP1dMmZY+l/p9leSqmOVsx5Fn2VtXbHlH+o2tz9mvXNOWO/YrW2U8rV8B2hqznNmv1q9fLyoGl75xKSkpQWRkJHJycjB06NCGcvPmzcP27duxc+dOAMBXX32FefPmwWq1YtSoUXj99ddttpWeno709HSYzWYUFha69I3LL12uAz4r1uHHC9fevoT6WvG7rmZ0bX7ySEREpHmqfOMi1tixYzF27Fi76qSkpCAlJQUmkwlBQUEYNWpUk1+8rq4OmzZtwsiRI+26YtueepMBfLW3BM+t24tzVwUsOeCFhxKikDqyO/x9bvwjF/NcW2Va+lzq91WSq2KWsx1HnmVvXbHlHek3tj5nv3JNW+7Yr2yV8bR+BWhrzHJmv8rKyhIVg0u3Q4eEhECv16OsrOyG35eVlSE8PNyVoThdUq/2WBBnxj0DO8BqBT7YeQp3LsnBjiPlSodGRETkthRJzk1ISMCSJUsAXEvOjY6Oxpw5czB//nzJballqagphysEfHpMh4s1AgAgob0F4ztZ0Np9/kFBRETkVGKXimSfuFRWVqKoqAgAMHDgQLz++utITExEcHAwoqOjkZGRgeTkZLzzzjtISEhAWloaVq9ejUOHDiEsLMzh9q8vFZWXlyu2VNRU+aqaevxzSxHe//4krFYgxN8bL9zVG7/pEcylol/R0mtXKXW5VCQNl4rkKc+losa0NGY5e6lIkRyXXbt2ITExseHn6zt+kpOTsWrVKkyePBnnz5/HwoULUVpairi4OGzYsEGWSYuatfbxwnN39MId/cKxYN1+HCuvwtxP92BU7/YYro4XQ0RERKrn1KUiV1LzUtGv1VmAjad12HxGgAUC/PRWTIixIL69FYKgdHRERESup9hSkdLUulTUlANnTZj/n304WFoJABjerR1evLsPOrZtZdezPO3Vq5Zeu0qpy6UiabhUJE95LhU1pqUxSw1LRaq7ZFFL+nQIRMaMQbgr2gxvvQ47ii7gjiXfYWXOCV7aSERE1ASPeePiTktFTSm7AmQc1ePo5WtrRdGtrx1cF9la4cCIiIhcgEtFbrBU9Otyer0XVuefwd+/LkRlTT28dAJmDo/BozdHIXvrN3ylr+J23PGVPpeK1N+WO/YrW2U8rV8B2hqz1LBUpMqTc+VgMBha/EOy9blc9cSWv15uyrDOGNWvA57P3IeNB8rw9vZifL2/DHeFAXc48J2kfl8luSpmOdtx5FnO7lvO+Jz9yjVtuWO/slXG0/oVoK0xy1n9SgyPnbjU1dU1eU22vVdyS60ntnxT5YJb6ZF+fyy+3l+GxV8cxLHyarxZ7oWyzH34w5ieCPBt/K+U5tqS+n2V5KqY5WzHkWe5sm/J9Tn7lWvacsd+ZauMp/UrQFtjlrP7lRges1Tk7jkuzamuBz4/oUPuuWt51EEGK+7tYkH/YI/4z0ZERASAOS5umePS0nrxW2s34/NSf5y8eAUAMKZvGBbe2QvtA3w8bs1YS+vFUuoyx0Ua5rjIU545Lo1pacxijosTuWuOS3O6B1mxftIwvJV9HO9mH8OG/WXIOXoBz93ZB+Njw2w+wx3XjLW0XiylLnNcpGGOizzlmePSmJbGLOa4OIE757g0V0YPC1J/0xVj+rTHHzP3Y3/JZcz77Cd8lt8GI9t4zpqxltaLpdRljos0zHGRpzxzXBrT0pjFHBcZeWqOS3PMVmD7WQFfntKhziLAIFgxNsqCERFW6HltABERuRnmuHhYjktzZU5crMZzmfvxffHPAIA+HQLw0vi+6BsRKCluNdDSerGUusxxkYY5LvKUZ45LY1oas5jj4kSeluPSXJluYUF4f9pgLHpvA9aX+ODA2cu4552dmDG8M576TY+G8u64Zqyl9WIpdZnjIg1zXOQpzxyXxrQ0ZimZ48K7ijyAIAgwhlqx4YmbcdeADjBbrHhn+zGMeSMbuccuKB0eERGRbDhx8SAh/j5Y+sBNWP7wYIQH+uLEhWo8vDIfHxfpcOmKeyW7ERERNcVjl4o8cVeR2N0fI7oH48u5w/DapiP4KO8Udp7XYfQb3+GFu3phTN8wCIK6s3e1lKEvpS53FUnDXUXylOeuosa0NGZxV5GMtLarSKyjJiDjmB5lV65NVvq3tWBSZwva+CgcGBER0S9wV5FGdhXZ+ryurg5ffb0JRT7dsOK7E6gzW+Hv44VnR3XH7wZ3hE6nvrcvWsrQl1KXu4qk4a4iecpzV1FjWhqzuKvIibSyq0jM5146IHVkD4y/KRp/+OwnFJyqwAv/PYgv9pbi5YkD0C3U3/YXUYCWMvSl1OWuImm4q0ie8txV1JiWxizuKiKX6BkegM9mD8ML4/rAz1uPH47/jDve2IElW46gtt6idHhEREQ2ceKiMXqdgGk3d8bGp2/FiJ7tUWu24LVNhRi35FvsPvmz0uERERG1iBMXjerY1g8rp8bjjd/FIbi1Nw6XXcY9b+dgYdY+XL7qXhn9RESkHZy4aJggCPhtXCS2pN6Ge27qCKsVeD/3BEa+no1NB8qUDo+IiKgRj03O1fI5LvY+399bwN8m9MG4AWF4PusATv18BTPf34XRfUKx8K7eCA1w7d5pLZ2JIKUuz3GRhue4yFOe57g0pqUxi+e4yIjnuMij1gxsOK3D1hIBFghopbfi7k4WDAm1QoU7p4mIyEPwHBee4yIp7usOnDXhuawD2HvGBAAY3KkN/vLbvujavrXoZ0ilpTMRpNTlOS7S8BwXecrzHJfGtDRm8RwXJ+I5LtLiuC42uh0yU27BqpzjeG3jYew6UYG703ORktgNs0d0hbeX89OjtHQmgpS6PMdFGp7jIk95nuPSmJbGLJ7jQqqk1wl45JYbt07/c3Mh7nxzB3Ydv6h0eEREpEGcuJBN17dOv3n/QIT4e+PIuUpMWpaL5zL3wsSt00RE5EKcuJAogiDg7tgIbE69DfcN7ggA+PD7kxj5+nZs2FeqcHRERKQVnLiQXdr4eeMfk2Lx8UwjYtr5ocxUg8c+zMesD3ahzHRV6fCIiMjDceJCkgzrGoINT92Kx0d0hZdOwNf7y5D02nZ8+P0JWCwetVGNiIhUhBMXkszXoMe8Mb3w37m3IDaqDS7X1OO5zH24751cHCm7rHR4RETkgVQ5cYmJicGAAQMQFxeHxMREpcMhG3p3CMR/fnHr9K4TP+OON3fgn5sKUVNvVjo8IiLyIKo9xyUnJwf+/v5Kh0EiXb91elTfcCzM3Icth87hjS1H8MVPJXh54gAkdA5WOkQiIvIAqnzjQu4rsk0rrEgejKUPDESIvw+Onq/Cfe/kYsF/9uLSFW6dJiIix8g+ccnOzsa4ceMQEREBQRCQmZnZqEx6ejpiYmLg6+sLo9GIvLy8Gz4XBAG33XYb4uPj8dFHH8kdIjmZIAi4a0AEtqTeht/FRwEAPsm7tnX6q71n4WG3TBARkQvJPnGpqqpCbGws0tPTm/w8IyMDqampeOGFF7B7927ExsZi9OjROHfuXEOZb7/9Fvn5+fj888/x0ksv4aeffpI7THKBID8D/nbPAHz66BB0CWmNc5drMPuj3Xj0g3ycvXRF6fCIiMgNyZ7jMnbsWIwdO7bZz19//XXMnDkT06ZNAwAsW7YM69evx7///W/Mnz8fABAZGQkA6NChA+644w7s3r0bAwYMaPJ5NTU1qKmpafjZZLp2KWBdXZ2s16ar9Yp4W5+r4Zr4QVGB+PzxIXhrezHe3VGMTQfKkHO0HM+M7I4H4qOg+9W101q6Il5KXbn6lrv3K3u5Mma52nLHfmWrjKf1K0BbY5az+5UYTr0dWhAErFu3DuPHjwcA1NbWws/PD2vXrm34HQAkJyejoqICWVlZqKqqgsViQUBAACorK3Hbbbdh2bJliI+Pb7KNRYsWYfHixY1+//HHH8PPz88ZX4scUFINZBzV43jltclKjL8Vk7uaEcH/VEREmlZdXa2+26HLy8thNpsRFhZ2w+/DwsJw6NAhAEBZWRkmTJgAADCbzZg5c2azkxYAWLBgAVJTUxt+NplMiIqKQmJiYpNfvL6+Hlu3bkViYiK8vMR/fXvriS0vppytMi19LvX7OtNDFisy8s/gn1uO4XilGa/tNWDGzdGYNbwTfLz0LotZznYceZZSfcvT+pUtroxZrrbcsV/ZKuNp/QpwXdxqGLOc2a/Wr18vKgaXvnEpKSlBZGQkcnJyMHTo0IZy8+bNw/bt27Fz507JbaWnpyM9PR1msxmFhYV84+IGKmqANcU67Pv5WqpVqO+1ty/dmp9oExGRh1LlG5eQkBDo9XqUlZXd8PuysjKEh4c79OyUlBSkpKTAZDIhKCiIb1wkxu1qk6xWbDp0Hn/56gjOVdZiyX4vTIwLR7zhNO4c6fn/epFSl29cpOEbF3nK841LY3zj4nh51b5xAQCj0YiEhAQsWbIEAGCxWBAdHY05c+Y0JOdKwTcu7q26HvjvCR1yzl17+xJgsOKezhbEBVshCDYqExGR2xP7xkX2iUtlZSWKiooAAAMHDsTrr7+OxMREBAcHIzo6GhkZGUhOTsY777yDhIQEpKWlYfXq1Th06FCj3Bcprr9xKS8vb/KL19XVYdOmTRg5ciQMBoPo59pbT2x5MeVslWnpc6nfVyl5xy/iucwDKL5QDQBI7BmCxeP6oEOQr+xtyfln48izlOpbWupXgGtjlqstd+xXtsp4Wr8CXBe3GsYsZ/arrKwsZZaKdu3adcP9QtcTZ5OTk7Fq1SpMnjwZ58+fx8KFC1FaWoq4uDhs2LBBlkkLub+EmGCsmxWPP7y/FVtK9Nh6uBx5xd8hdWR3PJgQBb2Or1+IiLTMqUtFrsSlIs9TWg18ekyP4svXJiud/K34XRczIlorHBgREclOsaUipXGpyP1fvf4yZr3eC5/sOo1XNx5BZU09vHQCZt4Sg5QRXeBj0MvWDpeKtNWvuFQkvTyXihrjUpHj5e1ZKuIli6RqOp2ABxOi8NUTwzCydyjqLVa8nV2Mu9Jz8f2xi0qHR0RELuYxb1y4VKQNey4IWFusg6nu2vKRsb0Fv+1kQWv3+ccZERE1gUtFXCqSFLca2Ir58tU6vLrpCD7OOw0AaNfaG8/d0RN39g+HYMfeaTW8dpVSl0tF0nCpSJ7yXCpqjEtFjpdXdFeRWhgMhhb/kGx9Llc9seXFlHPkO0n9vkpqLuZggwEvTYzFxJuiMP8/e1F0rhJPr9mL/+4tw5/H90Nkm1aytCNnzM6oK1ffYr9Sf1vu2K9slfG0fgW4Lm41jFnO6ldieOzEpa6Ot0PbE4eaiI05NjIAmbOH4N0dxXh7+zF8c+gcvn99O55O6oYpxmibW6fVcNOqlLpy9S32K/W35Y79ylYZT+tXgOviVsOY5ex+JYbHLBUxx0XbSquBjGN6HPvf1uno1lb8rqsZkdw6TUTkFpjjwhwXSXGrgdSYLRYrVuefwT82FuLy1XrodQJm3ByDOYld4NvE1mk1rBdLqcscF2mY4yJPeea4NMYcF8fLM8cFzHGRGoeaSIl5yrDOGN2vA174fD++2leKd3YUY8OBMrw0oT9u7hYiWzvNccdcBPYr9bfljv3KVhlP61cAc1yc0U5TeI4LeZzQQF+8/dAgvDtlEMIDfXHiQjUeXLETz6zZg5+rapUOj4iIHOCxb1yYnGtfHGoiV8yJPdrhy7nD8PrmI/go7xTW5p/GN4fK8NwdvXBX/3DU19fL0o6jMTM51zWYnCtPeSbnNsbkXMfLMzmXybn0K8WXgU+P6lF65Vrybu82FtzXxYJgH4UDIyIiAEzOZXKuxLjVwFkx19Zb8O6OYry1/RjqzFa0MugxJrIWLz70G/j6eDv0bHdMomS/Un9b7tivbJXxtH4FMDlXjvJMzgWTc6XGoSZyx2wwAE+P6oVxcR3xx//sRd7xi1h3XI+jK3fjb/fEok9E8/9DEd+G+yVRsl+pvy137Fe2ynhavwKYnOuMdprC5FzSnG6h/vj00SH489190EpvxU9nTBi39Fv8fcMhXK0zKx0eERG1gBMX0iSdTsDv4jtiQZwZo/uEwmyx4u1tRzEmLRs5ReVKh0dERM3w2KUi7iqyLw41cWWGfpA38M87+2J7UQQWfXEQxy9U44EVO3HPTRGYP7on2viJe7Xpjrs/2K/U35Y79itbZTytXwHcVSRHee4q4q4ikuBqPfDFSR2+LRNghQB/gxX3xFgwsJ0Vdlw6TUREEnBXEXcVSYpbDZTO0N99sgJ/ytyPovNVAIDbuodg8d29W7x12h13f7Bfqb8td+xXtsp4Wr8ClB+zXPks7ipyIu4qkhaHmiiVoW/s2h5fPnkrlm0/iqXfFGH7kXLcsSQHvx/VE1OHxbR467Q77v5gv1J/W+7Yr2yV8bR+BXBXkTPaaQqTc4ma4O2lwxO/6Y4vnxyOhJhgVNea8ecvDmDiW9/hQIlJ6fCIiDSLExeiFlzfOv3XCf0Q4OOFPacv4W5unSYiUgwnLkQ26HQCHjR2wubf34YxfcNRz63TRESK4cSFSKSwQF8smzII7/zv1unrW6efXbMHFdXutX2TiMhdeWxyLs9xsS8ONVH7mQi392iH+LnD8NqmI/j4h1NY879bp++KEJBUW+v0OHiOizQ8x0We8jzHpTG1j1lyPovnuMiI57iQEn5963SfNhbcy1uniYjsxnNceI6LpLjVwN3ORKipt2DZtiK8nV0Ms1WAn7ceTyd1wxRjdItbp6XGwXNcpOE5LvKU5zkujbnbmOXIs3iOixPxHBdpcaiJu5yJYDAATyb1QEBFETZVhGDXiQr89cvD+GJvGf42sT96dxB36zTPcXENnuMiT3me49KYu4xZcjyL57gQeYBwP+Cj6fH/v3X6VAXGLfkW/+DWaSIi2XDiQiSjprZOv3V96/RRbp0mInIUJy5ETvDLrdNhgT7Xtk4vv7512v6dR0REdA0nLkRONLpvODal3oaHhkQDANbkn0bS69vx+Z4SeFhePBGRS3DiQuRkgb4G/GV8f6x9bCi6hfqjvLIWT3zyIx55bxfOVFxROjwiIrfCiQuRiwyOCcb6J27BU0nd4a3X4ZtD5zDy9e1YlXsCFr58ISISRbUTl+rqanTq1AnPPPOM0qEQycbHS4+nknrgyydvweBObVFda8ZfvzyMtH16HCq9rHR4RESqp9qJy1//+lcMGTJE6TCInKJbaABWzxqKv4zvB38fL5yoFDDh7e+5dZqIyAZVTlyOHDmCQ4cOYezYsUqHQuQ0Op2Ah4Z0wldPDMOAYAu3ThMRiSD7xCU7Oxvjxo1DREQEBEFAZmZmozLp6emIiYmBr68vjEYj8vLybvj8mWeewcsvvyx3aESqFB7oi0d6WvDW/XE3bJ2et5Zbp4mIfk32I/+rqqoQGxuL6dOnY+LEiY0+z8jIQGpqKpYtWwaj0Yi0tDSMHj0ahw8fRmhoKLKystCjRw/06NEDOTk5NturqalBTU1Nw88mkwkAb4e2Nw410dJNq7+sM6J7WwyZOwyvbjqCj/NOY/Wu09hy8Byev7MX7ugXBkEQ7GqLt0PfiLdDy1Oet0M3pqUxy+NvhxYEAevWrcP48eMbfmc0GhEfH4+lS5cCACwWC6KiojB37lzMnz8fCxYswIcffgi9Xo/KykrU1dXh97//PRYuXNhkG4sWLcLixYsb/Z63Q5M7O2YCPj2mRxlvnSYijVDF7dC/nrjU1tbCz88Pa9euvWEyk5ycjIqKCmRlZd1Qf9WqVdi3bx9effXVZtto6o1LVFQUzp492+QXr6+vx9atW5GYmAgvL/EvnOytJ7a8mHK2yrT0udTvqyRXxSxnO448q7m6tfUWLP/uBN7ZcQL1Fiv8vPV4MrEL7hsYhuzt2xzuW+xX6m/LGf3K0fKOjlme1q8AbY1ZzuxX69evV9/EpaSkBJGRkcjJycHQoUMbys2bNw/bt2/Hzp07b6gvZuJyXXp6OtLT02E2m1FYWMg3LuQxSquvvX0pvnzt7Usnfyt+18WMiNYKB0ZEJCOxb1xUPaWdOnWq6LIpKSlISUmByWRCUFAQRo0a1eQXr6urw6ZNmzBy5Ei7rti2t57Y8mLK2SrT0udSv6+SXBWznO048iwxdadarPh012m8svEITlTW45W9esy4OQZP3N4NPga9pOeyX6m/LWf3KynlHR2zPK1fAdoas5zZr3696tIcl26HDgkJgV6vR1lZ2Q2/LysrQ3h4uCtDIXIrOp2ABxKi8NUTw5DUKwQWq4B3vz2Bu9JzsbP4otLhERG5jCLJuQkJCViyZAmAa8m50dHRmDNnDubPny+5LS4VkZbsuSBgbbEOprpry0dDQi34bScL/FT9DpWIqHmKLRVVVlaiqKio4efi4mIUFBQgODgY0dHRSE1NRXJyMgYPHoyEhASkpaWhqqoK06ZNc6hdLhV5zqtXLb12lVK3rq4O2LQJs8bfirStxfjkh9P4/pwORdW+eP6OXhj7v63TXCq6EZeK5CnPpaLGtDRmqWGpSPaJy65du5CYmNjwc2pqKoBrO4dWrVqFyZMn4/z581i4cCFKS0sRFxeHDRs2ICwsTO5QiDxagK8XXry7D+6O7YA/ZR7AsfIqPLn6J2TuCcGiu3qjfWu+fiEiz+PUpSJX4lIRaVm9Bdh0RodNZwSYrQJ8dFbcGW3B8HArdILS0RER2aaKc1yUcH2pqLy8nEtFEuJWAy29dpVSt6XyR85V4vmsA8g/WQHg2tbpNx5MQN+Obe1ul/1K+bbU0q/sLcelIvW3o8a/C7OyskRNXFR5ySIRSdM91B8fPxKPxeN6w99HjxOVAiYt34XXNx1BDW+dJiIP4DFvXLhURHSjihrgs+M6/HTx2r9PQn2t+F1XM7o2/w8ZIiLFcKmIS0WS4lYDLb12lVLX3r6FjnH4y1eFOF957abpBxOi8PuR3eGrt7JfqbwtNfcrLhX9Py2NWWpYKvLYbQcGg6HFPyRbn8tVT2x5MeUc+U5Sv6+SXBWznO048ixn9a07BkRgRO8IvPTlQWTsOoWP8k7hm8PnsXhcb1HPYb9Svi019itHxyxP61eAtsYsZ/UrMTx24lJXV9fkNdlir9h2tJ7Y8mLK2SrT0udSv6+SXBWznO048ixX9C0/gwF/+W1v3NEvFM9lHcCpn6/g0Q9/xKAQHQZVVCGsTeOLj9ivlG9L7f1KShlP61eAtsYsZ/crMTxmqYg5LkTi1JqBL0/psO2sACsEtPayYmKMBYNCrBC4dZqIFMIcF+a4SIpbDbS0Xiylrlx9a/fxC3jqk104W31ttnJbjxD8+e4+6BDka7M++5Vr2nLHfmWrjKf1K0BbYxZzXJyIOS7S4lATLa0XS6nraN+6KaYdnulvxin/Xnhr2zFsLyzH2De/wx/G9sJDxk6i2mG/ck1b7tSvxJbxtH4FaGvMYo6LEzDHxb441ERL68VS6srVt+rq6uClAx69OQqj+4TiT1kHsPtkBRZm7cfnBWew6M4ezdZnv3JNW+7Yr2yV8bTxCtDWmMUcFxkxx4XIMRYr8G2pgC9O6lBjEaAXrBjT0YLfRFih51GVRORkzHFhjoukuNVAS+vFUurK1bea+7yk4goWfn4Q24+UAwB6hrXGyxP6oX9kkOSY1YA5LvKUZ45LY1oas5jj4kTMcZEWh5poab1YSl25+tavP+/U3oBV0xPw2a6TeCFrLw6XVWHSOzsxY3gXPJ3UA6289ZJjVgPmuMhTnjkujWlpzFIyx4UvgImoEUEQ8Nu4CCyIM+Ou/uGwWIF3s49hzBvZyDlarnR4RKRhnLgQUbMCDMA/7xuAfyUPRnigL05cqMYDy3fiT5n7UV2vdHREpEUeu1TEXUX2xaEmWsrQl1JXzl1FYj+/tVswvpw7DK9uKsTHeaexOv8MNhj0aN25BGP7R4iKW2ncVSRPee4qakxLYxZ3FcmIu4qIXKPIBGQc1ePc1WsH18UFW3BPZwsCvRUOjIjcGncVcVeRpLjVQEsZ+lLqOntXkZjPK69cxbPvbcPWEj3MViuCWnlhwZiemDgwAoJK7w3griJ5ynNXUWNaGrO4q8iJuKtIWhxqoqUMfSl1nbWrSMzn/gDGRVsw5+5heO7zA9h3xoT56/Zj/b4yvDShP6KC1fu2k7uK5CnPXUWNaWnM4q4iInJLfSMCkfn4zZg/thd8vHTYcaQco/6ZjX99WwyzxaNe5hKRSnDiQkQO8dLr8NhtXbHhqVth7ByMK3Vm/PmLA7h3WQ6KzlUqHR4ReRhOXIhIFp1DWuOTmUPw1wn94O/jhd0nK3DHmzvw9rajqDdblA6PiDwEJy5EJBudTsCDxk7Y+PStuK1He9TWW/D3DYcw8e0cHCo1KR0eEXkATlyISHYRbVph1bR4vHpvLAJ9vfDT6UsYt+RbvLH5COr49oWIHOCxu4p4AJ19caiJlg5zklJXiQPopMbw2wFhGBIThBf+exBbDp3HPzcX4qt9Z/G3CX3RN6L57Y7OwAPo5CnPA+ga09KYxQPoZMQD6IjUy2oFdl8Q8FmxDlX1AnSwIinSitEdLfDie18iAg+g4wF0EuNWAy0d5iSlrhoOoJP6fS9U1mDxF4fw1f4yAEC39q3x8oS+iItqI/oZUvEAOnnK8wC6xrQ0ZvEAOifiAXTS4lATLR3mJKWukgfQ2RvDdeFtDXh7ymB8tfcsns/ah6LzVZi8PA8zhndB6sge8DXoRT9LKh5AJ095HkDXmJbGLB5AR0SaMrZ/B2x6+jaMj4uAxQq8m30MY9/YgR+OX1Q6NCJSOU5ciEgRbVt7I+13A7Hi4cEIC/RBcXkV7nsnF4s+34/q2nqlwyMileLEhYgUldQnDBufvg33De4IqxVYlXMco9OykVNUrnRoRKRCnLgQkeKCWhnwj0mxeH96AiLbtMKpi1fwwIqd+OO6vbh81b22xhKRc3HiQkSqcWuP9tjw1HA8NCQaAPDxzpMY/c9sbC88r3BkRKQWnLgQkaoE+Brwl/H98fFMI6KD/VBy6SqS/52HBf/5iW9fiEh9E5eKigoMHjwYcXFx6NevH5YvX650SESkgGFdQ7DhqeGYOiwGAPBJ3imMSduB75j7QqRpqpu4BAQEIDs7GwUFBdi5cydeeuklXLhwQemwiEgBft5eWHR3X3z66BBEBbfCmYoreHDFTjyfuQ9VNdx5RKRFqpu46PX6hqP6a2pqYLVa4WGH+xKRnYZ0aYcNT97akPvywfcnMPaNHcgr5rkvRFoj+8QlOzsb48aNQ0REBARBQGZmZqMy6enpiImJga+vL4xGI/Ly8m74vKKiArGxsejYsSOeffZZhISEyB0mEbmZ1j5e+Mv4/vjwESMignxx8mI1Jr+bixf/ewBX68xKh0dELiL7kf9VVVWIjY3F9OnTMXHixEafZ2RkIDU1FcuWLYPRaERaWhpGjx6Nw4cPIzQ0FADQpk0b7NmzB2VlZZg4cSImTZqEsLCwJturqalBTU1Nw88mkwkAb4e2Nw410dJNq1LqutPt0M5gjAnCF3OG4uUNhViTfwb//q4YWw+V4e/39MPAFu484u3Q8pTn7dCNaWnM8vjboQVBwLp16zB+/PiG3xmNRsTHx2Pp0qUAAIvFgqioKMydOxfz589v9IzHH38ct99+OyZNmtRkG4sWLcLixYsb/Z63QxN5vv0/C8g4qsOlOgECrLg9woqxURYYVLcITkS2iL0d2qWXLNbW1iI/Px8LFixo+J1Op0NSUhJyc3MBAGVlZfDz80NAQAAuXbqE7OxszJ49u9lnLliwAKmpqQ0/m0wmREVFITExsckvXl9fj61btyIxMRFeXuK/vr31xJYXU85WmZY+l/p9leSqmOVsx5FnKdW3PKFfjQQw7UodXv76CD7/qQxbSgQcr/PHy7/tjX4RN/7v35Uxy9WWO/YrW2XcoV/ZS0tjljP71fr160XF4NI3LiUlJYiMjEROTg6GDh3aUG7evHnYvn07du7ciby8PDz66KMNSbkpKSmYNWuWzbbS09ORnp4Os9mMwsJCvnEh0pi9FwV8ekyHyjoBOliRFGnF6I4WePHtC5FbUOUbFzESEhJQUFBgd72UlBSkpKTAZDIhKCiIb1wkxq0GWvrXi5S6fOPStJEAplXX4i9fHcFX+89h4xkBJ+oD8fJve6NXuD/fuMhUnm9cGtPSmKW5Ny61tbXw8/PD2rVrb8h7SU5ORkVFBbKysiS3xTcuRHTdjxcErDmmQ1W9AJ1gxZiOFiRFWqEXlI6MiJqjyjcu3t7eGDRoELZs2dIwcbFYLNiyZQvmzJnj0LN//cZl1KhRTX7xuro6bNq0CSNHjoTBYBD9fHvriS0vppytMi19LvX7KslVMcvZjiPPUqpveXK/ugPArMoaLPz8IDYdPIcvT+lxxhKIce0v4oG7nR+zXH8+7tivbJVx537VHC2NWc7sV2JfXsg+camsrERRUVHDz8XFxSgoKEBwcDCio6ORmpqK5ORkDB48GAkJCUhLS0NVVRWmTZsmdyhEpGEh/j5Ivz8Wn+85i8XrD2HPGRMOntVDH3kGDxo7QRD4+oXIHcm+VLRt2zYkJiY2+n1ycjJWrVoFAFi6dCleeeUVlJaWIi4uDm+++SaMRqND7XKpiIiac7EG+LhIhyOma5m6fdpY8LuuFgR5KxwYETUQu1Tk1BwXJVxfKiovL+dSkYS41UBLr12l1OVSkTQ1tbV4/oMtWH/KgFqzBW39DHjx7j4Y07fpwy0dwaUiLhWpvR01/l2YlZUlauLCjYJEpAk6QcCIDlasmTkYfToE4OfqOsz9dA+eXbsXl6+610mtRFrmMW9cuFRERGLVW4ANp3XYfEaAFQLaelvxYDcLugd5xHBI5Ja4VMSlIklxq4GWXrtKqculImmainn3yQo8+9lenLx4BQAwbVgn/D6pG3wMetnbcvVz1DpmeVq/ArQ1Zqlhqch9Tvixk8FgaPEPydbnctUTW15MOUe+k9TvqyRXxSxnO448S6m+peV+ZezaHl89eSv+sv4gPsk7iZU5J/Dd0Qt4/b449IsMkrUtpZ6j1jHL0/oVoK0xy1n9SgyPnbjU1fF2aHviUBMt3bQqpa5cfYv96hpvHfDiuF4Y0aMd/pS5H4VllZjw1neYm9gVjw7vDL3O/m3Tcv35uGO/slXG0/oVoK0xy9n9SgyPWSpijgsROaqyDsg4psNPF6/tW+gaYMVD3c0I9lE4MCINYI4Lc1wkxa0GWlovllKXOS7SiI3ZarViXUEJXvziEKpqzQjw9cKf7+6DO/uHy96WM5+j1jHL0/oVoK0xizkuTsQcF2lxqImW1oul1GWOizRiYp6cEIMhXdvjyU8LUHCqAk+t/gk7ii5i8W/7wt9H/LDJHBfmuKi9HbX9XSgGz3EhImpCp3atseaxoZh7ezfoBOCz3adx55s7UHCqQunQiDTNY9+4MDnXvjjUREuJblLqMjlXGqkxP5HYBUM7t8Uza/fixIVq3PN2Dp5I7IpZtzafuMvkXCbnqr0dtf5dKIbH5LgwOZeInKm6Hlh9TIcfLzBxl8gZmJzL5FxJcauBlhLdpNRlcq40csT868TdwP8l7t7xq8RdJucyOVft7ajx70Im5zI5V1IcaqKlRDcpdZmcK42jMU9OiIGxS3s8mVGAPacq8OTqn7Dj6EUsurtx4i6Tc5mcq/Z21PZ3oRhMziUislNMSGusfWwo5iR2gyAAa/OvJe7uYeIukdNx4kJEJIFBr8Mzo3vi05lDEBHkixMXqjFpWQ5W7DgGD1uBJ1IVj10q4q4i++JQEy1l6Eupy11F0jgr5puiAvF5ylD8MXM/Nh44h7+sP4gdhecwKoi7irirSL3tqPXvQjE8JjmXu4qISElWK/BtmYDM4zrUWwUEeVvxcHczujWfY0hEv8BdRdxVJCluNdBShr6UutxVJI2rYj5w1oQnM/bg+IUr0AnA3MSumH1bF8mXNbpbv7JVxtP6FaCtMYu7ipyIu4qkxaEmWsrQl1KXu4qkcXbMsdHtkDl7KGa+sxk/nNfhjW+O4ocTFUibHIfQQF9Jz3THfmWrjKf1K0BbYxZ3FREReZDWPl54qJsF/5jYD37eeuQcvYCxb+zAtsPnlA6NyO1x4kJE5CQTBkbgv3NvQe8OgbhQVYupK3/Ay18dRJ3ZonRoRG6LExciIifq2t4f6x4fhoeHdgIAvLP9GO57JxenLlYrHBmRe+LEhYjIyXwNerz4235Y9tBNCPT1wo8nK3DXkm+x5WCZ0qERuR2PTc7lOS72xaEmWjoTQUpdnuMijStjbq6t3/QMQdbjQ/HU6p+w5/QlPPLeLswa3hlP/aYrvPSN/x3pjv3KVhlP61eAtsYsnuMiI57jQkTuot4CZJ3QIbv02mSlW+C1M1+CvBUOjEhBPMeF57hIilsNtHQmgpS6PMdFGlfGLLatr/aVYkHmflTVmBHi741/3jsAQ7oEyxKzWscsT+tXgLbGLJ7j4kQ8x0VaHGqipTMRpNTlOS7SuDJmW23dPTAK/Tq2xeMf7cah0stIXrULvx/VE7Nv6wrdLw6sc8d+ZauMp/UrQFtjFs9xISLSqC7t/bHu8Ztx76COsFiBV74+jEfe+wE/V9UqHRqRKnHiQkSksFbeerxybyz+cc8A+HjpsPXwedy15FvsOX1J6dCIVIcTFyIilbgvPgrrHr8ZMe38cKbiCu5fkYfsswI8LBWRyCGcuBARqUifiEB8PvcWjO0XjjqzFZ8d1+PpNXtRXVuvdGhEqsCJCxGRygT6GvDWgzfhT3f0hE6wYv3eUkxIz8Hx8iqlQyNSHCcuREQqJAgCpg7thDl9zGjv743DZZcxbum3+OYQT9slbVPdxOXUqVMYMWIE+vTpgwEDBmDNmjVKh0REpJiugcC62UMwqFNbXL5aj+mrdiFtcyEsFua9kDapbuLi5eWFtLQ0HDhwABs3bsRTTz2Fqiq+HiUi7QoL9MUnM4c0XNSYtvkIZr6/C5euuNfR+ERyUN3EpUOHDoiLiwMAhIeHIyQkBBcvXlQ2KCIihXl76fDib/vh1Xtj4eOlw5ZD53D30m9xqNSkdGhELiX7xCU7Oxvjxo1DREQEBEFAZmZmozLp6emIiYmBr68vjEYj8vLymnxWfn4+zGYzoqKi5A6TiMgtTRrUEZ/NHobINq1w4kI1JqTn4PM9JUqHReQysh/5X1VVhdjYWEyfPh0TJ05s9HlGRgZSU1OxbNkyGI1GpKWlYfTo0Th8+DBCQ0Mbyl28eBEPP/wwli9f3mJ7NTU1qKmpafjZZLr2rw/eDm1fHGqipZtWpdTl7dDSqOF2aLme0zPUD+tmG/H06r347ugFPPHJjyg4cRHPjurecMu0WscsT+tXgLbGLI+/HVoQBKxbtw7jx49v+J3RaER8fDyWLl0KALBYLIiKisLcuXMxf/58ANcmIyNHjsTMmTMxZcqUFttYtGgRFi9e3Oj3vB2aiDydxQqsP6nD5pLrt0xbMLWHBQHud80PkejboV16yWJtbS3y8/OxYMGCht/pdDokJSUhNzcXAGC1WjF16lTcfvvtNictALBgwQKkpqY2/GwymRAVFYVRo0bxdmgJcauBlm5alVKXt0NLo8bboeV4zl0Avt5fhj/8Zx+KTMBbR/zw9oNx6B7SSpVjlqf1K0BbY5azb4cWw6XJueXl5TCbzQgLC7vh92FhYSgtLQUAfPfdd8jIyEBmZibi4uIQFxeHvXv3NvtMHx8fBAYG3vB/RERaMrpvGNbOMiKmnR9KLl3F5OV5+Gr/OaXDInIKly4VlZSUIDIyEjk5ORg6dGhDuXnz5mH79u3YuXOn5LbS09ORnp4Os9mMwsJCLhURkeZU1wPvFepw6NK1f5OOirRgbJQFOkHhwIhEUOVSUUhICPR6PcrKbjz5saysDOHh4Q49OyUlBSkpKTCZTAgKCkJiYmKTX7y+vh5bt25FYmIivLzEf31764ktL6acrTItfS71+yrJVTHL2Y4jz1Kqb7Ffqb8tKc+5a7QVr285ipW5p7DxjA41fiF4ZWJf+Pu0XN9VY5an9StAW2OWM8er9evXi4pBkeTchIQELFmyBMC15Nzo6GjMmTOnITlXCr5xISL6fz+cF/DpUR3qrQLCW1kxo6cZ7VspHRVR8xR741JZWYmioqKGn4uLi1FQUIDg4GBER0cjNTUVycnJGDx4MBISEpCWloaqqipMmzbNoXZ//caFybnS4lYDLSW6SanL5FxpPDU5tzkj6+oQ+p9N+PCEH0ov12LJYV+k3ReLW7q1c6gtJuc2pqUxyyOTc3ft2oWBAwdi4MCBAIDU1FQMHDgQCxcuBABMnjwZr776KhYuXIi4uDgUFBRgw4YNjRJ2iYjIMZ0CgDUzBiMuKgiXrtTjkffzsTLnBJz4op3I6Zy6VORKXCoiImpanQVYfUyHvPPX/q2a0N6CyV0s8FLdpS+kZWKXijxm4nLd9aWi8vJyLhVJiFsNtPTaVUpdLhVJo7Wlol/XtVqtWJV7En/bcBgWKxAXFYS37o9D+wAfu9riUlFjWhqznL1UJGbiwvk2EZEGCIKAacM64d/JgxDUygsFpy7hnnd24uDZy0qHRmQXj3njwqUiIiJxzl8Blh/Wo+yKAG+dFQ93t6B/sEf8VUBujEtFXCqSFLcaaOm1q5S6XCqSRutLRb9mulKHJzJ+wndHL0AQgNTbuyCqqhCjRnGpyF5aGrPUsFTkPif82MlgMLT4h2Trc7nqiS0vppwj30nq91WSq2KWsx1HnqVU32K/Un9bzuhX7QwGrJqegBf/ewAffH8Cr205BmN7HUYKepeMWZ7WrwBtjVnOGq/E8NiJS11dXZPXZNt7JbfUemLLiylnq0xLn0v9vkpyVcxytuPIs5TqW+xX6m/LFf1q4Z090bldK/zly0PYeV6H5JU/IP2BgQhu7S35uVL7jjv2K0BbY5azxysxPGapiDkuRETSHfxZwKojOlw1C2jnY8WjvcwI5xBKLsQcF+a4SIpbDbS0XiylLnNcpGGOi+3yH2Rtwgcn/HG64ioCfL3wxuQBGN4txO7nMsdF/e2o8e9C5rgwx0VSHGqipfViKXWZ4yINc1yaF+4HrH1sCOZ+ugc/HP8ZMz/4EQvv6oPkYTGSnsscF/W3o7a/C8Xw2IkLc1zsi0NNtLReLKUuc1ykYY6LuPKB3gJWJg/C858fwLofS/DC5/txpMyEP43tCS+9jjkuTdDSmMUcFxkxx4WISD5WK7ClRMAXJ3WwQkCvIAum9bDA12P/uUtKY44Lc1wkxa0GWlovllKXOS7SMMdFWvmNB8rw+7V7cbXOgl5h/njr/v7Yu3MHc1x+QUtjFnNcnIg5LtLiUBMtrRdLqcscF2mY42Jf+TtjOyKqnT+mr9qFQ2WVuP9fu/FwZ+a4NEVLY5aSOS68q4iIiFo0oGMbZKYMQ48wf5RdrsEb+/TYVnhe6bBIozhxISIimzq29cOax4ZhWJdg1FoEzPrwR3yQe1zpsEiDPHapiLuK7ItDTbSUoS+lLncVScNdRY6X9/MC3rq/P2a9uxU7z+vwfNZ+HC+vxLxRPaDTCaKe5Wn9CtDWmMVdRTLiriIiItewWoFNZwSsP6UHAAwItmBKNwu89QoHRm6Nu4q4q0hS3GqgpQx9KXW5q0ga7iqSp/wvy311oBzz1+1DndmK2I5BWPZgHEL8fbiryA3aUePfhdxVxF1FkuJQEy1l6Eupy11F0nBXkTzlDQYD7hkcjah2/nj0g13Yc/oS7n03D6umxaNTW1+bz/K0fgVoa8ziriIiInJLCZ2D8Z/Zw9CpnR9O/3wFE9/Kwc7ii0qHRR6MExciInJIl/b++M/sYbgpug1MV+sx7b18/HBesF2RSAJOXIiIyGHt/H3w8cwhuLN/B9SZrfiwSI9l24/Bw9IoSQU4cSEiIln4GvRYcv9APHJzJwDAa5uL8HzWPpgtnLyQfDhxISIi2eh0AuaP6YmJMWYIAvDh9ycx64N8XKk1Kx0aeQiP3VXEA+jsi0NNtHSYk5S6PIBOGh5AJ095sWPWbR2suDW+L+atO4jNB8tw//JcvPPgQAR4C83Wd8d+BWhrzOIBdDLiAXREROpz1ASsOKRHtVlAe18rHuttRoiv0lGRGvEAOh5AJyluNdDSYU5S6vIAOml4AJ085aWMWUXnKjHjg904U3EV7VobkNzlCmZM8Ix+BWhrzOIBdE7EA+ikxaEmWjrMSUpdHkAnDQ+gk6e8PWNW78i2WPf4zZi68gccOGvC0v169I6rwMi+EbLErBZaGrN4AB0REXm00EBfZMwagpu7tkOtRcBjHxUg44eTSodFbogTFyIicokAXwOWTxmI+PYWmC1W/OGzvfjnpkKe9UJ24cSFiIhcxqDX4cGuFsy+rTMA4I0tRzD/s72oM1sUjozcBScuRETkUoIApCZ1x18n9INOADJ2ncLM93ehqqZe6dDIDXDiQkREinjQ2AnvThkMX4MO2w6fx0P/3gVTrdJRkdpx4kJERIpJ6hOGT2YOQXBrb+wrMeGN/XqcvFitdFikYqqcuEyYMAFt27bFpEmTlA6FiIicbGB0W3w2exg6tm2F8qsCJi/Pw74zl5QOi1RKlROXJ598Eu+//77SYRARkYt0DmmNjJkJiPSzoryyFr9793vkFJUrHRapkConLiNGjEBAQIDSYRARkQuFBvhgbl8zjJ3borKmHlNX/oAvfipROixSGdknLtnZ2Rg3bhwiIiIgCAIyMzMblUlPT0dMTAx8fX1hNBqRl5cndxhEROSGWnkB/5pyE+7oH45aswVzP/kR7+ceVzosUhHZj/yvqqpCbGwspk+fjokTJzb6PCMjA6mpqVi2bBmMRiPS0tIwevRoHD58GKGhoXa3V1NTg5qamoafTSYTAN4ObW8caqKlm1al1OXt0NLwdmh5yjs6Zon5TAcLXp/UH21bGfBR3ikszNqP0ooreOo3XSEIgohv41paGrM8/nZoQRCwbt06jB8/vuF3RqMR8fHxWLp0KQDAYrEgKioKc+fOxfz58xvKbdu2DUuXLsXatWtbbGPRokVYvHhxo9/zdmgiIvdmtQIbzwj48pQeADA01IJ7u1igV9/chWQg9nZol16yWFtbi/z8fCxYsKDhdzqdDklJScjNzZX0zAULFiA1NbXhZ5PJhKioKCQmJjb5xevr67F161YkJibCy0v817e3ntjyYsrZKtPS51K/r5JcFbOc7TjyLKX6FvuV+ttyx35lq4y9n40CkJBfghe/PIzcczr4tgnFa/f0ga9Bb/M7uYqWxixn9qv169eLisGlb1xKSkoQGRmJnJwcDB06tKHcvHnzsH37duzcuRMAkJSUhD179qCqqgrBwcFYs2bNDeWbkp6ejvT0dJjNZhQWFvKNCxGRB/npooD3CnWotwroEmDFzF5m+LnHnJlEUuUbF7E2b95sd52UlBSkpKTAZDIhKCiIb1wkxq0GWvrXi5S6fOMiDd+4yFPe1W9crhsJ4LYTFUj5dC+OXa7Hv08E4d0HYhEW6GPzuzmblsYszb1xqa2thZ+fH9auXXtD3ktycjIqKiqQlZUluS2+cSEi8nwlVcCyg3pcqhPQ1tuK2X3MCGuldFQkB1W+cfH29sagQYOwZcuWhomLxWLBli1bMGfOHIee/es3LqNGjWryi9fV1WHTpk0YOXIkDAaD6OfbW09seTHlbJVp6XOp31dJropZznYceZZSfYv9Sv1tuWO/slVGjn41tuIKpr+Xj2Pl1XjrcCu8O2UgBka1sfkdnUVLY5Yz+5XYlxeyn+NSWVmJgoICFBQUAACKi4tRUFCAkydPAgBSU1OxfPlyvPfeezh48CBmz56NqqoqTJs2Te5QiIjIA0W2aYVPZiRgQMdAVFypw8Mrd2Fb4XmlwyIXkX2paNu2bUhMTGz0++TkZKxatQoAsHTpUrzyyisoLS1FXFwc3nzzTRiNRofa5VIREZG21JiBlYU6HKzQQQcr7u9qQUKo07IfyMnELhU5NcdFCdeXisrLy7lUJCFuNdDSa1cpdblUJA2XiuQpr4alohvqmC3447r9yNxzFgAwf0wPPHJzjKi6ctHSmOXspSIxExdV3lVEREQkhkGvw98n9sMjN3cCAPxtQyFe3XgEHvZvcvoFj3njwqUiIiJt23JGwOcn//+U3fu6WKDjKbtug0tFXCqSFLcaaOm1q5S6XCqShktF8pRX21LRr63edRrPf34AFiswuk8oXrt3AHy8nLu4oKUxSw1LRe5xcpQEBoOhxT8kW5/LVU9seTHlHPlOUr+vklwVs5ztOPIspfoW+5X623LHfmWrjLP61YNDO6NdgC+e+KQAXx84h8qPfsQ7UwbD38f5f91pacxyVr8Sw2MnLnV1vB3anjjUREs3rUqpK1ffYr9Sf1vu2K9slXFFv/pNzxCseHggZn9UgO+KLuD+d3OxYspNCG7t7dBzm6OlMcvZ/UoMj1kqYo4LERH90snKa6fsVtULCPW9dspusPI3BFAzmOPCHBdJcauBltaLpdRljos0zHGRp7zac1x+7ej5Kkx7Lx9nL11FeKAPViYPQrdQf1mefZ2WxizmuDgRc1ykxaEmWlovllKXOS7SMMdFnvJqzXH5tV4RbfDZ7GGY8q+dOHq+Cg/86wesmpaAWCdcEaClMUvJHBee40JERB4tok0rrHlsGGI7BuHn6jrcv/x7fHukXOmwSCKPfePC5Fz74lATLSW6SanL5FxpmJwrT3l3SM5tSoC3gFVTByHlkwLkHL2Iaavy8Nqk/hjbL9zhZ2tpzGJyroyYnEtERLbUW4APjuhQcFEHAVbc18WCYWEe8deg22NyLpNzJcWtBlpKdJNSl8m50jA5V57y7pac2xSzxYoX/nsQGbtOAwB+n9QNs27tDEGQdsyulsYsJuc6EZNzpcWhJlpKdJNSl8m50jA5V57y7pKc2+TzAfztngEICfBB+tajeG1zES5dNeOPd/SGzoE7ArQ0ZjE5l4iIyIUEQcCzo3vhuTt7AwBWfFuMZ9f+hHqzReHIyBZOXIiISLNmDO+C1+6NhV4n4LPdp/HYh7txtc6sdFjUAo9dKuKuIvviUBMtZehLqctdRdJwV5E85d11V1FL7h4QhtbesXgy4ydsPliGKf/aiXcejEOAr7ilDS2NWdxVJCPuKiIiIkcUmYDlh/S4ahYQ6WfFY73NCHTO9UbUBO4q4q4iSXGrgZYy9KXU5a4iabirSJ7ynrCrqCUHzpow/b3duFBVi07Bflg1dRA6tm3VYh0tjVncVeRE3FUkLQ410VKGvpS63FUkDXcVyVPenXcVtSQ2uh0+mz0MD/1rJ05crMb9K37AhzMS0C00wGZdLY1Z3FVERESkEjEhrbH2sWHoHuqPUtNV3LssFz+drlA6LPofTlyIiIh+JTzIF6tnDW243+iB5TuRe/SC0mEROHEhIiJqUtvW3vho5hAM7dIOlTX1SF6Zh80HypQOS/M4cSEiImqGv48XVk6Lx8g+Yaitt2DWh/nI/PGM0mFpmscm5/IcF/viUBMtnYkgpS7PcZGG57jIU94Tz3GxRQ/gzfv6Y8E6HTL3nMVTGQX4ueoqHjJGA9DWmMVzXGTEc1yIiMiZLFZg3XEdskuvLVbcGWXGyEgrJN7NSL/Cc1x4joukuNVAS2ciSKnLc1yk4Tku8pT39HNcbLFarXjzm6NYuu0YAOCRmzsh9fbO2Lx5sybGLJ7j4kQ8x0VaHGqipTMRpNTlOS7S8BwXecp76jkuYjwzpjfa+vviz18cwL++OwHT1XoMM2hrzOI5LkRERG7kkVs64x+TBkAnAGvyz2BVoQ419bxZ2hU4cSEiIpLgvsFReOvBm2DQC9hzUYfHPvoR1bX1Sofl8ThxISIikmhMvw5YPuUmeOus+LboAh5asROXqtW1K8rTcOJCRETkgJu7tkNKHzOCWnlh98kKTH43F+cuX1U6LI/FiQsREZGDYgKAj6bHo32ADw6VXsZ9y3Jx6mK10mF5JE5ciIiIZNAzPABrHxuKqOBWOH6hGvcuy0XRuctKh+VxVDlx+eKLL9CzZ090794dK1asUDocIiIiUTq1a401s3iztDOpbuJSX1+P1NRUfPPNN/jxxx/xyiuv4MIF3shJRETugTdLO5fqJi55eXno27cvIiMj4e/vj7Fjx2Ljxo1Kh0VERCQab5Z2HtknLtnZ2Rg3bhwiIiIgCAIyMzMblUlPT0dMTAx8fX1hNBqRl5fX8FlJSQkiIyMbfo6MjMSZM7yJk4iI3AtvlnYO2Y/8r6qqQmxsLKZPn46JEyc2+jwjIwOpqalYtmwZjEYj0tLSMHr0aBw+fBihoaF2t1dTU4OampqGn00mEwDeDm1vHGqipZtWpdTl7dDS8HZoecpr8XZoW1qK+/rN0n/M1GFdQeObpeVqx1XP8vjboQVBwLp16zB+/PiG3xmNRsTHx2Pp0qUAAIvFgqioKMydOxfz589HTk4OXnnlFaxbtw4A8NRTTyEhIQEPPPBAk20sWrQIixcvbvR73g5NRERqwZulbVPF7dC/nrjU1tbCz88Pa9euvWEyk5ycjIqKCmRlZaG+vh69e/fGtm3bEBQUhEGDBiEnJwft2rVrso2m3rhERUXxdmiJcasBb4eWpzxvh74Rb4eWp7zWb4duiti4rVYrlmw9iiVb//9m6T+M7gFB5OxFDWOW5m6HLi8vh9lsRlhY2A2/DwsLw6FDh64F5OWF1157DYmJibBYLJg3b16zkxYA8PHxgY+Pj1PjJiIicpQgCHji9m4IbGXAX7883HCz9J/v7gO9jq9exHLpG5fribc5OTkYOnRoQ7l58+Zh+/bt2Llzp+S20tPTkZ6eDrPZjMLCQi4VERGRau08J+CTozpYISA22IKHu1vgpbp9vq4ldqnIpW9cQkJCoNfrUVZ245awsrIyhIeHO/TslJQUpKSkwGQyISgoCImJiU1+8fr6emzduhWJiYnw8hL/9e2tJ7a8mHK2yrT0udTvqyRXxSxnO448S6m+xX6l/rbcsV/ZKuNp/QqQFvdIAEMOncfvP9uPPRd1+M/5dnjjvv7w89bL2o7cz3Jmv1q/fr2oGBRJzk1ISMCSJUsAXEvOjY6Oxpw5czB//nzJbfGNCxERuZvDFQJWHNah1iKgc4AVj/Yyw8995myyUuyNS2VlJYqKihp+Li4uRkFBAYKDgxEdHY3U1FQkJydj8ODBSEhIQFpaGqqqqjBt2jSH2v31G5dRo0YxOVdC3GrA5Fx5yjM590ZMzpWnPJNzG3Mk7jsAjDhVgZkf7Ebx5Xq8d6otVibfhBD/xrmbahiznJ2cK4bsK2q7du3CwIEDMXDgQABAamoqBg4ciIULFwIAJk+ejFdffRULFy5EXFwcCgoKsGHDhkYJu0RERFowMKoNPpwejxB/bxwqvYz7V/yAkoorSoelWk5dKnIlLhUREZE7O38FeOugHhdrBLTxtuLxPmaEtVI6KtdRxTkuSri+VMRzXKTFrQZcKpKnPJeKbsSlInnKc6moMTnjPnvpKqauysex8ioEtzZgZfIg9OkQKHs7avy7UOw5LhrffEVERKQeHYJ88fGMePSNCMDFqjo89O9dyD/xs9JhqYrHvHHhUhEREXmKK/XA8kN6HL0swFtnxSM9LejVxiP+um4Wl4q4VCQpbjXgUpE85blUdCMuFclTnktFjTkr7iu1Zsz9dA+2HymHQS/glYl9IJzeo/mlIo/dLW4wGFr8Q7L1uVz1xJYXU86R7yT1+yrJVTHL2Y4jz1Kqb7Ffqb8td+xXtsp4Wr8C5I/bYDBgeXI8nl5dgPU/nUXq2v2Y3EXAHSoYs5zVr8Tw2IlLXV1dk9dki71i29F6YsuLKWerTEufS/2+SnJVzHK248izlOpb7Ffqb8sd+5WtMp7WrwDnxi0AeO2efmht0GF1/hl8clSPmG+PYfotXRx6rlr/LhTDY5aKmONCRESeymoFsk7osPXstT01YzqaMaajFSIvlnYLzHFhjoukuNWAOS7ylGeOy42Y4yJPeea4NOaquGtra/Hsqi348tS1+4ymDeuEBWN6QJAwe1Hj34XMcWGOi6Q41IQ5LvKUZ47LjZjjIk955rg05oq4R3e0YtCAXvjz+kNYmXMCVbVmvDxxAPQ6aa9e1PZ3oRgeO3Fhjot9cagJc1zkKc8clxsxx0We8sxxaczVY9b9gzogwMcL89ftw+pdp2G6UodXJ/WHj5f4o9nU+nehGB6zVMQcFyIi0pI9FwS8d0QHs1VAryALpve0wEevdFTSMceFOS6S4lYD5rjIU545Ljdijos85Znj0piSY9Z3Ry9g9kc/4kqdBYOi2+DdhwYisJV8/70drcccFzswx0VaHGrCHBd5yjPH5UbMcZGnPHNcGlNizBrRKxwfzhiCaSvzkH+yAlNW5uP9RxIQ4u9j97OkxuDMdprCu4qIiIjc2KBObfHpo0MR4u+NA2dNuG9ZLkoqrigdltNw4kJEROTm+kQEYs1jwxDZphWOlVfh3mW5OHa+UumwnMJjl4q4q8i+ONSEu4rkKc9dRTfiriJ5ynNXUWNqGbM6BnnjkxnxmLpqF46VV+PeZblYmTwIvTsE2P0sqTFILc9dRdxVREREGnW5Dnj7gB5nqgW00lsxq7cZnRvPXVSHu4q4q0hS3GrAXUXylOeuohtxV5E85bmrqDE1jlmmK3V49MMfkX+yAq0MOrz1wEDc0q2dwzFzV5ETcVeRtDjUhLuK5CnPXUU34q4iecpzV1Fjahqz2hkM+GCGEbM+yMeOI+WY9eGPePP+gRjTL9zuZ0mNQY52msLkXCIiIg/k5+2FFcmDMbZfOGrNFjz+UT7W5p9WOiyHceJCRETkoXy89Fhy/0DcO6gjLFbgmTV7sOq7YqXDcojHLhURERER4KXX4e/3DECArwH//q4Yi/57ABXVtYhx0wxXvnEhIiLycDqdgOfv6o2nk3oAANK2FCHrhA7uuD+HExciIiINEAQBTyZ1x8K7+gAAtp7V4U9ZB2C2uNfkxWOXingAnX1xqIlaDnNy1bN4AJ1r8AA6ecrzALrG3G3MmmLsCB+dBc99fhBr8s+g8mo9Xp3UH95ett9l8AA6GfEAOiIiIvH2XBDw3hEdzFYBvdtYML2HBd565eLhAXQ8gE5S3GqgxsOcnPksHkDnGjyATp7yPICuMXces/y6DMITq/fiSp0Fgzu1wbsPDUSAr2P//aWU5wF04AF0UuNQEzUd5uSKZ/EAOtfgAXTylOcBdI2545g1olcYPpzhh6krf8CuExWYsjIf709PQDt/H1lj4AF0REREJItBnYLx6aND0K61N/aXmHDfO7k4e+mK0mE1ixMXIiIijesbEYTVjw1FRJAvjp6vwqS3c1FcXqV0WE3ixIWIiIjQtb0/1swehs4hrXGm4gruXZaLg2dNSofVCCcuREREBACIbNMKq2cNRe8OgSivrMHkd3Kx++TPSod1A05ciIiIqEH7AB98OnMIbopuA9PVejy0Yie+PVKudFgNOHEhIiKiGwT5GfDhDCOGdw9Bda0Z01f9gK/3lyodFgCVTlwmTJiAtm3bYtKkSUqHQkREpEl+3l5YkTwYY/qGo9ZsweMf7UZmQYnSYalz4vLkk0/i/fffVzoMIiIiTfPx0mPpAwMxaVBHmC1WPPvZPuwoFRSNSZUTlxEjRiAgIEDpMIiIiDTPS6/DP+4ZgKnDYgAAa4v1WPHtccXisXvikp2djXHjxiEiIgKCICAzM7NRmfT0dMTExMDX1xdGoxF5eXlyxEpEREQK0OkEvDCuD+YmdoFBsGJAx+aP5Hc2u4/8r6qqQmxsLKZPn46JEyc2+jwjIwOpqalYtmwZjEYj0tLSMHr0aBw+fBihoaEAgLi4ONTX1zequ3HjRkRERNgVT01NDWpqahp+Npmu7Tnn7dD2xaEm7nbTqqPP4u3QrsHboeUpz9uhG9PSmDV7eCe0rSjEwMgAUXVVdzu0IAhYt24dxo8f3/A7o9GI+Ph4LF26FABgsVgQFRWFuXPnYv78+aKfvW3bNixduhRr165tsdyiRYuwePHiRr/n7dBERETuQ+zt0LJeslhbW4v8/HwsWLCg4Xc6nQ5JSUnIzc2Vs6kGCxYsQGpqasPPJpMJUVFRSExMbPKL19fXY+vWrUhMTISXl/ivb289seXFlLNVpqXPpX5fJbkqZjnbceRZSvUt9iv1t+WO/cpWGU/rV4C2xixn9qv169eLikHWNy4lJSWIjIxETk4Ohg4d2lBu3rx52L59O3bu3CnquUlJSdizZw+qqqoQHByMNWvW3PC8pqSnpyM9PR1msxmFhYV840JERORGFHnjIpfNmzfbXSclJQUpKSkwmUwICgrCqFGjmvzidXV12LRpE0aOHGnXFdv21hNbXkw5W2Va+lzq91WSq2KWsx1HnqVU32K/Un9b7tivbJXxtH4FaGvMcma/ysrKEhWDrNuhQ0JCoNfrUVZWdsPvy8rKEB4eLmdTREREpEFOSc5NSEjAkiVLAFxLzo2OjsacOXPsSs61F5eKiIiI3JfTlooqKytRVFTU8HNxcTEKCgoQHByM6OhopKamIjk5GYMHD0ZCQgLS0tJQVVWFadOmSfsmIv16qYjJudLiVgMtJbpJqcvkXGmYnCtPeSbnNqalMcstk3O3bduGxMTERr9PTk7GqlWrAABLly7FK6+8gtLSUsTFxeHNN9+E0Wi0pxm78Y0LERGR+xL7xsWhpSI1uv7Gpby8nMm5EuJWAy0lukmpy+RcaZicK095Juc2pqUxy9nJuWImLqq8q4iIiIioKR7zxoVLRURERO6LS0VcKpIUtxpo6bWrlLpcKpKGS0XylOdSUWNaGrPUsFTkPmnbdjIYDC3+Idn6XK56YsuLKefId5L6fZXkqpjlbMeRZynVt9iv1N+WO/YrW2U8rV8B2hqznNWvxPDYiUtdHW+HticONdHSTatS6srVt9iv1N+WO/YrW2U8rV8B2hqznN2vxPCYpaLrOS719fU4cuQIVqxYwRwXIiIiN1FdXY0ZM2agoqICQUFBzZbzmInLdadPn0ZUVJTSYRAREZEEp06dQseOHZv93OMmLhaLBSUlJQgICIAgCE2WiY+Pxw8//GD3s+2tJ7a8mHK2yjT3uclkQlRUFE6dOtVispPaSP1vpGQ7jjxLqb7FfqX+ttyxX9kq42n9CtDWmOWsfjV48GB88803iIiIgE7X/GktHpfjotPpWpypAYBer5f0Pwp764ktL6acrTK2Pg8MDHSrgUDqfyMl23HkWUr1LfYr9bfljv3KVhlP61eAtsYsZ/UrLy8vm39/Axo9gC4lJcUl9cSWF1POVhmp30mtXPV95GzHkWcp1bfYr9Tfljv2K1tlPK1fAdoas5TsV4AHLhXRja6fa2NrXzyRPdivyBnYr0gMTb5x0RIfHx+88MIL8PHxUToU8iDsV+QM7FckBt+4EBERkdvgGxciIiJyG5y4EBERkdvgxIWIiIjcBicuRERE5DY4cSEiIiK3wYmLhk2YMAFt27bFpEmTlA6FPMSpU6cwYsQI9OnTBwMGDMCaNWuUDok8REVFBQYPHoy4uDj069cPy5cvVzokUgi3Q2vYtm3bcPnyZbz33ntYu3at0uGQBzh79izKysoQFxeH0tJSDBo0CIWFhWjdurXSoZGbM5vNqKmpgZ+fH6qqqtCvXz/s2rUL7dq1Uzo0cjG+cdGwESNGICAgQOkwyIN06NABcXFxAIDw8HCEhITg4sWLygZFHkGv18PPzw8AUFNTA6vVCv67W5s4cXFT2dnZGDduHCIiIiAIAjIzMxuVSU9PR0xMDHx9fWE0GpGXl+f6QMmtyNmv8vPzYTabERUV5eSoyR3I0bcqKioQGxuLjh074tlnn0VISIiLoic14cTFTVVVVSE2Nhbp6elNfp6RkYHU1FS88MIL2L17N2JjYzF69GicO3fOxZGSO5GrX128eBEPP/ww3n33XVeETW5Ajr7Vpk0b7NmzB8XFxfj4449RVlbmqvBJTazk9gBY161bd8PvEhISrCkpKQ0/m81ma0REhPXll1++odzWrVut99xzjyvCJDcjtV9dvXrVOnz4cOv777/vqlDJzTgyZl03e/Zs65o1a5wZJqkU37h4oNraWuTn5yMpKanhdzqdDklJScjNzVUwMnJnYvqV1WrF1KlTcfvtt2PKlClKhUpuRkzfKisrw+XLlwEAly5dQnZ2Nnr27KlIvKQsTlw8UHl5OcxmM8LCwm74fVhYGEpLSxt+TkpKwr333osvv/wSHTt25KSGWiSmX3333XfIyMhAZmYm4uLiEBcXh7179yoRLrkRMX3rxIkTGD58OGJjYzF8+HDMnTsX/fv3VyJcUpiX0gGQcjZv3qx0CORhbrnlFlgsFqXDIA+UkJCAgoICpcMgFeAbFw8UEhICvV7fKHGtrKwM4eHhCkVF7o79ipyFfYvswYmLB/L29sagQYOwZcuWht9ZLBZs2bIFQ4cOVTAycmfsV+Qs7FtkDy4VuanKykoUFRU1/FxcXIyCggIEBwcjOjoaqampSE5OxuDBg5GQkIC0tDRUVVVh2rRpCkZNasd+Rc7CvkWyUXpbE0mzdetWK4BG/5ecnNxQZsmSJdbo6Girt7e3NSEhwfr9998rFzC5BfYrchb2LZIL7yoiIiIit8EcFyIiInIbnLgQERGR2+DEhYiIiNwGJy5ERETkNjhxISIiIrfBiQsRERG5DU5ciIiIyG1w4kJERERugxMXIiIichucuBAREZHb4MSFiIiI3AYnLkREROQ2/g9XgEt2V6P2uAAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"result = virgo.model.run(FrequencyResponse(geomspace(5, 5000, 1000), [\"DARM_Fz\"], [\"B1p_56.I\"]))\n",
"\n",
"Figure = figure()\n",
"_ = Figure.gca().loglog(result.f, abs(result[\"B1p_56.I\", \"DARM_Fz\"]))\n",
"Figure.gca().grid(True, \"both\", \"both\")\n",
"show()"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "f40772e1-6e9d-462f-8a42-195cc772fb88",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([<Axes: title={'center': 'DARM TF'}, ylabel='Amplitude [$\\\\sqrt{W}$]'>,\n",
" <Axes: xlabel='f [Hz]', ylabel='Phase [deg]'>], dtype=object)"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAkIAAAHLCAYAAAAk8PeNAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAc75JREFUeJzt3XdcU+fiBvAng4S9CTOoqIATLQhVq6Ji0V61Vm2vHUrt+rW32kFtr7a22tt97bC2dt06OtTrtVY7tVYUte6FGxRB2QFEVphJzu8PIBpBZIQkkOf7+eQDOefkvC94ynn6riMSBEEAERERkRUSm7sCRERERObCIERERERWi0GIiIiIrBaDEBEREVktBiEiIiKyWgxCREREZLUYhIiIiMhqMQgRERGR1WIQIiIiIqvFIERERERWi0GIiExm9erVEIlE+petrS38/PwQGxuLZcuWoaysrNnPR0ZGQiQS4fPPP2/R+aVSKfz9/fHwww8jOzu70fHR0dEQiUTo3bt3k+f7888/9ef64YcfblqvhvPc6rV48WIAQPfu3W96TFVVVbO/AyIyLqm5K0BE1udf//oXevTogdraWuTl5SExMRHPPfccPvzwQ/z8888YOHBgo89cuHABhw8fRvfu3bFmzRo89dRTtzx/VVUVDhw4gNWrV+Ovv/7C6dOnYWtra3Csra0tUlNTcejQIURGRhrsW7NmDWxtbW8ZTl555RU89thj+veHDx/GsmXL8PLLL6NPnz767df/XIMGDcILL7zQ6FwymazZsojIuBiEiMjkJkyYgIiICP37BQsWYMeOHZg4cSImT56Mc+fOwc7OzuAz33//PRQKBT744ANMnz4dly5dQvfu3W95/sceewyenp5477338PPPP+O+++4zOLZnz57QaDRYt26dQRCqqqrCpk2b8Le//Q0bN25s9ucZN26cwXtbW1ssW7YM48aNQ3R0dJOf8ff3x0MPPdTseYmo47FrjIgswpgxY/Dqq6/i8uXL+P777xvtX7t2LaZPn46JEyfCxcUFa9eubfG5R4wYAQC4ePFik/vvv/9+rF+/HjqdTr/tl19+QUVFRaPgRERdC4MQEVmMmTNnAgC2bdtmsP3gwYNITU3F/fffD5lMhqlTp2LNmjUtPu+lS5cAAG5ubk3uf+CBB5Cbm4vExET9trVr12Ls2LFQKBSt+yFaqLa2FoWFhQavioqKDimLiG6OQYiILEZAQABcXFwatdx8//33UCqVGD58OABgxowZOHv2LJKSkpo8T0lJCQoLC5GVlYWNGzfi9ddfh1wux8SJE5s8vnfv3oiIiNC3MhUXF+P333/HAw88YLwf7gbbtm2Dl5eXwevf//53h5VHRE3jGCEisiiOjo4Gs8c0Gg3Wr1+PuLg4iEQiAHXdaAqFAmvWrMGgQYManSMmJsbgfffu3fH9998jICDgpuU+8MADeOONN/DZZ5/hhx9+gEQiwT333IOjR48a5we7QVRUFN58802DbUFBQR1SFhHdHIMQEVmU8vJyg+6obdu2oaCgAJGRkUhNTdVvHz16NNatW4f33nsPYrFh4/by5csRHByMkpISrFy5Ert374ZcLm+23BkzZmDevHnYsmUL1qxZg4kTJ8LJycm4P9x1PD09GwU2IjI9BiEishhZWVkoKSlBr1699NsaxgLdbNDyrl27MHr0aINtkZGR+lljU6ZMwR133IEHHngAKSkpcHR0bPI8vr6+iI6OxgcffIC9e/fecqYYEXUNHCNERBbju+++AwDExsYCANRqNX766Sf8/e9/x4YNGxq9fH19bzloWiKR4J133kFOTg4+/fTTZo994IEHsGfPHjg7O+Ouu+4yzg9FRBaNLUJEZBF27NiBN954Az169MCDDz4IANi0aRPUajWefvpp/RT4623btg0bNmzA8uXLm+36io6ORmRkJJYuXYrnnnuu0aKKDaZPn47MzEyEhIRwYUMiK8EWISIyuS1btuD777/H6tWr8d577yE2NhYxMTHw9vbGzz//rA8qa9asgYeHB4YNG9bkeSZPnozi4mL89ttvtyzzxRdfhEqlwurVq296jIuLCxYvXoz777+/TT8XEXU+bBEiIpN77bXXANQ9TsLd3R0DBgzA0qVLMXv2bP0A5fz8fGzfvh33338/JBJJk+cZO3Ys7O3t8f3332Pq1KnNljl16lT07NkT77//Ph5//PGbnpOIrItIEATB3JUgIiIiMgd2jREREZHVYhAiIiIiq8UgRERERFaLQYiIiIisFoMQERERWS0GISIiIrJaXEeoGTqdDjk5OXByctI/9ZqIiIgsmyAIKCsrg5+fX6OHMt+IQagZOTk5UCqV5q4GERERtUFmZiYCAgKaPYZBqBkNK9xmZmbC2dnZzLUhIiKiligtLYVSqdTfx5vDINSMhu4wZ2dnBiEiIqJOpiXDWjhYmoiIiKwWgxARERFZLXaNERF1ETqdDjU1NeauBpFJ2NjYQCKRtPs8DEJW5PClIrz52zlEB3vhkeE94GJvY+4qEZGR1NTUID09HTqdztxVITIZV1dX+Pj4tGuJGwYhK/Llros4kVmME5nFWPlXOpbcOxDj+/uau1pE1E6CICA3NxcSiQRKpfKW66YQdXaCIKCiogL5+fkAAF/ftt/LGISshEarw8G0IgBAgJsdsq5W4o1fzyG2X/uSNBGZn0ajQUVFBfz8/GBvb2/u6hCZhJ2dHQAgPz8fCoWizd1k/N8GK3EquwRl1Ro420qx9bmRkEnEyC6uRHqh2txVI6J20mq1AACZTGbmmhCZVkPwr62tbfM5GISsxL6LVwAAtwd5wFEuRXg3NwDAX6mF5qwWERkRW3fJ2hjjmmcQshL764PQsJ4eAIA7ensCAPZcYBAiIiLrxSBkBao1Why+VDc+aHivugA0oj4I7b94BbVazjIhIiLrxCBkBY5nFKNao4Onoxy9FI4AgH5+LnCzt0F5tQYnMovNW0EiIiIzYRCyAvvqxwEN6+mh70+ViEUYVt86tJvdY0RkBg8//DBEIhFEIhFsbGzg7e2NcePGYeXKlU2uhxQbGwuJRILDhw/f8lw9evTASy+9hKqqKoPjGo45cOCAwfbq6mp4eNT9jUxMTGxR/RvOdf3rjjvuaPkvoB0uXbpkUK5MJkOvXr3w5ptvQhAE/XFnzpzBtGnT0L17d4hEIixdutTodbn+d9/wGj9+fKPjfvvtN0RFRcHOzg5ubm6YMmVKs+eNjo7Gc889Z/T63ohBqIvT6QRsP1e3zkLD+KAGI+qD0F8XCkxeLyIiABg/fjxyc3Nx6dIlbNmyBaNHj8azzz6LiRMnQqPR6I/LyMjAvn37MGfOHKxcubLZc6WlpeGjjz7Cl19+iUWLFjU6TqlUYtWqVQbbNm3aBEdHx1bXf9WqVcjNzdW/fv7551afo0FWVpZBiGmJ7du3Izc3FxcuXMDrr7+Ot956y+D3U1FRgaCgILz77rvw8fFpc91upeF33/Bat26dwf6NGzdi5syZmD17Nk6cOIG9e/figQce6LD6tAaDUBe35uBlnM0thb1MgjF9FAb7GgZMJ2UWI7OowhzVI6IOIAgCKmo0Znm19kYul8vh4+MDf39/3HbbbXj55Zfx008/YcuWLVi9erX+uFWrVmHixIl46qmnsG7dOlRWVt70XEqlElOmTEFMTAz+/PPPRsfFxcXhv//9r8E5Vq5cibi4uFbVHbi2snHDy93dHQCwePHiJluMrv+ZbvTqq68iKCgIixYtQlpaWovK9/DwgI+PD7p164YHH3wQw4cPx7Fjx/T7hwwZgiVLlmDGjBmQy+VNniM6Ohpz5szBnDlz4OLiAk9PT7z66qut+rds+N03vNzc3PT7NBoNnn32WSxZsgRPPvkkgoOD0bdvX9x3330tPn9H4oKKXVh2cSXe3ZIMAPjn+FAonGwN9ge42WNEb0/suVCIZQkXsOTeMHNUk4iMrLJWi76v/WGWss/+Kxb2svbdWsaMGYOwsDD8+OOPeOyxxyAIAlatWoXly5cjNDQUvXr1wg8//ICZM2fe9BynT5/Gvn370K1bt0b7wsPD0b17d2zcuBEPPfQQMjIysHv3bixfvhxvvPFGu+reYN68eXjyySf179esWYPXXnsNERERN/3MsmXLsGHDBnz77bd48803MXz4cDz88MO499574eTkdMsyjxw5gqNHj2LWrFmtru8333yDRx99FIcOHcKRI0fwxBNPIDAwEI8//niLPp+YmAiFQgE3NzeMGTMGb775Jjw86nohjh07huzsbIjFYgwePBh5eXkYNGgQlixZgv79+7e6rsbGFqEu6lKhGvP+dwLqGi0iurlh5u2N/xgAQPy4YADAj8ezkVZQbsoqEhHdVGhoKC5dugSgrvunoqICsbGxAICHHnoIK1asaPSZX3/9FY6OjrC1tcWAAQOQn5+PF198scnzP/LII/oupNWrV+Ouu+6Cl5dXq+t5//33w9HRUf/avHkzAMDR0VHfOnLp0iUsXLgQq1atavbG7+TkhEceeQSJiYlIS0vDnXfeiffeew8+Pj546KGH8OeffzZqpRk2bBgcHR0hk8kwZMgQ3HfffW0KQkqlEh999BFCQkLw4IMPYu7cufjoo49a9Nnx48fj22+/RUJCAt577z3s2rULEyZM0C/02dC6tXjxYixcuBC//vor3NzcEB0djaKiolbX1djYItTF5BRX4p8bT+rXB5JJxHh32kCIxU0vOjU40A1jQxVISM7HxwkX8PGMwaasLhF1ADsbCc7+K9ZsZRuDIAj6yR0rV67E3//+d0ildbes+++/Hy+++CIuXryInj176j8zevRofP7551Cr1fjoo48glUoxbdq0Js//0EMPYf78+UhLS8Pq1auxbNmyNtXzo48+QkxMjP79jc+8ysjIwJQpUzBv3jx9V9CaNWvwf//3f/pjtmzZghEjRhh8rlu3bli4cCEWLlyIb775BnPmzMGaNWtw9epVuLq66o9bv349+vTpg9raWpw+fRpz586Fm5sb3n333Vb9HLfffrvB4oRDhw7FBx98AK1We8tHV8yYMUP//YABAzBw4ED07NkTiYmJGDt2rH7g+yuvvKL/91i1ahUCAgKwYcMGg9+FOTAIdSGF5dV4aMVBpBWoIRIBo4K98I/oXvop8zfz/LhgJCTn4+cTOXjsjiAMCHAxUY2JqCOIRKJ2d0+Z27lz59CjRw8UFRVh06ZNqK2txeeff67fr9VqsXLlSrz11lv6bQ4ODujVqxeAuvAUFhaGFStW4NFHH210fg8PD0ycOBGPPvooqqqqMGHCBJSVlbW6nj4+Pvoyb6RWqzF58mQMHToU//rXv/TbJ0+ejKioKP17f3//Rp8tLCzEunXr8N133yEpKQkTJkxAXFwcXFwM/z4rlUp9+X369MHFixfx6quvYvHixbC1tW10XlMICgqCp6cnUlNTMXbsWH047Nu3r/4YuVyOoKAgZGRkmKWO12PXWBdRUlmLWSsOIa1ADX9XO+x4IRqrZ0cisof7LT/b398Fk8L8IAjA3HXHUF6tueVniIg6yo4dO3Dq1ClMmzYNa9asQUBAAE6cOIGkpCT964MPPsDq1av13S83EovFePnll7Fw4cImB1YD0HdDzZo1q80P7LwZQRDw0EMPQafT4bvvvjNobXFyckKvXr30r4aHh1ZXV2PDhg2YPHky/Pz8sHLlSjz44IPIzs7GTz/9hKlTp97ykRISiQQajQY1NTWtqu/BgwcN3h84cAC9e/du0+8lKysLV65c0Qeg8PBwyOVypKSk6I+pra3FpUuXmhzDZWpdPgh99NFH6NevH/r27Ytnnnmm1TMaOoNjGVcx+dO/cDa3FJ6Ocnz/WBR6eDq06hxv3N0Pfi62uHSlAq9sOtUlf09EZHmqq6uRl5eH7OxsHDt2DG+//TbuvvtuTJw4EbNmzcKKFSswffp09O/f3+D16KOPorCwEFu3br3pue+9915IJBIsX768yf3jx49HQUGBQWuNsSxevBjbt2/Hl19+ifLycuTl5SEvL++moQwA/vGPf2Du3Lno3bs3jhw5guPHj+PZZ59tduzSlStXkJeXh6ysLGzZsgUff/wxRo8eDWdnZwBATU2NPjzW1NQgOzsbSUlJSE1NNThPRkYG4uPjkZKSgnXr1uGTTz7Bs88+e8ufs7y8HC+++CIOHDiAS5cuISEhAXfffTd69eqlH9Pl7OyMJ598EosWLcK2bduQkpKCp556CkDdv1GD0NBQbNq06ZZlGp3QheXn5wtBQUFCZWWloNFohGHDhgn79u1r8edLSkoEAEJJSUkH1rLtSiprhHe3nBOCFvwmdPvnr8LQt7cLZ3PaXtfD6Vf05/pmX7rxKkpEHaqyslI4e/asUFlZae6qtEpcXJwAQAAgSKVSwcvLS4iJiRFWrlwpaLVa4ciRIwIA4dChQ01+fsKECcI999yjP9fdd9/d6Jh33nlH8PLyEsrLywVBEAQAwqZNm5o839WrVwUAws6dO1tU/+bONWrUKP3Pdv1r1apVNz3fhQsXhNra2haVnZ6ebnBeiUQiBAQECI8//riQn59/0+MaXqNGjTKo6z/+8Q/hySefFJydnQU3Nzfh5ZdfFnQ63S3rUVFRIdx5552Cl5eXYGNjI3Tr1k14/PHHhby8PIPjampqhBdeeEFQKBSCk5OTEBMTI5w+fdrgmBt/P6NGjRKeffbZZsu/2bXfmvu3qL7wLqmgoAARERE4ceIEbG1tMXLkSKxbt85gcF1zSktL4eLigpKSEn26tgQllbX476EMfL7rIooragEAUwb54fW7+8PFzqZd5/488SLe25oMkQj48L4w3DM4wBhVJqIOVFVVhfT0dPTo0cNs40Ko84qOjsagQYM6ZNXpjnaza78192+L7hrbvXs3Jk2aBD8/P4hEIv20xOstX74c3bt3h62tLaKionDo0CH9Pi8vL8ybNw+BgYHw8/NDTExMi0OQpREEAUmZxVi4+RSGvpOAd7Yko7iiFr0UjvjPrAgsnTG43SEIAJ4cFYRZQ7tBEIB5G05iy6lcI9SeiIjIMll0EFKr1QgLC7tp/+769esRHx+PRYsW4dixYwgLC0NsbCzy8+seKXH16lX8+uuvuHTpErKzs7Fv3z7s3r37puVVV1ejtLTU4GVOWp2AI5eK8O6WZIx+PxFTlu/F9wcyUFGjRbC3I5ZMH4g/nhuJcX29jVamSCTC4kn9MD08AFqdgKfXHsN3+y8Z7fxERJ3B22+/bbA+0PWvCRMmmLt6JrNnz56b/h7a8kgSS9RpusZEIhE2bdpk8JC2qKgoDBkyBJ9++ikAQKfTQalUYu7cuZg/fz42bNiAxMREfZBasmQJBEHASy+91GQZixcvxuuvv95ou6m6xqo1WpzKKsGhS0U4nF6EI5euouy6GVx2NhKM6+uNGZFKDA3yuOXsgfbQ6gS8/OMprD+SCQB4YmQQ/jk+FJKbrEdERObDrjHjKyoquulif3Z2dk1Oee+KKisrkZ2dfdP9N1s6wFSM0TXWaReaqKmpwdGjR7FgwQL9NrFYjJiYGOzfvx9A3foK+/btQ1VVFWxsbJCYmIgnnnjipudcsGAB4uPj9e9LS0uhVCo7pP4lFbU4m1uK5LxSnMstxbncMqSoylCjMXzisrOtFKNDFRjbxxtjQxVwkJvmn0wiFuHdaQMQ4GaHD/48j692p+FUVgk+vn9Qo0d1EBF1Ne7u7vrnhlkzOzs7s4edjtZpg1BhYSG0Wi28vQ27hby9vZGcXPd8rdtvvx133XUXBg8eDLFYjLFjx2Ly5Mk3PadcLr/pQ+mMKae4EsPe3dHkPg8HGYZ0d0dkj7pXH19ns7XCiEQizB3bG908HTB/40nsT7uCuz7+C0umD8ToUMWtT0BEJtVJGviJjMYY13ybg9DPP//c6s+MGzdOv3CUqbz11lsGK49aAl8XWzjbSuFka4M+vs7o6+uEPr7O6OPrjG4e9h3a5dUWk8P80NfXGU+vOYYUVRlmrz6Me8MDsHBiX6MM0Cai9mlY9K6mpsbkf2OJzKmiogIAYGPT9ntRm4PQ9WN1WkIkEuHChQsICgpqa5EGPD09IZFIoFKpDLarVCr4+PgYpYyOIhKJcPDlGNjJjLuSaUfqpXDET3OG4/0/UrBibzo2HM3CjuR8vDQ+BPeGK2/6LDMi6nhSqRT29vYoKCiAjY0NxGKLngdD1G6CIKCiogL5+flwdXVt18rg7eoay8vLg0LRsi4SJyen9hTViEwmQ3h4OBISEvShTKfTISEhAXPmzDFqWR2hM4WgBrY2Eiyc2Bex/X0wf+NJXCxQ458bT+H7AxlYPLkfwru5mbuKRFZJJBLB19cX6enpuHz5srmrQ2Qyrq6u7W78aFUQkkgk+ue6xMXFtaoJ9qGHHmr1zKvy8nKDZcDT09ORlJQEd3d3BAYGIj4+HnFxcYiIiEBkZCSWLl0KtVqN2bNnt6ocap0h3d2x9bmR+GbfJXy8/QJOZZdg2uf7cM9gfzwfE4xAD3tzV5HI6shkMvTu3bvVz5gi6qxsbGyM8oy4Vk2fF4vF0Ol0tz7QSBITEzF69OhG2+Pi4rB69WoAwKeffoolS5YgLy8PgwYNwrJlywye6tselrqytCUpKKvGkj+S8b8jWQDqZptNvy0Ac8b0gtKdgYiIiEyvNffvVgWh61uEAECj0UAq7bQTz26JQajlTmYV44Nt57HrfAEAQCoW4b4hSjw1qicDERERmZTJHrERGBiI9957D8XFxe05DXUBAwNc8c0jkdj41FDc0csTGp2AtQczEP1+IuauO45TWSXmriIREVEj7QpCzz//PD777DMolUo8++yzSE9PN1a9qJMK7+aO7x+Lwv/+ry4QaXUCfjmRg0mf/oX7vzqAHckq6HRc64SIiCxDu7rGAECr1eJ///sfli5dimPHjuHuu+/Giy++aLRxOubErrH2O5NTgq/3pOOXEznQ1AcgpbsdHojshvsiAuDh2PELWBIRkXUx2RihG+3ZswcffvghfvnlF0RFRWHevHmYMmWKxS0Q2FIMQsaTU1yJVXvTsf5wJkqr6p6fJpOIcdcAHzx0ezeEd3PrtNcJERFZFrMFoQZpaWlYunQpvvnmGygUCly4cKGlRVgUBiHjq6zR4pcTOfj+4GWcvG7cUKiPEx6ICsTkMD+42svMWEMiIursTBaEFi1ahJKSkiZfxcXFKC4uhk6na1F4skQMQh3rZFYxvj9wGT+fyEFVbd2yDDKJGOP6emN6eABG9PaEVMIVcomIqHVMFoTEYjFsbW3x8MMPo0+fPnBxcYGzszOcnZ3137u4uLR49WlLwyBkGiUVtdh4LAv/O5KJ5Lwy/XYvJznuGeyP6eEBCPY27srkRETUdZksCO3cuRMffPABtm/fjhkzZmDevHno379/22tuYRiETO9MTgl+OJqFn5JyUKS+tkJuPz9nTA7zw8QwP/i78qGSRER0cyYfI5SSkoIPP/wQ33//PUaMGIEXX3wRY8eObX3NLQyDkPnUaHRITMnHD/UPd9VcN+U+opsbJg/yw10DfOHJWWdERHQDsw2WLigowPLly/HFF1/A19cX8+bNw4wZM4zyLBBzYBCyDEXqGvx+Khe/nMjBoUtFaLhiJWIRhvX0wOQwP8T294GzrY15K0pERBbB7LPGiouL8cknn2Dp0qVwdHTstE9DZhCyPLkllfjtZF0oOnHdrDOZRIyRwV6Y0N8HMX284WLPUEREZK1MFoSmTZvW5Iyx2tpaNJxWJBJx1hh1iEuFavxyIgc/n8jBhfxy/XapWIRhvTwxob8P7uzrzUUbiYisjMmC0MyZM+Hq6nrLl4eHR9t/GjNiEOocBEFAiqoMW07lYevpPKSors08E4uAId3dMaG/D8b394WPi60Za0pERKZg9q6xroJBqHNKKyjHltN1oehUtuHDXgcHutaFon6+CPSwN1MNiYioIzEIGQmDUOeXdbUCW+tD0dGMq7j+ag/1cUJMH2/E9PXGQH8XiMV8xAcRUVfQYUFILBZDp9O1u4KdBYNQ15JfWoU/zuRhy+k8HEwvgva6KfkKJznG9vHGuL4KDOvpCVubzjnTkYiIOjAIrVmzBiUlJQgJCUF0dHSnnRbfUgxCXVdxRQ0SUwrw5zkVdqUUoLxao99nZyPByGBPxPTxxphQBQdbExF1Mh0WhBqkpqZi586dqK2tRWBgIMaNGwe5vOvdLBiErEO1RouDaUX486wK28+pkFtSpd8nFgHh3dwQ08cb4/p6I8jL0Yw1JSKilujwIHS97OxsbNu2DZWVlVAoFBg/fjwcHbvGzYJByPoIgoAzOaX6UHQmp9Rgf5CXA8bVjysarHTlQ2GJiCyQSYPQ9QoLC7F161Z9BcaPH99pp84DDEIEZBdXIuGcCn+eVeFA2hXUaq/95+Jqb4PoYC+M6eONUb29uIgjEZGFMFsQAoATJ05g48aN2LZtG2677TZ88sknnXYsEYMQXa+0qha7zxfgz7MqJKYUoKSyVr9PIhYhvJsbxoYqMLaPAj29HCEScRYaEZE5mDwIXbhwAV9//TUOHjyIgQMHYtq0aRgxYgTE4s7dbcAgRDej0epwLKMYCckq7EzOx3lVucF+pbsdxobWDbaOCnKHXNo5/2eAiKgzMnkQSklJwebNm/HPf/6zvaeyKAxC1FKZRRXYkZyPhOR8HLh4BTXaa8tM2MskuKOXJ8b2UWB0iAIKZ65uTUTUkczSNVZbWwsbm641RoJBiNpCXa3B3tRC7EjOx47kfOSXVRvsHxjggtEhdV1o/f24kCMRkbGZJAidPHkS/fv3b3H315kzZxASEgKpVNqW4tosPT0djzzyCFQqFSQSCQ4cOAAHB4cWfZZBiNpLpxNwNrcUCefysSNZhRNZho/88HKSY0yIAqNDFRjR2xMOctP+90FE1BWZJAhJJBLk5eXBy8urRcc7OzsjKSkJQUFBbSmuzUaNGoU333wTI0aMQFFREZydnVscxhiEyNjyy6qQmFKAHefysedCAdQ11x5ZI5OIERXkjrGhCowJ9eaz0IiI2qg19+82/++nIAh49dVXYW/fsj/WNTU1bS2qzc6cOQMbGxuMGDECAODu7m7yOhBdT+Fki/silLgvQolqjRaH0ovqxhady0dGUQX2XCjEnguFWPzLWfRSOGJsaF1rUXg3N9hwzSIiIqNrc4tQdHR0q6cHr127Fr6+vi0+fvfu3ViyZAmOHj2K3NxcbNq0CVOmTDE4Zvny5ViyZAny8vIQFhaGTz75BJGRkQCAzZs3Y/Xq1dBqtcjOzsb06dPx8ssvt7h8tgiRqQiCgIsFauxMzkdCsgqHL101eBaak60UI3t7YVSIF6KDvTjgmoioGSZpEUpMTGzrR1tMrVYjLCwMjzzyCKZOndpo//r16xEfH48vvvgCUVFRWLp0KWJjY5GSkgKFQgGNRoM9e/YgKSlJv+r1kCFDMG7cuCbLq66uRnX1tYGtpaWlTR5HZGwikQi9FI7opXDE4yODUFJZt2bRjuR87DpfgCJ1DX47lYvfTuUCAPr7O2N0iALRIQoMUrpCwgHXRERtYvQFFTuKSCRq1CIUFRWFIUOG4NNPPwUA6HQ6KJVKzJ07F/Pnz8f+/fuxePFi/PHHHwCAJUuWAABefPHFJstYvHgxXn/99Ubb2SJE5qTVCTiZVYydKQVITMnHyRsGXLva22BUsBdGhygwMtgL7g4yM9WUiMgymHVl6Y5yYxCqqamBvb09fvjhB4NwFBcXh+LiYvz000/QaDQYMmQIduzYARcXF9x99934v//7P0ycOLHJMppqEVIqlQxCZFEKyqqx+3wBdqbkY/f5ApRWafT7RCIgLMAVo0MUGB3qxen5RGSVTNI1Zm6FhYXQarXw9vY22O7t7Y3k5GQAgFQqxdtvv42RI0dCEATceeedNw1BACCXyyGXyzu03kTt5eUkx7TwAEwLD4BGq8PxzGLsTM7HzpQCnMstRVJmMZIyi/HR9vPwdJRhVLAC0SFeGMnnoRERNdJpg1BLTZgwARMmTDB3NYg6hFQixpDu7hjS3R0vjQ9FXkkVdp3Px87kAvyVWojC8hpsPJaFjceyIBGLcFugK6JD6la47uPrxOehEZHV67RByNPTExKJBCqVymC7SqWCj4+PmWpFZF4+Lrb4+5BA/H1IIGo0Ohy5XITElALsTM7HhfxyHL50FYcvXcWSP1Lg7SyvH3DtheG9POFky9YiIrI+RgtCe/bswZdffomLFy/ihx9+gL+/P7777jv06NEDd9xxh7GK0ZPJZAgPD0dCQoJ+jJBOp0NCQgLmzJlj9PKIOhuZVIxhPT0xrKcnXr6rDzKLKpB4vgC7UvKxN/UKVKXV+O/hTPz3cCakYhGGdHfH6NC6Qde9FI5sLSIiq2CUILRx40bMnDkTDz74II4fP64fcFxSUoK3334bv//+e5vOW15ejtTUVP379PR0JCUlwd3dHYGBgYiPj0dcXBwiIiIQGRmJpUuXQq1WY/bs2cb4sYi6FKW7PWbe3g0zb++Gqtq6xRx3puQjMaUA6YVq7E+7gv1pV/D278nwd7XTh6KhPT1gL+u0jcdERM0yyqyxwYMH4/nnn8esWbPg5OSEEydOICgoCMePH8eECROQl5fXpvMmJiZi9OjRjbbHxcVh9erVAIBPP/1Uv6DioEGDsGzZMkRFRbXnx9HjgopkLdIL1UisD0X7066gRqPT75NJxYjq4Y7o+m60IE8HthYRkUUz+fR5e3t7nD17Ft27dzcIQmlpaejbty+qqqraW4RZMAiRNaqs0WJ/WiF2Jtct6JhdXGmwX+luh+hgBUYFe2FYL7YWEZHlMfn0eR8fH6SmpqJ79+4G2//66y+TP2SViNrHTibBmFBvjAn1xr8EARcLyrEzuQCJ5/NxKL0ImUWV+O7AZXx34DJkEjEie7hjVLAXokO8OLaIiDodowShxx9/HM8++yxWrlwJkUiEnJwc7N+/H/PmzcOrr75qjCKIyAzqHv3hhF4KJzw+Mgjqag32X7yCxPN13WhZVyvxV2oh/kotxFu/n4O/qx1G1oeiYT09OBONiCyeUbrGBEHA22+/jXfeeQcVFRUA6hYnnDdvHt544412V9Jc2DVGdHOCICCtUI3E+kd/HEwvMhhbJBWLENHdDdEhdd1ooT5ct4iITMNsj9ioqalBamoqysvL0bdvXzg6Ohrr1GbBIETUcpU1WhxIu4LElLoHxV66UmGw38fZFqOCvTCqft0iFzu2FhFRx+iSzxozBwYhora7VKjGrvN1rUX7066gqvZaa5FELEJ4oBtGhXhhVLAX+vk5s7WIiIzGJEEoPj6+xcd++OGHbSnC7BiEiIyjYd2ixJS6QddpBWqD/V5OcozsXTe2aERvT7jay8xUUyLqCkwShG5c3+fYsWPQaDQICQkBAJw/fx4SiQTh4eHYsWNHW4owOwYhoo5xbZXrAuy7WIiKGq1+n1gEDFK66tct6u/nArGYrUVE1HIm7xr78MMPkZiYiG+++QZubm4AgKtXr2L27NkYMWIEXnjhhfYWYRYMQkQdr1qjxZFLV/XdaOdV5Qb7PRxkGBlc14U2MtgL7g5sLSKi5pk8CPn7+2Pbtm3o16+fwfbTp0/jzjvvRE5OTnuLMAsGISLTyy6uxO76ULQ39QrKqzX6fSIRMDDAVb9uUViAKyRsLSKiG5h8QcXS0lIUFBQ02l5QUICysjJjFEFEVsLf1Q73Rwbi/shA1Gh0OJZxFYkpBdh1vgDncktxIrMYJzKLsSzhAlztbTCitxei61uLvJzk5q4+EXUyRmkRmjVrFvbs2YMPPvgAkZGRAICDBw/ixRdfxIgRI/DNN9+0u6LmwBYhIsuiKq3CrvpQtPtCAcqqNAb7+/s71z3+I8QLg5WukErEZqopEZmTybvGKioqMG/ePKxcuRK1tbUAAKlUikcffRRLliyBg4NDe4swCwYhIsul0epwPLMYu+pnop3OLjXY72wrxYjeXvq1i7ydbc1UUyIyNbOtI6RWq3Hx4kUAQM+ePTttAGrAIETUeeSXVWHP+UIkni/AngsFKK6oNdjfx9dZP7YovJsbbNhaRNRlcUFFI2EQIuqctDoBJ7KK9WOLTmYV4/q/dI5yKYb38tA//sPP1c58lSUiozN5EPrXv/7V7P7XXnutvUWYBYMQUddwpbwaey4U1o0tOl+AK+oag/3B3o76UBTR3Q1yqcRMNSUiYzB5EBo8eLDB+9raWqSnp0MqlaJnz544duxYe4swCwYhoq5HpxNwKrtEv25RUmYxdNf9FbSXSTCspydGhdTNRlO625uvskTUJhbRNVZaWoqHH34Y99xzD2bOnNkRRXQ4BiGirq+4ogZ7LhTqu9EKy6sN9vf0csCo4LpVriN7uMPWhq1FRJbOIoIQAJw6dQqTJk3CpUuXOqqIDsUgRGRddDoBZ3NLsav+8R9HM65Ce11zkZ2NBEN7eiA6xAvRwQoEerC1iMgSmXxBxZspKSlBSUlJRxZBRGQ0YrEI/f1d0N/fBU+P7oWSylrsTS1EYko+dp0vgKq0GjuS87EjOR/AGfTwdNDPRLs9yIOtRUSdkFFahJYtW2bwXhAE5Obm4rvvvsOoUaOwdu3a9hZhFmwRIqIGgiAgOa8MiSl1Y4uOXr4KzXWtRXKpGEN7etQHIwV6eHbu5UOIOjOTd4316NHD4L1YLIaXlxfGjBmDBQsWwMnJqb1FmAWDEBHdTFlVXWtR3aDrAuSWVBns7+Zhj+j6UHR7kAfsZGwtIjIVkwehy5cvQ6lUQiw2XKBMEARkZmYiMDCwvUWYBYMQEbWEIAg4rypHYko+ElMKcORyEWq11/60yqRiRPVwR3RI3aDrIE8HiER8WCxRRzF5EJJIJMjNzYVCoTDYfuXKFSgUCmi12vYWYRYMQkTUFuXVGuxLrVvleldKAbKLKw32K93t6rrQghUY1ssD9rIOHa5JZHVMHoTEYjHy8vIaBaHLly+jb9++UKvV7S3CLBiEiKi9BEHAxYLy+rFFBTiUXoQarU6/XyYRI7KHO6JD6p6L1kvhyNYionYyWRCKj48HAHz88cd4/PHHYW9/bSqpVqvFwYMHIZFIsHfv3rYWYVYMQkRkbOpqDfZfvILE83XdaFlXDVuL/F3t9Is5DuvlCUc5W4uIWstk0+ePHz8OoO7/eE6dOgWZTKbfJ5PJEBYWhnnz5rWnCKOoqKhAnz59cO+99+L99983d3WIyIo5yKWI6euNmL7eEAQBaYVq/Uy0g+lFyC6uxNqDGVh7MAM2EhEiutW1FkWHKBDszdYiImMzStfY7Nmz8fHHH1tsq8krr7yC1NRUKJXKVgUhtggRkSlV1mhxIO1K3aDr8wW4fKXCYL+vi61+3aLhvTzhZGtjppoSWTaTL6i4atUqY5ymQ1y4cAHJycmYNGkSTp8+be7qEBHdlJ1MgtGhCowOrRtvmV6oxq76ULT/4hXkllThv4cz8d/DmZCKRQjv5qafiRbq48TWIqI2aHOLUHx8PN544w04ODjoxwrdzIcfftimyu3evRtLlizB0aNHkZubi02bNmHKlCkGxyxfvhxLlixBXl4ewsLC8MknnyAyMlK//+6778aSJUuwb98+nD59mi1CRNQpVdU2tBbVPRMtvdBwEoq3s1y/mOPwXp5wsWNrEVkvk7QIHT9+HLW1tfrvb6Y9/4eiVqsRFhaGRx55BFOnTm20f/369YiPj8cXX3yBqKgoLF26FLGxsUhJSYFCocBPP/2E4OBgBAcHY9++fbcsr7q6GtXV1x64WFpa2ua6ExEZk62NpL71p6616PIVtX4xx30XC6Eqrcb/jmThf0eyIBGLEB7ohlH1M9H6+TmztYjoJjr0oavGJBKJGrUIRUVFYciQIfj0008BADqdDkqlEnPnzsX8+fOxYMECfP/995BIJCgvL0dtbS1eeOEFvPbaa02WsXjxYrz++uuNtrNFiIgsWVWtFocvFekHXV8sMGwt8nKS68cWjejlBRd7thZR12YxT583phuDUE1NDezt7fHDDz8YhKO4uDgUFxfjp59+Mvj86tWrb9k11lSLkFKpZBAiok4ls6hCv5jjvouFqKi5tqitWAQMDnTTP/6jn58zxGK2FlHXYpKusVuNC7peW8cINaewsBBarRbe3t4G2729vZGcnNymc8rlcsjlcmNUj4jIbJTu9ph5ezfMvL0bqjVaHLl0FYkp+dh1vgDnVeU4evkqjl6+ig/+PA9PRxlG9vbCqBAvjOztBTcH2a0LIOpC2jVGqCUspV/64YcfNncViIhMTi6VYHgvTwzv5YlX/gZkF1diV30X2t7UQhSW1+DH49n48Xg2xCIgTOmK6GAFRoV4YaC/C1uLqMtrcxDauXNnk9sbeto6OgB5enpCIpFApVIZbFepVPDx8enQsomIOit/Vzs8EBWIB6ICUaPR4ejlq0g8n49dKQVIzivD8YxiHM8oxkfbz8PdQYaRvT0RHaLAyGAvuLO1iLog8a0PaZkVK1agf//+sLW1ha2tLfr374+vv/7aWKdvRCaTITw8HAkJCfptOp0OCQkJGDp0aIeVS0TUVcikYgzt6YEFE/pg63MjsX/BGLw7dQAm9PeBk1yKInUNNifl4Ln1SQh/80/cvXwvPvzzPI5lXIVW1ymGlxLdklEWVHzttdfw4YcfYu7cufoQsn//fjz//PPIyMjAv/71rzadt7y8HKmpqfr36enpSEpKgru7OwIDAxEfH4+4uDhEREQgMjISS5cuhVqtxuzZs43xYxERWRVfFzvMiAzEjMhA1Gp1OHb5KhLrp+ifyy3FicxinMgsxrKEC3Czt8GI3nUz0UYGe8HTkeMrqXMyyqwxLy8vLFu2DPfff7/B9nXr1mHu3LkoLCxs03kTExMxevToRtvj4uKwevVqAMCnn36qX1Bx0KBBWLZsGaKiotpU3o24oCIRUR1VaRV21c9E232hAGVVGoP9A/xd6p+J5oVBSjdIOLaIzMjk0+ddXV1x+PBh9O7d22D7+fPnERkZieLi4vYWYRYMQkREjWm0OhzPLK4bdH0+H6ezDRefdbW3wahgL4wJVXAmGpmFyYPQ3LlzYWNj02ia/Lx581BZWYnly5e3twizYBAiIrq1/LIq7D5fiF3nC7D7fAFKKmv1+xrWLRoTqsDoEAX6+PKZaNTxzBKEvv32WyiVStx+++0AgIMHDyIjIwOzZs2Cjc21VUw7Yk2hjsIgRETUOg2tRTuS87EzOR/JeWUG+32cbTE61Auj65+J5iA3ylBVIgMmD0JNjeNpsjCRCDt27GhvcSbDIERE1D45xZXYmVIXivamXkFl7bVVrmUSMaKC3DE6RIExoQp093QwY02pK+mSj9gwBwYhIiLjqarV4mB6EXYm52NHcj4yiioM9gd5OmB0fRdaZA93yKRGW+GFrAyDkJEwCBERdQxBEHCxQI3ElLpQdCi9CJrr1iZykElwR29PjAlVIDpEAW9nWzPWljobswShqqoqnDx5Evn5+dDpdAb7Jk+ebIwiTI5BiIjINMqqavHXhcK6brSUAhSUVRvs7+fnXDfgOlSBsABXTs+nZpk8CG3duhWzZs1qcr0gkUgErVbbxKcsH4MQEZHp6XQCzuSUYkdyPnak5ONkVjGuv1O5O8gwKtgLo0MVGNXbCy72Njc/GVklkweh3r17484778Rrr73W6GnwnRmDEBGR+RWWV2NXSgF2pORj93nDxRwlYhHCA90wOlSBcX0V6OnlyOn5ZPog5OzsjOPHj6Nnz57tPZVFYRAiIrIstdq6B8U2zEQ7ryo32N/dwx5j+3gjpo83Irq7wUbCAdfWyORB6JFHHsHw4cPx6KOPtvdUFoVBiIjIsmUWVSAxJR/bz+Vj/8UrqNFeG6PqbCvF6FAFxvbxxqhgL7jYsQvNWpg8CFVUVODee++Fl5cXBgwYYLCAIgA888wz7S3CLBiEiIg6j/JqDf66UIA/z+ZjZ0o+itQ1+n1SsQiRPdwRU99aFOhhb8aaUkczeRBasWIFnnzySdja2sLDw8Ogf1YkEiEtLa29RZgFgxARUeek1Qk4nnEV28/lY/s5FVLzDbvQgr0d60JRX28MCnCFmLPQuhSTByEfHx8888wzmD9/PsTirtMfyyBERNQ1XCpUY/s5FbafU+HwpavQXrdmkaejDGPqu9BG9PaEvYyP/ejsTB6E3N3dcfjwYQ6WJiIii1dcUYNd5wvw51kVdqUUoKz62iw0mVSM4T09ENPXG+P6eEPBhRw7JZMHoeeffx5eXl54+eWX23sqi8IgRETUtdVodDh8qQh/nq1rLcq6WqnfJxIBg5WuiO3ng9h+PnwWWidi8iD0zDPP4Ntvv0VYWBgGDhzYaLB0Z3ri/PUYhIiIrIcgCDivKsf2cypsO6vCicxig/0h3k6I7eeNO/v5oJ+fM9crsmAW9fT5zvbE+esxCBERWa/ckkr8eVaFP87k4UBakcG4ogA3O9zZ1wex/bwR0d2dj/ywMBb10NXTp0+jf//+HVlEh2EQIiIioG5cUcK5fPxxJg+7LxSgqvbaekUeDjLE9PHGnf28MbyXJ2xtJGasKQEWEITKysqwbt06rFixAkeOHOGzxoiIqMuorNFi1/kCbDuTh+3nVCi97pEfDjIJokMUiO3vgzGhCjjKOQPNHMwWhHbv3o0VK1Zg48aN8PPzw9SpUzFt2jQMGTLEWEWYFIMQERE1p1arw8G0IvxxJg/bzuZBVVqt3yeXijEq2At/G+iLsX28GYpMyKRBKC8vD6tXr8aKFStQWlqK++67D1988QVOnDiBvn37tufUZscgRERELaXTCTiRVYw/zqiw9XQuLl2p0O+TNYSiAb4Y20cBJ1s+7qMjmSwITZo0Cbt378bf/vY3PPjggxg/fjwkEglsbGwYhIiIyGoJgoCzuaX4/VQufj+Vh/RCtX6fTCrGyN5e+NtAH8T08WYo6gAmC0JSqRTPPPMMnnrqKfTu3Vu/nUGIiIiojiAIOJdbVh+KcpF2fSiSiDEy2BN3DfBFTF9vODMUGUVr7t/t6rD866+/sGLFCoSHh6NPnz6YOXMmZsyY0Z5TEhERdSkikQh9/ZzR188ZL9wZjOS8ulD026lcpBWo65+Hlg+ZRIwRvetC0bh+DEWmYpTB0mq1GuvXr8fKlStx6NAhaLVafPjhh3jkkUfg5ORkjHq2WWZmJmbOnIn8/HxIpVK8+uqruPfee1v0WbYIERFRR2lYwPG3+pai6x8MK5OIER3ihcmD/DA21Bt2Mk7Jbw2zTp9PSUnBihUr8N1336G4uBjjxo3Dzz//bMwiWiU3NxcqlQqDBg1CXl4ewsPDcf78eTg43HqpdAYhIiIylfOqMvx2sq6l6PpQ5CCT4M5+Ppgc5oc7envCRtJ1Hm7eUcy+jhAAaLVa/PLLL1i5cqVZg9CNwsLC8Ouvv0KpVN7yWAYhIiIyNUEQkJxXhp9P5ODnpBxkF197/pmbvQ0mDPDF5DA/RHZ3h5grWjepNffvDouVEokEU6ZMaXcI2r17NyZNmgQ/Pz+IRCJs3ry50THLly9H9+7dYWtri6ioKBw6dKjJcx09ehRarbZFIYiIiMgcRCIR+vg645/jQ/HXP0dj41PDEDe0GzwdZbhaUYu1BzMw46sDGP7eDrz121mcyipBBz8kokuz+PY1tVqNsLAwLF++vMn969evR3x8PBYtWoRjx44hLCwMsbGxyM/PNziuqKgIs2bNwldffXXTsqqrq1FaWmrwIiIiMheRSITwbm54/e7+OLBgLL57NBL3hgfAyVaK3JIq/GdPOiZ9+hfGfrALH/553mCaPrVMhz9rzJhEIhE2bdqEKVOm6LdFRUVhyJAh+PTTTwEAOp0OSqUSc+fOxfz58wHUBZxx48bh8ccfx8yZM296/sWLF+P1119vtJ1dY0REZEmqarVITCnALydysP2cCtWaa88+G6R0xbTb/DFxoB/cHGRmrKX5WMQYoY5wYxCqqamBvb09fvjhB4NwFBcXh+LiYvz0008QBAEPPPAAQkJCsHjx4mbPX11djerqa8ujl5aWQqlUMggREZHFKq/WYNuZPPyUlIM9Fwqgq7+r20hEiA5RYNpt/hgdqoBcaj0zz0y2jpC5FRYWQqvVwtvb22C7t7c3kpOTAQB79+7F+vXrMXDgQP34ou+++w4DBgxodD65XA65XN7h9SYiIjIWR7kUU28LwNTbApBfVoWfk3Lw47FsnM0txZ9nVfjzrAoudjb420BfTLvNH7cFukEk4iDrBp06CLXEHXfcAZ1Od+sDiYiIOjmFky0eGxGEx0YEISWvDD8ez8JPx3OQV1qFtQczsPZgBrp52GPKIH9Mvc0f3TxuvZRMV9epg5CnpyckEglUKpXBdpVKBR8fHzPVioiIyPxCfJywYEIfvBQbiv0Xr+DH41nYejoPl69U4OOEC/g44QLCu7lhav14Ihc761zJ2uJnjTVHJpMhPDwcCQkJ+m06nQ4JCQkYOnSoGWtGRERkGSRiEe7o7YkP7xuEIwtj8NHfwzCityfEIuDo5at4ZdNpRL61Hc/+9zj+ulAIna7TDB02CotvESovL0dqaqr+fXp6OpKSkuDu7o7AwEDEx8cjLi4OERERiIyMxNKlS6FWqzF79mwz1pqIiMjy2MukuGdwAO4ZHABVaRV+SsrGD0ezcF5Vjp+ScvBTUg78Xe0wLTwA94YHQOlub+4qdziLnzWWmJiI0aNHN9oeFxeH1atXAwA+/fRTLFmyBHl5eRg0aBCWLVuGqKiodpfNlaWJiKirEwQBJ7NKsOFoJn5KykFZlUa/b2iQB+6NCMCE/r6d6nlnXXb6vKkxCBERkTWpqtXijzN5+OFoFv5KLURDQnCUSzEpzBfTw5W4LdDV4medMQgZCYMQERFZq+ziSmw8moUfjmYho6hCv72nlwPujVBi6mB/KJxtzVjDm2MQMhIGISIisnY6nYCD6UXYcDQTW07lobJWC6BuEHZ0sBfujVBibB8FbCSWM/+KQchIGISIiIiuKauqxW8nc7HhaBaOXr6q3+7pKMe9EQH4e4QS3T3NvzYRg5CRMAgRERE1LTW/HD/Ud50Vll97PNXQIA/MiFRifH8fsz3Wg0HISBiEiIiImler1SHhXD7+ezgDu84X6AdYu9rbYOrgANwfqURvbyeT1olByEgYhIiIiFouu7gS/zucif8dyURuSZV+e3g3N8wYosTEgX4mmYbPIGQkDEJEREStp9UJ2H2+AOsOZSAhOR/a+tWqneRS3D3YDzOGBKK/v0uHlc8gZCQMQkRERO2TX1qFDUezsP5wpsE0/AH+LpgRqcTkMD842Rr3OWcMQkbCIERERGQcOp2A/WlXsO5QBradUaFGqwMA2MskSHhhFHxd7IxWVmvu3xb/rDEiIiLq/MRiEYb38sTwXp4oUtfgx2NZWHcoA3YyiVFDUGuxRagZbBEiIiLqOIIgoEhdAw9HuVHP25r7t+UsA0lERERWRSQSGT0EtRaDEBEREVktBiEiIiKyWgxCREREZLUYhIiIiMhqcfp8Mxom1JWWlpq5JkRERNRSDfftlkyMZxBqRllZGQBAqVSauSZERETUWmVlZXBxaf5RHlxHqBk6nQ45OTlwcnKCSCRq9tghQ4bg8OHDrTp/az7TkmNvdUxz+5vaV1paCqVSiczMTIteR6ktv3tTn7+t52jp54xxXFv28Rox3vn5N8R8Ovr6MFYZvEZaThAElJWVwc/PD2Jx86OA2CLUDLFYjICAgBYdK5FIWv2P2JrPtOTYWx3T3P7m9jk7O1v0H7G2/O5Nff62nqOlnzPGcW3dB/AaMdc1wr8hxtHR14exyuA10jq3aglqwMHSRvL000936Gdacuytjmluf1vqbyk6uu7GOH9bz9HSzxnjuK56fQBd9xrh3xDjMEXdeY1Y7jXCrjG6KT5ihG6F1wg1h9cH3YolXCNsEaKbksvlWLRoEeRy8y5/TpaL1wg1h9cH3YolXCNsESIiIiKrxRYhIiIisloMQkRERGS1GISIiIjIajEIERERkdViECIiIiKrxSBEbfLrr78iJCQEvXv3xtdff23u6pAFuueee+Dm5obp06ebuypkgTIzMxEdHY2+ffti4MCB2LBhg7mrRBakuLgYERERGDRoEPr374///Oc/HVYWp89Tq2k0GvTt2xc7d+6Ei4sLwsPDsW/fPnh4eJi7amRBEhMTUVZWhm+++QY//PCDuatDFiY3NxcqlQqDBg1CXl4ewsPDcf78eTg4OJi7amQBtFotqqurYW9vD7Vajf79++PIkSMdcp9hixC12qFDh9CvXz/4+/vD0dEREyZMwLZt28xdLbIw0dHRcHJyMnc1yEL5+vpi0KBBAAAfHx94enqiqKjIvJUiiyGRSGBvbw8AqK6uhiAI6Kh2GwYhK7R7925MmjQJfn5+EIlE2Lx5c6Njli9fju7du8PW1hZRUVE4dOiQfl9OTg78/f317/39/ZGdnW2KqpOJtPcaoa7PmNfI0aNHodVqoVQqO7jWZCrGuD6Ki4sRFhaGgIAAvPjii/D09OyQujIIWSG1Wo2wsDAsX768yf3r169HfHw8Fi1ahGPHjiEsLAyxsbHIz883cU3JXHiN0K0Y6xopKirCrFmz8NVXX5mi2mQixrg+XF1dceLECaSnp2Pt2rVQqVQdU1mBrBoAYdOmTQbbIiMjhaefflr/XqvVCn5+fsI777wjCIIg7N27V5gyZYp+/7PPPiusWbPGJPUl02vLNdJg586dwrRp00xRTTKjtl4jVVVVwogRI4Rvv/3WVFUlM2jP35AGTz31lLBhw4YOqR9bhMhATU0Njh49ipiYGP02sViMmJgY7N+/HwAQGRmJ06dPIzs7G+Xl5diyZQtiY2PNVWUysZZcI2TdWnKNCIKAhx9+GGPGjMHMmTPNVVUyg5ZcHyqVCmVlZQCAkpIS7N69GyEhIR1SH2mHnJU6rcLCQmi1Wnh7exts9/b2RnJyMgBAKpXigw8+wOjRo6HT6fDSSy9xxpgVack1AgAxMTE4ceIE1Go1AgICsGHDBgwdOtTU1SUzaMk1snfvXqxfvx4DBw7Ujx/57rvvMGDAAFNXl0ysJdfH5cuX8cQTT+gHSc+dO7fDrg0GIWqTyZMnY/LkyeauBlmw7du3m7sKZMHuuOMO6HQ6c1eDLFRkZCSSkpJMUha7xsiAp6cnJBJJo0FpKpUKPj4+ZqoVWRJeI3QrvEaoOZZ2fTAIkQGZTIbw8HAkJCTot+l0OiQkJLBbgwDwGqFb4zVCzbG064NdY1aovLwcqamp+vfp6elISkqCu7s7AgMDER8fj7i4OERERCAyMhJLly6FWq3G7NmzzVhrMiVeI3QrvEaoOZ3q+uiQuWhk0Xbu3CkAaPSKi4vTH/PJJ58IgYGBgkwmEyIjI4UDBw6Yr8JkcrxG6FZ4jVBzOtP1wWeNERERkdXiGCEiIiKyWgxCREREZLUYhIiIiMhqMQgRERGR1WIQIiIiIqvFIERERERWiwsqNkOn0yEnJwdOTk4QiUTmrg4RERG1gCAIKCsrg5+fH8Ti5tt8GISakZOTA6VSae5qEBERURtkZmYiICCg2WMYhJrh5OQEoO4X6ezsbObaEBERUUuUlpZCqVTq7+PNYRBqRkN3mLOzM4MQERFRJ9OSYS0cLE1ERERWi0GIiIiIrBaDEBEREVktjhEi6gA6nQCNToBOEKBt+P6GbfrXje/rt+mu+1zDMUBdn7cIgEgEiCCq/wrghvci0fXf1x3Q8F4iFl17iUSG72+yTSwSQVr/PZeTIKKugkGI2kUQBFRrdCir0qCsqhZlVRqoqzWo1upQo6l71Wqvfa3W6FCj1aFWUxcIBKH+pi+g/j2g1Qn67xtCg06oK+v673U3fE533TahRcdc22dwTl3Tx2t0Ouh0uBZsmgk0XZ1IBEivC0fi+oAkFYthIxHBRnL912vfy6RiSMX126ViyOr3SSXXvreRiOvfX/d5qRg24mvnkEvFsLWRXPtqI4at1PCrXCqBRMzARkTNYxCiZpVXa3A2pxTpheXILq5CTnElcoorkVtShZLKWpRV1aJW2/Vv/MbU0OIiFgNSsRhiUUMLjRgSMepaYyQNx1xrnRHXt8IIqAttQH3gQ11ou36fUP/m+vfXH6u7Ibjd2PLUEDhvRhBQ/+8uoLrDflPtZyMRGQSjpgKTrVQCe5kE9nIJHGRS2MnqvtrL67fLpAbvHWRS/XZbGzFbx4g6OQYhMiAIAo5lFGPb2TzsOJeP1IJy/U22OSIR4CiTwslWCge5FHKbuv/Dl0nr/o9eXv9VVt8KIJWI61sU6rpwxCKRPhCI6r8X1wcB/ff1x0rEIn33jsGxN5xLXN81JK4PHXXvW3b8jV9v1mV0fXeRWCzSt5LcrJupoczOoKEFrqH1Sx+WrmsB02gN99VqBWh0da1/tVqh/qsONZrrtmuEulZBrQ4a7bXvG73X1H9eJ6C2oWWxvnWxuuFVq0W1Roeq675qrktwdXXQoKyD0ppIhOvCkwSOtlI429rASf/VBs520rqvtlL9e2dbG/1xTrZSSCUcrklkLgxCpHcs4yre+f0cDl+6arDdz8UWwT5O8HO1g7+rHfxcbeHjbAc3h7o/9E62UjjKpBCzG6JLEYlEkEpEkErMXZPW0WivBaWqJoJSo6+1WlTUNLw0UNdoUVmjhbpag4oaLdQ1mrr3NRpUVNcdV1mrBVDXMlZerUF5tQYF7aizvUwCZ1sbuNrbwN1BBjd7GdwcbOBuL4Orvaxum4MMbvY2cKt/by+TdJpQTWTJGIQIAPDulmR8sesiAEAuFePOfj64s683hvX0gIej3My1I2o5aX2Lo0MHXrZanYDKWi0qrgtL6motyqvrxsmVVmlQWtnwff3X+q7k0vrxdKWVGn2gaghieaVVLa6DTCKuC0sOcng5yeHlWP/VSQ6F07XvvZzkcJJLGZqIboJBiPC/w5n6EDQ9PAAv3BkMXxc7M9eKyHJJxCI4yqVwlLfvT2it9tpEg9JKDa5W1NS91DUoqqjFVXWNfluRuhbFFTW4oq6pm4ig1UFVWg1VaTXO5TZfjlwq1gckHxdb+LrYwdfFFn6u1756OcrZqktWiUHIyh3PuIqFm08DAOLHBeOZsb3NXCMi62EjEcPdoa6rq6UEoa416mp9UCosr0ZheQ3yy6pQUFZ97VVe97WsSoNqjQ5ZVyuRdbXypueVikXwdraFn2t9UHK1RYCbPbq526Obhz38XO1gw7FM1AUxCFkxjVaHOWuPo0arw519vTFndC9zV4mIbkEkEsFeJoW9TAp/11u33FbValFQVo38smrkl1Yht6QKuSWVyCmpQm79DFBVaRU0OgHZxZXILq4EcLXReSRiEfxcbdHN3QHK+nAU6F736uZhDydbmw74aYk6HoOQFUsrVCO7uBL2Mgk+/PsgNosTdUG2NhIo3e2hdLe/6TEarQ75ZdV1Aam4Sv8162oFLl+pQEZRBao1OmQWVSKzqOlWJQ8HGXp6OaKnwhE9vRzQS+GInl6O8He1498WsmgMQlbsXG4pAKCPr3O7xzoQUecllYjh52oHP1c7hHdrvF+nE1BQXq0PRRlX1Lhc1PB9Ba6oa+pfRTh0qcjgs3Y2EgR5OaCnlyN6KxwR6uuMUB8nBLjZcQA3WQTe/azY2ZyGIORk5poQkSUT148f8na2RWQP90b7y6pqcamwAhcLynGxoByp+XWvS1fUqKzV4kxOKc7U/71p4CSXItTXCX18nRHq44w+vk4I8XGCvYy3JTItXnFW7Gx9i1BfXxcz14SIOjMnWxsMCHDBgADDvyUarQ4ZRRW4WKBGan45zqvKcC63FBcLylFWrcHhS1cN1i0TiYDuHg4Y4O+CgQEuCFO6op+fM8MRdSheXVZKEAS2CBFRh5JKxAjyckSQlyPG9fXWb6/R6JBWWI5zuaVIzi3D2dxSJOeVoaCsGumFaqQXqvHziRwAgFgEBHs7YWCACwYGuGJggAtCfZwhk3IGGxkHg5CVKiirxhV1DcQiINTH2dzVISIrIpOKEepT1yWGwde2F5ZX40xOKU5mFuNEVglOZhUjv6wayXllSM4rw/+OZOk/Hxbggoju7hjS3Q3hge5wseesNWobBiEr1dAt1sPTAXayTvYMBSLqkjwd5RgV7IVRwV76bXklVTiZVYyTWSU4Uf+1pLJW3632ef1xId5OiOjuhiHd3RHR3Q3+rhyMTS3DIGSlzl43Y4yIyFL5uNjCx8UHd/bzAVDXrZ9eqMaRy1dx5FIRjly6irRCNVJUZUhRlWHNwQwAgK+LLYYGeWBYL08M7+XB1fLpphiErFTD+KC+fgxCRNR5iEQi/bij+yKUAOq61I5cqgtGhy9fxZnsEuSWVOHH49n48Xg2ACDIywHDe9aFotuDPOBq3/LVvKlrYxCyUuf0M8YYhIioc/N0lGN8fx+M71/XalRRo8Gxy8XYe7EQ+1ILcSq7BGkFaqQVqPHdgcsQiYD+fi4Y1ssD0cEKRHR34+NDrBiDkBWqqNEgrVANgEGIiLoee5kUd/T2xB29PQEAJRW1OJB+BXtTC7E3tRAXC9Q4lV2CU9kl+HJXGpzkUowI9kR0iALRIV5QONma+ScgU2IQskIpeWUQBMDTUQYvJ7m5q0NE1KFc7G0Q288HsfXjjPJKqrDvYiH2XCjErvMFKFLX4PdTefj9VB4AYIC/C0aHKjA6xAthAa58REgXZ7Ftgbt378akSZPg5+cHkUiEzZs3G+x/+OGHIRKJDF7jx483OKaoqAgPPvggnJ2d4erqikcffRTl5eUm/Cksk6q0GgCgdLfnrAoisjo+LraYelsAPvr7IBx+JQab/jEMz4ztjQH+dQtCnsouwbKEC7jns30Y8tZ2vLjhBBLOqVBVqzVzzakjWGyLkFqtRlhYGB555BFMnTq1yWPGjx+PVatW6d/L5YatGw8++CByc3Px559/ora2FrNnz8YTTzyBtWvXdmjdLV2tVgcAkHNBMiKychKxCIMD3TA40A3x44KRX1aFXSkF2JmSjz3nC3FFXYMNR7Ow4WgWHGQSRIcqENvPB6NDvOBky7WLugKLDUITJkzAhAkTmj1GLpfDx8enyX3nzp3D1q1bcfjwYURERAAAPvnkE9x11114//334efn1+gz1dXVqK6u1r8vLS1tdExX0BCEODiQiMiQwskW90YocW+EErVaHQ6nF+GPM3nYdlaF3JIq/HYyF7+dzIVMIsawXh4Y388HMX294enIYQadVae+EyYmJkKhUCAkJARPPfUUrly5ot+3f/9+uLq66kMQAMTExEAsFuPgwYNNnu+dd96Bi4uL/qVUKjv8ZzCHhiAkYxAiIropG4kYw3p54vW7+2Pf/DH46enheCq6J4K8HFCj1SExpQDzfzyFyLe24/6vDmDNwcsoUteYu9rUShbbInQr48ePx9SpU9GjRw9cvHgRL7/8MiZMmID9+/dDIpEgLy8PCoXC4DNSqRTu7u7Iy8tr8pwLFixAfHy8/n1paWmXDEO1WgEAIJVwfBARUUuIRCKEKV0RpnTFP8eHIjW/DH+cUWHr6Tycyi7B/rQr2J92Ba/9dAbDe3li4kBfxPbzgYsdu88sXacNQjNmzNB/P2DAAAwcOBA9e/ZEYmIixo4d26ZzyuXyRuOMuiJ2jRERtU8vhRN6KZzw9OheyLpagd9O5uKXkzk4nV2K3ecLsPt8AV7ZdAqjgr0wcaAfYvp6w1HeaW+5XVqX+VcJCgqCp6cnUlNTMXbsWPj4+CA/P9/gGI1Gg6KiopuOK7IW7BojIjKeADd7/N+onvi/UT2RXqjGrydy8OvJXKSoyrD9XD62n8uHXCpGTB9v3DPYHyODvSDjZBWL0WWCUFZWFq5cuQJfX18AwNChQ1FcXIyjR48iPDwcALBjxw7odDpERUWZs6pm19A1xhYhIiLj6uHpgLlje2Pu2N44ryrDrydy8MvJXKQXqvHbqVz8dioXbvY2mBTmh3sG+2OQ0pXLmJiZxQah8vJypKam6t+np6cjKSkJ7u7ucHd3x+uvv45p06bBx8cHFy9exEsvvYRevXohNjYWANCnTx+MHz8ejz/+OL744gvU1tZizpw5mDFjRpMzxqxJjaa+a0zK//iIiDpKsLcT4u8MwfPjgnE6uxSbjmfj5xM5KCyvxrf7L+Pb/ZfRw9MBUwb5Y8pgP3TzcDB3la2SSBAEwdyVaEpiYiJGjx7daHtcXBw+//xzTJkyBcePH0dxcTH8/Pxw55134o033oC3t7f+2KKiIsyZMwe//PILxGIxpk2bhmXLlsHR0bFFdSgtLYWLiwtKSkrg7Nx1HkXx763J+CzxImYP745Fk/qZuzpERFZDo9Xhr9RCbD6ejT/OqFB53SKN4d3cMGWwPyYO8IWbAx8K2x6tuX9bbBCyBF01CL39+zl8tTsN/zcyCAvu6mPu6hARWaXyag22ncnDpuPZ2JtaCF393dhGIsLoEAWm3haAMaEKjidqg9bcvy22a4w6TkPXGKfPExGZj6Nciqm3BWDqbQFQlVbh56QcbDqejbO5pdh2VoVtZ1VwtbfB5DA/TL0tAGEBLhxP1AEYhKwQp88TEVkWb2dbPD4yCI+PDEJKXhl+PJ6FzcezoSq9Np4oyMsB024LwJTB/vB3tTN3lbsMBiErxCBERGS5QnycsGBCH7wUG4q9qYX48VgWtp7JQ1qBGkv+SMH721Jwew8PTAsPwPj+PlyfqJ3427NCDdPnuY4QEZHlkohFGBnshZHBXiiv1mDLqVz8eCxbv4r1/rQreHXzaYzv74Opt/ljWE9PSMTsOmstBiErVKNvEeJ/MEREnYGjXKp/GGzW1QpsPp6NH49lI61QjU3Hs7HpeDa8neWYMtgf024LQLC3k7mr3GkwCFmhWv06QmwRIiLqbALc7DFnTG88PboXkjKL8eOxuvWJVKXV+HJXGr7clYYB/i6Yeps/Jof5wcOx6z86qj0YhKyQpn6Opo2YQYiIqLMSiUQYHOiGwYFuWDixD3Ym52PjsWzsTM7HqewSnMouwVu/nUN0iJd+Kr6tjcTc1bY47Q5Ct912W6uOF4lE+Pnnn+Hv79/eoqmN9IOlubI0EVGXIJdKML6/L8b390WRuga/nMjBj8eycCKrRP+8M2dbKSbVT8W/LZCP9mjQ7iCUlJSEF154oUWrNQuCgHfffRfV1dXtLZbaQf+IDQ6WJiLqctwdZIgb1h1xw7ojNb8MPx6rG0OUW1KFNQczsOZgBrp72GPqbQG4Z7A/lO725q6yWbV7ZWmxWIy8vDwoFIoWHe/k5IQTJ04gKCioPcWaRFddWXrqZ3txLKMYX84MR2w/H3NXh4iIOphWJ+BA2hVsPJaFrafzUFFz7dEekT3cMe02f0wY4AtnWxsz1tJ4TLqydHp6Ory8vFp8/NmzZ63+oafmxunzRETWRSIWYXgvTwzv5Yk37tbgjzN5+PFYNvZeLMSh9CIcSi/Caz+dQWy/uqn4d/TyhNRK7hHtDkLdunVr1fFKpbK9RVI7cUFFIiLr5XDdoz1ySyqx+XgONh7LQmp+OX4+kYOfT+TAy0mOKYP8MGWwP/r6Onfp8URGnTV28uTJJreLRCLY2toiMDAQcjmn8ZlbLdcRIiIiAL4udngquieeHBWEU9kl+PFYNn5KykZBWTX+sycd/9mTjiBPB/xtoC/+NtAXId5OXS4UGfXp82KxuNlfkI2NDf7+97/jyy+/hK2trbGK7TBddYzQyH/vREZRBTY+NQzh3dzMXR0iIrIgNRoddp0vwI/HsrAjOR/V9RNsAKCnlwP+NtAPEwf6WvSija25fxu1b2TTpk3o3bs3vvrqKyQlJSEpKQlfffUVQkJCsHbtWqxYsQI7duzAwoULjVkstVJDixDHCBER0Y1kUjHG9fXG5w+F4+ir4/DxjEEY19cbMokYFwvUWJZwAXd+tBvjPtyFpdvP47yqDEZsUzE5o7YIRUZG4o033kBsbKzB9j/++AOvvvoqDh06hM2bN+OFF17AxYsXjVVsh+mqLUIRb/6JwvIabH1uBEJ9us7PRUREHaesqhbbz6nw28lc7D5fqH9cEwD08HTAuL7euLOvNwYHupn9mWcmnTV2vVOnTjU5eLpbt244deoUAGDQoEHIzc01ZrHUSlxHiIiIWsvJ1gb3DA7APYMDUFJZi+1nVfjtVC7+ulCI9EI1vtqdhq92p8HTUYaYPt4Y19cbw3t5Wvxq1kYNQqGhoXj33Xfx1VdfQSaTAQBqa2vx7rvvIjQ0FACQnZ0Nb29vYxZLrcTp80RE1B4udjaYFh6AaeEBKK/WYFdKAbadzcOO5HwUltfgv4cz8d/DmbCXSTAq2AsxfbwxKsQLnhb43DOjBqHly5dj8uTJCAgIwMCBAwHUtRJptVr8+uuvAIC0tDT84x//MGax1EqcPk9ERMbiKJfqZ5XVanU4mFaEbWfz8OdZFXJLqrDldB62nM6DSAQM9HdBdIgCo0MVGOjvArGZu9AAI48RAoCysjKsWbMG58+fBwCEhITggQcegJOT5Y4uv5muOEZIEAT0WPA7AODIwhiLTOdERNT5CYKA09ml+paiMzmlBvs9HGQYGeyF6BAvxPbzMWoXWmvu30YPQl1JVwxCNRodghduAQCcWHQnXOy6xnLqRERk2VSlVdiVUoCdKfn460Ihyqo1AOrWtDv+2p1wlBuvk8psg6UB4LvvvsOXX36JtLQ07N+/H926dcNHH32EoKAg3H333cYujlqp9rpR/hwjREREpuLtbIv7hihx3xAlarU6HL18FTtT8lFaqTFqCGoto94JP//8c8THx2PChAm4evUqtNq6h7q5ublh6dKlxiyK2uj6IMSVpYmIyBxsJGLcHuSBBRP64J2pA8xaF6MGoU8++QT/+c9/8Morr0AqvZbuIiIi9NPnybwa1n0QiWD2dR6IiIjMzahBKD09HYMHD260XS6XQ61WG7MoaqOGqfM2kuYfh0JERGQNjBqEevTogaSkpEbbt27dij59+hizKGojDR+vQUREpGfU0Unx8fF4+umnUVVVBUEQcOjQIaxbtw7vvPMOvv76a2MWRW3UMEZIyvFBRERExm0Reuyxx/Dee+9h4cKFqKiowAMPPIDPP/8cH3/8MWbMmNGqc+3evRuTJk2Cn58fRCIRNm/ebLBfEAS89tpr8PX1hZ2dHWJiYnDhwgWDY4qKivDggw/C2dkZrq6uePTRR1FeXt7eH7NTq9Fc6xojIiKydka/Gz744IO4cOECysvLkZeXh6ysLDz66KOtPo9arUZYWBiWL1/e5P5///vfWLZsGb744gscPHgQDg4OiI2NRVVVlUFdzpw5gz///BO//vordu/ejSeeeKLNP1tXwCfPExERXdNhE/ft7e1hb2/f5s9PmDABEyZMaHKfIAhYunQpFi5cqF+b6Ntvv4W3tzc2b96MGTNm4Ny5c9i6dSsOHz6MiIgIAHWz2u666y68//778PPza3PdOrNrj9dg1xgREVG7g9DgwYNbPPvo2LFj7S0OQN3stLy8PMTExOi3ubi4ICoqCvv378eMGTOwf/9+uLq66kMQAMTExEAsFuPgwYO45557Gp23uroa1dXV+velpaWNjunsavicMSIiIr12B6EpU6bov6+qqsJnn32Gvn37YujQoQCAAwcO4MyZM0Z90GpeXh4ANHqKvbe3t35fXl4eFAqFwX6pVAp3d3f9MTd655138Prrrxutnpbo+unzRERE1q7dQWjRokX67x977DE888wzeOONNxodk5mZ2d6iOtyCBQsQHx+vf19aWgqlUmnGGhmfhl1jREREekZtFtiwYQNmzZrVaPtDDz2EjRs3Gq0cHx8fAIBKpTLYrlKp9Pt8fHyQn59vsF+j0aCoqEh/zI3kcjmcnZ0NXl1NLbvGiIiI9Ix6N7Szs8PevXsbbd+7dy9sbW2NVk6PHj3g4+ODhIQE/bbS0lIcPHhQ3yU3dOhQFBcX4+jRo/pjduzYAZ1Oh6ioKKPVpbOpYdcYERGRnlFnjT333HN46qmncOzYMURGRgIADh48iJUrV+LVV19t1bnKy8uRmpqqf5+eno6kpCS4u7sjMDAQzz33HN5880307t0bPXr0wKuvvgo/Pz/9mKU+ffpg/PjxePzxx/HFF1+gtrYWc+bMwYwZM6x2xhgA1GrqW4SkDEJERERGDULz589HUFAQPv74Y3z//fcA6gLJqlWrcN9997XqXEeOHMHo0aP17xvG7sTFxWH16tV46aWXoFar8cQTT6C4uBh33HEHtm7datDytGbNGsyZMwdjx46FWCzGtGnTsGzZMiP8pJ3XtXWEOEaIiIhIJAiCYO5KWKrS0lK4uLigpKSky4wX+m7/Jbz60xlM6O+Dzx8KN3d1iIiIjK4192/2j1gZTp8nIiK6pt13Q3d3dxQWFrb4+MDAQFy+fLm9xVIb8aGrRERE17R7jFBxcTG2bNkCFxeXFh1/5coVaLXa9hZLbcRnjREREV1jlMHScXFxxjgNmQCnzxMREV3T7iCk0+mMUQ8yES6oSEREdA3vhlbm2jpCHCNERETEIGRlNLq6rjGOESIiImIQsjo1DbPGxPynJyIi4t3QyrBrjIiI6BoGISvD6fNERETXGP1uePHiRSxcuBD3338/8vPzAQBbtmzBmTNnjF0UtQFXliYiIrrGqHfDXbt2YcCAATh48CB+/PFHlJeXAwBOnDiBRYsWGbMoaqMaTp8nIiLSM+rdcP78+XjzzTfx559/QiaT6bePGTMGBw4cMGZR1EbX1hHiGCEiIiKjBqFTp07hnnvuabRdoVC06nlk1HE09V1jMilbhIiIiIx6N3R1dUVubm6j7cePH4e/v78xi6I24vR5IiKia4x6N5wxYwb++c9/Ii8vDyKRCDqdDnv37sW8efMwa9YsYxZFbcSuMSIiomuMGoTefvtthIaGQqlUory8HH379sXIkSMxbNgwLFy40JhFURvpgxC7xoiIiIzz9PkGMpkM//nPf/Daa6/h1KlTKC8vx+DBg9G7d29jFkPtUKvhIzaIiIgaGDUINVAqlVAqldBqtTh16hSuXr0KNze3jiiKWolPnyciIrrGqHfD5557DitWrAAAaLVajBo1CrfddhuUSiUSExONWRS1Ua2OY4SIiIgaGDUI/fDDDwgLCwMA/PLLL0hLS0NycjKef/55vPLKK8YsitqooWuMLUJERERGDkKFhYXw8fEBAPz++++47777EBwcjEceeQSnTp0yZlHURuwaIyIiusaod0Nvb2+cPXsWWq0WW7duxbhx4wAAFRUVkEgkxiyK2qiG0+eJiIj0jDpYevbs2bjvvvvg6+sLkUiEmJgYAMDBgwcRGhpqzKKojdgiREREdI1Rg9DixYvRv39/ZGZm4t5774VcLgcASCQSzJ8/35hFURvV8hEbREREekafPj99+vRG2+Li4oxdDLWBTidAq+NgaSIiogZGD0JqtRq7du1CRkYGampqDPY988wzxi6OWqFh6jzAMUJERESAkYPQ8ePHcdddd6GiogJqtRru7u4oLCyEvb09FAoFg5CZNXSLAWwRIiIiAow8a+z555/HpEmTcPXqVdjZ2eHAgQO4fPkywsPD8f777xuzKCxevBgikcjgdf2A7KqqKjz99NPw8PCAo6Mjpk2bBpVKZdQ6dDa1mutbhBiEiIiIjHo3TEpKwgsvvACxWAyJRILq6moolUr8+9//xssvv2zMogAA/fr1Q25urv71119/6fc9//zz+OWXX7Bhwwbs2rULOTk5mDp1qtHr0Jk0zBgTiwCJmF1jRERERu0as7GxgVhcl60UCgUyMjLQp08fuLi4IDMz05hFAQCkUql+AcfrlZSUYMWKFVi7di3GjBkDAFi1ahX69OmDAwcO4Pbbbzd6XTqDGk6dJyIiMmDUO+LgwYNx+PBhAMCoUaPw2muvYc2aNXjuuefQv39/YxYFALhw4QL8/PwQFBSEBx98EBkZGQCAo0ePora2Vr+OEQCEhoYiMDAQ+/fvv+n5qqurUVpaavDqSvRT5xmEiIiIABg5CL399tvw9fUFALz11ltwc3PDU089hYKCAnz11VfGLApRUVFYvXo1tm7dis8//xzp6ekYMWIEysrKkJeXB5lMBldXV4PPeHt7Iy8v76bnfOedd+Di4qJ/KZVKo9bZ3DQNLUJcQ4iIiAiAkbvGIiIi9N8rFAps3brVmKc3MGHCBP33AwcORFRUFLp164b//e9/sLOza9M5FyxYgPj4eP370tLSLhWG+HgNIiIiQ12macDV1RXBwcFITU2Fj48PampqUFxcbHCMSqVqckxRA7lcDmdnZ4NXV9LQNSYVd5l/diIionYx6h1RpVJh5syZ8PPzg1QqhUQiMXh1pPLycly8eBG+vr4IDw+HjY0NEhIS9PtTUlKQkZGBoUOHdmg9LFnDrDE+XoOIiKiOUbvGHn74YWRkZODVV1/VP3i1o8ybNw+TJk1Ct27dkJOTg0WLFkEikeD++++Hi4sLHn30UcTHx8Pd3R3Ozs6YO3cuhg4darUzxoBr6wixa4yIiKiOUYPQX3/9hT179mDQoEHGPG2TsrKycP/99+PKlSvw8vLCHXfcgQMHDsDLywsA8NFHH0EsFmPatGmorq5GbGwsPvvssw6vlyXj9HkiIiJDRg1CSqUSgiDc+kAj+O9//9vsfltbWyxfvhzLly83SX06g4YxQgxCREREdYx6R1y6dCnmz5+PS5cuGfO0ZCQN0+e5jhAREVGddrcIubm5GYwFUqvV6NmzJ+zt7WFjY2NwbFFRUXuLo3Zo6BqTcowQERERACMEoaVLlxqhGmQK7BojIiIy1O4gFBcXZ4x6kAnUcrA0ERGRAaPcEXU6Hd577z0MHz4cQ4YMwfz581FZWWmMU5MRXVtHiF1jREREgJGC0FtvvYWXX34Zjo6O8Pf3x8cff4ynn37aGKcmI6rRsEWIiIjoeka5I3777bf47LPP8Mcff2Dz5s345ZdfsGbNGuh0OmOcnoxEo+MYISIiousZ5Y6YkZGBu+66S/8+JiYGIpEIOTk5xjg9GUktW4SIiIgMGOWOqNFoYGtra7DNxsYGtbW1xjg9GUktnz5PRERkwCgrSwuCgIcffhhyuVy/raqqCk8++SQcHBz023788UdjFEdtVMPp80RERAaMEoSamkL/0EMPGePUZEScPk9ERGTIKEFo1apVxjgNdTB1tQYAYGcjMXNNiIiILAObBqzIpStqAECgh52Za0JERGQZGISsyKXCCgBAdw+HWxxJRERkHRiErERFjQZ5pVUAgB6eDEJEREQAg5DVaGgNcrW3gau9zMy1ISIisgwMQlaiYXwQu8WIiIiuYRCyEumFdUGI3WJERETXMAhZiUuFbBEiIiK6EYOQldB3jXnam7kmREREloNByEqk1w+WDvJ0NHNNiIiILAeDkBUoq6pFYXk1ALYIERERXY9ByAo0TJ33dJTBydbGzLUhIiKyHAxCViCdU+eJiIiaxCBkBfQzxjh1noiIyACDkBW4xDWEiIiImsQgZAXYNUZERNQ0BiErkFlUN1i6mwdnjBEREV1Pau4KmMLy5cuxZMkS5OXlISwsDJ988gkiIyPNXa12qdHooCqtQnZxJTKuVOBiQTnyy6rh62KL7p4OiO3nAxc7G1TValFYXgMACHCzM3OtiYiILEuXD0Lr169HfHw8vvjiC0RFRWHp0qWIjY1FSkoKFAqFWeqk0wm4WlGDylotqmq1qKzRobJWW/eqqd9W/315tQZXK2pQUlGL4spaFKlrkFtSifyyagjCzcvYm1qIj2cMRnZxJQDAQSaBix2nzhMREV2vywehDz/8EI8//jhmz54NAPjiiy/w22+/YeXKlZg/f75Z6lRZq0X4m9vbfR6ZVAx/VzsEuNmhp5cjFM5ynMkpxW8nc3E6uwQAkHW1LggFuNlDJBK1u0wiIqKupEsHoZqaGhw9ehQLFizQbxOLxYiJicH+/fsbHV9dXY3q6mr9+9LS0g6pl62NBAAgl4phJ5PAzqbuZWsj0b9v+N5BJoGrvQxu9jZwtbeBq70MPs628Hezg4eDrFG4ybpagd9O5iKzqBJanYDs+iDkz24xIiKiRrp0ECosLIRWq4W3t7fBdm9vbyQnJzc6/p133sHrr7/e4fWSiEVIe/suiMXGb6HxdbGDTCJGjVaHnOJKZBfXDZT2d2UQIiIiuhFnjV1nwYIFKCkp0b8yMzM7rKyOCEFAXcgKrJ8ddumK+rquMQYhIiKiG3XpFiFPT09IJBKoVCqD7SqVCj4+Po2Ol8vlkMvlpqpeh+nu4YDU/HJcKlSza4yIiKgZXbpFSCaTITw8HAkJCfptOp0OCQkJGDp0qBlr1rG617cIpRdW6GeNBbhxDSEiIqIbdekWIQCIj49HXFwcIiIiEBkZiaVLl0KtVutnkXVFDc8Uu5BfhrzSKgAcI0RERNSULh+E/v73v6OgoACvvfYa8vLyMGjQIGzdurXRAOqupOGZYkcuXYUg1M1O83SUmblWRERElqfLByEAmDNnDubMmWPuaphMw6M0Kmu1AOrGB3ENISIiosa69Bgha+XnYgeZ9No/LbvFiIiImsYg1AWJxSJ0c782OJoDpYmIiJrGINRFdfNw0H/PNYSIiIiaxiDURfXwvNYKxK4xIiKipjEIdVENU+gBtggRERHdDINQF9Xjuq4xripNRETUNKuYPm+NeikcIRIB9jYSKJxszV0dIiIii8Qg1EUpnG2x/IHb4GJnA0kHPeCViIios2MQ6sLuGuBr7ioQERFZNI4RIiIiIqvFIERERERWi0GIiIiIrBaDEBEREVktDpZuhiAIAIDS0lIz14SIiIhaquG+3XAfbw6DUDPKysoAAEql0sw1ISIiotYqKyuDi4tLs8eIhJbEJSul0+mQk5MDJycniETNr8UzZMgQHD58uFXnb81nWnLsrY5pbn9T+0pLS6FUKpGZmQlnZ+cW1dMc2vK7N/X523qOln7OGMe1ZR+vEeOdn39DzKejrw9jlcFrpOUEQUBZWRn8/PwgFjc/CogtQs0Qi8UICAho0bESiaTV/4it+UxLjr3VMc3tb26fs7OzRf8Ra8vv3tTnb+s5Wvo5YxzX1n0ArxFzXSP8G2IcHX19GKsMXiOtc6uWoAYcLG0kTz/9dId+piXH3uqY5va3pf6WoqPrbozzt/UcLf2cMY7rqtcH0HWvEf4NMQ5T1J3XiOVeI+wao5sqLS2Fi4sLSkpKLPr/5sh8eI1Qc3h90K1YwjXCFiG6KblcjkWLFkEul5u7KmSheI1Qc3h90K1YwjXCFiEiIiKyWmwRIiIiIqvFIERERERWi0GIiIiIrBaDEBEREVktBiEiIiKyWgxC1Ca//vorQkJC0Lt3b3z99dfmrg5ZoHvuuQdubm6YPn26uatCFigzMxPR0dHo27cvBg4ciA0bNpi7SmRBiouLERERgUGDBqF///74z3/+02Flcfo8tZpGo0Hfvn2xc+dOuLi4IDw8HPv27YOHh4e5q0YWJDExEWVlZfjmm2/www8/mLs6ZGFyc3OhUqkwaNAg5OXlITw8HOfPn4eDg4O5q0YWQKvVorq6Gvb29lCr1ejfvz+OHDnSIfcZtghRqx06dAj9+vWDv78/HB0dMWHCBGzbts3c1SILEx0dDScnJ3NXgyyUr68vBg0aBADw8fGBp6cnioqKzFspshgSiQT29vYAgOrqagiCgI5qt2EQskK7d+/GpEmT4OfnB5FIhM2bNzc6Zvny5ejevTtsbW0RFRWFQ4cO6ffl5OTA399f/97f3x/Z2dmmqDqZSHuvEer6jHmNHD16FFqtFkqlsoNrTaZijOujuLgYYWFhCAgIwIsvvghPT88OqSuDkBVSq9UICwvD8uXLm9y/fv16xMfHY9GiRTh27BjCwsIQGxuL/Px8E9eUzIXXCN2Ksa6RoqIizJo1C1999ZUpqk0mYozrw9XVFSdOnEB6ejrWrl0LlUrVMZUVyKoBEDZt2mSwLTIyUnj66af177VareDn5ye88847giAIwt69e4UpU6bo9z/77LPCmjVrTFJfMr22XCMNdu7cKUybNs0U1SQzaus1UlVVJYwYMUL49ttvTVVVMoP2/A1p8NRTTwkbNmzokPqxRYgM1NTU4OjRo4iJidFvE4vFiImJwf79+wEAkZGROH36NLKzs1FeXo4tW7YgNjbWXFUmE2vJNULWrSXXiCAIePjhhzFmzBjMnDnTXFUlM2jJ9aFSqVBWVgYAKCkpwe7duxESEtIh9ZF2yFmp0yosLIRWq4W3t7fBdm9vbyQnJwMApFIpPvjgA4wePRo6nQ4vvfQSZ4xZkZZcIwAQExODEydOQK1WIyAgABs2bMDQoUNNXV0yg5ZcI3v37sX69esxcOBA/fiR7777DgMGDDB1dcnEWnJ9XL58GU888YR+kPTcuXM77NpgEKI2mTx5MiZPnmzuapAF2759u7mrQBbsjjvugE6nM3c1yEJFRkYiKSnJJGWxa4wMeHp6QiKRNBqUplKp4OPjY6ZakSXhNUK3wmuEmmNp1weDEBmQyWQIDw9HQkKCfptOp0NCQgK7NQgArxG6NV4j1BxLuz7YNWaFysvLkZqaqn+fnp6OpKQkuLu7IzAwEPHx8YiLi0NERAQiIyOxdOlSqNVqzJ4924y1JlPiNUK3wmuEmtOpro8OmYtGFm3nzp0CgEavuLg4/TGffPKJEBgYKMhkMiEyMlI4cOCA+SpMJsdrhG6F1wg1pzNdH3zWGBEREVktjhEiIiIiq8UgRERERFaLQYiIiIisFoMQERERWS0GISIiIrJaDEJERERktRiEiIiIyGoxCBEREZHVYhAiIiIiq8UgRESdjiAIeOKJJ+Du7g6RSISkpKRGx6xevRoikQgikQjPPfec0evQcG5XV1ejn5uITIdBiIg6na1bt2L16tX49ddfkZubi/79+zd5nLOzM3Jzc/HGG2/ot0VHRzcZjFavXt2qUJObm4ulS5e2suZEZGn49Hki6nQuXrwIX19fDBs2rNnjRCIRfHx8OqQOPj4+cHFx6ZBzE5HpsEWIiDqVhx9+GHPnzkVGRgZEIhG6d+/eIeV0795d3/11/YuIuha2CBFRp/Lxxx+jZ8+e+Oqrr3D48GFIJJIOKefw4cPQarUAAK1Wi+nTp8PGxqZDyiIi82EQIqJOxcXFBU5OTpBIJG3u9vrss8/w9ddfG2zTaDSwtbXVv/fy8tJ//+yzzyI3NxeHDx9uW6WJyGIxCBGR1XnwwQfxyiuvGGz78ccf8fbbbzc69quvvsKKFSuwb98+g3BERF0DgxARWR0XFxf06tXLYJtCoWh03M6dOzF37lysW7cOAwcONFX1iMiEOFiaiKgJqampmD59Ol5++WVMnTrV3NUhog7CFiEiohtUVlZi0qRJGDx4MJ544gnk5eXp93XUdHwiMg8GISKiG6hUKiQnJyM5ORl+fn4G+wRBMFOtiKgjiAT+V01EXdDq1avx3HPPobi4uFOXQUQdi2OEiKjLKikpgaOjI/75z38a/dyOjo548sknjX5eIjIttggRUZdUVlYGlUoFAHB1dYWnp6dRz5+amgoAkEgk6NGjh1HPTUSmwyBEREREVotdY0RERGS1GISIiIjIajEIERERkdViECIiIiKrxSBEREREVotBiIiIiKwWgxARERFZLQYhIiIislr/D8QPAAOBzI6BAAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"virgo.plot_DARM()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "870833dc-4b69-403d-9d50-2c187e66750d",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"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.2"
}
},
"nbformat": 4,
"nbformat_minor": 5
}