{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Introduction to MulensModel\n", "\n", "How to create and plot a model and then add some data and fit for the source and blend fluxes.\n", "\n", "This example shows OGLE-2003-BLG-235/MOA-2003-BLG-53, the first microlensing planet. See [Bond et al. 2004](http://adsabs.harvard.edu/abs/2004ApJ...606L.155B). The data were downloaded from the [NASA Exoplanet Archive](https://exoplanetarchive.ipac.caltech.edu/cgi-bin/DisplayOverview/nph-DisplayOverview?objname=OGLE-2003-BLG-235L+b&type=CONFIRMED_PLANET)." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Import basic packages\n", "import MulensModel\n", "import matplotlib.pyplot as plt # MulensModel uses matplotlib for plotting.\n", "import os.path" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Define a point lens model:\n", "my_pspl_model = MulensModel.Model(\n", " {'t_0': 2452848.06, 'u_0': 0.133, 't_E': 61.5})\n", "\n", "# Or a model with 2-bodies:\n", "my_1S2L_model = MulensModel.Model(\n", " {'t_0': 2452848.06, 'u_0': 0.133, 't_E': 61.5, 'rho': 0.00096, \n", " 'q': 0.0039, 's': 1.120, 'alpha': 223.8})\n", "# Since rho is set, define a time range and method to apply \n", "# finite source effects:\n", "my_1S2L_model.set_magnification_methods(\n", " [2452833., 'VBBL', 2452845.])" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deZyNdf/H8dfHzDBkjamxjbFkHWOK\nyc4gEUopS7eESD9bN0ISlaQid7m1WJLKEtFKkSVrdyT7NolkjRg7Y5nl+/vjOjONMcw5M3PmOmfO\n5/l4nMec6zrX8pljvM/3fK/r+l5ijEEppZTvyGV3AUoppbKXBr9SSvkYDX6llPIxGvxKKeVjNPiV\nUsrH+NtdgDOKFStmQkND7S5DKaW8yqZNm2KMMUGp53tF8IeGhrJx40a7y1BKKa8iIgfTmq9dPUop\n5WM0+JVSysdo8CullI/xij5+pZRyRVxcHEeOHOHKlSt2l5ItAgMDKVWqFAEBAU4tr8GvlMpxjhw5\nQoECBQgNDUVE7C7HrYwxnDp1iiNHjlC2bFmn1tGuHqVUjnPlyhWKFi2a40MfQEQoWrSoS99uNPiV\nUjmSL4R+Eld/Vw1+5fW2bt3KL7/8YncZSnkNDX7l9e6++27q1KljdxlKXcfPz4+IiAjCwsJo3749\nsbGxAIwZM4Zq1aoRHh5OREREcqMlKiqKSpUqUaNGDerXr8+ePXuS52f1Bawa/Eop5QZ58+Zl69at\n7Ny5k9y5czN58mTWrVvHd999x+bNm9m+fTvLly+ndOnSyevMnj2bbdu20bVrV4YMGeK22jT4lVLK\nzRo2bMi+ffs4duwYxYoVI0+ePAAUK1aMEiVK3LB8o0aN2Ldvn9vq0eBXSuV8UVHwySfW87g4a3rW\nLGs6Ntaa/vxza/rcOWv6q6+s6ZgYa3rhQmv6+HGXdh0fH8/ixYupXr06999/P4cPH6ZixYr06dOH\n1atXp7nOwoULqV69ukv7cYXbgl9EpovICRHZmWLeWyLym4hsF5GvRaSwu/avlFJ2unz5MhEREdSq\nVYuQkBB69OhB/vz52bRpE1OnTiUoKIiOHTvySdIHEtC5c2ciIiL43//+x/jx491Wmzsv4PoEeA+Y\nkWLeMuAFY0y8iIwFXgCed2MNSikFq1b98zwg4PrpfPmuny5U6PrpYsWunw4OdmqXSX38qfn5+REV\nFUVUVBTVq1fn008/pVu3boDVx1+rVi2ntp8ZbmvxG2PWAKdTzVtqjIl3TK4HSrlr/0op5Wn27NnD\n3r17k6e3bt1KmTJlsr0OO4dseAr43Mb9K6VUtrp48SL9+/fn7Nmz+Pv7U6FCBaZOnZrueq1bt04e\nh6du3brMnz8/U3XYEvwi8iIQD8y+xTK9gF4AISEh2VSZUkpljYsXL94wr2bNmvz8889pLr8qZXeS\nE/MzI9vP6hGRrkAboLMxxtxsOWPMVGNMLWNMraCgG+4cppRSKoOytcUvIi2xDuY2NsbEZue+lVJK\nWdx5OuccYB1QSUSOiEgPrLN8CgDLRGSriEx21/6VUkqlzW0tfmPM42nM/shd+1NKKeUcvXJXKaV8\njAa/Ukr5GA1+pZRyg6eeeoo77riDsLCw5Hnr16+ndu3aREREUKVKFV555RXAumI3PDyc8PBw6tWr\nx7Zt25LXyZ8/f5bXpvfcVUopN+jWrRv9+vXjySefTJ7XtWtX5s2bR40aNUhISEgec79s2bKsXr2a\nIkWKsHjxYnr16uXWmwtp8CullBs0atSIAwcOXDfvxIkTFC9eHLDG7KlatSoA9erVS16mTp06HDly\nxK21afArpXK0AQMGpDlYWmZEREQwYcIEl9cbOHAglSpVIioqipYtW9K1a1cCAwOvW+ajjz7igQce\nyKpS06R9/EoplU1eeuklNm7cyP33389nn31Gy5Ytr3t95cqVfPTRR4wdO9atdWiLXymVo2WkZe5O\n5cuXp3fv3jz99NMEBQVx6tQpihYtyvbt2+nZsyeLFy+maNGibq1BW/xKKZVNvv/+e5KGKNu7dy9+\nfn4ULlyYQ4cO0a5dO2bOnEnFihXdXoe2+JVSyg0ef/xxVq1aRUxMDKVKlWLUqFEsW7aMgQMHki9f\nPvz9/Zk9ezZ+fn68+uqrnDp1ij59+gDg7+/Pxo0bAYiNjaVUqX9uXTJo0CAGDRqUqdrkFgNkeoxa\ntWqZpDdBqdREBABv+FtW2SM6OpoqVarYXUa2Sut3FpFNxpgbbumlXT1KKeVjNPiVUsrHaPArpXIk\nX+r6c/V31eBXSuU4gYGBnDp1yifC3xjDqVOnbrgQ7Fb0rB6lVI5TqlQpjhw5wsmTJ+0uJVsEBgZe\nd+ZPejT4lVI5TkBAAGXLlrW7DI+lXT1KKeVjNPiVUsrHaPArpZSP0eBXSikfo8GvlFI+RoNfKaV8\njAa/Ukr5GA1+pZTyMRr8yqv5wiX5SmU1DX7l1RITE+0uQSmvo8GvvJoGv1Kuc1vwi8h0ETkhIjtT\nzLtdRJaJyF7HzyLu2r/yDQkJCXaXoJTXcWeL/xOgZap5w4AfjTF3AT86ppXKMA1+pVzntuA3xqwB\nTqea3Rb41PH8U+Bhd+1f+Yb4+Hi7S1DK62R3H/+dxphjAI6fd2Tz/lUOoy1+pVznsQd3RaSXiGwU\nkY2+cjMF5ToNfqVcl93B/7eIFAdw/DxxswWNMVONMbWMMbWCgoKyrUDlXTT4lXJddgf/AqCr43lX\n4Nts3r/KYTT4lXKdO0/nnAOsAyqJyBER6QG8CTQXkb1Ac8e0Uhmmwa+U69x2z11jzOM3eamZu/ap\nfI8Gv1Ku89iDu0o5Q4NfKddp8CuvpsGvlOs0+JVX0wu4lHKdBr/yatriV8p1GvzKq2nwK+U6DX7l\n1TT4lXKdBr/yahr8SrlOg195NQ1+pVynwa+8mga/Uq7T4FdeTYNfKddp8CuvpsGvlOs0+JVX0+BX\nynUa/Mqr6ZW7SrlOg195NW3xK+U6DX7l1TT4lXKdBr/yahr8SrlOg195tbi4OLtLUMrraPArt4mJ\nieH333936z6uXLni1u0rlRO57daLSlWsWJEzZ85gjHHbPq5eveq2bSuVUzkV/CJSERgClEm5jjGm\nqZvqUjnAmTNn3L4PDX6lXOdsi38+MBn4ENCjacpjaPAr5Tpngz/eGDPJrZUolQFJwZ8rlx6uUspZ\nzv5vWSgifUSkuIjcnvRwa2VKOSEp+P399XCVUs5y9n9LV8fPISnmGaBc1pajlGuSgt+dB5CVymmc\nCn5jTFl3F6JURmgfv1Kuc/asngCgN9DIMWsVMMUYo1fPKFtpi18p1znb1TMJCAA+cEx3cczr6Y6i\nlHKWBr9SrnM2+CONMTVSTK8QkW3uKEgpV8TGxtpdglJex9mzehJEpHzShIiUIxPn84vIQBHZJSI7\nRWSOiARmdFvKt128eBHQFr9SrnC2xT8EWCki+wHBuoK3e0Z2KCIlgWeBqsaYyyIyD+gEfJKR7SnP\nl5iY6Lbz7C9dugRo8CvlCmfP6vlRRO4CKmEF/2/GmMycTuEP5BWROCAf8FcmtqU8XEJCgtuCP6nF\nr5Ry3i2DX0SaGmNWiEi7VC+VFxGMMV+5ukNjzFERGQ8cAi4DS40xS9PYdy+gF0BISIiru1EeJCEh\ngYCAALdsW1v8SrkuvWZYY8fPB9N4tMnIDkWkCNAWKAuUAG4TkSdSL2eMmWqMqWWMqRUUFJSRXSkP\n4c6bpWiLXynX3bLFb4x52fH0VWPMnylfE5GMXtR1H/CnMeakYztfAfWAWRncnvJwiYmJbtu2Br9S\nrnO24/XLNOZ9kcF9HgLqiEg+ERGgGRCdwW0pL+CuFn9iYiLnzp1LntbuHqWck14ff2WgGlAoVT9/\nQSBDp2AaY34RkS+AzUA8sAWYmpFtKe/gruC/cOECxhgKFizI+fPnMcZgtSWUUreS3lk9lbD68gtj\n9esnuQA8ndGdOrqQXk53QZUjuCv4z549C0CRIkU4f/68W/ahVE6UXh//t8C3IlLXGLMum2pSOYy7\ngj/pDl9FihTh4MGD2tWjlJOcvYBri4j0xer2Se7iMcY85ZaqVI7iruA/ffo0ALffbt0aIsuDf8cO\n+PtvuO++rN2uUjZz9uDuTCAYaAGsBkphdfcolS53ndVz/PhxAIoXL551G71y5Z/nY8dCr17/TA8e\nDCNHZt2+lLKJs8FfwRgzErhkjPkUaA1Ud19ZKidxV4v/2LFjwD/Bn+kW/5dfQunSsH+/NT1qFCxc\n+M/rly6B48NGKW/mbFdP0rj7Z0UkDDgOhLqlIpXjuCv4jx8/TmBgIIULFwayIPjr1IGWLcHPz5ou\nX/761ydNgqRvL3v3wsSJMH485MmTuf0qlc2cbfFPdVxxOwJYAOwGxrmtKpWjuDP4g4ODM3cK5/nz\n8O67YAyULAkzZ0KZMjdfPmnMoWXLYP58OHgw4/tWyibODtI2zfF0DXqfXeUid3b1BAcHJ09nqMX/\n6acwaBA0agQ1atzwsjGGo0ePcvDgQfz8/ChfvjxBQUHQpw906QIFClgLJib+86GglIdz6i9VRF4X\nkcIppouIyGvuK0vlJO46uJsU/Ekt/gwFf79+sHHjDaF/4cIFxowZQ4UKFShdujQNGjSgbt263HHH\nHdxzzz1MmjSJq7lzWwtPnAidOv3TDaSUh3O2ifKAMeZs0oQx5gzQyj0lqZzGHS1+YwwHDhwgNDQ0\nY109b70FR4+CyA2hv3DhQipUqMCIESMoW7YsEydOZPHixXz33Xe8+eab+Pn50adPHypVqsTy5cvh\n6lUr9K9dy6LfTin3cjb4/UQk+QiWiOQF9IiWckp8fHyWb/P48ePExsZSoUKF5HlOt/gPHoRXX4Xp\n06+bbYzhhRde4KGHHqJkyZJs2LCB5cuX079/f1q2bEnr1q15/vnn2bBhA0uXLiUwMJDmzZszLCaG\nxLlzIVBvJKe8g7PBPwv4UUR6iMhTwDLgU/eVpXKCpJa4O4J/3759AJQvX971rp4yZWDbNnjhheRZ\niYmJ9O7dmzfffJOnn36an3/+mcjIyDRXFxGaN2/Oli1b6NWrF2PHjaPj449z7dAhaN9eD/gqj+fs\nwd1xIrIDayRNAUYbY5a4tTLl9fz9/YmLiyMuLi79hV30xx9/AFChQgW2b9/u3EqXL8OaNdCiBZS7\n/hyFoUOHMmXKFIYNG8brr7/uVPdR3rx5mTx5MhUrVmTw4MHEnzvHvG3bCNi+/dZnBillM2fP48cY\nsxhY7MZaVA7jzuDfvXs3uXPnpkyKgE23xT9hArz4IuzeDZUrJ8/+8MMP+c9//kPfvn2dDv0kIsJz\nzz1Hnjx56N+/P0+2b8/s1q2d/iqtlB3SG5b5J2NMAxG5AKT8XyWAMcYUdGt1yqv5OS6EckdXz+bN\nm6levToBAQHOd/UMGgRhYdeF/rp16+jTpw8tW7ZkwoQJGb4moF+/fsTGxvL8889ToVIlRjdoAMWL\nQ3h4hranlDul1+J/EsAYUyAbalE5jL+/9eeV1S1+Ywxbtmzh0UcfBUg/rK9csc6xz5MHHvxndPHz\n58/TuXNnSpUqxdy5c5PrzaghQ4awd+9eXnvtNaoUKsS/WraEuXMztU2l3CG9b6TzAUTkx2yoReUw\n7gr+Q4cOcfr0ae6+++7r5t+0xT9yJERGQmzsdbP79+/PwYMHmTVrFoUKFcp0XSLC+++/T+PGjXnq\n8mW2DhyY6W0q5Q7pNXFyicjLQEURGZT6RWPM2+4pS+UEScGf1V09mzZtAuCee+4BSL+rp3Fj61TL\nfPmSZ/3www/MmDGDkSNHUr9+/SyrLXfu3MyfP5+IiAg6Pvkkm9avJ39cHNxxR5btQ6nMSq/F3wm4\ngvUBUSCNh1I35a4W/9q1awkMDKSG48KrdLt62rSB0aOTJy9fvkzfvn2pVKkSL774YpbWBhAUFMRn\nn33Gvn376FulCnToYI0FpJSHSO8OXHuAsSKy3XFWj1JOSzq4m9XBv2LFCurXr09gqgumbmjxf/yx\ndVVtr17XjaPz+uuvs3//flasWEEeN42s2bhxY1566SVeeeUVmnfqxBNu2YtSGXPLFr+IJP29VhWR\nQakf2VCf8mLu6Oo5efIk27dvp1mzZsnzbtrVs3ChNcZ+im8E+/btY+zYsXTp0oUmTZpkWV1pGTFi\nBA0aNKDfJ59w5OhRt+5LKVek19Vzm+NnfrSrR7nIHV09S5cuBaBp06bJ827a1fPllzcE/4gRIwgI\nCGDs2LFZVtPN+Pn58cknnxAXF8fTDzyASXGlsFJ2Sq+rZ4rj56jsKUflJO7o6vn6668pXrx4msMp\nJLf4f/sNgoOhcGEo+M+lJr/++iuff/45I0eOzNrbNd5C+fLlGTduHP369eOj+Hh6jh4NmTxtVKnM\ncuovUESCgKex7rqVvI7ebF3dSlZ39Vy+fJnFixfTrVs3cqXos7+uq8cYa4jkPHlg/frk1r4xhuef\nf56goCAGDx6cJfU4q3fv3nz1xRcM2rSJ5kePXne1sVJ2cPbK8m+BQsBy4PsUD6VuKqu7ehYsWEBs\nbGzyhVtJruvqEbFG3Rw79rouniVLlrBy5UpGjhxJwYLZe8F5rly5+OjjjzHG0OuppzBr12br/pVK\nzdnvnPmMMc+7tRKV42R18E+fPp2QkBCioqLSfN0k3QjFcX5/8nxjeOmllwgNDeWZZ57JklpcFRoa\nyhtvvEH//v35bMMGOp84AXnz2lKLUs62+L8TEb3xinJJVgb/oUOHWLZsGd27d7+umwdSdPX06wcv\nv3zDusuWLePXX39l+PDh5E66a5YNevfuTZ0aNRjg70/MpUu21aGUs8H/b6zwvywi50Xkgoicd2dh\nyvslBfKVK1cyva0PP/wQgK5du950P/j5Xde9k+S1116jVKlSPPnkk5muIzP8/PyYOnMmZy9etI4z\nuOlexEqlx9nx+PXUTeWypLNsLl68mKntXLp0iQ8++IC2bdtStmzZm+9vwgQoVuy6eatXr2bt2rW8\n++67brtYyxXVq1dn6NChvP7663TZsIFmO3ZYH1hKZSNnb7Z+TxqP8iKSofPSRKSwiHwhIr+JSLSI\n1M3IdpRnSwr+S5ns1pg+fTqnT59myJAhab4uu3b9s79ULf7XXnuNO++8kx49emSqhqw0cuRI7goO\n5pkjR7h8+rTd5Sgf5GxXzwfAeuBDx2M9MBf4XUTuz8B+/wv8YIypDNQAojOwDeXhsiL4r127xttv\nv029evWoV6/ejQvExSE3Gfp48+bNLF++nOeee468HnQgNTAwkCmffcYfFy7w6ts6zqHKfs4G/wHg\nbmNMTWNMTSAC2AncB4xzZYciUhBoBHwEYIy5Zow568o2lHfIiq6eadOmceDAAUaMGJH2AgEByNCh\ngHXf3JTeeecd8ufPT69evTK8f3dp0qQJ3bt356233mKbXtGrspmzwV/ZGLMracIYsxvrg2B/BvZZ\nDjgJfCwiW0RkmojclnohEeklIhtFZOPJkyczsBvlKTLa4r948SKvvvoqjRs3pmXLlmltGIDcjiGP\nU549dPToUebOnUuPHj2yZKx9dxg/fjxF8+Sh57hxJBw7Znc5yoc4G/x7RGSSiDR2PD7A6ubJA7h6\nrp4/cA8wyRhzN3AJGJZ6IWPMVGNMLWNMraCgIBd3oTxBUov/woULGVp//Pjx/P3337zxxhs3jseT\nkABRUdC/PwEBAYDVLZTk/fffJyEhgWeffTZD+84Ot99+O/99/302Jiby7uef212O8iHOBn83YB8w\nABgI7HfMiwNcHeLwCHDEGPOLY/oLrA8ClcMkBX9GvrHt27ePN998k44dO1K3bhrH/hMToW1baNgw\n+dz8pOC/dOkSkydP5pFHHqFcuXIZ/wWyQceuXWnVqhUvvvgiB3791e5ylI9wKviNMZeNMf8xxjxi\njHnYGDPeGBNrjEk0xrjUgWuMOQ4cFpFKjlnNgN0u1q28QFLwHz9+PP0boadar1+/fuTOnZu3b3bw\nMyAARoyADh1uCP4ZM2Zw5swZBnrBrQ9FhEmTJiHx8fxfvXqYv/+2uyTlA5w9nfMux+mXu0Vkf9Ij\nE/vtD8wWke1YB4pfz8S2lIdKCvtr165x5swZp9ebM2cOS5Ys4bXXXqNEiRI3LjB+PPz0U/JkyuA3\nxjBx4kRq1aqVpbdUdKeQkBDeGDKEJfHxfLZokd3lKB/gbFfPx8AkIB6ra2cGMDOjOzXGbHX034c7\nvkE4nwrKa6Rs5R8/ftypdQ4dOkSfPn2oV68effr0uXGBS5fg/fdh/vzkWSmDf/Xq1fz222/07ds3\n/VsyepA+o0ZRu3ZtBgwdSkxMjN3lqBzO2eDPa4z5ERBjzEFjzCtA03TWUT7OGEP+/PkBOObEWSsJ\nCQl069aNhIQEZs6cmTzWz3Vuuw127IAxY5JnpQz+SZMmUaRIETp27Jg1v0Q28fPzY9q0aZw9c4ZB\n994LWXyDeqVScjb4r4hILmCviPQTkUeAO9xYl8oBjDHJNzxxJvhffvllVq5cycSJE9M+KLttm3VQ\nN39+6+GQFPyHDx/mq6++olu3bh51wZazwsLCeOHhh5n5558smZnhL9RKpcvZ4B8A5AOeBWoCXYAb\nR8tSKpUyZcqQK1cu9uzZc8vlvvzyS8aMGUPPnj3p1q3bjQscOwb16kEaFzslBf/kyZOJj4+3bejl\nrDB85kwqVazI/736aqaHulDqZpw9q+dXY8xFY8wRY0x3Y0w7Y8x6dxenvJsxhnz58lGtWjV+vcWp\niv/73/948sknqVOnDu+9917affPBwTBlCqTR758U/OvWraNp06ZUqlTphmW8RWDevHzouFr5pe7d\n7S5H5VC3HGRNRBbc6nVjzENZW47KSYwxiAiRkZF88803JCYm3jCW/tatW2ndujWlSpXim2++SXsE\nzaTB1554Is39pBxjv3fv3ln6O9ihYcOG/F+DBkyYP59OM2YQafNw0irnSa/FXxcoBawFxgP/SfVQ\n6qaSgr9Zs2acPn2aVatWXff6ypUradKkCQULFmTZsmXceeedN27k7FmoVQsWL77pfpKCPzg4mLZt\n22blr2CbN+fPp/jtt9P1zTe5fPmy3eWoHCa94A8GhgNhWCNqNgdijDGrjTGr3V2c8m5Jwf/II49Q\npEgRRo8eTVxcHFevXmXMmDG0aNGCEiVKsGbNGkJCQtLeyJkz1sVaqcbZTynpzKGePXsmD9/g7QoF\nBzN9zhyio6MZqYO4qaxmjHHqAeTBGqbhJNDf2fWy4lGzZk2jvE+1atVMu3btjDHGTJs2zQCmdOnS\npnDhwgYw7du3N2fOnEl/Q4mJ6S6ycuVKc+XKlcyW7HF6P/aYETCr33vP7lKUFwI2mjQyNd2DuyKS\nR0TaAbOAvsBE4Cs3fQ6pHMQ4WvwAPXr04IsvviAyMpJ27dqxYsUK5s2bR+HChdNe+dw5GDUKLl9O\n83aKqUVFRXnEHbay2riJEymbLx/d3ngjw4PdKZVaegd3P8Xq5lkMjDLG7MyWqlSOYFKNz/Poo4/y\n6KOPOrfy99/D6NHQurXVx++j8hcvzqdLltCoUSOGDBnC5MmT7S5J5QDptfi7ABWxbrb+s+NG63qz\ndeW0DA+b8K9/wd69Ph36SRo0aMBzAwYwZcoUftA7dqkscMvgN8bkMsYUcDwKpngUMMYUzK4ilXdK\n2dXjtJgY+O036/ktbqzua0YPG0a1gAC6v/wyf+sIniqTnL1yVymXZSj4Bw+GunXhvH6hTCnwjjuY\ns2oVZ+Pj6dq16w23mVTKFRr8ym0yFPyvvw4ffQQF9QtlatXr1eOdd95hyZIlvJ0DLlRT9tHgV27j\nUvCfO2ddoVuiBLRr597CvNgzzzzDo0FBvDB1Kht+/tnucpSX0uBXbuN08F+7Bk2aQL9+7i/Ky4kI\nHy5dSomSJen0xBOcO3fO7pKUF9LgV26T+nTOm/L3h8ceg+bN3VtQDlEkIoI58+Zx6NAhejz6qEu3\ntVQKNPiVm6Xb4k9MhFy5YPhwePjh7CkqB6hXrx5vPPYYX/74I+P79rW7HOVlNPiV26Tb1bN7N9So\nYd1RS7ls8PTptI+IYNiUKSxfvtzucpQX0eBXbpNu8F++DHnzQtGi2VdUDiL58jF97VoqV65Mp44d\nOfj773aXpLyEBr9ym3SDv2ZN+OUX60welSH58+fn69mziTt7lnb16+sQzsopGvzKbW4a/KNHw9ix\n/9xgRWVKxYgIZnXpwpZTp+jevbte3KXSdctB2pTKjDSDPzHR6tvPgSNp2unBTz5hbLVqDB06lPKl\nSzPmrbfsLkl5MA1+5TZpBn+uXPDZZxAfr639LDZ48GD2LlnC6+PHU6FIEboPH253ScpDaVePcpvr\nzi8/fBjat4eTJ63AzyF3yvIkIsL7U6bQvHhxer38MitXrrS7JOWhNPiVWyW3+KOjYc0aOHHC3oJy\nuIDy5ZkfHU3FihV55JFH2LZ+vd0lKQ+kwa/c5rqunvvvh/37oVo1e4vyAYUKFWLxokUUuHqVFo0b\ns09P81SpaPArtzHGICtXwuLF1ozbbrO3IB8SUqYMy0aOJCF3bu67/36OHj1qd0nKg9gW/CLiJyJb\nROQ7u2pQ7mUSEpDz563+fZXtKg8fzg+rVnH69GmaN2lCzMmTdpekPISdLf5/A9E27l+5S2IiJCZi\nRJC2baFXL7sr8lk1a9ZkwbRp7N+7l/tr1ODUqVN2l6Q8gC3BLyKlgNbANDv2r9zIGOjTB3r2tLp6\n/PWMYbtFtW/P1507s/vUKZo1a0ZMTIzdJSmb2dXinwAMBW56iaGI9BKRjSKy8aR+RfUeIhAcDMHB\nOlywpxDhgVmz+HbBAvbs2UPT2rU5qWdX+bRsD34RaQOcMMZsutVyxpipxphaxphaQUFB2VSdyrCE\nBPjrL+v5K6/AmDEZu/WicpsWLVqwcMQI9u7fT9PISI4fP253ScomdrT46wMPicgBYC7QVERm2VCH\nykrPPmvdJD3pjlCOwNfg9yz3vfAC3w8Zwv6YGOrXr8++ffvsLknZINuD3xjzgjGmlDEmFOgErDDG\nPJHddags1qMHDBwIhQphjGHixImcOnWKQoUK2V2ZSilXLpqOG8eKFSs4d/Ys9WvUYLNe5OVz9Dx+\nlXExMfD559bze+6BAQOIiUe5720AABV4SURBVInhoYce4t///jetWrVi6NCh9tao0lS7dm1+evFF\nAmNjadysGT/++KPdJalsZGvwG2NWGWPa2FmDyoTXX4fu3eHYMQBWrlxJeHg4S5cuZeLEiSxYsIDb\nb7/d5iLVzVQeNIifN2wgtFw5WrZsybQPPrC7JJVNtMWvXJc03vtrr8HatcQVK8aIESNo1qwZBQsW\n5JdffqF///7av+8FSkZGsnbtWprWqMHTffsy8IknSEhIsLss5WYa/Mo1U6ZAs2Zw9Srky8eBokVp\n3LgxY8aMoXv37mzatImIiAi7q1QuKFy4MN/PmMGz5coxYfZs2rRpw7mkg/QqR9LgV64pXBgKFoS4\nOObMmUONGjXYtWsXc+bM4aOPPuI2HY/HK/lXrcp///iDKY4bt9euWJGdW7bYXZZyEw1+lb79+2H5\ncut5x45cmDmTrn378q9//YuwsDC2bdtGp06d7K1RZYlevXqx/NVXOXviBPfWqcOMGTPsLkm5gQa/\nSl/v3vD00xAXx4YNG7j7nnuYNWsWL730EqtXryY0NNTuClUWavzCC2xdtIja9erRtWtXej75pN7E\nPYfR4Fdpu3wZLl2ynn/4IVcXL+al0aOpX78+165dY9WqVYwaNQp/HYsnRwp+4AGWLVvGi3378tHM\nmdxboQLbtm2zuyyVRTT41Y2uXIHISHjuOQB+/ftvaj72GKNHj6ZTp05s27aNhg0b2lykcjd/f39e\ne/NNfnjwQWLi4oiMjGTs2LF61k8OoMGvbhQYCF26cOXBBxk2bBh16tThzJkzLFy4kJkzZ1KkSBG7\nK1TZJX9+WixYwI7du2nbti3Dhg0jqlIl9u/fb3dlKhM0+JXljz+gcWPYsQOApXffTfjAgYwdO5bu\n3buza9cu2rTRa+18VbFixZg3ezYzw8PZceQI4eHhTJgwQVv/XkqDX1kKFYJTpzi0ZQuPPfYYLVq0\nwBjD0qVLmTZtGoULF7a7QmUzyZ2bJ7ZuZcfu3TRu3JiBAwdSOySEzb/8YndpykUa/L7sxx+tM3aM\n4WqBArzZuTNVevdm0aJFjBkzhp07d9K8eXO7q1SeRITS5crx3Xff8fljj3H0+HEi69Xjueee48KF\nC3ZXp5ykwe/Ldu4kcflyPpsyhcqVK/PC8OG0aNGC6Ohohg8fTp48eeyuUHkoEaHD/PlER0fz9NNP\n8/bbb1OxRAk+fustEhNven8l5SE0+H1JbCwMGQKLFwOwokoVIgsUoHPv3hQuXJilS5fy1VdfUaZM\nGZsLVd6icMWKTJ48mV8+/5yysbE8NXQokY7xf5Tn0uD3Jf7+8P33bPz6a1q1akWzFi04GRPDjBkz\n2LRpk3brqAy7t0MH/vfnn3w2axYnTpygUaNGPFazJtGOkwWUZ9Hgz+kWLIDmza2rbrdupXVICJEf\nfsj69esZN24cv//+O126dCFXLv1TUJkjISE83rkze/bsYVSNGizZupWwiAi6du2qp396GP3fnhMl\nJkJcXPLk+qNHeeD++6lduzbrf/2VMWPGcODAAYYMGUJgYKCNhaqcKF++fLy0ZQv7d+xg0KBBzJs3\nj0oVKvB/bdty+PBhu8tTgBhj7K4hXbVq1TIbN260uwzvcPo01K1LYu/efFeuHG+//TarV6+maNGi\nDB48mL59+1KgQAG7q1Q+5K+1a3m9TRumxsaCCE88/jhDBg+mSvXqdpeW44nIJmNMrdTztcWfE5w/\nDz/9BEBsYCCT7ryTym+9Rdu2bdm/fz/jx4/nwIEDDBs2TENfZbsSDRvy3qlT7N27l2eeeYa5c+ZQ\nNTycR9q0Yb3e79cWGvw5Qe/eHGjdmheHDiUkJIQ+a9dSqEQJ5syZwx9//MFzzz1H/vz57a5S+TJ/\nf8qEhvLuu+9ycNYsXqpTh9U//0zdunVpFBbG/EmTiEvRPancS4PfG+3eDQ8/TPzhw3z77be0OnyY\nchcu8OZ//kODBg1Ys2YNGzZsoFOnTgQEBNhdrVLXCerQgVHr1nHo0CHeeestjkRH06FPH0JDQxk1\nahTHHPdwVu6jffze4sQJuHYNSpXi0Jo1TG/Viml583I0JoYSJUrQo0cPevbsSUhIiN2VKuWShP37\nWbxiBe9/+SU//PAD/sCjzZrRa/hwoqKi9IyzTLhZH78Gvze4epWzd9zBl3ffzUxg9erViAgtWrTg\nmWeeoU2bNjouvsoR9n35JZP69WN6bCxnz58npHhxujZpQteXX6Z8xYp2l+d1NPi9zdixXNuyhR/+\n9S9mzZrFgm++4WpcHBUrVqRLly488cQTeucrlWNdvnyZb7/9lo+HDmXZ4cMYoFGjRnTr0IF2nTtT\nSAcNdIoGv6fbuxe++IIrAwawdNkyvnzlFRbs3MnZuDiCgoLo1KkTXbp0oVatWoiI3dUqlT2uXuXI\n8uXM3L6djz/+mL1795I7Vy5atG5Nhw4deOihhyhYsKDdVXosDX5PY4w19n25clwEFj//PF9+8AHf\n58vHxdhYChcuzEMPPUT79u1p0aKFHqRVPs8Yw4YXX2Teli3M27mTI0eOkCdXLh6oVYsOAwbQqlUr\nChUqZHeZHkWD3xNcvWrd1rBQIfZ/8w2LHnmEReHhrPz9d65cuUJQ0aI83K4djz76KE2aNCF37tx2\nV6yUR0pMTGT9okXMe/pp5l+9yl9nzuDv70+j4sV5sEsXHnzqKcqXL293mbbT4LdLXBwEBHD1zBnW\nli7NorAwFp09y549ewCoWL48D7Rpw8MPP0zDhg3x8/OzuWClvEtiYiLr1q1j4XvvsfDzz9ntyLQq\n5cvTpmJF2jz7LHWbNfPJb80a/NklIQH8/EhISGBzvXqsvHKFFcWL89NPP3Hp0iXy5M5Nk6ZNadWq\nFQ888AAVKlSwu2Klco6LF9l//DjfLVrEwnffZfW+fcQB+fPnp3GNGtwXFsZ9vXtTLTzcJ46VeUzw\ni0hpYAYQDCQCU40x/73VOh4d/PHx4O9PYmIiO7t0YcXatayIiGDNmjWcO3cOgGrVqtGkSRNatmxJ\nkyZNyJcvn81FK+UDEhM5v3kzPx4+zPLly1k+cya/O+4SFhwcTLOqVWl277006tmTcuXK5cgPAk8K\n/uJAcWPMZhEpAGwCHjbG7L7ZOh4V/KdOQdGiXLx4kQ2DB/PzjBmsi4pi3fr1nDlzBoAKFSrQtGlT\nmjRpQlRUFMHBwTYXrZTi0iUOrV3L8r/+sj4I5s/nZHw8AMWLF6dB6dI0qFuXhl27Eh4eniO6XT0m\n+G8oQORb4D1jzLKbLWNb8MfHQ3Q0pmxZ/jxxgnXjxvHzlCmsq1aNbdHRybeYq1alCnXr16dBgwY0\nbdqU0qVLZ3+tSimXJJ46xe516/jpyBHWrl3LT3Pncsjxf7pAgQLULVqUeo0bE9mhA5GRkQQFBdlc\nses8MvhFJBRYA4QZY86neq0X0AsgJCSk5sGDB91f0PnzJH77LX+ULMmmEyfY/PXXbJo3j83583P2\n4kUA8ufOTZ06dagXFUXdunWpXbs2RYoUcX9tSin3uniRQ7/9xk+//85PK1aw9pNP2JWYSFJGlgkM\nJDIyksg2bbg3MpKaNWtSwMOvIfC44BeR/MBqYIwx5qtbLZvlLX5jQIQrMTH8NmAAO0qWZGt8PJt+\n/pkt69eT9AmUO3duwkuWpGbDhtxTvz61a9cmLCwsR3wFVEqlwxgunDrF5l27+HXJEn6dPJlfc+fm\nz7//BkCAiqVKUaNePWpUqkR4sWLUaNOGUmXLeszxAo8KfhEJAL4Dlhhj3k5v+UwF/4EDxF++zD4R\ndm7fzs5+/dhZrBg7ExPZu3dvcndNYGAgNWrU4J7QUGo2bUrNe++latWqei69Uuo6MTEx/PrNN/w6\ndSqbCxVi+x9/8Oeffya/XqRIEcLLlaNGnjyEd+xI9Tp1qFy5si1XGHtM8Iv1UfgpcNoYM8CZdTIa\n/P/973/5ZOhQouPjueoIeAEqBAcTVrcuYWFhhFWtSlh4OHfddZdPnuerlMq88/v2sWP2bLYVKsS2\n6Gi2L1vGjj//5FKKZUoEBFClXj0qh4VR5c47qRISQuX77qN4iRJu+4Zws+C3Y0jH+kAXYIeIbHXM\nG26MWZTVO0pMTCQ4PJz7qlYlrGlTwsLCqFKlip5OqZTKUgUrVKD+yy9TP8W8xDNn+CMmhl27dvHb\nF18QvWYN0bGxzJgxgwuO00oBChUqROWiRakYGEiFjh2pUKECdwUHU6F6dYq46YCy7Wf1OMOjTudU\nSqlMMMbw1+rV/LZmDdG33050dDTRCxey7++/OXzt2nXL3n777cydO5fmzZtnaF+e1OJXSimfJSKU\njIqiZFQUzZJmvv8+YA1HvX//fvbOmsW+/fvZV6QIZcqUyfoatMWvlFI5081a/HpPM6WU8jEa/Eop\n5WM0+JVSysdo8CullI/R4FdKKR+jwa+UUj5Gg18ppXyMBr9SSvkYr7iAS0ROAhkdkL8YEJOF5WQV\nrcs1WpdrtC7XeGpdkLnayhhjbhjwxyuCPzNEZGNaV67ZTetyjdblGq3LNZ5aF7inNu3qUUopH6PB\nr5RSPsYXgn+q3QXchNblGq3LNVqXazy1LnBDbTm+j18ppdT1fKHFr5RSKgUNfqWU8jFeF/wiUlpE\nVopItIjsEpF/O+ZHiMh6EdkqIhtF5F7H/KKO5S+KyHuptlVTRHaIyD4RmSiZuONxBupqLiKbHPvf\nJCJNPaSuex3ztorINhF5xBPqSrFeiOPfcrAn1CUioSJyOcV7NtkT6nK8Fi4i6xzL7xCRQLvrEpHO\nKd6rrSKSKCIRHlBXgIh86th/tIi8kGJbWVZXBmvLLSIfO2rYJiJRma7NGONVD6A4cI/jeQHgd6Aq\nsBR4wDG/FbDK8fw2oAHwf8B7qba1AagLCLA4af1squtuoITjeRhw1EPqygf4p1j3RIpp2+pKsd6X\nwHxgsIe8X6HAzptsy866/IHtQA3HdFHAz+66Uq1bHdjvIe/Xv4C5Kf4PHABCs7quDNbWF/jY8fwO\nYBOQKzO1eV2L3xhzzBiz2fH8AhANlAQMUNCxWCHgL8cyl4wxPwFXUm5HRIoDBY0x64z1Ds4AHs7G\nurYYY/5yzN8FBIpIHg+oK9YYE++YH+hYzvb3y1HDw8B+rPcraZ7tdaXFA+q6H9hujNnmWOeUMSbB\nA+pK6XFgDnjE+2WA20TEH8gLXAPOZ3VdGaytKvCjY/kTwFmgVqZqy8wnl90PrNbWIcebVcXx/DBw\nFOtS5ZTLdiNFix+oBSxPMd0Q+C6763Is/1hSLZ5QF1AbK1wvAo94Ql1Y39zWAfmBV3C0+D2grlDg\nErAFWA009JC6BgAzgSXAZmCoJ9SVavk/gDBPqAsIAOYCJx3/nr3cXZcLtfXC+pbrD5TFCv5HM1Ob\n17X4k4hIfqyv/QOMMeeB3sBAY0xpYCDwUXqbSGNeps9tdbUuEakGjAWe8ZS6jDG/GGOqAZHAC46+\nYbvrGgW8Y4y5mHoTNtd1DAgxxtwNDAI+E5GCHlCXP1YXZ2fHz0dEpJkH1JW0fG0g1hizM2mWzXXd\nCyQAJbDC9TkRKeeuulysbTpwBNgITAB+BuIzVVtWfXJl5wPr03kJMCjFvHP8c12CAOdTrdON61v8\nxYHfUkw/DkzJzrqAUlj9e/U9qa5U667EalnYWhewFqvf9QBWi+c00M/uutJYd5WHvF+dgE9SLDcS\nGGJ3XSlefwcY7il/98D7QJcUy00HOrijriz4G/sZq/snw7V5XYvfcdT6IyDaGPN2ipf+Aho7njcF\n9t5qO8aYY8AFEanj2OaTwLfZVZeIFAa+B14wxvzPg+oq6+jnRETKAJWAA3bXZYxpaIwJNcaEYrV6\nXjfGvGd3XSISJCJ+juflgLuwDljaWhdWqISLSD7Hv2djYLcH1IWI5ALaY3WtAPb/3WN1sTQVy21A\nHaxQzdK6MlKb49/wNsfz5kC8MSZz/5aZ/eTK7gfW11aDdcbCVsejlWP+JmAb8AtQM8U6B7BaiBex\nvjJVdcyvBezE6mt8D8enbXbUBYzA6kvcmuJxhwfU1QWrf38rVt/wwym2ZVtdqdZ9hevP6rHz/XrU\n8X5tc7xfD3pCXY51nnDUthMY50F1RQHr09iWnf+O+bH60XcBu4Eh7qgrg7WFAnuwDgIv5/rjcRmq\nTYdsUEopH+N1XT1KKaUyR4NfKaV8jAa/Ukr5GA1+pZTyMRr8SinlYzT4lccTa4TVpNEcj4vIUcfz\niyLyQTbXMltE9ojIThGZLiIBqV6PFJEEEXksxbwDjhEUt4rIxhTzbxeRZSKy1/GzSIrXXnCMuLhH\nRFqkmJ/maIyOcZ4+d8z/RURC3fk+KO+mwa88nrEGGIswxkQAk7GGbogwxuQ3xvTJ5nJmA5WxRpbM\nC/RMesFxIddYrIunUmviqLlWinnDgB+NMXdhDcI1zLGdqlhX3lYDWgIfJF0kBkzCGrvlLsejpWN+\nD+CMMaYC1lWxYzP/q6qcSoNfeS0RiRKR7xzPXxFrPPWljhZ2OxEZ52gd/5DUMne0mFeLdQ+EJWKN\ncOg0Y8wi44A1JG6pFC/3xxp75YSTm2sLfOp4/in/jKzYFmuI4KvGmD+BfcC9cuvRGFNu6wugWdK3\nAaVS0+BXOUl5oDVWCM4CVhpjqgOXgdaO8H8XeMwYUxNrPJYxGdmRY1tdgB8c0yWBR7C+kaRmgKWO\nD5teKebfaazL7nH8vMMxvyTWCI1JjjjmlXQ8Tz3/unWMNaz2Oawx+JW6gb/dBSiVhRYbY+JEZAfg\nhyOUgR1Yl71XwrrpzTJHY9gPa3TNjPgAWGOMWeuYngA8b6wx71MvW98Y85eI3OHY92/GmDW32PbN\nRl281WiMbhtFUuU8GvwqJ7kKYIxJFJE48894JIlYf+sC7DLG1L3ZBhx96ZsckwuMMS+lsczLQBD/\nDKUN1pgpcx2hXwxoJSLxxphvjOOGO8aYEyLyNdYQwGuAv0WkuDHmmKMbJ6mL6AhQOsW2S2EN4HWE\n67uWkuanXOeIY1C2QljjUyl1A+3qUb5kDxAkInUh+T6r1VIuYIxJSDqQfJPQ7wm0AB43xiSmWK+s\n+Wf00C+APsaYb0TkNhEp4Fj3Nqw7YyWNQb8A6Op43pV/RlZcAHRynKlTFusg7gZz69EYU27rMWBF\nig8+pa6jLX7lM4wx1xynWU4UkUJYf/8TSHErRydMBg4C6xyt+6+MMa/eYvk7ga8dy/oDnxljkrqg\n3gTmiUgPrGGB2zvq3CUi87BGiYwH+hpjEhzr9AY+wTqjaLHjAdYwvzNFZB9WS7+TC7+T8jE6OqdS\nSvkY7epRSikfo8GvlFI+RoNfKaV8jAa/Ukr5GA1+pZTyMRr8SinlYzT4lVLKx/w/6fpY1DbNiuYA\nAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Plot those models:\n", "my_pspl_model.plot_magnification(\n", " t_range=[2452810, 2452890], subtract_2450000=True, color='red', \n", " linestyle=':', label='PSPL')\n", "my_1S2L_model.plot_magnification(\n", " t_range=[2452810, 2452890], subtract_2450000=True, color='black', \n", " label='1S2L')\n", "plt.legend(loc='best')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Suppose you also had some data you want to import:\n", "\n", "OGLE_data = MulensModel.MulensData(\n", " file_name=os.path.join(\n", " MulensModel.DATA_PATH, \"photometry_files\", \"OB03235\", \"OB03235_OGLE.tbl.txt\"),\n", " comments=['\\\\','|'], plot_properties={'label': 'OGLE', 'color': 'orange'})\n", "\n", "MOA_data = MulensModel.MulensData(\n", " file_name=os.path.join(\n", " MulensModel.DATA_PATH, \"photometry_files\", \"OB03235\", \"OB03235_MOA.tbl.txt\"),\n", " phot_fmt='flux', comments=['\\\\','|'], plot_properties={'label': 'MOA', 'color': 'cyan'})" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/jyee/MulensModel/source/MulensModel/utils.py:149: UserWarning: Flux to magnitude conversion approached negative flux\n UserWarning)\n/Users/jyee/MulensModel/source/MulensModel/utils.py:150: RuntimeWarning: invalid value encountered in log10\n mag = zeropoint - 2.5 * np.log10(flux)\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOy9eVxU5734/35mANlxQWAIuMAwIAqa\nVZsmCopJs5k2baK5XRJrXJKYWJv87m2b2+/NvbfNvW3SW7MaE26StrdVm6RJzB4XRmPSpMYVERgG\nXEAGcGdAEJzz/P44ZyYDDItsKj7v1+u8ZuY8Z/k855x5Puf5fD7P5xFSShQKhUKh6A9M51sAhUKh\nUAwdlFJRKBQKRb+hlIpCoVAo+g2lVBQKhULRbyilolAoFIp+QykVhUKhUPQbSqkoFAqFot9QSkWh\nUCgU/YZSKhcoQoh7hRCFQojTQogaIcRKIcRwv/I0IcQaIcQRIUS9EKJMCPGsECLJKM8RQlR1cuzX\nhBAtQogGv2V3F7I8ZRzfLYQoEUL8qF35FCHEdkPW7UKIKe3Klxt1OCWEeEUIMcyv7P+EEC6jDg4h\nxH1dyJEjhND8ZD4shPj3dttIIYS1k/0tQoiXhRDVxv4VxrXI6E3dhRCxQojPhBDHhBAnhRB/F0J8\n06/8XiGEp911zml3/ETvfRJCHDDuS2y7bXYZ9Rrnt+5aIcQmQ65TQoh3hRCZAeS/19j3rs7qaGxn\nE0K8YzxPx4UQHwsh0ttt09V9HCmEeEsI0SiEOCiE+Ce/skwhxFdCiBPGsiGQrH7b+z+fbuOZmtGu\nTlu72H+2EKLA2PeYcf3+RQgR2sU+0pDde5/y/crmCSFKjXrXCSH+IISI7up6XtJIKdVygS3AI0At\n8C0gGBgHfABsA0IAK3Ac+B8gydgnDvgJMM/4nQNUdXL814BfnYM8/w5koL+ETAVOANcaZSHAQWA5\nMAx42PgdYpTfaNRlIjACsAP/7XfsicAw43sGUANc2YkcbeoEjAeqgG/7rZOANcC+o4D9wJ+BVEAA\nw4H5wEO9rHsokG6UCeDbxn0JMsrvBbZ2c23vA/KN7weAUn95gCxjnQTGGeu+ATQAy4AoYCTwK0O2\nlHbHLwCOAe93I8c1wALjWMHAfwIlfuXd3cfVwFogErgOOAVMNMqGoz/DAjAbz8ieLmR5DeP5NK7t\nfUYdzN1dV+BO49wLgZHGunTgWSCti3MGfG6MsmQg1vgeaTxDz5zvduJCXc67AGppd0Mg2mgw7mq3\nPhKoA34M/B/wbjfHyaGflEqA/dcBjxjfbwAOA8Kv/BDwLeP7X4An/MpmATWdHDcdcLWve1d1Av4K\n/MLvd2dK5VfAbsDUx/vjq3u79SbgNuP8cca6Ths/v/3+BtxhfD8A/Cuwza/8KeAx2iqVT4EXAhzr\nQ+CPfr/HAhrwXeAsEH8O9RxpnHNUd/cRiABaAJtf+Z/wUzp+64OAB4HTXZy7zfMJhBuyJHZ1XdGV\nVmWg+9OD+naqVNptFwn8EfigL8/RUF6U+evC41r0N+C/+a+UUjagNxqzgTzgzcEXDYQQYcDVQJGx\naiL6W6d/Erk9xnpvub9pbTcQL4QY5XfMF4QQp4ESdKXyQQ9lSQO+CXzRg83zgLeklFpPjt3J+drX\n3bt+D9CMrnDypZR1fsWXCyGOGqa9Xwohgvz2CwamA+v9tv8CiBZCTBBCmIG56C8R3n3C0Z+R1wOI\n+Ff058PLj4CvpJRvAsXA98+hutPRlcYx43dX99EGeKSUjnblE/1+I4Q4iX6dngWe6IkQxjX4EXov\ns7abzdOBJHr/39himPf+5m9qNOS4TghxCnCjK+kVvTzHkEcplQuPWOColPJsgDKXUR6LbiYCQAix\n1LDpNwghXu7heR419vEuf+jhfi+iNxgfG78j0c0N/pxCN8sEKvd+95YjpXzA+H09ujI908X5Ew15\n6wEH8CXQqX3dj/bXbI5xHLcQ4pMe7A8d6+6VPxu9h/lP7WTZAkxCN01+F7gb+P/8yqcDu6WU7nbn\n+RN6QzobXdEe9isbif6/dQWQz/t8ePkReg8D4/OeLmtnIHS/3PPAT/1Wd3Ufu3sGAJBSDgdigKXA\nzm7EeNRQQo3oDfgvpZSebvbx1t3/Pq8x7vNpIcQPu9h3BrqJLgOoBt7zfwGQUm6VUsagK60n0XuV\nigAopXLhcRSI9X+g/bAY5ceM7wBIKZ8z/rAr0O3hPeEpKeVwv+UeACHEi37Oyl/47yCEeBK9kbzL\nr2fSgN6g+hON/kYXqNz7vU1DKqX0SCm3ov9p7zfO96GfLN637GpD3mh0W30T0BOF2P6arTOu2XJ0\nv1Bv6u4vf7OUcjXwMyHEZGNdhZRyv5RSk1IWAv8BfM9vt5sJ3Cv7E7qCuhfd1OLPCXSTloWOeJ8P\nhB4wMB5YY5T9BcgSRhCFaBs8MMavnqOBT9DNa6v9jt3VfezuGfAhpWxEV85/FELECSHG+Mvit+lT\nxv0JA64CnhRC3BSgzv54e1X+93mecZwd6P6cgHWXUm6RUrZIKU+i+6rGAxMCyH8Y+Iivr6uiHUqp\nXHj8Hf1N/Q7/lUKICOAmYKOx3NFx174jpVwipYw0Fp+JQuhRVjcBN0gp6/12KQKyhRDCb102X5uI\nioDJfmWTgVo/s0p7gtAd6Ugpb/KT5c8BZD2F3lje1oOqbQS+LYTo9JnvRd0DEQykdHYKdLu/l5uB\n9wPIcRDd3HMzHc2gjejPyJ0Bjn8Xej1B75UIYJcQoga9Rwd67wW/ekZKKQ8Z9RyBrlDWSSl/3e7Y\nXd1HBxBkmCP9y9uYCf0woftJLpNSHvKXJcC1kFLKvcBnwC2dHM+Lt1fX5X8jUN0DbUbbe+WP7xlV\nBOB8O3XU0nEB/pnA0V870COsbMBJ9Oivy4x9YtFtya8Zv3PQI6NC2y2Cc4/++jlQBlgClHmjv5YZ\nsi2lbfTXt9DNEZnoUUObMBy46GaheejmEzN6hFEjcHsncuTQNvorEv2t/ku/ddI4l3+dzcb1OUTb\n6K8o4DnA3su6T0OPdApBf6P+F/S3c69D+SYM5zi6WWUv8G/G7/FARbvjHQDyjO+pwFXG9yDaOuqv\nM67Tw0YdRqAHIpwE0ow6n0SP5krwWx40nqugAHWJBv4BPNfJdej0Phrla9AjwCLQ/Vz+0V+zgcuN\n+xANPINuYgrt5Fyv0dZRn4Fu2rvf+H0vupJp82wbZXOBevTorxHGfU4z7uG9nZxvIjDFkC8Svcdf\nCgQb5d8HxhjHGgtsBv52vtuJC3U57wKopZMbozcIe9HNO7XAKmCEX3kGumP2qNGQlaI7QJON8hyj\nIWq/WI0/bQu62cK7HO1CFonee/Lf3j/i6nJguyHrDuDydvv/1KhDPfAqX4cQjzb+oCeNskJgYRdy\n5KCbfrwyHEN/07e2k7X9cp9Rlgj8r9FANQDl6KazCb2pO7odfrdx/Y8bdZnut+9TRr0bgQp085e3\noVpKuwYcP6XSbn0bpWKsuw49rLfBuHbvA5OMsnlGHYPbHSfUeF5uDXCOe4xzNLar65ju7qNRNhJ4\n29j/EPBPfmV3ovciGoAj6C9I2V1c89f4+vn0Hu8JjMg9dKUS6D57Q7m/ZdwL7zOyE92XFdHJ+Wai\n/38a0SMs38Yv/Bj4NfoLWqPx+RJGVJxaOi7CuGgKhWIQEUJ8gK5UehTpplBcLAyYT8UYcVsnhNjb\nbv1DxujUIiHEbzvZ94DQR5PvEkJ8NVAyKhTnETv6wESFYkgxYD0VIcR09O7nH6WUk4x1uegDuW6R\nUp4RQsTJtjH93n0PoNuTjw6IcAqFQqEYEAaspyKl3IJuZ/bnfnTn3hljmw4KRaFQKBQXL4HGQgwk\nNuB6IcSv0UfWPiql3BZgOwl8IoSQwCop5UudHVAIsQhYBBAREXFlRkanuQEVCoVC0Y7t27cflVKO\n7q/jDbZSCUIP85uGnu7ir0KIFNnRBvdNKWW1ECIOWC+EKDF6Ph0wFM5LAFdddZX86ivlglEoFIqe\nIoQ42J/HG+zBj1Xo8d1SSvkP9PDQ2PYbSSmrjc864C30DKoKhUKhuMAZbKXyNnpMOEIIG/qgsTbO\neCFEhBAiyvsdPQvuXhQKhUJxwTOQIcWr0dNJpAshqoQQC4BXgBQjzHgNcI+UUgp9oiJvvH48sFXo\nk0b9A30eiI8GSk6FQqFQ9B8D5lORUt7dSdEPAmxbjZ7nCCllBW1zDCkUCkW/0draSlVVFc3Nzedb\nlEElNDSUpKQkgoN7mnO2dwy2o16hUCjOK1VVVURFRTFu3Dja5kEdukgpOXbsGFVVVYwfP35Az6Wy\nFCsUikuK5uZmRo0adckoFAAhBKNGjRqU3plSKgqF4pLjUlIoXgarzkqpKBQKRTfkGIuie5RSUSgU\nikFGCMEPf/j17MZnz55l9OjR3Hrrrb51b7/9NtnZ2WRkZJCVlcXbb7/d5hhHjhwhODiYVatWDZrc\nPUEpFYVCoegCTdMY5XAwdvNmHA4Hmqb1+ZgRERHs3buXpqYmANavX89ll13mK9+9ezePPvoo77zz\nDiUlJaxbt45HH32UPXv2+LZ5/fXXmTZtGqtXr+5w/POJUioKhULRCZqm8eratcTb7VS3trLKbufV\ntWv7RbHcdNNNvP++Ppv06tWrufvur0dhPPXUU/ziF7/wRWqNHz+en//85zz55JO+bVavXs3vfvc7\nqqqqOHz4cJ/l6S+UUlEoFIpOcDqd7HO7WbVgARvy8lixYAHF9fU4nc4+H3vevHmsWbOG5uZm9uzZ\nw9SpU31lRUVFXHnllW22v+qqqygqKgKgsrKSmpoarrnmGu666y7Wrl3bZ3n6C6VUFAqFohNcLheF\nKSloZjMAmtnMntRUampq+nzs7OxsDhw4wOrVq7n55pvblEkpO0Rr+a9bs2YNd911F6ArpwvJBKaU\nikKhUHSCxWIhq6ICk8cDgMnjIbu8nISEhH45/pw5c3j00UfbmL4AJk6cSPuM6zt27CAzMxPQTV+v\nvfYa48aNY86cOezevZuysrJ+kamvqBH1CoVC0QlWq5XM7dtZnJ9PWWoqk8vLmRAdjdVq7Zfj//jH\nPyYmJoasrCzsdrtv/aOPPsqdd97JzJkzGTduHAcOHOCJJ57gjTfeoLS0lMbGxjZ+lH/7t39jzZo1\n/PKXv+wXufqC6qkoFApFJ5hMJubPnUttbi6JISEsys1l/ty5mEz903QmJSWxbNmyDuunTJnCb37z\nG2677TYyMjK47bbb+O1vf8uUKVNYvXo13/nOd9ps/93vfveCMYEN2Bz15wM1SZdCoeiO4uJiJkyY\ncE775Bif9v4WZpAJVHchxHYp5VX9dQ5l/lIoFIpusJ9vAS4ilPlLoVAoFP2GUioKhUKh6DeUUlEo\nFApFv6GUikKhUCj6DaVUFEOaHFTKckU/sCFHXxTdopSKYsgyENllFYr+oqqqittvv520tDRSU1NZ\ntmwZLS0tAPzjH/8gJyeHtLQ0rrjiCm655RYKCwsBePzxx3nqqac6HM9sNjNlyhTf8t///d+DWh8v\nKqRYMSTxZZd1uylLSWGV3U7m9u39OnBNcWmgaRrO46NwNUZicTiwWq19foaklNxxxx3cf//9vPPO\nO3g8HhYtWsRjjz3Go48+yl133cVf/vIXrr32WgC2bt1KeXk5WVlZnR4zLCyMXbt29Umu/kApFcWQ\nxD+7rGY2s8njYXl+Pk6nE5vNdr7FU1wkaJrG2j+/irs6npSwMuzrVrE9PpO535/fJ8WyadMmQkND\nmT9/PqD3Mn7/+9/7Ut3fc889PoUCcN111/WtIoOIemVTDEkGMrus4tLB6XTirt3HgqRV5MVuYIFl\nBfW1xX1OfR8otX10dDRjxoyhvLycK6644pyP2dTU1Mb8db7S4aueimJIYrFYyLLb2ejxoJnNX2eX\nzc0936IpLiJcLhcpIYWYhe6PMwuN1JA91NTU9KnHGyi1vXd9+9RZU6dOpb6+nhtuuIGnn36602Ne\nKOYv1VNRDEmsVitjhGDhihXMWr+e5fn5/ZpdVnFpYLFYqGjJwiP1ptIjTZS3ZPc59X2g1Pb19fVU\nVlZitVrZsWOHb/2XX37Jf/7nf3Lq1Kk+nXOwUEpFMSQxmUxMHD+egpdeIunPf+737LKKSwOr1UpU\nfCb5VYvZcDSPfNdyouMn9PnlZNasWZw+fZo//vGPAHg8Hh555BHuvfdeHnnkEV577TU+//xz3/an\nT5/u0/kGE/UPUwxZNE3D4XBQXV2NzWZTCkVxzphMJuZ+fz6542sJiUokd86iPjvpAYQQvPXWW7z+\n+uukpaVhs9kIDQ3liSeeICEhgbVr1/Lzn/8cq9XKtddeyxtvvMHSpUt9+//qV78iKSnJt0BHn8rP\nfvazPsnYW5RPRTFkOXv2LABBQeoxV/Qek8mEbeQxbCOPQT9GDiYnJ/Puu+8GLJs2bRqbN28OWPb4\n44/z+OOPd1jvMWanPN+of5tiyOJVKmYjAkyh6DV59vMtwUWDsgcohizeNzelVBSKwUMpFcWQRZm/\nFJ0xlGa87SmDVWelVBRDFtVTUQQiNDSUY8eOXVKKRUrJsWPHCA0NHfBzqVc4xZDFq1RUT0XhT1JS\nElVVVRw5cuR8izKohIaG+iLFBhL1b1MMWZT5SxGI4OBgX44tRf+jzF+KIYtSKgrF4DNgSkUI8YoQ\nok4Isddv3VohxC5jOSCECJioRgjxLSFEqRDCKYQ4PyN4FBc9yqeiUAw+A9lTeQ34lv8KKeVcKeUU\nKeUU4E3gb+13EkKYgeeBm4BM4G4hROYAyqkYoiilolAMPgOmVKSUW4DjgcqEnp7zLmB1gOJrAKeU\nskJK2QKsAW4fKDkVQxdl/lIoBp/z5VO5HqiVUpYFKLsMqPT7XWWsC4gQYpEQ4ishxFeXWjSHomvU\niHqFYvA5X0rlbgL3UgA6TjIAnQaUSylfklJeJaW8avTo0f0inGJooEKKFYrBZ9D/bUKIIOAO4MpO\nNqkCkv1+JwHVAy2XYuiheioKxeBzPnoqeUCJlLKqk/JtQJoQYrwQIgSYB6wbNOkUQwblqFcoBp+B\nDCleDfwdSBdCVAkhFhhF82hn+hJCJAohPgCQUp4FlgIfA8XAX6WURQMlp2Loohz1CsXgM2D/Ninl\n3Z2svzfAumrgZr/fHwAfDJRsikuD1tZWQB9BrVAoBgc1ol4xZPEqFTXjo0IxeKh/m2LI4lUqmqad\nZ0kUiksHpVQUQ5aWlhbg0pw7Q6E4Xyilohiy/M3oqSilolAMHiosRjEk0TSN5KAg0qZPR9M0NE1T\nvhWFYhBQ/zLFkEPTNF5du5bL09MJfuwxjsfG8uratcq3olAMAkqpKIYcTqeTfW43Lz/4IBtvuIEX\nHniA4vp6nE7n+RZNoRjyKKWiGHK4XC4KU1LQjJH0mtnMntRUampqzrNkCsXQRykVxZDDYrGQVVGB\nyUjTYvJ4yC4vJyEh4TxLplAMfZSjXjHksFqtZG7fzsJnn8WZmcmE0lImxMZitVrPt2gKxZBH9VQU\nQw6TycT8uXP5oqCA1l//mh0uF/PnzlXRXwrFIKD+ZYohiclkwlRZSU1NDSOFwOl0qugvhWIQUEpF\nMWj88z//M2lpaYNyLk3TuPyb3yRn0SKacnJYZbcPWlhxjrEoFJciyqeiGDSefPLJQTuX0+kkePx4\nXl62DM1spsDjYXl+Pk6nE5vNNmhyKBSXGqqnohiSuFwunJMmqbBihWKQUUpFMSSJi4sjtbBw0MOK\nczSNWoeDYZs343A4lB9HccmhzF+KIUlcXBzu4mIWrliBMzubKRUVTIiOHtCwYk3TyFy7lnS3m/KU\nFFbZ7WRu364izxSXFEqpKAYdKSVCiAE9R0NDA2teeYW0tDRs3/42ixYswGq1Dmjj7nQ6we0mf8EC\n5cdRXLKo1yfFoOOdO34gOXbsGFJKHA4H406fxmazDXhvweVyUdbD9DA5qAgxxdBEKRXFoOOdPGsg\nOXLkiO/7QM+nkmMsFouFNJUeRnGJo8xfikGnpaWFiIiIAT2Hv1IZLGe51WqF7dtZmJ+PMzWVKeXl\nnI6OZpHVit1vuxxgFzBlUKRSKAYX1VNRDDqD2lMxm3l7kGZ+NJlMrJ47l4LcXEwhISzKzWXf3LnQ\n3uymaYx3OBirIsQUQxDVU1EMOoOhVKqrqwkNDcUTHQ3dKJUc49Pei/NomsYop5NIl4uS+HjigITa\nWiotFlJSUnxlDovFF3mWuXYtmooQUwxRlFJRDDoDqVRyjM+EykqSk5MJHjGC8REROByOfo/+8s4w\nGe92UzZ+PK++9x55oaEciYvjqm3b+BUQHx1NmZ/y+Obll4PbzctGhNgmj4fF+fl84HTypooQUwwB\nlFJRDDqtra0Dfo7Kykpm3347Zy0WyrOyBqRH4J1hctWCBVjLy0kpL8cdFUVCTQ1N4eGc0TRW+SmP\n5fn5jNi3D0e7CLGy1FQSa2pAKRXFEED1txWDhrcxH6ieiqZpjDJ8FcOHD0ckJ5O/bBkbZ89mxYIF\n/T6lsP8Mk5bDh2kNCmLcwYMANIWGUpqe3iG8GCC1XYRYWnk5DSpCTDFEUD0VxaARFBRES0vLgCiV\nNqaolBSyr7uOPRkZWMvLSXC5qLFY2JOSQk1NTZuBiP4+Ea/fY6ah/OzdnNNisZBlt7OptZWrC79i\nTEMlGRHF7GnJ5mDrOFLKy0lPTia+tpba+HiynU4ypk/n048/5v6VKylNTyejtJTQ1laOpaT0+zVR\nKM4HSqkoBo3g4OABUyr+pijNbKYqMZHb3n+flAMHKE9JIcduZ0R9PXG33urbp70i8prICBSxFQDv\nDJNLX3iB5NNVPDD2ecxCY6bcxO+rf8LRs7HMWbeOA2PHklVYSExrK1JKmkNDKcjNJb6ujvV5eeTY\n7ZyqqFDmL8WQQCkVxaARHBwMDIz5y98UBaCZTJwOC2uTMmXpiy+22ae9Itrk8bD85ZcZt3kz0mTi\nuxYLm6xWJptMbXotmqbhfONOyo6E0ei5hiNVVdyStA+z0EODBZKw1jNchov00BKKDk7icGgixzDz\n3scf48zKojQjg9KMDACSqqqI6QefSo7xae9iG4VioFE+FcWg0RelkkPXaU0sFgtZfr4KS01NB59G\ncXo6dXV1vn3aKyIpBE0eD+ElJVS3tJD80Uf86IUXGFVS4htLomkaq174PR+WjCa4tYbW40WMbG6g\nyG3DI/W/k6PRhgQeGPs8N8Su5+HLniaxsYrw06dpbGjg8p07CTpzBtB9KrayMmqVT0UxRFA9FcWg\nMZA9Fa8panF+Po7x48navRstPJxNM2eimc1fp0zJzfXt4/WJbPR40Mxm0ktKEFKy6r77uHfNHxjb\nvJ9Jw/ZS8lYZr4xMZcbMm3j88ce5Ii2Sn6TlYxYas0YX8EzEck7JGJ6rXEpm2D62N1zFlZFf+Xou\nUoLnbAiJZhfpMaWUNqbz+JP/QcHV08kuKiK4tRU0DU3T1FgVxUWPUiqKQSMoSH/cBiKk2GQyMX/u\nXN4rLKR12TKipk9nzGWXsfCZZ3BOmhQw9b1XES3Mz8eZksK0bdv48uqrSS0vJ/twIQ8mP+fzkTx/\naClPPfssR44cIfuaOp/CMAuNScP28HHQt2k+aaauIY6zZjPFpzOZKTdhFhpbjs8gzHya+8eu1I8X\nu4kXDjzAN3d8wUhxnMywfVz2VjUvRiRx35JlhISEdFvfHOPT3u9XUqHoG+q1SDFoeJVKoJ5KDn3P\n2jvTZOKTgwfZsnkz06dP58fz5vHZO+8Qsm4di3JzO4xR8SqiL3NzGeV20xIRQcr+/UwqKmJCWFsf\nSXxwDaPi4/nBD37A/tZsn6nLI01UtGQzdfRogs6e5WRMDGdCQzlijuWFgw+y/kgen5+6lvSI0jaK\nKG5YLbHyKEuTn2P26A08OOY5zjZW8l8rVgxKFmeFYqBQSkUxaPTW/OU//iRQrixN03AY5ZPWrycm\nJobp06djMpkor67mq6NHO019bzKZ2G+zcXLUKHZlZuKOisLqdFJ8OhOPNKFJwZrqeRxpiSM3ajMN\n1TtolMPJr1rMhqN55LuWEx0/gdjYWI6PGEFEUxNN4eH8/fJpVJkS+aTlBtzBkZQ2prdRRAebxpHR\nTtFMCC+m5exZtmzZ0pvLq1BcEAyY+UsI8QpwK1AnpZxkrFsLpBubDAdOSik7JGsVQhwA3IAHOCul\nvGqg5FQMHr1RKp2F/Xp7HZqmsfbPr+Ku3s4PwsrYE5vKsUX3+npFCNFl7q8c4BRwKD6eGzZu5PlF\ni7j+00/J+WIzzx9aSnxwDcdbR3H/2BcMU9hGXq5eju2aOzGbzeQmJGC1WnE6nQz/8ktOREf7Is6q\nkpPJW7+eTTNnMm/dWlYefID0iBJKGzNoCI6k+HQGM+VGzELDI03sac3mdFgYDoeDmTNndnlN2o+t\nUb4YxYVCj55EofMDIcT/M36PEUJc081urwHf8l8hpZwrpZxiKJI3gb91sX+usa1SKEOE3igV/7Df\nDXl5HUbGO51O3LX7WJC0itmxG/iJ9WUuGyl95cJk8imV4cYSCAGENjVx3yuvYPZ4aDCHU6PFsedM\nNumRbXsU1mF7MJvNTJ8+3dcDslqtxAQHU5aW5osmc9hsSCG48eOPqYmKo641FntDDn+b9h2e+OnP\nOGIezbNVD/HJ0dmscC2nOD6TvRMnEhoa2un18ClZu53q1lZW2e28unZtnzMd56AmDVP0Dz19vXkB\n+AZwt/HbDTzf1Q5Syi3A8UBlQp9L9i5gdQ/PrxgCdOVTaU+OsXQYf9JuNkWXy0VKSGGbRj89vPjr\n2RaFQHbT4JqB5Npadk6eTEFuLq3DhrH96qs5FR1Nqc1GceOENqar8pbsgBNvhUVFkVFa2iYFS8Tp\n0wR5PGAy8fdp0/CYzUzct4/czZtpFKEcCh7DE9m/4KU5i3ht3o/IKC3l6quv7lTW7pSsQnG+6alS\nmSqlfBBoBpBSngC6D1HpnOuBWillWSflEvhECLFdCLGoqwMJIRYJIb4SQnzlPzGT4sLjXJSKl/bj\nT9rPpmixWKhoyWrT6DuaJjJ69Gg0TSM1MZGrRo/G4XBQr2mc8jt2DvpkWQB1FgtZ+/fjTE1ly/Tp\nSACTiVHHj1MVkcTzh5ay/uhy09gAACAASURBVEgeq6oeJjp+QpsoMtAb+wNScmzkSBa99BJ569fz\n8LPP4gkOpj4mhlWLF7Phxhv53fLlhDY1kVpezjfGjSPY42FCcTFjDh5k4csvE+F2s/uTZzqdZ6U7\nJatQnG966lNpFUKY0Rt7hBCjgb70t++m617KN6WU1UKIOGC9EKLE6Pl0QEr5EvASwFVXXTU4szEp\neoXX/HUuIcX+40/KUlOZXF7eJjTYarXyVdwEXjj4IBkRxRQ1T+SkNpy/fvABl+3cyTduuYXySZNY\nVVDAgs2bKUtL47uJiaxNSWFURQVTDx8m+tgxot1uhh075svJddVXX+GOiiL/vvuQQmBzOJizbh0T\n4kYz9/v3dPBhuFwuClNT2TBzJmlOJwk1NZwZNoymsDAOjRnjUwJnhw1j15QptAwbxt1SsnP4cA6O\nHUuCy0Ww2UOCuY7RpgoK3lnFV/ETuOKa66itrcVi+E7aj63xKtm4GTMY5XAoP4vivNNTpfIM8BYQ\nJ4T4NfA94F97c0IhRBBwB3BlZ9tIKauNzzohxFvANYAKibnI6U1PxRv2+4HTSWJNDYtyc9s0mCaT\niZFxSXx1tJ43J9+BKzGR8vHjefjZZzlx6hT5Dz+MFIK4I0cYeeIEQWfPEl9QwH+9/z7xYWGkuN00\nhYVRN3o0YZGRbMjLI76ujpMxMTj9/COlGRnsrKxkYnV1wMba19jn5uKw2XCmpjJt505C3G5SKyoo\n8FMCaU4nf581iziPh/SNG9loOOXnbFvHQ0nPYhYarVoBKw48Qn31TqzhDuzbstken8mdd9/z9dga\nY8rijKgotu7YQXxDQ8BgBoViMOmRUpFS/lkIsR2Yhe7T/LaUsriX58wDSqSUVYEKhRARgElK6Ta+\n3wD8Ry/PpbgAOdeQYpPJxDGbjY02GwfpOOBv27ZtFE6bxuacHN+6k8OHU5WcjGY2Y3M4iGpo4MVF\ni0jdv5+YU6cwu1yUWq1MKC1l5ZIlXLd1K8dHjfLl5DJ5PEwsKmozIj+jtJQxEycGlLHNQMrUVGzl\n5VhjY9l7+jTRJ0+yZNUqHGlpZDoctI4YQYzVistuJ6ypiYX5+ZwJCSE7eI/PN7T/dCoRZjcLk1cZ\n0WGbyHctp6Kigvlz5/Jboze0KDcXTdN4ecuWDnO3OJ3ONhmZFYrBoMvXGCHESO8C1KGbrP4C1Brr\nutp3NfB3IF0IUSWEWGAUzaOd6UsIkSiE+MD4GQ9sFULsBv4BvC+l/OhcK6a48JBGFFZ/pmmpq6vj\ns88+I72dg3z4yZNkOp2YPB4s1dVUjB3LPX/8I7M2biS6vp4Sm40rd+yg1GZDM5upsVjazHOy9dpr\niWpo4IGVK8lbv54HVq5kdGsr119/fUA5vD2qI8b89LUzZiCEoFXTODhmDJFuN1cWFrJk5kz2zZsH\nJhOVlZXszM6m1GYj9MwZipon+nxD1WcspIWXtQlASA3Z09Z3YlzPmpqaDn6W3amp/Fr5WRTnge56\nKtvR/SgCGAOcML4PBw4B4zvbUUp5dyfr7w2wrhq42fheAUzuXnTFxYZXqTQ2Nna5XaBxGHaTKWDI\n6y9/+UsKCgqYNnt2mzlKhgOJMTEsfOYZWqKiiD9yBHd0NCuXLMFaXs6NH32Ex2Qi3eFg08yZlFmt\nXLltG8uefpqTw4czur6e0OBgzjY3M2nvXkI9HhISEro0J3l7VLtsNmY5HBQ3NLDqvvt8PZ3F+fmY\nzWZfWv3k5GQu376dUzExlKWmEre9zucb2td0BWbtNDNjN/nGsZS3ZDMjLo5X1qzhxqNHcUdF8b+7\nd3MmNJQUdGXqM7GVl1Prl+dMoRgsulQqUsrxAEKIF4F1UsoPjN83oZuxFIoe441mOnXqVJfbdDbY\nsf0cJ19++SUvv/wyy5Yt41+XLePWTz8lubKSOyZO5Prrr8dkMvHkrbcyc9IkPLGxlNhs+qRd1dWY\nPR72ZGUxad8+XRnZbCQdPkxTeDiVycnEFhXRqGnUx8RQnpJCakUFsqYGh8NBhpGyvjOmAA+7XPw6\nwLTBNX4p7i0WC01+6fntOTkse/ZZEpNmcOPEiezYtpUX9i8iPbKM0tOZxCZPQNM0yqqrORMdTWVy\nMiEVFUSfOkXDiBEszM+n3AhmqIuOZlO7CDWFYjDoqaP+ainlEu8PKeWHQoj/HCCZFEMUb0+lvVLx\n75ls1rSOc5wY/gH/+UZOnz7Nj370I5KSknj88ccJCgriC+PN/BW/YzuOHOEbHg/Hxoxhyu7dpO7f\nT3lKCh6TiYyyMp574AGu+/xzMkpKaAoP54X770czm4murye+trbNfCxLVq1i37593SoV6JgB2dt7\n8GZJ3gX8d10drowMn+LxBAdTmJ3N7SEhZGRkYLPZKCu7npUrV/Lhhx/yz/88meLiYtzh4R3kOj5i\nBEVZWUwz/CxTrVbkOTjpd3W/iULRI3r61B0VQvyrEGKcEGKsEOIx4NhACqYYeniVysmTJ33r2o8Q\n/2THDgrHj+92HMZjjz2Gw+HglVdeISYmBtB7CP45fzRNIyM+ntGhoSS6XL5ewca8PF544AHCTp/m\nvpdewmwMUvSffwUh2oyO18xmHGlpPa6r1WolMyqKxfn55G3YwOL8fGiXJbmhmzE4JpOJ9PR0nnzy\nScaNG8eSJUuoq6sLLJcQOGw2Dhqj/HuiUHKMRdM0xjscTO0kt5pCcS70VKncDYxGDyt+G4jj69H1\nCkWPCNRTaT9C/K1bbiHd4ei0oQX46KOPWLFiBQ8++CB5eYGtsJqm8cqaNcyeMYPgsDAimppwGE55\n0BvjA2PGgNErKbr8crL27/edd19GRofR8ZkOB5mZmT2qq9dxX5ubS2JICLW5uewzQnzt6MrvmKF4\nFubnM2vDBpbn53dIzw/6+J61a9dis9l49913uWLXLn3uF03D5PGQ7nBQ5CdXDtDQIykBQ6lPtdvx\n9GPaF8WlS09Dio8DywZYFsUQx6tUjh//OntP+xHiDpuNvI0buX/VKkpsNmzGOAxt58+5/7NI/i91\nAd///vfJzs7mt7/9bafncjgcuu8hJgZzfT2hZ86QUVJCZXIy8TU1pDmdhDY3U5KRQWpFBU3NzUyI\njPx6kKXTiWhtZdFLL+G0WslyOskcMeKcQnS9jvtjne3jNwbnsgBjcPyJjo7mkWX3c9JVSFbIJyS9\nVUmdeTSN5nDMgJCS6Zs3M8pi4ZjV2sH/1BmjDKX+sp85TYUjK/pCj5SKEKIAYzS9P1LKzlOpKhTt\n8CqVuro6PB4PZrO5g+9BSEmo2czR9HRMZjO1M2YQvm0rWyrjGR9WRq59NZd95xYe+ZdfEh4e3um5\nioqKfL6HpU8/TfOwYUQ0NpK3cSN1sbGENjfz4uLFpO7fz9mgIKbs3s21s2fzYVCQPshy5kwWpaRQ\nW1FBck0NS2bN6tEodfs5rvdXPF014U6nk7P1+1lmzfdNHPZs1UO4zVGENDXxrU8+4VR0NDG7dxM/\nahSuK68ksba229H1kV2kfVFKRdEbeuqof9TveyjwXUDNJKQ4J7wmFY/HQ11dnS/1iG8aYGPQ4ISY\nGP50/fWcqqggtbCQetc+FibpgwBnjtrEqpCHfArKH7vfd+HnEwlraUFomi+k+LqtWzk+YgTzX3uN\n6Pp6KpOTaQkK4oPNmzm2eDHHbDasRvCAcLk40oe0J/YA69qHTNODnkWgxJmZYfvYmngdKfv30xQW\npkeDlZczorKSWadOUZaWxoubNpH51Vf8eN68gPI3dJL2JUGFIyt6SY/+JVLK7X7LZ1LKnwJTB1g2\nxRDDXxFUVekJFfx9DybD93DPnXeS+frrTLXbOVNfT+qwto2pLXwfLpery3NlZmaSafhmzprNNIeF\n+RzxtfHxXP/V52QdLeTG4I+ZeuhLWqKHUS0lLzmdbPILHjC3tjK1H/0MgVLXZ65dC90cO1DizD1n\nskkrK0MCm6dPZ9PMmWy5/noaIyP1BJazZ/P0woUUnjiBw+EIeNye+nUUip7SU/OX/+h5E3rero65\nvxWKLpBSMnbsWA4ePEhxcbEvxbv/oMEpQIXDAYad31pezsR1e33zvXukiX2NEzA5HL6xKIGw2WxM\n2r6dJc8/j8dkIvT0aTJKS9k0cyYSGCGO8+CY5zELjVlyIytcy/nHmGt8UWb+Yc0FHg+R/eRn8A9M\n8A+ZfqldyHR7rFYr2+MzebHyIWxhRexuyabRE0midJExrISE92oojs/kRPhwSvyi2DSzmX02W4dQ\n6Bz0MOIphlJ/zulkdDd+HYWiJ/TU/OU/sv4ssB9Y0OUeCkU7pJSkp6cTGRFK8Rev45g27esGTNMY\n73Qy1uVi7941OG136AMGrVYORo3j+UNLyYgoprAlm+KkDIY3nu2ykdc0jeraWoYBZ0JCiGxoaJNn\na0J4cZveT1bIHo5VWkiYMqXL9PJ9VSq9PbbJZGLu9+fzxhtvsNkRS0tICEnyMA+Mfd43I+WzVQ9T\nGpVGfG1tm5xltrIyaDf/y25Nw2Jcb6fFwgGrlf3d+HUUip7QU6UyQUrZ7L9CCDFsAORRDGGklEzJ\nsnHTteOZONyJfd2qrzPvvv46mttN+fjxbGscy+U7d3IoORmHzcYe2yRq6uJ40/JdahISKLNaydu0\nqcuG2G63czQ4mPWzZxN/+DAzPvuMnZMnc2jcOCbu3cveI5Pa9H5KGieQGPu12Weg/Aydpa7vybFN\nJhONjY0cHzGC2oQEZh/e0EYxTgwrovZUHB6zmaXPPUd9dDTR9fWEtbS0CYXWNI3b1q4lzO2mwsha\ncNv27bwbIGuBQnGu9FSpfA5c0W7d3wOsUyg6JS4ujsQYjaUpLxuNeQH5ruV8+umn4HaTP38+9675\nI7GmI2QH7+Gyt6qoCk9GyGCag4N583vf63FDXFhYiBYSwowtWwhqaaF52DDSy8rYNGsWZVYrI/98\nnOcPPciEiBL2NU0kcmQK98yf75sauKs5XPpCX4/d3NxMmc3GobFj2bM/q+0c9y3ZvD1nDre/+y4e\ns5mqpCQySksJDwpqc3yn00m4XxjxJo+Hhfn5jOvGBKdQ9IQulYoQIgG4DAgTQlyObv4CiAY6j+dU\nKALQGh3NhGhHh8y7hw5lUZaSQur+/Uyo3cdPLCt8YbPPVT1M1tW3UVFTc04N8ZkzZxDBwXpI8bPP\n0hgeTpTbzQMrV1KSnk7YySZqtXjez76DyMREXm83R0tXc7j0hb4eOz09nYx9+yjIyaE4PtOnGPe0\nZFMcPwEpBO6oKF5cskRXGDNnsjw/n4qKCl+vrrq6Gmc7E5wzNZUYv7xkCkVv6a6nciNwL5AE/I/f\nejfwiwGSSTFEqT15kr31Nma1ybybhW3MGNLKymgNCiKrQ9hsEUFB3+HH8+bx4Tk0xKGhob5or/qY\nGCIaG/nd8uVc9/nnjDl0iLCmJlISEvgsJ4czdAyD7On4kd7Ql2NPnz6drbt3+waHJhyopqo5iXVz\n5lBmtepK0y+fWHufjaZp7C0rI6OlhYLcXF/Pz1pezmYVRqzoB7rLUvwH4A9CiO9KKd8cJJkUQ5T9\nLhdl9TbyXctJCdnDnuMpVB1vYEfZGhpDQph84gT7yGSWn0mnqCWbm42U8+fSEGdlZeEsLmbTzJkc\nGDOGIE3j7LBh2I2GM2/9er7lzfN1EREUFMTPHnyQf1q3jsziYhJGjKCqoYG8DRvI3rWLsKYmsvfs\nodLwRwkpmVRWRsKsWYBu+trv8dA8bBiLV63CYbMxobgYKQQJ1dU4wKewc4xz2s9TXRUXJ92Zv34g\npfw/YJwQ4qfty6WU/xNgN4UiIFJK3nBWsnzhb6ipqSGqpIT1bzyDy+VCTplC3OTJ3HrazEuVS0gL\nd7CnZTLV8ROwWq1fh8D28FwzZszgsz17eGDlShrDwghvamoTEZVeWoroZBbHCxlN0/jTm28yyu2m\naOJEwnftwh0eTqTbTfaJGqJCGsgM3sfYtw5SFZFMozmcRMNPBLrp64zHg1kIjsbGMnnXLoI8HtzR\n0ZjPnu1yqgGFoid0Z/6KMD4jB1oQxSWAlCAENpsNm83G9OnTWbRoEaCPm9gMlGsaRU4nv66poSEh\noY2v41wwmUxY4uI45XIR4Xbjjo72TfVrLS9HCnFex2LYe7mfd5yLN/X94cREbnv/fc6Eh5PYVMPS\npOeMzAMbeabyYWqIZ3xaGk6nE6vV6huAmm9MHpZeUkLehg2sWrSoQ+4v5V9R9IbuzF+rjM9/Hxxx\nFEOVHKBRSs4IEbDcO5+HNN6qG9DzUjnRzTHn+tbsdDopaWxkxfLlWMvLySkoYPOMGcTX1WGfMYOZ\ndjsWi6UPNTo/dBjnYjJxOiyMvZMmcfOe933+KIEkRGshicPEHN6IfX822+MzSUhKaZOtOb62tk3K\n/zZTDVwgSiXH+LSfRxkUPadH/1QhxGghxC+EEC8JIV7xLgMtnGKIYfRUukIESGPyag/SmLTHv/Et\ns1pxR0cza9MmQs6cIWfzZoiJ0acppuskkJ2VnS8s7eZgsdTUUJqejisxkUK/NC6ORhsSwQNjXyAv\ndgMLLCuory0GaJPivzY+noxuphpQKM6Fno5TeQf4FNgAeAZOHMWQRtM6VSpTNI1ap5PJhYUUnTrF\nqoUL26QxGXWO5hiLxUJWQQGHkpKIr61l5+TJ3LB+PVH19Ryw2YgXwmcSuphSkrQf53JFURHNQUG+\nEONnDi9j0rBCdriv4IqonR3Ct02m23y5vpypqVidTuJaWny/pxjh2ousVrai7N6Kc6enSiVcSvkv\nAyqJYugjJSJAA65pGplr15LudtMSHExhamoHc4zlHM0xKSkpxLz3Hjd/+CEnY2K4cscOpKYRduYM\nCWVllBkjyb1O6YtFsbQf5zL/xhvZumMH973yCs6UFBIOH6GgMZcaSzwhx1raDI4sbZzAbIuF66+/\nnuecTqJratg8cyaFKSlcXVHRJvfXny6S66G48OipUnlPCHGzlPKDAZVGMbTpxPzlcDgQp06Rv3Ch\n7v+w2zH5pzFxOjmans7UzZuJtFjQetC7cDqdNGoajdHRHBozhuCKCkYcP050fX2HXtDFNiGVf3h1\nBnryzCedThJqaphitfL5oUOMPH6ck9oInq16iImhRRQ3TmC/q5nZxv77bTZO2WyY0RuB/Tabyv2l\n6Bd6+jqyDF2xNAkh6oUQbiFE/UAKphhiaBq2xESui41tMw+6pml8uGEDDqN3Uma14o6KYsmqVeSt\nX8/yl18mpqmJ8LIyPK2txPcwDb3/JF0b8/L43/nzAXAG6AV5MxNfrJhMJhw2G1umTyc2Npaqyy6j\nKTycXz/yc1Z+5wF+PfkxDoaO49O/b+f2229vM51zf5u3cvjasa64NOnpfCpRUkqTlDJMShlt/I4e\naOEUQwOveStnzhxMt9/eZh50p9NJjaaRYjifpcnEX7/3PUJaWrC63dySns7JiAhWGcph1YIFFNfX\n6yGvXeA/SZfQNO564w2ahw0jraxsSDulExMTGVtVRWl6Op7gYBw2G5tzcijMzuahhx6ivLycJUuW\nMHXTJnIKCphqt+s9RU2jgc4VQk4XZQqFPz2dTyVQ4shTwEEppZoBUtElTqdTTxj58MMdxkK4XC52\nZmQwsaiI+1eupDQ9nfTSUkJaW9k5Zw6mzz7rVar4zMxMMjduZNPMmVjLy4lyu3lu6VLueuMNFubn\nU56SQpbTyYThw4fEhFQxxqfVaqU1OpqM0tI2aViyy8uZOmMGT/zHvyKaD5Nd+F8Un55AVUQyq0qj\nuW37dlarAY+KfqCnT9ALwBfAy8byBbAGcAghbhgg2RRDBJfLRVknisFisXB5SQmnw8PZkJdHy7Bh\nbMjLoyEiglEVFR1CaHvau7DZbGSNGMGil17im1u3UjF+PJ7gYNbMnUtBbi7JVVVcFx9/UTnp/bHT\nNtx5irGYTCZ23XMPmhDc/+KL5K1fz0+M2RwBooNO8ZO0fGaP3sCDY54nsamaXampjKir4/a332aU\n1zSpaaQ5HFy+eTPfdTjOOaRbcenS03/TAeByKeWVUsor0Z/fvUAe8NsBkk1xgZLDuZlCLBYLaZ0o\nBqvVSoLJREVKCqUZGWyZPp3SjAwcNhuRNTV6CG1UFIvz88k7h+luTSYTP543jy9mzaJx+HCyyst9\n5jVnaiqhLS1MmjTpolQoXaFpGj988008ZjN1o0czec8exptM3HPnndTW1nac5z60iG988QXDzpzB\n1NpKfEEBr6xZQ+aaNUw3plOO7+GUxwoF9Dz6K0NKWeT9IaXcJ4S4XEpZIboZzKZQpKSkEPH++z7z\nVkZpKTGtraSkpGAymbgpLw+H3d4mN1daebk+b30fUsV7o5wOWK08tXatbyyGrbwchtg87Hbj0+FN\n42KkYTEZpsa5FRVgsXDjtiw8Rphxq2Zmp/ty4s11pIeWUFiZRXH8BMrq6hChoeQbUXIFHg+LezFW\nSHFp0lOlUiqEWIlu8gKYi276Gga0DohkCmBopKioqKigMTycTUaalPV5edxqt/vm+LDZbDTu3Nlm\nAN7p6GiOGY1+X9PQS0MxPeV0El9TQ21uLscuskGPPaWz6YotNTUcvO46ouIzWXXoftIiHOw4fTnh\n5iYeGKNPSTxLbmSFazkV4Sm+aDzvMcpSUxE1NeTYbBf1s6gYeHr6r7oXcAI/AZYDFca6VkBNwqDo\nEpfLRVlqahvz1h6r1RfKazKZeNfwdZhCQliUm8u+fnYam0wmXDYbu6ZP55jNNmQd0p35oBoSEsCY\n5z5j6h183HID/0i7hvSIkjbmsKyQPQx3u7EZ5kLvMVLLy6kcQlFyioGjRz0VKWUT8DtjaU9Dv0qk\nGHJYLBbS7HY2dTEvuzTGWtQavRH7AMjR07T5FzPeNC7+GZknREfzlZ+pr6KmBgGMOXSI4pa289eU\nNE4gYdQITkRGtjEXNkRH47Raie/i3JqmccrpJM7lwmGxnHMKnBzj096LeisuHHoaUpwG/BeQCYR6\n10spUwZILsUQwmq10rB9OwuffhpnVhZTKiraONs1TWOW08lxl4uTPRwxrwiM1wflTcPyZW4uL/il\nXXE6nRQ3NPD0smVYy8u59ZP3eOHQg2SEF1NUn0aVKYYn7r2XuysqkEVFRNfXUztjBptsNmQX90Qz\nEoFOdbtxnkMKnBzj095/l0BxnumpT+VV4N+A36Obu+bz9Xz1CkWXeM1bcddeS2JVFYuWLPG9xXob\no3i3m1MpKXzDbufVfszHFah3Yu/zUS9svAEKDTYb19HWxu1yuSgcP57U/fuJr63lvRtuZfLu3exs\nvpLPdmzC+eeXGBkxkviwMI6OGEFyZSWeo0dJkZLIujpOdaL0vfO8vGzM81JwkabAUfSdnv5rw6SU\nGwEhpTwopXwcmDlwYimGGtJkwnHwIJ81NGCz2XyNkrcx6m7EvJ3eKYPe7jdUiY+P5/Ldu8kpKCC+\ntpZb3n+f1P37OXX11dS8+CJZt93GwbNnORkVxYiTJ9mdnU3IqVN8Y9MmgltbmdpJmpzOAgQu9hQ4\ninOnp0qlWQhhAsqEEEuFEN8B4gZQLsUQQtM0xjscTL/iCtJCQ9s0SIPZGNlRCgagOTQUt1dpTJ5M\nfWQkSZ9/jmnYMLKWLOFEcjJRDQ3kL1hA5dixnIqJYeWSJWzMy+PlTpR+bwepKoYePVUqPwHCgYeB\nK4EfAvd0tYMxkVedEGKv37opQogvhBC7hBBfCSGu6WTfbwkhSoUQTiHEz3oooyIAOZx7zqbe7NMZ\nPlu73U7wsmVcZ7W2edNVjdHgUltbS01cHFFuN5unT6c1OJgNeXl4zpxhnNPJ6fHjSa2spNxQ9Aku\nl+87dK70vYNUF+bnM+scBql2h6ZpjHI4GLt5c5tEpIoLl54mlNwmpWyQUlZJKedLKe+QUn7RzW6v\nAd9qt+63wL9LKacA/48Ao/GFEGbgeeAm9MCAu4UQmT2RU3Hu5NB7BdKTff1t7RtvuIH8Bx9s86bb\n2xHziq4RRpoVb2OMnxJPOXSIoLNnmbFlC8GtrczYsgUJjHa5OGa1cllMDBmlpZg8HmosFlJ7oPS9\nAQJf+oWF99Uv5vO3tZsFVCmWC5suHfVCiHVdlUsp53RRtkUIMa79asCb3TgGqA6w6zWAU0pZYciw\nBrgd2NeVLIoLk67MW17fSm9HzCsCo2kat61dS5jbTYVfJNa+uXNJSUkhtLUVd2go+X5O9ftXrsSk\naWAysWT+fB75wx9Y/OKLOKxWok+c4P4XX6QkPd0XohxI6c80mSi02Yjsp3lZ/P1tF/P8N5ca3UV/\nfQOoBFYDX9L3iK+fAB8LIZ5C7yVdG2Cby4xzeqkCpnZ2QCHEImARwJgxY/oo3oWFpmmMcjqJ7GXc\n/0DL1GCx6KPeu5DJYrGQZbezsYsxKn0dMa9oi9PpJNKvMS5obWXpCy8w4e23+XTECJpDQihNT2+j\n6EvT033PVlBQELvnz8fldJJQUMC7f/gDIjOTuIwMX4jyYDyH3b2QKC5MunsyEoBfAJOAp4HZwFEp\n5WYp5eZenO9+YLmUMhl9ZP7/BtgmkOKSnR1QSvmSlPIqKeVVo0eP7oVIFyYXYtffX6bDra2MNhIN\ndiVTG1v7J5+w6PnnlXlrgPHPCu2dS6YlOBhnVBSf7NhBRVISKfv3tzFppZeVUWexfH0Q78Rfixdj\n+vnP0QoLYc0a5FljposNOfoygCh/28VJlz0VKaUH+Aj4yMjzdTdgF0L8h5Ty2V6c7x70WSQBXgfy\nA2xTBST7/U4isJlsSHMhdv3by+RNNNiVTG0G4z30ELsnTeL5pUvPe49rKOOfwcA7l8yqxYvRzGYq\nk5OZvWEDx0aO9I2Yn1BSgiksjANWK5O9B9E0bE4nCdXVTGlpofnhhylPS2Pq+vW8smsXPx4NpgEe\nqebNDrA4P5+y1FQmd2F6uxDJMT7t51GG80G3gx8NZXILukIZBzwD/K2X56sGZqBf55lAWYBttgFp\nQojxwGFgHvBPvTzfVHwuSAAAIABJREFURcuF2PUPJFNZD2TyzYn+5Zc8nJGhFMoAY7VaudJojJtD\nQqgYP953zxw2G3kbNzL8xAmOjh7NlN27GRcdzRv33OMbMe+dqTPd7ebE8OF4mpvJX7JEf5GYNYvg\n557D3jQWswksDgfWg4t0BZNn79d6KH9b38kxPu2DeM7uHPV/QDd9fYgetbW3q+3b7bsavU6xQogq\n9BH5C4GnhRBBQDOGL0QIkQjkSylvllKeFUIsBT4GzMAr/mn3LxV64osYSHah3zx7NzKlnYtMUqKm\nShh4/BvjuL17ySov900rIKQk1Gxmu+FTKZsyhe1WK2/4NdS+mToXLOC6rVs5PnIkSElOQQFjDh5k\nRHM9RQ2RTIhyUPDOKjaLKaSNPMEohwPRjY+tN3VR/raLi+56Kj8EGgEb8LBfgyAA2dU89VLKuzsp\nujLAttXAzX6/PwA+6Ea2Ic2F2PX3l8lhJCvsybwkHo+H5sceg/p64uO7Skmo6C98jbHVyg/Xrm3z\nHGVER/O+xUJcbW3Aff19MjUWCzmbNpFVWEhrSAh1sbGM1o7wwLgXEUjqqkdzomUkZxuquX7dKkbG\nZ/LG9+cP2SzQiu7p8s5LKU1SyihjifZborpSKIq+433brM3NJbGXcf/9PXDMXyZTSAhfGinqu5LJ\n7XYzZ84czvzmN4QsXsxPf/rTPsmgOEfaPUf3zZgBwHRjjMo3jACQTZrmy5PmP1NnmdVKSEsLrcHB\nrFyyhNqEBDKMdPnOxjQaPFEsGbuSvNgNPGRZQUZtMePajbZXXFr0NKGk4jzQl66/f6LGsnPIGusN\nGZ7qchEZIHmgV6Zdhg+lq3Tyhw8f5tZbb6WwsJDQF19k2OLFDOtkW/s51k/Rc/yfI5PDQbGRgsU/\n8aPD4WCUyUSky4UWHw9+qe+jGhrYdvXVvp5L4bYsZsmNuM4kkBJW3mE+luyaGjVD5CWM6qMOUfwj\ntTbk5bGik5xN/viHDHuMucl7G8a8Z88epk2bhtPp5P3332fY4sV9qY6inwgYAJKSwnsFBb7w9Ze3\nbAEpKU1L47LDhzn9/7d35uFRlVn+/7w3G1kRAkkqECBJpbKwu9vIEgguiG27ICjjKDuoKPRuM/Pr\nnsV2mZkeRBSjEXRoARVb29ZulwAVwF3ZE5JKJWwhlbAIZIGQpd7fH3WrqCoSyFIJRXg/z5MnVTf3\n3jpVqXvP+55z3u8xGEi1WNCamrAmJXEgciAvH3iMirOxFNam0SQdt5EmqbGnfpijIVgbUXIs3Qc1\nU+mmtKV6zJmUf7WZkuH2lDF/+umnTJkyhaioKLZu3crw4cMvfpCiS2iu2GJkQQHlISEeiyUXv/AC\naVVVHI+OJubAAUKbmljw8suE2evoV3uImOAKDtUNhMBwXiqdS1pkMXtq07EMSHO1gW417ZxVK/wT\n9R/rprRn4ZgvFINzcnK44447SExM5Ouvv1YOxc9oTmstLiCAPUaj6/+evG8ftaGhnOrZk14nT7J9\n5Eiqw8MJO3OG/qfLeHzgS0yNf5fFiX8iIqCGIaOm8PnhkWT/32dszbcycMuWNs02otsxq1b4L8qp\ndFPaI9TYkRXMdrudJUuWMGfOHLKystiyZQv9+/f32ftR+IbmCkBuz8ry+L8byss51qePS/5+Q1YW\nKx59FCkEaeF7PXIoxpBdBAYG8uyzz3LTAw8wqm9fDtfXt0kBIkL1YulWqPBXN6U1C8e8k/JJSUke\nJcOmVpYxnz17lkceeYR169Yxd+5cli9fTlBQUGe/RUU78S4AsdvtZGzf7io7vjo/H+rr+WHkSI8b\n/YGBA9m7z7OnfUn9MDLj4igpKSEsPt4jdLqolaHTmjasyRqn/zZ38DPw1Xlawh91+7oK5VS6MS1V\nj40DsNt5yK2Nb6zZzJs//MDDU6bw99JStIoKKjMzmXGRi6Hhxx+55e672bx5M8888wy/+c1v1AJH\nP8Hcyv28ByAzbr2VzzdvpraoiE2Zma4bfe8ff2TfVcnklM0jObSYkvrhRMWmYzQa2bJli2ttC5wr\nACgrK8NoNGJdPwVbbQSGUUvOu8EeNxqZ5GdrsjpCeysvuwvKqXQzxum/zRfZL7qFpHxpaamrZHgE\nF46P2vftY/vtt9O0bx9r1qzhgQeaX+9qt9tJtFqJuQJHbZcL7gOQNBzh0/lvvsm87GyKTSZSSko4\n2bMn70yZwtd/mUpFbTyZo+a6/pfuemNOJ5S8axevfPMN5QeLOV0ZS1JoMeYPs/khNoOp02ec+w5c\nYFY9DkchyYVK1/0Nf9Tt60qUU7nCGIfjIr3rQnFso5FEq5WBF3ACjd99x+nJkwloaGDD558zZsyY\nZl/P1fmxuhrrFThqu1zRNI0vR40iLT8fY1UVlWPH8qHJRLimYep9HFPv4x5rUYxGI/zwA3NycijR\nZxtXCUHe7t38ZEgfFiZn6yGzjeTYFmO1WjEemMtvT0SzNvYJrHobhebWZNXg+M76inFc3FGN03+b\n23F+f9Tt60rUVe3nmGn7F9tZItwSNcCuFpLyMTExZOjtf1uS3D/24YfUjh0LYWGM/PLLFh0KeHV+\nVJU9lwwzrf8eOQcCozdvpqpnT8SxY0zavp2IixxXOXIkJ/v0wVhVxeyxY/ndokU8+eSTZERZPJL7\nycG7KN/yNC/vHI75wHWcqawke9MmMt5+29Wh8nLmSpfsVzOVK5T9XjpezqQ8ALoTaG7qvnz5cgoW\nLSLttttIePRRIjUNu93e4qzjSh+1XY44BwLeq+4HWa3sM5nOUyN2zyFYkpLoXVrKF9u3YzKZGD9+\nPBveL6RJbnIl94vrhlFfK9BkHeMjNtL34FH2xqYTfOoU0VYrmEzY7Xas66ewoDaCylFLqPCxUGVn\n4o+6fV3J5fFfUvgc6aXjValri1VWVp6fcE1O5vDhwzz22GM88cQTPP7LX5J5553Yg4Mvuur+Sh+1\nXY60NBDo20KJr3sOYUNWFtmzZrH3RCXW9VMwGo1cFT+UVw/O47OjE1haPIdT9WEEyjM8OvBlbunz\nOYsMS0mvLORY795EVFRgt9t5+61VmPfFEl1dzuwPs5n61irsdjs78G0oDByz+nE+PJ8vdPsuZ9RM\nxY8Yp/82d9Hruet4OZPyzSVch1qtvJSby3vvvcdvf/tb6hMTeamVq+6v9FGbv2K+wN9aaruwoYUW\nB806oZTBVJRtx6RpTJ0+A6t1NBaLhe/+/GeGp4WT2qvwPM2wikMGSkeMwGq1Ul1ZwKz+jjxMptzI\nC3oeBpOJGs5vy+BrOuq4rmTJfuVUFB54J1yHFhdzoqCAjz76iDfffJOBAwfydEMDUghMFgtxNhuH\nevXCZrM161RUo6XLD+dAwCko6QyNfm80Npvc9nZCAQ0NjNyzg2PBYVgsFoxGIyaTCZPJRHBwMCd2\n5VB6JtkjJFZUk0ZD3yiOG43YtmwhKXi3h9MZErzLUUTSzHdsnP7b3FkfiKJNKKei8EDTNAqmTuWU\n1YoxL4+N2dlUVlZiNpu58cYbsVgsDN20iT5HjhBZU0NJUhIxR4+yq7qa0aNHN+ssruRR2+WIewvo\nnm7rlVa3MBDw6LOTlMSoH74kRh4jKsB6XgnxoEGDKN2ZQFSAjZxDc0gOtVJYm05ARAI7Zjj6sBgM\nBszfDaXJbZFlQf0wJsXFUQM04fsQmMJ3KKfSzTE3s60JRwVYS0ghKPzgA3b89rdcPXIk3377LQkJ\nCYDjBpKYl8fZEydY4Wwxm5l5RdXhXwk4W0DXmEzczIWTr+6z0eg9e+itnWJ+v5fPKyE2mUwYjUZ+\niMvg1OEz9Ao8xncnRmCraeBXD00j6S9TOVUbgf2mp4iITXctstxVP5xyfZGlL9iB5/dfOSjfopyK\nwoPq6mpiZs6kbv16gu6/ny2rVhEWFub6u6ZpDElJ4a+NjaqiS+HCORsNsdnIKNtzXgmx87uhOXMs\n6/9ORW0Er0UPJm/hQuw1VaQPjCE13Mrmj14jIiadsXc9wTNHjrArLo6eukNJsViIsdmoaqbXT1tQ\ns53OQzmVFhin/zZfQhu6moKCAu655x6sVis9/uu/CP7FLwhrRnIlPj7+vGR+S1pNiiuLIwYD+d8N\nZXwz+mBONE3DdP97mIBGYGJcHMnf/pX5A151zW5eK1+M7YutTNIku+KWOOTx332XMdXVlCQlccOm\nTbz88esM6VlB9M1LHHL7Kk/nFyin4seM03+b23hce1YgH1m3jutnzyYiIoINGzZwl952tjmaWz2t\nKrq6P+ZW7LPfaKQsNqNZfbCW6BsWxtBepa7ZjUDS1FBH4ZEIUsKLeOjDbI5FDqKAYHJmz0YKwSDb\nfuSPZ2gMKGf0h9mUxWZgd5d+uQDeod8LhYIVbUc5FT/DuRre3I5jncqoo202TrYyPGCvr+fMr37F\n3mXLGDVqFO+88w7x8fEXPMY9mT9MVXQp3JCaxvfTZzB3/d/P0wdrjh3AKYOBUrfEvKXWhATmDljh\nKileXvYEXybehD0gAJPFQnrlXh4d8JLr7y+65W0UlxblVLoJrlXNVVUE9OrF1d9/z8tRUcyfMYPA\nwOb/zQcPHmTnAw9Q/+WX9Fu0iE3PP996yXo9kdtTVXR1W0bQzrxDC/pgLbHfaCQyNoNXy+aRElrM\ntqqrGRm5zSMvkxGaT8Wh/qQWFnLDN98w1KvkeLBb3kZxaVFOpZtgtVopqKriZGQkvU6eZOfw4dQX\nFbHizTd5bIYjLODsA353fj7i0CHuXr2a08XFhK5bh3HqVFQHFIXP8JJyuRBST95PsY7m64oKHmlq\novDbs4yXG115mb216QTKs2Rt2MCRPn0oLE9jQrRbyfGZwQxrarqgZBA4Bl8mi4W0/HyChaAgLY0U\nTSOuslIpaPsI5VS6CTabjbJevbjq5Mlzmk2ZmTyRne1QhDUaWbluHUnl5dSGhVE8eDA3PfII14WF\n8faUKZfafMUVjvtaprF2O6WWPbx08DHSwgvZXT+MA70HYG/QWDF/PlIIot6q5qWDj5Eetpe9pzM4\nHB7PQYuF0oqKFiVR7HY7K9etY1J5OafDwig2Gnnj448ZHxZGUWqqh4K2L5L+5g6f4fJEOZUL0JH8\nRlcxTv/9qsFA0vff893w4R6lvvkmk6sta9HRo1RFRbmczsbx45n/6qsM0uUvFAp/QNM0klOHsqIy\nhvcM91ERF4fh8GGC3MrYV06fyYxVqyj/sR8f3H0XFpMJIeUF10tZrVaKjx2jWr8GjCUlJO3ff55w\nplVdDx1CzfO6CUajkX5RUaQVFZ0n3hgbG8tHH33EschISrw0moqNRoapXuAKP8DMuQFcfHw8/U+d\nYuuoUVhMJjQpSXX7bgspiT5+nG0jrqYoLQ2paa5uk67e9rnjHD86NpuNo27XQJzNdt714OoppGg3\nyql0EzRNY/6MGQT06MG87Gwm5OayOCeH5B49+PWvf012djY9jx/3uDC1pibSioqoiYk573wjuLy6\n7Sm6F0ajkYzISObl5DAhN5drfviBoIYGFi9dypR33mHOa6/RGBREhtXqOYgqzicuLg673Y7lx2jy\nDg3EYrFgt9sxGAz0ra4mWVfNrjAYSC4pIbWwkDF5eaQWFjLMalUK2h1Ehb+6EYGBgeycMQOb1Yqh\nooIhISH8etFcTtU2MODZ/yYpPJyyqiqXUKCxpISwM2eILC8nQiUqFX6Eh/6YzUaPEIip3UdqeCEF\n+zIovSqJM01NJDY1Ob7PSUlcs2cHcfWVNG57inXf3E5NuWcL4ykPPExKnz4UlZczPzsbS3IyvU6c\nICs3l6LUVCZu2EDP+nqSkpIu9du/rFF3jw4yDt/2YugwmoYlJobNOTnMnDmTa0eksfS5aQybPJmh\nJhPbR4xgU2YmDcHBmMeOpSYsjIiioha7PCoUbWEEvstBOvXHKg0GousreHTAS9zS53MW9n8R43Er\nIZWV/NtvfoP51ClGFO5kUN1++gUc5G/FKZwo28Ws/tlk9clllmEpVZV7KS0tZea0aXw0aRI1YWEM\nKSykLiqKFQsWkDtxIi/Pn8/JsDBKS0t99A6uTJRT6Wac2LABhg2DtWv5j98/xV23jSIl4CijP8ym\n1LKHlNJSrMnJbNZbAAtgxdy55KpWvwo/Jc5mY7DXupT0iEJC4+KYuXAh8vXXiTtRSs/AE5xs7EVM\ncCWmsL3N6o/Z7XZ+8tVXhJ85w8nISPJTUjxyKruTkrDZbJfsvXYHlFNpAW8lU3+nqqqKxx57jF1Z\nWRAWhun99+kbXsdsfbS20LCUxuoDBGsac/Q49T0ff0yRyXReotJmsxFtsTAwL88Vj1Z0f8z4Z6Vj\nhcFAfv1QmqTjdtUkNXbVD+P9O+4gJDGRIbNnkxB2jJqmSGYl5DCq15fsO5PksX9J/TBiYmJ4ZdUq\nQuvqyJ47l/2JiaRaLB45mdSiIux2OzX43/U/jrZFRZzr0rr6OlY5lW7A8Y8/ZvD8+Rw+fJiQRYsY\nOHMmN27bdl6jo+TgXXxpupOvDQYMFRVMvPpqzhYXe3b4s1rZFRREbFERxUlJHrX7KteiuBRYdT2x\n7IMLSAm3sKt+GHtj07GYTOwqKyOyqorSwh+5Lvw7AoQdY3gxP1RdQ/bB+aSEWShtcOiPARyuqqJQ\nL7u3C4EEjxyjFKJbfM9dChvV1V1+HV/+n14n4WwG5M8cPXqUvdOns2fyZHr27MkXX3zB/ddfz7iv\nvuJ0aCh7a9POG63VxMYCjpJMg8FAekQE83JyyNKrxQYFBLC/qYnsWbNUSEzRphyJuQ37tgWnnlhq\nTA25Zyfw2uTZrJw+AyElw0pK2JuRQXnvfhTp33dNSO6Le4ez9lCqExaQ+dO5TJ0+g8rKSkoTEs5V\nf8XH0xgYiHnsWFeOsUdAAMtiY0m2WBh1Gc/UrVYrBdXVl+Q6Vk6lA4zDdz0ZnFPVkXl5nLrIF1lK\nyZo1a8jIyODou+8y8A9/YNu2bURHRxNWXU3OrFmsv+8+dvYfxksHH+PzY1m8aFtMZGwasdu2Mc5s\nJqChgdc2bwagcuxY4oODmZuZyZCUlPP7jTdTu29GlRx3d8z4UTjMbueLU8OQaEzYuJEJGzbw2IoV\npEVGUmwykT17Pj1CAsg+OJ/Pj2axvHQu+cXlvPHGG4SFhaHpHSX7nzhBdUQEc3JySDhwgJ5VVWTl\n5hJUX8/kvDzSoqJc10jQZVy8YrPZWnUddwad5lSEECuFEEeEEHvcto0QQnwthNghhPheCHF9C8fu\nF0Lsdu7XWTZeiK6cpbimqvrN/gb9i0wzX+QbSkvpc8cdTJ8+neTkZK7Zvp1Bv/89wcHB2Gw2rImJ\nGEtKGL1lC1uvu5kvkkZhsWew5adzufq6m6GmhpxZs9igj14Kq6tB0zgwZgwmk4n4+HiG6iM5OLeA\nUtXuKy4lg7Zs4UhICH9atIjcrCwagoIIrK8nKS7OsfAxMJAZw3aSlVxByMh/5fb7nyB96DXk5uaS\nkZHBK6+8QlJSEumRkfQ8eZIfr7qKETt30j8qitwJE2jQB1U3X331edfI5ThTNxgMl+w67sycyhvA\ncuD/3LY9D/yblPIfQohJ+vNxLRyfKaU85itjnC9i9tUJfYj7VNUpFxGRk0O01cpxXS6irq6O5557\nju+eeQYtKIilS5fy+OOPM0EfiQDExsYy8qOPSN63j5KkJMZu3kyvUydIGlDGuyYTlXl5FDczejFU\nVLhex73feLEPe6WYO3S0ojswTv9tbsMxwm4nxWoldtcutmVk0BQUhMVkwmIyEdTQwIiyMte+msCh\njqxXNqalpXH33Xczb948FixYwJo1a7j3n/+ZoMZGwk+fpj44mLDISCypqRSlp2MC8lq4Ri43BeTO\nuo5bQ6fNVKSUm4EfvTcDUfrjnkB5Z71+RxG6mmlXVE60NFWN0Keq//jHPxgyZAh/+MMf6POzn3Fd\nYSFPPvkkQgiP6g4pJadDQ12jrJxZs6gJi2DJ8KcBx+glpZnRS41XV74ZU6dSmZnpCompJL3iUmC3\n27nz7bcZYzZzNCKiWTWIhIQEIoAIcCgje6kjJyUl8dlnn7Fq1Spqa2spOHWKlx57jDdmzOClxx6j\nsLqaFLdZiKFyGUMsey77mfqlvI67uvprEfCpEOK/cTi0n7SwnwQ+E0JIIFtK+WpLJxRCzAXmAgwY\nMMAnRtrtdqa9/TaR1dWUtrdywqk51AoJcIPBwFCz2bMKq6SEoiFDyL/nHia9/z6pqank5ubyHxMm\nuGz0ru7oX19PUVqah3MqSk2lz5Ej1KeleXRstCYnM0IfvXzvNXpxV4y9fMZmisuRC+XlrFYrYdXV\nvDZrFkjJEy++yKMrVlCYmkpqURExDQ2MHj36oq8hhOCRRx6hb9++LA0OPm/wFldRgcU5U+91nIFH\ny867Ri7HrqaX6jru6uHnAmCxlDIBWAy83sJ+o6SUVwO3A48JIca0dEIp5atSymullNf27dvXJ0Za\nrVYi9YR3eyonmtMduhDeOkdzcnLo8eOP/OPWW/nx00955pln2LVrFxN0h+K00bu6o6KpiVQvLaSU\nkhKO6qMsZ8fGTZmZNLqNXlRvb4U/YrPZsOozeHtgIMsWLqS2Rw+G79zJmYgI7rvttjaNvFNSUhi2\nb5/n9ZGfT4UQgOO6tZ6IJi70KEUmEw1BQWqm3g66eqbyMPCk/vhdIKe5naSU5frvI0KI94Hrgc1d\nYiGOL3NL6qUXi6va7XbefmsV1V66Q1Mv0D/bOVX9uLiYhs8/Z9O6dRR/+SXRd99N8p/+xG8HDnTt\n66w2ay5ktj09neF79zIvOxuLycSIkhJOR0Wx32hk+LkXw2IyUaJmIQo/x2AwYDSb2aTP4IWURJ88\nSU1EBAf69ydnyxYyduxATJ2KbMVN35lncK1LKSzkVEEBxU89hXjqKdYOiKO2wnHd/vS7v3E4NB57\n7GTPk7QhAnEh9LN0y1xjVzuVcmAsjs9yPFDsvYMQIhzQpJTV+uNbgH/vSiMNBgPJbl9mV1w1M/Oi\nx1qtVqorC5jVP1vvSreRnFb0z/7222/JXbyY6q+/Rhs5kg0bNrC8Xz9q9u/n3rNneddL6NE7ZBbQ\n0MDVu3ZRHR7OsehoRuzcyaCoKNZPmeJxwZlRK14VlwdGoxFNdwIlyclcnZ9PTWgor86Z49H/ZJDV\nyr5WJNGdg7fnrVbiKirIu+UW/n7HHXzX0MDg7ds5ShwLk7LP9b0vW8iqzz7DtH27Y0YPWH+MxlYb\ngcFiUeKrLdBp9xchxFocDrmPEKIM+D0wB3hBCBEI1KHnQoQQ8UCOlHISEAu8LxxT0kBgjZTyk86y\nszmMRiPVbl/mtlRO2Gy2ZleytzTLOXjwIE899RRr1qyhT+9gHv1FKq888w2l69cTW1zMqaQkYs1m\nVnl1pPOu7rg6P58zoaG8MneuyxEuzskhurS00xoOmTvlrAqFA2e49pTVyrCKCm5ISODlqKjzIgh9\nKypa5VSc53RWjwXa7Wx8+23G3Xgj9YGBDK5e63HdZoQW8MGwn9FUVITFYmH7d1+0KQJxpdJpTkVK\n+UALf7qmmX3LgUn641I4F625FGiaxrqpU0mxWrmxooK5mZnnjUrsdjvRVis32GxEGAzY9b8bDAbM\n3w2lSZ7rn11SP4zMuDiPqfOJEyd47rnneOGFFwBYsmQJ16Xlc7Ihihu3bqWgqors2bNb7EjnHHX9\n3WolXr/gXmnmgjNUVJznVCI6/yNUKHyDrlTc02RiiMXiEQ5zRhA2tCKC4E4AjnVoRj0v6ewCOeTD\nfMbLja7rdnf9MGzx8exqbKSgoKBdEYgrEeVinbh1ibPrtfFxNhs1cXHNOpRVb79NX7OZpoYGx0xC\nX3VrNBqJjM0gp2weuceyyLEtJio23TXLqT3TxDPPPENiYiLPP/889913H3v37mVwaiLlhw3EV5dz\n+7d/JeBsFdKZQGxhNayzuuPAmDEMGTKk2cVONTExJCpxSEU3wGg0cjoy0iWIujgnh3Q9Z9geBrjl\nJYuNRg5EDnIoUBzNYqltEXtj0ylJTCS1YA8NDQ3NRyC+eNqju2RnsgPfKXh0Jiq87oXdbmflunXc\neuwY1ZGR9Ny5k5V9+vDnadNA0zBzrvLqtWZ6W5tMJqZOn4F1/d+pqI0nc9RcjEYj9fX1rPgH5OYH\nYd33BmPGjOHpp59m6NChWCwWqisLmKuPgsbLjbx08HFMFgtFaWkuB/FCZiY1ND/TcA+HWXRxvLTI\nSE5v20bfmhpK3EqjW5vYVCj8CU3T+NvUqcRYrQx0iyD8op3f5YMGA5Pc8pK7TEOoOBKD5aiJAJrY\nH5vAwuXLCTx6lNWf5XLrT9Lo3+MglWdjiQ2ppLhuGGnhx8g7NPDCOZYWkvs7cOQHPLe2jL+pJreE\nciropYR6Aq7JbKb48GHO9uzJoYQEgktLsZaXE22xcDwtDbiwro7JZELTNEz3v4cJaGpqYu3atXy9\ndQODYlNZdHsJ1rpJXBU/hMGDB7vO5z0KSg/fy50ff0z/sjJMLawnccc9HNZYUcE3mZksstvZtnmz\ny/ltbGNiU6HwN6SeE6nsQPWiU4LJ6lUNdk1+PnUBGi/Nm4+xpITJn33EwLqDpMXuxXBdClX1oeQe\nyyI1vIjco1nUyh4UHwtTORYvrux3z7kSYPO+WBqqyyn89q/0aDzN6zNmuFalV4eFEVNQ4DqmNbo6\nTmcydOhQ/v3f/50kQwiLUnKY2HcDc/svo6pyr2vdi8FgoNSrX0Rp/TCOXn01WnAwla1cT+IMh20d\nM8bRMa+yslnJib5dICqnUHQUp8jqDXl5RHdC6DZcH4h9k5mJFhzMjFtvZWjvaGavXMmwXbsYUHuQ\nRwe8xMS+uUwybKRncC0LBq5gYt9csvrmEq6dPK+7ZLNr2U7scPxcIVzxTsW9BDirTy7zEpbTSztB\n8r59gONGbElJ8TjGuVjRGdudp8d2jUYjjY2N/PnPf2bw4ME8+OCDaJrGE088wdCrSpqtCHOeLzI2\ng1fL5vH5sSwbdKdEAAAeD0lEQVRe0PMw+8eO5dsxYziuz37aSkuyLJFxcUphWOEXtJQncBdZdc9b\nCi/HYrfbSdQdT2tzhgFuj50ti78dM4a0tDRmTpvGN5mZhNbVkRZe6LpmK8/Gkhpe1OJz72va3T5L\ndSJ5R2/wbU7TLQfsb1zxTqW50NPgHvkYyh2yZFpTE+lFRRxJT3cdo7mNcBr1mcRD997L6tWrycjI\n4KGHHiI4OJh3332XXbt2MXHixPNmIiX1w1wzG03TmDp9BomJlZRHxrNa7//Q0ZXuRqMRdOeX5ZbY\nPH4ZSk4orizcFSM2ZGWRratajLda6YlD3sXpeG7QHU9rZeqdun5OR+TuqJxO5usbb2S32zUbG1KJ\nxa0/UWxIJUW1qR7XdNHpDGJiYlznstvtrHtrJRsrb6S+wc7GD15h3Vsru32xzBWfU2muBLiwNp3h\nO3cS2NiIsaSE0Lo60CuxnDi/fNUDB2JcvZr022+ntLSUESNG8Je//IW77rrLNbswGo38oFeEJYcW\nU1I/3KMizHm+uN7HCe19nH0mU4e8fQTnNJUqR44kKD8fY1UVs8eOxWQysboVzsrcgddXKDpKS3lL\n9xL5ixXMNIfdbmeKruvnLF6584cf+JtXeLnYaGSv2zVrPTuc+qCeHs8bgnqSc2geSaEW9pxKIb/0\nCC/fey9Llixh6tSpFBcXc+LwbuYPeNVVgPPKocexWCyg52e7I1fmTMVt6uhdApxd9gSHwvvz6a23\n0hAczKbMTLYPG0ZNZaXHNP3kyZPUPfss9kGDsMyZw1VXXcUHH3zAtm3buPvuuz3CVc6ZSGZiJcGR\n8a5OdN4hrUVZZiZ0UP7BhTN8sHkzx3r2RBw7xhfbt/vm3ApFJ9NS3tJdUbs9jajcdf2c/VLCqqoY\n5JULkZrGSrdrdvxdc1m4+KnznydVEhLVj9vvf5Ks236KpmlMnz6d9PR0/va3v2EKLfCIgphCCyhw\ny8+2FrvdTorFwui8PAoLCyk83nptwa7mypip5I5jKY6btjfOG76zBDgmfhQFFRUUp6S4ynlv+ewz\nigcPpga4qayMUUuXkp2dzdmaGpg4kWGrV/P9hAkIr9mM9+s4K8K6guhmerR4L6BUKPwV7xJ5ZwXk\n80Yj4/V9WlL3vpCcUnO6flZn8Yp+XdTa7Zj0dWqMeIab3UqFTb2PY+p93LWv+zWdlpbG/fffzwcf\nfMB//Md/8NVXXxF+XYrHgsri0ybaKqLvDPONqa6mJDGRdz54h2gZT2p4oV9WnV0ZTuUiaJrm+rKY\nA/tAZeU50bmSEqQQNB07hv2RR/jmrbf4Tkruv/9+/varX3Fm5Eh6AS27kzagL7ocYLNhMRjAaGx3\nXiXiQqM45VQUfo57ibxWUeGogPRaB+ItENkamfrmdP2MJSV8ozui5sJjHm0vLhJJ0DSNe+65h7vv\nvps1a9ZQVrKDnENzSA6zUnLaSG1TBBkZGW36LKxeK//Dd37CfMPLfruy3z9cWxdyMVn6+Ph4QgID\nMY8dS0NwMOaxY5H19Rz8l39Bvvsu8QsWUFxczJo1awgYOdKndmWsW8fETz4hvrSU1z/5hIx1686r\ndvHGjGf+wwycBJZcwnaiCkVruFjllrNEvqUKSM2rJLg1MvVGo5Fqr1X5p91W5TcXHmtPO2EhBA88\n8AAJg1I50xjCwTMJnKwLwVJSxjPPPEPN7t0XPoFbiN49zBdnszG0BW1Bf6HbOxWnEyk6NJDowkLW\nvbXStSbF/GE2b7+1yvFl1rvGGY1GTGFhjMvNJaiujnEff8xpq5VxNzfy19UjMC5bRmJiIqC3OvWR\nBIrFYiG8vJyzPXpwKCGBkz16EF5eziCLpV3n8+7R4l72bEYl4hWXlvZWbnnjXhJsamXp/baRIznS\npw/RevHK39wUJi7U9qLV6A5B0zSm/dNMJsVvxhh5mJ9MuI+wqD689957VN93HyNnziRo9WoaGho8\nDvce+MbGxroGiBUGg0dVmnclqT/QrcNf3r1NbvxgFSeoY35C86JwBQUFvPjii6xevZp+/foRd8MN\n3JCZyfSFC/l2c75rJa7z3He+/TahHekO6UZ+fj61YWHkuOVA5mdnYyoooNJZKeIVHruQ9HZrwgcK\nxaWiPZVbHcXpyMbqoS1jaamjeMXt9TrS9qI5NE3DFLkPU+Q+GD2aUaNGkX7NNeyoqsJqMpG0ezfz\nFi8mbcAAZs2aRa9evc7rxxQRm056RIQjzJeURG1jFCsOPkpaeGGzlaSXmm7tVLx7mwQeb6DBHnTe\n1PHzzz9n/vz5bNq0iZCQEKZPn87KhQuxjBhBHo4v47EfozlVG+FY2Ws0erQ6dZdAae9FIYSgOCXF\nY4RkSUmhd3W147ndTsbbb5PaUqy3GZzhgx0mEyO4AqalisuGi0kddQbuuQmnI4vwki0yGo2caWOe\nxh13ySeXHtiUkx42FNfV8dqCBQ4bJkxg3osv8vqKFfz+979n5syZDBkY7NIBdA58xybso1cv2Bvy\nJPf/bCrsXMKR0+e0Bf1psOg/lnQC3gsb40NsFJ9O8Zg67jqRxLJly9i3bx9//OMfKSsr4/XXX4cR\njpUeztmObV8s8dXljNJDZocPH+74NNmNjIwMUi0WjxxIalERVVddxcC8PPLy8sAHsd7mMKPCYYqu\npTVSR77G6cikEJgsFm7eupVDvXrR12Zz7eMUrdyUmUlDK/M0TrwlnzzC6142eFSfDRnCf/7nf/Lw\nww9TVFREcjM5kyOnI4jufZzN+sr/tOjjjEk40OqQX1fSrWcq3gsbE8NKqD1yF8tL55EeaWH3SSPH\nT4ewbNkyJk6c2Ow/xznbcVcQzrEtxmazkhrWj42ZmT6ZJptMJuzbtjH31VexGo0MLS4moKGBgVYr\nJUlJfLZzJ8VDh2IsKSHOZqPCYGBXUlKnjuwUis6iPZVbHcVgMDB00yb6HDlCpK7cHXP0KD2rq/lu\n9GhXpaVTtNLSRtHK1nR9bakMenhmJlOmTGHbtm18+ekbTJCbXGXI+VWpxNQfp56BmPRIiX+5EU+6\ntVNxrGRP59WD8zCGWdh10kjRvjKKi3Zzy03JTJw2k2uvvfaCnr6lTo4bzo7nx969XRdFemEhA0ND\n231RaJpGwbRpVFqtJFRUcGd6Ou8WF7vCa2Xx8dz58cck7d9PSVIS48xmelVVETN58sVPrlD4Gc6c\n33Krlb4tNMJrDXa7nUSrlZhW5BmNRiOJeXmcPXGCFfPnO8JPmZnM0UNgHS21b03XV6czfSL7RfJT\nBmMsLfVwpiNGjKAof7hrpX5+VRo/ngmksTGSoUElzH1vOW/FDWb6wxv9bobixD+t8gHl5eX8z//8\nD398/gWWvrGBP37QwLoD9Uz7p5ls+rfeLJl0kuuvv/6i/5jmFIQtZwZTOiCRddOmuabJR2JiGNbB\n2KZZ07DplSyapnkqDGsap0NDXeGvnFmzOBMa2u7XUiguNe2p3HKnrRVkmqYxJCWFwtTU5hc/dpDm\n7hXelVmapvHwlClcH7KDyfkfcdpk4uHobLSN411/nzp9BplJlfwY1Y8dN99Lv16CRSk53NJ3A4sG\nvMSRA9tZuHAheXl5freaHrqZU7Hb7axdu5bbbruNhIQEfv3rXxMREcHi2yR/WNhI+MqVZGVloWmt\nX6rolHHJdlMQDu6ZyKDyA4zZvBmAL266iYQTJzAYDD57L94Kw4aKCoq8Loa9qakcOXLEZ6+pUFxO\nuFeQtTbPGBcXR4rV6rquAhoaGJyfT+9jx85bFhDQ0kla4GJdX8Fxj3p37ZucPBZGulbAiB8+4t3C\nDOzy3Hmci7FTEw4gNY3kkD0es59hvUqwWCyMGzeOpKQklixZQmFhYRut7Ty6lVPZuXMnDz74IHv3\n7uV3v/sdRUVFfPXVV8y/M56eUUHndtTXpLQG58gh3qkgPHk2YWFhDKg7xJKdT7Pg/ZdZ8j9/JC0y\n0qfxYG+F4Rt27jwvka8WMyquZNqq/WW329m6bRvhNTXMz84m67PP+PnSpTQFBnI8Kso107nYguOW\naI3Gn3uOdmKfXBYallJ1Bqwnoj1PlmVmUZaZmmZmP9azQ/nv//5v3nrrLdLT03n22WdJT0/nuuuu\nY9myZZd8oNmtciq9evVi/fr1jNHDRy6yzCzqwHk1TaNP7+Nc1fs4aBo1R/Yyf8DLrsT98sNPcvPV\nV/s0xqlpGgVTp3LKamVYRQUzbr2VZdu2dWliU6HwZ9qq/WW1WtlbU8MLTz5J8r59DNm9mzOhoa78\nintn1B3tzK+4Sz41l6NpNu8SWkxFbbxnUYBTr3D8RiJi01lx8FFSwwrJrxtMlbyK7/LzmTltGg8+\n+CAVFRWsXbuW1atX8+STT/Lzn/+cW2+9lenTp9N4550ERka26720l27lVAYNGsS4ceM6fB67m6Cc\nM/k3OcvMKeCWvLzzvhQZPfZw5MgR0nwtZ63HnHuaTKQBBSaTK5Hf3sSmQtFdaGsFmXNm0xQUhMVk\nIs5moyoqqvnOqB1J2l8gCuJdkdpgDyC/ZjAJAWexNNfnXtO4+rqbWXXsBOuH34stPp6SxEQWrVzp\nqiqLi4tj8eLFLF68mPz8fP785z/z1ltvMX36dERICL1vv52199/P5MmTiewCB9Ot7kh1dXUdTlzZ\n7XZWrlvHrZ98QpKuwfX62rUkFRYyJi8PYbefNx3N7yqZBE2juB2JzRGodSiK7kdbtb+818ZUxsaS\n1kxI+WgnXsvuXV4/O5rFnw78kkDRSFST1bWupbGx8Zy0lMVCRUUF24eMIG/cOCwmE01BQS2G+QYP\nHswzzzzD/v372bp1K/Hz5lH9zTc8+OCDxMTEcM8997Bu3Tpqamo67T12q5lK9aljvP3Wqg7JQFss\nForLyzkbFcWhhASCS0ooPXSIG6qrKTYaGVRSQjW9ySmbR1JoMbvqh2PzM5kEheJKwVlBtq8Va0q8\nZzZGq5WY+nrm5eRQnJzM8JIS0iIj2Wi3MyYvjyMGg2NNiI/D2lOnz+Cz9X+n9JiRHkH1zO+3wrWu\n5bXyxbz5+gqaTjpkWkZ/mE1J5CCGNjayYfz4Vq+J0zSNUaNGYRw1iuT//V+e/vJL3nnnHdavX8/7\n779Pjx49mDRpEvfdd5/P3pvrtX1+xktI36BKqir3dmiVeX5+PtW6BteGrCzyxo6lOjycTePG0RAU\nhHnsWGzhvTH1OU1FZDw5P53L953Uy8AMqpe8QuEjnDObPH0ZwDfjx/PUwoVUZmYSHxzM7LFjAbh+\n82aCGhoY006Ry9bY8ex972KOm0hGD8/Krt7iEPUnS5jVP5ssPZHfWH2AQUHhPL58OY+sWsUvly9v\nU2GQ0DRuvvlmli1bRllZGZs3b2b27Nl8+eWXPPjggz59b9DNnArgExlodw2u2IoKpKYxVv+ijd28\nmTq7HU1IEhMOUGwydbiXfGdiRoW+FJc33uFbM+3/Tmt6CHnzmDHsM5kIDAzkuMnEAb24Z29NjatE\nOceHUkjuOHX8YsvKKDiT4RFKP1SfTIpXt8jk4F001NdTHxJCWf/+NISEtPq1zHh+VpqmMXr0aF58\n8UUOHz7MF1984bP35XoNn5/xEtNRGehevXqRVlR0ro7dbkdIyeszZnBowACKTCaCz56lye6TtlwK\nheIS4R0FaE974vZgtVqhupqXFixgd/xQltoW8fnRLF49/CRRvePPy9kW1w2h0m7n1TlzyJ04kRfm\nzKGwurrDzk7TNH7yk5/44i150K1yKkcbYjssA61pGlIIlwbX8B072Dl8OPevX+/qBlfXowfb6q/n\n6fErfWi9QqG4lLSnPXF7sNlsFOtVaCunzyDFamXU1q1MTriKn/70p7y79k1yyuaRHFrMnvrhBIfG\nsT0tvUsVnTtCt5qpRPbs0+Fezc7Oj5syM6kPCWHbNdcwND/foxtc9rx5lAvh0AvqIux2O9EWC6Pz\n8kjsYEMwheJywUzXhW+dje2cXSHnuDW28yXuahlS07AmJ9O7roohQ4YQGBjoWkB5PDKeLT+dy/iJ\nkxm6b99ls/C5W81UevTo0eGEubNCpDYvz1H7brUSdfYsuwcP9lAI3mM0dryevbXoGkex1dWcSkri\nBrOZVR1oCKZQKM7HXeQyqqKCbzIzebkT1oIZjUbQq9BKkpMZbtlDuixzOS9N0zDd/x5z9f1NdjsZ\n27dfNgufu5VT8QWapvHQvffyjw8/ZHB+Pjf260fqqFFUf/LJeQrBn3aRQnC0rnGU7dVcqDO75CkU\nlwu+rJB0liifMpnoSeeEcrzVMuZmTWpVF9fnrVbiLoOFz8qpeNHY2MhzL71EXHAwBRkZNFks5JeW\ncjo83KNj3OOvvEJXpeojLkGXPIVC0Ym4qWW05grW2tnj5VLgn67uErJlyxYqg4NZMX8+uRMnsmL+\nfGoBi8l0nkJwny4Sbqu5BF3yFApF52Gm+65BU07Fi4MHD1Lk5UD2DxzYrEJwZ8o5uHNcTyDO6+QE\nokJxpRBA8wUAdrudRIuFMX5YEBNA2+X4LwXKqXgxYMCA8xxI7JEjIISrKmSxflPf3wk39RE0M4LR\nY6qVmZk0BgfzTRv6ZisUitbh3vQrqKGBmzppRX13p9NyKkKIlcBk4IiUcoi+bTjwChAB7AemSymr\nmjn2NuAFHI45R0r5bGfZ6c3o0aPZumMHC1asoCg1lbSiImIaG1k6bx599u/H4JYoW61p7OgiuzRN\n47jJxA6TiRGo0YBC0R7sdjspXgrkTkUM96ZfztzpYlUQ02Y68970BnCb17Yc4LdSyqHA+8CvvA8S\nQgQALwG3AxnAA0KIjE6004PAwECeWriQ/MGDibfZuGfwYJ5auJCm4GAsJhNb29n6VKFQXFqcM5Ex\nZjMBbu2H0WciXbWivrvTaTMVKeVmIcQgr82pwGb98efAp8C/eu1zPWCVUpYCCCHWAXcBBZ1lqzeB\ngYGY9VW0b3bViyoUik7FORPJ8ZqJRFutHDeZumxFfXdHSCkvvld7T+5wKh+5hb++BJ6TUv5VCPFz\n4N+klJFex9wH3CalnK0/fwi4QUr5eAuvMRdc64SGAHs64734kD7AsUttRCtQdvoWZadvabOdUVFR\nBhkbG18dFeXaFllVhaisLK+qqrIB9O7b1xQQEhJeFxqq9Thzxt509mztj0ePWrrSzktAqvd9uCN0\n9TqVmcAyIcT/Az4E6pvZp7nlHy16Pinlq8CrAEKI76WU1/rC0M7icrARlJ2+RtnpW5SdvkMI8b0v\nz9elTkVKWQjcAiCEMAF3NLNbGZDg9rw/UN751ikUCoWio3RptlkIEaP/1oB/wVEJ5s13QIoQIlEI\nEQxMwzGrUSgUCoWf02lORQixFvgKSBVClAkhZuGo5LIAhThmH6v0feOFEH8HkFI2Ao/jSOLvBd6R\nUua38mVf9fHb6AwuBxtB2elrlJ2+RdnpO3xqY6cm6hUKhUJxZaEWWygUCoXCZyinolAoFAqf4ddO\nRQiRIITYJITYK4TIF0I8qW8fIYT4WgixQwjxvRDien17tL5/jRBiude5rhFC7BZCWIUQy4QQPlGu\nb4eNE4UQP+i2/CCEGN/ZNrbTzuv1bTuEEDuFEHf7o51uxw3Q/++/9Ec7hRCDhBBn3D7TV9zO5Td2\n6n8bJoT4St9/txCih7/ZKYSY7vZZ7hBC2IUQI/zQziAhxJu6PXuFEE+5nctf7knBQohVui07hRDj\nOmSjlNJvfwADcLX+OBKw4JBu+Qy4Xd8+CTDrj8OBm4H5wHKvc30L3IRjHcw/nMdfAhtHAvH64yHA\n4c62sZ12hgGBbscecXvuN3a6Hfce8C7wSz/9PAcBe1o4lz/ZGQjsAobrz6OBAH+z0+vYoUCpn36e\nDwLr3K6p/cCgzrSzHTY+BqzSH8cAPwBae23065mKlNImpdymP67GUQ3WD8diSOey2J7o61iklLVS\nyq1Anft5hBAGIEpK+ZV0fFL/B/zsEtm4XUrpXHeTD/QQQoR0po3ttPO0dFTiAfTQ9+vUz7I9duo2\n/QwoxfF5Orf5nZ3N4Yd23gLsklLu1I85LqVs8kM73XkAWAt++XlKIFwIEQiE4ljwXeVP9yQcDmeD\nvv8R4CRwbbtt9JUH7+wfHCO9g/qHkq4/PgQcBgZ67fsIbjMV4Fog1+35aBzyMZfMRn3/+5x2dZWN\nbbETuAHHjboGuNsf7cQxO/0Kh/L1H9BnKn5o5yCgFtgO5AGj/dTORcBqHCX924Bf+6OdXvuXAEP8\n0U4gCFgHHNX//3O70s5W2jgXxyw/EEjE4VTuba+Nfj1TcSKEiMAR3lgkHVL5C4DFUsoEYDHw+sVO\n0cw2n9ZSt9VGIcRg4DlgXlfZ2FY7pZTfSCkHA9cBT+mxdX+z89+A/5VS1nifws/stAEDpJQjgZ8D\na4QQUX5oZyCOEPJ0/ffdQogJfminc/8bgNNSSqfmn7/ZeT3QBMTjuGH/QgiR1BV2tsHGlTiUTL4H\nlgJfAo3ttrEzPLiPPW0QjlHTz922neLcGhsBVHkd8wieMxUDUOj2/AEg+1LZiEN6xgKM6iob2/tZ\nuu23CcfIxa/sBLbgiFPvxzHC+hHH4lm/srOZY81++nlOA95w2+9fcbSo8Cs73f7+v8Dv3J77lZ04\n2ng85LbfSuD+zrazg9/NL3GExNplo1/PVPRKg9eBvVLKP7n9qRwYqz8eDxRf6DxSShtQLYS4UT/n\nPwN/vRQ2CiGuAj4GnpJSftEVNrbTzkQ9DowQYiCOtgX7/c1OKeVoKeUgKeUgHKOsP0opl/ubnUKI\nvsLRKwh9pJqCI7nsV3biuBENE0KE6f//sUCBH9rplHuagiO8BPjfdYQj3DReOAgHbsRxo/ane1KY\nbhtCiIlAo5Sy/f9zX3pwX//gmH5LHNUoO/SfSfr2H4CdwDfANW7H7McxWq3BMaXL0Ldfi0MWvwRY\nju6xu9pGHJpntW777gBiOtPGdtr5EI58yg4csfWfuZ3Lb+z0OvYPeFZ/+Y2dOGLU+fr2bcCd/min\nfsw/6bbuAZ73YzvHAV83cy6/sRNHru9d/fMsAH7V2Xa2w8ZBQBGOhH4unnnVNtuoZFoUCoVC4TP8\nOvylUCgUissL5VQUCoVC4TOUU1EoFAqFz1BORaFQKBQ+QzkVhUKhUPgM5VQU3QLhUKh2KtZWCCEO\n649rhBAvd7EtbwkhioQQe4QQK4UQQV5/v04I0SSEuM9t235dDXaHEOJ7t+29hRCfCyGK9d+93P72\nlK4eWySEuNVte7PKsrrG3Nv69m+EEIM683NQXJkop6LoFkiH8OEIKeUI4BUc0i0jpJQRUspHu9ic\nt4A0HOq5ocBs5x/0BZDP4Vhk6E2mbvO1btt+C2yQUqbgEP37rX6eDByr3wcDtwEvOxdXAitw6Dml\n6D+36dtnASeklEYcK9Gf6/hbVSg8UU5F0a0RQowTQnykP/6DcPS2+EyfGdwjhHheH9V/4pxR6CP9\nPOHod/OpcKi1thop5d+lDg7p8P5uf16IQ4/pSCtPdxfwpv74Tc6pxN6FQ1L9rJRyH2AFrhcXVpZ1\nP9d6YIJzFqNQ+ArlVBRXGsnAHThusH8GNkkphwJngDt0x/IicJ+U8hocWk1Pt+eF9HM9BHyiP+8H\n3I1jJuWNBD7THdlct+2x0iGXgf47Rt/eD4farJMyfVs//bH3do9jpKOtwSkc/VIUCp8ReKkNUCi6\nmH9IKRuEELuBAPQbPrAbh1xFKo7maZ/rg/gAHArD7eFlYLOUcov+fCnwG+noT+K97ygpZbkQIkZ/\n7UIp5eYLnLslBdkLKct2iYKv4spGORXFlcZZACmlXQjRIM/pFNlxXA8CyJdS3tTSCfTcxQ/60w+l\nlP+vmX1+D/TlXGsDcOgordMdSh9gkhCiUUr5gdQbt0kpjwgh3schmb4ZqBRCGKSUNj205QyblQEJ\nbufuj0MwsAzPcJtzu/sxZbpYZE8cOnkKhc9Q4S+FwpMioK8Q4iZw9Rgf7L6DlLLJWRTQgkOZDdwK\nPCCltLsdlyjPKSqvBx6VUn4ghAgXQkTqx4bj6L7o7A/yIfCw/vhhzqnEfghM0yu6EnEk5L+VF1aW\ndT/XfcBGN6eqUPgENVNRKNyQUtbrpb7LhBA9cVwjS3FrVdwKXgEOAF/ps5K/SCn//QL7xwLv6/sG\nAmuklM6w3LPAO0KIWThk1KfoduYLId7BoXzbCDwmpWzSj1kAvIGj8uwf+g845NBXCyGsOGYo09rw\nnhSKVqFUihUKhULhM1T4S6FQKBQ+QzkVhUKhUPgM5VQUCoVC4TOUU1EoFAqFz1BORaFQKBQ+QzkV\nhUKhUPgM5VQUCoVC4TP+Pzy+NuDT0qgCAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Now suppose you wanted to combine the two together:\n", "my_event = MulensModel.Event(\n", " datasets=[MOA_data, OGLE_data], model=my_1S2L_model)\n", "\n", "# And you wanted to plot the result:\n", "my_event.plot_model(\n", " t_range=[2452810,2452890], subtract_2450000=True, color='black',\n", " data_ref=1)\n", "my_event.plot_data(\n", " subtract_2450000=True, data_ref=1, s=5, markeredgecolor='gray')\n", "# MulensModel automatically fits for the source and blend flux for the \n", "# given model.\n", "\n", "# Customize the output\n", "plt.legend(loc='best')\n", "plt.title('OGLE-2003-BLG-235/MOA-2003-BLG-53')\n", "plt.ylim(19., 16.5)\n", "plt.xlim(2810,2890)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Chi2 of the fit: 1774.42\n" ] } ], "source": [ "# If you want to see how good the fit is, output the chi2:\n", "print('Chi2 of the fit: {0:8.2f}'.format(my_event.get_chi2()))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If you want to optimize the chi2, we leave it up to you to determine the best method for doing this." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.14" } }, "nbformat": 4, "nbformat_minor": 2 }