diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb b/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb index d82127068..7f6a7e57d 100644 --- a/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb +++ b/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb @@ -409,311 +409,147 @@ " fill: currentColor;\n", "}\n", "
<xarray.Dataset>\n",
-       "Dimensions:  (id: 19, time: 5)\n",
+       "Dimensions:  (id: 209, time: 3651)\n",
        "Coordinates:\n",
-       "  * id       (id) int64 0 2 3 4 5 6 ... 6000205 6000206 6000207 6000208 6000209\n",
-       "  * time     (time) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n",
+       "  * id       (id) int64 0 2 3 4 5 6 ... 6000395 6000396 6000397 6000398 6000399\n",
+       "  * time     (time) float64 0.0 1.0 2.0 3.0 ... 3.648e+03 3.649e+03 3.65e+03\n",
        "Data variables:\n",
        "    Mass     (time, id) float64 0.0002959 4.913e-11 7.243e-10 ... nan nan nan\n",
        "    Radius   (time, id) float64 0.00465 1.631e-05 4.045e-05 ... nan nan nan\n",
        "    J_2      (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n",
        "    J_4      (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n",
-       "    px       (time, id) float64 nan 0.3347 -0.4641 ... -0.9303 0.1796 1.422\n",
-       "    py       (time, id) float64 nan -0.2106 0.5473 ... -1.162 -1.411 -0.5606\n",
-       "    pz       (time, id) float64 nan -0.04792 0.03429 ... -0.02492 0.04703\n",
-       "    vx       (time, id) float64 nan 0.009457 -0.0155 ... 0.01503 0.007312\n",
-       "    vy       (time, id) float64 nan 0.02511 -0.01319 ... 0.001562 0.0129\n",
-       "    vz       (time, id) float64 nan 0.001184 0.0007137 ... 0.0005511 0.0005797
  • " ], "text/plain": [ "\n", - "Dimensions: (id: 19, time: 5)\n", + "Dimensions: (id: 209, time: 3651)\n", "Coordinates:\n", - " * id (id) int64 0 2 3 4 5 6 ... 6000205 6000206 6000207 6000208 6000209\n", - " * time (time) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n", + " * id (id) int64 0 2 3 4 5 6 ... 6000395 6000396 6000397 6000398 6000399\n", + " * time (time) float64 0.0 1.0 2.0 3.0 ... 3.648e+03 3.649e+03 3.65e+03\n", "Data variables:\n", " Mass (time, id) float64 0.0002959 4.913e-11 7.243e-10 ... nan nan nan\n", " Radius (time, id) float64 0.00465 1.631e-05 4.045e-05 ... nan nan nan\n", " J_2 (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n", " J_4 (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n", - " px (time, id) float64 nan 0.3347 -0.4641 ... -0.9303 0.1796 1.422\n", - " py (time, id) float64 nan -0.2106 0.5473 ... -1.162 -1.411 -0.5606\n", - " pz (time, id) float64 nan -0.04792 0.03429 ... -0.02492 0.04703\n", - " vx (time, id) float64 nan 0.009457 -0.0155 ... 0.01503 0.007312\n", - " vy (time, id) float64 nan 0.02511 -0.01319 ... 0.001562 0.0129\n", - " vz (time, id) float64 nan 0.001184 0.0007137 ... 0.0005511 0.0005797" + " px (time, id) float64 nan 0.3347 -0.4641 ... 1.069 0.6187 -0.5527\n", + " py (time, id) float64 nan -0.2106 0.5473 0.6083 ... 0.4382 1.45 2.533\n", + " pz (time, id) float64 nan -0.04792 0.03429 ... -0.07312 0.06803\n", + " vx (time, id) float64 nan 0.009457 -0.0155 ... -0.01209 -0.009084\n", + " vy (time, id) float64 nan 0.02511 -0.01319 ... 0.005512 -0.0004694\n", + " vz (time, id) float64 nan 0.001184 0.0007137 ... 0.001462 0.0007843" ] }, "execution_count": 4, @@ -740,44 +576,54 @@ "metadata": {}, "outputs": [], "source": [ - "swiftdiff = swiftdiff.rename({'time' : 'time (y)'})" + "swiftdiff = swiftdiff.rename({'time' : 'time (d)'})" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, + "outputs": [], + "source": [ + "swiftdiff['rmag'] = np.sqrt(swiftdiff['px']**2 + swiftdiff['py']**2 + swiftdiff['pz']**2)\n", + "swiftdiff['vmag'] = np.sqrt(swiftdiff['vx']**2 + swiftdiff['vy']**2 + swiftdiff['vz']**2)" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [], + "source": [ + "plidx = swiftdiff.id.values[swiftdiff.id.values < 10]" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "[,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ,\n", - " ]" + "array([2, 3, 4, 5, 6, 7, 8, 9])" ] }, - "execution_count": 7, + "execution_count": 20, "metadata": {}, "output_type": "execute_result" - }, + } + ], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 21, + "metadata": {}, + "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEzCAYAAADAeS+TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABtvElEQVR4nO3dd3hUVfrA8e+Znt5DAgkJIaEGooJSVFRQxN4Ae8Wyq2vbXcsW9WfvrljYVbFXii72joKLIEUIRHqVkN7b1HvP74+ZhCEEEkgmM4HzeZ55ZubOLecO5L5zyn2PkFKiKIqiKPtjCHYBFEVRlNCngoWiKIrSLhUsFEVRlHapYKEoiqK0SwULRVEUpV0qWCiKoijtMgW7AAciMTFRZmZmBrsYiqIoPcqKFSsqpJRJndlHjwoWmZmZLF++PNjFUBRF6VGEEDs6uw/VDKUoiqK0SwULRVEUpV0qWCiKoijt6lF9FoqiHNrcbjeFhYU4HI5gF6VHstlspKWlYTabu3zfKlgoihIyCgsLiYqKIjMzEyFEsIvTo0gpqayspLCwkH79+nX5/lUzlKIoIcPhcJCQkKACxUEQQpCQkBCwWpkKFoqihBQVKPZmd7k7tF4gvzsVLBRFUfZj7NixbS6/6qqrmDt3bsCPX2O3U1VRTmVDY8CPtT8qWCiKouzHzz//HLRja7pOfW0dUhiICQ8LWjlAdXAriqLsV2RkJA0NDUgpufnmm5k/fz79+vWjO2YZLa6rx6hrhMXEYDIE97e9qlkoiqJ0wH//+182bNjAmjVreOWVVwJe46hzuaGpCWkyExseHtBjdYQKFoqiKB2wcOFCLr74YoxGI71792b8+PEBO5YmJZW1tQgkSXGxIdHpr5qhFEVROqi7LtpFDU1Y3S6s4eFYAnCD3cFQNQtFUZQOGDduHB988AGaplFcXMwPP/wQkOPUuj3ojQ0gDMRFRwfkGAdD1SwURVE64LzzzmP+/PkMGzaMAQMGcMIJJ3T5MTy6pKy+gQhdIyY2FkOQO7X9qWChKIqyHw0NDYC3CeqFF14I6LF2OZyEOewYzWbCw4I7VLa10AlbiqIoh7EatwdPYwMGJPGxodGp7U8FC0VRlCBz6zrFTXZsbhfh4eEByRrbWSpYKIqiBJGUkl0ON+H2JgxCEB1Cndr+VLBQFEUJohqPhsNhx6xrREdHh1Sntr/QLJWiKMphwK3r7LI7iXTaMZvNhIfAndr7ooKFoihKEEgp2elwE+Z0YJCSmJiYkOvU9qeChaIoip+dO3dy0kknMXjwYIYOHcr06dMDcpwqt0aTy0WY20l4eDgWiyUgx+kq6j4LRVEUPyaTiaeffpqjjjqK+vp6RowYwSmnnMKQIUO67BguXafI6SLWZUcIQVRUVJftO1BUzUJRFMVPamoqRx11FABRUVEMHjyYXbt2ddn+vc1PLixuN0aPh+joaIxGY5ftP1BUzUJRlJB0/6e/sbaorkv3OaR3NPedNbTD62/fvp2VK1cyatSoLitDpdtDg0cjyWXHZDKFdKe2P1WzUBRFaUNDQwMXXHABzz77bJfd++DUdYqcbuLcTtD1kO/U9qdqFoqihKQDqQF0NbfbzQUXXMCll17K+eef3yX7lFKy0+7CpOuYnA7CwsKwWq1dsu/uoGoWiqIofqSUTJs2jcGDB/PnP/+5y/Zb7vLQqOnEux2IEL5Te19UsFAURfGzaNEi3n77bebPn88RRxzBEUccwRdffNGpfTo0nRKXm1ipobtcREVF9YhObX9Ba4YSQqQDbwEpgA68LKUMzIBmRVGUDjruuOOQUnbZ/qSU/O5wYQCsjiYMJhMRERFdtv/uEsw+Cw/wFynlr0KIKGCFEOJbKeXaIJZJURSlS5W5PNg1nVTpxqNpxMbH95hObX9Ba4aSUhZLKX/1va4H1gF9glUeRVGUrmbXdEpdbmIN4GlsxGazYbPZgl2sgxISfRZCiEzgSOCXNj67XgixXAixvLy8vNvLpiiKcjB0X/OTEUGk0wHQ4zq1/QU9WAghIoEPgduklHvdgSOlfFlKOVJKOTIpKan7C6goinIQylxuHJpOqkHidDqIiorCZOq5dysENVgIIcx4A8W7UsqPglkWRVGUrtKkaZS6PMSZjHgaGzAajT2yU9tf0IKF8PbwvAqsk1I+E6xyKIqidKXm5ieTEMRqLjweDzExMSE7qVFHBbP0xwKXA+OFEKt8j9ODWB5FURQcDgfHHHMMeXl5DB06lPvuu++Ati9xunFqkjSzkcaGBqxWa4/t1PYXtAY0KeX/gJ43fkxRlEOa1Wpl/vz5REZG4na7Oe644zjttNMYPXp0u9s2ejTKXR7izUb0pkakb1KjQ0HPrhcpiqJ0MSEEkZGRgDdHlNvt7tB9EZqv+clsECQKsNvtREZG9uhObX+HxlkoinLo+fJuKFnTtftMGQanPdbuapqmMWLECDZv3sxNN93UoRTlJU43Ll2SFW6hvqoKo9HYEnQOBapmoSiK0orRaGTVqlUUFhaydOlSCgoK9rt+g0ejwuUhwWJCOBx4fJMa9fRObX+qZqEoSmjqQA0g0GJjYznxxBP56quvyM3NbXOd5uYni0HQy2ykoroei8VySHRq+zt0wp6iKEoXKC8vp6amBvD2O3z33XcMGjRon+sXO924dUlfm4XG+vqWTu2emP9pf1TNQlEUxU9xcTFXXnklmqah6zpTp07lzDPPbHPdeo9GpctDksWEWdeobWoiIiICs9nczaUOPBUsFEVR/AwfPpyVK1e2u55Hl+x0uLAaBL2sZqoqKjAYDERFRXVDKbufaoZSFEU5CEVOF24pSQ+z4Ghqwu12H3Kd2v4OzbNSFEUJoFq3h2q3RrLFhA2oq6vDYrEQFhYW7KIFjAoWiqIoB8CjSwqdbmxGQS+LmfpDuFPbnwoWiqIoB2CX04VHekc/edxumg7hTm1/KlgoiqJ0UI3bQ41bo5fFjM1goLa29pDu1PangoWiKEoHuHWdQoeLMKOBZIsJu92O2+0mKirqkO3U9nfon6GiKMpB0DSNI488kjPPPBMpJYUONzrQ12YBKamrq8NsNhMeHh7sonYLFSwURVHaMH36dAYPHgxAtUejzqORYjFjMxqor69H1/VDvlPbnwoWiqIorRQWFvL5559z7bXXoktJkcNFuNFAksWE2+2msbGR8PBwLBZLsIvabdQd3IqihKTHlz7O+qr1XbrPQfGDuOuYu9pd77bbbuOJJ56grq4Ohy53Nz8BtbW1CCEOi05tf6pmoSiK4uezzz4jOTmZESNGUOfR8UhJqtWM1WjA4XDgcrmIjo7GaDQGu6jdStUsFEUJSR2pAQTCokWL+OSTT/jiiy9osNtprK/n9muu5q233qK2thaTyXTYdGr7UzULRVEUP48++ig7d+7k+7XreeL1NznppJN45513aGhoOOw6tf2pYKEoitJKpdtDg0cnwWzCIAQej4eGhgbCwsKwWq3BLl5QqGChKIrix6npFDndRJkMnH3yBD777LOWTu3o6OhgFy9oVLBQFEXxkb4pUgWQZrMghMDhcOB0OomKijrsOrX9qWChKIriU+7y0KTp9LFZsBgM6Lre0qkdERER7OIFlQoWiqIogEPTKXG5iTYZiTN5axCNjY1omnbYdmr7U8FCUZTDXnPzkwFIs5kRvk7t+vp6bDbbYdup7U8FC0VRDntlLg92TSfNZsHsyyBbV1d32Hdq+1PBQlGUw5pd0yl1uYkxG4k1e+9TdjgcOBwOIiMjMZnUvcsQ5Du4hRCvAWcCZVLK3GCWRVGUw48uJb87nBgRpFm9uZ+klAwYMICIiAisVismk4nly5cHuaTBF+yQ+QbwAvBWkMuhKEo3cTs1yn+vJyYpjPAYS1A7jktdbhyaJDPMgsngLUdjYyNSSr799lvS0tKCVrZQE9RgIaVcKITIDGYZFEUJPKlLijbVsH5JMVt+Lcft1AAIizKTlB5FYnokielR6DYdKWW3BJAmTaPM5SHObCTG1/ykaRr19fUIIbDZbAEvQ08S7JpFu4QQ1wPXA/Tt2zfIpVEU5UDUlDWxYUkJG34pob7SgdlmJHtkMpnDEqmvclBR2EDFznpWfbcTXZMcfVkcFTsbMFkMNLz4DO4tGxEGEELQFeHDOngQKX//u7f5ye7CLAS9rbvnpKirq0NKidFoZOLEiQghuOGGG7j++uu74Og9W8gHCynly8DLACNHjpRBLo6iKO1w2j1sXl7KhiUlFG+pRQhIGxzP6HOy6HdEEmbL3ndBa26dquJGiqt2YIs043FpeDw6mkdvWUcIX9AwCO9rw8EHkBKnG6cu6RdubWl+cjqd2O12IiMjWbRoEb1796asrIxTTjmFQYMGMW7cuIM82qEh5IOFoiihT9d0dq6vZsPiYrbmV6C5deJSwhlzXn8GHJNCZNz+71Mwmg0k9Y2iotFIVLy3+Sf2kXvRPDoel+4NHi4dj1tD13b/ZjSaDJgsBkwWIyaz99lo2v8gz0aPRrnLQ7zZRLTv5jspJbW1tRiNRiIjI1uGyyYnJ3PeeeexdOlSFSyCXQBFUXquyl0NrF9SwsalJTTVurBGmBgyNpWBY1JJzojqVN+DEAKT2YjJbIQIc8tybwDR8Lh3BxFnk6flc4PRt11zELEYMJoMCCHQfDffmQ2C3rbd+2xsbMTj8RAXF4fdbkfXdaKiomhsbOSbb77h3nvvPejzOFQEe+js+8CJQKIQohC4T0r5ajDLpCjK/tkbXGxc6m1mKv+9HoNBkDEsgUGjU8nITcBoDuztW0aT9+LvX1fRdbm79uF7bqpztXwuhMBkMeAwgFFIekVaMUhA7O7Utlgs2Gw2tm3bxnnnnQeAx+PhkksuYdKkSQE9p54g2KOhLg7m8RVF6RjNo7NjTSXrlxSzY00lui5J6hvFcVNzGHB0L8KiLO3vBJCaRuPPP9O4ZAlhQ4cSPno0pvj4TpfPYBBYbCYsfgOYpC69tQ+3r/bh1BAOjVgJziY75YDJbMRjsCOlJCIsCikhKyuL/Pz8TpfpUKOaoRRFaZOUkvLf61m/uIRNy0pxNLoJj7YwfEI6g0ankNAnssP7cpeWUvPhh9TMnYunqNjbWy29fQ/WIYOJHDuWiLFjkbGxXVZ+YRCYrUbMViOalOxsdCAw0N9iQfc1YbmcLjy6C6GZqa9wUo9zz34QX19Ie/0ghwMVLBRF2UNDtZONS0tYv6SE6uJGjCYD/Y5IZNDoVNIHx2EwduzCKTWNhoULqZk9h4YFC0DXiRg7hl533knkiSfi3LCBxsWLaVz0M5VvvkXlzFfxvPgizvBwjBERGCIjETZbl9xzUeR049Yl2eFWLCYjWIzIcBP2inoMGEhMikdz6wfVD3K4UMFCURTcLo1tq8pZv6SEwnVVSAkpWTGceOlAskckYw03t7+T5n0VFVHz4UfUfPghnpISjImJJFx7LbGTL8Did69UWF4eYXl5JP7hD+iNjTQtX852qxXcbtylpVBaijCaMER6A4chIgKDpWPNXf7qPBpVLg9JFhMRpt3DdpuamnC73cTGxrZ0pB9MP0hz8GgOJodqAFHBQlEOU1JKijfXeu+qXlGGy6ERGW9lxGmZDByVQmyv8I7vy+OhYcECqmfPpnHhTwBEHHccvf7xd6JOPBFh3n+wMUREEHnCCRjXrcOak4PudqM3NqI3NKA3NKDV1gIgLBaMkZEtwUO0M3OdR5cUOlxYDYIUq9+IKk2jrq4Oi8VCWFhY22XaXz+I31BeR4MbKXcP5/UfxtscTAyGzgUQu8uDzWwMaiBSwUJRDjN1FXbWLylhw5Ji6iocmKxGso9KYtDoVHrnxCIO4MLmKiykZu5caj/8CE95OabkZBL/+Adizr8AS1qfgy6jwWzGEBsLsbFIKZFOpzdoNDbiqamBqipAYAgP89Y4IiMxhIUhDHs2kRU5XbilJCfcisHvQltfX4+U8oAnNfLvB2kmpdzrfhCX3YOj0d2yjrcfxLjHPSEd6QdxazrFtQ5qmlxkJIQTE3bgNauuooKFohwGXHYPm38tY8OSEoo21YCAtIFxHHNmP7KOTN7j4tce6XZTP/8HambPpvHnn0EIIseNI3bqFCLHjUN0cUpvIQTCZsNgs2FKTETqOrrd3lLr8JSXQ3k5wmDYHTgiIqg3mqh2ayRbTYT71UBcLhdNTU1ERERgbqfG09Hytb4fREqJrsk9aiAel4azaXcA2V8/iC4llQ0uSuscSCA52kaUtfNl7QwVLBTlEKXrkl3rq1m/pJitK8vxuHVie4Uz6pwsBo5KablTuqNcO3ZQM3cuNf+dh1ZRgSklhcSbbiL2gvMxp6YG6Cz2JgwGjBERGCMioFcvpMeD3tiI5mu20urrATAYTfQJCyM2NgZpNCBMppY7tQ0GA1FRUfs8Rk1NDddeey0FBQUIIXjttdcYM2ZMx8soBEaT8N4P4teap2t7dqJ7XBpNdZ49thMmgV3XcUhJlNVISlwYVnPHg3mgqGChKIeYquJGNiwpZsMvpTTWOLGGmxg4JpVBo1Po1S/6gJpddJeLhu++o3rOHJoWLwGjkcgTTyRu6hQijjuu3T6D7iBMJowxMRhjYgBvmcuraxCNjUTZm3A31OMGDDYbrogI3FISEx2NwbDvZqBbb72VSZMmMXfu3JaaSFcwGA1YjAYstt2X3uZ+EKfDQ12jG92tYwGsCLDr1DoaMZmNRMZZ99iuu6lgoSiHAEeDm03LS1m/uJiyHfUIg6Dv0HiOm5JD5vAEbxPJAXBu20bNnLnU/ve/aNXVmHv3Jum2W4k573zMvZIDdBZdo04YKAmLJCU2DpvFhLQ70Bob0BoaaNA0jLqOYccOnOHhGCIjMUZEIsJ2D9Gtq6tj4cKFvPHGGwBYLBYsBzEKq6OkgBqXm9IGp7fJKc5KYqQVqe3ZD9LZTvLOUsFCUXooTdP5vaCS9UtK2L66Al2TJPSJ5NjJ2Qw4JoXw6AO7wOlOJ/XffEvN7Nk0LVsGJhNR48cTO2UKEceO3avzONB+mr2Rip0NB7SNRNKk6RiEYGur8no0D1FJZk48rx9GswW9sQFPaSkeShFGI4aISAyREWzesoWkpCSuvvpq8vPzGTFiBNOnTyciIqIrTw+ABoebXTUOnB6NaJuZ1Fgb1ubhvYa982IFkwoWitKDSCmp2NnA+iXFbFpWir3eTViUmWEnpjFoTAqJaftuh98X5+bN1MyZQ+28j9FqazGnp5P05z8Te965mJKSAnAWgePUfXeFt/oV7u1w1jCZwrDFxkKsb7nb3dLXoTc0otXV0rR1K7/++ivP3HsvY6ZP5/Z77uGxxx7jwQcf7LJyuj06xbV2auxuLCYDmQkRRIeFRlDYFxUsFKUHaKx1+pL3FVO5qxGDSdBvuO+u6qHxGDt4V3Uz3eGg7quvqJk9B/uvv4LZTNTJE4ibOpXwUaO6vRbRluOnDjig9avcHnbaXaRazST7jRySUlJZWYnb7SY5ec8mNGE2Y2o1RDdTGEhLSWFkVhaunTs5e8QInn7jDdylpfscottRupRUNDgpq3MC0CvaRlKkNehNTB2hgoWihCiPW2NbfgXrF5ewc20lUkKvftGccPEAskf2wnYQzROODRupmT2b2k8/Ra+rw5KRQfIddxBz3rldktAvWFy6zi6Hi3CjgSTLnpc1h8OBy+UiJiYG43465JuH6KYNHUJ6v35sF4Kcfv1Y8MYbDM7OxlNeAeXl4Bui25KSxGrt0KCB/TY59QAqWChKCJFSUrK1jvVLitm8vAyX3UNknJUjT81g0OgU4lIOvN1cb2qi7ssvvbWI/HyE2UzUqacSO2UK4ccc3ePTU0jpvUtbAn1tlj3OR9d1amtrMZlMhId3/I70559/nssuvxyXy0VWVhavv/46tujolrvKtYYG3L4husJk8tY4mlOStLp3w+VrcqrtQU1ObVHBQlFCQF2lnY2/eJP31ZbZMVkM9D8ymYGjU+gzMO6gmikc69ZRPXs2dZ9+ht7QgKV/f5LvvouYc87BFBcXgLMIjiq3Rr1Hp7fNjLVVc1xDQwO6rhMXF3dAQfGII45g+fLley03RkdjjI7GjHeIrt7Q4A0g9fVoNTUAGKxWb40jIpJqaaK0wZtPqic1ObVFBQtFCRKXw8PWleWsX1LMrg01APQZEMuISZn0PyrpoMbUaw2N1H3xOTWz5+AoKEBYrURPOpXYqVMJO+qoHl+LaM2p6xQ5XUSYDCSa9/y+PB4PDQ0NhIWFYbXuf1rXg2GwWDDEx0N8vLe/w+HYnZKkqgoqK4lAkG6xYo2OwmIy05O/fhUsFKUbSV2ya2M165eUsGVlOR6nRnRSGMec1Y+Bo1KITmw7qd1+9ykljoLfqJk9m7rPP0dvasKak0Ovf/yDmLPParlZ7VAjpWSnw/urPb1V8xNAbW0tQoiW+bQDSQiBCAvDY7ZSYginzhhDtHSTaPBgczShV5TjrCj3DdH1y6IbgCAWKCpYKEo3qCltYv3iYjb8UkJDtROLzciAo3sxaHQKKf0PLJldM62+nrrPPqN6zhyca9chbDaiTz+duKlTsOXlHXK1iNYq3B4aPTppNgvWVqOTHA4HTqeT6Ojo/XZqd5XWo5ySY8JIioxtaXKSHg9ac5NVQwNaXR0AwmzBEBnhzaQbEdHlebW6UuiWTFF6OEejm80ryli/uJjSbXUIAelD4hl7fjb98hIxWQ78IialxLF6tbcv4osvkXY71sGDSbnvXqLPPBPjfvIdHUqcmk6x002UyUB8q7vT/Tu1A3EjXWv1DjdFfqOcesfavBMs+REm055DdJv7Oxoa0Wtr0aqrATCEhbXcHGgIDw+JIczNVLBQlC6kazq/r61i/WLvXdWaRye+dwRjzu/PwGNSiIg9uGYHra6O2k8+pWb2bJwbNyLCw4k580xip07Fljv0kK9F+JNS8rvDhQFIa6P5qbGxEU3TSEhICOj34j/KyWoy0C8xgihb+6OchBAIq9XbBJWQ4L1hsMmO3ujLoltZARXlIAwYIsJb5u/o6BDdQFHBQlE6SeqSyqIG1i8pYePSUux1LmyRZoYe35tBY1JJTI88qD9yKSX2lSupmTWbuq++Qjqd2HJzSbn/fqLPOANjZOB/NYeicpeHJk0nPcyCpXVKD4+H+vp6bDZbQDq1wdfkVO+krN7b5JQSbSMxas/5Mg6EEAJjRDjGiHBITkZqmq+5qhGtsQF3SQkAlvT0oPY/qWChKO2QuqSp3kV9pYO6Srvv2UG930Pz6BiMgsxhiQwcnUJGbkKHJrdpi6e6mrpPPqF6zhxcm7dgiIgg5rxziZ0yhbChQ7v47HoWh6ZT4nITbTIS18YNbXV1dZ3u1N6wYQMXXnhhy/utW7fywAMPcNttt+3R5BQTZiY1Zu8mp84SRuOeQ3TdbvSGBgzd0KS2PypYKIc9qUua6ly+AGD3PlftHQz8hUWZiYq3kdAnkn7DE4lNCadfXiJhkQeXnVRKSdOyZdTMmUv9118jXS5secNJffghoidNCvqFIhTs2fxk3qu25nA4cDgcREVFYepER/HAgQNZtWoV4J1+tU+fPpxx1jnsqGz0NTkZO9zk1BUMZjOGELgvRgUL5ZDXZjBoFRh0j9xjm7AoM1EJYSSmRdIvL5HoBBtRCWFEJdiIircd0Mxy++OpqqJ23sfUzJmDa9s2DFFRxE6ZQuzUKdgGDuySYxwqylwe7JpORpgFc6vmJykldXV1GI1GIiMju+yY3373HX0z++G0xeN0eDrd5NSTqWCh9HidCQZJ6VFkHZEUsGDQdnl1mpYu9d4X8e134HYTdtRRpF5/PdGTTsUQduD3WhyKfnjjZcp2bAVAl2DXdIxCsN6494Va0zQ0j4bJbNrvpEbJGVmcdNX1HTp+vcPNK2+8w4QzzyPKZiI1JgzLQTYtHgpUsFBCXpcHgwQb5oMYttpZnooKav77X2rmzMX9++8YYmKIu/gi4qZMwZqT0+3l6TEkOHQdxN6pxwGkBM2jeefg7oKhpi6PRlGNg8q6RuZ//QVPPPYoGQmqGVAFCyXopC5prHW1HQgqHdRX94xg0Bap6zT+vJia2bOpnz8fPB7CR44k6eY/ETVxYo+6g7e7NdcAip0uypweMsMsxJj3vmRVV1djt9tJTk7uVF+FrkvKG5yU+0Y5FSxdwMiRI+ifkXbQ+zyUdOibFUJMk1K+6vfeCPxTSnl/wEqmHDI6FQz6RpF1ZOgGg31xl5ZR+9+PqJn7Ie7CQoyxscRffjmxUyZjzcoKdvF6jCZNo8zlIc5sbDNQOJ1O7HY7kZGRnQoUdQ43RTV2XB7dN8opjIfmfcglF1/cmeIfUjr67U4QQlwATAMSgNeBBZ09uBBiEjAdMAIzpZSPdXafSvc7qGAQbSE6wUZSRhT9j0pqCQTRCTYi40M/GLRFahqNixZRPXs2DT/8CJpG+OjRJP/5diJPPhlDAOdxPhTpUvK73YVZCHpb9/7upJTU1tZ2qlO7ucmpzrHnKKempia+/fZbXnrppc6exl6klCABXXpf64CUSH33cu97v9cSjNEWDEH8u+hQsJBSXiKEuBBYAzQBF0spF3XmwL7ayYvAKUAhsEwI8YmUcm1n9qt0PV2XNNU62w4EzX0GWvcGA+8fmQ6a5n2taUhdB11Hapr3M133LtM07x+irrX/WfO+NB2k7n3WW++7ef3mdTTcOwup+e9HeIqKMSYkkHDN1cROnowlI6NT53k4K3G6ceqSfuFWTG30VTQ2NuLxeIiLizvgvorWTU4pMTYSI3ePcgoPD6eysnKPbQ7oIi+lbz3/5bu36RAhwOC9aQ+D8B4jiDraDJUD3Ap8CAwGLhdCrJRSNnXi2McAm6WUW33H+AA4B+jyYPHq7fejW0JzhIkgNIfgmaSRFFcvImQ8urAhES3/VyUgdCexup0E3Y5Bb8Sk2TFqjZh8D2OdG7FTIqT3giqkxK1LKqWkUuotf0D4Pt/9R6W3PEvfH5x3PV9waHnG94tL95ZI+j7b67VvHe8B9l7H936vddraVzsM6RqmU3UM/YpoML5Iwycv7rmV70Lkv0zu57OW962Gae6x/b4+81u+v30e0PH2V/4O7tP/M81gQBei5VkXgphBT1JXCQ5hodKYTJRsQq/7nZpWfyZSCuyOSAwGDad9B04H++YtAEIaAIFLs9LojkBKIwkGNxFGJ8amBuxNBoQUIA3ev0vfc/N23uf2SSQIHSm8/5+kkEh0EBJp8v5/lWLv59bLWl8aZL2JOGt2h8oQCB1thvoUuElK+b3w3glzO7AM6MztpH2AnX7vC4FRrVcSQlwPXA/Qt2/fgzqQNIVRYXMd1LaBFOQfCvvlxsP28DoStEj6e3rR35OEDQtC+AKcsAI2hEgIdlG7T3NgQvc9ZMt7IXRAQwgNSj0INMD3LJrfa3sux+Nd3/d+92ce37q+5cLT9rZ77dtvO+G/bqt9iz3XBX33+yD/dlnX302000E0DpJkHd4rpgGJEYkBhAGJgTq8AwNidDA6IkEagOaH2P26ZfluYUBL0gzd6n3gd5FvvrALHYlEN7R9ce/oRf5Q0dFgcQxwrRDiJrx/Hf8DLurksdv6Sve6fkopXwZeBhg5cuRBXV+vffLOg9nssCWl5B+vfIFrzQ8YwhtZGtbAMttW+mZlkTlkMLEZfXDgodHjosntoNHjxOFx4XC7cLidODQnTo8bt9uFU3Ph8rjxaC7cmgdNc6PpbjyaB01zoesaUveg6R6k7gHd42si8iCkhtA1kBoGqWPQdYTUMEoNAxKD75efwferz4DAgPfXoe83IQbfcu/r1st86woTRswYhAkDJozCiAETJkwYhQmj79mEEaMwYhImTMKIESNmvK/DjOH0j84ixmDCLEFqvmYITYKme5sgNNnyrOu+5X7LWp416Wvyaq5ZdeM/vgGEQYBReJ9b3rP7vdFvuf+z/zoGWq3DHs8gkU4P0qmhOzR0p450STQRi0P2BwmGti4REtxoOAxNhEkLRmn1LhS6r+7rDeTNwY/mIN4S4PV9v259OP+qNMJbeOF7lr7n5sAkBEizrwblt17r7Vp/fiDrBVlHg8VMoB543vf+YmAMMLUTxy4E0v3epwFFndif0kWEEPzfNadx2avxlG0oYIpWQK3DQaGusWPLFmw2G0OHDiUvL4/09KxuyYQppcQjwaXrOKXErnlo8rio97hocttp9Lho9Dhp8riwe5w0aU4cmguHx4ldc+LUXL6H97Vbc+PSnbh0N27NhUdvfjSi6W7fw4Wuu9GlG6m7kNKDlG6QHoR0+Z49CNwAuPV+1MXdRFx4H/qFWcgMs9IvzEpWuNX32kJsGyN62j13fR+BpznA7OMz9OYApLcKXtLXX+P9fPeyvT/bY1+a33Lfds1lkG4dHK3L47evPcrjvQoLswFhM2GwWRFWI4YYE7pRp85kwGo0EGE2eoONAAzC13YPtbVVGHQDMYlxCKNhv///6uzeUU4eTSMmzERKlAWzgd1NoHs0Ser7WNa66dK/ybK5aVQD6WF3M6a+5zE6TUB8FtgCP5HTPksgO9DZIoTIl1LmtbfsgA4shAnYCEwAduFt1rpESvnbvrYZOXKkbGteXCUwqhpdnPPi/3C4NKYfIyn45AOqmxxYM3NoMJhbOhfz8vIYPnw48fHxwS5yt9CkxKlLXLqOS5fYNY35O+czY/mDaFJncOafcESMYZvdSZHTvce2cSYj/cK9QSQzzEI/X0DpF24lzmQ8LFKNN3cUi1ad1vUejWUFv5GancOACBvGNr6LxsZGamtriY2NJTw8fJ/HcHo0iv1GOfWJtRHZTbmc9uLXf7ZXENlfAGrdjxYeDyZbu4dbt24dgwcP3mOZEGKFlHJkZ06jo8HiDeA/UsolvvejgCullDd26uBCnA48i3fo7GtSyof3t74KFt1vY2k958/4mb7x4XwwbQQbfviKX/47C6fDSfLRY3FHxfN7YSHg7VPKy8tj6NCh2Gzt/6c+1BQ1FHHXwrtYVb6Kc7PP5W/H/A1hsLHD4WR7k4utdifb7U62+R67HO49fnPGmIx7BZDmoJJoNh3ygeQv639nkr2ascNyiWgjk6umaZSVlWE2m/c5V0XzKKeyeicC6BVtJSHy8MrlFOxgsQ4YCPzuW9QXWEfzQDAph3emEB2lgkVw/LC+jGlvLmPikBRmXHoUzqYGfvloFqu+/gxhMDJ04hkY+2Ty29q1VFRUYDKZGDhwIHl5efTv379bprUMFR7dw4xVM5i5ZiaZMZk8Oe5JBsa3nRDQqevssLv8AoiLbU3e14UOF/55bqOMBm/g8AsgWb6gkmTp+YHk+8o6Ll29lf8mmBgzPLfNdWpqamhqaiIpKQmzee9aQp3dTVGt98a62DALqTE2zAeZy+lf//oXM2fORAjBsGHDeP3113vMD6BgB4uM/X0updzRmUJ0lAoWwTPzp6089Pk6bh6fzV8mei9+tWUl/O+Dt1m/aAFh0TGMueBiEocMZ01BAWvWrMFutxMREcGwYcPIy8sjJSWlx1/UOuqX4l+4+6e7qXPWccfRd3DhwAsP6Nxdus5Oh2uPALLNVzP53eFqbvYHINxooJ9/jSTM20eSFW6lVw8IJDVuDycu3UC0ych/onSGDBmy1zoul4uKigoiIiKIaTUBkH+Tk81kpHcnm5x27drFcccdx9q1awkLC2Pq1KmcfvrpXHXVVQe9z+4UqGDR0ZvyuiUYKKFr2nH92Fhaz/PzN5OdHMk5R/QhJjmFM265gxGnn8OCd19j/uv/IS61D8dfciUT//IXNm/eTH5+PkuXLmXJkiUkJyeTl5fHsGHDOjU5TU8wKnUUc8+ayz8X/ZOHf3mYJcVLuH/s/cRYY9rfGLAYDPQPt9E/3ObNmeDHrUsKHa6WALLN7mRbk4t1DQ6+qqjF/2b5MINhd9NWq76SVKs5JJpn/rlpF+VuN28O74co3PtS03yntsFgIMpvjnFdl5TVOylv8DY5pcaEkRBp6ZJz8ng82O12zGbv3dy9e/fu9D57ug7VLEKFqlkEl8ujc9nMX1hVWMPsG8ZwRHpsy2dSSrb+upSF775B1a6d9Bk0hHGXXkPvAYNoamrit99+Iz8/n8LCQoQQZGVlkZeXx6BBg7AcwmkwdKnz9tq3efbXZ0kMS+SJcU9wZPKRATueR5fscrpamrW2+9VKdthduPz+3m0GQYZvlFbLyC1fU1dvq7nNDuau9mV5DVcXbOf2jF7clZW6x6/imk+34CpqRNc1PB4PJpMJg8EISDRd4vLo6BJMRoHFZOjQDa6W3hHEntW/3fWmT5/OP/7xD8LCwpg4cSLvvvtuZ0+12wS1GSpUqGARfJUNTs55cRFOj84nfzqW1Jg974zXNY2CH75l0ex3aKqtYcCoYznukiuJS/H+MquoqGD16tXk5+dTW1uLxWJhyJAh5OXlkZGR0SUppkNRQUUBdyy4g+LGYm484kam5U7DaOjevhxNSoqc7pYAsrvD3cUOuxOHvvtaYBGCDL+mLW9fifd9H6ulzfQbB6rS5eGEpetJsZr5YkQOFoOhjWDRgMvlRggwm81ICU6PjqZLDAawGI0YD6AsHQkW1dXVXHDBBcyaNYvY2FimTJnC5MmTueyyyzp1vt1FBQtUsAgVG0rqueDfP5OREM6cP4wh3LJ3a6bLYWf5p/9l+acfoXk85E08jdHnX0R4tLcZRtd1fv/9d/Lz8/ntt99wuVzExMQwfPhw8vLySExM7O7TCrgGVwMPLHmAL7d9yaiUUTxy/CMkhycHu1iAN2lfidPtCyCulv6RrU3eZ7tfIDELQV+brzYSvmdfSZrNgrmDF+/rCrbzVUUt34wcwOBI74+O1he62tpaGhsbSUhIpNqhU97gxAAkR9u6rMmptTlz5vDVV1/x6qveRNtvvfUWS5YsYcaMGV1+rEBQwQIVLELJ/PWlTHtzOZOGpvDiJUdh2McFoqG6isVz3mPN/G8w22yMOm8qR552FmbL7nkcXC4XGzZsID8/ny1btiClpE+fPuTl5ZGbm7vf8fQ9jZSSeZvn8ejSR7EZbTx83MMcn3Z8sIu1X1JKSl0eX9+If2e7dzhwk7Z73JZRQLqtdWe7haxwK+k2CxZfzXFeaTV/WLuDv/VL5dbMXi3b+1/o3G435eXlmCw2qtwmXJpObLhvlJMxcDXQX375hWuuuYZly5YRFhbGVVddxciRI7n55psDdsyupIIFKliEmlcWbuXhL9Zxy/hs/jxx//NFVxb+zsL33mDriqVEJSRx7IWXMeT4kxCtmp3q6+tZs2YN+fn5lJaWYjAYGDBgAHl5eeTk5HRqzoJQsrVmK3csvION1Ru5csiV3HrUrZiNQbpprBOklJQ3BxK/ALK9ydvM1eAXSAxAmi+Q5Nc30S/MyqdH5ezRpNV8oZNSUl5Ridvtplq3YTGb6B0bRqS1e/7977vvPmbNmoXJZOLII49k5syZWHvIRFUqWKCCRaiRUnLn3NXMWVHIcxcfydl57Y8Y2fnbaha88zqlWzeRlNGPcZddQ+bwtjt8S0pKyM/PZ/Xq1TQ2NhIWFkZubi55eXn06dMn5IeEtsfhcfDU8qeYtWEWQxOG8uS4J0mPTm9/wx5CSkmlW/M2Z/lqJc19JHUejTeH9SMnYs97F9atW8fAgYMora5HOhtolBZio6MC1uR0KFLBAhUsQpHTo3HZzF9YXVjLrFYjpPZF6jrrF//E/95/i7ryUjLzjmLcpVeTlNGvzfU1TWPr1q3k5+ezfv16PB4PCQkJLWlGYmPbP2Yo+37H99zz8z3oUufe0fdyetbpwS5SUEgpWbn6N8KT+xKhN2IwGEhOSsTcxt3cyr6pYIEKFqGqeYSUy6PzcRsjpPbF43az6uvP+OWjWTiaGhl6wgSOnXoZUQn77tx2OBysXbuW/Px8duzwjsnPzMwkLy+PwYMH95i7bFvzTxVyXvZ53H3M3YSbD52+mvZsLW/g/z5dy+WDTOT064tFukhISOgxTT+hRAULVLAIZRtK6jl/xiL6JUUw+4a2R0jti6OhgV/mzWbll58gDEZGnHEOR589GWs7HdvV1dUtw3CrqqowmUwMHjyYvLw8srKyetww3ANJFXKoaHJ5ePGHzbyycBtWk4FXz0khLSWJsLAw4uLigl28HkkFC1SwCHXfryvl2reWc1puCi9cvO8RUvtSW1bK/z54a3f6kMkXM3zCJIztdGpLKSksLCQ/P5+CggIcDgeRkZEtw3B79eq13+1DTWdThfQEUkq+Kijhwc/WUlTr4Pyj+nDXpIH8vqGAtLQ0kpOTD6ucYl1JBQtUsOgJXl64hUe+WM8tE3L48ykDDmofJVs2sfCd19i5dg1xqb05/uKryD5mTIcumB6Ph40bN5Kfn8+mTZvQdZ2UlJSWNCORkZEHVabuVmmv5J+L/sn/dv2PCX0nHFCqkFBXVu/gL7Pz+WlTBYNSonjw3FyOzoxnw4YNVFVV9ah/p1CkggUqWPQE/iOknr/4SM7qwAipfe1n28rlLHz3dSoLf6f3gMGMu+wa+gwc3P7GPo2NjRQUFJCfn09RURFCCLKzs8nLy2PgwIFtZi4NJS2pQlY8S2J44FOFdIfC6iYum/kLZfVO7jx1IJeNzsBkNOB2u3nxxRcZO3YsRx999CFXk+pOKliggkVP4T9CavYNY8jrwAipfdE1jYIfv+PnOe/SWF1FzqixHH/xlcSl9jmg/ZSVlbF69WpWr15NXV0dVqu1Zba/vn37hvTFKRRShXSFbRWNXPrKEhqcHt645hiO6ru7T2LBggX88MMPXHDBBQwbNiyIpfSaPn06r7zyClJKrrvuOm677bZgF6nDVLBABYuepKLByTkvLMKt6Xzyp+NIiencKCW3w8Hyz/7Lsk8+RPO4yTvldEZfsDt9SEfpus727dvJz89n7dq1uN1uYmNjW4bhJiQktL+TIGhwNfDA4gf4cnvopQrpiPUldVw2cylSSt6eNoohvXdnHa6urubFF19kwIABDBs2bK8LXXcrKCjgoosuYunSpVgsFiZNmsS///1vcnJyglqujlLBAhUsepr1JXVcMONnspIimX3DGMIsnf813FhTzeK577H6+68xW20cc85kjjrjnD3Sh3SUy+Vi3bp15Ofns3XrVgDS09NbZvsLC+vYEODu0hNThQCsLqzhiteWYjMZeefaUWQn79kfMWvWLDZv3sxNN91EcXFx0IPFnDlz+Prrr5k5cyYADz74IFarlTvvvDOo5eooFSxQwaInah4hdXpuKs9ffOQBj5Dal8rCnfz0/htsWf4LkQmJHDv1MoaMO8mXwvrA1dbWtqQZKS8vx2g0tsz2l52dHVIjc3pSqpCl26q45o1lxEWYee/a0aTH7zkcevPmzbzzzjuMHz+ecePG7XGh+/LLLykpKenS8qSkpHDaaaftd51169ZxzjnnsHjxYsLCwpgwYQIjR47k+eef79KyBEpQJz9SlIM1YXAv7p40iEe/XE92ciS3H+QIqdYS0tI59457KFxbwIJ3XuXrfz/Lr198zLhLryYz76gD3l9MTAzHHXccxx57LMXFxeTn57NmzRrWrl1LeHh4y2x/qampQe/fyIrN4t3T3+Wp5U/x5to3WV66PCRThSzcWM71by+nT2wY7147eq+mSI/Hw5dffkl8fDxjx44NUin3NnjwYO666y5OOeUUIiMjycvLO2RyknWGqlkoASel5K9zVvPhr4W8cMmRnDm8a2cdk7rOhiX/43/vv0ltWSkZw49k3KVXk5yZ1an9aprWMtvfhg0b0DSNpKSklmG4raf3DIZQTRXy9W8l3PzeSrKTI3lr2jEkRu7dTLho0SK+/fZbLrnkEgYM8P6IaOtXcbD9/e9/Jy0tjRtvvDHYRekQ1QyFChY9mdOjcekrv7BmVy1z/jCG4WmxXX4Mj9tN/jdfsOSjD3A0NjB03HiOvfDy/aYP6Si73d4y29/OnTsB9pjtL5hpKUItVci8lbv4y5x8hqfF8MbVxxATtncTWV1dHS+88AKZmZlccsklLctDJViUlZWRnJzM77//zsSJE1m8eHGPuaNcBQtUsOjpmkdIeXSdj2/q/AipfWlJH/LVpwgER51xDsecMxlreESX7L+ysrIlzUhNTQ1ms5ns7GwsFgtCCAwGA0KIPR6BXiaRfLPjG777/TuSIpK4Ovdq0qLS9lov0OWYtayQf35SwOh+icy8ciQR+0gp/uGHH7J27Vpuuukm4uPjW5aHSrA4/vjjqaysxGw288wzzzBhwoRgF6nDVLBABYtDwbriOib/u2tHSO1LXXkZ/5v1Nut++oGwqGhGX3AxeadMwmjqms5gKWXLbH/btm1D13WklEgp93i9v2U96e/vQO0vqDgcDsaNG8f48eP32CZUgkVPpoIFKlgcKr5bW8p1by/n9GGpvHDxkQHvMC7dupmF777G7wWriU1J5fiLryRn1LFB76gG2gweHQ00+1pW46hhxqoZ5Jflc3Svo5mWO41wU/h+t+3sMRdtquB/m8sZ1CuSSbkpCPZ/XuHh4Rx77LF73UWvgkXnqWCBChaHkv8s2MJjX67ntpNzuO3krhkhtT9SSravWsHCd1+nYucOUgcM4oTLph1Q+pCexD9VSFJ4Eo+PezwgqUKklDz65XpeXriVySPSePyC4Rg7MTxaBYvOC1Sw6Fk5nJVDxg3jsjj/qD48+90mPl9dHPDjCSHod+RILn/iOSbecAt15WV8cO8dfPL0I1QV7Qr48bubQRi4cuiVvH362xiFkau/uppXVr+Cpmtddgxdl/xzXgEvL9zKlWMyeKKTgUIJbSpYKEEhhODR84cxIiOOv8xZxZrC2m45rsFgZNj4iUx79mXGTr2U7atX8uZfb+T71/5NU21Nt5ShO+Um5jLnrDlMzJjIcyuf44Zvb6CsqazT+/VoOn+dk8+7v/zOH0/sz/+dPbTLbrhU9qTrOvVVlei63v7KAaSChRI0VpORly4fQUKElWvfWkZpnaPbjm222RhzwcVMm/4yw8afSv63X/Lqrdex5KNZuJ3dV47uEGmJ5PFxj/PA2AfIL89n8ieT+anwp4Pen9Oj8af3VvLRyl3ccepA7po0KCT6fw5FmsdDdfEuGqurcNmbglqWoAQLIcQUIcRvQghdCNGpdjSlZ0uMtDLzypHUOzxc99Zy7K6uaybpiIjYOE6+9kaufOpF0ofmsWjW27x22w0U/PAtehc22QSbEILzcs5j1pmzSAxP5Mbvb+SpZU/h1twHtB+7S+P6t1bw1W8l3HvmEG46KTtAJVY8LidVRYV4nE5ie6ViiwjuHB/BqlkUAOcDC4N0fCWEDE6NZvpFR7JmVy13zM0PynDShD7pnHvHP7nw/x4jKj6Rr/8znbfvupVtq1YcUsNbs2KzeO/097hw4IW8ufZNLv/ycnbW7ezQtg1OD1e9vpSFm8p54oLhXHNcvwCXNnhqamqYPHkygwYNYvDgwSxevJiqqipOOeUUcnJyOOWUU6iurm5Z/9FHHyU7O5uBAwfy9ddftyxfsWIFw4YNIzs7m1tuuaXl/9IzzzzDkCFDGD58OBMmTGiZTx7gzTffJDs7mwEDBvLBnDnE907DFhnJtm3bGDVqFDk5OVx44YW4XK7u+0Kg7aF73fUAfgRGdnT9ESNGSOXQNeOHzTLjrs/ks99uDGo5dF2X639eKF+5eZp8auoZcvaD/5Cl27YEtUyB8N327+SY98bIUe+Okp9v+Xy/61Y3OuXZL/xP9v/b5/LjVbsCVqa1a9cGbN8H4oorrpCvvPKKlFJKp9Mpq6ur5R133CEfffRRKaWUjz76qLzzzjullFL+9ttvcvjw4dLhcMitW7fKrKws6fF4pJRSHn300fLnn3+Wuq7LSZMmyS+++EJKKeX8+fNlY2OjlFLKGTNmyKlTp0oppaysrJSZGRly7YqlctPq1bJfv36yqqpKSinllClT5Pvvvy+llPKGG26QM2bMaLPsbX2HwHLZ2et1Z3fQqYN3IFgA1wPLgeV9+/Zt88tRDg26rsvbZ62UGXd9Jj/LLwp2caTH7ZIrPp8nX7jmIvnUhWfKL154WtaWlwa7WF1qV/0uednnl8ncN3LlPf+7Rza6Gvdap6zOIU/91wKZ8/cv5De/lQS0PKEQLGpra2VmZqbUdX2P5QMGDJBFRd7/l0VFRXLAgAFSSikfeeQR+cgjj7SsN3HiRPnzzz/LoqIiOXDgwJbl7733nrz++uv3Ot6vv/4qx44dK3Vdl6++9JK87KILZVVRodQ8Hnn99dfL9957T+q6LhMSEqTb7ZZSSvnzzz/LiRMntln+QAWLgKVSFEJ8B6S08dE/pJQfd3Q/UsqXgZfBe59FFxVPCUFCCB45bxjbKxr5y5xV9I0PZ1ha8JL1GU1mjjr9HIacMIGl8+bw65efsGHxTxx1+jmMOndKl6UPCabekb15fdLrzFg1g5lrZrKqfBVPjnuSgfEDASiqsXPZzF8ornXw2lVHc1xO5/NsddTGjQ9S37CuS/cZFTmYAQPu2e86W7duJSkpiauvvpr8/HxGjBjB9OnTKS0tJTU1FYDU1FTKyryjynbt2sXo0aNbtk9LS2PXrl2YzWbS0tL2Wt7aq6++yqRJk6gtLWHHtq1kZmYSm9IbIUTLNpWVlcTGxrZkv93XvgIpYH0WUsqTpZS5bTw6HCiUw4/NbOSly0eSEGHlureWd+sIqX2WKSKScZdezTXPvsTA0cex7JMPmXnLdfz6xcdongPrIA5FJoOJW466hVcmvkK9q55LPr+ED9Z/wPaKBqb8ZzHl9U7ennZMtwaKYPJ4PPz666/88Y9/ZOXKlURERPDYY4/tc33ZRp+WEGKfy/298847LFu2jGmXXIyjsQFLWDiW8Ig91uvovgJNJWlXQk5SlHeE1AX//pnr31rOrBvGYDMHf/Kh6MRkTvvTXzjqjHNZ+M5r/PDmK6z86jOOu/hKBowOjfQhnTEqdRRzz5rLPxf9k4d/eRiD/ROEZyrvX38SuX26v4bXXg0gUNLS0khLS2PUqFEATJ48mccee4xevXpRXFxMamoqxcXFJCcnt6zfnIkYoLCwkN69e5OWlkZhYeFey5t99913PPTQQ3z07tsYDRCTnEpWTg4//vjjHtuceOKJJCYmUlNTg8fjwWQy7bWv7hCsobPnCSEKgTHA50KIr9vbRjm8DE6N5tkLj2D1rlrumLs6pEYk9erXn8n/fIjz/3Y/JouFz559jPfvvYPq4p5/J3hCWAJ/HPQIhuqz0W2/EZfzAm7z1mAXq1ulpKSQnp7Ohg0bAPj+++8ZMmQIZ599Nm+++SbgHbF0zjnnAHD22WfzwQcf4HQ62bZtG5s2beKYY44hNTWVqKgolixZgpSSt956q2WblStXcv311/P6jBdJSkz0jniKiOTUU0/lm2++obq6murqar755htOPfVUhBCcdNJJzJ07d6/jd5vOdnp050ONhjr8vPjDJplx12dy+nfBHSG1L5rmkavnfy1fuOYiOf2KyXLtTz8Eu0idsnx7pcy97ys59tHv5deblspJcyfJvDfz5Mv5L0uP5gn48UOhg1tKKVeuXClHjBghhw0bJs855xxZVVUlKyoq5Pjx42V2drYcP368rKysbFn/oYcekllZWXLAgAEtI56klHLZsmVy6NChMisrS950000tneYnnnCCTExIkEOHDJF5w4fLs846q2WbV199Vfbv31/2799fvvbaay3Lt2zZIo8++mjZv39/OXnyZOlwONose6A6uFUiQSWkSSn5y+x8Plq5i39fehSnDUsNdpHaVFdRzufPPUnRhrXknjSR8Vdfj9kamPk6AmXR5gqufXM5qTE23rl2FL1jw2hwNfDA4gf4cvuXjEoZxaPHP0pSeFLAynCoJxKUUtJQWUFjbQ3W8AhievU66Hnj90UlElQOS0IIHjl/GEf2jeX22aso2NU9OaQOVHRiEhfe9yijzruQgh+/5d2//5mKnTva3zBEfLe2lKvfWEZGQjizbhhD79gwoI1UIZ92LlXI4UzXdWpKi2msrSE8JpbYlNQuDxSBpIKFEvJsZiMvXz6S+HAL1765nLIQGCHVFoPRyHEXXc4Ff38Ae30d7/79z6z54ZuQ6m9py6f5RfzhnRUMTonig+tHkxS15xSx/qlCEsISuPH7G3l6+dMHnCrkcKZ5PFQXFeJsbCQqMYnoxKQeNyBCBQulR/COkDqaOoeb695ajsMdunmbMocfyRVPPE/vAYP45j/P8cXzTwU9Cdy+zF62k1s/WMlRGXG8c+0oYsMt+1zXP1XIG7+9wRVfXsHO+o6lCjmcuZ1OqnbtxON2E5fSm4iY2GAX6aCoYKH0GEN6R/OvC48gvzD0Rki1FhEbxwX/eIBjL7ycDT//xNt330rp1s3BLtYeXl+0jTs/XM1xOUm8efUxRNnan27WZrLxz9H/5F8n/osd9TuY8ukUvtj6RTeUtmdyNjZSVeQdPhvfOw1rRM+9kVMFC6VHOXVoCnecOpBP84t4YX5oXXxbMxiMjD7/Qqbe9wgel4v37/krK7/6NCSC3Is/bOb+T9dy6tBevHLFiAOeC/3kjJOZe9ZccmJzuOunu7h30b00uUOz9hQsjbU1VJcUYTKbie+TjtlqbX+jEKaChdLj3Hhif847sg9Pf7uRL9cEfpa9zkobnMvljz9HxvAjmf/6S3zy9CM4GhqCUhYpJU98tZ4nv97AeUf24cVLjsJqOrhO1uZUIdcNu455m+dx0ecXsaFqQxeXuOeRUlJXUU59RTnWiAjieqdhNPX8+59VsFB6nOZZ9o7sG8ufZ+eH7Agpf+HRMZx7572ccPk0tv66jLfvvoWijeu7tQy6Lrn/07XM+HELl4zqy9NT8jAZO3cJ2FeqkFCoPXVGZ1KUf/jB+zTV1hARE8uWwiLy8vIOOEV5Tk4OOTk5LTcBArzwwgtkZ2cjhKCioqL7voxmnb1Rozsf6qY8xV9pnV2OeeQ7OfqR72RprT3Yxemwok3r5St/ukY+c/HZ8pd5c6SuaQE/pkfT5V9mr5IZd30mH/rst70yqnaFiqYK+Ydv/yBz38iVt86/VdY4ag54H6FyU97BpChvbKiXyxf9JDP6pss63w17B5OivF+/frKyslJWVVXtkaL8119/ldu2bZMZGRmyvLx8n2UP1E15qmah9FjJUTZeuXIkNU1urn97RUiPkPKXmj2Qyx6bTvbI0fz03ht89Pj9NNUFrnbk8ujc8sFK5q4o5LaTc/j76YMDMmwzISyBFye8yF9H/pUFOxcw5dMprCxb2eXHCbS6ujoWLlzItGnTALBYLMTGxvLxxx9z5ZVXAnDllVcyb948AD7++GOmTJ5MY0U5aamp5OQMoGDDBoqLi6mrq2PMmDEIIbjiiitatjnppJMIDw8HYPTo0S05pL7++mtOOeUU4uPjiYuL45RTTuGrr74C4MgjjyQzM7P7vohWen5DmnJYG9o7hn9deAR/eGcFd85dzfSLjugR49dtEZGcefvdrP7uS3548xXeuvNmzrj5r6QPHd6lx3G4NW5691e+X1/GP04fzHXjsrp0/60ZhIErh17JiF4juGPBHVz91dXcdMRNXJN7DcYDvAHtnk2FFDTYu7R8uZFhPJiTtt91DjRF+Y7t2xk2aADgHfHUNyPjgFOUn3baaYA33Xl6enq72wSDqlkoPd6kXO8IqU/yi3jxh9AeIeVPCEHeKadzyUNPYwkLZ86D/+TnOe912dzfjU4P17yxjPkbynj4vNyABwp/uYm5zDlrDhMzJvLcyue44dsbKG8q77bjd8aBpChvrKnB2diAwWjaY8TTgaQoX758OXfccQew73TnoUDVLJRDwo0n9mdTaT1PfbOR7ORIJuWGZg6ptiRnZnHZo//i+1f/zeK571G4dg2n3/xXIuMTDnqftXY3V7++lPzCWp6Zmsd5R+7/13QgNKcKGdN7DI/88giTP53MQ8c+xPFpx3do+/ZqAIHSkRTlRUVFJCYkUF9ZTnrfDKobGltGPB1IivKHH36YBQsWYPUFmbS0tDZTlIcCVbNQDglCCB67YDhHpMdy+6yeMULKn8UWxmk3/ZlJN95O8ZaNvHXXLWxbteKg9lXZ4OTil5ewZlctL15yVFACRbOemCqkvRTluqbx0gvPM3H8SUTExjH1kkuYNWvWAacov+GGG/jkk09a5sUA9pmiPCR0toe8Ox9qNJTSnj1GSNX1nBFS/ip2/i7f+OtN8qmpZ8gF77wmPb55lzuiuMYuJzz9oxz4zy/kjxvKAljKA2d32+WDix+UuW/kyvM/Pl8+uPhBOWPlDPnBug/kt9u/lStKVsjVBaulR/MEZLTWgdhXivKTTjpJ9svMlMeNGS13btvWsv6BpiifMGGCTE5Olnl5eTIvL69DKcqnT58u+/TpI41Go0xNTZXTpk1rs+wqRTkqRbnSMb8V1TL534sZ6EuMFwqz7B0ot8vJj2++wurvviI1ZyBn3noX0UnJ+91mZ1UTl878hapGF69eOZJRWQffjBVI3+34jn/n/5uypjJqnDV7fPbskGdJ6ZeCEAKTMGE0GDEZTN6H8D4bDcaW1yaDCaMwdku7vtvhoLq0CKlLYnulYvWNZgo1gUpRroKFckj6qqCYP7zzK+cc0ZtnL+wZI6TasmHxT3zz0vMIg+DUP95GztFj2lxvc1kDl838Bbtb461rjiEvPbZ7C3qQ3LqbGkcNlY5KKu2V2CptZORk4NE9aLqGR3rw6J6W95K2r1ctQaVVEGkr0BjEgbe+OxobqC0twWAyEZuSitkSuqk7AhUsVAe3ckialJvKXycO4KlvNjKgVxQ3nZQd7CIdlIFjjqdXVg6fPfs4nzz1MEdOOotxl12Dybw76d/aojouf/UXhIAPrh/N4NToIJb4wJgNZpLCk1omVFpXt47EsMQ215VSokltryDikb7A4nvd5GlC0zV0qbe5H4Mw7BFEWgca//cGYaCptob6ygrMNhuxvVIPidQdB+PwPGvlsHDTSdlsKmvgya830D8pkkm5KcEu0kGJ7ZXCxQ8+wU/vvcGKzz9m14a1nHnrncSl9mHl79Vc+dpSIq0m3rl2FFlJkcEubsA0N02ZDCastP/LXtM1NKntrpn4vW4OME7Nicfjfd2WSLsJq8uAZhW4InUcjtKgN4cFiwoWyiFLCMHjFwxnR2UTt89aRXr8GIb2jgl2sQ6K0WTmxCuuI33ocL6a8Sxv330bmWdfyd/yTSRGWXn32lGkxYVmG3qwGA1GjBixGPc9R0czXeotNRZN13BrLlyVdUiXGz3CjCdc4NE92F0ONKntM/dVoJvDgkn1WSiHvLJ6B+e8sAgBzPvTsSRH9ay5sVurqyjnvccfpvH3zexMGsad991N76SeGQRbC4U5uD1uNzUlRWhuN1GJSYRH7/nd7qs5rHXNpbmfpauawzpaa1F9FopykJKjbLxyxUim/GcxN7y9gvev65kjpJotKvbwjPlkTk1LInPXEr557G+cedtdJKZnBLtoPZ7L4aCmpAiQxKb2xhq2d23tQJvDdKnvDiL+QcUvoLTXHCaEoG9UXyItwWtm7Fn1IEU5SLl9YvjXhXms/L2Guz8M7Vn29ufDFYXc9N6vDE+P5+GH7+xx832HMkdDPdVFhQiDAUN4FJdefsVBpSj/+uuvW5avWLGCvOF5DBk4hLv/cjdRlijibHG8/Z+3mXD0BM449gymnT8Nc52ZQfGDGJIwhCWfLuGcMedwzuhzWDhvIb0iehFvi+faq65l4MCB5Obmcs011+B2d/ONjZ29UaM7H+qmPKWznvtuo8y46zP5wvxNwS7KAXtr8XaZcddn8tJXlshG5+4b9Rqqq+TsB/4un5p6hvxs+hPS2dQYxFJ2TjBSlOu6LuurKmXx5o2yovB36fG4DypFucPhkFu3bpVZWVnS4/FIKbs2Rfnnn38udV2Xuq7Liy66SM6YMaPN81EpyhWlC/xpfDZn5/Xmya838PVvJcEuToe9tGAL98wr4OTBvZh55UjCLbtbkHvCfN+hSkqduvIyGqoqsUVGEZ/ah8bGpgNOUX7RRRdhtVrp168f2dnZLF26tMtTlJ9++ukIIRBCcMwxx+yRd6o7qD4L5bAihOCJycPZUeUdITX3D2MZ0jt070uQUvKvbzfy3PzNnJXXm2em5mFuY3a75vm+0wYP5fPnnuT9e/7KCZdP44hTz+yxwznv//Q31hbVdek+h/SO5r6zhgKgaxo1pcW47HYi4+KJiItHCHHAKcp37drF6NGjW47RnFY8UCnK3W43b7/9NtOnT++ib6VjVM1COezYzEZeuXwEMWFmrn1zGeX1zmAXqU1SSh76fB3Pzd/MhSPTefbCI9oMFP5Cab7vUOZxu6kqKsTtcBCT3IvI+ISWoHogKcph32nF97Xc38GkKL/xxhsZN24cxx/fsey9XSUoNQshxJPAWYAL2AJcLaWsCUZZlMNTcrR3hNTk//zMDW8v570QGyGl6ZJ/zlvD+0t3cvWxmdxzxhAMho7VEJrn+17x+Tx+eu9N3r77Fs645U56DxgU4FJ3reYaQFdzOezUlBQDkrjU3lhajXjqSIry4uLilmyxaWlp7Ny5s2X7QKYov//++ykvL+ell17qqq+jw4JVs/gWyJVSDgc2An8LUjmUw1hunxj+NfUIfv29hr99tCZkRhK5NZ0/z17F+0t3cvP4bO49s+OBopkQgpFnnsdFDzyOEIJZ/3cXSz+ei9TbHvN/uLA31FNdtAthMBDfO32vQAHtpygHePPNN1vSjZ999tl88MEHAU9RPnPmTL7++mvef/99DIYgXLo720Pe2QdwHvBuR9ZVo6GUQJjuGyH14g/BHyHlcHvktW8ukxl3fSZn/LC5S/Zpb6iXnzz9iHxq6hly7iP3ysbami7ZbyAEajSU/4inysKd0uPZf9r3faUoHz9+vMzOzpbjx4+XlZWVLet3R4pyo9Eos7KyWra5//772yz7IZuiXAjxKTBLSvlOe+uqO7iVQJBScssHq/hsdREvXTaCiUODk0OqyeXhhrdX8NOmCh44ZyhXjMnssn1LKVvm+7ZFRgVkvu+uEIg7uL0jnsqx19cRFhVFdGIyIhi/zLtJoO7gDtg3JoT4TghR0MbjHL91/gF4gHf3s5/rhRDLhRDLy8t7xhy+Ss8ihODJycMZ3ieG22at6vIROB1R53Bz5WtLWbS5gicnD+/SQAGt5vu2hXX5fN+hStc0qouLsNfXERkXT3RSr0M6UARS0GoWQogrgT8AE6SUTR3ZRtUslEAqq3Nw9guLMBoE8246lqSo7pmzoLrRxRWvLWVdcR3TLzqSM4YHdv5wl8PO9zNnsPanH0gfMqzT8313pa6sWXjcLmqKi9E8bqKTkgmLCt0h0l2px9Us9kcIMQm4Czi7o4FCUQItOdrGzCtHUtno5A/vrMDhDvyv7rI6Bxe+vJgNpfW8fMWIgAcK8M33/ae/dMl836HKZbdTtasQXdeIS+1z2ASKQApWfewFIAr4VgixSgjxnyCVQ1H2kNsnhmemHsGKHdX8PcAjpAqrm5j60mIKq+28cfXRjB/UK2DHasvQEyZw2SPPEhEbx0eP3sfCd19H83i6tQyBYK+vo7p4FwaDkfg+aVjCwoJdpENCUIKFlDJbSpkupTzC9/hDMMqhKG05fVgqt588gI9W7uI/C7YG5BjbKhqZ+p/FVDa6eHvaKMb2b3t2uEBLSEvnkoefZvjJk1j2yYfM+r+7qCsvC0pZOktKSUNVJbVlpZhtNuL7pGEytz+XhdIxqqdHUdpwy4RszhyeyhNfr+fbtaVduu8NJfVM+c9inB6dD64fzYiMuC7d/4EyW6ycct2fOPO2u6gs3Mlbd93MpmWLg1qmAyV1nbryUhqqqwiLiiYutTcGY+jcZHkoUMFCUdoghOCpKXkM7xPDrR+sZF1x14yQWl1Yw4UvL8ZogFk3hNbMfQPHHM/lj00ntlcqnzz1MPNffwlPd6fBPgia5vGNeKonMj6B6KRkRCdnoaupqWHy5MldkqJ82LBhZGdnc8stt7Q0az7zzDMMGTKE4cOHM2HCBHbs2NGyzZtvvklOTg45OTktNwECTJs2jby8PIYPH87kyZNp6O40Lp29UaM7H+qmPKW7ldTa5aiHv5NjH/1eltc7OrWvX7ZWyqH3fiWPe/x7uaMidNOIu10uOf+Nl+VTU8+Qb911i6wqKuy2Yx/oTXlup1OW7dgmS7Zskk31dV1WjlBMUV5bW9tSvttvv72lLK2pFOWKEgS9fDmkKhud3PD2CpyegxshtXBjOVe89gvJ0VZm3zCGvgmhO1+2yWzmpCuv49w776GuvIy3776NdYsWBLtYe3HZm6gq2onUdeJ69yEsMqpL9ltXVxeSKcqjo70juqSU2O32bs8mrFKUK0o7hqXF8PSUI7jpvV/520dreHpK3gH9oX79Wwk3v7eS/smRvD3tGBIju+f+jc7qP2IUlz/+HJ8/9yRfPPckOwvyOemq6zFbu2kO8y/vhpI1bX6kaR5wu4kTApPF0vFmp5RhcNq+M8gCIZ2i/Oqrr+aLL75gyJAhPP300x075y6iahaK0gFnDE/ltpNz+OjXXby0sOMjpD5etYsb3/2VoX2i+eC60T0mUDSLTkziwvseZdR5U1nzw7e8+/c/U7FzR/sbBohEonncaG4XwmDAZLV2un+itVBOUf76669TVFTE4MGDmTVrVofPqSuomoWidNCtE3LYVNbA41+tp39SJKcM2f99Ee8v/Z2//3cNo/rFM/PKo4m09sw/N4PRyHEXXUHakGF8+cLTvPv3PzP+mhvIPfGUwDaFtKoBSF2ntrwMR0M9YVHRRCcldXmggNBOUQ5gNBq58MILefLJJ7n66qu7+vT3SdUsFKWDhBA8NTmPYX1iuO2Dlawv2fcIqZk/beVvH63hxAFJvHH1MT02UPjLHH4kVzzxPL0HDOKb/zzHF88/hcvePQkYNM1DVfEuHA31RCV0zYinfQnFFOVSSjZv9k6VK6Xk008/ZdCgbp6fpLM95N35UKOhlFBQUmuXxzz8bZsjpHRdb0l5/sd3lkunWwtSKQNH0zxy8YcfyKcvPEvOvOVaWbK1a1KpS9n2SB7/EU/2LhzxtD+hlqJc0zQ5duxYmZubK4cOHSovueSSPUZH+TtkU5QfCJVIUAkVqwtrmPrSYnJ7x/DudaOwmoxIKXnsy/W8tHArFxyVxuMXDMPUzjSoPVnhugI+f+5J7HW1XTbfd+skeE57E7UlxSAEsSm9sdi6qXO9BzukEgkqSk83PC2Wp6bksXxHNX//qABdl9zzcQEvLdzKFWMyeHLy8EM6UEDg5/tuqqulprgIg8lEfJ90FSiCrOc3pCpKkJw5vDebShuY/v0m1hbXsa64jj+c0J+7Jg3s9jHwwbLnfN9vdMl831J6czw11lRjCQ8nNjlFpe4IAYf2Tx9FCbBbJ+RwxvBU1hXX8deJAw6rQNFs93zfT3R6vm8pJbVlJTTWVBMeHUNcisrxFCpUn4WidJJb09lc1sDgVDVngqOxgW9fep6Nvywi84gRnHbTnwmP7lj+q8aaatatW0dachJRCYmEx8QedoG3K6g+C0UJUWajQQUKH1tEJGfefjcTpt3Izt9W8/adN7Pzt9Xtblexcwfv/fMvaB4PsSmpRMTGqUARYlSwUBSlSwkhOGKid75vcwfm+96+eiXv33MHmttNREwstojIbi6x0hEqWCiKEhDJmVlc9tizDD7uBBbPfY+5D/6ThqrKPdZZ/f1XfPTofUQnJXPJw09jNJuDVNo9hWKK8mY333wzkZFBCKidvVGjOx/qpjxF6ZkKfvxOPnv5+fLFay+RW1cul7qmyR/fflU+NfUMOfeRe6XDl677QFOUB0oopiiX0nuT32WXXSYjIiL2WXaVolxRlB6r9Xzfb//tNpZ/+hF5E8/gvDvvxRoeOinbQzVFuaZp3HHHHTzxxBPd+G3spu6zUBSlWzTP9/3jm6+wZv43nHjFdRx1+tn77Mh+fOnjrK9a36VlGBQ/iLuOuWu/64RqivIXXniBs88+u6UM3U0FC0VRuk3zfN8nXnkdZktopmtvTlH+/PPPM2rUKG699daApyhfsGDBfvdVVFTEnDlz9shI291UsFAUpdt1JFC0VwMIlFBMUb5y5Uo2b95MdnY2AE1NTWRnZ7dkou0Oqs9CURTFTyimKD/jjDMoKSlh+/btbN++nfDw8G4NFKBqFoqiKHt5/vnnufTSS3G5XGRlZfH666+j6zpTp07l1VdfpW/fvsyZMweAoUOHMnXqVIYMGYLJZOLFF1/E6EtR8u9//5urrroKu93Oaaed1tI3cccdd9DQ0MCUKVMA6Nu3L5988gnx8fHcc889HH300QDce++9xMfHB+Eb2JtK96EoSshoK1WFcmBUug9FURQlaFSwUBRFUdqlgoWiKIrSLhUsFEVRlHb1qNFQK1asqBBC7Gh/zTYlAhVdWZ4gU+cTug6lc4FuPJ9vv/12mKZpnkAeQ9M0k9FoDOgxulPr8ykpKTENGTJkTavVMjp7nB4VLKSUSQe7rRBieWdHA4QSdT6h61A6F+je88nPz9+em5sb0MBUUFAwODc3d10gj9GdWp+PpmmJgfj3Us1QiqIorVRUVBgnTZqU1a9fv6FZWVlDv/vuu4jS0lLj2LFjczIyMnLHjh2bU15e3jLf69/+9reUvn375mZmZuZ++OGHLTNh/fTTT+EDBgwY0rdv39yrrroqXfdNNft///d/vfr37z90wIABQ8aMGTNg48aNluZtnn/++YSMjIzcjIyM3Oeffz6hefkFF1yQ2adPn2GDBg0aMmjQoCE///xzWDd9HYAKFoqiKHu5/vrr0ydOnFi3bdu239auXbv2iCOOcNx3332pJ554Yv2OHTsKTjzxxPp77703BWDFihW2jz76KH7Dhg2/ffXVVxtvu+22vh6Pt1XoxhtvzJgxY8aO7du3F2zdutU2d+7caIARI0Y0rVq1at3GjRvXnnvuudW33357GkBpaanx8ccf77106dJ1y5cvX/f444/39g9KDz30UOH69evXrl+/fu3YsWPt3fmdHE7B4uVgF6CLqfMJXYfSucAhdj6JiYnl+/u8qqrK8Msvv0TddtttFQA2m00mJiZqX331VewNN9xQCXDDDTdUfvnll3EAc+fOjT3//POrwsLC5KBBg1wZGRnOH3/8MWLHjh3mhoYGw8knn9xoMBi49NJLK+fNmxcHcNZZZ9VHRUXpAMcdd1xDcXGxBWDevHkx48aNq+vVq5eWlJSkjRs3ru6jjz7a7yTm7Z1PV+lRfRadIaU8pP7Dq/MJXYfSuUDwzqfo7/9Id27a1KUTXVhzcpp6P/Lwzv2ts379emt8fLxnypQpmWvXrg0fPnx44yuvvLKzsrLSlJGR4QbIyMhwV1VVmQB27dplGT16dEPz9r1793bt3LnTYrFYZGpqqrt5eUZGhqu4uHivqQBfeumlpJNPPrnWty9zWlqaq/mzPn36uHbt2tWyzf3339/n0UcfTT3++OPrX3jhhcKwsDCZkpLSLYMPDqeahaIoSrs8Ho9Yt25d+E033VS+bt26teHh4fo999yTsq/195FWXHYkRfmMGTPi8/Pzw++///6S/ewLgGeeeWbX1q1bC/Lz89dVV1cb91emQDhsahaKovQs7dUAAiUzM9PVq1cv1/jx4xsBLrzwwurHHnssJSEhwbNjxw5zRkaGe8eOHeb4+HgPQFpammvnzp0tHdRFRUWWtLQ0d2Zmptu/JrFjxw5LSkpKS01j3rx5UU899VTqTz/9tCEsLEz69uVesGBBVPM6u3btspxwwgn14K3NAISFhclrrrmm8umnn+4V6O/C32FRsxBCTBJCbBBCbBZC3B3ksrwmhCgTQhT4LYsXQnwrhNjke47z++xvvnJvEEKc6rd8hBBije+z54Tv54cQwiqEmOVb/osQItNvmyt9x9gkhLiyC84lXQjxgxBinRDiNyHErT38fGxCiKVCiHzf+dzfk8/Hb79GIcRKIcRnPf188vPzh61Zs2ZIQUHBkIKCgsEAbrfbuG7dupzVq1fnrlu3Lsftdrd0CBcWFqasXr06d/Xq1bnV1dUto5Tq6+vD16xZM2T16tW527ZtS2/+Ra/runA6nX2Tk5PNn3766WCHw2H55ptvogcOHOiYMGGC4+mnnx6yevXq3Oeffz590qRJNQAXXHBBzUcffRRvt9vF+vXrLdu3b7edeOKJjRkZGe6IiAj9+++/j9B1nXfffTfhnHPOqQFYtGhR2M0335w5Y8YMZ2Vl5cA1a9YMraurizjzzDMbFixYkPDTTz/lLl68OGfBggXR5557bi3AsmXLeq9evTo3Pz8/d86cOcmDBw+2t3Uuft9/1/7bdHYS71B/AEZgC5AFWIB8YEgQyzMOOAoo8Fv2BHC37/XdwOO+10N85bUC/XznYfR9thQYAwjgS+A03/Ibgf/4Xl8EzPK9jge2+p7jfK/jOnkuqcBRvtdRwEZfmXvq+Qgg0vfaDPwCjO6p5+N3Xn8G3gM+C/X/b6tWrdoupVy+r8eqVaucLpdrlf+y7du3l+zcubNQSrl8586dhdu3by+WUi5vaGgoWLNmTZOmaSuamppW5+fnO3RdXy6lXF5QUNBYW1u7Ttf15evWrautrKzcKKVcXlRUtGPLli1lixYt+m3IkCGO7Oxsz4QJE6p37dq1+scff3SOGjWqvm/fvo6jjz5aLywsXN1chjvvvLMwLS3NkZGR4Zg1a9bG5uULFixYm52dbU9LS3NcfvnlZZqmLZdSLh8zZkxdfHy8PmDAANfAgQObTjrppBq3271y+/btJQ8//HBVenq6Iy0tzfPwww/XNJ/LyJEjtezsbHv//v0dp512mlZdXf1rW+eyZMmS0kD82wT9Yh7oh+8/+Nd+7/8G/C3IZcpkz2CxAUj1vU4FNrRVVuBr3/mkAuv9ll8MvOS/ju+1Ce+dt8J/Hd9nLwEXd/F5fQycciicDxAO/AqM6snnA6QB3wPj2R0sQvZ8DiZY5OfnOxwOR76UcrnD4cjPz893SF/gaA4iUsrl69atq62trV3nW8fevLysrGzrli1byv3XkVIu1zRt+a+//urWdX2PdaSUy7ds2VJeVla2dX9l3d/D7Xb/umrVKmdz8OrKc/nxxx/rA/Fvczj0WfQB/Ns+C/FeAEJJLyllMYCUslgI0Tx1Vh9gid96hb5lbt/r1subt9np25dHCFELJND299CHLuKr4h6J99d4jz0fIYQRWAFkAy9KKX8RQvTY8wGeBe7EW/Nr1pPPhw0bNuSAd8hoSkpKhcfjMVmtVjeA1Wp1ezweE4Db7bZERES0jFIym80ul8tlEUJIs9nc0ndgsVhcbrfb3LyN1Wp1ARgMBoxGo+bxeEwul8tsNptdrfa118imjnI4HFaTyeTZsmVLpsPhCA8LC2vMzMzc2RXnomlaczNcl/7bHA7BQrSxrKfM+LSvsu/vnA5mm04RQkQCHwK3SSnrhGjrUAddtm49HymlBhwhhIgF/iuEyN3P6iF9PkKIM4EyKeUKIcSJHdlkH2UIifMBGDRo0Hqr1ep2uVymjRs3DggLC3Mc4C4Opgxdfr2QUgq73R6enp7+e3R0dOO2bdvSd+3adaCjm9osl9/fX5f+2xwOHdyFQLrf+zSgKEhl2ZdSIUQqgO+5zLd8X2Uv9L1uvXyPbYQQJiAGqNrPvjpFCGHGGyjelVJ+1NPPp5mUsgb4EZjUg8/nWOBsIcR24ANgvBDinR58PjT/6rZYLJ6YmJiahoaGCJPJ5HE6nWYAp9NpNplMHtj967t5W7fbbbFYLG6LxeJurkkAuFwuS/Ovc7PZ7HI6nRYAXdfRNM1oMpk03zZ77asT5+Eym82u6OjoRoD4+Phqu90e3hXnYjAYNN/bLv23ORyCxTIgRwjRTwhhwdvR80mQy9TaJ0DziIQr8bb9Ny+/yDeqoR+QAyz1NSHUCyFG+0alXNFqm+Z9TQbmS2/D5NfARCFEnG/0y0TfsoPmO/arwDop5TOHwPkk+WoUCCHCgJOB9T31fKSUf5NSpkkpM/H+v58vpbysp56PpmkGj8djaH5dX18fHRYWZo+Ojq4pLy9PACgvL0+IiYmpAYiLi6upqamJ13Vd2O12i9PptEVFRTVarVa30WjU6+rqIqSUVFZWJsTGxtYAxMTE1FRUVCQAVFZWxkVGRtYLIYiNja2tr6+PdrvdRrfbbayvr4+OjY2tPdhzsVgsHrPZ7GpqarIC1NXVRVutVkdXnIvNZmvyHaZL/20O+WYoX1vdn/B+GUbgNSnlb8EqjxDifeBEIFEIUQjcBzwGzBZCTAN+B6YASCl/E0LMBtYCHuAmXzMJwB+BN4AwvKNTvvQtfxV4WwixGe+viIt8+6oSQjyIN3gCPCClrOrk6RwLXA6sEUKs8i37ew8+n1TgTV+/hQGYLaX8TAixuIeez770yH8fl8tl2rJlS7Zv/yIuLq4yPj6+LioqqnHz5s39V69enWg2m13Z2dlbACIiIhyxsbFVBQUFQwHS09N3NDfR9O3bd8f27dv7SSlFVFRUXVxcXC1AcnJyxZYtW/qtXr0612g0allZWVsAzGazlpKSUrRu3brBACkpKUVms1lro5gd1rdv39+3bt2aJaUUFovFmZWVtR2gs+dis9mac0Z16b+N8PWGK4qiBF1+fv72vLy8Q2kukG6Xn5+fmJeXl9nV+z0cmqEURVEOSCimKNd1nZtvvrlPZmZmblZW1tCHHnqoeRRbt1DBQlEUpZVQTFH+/PPPJxQWFpq3bNlSsHXr1t+uvvrqQDVTtkkFC0VRFD+hmqJ85syZyQ8++GCx0eit0PTp06dbp4Y95Du4FUXpmb5/a1161a6GLk1RHt8nsmnCFYN7ZIrynTt3Wt9+++24zz//PC4+Pt7z4osv/j5s2DBnZ7+TjlI1C0VRFD+hmqLc5XIJm80mCwoK1k2bNq38qquuyjywM+scVbNQDlu+eyoukVLO8L3vDTwnpZwcgGOdCwyXUj6wn3WeAr6QUs7v6uP3RO3VAAIlVFOU9+rVy3XJJZdUA1x++eU1f/rTnzID/FXsQdUslMNZLN7MnABIKYsCESh87gRmtLPO83izwCpB1LdvX09KSoorPz/fCtCcovzUU0+teemllxIAXnrppYQuSFGe8fHHH2/273s499xzaxcsWBBdXl5uLC8vN/qnKD/ttNNqvvzyyyiAL774IiojI6PbmqBA1SyUw9tjQH/fDYXfAi/izcyaK4S4CjgX742cucDTeFPcXw44gdN9Nzf1922XBDQB10kp1/sfRAgxAHBKKSuEEFHAamCAlNIthIj2vc+RUu4QQiQIIVKklCWBPnll355//vnfL7300iyXyyX69u3rfP/997drmsZ5553XPyMjI7F3796uefPmbQEYOXKk49xzz60aMGDAUKPRyDPPPLPDZPJeWmfMmLFj2rRp/RwOhzjppJPqpkyZUgtwxx13pDc1NRmnTJnSH7z9HPPnz9/cq1cv7Y477igaMWLEYIA777yzqFevXhrAAw88UDJ58uR+M2bM6BUeHq6/8sor27vzO1E35SmHLV+m3M+klLmt3/uCxT/xZtK1AZuBu6SU/xFC/AvYIaV8VgjxPfAHKeUmIcQo4FEp5fhWx7kayJVS/sX3/nXgYynlPCHE9cBAv89eAb6SUn4Y8C8gBKmb8jovUDflqZqFouzbD1LKerx5kWqBT33L1wDDfdl2xwJz/DourW3sJxUo93s/E2+z1DzgauA6v8/KgN5ddQKK0lVUsFCUffNvE9b93ut4/3YMQI2U8oh29mPHm/ETACnlIiFEphDiBLwz0RX4rWvzra8oIUV1cCuHs3r2nBTogEgp64BtQogp4M3CK4TIa2PVdXgnU/L3FvA+8Hqr5QOAAhQlxKhgoRy2pJSVwCIhRIEQ4smD3M2lwDQhRD7wG3BOG+ssBI4Uew6yfxfv/MfvNy8Q3rlBsoHlB1kWRQkY1QylHNaklJe0WpTrW/4G3pTczetl+r1u+UxKuQ3vBEn7O0aTEOI7YALwnW/xccBc3yRLzc70LevWNA6K0hEqWChK93gE39zvQojngdOA01utY8I7RFdRQo5qhlKUbiClLJVSfuJ7fbOUMltKubHVOnNa1TSUIAnFFOUjRowYOGjQoCGDBg0akpycPPzkk0/u301fB6CChaIoyl5CMUX5ihUrNqxfv37t+vXr1x555JGN5557bk13ficqWCiKovgJ1RTlzaqrqw2LFy+Oas4T1V1Un4WiKCHp638/m16xc0eXpihPTM9oOvWPt/XIFOXN3n333bixY8fWxcfH6wf7PRwMVbNQFEXxE6opypvNnj07/qKLLurWWfJA1SwURQlR7dUAAiVUU5QDlJSUGFevXh0xderUzYH9FvamahaKoih+QjVFOcBbb70VP378+Jrw8PBuzwCrahaKoiithGKKcoC5c+fG33nnncXd/oWgUpQrihJCVIryzgtUinLVDKUoiqK0SwULRVEUpV0qWCiKoijtUsFCURRFaZcKFoqiKEq7VLBQFEVR2qWChaIoSiuhmKL8448/jhoyZMjgQYMGDRkxYsTAgoICazd9HYAKFoqiKHsJxRTlt956a8Y777yzbf369WunTJlSdd9996V253eigoWiKIqfUE5RXlNTYwSora01+me07Q4q3YeiKCGpau7GdHdJY5emKDenRDTFTx7QI1OU/+c//9l+/vnn51itVj0yMlJbtmzZus5+HwdC1SwURVH8hGqK8meeeabXRx99tKm0tHT1JZdcUvHHP/4x/QBPrVNUzUJRlJDUXg0gUEIxRXlRUZFp3bp1Yc1luuKKK6onTZqUE/hvYzdVs1AURfETiinKk5KSPA0NDcbVq1dbAT777LPo7OxsR3d+L6pmoSiK0koopiifPn36jsmTJ/cXQhATE6O98cYb27rzO1EpyhVFCRkqRXnnqRTliqIoStCoYKEoiqK0SwULRVEUpV0qWCiKoijtUsFCURRFaZcKFoqiKEq7VLBQFEVpJRRTlH/yySdRQ4YMGZyTkzP0/PPPz3S7uzWPoAoWiqIorYVainJN07j++uv7ffDBB1s3bdr0W9++fV0vvPBCYnd+JypYKIqi+AnFFOWlpaUmi8WiDx8+3AkwadKkunnz5sV25/ei0n0oihKS5s2bl15WVtalKcqTk5Obzj333B6XojwlJcXj8XjEwoULw8eNG9c0a9asuOYA011UzUJRFMVPKKYoNxgMvPXWW1tvv/329GHDhg2OiorSjEbjXusGkqpZKIoSktqrAQRKKKYoBzj55JMbV6xYsQHgo48+it68ebMt0N+FP1WzUBRF8ROKKcoBdu3aZQKw2+3iySefTPnDH/5Q3p3fi6pZKIqitBKKKcofeOCBlG+//TZG13VxzTXXlJ199tn13fmdqBTliqKEDJWivPNUinJFURQlaFSwUBRFUdqlgoWiKIrSLhUsFEUJJbqu66L91ZS2+L47PRD7VsFCUZRQUlBeXh6jAsaB03VdlJeXxwAFgdi/GjqrKErI8Hg815aUlMwsKSnJRf2YPVA6UODxeK4NxM7V0FlFURSlXSpyK4qiKO1SwUJRFEVplwoWiqIoSrtUsFAURVHapYKFoiiK0q7/B674n/wEiAtvAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAERCAYAAAB/4wAeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABA/UlEQVR4nO2deZhcZZXwf6equ9OdnWxkJwlbErZsJCCKyAwxgA7DpqCO4hcHQZnBYfiQEQdFHgWdz5EgjAiKiDriKAoRExcWB4hASEgChBCyEekQ6NDZSa9V5/uj7q3cqr5Vde/t2rpyfs/TT1fd+973nqpO3nPP8p4jqophGIZhuMQqLYBhGIZRXZhiMAzDMDIwxWAYhmFkYIrBMAzDyMAUg2EYhpGBKQbDMAwjg5pQDCJyr4i0iMjLRZhrhog8IyJrReRFEfmo59xkEXlORDaIyC9EpKG39zMMw6g2akIxAPcBC4o01wHgk6p6nDPnbSIy1Dn3TeA7qno0sAtYWKR7GoZhVA01oRhU9Ulgp/eYiBwpIr8XkZUi8pSITA0412uqusF5/SbQAowUEQHOBH7lDP0x8PfF+gyGYRjVQl2lBSghdwNXqOoGEZkH/BephT0wIjIXaAA2AcOB3ara7ZxuBsYVUV7DMIyqoCYVg4gMBN4D/DL1oA9AP+fcBcDXfC7bpqof9MwxBvgJ8ClVTYpnIg9WT8QwjJqjJhUDKRfZblWdkX1CVX8N/DrfxSIyGPgd8GVVfdY5/A4wVETqHKthPPBmUaU2DMOoAmoixpCNqu4FtojIxQCS4qQg1zqZRr8B7lfVX3rmVOAJ4CLn0KeAh4squGEYRhUgtVBdVUR+DpwBjADeBr4CPA58DxgD1AMPqKqfCyl7rk8APwLWeg5fpqqrRWQK8AAwDFgFfEJVO4r4UQzDMCpOTSgGwzAMo3jUpCvJMAzDiE6fDz6PGDFCJ02aVGkxDMMw+hQrV658R1VH+p3r84ph0qRJrFixotJiGIZh9ClEZGuuc+ZKMgzDMDIwxWAYhmFkYIrBMAzDyKDPxxj86Orqorm5mfb29kqLkpfGxkbGjx9PfX19pUUxDMNIU5OKobm5mUGDBjFp0iT8SxxVHlWltbWV5uZmJk+eXGlxDMMw0tSkK6m9vZ3hw4dXrVIAEBGGDx9e9VaNYRiHHjWpGICqVgoufUFGwzAOPWrSlWQYhlGLqCr//ep/s6t9FwCzRs3iPePeU/T7lN1iEJG4iKwSkUd8zomI3C4iG51+y7PKLZ+X97zH/wu/7LLL+NWvfuV7zjAMo1Rs27+NW5ffyvdf/D53v3g3y99aXpL7VMJiuBpYBwz2OXc2cLTzM49UddR55RMtk7/85S+VurVhGEYPEpoA4Jb33cKHpnyoZPcpq8UgIuOBc4Ef5BhyHqk+COo0yBnqdFKrCAMHDgRS5ttVV13F9OnTOffcc2lpaamUSIZhHMIkNQmAUNr4ZLldSbcB1wHJHOfHAW943vv2VRaRy0VkhYis2LFjR9GFzOY3v/kN69ev56WXXuKee+4xS8IwjIqgTjfhmlEMIvIhoEVVV+Yb5nOsR8MIVb1bVeeo6pyRI32LAxaVJ598kksvvZR4PM7YsWM588wzS35PwzCMHjirYUxKu3SX02I4Dfg7EXmdVBe0M0Xkp1ljmoEJnvdV01fZUksNw6g0riupxAZD+RSDqv6bqo5X1UnAJcDjqvqJrGGLgU862UmnAHtUdXu5ZMzF6aefzgMPPEAikWD79u088cQTlRbJMIxDkHK5kiq+j0FErgBQ1buAJcA5wEbgAPDpCoqW5vzzz+fxxx/nhBNO4JhjjuH9739/pUUyDOMQxFUMpXYlVUQxqOqfgT87r+/yHFfg85WQyY/9+/cDKTfSHXfcUWFpDMM41EktkTUUfDYMwzB6R81lJRmGYRi9I20xlDgZxhSDYRhGHyFJbW5wMwzDMKLi7GMwi8EwDMMALMZgGIZhZGExhj7OG2+8wQc+8AGmTZvGcccdx6JFiyotkmEYfZxyxRgqvsGtVqmrq+Pb3/42s2bNYt++fcyePZuzzjqL6dOnV1o0wzD6KGYx9HHGjBnDrFmpPkODBg1i2rRpbNu2rcJSGYZRC8RKvHTXvMVw02/X8sqbe4s65/Sxg/nKh48LPP71119n1apVzJtXsZ5DhmHUADVXRO9QZf/+/Vx44YXcdtttDB7s17TOMAwjGIdMEb1SE+bJvth0dXVx4YUX8vGPf5wLLrigYnIYhlEbuDGGWurHcEihqixcuJBp06ZxzTXXVFocwzBqANvH0MdZtmwZP/nJT3j88ceZMWMGM2bMYMmSJZUWyzCMPky5spJq3pVUKd773vem/4iGYRjFwCwGwzAMIwM3K6lm9jGISKOILBeRNSKyVkRu8hlzhojsEZHVzs+N5ZLPMAyj2qnFrKQO4ExV3S8i9cDTIrJUVZ/NGveUqn6ojHIZhmH0DcpUXbVsisFp27nfeVvv/JgT3jAMIyA1GWMQkbiIrAZagD+p6nM+w0513E1LRcR3E4KIXC4iK0RkxY4dO0opsmEYRtVQczEGAFVNqOoMYDwwV0SOzxryAnCEqp4EfBd4KMc8d6vqHFWdM3LkyFKKbBiGUTXUpMXgoqq7gT8DC7KO71XV/c7rJUC9iIwou4BFoL29nblz53LSSSdx3HHH8ZWvfKXSIhmGUSPUzM5nERkpIkOd103A3wKvZo0ZLY6NJCJzHflayyVjMenXrx+PP/44a9asYfXq1fz+97/n2Wez4+yGYRjBSbuSaigraQzwYxGJk1rw/0dVHxGRKwBU9S7gIuBKEekG2oBLtI/uEhMRBg4cCKRqJnV1dZXcL2gYRm3zSusrqRclXkrKmZX0IjDT5/hdntd3AHcU9cZLr4e3XirqlIw+Ac6+teCwRCLB7Nmz2bhxI5///Oet7LZhGL3i9T2vAzC6/+iS3sd2PpeQeDzO6tWraW5uZvny5bz88suVFskwjL6MwMRBExneNLykt6n9WkkBnuxLzdChQznjjDP4/e9/z/HHZydiGYZhBERLn6oKZjGUjB07drB7924A2traePTRR5k6dWplhTIMo0+jaMkDz3AoWAwVYvv27XzqU58ikUiQTCb5yEc+woc+ZJU+DMOIjqJlsRhMMZSIE088kVWrVlVaDMMwaoikJstiMZgryTAMow9hisEwDMNIo1oeV5IpBsMwjD5CUpOmGAzDMIyDlCsryRSDYRhGH0HRkhfQA1MMhmEYfQZVsxhqgkQiwcyZM20Pg2EYvUbL1PTSFEOJWbRoEdOmTau0GIZh1ACq5krq8zQ3N/O73/2Oz3zmM5UWxTCMGiBJeTa41fzO528u/yav7ny18MAQTB02lS/O/WLBcV/4whf41re+xb59+4p6f8MwDlFqrYieiDSKyHIRWSMia0XkJp8xIiK3i8hGEXlRRGaVS75i88gjjzBq1Chmz55daVEMw6gRarGIXgdwpqruF5F64GkRWaqq3n6XZwNHOz/zgO85vyMT5Mm+FCxbtozFixezZMkS2tvb2bt3L5/4xCf46U9/WhF5DMPo+9TczmdNsd95W+/8ZIfYzwPud8Y+CwwVkTHlkrGY3HLLLTQ3N/P666/zwAMPcOaZZ5pSMAyjV5QrxlDW4LOIxEVkNdAC/ElVn8saMg54w/O+2TmWPc/lIrJCRFbs2LGjZPIahmFUFbUWYwBQ1YSqzgDGA3NFJLudmd8n7pG4q6p3q+ocVZ0zcuTIEkhaXM444wweeeSRSothGEYfp6ZLYqjqbuDPwIKsU83ABM/78cCb5ZHKMAyjuqm5InoiMlJEhjqvm4C/BbLzSBcDn3Syk04B9qjq9nLJaBiGUc3UYlbSGODHIhInpZD+R1UfEZErAFT1LmAJcA6wETgAfLqM8hmGYVQ1qko8Fi/5fcqmGFT1RWCmz/G7PK8V+Hy5ZDIMw+hL1HSMwTAMwwiPqvqn6BQZUwyGYRh9hFqMMRxyTJo0iUGDBhGPx6mrq2PFihWVFskwjD5MuaqrFlQMIjIx4Fy7VXVvL+WpOZ544glGjBhRaTEMw6gBqsli+DGpTWb5pFHgPuD+IshkGIZh+FCuGENBxaCqH8g+JiKjVfWt0ohUXN76xjfoWFfcstv9pk1l9Je+VHCciDB//nxEhM9+9rNcfvnlRZXDMIxDC0WJlSE0HDXG8EngW8UUpBZZtmwZY8eOpaWlhbPOOoupU6dy+umnV1oswzD6KOXa+RxVMZwnIgdIFcJbX0yBik2QJ/tSMXbsWABGjRrF+eefz/Lly00xGEYNsKt9F3euvpOOREf6mCBcdMxFnDjyxNDzPfjag8w8fCZThkzJO27dznWM6j8q9PxhiaoYLiC1We18ETlKVa13ZRbvvvsuyWSSQYMG8e677/LHP/6RG2+8sdJiGYZRBJ5/63l+sf4XDG8cTn28HoCWAy3EY/FIiuGrz3yV+lg9L/zDCwXH7unYE3r+sERSDKr6NvB758fw4e233+b8888HoLu7m4997GMsWJBdM9AwjL5IUpMA/PCDP+TIoUcC8De//Jv08Sh0JbsKjqmTOk4efXLkewQlkmIQkTuBAap6mYjMV9U/FlmuPs+UKVNYs2ZNpcUwDKMEJDQBZPZGiEmMRDJR0vsmqe7qqp3AZuf1mUWSxTAMo0/gWgZxOVjQLi7xXlkMQe9bjg1uUe9wABji9G4OugHOMAyjJnAVgHeRjkmMJKVTDKkao1R1uupOoA24E1hWPHEMwzCqn5yKIVk6xeDes+pcSSIyVER+BFzoHLofmFN0qQzDMKoYP1dSqS0Gd+6qcyU5LTlvBW4CngOOBn4d5FoRmSAiT4jIOhFZKyJX+4w5Q0T2iMhq58fyOw3DqDrSwWdPfYoYsZLGGNKupGoooufDQmCLqv4BWBnium7gX1X1BREZBKwUkT+p6itZ455S1Q9FkMswDKMsuIu0t5taLFbarKS0K6lKG/XsAq4QkdtE5NMi0qMrmx+qul1VX3Be7wPWAeMi3L/PsHv3bi666CKmTp3KtGnTeOaZZyotkmEYRcC1GLxP73GJR3IluUqmEH5xjVIR2mJQ1VtE5DHgNWAGcDqwKswcIjKJ1M7p53xOnyoia4A3gWtVda3P9ZcDlwNMnFi9SVFXX301CxYs4Fe/+hWdnZ0cOHCg0iIZhlEElJ4ZQoJEciW5cwW+ZzUqBhH5GhAHVgOrVfXPIa8fCDwIfMGnf8MLwBGqul9EzgEeIhXHyEBV7wbuBpgzZ06wb7XM7N27lyeffJL77rsPgIaGBhoaGiorlGEYRcF1GcViWRZDCWMM5XQlRbEYbhSRw0k98V8oIkeq6j8GudbZ9/Ag8DNV7RG09ioKVV0iIv8lIiNU9Z2wcro89T+v8c4b+6Ne7suICQN530eOyTtm8+bNjBw5kk9/+tOsWbOG2bNns2jRIgYMGFBUWQzDKD9+FkNMogWfw7qSqi5d1cNngVWqemsIpSDAD4F1qvqfOcaMdsYhInMd+VojylhRuru7eeGFF7jyyitZtWoVAwYM4NZbb620WIZhFAG/GENMYunjYQjsSqryrCSAe4ErRWQAqaf/1QGuOQ34B+AlEXHHfwln57Sq3gVc5MzbTWoD3SUaVJ3moNCTfakYP34848ePZ968eQBcdNFFphgMo0ZI72OIZe5j6OVylRdXgVSlK8nhn0nVS6oDbicVgM6Lqj5NgaZ0qnoHcEdEmaqK0aNHM2HCBNavX8+xxx7LY489xvTp0ystlmHkJalJtu3blpFdE5c44waOK4sLo6/QnewGerqS2rvbQ88V1GKo6qwkh02kgsIPq+q/FFGemuK73/0uH//4x+ns7GTKlCn86Ec/qrRIhpGX+9bex3dWfqfH8ZvecxMXHH1BBSSqTl7b9RqQuUh3JjrZuHtj+MkCGhlVnZXksBZ4A1goIv+hqqUvEN4HmTFjBitWrKi0GIYRmNa2Vupj9XzttK8BqafUG56+gR0HdlRYsupiWOMwINOV1FjXmD5eCsoZfI6qGI4BdpBKGd1VPHEMw6gkSU3SEG/gQ1NSxQdUlRuevoFu7a6wZNVFQhMMbxyecWxU/1Fs278t9FyhXUllqK4a9Q5TSW1quxZno5lhGH2f7Hr/IkKd1KV96kYKVc0ooAdlCD6XMSsp6h2GAl8ErgPCR1sMw6hKEprosfDUxUwxZJPQRA+XTqnTVd2EgGp2JX0NmKqq60WktC2LDMMoG35PwvWx+kD9iA8lkprs8T1F3flcjbWSAt1BROIisl1EPgOgqs2q+qjz+vpSCmgYRvkwiyEYSe3ZeznqzueguAqkaqqrqmoCeBk4srTiGIZRSRTtEdw0xdATP4shckmMKtzHEOYO/YHrRGSFiCx2fh4ulWB9nfXr1zNjxoz0z+DBg7ntttsqLZZh5CWRTGQUhoOUT/vBDQ/yTlvkkmU1R3aQHqJ3cAvsSipjB7cwMYZTnd+znB8IvDXj0OPYY49l9erVACQSCcaNG8f5559fWaEMowB+FsMpY05h8abFNO9rZkTTiApJVl3kVAwl7PmcdiVVWRG9yT4/U0ohVK3x2GOPceSRR3LEEUdUWhTDyItfjOG8I88DsAC0h1yKIUpWUph7Qnn2MQS2GFR1aykFKRVP3Hc3LVs3F3XOUUdM4QOXBd++8cADD3DppZcWVQbDKAXJZDJjNy+kYgxgisGLn2KISzxwvMBLX48xGBHo7Oxk8eLFXHzxxZUWxTAKkiTZI+vFVQwWgD6IX/BZREra8zldXbWK9zH0GcI82ZeCpUuXMmvWLA4//PCKymEYQfBb8Ewx9MRvg1vZej5XY0kMEflwKQSpVX7+85+bG8noM/jl59fH6gFTDF6S+FgMpe75XKXBZ5evF12KGuXAgQP86U9/4oILrFyx0TdIaCKnxWAxhoMkkz0VaDyW2vlcqnpJ1V52O5K6EpEJwP3AaCAJ3K2qi7LGCLAIOAc4AFymqi9EuV810L9/f1pb+2RnUuMQZVXLKsYOGJtxzFUMS7cs5dwp51ZCrJwkNcm3nv8WLQdaABjSbwhfmveltJUThH2d+/jm8m9yoPtA+lhDvIFrZl/DqP6jeoy/c/WdrN+1niMGZ2YZugv2tf97LQtPWMj04bkbczXva+b2VbfTnezOULjXP3U9X5j1BUYPGN3jmme3PwtU0c7nLKKqw27gX1V1GnAK8HkRyf7mzibVAOhoUlVbvxfxXoZhRKAx3sjezr0Zxw7vn4qPbX93eyVEyktrWys/W/czVres5sUdL/Kr137FX/f+NdQcr+58lYc3Pcwrra+wZc8W1u9cz+82/46Vb6/0Hf/9Nd9HVTl9fGbjytmjZnPMYcfwx61/5NGtj+a95zPbn2HplqVs2LWB5n3N6eO/2/w7nn/red9rNuzaAMBRQ48K8/EiUbbgs6puB7Y7r/eJyDpgHPCKZ9h5wP1On+dnRWSoiIxxrjUMo8QoytzRczOONcQbmH/EfDbs3lAhqXLj7hv43IzPMbhhMP/6v/8aei+BGxe4+bSbOXn0yby+53U+/NCHfeMFqoqiXDr1Uj5zwmcyzs0dM5cH/+5BZv5kZsG4getuuveD9zKy/0ggZUWc/euzc8Ypkppk0uBJjBk4JtTni0JF0lVFZBIwE3gu69Q4Up3hXJqdY9nXX+6U5lixY4d1ljKMYqGqvj7shngDnYnOCkiUH3eBjUs8HRsJ6+NPp4E6Lhr38/sqhgApozEK10zyCyTnu697vFx9t6Mohrd7c0MRGQg8CHxBVfdmn/a5pMdfWVXvVtU5qjpn5MiRvRHHMAwPfjufwSm9nai+4LNrHYhIetEMazFkN8Bxf/vNE2STWZCGPdnKqNB93XtnJwaUitCKQVXPinozEaknpRR+pqq/9hnSDEzwvB8PvBn1foZhhCOfxVCNWUlFsRiynt7zzROki5pI4bRVv/7N7n3zWQzlyEiCMrqSnIyjHwLrVPU/cwxbDHxSUpwC7LH4gmGUj3wWQ2ey+lxJRbEYsp7e882T7qKWJzMoSJVVP4vBve8hpRiA04B/AM4UkdXOzzkicoWIXOGMWQJsBjYC9wCfK6N8Rec73/kOxx13HMcffzyXXnop7e3WBdWobnK5K6o1xuAuol6LIewms+y4Qb55ArmSCN772TtPIflzKe1SECkrSUSucZ/6ReRYVV1f6BpVfZoCeyCcbKTPR5Gp2ti2bRu33347r7zyCk1NTXzkIx/hgQce4LLLLqu0aIaRk1xPpa4rSVXLFgANgnehLhS8LTRHtsWQKyvJvV8uwriSvBQMPvvsti4VoRSDiAwFvgNMFZF24EVgIfDp4ovW9+nu7qatrY36+noOHDjA2LFjC19kGBUkp2KINQCp3c8N8YZyi5UT193jVQxRS1+7iiGIxZB3ngCKIVJWks9u61IRSjGo6m7g0yJyLvAWMB/wCyJXDbt/u4nON98t6pwNYwcw9MP5u5yOGzeOa6+9lokTJ9LU1MT8+fOZP39+UeUwjGKTz2IA6Ex0VpVi8JaJcOXubfA53wIdpItajFjhfQx5spKqwWKI6rB6P6m01VOAyFlKtcyuXbt4+OGH2bJlC2+++SbvvvsuP/3pTystlmHkJZdicEtMVFsAOm0xEN1iyI4x5N3HUCRXkot3nvR9cwSuk9qzJHqpiLrzeSjwReA6Uq6kqqXQk32pePTRR5k8eTLuPosLLriAv/zlL3ziE5+oiDyGEYRcAU7XSni3612GNQ4rt1g5cVtpxmPRg8/ZMYa8FoMGzEoKmq4awmJIJBM9miiViqgWw9eAh52gc+manPZhJk6cyLPPPsuBAwdQVR577DGmTZtWabEMIyduuQc/d0W/eD8AntueXaygsnhdO1GDzzl3PvssbYGzkgK6kjKuKyC/4r/HpBREshhUtZnUZjRU9fqiSlQjzJs3j4suuohZs2ZRV1fHzJkzufzyyjYNMox8+G26cjl17KlA9fVk8DaviaoY3DU6e+dzvpIYvXUl+bmk3AY8uVxhCU2UpUkPRE9XvRMYoKqXich8Vf1jkeWqCW666SZuuummSothGIFwn5D9LIZ0jKHK9jKkFUOsF+mqWZZBOlbh06Yzn/L0Xl9QMfjUXCpoMagSi1WxYgA6OVgz6UzAFINhVBBVZcPuDT3KVoxoHMHhA4K1lc3nJnFjDOUsi9HW3cbmPZs5auhRaVeWly17trB5z2Yg02LYuncra1vXMqppVLpyaT5ylcToSHT0GBukvWY+a0JVefvA2wfv6bPzeeverb7Xvrn/TYb2G5pz7mISVTEcAIY4tY8mFlEewzAi8OhfH+WaP1/T43hTXRPLLllGfbxw45oXd7wI+C/+7j6GcmYlff3Zr/Pwpoe5+JiLufHUGzPOrW1dyyWPXJJ+P6B+AAPqBgDw7ZXfBmBQwyD+culfCt4n7R5yFnt3sd6yZ0vPsUGykvK0+PzvV/+bW5ffytmTz06N9bE8lm5Zyj/N/CcmDDpYNm7T7k20trf6KqtSEFUx7ATagDuBZcUTxzCMKOzp2AOkegq4T5WP/fUxHtr4EF3JrkCK4d2u1H6feaPn9TjnZsO0d5evrMueztRncj9bxjnn2BdmfYETR57I9OHTERF+vODH7O3cy5LNS1j6+tJAO7WzLQYRoS5Wx+B+g3uMTddKiuhKemrbU0Dq6R96Zjf988x/5vZVt7Ovc5/v5z3/6PPzfpZiEXXn87HAT0m16qzqdFXDOBRwF6LTxp6Wdp+4Lomgef3uuKa6Jt/zMYnxSusrvudKgfuZ8pW/nn34bGaMmpE+PuvwWQCs27kOSFkDhXL//TabjWoa5RtoL2XZbYBjDjsmdT7revc7GNlUnjYDoXc+i8itwCTgHeBEqnzns2EcCvgtWGHz+gu5SQY3DGZA/YDeiBkKdzHMFwTOtRPYdQsFqUjqV56iLlbn61Lziw1kIyIFq6tmZ0J5r3Xl9rtvNaerLgS2qOofAP+mqIZhlBV3EfUulIXKOOeaI9fiM3nIZPZ37e+NmKFwF8NuzfPkniNLJ0x5DL+n97pYXXSLIUAHt1zZTbn2UBT62xSbKHfZBVwhIreJyKdFZGaxhaoVFi1axPHHH89xxx3HbbfdVmlxjBrG76k3qsWQ6ym8f31/3u0sbt2xfOSzGLylMPxIK8UA+2/9FulcFkOQWkkiEtiVlI37ebKvL/S3KTZROrjdAvwj8FVgC3B6kWWqCV5++WXuueceli9fzpo1a3jkkUfYsKH6mqkbtYGfxRA2r9/b9MaPgfUDK2Ix+MUYCrlWouxpCGIx+Clgv3sH2cfg547KZeUV+tsUm9CKQUS+BpxHqnjeNlVdFPC6e0WkRUReznH+DBHZ42nic6PfuL7CunXrOOWUU+jfvz91dXW8//3v5ze/+U2lxTJqFD8XR9jFsZDffmD9wHTmUlge2/oYc382l79sK5w+6uIuhn4LdCHXSq4nbz/8NpsVdCUV2MdQsINbjhaquf5mhf42xSZ0jEFVb3QW7RhwoYgcqar/GODS+4A7SGUy5eIpVf1QWJnysXTpUt56661iTsno0aM5++yz8445/vjjueGGG2htbaWpqYklS5YwZ86cosphGC7FCD4X2tU7oH5AZIvh2ievpTvZzece+xwTB0/kwb97ML2bupA8L7S8QHt3O411jT3O5Voow8RX/Ara9Yv347ntz3Hdk9fxrdO/1WNsoX0MhRRSrkqp6dhIlqspyI7rYhI1knEvMA0YDvxXkAtU9UlS+x8OCaZNm8YXv/hFzjrrLBYsWMBJJ51EXV3UbSOGkR+/hTJsD+RCi+1hjYfR1t0WaS+Du1DOGT2HLXu2sL+zsILxLuo7DuzwPVfQlRQgxuDnHrrypCuZMmQKL7+T6eDwsy787h1IGftM4SqLPmcxOPwzqbIYdcAiihdnOFVE1gBvAteq6lq/QSJyOXA5pKqY5qPQk30pWbhwIQsXprZ5fOlLX2L8+PEVk8WobfyeKN1FJGjjmkKL7fDG4QC0trcybuC4wLKpKglNcMVJVzCyaSTPbX8uUDG+fItrIVndBTZM0x6ve+jk0Sdz0qiTeHrb0/737aUrKUnSd46crqQAQe9iEvUum4BGUqW3i6UUXgCOUNWTgO8CD+UaqKp3q+ocVZ3j9juoRlpaWgD461//yq9//WsuvfTSCktk1CrFsBgK+e2HNzmKoa01lGxudk9DrIG6WOpZtOSKIYorKcsKaIg19CgaGLRRTyGFlEgmfK2OXGm2hbKwik1Ui2Et8AawUET+Q1VP7q0gqrrX83qJiPyXiIxQ1Xd6O3eluPDCC2ltbaW+vp4777yTww47rNIiGTWKX6mGdIwhYMuUQoueqxjeaQv3X9JdXBvi0RVD9mcI7EoKoBhypY42xBt61CYKVF014D6GvFlJWZ83/bep8uqqR5Laz3C387vXiMho4G1VVRGZS8qaCfdoUmU89dRTlRbBOERIJBM9/M/pRSZZnA1uXldSGNzFtV+8Xzrg3KWFq7R6LZ1sRRI4K6lAwxzvmOy5+sX70ZXIlDNQz+cAJTG6kl3+FgP+Cq2vWAxvqOrjIjIGaAlygYj8HDgDGCEizcBXgHoAVb0LuAi4UkS6SRXou0TDdvU2jEMUv+5eoS0G8m+iSiuGkK4k12LoF++XthiyF1w/vItjtmIoZlZSrjIXDfEGurWb7mR3Wu70k3ueBTpISYyEJvJnJVV4g1tUxbBARF4jVV11K6lgdF5UNa+DXVXvIJXOahhGSPx6Ned6+sw5RzL/Jqr6eD1D+g0J7UpyLYb6eD114riSfMpcZBNEMeSSNZQrKcemNbcHRWeiM60Y0vPlyRrNV3bbJZlM+svuHKr0BreoimEo8EXgOuAzRZOmiAQpt1tpzCAyisUv1/+yR5DZXRxvfvZmBtYP5MwJZ/LRqR9Nn3+n7R1ufuZm2hOp9FO3FHS+p9LhjcP57abf8uVTvhxMrtd+yeKNi4FMi+Ebz36DgQ0D03JeNfMqjht+XPo6VU034QFY/tZyjhuROr+6ZTW/WP+LvLLmevLuSHTw70//e7qkN8D2d7f7zuE2B1q3cx2zD59NUpPc/OzNGfPnuvfyt5an3HtOufLORCf/vuzfWbYt1aVgZ8dOGuONPa91lXmWxbHy7ZV5P2+xieqw+hqpjKT1QLCUhzLS2NhIa2trVS+8qkprayuNjT3/cRhGWPZ37e/R5Wza8GnMG5PqrfDSOy/x4IYHM86/0voKj7/xOC0HWtjftZ/B/QZzzuRzcpbdhtRieaD7QOD/Ww9teIhNuzcxb/Q8jh9+PFOHTeWUMaeApGTe37mfp7c9zdPNmWmh2RbCtv3b0q+fbH6Sjbs3Mv+I+Qxu6NkzATzpqlkxhq17t7L09aVs278tdf+u/QxqGMRZR5zFYf0yk0OmDZsGkN7LsLtjNxt3bwTgqKFH5fzM7t9hR9vBvRfN+5tZsmUJh/dPddObPnw65x11Xo9r0won6+t13XEj+o/Ied9iEshiEJE40Az8u6r+QFWbnfeo6vUllC8S48ePp7m5mR07dhQeXEEaGxttb4NRFOIS52NTP5ZxbPSA0fxg/g8AuPrxq3lj/xsZ5113xdff+3WmD58e6D5nHXEW63auoyvZlXa15COhCU4adRLf+9vvpY/dM/+eDBlOuv+knK6Tq2ddzQOvPpCRNqoodbE6vn3Gt3Pet1BpiX+Z9S/8zRF/k1f2E0acABx0hbnXfnnelxnRlHuBnj9pPk9teypDubkJANeefC0LJi3IeW2urKSkJpk0eFLB3eLFIpBiUNWEU+PoyBLLUxTq6+uZPHlypcUwjLJRqO+AX+0fv1IQhXCtibbutkCKIanJvO6PXDuUvcHl7LTRpPpvDvObN9uyCVLSwsWNibR1t2VeWyBl1C8l1/18hVxBueJCQfpKFJMwMYb+wHUichapnckAqqo97SHDMMqGqvpmJXnxUwxRmr+49Yre2PcGQ/oNKTg+qTmCrB5iEutRWtu7gPeL98u0GHIUoPOS78nbnTcIjXWN6RIgQXY9Qw7FEPDafAqtWhXDqc7vWc4P9PCEGYZRboIsdr6KIUDdn2zGDBgDpAKyx484vrBs5LcYwL+2kHe/gJ/FUFDZFNoPEEIxuBZD0GvrJeXu8WZeBb02X9ntalUM5psxjCokiGKoj9X3aDwT9CnWy7ThqYBskH0IkPKtF1rQ4hLv+WSfzLIYkgcthiSF5yzWfoDGeGM6ayuotVEfdzbxeb6j9H1jhZUk+O98rirFICJulTpf68Bzfre3rIVhGOUhyNOob4zBp4xGIRpiqbiCX3ezXLIFWcSzd2d7P1NDvKGnK6mAMitWwxtfV1KBz+O3VyN93wLxnFx9JKrRYvgxKaWQ7xMpqX4L+XotGIZRAoKUZK6P1bO7Yzfb9m9LV0YN0o0sG++mryAUin2AfzXSdJkKYjTEGnil9ZX0uaBxC8gf1A6Cn8VQ6Fo3xrDy7ZXMHJXqfBzUUsml0IK45IpJQcWgqh8ohyCGYUQjyJPslKFTAPjF+l9wzexrgGDlHbJJ1zoKajEkA1gMPkXn3GB0LBZLbUbzPEAHCcTmevIO2/AmisUwZmAqDvPY1sf4zAmp/b9BLZWcaba5dkqXiPLZJoZhlIQgrqSLj7mYQfWDMprsRHEliQj1sfrAFkOhdFVILf45s5KIceKIEzM2qgWxQgq5kgJbDHWNdHR3ZFxb6N4TBk3g1DGnZnyvQa2NfB3cymkxmGIwjD5O0LRTP189hK/Y2RBvyAgG5yNIoDgu8Z4LYVZWklfuXCWrveQL4nrPF8LrSgoTuK6P12eUKAlqbeTs4BbgeywmphgMo48T9Em2Id6Q4QKK2kfYr4FNLoJkJQnSo86TNyvJVUTuwhwmxpCz4U0v0lWDfF9xiWdYQUHvm8uVFMQlV0xMMRhGH6dQuWyXHhZDjj4EhfBLfc1FYIshOxZAZroqkLZSlABZSTmevEOnq9Y1pvdQhLm2LlaXoeyCWiqu0umRZhvAfVZMTDEYRh+nULlsl+zYQK4+BIWoj4eLMQRJLc22GLyxADdF1luzKGoQN3S6avxg8DmsxeBNDw5tMdBT7nI16YEyKgYRuVdEWpyaS37nRURuF5GNIvKiiMzyG2cYRiZBA5vZsYHIrqR4cFdSkPz7uMR9s3Agy2Jw7hkkKyn95O0TxHXvGYTGuvDpqpDayOYXY4hcKymZLFtbz5Qc5eM+IHdZQTgbONr5uRz4Xp6xhmGQeup3S0EXjDHEsmIMAdpU5ppnR9sOnn/red7tetd3zKs7X+X5t56nK9FVcLeviNC8rzn9fm/nXl5656W0bNl7J4LsAnbPe5/aVTVj3iA0xhvpTnbTlexKtzQNcm12jOFA14FA17oKzfu9tne383Lry9W1j6FYqOqTIjIpz5DzgPuddp7PishQERmjqv5dNAzDYNPuTXzusc8BMKh+UN6xDfEG3zINYV1JQ/sN5bm3nuP//OH/8NFjP9qjac/WvVu5+LcXp98PrB+Yd779XfvZ1nKw38Itz93CI5sfSV3bMDCtzNJP7hTOSnItDrf5EMDa1rXc+/K9gWRycYsGdnR38MOXfgjAoIb83zM4riTPzudlby4LdF/XOlr7ztr0sR+89INQMheDaooxjAO8BeObnWM9EJHLRWSFiKyo9p4LhlFK9nftB+D6uddzxoQz8o7Njg1E2fkM8M3Tv8m9H7yXUf1Hsa9zX0+ZOlMyXT3ran70wR+lN3nl4swJZ/a4fuKgifzsnJ9xyphTMhZnV+5CT96jB44GyOhfsLcjVbHny/O+nN6EVgi3y1p7op2GeAPDGodldJrLRTyW6R7rX9cfgAmDJ+S/X10jjfHGdHc7gD0dqW5zN5xyQyCZi0HZLIYA+P3r9K3PpKp3A3cDzJkzxyq8Gocsrh97ypApBV02DbGGHsXoILwraXjTcIY3Dad/XX/f3sauTMccdgxzRs8pON/hA1Jdzdx2vAlNMLBhICeOPBGApniqB4Q3OyhKhVP3ddCmRHDQYmjvbiepSY4eenSkdNWkJhnWOCzQPUf1H5Xh8uvWboY3Ds/Zra4UVJPF0Ax41el4DvZ9MAzDhzAB0YZ4A7s7dve4Nmq2S0xiPbKJwsrkzuO9LrsukLs4pxvmUDgrya/CqfvaPReEfnUp146rGIIq0bpYXYZSCrNBrT5WnxEb6Up0hZK5GFSTYlgMfNLJTjoF2GPxBcPIT5jGM6pKy4EWdrfvzjgetQZPTGK+vZ/DNsNxlYBrwWTXBfIuzu78UZrldGlKMbjVT4PgtVbCKAY/iyGoAs6uhNut3aFkLgblTFf9OfAMcKyINIvIQhG5QkSucIYsATYDG4F7gM+VSzbD6KuE2ck76/BUBnhLWwsQPV3VJZfFEHZ3cQ+LIasukLs4e0tTFJLZr/S1u9i6SiMIrlJq624LVfo6O101kUwETjeti9VlupKS3aFkLgblzEq6tMB5BT5fJnEMoyYIU/tnypBUhVU3YNxbV5LfjuWwMnnHeV1J3mu9fn53XMGeCD4WQxTF4AafOxIdqGpg91id1GVYDEqIa7MthgoohmpyJRmGEZIwT+dumqWbNeQStdSC347lsDJBz01d2XWB0orBk64auIua58nbzcjyZioVoqnOsVa620loIrB1FY+l0lVdJZnQROC04OySI+2J9lAyFwNTDIbRhwkT6HWzWtxNXn4ZRWHw27EcViY4qEBchZJdFyidMtrtcSUVWGRdV9K61nUA3P7C7dz87M3Awb0CQXDHtiXaQpW+dhXC5j2bAcc9ViBrzMvKt1cC8Oz2Z1m2bVlZN7eBKQbD6NOECfROGJRK+nN31YYNEmcTk54NdqLM6y6Y3qdrr3sriivJnXNIvyEAbNi9gRFNI7j5tJsZ2jg0kFzee3d0dwSq0eQybViqN/a2/dvSMge1GPrX90+7jrbu2QrAwhMWBpa5GJhiMIw+TJhFWESYPGRyejdw1OqqLsVSDO6CmbYYVDMCtTFJtfdsSxxMVw0y96j+o9K++q5EF2MGjOHvj/r7QDK5pF1JifZQFsP4QeOBzPpOQa89eujR6deuS+nk0ScHlrkYmGIwjD5M2EX4qKFH8dqu14DoJTFc/Ho1R5Epna7qxhh8Kol6W2wGyUqCzP0AncnOSH5615Xk7mMIajFk13cKc60bfFbVtPwWYzAMIzBh/fnHDT+O5v3N7OnY0+t01exc/agyudaBe51f9o+3L0LQPQHeIG5nojO9WIchrRhCWgzuvbwyB73WVQLd2p2W3xSDYRiBCdtfwC0/sat9V+SSGC4ikrckRlCZemQl+WT/NMYPdlILUhIDMtM+oyoGEaGprin0zucopcK9cgPpqq7eY+XCFINh9GFC9xfIysuH6K6kuMT9XUn0siSGz9O115UUpCQGZCqGrmRXuuFPWPrF+6UUQ4iyFn4WQxTF0J1M7XqOatVFxRSDYfRhwrqDvAuWG3yOuuiISLq8dYZMyXAyZaer+vnjMxRDwEW2TurSZTA6E52R6w25zXqC9K92yW5H2huLodx1kqC6qqsahhGC21bexh9e/wMQzWIIs1j5EZc4m/Zs4raVt/GZEz7DtU9ey77Ofexq3xVKJleGa//3WuaNnsfmPZs55rBjMsY0xZvSG9zW71wfqEJqfbyeZdtSfRA6k52RLQa3vWdLW0twi8G51/qd62nvbmfZm8vS1WILyu3EE/74+h/pSnaVvU4SmMVgGH2WxZsW05ns5O+O/DtG9R8V6JoMi0G1V32EPzzlwwxqGMRvN/+WrXu3smzbMroSXUwYNIFzp5zL2AFjA80z6/BZnDH+DF5pfYUfrf0RAGdPPjtjjNdiOKzxMHZ17Co4r+siS2oycowBHFeSo5RydazrcW/H4mmqa0rvZRjVFOxvNHf0XADW7FhDd7K7IhaDKQbD6KN0J7v5wIQP8PX3fj1wcDK9YSsRbsOWHwsmL+B9496HqqbdQFfNvIrvn/V9bn3frYEXtHEDx7HozEXp96MHjObMiZnNe7yLs6r2sCj8eO+49wKOSybRFVkx1Mfq04Fv7x6DQkwcNJGOREc6znHulHMDXTdpyCSOHHIk7Yn2lMVQ5sAzmGIwjD5LlOJqru+7ozsVY4gaeHaJSQxFQwfBs/HK4TdHxj6GgAXpvL763riS6uMHFUOY77shnuqxHal4X10qC6sr0VX2VFUwxWAYfZZu7Q69aLiK4c7VdwZO+8yHu/u5t3siRCStHPyUlZsyCsEL0rkLcVeyq1fB5/pYPQe6DmTMGYR+8X50JDoipZw21TWxfPvyVPDZFINhGEHpSoR3M7iK4a/7/tprVxKkFnHV3lsMcDAI7VdsrjHeeLC6ajJYQTp3QT3QdQBF04H3sHhdSWEW6YZ4Q4ZiCHNtV7KLpCYrUnIbyqwYRGSBiKwXkY0icr3P+TNEZI+IrHZ+biynfIbRV1DVVGevkIuGG2OA1OJTDFdSkmToMhh+uErKTybXlaSqgfcTuN/N3s69QKo4XRR640rqTHRGciWdOvZUurWbjkRHRSyGsqkiEYkDdwJnkerv/LyILFbVV7KGPqWqHyqXXIbRF3E7k4VNZfQGYNu624rnSurlLmo4uAM6V4xBUTqTnYFLYrgL8Z6OPQD0r4uoGGLRFEO/eD/2de6LpBjc4n3vdr1b866kucBGVd2sqp3AA8B5Zby/YdQEqsp9L98HhC+V4F1kXmh5ofeuJHFcScniuZL8ZHLdQBt2bWBn+85QFoObLjqgfkAkuXqjGDoTnSzdsjT0te7n3du5t+ZdSeOANzzvm51j2ZwqImtEZKmIHFce0Qyj7/DGvje4fdXtNNU1BUrbzMXWvVt7bTEIqXpJYesj+c7lXJvLYgC47snrANjXta/gfO5T9/pd64FeuJI8ynRE04jA17kxhv9t/l8AxgwYE/haV/a9nXtrfh+D37+Y7IaxLwBHqOpJwHeBh3wnErlcRFaIyIodO3YUV0rDqHLcTVa3vO8W3jf+faGvf+lTLzF+YKpfQLHSVd3yGsWwGPyUlasY3tiXerb0K8WRzeQhk4GDrUyjWgxe99uxhx0b+LrGeCMd3alaSZdOvZRhjcMCX+sqhn2d+2reYmgGJnjejwfe9A5Q1b2qut95vQSoF5EeKlpV71bVOao6Z+TIkaWU2TCqDjc7pyneFHmOfItwGNwKq2757WLEGPzmyP6sfr2ms3HjL/s6U9ZFb2IMLmEUaWNdI22JtkiZRa5109bdVvMxhueBo0Vksog0AJcAi70DRGS0OPakiMx15Gsto4yGUfW0tqX+SzTV914x9NpioHjB53yupH51mX2atYezoSduSqvrdiqGKymMq6xfvB8d3R2RFINrMWTfv1yUzUZR1W4RuQr4AxAH7lXVtSJyhXP+LuAi4EoR6QbagEvUrQ1sGAYAO9pS7tOorhEonsUQk1jGPoZeWQwBgs8ufg2CsnEVzN4OJ101qsXg8fGH+XyNdY3p6qphs8e8iqESrqSy3tFxDy3JOnaX5/UdwB3llMkw+hrN+5oBOLz/4ZHnKJbFICKZ+xh64YRwZfGzGLwLJeDbICgbd0F1XUm9yUrKljEIXmUWNoBcaYvBdj4bRh/j/lfuB3oulmHI93Qedh5vSQy3TWcU0hvc/CyGukyLwa9BUDZpi6FzLzGJpXd9h8UbfA5rMbiEXdy91o0pBsMwAtObBSPIE3cQXAthVcsqoEglMfKkq7oECT67MYb9XfsZUDcgshKMGmPwWgx9zZVkisEw+ii9edrfuHsj0PtFx5Xh56/+HICB9QMjz+UqmZ3tO3ucy37aP3FE4aY33sW4N4F673cUNispKl55zWIwDCMv7pP++UedX5T5Fn1gUeFBefC6Vq486UqGNw2PPJerZGaMmtHjnPcJelT/UXz2xM8WnM9baK83gXrvwhzGleRVZmH3mzTEGtL3qvUNboZh9BK37s+xw4JvtMpHbxZyyFwoxw8aX5S5/KwO7yI7fdj0QNVVvbJFzUiCXriSPBZD2CZBIpKW2Vp7GoaRl3fa3gFgeGPvFnSX3sQEINO1UozUV/BffOtidekFOqj7y/vZemMxZASfQyyZ3hhDlCZBrpVkFoNhGHlpbU9tbuvtk75LkCfvfHiVQbGUTK553IU2qGKISSwtX9TNbZAZqI9qMURZ3NOKwWIMhmHkw931XCzF0Fs3hddiKEbqa/acXtwn9zABc1fJ9MZimDDoYCWfqPsYemUxmGIwDCMfacVQLFdSLy0GrzLotcXglsTIIZMbZwizULpKpDfZUt5ry2kxuFaOpasahpGXd9rfoT5Wz+CGwUWZr7eLudeV1Jtdz3DwaTxXrCKKxeAyqGFQZLmiprp6A+ZRnvrNYjAMIxA723YyrHFYr902LkVVDL0MPruF8XIpmNf3vg6E25znNtgplsUQBm+KbZTvxhSDYRiB2Nu5lyH9hhRtvl53cCtiVpJLIffWqP6jQs/ZG4sh6vcdtQSHi6sYzJVkGEZe9nXu69UiV2yKajE4hZQLBXijFA/sTV2pqPQ2fuPuY7B0VcMw8rKvcx+D6mtUMRShC1wuihWsLydpV5KYYjAMIw/rd63vVU5+sfG6oorlSirWPF6GNBbP/VYu3KC3WQyGYeTEDbpWwuecC2+guFiupFzzTB02NfLcUfYRVJpDJsYgIgtEZL2IbBSR633Oi4jc7px/UURmlVM+w6hmOhKpxvJHDj2ywpIcpCRZSTnmmTZsWsa4MIStVVQNpGMMtZyVJCJx4E7gbGA6cKmITM8adjZwtPNzOfC9cslnGNWOm3qZ3eaykhRzg5tLKVxJfdliqIRikHK1VBaRU4GvquoHnff/BqCqt3jGfB/4s6r+3Hm/HjhDVbfnmnfOnDm6YsWK0PL89N++QUudedIMw+i7jGgTPvn/vhjpWhFZqapz/M6V03k1DnjD874ZmBdgzDggQzGIyOWkLAomTpwYSZj6fnU0dRRnk5BhlJNq+1frPloWQy4tME+h88W6pq9QL6V5sC+nYvD722R/qiBjUNW7gbshZTFEEeajX70uymWGYRg1Tzl9Kc3ABM/78cCbEcYYhmEYJaSciuF54GgRmSwiDcAlwOKsMYuBTzrZSacAe/LFFwzDMIziUzZXkqp2i8hVwB+AOHCvqq4VkSuc83cBS4BzgI3AAeDT5ZLPMAzDSFHWnROquoTU4u89dpfntQKfL6dMhmEYRiaWr2kYhmFkYIrBMAzDyMAUg2EYhpGBKQbDMAwjg7KVxCgVIrID2Brx8hHAO0UUp1SYnMWjL8gIJmex6QtyllvGI1R1pN+JPq8YeoOIrMhVK6SaMDmLR1+QEUzOYtMX5KwmGc2VZBiGYWRgisEwDMPI4FBXDHdXWoCAmJzFoy/ICCZnsekLclaNjId0jMEwDMPoyaFuMRiGYRhZmGIwDMMwMjhkFYOILBCR9SKyUUSur7Asr4vISyKyWkRWOMeGicifRGSD8/swz/h/c+ReLyIfLKFc94pIi4i87DkWWi4Rme18vo0icrt4GwWXTs6visg25ztdLSLnVFJOEZkgIk+IyDoRWSsiVzvHq+r7zCNntX2fjSKyXETWOHLe5Byvmu8zj4xV9V36oqqH3A+pst+bgClAA7AGmF5BeV4HRmQd+xZwvfP6euCbzuvpjrz9gMnO54iXSK7TgVnAy72RC1gOnEqqQ99S4OwyyPlV4FqfsRWRExgDzHJeDwJec2Spqu8zj5zV9n0KMNB5XQ88B5xSTd9nHhmr6rv0+zlULYa5wEZV3ayqncADwHkVlimb84AfO69/DPy95/gDqtqhqltI9a6YWwoBVPVJYGdv5BKRMcBgVX1GU//C7/dcU0o5c1EROVV1u6q+4LzeB6wj1c+8qr7PPHLmolJyqqrud97WOz9KFX2feWTMRcX+D2VzqCqGccAbnvfN5P/HX2oU+KOIrBSRy51jh6vTvc75Pco5XmnZw8o1znmdfbwcXCUiLzquJtelUHE5RWQSMJPUE2TVfp9ZckKVfZ8iEheR1UAL8CdVrbrvM4eMUGXfZTaHqmLw889VMm/3NFWdBZwNfF5ETs8zttpkd8klV6Xk/R5wJDAD2A582zleUTlFZCDwIPAFVd2bb2gOeSolZ9V9n6qaUNUZpHrDzxWR4/MMr4icOWSsuu8ym0NVMTQDEzzvxwNvVkgWVPVN53cL8BtSrqG3HRMS53eLM7zSsoeVq9l5nX28pKjq285/yiRwDwfdbRWTU0TqSS22P1PVXzuHq+779JOzGr9PF1XdDfwZWEAVfp/ZMlbzd+lyqCqG54GjRWSyiDQAlwCLKyGIiAwQkUHua2A+8LIjz6ecYZ8CHnZeLwYuEZF+IjIZOJpUYKpchJLLMef3icgpTibFJz3XlAx3cXA4n9R3WjE5nTl/CKxT1f/0nKqq7zOXnFX4fY4UkaHO6ybgb4FXqaLvM5eM1fZd+lLKyHY1/wDnkMq42ATcUEE5ppDKRFgDrHVlAYYDjwEbnN/DPNfc4Mi9nhJmJwA/J2XqdpF6alkYRS5gDql//JuAO3B23JdYzp8ALwEvkvoPN6aScgLvJWX+vwisdn7OqbbvM4+c1fZ9ngiscuR5Gbgx6v+bUsmZR8aq+i79fqwkhmEYhpHBoepKMgzDMHJgisEwDMPIwBSDYRiGkYEpBsMwDCMDUwyGYRhGBqYYDMODiAwVkc953o8VkV+V6F5/LyI35ji33/k9UkR+X4r7G0YuTDEYRiZDgbRiUNU3VfWiEt3rOuC/8g1Q1R3AdhE5rUQyGEYPTDEYRia3Akc6dfL/Q0QmidPnQUQuE5GHROS3IrJFRK4SkWtEZJWIPCsiw5xxR4rI752iiE+JyNTsm4jIMUCHqr7jvJ8sIs+IyPMicnPW8IeAj5f0UxuGB1MMhpHJ9cAmVZ2hqv/X5/zxwMdI1bf5OnBAVWcCz5AqVQCppu7/pKqzgWvxtwpOA17wvF8EfE9VTwbeyhq7AnhfxM9jGKGpq7QAhtHHeEJTfQr2icge4LfO8ZeAE52qpO8BfulpstXPZ54xwA7P+9OAC53XPwG+6TnXAowtjviGURhTDIYRjg7P66TnfZLU/6cYsFtTpZbz0QYMyTqWqz5NozPeMMqCuZIMI5N9pFpaRkJTvQu2iMjFkKpWKiIn+QxdBxzleb+MVJVf6BlPOIaDFTgNo+SYYjAMD6raCiwTkZdF5D8iTvNxYKGIuBVz/drGPgnMlIP+pqtJNWl6np6WxAeA30WUxTBCY9VVDaNCiMgi4Leq+miBcU8C56nqrvJIZhzqmMVgGJXjG0D/fANEZCTwn6YUjHJiFoNhGIaRgVkMhmEYRgamGAzDMIwMTDEYhmEYGZhiMAzDMDIwxWAYhmFk8P8B3zCAXvNCz1sAAAAASUVORK5CYII=\n", "text/plain": [ "
    " ] @@ -789,606 +635,76 @@ } ], "source": [ - "swiftdiff['px'].plot.line(x=\"time (y)\")" + "fig, ax = plt.subplots()\n", + "swiftdiff['rmag'].sel(id=plidx).plot.line(ax=ax, x=\"time (d)\")\n", + "ax.set_ylabel(\"$|\\mathbf{r}_{swiftest} - \\mathbf{r}_{swifter}|$\")\n", + "#legend = ax.legend()\n", + "#legend.remove()\n", + "fig.savefig(\"rmvs_swifter_comparison-mars_ejecta-rmag.png\", facecolor='white', transparent=False, dpi=300)" ] }, { "cell_type": "code", - "execution_count": 8, + "execution_count": 23, "metadata": {}, "outputs": [ { "data": { - "text/html": [ - "
    \n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "
    <xarray.Dataset>\n",
    -       "Dimensions:   (id: 18, time (y): 5)\n",
    -       "Coordinates:\n",
    -       "  * id        (id) int64 2 3 4 5 6 7 ... 6000205 6000206 6000207 6000208 6000209\n",
    -       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n",
    -       "Data variables:\n",
    -       "    Mass      (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n",
    -       "    Radius    (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n",
    -       "    px        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.8965 0.0579 2.692\n",
    -       "    py        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -2.79 0.004403 -2.209\n",
    -       "    pz        (time (y), id) float64 0.0 0.0 0.0 ... -0.1535 0.002265 0.06436\n",
    -       "    vx        (time (y), id) float64 0.0 0.0 0.0 ... 0.02465 -6.413e-05 0.01524\n",
    -       "    vy        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.00847 0.00055 0.0203\n",
    -       "    vz        (time (y), id) float64 0.0 0.0 0.0 ... -1.248e-06 0.001067
    " - ], + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAERCAYAAAB4jRxOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABCvElEQVR4nO29eZxcVbW3/6waekg6IZCBJISQgEBCAoQkTMJlUoYADkDQi+CACOqrXvgpzu8VUa+ivigRRQURFXzlvlfDpAQZEmQmJiSBIAQIBNKEzHN6rtq/P6pO9amqM9fYXevhkw/VVefsveokvc4637X22mKMQVEURRn8xGptgKIoilId1OEriqI0COrwFUVRGgR1+IqiKA2COnxFUZQGQR2+oihKg1D3Dl9EfisiG0VkZRnGmiEiT4vIiyLyvIh82PbZH0VklYiszM6ZLHU+RVGUeqLuHT7wO+CsMo3VAXzMGDMtO+YNIjIi+9kfgSnA4UAr8KkyzakoilIX1L3DN8Y8Bmy1vyciB4nIAyKyVEQeF5EpAcd6xRjzavb1OmAjMDr78/0mC7AYmFDWL6IoilJj6t7hu3Az8AVjzCzgauCmsAOIyDFAE7C64P0k8FHggTLYqSiKUjckam1AWESkDXg38D8iYr3dnP3sfOA7Dqe9bYw50zbGOOB24OPGmHTBsTcBjxljHi+37YqiKLVkwDl8Mk8l240xMwo/MMbMB+Z7nSwiw4G/Af/bGPNMwWfXkJF4Pl02axVFUeqEASfpGGN2Am+IyIUAkuHIIOeKSBNwF/AHY8z/FHz2KeBM4CKHqF9RFGXAI/XeLVNE/gScAowCNgDXAAuBXwLjgCRwpzHGScopHOsS4DbgRdvbnzDGLBeRPuBNYFf2/flBxlQURRko1L3DVxRFUcrDgJN0FEVRlGjUddJ21KhRZtKkSbU2Q1EUZcCwdOnSzcaY0U6f1bXDnzRpEkuWLKm1GYqiKAMGEXnT7TOVdBRFURoEdfiKoigNgjp8RVGUBqGuNXwnent7aW9vp6urq9amuNLS0sKECRNIJrXDsqIo9cOAc/jt7e0MGzaMSZMmYeulUzcYY9iyZQvt7e1Mnjy51uYoiqLkGHCSTldXFyNHjqxLZw8gIowcObKun0AURWlMBpzDB+rW2VvUu32KojQmA07SURRFKQedfZ3835f+L519nQC0Jdu4+LCLScYGb+6tIR3+u9/9bp566qmi9z/xiU9w7rnnMnfu3BpYpShKNXluw3Pc8NwNee/N2ncWh48+vDYGVYEBKemUipOzVxSlsUiZFAB/OudP/Pq9v857b7DSkBF+W1sbu3fvxhjDF77wBRYuXMjkyZPRzqGK0jhYv++C5PJu6UG+FUZDRvgWd911F6tWreKFF17glltu0chfURoIQzbAE4hJxhWqwx/EPPbYY1x00UXE43HGjx/PaaedVmuTFEWpEpZzF0QdfqOgJZSK0phYEb5k/wNIow5/0HLSSSdx5513kkqleOedd1i0aFGtTVIUpVpYio40ToTfkElbi/POO4+FCxdy+OGHc8ghh3DyySfX2iRFUaqEPcK3HP5gL9yoqsMXkTVkNglPAX3GmNnVnN9i9+7dlj38/Oc/r4UJiqLUmFzSlsZJ2tYiwj/VGLO5BvMqiqLkyJVl2iQd+01gMNLQGr6iKI1LXtJW6/ArggEeFJGlInKF0wEicoWILBGRJZs2baqyeYqiNAp5Gj6NIelU2+GfYIyZCcwBPiciJxUeYIy52Rgz2xgze/Rox43XFUVRSqcBq3Sq6vCNMeuy/98I3AUcU835FUVRLFTSqSAiMlREhlmvgTOAldWaX1EUxU6uBFMgLnFAF16Vk32BJ0RkBbAY+Jsx5oEqzl8W1q5dy6mnnsrUqVOZNm0a8+bNq7VJiqJEwCnC1zr8MmGMeR04slrzVYpEIsH111/PzJkz2bVrF7NmzeL000/nsMMOq7VpiqKEwN5LJ9daQSUdxc64ceOYOXMmAMOGDWPq1Km8/fbbNbZKUZSoNFLSdkC3Vrj2vhf517qdZR3zsPHDueZ90wIdu2bNGpYtW8axxx5bVhsURak8lqQTI0Y2wNeFV4ozu3fv5oILLuCGG25g+PDhtTZHUZSQ2JO2VoT/6xW/rqFFlWdAR/hBI/Fy09vbywUXXMDFF1/M+eefXxMbFEUpDXvSdsyQMQC8teutWppUcTTCD4kxhssuu4ypU6fyxS9+sdbmKIoSEXsvnWQsyacO/xSJ2ICOgX1Rhx+SJ598kttvv52FCxcyY8YMZsyYwf33319rsxRFiYhVoSOIlmUq+Zx44omD/h+FojQCdkkHMpG+Jm0VRVEGIXZJBzKJ28FelqkOX1GUhqQwmrci/cH8BK8OX1GUhsRJ0rG/PxhRh68oSkNSKOk0QnsFdfiKojQkOYdPv4YPGuEriqIMOnKSjhQ4fNXwFYuuri6OOeYYjjzySKZNm8Y111xTa5MURYmAWyQ/mCUdrcMPSXNzMwsXLqStrY3e3l5OPPFE5syZw3HHHVdr0xRFCYFKOoovIkJbWxuQ6anT29ubeyRUFGXgUCTpMPglnYEd4S/4Gqx/obxjjj0c5lzneUgqlWLWrFm89tprfO5zn9P2yIoygCksyxzMko5G+BGIx+MsX76c9vZ2Fi9ezMqVujWvogw0CiWd3MKrQSzpDOwI3ycSrzQjRozglFNO4YEHHmD69Ok1tUVRlHAUSjoa4StFbNq0ie3btwPQ2dnJww8/zJQpU2prlKIooXHqpTPYGdgRfg145513+PjHP04qlSKdTvOhD32Ic889t9ZmKYoSkqLWCg2w0lYdfkiOOOIIli1bVmszFEUpE42UtFWHryhKzfjtyt/y82U/z3vv8FGH8/s5v6/YnJs6NjH3vrls794O9Es5ubJMl6Ttfavv49qnry26IbQkWrjtzNs4dJ9DA83/vWe+x/xX5+e999kjP8vlR1we5mtEQh2+oig1Y9XWVbQkWvjQIR8CYPH6xby09aWKzrmhYwNbu7Zy+gGnc9y44xiSHALYumW61OG/uv1V+tJ9fGLaJ3Lvbenawt2v3U377vbADv/lrS8zunU0cybPAeDPr/6ZV7a9UsI3Co46fEVRaobBsHfz3lw16yoAfrLkJ7y67dXKzpl16B981wc5acJJufd92yMbSMaSOVsh47zvfu1uwlRyGgwHDD8gN86itYuqVgo6+NPSiqLUL4a8leoiUnEN3c25WpKO2/xpky5aVR+pdr/wO1P572yhDl9RlJphMDmnCdltBqmOw7fPC/6STpp00TmFYwad3z6OSPU2T1eHryhKzXDaZrDiEb6Lc/WL1o0xRbX6fjcJ1/lt942YxAavpCMicRFZJiJ/rfbc5SSVSnHUUUdpDb6ilIAxJk/eqOZG4oXyjOXM3eYvjMwhmqRTFOEPcknnSqCyafgqMG/ePKZOnVprMxRlQGMwOe0cqrsJSWhJx0HDj9JS2ZD/pBCT2OCUdERkAnAO8Jtqzltu2tvb+dvf/sanPvWpWpuiKAOewqQtVHbxk6uGH0XSscYIk7M1xU8K1ZJ0ql2WeQPwFWBYOQb74eIf8vLWl8sxVI4p+0zhq8d81fOYq666ih/96Efs2rWrrHMrSqNRGNnGJQ5kEqRx4pWdsyD/6nezcZJ0+v19OIddlKgebJKOiJwLbDTGLPU57goRWSIiSzZt2lQl64Lz17/+lTFjxjBr1qxam6IoA4L1e9Zz+O8P52+v/63oM0Oxhg+VlXT8yjK9Iny3sswwDttQu6RtNSP8E4D3i8jZQAswXETuMMZcYj/IGHMzcDPA7NmzPa+CXyReCZ588knuvfde7r//frq6uti5cyeXXHIJd9xxR9VtUZSBgLWK9L7X7+OcA8/J+yxt0kUJTOv9SlHYB9/C72aTJu0q6YRx2E7fedBp+MaYrxtjJhhjJgH/DiwsdPYDgR/84Ae0t7ezZs0a7rzzTk477TR19ooSEac6fKhOA7OirUmzP7pKOg7ae6SyTIc6/EEn6SiKohRRsOq0Gg7fLWnrK+ngLumEmr9AGhJkUEo6OYwxjwKP1mLucnLKKadwyimn1NoMRalrvKJfp5p0oOKrbe1zWfjdbAqlGPsYpSZtB52koyhKY/LX1zNrLJtjzUWfFTrKeCxTmVPRpG3BTlcWfs47bYo1/FyVTsiVtkX9g6pwgwN1+IqiVBjLiR8++vCiz5zkDaiOpFNIED2+HElbp3YSGuErijIoSKVTrp/VMmlbiJ+c5GRTOZK2g7IOX1GUxiRlMg6/L91X9Jmbw69kEtM1aWtF767t8D1W2oaZ3+Wpphqow1cUpaJs6dwCZBZg7eopWJ1u8mUSyxE63RzKhauGn/15R/cOx/OcNPyoNyh7/6C6KssUkYkB/wyvhsGKogwslm1cBsBfXv0LX/7Hl/M+K1x12hRrAuDBNQ9WzB63CH9IIrPV4Q//+UPn8xzq8O2fhZq/YKVttZK2Qcoyf0/mIcfrucMAvwP+UAab6p5JkyYxbNgw4vE4iUSCJUuW1NokRalLuvq6MBiOHns0u3t2s7NnZ97nhU70lP1PAfploIrg4ptn7juTfVr2IRlLupzmXocfKmlrHEpRq5Oz9Xf4xphTC98TkbHGmPWVMWlgsGjRIkaNGlVrMxSlrlm/J+MmznvXedz/xv1s69qW93mhht8cz5Ru1mKlbUxizBg9g/bd7Y7HO9bhR+juWXjjqCtJx4WPldUKRVEGJes7Mg5/7NCxjtUoThugQG2Sttb8Xs63HElbx146db7S9gMi0gE8ZIxZVU6DwrD++9+n+6XytkdunjqFsd/4hucxIsIZZ5yBiPDpT3+aK664oqw2KMpgYeXmlQCMHTKWGMVdIZ36ykBt6vCt+cNsgBKlLBNq1x45qsM/HzgKOE9E3mWMaaidQJ588knGjx/Pxo0bOf3005kyZQonnXRSrc1SlLrj7tfuBmDM0DGO0kWh87WqV6rSLbOweRreCVQnSSc3ZlgNv0DSqesI3xizAXgg+6dm+EXilWL8+PEAjBkzhvPOO4/Fixerw1cUBwThiNFH0BxvJi7x4JJODfrhQ+aG47UBSkOutBWRX4jI77KvzyirRXXOnj17cjtd7dmzhwcffJDp06fX2CpFqU+6Ul1MGj4JcJZLnCQdoTq9ZRw1/JiHw69Qe+R6K8t0ogfYkH19GlC5otk6Y8OGDZx33nkA9PX18ZGPfISzzjqrxlYpSn3Sk+qhJd4CODs2YwyxWPFipqpU6Tg5fI8I31HDj1qWWdgeuUoRflSH3wHsJSJJYGIZ7al7DjzwQFasWFFrMxSl7jHGsLVrK82JTKmlmzN1ipprpeF7JW0NJm+FLERsreDwVFPvDn8r0An8AniyfOYoijJYeGnrS0B/8zS3pK3TRiS12AAFvOUVxz1to1TpOGz6UpftkUVkhIjcBlyQfesPwOyyW6UoyoDH6pvznonvAZylGmPy2wxYx1WjH75TcO4lJzmttLV/Fnh+h01f6jLCN8ZsF5HrgEnAZuAIYH4F7FIUZYDTm+4F6Jd0HBy5o0xSaUnHqw4f97nTJl1ka5SqIqeVtvVclnkZ8IYx5u/A0jLboyjKIKEn1QOQ603jmLR1iJrjEq9ZlY5T6aiFo6RTpl469bzwahvwGRE5FFgBLDfGLCuvWYqiDHSsCD/P4Rc6NuOctK1GHb6Tw/dL2laqLLMuJR0AY8wPROQR4BVgBnASoA5fUZQ8rAi/KZ5peeykVadN2lHDr2jEa0n4UVbaVqAsMybFLScqReiFVyLyHeADwOnA28aYeWW3qs7Zvn07c+fOZcqUKUydOpWnn3661iYpSt1hbWLiFeHXZZWOl4ZfuIl5wZhh5i8cuxqEdvjGmG8B3dlzLxCRW8puVZ1z5ZVXctZZZ/Hyyy+zYsUKpk6dWmuTFKXusCQdK8IP6vBrnbT13MTcIcEMZZB06jXCz/JbYCowEripfObUPzt37uSxxx7jsssuA6CpqYkRI0bU1ihFqUMKk7ZOLROcEqHVrEsvJCYx181XnOSncq20reekLcB/kGmvkADmkdHxq87j/+8VNq/dXdYxR+3fxr996BDXz19//XVGjx7NpZdeyooVK5g1axbz5s1j6NChZbVDUcLyxo43eOLtJ/LeO2TvQzh23LGBzu/q6+Le1ffSnerOe3//YfvndqIKyivbXuHZ9c8C/Q4/Houzo3sHvalekvEkT7z9BFu6tjBmyJii8+997V6unn01Q5PFv1epdIp7V9/L7t7M735c4pw9+WxGtIwIZJvXStu4xOns62TtzrXsP3z/3Purt69m+ablHD/u+LzjrTH+e9V/s++QfTl78tnEY3HXuR9d+yg9qZ6iCH9jx0bHm1+5ierwVwMHA/cYY/6/MtpT9/T19fHcc89x4403cuyxx3LllVdy3XXX8d3vfrfWpikNzi+X/5IFaxbkvTeqdRSLPrQo0PlPrXuK7z5T/O84LnGWXLKERCy4u/jp0p/yxNtPMLp1dG4XK8vJLd24lOPGHccXH/0inX2duYVZFm3JNjZ3bubRtY9yzoHnFI29atsqvvXUt/LeS5kUHz3so4Fs89Lw92vbD4BfP/9rvnfi93Lv37jsRgDGt43POz4hmWuyfs96vvHENxjeNJyT9z/Zcd6+dB9XLrqStEnnjdPR2wHA2l1rmTi8sp1qojr8F4G1wGUi8mNjzNFltCkwXpF4pZgwYQITJkzg2GMzUdPcuXO57rrrqm6HohTSm+5l8l6TuePsO4CM033ozYdCnQ9w+5zbOXDEgZnX/7qdX634FSmTIhHCXfSmejli9BH87szf5SLeOZPncMdLd9Cb6s0d87HDPsbVs6/OO/em997E2fPPzklCbnZef/L1zBgzg/f8z3ty7wXBy+F/eMqH+fXzv6YnnT93T6qHA4YfwDXHX5P3fjwW5xfv+QWfe+RzAOzo2eE+rzGkTZrPHPkZPj7t47n3TzvgNBasWRDqO0QlqsM/iEw9/s3Z/zcMY8eOZf/992fVqlUceuihPPLIIxx22GG1NktRMBgSsQTDm4YD0BJvCb3XKsCwpmG5MYYmMpJKKp0Cd6WiiDRpmmJNJOP9G4IXbl9oMDTHm4tkDEsCctPSLUlmaHJozs6wvWzAWdKxxi28bmnSDG8a7niO9QQD3s3UrO9tPx6qs+mLRVSHv9YYs1BExgEby2nQQODGG2/k4osvpqenhwMPPJDbbrut1iYpSvEKzrALmHItZvrHsGScsNFnKp0qkoAKty90600Tl3jecW4I0l8lU8YqlyAbtdhxK9UsxM1G6/x6dvhnicgrZLplvkkmieuJiLQAjwHN2Xn/bIy5xvus+mTGjBksWbKk1mYoSh5F5X64V5y4nQ/kVaJYckyYcXK2FFbfZCNZayy3LQP95sw5RonenjhzuosDj8VyHT7tcxaWZOaOtzl8r6RrLlnstlq3CqWZUcsyRwBfBb5CpiY/CN3AacaYI8ms0D1LRI6LOL+iKAUUOtmwS/YtR2p3SFa0XegAg4wVpNGYV4TvNqfdYVu2hpKuPLplWvMXSToei66CRvi56+tzI6wkUSP87wBTjDGrRCSQlSZzla0aymT2T3VWGyhKA+Ak6UTR8O1jWHq6tWo2jC2FjtAuXbhFu/bj/DR8Eck57bALn9zmtuYvnNvL4Yd9yig83nqiqUY/ncARvogcab02xrQbYx7Ovv5aiDHiIrKcjO7/kDHmWYdjrhCRJSKyZNOmTUGHVpSGp1ASCLuAyak+3XJGfSacw3fqO5Nz+KQ9na6fhu8U4UdpbeDmqEuJ8IMkbYsknQhPKVEJI+ksE5HnReQrIrK//+HFGGNSxpgZwATgGBEp2v3bGHOzMWa2MWb26NGjo0yjKA1JufZKLYukQ7GDtLch8JJVguYNIjt8n2vi2gKi1KSty4KvaiZtwzj864GhwHXAGyKySEQ+GWVSY8x24FFAd/9WlDJRmLSNx9x7u7udDxSNAfD6jtdD2eIUEdvLD/0amIGHhm9znDnnGUEN8XLgjhG+i7ssvMm64XZTqsukrTHmy8aYg8hsafgbMu0Ubg56voiMFpER2detwHuBl0NZqyiKK05VOgYTOMp3ironD58MwPJNy0PZ4pW09XP4QSUd+/lRnGUYDd8pJ1FoL/hU6bh857qM8EVkpIh8Cvg+cCmZfxZrQ8w1DlgkIs8D/ySj4f81jLH1wKpVq5gxY0buz/Dhw7nhhhtqbZaiFG2OXVj37n96sUM6dJ9DGZocGjpp6xTh59nj05MePJK2dg0/QnQcVcN3c+ZBk7aukk6dLrxaT+YGsQ24DbjDGPOE9yn9GGOeB44KZ179ceihh7J8+XIAUqkU++23H+edd15tjVIUnNvuQkZPjwdYJuvmkJrjzXT3Ba2+7h/LS6v2c7oJSXhuNVhoZ6gqHZ+yTMcqHYechP14iyDOvxy7ZkUljMO/C7gDWGCMqXzThwHAI488wkEHHcQBBxxQa1MUJeMwbL4kygbbUOyQmuPNdKW6Qo3h5CCtSNZk/wNvHd2tMsg61xpfCLcJeOQqnQALrzzn9Una1lUdvjHmQ5U0JAqLfnczG98Ml0zyY8wBB3LqJ64IdOydd97JRRddVNb5FSUqhRF+2HI/N0fYHG9m1dZV4WwxxlvD96jDh2zCOe1id0ELiHLvgRu6Dj9gS2O362vlAOqqDl/Jp6enh3vvvZcLL7yw1qYoCuBQpROwJ03ufJcIdFfPrtyuVUFJmRSxWICkrUeE76fhW4Rd+OTVDx8y181p711XScfuRj1McfvOudxGds3Eso3LWPBGfpvrchF6pa2IvM8Yc18ljAlL0Ei8EixYsICZM2ey77771swGRcnDJWkbVO5wO27GmBm8ufPNUKY4SSBOdfhelTJee8vaxwsr6eTscZlbRByrdLxuEH5jWmM4HVOYtL3ntXv4R/s/mDN5js83CE+UCP+/ym7FAORPf/qTyjlKXeGWtA2qDTv10gEYkhhCZ19nOFs8krYpkwqUtA1SpZN9UfJ6AzuOGr5H0jYveexx4/Ery7RuCClT3Gm0XERx+OHb0w0yOjo6eOihhzj//PNrbYqi5HBqrQARkrYFjro10Rra4fsmbX1kFS9JJ2cn0SL8KHM75STsx/cfGH7ewvLZVDqV99RQTqI4/IZveDZkyBC2bNnCXnvtVWtTFCVH4eIgyyF++6lv89bOt4qO39y5ma8//nUea38sdz4UV50MSYaL8Nt3tbN+z3rXOvz5r87nNy/8xnMMEWH5xuWOnxU6zpjEyuqVCiP8Hd07WLNzTaDWCl5PGn798BeuXQhkGtXVU4SvKEodUtjv5cjRR3LI3ofw8FsP8/jbjxcdv3TDUv76+l+54bkbcuc7YUX4QWWTp9Y9BWQ2ULfTlmzj3/b7NzZ0bODWlbcC7rLKrp5dRTtDWRRKI1HLMt0o7If/wuYXABjd6tzba2TLSMYOHes7tluEP6FtAgDrdq8DMo3qrL1yy406fEUZJBS2R542ahq/O+t3gHNfGssBWXvMumnMQxJDAOjqC1aLb4175qQz896Px+Lc9N6b+MBBH8i95xY1HzfuON+FVzkJP2RZZlgN33pd+H0skvEkt55xa96xYeZta2pj1r6zcp/3pftyPYzKTRSHv6HsViiKUjKFSVuwdbt00MOdtvGDYic8JJlx+Lt6dgWywyovDNI/3jVpG0v4Lryyn1tJDd/veLstQexw+s5C/94FdZW0NcacXglDFEUpjcKVtuC9gUlhr3y3CHTc0HEArNuzLpAdluMqZUvAuMRd+/eUKunk5vaI8O1PDH5rBuyfeT1peN047LuT9aVV0lEUxQfHCN9jAxM351Q4xv7DMttfrN0VrFeib0QcoM4vHou79+AvaLwWdaWtazO0gjp8t3JVp7GiSDrWe9bnqXSqriQdRVHqECeHH5MYgjg6z6CSzughmWTlls4tgezIRfgBJB23Y5KxZOA6/Kgrbd0o1PALe/c44fY0kzevx5OC/abVZ+qsSkdEvmh7fWj5zBkY/PSnP2XatGlMnz6diy66iK6ucI2lFKUiGGdnkoglnCUdh12dnGhNtAKwp3dPIDP8HH6Q7pJxidObdu7RWOg4y16l46bhlxrh++zja0lsfem++qjDF5ERInIbcKGI/C8ROREIvKftYODtt9/mZz/7GUuWLGHlypWkUinuvPPOWpulKI4RPmQcfilJ25jEGJIYEtzhlylp6ybpFEksErI9csQqnSBJW689hD0jfPoj/FS6cknbUKNmtya8VETOBDYDRwDzK2BXXdPX10dnZyfJZJKOjg7Gjx9fa5MUxTFpC+4J0KBJW4ChyaF09HUEssMvwrfjlbQNI+lUtErHknQ84uNAq5oLunzaKZJ0KpS0jTpqrzFmqYisAzaW06AwbL9vNT3rgkUdQWkaP5QR7zvI9fP99tuPq6++mokTJ9La2soZZ5zBGWecUVYbFCUKbhG+INzx0h185eiveG4a4iU5bOrcxPxX5/OuEe/io4d91NsO4+0gg9wIErEE27u3Ozctc0jalruXTmdfJ92pbprjza4rkO0EqdLxiv4LJZ260vCBs0RkAvAr4KdltKfu2bZtG/fccw9vvPEG69atY8+ePdxxxx21NktRXB3+8ObhQHEtvpuG7yVd/OifP/K1w08CceroWYj1NOG08UqpSVs/LGf7zu53ANt18pgmkKTj8WSRJ+mYylXpRL2NjAC+CnwF+FTZrAlrhEckXikefvhhJk+ezOjRmcqF888/n6eeeopLLrmk6rYoih23Fr4XHnIhNzx3A73p3rzI0dXhl+hAy6HhW5unO0XuRSttiViW6TL3lH2mANCT7snMVyZJx7MsU/plqbpJ2tr4DnC3MWYVUPl9ueqIiRMn8swzz9DR0YExhkceeYSpU6fW2ixFAZydidviK7ekban41a3nRfgux3i1di6K8KW8Gr612YtVJRQmaRulPbL1nr1bZr1JOl8HLCFvUZlsGRAce+yxzJ07l5kzZ3L44YeTTqe54orabcSiKBZukk4ynnH4hWWOUSSdIFi7Q7lKOvhLOpak4bTNYcmtFXyOtW6QVo8hC68IP9BG5AW5h7yxJZYX4ddFlY6NHvp76pwK/L085gwMrr32Wq699tpam6EoebhV6VjOo9CBuTm+UiUdr97xheN7JU7Be4VwXh1+mLLMCkT49u0bXef1iPDzWivUYbfMDmAvEUkCE8toj6IoEXGN8C1JxwSTdJwc23dP+C5Arg2wF2mTDuQcvQjiQO2SThTcbjbW9epJ9eTZULKkU5B7KMSap966ZQJcA6wGfgH8sXzmKIoSFT+HXxjhh0nafvBdH+T9B70/UAsBr+0AIViVjtcG7JXuh5+7Xungkk6pSVu7pJMyldvxKupzw38YY34CjdlaQVHqkcIdryxykk7aW9Lxkzq8WhYHscMiiKTjmbQtsR++39xWziNUhG+1Vgiy0tYlaWtfaWvddMpNKIcvIiOAXwIHiEgXsIJMWeal5TdNUZQwuJVlukWsha0LgixIcu1gaR/XpAL1nQF3ece6STlu3FLQzKzcK20Lr5flxD2/k1WHH6CXjtN3ti8eq6SkE7q1goi0A48BzwJH0oCtFRSlnBhj+OTfP8manWscP5+17yz+z8n/x3OMd3a/w+odq5m81+Siz9zKMpduWArAmp1rOPX/nZrrleOlbW/t2so588/h5jNuZr+2/YqOeX3769z+r9tdtyf0Gt+OV4T//We/nzkmK7HEJMaCNxbw9Lqn+e4J3+X48cd7jv1o+6OedlhJ29XbVwPejrrQ3nnPzePosUdz5Ogj8z6/atFVuevtVaVjjKlot8woo24BPgMcSibCby+rRYrSYKRMiiUbljBt5DSmjsxf07FswzIWv7PYd4y3dmU2KS88H9zLMtua2gCYe8jc3HuThk9yjS7PO/g81u9Zz8K1C3lz55vODn/H60BG83fD7jiPHnu04zFeGn5nXyct8ZZcAvnzR32eJeuXcNdrd/HS1pd8Hb7V/dPayauQMa1jAGhJtADBdrxqijfx0cM+yu3/up1VW1cVOfxn33mWsUPHcu6B5zJ739lF51t1+NYNrm40fGPMdSKyEHgFmAH8G7DM7zwR2R/4AzAWSAM3G2PmhZ2/Hpg3bx633HILxhguv/xyrrrqqlqbpAxgLDnitImnccUR+Ws6rlt8Hfe+dm/gMWaOmVn0mVXi59Ru+F0j3sU1x18TyM4p+0zh0umXsnDtQlfN3HLQHz70w4HGHNU6yvF9rwg/IQk+ethHcw74/Qe9nzMnncldr90VqKeOwTB95HTXzwtzHkFXIH9y+ie5/V+3O16btEnz7vHv5stHf9nx3MLvWzcRvoh8B4gDy4HlxphHA57aB3zJGPOciAwDlorIQ8aYf4W1oZasXLmSW265hcWLF9PU1MRZZ53FOeecw8EHH1xr05SBikcXxbjEAydKwTkKzUX4PlU6QfCKvMG/rQIE2wAlt/DKYZ40xWWfgbpVWrjsG5CzT4RkLJm7XkG7f/qtDvaK2nMRfjZnUTd1+MaYbwE/A3YBF4jILQHPe8cY81z29S7gJaD4mbDOeemllzjuuOMYMmQIiUSCk08+mbvuuqvWZikDGK+Eo1dfeKcxwrRWcEvyepFzrC42Wytjg9bh+y28ctupq9D5Wnq+W0tlO27lq3YSsUToCN+y2enapNIp35tM2qRzc9ZF0tbGp4FfG2MeiHKyiEwCjiKT+C387ArgCsj0rfFiwYIFrF+/PooJrowdO5Y5c+a4fj59+nS++c1vsmXLFlpbW7n//vuZPbtYk1OUoHi1NPDqC18wiOsYuYVE2WZg9nmD1NXbscZ33ZwkG+F7RrMB6vDdomW3DdLDRPhuK5LtJGPJfocfQMO3f+54k/JZm2AlbetO0snyW+CzIjIU+KMxZnnQE0WkDfgLcJUxZmfh58aYm4GbAWbPnl2ebk5lZOrUqXz1q1/l9NNPp62tjSOPPJJEojJ/OYpi7VblF417RaFW8rGrL7/VcCkRvlu9uV8v/KC4SUduNfFB6uBzNgaI8JOxZO6JKEgdvt1mpwjfGO85rTp8a866Sdpm+Q8y/XQSZOSdk4KclG3F8BcyN4mSyzm9IvFKctlll3HZZZcB8I1vfIMJEybUxA5lcOC18YgV6fWZPpLivhjHSxYakshUoxRuURjE8RWSayHgk7QttbWCJWkURvjWvE4OMSaxwElbX4cfTxZJOn43MWtMt+0k/VYfG0zu6aBSC6+i3oZXAy3APcaYoM5egFuBl6xVugOVjRszm3y99dZbzJ8/n4suuqjGFikDGT9JB9wllNwYHrLD0ORQgKItCv163jjh1+MmSIIzyE3G7XvnFkE51bITK0vSFpwlHb8blVui2RiTkc+8JB0yN6tchF9nGv6LwFrgMhH5sTHGuZg2nxPItFR+QUSWZ9/7hjHm/og21IwLLriALVu2kEwm+cUvfsHee+9da5OUAYyXk8qtOPXR8b0kneZ4M3GJ09HbUXROWOnFV9LB3zmWsvDK64YSdKvDwEnbVP5K26A2R2k7bbWHsCqy6k3SOQjYRkZr3xbkBGPME/imSgYGjz/+eK1NUAYhXk7oxc0vcsy4Y3zHcHOEQxJDiiWdUqp0XG5SloMO2lrBj7W71nI8/QupvHIEcYkHd/gViPDdKoWCLKay2kPkyjLrLGm71hizUETGUcNNzBVlMOHkJKePyiwQ2tq91fNcrzwAZFaVlkPD96o1t9vhJUkEcfjWgqzCFg25G4rTvrBBI3yfBCoUOPwATy32z902h/er0kmbdC7Cr5s6/CwNu4m5opQbr0d+y/H5afh+G20PTQ4t0vAjRfh4R/h+2xv6fWZhRbiFyWgvSScmsUDyS5DchWOVjo/dIoIgrjKU31OPMZWP8KM6/BH0b2LeXTZrAlKuvTcrRb3bp9QXXv9ecjs/pb1X2/otDmpJtBSVZaZJh6/SscofS0jaBqnScasG8oqWAydt8XfeeVU6AevwwVlWCprINpiKL7wK7PBFxN4N6DtkKnSqvol5S0sLW7ZsqVunaoxhy5YttLS01NoUZYDgt9IWSkvaArTEW+hOFcRmJpjzteNXpVPupG1RhO/RqrisSVtJhGqPbLfbLWnrV5Zpb55WDxr+MhFZCdwB/MkY8zCAMeZrFbHMhQkTJtDe3s6mTZuqOW0oWlpatDZfCYyXs87V4ftE+F4rbSET4e/u2Z33XtDqEzu+rRUCLlIKilu0XFIdfpCVtvEkvd3Z3kOmf3w/nGwIIunEyKy0raeFV9cD5wPXAd8XkceB240xv62IZS4kk0kmTy7u+a0oAx2vOvySJZ14C5tS+UGS385UTgSuww+wHaAXbnvEeskr5Vx41RRrKo7wA9zEvBy+l0xjbWJeNxq+MebLxpiDgNnAb8isrr25IlYpSgPhVWHjtuK0aAyfHZ+aE8109+VLOjVbeBXQcUKxhl+OpG0ghx9vytviMGiuI2qEj1CVKp3Ao4rISOA8YC5wasZE3qqIVYrSgDg5QusX3zfC90kstiZai5K2EL7njd9WfuVaaWt9j1AOv4wrbZviTbmcR5gnoZjEQi0Wy51XKOnUwUrb9WSeCLYBtwF3ZBdTKYpSAkFW2jo567wxAkg6nanOvPfSJh16KaTltLZ0baGzrzO3e5SFVetfaoRfKOkYY3h799u8s+edvM8L2dWzy3fssJJOb7o3WMdSnCN8axy/a9Kb7s05/JpLOsBdZCL8ccaYz6izV5Ty4CXHWBr+A2u8O5H7OXwnSSfKwqumeBOC8KsVv+KKB68o+vzWlbcC5YvwLed568pbmTN/Dp/8+yeB/g6gdnrSPTz45oO+YwdJ2toj/MXr/beYtHBy+NcvuR4oXkRWOB/AnavuBGoo6YiI1ZT+6uz/x7ncobc7tTtWFMUbL2cdj8XZp2UfRjSP8BsEcHe0rfFWetI9pNKpnFzg19DLiaHJodxyxi3cuOxGtnYVr/4d0TyC0a2jPSPUQA6/IMLf3LmZ1kQr/3ncf5KMJzl1/1OLzpk2ahrPrHvGd+wgPYTsGv7IlpEMaxrmOy441+F3pTJPZ2dMOsP1vIsOvYiblt/Exo5M44JaSjq/J/fPyfVvygC/I7NnraIoEXCTOg7Z+xA6+zodP7PIPSW4/IZaEXF3qpshsUy75CAtBpw4dtyx3P3a3Wzu3Oz4eeEG3oUELW+005fuoyXewvsOep/rOYfufWigDd+D9NJpjjfTm+7FGEOaNAcMO8B3XHBeC5A2aabuM7VI/rKzV/NeQP82lDWrwzfGFN9KFUUpG359cIJsc+g3hiUndKW6GJK0OfyI9fJWGWEhfn3fs0YGxnKeKZPyjXrjEg+841WQKh3IyETGGGKxYE9Cbitt/erqrb+H3ErbCtXhl7YtjaIoJeOnvyck4buRud9qUCu6tCd/o7RWsBDEsQQyiMMPtdI268D70n2+Ua+Is02FBN3xCqAn1ZP5TgFdpeNK24BVPjGJ5Rx+PSRtFUWpAYlYouSyTKdtDqMsvLJwW+QUZMwoC69S6ZRv1FvOlbbWE1F3qrssZZlBK5NyDr/OumUqilIm/Jx1PBZwI3OCSTq5eSNU6Vi4Sjr4O7dIEb7xj/D9OnlaBF14BRlNPch3snCSlfw2MLcQkYpr+OrwFaXOiUu8LN0ygbzkb5BI1w23aDqI/BHmJmNJNH3pPt+o16+Tp0UYh29VNgWN8EWK2yMHTY7bI/yad8tUFKUy+G2BFyZp6+ZXnDT8KFscWng6fJ8EZyB5o2Clrb2c1A1L8vF1+AGS1U2xjMPvTnWHKl+NS7xoXUWQhDNkn5qy56qkoyiDlCBVOuv2rPOM8v0i/CGJTGWOPcKP0kvHwjNpW4YIv7Ar57bubYErXYIkbkNJOiGStiJSdHM2JtiN1X5T0QhfUQYpfo3PrIqRBW8scB/DJw/g5PAhfD98+3lRk7ahWitkv9fGjo3s7t3tdYpvYze7jUF66UD4pK1bWWaYG6sgkf9e/FCHryh1glvUefnhlwOwo3tH5DFakxlJp6O3f5vDIBUtbpSStA2awIT+m2Fbso0xQ8Z4jxsiaeuHJen0pHtCOWynjp1Bk7bWMZWK7kEdvqLUDW5OxVrW35PucT03iqQTpbWC3VZXDb8M0WlhV840ad/2EuVM2lpVTT2pnsAOGzI3HcekbZCW0Fl3XCn9PjOHoig1xU/Dt+vJvmP41OHbNzKP2loBnJOTuf1my6DhF9bhp9P+Tjcn6fho+EEkmqgafixW/OQT+PzsZamUnAPq8BWl5vhJDHGJI0hJEX5MYjTFmnLtha1zIq+0dUhOWpFtORZe5Y4x/WMHdvjp0pO2uZW2WUmnlAg/SGsFsEk6FWqrAOrwFaXm+Dl8EaEp3t+f3WsML+mgJ93Dpo7+bQ5L6qVDsVadi/DL0FqhsOImiPwUOMInwJ62WYffm+4N5/BddrwKk6gO2rcnCurwFaXG+MkxkHFAgSQdD082cdhEkvFk7udSyjKdkraB934NOKUgeXX4vg6fYFU6QbYstK5Tb6o31I0xHnOo0gmbtNUIX1EGP15OKBlLekf4AW4aw5uG541RysIrtzbAUB5Jp3COMB0nfat0AjjwPEknZNI2avM06+9fHb6iDGKCyDHJeDK3IYfXGF4UNmEruT0yJs+5ljNpWzhOKA2/DFU69gi/1G6ZQTdBt/4uKpm0rVz9TwEi8lvgXGCjMWZ6teZVlLon1xXBO8J/Y8cbRe8v3bCUt3e/zfObnvcdw3L4O3t28vKWl0tungZwz+p7aEm0cPCIg1m6YWnGhjI0T4NMw7RH2x/lC0d9IVQXzr+v+TtTR07l6LFH533+xNtPsLVrKx29Hf4rbbN1+L3pXtbvWc+0kdMC2Wx3+B29HTy69lF29uwMJulQeUmnag6fzI5YP0d3xVKUPIJE53t69zhGjpc/eHlOpmlNtHruqpSIJejs6+S6Z6/jvtfvozXRGjnCH9U6CoD/fPI/XT9zI8ycr257lVe3vxoowrfm/fGSHxOXOM985JlcOeqmjk189uHP5o7dp2Ufz7Hskk5PqoctnVsC2Wtvj3zf6vv43rPfy7PNkyqUZVbN4RtjHhORSdWaT1EGCkEc/on7ncjT657Oey9t0vSme7lk6iV8ZOpHGN40PLeblRNWhP/UuqeATNuAqBH+3IPncsL4E3hr11tc/mBmJfBhIw9j3qnzGDt0rOe51pyWU3Xj+pOv50v/+BIdvR2BNPyTJpzEw3Mf5r9X/Te3vHALveleWsjuA5BtC/2lWV/iPQe8h/FDx3uOFY/FMxuSpHpJxBIcss8hnsdb2JPZ1iK3v7z/Lxy010GBzoXKtUaG6kb4iqJ44BX5jm4dzZauLezs2cnwpuFAv2a+d8ve7D9sf9/xrcSvfUVqKSttx7eNz8sJtMRbfJ09BHdse7fsDYRb7brv0H1z0bv9ici6ViNbRwa6VpCRdfrSfcQkFtgJ2yN86/8Th00M1i0zK+k01MIrEblCRJaIyJJNmzb5n6AoA5wgJZUjW0cCsHbX2tx7VhlkUAcRZOessNhtDirVBD0urxY+wEpbC+s4+wIoy/mH0ceTsWT/wquArtLePC3snNVI2tadwzfG3GyMmW2MmT169Oham6MoFcdvlSxk5BKAPT17cu8FLYO0sBy+fZ6SnYvYXwZ0+AGPsyplelI9gTR8C6e++GGvlTV/rkonQilp0JXHuXO1LFNRBj9BNPyhyaEA7OrdlXsv58QC/honY8kihx9Vw3c6v5IRvsEEdoTWSlV76wfrWoVJGFsSWBiHb4/wwzr8QdUtU0T+BDwNHCoi7SJyWbXmVpS6xirL9HBGw5KZjpl7eosj/KBOzEnSiVql43R+uSN8e2lkKp0K/j2z3SbzJJ2Q8hf0SzphuormRfjZjdfD3ggHRVmmMeaias2lKAMRL0fY1tQGwK4ehwg/qKQjiaLNwMsa4VdQ0gnagAycNfygi8IK57cWu4XZxNyaN2xSPNdLp5E0fEVpNEz/yitX2pIZh2+P8IM2K7OwNGk7pToXyRfxAxF0zrwIP4SG77TiNoqk0xRrojvVnRkzxErbXP8fkwoVrVcjwleHryg1JkiVTjKeJC7xvE3Ic83KAnpaK8KPoru7EUXSCXpYVA3fOq4cko79BhuEmMRYt2cdKzevDHWTAttK28Gg4SuK4kyQpC1kNuWw99MJW/aXiCWKGrCVKulEGSuwhp/dhMRawBSmYyUUJG3T0ap0rHYRXvsJ2zlwrwMB+NlzP8vIUCGct/X97C2sy406fEWpMUHKMiGz7Z4lMUD0pG1eVF5qhB/hacFyun5dLQsdftgIP0/SiRDhW5ISwPbu7YHOuXT6pRw99mg6+zpzSdugWNdv+qjKtRpTh68odYKfw2yKN+XtehWlDh/yd4QqOWlbwSqdRCxBQhI5h1/KwqsoSdtE3JbcDnFjHJIYQneqO7KkU8nWCurwFaXGBNHwoTjCD520tWniFrVI2oahOdEc2uF7LbwKW4dvEeZGYf09haksAt0ARVEaiiAOP0/DD5u0zUaOfaa/Fr+cEX6YShYIlrtojjfT0ZvZeD2wpJPVze1rDqKstLVLOmEuU0uiJRfhR0nAaoSvKIOYMElbJw0/rKRjd4S10PDD3GSa48109GUcfillmdY1Dpu0jUJzvJmuvq7QGn41umWqw1eUGpNLXvr4wXJJOvbqlZpo+CFuMs3x6JKOXcMP2+YACiSdEOdZf09hNXztpaMoDYTvLkwuZZlhI3y7IyxrhF8Bh9+SaImctC21eZpd0gmj4bcmWunq64qs4WuEryiDmFLLMkty+GVsrRB83VWI1a7xptAafu57pourdMLMbZd0wj6V9Jk+ulPddbfSVjdAUZQaE2QTc+hP2r64+UX+/Oqf2dG9I9B5Fk47TNWyeZr95uVGS7wlt71g2Aj/+qXXM2vsLHb17OIzD38m1BjgvyOXG9a2iovWLuLgvQ8OfJ51g9IdrxRlMBMsZ5tL2s5/dT5/eeUvjGodxYS2CYG2zwNnhzx9ZPkW+QR1+FYjuCBE0fAntE0A4JVtr/Di5hd5bftruc/Gt3lvbWjH7vCvmnVV4PNa4i2513YJzo+VW1YC+VJUuVFJR1FqTFhJpzfdy5ghY1j4oYUsuGBB4CjS2n5wZMvI3HunTjw1otUZoqzatTtEP6I4/JGtI7nznDsB2NG9I8+BBtpMPItd0jlm7DGBz7MifIALD7kw8HnW1ovTRk4LfE5Y1OErSp3gu9I21pTb/SnKY791TmE/nVKIkrQNVfFiW3gVRtu2toT8R/s/Ap9TiD3CDzN3c6I59zrK35M2T1OUQYxfTxkLK8LvS/dFSuxZ55RzX9tydt50wv40EGZ862mmFHnEXqUT5nq3xltzr6PkAbQfvqIMYoJKOlZZZl+6r6TIsawRfgQnH6onfTya0wU4eO+DAzc9c8Iu6YR9KsmNoQ5fURQ7Yap0UiZFT6qnpCX7tZZ0wpRG2iP8sI5wRPMIdnTvCLySuZA8SSfE9bbbHEl6E63DV5RBSxhJB2BP355ITqESjiSvl06IfV+DYo+Ww0b4ezXtVZS0DUPUpwt70jZKewaN8BWlAQgi6QB09HbUTTKw0hG+dZOD8PLRXs17sb17e2SHHzlpa7M5KeEdviZtFaUBCFKWCdDR1xEpaVvpCL8SK22jOl3ol3TK4fDD3GxaE7akbYQIX3vpKMogJsgm5tAf4b+x441IUWC9RPhhujnYn2SiaPh9po9dPbtCnWdhl3TCYI/wwzyJafM0RWkAgvZ5sTuPMw44I/Q8ezfvHfocPyrdHtkeZdsXjAVhr+a9ANjWtS3UeRZRHX6ehh+iSidKC+ewqMNXlBoTtIrEqi0f0zqGj0z9SOh54rE4+w7ZN/R5nuQpOpV1+PsODWf7iOYRAGzrjubww95gLPI0fC3LVBTFTi7C94mQrUf9UhY4lbsxV6XbI9vtDZuDGNEyAoCtXVtDnWcRVQKzO+woDl/bIytKA+DnMC0HVIrGW259uOKSTjxaLTyULumUI8kdRcOvZISv3TIbkFVbV/H9Z79Pn+mjOd7Mt4//NhOHT6y1WQ1L0JW2lgMqJcIvd+K20ittS63SAXh9x+uhzrMoR6QdRcPXpK1SVpZsWMJzG58D4J/r/8mKTStqbJEC/o7QckClRID1IOmEwb5qNaw8UmqSuhyOd8KwCSWPUU7U4TcgXX1dAPzXCf8FhOvZrZSfoCtty7EjUrlr8aO0Rw6DPckc9ulERDhz0pmR5y7H01CYSh/rhlmJG6dFVR2+iJwlIqtE5DUR+Vo151b6sXYaGtY0DICuVFctzWl4glbpBE3uelFJuaASjsreWiHKk429YiYsUXe8qmeq5vBFJA78ApgDHAZcJCKHVWt+pZ/uVDfJWDK3IjDIVnNK5fFzmLk9bEv4ta3E4ivLEVci2Rhms5Ryn1/Jm6MT1o0/arO3IEjQx8mSJxI5Hvi2MebM7M9fBzDG/MDtnNmzZ5slS5aEnuuG//1D0lKd76UoilJukibGF773lUjnishSY8xsp8+qWaWzH7DW9nM7cGzhQSJyBXAFwMSJ0SpHWvtipCsngw0K7JdHb431QbX+yVp/3+Wcz5R5vMKxKWH8SnzfShOv0La21XT4Tte7yNcYY24GboZMhB9lok9f9+UopymKogxqqpm0bQf2t/08AVhXxfkVRVEammo6/H8CB4vIZBFpAv4duLeK8yuKojQ0VZN0jDF9IvJ54O9AHPitMebFas2vKIrS6FS1tYIx5n7g/mrOqSiKomTQlbaKoigNgjp8RVGUBkEdvqIoSoOgDl9RFKVBqFprhSiIyCbgzYinjwI2l9GcSjAQbAS1s9wMBDsHgo2gdjpxgDFmtNMHde3wS0FElrj1k6gXBoKNoHaWm4Fg50CwEdTOsKikoyiK0iCow1cURWkQBrPDv7nWBgRgINgIame5GQh2DgQbQe0MxaDV8BVFUZR8BnOEryiKothQh68oitIgDDqHX28bpYvIGhF5QUSWi8iS7Hv7iMhDIvJq9v97247/etb2VSJyZgXt+q2IbBSRlbb3QtslIrOy3+81EfmZlLLDdjAbvy0ib2ev53IRObuWNmbH319EFonISyLyoohcmX2/bq6nh411dT1FpEVEFovIiqyd12bfr5tr6WNnXV3PIowxg+YPmbbLq4EDgSZgBXBYjW1aA4wqeO9HwNeyr78G/DD7+rCszc3A5Ox3iVfIrpOAmcDKUuwCFgPHk9nRbAEwp8I2fhu42uHYmtiYHX8cMDP7ehjwStaeurmeHjbW1fXMjtmWfZ0EngWOq6dr6WNnXV3Pwj+DLcI/BnjNGPO6MaYHuBP4QI1tcuIDwO+zr38PfND2/p3GmG5jzBvAa2S+U9kxxjwGbC3FLhEZBww3xjxtMv9y/2A7p1I2ulETG7N2vmOMeS77ehfwEpk9nOvmenrY6Eat/s6NMWZ39sdk9o+hjq6lj51u1Ozfp53B5vCdNkr3+kddDQzwoIgslcwG7QD7GmPegcwvIjAm+36t7Q9r137Z14XvV5rPi8jzWcnHerSvCxtFZBJwFJmIry6vZ4GNUGfXU0TiIrIc2Ag8ZIypy2vpYifU2fW0M9gcfqCN0qvMCcaYmcAc4HMicpLHsfVoP7jbVQt7fwkcBMwA3gGuz75fcxtFpA34C3CVMWan16EuNlXcVgcb6+56GmNSxpgZZPa9PkZEpnscXm921t31tDPYHH7dbZRujFmX/f9G4C4yEs2G7KMc2f9vzB5ea/vD2tWefV34fsUwxmzI/qKlgVvol7xqaqOIJMk40j8aY+Zn366r6+lkY71ez6xt24FHgbOos2vpZmc9X08YfA6/rjZKF5GhIjLMeg2cAazM2vTx7GEfB+7Jvr4X+HcRaRaRycDBZBI61SKUXdlH610icly2suBjtnMqgvVLn+U8MtezpjZmx70VeMkY8xPbR3VzPd1srLfrKSKjRWRE9nUr8F7gZeroWnrZWW/Xs4hKZYNr9Qc4m0wFwmrgmzW25UAymfkVwIuWPcBI4BHg1ez/97Gd882s7auoYLYe+BOZR85eMlHGZVHsAmaT+Ue9Gvg52dXbFbTxduAF4Hkyv0TjamljdvwTyTyGPw8sz/45u56up4eNdXU9gSOAZVl7VgLfivo7UyM76+p6Fv7R1gqKoigNwmCTdBRFURQX1OEriqI0COrwFUVRGgR1+IqiKA2COnxFUZQGQR2+0hCIyAgR+V+2n8eLyJ8rNNcHReRbLp/tzv5/tIg8UIn5FcUNdfhKozACyDl8Y8w6Y8zcCs31FeAmrwOMMZuAd0TkhArZoChFqMNXGoXrgIOyPcp/LCKTJNtnX0Q+ISJ3i8h9IvKGiHxeRL4oIstE5BkR2Sd73EEi8kC2Ed7jIjKlcBIROQToNsZszv48WUSeFpF/ish3Cw6/G7i4ot9aUWyow1caha8Bq40xM4wxX3b4fDrwETK9T/4L6DDGHAU8TWa5O2Q2ov6CMWYWcDXOUfwJwHO2n+cBvzTGHA2sLzh2CfBvEb+PooQmUWsDFKVOWGQyfeJ3icgO4L7s+y8AR2S7TL4b+B/bhkTNDuOMAzbZfj4BuCD7+nbgh7bPNgLjy2O+ovijDl9RMnTbXqdtP6fJ/J7EgO0m0w7Xi05gr4L33PqXtGSPV5SqoJKO0ijsIrO1XyRMpnf8GyJyIWS6T4rIkQ6HvgS8y/bzk2S6tkKxXn8I/d0UFaXiqMNXGgJjzBbgSRFZKSI/jjjMxcBlImJ1P3XaPvMx4Cjp132uJLPxzT8pjvxPBf4W0RZFCY12y1SUMiMi84D7jDEP+xz3GPABY8y26limNDoa4StK+fk+MMTrABEZDfxEnb1STTTCVxRFaRA0wlcURWkQ1OEriqI0COrwFUVRGgR1+IqiKA2COnxFUZQG4f8HDZ6fvFVlCzYAAAAASUVORK5CYII=\n", "text/plain": [ - "\n", - "Dimensions: (id: 18, time (y): 5)\n", - "Coordinates:\n", - " * id (id) int64 2 3 4 5 6 7 ... 6000205 6000206 6000207 6000208 6000209\n", - " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n", - "Data variables:\n", - " Mass (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n", - " Radius (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n", - " px (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.8965 0.0579 2.692\n", - " py (time (y), id) float64 0.0 0.0 0.0 0.0 ... -2.79 0.004403 -2.209\n", - " pz (time (y), id) float64 0.0 0.0 0.0 ... -0.1535 0.002265 0.06436\n", - " vx (time (y), id) float64 0.0 0.0 0.0 ... 0.02465 -6.413e-05 0.01524\n", - " vy (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.00847 0.00055 0.0203\n", - " vz (time (y), id) float64 0.0 0.0 0.0 ... -1.248e-06 0.001067" + "
    " ] }, - "execution_count": 8, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots()\n", + "swiftdiff['vmag'].sel(id=plidx).plot.line(ax=ax, x=\"time (d)\")\n", + "ax.set_ylabel(\"$|\\mathbf{v}_{swiftest} - \\mathbf{v}_{swifter}|$\")\n", + "#legend = ax.legend()\n", + "#legend.remove\n", + "fig.savefig(\"rmvs_swifter_comparison-mars_ejecta-vmag.png\", facecolor='white', transparent=False, dpi=300)" + ] + }, + { + "cell_type": "code", + "execution_count": 26, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 26, "metadata": {}, "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1iElEQVR4nO2deZwkdXnwv08fM7OzO7P3fcNyw3KtiEJACArimZhDxYMkShKP6JsQo69vjDl81eSjryaKSjQqSuTjrXgQRVxA5VpYhIXdZVnYXZa9752d2Tm6nvePquru7emjqrqqu6f7+X4+85np6ur6PV2wv6eeW1QVwzAMo/NINVsAwzAMozmYAjAMw+hQTAEYhmF0KKYADMMwOhRTAIZhGB2KKQDDMIwOxRSA0fGIyBMi8pIK771ERLY3ViLDaAymAIyOR1XPUtXVzVhbRN4qIioib2vG+kZnYwrAMJqEiEwHPgA80WxZjM7EFIDR8YjIFhG5yvt7koh8RUQOisiTwAsSXPqjwL8D+xJcwzAqYgrAME7kH4CTvZ+rgbdWO1lEHhORQxV+bqryuYuAVcDn4xTeMMIw4RSAiPyXiOwRkXUxXOs8EbnPCwI+JiJ/XPTerSKyUUTWeWtm613PmBD8EfARVT2gqs/hPqFXRFVXquq0Cj/vKPcZEUkDNwHvVlUn/q9gGMGYcAoA+ApwTUzXGgTeoqpnedf8lIhM8967FTgdOAeYBFiQrjNYADxX9HprAmu8A3hMVe9L4NqGEZgJpwBU9R7gQPExETlZRO4QkYdF5F4ROT3gtZ5S1U3e3zuAPcBs7/VP1AN4EFgU6xcxWpWdwOKi10uqnexZjwMVfiq5d34X+D0R2SUiu4AXA58Qkc/E8xUMIxiZZgsQEzcDf6Gqm0Tkhbjm9ZVhLuD5ZLuAzSXHs8CbgffEJKvR2nwT+ICIPABMBt5d7WTPegzL9UBP0evvAt8GvhThWoYRmQmvAERkCu4T1LdExD/c7b33+8A/lfnY86p6ddE15gNfA95axid7E3CPqt4bt+xGS/KPuIHZZ4EdwJeJWfmr6qHi1yIyAhxR1cNxrmMYtZCJOBBGRJYBP1LVs0WkH9ioqvMjXqsfWA18VFW/VfLePwDnA79vwTrDMNqNCRcDKEVVjwDPisgfAojLuUE+KyJdwPeAW8ps/m/DTQN8g23+hmG0IxPOAhCRbwAvAWYBu3Hztu8CPgfMB7LAbapazvVTeq034Zr4xZWY16vqoyIyhpsBctQ7/t0g1zQMw5goTDgFYBiGYcTDhHcBGYZhGNGYUFlAs2bN0mXLljVbDMMwjAnFww8/vE9VZ5cen1AKYNmyZaxZs6bZYhiGYUwoRKRsRbu5gAzDMDoUUwCGYRgdiikAwzCMDsUUgGEYRodiCsAwDKNDMQVgGIbRoZgCMAzD6FAmVB2AYRhGkvzPE7t44nmvK7cIrz53ASvmTGmuUAliCsAwDMPj/3x/HXuPDiMCqnDg2DD/8tpzmi1WYpgLyDAMw2Ms5/Dmi5fy7Edfwbz+HkbH2rtZpikAwzAMD0ch5Q0WTKeEMccUgGEYRkfgqOKPls2khZzT3rOgTAEYhmF4qELKUwBmARiGYXQQjmreBZRJCTlTAIZhGJ2Bo0rK0wApMQvAMAyjY3AUxLcA0oJjCsAwDKNDOCEGkDILwDAMo1OwGIBhGEaH4iqAQhbQr57ex+Gh0SZLlRymAAzDMDzcGICrAGZN6QLgnqf2NlOkRDEFYBiGAai67h7fBfS3V58OwFgbF4OZAjAMw8B9+oeiILD3O9e++78pAMMwDHD9/1CwAFLe7tjOqaCmAAzDMCgoAD8G4FsC/vF2xBSAYRgGbh8gKBSCpT1TIGcKwDAMo70puIBKLABzARmGYbQ3hSCw+ztvAZgCMAzDaG9KLYB8FlD77v+mAAzDMADUS/fMB4E7IAvIhsIbhtEyDI3k+MI9mxkcyeWP9WRSvO2yk+jvySa27i/W7+Zur+J3nAso4SDwvoFhvvzrZxktMjWWzuzluhcuTXRdMAVgGEYLsXbbQT515ya6MinSIjiqDI85nDG/n5efMz+xdf/1jo1s2nOUvp4Mp8zpAwquoKRjAL9Yv5vP/nIz3ZkUKRFGcw5jjvIHFy6iO5NOdG1TAIZhtAx+++VvvP1iLlw6naf3HOWqT97DaMKb8Jjj8PKz5/PZ6y7IH/MtgKRdQP53vvd9VzCnv4cv3L2Zj/50Q0OCz02PAYhIWkTWisiPmi2LYRjNpVCM5b5Oe474pDdhBZATjxWCwMmu7eTrD05MP+0IBQC8B1jfbCEMw2g+WtKPJ+M9hSc9mKV4GLxPqkEWQGkTusK6iS7rrpX8EpURkUXAK4AvNlMOwzBag9J+PIVc/GR3Q1UtNQDy6yf9IO4rmEL6qXe8ARXIzbYAPgW8D6j4X1dEbhCRNSKyZu/e9u3LbRjG+I6cDbMAKCidYlLSSBeQ+7uRLSiapgBE5JXAHlV9uNp5qnqzqq5S1VWzZ89ukHSGYTSD8TGAxvjDHdW8D76YlCQ/GH5cE7oGuZ6guRbAJcCrRWQLcBtwpYh8vYnyGIbRZLS0Gte3ABIux1UtKJ1i0g2YC1yIe3hrNij4DE1UAKr6AVVdpKrLgNcDd6nqm5olj2EYzWfcUBb/aTjhzVAVpEwUIC3SABdQSRO6BvYganYMwDAMI09pEDjjpYEmnwWk5WMAqeRdQP7VS3sQtX0WkI+qrlbVVzZbDsMwmktpTrxvASQ9mN2p4AICuOX+rYk+jZfGPfweRG3tAjIMwyilNCc+6+VE7jpyPNl10XF1AAA92RSqsGnP0eTWLnF7dVohmGEYBjDeHy4ivHJlcj2ACuuWtwD+5bXnAMkGoQt1AO7rRsU9wBSAYRgthO/3Ln4aTzfCD6+UTQNNeztkkjGIcYFvswAMw+hESv3h4CqARgSBy1cCu1tkkpXI42MApgAMw+hA8v7wopScTCNy8SnvAsrkN+ME11ZFpCjwLeYCMgyjA8k/DRcdS6dSiVsAjpYPAucL0RK1AMa7vPzjSWMKwDCMlqHUHw4NsgB0XDfo/NqQrDvGKalBMBeQYRgdSWkhGHgxgCR9MPhumGoWQLJB4OK1/e9uLiDDMDoKLWmMBu5TeNIPw5V6AfmVyLkE00BLq5AtC8gwjI6k4AIqHEunJVEfPPjtoJtlAZwYf2hkN1CbCWwYRkW+8eA21u88csKxbDrFn19+EnP6emJda93zh/nx4zuB8TGA46MO9zy1l8tOTaYlvFMhDTTjVSJ/9TdbuObsebGv+8i2g/z66f1lg8Cb9gzw4hWzYl+zGFMAhmFU5F9+9CRjjtLblQbcJ/TDQ6OsmDOFN1y0JNa1bn1gGw9tOcBZC/qZ0lPYms5ZOBWAm1Y/nZgCUD0x9dRn3lRXyd33zH5Gcw7ZdLxOk1t+s4UNu45wedH3WjKjF4DtBwdjXasc5gIyDKMiOVWuf/Ey1n7oZaz90MtYfeNLABgezcW+luMoc/t6+PFf/c4JG+01Z8/n0hWzGE2yHUMFC6C/J8vfXXM6kEw7iJzCslmT+fKfXJQ/Nre/h96uNA2IAZsCMAyjMqUZKt1Zd8sYHovfJ1+aDllMNi2MJpgJ5BaClV/cb0g3mkAcolr9QdK1D2AKwDCMKvhVqj5d6SQVQLVNOMVIAmv6lH7PYvJziROwACq1oHAzn0wBGIbRRNwq1cLrjKcABobHYl9LVfO98EvJZlLJWgAVCsGg8J2TqEVwnErZR8lXP4MpAMMwqlDORZFOCV+895mGrOWTTUniMYCKa+ddQPGv71SwPDIpSbT2wMcUgGEYZVHVsm2Sz100lUnZdOzrlfbEKSabTtgCoPJEMD8gPZqQ28tiAIZhtBxapigL4NJTZnNsJBd7oVKlp2HwXUBJVuNWjj/kXUAJBIErub3SKUm0BbWPKQDDMMpSOp3Lx3+17UC8eerV/PBdCVoA+fYTFd7PehowCQVUyfWUSQkN8ACZAjAMozzl2jIAnL9kGgD7BoZjXq+6Hz45BeD+rrR2IQichAKoNInMLADDMBrM1+7fyk2rnwaKJ1WduEHNnNwNwIFjI7GuXV0BpBgcySXSIK3cFLJiMonXAYw/nk4JP3l8V946SQpTAIZh5Pn776/jX+/YyNHjo/ljpZvyjCldQBIKoPIm7G+Dz+wdiHXN4mtXLEJLJWcBuOuOX9hXujsPH09kzfzaiV7dMIwJSc7Rsr35AWZO7kIk/s1Jq1gAFyyZDiRXgQwBKoGTqAOoYAH86SXLgORbQpsCMAxjHGOOlp3OBdCTTXPa3D4e2XYw1jUdpWIhWHfGr0COvweR72Wp7ALy0kATKgSrNojGFIBhGA2hOK2z2AIotzGeOrePrfvjzQKqFgPo8hXAaBKpmO5vqZAH5FsAyQSBK8cAwG3GlySmAAzDAE7cbMYcRb29ttymPLuvm20HBjlSFCuol2q9gAoWQAIKgPKuLh9/KlgydQDl769/LOmhMKYADMMATnQ3OFViAAA7Dw8BsPLDP4stGFw6GrEY3wLYuPtoLGsV49RwAfmzEO7bvD+BtSt3AwWzAAzDaBDFCmDsBBfQ+A1q79FCDcD+mOoBKj0NA8zrdwez9GTi37IKhWDl11460x3QEvcwGKhc/VywAGJf8sR1kr28YRgThbETYgBOxUIwgAXTJuX/jsstU20egG8BJFGN61+xkgUgIsya0s2xkQSG4FRwe/n3IemW0KYADMMAxlsAWsUCeNXKBfm/48rMcZ+GKxeCQTLFWH6so9La4LqBhkYSaoFdLQhsMQDDMBpBcZBzLFc5DRTgqjPn8pk3ng/El5lTOnugmEJHziQsgOpBYHAVQFIWQNkgcINiAE0bCi8ii4FbgHmAA9ysqp9uljyG0Uo8+OwBvrXmubLv9Xal+euXncbUSdlY1zx4rJDR86k7n6Lba/lcaWNc6LmBPnf3ZjLpFBctn1HX+g8+e4AXnzyz7HvplJBOJdMPaN+AG8Susv8zuTvD49sPx7bmuucPc8t9W3ju4CBz+7vHvZ9uUBZQ0xQAMAb8jao+IiJ9wMMi8nNVfbKJMhlGS/DfD2zl9sd2MrfvxM1hJKfsGxjm8tNmc+Xpc2Ndc/vBQl7/kzuOAG4A9Iz5/WXPXzZzMqfP6+M3m/czY3JXXQrAH/c4WOUpO5OQAtiwy/2u/VUU6sDxMbKZaioiHN95ZDvfeng78/t7eOHy8UqvUS6gpikAVd0J7PT+Pioi64GFgCkAo+NxFBZPn8Tqv73ihOMbdh3hmk/dm2hB1O3vupRzFk2tef70yV3c8d7L+N1PrK57eIm/0V1z9ryK53SlU4wkoAD8772yync+f8k07tqwJ7Y1x3LK9N4ufvOB3y37vu8W6og0UBFZBpwPPNBkUQyjJVDK+4a7M65bJomCqFyNrpiVyKRSdY8v9OMP6SqLZzOpxKpxoXoQeFJXmqHR+GIAY47mn/LL4b/X9mmgIjIF+A7wXlU9Uub9G0RkjYis2bt3b+MFNIwm4FSYjpJsTxx3I6y2MZUjjvGF/kZXbe2kZgLUmgcAMCmbZijGIHDOcchUVQDeee1sAYhIFnfzv1VVv1vuHFW9WVVXqeqq2bNnN1ZAw2gWFbJD/Hz4ar7yqOSqtH6oRiZd//AS3wLwe++XXSeVSuR7V6t49untSjPmaGwKqJYF0PatIMS1t74ErFfVTzZLDsNoRRzVslkpfluCHzy6I/Y1c3kLINznUlK/BeDHAKptijlHeWJHfJk4PtXSXX16vIyouNxAuYAuoHauA7gEeDNwpYg86v1c20R5DKNl0ArDUXq7MkzpztDXE3/+huafhENaACmpu2LVVyDV3CLd2RSzpoxPmayXWhPBwI0BALG5gYJaAG1bB6Cqv6J66q1hdCzVWiOfvbA/mSCwE00BpFNSd3C2YAFUfiZdMHVSQh05a3/v3pgVgONoVWXny2IjIQ2jA6n2z747k05oMpb7O2wQ2I0BxGMBVHM/pVNCEh6RIC6gSVn3WTmuGIRrAVT+sgUXUCzLVcQUgGG0INU6Y3ZnUgzHmJLo4wccw6aBplOpGGMAlbckkWR84kGCwFO6XQVwLKZ+QLkaFkBHZAEZhlEerdAmGNxMoA27jsaeIZKLmAaaSQlb9h/jQz9YF9lFcnhoJH+tSqRTkohLpDAPoPLak7tdF9DA8foVwIFjI9y1YU++3085fOX/3Ue2171eNUwBGEYLUqkQDAqb5POHhmJd038SrlaMVY5LVswim05xy31bI2fpDAzXVhwpkUSeiDWEBTAwXL8CePDZAwAsmdFb8Zx5U935B5t2D9S9XjVMARhGC1JpUAjAlWe4PYDiLgYruIDCKYA/u3Q5n379eQCRXUH+JuxvfOVIiSRSGRukEGxKT3wKwHdjvfvKFRXP6e3K8LoLFtW9Vi1MARhGC6JNmI8bNQgMhbm5UX30QTbhlCQzIMUJkAU02Y8BxKAA8m0vatzndKq96wAMw6hApUIwSE4BFNJAw3/WD1pGtQCCrJ2Ood6gHPkYQJXdcHKXqwCOxhADyAWoeQD3+1oQ2DA6lEr7g98Q7sjQaPkTIpJ/Eo6gAdJ5CyCaUgryFJ6S+tNNyxGkDiCdEncoTCwWQLBgu+vyMgVgGB1HtfGIfhXwR368PvY1IXwhGBSeZqMWhAVSAInVAQSzfCZ3Z2JJAw3S9sJ/3ywAw+hAKjQDBeCsBf1M782GztevhV90FDYLCOrvXRMk/pBcDMD9LTUaE0zpzsTiAgpjAVgMwDA6kGqtIESEK06bE3tnzIILKPxnfQsg6hNroBhAQhtikF5A4NYCxHHPnXwMoPqNduse6l6uKqYADKMFUaVqp6zemDajYhwnugsoVbcFUDv+kEpoQwySgQTQm80wGIMLKLgFYFlAhtGRuIVgld/v7YpnMyom74ZpQgwgaBpoIhZAwOynuJSuHyivlQWUakAMoJlD4Q3DqIDbCqLy81lvV5rjow7HhsfyOer1EnUkJBSeZm+5fyurn9rLK1fO5+qzKs/3Hbd2wDTQXUeO8/DWg1y4dHp4IUv4yeM7+em6XTy166i3dg0LoCvN9oPRFYCq8vE7NrJ6oztbuGYQ2LKADKMzUa3ui585uQuA+zbvj21Nx1FSEr4SGGB2XzcXLZ/B0aFRfvbELr5+/9ZwawfIArr4pJkAfCem/jhf+c0Wfv7kLkZzDi89c25NxTcpm2GwjjTQw0OjfP7uzew9OsyVp8/J13NUwrKADKNDcQvBKu9ILz3Tfbq++6n45mRXCzzXojuT5pt//iLuuvElrFw0NXS2TpAYwGvOW8iSGb2x9eTPOcqqpTO468aX8J9vWVVT8U3uTjNYRxfWES/N6n+99FT+6/oX1FwvJW7MI8mZAKYADKMFUaq7YuZN7aGvO8OR4/EVg+VUIxWBlSIRevYEjT/EVYwFtadylTKpq74YwKgXH+kKOHOzEWMhTQEYRgviVOkF5HP6/D52Hzke25qq0QLApUTJ1w9ajNVb5yZcTM5xagZii5nclWFkzGEs4pSWUa91RzYTbM10nam1QTAFYBitiGrNzXBOXw97jgzHtmTOqb1mEHzXRRiCdiKd3B1f9tNYLpwF4I+FjOoG8pvAZQNaAL47LokOqPk1kru0YRhRcapUAvtM681yOMZ+QE5sLqAoFoD7u9aG7LqA4osBZNLhXEAQfS7wyJj7JYMqgEZMBTMFYBgtiFI7INvXk42lNYGPmwUUhwsofNfOoJ1I+3qysfTk99esNoKyFL8jaNQYxKjnOgoaA/D/W1gMwDA6DMepnY/f15NhJOdwPKb5wDkN5xKphEj4pm1BO5H29WRi64I65ighDIC8BRA1BuErgKBWh//fIslagJoKQES+JCLnlRz7cFICGYbhZgHVcgL5XUHjsgKcKoPow5CS8KmLQdsx9PdkGRgZi2VTjGoBDEVUuH4aaOgYQJNdQFcDXxGRtxQde3VC8hiGgbuB1naH+AognidiJ8YgcNj9ORcwC6ivJ4MqHI1pNGOYLCDfAojqAvLbZARWAC2SBbQHuAz4QxH5rIhkqB2fMgyjDtyRkNXP6evOAvFYAAePjXDbQ8/VfR2oNw20hgUwyf/O9Sm9j/10A7uOHCcdwgc0ubu+IPCvnt4HhKgD8O5FXIVv5QgkiaoeUdVXAXuBu4GpiUlkGEbAIHB8LqC1zx0EYMWcKXVfK1IMIGAn0n7vOx8Ziv6dHUf5/N2bAciGSQPNekHgiBvy9oODACyZ2RvofF8pHhqMd/JbMUEUwAoRWQigqh8GPgpsSUwiwzC8QrDq5/T1xPM0DIVhMB94+Rl1XytKDCBoGmh/DN+52KXy4hWzAn+ukAYaTfkMjeQ4Z+FUpnpWTC0Wz3AVRdQ5y0EI0kbwk8DPROQAcBvwbVW9MjGJDMPwuoE2zgLIp2HGkBcYpRAsTBoowJE6vnNxWmWoSuDu+rKAhkZzTMqmA5/vWydRK4+DUPM/t6p+WFXPAt4JLADuFpE7E5PIMIyqIyF9+vObYf0WgO9uiCMNNEodgKvwalcC90/yXUDRv3PxE3WYwreejBcEjqwAHHq6giuAjBcrSNICCKPv9wC7gP3AnGTEMQwD/IEwNWbUJmABxNELKGolcJAU1Cne7IN6isFyuWgWQColdGdSHB4cibTu0MgYvSEsAF8ZjzbTAhCRvxSR1cAvgFnA21V1ZWISGYbhtoOusTelU8LkrnQsCiBoIVYQIrmAAqS9gjsJDaLn4kOhJw+Et3hE4MmdRyKtOzSay8cRgpBNJ18JHCQGsBR4r6o+mpgUhmGcgAZ8InbbQcQRBI7PAoiaBhrk+/pDVOpJjTwxBhAu6DGjtytvhYRlaMShJ4QF4Ms2GnHMZqA1ap2gqu9PbHWjZfjZE7vYfnAIgEtWzOK0eX1NlqizcQfC1KavJxOvCyi2GEC4zwTtQ5RKCV3pFM/sOxZRuhN96mG/74JpkxgeC++SUVX2DQznO4oGwbcAxhJsB9rUmcAicg3waSANfFFVP9ZMeTqVY8Nj/PnXH86b7ZefOpuv/ulFzRWqw1ElULllX0+Go8PxBYFjGwgTIQYQdDMeyTl1jWaMmgUE0J1NcXw0/Ib8rKewwtyXfBA4QQugac3gRCQNfBZ4OXAm8AYRObNZ8nQyh4dGUYUPvfJMLl0xi0MRg1xGfGhAl0hcHUH9PTGOVhAihJ8HECDm4XPu4mmRnsJ9cnVYAN2ZNMNj4d1PB465/6YuP3V24M9kGhAEbqYFcBHwtKo+AyAitwGvAZ6Me6FbH9jK6o3xzU5tN/wBG3P7e5jWm2XttoO8/ZY1TZaqs9l/bCSwC2jbgcHI6zx3YJCP3bGBHz+2E4gvBrBvYJhP/Gwjf/Oy0wJ9Jkwr6pS4bRV2HznO3P6eULINj+W48Vu/zb8OrwBSrHv+iCtviM/6qbpBi8Cg0DW02YVgSbEQKG4+sh14YelJInIDcAPAkiVLIi10aHA07982ynPBkmmsXDSVMcdh895jdr+azMmzp3BZgKfFeoPAv9m8L7/5QzwuoMtPncNdG/byudWbgyuAEC6gtdsOAXDr/Vv564DX99m0e4A1W922FwunTWLJjGBtGXx8t8yuI8dZMG1S4M/5g3tCKYCU7wJqTwug3H/tcapOVW8GbgZYtWpVJFX4zitW8M4rVkT5aMexeEYvrzlvYbPFMALS35Opqyq21JUShwXwipXz2bDrCDet3hz4M0HTQMGtBRgYHotUkOUXkH3j7RfzopNnhv78FafN5vbf7mAkpAvq8GB4BeAHgZPMAmrmQJjtwOKi14uAHU2SxTAmJH097qDyKH5pYNxGFocFAOEDwUFjHgDzp7punzDjHH1ufXAbUKgoDku3Vw08EvKp/LDXvK4/lAvI3Z7bdSLYQ8ApIrJcRLqA1wM/bKI8hjHhKDSEi2YFjLMAYlIAKS8QHLQpXC5EDOBzb7oQgHkh/f8Auw8fB+D0ef2hPwvQ5dUhDIfMBDo8NMrkrnTgWQBQFAROMA20aQpAVceAdwH/A6wHvqmqTzRLHsOYiPgN4TbsPBrp86UZX3G4gKBQxBbUCHBbQQQ7d8E0d+OP0hJ624FB/mjVosiKzi9EC2txHR4aDfX0DwUF0JZpoACq+hNVPVVVT1bVjzRTFsOYiMya0g3Am770QD7VMAwPbjl4wus4uoFCYTMP6gZyNHhWje+GeeDZ/aFkGhrJsefocOjAbzH+pvz0noFQnztyfDSU/x8K1lhTu4EahtG6XLJiFn+0ahEw/mk+CFO6T6xMjcsC8Lt6BnVfh0kDTaeE/p5M3voJij+QZXEdCmDhdDfzJ4wrB6JZACJCNi2MtmkMwDCMOkmnhCtOc5vzRimOKvVlxxUDkNAWQLi1F03vDR0c3brfVQBLZ04O9bli/BhA2OKsI0PhLQBwU0HNAjAMoyLdWd8vHUEBjDkntEOo1Y8/KGFjALkQlcDgpkiGTY/cst9tx1CfCyiaAjgcVQFE+J5hMAVgGBOcfGpiJAUQrkVxUMLGAFQ1lPspk06FbpL27Ye305VJMb03/Ebs4w90Hwm5KUdVANl0qm3TQA3DiAE/M+UjP1nPn3z5QW7/bfBymqd2D4QaUxiUVD4GEGzz2rDzaCgLIJMK92T8lV8/y4ZdRzllzpS6rJxsxv3sdx/ZHvgzozmHwZFcJAWQEvjF+t2hPxf4+old2TCMhnDKnD4uWTETVeX+Zw7w7YeDb04i0VxHta8bLgg8rTfL/oHgQexsOpxv/Jb7tgLw7itPCfyZSusCPLEj+FCYIxHaQPiM5nR8e4QYMQVgGBOcqb1Zbn3bxfzwXZdyxvy+wE/dOUdRhZWLpgKwYGr4wqpK+C6goIVgjsJZC6cGvn4mLaGapI06Dr9//kKuOXte4M+UXTdCkDxKHyCfV6ycbzEAwzCC4WaNBNswjntjFSd7YxbjCgBDsQso2PlOiF5A4H7PMBvj6JiGTt0sR5R75CuAKO0nutKpyG0+gmAKwDDaiHRKAgcN8wog4ojDakQpBAsTBM6mJZQLaDTnROodFAd+s74oFkB3NpWIi87HFIBhtBGuayTYhjGUVwBuEDhGA6AoBhDUHRXu6TqMogO3eVscFkBYVJW129xqa79vUxi6M2lGxpzArrSwmAIwjDYiJUJQz4g/2rA37wKKVw4IXgegIV1A2XQqVJO00ZyTL+JqJBt3H+VTd24CYObkrtCf91N7o7T5CIIpAMNoIzIpIRdwY/RdQH47CAk0gywYkVxAITRAJiWhmqSN5TTfX79eXrh8RuDN/OAx1///wWvPYKbXtykMftHaseFk4gCmAAyjjUiH2Bh9BZCkBRDUSxOmHTS4hWBBg8COo4w58QSBwW1D0ROwdsJ3x523ZFqktaZ4/Y7Czh8IiikAw2gjMukwQWB3U8nHAGKUI98LKKAsquGG0WRDxDp8V1FcCiATIv7gK+Mo6aMAXfmpYKYADMOoQToVvHWAbwH4rSSSSAMN0wsobBpoUEvHtxTicgGlUsFrEPzzoiqfqL2HAl8/kasahtEU0gI7Dg/x5i89QCYl3Hj1aZy1oHyB1VBeAbibTNh2xdXwdcnffOtRXnHOfK6/ZHnV86OkgQbZFP/x9ifyw3LitAD2DQyjqjWVpp+qGrXLajZi99GgmAVgGG3Ey86axxnz+zkyNMovN+7l3k37Kp7rWwCnzuvj2nPm8eFXnRmbHOcvmc7FJ83gqd0D3PbQczXPd0KmgQapBHYc5cu/3sKW/ce4aNkMLlo+I/D1qzEw7Ob1B7ECRp36rI+kB8ObBWAYbcS158zn2nPmc3w0x+l/f0dVF8xxL8WwryfDTdddGKscy2dN5rYbXsQ7bn2YTbtrT88KWwnsu7qqPYX7G/SbLl7KO69YEfziNThl7hTADVzXigX7GVmZiKPW/O6jZgEYhhGYIN04120/DBSygJIgaEwibBpo1h+XWOXa/rpxDbnx8QO6Qb6X/+Qe2QVkCsAwjLAEacbmP5ROSaAVhE8mYMA0bCVwxtsYqwWCx/JP3/EqgLR344J8r7FcfUFg/3MjY1YJbBhGQILk4e8fGOG0uX2JyhG0ZYOqEmaPzPvGq6SCtoIFkHcB1R0DMAvAMIyABJnJe3BwhOmT48v8KYdrAdTevNw00OCbpN8gbd/R4Yrn+E/o8VsAvvup9vcarbMOoGABmAIwDCMgIoJIDQvg2AgzJ4dvTxCGoBaAE7ISeJk32N1PZS1HwQKId5sLYwHk3VARXUD+vGerBDYMIxRC9RjA4cFRptUxHzcIQWMAjhJKAfjVy9VaJSduAQRIzaxXhnrmPQfBFIBhtCkpkYouIFXlUMRB5WFIp1LkAmyUYdNA/Y1xeLRKDKDODJxK+P78QBZAnS4gv0gvqaEwpgAMo01xFUD5946N5Mg5mrwFEHB0Y9g00K4AG+NYnQHYSoTLAqqvEjivAKoounowBWAYbYobAyi/SR0adPvLJ20BpEQYGs3x9J7qxWBhK4H9jfFjP93AG26+n288uG3cOY9sO5SXIU78lhXPHRysea7bhVQi91nKpFOkU5LYVDBTAIbRpqREKlYCFwaVhx9SEoazF/YD8MCz+6ueF9YFdPLsKVx1xlz6e7Kse/4w31v7/Lhztu4/BsB5i6cFv3AA/MypI949rMaYE86yKUd3Jrm5wNYKwjDalJRUbsd8eNBXAMlaABefNBOo3RY6F9IFNKkrzRffugqA6754f1kXSc57+l7sDVWJiwVTJ+WvX4vRnEO2ziwkVwGYBWAYRgiqxQB8C6ARWUBQ3V+uqqhGb0edTqXKXj8Xw9N3+fXCFIJp3TGI7kzaYgCGYYSjagxgqDEWQJDN0n8rTDvoYioNaBlzNHITtmr4g2uCjLsczWnddQjd2eRcQKYADKNNSaWkYh3AocFGWQC1M2b8jTTqw3q6Qq1BYhaA+Eqt9rljOafuQTRt5wISkX8TkQ0i8piIfE9EpjVDDsNoZ6q5gPYPDDMpm060EygEtQA8BRBxs06L5HvuFDPmOLEXgUGhiV4ugAUQmwuonRQA8HPgbFVdCTwFfKBJchhG25Kq4gLaNzDMrL5kM4CgKAZQpRjM37ujpmumK8xBTtoCCDLveDQGN1SSWUBNUQCq+jNVHfNe3g8saoYchtHOSDULoAF9gMB9qheh7BO6T70uoIoxgJwmYgGECQKP5eq3QrqzqbYOAv8p8NNmC2EY7UZKKvcCOjg4wvSE/f8+tfoB+a6UqE/rVWMAMVcBQ7ggcDx1AMm5gBJzAIrIncC8Mm99UFV/4J3zQWAMuLXKdW4AbgBYsmRJApIaRntSrRfQ0EiO3gQHwRRTqyOoentb5GrZBmcBFYLAwSyAeofRT0gXkKpepapnl/nxN/+3Aq8ErtMqLQtV9WZVXaWqq2bPnp2UuIbRlnxzzfayVsDQSI5JtQbaxoQqfOGeZ3j+0FDZ9/MWQB1ZQDsPH2fNlgMAPLTlAH/8hfu4d9PeyG6lWutB7SDwjx7bwS837o2pEriNXEAicg3wd8CrVbV2Qw3DMELT2+Vu8IMj458eh0YbpwB+55RZADy+/VDZ9+vNArr6LNfRcO+mfQDc89ReHnj2AGfM7+ePX7A40jWr4Qera3mA7li3C4A/uLC+EGc7FoJ9BugDfi4ij4rI55skh2G0La9/gesyLbdPDY7k8goiad53zelA5VoAP5smqgvoJafNAQrfU9V9Sv/vt1/MDZedHOma1fD1VC0XkKPKijlTeNPFS+taL8lCsKb0AlLVFc1Y1zA6iUpjIR1HGR5z6GmQBVAra6beSmA4MeAdtrFcWIJmAeUcres7+XSl28wFZBhG8uRdFSV7x3HvabJRFkCtEYr1poG6ny0EvJ06+goFoTBus5YCiO7WKsa1AEwBGIYRglQFC2D/gDsLYFKDFEBhiHr5DdNXDPVslsVVz5qwBQB+9XFtF1CdCUCAGwPIOZofLhMnpgAMo02plK/+tfu3AjSkEAwK/YAqbZi+ePUMbil+InddQMlqgFRKamYBxeUCKoyFNAVgGEZAfDdI6b47POq6gF5+drkynfipaQHE5ALy9+OwA+ajkBap2QrCUY3HBWQKwDCMsPh7T2kdwEjOYU5fdyybUxDyMYAKLgynzkpgOHH4jaNKwvu/V9xW/ZzYLAAvWJ9EJpApAMNoU1IVLQCH7mzj/un77RiSSgOF0hgAJK3aqjXa88k5MVsACdQCmAIwjDalUhB4eMyhO9OYADAEyQJyf9fztDwuBpCwdVOrvYUvRzwxAN8CMAVgGEZACjEAd6NSVd5z21p+9fQ+uuJITwmI79r5zF1Pj3vv2w9v5z23rQXqjAEUDb9pRBAY3GB6pS42X79/Kw9tOci2A/U3OijEAMwFZBhGQPwt0N+jjo3k+MGjO5g5uYs3vLBxjRW70im6MymODo+Ne++OdTvZfnCIl545l/OXTI+8RrELyA0CR75UIPp63E6qoxXmHHxu9WaAiv2PwuC768wCMAwjMKkSCyDnbVbXXbyUN9fZniAMIsJfXO62ZCh9Ys45ykmzJ/Ofb1nFvKk9kdco9smraqKFYABvuMhVoJXiAEFaRQcl7wKyGIBhGEHxOyH7T8Z+umUSQ1JqUSkOEEe/fDhx+I3jJG8BZGqktsa5/LTeLOctnpZI4V5TegEZhpE8pRbAmDeVK4kxibUozgQqjj/nnHimdo3vBZR8EBgKVlWSnDq3j++/85JErm0WgGG0Kb4bxN8Y/afvdrQASnsBJa0AMnmFlkyPnkZhCsAw2pRCIZj72x/M3hQLwPNHlbpMcjFN7SquBFYaUwgGwaaCtTKmAAyjTSktBMtbAAnMya1FNQsgjpx9tw7A/Vsb1AoCKscAJgqmAAyjTSktBPM3q3QCc3JrUegHdKLLxIktBlBaB1D3JatiFoBhGC1NaSFYe8cAaEoMwBSAYRgtSens2qZmAfkWQK40BuDEZgEUCsEaEQMoH9PwmShqwRSAYbQp/r76ntvW8qr/+BV/9Q235UIc/WnC4j8x/+WtD/PkjiMADAyP8dTugdhiAKs37uFV//Er7n1qb+KFYL7SuuFraxgtaQv6/bXPs/PwcQAmNWjsZlSsDsAw2pRzF0/jFSvnMzTi95Dp5vR5/Zy3ZFrDZXnBshlcdcZc7ly/m4e2HODMBf1s2XcMgNlT6h9M8+aLl3LPpn3u9fq6ufzU2XVfsxoXLnXbVjyz9xj7B0ZOqGK+c/1uAM5dNJW/vfr0ROWoF1MAhtGmzJrSzWffeEGzxQBg0fRePvPG8zn97+/g2IjbE8h3n1x26qy6r3/9Jcu5/pLldV8nKHP7e/j4687h777z+LjJYKpw8uzJ/OBdlzZMnqiYC8gwjIbQnUkhAoPDrkWSy8ckJuY2lA+yl8QBlMZ0I42DiXnnDcOYcIgIk7syDI74CsA93oyspDhIl2RZ+bi9iCbGdzIFYBhGw+jJptmwyw0CNzMrKQ4q1QI0IgspLkwBGIbRMMYch7XbDgHNrUuIAz97aZwF0IA6hLgwBWAYRsO4+sx5+Sd/PwjcqOH0ceO7gEqHw6sqEyWsMUHENAyjHVg8YxKjOWX7wcF8K+WJagH4UzXLuYDMAjAMwyjhrIVTAfjeI88X9SaaGJtlKaXzFnwcJfFCtLgwBWAYRsO44rQ5LJvZyxM7jhTFACbmNlQtCDxRdJoVghmG0VDOXTyNHzy6gwe3HAAKrpSJhh+7KFcIZi4gwzCMMqyYPQWAA8dGeP0LFrN05uQmSxSNdIVCMLMADMMwKnCB10cH4GOvW9lESeqjeh3AxNAAZgEYhtFQpvd2NVuEWPD3+FIXkFsH0ASBItBUBSAiN4qIikj93aAMw5gQnDG/r9kixEK6ZN6Cj1oaaG1EZDHwUmBbs2QwDKPxTBT3SC0qu4AsCByE/we8j4kzPMcwjJhICfT1TOwQZKUsoInUC6gp/wVE5NXA86r621pPAyJyA3ADwJIlSxognWEYSfPkP13TbBHqplIW0ERKA01MAYjIncC8Mm99EPjfwMuCXEdVbwZuBli1apVZC4bRBvS0+KjEIFRyAamlgYKqXlXuuIicAywH/Kf/RcAjInKRqu5KSh7DMIw4qdYKouMtgEqo6uPAHP+1iGwBVqnqvkbLYhiGEZWCBXDi8YlUBzCxozCGYRhNwm9h8c8/epJP3flU/vjWA4MsmDapSVKFo+kKQFWXNVsGwzCMsCydOZnrXriEg4MjJxw/Ze4UXnfBoiZJFY6mKwDDMIyJSDad4iO/d06zxagLawVhGIbRoZgCMAzD6FBMARiGYXQopgAMwzA6FFMAhmEYHYopAMMwjA7FFIBhGEaHYgrAMAyjQxEtHWfTwojIXmBrxI/PAiZCvyGTMz4mgoxgcsbJRJARGi/nUlWdXXpwQimAehCRNaq6qtly1MLkjI+JICOYnHEyEWSE1pHTXECGYRgdiikAwzCMDqWTFMDNzRYgICZnfEwEGcHkjJOJICO0iJwdEwMwDMMwTqSTLADDMAyjCFMAhmEYHUpHKAARuUZENorI0yLy/ibLskVEHheRR0VkjXdshoj8XEQ2eb+nF53/AU/ujSJydYJy/ZeI7BGRdUXHQsslIhd63+9pEfl3iXk4agU5Pywiz3v39FERubaZcorIYhH5pYisF5EnROQ93vGWup9V5GyZ+ykiPSLyoIj81pPxH73jrXYvK8nZMveyLKra1j9AGtgMnAR0Ab8FzmyiPFuAWSXH/hV4v/f3+4GPe3+f6cnbDSz3vkc6IbkuAy4A1tUjF/Ag8CJAgJ8CL2+AnB8GbixzblPkBOYDF3h/9wFPebK01P2sImfL3E/velO8v7PAA8DFLXgvK8nZMvey3E8nWAAXAU+r6jOqOgLcBrymyTKV8hrgq97fXwVeW3T8NlUdVtVngadxv0/sqOo9wIF65BKR+UC/qt6n7v/JtxR9Jkk5K9EUOVV1p6o+4v19FFgPLKTF7mcVOSvRcDnVZcB7mfV+lNa7l5XkrETT/g0V0wkKYCHwXNHr7VT/nzxpFPiZiDwsIjd4x+aq6k5w/1ECc7zjzZY9rFwLvb9LjzeCd4nIY56LyHcHNF1OEVkGnI/7RNiy97NETmih+ykiaRF5FNgD/FxVW/JeVpATWuheltIJCqCc/6yZua+XqOoFwMuBd4rIZVXObTXZfSrJ1Sx5PwecDJwH7AQ+4R1vqpwiMgX4DvBeVT1S7dQK8jRLzpa6n6qaU9XzgEW4T8lnVzm9afeygpwtdS9L6QQFsB1YXPR6EbCjSbKgqju833uA7+G6dHZ7ph/e7z3e6c2WPaxc272/S48niqru9v7xOcB/UnCTNU1OEcnibqq3qup3vcMtdz/LydmK99OT6xCwGriGFryX5eRs1Xvp0wkK4CHgFBFZLiJdwOuBHzZDEBGZLCJ9/t/Ay4B1njxv9U57K/AD7+8fAq8XkW4RWQ6cghsgahSh5PJM8aMicrGXufCWos8khr8RePwe7j1tmpzeNb8ErFfVTxa91VL3s5KcrXQ/RWS2iEzz/p4EXAVsoPXuZVk5W+leliWp6HIr/QDX4mY4bAY+2EQ5TsKN/P8WeMKXBZgJ/ALY5P2eUfSZD3pybyTBbADgG7gm6ijuU8ifRZELWIX7P/lm4DN41eYJy/k14HHgMdx/WPObKSdwKa7Z/hjwqPdzbavdzypytsz9BFYCaz1Z1gEfivpvJuF7WUnOlrmX5X6sFYRhGEaH0gkuIMMwDKMMpgAMwzA6FFMAhmEYHYopAMMwjA7FFIBhGEaHYgrA6EhEZJqIvKPo9QIR+XZCa71WRD5U4b0B7/dsEbkjifUNoxKmAIxOZRqQVwCqukNV/yChtd4H3FTtBFXdC+wUkUsSksEwxmEKwOhUPgac7PVo/zcRWSbejAERuV5Evi8it4vIsyLyLhH5axFZKyL3i8gM77yTReQOr7HfvSJyeukiInIqMKyq+7zXy0XkPhF5SET+ueT07wPXJfqtDaMIUwBGp/J+YLOqnqeqf1vm/bOBN+L2bvkIMKiq5wP34ZbngzvY+92qeiFwI+Wf8i8BHil6/Wngc6r6AmBXyblrgN+J+H0MIzSZZgtgGC3KL9XtkX9URA4Dt3vHHwdWeh00Xwx8q2hgU3eZ68wH9ha9vgR4nff314CPF723B1gQj/iGURtTAIZRnuGiv52i1w7uv5sUcEjd9r/VGAKmlhyr1H+lxzvfMBqCuYCMTuUo7hjESKjbN/9ZEflDcDtrisi5ZU5dD6woev1r3I60MN7ffyqFbpGGkTimAIyORFX3A78WkXUi8m8RL3Md8Gci4nd3LTdq9B7gfCn4id6DOwjoIcZbBlcAP44oi2GExrqBGkbCiMingdtV9c4a590DvEZVDzZGMqPTMQvAMJLn/wK91U4QkdnAJ23zNxqJWQCGYRgdilkAhmEYHYopAMMwjA7FFIBhGEaHYgrAMAyjQzEFYBiG0aH8fxsaw/MCnQ8tAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
    " + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" } ], "source": [ - "swiftdiff" + "swiftdiff['vz'].sel(id=4).plot.line(x=\"time (d)\")" ] }, { "cell_type": "code", - "execution_count": 9, + "execution_count": 27, "metadata": {}, "outputs": [ { @@ -1745,35 +1061,35 @@ " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", - "
    <xarray.DataArray 'vx' (time (y): 5)>\n",
    -       "array([ 0.00000000e+00, -4.35675281e-11,  2.08934259e-12, -1.73709378e-10,\n",
    -       "       -2.60294819e-10])\n",
    +       "
    <xarray.DataArray 'vz' (time (d): 3651)>\n",
    +       "array([ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00, ...,\n",
    +       "       -2.48154184e-22, -2.60768688e-22, -2.74623963e-22])\n",
            "Coordinates:\n",
    -       "    id        int64 6000200\n",
    -       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05
    " + " id int64 4\n", + " * time (d) (time (d)) float64 0.0 1.0 2.0 ... 3.648e+03 3.649e+03 3.65e+03
    " ], "text/plain": [ - "\n", - "array([ 0.00000000e+00, -4.35675281e-11, 2.08934259e-12, -1.73709378e-10,\n", - " -2.60294819e-10])\n", + "\n", + "array([ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, ...,\n", + " -2.48154184e-22, -2.60768688e-22, -2.74623963e-22])\n", "Coordinates:\n", - " id int64 6000200\n", - " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05" + " id int64 4\n", + " * time (d) (time (d)) float64 0.0 1.0 2.0 ... 3.648e+03 3.649e+03 3.65e+03" ] }, - "execution_count": 9, + "execution_count": 27, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "swiftdiff['vx'].sel(id=6000200)" + "swiftdiff['vz'].sel(id=4)" ] }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 12, "metadata": {}, "outputs": [ { @@ -2130,30 +1446,27 @@ " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", - "
    <xarray.DataArray 'px' (time (y): 5)>\n",
    -       "array([ 0.00000000e+00,  7.07309283e-09, -5.01222758e-09,  1.41219044e-08,\n",
    -       "        4.51840986e-08])\n",
    +       "
    <xarray.DataArray 'pz' (time (d): 3651)>\n",
    +       "array([0., 0., 0., ..., 0., 0., 0.])\n",
            "Coordinates:\n",
    -       "    id        int64 6000200\n",
    -       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05
    " + " id int64 6000204\n", + " * time (d) (time (d)) float64 0.0 1.0 2.0 ... 3.648e+03 3.649e+03 3.65e+03
    " ], "text/plain": [ - "\n", - "array([ 0.00000000e+00, 7.07309283e-09, -5.01222758e-09, 1.41219044e-08,\n", - " 4.51840986e-08])\n", + "\n", + "array([0., 0., 0., ..., 0., 0., 0.])\n", "Coordinates:\n", - " id int64 6000200\n", - " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05" + " id int64 6000204\n", + " * time (d) (time (d)) float64 0.0 1.0 2.0 ... 3.648e+03 3.649e+03 3.65e+03" ] }, - "execution_count": 10, + "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "swiftdiff['px'].sel(id=6000200)" + "swiftdiff['pz'].sel(id=6000204)" ] }, {