From e219e57290937ac5f1de8e9cf36ab3a70a8f82be Mon Sep 17 00:00:00 2001 From: David A Minton Date: Thu, 23 Sep 2021 10:15:25 -0400 Subject: [PATCH] Fixed initial conditions generator and Notebook for Helio Swifter comparison --- .../helio_swifter_comparison/init_cond.py | 10 +-- .../param.swiftest.in | 2 +- .../swiftest_vs_swifter.ipynb | 62 ++++++++----------- .../helio_swifter_comparison/tp.swifter.in | 8 +-- 4 files changed, 36 insertions(+), 46 deletions(-) diff --git a/examples/helio_swifter_comparison/init_cond.py b/examples/helio_swifter_comparison/init_cond.py index b8b9b4369..c3f4be742 100755 --- a/examples/helio_swifter_comparison/init_cond.py +++ b/examples/helio_swifter_comparison/init_cond.py @@ -32,10 +32,10 @@ "Saturn": 6, "Uranus": 7, "Neptune": 8, - "Ceres": 101, - "Pallas": 102, - "Juno": 103, - "Vesta": 104 + "Ceres": 9, + "Pallas": 10, + "Juno": 11, + "Vesta": 12 } for name, id in bodyid.items(): @@ -44,7 +44,7 @@ sim.param['PL_IN'] = "pl.swiftest.in" sim.param['TP_IN'] = "tp.swiftest.in" sim.param['CB_IN'] = "cb.swiftest.in" -sim.param['BIN_OUT'] = "bin.swiftest.dat" +sim.param['BIN_OUT'] = "bin.swiftest.nc" sim.param['ENC_OUT'] = "enc.swiftest.dat" sim.save("param.swiftest.in") sim.param['PL_IN'] = "pl.swifter.in" diff --git a/examples/helio_swifter_comparison/param.swiftest.in b/examples/helio_swifter_comparison/param.swiftest.in index 520d03553..3d20ed4e3 100644 --- a/examples/helio_swifter_comparison/param.swiftest.in +++ b/examples/helio_swifter_comparison/param.swiftest.in @@ -11,7 +11,7 @@ IN_TYPE ASCII PL_IN pl.swiftest.in TP_IN tp.swiftest.in CB_IN cb.swiftest.in -BIN_OUT bin.swiftest.dat +BIN_OUT bin.swiftest.nc CHK_QMIN 0.004650467260962157 CHK_RMIN 0.004650467260962157 CHK_RMAX 1000.0 diff --git a/examples/helio_swifter_comparison/swiftest_vs_swifter.ipynb b/examples/helio_swifter_comparison/swiftest_vs_swifter.ipynb index 709b6cd44..26e1b3f79 100644 --- a/examples/helio_swifter_comparison/swiftest_vs_swifter.ipynb +++ b/examples/helio_swifter_comparison/swiftest_vs_swifter.ipynb @@ -43,8 +43,8 @@ "output_type": "stream", "text": [ "Reading Swiftest file param.swiftest.in\n", - "Reading in time 1.000e+00\n", - "Creating Dataset\n", + "\n", + "Creating Dataset from NetCDF file\n", "Successfully converted 1462 output frames.\n", "Swiftest simulation data stored as xarray DataSet .ds\n" ] @@ -79,8 +79,8 @@ "metadata": {}, "outputs": [], "source": [ - "swiftdiff['dr'] = np.sqrt(swiftdiff['px']**2 + swiftdiff['py']**2 + swiftdiff['pz']**2)\n", - "swiftdiff['dv'] = np.sqrt(swiftdiff['vx']**2 + swiftdiff['vy']**2 + swiftdiff['vz']**2)" + "swiftdiff['dr'] = np.sqrt(swiftdiff['xhx']**2 + swiftdiff['xhy']**2 + swiftdiff['xhz']**2)\n", + "swiftdiff['dv'] = np.sqrt(swiftdiff['vhx']**2 + swiftdiff['vhy']**2 + swiftdiff['vhz']**2)" ] }, { @@ -100,7 +100,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAElCAYAAADDUxRwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAiXUlEQVR4nO3dfbxVZZ338c+3A4QKigrIwwFBQwFREQjUzNQGBswyFBtRS80iTZsa9VaqmUznnrSa8iEthsyn7JYcK0XDJ0BvC2UUBVQkEhHjCAiixJMEHH73H2vhvT3tc84+e6/zsNnf9+t1Xuy91rWu9Vug+3uutde6liICMzOzUn2otQswM7PdgwPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDHLQ9J3Jd2dvu4raZOkqiL6mSLp37Kv0KztcaDYbknSckn/UGfZeZL+2NS+IuIvEdEpImqL2PbCiPj3QtpKukPS/27qPrJS7N+P2S4OFLMKIKnd7rAPa9scKFaxJPWS9BtJayW9Lumf62nXT1Ls+sBMt5su6R1JSyV9uYF9vD/qkHSCpBpJl0laI2mVpPPTdZOAs4Er0tNrDzZWo6Q9JN0p6V1JiyVdIakmZ/1ySVdKehHYLKmdpMmSXpO0UdIrksanbQcBU4Bj0v2vT5fvI+mudP9vSPpXSR9K150naY6k6yW9A3y32H8L2z34NwqrSOmH4oPAA8BEoBqYKWlJRDzayOb3AIuAXsBA4HFJyyJiVgG77gHsA/QGRgP3Sbo/IqZKOhaoiYh/LbDGq4B+wEHAXsCMPPubCHwKeDsidkh6Dfg4sBo4A7hb0kciYrGkC4EvRcRxOdv/JK33IGB/4DFgFfCLdP0oYBrQHWhfwPHbbswjFNud3S9p/a4f4Kc56z4KdIuIayJiW0QsA34OnNlQh5L6AMcBV0bE1ohYANwKfL7AmrYD10TE9oiYAWwCDq2nbWM1fg74XkS8GxE1wE15+rgpIlZExHsAEfHfEbEyInZGxK+BV4GR9RxrFfBPwDcjYmNELAd+VOdYV0bETyJix659WOXyCMV2Z5+NiJm73kg6D/hS+vZAoNeuUzupKuAPjfTZC3gnIjbmLHsDGFFgTesiYkfO+y1Ap3raNlZjL2BFzrrc13mXSfoCcCnJyIZ0313r2X9XoAPJ8e3yBsnoqqF9WoVyoFilWgG8HhEDmrjdSmA/SZ1zQqUv8GYGNdWd+ruxGleRnAZ7JX3fp6E+JR1IMsL5JPBMRNRKWgConv2/TTKiOjBnH3WP1dOV2/t8yssq1bPAhvRL6z0kVUkaIumjDW0UESuAp4FrJXWUdARwAfCrDGp6i+S7ikJrvBf4pqR9JfUGLmmk/71IAmAtQHpBwJA6+6+W1AEgvUz6XuA/JHVOA+lS4O7SDtN2Vw4Uq0jph+WngaHA6yS/jd9K8gV0YyaSnDJaCfwOuCoiHs+grF8Ag9PvfO4voMZrgJp03UzgPuBv9XUeEa+QfAfyDEl4HA7MyWkym+Rig9WS3k6XfQ3YDCwD/gj8H+C2Ug/Udk/yA7bMdg+SLgLOjIhPtHYtVpk8QjErU5J6SvqYpA9JOhS4jGTEZNYq/KW8WfnqAPwX0B9YT3I/yE8b2sCsOfmUl5mZZcKnvMzMLBMOFLMmyDeL8e6i7pxlZk3lQDGrI/1Q3ZxOkvimpB+riGehZFDDR1pyn2alcqCY5XdkRHQiuav8LKDeGYXNLOFAMWtARPyJZO6sIXXXSRop6Zn0RsRVkm7edZd5uj4kXSjp1XSK+VskKWf9F9Np59+V9Gh6JzqSnkqbLExHSf8kqaukh9J9vSPpD7umkc9T17GSnpP01/TPY3PWPSnp39Np5zdKekzS383lJekMSc/XWXaZpPub9jdolcSBYtYASYNJpnufn2d1LfAvJJMoHkMymvlqnTankMwafCTJ7MD/mPb7WeBbwGlAN5LQugcgIo5Ptz0yfVLkr0nuMalJ2x6Qbvt3l2hK2g/4PcnMw/sDPwZ+L2n/nGZnAeeTTDnfAbg8z7FNB/qnz0nZ5Rzgl3namgEOFLP6vCDpXZLnkdwK3F63QUQ8HxFz06nbl5PcE1L3LvXrImJ9RPwFeIJkGhWArwDXRsTidPbh7wFDd41S8tgO9AQOTKe+/0Pkv+b/U8CrEfHLtK57gD+RTOGyy+0R8ed0uvl7c2rKPba/Ab8mCREkHUYy3cxD9dRn5kAxq8ewiNg3Ig6OiH+NiJ11G0g6JD0NtVrSBpJQqHv6aHXO69yp6g8Ebsx5Vss7JLP+9ia/HwJLgcckLZM0uZ52vfjgdPPw91PO11dTXXcCZ6Wn6T4P3JsGjVleDhSz4v2M5Lf/ARGxN8lpKDW8yftWAF+JiC45P3tExNP5GqcPuLosIg4iGW1cKumTeZquJAmrXEVNrx8Rc4FtJKf8zsKnu6wRDhSz4nUGNgCbJA0ELmrCtlNIpp4/DN5/dvsZOes/MJW9pFMkfSQdLWwg+f6mNk+/M4BDJJ2l5Bny/wQMpvhTVXcBNwM7IuKPRfZhFcKBYla8y0l+c99I8uCqXxe6YUT8Dvg+MC09XfYyMC6nyXeBO9NTYp8DBpBMUb+JZPr5n0bEk3n6XUdyIcBlwDrgCuCUiHi7btsC/ZLkCjePTqxRnsvLzOolaQ9gDcl3Sq+2dj3WtnmEYmYNuQh4zmFihfCcPWaWl6TlJBcZfLZ1K7Fy4VNeZmaWCZ/yMjOzTFT0Ka+uXbtGv379WrsMM7Oy8vzzz78dEd3qLq/oQOnXrx/z5s1r7TLMzMqKpLqzMQA+5WVmZhlxoJiZWSYcKGZmlomK/g7FzKw1bN++nZqaGrZu3drapTSoY8eOVFdX0759+4LaO1DMzFpYTU0NnTt3pl+/fuQ8xLNNiQjWrVtHTU0N/fv3L2gbn/IyM2thW7duZf/992+zYQIgif33379JoygHiplZK2jLYbJLU2t0oJiZWSYcKGZmZerYY4/Nu/y8887jvvvua+FqHChmZmXr6afzPjG61fgqLzOzMtWpUyc2bdpERPC1r32N2bNn079/f1prFnmPUMzMytzvfvc7lixZwksvvcTPf/7zVhu5OFDMzMrcU089xcSJE6mqqqJXr16cdNJJrVKHA8XMbDfQFi5DdqCYmZW5448/nmnTplFbW8uqVat44oknWqUOfylvZlbmxo8fz+zZszn88MM55JBD+MQnPtEqdThQzMzK1KZNm4DkdNfNN9/cytX4lJeZmWXEgWJmZplwoJiZWSYcKGZmlgkHipmZZcKBYmZmmXCgmJlVqC9+8Yt0796dIUOGZNKfA8XMrEKdd955PPLII5n116YCRdJYSUskLZU0Oc96SbopXf+ipGF11ldJmi/poZar2sysPB1//PHst99+mfXXZu6Ul1QF3AKMBmqA5yRNj4hXcpqNAwakP6OAn6V/7vJ1YDGwd4sUbWZWoqsfXMQrKzdk2ufgXntz1acPy7TPQrSlEcpIYGlELIuIbcA04NQ6bU4F7orEXKCLpJ4AkqqBTwG3tmTRZmaWaDMjFKA3sCLnfQ0fHH3U16Y3sAq4AbgC6NzQTiRNAiYB9O3bt6SCzcxK1RojiebSlkYo+Sbzr/scy7xtJJ0CrImI5xvbSURMjYgRETGiW7duxdRpZmZ5tKVAqQH65LyvBlYW2OZjwGckLSc5VXaSpLubr1Qzs/I3ceJEjjnmGJYsWUJ1dTW/+MUvSuqvLZ3yeg4YIKk/8CZwJnBWnTbTgUskTSM5HfbXiFgFfDP9QdIJwOURcU4L1W1mVpbuueeeTPtrM4ESETskXQI8ClQBt0XEIkkXpuunADOAk4GlwBbg/Naq18zMPqjNBApARMwgCY3cZVNyXgdwcSN9PAk82QzlmZlZA9rSdyhmZlbGHChmZpYJB4qZmWXCgWJmZplwoJiZVaAVK1Zw4oknMmjQIA477DBuvPHGkvtsU1d5mZlZy2jXrh0/+tGPGDZsGBs3bmT48OGMHj2awYMHF92nRyhmZhWoZ8+eDBuWPAGkc+fODBo0iDfffLOkPj1CMTNrTQ9PhtUvZdtnj8Nh3HUFN1++fDnz589n1Ki68/E2jUcoZmYVbNOmTZx++unccMMN7L13aY+S8gjFzKw1NWEkkbXt27dz+umnc/bZZ3PaaaeV3J9HKGZmFSgiuOCCCxg0aBCXXnppJn06UMzMKtCcOXP45S9/yezZsxk6dChDhw5lxowZjW/YAJ/yMjOrQMcddxzJfLvZ8QjFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMKtDWrVsZOXIkRx55JIcddhhXXXVVyX36PhQzswr04Q9/mNmzZ9OpUye2b9/Occcdx7hx4zj66KOL7tMjFDOzCiSJTp06AcmcXtu3b0dSSX16hGJm1oq+/+z3+dM7f8q0z4H7DeTKkVc22q62tpbhw4ezdOlSLr74Yk9fb2ZmxamqqmLBggXU1NTw7LPP8vLLL5fUn0coZmatqJCRRHPr0qULJ5xwAo888ghDhgwpuh+PUMzMKtDatWtZv349AO+99x4zZ85k4MCBJfXpEYqZWQVatWoV5557LrW1tezcuZPPfe5znHLKKSX16UAxM6tARxxxBPPnz8+0T5/yMjOzTDhQzMwsE20qUCSNlbRE0lJJk/Osl6Sb0vUvShqWLu8j6QlJiyUtkvT1lq/ezKyytZlAkVQF3AKMAwYDEyUNrtNsHDAg/ZkE/CxdvgO4LCIGAUcDF+fZ1szMmlGbCRRgJLA0IpZFxDZgGnBqnTanAndFYi7QRVLPiFgVES8ARMRGYDHQuyWLNzOrdG0pUHoDK3Le1/D3odBoG0n9gKOA/8m+RDMzq09bCpR8s5JFU9pI6gT8BvhGRGzIuxNpkqR5kuatXbu26GLNzHYHtbW1HHXUUSXfgwIF3IciqW+Bfa2v70O8QDVAn5z31cDKQttIak8SJr+KiN/Wt5OImApMBRgxYkTdwDIzqyg33ngjgwYNYsOGUj6+E4Xc2HgnySigoXmNA7gDuKuEWp4DBkjqD7wJnAmcVafNdOASSdOAUcBfI2KVkjmXfwEsjogfl1CDmVnFqKmp4fe//z3f/va3+fGPS//obDRQIuLEussk9YiI1SXv/YP72SHpEuBRoAq4LSIWSbowXT8FmAGcDCwFtgDnp5t/DPg88JKkBemyb0XEjCxrNDPL2urvfY+/Lc52+voPDxpIj299q9F23/jGN/jBD37Axo0bM9lvsVOvfAH4QSYV5EgDYEadZVNyXgdwcZ7t/kjDIygzM8vx0EMP0b17d4YPH86TTz6ZSZ/FBsqpkrYAj0fEkkwqMTOrQIWMJJrDnDlzmD59OjNmzGDr1q1s2LCBc845h7vvvrvoPou9yus0ktNO4yXdWvTezcysVVx77bXU1NSwfPlypk2bxkknnVRSmECRI5SIeAt4JP0xMzMrboQi6RZJd6Svx2RakZmZtagTTjiBhx56qOR+ij3ltQ1Ylr4+qeQqzMys7BUbKFuAfdKbCQu98dHMzHZjxV7l9Q7wHsnswHOyK8fMzMpVk0YokrpIuh04PV10FzAi86rMzKzsNGmEEhHrJV0H9APeBo4A6p03y8zMKkcxp7wuAF6PiEeB5zOux8zMylQxgfIucKGkQ4GFwIKImJ9tWWZm1tz69etH586dqaqqol27dsybN6+k/pocKBFxraRZwJ+BocDxgAPFzKwMPfHEE3Tt2jWTvpocKJKuIZkNeAHJ6OTJTCoxM7OyVswI5TuSDiB5zO7pkg6OiC9nX5qZ2e7vD/f+mbdXbMq0z659OvHxzx3SaDtJjBkzBkl85StfYdKkSSXtt9j7UL4C/FdEeC4vM7MyNWfOHHr16sWaNWsYPXo0AwcO5Pjjjy+6v2ID5TbgIkl7kTxyd0HRFZiZVbBCRhLNpVevXgB0796d8ePH8+yzz5YUKMVOvfLPJGHUDrip6L2bmVmr2Lx58/tPaty8eTOPPfYYQ4YMKanPYkcorwEDgAci4l9KqsDMzFrcW2+9xfjx4wHYsWMHZ511FmPHji2pz2IDZRGwArhA0g8j4qMlVWFmZi3qoIMOYuHChZn2WWygHAKsBaaS3OhoZmYVrtjvUAaS3Mx4OVDadWZmZrZbKDZQugBXAlcAWzOrxszMylaxp7yuAQZGxBJJO7MsyMzMylNBIxRJVZJWSfoSQETURMTM9PXk5izQzMzKQ0GBEhG1wMvAwc1bjpmZlaumfIeyJ3CFpHmSpqc/DzRXYWZm1rzWr1/PhAkTGDhwIIMGDeKZZ54pqb+mfIdyTPrnsPQHIErau5mZtZqvf/3rjB07lvvuu49t27axZcuWkvprSqD0L2lPZmbWZmzYsIGnnnqKO+64A4AOHTrQoUOHkvosOFAi4o2S9mRmZn/niTumsuaNZZn22f3AgzjxvIZvEVy2bBndunXj/PPPZ+HChQwfPpwbb7yRvfbaq+j9FnsfipmZlbEdO3bwwgsvcNFFFzF//nz22msvrrvuupL6LPY+FDMzy0BjI4nmUl1dTXV1NaNGjQJgwoQJJQdKk0cokj5d0h4b7nuspCWSlkr6u/tblLgpXf+ipGGFbmtmZv9fjx496NOnD0uWLAFg1qxZDB48uKQ+ixmh/AfwYEl7zUNSFXALMBqoAZ6TND0iXslpNo5k2vwBwCjgZ8CoArc1M7McP/nJTzj77LPZtm0bBx10ELfffntJ/RUTKCppj/UbCSyNiGUAkqYBpwK5oXAqcFdEBDBXUhdJPYF+BWybmTsu+x7vdWjfHF2bWQUY/umPs6ZmdavW0L5KDB06lHnz5mXWZzGB0lz3nvQmecbKLjUko5DG2vQucFsAJE0inSG5b9++RRW6U1W81662qG3NzEKwU617G1/szH7/belL+Xwjn7pHXF+bQrZNFkZMJXmOCyNGjCjqb/SL/3llMZuZmQGwePFievTu2dplZK4tBUoN0CfnfTWwssA2HQrY1szMmlEx96G8lXkVieeAAZL6S+oAnAlMr9NmOvCF9Gqvo4G/RsSqArc1M7Nm1OQRSkSMbo5CImKHpEuAR4Eq4LaIWCTpwnT9FGAGcDKwFNgCnN/Qts1Rp5mZ5deWTnkRETNIQiN32ZSc1wFcXOi2ZmbWcjz1iplZBVqyZAlDhw59/2fvvffmhhtuKKnPokYoki6NiB+nrw+NiCUlVWFmZi3q0EMPZcGCBQDU1tbSu3dvxo8fX1KfTQoUSV2A64GBkrYCLwIXkH6XYWZm5WfWrFkcfPDBHHjggSX106RAiYj1wPmSPgWsBsYAvy2pAjOzCrb+wdfYtnJzpn126LUXXT5d+BPbp02bxsSJE0veb7HfoXyC5PLho0nmzzIzszK0bds2pk+fzhlnnFFyX8Ve5dUFuBK4guSUl5mZFaEpI4nm8PDDDzNs2DAOOOCAkvsqNlCuAQZGxBJJO0uuwszMWsU999yTyekuKPKUV0TURMTM9LWfPWJmVoa2bNnC448/zmmnnZZJf0UFiqRbJN2Rvh6TSSVmZtai9txzT9atW8c+++yTSX/Ffim/DViWvj4pk0rMzKysFRsoW4B9JLUHinuoiJmZ7VaK/VL+HeA9ksfuzsmuHDMzK1dNGqGkj9y9HTg9XXQXMCLzqszMrOw0+U55SdeRPMP9beAIfKe8mZlR3CmvC4DXI+JR4PmM6zEzszJVzJfy7wIXSrpB0vmSjsq6KDMza37XX389hx12GEOGDGHixIls3bq1pP6aHCgRcS3wZeC7wOvA8SVVYGZmLe7NN9/kpptuYt68ebz88svU1tYybdq0kvps8ikvSdeQPGZ3AbAgIp4sqQIzM2sVO3bs4L333qN9+/Zs2bKFXr16ldRfMc+U/46k75CMbk6XdHBEfLmkKszMKtTDDz/M6tWrM+2zR48ejBs3rsE2vXv35vLLL6dv377ssccejBkzhjFjSpv4pNgbG28DBgH7Az8tqQIzM2tx7777Lg888ACvv/46K1euZPPmzdx9990l9VnsjY3/TDL9SjvgRvw9iplZURobSTSXmTNn0r9/f7p16wbAaaedxtNPP80555xTdJ/FjlBeAzoCD0SEw8TMrMz07duXuXPnsmXLFiKCWbNmMWjQoJL6LDZQFgGzgQskPVdSBWZm1uJGjRrFhAkTGDZsGIcffjg7d+5k0qRJJfVZ7Cmvg0nuR5ma/mlmZmXm6quv5uqrr86sv2IDZUVEzJbUE1iTWTVmZla2ij3lNVZSNTAFuD7DeszMrEwVGyhdgCuBK4C/ZVaNmVmFiIjWLqFRTa2x2EC5huQKryVAbZF9mJlVpI4dO7Ju3bo2HSoRwbp16+jYsWPB2xT0HYqkKqAG+LeIuDUiatL3RMTkYoo1M6tU1dXV1NTUsHbt2tYupUEdO3akurq64PYFBUpE1Ep6meTqLjMzK0H79u3p379/a5eRuaac8toTuELSPEnT058HsihC0n6SHpf0avrnvvW0GytpiaSlkibnLP+hpD9JelHS7yR1yaIuMzMrXFMC5RhAwDDglJyfLEwGZkXEAGBW+v4D0tNutwDjgMHAREmD09WPA0Mi4gjgz8A3M6rLzMwK1JT7UJpzfHYqcEL6+k7gSZKryHKNBJZGxDIASdPS7V6JiMdy2s0FJjRjrWZmlkejgSKpb/oy7+UIOevXR8SGIus4ICJWAUTEKknd87TpDazIeV8DjMrT7ovAr4usw8zMilTICOVOkjBRA20CuAO4q74GkmYCPfKs+nYBNVDP/j8QcpK+DewAftVAHZOASZBMjmZmZtloNFAi4sQsdhQR/1DfOklvSeqZjk7qm86lBuiT874aWJnTx7kk3+l8Mhq4uDsippLMQcaIESPa7kXgZmZlptgbG7M2HTg3fX0ukO/qseeAAZL6S+oAnJluh6SxJN+5fCYitrRAvWZmVkdbCZTrgNGSXgVGp++R1EvSDICI2AFcAjwKLAbujYhF6fY3A52BxyUtkDSlpQ/AzKzSFTvbcKYiYh3wyTzLVwIn57yfAczI0+4jzVqgmZk1qq2MUMzMrMw5UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLRJsIFEn7SXpc0qvpn/vW026spCWSlkqanGf95ZJCUtfmr9rMzHK1iUABJgOzImIAMCt9/wGSqoBbgHHAYGCipME56/sAo4G/tEjFZmb2AW0lUE4F7kxf3wl8Nk+bkcDSiFgWEduAael2u1wPXAFEM9ZpZmb1aCuBckBErAJI/+yep01vYEXO+5p0GZI+A7wZEQsb25GkSZLmSZq3du3a0is3MzMA2rXUjiTNBHrkWfXtQrvIsywk7Zn2MaaQTiJiKjAVYMSIER7NmJllpMUCJSL+ob51kt6S1DMiVknqCazJ06wG6JPzvhpYCRwM9AcWStq1/AVJIyNidWYHYGZmDWorp7ymA+emr88FHsjT5jlggKT+kjoAZwLTI+KliOgeEf0ioh9J8AxzmJiZtay2EijXAaMlvUpypdZ1AJJ6SZoBEBE7gEuAR4HFwL0RsaiV6jUzszpa7JRXQyJiHfDJPMtXAifnvJ8BzGikr35Z12dmZo1rKyMUMzMrcw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTCgiWruGViNpLfBGkZt3Bd7OsJxy4GOuDD7mylDKMR8YEd3qLqzoQCmFpHkRMaK162hJPubK4GOuDM1xzD7lZWZmmXCgmJlZJhwoxZva2gW0Ah9zZfAxV4bMj9nfoZiZWSY8QjEzs0w4UMzMLBMOlEZIGitpiaSlkibnWS9JN6XrX5Q0rDXqzFIBx3x2eqwvSnpa0pGtUWeWGjvmnHYflVQraUJL1pe1Qo5X0gmSFkhaJOn/tnSNWSvgv+t9JD0oaWF6zOe3Rp1ZknSbpDWSXq5nfbafXxHhn3p+gCrgNeAgoAOwEBhcp83JwMOAgKOB/2ntulvgmI8F9k1fj6uEY85pNxuYAUxo7bqb+d+4C/AK0Dd93721626BY/4W8P30dTfgHaBDa9de4nEfDwwDXq5nfaafXx6hNGwksDQilkXENmAacGqdNqcCd0ViLtBFUs+WLjRDjR5zRDwdEe+mb+cC1S1cY9YK+XcG+BrwG2BNSxbXDAo53rOA30bEXwAiohKOOYDOkgR0IgmUHS1bZrYi4imS46hPpp9fDpSG9QZW5LyvSZc1tU05aerxXEDyG045a/SYJfUGxgNTWrCu5lLIv/EhwL6SnpT0vKQvtFh1zaOQY74ZGASsBF4Cvh4RO1umvFaT6edXu5LL2b0pz7K611kX0qacFHw8kk4kCZTjmrWi5lfIMd8AXBkRtckvsGWtkONtBwwHPgnsATwjaW5E/Lm5i2smhRzzPwILgJOAg4HHJf0hIjY0c22tKdPPLwdKw2qAPjnvq0l+e2lqm3JS0PFIOgK4FRgXEetaqLbmUsgxjwCmpWHSFThZ0o6IuL9FKsxWof9dvx0Rm4HNkp4CjgTKNVAKOebzgesi+XJhqaTXgYHAsy1TYqvI9PPLp7wa9hwwQFJ/SR2AM4HpddpMB76QXi1xNPDXiFjV0oVmqNFjltQX+C3w+TL+jTVXo8ccEf0jol9E9APuA75apmEChf13/QDwcUntJO0JjAIWt3CdWSrkmP9CMiJD0gHAocCyFq2y5WX6+eURSgMiYoekS4BHSa4SuS0iFkm6MF0/heSKn5OBpcAWkt9yylaBx/wdYH/gp+lv7DuijGdqLfCYdxuFHG9ELJb0CPAisBO4NSLyXnpaDgr8N/534A5JL5GcCroyIsp6SntJ9wAnAF0l1QBXAe2heT6/PPWKmZllwqe8zMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhSzjEjqIumrOe97Sbqvmfb1WUnfaaTNf0o6qTn2b5aPLxs2y4ikfsBDETGkBfb1NPCZhu6TkHQg8POIGNPc9ZiBRyhmWboOODh9hsgPJfXb9RwKSedJuj993sbrki6RdKmk+ZLmStovbXewpEfSCRn/IGlg3Z1IOgT4W0S8Lalz2l/7dN3ekpZLah8RbwD7S+rRgn8HVsEcKGbZmQy8FhFDI+J/5Vk/hGRa+JHAfwBbIuIo4Blg12y+U4GvRcRw4HLgp3n6+RjwAkBEbASeBD6VrjsT+E1EbE/fv5C2N2t2nnrFrOU8kQbARkl/BR5Ml78EHCGpE8nDy/47Z0bjD+fppyewNuf9rcAVwP0kU2d8OWfdGqBXVgdg1hAHilnL+VvO650573eS/L/4IWB9RAxtpJ/3gH12vYmIOenptU8AVXXm3OqYtjdrdj7lZZadjUDnYjdOn7vxuqQz4P3nfR+Zp+li4CN1lt0F3APcXmf5IUDZTupo5cWBYpaR9LkwcyS9LOmHRXZzNnCBpIXAIvI/ivgp4Ch98ElfvwL2JQkVANIv6j8CzCuyFrMm8WXDZmVI0o3AgxExM30/ATg1Ij6f02Y8MCwi/q2VyrQK4+9QzMrT90geeoWknwDjSJ5rkasd8KMWrssqmEcoZmaWCX+HYmZmmXCgmJlZJhwoZmaWCQeKmZllwoFiZmaZ+H+yd8uFZ3nA5gAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAElCAYAAADgCEWlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABqDUlEQVR4nO3ddZgV1RvA8e97Y7uAXdhi6e4QxEAsbLE7sLD4iWKhYhe22CKiIioqBiGolFIi3SENuwubbNeN8/tjLrC7LGxwY+N8nuc+996ZuTPvEOedOefMOaKUQtM0TdMOMfk6AE3TNK120YlB0zRNK0MnBk3TNK0MnRg0TdO0MnRi0DRN08rQiUHTNE0rQycGrd4SkedEZJLrc4KI5ImIuQb7+UREnnZ/hJpWO+nEoNVaIrJbRM4pt2yoiCyq7r6UUnuVUiFKKUcNfnuPUurFqmwrIl+KyEvVPYa71PTPR9NK04lB0+oIEbHUh2NotZ9ODFqdJiKxIvKTiKSJyC4ReeAY27UUEXWo4HP9bpqIZIrIdhG56zjHOHwXICKDRCRRRB4WkVQR2S8it7nWDQNuBB5zVVtNryxGEQkUka9E5KCIbBaRx0QksdT63SLyuIisA/JFxCIio0Rkh4jkisgmEbnctW0n4BNggOv4Wa7l4SIy0XX8PSIyWkRMrnVDRWSxiLwjIpnAczX9u9DqD311oNVZrsJtOjAVuB6IB+aIyFal1B+V/Pw7YCMQC3QEZovITqXU3CocOhoIB+KAc4EpIvKrUmqciJwCJCqlRlcxxmeBlkBrIBiYWcHxrgcuAtKVUnYR2QGcDhwArgYmiUhbpdRmEbkHuFMpdVqp37/virc10AT4E9gPfO5a3x+YDDQFrFU4f62e03cMWm33q4hkHXoBH5VadxIQpZR6QSlVopTaCXwGXHe8HYpIc+A04HGlVJFSag0wHri5ijHZgBeUUjal1EwgD+hwjG0ri/Ea4BWl1EGlVCLwXgX7eE8ptU8pVQiglPpRKZWslHIqpb4HtgH9jnGuZuBa4AmlVK5SajfwVrlzTVZKva+Ush86htaw6TsGrba7TCk159AXERkK3On62gKIPVRl4mIGFlayz1ggUymVW2rZHqBvFWPKUErZS30vAEKOsW1lMcYC+0qtK/25wmUicgswEuNOA9exI49x/EjAD+P8DtmDcbdzvGNqDZhODFpdtg/YpZRqV83fJQONRSS0VHJIAJLcEFP54Yori3E/RvXSJtf35sfbp4i0wLjjOBv4RynlEJE1gBzj+OkYdzgtSh2j/LnqIZa1MnRVklaXLQNyXI2zgSJiFpGuInLS8X6klNoHLAFeFZEAEekO3AF844aYUjDq8qsa4w/AEyLSSETigOGV7D8YoyBPA3A1fHctd/x4EfEDcHXP/QF4WURCXYllJDDpxE5Tq890YtDqLFehdwnQE9iFcXU8HqOhtTLXY1TFJAO/AM8qpWa7IazPgc6uNpFfqxDjC0Cia90cYApQfKydK6U2YbQR/IORBLoBi0ttMg+jUf2AiKS7lv0PyAd2AouAb4EJJ3qiWv0leqIeTas9RORe4Dql1Bm+jkVruPQdg6b5kIjEiMipImISkQ7Awxh3MJrmM7rxWdN8yw/4FGgFZGE8T/DR8X6gaZ6mq5I0TdO0MnRVkqZpmlaGTgxag1PRqK31RfkxoTStJnRi0OolV+GY7xpMLklE3pYazMXghhjaevOYmuYOOjFo9VkPpVQIxlPCNwDHHEFV07QjdGLQ6j2l1BaMsYm6ll8nIv1E5B/XA2n7ReSDQ08Nu9YrEblHRLa5hsb+UESk1PrbXcNlHxSRP1xPFiMiC1ybrHXdtVwrIpEiMsN1rEwRWXho+OsK4jpFRJaLSLbr/ZRS6/4SkRddw2XnisifInLUWEkicrWIrCy37GER+bV6f4JaQ6MTg1bviUhnjGGqV1ew2gE8hDHY3ACMu4v7ym1zMcYoqT0wRkM9z7Xfy4AngSuAKIzk8x2AUmqg67c9XDPHfY/xjEKia9tmrt8e1S1QRBoDv2GMtNoEeBv4TUSalNrsBuA2jKGy/YBHKji3aUAr1zwNh9wEfF3Btpp2WL1IDCIyQYxJUza4aX+vicgG1+tad+xT84lVInIQYz6E8cAX5TdQSq1USi11DTm9G+OZgvJPHY9RSmUppfYC8zGGtwC4G3hVKbXZNdrqK0DPQ3cNFbABMUAL15DdC1XF/cUvArYppb52xfUdsAVjaI1DvlBK/ecaJvuHUjGVPrdi4HuMZICIdMEYBmTGMeLTNKCeJAbgS+B8d+xIRC4CemP8R+sPPCoiYe7Yt+Z1vZVSjZRSbZRSo5VSzvIbiEh7V/XOARHJwSjcy1fLHCj1ufQQ2y2AsaXmisjEGOU0joq9AWwH/hSRnSIy6hjbxVJ2mGw4eqjsY8VU3lfADa7qr5uBH1wJQ9OOqV4kBqXUAoz/lIeJSBsR+V1EVrrqcjtWcXedgb9dV2r5wFrclHS0WuljjKvxdkqpMIzqHTn+Tw7bB9ytlIoo9QpUSi2paGPXRDkPK6VaY1z9jxSRsyvYNBkj6ZRWo2HBlVJLgRKMqrQb0NVIWhXUi8RwDOOA/yml+mDUv1Z1mIG1wAUiEuRq0DuTisfI1+qHUCAHyHNdPNxbjd9+gjFkdhc4PLfy1aXWlxmCW0QuFpG2rqv3HIz2DUcF+50JtBeRG8SY4/lajAuWmlYBTQQ+AOxKqUU13IfWgNTLh2BEJAQ4BfixVAcSf9e6KzCGOi4vSSl1nlLqT9dY+Uswxrz/B7BXsL1WPzyCcRHxGEbj9PfAWVX5oVLqF9e/tcmudoVsYDbwo2uT54CvRCQQGIZRFfQBRuPzQeAjpdRfFew3Q0QuBsZi3NFsBy5WSqWX37aKvgZedL00rVL1ZqwkEWkJzFBKdXW1CWxVSsW4Yb/fApNcc/tqWp3jSkypGG0u23wdj1b71cuqJKVUDrDr0G29GHpU5bdizLDVxPW5O9Ad+NNjwWqa590LLNdJQauqelGVJCLfAYOASBFJBJ4FbgQ+FpHRgBVjOOO1VdidFVjoqoLKAW4qN/G7ptUZIrIbozH9Mt9GotUl9aYqSdM0TXOPelmVpGmaptVcna9KioyMVC1btvR1GJqmaXXKypUr05VSURWtq/OJoWXLlqxYscLXYWiaptUpIlL+6frDdFWSpmmaVoZODJqmaVoZOjFomqZpZXi9jcE1veIKjCEoLi63TjCGAbgQY8TIoUqpVdU9hs1mIzExkaKiIneE7DEBAQHEx8djtVp9HYqmadphvmh8HgFsBioayvoCoJ3r1R9jnJj+1T1AYmIioaGhtGzZklJjJdUqSikyMjJITEykVatWvg5H0zTtMK9WJYlIPMYkJOOPsckQYKIyLAUiRKTa4x0VFRXRpEmTWpsUAESEJk2a1Pq7Gk3TGh5vtzG8izGK5VETprjEYYxxf0giFUx6IiLDRGSFiKxIS0urcEe1OSkcUhdi1DSt4fFaYnANI5yqlFp5vM0qWHbUmB1KqXFKqb5Kqb5RURU+n6FpmlavjZ2zjSXbazoS+/F5847hVOBS16Bek4GzRGRSuW0SKTspTjzGbFY+ccopp1S4fOjQoUyZMsXL0Wiaphky80t4d+5/LN990CP791piUEo9oZSKV0q1BK4D5imlbiq32TTgFtcw2ScD2Uqp/d6KsbwlSyqcoVHTNM2nFm5LQykY2L789OTu4fMhMUTkHgCl1CcYUxpeiDFjVQFwmw9DIyQkhLy8PJRS/O9//2PevHm0atUKPSKtpmm+tOC/dCKCrHSPj/DI/n2SGFzTGf7l+vxJqeUKuN8XMR3PL7/8wtatW1m/fj0pKSl07tyZ22+/3ddhaZrWACmlWLgtjdPaRmI2eaYDi37yuQoWLFjA9ddfj9lsJjY2lrPOqtKUwJqmaW635UAuqbnFDGzvuY43OjFUke5aqmlabbDgP6OL/sB2OjH41MCBA5k8eTIOh4P9+/czf/58X4ekaVoDtWBbGh2ahRIdHuCxY/i88bkuuPzyy5k3bx7dunWjffv2nHHGGb4OSdO0Bii3yMayXZncdqpnh9HRieE48vLyAKMa6YMPPvBxNJqmNXQLt6VjcyjO7tgU8tMhqAl4oJpbVyVpmqbVEXM3pxIeaKVPRB6MGwTzXvTIcXRi0DRNqwMcTsX8ralc0saM5ZsroCgHOg/xyLF0VZKmaVodsGZfFrb8LB5NexMKkuHmXyCmh0eOpRODpmlaHfDXxn185vc2Ybk74IbvIeFkjx1LJwZN07Tazumkz+onOdm0CS77DNqe7dHD6TYGTdO0Wi532uMMsi1kaZsR0P0ajx9PJwYPuf3222natCldu3b1dSiaptVlS94ndM04vrCfR+wFj3vlkDoxeMjQoUP5/ffffR2Gpml12fop8OdoFvufxs9N7ychMtgrh9WJwUMGDhxI48aNfR2Gpml11b7l8Ot9FMedzO3Zd3JB96NmOfaYet/4/Pz0jWxKznHrPjvHhvHsJV3cuk9N006Mw6lYtiuTfQcLaN4oiJNaNsJirqPXvln7YPINEBbDD61foXjHAS7sGuO1w9f7xKBpWv23KTmHEZNXsy017/CyVpHBvDCkC6d7cBRSjyjOg++uB3sR3DqdX6Zk0DkmjJZeqkYCLyYGEQkAFgD+ruNOUUo9W26bQcBUYJdr0c9KqRdO5Lj6yl7T6rcNSdlc++k/hAZYGXtdT3o2j2BdYjbvzPmPWyYs4+mLOnP7aZ4ddM5tnE74+S5I3Qg3/kiyXwtW7d3Bo+d18GoY3rxjKAbOUkrliYgVWCQis5RSS8ttt1ApdbEX49I0rY7KzC/h9i+XExHkx0/3nnJ4KOoWTYI5u1NTRn6/lhdmbEIEj49I6hZzn4etM+GC16HtOUz9awcAF3f3XjUSeLHxWRkO3edZXa96O3ny9ddfz4ABA9i6dSvx8fF8/vnnvg5J0+qdF6ZvJDO/hHG39DlqfoIgPwvv39CL87o04/npm5i9KcVHUVbRmm9h8bvQ93boNwylFD+tSqRvi0a0aOK9aiTwcq8kETGLyBogFZitlPq3gs0GiMhaEZklIhXWA4nIMBFZISIr0tLSPBlyjX333Xfs378fm81GYmIid9xxh69D0rR6ZfH2dH5dk8x9Z7alS2x4hdtYzSbGXteLbnHhPPT9GraXaoOoVZJXw/QHodVA425BhPVJ2WxPzePKPvFeD8eriUEp5VBK9QTigX4iUv7pr1VAC6VUD+B94Ndj7GecUqqvUqpvVFQda1jSNO2EKaV4/Y+txEUEcv+ZbY67bYDVzKc398HfYuL+b1ZRZHN4KcoqKsiEH26B4Ci46kswWwH4aWUifhYTF3bzbjUS+Og5BqVUFvAXcH655TmHqpuUUjMBq4hEej1ATdNqtbmbU1m7L4sHzm6Lv8Vc6faxEYG8eU0Ptqbk8sYfW70QYRU5nfDL3ZCzH675CoKbAFBidzJtbTKDOzcjPNDq9bC8lhhEJEpEIlyfA4FzgC3ltokWMaYjEpF+rvgyvBWjpml1w6cLdtC8cSBX9K56NcuZHZpyy4AWfL5oF4u2pXswumpY9BZs+xPOfxXi+x5ePGvDfg4W2Li6b3OfhOXNO4YYYL6IrAOWY7QxzBCRe0TkHtc2VwEbRGQt8B5wnVKq3jZQa5pWfRuSslm++yC3DmiJtZoPsD1xQSfaRAXz8I9ryCoo8VCEVbRjHsx7GbpdAyfdWWbVxH/20LJJEKe39U2FiTd7Ja1TSvVSSnVXSnU99HyCUuoTpdQnrs8fKKW6KKV6KKVOVkot8VZ8mqbVDV8u2U2Qn7lGV9OBfmbGXteLjLwSRv+6wQPRVVF2Ivx0J0R1hEveLTNv88bkbFbuOchNJ7fAZHL/fM5VUUefF9c0rSHKyCtm2tpkruwdX+O6965x4Tx4TjtmrNvP1DVJbo6wCuwl8MOtYC+Ga78Gv7JdUb/+Zw8BVhNX9/FNNRLoxOAx+/bt48wzz6RTp0506dKFsWPH+jokTavzpq1NpsTu5KaTW5zQfu45ow29EiJ4+tcN7M8udFN0VfTnU5C0AoZ8CJHtyqzKyCvm1zVJDOkRR3iQ9xudD9GJwUMsFgtvvfUWmzdvZunSpXz44Yds2rTJ12FpWp02dU0ynWLC6BAdekL7sZhNvH1NT2wOxWNT1uG1psz1U2DZOBgwHLpcdtTqL5fsptju5K6Bvn1KWycGD4mJiaF3794AhIaG0qlTJ5KSfHDbqmn1xO70fNbsy+KynrFu2V+ryGCevKgTC7el8/XSPW7Z53GlboZp/4OEAXDOc0etzi2y8eWS3Qzu3Iy2TU8s8Z2o+j+66qxRcGC9e/cZ3Q0uGFPlzXfv3s3q1avp37+/e+PQtAZk6ppkROCSHu5JDAA39U9gzqYUXpm5mVPbRtImKsRt+y6jOBe+vxn8QuCqLw4/xFbapKV7yS2yc9+gtp6JoRr0HYOH5eXlceWVV/Luu+8SFhbm63A0rU5SSjF1TRL9WjYmNiLQbfsVEV6/qjsBVjMjf1iL3eF0274PUwqmDofMHXDVBAg7+knm7EIbn/y9g4Hto+jRPML9MVRT/b9jqMaVvbvZbDauvPJKbrzxRq644gqfxaFpdd2GpBx2pudz18DWbt93s7AAXrqsK8O/Xc2H83cw4px2lf+oOv79BDb9Cuc8D61Or3CTj//aQU6RjcfP9+7w2sei7xg8RCnFHXfcQadOnRg5cqSvw9G0Ou3XNUn4mU0em8Xs4u6xDOkZy3vztvHPDjcOtrB3Kfw5GjpeDKeOqHCT5KxCvli8i8t6xh1zMEBv04nBQxYvXszXX3/NvHnz6NmzJz179mTmzJm+DkvT6hyHUzF9bTKDOkR5tAvnS5d1pWWTIIZ/u8o9XVjzUuHHoRCRAJd9VOYhttKen74RgJHntj/xY7pJ/a9K8pHTTjvNe13gNK0e+2dHBqm5xQzpGefR44QGWPn05r5c9uFi7pm0ih/uPrlKA/RVyGGHKbdD4UG48UcIqPhO4M+NB/hjYwqPnd+B5o2DTiB699J3DJqm1Wq/rkkixN/C2Z2aevxYbZuG8ObVPVi7L4tRP63H6azhxd38l2D3Qrj4HaMXYwWyCkp4dtpGOjQL5a7T3d92ciJ0YtA0rdYqsjn4fcMBzu8aTYC1hlfv1XR+12geGdyeX1Yn8dofWyr/QXlbZsKid6DPUOh5Q4WbKGU8WJeeV8zrV3Wv9mCAnqarkjRNq7XmbUklr9jOZR6uRirv/jPbkpJTzKd/7yQi0I97Bx1/MqDDMnfCL/dATE84/7VjbvbF4t38uSmF0Rd1qhXdU8vTiUHTtFrr19VJRIX6M6BNE68eV0R47tIuZBXaeO33LRTZHDx4TjvkGA3IANgK4ftbjEbmayaCNaDCzeZsSuGl3zZxbudm3HGab4e+OBadGDRNq5WyC2z8tTWNm05ugdkHw0+bTcK71/YkwGJi7NxtHCwo4emLO1dc7aMU/PYwpKyHG36ERhUP8rdoWzrDv1tF17hwxl7X8/iJxod0YtA0rVaatWE/JQ4nl/Vy3xAY1WU2Ca9d2Z2IICufLdzFlgO5fHB9L5qGlbsbWDUR1nwDAx+D9oMr3Ndv6/bz0PdraB0VzIShJxHkV3uLX29O7RkgIstEZK2IbBSR5yvYRkTkPRHZLiLrRKS3t+Jzt6KiIvr160ePHj3o0qULzz77rK9D0rQ65dc1SbSODKZbnG8f+jKZhKcu6sw71/ZgXWIW57z9N98v33ukx1LyGpj5KLQ+EwaNOur3BSV2np++kfu/XUXXuDAmDzuZyBB/755ENXkzZRUDZyml8kTECiwSkVlKqaWltrkAaOd69Qc+dr3XOf7+/sybN4+QkBBsNhunnXYaF1xwASeffLKvQ9O0Wm9/diH/7spkxNmV1Ot70eW94ukeH8ETP6/n8Z/WM2HRboYPaMxF/9yMKTgKrvwcTEd6TuUX2/l5dRIfzd/O/uwihp7Skicv7ISfpXb1QKqI1xKDa+7mPNdXq+tVvpPwEGCia9ulIhIhIjFKqf3eitNdRISQEGOkRpvNhs1mqzX/wDWttpu2Jhml8PhDbdXVJiqEyXedzLS1yXw07z+CfhuO3ZTES03fxjT3AKEB6eQW2dmemsfy3ZkU2530iA/n/et70bdlY1+HX2VereQSETOwEmgLfKiU+rfcJnHAvlLfE13LyiQGERkGDANISEg47jFfW/YaWzJr0Bf5ODo27sjj/R6vdDuHw0GfPn3Yvn07999/vx52W9Oq6JfVSfRsHkGryODKN/Yyk0m4rFccl+Z8i2n+an5r/jCrC9qya2Ui+SV2QvwsxDcO4ob+CVzcPYbeCY3q3EWhVxODUsoB9BSRCOAXEemqlCo9I3dFf3pHPXqolBoHjAPo27dvrR13wmw2s2bNGrKysrj88svZsGEDXbt29XVYmlarbUrOYcuBXF4Y0sUnx7c5bWQVZXGw+CD5tnyK7EUUO4opshdR5CiixFECaf8hKz5BOp6B9O7MbZKK1WTF3+xPoCXQ9YIAy0EO5BcT7h9OoCWwziQInzSLK6WyROQv4HygdGJIBErPgB0PJJ/IsapyZe9pERERDBo0iN9//10nBk2rxK9rkrCYhIu7u783ksPpIDEvkX25+0jOSz7yyk8msyiTrKIscm25VdtZZCMo3gX/HNWPpkJWk5Vw/3Ai/CMI8wsjwj+CRgGNaBLYhKaBTYkMiiQqMIqowCgiAyOxVjCZj7d4LTGISBRgcyWFQOAcoPyjgdOA4SIyGaPRObsuti8ApKWlYbVaiYiIoLCwkDlz5vD4475PUppWmzmcxoQ8gzo0pXGw3wnt62DRQdanr2dr5la2Z21nR9YOdmXvosRZcngbi8lCdFA0sSGxdIvsRqOARkaB7d+IiIAIQqwhBFgCCDAH4G/2J0BM+P0wFDJ3oW6agnI9r6CUwua0UWgvpMhRRJG9yPhsLyLflk92STZZxVnkFOeQVZxFdnE2e3P3sjZtLQeLD+JUR08Q1Mi/EZFBkTQNbEqz4GZEB0XTLLgZzYKaER0cTbOgZoT4eWbGOW/eMcQAX7naGUzAD0qpGSJyD4BS6hNgJnAhsB0oAG7zYnxutX//fm699VYcDgdOp5NrrrmGiy++2NdhaVqttmRHOik5xTx7SfUanZ3KybaD21iRsoJ1aetYn76efblHmitjgmNoE9GGAbEDaBPRhhZhLYgJjiEqMAqzqRpjME1/EJLXwrWTIN49bYZ2p53MokzSCtNIL0gnrTCNtII0470wjdSCVLYe3Ep6YfpRv72ty22M7Ov++V682StpHdCrguWflPqsgPu9FZMnde/endWrV/s6DE2rU35ZlURogIWzOlY+kmpyXjJL9y9lafJS/j3wL5lFmQBEBUbRPao7V7W/im6R3ejUuJN7rqxXT4KVX8BpD0GnS058fy4Wk4WmQU1pGtQUjjPyh81hI7UwlZT8FA7kHyClIIXOTTq7LY4yMXlkr5qmadVUUGLn940HuLRHbIUjqSql2JK5hbl75zJ371y2Z20HIDIwklNiT+HkmJPpF92P6OBo9zfyJq+BGSOh1Rlw5mj37ruKrGYrcSFxxIV4vguvTgyaptUKM9btp6DEwZV94g8vO5QMZuycwdy9c0nKS8IkJno17cWjfR/llNhTaBPRxrO9fQoy4YebITgSrpoA5vpfbNb/M9Q0rU749t+9tG0aQt8WjUjJT+G3Xb8xfcd0tmdtx2KycGrsqdzd/W7OaH4GjQO89LCY0wk/3wU5++H2343k0ADoxKBpms9tSs5hzb5MbjqzkPvn3s+ipEUoFD2ievD0yU9zXsvzCPf3wZhJf4+B7XPgorchvq/3j+8jOjFomuZTB4sO8sLCjwlp+wdTD2QSGRjJXd3v4tI2l9IirOLhq71i0zT4+zXoeSP0vd13cfiATgyapvnEvtx9fLXxK37d/ivFjmIa+XVk9MAnOSvhLKwm3z3cBcCBDfDL3RDX17hbqCNPLLuLTgwe5nA46Nu3L3FxccyYMcPX4Wiaz23N3MrnGz7nj91/YBYzXcPOYtHKznx+22X0qw0DzeWnw3fXQ0A4XPfNMWdiq88qTQwicvxR6o7IUkrlnGA89c7YsWPp1KkTOTn6j0Zr2DZlbOKD1R+wMGkhwdZgbu18Kzd2upFbP9tKu0ZwUstGvg4RHDb44VbIS4HbZ0FotK8j8omq3DF8hTGQ3fHupRTwJTDRDTHVG4mJifz222889dRTvP32274OR9N8YkfWDj5c8yGz98wm3D+cB3o9wLUdryXML4wF/6Wx5UAub1zVvXYMMDfrcdizCK74DOL6+Doan6k0MSilziy/TESilVIHPBOSex145RWKN7t32G3/Th2JfvLJSrd78MEHef3118nNreKgXJpWjyTlJfHRmo+YsXMGgZZA7u1xLzd3vplQv9DD24xftIuoUH8u7em76TsPW/45rPgcTh0B3a/xdTQ+VdM2hluA190ZSH0zY8YMmjZtSp8+ffjrr798HY6meU2+LZ/x68czceNERIRbOt/C7V1vp1FA2aqirQdyWfBfGo8Mbo+/pRrjFXnC9rnG9JztBsPZehremiaGISJSAMxWSm11Z0DuVpUre09YvHgx06ZNY+bMmRQVFZGTk8NNN93EpEmTfBKPpnmaUzmZun0q761+j/TCdC5pfQkP9H6A6OCK6+k//XsHAVYTN/b3YZdUgJSNRrtC005HTc/ZUNU0MVyBMSDe5SLSVil1pxtjqhdeffVVXn31VQD++usv3nzzTZ0UtHprTeoaXl32KpsyNtE9qjtjzxxL96jux9x+R1oev65J4o7TWtHoBIfXPiE5++Gbq8E/BG74AQLCfBdLLVKjxKCUSgF+d700TWugckpyeHflu/z43480DWrKq6e/ykWtLqq0Ifm9udvwt5i5+4w2Xoq0AsV58O01UJQNt82C8No1v7Qv1SgxiMiHQLBSaqiIDFZK/enmuOqVQYMGMWjQIF+HoWluo5Tijz1/8Nqy18gsyuTmzjczvOdwgqxBlf52e2ou09YmM2xgayJD/L0QbQUcdphyu1GNdMP3EHPsu5uGqKZVSSVAiuvzWYBODJrWQCTnJfPS0pdYmLSQzk068+HZH1ZrXoBXZ24h2M/CsNNbezDK41AKfhsJ2/6Ai9+Bduf6Jo5arKaJoQAIFxErUKUH4ESkOcZzDtGAExinlBpbbptBwFRgl2vRz0qpF2oYo6ZpbqSU4qdtP/HG8jdQKB476TGu73g9FlPVi5G//0tj7pZUnrigI018dbcw5zlY9RWc/kidHgNJKeWxZz9qmhgygULgQ2BxFX9jBx5WSq0SkVBgpYjMVkptKrfdQqWUngNT02qR1IJUnl3yLIuSFtEvuh8vnvoisSHVe/bA5nDy4oxNtGgSxNBTW3om0MosehcWv2skhLN8M+GOO9hSUznw9DOED7mUsAsvdPv+TdXZWEQiROQL4ErXoolAlcaiVUrtV0qtcn3OBTYDurVH02oxpRS/7fyNy6dezooDK3ii3xN8NvizaicFgM8X7WJ7ah5PXdjJN88trPwS5jwLXa+EC9+sswPj5fz+O7suuZT8pUtx5Od75BjVumNQSmWJyBigJZAOdAd+ru5BRaQlRnfXfytYPUBE1gLJwCNKqY0V/H4YMAwgIaGqQzlpWt23Oz2fv/9LIzm7kACLmQ7RoQxo3cQjXT6zi7N54Z8X+HPPn/SI6sHLp71c42Gwt6fm8vbs/zi/SzTndm7m5kirYOMvMOMhaHsOXPZJnXxWwZGdzYEXXyJnxgwCunUj9rUx+Lf2TDtNTaqS7gB2KaX+AFZW98ciEgL8BDxYwaB7q4AWSqk8EbkQ+BVoV34fSqlxwDiAvn37qurGoGl1TXJWIc9P38gfG40+H1azYHOow58Hd47m3kFt6Brnnsls1qSu4bEFj5FWkMaI3iO4rcttmGtYmDqcikenrCPIz8yLl3X1/phIm6bClDsgvh9c8zVYfPjcRA3lLVzE/qeewp6ZSeQD/yNy2DDE4rnBsWuy54PAPSLSAVgLrFFKra7KD12N1T8B3yiljrrTKJ0olFIzReQjEYlUSqXXIE6fa9myJaGhoZjNZiwWCytWrPB1SFodtHLPQe7+egUFJQ4ePKcdV/aOJ75RIDaHYn1SNr+t28+Ulfv4bf1+LuwWzZMXdiK+UeXdRivicDqYsGECH675kJjgGCZeMJFuUd1OKP53Zv/H6r1ZjL2uJ1GhXm5w3jQVfrzNmH3tpingV7M/F19xFhWR+sabHPzmG/zatqHlRx8R2LWLx49b7cSglHpVROYC/wE9gYFApYlBjMuEz4HNSqkKhxoVkWggRSmlRKQfRhtIRnVjrE3mz59PZGTDmCdWc7/1idnc8vm/RIX6M3nYANo2DTm8zs8i9GnRiD4tGvHgue34fOEuxi3YyfwtaTxwdjvuOK0VfpaqNyOmFaTxxKIn+Hf/v5zf8nyeGfBMmQHvamLelhQ+mL+da/s2Z0hPLzcpHk4KJxlJwf/EzsXbirZsIemRRyjZvoPGt95C1MiRmPy9k1irnRhE5AXADKzBuFv4q4o/PRW4GVgvImtcy57E1d1VKfUJcBVwr4jYMXo9XaeU0lVFWoOUmlvEbV8uJyLIj8nDBhAdfuwJY8ICrDx0bnuuOak5L0zfyGu/b+GnVYm8cnk3+rWqfPKbpfuX8viCxymwFfD8Kc9zedvLT7jKZ/P+HEZ8t4bOMWE8P8TzV7llrJ8CPw+rk0lBOZ1kfvkVae+8gykinObjxxNy2qlejaEmdwzPiEgzjMbjK0WkjVLqrir8bhHHn9MBpdQHwAfVjel4Fv7wH+n78ty5SyKbh3D6Ne0r3U5EGDx4MCLC3XffzbBhw9wah1Z/OZ2Kh39YS26RjW+Gn3bcpFBaXEQgn97cl7mbU3hm6kau+fQfru3bnCcu7EhE0NF160opJmyYwHur36NVWCsmnDeBNhEnPkxF4sEChn6xjGB/C5/d2pcAqxcbe5d9ZoyU2uIU46nmOpQUbCkpJI8aRcE/Swk552xiXnwRSyPvT2BU09aLu4FPlVJ6rKTjWLx4MbGxsaSmpnLuuefSsWNHBg4c6OuwtDpg0r97WLgtnZcu60qH6OoXbGd3asaANk0YO2cb4xftYs7mFEZf3InLesYdvhPIK8nj6cVPM2fvHM5reR4vnPJClYa0qMzu9HxuHP8vBSUOfrxnAHERgSe8zypRCv4aA3+PgQ4XwVUT6tS0nDl//MmBZ57BWVJC9IsvEHHVVT6bvKimiWECRpVPMEZD8hr3heReVbmy95TYWKOvd9OmTbn88stZtmyZTgxapdJyi3njj62c1jaSG/vXvDt2kJ+FJy7sxJCecTz5y3oe+n4tU1Ym8uKQrmBNZcT8EezL3ccjfR/hls63uKUQWrsvi7smrsDmcPLtnSfTMdpLo5U67DDrMWOinZ43wSVjwVw3prR35udz4JVXyP7pZwK6diX2jdfxb9XKpzFV6wG3Uh7ASCoW4D33hVN/5OfnH565LT8/nz///JOuXbv6OCqtLnjt9y0U2Rw8P6SLWwrrzrFh/HzvKbx0WVfWJWZz/mcfcsXUa8kqzuGzwZ9xa5dbT/g4SikmLd3D1Z/8g9Vs4vu7B9At3j1dZytVlG2Mknpo9rUhH9SZpFC4bh07r7iC7J9/ocndd9Pyu299nhSg5ncMOzCeL5iqlHrIjfHUGykpKVx++eUA2O12brjhBs4//3wfR6XVdhuTs5myMpF7zmhDm6iQyn9QRSaTcGP/BPbLVCZumYStMIGUXTfxW1AQTQbk0foEjrU9NZdnpm5kyY4MTm8XyXvX9fLeHAuZO+Hb6yBzh3GX0Geod457gpTTSeaECaS+OxZL0yhaTPyKoJNO8nVYh9U0MWwE9gF3iMgbSqnac0a1ROvWrVm7dq2vw9DqmA/mbSc0wMK9g9w7T0GRvYjRi0fzx+4/uLTNpdzR4VE+/msP3/y7hy+X7Ob0dpFc2iOWczo1q1Kh7nQqlu/O5Jt/9zJ9XTIhfhZevrwr15+UgMnkpXrxHfOMB9dQcPOv0Op07xz3BNkzMkh+fBT5ixYRet55xLz4Auaw2jVBUE0TQ3sgDePp44PuC0fTGq6tB3KZteEAD5zdjvBAq9v2m1qQygPzHmBTxiZG9hnJ0C5DERHeuDqCx87vyPfL9/Ldsn08OmUdZpPQLS6cns0jaNcshOiwAMIDrSggq8BG0sEC1iVms3RnBsnZRQT7mbl7YBvuOr2V90ZLdTqMRuYFb0BUB7juW2jiwwl/qiH/n39IeuwxnDm5RD/3HBHXXuOzBubjqWli6IgxB8OHwB6MNgdN007Ae/O2EeJv4XY3jjy6MWMjD8x9gFxbLmPPHMuZCWeWWR8V6s/ws9px/5ltWZ+UzZ8bU1i2O5Pvl++j0OaocJ+RIX6c1LIxj57fjPO6RBPk58X6/NwU+OkO2L0Qet4IF74BfsHeO34NKbudtA8+IOPTcfi1bk3C+M8J6OC7jjGVqenfaATwOPAYxthJmqadgG0pucxcv5/7BrWp8HmDmpi9ZzZPLnySRgGN+PqCr+nQuMMxtxURusdH0D0+AjDGN0rJKeJAThF5RXZEIDTASmxEAFEh/t6/ylUKNvxkPJ9gK4QhH0GvG70bQw3ZkpNJeuRRCletIvzKK4h+6ilMQbV7aI6aJoYXgI5Kqa0i4nRnQJrWEH0wfzuBVjN3nHbio2UqpZi4aSJvrniTHlE9ePfMd4kMrN6wLGaTEBsRSKy3nkE4nrw0Y8a1zdMgri9c9jFE1d6r7dJy58wh+anRYLMR+8YbhF9SN6aaqVJiEBEzkAg8rZQar5RKdH1HKTXKg/FpWr23Iy2P6WuTuWtgaxqfYG8eh9PBmyveZNLmSZzb4lxePf1V/M0+mintRDnsRhfU+S8bdwnnPA8DhteJrqjOkhJSX3+Dg5MmEdC5M3HvvI1fi5oNWe4LVfoTVko5RGQDUDdaeDStDvlw/nb8LCbuOsE5kIvsRTy56Elm75nNTZ1u4tGTHsUkNX1Uycd2LYBZoyB1I7Q6w2hLiDp2VVhtUrxrF0kPP0zxps3G4HcPP4zJr24N9V2d1BsEPCYi52JMogOglFJD3B9W/ZCVlcWdd97Jhg0bEBEmTJjAgAEDfB2W5m25KZD+H2Rsg4wdkJMM+WlQkIE9/yBP5eXzip+TgHcdYLIYwzhYAo0hooOaQHAUhDQ1XuEJ0LgVNG5trHPV9WcVZfHA/AdYk7qGR/s+yi1dbvHxSdfQnn/gr1eMxBCeYMyf0OmSOjPbWvZvv3Hg6WcQq5X4jz4i9KwzK/9RLVSdxHCoROvtegHokU+PY8SIEZx//vlMmTKFkpISCgoKfB2S5mlKwYH1RsGWuNx45SQdWW8JgLA4o7Bv3JrVJU62KRuX9WgJAQFG9Ym9EGxFYMuH/AxI2QA70qA4u+yx/EKhcSuSIltxT8lOku35vNHvSc7reK1XT/mEOZ2wfTYs/Qh2/gXBTeG8V6HvbWCtBW0cVeAsKSF1zGsc/PZbAnv3Ju7tt7BGR/s6rBqrTmLw/XPadUhOTg4LFizgyy+/BMDPzw+/OnY7qVWRw248bLVlOmybDbn7jeURCZBwstFg2rQjNGlnJAWTUb2zL7OA69/8i5sHtCDokioMS20rgqy9cHCX8cRv5i62pG/k3txVFCvFuJQ0+uy8GwJHQdMu0KwLRHeFZl2haafaV8jm7IcNU2D5eDi4G0Jj4NwX4aQ769SEOrakJBIffIii9etpfNttNB35EGJ133MovlDlxKCU2uPJQDxl/pfjSN2z0637bNqiNWcOPf4Q2jt37iQqKorbbruNtWvX0qdPH8aOHUtwcO3vc61V0YH1sPobo3DLTwP/MGhzJrQbDG3OhrCY4/78o792YDIJ95xRxaY7a4DRG8fVI2dlykqGzx1OSFAM4097jTZFBZCyyaiXT9kIqycZdx0AYoImbY0k0awLRHcz3sPivFtNk7UPtv1pzMG8exGgIGEAnP2sUWVkrlsFat7ff5P02OPgcBD33ljCBg/2dUhuUfub9+sou93OqlWreP/99+nfvz8jRoxgzJgxvPjii74OTTsRTgdsnQX/fmI8ZGX2g/bnQ4/roO25VZ5POCmrkCkr93F9vwSahVV/aOgFiQsY+ddIYkNiGXfuOKKDXdUWLU8rFavTuLtI2WhURx3YAEkrYWOpWXUDGx1JFs26QlRH404npOmJJwyHzWhb2b8WElcY1USZO4x1TdrBoFHQ5Yo60/W0NOVwkPb++2R88in+HTsSP/bdOtXrqDJeSwwi0hyYCEQDTmCcUmpsuW0EGAtcCBQAQ5VSq07kuJVd2XtKfHw88fHx9O/fH4CrrrqKMWPG+CQWzQ0cdlj/A/z9ulHYhjc3qj163QRBlc+QVt5H87cDVP1uoZTfdv7G6EWjad+4PR+f8zGNA45xfJPJGCqiSRvofOmR5UU5kLrJuOM5lDRWfX3k7gKMtpDw5kaSCI02EkhghPFuCQAxg8ls3ImU5ENxLpTkQe4ByN5n3Bkc3A2OYmN/fiHQ4lTodxe0HmQkoDrSoFyePT2dpEcepWDpUsKvupLo0aMxBdSdeR+qoiZTe16ilJpeg2PZgYeVUqtEJBRYKSKzlVKbSm1zAcaore2A/sDHrvc6Jzo6mubNm7N161Y6dOjA3Llz6dy5s6/D0qrL6YRNv8D8V41eRTE94OqvoOPFNe5Pn3iwgB9W7OPak5pX+wGyyVsm88q/r9CnWR/eP+t9QvxqMCpqQJjR9pFw8pFlh+4uMrYb7RhZe1zve42r/oLMsonjWIKaGAklqj20HwzRPSCmu1GNZfLiLG4eUrB8OUkjH8aRk0PMK68QccXlvg7JI2ryL/tloNqJQSm1H9jv+pwrIpuBOKB0YhgCTHTN87xURCJEJMb12zrn/fff58Ybb6SkpITWrVvzxRdf+DokrTqSVsHMR4zql6hOcO0kIyGc4JXuh/N3IAj3n9m2yr9RSjFu3Tg+WPMBg+IH8cYZbxBgceNVaum7i2OxF0PhQeNdOYxkopzGWEX+IcZdQT0o/CuilDKGyX77Hfzi42k+/jMCOtSN5ypqoiaJ4YTv/0SkJcac0f+WWxWHMZz3IYmuZWUSg4gMA4YBJCTUfIYrT+vZsycrVqzwdRhadeVnwNznYdVEo1vpZR9D92vdUujtyyzgxxX7uKF/AjHhVbtbcConb654k683fc3FrS/mhVNfwGryQSOtxd+oVmpgHDk5JD/xJHlz5xrDZL/8EuYQ982VURvVJDGc0LMLIhIC/AQ8qJTKKb+6KsdTSo3DGPKbvn376mcpNPfZNBVmPASFWTDgfjjjcaPqxU0+nL8dkwj3Dara3YJTOXluyXP8sv0Xbuh4A4/3e7zuPs1cBxVu3EjSiAexHThAsyefoNHNN9fKYbLdzau9kkTEipEUvlFK/VzBJolA81Lf4znylLWmeU5BpjFy54YpENMTbp0BzdzbJrQzLY8pKxO5sX8C0eGVVwM5nA6eXfIsU3dMZVj3YQzvObxBFEq1gVKKrB9+JOXllzE3bkyLrycS1KuXr8PyGm/2ShLgc2CzUurtY2w2DRguIpMxGp2z62r7glaH7JgHv9xjJIczR8NpD3qkP/2YWVvwt5gYfla7SrctnRTu7XEv9/W8z+3xaBVzFhRw4PnnyZ46jeBTTyX2zTewNGrk67C8qiaJIaWGxzoVuBlYLyJrXMueBBIAlFKfADMxuqpux+iuelsNj6VplSs/E9iNU4weNB6wdGcGf25K4ZHB7YkKPf5opw6ng2eWPMO0HdN0UvCy4p27SBrxAMXbdxA5fDiR996DmOtng/rxVDsxKKXOrcmBlFKLqKTh2tUb6f6a7F/TqsWLM4E5nYqXf9tMTHhApfMtlE4K9/W4j3t73uuRmLSj5cyaxf6nRiP+/jQf/xkhp57q65B8Rj/5rDU8SSvhuxugKNsrM4FNWZXI+qRs3r6mB4F+x776LJMUet7HvT10UvAGVVJCimvuhMCePYl79506PQCeO+juDR6ydetWevbsefgVFhbGu+++6+uwtA0/wRcXGkNX3DXX40khI6+YV2Zupk+LRlzWM+6Y2zmcDp5e/DTTdkzj/p7366TgJbbkZHbffDMHJ02i8a230uLriQ0+KUAN7xhEZOShBmQR6aCU2uresOq+Dh06sGbNGgAcDgdxcXFcfnn9fEqyTnA64e8x8PdrkHAKXPs1BFdvusuaeOm3zeQX23n1im6YTBXXpB5KCtN3Tuf+nvdzT497PB6XBnkLF5L8yKMou524sWMJO69+DIDnDtVKDCISAbwDdBSRImAdcAe6kfi45s6dS5s2bWhRjwbZqlPsJfDrvUZX1J43wcXvVHmwuxPx939p/LI6if+d1Zb2zUIr3MbhdDB68Whm7JzB8J7DubvH3R6Pq6FTDgdpH3xgDIDXvr0xAF7Llr4Oq1apVmJQSmUBt4nIRcABYDBQ0fMItUbW9B2UJFdhjJdq8IsNJuKSqg9+NnnyZK6//nq3xqBVUXEe/HCz0SX17GfhtIe8MnhbRl4xj/y4lrZNQ4459IXD6eCpxU/x287f+F+v/zGsu28GfGxI7BkZJD3yCAX/LCX8iiuIfubpejcAnjvUtI3hDIxuqycDNeql1FCUlJQwbdo0rr76al+H0vDkp8NXl8DOv2HIh3D6SK8kBaUUj01ZR3aBjfeu60WA9egGZ50UvK9g5Up2XX4FhatWE/PyS8S+8rJOCsdQ015JEcDjwGMYVUm1VnWu7D1h1qxZ9O7dm2bNmvk0jgYnay98fTlkJ7oGv7vQa4f+bOFO5m5J5ZmLO9M59ujhNOxOO08teoqZu2byQK8HuKv7XV6LrSFSSpH5xZekvvUW1vg4Wo77lICOHX0dVq1W08TwAtBRKbVVRJzuDKi++e6773Q1krcd3A1fXmzMO3DzL9DiFK8dev6WVF6dtYWLusVw26ktj1pfOimM6D2CO7vd6bXYGiJHTg77n3qK3NlzCD33XGJeeRlzaMXtPdoRNUoMSqlEjHGNUEqNcmtE9UhBQQGzZ8/m008/9XUoDUfmTvjyEmPSmFunQWxPrx16Y3I2D3y3ms4xYbx5dY+jxjWyO+08uehJZu2apZOCFxRt3kziiAexJSfTdNTjNL71Vj3WVBXVqI1BRD4UkS9dn3Ufr2MICgoiIyOD8PBwX4fSMGTsMO4UbPleTwrbU3O5+fNlhAZY+OyWvkc9yGZ32nlyoU4K3qCU4uCPP7L72utQxcW0mPgVTYYO1UmhGmpalVTCkTGTzgL+dE84mlZD6dvhq4vBUQK3Tjcmu/eSLQdyuOXzZZhNwjd3nXzUrGyHk8LuWTzY+0Hu6Farm+XqNGMAvBfInjqV4FNOMQbAa1z9qVcbupomhgIg3DWMdu2dKUdrGDJ3uZKCzSPDZR/PPzsyGDZxBUH+Zibd0Z9WkWXHWyqdFB7q8xC3d73da7E1NHoAPPepaWLIBAqBD4HF7gtH06opJxkmDgF7EQyd6bWkoJRi/MJdvPb7FlpFBvPl7f2Iq+BO4YmFT/D77t8Z2Wckt3XVz4F6Ss7Mmewf/bQxAN5nnxFyWsMdAM8davrkcwdgEjCRWt5dVavH8jNg4mVQkGG0KXgpKSRlFTL6l/XM35rGeV2a8fpVPQgPLDt/g91pZ9TCUfyx+w+dFDzIWVJC6muvc/Cbbwjs1Yu4d97WYx25QbWffBaRMUBLIB3oTi1/8lmrp4qyYdIVkLXHmEchro/HD1lY4uCrf3bz3txtKAXPXtKZoae0rLD30aGk8HCfhxnadajHY2uIbElJJD74EEXr19N46FCaPjwSsfpgLux6qCZVSXcAu5RSfwAr3RyPplWupAC+vQ5SNsB130Kr0z16uIP5Jfy4ch/jFuwiPa+Yszs25fkhXYhvFHTUtjanjVELRvHnnj95pO8j3NrlVo/G1lDl/vUXyY+PAoeDuPfGEjZYd450p5okhoPAPSLSAVgLrFFKra7sRyIyAbgYSFVKda1g/SBgKrDLtehnpdQLNYiv1njnnXcYP348IkK3bt344osvCNCP4J8Yhx1+HAp7/4GrPof253nkMPnFdhZuS2fWhv3M2nCAEruTU9o04eObenNSy4p7udicNh5f8Diz98zWScFDlN1O2nvvkzFuHP6dOhH/7jv46cEp3a4mM7i9KiJzgf+AnsBAoNLEAHwJfIDRLnEsC5VSF1c3ptooKSmJ9957j02bNhEYGMg111zD5MmTGTp0qK9Dq7uUgt8egm1/wEVvQ9cr3bLbEruTpKxCth7IYc2+bNbsO8iqPVmUOJyEBVi47qTmXHdSQoXDWxyik4Ln2dPSSBr5MAXLlxNx9dU0e+pJPdaRh1Q7MYjIC4AZWINxt/BXVX6nlFogIi2re7y6zG63U1hYiNVqpaCggNjYWF+HVLcteANWTYTTH4GTju7zoJQir9hOdqGNnEI7OUU2cgpt5BTZyS0quyy3yM7BghISDxayP7sQpzL2YTULnWPCuGVAC87u1Iy+LRthNR//OdDSSeHRvo9yS5dbPHH2DVr+v8tIeuRhnLl5xIx5lYjLLvN1SPVaTe4YnhGRZzCemr5SRNoopdw1CtgAEVkLJAOPKKU2VrSRiAwDhgEkJBz/MYpZs2Zx4MABN4VniI6O5oILLjjuNnFxcTzyyCMkJCQQGBjI4MGDGazrQWtu9SSY/zK2rteyutV9rF+0i93p+SQeLCApq5C03GJyiuw4DpXwxxDsZyY0wEpYoIWwACsntWxEQuM4EpoE0yYqmM6xYfhbqt73XScFz1JOJxmfjSdt7Fj8WrQg4fPPCWjf3tdh1Xs1fY5hAnAnEAx85KZYVgEtlFJ5InIh8CvQrqINlVLjgHEAffv2PX5J4CMHDx5k6tSp7Nq1i4iICK6++momTZrETTfd5OvQ6py01TNoMvUB1vv14pqVF1G8YikA4YFW4hsF0rJJMP1aNSY80FrmFRZgJcz1HhpgITTAgqWSq//qsDltPPb3Y8zZO4fHTnqMmzvf7LZ9a2BPTyf58VHkL15M2IUXEP3Ci5hDgiv/oXbCapoYHsAYFsMCjMVoZzghSqmcUp9nishHIhKplEo/kf1WdmXvKXPmzKFVq1ZERUUBcMUVV7BkyRKdGKph7b4sps6aycNJD7JZxfNS8BMM69uSXgkRdI0Lp2mo7+qXdVLwrPylS0l69FGcOblEP/88Eddcrcc68qKaJoYdGFfzU5VSD7kjEBGJBlKUUkpE+mFUVWW4Y9++kJCQwNKlSykoKCAwMJC5c+fSt29fX4dVJ6TmFvHijM2sXreWX/2fxeHfiEY3/8qPzVv7OjQAbA4bjy54lLl75/L4SY9zU2ed7N1FORykf/gR6R9/jF+rViSMH09Ahw6+DqvBqWli2AjsA+4QkTeUUidV9gMR+Q4YBESKSCLwLGAFUEp9AlwF3CsidozhNq5TStXKaqKq6N+/P1dddRW9e/fGYrHQq1cvhg3Ts3RVZvamFB7/aR0U5/JHo/do4nAid04jLKr2JYVR/UZxY6cbfR1SvWFLSSH5kUcpWL6c8MsuI/rp0ZiCddWRL0hNyl4RuRXjeYaDwEGl1AZ3B1ZVffv2VStWrCizbPPmzXTq1MlHEVVPXYrVk5RSfPL3Tl77fQvdYoL5Lux9QvbOh5umQJuzfB0eYCSFR/5+hHn75umk4GZ5CxaQ/PgonMXFRD/ztO515AUislIpVWE1Rk3vGPYppeaJSAyQWvPQNM1ICi/O2MyExbu4pEcs7zT6CcvSOXDhm7UqKTz898PM3zdfJwU3UjYbaWPHkjH+c/w7dCDunbfxb1077g4bspomhvNF5D+M0VX3YDRGa1q1KaUY8/sWJizexW2ntuTpuNWYpr0Pfe+AfrVjLuTSSeGJfk9wQ6cbfB1SvWBLSiJp5MMUrl1LxHXX0mzUKP3AWi1R08QQATwOPIbRbbXWUUrV+l4MdbgJxW3GLdjJp3/v5KaTE3imew4y8SFodQZc8JqvQwOMpDDy75H8te8vnuz/JNd31PN3u0POn3+yf/TT4HQS987bhPmo96BWsZomhheAjkqprSLicGdA7hAQEEBGRgZNmjSptclBKUVGRkaDHjtp3pYUxvy+hYu6x/DCwFBk/GUQHg9Xfwlm34+SWeIo4eG/H9ZJwY2cBQWkvPoqWT9OIaBrV+Lefgu/Sh5S1byvSolBRMxAIvC0Umq8UirR9R2l1CgPxlcj8fHxJCYmkpaW5utQjisgIID4+Hhfh+ET21NzeeC7NXSOCePNS1pj+vpCYwa267+HIN9PxVjsKGbkXyNZkLiAp/o/xXUdr/N1SHVe4YaNJD/yCCV79tBk2DCi/jdcD5NdS1UpMSilHCKyAWjj4Xjcwmq10qpVK1+HoR1Dkc3B/d+sxt9i4rObexM4405I2wI3/ghRvh/uoNhRzIj5I1ictJinT36aazpc4+uQ6jTldJI5YQKpY9/D0rgxCV9+SXD/fr4OSzuO6lQlBQGPici5GGMZASil1BD3h6XVZ6/O3MzWlFy+vO0kYle8DltnwgVvQNuzfR0aRfYiRswfwT/J//DcgOe4sr17RnBtqGwpKSSPGkXBP0sJHTyYmBeexxwR4euwtEpUJzEMcL33dr0AdOupVi1zN6fw1T97uP3UVgwqnAOL34W+t9eKHkiF9kIemPcA/+7/l+dPeZ7L213u65DqtNw5c9j/1GicJSXEvPQi4VdeWWvb/LSyqpMYdN2MdkJSc4p4dMo6OsWEMaprFkwaAa0GwgWvg48LjAJbAQ/Me4BlB5bx4qkvMqStvhGuKWdhISljXiPr++8J6NKF2DffwF9X7dYplSYGETnUZaDCu4NS67NKD4SnaaU5nYqRP6yloMTOxxdH4vfjxa4eSF/5vAdSga2A4fOGszJlJS+f9jKXtLnEp/HUZYVr15L8+CijgfnOO4h64AHEz8/XYWnVVJU7hq8wksLxLukUxgxtx5udTWvAPl+0i0Xb03ntkta0/PPOWtMDKd+Wz31z7mNN2hpeOe0VLmp9kU/jqatUSQlpH31ExrjPsEQ3I+GLLwg+ub+vw9JqqNLEoJQ60xuBaPXXxuRsXv9jC+d1iuKavS9A6ia4cYrPeyDlleRx39z7WJe2jtdOf43zW53v03jqqqKt/5E8ahTFmzcTfsUVNHvyCcwhIb4OSzsBNX3ATdOqpLDEwYjJa2gU5Me7Tacj/8402hR83AMppySHe+fcy8b0jbw28DXOa3meT+Opi5TDQeYXX5A29j1MYWHEf/QhoWfVjrGttBOjE4PmUa/O2sz21Dz+ODOJwH/egz63QT/fDj+eWZTJ3bPvZnvWdt464y3ObuH7brJ1TcnevSSPeoLCVasIPfdcop9/Dktj3z+YqLmHTgyax8zdnMLEf/bwXM9cOix7ClqeDhe+4dMeSCn5Kdw1+y6S85J5/6z3OS3uNJ/FUhcpp5OsH34g5fU3ELOZ2NdfI+ySS3Q31HpGJwbNI/ZlFjDyh7UMbFbErftGQ1gcXDPRpz2QEnMTufPPO8kqzuKTcz6hb7SeUa86SvbuZf/opylYtozgU04h5uWXsMbE+DoszQPcNzN6JURkgoikuobWqGi9iMh7IrJdRNaJSO+KttNqvyKbg3smrSRAFTLe8iZiL4YbfNsDaWf2Tm79/VZyS3IZP3i8TgrVoBwOMr74kp2XDqFo0yaiX3yB5p+P10mhHvPmHcOXwAccu0vrBRjzSLcD+gMfu961OkQpxVO/bGBTchbL2nyFX/KhMZB8N2/vlswt3D37bgAmnDeBDo31HMJVVbxtG8mjR1O0dh0hZ55J9HPPYm3WzNdhaR7mtcSglFogIi2Ps8kQYKJrnuelIhIhIjFKqf3eiVBzh/fmbuenVYn82OZPopLmusZAOsdn8axNW8u9c+4lyBLE+MHjaRne0mex1CXKZiNj/HjSP/oYU3AwsW++SdhFF+q2hAaiNrUxxAH7Sn1PdC07KjGIyDBgGECCHsu91vh++V7emfMfY1qv56SkiT6fhW1J8hIenP8gkYGRfDb4M+JC4nwWS11SuGEj+0ePpnjLFsIuvIBmTz2FpUkTX4eleVFtSgwVXYpUOAyHUmocMA6gb9++eiC/WuCnlYk88fN67kxI4doDbx6Zhc1HV5izds3iyUVP0iq8FZ+c8wlNg5r6JI66xJGXR9q7Yzn47bdYmjQh/sMPCD1bd+VtiGpTYkgEmpf6Hs+R4b21Wuy7ZXt58pf1DGlh46nsl5BGLeAa342B9M3mbxizbAy9m/bm/bPfJ8wvzCdx1BVKKXL/+IOUl1/Bnp5Oo+uvJ+qhBzGHhvo6NM1HalNimAYMF5HJGI3O2bp9oXZzOBWv/b6FcQt2cl7bIN4uHoUouzEGUmAjr8ejlOL91e/z2frPOLP5mbw+8HUCLA136tSqKNm3jwMvvkj+goX4d+5E/EcfEtitm6/D0nzMa4lBRL4DBgGRIpIIPAtYAZRSnwAzgQuB7UABcJu3YtOq70B2EY9OWcvCbenc0j+e5/JfxJS0DW7+GSLbej0eu9POi0tf5OdtP3NluysZffJoLKbadN1Tu6iSEjK++JL0jz5CzGaaPfkEjW64AbHoPzPNu72SjjuTuqs30v1eCkerIadTMWVVIi//tpkSu5NXr+jG9Zkfw9rZcPE70HqQ12Mqshfx2ILHmL9vPsO6D2N4z+G698xx5C1eTMorr1KyYwehgwfT7KkndRdUrQx9eaBViVKKv/9L440/trIxOYfeCRG8dU1PWu2YBEs/gv73GDOxeVlmUSYj5o1gbdpaRvUbxY2dbvR6DHVFyd69pLz2Onlz52JNSCD+k48JHTTI12FptZBODNpxZeaXMHP9fr5aspttqXnERQQy9rqeXNI9FtN/M2HW49DhIjjvFa/Htit7F/fNuY+0wjTePONNBrcc7PUY6gJnfj7pn44j84svwGol6uGRNL71Vkx6Ah3tGHRi0MpwOBWb9+ewdGcGf/+XxpIdGTicis4xYbx9TQ8u7h6Ln8UEiSthyh0Q1xuuHA8ms1fjXH5gOQ/OfxCLycLn531Oj6geXj1+XaCUImf6dFLffAt7airhQy4lauTDWJvprrva8enE0IAdzC9ha0ou21Jy2ZqSy38peWzen0NukR2A1pHB3D2wNRd1j6FzTNiRevvMXfDtNRDS1OiB5Bfk1bin75jOM0ueoXlocz48+0Oahzav/EcNTMGq1aS+/jqFa9YQ0KULcWPfJahXL1+HpdUROjE0ADlFNra5Cv6tB3LZlprL1gN5pOcVH94mNMBCh2ahXNojln6tGtO/VROiwyvo6lmQCd9cDcoBN/0EIVFeOw+lFB+v/ZiP135Mv+h+vD3obcL9w712/LqgeOcu0t55m9zZczBHRRLz0ouEX3EFYvLaeJlaPaATQz2jlGLLgVyW785k1Z6DrNqbxd7MgsPrA61m2jcLYVCHKDo0C6V9dCjtm4UQHRZQeU8eWxFMvgGy9sItUyGynYfP5ohCeyHPLn6WWbtncWmbS3luwHNYfTiEd21jT08n7cMPyfrhR0z+/kQ+8D+aDB2KKci7d3Na/aATQz1gczhZuC2NPzemMH9rKik5xp1A01B/eic04tqTmtMxOpT2zUKJiwjEZKpBV06nA36+C/b+A1dNgBYD3HwWx7Y/bz8j5o9gS+YWRvQewR1d79DdUV2c+flkfPklmZ9PwFlSQqNrryXy/vv02EbaCdGJoQ7bnprHd8v2MnVNEul5JYT6Wzi9fSSDOjTllDZNiIsIdE8BqhTMeBA2T4PzXoWuV574PqtoxYEVPPz3w5Q4Svjg7A8YGD/Qa8euzZxFRWR9/z3pn43HkZ5O6ODBRD30IP6tWvk6NK0e0ImhDlq19yCf/LWD2ZtTsJiEszs246o+8QxsH2X0GHK3Oc/Bqokw8FEYcJ/7938M32/5njHLxhAfGs97Z71Hq3Bd6DmLi8n6cQoZn36KPS2NoAEnE/X+e7phWXMrnRjqkN3p+YyZtYXfNx4gPNDK/85syy2ntCQyxN9zB108Fha/azy8duZTnjtOKcWOYl7991V+2vYTp8edzpiBYxr8QHiqpISsn38m/ZNPsR84QFDfvsS+9SbB/fr5OjStHtKJoQ4otjt4f+52Pl2wA6vZxMPntuf201oR7O/hv75VX8PsZ6DL5XDhm14ZQntfzj5G/j2SLZlbuLPbnQzvORyzl5+RqE2cxcVk//IL6ePGYU/eT2CvXsSOeZWg/v11O4vmMTox1HIbkrJ55Me1bDmQyxW94xh1fkeahnlhxNANP8P0B6DNWXD5OK88wDZnzxyeXvw0JjHxwVkfcEbzMzx+zNrKkZfHwe++I/OriTjS0wno0Z2Y518g+LRTdULQPE4nhlps8rK9PDN1IxFBVj6/tS9nd/LSQGebpsJPd0Lz/nDN12Dx7NAJNoeNt1e+zaTNk+gW2Y03z3iT2JBYjx6ztrKnp5M58WsOfvcdztxcgk85hSZvvklQ/346IWheoxNDLVRid/LstA18t2wfp7eL5L3retEo2Evj2myeAVNuh/i+cOOP4B/i0cPty9nHqIWjWJe+jhs73cjDfR5ukM8nFO/cRebXE8n++RdUSQmhgwfT5K67COzaxdehaQ2QTgy1TEGJnXsmrWLBf2ncN6gNDw/ugLkmzx3UxNZZ8ONQiOkJN04Bf8/N4KWU4tftvzJm2RjMYuatM95qcIPgKaXIX7SYzIkTyV+4ELFaCbv0EprccSf+rXUPLM13dGKoRQ7ml3Dbl8tZl5jFa1d249qTErx38P/+gO9vhuhuxmQ7AZ7rBZRVlMXz/zzPnL1zOCn6JF4+9WViQmI8drzaxllQQPbUqWR+PYmSnTsxR0US+b/hNLr2WiyRkb4OT9O8mxhE5HxgLGAGxiulxpRbPwiYCuxyLfpZKfWCN2P0lZwiGzdP+Jf/UvL4+KY+nNcl2nsH3/AT/DwMmnWFm3+BAM+NP7QkaQmjF4/mYPFBRvYZya1dbsUkDWMcn+Jt2zj4449k/zoVZ04OAV27Evv6a4Sdfz6ih8DWahFvTu1pBj4EzgUSgeUiMk0ptancpguVUhd7K67aoKDEzu1fLGfrgVzG3dyXMzt6cVjklV/C9AchYQDcMNljSSGnJIe3V7zNT9t+ok14Gz465yM6Nu7okWPVJs7CQnJm/U7WDz9QuGYNYrUSeu65NLrpJgJ79dQNylqt5M07hn7AdqXUTgARmQwMAconhgbF5nBy99crWbX3IB/c0Nu7SWHJ+/DnaGh7Llwz0WPDZ8/bO4+Xlr5EZlEmt3e9nXt73EuAxQtdbn1EKUXRpk1k//QT2dOm48zLw69VK5o+/jjhlw3B0qiRr0PUtOPyZmKIA/aV+p4I9K9guwEishZIBh5RSm0sv4GIDAOGASQkeLEe3s2UUjwzdSMLt6Xz+pXdubCbl+rZnU6Y9yIsehs6XwZXfOaRLqkZhRmMWTaG33f/TodGHXj/7Pfp0qT+9rKxJSWRPeM3sqdPo2T7DsTfn7DzzyPi6qsJ7NNH3x1odYY3E0NF/ytUue+rgBZKqTwRuRD4FThqbGel1DhgHEDfvn3L76PO+HLJbr5btpd7B7XhmpO8NNmMrQim3me0K/QZChe97faH1+xOOz/+9yPvr36fInsR/+v1P27rehtWU/3rhurIzibn9z/Inj6NwhUrAQjs04fo554j7ILzMYfr+SK0usebiSERKF36xWPcFRymlMop9XmmiHwkIpFKqXQvxeg1f/+XxoszNjG4czMeHdzBOwctyDTmU9j7D5zzHJz6oNuHuVidupqXl77M1oNb6R/Tnyf7PUnriNZuPYav2Q8eJG/ePHL/nE3+kiUomw2/Vq2IenAEYRdfjF98vK9D1LQT4s3EsBxoJyKtgCTgOuCG0huISDSQopRSItIPMAEZXozRK5KzChkxeTXtm4XyzrU9azY/QnWlbzem48xONOZTcPPQ2Sn5KYxdNZbpO6cTHRzNW2e8xbktzq031Se21FRy58whd/ZsCpYtB4cDa2wsjW64gbBLLiGgS+d6c66a5rXEoJSyi8hw4A+M7qoTlFIbReQe1/pPgKuAe0XEDhQC1yml6mxVUUXsDicPfLcam93Jxzf18fxAeABbfoNf7gGzFW6dBgknu23XOSU5TFg/gUmbJ+FUTu7qdhd3druTIGvdnjlMOZ0UbdpM3oK/yf97AYXr1oFS+LVqRZM77yR08LkEdNbJQPMepRS2YgfFBXZKCu0UF9oJCvMjoqn7/6959TkGpdRMYGa5ZZ+U+vwB8IE3Y/K2d+b8x4o9Bxl7XU9aRQZ79mBOB/z1Kix4A2J7GeMeRbinLaPYUczkLZMZt24cOSU5XNT6Iob3HE58aN2tRnFkZ5O/ZAl5fy8gb+FCHBkZIEJAt25E/m84YYMH49+2ra/D1Ooop8NJSaGD4kKb8V5go7jQVcgX2A9/Lin1ufT6kkI75S+Te5+XwIDL3f9vUj/57EUL/kvjo792cN1JzRnSM86zB8tNgV/vhR1zodfNxrDZ1hPvIlpkL+LnbT8zYcMEUgpSODX2VEb0HkGnJp3cELR3OQsKKFi5ioJl/5K/bBlFGzaCw4EpPJyQ004j5IyBBJ92GpbGjX0dquZjSikcNmeFBXnpK/iScstLF/D2Ykelx/ELMOMXZME/0IpfoJmQRgH4x1rwC7TgH+R6DzzyPTwq0CPnqxODl6TmFDHyhzW0axrCs5d4uMvm1t9h6v1QkgeXjDV6H52gfFs+P2z9ga82fkVGUQa9mvbipdNe4uQY91VLeZozP5/CdevI//dfCv5dRuH69WC3g8VCYLduNBl2FyGnDySwR3fE3HDngKiPlFNRUmxcpZcUOigptFFc6KDkGFftZQp713en4/i12iaTuAr1IwV3cHgQfoGWo5aX/nyosLcGWLzT3lgFOjF4gcOpGDF5DXnFdr6762QC/TxU6JQUGBPrLP8MmnWDK8dD0xN7uji1IJUftv7A5K2TyS7O5uSYk3mj+xv0bda3VtevK6Uo2b2bwrVrKVyzhsK16yjeutV4hsNsJqBrF5rcdhtB/fsT1LsXpqC63SZS3znszgqrVqp8BV9kP7pzfDkWP5NRYAdZ8Q80ExDiR3hUoOu7Bb9AM/6HP5e7gg+yYLGaavX/ierQicELPpi3nX92ZvD6Vd1p18xDI5bumAczHoKDu2HAcDj7GbDUbMpPpRRr09by7eZvmb1nNg7l4Iz4M7iz+530iOrh3rjdQCmFPTmZoi1bKNq0mcIN6ylasxZHdjYAppAQArt3J/Seewjs2YPA3n0wh3i4fUcrw+FwGoV0QdkC+1A9e3FB6WV244q+VGFvL3Ee/wDC4QL7UGEd1iQAv/iKr87LX8H7BVowmxvGmF1VoRODh/2zI4Oxc//jsp6xXN3HAw2z+Rnw51Ow9jto0hZunQGtTq/RrrKLs/l91+/8vP1nNmVsItQayvWdruf6DtfTPMxLD+BVQtlsFO/aRbErCRRt2ULx5s2HkwAi+LVpTci55xDYowdBPXvi16YNYtL/6U9E6YbTowt3o4A//L38uirUr5euhjlUUAdH+Je5gvcLdL2Xu2r3D7Rg9TcjtaQapj7QicGD0nKLeWDyalpGBvPy5d3ce5tpLzGqjP5+DUry4fRHYOCj1W5gdjgdLElewtQdU5m/dz4lzhLaNWrH6P6juaTNJT7rduosLqZk1y6Kd+ygZMcOirfvoHjnDkp27zHaBQDx88O/fXtCBw8moHMnAjp1wr99e10tVAGnUxk9XgrLX7W7CvoKruRLb2srOn7BLkKpgt2KX6CFRs2CyhX21iP160GHCnVjmcWv/lTDVIdSCuxOnCVOVInDeNmMz0eWOVE2491Z4jj8WZU4COjQmKAeUW6PSycGD3E6FSN/WENOoY2Jt/dz3/MKSsHmaTDnOcjcCW3OhvNehqZV7xVkc9pYtn8Zs/fMZv6++WQWZRLhH8HVHa5mSJshdGzc0Sv/SZXdji05mZK9+7Dt20vJ3n2U7N5N8c4d2PYlGu0BACYT1ubx+LdpS+iZZ+Hftg3+nTrh37o1YmkY/4TLNp7aDxfmJaWu2Cu6Ui85dDVfScFeuirmUOEdHlW+YD+yzj/oSALwD3JdsdfDgl0pBQ51pKAuVSgfKsSdJaWWHS7UyxboqsSB01Z2m0OfK2v7KE+sJsTPhFjNWJt55iKoYfyv8oEP529n4bZ0Xr2iG51i3DDpjdNhJIQFb0LKBojqBDf+BO3OqdLPs4uz+Xf/v/yd+Dfz980ntySXIEsQA+MHcl7L8xgYPxA/s/sH0nPk5WNLTsK2bx8le/ZSsm8vtr37KNm3D1tSEjiOFFji749fQnMCOnUm/KKL8W/bBr82bfBr2RKTf83aS2qLox5OOlxw2yg6atmRK/nSyysrQEp3dfQPMurY/eMtR13J+wcdefkdWl6Hq2KUw1nBVXbpAvjQlXapK/AyhXq5bW2OMlfqVNK8cRSLIFYzJj8T4mc2XlYTpkALpjA/1zKjYBfXNibXNofXHWOZWExe+XvSicED/th4gLfnGO0K153o4HjFebD+R1j6EaT/B03aweWfQterwHzsvz6b08bG9I0sTl7MkuQlbEjfgFM5CfUL5czmZ3JOwjmcEncK/uaaF7jK4cCeloYteT+2/cnY9+83PicnY9u/H9v+/Thzcsr8xhQWhl9CAoFduxB2wQX4JTTHLyEBa0IClqioWt0WYLcZBXtxftlG06MaS8tdtR/qIqmcxy/Zrf7mMlfmh/uwV3TV7qpnP7TML7D2dHUsTzlV2Svt0gVy+WqSkqOvqo+6Ui+dBGwOqKQb6VFMUnHh62/GFOKHxc90dKFsdW17qIC2Hin0TaULeasZMdfOv4fq0InBzTYmZ/Pg5DV0j49gzJXda357nboZVnxhNCoX5xhTbl71BXQeUuFoqNnF2axNW8ua1DWsTl3NhvQNFDmKMImJrpFdGdZ9GKfGnkrXyK5YTJX/tSulcObmYjtwAPuBA67C30gAtuRk7Mn7saWmHq7vP8QUHo41JgZrbCxBffpgjY3BEhODX0ICfs2bY46IqNmfhxuUvmo/XP1S0edDhXx+qeWFdhy24186lu7uGBBkISjcj0bRQaV6wViPFO7BpapkXA8zmXzUK0Y5Fcp+rPrsUlfgrnXHrQ+3Octt6wB7detKKFP4mkoVxqYg/yOfy1x5l7rSLlegy6GC/tA+LbX34qO20InBjXan53P7l8uJCLLy2c19CLBW83mFg3uM4bA3/Awp68HsB10uh5PuhPiTDo+Eml2czdbMrWzO3MyWzC1sztjMjuwdAJjFTMfGHbmq/VX0atqL/jH9CfcvO/SzUgpHVpZR4B84gD0lxZUAUrClHHpPQRUUlI3PYsHarBnWmBgC+/YhLCYWa2ws1tgYrDExWKJjPN4NVDkVJUVHrtSLDvWGcX0+fAXv+lxUqtAvKbDjPN5Vu6ue/fAVeZCF4PDgI9+Dy16pH9rmcHdHDxU4RgOlOqrhsaJqj8P12dWoD1eVJLyKlL+aPlT4mkL9jq4mOWpb1/altild2GORetleURmn04nD4TjqZbfbK1zucDho1KgRTZu6f3Ivqetj1PXt21etWLHC12GwN6OA68b9Q5HdyXd3nUyH6Co8r+CwQeJy2DYbts+GA+uN5fH9oOuVFHS8kD2OfHZl72J3zm7+O/gfWzK3kJSXdHgXUYFRdGzckR5RPejVtBddGnfGP7cI24EU7Af2G+8pB1zfD2BLMd5VSUnZWMxmLE2bYm3WDEt09JH36GZYmkVjjYvFEhnplieCnQ5jaAHjirzsVXlxge2oK/Xy2xyvrl1MUqoO3bhy9ytXiAccqmsPNj4fLtwDLDWuvz10xe20HX3lXXZ5RdUnldSH2xzVbqDEYipbxVGmSqRsYVz+qrx8QV2+oPdWPbc7KaWqVeD6Yl1NyuJTTz2Vc889t0Z/JiKyUinVt8J1OjGcuJV7Mhk2cSUOpfj2zpPpHHuMxuaSAkheBfv+hb3/wt6lOIuzSbH4sSuuG7uatmF3SBN2F2ewK3sXKQUph38qCAnB8fQyt6KLiqZ1cRjRBX74Z+SVKfBtqalgs5U9rtXqKuibYW0WXfY9OhpLs2gskU2qVegfrm8/brWMrcJtbJX0aTdbTGUK94o/l/5+ZPmxescohyp3RV1xtUf5Av2YV9nlqk+opP3g6JMs10BZvuGxdAF9VJVIJQ2XVu82JCulylzt1rYC1+Fw4HRW/66oMiKC2Ww+/LJYLGW+e2NdSEgIYWE169yiE4OH2BxOPv5rB+/N3UbzxkF8fmtfWkeFGCuLcyFlk1EldGADBQfWsjtjM7vMwm6rld0hjdgVEMQeVUyR04YoRVg+JBQE0rEkkpZFIcTkW2mcowjKKsScloUjPb1MLx4wevJYo6PLXN0fKfCNd3OjRkc16h6qby/d3bEo/0gXyKJyhXtJuWqZSuvb/c0ElCrA/QItru+HqmSO9JQ5XO/uZ8JqFsyK4xfIFVyVV1YfXv0GSspVhxzd8Fhhr5FjVJEcdVVejfaEQ4VubSxwD708wZcFblXWmWpxR4mq0InBzfKK7Uxfm8yH87eTeLCAoV3MPNLdRlDmJlJSVrErYyu7itPZZbWy22pll58fB51CZA5E5iia5ggtC0NoXhBIVK4QllmMX0YOYivbkCtBQVjLF/ilCn5zVFMcASGuIXyPfYV+rGqZygYFK9O1MdBCQKCFwAAz/v5m/P3MBFhNWK2C1WLCahKsJsEixuxK2I/uiXLM/t0lTrBX84pOOKrao0oNj36llllNYDXhNIPTonCaQJkUDpPCgROn0+nzAtfhcNSoiqEyJpOp1ha4hwrdhtjO4E06MZyg/GI765Oy2bB7P6nbV5OTuIrGlu00Cj5AsfUge01OUm1WCgrNhOW6EkC2IjrXREyelcY5TgJyy9Xpm0xYmjbFHBML0c1RUXE4mzTDGRaJI7gRdv9QSpyWsg8wlSvwy4/PbgLMYsyCZBGwmISAAAuB/iajMLea8LOa8bMIVrMJixksIod/Y1IKkxPE4TxSsJduoKzuPxWroCyC8hOcfuC0CsosKKvCacFVILvezQplMt4donCaXC8xXg5x4sRVYOPEoY5uqKtJYeyJKgaTyVRrC9xDL13oasdLDF7tlSQi5wNjMcqh8UqpMeXWi2v9hUABMFQptcqbMTqcim0pOWzdsoHEXYvJPLgaZ8kerCUHcRbZKSmwEJQvhGRDeA5E5lg5LdcCpiDslkDsliBsliDsQeGoJjE4m0WRGtIIR0AodkswTvHH4TTjsIGj2AElTsw2MO8Hy35XwS7ZWMjG7CrcA60mws2C1SxYTILZAqYwM6YwEzgVTqcdp8OBUx0qQMsXogpnkRNHkaIIJ/k4jcLWpFCuQtlpxiiIzcaVs1OcOP3BGeDEgcIprnecxn6V8e5wOnEoo4A9XFg7S1UxKKDY9XKTigq7ipZZrVafFMZ1vYpB07yWGETEDHwInAskAstFZJpSalOpzS4A2rle/YGPXe8ek5SSzIrFv5CyeRm2tDSs+YJ/cTBWWxixjlDaqLaYTd1x+AWiLP44Lf44/fxwxlqhuZkck4kcUYiA4ARRIE7AeFdy6D0H5cw6XLg6LU6cVuNzcfmCXA4VwAqnybXOdcV8+Gr50OOYJtfrRP4mHa4XVLEA9MPPR1e/utDVNM/z5h1DP2C7UmongIhMBoYApRPDEGCia57npSISISIxSqn97g5m0qiXSbGaUIJR4Ep7nOHtcIQbBbCS0vUm7rnkFQWCUY4f+my8K0CBUsZnp8J4Dl8hduX6jcLk2lZQGOGpw78XlJGLUK79Hr3sUL2T8XvXslL7MNZVTgF216tOqOPVpTXT8M65rleL10T3s8/npEuvdPt+vZkY4oB9pb4ncvTdQEXbxAFlEoOIDAOGASQkJNQoGL+gAIILnUcKTwGTCCaTICYx7gAETCZl9FA5XIofuitwFdAmVwGPYJJDn42ubCZw3UkcWeZrvo/AR2rBn7231YZ/b5pnhTaJ9Mh+vZkYKvpXWj7FV2UblFLjgHFgND7XJJhrnnm4Jj/TNE2r97xZYZsIlB5RLh5IrsE2mqZpmgd5MzEsB9qJSCsR8QOuA6aV22YacIsYTgayPdG+oGmaph2b16qSlFJ2ERkO/IHRXXWCUmqjiNzjWv8JMBOjq+p2jO6qt3krPk3TNM3g1ecYlFIzMQr/0ss+KfVZAfd7MyZN0zStLN0pXNM0TStDJwZN0zStDJ0YNE3TtDJ0YtA0TdPKqPOjq4pIGrCnhj+PBNLdGE5doM+5YdDn3DCcyDm3UEpFVbSizieGEyEiK4417Gx9pc+5YdDn3DB46px1VZKmaZpWhk4MmqZpWhkNPTGM83UAPqDPuWHQ59wweOScG3Qbg6Zpmna0hn7HoGmappWjE4OmaZpWRoNIDCJyvohsFZHtIjKqgvUiIu+51q8Tkd6+iNOdqnDON7rOdZ2ILBGRHr6I050qO+dS250kIg4Rucqb8XlCVc5ZRAaJyBoR2Sgif3s7Rnerwr/tcBGZLiJrXedcp0dpFpEJIpIqIhuOsd795ZdSql6/MIb43gG0BvyAtUDncttcCMzCmEHuZOBfX8fthXM+BWjk+nxBQzjnUtvNwxjl9ypfx+2Fv+cIjHnVE1zfm/o6bi+c85PAa67PUUAm4Ofr2E/gnAcCvYENx1jv9vKrIdwx9AO2K6V2KqVKgMnAkHLbDAEmKsNSIEJEYrwdqBtVes5KqSVKqYOur0sxZsury6ry9wzwP+AnINWbwXlIVc75BuBnpdReAKVUXT/vqpyzAkLFmPQ6BCMx2L0bpvsopRZgnMOxuL38agiJIQ7YV+p7omtZdbepS6p7PndgXHHUZZWes4jEAZcDn1A/VOXvuT3QSET+EpGVInKL16LzjKqc8wdAJ4xpgdcDI5RSTu+E5xNuL7+8OlGPj0gFy8r30a3KNnVJlc9HRM7ESAyneTQiz6vKOb8LPK6UchgXk3VeVc7ZAvQBzgYCgX9EZKlS6j9PB+chVTnn84A1wFlAG2C2iCxUSuV4ODZfcXv51RASQyLQvNT3eIwriepuU5dU6XxEpDswHrhAKZXhpdg8pSrn3BeY7EoKkcCFImJXSv3qlQjdr6r/ttOVUvlAvogsAHoAdTUxVOWcbwPGKKMCfruI7AI6Asu8E6LXub38aghVScuBdiLSSkT8gOuAaeW2mQbc4mrdPxnIVkrt93agblTpOYtIAvAzcHMdvnosrdJzVkq1Ukq1VEq1BKYA99XhpABV+7c9FThdRCwiEgT0BzZ7OU53qso578W4Q0JEmgEdgJ1ejdK73F5+1fs7BqWUXUSGA39g9GiYoJTaKCL3uNZ/gtFD5UJgO1CAccVRZ1XxnJ8BmgAfua6g7aoOj0xZxXOuV6pyzkqpzSLyO7AOcALjlVIVdnusC6r49/wi8KWIrMeoZnlcKVVnh+MWke+AQUCkiCQCzwJW8Fz5pYfE0DRN08poCFVJmqZpWjXoxKBpmqaVoRODpmmaVoZODJqmaVoZOjFomqZpZejEoGmliEiEiNxX6nusiEzx0LEuE5FnKtnmTRE5yxPH17Rj0d1VNa0UEWkJzFBKdfXCsZYAlx6vj72ItAA+U0oN9nQ8mnaIvmPQtLLGAG1c8xe8ISItD42DLyJDReRX11j/u0RkuIiMFJHVIrJURBq7tmsjIr+7Bq1bKCIdyx9ERNoDxUqpdBEJde3P6loXJiK7RcSqlNoDNBGRaC/+GWgNnE4MmlbWKGCHUqqnUurRCtZ3xRjKuh/wMlCglOoF/AMcGrl0HPA/pVQf4BHgowr2cyqwCkAplQv8BVzkWncd8JNSyub6vsq1vaZ5Rb0fEkPT3Gy+qyDPFZFsYLpr+Xqgu4iEYEyC9GOpEVz9K9hPDJBW6vt44DHgV4whDe4qtS4ViHXXCWhaZXRi0LTqKS712VnquxPj/5MJyFJK9axkP4VA+KEvSqnFrmqrMwBzufGMAlzba5pX6KokTSsrFwit6Y9dY/7vEpGr4fB8vBXNp70ZaFtu2UTgO+CLcsvbA3V24Dut7tGJQdNKcc1LsVhENojIGzXczY3AHSKyFthIxVOMLgB6SdkZg74BGmEkBwBcDdJtgRU1jEXTqk13V9U0HxGRscB0pdQc1/ergCFKqZtLbXM50Fsp9bSPwtQaIN3GoGm+8wrGxDmIyPvABRjj6pdmAd7yclxaA6fvGDRN07QydBuDpmmaVoZODJqmaVoZOjFomqZpZejEoGmappWhE4OmaZpWxv8BxlSK+txCmxMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -125,9 +125,25 @@ "execution_count": 9, "metadata": {}, "outputs": [ + { + "ename": "AttributeError", + "evalue": "'numpy.ndarray' object has no attribute 'mid'", + "output_type": "error", + "traceback": [ + "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", + "\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)", + "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0mfig\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0max\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubplots\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mswiftdiff\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'dr'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msel\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mid\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mtpidx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mline\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"time (y)\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0max\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0max\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 3\u001b[0m \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_ylabel\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"$|\\mathbf{r}_{swiftest} - \\mathbf{r}_{swifter}|$\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 4\u001b[0m \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_title\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Helio integrator \\n Test Particles only\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 5\u001b[0m \u001b[0mfig\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msavefig\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"helio_swifter_comparison-tp-rmag.png\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfacecolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'white'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtransparent\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdpi\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m300\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;32m~/.conda/envs/cent7/2020.02-py37/swiftestOOF/lib/python3.7/site-packages/xarray/plot/plot.py\u001b[0m in \u001b[0;36mline\u001b[0;34m(self, *args, **kwargs)\u001b[0m\n\u001b[1;32m 457\u001b[0m \u001b[0;34m@\u001b[0m\u001b[0mfunctools\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mwraps\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mline\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 458\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mline\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 459\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mline\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_da\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 460\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 461\u001b[0m \u001b[0;34m@\u001b[0m\u001b[0mfunctools\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mwraps\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstep\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;32m~/.conda/envs/cent7/2020.02-py37/swiftestOOF/lib/python3.7/site-packages/xarray/plot/plot.py\u001b[0m in \u001b[0;36mline\u001b[0;34m(darray, row, col, figsize, aspect, size, ax, hue, x, y, xincrease, yincrease, xscale, yscale, xticks, yticks, xlim, ylim, add_legend, _labels, *args, **kwargs)\u001b[0m\n\u001b[1;32m 296\u001b[0m \u001b[0;31m# Remove pd.Intervals if contained in xplt.values and/or yplt.values.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 297\u001b[0m xplt_val, yplt_val, x_suffix, y_suffix, kwargs = _resolve_intervals_1dplot(\n\u001b[0;32m--> 298\u001b[0;31m \u001b[0mxplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0myplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkwargs\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 299\u001b[0m )\n\u001b[1;32m 300\u001b[0m \u001b[0mxlabel\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlabel_from_attrs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mxplt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mextra\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mx_suffix\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;32m~/.conda/envs/cent7/2020.02-py37/swiftestOOF/lib/python3.7/site-packages/xarray/plot/utils.py\u001b[0m in \u001b[0;36m_resolve_intervals_1dplot\u001b[0;34m(xval, yval, kwargs)\u001b[0m\n\u001b[1;32m 544\u001b[0m \u001b[0mx_suffix\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m\"_center\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 545\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0m_valid_other_type\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0myval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mpd\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mInterval\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 546\u001b[0;31m \u001b[0myval\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_interval_to_mid_points\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0myval\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 547\u001b[0m \u001b[0my_suffix\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m\"_center\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 548\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;32m~/.conda/envs/cent7/2020.02-py37/swiftestOOF/lib/python3.7/site-packages/xarray/plot/utils.py\u001b[0m in \u001b[0;36m_interval_to_mid_points\u001b[0;34m(array)\u001b[0m\n\u001b[1;32m 478\u001b[0m \"\"\"\n\u001b[1;32m 479\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 480\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marray\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmid\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mx\u001b[0m \u001b[0;32min\u001b[0m \u001b[0marray\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 481\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 482\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;32m~/.conda/envs/cent7/2020.02-py37/swiftestOOF/lib/python3.7/site-packages/xarray/plot/utils.py\u001b[0m in \u001b[0;36m\u001b[0;34m(.0)\u001b[0m\n\u001b[1;32m 478\u001b[0m \"\"\"\n\u001b[1;32m 479\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 480\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marray\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmid\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mx\u001b[0m \u001b[0;32min\u001b[0m \u001b[0marray\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 481\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 482\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n", + "\u001b[0;31mAttributeError\u001b[0m: 'numpy.ndarray' object has no attribute 'mid'" + ] + }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAElCAYAAADDUxRwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAiMElEQVR4nO3dfZxWdZ3/8ddbRCeBRAUVHHFQMSVN1FlvysW7KNGU/FX+pDLvyqylbS0XJmvzrhI3+9W66rqablCtbGoitqQixE/TNLFQRGJFEJkYFVC8DRX97B/nYBfjNTPXdc33mplr5v18PObBdc75nnM+B/R6z/fcfI8iAjMzs87aorsLMDOz3sGBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8WsCEkXSvpZ/nmEpFck9atgO9dI+qf0FZr1PA4U65UkPSXpw63mnS7pt+VuKyKejoiBEfFWBeueExGXlNJW0k8kfafcfaRS6d+P2SYOFLM+QNKWvWEf1rM5UKzPkjRc0i2S1khaIenv22jXICk2fWHm682S9LykZZK+0M4+3ul1SDpSUrOkr0t6TlKLpDPyZWcDnwEm56fXbu+oRknvkTRN0guSlkiaLKm5YPlTkqZIehR4VdKWkpokPSnpZUmPSzopb7sPcA1wWL7/9fn8bSVNz/e/UtK3JG2RLztd0n2SfijpeeDCSv8trHfwbxTWJ+VfircDtwETgXrgbklLI+LODla/EVgMDAf2BuZIWh4Rc0vY9c7AtsAuwDjgZkkzI+JaSR8EmiPiWyXWeAHQAOwODABmF9nfROB4YG1EbJT0JPC3wDPAp4CfSdozIpZIOgf4fEQcXrD+v+b17g7sANwFtADX58sPAWYAOwL9Szh+68XcQ7HebKak9Zt+gKsLlv0NMDQiLo6INyJiOXAdcEp7G5S0K3A4MCUiNkTEQuDHwKkl1vQmcHFEvBkRs4FXgPe10bajGk8GvhcRL0REM3BFkW1cERGrIuIvABFxU0Ssjoi3I+K/gCeAg9s41n7A/wW+EREvR8RTwA9aHevqiPjXiNi4aR/Wd7mHYr3ZxyPi7k0Tkk4HPp9P7gYM33RqJ9cPuLeDbQ4Hno+IlwvmrQQaS6xpXURsLJh+DRjYRtuOahwOrCpYVvi56DxJnwO+RtazId/3kDb2PwTYiuz4NllJ1rtqb5/WRzlQrK9aBayIiFFlrrca2F7SoIJQGQH8OUFNrYf+7qjGFrLTYI/n07u2t01Ju5H1cI4BfhcRb0laCKiN/a8l61HtVrCP1sfq4crtHT7lZX3V74GX8ovW75HUT9K+kv6mvZUiYhVwP3CppDpJHwDOAn6eoKZnya5VlFrjL4BvSNpO0i7ApA62P4AsANYA5DcE7Ntq//WStgLIb5P+BfBdSYPyQPoa8LPOHab1Vg4U65PyL8sTgDHACrLfxn9MdgG6IxPJThmtBm4FLoiIOQnKuh4YnV/zmVlCjRcDzfmyu4Gbgdfb2nhEPE52DeR3ZOGxH3BfQZN5ZDcbPCNpbT7vK8CrwHLgt8B/Ajd09kCtd5JfsGXWO0j6EnBKRBzR3bVY3+QeilmNkjRM0ockbSHpfcDXyXpMZt3CF+XNatdWwL8DI4H1ZM+DXN3eCmbV5FNeZmaWhE95mZlZEg4Usx5I0mck3VVCu3eG2e8JunvEZOteDhSrefrr+0o2/YSkVwum/7aCbb5r+PtWy4+U9Ha+/ZclLd000GMF+9ps8EmAiPh5RHykku2ZdRdflLeaFxFPUzB8iaQA9o+IZVXe9eqIqJckYALZQI8P5s97lEQe8t16EfdQrFeTtLWkyyU9LelZZW9QfE++bIikX+UPEj4v6d78Ftyfkg0xcnveA5nc3j4iMxN4gezBxOMl/VHSS5JWSbqwoJ5NvZGzJD1N9jDhPfni9fn+DlOrl11Jer+kOXmdz0o6v43jPVTS/fkxPSLpyIJlp0tanveoVkj6TDt/Zz+StDr/+ZGkrfNlbQ7BX2Q7j0k6oWC6v6S1ksa09/dptcuBYr3dZcBeZE+b70k2sOG382VfJ3vSfCiwE3A+WT6cCjwNnJC/qfGf29tBHkInAYOBRWRPln8unz4e+JKkj7da7QhgH+CjwNh83uB8f79rtf1BZE/C30E2IOSewLuGys+HX/lv4DvA9sB5wC2ShkoaQDYa8fiIGAR8EFjYxiF9EziU7O9sf7LRiL9VsLxwCP6zgKskbVdkO9OBzxZMHwe05CM0Wy/kQLFeKz8V9QXg3IjYNELw9/jr8O9vAsOA3fLh5O+N8u6j3zQS8Fqyd5OcGhFLI2J+RCzKh4h/lOz9Ka2fXr8wIl4tccj3jwHPRMQP8iHzX46IB4u0+ywwOyJm5/ueAywg+yIHeBvYV9J7IqIlIha3sb/PkA2x/1xErAEuYvMh60sdgv9nwHGS3ptPnwr8tITjtRrlQLHebCiwDfCw/vpOlDvy+QDfB5YBd+WngprK3P7qiBgcEdtHxJiImAEg6RBJv1H2lsMXgXN49xDx5Qz7vivwZAntdgM+pc3fAXM4MCwiXiV7t8k5QIuk/5a0dxvbGc67h6wfXjBd0hD8EbGabKywT0gaDIwnzSCa1kM5UKw3Wwv8BXh//sU/OCK2jYiBAPlv+l+PiN3JBmH8mqRj8nU788TvfwKzgF0jYluyV+uqVZto43Mxq4A9StjvKuCnBcc6OCIGRMRUgIi4MyLGkfXK/kQ2lH0xq8nCaZMR+bxKTCPrOX2KbMj8FMP8Ww/lQLFeKyLeJvvS/KGkHSG7ziDpo/nnj0naMz819hLwVv4D7x5KvhyDyF7CtUHSwcCnO2i/hux0VFv7+xWws6R/yC+YD5J0SJF2PwNOkPRRZUPd1+UX0esl7STpxPxayutkp6neKrINyE7RfSu/9jKE7JpTpc+6zAQOBL5Kdk3FejEHivV2U8hOaz0g6SWyi9ubzvePyqdfIRvS/eqImJ8vu5TsS3W9pPPK3OeXgYslvUz2ZfyL9hpHxGvAd4H78v0d2mr5y2Tvnz+B7F3wTwBHFdnOKrLbl88nC6lVwD+S/X++BdlNCKuB58mu6Xy5jZK+Q3bt5VGymwz+kM8rW36N6Bay8cZ+Wck2rHZ4LC8zqypJ3wb2iojPdtjYapofqjKzqpG0Pdmtxad21NZqn095mVlVSPoC2Wm3X0fEPR21t9rnU15mZpaEeyhmZpZEn76GMmTIkGhoaOjuMszMasrDDz+8NiKGtp7fpwOloaGBBQsWdHcZZmY1RdLKYvN9ysvMzJJwoJiZWRIOFDMzS6JPX0MxM0vhzTffpLm5mQ0bNnR3KUnV1dVRX19P//79S2rvQDEz66Tm5mYGDRpEQ0MD2VijtS8iWLduHc3NzYwcObKkdXzKy8yskzZs2MAOO+zQa8IEQBI77LBDWb0uB4qZWQK9KUw2KfeYHChmZpaEA8XMrIf74Ac/WHT+6aefzs0339zF1bTNgWJm1sPdf//93V1CSXyXl5lZDzdw4EBeeeUVIoKvfOUrzJs3j5EjR9LTRot3D8XMrEbceuutLF26lEWLFnHdddf1uJ6LA8XMrEbcc889TJw4kX79+jF8+HCOPvro7i5pMw4UM7Ma0pNvT3agmJnViLFjxzJjxgzeeustWlpa+M1vftPdJW3GF+XNzGrESSedxLx589hvv/3Ya6+9OOKII7q7pM04UMzMerhXXnkFyE53XXnlld1cTdt8ysvMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM+sFzjzzTHbccUf23Xffd+Y9//zzjBs3jlGjRjFu3DheeOEFANatW8dRRx3FwIEDmTRpUrIaHChmZr3A6aefzh133LHZvKlTp3LMMcfwxBNPcMwxxzB16lQA6urquOSSS7j88suT1tCjAkXSsZKWSlomqanIckm6Il/+qKQDWy3vJ+mPkn7VdVWbmXW/sWPHsv32228277bbbuO0004D4LTTTmPmzJkADBgwgMMPP5y6urqkNfSYJ+Ul9QOuAsYBzcBDkmZFxOMFzcYDo/KfQ4B/y//c5KvAEuC9XVK0mVkrF92+mMdXv5R0m6OHv5cLTnh/2es9++yzDBs2DIBhw4bx3HPPJa2rtZ7UQzkYWBYRyyPiDWAGMKFVmwnA9Mg8AAyWNAxAUj1wPPDjrizazMwyPaaHAuwCrCqYbmbz3kdbbXYBWoAfAZOBQe3tRNLZwNkAI0aM6FTBZmatVdKTqJaddtqJlpYWhg0bRktLCzvuuGNV99eTeijFBvlv/X7Lom0kfQx4LiIe7mgnEXFtRDRGROPQoUMrqdPMrCaceOKJTJs2DYBp06YxYULrkz5p9aQeSjOwa8F0PbC6xDafBE6UdBxQB7xX0s8i4rNVrNfMrMeYOHEi8+fPZ+3atdTX13PRRRfR1NTEySefzPXXX8+IESO46aab3mnf0NDASy+9xBtvvMHMmTO56667GD16dKdq6EmB8hAwStJI4M/AKcCnW7WZBUySNIPsdNiLEdECfCP/QdKRwHkOEzPrS2688cai8+fOnVt0/lNPPZW8hh4TKBGxUdIk4E6gH3BDRCyWdE6+/BpgNnAcsAx4DTiju+o1M7PN9ZhAAYiI2WShUTjvmoLPAfxdB9uYD8yvQnlmZtaOnnRR3szMapgDxczMknCgmJlZEg4UMzNLwoFiZtYLlDN8/Zw5czjooIPYb7/9OOigg5g3b16SGhwoZma9QDnD1w8ZMoTbb7+dRYsWMW3aNE499dQkNThQzMx6gXKGrz/ggAMYPnw4AO9///vZsGEDr7/+eqdr6FHPoZiZ1bxfN8Ezi9Juc+f9YPzUslcrZfj6W265hQMOOICtt96602U6UMzM+qjFixczZcoU7rrrriTbc6CYmaVUQU+iWtobvr65uZmTTjqJ6dOns8ceeyTZn6+hmJn1Um0NX79+/XqOP/54Lr30Uj70oQ8l258DxcysF5g4cSKHHXYYS5cupb6+nuuvv56mpibmzJnDqFGjmDNnDk1NTQBceeWVLFu2jEsuuYQxY8YwZsyYJK8HVjbeYt/U2NgYCxYs6O4yzKzGLVmyhH322ae7y6iKYscm6eGIaGzd1j0UMzNLwoFiZmZJOFDMzCwJB4qZmSXhQDEzsyQcKGZmloQDxcysFyhn+Prf//737zx/sv/++3PrrbcmqcGBYmbWC5QzfP2+++7LggULWLhwIXfccQdf/OIX2bhxY6drcKCYmfUC5Qxfv80227DlltlQjhs2bEBSkho8OKSZWUKX/f4y/vT8n5Juc+/t92bKwVPKXq+94esffPBBzjzzTFauXMlPf/rTdwKmM9xDMTPrgw455BAWL17MQw89xKWXXsqGDRs6vU33UMzMEqqkJ1Et7Q1fv8k+++zDgAEDeOyxx2hsfNfwXGVxD8XMrJdqa/j6FStWvHMRfuXKlSxdupSGhoZO7889FDOzXmDixInMnz+ftWvXUl9fz0UXXURTUxMnn3wy119/PSNGjOCmm24C4Le//S1Tp06lf//+bLHFFlx99dUMGTKk0zV4+HoPX29mneTh6zM+5WVmZkk4UMzMLIkeFSiSjpW0VNIySU1FlkvSFfnyRyUdmM/fVdJvJC2RtFjSV7u+ejOzvq3HBIqkfsBVwHhgNDBR0uhWzcYDo/Kfs4F/y+dvBL4eEfsAhwJ/V2RdMzOroh4TKMDBwLKIWB4RbwAzgAmt2kwApkfmAWCwpGER0RIRfwCIiJeBJcAuXVm8mVlf15MCZRdgVcF0M+8OhQ7bSGoADgAeTF+imZm1pScFSrHRyVrf09xuG0kDgVuAf4iIl4ruRDpb0gJJC9asWVNxsWZmPUk5w9dv8vTTTzNw4EAuv/zyJDV0GCiSRpT4895O1tIM7FowXQ+sLrWNpP5kYfLziPhlWzuJiGsjojEiGocOHdrJks3MeoZyhq/f5Nxzz2X8+PHJaijlSflpZL2A9sY3DuAnwPRO1PIQMErSSODPwCnAp1u1mQVMkjQDOAR4MSJalI29fD2wJCL+XydqMDOrSWPHjuWpp57abN5tt93G/PnzgWz4+iOPPJLLLrsMgJkzZ7L77rszYMCAZDV0GCgRcVTreZJ2johnklWR7WejpEnAnUA/4IaIWCzpnHz5NcBs4DhgGfAacEa++oeAU4FFkhbm886PiNkpazQz68gz3/sery9JO3z91vvszc7nn1/2em0NX//qq69y2WWXMWfOnGSnu6Dysbw+B/xzsipyeQDMbjXvmoLPAfxdkfV+S/s9KDMzy11wwQWce+65DBw4MOl2Kw2UCZJeA+ZExNKUBZmZ1bJKehLV0tbw9Q8++CA333wzkydPZv369WyxxRbU1dUxadKkTu2v0kD5P2S35p4kac+I+HynqjAzs+Q2DV/f1NS02fD199577zttLrzwQgYOHNjpMIEKAyUingXuyH/MzKyblTN8fbVUFCiSrgIGRMTpkj4SEXclrsvMzMpw4403Fp0/d+7cdte78MILk9VQ6YONbwDL889HJ6rFzMxqWKWB8hqwbf4w4YiE9ZiZWY2q9KL888BfyEYHvi9dOWZmtSkiyJ6x7j3KfaNvWT0USYMl/QfwiXzWdOBdr4E0M+tL6urqWLduXdlfwD1ZRLBu3Trq6upKXqesHkpErJc0FWgA1gIfANocN8vMrC+or6+nubmZ3jbgbF1dHfX19SW3r+SU11nAioi4E3i4gvXNzHqV/v37M3LkyO4uo9tVEigvAOdIeh/wCLAwIv6YtiwzM6s1ZQdKRFwqaS7wP8AYYCzgQDEz6+PKDhRJF5ONBryQrHcyP3FNZmZWgyrpoXxb0k5kY3l9QtIeEfGF9KWZmVktqfQ5lC8C/x4RHsvLzMyAygPlBuBLkgaQvXJ3YbqSzMysFlU69Mrfk4XRlsAV6coxM7NaVWmgPAnUAbdFxNiE9ZiZWY2qNFAWA/OAsyQ9lLAeMzOrUZVeQ9kLWANcS/ago5mZ9XGV9lD2JnuY8Tzg7HTlmJlZrao0UAYDU4DJwIZk1ZiZWc2q9JTXxcDeEbFU0tspCzIzs9pUUg9FUj9JLZI+DxARzRFxd/65qZoFmplZbSgpUCLiLeAxYI/qlmNmZrWqnFNe2wCTJY0DVufzIiImpC/LzMxqTTmBclj+54H5D0Dved+lmZl1SjmB4teRmZlZm0oOlIhYWc1CzMystlX6HIqZmdlmHChmZpZE2YEi6YRqFJJv+1hJSyUtk/Su51uUuSJf/qikA0td18zMqquSHsp3k1dB9vAkcBUwHhgNTJQ0ulWz8cCo/Ods4N/KWNfMzKqokqFXlLyKzMHAsohYDiBpBjABeLygzQRgekQE8ICkwZKGAQ0lrJvMf513Ilv/aUU1Nm1m1iW2/ewZHHXK15Jus5JAqdazJ7sAqwqmm4FDSmizS4nrAiDpbPIRkkeMGFFRofHierZ9YWNF65qZ9QSvv5z+zSOVDg5ZDcV6Pq3Dq602paybzYy4luw9LjQ2NlYUjqdcd08lq5mZ9Wo9KVCagV0Lpuv56xAvHbXZqoR1zcysiiq5KP9s8ioyDwGjJI2UtBVwCjCrVZtZwOfyu70OBV6MiJYS1zUzsyoqu4cSEeOqUUhEbJQ0CbgT6AfcEBGLJZ2TL78GmA0cBywDXgPOaG/datRpZmbFKbthqm9qbGyMBQsWdHcZZmY1RdLDEdHYer6flDczsyQqChRJXyv4/L505ZiZWa0q6xqKpMHAD4G9JW0AHgXOIr+WYWZmfVdZgRIR64EzJB0PPAN8BPhlFeoyM7MaU+k1lCPIbh8+FKjKXV9mZlZbKg2UwcAUYDKwIVk1ZmZWsyp9Uv5iYO+IWCrp7ZQFmZlZbaooUCKimWwYFCLC7x4xM7OKbxu+StJP8s8fSVqRmZnVpEqvobwBLM8/H52oFjMzq2GVBsprwLaS+gOVvVTEzMx6lUovyj8P/IXstbv3pSvHzMxqVVk9lPyVu/8BfCKfNR141wBhZmbW95T9pLykqWTvcF8LfAA/KW9mZlR2yussYEVE3Ak8nLgeMzOrUZUEygvAOfkow48ACyPij2nLMjOzWlPJGxsvlTQX+B9gDDAWcKCYmfVxZQeKpIvJXrO7kKx3Mj9xTWZmVoPKfg4lIr4NvJ6v+wlJ1yWvyszMak6lDzbeAOwD7ABcna4cMzOrVZUGyt+TnS7bEviXdOWYmVmtqjRQngTqgNsiYmzCeszMrEZVGiiLgXnAWZIeSliPmZnVqErH8tqD7HmUa/M/zcysj6s0UFZFxDxJw4DnUhZkZma1qdJTXsdKqgeuAX6YsB4zM6tRlQbKYGAKMJnsmRQzM+vjKj3ldTGwd0QslfRWyoLMzKw2ldRDkdRPUoukzwNERHNE3J1/bqpmgWZmVhtKCpSIeAt4jOzuLjMzs3cp5xrKNsBkSQskzcp/bktRhKTtJc2R9ET+53ZttDtW0lJJyyQ1Fcz/vqQ/SXpU0q2SBqeoy8zMSldOoBwGCDgQ+FjBTwpNwNyIGAXMzac3I6kf2TvsxwOjgYmSRueL5wD7RsQHyIbV/0aiuszMrETlXJQfWbUqYAJwZP55GjCf7C6yQgcDyyJiOYCkGfl6j0fEXQXtHgA+WcVazcysiA4DRdKI/GN0sHx9RLxUYR07RUQLQES0SNqxSJtdgFUF083AIUXanQn8V4V1mJlZhUrpoUwjCxO10yaAnwDT22og6W5g5yKLvllCDbSx/81CTtI3gY3Az9up42zgbIARI0a01czMzMrUYaBExFEpdhQRH25rmaRnJQ3LeydtDefSDOxaMF0PrC7Yxmlk13SOiYiivam8jmvJxiCjsbGxzXZmZlaeSp+UT20WcFr++TSg2N1jDwGjJI2UtBVwSr4eko4lu+ZyYkS81gX1mplZKz0lUKYC4yQ9AYzLp5E0XNJsgIjYCEwC7gSWAL+IiMX5+lcCg4A5khZKuqarD8DMrK+rdOiVpCJiHXBMkfmrgeMKpmcDs4u027OqBZqZWYd6Sg/FzMxqnAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJdEjAkXS9pLmSHoi/3O7NtodK2mppGWSmoosP09SSBpS/arNzKxQjwgUoAmYGxGjgLn59GYk9QOuAsYDo4GJkkYXLN8VGAc83SUVm5nZZnpKoEwApuWfpwEfL9LmYGBZRCyPiDeAGfl6m/wQmAxEFes0M7M29JRA2SkiWgDyP3cs0mYXYFXBdHM+D0knAn+OiEc62pGksyUtkLRgzZo1na/czMwA2LKrdiTpbmDnIou+WeomiswLSdvk2/hIKRuJiGuBawEaGxvdmzEzS6TLAiUiPtzWMknPShoWES2ShgHPFWnWDOxaMF0PrAb2AEYCj0jaNP8Pkg6OiGeSHYCZmbWrp5zymgWcln8+DbitSJuHgFGSRkraCjgFmBURiyJix4hoiIgGsuA50GFiZta1ekqgTAXGSXqC7E6tqQCShkuaDRARG4FJwJ3AEuAXEbG4m+o1M7NWuuyUV3siYh1wTJH5q4HjCqZnA7M72FZD6vrMzKxjPaWHYmZmNc6BYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRKKiO6uodtIWgOsrHD1IcDahOXUAh9z3+Bj7hs6c8y7RcTQ1jP7dKB0hqQFEdHY3XV0JR9z3+Bj7huqccw+5WVmZkk4UMzMLAkHSuWu7e4CuoGPuW/wMfcNyY/Z11DMzCwJ91DMzCwJB4qZmSXhQOmApGMlLZW0TFJTkeWSdEW+/FFJB3ZHnSmVcMyfyY/1UUn3S9q/O+pMqaNjLmj3N5LekvTJrqwvtVKOV9KRkhZKWizp/3d1jamV8N/1tpJul/RIfsxndEedKUm6QdJzkh5rY3na76+I8E8bP0A/4Elgd2Ar4BFgdKs2xwG/BgQcCjzY3XV3wTF/ENgu/zy+LxxzQbt5wGzgk91dd5X/jQcDjwMj8ukdu7vuLjjm84HL8s9DgeeBrbq79k4e91jgQOCxNpYn/f5yD6V9BwPLImJ5RLwBzAAmtGozAZgemQeAwZKGdXWhCXV4zBFxf0S8kE8+ANR3cY2plfLvDPAV4Bbgua4srgpKOd5PA7+MiKcBIqIvHHMAgyQJGEgWKBu7tsy0IuIesuNoS9LvLwdK+3YBVhVMN+fzym1TS8o9nrPIfsOpZR0es6RdgJOAa7qwrmop5d94L2A7SfMlPSzpc11WXXWUcsxXAvsAq4FFwFcj4u2uKa/bJP3+2rLT5fRuKjKv9X3WpbSpJSUfj6SjyALl8KpWVH2lHPOPgCkR8Vb2C2xNK+V4twQOAo4B3gP8TtIDEfE/1S6uSko55o8CC4GjgT2AOZLujYiXqlxbd0r6/eVAaV8zsGvBdD3Zby/ltqklJR2PpA8APwbGR8S6LqqtWko55kZgRh4mQ4DjJG2MiJldUmFapf53vTYiXgVelXQPsD9Qq4FSyjGfAUyN7OLCMkkrgL2B33dNid0i6feXT3m17yFglKSRkrYCTgFmtWozC/hcfrfEocCLEdHS1YUm1OExSxoB/BI4tYZ/Yy3U4TFHxMiIaIiIBuBm4Ms1GiZQ2n/XtwF/K2lLSdsAhwBLurjOlEo55qfJemRI2gl4H7C8S6vsekm/v9xDaUdEbJQ0CbiT7C6RGyJisaRz8uXXkN3xcxywDHiN7LecmlXiMX8b2AG4Ov+NfWPU8EitJR5zr1HK8UbEEkl3AI8CbwM/joiit57WghL/jS8BfiJpEdmpoCkRUdND2ku6ETgSGCKpGbgA6A/V+f7y0CtmZpaET3mZmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMUtE0mBJXy6YHi7p5irt6+OSvt1Bm8slHV2N/ZsV49uGzRKR1AD8KiL27YJ93Q+c2N5zEpJ2A66LiI9Uux4zcA/FLKWpwB75O0S+L6lh03soJJ0uaWb+vo0VkiZJ+pqkP0p6QNL2ebs9JN2RD8h4r6S9W+9E0l7A6xGxVtKgfHv982XvlfSUpP4RsRLYQdLOXfh3YH2YA8UsnSbgyYgYExH/WGT5vmTDwh8MfBd4LSIOAH4HbBrN91rgKxFxEHAecHWR7XwI+ANARLwMzAeOz5edAtwSEW/m03/I25tVnYdeMes6v8kD4GVJLwK35/MXAR+QNJDs5WU3FYxovHWR7QwD1hRM/xiYDMwkGzrjCwXLngOGpzoAs/Y4UMy6zusFn98umH6b7P/FLYD1ETGmg+38Bdh200RE3JefXjsC6NdqzK26vL1Z1fmUl1k6LwODKl05f+/GCkmfgnfe971/kaZLgD1bzZsO3Aj8R6v5ewE1O6ij1RYHilki+Xth7pP0mKTvV7iZzwBnSXoEWEzxVxHfAxygzd/09XNgO7JQASC/UL8nsKDCWszK4tuGzWqQpH8Bbo+Iu/PpTwITIuLUgjYnAQdGxD91U5nWx/gaillt+h7ZS6+Q9K/AeLL3WhTaEvhBF9dlfZh7KGZmloSvoZiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl8b+j0m08j6dtgAAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAANQklEQVR4nO3cX4il9X3H8fenuxEak0aJk5DurmRb1pi90KITI6VpTUObXXuxBLxQQ6QSWKQx5FIpNLnwprkohKBmWWSR3GQvGkk2ZRMplMSCNd1Z8N8qynSlOl3BNYYUDFRWv704p51hnHWenXNmZp3v+wUD85znNzPf+TH73mfPznlSVUiStr7f2ewBJEkbw+BLUhMGX5KaMPiS1ITBl6QmDL4kNbFq8JMcSfJakmfPcz5JvptkPsnTSa6b/piSpEkNucJ/GNj3Huf3A3vGbweB700+liRp2lYNflU9BrzxHksOAN+vkSeAy5J8YloDSpKmY/sUPscO4JUlxwvjx15dvjDJQUb/CuDSSy+9/uqrr57Cl5ekPk6ePPl6Vc2s5WOnEfys8NiK92uoqsPAYYDZ2dmam5ubwpeXpD6S/OdaP3Yav6WzAOxacrwTODOFzytJmqJpBP8YcMf4t3VuBH5TVe96OkeStLlWfUonyQ+Am4ArkiwA3wI+AFBVh4DjwM3APPBb4M71GlaStHarBr+qblvlfAFfm9pEkqR14SttJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJamJQ8JPsS/JCkvkk965w/iNJfpLkqSSnktw5/VElSZNYNfhJtgEPAPuBvcBtSfYuW/Y14Lmquha4CfiHJJdMeVZJ0gSGXOHfAMxX1emqegs4ChxYtqaADycJ8CHgDeDcVCeVJE1kSPB3AK8sOV4YP7bU/cCngTPAM8A3quqd5Z8oycEkc0nmzp49u8aRJUlrMST4WeGxWnb8ReBJ4PeBPwLuT/J77/qgqsNVNVtVszMzMxc4qiRpEkOCvwDsWnK8k9GV/FJ3Ao/UyDzwEnD1dEaUJE3DkOCfAPYk2T3+j9hbgWPL1rwMfAEgyceBTwGnpzmoJGky21dbUFXnktwNPApsA45U1akkd43PHwLuAx5O8gyjp4DuqarX13FuSdIFWjX4AFV1HDi+7LFDS94/A/zldEeTJE2Tr7SVpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDUxKPhJ9iV5Icl8knvPs+amJE8mOZXkF9MdU5I0qe2rLUiyDXgA+AtgATiR5FhVPbdkzWXAg8C+qno5ycfWaV5J0hoNucK/AZivqtNV9RZwFDiwbM3twCNV9TJAVb023TElSZMaEvwdwCtLjhfGjy11FXB5kp8nOZnkjpU+UZKDSeaSzJ09e3ZtE0uS1mRI8LPCY7XseDtwPfBXwBeBv0ty1bs+qOpwVc1W1ezMzMwFDytJWrtVn8NndEW/a8nxTuDMCmter6o3gTeTPAZcC7w4lSklSRMbcoV/AtiTZHeSS4BbgWPL1vwY+FyS7Uk+CHwWeH66o0qSJrHqFX5VnUtyN/AosA04UlWnktw1Pn+oqp5P8jPgaeAd4KGqenY9B5ckXZhULX86fmPMzs7W3NzcpnxtSXq/SnKyqmbX8rG+0laSmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmBgU/yb4kLySZT3Lve6z7TJK3k9wyvRElSdOwavCTbAMeAPYDe4Hbkuw9z7pvA49Oe0hJ0uSGXOHfAMxX1emqegs4ChxYYd3XgR8Cr01xPknSlAwJ/g7glSXHC+PH/l+SHcCXgEPv9YmSHEwyl2Tu7NmzFzqrJGkCQ4KfFR6rZcffAe6pqrff6xNV1eGqmq2q2ZmZmYEjSpKmYfuANQvAriXHO4Ezy9bMAkeTAFwB3JzkXFX9aBpDSpImNyT4J4A9SXYD/wXcCty+dEFV7f6/95M8DPyTsZeki8uqwa+qc0nuZvTbN9uAI1V1Ksld4/Pv+by9JOniMOQKn6o6Dhxf9tiKoa+qv558LEnStPlKW0lqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSE4OCn2RfkheSzCe5d4XzX07y9Pjt8STXTn9USdIkVg1+km3AA8B+YC9wW5K9y5a9BPxZVV0D3AccnvagkqTJDLnCvwGYr6rTVfUWcBQ4sHRBVT1eVb8eHz4B7JzumJKkSQ0J/g7glSXHC+PHzuerwE9XOpHkYJK5JHNnz54dPqUkaWJDgp8VHqsVFyafZxT8e1Y6X1WHq2q2qmZnZmaGTylJmtj2AWsWgF1LjncCZ5YvSnIN8BCwv6p+NZ3xJEnTMuQK/wSwJ8nuJJcAtwLHli5IciXwCPCVqnpx+mNKkia16hV+VZ1LcjfwKLANOFJVp5LcNT5/CPgm8FHgwSQA56pqdv3GliRdqFSt+HT8upudna25ublN+dqS9H6V5ORaL6h9pa0kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNGHxJasLgS1ITBl+SmjD4ktSEwZekJgy+JDVh8CWpCYMvSU0YfElqwuBLUhMGX5KaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrC4EtSEwZfkpow+JLUhMGXpCYMviQ1YfAlqQmDL0lNDAp+kn1JXkgyn+TeFc4nyXfH559Oct30R5UkTWLV4CfZBjwA7Af2Arcl2bts2X5gz/jtIPC9Kc8pSZrQkCv8G4D5qjpdVW8BR4EDy9YcAL5fI08AlyX5xJRnlSRNYPuANTuAV5YcLwCfHbBmB/Dq0kVJDjL6FwDA/yR59oKm3bquAF7f7CEuEu7FIvdikXux6FNr/cAhwc8Kj9Ua1lBVh4HDAEnmqmp2wNff8tyLRe7FIvdikXuxKMncWj92yFM6C8CuJcc7gTNrWCNJ2kRDgn8C2JNkd5JLgFuBY8vWHAPuGP+2zo3Ab6rq1eWfSJK0eVZ9SqeqziW5G3gU2AYcqapTSe4anz8EHAduBuaB3wJ3Dvjah9c89dbjXixyLxa5F4vci0Vr3otUveupdknSFuQrbSWpCYMvSU2se/C9LcOiAXvx5fEePJ3k8STXbsacG2G1vViy7jNJ3k5yy0bOt5GG7EWSm5I8meRUkl9s9IwbZcCfkY8k+UmSp8Z7MeT/C993khxJ8tr5Xqu05m5W1bq9MfpP3v8A/gC4BHgK2Ltszc3ATxn9Lv+NwC/Xc6bNehu4F38MXD5+f3/nvViy7l8Y/VLALZs99yb+XFwGPAdcOT7+2GbPvYl78bfAt8fvzwBvAJds9uzrsBd/ClwHPHue82vq5npf4XtbhkWr7kVVPV5Vvx4fPsHo9Qxb0ZCfC4CvAz8EXtvI4TbYkL24HXikql4GqKqtuh9D9qKADycJ8CFGwT+3sWOuv6p6jNH3dj5r6uZ6B/98t1y40DVbwYV+n19l9Df4VrTqXiTZAXwJOLSBc22GIT8XVwGXJ/l5kpNJ7tiw6TbWkL24H/g0oxd2PgN8o6re2ZjxLipr6uaQWytMYmq3ZdgCBn+fST7PKPh/sq4TbZ4he/Ed4J6qent0MbdlDdmL7cD1wBeA3wX+LckTVfXieg+3wYbsxReBJ4E/B/4Q+Ock/1pV/73Os11s1tTN9Q6+t2VYNOj7THIN8BCwv6p+tUGzbbQhezELHB3H/grg5iTnqupHGzLhxhn6Z+T1qnoTeDPJY8C1wFYL/pC9uBP4+xo9kT2f5CXgauDfN2bEi8aaurneT+l4W4ZFq+5FkiuBR4CvbMGrt6VW3Yuq2l1Vn6yqTwL/CPzNFow9DPsz8mPgc0m2J/kgo7vVPr/Bc26EIXvxMqN/6ZDk44zuHHl6Q6e8OKypm+t6hV/rd1uG952Be/FN4KPAg+Mr23O1Be8QOHAvWhiyF1X1fJKfAU8D7wAPVdWWu7X4wJ+L+4CHkzzD6GmNe6pqy902OckPgJuAK5IsAN8CPgCTddNbK0hSE77SVpKaMPiS1ITBl6QmDL4kNWHwJakJgy9JTRh8SWrifwHXe3WluIZOawAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -148,22 +164,9 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAElCAYAAADDUxRwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAiXUlEQVR4nO3dfbxVZZ338c+3A4QKigrIwwFBQwFREQjUzNQGBswyFBtRS80iTZsa9VaqmUznnrSa8iEthsyn7JYcK0XDJ0BvC2UUBVQkEhHjCAiixJMEHH73H2vhvT3tc84+e6/zsNnf9+t1Xuy91rWu9Vug+3uutde6liICMzOzUn2otQswM7PdgwPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDHLQ9J3Jd2dvu4raZOkqiL6mSLp37Kv0KztcaDYbknSckn/UGfZeZL+2NS+IuIvEdEpImqL2PbCiPj3QtpKukPS/27qPrJS7N+P2S4OFLMKIKnd7rAPa9scKFaxJPWS9BtJayW9Lumf62nXT1Ls+sBMt5su6R1JSyV9uYF9vD/qkHSCpBpJl0laI2mVpPPTdZOAs4Er0tNrDzZWo6Q9JN0p6V1JiyVdIakmZ/1ySVdKehHYLKmdpMmSXpO0UdIrksanbQcBU4Bj0v2vT5fvI+mudP9vSPpXSR9K150naY6k6yW9A3y32H8L2z34NwqrSOmH4oPAA8BEoBqYKWlJRDzayOb3AIuAXsBA4HFJyyJiVgG77gHsA/QGRgP3Sbo/IqZKOhaoiYh/LbDGq4B+wEHAXsCMPPubCHwKeDsidkh6Dfg4sBo4A7hb0kciYrGkC4EvRcRxOdv/JK33IGB/4DFgFfCLdP0oYBrQHWhfwPHbbswjFNud3S9p/a4f4Kc56z4KdIuIayJiW0QsA34OnNlQh5L6AMcBV0bE1ohYANwKfL7AmrYD10TE9oiYAWwCDq2nbWM1fg74XkS8GxE1wE15+rgpIlZExHsAEfHfEbEyInZGxK+BV4GR9RxrFfBPwDcjYmNELAd+VOdYV0bETyJix659WOXyCMV2Z5+NiJm73kg6D/hS+vZAoNeuUzupKuAPjfTZC3gnIjbmLHsDGFFgTesiYkfO+y1Ap3raNlZjL2BFzrrc13mXSfoCcCnJyIZ0313r2X9XoAPJ8e3yBsnoqqF9WoVyoFilWgG8HhEDmrjdSmA/SZ1zQqUv8GYGNdWd+ruxGleRnAZ7JX3fp6E+JR1IMsL5JPBMRNRKWgConv2/TTKiOjBnH3WP1dOV2/t8yssq1bPAhvRL6z0kVUkaIumjDW0UESuAp4FrJXWUdARwAfCrDGp6i+S7ikJrvBf4pqR9JfUGLmmk/71IAmAtQHpBwJA6+6+W1AEgvUz6XuA/JHVOA+lS4O7SDtN2Vw4Uq0jph+WngaHA6yS/jd9K8gV0YyaSnDJaCfwOuCoiHs+grF8Ag9PvfO4voMZrgJp03UzgPuBv9XUeEa+QfAfyDEl4HA7MyWkym+Rig9WS3k6XfQ3YDCwD/gj8H+C2Ug/Udk/yA7bMdg+SLgLOjIhPtHYtVpk8QjErU5J6SvqYpA9JOhS4jGTEZNYq/KW8WfnqAPwX0B9YT3I/yE8b2sCsOfmUl5mZZcKnvMzMLBMOFLMmyDeL8e6i7pxlZk3lQDGrI/1Q3ZxOkvimpB+riGehZFDDR1pyn2alcqCY5XdkRHQiuav8LKDeGYXNLOFAMWtARPyJZO6sIXXXSRop6Zn0RsRVkm7edZd5uj4kXSjp1XSK+VskKWf9F9Np59+V9Gh6JzqSnkqbLExHSf8kqaukh9J9vSPpD7umkc9T17GSnpP01/TPY3PWPSnp39Np5zdKekzS383lJekMSc/XWXaZpPub9jdolcSBYtYASYNJpnufn2d1LfAvJJMoHkMymvlqnTankMwafCTJ7MD/mPb7WeBbwGlAN5LQugcgIo5Ptz0yfVLkr0nuMalJ2x6Qbvt3l2hK2g/4PcnMw/sDPwZ+L2n/nGZnAeeTTDnfAbg8z7FNB/qnz0nZ5Rzgl3namgEOFLP6vCDpXZLnkdwK3F63QUQ8HxFz06nbl5PcE1L3LvXrImJ9RPwFeIJkGhWArwDXRsTidPbh7wFDd41S8tgO9AQOTKe+/0Pkv+b/U8CrEfHLtK57gD+RTOGyy+0R8ed0uvl7c2rKPba/Ab8mCREkHUYy3cxD9dRn5kAxq8ewiNg3Ig6OiH+NiJ11G0g6JD0NtVrSBpJQqHv6aHXO69yp6g8Ebsx5Vss7JLP+9ia/HwJLgcckLZM0uZ52vfjgdPPw91PO11dTXXcCZ6Wn6T4P3JsGjVleDhSz4v2M5Lf/ARGxN8lpKDW8yftWAF+JiC45P3tExNP5GqcPuLosIg4iGW1cKumTeZquJAmrXEVNrx8Rc4FtJKf8zsKnu6wRDhSz4nUGNgCbJA0ELmrCtlNIpp4/DN5/dvsZOes/MJW9pFMkfSQdLWwg+f6mNk+/M4BDJJ2l5Bny/wQMpvhTVXcBNwM7IuKPRfZhFcKBYla8y0l+c99I8uCqXxe6YUT8Dvg+MC09XfYyMC6nyXeBO9NTYp8DBpBMUb+JZPr5n0bEk3n6XUdyIcBlwDrgCuCUiHi7btsC/ZLkCjePTqxRnsvLzOolaQ9gDcl3Sq+2dj3WtnmEYmYNuQh4zmFihfCcPWaWl6TlJBcZfLZ1K7Fy4VNeZmaWCZ/yMjOzTFT0Ka+uXbtGv379WrsMM7Oy8vzzz78dEd3qLq/oQOnXrx/z5s1r7TLMzMqKpLqzMQA+5WVmZhlxoJiZWSYcKGZmlomK/g7FzKw1bN++nZqaGrZu3drapTSoY8eOVFdX0759+4LaO1DMzFpYTU0NnTt3pl+/fuQ8xLNNiQjWrVtHTU0N/fv3L2gbn/IyM2thW7duZf/992+zYQIgif33379JoygHiplZK2jLYbJLU2t0oJiZWSYcKGZmZerYY4/Nu/y8887jvvvua+FqHChmZmXr6afzPjG61fgqLzOzMtWpUyc2bdpERPC1r32N2bNn079/f1prFnmPUMzMytzvfvc7lixZwksvvcTPf/7zVhu5OFDMzMrcU089xcSJE6mqqqJXr16cdNJJrVKHA8XMbDfQFi5DdqCYmZW5448/nmnTplFbW8uqVat44oknWqUOfylvZlbmxo8fz+zZszn88MM55JBD+MQnPtEqdThQzMzK1KZNm4DkdNfNN9/cytX4lJeZmWXEgWJmZplwoJiZWSYcKGZmlgkHipmZZcKBYmZmmXCgmJlVqC9+8Yt0796dIUOGZNKfA8XMrEKdd955PPLII5n116YCRdJYSUskLZU0Oc96SbopXf+ipGF11ldJmi/poZar2sysPB1//PHst99+mfXXZu6Ul1QF3AKMBmqA5yRNj4hXcpqNAwakP6OAn6V/7vJ1YDGwd4sUbWZWoqsfXMQrKzdk2ufgXntz1acPy7TPQrSlEcpIYGlELIuIbcA04NQ6bU4F7orEXKCLpJ4AkqqBTwG3tmTRZmaWaDMjFKA3sCLnfQ0fHH3U16Y3sAq4AbgC6NzQTiRNAiYB9O3bt6SCzcxK1RojiebSlkYo+Sbzr/scy7xtJJ0CrImI5xvbSURMjYgRETGiW7duxdRpZmZ5tKVAqQH65LyvBlYW2OZjwGckLSc5VXaSpLubr1Qzs/I3ceJEjjnmGJYsWUJ1dTW/+MUvSuqvLZ3yeg4YIKk/8CZwJnBWnTbTgUskTSM5HfbXiFgFfDP9QdIJwOURcU4L1W1mVpbuueeeTPtrM4ESETskXQI8ClQBt0XEIkkXpuunADOAk4GlwBbg/Naq18zMPqjNBApARMwgCY3cZVNyXgdwcSN9PAk82QzlmZlZA9rSdyhmZlbGHChmZpYJB4qZmWXCgWJmZplwoJiZVaAVK1Zw4oknMmjQIA477DBuvPHGkvtsU1d5mZlZy2jXrh0/+tGPGDZsGBs3bmT48OGMHj2awYMHF92nRyhmZhWoZ8+eDBuWPAGkc+fODBo0iDfffLOkPj1CMTNrTQ9PhtUvZdtnj8Nh3HUFN1++fDnz589n1Ki68/E2jUcoZmYVbNOmTZx++unccMMN7L13aY+S8gjFzKw1NWEkkbXt27dz+umnc/bZZ3PaaaeV3J9HKGZmFSgiuOCCCxg0aBCXXnppJn06UMzMKtCcOXP45S9/yezZsxk6dChDhw5lxowZjW/YAJ/yMjOrQMcddxzJfLvZ8QjFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMKtDWrVsZOXIkRx55JIcddhhXXXVVyX36PhQzswr04Q9/mNmzZ9OpUye2b9/Occcdx7hx4zj66KOL7tMjFDOzCiSJTp06AcmcXtu3b0dSSX16hGJm1oq+/+z3+dM7f8q0z4H7DeTKkVc22q62tpbhw4ezdOlSLr74Yk9fb2ZmxamqqmLBggXU1NTw7LPP8vLLL5fUn0coZmatqJCRRHPr0qULJ5xwAo888ghDhgwpuh+PUMzMKtDatWtZv349AO+99x4zZ85k4MCBJfXpEYqZWQVatWoV5557LrW1tezcuZPPfe5znHLKKSX16UAxM6tARxxxBPPnz8+0T5/yMjOzTDhQzMwsE20qUCSNlbRE0lJJk/Osl6Sb0vUvShqWLu8j6QlJiyUtkvT1lq/ezKyytZlAkVQF3AKMAwYDEyUNrtNsHDAg/ZkE/CxdvgO4LCIGAUcDF+fZ1szMmlGbCRRgJLA0IpZFxDZgGnBqnTanAndFYi7QRVLPiFgVES8ARMRGYDHQuyWLNzOrdG0pUHoDK3Le1/D3odBoG0n9gKOA/8m+RDMzq09bCpR8s5JFU9pI6gT8BvhGRGzIuxNpkqR5kuatXbu26GLNzHYHtbW1HHXUUSXfgwIF3IciqW+Bfa2v70O8QDVAn5z31cDKQttIak8SJr+KiN/Wt5OImApMBRgxYkTdwDIzqyg33ngjgwYNYsOGUj6+E4Xc2HgnySigoXmNA7gDuKuEWp4DBkjqD7wJnAmcVafNdOASSdOAUcBfI2KVkjmXfwEsjogfl1CDmVnFqKmp4fe//z3f/va3+fGPS//obDRQIuLEussk9YiI1SXv/YP72SHpEuBRoAq4LSIWSbowXT8FmAGcDCwFtgDnp5t/DPg88JKkBemyb0XEjCxrNDPL2urvfY+/Lc52+voPDxpIj299q9F23/jGN/jBD37Axo0bM9lvsVOvfAH4QSYV5EgDYEadZVNyXgdwcZ7t/kjDIygzM8vx0EMP0b17d4YPH86TTz6ZSZ/FBsqpkrYAj0fEkkwqMTOrQIWMJJrDnDlzmD59OjNmzGDr1q1s2LCBc845h7vvvrvoPou9yus0ktNO4yXdWvTezcysVVx77bXU1NSwfPlypk2bxkknnVRSmECRI5SIeAt4JP0xMzMrboQi6RZJd6Svx2RakZmZtagTTjiBhx56qOR+ij3ltQ1Ylr4+qeQqzMys7BUbKFuAfdKbCQu98dHMzHZjxV7l9Q7wHsnswHOyK8fMzMpVk0YokrpIuh04PV10FzAi86rMzKzsNGmEEhHrJV0H9APeBo4A6p03y8zMKkcxp7wuAF6PiEeB5zOux8zMylQxgfIucKGkQ4GFwIKImJ9tWWZm1tz69etH586dqaqqol27dsybN6+k/pocKBFxraRZwJ+BocDxgAPFzKwMPfHEE3Tt2jWTvpocKJKuIZkNeAHJ6OTJTCoxM7OyVswI5TuSDiB5zO7pkg6OiC9nX5qZ2e7vD/f+mbdXbMq0z659OvHxzx3SaDtJjBkzBkl85StfYdKkSSXtt9j7UL4C/FdEeC4vM7MyNWfOHHr16sWaNWsYPXo0AwcO5Pjjjy+6v2ID5TbgIkl7kTxyd0HRFZiZVbBCRhLNpVevXgB0796d8ePH8+yzz5YUKMVOvfLPJGHUDrip6L2bmVmr2Lx58/tPaty8eTOPPfYYQ4YMKanPYkcorwEDgAci4l9KqsDMzFrcW2+9xfjx4wHYsWMHZ511FmPHji2pz2IDZRGwArhA0g8j4qMlVWFmZi3qoIMOYuHChZn2WWygHAKsBaaS3OhoZmYVrtjvUAaS3Mx4OVDadWZmZrZbKDZQugBXAlcAWzOrxszMylaxp7yuAQZGxBJJO7MsyMzMylNBIxRJVZJWSfoSQETURMTM9PXk5izQzMzKQ0GBEhG1wMvAwc1bjpmZlaumfIeyJ3CFpHmSpqc/DzRXYWZm1rzWr1/PhAkTGDhwIIMGDeKZZ54pqb+mfIdyTPrnsPQHIErau5mZtZqvf/3rjB07lvvuu49t27axZcuWkvprSqD0L2lPZmbWZmzYsIGnnnqKO+64A4AOHTrQoUOHkvosOFAi4o2S9mRmZn/niTumsuaNZZn22f3AgzjxvIZvEVy2bBndunXj/PPPZ+HChQwfPpwbb7yRvfbaq+j9FnsfipmZlbEdO3bwwgsvcNFFFzF//nz22msvrrvuupL6LPY+FDMzy0BjI4nmUl1dTXV1NaNGjQJgwoQJJQdKk0cokj5d0h4b7nuspCWSlkr6u/tblLgpXf+ipGGFbmtmZv9fjx496NOnD0uWLAFg1qxZDB48uKQ+ixmh/AfwYEl7zUNSFXALMBqoAZ6TND0iXslpNo5k2vwBwCjgZ8CoArc1M7McP/nJTzj77LPZtm0bBx10ELfffntJ/RUTKCppj/UbCSyNiGUAkqYBpwK5oXAqcFdEBDBXUhdJPYF+BWybmTsu+x7vdWjfHF2bWQUY/umPs6ZmdavW0L5KDB06lHnz5mXWZzGB0lz3nvQmecbKLjUko5DG2vQucFsAJE0inSG5b9++RRW6U1W81662qG3NzEKwU617G1/szH7/belL+Xwjn7pHXF+bQrZNFkZMJXmOCyNGjCjqb/SL/3llMZuZmQGwePFievTu2dplZK4tBUoN0CfnfTWwssA2HQrY1szMmlEx96G8lXkVieeAAZL6S+oAnAlMr9NmOvCF9Gqvo4G/RsSqArc1M7Nm1OQRSkSMbo5CImKHpEuAR4Eq4LaIWCTpwnT9FGAGcDKwFNgCnN/Qts1Rp5mZ5deWTnkRETNIQiN32ZSc1wFcXOi2ZmbWcjz1iplZBVqyZAlDhw59/2fvvffmhhtuKKnPokYoki6NiB+nrw+NiCUlVWFmZi3q0EMPZcGCBQDU1tbSu3dvxo8fX1KfTQoUSV2A64GBkrYCLwIXkH6XYWZm5WfWrFkcfPDBHHjggSX106RAiYj1wPmSPgWsBsYAvy2pAjOzCrb+wdfYtnJzpn126LUXXT5d+BPbp02bxsSJE0veb7HfoXyC5PLho0nmzzIzszK0bds2pk+fzhlnnFFyX8Ve5dUFuBK4guSUl5mZFaEpI4nm8PDDDzNs2DAOOOCAkvsqNlCuAQZGxBJJO0uuwszMWsU999yTyekuKPKUV0TURMTM9LWfPWJmVoa2bNnC448/zmmnnZZJf0UFiqRbJN2Rvh6TSSVmZtai9txzT9atW8c+++yTSX/Ffim/DViWvj4pk0rMzKysFRsoW4B9JLUHinuoiJmZ7VaK/VL+HeA9ksfuzsmuHDMzK1dNGqGkj9y9HTg9XXQXMCLzqszMrOw0+U55SdeRPMP9beAIfKe8mZlR3CmvC4DXI+JR4PmM6zEzszJVzJfy7wIXSrpB0vmSjsq6KDMza37XX389hx12GEOGDGHixIls3bq1pP6aHCgRcS3wZeC7wOvA8SVVYGZmLe7NN9/kpptuYt68ebz88svU1tYybdq0kvps8ikvSdeQPGZ3AbAgIp4sqQIzM2sVO3bs4L333qN9+/Zs2bKFXr16ldRfMc+U/46k75CMbk6XdHBEfLmkKszMKtTDDz/M6tWrM+2zR48ejBs3rsE2vXv35vLLL6dv377ssccejBkzhjFjSpv4pNgbG28DBgH7Az8tqQIzM2tx7777Lg888ACvv/46K1euZPPmzdx9990l9VnsjY3/TDL9SjvgRvw9iplZURobSTSXmTNn0r9/f7p16wbAaaedxtNPP80555xTdJ/FjlBeAzoCD0SEw8TMrMz07duXuXPnsmXLFiKCWbNmMWjQoJL6LDZQFgGzgQskPVdSBWZm1uJGjRrFhAkTGDZsGIcffjg7d+5k0qRJJfVZ7Cmvg0nuR5ma/mlmZmXm6quv5uqrr86sv2IDZUVEzJbUE1iTWTVmZla2ij3lNVZSNTAFuD7DeszMrEwVGyhdgCuBK4C/ZVaNmVmFiIjWLqFRTa2x2EC5huQKryVAbZF9mJlVpI4dO7Ju3bo2HSoRwbp16+jYsWPB2xT0HYqkKqAG+LeIuDUiatL3RMTkYoo1M6tU1dXV1NTUsHbt2tYupUEdO3akurq64PYFBUpE1Ep6meTqLjMzK0H79u3p379/a5eRuaac8toTuELSPEnT058HsihC0n6SHpf0avrnvvW0GytpiaSlkibnLP+hpD9JelHS7yR1yaIuMzMrXFMC5RhAwDDglJyfLEwGZkXEAGBW+v4D0tNutwDjgMHAREmD09WPA0Mi4gjgz8A3M6rLzMwK1JT7UJpzfHYqcEL6+k7gSZKryHKNBJZGxDIASdPS7V6JiMdy2s0FJjRjrWZmlkejgSKpb/oy7+UIOevXR8SGIus4ICJWAUTEKknd87TpDazIeV8DjMrT7ovAr4usw8zMilTICOVOkjBRA20CuAO4q74GkmYCPfKs+nYBNVDP/j8QcpK+DewAftVAHZOASZBMjmZmZtloNFAi4sQsdhQR/1DfOklvSeqZjk7qm86lBuiT874aWJnTx7kk3+l8Mhq4uDsippLMQcaIESPa7kXgZmZlptgbG7M2HTg3fX0ukO/qseeAAZL6S+oAnJluh6SxJN+5fCYitrRAvWZmVkdbCZTrgNGSXgVGp++R1EvSDICI2AFcAjwKLAbujYhF6fY3A52BxyUtkDSlpQ/AzKzSFTvbcKYiYh3wyTzLVwIn57yfAczI0+4jzVqgmZk1qq2MUMzMrMw5UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLRJsIFEn7SXpc0qvpn/vW026spCWSlkqanGf95ZJCUtfmr9rMzHK1iUABJgOzImIAMCt9/wGSqoBbgHHAYGCipME56/sAo4G/tEjFZmb2AW0lUE4F7kxf3wl8Nk+bkcDSiFgWEduAael2u1wPXAFEM9ZpZmb1aCuBckBErAJI/+yep01vYEXO+5p0GZI+A7wZEQsb25GkSZLmSZq3du3a0is3MzMA2rXUjiTNBHrkWfXtQrvIsywk7Zn2MaaQTiJiKjAVYMSIER7NmJllpMUCJSL+ob51kt6S1DMiVknqCazJ06wG6JPzvhpYCRwM9AcWStq1/AVJIyNidWYHYGZmDWorp7ymA+emr88FHsjT5jlggKT+kjoAZwLTI+KliOgeEf0ioh9J8AxzmJiZtay2EijXAaMlvUpypdZ1AJJ6SZoBEBE7gEuAR4HFwL0RsaiV6jUzszpa7JRXQyJiHfDJPMtXAifnvJ8BzGikr35Z12dmZo1rKyMUMzMrcw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTDhQzMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhQzM8uEA8XMzDLhQDEzs0w4UMzMLBMOFDMzy4QDxczMMuFAMTOzTCgiWruGViNpLfBGkZt3Bd7OsJxy4GOuDD7mylDKMR8YEd3qLqzoQCmFpHkRMaK162hJPubK4GOuDM1xzD7lZWZmmXCgmJlZJhwoxZva2gW0Ah9zZfAxV4bMj9nfoZiZWSY8QjEzs0w4UMzMLBMOlEZIGitpiaSlkibnWS9JN6XrX5Q0rDXqzFIBx3x2eqwvSnpa0pGtUWeWGjvmnHYflVQraUJL1pe1Qo5X0gmSFkhaJOn/tnSNWSvgv+t9JD0oaWF6zOe3Rp1ZknSbpDWSXq5nfbafXxHhn3p+gCrgNeAgoAOwEBhcp83JwMOAgKOB/2ntulvgmI8F9k1fj6uEY85pNxuYAUxo7bqb+d+4C/AK0Dd93721626BY/4W8P30dTfgHaBDa9de4nEfDwwDXq5nfaafXx6hNGwksDQilkXENmAacGqdNqcCd0ViLtBFUs+WLjRDjR5zRDwdEe+mb+cC1S1cY9YK+XcG+BrwG2BNSxbXDAo53rOA30bEXwAiohKOOYDOkgR0IgmUHS1bZrYi4imS46hPpp9fDpSG9QZW5LyvSZc1tU05aerxXEDyG045a/SYJfUGxgNTWrCu5lLIv/EhwL6SnpT0vKQvtFh1zaOQY74ZGASsBF4Cvh4RO1umvFaT6edXu5LL2b0pz7K611kX0qacFHw8kk4kCZTjmrWi5lfIMd8AXBkRtckvsGWtkONtBwwHPgnsATwjaW5E/Lm5i2smhRzzPwILgJOAg4HHJf0hIjY0c22tKdPPLwdKw2qAPjnvq0l+e2lqm3JS0PFIOgK4FRgXEetaqLbmUsgxjwCmpWHSFThZ0o6IuL9FKsxWof9dvx0Rm4HNkp4CjgTKNVAKOebzgesi+XJhqaTXgYHAsy1TYqvI9PPLp7wa9hwwQFJ/SR2AM4HpddpMB76QXi1xNPDXiFjV0oVmqNFjltQX+C3w+TL+jTVXo8ccEf0jol9E9APuA75apmEChf13/QDwcUntJO0JjAIWt3CdWSrkmP9CMiJD0gHAocCyFq2y5WX6+eURSgMiYoekS4BHSa4SuS0iFkm6MF0/heSKn5OBpcAWkt9yylaBx/wdYH/gp+lv7DuijGdqLfCYdxuFHG9ELJb0CPAisBO4NSLyXnpaDgr8N/534A5JL5GcCroyIsp6SntJ9wAnAF0l1QBXAe2heT6/PPWKmZllwqe8zMwsEw4UMzPLhAPFzMwy4UAxM7NMOFDMzCwTDhSzjEjqIumrOe97Sbqvmfb1WUnfaaTNf0o6qTn2b5aPLxs2y4ikfsBDETGkBfb1NPCZhu6TkHQg8POIGNPc9ZiBRyhmWboOODh9hsgPJfXb9RwKSedJuj993sbrki6RdKmk+ZLmStovbXewpEfSCRn/IGlg3Z1IOgT4W0S8Lalz2l/7dN3ekpZLah8RbwD7S+rRgn8HVsEcKGbZmQy8FhFDI+J/5Vk/hGRa+JHAfwBbIuIo4Blg12y+U4GvRcRw4HLgp3n6+RjwAkBEbASeBD6VrjsT+E1EbE/fv5C2N2t2nnrFrOU8kQbARkl/BR5Ml78EHCGpE8nDy/47Z0bjD+fppyewNuf9rcAVwP0kU2d8OWfdGqBXVgdg1hAHilnL+VvO650573eS/L/4IWB9RAxtpJ/3gH12vYmIOenptU8AVXXm3OqYtjdrdj7lZZadjUDnYjdOn7vxuqQz4P3nfR+Zp+li4CN1lt0F3APcXmf5IUDZTupo5cWBYpaR9LkwcyS9LOmHRXZzNnCBpIXAIvI/ivgp4Ch98ElfvwL2JQkVANIv6j8CzCuyFrMm8WXDZmVI0o3AgxExM30/ATg1Ij6f02Y8MCwi/q2VyrQK4+9QzMrT90geeoWknwDjSJ5rkasd8KMWrssqmEcoZmaWCX+HYmZmmXCgmJlZJhwoZmaWCQeKmZllwoFiZmaZ+H+yd8uFZ3nA5gAAAABJRU5ErkJggg==\n", - "text/plain": [ - "
" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "fig, ax = plt.subplots()\n", "swiftdiff['dv'].sel(id=plidx).plot.line(x=\"time (y)\", ax=ax)\n", @@ -174,22 +177,9 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAElCAYAAADDUxRwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAiMElEQVR4nO3dfZxWdZ3/8ddbRCeBRAUVHHFQMSVN1FlvysW7KNGU/FX+pDLvyqylbS0XJmvzrhI3+9W66rqablCtbGoitqQixE/TNLFQRGJFEJkYFVC8DRX97B/nYBfjNTPXdc33mplr5v18PObBdc75nnM+B/R6z/fcfI8iAjMzs87aorsLMDOz3sGBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8WsCEkXSvpZ/nmEpFck9atgO9dI+qf0FZr1PA4U65UkPSXpw63mnS7pt+VuKyKejoiBEfFWBeueExGXlNJW0k8kfafcfaRS6d+P2SYOFLM+QNKWvWEf1rM5UKzPkjRc0i2S1khaIenv22jXICk2fWHm682S9LykZZK+0M4+3ul1SDpSUrOkr0t6TlKLpDPyZWcDnwEm56fXbu+oRknvkTRN0guSlkiaLKm5YPlTkqZIehR4VdKWkpokPSnpZUmPSzopb7sPcA1wWL7/9fn8bSVNz/e/UtK3JG2RLztd0n2SfijpeeDCSv8trHfwbxTWJ+VfircDtwETgXrgbklLI+LODla/EVgMDAf2BuZIWh4Rc0vY9c7AtsAuwDjgZkkzI+JaSR8EmiPiWyXWeAHQAOwODABmF9nfROB4YG1EbJT0JPC3wDPAp4CfSdozIpZIOgf4fEQcXrD+v+b17g7sANwFtADX58sPAWYAOwL9Szh+68XcQ7HebKak9Zt+gKsLlv0NMDQiLo6INyJiOXAdcEp7G5S0K3A4MCUiNkTEQuDHwKkl1vQmcHFEvBkRs4FXgPe10bajGk8GvhcRL0REM3BFkW1cERGrIuIvABFxU0Ssjoi3I+K/gCeAg9s41n7A/wW+EREvR8RTwA9aHevqiPjXiNi4aR/Wd7mHYr3ZxyPi7k0Tkk4HPp9P7gYM33RqJ9cPuLeDbQ4Hno+IlwvmrQQaS6xpXURsLJh+DRjYRtuOahwOrCpYVvi56DxJnwO+RtazId/3kDb2PwTYiuz4NllJ1rtqb5/WRzlQrK9aBayIiFFlrrca2F7SoIJQGQH8OUFNrYf+7qjGFrLTYI/n07u2t01Ju5H1cI4BfhcRb0laCKiN/a8l61HtVrCP1sfq4crtHT7lZX3V74GX8ovW75HUT9K+kv6mvZUiYhVwP3CppDpJHwDOAn6eoKZnya5VlFrjL4BvSNpO0i7ApA62P4AsANYA5DcE7Ntq//WStgLIb5P+BfBdSYPyQPoa8LPOHab1Vg4U65PyL8sTgDHACrLfxn9MdgG6IxPJThmtBm4FLoiIOQnKuh4YnV/zmVlCjRcDzfmyu4Gbgdfb2nhEPE52DeR3ZOGxH3BfQZN5ZDcbPCNpbT7vK8CrwHLgt8B/Ajd09kCtd5JfsGXWO0j6EnBKRBzR3bVY3+QeilmNkjRM0ockbSHpfcDXyXpMZt3CF+XNatdWwL8DI4H1ZM+DXN3eCmbV5FNeZmaWhE95mZlZEg4Usx5I0mck3VVCu3eG2e8JunvEZOteDhSrefrr+0o2/YSkVwum/7aCbb5r+PtWy4+U9Ha+/ZclLd000GMF+9ps8EmAiPh5RHykku2ZdRdflLeaFxFPUzB8iaQA9o+IZVXe9eqIqJckYALZQI8P5s97lEQe8t16EfdQrFeTtLWkyyU9LelZZW9QfE++bIikX+UPEj4v6d78Ftyfkg0xcnveA5nc3j4iMxN4gezBxOMl/VHSS5JWSbqwoJ5NvZGzJD1N9jDhPfni9fn+DlOrl11Jer+kOXmdz0o6v43jPVTS/fkxPSLpyIJlp0tanveoVkj6TDt/Zz+StDr/+ZGkrfNlbQ7BX2Q7j0k6oWC6v6S1ksa09/dptcuBYr3dZcBeZE+b70k2sOG382VfJ3vSfCiwE3A+WT6cCjwNnJC/qfGf29tBHkInAYOBRWRPln8unz4e+JKkj7da7QhgH+CjwNh83uB8f79rtf1BZE/C30E2IOSewLuGys+HX/lv4DvA9sB5wC2ShkoaQDYa8fiIGAR8EFjYxiF9EziU7O9sf7LRiL9VsLxwCP6zgKskbVdkO9OBzxZMHwe05CM0Wy/kQLFeKz8V9QXg3IjYNELw9/jr8O9vAsOA3fLh5O+N8u6j3zQS8Fqyd5OcGhFLI2J+RCzKh4h/lOz9Ka2fXr8wIl4tccj3jwHPRMQP8iHzX46IB4u0+ywwOyJm5/ueAywg+yIHeBvYV9J7IqIlIha3sb/PkA2x/1xErAEuYvMh60sdgv9nwHGS3ptPnwr8tITjtRrlQLHebCiwDfCw/vpOlDvy+QDfB5YBd+WngprK3P7qiBgcEdtHxJiImAEg6RBJv1H2lsMXgXN49xDx5Qz7vivwZAntdgM+pc3fAXM4MCwiXiV7t8k5QIuk/5a0dxvbGc67h6wfXjBd0hD8EbGabKywT0gaDIwnzSCa1kM5UKw3Wwv8BXh//sU/OCK2jYiBAPlv+l+PiN3JBmH8mqRj8nU788TvfwKzgF0jYluyV+uqVZto43Mxq4A9StjvKuCnBcc6OCIGRMRUgIi4MyLGkfXK/kQ2lH0xq8nCaZMR+bxKTCPrOX2KbMj8FMP8Ww/lQLFeKyLeJvvS/KGkHSG7ziDpo/nnj0naMz819hLwVv4D7x5KvhyDyF7CtUHSwcCnO2i/hux0VFv7+xWws6R/yC+YD5J0SJF2PwNOkPRRZUPd1+UX0esl7STpxPxayutkp6neKrINyE7RfSu/9jKE7JpTpc+6zAQOBL5Kdk3FejEHivV2U8hOaz0g6SWyi9ubzvePyqdfIRvS/eqImJ8vu5TsS3W9pPPK3OeXgYslvUz2ZfyL9hpHxGvAd4H78v0d2mr5y2Tvnz+B7F3wTwBHFdnOKrLbl88nC6lVwD+S/X++BdlNCKuB58mu6Xy5jZK+Q3bt5VGymwz+kM8rW36N6Bay8cZ+Wck2rHZ4LC8zqypJ3wb2iojPdtjYapofqjKzqpG0Pdmtxad21NZqn095mVlVSPoC2Wm3X0fEPR21t9rnU15mZpaEeyhmZpZEn76GMmTIkGhoaOjuMszMasrDDz+8NiKGtp7fpwOloaGBBQsWdHcZZmY1RdLKYvN9ysvMzJJwoJiZWRIOFDMzS6JPX0MxM0vhzTffpLm5mQ0bNnR3KUnV1dVRX19P//79S2rvQDEz66Tm5mYGDRpEQ0MD2VijtS8iWLduHc3NzYwcObKkdXzKy8yskzZs2MAOO+zQa8IEQBI77LBDWb0uB4qZWQK9KUw2KfeYHChmZpaEA8XMrIf74Ac/WHT+6aefzs0339zF1bTNgWJm1sPdf//93V1CSXyXl5lZDzdw4EBeeeUVIoKvfOUrzJs3j5EjR9LTRot3D8XMrEbceuutLF26lEWLFnHdddf1uJ6LA8XMrEbcc889TJw4kX79+jF8+HCOPvro7i5pMw4UM7Ma0pNvT3agmJnViLFjxzJjxgzeeustWlpa+M1vftPdJW3GF+XNzGrESSedxLx589hvv/3Ya6+9OOKII7q7pM04UMzMerhXXnkFyE53XXnlld1cTdt8ysvMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM+sFzjzzTHbccUf23Xffd+Y9//zzjBs3jlGjRjFu3DheeOEFANatW8dRRx3FwIEDmTRpUrIaHChmZr3A6aefzh133LHZvKlTp3LMMcfwxBNPcMwxxzB16lQA6urquOSSS7j88suT1tCjAkXSsZKWSlomqanIckm6Il/+qKQDWy3vJ+mPkn7VdVWbmXW/sWPHsv32228277bbbuO0004D4LTTTmPmzJkADBgwgMMPP5y6urqkNfSYJ+Ul9QOuAsYBzcBDkmZFxOMFzcYDo/KfQ4B/y//c5KvAEuC9XVK0mVkrF92+mMdXv5R0m6OHv5cLTnh/2es9++yzDBs2DIBhw4bx3HPPJa2rtZ7UQzkYWBYRyyPiDWAGMKFVmwnA9Mg8AAyWNAxAUj1wPPDjrizazMwyPaaHAuwCrCqYbmbz3kdbbXYBWoAfAZOBQe3tRNLZwNkAI0aM6FTBZmatVdKTqJaddtqJlpYWhg0bRktLCzvuuGNV99eTeijFBvlv/X7Lom0kfQx4LiIe7mgnEXFtRDRGROPQoUMrqdPMrCaceOKJTJs2DYBp06YxYULrkz5p9aQeSjOwa8F0PbC6xDafBE6UdBxQB7xX0s8i4rNVrNfMrMeYOHEi8+fPZ+3atdTX13PRRRfR1NTEySefzPXXX8+IESO46aab3mnf0NDASy+9xBtvvMHMmTO56667GD16dKdq6EmB8hAwStJI4M/AKcCnW7WZBUySNIPsdNiLEdECfCP/QdKRwHkOEzPrS2688cai8+fOnVt0/lNPPZW8hh4TKBGxUdIk4E6gH3BDRCyWdE6+/BpgNnAcsAx4DTiju+o1M7PN9ZhAAYiI2WShUTjvmoLPAfxdB9uYD8yvQnlmZtaOnnRR3szMapgDxczMknCgmJlZEg4UMzNLwoFiZtYLlDN8/Zw5czjooIPYb7/9OOigg5g3b16SGhwoZma9QDnD1w8ZMoTbb7+dRYsWMW3aNE499dQkNThQzMx6gXKGrz/ggAMYPnw4AO9///vZsGEDr7/+eqdr6FHPoZiZ1bxfN8Ezi9Juc+f9YPzUslcrZfj6W265hQMOOICtt96602U6UMzM+qjFixczZcoU7rrrriTbc6CYmaVUQU+iWtobvr65uZmTTjqJ6dOns8ceeyTZn6+hmJn1Um0NX79+/XqOP/54Lr30Uj70oQ8l258DxcysF5g4cSKHHXYYS5cupb6+nuuvv56mpibmzJnDqFGjmDNnDk1NTQBceeWVLFu2jEsuuYQxY8YwZsyYJK8HVjbeYt/U2NgYCxYs6O4yzKzGLVmyhH322ae7y6iKYscm6eGIaGzd1j0UMzNLwoFiZmZJOFDMzCwJB4qZmSXhQDEzsyQcKGZmloQDxcysFyhn+Prf//737zx/sv/++3PrrbcmqcGBYmbWC5QzfP2+++7LggULWLhwIXfccQdf/OIX2bhxY6drcKCYmfUC5Qxfv80227DlltlQjhs2bEBSkho8OKSZWUKX/f4y/vT8n5Juc+/t92bKwVPKXq+94esffPBBzjzzTFauXMlPf/rTdwKmM9xDMTPrgw455BAWL17MQw89xKWXXsqGDRs6vU33UMzMEqqkJ1Et7Q1fv8k+++zDgAEDeOyxx2hsfNfwXGVxD8XMrJdqa/j6FStWvHMRfuXKlSxdupSGhoZO7889FDOzXmDixInMnz+ftWvXUl9fz0UXXURTUxMnn3wy119/PSNGjOCmm24C4Le//S1Tp06lf//+bLHFFlx99dUMGTKk0zV4+HoPX29mneTh6zM+5WVmZkk4UMzMLIkeFSiSjpW0VNIySU1FlkvSFfnyRyUdmM/fVdJvJC2RtFjSV7u+ejOzvq3HBIqkfsBVwHhgNDBR0uhWzcYDo/Kfs4F/y+dvBL4eEfsAhwJ/V2RdMzOroh4TKMDBwLKIWB4RbwAzgAmt2kwApkfmAWCwpGER0RIRfwCIiJeBJcAuXVm8mVlf15MCZRdgVcF0M+8OhQ7bSGoADgAeTF+imZm1pScFSrHRyVrf09xuG0kDgVuAf4iIl4ruRDpb0gJJC9asWVNxsWZmPUk5w9dv8vTTTzNw4EAuv/zyJDV0GCiSRpT4895O1tIM7FowXQ+sLrWNpP5kYfLziPhlWzuJiGsjojEiGocOHdrJks3MeoZyhq/f5Nxzz2X8+PHJaijlSflpZL2A9sY3DuAnwPRO1PIQMErSSODPwCnAp1u1mQVMkjQDOAR4MSJalI29fD2wJCL+XydqMDOrSWPHjuWpp57abN5tt93G/PnzgWz4+iOPPJLLLrsMgJkzZ7L77rszYMCAZDV0GCgRcVTreZJ2johnklWR7WejpEnAnUA/4IaIWCzpnHz5NcBs4DhgGfAacEa++oeAU4FFkhbm886PiNkpazQz68gz3/sery9JO3z91vvszc7nn1/2em0NX//qq69y2WWXMWfOnGSnu6Dysbw+B/xzsipyeQDMbjXvmoLPAfxdkfV+S/s9KDMzy11wwQWce+65DBw4MOl2Kw2UCZJeA+ZExNKUBZmZ1bJKehLV0tbw9Q8++CA333wzkydPZv369WyxxRbU1dUxadKkTu2v0kD5P2S35p4kac+I+HynqjAzs+Q2DV/f1NS02fD199577zttLrzwQgYOHNjpMIEKAyUingXuyH/MzKyblTN8fbVUFCiSrgIGRMTpkj4SEXclrsvMzMpw4403Fp0/d+7cdte78MILk9VQ6YONbwDL889HJ6rFzMxqWKWB8hqwbf4w4YiE9ZiZWY2q9KL888BfyEYHvi9dOWZmtSkiyJ6x7j3KfaNvWT0USYMl/QfwiXzWdOBdr4E0M+tL6urqWLduXdlfwD1ZRLBu3Trq6upKXqesHkpErJc0FWgA1gIfANocN8vMrC+or6+nubmZ3jbgbF1dHfX19SW3r+SU11nAioi4E3i4gvXNzHqV/v37M3LkyO4uo9tVEigvAOdIeh/wCLAwIv6YtiwzM6s1ZQdKRFwqaS7wP8AYYCzgQDEz6+PKDhRJF5ONBryQrHcyP3FNZmZWgyrpoXxb0k5kY3l9QtIeEfGF9KWZmVktqfQ5lC8C/x4RHsvLzMyAygPlBuBLkgaQvXJ3YbqSzMysFlU69Mrfk4XRlsAV6coxM7NaVWmgPAnUAbdFxNiE9ZiZWY2qNFAWA/OAsyQ9lLAeMzOrUZVeQ9kLWANcS/ago5mZ9XGV9lD2JnuY8Tzg7HTlmJlZrao0UAYDU4DJwIZk1ZiZWc2q9JTXxcDeEbFU0tspCzIzs9pUUg9FUj9JLZI+DxARzRFxd/65qZoFmplZbSgpUCLiLeAxYI/qlmNmZrWqnFNe2wCTJY0DVufzIiImpC/LzMxqTTmBclj+54H5D0Dved+lmZl1SjmB4teRmZlZm0oOlIhYWc1CzMystlX6HIqZmdlmHChmZpZE2YEi6YRqFJJv+1hJSyUtk/Su51uUuSJf/qikA0td18zMqquSHsp3k1dB9vAkcBUwHhgNTJQ0ulWz8cCo/Ods4N/KWNfMzKqokqFXlLyKzMHAsohYDiBpBjABeLygzQRgekQE8ICkwZKGAQ0lrJvMf513Ilv/aUU1Nm1m1iW2/ewZHHXK15Jus5JAqdazJ7sAqwqmm4FDSmizS4nrAiDpbPIRkkeMGFFRofHierZ9YWNF65qZ9QSvv5z+zSOVDg5ZDcV6Pq3Dq602paybzYy4luw9LjQ2NlYUjqdcd08lq5mZ9Wo9KVCagV0Lpuv56xAvHbXZqoR1zcysiiq5KP9s8ioyDwGjJI2UtBVwCjCrVZtZwOfyu70OBV6MiJYS1zUzsyoqu4cSEeOqUUhEbJQ0CbgT6AfcEBGLJZ2TL78GmA0cBywDXgPOaG/datRpZmbFKbthqm9qbGyMBQsWdHcZZmY1RdLDEdHYer6flDczsyQqChRJXyv4/L505ZiZWa0q6xqKpMHAD4G9JW0AHgXOIr+WYWZmfVdZgRIR64EzJB0PPAN8BPhlFeoyM7MaU+k1lCPIbh8+FKjKXV9mZlZbKg2UwcAUYDKwIVk1ZmZWsyp9Uv5iYO+IWCrp7ZQFmZlZbaooUCKimWwYFCLC7x4xM7OKbxu+StJP8s8fSVqRmZnVpEqvobwBLM8/H52oFjMzq2GVBsprwLaS+gOVvVTEzMx6lUovyj8P/IXstbv3pSvHzMxqVVk9lPyVu/8BfCKfNR141wBhZmbW95T9pLykqWTvcF8LfAA/KW9mZlR2yussYEVE3Ak8nLgeMzOrUZUEygvAOfkow48ACyPij2nLMjOzWlPJGxsvlTQX+B9gDDAWcKCYmfVxZQeKpIvJXrO7kKx3Mj9xTWZmVoPKfg4lIr4NvJ6v+wlJ1yWvyszMak6lDzbeAOwD7ABcna4cMzOrVZUGyt+TnS7bEviXdOWYmVmtqjRQngTqgNsiYmzCeszMrEZVGiiLgXnAWZIeSliPmZnVqErH8tqD7HmUa/M/zcysj6s0UFZFxDxJw4DnUhZkZma1qdJTXsdKqgeuAX6YsB4zM6tRlQbKYGAKMJnsmRQzM+vjKj3ldTGwd0QslfRWyoLMzKw2ldRDkdRPUoukzwNERHNE3J1/bqpmgWZmVhtKCpSIeAt4jOzuLjMzs3cp5xrKNsBkSQskzcp/bktRhKTtJc2R9ET+53ZttDtW0lJJyyQ1Fcz/vqQ/SXpU0q2SBqeoy8zMSldOoBwGCDgQ+FjBTwpNwNyIGAXMzac3I6kf2TvsxwOjgYmSRueL5wD7RsQHyIbV/0aiuszMrETlXJQfWbUqYAJwZP55GjCf7C6yQgcDyyJiOYCkGfl6j0fEXQXtHgA+WcVazcysiA4DRdKI/GN0sHx9RLxUYR07RUQLQES0SNqxSJtdgFUF083AIUXanQn8V4V1mJlZhUrpoUwjCxO10yaAnwDT22og6W5g5yKLvllCDbSx/81CTtI3gY3Az9up42zgbIARI0a01czMzMrUYaBExFEpdhQRH25rmaRnJQ3LeydtDefSDOxaMF0PrC7Yxmlk13SOiYiivam8jmvJxiCjsbGxzXZmZlaeSp+UT20WcFr++TSg2N1jDwGjJI2UtBVwSr4eko4lu+ZyYkS81gX1mplZKz0lUKYC4yQ9AYzLp5E0XNJsgIjYCEwC7gSWAL+IiMX5+lcCg4A5khZKuqarD8DMrK+rdOiVpCJiHXBMkfmrgeMKpmcDs4u027OqBZqZWYd6Sg/FzMxqnAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMTOzJBwoZmaWhAPFzMyScKCYmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJdEjAkXS9pLmSHoi/3O7NtodK2mppGWSmoosP09SSBpS/arNzKxQjwgUoAmYGxGjgLn59GYk9QOuAsYDo4GJkkYXLN8VGAc83SUVm5nZZnpKoEwApuWfpwEfL9LmYGBZRCyPiDeAGfl6m/wQmAxEFes0M7M29JRA2SkiWgDyP3cs0mYXYFXBdHM+D0knAn+OiEc62pGksyUtkLRgzZo1na/czMwA2LKrdiTpbmDnIou+WeomiswLSdvk2/hIKRuJiGuBawEaGxvdmzEzS6TLAiUiPtzWMknPShoWES2ShgHPFWnWDOxaMF0PrAb2AEYCj0jaNP8Pkg6OiGeSHYCZmbWrp5zymgWcln8+DbitSJuHgFGSRkraCjgFmBURiyJix4hoiIgGsuA50GFiZta1ekqgTAXGSXqC7E6tqQCShkuaDRARG4FJwJ3AEuAXEbG4m+o1M7NWuuyUV3siYh1wTJH5q4HjCqZnA7M72FZD6vrMzKxjPaWHYmZmNc6BYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl4UAxM7MkHChmZpaEA8XMzJJwoJiZWRKKiO6uodtIWgOsrHD1IcDahOXUAh9z3+Bj7hs6c8y7RcTQ1jP7dKB0hqQFEdHY3XV0JR9z3+Bj7huqccw+5WVmZkk4UMzMLAkHSuWu7e4CuoGPuW/wMfcNyY/Z11DMzCwJ91DMzCwJB4qZmSXhQOmApGMlLZW0TFJTkeWSdEW+/FFJB3ZHnSmVcMyfyY/1UUn3S9q/O+pMqaNjLmj3N5LekvTJrqwvtVKOV9KRkhZKWizp/3d1jamV8N/1tpJul/RIfsxndEedKUm6QdJzkh5rY3na76+I8E8bP0A/4Elgd2Ar4BFgdKs2xwG/BgQcCjzY3XV3wTF/ENgu/zy+LxxzQbt5wGzgk91dd5X/jQcDjwMj8ukdu7vuLjjm84HL8s9DgeeBrbq79k4e91jgQOCxNpYn/f5yD6V9BwPLImJ5RLwBzAAmtGozAZgemQeAwZKGdXWhCXV4zBFxf0S8kE8+ANR3cY2plfLvDPAV4Bbgua4srgpKOd5PA7+MiKcBIqIvHHMAgyQJGEgWKBu7tsy0IuIesuNoS9LvLwdK+3YBVhVMN+fzym1TS8o9nrPIfsOpZR0es6RdgJOAa7qwrmop5d94L2A7SfMlPSzpc11WXXWUcsxXAvsAq4FFwFcj4u2uKa/bJP3+2rLT5fRuKjKv9X3WpbSpJSUfj6SjyALl8KpWVH2lHPOPgCkR8Vb2C2xNK+V4twQOAo4B3gP8TtIDEfE/1S6uSko55o8CC4GjgT2AOZLujYiXqlxbd0r6/eVAaV8zsGvBdD3Zby/ltqklJR2PpA8APwbGR8S6LqqtWko55kZgRh4mQ4DjJG2MiJldUmFapf53vTYiXgVelXQPsD9Qq4FSyjGfAUyN7OLCMkkrgL2B33dNid0i6feXT3m17yFglKSRkrYCTgFmtWozC/hcfrfEocCLEdHS1YUm1OExSxoB/BI4tYZ/Yy3U4TFHxMiIaIiIBuBm4Ms1GiZQ2n/XtwF/K2lLSdsAhwBLurjOlEo55qfJemRI2gl4H7C8S6vsekm/v9xDaUdEbJQ0CbiT7C6RGyJisaRz8uXXkN3xcxywDHiN7LecmlXiMX8b2AG4Ov+NfWPU8EitJR5zr1HK8UbEEkl3AI8CbwM/joiit57WghL/jS8BfiJpEdmpoCkRUdND2ku6ETgSGCKpGbgA6A/V+f7y0CtmZpaET3mZmVkSDhQzM0vCgWJmZkk4UMzMLAkHipmZJeFAMUtE0mBJXy6YHi7p5irt6+OSvt1Bm8slHV2N/ZsV49uGzRKR1AD8KiL27YJ93Q+c2N5zEpJ2A66LiI9Uux4zcA/FLKWpwB75O0S+L6lh03soJJ0uaWb+vo0VkiZJ+pqkP0p6QNL2ebs9JN2RD8h4r6S9W+9E0l7A6xGxVtKgfHv982XvlfSUpP4RsRLYQdLOXfh3YH2YA8UsnSbgyYgYExH/WGT5vmTDwh8MfBd4LSIOAH4HbBrN91rgKxFxEHAecHWR7XwI+ANARLwMzAeOz5edAtwSEW/m03/I25tVnYdeMes6v8kD4GVJLwK35/MXAR+QNJDs5WU3FYxovHWR7QwD1hRM/xiYDMwkGzrjCwXLngOGpzoAs/Y4UMy6zusFn98umH6b7P/FLYD1ETGmg+38Bdh200RE3JefXjsC6NdqzK26vL1Z1fmUl1k6LwODKl05f+/GCkmfgnfe971/kaZLgD1bzZsO3Aj8R6v5ewE1O6ij1RYHilki+Xth7pP0mKTvV7iZzwBnSXoEWEzxVxHfAxygzd/09XNgO7JQASC/UL8nsKDCWszK4tuGzWqQpH8Bbo+Iu/PpTwITIuLUgjYnAQdGxD91U5nWx/gaillt+h7ZS6+Q9K/AeLL3WhTaEvhBF9dlfZh7KGZmloSvoZiZWRIOFDMzS8KBYmZmSThQzMwsCQeKmZkl8b+j0m08j6dtgAAAAABJRU5ErkJggg==\n", - "text/plain": [ - "
" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "fig, ax = plt.subplots()\n", "swiftdiff['dv'].sel(id=tpidx).plot.line(x=\"time (y)\", ax=ax)\n", diff --git a/examples/helio_swifter_comparison/tp.swifter.in b/examples/helio_swifter_comparison/tp.swifter.in index 58c816e2f..180f3efdd 100644 --- a/examples/helio_swifter_comparison/tp.swifter.in +++ b/examples/helio_swifter_comparison/tp.swifter.in @@ -1,13 +1,13 @@ 4 -101 +9 1.7496059999633410964 2.170163391141847864 -0.2537726760879844834 -3.0064589998644978604 2.1233488530690124423 0.6210068204130407379 -102 +10 3.0772345391474851262 -0.5509101822792066283 0.11666058691376969547 -0.08868603569822111026 2.7292630488987525612 -1.882742859645719835 -103 +11 0.13917497353384339354 -3.081984978409241016 0.69426813140927812196 3.105664664763373206 0.67090307556352112164 -0.2786153399455880027 -104 +12 -1.5389664718057010084 -1.5223401530194009545 0.23276670506845731357 3.2083305098906780644 -3.027143636331455024 -0.2998683055841925141