[15:50:25] WARNING: Simulation final field decay value of 0.000136 is greater than the simulation shutoff log.py:50\n",
"threshold of 1e-05. Consider simulation again with large run_time duration for more accurate \n",
"results. \n",
"
\n"
],
"text/plain": [
"\u001b[2;36m[15:50:25]\u001b[0m\u001b[2;36m \u001b[0mWARNING: Simulation final field decay value of \u001b[1;36m0.000136\u001b[0m is greater than the simulation shutoff \u001b]8;id=966988;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/log.py\u001b\\\u001b[2mlog.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=288417;file:///Users/twhughes/Documents/Flexcompute/tidy3d-docs/tidy3d/tidy3d/log.py#50\u001b\\\u001b[2m50\u001b[0m\u001b]8;;\u001b\\\n",
"\u001b[2;36m \u001b[0mthreshold of \u001b[1;36m1e-05\u001b[0m. Consider simulation again with large run_time duration for more accurate \u001b[2m \u001b[0m\n",
"\u001b[2;36m \u001b[0mresults. \u001b[2m \u001b[0m\n"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"batch_data = batch.load(path_dir=\"data\")\n",
"flux_norm = batch_data[\"normalization\"][\"flux\"].flux\n",
"trans_g0 = batch_data[\"Si-resonator-delta-0\"][\"flux\"].flux / flux_norm\n",
"trans_g20 = batch_data[\"Si-resonator-delta-20\"][\"flux\"].flux / flux_norm\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The normalizing run computes the transmitted flux for an air -> SiO2 interface, which is just below unity due to some reflection.\n",
"\n",
"While not technically necessary for this example, since this transmission can be computed analytically, it is often a good idea to run a normalizing run so you can accurately measure the *change* in output when the structure is added. For example, for multilayer structures, the normalizing run displays frequency dependence, which would make it prudent to include in the calculation."
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"tags": []
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAikAAAGjCAYAAADpUIU9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAACDk0lEQVR4nO3deZxN9f/A8ded3VhmsYwla4pE9iQVMVHJ+vtGZAkplawJ2WVXaBFR2gulJEohRERhxpIo20jG2GYwYxYz5/fHZ+42c2fm3nHvPffeeT8fj/s455577r2fOXPvOe/7/mwGTdM0hBBCCCE8jJ/eBRBCCCGEsEWCFCGEEEJ4JAlShBBCCOGRJEgRQgghhEeSIEUIIYQQHkmCFCGEEEJ4JAlShBBCCOGRJEgRQgghhEeSIEUIIYQQHkmCFCGEEEJ4JF2DlF9++YUOHTpQsWJFDAYDq1evLvA5W7ZsoVGjRgQHB1OzZk0+/PBDl5dTCCGEEO6na5CSnJxM/fr1WbhwoV37nzhxgvbt2/Pggw8SExPDsGHDePrpp/nxxx9dXFIhhBBCuJvBUyYYNBgMfPPNN3Tu3DnPfUaPHs26des4ePCgadsTTzxBYmIi69evd0MphRBCCOEuAXoXwBE7d+4kOjraalu7du0YNmxYns9JS0sjLS3NdD8rK4tLly5RunRpDAaDq4oqhBBC+BxN07h69SoVK1bEz8/1lTFeFaTEx8cTFRVltS0qKoorV65w/fp1ihUrlus5M2fOZMqUKe4qohBCCOHzTp8+zS233OLy9/GqIKUwxo4dy4gRI0z3k5KSqFKlCqdPn6ZUqVI6lsy7vPUWjB+v1v/7D4oX17c8Pumee+DwYbV+yy0cW3OIRo3MDz/wAHz3nT5FE0IIgCtXrlC5cmVKlizplvfzqiClfPnynDt3zmrbuXPnKFWqlM0sCkBwcDDBwcG5tpcqVUqCFAdYZvVKloQSJfQri8+6csW8npqKv7/15zM9HeQjK4TwBO5qLuFV46Q0b96cTZs2WW3bsGEDzZs316lERUdGhnk9K0u/cvi0xETzekoKKSnWD+e8L4QQvk7XIOXatWvExMQQExMDqC7GMTExxMXFAaqqpk+fPqb9Bw0axPHjx3n55Zf566+/eOedd1i5ciXDhw/Xo/hFSnq6ed0z+oP5mIwMSE42309J4XqK9YG+ft3NZRJCCJ3pGqT88ccfNGzYkIYNGwIwYsQIGjZsyMSJEwE4e/asKWABqF69OuvWrWPDhg3Ur1+f119/nffee4927drpUv6iRIIUF7PMomRLSVQH3VjVJpkUIURRo2ublFatWpHfMC22RpNt1aoV+/btc2GphC1S3eNiSUlqGRpqikZUkBJMmTKQkCBBihCi6PGqNilCP5JJcbGrV9UyPByCggBzJiU8XD1k+T8QQoiiQIIUYZfMTPO6BCkucO2aWpYoobIpwPWrKn0lQYoQoqiSIEXY5cYN87pU97iAMZNSsqQpSEm5og66MUjJyJAAUQhRtEiQIuxi2SZFLpQuYCOTknJFpa/Cwsy7Wf4fhBDC10mQIuwimRQXs5VJuaYOtGWQIlU+QoiiRIIUYRfJpLiYZSYle4Rk4zgpxuoekCBFCFG0SJAi7CJBiosZMykWQUpK9uBtJUuCcQRqqe4RQhQlEqQIu8g4KS5mzKSULGnugnxdRSbFi5s2SSZFCFGkSJAi7GLZJkUyKS5gWd1jClLU17NYMQlShBBFkwQpwi5S3eNilg1njW1S0lQmJTQUAgPVwxKkCCGKEglShF2kusfFbGVSUv0BFaRIJkUIURRJkCLsItU9LmbZcNYYpKRJkCKEKNqKbJCSeFmutI6Q6h4Xs2w4a+zdkx2kWLZJkd49QoiipMgGKVWrGZj6wjmrOWlE3qS6x8VsZFKuZ0gmRQhRtBXZIAVg0jtRBATAO6/8y+XLepfGs0l1j4vZ6oKcHgBIkCKEKLqKdJBi9MLMW6hR7io/LDsrF+A8SCbFxWwN5pahuvRIkCKEKKqKbJDy1Vdwx23mK2/ijZI8OqACHasfIOX0RR1L5pmkTYqL2erdc0MtixWTLshCiKKpyAYpDz0Efx4N5NAh6+1rT9WjeJXSfPjiHjIzJGVgJEGKC2VkQFqaWs+u7skggBtZ0iZFCFG0FdkgxahOHdXeYu0a64Ck39uNCQrSuPbfFZ1K5llkFmQXMmZRwFTdc51ipk2WQYr07hFCFCVFPkgB8PeH9h38uHQJevcwXwWy8KdkpVJM6RJT5LMHkklxoeRktQwMVNFIUBAphAJqYsHgYMmkCCGKJglSLEREwMefB3LsmPX2yasbcFeZM/y5/ZI+BfMAEqS4kDGTUry4WloEKcWKqUBF2qQIIYoiCVJsqFFDXZTfecs8iMrBS5W48/5I4lb9rmPJ9CPVPS5kzKSUKKGWFtU9oSpWIUD1RpZxfYQQRYoEKXkICIDnBvvz009QPMR8ha76v6Z0r3+4yDWqlUyKC+WTSckZpFgGi0II4W6LF7v3/SRIKcBDD8HFxAB6djNfpVfuv4Oet/9OSlLRaMWoaRKkuFTOTIoEKUIID3PuHLRrB6NHu/d9JUixQ3AwfLYikF+2mq/OK082o3h4IPvXnNSvYG6Ss3pHqnuczHKMFIDAQKs2KaAad4MEKUII99I0ePddqFwZfvrJ/e8vQYoD7n/AwL//wt23m8fQr9+pGke++VPHUrlezm6vkklxspzVPQEBebZJkSBFCOEu58/Dww/DoEHqOlCyeCarKr3o1jJIkOKgSpXg10MRVL3F3IKxdtc6TO6wR8dSuVZKivV9yaQ4Wc7qHotMigQpQgg9rFgB1aubsye96u/n7PUIos987NZySJBSCAEBcPK0P8uXma/eU9Y2ZtH/NupYKtfJGaRIJsXJbGRS8gpSpHePEMKVUlOhb1944gn1+ykyPJN1dUfzSWx9imddhS5d3FoeCVJuQvd+oXz1pfmK/fyqaG4tlUDGdd/6uWv8oW8kQYqT5cykWAQpxjYpkkkRQrja3r1w++3wcXay5PF7z3CcW3n04Bx1fvr4Y/jwQ7eWSYKUm/R//zMQG2O+ah+/Wo6g0ACfClRyBilS3eNkOTMpgYHSJkUI4TaaBvPmwd13w+nTEBKi8UHrT1i54xbCEk9Bkyawbx/07u32skmQ4gR31TeQmGi97fEqu0g9e9nm/t5GqntcLJ9MigQpQghXSkqC//s/GDlSVSc/1vIqR6s8xFM/91E7jBoFv/4KNWvqUj4JUpwkLAz++898/9sLLShWMYJjh1L1K5STSHWPiznQJkWCFCGEs+zbp5Ik33wDgYEaC3tsZ81v5ah8dBNERcGPP8KcOebJw3QgQYoTVaigqkLmv3TGtK1m3RBO/n5ex1LdPKnucTE72qTIOClCCGf64ANo3hz++QeqVNbY/vB0nv/ifgxpqWrUtthYaNtW72JKkOJsBgMMmV2JZnWumLY1vxcunUjSsVQ3R6p7XEzapAgh3CQ9Hfr3V7e0NGjfJpW95R7m7u8mgJ8fzJoF33+vMikeQIIUF/Dzg1/3l6JFo+sAxN8oS+kaYfz5q3e2UZHqHheTNilCCDeIi1PVOx98oO6P6X2G1QdrUnrPTxARAT/8oMa99/Oc0MBzSuJj/P1h+55ifDPvhGnbnfdF8Mt33pdRkeoeF8s5LL4EKUIIJ/v1V2jUCA4cUE1MVjz7MzNX1CDg3BmoWxd+/90jqndykiDFxToPr847E8+a7rfsGMaRvcn5PMPzSCbFxYwH2KK6J69xUmQwNyGEoz79FFq3hosXoXJljd87T6fbu21U3U/XrrBzJ9x6q97FtEmCFDd4bkoFnut+yXS/+d03uPBPon4FcpDxh76RZFKczEYmRdqkCCFulqbBq6+q4U3S0+G+Zun8UbETd60crxpQvvoqfPml+dzjgSRIcZN3lkcypFs8AJczw3i0cTxJF73jinP1qvV9yaQ4Wc5MimV1TzHNuAmQIEUIYZ+MDNU4duJEdX9Qt4ts+q8O5XZ9B6VKwbffwvjxHtX+xBbPLp2PeWNFeb6ae4IQrvP7ldqElwng1D8ZBT9RZxKkuFBmpposA2xPMBis6nekC7IQwl6XL0ObNmoEe39/ePvFIyxaX4Og08fUuPe7dkGHDnoX0y4SpLjZ/71UnTWv7jfdr3ZbIEePePZVP2eQItU9TmTZ4MdGJqVYUKZxEyBBihAif8ePqx4827ap3z2rh/zMC4vqwpUrcP/9qv1J7dp6F9NuEqTo4KHxzZjW+4jpfq3aBv79V8cCFUAyKS5kbI/i7w/BwWrdsk1K0A3jJkCCFCFE3vbuhWbNVKASFQWb+3/CY/PbqBPHE0/ATz9BZKTexXSIBCk6GfdxLZ5o8o/pfuXKuQdN8xRXrljflyDFiSzboxgMat2yuidIVQdKkCKEyE9MDDz4IFy4ADVqaPzy2ByavJk9/86IEfDZZxASomsZC0OCFB198XtNWlQ9bbrfrvF5j+xiKtU9LpSzZw+QiT9pqJNJaKBkUoQQ+Vu/XtXkXLkCjRtlsb3xMG5/f7R6cM4ceP11j28gmxfvLLUP2X6yMq2r/K3W/yrLxH6nC3iG+xmDFOMcU5JJcaKcPXuA62nmr2WxHEGKJwaxQgj9rFwJHTuq3zutW2aysUwPKnz5pqpC/uADNYuxF5MgxQNs/Kc6PStsBmDGJ5VZNM6zGqgYg5SwMLWUIMWJbGRSrl83P1zMPx2QTIoQIreFC6F7d9Xd+PFOafyQ+iDhP61U1TrffANPPaV3EW+aBCkewBAYwKd/N2Nopa8AeHFGeb768GoBz3IPTTNfR41BilT3OJGNTIqxbVII1/HLkuoeIURuc+fC4MFqvX/PVD4/3pygXdsgPBw2bvSaLsYFkSDFQxiKhzI/tjV9i39FJgH06h/Ej9+m6l0skpPNmZNSpdRSMilOZCOTYgxSQkkxRSUyTooQwui11+Dll9X6c/2us2T/PQQc2Afly6u+xy1a6FtAJ5IgxYMYSkfy/vZatAnYQpoWzMOdQ/h9t74RgbGqx8/P/GPfVzIpFy+qNmX//adjIXLOgIy5uqcY11UeF8mkCCHUD8RZs8zNTEY8d52Ff9yD/8FYFaBs3qwmC/QhEqR4GP8G9fh2fQi3cxSAu5sZ2LtXv/IYg5QSJcyNw30lk/LUU2pW8jZtdCyEMZNio7rHMpMiQYoQRZsxQBk7Vt2f/HIKr/3aHMOB/dmDomz2qkHa7CVBigcq3uYeNs/8jUBUo8nGjdXgPHowjpFSsqTvBSk//aSWf/2lYyFsZFIkSBFCWNI0mDkTXnlF3X91bAqTfrwXw/5YFaBs2eKTAQpIkOKxKo7uzZbui033a9bUOH/e/eUwZlJKljSPNeYr1T3FiuldAiSTIoQo0BtvwLhxan3M0OuM+74FxMb6dAbFSIIUT2UwcO/ng/mk8QIANM3AXXUz3T4qra0gxVcyKaGhepeAfBvOSpsUIcSSJTB8uFp/afB1ZmxtgSE2BsqVg59/hjvu0LV8riZBiifz86PXxqeYGP4mAPEJ/rSNznTrgF6WQYqvVfd4RJBiowuyaRPJuTIpMpibEEXHp5/Cs8+q9ZGDU5mzvQWGmH0qQNm8GerU0beAbiBBiqcLD2fKrod5MuhLAH7d6c+YMe57e2OQUqqUVPe4RD6ZFFtBimRShCgavv4a+vZV60Oev8HcfQ+ZA5Sffy4SAQpIkOIdbr+dj38oy4P8DKg+8h995J63PntWLcuWleoelygok5Jd3SPjpAhRdGzaBD16qB+ET3TPYt7pxzH8ul2NqLlhA9x5p95FdBsJUryEX+tWfD9uB435A4CnB2Tx22+uf99/sidqrlnT96p7PDWTYgxSpOGsEEXP3r3QtSukp8Ojj2p8FDII/+9WQ3AwrFkDd92ldxHdSoIULxLy6ji2PP4OTdnNjUw/Hn0ki7g4176nZZDia9U9XpFJkSBFiCLjyBE1mv2VK2pW46/qTCLoo6XqF+KKFfDAA3oX0e10D1IWLlxItWrVCAkJoVmzZuzevTvf/RcsWECtWrUoVqwYlStXZvjw4aSm6j98vFsYDJT4+B2+q/0ydTjE5UQ/2rXVuHzZdW958qRaVq8umRSXKKhNivTuEaJI+O8/aN9eLevUge8eXUSx115VD777LnTqpG8BdaJrkLJixQpGjBjBpEmT2Lt3L/Xr16ddu3YkJCTY3P/zzz9nzJgxTJo0icOHD/P++++zYsUKXjGOcFMUhIQQ9dVCvgt9gls4zV9HDPToobmk10dmJly4oNYrVJBMikvkk0mxVd2Tmek7QaIQQklOVjHIsWNQpQpsHryKsLHPqwenT4enn9a3gDrSNUiZN28eAwcOpF+/ftSpU4fFixcTGhrKsmXLbO6/Y8cOWrRoQc+ePalWrRpt27alR48eBWZffM6dd1Jj5Sy+M3SiGCn8+KOBAQNMP7qd5uJFFZAYDFCmjG83nNWta69l96ls+VX3gHRDFsKXZGbC44/DH39ARARsmryNckOeUA8OGWIeB7+I0i1ISU9PZ8+ePURHR5sL4+dHdHQ0O3futPmce++9lz179piCkuPHj/P999/z6KOPuqXMHqV9exq8+n98yFOA6u3z2GPOvYCdO6eWpUuri6SvVfeEhJjXjZP6uVVGhvmNLYKU/Logg1T5COErNA0GDYIfflDVz+vm/knNFx9RX/IePWD+fPOvwyJKtyDlwoULZGZmEhUVZbU9KiqK+Ph4m8/p2bMnU6dO5b777iMwMJBbb72VVq1a5Vvdk5aWxpUrV6xuPmPMGLpFX2YuLwFqLpr27Z0XRBiDFOO/yNeqe3TPThizKKBGy8tmVd2TowsySCZFCF8xeza89576fn/0+nmav/KgOgG0bQsffmj+ZViEedUR2LJlCzNmzOCdd95h7969fP3116xbt45XX301z+fMnDmTsLAw061y5cpuLLGL+fvD558zsuoqHmUdAD/+CM8845yX/+8/tcwZpPhKJsXywq9LdsIYMBcrZhUxFVTdI5kUIbzft9+aa3JmTk7j8XcfgoQEqF8fVq2CoCB9C+ghdAtSypQpg7+/P+eMP9eznTt3jvLly9t8zoQJE+jduzdPP/009erVo0uXLsyYMYOZM2eSlcfP+7Fjx5KUlGS6nT592ul/i67KlsWw6iu+DuhOWVSD4/feg7feuvmXjolRy7p11dLXqnsss6i6ZCeMQYpFVQ9IkCKEr9u/X9XmAAx6VmPUnifMEwauWWPV26+o0y1ICQoKonHjxmzatMm0LSsri02bNtG8eXObz0lJScEvR/rLP/vnsJbHlTM4OJhSpUpZ3XxO48YEz5zMUW43bRoyBPJo2mO3P/4wvTzge9U9nhqk2OqCbPmxlyBFCO915gw8/LBqjta2LbwVPgFWr1aDtX3zjereI0x0re4ZMWIES5cu5aOPPuLw4cM899xzJCcn069fPwD69OnDWIuWzR06dGDRokUsX76cEydOsGHDBiZMmECHDh1MwUqRNWIE4V3b8B8VKO13CVBfhMOHC/dymZlq5EOAJk3U0tcyKZZ/h67VPXlkUiy7IBsMMlaKEN4uLQ169lTTjdx5J3zReQUBs6erB99/H/L4gV6UBRS8i+t0796d8+fPM3HiROLj42nQoAHr1683NaaNi4uzypyMHz8eg8HA+PHjOXPmDGXLlqVDhw5Mnz5drz/Bc/j5wUcfUSG2AQeP1aF9xA72Xq7BQw+puahuv73gl7B05Ii6WBYvDrVqqW2+lkmx5EmZFFvVPaCClBs3JEgRwhtpGgwcCL/8otrJfzVuH5FP9VEPjh0LTz6pbwE9lK5BCsDgwYMZPHiwzce2bNlidT8gIIBJkyYxadIkN5TMC5UoAZ99RvkWLVh/uRktyh7l7zMRNGsGu3fDbbfZ/1LGqp6GDc0NTH2t4azl3+EpQUpGhjkIsazuAcmkCOHN3noLPvlEnU+/WniO2sMeVhP0dO4M06bpXTyP5VW9e4QdmjWD2bMpywU2XGpCqeKZJCaqTIojVT/GIMVY1QO+V91jyVOqe4xZFLCu7gEJUoTwVjt2wMiRav31GWm0fa2tuSfPJ59IV+N8yJHxRSNGwP/+R9XM42wv28W0uWlTVY1TkKtX4csvzc8x8rXqHk/MpBiDFH9DJkGkW0UkxoyWjJMihPdISIBu3dRXudvjGkN+7a6690hPHrtIkOKLDAZYvBgqVKDeye841HkcpUtrJCermTW3bcv/6WvXQny8amTeubP1y/oqTwlSTD17AjMwgGRShPBimZmqoeyZM6pt33u15mJY863qybN6tfTksYMEKb6qdGlTGrHO6hnsH7uc+vXh/Hk12/frr9uuttE0+Owztd6tm+1J+KS6x0nyyaSEBma3RZE2KUJ4rcmTYdMmdR5d9dJOSk4fox5YtAjuuUfXsnkLCVJ8WZs2MGMGABUnP8OOL06Zeri99JIKQo4ds37K8uWwbp26IPbqZf2Yr2VSPKa6x8aQ+MWDsoMTyaQI4ZW+/97cHnbprIvcObajOukMGADZw2yIgkmQ4uteeglatIBr1wgd+CTbf05n6FD10FdfqQa13bvDxx+rpiw9e6rHxo9XbbqKCk+p7jEFKcESpAjhrU6eNP/Ie/7ZTHp+0QEuXIAGDZwzHHgRIkGKr/P3V1MklyoFv/6K3+hRLFgAv/4Kjz6qGsGuXAl9+6oJNwGioyGfORt9prpH90yKcYJBiyDFuKlksFT3COGNUlPhf/+Dy5fh7rthXuBoNfx3WJiak6dYMb2L6FUkSCkKbr3V3NDkzTfhxx+5915VrbN/Pzz/vEq29OoFK1ao2ZQDA3O/jK9V91jylDYppk3F0tWKZFKE8CrDhsGePRAZCSv7riP47dfVAx9/DDVq6Fo2b6T7YG7CTR57DF58UaUa+/WDAwegdGnq1YOFCx17KcmkOImNIMWUXCmWHYnY6IIsQYoQnumTT+Ddd9UPus9mn6bq8CfUA6NHQ8eO+hbOS0kmpSiZPRvuuENNHDFggMPRhi9nUjwlSDG1pTUGKVLdI4RXOHAAnn1WrU8cm8HDCx6Ba9egVSsZUfYmSJBSlBQrpqp9goLg229VyF8IvpJJseT2C39WVr5tUkqF5s6kGIMUGcxNCM9y5Qr83/+ZZzaecH4IHDoE5cvDF1+Yv7zCYRKkFDUNG8KsWWp95Ej4+2+7n+prmRRdq3uSk80FsJVJCc3MVTDJpAjheTQN+vdXp9LKleGzHmvxX7o4u87nMxWoiEKTIKUoGjoUWrdWw5v26mVVpWAPX8ykuD1IMUYjAQEQEpJrc36ZFAf/XUIIF1qwQHXaCQyELxcmUGZE9szGo0er86y4KRKkFEV+fvDhhxAerqZHHj/erqf5cibF7dkJy4HcLA6sqbqneGaugsncPUJ4lpgYFYsAzHsti2bzuqu+x02bwtSpupbNV0iQUlRVrgzvv6/W586F7dvtfqpkUpzARqNZy80li2fP4mijukeCFCH0l5pqTkR37gwvXJsNW7ZA8eLw+ee2x3EQDpMgpSjr2hWeekpFHU89pVqi58PXMimWPCVIMWVSSmQHKdIFWQiPNHq0ahsbFQVLnt2DYeIE9cDChVCzpr6F8yESpBR1CxaorMqxY2pUNzv4SiZF14az9mZSpHePEB7nhx/UuJgAHyxMoewL3dQX84knoE8ffQvnYyRIKerCwlQLdH9/NRLRl1/muasvZ1Lcnp24fFktw8OtNpsyKSWzIyiLiEQyKULoLyFBJZ4BhgyBR74dBMePQ7VqsHixb58odSBBioD77zdP1vPCC3D+vL7lcRNdMynGICUiwmpzrgSLZFKE8BiapgbsTkiAunVhdoMv1I87Pz/1Yy8sTO8i+hwJUoQyfjzUq6cClAJGo/WV6h5Lbr/wJyaqpUWQcuOGGgwKoGQJzbwxm2RShNDXokXw/fcQHAyfzzlNyNDsIWYnTYJ779W3cD5KghShBAWpCbCCg+G778wVrhZ8LYupaxdkG9U9xqoeUD2TAZvVPZJJEcL9/v4bRo1S67NnZFJvajf1pb3vvvynjRc3RYIUYdagAbyePWPnmDFw5IjN3SST4gQ2qnuMVT3BwRBULHfaRKp7hNBHZqaq5klJUeOzvXh9Dvz2m7lNnwx77zISpAhrzz8PDz2kBgHo29fqIulrmRRLnlDdYzWVj426HanuEUIfb70Fv/6qMpzLRh/Bb+pk8wNVquhYMt8nQYqwZjCoQd7CwmDXLjVzcg6+kknxiOoeG5mUUqWwOVGPZFKEcL+//zbX5rw2O5Oq43pBejo89pgazU24lAQpIrfKlc1tUiZPVmM/I5kUp8qnTUrJktiMSCSTIoR7Gat5rl+H6GgYeGk2/PGH+t6++65vnxQ9hAQpwrbevaFLF3VF7NsX0tJMD/liJsUTqnusMik2IhLJpAjhXm++aa7meW/UEQxTJpsfqFhR17IVFRKkCNsMBjUwUZkysH8/TJni0z8a3Hrh17R8q3usMinSJkUIXRw9alnNc4OqrzypJurp0EGqedxIghSRt3LlVEoTYNYsOHsW8M1Milsv/Ckp6mQHNqt7rNqkyASDQridsZonNVX1Ixh4cTbs2SPVPDqQIEXkr2tXVd2jaRg2/6x3aVzGrUGKsarH3x9KlDBtLqi6RzIpQrjHW2/Bjh0W1TxTp6gH3nwTKlTQt3BFjAQpomBvvAG33AJXkvQuicsYExtuYVnVY/GLzGbDWWmTIoRbnTqlBuAGmDsrkypjpZpHTxKkiIKFhcHSpaa72j/HdCyM81hW9+gSpOSYXNBmF2Tp3SOE22gaPPccJCerKc0GXpylqnkiIqSaRycSpAj7PPwwhjvvVOurV5urLHxEerob3+ziRbUsXdpqs1UmJZ/qHsmkCOEaK1fCDz+oWUKWvPwPfq9KNY/eJEgR9muhJtDSrl0zT2LhxXTLpBiDlDJlrDYb477wcMyZlKwsU0GlukcI17l8GYYMUevjXsmi9sy+6sTw2GPw5JP6Fq4IkyBF2M0QGGi+8957sGGDfoVxMrdmUi5cUEt7ghQwRSVS3SOE67z8MiQkwB13wOjI91TL2RIl4J13pJpHRxKkCIdpTe5WK/36mdtXeDmPDlKyoxLJpAjhGlu3qt9dAEumnyd4fHameMYMNQK30I0EKcJuph8TbdrAbbfBmTPwwgu6lulm6FbdY0+QYkybgClIkUyKEM6XmgrPPqvWn30W7vt0kGrF3qyZmnBV6EqCFOEwLTAIPvkE/Pzgiy/gyy/1LtJNk0yKEEXTjBlw5IhqFzvr/nXw9dfqy7ZkifWPBaELCVKE3ayqZZs1g7Fj1fpzz8G5c7qU6Wbonkmx6N2TlqYmMQNpkyKEuxw6pAbTBnhrdgrho7NTKqNGwV136VcwYSJBinCY6eI+cSLUr696qzzzjFePl693JsWYRTEYssdJ8bP4akomRQiny8pSp62MDOjYEbr+9rKqwq5ZEyZM0Lt4IpsEKcJuuRq4BwXBxx+r5Zo18NlnupSrsCxjKk8JUsLCLOKTHKPOSiZFCOf54ANzB563B+zDsOgd9cC770KxYvoWTphIkCIcZpUwuesumDRJrQ8ZYpqE0Nu4rbrnxg1zRGIjSLEahDZH6kQyKUI4x8WLMHq0Wp8yKYvKU55WJ7a+faF1a30LJ6xIkCLsludQAS+/DI0bq+7I/ft7ZbWP2zIply6ppcGghtrOZjNIyZE6kRFnhXCOV15RgUrduvBiyFLYu1elMefM0btoIgcJUsTNCwhQ1T4hIbB+vUqXegFdGs6eP6+W4eFWjWPzzaRIdY8QTrNrl3kqsnemXyZwwhh1Z/p0KFdOv4IJmyRIEQ6zmSipUwdmzlTrL78Mp0+7tUw3y22ZlIQEtYyKstpsc85Bqe4RwqkyM9XQJ5oGffrA/atHql8IDRvCoEF6F0/YIEGKsFuBI0O/+CI0b65myvOC3j66ZFLi49UyR5DiSHWPZFKEKJwlSyxqdh7/XbWeBTX0vYyJ4pEkSBEOyzP28PeH99+H4GBV7bNkiVvLdTPclkkxBinly1tttqe6RzIpQhReQoJqiwIwbUomUROeUXcGDIB77tGvYCJfEqQIu9k1x9Ydd5irfUaMgH/+cWmZboYuXZCNg97lEaRYtKXNFZVIJkWIwhszxlyz8xyLICZGfeGM5yvhkSRIEQ4rsBZn6FBo1QpSUtQkhF7w098bqnskkyJE4ezcaa7ZWfjqJfwnjlN3ZsyAsmX1K5gokAQpwm52z1bu56fOCCVKwPbtMH++S8vlDJ6SSZHePUI4V2ameR7Ufv2g+aqX1ASCjRvDwIH6Fk4USIIU4TC72sNWq2YOTsaNg4MHXVmkQtGluiePTEq+vXskkyJEoS1dCvv2qe/WrO77zCmVt9+WxrJeQIIUYTe7MylGAwZA+/YqAnjqKY9OAaSnu6kzUmEyKTnapEiQIoR9EhNh/Hi1PnVKFuUmP6/u9O4tjWW9hAQpwmF2X8wNBtXDJzwc9uzxuNEcc/4dqakufsOsrDzHSXGkTYoHx3pCeJRXX1Ujy9apA8+V+hx++w2KFzdPfSw8ngQpwm4OZ1IAKlaEN95Q65Mnq0EKPJTLg5SLF81pEIuRLTXNsTYpkkkRomBHj8Kbb6r1edNSCHjlZXVn/Hh1XhJeQYIU4Xq9e0OXLqoLTa9ekJamd4lsun7dxW9w5oxalisHgYGmzamp5jYx9lT3SCZFiIKNGqW+K48+Cu12v6omP731Vhg+XO+iCQcEFLxLbomJiezevZuEhASysrKsHuvTp49TCiY8l8NtNwwG1Xptxw44fFjlYKdNc0nZHOH26h5jkFKpktXmCxfUMjAQSpa0eEC6IAtRKBs3wpo16iv0+pBT0HGeemDePDXYpPAaDgcp3333HU8++STXrl2jVKlSGCzqAAwGgwQpPqxQ1T1GpUvDwoXwv/+p+uDOnaFJE2cVzSn0ClKMzVTKls1xjKULshAOy8xU40iC6npc+50hKlXZti106KBv4YTDHK7uGTlyJP379+fatWskJiZy+fJl0+2ScRp64dMK3Qvm//4PunVTZ5Enn4TkZKeWy1E5/w6XV/f8+69a3nKL1WbjxMi5xpSSCQaFcNhrr8GBA2ow2Un3/6xSKgEBsGDBTf7SEnpwOEg5c+YMQ4YMITQ01BXlER7MKd/vRYtUo7WjR9VsyR5E70xKrlniZYJBIRySmGge5f71OZlEThys7rz4opqyQ3gdh4OUdu3a8ccff7iiLMJL3NR4IpGR8NFHav2dd2DDBqeUqTDcnkkxBil5ZFJyBSkymJsQDhk9GpKSVJfjvulLVRu4MmVg4kS9iyYKyeEgpX379owaNYrJkyezatUq1qxZY3Vz1MKFC6lWrRohISE0a9aM3bt357t/YmIiL7zwAhUqVCA4OJjbb7+d77//3uH3FY5zWqY0Oto8TnX//ub+tzpzeSbFWN2TT5sUKzKYmxB2O3TIPPH6O6+l4DdlkrozcWKObnPCmzjccHZg9lwHU6dOzfWYwWAg04Ez6IoVKxgxYgSLFy+mWbNmLFiwgHbt2nHkyBHK5fpZCenp6Tz00EOUK1eOr776ikqVKnHq1CnC5QPoVk4ZmXX2bPjxRzVL8tCh5uyKjtyWSSlkdY8M5iZE3ozJkq5doeVvs9UXq2ZNePZZfQsmborDmZSsrKw8b44EKADz5s1j4MCB9OvXjzp16rB48WJCQ0NZtmyZzf2XLVvGpUuXWL16NS1atKBatWq0bNmS+vXrO/pniEJwapuz4sVVYOLnBx9/DF9/7cQXt49buyAnJ5szRo42nJXB3ITI17ffqlNIQABMHJSgWs+CaqASFKRv4cRN0W0wt/T0dPbs2UN0dLS5MH5+REdHs3PnTpvPWbNmDc2bN+eFF14gKiqKunXrMmPGjHyDo7S0NK5cuWJ1EzfHaXPc3HuvufHswIHmTINOXBqknDihlhERUKqU1UN5ZlLyCFI0TY2wL4RQ3wfjsEsjR0L9L8dDSgo0b656FAqvVqggZevWrXTo0IGaNWtSs2ZNOnbsyLZt2xx6jQsXLpCZmUlUjjlMoqKiiDfOFJvD8ePH+eqrr8jMzOT7779nwoQJvP7660zLZ2CwmTNnEhYWZrpVrlzZoXIKM5f03psyRU2ZfukS9O3r1quvWxvOGoOU6tVzPZRnJiVH6iTAonJWsilCKO+9B3/8AcWKwYhH/4L331cPvPaadDn2AQ4HKZ9++inR0dGEhoYyZMgQhgwZQrFixWjTpg2ff/65K8pokpWVRbly5ViyZAmNGzeme/fujBs3jsWLF+f5nLFjx5KUlGS6nT592qVlFA4KCoLPP4fQUNi0Cd56S7eiuHTYljyCFE1zPJNisUmIIi0lBcaMUevTpkG5OS+pHzpdu6pMrfB6DjecnT59OnPmzGG4xfwHQ4YMYd68ebz66qv07NnTrtcpU6YM/v7+nDNOXZ/t3LlzlM8xjb1RhQoVCAwMxN/ibH3HHXcQHx9Peno6QTbqHoODgwmWYZCdqqDqHk1z8AfM7berXz3PP6/OOG3bumVMg5x/x9WrLnwzY5BSrZrV5qQkdaIFyPWxz6MLMkgmRQhQSZNLl9TXauhdm2HkOvVFMQ6WIryew5mU48eP08HG0MIdO3bkhPFEbIegoCAaN27Mpk2bTNuysrLYtGkTzZs3t/mcFi1a8M8//1jNF3T06FEqVKhgM0ARzlVQ4JGZCRMmqBHwDQb1Y2bLFjvbsAwaBA8/rBqG9OplnnHPDYoXV0u3BCk5MilxcWpZurS5HCY5qnskkyKE2fnz5h49L7+Uhf+YUerOs8+qHz7CJzgcpFSuXNkqsDDauHGjw+09RowYwdKlS/noo484fPgwzz33HMnJyfTr1w9QkxWOHTvWtP9zzz3HpUuXGDp0KEePHmXdunXMmDGDF4xjbgi3yCvoWL5cpVwvX1b3v/kGHnxQtRX9/fcCXtRgUD+LIiNh7141CaGbGCf1u3bNhW+SR5BirH2sUsXGcySTIkSeXnlFdZhr0ACeCVsBe/aoL7MM3OZTHK7uGTlyJEOGDCEmJoZ7s+v8fv31Vz788EPeeOMNh16re/funD9/nokTJxIfH0+DBg1Yv369qTFtXFwcfn7mOKpy5cr8+OOPDB8+nLvuuotKlSoxdOhQRo8e7eifIQohv0zKihUqAWJLUhLcfbcaYHbQoHxep2JFWLxYze8zY4aaYz2PrJozGIOtkiUhPt6FmRRNKzCTYk+QYvFVkCBFFGm7d5vbx749Lx3/ftk/ZkePttG4S3gzh4OU5557jvLly/P666+zcuVKQLULWbFiBZ06dXK4AIMHD2bw4ME2H9uyZUuubc2bN+e3335z+H2E89jKpLz5pnm9Th146inV4j77IwKoJifFiqnH8vT442rywc8+g549Yd8+l48W6fJMijEC8vODGjWsHjIGKTaTkDlGnDUYVJVPZqZU94iibdw4dR7q0wda7H0LTp1SgyRatJUUvsHhIAWgS5cudOnSxdllER4urwyIpsGOHWp92DCYP9/82FtvwdNPw3ffqfv9+sH+/TBvXj5vtHAh7NwJx4+rJ3/5pUu6ElpmUsCFmZS//lLLGjUgRyPufDMpNmYUNAYpkkkRRdW2bbBxozolTBmeCA9mD0Hx6quql6DwKboN5ia8V85MyowZ5vWc1cHlyqmZ0o2pWVBBzJQp+TSoDQtTDVwCAmDVKlUF5EJuC1Jq1871kCNtUvLYJESRkZGhpvsClWit9uk01TClXj2VVhE+x64gJTIykgsXLgAQERFBZGRknjfhu2wlMzRNjWoP6mIfEWH7uf37q2DFaPJk2Lo1nzdr2lTN7wMqhRsbW5gi5ytnJsVl1T35BCmOtEkBGRpfFG1r16rpvsqUgXdePmkeV2nOHOvub8Jn2FXdM3/+fEpmn8nnz5+PQUbxK9IsMyAXLsDRo2rdOMlvXjp0gJ9/htat1f0HH4QfflA9j20aPhw2b1Znpm7dVOv9EiVuuvw5hYWpZVKS019aMQYpOcZ+ycw0H7N8q3ssIhIbNUBCFAlJSTBkiFrv2xdKzXpFDVUQHQ3t2ulbOOEydgUpffv2Na0/lW+rR+HLbMWmJ0+qZYUKuaaksenBB1WblLvuUvd791bZhGLF8njDDz5QfQyPHoUXXnDJbMmlS6vl5csqHnD6D7LDh9UyRybl7Fn1fgEBNgZyg3yreySTIoqaN95QQf2tt8Ir7WOh9RfqHDFnjgx/78McbpOyd+9eDhw4YLr/7bff0rlzZ1555RXS3TgAl9CPZSbl1Cm1rFrV/ufXq6dOOKAyMSVKQFpaHjuXKQNffGGeLdmJQYrx7zAGKZpmHuPFaa5dMzc8qVXL6iHj5kqV8giM8qnukUyK8FaaBn//7dg0XZcvw5Ilan3SJIickz0Wfs+e0LCh8wspPIbDQcqzzz7L0ez8/vHjx+nevTuhoaF8+eWXvGyc0Vb4JFs/Voy9dBwJUkClbY1j8GVlmSdDtun++1VLW1D9mI3VJ04SGGju5Zzd9Mp5jHVhZcuao6Fs+bZHAZvVPZJJEd7u7bfVgLDPPWf/c4YOVZOkV6gAXaJ2wPr16stgPC8In+VwkHL06FEaNGgAwJdffknLli35/PPP+fDDD1m1apWzyyc82I0balAlgPr1HX/+a6+BsSf7m2/CokX57Dx2LLRpoya66dbNKVMWW2aEypRRS6cHKXlU9YAdQYo0nBU+yNgD0JgZKciFC2oOUoBvvtYoMT174Lb+/VXdj/BpDgcpmqaZ5s7ZuHEjjz76KKBGg73g9DO88ETGi/v58+aLZWGSaCEhakgUo+efhytX8tjZ3x8+/VT1aT5wwKmDNhkM5iTHxYtOe1klJkYtjY1wLBQmSJEuyMLb2Wx/lgdNU3OOZmaqWp1mVzbAL7+o8YYmTHBdIYXHcDhIadKkCdOmTeOTTz5h69attG/fHoATJ06YhrMXvilndU98vFqWL1/4xqYVKqg5BY3dgG+9VQ17YFP58ipQMRjg3Xeth7O9SWXLqmWOSblv3r59ammj3tzYJiXPKa/y6d0jmRThrRwJUn76SY2x5OcH06dpaqhZUHVFt9zimgIKj+JwkLJgwQL27t3L4MGDGTduHDVr1gTgq6++Ms3lI3ybMZNiGaTcjOBg84BwFy6o7EqeA7099JCq+gEYOFCNSltIlu9hPN8ZAwen0DQ1WSLYDFIkkyKKInuDFE0zt5N/5hl4JP1bNddG8eLmc4DweQ4Pi3/XXXdZ9e4xmjt3Lv4ymI5Pyy+TcrMGD1YBypQpMH68qgoaOTKPnadMUSPB/fordO+ulkFBhX5vg8GczShorBeHxMWpbgkBAXDnnbkeNnbfzrPRsbRJET7I3pHr33pLdewD6Nk9E17Mrt4ZOlQmESxCHM6knD59mn8tzuS7d+9m2LBhfPzxxwQGBjq1cMIzOTuTYjRsmJoIGdTQ+Xm2jQ0IUGevyEj1y2rMmEK9n2UmxRikODWTsmePWt55Z645exITzd2dq1XL4/n59O6RTIrwVvZmUt5+Wy2nTYP7/1sBBw+qkRdfesl1hRMex+EgpWfPnmzevBmA+Ph4HnroIXbv3s24ceOYOnWq0wsoPEfOTIqx/YazgpTwcPjzT5VFOXMGqlfP52JcubIa6A1URGM55r6DDAYXVff8+qta3nNProdOnFDLcuXyGURXMinCB4WEFLzPvHlqLJWAAHhxUIYaHAVg1Ki8594QPsnhIOXgwYPcfffdAKxcuZK6deuyY8cOPvvsMz788ENnl094oJyZFGe2lw4Lg5kz1fq5c7BuXT47d+yo0i+gJhczjklip7wyKXm2h3GUMUi5775cDxmb0lSvns/zZTA34YMsMym2vmvHjpmregcOhFLffKQm7ClbVlX1iCLF4SAlIyOD4OzU9caNG+nYsSMAtWvX5uzZs84tnfAormyTYmnYMPMcHZ07m2tNbJo9WwUBSUnQqVOhJuCxzKRcv+6kUWdTUswFb9Ei18PGTIqjQYoM5ia8nWWblJxjpWgazJ2r1qtWhbdeSwNjhn7sWJfM3SU8m8NByp133snixYvZtm0bGzZs4OHs2eH+++8/SucYUVP4Jle1SbFkDFIAnn02nx2DguCrr1SU8ddfajIgR8bbzhYSYu6G7JQqn127VHBRoYLNRid2BSkywaDwQZaZlEGDrB97/301ugCoxvT+772rvpCVKjk2RK3wGQ4HKbNnz+bdd9+lVatW9OjRg/rZQ42uWbPGVA0kfJNlJkXTVLsRUNdhZ7v1VnMiYs8eeO+9fHaOioJvvlGNU7/7zjykZQFyppqd2i5lwwa1bNPG5nwCxuqeGjXyeQ3JpAgfZBmkWLYn1zT47DPz/ef7JsP06erOhAn2NWYRPsfhIKVVq1ZcuHCBCxcusGzZMtP2Z555hsWLFzu1cMJzJSSoufMMhnx6p9ykRo3Mw+YPHAg7d+azc5MmsHSpWp8+HVassPt9jDGEU3v4/PSTWrZta/NhYy/+22/P5zWk4azwQZZBiuXM6fffD1u2qPVDhyD0/bfUiaZGDTUEviiSHA5SAPz9/YnI0cK6WrVqlJO+60WCpql2bKAGIsvRu9apPvoIjAm6SZMKuDj37m3unvjUU7BjR76vnTOT4rQg5cIF8yBu0dG5Hj5zRt38/FQglifpgix8kOWQRsYfOP/9Z25nfvvtcEeFRJgzR22YPFnNAiqKJLsGc2vUqBGbNm0iIiKChg0bYrA1HW62vcaTs/A5lv92Y5uKfKsrnKBkSTUabXS0qkF5+WV4/fV8njBrFhw5oqp9OnZUZ75atfJ9j5yZlJse0G3TJhUB1atnsy5s1y61rFu3gHaAkkkRPu6PP2DjRjWQtNHKlWCYP0+1YK9TB3r21K+AQnd2BSmdOnUy9ejp3LmzK8sjvICmmTvRREa6/v1at4ZHHoEfflBD5r/yinlCwFz8/dVAbw8+CL//rs5+27fbHHs+r0zKqVM3WeDvv1fLPKp6jDNHF9iES4bFFz7I8nunadYBysiRUL/ieTX2EaiePTKSeZFmV5AyyTiQTo51UbRYZlJSUtSyeHH3vO+6dapqJCYGBgyAL7/MJwNcvLgKFO6/X/X4iY6GbdvyHNDF+HcZEy6HDqmTZz4Jw7ylpsLq1Wq9UyebuxgzKc2aFfBaMsGg8EH5jUPUvz9qWIFr19QXvmtXt5VLeKZCtUkxunbtGleuXLG6Cd+naZCcrNbdEaSAChiMPRC//RbeeKOAJ5Qpo+qHqlZVQ1e2a1fgACh33qmCgIsXVR15oXz/PVy5otIyNsZHycxUKW6wI0iRTIooQv75B+qE/6fSpaDGwy/ULwXhSxwOUk6cOEH79u0pXrw4YWFhREREEBERQXh4eK7GtMK3WJ4v3B2kgJoJdfx4tf7RR3aMDHvLLarCOyoKYmOhfXtzwcn9/JAQuOMOtf7bb4Us5Pvvq2WPHqplbA6HD6sfiSVKqOr2fEmbFOGDbH1vQ0PVsANMm6aykS1aQPYYXKJoc3gW5F69eqFpGsuWLSMqKirfRrTCN1lmUuyd0dRZhg5VI1IePKjaqaxfX8ATatZU3YFbtlR9mLt0UY1qLbokWX6EW7VSr/3zz/B//+dg4Y4fVw1nAJ5+2uYuxqqeJk3sqGqXwdyED8oZpNxxR3YN6YkT5gGRpk+XLIoAChGkxMbGsmfPHmoV0GNC+B692qRYKlNGjes0fjz8+KPqvGOjVsXaXXep4MHYRahnT1ixAk3L/fFv00bNvvrzz4Uo3OLF6gzcrh3cdpvNXYxBil3jHspgbsIHGYOUHj3g3nuhb1/Vi4+np0NGhmpJ27KlrmUUnsPh6p6mTZty2qlTxQpvo0ebFEvjxpl7Jb73np0TAt5zj/q5FhQEX38NTz5pGj7fMvhq2VLV0vz1l4PtUq5cMf8KfOGFPHezu9EsSHWP8GkVK6qh70uWRGVRPvpIPTBliq7lEp7F4SDlvffeY/bs2Xz00Ufs2bOH/fv3W92E77LVJsXd1T1GTz2llh9+qMZ6skt0tLlb0MqV5hasFiIiIHumB7Zvd6BA87LHdahdGx591OYuycmqKgnsDFJs1O1Iw1nh7Wz+qJgxQ32oH3oImjd3e5mE53I4SDl//jzHjh2jX79+NG3alAYNGtCwYUPTUvg+vTMpoOINY3Ayf765+qlAHTuqeX5CQtASEtS2xESrXYzVR8YRMAt05Ih5dMx8xnX4/XeVvKlUSd0KZKNuRzIpwtsZgxTTj56TJ9WvDVDDSgthweEgpX///jRs2JCdO3dy/PhxTpw4YbUUvssT2qRYlmXCBDWL8NWruad8z1f79mqSkCDVeNYw7VXzXDuYgxS7MikZGWo4/uvX1a/A//0vz12Ncw/de6+d5ZQuyMKHmc4nM2eqD3SbNnY0MBNFjcMNZ0+dOsWaNWuoWbOmK8ojvITemRRQbUf69FFV2CNHqrHbGje288nNmqHdmwJbgMTLqrHr00/D669z331q1rOYGLh0qYBRdSdNUimS8HD44IN8eyQYpxKyO5ttmTbJHl1OMinC21llUuLi1PcGJIsibHI4k9K6dWtiY2NdURbhJfTsgpzTqFHqop+Vpaq1HZJdeINxEsD33oO6dbklZi3162tkZcE77+Tx3KwsVd80c6a6/847+dbhaJp57BW7g5QAi98Q2Y18JZMivJ1Vm5SZM1U2snVr9StDiBwczqR06NCB4cOHc+DAAerVq0dgjrHJO3bs6LTCCc+i92ButhQvDm++CU2bqi7JxoHS7GE6WfbsCeMrqTG5jx+HDh14uc40nmQcs2ZpvPiigbAwiyfFxsLo0eYqoilTVH/KfPzxh5ocuXhxsLvplmWQcuMG+PtLJkV4PVMm5eoV8+CHkkUReXA4SBk0aBAAU6dOzfWYwWAgU86ePk/T9G+TYqlRIzVm2z//qLHafvrJsXGgDAZU3+P9+1XD1zfeoMef45lOZ/5MvpPlzebzbO2tasc//4SjR9UTixVTQ3j361fgexh7J7dvbzWOXP5yBinBwZJJET7D8NtOlUVp1QoeeEDv4ggP5XB1T1ZWVp43CVB8m/HCn5Wl2omCZwQpfn7w1ltqfeNGqzawjileXE1uduQIhuHD6RvxHQBfH6mjJgxavVoFKCEhauKz33+3K0A5cQKWLVPrgwc7UB7LXkLZ3y3JpAhvZ8qkxMaoFcmiiHzc1ASDRok5unAK32bZ3VfvNilGDz9sHkNtxQr7npPnIHBVq8K8eXTeORqAn/0f4r/Zn6gRZVeuhIQEWLVKzUhoh4ULzUNAOFTtnjOTggQpwvuZvndZmSqD0qqVnsURHs7hIGX27NmssLgKPP7440RGRlKpUiVpUOvjjJkUizn6KFZMn7LY0rmzWn77rWr/Ya+8qoZur2WgRQu4kelH1697cbXns/D449lDZNonIwM++UStv/ii/WUCrDMp2UGKVPcIb6ddvQaAAU2yKKJADgcpixcvpnLlygBs2LCBjRs3sn79eh555BFGjRrl9AIKz2PZs8fGRL+6adlSDfh66RJ061bw/vYMp//226oL8q5daoJWR33/vUq8REWpCREdYjCYD7BU9whfsWcPAIbKleHBB3UujPB0Dl9i4uPjTUHK2rVr6datG23btuXll1/m999/d3oBhecwZhyuqR9CHtEexVJgIHzxhbqub94Mhw/b97z8Gtk2aGAexmHRIjVwnCOMbVF697auvbFbjtSJZFKEV/vvP7SDhwAwPHC/zHQsCuRwkBIREWGaYHD9+vVEZ48xoWmaNJwtIjxljBRbGjRQPWjAXM2SF7smJgQ6dIDbb1cByuOPg73za65bB9+ptrf072/fc3LJEZVIJkV4tTlz0LLH/KFaNV2LIryDw0FK165d6dmzJw899BAXL17kkewc9r59+2QUWh+Xs02Kp2VSjHr3VstPPjGNgZavgn7MGQzmMdt+/FG18ysoo7J0qQpuNA2eeALuuKPgctiUY5JByaQIb6Rp8MU7lzmy6Gc01BfO4CdZFFEwh4OU+fPnM3jwYOrUqcOGDRsokT1y1tmzZ3n++eedXkDheTw5kwIqOAgLg3//heXLnfOaXbvC1q1QubIa7+3JJ2HfPtv7njkDw4erE/Mzz5jnTiuUHJMMSiZFeKOvvoKeL0TQOH0HRJUHpKZH2MfhWvLAwEBeeumlXNuHDx/ulAIJz2ccI8WTevZYCgmBYcPUQLCvv64GlLXF3uoeowcegHffhUcfVdU4330Hr72m2v5FRUFSEsybp8ZqSU5Ww98vWnSTjYulTYrwcllZ8Nbr6UAQyZRgZWoHQIIUYZ/CNOXj77//ZvPmzSQkJJCVI58+ceJEpxRMeB7jSSUjQy3tHjlVB88+q4KUffvg4kUoXTrvfR05WT7yiKr6eecd1TbFRrwOqEa877zjhN5POap7JJMivIlxfKBtu4JM2y4kBeXzDCGsORykLF26lOeee44yZcpQvnx5DBZneIPBIEFKEeLJQUqFCmqstUOHVE+f//0v9z6OZlKMxoxRU/c8+6xqe2LpjjtUdVPnzqoR703LUd2T464QHi0mBrZssf2YZFKEPRwOUqZNm8b06dMZPXq0K8ojPFjOk0qQh/8gio5WQcr69baDFKPCnCwNBlX1M2IElC2rBqI9e1bdDw8vdJFzy6N3j1T3CE/3669w3315Py5BirCHw8noy5cv8/jjj7uiLMLLeHImBcwj0C5fDq6YucFgUIPHlS4Nzz2n5iZ0aoACUt0jvNaUKfk/LkGKsIfDQcrjjz/OT4WewU14s5wnFU8PUlq2hLp1VSNW44Bslgpb3eNWeVT3SCZFeLpjx8zrG6sPpFNHb/jCCU/jcHVPzZo1mTBhAr/99hv16tUjMDDQ6vEhQ4Y4rXDCs3l6kGIwwMCBMHQo/PCD6hac134eSwZzE15G06BZM9VVHyCBspSd+TZ1Wxn4do15P4/+3gmP4XCQsmTJEkqUKMHWrVvZunWr1WMGg0GCFB/mbW1SQGVTAH77TV3YLefs86pMinRBFl5i+3YwzpByG0cpW6s0/O9/RPnDY4/B2rXqMQlShD0cDlJOnDjhinIIL+TpmRRQ1T0lS6oRYg8ehPr1c+/j0SdLaZMivMyPP5rXl9EfXnnF9MENCTE/5tHfO+ExPGgOW+HpvK1NCqhz4z33qPVff7V+zKsyKdImRXiB9eth+nS1vphnua/6f9Cjh+lxbzhnCM9SqMHc/v33X9asWUNcXBzp6elWj82bN88pBROez1tOOPfeCxs2wI4dYGvmBo/+RSdtUoQXsaztb83PMHasGtkwm2RShKMcDlI2bdpEx44dqVGjBn/99Rd169bl5MmTaJpGo0aNXFFG4aG8oU0KQIsWapkzk+IV8phgUIIU4Wn++w/+/lut/0YzbrslFfr0sdpHghThKIere8aOHctLL73EgQMHCAkJYdWqVZw+fZqWLVvK+Ck+zhure0D1NPDzg5Mn1YnUyBure2QwN+GJbtyAxo3VepOgWJqxG15+OddJwvKuBCnCHg4HKYcPH6ZPdnQcEBDA9evXKVGiBFOnTmX27NlOL6DwXN4SpJQqBfXqqfUdO3I/7tEnyzx690gmRXiSL76A+Hi13jP9QzXj5tNP59rPMpMihD0cDlKKFy9uaodSoUIFjlmM2HPhwgXnlUx4HG/NpAA0baqWsbHmbV6RScmjd49kUoQn+fRTtSzjd4kXWKhm3rQxTbpU9whHORyk3HPPPWzfvh2ARx99lJEjRzJ9+nT69+/PPcZuFA5auHAh1apVIyQkhGbNmrF79267nrd8+XIMBgOdjeOfC7fyljYpADVrqqWtHvQefbLMo7pHMinCUyQkwC+/qPWNWQ8SVLoUDBpkc1+p7hGOcjhImTdvHs2aNQNgypQptGnThhUrVlCtWjXef/99hwuwYsUKRowYwaRJk9i7dy/169enXbt2JCQk5Pu8kydP8tJLL3H//fc7/J6icLw5k1K9ulpaBilekUmRwdyEB0tLg1atIDUVKgXEU48DamjnEiVs7i+ZFOEoh4KUzMxM/v33X6pUqQKoqp/Fixezf/9+Vq1aRdWqVR0uwLx58xg4cCD9+vWjTp06LF68mNDQUJYtW5ZvOZ588kmmTJlCjRo1HH5P4RzelEkxfkyMQ3Vb8uiTpQzmJjzY5s1w+LBaX33jMfzCSsHgwXnun2MWFSEK5FCQ4u/vT9u2bbl8+bJT3jw9PZ09e/YQHR1tLpCfH9HR0ezcuTPP502dOpVy5coxYMCAAt8jLS2NK1euWN1E4eS8mAcUapQdfRiDlPh4uH5d37I4RAZzEx7sm2/UclDpL2nCHnjxRQgLy3N/y3OIR/84EB7D4eqeunXrctzWz9FCuHDhApmZmURFRVltj4qKIt7YVDyH7du38/7777N06VK73mPmzJmEhYWZbpUrV77pcgvFz4vGK46IUL18QHVFBu+s7pFMivAUW7bAkiVqvdPF96F4cRg2LN/nWJ4zJEgR9nD4MjNt2jReeukl1q5dy9mzZ92apbh69Sq9e/dm6dKllClTxq7njB07lqSkJNPt9OnTLi2jL8t5UrGcrM/TGQx5V/l49MlSBnMTHuq559SyTrHjanTZ55+H0qXzfY5kUoSj7E7YT506lZEjR/Loo48C0LFjRwwWnzJN0zAYDGQ6cPYsU6YM/v7+nDt3zmr7uXPnKF++fK79jx07xsmTJ+nQoYNpW1ZWlvpDAgI4cuQIt956q9VzgoODCfamFp5exJuCFFCNZ2NifCOTItU9Qk/Hj8Nff0GAfxa/Xm9EUIg/jBhR4PMkMBGOsjtImTJlCoMGDWLz5s1Oe/OgoCAaN27Mpk2bTN2Is7Ky2LRpE4NtNL6qXbs2Bw4csNo2fvx4rl69yhtvvCFVOS7mzZkUAOPHI2cyzaNPnHm0SZFMitDT3LlqeW+J/YQnJcEzQ8DGD8ucJJMiHGV3kKJl/+xs2bKlUwswYsQI+vbtS5MmTbj77rtZsGABycnJ9OvXD4A+ffpQqVIlZs6cSUhICHXr1rV6fnh4OECu7cL1vKlNCpiDlLg4fcvhkDx692RlqUyQnOiFu337LSxerNafS5qluuyMGmXXcyVIEY5yqH+GwQWfqu7du3P+/HkmTpxIfHw8DRo0YP369abGtHFxcfh529WwiPD2TIo3VvdY9qjKzPSuHlbCN7z7rlp2i9pK93MroN8zcMstdj1XghThKIdOcbfffnuBgcqlS5ccLsTgwYNtVu8AbNmyJd/nfvjhhw6/nygcqe7RQR4jzoKKWyRIEe507Rps2qTWJ517DoO/P4wZY/fzPfq7JjySQ6e4KVOmEJZPH3hRtHhbgssYpJw5o675XpVJycgArIMUaZci3EnToHdvSE+HmsX/447kw9Crr3k4ZztIJkU4yqEg5YknnqBcuXKuKovwcN6eSalQQZX5xg2w7FDm0SdL4xCd2UGKZeZEevgId9q6FVavhuCgLJYm91RZ9bFjHXoNCVKEo+z+LeyK9ijCu3lbkBIQABUrqvW4OC/JpOQIUiSTIvSycqVa9qq0mVZshe7doVYth15DBnMTjrI7SNG84owuXMnbMylgbt935ox5m0efLPMJUiSTItzl4kX44gu13u3EHLXyyisOv45Hf9eER7K7usc4aJoQRt7WJgXMA2I6afop18sRpBgM6rhnZUkmRbhPv36QmAiVQi/xYMrP0Lkz1Kvn8OtIdY9wlBdeZoRefCGTEhGhlpcve2d1D8iAbsK9/vkHvvsO/P011l6PJpAbMH58oV5LghThKAlSRKF5Y5CSPfafVSbFo0+WOXr3gAyNL9zLONNxq6jDNND2wSOPQOPGhXotCVKEoyRIEXaTTIoOjJkUi4hEMinCXU6fNg+B3z3+TbVSyCwKSGAiHCdBiig0b2yTYgxSEhPN2zz6xGmjukcyKcJdRo+G8+fh1lLn6Zn1CbRuDffeW+jXk0yKcJQXXmaEp5BMihtImxShkytX4Kuv1PqK6x0pTspNZVFAghThOAlShN18rbrHyKNPlpJJETr5+mv1sbs9IoHGGb+pDEqrVjf1mjJOinCUBCmi0LyxuqdkSbW8dk3fctgtn0yKBCnCVeLiYMQItd7v6ltqZcKEm44sJDARjvLCy4zQiy9kUkqUUMvkZO+v7pEgRbjK66+rbGO9smd58cY8aNIE2rW76deV6h7hKAlSRKF5c5BimUnx6JOljSDFxiYhnCYz0zy67Jyrz6u2KBMnOuWLIkGKcJQEKcJuvpRJuXbNezMpEqQIV/rmG9WjJyIkhTapa6FBA3jsMae8tgQpwlESpIhC88Y2KcWLq2Vamvki79EnSwlShBtduABPP63Wn89aaB5d1klfEo/+rgmP5IWXGaEXX8qkgGqX4vEkSBFu9NlnkJQEd5Q9z7j0iXDnndCli9NeXzIpwlESpIhC88YgJSjIfJG/elXfsthFghThRh9+qJYvJM+lGKkqi+LElKkEKcJREqQIu+U8qXhjdQ/kbjzr0SdLCVKEm3z7LcTEQJD/DZ5IeR9q1YLHH3fqe0iQIhzlpZcZ4Qm8MZMC1t2QPZ4EKcINkpOhXz+1PiToXUpzCV55xelfcm/9YSP0Ix8ZYTdfyaQYG88aefQvOglShBusWqXGRakWmcS06yOgRg3o2dPp7yOZFOEoL73MCL35+XnvSSY0VO8SOECCFOEGy5ap5YD0xQSTrrIoxlEDnUiCFOEoCVJEoXhrFgUgJMT6vkefLCVIES62Zg1s3Qr+fln0ubYQqlSB3r1d8l4SpAhHOT9UFj7L8qTire1RAIoV07sEDrCMSDQNDAYJUoTT3LgBgwap9RdDl1Hl2mkYu0h1g3MBCUyEo7z497DQkwQpbmKMSECNV45kUoTz/PADnD0LZUteZ+a1wVCpkrkFrQtIJkU4SoIUYTdfzaR49MnSMkjJjkokSBHO8u67atnH8CkhpMHLL0NwsMveT4IU4SgJUkSheHObFK/NpEiQIpxo3Tp18zNkMfDKa1C+PAwc6NL3lCBFOMqLLzXC3SSTogMJUoQLZGTA4MFq/YUSH1GLozBmjMsjeI/+rgmPJEGKKBRfClI8muWBliBFOMmPP8LJk6otyuyrz6ssyjPPuPx9LTOwErAIe0iQIuxmeVLx5uoer+qCbDDkikokSBE364MP1LInn6s5etyQRQGp7hGO8+JLjdCTZFLcyBiV3LhhdVeCFFEYP/wAX3+t1gdenQcVKrgliwISpAjHSZAi7OarbVI8nmRShJNkZsLQoWr9hRIfcSd/ui2LAhKYCMdJkCIKxZeCFI8/cUqQIpxk61b4+28ID01T46JUqODyHj2WJJMiHCVBiigUb26TIpkUURRpGsyYoda7G76kJNdg7Fi3fiEkSBGO8uJLjXA3X6nuyTlWlcefLCVIEU6wZg1s2gTBATd4OXkiVKzo1iwKSJAiHCdBiigUbw5SXDQtietIkCJuUlaWGkwWYGSxRdTghMqi5Ozq5mISmAhHSZAi7OYrXZAlSBFFzZYtcPQolApJY8zVV1QW5emn3V4OyaQIR3nxpUboyZszKZaDuIIXnCwlSBE3yThHT0//lea2KG7OooAM5iYcJ0GKsJuvtEmRTIooSn78EVauVOvPJM9TMx3rkEUByaQIxwXoXQDhnbw5SJFMiigqNA2GD1frLxZfRsPkGBj7ti5ZFPCC75rwOJJJEXaTNik6kSBFFNK2bXD4MBQPymBa8jCVRRkwQLfySCZFOMqLLzVCT5JJcSMJUkQhvfWWWvYMWEkprsIrr+iWRQEJUoTjJEgRdpM2KTqRIEUUwg8/wFdfgcGg8WLKLKhcWdcsCkiQIhwnbVJEoXhzkJIzk+LxJEgRDtI0GDVKrQ8Ofo96qQdh0nu5RzJ0MwlMhKMkkyIKxZfapHj8iVOCFOGgjRvh0CEoHpTO1NRRcNtt0Lev3sWSTIpwmGRShN2kukcnEqQIB82bp5YDeJ9wkmDqYgjQ/3Qv46QIR3nx72GhJ28OUqThrPBla9fC+vXgZ8hiaPpcuOsu6NZN72IBkkkRjtM/tBZeQ7og68T4C1iCFFGA9HQYNEitD/N/ixo3TsC0NR7zhZXARDjKMz65wutIJsWNjFFVejogQYrI29dfw5kzUCE0iWk3xsA998Bjj+ldLBPJpAhHSZAi7CZtUnRi7JGRlgZIkCJs0zR48021/kzqmxQjFaZP96hoQIIU4SgJUkSheHOQ4nVdkI2Db0mQIvLx6aewcyeE+KfzTNYiaNMGWrfWu1hWJEgRjpI2KcJuvtImxd9flT8rS933+JOlZFJEAa5fh5Ej1fq4rFepyFmY/o2+hbLB479rwuN48aVG6MmbMyngZdkUCVJEAT78EM6fh6qhCYzRZkLHjtCsmd7FykUyKcJREqQIu/lKmxSwbpfi8SfLPIIUTYPMTJ3KJDxGZia89ppaH5nyKgGGLHj1VX0LlQcJUoSjJEgRheLN1T3gpZmU1FTAuuySTRFvvAHHj0PpwCT6swyeeEKNjeKBZDA34Sgvv9QId5JMik7yaDgLEqQUdadPq4mNAWZmvERx/zSYMkXfQuXD479rwuNIkCIKxduDFK8qfx7VPSBBSlH31lvqY/FAyX08zXvQv7+ap8dDSXWPcJQEKaJQvOoib4NXlT9HkGJZdglSiq6LF2HJErU+6uoEDKGhHp1FAQlShOM8IkhZuHAh1apVIyQkhGbNmrF79+489126dCn3338/ERERREREEB0dne/+wnl8pQsyWF/oPf5kmSNIMRikh4+A0aMhKQnqBh3hUb6Hl16CChX0Lla+JEgRjtL9UrNixQpGjBjBpEmT2Lt3L/Xr16ddu3YkJCTY3H/Lli306NGDzZs3s3PnTipXrkzbtm05c+aMm0tetHlVJsIGryp/joazkGukfFHE7NgB77+v1hen98cvqhyMGqVvoewggYlwlO5Byrx58xg4cCD9+vWjTp06LF68mNDQUJYtW2Zz/88++4znn3+eBg0aULt2bd577z2ysrLYtGmTm0te9PhSw1lvzqTksUkUIcYexv2CP6MFO1Q1T4kS+hbKDpJJEY7SNUhJT09nz549REdHm7b5+fkRHR3Nzp077XqNlJQUMjIyiIyMtPl4WloaV65csbqJm+dLQYrHy9G7x3KTRXJFFBGHDsH69WAgi3FpE+GOO2DAAL2LZRcJUoSjdA1SLly4QGZmJlFRUVbbo6KiiI+Pt+s1Ro8eTcWKFa0CHUszZ84kLCzMdKtcufJNl7uo8tU2KR7PRtpEgpSiSdNg6FC13sVvDbdyHGbPhgDvmOFEAhPhKK++1MyaNYvly5fzzTffEGI8a+cwduxYkpKSTLfTp0+7uZS+yasu8jZ4e3WPjeSKKAI++gg2bYJi/mnMyRoJDzwAjz2md7HsJoO5CUfpGn6XKVMGf39/zp07Z7X93LlzlC9fPt/nvvbaa8yaNYuNGzdyVz6jKwYHBxNsPMmLm+JLbVK8KhNko+GsjU3Cx6Wnmwdum5w5UWVRXlvuVVd7qe4RjtL1VB0UFETjxo2tGr0aG8E2b948z+fNmTOHV199lfXr19OkSRN3FFXk4FUXeRt8JZMiQUrRsWgRnD0LFYMvMIz5avj7pk31LpZDJEgRjtK9InPEiBH07duXJk2acPfdd7NgwQKSk5Pp168fAH369KFSpUrMnDkTgNmzZzNx4kQ+//xzqlWrZmq7UqJECUp4Qet2b+ZLmRSvKr80nC3yEhJg0iS1PintFYKCDDBjhr6FKgQJTISjdA9Sunfvzvnz55k4cSLx8fE0aNCA9evXmxrTxsXF4Wfxs33RokWkp6fzv//9z+p1Jk2axOTJk91Z9CLNqy7yNnhlJiUrC27cgIAAaZNSxIwdqwZuaxTyJwNS34fBw6B6db2L5TCP/64Jj6N7kAIwePBgBg8ebPOxLVu2WN0/efKk6wskCuRLQYrHs2xTlZYGAQHSJqUIOXgQjMNGvZ06AP/IcBg3TtcyFZYEKcJRXt6yQLiTdEHWSc4gBanuKUpmz1bL/wtcQ3N+UyO55TEulKezPIdomn7lEN7Dyy81Qi9edZG3wauqewICzFFhdlQiQUrRsG8ffPGFWh+bMQXq1YNnntG3UDfB479rwuN4RHWP8A7ScFZHISGQkpIrkyJtUnxXejo89RRkZsLjhi9prO2FBZu8ZuA2WyRIEY6STIooFF+q7vGKE2eObsiSSfF9b74J+/dDmcBE3tZegK5doXVrvYt1UyzPG1LdI+zh5Zca4U6SSdFRjiBFGs76tpQUmDtXrc/OGEG54CvmDV5M2qQIR0mQIgrF6y7yOUgmRXiy115TY6NUCzhNbz6BkSOhRg29i3XTvOK7JjyKBCnCbpJJ0VGO1IkEKb7r8GGYNk2tz7jxMoGVotRAKT5AghThKO9tgSV05UttUrxCjpay0nDWd40YARkZ8JhhHU9oy2HBl+Ajo2lLkCIc5eWXGuFOvppJ8YoTp7RJKRLWr1e3QEMG87WhGNq1g//7P72L5TTSJkU4SjIpecjMzCQjI0PvYngUPz+oWlWth4Z63wUyMDAQ/+zoxOuCLGmT4vNu3ICXXlLrL2pvUjP4X3h7vZdE0fbxoT9FuIkEKTlomkZ8fDyJiYl6F8XjlCwJixer9dKl4cQJfctTGOHh4ZQvXx5/f/PZ0itOnBKk+LyJE+HQIYg0XGa8Ng3GjIGaNfUullN5xXdNeBQJUnIwBijlypUjNDQUg3yrTJKT1cBSAJUqQUSEvuVxhKZppKSkkJCQAIC/fwWdS+SgPBrOSpsU37BnD8yapdYXac8ScWtpGD1a30K5gFT3CEdJkGIhMzPTFKCULl1a7+J4HMvar6Ag84XSWxQrVgyAhIQEQkLKAV5U55MjdSJtUnyHpqlqHk2DHoYv6KZ9CW99D9mfV18ig7kJR0nDWQvGNiihoaE6l8TzeWuCyfi/LVHCHHF5xd9SvLhaJicDUt3jS9asgS1bINiQxkxtjGoo+8gjehfLJbziuyY8igQpNkgVj+8y/m+9bvoTYxdUCVJ8yrlz5vkCh2nzqRp+Bd5+W99CuZCcWoWjvO1ULTyEt59sLMvvFX+LZFJ8jqZB//5qZNl6hgNM1ibDawuhfHm9i+YyXvFdEx5FMimiULz9ZON1XZCNQcq1a4A0nPUFP/wA338PQYZ0lmvdCWndQkUtPkwazgpHSZDiY+Lj4+nZsyfly5cnKCiIihUr8tprrzn1PVauXMhdd1UjJCSEZs2asXv3bqe+vjtYNuDzioArRyZFGs56t+RkGD5crQ/R3qBOsZOwZImXfBiFcB+p7vExzz77LBkZGWzcuJGIiAjOnTvn1DFffvppBQsWjGD+/MU8+GAzFixYQLt27Thy5AjlypVz2vu4mtdmUqS6xye89BIcPQoVDf8xTpsOU6fCrbfqXSwhPI4EKQXRNDVvuh5CQx3+ZZWWlsapU6fYuXMn0dHRNGrUyGnFMRjg88/n0bnzQHr16kdYGCxevJh169axbNkyxowZ47T3cjWvC1LyaDh744Yau8br/p4ibO1a86CIH2l9CG9cE4YN07VMQngqCVIKkpKi3+Re166Zf0Hb4caNGzz88MM8+OCDhIWF8fbbb3P06FE+//xzSuT4G2bMmMGMGTPyfb0///yTKlWqmO6np6fz1197eOqpsabYyc/Pj+joaHbu3Gn/3+UBvL26x3KMmrQ0Fc8Kz3f9Ojz/vFofzjyiA7bCe797YXezmydtUoQ9it43w4cNHTqU1q1bU79+fQBee+01qlWrxqJFixg1apTVvoMGDaJbt275vl7FihWt7l+8eIHMzEwiI6OstkdFRfHXX3854S9wH6/LPOQTpCQnS5DiLWbPhtOnobLhNNO1cTBpEjRooHexdCFBirCHBCkFCQ019ajQ5b3tFBMTw6effsqCBQustoeFhXH27Nlc+0dGRhIZGVnoonlF9iEfXptJyf4s+vurQCU1VQUpZcvqWDZhlz/+gOnTNcDAa9pIit19l5qfRwiRJwlSCmIwOFTlopdVq1Zx++23ExgYaNqWnJzM0aNHGTJkSK79C1PdU7p0Gfz9/bl06ZzVfufOnaO8l43t4O2ZFOMmY5AiPNv169CnD9y4YaAbK3g8ZC18vK9IVvMI4Qj5hviIy5cvk5zjarVkyRIAunbtmmv/wlT3BAUFUbt2Y37/fRODBnUGICsri02bNjF48OCbKL37eV0mJUfDWeOmixclSPEGEybA4cNQnnje4XkMc+dArVp6F0tXUt0j7CFBio9o1qwZCxcuZP78+Tz22GOsX7+esWPH8s477xBhY7riwlb39Ow5gilT+tK6dRMeeOBuFixYQHJyMv369XPGn+E2vpJJAf1qI4V9tm2DefNUNc9SnqZ0dCNz61khRL4kSPERvXr1Ii4ujjfffJNJkyZRr149vv76ax599FGnvk/btt1JTDzP9OkTGTIkngYNGrB+/XqioqIKfrIH8fO2YQyNEcmNG5CeDkFBtpIrwsMkJECvXqBpBvrzPo9F7oQPYr3wAyiEPiRI8REGg4Fx48Yxbtw4l79Xt26DmTx5sFf3KLFsCuAV1T2W7aKSkyEoSDIpHu7GDTWhcVwc1ORv5jMcPvwMbrlF76IJ4TUknBdFktdV9wQGqhuYohLJpHi2qVNh+3YoZbjCWh6j1ND+0KGD3sXyGNImRdhDghRRKF6RfciH12VSIFe7FMmkeK6ff4Zp09RVeLH2LLUaFleDpAjuuw/Cw+H++/UuifAGUt0j7OY1F3M7eGXPzxIlIDHRFKRIJsUznTsHTz6p2qE8zVJ6lFgLK/aaZ4Us4n75RVWFWYyWIESevPFULTyAtwcsXlfdA7kyKcYg5epVncojcsnKgt69IT4e7uQgbzAUlrwPt92md9E8hsEgAYqwn1T3iCLJMkjxmoArR5BSqpS6e+WKTuURucyaBRs2QDFSWEk3Qoc+Az166F0sIbyWBCmiULzmwp4Hr/wlZ0ydZDdCCQtTd5OSdCqPsLJtG0yYoNqhLOQF6txXGubO1blUQng3qe4RRZJXZlLCw9Xy8mXAHKQkJupSGmHh+HHo3l0jK8tALz7hqaj1sHKvl0bDQngOCVJEoXjNhT0PXtkmpXRptbx0CZBMiqdITISHH4azZw3U5QDvBA7DsGoNVKigd9GE8HpS3SOKJK/sgmycxkCCFI+Rmal68vz9N1ThFD/RlpLvzYcWLfQumhA+QYIUHxMfH0/Pnj0pX748QUFBVKxYkddee81pr//BBzPp06cppUuXpFy5cnTu3JkjR47k2m/hwoVUq1aNkJAQmjVrxu7du51WBmfwyi7IOYIUY+2PBCn6GTUKvv8eQrjON3Shwpin1HTHQginkCDFxzz77LMkJiayceNGTpw4wdq1a2nUqJHTXn/v3q08/vgLbN/+Gxs2bCAjI4O2bdtazcC8YsUKRowYwaRJk9i7dy/169enXbt2JCQkOK0cN8srq3skk+JR3noL5s9X6x/Qj0ZdqsH06bqWSQhfI0GKj0lLS+PEiRPs3LmT9PR0GjVqROvWrZ3y2poGb721ng4dnuLOO++kfv36fPjhh8TFxbFnzx7TfvPmzWPgwIH069ePOnXqsHjxYkJDQ1m2bFmer33s2DEMBgNr166lTZs2hIaGUqtWLXbt2mXa5+TJkxgMBlatWsUDDzxAsWLFaNq0KXFxcWzbto177rmH0NBQ2rRpQ2IBrUl9obrHmEm5ckUNjiXc59tvYehQ1ZNnJmN4oulx+OQTmThQCCeTb1QBNE0NS6HHzdG5LW7cuMHDDz/M8uXLeeihh1i4cCEdO3bkmo1x02fMmEGJEiXyvcXFxeX5XsYLe1L2z/jI7Atoeno6e/bsITo62rSvn58f0dHR7Ny5M8/Xi42NxWAwMG/ePCZMmEBsbCxVqlRhzJgxVvsALFq0iBkzZrBjxw7OnTtHr169mDVrFm+//TabN28mNjaWDz74IN9j5dWZlIsXAdWO1nhNvHBBpzIVQbt3Q48eGppm4BneZfRt38C6ddaTQAohnMIba+bdKiXFPDyFu1275th5b+jQobRu3Zr69esD8Nprr1GtWjUWLVrEqFGjrPYdNGgQ3bp1y/f1KlasaHU/Z9CUlZXFsGHDaNGiBXXr1gXgwoULZGZmEhUVZbVvVFQUf/31V57vFRsbS3h4OCtWrKBs2bIAdOzYkXfffde0T0xMDJGRkaxYsYLS2T1dWrZsyfbt2zl06BCh2dMyN23alPj4+Hz/Nl/IpPj7Q5kykJCghmIvX17HshUR//wDHR7TuH7dwCN8z8Ly0zD8tA2yP7NCCOeSIMVHxMTE8Omnn7JgwQKr7WFhYZw9ezbX/pGRkabsR2EYDPDCCy9w8OBBtm/fXujXMYqNjaVTp06mAAXgxIkT1KxZ02qfLl26mAIUgLi4OLp3724KUIzbOnXqlO/7+ULDWYBy5VSQ4kHNfXzW4cPQ+kGNhPMGGrKXFaWeIeDH76FaNb2LJoTPkuqeAoSGqoyGHjeL626BVq1axe23306gxeBRycnJHD16lDvvvDPX/jdb3fPii4NZu3Ytmzdv5pZbbjFtL1OmDP7+/pw7d85q/3PnzlE+n5/6sbGxNG/e3GpbTEwMDRo0sLrfrFmzXM+75557TPdTU1M5cuSIKZuUF68czM0YnKWmwvXrABgTVjkOt3Cy48chuo1G/DkD9djP98W7UfKnVXDXXXoXTQif5o2/J93KYPCOqubLly9b9bABWLJkCQBdu3bNtX9hq3s0TWPu3Bf59ddv2LJlC9WrV7faJygoiMaNG7Np0yY6d+4MqGqhTZs2MXjwYJvvk5SUxMmTJ2nYsKHV9piYGIYMGQLAlStXcu1z4sQJkpKSrLYdOHAATdOoV69evn+bV2ZSSpRQBb9xQ2VTKlWiXDn1kGRSXOfff6FN6yz+O+vHnRzk5+IdKfPT55AjYBZCOJ83nqqFDc2aNWPhwoXMnz+fxx57jPXr1zN27FjeeecdIiIicu1f2Oqe2bNf4McfP+e7776lZMmSprYfYWFhFCtWDIARI0bQt29fmjRpwt13382CBQtITk6mX79+Nl9z//79BAQEWAUWp06d4vLly6ZMSmxsLP7+/qa2L2Buo1K1alWrbbfeeislCmhI5JVBisGgqnwSEkxBimRSXOvwYXikXRanTvtRk7/ZWKwjZdZ/Cvfeq3fRhCgSvPFULWzo1asXcXFxvPnmm0yaNIl69erx9ddf8+ijjzrtPTQNVq1aBMCDD7ayeuyDDz7gqaeeAqB79+6cP3+eiRMnEh8fT4MGDVi/fn2uxrRGsbGx1KpVi5CQENO2ffv2ER4eTrXs+n5b+8TGxubKvsTGxhZY1QNeWt0D5iAluzuPMZMiQYrzbd0KnTtlkZjkx20cZWOJLpT//iO47z69iyZEkWHQNEc7unq3K1euEBYWRlJSEqWMc91nS01N5cSJE1SvXt3qYiiUpCQ1/DdAkyb6lqWwjP/jEiWqU6WK+h8fPw45aq08V5s28PPPakyOXr1YtgwGDIBHHlEjnwrn+OQTGDBAIyPDQHN2sKbMAMps+AIs2kgJURTldw11BWk4K+zmS+GsZXWPV/1dVaqoZXajZqnucS5Ng0mT1Mj2GRkGHmclm6r0p8zO7yRAEUIHUt0jiiTLIMWrRms1tr85dQqQ6h5nSkuD/v3h88/V/bHMYFrDr/Fb+zPkaEQuhHAPCVKE3bwq41AAyzYpXhWk5MikGO/+95/qmSy1lIVz4QJ06ZTF9h1+BJDBYgYxoGsifLzVO7r3CeGjpLpHFEm+lEkJC1MB5LFjOpbLix0+DM2bZLB9hx9hJLKehxkwrgJ8+aUEKELoTIIUYTdfyqR4bZBimUnRNAwGuP12tenoUf2K5Y2ysmDBfI1G9W/wz6lAqnGCHWGP0mbtCJg2TSYLFMIDyLfQhiLW4cluXjm2SA7G/63l9SczU6fCFEblymqZnAyXLwPmIOXIEZ3K5IWOH4cHW6QzfISB1IwA2rGe3xoPps7+5dC+vd7FE0JkkyDFgnFI+ZSUFJ1L4plKllS9Sbx5qhLj/9Zy+gCLVc8XEmKeSfD4cUAyKY5ITYUZ0zXq3XGDX34LojjXWOT3Aj9M/I2onavNmSohhEfwgd/GzuPv7094eDgJ2WOMh4aGYvCqkb5czzj/X2qqvuVwlKZppKSkkJCQQHh4OP7+/kydqmpN7Bj7zbPUqQPx8XDgADRp4ngmJT4eYmNVu5YLF9TotcHBUKqUulWoALVrw623elkElzc1ECGMGpLGybPBQAAt2cIHtedQ/YsZ0r1YCA8lQUoOxknwEmQyFJ8UHh5u+h9PmKBzYQrrrrvUgG6xsQAYZwqIjYWMjDziir//hqVLYf16Fdzk4Qb+ZOFHIBkYAgJUoFK7tpqn5v77oWlTFdB4kX37YNgzyfzyR3EgmEr8y+zgifSYcBt+L3/rM4GYEL5IgpQcDAYDFSpUoFy5cmRkZOhdHOFEgYGB+Fv2PfZWjRqp5W+/ASqxEhmpEiJ79+aY927vXpgyBb77ztzy2WDgUPXH+DqgG9uTG/D3tQokpoaQnBFEepb5gh1wI4PSRy5S8ch/VP32FHXZTL2ARTSor3Hbw7dieLCVmsMme84mT5KVBRs2wNszklj3S0k0ilOMFF5mDqN6nKH4a9Nk7BMhvIHmAd5++22tatWqWnBwsHb33Xdru3btynf/lStXarVq1dKCg4O1unXrauvWrbP7vZKSkjRAS0pKutliC6GPkyfVhNT+/pp29aqmaZrWqZPaNHt29j7JyZo2cqSm+fmpB0DT2rfX/pi+XuvQLtW0qbC3CC5qbVmvTfCbpn1/50vapVEzNO3nnzUtJUW3w6Jpmnb5sqYtmJuu3VY+yaq8PfhMi3vkGU2LidG1fEJ4O3dfQ3UPUpYvX64FBQVpy5Yt0w4dOqQNHDhQCw8P186dO2dz/19//VXz9/fX5syZo/3555/a+PHjtcDAQO3AgQN2vZ8EKcIn1Kyprr7Ll2uapmmvv67uPvywpmkbN2pajRqmK3RWt+7a5o9OaY88Yr5o+/lpWocOmrZokab98oum/fmnpp06pWkXLqgLfUKCpv37r7qmr1unaW+8oWlPD8jS7r4rRQsOyLAZuNzBIW2A3zJtye1ztL3PLNLSftikaZcuufxQpKVp2rpv0rReLeO04gHXTeUpRaI2zLBAO/LoME3bv9/l5RCiKHD3NVT3CQabNWtG06ZNefvttwHIysqicuXKvPjii4wZMybX/t27dyc5OZm1a9eatt1zzz00aNCAxYsXF/h+7p4cSQiXGDcOZsyABx6ALVs4/JeBOnUgwHCDf7VKRJFAVqXKrOy5mtc3N+KPP9TT/PzgySfV02vVKtxbp6fDwYOw6zeN3zZeZeevGn8nhOXaL4g06hNLk5JHqV8zmdqNQrmjTUXKtroTQ/moQk8/feEC7Nlyle1rLrJ9G+yKK8/1LPNQu3U5wOCIz3nyhXBKPN9HNQQWQjiFu6+hugYp6enphIaG8tVXX9G5c2fT9r59+5KYmMi3336b6zlVqlRhxIgRDBs2zLRt0qRJrF69mtjshoT5kSBF+IR//1WNWtPTYcQICA6m+ZzO/JZ5N334iE4Pp/Pqf/2J2a/a4AQHq3lpRoyAmjWdX5zz5+G3nRo71l3mj20p/PFPOIkZJWzuG8lFavodp0JoEuXDUqlQ9gZlKwQQEhlKUJBqx5qWGUBKZhApqf4kXcrk7PkA/j0fzP7z5TmTVjbXa0YRz+PFf6DHY1dpPqQphnuayWBsQriAu6+hujacvXDhApmZmUQZp3LNFhUVxV9//WXzOfHx8Tb3j4+Pt7l/WloaaWlppvtJSUmAOtBCeK1SpVT3pAkTYN48AMbyK534jo/pwsfrAZIpWRKefhp691YxDYArPvrBwdCyFbRsFQCUQtOyOHHiCjG/XmPfz5c5fOAGR8+U4FRKaS4RyO6sWnANdTsDxDjybleowTHuDj1I89sv0rxdKW7rehd+tbqCwcBVgGvXnP0nCiEwXzvdld/w+d49M2fOZMqUKbm2VzaO3CmEz/gFsK52uXoV5s9XN19yHDieAstjUAHObF2LI0SRc/HiRcLCclfzOpuuQUqZMmXw9/fnXI555s+dO2cayyKn8uXLO7T/2LFjGTFihOl+YmIiVatWJS4uzi0H2JNduXKFypUrc/r06SJd9SXHwUyOhSLHQZHjYCbHQklKSqJKlSpERka65f10DVKCgoJo3LgxmzZtMrVJycrKYtOmTQwePNjmc5o3b86mTZus2qRs2LCB5s2b29w/ODiYYBuDT4WFhRXpD5qlUqVKybFAjoMlORaKHAdFjoOZHAvFz01tvnSv7hkxYgR9+/alSZMm3H333SxYsIDk5GT69esHQJ8+fahUqRIzZ84EYOjQobRs2ZLXX3+d9u3bs3z5cv744w+WLFmi558hhBBCCCfTPUjp3r0758+fZ+LEicTHx9OgQQPWr19vahwbFxdnFbHde++9fP7554wfP55XXnmF2267jdWrV1PXODa4EEIIIXyC7kEKwODBg/Os3tmyZUuubY8//jiPP/54od4rODiYSZMm2awCKmrkWChyHMzkWChyHBQ5DmZyLBR3HwfdB3MTQgghhLBFRjsSQgghhEeSIEUIIYQQHkmCFCGEEEJ4JAlShBBCCOGRvDJI+eWXX+jQoQMVK1bEYDCwevVqq8c1TWPixIlUqFCBYsWKER0dzd9//221T7Vq1TAYDFa3WbNmWe2zf/9+7r//fkJCQqhcuTJz5sxx9Z/msIKOxddff03btm0pXbo0BoOBmJiYXK+RmprKCy+8QOnSpSlRogT/93//l2tU37i4ONq3b09oaCjlypVj1KhR3Lhxw4V/mWOccRxatWqV6zMxaNAgq308/ThA/sciIyOD0aNHU69ePYoXL07FihXp06cP//33n9VrXLp0iSeffJJSpUoRHh7OgAEDuJZjPhxP/3444zj4wnmioO/G5MmTqV27NsWLFyciIoLo6Gh27dpltY8vfB7AOceiKHwmLA0aNAiDwcCCBQustrvrM+GVQUpycjL169dn4cKFNh+fM2cOb775JosXL2bXrl0UL16cdu3akZqaarXf1KlTOXv2rOn24osvmh67cuUKbdu2pWrVquzZs4e5c+cyefJkjxs0rqBjkZyczH333cfs2XlPbjJ8+HC+++47vvzyS7Zu3cp///1H165dTY9nZmbSvn170tPT2bFjBx999BEffvghEydOdPrfU1jOOA4AAwcOtPpMWH6pvOE4QP7HIiUlhb179zJhwgT27t3L119/zZEjR+jYsaPVfk8++SSHDh1iw4YNrF27ll9++YVnnnnG9Lg3fD+ccRzA+88TBX03br/9dt5++20OHDjA9u3bqVatGm3btuX8+fOmfXzh8wDOORbg+58Jo2+++YbffvuNihUr5nrMbZ8JzcsB2jfffGO6n5WVpZUvX16bO3euaVtiYqIWHBysffHFF6ZtVatW1ebPn5/n677zzjtaRESElpaWZto2evRorVatWk4tvzPlPBaWTpw4oQHavn37rLYnJiZqgYGB2pdffmnadvjwYQ3Qdu7cqWmapn3//fean5+fFh8fb9pn0aJFWqlSpayOj6cozHHQNE1r2bKlNnTo0Dxf19uOg6blfyyMdu/erQHaqVOnNE3TtD///FMDtN9//920zw8//KAZDAbtzJkzmqZ53/ejMMdB03zvPGHPcUhKStIAbePGjZqm+ebnQdMKdyw0reh8Jv7991+tUqVK2sGDB3P9ze78THhlJiU/J06cID4+nujoaNO2sLAwmjVrxs6dO632nTVrFqVLl6Zhw4bMnTvXKm2/c+dOHnjgAYKCgkzb2rVrx5EjR7h8+bLr/xA32bNnDxkZGVbHq3bt2lSpUsV0vHbu3Em9evVMowCDOhZXrlzh0KFDbi+zK3322WeUKVOGunXrMnbsWFJSUkyP+epxSEpKwmAwEB4eDqi/Mzw8nCZNmpj2iY6Oxs/Pz5T69sXvR87jYFSUzhPp6eksWbKEsLAw6tevDxTdz4OtY2Hk65+JrKwsevfuzahRo7jzzjtzPe7Oz4RHjDjrTPHx8QBWFxLjfeNjAEOGDKFRo0ZERkayY8cOxo4dy9mzZ5k3b57pdapXr57rNYyPRUREuPLPcJv4+HiCgoJynZgtj1d8fLzN42l8zFf07NmTqlWrUrFiRfbv38/o0aM5cuQIX3/9NeCbxyE1NZXRo0fTo0cP06Rp8fHxlCtXzmq/gIAAIiMjrT4TvvT9sHUcoOicJ9auXcsTTzxBSkoKFSpUYMOGDZQpUwYoep+H/I4FFI3PxOzZswkICGDIkCE2H3fnZ8LnghR7jRgxwrR+1113ERQUxLPPPsvMmTOL/LDHRZVlfWq9evWoUKECbdq04dixY9x66606lsw1MjIy6NatG5qmsWjRIr2Lo5v8jkNROU88+OCDxMTEcOHCBZYuXUq3bt3YtWtXrgtRUVDQsfD1z8SePXt444032Lt3LwaDQe/ieGfD2fyUL18eIFfvlHPnzpkes6VZs2bcuHGDkydPml7H1mtYvocvKF++POnp6SQmJlpttzxeReVY5NSsWTMA/vnnH8C3joPxwnzq1Ck2bNhglT0oX748CQkJVvvfuHGDS5cu+dxnIr/jYIuvnieKFy9OzZo1ueeee3j//fcJCAjg/fffB4rW5wHyPxa2+NpnYtu2bSQkJFClShUCAgIICAjg1KlTjBw5kmrVqgHu/Uz4XJBSvXp1ypcvz6ZNm0zbrly5wq5du2jevHmez4uJicHPz88ULTdv3pxffvmFjIwM0z4bNmygVq1aXpGus1fjxo0JDAy0Ol5HjhwhLi7OdLyaN2/OgQMHrD6UxhN6nTp13F5mdzF2U65QoQLgO8fBeGH++++/2bhxI6VLl7Z6vHnz5iQmJrJnzx7Ttp9//pmsrCxT4OYL34+CjoMtReU8kZWVRVpaGlB0Pg95sTwWtvjaZ6J3797s37+fmJgY061ixYqMGjWKH3/8EXDzZ8KhZrYe4urVq9q+ffu0ffv2aYA2b948bd++faZW+bNmzdLCw8O1b7/9Vtu/f7/WqVMnrXr16tr169c1TdO0HTt2aPPnz9diYmK0Y8eOaZ9++qlWtmxZrU+fPqb3SExM1KKiorTevXtrBw8e1JYvX66FhoZq7777ri5/c14KOhYXL17U9u3bp61bt04DtOXLl2v79u3Tzp49a3qNQYMGaVWqVNF+/vln7Y8//tCaN2+uNW/e3PT4jRs3tLp162pt27bVYmJitPXr12tly5bVxo4d6/a/Ny83exz++ecfberUqdoff/yhnThxQvv222+1GjVqaA888IDpPbzhOGha/sciPT1d69ixo3bLLbdoMTEx2tmzZ003y1b4Dz/8sNawYUNt165d2vbt27XbbrtN69Gjh+lxb/h+3Oxx8JXzRH7H4dq1a9rYsWO1nTt3aidPntT++OMPrV+/flpwcLB28OBB02v4wudB027+WBSFz4Qttno0uesz4ZVByubNmzUg161v376apqluyBMmTNCioqK04OBgrU2bNtqRI0dMz9+zZ4/WrFkzLSwsTAsJCdHuuOMObcaMGVpqaqrV+8TGxmr33XefFhwcrFWqVEmbNWuWO/9MuxR0LD744AObj0+aNMn0GtevX9eef/55LSIiQgsNDdW6dOliFcRomqadPHlSe+SRR7RixYppZcqU0UaOHKllZGS48S/N380eh7i4OO2BBx7QIiMjteDgYK1mzZraqFGjtKSkJKv38fTjoGn5HwtjF2xbt82bN5te4+LFi1qPHj20EiVKaKVKldL69eunXb161ep9PP37cbPHwVfOE/kdh+vXr2tdunTRKlasqAUFBWkVKlTQOnbsqO3evdvqNXzh86BpN38sisJnwhZbQYq7PhMGTdM0+/MuQgghhBDu4XNtUoQQQgjhGyRIEUIIIYRHkiBFCCGEEB5JghQhhBBCeCQJUoQQQgjhkSRIEUIIIYRHkiBFCCGEEB5JghQhhBBCeCQJUoQQQgjhkSRIEUIUWqtWrRg2bJjexTApbHkuXrxIuXLlTDPZusoTTzzB66+/7tL3EMKXSJAihPA6zg6Opk+fTqdOnUxT0bvK+PHjmT59OklJSS59HyF8hQQpQogiLSUlhffff58BAwa4/L3q1q3Lrbfeyqeffury9xLCF0iQIoSXWLt2LeHh4WRmZgIQExODwWBgzJgxpn2efvppevXqBcD69eu57777CA8Pp3Tp0jz22GMcO3YMgCVLllCxYkWysrKs3qNTp07079/fdD8rK4uZM2dSvXp1ihUrRv369fnqq6/yLGNB+7dq1YohQ4bw8ssvExkZSfny5Zk8ebLVa1y9epUnn3yS4sWLU6FCBebPn2+VOXnqqafYunUrb7zxBgaDAYPBYFVNk5WVle/r5/T9998THBzMPffcY9p27NgxDAYDa9eupU2bNoSGhlKrVi127dpl2ufkyZMYDAZWrVrFAw88QLFixWjatClxcXFs27aNe+65h9DQUNq0aUNiYqLpeR06dGD58uX5lkkIkc3heZOFELpITEzU/Pz8tN9//13TNE1bsGCBVqZMGa1Zs2amfWrWrKktXbpU0zRN++qrr7RVq1Zpf//9t7Zv3z6tQ4cOWr169bTMzEzt0qVLWlBQkLZx40bTcy9evJhr27Rp07TatWtr69ev144dO6Z98MEHWnBwsLZlyxZN0zStZcuW2tChQx3av1SpUtrkyZO1o0ePah999JFmMBi0n376yfQaTz/9tFa1alVt48aN2oEDB7QuXbpoJUuWNL1PYmKi1rx5c23gwIHa2bNntbNnz2o3btyw+/VzGjJkiPbwww9bbVu1apVmMBi0Bx98UNu8ebN29OhRLTo6WmvVqpVpn9WrV2uA1qZNG23btm3a3r17tcqVK2v333+/9uijj2q///679ttvv2mlS5fW5s2bZ3reDz/8oAUFBWmpqan5/LeFEJqmaRKkCOFFGjVqpM2dO1fTNE3r3LmzNn36dC0oKEi7evWq9u+//2qAdvToUZvPPX/+vAZoBw4c0DRN0zp16qT179/f9Pi7776rVaxYUcvMzNQ0TdNSU1O10NBQbceOHVavM2DAAK1Hjx6aplkHKfbuf99991k93rRpU2306NGapmnalStXtMDAQO3LL780PZ6YmKiFhoZaBUM5gyPL7fm9vi05j4OmadrEiRO1iIgILSEhwbTtzTff1O68807T/cmTJ2uRkZHahQsXTNt69eqlVatWTUtOTjZte/jhh7WXX37ZdD82NlYDtJMnT+ZZJiGEItU9QniRli1bsmXLFjRNY9u2bXTt2pU77riD7du3s3XrVipWrMhtt90GwN9//02PHj2oUaMGpUqVMjUKjYuLA+DJJ59k1apVpKWlAfDZZ5/xxBNP4OenTgv//PMPKSkpPPTQQ5QoUcJ0+/jjj03VRpbs3f+uu+6yel6FChVISEgA4Pjx42RkZHD33XebHg8LC6NWrVp2H6P8Xt+W69evExISYrUtNjaWTp06UbZsWdO2EydOULNmTat9unTpQunSpU3b4uLi6N69O6GhoVbbqlevbrpfrFgxQLWFEULkL0DvAggh7NeqVSuWLVtGbGwsgYGB1K5dm1atWrFlyxYuX75My5YtTft26NCBqlWrsnTpUlP7k7p165Kenm56XNM01q1bR9OmTdm2bRvz5883Pf/atWsArFu3jkqVKlmVIzg4OFfZ7N0/MDDQ6jGDwZCrbczNcPT1y5Qpw+XLl622xcbGMnbsWKttMTExPPDAA1b3c+4TGxvL8OHDTfdTU1M5cuQI9evXN227dOkSgFUAJISwTYIUIbzI/fffz9WrV5k/f74pIGnVqhWzZs3i8uXLjBw5ElDjfhw5coSlS5dy//33A7B9+3ar1woJCaFr16589tln/PPPP9SqVYtGjRqZHq9Tpw7BwcHExcVZBT95cXR/W2rUqEFgYCC///47VapUASApKYmjR49aBQhBQUGmBsQ3q2HDhla9bZKSkjh58iQNGza02i8mJoYhQ4YAcOXKlVz7nDhxgqSkJKttBw4cQNM06tWrZ9p28OBBbrnlFsqUKeOU8gvhyyRIEcKLREREcNddd/HZZ5/x9ttvA/DAAw/QrVs3MjIyTMFBREQEpUuXZsmSJVSoUIG4uDirXkBGTz75JI899hiHDh0y9QoyKlmyJC+99BLDhw8nKyuL++67j6SkJH799VdKlSpF3759b2p/W0qWLEnfvn0ZNWoUkZGRlCtXjkmTJuHn54fBYDDtV61aNXbt2sXJkycpUaIEkZGRpmoqR7Vr146xY8dy+fJlIiIi2L9/PwEBAVaBxalTp7h8+TINGjQAVMbE39+funXrmvaJiYkhMjKSqlWrWm279dZbKVGihGnbtm3baNu2baHKKkRRI21ShPAyLVu2JDMzk1atWgEQGRlJnTp1KF++vKnthp+fH8uXL2fPnj3UrVuX4cOHM3fu3Fyv1bp1ayIjIzly5Ag9e/bM9firr77KhAkTmDlzJnfccQcPP/ww69ats2pjcTP72zJv3jyaN2/OY489RnR0NC1atOCOO+6wajfy0ksv4e/vT506dShbtqypnU1h1KtXj0aNGrFy5UpABSC1atWyer99+/YRHh5uatdja5/Y2Nhc2ZfY2Firqp7U1FRWr17NwIEDC11eIYoSg6Zpmt6FEEKIvCQnJ1OpUiVef/11lw24tm7dOkaNGsXBgwcLnZGxx6JFi/jmm2/46aefXPYeQvgSqe4RQniUffv28ddff3H33XeTlJTE1KlTATXQnKu0b9+ev//+mzNnzlC5cmWXvU9gYCBvvfWWy15fCF8jmRQhhEfZt28fTz/9NEeOHCEoKIjGjRszb948qzYiQoiiQYIUIYQQQngkaTgrhBBCCI8kQYoQQgghPJIEKUIIIYTwSBKkCCGEEMIjSZAihBBCCI8kQYoQQgghPJIEKUIIIYTwSBKkCCGEEMIjSZAihBBCCI8kQYoQQgghPJIEKUIIIYTwSP8P76hRW0gqnX4AAAAASUVORK5CYII=\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# plot transmission, compare to paper results, look similar\n",
"fig, ax = plt.subplots(1, 1, figsize=(6, 4.5))\n",
"wavelengths_nm = td.C_0 / trans_g0.f / nm\n",
"plt.plot(wavelengths_nm, trans_g0.values, color=\"red\", label=\"$\\delta=0$\")\n",
"plt.plot(wavelengths_nm, trans_g20.values, color=\"blue\", label=\"$\\delta=20~nm$\")\n",
"plt.xlabel(\"wavelength ($nm$)\")\n",
"plt.ylabel(\"Transmission\")\n",
"plt.xlim([1050, 1400])\n",
"plt.ylim([0, 1])\n",
"plt.legend()\n",
"plt.show()\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Results Comparison\n",
"\n",
"Compare this plot to published results:\n",
"\n",
"\n",
"\n",
"(Citation: Opt. Lett. 43, 1842-1845 (2018). With permission from the Optical Society)"
]
},
{
"cell_type": "code",
"execution_count": null,
"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.10.9"
},
"widgets": {
"application/vnd.jupyter.widget-state+json": {
"state": {
"13181b6ba8ff496bb01cd4687b3709f2": {
"model_module": "@jupyter-widgets/output",
"model_module_version": "1.0.0",
"model_name": "OutputModel",
"state": {
"_dom_classes": [],
"_model_module": "@jupyter-widgets/output",
"_model_module_version": "1.0.0",
"_model_name": "OutputModel",
"_view_count": null,
"_view_module": "@jupyter-widgets/output",
"_view_module_version": "1.0.0",
"_view_name": "OutputView",
"layout": "IPY_MODEL_bd861d41cf8b4a2cb2b41312faa9f441",
"msg_id": "",
"outputs": [
{
"data": {
"text/html": "
normalization: status = success \u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501100%0:00:00\nSi-resonator-delta-0: status = success \u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501100%0:00:00\nSi-resonator-delta-20: status = success \u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501\u2501100%0:00:00\n