{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "To obtain the data to run this tutorial, download the tar.gz file from https://zenodo.org/records/10668525 and unzip it. Then run the creating_h5 tutorial to create a .h5 file to use with this tutorial. Then set example_h5_filepath to be the path to that .h5 file.\n", "\n", "Once you have run the tutorial on the sample simulation, you can also try downloading a file from the MAYA catalog in the MAYA format (https://cgp.ph.utexas.edu/waveform) and using it with the tutorial." ] }, { "cell_type": "markdown", "metadata": { "collapsed": true, "jupyter": { "outputs_hidden": true }, "pycharm": { "name": "#%% md\n" } }, "source": [ "# Compact Objects" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "from mayawaves.coalescence import Coalescence\n", "from mayawaves.compactobject import CompactObject\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Create a Coalescence object using the simulation h5 file" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "example_h5_filepath = \"D11_q5_a1_-0.362_-0.0548_-0.64_a2_-0.0013_0.001_-0.0838_m533.33.h5\"" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "coalescence = Coalescence(example_h5_filepath)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Obtain the CompactObject objects associated with the larger object (primary) and smaller object (secondary)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "primary_compact_object = coalescence.primary_compact_object\n", "secondary_compact_object = coalescence.secondary_compact_object" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Obtain trajectory data" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "time, primary_position = primary_compact_object.position_vector\n", "time, secondary_position = secondary_compact_object.position_vector" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASAAAAEKCAYAAACytIjQAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAACD/klEQVR4nO1dZXgUZxc9bzwkhOAS3N3dXUuRlhZa6t5Sd+9Xd6WlpUYLtJQWKS3ursEdAiEJMeIuu/t+P85MVzK72U02yQbmPM8+m+zOzs7szNy5cu65QkoJHTp06KgIeFX0BujQoePahW6AdOjQUWHQDZAOHToqDLoB0qFDR4VBN0A6dOioMOgGSIcOHRWGCjVAQoifhBCJQojjFq/VEEKsF0KcU56r2/nsWCHEGSHEeSHEC+W31Tp06HAXKtoDmgdgrM1rLwDYKKVsBWCj8r8VhBDeAL4GMA5AewAzhBDty3ZTdejQ4W5UqAGSUm4DkGLz8iQAvyh//wJgssZHewM4L6W8IKUsALBI+ZwOHToqEXwqegM0UFdKGQcAUso4IUQdjWXCAERb/B8DoE9xK65Vq5Zs2rSpWzZShw4d2ggPD0+SUtZ2ZllPNEDOQGi8ptlTIoS4H8D9ANC4cWMcOHCgLLdLh45rHkKIS84uW9E5IC0kCCHqA4DynKixTAyARhb/NwQQq7UyKeVcKWVPKWXP2rWdMso6dOgoJ3iiAVoB4A7l7zsA/K2xzH4ArYQQzYQQfgCmK5/ToUNHJUJFl+F/B7AbQBshRIwQ4h4A7wMYJYQ4B2CU8j+EEA2EEKsAQEppADALwFoApwAsllKeqIh90KFDR8lRoTkgKeUMO2+N0Fg2FsB4i/9XAVhVRpumQ4eOcoAnhmA6dOi4RqAbIB06dFQYdAOkQ4eOCkNl5QHpqOwwFgJZCUBuKpCbBuSl8bkgGzAZAFMhYDTwby9vwNsP8AkAfPwAn0AgqDZQtS4QXA+oUhPw0u+llRG6AdJRdjAZgeQIIP4oEHcESLsEpF8GMi4DmfGwwx11HV4+QFAds0GqWheo2gCo3Rqo0x6o0Rzw9nXPd+lwK3QDpMM9MBQAV06bjU3cESD+OFCYzfe9/YDQJkC1MKDFCD5XrQ9UqQEEhAKBoXz2Cwa8fQAvXxoNLx8aMmM+YFAehTlA9hUasawE6+f0aCBmP5CTZN42L1+gViugbkegYS+gYQ+gbid6UzoqFLoB0lEyZMQBFzYDUXtobBJPAsYCvucXDNTrBHS/DajXGajfBajdxuyFGApoMLISgMw4IOksQ7H8TKAgC8jPYigmjYCUACSfvbwB30DANwjwCwL8qgCB1bnupgOB4Dr0gIJqAYW5XO+V00DiKT5f3AYcW8xt8PbndjUfArQcBYT1oOHTUa4Q19JYnp49e0q9F6yEMBYCl3YC5zcA5zcBiQrvM7A6L2TV0NTvypDHywvISQGO/UXjkpUAJJ8HUi7QS9GCtz/gH0wD5hdMgyMEAMFnk5GGpTCHhqogh7kiW/iHADVbADVbATVbArVa8rlGcyAvHYg5QC8pei9wORyQJnpfLYYBLUfyUbVeGf2QVz+EEOFSyp5OLasbIB12YSgALm4FTi4HTq+kIfH2Axr3ZRjVYjjDGi8vei1xR4DYg8Dlg7yw0yx6Ev2r0RDUaAHUaAaENGAIFlzXHIq5mqeRkl5T9hUgK9HsUSVH0Nglny9q7Op3ZRjWqDefA6oBF7YA5zfSuGbFWyzbBZgyF6jTtoQ/4LUJ3QDZgW6AnICU9AwOLQBOrgDy0+lRtBkHtJ8ENB/K8MeQz+UiNjMUiztCTwIAqjUGwroBDbrT82jYi+GR0BIxKGMU5tLrurCFBsZYSAOp5qaEF7d75BtAu+vpXZ1bD2z8n3kddTsBnacBHW8AqjUs/32oZNANkB3oBsgBMuKAI78DhxfSc/ANAtpfD3SYQqPj48+Q6vRKYMUs8+e8fGhgmg0GwnoCDboBwR6uOmAyMmcVvQ84+Te9PBU1mjMEazMOqNUGOPUPcOxP4PIBAAJoOQLoeTfQaoyeM7ID3QDZgW6ANBBzANjzDS9EkwFo3B/odivQfjLzMTkpwOl/gRPLeaGaDPxc/S7A0BeZ/PWvWpF74B4kRyhh2Hrg4nbAkEt+UdvraISrNWQ+6+CvQGYsEBIGdL8D6HGHni+ygW6A7EA3QAqMBuDU38CeOUzG+ocA3W/nnb1mC4ZXp/4BDv9mNjqhTYAOk3kx1u9aNuFUfhaQn8GwqSCb4ZAhD0xCe5kfPv7M3fhX5cMnwL3bU5hLY3RiGXBmNcO1oNpAl+lAl1sY0h34CYjYyMR599uAAY8DoY3dtw2VGLoBsoNr3gAZC4Gji4HtH/MiqtEc6PMQ0HUGL+Skc0D4PBqe3BTmcjpOpeEprdExFgKpl4CMGDMZMT1Geb4MZMQy31QSePkymR1Sn95I1QZAaCOgVmvmoEKblDxcKsxlTujoH8DZNTTGjfrQYDfoBuz9jr8XJNB5OjDwSSbbr2HoBsgOrlkDZCxkfmf7J0BqJEvmQ54D2kygUbmwBZg/mct6+QBtxgM97gSaDytZi0NBDqthcUeB+GNAwjHgyhkzT0hFlVokJIY0VJ4b0LPxVTg+vlXo3ag8IGkiN8iQz+pXfgafc9NYAcuIJRkxM47vqfD2Y/WtQVegYU/mqup2cL3qlpXI3/HgfCD5HEO0nnfz9zr6B423sYDh66CngXodXf/trgLoBsgOrjkDJCWTxutfA1IieMce8jzQeiwv5lP/ADs+A+IOc/mwHsCMRaxYuYKcFHKEovYAUbtZEVNzRUF1SEqs1xGo3Za5lJAwPnwD3Lq7RbYp6RwNRdJZGsDL4SzZA+wna9CVFIAG3WhwA0KcW7eU3N/d3wBnVtGQdbyRubNz64H9P5CnBACPHSbt4BqCboDs4JoyQJfDgbWvAFG7WM0Z9T8aHmMB7+I7v6RRqtGc+YvO0503CCYTEH8EOLeBSduY/TRo3v40Yo378tGgm+vGrCwhJZAWxYpWjPrYx/e8fIBGfYGWw1kFq9fZuZAzOYK5tMMLma/qdBPQ5wFg9fNct7cf/x/0DNtNrgHoBsgOrgkDlHUFWP8qjUxQbWDYS0C323kxHf0D2PQO8zD1uzJf0W4iGcfFwWTiBXV8KStmKmGvQTe2MrQYDoR1Z4K4MqEwl+V4lScUf5SvhzZhwt3Z/Fd2ErDzc2Df9wx5u81knij8Z+DQQjLGh78M9LjLud+7EqPSGyAhRBsAf1i81BzAa1LKzy2WGQoK1l9UXloqpXzT0XqvagNkMgGHfgXWv84KUv9ZwMCnmFw+tw7Y8Aa5L/W7AiNfZ37HmTt8wgkmWU8sp+Hy9gdajWJ5uuUIz/Jw3IHMBP5eJ5fTKJkMQPWmQOebaVSKq3RlxgPbP6XhAYDAGsCk2cDOL4DI7cw/Tfziqs4PVXoDZAllDPNlAH2klJcsXh8K4Bkp5XXOruuqNUCJp4B/HiczuclA4LpP2aB55SzwdS8uU70ZMOI1JkiLSyznZwEnlgLhvzBc8fJlWNJxKsM4Z3MllR0q8fL4EhojgP1i3W9nAt9RN31aNLD2JeCUMqxl3EfkVa17hf1o/WYxH+dXpcx3o7xxtRmg0QBel1IOsHl9KK51A2QyArtnA5vepqcz+m2gywzmIrZ/Auz4nM2abSYA0+YVLz+RcoGJ1SOLgIJM5o563En+S5UaZbsfGbFAdiKQncxEcU4SL9TCPJICC/MoySG8AOHNnI2XFytlAdUsHqEsx4c0ZAjqLqGytCiGUocW0BMMqgP0vh/odY/j3+bKGWD1czRgdToAw15kOf/QAlbmbviBoetVhKvNAP0E4KCUcrbN60MBLAGHFMaCxsjhaJ6rygClXASWP8SqU9vr6NYH1SKBbuXTQOpFhg2j3y4+TIreB+z6ilUxLx/2PPW8i3wXdxL8ctOAhOMM666cISUg9SK9Ba2uduHFapVvAA2Ntx8ASYNlMjI8MuQCeRnQFDfz9mOVq0ZzeoS127ASV7ttyQ2qyQhEbAL2fsuckU8g0PUWoN8jJHFqQUr+tmtfBtKjGMq1HstEdVYCMPxVoP9jV42q41VjgJShg7EAOkgpE2zeCwFgklJmCSHGA/hCStlKYx2Wo5l7XLrk9NRYz8WhhcCqZ5nMHP8RDU32FZ7QJ5aSfDfhU2rdOMLFbcDm91gpCwglp6X3/fQgSgtDAcv7UbuBqL1M7lp2pgeEsjxdvSkfoU1IIqxSi4Y0qBYlOZwxgCYTeT956SRQZsRZEx2TI2jw1AZUgCGpyglq2Auo39l1XlDiKXqgRxfTGHaeDgx5lgZPCwU5wLYPWYEMrgOMeB04u5pJ/WaDgSnfkQtVyXE1GaBJAB6RUo52YtlIAD2llEn2lqn0HlBBDg3P4QVA00HA5Dlk/J5eBax4lKS8QU8DA59wXI2K2sOwLXI7WcMDHudd2T+45NtmMtHgnN/AcONyuNJGAYUE2I2J17oKJyi4bvl2x5tMNEZXztAYXg7nIzOO7/tVBZoNYuNt82FUUHR2+zITgF1fkv9jLKRHNOQ5+wnryweBv2dRU6nTNPbVbX6XpMtJs4G2E9yyyxWFq8kALQKwVkr5s8Z79QAkSCmlEKI3gL8ANJEOdqhSG6Ckc8DiO1jJGvwsMPQFlpDXvgQc/IVkv6k/ONausUxKB9UBBj3FsnBJCYEFOcC5tcCZNeyLyr4CQPCCajqQXKBGfTy7UpZ+mfSCC1tpOFOVompoE1IU2k+il+RMeGRbAev3CG8IWs26hgLm6bZ9RFLm0BeAfd+RxBlUB3jyeOWjNCi4KgyQEKIKgGgAzaWU6cprDwKAlPJbIcQsAA8BMADIBfCUlHKXo3VWWgN0ehWw9D6ekFPnsiIVuROYpwyKHfAEMOxl+0nmvAxg6wfMW5gMlJKY9jN1fVxFYR7Jh8eXMplamMNSc8sR5AO1HMHwqbIiNZI5ntOrlDJ8IfNI7SfTS3SmfJ4eA2x8Czi6iJ7eiNdZHNAyYtH7gSX38DODn6EHeTmchM7pv1Ngv5LhqjBAZYFKZ4CkJH9kwxsMYW5ewJ6po4tpkAC2TrQZp/15k4kXwfrX6Z10vw0Y/lrJ9HrijtLTOvonm0ar1KSAV8epQJMBVye5Li8dOKtwgs6tI4u8QXf+jh1vYNXNEWIOAGteIFM8rAcw8UttA5aXDvz7JMv9zYZQh2ndqyQvTl/IY1+JoBsgO6hUBsiQz5Py8EKgw1Rg8jesUK17Fdg7hxf9tHn2w5v448C/T/Dkb9gLGPcBLwJXUJBD9vTBX4DYQyQhtp/EsnyzIdeWIFdOivJbzDfrYQPA/VvZU2YPUvKGsfYlzj4b8ARDaNuwV0qW5lc+Ta9p6AvAlvd445j8DQ1eJYFugOyg0higvHTgfSWBOfRFEtayrwB/3skmyL4PA6Pe1K7aGAvZYLr1Q/YejXqLVTJXSryZCcC+ucCBH6kDXacDhbc638S78rUMKdnpv/51JvEBtqH0f4wJbHuJ65wUluGP/MYq5aSvmSOzxeVw4I/buPzI11khi9oNjHmXOaVKAN0A2UGlMEDZycCCKUxG9n0YGPseqyaLbqUxuP5LGgItxB8H/n6Yn+14I0v0rvBdkiNovI7+QUPWdgIZu437VoyeswopWVEryGalz2S0fl8Iluz9g8kXKq9tzU2lMNne78jnqdeZvXetx9rfhohNwD9PkJIw6BneXGw9yaxEFhyidtGwpV4kj6j/Y8DI/3k8X0g3QHbg8QYoI466PKmRwE3zgdajKe+w+HbyY6YvJF/FFlKyI3v9a/R6rvuMFRxnkXqJ/JTDv9Or6nqrY2KdO6GWx9UpFmlRyqDBeOU5kUZHGotfF0Dyol9VJsJDGigCZfUorlarJUXKQsLca6QM+QyzdnxKNnlYT2DEq/SItJCXQc7Wkd8YHk+dW5Q7ZCgAVj4FHJpPfpGPP0PhztPpPXlw+KsbIDvwaAOUegn49Xp2Vc9YRE7KoQXAisconnXrX9oVkZwUYPnDJLS1GQ9cPxsIquncd2YmAFvfZ15DeJGIOPDJsqu8FOaRBR17kJyhuCNA0nmymVV4+ylGQxnZE1yXvWd+QTQsfkEKI9oCJoMyJ0wZapifybYOVZwsI45tHCp8q5Dn00CZ3BHWHajdrvQXtbGQjbtbP2S7RrPBwNj3efy0cHwJ8M+TlDKZ/DXza5aQkuva8i65SfU7syjRfhJww48eO25aN0B24LEG6MpZ4NdJZOrOXMpk8a+TqMfcfCi9Ia0G0Kg9wF93Mz806i3qzjhzZzfkU4h+28f8u/ttDAeqhbl3vwqyuY2R24HIHUxkq0JlVWqyM792W/PgwJotaXjcHUJJyRApSREnSz5PPlXsIebbABqlxv3YbNpiOGfKl3Q7CvPIBdr6Ab2dXveyB0wrf5YWzdze5QM0/sNfLVpRVG9E9ToCrUaTO9RuInDDTx45Xlo3QHbgkQYo/jjw7QA2WD6wjSfZ5vfomXj5AC/FaZ9k+39kk2O1RuT0OFOqVRUS173MMK/NePaKuTPUSo6gSuCZ1ezONxm4H2E9gCb96XE06EZlxIrMKwH8PVIuMMemkhGTzvC94LpmIfqe95Ssaz0nBdj8DvNEAaEsHHSbWXS/DfkMycJ/pqdz409Fc3dn1wF/zKTn1mYcjZCzTcblDN0A2YHHGaD0y8DXvRk6PLKPzZJbPqDL3XUmcP1XRROOxkKerAd+JPHvhh+cU9pLjwFWPsNQrXY7YOy7vNO7A/HHgeN/0bglneVrdTuRlNhsEJUGS9PmUZ5Iv8xBixGbuU8ANarbjmcpvMVw1xnK8cfYQhO1m5+f+IV2m8bBX3mMQhow5LYVt4/YDPw+A6jehA3I2z8mF2vaPI/iYekGyA48ygDlpgI/jaNhuGsV4/utHwGb3+bol0lfFzU+aVHA553494DHybAt7sQzmWisNvyPidxhLwN9Hix9viP9Mi/QI3+QF+Plw/aLNuN5h74aRtSYjKQ9HF/CcnhuKl+vUgu4/W/XRMXU47D+dXpAo/4H9Li76DGO3g/8fjO9sxmLgMZ9rN+P3AEsvIkNw20nMCfU407gus8r3qNUoBsgO/AYA1SYC8yfQs7HzCVMVm7/BNj4Jqsck78paljSY4DPlGTm9V9RFKs4pEYCyx7knbf5MGDi5+w8LylMJvZ87fuezGBIVnE630yypLPJ78oIYyFbMxbeaH6tUV/md9pPcj4MSr0E/PMY19VqNBuKbVtXUi4AC27kMZ86l7KwlojaA8yfymPZuA9DvEHPsPLmAdANkB14hAEyGoA/72C4cuNPbGXY8Tmw4XUKmk/5tqjxSToH/DqZkhMzfqenURyO/EFWrRBkQXeZUfI7ZG4qK2UHfqRRC6pDYmKXGeVTqvc05KSw2nXgRxqLkDCg/6O8KTjTXycljfi6l5mMn/o9Q1VLZCcDi2ZQq+n6L4vecCI2A7/dxMbf6s2AY4uBsR8AfR90336WELoBsoMKN0BSsj0ifB4w7kNWrQ7O56z1jjfybmdrfGIPAQtuYJl85hKecI6Ql07Dc+xPVnWmzi15OJSVCOz+mgnvgkyur9e9zDt4WOKzQqB6hDs+By7toDHp+xCHPTqT84o7Cvx1FxP3Q19ki4ZlSFaYSwJqxEZg/MdA7/usP39yBW9mTQcxL3V+AzDjD/LHKhC6AbKDCjdAao5n4JPAyDd4wiy8iSHYrX8W5XVc3MakY5UawG3Li/c2Ek4Ci25hrmjoCxSlL0muJz2GF9Wh+eZBewOfKN74XcuI2kMpjnNr6SEOfxnodlvxObr8LBIOj/7B0vrkb62NlyEf+PMu4MxKViz7P2r9+UMLgL8fYRiYcpGPe9cDddq5fx+dhG6A7KBCDdCZNUwudr6Zynfxx4Cfx9F9vmtVUZ5PxGbgt5tpdGYuLV6l8PgSilz5VwVu+lW7z6g45KSQzbt3LslxXWewefJaDLNKipgD7PmK3kMu0ai3WA10FP5KSV7Wulf4mem/sdKlwlhI9YMTy4Ax7wH9Hrb+/Jb32bja6z6K4PsEAPdtqjBZFN0A2UGFGaDkCGDuMKBGU+DutQyT5g5lWHXvhqIynJE7GXbVbAHc8Y/jfi6TkfmjXV8xKXrTL2QSu4LCXGoFbf+MeaautzAkCG3k6p7qAMwa0BteZ44IcG5C6vkNJJZ6+QK3LrZWLzAaGK6dWsEKabeZ1t+35F5WJQc+ybachr1YqauA8rxugOygQgxQfhbw4yi2BTywFQiuB/xyHT2ge9YXLeVG72c/WEgYcOdKx9o9BTkUszqzine/Me+6nps5s5qExrQoNlGOeM1+60BFwViotFTE8jkvnQzjvHSlOdVgnhsvTWzV8K2itG8EUbcnuB5bTILrkWBYHr1UhgIajdP/cnucoUAknQcWTGVLzvQF1lwtQz7w+3RW0Kb9Qt0gFYV5wC8TeV71uJOSLYOfBYa/UlZ7Zxe6AbKDCjFAS+7jnWnmEp5MK5+mdvCNP7MCZomEE+QGVakB3LXacdiVdYUhXewhJrRtE5TFITUSWP0CiYm12rBzvjgR+7JGZjxbJK6cAa6c5nPKRbZRaE29UBtPvX2VcT3Kw1hApcbCHO3v8fJh6FurFds/archO7t227LxGNJjeNzPruH3TP7WsXRuZjzL8FdOsyrayaL0X5CtUDgOArcts66eZSUC3w3hTaheJ3phty4BWo10/z45gG6A7KDcDZCqXDj0JWDo8+aEYf9HmVC0RPpl4IeRACRwzzrHlSv1LpmVCNz4o2si5iYjXfRNb/NiHfo8qzblXdUy5JMHFbOfeZPL4eyKVxFYnYztms054yukAXvVqtZnW0NASPFTM0wmGqG8NDbeqh326dHsB0s6D6RE0GABZDzX78JpGc2HAI37u29woJTM4ax6hp7r2Hepx21v+/PSWYC4tLMo7ys3DfhpDPfl3o3WjOmovZTqbTaY+5wZBzy4na0v5YSrwgApUy4yARgBGGx3SAghAHwBYDyAHAB3SikPOlpnuRqg1EhgzkCGM3cqLQrfD6NI+8yl1m54XgYT0qmXgLtX8+5lD4mngF+uZ7hxy2JeLM4iOYIGMGo30HocMOET9zeg2oOU9G4iNrPV4dIus4cS2sQ8IqdeJ3oiQbXKh9lrMpr7wS6Hs1M/9jC1oL39eLxajSL1oLgcjjPIjCc59MJmVr0mfmk/x1eYB/xxK2e9Tfoa6Har+b3USOD7EQwv791gvY7d3wBrXwS638HiRINuwO0ryk1H6GoyQHbH7CizwB4FDVAfcC5YH61lVZSbATIaaFCunAYe3EHZ1O+Hs2v9oV3WMqpGA/DbNDZC3vonKyb2EHeU+SEvX+COFQwdnIGUDPvWv8bPjv+Q1biyvsBNJjZ5nlzBcCA9iq/XbKWMvxmqTM0ogUZ1WaIgm0b6whYgYguQcIyv1+vMcnenadZVKldhMnGe2MY36dnN+N1+3q0wz5z3mTyHlUkVUXuZ92nYC7h9uZnGISU1pE6vpNE6+CtlQfo+VPJtdgHXigH6DsAWKeXvyv9nAAyVUsbZW2e5GSC1reKGHxm/r3qOI1e04vE1LwF7vi6+veLyQcb+fsE0Ps6WxnNTWZ4//S+naVz/VdkPv4s7SqbwiWUMe7z92ArSdgLzYJWtupZ6iQb05N80qBA0nt1v5z6VdHxO9H52uOdnMtdjmVS2RGEuWc8Xt7PKaakbpIb5fR9hWKciLx34pj+3LaQBQ90Htjl/0yoFXDFAkFJ65APARQAHAYQDuF/j/X8BDLT4fyNosOyus0ePHrLMceWslG/WlvKP2/j/mbVSvh4i5eoXii57ZDHfW/Wc43XGH5fyvcZSftZRypRI57cl5gA/878aUu6aLaXJ5PxnXUV2spS750g5ZwD36c1aUi66Vcqjf0qZm15231veSI2ScvN7Un7Snvv5QXMpN7/P/S8J0mOlnDuc69rygf1jlJ8t5fcjeG5F7rJ+b+Wz/PyJ5davR2zh64tmSvl+Uym/GyKlobBk2+kCAByQzl7nzi5Y3g8ADZTnOgCOABhs8/5KDQPUQ2M99wM4AOBA48aN3fUba8NolPLHMTQWGfFS5qRK+VFrKb/uJ2VBrvWycUelfKuulD+Nk9JQYH+dyRFSftRKyo/buGZ89v8o5f9qSvlpBymj9pVod5xCwkkp/35Uyrfq8GT/drCUe+eW/IKsLDAapDy7XsoF07jfb9fjjSQt2vV1FeZJufQBruefJ7luLWQlSflld55fiactPp9PI/ZOmJRJ560/s/oFrnf5I3ze+ZXr2+ciXDFAHhuCWUII8QaALCnlxxaveV4ItukdaiurRLEVj7Lydd8ma8GwvAzgu0HkiTyw1f5onYxYVjvys1iWd1S6VWE0MAG5by71gqbOdU2Y3hlIyf6kXbOZTPUJYE6p932OE+hXKxJOkgh6bDGF5Xrfx4morvzuUnL+287PmfCe+r32xNrUSOCHUQyt7t9iZjunRQPfDWZR4d5N5qpmYS5L8wVZ7J6PPQw8srdMw2BXQjCPlNcXQgQJIaqqfwMYDeC4zWIrANwuiL4A0h0ZnzJHdhKND0BR9wtbmPzr/2hRtcJVz5L4N+1n+8YnP5NckJxUcoicMT65aUxo75vLaRa3/OFe4yMlRfJ/GEGm9pXTJC4+eZId29ei8QGAuu2BKXOAxw5RtGz318AXXTlhxFDg3DpUjaAx75HtvGgGE9C2qN6ULOnsKyQ5GhWJ29BGvPHFH6OapgrfQOb9Mi7zXJAmCtp5Cpx1lcrzAaA5GHYdAXACwMvK6w8CeFD5WwD4GkAEgGMoJv8jyzoH9O9TdHHjjjFe/7yzlF90lbIgx3o5Ne+z+T376zIa6Nq/UV3K8xud+/60aCm/6sWwK/zXku+HPZzbIOV3Q7ntn3aU8sDPdP11FEX8CSkX3sTf6qteUl7c7trnD86X8vVqUs6fyvBMc5kFXP+6V61fX/6wlG+ESnlpj/Xrfz/K8+nPu/m5Uytd2yYXgKshB1QWjzIzQImneXBXPsP/N73Lg3xhq/VyKRelfLehlD+MdpwMXP0iP7/vBye//wyTou82dP1kLw4JJ3kh/Gd45umGx1mcWcMiwOshUi57iDlBZ3FgHj+3YJp9I/TPE0WTz7np/M7Pu0iZl2l+PTtZyg9bSPntICm/6smbYxkdR1cMkEeGYJUO615leXzIC4zFd34OdJhCNqoKKVkOB5iXsdcPFP4Ly/J9HgR63VP8d8ceAn4ey7Ezd650TqzMGWQncTT0nP4sF49+B3j0AIXIdC0g59B6DPDwXsqiHFkEfDuQkqrOoMcdlFk9t5bkUZOp6DJj3yd5c/HtlAMByBCf/C2HGW55z7xslRrsFYw7QqmOlAtUUqxg6AaotIjYzJNk8DOUJF3/KgBBGQZLHPyV42lGvWmfxBZ7mPmhFsN5wReH6H3AvIlsIbh7rfbQQlchJTk8s3vSGPa6j7mN/rNKzne5luFXhSOW71lHouC865hsVnM3jtDzLubYjv0JbHqz6Ps+/pReAVisUPNNTQeQBb1nDnNCKjpNo8GK2gs07M1cUW5aafewVNANUGlgMlLDJbQJ1Q0jd5J8N/AJ6ypDRhy9pCYDeWJoIS+d6nZVagJTfyi+W/tyOBPBwbWBe9a6R7MnOYLDEZc/xAmiD+0ka/pq1nouLzTsCTywnTPYdnymdLwnF/+5gU9xYOSOz8hmt0W1MBJeAevk88g32E/3zxNm70kIYMw7JIcG16Hx2f4xKhK6ASoNVj4FJBxn9cLbj3e2qg04w9sSq59liHT9l9r9OFLSzU6P4YiV4i74uCNkRQdWp15QaZnNJhPnm8/pD8Qe4Wjnu9ZUqKreVQn/YFakJn3NkGnuUB5LRxACGPcRpVJWPacdwnW6kbSPHZ/RKwaUkOsdDjw8OM+8bOO+LPNHbOYQxn3fsz+tgqAboJLCZKS2M0DJ0nPrSdMf8px1B3XEJtL4Bz9r30sJ/5nLjHi96BgWW6gC9f4hND6l7XLOjOekh9XPAc2GALP28Y5bTo2L1yS6zWTTsTQCP4/nOeII3j7kBdVozimqGbFFlxnzHlUDlj1oLt93vhloMoD8tLwM87Ij3+ANUXhRa2nnF+7aM5ehn2UlxakVfL7pV3owm94iR8NSqc5oANa8SO0ZWy1fFSkXgLWvsFfK3jIqshIZdnl5U+2uNA2RAHB6FfBNP3amT/iEvCFX1RR1lAxhPSilEdqEuuDH/nK8fEAIcPMCSnksvr0ovyggBJj0FeVFdn3J14Sg7EtOEgsjKmq2IFft4jbqCR34idIdFYBykIW7CiEl3d0aLTih8vQ/QPxRaj1bCssf+IlkvZsXaidwTUZg+SMUyJo023F3ekE2GxKzrwB3/lu6nI/JyJHB2z9hh/cNPwK1W5d8feUNo4FVnitn+JwZT92bzHhW7wpzyPwtyFFGQ3uToezlwzAosDofVWrQAFRvyketVvy/vAb8hdSnHvjvM6hsWZDN6pc91GkLTP6aXtDmdxj6W6L5UFZft38CdL6J+xTWncnn3V/Ts1U95kFPAYcXUh3BWECjNcaJwoeboRugkuDCFsbuE7+gG7vtI0pMdJpmXiY3lSOWmw2xLxi2Zw4QtYtlU0ehlMkEvKvkeWb8Ya0V7CpyUniyR2xiN/e4j7Qp/54CQz5HP18+QOGyuCP0Gk2F5mV8AnkxV63Pbm+/YEWOtQqNjsnIcMdkJMM8N5UJ2CtngXMbAEOueV0BoRQla9AVaDoYaNLPuVlfJUVgKHDbUuCP24B/HucNrOst9pfvMIXHbucXLPM36W/9/uh3OEd+zUvAjN/42vBX2cm/+V0OvQRonLpMB47+SbrIgZ+YJnBmzLcboRugkmDn50BwXU4xvbCFpc7rv7KW89z1FU/yMe9o31EzYnlCtB7HE8ERtn7A59bjgDZjS77diafZqpEZT+PZ486Sr6usYDLSyERsYqI0Zp9ZsbBqfba1tB1PGdlaramYGBBacq9FSkq+plyktxp3mN+/Zw4vci9f6u20Hc9cX1n0UPkGMrz6/WYWI7z9rGVYbTHmXYZPyx4AHtxpPVGlWhgw5FkWRM6tp5ha9Sac57b3OxoZVVht0NPA4d/ZnlGYQ/qF7cSNMkalaEZ1F9zSjHo5nOJiI9/gBIL5U6jl/MQxc5iVnQR83pl3qGk/a6/nr3uYeJ61z/G45LNrGXp1vZXVk5JeaJE7ODPMJ4BjX1xRUixrGPJpbE4u5/7mpvD1ep15d27Um/yV8lJvBBi+Re9RRMk2mfk0DXvxhtH5Zo5Acvd3LryRlazblzsmlUbvI/enx52sWlrCUAC8rYi8vZpEryojDviiC9DlZt4sVfx5J3B+Ew1rYQ4wK7zUBYhK34zq0VisxOg972ZIELGJHCDLHM/Oz+nWD31Rex2ROylUP+Bxx8Yn5QLFpup1YpK4pMbn2F80lMH1OInDE4yPycTt+l8N4KOWvPufWcV56VN/AJ45Ty3jMe9QgKs8jQ/A8K3FcBJHH9xBMuaI19ldvvJp4JN2wMpnyJ1y53dOX0gPZdGt1Ky2h0a9yZY/8DMQE279no8fb5AA8zwAQ9Tut9HjSYs2L9v3YSA/ncY05QJVDsoRugfkCvLSgfcbk0V673rgzVrMRTwfyaQmwErV54p059Tviq7DaADmDuG6HtlnX/S8MI8i9enRlOxwZKgcYd/3FEJvMpBjXtTtrChkxJJhfXgh9w1gKNvpRubLPL3NQ0reePb/AJxYyiR355sd0yxcRcpFKg4EhAL3b6busxbyM4HZvThm6L7N1uRVKTkOKiMWePQg83xp0cCXXXnzHP+RebnvhwM5yTSuDbpSGrgU0D2gssKJZXwe9z41erx9WQWzvKj3zQUMeTwhtRD+M8mLo992PHFh01vUIp46t+TGZ/fXND5tJjDRWZHG53I4h+d93ok5rVqtgBt/Al5OoKFuNcrzjQ9AL7RRL27zkyfoQZxYztaVlc8wwV1a1GjGymlqJLDiMRoJLfhXBca+xwqsbV+XEJwJlnEZOPgLXwttxELJoYW8AarL9X0YSLvESSzn1rvXqysGugFyBYd/Y/KzQXfmKwpzrFnPBdm8M7adYD0qRUVuGsfhNBtsretri8gdStn0HuaRSoIdnwFrX2Li9KZfKqaPS0rmud6oxrvsmTVA7weAxw9zplXHGzy7AlccguswRHziKI/VgR/p/a58xr7RcBZN+rEP7ORy7RYMFe0n83za+r412RCgR9l0EMvyhny+1ucBoDCbRui/dUxiC5CXNwBJcmQ5QTdAziI5Aojey6kEQlDpsGYrxuIqDv/GO6A9QuGebzijarSdyhjAz8+bwLvg6Le0lykOe+awCtJpGjk+ltyk8oBUhMu+H07RdYC/yVMnKZxeUo/OUxFcB5jwMXD/Vo623v89CaNZV0q33v6PMSe29iVWMLUgBPM9OcmctGH73qCnWOU7poRVDbpxEsm+ueYeMR8/hpGXlalWWfHONcu6AboBchZHfifnp/PNTA5G7SbrWTUkJhMNTFhPHmBb5KTQMLS73nHX+ldK6Dzlu5LxT478Aax5gTOnpnxXPiOILRF3lA2tC28kA3fiF8ArVxhyWpaLr0bU7wy8lkLeTeQOYE4/VvdKCi8vYNI35DWtmEWKghbCetAT2jWbOUhLNB8G1O1Ij1r1yvo8QALn+fXm5breynxmU2XS6sUtJd9uF6AbIGdgMlHPpfkwNn4emk9mrSV/5+IWVhH6PqTt3ez+mnfHoS/Y/57L4bxoO0yx9qycxdl1wN8P8ySa+kPZjBm2h4w4YNlD1CWOP85x0bPCWSauDLkdd8HLm9Is928GqtSiJ7Tv+5KvL7g2MO4DjtXZ+6395Ya/ysqrlhfUb5YyFFKpcLW7ntt2aIF5uXodScDMSgT8q5GgWA7QDZAziNzOik3XW+iaHvmdrrFl39TBX5nkbTex6Ofz0kkCaz/Z/gA6k5Hl3eB6nJbpKmIPsUeobgfyfMort2Iy8QL7ujepBf0fZcm6zwPXluGxRd0OrJS2Gs1CwKrntEXFnEGnaeyG3/Q2Db0WarVkLufAz+YEs4qON/C82qtUZb19uc6za+iZq+hyC5B0hioIp/8lL6mM4ZEGSAjRSAixWQhxSghxQgjxuMYyQ4UQ6UKIw8rjtTLboBNL6Qa3ncDWiawEa+8nOwk49S/QZYZ2svfgfKAgk8RFewj/mUZkzDuuhyqZCeSNBNUCbv2r/EKdhJMkw616hrmFh/cwb1XOdH6PhX9V8nr6PsLBlP8+bj+McgQhWO0yGWiE7GHAE/SyD9iQX338mC44v8HcSd91BhnmJ5aal1MHI/r4s5fOMkQrI3ikAQJgAPC0lLIdgL4AHhFCtNdYbruUsqvy0JCMcwOkJDu3xXBS5k+vJJu41SjzMkcWMX7WmmxqNPDO02QAORZayEnhJNWmg3i3cgWGfGDxbVzH9N/sT9lwJ0wm5hvmDgGSzzPXdPvf7uPBXE3w8uZNZfCz9JL/cVBWd4QazelVHl5oX0OoQVemCfbMocyGJbrNZMuFSkys1xmo04HERBUhDcj0zkkmB+nMGte300V4pAGSUsZJKQ8qf2cCOAWgnKmwCuIOs9O6zTieOKdX0hhZJoiPKg2iWgJep5WZ6H0d9Njs/Jwl1HEfuM52Xv08q3OTv3GPJGtxyIgDFkwB1r3MuWOz9tMbLK8O8soIlZMz+FnmXba8X/xntDD4WYb5mxx0rfd9iFWsszbGo0Yz3uAOLeANRAig8zQ2+abHmJdrdz15arXbUGq4JB6bC/BIA2QJIURTAN0A7NV4u58Q4ogQYrUQwk5ypZQ4uxaAYCwfd4S5IMvu9uQIEsHseS77fqDEQ5tx2u9nxNFD6nyT/fyQPZxYxtBtwONAx6mufbYkiNikiNTvY3Vr+kLzYDwdxWPYy0DXmeTsHFnk+ucDqvFGdm4t+w+10HIkEBJmFsuzRLfbSG6MVi6ltkq+8vRK8zJqDtNkpCcUU4aDPOHhBkgIEQxgCYAnpJQ2LCscBNBEStkFwFcAlttZx/1CiANCiANXrpSAl3F2Dd3SoFrKXUWwK13FSeVrtYiFKReASzsYmtmrSG37kLG9vb4xe0iLAlY8Ts9r+KuufdZVSMnO8AU3MPH+wDZWt3SvxzUIAUz8nG0x/z5JORBX0fte5iN3fK79vpe3ku/ZyHPEEm3GAd7+ZjG9Wi1JrD39r3mZGs2oMiC8KGVydrXr2+gCPNYACSF8QeOzUEq51PZ9KWWGlDJL+XsVAF8hRJHbsZRyrpSyp5SyZ+3atV3biNw0TqpoMZz/X9jCONtSs/nEMvaGaen5HP6NB7LLDO31p19mgrr7HWaJBGdgMgJL7mNMX9ZEw4Ic6getf413x3vWs41CR8ng7Qvc8APziV/3Yj+XKwisTuN/bDGLAFroNhOAZGrAEgEhPJdPrjDnodpOYHO0ZTWs+TB2/9fvAlzY6tr2uQiPNEBCCAHgRwCnpJSf2lmmnrIchBC9wX1xYsyAC4jaDUBSFiE/i1yMZkPM76de4oHS8n5MRhqgFsPtd3LvnUMjMqBIkc8x9n1PqYgJH7tmuFxFTgowfzJwfCk7waf9QkVBHaVDSH1qOAPsj3MVaj5x2QPa74cqDdMn/y76XvtJQEaMmfXcZjzF2ix1qVsMJ6fIy5c5UNsWDzfCIw0QgAEAbgMw3KLMPl4I8aAQ4kFlmRsBHBdCHAHwJYDp0t2t/ZE76LI27EVjZDIAzS0M0PkNfNbq17q0i42A9tTt8tKBA/NIOnRF2zktmhWzliPJyi4rpMcAP48jNWDaPFL69ZDLfehyMz3K8xs5aMAVVAtjTjIrwX7LRIfJvDnaNpa2GcfQ6oyS92nQjQMOLm4zL9N0oCLVms8bpDr0sAzgkQZISrlDSimklJ0tyuyrpJTfSim/VZaZLaXsIKXsIqXsK6Xc5fYNidxB4+MbAFzcSqW6Rn3N75/fCFRrDNTUaDw9s5rLt7LTTBo+j9ygAY9pv68FKcm5gSydPlBxSLkA/DianJGZS3ky63A/JnxKRYQ1Dtjx9tDjThog9SZoi3YKp8fWCwoMZbuQ2iLi7UODc9Ei1PIPZm7RaOA5HLnd9e1zEh5pgDwCBdmsbjXpx/+j9wP1u5olNAwFvGu0HFHUEEhJca1mQ7RDFpOJ8glNBjLOdhan/2UifNhLZdfQmRrJaauFuRz13GxQ2XyPDnK2Bj5FI6LO83IWrUazncI2z6MitBG5PloGqsUwerZq3qfZEB731EvmZRr1ApLOUgzvkvvv7Sp0A2QP8cfpfjboTlJX3BFrJcHLB+jBtBxR9LNJZ9nsZ6/0HrmNB9wVTWZjIRPBtdsCfR5yZU+cR1o08MtEsmBv/7t8eEXXOnrfR0NiOcfdGXj7Uqf6/IaiI3pUtBhGw5afZf1686EApNmzaTaYz5d2mpdp2IshmJcveUFl1B2vGyB7iDvM5wZdgcRTTMpZTqNQ4+LG/W0/Se8HYP+OFsJ/sd83Zg/h8xgajXqzbDrcc1Io25qbTj1i3fiUD/yCKAQfscn1snybCWy9sBciNR9Khr6tBxPWA/Cras771G7L/y9bSLs27MVnUyEF9q7YkQMpJXQDZA9xRyh1WbW++cBYGqDofdQD0hqjfGY1wzWt6ldOCkOpztOdbxjNyyB7tukgut7uhiGfuj1pl4BbFjExWUEwmiSiknOw41wS/jkSi0/WncHb/57EptMJOBKdhtyCsmXmVgi63U5PI9zOAAN7aD4E8K1SlPWsonE/FlEu2pTSvX2BsG5mkqGXF2+0amUMYFtGSJiZCa3ekN0MfSyPPcQephERgrmggGrmvIuUZJO20VCOy06mcbInu3FmNZsAO9/k/LYc+JEyHaP+5/7Es5TA37Poft/wY9E5U+WA3AIjFu69hD0XUrDhlPaEzh92XPzv76oBPnhgcHNM6d4QYaGB5bWZZYfg2kC766iyMOot5z1c30BqT9nL0fgGMocTe6joe2E9ODqqMI83wgbdKPdhKDCrGNTtSDKjX1Wuw3Lqr5uge0BaMBmZx1F7u66cAWq3M1/8yREcHaOl2RO9B4C05gtZ4uRyVs6c9TIKc4Hd35CbUZqBhPawby5JbcNecTyLqgyQml2A91adQrvX1uDtlafsGh9bZOYZ8PG6sxj4wSbc+8sBHI1JK9sNLQ90mEo1zKjdrn2ucT+2ZeSmab/foCtF4mylQMJ6klYSf1T5vztvjAnHzcvUacdm49ptmIYoA+gGSAtpUYx9VcbvlTPWo4vVg6bV3R61h6VLLQOTm8byZ/vrnfdkDi8EshNZLXE3YsKBtS8zVzXoafev3w5MJon5ey6h21vr8d22CyVej5TAhlMJuH72Tjw4PxyJGXlu3MpyRovhDJcs2yKcQZN+AKS5v8sW9buwWJJ60fr1sO58VsOuuh35fOWMeZk67XgdCC8aojKAboC0oJK3arZkSJWTxJ4ZFYmneFAsX1MRtYehm1Z+5/wGHlBHgvSWkJLeT1hPx0PqSoK8DA6lq1ofmDyn1MPonEV6TiGmz92DV5cfL35hF7DmRDx6v7sRa47bEezydPgHA00HWBMCnUEDxZCoN0Vb1O/KZ9swrGp9phWSlMR39aYkKCZZJMLVCMBUSM5RGTCidQOkBdXa12xpPiC1LQ3QSaBGi6JGxpDPZF1jDU1ogN5PQKjzodSlnUBKBEu17s79rHuZlPwbfwSq1HDvuu0gOiUHPd5ej32RKXaX8ffxQv8WNdGjSXU0qxWEjmEhaFE7CD5ejve/qj/zJg8uOIivNp5DpZx316gvb262ioaO4B/M1gt7IVLttnxe/Zz160LwBqqe396+1ByyNEA1bXr+ysAL0g2QFlIimHgLqm3uKA61aJdIPAnU1dBHS45gHF1Pg1woJXBhMysXzmo1H/yV+rwqq9VdiNjEdfd/rGTa0yVAbFouBn24GQaTtmHo0rAapnQLg5+PF3ZFJCP8UiouJmXj+OUMRFzJtvs5FZn55Kn0blYDn6w/ixeXHqt8RqhRLwDSuhrlDGq3sz81Q00oV21Q9L1ara0NTq3W1m0hflXIUfJRbrRlYID0KpgWMmJZQheC/VyAuaRuKCCJsKNGwlY9QFozwZLOcV3Nn3FuG3JTSaPvNtPxAENXYSigPnGN5q5LgJQQeYVG9H9/k+Z7jWoEonfTmjgUnYplh/hbdwqrhsY1qqCKnzf8fb2QlWfA2YQsCAHEpOYiwNcLCRn5mus7HJ2GQa1qYdH+aHh7CbwzpVOZ7ZfbUUfRg0o6RxKhs6jdmkoNUmp7yi1HMY9oi1qtgMML6HGpVd7zNqOZQxuZO+fVa8GN0A2QFjLjzYLzGZcZNqkKiBmXyZAObVz0c6oBqqEhTaqyTO1Vx2xxYhkJYN1uc2nTi8Xeb4Hkc8Atf5abcP3Tf2pLiHZtFIpCowlLDsagdd1gvDy+HZKy87HmeDxWHmMux9tLICTAB/kGE3IUDlB6rv3vKjCYkJxVgAbVArBwbxRGd6iHIa1dlGGpKATXAXyDSDh1BSFhZC3npmqH06GNydwv8nojPqdfpgGqUpOE24Ic802vWiOSEP1DeF24GboB0kJWgrkClhHLA6xCnWduzwAF19MWhY89RENWo7lz23B6JZd1pVesOOSkAFs/ZINs6zIgNGrgYFQqVh4tmhge3Lo2tp2lQNwzo1vjTEIW3lnFPMagVrVwV/+m6NSwGjo3DIWPl0BOgRFeQiAjrxA7ziXhi43nEJWiPbXhZFwGnhzZGquOxeHZP49g/VNDUC2wnIczlgRCUBkh7VLxy1pCvVlmxtk3QLmp1B7yr2rxOSUsy4xlSqGKQqrNSTYbIHVcc2gjs6C9G6EbIFtISUsfXJf/W3pDgEVOqFHRzyaf1+6MB4DYgyzNO5NMzs9iNaSXm5PPe79lSXbk6+5bpwNIKTH1m6Ikue6NQ3EkOg3eXgLPjG6DBXsu4XJaLsZ3qodnx7RF4xpV8N22CPy0M7KI8fLz9kLz2kEY0a4O+jSriVXH4rDiSNEL47MNZ/HL3b1xx0/78PTiw/jhjl5ltp9uRWAN+5wee6han8+Zcdqyvqpsbk6KjQFSDVe8zXJJ5vM7qDa9ooBqugdULshLY9lRnS6Rl24t+pWuxMEhGgqISefMo00sUZBNbZYOU5zbhgtbmMxuY6eXrCTITQP2fMv+M1e1p0uIfReLVrvCQgMRnZoLfx8vjOvYEB+sYfJ06cP90SmsGl5dfhyL9kf/t3zNID8UGEx4eBgNe2pOAU7HZ2Lh3ij8vDMSg1vXxmPDW+LLTechYEJ7EYV2XpfQVMQjeMNS/FgrH8fP+eLKntGo3X4oxcA8GQEhRaVUi/1MNT7bU1cMrM7n3FRr7SnVcKmzxiw9INvPevtZi9e7CboBsoXaOazeKdQEnYqcZP5vO3QvJ4XsaC0PSOUV1XOywfPsGla/Gvdzbdsd4dACID+dkxXKCTfPLSpk5e0lkJ5ViJt6NsSCPVEI9PXG9ueHISEjD+O/2I5zifz9v7m1O4a1qYNAP+2KYWZeIX7fF4WvNp5HXlQ4fq+3Ey1Tt6G2YAm7UHojKyEQIYG+GOadBq81S4E1AoAE7lxFzo0nwjcQKHRxIKCv0o5SaCc5FhDK59xUm88FUF9afb2KhaekQg3pTAbX5WOdgG6AbFGgGCC/YIZjtgYoN9V8V7CEyjTVyvGoMb2zlY3IHdThcZfWs5Tspm/Ux705JQcwGItOAQ2t4ouolBxM60Hj07B6INY/OQSrjsXhxaXHUD3IF3Nv64FR7etCaISealldCIGqAb64v0U6bj/zOQIu70JOqj8OBfbBuxkdcFi2xCVZFyZ44Yebe2LNsWhEn9yHhYNS4LPjI2DeeGDQMxyV42kqjyYjCYGuwFfJ19gzXOr5q8Uv8gtmWA6YjU12kvl99Vw3FtIA2au0lRC6AbKFpQdUmMtwzN8iqWzPAKl3jSCNikuKYpycERHLTqIxc0UrqDhc2sXK16A57ltnMTii0Z+VllOIKn7eOJOQiepVfPHng/3w9sqTWLg3Cv1b1MTsW7qjRpC1Z3k4Og1LD8ZgV0QyolNyYDBJ1PQz4ONqf2Jw+goEVKmFhH6vYtTmxgirUR+n0qzZuksOxmBm36a49VAC1tWdgvEvPk4Fwu0fkz18a/nMQHcaJgM7412BarBshxGqUG9kJg1NH/9g8zkfEAoIb+sQTDVe0qRIc+S7tXrqsUREIcRYIcQZIcR5IUSR1nJBfKm8f1QI0d0tX5yuxN9+QSxtAmYXFyjeAAVqVCFSL/Lgan3OFjH7+exOguCJZYBPIGfTlxMOR2uzeb2EwNGYdDw7pi3i0vOwcC9/75/v6mVlfM4lZOKW7/dg8tc78eeBGDSsHog7+jfFC3388KfXSxicvgI/GsZh+/h1qDvmGTw8rgdOxWUUMWCrj8ejT7Ma8PfxwsFLqbyxXD+bofK5dWU+9cFlFGS7zvtSDY89j1koYazUmE3vF2T2+r286AXlWHhAKglRfVaXdRM80gAJIbwBfA1gHID2AGZojGYeB6CV8rgfgHtu76o7aywwa6EIizxEYY55GUvkKgZIqwyactH56RVxRwAI94VKqjxsyxHuJTQWA63KFABkKYzlG3qE4Zk/j6Bh9UAceW00/H3Mv/GWM4kY9dk27IpIxisT2mHfyyMw767eeKlrHu47+yCa+GchZepi/FHzIdy24DQy8wpxQ3cWBVKytdUB29arihOxinckBPCgwsv61c0s89IiK8FcgXUW6o3S21/7fZV5rzXl1CeQfDMVfsE0gipUo6auw9X8VDEo1gAJIWYJIZy4dbsVvQGcl1JekFIWAFgEwLaDcxKAXyWxB0CoEKL0JY4sRRLCJ9B8wCxbJ4yF2ncatXRqmS9SkRlnzSVyhKSzLIFajn4uDRKOkzyppV1UhjgSnWb3vXb1Q/D3oVhcuJKN1yd2QLUq5t9zf2QK7v81HO3rh2D7c8Nw76DmMGYlI33fb8C863hc7l6HGp3H4Nkx7HPacS4Jtav6o15IANrWq1rk+7LzjWjfoBpOxlmEZ74B5mm2mc7JgJQLMuPNFVhnoUqyevs5Xs4uLHI6QljPrlfXqYZvbm5vccYDqgdgvxBisRIWlUfWLgxAtMX/MSg6G96ZZVzHP8qMLpPB/KNbGiBToXaMbixgLK7V55WbxmkEziDpXNEmwNLgkqIv40Hi8m9P7ojnlrB7e2Q764vtnZWnUD3IF7/d1weNalTBzp1bEDq7DaqtegiyWkPgnrX/SaO0q09jc0oxLIF+3vD3Lfr75xmMqBcSgPTcQuvkeB9lwlOsi71XZYW8DHrSWkMuHUH1vu2F+KrXYplK+A+2BsXm8lYNkBrmaYVxpUCxSWgp5StCiFcBjAZwF4DZQojFAH6UUkY4/nSJoWXkivmlNJeBEOJ+MERD48Ya7GV7MBksVmfxVUaDtgdkzzABSiUttPjvlJI0fHeqEsbsI+O1mgZxsoLQtCZDwb7Na1hVu9JzC3E4Og1PjWqN0Cp+OH45Hf5rn/3vNinuWm0V4h6IZPm4ee1gSClxMckidLBAgI83RrWvi0Y1bC5A1SstA4JdiaDqLtfRaHR2hGxl5HiwnZYTNaTyszNUsohPoeHl/LeMez0gp6pgUkophIgHEA/AAKA6gL+EEOullM85/nSJEAPA8oppCMA2qeDMMpBSzgUwFwB69uzp/K9nMpgTbwaLxkdp0i5D2iufGvLJJHXGAyrIZpKvqhvJclfOkGZfzuVmHy9ht4M9PZd30wEtrCdpB/jS0ny6/iw+Xc8u7RbiPnT3OoePXnu9yIijJ/44DADo1awGziTY56gEB/igfYMQtG9g0yKj3s2dVScoayQqo5ZdNUBZSqOpVgUWsDBAGjlAk5HaVvZgVMI79Voo7xBMCPGYECIcwIcAdgLoJKV8CEAPADe4dWvM2A+glRCimRDCD8B0ACtsllkB4HalGtYXQLqU0n1qVCajhQGyIHj5+GuPQTEZtEW91ASfjxPaxWr1IaiW4+WchepR2WsPKUM8OESjIVdBvoEXflKWdUe7v483Zg2z3tYePfrgo3c+KmJ8dp3nb9W7WQ2EhQbirX958Q5sWfS387anJaR6PlXc9HuXFpfDmUN01VvNiKURCbKTO1L5P/5F82MoyLJ+3ZBnHaoZbc51d3HTFDjjAdUCMFVKadUhJ6U0CSGuc+vWmNdtEELMArAWgDeAn6SUJ9SxzMp01FUAxgM4DyAHDA9LDy9fhlPSaMEwtagS+ARYVw1UePtp8zDUO4YzHohKAHPXBZGXzhOsAsKvYW3rYPZmbf2YozHsA7ugETI9M6YNnhmjoTRpgUNRqbj31wNoXisIP9/ZC7FpuTgQmYqR7eog/FKqw89aQZ304CkjiCJ3AE0GuK5OmXyeTaO27HwVamFFy7O2bVAtyLIO1dRzWs2H+rqpOKLAmRzQaw7eKxulaq57FWhkLF/71uJvCeARt39xYHVqp5gMdM29fDU8IA0D5BtI4mIRpqhqgJw4qdRkobsqYPlK1UerMlfG6NLQ/nd+uv4snhzZCh+vO4u/D1/GpK7O1w5WHYvDwwuZNP7l7t4I9PXGDXN2Id9gQpC/D1JzrG8CL4xra39lZ9fywnW2QlmWyIilt9rzHtc/m1xM4SIzjukBLY5afibF9wCeuwXZ1uefmn74zwC5dwqJR/KAKhRqCbRAMQaBodbdyT4B2j03PgEAZFGX1ZWY+b+qm5sI6o5i/zKGj7cXagZp35ETMvLRp3lNNK1ZBc/+dRSbT2uIZdngVFwG7pm3/z/js2LWADSsHog3/z2J0/GZmNItDH8fLso9mtLNjnHJTmLTbzsXBgSUJc6t43NzJ/WiVJhM7DV0FGZnJlAmxtazMhmtPSBjAc9BSwOUl8Znbz8Awu0GSG/FsEX1puTOqHT0KrXMVQaABkntiLeEZT+OjwUhTD2YzhC4/uMduemw/FdCLZuxusVh9i3dMeP7og2pADDt291Y8lA/PPrbIdw1bz8GtKyJSV3C0DGsGqoH+cJglDgYlYqY1FxsOJWAQ1FpqBrgg+fHtsW9g5rBJCWeWnwEyw5dxp39m+JkbIbm99QNsdM2sG8uQ+3ud7hrd0uHkyt47qnTKZxF8jmeW/UcfC41UntIZlYiAAlUVYiPalOqZeuRyvD39uW57GZjrRsgW6j5EsuEsGVvTFAtha1sg/+kDFKs+Rg+/jQEznQSqwdXq2enJPgvh6Vdni5r9G3uWOz+hjm7senpIVh3MgHzd1/6jxtki7b1quKVCe1wQ/eGqB7khxOx6Zjw5Q4AwKPDW2Lpwcu4nFbUK13ykB06Q04KsPc7kjMtxy1VFHJTOb2078OuX+DqdNOwntrvSwkknuDcMVtkKh7jf8JkSlLeUv9KNUomk/kcdyN0A2QL9UdW1d+CagHxFiNkgmrTfbfN9agcjKwEoKZNBci/qnMGSL3zuKvfpkpNAKLCmL5CCCx5qD9umGNncieA4Z9sxbiO9bD8kQFIz6XWT0au4b8pGANa1UJYaCCklDgYlYaflh/HquNxqOrvg+fGtXU43qdHEzvEvE1v83gMf7VU++c2HF/Km06Hya5/9vIBnje17BjSjFgWI7Q0oNRzXNVIUpPVwZYGSPGATAb7Zf5SQDdAtlDd0YQTfA6qwwOjGpwqtei652dYJ3fV/p0sjXyGf0hRLRYtqLF4vnY44TJ8/JlgtR1KV46wawQssPp4PFYfj8eo9nXRrXEomtYMQmJGHuqGBGD2pnOITMrB7gv0QoP9ffDQkBaoFujr0PjsfGG49huXdnMGe+/7tSeblDekBA78zBHK6owvVxC1h2Oe7FXO1PPYkQEq4gFZ9KJlXWF6oSBLW4a4lNANkC1UQ6IOeqvehAZBFfz+T6w+1toAqRyMLA1vIyTMOT1dlS1tKQhVWtRpSzXGCsTZt8eh9SurHS7TrFYQTsdnYP1J+97aKxPaISTA126opuL723tqz4zPTQWW3MsRS8NfcWrbyxyXw4GEY8B1n7kefmXEkrzYZbr9ZRIVA6QOGbRE0jlWwNTCS6ZCo7Nshk2PYmtIVmKZjAbXDZAtbDuRVQ2f1Is0QNWVrvaUi9YHNagW7xSpkUXXWS2MmjzOfLeXr3ulLxv3ZciRk1JuAwht4efjhfBXRqLH2xvsLnMxKRveXgJjO9RDSKAP0nIK0bB6FeQWGpFfaMS6kwl4e2XxrI8XxrXFqPYa3eRGA7D0fiArHrh7nTYpryKwZw55N52muf7ZCGXUUcuR9pdJOMkboFafWOIpnsOq4Uu5QKlhyyJKWjTPy6Rz1rkhN0E3QLaw/JHzMiwMUCTvAKqshm1YIwRzP1rD21QPyGR0TPv38mInvKuawI7QdDCfz28AOt/kvvW6iJrB/jjwykj0dGCEjCaJNSdK3pf1xMhW9hnYa19kqXvCp0BD99/JS4Sk88CJpUD/R0tmEM9vILnQUetGwnHt96Wk99Ruovm15PNF85fp0crnpfnm60boPCBbWCba0i6ZJ6KqqoaB1anXnKKRV6nRwqz/bPV6MzKrtbwjW1Rvysms7kLDXryrHfvLfessIWoF++P0W2PLZETO/Ht644mRdhKxWz9k2b3fLKBXCYh+ZYUdn1LDp98s1z9bkAOcXQe0Gm0/dMtJoZejJW6XlcgEs2qcpCw61aUgmxVgtSrrrKaVC9ANkC0sPZTkCPYghYSZO5WF4IHQ8nRqtqSRse0VUxOAakLQEep14kmj1W9WEnh5AZ2n8W7pjAEsYwT4euPI66Px1iT3TObw8/HC0TdGY1ArOxWaLe8Dm98ButwCjHrTLd/pFiSdB47+AfS4w3X9HwA4t5b0ik4aE3pVRO4AIIGmGlIsauOrmojPSWa1zNIA/TemWTFwugdUTmgykM+JSs6hbkfrUnzdjkzs2rKc63agp6MeXBW12wEQRV/XQv0uZKQ6s6yz6H0/t+vXye5bZylxW7+mOP/OODw/1kGrhAO0rVcVW54ZirNvj0NIgJY8ihFY8xKw5T2g60xg0mzP6XoHgHWvsEF50NMl+/zxJczNNHEw3SNyO/OSWslj1SNWPaAE5fyuY3E81JuubwDXUxJDWQx0A6QF9a6gVhDqdaRSodqUWq8TiYq2Fa8G3fgce8j6db8qnJbhTDVKXYfWKN2SIkQps6ZepDyHh8DH2wsPDW2ByPcnYOuzQ3H3AMd32D7NauCL6V1x5u2xWPPEYDStZadnriAb+OM2YM/XFB27/ivPMj4Rm4Czq4HBz5Tsos5JYfjVYYrj/bq4jUUIrSbV8+v5rCovxB7mc/2u5mUST5JEW5gH1G5bJi0rehJaC+ponYvb+Vy3Iz2IhBNMYKrd03FHrZPW1ZsyRxR7CEWa8xv2AiI2Fj/WpHozsrEvbAF63eumHQLwzHlgdg9g5dPA7Stc77guYzSpGYTXJrbHaxNLyc2JPwZ8q3iw4z4E+jxQ+o1zJwz5wJoXea70fahk6zi8kDrQ3W+3v0xyBD0YrekqUlLn3JIdHXeYPB/LSmniKTa5Xjllnax2IzzrLPQU1FAqAXlpvJuqc+IPL+Cz2q9j25IhBD2YyxoSn036sadMK3dku44Ww4AL29zbwxVcGxj1Ft3yHZ+6b72eAimB1S+Yjc/kOZ5nfABg28c0DOM+si53OwuTCTjwE4dWOppwe3oln9tOKPpe2iW2YVgqb8YesvZ+AN5wA0PJn3J2qKaL0A2QFiwb++KO0OBUqWlupwgIYV4nWqPRsnE/pZnVhkyoxuqXdhb//S1GcIpplBPcIVfQ/XYKsW9+hy781YLMBOB/ocBeZTDKE8eBrrdU6CZpIu4ojX+XGUDr0SVbx4XNzsl2nP6XRkOLvaxy0lQDlHWFBQo1/AdIG8m4bBYgq9epZNtbDHQDpIWq9fFf5j96L72S5kNZVVATz036A1F7i3opzYcCkIy/LVGzJXtsIjYX//2tRrN948ii0u2HLYQAJn7Bk2nxbc6RIz0ZJhOw/0fgM8UTaD8ZeC2FXCpPQ2EesPxhavKMebfk69n5BZPP7R2ME8qIBaL3AW3t6AVe2MpUQW2FSKsSGpsPNS+jzqdTZTgceVulgG6AtCCEWZcleh+fmw4kVT3lAv9v0p8jbRNsEssNupPeftFm4J0QQJuxLIdbakxrwa8KE4wnlpunVroL/lWBmUuZZ5o/FTi9qvjPeCIu7QJ+GAGsfIqJ1kf2Azf94lnJZkuse5nnyqTZJWekR+/nedX/Ucfh25FFACTpF7YwGZmAbjnKnAeM2EgP3zIEi95HjlJBNiOAMmKO6wbIHlSXM3InvR6VSxGpJKbVkCrSJqTy9qGxithUtEzfZjyb+tR1OELXW8jzOGUrhe0GBNUC7lrFkusftwI7PqM3URkQEw68UQ34eRybJ6fMBe74xzNkNezhxHJg/w8kHLYeU/L1bP+YnksPB+rDUgKHfwMa9zcXUyxxOZycH3U7TCaeqy2GWxcmYvYzFRF7mAa+jOBxBkgI8ZEQ4rQybnmZECLUznKRQohjQojDQgg31qwVqNyJ/HSWrmu2pOsbSR0ahNQHarWhR2OLNmMZU9sSD5sNJp/i1D/Ff3+jPjyBDv5aqt2wi+A6wJ0rqQi44Q1gwRQg9VKxH6sQSCWk/XUy8IPS5d6gO/BoONDlZs9QNLSHxFPA37Oo1zPi9ZKv5/JB4OwaagbZCPRbIeYARcrs5cDOrmEFrOUI/h9/lMWRFiPMy+Rn8fv8gnkTvJYMEID1ADpKKTsDOAvgRQfLDpNSdpVS2lFjKgUaW1QILmzmSd50oHUeqPUY/m+r9dNmAjWgbb0X30DG5SeWWQvda0EIoNd9QNRu5prKAn5BwLR5wMQv6XJ/3RvY/J77w76SIjcN2PMt8HUf4JeJTO6PehN48TJw/+YKkZp1CdlJwG8387hPm2dfNL44SAmse5VtQsWV7sPn8SbX3naQsLKe06toUNTm1IiNfG5hIV9yaSclZ1RBu8b9SrbdTsDjDJCUcp2UUs3s7gHnfZU/qtY1l+PPKwep2WDmgVSWcuuxPFC2ieXg2jRgJzXCp663kPJ+xoncS487GJtv/7jk+1EchOD3zNrPEHHr+8DnHdnCoE7pKE8UZNNAL74d+KQtsOZ53vEnfQM8cQwY8LhjD8BTYMgH/pjJMHHG76VLjJ9ZDVzaAQx9wXEuJusKcGwxz7GAkKLvJ5wgp6fDFPNr5zcx3WCpARSxiSztwlx64Vpyrm6CxxkgG9wNwJ6QjASwTggRrkw/dT/UMmXkdnosrccBEMCpf/l6oz7UBDq7puhn21/Pg622c6hoNpi9ZYd/K/77/YLocp9bZ06GlxWqNQSm/Qzcs4F3vC3vAZ+04Tz2DW+UnVckJfui9n0P/H4L8FFL4M87mWTuditw/1bgvk38282C6GUGkxFY9gC918nfAA1L4aAbCoD1r1HxsDj96gM/sY1HHTlti2N/Um9cJSDmZZBKYhl+AbyhhvXgOddyVMm33QlUCBNaCLEBnDlvi5ellH8ry7wMTmFdaGc1A6SUsUKIOgDWCyFOSym32S5U4tHMAI3FofkcwxO1i25qoz7A6X+Aoc8z4dx6HDkXhs+sKxMdbwDWvgwcWgCMecf8upc30G0mO7STI4rKH9iiz4O8ONe8QONQ1gzmRr14x048TeLlrq9ogHfNpkFu1IcnZ/3OpCu4kn8xFlLrKPk8k5txh5lrULWJQxtTXKvDFGU+lodWtBxBSuCfx+nFjX7bcbOoM9j1BXM6t/zpeCigIR848CMpHCpx1hImE/u/WgwHghTZ4cjt7HS31BNKvQQknWH7jiEXaHUVGiAppQMFJUAIcQeA6wCMUOZ/aa0jVnlOFEIsA9AbQBEDVOLRzID1gTm9kgev3XVsJEyNJJ2+0zTg6CImoy1Zp0G1gDbjWBId8bp1/N/zHlae9swBJhQTXvkHAyNfB5Y/RPfakfqdO1GnLS+g4a8xBxa5g3fG7R+bRxr7BJpnpnWYSm/Qyxv/cagKspjHyUtj6JoWzZYWFTVb0qg16U/2d/Vmnp1QLg5SKjed+cDg51guLw2SI4CtHzGfUxxx8dAC9ib2fVj7/Us7gIwYnksqTq6gtEyjPubX1AKJjz9HTTUdWLp9KAYe1wsmhBgL4HkAQ6SUmrNshBBBALyklJnK36MBuF9roUoNHpzovTww4z5kEnndKyytDnyC5K0qtYCji4vS3rvdxkT02TXWxLGqdWm4Di8Ehr1UPC+k83QS7ta+RHdZFcAvD/j4Mdmulm0LssnoTThOI7x7Nl+PP8rclsmI/4Yx+gVTZjYwlF5TxxtptGs0Z95BK09RWSElj8+eb+i1Dnup9Ov79wkagrEfOF7WkA9s/xRo2NuaTGiJAz/xWKjkxIJsntOdbrS+OZ5aAdTtxHaRpoPKPOz1OAMEYDYAfzCsAoA9UsoHhRANAPwgpRwPoC6AZcr7PgB+k1JqJGLcgNZjaICyEvjcpD8bS48sYkLU2wfoOJXl8rwM64uq5Qjme/Z/X5S52vdhGqC93xZ/snp5kcD23WAS7276teI8Bb8g9rU1USojluHltQqTkcbi4K9An4eAse+V/viE/0zqwYRPzFMr7OHQfHo3k77S/t7MeBqb3g+YK4en/mWJvfPN5uUy4niOt5kAnFkJDHqmdPvgBDwuCS2lbCmlbKSU17tKKdV58LGK8YGU8oKUsovy6CClLLuroPU4898nlvO56y1MMKuzxbtMZ57o2GLrz3p5syHy4jZ6DZao15EcnN1fA9nJKBZ12tFQnVrh/hYNHSWHIR94swaNz+Bn3WN8rpylllGL4UCPux0vW5hL76dRX6D5MO1lDs5nrqenxbqOLgKqNbYusZ9WiismA5PVWo2sbobHGSCPQ5125p6ZE8uYSO0wlTR1tZLVoDuFxPb/VJT93P0OhiJqqGKJYS/TFd75mXPb0v8xiqX9+2SFT7rQAd44flX4Ns2GcNJGaY2PoQBYcg9Dn8lzii867J7NplF7320ooDfVfChQS1E7TIum3Evnm6zXf+R3nuuJJ2n8ymGIgW6AioMQ5kpGdiKTzYGhQNvxLGsW5nGZnvdQwCzKpkM+MJRG6PiSotMu6rSlC7zve+dYyF7eLJUHhFB2wp3i9Tpcw6XdwEfNWcW78SfgDje1zGx4g/m0SbOLn0KREQds/4x5nWYasqsAvfKMy9a60+E/87mHRVk/8TTbNGq3oRC91iTVMoBugJyBZSn1kKIJ1ONO6qScWGpexr8a8z226PsQmdHbPir63vBX+N5aJ5OWwXWA6xVv6ttBFCfXUb44vRL4eSz/vmMFKRfuwNE/qeLY+37nwp+N/yMRdvRb2u+bjMCOz5nwVyu6hnwg/BeSaC2lOg4vZNhlLCSTuhzCL0A3QM6helNWGACyUrMS6XLXbstSupRMznabyTyRrTcT2ogG69ACcze95XuDn2X8fc7+yBortB5Nan9eOvDXXe4VLtNhHyYjPZRFt1Br57HD7uuTijsKrHiUDHpn5Dqi9jJk6vuQdtMpwHMq+Rww8ElzeHZiOeWELdU2jQYK5Dfpz277DlPLrUKpGyBnocpfSiNzP0IwwRx/lJUDAOj3CL2ZXV8V/fygZzh0cMv7Rd/rN4vSl6ueYU7IGXSYQg7R2TXA348o5W8dZYa0KCabd3zGm8k96903piY7iaoEgdUpKeKIcAjQi1nxKMctDX5WexmTiR53jebUSQJ4o9z7LVuMLBPWZ1ayyusXTO6WI6lXN0M3QM6i41SGWABzNkYD8zcB1cj9ANgz0+VmlkVtZ8RXrQv0uZ98IduKmI8fMPFz8mrWu9Ax3etehnBHF9ETctcoHx1mSMmq4xxFfqXPQxR18w1wz/oLsoHfbuL5cvN850Tqt31MtvLEz+33hp1YykLFkBfMjPILW4DYgyRIWiaf93zLcCwznl691hyxMoJugJyFXxDQdQb/zohhO4ZfEEubJ1eYZygNeIJ3KNUoWWLgU2wuXfVs0WpZ04HkBu3/3tz86gwGP0uX/eTfwKIZRTvzdZQcWYnsS1v2ABUBHz8KjNPwYEsKowH4627qMd/4k3M9YwknKOva+Wb7bRKGAmDTWyQUWo583v4JVTktpTrijrDNKKwnjVP3O8qVY6YbIFdgKQS1R9Ef7vsIKevbP+H/tVoxKbn3O95RLBEYSip89B56QrYY8So1hv5+xLVO9H6PcPRMxCbgx9EeMYCwUsNkAg78DMzuSdWCEa9RO6l6E/d9h5QklZ5dA4z/yLmkb2EesOQ+MprHvGd/ufB5PAdGvm72dKL3sffLVk1xz7eAbxBFygKqAd1vK8VOuQ7dALmCOm3N1YTovRR/Cq7Ncb9HF5sTzMNeYlfyVg0KfdeZ5A2tf5VJZEv4BgI3/MDq2l93u5bX6X47MHMJS65zh9Hd1uE64o9TbfHfJ5hofmgXhwe6szFWSmD188DBX7huZ8cvrX+NVI8p35obSm2Rk0JJlaaDrHsZt35IPWrLMT0ZscDxvxhyXdzG7Sgj6VV70A2Qqxj4lPlvtaze/1GWMFUvqGYLekvhv1BqwhJeXqTXZ19hT5kt6ncGJnzKasSmt13bthbDgfs2U7jq18ms2Oh5IeeQEUflwu8GcQjl5DmUetXqLC8N1J6xfd+x+DD8Vec+d2YNP9P3Yccd6pve5g3MkpF9cTt1oG21lHZ+Yb7JefvZl/EoQ+gGyFU06W8uyZ9dw/i9aj3mgg7/RkIXAAx5jqHZBo2kclh3spoP/qpdeu92K+9UOz4lgdEV1GxBtcDut7Ni89PooppEOszIy6AK5FfdWYru+zDw2EHmSdydC5GSN5093zCZPfpt574jLYpqCHU7ASPfsL9c7CE2nfZ+wKxpLiXPwZAw6zlpmQkM1VoMowJit1vLZPRycdANkKsQwjqJt1XxggY/yzKmanCC6wCDnrLP7xn6IisOKx6lZIUtxn1ITZxlDxYVvi8OfkHA9V8CN80nJ+nbgfSGdNKiGblpDEs+78SQpfUY4JF9bK5V5UrdCaMBWDGLrRO9H3C+Z6wgB1h0K/uzps2zPw3DZOTU26DawDALFeNTK8hwHvaSdWf7ri+ZJvAL4rotmdLlCN0AlQQ97jSfpGdWspIQVJMG5+wa80jn/o9S82b1s0U1oH0DqJaXFc/eLtuqmI8/cPMCkiAXzTB7Vq6g/fXArAOsmOz4DPimD0WpKssEjLJAZjzwRRfggyYc0NikP8PWafPcx+uxRWEu57AdWsCy+LgPnDM+UgL/PMZy+g0/mHu5tLDnGxqaMe8wmQywGrvhDfZ3dZlhXjYznp5S82GsntbtULwwXhlBN0AlgRDALRZVrPWv8WTp8yDnba17hRe5jz8rHCkXqJJni7AebEg9sZQnhC2q1ABu/YuNr/OnUKDKVQTVpKG7cxWHHS65B5g7hD1t2lpvVydiwoEl93KIoVolfGA71R/Dupfd92Yn8didWQ2M/5jeibOh3a4v2W84/BXH43wSTgIb32RPmGXZfdeXPPfGvGOdRN/8Llsu/IL4v7N5qDKAboBKika9zQ17F7aQu+MbyBg97rC54U+dNnB0kXYuZuBTrFasebHorHmApd/blwPGfOozl8QIAUDTAbzgpn7P6tuCG4DvhzHHdLW2cmTEUu5k7lCO8zmzhpNGHj0IvJHOhH9ZIu4oK5Jqw2rv+5z/7LG/eGPrMIWVMnswFADL7qfXM/ELs3FLvQRs+4RqiuoIHoDn4KH5VOs8tw7odFPpZpWVEroBKg1GvoH/5EfXvcILueMN1JLe8D8zG/rps/Rilt5XtCrl5cXhelVqAovvKDpTHqCLfPsKag79MrHkRsjLixIMsw4A131O0uJfdwNv1QT+vIsXbGVHdhKHAP48Hvi0PStO0kRVwadOkkhYHuHGiWXAT2OYX7l7DZn0zuLiNiadmwwAJn/r2GPa+j5DtIlfUgZYxZoX2BZkyxda/zon90oTH8M1KrHlCN0AlQbVm7C0CVCg7NB8niwTPqVWstrhXrUuZTTijwFbNBoNg2pS5TAjFvjjNu3Seb2O7Lw25AE/jmK8X1L4+AE97+I44+mKptGJpQxP5k8BNr7lObPBnEFeOnBoIUdNf9yaydjsK0z0zzoAPLAN6Ptg+TRYmowMh/68k5Wo+7e4FuLFH2fSuUYLYPpCxy0f5zaQ+tFtJuVhVJxcQQLlkOesR+pEbALOrSXp8cwq8n7cSa4sAYQdzferEj179pQHDrh5iGpeOvBlNzJJA6sDs8JpUDa/SyLizKVmF/jvWUxE3rXaLGlqiaN/Akvv5Ql1/WztO1/SeU4xzU5ilauVQ31/55EcQRqB5QyyliPpqjcfVmFJSk0Y8jk6+MJW8qUuh9PTCG1MD7TjjfQay1u2Ni0aWHo/Wxu63Ua+l6MZ7rZIPMUw28cfuGcdRyXZQ04K8KGSNH8h2mxcs5M4yLFaGHDvRnNja2EeMKc/AMlK2ZUzDEXtERpLASFEuLPDQj3OAAkh3gBwH4AryksvSSmLTPFTxOu/AOANakUX26RTJgYIoMu/UonTu97KpG9hHklt+VnAw7tonPIzWRKXJuDBndp3ZNVwDX8VGGxHkzczAVh4I/uCxrxLfoe7LjaTETi5nHmLM6vM7O7gulTVa9yXQv21Whffte2u7Uk8ySpg7EFuV9wRepjCC2jQjSFv2+uY1K8oreyTf5NSYTLS8Lg6veTKGWDeBBJa71zp2OCbjMDCaZxqettycnkAFhX+vIMJ7/u3AnXbmz+z5X3Oeus6k+OWJnxKBn8Z4GowQFlSSrvzaoQQ3uDY5lEAYgDsBzBDSnnS0brLzAAZDTQsV5Qk850r2Vx6+SDww0iKlU2dy/ei9pDq32kaMOW7oheMlLyLHlvMqom9xGVeBjlCZ1ZyasbEz90/wUBKXvyb3qH4/qXdVIUEKC1SqxUla2u24h23WkNKRATXZp7B24mZB4YC89iezHiGoekx1LFJOmeeQgtwDFD9zjQ0TQcyRxIY6t59dhW5acA6Zf5bg+4sl7vqLV45C/yiTKu4c2Xx7OtNb5OFf93nDKVVHPuLVc4Rr5MSoiI5AvimH9B8CBPjIfXpHZXR3DVXDJAnTsVwBr0BnJdSXgAAIcQiAJMAODRAZQZvHxL/flQo8v8+yYpTWHfG4Vve49jjDpPpQQx5gbmgsJ6U6LCEEPSgCrKoD+RfVftuGhBCntD2j+k1JZ5kpcWdrQNCMJSZoeSJpARSL7KknXiSIUPMAeD4Uvw3iscSvkHcTh9/AEIxtoIqfvlZ3EejRr5LeHFGWK3WDAPzM8m9qtvROaNWXjj1rznfNPBJYOhLrs9/jwkHfpsGCG/gzn+LP36n/qXx6TbTuq8rNZLNrWE9ybJXYTKZx/v4V6Xuz4zfPGboowcdTSvMEkLcDuAAgKellKk274cBiLb4PwZAH2igVJNRXUGj3uyM3/M1e4k2v0OpzEFPk5z47xOUW6imiEjFHgLWvsjksjoCWoW3L3Djzzwxlz/MPh2tKoqXFw1cvc7A8gc5tmfMuzwxyyIUEYICV7YKfIYCejAZl+m9ZCfRaORnMEdmyAcgFd6RZJjhF8y+JL+qLCFXrcdpnFXrMdwrj/CupMiIY5Xp5HImmm/5A2jQ1fX1nFsPLL6drPmZS4v3nC4rXKYG3YHxn5iPsaGAVUwJ4MYfrY30/u9ZVes6k7KrPe+iB+khqJAQzNFoZgB7ACSBP+dbAOpLKa1mkwghpgEYI6W8V/n/NgC9pZQOR1GWWQimoiCHib7Ui/z/rtU0LskRwHdD2E1/5yreJfPSge+H8/n+rdbVChX5Wcz1RO8FJn1t3QJii4xYlm4vbKG3dd1nxYua63ANhXm8wWz7hF7ckOeV2XAlMJaHf2NRom4Hkk2r1nW8fOolhvO+AQyfLPu21rzE7bppvvX8uaRz1A1v2JPnh7GQ+cgy7nh3JQSrkDK8lHKklLKjxuNvKWWClNIopTQB+B4Mt2wRA6CRxf8NAVQ8icWvCqcZqFj2ID2Bmi34esx+c69YQDWWwAvzqIhnK80B0EOYuYRJ1uUPUYnRHkIaADOX0QM6vxGY3ZvJ8Wu57cJdkJKhzzd9WGJvMQx4ZC+LBK4aHymB32/h8Ww2CLhrVfHGJzeVSWdjPo2VpfE5uUIRsn/A2vgYDRRS8/Hn8ikRHFxYznIbxcHjeEBCCMsxkFMAHNdYbD+AVkKIZkIIPwDTAbhpLkop0XSgmRuUdglY/QL/7jCZHdB7viFJDeAIlJvncwyuPf6PXxAw4w96NaueYee2Pa/Vy4viZA/tAhp0YX7ip9FA7GF37+W1ASmBiM30PP64leoGty0nP8eeELwjFOQwhDqzkjyfW/4s3iDkZ9L4pFwAbl7Ic0ZFTDgLFmE9i07G2PoBQ7aONzBH1/Me+2ObKxJSSo96AJgP4BiAo6BRqa+83gDAKovlxoOVsAgALzuz7h49eshyQWG+lN8NlfL1ED4OLTS//v0IKd+uJ2XsYfPyh3/nckvuk9Jk0l6noUDKZQ9xub/ulbIwz/E2mExSHl4k5QfN+ZnFd0qZdN49+3ctIHKnlD+N52/3SXsp9//EY1BSpFyU8psBUr5eTcqtH9o/zpYoyJHy5wlSvlFdypP/FF3fhy2k/KyTlJmJ1u+d28DvWXizlB+3lfKLrlLmZZZ8210EgAPSyevd48rwZYkyzwFZIuUC8E1/8lW8/cy8jMwE5n6kCbhvk3nu97aPWF7t/YD9bmkpyXzd9BZH6t68wJp+r4W8dGDXbPZEGfOpEzTkeT0/pAWTkTO/dn0FxOxjMnzQMxzg5wqh0Bbn1rMNR5qAG350LCimwpDP8T/nN7K0bzmbLjcV+HEMK1r3rAdqtza/lxFLSkhQHZ4b0Xu5TEmS5CWEx+eArgnUaM7mQICl5j/vYFK5al3glkU0DL9PN2v0DHqGmiz7vmNfmdaNQQjmHW78mRyj7wYD0fsdb0dANWD4y8Djh6nSePBX4IuuwKrngJSL7tzjyov8TObXvupB2YysBPaOPXaYNImSGp/CPIbgC2+kINh9m50zPoW5wB8zqVhw/ZfWxsdQwHA95QJDQUvjYywE/rqH39uoNzWgR79TrsbHVegeUFlj1bPAPoWE2G4iMO1X5mrOrAZ+n8Hczk2/snSqagXv+468khGv2y+nxx5i82pGLOUWet/vXOk95SK9raOLOeOs3UQavoa9Ko5FXBGQkhymg/OA48uAwmyFQ/Mof5PS8mQST5MUmHCcMi0j/+fcKJ/8TOC36VQpvO4za6KhlCxsHF1EVYPON1l/9t+ngAM/Mt9zUJl+evOCcj+ulZoJXZaoEANkLAR+uZ79QQA9nRGK/sreuRQr63ILy+xeXjzJ/n2Sch79ZjmW7cxNBZY9BJxdDbS7nsxYZ3t7MuJo6A78RG+sdjtOROg8vUz6gzwGqZfI3znyBwXefYPIsepxp3taOYwGqh5ueY9cp8nfOC93kZNCbyn2MJnzlp4PYG6nGPYy+V+W2Pc9ixQdb2CPXEA14L6NZaPuWAx0A2QHFWKAAMpyfDcEyFSYAlPmcoAhAGz5gKzoPg+ZZTpNJhLd9n3H3rKJX9pnAJtMwO6v2MEeWJ3jedqMdX7b8rM4GeHgfODyAbZYtBpNHZk2Y83qepUZKRcpjXtimVlFIKwnDW7HG9xXmo47StnVuCPsTZvwafEldhWpl0jHSLkATPvFursdAPb/SKZz11t5s7I0lBGbgAU3sgKbk0wN6Xs3Wodn5QjdANlBhRkggCXTH0cx7AGAu9awI16dkrDnGyaHhykSHlKylLrlPaB+F0XRMNj++uOP0T1POE6a/uh3XO+TSjjJnqYTy2gsvXzJeWk9liXcGs0rR5hWkMMQ5vwGPpKVyST1u9LbaT/ZvTIU+ZmcVrrrK+o6jf+IBtzZ3yrmAPOBxgKGTM0GW7+vNju3GsP3Lds9Ek9RdyhE6cU7t55qna1Hu2//XIRugOygQg0QAJz6h8lFgC0I96wlE9ZkAv55lBf/oGcoEqWevHvm0BsCgKfPOK5eGfJptHZ8xgth9DvME7hqNEwmekMn/+Y2p13i66GNaYiaDmLOqHpTzzBIWYls8o3ey0fsYTKVfQLpFbQcyQuyJNwdRzCZOEljwxvU9u56K0PmKjWcX8eJ5SQMBtclydDWa1FDq9bjODfeMiGeeonGByDb+dQ/zDUNfKKUO1Y66AbIDircAAGcRLnmef4dXI9qeTWaKU2Dj7NK1f9RYNRb5ot79zc8yavUYOWjuF6euCNMSF4+ADQZCEz4mF3rJYGUbCW5sJltHhe3A/kKazsgFKjTnvSCOu1pTGu3LZsOdSkZXiRHMExJPkfxrvhj5tDW25/yHI37AM2GsA3G3QoBKqL2spfvcjiPx9gPgEa9nP+8yciczrYPKW8y/beilAr1XNEyPllXaHxykoCWoxhG93kQGPt+hd8UdANkBx5hgABz7w7Aru+71zJXYDIBq59jA2Hv+3lSq6N1449zOkZmAvM8ag7JHkwm4NCvlODMz2BYNvRFtmyUBqo+T/RebpPaFZ+fYV7GP4QhQUgD9rgF1mBXvH8IjZZvADvehRe7wKWJfKmCHKAwByjI5oWVlWh+pEdbf4fwJiu4Xic+GvVhqFoavo4ziD1EeZLz63kDGfU/6ip7ucBoyU5mhezCZnpNEz61rpBJSU7Y9o+ZS7rxJ+v9ysugfMeVs2TYH/mdk08mf+vadpQRdANkBx5jgEwm4O+HeeIA9Bru+Ic9O+rwut2zKXo/eY755MxOZvf0pR30kkb+r/hycXYyyYv7v+cF3+dBtoq4EiYUBynZBZ94ksJaald8xmXSBHJTtWU3HMEngL9HUB0+h4Sxp65Gc7YxhDZ2XfqiNIg7yvD29L9M9g94nDcJdbKEs4g5QPpE9hV6pt1vt37faKBywqH5QPc7aJwsCxB56ZSejTtMTakji1g0mL7QYxQEdANkBx5jgACeaEvuYUkYoPbNHf8wxyMlR6qsf42iW9MXmsupxkLmhPb/wPBq6neOpTtVpF6ibtDRP3jR9LyLZf7yYkQX5lnIc+TRk1KF0QHAtwqbeX2Vh19QhYcS7AXbyOTyhS304PrNAvo+5Lq+tNHA3NzW9+kZ3jS/KEEwPxNYch9pFYOfY0HC8jfITQMWTKUx7H47w/WGPSnl4VellDvrPugGyA48ygABNCaLb6f0KcAhhnf8a27POPYXu6arNwNu/dO6cnP4N2Cl0o19/ZesujiDhBO8EI4voS5P11s4jtiyyfFaR0E2Gzj3zCFXKLgeRe0tB1K6guQIJppj9tNrGf9R0fUkR7D1IukcW3FslTBzUzkwIP44j9nBX5jzeuZsxatC2kA3QHbgcQYIYOXqj5mc0QTQ2Ny21FyxidzBE1N4c7KGZUdzcgS7q2MPUgR9zLvO35lTLtLLOrSA4VGTgUCvu4G2E8s3tPEkxB7mhX3sL3prdToA/WdR5L4kv4nJRGby+tcZRk34tCi5EGDp/K97GE5Pm0fpVEtkxJGgmHSWnk/4PIrwP3rQs4YFKNANkB14pAEC2N+z9D5zOBZUmyVZ1UX/7+54Fhj1JsMA1TU3FjK02vk579QTPilKYnOE7CTmGw78zHJ7UB3epTtPI2+mosOgskbqJfKeji8B4o8y99RhCvMvjfuWfP8TTnAabsx+3jQmfVNUdO6/5uK3qYx588Ki/KSkc8z55KaQUnFwPquOM5d5LGNdN0B24LEGCGBO5J/H6JEAzINMX2ierJqfSXnWUyuYnJ74hbW3ExPOqQyJJ3gBjf3AeRYuwLt1xEbeXc+uJY+mZksao7bXVcyYm7JCcgR78U4stWBG9+D89E7TShfSFOQwWb17NlnkY95lhcr2t8tOBr7uzWpfxxtZ2bTN48QcoBaQlzeP+f4fuJ23/ulxYZcldANkBx5tgAAagXWvmEv0AOeDdb+Nf0vJ/M2mt1gFuuEnoKEFJ8hQAOz6Atj6Ie/kQ56jvIer4UNuKpX2jv3JEBCSVajWY8jGbTrA45T1HKIgG4jcydL5ufVmydx6ncmM7jCFpMrSQCUlbnyTvKRuM8nl0qo2Ru2lUBzAAQVDXyhqoI4v5Q2nal1SDI7+QU/q5oWOGfEeAN0A2YHHGyAVlmRFgKXz0e+Yy7GXdjNky4xjpWTAE9bl+KRzbO84t44l6zHvsJ2iJB5MZgLXc3YNK0EFWcxH1e8MNO5Psl+j3tYyoRUJKfm7RO9T2NF7WDWSRjKjmw2mJEbLkSSAugMXt3M0T9wRCsaPeVd78KSxkDeH7R+zcjnpG8qyWsJk4kCD7R9zXVVqsJ2k661sNq4E+TndANlBpTFAAHB6FTWEVP5MsyFMUKp31Nw0ds2fWEpDcP1XQK2W1us4t4Fs3aSzbJ8Y9rL2heEsDPnApV3ss7q0iyGCMZ/vBdc1kwLrduTFXb0Zqz1lEbpJSd2e1EiGVAnH+Yg/znwJQIPTsCc9iCb9SWlwRhLDWUTvZ1n9/AbOQxv5OsMpLTJg4ing70cY8nWZAYz7sGjBID8TWPoAJVvbXkfmd9RuHrfBz1aaEFg3QHZQqQwQQNbtb9PZZwQAoU1YCVNbMaQkEW3N8+TZDHsR6PeoNXHNWEjJjW0fc6hgi+GcX+VK24A9GPIpjBZ7iC0R8cc4nNFkMC/jX42J1ar12WoQVItJ9oBQtkn4+NNQ+Phxf0yFzIflZ5LM6BtImYrcFPNz+mUmzA155u/xCWS7Sb2ONIBhPWkMy8JjsDQ8gTWAAY/RS9Vq+zDk87ff8RnD1gkfswPfFvHHOU8+5QINVOQ2Dmq8/ivXp6xWMCq1ARJC/AFAJaWEAkiTUnbVWC4SQCYAIwCDMztc6QwQwBDor7vodagY8y65O+odMTOe3dKn/2Ve4/ov2RNliYIcloR3fM7EZ/NhLDG3GOHeO6uhgN3nqZHMtaRc5N/Ziay4ZV9xnRUNMKEbWIMeYEgDGuPqTZVHM3pcZTlsz2QCzq2ltG3kdjb79n8U6HWf/ZxM5E6ympPOsl1j7HtF+72kpPbT6heYWFb7uqrU5MACD5rh5SwqtQGyhBDiEwDpUso3Nd6LBNBTSpnk7PoqpQEC6MVseIOVFRWtxzKHYFmKPfk3yYnZV5i4Hv4axyRbIj+L1ZQ9c+hZ1W4H9HuYF4g7wxN7kJLeTV4avYPCXD4b8tgq4uVjVof0rcLcUkBoxU1Ezc8k6XPvt/ROQsKAPg9QddCe4UmO4PE6tYLFgus+Y87JFrlpNFAnljGMDqrFzzQdRNld22NXSXBVGCAhhAAQBWC4lPKcxvuRuFYMkIoTy4C/HwUKMvl/UB3OhG87wbxMbhoTnfu+4wU85Hn2LNmGIoYCcl92fw0kHGOupvN0Gq66HcprjzwTUjKBrWojFWYDDXuzBaPdRPs9VzkplLvd9z0HEQx4nF6mVr/YuQ2kTWQlsAoXf5SeUv9HgRFveNYIahdxtRigwQA+tbcjQoiLAFLBCarfSSnn2lnOcjRzj0uXLpXRFpcTUiOZqIzeY36t002k71uWfK+cZQL6/AaGKUNfJMfFNkyRkqN7w3/mRAhjAasv3W4F2k2qtHfhEiE5ggbn8G8c5OcXTOPQ404ms+3BkE+js+1DekzdZjJxrNVnl5fBitnBX4FabdjBf/Jvhl+T5wAtR5TV3pUbPN4AORrNLKX8W1lmDoDzUspP7KyjgZQyVghRB8B6AI9KKbc5+t5K7wGpUBsbN79tfi24Lisrtkp859YDG//HBHGtNkxUt5ukXanJTgaOLSbbNvEEQ6ImAyj50O56zym1uxNXzpDzdPJveoIA97nbTP6WjrrdC/OAwwuYfwMYZo16074HeWY1hxRkXObvmRnPG0nrsZRZLW7EUiWBxxug4iCE8AFwGUAPKWWME8u/ASBLSvmxo+WuGgOkIiacsh5XTptfazEcGP+xdY+QyQSc/octG1dOM+/T/1F6RFpVIikprXFiOdtDks4CEExstxjOu3TDXh4j/+AS8tLp8UVs4iM1kq836kuD024iENrI4SqQn0WPcddsc4Vy4hf0lLSQGslpJ2fXADVbsVJ3Zg3zXaPf5LikSlJidwZXgwEaC+BFKeUQO+8HAfCSUmYqf68H8KaUco2j9V51BghgLkeVeVClLbx8SE4c9JT1HdxkZN5nx+f0cKrWZ0K1x132qf1S0mid+odD8mL2k9TnV5W9Uo168xHWwzPZ0ekxJCXG7Odz7CFl+4NJSmwxnDk0Z4TaUiMpDn/wVybRmw2mhG6zwdoGpDCXUh7bPyF5s911/C3jjtDrmfBp0f6wqwBXgwGaB2CPlPJbi9caAPhBSjleCNEcwDLlLR8Av0kp3yluvVelAVJx5Qyw4jHr3FBwPdL8u91mndSUknf/XV+S3exbRRlNc1fxo2ly0xQPYiMTtar3JbworFa3A/k4ddrzuVqjsi2PqyjIoUxr4imzSmP8cbNcq08gPbgm/V3z4EwmKhfu+54ejPCiwer/mH0ulcnIPNLmd/n9TQfxNz6/nsn+cR+SC3QVeT2WqPQGqKxwVRsgoKhIuoqarcjSbXtd0ZM+7ij5Qcf+YptFnQ4cRdxhqnMJ6Nw0ak+r3kXiKcqnqvDyYdtBtUbk7qgSrYHVmTQPrEECn7cfDYK3Hy9yUyET4sZCJnnzM6wJier4mfRoIC2a3CYV3n7Md9Vpx+Rxw14kJboSMl45Q5Ln0cVARgzJkz3upJG257VIyUbeDW+QkFm/C/c7cjvDtt73A0Ofr5BZXeUJ3QDZwVVvgFTkZzEs2/WVuVUCoLzGoKdpiGyT0PmZNELh8yj3KbzY/tFxKvMirlw0eem8gBNPMmxJizYbi8x4sHBZSvgEKEatEbk21RpRQ6luB/a/laSMnXKRlcDjf9GYCm+GaF2m8zewpzdtMlHFcNtH/Fz1ZtyOhBMkYzYfRrH4Om1Lt8+VBLoBsoNrxgCpSI0ENr/HUb6WqN2WhqjDVO0LNeEEu7GPL+EF5OXLDvhWo/mo2bLk4YPJSAOVm2r2aAx59HSMBfR2pFHxiCy8Iv+qZiZ0YHWGNKUNYUwmIP4I++7OrGIvGUA2eZfp7OtyJGliMjJJv+0T5tSqN6XnpU7tqNsRGP4qVQSu0nBLC7oBsoNrzgCpSDzNAYeq4JmKao2AXvcA3W7XFreSkt7Q8aXsiFfzPdWbsoWjidINX9pJG+UFKdkmcnErxxdHbqchFF5A437M7bQZX3yXfE4KE9H7fwTSo5S2kCYclXPlFA30sJeA9lM8YkpFeUM3QHZwzRogFXFHSZY79Y/16z4BvNv3uoeJWnt369RLTKSeXcfetIIsvh7amK0EDbqZm0E9QTArK5Eh0eWD7EKPPWTOFYU0pPRpsyHk7xSnLiglq1f7f6BOkiGPIZ+XD78nL435s36PUICsEjOZSwvdANnBNW+AVCRHsK/s8G/WHeUAq1ddZlD+09HEDKOBxL1LuykZEbWHDacqqjWmdKjaKKo+V63vvokXJhP73jJjOf4n5QI5S1fO8lmV5VArdA26MyndfAi3x5ltyIxnIvrI78xp+QTS4EoTvw+SXlOfBzmF9RoKtexBN0B2oBsgG2RdYQUs/BdzuVqF8GICtv1kXmDO6A9nJiiyHEeZT0k8zRxSYY71cj6BrLAF1WbXtzqKxzeAz17eytgeyXyQyUhvKy/DPNonN5XGwVRove6g2hxxVKsVn+t3ZTXKFRXB7CTmhE7+TbqCNAF1FWmPrCsMuwKqUTe6173unTN/FUA3QHagGyA7MBqY4wmfp0znsDknhBdDrHbXkUDnipKglPRSUi7SGGXG839VmiMniYS9wjwaKjUh7eXNKpSXMkHVr6rFdNVqDPGq1mf+SX2ENin5wMWUiyyhn/oHiNpFo1O1AZPQ3n7Mf+Wlk0Xe5wF6iK4OJbxGoBsgO9ANkBNIiyaX6PgShhxaCG1CfeLmQ5lD8dDpDA6Rnawko7fwkaY0KVdvRiMmvGiUcpJo/Npdx8pYsyF6mFUMdANkB7oBchEJJ2mITixV8h12ULMlFQgb9iSTum5Hz9IuNhSwTB5zQHnsZ7c7QI8qtAk9HlMhWzcKcxgmthlLxnLLUeWjlXSVQDdAdqAboBJCSgrdn1vLMCVqt7Xsqi28fEgGrNOWyd/abVi6r9aIOZqy0ojOvqKQHiMVIuQpPqdEmLfXy1dpDRFkYEsTK1gAk8utxpDr1GyQtsSqjmKhGyA70A2Qm5CbxnE9qkB9/FFzI2xx8PZna0ZIAxIK/3uEskHU25dGwtuPpWyTUWnJUNoyCrKZgM5NUciMqZyCkR5jzfoWXswRqZ8x5gMQFklrQcMY1gMI687JsLXb6OGVG+CKAbp2yQo6So7AUOZE2l3H//MygOi95NnEHeHDsh/MEsZ8eiRqCFQSePnQQEkJGHJt3hQ0Pr6B1N1RUbUBDU1YD4aK9bs6P8ZaR5lBN0A6So+AEM7aajXK/FpOCkvyyeeV1gTF6KRFFeUeFQuhlOaVMMpk0A4BA6qRYFgtTKmMNWQYGNaj8rC1rzHoBkhH2aBKDRL+mttIOknJcnZWAkOnzAT+X5DJJtqCLIZM0sRlIfksvDi62LcKy9+plyi61rAnDU1IA4+fGKqjKHQDpKN8IQRDuMBQ5lx0XNO49jrldOjQ4THQDZAOHToqDBVigIQQ04QQJ4QQJiFET5v3XhRCnBdCnBFCjLHz+RpCiPVCiHPK89UtMadDx1WKivKAjgOYCsBqjI4Qoj2A6QA6ABgL4BshhJag8AsANkopWwHYqPyvQ4eOSoYKMUBSylNSyjMab00CsEhKmS+lvAjgPIDedpb7Rfn7FwCTy2RDdejQUabwtBxQGABLBluM8pot6kop4wBAeb4KJ+bp0HH1o8zK8M5MP9X6mMZrpeoVsRnNXJpV6dChw80oMwMkpRxZgo/FALAcS9kQQKzGcglCiPpSyjghRH0AiRrLqNsxF8BcgL1gJdgmHTp0lBE8jYi4AsBvQohPATQA0ArAPjvL3QHgfeXZnkdlhfDw8CQhRDaApGIXrnyoBX2/Kguuxn0CzPvltERkhXTDCyGmAPgKQG0AaQAOSynHKO+9DOBuAAYAT0gpVyuv/wDgWynlASFETQCLATQGEAVgmpQyxcnvPuBsp25lgr5flQdX4z4BJduva0qOA9APfmXD1bhfV+M+ASXbL0+rgunQoeMawrVogOZW9AaUEfT9qjy4GvcJKMF+XXMhmA4dOjwH16IHpEOHDg/BNWGAStv8WhkghHhDCHFZCHFYeYyv6G0qDYQQY5Vjcl4IcdX0+gkhIoUQx5RjVGkFyoUQPwkhEoUQxy1ec7lJ/JowQCh982tlwWdSyq7KY1VFb0xJoRyDrwGMA9AewAzlWF0tGKYco8pcCZsHXjOWcLlJ/JowQG5oftVRvugN4LyU8oKUsgDAIvBY6fAQSCm3AbDl3rncJH5NGCAHcLb5tbJglhDiqOIeV2aNpKvtuFhCAlgnhAhX+hSvJrjcJO5prRglhqc0v5YlHO0jgDkA3gK3/y0An4CM8sqISnVcXMQAKWWsEKIOgPVCiNOKN3FN4qoxQGXc/OoRcHYfhRDfA/i3jDenLFGpjosrkFLGKs+JQohlYLh5tRggp5vEVVzrIdgKANOFEP5CiGaw3/zq8VAOuIopYOK9smI/gFZCiGZCCD+wULCigrep1BBCBAkhqqp/AxiNyn2cbKE2iQNONolfNR6QI9g0v64UQhyWUo6RUp4QQiwGcBJsfn1ESmmsyG0tBT4UQnQFQ5VIAA9U6NaUAlJKgxBiFoC1ALwB/CSlPFHBm+UO1AWwTHD8sw+A36SUayp2k0oGIcTvAIYCqCWEiAHwOqhOsVgIcQ+UJvFi16MzoXXo0FFRuNZDMB06dFQgdAOkQ4eOCoNugHTo0FFh0A2QDh06Kgy6AdKhQ0eFQTdAOnToqDDoBkiHDh0VBt0A6fBICCF6KY21AQqD+IQQomNFb5cO90InIurwWAgh3gYQACAQQIyU8r0K3iQdboZugHR4LJQ+sP0A8gD0r8RtMjrsQA/BdHgyagAIBlAV9IR0XGXQPSAdHgshxApQDbEZgPpSylkVvEk63IxrohteR+WDEOJ2AAYp5W+KRvQuIcRwKeWmit42He6D7gHp0KGjwqDngHTo0FFh0A2QDh06Kgy6AdKhQ0eFQTdAOnToqDDoBkiHDh0VBt0A6dCho8KgGyAdOnRUGHQDpEOHjgrD/wHdEVrOFRsKzQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(primary_position[:, 0], primary_position[:, 1])\n", "plt.plot(secondary_position[:, 0], secondary_position[:, 1])\n", "plt.gca().set_aspect('equal', adjustable='box')\n", "plt.xlabel('x')\n", "plt.ylabel('y')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Obtain initial data such as masses and spins" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "initial dimensionless spin of primary BH: [-0.36451101 -0.05518085 -0.64443774]\n", "initial dimensionless spin of secondary BH: [-0.00129516 0.00100747 -0.08374188]\n", "initial dimensional spin of primary BH: [-0.25138458 -0.0380554 -0.4444357 ]\n", "initial dimensional spin of secondary BH: [-3.59770435e-05 2.79857020e-05 -2.32619237e-03]\n" ] } ], "source": [ "print(f\"initial dimensionless spin of primary BH: {primary_compact_object.initial_dimensionless_spin}\")\n", "print(f\"initial dimensionless spin of secondary BH: {secondary_compact_object.initial_dimensionless_spin}\")\n", "print(f\"initial dimensional spin of primary BH: {primary_compact_object.initial_dimensional_spin}\")\n", "print(f\"initial dimensional spin of secondary BH: {secondary_compact_object.initial_dimensional_spin}\")" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "initial horizon mass of primary BH: 0.8304509510300787\n", "initial horizon mass of secondary BH: 0.16666770757317356\n", "initial irreducible mass of primary BH: 0.7588330514\n", "initial irreducible mass of secondary BH: 0.166521231\n" ] } ], "source": [ "print(f\"initial horizon mass of primary BH: {primary_compact_object.initial_horizon_mass}\")\n", "print(f\"initial horizon mass of secondary BH: {secondary_compact_object.initial_horizon_mass}\")\n", "print(f\"initial irreducible mass of primary BH: {primary_compact_object.initial_irreducible_mass}\")\n", "print(f\"initial irreducible mass of secondary BH: {secondary_compact_object.initial_irreducible_mass}\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Obtain timeseries data for parameters such as spin\n", "### Note that for many systems, the initial horizons are tracked all the way through merger after which the values are no longer reliable." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "time, primary_spin = primary_compact_object.dimensionless_spin_vector" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAEGCAYAAAC3lehYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAAA0vklEQVR4nO3deXyU1fX48c9JJvsK2QhL2GVfEiPuKALuirZurbVUsX7tYm1rW+3XVru4ddF+bWv1Z0XFahVFrdZaF1yq1pV9B5HNQAghQEISss75/XEnIWACCZnJM5Oc9+v1vGZ7ZuZclpzce89zr6gqxhhjTChEeR2AMcaY7suSjDHGmJCxJGOMMSZkLMkYY4wJGUsyxhhjQsbndQBdITMzUwcNGuR1GMYYEzEWLly4U1WzOvs5PSLJDBo0iAULFngdhjHGRAwR2RyMz7HhMmOMMSFjScYYY0zIWJIxxhgTMj1iTsYYYzqjvr6eoqIiampqvA4l6OLj4+nfvz8xMTEh+XxLMsYYcxhFRUWkpKQwaNAgRMTrcIJGVSkrK6OoqIjBgweH5DtsuMwYYw6jpqaGjIyMbpVgAESEjIyMkPbQLMkYY0w7dLcE0yTU7bLhskN59x7QRvAlQEqfwJHrjthEr6MzxpiwZ0nmUP7zW2jY1/pr8WmQ0tclntS+LvGk5kKvwZAxDNL6Q1R018ZrjDFhxpLMofzvVvA3Qn01VJZAxTbYW+yOiuL99z9b617Xxv3vjY6D3oGEkzHU3WaPgZzREJPgXZuMMRFr3rx5/P73v2ffvn2kpKTw/PPPk5XV6ZVfQsqSzKFERbvDFwsJ6ZA1ou1z/Y1QuQN2bYCy9YHjM9j5Kax7Ffz17jyJhqyRkDsB+k6EvOMhZyxE2fSYMebQpkyZwkUXXQTAL3/5S55++mm+853veBzVoVmSCZaoaDdclpoLg0488DV/I+zZDNtXQPFSd6x/HZb+3b0enwYDT3THiLNcz8cYE5Z++c+VrNpWEdTPHN03lVvPG3PY8x599FHmzp1LbW0t27dv54477mDDhg3cfvvtlJeXM2/evKDGFQz263NXiIqG3kNg9Pkw9efwtXnwo0/h+yvgwgdh1PlQugZeuxn+VACPXQCrX4LGBq8jN8aEiccee4yPP/6YN998k6VLlzJixAjGjBnDkCFDmD17dpvvKy0t5corr6SoqIirrrqK+vr6LozaejLeEYH0AZB+KUy41D2353NY+iQsfBTmXg6p/eDob0DB112BgTHGc+3pcYTC8uXLOeGEE0hOTubZZ5/l/fffZ9y4cYd9X1ZWFnl5edxwww3Mnj07ZFf2t8V6MuEkfQCc8hO4fhlc+oSbA3rrdvjDGHh6Jmx8B1S9jtIY44GZM2fyxz/+kZNPPpl169YxZMgQkpKSDvu+yspKNmzYgM/nIzk5uQsiPZAlmXAU7YNR58IVz8N1i+DYa2HD2zDnPLjvWPjwAdi3x+sojTFdaOzYsXz22We8++67/PSnP2XdunUAlJWVce2117J48WLuvPNOSkpKmofPGhoa+N73vsdtt93GxIkTefvtt7s8btEe8JtxYWGhRvymZfX7YOXz8MlDsHUhxCTCuIugcJarUjPGhMzq1asZNWqU12G0yyuvvEJsbCynnXZau9/TWvtEZKGqFnY2HpuTiRQxCTDxq+7Ythg+mQ3LnoFFj0G/QjjmahhzIcTEex2pMcZDZ555ptchHMCGyyJR33yY8We4YTWceRfUlMM/roV7RsJrP3PX6hhjTBiwJBPJEnrBcd+C734CX38RBk+GD/4Cf8yHv30J1vzLyqCNMZ4Ki+EyEekNzAUGAZuAS1R1dyvnPQycC+xQ1bFdGWNYE4Ehp7ijotgNoS18BJ76KqT2b1EGneN1pMaYHiZcejI3AW+o6nDgjcDj1jwKhNeAY7hJzYVTb3QXel76OGQOg7dugz+Mhme+AZveszJoY0yXCYueDDADODVwfw7wNnDjwSep6jsiMqjLoopk0T4YdZ47dq6HBQ/DkidchVrWSFeVNuFSt6SNMcaESLj0ZHJUtRggcJvd2Q8UkWtEZIGILCgtLe10gBEtcxiceQf8cDXMuM+VP//7x3D3KPjn9VC8zOsIjTHdVJf1ZERkPtDa2ig3h+L7VPVB4EFw18mE4jsiTmwi5H/NHVsXwYLZsPQpt4xN/0lwzCwYfYGVQRtjgqbLkoyqTmvrNREpEZFcVS0WkVxgR1fF1WP1K3DH6bfBkiddwnn+f+CVn7okVHilW9TTGGM6IVyGy14EZgbuzwRe8DCWniWhFxz/bfjuAvj6C26bgg/ugz8WwONfhrX/dlsVGGPMEQiXJHMXMF1EPgWmBx4jIn1F5OWmk0TkSeADYISIFInILE+i7Y5EYMipriLtByvglBuhZCU8eRncOwHe+b3blM0YYzrA1i4zbWush7UvuyVsNv4HomLcnjiFs2DgCS4xGdMDhMvaZaHaftnWLjPeiI6B0TPcsfNTVwa9+AlY8SxkjXKFAuMvhfhUryM1puv8+ybYvjy4n9lnHJx112FPi8Ttl8NluMyEu8zhcOadcMMaOP/P4IuDl38Ed4+Ef34/+P/pjDFf8OijjzJp0iQmTJjAX/7yF+Lj49mwYQOzZs1qTj7hxnoypmNiE6HgCndsXQifPBzYzfMRGHCsWw169AyXhIzpjtrR4wiFltsvJycnM3ny5AO2X24ryTzyyCP06dOHM888k1mzZnHfffeRkJDQZXFbT8YcuX5HwwX3uYs8T78dqnbCc9+Ee0bB67fC7k1eR2hMt3Gk2y9PnjyZ9957j9mzZ3PppZd2aYIBSzImGBJ7wwnfdWXQVzwPecfD+3+EeyfCExfD2lesDNqYTjrS7ZeHDh3K4sWLWbJkCWeccUYXRHogSzImeKKiYOhpcNkTboHOU37ilqx58lKXcN69Gyp7+BI/xhyhI9l+uYnP5+OWW27xImybkzEhktYPpvwvTP6x29dmwWx441fw1p1uzuaYWa7HY2XQxnRKRkYGDzzwQPPjV155hcGDBwNQXl7OzTffzMyZM8nO7vSSkEfErpMxXad0XWA16L9DbTlkj95fBh2X4nV0xrQpXK6TCZVQXidjw2Wm62Qd5SpzblgN5/8JonzwrxtcGfRLP3QrDBhjuhUbLjNdLzbJ7dSZ31QGPdvtdbNgthtCK5zlVhawMmhjIp71ZIx3RKB/IVx4f6AM+jaoLIHnroZ7RsP8X8DuzV5HaYzpBEsyJjwk9oYTroPvLoSvPQd5x8F/73WLcz5xCax71cqgjYlANlxmwktUFAyb6o7yIlg4BxbNgb9fAul5UHiVG2ZLyvQ6UmNMO1hPxoSvtP5w2s3wg5Vw8aOQPtANod0zCp69GrZ8CD2gOtKYSGY9GRP+omNgzIXu2LHGlUEvfRKWPwPZYwJl0JdYGbQxYch6MiayZI+Es3/rVoM+7143vPavH8Ldo1w5dMkqryM0xrRgScZEptgkOPob8D/vwqz5MOpcWPQ3uP94ePgsWD4PGuq8jtKYoJo3bx7HHXccEyZM4KSTTqK0NPyXabIkYyKbCAw4Bi58wJVBT/817N0Gz86CP4x2S9ns2eJ1lMYExZQpU/jwww9ZunQp06dP5+mnn/Y6pMOyORnTfSRlwInfg+O/CxvedBd5vvcHdww/3e11M3SqG2Iz5gj95uPfsGbXmqB+5sjeI7lx0o2HPe/RRx9l7ty51NbWsn37du644w42bNjA7bffTnl5OfPmzQtqXMFgScZ0P1FRMGyaO/Z8DgsfhUWPwbqLXIVacxl0hteRGtNuR7pp2ebNm/nrX//Kxo0bEREef/zxLo3bkozp3tIHwNSfwyk3wpp/up08598Kb93uqtUKZ8GASbYatGm39vQ4QuFINy0bOHAgs2bN4tZbb+X+++/vgkgPZEnG9Ay+WBj7ZXfsWB1YDfpJWDYXcsbBMVfBuEsgLtnrSI1p1cyZM5kxYwbz5s3j7LPPbvemZZs2beIXv/gF999/f7vODzYbnDY9T/YoOPt3rgz63P9zz730A7ca9L9+5JKQMWHmSDctO/vss+nduzd33nknu3bt6vK4rSdjeq64ZCi80pVCF30CnzzklrD55K8w8ER3kefI81wvyJgwdahNywBWrfL22jFLMsaIuHmZAZPgjDtg8eNuOG3eVZCU7bYlOPobbn7HmDB35plneh3CAcJiuExEeovI6yLyaeC2VyvnDBCRt0RktYisFJHrvYjVdHNJmXDS9+F7S+DyedDvaHjvHrh3PDz5Ffh0Pvj9XkdpTMQIiyQD3AS8oarDgTcCjw/WANygqqOA44DviMjoLozR9CRRUTB8Onz1Kbh+KZz0Azek9sSX4U8FbhuCqjKvozQm7IVLkpkBzAncnwNccPAJqlqsqosC9/cCq4F+XRWg6cHS82DqLfCDVfDl2ZDaF16/xa0G/dz/wOef2GrQxrQhXOZkclS1GFwyEZHsQ50sIoOAfOCjQ5xzDXANQF5eXvAiNT2XLxbGXeSOklVuu+ilc2HZU9BnnFtRYNzFbl01YwzQhT0ZEZkvIitaOWZ08HOSgWeB76tqRVvnqeqDqlqoqoVZWVmdDd+YA+WMhnPuhhtWwzn3uJ7MP693ZdAv/wRK13odoQky7aa91VC3q8t6Mqo6ra3XRKRERHIDvZhcYEcb58XgEswTqvpciEI1pv3iUlypc+FV8PnHrgx64SPw8f+DQSe750eea2XQES4+Pp6ysjIyMjKQbrQ6hKpSVlZGfHx8yL5DwiE7i8jvgDJVvUtEbgJ6q+pPDjpHcPM1u1T1+x35/MLCQl2wYEHQ4jXmkKp2wuK/uTLoPVsgOWd/GXRaf6+jM0egvr6eoqIiampqvA4l6OLj4+nfvz8xMTEHPC8iC1W1sLOfHy5JJgN4GsgDtgAXq+ouEekLPKSqZ4vIScC7wHKgqYb0f1X15cN9viUZ4wl/I6x/w83drHvVXY9z1Fmu5zNkiq0GbcJat0oyoWZJxnhu9+b9q0FX74Reg12ymXg5JPb2OjpjvsCSTAdYkjFho6EWVv/Tzd1s+QCi49yincfMchd+dqPxfhPZgpVkwqWE2ZiewRfXogx6pdtYbdlcWPp3yJ3gth4Yd5GVQZtuw3oyxnitdq9LNJ88DDtWQlwaTPyKSzhZR3kdnemhbLisAyzJmIigCls+dIUCK/8B/npXBn3M1TDyHIiOOexHGBMslmQ6wJKMiTiVpYEy6EegfAsk94GjZ0LBTEiz1ZRM6FmS6QBLMiZi+Rth/XxXKPDp6yBRMCJQBj34VCuDNiFjE//G9ARR0XDUGe7Yvcn1bBb/Dda8BJkj4ITrYPwlrqDAmIBGv7K6uIKi3fs4c2wfT2OxnowxkaahFlY+D+//GUqWu6G0466Fo6+EhHSvozMeqWvw8976Ul5aVsz8VSVU1DSQEudjya2nEx3V8dJ4Gy7rAEsypltShQ1vub1tNrwNsclu6ZrjvmXL1/QQ9Y1+/rt+J/9aVsyrK7dTUdNAaryP6aP7MPmoTI4dnEGftCNbl6zLh8tE5GLgFVXdKyI/AwqA25r2eDHGdDERGHqaO4qXwvt/gg/vh48egLEXwYnfg5wxXkdpgqyh0c/7n5W5xLJqO3uq60mJ8zF9TA7njs/lpGFZxPrCZ66u3T0ZEVmmquMDa4jdCfwet3bYsaEMMBisJ2N6jD1b4IO/uOVr6qtg+Olw4vUw8ERbTSCC1dQ38v5nO3ltZQmvrtzO7up6kmKjmT46h3PG9+Xk4ZnEx0QH9Tu7fLhMRBarar6I3AksV9W/Nz3X2SBCzZKM6XGqd7nVBD56wK2V1q8QTvo+jDjHKtIiQEVNPSuKyllStIfFW/bw3/U7qa5rJDnOx5SR2ZwzLpdTR2QFPbG05EWSeQnYCkwDjgb2AR+r6oTOBhFqlmRMj1W/D5Y84YbSdm+CrJFw6k0waoYlmzBRU9/I6uIKlhWVs/TzPSwt2sNnpVXNrw/MSOSkYZlMH53D8UMziPOFLrG05EWSSQTOxPViPg1sLjZOVV/rbBChZknG9HiNDbDqH/Cf38LOtZAzFk79qVtJwIbRuoyqsr2ihkWb97Bw824Wbt7FquIK6hvdz+GslDgm9E9nQv80JgxIZ3z/NNITvdnwzosk8xtVvfFwz4UjSzLGBPgbYcWz8PZdsOszyJ0Ip/0chk21ZBNkqsq28hqWF5WzYms5y7eWs3JbOTsr6wCIj4liQv908vN6MXGASyp9UuPDZudNL5LMIlUtOOi5Zao6vrNBhJolGWMO0tjgFuX8z29gz2a3Rtq0X0D/Tv9M6ZFUla179jUnk+VbK1ixtZxdVS6hREcJw7OTGdcvjbH90sjPS2dUbiox0eE7ZNllJcwi8i3g28AQEVnW4qUU4P3OBmCM8UC0D/Ivh3EXw8JH3DDaQ1Nh1Hlw2i22+vMhqCpFu/cFkonrpazYWs7u6noAfFHC8JwUpo3Kbk4qo3JTQzpJH84O25MRkTSgF65s+aYWL+1V1V0hjC1orCdjzGHU7oUP7nMFAvX7XAI69aeQ2tfryDylqny+66CEsq2cPS0Syog+KYztm8bY/mmM65fGyD4p3SKheDFcFg98CRhEix6Qqv6qs0GEmiUZY9qpshTe/b0rf46KhmOvdaXPCb28jizkVJXNZdWs2Nayh1JB+T6XUGKiXUJp6p2M65fGiD4pXVbt1dW8SDKvAOXAQqCx6XlVvbuzQYSaJRljOmj3JnjrDlj2NMSnweQfw6RvdpuFOP1+ZfOu6uZksrzI9VD21jQAEBsd5XoogWQyrl8aR/VJ7rYJpTVeJJkVqjq2s1/oBUsyxhyh7cvh9VvgszchfaArDhhzYURVovn9yqayqv0JZWs5K7dWsLd2f0IZlZvCmJYJJSclrJZm8YIXS/2/LyLjVHV5Z7/UGBMh+oyDK553e9q8dgvMuxI+/Aucfjvkhd+KUn6/smFnVYsqr3JWbaugsimh+KIYlZvKjPy+zcNeR+WkhHWVV6TrSE9mFTAc2ADUAgKolTAb00P4G2HJ3+HN26ByO4w63/VsMoZ6Ek6jX9m4s9IlkyJXMrxyWzlVdW40Py6QUJp6J2P7pTE8J9kSSjt5MVw2sLXnVXVzZ4MINUsyxgRRXZXby+a/90JjHRxzNZzyE0jsHbKvbPQrn5VWNvdQXEKpoDqQUOJjohgdSChjmxJKdjI+SyhHzIskI8DlwBBV/ZWI5AF9VPXjTgch0huYi6tc2wRcoqq7DzonHngHiMMN881T1Vvb8/mWZIwJgb3bXXHA4r9BbApM/hFMugZijmz/kiYNjX4+Kz1wDmXVtgr21buEkhATzei+qQdUeQ3NSrKEEmReJJn7AT9wmqqOEpFewGuqekyngxD5LbBLVe8SkZuAXq0sYSNAkqpWikgM8B5wvap+eLjPtyRjTAiVrHLFAetfh/Q8mHorjP1yu4oDGhr9rC+tPGDplVXFFdTU+wGXUMb0Td1f5dU/jaFZyUe006PpGC8m/o9V1QIRWQygqrtFJFgrt80ATg3cnwO8DRyQZNRlw8rAw5jA0f239TQm3OWMhq/NcxVor/0cnp0VKA64DQae0HxafaOfT0taDHltK2d1i4SSFBvNmL5pfHXSQMb1T2Vs3zSGWEKJeB1JMvUiEk3gB7uIZOF6NsGQo6rFAKpaLCLZrZ0U+P6FwDDgPlX9qK0PFJFrgGsA8vLyghSmMaZNQ0+D/zkFlj6Fvvlr5JGzKMqZyjO9ZvF2WTpriiuobXA/MpLjfIzum8rlxw5sHvYanJlkCaUb6shw2eXApbhtl+cAFwE/U9Vn2vn++UCfVl66GZijquktzt2tqm1eYiwi6cDzwHWquuJw323DZcaETl2Dn3Ule5tLhlduLWfj9p18XV/iWt8/iaOe+UnnsGr4txg2eKBLKBlJRFlCCWtdPlymqk+IyEJgKq58+QJVXd2B909r6zURKRGR3EAvJhfYcZjP2iMib+P2tzlskjHGBEdtQyPrtlcesJbX2u17qWt0PZSUeB9j+6Zx2QkjOKrfJMp6/Yy8pfdy1qJHOWvN25D1Q0j/FliC6THa3ZMJaRAivwPKWkz891bVnxx0ThZQH0gwCcBrwG9U9aXDfb71ZIzpuJr6RtZu38uKbfsn5ddu39u8wVZqvK95Qr7pNq93Yus9lB1rXHHAp69C2gCYeguMvch25wxjXlSXxQFfJgQLZIpIBvA0kAdsAS5W1V0i0hd4SFXPFpHxuGG6aCAKeLq9321JxpjW+f1KWVUd2/bsY+POKtbvqHRHaSUbd1bR6Hc/H9ISYhjbL/WAtbzyeid2fIOtDf+B134G25e5DdPOuB0GnRT8hplOswUyO8CSjOlpVJU91fUUl9dQWllL6d5aduytCdzWUlJew/aKGkoqapp7JuA21xrYO5Gh2cmMyElpLh/u3ysheDs2+v2w/Gl441dQsRVGnA3Tfml72IQZWyCzAyzJmO5qX10j63dUsrZkL5+W7GVzWTVbdlXz+a7q5gUgW0qO85GVEkdOahy5aQn0SYunT2o8fdLiGZKZxMCMpK5bGLJ+nyt1fvcPUF8NhVfCKTdBclbXfL85JFsg05huTFWpqmtk595adlY2HXXN90sqal1S2VVN0++Jsb4o8nonktc7kUmDezOgdyK5afFkp8SRnRJPZkosibEd+S8fYjEJcPINkP91+M9dsOARWDoXTv4BHPdt97qJeB1dIHMYsBFbINOYDvP7lfJ99eysrKU0kDTKmhLI3kACqaprTixN15QcrFdiDFkpcQzLTuaonBRG9knhqJwUBmZE+HUmpetg/q2w9mVI7Q9Tfw7jLrHiAI/YApkdYEnGhEpDo59dVXWUVtZS1qKnsb/XsT9p7Kqqo8H/xf9v0VFCRlIsGclxZCbHkpUcR2ZKHBlJsWQG7mcmu/u9k2K7/yrCG991xQHFS6DPeLdywJBTvI6qx+nyJBP40gnAyYGH76rq0s4G0BUsyZiOqKlvZGc7kkZZVR27q+to7b9QrC/KJYtAcshMjiMj+YtJIzM5jvSEGLsw8WB+P6x4Ft74JZR/DsPPgOm/guyRXkfWY3jRk7ke+CbwXOCpC4EHVfVPnQ0i1CzJ9GyqSmVtwwFJo7SVoaqywFBVaxPm4CbN20oaWclNPRGXQJLjfMGrxurJ6mvgowfg3buhrhIKZsKU/4XkVleeMkHkRZJZBhyvqlWBx0nABzYnY7wQzPmNAxJGchxZbQxVxcf0nP3dw07VTvjPb2HBbPDFw4nfh+O/A7GJXkfWbXlRXSa0uD4mcN9+VTNB03J+44Ck0TREVdXx+Y2hmUk9e36ju0jKhLN/6/armX8rvHUbLHgYTvsZTLgMouwXgHDVkSTzCPCRiDyPSy4zgNkhicp0G03zGwcnjdK9+4enmoawdlfXt/oZLec3ctPiGdcvzeY3eqrMYXDZE7D5fXj1Znjh2/Dh/XD6r2HoFK+jM63o6MR/AdC0BsS7qro4JFEFmQ2XBU/T/EbLpFEa6GmUVbUYqgpMnNv8hgkZvx9WPgfzfwnlW2DYdFcckDPa68i6hS4fLgtsf3wqrrrMD0SLyGpVrelsEMZbwZ7fyEyOY2y/tOb5jczkWDKSbH7DBFlUFIy7CEaeCx8/CO/8Hh44EfKvcMUBKa3tLGK6Wkcm/p8G9gKPB576Cm6b5ItDFFvQ9MSeTCjmN9q6fiMrOY5eNr9hvFa9yxUHfPIQRMfCidfDCd+F2CSvI4tIXlSXLVXVCYd7Lhx1lyQT7PkNu37DdEtln8H8X8DqFyG5jysOmPhVKw7oIC+qyxaLyHGq+mEggGOB/3Y2gO5IValt8FNV20B1XSO1DY3UNSh1jX7qGvzUB27rGt19v7r3qMLe2oZOzW8MzUpm0uDeNr9heq6MoXDp32DLh6444MXv7i8OGDbV6+h6nI70ZFYDI3D7vYDb+2U1bn4mrNcw60xPZmdlLcuK9rClrJqyqjr21TVS1+hnX10j1XWNVNU1UFnTwN6aBiprG9hbU091XWOrw08d0XJ+o6mnYfMbxnSQKqz6B7x+K+zZDENPg+m/hj4RuaB8l/KiJ3NmZ78skvj9ytWPLeDDDWVU17nLg6IE4nzRxMVEkRATTWJsNImxPpLjfAzMSCQ53kdKnI+kwJEc5yMxNpr4mGhioqOI80UREx1FrC+KmGgJ3EYRJSAiCJAU57PrN4wJFhEYc6Hbs+aTh9yczQMnQf7lMOVnkJrrdYTdXruTTCQshBlMUVFCXYOf00fncNmkPIZkJZGVHGdDTcZEIl+cWyFgwlfcEjUf/T9Y8RyccB2c8D2IS/Y6wm6rI8NlhcDNwEBccrKl/o0xkWnXBrcz58rnITnHlTxP/BpEh9F+Ox7zorpsLfBjYDluHgaIjB6OJRljTKs+/9htK/D5R5A1KlAcMM0Ns/VwwUoyHRn4L1XVF1V1o6pubjo6G4AxxnhmwCS46lW45DFoqIEnLoK/XQDFy7yOrNvoSE9mKu4CzDdwO2MCoKrPtfmmMGE9GWPMYTXUuVWe//Mb2LfHXVsz5WZI6+d1ZJ7worrsSmAkEMP+4TJl//4yxhgTuXyxcNy33KrOBxQHfNetHhCX4nWEEakjPZnlqjouxPGEhPVkjDEdtnsTvPFrWDEPkrJccUD+13tMcYAXczIfiogtb2qM6Rl6DYKLZsPVb0LGcHjpB3D/CbDuVVrdc9u0qiNJ5iRgiYisFZFlIrI8sFumMcZ0X/2PhitfhkufAH8D/P0SmHMebFvidWQRISyu+BeR3sBcYBCwCbhEVXe3cW40sADYqqrnhiomY4xpJgKjzoWjzoAFj8Dbd8KDp8D4y2DqzyGtv9cRhq1292QC5crpwHmBIz2IJcw3AW+o6nBc9dpNhzj3etyaacYY07WiY+DYa+D6JXDi993FnH862m2cVlPhdXRhqd1JRkSuB54AsgPH4yJyXZDimAHMCdyfA1zQRgz9gXOAh4L0vcYY03HxaTD9l3DdAhh1Prx3D/wxHz7+KzS2vs1GT9WR6rJlwPGqWhV4nAR8EIxlZURkj6qmt3i8W1V7tXLePOBOIAX40aGGy0TkGuAagLy8vKM3b7brRo0xIbJ1Ebz2c9j8nisSmP4rGHFWRK8c4EV1mQCNLR43Bp5r35tF5ovIilaOGe18/7nADlVd2J7zVfVBVS1U1cKsrKz2hmmMMR3XrwC+8RJc9qR7/NRX4NFzXfLp4Toy8f8I8JGIPB94fAHwcHvfrKrT2npNREpEJFdVi0UkF9jRymknAueLyNlAPJAqIo+r6tfa3QJjjAkVERh5NgyfDovmwFt3wl+nwLhLXHFAep7XEXqi3cNlACJSgCtlFuAdVV0clCBEfgeUqepdInIT0FtVf3KI80/lMMNlLdnFmMaYLldTAf/9P/jgPnddzXHfgpN/6OZzIkCXD5eJyBxgg6r+UVXvBTaJSLt7ModxFzBdRD4FpgceIyJ9ReTlIH2HMcZ0nfhUmHoLXLcQxn7JJZx7J7rlanpQcUBHJv4Xq2r+4Z4LR9aTMcZ4btsSt63Apneh91BXnTby3LAtDvBi4j9KRJorvgIXUPaMRXyMMaaz+k6Emf+Erz4NUT6Y+zV45CwoalctU8TqSJK5G3hfRH4tIr8C3gd+G5qwjDGmGxJxqwZ863049w9Qth4eOg3mXeUW5OyGOnLF/2PAl4ESoBT4kqr+LVSBGWNMtxXtg8Kr4HuLYfKPYc3LcN+x8N97obHB6+iCqkPVZZHK5mSMMWGtfCv8+yew5iXoMx5m/BlyJ3gakhdzMsYYY0IhrR9c+rjbBrqyBB6cAq/fAvX7vI6s0yzJGGNMOBCB0TPgOx9B/uVu6OyBk2DLR15H1imWZIwxJpwk9ILz/wRffwEa6uDhM+DVmyO2V9PpJCMiNwYjEGOMMS0MORW+/T4UXgkf/Nn1aj7/2OuoOqzDSUZEnm5xPANcHYK4jDHGxKW4UueWvZrXfhZRvZoj6clUqOolgeNiYH6wgzLGGNNCU6+mYCa8/yd44OSI6dUcSZK5/aDHNwcjEGOMMYcQlwLn/R9c8Q9oqImYXk1HFsjMEJFvAaeKyCQRSQBQ1V0hi84YY8yBhk6Bb39wUK/mE6+jalNHejLPA1nAHcDvgHIRWROSqIwxxrTtC72a0+HfN0JNudeRfUFHkkyKqv4KKFHVU4Cv4DYyM8YY44WhU9w6aIVXuS0E/nwMLHvG7V8TJjqSZGoCt7UikqCqzwJnhyAmY4wx7RWfCufcDd98E1L7wXNXw5zzoHSt15EBHUsyvw8s7z8XeFhErgP6hSYsY4wxHdKvAK6eD+fcA9uXw/0nBJamqTn8e0OoI6swP6uqu1T1HuBlYAAwI2SRGWOM6ZioaDhmltuNc/xlsOFtiI7xNKQj2nTMlvg3xpgwlpQJF9znypujoj0NxdYuM8aY7iomwesILMkYY4wJHUsyxhhjQsaSjDHGmJCxJGOMMSZkLMkYY4wJmSMqYQ62Fhd5DgI2AZeo6u5WztsE7AUagQZVLey6KI0xxnRUuPRkbgLeUNXhwBuBx22ZoqoTLcEYY0z4C5ckMwOYE7g/B7jAu1CMMcYES7gkmRxVLQYI3Ga3cZ4Cr4nIQhG55lAfKCLXiMgCEVlQWloa5HCNMca0R5fNyYjIfKBPKy91ZGfNE1V1m4hkA6+LyBpVfae1E1X1QeBBgMLCwvBZ99oYY3qQLksyqjqtrddEpEREclW1WERygR1tfMa2wO0OEXkemAS0mmSMMcZ4L1yGy14EZgbuzwReOPgEEUkSkZSm+8DpwIoui9AYY0yHhUuSuQuYLiKfAtMDjxGRviLycuCcHOA9EVkKfAz8S1Vf8SRaY4wx7RIW18moahkwtZXntxHYfVNVNwATujg0Y4yJWNX11STGJHoaQ1gkGWOMMcHR6G/knaJ3eGLNE5RWl/KPGf9ARDyLx5KMMcZ0A1srt/LC+hd4Yf0LbKvaRk5iDpeNvIx6fz2x0bGexWVJxhhjIlRNQw1vbnmT59Y/x0fFHyEIx+Uexw8Lf8hpeacRE+Xt1stgScYYYyLO2l1reWbdM7y88WX21u2lX3I/vj3x28wYOoO+yX29Du8AlmSMMSYC1DXW8frm15m7di6LdywmNiqWaQOn8aXhX+KYPscQJeFSLHwgSzLGGBPGiiuLeWbdMzz76bPsqtlFXkoePyr8ERcMu4C0uDSvwzssSzLGGBNmVJWPtn/Ek6uf5O2itwGY3H8yl424jOP7Hh+2vZbWWJIxxpgwUVVfxYufvciTa55kY/lGesX14soxV3LJiEvCbq6lvSzJGGOMh1SVJaVLeGH9C7yy6RWq6qsYlzmO20+6nTMGnUFcdJzXIXaKJRljjOliqsq63euYv2U+/974bzZXbCbBl8D0gdO5bMRljMsa53WIQWNJxhhjusDOfTtZsmMJC0sW8k7RO2zZuwVBKOxTyNXjrmb6wOkkxSR5HWbQWZIxxpggU1U2VWxi8Y7Fzcfmis0AxEXHcXTO0Xxj7DeYMmAKmQmZHkcbWpZkjDGmk+ob61m1axWLS/Ynld21uwFIj0snPzufi4ZfRH5OPqN7jyYm2vsr8buKJRljjOmgiroKlu5Y2pxQlu9cTm1jLQB5KXlM7j+ZgpwCJmZPZHDqYE8XqPSaJRljjDmM4spiFu1Y1JxUPt39KYoSLdGM6j2KS0ZcQn52PvnZ+d1++KujLMkYY0wLjf5G1u9Zz+Idi5sTy/aq7QAk+hKZmD2RaQOnUZBdwLjMcZ7v1xLuLMkYY3q0fQ37WLFzRXNSWbpjKZX1lQBkJ2S7Ya8xEynILmB4r+H4ouzHZkfYn5YxpkfZVbPLDXsFJulXla2iQRsAGJY+jLMGn0V+dj4FOQX0Terbo+dTgsGSjDGm21JVtuzdwqKS/fMpmyo2ARAbFcvYzLHMHDOTgpwCJmRNiIgFJyONJRljTLdR769nTdmaAybpd9XsAiAtLo38rHwuHH4hBdkFjM4Y7emOkT2FJRljTMSqrKtkaelSFu1YxJIdS1hWuoyaxhoA+if356R+J7mhr+wCBqUNiqjVi7sLSzLGmIixvWq7m6AvWcSS0iWs270Ov/qJkihG9h7JRUdd1FxKnJWY5XW4Bksyxpgw5Ve/KyUuWcziUjdRv61qGwAJvgQmZE3g2vHXMjF7IuOzxnfLdb+6A0syxpiwUNNQw4qdK1hSuqS5p7K3bi8AmQmZFGQXcMXoK8jPyWdErxFWShwh7G/JGOOJ3TW7WbJjSfP1KSvLVtLgd6XEQ9OGcvrA0ynIKSA/O5/+yf2tlDhChUWSEZHewFxgELAJuERVd7dyXjrwEDAWUOAqVf2gywI1xhwRVaVob1Fz1deiHYvYWL4RgJioGMZkjOGK0VdQkF3AxKyJpMenexuwCZqwSDLATcAbqnqXiNwUeHxjK+fdC7yiqheJSCxg6zkYE4Ya/A2s3bV2f1IpWURZTRkAqbGpTMyeyPlDz6cgu4AxmWMifvdH07ZwSTIzgFMD9+cAb3NQkhGRVGAy8A0AVa0D6roqQGNM26rqq1haurT5SvplO5exr2EfAP2S+3F83+ObS4mHpA+xUuIeJFySTI6qFgOoarGIZLdyzhCgFHhERCYAC4HrVbWqtQ8UkWuAawDy8vJCE7UxPVRJVUlzxdfiHYtZu3ttcynxiF4juHDYheTn5JOflU9OUo7X4RoPiap2zReJzAf6tPLSzcAcVU1vce5uVe110PsLgQ+BE1X1IxG5F6hQ1Z8f7rsLCwt1wYIFnYrfmJ7Kr3427NlwwFX0Wyu3Aq6UeHzmeJdQsvMZnzme5NhkjyM2wSAiC1W1sLOf02U9GVWd1tZrIlIiIrmBXkwusKOV04qAIlX9KPB4Hm7uxhgTBH71U1xVzIY9G9hYvpEN5e52/Z71VNRVAJARn0FBTgGXj7qc/Ox8RvQeQUxUz9nl0XRcuAyXvQjMBO4K3L5w8Amqul1EPheREaq6FpgKrOraMI2JfLWNtWwq38TGio1s3LOxOaFsrtjcvCQLuG2DB6cNZvrA6UzMdkvdD0gZYKXEpkPCJcncBTwtIrOALcDFACLSF3hIVc8OnHcd8ESgsmwDcKUXwRoTzqrrq9levZ3iymKKq4rZVrmN7VXb2Va1rfk5xQ2TC0Lf5L4MThvMpNxJDEkbwuC0wQxJG0Kv+F6H+SZjDq/L5mS8ZHMyJpLVNdZRVV91wFFZX0lVfRU7qnewvWr7Aclkd+2Bl5hFSRQ5iTnkJuXSJ6kPeal5zclkYOpAEnwJHrXMhLOIm5MxpqdQVWoba6mqr6K6vro5IVQ3VFNZV0lVQxVVdVXutv7wR72//pDfl+BLoG9SX3KTcxmbOZbcpFxyk3PJTcqlb1JfshKzbAkW4xn7l2cMLjHsa9hHdUN1c0+huv7A+83Joo37TUd1fXXzTouHk+hLJCkmqflIjkkmPTn9gOfaPHxJZCVmkRqbavMkJmxZkjERy69+9jXsO2QvoLVkccA5LXoVfvUf9jsFafUHfmZCZruSQlJs4DYmicSYRLso0XR7lmRMl/Krv/m3/zZ7B03DSm0kg5Y9hqYJ7EOJlmgSYxJJjklu/uGeHJtMTlIOiT53v+k2yZd0wLkHHwm+BOs1GNMBlmTMYTX4G9wwUt1BvYOGKirrKr8wxHRw4mjZk2haauRwfFG+5uGjxJhEknxJpMWn0dfXtzkpHPB6G/eTY5KJi46zxGCMRyzJdFP1/vpW5wuO5Gh57cShxEbFfuE3/4z4DPJS8g47lHRwgrC9143pHizJtEODv4GS6hKiiGr+jThaoomOiiZaoomJiiHeF3/E4+uqSoO/gXp/PbWNtZ2bdA4MNdX527d2aHx0/Bd+4GcnZn9xeKmN4aPmuYaYJGKi7cpvY8yBLMkcxpyVc7h7wd3tGvv3RfmIiYrBJz6io6KJkigEab5VlEZtpN5fT6O/kQZ/Aw3a0K4J55YSfAlfSAC5yblfGCZKiklqc66hqddgpa3GmFCynzBtqPfXc9lLl7Fu9zp6x/fmqrFXuWqgQG+mwd9Ag7+hOWnUNNRQ11jnEog20uhvpFEbURRVRVEEwRflI1qi8UX59h+y/35cdNwh5xcSfAlER0V7/cdjjDHtYkmmDTFRMYzsPZJpA6dx9birbRFAY4w5ApZkDuH2k273OgRjjIlodiWYMcaYkLEkY4wxJmQsyRhjjAkZSzLGGGNCxpKMMcaYkLEkY4wxJmQsyRhjjAkZSzLGGGNCRlQPvyZXpBORUmDzEb49E9gZxHDCgbUp/HW39oC1KVI0tWmgqmZ19sN6RJLpDBFZoKqFXscRTNam8Nfd2gPWpkgR7DbZcJkxxpiQsSRjjDEmZCzJHN6DXgcQAtam8Nfd2gPWpkgR1DbZnIwxxpiQsZ6MMcaYkLEkY4wxJmQsybRBRM4UkbUisl5EbvI6nvYSkQEi8paIrBaRlSJyfeD53iLyuoh8Grjt1eI9Pw20c62InOFd9G0TkWgRWSwiLwUeR3R7AEQkXUTmiciawN/X8ZHcLhH5QeDf3AoReVJE4iOxPSLysIjsEJEVLZ7rcDtE5GgRWR547Y8iIl3dlkAcrbXnd4F/d8tE5HkRSW/xWnDbo6p2HHQA0cBnwBAgFlgKjPY6rnbGngsUBO6nAOuA0cBvgZsCz98E/CZwf3SgfXHA4EC7o71uRyvt+iHwd+ClwOOIbk8g1jnA1YH7sUB6pLYL6AdsBBICj58GvhGJ7QEmAwXAihbPdbgdwMfA8YAA/wbOCqP2nA74Avd/E8r2WE+mdZOA9aq6QVXrgKeAGR7H1C6qWqyqiwL39wKrcT8AZuB+qBG4vSBwfwbwlKrWqupGYD2u/WFDRPoD5wAPtXg6YtsDICKpuP/8swFUtU5V9xDZ7fIBCSLiAxKBbURge1T1HWDXQU93qB0ikgukquoH6n5CP9biPV2qtfao6muq2hB4+CHQP3A/6O2xJNO6fsDnLR4XBZ6LKCIyCMgHPgJyVLUYXCICsgOnRUJb/w/4CeBv8VwktwdcL7kUeCQwDPiQiCQRoe1S1a3A74EtQDFQrqqvEaHtaUVH29EvcP/g58PRVbieCYSgPZZkWtfaWGNE1XqLSDLwLPB9Va041KmtPBc2bRWRc4EdqrqwvW9p5bmwaU8LPtwQxv2qmg9U4YZh2hLW7QrMUczADbH0BZJE5GuHeksrz4VNezqgrXZERPtE5GagAXii6alWTutUeyzJtK4IGNDicX9c1z8iiEgMLsE8oarPBZ4uCXR5CdzuCDwf7m09EThfRDbhhi1PE5HHidz2NCkCilT1o8DjebikE6ntmgZsVNVSVa0HngNOIHLbc7COtqOI/UNQLZ8PGyIyEzgXuDwwBAYhaI8lmdZ9AgwXkcEiEgtcBrzocUztEqj4mA2sVtV7Wrz0IjAzcH8m8EKL5y8TkTgRGQwMx03whQVV/amq9lfVQbi/hzdV9WtEaHuaqOp24HMRGRF4aiqwisht1xbgOBFJDPwbnIqbD4zU9hysQ+0IDKntFZHjAn8eX2/xHs+JyJnAjcD5qlrd4qXgt8eLaodIOICzcZVZnwE3ex1PB+I+CdeNXQYsCRxnAxnAG8CngdveLd5zc6Cda/GoAqadbTuV/dVl3aE9E4EFgb+rfwC9IrldwC+BNcAK4G+4CqWIaw/wJG5eqR73G/ysI2kHUBj4s/gM+DOBFVbCpD3rcXMvTT8jHghVe2xZGWOMMSFjw2XGGGNCxpKMMcaYkLEkY4wxJmQsyRhjjAkZSzLGGGNCxpKMMSEUWGn52wc994qI9BORt0VkS8vVbEXkHyJS2fWRGhMalmSMCa10oDnJiEgC7hqLrYGn9uBWNSCw3Hpu14ZnTGhZkjEmtO4ChorIEhH5He6C0rdbvP4UbiUDgC/hlmMxptuwJGNMaN0EfKaqE1X1x8BZwCstXn8DmCwi0bhkM9eDGI0JGUsyxnStE4H3WjxuDDy+FLfh1yYvgjImVCzJGNNFRGQI8Lm6jfBaegr4E243SWO6FUsyxoTWXtw22PDFobIm7wJ34hYyNKZb8XkdgDHdmaqWich/RWQFbiXcb7ZyjuJ2lTSm27FVmI3pAiISB/xXVQu9jsWYrmRJxhhjTMjYnIwxxpiQsSRjjDEmZCzJGGOMCRlLMsYYY0LGkowxxpiQsSRjjDEmZP4/6s3cVs6XrHMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(time, primary_spin[:, 0], label=r'$a_{1, x}$')\n", "plt.plot(time, primary_spin[:, 1], label=r'$a_{1, y}$')\n", "plt.plot(time, primary_spin[:, 2], label=r'$a_{1, z}$')\n", "plt.legend()\n", "plt.xlabel('t/M')\n", "plt.ylabel(f'$a_{1}$ components')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "time, secondary_spin = secondary_compact_object.dimensionless_spin_vector" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEGCAYAAABCa2PoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAAAw4ElEQVR4nO3deZxddZnv+8+zh5qrUkmqiiQVQhKIYZDRNKAiHkQw0LbBvj3osVtErlxP4wv7HLWFSx+HPrbiaYdubYWLEgGHBi4IcgVRRLFRRAkzhJCEBEhlqspUVUlNe3juH2vtql2VXaF2slftXVXf9+u1Xmte+9kh7Ce/Yf1+5u6IiIgcqVi5AxARkelBCUVEREpCCUVEREpCCUVEREpCCUVEREoiUe4AyqmlpcUXL15c7jBERKaUJ554Ype7t449PqMTyuLFi1mzZk25wxARmVLM7NVCx1XlJSIiJaGEIiIiJaGEIiIiJTGj21AKSaVSdHR0MDAwUO5QSq6mpoaFCxeSTCbLHYqITEMVlVDMbCXwb0Ac+K67Xzfm/PHA94AzgGvd/St5514BeoEMkHb3FYcTQ0dHB42NjSxevBgzO7wvUoHcnd27d9PR0cGSJUvKHY6ITEMVk1DMLA58C7gA6AAeN7N73X1t3mV7gKuAS8Z5zHnuvutI4hgYGJh2yQTAzJg7dy5dXV3lDkVEpqlKakM5E9jo7pvcfQi4DViVf4G7d7r740AqykCmWzLJma7fS0QqQyUllHZgS95+R3hsohz4hZk9YWZXjHeRmV1hZmvMbI3+tS4iM822ff189RcvsXnXgZI/u5ISSqF/PhczWctb3f0M4CLgSjM7t9BF7n6ju69w9xWtrQe96CkiMq1t7+7nm7/ayJY9fSV/diUllA7g6Lz9hcC2id7s7tvCdSdwN0EV2pR15513cvbZZ3PqqadyzjnnqO1DREoilQn+nZ6Ilb4KvJISyuPAMjNbYmZVwPuAeydyo5nVm1ljbhu4EHg+skgnwXnnncdjjz3GM888wwUXXMAdd9xR7pBEZBrIZIOEEo8goVRMLy93T5vZx4CfE3QbXu3uL5jZR8PzN5jZPGAN0ARkzezvgROBFuDusNE5AfzI3R840pg+//+9wNptPUf6mFFOXNDEZ//spNe97uabb+b2229ncHCQHTt28MUvfpF77rmH++67j87OTq688kouvPDCksYmItNfOkwoiXjpyxOVVELB3e939ze4+7Hu/s/hsRvc/YZwe4e7L3T3JndvDrd7wp5hp4bLSbl7p6pbb72VP/7xj/zqV7/imWeeYfny5Zx00klccsklfOc73xlONmN1dXVx2WWX0dHRwYc//GFSqUg7w4nIFJTOZIFoqrwqpoRSiSZSkojCc889x1ve8hYaGhq46667ePTRRzn55JOHz3/hC1/gyiuvPOi+1tZWFi1axCc+8QluuukmvREvIgdJR1jlVVElFAlceumlfOMb3+Btb3sb69evZ+nSpdTX1+PufPrTn+aiiy7ijDPOOOi+/fv3s2nTJhKJBA0NDWWIXEQqXa4NJTndq7wk8MY3vpGXX36ZRx55hGuuuYb169cD8M1vfpNf/vKX3Hnnndxwww3s3LmTm266CYB0Os1VV13FF77wBU477TQefvjhMn4DEalUqbDKa1o3ysvru+qqq7jqqquG9x944IHhcbkSiQSrV68G4FOf+lRZ4hORyjdSQlFCkTwrV64sdwgiMsWoDUVEREoiPfxio9pQRETkCGSyYbfhCKq8lFBERGaQmTL0ioiIRCzKoVeUUEREZpC03kMREZFSSEf4HooSiojIDDI8OKQSioiIHIlM1onHLJIpwZVQRERmkFQ2G0l1FyihiIjMKJmMk1RCmVk0BbCIRCEdVnlFQWN5HcrProYdz5X2mfNOhouue93LzjvvPP7iL/4CgM9//vPccccdBedAEREpRjqbjWS2RlAJpWLdfPPNnHnmmZx66ql8+9vfpqamhnvuuYePfOQjrFq1il/84hflDlFEpqBM1iPp4QUVVkIxs5XAvxHMKf9dd79uzPnjge8BZwDXuvtXJnrvYZlASSIK+VMANzQ0cO6553LSSSdx9tlnc8kll7B3714++clPHjSn/Pe+9z3mzZvHypUrufzyy/nWt75FbW1tWb6DiFSmVGYGJBQziwPfAi4AOoDHzexed1+bd9ke4CrgksO4d8o43CmAzz33XFavXs3WrVv567/+ayUTETlIJuvEIxgYEiqryutMYKO7b3L3IeA2YFX+Be7e6e6PA6li751KDncK4GOPPZannnqKp59+mne9611liFxEKl066yQjGLoeKqiEArQDW/L2O4CzJuHeipObAjjnmmuuAUamAO7u7mbjxo28973v5ac//SmXX3758LWJRILPfOYzkx6ziEwN6Ux076FUUkIp9A291Pea2RXAFQCLFi2a4OMrw6GmAO7u7ubaa6/l0ksvpa2trVwhikiFS2c9sl5elZRQOoCj8/YXAttKfa+73wjcCLBixYqJJqyKlD8F8KxZs/j3f//3MkYjIlNBlL28KqkN5XFgmZktMbMq4H3AvZNwr4jIjJGaCVVe7p42s48BPyfo+rva3V8ws4+G528ws3nAGqAJyJrZ3wMnuntPoXvL8kVERCpYJuskI+rlVTEJBcDd7wfuH3PshrztHQTVWRO6V0RERoty6JVKqvISEZGIpTNZEhF1G1ZCERGZQTIqoYiISCkMZTyS+eRBCUVEZEZJZ7JUJVRCERGRI9Rnm/hd+iM8uvXRkj9bCaVCaYItEYlCOpslSwqf8EAkE6eEUqHOO+88HnvsMZ555hkuuOAC7rjjjnKHJCLTQDqbAcCs9NVeFfUeSqX58h+/zLo960r6zOPnHM+nz/z061538803c/vttzM4OMiOHTv44he/yD333MN9991HZ2cnV1555UHzoYiIvJ5MNguAFRwC8cgooVSgw51g69VXX+U73/kOmzdvxsz4wQ9+UKZvICKVKpdQYlb6CiollEOYSEkiCoc7wdYxxxzD5Zdfzmc/+1muv/76yQxZRKaIVDZLkmhKKGpDqUCHO8HWK6+8wuc+9zmuv/566uvryxC5iFS6dK7KK4I2FCWUCpSbYOuRRx7hmmuuYf369cDIBFt33nknN9xwAzt37uSmm24avu/iiy9mzpw5fOlLX2LPnj3lCl9EKlQm67gHvbvUhjLDHWqCLYC1a9eWIywRmSJSmSy5uQfVhiKj5E+wJSLyetJZBwtLKKryEhGRw5VKZ4e31SgvIiKHLb/KSyUUERE5bKmsM9yGEsHPvxKKiMgMkUpn1YYy2XLd6qab6fq9RGRiZlSVl5mtNLOXzGyjmV1d4LyZ2TfC88+a2Rl5514xs+fM7GkzW3O4MdTU1LB79+5p9+Pr7uzevZuamppyhyIiZZLKOMYMeA/FzOLAt4ALgA7gcTO7193zX664CFgWLmcB14frnPPcfdeRxLFw4UI6Ojqm5XDxNTU1LFy4sNxhiEiZpDIjVV7T/T2UM4GN7r4JwMxuA1YB+QllFXCrB8WHx8ys2czmu/v2UgWRTCZHvSwoIjJd5IZdgenfbbgd2JK33xEem+g1DvzCzJ4wsyvG+xAzu8LM1pjZmulYChERGc9QeqSX13RvQyn07cY2ZBzqmre6+xkE1WJXmtm5hT7E3W909xXuvqK1tfXwoxURmWIG0xlmSrfhDuDovP2FwLaJXuPuuXUncDdBFZqIiIQGUjOnl9fjwDIzW2JmVcD7gHvHXHMv8MGwt9fZQLe7bzezejNrBDCzeuBC4PnJDF5EpNINpjPD9TzTegpgd0+b2ceAnwNxYLW7v2BmHw3P3wDcD1wMbAT6gMvC248C7g7/gBLAj9z9gUn+CiIiFW0gNVLlNa27DQO4+/0ESSP/2A152w4cNFVh2DPs1MgDFBGZwvKrvKLoNlxJVV4iIhKhqEsoSigiIjPEQEpjeYmISAkMpDMkwl99lVBEROSwDaQyJGNBIilrCcXM/jKva+4/mtmP8wdnFBGRyjaQypIIu2LFLV7y5xdTQvmf7t5rZucA7wJuIRicUUREpoD+oTTJRFhCKXOVVyZc/ylwvbv/BKgqeUQiIhKJ3oE0NckgkcRj5S2hbDWz/wf4K+B+M6su8n4RESmjnoEUNWGrfLlLKH9F8Bb7SnffB8wBPlXyiEREJBK9A2mqwxJKuV9s/Ky7/9jdNwCEc5CcX/KIREQkEpWUUC4ocOyiUgUiIiLRcXf29g1Rkwx+9ssyY6OZ/Tfg74ClZvZs3qlG4NGSRyQiMgOkMlle3N7D/c/t4A+bd5POOJmsk3XHHTI+st1Uk+Ck9lk0Vid4dXcf27v7ScZjwZKIURW3kf14jKrE2P0Yb39DK31DGRprYtBXvimAfwT8DPgScHXe8V5331PyiEREpoGhdJbt3f1s3dtPx95+Ovb2hetge0fPAFmHRMw4Y9FsZtcliBnEzIIlNrLd1TvIfc9upz+Vob25loWza8m6k0o7/f0pUplsuDhD6ezo/UyWoXSWbzy0AYDGmuBnvywJxd27gW4zuwz4c2Bx7j4zw93/qeRRiYhMAV29g7y0o5eOvX1s3Tc6cezoGcDz5pyNGcyfVUt7cy1nL53Lwtm1LGmt57zlbTTXRfsGxvNbu3n3N38LwLxZVbCtfCWUnHuAbuAJYLDkkYiIVIhUJsuu/YN09gzS2TtIZ+8AO3sG6eodOOhYTjxmzGuqYeHsWt5ybAsLZwclifbZtRw9u455s2pIxsvzpsUb22ex+kMr6B1IszO2FYhmCuBiEspCd19Z8ghERMpgX98QGzr3s2HnfjZ09rKp6wA7ewbo6h1kT9/QqNJFztz6Klobq2lrqmH5vEaWttZz+tGzOXpOLfOaakiUKWFMxDuOPwqAG57JAuUvoTxqZie7+3Mlj0JEJCK79w8GiaNzPxt39rKhcz/rd+5n1/6R0kVtMs7S1noWzq7l9EWzaWuspq2pmrbGmuHtlobqspUwSsk9ugm2ikko5wCXmdkmgiovC2LzU0oelYhIEdydrt4wcYRJY0PnfjZ27mfPgaHh6xqqExzX1sB5y1tZdlQDy9oaOa6tgfbmWmKx0r85XokyHoyiVe455SN/58TMVgL/RjCn/Hfd/box5y08fzHBnPIfcvcnJ3KviEwvG3b28tNnt/PbjbvYsLOXnoH08LmmmgTLjmrkXScdxXFtjSxra2DZUQ3Ma6qJ5Id0Ksl6NpKRhqG4hPIa8AFgqbv/k5ktAuYBr5YiEDOLA98ieIGyA3jczO5197V5l10ELAuXswhGOz5rgveWzq4N0LMNLBYssThYHGLhvsVHjg2fjxU4lrtn7LFwf4b/xZeZrX8oE/acCnpQ5brf5o7t7BnEDE47upk/O3VBmDSC5NHaWD3jE8d4HI/sz6aYhPJtIAu8A/gnoBe4C/iTEsVyJrDR3TcBmNltwCogPymsAm71oBLwMTNrNrP5BF2ZX+/ekvHHrsfW3BTFo8ewAkkmLwnF4hCvhkTVmHW4xKvCdaFraqC2GermQn3ryFI7O3i+SMR6B1IFE8XWcHvX/qFR1ydixvzmGtqba3nbslZOWTiLlSfNo62ppkzfYGrKZDMVUUI5y93PMLOnANx9r5mVsvN0O7Alb7+DoBTyete0T/BeAMzsCuAKgEWLFh1WoP964EJ+P7iAhGVJxqE6DtUxpyoGVbntuFMVc5IxqArPJWNOMuZUWbBOWHBNItxOmpOwLAkjXAf7cQvOx8kStyxxnLg5Mc9AZhDSg5AeCrZTg9jAfiyzGzJDxDJDWHYIywwRzwwSyw4Ryw5hFOjCAkHSqmuB+twSJpq6FqhugGQdVNWPLMl6qKobvZ2oUelKGExn2L5vgG37ggSxdV8/2/b1s37nfjbvOkB3f2rU9VWJGAubg262Jy5oCl/gq6N9dvDuxlFNNcRnSDtHlFLZFMlYMpJnF5NQUmHVkgOYWStBiaVUCv1NGfurN941E7k3OOh+I3AjwIoVK8b5VT200049A+YsJZN1Upksg+GbqUPpLH2ZLPvSwfbQmHNDQ9nht1Zz53PrQl0Uo+MkyTCLA8yxHuZaD22xHhZX97Gw+gDz4/tpGeph9mA3DZ2vUp3aQ2Kod+KPt1jhRFNVDzXNMKsdZh0dLguDpXa2ktAU4u5094+UMLYNJ4yB4eTR1Tv6dTUzaGusZvHcet5z6gLac+9phImjpaFK1VSToFISyjeAu4E2M/tn4C+AfyxhLB3A0Xn7C4FtE7ymagL3lsx5x7dx3vFtJXueu5PO+nCiySWpscknlc4ymHcsnc0Sj8VIxox4zEjEjUQsRmJ4f2Q7GY8Fx8Lr4jFjX1+K7d3Bj8CO7gG2dPfzePcA27sH2L63n4HUyL8XkqSpY4A6BqmzcJ233RAbpMEGqbdgu25okLrUIPV9wfla+qljL7N4kdbsLqoY/a/TQathX9VR9FTNY3/NPA7UzKO/bj6Dde0MNSwg0zCfZFUN1YkY1ckYNYk41ckY1Yk41YkYDTUJjmqsmTE9daKWzmTZ0TPAtjEljFzy2LavnwNDmVH3VCdiw6WJdyxvo312LQuag/325lrmzaqhKqHq1HJLZ9MkYsX89E/chJ/q7j80sycIhqw34BJ3f7GEsTwOLDOzJcBW4H3Afx1zzb3Ax8I2krOAbnffbmZdE7i3YpkZyXBwt/rqyfvc+bNqOWF+U8Fz7h4mnAF29PSzs2cwTGJONhskwEw2SyYLmWxwPONOX8bp9WCQu9HX+kiJLpWhemgPTYM7mDW0k9npncxNd9Ey1EnbYBcLe9bRat2j4sm60cUstnkLW30uG72Frd7Cdp9DL3UMehLiSZobGxlsaCdW00h9VYK66vjodVWc+upwnXe8vjpOXVVi+Nh0eN8AIJt1Utks6YwPj+2UzmaDMaBSmbySRd46HDYkO6bUPLe+igXNtSxtredty1pZ0By8FZ5LGnPqVcKYCiqlhIK7rwPWRRGIu6fN7GMEk3jFgdXu/oKZfTQ8fwNwP0GX4Y0E3YYvO9S9UcQ5U5gZs+urmF1fxYkLCiedKHmqn/S+raT3vEp6zxa8ewu13R0s69nKCfu3UrX/aWKZAiMA9QdLZ6yVV2KL2MhCXsos5InUfNamF9DHxBpwq+KxUUmorjpBfVWYdMLk01A9ej8/KSXjxr7+FAOpDOncj3j4o57OBD/yqfTI8XRuML9sbnvk2nQ2y1B4TToc7C89NkmE9w4/P5MdTuQTkd/gffaxc4dLFQvCNo0Fs2qprYqmIVcmVyqbKn8JJZzy9/8gb3BIoKSDQ7r7/QRJI//YDXnbDlw50Xtl6rJkLcnW40i2Hlf4Anfo2w3dW2DoQNAxITMEqT7Ys5m2rnW0da7jzF0PgA8Gf2MTkG06mqG5y+lvXsb+xuPY13Asu+sW05uppm8ozYHBTLAeytA3GK7D4wcG0+zt6x91Xd+Yap+iv6dBMhYjMTz8eFhtGTeq4rHhasxkIqjaTMSNxmSCZFidOXxPuA6O546NDGE+9tpEzKitig8PVtjaWK0G7xkinU1XRAnlJ2hwSKkUZiM90Q4lk4a9r0DXOuh6kVjnOmq61lHz2n8yOzM00vDWvAhaT4C240fWLW8IOhIcQjYbVB0dGErTNxiuhzIMpbPMqk1Skwx/0ONBQkjGR5JHrn1L1UQymSqihIIGh5SpKJ6AluOC5YR3jxzPpGHv5iDRdAbJhs51sOnXQUkHAAsSTdsJ0Hr8yLrlDUGvNSAWM+qrE9RXJ4Ip50QqXKW0oWhwSJk+4gloWRYsJ/zZyPFMGvZsCks066DzxWC98SHI5nqmGcw+pnCJJllblq8jMlHpbJpkvPwJ5RzgQ2a2GQ0OKdNVPAGtbwgW3jNyPJMaSTT5JZqNvxyTaBYHJZl5p8DSt8P8U1+32kxkMqUylVHlFfngkCIVK56E1uXBcuKqkeOZFOx++eASzfoH4Dfh+KS59pnW5WG12XJoWR6MPCAyydKepiYWzXA1xbyH8qqZnQq8LTz0iLs/E0lUIlNFPBlUd7UdP/p4/1549VHYuTYozXS9NKZ9Bpi1KKwyWx4mnHBbiUYilMqkSFSXv9vwx4GPAD8OD/3AzG50929GEpnIVFY7G47/02DJGe5x9mJe1dlLsOk3wThsObOOHkkubfmJRq3+cuQqpVH+coIBIg8AmNmXgd8DSigiEzGqx9mYjgD7Xh2pLsstrzwC6YGR65oWhiWa/GU51Ez+i6cydVVKQjEg/y2uDIUHZRSRYsQTMPfYYMnv2pzN5L1DkyvRrINXfjsm0bSPJJi2/EQza9K/ilS+/UP7qU9G01GkmITyPeAPZnY3QSJZBUzGpCAiM1MsPpJo8qvOspmwRLNudIlmzWpI949c17hgpFtzfocAJZoZ7UDqAI1V0VSfFtMo/zUze5ig+zDAZe7+VCRRicj4YnGYszRYjr945Hg2A/teO7hEUyjRtC6HthNhwenQHkzHoOkDpr9UNsVAZqD8JRQzqwH+C0EvrywQN7MX3X3gkDeKyOSIxWHOkmBZntfLP5uF7tcOXaKpaR5JLu1vggVnQNP8snwNic7a3cEktosaD29ywddTTJXXrQTT/n4j3H8/8H3gL0sdlIiUUCwWvHA5ezEszxs9KZMOEsvWJ2Dbk7D1Sfjtv4KHTaWNC4IEs+D0MMmcHkwbLVPWIx2PELc4b21/ayTPLyahLHf3U/P2f21meg9FZKqKJ2DeG4PlTZcGx1L9sOO5IMlsfTJINOt+OnLPnGOD5NJ+RlCKmX+KhpuZQtbvXc/ipsXMqo6mHa2YhPKUmZ3t7o8BmNlZwO8iiUpEyiNZC0efGSw5/Xth29NhSeapoJfZc3cE52KJoLE/V03WfkbQCSAezYtzcmQ2d2/muOZxpoQogWL+q58FfNDMXgv3FwEvmtlzaEwvkemrdjYce16w5PRsH6km2/oEvHA3PHFzcC5RG4xhNtwec7oa/StAKpuio7eDdx7zzsg+o5iEoqHrRSTQNB+a8kYCcA8Gz8xVk219Imj0f+zbwfma5pFqsvYzgh5mzYuCjgQyKTp6O0h7miWzlkT2GUWN5RVZFCIytZmNvDNzSthPJ5MK3v4fLsk8Cb/9+kijf7w6KLnkphGYm1sfp8b/CGzq3gTA4qbFkX1GMd2GVwDXAseE95Vs+HozmwPcTjC98CvAX7n73gLXrQT+jWDe+O+6+3Xh8c8RjDPWFV76f4dTAotIucSTQaP9/FPgTR8Kjg31BY3+u16CXRtg90boXAvr7htJNAD1bYUTTfMxap85TJu7NwOwdNbSyD6jmP8yPwQ+BTxH8B5KKV0NPOTu15nZ1eH+p/MvMLM48C3gAqADeNzM7nX3teElX3f3r5Q4LhEppao6WHRWsOTLpIJhZnatDxPNBti1EdbeC/17Rq6LVwWlmrnHjSSbWe1BF+em+RpA8xBe3vcy8+rn0VAV3WjWxSSULne/N6I4VhG8NAlwC/AwYxIKcCaw0d03AZjZbeF9axGRqS2eHCmRjNW3Jy/JrA8Sza71wZwz2fToa6uboGnByNKY224fOVY7e0Z2EFi7ey3Lmgv8+ZZQMQnls2b2XeAhghkbAXD3H49/y4Qd5e7bw+dtN7O2Ate0A1vy9jsIep7lfMzMPgisAT5RqMoMwMyuAK4AWLQomrdFRaSE6uYcolTzKuzfEfQ669kKPdugd1uw7nwRencAPvq+RG3YqSAvyTS1Q+O8ICFVNwTrqoZgO1kfvBw6hW3p3cKm7k1cctwlkX5OMQnlMuB4IMlIlZczMj/KIZnZL4F5BU5dO8HPL/RPitzflOuB/xXu/y/gq8CHCz3E3W8EbgRYsWKFF7pGRKaAeHJkOoDxZNKwf2eQYHIJp2cr9G4Ptl/7fZCMhqdxLsRGkktVQ1Ctlks6rccH1W/JmiBRJfOWRG1wPFkHiZpgmeTE1D3YzYOvPsiP1v2IRCzBRUuinXi3mIRyqruffLgf5O7jdn42s51mNj8sncwHOgtc1gEcnbe/ENgWPntn3rO+A/wUEZF4ImhjmdUO/Enha7JZ6NsdJJnB3mAZ2g+DPTC4P9zuHXNuP/RugJd+NrozwevJJZZkLcxaCEefBfNODqrhMqlgorV0uGSGRrYTVXD8u4NedK+jL9XHw1se5mebf8Zvt/2WdDbN4qbFXPe265hXX+jf9KVTTEJ5zMxOzGsEL6V7gUuB68L1Twpc8ziwzMyWAFuB9wH/FSCXjMLr3gs8H0GMIjIdxWLQ0BosxRrqC6rcUgPBQJup/jHb/cHcNam+Mcf7YPfL8MfvjJ6t81Ae/AwsejOc9gE46ZJRHRBSmRS/2/Y77t98Pw9veZj+dD9tdW184PgPcPHSizlhzgnYJLQbmfvEan3M7EXgWGAzQRtKKbsNzwXuIHj7/jXgL919j5ktIOgefHF43cXAvxJ0G17t7v8cHv8+cBpBldcrwP+Vl2DGtWLFCl+zZs2Rhi8icnjSQ8HcNgPdQRVevBoS4RKvDkomiRo4sAuevQ2e+mHQOSFZT+akS3hyyVncd2AzD776ID1DPcyqnsWFx1zIRUsu4k1HvYmYRVPFZmZPuPuKg44XkVCOKXR8Kr/wqIQiIlOJZ7OsXXsH9z9/Mw/0baEzHqPW4R2Nx3LxKZfz5mNXRja9b77xEkpRb8qb2akE86EAPOLuGm1YROQwuDvbDmzjuV3PsX3/dqriVdQmaqlN1FITr6E2Ga4TtewZ2MPjOx7ngVceYEvvFpKxJOcsejufijXz9k1/pPa5X8MLj8AbVsLpfwvHvbMsL4AW86b8xwneRs/16vqBmd3o7t+MJDIRkQo3kB6gq7+LvlQfB1IHOJA6QF+6b3i/Lz1yvD/dP3JNqo+O/R3sGdjz+h8Silucs+afxUdO/gjnH3M+TVVNwYl3EEye9tT34ZnbgukGGubBae8PkssEGvJLpZgqr2eBN7v7gXC/Hvj9VB5lWFVeIpLOptk7sJc9A3uGf/APpA/Qn+ofNyn0pfrYM7CHl/e9TNrTh3x+MpakLllHfaKeumTd8HZrXSsnt5zMya0ns7hpMalMiv50P/2ZfvrT/QykB4bXdck63tjyxpEkMu6XGYINP4cnvw8bHwTPwjFvDRLLSe8NujGXwBFXeRE0wuf3j8tQ+N0QEZGycnd6hnrY1b9r1LK7fze7B3aPOrZ3YC8+9uXHMZKxJPXJeuoSdcNJobWulXMXnsuipkU0JBuGz9Un64evrU/Wk4xPsE0jCc00H9kXT1TBCX8WLD3b4OkfwVM/gHs+Cs/9v/A3d0U6SkAxCeV7wB/M7O5w/xJgdckjEhEZR1+qb3SCCJPD7v7dBx1Pjx2WBaiKVdFS20JLbQsLGhZwSuspwX5NC3Nq5wwng1xpIpcYJpwUKknTAjj3k3DO/4DffBl+c10wQVr7GZF9ZDGN8l8zs4eBcwhKJpe5+1NRBSYiM8NQZqhgyaFQiaI/3X/Q/TGLMadmDi21LcytnctxzccNJ43csbm1c2mpbaEx2Tgp72NUlFgMzv5okFTWP1AZCcXMbgE+7u5PhvuzzWy1uxcc4kREZq5MNsPewb0FSw5jSxQ9Qz0FnzGrehYtNUFSOLnl5OHkkCtR5Labq5uJa6KuQ6udHUxq1hFtm3ExVV6nuPu+3I677zWz00sfkohUoly7RC4ZjC1R5CeJvYN7yfrBs1zUJeqGSw7HNh/LmfPOHFWayCWNOTVzqIpXleFbTmMLToOND0X6EcUklJiZzc6N4htOiqWZbkSmuL5UH7sHdh9UmsgljfzjqQKDKCZiieFSw/z6+byx5Y2jkkN+iaIuWVeGbygAzFkSDBMz1BfMSxOBYhLCV4FHzexOgiFO/gr450iiEpFDynp2VBfWA+kDI+8+jNnP7/qa/35E71Avu/t305fuO+j5hg23S7TUtrBk1pIgOdQcXJpoqmqaee0SU9HscC75va/AUSdG8hHFNMrfamZrCF6jMeDPIxooUmTayU8AY3/Y848dcj8vKRRKAuOpTdQOd2GtT9ZTm6hlbs1cjmk8ZqQEMaZE0VzdTCKmCohppTmc/6m7o/wJBSBMIDM+iXQPdpPOpplbO7fcocgEuTvpbJqh7BCDmUGGMkMMZcLt7NDo/bHb4T2pTIrBzGCwnU2N+5z847kkUKh30nhqE7Wj3mPIvfOQv5/ftTX3olzuXP5+baJWDdYSqA9HUz5QaHaQ0tA/QQ7D15/4OndtuIs5NXNIWIJYLEbc4iRiCeIWJx6LB+twO2EJYhYb3h57PndvzIp7Tv5IooVezMofBSF3Pv+6QucPdb+7j1pnPVvweP46S7bgfbnn584Dw88rdH6854+XJHI/+PmJ4PVeXpuIRCxBdbyaqlgVVfGqYDsebFfFgv2GZMPwsVEvuxVIAMNJIEwWSgASmeGE0hXZRyihHIYPnvhB6pJ1DKQHyHqWdDZNxjPBkh2zDrfTnmYwPUi/95P29EHnM54Zfs6oZ+ZdV6jXTLkZRsxiGAYGMWKYGYYNr0edD7fzz5sZMWKHPp97BnnXmJGIJYZ/zJuqmoZ/yKvj1SRjSarj1cF2PDluIsg/Xuhc7jlV8arIhgMXiVx1QzB75H4llIqytHkp//An/zDpn5v17HCSyXqWtKeHf2SB4e38BtL888PHxjlf6Pio+wskDBGZQupaoG9XZI8vOqGY2QUEPby+5e5Pm9kV4TztErGYxYhZbFLmOxCRaai6MZi+OCKHU0L5O+Ay4B/Dd1FOK2lEIiISjeoGGOqN7PGHUyHc5e773P2TwIXAn5Q4JhERiUJVQ6QllMNJKPflNtz9auDWIw3CzOaY2YNmtiFczx7nutVm1mlmzx/O/SIiM1p1IwxVQEIxs7lm9t+AOWZ2ppnVApRoxsargYfcfRnwULhfyM3AyiO4X0Rk5qpugMHKqPK6G2gFvgj8C9BtZutKFMcq4JZw+xaCuVYO4u7/CRSaM3NC94uIzGhV0TbKF5NQGt39n4Cd7v524P0Ek26VwlHuvh0gXLdN8v0iItNfdUNQ5TXBqd+LVUwvr4FwPWhmte5+l5n9BvjyRG42s18C8wqcuraIGI6YmV0BXAGwaNGiyfxoEZHyilcDDtk0RDALZTEJ5SthN+HbgdVm9ijQPtGb3f2d450zs51mNt/dt5vZfKDYwWYmfH/4zsyNACtWrIgmTYuIVKJEdbBOD0SSUCZc5eXud7n7Hnf/GnA/cDRB20Up3AtcGm5fCvxkku8XEZn+hhPKUCSPP6yBidz9++7+D+7+QoniuA64wMw2ABeE+5jZAjO7P3eRmf0H8HtguZl1mNnlh7pfRETy5JdQonh8JE8tkrvvBs4vcHwbcHHe/vuLuV9ERPLEw4SSGYzk8Ro6VURkphguoSihiIjIkVBCERGRklBCERGRklAbioiIlESiJlirhCIiIkck9zJjJhXJ45VQRERmilj4pkhWCUVERI5EroSSTUfyeCUUEZGZIldCySihiIjIkVCVl4iIlMRwQlEJRUREjoR6eYmISEkMl1Ay0Tw+kqeKiEjlURuKiIiUhLoNi4hISajbsIiIlIR6eYmISEmYgcXVhiIiIiUQT07vbsNmNsfMHjSzDeF69jjXrTazTjN7fszxz5nZVjN7OlwuLnS/iMiMF0tM+27DVwMPufsy4KFwv5CbgZXjnPu6u58WLvdHEKOIyNQXS0z7Kq9VwC3h9i3AJYUucvf/BPZMUkwiItNPPDntG+WPcvftAOG67TCe8TEzezasFitYZQZgZleY2RozW9PV1XW48YqITE2xxNRvQzGzX5rZ8wWWVSV4/PXAscBpwHbgq+Nd6O43uvsKd1/R2tpago8WEZlCYsnI2lASkTy1AHd/53jnzGynmc139+1mNh/oLPLZO/Oe9R3gp4cfqYjINBab/t2G7wUuDbcvBX5SzM1hEsp5L/D8eNeKiMxoM6AN5TrgAjPbAFwQ7mNmC8xsuMeWmf0H8HtguZl1mNnl4an/bWbPmdmzwHnAf5/c8EVEpogI21AmrcrrUNx9N3B+gePbgIvz9t8/zv1/G110IiLTSCwx7UsoIiIyGZRQRESkJGbAm/IiIjIZVEIREZGSiMVVQhERkRKIxVVCERGREoglwFVCERGRI6U2FBERKQlTG4qIiJSCGuVFRKQkVOUlIiIloV5eIiJSEurlJSIiJaE2FBERKQlTlZeIiJSCBocUEZGSUC8vEREpCbWhiIhISUz3bsNmNsfMHjSzDeF6doFrjjazX5vZi2b2gpl9vJj7RUSEGdFt+GrgIXdfBjwU7o+VBj7h7icAZwNXmtmJRdwvIiIzoA1lFXBLuH0LcMnYC9x9u7s/GW73Ai8C7RO9X0RECLoNexbcS/7oSkkoR7n7dggSB9B2qIvNbDFwOvCHYu83syvMbI2Zrenq6ipF7CIiU0csEawjaJhPlPyJ4zCzXwLzCpy6tsjnNAB3AX/v7j3FxuHuNwI3AqxYsaL0KVpEpJINdgfrnq0w+5iSPnrSEoq7v3O8c2a208zmu/t2M5sPdI5zXZIgmfzQ3X+cd2pC94uIzHh7XwnW258ueUKplCqve4FLw+1LgZ+MvcDMDLgJeNHdv1bs/SIiEq1KSSjXAReY2QbggnAfM1tgZveH17wV+FvgHWb2dLhcfKj7RURkHBE0yk9aldehuPtu4PwCx7cBF4fbvwWsmPtFRGTyVEoJRUREJkM2G6yt9D//SigiIjOJhRU9NU0lf3RFVHmJiMgk+dOvQuvxsPhtJX+0EoqIyEzSOA/O/5+RPFpVXiIiUhJKKCIiUhJKKCIiUhJKKCIiUhJKKCIiUhJKKCIiUhJKKCIiUhJKKCIiUhLmEYw4OVWYWRfw6mHe3gLsKmE4UVKs0ZlK8SrWaMzEWI9x99axB2d0QjkSZrbG3VeUO46JUKzRmUrxKtZoKNYRqvISEZGSUEIREZGSUEI5fDeWO4AiKNboTKV4FWs0FGtIbSgiIlISKqGIiEhJKKGIiEhJKKEcBjNbaWYvmdlGM7u6AuI52sx+bWYvmtkLZvbx8PgcM3vQzDaE69l591wTxv+Smb1rkuONm9lTZvbTSo4z/PxmM7vTzNaFf75vrtR4zey/h//9nzez/zCzmkqJ1cxWm1mnmT2fd6zo2MzsTWb2XHjuG2a5+Wwjj/Vfwr8Dz5rZ3WbWXAmxjhdv3rlPmpmbWcukxOvuWopYgDjwMrAUqAKeAU4sc0zzgTPC7UZgPXAi8L+Bq8PjVwNfDrdPDOOuBpaE3yc+ifH+D+BHwE/D/YqMM4zhFuD/DLergOZKjBdoBzYDteH+HcCHKiVW4FzgDOD5vGNFxwb8EXgzYMDPgIsmKdYLgUS4/eVKiXW8eMPjRwM/J3h5u2Uy4lUJpXhnAhvdfZO7DwG3AavKGZC7b3f3J8PtXuBFgh+YVQQ/iITrS8LtVcBt7j7o7puBjQTfK3JmthD4U+C7eYcrLk4AM2si+J/1JgB3H3L3fZUaL8GU3rVmlgDqgG2VEqu7/yewZ8zhomIzs/lAk7v/3oNfwFvz7ok0Vnf/hbunw93HgIWVEOt48Ya+DvwDkN/zKtJ4lVCK1w5sydvvCI9VBDNbDJwO/AE4yt23Q5B0gLbwsnJ+h38l+EuezTtWiXFCUArtAr4XVtF918zqKzFed98KfAV4DdgOdLv7Lyox1jzFxtYebo89Ptk+TPAveKjQWM3sPcBWd39mzKlI41VCKV6hesWK6HttZg3AXcDfu3vPoS4tcCzy72Bm7wY63f2Jid5S4Nhk/lknCKoSrnf304EDBFUz4ylbvGH7wyqCaowFQL2Z/c2hbilwrCL+HjN+bGWP2cyuBdLAD3OHClxW1ljNrA64FvhModMFjpUsXiWU4nUQ1E3mLCSoWigrM0sSJJMfuvuPw8M7w6Is4bozPF6u7/BW4D1m9gpBVeE7zOwHFRhnTgfQ4e5/CPfvJEgwlRjvO4HN7t7l7ingx8BbKjTWnGJj62Ckqin/+KQws0uBdwMfCKuFoDJjPZbgHxbPhP+vLQSeNLN5RByvEkrxHgeWmdkSM6sC3gfcW86Awt4YNwEvuvvX8k7dC1wabl8K/CTv+PvMrNrMlgDLCBrkIuXu17j7QndfTPDn9it3/5tKizMv3h3AFjNbHh46H1hbofG+BpxtZnXh34fzCdrSKjHWnKJiC6vFes3s7PA7fjDvnkiZ2Urg08B73L1vzHeoqFjd/Tl3b3P3xeH/ax0EnXZ2RB5vFL0OpvsCXEzQk+pl4NoKiOccguLps8DT4XIxMBd4CNgQrufk3XNtGP9LRNT75HVi/i+M9PKq5DhPA9aEf7b3ALMrNV7g88A64Hng+wQ9eSoiVuA/CNp2UgQ/cJcfTmzAivD7vQz8O+FoH5MQ60aCtofc/183VEKs48U75vwrhL28oo5XQ6+IiEhJqMpLRERKQglFRERKQglFRERKQglFRERKQglFRERKQglFpEwsGMn478Yce8DM2s3sYTN7LX/EVzO7x8z2T36kIhOjhCJSPs3AcEIxs1qCdzG2hof2EYwuQDhc+vzJDU+kOEooIuVzHXCsmT1tZv9C8LLnw3nnbyMYUQDgzwmGUxGpWEooIuVzNfCyu5/m7p8CLgIeyDv/EHCumcUJEsvtZYhRZMKUUEQqx1uB3+btZ8L9vyaYOOuVcgQlMlFKKCIVwMyWAls8mLQt323ANwlmYBSpaEooIuXTSzBlMxxc3ZXzCPAlggEARSpaotwBiMxU7r7bzH5nZs8TjBb7kQLXOMFMjCIVT6MNi5SZmVUDv3P3FeWOReRIKKGIiEhJqA1FRERKQglFRERKQglFRERKQglFRERKQglFRERKQglFRERK4v8H7CfQ9zlWIAEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(time, secondary_spin[:, 0], label=r'$a_{2, x}$')\n", "plt.plot(time, secondary_spin[:, 1], label=r'$a_{2, y}$')\n", "plt.plot(time, secondary_spin[:, 2], label=r'$a_{2, z}$')\n", "plt.legend()\n", "plt.xlabel('t/M')\n", "plt.ylabel(f'$a_{2}$ components')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## You can access any data that is available for the given compact object" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ]\n" ] } ], "source": [ "available_columns = primary_compact_object.available_data_columns\n", "print(available_columns)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEGCAYAAABy53LJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAAsTAAALEwEAmpwYAAAuQElEQVR4nO3dd3xc9Znv8c93Rr1Ykm3JknvBYIxjAxZgIJheA3HCQhYSsskG4pDABrh7E9i79ya72bs3BLKkAAkhgYUklBAIJQk1pndEccM27t1Y7kWy1Z77xzmSx0KyRuVoRvLzfr3mNTO/c86cR0Keh1+XmeGcc84lK5bqAJxzzvUtnjicc851iicO55xzneKJwznnXKd44nDOOdcpGakOoDcMHjzYRo8eneownHOuT3n33Xc3mVlp6/KDInGMHj2aqqqqVIfhnHN9iqSVbZV7U5VzzrlO8cThnHOuUyJLHJLulrRR0ryEsoGSnpO0OHwuaeO6EZJekLRA0nxJ1yQc+zdJayV9ED7Oiyp+55xzbYuyxnEPcE6rshuAWWY2HpgVvm+tAfhnMzscmAZcJWliwvGfmNmR4ePJCOJ2zjl3AJElDjN7GdjSqngGcG/4+l7gc21ct97M3gtf7wQWAMOiitM551zn9HYfxxAzWw9BggDKDnSypNHAUcBbCcVXS5oTNoV9oqkr4dqZkqokVVVXV/dA6M455yCNO8clFQCPANea2Y6w+JfAOOBIYD3wX+1db2Z3mlmlmVWWln5iGLJzzrku6u3E8bGkCoDweWNbJ0nKJEga95nZn5rLzexjM2s0sybg18CxkQe8Yw9Pz1sf9W2cc67P6O3E8QTwlfD1V4DHW58gScBdwAIzu6XVsYqEt58H5hGxi+94gyt//x57GxqjvpVzzvUJUQ7HfQB4AzhM0hpJlwM3AmdKWgycGb5H0lBJzSOkTgS+DJzWxrDbmyTNlTQHOBW4Lqr4m63aUgPAtpr6qG/lnHN9QmRLjpjZpe0cOr2Nc9cB54WvXwXUzmd+uccCTFJmXNQ3Grv3NvT2rZ1zLi2lbed4usjOiANQU+dNVc45B544OpSdEfyKPHE451zAE0cH9iUOb6pyzjnwxNGhnExvqnLOuUSeODqQ5U1Vzjm3H08cHchuqXF4U5VzzoEnjg5557hzzu3PE0cHPHE459z+PHF0oMkMgJq9DZgZFr53zrmDlSeODjQ0homjvpGf/G0xp9/yEvWNTSmOyjnnUscTRwcSaxw/n7WYZdW7eXfl1hRH5ZxzqRPZWlX9RWNTkDh2J/RxvLF0M6s217Bp917+8YQxLNu0i/Xb9lCYk8GuvQ3U1DUSj4mde+oRIh4TGXGRnREjJzNObmacnJZHjOyM4DknM05MQoKYRJMZjU1GQ5PR2Gg0NDXR0Op9y/Emo76xiayMGEW5mRTlZlKQnUGw2LBzzvUcTxwdCFuqWB2ukgvwuzdXsmV3HQC3P79kv6SSTuIxUZybycD8LErysijOy2x5zs/OIDsjRlZGDDNaklSjGU1NRmMTLa8bmozsjBiFORnhI/jMssJsSguzPUE5d5DxxNGBxqagP2Phhp0ADCvOZe22WrIyYtx26VE8PnsdE4YUcsIhg6ita6IgJ4P8rDiNZhTmZALQ0NhEfaNR19BEbX0je1oeTcFzQyN765vY09AYfIk3GU0G8RjEYzEyYvtqLRmt3sdjCt8H5XsbGtleW9/y2FpTz5ZddWytqWPl5ho+WL2NbTX11CXZTxOPibh0wPOzMmIMzMuiJD+LgflBchqYv+/R/L4kL4uivEwyYvtqVTEFP0NBdgbxmCcf5/oCTxwdaP19+cXjRnLzM4v47tmHcdYR5Zx1RHlqAuumhsYm9jY0UdfQFHyJhwkiHtv3ZR4TLTWJxiZjd10DO/c0sKO2ni2769i4cw/VO/eyZXc9W3fXsXl3kKA+XLeDzbvr2F6b/B4mEhRmZ1Ac1oiam9sG5mcxvCSX4SV5jCjJY3hJLsV5mV7DcS6FPHF0oKlp/+G3F00dzmXHjaIoLzNFEfWMjHiMjHiM/Ozkzo/HxICcTAbkZDKsODepaxoam9hWm5BUwmTSFDaNNSU0he3Y08D2muD4ttp6ttXUs2ZrLZt27mVnq71QCrIzWpJJ8JzLiIHB69LCbLLjcbLCZjivxTjX8zxxdKChaV+VIzMuSguyifmXUVIy4jEGF2QzuCCb8d34nO219azZWsPqLbWs2VrDmq3NzzW8sXTTAfuYMmIiKyNGbmacAbmZlORlUlaYw8hBeYwcmMf4sgImlA/o8/8j4FxviixxSLobOB/YaGaTwrKBwB+A0cAK4Atm9omxrZLOAX4GxIHfmNmNnbm+JyVWOIpyszxppEDQbFXEEUOLPnHMzFpqJ6u31rB5196gCa4xaIZrftTWN7JjTwNbd9expHoXzy/aSF3Dvv8pqCjK4bDyQiaUD2BCeSETKgoZO7igZZFL59w+UdY47gFuA36bUHYDMMvMbpR0Q/j++sSLJMWB2wn2JF8DvCPpCTP7MJnre1pjk5EVj1HX2ERm3JNGupFESX7QMf+p4Z9MLO1pajI27NjDRx/vZOGGnSzasJMF63fw2pJN1IdD6TJiYlxpARMqCjmsvJBDywoZWJBFUW7QZFeQnUFOZsz7W9xBJ8o9x1+WNLpV8QzglPD1vcCLfPKL/1hgiZktA5D0YHjdh0le36Mam4xxZQUsWL+DI0cUR3kr14tiMTG0OJehxbmcclhZS3l9YxPLqnezcMOOloRStWIrj3+wrs3Pyc+KU16UQ0VRLhVFOVQU5VBelEtFcQ6Hlw+gvCint34k53pNb/dxDDGz9QBmtl5SWRvnDANWJ7xfAxzXiesBkDQTmAkwcuTILgfc2GRMGjqAr504mrMm9s0RVC55mfEYh5UHNYwZCeXba+tZWr2L7TX17NgTDHXetbeBTTvrWL+9lvXb9/Dy4mo27txL4nJmw0tyOWb0QCpHl3DM6IEcUlrgzZ2uz0vHzvG2/lV1emVBM7sTuBOgsrKyyysTNpqRERcXV47o6ke4fqAoN5OjR5Z0eF59YxMbd+5l3bZa5qzZTtWKLbyyeBOPvr+25XMqR5VwzJiBHDO6hEnDisjOiEcdvnM9qrcTx8eSKsLaQgWwsY1z1gCJ39LDgeZ2gmSu71FNTUbM27BdkjLjMYYV5zKsOKhpXP7pMZgZKzfX8M6KLVSt2Mo7K7cwa2Hwp5uVEeOIoQMYMyifYSW5lBVmUzYgh7LCbIYMyKG0MJvMuHfQu/TS24njCeArwI3h8+NtnPMOMF7SGGAtcAnwxU5c36MazXwugOsWSYwenM/owfktNddNu/by7sqtVK3Ywty123lz2WY27NhDUxt140H5WZSGieSQsgKOGV3C1FEDKS1MchKOcz0syuG4DxB0ZA+WtAb4PsEX/kOSLgdWAReH5w4lGHZ7npk1SLoaeIZgOO7dZjY//Ng2r49SY6MnDtfzBhdkc/YR5ZydsPJAY5OxefdeNu7Yy8c79rBx577njTv2sGHHHt54czN3vbocgAnlhZw2oYwzJw5hyvBi7ztxvSbKUVWXtnPo9DbOXQecl/D+SeDJNs7b3Nb1UWo0I+5NVa4XxGOirDCHssIcJg1re2jx3oZG5q3dwTsrtvDioo386uVl/OLFpZQPyOH0w8s45bAyThg3iPzsdOy+dP2F/3V1oLHJaxwufWRnxJk6qoSpo0q48uRxbK+pZ9bCj3lm/gYee38t9721iqx4jKNGFnPE0CImVBRyePkAxg8pICfTO+Fdz/DE0YEmM28CcGmrKC+TC48ezoVHD6euoYmqFVt4YdFG3lq+hfvfXsme+mB2fExQWpjNsOJcRg4MllsZOSifkQPzGFuaz+AC7y9xyfPE0YGGJiPDE4frA7IyYpxwyGBOOGQwENSWV27ezcINwez4ddtqWbu1lndWbOWJ2ev264gfV5rPiYcM5syJQzh+7CAyfCSXOwBPHAdgZpjhw3FdnxSPibGlBYwtLeC8T1Xsd6yuoYm122pZtaWGhet38PrSzfyxag2/fWMlgwuyuWBKBRdPHcHEoQNSFL1LZ544DqB521jv43D9TVZGjDGD8xkzOJ+TDy3lGyePY099Iy8s3MjjH6zjvjdX8d+vrWBCeSHnTqrgM5PLOaSsMNVhuzThieMAGs0Thzt45GTGOfdTFZz7qQq21dTx2Ptr+evc9fx01kf85G8fceiQoOZyxuFDmFBe6M1ZBzFPHAfgNQ53sCrOy+KrJ47hqyeO4eMde3hq7nqenLeBn81azE//tpiczBiHlBUwrjThUZbP6EH5PnrrIJB04pCUb2a7owwm3bQkDu/jcAexIQNyWpLIxh17eGPZZuas2c7ijbs+sXKwBBPKB3DWxCGcdcQQJlYM8GXn+6EOE4ekE4DfAAXASElTgG+Y2beiDi7Vmjf/8+G4zgXKBuQw48hhzDhyWEtZbV0jyzbtYmn1bpZs3MWbSzdz6/OL+dmsxQwrzuX8yRVcNm0UIwbmpTBy15OSqXH8BDibYJ0ozGy2pOmRRpUmmvs4fDiuc+3LzYpzxNCEHRrPDNbien7BRp6Zv4HfvLqcX7+yjLMmlvP16WOYOmpgagN23ZZUU5WZrW5V3Wx/k+d+pHm/ca9xONc5gwuy+cIxI/jCMSNYv72W376xkvvfWsXT8zdw1Mhivn7SWM6cOMRX/u2jkkkcq8PmKpOUBXwbWBBtWOmhuanK+zic67qKolyuP2cC/3TaIfyxag13vbqcb933HsV5mZw9sZzPTK5g2thBvr97H5JM4rgS+BnBznxrgGeBq6IMKl3sG46b4kCc6wfysjL4ygmjuWzaKF5ctJG/zFnPX+eu5w9Vq8mMi0OHFDJpaBGThg3giGFFTCgvJC/LB36mow7/q5jZJuBLvRBL2mlqGY7rmcO5nhKPidMPH8Lphw9hT30jryzexLsrtzJ/3Xae/XADf6gKdo6OCaaOKuH0w4dwxuFljCst8BFaaSKZUVWlwNeB0Ynnm9nXogsrPTQ0eY3DuSjlZMY5c+IQzpw4BAiW+Vm/fQ/z1m5n9pptvLCwmhufWsiNTy1kbGk+M08ay99NHe59IymWTD3wceAV4G8cJJ3izZrncfhaVc71DkkMLc5laHEuZx1RznfOnsC6bbXMWriRP1at5oY/zeXW55dw1amHcNHU4d4vkiLJJI48M7s+8kjSUJMvOeJcyg0tzuXL00Zx2XEjeXFRNT+dtZj/9ehcbn9hCd88ZRwXVw4nO8Nnq/emZNL1XySd1/FpyZN0jaR5kuZLuraN49+R9EH4mCepUdLA8NgKSXPDY1U9GVdrzTUOn8fhXOpJ4tQJZTz2rRO45x+PoWxANv/7sXmcevOL3PXqcnbuqU91iAeNZGoc1wD/S9JeoB4QYGbWpfWWJU0i6DM5FqgDnpb0VzNb3HyOmd0M3ByefwFwnZltSfiYU8NO+0h5U5Vz6UcSpxxWxsmHlvLK4k3c+vxi/uMvH/LjZxZxxsQhnB8O7y3KzUx1qP1WMqOqenot5cOBN82sBkDSS8DngZvaOf9S4IEejiEpvsihc+lLEtMPLWX6oaXMXr2Nh6pW8+Tc9fx59jokmFgxgMpRJRw9qoTjxgyivCgn1SH3G0kNkpZUAowHWn7zZvZyF+85D/hPSYOAWuA8oM0mJ0l5wDnA1QnFBjwryYBfmdmdXYyjQ83zOHzmuHPpbcqIYqaMKObfPnsEVSu28tbyzby1bAt/fHcN976xEoAJ5YV87qhhXHjUMMoGeBLpjmSG415B0Fw1HPgAmAa8AZzWlRua2QJJPwKeA3YBs4GGdk6/AHitVTPViWa2TlIZ8JykhW0lMUkzgZkAI0eO7EqoLTWOTJ/H4VyfkBmPcfy4QRw/bhAQ/BtesH4Hry/dxDPzP+bGpxbyX88u4kvHjeLq0w7xvda7KJlvxGuAY4CVZnYqcBRQ3Z2bmtldZna0mU0HtgCL2zn1Elo1U5nZuvB5I/AoQV9JW/e408wqzayytLS0S3E2NDbXOLp0uXMuxeIxMWlYETOnj+ORb57ArH8+mYumDud3b67k5Jte4CfPfcTuve39f6trTzJfiXvMbA+ApGwzWwgc1p2bhrUFJI0ELqSNPgxJRcDJBPNImsvyJRU2vwbOImj6isS+UVWeOZzrD8aVFvDDCyfz7HXTmX5oKT+btZgzbnmJv8xZh4VN065jyfRxrJFUDDxG0DS0FVh3wCs69kjYx1EPXGVmWyVdCWBmd4TnfB54ttXmUUOAR8NlBzKA+83s6W7G0i7fOta5/mlcaQG/vGwq767cwv95bD5X3/8+945ewTWnH8qJhwzypU06oM5kWUknA0XA02ZWF1lUPayystKqqjo/5eP5hR/ztXuqePyqE5kyorjnA3POpVxjk/HA26u47fklbNixh6mjSrjq1HFMH1960O+rLuldM6tsXd5ujaN5wl0rc8PnAoK+iX6tuY/DaxzO9V/xmLhs2igurhzOQ1Vr+OULS/jaPVUMys/ivE9V8NkjhzJ1ZImPrkxwoKaqdwmGvrb12zJgbCQRpZGWPo64/8E4199lZ8T58rRR/H3lCJ5fuJE/z17HQ1Wr+d2bK6koyuH8yRV8dsowJg3zfdTbTRxmNqY3A0lHDb7kiHMHnayMGOdMKuecSeXs2tvArAUf88QH67jn9RX8+pXljBmczzemj+WiqcMP2qasZCcAXgh8mqCm8YqZPRZlUOmi0ffjcO6gVpCdwYwjhzHjyGFsq6njmfkbuP/tYJXe37+1khsvnMykYUWpDrPXdfiNKOkXBLsAziUY+nqlpNujDiwdtOzHcZBXS51zUJyXxd8fM5LHvnUCt33xKDZs38uM21/jxqcWsrfhoNpxIqkax8nAJAuHX0m6l32d5P1aY7jpeNz7OJxzIUmcP3koJx1Syn8++SF3vLSUFxdt5McXTzloah/JtMEsAhLX7BgBzIkmnPTSGOQN7+Nwzn1CUV4mN100hbu/Wsnm3XV87vbX+MGfP2Tzrr2pDi1yySSOQcACSS9KehH4ECiV9ISkJyKNLsVaahyeOJxz7ThtwhCeu246F00dzj2vL2d6uJRJf94fJJmmqu9FHkWaqmv0RQ6dcx0rzsvixr+bzBUnjeWW5xbxs1mLufeNFXzpuJF8edrofrekezL7cbzUG4Gko5pw8bPcLN+W0jnXsUPKCvjFl6YyZ802bnt+Cb94cSm/emkZ532qgm+dOo4J5V3a/y7tJDUc92C1rbae7IwYWRle43DOJW/y8GLu/IdKVm2u4d43VvDQO6v569z1fO3E0Vx7xqHkZ/ftr17/RjyAylElXP7pg34epHOui0YOyuP/nD+RV64/lS9UDufXryznzFte4vWlke98HakOFzmUdD7wpJk19U5IPa+rixw651xPenflFr7z8BxWbNrNd86ewJUnj03r5UvaW+QwmRrHJcBiSTdJOrznQ3POuYPD1FEDeeLqT3PupAp+9PRCvv7bKrbV9JmFxlt0mDjM7DKCXf+WAv8t6Q1JM5s3VHLOOZe8guwMbvviUfzbBRN56aNqPvPzV3lv1dZUh9UpSfVxmNkO4BHgQaCCYJOl9yT9U4SxOedcvySJr544hoevPIFYDC6+4w1++OQCauv6xtIlyaxVdYGkR4HngUzgWDM7F5gC/M+I43POuX5ryohi/vJPJ3Hx1OH86uVlnPXTl3h1cfp3nCdT47gY+ImZTTazm81sI4CZ1QBfizQ655zr54pyM7nx7ybzh5nTyIzFuOyut7j+4Tlsr03fmefJJI4bzezlxAJJpwCY2ayu3FTSNZLmSZov6do2jp8iabukD8LH9xKOnSNpkaQlkm7oyv2dcy7dHDd2EE9ecxLfPGUcD7+3hjNveYk5a7alOqw2JZM4HpL0XQVyJd0K/LCrN5Q0Cfg6cCxBc9f5ksa3ceorZnZk+PhBeG0cuB04F5gIXCppYldjcc65dJKTGef6cybw2LdOJCsjxhd//RZvL0+/XbqTSRzHEayO+zrwDrAOOLEb9zwceNPMasysAXiJoLM9GccCS8xsmZnVEXTWz+hGLM45l3Y+NbyIP155PGUDsvmHu9/izWWbUx3SfpJJHPVALZAL5ADLuzkZcB4wXdIgSXnAeQRLtbd2vKTZkp6SdERYNgxYnXDOmrDsE8Ihw1WSqqqrq7sRrnPO9b6Kolwe+sbxDC/J44p7q5i9eluqQ2qRTOJ4hyBxHEOwfeylkh7u6g3NbAHwI+A54GlgNtDQ6rT3gFFmNgW4FXgsLG9rimWbU9/N7E4zqzSzytLS0q6G65xzKTO4IJvfX34cJfmZfOW/32bRhp2pDglILnFcbmbfM7N6M9tgZjOAx7tzUzO7y8yONrPpwBZgcavjO8xsV/j6SSBT0mCCGkZi7WQ4QdOZc871S+VFOdx3+TSy4sGIq7lrtqc6pKRmjn9ikScz+113biqpLHweCVwIPNDqeLnCBVwkHRvGuZmg9jNe0hhJWQTLofTrzaScc27koDzuu+I4YoLP/eI1fvjUAvbUp26yYKpWx31E0ofAn4GrzGyrpCslXRkevwiYJ2k28HPgEgs0AFcDzwALgIfMbH4qfgDnnOtN44cU8ux1JweTBV9axmd+/grvp2ipkg5Xx+0PfHVc51x/8vJH1Vz/yBw27tzLfVccx7SxgyK5T3dWx3XOOZdGph9aytPXTmfkwDz++aHZvd5slcxaVRdKWhzO5N4haaekHb0RnHPOubYV5WbygxlHsHZbLU/M7t0xQsnUOG4CPmtmRWY2wMwKzax/bJzrnHN92KcPGcxhQwp54O1V+5XX1DUwf110o6+SSRwfh3MvnHPOpRFJnD+5gvdXbWPjjj0t5dc/MpfP/PzVyBZKTCZxVEn6g6RLw2arCyVdGEk0zjnnOuXsSeUAPLfg45ayFxZuBGDnntQljgFADXAWcEH4OD+SaJxzznXK+LICKopyeGPpvvWsasPO8qg6zTM6OsHM/jGSOzvnnOs2SRwzeiBvLd+MmSGJxqZgmkVtXXeWFWxfu4lD0nfN7KZwGfVPTPYws29HEpFzzrlOOWbMQJ6YvY5VW2oYNSi/pbw2BTWO5g5xnznnnHNp7NjRAwF4e/mW/RJHQ2Mv1zjM7M/h872R3Nk551yPGF9WQHFeJm8v38JFU4e3lDc0RbMyiM8cd865Pi4WC/o53l6xhZ179+1S0eiJwznnXHuOGzOQlZtr9tvwyWsczjnn2nXsmKCf47H39y0/ElUfRzJrVR0qaZakeeH7yZL+dyTROOec65KJFQMoyM7gkffWtJSlssbxa+BfCPYex8zmEGyg5JxzLk1kxGOcNH7wfmWp7OPIM7O3W5W13iPcOedcin3j5HHkZsa5YMpQILoaR4czx4FNksYRTgKUdBGwPpJonHPOddmRI4r54PtnsmV3HX+evS51fRzAVcCvgAmS1gLXAt/szk0lXSNpnqT5kq5t4/iXJM0JH69LmpJwbIWkuZI+kOSTE51zLkF2Rpx4TEAKaxxmtgw4Q1I+EDOznd25oaRJwNeBY4E64GlJfzWzxQmnLQdODvciPxe4Ezgu4fipZrapO3E451x/lREL6gQp6+OQ9P8kFZvZbjPbKalE0v/txj0PB940sxozawBeAj6feIKZvW5mzbuwvwkMxznnXFKirnEk01R1rplta34TfqGf1417zgOmSxokKS/8rBEHOP9y4KmE9wY8K+ldSTPbu0jSTElVkqqqq6u7Ea5zzvUtGc2Jo7fXqkoQl5RtZnsBJOUC2V29oZktkPQj4DlgFzCbdkZpSTqVIHF8OqH4RDNbJ6kMeE7SQjN7uY373EnQxEVlZWU0adc559JQTEHiiOqLL5kax++BWZIul/Q1gi/8bi18aGZ3mdnRZjYd2AIsbn2OpMnAb4AZZrY54dp14fNG4FGCvhLnnHOhMG/QZKnrHL9J0lzgdEDAf5jZM925qaQyM9soaSRwIXB8q+MjgT8BXzazjxLKWzrow9dnAT/oTizOOdffNNc4mlI4jwMze4r9+xm66xFJgwhmo18Vjp66MrzXHcD3gEHALxT8AhrMrBIYAjwalmUA95vZ0z0Yl3PO9XnNneMR5Y2OE4ekC4EfAWUENQ4BZmYDunpTMzupjbI7El5fAVzRxjnLgCmty51zzu0TS3VTFXATcIGZLejwTOeccyknRVvjSKZz/GNPGs4517fElNo+jipJfwAeA/Y2F5rZnyKJyDnnXLfFY0ppU9UAoIZgBFMzIxj15JxzLg1JSl3nuJn9YzS3ds45F5WYwCKqcfgOgM451w/FFF1Tle8A6Jxz/VBMIqKlqnwHQOec649iim4eRzKJw3cAdM65PiYWU2R9HMmMqrqKYJXZ5h0AlwNfiiQa55xzPSKWqlFVkuLAN82sx3YAdM45F70om6oOmDjMrFHS1PD17kgicM451+MU4aiqZJqq3pf0BPBHoCV5+Mxx55xLX3GJpohGVSWTOAYCm4HTEsp85rhzzqWxlDVVgc8cd865vijKJUd85rhzzvVDsVhq53H4zHHnnOtjUr3kSI/PHJd0jaR5kuZLuraN45L0c0lLJM2RdHTCsXMkLQqP3dCdOJxzrr+Kp7Kpih6eOS5pEvB14FiCbWDPlzS+1WnnAuPDx0zgl+G1ceD28PhE4FJJE7sai3PO9VdK8ZIjVwG/Yt/M8WuBK7txz8OBN82sxswagJeAz7c6ZwbwWwu8CRRLqiBINkvMbJmZ1QEPhuc655xLEFN0S460mzgkXRO+rDCzM4BSYIKZfdrMVnbjnvOA6ZIGScoDzgNGtDpnGLA64f2asKy98rbinympSlJVdXV1N8J1zrm+J1gdt/drHM3DcG+FYOZ4Tyw3Eu5f/iPgOeBpYDaf7DNRW5ceoLyt+9xpZpVmVllaWtqNiJ1zru+JxVKzVtUCSSuAUklzEsoFmJlN7upNzewu4C4ASf+PoOaQaA3710KGA+uArHbKnXPOJYhyB8B2E4eZXSqpHHgG+GxP3lRSmZltlDQSuBA4vtUpTwBXS3oQOA7YbmbrJVUD4yWNAdYSDAv+Yk/G5pxz/UHKVsc1sw0EI5962iOSBhHMDbnKzLZKujK85x3AkwR9H0uAGsJmMzNrkHQ1QTKLA3eb2fwI4nPOuT4tpUuOSDof+A9gVHh+c1PVgK7e1MxOaqPsjoTXRjCaq61rnyRILM4559qhCDvHk1nk8KcEzUlzLaoGM+eccz0qHhNRfWMnM49jNTDPk4ZzzvUdKW2qAr4LPCnpJWBvc6GZ3RJJRM4557ot1Rs5/SewC8ghGA7rnHMuzcVEajdyMrOzorm9c865KMQkGiyazJFMH8ffJHnicM65PiQeS+2y6lcBT0uqlbRD0k5JOyKJxjnnXI+IcgfAZLaOLYzm1s4556KSkiVHmkma3la5mb3c8+E455zrCTGJxhSOqvpOwuscgj0x3gVOiyQi55xz3RYsqx7NZyfTVHVB4ntJI4CbognHOedcT4iyqSqZzvHW1gCTejoQ55xzPSfYATCaz06mj+NW9m2WFAOOJNh8yTnnXJqKxVK75EhVwusG4AEzey2SaJxzzvWIVC858jCwx8waw2DikvLMrCaSiJxzznVblE1VyfRxzAJyE97nAn+LJhznnHM9IcrVcZNJHDlmtqv5Tfg6L5JonHPO9Ygot45NJnHslnR08xtJU4Ha7txU0nWS5kuaJ+kBSTmtjn9H0gfhY56kRkkDw2MrJM0Nj1W1fQfnnDu4KcX7cVwD/FHSuvB9BfD3Xb2hpGHAt4GJZlYr6SHgEuCe5nPM7Gbg5vD8C4DrzGxLwsecamabuhqDc871dykbjispDpwETAAOI9hvfKGZ1ffAfXMl1RM0e607wLmXAg90837OOXdQSVkfRziSaoaZ1ZvZPDOb292kYWZrgR8Dq4D1wHYze7atcyXlAecAjyR+BPCspHclzWzvPpJmSqqSVFVdXd2dkJ1zrs+JRTgcN5k+jtck3SbpJElHNz+6ekNJJcAMYAwwFMiXdFk7p18AvNaqmepEMzsaOBe46gCLMN5pZpVmVllaWtrVcJ1zrk8K+jii+exk+jhOCJ9/kFBmdH2RwzOA5WZWDSDpT+E9ft/GuZfQqpnKzNaFzxslPUqw6KKv1OuccwkkpW5ZdTM7tYfvuQqYFjZD1QKns//sdAAkFQEnA5cllOUDMTPbGb4+i/0TmnPOOZoXOYzms9tNHJIuM7PfS/ofbR03s1u6ckMze0vSw8B7BEuYvA/cKenK8Pgd4amfB541s90Jlw8BHpXUHPv9ZvZ0V+Jwzrn+LMo+jgPVOPLD5x7fAdDMvg98v1XxHa3OuYeEIbph2TJgSk/H45xz/U2UEwDbTRxm9qvw+d+jubVzzrmoRDkBsMNRVZLulVSc8L5E0t2RROOcc65HpHqRw8lmtq35jZltBY6KJhznnHM9IdWLHMbCuRcAhGtGJTOM1znnXIqkqnO82X8Br4cjoQz4AvCfkUTjnHOuRygVnePNzOy34Sq0pxGsVXWhmX0YTTjOOed6QjCPI3U1DsJE4cnCOef6iFTvx+Gcc66PSXXnuHPOuT5G4XDcKJqrPHE451w/FAuWZopkLocnDuec64diQd6IpLnKE4dzzvVDsTBzRNFB7onDOef6IXmNwznnXGd4H4dzzrlO8T4O55xzndJc4/DE4ZxzLilSP+scl3SdpPmS5kl6QFJOq+OnSNou6YPw8b2EY+dIWiRpiaQbej9655xLf81NVf1iAqCkYcC3gUozmwTEgUvaOPUVMzsyfPwgvDYO3A6cC0wELpU0sZdCd865PiPW32ocBIsr5krKAPKAdUledyywxMyWmVkd8CAwI6IYnXOuz+pXneNmthb4MbAKWA9sN7Nn2zj1eEmzJT0l6YiwbBiwOuGcNWHZJ0iaKalKUlV1dXUP/gTOOZf+1J86x8PdBGcAY4ChQL6ky1qd9h4wysymALcCjzVf3sZHtvlbMbM7zazSzCpLS0t7JHbnnOsr+ts8jjOA5WZWbWb1wJ+AExJPMLMdZrYrfP0kkClpMEENY0TCqcNJvpnLOecOGv2qqYqgiWqapDwFdanTgQWJJ0gqD48h6ViCODcD7wDjJY2RlEXQqf5Er0bvnHN9QJSd40ntANiTzOytcP/y94AG4H3gTklXhsfvAC4CvimpAagFLrFgTFmDpKuBZwhGY91tZvN7+2dwzrl017JWVQSZo9cTB4CZfR/4fqviOxKO3wbc1s61TwJPRhedc871ff2tj8M551zEYuG3e3/p43DOORcxX6vKOedcp/S7taqcc85Fq1+tVeWccy56/XGtKueccxHqbxMAnXPORaxfrVXlnHMuej6PwznnXKd4U5VzzrlOWb2lBoD563b0+Gd74nDOuX5o8+46AKp37u3xz/bE4Zxz/VBz53huZrzHPzslixw655yL1jemj2VvfSNfPn5Uj3+2Jw7nnOuH8rMz+JfzDo/ks72pyjnnXKd44nDOOdcpnjicc851SkoSh6TrJM2XNE/SA5JyWh3/kqQ54eN1SVMSjq2QNFfSB5Kqej9655w7uPV64pA0DPg2UGlmkwj2Dr+k1WnLgZPNbDLwH8CdrY6famZHmlll5AE755zbT6pGVWUAuZLqgTxgXeJBM3s94e2bwPBejM0559wB9HqNw8zWAj8GVgHrge1m9uwBLrkceCrxI4BnJb0raWZ0kTrnnGtLKpqqSoAZwBhgKJAv6bJ2zj2VIHFcn1B8opkdDZwLXCVpejvXzpRUJamqurq6R38G55w7mCmKbQUPeEPpYuAcM7s8fP8PwDQz+1ar8yYDjwLnmtlH7XzWvwG7zOzHHdyzGljZxZAHA5u6eG1v81ij05fi9VijcTDGOsrMSlsXpqKPYxUwTVIeUAucDuw3OkrSSOBPwJcTk4akfCBmZjvD12cBP+johm394MmSVNVXOuE91uj0pXg91mh4rPv0euIws7ckPQy8BzQA7wN3SroyPH4H8D1gEPCLcKGuhvCXMAR4NCzLAO43s6d7+2dwzrmDWUpGVZnZ94Hvtyq+I+H4FcAVbVy3DJjSutw551zv8ZnjHWs9hySdeazR6UvxeqzR8FhDvd457pxzrm/zGodzzrlO8cThnHOuUzxxHICkcyQtkrRE0g1pEM8ISS9IWhAuEnlNWD5Q0nOSFofPJQnX/EsY/yJJZ/dyvHFJ70v6SzrHGd6/WNLDkhaGv9/j0zXethYJTZdYJd0taaOkeQllnY5N0tRwMdMlkn6u5n1Qeyfem8O/gzmSHpVUnA7xthVrwrH/KckkDe6VWM3MH208CBZfXAqMBbKA2cDEFMdUARwdvi4EPgImAjcBN4TlNwA/Cl9PDOPOJpipvxSI92K8/wO4H/hL+D4t4wxjuBe4InydBRSnY7zAMIJFQHPD9w8BX02XWIHpwNHAvISyTscGvA0cD4hgyaFzezHes4CM8PWP0iXetmINy0cAzxBMch7cG7F6jaN9xwJLzGyZmdUBDxIslZIyZrbezN4LX+8EFhB8kcwg+OIjfP5c+HoG8KCZ7TWz5cASgp8rcpKGA58BfpNQnHZxAkgaQPCP8i4AM6szs23pGi/7FgnNYN8ioWkRq5m9DGxpVdyp2CRVAAPM7A0Lvul+m3BN5PGa2bNm1hC+TVxkNaXxtvO7BfgJ8F2CdfyaRRqrJ472DQNWJ7xfE5alBUmjgaOAt4AhZrYeguQClIWnpfJn+CnBH3NTQlk6xglBrbIa+O+wae03ClYmSLt4rf1FQtMu1gSdjW1Y+Lp1eSp8jX2LrKZdvJI+C6w1s9mtDkUaqyeO9rXV7pcWY5clFQCPANea2Y4DndpGWeQ/g6TzgY1m9m6yl7RR1pu/6wyCJoBfmtlRwG6CJpX2pCxedWKR0OZL2ihLi79j2o8tLWKW9K8Eq1vc11zUxmkpi1fBsk3/SrDSxicOt1HWY7F64mjfGoK2w2bDabVvSCpIyiRIGveZ2Z/C4o/DKijh88awPFU/w4nAZyWtIGjiO03S79MwzmZrgDVm9lb4/mGCRJKO8Z4BLDezajOrJ1jT7YQ0jbVZZ2Nbw/578PR6zJK+ApwPfCls0oH0i3ccwf9AzA7/rQ0H3pNUHnWsnjja9w4wXtIYSVkEuxQ+kcqAwtEPdwELzOyWhENPAF8JX38FeDyh/BJJ2ZLGAOMJOsYiZWb/YmbDzWw0we/teTO7LN3iTIh3A7Ba0mFh0enAh2kab8sioeHfw+kEfV3pGGuzTsUWNmftlDQt/Bn/IeGayEk6h2Arh8+aWU3CobSK18zmmlmZmY0O/62tIRg8syHyWHu6578/PYDzCEYuLQX+NQ3i+TRBtXIO8EH4OI9gQchZwOLweWDCNf8axr+IiEamdBDzKewbVZXOcR5JsErzHOAxoCRd4wX+HVgIzAN+RzByJi1iBR4g6HupJ/giu7wrsQGV4c+3FLiNcJWLXop3CUH/QPO/sTvSId62Ym11fAXhqKqoY/UlR5xzznWKN1U555zrFE8czjnnOsUTh3POuU7xxOGcc65TPHE455zrFE8czkVMwcq732pV9rSkYZJelLQqcYVSSY9J2tX7kTqXHE8czkWvGGhJHJJyCeYyrA2LthHMtidcwruid8NzrnM8cTgXvRuBcZI+kHQzwaTIFxOOP0gwwx7gQoJlRJxLW544nIveDcBSMzvSzL4DnAs8nXB8FjBdUpwggfwhBTE6lzRPHM71vhOBVxPeN4bv/55gg6YVqQjKuWR54nCuF0kaC6y2YHOwRA8CtxLs6OdcWvPE4Vz0dhJs9QufbKZq9grwQ4KF7JxLaxmpDsC5/s7MNkt6TdI8gtVNv97GOUaws59zac9Xx3Wul0jKBl4zs8pUx+Jcd3jicM451ynex+Gcc65TPHE455zrFE8czjnnOsUTh3POuU7xxOGcc65TPHE455zrlP8PnlJO+UzQuxEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "data = primary_compact_object.get_data_from_columns([CompactObject.Column.TIME, CompactObject.Column.XY_PLANE_CIRCUMFERENCE])\n", "time = data[:, 0]\n", "circumference = data[:, 1]\n", "plt.plot(time, circumference)\n", "plt.xlabel('t/M')\n", "plt.ylabel('circumference in xy plane')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Close the Coalescence object to close the associated h5 file" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "pycharm": { "name": "#%%\n" } }, "outputs": [], "source": [ "coalescence.close()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "mayawaves_dev", "language": "python", "name": "mayawaves_dev" }, "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.8.15" } }, "nbformat": 4, "nbformat_minor": 4 }