From 28a93c2e4f475bd23a0318be79b6a2f2d7be49a5 Mon Sep 17 00:00:00 2001 From: Ben Gyori Date: Tue, 26 Mar 2024 11:30:20 -0400 Subject: [PATCH] Update scenario 5 notebook --- .../scenario5/epi_scenario5.ipynb | 82 +------------------ 1 file changed, 2 insertions(+), 80 deletions(-) diff --git a/notebooks/evaluation_2024.03/scenario5/epi_scenario5.ipynb b/notebooks/evaluation_2024.03/scenario5/epi_scenario5.ipynb index 5a3cae04..da709e51 100644 --- a/notebooks/evaluation_2024.03/scenario5/epi_scenario5.ipynb +++ b/notebooks/evaluation_2024.03/scenario5/epi_scenario5.ipynb @@ -214,43 +214,6 @@ { "cell_type": "code", "execution_count": 9, - "id": "6148a6ab", - "metadata": { - "scrolled": false - }, - "outputs": [ - { - "data": { - "text/plain": [ - "Text(0, 0.5, '# Molecules / cell')" - ] - }, - "execution_count": 9, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAArUAAAGwCAYAAAC3hZ1vAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAByuklEQVR4nO3deXQUVcI28Keq9+z7BkmIskMEBgQCo8MIEhEdEQaXTwU0wxzZZBmVRQRREcVRUQH19UWRd1AQjziKDgioqOwgKJvgyBKWLED2pde63x+9JE0nkIRe0uH5ndOnq2/drrp1s/STm1tVkhBCgIiIiIgoiMmBbgARERER0dViqCUiIiKioMdQS0RERERBj6GWiIiIiIIeQy0RERERBT2GWiIiIiIKegy1RERERBT01IFuQHOgKArOnTuH8PBwSJIU6OYQERFRAwghUF5ejpSUFMgyx+mudQy1AM6dO4fU1NRAN4OIiIia4PTp02jdunWgm0EBxlALIDw8HID9hyIiIsJr21UUBcXFxYiOjuZfkD7EfvYf9rV/sJ/9g/3sH77s57KyMqSmpro+x+naxlALuKYcREREeD3UWq1WRERE8BemD7Gf/Yd97R/sZ/9gP/uHP/qZUwcJ4IliRERERNQCMNQSERERUdBjqCUiIiKioMdQS0RERERBj6GWiIiIiIIeQy0RERERBT2GWiIiIiIKegy1RERERBT0GGqJiIiIKOgx1BIRERFR0AtoqH3mmWcgSZLbo2PHjq71RqMREyZMQGxsLMLCwjBixAgUFBS4bSM3NxdDhw5FSEgIEhIS8MQTT8Bqtfr7UIiIiIgogNSBbkCXLl2wadMm12u1uqZJU6dOxZdffok1a9YgMjISEydOxPDhw7F161YAgM1mw9ChQ5GUlIRt27YhLy8Po0aNgkajwQsvvOD3YyEiIiKiwAh4qFWr1UhKSvIoLy0txbJly/Dhhx/illtuAQC8//776NSpE3bs2IG+ffvi66+/xuHDh7Fp0yYkJiaie/fueO655zB9+nQ888wz0Gq1/j4cN9bz52ErKIASGgrZYAhoW4iIiIhasoCH2t9++w0pKSnQ6/XIysrCggULkJaWhr1798JisWDQoEGuuh07dkRaWhq2b9+Ovn37Yvv27cjMzERiYqKrTnZ2NsaNG4dDhw6hR48ede7TZDLBZDK5XpeVlQEAFEWBoiheO7bcUaNhOXUKIStWILRXT69tl9wpigIhhFe/dlQ39rV/sJ/9g/3sH77sZ37tqLaAhto+ffpg+fLl6NChA/Ly8jBv3jzcdNNNOHjwIPLz86HVahEVFeX2nsTEROTn5wMA8vPz3QKtc71zXX0WLFiAefPmeZQXFxd7dT6u80etrKQYpqIir22X3CmKgvLycgghIMs899GX2Nf+wX72D/azf/iyn8vLy726PQpuAQ21Q4YMcS3fcMMN6NOnD9LT0/Hxxx/D4MN/18+cORPTpk1zvS4rK0Nqaiqio6MRERHhtf2UarWwAQgPCUFYTIzXtkvuFEWBJEmIjo7mB5OPsa/9g/3sH+xn//BlP9c+D4eoWX03REVFoX379vjvf/+LW2+9FWazGSUlJW6jtQUFBa45uElJSdi1a5fbNpxXR6hrnq6TTqeDTqfzKJdl2as/cJLjh01SFP7C9DFJkrz+9aO6sa/9g/3sH+xn//BVP/PrRrU1q++GiooK/P7770hOTkbPnj2h0WiwefNm1/qjR48iNzcXWVlZAICsrCwcOHAAhYWFrjobN25EREQEOnfu7Pf2X0pSqwAAwmYLcEuIiIiIWraAjtQ+/vjjuPPOO5Geno5z585h7ty5UKlUuP/++xEZGYmcnBxMmzYNMTExiIiIwKRJk5CVlYW+ffsCAAYPHozOnTvjoYcewsKFC5Gfn4/Zs2djwoQJdY7E+p3K3r3CylBLRERE5EsBDbVnzpzB/fffj4sXLyI+Ph5//OMfsWPHDsTHxwMAXnvtNciyjBEjRsBkMiE7OxtLly51vV+lUmHdunUYN24csrKyEBoaitGjR+PZZ58N1CG5kVT2kVrYeDMIIiIiIl8KaKhdtWrVZdfr9XosWbIES5YsqbdOeno6vvrqK283zSs4/YCIiIjIP5rVnNoWxzH9AJx+QERERORTDLU+5Jx+IDj9gIiIiMinGGp9idMPiIiIiPyCodaHJE4/ICIiIvILhlof4vQDIiIiIv9gqPUltfOSXhypJSIiIvIlhlofknjzBSIiIiK/YKj1JU4/ICIiIvILhlofct58gSeKEREREfkWQ60P1ZwoxlBLRERE5EsMtb7kvKQXpx8QERER+RRDrQ85px/wRDEiIiIi32Ko9SVOPyAiIiLyC4ZaH5I4/YCIiIjILxhqfahm+gFDLREREZEvMdT6kopzaomIiIj8gaHWh2qmHzDUEhEREfkSQ60vqXlHMSIiIiJ/YKj1IYnTD4iIiIj8gqHWhzj9gIiIiMg/GGp9idMPiIiIiPxCHegGtGSu6QcWi9e3LYQAFMV+YwerFZBlSLJsv+KCLEOSJK/vk4iIiKi5Yqj1IXVSMgDAcirXVSZsNtiKimAtKoat6CKsRUWwFRXDVlwEa3ExlIpKKFVVUCprPVdWQlRXQ9hsEFarPcheKSg7Aq4ky5A0GkgGPWS9AbJeD8ngfNZDNoRA1ushh4VBFREOOSICqvAIqCIjIIdHQBURDlVEBOSICMhhYQzLRERE1Cwx1PqQvlNHAID5+HEcv2sYbCUlsF644J85topiH8kFIEwmoKICV71XjQbq6Gio4mKhjomFOjYGqrqe4+OgjouDpOa3FxEREfkHU4cPqePjoUpLgy03F6ajR2tWSBJU0dFQxURDHR0DVUyMPRBGRUMOC4McGgo5JMT+7Fw26CFpNIBKDUmjtk9tUKnso7CyDKEIQLHZR3FtNgib4nitQFjMECYTlOpqCKMRSrURwlgNpdoIxVgNUV0NW0UFlLIy2MrKYSsrhVJWDltZGWzlZVBKyyDMZsBigbWwENbCQpiudPCyDHVsLNSJiVAnJUKTkGhfTkyAJjER6sQkaBITIIeG+vJLQERERNcIhlofi1qyGJqDB6EOD7eHuoREqONiXfNtg4ViNMJWXAzrxSL7tInazxcvwnrxIqxFF2G7WATrxYuA1Qrr+fOwnj8PHDxY73bl8HBoUlKgadWq5rlVCjQp9mdVVBSnPBAREdEVMdT6mBwdjci77oIsB/eFJmS9HnJyMjTJyVesKxQFtosXYSkohLWwAJb8fFgLCmEtKHC8LoC1oMA+X7i8HKajR91HsmvvNySkVsh1BN5WraBJTYU2LQ2q8HBvHyoREREFIYZa8jpJlqGOj4c6Ph5Al3rr2SoqYc3Pg+XcOfvj7FlYztqfzefOwnb+ApSqKph++y9Mv/23zm2ooqKgSUuDSEyE0vZ6aNPToU1Lswfe2FiO8hIREV0jGGopYFRhoVC1bQtd27Z1rldMJs/Ae+4cLGfOwHz6NGwXLsBWUgJbSQkAwLRxo9v7pZAQaB0jupq0VGhT06BNT4MmNQ2a5KSgmwJCRERE9WOopWZL1umgy8iALiOjzvVKZSXMp0/DdOoUSn/9FeoLF2E5cxqWU7mw5OdDVFXVO7VB0mqhTU+Dtk2bmkdGBrRt2kAVHc0RXiIioiDDUEtBSw4Nhb5jR2jbt4elZ0/ExMS45i4rZjMsZ87CcjoX5tzTMOfmwpKba38+cwbCbK53WoMcEQFtmzbQZbRxD73p6ZBDQvx9mERERNQADLXUIslaLXTXZUB3necor7DZYMnLg/nESZhPnoT5xAn788mTsOTlQSkrg/GXX2D85ReP96qTkhwhNx06x8iutk0baFq14nV5iYiIAoifwnTNkVQqaFu3hrZ1a+CmP7qtU4xGmHNzawLvyZrgayspgTU/H9b8fFTt2OG+UY3GPn+3TRtoM9rYA6/jwekMREREvsdQS1SLrNdD37499O3be6yzFhfDcuoUTCdPuofeU6cgjEaYjx+H+fhxz21GRtpHdtvUBF1thmM6g07nj8MiIiJq8RhqiRpIHR0NdXQ0DN27u5ULRYG1oADmEydgOnHCHnhPnID5xAn7dIbSUhh//gXGny+ZziBJ0LRq5Qq5rukMGRlQJyZydJeIiKgRGGqJrpIky9A4bkwR2q+f2zrFaIT5VK5j3u4JR/C1h16lvByWM2dgOXMGlT/84L7NkJA6R3d1bdrw1sJERER1YKgl8iFZr4e+Q3voO7hPZxBCwHbxYs3o7slTrtFd8+nT9suRHT4C0+EjHttUJyS4j+46HpqUFF57l4iIrlkMtUQBIEkS1HFxUMfFIeTGG93WCYsF5tNnao3snnDN4bVdvAhrYSGshYWo2rnTfZuua+9muK656wy+qqgoPx4dERGR/zHUEjUzkkZT7+XIbKWlMJ886TF313zq1GWvvauKjq5zdFfbujUkrdYfh0VERORTDLVEQUQVGQlDt24wdOvmVl5z7V3nqG7NCK81Px+24mJUFxej+qefLtmgCprWrTzn7mZkQBUXx5PViIgoaDDUErUA7tfevcltnVJV5T66e7JmhFepqrLfVvhULrBli9v75LAwj2kM6vR0iIgIfx4aERFRgzDUErVwckgI9J07Q9+5s1u5EALWwvOXXJnBHnwtZ89CqaiA8cABGA8c8NhmSXKy+zQGx5UZ1MnJkBy3KiYiIvInhlqia5QkSdAkJkCTmIDQvn3c1ilmMyy5uR5zd00nTkApLYU1Lw/WvDxUbtvmvk29Htr09Drn76rCwvx5eEREAABFUWA2mwPdDGoCjUYDVSOu6sNQS0QeZK0WurZtoWvb1q1cURRcOH4cISUlsJ7KdZu7a87NhTAaYTp6FKajRz22qYqPgy69jcfcXU3r1pDU/FVERN5nNptx4sQJKIoS6KZQE0VFRSEpKalB53jwk4SIGkWOikLIdddB7tXLrVxYrbCcPet5ZYaTJ2E9fx628xdQdf4Cqvbscd+gRgNtaqp97m5aGrTpadCkpUGblgZNcjIDLxE1iRACeXl5UKlUSE1NhcypUUFFCIGqqioUFhYCAJKTk6/4Hn5aEJFXSGq1fepBejowwH2draLCdVWGS6+9K4xGmI8fh/n4cc+NqtXQtmrlCrna9DRoUlOhTU+HpnVryLwcGRHVw2q1oqqqCikpKQgJCQl0c6gJDAYDAKCwsBAJCQlXnIrAUEtEPqcKC4MhsysMmV3dyoWiwFpQYA+6J0/Cknsa5txcmHNPwZJ7GsJshvnUKZhPnULlpRuVJKiTk6BNS7cH3rRUe/hNT4c2NRUyP8SIrmk2mw0AoOUfv0HN+QeJxWJhqCWi5kuSZWiSk6FJTkZov35u64SiwFpYCPOpmpBrD7y5sOTmQqmshPVcHqzn8lC1Y4fHtlXxcdCm1h7htT9r09Kgioz01yESUYDxetvBrTFfP4ZaImqWJFmGJikJmqQkhPbp7bZOCAFbUZEr4JpP5cJ8+rQ9/J7Kha2kBLbzF1B9/oLnDScAyOHh0KS2hrZVa2hSU6Fp3Qra1q2haZ0KTasUyDqdvw6TiIi8hKGWiIKOJElQx8ZCHRsL9Ojhsd5WVgZz7mlYck/BfMmUBmthIZTycpgOH4Hp8JE6t69OTISmdWtoW7eyB93WraFNbQ1N69ZQJyTwWrxERM0QQy0RtTiqiAgYunaBoWsXj3VKdTUsZ8/CfPo0LGfOwnLmDMxnzsBy5gwsp09DqaqCtaAA1oICVO/d6/F+SaOBplUraFq3to/2tm4NTauaZU5tIKKrMWDAAHTv3h2LFi0KdFOCDkMtEV1TZIOhzmvwAo5pDSUlsJw+7Qi7Z+3LZ8/AfPoMLHl5EBaL/VbDJ0/Wvf2ICPt0hlatoElJgSYlBeqUFNdrOTKSc/yIqF6ffvopNBpNoJsRlBhqiYgcJEmCOjoa6uhoGG64wWO9sFphyS+wj+qePVMz2nv6NMxnz8J24QKUsjKYDpfVO7VBDgmBppU96DpDrzPwqlNSoI6L8/VhElEzFhMTE+gmBC2GWiKiBpLUamhbt4K2dSsAfTzWK1VV9qkNZ87AcvYcLOfcH7YLF6BUVcH0239h+u2/de9Do4E6ORmIj4c5PR2a1jUjvpqUVtAkJfKGFERNIIRAtcUWkH0bNKoG/4eG0w+ajr8ZiYi8RA4Jga5dO+jatatzvWI0wnIuzxFyz3qEXmt+AYTFAktuLpCbi9I65vRClu0nsjmDbnIyNMlJUCcl2ZeTkjjFgagO1RYbOs/ZEJB9H342GyFaRi5fYw8TEfmJrNdDd10GdNdl1LleWCywFBTCfPYMSo79Bm1ZKax5ebCeO2cf+c3LgzCb7WV5eXWeyAYAksEATVIS1EmJ0CR5hl51cjJUYWG+PFQiIr9jqCUiaiYkjQba1q2gTkmG8frrERMT43a/eqEosF64YA+5zkdePiz5+bDm5cGSnw9bURFEdTXMJ+y3JK6PHBbmCLvOoJvkEYBlvd4fh03kFwaNCoefzQ7Yvsn3mk2offHFFzFz5kxMnjzZNY/EaDTiH//4B1atWgWTyYTs7GwsXboUiYmJrvfl5uZi3Lhx+PbbbxEWFobRo0djwYIFUHPOGRG1MJIsQ5OQAE1CAgzdu9dZRzGZYM3Pd4TdPPdlRwBWysqgVFRcdm4vAKiioqB2ju4mJkCTmAh1QqJ9+kNiAtSJiZDDwznVgYKCJEmcAtDCNYuv7u7du/HOO+/ghkvONp46dSq+/PJLrFmzBpGRkZg4cSKGDx+OrVu3ArDf13no0KFISkrCtm3bkJeXh1GjRkGj0eCFF14IxKEQEQWUrNNBm54ObXp6vXWUykpYHGHXWuAZeq15eVCqqux3ZispgelI3VdyAOxTHdQJ8dA4wm5N+E2AOsERfuPjIWm1vjhcIiKXgIfaiooKPPDAA3j33Xfx/PPPu8pLS0uxbNkyfPjhh7jlllsAAO+//z46deqEHTt2oG/fvvj6669x+PBhbNq0CYmJiejevTuee+45TJ8+Hc888wy09fwSNZlMMJlMrtdlZWUAAEVRoCiK145NURQIIby6TfLEfvYf9rV/+LyfDQZoMjKgyahnbq8QUMrLXaO81vx8WAsLYC0ohMXxbC0ogFJWBlFdDcupXFhO5V52l6rYWHvQTUyAOt7xnJgITYL9WZ2Q4PcT3Pj97B++7Gd+7ai2gIfaCRMmYOjQoRg0aJBbqN27dy8sFgsGDRrkKuvYsSPS0tKwfft29O3bF9u3b0dmZqbbdITs7GyMGzcOhw4dQo86bp8JAAsWLMC8efM8youLi2G1Wr12bIqioLy8HEIIt3lx5F3sZ/9hX/tHs+nnuDj7I7MrZABax8NJmExQLlyA7fx5KOcvQDl/HsqF81AuXLSXXbgA5cIFwGKB7eJF2C5evOyoL7RayHFxUMXHQ46NtT9iYiDHxda8jo2F5KXw22z6uYXzZT+Xl5d7dXsU3AIaaletWoWffvoJu3fv9liXn58PrVaLqKgot/LExETk5+e76tQOtM71znX1mTlzJqZNm+Z6XVZWhtTUVERHRyMiIqKph+NBURRIkoTo6Gj+wvQh9rP/sK/9I6j6OTn5squFELAVF9tvPVx43jXi6/ZceB624mLAbIZy7hyUc+cuv0+1Guq4OKjj46COi4cqPh7q+HhXmet1bCyky9yZKaj6OYj5sp9b4vkz3333XaCbELQC9t1w+vRpTJ48GRs3boTez2fY6nQ66HQ6j3JZlr3+AydJkk+2S+7Yz/7DvvaPltTPqrg4aOPigC7111HMZlgLCx3htxDW8+cdjwu1lh3h12q1T4m4zOAFAECSoIqOdoTdePdHQjzkmFgoWg1gMEDmJc58ylffzy3h54O8J2Chdu/evSgsLMQf/vAHV5nNZsP333+PxYsXY8OGDTCbzSgpKXEbrS0oKEBSUhIAICkpCbt27XLbbkFBgWsdEREFB1mrhbZ1a2hbt75sPWE2w3rxoj3kXrhgH/09f8njwgVYL1wArFbYiopgKyqC6dixerdZBPuNM1RxcVDHxkIdF2ufAxwbV7PsWKeKjYMcGsIrPhA1QwELtQMHDsSBAwfcyh5++GF07NgR06dPR2pqKjQaDTZv3owRI0YAAI4ePYrc3FxkZWUBALKysjB//nwUFhYiISEBALBx40ZERESgc+fO/j0gIiLyOUmrddxF7QrTHhQFtpISe8gtvCTwOpcLC2E5fx4wGqFUVUHJzbXfze1KbdDrHQE39oohWI6IYAAm8pOAhdrw8HB07drVrSw0NBSxsbGu8pycHEybNg0xMTGIiIjApEmTkJWVhb59+wIABg8ejM6dO+Ohhx7CwoULkZ+fj9mzZ2PChAl1Ti8gIqJrgyTLUMfEQB0TA3ToUGcdRVFw8eJFROl0UC5ehK2oCNYLF2G9eAG2CxftI8JuyxchqqogjEZYzp6F5ezZK7dDo3GFX1WcI/g6g3BMreeYaKiioyG1wDmiRP7SrH96XnvtNciyjBEjRrjdfMFJpVJh3bp1GDduHLKyshAaGorRo0fj2WefDWCriYgoWEiSBFVYGDQREUA9lzirTamqsgfcCxdgu3ixJgRfLPIIwUp5OYTF0rD5vw6qyEioYmKgiolxBN0YqGKioY6JrSmLiYEqOgbq6Che/5eolmYVai8940+v12PJkiVYsmRJve9JT0/HV1995eOWERER2efeakNCoE1NvWJdxWSyB1+PEHwRtosXapaLimArKQGEgK20FLbSUuAytzh2a094ONSOEOwZhGNqRoEd62WGYGrBmlWoJSIiailknQ5ySgo0KSlXrCtsNnugvXgR1qJi2IqLYC0qgq2o2D4twnHCm7282H4VCEWBUl4Oc3k5cOpUw9oUGuoIuM7RX0f4jYqGKioKquhoqKKjoI62v5YjIiDxCgMUJBhqiYiIAkxSqVxzgBtyRohQFHsIdoRdzyDsCMdFRbAWF8FWXAJYrVAqK+23ST59umENk+WasBsV5RZ4VVHRrhCsinKUR0dDDg9nEKaAYKglIiIKMpIsQx0dDXV0NHD99VesL4SAUlYG68Ui9/DrXC4pga24BLbiYsdyMZTKSkBRXMG5wWoH4egoqCKjYA0xIHzyZOiucMk2oqvBUEtERNTCSZJkPwktMhLAlU+IA+w3xLAVl7hCrq2kJvBai4vd1znC8OWCsBg33gdHRlSDoZaIiIg8yFot5MQEaBITGvyemiBc7Hq2FBWh4lwe1PFxPmzttatNmzaYMmUKpkyZEuimBBxDLREREXlFXUFYURSIoiLIoaEBbFnwMpvN0PKqFQ3CUEtEREQtnxCApSow+9aEAA28s9yAAQPQtWtXqNVq/Otf/0JmZia+/fZbHzewZWCoJSIiopbPUgW8cOXLq/nErHOAtuEj1R988AHGjRuHrVu3+rBRLQ9DLREREVEz0q5dOyxcuDDQzQg6DLVERETU8mlC7COmgdp3I/Ts2dNHDWnZGGqJiIio5ZOkRk0BCKRQnlTXJLzlBxEREREFPYZaIiIiIgp6DLVEREREFPQ4p5aIiIiomfjuu+8aVf/kyZM+aUcw4kgtEREREQU9hloiIiKiZmjlypUICwur89GlS5dAN6/Z4fQDIiIiomboL3/5C/r06VPnOo1G4+fWNH8MtURERETNUHh4OMLDwwPdjKDB6QdEREREFPQYaomIiIgo6DHUEhEREVHQY6glIiIioqDHUEtEREREQY+hloiIiIiCHkMtERERUTMxZswYDBs2LNDNCEoMtUREREQU9HjzBSIiImrxhBCotlYHZN8GtQGSJAVk39cShloiIiJq8aqt1ejzYd23nPW1nf9vJ0I0IQHZ97WE0w+IiIiIKOhxpJaIiIhaPIPagJ3/b2fA9k2+x1BLRERELZ4kSZwC0MJx+gERERERBT2GWiIiIiIKegy1RERERBT0GjyntkePHg2+xtpPP/3U5AYRERERXauWL18e6CYErQaHWt6yjYiIiIiaqwaH2rlz5/qyHURERERETcY5tUREREQU9Bo8UhsdHd3gObVFRUVNbhARERERUWM1ONQuWrTIh80gIiIiImq6Bofa0aNH+7IdRERERERN1uQ5tb///jtmz56N+++/H4WFhQCA//znPzh06JDXGkdERERE1BBNCrVbtmxBZmYmdu7ciU8//RQVFRUAgJ9//plXSSAiIiIiv2tSqJ0xYwaef/55bNy4EVqt1lV+yy23YMeOHV5rHBERERFRQzQp1B44cAB33323R3lCQgIuXLhw1Y0iIiIiImqMJoXaqKgo5OXleZTv27cPrVq1uupGEREREV2LxowZ47O7uA4YMABTpkzxybabgyaF2vvuuw/Tp09Hfn4+JEmCoijYunUrHn/8cYwaNcrbbSQiIiIiuqwmhdoXXngBHTt2RGpqKioqKtC5c2fcfPPN6NevH2bPnu3tNhIRERFdFSEElKqqgDyEEF45hldffRWZmZkIDQ1Famoqxo8f7zpZ32nr1q0YMGAAQkJCEB0djezsbBQXF2PMmDHYsmULXn/9dUiSBEmScPLkSQD2CwD07t0bOp0OycnJmDFjBqxWq2ubn3zyCTIzM2EwGBAbG4tBgwahsrLStf5///d/0alTJ+j1enTs2BFLly71yvE2VoOvU1ubVqvFu+++izlz5uDAgQOoqKhAjx490K5dO2+3j4iIiOiqiepqHP1Dz4Dsu8NPeyGFhFz1dmRZxhtvvIGMjAwcP34c48ePx5NPPukKkfv378fAgQPxyCOP4PXXX4darca3334Lm82G119/HceOHUPXrl3x7LPPAgDi4+Nx9uxZ3H777RgzZgxWrFiBX3/9FWPHjoVer8czzzyDvLw83H///Vi4cCHuvvtulJeX44cffnAF9ZUrV2LOnDlYvHgxevTogX379mHs2LEIDQ31+z0OmhRqnVJTU5GamuqtthARERFRPWrPh23Tpg2ef/55PProo65Qu3DhQvTq1cttpLRLly6uZa1Wi5CQECQlJbnKli5ditTUVCxevBiSJKFjx444d+4cpk+fjjlz5iAvLw9WqxXDhw9Heno6ACAzM9P1/rlz5+KVV17B8OHDAQAZGRk4fPgw3nnnneAItSNGjEDv3r0xffp0t/KFCxdi9+7dWLNmjVcaR0REROQNksGADj/tDdi+vWHTpk1YsGABfv31V5SVlcFqtcJoNKKqqgohISHYv38/Ro4c2ahtHjlyBFlZWZAkyVXWv39/VFRU4MyZM+jWrRsGDhyIzMxMZGdnY/DgwfjrX/+K6OhoVFZW4vfff0dOTg7Gjh3rer/VakVkZKRXjrkxmhRqv//+ezzzzDMe5UOGDMErr7xytW0iIiIi8ipJkrwyBSBQTp48iTvuuAPjxo3D/PnzERMTgx9//BE5OTkwm80ICQmBwUvhuTaVSoWNGzdi27Zt+Prrr/Hmm2/iqaeews6dOxHi6M93330Xffr08XifvzXpRLGKigq3my44aTQalJWVXXWjiIiIiKjG3r17oSgKXnnlFfTt2xft27fHuXPn3OrccMMN2Lx5c73b0Gq1sNlsbmWdOnXC9u3b3U5m27p1K8LDw9G6dWsA9j8I+vfvj3nz5mHfvn3QarVYu3YtEhMTkZKSguPHj6Nt27Zuj4yMDC8efcM0aaQ2MzMTq1evxpw5c9zKV61ahc6dO3ulYURERETXotLSUuzfv9+tLC4uDhaLBW+++SbuvPNObN26FW+//bZbnZkzZyIzMxPjx4/Ho48+Cq1Wi2+//RYjR45EXFwc2rRpg507d+LkyZMICwtDTEwMxo8fj0WLFmHSpEmYOHEijh49irlz52LatGmQZRk7d+7E5s2bMXjwYCQkJGDnzp04f/48OnXqBACYN28eHnvsMURGRuK2226DyWTCnj17UFxcjGnTpvmry+xEE3z++edCrVaLUaNGieXLl4vly5eLhx56SKjVarF27dqmbDKgSktLBQBRWlrq1e3abDZx/vx5YbPZvLpdcsd+9h/2tX+wn/2D/ewfvuzny31+V1dXi8OHD4vq6mqv79eXRo8eLQB4PHJycsSrr74qkpOThcFgENnZ2WLFihUCgCguLna9/7vvvhP9+vUTOp1OREVFiezsbNf6o0ePir59+wqDwSAAiBMnTrjec+ONNwqtViuSkpLE9OnThcViEUIIcfjwYZGdnS3i4+OFTqcT7du3F2+++aZbm1euXCm6d+8utFqtiI6OFjfffLP49NNPvdIfjfk6NinUCiHEunXrRL9+/URISIiIjY0Vf/7zn8V3333XqG0sXbpUZGZmivDwcBEeHi769u0rvvrqK9f66upqMX78eBETEyNCQ0PF8OHDRX5+vts2Tp06JW6//XZhMBhEfHy8ePzxx11fiIZiqA1u7Gf/YV/7B/vZP9jP/sFQS1ejMV/HJl/Sa+jQoRg6dOhVjRK3bt0aL774Itq1awchBD744APcdddd2LdvH7p06YKpU6fiyy+/xJo1axAZGYmJEydi+PDh2Lp1KwDAZrNh6NChSEpKwrZt25CXl4dRo0ZBo9HghRdeuKq2EREREVHwuKrr1F6tO++80+31/Pnz8dZbb2HHjh1o3bo1li1bhg8//BC33HILAOD9999Hp06dsGPHDvTt2xdff/01Dh8+jE2bNiExMRHdu3fHc889h+nTp+OZZ56p82Q2ADCZTDCZTK7XzpPbFEWBoiheOz5FUex3MPHiNskT+9l/2Nf+wX72D/azf/iyn/m1o9oCGmprs9lsWLNmDSorK5GVlYW9e/fCYrFg0KBBrjodO3ZEWloatm/fjr59+2L79u3IzMxEYmKiq052djbGjRuHQ4cOoUePHnXua8GCBZg3b55HeXFxsdtt4a6WoigoLy+HEAKy3KQLTVADsJ/9h33tH+xn/2A/+4cv+7m8vNyr26PgFvBQe+DAAWRlZcFoNCIsLAxr165F586dsX//fmi1WkRFRbnVT0xMRH5+PgAgPz/fLdA61zvX1WfmzJluZ+SVlZUhNTUV0dHRiIiI8NKR2X+QJUlCdHQ0f2H6EPvZf9jX/sF+9g/2s3/4sp/V6oDHGGpGAv7d0KFDB+zfvx+lpaX45JNPMHr0aGzZssWn+9TpdNDpdB7lsix7/QdOkiSfbJfcsZ/9h33tH+xn/2A/+4ev+plfN6qtUd8NN910E/75z3/i2LFjXmuAVqtF27Zt0bNnTyxYsADdunXD66+/jqSkJJjNZpSUlLjVLygocN2zOCkpCQUFBR7rneuIiIiI6NrQqFA7duxYbN++HT179kSnTp0wffp0bN261e0uFFdLURSYTCb07NkTGo3G7c4YR48eRW5uLrKysgAAWVlZOHDgAAoLC111Nm7ciIiICN4EgoiIiOga0qjpB6NGjcKoUaNgMpmwefNm/Pvf/8bIkSNdl9b6y1/+guzs7Abfe3jmzJkYMmQI0tLSUF5ejg8//BDfffcdNmzYgMjISOTk5GDatGmIiYlBREQEJk2ahKysLPTt2xcAMHjwYHTu3BkPPfQQFi5ciPz8fMyePRsTJkyoc3oBEREREbVMTZqMotPpcPvtt+Odd97BuXPn8PnnnyM5ORlPP/00YmNjcccdd7iuJXs5hYWFGDVqFDp06ICBAwdi9+7d2LBhA2699VYAwGuvvYY77rgDI0aMwM0334ykpCR8+umnrverVCqsW7cOKpUKWVlZePDBBzFq1Cg8++yzTTksIiIiIgpSkvDm3AEAv//+Oz7//HOkpqbir3/9qzc37TNlZWWIjIxEaWmp169+UFRUhJiYGE5m9yH2s/+wr/2D/ewf7Gf/8GU/X+7z22g04sSJE8jIyIBer/fqfsl/GvN19PrVD66//npMnTrV25slIiIiavHGjBmDkpISfPbZZ16p11QDBgxA9+7dsWjRIp9s3xf4pykRERERBT2GWiIiImrxhBCwmGwBeXh5pmeDvPrqq8jMzERoaChSU1Mxfvx4VFRUuNXZunUrBgwYgJCQEERHRyM7OxvFxcUYM2YMtmzZgtdffx2SJEGSJJw8eRIAsGXLFvTu3Rs6nQ7JycmYMWOG291YP/nkE2RmZsJgMCA2NhaDBg1CZWWla/3//u//olOnTtDr9ejYsSOWLl3qtWMO+M0XiIiIiHzNalbwP5N9e3On+vz99T9Bo1P5dZ+yLOONN95ARkYGjh8/jvHjx+PJJ590hcj9+/dj4MCBeOSRR/D6669DrVbj22+/hc1mw+uvv45jx46ha9eurpPv4+PjcfbsWdx+++0YM2YMVqxYgV9//RVjx46FXq/HM888g7y8PNx///1YuHAh7r77bpSXl+OHH35whfqVK1dizpw5WLx4MXr06IF9+/Zh7NixCA0NxejRo6/6mL0Sam02Gw4cOID09HRER0d7Y5NERERE1ERTpkxxLbdp0wbPP/88Hn30UVeoXbhwIXr16uU2UtqlSxfXslarRUhIiNvNrJYuXYrU1FQsXrwYkiShY8eOOHfuHKZPn445c+YgLy8PVqsVw4cPR3p6OgAgMzPT9f65c+filVdewfDhwwEAGRkZOHz4MN55553AhdopU6YgMzMTOTk5sNls+NOf/oRt27YhJCQE69atw4ABA666YURERETeotbK+PvrfwrYvv1t06ZNWLBgAX799VeUlZXBarXCaDSiqqoKISEh2L9/P0aOHNmobR45cgRZWVmQJMlV1r9/f1RUVODMmTPo1q0bBg4ciMzMTGRnZ2Pw4MH461//iujoaFRWVuL3339HTk4Oxo4d63q/1WpFZGSkV465SaH2k08+wYMPPggA+OKLL3DixAn8+uuv+L//+z889dRTDbpGLREREZG/SJLk9ykAgXLy5EnccccdGDduHObPn4+YmBj8+OOPyMnJgdlsRkhISINvlNUYKpUKGzduxLZt2/D111/jzTffxFNPPYWdO3ciJCQEAPDuu++iT58+Hu/zhib96XDhwgXXcPRXX32FkSNHon379njkkUdw4MABrzSMiIiIiBpv7969UBQFr7zyCvr27Yv27dvj3LlzbnVuuOEGbN68ud5taLVa2Gw2t7JOnTph+/btbie+bd26FeHh4WjdujUA+x8P/fv3x7x587Bv3z5otVqsXbsWiYmJSElJwfHjx9G2bVu3R0ZGhleOu0kjtYmJiTh8+DCSk5Oxfv16vPXWWwCAqqoqr6VtIiIiIqpfaWkp9u/f71YWGxuLtm3bwmKx4M0338Sdd96JrVu34u2333arN3PmTGRmZmL8+PF49NFHodVq8e2332LkyJGIi4tDmzZtsHPnTpw8eRJhYWGIiYnB+PHjsWjRIkyaNAkTJ07E0aNHMXfuXEybNg2yLGPnzp3YvHkzBg8ejISEBOzcuRPnz59Hp06dAADz5s3DY489hsjISNx2220wmUzYs2cPiouLMW3atKvujyaN1D788MO455570LVrV0iShEGDBgEAdu7ciY4dO151o4iIiIjo8r777jv06NHD7TFv3jx069YNr776Kl566SV07doVK1euxIIFC9ze2759e3z99df4+eef0bt3b2RlZeHf//431Gr7eOfjjz8OlUqFzp07Iz4+Hrm5uWjVqhW++uor7Nq1C926dcOjjz6KnJwczJ49GwAQERGB77//Hrfffjvat2+P2bNn45VXXsGQIUMAAH/729/wv//7v3j//feRmZmJP/3pT1i+fLnXRmqbfJvcTz75BKdPn8bIkSNdQ84ffPABoqKicNddd3mlcf7C2+QGN/az/7Cv/YP97B/sZ//gbXLpavjlNrl//etfXTtz8sblGIiIiIiIGqtJfzLZbDY899xzaNWqFcLCwnD8+HEAwNNPP41ly5Z5tYFERERE15Lc3FyEhYXV+8jNzQ10E5ulJo3Uzp8/Hx988AEWLlzodq2xrl27YtGiRcjJyfFaA4mIiIiuJSkpKR4ngF26njw1KdSuWLEC//M//4OBAwfi0UcfdZV369YNv/76q9caR0RERHStUavVaNu2baCbEXSaNP3g7NmzdXa2oiiwWCxX3SgiIiIiosZoUqjt3LkzfvjhB4/yTz75BD169LjqRhERERERNUaTph/MmTMHo0ePxtmzZ6EoCj799FMcPXoUK1aswLp167zdRiIiIiKiy2rSSO1dd92FL774Aps2bUJoaCjmzJmDI0eO4IsvvsCtt97q7TYSEREREV1Wk69Te9NNN2Hjxo3ebAsRERERUZPwFipEREREFPQaHGqjo6MRExPToAcRERERNd3p06fxyCOPICUlBVqtFunp6Zg8eTIuXrwY6KY1Ww2efrBo0SIfNoOIiIiIAOD48ePIyspC+/bt8dFHHyEjIwOHDh3CE088gf/85z/YsWMHBxHr0OBQO3r0aF+2g4iIiMhnhBCwmkwB2bdap4MkSQ2uP2HCBGi1Wnz99dcwGAwAgLS0NPTo0QPXX389nnrqKbz11lto06YNcnJycPjwYXz++eeIiorCrFmzMGHCBF8dSrPWpBPFrnTP4bS0tCY1hoiIiMgXrCYT3hj914Ds+7EPPoFGr29Q3aKiImzYsAHz5893BVqnpKQkPPDAA1i9ejWWLl0KAHj55Zcxa9YszJs3Dxs2bMDkyZPRvn37a/JqVE0KtW3atLnsXxw2m63JDSIiIiK6Vv32228QQqBTp051ru/UqROKi4tx/vx5AED//v0xY8YMAED79u2xdetWvPbaawy1DbVv3z631xaLBfv27cOrr76K+fPne6VhRERERN6i1unw2AefBGzfjSWEaFC9rKwsj9fX6nlQTQq13bp18yjr1asXUlJS8PLLL2P48OFX3TAiIiIib5EkqcFTAAKpbdu2kCQJR44cwd133+2x/siRI4iOjkZ8fHwAWte8efU6tR06dMDu3bu9uUkiIiKia0ZsbCxuvfVWLF26FNXV1W7r8vPzsXLlStx7772uaaA7duxwq7Njx456py60dE0KtWVlZW6P0tJS/Prrr5g9ezbatWvn7TYSERERXTMWL14Mk8mE7OxsfP/99zh9+jTWr1+PW2+9Fa1atXKb6rl161YsXLgQx44dw5IlS7BmzRpMnjw5gK0PnCZNP4iKivI4UUwIgdTUVKxatcorDSMiIiK6FrVr1w579uzB3Llzcc8996CoqAhJSUkYNmwY5s6d63aN2n/84x/Ys2cP5s2bh4iICLz66qvIzs4OYOsDp0mh9ptvvnELtbIsIz4+Hm3btoVa3aRNEhEREZFDeno6li9ffsV6ERER+Pjjj33foCDQpAQ6YMAALzeDiIiIiKjpmjSndsGCBXjvvfc8yt977z289NJLV90oIiIiIqLGaFKofeedd9CxY0eP8i5duuDtt9++6kYRERER0eWdPHkSU6ZMCXQzmo0mhdr8/HwkJyd7lMfHxyMvL++qG0VERERE1BhNCrWpqanYunWrR/nWrVuRkpJy1Y0iIiIiImqMJp0oNnbsWEyZMgUWiwW33HILAGDz5s148skn8Y9//MOrDSQiIiIiupImhdonnngCFy9exPjx42E2mwEAer0e06dPx8yZM73aQCIiIiKiK2lSqJUkCS+99BKefvppHDlyBAaDAe3atYNOp/N2+4iIiIiIrqhJc2qd8vPzUVRUhOuvvx46nQ5CCG+1i4iIiIiowZoUai9evIiBAweiffv2uP32211XPMjJyeGcWiIiIiLyuyaF2qlTp0Kj0SA3NxchISGu8nvvvRfr16/3WuOIiIiIriVjxozBsGHDGlT39OnTeOSRR5CSkgKtVov09HRMnjwZFy9e9G0jm6kmhdqvv/4aL730Elq3bu1W3q5dO5w6dcorDSMiIiKiuh0/fhy9evXCb7/9ho8++gj//e9/8fbbb2Pz5s3IyspCUVFRoJvod006UayystJthNapqKiIJ4sRERFRsyOEgLAoAdm3pJEhSZJXtzlhwgRotVp8/fXXMBgMAIC0tDT06NED119/PZ566im89dZbAIA2bdogJycHhw8fxueff46oqCjMmjULEyZM8GqbAq1Jofamm27CihUr8NxzzwGwXw1BURQsXLgQf/7zn73aQCIiIqKrJSwKzs3ZFpB9pzzbD5JW5bXtFRUVYcOGDZg/f74r0DolJSXhgQcewOrVq7F06VJXmH755Zcxa9YszJs3Dxs2bMDkyZPRvn173HrrrV5rV6A1KdQuXLgQAwcOxJ49e2A2m/Hkk0/i0KFDKCoqqvNOY0RERETkHb/99huEEOjUqVOd6zt16oTi4mKcP38eCQkJAID+/ftjxowZAID27dtj69ateO211xhqu3btimPHjmHx4sUIDw9HRUUFhg8fjgkTJiA5OdnbbSQiIiK6KpJGRsqz/QK2b19ozKVUs7KyPF4vWrTIyy0KrCaFWgCIjIzEU0895c22EBEREfmEJElenQIQSG3btoUkSThy5Ajuvvtuj/VHjhxBdHQ04uPjA9C6wGlwqP3ll18avNEbbrihSY0hIiIiosuLjY3FrbfeiqVLl2Lq1Klu82rz8/OxcuVKjBo1yu3ktB07drhtY8eOHfVOXwhWDQ613bt3hyRJVxzqliQJNpvtqhtGRERERHVbvHgx+vXrh+zsbDz//PPIyMjAoUOH8MQTT6BVq1aYP3++W/2tW7di4cKFGDZsGDZu3Ig1a9bgyy+/DFDrfaPBofbEiRO+bAcRERERNVC7du2wZ88ezJ07F/fccw+KioqQlJSEYcOGYe7cuYiJiXGr/49//AN79uzBvHnzEBERgVdffRXZ2dkBar1vNDjUpqen+7IdRERERNe85cuXN7huenp6g+tHRETg448/blqjgkSTTxT7/fffsWjRIhw5cgQA0LlzZ0yePBnXX3+91xpHRERERNQQTbrGxIYNG9C5c2fs2rULN9xwA2644Qbs3LkTXbp0wcaNG73dRiIiIqJrRm5uLsLCwup95ObmBrqJzVKTRmpnzJiBqVOn4sUXX/Qonz59eou6kC8RERGRP6WkpGD//v2XXd8YJ0+evLoGBYkmjdQeOXIEOTk5HuWPPPIIDh8+3ODtLFiwADfeeCPCw8ORkJCAYcOG4ejRo251jEYjJkyYgNjYWISFhWHEiBEoKChwq5Obm4uhQ4ciJCQECQkJeOKJJ2C1WptyaEREREQBpVar0bZt23ofanWTZ4+2aE0KtfHx8XX+BbF//37X7dgaYsuWLZgwYQJ27NiBjRs3wmKxYPDgwaisrHTVmTp1Kr744gusWbMGW7Zswblz5zB8+HDXepvNhqFDh8JsNmPbtm344IMPsHz5csyZM6cph0ZEREREQahJUX/s2LH4+9//juPHj6NfP/st57Zu3YqXXnoJ06ZNa/B21q9f7/Z6+fLlSEhIwN69e3HzzTejtLQUy5Ytw4cffohbbrkFAPD++++jU6dO2LFjB/r27Yuvv/4ahw8fxqZNm5CYmIju3bvjueeew/Tp0/HMM89Aq9U25RCJiIiIKIg0KdQ+/fTTCA8PxyuvvIKZM2cCsM/veOaZZ/DYY481uTGlpaUA4Lq22t69e2GxWDBo0CBXnY4dOyItLQ3bt29H3759sX37dmRmZiIxMdFVJzs7G+PGjcOhQ4fQo0cPj/2YTCaYTCbX67KyMgCAoihQFKXJ7b+UoigQQnh1m+SJ/ew/7Gv/YD/7B/vZP3zZz/zaUW1NCrWSJGHq1KmYOnUqysvLAQDh4eFX1RBFUTBlyhT0798fXbt2BWC/1ZtWq0VUVJRb3cTEROTn57vq1A60zvXOdXVZsGAB5s2b51FeXFzs1bm4iqKgvLwcQgjIcpNmelADsJ/9h33tH+xn/2A/+4cv+9mZQYiAq7hOrdPVhlmnCRMm4ODBg/jxxx+9sr3LmTlzpts0ibKyMqSmpiI6OhoRERFe24+iKJAkCdHR0fyF6UPsZ/9hX/sH+9k/2M/+4ct+5glTVFujvhuc81qv5JtvvmlUIyZOnIh169bh+++/R+vWrV3lSUlJMJvNKCkpcRutLSgoQFJSkqvOrl273LbnvDqCs86ldDoddDqdR7ksy17/gZMkySfbJXfsZ/9hX/sH+9k/2M/+4at+5teNamvUd8N3332HEydOoHPnzujWrVu9j4YSQmDixIlYu3YtvvnmG2RkZLit79mzJzQaDTZv3uwqO3r0KHJzc5GVlQUAyMrKwoEDB1BYWOiqs3HjRkRERKBz586NOTwiIiKigBozZgwkSYIkSdBoNMjIyMCTTz4Jo9HoqiNJEvR6PU6dOuX23mHDhmHMmDEe29y+fTtUKhWGDh3qse7kyZOu/UmShNjYWAwePBj79u3z+rH5WqNGal966SW8//77WLNmDR544AE88sgjrvmvTTFhwgR8+OGH+Pe//43w8HDXHNjIyEgYDAZERkYiJycH06ZNQ0xMDCIiIjBp0iRkZWWhb9++AIDBgwejc+fOeOihh7Bw4ULk5+dj9uzZmDBhQp2jsURERETN2W233Yb3338fFosFe/fuxejRoyFJEl566SVXHUmSMGfOHHzwwQdX3N6yZcswadIkLFu2DOfOnavz5g2bNm1Cly5dcObMGTz22GMYMmQIfv31V4/zmpqzRo3UPvHEEzh8+DA+++wzlJeXo3///ujduzfefvtt1xUEGuOtt95CaWkpBgwYgOTkZNdj9erVrjqvvfYa7rjjDowYMQI333wzkpKS8Omnn7rWq1QqrFu3DiqVCllZWXjwwQcxatQoPPvss41uDxEREbVMQgiYzeaAPIQQjWqrTqdDUlISUlNTMWzYMAwaNAgbN250qzNx4kT861//wsGDBy+7rYqKCqxevRrjxo3D0KFDsXz58jrrxcbGIikpCb169cI///lPFBQUYOfOnY1qd6A1aYZ1VlYWsrKy8Prrr2PNmjVYsmQJHn/8cZw7d65RJ1o15Ius1+uxZMkSLFmypN466enp+Oqrrxq8XyIiIrq2WCwWvPDCCwHZ96xZs5p83fyDBw9i27ZtSE9Pdyvv378/jh07hhkzZmDdunX1vv/jjz9Gx44d0aFDBzz44IOYMmUKZs6cCUmS6n2PwWAAAJjN5ia1OVCuaob1Tz/9hC1btuDIkSPo2rUrNBqNt9pFREREdE1at24dwsLCoNfrkZmZicLCQjzxxBMe9RYsWID169fjhx9+qHdby5Ytw4MPPgjAPq2htLQUW7Zsqbd+SUkJnnvuOYSFhaF3795XfzB+1OiR2nPnzmH58uVYvnw5ysrK8OCDD2Lnzp08KYuIiIiaLY1Gg1mzZgVs343x5z//GW+99RYqKyvx2muvQa1WY8SIER71OnfujFGjRmHGjBnYunWrx/qjR49i165dWLt2LQD7JdDuvfdeLFu2DAMGDHCr269fP8iyjMrKSlx33XVYvXq1x30AmrtGhdrbb78d3377LQYPHoyXX34ZQ4cO5TXiiIiIqNmTJKnJUwD8LTQ0FG3btgUAvPfee+jWrRuWLVuGnJwcj7rz5s1D+/bt8dlnn3msW7ZsGaxWq9uJYUII6HQ6LF68GJGRka7y1atXo3PnzoiNjQ2qk8Nqa9T0g/Xr1yMmJga5ubmYN28eevfujT/84Q8eDyIiIiK6erIsY9asWZg9ezaqq6s91qempmLixImYNWsWbDabq9xqtWLFihV45ZVXsH//ftfj559/RkpKCj766COP7Vx//fVBG2iBRo7Uzp0711ftICIiIqI6jBw5Ek888YTrxPxLzZw5E++++y5OnDiBe++9F4B9Xm5xcTFycnLcRmQBYMSIEVi2bBkeffRRv7TfXxhqiYiIiJoxtVqNiRMnYuHChRg3bpzH+piYGEyfPt1tzvCyZcswaNAgj0AL2EPtwoUL8csvvzTqqlXNnSQae/G0FqisrAyRkZEoLS316hdXURQUFRUhJiaGt/LzIfaz/7Cv/YP97B/sZ//wZT9f7vPbaDTixIkTyMjIgF6v9+p+yX8a83XkTzERERERBT2GWiIiIiIKegy1RERERBT0GGqJiIiIKOg1OdROnDgRRUVF3mwLEREREVGTNCrUnjlzxrX84YcfoqKiAgCQmZmJ06dPe7dlREREREQN1Kjr1Hbs2BGxsbHo378/jEYjTp8+jbS0NJw8eRIWi8VXbSQiIiIiuqxGjdSWlJRgzZo16NmzJxRFwe2334727dvDZDJhw4YNKCgo8FU7iYiIiIjq1ahQa7FY0Lt3b/zjH/+AwWDAvn378P7770OlUuG9995DRkYGOnTo4Ku2EhERERHVqVHTD6KiotC9e3f0798fZrMZ1dXV6N+/P9RqNVavXo1WrVph9+7dvmorEREREVGdGjVSe/bsWcyePRs6nQ5WqxU9e/bETTfdBLPZjJ9++gmSJOGPf/yjr9pKRERE1KKNGTMGkiRBkiRoNBpkZGTgySefhNFodNWRJAl6vR6nTp1ye++wYcMwZswYj21u374dKpUKQ4cO9Vh38uRJ1/4kSUJsbCwGDx6Mffv2NbjNhw4dwj333IP4+HjodDq0b98ec+bMQVVVVcMP3AsaFWrj4uJw5513YsGCBQgJCcHu3bsxadIkSJKExx9/HJGRkfjTn/7kq7YSERERtXi33XYb8vLycPz4cbz22mt45513MHfuXLc6kiRhzpw5DdresmXLMGnSJHz//fc4d+5cnXU2bdqEvLw8bNiwARUVFRgyZAhKSkquuO0dO3agT58+MJvN+PLLL3Hs2DHMnz8fy5cvx6233gqz2dygNnrDVd18ITIyEvfccw80Gg2++eYbnDhxAuPHj/dW24iIiIi8QggBm60qIA8hRKPaqtPpkJSUhNTUVAwbNgyDBg3Cxo0b3epMnDgR//rXv3Dw4MHLbquiogKrV6/GuHHjMHToUCxfvrzOerGxsUhKSkKvXr3wz3/+EwUFBdi5c+cV+zQnJwedOnXCp59+it69eyM9PR0jR47EF198ge3bt+O1115z1ZckCW+99RaGDBkCg8GA6667Dp988knDOqUBGjWntrZffvkFrVq1AgCkp6dDo9EgKSkJ9957r9caR0REROQNilKN77ZkBmTfA/50ACpVSJPee/DgQWzbtg3p6elu5f3798exY8cwY8YMrFu3rt73f/zxx+jYsSM6dOiABx98EFOmTMHMmTMhSVK97zEYDABwxVHW/fv34/Dhw/jwww8hy+7jpN26dcOgQYPw0UcfYfr06a7yp59+Gi+++CJef/11/N///R/uu+8+HDhwAJ06dbrsvhqiySO1qamprgM4ePAgUlNTr7oxRERERNe6devWISwsDHq9HpmZmSgsLMQTTzzhUW/BggVYv349fvjhh3q3tWzZMjz44IMA7NMaSktLsWXLlnrrl5SU4LnnnkNYWBh69+592XYeO3YMAOoNpJ06dXLVcRo5ciT+9re/oX379njuuefQq1cvvPnmm5fdT0M1eaSWiIiIKFjIsgED/nQgYPtujD//+c946623UFlZiddeew1qtRojRozwqNe5c2eMGjUKM2bMwNatWz3WHz16FLt27cLatWsBAGq1Gvfeey+WLVuGAQMGuNXt168fZFlGZWUlrrvuOqxevRqJiYkNam9jpldkZWV5vN6/f3+D3385DLVERETU4kmS1OQpAP4WGhqKtm3bAgDee+89dOvWDcuWLUNOTo5H3Xnz5qF9+/b47LPPPNYtW7YMVqsVKSkprjIhBHQ6HRYvXozIyEhX+erVq9G5c2fExsYiKiqqQe1s3749AODIkSPo0aOHx/ojR4646vjDVZ0oRkRERES+I8syZs2ahdmzZ6O6utpjfWpqKiZOnIhZs2bBZrO5yq1WK1asWIFXXnkF+/fvdz1+/vlnpKSk4KOPPvLYzvXXX9/gQAsA3bt3R8eOHfHaa69BURS3dT///DM2bdqE+++/3618x44dHq+9MZ8WYKglIiIiatZGjhwJlUqFJUuW1Ll+5syZOHfuHDZt2uQqW7duHYqLi5GTk4OuXbu6PUaMGIFly5ZddbskScKyZctw+PBhjBgxArt27UJubi7WrFmDO++8E1lZWZgyZYrbe9asWYP33nsPx44dw9y5c7Fr1y5MnDjxqtsCMNQSERERNWtqtRoTJ07EwoULUVlZ6bE+JiYG06dPd7tBw7JlyzBo0CC3KQZOI0aMwJ49e/DLL79cddv69euHHTt2QKVSYciQIWjbti1mzpyJ0aNHY+PGjdDpdG71582bh1WrVuGGG27AihUr8NFHH6Fz585X3Q4AkERjL57WApWVlSEyMhKlpaWIiIjw2nYVRUFRURFiYmI8LnVB3sN+9h/2tX+wn/2D/ewfvuzny31+G41GnDhxAhkZGdDr9V7dLzWNJElYu3Ythg0b1uD3NObryJ9iIiIiIgp6DLVERERE5OGHH35AWFhYvY/mhpf0IiIiIiIPvXr18to1ZIHGXc+2KRhqiYiIiMiDwWBwXS83GHD6AREREREFPYZaIiIiIgp6DLVEREREFPQYaomIiIgo6DHUEhEREVHQY6glIiIiClKSJOGzzz4LdDOaBYZaIiIiIgp6DLVEREREzZTZbA50E4IGQy0RERG1eEIIVNpsAXk05k5aAwYMwMSJEzFlyhTExcUhOzu7Ucc5ffp0tG/fHiEhIbjuuuvw9NNPw2KxuNX54osvcOONN0Kv1yMuLg533313o/bRXPGOYkRERNTiVSkKrv/+QED2/fvNmQhVqRpc/4MPPsC4ceOwdevWRu8rPDwcy5cvR0pKCg4cOICxY8ciPDwcTz75JADgyy+/xN13342nnnoKK1asgNlsxldffdXo/TRHDLVEREREzUi7du2wcOHCJr139uzZruU2bdrg8ccfx6pVq1yhdv78+bjvvvswb948V71u3bpdXYObCYZaIiIiavFCZBm/35wZsH03Rs+ePZu8r9WrV+ONN97A77//joqKClitVkRERLjW79+/H2PHjm3y9pszhloiIiJq8SRJatQUgEAKDQ1t0vu2b9+OBx54APPmzUN2djYiIyOxatUqvPLKK646BoPBW81sdniiGBEREVELsG3bNqSnp+Opp55Cr1690K5dO5w6dcqtzg033IDNmzcHqIW+xZFaIiIiohagXbt2yM3NxapVq3DjjTfiyy+/xNq1a93qzJ07FwMHDsT111+P++67D1arFV999RWmT58eoFZ7D0dqiYiIiFqAv/zlL5g6dSomTpyI7t27Y9u2bXj66afd6gwYMABr1qzB559/ju7du+OWW27Brl27AtRi75JEYy6e1kKVlZUhMjISpaWlbpOpr5aiKCgqKkJMTAzkRk4Sp4ZjP/sP+9o/2M/+wX72D1/28+U+v41GI06cOIGMjAzo9Xqv7pf8pzFfR/4UExEREVHQY6glIiIiaoZWrlyJsLCwOh9dunQJdPOaHZ4oRkRERNQM/eUvf0GfPn3qXKfRaPzcmuaPoZaIiIioGQoPD0d4eHigmxE0OP2AiIiIiIIeQy0RERERBT2GWiIiIiIKegENtd9//z3uvPNOpKSkQJIkfPbZZ27rhRCYM2cOkpOTYTAYMGjQIPz2229udYqKivDAAw8gIiICUVFRyMnJQUVFhR+PgoiIiIgCLaChtrKyEt26dcOSJUvqXL9w4UK88cYbePvtt7Fz506EhoYiOzsbRqPRVeeBBx7AoUOHsHHjRqxbtw7ff/89/v73v/vrEIiIiIioGQjo1Q+GDBmCIUOG1LlOCIFFixZh9uzZuOuuuwAAK1asQGJiIj777DPcd999OHLkCNavX4/du3ejV69eAIA333wTt99+O/75z38iJSXFb8dSl2PHjqGiogJGoxGyLEOSJACAJEn1Ll9pfWPqXm79peor92a9hr6XWoam3KzQXzc4bO77URQFNpsNVqu1wXdgau7H1Bz3oygKzGYzTCaTz38/taR+a+x+hBCorq6Goii8cxv5VLO9pNeJEyeQn5+PQYMGucoiIyPRp08fbN++Hffddx+2b9+OqKgoV6AFgEGDBkGWZezcuRN33313nds2mUwwmUyu12VlZQDsv+AURfHaMWzevBnnz5/32vauZQ35wGlOobk5frAQEQXSuHHjEB8f79VtevMzu7n47rvv8Oc//7ne9QMGDMC3337rxxYFj2YbavPz8wEAiYmJbuWJiYmudfn5+UhISHBbr1arERMT46pTlwULFmDevHke5cXFxbBarVfbdJfo6GhYrVaYzWZXmTOEXO65dlCp/bq+52tBQ471WuoPIqJgU1paCpVK5dVtlpeXe3V7zUG/fv2Ql5fnUf7555/j0Ucfxfjx4wPQquDQbEOtL82cORPTpk1zvS4rK0Nqaiqio6MRERHhtf3cc889KC4uRnR0tE//5XK5YFzfckO3GQz1FEVBaWkpIiMj3frZmyHXX6PAzX0/9fW1t/fTWC1tP0IIlJSUICoqyuf/rm1pfdeY/SiK0uTf0c3xeJrrfmr3s7dDrVrd8BgjhEC1xebV/TeUQaNqcB9rtVokJSW5lR05cgSPP/44Zs2ahZEjR/qiiS1Csw21zi9oQUEBkpOTXeUFBQXo3r27q05hYaHb+6xWK4qKijy+IWrT6XTQ6XQe5bIse/0DRJIkn2yXajjnHzY0aFHTKYoCi8WC8PBw9rUPKYqC6upqGAwG9rMPKYoCjUYDrVbLfvYh5+egSqXyej83ZnvVFhs6z9ng1f031OFnsxGibVrkKikpwV133YUBAwbgueee83LLWpZm+1OckZGBpKQkbN682VVWVlaGnTt3IisrCwCQlZWFkpIS7N2711Xnm2++gaIo9d4rmYiIiCgYKIqC//f//h/UajVWrlzZrM4daY4COlJbUVGB//73v67XJ06cwP79+xETE4O0tDRMmTIFzz//PNq1a4eMjAw8/fTTSElJwbBhwwAAnTp1wm233YaxY8fi7bffhsViwcSJE3HfffcF/MoHRERE1HwYNCocfjY7YPtuilmzZmH79u3YtWsXwsPDvdyqliegoXbPnj1uZ/g557mOHj0ay5cvx5NPPonKykr8/e9/R0lJCf74xz9i/fr10Ov1rvesXLkSEydOxMCBAyHLMkaMGIE33njD78dCREREzZckSU2eAhAIq1atwj//+U98+eWXaNeuXaCbExQC+tUdMGDAZU/mkSQJzz77LJ599tl668TExODDDz/0RfOIiIiI/G7//v3IycnBiy++iOzswIwuB6Pg+ZOFiIiIqIW7cOEChg0bhgEDBuDBBx/0uESpSqXy+vV+WwqGWiIiImoUoQjYrAqsFgU2qwKb49lqcSxbFFgd5RazFWUl5Qj7YwT0IdpAN73Z+/LLL3Hq1CmcOnXK7epPTunp6Th58qT/GxYEGGqJiIiCjBCOUGl2PCw2V7CsHShtFvfg6Qqdl61nq1XPGV5tbmFVsTb+OuDX35DCUNsAo0ePxujRowPdjKDEUEtEROQFiiJgNdtcAdFqtjkCp2O5VpnNYoPF8Vw7mLrXrVm2uW3PXobmchNFCVBrZKg0MtRq+7PK8azWyJBVMhTYoFI326uIUgvBUEtERNcEm02B1WSDxWQPiGajBRcvVKIyH7BZBCxmGywmG6xuzwosZpv9fc5nUx3h02Jr0uilN0iyBLVGhlrrHiZVtQKmM3TWBE9VrXoSVOpar+t6j3N7tbfveMiydNnrpyqKgqKiIoTH6OutQ+QNDLVERNRsCCHs8zBNNpiN9gBpMVphNtlgMbqHzpqwqbiHTlcoVdzqKzb/hU6V2h4y1RoZKq3KETpVrvCprqfMHhhVjte1li95VmlkaLQqqLQyVCqOgBIBDLVERHQVhBCwmhWYjVZX6LSYrLUCqc2+zrnsWG8vrwmtFscIqNlog1B8Gz4lWYJGK0OtU0FWAzqDBhqtGhqdPVhqdCqodSpoNCqodbL9taNco7Wvcw+ajnDqKHOOXhKRfzHUEhFdgxRFwGK0wlRtD5imaivM1VaYjVaYq201y0bHcq11bqOnJpvP5naqdSpodY4wqVdBq1c7wqXsCpe1nzU62aNMrfUMpbLa/u9y57/FY2JiIMsc7WypLnc9fGr+GvP1Y6glIgoiQghYLQrMVfZAag+ajiDqWq4VRi8NptU1o6heJaFWAFVDq69ZdoXSWq+1enuZRqd2lKsc5Wr7slYFiaOddBVUKvutac1mMwwGQ4BbQ01VVVUFANBoNFesy1BLRORnNpsCc7UVpkpHMK2ywlhlsZdVWWGstKCspAKSLQ/mavsoqqnKClOVBaZqq1dPSFJpZGj1KmgNamj1amgNKsez41Frnc6gdo2Y1h491Tj+HX+5k4WI/E2tViMkJATnz5+HRqPhaHyQEUKgqqoKhYWFiIqKcv2RcjkMtURETWCzKDBWWWCssMBUZYGx0h48zdWOgOoYSXWFUcc6U5XVK6OkkizVBFBnGHUFUzV0BvuoqO6SYGoPq47XOjVUGn7QU8skSRKSk5Nx4sQJnDp1KtDNoSaKiopCUlJSg+oy1BLRNc1mVeyjoxUWj5BqrLTAVGmBsfKS11VWWL0QTDV6FXQGNXQhGuhC1PaHQQ2NQQVFsiIqJhz60FrrQjT2wBpiHx3lyCjR5Wm1WrRr1w5msznQTaEm0Gg0DRqhdWKoJaIWQQgBc7UV1RX2YFpdUTuQOkZSa702OULq1YyaShKgC9V4BE+dQQ1trZDqFlpD1NAZNNAaVJDruRQTT2Ai8h5ZlqHX8xq51wKGWiJqlmxWxRFOzfagWm5xBFZzreBqtj+X218rTb0UlAToQtTQuwKqBvowNfQhGujDHK9D7evtIda+rNWreTITEVEzwVBLRD4nhIDZaEN1mRnGSguqy81uI6rGcjOqKx2vy+1B1Wxs2giqRqeCIdweTl2BNEwDfYi6VkB1hFdnUDUwnBIRBTuGWiJqEkURrhBaVWZ/VJebXa+ryy21yiywWZVG70OSJehD1TCEa6EP1djDapgWhjB7UDWEa2AI1UIfrnGVqTUNn39FREQtB0MtEblYLTa3MOoKpWUWVJWbUVVmQkWJEaZKG4yVlkZfdF+tU8EQ5gygWkdI1TjKtDXL4fZljqASEVFDMdQStXCKTUF1uQWVpSZUlZrtz2VmVJWaUVVuRnVZzXOj/+UvAfpQDUIi7AE1JFwLQ7gWhgito0zrKNPAEKGFRstRVCIi8g2GWqIgZbMoqCyrFVTdnu2jqpWlZhjLzWjMXSJlleQKpPZQqnGFVH2oGlbJhISUWIRG6mAI09R7Bj8REZE/MdQSNTNmo9U+iuoIpR6htcz+bKq0NnibkgQYIrQIjdQhJFKL0AgtQiJ1NaOpERrHsxZag7re65/WXGoqjJeaIiKiZoWhlsgPhBAwVVld4bSqtFZgvWS0tTHXTZXVEkIjHEHVGVgjawKrs8wQroXMualERNSCMdQSXQVFEfYTqtzmqnqOsFaVmRt19r9ap3KMptYOq57PupD6R1WJiIiuJQy1RHWwWRXXv/ndR1ZNqCyrCazV5RaIRlzwXxeiRkikzjGaqq1jlNX+rNXzR5OIiKgx+MlJ15S65qvWLNeUGSstDd+oBBjCHf/2j6gVWGsHVceoK6+hSkRE5BsMtRT0hCJgrLTiYlWF49JV9qBa+yoA9rmrZlgbM1/VcRWAmpFVx3OEe2A1hPMKAERERIHGUOtDz2x/BqdLTmNm1ky0j2kf6OYEldq3VXXekaqq3H06gH16gP36qkojpgBodCrXyGntUdRLpwDoQzS88D8REVGQYKj1oX2F+3Cy7CTKTGWBbkqz4LxblTOkVpc7LvrvKnMvV6yNu12VPkxTM6LqCKohEZ5TADhflYiIqOXhp7sPqSQVYipTcGpTFaT4kwiN0tkfkTqEhGvtZ64H2UigEAI2iwKz0QZjhQXGKgtMlRYYKy0wVlrty1VWGCssMFU5yy0wVVobdakqJ41eVXMDgPDal6qqmQ6gD1PDaK1EXHwsr51KRER0jWKo9SG1pMbQI48i32JFPo57rJdkye1e94YwDbR6FTR69SXPKmi0KsgqGbJagkol2ZdrPUsSIIQ9dMLxLARcy4rNHkatVgU2i/1htSiwWWz2Z6sCq1mBxWiD2WSFudoGyyXPZqMVFqOtUf/qv5Sskhx3qqp1S1VHYHWVO+9mFaaBugG3VVUUBeaiqia3iYiIiIIfQ60PhRgjEWqJBAC0uzER1eVmVJbYz7A3V1shFGGfM1pmBlAZ2MY2lmS/PJU+RANdqAb6UA30oWr7coga+jANdCH2cl2o2rFew+uqEhERkU8w1PpQRHkCAECXCAzO6eK2zmZVYKywoLrCjOoyx3O5xT4qarR5jJhaTAoUmwLFJhwP+7LNJqA4L+ovAZJkH7V1PTumN8gqCSqNCmqNDJVahkoj25cdr53LtUeItXo1NI5nt2WDfeQ42KZOEBERUcvFUOtDOnMIAEAd4fnvepVads2xJSIiIqKrw7NqfEjl/JtBbvocVCIiIiK6MoZaH1LBcZKTpAS2IUREREQtHEOtD8nCHmoFR2qJiIiIfIqh1odUwj79QHCkloiIiMinGGp9SHZ0r5A4UktERETkSwy1PuSafsCRWiIiIiKfYqj1IYZaIiIiIv9gqPWhmlDL6QdEREREvsRQ60M1c2o5UktERETkSwy1PuQcqVUYaomIiIh8iqHWh2TBkVoiIiIif2Co9SGeKEZERETkHwy1PiQ5upfTD4iIiIh8i6HWhyThDLW2ALeEiIiIqGVjqPUhmaGWiIiIyC8Yan3INVILTj8gIiIi8iWGWh/i9AMiIiIi/2Co9SFJkQAw1BIRERH5GkOtDzlHam1gqCUiIiLyJYZaH5IER2qJiIiI/EEd6Aa0aI5Qy5FaaihFEbAoCiw2AatNgdlmX7ZYFVgVBWargMWmwOJYZ1OE20MRAtZayzYFsCmK/VkI2GwKbMK+H6urTq2HEK51NkVACAEBQAhAEfb6RqMJOt05CAEICMezvU5NffuzUqsMjm3Ufp/iWFGzDeHWH5Ik1Sy7ytxq1KrrWepeVk9dyXN9XYuSZK8hSYDsWEatZVe5Y4Us2csk17J9S846EgBZrns7EgBjtRGGkItQyRIkSI7t17Md137r276jDbW3A+ky7a+9nZpjQa1jsW+/7u043wO3fda0S5Jq96fkcQzu9WvaDHges3O9+zbdv0bSJe121hVCQYXJCq3JCrVKdm3r0n6UJPevIRE1Twy1PiQpPFEs2FhsCqpMNlRZrDBaFBgtNpis9T+b6iq3KDBaPZ8tNgUWqzO0OpYd4dRqsy9bFXHlRhJRQNUOzR7BunaIBxyhuuYPDLcwXyvky7XKgZo/GGqHeFwasusK4W5/VHiG/Mu21/VHhXP7Nct1/kHhti37m2v/YeY8VgAwm014+i+hSIw0+PJLQ9c4hlofkjhS6zNCCFSabagwWlFWbca58xVQFSmotiiotthQZbah0mRFtdmGSrMN1WYrqsw2x8PqKLOh0myv4yy32JpXqJQkQKOSoVXJ0KgkaFSy4yFBrZKhliWoVRJUkgRZtj+r5EsejnVquY46tdapZAmyZN+eLElQybU/mB0fbELAaDQixGCwf0jXMZpW+4P20g9HyWPZ80PbRdRetL8Qos7Vl5Rfvm7tFcKzyG20uPb7nKPOzvrKJaPYzvfa69nboTgq2Eeta0akFceCqOM9wjGSXm00QqvTudohnO+/ZN+1R8eVWtts0nvcRtHdR9DFpe11dIRyych73X1U175r9uHWd/WM4rvac7ltOTrL2fe113uTEPb/fNQq8e4OWqgp2VYkBroR1KIx1PqS4+oHVlgC3JDmQwgBo0VBudGCcpMVFUYryo1WVJgsKHctW1FutDiea8rsdR3vM1m9/kFVm1qWYNCooNOooFPL0Gtk6NSqyz9rVNCrZdd7Ln2tVTvDqWdA1ajs69WyBE2teiq5ef27U1EUFBUVISYmBrLMKfm+wn72vrqCuU2xoaioGFFR0ZBkySNA1xeQ3ZZrr3f8Trr0D4LaAd2mCMd77NN7bLX/mHBMAXJuw1nHuXxpuTNcO7fnnOLjXK69rib41/xxULuu4vHHjHD7Q0KBgFBqHVOt7TjbfOkfTK7tKPY/hsN0fvyC0zWpxYTaJUuW4OWXX0Z+fj66deuGN998E7179w5omzSSBlYA542FAW2Ht5itiitwOsOmc9lV7hZU7ctljoDqfO3Nf7GrZQnhejUMahlhBg1CdWqEaFUwaNQI1akQolUhROso06oQqlW7ni8tC6m1rFUHZ5AQrg8tx3Pt5fqencsNeI+iKCiqKkK1tto+J7Gp24PzA1hxqwPUKnPUravs0mdXnUu3jzrKrrBPwPmBrjSsXbVe2z/s6yirtc+aoHSZfSoC1cZqaHVaQIJnnXr26axTe5+1vy9cZc4Q5ShrSB23Msc+a9dx7rN2Gy77/VH7++qSNtT3fue2nXVd+7p0vXP0t651l/ysKIpS871c80a317X7+9LXbt9ftd57aT8QkCP+DSA80M2gFqxFhNrVq1dj2rRpePvtt9GnTx8sWrQI2dnZOHr0KBISEgLSJiEE1NBBiAqcqz6HF3a8CL3aAFmS3U9IcairzLGiQXXrfT8kWG0SLFYJZqsEi1WG2QKYrBKMFgkmswyzRYLJIsPkeDZbZPuzVXa9Nltk2BTvBT0JAmq1FRqVBRq11f5QWaBWWaBRO55VZqhVVqjVZqhli/3ZUaZS2ZdlyQohCVgtFqjUKtdHiSIEqiFQBQFYAMWiQFQ4P2wu/YB1//Bz1HB9ONpq1QEkKKj5N6n9SOxLihAQrqAn1doOPLbrLAMABXCtcwsJjr5ynXRVU+uSdsOtfkN6P3jUdVQNO1KpwT3SgH3Uuammt83b7w2mrygFlv38vZr5rzXljjm3zteuObK15srWqotL6tbeLi6pKwGAAGRjuTcPhciDJC493TgI9enTBzfeeCMWL14MwP6vu9TUVEyaNAkzZsy44vvLysoQGRmJ0tJSREREeKVNZRVVWPDUEqRFqlEdBRjVFtgkwCYkKJIEBZL9NQAbJCgAFEmGRUiwQoJVyPYHZFgVCVbIsAkJNiHD6nqWYVU8HzZFhtUmw6ZIjiDq3Y88SRaQ1QKyyvkMSGpAVglIavuypKr9LAFqCVBJgEaCpJIgVICQZQhINQ9JggLZ7XXNegCSM6bYj8e5XPtZSM51NccsIHk+S5dZB0BIwTlSS1Qfqc5f9Q0LyHX9gVDnb5U69tHQ7TV1n3Uf11Vsr662iSvXuap+q6teHc1o+H6v/L66iq6qjxpQb3mrEHTr2qdB720oX3x+U/AK+pFas9mMvXv3YubMma4yWZYxaNAgbN++vc73mEwmmEwm1+uysjIA9jCsKIpX2mU1WbA6tjOUKgDNZPaBUNmDpVBJgFoG1BJEXc8a2V5HI0Oo7XWdz1A7TnkloqAi6vy59Sy7qlEO/mqgyygqPuO1z1gnb2+PglvQh9oLFy7AZrMhMdH9nMrExET8+uuvdb5nwYIFmDdvnkd5cXExrFarV9qlmGyQZcDmGCh1XuOx9v9wJOef4rVyoiQ7H8LxsL+WJeEqkyUBWRaQZUClUiDLCtQqBSqVgEpWoHK8VqsUqGUFapUNakmBShKQhIAM4RoPlYXb2CgkAcjO8VLHsmQTkG0CkklAFs56tbdR834ZiuM1XK9ddYWzjvP9iv2sd7d1ittrZ30IOF4DzvHUmpEMUbtbHe+tXV7zMe3ahnBfV/vZ9blcRx33/dexXdc+lVptkdCgqNCkQNDUCNLY9zVhP3UNNaHOwbxa76n3Rb2tqHNzdfSlaOD2Ln2zq06t46lrW3Wpu711zinyoMBzUtFltyeJeus0vC+v4rga+Mfupfuot72XfP807LhEk7/OAJyXFr/kvVf+el1dXzb9D4s63+tRdBXbr+Nr2tTvLwEgOnUUioqKGrj3hikv55QGqhH0obYpZs6ciWnTprlel5WVITU1FdHR0V7998Whx/+EklOliL0+kWcw+5CiKCguLkZ0dDT72cfY1/7BfvYP9rN/+LKf1eprMsZQPYL+uyEuLg4qlQoFBQVu5QUFBUhKSqrzPTqdDjqd57VFZFn26g+cRhsCVazR69slT5IksZ/9hH3tH+xn/2A/+4ev+plfN6ot6L8btFotevbsic2bN7vKFEXB5s2bkZWVFcCWEREREZG/BP1ILQBMmzYNo0ePRq9evdC7d28sWrQIlZWVePjhhwPdNCIiIiLygxYRau+9916cP38ec+bMQX5+Prp3747169d7nDxGRERERC1Tiwi1ADBx4kRMnDgx0M0gIiIiogAI+jm1REREREQMtUREREQU9BhqiYiIiCjoMdQSERERUdBjqCUiIiKioMdQS0RERERBj6GWiIiIiIIeQy0RERERBT2GWiIiIiIKei3mjmJXQwgBACgrK/PqdhVFQXl5OdRqNWSZfz/4CvvZf9jX/sF+9g/2s3/4sp+dn9vOz3G6tjHUAigvLwcApKamBrglRERE1Fjl5eWIjIwMdDMowCTBP2+gKArOnTuH8PBwSJLkte2WlZUhNTUVp0+fRkREhNe2S+7Yz/7DvvYP9rN/sJ/9w5f9LIRAeXk5UlJSONpOHKkFAFmW0bp1a59tPyIigr8w/YD97D/sa/9gP/sH+9k/fNXPHKElJ/5ZQ0RERERBj6GWiIiIiIIeQ60P6XQ6zJ07FzqdLtBNadHYz/7DvvYP9rN/sJ/9g/1M/sITxYiIiIgo6HGkloiIiIiCHkMtEREREQU9hloiIiIiCnoMtUREREQU9BhqfWjJkiVo06YN9Ho9+vTpg127dgW6SUFjwYIFuPHGGxEeHo6EhAQMGzYMR48edatjNBoxYcIExMbGIiwsDCNGjEBBQYFbndzcXAwdOhQhISFISEjAE088AavV6s9DCSovvvgiJEnClClTXGXsZ+85e/YsHnzwQcTGxsJgMCAzMxN79uxxrRdCYM6cOUhOTobBYMCgQYPw22+/uW2jqKgIDzzwACIiIhAVFYWcnBxUVFT4+1CaLZvNhqeffhoZGRkwGAy4/vrr8dxzz6H2OdHs58b7/vvvceeddyIlJQWSJOGzzz5zW++tPv3ll19w0003Qa/XIzU1FQsXLvT1oVFLIsgnVq1aJbRarXjvvffEoUOHxNixY0VUVJQoKCgIdNOCQnZ2tnj//ffFwYMHxf79+8Xtt98u0tLSREVFhavOo48+KlJTU8XmzZvFnj17RN++fUW/fv1c661Wq+jatasYNGiQ2Ldvn/jqq69EXFycmDlzZiAOqdnbtWuXaNOmjbjhhhvE5MmTXeXsZ+8oKioS6enpYsyYMWLnzp3i+PHjYsOGDeK///2vq86LL74oIiMjxWeffSZ+/vln8Ze//EVkZGSI6upqV53bbrtNdOvWTezYsUP88MMPom3btuL+++8PxCE1S/PnzxexsbFi3bp14sSJE2LNmjUiLCxMvP7666467OfG++qrr8RTTz0lPv30UwFArF271m29N/q0tLRUJCYmigceeEAcPHhQfPTRR8JgMIh33nnHX4dJQY6h1kd69+4tJkyY4Hpts9lESkqKWLBgQQBbFbwKCwsFALFlyxYhhBAlJSVCo9GINWvWuOocOXJEABDbt28XQth/CcuyLPLz81113nrrLRERESFMJpN/D6CZKy8vF+3atRMbN24Uf/rTn1yhlv3sPdOnTxd//OMf612vKIpISkoSL7/8squspKRE6HQ68dFHHwkhhDh8+LAAIHbv3u2q85///EdIkiTOnj3ru8YHkaFDh4pHHnnErWz48OHigQceEEKwn73h0lDrrT5dunSpiI6Odvu9MX36dNGhQwcfHxG1FJx+4ANmsxl79+7FoEGDXGWyLGPQoEHYvn17AFsWvEpLSwEAMTExAIC9e/fCYrG49XHHjh2Rlpbm6uPt27cjMzMTiYmJrjrZ2dkoKyvDoUOH/Nj65m/ChAkYOnSoW38C7Gdv+vzzz9GrVy+MHDkSCQkJ6NGjB959913X+hMnTiA/P9+tryMjI9GnTx+3vo6KikKvXr1cdQYNGgRZlrFz507/HUwz1q9fP2zevBnHjh0DAPz888/48ccfMWTIEADsZ1/wVp9u374dN998M7RaratOdnY2jh49iuLiYj8dDQUzdaAb0BJduHABNpvN7UMeABITE/Hrr78GqFXBS1EUTJkyBf3790fXrl0BAPn5+dBqtYiKinKrm5iYiPz8fFedur4GznVkt2rVKvz000/YvXu3xzr2s/ccP34cb731FqZNm4ZZs2Zh9+7deOyxx6DVajF69GhXX9XVl7X7OiEhwW29Wq1GTEwM+9phxowZKCsrQ8eOHaFSqWCz2TB//nw88MADAMB+9gFv9Wl+fj4yMjI8tuFcFx0d7ZP2U8vBUEvN3oQJE3Dw4EH8+OOPgW5Ki3P69GlMnjwZGzduhF6vD3RzWjRFUdCrVy+88MILAIAePXrg4MGDePvttzF69OgAt67l+Pjjj7Fy5Up8+OGH6NKlC/bv348pU6YgJSWF/UzUwnH6gQ/ExcVBpVJ5nCFeUFCApKSkALUqOE2cOBHr1q3Dt99+i9atW7vKk5KSYDabUVJS4la/dh8nJSXV+TVwriP79ILCwkL84Q9/gFqthlqtxpYtW/DGG29ArVYjMTGR/ewlycnJ6Ny5s1tZp06dkJubC6Cmry73eyMpKQmFhYVu661WK4qKitjXDk888QRmzJiB++67D5mZmXjooYcwdepULFiwAAD72Re81af8XUJXi6HWB7RaLXr27InNmze7yhRFwebNm5GVlRXAlgUPIQQmTpyItWvX4ptvvvH4l1TPnj2h0Wjc+vjo0aPIzc119XFWVhYOHDjg9ot048aNiIiI8AgX16qBAwfiwIED2L9/v+vRq1cvPPDAA65l9rN39O/f3+OydMeOHUN6ejoAICMjA0lJSW59XVZWhp07d7r1dUlJCfbu3euq880330BRFPTp08cPR9H8VVVVQZbdP9pUKhUURQHAfvYFb/VpVlYWvv/+e1gsFledjRs3okOHDpx6QA0T6DPVWqpVq1YJnU4nli9fLg4fPiz+/ve/i6ioKLczxKl+48aNE5GRkeK7774TeXl5rkdVVZWrzqOPPirS0tLEN998I/bs2SOysrJEVlaWa73zUlODBw8W+/fvF+vXrxfx8fG81NQV1L76gRDsZ2/ZtWuXUKvVYv78+eK3334TK1euFCEhIeJf//qXq86LL74ooqKixL///W/xyy+/iLvuuqvOyyL16NFD7Ny5U/z444+iXbt21/Slpi41evRo0apVK9clvT799FMRFxcnnnzySVcd9nPjlZeXi3379ol9+/YJAOLVV18V+/btE6dOnRJCeKdPS0pKRGJionjooYfEwYMHxapVq0RISAgv6UUNxlDrQ2+++aZIS0sTWq1W9O7dW+zYsSPQTQoaAOp8vP/++6461dXVYvz48SI6OlqEhISIu+++W+Tl5blt5+TJk2LIkCHCYDCIuLg48Y9//ENYLBY/H01wuTTUsp+954svvhBdu3YVOp1OdOzYUfzP//yP23pFUcTTTz8tEhMThU6nEwMHDhRHjx51q3Px4kVx//33i7CwMBERESEefvhhUV5e7s/DaNbKysrE5MmTRVpamtDr9eK6664TTz31lNtlotjPjfftt9/W+Tt59OjRQgjv9enPP/8s/vjHPwqdTidatWolXnzxRX8dIrUAkhC1brNCRERERBSEOKeWiIiIiIIeQy0RERERBT2GWiIiIiIKegy1RERERBT0GGqJiIiIKOgx1BIRERFR0GOoJSIiIqKgx1BLREREREGPoZaImp0xY8Zg2LBhAdv/Qw89hBdeeKFBde+77z688sorPm4RERFdCe8oRkR+JUnSZdfPnTsXU6dOhRACUVFR/mlULT///DNuueUWnDp1CmFhYVesf/DgQdx88804ceIEIiMj/dBCIiKqC0MtEflVfn6+a3n16tWYM2cOjh496ioLCwtrUJj0lb/97W9Qq9V4++23G/yeG2+8EWPGjMGECRN82DIiIrocTj8gIr9KSkpyPSIjIyFJkltZWFiYx/SDAQMGYNKkSZgyZQqio6ORmJiId999F5WVlXj44YcRHh6Otm3b4j//+Y/bvg4ePIghQ4YgLCwMiYmJeOihh3DhwoV622az2fDJJ5/gzjvvdCtfunQp2rVrB71ej8TERPz1r391W3/nnXdi1apVV985RETUZAy1RBQUPvjgA8TFxWHXrl2YNGkSxo0bh5EjR6Jfv3746aefMHjwYDz00EOoqqoCAJSUlOCWW25Bjx49sGfPHqxfvx4FBQW455576t3HL7/8gtLSUvTq1ctVtmfPHjz22GN49tlncfToUaxfvx4333yz2/t69+6NXbt2wWQy+ebgiYjoihhqiSgodOvWDbNnz0a7du0wc+ZM6PV6xMXFYezYsWjXrh3mzJmDixcv4pdffgEALF68GD169MALL7yAjh07okePHnjvvffw7bff4tixY3Xu49SpU1CpVEhISHCV5ebmIjQ0FHfccQfS09PRo0cPPPbYY27vS0lJgdlsdptaQURE/sVQS0RB4YYbbnAtq1QqxMbGIjMz01WWmJgIACgsLARgP+Hr22+/dc3RDQsLQ8eOHQEAv//+e537qK6uhk6nczuZ7dZbb0V6ejquu+46PPTQQ1i5cqVrNNjJYDAAgEc5ERH5D0MtEQUFjUbj9lqSJLcyZxBVFAUAUFFRgTvvvBP79+93e/z2228e0wec4uLiUFVVBbPZ7CoLDw/HTz/9hI8++gjJycmYM2cOunXrhpKSEledoqIiAEB8fLxXjpWIiBqPoZaIWqQ//OEPOHToENq0aYO2bdu6PUJDQ+t8T/fu3QEAhw8fditXq9UYNGgQFi5ciF9++QUnT57EN99841p/8OBBtG7dGnFxcT47HiIiujyGWiJqkSZMmICioiLcf//92L17N37//Xds2LABDz/8MGw2W53viY+Pxx/+8Af8+OOPrrJ169bhjTfewP79+3Hq1CmsWLECiqKgQ4cOrjo//PADBg8e7PNjIiKi+jHUElGLlJKSgq1bt8Jms2Hw4MHIzMzElClTEBUVBVmu/1ff3/72N6xcudL1OioqCp9++iluueUWdOrUCW+//TY++ugjdOnSBQBgNBrx2WefYezYsT4/JiIiqh9vvkBEVEt1dTU6dOiA1atXIysr64r133rrLaxduxZff/21H1pHRET14UgtEVEtBoMBK1asuOxNGmrTaDR48803fdwqIiK6Eo7UEhEREVHQ40gtEREREQU9hloiIiIiCnoMtUREREQU9BhqiYiIiCjoMdQSERERUdBjqCUiIiKioMdQS0RERERBj6GWiIiIiIIeQy0RERERBb3/D0tkWmVK3LwMAAAAAElFTkSuQmCC", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "plt.plot(ts, res, label=om.vmap.keys())\n", - "plt.grid(alpha=0.25)\n", - "plt.legend(loc=\"center left\", bbox_to_anchor=(1, 0.5))\n", - "plt.xlabel('Time (s)')\n", - "plt.ylabel('# Molecules / cell')" - ] - }, - { - "cell_type": "code", - "execution_count": 10, "id": "f8742e84", "metadata": {}, "outputs": [ @@ -260,7 +223,7 @@ "Text(0, 0.5, '# Molecules / cell')" ] }, - "execution_count": 10, + "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, @@ -284,39 +247,6 @@ "plt.ylabel('# Molecules / cell')" ] }, - { - "cell_type": "code", - "execution_count": 11, - "id": "e29ef3bd", - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([[1.00000000e+00, 0.00000000e+00, 5.00000000e+01, ...,\n", - " 0.00000000e+00, 0.00000000e+00, 0.00000000e+00],\n", - " [1.00000000e+00, 1.99201532e-02, 6.07414794e+00, ...,\n", - " 1.72803021e-01, 5.07038542e-03, 2.48748326e-04],\n", - " [1.00000000e+00, 3.96419002e-02, 2.44132100e+00, ...,\n", - " 1.76855645e-01, 1.02367562e-02, 1.00575004e-03],\n", - " ...,\n", - " [1.00000000e+00, 1.99990712e+00, 4.54120600e+00, ...,\n", - " 3.58319690e-01, 1.10834085e+00, 8.19104745e+00],\n", - " [1.00000000e+00, 1.99990804e+00, 4.54279162e+00, ...,\n", - " 3.58230823e-01, 1.10800706e+00, 8.18857504e+00],\n", - " [1.00000000e+00, 1.99990896e+00, 4.54437502e+00, ...,\n", - " 3.58142129e-01, 1.10767393e+00, 8.18610758e+00]])" - ] - }, - "execution_count": 11, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "res" - ] - }, { "cell_type": "markdown", "id": "08c1353b", @@ -327,7 +257,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 10, "id": "a1d9cfef", "metadata": {}, "outputs": [], @@ -338,14 +268,6 @@ " with open(f'scenario5_{fkey}.json', 'w') as fh:\n", " json.dump(fun(tm), fh, indent=1)" ] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "d9800a23", - "metadata": {}, - "outputs": [], - "source": [] } ], "metadata": {