diff --git a/20180708Fields/Mercury Interconnect Sept 23 2018.ipynb b/20180708Fields/Mercury Interconnect Sept 23 2018.ipynb index a5d3f0031..7302f63fd 100644 --- a/20180708Fields/Mercury Interconnect Sept 23 2018.ipynb +++ b/20180708Fields/Mercury Interconnect Sept 23 2018.ipynb @@ -49,7 +49,7 @@ }, { "cell_type": "code", - "execution_count": 134, + "execution_count": 235, "metadata": {}, "outputs": [], "source": [ @@ -68,7 +68,7 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": 236, "metadata": {}, "outputs": [], "source": [ @@ -89,7 +89,7 @@ }, { "cell_type": "code", - "execution_count": 54, + "execution_count": 237, "metadata": {}, "outputs": [], "source": [ @@ -224,7 +224,7 @@ }, { "cell_type": "code", - "execution_count": 135, + "execution_count": 238, "metadata": {}, "outputs": [], "source": [ @@ -240,7 +240,7 @@ }, { "cell_type": "code", - "execution_count": 136, + "execution_count": 239, "metadata": {}, "outputs": [], "source": [ @@ -277,12 +277,12 @@ }, { "cell_type": "code", - "execution_count": 124, + "execution_count": 339, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbsAAAGfCAYAAADGVHw+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXusbNld3/n9rb3rcc65j36624+23cYG26AkHS62edlgG2wkhJ0AwUCISZixIg3KH6MZ0SNmgmQUyWgihhmFBEyIMQyWeQyYnhGTlsFJFCEMvh0zhDbY/bBxP3D3fZ9z7jlVtfdav/ljrbX32rv2rlPnnjr3Uff7QdVV+1FVu+qa+pzfY60lqgpCCCFknTE3+gIIIYSQ44ayI4QQsvZQdoQQQtYeyo4QQsjaQ9kRQghZeyg7Qggha89KZCci7xaRz4vIkyLycMfxt4rIfxGRUkS+r3XMisifhdsjq7geQgghJEWOOs5ORDIAXwDwHQCeBfAZAD+oqp9Lznk1gFMA/gcAj6jqbyfHdlX1xJEughBCCFlAvoLXeBOAJ1X1aQAQkY8DeA+ASnaq+qVwzK3g/QghhJBDsQrZvRzAM8n2swDefIjnj0XkLIASwIdU9RNdJ4nIBwB8AAC2tra+/vWvf/01Xi4hhJBbhccee+y8qt571NdZheyOyqtU9TkReQ2AT4nIf1XVp9onqeqHAXwYAM6cOaNnz5693tdJCCHkOiMif72K11lFg8pzAB5Itl8R9i2Fqj4X7p8G8B8BPLSCayKEEEIqViG7zwB4nYg8KCJDAO8DsFRXpYjcKSKj8PgeAN+MpNZHCCGErIIjy05VSwA/DuBRAH8J4DdV9XER+aCIfA8AiMg3iMizAL4fwC+KyOPh6W8AcFZE/j8A/wG+ZkfZEUIIWSlHHnpwI2DNjhBCbg9E5DFVPXPU1+EMKoQQQtYeyo4QQsjaQ9kRQghZeyg7Qgghaw9lRwghZO2h7AghhKw9lB0hhJC1h7IjhBCy9lB2hBBC1h7KjhBCyNpD2RFCCFl7KDtCCCFrD2VHCCFk7aHsCCGErD2UHSGEkLWHsiOEELL2UHaEEELWHsqOEELI2kPZEUIIWXsoO0IIIWsPZUcIIWTtoewIIYSsPZQdIYSQtYeyI4QQsvZQdoQQQtYeyo4QQsjaQ9kRQghZeyg7Qgghaw9lRwghZO2h7AghhKw9lB0hhJC1h7IjhBCy9lB2hBBC1h7KjhBCyNpD2RFCCFl7KDtCCCFrD2VHCCFk7aHsCCGErD2UHSGEkLWHsiOEELL2UHaEEELWHsqOEELI2kPZEUIIWXsoO0IIIWsPZUcIIWTtoewIIYSsPZQdIYSQtYeyI4QQsvZQdoQQQtYeyo4QQsjaQ9kRQghZeyg7Qgghaw9lRwghZO2h7AghhKw9lB0hhJC1ZyWyE5F3i8jnReRJEXm44/hbReS/iEgpIt/XOvZ+EXki3N6/iushhBBCUo4sOxHJAPw8gO8C8EYAPygib2yd9mUAPwrgY63n3gXgpwC8GcCbAPyUiNx51GsihBBCUlYR2b0JwJOq+rSqzgB8HMB70hNU9Uuq+ucAXOu57wLwSVW9qKqXAHwSwLtXcE2EEEJIxSpk93IAzyTbz4Z9K32uiHxARM6KyNlz585d04USQgi5PbllGlRU9cOqekZVz9x77703+nIIIYTcQqxCds8BeCDZfkXYd9zPJYQQQpZiFbL7DIDXiciDIjIE8D4Ajyz53EcBfKeI3BkaU74z7COEEEJWxpFlp6olgB+Hl9RfAvhNVX1cRD4oIt8DACLyDSLyLIDvB/CLIvJ4eO5FAD8NL8zPAPhg2EcIIYSsDFHVG30Nh+bMmTN69uzZG30ZhBBCjhkReUxVzxz1dW6ZBhVCCCHkWqHsCCGErD2UHSGEkLWHsiOEELL2UHaEEELWHsqOEELI2kPZEUIIWXsoO0IIIWsPZUcIIWTtoewIIYSsPZQdIYSQtYeyI4QQsvZQdoQQQtYeyo4QQsjaQ9kRQghZeyg7Qgghaw9lRwghZO2h7AghhKw9lB0hhJC1h7IjhBCy9lB2hBBC1h7KjhBCyNpD2RFCCFl7KDtCCCFrD2VHCCFk7aHsCCGErD2UHSGEkLWHsiOEELL2UHaEEELWHsqOEELI2kPZEUIIWXsoO0IIIWsPZUcIIWTtoewIIYSsPZQdIYSQtYeyI4QQsvZQdoQQQtYeyo4QQsjaQ9kRQghZeyg7Qgghaw9lRwghZO2h7AghhKw9lB0hhJC1h7IjhBCy9lB2hBBC1h7KjhBCyNpD2RFCCFl7KDtCCCFrD2VHCCFk7aHsCCGErD2UHSGEkLWHsiOEELL2UHaEEELWHsqOEELI2rMS2YnIu0Xk8yLypIg83HF8JCK/EY7/iYi8Oux/tYjsi8ifhdsvrOJ6CCGEkJT8qC8gIhmAnwfwHQCeBfAZEXlEVT+XnPZjAC6p6mtF5H0AfgbAD4RjT6nq3znqdRBCCCF9rCKyexOAJ1X1aVWdAfg4gPe0znkPgI+Gx78N4B0iIit4b0IIIeRAViG7lwN4Jtl+NuzrPEdVSwBXANwdjj0oIp8Vkf8kIt/a9yYi8gEROSsiZ8+dO7eCyyaEEHK7cOQ05hH5GwCvVNULIvL1AD4hIl+rqtvtE1X1wwA+DABnzpzR63ydhNwUWGtRlg6lVVinKJyDg8I5//8SqkBMmuRGYESQZwZGgMwYGCPIjMAYAZMr5HZiFbJ7DsADyfYrwr6uc54VkRzAaQAXVFUBTAFAVR8TkacAfDWAsyu4LkJuefb2ptibWezOSkzLAmXhUFoH5xQOgP9/IY+IIIuiy8zcffs2CLc8M8iMIMvYnE3Wl1XI7jMAXiciD8JL7X0Afqh1ziMA3g/gjwF8H4BPqaqKyL0ALqqqFZHXAHgdgKdXcE2E3JLMZjNc3pvhwtUCu3sTTK1DUSqsKqw25ZYiIjAS7gEYAYz4CC4XeJmJF1pmguxyfz/MMwzC/SjsY+RH1o0jy05VSxH5cQCPAsgA/DtVfVxEPgjgrKo+AuCXAfyaiDwJ4CK8EAHgrQA+KCIFAAfgn6rqxaNeEyG3EmVZ4itXdvGVKyUu7e5hFlKU1ikUgHMKVaBLc8YIpONxlaoMjzPjJZgbad4yg0GQXZ5nGOQ+2hsPcoxyg+Ego/TIWiB9fynezJw5c0bPnmWmk9zavHDhMp46P8WFq/uYFq4SnAsRnOsRXEQASIjm0scmuY/RXXwcxTcnPSM+uhvklfCGQXjjQY5hbpjmJDcEEXlMVc8c9XVudIMKIbcV0+kUX3xxF198YRuXZ77+Zp1PUboQyan2R3IRE4KtGHWZDtGl91kQXhaaVmKjShTfwAgGxmFQ2irSG8xKjIY59mclxsO8ivYoPXIrQtkRch24evUq/uL5XXzpwg72pw6FrdOU1mlDcO1siwJIE4neb+IjObSjuy7R+aiuLb14K42giNJzBrlVDIxgmBnMrMMwMyisw7SwGA0ybDDSI7cglB0hx8jlK7v4ixd28MUXtzGZalWPi7U4pzGaA4DDR3RxXyo7kTTSCylM03wcI7xUepkRzKximHn5FU4xtA5FkN0wzzArM8xKh1FusDHMMWJNj9wiUHaEHAO7u7v4ixd28dRXrmB7zzWbTrROV4bhcQ3Z9SEN0fkTo9QATR77yE2gVb1OXJSbf9yUXRrpGVgnGGSC0inKzCB3itI5FNanXUvrUA5zf9w6jId5NcSBkJsVyo6QFbK/v4/HX9jGF7+yjfO7BfZnDtYhNJ6EyE3jAPB5uS0aWpCcBYQUpguik0R6mRFYpz51Ge8l7vcpTRvkZl1bdloJLzNeeF58BqMgt9I62Cg86yXIKI/c7FB2hKwAay2e/so2/ur5i3hhd4adfYtpmURyleiaQqvEd+A7+DMk/LdyiigkbosXm3Ma0pnaSGc61ZDKVBhXC8+pl6CXHpAZV8lP1R8vDaAq1Vg/qyFSHThYVZQhat0Y5jCGwiM3H5QdIUfk/OVdPP7MRXz5ygTbuzPszPx0XtX4OK1TlwAaj5HsA+alJ8mjOmjSENzFfVrV61xIXfpoT+AkNrAoMvWiq5pWVKsILkqvfYuyyyrZaSPKsxo+59BHrk698JjWJDcblB0h18hsNsPnvnIFX/rKNi7sFbi8M8VeqXC2Kbj6Hkg11yW9fqLgAKSSS+bCFA2Csz6Ci8IT8QJMI7tacvFapSG+OCRBs6bs2lEeUEBDitZq3V26OcwxHGSr+aIJWQGUHSHXwIXLu3j8yxfx3O4UV67OcHm/wGSm0NhpGSK1eL8osquaU3reay512Y7stE5lVsILYkO4j9sqCicCY/z7uuo4gtBq+Xm5GeQZ4NRVssszfz0uXJdq6UWn4fMP69QtZ2AhNwuUHSGHwDmHJ164hKee38aLV2fY2Stweb9AUdSRnGtEdbXG2sIDlqnVRdLIrm5GqSO7+nFbeBq2FT6lWYnOKIz6bk3fGRpre15+/ubg1HdnqprwmZx/nNzb8GFjObL6/AAbV8hNAWVHyJJMJhM8/vw2vnR+F5f3ZtjeK3Blv0BZ9omueQ90RXXLSa8pOYRoLmwnj6Po4r1W275uZyBQqxAjMCpQ8dJzKshCdBdevrpFnDr4SM7Ar7kgYehEjPIAoKw/W3L9FB650VB2hCzBxStX8blnL+H5kLa8sldge79AaZcT3bINKkBTEnMNKlp3ZR4kOukSXdhnXJiZJUR3aEjO1/jCK8MLzt8A05DeIHMAfHRXf5Yyuer6Nhry54bcOPi/PkIO4PkLu/jC8xfxN7szbF+dYXtSYme/QFG6agaUVHTtOl2f8AC0orz2nohvMomHJNnXlF4cblDX6qLwGqJDXOTVD0GA1JJzoSszSyK7mNIMVchKel5wBn7BkngfrrFIhBemKhMRNq2QGwZlR8gCnnrhEr74lW28uDvDzr6P5nYmJaaFrcTWJbq++4MjPKBrAEKaAUwbVirZIcx/qX6geRxwHoVXPwZUagHG6zDiIzmjqKK8LuoAtEt0/r66VLFh2EOsFYb19TgsgdwAKDtCOrDW4okXLuPLL+7i/N7MS26/xM6kxKSwnZJbLp2ZNqwA7S7MPtkByXRh4b+dskOcJ7OO8trCU8ToTqEh4lIAxmcnG/W69i3ia3UdohMDsV54WZzVRWxlPAGwNRpw4Dm57lB2hLSw1uILL17BM+d2cWHfd1xuT0rsTgrsz+o2+74UZvOx11lXV2a6D2iKbr5ut0h4zXtTyQ617ODH16lIJT0TJZe8v8Y0KPpu7ajMi04kFZ3xKVbrIFLCAH6VdBFMwvtujnI2rJDrCmVHSIK1Fk+8eAXPvrCDi/sFtq+W2N4vcHVSYjIrUTrXSF+25dYV1aXSA/plBzSjpyahMaWVzqwGlLdEl04T5relGlDux9vN32LUpwA0NLCkN0+M4jwmzM4Sj4WYESIOJkgPsBAj1c1UC8WyfkeuH5QdIQFVxdPnruDZczs4v19gZ6/AzqTA1anFpLAoXJgeK5Ga6xFdX4Tn3+co15g8buzQzujOSHqLkZ4JjxHWtquF1ojbHOZsp9UB0zrRVOITcShsKj0gK2wV3WUimGSG9TtyXaHsCAl88dxlPHtuFxdDJHd16m+TWYlZWaJMRNeWXC0/L7V2CrOL5aO6edoZQIHUQlXAhWirKbs6usvEJKnMsNxPEuXBJBJy7ffyO0VMeI80uqvvTZCen3jawVgHKS2yzGAa5WeE6UxyXaDsCAHw7PltPPviLi5OSuzsW+xNS+xOS0xnFjNrUbRE1751ia7NspHdUjNmalMQ8TmVN8JxqwoXVkaIcstE4EwtPEUUnamiPDgXxGfaowoqpBKqvwL/2EvQiENpTRjeoDBWkYlFJsAsRHVMZ5LrCWVHbnvOXdnDX7+4jYv7vgnFi85iMikxLUtMrW2kL+cl115x3LNIbvONKYejLcRkYhUf5SXy8/t8FCoCWBFkairhqQo0LPIaxQdj6pSlC5FeFF7VBZoKNkoOydJCXnSFVWTiUNjQqGJ8VJcbn87Mg/wIOU4oO3Jbs7s7wV+/uI0LkxI70wJX9y32ZiX2pyVmzmFqLUqnYY7Ivsiu+Zp9kltWfsnOflpuiJKLr9WO8Kp9YZkeJxZODZwonDFwEGRBdip+OrAsRnUA4FwtPAEs0tRlLb+27ErbTGdmziGzDrl1KEqHWWYxzQw2OLsKOWb4vzBy21IUBb54fgcXrk6xNy2xP7O4OiuxP7WYFRazspxrSokrjvtITuciuS7B9UkvbPSyMJ0ZDkliPa13dkZ4UKmk58fNOVgBclUvPqN+3ksTXsSgjuqAhvAEgE1kNxfRVfd1OjNzisJpJbsovCJ3GFjHZhVyrFB25LblmUu7OLe9j+3QiLIXui5nhcXUliGqqyO6Zidm/Tpd0gO6U5ptfy0rvj4az6+iuaboavmlaU1BmE4FRYhao+gcFIqwwoGEgl0QnqhWskvTmNZJj+wQGlT8auaZU5TOoUiENy0sBplhswo5Vig7clvywuU9vHBxHzszi/3CYn9qsT8rMZlaTK1FYeuhBraRxqxf46DIrktwfdJrHKt3LEciOb/ZFF0qP02OV0PFBbAKOGeRB+nF4QY+ykuEF6I7cfUbOvHdn9Z1RXiCzPjvMMquiM0qYehBbh1mpRdenlF25Hig7Mhtx/7+FM9d2MGlSVGlLfcLi8nMeckF0cVanW1Fc0tJ7jBR3bI1vD76TpXWsXQ7jfaSW5HUJpuDyWvhiXpRwsWGl2ZnZio7LzqE1dGD8MTfl9ZVtyLcmMokxwVlR247vnRpD1f2Zrg6s9ibWuyHel1R+mEG03BfR3VJ5yQWS29OckeU3oLdc9QB3uHyoYpkMrAkygPQWHoIcXaUKq3pIM4E4XnxORVYp8iMX2XBumYq0zpFGWRnncJG4TlFYR1GTtmZSY4Fyo7cVlzY3selK3vYnvqOy2lhsV8635BSRXW+A9Oqq9KWi6K5pSS3lPT6Ng7cPRfEpTsECwQYz5O6xhjFN3O2sVpDdUzq6E5UIc5HcS5IzwgqybVTmdYprInNPm4uuiutw9Bw3B1ZPZQduW0oigLPX9zF7swLblo6TAqH6cyGNJpPXxbWouxJWy4rubbQuut37Qcdm0sGaYeP5boJIwuqKLFQF2p9WVhNAfVqCzGlGWRnWynMdiozMwrn/C2Kr7Sull6I7ga5YaMKWTmUHblteP7yHq7sTXG1sJjMLKYzi0nh29/LI4juMJLrEtzCfb07FjAX5i0gFu5iLjO8mSZrLZQO8CPrAEFcrUD9PJgxwgvpR1/vq+9jKjMzUgtPQ80uRM+xfuesQ+kcrFM2qpCVQ9mR24LJZIaLO1Nsz3zX5TR0Yc4K24joClenLtNaXZXOOyia6436cGjhdWxe0/CE5agvTkPkFqM8AChcTGPWEvSTTSvE+HSmqxpWQipTmrU7H+FJlcq0qigdUFoHjQ0rIfIDM5lkxVB25Lbg+e097ExmmBQW08L58XRliOpS0YUf9VRuS0VznY/ReKCNnT3H5w72+O2QzZrS3o51OrQDO5070fm1FFA4B6BsrpMngFGTyE46JNcUnm9i8Y0wVrUSXWwIYiqTHAeUHVl79vamuLwzDUv1+DrdzDrMCofCJWPqnMJBO5tSuoR2LZLrFVyy/bnzh/t8b7hryRPb1mvvT6I6n9706cxGDc85oCzrBWPDKgoShiuYkAKel5x/HOt1jfqd1qKzQXqq8ys7EHIUKDuy9rx4dd+vMl5aTIs4FVis1blKdHHlAqBbdEtHc8s8bm0fVnApf3mxfvz6u9I3qDlUdJdIT5IXi/NAl85han2EVy8VFGaXCdGdmq6oTvzKEGECbasKq74m6JzCBdHFGh67MskqoezIWrO/P8Olbd+UMiudl5zTKqornfWiQ0ipoSW6jlpdt/xwTcI7iuS6+KtEfF9zZ/NYKrhqu3EQregu7BMNN4EL04bNrIURQW7qdfGi8EyrOSUOPXBO4UxsYEG4BfHZ+eiOkFVC2ZG15sLeFHuzsk5fFn4uxqK0oU7kz7POhZXF+0V3UDR3IyXXxecvAV99R3NfW3hRaHF/Kro6uqujPOMcXFgNYWptYxHYSnguik2SJZCaUV1jCEJSu9M48XaItFm3I6uCsiNrS1mWuLQ7wd4szIpS2KpWZ0MazafTfILOHVJ0XUI7SHifu3C8n7nNFy77+68+Xe+rrkua8quCu1Avq/YLQlOKnwNTjAOcgYXDpCz9HJfG+DXyEpmlacy29KobmmlNp4A6Zd2OrBxOREfWlu39Elf3Cz9DSukwKxWTwsFaP/O+at2Q0pgtZZHotCk6hZfDIvnF49dbdClfuNKxM163Aur8vXP1dhSSurAd9zmFC3XOwjlMyxJFGBzuXL2obUxjVq+TiFDTc6L0XLK6RGuaNkKOCmVH1pbzu/vYL0JUVypmcbYOZ/2PKeof3jjkoE90Pr0Zj2Gh8DQKJDl+I0UX6RReSvwMrpacS6TnggydS4XnF7idlGWyHFIizSA5l0qvfqu5hhWfyqwjO0JWBWVH1pLJpPBRnfURXZx3sSiSQeOhMaWe4b9fdHVkhznR9R4Lt7+8CUQXOVB4keT6K+mFm7aEV1ifzpxZizKujYc6OqsEl4itvoWJp5P0ZXUJtB1ZIZQdWUt2piX2p4VfCds6lKXDtKiXr3HBUC7U62JjyoGiQ39k15XSTIcF3Cx84Uoz+oy3TjqEZ5N7G4Q3sxaTsvDTfmkzlZmmLhtRYprajNfimMYkxwNlR9aSS1cnmIbJnkvrQgpTwzADrdKYql54S4uuJ7LrSmnejKKLPLEdHiS2awtwDm3V7xzgLGCtwjqfyiyc9bU71SqV6Vqvqa338vvm5UbVkVVC2ZG1YzabYW9aYOLUj6krHQqb1uX8eXUaE1XjiT+hea9d9wcI72YWXaQSXptUfpiXTmxkcTZGd0BZKqalxX5RhOafGC333xD/2EheO9btYnRHyKrg0AOyduxO/YoGs8KitIrShq5Ap9UPcExjxjpdrBXp3D0W34f/tO+vJ/mpz2L00t+CiGvst1e/CvvP/LcLn/vENvC6UwtOSD9nGAZQDdELkV5tRIfd2QxbwyGGmiFTmW9waURz/pWj8LpqdBxrR1YFZUfWjt1pgf2iROE0DB5XWFvLK63PRQEC6LhH9ZyF9+E/8f6vLh3rx6sY3fcJDO78NIDu8WjZ1lPYeOCXDhSeKubmzOzUS0t8cfC50/q7mBUWV2cFxnkONVnVhJLe4mtUKV8k4xtbTSqErArKjqwdu/slZlZRlKFOFxsmUHcHAvPia9brUN9jCeGF++MWXX7qsxjd/zsQUwBYPOhaxAvvIJ7cAV57srmv0k0rmktPSKWnCljr7y9Ppjg9HiWiSwqD6JBeeGy1+fqErBLKjqwVRVFgMp1hFqeisn7YAVDLDQhNKah/aHtF1yO8G5G+3HzwZ2FGL17fWUW0cTcvv0R6UF+/m04sdqb72MgH9R8T1a0WXyo92o0cN2xQIWvF/qTELK5mEOp19Y8s6h/vRHya/KL3SawtvK5zPn/5eD5Tfuqz2Pqah49VdE/uLHlicFXazxP3x+/BWuD8zqQxMLyO69CSXA8s05EVw8iOrBX71UwpYSHQpNEklZuIVNFcemxZ4fUdWzUbD/wSsq2nrklyqr5J5TDnVyxZw0tTmen2ZApc2N/Dy06crNKZ4ZTmC3RAz5HjgJEdWSuulgUKl3RfarLqeIfcBNJMXwKLhXcdU5ibD/7sNYkuRk7LdGOmPLWbvkjztnAMXjvSC49f3L6KSVlUX1C38CL+Q0pYPSE+ZicmWRWM7MhaMZmUYdaU5rg6TX5ipbWv6kZ0Swiv59jS03AtycYDv3SotGWdXs0x+ZvvRbn90GovqHqj5G5B/U4UsCVwcTLB5iDHQIfN62w+EwCQJUuiG3DFA7JaKDuyNlhrURS2iuys1sMK0g6LWnz+salSmjg4wrsOEd3ovk8sHdFVnYyHjOJ6X6+1vfAS+sQXvpuyBKazKS7PRhgOcmBw8MrjAkCMj+iYdiKrhLIja4O1tpqb0bpmt18zugutg+39ywivtW/V5Kc+i8Gdnz5QdBoMM3n+H6w2itOFmw37SceJaQ1PFdgrSpwoZrg6zbE5zueeZ8SnKzNppjCNYQqTrBbKjqwNs5mf9DmuuQYA0BAtiHQKz4T96crdC4V3jKIDgNF9//dSoltVJHdotOPhgnRmUQBFabFfWuxNLIYnfHTn63HzL+/FJxB44RGyKlaSKRCRd4vI50XkSRF5uOP4SER+Ixz/ExF5dXLsfwr7Py8i71rF9ZDbE6uKUusZ94FQ9wm5tbQDE6ijukYa05/aEF6X/I4LyfYWHlcFJs//wLGJ7umr1/CktIkFze/JOr8wbmkLXJ0VKEpbPS3+EZJJLTmIT2OaJMojZBUcWXYikgH4eQDfBeCNAH5QRN7YOu3HAFxS1dcC+N8A/Ex47hsBvA/A1wJ4N4B/HV6PkENTKuqVDQT1/JfJr28qvDTakwXCa0d2bVbVnDK67xOLT1DFh148jy9M/tfVvGH/28zfMH/rfjIa4rMKlM6vDF9Yi+1JCWcV7aAt1uqyIL3MdEd+hFwrq4js3gTgSVV9WlVnAD4O4D2tc94D4KPh8W8DeIf4hPx7AHxcVaeq+kUAT4bXI+TQ2DBJY4jpADSjuLhd14ZCVBcaVcKueeG1Hx8Tgzv/pP8HXhU/sL2D797bgwjw9PCHjvlq2u8/f+sSYfs56oCZK3wd1ZYoCoftWQnn4r9DiOqC6IwRLz1jWLMjK2UVsns5gGeS7WfDvs5zVLUEcAXA3Us+l5CliFGdi3W68F+EH9S4B/DNKyJBcojj7QId6cvjrNPVLH6D//min6JF5CZry28LEK1UpqpfjRzw694VDvszCwl/eGQShhoEyWVhHyGr5Jbp7hWRD4jIWRE5e+7cuRt9OeQmJDZFqEtiuyA+iMZhy0F8Go7F6T/qHs2u9OXxi24BIaq77m97wO2gJ/oB/TEC9OvUOefgSoeJdZiVrhKeiZGdwEd52S3z00RuEVbxv6jnADyQbL8i7Os8R0RyAKcBXFjWflexAAAgAElEQVTyuQAAVf2wqp5R1TP33nvvCi6brBsaIojKcvARkFbi0qolXqSWHKrnVE9bOn15vRwYo7rrQTVk44Dboppe9Vrh9Rx881BcVmlWOOyXDmWpMECVwswygzwzyNiJSVbMKmT3GQCvE5EHRWQI33DySOucRwC8Pzz+PgCfUt8u9wiA94VuzQcBvA7An67gmshtiMAPRPYRns+NxfqcINTpBFX3pcTTwvGY9WxEdUCv0a5ParPjfXsaZW4oXSJE/R3FNLFfgdzCOsWsVMzUD/4XAbLMIMsMBhnrdWT1HFl2oQb34wAeBfCXAH5TVR8XkQ+KyPeE034ZwN0i8iSA/x7Aw+G5jwP4TQCfA/DvAfx3qmrb70HIMih8zQ4I3lKENJnf0W5SqVKarXFiSzelhIMLV/peEf/P5mYjknrN7GPH8j6v3ljhi8XvW/3wAx9Z+w5YW1g4VZROMXN+lYo8iC43TGGS1bOSQeWq+vsAfr+1758njycAvr/nuf8CwL9YxXWQ2xsjvpOvigmkruO1IzzpiPDas/cvityua2Algp+4+5X48Yv//OBzV0DnZzvoA3cNEA//EcS/qsVHdlDAeAmW1qF0foiCsw65YQqTHA+cQYWsDVVXX5gtvxIb4MO8SGxaCQ+qrk0giNG3zPdxI9KXBw02XxXXnB7teJ4ifI/iOzGRdL46B7jSwtncd9AiLM1UOs6cQo4F5gvI2jDIBBrHbiHpxhSBgVQt+2kKM9bx4vlVyjM8t5NjEp26wcLj+anPHs8bHyP1D0z9pTkoxPi/P1T9avK+sUhQqsNkVt6AKyXrDmVH1gYjggyoWtnFhHugDtySFGb8v5jWFGnEeJ30eW4VdbvpV/5+b2Ql4ufNvFEs0ZzZeTNSNwHVU4AJ1IUZbtRHchrepVSHq7MS04Kle7JaKDuyNoyHBsZINRVVXBNNErHVItNG9BcjvJjG7BXeMdbwDlq94LhTma8a+/uGsJYchtB3S7PHAkDFCzAOMxAA6lz1b+AcMLNeeDMKj6wQyo6sDcYY5EFYxgAwMaJA1aRSR3JSiS2N8iQRXpuFWc0VpTZVF/+/5HGmMqPclhs5vhzZAOEPEKkiZxOGIQjq+TBFAGt9V2ZpLWalxc60oPDIyqDsyNowGAwwzDMf0SFOQYU6WmtFed3Sq2t3c8I7oIb32pNH/wzTv/n+G5LKfOUI1yS3g4K7PA8zooip0pgigjykl7MwJsQAYao3X8MrrMOstNhlSpOsCMqOrBWj8aCu1YXVrg3iYqBRZDJXs4t1u0p8pim769WrUm4/hOLSW/qFd526Mtt0pTYPSm8ajWP70/Rx8odEBuS5wSA3MPDDEGaFhYsrzTvFpLC4OmOER44OZUfWio1hhtxI1RiBJH1WLxiqjYaJdqQnjUgwvPCStltJdPfCexceP3ApoEPyytH8voU1uyUxAmQZkBsDYwwy4weMx1lScjG+MxMIq8u7KrorS4vCOlhrMS0sdltr4RFyWCg7slZsDQbIK2lJFdWlHZdG6iV90k7BtG6XRneHjeq+agXCU9s9lYkIMLjz08dSu0sFd5SaXZXCHNXNKHk10bMgywR5bjDMvfSyzKBUoLSKWZCcah3dFVYxKx22J4zwyLVD2ZG1YmtoMMjC+mghnVkJL9SJ0qguTV+mNaW28Nr0NaWsqlll+sL3XJfa3QOjaxNcX1ozvo5RX68zRpCLgRFTyS7WRrNMkGcShiD4MXfOKaz1EV5ZWljr4JxDUVrMSoddpjTJNULZkbViczzAODPIgFp2Ia1pJKnlybz4qiaVVirTSLfwOgmy+KoTR/sc5fZDvdEd4Gt3R01nPhCbUpYQ3GHTmgaAyYFBlmFgMuRikJkMeeYjutHAR3VGDKwDCqeYhUiutPX0NWXYLqxPczKlSa4Vyo6sFYPBABujgV/pGsmYrmQKsXrEXbNBJd6yUGMyxndzLtuZ2d51VOEdFN0N7vz0NQvvgSEOlFyn3JZkMPKR3TAILjcGuREMB1mQXYZBWLPOv4+P7mLdrigtSuvgnIaIzwXp+XXwdpjSJIeEsiNrx+nNEQYhRRbH2VXRXUcqM/2/uFJ2GuEZ00xn9v7mdxx4zRGEd2Bn5jHU7+bqdod9HpopTGPg05hGqkme88zX8JzWww3KsNSPg5dbfGvr4tg7V8lvVlpMSsuUJjkUlB1ZO+7aGiDP6ppdtTBoOoWYdEd0jbpS1cXZqt8tEdWlB16zde2fZfrCexenMw9Zv3tgGKK6Foep23XV69KdeQ7kA2CYZRhlOQZZhqHJMMgNxgOD8SBDlizjo86vaO4UsNbPoFKUoTPTpeLTID2/f1Jw4DlZHsqOrB0nN4bYyEKre1aLzgutWaNL63JRcJUAq5pfEt31/X/MAQI8mvD605nA8vW7AyW3gM5OzY7nZACGIz9zyijLMAy3QajTbQwzDPOs8z1c8iFVQ/pSFUXpo7rYrOKcH3Belj6tuUvhkSWg7MjaMRgMcOLkuOrKjB2ZWRx/Z+YnJ86M6YzyusbdtdcWXTatea3CWzaduUh4bdEdSnKHqNkNQlSXZwajLKsWY41NKXlmqiEF1vnUZXxcupDOdL4ZRcMHjvdO/TCEmM50qpgWdUqTM62QRXA9O7KW3D0e47zZw8QYZMbV0pOwwGtXJBduMAZwDjbU8DxhRTZF9Seii02Dh0hrPrjp7794yIlQ4kDzwZ2f7uwMjcKz+69qTCjdF80tEpfOPVhMPC1TYDAE8gwY5TnGeY5xlmM0zLExzDAaZHN/PLRxTqEisKGel15EHKMH+NlWVNX/UWIE+zO/8rlgiOGgO3IktzeM7Mhacu/JHOM8NEPEVKagOf6uQ3R5EuHNNask6cx4O0yzSrorSu8wLFO/G9//u9X2YaO56tA1pjVHIy+74TDH5mDgZTcYYGOYYWOUN+p0C18zbiehrKtSmhZFaUNkhyq6s6HOxxoe6YOyI2vJeDzGya0NDMPsHbFud5Ds0npdlkrPmLlU5sIB5137Wgce3ARevQm8ut9fcxxYvzMz3PWyT3SK7sBo7gDJddbtAgMAwyGQDwTjQY6NENltjf3NGFM1/MxfdN8FJZvaPBTFV8Y6niqmpR+Htz2ZcQFYMgdlR9aWe0+OMcrFj/PKWrJLuy1b4+3S6C4P0Z2v6ZmG5HoHnB8m/RfOffVGfVvEQfU7CFCc+jSu3vOJ6j0OjOaWkdyCup0BMN70Ud3GcIDN4QBbgwG2NgbYGg8aEV2v8DDvvL5Laogv1PhmpYWzDtY51vBIJ6zZkbXlnhNDnBgNsTebYBqiuTwTZKUXnjWCTH3DRBRbGQpxMaJT5/yPtXNw7fodEFZG9/W7sAbpHL1RVU+qMy6i2vsal9+LS6cfA7Ki+4NH4QHYPNc9qfRBkd6y5wiA8QgYjYHxMMfWcIATQy+5k+NBHdElJktXhE9vC69FtTorfewUyPwqQVV0N8wMZuKwrQVOKTAasoZHGNmRNWY4HOL0ySE2BmHS4Xz5VGaM7qooL6Q109lVum7LcuhaX2v/5vm/v1hEQXizE/MDzpdJaS47qHwAYDwGRgODreEQW6MhTo4GODEeIs+ySmJxvGJ1eSGyPvAa+o63DsahCgirJ8xKP5/mzpQpTeKh7Mhac//WBsZ5hsyIn2E/RneNeTPnhyAI5oWXpSlNY8KMLK3bEvWnvn2HEeDo6kMYbL/lQOFN7n50/qWOkNJMGQLYOgUMR4IToxFObYxwajTEyfEYg6wZ0aUD9BuXWA3vkIMFmI7DS/bVQxP8Phs+SJmkNPenFN7tDtOYZK05fXIDp06McXVWorCKLPPDEPzUVQqXyC6iyXZMbebGAOogIaUpxjR+aIH5tGaY7nGOw6Q6F9Wtts6/F1fhI7i+PKDmlxvP6WOZlGXKEMDmSd+Becd4jBPjEU6OhtgcDavIOdIlvD6nSeu+93rrTCbiwyq9qerTms6no6eFg6qX3caIP3m3K/yXJ2vPS09vYGdngr3CYWAMiswhsz7Csy5EdxrqQfCCS2fzyI2BOj+mS8J+qw5O43lhwJ1Uj/xrqZ/+ao7DRHoHCPAg4Ul5x4HR3GEkBwTRnQA2NgSnx2OcGo1wamOM0TD369YZCbW6+YgubVCZi/IOeF/tOachPlWISFXLs06RGaC0FrszharD5rhj8CFZeyg7svbcdWKMra0RNguLmXUYOgNrFdYJSuMHKjtX/4xaoNGsAgCZGIjxqbE4wE4VsJUUnTdd0qlpnZ8UuZo/EstHen0HuuTUKzyXQWWK7a96GFLegdH5d2G4+9DC1zqIIYCtk8DWpsHp8Rinx2NsjkfIc1OnisPA71pqYbtLeGhesjSed7hrA9KGFW/AeB8nnd4rHXR/hq0NCu92g7Ija48xBi+7YxO7VyfYLwyKTMOsKgZ5piitj+7SsExDrc7F2fiDreI+depTmyGlaQUAXJXCBACI/+GOKU3ncKT63aJIb/PcezHbfxUmdz8KzS9D7CbUTIB8358zuIzJfb8DABjuPnRNohsBOHESOLGV4c6NDZwYjzHM82oM4zDURE0VzdU1O9MnvGC3uY7MA0TX152Z4tR32TpVGPh/g9I5TASQaYlNpjRvK/ivTW4L7j45xqkTG9grHabWYTgwYT5GPw7PuWaBzbTCCg2Si9FejPQA/6MqVQTnOzxEAFhAMv8DH3+PnZ0fIN3msEMV4v7h7kNV5Lb9qg8BpjUnmSkwvedRDHYewmEZAzh9J3Dn1ggnx2OMhkPfuGMEWViQNQtRcltyBwkvld5S4xA6iOlLDfU6/9i/llPfOGRVkavAOsWk9P8Qm+PB4d+M3JJQduS2wBiD+09vYHdvgv3CorSKWeaQW9+oYoxAM/jozs3LLm1YiVFevDeqgAPUBPGJb3/P4KO5KtqT2EhRR3vHFemljSmNYz37+xjCN93cfy9wcmsLm+NxWIjVT+qc58ZPuB3mrYyL5jY6MYHFwmuNuwOWc14UXPzM8fy+hhUJac04zdjMOci0wHiYL+wCJesBZUduG+4+Oca5q5u4OrOYWsXQ1rU76wz8z2QYVK5Jdya02a2ZbFdpzXAszuEoAtiQ1vSRBiCulp2EiOO4Ij0p74AOusQm2Hmdr+ENz78LwwVR3hDA3XcAd5wcY3M8xiDLQhdrMvYw3rK6AzOVnN9uzppSCS+9qjjsINTbnPohBPH7AvwfFY1/B62HKhg1yEz9miarU5xpw4pTP2G1T08DM6eQwlaTVJP1hbIjtw3GGLz81Aa2dyeYhIVAC6sYOL90jP9RNRAJkYCq775UQKC1yFrbVn1niqCWHRwgma8XOectFxtXVH2Pi7U+vbaqml56YHT+Xb5GZ4rmMfEn6OAyprGG1xJeDmA8AB542SaGeV5JrhqLKII895HdIDMYhls7qvMrSvjXTB8DdXozig0uisqvSqEhtls0XKLxsUP0pgDUAfGPFqcGuQFEFJnJkoYVhVF/XziFlA7D3FB4awxlR24rTm4OcfepMWZFiaIwfmHQTDDMmj+szYYVaYQhCoE2fhP9cQlRXprmdKpwIrDqoKKV7FyM7mJzp/ruTU3Fh2sfkxdrd9N7Hg2pS6lEV3/IArN7Hq1kZwBsjYCX3reBQZZVkksH18c6XZx6rR5q0Izq6nRlc2A5pJYbTCLDjlubZTXUbl6pVgoqbRgS4ZuT4h8lRhSFdRAAAwpvbaHsyG2FMQYvPbWFvf0yNKv4iYR9a7ofFCcA1IQ6kvN7NPnBno/sAECrSFBCpFfJLvyoOqMhagwNLVItnVelNtXUQxXS+TYbEc6SQxKGuw9VzSg7r3u48/vQ/DI2xsAdJ4DN8QjDPK+klrVElwXR+Tpd3ZCSNqU0m1PivUChsM4vsRSj3Or7jF+6LJbdtdCeR1NUUVhA1UHUd5Bap4ABSlWIdRj0rKRObm0oO3LbcWJziHtOj7E3KzArwyrZVlHmBijD8IG0Zqf+51eqqK2Wn+mRX9xWAE5d1QJfneN8pBdTmk4B42rRxehPe/b560I9S8sBkV5/DQ+49IqfwWb53bhL3lxJrV7pIZ0uzXdexinX0vRl/C7Se2lJrtGhiWZNz7Sk1zlA/sCpxKSxXTWvxE7N+iDiXAAlnJcgBDZcg3HamAGGrAeUHbktuXtzjJ2TBSazEqUzKK3/C19gIK5ORYYwD1BNfoPFp+A0Cs3LLwpNFLBQ343Zkp+/ORjJqm0NEWEUaBrZqWnui7f2vjlilBSiKLP7Lmyf/h2oKebOK+QSvjz4OEaa46Xy5kpy6eK1eWaq1OUw1OuyrE5f1rW5+MI9kmt0YbZSl4n0VqmaSoNpejM8itHezAJOLSQ3KK1AYKrPRtYDyo7clozHA9x7YgPTaYmJ28fM+uhO4Ws3mvnuTBH1Y/CMSwaLByFWr+bH3MVBzHDo3DZVpOfHegFoRnomNlr482KkF4m1p7nuzWRboim0XnfPGOCU/F1sTgzOj/89Crk0ZxMnBb4kv4dXyDfOr9YeanKDKDkjGITZUrrSl/Fz9421awuvnb7smhC6s4a3KNJLnpMGfWl3ppG6scU/RzC1LkworRgKOzTXCcqO3LbcsTXEzmSMvVmJIqQzfe1MoYhdmH76MJeFWpMKpOqu9KnNaMEY2WXG1FFfyFOmTSvWASaLzRPaGMoQz0vlB4ThDq55/anoYgQXaa/kkIlgS96MlxffiLPDf9b5fUxxEX+oH8DY3o3XZ9+LV2XfVDWjDPJacsMwvi7W2+pZSvwFZFW0t5zw4naM9I5zKZY0yvP/1pJEef4zxHXxFILxgPW7dYGyI7ctWZbhrhNjTGYFitKidA4utKEbZD7KcwrjfKRlxU8TFnKWPmUJQEMkolUk50Kzoe/ClNDiDqB+DlBFcNY5ZFktOoUiS2pMGoSoLQukdal0sHa66rok0oviG+EuTHGx93uZ4AL+q/0V5Ebw6sG3VEsiRcnFKcHi+8XP5K+jfb9YeI1oDklk2uKwAVarhNe5vyG5KGwR2PC9ltZhCnBIwppA2ZHbmlObQ+xOhpiEJYDKkF4sRDFQP6eiAeBCd6YTAGIAq5Aw93MWl/sRgVPXiOxE6pXQ42BzgWkMQjdBlPU+1xraUIswUkWeUTrhP3H4diq4tvAexHvwBfw6HGa934vFDJ8rfxuv3fxWP54udmGGJpV4FU6bacyuMXZzkkNzfzt92Z69pvp86XY7zbmkjBpNrdF2Wh+M0Z2qT1NbVcxKV6Vvya0LZUdue+49OcbVqcWstJg6hYiFFH6Edyl+0LEKUIZhCn5uMC8+UQGsgwZhSVLPqyI746M0EyI0m/zkWnX18kJJTS9uAyGFqTGSah6rIyKpIr12s0dbeC+Xb0SuBk/idzFZEOHt6YWFogO8rNpdmKYlvnQf4Mcwtut2Iv6Pij6dHCWyanRqtppUIvUfHf5wPQ7SD0swqlCroeGG0rsVoezIbc9gMMBLTm36geZOQ1cmEJKNYbS3H2JgUEcg1oW0Jvx6d40x6In8nIYf3ER+AKrILwZtVf2uFda1ozp/XXVOM45jgzYjI9PxOG6/3HwjHpBvghHBp4r/Efu40Pnd/O7lf4Y3nXgf3jB8ayW6qgEFqeDQul+ubhe/z7Ru16a7OaXzcpeiSv+m4gsPmzOs+HqtgYR18QSlU2QCdmreglB2hAA4vTXE7myIaWlx0rqq5R6wML5EByM+zSkOgKqf6ssIULoqyhMXIrkkrRnrev6H1ITIoU5tptuKujsQQKf80lodkAx5EMzJLt53SS+e+8b8+/Bn5UdgO9KaV915/OedD2OQG3z11rf6cYFpfQ3+u+lOX9bXd6Dw0Lz2lGUFeC3UY/BS8bXm06yiO3+NVtXPvELh3VJQdoQE7ju5iav7FkVpofA/2k4VhYmTHMeGFf/YhShQM4GoQF1IY1ptyC1W49TnxQDEaAyNbQCV/CLtbS/DOlUJhBSgzgsu3jejrGZ60xjBq/NvRpYJ/mL2W9jT+Qiv1Bn++MrH8NWb39pshGmlK9vpy4V1O0FSp6vTmF10R3btml37eM+LoTUsIT5uDfHIBJUA6y5Nvx/wtTx1YB3vFoKyIySQ5znuP72JWWFRuClEckCA/ZnDLNZxrAvjswQ2SWc5VTgjsNb5SMH5weZpPU/j4DfUHZYStxFmWlkQs7Qjv1SQxjQbP+J9OjwgfRxFJ0Hkrx18K75m6634tYs/1PneO/Y8/tWz/wAns7vxTad/CG84+baG4NppS7OE8ExHGrNNKvVqX8eJcQ7MduNP5/eoyR8AKtUMMALn/z3Fb2mSfo4NKzG6i+8BCu+WgbIjJOHU1gh3zQpMixJ7hcXGIIfAIiv9j/Gs9C3pAGBjfU/8UkFV96VVqJFGPQ8icFCIq6M8RRiwHjBJ/Q7oTmE6bf64KhS5NKfrSptUukQnYeycCRMx++m/fLfhCXMPdt35nm9HsWPP4w8v/SLECL725Nsa75uOr4v3XcJLJSchRdwltUW4MBwkNuqkM7Wk13AQldDglxRSB8Ap1CiMMYACg7y1RFA78qbwbgmOc/wmIbck95/exOmTG35cWRDBxjDD1jDDeGAwGmQYDfzq3MPML146CJMjx8d56GDMjCDLDUxc4DSdezLdNqYxPm5uUHichLk9u4nUU3ulA8jb81v6xVUNstxP85VngkGY+msYOi4HRvCWUz+IXIYLv59Sp/ijSx8DkAoOjfv2wPIovDR1aZDWELvfy7TEYsNsM9WA/cAqOiTbM9OoKkp1KKzDrLSw1sE6N/c8P1FAu4mI3GwwsiOkhTEGrzi9gdm0wM6kwDgXTEtgNMgBWGTiQkTiYI1CrG/YdE6ROaAofWTnYs1OfZ3NmRDZqZ+ayjlFMj1xldKMzNfvUHVyCkJNSerEZzuia9TSQuQTF1v1YqxveYj03nDybTBG8MdXPoYdewGdM0zDpzVT0VUyk/4GFQmp0yp1aZq1uzYhm+jrY8lltKOozucukF97YHm1v+ucJJVZapjFprDVXKERp75xiU0rNy+UHSEdjMdDvOyuE/jr8zuYzEpY42BVMcjDem4WyI2gsApjgNIiNKzAT/oMhDSmg6pUM6U4kSpqiNryAxwU4gSNsXVojukyrUimPdYOCOk7ad5Lkq6MkmusMp7MiAIAbzjxVnzN1lsBAL/y/D/Fju1Kayp+4a8/gLfd9cP4utPf1tmY0lWvE6BKn9b1w/lXr2aTcS2xHfxPd2CNr9eDjTF4zfF4groOaNVP3G2tVtfPLs2bH8qOkB7uODHG1WmBc9v7UIR10ABAgVGew4hDJg6ZEczE+ZlSnEDgEtmF5hOnMGLgqsmeBWVcNVzrVbbT9vfKbRKjC1OdoaKAmnrmlBjNmfo+FU4UmqnkVu9rdHVKPQ2aCPDNd/ww/uDiL6DU6dz3s12ew/977t9AjOBvBeHNpS1lXmzp4xjlRVxoNAH8HxNtukR1LG5JU5rVGLz6UD08QatJpTMTRbiatCpZLZQdIQt46Z1b2C8dsAdASx8BqY/ynBGYQYbM+Q7N2LASOzedi2PovPhsrO2owMF3dnr8oGWXZDFT2UncjmeL/+VVJFFdh2CieLIkmotpRJPMBBLFFn+8gTo6ef2Jt8KI4D9f+nXs2HNz30+pU/yn8/8n/vYd39bbiTmXvmw1qABNyQH9Auseh7dk5NZF3ySa6SlJlFePj5xvVHG2HmyfZ5TdzQZlR8gCjDF41V2beMo6qAKz0voVEBwA8VOIwRiMBpJ0Z4bZVcJAbw3iq2Xn17szqIcdGKPNzsCkD0JEqq7NWg7SkFSj2zGKJdynokvrY+nzEa6kcTyI6w0n34o3nnor/uXT34uuGt6V8hx++q/ei9P5PXjnff8IfycRXyOyM7WE42MIqvlIG997Xw2vvW/J83ppPb9r0dfuKK9bkk4VM+v/LQc5+/9uJig7Qg5gMBjg5ac28BwATAGIAiWQqUDVQULkZkRgMmlEecaIXycviC6VXZwjU+GbMFxcUQGANfUvrF/xIBmBJwrnWoO4UQslRnPSEcUBXZJDkHI6Dq4+Hh+fyu/Bdjkf3cWrvFKewyPP/ysYAR6689u7hZc8BppSjyzTmZle9/y+xZHetcRc7ZlW0lRmF4VzEAvkGYV3s0DZEbIEJ0+McR8EL27v+bSiK2E1LMSqgEBQAtU0YTCCzACZc8iNT9FZ5/w0m+rH5JVSj4IuwtguhNfIksYMl6TRqpSm1JFb1YiCWnZpLS7SJ7lq2rJku+vxt939D/H7L/7rzvpdpNApPvnCr+Lr73579fws3ifX6+uC3VFZX1TXJZb2uUeO9BK6ZlpJu0Lbqcw2hfX/cpkIm1ZuAig7QpbkzhMjFLas6mmz0gIimJV+ALJTv+AnMvjlf8Jf/pkRZPDpxDgY2i8FVM/8EdfMizOrpMO5XHi/KDQbaoEiTdkBHVFN+E9bFpLsi9XAtKW/K7L7utNvgxHBfzj/a7hSnkffsITLxTl86HM/iu962Y/izN3f3ojsYqNH36DvvnrbMhFcH9e6HFAfVZR3QL3Pd+AqYMCpxW4CKDtCDsFLTm/BqTeRCGBKQI3zDSkQvwo5BEWc+suEmp3W9TB/X6clY4ozqiOuUo5qu166B/BiBPp/tGME13VOGq3FbSTpy0WRnQnC+1t3vA0igp974sdwpSetebk4h9/+8v8BI8Cb73l7NRTCz5jWfd1xkHnX5+na3/Uy7X0HDkPovBJPYzxea2xe+CgHEldPcADUKtOaN5AjffMicpeIfFJEngj3d/ac9/5wzhMi8v5k/38Ukc+LyJ+F20uOcj2EXA/uPbmJE5tDbI4HGA5yjPIMuREMMz8LSSZ1B2RmalHUY9piOtG/nkhzvbl0VpVq9pSYlpT6fKAWm5G6/T3W6JqNJmiMpxNpTrElreOxkSR9XD/Hf4Z33vePMJBR7/dU6BS///xH6+hxgeiA/mN9UV2XyI7UmbkkjfipMQ4AACAASURBVJlW9GDlxdlVFEBp3VLPIavnqJHdwwD+UFU/JCIPh+2fSE8QkbsA/BSAM/D/3o+JyCOqeimc8sOqevaI10HIdSPLMtx/cgMvIMyaYgQyK2EsgDDOCvBzLYqaMEYupCOjrFR8o0vyY1wPKG/S+MHWZoqxi7rJZP6Hv53OVO1PX0pLeHF/PP+hO78NRoBPvvCruFx0R3iXZuegoV1/URqv71BvDW9JibVPm2tWWfJ1FN0pyzC5ygHPTRtcfPdpbo6eTiWH46gx9XsAfDQ8/iiA93ac8y4An1TVi0FwnwTw7iO+LyE3lDzPcO/WGCc2BxgOcwwHOQZ5hkEuVZSXG3+rZyox1Tp5XVFeO8JLx63V263OSmlGbtXg8URYadTWrMOhI9KT1vnz++vnC/7uXd+Oh9/4K7hjcG/PN6X4qT9/P85e+FTvd9mXvozv073/4LRmZ6TXdYWhE9Y6Xz8trVt4K6xFaf0kAqXzzz0oWmvPnVku8RyyWo4qu/tU9W/C468AuK/jnJcDeCbZfjbsi3wkpDD/F1nwp46IfEBEzorI2XPn+tqfCbl+DIc57t3cwOZogM1RjtEgwzjPwsTKppLdIAxHiF15qeSMSdOXh7+lEgTm5dcWn0nkBTTPu5ZUZrz/rpf9aG9K89LsRXzsS/87/vT8vPAEC9KX6OnA7Di/S2x9RLGV1oXJnbWaWHoZASnikkwKGyYYKJ2iCDK0HQKM0V0KhXd9OTCNKSJ/AOD+jkM/mW6oqorIYf/lflhVnxORkwD+LwA/AuBXu05U1Q8D+DAAnDlzhv8LITcFo1GOewW4IH7IgTECRQnjJ8lE6bwoZkDV2eBcODeZDgxYPo2Zbkt1Tn+qrzOdKQDcfLrSn784ldneZ0TwDXd/O4wAv//8R3Fp9uLctczcFL/3zEfwpnve3ti/KH3Zl/ZcNq1pku+1b1GCxve2bL1vwa9PQ2oKCLT6g8CpzA1VYErz+nFgZKeq71TVr+u4/R6AF0TkpQAQ7uf/Vw48B+CBZPsVYR9UNd7vAPgYgDcd7eMQcv0ZDXPctTHCyY0BhnmG8cCnNeMSQQPjl9BppDXFL7MzF+WZVhqzZ9mfdlpzmXRmuj/WmtqR3EGpzHaEV42fE8Gb73k7fvpv/yr6ehwvzl7ET372R6oIL5P+H/ne9OUhBKghJZnOZpOyCsEcFJjFCLB0Wi0T1IbLA10fjprGfARA7K58P4Df6zjnUQDfKSJ3hm7N7wTwqIjkInIPAIjIAMB3A/iLI14PITeE8WiAO8ZDnNgYYGM8wHCQYTxMhJfVKc0skZ4JkssPSGealsgWpTP7xNeu39U1wcWS6xNelQJFIk0B7hz21e+88H79iz+Hz5z/VK9s+tKXsZ7ZJhVgOlONa5noILW1j6861lL4yQOKRMBxf9k1lQxZKUeV3YcAfIeIPAHgnWEbInJGRP4tAKjqRQA/DeAz4fbBsG8EL70/B/Bn8NHeLx3xegi5YYwGGU6Ph9ga5dgY+YaV8TDHIDcYZRnyTDBManl5j7hSuVXiMsn6b+n+DglWkZ9pyizuV2hzn2mtgoCex2k0lwivnuTZN5i+94F/jKHpH5Iwc1N84pmPdB5bnL7sfj0fqWojimuf21XTmz+n95IXcpi6WxRwjPbi9VYD0MmxIbdigfTMmTN69ixHK5CbD1XFZFZityhxdb9AaR2mhe/emxUlnPp5E53z04AVYdJgm9SV/Erc/nGcPcUvAdScWQXoaKVHf/0uHos/qum+uF2lPeNxkUYTSZXGTB4b8U9IZ0b50/Ofwu898xFc7KjfpVd71/Be/L1X/mO85d53AOhPbcaZaLo+FyAN4XSd27WvPcA7bdxpH28/Pz2Wm+WbY/rOj99jbgynFmshIo+p6pmjvg5nUCFkhYgIRoPMz5u4KdibloAIstJCRFBaC2MFhYRFXdXPn2jUN7Oo+lXOJZGaAl4mWCIVl/ZboCkOgV8HL/3RTqUWn98lvDihNJBILuwT1EvfRN50z9vxpnvejp/87I8sEJ7i4uxF/NrTPwcA+OaXvONQogszcaHdMdLliq5Ir+t9ko2eaz46XXNqxiWOVB2GYtiwcgxw7hpCVowxBuNBhkG43xhlGA5zbAwzDPMMw0FWNawMc5/WHLaGKnStJn6YoQhdNTkvpjoqq5b+kVYaM6ZHjW+kyWKqUnxU0t6XyrLNew5IaQJ1WnPZZpWYsuxqi7zmmVXaxxcfPhKLsmlWFdPCckjCMcDIjpBjIMsMNod5iMoERhwK42fQ8B15OaxzMFZRwlUrlbswj2JcEMEBVerSqS4M7TqjlfCfKspz6XYr8ktSlGmE1xXVRWFCupfpibzpnrfDAPjdZz6Ci7Nz6OvbvzDtHjvbXp6o7qzsaVbp+H66RNV+7qHqd0cUXzqjShcOwKx0GGRMaa4Syo6QYyLLDLaGOTAtAQTBQDArLUQsrDMojU9rWucgElKZEEhoWnAA/Cxk/T+OkXh0kfRMmrKs9jXTlp3C69nftfBq+p5GgLe85B14y0t8Te4nHvuHnWnNu0fzHZyp0GKXZaS7ftddO7uWZYAO2j4qB00zZlUhzsGocPLoFcFvkZBjJEZ4sSlhkBuMBhk2RoPwOEeeGQzCzc++4sfgDcyC7k2Zv9Xj9HqGKSSpx5jCzCVJTyYpyvg4TV2mqU4Ac639KXE4Qlssf++V82nNoRnhe1/5Txr70jpdW3T9Y+069vU2tTR3dKU++7a71uE7LMutmMDJo1cJIztCjpk8N9jUHCIW08J6WRhgNAipTBE/9spaFDa2ZFqoq6em0tilEn8mF/32SeOu+qGusqBpWhOYWxOvGkbQEcmldAV1MZrri0Jj1+UnnvkILkzP4e7RvfjeV/4TfGOI/CLRUa41Xq4vfdlXq+uK6ub2zb/Y/IdaMQethQfUywNBBKXT6o8Vcm1QdoRcB4aDDCKoVkCYFRZGAA3RHiS0sxelF5sqrGg1Q4drOy5dZiZ5n8ZPYSI9Eb+iOq5BeF10RXVpbW8R3/ySd+Bb7ntn7/EYEaYDr4FFww8WDUs4eN9cva79HLSFOf/+h2XZOM06RR5aN636gi6Fd21QdoRcJwZ5hk0FVEsgR5hJw9bRhtTRXmYERWn9cAOrYX5HnZNd34+mtB5UVbkouXqzjuCwWHApaVR3UDR3mPP6RAcsGlTe8V69DSxzpjr0YPOuy6hmQ1nij5D4esvUYdvNLBTetUPZEXIdGQ4yQARXp4VfvTw3KKwfo+bCzB8qApPV9SojSc0qCNAfbNzN0RYeXP2j2SW8ZYlR3aIUZ9e1LFrPDlgsuq76H9DflNId1XXJb/6Jc+nQ5LFfGUGqVRKuKcVZ/bs5ZOH9Fv2hkUZ3AIV3rVB2hFxnhrmBcxkAxdQCeZi+a1bGHv44w4n/MUub8YZ5BudcJbhl+xbib+hhxdaFnyllOckBB5+bLhfUtTZce/hB+ry+tGbX+3Wde3D9zn/eVG6SplaPUNBTBVTqqND/o2q1YG2UX9dQBQrv8FB2hNwAxsM8/LhZTNVCVTDITJWqdE4hcX1sEVQCRJ12rHYsgZ8P82g/jP69tT+f2HH+QdFcKqx2MwqAatB7F30v3Sm1Hin2vWa8FgnNRNVzDqjfrQJFkFlYIsiIwMn88kBW1S92fxwXsYZQdoTcIMajHJA6LViKgysVuRFYAEZMEu15Eu9dF9qpSrvEey+b3mwPLziM6PrSmn1SWybSi4P6NZmA9Frqd8uiS0wAF8VnS4tRns3VPLke3vJQdoTcQMbDHM75hV5RAsjCQO3QuFL9IGcGZai5WWcP/T5LpzvjfYesFLrwdQ5Tw2uLrD3j/8KIbkFac9mFXdPnx8VdRboWzl2c5ryekimtC1O+abVaBTBf0yPdUHaE3GA2x0PItMSelmGPm29cMYIsHNXMVMvCNLr/DjnwOBXbMrW8vpc/VKNKUp+LtNdy6xPdH73wB/itL/0yLkzP4Z7RS/ADr/lv8C33v7N63cPW76Lk4veWGTP33PYHbU4WvfizHsRhx4k7VRhIFe05Gyb1Dk09B6WMb3coO0JuAjZG/v8V91CGdJ7xKyI0UpnNxhWp/hOpU4JtxCH8TF57NNIeRH4YyQHdEnPJi3aJMPJHL/wB/t0TP4uZmwIAzk9fwC99/l8CAL7l/nceun7X7vjsSn/Oj79rR3nXVy7tf1WFzwIIvOgEnEtzEZwujJCbhI1Rjq1hjmGe+Ugu8+m6PDNhei9TCWERsR2/73YUJNzilGLLRnNd56Z1unhO3/X91pd+uRJdZOam+I2n/+2h6ncA5js+O6K/ru/qUJNFHxOuY9qaKL0ZpxVbCCM7Qv7/9s411raquuP/Mdfa55zro5U3VKiApRqaNDS9bUxbTcvDS5sqJiWpkehNsRroN5s2Skj7wbQR6wfSxARCKAppG6k2UWKaS/Cq1aapLVoUJNJ7RRtA4F5QWqCXc85ec/TDnHPvuddej7nWmmvt1/glh7sf63n2Yf73GHPM8V8iJnN4YIwzArOxe2EijBLjMKBIgZ1TwgCDG5FJ+ZUZq1btV7WIfGoiW5yCdChFpa4Iz++eCJ6/yzQXXk9hT825ybmi/YZXu6pPW1t7oO1RIgUrBYjYCcKSMZnDw9guKifwWE9E0IiDmvRXzM/fdRHASYo0N49XZeNTdIy6Tilac9B2LvV55vbZeG732bn3z9g5u/D8vni6RtKFkW3RnF5hVBexDLMDmhlJxcnFHqgcSWMKwhJyYDvFq0YpRokCKWPomiYKaaJm19lhGsU4dwO33cTQteKndJ/cgB8ioORdQ0gbrKrt8qnP37/4DwvdEt7zxg8U7ufQeip0xYvK589dVNg473+3OCGp+ywyZvNTYb+0iUhkJwhLyoHtFMwaIMIeAA2NFAqAnlmDVzb4lReyxMNES+GD/yQlWhF1FAmTq7q89/E78fzuCZyxczbe88YP4K3nXTWznb+bX4QS3D+zIKorXnze/Rfa9hB1XniAndtTANuKTUlritgJwlLzqp0t4JV9gM0AuzfOoDRBw7gkmF6NUxugcuFD68XoeecBIKzx8+z5p4U1VQvTq9bX/ca5V+JtP3NV0EJzfzlD2bq8IsEtjupmny9aOEI+xsxbpiALzw0idoKw5LxqZwQiwsv7Y2jNUCMjegAwznia1rT/iTWH5w+OScvBMj8vV5Zaq5u/q3u/TOhK3Q9K1uQVRXUhBS1taLt0IcQLD8DM2jsRPBE7QVgJ3Do8ANjdz5AqBU3O+mc++ppNXdoHViDrXBLyA2KT4hT/GvLixDzf5BmojubcsarSnlVC171RdH7n+IJRZA+Ux52yybnzxSybLngidoKwIkzm8MDYJ2B3rJEmys5NhaQyvTm8Hqibv5tblF4jYiHbNBU6lBWqBEZ6XebqZlPOpiuO5sDj5X2dMP1SQxUCrDXP3NsmC56InSCsEG5ZwktsHM0zAjLWM/50gPNbM/t0XYtXJ56h83fuOCFLDpqmNYOEDsVzcsGRXoNCHACFUXfucK1xgjl5UmANBExbjPlsquCJ2AnCiuFSmi9bP7xtmMHeH1id87nBzeOpmQ4cocsJ/EExv/4uFOb6dKV/zibVmlmu5Vip0DUoVCluHxYm6FUC59NFbIpczousgYiKt93E5tEidoKwghjBY2A/wx5sMYIX1elcl2gi44dGMytrZ+fT/G19ODTVVoBLbRp/vvqUZUjEl19H50eMZULXpCKzeElCtdgVefFV0rPO+MIHMEZJbv4U5gtSmmzOUmsRO0FYUQ5smyrNl9iMnfuZNnM0VjTcWOeEr2r5QcyUVn7uri7KaZvWLOqtWYRShK/96AHc89jteO6VEzhz52y870034PILDjVYklB87HxD6VC6NpEOq8c0jLUGkZqzBmJgo9wSROwEYYVxjucv7RtRU8STAdgVMKiZiK96IXpTytqL+ZRaA9UUtPjbFQ3ILn1Z5333tR89gE8+fAt2bSPpk688i08+cguUIlx+/qG57Yuvc/b1xpHc3DFb79oKV6jirIGc6GlmkC5vmr1ObE4MKwhrys52itdspdixDYBd+y9lvc4ciqYOCkXtxNygXvXjtxmrai/mkxcFv61YyJKDKqGrOoYTwXseu30idI7dbBd3f+/2ufMVpS/zaVMzP9rty8LQ7cY49zhjcx9sW4ttAhLZCcIasLNl/ldmMPbHplhFAXPpTEU86WAytxShZPztMmc3OXRgFBey7dRstVxg/feee+VE4TYnT3mNpak6femaSUdhAUFUUQNp111Fkfm7WPd0pkR2grAm7Gyl+KntLYzSZBoR0TRiSb2Ij4h6zaVNIkHC1I8vIF2pAiK+yibSNN+I+swCZwQAOOvAOXYfIFXzQ6G7B9dMOhaLaiJdlrrWzNgbZ8jadA9YIUTsBGGN2Bol+KntEVJlzF5TpWYGf5MatKlM5RvDegLYYDDOpzjz6c2QMCY0remuq3QbK+h5EXzfm27AdjLrmLCd7ODwm28AUFJ8YiO9GCnLgkMvhKrbYAB7mV5rpwRJYwrCmrE1SvBabOHlvX1k2jReThiYFqnwzMA3315s+sR0/Sgu2gihtDXZRLgCjmWXTQDlTaSrilR+6/xDUIpw9/dux8lTz+KsA+fg8JtvwOXnHyp3OcfsQvVohNxv2GEaUydjbjkCM63lkgQRO0FYQ7ZHZobm5f0xyFZnuk74W2kyiVjq3M6LKhGbMNezs6Jqs+Dkdj7JbFsYdbj1gyXHc1Hj5ecfmqu8LBM6wnxrs1gs3gev+vzub2Sc6bWzBhKxE4Q1ZWuUAABO7Y8xZoCYpqLnBjEb0vlOCWYdcsfR3i17UBQubt6+eQEraiJd15Gl6v2ygbzvJF4ssWsjQqH35tbejW2jgnVZliBiJwhrzNYogSLC/+2PodlUYhITmBh7Yz3TfWQ2lUnTbvxQlZ35Z3QT04GYmcEN5//K0pqaw7ZzF1QV7VUJXZ/D+jJESUWtw/L4lZsZM7Am6/BE7ARhzUlThVdTilP7GcjO2Wn2zV9tP0fwzNd/NygmBExmrxqMebXBobfovWoAnkR1NSLmrrm0hJ5MMUqZmPY9nkfzwRtAd3y3hHURPBE7QdgAkkThAIDdcYaxtt3xobGXmZ6Vrq/mbDrTNhUmxMnvkWdJgyZFLlRvHhs4d1fEEELXde5zBp7aBdmn5eednN9uyxQklnm3hHUQPBE7QdgQkkRhhwh7mcY40yClABAyPXVMmE1neulIuDZj5lhBA6wCFDcTtulBpsUppiqy3dxc0PsDjN9tF2znPfAAcz/ctFB08rlpJEyTz5lKviAw5r3wMmYQL0c6tg0idoKwQShF2Cbznd0UqjB2YQYxtzQhX6E500aryTinGwaEBWlNXVEWWWsZVJG2tG/3XpACNBcHrXlG3OaO1+FamAGmeS885RUUTa5jzbzwROwEYcMgImylxuF8TGzbRunp4OalNN1cHhGVDr5lBG3vik1QVpgyf4w6vzugRgjJrKMLWV5w9IkjuOvR2ybr866/9EZcccHV9Tt6hER1TRpL9yE07txZNjWBLUtZrqoXnoidIGwgRIQ0IRCZwQ3M0DQVOG1Tmm4uT4GRaT2ZzwPsnJ55UHcy84/3tM1yBH/NXdV91Tmcu84odRx94ghufehj2M1eAQCcOPUMbn3oYwAQLHhVc1xNjF4nDKAxrlF0lpn8qfk7odn3V9AaSMROEDaYRBEICkTA3lhDw4iRmjSPNmlNN7eTX57QJ05ww7qsVKcsgWm0F9oZ5a5Hb5sInWM3ewV3PXpbkNgVuiigpchZhl6UPtYaIGW/aHjp5RW0BhKxE4QNRynClkoAkElnalOOwpgOrkb04veJnIMwmT9yi+Drtq8TOT/tWTUHmGfGFSHg9Txl1kRdvASHlhZ3pW65iqLpfa1awYqInSAIAICtVIHGjDEABZrMIznRSxM1EbuiJQqNmNS7lLcQY10egQVFfFY4fcf0JmJ91oFzcOLUM4Wv15GPeLqKHIBofTWb4ldlas/8NVG0UvN369ftUxCE1qSJMq4FNHUwSG36zx/TXIrOuCeo6X6KKn/cdqnnyFBk/FooDGRExJ2nLmWZKjWT9mva0f/6S2/EdrIz89p2soPrL72xcj/KCawzSe0KLcgvoejKNU9NbJtEy4tEIjtBECZMilKUF93BzNklKgHGGTSmc3lF+zuUUq0HebdbE9NXc87iebI2A7Kbl2tUjemn+WJEcx6Lmh4rayDtzF81a7OcZcnTmSJ2giDM4BYtazbpKz+lmSgFeI2k8+nM1mN7Lq0JYG6NV9W+frqyiLZzjVdccHX4UgMyJrBRXc29Y0exBmpxiLo7Meav2qTBl1jwROwEQZgjUQToqUg40UsUY9c2kGbkKzS98nRm4wlXMlJOlyEUD466boglI4whlZpDpdkSu3SjDwPURVoDAfUNpDNm7Gcao2R5BU/EThCEQhJF4IxnZIeIMLKFKq6VlStiQW47oshFAX6lZoMBtfcKUtjfFfd0rgZp3IBDtcItP6lCM5baFkjEThCEUtJEmQIL7zVFZmCbCpqX0vT6OHZqIE3VlZqhDBHVJYpar5sLIWZU12fUpZmhGMgAsG7fD7QvROwEQagkL3iu7D8fxTjxs8/Mf71KxLDu/G5tXZyYsO+orm+hAxafwgTCv7NotpZQ1iVhmQRPxE4QhFrSxPTSnMzhUVhvSTOVN9subCj6juoGEbqIYjGEZvrGr8smeJ2+PhHR6UT0ABEds/+eVrLdESJ6gYi+mHv9IiL6BhEdJ6J7iWiry/UIgtAfiVfWX9YKK0+btFks7egtqqNhhC7mXJ05XPtjNblP/0uG7qlgpw1dcwUfAXCUmS8BcNQ+L+ITAN5b8PrHAdzKzD8H4CcA3t/xegRB6JFEERI7AId8Y1/Ud/reRMguL6gTuqNPHMF191+Dt3/+Lbju/mtw9IkjjU8VS+jYLvx2wjO2fob7FT9um0ybfXSDNYP5LxnLInhdxe4aAHfbx3cDeFfRRsx8FMCL/mtkvvJdDuBzdfsLgrA8mO4kJk5IagbkPtZ1hdDH2EpkurLUDfzOLeHEqWfA4IlbQiPB6xDVcU7U/PSzq5yt+/W4bbRtsZaBMdY8EcKqBfNFXnzL0Gmlq9idw8xP28fPAKhvGjflDAAvMPPYPn8SwOvLNiaiDxLRg0T04MmTJ9tdrSAIUSCatv6qGpIn6/AGJmpkZ9OWZnlBfX/NKreEUOq+ROTReipwhULU8UPwD+dE0Imfi/xmrqfgV5Rx3I4yTaktUCGiLwE4t+Ctm/0nzMxE1NudMPMdAO4AgIMHDy4+JhYEAYkibKfKLDQv2aaN8WsXYp7Ld0wIXTDe1S1B1fT9dDSxCuqzotM3flW2Ute8Uuz6sKjG0bVix8xXlr1HRM8S0XnM/DQRnQfgRINzPw/gdUSU2ujufABPNdhfEIQlQCmF7ZSwlxVbAIVWbsYiitbRvHVQ6LxTF7eEkPRlE1dz77CD4KyACMYUOElmk4cMYJxppMnwHgRdz3gfgMP28WEAXwjdkc3Xka8AuLbN/oIgLA9Kmc4qo0TNDdZDt4/qVIXpnBXUbNurJvNNbd0SXPFLGVpPnQYaEaGvZtPdXZPo/UzP/e4Yw7Vw8+kqdrcAuIqIjgG40j4HER0kojvdRkT0dQCfBXAFET1JRIfsWx8G8MdEdBxmDu9vOl6PIAgLInFFK7aAxRe9Jmm0hSS5PJHLX2tTH7wrLrgaH7rsJpx94FwQCGcfOBcfuuym2obSZfN0rUXOsqhF6ZP0prMD8gRuEfN3tMgJw7YcPHiQH3zwwUVfhiAIOdgWLvjPmYFMa+wHfpvXmpG1tgZqVuYeYgI7zspNZGNRZE0Uq6l0jJQh1USdZYxy53ZfhlzhUsi1EdE3mflg45PnEPNWQRCiYcrzaea5UoRRmiANXoje/vxBGknTwpM6E9gh0m1FQufWt3UlVgq5bU1tPphy6c1M8+DpTGkXJghCVIzgYSbCA2zLMdtOSrvBLnJmqfRoS+qYkO9EE9siaNGdusrcEjQzOGNAEZRKBrkWETtBEKLjHM/9dKRSBG0tg1zPxwQ055bQxS7B7L8ajgn+sgZ3vqjiGqEwZXKoloepuhsGsK8ZNNZI0/6TjCJ2giD0glIE4tkITxHNzccVuiVgGlW5refWSXsDsFcKE+HKDX1GdXmhq+pI0pamC9Or6DNA3NcapKn3htEidoIg9EY+pelHd3X7TWx/kH9QRRzB6DOqy8/R9VIAEzGqA9ofywh49b5u7R0z9br+TgpUBEHoFVe04oa8oCKVFueJGRn1FdX5Qse2JL8PokZ1A1kDTUWvp999L0cVBEHwcL00la3OrBs/F+lX2stgS7MWSbELUWZORWHtxoKP18UaqMF2bAWvj5QuIGlMQRAGJFEE0gArmqvW9KEORSpdia5BuTVqVUJ39IkjuOvR23Dy1LM468A5uP7SG2sXo+eJPfU11BcP53LuBC9RcVOxInaCIAyKUoQRKYA1xhXf4JvKXYxggCN39ihbQ1eEswZyjgnOGghAsOCFNpFuwpB9NRXbxuGI3zRa0piCIAwOEWFrlFTOLS2izVU0nculLR1Vc3SdrYEiO5sDJqrzxdN9GdB6/ifGF4W8lVDMOU2J7ARBWBhbqZp4ouWHyUVkMjsXpljBKRKdujm6rtZAsYpS/HWPBABcbt80u+PsA7fGMSGeE80yNDN8h8SYH79EdoIgLAzn/p0maqZi07234AYgjShrJA2g1t0cKLcACrEG6pq+zJu/ahulEbUXHNchJ+NmLud9IWInCMJCcfNarmIz8dfYDZjKbLu2TqlppWkRoa4Jba2BiuYF63CpyLEVIF2Ugoy8Vg8IdDnvqUpV0piCICycRBHYayWmbBsxArCXhQ1+XdOejXatSFfmCV1i4IpQmlRj5ntr1rEs7uaOOtkSTAAADPZJREFUQpfzns4rYicIwlKQJsosKrbPTaRnBr5Jaq3H89cKgO25WWcJ5NN0Ld0VF1wdvtTAFsGE0KbvZgyxa3II53KudQY1SqJHlSJ2giAsDYmiuWKVRJlSdFe4UOaY0MvavJZuCUD8ZQz56wrxl2vbXHpod3kfDdNejsBR+2WK2AmCsDSYBsmYETxTqDJ9XumY0JKYbgmOvjqkhAhd1w4tsTSmkw8eVTceaIqInSAIS0WR4BW5Jbht844JvvBVDZV+EYzWDIoYhS1S6Dq32+qhMKUpZT54XRCxEwRh6cgLXhO3BFMyP7UKCiG2NPXVX7NK6GL54S1DE+k+virI0gNBEJYSv3k00J9bAhBXnPqI6tx6xKpzxhC6+E2k2yGNoAVB2DgSRaZqQaE2umtboxJraO2jKKVqHV1s94T4TaSXpy2ARHaCICw9iSIkVO9mvejBNXZQV9Rf0xFb6KJHdR0O1UcaUyI7QRBWAueWwKwLi1Uczd0S4g2tsY5FNcIea35uesLwNXvhh1yeqA4QsRMEYYUgImyPEuztZ6X2QK5IZWhitbmqa/9VVG3Z1QcvZlGKo+sh/UKjGIjYCYKwcoxSBWTadNzIDfyKejBgDaBrpFUXzQHlQtfFBy92+tKxCIumKmTOThCElcMIg0KiCGkuElrEvF2XSNKJXBuhAzr64PWQvgSGczdvgkR2giCsJK55NIiQEKDYW5dHFHdOq4Y2kSRReJ/NKougLj54faQvgeWL6gCJ7ARBWGHSRM10QnG+eH0N4mU0iez8SC5E6Oosgtr64IWevw3LJ3UidoIgrDj5NBwRYZQqjHLeeH0RJHQ09b1rIjIhywva+OB1NXutwkV1zjPPGcOW/bhtem2cDUljCoKw4pjuIphrGqyIAFXvltCVssO5vp1dUnoh6+ia+uA19cCrI9+MO1EE3eb4PH1AABQlUfuVitgJgrDyFAleviqz1C0BQJdlzLoHxwSgWduxUB+8kIrPEEpNYCM1ke7jS4mInSAIa0Fe8PLWQPltfbcEt02IW4LZY3qcPqgqSGlLV6ELcTmPFTH28WsVsRMEYW3IC16iwjzRXNeVxm4JPc0xxY5qughdk24ty1iF6RCxEwRhrXCCl2mzLMGkM/sRpT4OG9s1oa4jSxlNW5KpiOv1+igrkmpMQRDWDrcMwRQ6hG3fhthaFzt9WdVIuuoaxplu9gWhYyHO3OEkjSkIghBOmijr9VY9eC+L71q0CJTMgvEmIt7FRSF2+rKPZKiInSAIa02iCGDCblYhdi198GISq5F0m/m5slZkYSfsQex6CO1E7ARBWHuSRGGHCHsl6TlqqXYx9bFzVDdwNOcYultNW0TsBEHYCJQipEzINAr98BYZ3HWN6gaP5ix9dGLpq6JTxE4QhI0hUWZNnQLZuTx/ETpVmsLmiTlf1zaqa9JM2mec6Vbnmz15P8LUV6AoYicIwsZA1iEhs22tEk/0FjVv1yqqa5GyBObTll1MX1fNMUHEThCEjUIpAutpNJUogmJAMyHT2eB61ySqaxvJAcVC19b0ta9G0n1O/8k6O0EQNo5EzS5bdnNeW2kyiFOCo4ljQhNboKLz5AtR2pq+xm4k7dNnBxaJ7ARB2EjSRGGc6ZlITllhUV6zaM3F/TVjUJrBtI2l20ZxPmUVl21MX2M1ki49fm9HlshOEIQNpsgLz2/y7DzonDeeihz1TSI7mp4vUYRUNfO9qzp+2dKCxqavNP/7iglFckwoQ8ROEISNxfTRnB1gy1JpEyGy4pd6qUXlhJDK553ce25bZcXDuKuryXFiDfh1a+gamb4SkKp+5aLvJtKSxhQEYaPJOyUoRcgquq04lFJQbEv4PbugULTutyCjbrF4sOnrAELX1eQ2BBE7QRA2njnBI+rNKcHR5/FD19HVmr4OIHTAMNZAInaCIAiYtQZKFEEHRHddlub1JXWxLIL6LkbxUTQtCAIwZ6YbAxE7QRAEi7EGMgvNk4YdVZrQl+mrEYzux27rgReKEzaGqYDdj9HRpQYRO0EQhBwmolHQ4+pF5mQjkqb0lcGMEdX1JXRamyUcM78vAoiGqZMUsRMEQSggUYStRGEvtxbPp7UPXtuLqqCz0LVsQVYFM0NzeSQ7xFydQ8ROEAShhCRRGMEUrhRbA6GVcsVOY3ZNX8aO5nTJ72uGASowfUTsBEEQKpg4JbC1xfHeWwYfPKBDVBc5mgsSOcvQPngidoIgCBVMnBIApAlBa54pXFm0yXlb1wS3uD3WNTRZStFXI+nKc3bZmYhOJ6IHiOiY/fe0ku2OENELRPTF3OufJqIfENFD9ueyLtcjCILQB8prHK0UYZSoiVC0sdmJSaP1erb3Z6pUFKFjZoxL3N8rr2EB7uZdy2A+AuAoM18C4Kh9XsQnALy35L0/ZebL7M9DHa9HEAShF/JrzhIreumALgl5gtOXkUXOnbtN+nTo9KWjq9hdA+Bu+/huAO8q2oiZjwJ4seO5BEEQFoZJZ84P1GmqkNpG0SHDeMzAri5KdAvDY4ocM+P+H/4T3nPknXj759+C6+6/BkefOBK0b4zm1m3pOmd3DjM/bR8/A6CkXXYlf0lEfw4bGTLzbtFGRPRBAB+0T3eJ6JEW51oHzgTw3KIvYoFs8v1v8r0Dm33/S3PvyWuT00enj94AMsHSMRzDv/K/6P0f7/939mL24x5O+aYYB6kVOyL6EoBzC9662X/CzExETb+z3AQjklsA7gDwYQAfLdqQme+w24CIHmTmgw3PtRZs8r0Dm33/m3zvwGbf/6bfe4zj1IodM19ZcRHPEtF5zPw0EZ0H4ESTk3tR4S4RfQrAnzTZXxAEQRBC6Dpndx+Aw/bxYQBfaLKzFUiQSeK+C8CmpiYFQRCEHukqdrcAuIqIjgG40j4HER0kojvdRkT0dQCfBXAFET1JRIfsW39HRA8DeBgmJ/0Xgee9o+N1rzKbfO/AZt//Jt87sNn3L/feEeqr+7YgCIIgLAvDtJsWBEEQhAUiYicIgiCsPUsrdpvciizCvV9ERN8gouNEdC8RbQ1z5XFocP+H7TbHiOiw9/pXiegx77M/e7irbwcRXW2v+TgRzXUiIqJt+1ket5/thd57N9nXH/Pmw1eGtvdORBcS0Snvc7596GvvSsC9v42IvkVEYyK6Nvde4d//KtHx/jPvs7+v9mTOGmLZfgD8Fcwic8C0Ift4yXZXAHgHgC/mXv80gGsXfR8Luvd/APBu+/h2ADcu+p5i3z+A0wE8bv89zT4+zb73VQAHF30fDe43AfB9ABfDrDn9NoBLc9v8EYDb7eN3A7jXPr7Ubr8N4CJ7nGTR9zTQvV8I4JFF30PP934hgF8EcI8/nlX9/a/KT5f7t++91OR8SxvZYbNbkbW+d7uM43IAn6vbf4kJuf9DAB5g5h8z808APADg6oGuLza/CuA4Mz/OzHsAPgPzO/Dxfyefg6lsJvv6Z5h5l5l/AOC4Pd6q0OXeV53ae2fmHzLzdwDo3L7r8Pff5f4bs8xiF6sV2XeI6FYi2o54bX3T5d7PAPACM4/t8ycBvD7mxQ1AyP2/HsAT3vP8fX7Kpjf+bAUGxrp7mdnGfrb/A/NZh+y7zHS5dwC4iIj+k4j+mYje2vfFRqbLZ7fqnzvQ/R52iOhBIvo3Iqr9Qr9QPztaklZki6Dne196er7/65j5KSJ6LYB/hHHcuKfdlQpLzNMAfpaZnyeiXwbweSL6BWb+30VfmDAIb7D/n18M4MtE9DAzf79s44WKHW9wK7Ie7/15AK8jotR+Cz4fwFMdLzc6Ee7/KQC/6T0/H2auDsz8lP33RSL6e5h0yTKL3VMALvCeF31mbpsniSgF8NMwn3XIvstM63tnM3GzCwDM/E0i+j6AnwcQpZfiAHT57Er//leITn+73v/njxPRVwH8EswcYCHLnMbc5FZkre/dDgBfAeAqlxr/7paAkPu/H8Dbieg0W635dgD3E1FKRGcCABGNAPwulv+z/w8Al9gq2i2YIox8dZn/O7kWwJftZ30fgHfbisWLAFwC4N8Huu4YtL53IjqLiBIAsN/uL4Ep1FgVQu69jMK//56usy9a37+97237+EwAvw7g0cqdFl2RU1GpcwaM7c8xAF8CcLp9/SCAO73tvg7gJIBTMDnfQ/b1L8O0IXsEwN8CeM2i72nAe78YZsA7DtOmbXvR99TT/V9v7/E4gD+wr70awDcBfAfAdwH8NVagOhHA7wD4L5hvpjfb1z4K4J328Y79LI/bz/Zib9+b7X6PAfjtRd/LUPcO4PfsZ/wQgG8BeMei76WHe/8V+//2yzCR/He9fef+/lftp+39A/g1O75/2/77/rpzSbswQRAEYe1Z5jSmIAiCIERBxE4QBEFYe0TsBEEQhLVHxE4QBEFYe0TsBEEQhLVHxE4QBEFYe0TsBEEQhLXn/wEiwuUYyUMLWAAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbsAAAGfCAYAAADGVHw+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXuwLEl93/n9ZVZVd5/HPfcxlwGGpyRkQArbWFcDCAskQAJFKAQ2YCHJMrIVSyhiCf/h2A2xIa+0gcIRKHZjV96QvBayF4FWBJKRBLMbrAks7A3FhpC4Y7QyAwKG4XFnhrnve8+ru6sq87d/ZGZVVnV1n77n9H3V+X0maurR1d3VfWbq079f/jKTmBmCIAiC0GfU3b4AQRAEQbjdiOwEQRCE3iOyEwRBEHqPyE4QBEHoPSI7QRAEofeI7ARBEITesxLZEdGbiejLRPQ4Eb234/HXEtF/JqKSiN7eeswQ0V/65ZFVXI8gCIIgxNBR+9kRkQbwFQA/AuBJAJ8D8FPM/MXonBcBOAHgvwHwCDN/LHpsl5k3jnQRgiAIgrCAZAWv8TCAx5n5CQAgoo8CeAuASnbM/A3/mF3B+wmCIAjCLbEK2T0E4EK0/ySAV97C84dEdB5ACeD9zPzxrpOI6N0A3g0A6+vr3/fSl770kJcrCIIg3C88+uijV5j57FFfZxWyOyovZOaniOg7AHyGiP4LM3+tfRIzfwDABwDg3LlzfP78+Tt9nYIgCMIdhoi+uYrXWUWBylMAnh/tP88fWwpmfsqvnwDwnwC8YgXXJAiCIAgVq5Dd5wC8hIheTEQZgHcCWKqqkohOEdHAbz8A4DWI2voEQRAEYRUcWXbMXAJ4D4BPAfgSgD9g5seI6H1E9BMAQETfT0RPAngHgN8iosf8018G4DwR/X8A/iNcm53IThAEQVgpR+56cDeQNjtBEITjARE9ysznjvo6MoKKIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9ZyWyI6I3E9GXiehxInpvx+OvJaL/TEQlEb299di7iOirfnnXKq5HEARBEGKOLDsi0gB+E8CPAXg5gJ8iope3TvsWgJ8D8JHWc08D+BUArwTwMIBfIaJTR70mQRAEQYhZRWT3MIDHmfkJZs4BfBTAW+ITmPkbzPxXAGzruW8C8GlmvsbM1wF8GsCbV3BNgiAIglCxCtk9BOBCtP+kP7bS5xLRu4noPBGdv3z58qEuVBAEQTie3DcFKsz8AWY+x8znzp49e7cvRxAEQbiPWIXsngLw/Gj/ef7Y7X6uIAiCICzFKmT3OQAvIaIXE1EG4J0AHlnyuZ8C8KNEdMoXpvyoPyYIgiAIK+PIsmPmEsB74CT1JQB/wMyPEdH7iOgnAICIvp+IngTwDgC/RUSP+edeA/CrcML8HID3+WOCIAiCsDKIme/2Ndwy586d4/Pnz9/tyxAEQRBuM0T0KDOfO+rr3DcFKoIgCIJwWER2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9B6RnSAIgtB7RHaCIAhC7xHZCYIgCL1HZCcIgiD0HpGdIAiC0HtEdoIgCELvEdkJgiAIvUdkJwiCIPQekZ0gCILQe0R2giAIQu8R2QmCIAi9R2QnCIIg9J7kbl+AIAj3L8aYxj4RVYsg3EuI7ARBmIGZcfHiRVze38c1w9i5BjAD+yWBQDAaUARoRUgSQkKEVBMGicIwURgmCYZZgmGmsZZqDNMEg1QhTTWUkoSScOdZieyI6M0A/iUADeDfMPP7W48PAHwYwPcBuArgJ5n5G0T0IgBfAvBlf+pnmfkXVnFNgiAsZm9vD1du3MCF3T1885kJru/luHS9wDZbjKcWTAABUIqgQVCKkJCCIkKilNtXCokipKnCINXIEoVBopAlCqPM76cKa6nGKFFYH6RYH6XYGGZYz1IMMpGfcGc4suyISAP4TQA/AuBJAJ8jokeY+YvRaT8P4DozfxcRvRPArwH4Sf/Y15j5bx/1OgRB6ObmzZt4+uY2vn59H5cuT/H0jRyX9vdxc3uC7bxAXgAlATRnUYqglIImglbKC657SZVCqjUSL7xhpjHw4humGsM0xzBTWM80NlKNjVHmlkGGYaaRJPpuf11CT1lFZPcwgMeZ+QkAIKKPAngLgFh2bwHwP/jtjwH4DZKkviCsnMuXL+Nbu2M88a19PLMzxYWbu7i5O8XudIJxwdjLgaIEyhIwZp7cmttKMYiM3w6POwEmSiHVQXQaqXbCS5RG5o9niZNfLTyNYaYxyhRGezk2Mo2NLMHmKMPm2gCbwwyDTEu7n7BSViG7hwBciPafBPDKeecwc0lENwGc8Y+9mIg+D2AbwD9n5j/tehMiejeAdwPAC17wghVctiDc3+zu7uLCjW185cI+ntmd4MLVHVzfzbEznWBclNjLgbwAigLIc6C0QOmfS9Fa+TWhKbt43d5WiqGUQaEMcg0X+XnhJUoh07oSYVpoZFpjPHXiG3jZDVMnwJ1MY5SVWB8X2NidYGOYYXMtw+Yww/owhdaS5hSOzt0uUPk2gBcw81Ui+j4AHyei72Hm7faJzPwBAB8AgHPnzvEdvk5BuOvs7e3hm9du4KtPTvDEjV1curmHnd0c23mO/WmBcQFMvdyKwkVuBZzgSgAWkdQ6FgWAGFAm2p8ru/ZiobVFoUqoKsKbXWelxqRQ2J9oZJnGKFr2vPjWpwab4xyboxQbowFOrg2wNkykbU84EquQ3VMAnh/tP88f6zrnSSJKAGwBuMrMDGAKAMz8KBF9DcB3Azi/gusShPuebz5zGV/59hhPXLmJZ27u49pujr3JFHtFiUlhMI7kVhQueguCC+uDJNeQXXsd5Gdc9eV82Tm5KgVobWCMQa4Mkg7pZdpFernRyHONSaYxjqQ3yTT2M429wmBjUmBvUmBzlGFrLcPaMJX0pnAoViG7zwF4CRG9GE5q7wTw061zHgHwLgB/BuDtAD7DzExEZwFcY2ZDRN8B4CUAnljBNQnCfcl0OsWFa9v48lO7+Nrlm7iyO8WNnRx7RY69vMC4MI3orSyBwkdwBQATrWPaqZB5opsrPQDKi0+ZxdKzto72rLUolUFhNTJtUViF0lqU1iJVFoV1j+WFxjRaJn7ZyzTGhcV4WmB/WmBrfYCttQGyVApZhFvjyLLzbXDvAfApuK4H/zszP0ZE7wNwnpkfAfBvAfwuET0O4BqcEAHgtQDeR0QF3A/QX2Dma0e9JkG4nyjLEheuXMeXntrDE5dv4upejhu7OfbzHHtFgXFeYD+P5BbWqFOUYXvZ/D63zm2026Epuva2AkAWULYZ7WndlJ3WbltrhrUljLXQqpZdogwyq6v90loUpUXeWkrDmBqLSWlQlAbTwuDEKMPGKJXUprA0K2mzY+ZPAvhk69gvR9sTAO/oeN4fAvjDVVyDINxvXLx6A499e9cJ7uYUV3en2Jvm2C9KjIsc48JiGgmuKIDSNNOUXVHcYbB+HWRn4aQW1nPFF6I9G8TWjPBq4QHWWlhtYSyjtBaZ1rDMMJZhrIXxwjPGovBLadzjhbGYlhaFYUwLg0leYlpKlCcsz90uUBGEY8V0OsXXL+3ir5+8iie3p7i2O8WNvRzjosR+nmNcFNgvLYq8FcWZZgS3Ksm1iSO+tvC4tR3vB+kxN1OaYb+5NmC2lehKbWFYwzDDMLvjzDCGURonuhDhGeuivNwYlMYiLw1OrQ8xGiTSlicsRGQnCHeAvb09fOHpXXz14nVc285xZS/H7n6OcVliL88xKQuMc4uidN0FSt8XrihnU5Xl4rdaGQwnVMKs6CrJoRYfA+CQ3vRpzG7ZueHImEsYxTCsYL3kZhYbIj8GM+pIkBmWgdJYMAObJsOJtUyEJ8xFZCcIt5EbN3fxhYs7eOLiTdzYyXF1L8fuuHARXFFgv8gxKW0VwTVEx4dvk1slbelxa7tzsU5sQXjxdvOY8eJDJbOuxfCs7JiB0qdAS+PSoFtrA+mXJ3QishOE28Du7i6+cHEXX3vmJq5t57g+LrC959KU47LApCwxzgtMvdjKUGFZNrsPBNHZxW93RwjSY3REdF0LA1wCSgO6S3gI2xYAg8GwrCvxtRe0ZWfhtwHmCUprYRk4tS7CE2YR2QnCChmPx3js4ja+/sw2Lu/kuLGfY3vPVVSOywLjonQpy8I6wZVRJFfUgrvb0dwinJqWFB5QGZL9B2msddhmMBtYzZj7Tx6iO38hXH83HLb8amstQyrjbAoRIjtBWAHGGDzxzDb++ulruLib4+ZegRvjAvuTwldW+miuKJFHgqsiOjsruttRgLIqQpSnov1O0YXH7RzZoY70wqtO50V24bX8TnysgW+2O7k+QCIRnuAR2QnCEblyYxePXbiGb92cYHs3x41xgR3fLjcuSpe2LApMSlsLLqq0LFttcwXuvWhuHnGUN49KfN6Q86RXnc8G+QHvGwpR2vUoBGo8LilNISCyE4RDkuc5vvjMTXzjmW1c3S9wY2eKm5MC+5OQqiwwLktMyxLTkjvTlnHbXFjuF9EFQnVmvN8V3QFebEs1QHYLb16xJRFBkVtrAkiRn5KIsLWWSedzQWQnCIfh6o1dPPata3hqd4qbey6a294vMMmd4PaLAoUxbh1Fc1VkZ2YlV9zdj3QkgsPmRXncdfIBELWER26MTqr+MSA/WzqRqWSnCFDkZKeUE97mSLolHHdEdoJwC1hr8dWL1/G1p7dxaS/Hzr5rm9sbu1TluCxdNFcUGJdmRnKl8QtWPwrK3aYrklt4cofwZn1kURJBEUEZgvKSC+eqgjD2ooulp7z0tCJoP6/e2jA99GcT7n9EdoKwJJPJBI89vY1vXNl1VZb7BW76IpRJUWBijCtEKQpMTVRtGQvPzqYu+yC6mKW7SQThecHFoqsnk2WUFH1DYVJZLz+CGz8tjuoq6SkgUQWUVlVKc5DJLe+4In95QViCazf38MUnr+Npn7a8uV9ge1xgkpczEd2ktHUUFwuvo5P4vdB/7naw6HPFwRt54VH0QBBdDcOECK8rujOEvCAoZSLREXROSJVFMincbOp+kYKV44nIThAO4Omru/jK09fw7d0c23s5tielq7acuu4EE99G5wpR6tFQGsLj2SG/+iq6QNfno9Zx8icG2YXZ0qvHqwjPyY78oqLUJpGphJd44SWFQaIIY79OJgUSrTBINTZkTrxjichOEBbwtYvX8fVntnFpN8fO2EVzO5OyEt24LKvUZW64Fl3pZygoXNeC0C5ncDxEF2hXZrY/d1AO2Wi/FeG1hVcVn0TCK8hAlYSJImhtnfCC+JRBqghpqpFNCmRaSTrzGCJ/cUHowBiDr168gW9d2sWV/dxJblw60eXtaK5EblzXAmPQ6EtnuG6X63vqsotFhZdxlEf+ZLJObBZdsuOG8FRLeIosVEmY+mhOR+JLtUXm05lpoiWdeQwR2QlCC2MMvnLpJi5c3sXVsau43J6U2J242bInVd+5ElNjMI3a6IoiSl+2Upf3U2fxVTIvorOtbWI39x35SfWC+NrRnSHqFp410K10ZhrSmaVFpg3SvESWlhilGmsiu2OFyE4QIowx+Oqlm3jy4g6ujQts75XYHhfYm5R+wtASuTFOeMZgUhiYOGXZMWPB/dpZfJXMS2cujO4WCK8rnamtQumjuzxRSEqL1K+z0mKiCYO8xCQtsZ9qDFIt0d0xQmQnCB5mxhOXb+LJyzu4Mi6ws19gZ1Jgb+qkFiK5MM7l1IsuHhWlLJrDfx2XYpRlaBemdEmP4AZ6XiadGUd4Wik31Y/vSF6WFrkiZH6dlxZTTZhoiywvMUk1xqnGxii7s1+CcNcQ2QmC5+uXb+DJy7u45iO5valbJnnoUlCnL/PS1B3EvehM1EYnouumK7KLu9tV6cxIeNzan4nuYtGFCK+0KLSTXOLXqZfdoCgxyTX2U41hlshg0ccE+SsLAoAnr2zjyUu7uDYpsTM22J+W2J2WmOYGuTGubc7U26WJilG86Mpo5oJQeSmia9KVzrTRdrW2teSsrfebi5/Y1VqUfjHxvmHkpUXh17lhTErrlsJgWhhM8js177twtxHZCceeyzf38c1L27g2DkUoJXanBpOJi+ImVSGK225LzkZDgAXBWdQTnQpNYsHF0muvY8kF6TWXWnJt0RlmmNKitIwiWhvLyI1FXrjofFpaGCM/SY4DIjvhWLO7O8E3L23j6qTEzrTA3thgPy8xnpbIrfWSM8h9YUpZsktXhmguSA/1TToIT0Q3n1hsXdEeY35kZ/0SR3chwjOt7bIlvBDh5ZaRFwbTokQhsjsWiOyEY0tRFPj6lR1c3Ztif1pinBvs5SXGU4O8qAU3NS6yK4yt5eajOWvrm3NIXS49i80xZkZsmBPt2abgKvFF0Z1hbkiuZFsLryOyK4x1BSuF+ztPClNNCCv0F5GdcGy5cH0Xl7fH2PaFKPu+6jIvasFNyxKlMSiMcZKzfjHRDRe15IL0RHYHs0h47eiOeV67nRdeFNVZy7XwzGxkVxqfyixd1D4pDIpS/mJ9R2QnHEsu3tjHxWtj7OQG48JgPDVuZJSp8VGcK0YprMWkdKKzpl6YnfTi9qc4shMOpiuF2bksEF1XoYppC89HdKXxa8tuaLcgwdJIKvMYILITjh3j8RRPXd3B9UlRpS3HhcEktyi86ArrfvnnZQlj2YktLBxFdpiN7OS2uTxLy64dzc1pt2sslmvxWUZpeEZ6pbE+rSmpzL4jshOOHd+4vo+b+zn2coP9qcHYt9cVpe9aULr0ZWHdL/4qovOiC5Fdu6ow7AvLc5DoLBZJbgnhVVEfV4IL24Vl98OmNMi9AIX+IrITjhVXt8e4fnMf21NXcTktDMaldQUpVVTnJFcY4woiOqoCLc9KTiowD8dho7uuVObcJZJcFeFF0V1prLTb9RyRnXBsKIoCT1/bxW7uBDctLSaFxTSv5VZYt86tgbFci60tvTmRncjucCwru3Ag3l8qsvORm2lIz3U4L7zoCmMlldljRHbCseHpG/u4uT/FXmEwyQ2mucGkcDe6MhJd4UvWu0THHaKTfnVHI5bawn2eL72DhOcKipqyswzXbcHYqrhFXNdfRHbCsWAyyXFtZ4rt3FVdTn0VZl6YRkRXi46bacvQt6sjhdm+OQu3TjuSQ9f2IVKZjWPWLQwnPlfAUguvNK56U+gnIjvhWPD09j52JrkfE9GNjZiXPqqLRGcr0dWRnW1HdujuAC0cnq6o7pbb7sD+HC86oFozNyUXF62YOOqTLgi9RWQn9J79/Slu7Ez9VD2unc6Nj2h9MUod0Rmu01lx2pLR2sZse51wNNqR3MwxniO96AVi0TEzLOrojiPJuQi9mdYMHdCFfiKyE3rPpb0xdicFxqWphoiqojpjG6nLUH4ey60R0cU3XIjoVslBbXazO/WxOJXZiO64ju5sJDk7I7wwkLTIrq+I7IReMx7nuL7tilLy0o2JWFiuoroyVF3Ctd1UqS9vsFAMEYQHSER3O2m7rL29KI1ZndeK7uJ1nZ6OUtXhdW3dzif0D5Gd0Guu7k+xn5d1+rJwAwAXpXEjafj7mrHW3yRtsyjFiy4IL8gtrr4U4a2GGbF1bLcfbzwniuzc/qzwAFQRX/U6zLDWVoUsQj8R2Qm9pSxLXN+dYD/30/QUpmqrCyPlM7t2OgC+fcc9t3EztLXw3Hn1YyK61TJPcu22uth2XZFdLLx4HUQXIviQzqxS1OyiP6F/iOyE3rI9LrE3LtwIKaVFXjImhYUx7ItQXAGD62/V6nLQaqvrEl1XdCGshnmpzIXPmSM8i7odFqHNtRXdhfSlsfLzpa8kd/sCBOF2cWV3jHERJl9lP/6hRWmN+0UPbnQ1AJoRXbgxhnWX6I4Do+f/NvT615Y+v7j+KkwvvvVQ73Vgm11YQkTXfn5HGrMhONSiiyXXKG451JUL9zoiO6GXTCaFi+qMi+iKMP5hYRGyVKGDcd0ny9Y30HCDtM3ilLbo+hIHDB78ONJTn537ONHyrxVe57DCC8yN7loPHJTGjNvr4oiuPr/jPYTeIbITesnOtMR4WlRjHpalxbSoixBCast6izHqX/ux4OIUZl9ENy9SuxWhLYIISE/9+aFkx611e3vhc+cI71bet/06Qn8Q2Qm95PreBFM/2HNpbDWFi+tmwFUa07XNhYiOZ1KXlei4Kbqum/K9yLyIbVVim89qvpmuVzlsGnkmJTrHarf/uxHuBiI7oXfkeY79aYGJn7MsL93o9u12uTqN6SOBuFEuWrfTXPey6JITn8fgOf8ORHXceXdu3od/03nf7zKCW/ZvQgAo+mLEb/1HZCf0jt2pm9EgLwxKU89fZn1Ux1EaMyS7DDdF2NUFoWt9L9AW3N2OTJiB4vorb+MboJFajln00ec9Fr4vihahf4jshN6xOy0wLko3E3XphGdMLbO6fa4WIBBV+fnXqYscWvut9d0iTlHeKcEt0551lGrMue875/itfmwiasqNCEQERYAigrrbvxSE24bITugdu+MSuWEUZRjcOQwfFaZ7cee1xdfVXheLjTC7f6dJTnweg2f/EUgVAG6f5OZJzU6fhf2v/7Pb86bLskz4Rf4UaqYsa8G5747gOhsTEUgRlF+E/iGyE3pFURSYTHPk1Uj2rtsBUMsN8EUpiG7q3BLaAcI7TIHEUVl78f8MNbi0UsF1So0VJt9+B8rtV6zujY5I10duH4u/F+UfVd52lfCC4KKIjqKoTiK7/iKyE3rFeFIiD7MZ+Pa60GEYQKu6stluNzey6xBe9FK3Hdcm9/tVNHIU2nIze9+J8YX/6mgveoeg1tJ4jOJtaq5bwmunMYP44m2hf4jshF4xrkZKcZFdERWexHIjIi+xOo1JPmxbJLp4fScIfeIOe/9tyI0TTL79tnsqYlsFbdFR1xrNqE6pOpWpokVLCrO3iOyEXrFXFihsVH3JqEayn2mjg7sJNtrrgCp8i9t8/OYdTV8eNm1Zp2bvb7kd+LGpLbrZ58bCU6glFxekaCJoraC1QiKy6y0iO6FXTCalHzWl2a8uHk+DWsc4aoxjoDvCQ53CjAtVbhej5//2LYmuL4JrMy9lSdQ83miv8w9W6Ul/jFSz6lKreFFu0TI2fl8R2Qm9wRiDojBVZGe47lZQ24ki8bltVaU0WyKbI7zbXYk5ePDjS6cug+Tup7a3g5hXjLJIfPW2b3tDS3QUIrpacCqSXeIXKVDpLyI7oTcY4yZkteyHBeNm5/BYcsBs1NeuxKyiBwKI62feztthcuLzSE999kDRuc9FmDz9D3oTxcUs/PjUJbmwHQuuW3She4Em+MWlMVOt7nqHfOH2IbITekOeu0GfLUcTcHLdbtMlPOWPN0osaTbCQ+uU+hVWy+DB/3Mp0fUpkouh1vbMErehdgkvElwV3VHdXhckpwhQIaLzonOyE9v1lZUkqInozUT0ZSJ6nIje2/H4gIh+3z/+50T0ouix/84f/zIRvWkV1yMcTwwzytCB3Id04QaIVgUmUEd1IY1ZFTy0ClTC3TWW3O2K8EjvL3ycGZg8/ZO9FF1g0fdatdfNiera6UvljzXa6HwhSqIIqSKkWjnhJfp2fzThLnJk2RGRBvCbAH4MwMsB/BQRvbx12s8DuM7M3wXgfwHwa/65LwfwTgDfA+DNAP6Vfz1BuGVKRj2zAUUzVEchWiy8ONqrjzcLVGKpNSr8bsP1Dx78+OITmPH+S1fwlcn/eBve/d6iM6qLH18Q1c0sUftckJz2gku0QpJoZKlGoiWq6zOriOweBvA4Mz/BzDmAjwJ4S+uctwD4kN/+GIA3kMsXvAXAR5l5ysxfB/C4fz1BuGWMm74gaoFrVVyiFpuiOo0ZClVCiixur2t0U8Os/FZ5e0xP/fn8FCYzfnJ7Bz++vw8i4Insp1f4zvcOc+UWlg7J1f3lVFNwoVsBKSjdrL6Mo7o0Ucgkhdl7ViG7hwBciPaf9Mc6z2HmEsBNAGeWfK4gLEWI6mxopwu3SF+lF47A97Uj8pLzj8xILk5lYs7Nd6WfYHEL4D+/dsO9d3TD7xvtj9X4cUEHLbXg2tGda5tzSxBekF2WaGSSwuw9902nEiJ6NxGdJ6Lzly9fvtuXI9yDBEGxjWI7DsJif9MMw0Gxf4waMmtILmq7W7Tc/g/morrjwtzvegnR1R3FVR3dxWlLRUg1IQuLF12a3De3QuGQrOIv/BSA50f7z/PHOs8hogTAFoCrSz4XAMDMH2Dmc8x87uzZsyu4bKFvMDNMVVlZF6iw7xwXorgQ1bWoFHA5AAAgAElEQVQrMIFIXiF6OkB4d4oQ1fWZA39UHCA73ZG+1IqgtBNcqLxMtMIgURikCYZZglGqJYV5DFiF7D4H4CVE9GIiyuAKTh5pnfMIgHf57bcD+Ay7crlHALzTV2u+GMBLAPzFCq5JOIYQCAohOiNA1b/4yVfkhSKGuksCqsejrGfz5or5orsbt8i4/2CfOEh2Dbmp5lop116XkBsJxXUxcMIL3QqyJFq0Qpa6wpRBKinM48CRZefb4N4D4FMAvgTgD5j5MSJ6HxH9hD/t3wI4Q0SPA/hnAN7rn/sYgD8A8EUA/x7Af83M5qjXJBxPGK7NDvA3RwYoDPBMmClSqVKaHVFdtQrbC6R3JxJg/9faWiU5ZuA78o/cgXe9sxwkOhXJTsXrILZIcgkpJEpBJwqpj+zCMtCEgVYYJBrDNJEhwo4JK+lUzsyfBPDJ1rFfjrYnAN4x57n/AsC/WMV1CMebMP5hLK2qDa4V4YWILo7w2lFdtR9usjwruTB8mL2dH4wIv3jmBXjPtV8++Nz7mLjadeYHhYpEp6K1clWYWqlKcm6cS191qevUZZzCHGYuhTmUqO7YID9phN6gKZp1OhYbgGq8L6AuWvEbVCuumbps3VwXtifdZg7qbN4H4nFHG6KjOaKjeroeTeQiudC9wEuvnb5cSxVGiXZtdVmCTGR3bBDZCb0h1QQmP/s0UNWp1CNr1MUMIYUZ2vHidFlnVBdutrh9smObLnw8OfH5Fb3TvUtbeO7HyxzRKZfCrFKWygkvVaGzeC25YeqXKKpbHy7+voV+IbITeoMigoaXm/JTuoS0ZtX2Fg0QjGZas57oc1Zycd+2Ltmt4rY5febvzy08IXLjZh4HwldQ/bDoSl36RSsnOrfoajvVkeBSL7xEYeijuvVhikTa6o4V8tcWesMwC32r3H4dFdRiq6MwrqMyrudA60pjNiIL1BFHvKyCg2YvOA6pzEAlPNWxUNiORddc0qRVfellNxqmWBumWMtkDPzjhshO6A1KKSTeRspbSPm0JiLphVHxg9jiKI9i4bWjugWpzFX9j8S8+JX6nsqM08Kx8GhGetSK6nz60qcwG6JLFdZTjY1BivVBio1BKhWYxxD5iwu9IU1TZImuZKTgbpJNebXTl23pdRenhKiuK7pbaSrz2+841qnMdjtoGMKtjuaaUV2q6n50bl9jmGmMomUt1RhlCuvDFJujDKOBRHXHEZGd0CsGw7RuqyOqU42KovQkzbTZhXY7tCQXttESX5AcgJWmM8vtV6C4/qr5wutxKpOidRUtUyS80E6n/QDOSjcXrZBlrpP4INUYpk52G5nGxiDDaJhiY5jKaCnHFJGd0CtGmUbi2+3Ih2hV4UkQmp/tQHVIL5wzs2B2P47swna2gs8wvfjWhY8fOBXQfUqjAjNsRwasi1JqucVLluhGUcogq9OXG6MMJ4aZpC+PMfKXF3rFepoiqaRFVVQXV1wqiqb0IcxIr53CbEd1VbSHbuGtIknGZtR5nAhIT322l2137TbQOJqjKsKjWm5aI9NOfJmeTV+uZxqbgwQbo0zSl4LITugX65kbHip0LlYqEp5qzWDdSl/W+7PCi2+8VcEEmssqi1WmF3/iWLXdtb/Hrv51RO5710SN9OUg0RjEohtobAwSbHnRnVhzi6QvjzciO6FXrA1TDLWCBmrZ+bRmGG2jbrPDjPS6UpmqQ3qNopWO5ajFKuX2K+ZGd4Bru+tbOrPxXRKgdaudVAFECuTnoRskTnSZb5tbG2isDxJsDBJsDjQ2RxlObgyxtTYQ0QkiO6FfpGmK0SCFVqpq8wnDSTX727k4LC5QqeZCU8qNou8juHZkFwRYrTGbylxFOvOg6C499dneCE8BbkAAhIKiuCCl7m6QKAVNbogvTa6dbm2gsTFMsD5I3HamsTXKsLU+xImRtNMJDvmvQOgdW2sDN3+ZDgJDHd11pDLjf8KAwnGEF0supNIaZfDoXrRfDsuBlZk9ab+LfySE70zraFFhcRWYiSJoraE0YZC6aG6UufXWMMFJH9Gd3hjK2JdChchO6B2n11Mkum6zC5N66ngIMeqO6HQswaqKs9Vm19WGh1nRraJLwvTiWxenM3vQfjfT9ukj5kpyCdxUPYlGqjUynWDgC1N0QhimLrILojtzYoQzmyMRndBAZCf0js1RhpF2o94nuhadE1qzja7RLkctCVZtfs3ClLboqiGs0F1oEcR3WBalM4H7u/2uEc3BpX6TEM3pOoWZaIWB1hj6drpU+ba6QYIkVdgYaJxaE9EJ8xHZCb0jTVNsbA6rqsxQkalD/7uq3a3VTtcR5XX1u2uLrkprYjaqW0WEt2w6834UXiy70FYXpzCTxIlumKTIdFJFdcMswcYowSjVWBskGGQaJ9YzEZ0wF5Gd0EvODIcYqOYcZ1V/O9WsvJyXzgxteKFYpTkIcbNYpSqkQB3VhRt4fDM/bE3g9OJbe9d+FyK5OKrT2kV2SRJERxhqjYF289CNkgSjLMW6L0g5MUpwcpRgazTAQPrRCQsQ2Qm95OxmgmESZqoObXeYEV97SaIIb6ZYpZ26VLPFKkF4sejaEd5RhHdQ+93w2X98yFe/s4Rq1QSum0YKJ7c0AZLULWkCjNIEa2mK9TRxoks01gcaJ0YpTq6n2FpLcXKUYXOYwlpgZ1ogL8xd/WzCvYnITuglw+EQm+sjZIrc8GGq7mi+SHZxe52OpafUwlRmY+kQXoj0dLR/GA5qv4PK7/l0JqEWXbUkQJrWS5bCSy7DWpZiLUmxkWXYXEuxtZ7h5FqCzUGC06MU64MEDGBaWkwLg+1Jjkle3s2PKNyDSNwv9Jazm0Nc3dnHuFRItEUZyy6utoz62wFAohTYchXpwVpAKVhmQFkA7lBIYbZ/MjL7Y+5UEOrZtxGtDeppbJal3H4FitE3kZ76LLr6SYd0JnDwGJt3gxDRpfGSAFnm1mnmZLeWpljPMqylKUZJgrUsw4n1FFujFCfWM5z224M0gbEMLg0SpWAImJSMMKXTQNrvBI/ITugtD2xk2Bhk2M8nmPpoLtEEXTrhGUXQrGAisZXWGSpEdGwttBeeDQNieosxfJQWRXwMQEUGIy+8OJqLHXUY4U0vvhXJ1qMgXXQ+fq8KTyFKWfp1FkVyYb0xyLARZJdlWM8SnFhLcXIjw9Yoxck1N1UPEcEyoMn9wLDMmJYWmVbIyWKbC5xgYJCJ8ASRndBjsizD1maGnUmOqVEojIWxrVSmrVOXAJrRHdtKgKwImqMRoGHBmBVVkFpDgL4JqUt27pWqIHBpps/8fQyf+/ud0R1QC8+MX3jgDOh3glCAEmSXoY7istRtZymwkWbYGAywkbq05Ym1FCfWXBR3ZiPDxihFomflxcxgF9DBWAtbAmBgBwxGiqHMTH7skf8ChF7z7PURrlwfQyuDVCsUil10Z+NxM2vZVVEc15FeSGvqKKUJAMwW7YnFqxRmhGYX4XXJjuBEZ/yyLAelMwHf4fzsp+667GbSluTTlkF2GTDMNNbSFJtphs3MTcdzZiPDqY0MW2spNn3KMoxxGdLCleSiKM8wkDCjtBbWpzSZIbMeHHPkry/0mq3NEU5sDLGXlygMQ2sLXXVJYNiW7ACAqdl+F4QHtiAvPFLK32jr2K4rrRnMRuSFx0AonYhFF7BYPq0ZUpQLhZfeWPLVbg9VFBfWvl0uy+qIbj3LsO7TlifSDCfXMzywkeH0RoaTGxmyZPY2xVEjaCy+YDZm/7uECNPCgtl96yK844v85YXe85ytEXZ2JtgvLFKlUGgLbVz7nQkRHtcRQxKKUTxxZEf+uGELy+G8uhIlTkcGATXSmj5nGSQX+bAhv2XTmgcJj4uTS77SatFoiU7VUVyI6tws4hnWMxfNba1nOLOe4ezmACdHGbI0gVbRjxB0V7E2xMfcjPIsQyugNAa7OYPZYm24iil2hfsNkZ3Qe05vDLG+PsBaYZAbi8wqGMMwllyFpiZYW99GDdAoVgHcCPukXGoMSgHsbrKmkqI3WDRgtLGuWKUz0rOzkV1ouwuFnL7Z6UDmCY+tBtQUGy99L7g4ienlN92RlGYGYAAvOaqjuGzguxQMNDay1EtugJMbGU6dGODUKMPmIEOq3UAA86LVRdQFK86AYW3ZiW+/tOBxjvWRCO+4IbITeo9SCs89uYbdvQnGhUKhGVpZP3YmozQuuovDKfbVmZYZhrmK9MIxtlx1SyBf8g5YELn0GQBX/o66m8JMpOfTmuGYjZYQ9YUo7yDpTS++FWb8QgzOfgqU3gCbNZCaQCVj9/rZDQyf80eYALdNeArACE50A+UlNwAGfn0i09gcZNgcDHByfYCTGwOcGA6wPhxAk+sPSfGXtACOwrl4O8ayq7K1zFBwP05KazEhgKYl1iSleayQv7ZwLDizOcSJjRH2S4upschS1+WgtIREK1jLDaOoVljBXnIh2guRHuBuquQjPXfjZSc2A5D2FZnR/Th0SrfWRX9kukUXR3nh+CLplduvqES2/p3vByX7jcdJFbetYGUIJ7pRBowGwGDolpMjha1BhtPDIU4PBxgNB9hYX/Md9amaD7DK58Z53VsgpC/jgpXwnVt2/SENMxImGMuYlAZgxtrwqNPsCvcLIjvhWKCUwrO3Rtjdn2BcGJSGkWuLxLhCFaUIrFGZpS27uGAlRHlhrZgBC7Dy4iOGsRYadVRnLapIL4guvIXyj1s7K7hYdO3ji5hXmLLKghUFYAPASQAPZMDGBvDABrCxPsDGxjoeeOBMaz5BqtKTVQDX6NzfbL88yHlBcECzPW9ewQr5tKb1ffJya0HTAsOsrvIU+ovITjg2nNkc4vLeGvZyg6lhZKZuuzNWwd0mfadyjgsjuFmtGe1XaU3/mPUpTyLA+LQmh3Y7G27CqCaBDY9XsmPARpEeoym6cIVxIUtXtMfFSVDWJTY6VBveGpzYtgBsDoEHz6Z48MGzSBOFLFrc7BGoBAY0t4H6h0OjC4aXHnx7m2XXhSB8X4D7UdH4O3Cd9lSsoFX9mkrXKc64YMWy6wri0tNAbhlUGAxSLcLrOSI74diglMJDJ0bY3p1gYhmlsSgMI7WMwrC/qSoQ+UiA2VVfMkDgWmStfcMuHiPUsoMFSLv2ImutH76qlptSgDF1VBeiPWsBGwmQ2YkullwQoEZTgOE4AEwvvwnD5/wRSNWjrDjJOnMs04b3bADPOQE8/6FT1QDZiVJuSRQSrZBqhcwvWinfHhmiNxfVAc1tIMw4UYsNVT9ENysF+9hu4TigESF9yQDYAuFHi2WFRLnPrZWOClYYit26sAwqLbJEifB6jMhOOFZsrmU4c2KIvChRFApFaVFqQqabN9ZmwQo1whAGgRv3RPc4+SgvTnNaZlgiGLZg4jqS4zq9qZSv7LRu6M1YdJYBZertWHRhBJdYdOEYb78CBYDEF6wAVImuuuqoDS8M5fUsAN/54hFS7WYFD3KLRaeiodeSaDQaIE7Nhv26qtJFbrXcoCIZdixtltVQu3jFhi+lNFBKgcgVJ4UfJYoYhbEgAKkIr7eI7IRjhVIKzzmxjv1x6YtVGMb6tJl1cREBYOXbkaw7wtENezayAwCuIkHykV4lO39TtYp91OgLWqJClZDaZBUKXepILz7WOI5Z4QFN4bGP2vil7+3+PtIb+M514OQGsDYcIEuSSmq6JTrtRZckquoeoLVb4mgubot0awKDYaybYilEudX3Gb50OlKNSiex+KyP1gvjRr8hJmTaFSpBASUzyFikiYyl2UdEdsKxY2MtwwNbQ+znBfLSya40jDJRQOm7D8Rtduxuv1RFbbX81Bz5hX0GYNlWJfDVOdZFeiGlabmO6iqpsUvJdR1z14VG6i6OTKt7vJfoXnkSPKc45frzfg1r5Y/jNL2ykpqmevb2IDwnNye5RFMjfRm+i3hNLcnFUV+jGMWPWDNTldliUcTV7n7QKF4JlZrRlxPGAihhnQRBMP4alOVGZ3ahH4jshGPJmbUhdjYLTPISpVUojfuFT1AgW6cifZgHMEf3YHIpOA5Cc/ILQiMGDNi1qbXk5xYLRbraZx8RBoHGEVwYezMWXri09n6DECX5KErvvgnbW38EVsXMeQVdx7fSj2LACZ5Dr6wkF09em2hVpS4z316ndZ2+rNvmwgvPkVyjCrOVuoykt0rVVBqM05t+K0R7uQEsG1CiUBoCQVWfTegHIjvhWDIcpji7McJ0WmJix8iNi+4Yru2GtSsLIWLXB0/ZqLO4F2L1aq7PXejEDIvOfVVFeq6vF4BmpKdCoYU7T7WitdD21Ijgqn/5Swum8BFj6NN3gv4O1iYKV4b/HgVdn7GJpQLfoE/gefTq2dnafZtcGiSnCKkipFp1pi/D5+6SXJfw2unLqjIzorMNb1GkFz0nDvri6kxFdWGLew5halykBzAykgrNPiGyE44tJ9cz7EyG2M9LFD6dGSoWGaEK05f4a9/WxASqqitdajNYMER2Wqk66vN5yrhoxVg3mzlXkWFdUh/Oi+UH+O4Orc51sejivmvhWhrCIsI6vRIPFa/G+eyfdn4fU1zDn/C7MTRn8FL9NrxQ/0BVjJImteQy7aoxQ3tbPUqJuwBdRXvLCS/sh0ivNWnESomjPPe3pijKc58hzIvHIAxl8tfeILITji1aa5zeGGKSFyhK46aE8WXoCtpFeZahrIu0DLlhwnzO0qUsAbCPRLiK5KwvNnRVmORL3AHUzwGqCM5YC61r0TEYOmpjYi/E2emE6napup3MtT+pSiI0I74BTmOKa3O/lwmu4r+Y30GiCC9K/66bEimSXJjpPbxf+EzuOtrrxcJrRHOIItMWtxpgtZrwOo83JBeETQTjv9fSWEwB6ZLQE0R2wrHmxFqG3UmGiZ8CqPTpxYIYKbsxFRUA66szLQEgBRgG+RFPdJjuhwiWbSOyI6pnQg+dzQmq0QldeVHWx2yra0MtwkAVeQbp+H+F7tux4NrCezHegq/g92CRz/1eDHJ8sfwYvmvtB11/ulCF6YtUwlVYbqYxu/rYzUgOzePt9GV79Jrq88X77TTnkjKKv8W4iCc8GKI7ZpemNszIS1ulb4X7F5GdcOw5uznE3tQgLw2mlkFkQIUr5i/JdTpmAkrfTcGNDebER0yAsWAvLIra86rITrkoTfkIzUS3XMO2nl4oatML+4BPYXKIpJqP1RERVZFeu9ijLbyH6NVIWOFx/DEmCyK8fb66UHSAk1W7ClO1xBcfA1wfxna7HZH7UTFPJ0eJrBqVmq0ilUD9o8M9XPeDdN0SFDPYsC+4Eendj4jshGNPmqZ41ok119Hcsq/KBKoea8YC7LoYKNQRiLE+rQk3312jD3okP8v+hhvJD0AV+YWgrWq/a4V17ajOXVed0wz92MDNyEh1bIf9h9Sr8Xz6ASgifKb4bzHG1c7v5o9v/FM8vPFOvCx7bSW6qgAFseDQWi/Xbhe+z7jdrk13cUrn5S5Flf6Nxec3myOshKHayM+LRygtQxOkUvM+RGQnCAC21jPs5hmmpcGmsVXJPWDcQMbkRtoorUtfgtkN9aUIKG0V5ZH1kVyU1gzteu5GqnzkUKc2431GXR0IoFN+cVsdEHV5IMzILqy7pBfOfXnydvxl+UGYjrTmnr2CP935ANJE4bvXf9D1C4zb14B6LMyZ9GV9fQcKD81rj1lWgIeh7oMXi681nmYV3blrNMxu5BUR3n2FyE4QPA9urmFvbFCUBgx307bMKPyQWKgKVty29VEgawIxga1PYxpuyC20xrHLiwEI0Rga+wAq+QXa+06GdaoS8ClAnhVcWDejrGZ6UynCi5LXQGvCF/J/h32ejfBKzvFnNz+C7177wWYhTCtd2U5fLmy3I0TtdHUas4vuyK7dZtd+fM6LodUtIWy3unhoQiXAukrTHQdcWx5bSDvefYTIThA8SZLg2VtryAuDwk5BlAAEjHOLPLTjGOv7ZxFMlM6yzLCKYIx1kYJ1nc3j9jwOnd9QV1hS2IcfaWVBzNKO/GJBKtUs/AjruHtAvB1ER17k35X+IP7G+mvxu9d+uvO9d8wV/MaT/wCb+gx+YOun8bLN1zUE105bqiWEpzrSmG1iqVfHOk4MY2C2C386v0eOfgAwVSPAEKz7e5Lb4yj9HApWQnQX3gMivPsGkZ0gRJxYH+B0XmBalNgvDEZpAoKBLt3NOC9dSToAmNC+R26qoKr60jBYUaM9D0SwYJCtozyG77DuUVH7HdCdwrTcvLkyGAk1h+uKi1S6REe+75zyAzG74b9cteGGegC79sqcb4exY67gT67/FkgRvmfzdY33jfvXhXWX8GLJkU8Rd0ltEdZ3BwmFOvFILfE1HEQlNLgphdgCsAxWDKUUwECatKYIakfeIrz7gtvZf1MQ7kuevbWGrc2R61fmRTDKNNYzjWGqMEg1Bqmfw00rpNp1uk6j7cRXMGpF0ImC0lSNN1mNPRnvK9XoHzfTKTwMwtwe3YTqob3iDuTt8S01+XEtEzfMV6IJqR/6K/MVl6kivOrETyGhbOH3U/IU/+/1jwCIBYfGut2xPAgvTl0qxG2I3e+lWmIxfrSZqsO+ZxUVku2RaZgZJVsUxiIvDYyxMNbOPM8NFNAuIhLuNSSyE4QWSik8b2uEfFpgZ1JgmBCmJTBIEwAGmqyPSCyMYpBxBZvWMrQFitJFdja02bFrZ7PKR3bshqaylhENT1ylNAOz7XeoKjkJvk2J6sRnO6JrtKX5yEeHKXmIqm2t3FQ9ShFetvk6KEX4s5sfwY65imaBfs2OudIQXSUzml+gQj51WqUuVbPtro3PJrr2segy2lFU53MXyK/dsbw63nVOlMos2Y9iU5hqrNCAZVe4JEUr9y4iO0HoYDjM8NzTG/jmlR1M8hJGWRhmpImfz80AiSIUhqEUUBr4ghW4QZ8Bn8a0YKZqpBRLVEUNQVuugwODLKHRtw7NPl2qFcm0+9oBPn1HzTVF6cogOaWasovTfi/beC3+xvprAQC/8/QvYMd0pTUZ//qb78brTv8MvnfrhzoLU7ra6wio0qd1++Hsq1ejydiW2A7+0x3YxjfXg40+eM3+eIS6HdCwG7jbGK6uX6o0731EdoIwh5MbQ+xNC1zeHoPh50EDAAYGSQJFFpostCLkZN1IKZZAsJHsfPGJZShSsNVgz4TST6Yaz7Idl79XbqMQXajqDCYGWNUjp4RoTtXrWDhBaKqSW32sUdVJ9TBoRMBrTv4M/sO1f42SpzPfz3Z5Gf/35f8NpAh/0wtvJm1Js2KLt0OUF7C+0ARwPybadInqtrglTmlWffDqh+ruCVwNKq1VEOFq0qrCahHZCcICnnNqHePSAvsAuHQRELsozyqCSjW0dRWaoWAlVG5aG/rQOfGZ0LbDBAtX2elwnZZtlMWMZUdhP5xN7s7LiKK6DsEE8egomgtpRBWNBBLEFm7eQB2dvHTjtVBE+NPrv4cdc3nm+yl5iv/nyv+Bv3Xyh+ZWYs6kL1sFKkBTcsB8gXX3w1sycuti3iCa8SlRlFf3j5wtVLGm7myfaJHdvYbIThAWoJTCC0+v4WvGghnIS+NmQLAAyA0hBqUwSCmqzvSjq/iO3uzFV8vOzXenUHc7UIqblYFRHQQRVVWbtRyoIalGtWMQi1/Hoovbx+Lnw19J43EvrpdtvhYvP/Fa/E9PvA1dbXg3y8v41b9+K7aSB/DGB/8R/nYkvkZkp2oJh20QqvFIG9/7vDa89rElz5tL6/ldk752R3ndkrTMyI37W6aJ1P/dS4jsBOEA0jTFQydGeAoApgCIgRLQTGC2IB+5KSIoTY0oTyly8+R50cWyC2NkMlwRhg0zKgAwqr7DuhkPoh54xLC21YkbtVBCNEcdURzQJTl4Kcf94OrHw/aJ5AFsl7PRXbjKm+VlPPL0b0AR8IpTP9wtvGgbaEo9sExlZnzds8cWR3qHibnaI63EqcwuCmtBBki0CO9eQWQnCEuwuTHEgyBc2t53aUVbwrCfiJUBAqEEqmHCoAhaAdpaJMql6Iy1bphNdn3ySqp7QRe+bxf8a+ioMMNGabQqpUl15FYVoqCWXdwWF5gnuWrYsmi/a/uHzvxDfPLSv+psvwsUPMWnL34Y33fm9dXzdVhH1+vaBbujsnlRXZdY2uceOdKL6BppJa4Kbacy2xTG/eU0kRSt3AOI7ARhSU5tDFCYsmpPy0sDECEvXQdky27CT2i46X/8L3+tCBounRg6Q7upgOqRP8KceWFklbg7l/XvF4RmfFsgUVN2QEdU4//VlgVFx0JrYFzS3xXZfe/W66CI8B+v/C5ullcwr1vCjeIy3v/Fn8OPPffncO7MDzciu1DoMa/T97z2tmUiuHkcdjqgeVRR3gHtfa4ClwEFGVrsHkBkJwi3wLO21mHZmYgIUCXAyrqCFJCbhRyEIgz9pXybHdftYW5dpyVDijOoI8xSjmq/nroHcGIE5t+0QwTXdU4crYV9ROnLRZGd8sL7mydfByLCr3/153FzTlrzRnEZH/vW/wpFwCsfeH3VFcKNmNZ93aGTedfn6Tre9TLtYwd2Q+i8EkejP16rb57/KAcSZk+wANiwpDXvIkf65onoNBF9moi+6ten5pz3Ln/OV4noXdHx/0REXyaiv/TLs45yPYJwJzi7uYaNtQxrwxRZmmCQaCSKkGk3CommugJSq1oUdZ+2kE50r0fUnG8uHlWlGj0lpCWpPh+oxaaoLn8PbXTNQhM0+tMRNYfYotbjoZAk3q6f4z7DGx/8R0hpMPd7KniKTz79oTp6XCA6YP5j86K6LpEdqTJzSRojrfDBygujqzCA0tilniOsnqNGdu8F8CfM/H4ieq/f/8X4BCI6DeBXAJyD+3s/SkSPMPN1f8rPMPP5I16HINwxtNZ49uYIF+FHTVEEyksoA8D3swLcWIvEyveR8+nIICsmV+gS3YzrDuVNGjdsbqYYu6iLTP25SLcAACAASURBVGZv/O10JvP89CW1hBeOh/NfceqHoAj49MUP40bRHeFdzy+Dfbn+ojTevIfmtuEtKbH2aTPFKku+DqM7ZekHVznguXGBi6s+TdTR06nCrXHUmPotAD7ktz8E4K0d57wJwKeZ+ZoX3KcBvPmI7ysId5Uk0Ti7PsTGWoosS5ClCdJEI02oivIS5ZZ6pBJVzZPXFeW1I7y431q936qspGbkVnUej4QVR23Ndjh0RHrUOn/2eP18wt85/cN478t/ByfTs3O+Kcav/NW7cP7qZ+Z+l/PSl+F9uo8fnNbsjPS6rtBXwhrr2k9LYxcuhTEojRtEoLTuuQdFa+2xM8slniOslqPK7kFm/rbffgbAgx3nPATgQrT/pD8W+KBPYf73tOCnDhG9m4jOE9H5y5fnlT8Lwp0jyxKcXRthbZBibZBgkGoME+0HVlaV7FLfHSFU5cWSUypOX976EksQmJVfW3wqkhfQPO8wqcyw/rHn/tzclOb1/BI+8o1/ib+4Mis8woL0JeZUYHac3yW2eQSxlcb6wZ25Glh6GQExwpRMDOMHGCgto/AyNB0CDNFdjAjvznJgGpOI/gOAZ3c89EvxDjMzEd3qX+5nmPkpItoE8IcAfhbAh7tOZOYPAPgAAJw7d07+CxHuCQaDBGcJuEquy4FSBEYJ5QbJRGmdKHKgqmyw1p8bDQcGLJ/GjPepOmd+qq8znUkA7Gy60p2/OJXZPqaI8P1nfhiKgE8+/SFczy/NXEtup/jEhQ/i4Qde3zi+KH05L+25bFpTRd/rvEkJGt/bsu19C+4+DakxQODqB4FlmumqICnNO8eBkR0zv5GZv7dj+QSAi0T0HADw69n/yoGnADw/2n+ePwZmDusdAB8B8PDRPo4g3HkGWYLTowE2RymyRGOYurRmmCIoVW4KnUZak9w0OzNRnmqlMedM+9NOay6TzoyPh7amdiR3UCqzHeFV/eeI8MoHXo9f/Vsfxrwax2v5JfzS53+2ivA0zb/Jz01f3oIA2ack49FsYlYhmIMCsxABlparaYLayPRAd4ajpjEfARCqK98F4BMd53wKwI8S0SlfrfmjAD5FRAkRPQAARJQC+HEAXzji9QjCXWE4SHFymGFjlGI0TJGlGsMsEp6uU5o6kp7ykksOSGeqlsgWpTPnia/dfle3CS6W3DzhVSlQRNIk4FQ2r/3OCe/3vv7r+NyVz8yVzbz0ZWjPbBMLMB6pxrZMdJDa2o+vOtZiuMEDikjA4XjZNZSMsFKOKrv3A/gRIvoqgDf6fRDROSL6NwDAzNcA/CqAz/nlff7YAE56fwXgL+Givd8+4vUIwl1jkGpsDTOsDxKMBq5gZZglSBOFgdZINCGL2vKSOeKK5VaJS0Xzv8XHOyRYRX6qKbNwnMHNY6o1CwLmbMfRXCS8epBnV2D61uf/Y2RqfpeE3E7x8Qsf7Hxscfqy+/VcpMqNKK59bleb3uw5cy95IbfS7hYEHKK9cL1VB3ThtkH3YwPpuXPn+Px56a0g3HswMyZ5id2ixN64QGkspoWr3suLEpbduInWumHACj9osInaldxM3G47jJ7ipgBqjqwCdJTSY377XXgs3FTjY2G/SnuGx4kaRSRVGjPaVuSeEI+M8hdXPoNPXPggrnW038VXezo7i7/3gn+MV519A4D5qc0wEk3X5wKoIZyuc7uOtTt4x4U77cfbz48fS9TyxTHzzg/fY6KUDC3WgogeZeZzR30dGUFFEFYIEWGQajdu4hphf1oCRNClARGhNAbKEAryk7qyGz9RsStmYXaznFMkNQacTLBEKi6ut0BTHAQ3D158046lFp7fJbwwoDQQSc4fI9RT3wQefuD1ePiB1+OXPv+zC4THuJZfwu8+8esAgNc86w23JDo/EhfaFSNdruiK9LreJ9qZc81Hp2tMzTDFEbNFRkoKVm4DMnaNIKwYpRSGqUbq16OBRpYlGGUaWaKRpboqWMkSl9bMWl0VumYTv5WuCF1tck5MdVRWTf1DrTRmSI8qV0ijQ6qSXFTSPhbLss1bDkhpAnVac9lilZCy7CqLPPTIKu3HFz98JBZl0wwzpoWRLgm3AYnsBOE2oLXCWpb4qIygyKJQbgQNV5GXwFgLZRglbDVTufXjKIYJESxQpS4t88LQrjNa8f+qojwb77civyhFGUd4XVFdECaoe5qewMMPvB4KwB9f+CCu5Zcxr27/6rS772x7eqK6snJOsUrH99MlqvZzb6n97ojii0dU6cICyEuLVEtKc5WI7AThNqG1wnqWANMSgBcMCHlpQGRgrEKpXFrTWAsin8oEgXzRggXgRiGbf3MMhEcXSU/FKcvqWDNt2Sm8Oce7Jl6N31MR8KpnvQGvepZrk/vFR/9hZ1rzzGC2gjMWWqiyDHS333W3nR1mGqCD9o/KQcOMGWaQtVBMMnj0ipBvURBuIyHCC0UJaaIwSDVGg9RvJ0i0QuoXN/qK64OXqgXVmzS71P305nRTiFKPIYWZUJSejFKUYTtOXcapTgAzpf0xoTtCWyx/7wWzac1MDfC2F/yTxrG4na4tuvl97TqOzS1qaR7oSn3O2++ah+9WWW7GBBk8epVIZCcIt5kkUVjjBEQG08I4WShgkPpUJpHre2UMChNKMg3Y1kNTcahSCbfJRfc+aqyqG3WVBY3TmsDMnHhVN4KOSC6mK6gL0dy8KDRUXX78wgdxdXoZZwZn8bYX/BO82kd+geAo2+ovNy99Oa+triuqmzk2+2KzH2rFHDQXHlBPDwQilJarHyvC4RDZCcIdIEs1iFDNgJAXBooA9tEeyJezF6UTGzMMcTVCh207Lp5mJnqfxq0wkh6Rm1EdhxBeF11RXdy2t4jXPOsN+LsPvnHu4yEijDteA4u6HyzqlnDwsZn2uvZz0Bbm7PvfKsvGacYyEl+6adg16IrwDofIThDuEGmiscYAcwkk8CNpmDraoDra04pQlMZ1NzDsx3fkGdnNu2lSa6NqlQuSq3frCA6LBRcTR3UHRXO3ct480QGLOpV3vNfcApYZU91yZ/Ouy6hGQ1niR0h4vWXaYdvFLCK8wyOyE4Q7SJZqgAh708LNXp4oFMb1UbN+5A8mgtJ1e5WiqM3KC9A92FjN0BYebH3T7BLesoSoblGKs+taFs1nBywWXVf7HzC/KKU7quuS3+wTZ9Kh0babGYGqWRIOleKs/m4W2r/foh8acXQHiPAOi8hOEO4wWaJgrQbAmBog8cN35WWo4Q8jnLibWVyMlyUa1tpKcMvWLYR76K2KrQs3UspykgMOPjeeLqhrbrh294P4efPSml3v13Xuwe137vPGcqM4tXqEBj1mgKmOCt0flasJa4P8uroqiPBuHZGdINwFhlnib24GUzZgJqRaValKaxkU5scmQiVA1GnH6sASuPEwj3ZjdO/N8/OJHecfFM3FwmoXowCoOr13Me+lO6U2R4rzXjNcC/liouo5B7TfrQKGl5mfIkgRwdLs9ECG2U12fzsuooeI7AThLjEcJADVacGSLGzJSBTBAFCkomjPEXnvjtBOVZol3nvZ9Ga7e8GtiG5eWnOe1JaJ9EKnfo4GID1M+92y8BIDwAXxmdJgkOiZNk+ZD295RHaCcBcZZgmsdRO9ogSgfUdtX7hS3ZC1Qunb3Iw1t/w+S6c7w7pDVgxe+Dq30obXFll7xP+FEd2CtOayE7vGzw+TuxJ1TZy7OM15JyVTGuuHfONqtgpgtk1P6EZkJwh3mbVhBpqW2OfSH7H/f3vnHmtZVd/x72/tfe6946vyGqBCBSzV0KSh6bQxrZqW19CmiklJNBKdFB+B/mfTRglp/zBtxNqUNDGBUEQhbSNiEyWmGQKjVmtTW7QoSKQzIs2AMDOgtIDDvffs9esfa61z1tlnP9bee+3z/H2Sy5zHft5zWd/z+63f+n2nC1cUIbHvcqJGtjAT1X8NFx77whYyl1d2+EaFKt78nCPv5VYmdN889gDueeLTeG77BE7f3It3XfABvOWsy0bHbTp/50TO/d4Spab2zd/oZLPo6nuto+k6cc0MBRpFezqzTb1tUU9dynjdEbEThAVgz6b5X/FnGNp0njKOCBOpzMnCFRr9xzFOCeYhDTtMto9G8ovIm4gcUCxi2jtokRA6vnnsAdxx+G+wo7cBAM9uH8PfPfbXAIC3nHVZ4/m7fMVnUfpzev1dPsqbrbjkP1WGyQIQjNARpJdmFdIuTBAWhD2bKV65kWIjTUwkl5h0XZoo295LjQShCleOX/bTBbI/rqVYaDRXtK0/T+e2Kbu+e5749EjoHDt6G3c/fnuj+TsA0xWfBdFf0e+qUbPontAFbWuc6O1IW7FKJLIThAViNIcHxjAjMBu7FybCIDEOA4oU2DklzGBwIzIpvzJj1ar9qhaRj01ki1OQDqWo1BXhue3jwfN3mebC6ynsqTk1OVe03+zVrurT1tYeaHOQSMFKASJ2grBgjObwMLSLygk81CMRNOKgRv0V8/N3XQRwlCLNzeNV2fgUHaOuU4rWHLSdS32evrkXz24fm3r/tK29hef3xdM1ki6MbIvm9AqjuohlmB3QzEgqTi72QOVIGlMQFpA9myleMUgxSBRIGUPXNFFIEzW5zg7jKMa5G7jtRoauFT+l++QG/BABJe8aQtpgVW2XT32+64IPFLolvOcNHyzcz6H1WOiKF5VPn7uosHHa/25+QlL3WWTM5qfCfmkdkchOEBaUPZspmDVAhB0AGhopFAA9sQavbPArL2SJh4mWwgf/UUq0IuooEiZXdXn347fjue3jOG1rL97zhg/irWdfPrGdv5tfhBLcP7MgqitefN79F9r2EHVeeICd21MA24pNSWuK2AnCQvOKrQ3g5V2AzQC7M8ygNEHDuCSYXo1jG6By4UPrxeh55wEgrPHz5PnHhTVVC9Or1te95azL8Lafvzxoobm/nKFsXV6R4BZHdZPP5y0cIR9j5i1TkIXnBhE7QVhwXrE1ABHhpd0htGaogRE9ABhmPE5r2v/EmsPzB8ek5WCZn5crS63Vzd/VvV8mdKXuByVr8oqiupCClja0XboQ4oUHYGLtnQieiJ0gLAVuHR4AbO9mSJWCJmf9Mx19TaYu7QMrkHUuCfkBsUlxin8NeXFinm7yDFRHc+5YVWnPKqHr3ig6v3N8wSiyB8rjTtnk3PlilnUXPBE7QVgSRnN4YOwSsD3USBNl56ZCUpneHF4P1M3fTS1KrxGxkG2aCh3KClUCI70uc3WTKWfTFUdz4PHyvk4Yf6mhCgHWmifubZ0FT8ROEJYItyzhRTaO5hkBGesJfzrA+a2ZfbquxasTz9D5O3eckCUHTdOaQUKH4jm54EivQSEOgMKoO3e41jjBHD0psAYCxi3GfNZV8ETsBGHJcCnNl6wf3ibMYO8PrM753ODm8dREB47Q5QT+oJhffxcKc3260j9nk2rNLNdyrFToGhSqFLcPCxP0KoHz6SI2RS7nRdZARMXbrmPzaBE7QVhCjOAxsJthB7YYwYvqdK5LNJHxQ6OJlbWT82n+tj4cmmorwKU2jT9ffcoyJOLLr6PzI8YyoWtSkVm8JKFa7Iq8+CrpWWd84QMYgyQ3fwrzBSlN1meptYidICwpezZNleaLbMbO3UybORorGm6sc8JXtfwgZkorP3dXF+W0TWsW9dYsQinC1398P+567FY8+/JxnL61F+9743W45Nz9DZYkFB8731A6lK5NpMPqMQ1DrUGkpqyBGFgrtwQRO0FYYpzj+Yu7RtQU8WgAdgUMaiLiq16I3pSy9mI+pdZANQUt/nZFA7JLX9Z53339x/fjUw/fhG3bSPrEy8fwqUduglKES87ZP7V98XVOvt44kps6ZutdW+EKVZw1kBM9zQzS5U2zV4n1iWEFYUXZ2kzxqo0UW7YBsGv/pazXmUPR2EGhqJ2YG9Srfvw2Y1XtxXzyouC3FQtZclAldFXHcCJ412O3joTOsZ1t484f3Dp1vqL0ZT5tauZHu31ZmHW7Mc49ztjcB9vWYuuARHaCsAJsbZj/lRmM3aEpVlHAVDpTEY86mEwtRSgZf7vM2Y0OHRjFhWw7NlstF1j/vWdfPl64zYmTXmNpqk5fumbSUZhDEFXUQNp1V1Fk/i5WPZ0pkZ0grAhbGyles7mBQZqMIyIaRyypF/ERUa+5tFEkSBj78QWkK1VAxFfZRJqmG1GfXuCMAABn7DnT7gOkanoodPfgmknHYl5NpMtS15oZO8MMWZvuAUuEiJ0grBAbgwSv2RwgVcbsNVVqYvA3qUGbylS+MawngA0G43yKM5/eDAljQtOa7rpKt7GCnhfB973xOmwmk44Jm8kWDrzpOgAlxSc20ouRsiw49Fyoug0GsJPplXZKkDSmIKwYG4MEr8YGXtrZRaZN4+WEgXGRCk8MfNPtxcZPTNeP4qKNEEpbk42EK+BYdtkEUN5EuqpI5XfO2Q+lCHf+4FacOHkMZ+w5EwfedB0uOWd/ucs5JheqRyPkfsMO05g6GXPLEZhpJZckiNgJwgqyOTAzNC/tDkG2OtN1wt9Ik1HEUud2XlSJ2ISpnp0VVZsFJ7fzSWbbwqjDrR8sOZ6LGi85Z/9U5WWZ0BGmW5vFYv4+eNXnd38jw0yvnDWQiJ0grCgbgwQAcHJ3iCEDxDQWPTeI2ZDOd0ow65A7jvZu2YOicHHz9s0LWFET6bqOLFXvlw3kfSfxYoldGxEKvTe39m5oGxWsyrIEETtBWGE2BgkUEX62O4RmU4lJTGBi7Az1RPeRyVQmjbvxQ1V25p/QTYwHYmYGN5z/K0trag7bzl1QVbRXJXR9DuuLECUVtQ7L41duZszAiqzDE7EThBUnTRVeSSlO7mYgO2en2Td/tf0cwRNf/92gmBAwmr1qMObVBofeoveqAXgU1dWImLvm0hJ6MsUoZWLa93gezQdvBrrjuyWsiuCJ2AnCGpAkCnsAbA8zDLXtjg+Nncz0rHR9NSfTmbapMCFOfo88Sxo0KXKhevPYwLm7ImYhdF3nPifgsV2QfVp+3tH57bZMQWKZd0tYBcETsROENSFJFLaIsJNpDDMNUgoAIdNjx4TJdKaXjoRrM2aOFTTAKkBxM2EbH2RcnGKqItvNzQW9P4Pxu+2C7bwHHmDuh5sWio4+N42EafQ5U8kXBMa0F17GDOLFSMe2QcROENYIpQibZL6zm0IVxjbMIOaWJuQrNCfaaDUZ53TDgLAgrakryiJrLYMq0pb27d4LUoDm4qA1T4jb1PE6XAszwDTthae8gqLRdayYF56InSCsGUSEjdQ4nA+JbdsoPR7cvJSmm8sjotLBt4yg7V2xCcoKU6aPUed3B9QIIZl1dCHLCw4dPYg7Hr1ltD7v2ouux6XnXlm/o0dIVNeksXQfQuPOnWVjE9iylOWyeuGJ2AnCGkJESBMCkRncwAxNY4HTNqXp5vIUGJnWo/k8wM7pmQd1JzP/eE/bLEfw19xV3Vedw7nrjFLHoaMHcfNDH8d29jIA4PjJZ3DzQx8HgGDBq5rjamL0OmIGGuMaRWeZyZ+avxOafH8JrYFE7ARhjUkUgaBABOwMNTSMGKlR82iT1nRzO/nlCX3iBDesy0p1yhIYR3uhnVHuePSWkdA5trOXccejtwSJXaGLAlqKnGXWi9KHWgOk7BcNL728hNZAInaCsOYoRdhQCQAy6UxtylEY48HViF78PpFTEEbzR24RfN32dSLnpz2r5gDzTLgiBLyep8yaqIuX4KylxV2pW66iaHxfy1awImInCAIAYCNVoCFjCECBRvNITvTSRI3ErmiJQiNG9S7lLcRYl0dgQRGfFU7fMb2JWJ+x50wcP/lM4et15COeriIHIFpfzab4VZnaM39NFC3V/N3qdfsUBKE1aaKMawGNHQxSm/7zxzSXojPuCWq8n6LKH7dd6jkyFBm/FgoDGRFx56lLWaZKTaT9mnb0v/ai67GZbE28tpls4dqLrq/cj3IC60xSu0Jz8ksounLNYxPbJtHyPJHIThCEEaOiFOVFdzBzdolKgGEGjfFcXtH+DqVU60He7dbE9NWcs3ierM2A7OblGlVj+mm+GNGcx7ymx8oaSDvzV83aLGdZ8HSmiJ0gCBO4RcuaTfrKT2kmSgFeI+l8OrP12J5LawKYWuNVta+friyi7VzjpedeGb7UgIwJbFRXc+/YUayBWhyi7k6M+as2afAFFjwRO0EQpkgUAXosEk70EsXYtg2kGfkKTa88ndl4wpWMlONlCMWDo64bYskIY0il5qzSbIldutGHAeo8rYGA+gbSGTN2M41BsriCJ2InCEIhiSJwxhOyQ0QY2EIV18rKFbEgtx1R5KIAv1KzwYDaewUp7O+KezpXgzRuwKFa4ZafVKEZC20LJGInCEIpaaJMgYX3miIzsI0FzUtpen0cOzWQpupKzVBmEdUlilqvmwshZlTXZ9SlmaEYyACwbt8PtC9E7ARBqCQveK7sPx/FOPGzz8x/vUrEsO78bm1dnJiw76iub6ED5p/CBMK/s2i2llDWJWGRBE/EThCEWtLE9NIczeFRWG9JM5U32S5sVvQd1c1E6CKKxSw00zd+XTTB6/T1iYhOJaL7ieiw/feUku0OEtHzRPTl3OvnE9G3iOgIEd1NRBtdrkcQhP5IvLL+slZYedqkzWJpR29RHc1G6GLO1ZnDtT9Wk/v0v2Tongp22tA1V/BRAIeY+UIAh+zzIj4J4L0Fr38CwM3M/IsAfgrg/R2vRxCEHkkUIbEDcMg39nl9p+9NhOzygjqhO3T0IK657ypc8cU345r7rsKhowcbnyqW0LFd+O2EZ2j9DHcrftw2mTb76AZrBvNfMhZF8LqK3VUA7rSP7wTwzqKNmPkQgBf818h85bsEwBfq9hcEYXEw3UlMnJDUDMh9rOsKoY+xlch0Zakb+J1bwvGTz4DBI7eERoLXIarjnKj56WdXOVv363HbaNtiLQNjqHkkhFUL5ou8+Bah00pXsTuTmZ+2j58BUN80bsxpAJ5n5qF9/iSA15VtTEQfIqIHiejBEydOtLtaQRCiQDRu/VU1JI/W4c2YqJGdTVua5QX1/TWr3BJCqfsSkUfrscAVClHHD8E/nBNBJ34u8pu4noJfUcZxO8o0pbZAhYgeAHBWwVs3+k+YmYmotzth5tsA3AYA+/btm39MLAgCEkXYTJVZaF6yTRvj1y7EPJfvmBC6YLyrW4Kq6fvpaGIV1GdFp2/8qmylrnml2PVhXo2ja8WOmS8re4+IjhHR2cz8NBGdDeB4g3M/B+C1RJTa6O4cAE812F8QhAVAKYXNlLCTFVsAhVZuxiKK1tG0dVDovFMXt4SQ9GUTV3PvsDPBWQERjClwkkwmDxnAMNNIk9l7EHQ9470ADtjHBwB8KXRHNl9Hvgrg6jb7C4KwOChlOqsMEjU1WM+6fVSnKkznrKAm2141mW9q65bgil/K0HrsNNCICH01m+7umkTvZnrqd8eYXQs3n65idxOAy4noMIDL7HMQ0T4iut1tRETfAHAPgEuJ6Eki2m/f+giAPyaiIzBzeJ/ueD2CIMyJxBWt2AIWX/SapNHmkuTyRC5/rU198C4990p8+OIbsHfPWSAQ9u45Cx+++IbahtJl83StRc4yr0Xpo/SmswPyBG4e83c0zwnDtuzbt48ffPDBeV+GIAg52BYu+M+ZgUxr7AZ+m9eakbW2BmpW5h5iAjvMyk1kY1FkTRSrqXSMlCHVRJ1lDHLndl+GXOFSyLUR0beZeV/jk+cQ81ZBEKJhyvNp4rlShEGaIA1eiN7+/EEaSePCkzoT2Fmk24qEzq1v60qsFHLbmtp8MOXSm5nmmaczpV2YIAhRMYKHiQgPsC3HbDsp7Qa7yJml0qMtqGNCvhNNbIugeXfqKnNL0MzgjAFFUCqZybWI2AmCEB3neO6nI5UiaGsZ5Ho+JqApt4Qudglm/+VwTPCXNbjzRRXXCIUpo0O1PEzV3TCAXc2goUaa9p9kFLETBKEXlCIQT0Z4imhqPq7QLQHjqMptPbVO2huAvVKYCFdu6DOqywtdVUeStjRdmF5FnwHirtYgTb03jBaxEwShN/IpTT+6q9tvZPuD/IMq4ghGn1Fdfo6ulwKYiFEd0P5YRsCr93Vr75ip1/V3UqAiCEKvuKIVN+QFFam0OE/MyKivqM4XOrYl+X0QNaqbkTXQWPR6+t33clRBEAQP10tT2erMuvFznn6lvQy2NGmRFLsQZeJUFNZuLPh4XayBGmzHVvD6SOkCksYUBGGGJIpAGmBFU9WaPtShSKUr0TUot0atSugOHT2IOx69BSdOHsMZe87EtRddX7sYPU/sqa9ZffFwLudO8BIVNxUrYicIwkxRijAgBbDGsOIbfFO5ixEMcOTOHmVr6Ipw1kDOMcFZAwEIFrzQJtJNmGVfTcW2cTjiN42WNKYgCDOHiLAxSCrnlubR5iqazuXSlo6qObrO1kCRnc0BE9X54um+DGg9/RPji0LeSijmnKZEdoIgzI2NVI080fLD5DwymZ0LU6zgFIlO3RxdV2ugWEUp/rpHAgAut2+a3HHygVvjmBBPiWYZmhm+Q2LMj18iO0EQ5oZz/04TNVGx6d6bcwOQRpQ1kgZQ624OlFsAhVgDdU1f5s1ftY3SiNoLjuuQk3Ezl/O+ELETBGGuuHktV7GZ+GvsZpjKbLu2TqlxpWkRoa4Jba2BiuYF63CpyKEVIF2Ugoy8Vg8IdDnvqUpV0piCIMydRBHYayWmbBsxArCThQ1+XdOejXatSFfmCV1i4IpQmlRj5ntr1rEo7uaOQpfzns4rYicIwkKQJsosKrbPTaRnBr5Raq3H89cKgO25WWcJ5NN0Ld2l514ZvtTAFsGE0KbvZgyxa3II53KudQY1SKJHlSJ2giAsDImiqWKVRJlSdFe4UOaY0MvavJZuCUD8ZQz56wrxl2vbXHrW7vI+Gqa9HIGj9ssUsRMEt+CrhQAADHpJREFUYWEwDZIxIXimUGX8vNIxoSUx3RIcfXVICRG6rh1aYmlMJx88qm480BQRO0EQFooiwStyS3Db5h0TfOGrGir9IhitGRQxCpun0HVut9VDYUpTynzwuiBiJwjCwpEXvCZuCaZkfmwVFEJsaeqrv2aV0MXyw1uEJtJ9fFWQpQeCICwkfvNooD+3BCCuOPUR1bn1iFXnjCF08ZtIt0MaQQuCsHYkikzVgkJtdNe2RiXW0NpHUUrVOrrY7gnxm0gvTlsAiewEQVh4EkVIqN7Net6Da+ygrqi/piO20EWP6jocqo80pkR2giAsBc4tgVkXFqs4mrslxBtaYx2LaoQ91vzc+ITha/bCD7k4UR0gYicIwhJBRNgcJNjZzUrtgVyRyqyJ1eaqrv1XUbVlVx+8mEUpjq6H9AuNYiBiJwjC0jFIFZBp03EjN/Ar6sGANYCukVZdNAeUC10XH7zY6UvHPCyaqpA5O0EQlg4jDAqJIqS5SGge83ZdIkkncm2EDujog9dD+hKYnbt5EySyEwRhKXHNo0GEhADF3ro8orhzWjW0iSSJwvtsVlkEdfHB6yN9CSxeVAdIZCcIwhKTJmqiE4rzxetrEC+jSWTnR3IhQldnEdTWBy/0/G1YPKkTsRMEYcnJp+GICINUYZDzxuuLIKGjse9dE5EJWV7Qxgevq9lrFS6qc555zhi27Mdt02vjbEgaUxCEJcd0F8FU02BFBKh6t4SulB3O9e3sktILWUfX1AevqQdeHflm3Iki6DbH5/EDAqAoidqvVMROEISlp0jw8lWZpW4JALosY9Y9OCYAzdqOhfrghVR8hlBqAhupiXQfX0pE7ARBWAnygpe3Bspv67sluG1C3BLMHuPj9EFVQUpbugpdiMt5rIixj1+riJ0gCCtDXvASFeaJ5rquNHZL6GmOKXZU00XomnRrWcQqTIeInSAIK4UTvEybZQkmndmPKPVx2NiuCXUdWcpo2pJMRVyv10dZkVRjCoKwcrhlCKbQIWz7NsTWutjpy6pG0lXXMMx0sy8IHQtxpg4naUxBEIRw0kRZr7fqwXtRfNeiRaBkFow3EfEuLgqx05d9JENF7ARBWGkSRQATtrMKsWvpgxeTWI2k28zPlbUiCzthD2LXQ2gnYicIwsqTJApbRNgpSc9RS7WLqY+do7oZR3OOWXeraYuInSAIa4FShJQJmUahH948g7uuUd3MozlLH51Y+qroFLETBGFtSJRZU6dAdi7PX4ROlaaweWLO17WN6po0k/YZZrrV+SZP3o8w9RUoitgJgrA2kHVIyGxbq8QTvXnN27WK6lqkLIHptGUX09dlc0wQsRMEYa1QisB6HE0liqAY0EzIdDZzvWsS1bWN5IBioWtr+tpXI+k+p/9knZ0gCGtHoiaXLbs5r400mYlTgqOJY0ITW6Ci8+QLUdqavsZuJO3TZwcWiewEQVhL0kRhmOmJSE5ZYVFes2jNxf01Y1CawbSNpdtGcT5lFZdtTF9jNZIuPX5vR5bIThCENabIC89v8uw86Jw3nooc9Y0iOxqfL1GEVDXzvas6ftnSgsamrzT9+4oJRXJMKEPEThCEtcX00ZwcYMtSaSMhsuKXeqlF5YSQyued3HtuW2XFw7irq9FxYg34dWvoGpm+EpCqfuWi7ybSksYUBGGtyTslKEXIKrqtOJRSUGxL+D27oFC07rcgo26xeLDp6wyErqvJbQgidoIgrD1TgkfUm1OCo8/jh66jqzV9nYHQAbOxBhKxEwRBwKQ1UKIIOiC667I0ry+pi2UR1Hcxio+icUEQgCkz3RiI2AmCIFiMNZBZaJ407KjShL5MX41gdD92Ww+8UJywMUwF7G6Mji41iNgJgiDkMBGNgh5WLzInG5E0pa8MZoyori+h09os4Zj4fRFANJs6SRE7QRCEAhJF2EgUdnJr8Xxa++C1vagKOgtdyxZkVTAzNJdHsrOYq3OI2AmCIJSQJAoDmMKVYmsgtFKu2GnMrunL2NGcLvl9TTCDCkwfETtBEIQKRk4JbG1xvPcWwQcP6BDVRY7mgkTOMmsfPBE7QRCECkZOCQDShKA1TxSuzNvkvK1rglvcHusamiyl6KuRdOU5u+xMRKcS0f1EdNj+e0rJdgeJ6Hki+nLu9c8S0Y+I6CH7c3GX6xEEQegD5TWOVoowSNRIKNrY7MSk0Xo92/szVSqK0DEzhiXu75XXMAd3865lMB8FcIiZLwRwyD4v4pMA3lvy3p8y88X256GO1yMIgtAL+TVniRW9dIYuCXmC05eRRc6du036dNbpS0dXsbsKwJ328Z0A3lm0ETMfAvBCx3MJgiDMDZPOnB6o01QhtY2iQ4bxmIFdXZToFobHFDlmxn1P/DPec/AduOKLb8Y1912FQ0cPBu0bo7l1W7rO2Z3JzE/bx88AKGmXXclfEtGfw0aGzLxdtBERfQjAh+zTbSJ6pMW5VoHTATw774uYI+t8/+t878B63//C3Hvy6uTUwamD14NMsHQYh/Fv/K969ye7/5O9kP2kh1O+McZBasWOiB4AcFbBWzf6T5iZiajpd5YbYERyA8BtAD4C4GNFGzLzbXYbENGDzLyv4blWgnW+d2C973+d7x1Y7/tf93uPcZxasWPmyyou4hgRnc3MTxPR2QCONzm5FxVuE9FnAPxJk/0FQRAEIYSuc3b3AjhgHx8A8KUmO1uBBJkk7jsBrGtqUhAEQeiRrmJ3E4DLiegwgMvscxDRPiK63W1ERN8AcA+AS4noSSLab9/6ByJ6GMDDMDnpvwg8720dr3uZWed7B9b7/tf53oH1vn+5945QX923BUEQBGFRmE27aUEQBEGYIyJ2giAIwsqzsGK3zq3IItz7+UT0LSI6QkR3E9HGbK48Dg3u/4Dd5jARHfBe/xoRPeZ99ntnd/XtIKIr7TUfIaKpTkREtGk/yyP2sz3Pe+8G+/pj3nz40tD23onoPCI66X3Ot8762rsScO9vI6LvENGQiK7OvVf4979MdLz/zPvs7609mbOGWLQfAH8Fs8gcMG3IPlGy3aUA3g7gy7nXPwvg6nnfx5zu/fMA3m0f3wrg+nnfU+z7B3AqgMftv6fYx6fY974GYN+876PB/SYAfgjgApg1p98FcFFumz8CcKt9/G4Ad9vHF9ntNwGcb4+TzPueZnTv5wF4ZN730PO9nwfgVwDc5Y9nVX//y/LT5f7tey82Od/CRnZY71Zkre/dLuO4BMAX6vZfYELufz+A+5n5J8z8UwD3A7hyRtcXm98AcISZH2fmHQCfg/kd+Pi/ky/AVDaTff1zzLzNzD8CcMQeb1nocu/LTu29M/MTzPw9ADq37yr8/Xe5/8YsstjFakX2PSK6mYg2I15b33S599MAPM/MQ/v8SQCvi3lxMyDk/l8H4Kj3PH+fn7HpjT9bgoGx7l4mtrGf7f/CfNYh+y4yXe4dAM4nov8ion8horf2fbGR6fLZLfvnDnS/hy0iepCI/p2Iar/Qz9XPjhakFdk86PneF56e7/8aZn6KiF4N4J9gHDfuanelwgLzNIBfYObniOjXAHyRiH6Zmf9v3hcmzITX2//PLwDwFSJ6mJl/WLbxXMWO17gVWY/3/hyA1xJRar8FnwPgqY6XG50I9/8UgN/2np8DM1cHZn7K/vsCEf0jTLpkkcXuKQDnes+LPjO3zZNElAL4OZjPOmTfRab1vbOZuNkGAGb+NhH9EMAvAYjSS3EGdPnsSv/+l4hOf7ve/+ePE9HXAPwqzBxgIYucxlznVmSt790OAF8F4CqXGv/uFoCQ+78PwBVEdIqt1rwCwH1ElBLR6QBARAMAv4/F/+z/E8CFtop2A6YII19d5v9OrgbwFftZ3wvg3bZi8XwAFwL4jxlddwxa3zsRnUFECQDYb/cXwhRqLAsh915G4d9/T9fZF63v3973pn18OoDfAvBo5U7zrsipqNQ5Dcb25zCABwCcal/fB+B2b7tvADgB4CRMzne/ff0rMG3IHgHw9wBeNe97muG9XwAz4B2BadO2Oe976un+r7X3eATAH9rXXgng2wC+B+D7AP4WS1CdCOD3APw3zDfTG+1rHwPwDvt4y36WR+xne4G37412v8cA/O6872VW9w7gD+xn/BCA7wB4+7zvpYd7/3X7//ZLMJH89719p/7+l+2n7f0D+E07vn/X/vv+unNJuzBBEARh5VnkNKYgCIIgREHEThAEQVh5ROwEQRCElUfEThAEQVh5ROwEQRCElUfEThAEQVh5ROwEQRCElef/AaX4RrlTKOa5AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -297,7 +297,7 @@ "\n", "sim = setup_sim()\n", "patches = [[] for i in range(sim.N)]\n", - "steps = 200\n", + "steps = 300\n", "xyz = np.zeros((steps,sim.N,3))\n", "for k in range(steps):\n", " sim.step_manual(\"ias15\")\n", @@ -322,7 +322,7 @@ }, { "cell_type": "code", - "execution_count": 212, + "execution_count": 340, "metadata": {}, "outputs": [], "source": [ @@ -332,11 +332,11 @@ " Omega2 = np.concatenate((-np.identity(Ndim//2),np.zeros((Ndim//2,Ndim//2))),axis=1)\n", " Omega = np.concatenate((Omega1,Omega2))\n", " J = np.zeros((Ndim,Ndim))\n", - " delta = 1e-9\n", + " delta = 1e-6\n", " sim_0 = sim.copy()\n", - " sim_0.move_to_hel()\n", + " #sim_0.move_to_hel()\n", " sim_0.step_manual(algo)\n", - " sim_0.move_to_com()\n", + " #sim_0.move_to_com()\n", " for i in range(Ndim):\n", " sim_p = sim.copy()\n", " if i" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ - "fig, ax = plt.subplots(2,1,figsize=(7,10))\n", + "fig, ax = plt.subplots(3,1,figsize=(7,15))\n", "ax[0].set_yscale(\"log\")\n", "for i,algo in enumerate(algos):\n", " ax[0].plot(times,np.maximum(energy[i,:],1e-16),label=algo)\n", @@ -610,7 +433,12 @@ " if i>0:\n", " ax[1].plot(times,simdiff[i-1],label=algos[0]+\"-\"+algo)\n", "ax[1].legend();\n", - "ax[1].set_ylabel(\"difference between integrators\");" + "ax[1].set_ylabel(\"difference between integrators\");\n", + "ax[2].set_yscale(\"log\")\n", + "for i,algo in enumerate(algos):\n", + " ax[2].plot(times,symple[i],label=algo)\n", + "ax[2].legend();\n", + "ax[2].set_ylabel(\"symplecticity error\");" ] }, {