finesse-simulation-O4/generate_from_finesse-virgo.ipynb
2025-03-27 14:49:09 +01:00

104 lines
28 KiB
Text

{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"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"
]
},
{
"cell_type": "code",
"execution_count": 2,
"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": 3,
"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": 4,
"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",
")"
]
}
],
"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
}