{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'vel_true')" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgcAAAGzCAYAAAC7ErTFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABHgElEQVR4nO3de3RU5b3/8c+ea8IlQSwmUAPFesErpaAYsb9apaVUWd6OtdYeqUfbVU+0CnXZ0qpY2xJti6I9CGoVejkYL5WqPVWqtODxCF5iOUtrRal4iEqgegoJgSQzs/fvD2vqHDLPd5KdZGbI++Watcw8s2c/s2cmPHn2/jxfLwiCQAAAAH8XKXQHAABAcWFwAAAAsjA4AAAAWRgcAACALAwOAABAFgYHAAAgC4MDAACQhcEBAADIwuAAAABkYXAA7EPWrFkjz/O0Zs2aQncFQAljcAAAALIwOAAAAFkYHAAAgCwMDoACeuCBB+R5ntauXbtX2+233y7P8/TSSy9Jkl555RX90z/9k0aOHKmysjJNmTJFDz/88EB3GcAgwOAAKKBTTz1Vw4YN03333bdX27333qsjjzxSRx11lP70pz/p+OOP15///Gd961vf0sKFCzV06FCdccYZWrlyZQF6DmBfxuAAKKDy8nLNmjVLDzzwgDKZTNf9zc3NWrt2rc4991xJ0uWXX66xY8fqhRde0FVXXaW6ujqtWbNGtbW1+uY3v1mo7gPYRzE4AArs3HPP1fbt27Pihw888IB839e5556r//3f/9Xvf/97ff7zn1dra6veeecdvfPOO3r33Xc1Y8YMvfbaa3rrrbcK9wIA7HMYHAAF9tnPflaVlZW69957u+6799579bGPfUyHHnqoNm3apCAIdM0112jUqFFZt/nz50uStm/fXqjuA9gHxQrdAWCwSyaTXdcO3Hbbbdq2bZv+67/+SwsWLJAk+b4vSbryyis1Y8aMbp/j4IMPHrD+Atj3MTgAisC5556rn/3sZ1q9erX+/Oc/KwiCrusNDjroIElSPB7X9OnTC9lNAIMEpxWAIjB9+nSNHDlS9957r+69914dd9xxGj9+vCTpgAMO0EknnaTbb79dW7du3Wvbv/71rwPdXQD7OGYOgCIQj8d11llnqaGhQW1tbfrxj3+c1b548WKdeOKJOvroo/WVr3xFBx10kLZt26Z169bpzTff1H//938XqOcA9kXMHABF4txzz9WuXbskSZ///Oez2o444gg9//zzOvXUU7V8+XLV1dVp6dKlikQiuvbaawvRXQD7MC8IgqDQnQAAAMWDmQMAAJCFwQEAAMjC4AAAAGRhcAAAALIwOAAAAFkYHAAAgCxFtwiS7/t6++23NXz4cHmeV+juAACKWBAEam1t1ZgxYxSJ9N/fu+3t7ers7Az9PIlEQmVlZX3Qo/5VdIODt99+WzU1NYXuBgCghDQ1NenAAw/sl+dub2/X+HHD1Lw9E/q5qqurtXnz5qIfIBTd4GD48OGSpAl3f13RIckC9wbI5nmsGQb0VBD03yxwZneHXvmXW7v+7egPnZ2dat6e0ebGcaoY3vvZiZZWX+Mn/486OzsZHPTU+6cSokOSDA5QdBgcAD3Xn4OD9w3EaeiK4ZFQg4PeeOutt/TNb35Tjz76qHbv3q2DDz5Yy5Yt05QpU7p9/IMPPqglS5Zow4YN6ujo0JFHHqnrrrsuZ7n3XLggEQCAPGQCP/StJ/72t79p2rRpisfjevTRR/Xyyy9r4cKF2m+//XJu8+STT+rTn/60fvvb36qxsVGf+tSnNGvWLP3xj3/s0b6LbuYAAIBi5CuQr97PHvZ02xtvvFE1NTVatmxZ133vl3LPZdGiRVk/L1iwQA899JAeeeQRTZo0Ke99M3MAAEAe/D74T5JaWlqybh0dHd3u7+GHH9aUKVN0zjnn6IADDtCkSZN055139qzPvq/W1laNHDmyR9vtkzMHfsjzW4U+qzwQ5+d6q9Dn3Pv7yESM1xf2vSn08XMp5r4NBsX8vQ8r7GdrXzs2/zeRN3/+fF133XV7Pe7111/XkiVLNHfuXH3729/Wc889p69//etKJBKaPXt2Xvv68Y9/rF27du1VBt6yTw4OAADoa5kgUCbo/UDn/W2bmppUUVHRdX8y2f3F977va8qUKVqwYIEkadKkSXrppZe0dOnSvAYHK1as0He/+1099NBDOuCAA3rUV04rAACQh/evOQhzk6SKioqsW67BwejRo3XEEUdk3Xf44Ydry5YtZl8bGhp08cUX67777tP06dN7/FoZHAAAUISmTZumjRs3Zt336quvaty4cc7t7rnnHl144YW65557dOqpp/Zq3wwOAADIg69AmRC3nqYV5syZo/Xr12vBggXatGmTVqxYoTvuuEN1dXVdj5k3b54uuOCCrp9XrFihCy64QAsXLtTUqVPV3Nys5uZm7dy5s0f7ZnAAAEAe+uq0Qr6OPfZYrVy5Uvfcc4+OOuoofe9739OiRYt0/vnndz1m69atWacZ7rjjDqXTadXV1Wn06NFdt8svv7xH++aCRAAAitRpp52m0047LWf78uXLs35es2ZNn+yXwQEAAHnoq7RCKSjJwUHGyLxamVi73b1/33efjQn7/rv619+fLWt58rB55bDPb7eHe37r1Vlpa2udBEuY4xvp5/fO3n/p/OIrBHv9lf5dY8PF78ffWfsS/++3MNuXCq45AAAAWUpy5gAAgIH2fuogzPalgsEBAAB5yATv3cJsXyoYHAAAkAeuOQAAAIMWMwcAAOTBl6dMiNqwfr/Xle07RTs48ANPXo54jBUl9H33G2C2Z6znN6KMxtxRYDx/4Opf2HkpM4sXLipob2+0G3NZ9vbuAxSNurePGNtHjNdniYaI+0WNfUetvlvHrsc96tnz7+tRRyuq2J+l5K3ntqKGuX7X5ru9jPc2Y/zOLRV+EC72GTYyOpA4rQAAALIU7cwBAADFJBPytEKYbQcagwMAAPIwmAYHnFYAAABZmDkAACAPfuCFurA07EWpA4nBAQAAeeC0AgAAGLSKdubAvc6Be/SVMdYRyKSj7n2njHUIOo0xVdrIFBvPH007GsPmha2Sxe5Do8BYJ8Aabvrm9ka7sb0XM/LWMWsdBKM95PaWRCyT+7mtNRyM9zZmbW+0xzx3u7UGRamvgxB2HQNrrQBrgZy0Y30V67ld20p2qfKMsb3VbpZKt9ZhcGzf36XIPyijiDIh/qbO/e0uPkU7OAAAoJgEIa85MBeTKiIMDgAAyAPXHAAAgEGLmQMAAPKQCSLKBCGuOSjuy2qyMDgAACAPvjz5ISbcffPSz+LBaQUAAJClaGcOgsDLeWWnFWVMp4yoYoe73Wt3t0d3G1HEdmezoh3u/kccUUYjTWZnkoyay74VZTSGk4HxifKNqGEQD7e9nwjb7j7AmYQ7jBSLu9s9a/8h/rCwoopJZ0ZWikfdfU+4Pph57N+KQtpRx7D1ysPxjQ9/2mo34n5p48vnen7ruVPGc1tRRIt5BX8JXaXvMpguSCzawQEAAMUk/DUHnFYAAAAlipkDAADy8N4FiSEKL3FaAQCAfYsfcvlk0goAAKBkMXMAAEAeBtMFiUU7OPB9ycsRWfSNqotBxjivY1RF9DqNqGHK/fTRdvf2MTPqmLvNSxtVCUOmvcyoYtSIQlpRxrixvRVlTBgVOZNGVNBqNz47vnF83WG/cJUL40bFR+u5raqLVlSxzIhCWlFJ6/nNqpJeYWvaWVHGlNHeaXw5rLhhp6O9M+N+bismmvKM+LcRRcxYVReNCHWprBzoKzJoFkEq2sEBAADFJBN4yoRYsyHMtgONaw4AAEAWZg4AAMhDJmRaIcNpBQAA9i1+EDGvPXFvXzqDA04rAACALMwcAACQB04rAACALL7CJQ4KW1e0Z4p2cOD7Xu51DozyooGxDoJnZNkjKXe7VXLZtU6BJEX3uNtje3KPLqOdRk7eCNobcWSzsmoQNdYJiBl5ZnMdA3d7JmmU6y439m+8txkjcG3lsY1q4kpHjLUIHGsZWGV1c5U4z5dVcrk86l7gw27vdLbHjXUM4hF3e6Sf/ypLBe61AKx2q6zyHuPD3+FYy6DDODbtxjoI7cYX0+p71PjsmSWdS6jmwGBRtIMDAACKSfhFkErnMj8GBwAA5CH88smlMzgonZ4CAIABwcwBAAB58OXJD3F9RJhtBxqDAwAA8jCYTiswOAAAIA/h1zlgcBBakInkLM1sldVVOly7Z8QBI+5ElhlldEUVJSm+O3d7rN0dN4tYUceQtVHNqKMVZUy4vxxmVLHMvf+IUW47bUU9rSyi8eW2ooqZiBGHc8QJ0zF3XM2Km1liRr1vqyTz0Jj7gz886q5VnjRqoZcZX8xI2HrlBmvZ3PbA/eu0w6hHXu67X/+eSO7trRhkxLPajZLOGffn1orZWiWdrXLjYWO66LmiHRwAAFBM/MDLY80G9/algsEBAAB58EOeViildQ5Kp6cAAGBAMHMAAEAewpdsLp2/xxkcAACQh4w8ZUKsVRBm24FWOsMYAAAGmbfeektf+tKXtP/++6u8vFxHH320nn/++ZyP37p1q774xS/q0EMPVSQS0RVXXNGr/TI4AAAgD++fVghz64m//e1vmjZtmuLxuB599FG9/PLLWrhwofbbb7+c23R0dGjUqFG6+uqrNXHixF6/1lCnFW644QbNmzdPl19+uRYtWiRJam9v1ze+8Q01NDSoo6NDM2bM0G233aaqqqoePXcmE8lZejlIGwfYWMfAysJHjXYjjm2WVY51WOsc5M5rx3a7s+7Rdne7lzLaM0ZW3HMfmyDmfm98oz0zxP2RTA9xbx9JGeW6zXUMjNcXMdqjxuszPrqZWO7j3xl3Z82TcePYGzEqK+ueMOqBDzEWABlmrHMwPOJuLzO+eHFjHYSoUdLZmvK11zlwr2PQbqxzsNtPutsjudcqiBnf27DlrHdH3eskdPruz6YXcg2OYpFRuFMD7t++e7vxxhtVU1OjZcuWdd03fvx45zYf+chHdMstt0iS7r777p52sUuv37HnnntOt99+u4455pis++fMmaNHHnlE999/v9auXau3335bZ511Vq87CADAvqSlpSXr1tHR/QJiDz/8sKZMmaJzzjlHBxxwgCZNmqQ777xzQPrYq8HBrl27dP755+vOO+/Mmt7YuXOn7rrrLt100006+eSTNXnyZC1btkxPP/201q9f32edBgBgoPXVaYWamhpVVlZ23err67vd3+uvv64lS5bokEMO0apVq3TJJZfo61//un72s5/1+2vt1WmFuro6nXrqqZo+fbq+//3vd93f2NioVCql6dOnd903YcIEjR07VuvWrdPxxx+/13N1dHRkjZpaWlp60yUAAPpVXxVeampqUkVFRdf9yWT3p5R839eUKVO0YMECSdKkSZP00ksvaenSpZo9e3av+5GPHg8OGhoa9MILL+i5557bq625uVmJREIjRozIur+qqkrNzc3dPl99fb2++93v9rQbAAAMqCBkyebg79tWVFRkDQ5yGT16tI444ois+w4//HD96le/6nUf8tWjIVBTU5Muv/xy/fu//7vKyowKOHmaN2+edu7c2XVramrqk+cFAKCUTZs2TRs3bsy679VXX9W4ceP6fd89mjlobGzU9u3b9fGPf7zrvkwmoyeffFL/9m//plWrVqmzs1M7duzImj3Ytm2bqquru33OZDKZc0oFAIBi0VenFfI1Z84cnXDCCVqwYIE+//nP69lnn9Udd9yhO+64o+sx8+bN01tvvaWf//znXfdt2LBB0nvXB/71r3/Vhg0blEgk9pqFcOnR4OCUU07Riy++mHXfhRdeqAkTJuib3/ymampqFI/HtXr1ap199tmSpI0bN2rLli2qra3tya4U+J6CXKWZU1ZU0Yi7hS3JbLTb27tjRdGO3LGk6B4jrtXm3rnX7o6DeSnj4ARGaVUjyhdNuONc0XajtGynu11G2dzAKJkcuJvlG98YP27E4cz23B1IGVHGlFEOO2zZ27jnDmINMWqVj4judrZbUcahEffzh40yWqwImx1ldH92W42SzfFMec62sFFFa6p8V8r9B1xHxP3FsGKypWKgqzIee+yxWrlypebNm6frr79e48eP16JFi3T++ed3PWbr1q3asmVL1naTJk3q+v/GxkatWLFC48aN0xtvvJH3vns0OBg+fLiOOuqorPuGDh2q/fffv+v+iy66SHPnztXIkSNVUVGhyy67TLW1td1ejAgAAHI77bTTdNppp+VsX758+V73BcYfcfno89oKN998syKRiM4+++ysRZAAAChlmZAlm8NsO9BCDw7WrFmT9XNZWZkWL16sxYsXh31qAACKxkCfViik0hnGAACAAUHJZgAA8uArIj/E39Rhth1oDA4AAMhDJvCUCXFqIMy2A61oBwdBOpKz+qJnVGX0zKij0W6k+YzElKIp95WiEau9M3eUMWJEGb3d7riXt9sdFwtyFADpYlVtNKKMnrGmRSTlXlzL+sAG0XBRwowRB4wmjedPhGt37T+TdkcZ0xl3uxVXi3ju99aMMhoZ3iFGFHFEtM3ZPtRzP7/VPyvKaMXtrPPF7UaMts0LV1UyIuO752C99x1GRjcRdR9b63ea9U+iUexUmX0jCVlSinZwAABAMRlMFyQyOAAAIA/BByor9nb7UsHgAACAPGTkmStlWtuXitIZxgAAgAHBzAEAAHnwg3DXDfgldGElgwMAAPLgh7zmIMy2A610egoAAAZE8c4cZLz3bt1JG3ltdyQ39DoGkbSR6TX3b2yfzp1n9jLuJ/c63VnqoN1Y56DNXVY3sNY5MALLEaN/Zh7aKPkcKXN/pKMd7rUAou7uKRK63VgHwfHZzlnC/O8yvnusHzZGFTe+WGVGjt8quWytYzDEOLhlxhcvHrKssfXhTAbG/o12ax0G1+p6KWONhd1GuWhrjYiYsQaGF7Ikc6lMt/vyzDUjrO1LRfEODgAAKCKDaYVETisAAIAszBwAAJCHwXRBIoMDAADy4Cvk8skldM1B6QxjAADAgGDmAACAPAQh0wpBCc0cFO/gwPfeu3XDSNWYUUKz3Xr+QsZurKqtVtQw5c5pBulw7fLcH/7AaJf1/GkjymnUdrXarZipGYM14obmZ88RZfQz7om+wPhcWtOhUeODHTGidlbJYSvqWGZub0Up3f0rMz571jSq9dWLBsYjfPfrS5nHJ3e7We7Z+KVmlasOq0SSiiaqMgIAgCyD6YLE0ukpAAAYEMwcAACQB04rAACALINp+WROKwAAgCzMHAAAkAdOKwAAgCwMDopdyNCs9f5YaZPw7UYHoo52V5skRY2dW/vubxF3/zxrnYS+7EsvdmCugRGy3RmmN7YN+vkXT9TsfDhW1t4quewu5i3FjfO9EWsNDpP7+GSM47fbWMchYq600HvWP1pWe39/9jDwSnNwAADAAGPmAAAAZBlMgwPSCgAAIAszBwAA5CFQuLUKSqnGBIMDAADyMJhOKzA4AAAgDwwOioEX5K6NbEUFo0a78aqtdj9mxH6MTFUm4d4+E8/9AqNx94sL4kbnk0lns2eUdFbUvX8riugljIMTM/ofM95cK6ppfTfDfneteUPzl0PvJx6tOJnVnjHb+/cSJfMXp5kADhdVjMv4bFmM/rUbUcV4P0ZFrWqA1rFPW9sbL76U/lHEe4p3cAAAQBFh5gAAAGQZTIMDoowAACALMwcAAOQhCLxQS0WX0jLTDA4AAMiDLy/UOgdhth1onFYAAABZinfmIBq8d+tGEHfHvXx3YkhGokiZjBHHM7aPGGnAdMr9/NHO3JGqTIc7CuilypztkUy4uJSXNl6cFWU0opRBubvdL3N/ZP2EEbkyYqhmRc3+jkL2I7vyXri4WsrIAFtRyIwVh3O22qLG85tRSONvKT9w99CqGhmGb/TNOvYpI/9N1cb3DKYLEot3cAAAQBEZTNcccFoBAABkYeYAAIA8cFoBAABkGUynFRgcAACQhyDkzEEpDQ645gAAAGRh5gAAgDwEkoLeF04NUXN14BXt4MCL+/Li3eeGAz9cXjv0O2RMDXlG/6zKrJFM7v57GavetHudgFjYksoZY5EHg1VS2i9LuHdf7u5fusyd17bKZZvrIETDrZOQswx5H7B+aYXNqqd897G1svLtgfu9s7bvNA5up7HOQNI49v25DkF/C7+OQbiSzv19oZ1rOn4gp+p9efJYIREAABTSW2+9pS996Uvaf//9VV5erqOPPlrPP/+8c5s1a9bo4x//uJLJpA4++GAtX768x/tlcAAAQB7eTyuEufXE3/72N02bNk3xeFyPPvqoXn75ZS1cuFD77bdfzm02b96sU089VZ/61Ke0YcMGXXHFFbr44ou1atWqHu27aE8rAABQTPzAkzeA6xzceOONqqmp0bJly7ruGz9+vHObpUuXavz48Vq4cKEk6fDDD9dTTz2lm2++WTNmzMh738wcAAAwgFpaWrJuHR0d3T7u4Ycf1pQpU3TOOefogAMO0KRJk3TnnXc6n3vdunWaPn161n0zZszQunXretRHBgcAAOQhCMLfJKmmpkaVlZVdt/r6+m739/rrr2vJkiU65JBDtGrVKl1yySX6+te/rp/97Gc5+9jc3Kyqqqqs+6qqqtTS0qI9e/bk/Vo5rQAAQB76aoXEpqYmVVRUdN2fzFGt1vd9TZkyRQsWLJAkTZo0SS+99JKWLl2q2bNn97of+SjawUEk5iuSI8polW613jvfMyZMQqdNjCewOuhMXLkjSUHEiBwl3NtH292BrkgqXOFcP25Epoz+pcuN9iFG6dqk+/hk3ElK+UbezUzRWnN1/TiXZ4UoO42o4m7ffXB2++4YrdXeFul+avV98cAdo40GKXe7VcvdOPZWyeeUEaV09y6PKKnju5+xykmbfbeijO72/o4Teo4YqqutWFVUVGQNDnIZPXq0jjjiiKz7Dj/8cP3qV7/KuU11dbW2bduWdd+2bdtUUVGh8vLyvPtYtIMDAACKyUDXVpg2bZo2btyYdd+rr76qcePG5dymtrZWv/3tb7Pue/zxx1VbW9ujfXPNAQAAeXi/KmOYW0/MmTNH69ev14IFC7Rp0yatWLFCd9xxh+rq6roeM2/ePF1wwQVdP3/ta1/T66+/rquuukqvvPKKbrvtNt13332aM2dOj/bN4AAAgDz01QWJ+Tr22GO1cuVK3XPPPTrqqKP0ve99T4sWLdL555/f9ZitW7dqy5YtXT+PHz9e//Ef/6HHH39cEydO1MKFC/XTn/60RzFGidMKAAAUrdNOO02nnXZazvbuVj886aST9Mc//jHUfhkcAACQh/f++g9zzUEfdqafMTgAACAPA31BYiFxzQEAAMjSo5mDJUuWaMmSJXrjjTckSUceeaSuvfZazZw5U5LU3t6ub3zjG2poaFBHR4dmzJih2267ba/VmvIRjfqKRLvPDVtjr4y1jEDEPbeTiRiZXqvd6qBRNtm1VoFReVWZuPsB8TL3vqMd7u29jDEvZjUbJY/9uLUOgdFurWPgjtqb21vrHJjrIBjvn/OzY32ujQ9exjdKHvvuXwd7Mu4XtzPjzlAnI+6kf8SoZW6VJc5ot7NdkU739sY6CAnje9tpzBm3G+9Pe+A+/ilHu1WS2Sq3ba1jkDYX6HDr75LOAyWQvV6ItX2p6NE7fuCBB+qGG25QY2Ojnn/+eZ188sk6/fTT9ac//UnSe7GLRx55RPfff7/Wrl2rt99+W2eddVa/dBwAgIE00FUZC6lHMwezZs3K+vkHP/iBlixZovXr1+vAAw/UXXfdpRUrVujkk0+WJC1btkyHH3641q9fr+OPP77veg0AAPpNr+eKMpmMGhoa1NbWptraWjU2NiqVSmVVg5owYYLGjh3rrAbV0dGxV4UqAACKTtAHtxLR48HBiy++qGHDhimZTOprX/uaVq5cqSOOOELNzc1KJBIaMWJE1uOrqqrU3Nyc8/nq6+uzqlPV1NT0+EUAANDvwp5SKKHTCj0eHBx22GHasGGDnnnmGV1yySWaPXu2Xn755V53YN68edq5c2fXrampqdfPBQBAfxnoFRILqcfrHCQSCR188MGSpMmTJ+u5557TLbfconPPPVednZ3asWNH1uzBtm3bVF1dnfP5kslkznKVAABg4IVeBMn3fXV0dGjy5MmKx+NavXq1zj77bEnSxo0btWXLlh5Xg5KkaMxXNNZ9tMg3oohexIhERYxYj1ny2d1uRh2tsr6OuJ8Z9TPGWely986NtJkiRtVbM8povnZ3u5G2M4+PGUU0SjbbJZ3dB8B6fYo6tjejjO72VMZ98Pekjahiyooqpt0dMFhxvPaou3/tgdEedUcdKzyjZLQRtbS0GVHFFr/Mvb2j5LVVTts6th1m1NH4XpXQdHkYg2kRpB4NDubNm6eZM2dq7Nixam1t1YoVK7RmzRqtWrVKlZWVuuiiizR37lyNHDlSFRUVuuyyy1RbW0tSAQBQ+sJeN7CvDg62b9+uCy64QFu3blVlZaWOOeYYrVq1Sp/+9KclSTfffLMikYjOPvvsrEWQAABA6ejR4OCuu+5ytpeVlWnx4sVavHhxqE4BAFBswl5UuE9fkAgAwKA0iNZPpvASAADIwswBAAB5IK1QBKJRX9EcVRkjRpQx4qhqKJlFEZUxnt+PWlFFdywoY23viOPZVQfd7ZGU1e5slpXm8qwooxkDdbdbc11GIsuOUppRSeOzEXL7wBllDDcnmc64D86elDsK2BJxR+0ixpxph3Fwdhs50d0xd063NeaOWrb67vbhkT3O9jLP/eWIGl8OK2rZalS13JEZkrNtV8b93uwxjm06ZJRxUCmhUwNhcFoBAABkKdqZAwAAigmnFQAAQLZBlFZgcAAAQF48meuYm9uXBq45AAAAWZg5AAAgH5xWAAAAWRgcFF7UCxTNkeu2Cqd6rqy47HUOcq2v8L60UZI5Y+zfjxvbO9r9pLFtWbh1DjyjJHPIqrW2kKfkzHUSrHUWQq7DEMSMdQyMdQ7kWGPD+twGvrtzaatkc8r6deDO0qeMrPyutHudgp3GOgU7Yrlz/pJUEXOvUzAs5i7JXBl1bz8k4t4+YqxDYZVNtsouu9Yy+FvKfWxajGPfnnG/9xnjgx92HYRSuop/sCjawQEAAEWFks0AAOCDBlNVRtIKAAAgCzMHAADkgwsSAQBAlkF0zQGnFQAAQJainTnwvEBejmhQrohjvgJje98o+WxFHTMxo92KMjriin7aiBSljPaMMXL1rbycu9lq98KOnMNOy4Xd3jo8RozVFVWUJIUo2Wy9tHTaKNlsvDcpo+Tz7k53SeKY8b1JRN052mQs7WwfGu90tg+JuduHG1HH8qh7+7iRA84Yf4tZZZP3ZHIfX1ebJO1Ou2OSe9Lu7a33PmN8dsJGFV3bD2QM0gvCVU4P+U/XgCrawQEAAEWFaw4AAEAWrjkAAACDFTMHAADkg9MKAAAgyyAaHHBaAQAAZGHmAACAfAyimYOiHRy41jkwliHI47nD1R22crW+cVSt7TOO0ruuNknyjXUKMkbZXrMwSMg8s/n0Rv+tdmudBnP7kOs4mMyS0r3fQWCsYZExSgZnjHUQOjvdH2zra5nr+9zVbqwBEYu51xFIGO3JuHudhDJjHYVk1N0ejbh/r1glna2yx652q1y29XsjbbRb6xiYv5dK6B9FJ9IKAABgsCramQMAAIoJKyQCAIBsg+iaA04rAACALAwOAABAFk4rAACQB08hrznos570v5IcHFiRKCsyZG1vlYS2nt9qt7giS1bcydqztb3FjHGakSdj+5CRK6ssccaI+1n7DxuFNMvLhvjoBEZM1Xcn8exy3Wa78fwW4+k7jajjnpjxvTdKqVul2KMht48YUcdi/ofDimkOZNnkgiLKCAAABquSnDkAAGDADaK0AoMDAADyMYgGB5xWAACgCF133XXyPC/rNmHChJyPT6VSuv766/XRj35UZWVlmjhxoh577LFe7ZuZAwAA8lCIFRKPPPJIPfHEE10/x2K5/9m++uqr9ctf/lJ33nmnJkyYoFWrVunMM8/U008/rUmTJvVovwwOAADIRwFOK8RiMVVXV+f12F/84hf6zne+o8997nOSpEsuuURPPPGEFi5cqF/+8pc922+PezpAIl7u6otWVNCK3VhRxXjUXd0tZlR1tPYfs6q3hfj0RUJWnLT4gRGXMwJZVhTRqi7XkXZ/ZNuN9pQRdUxl3O3ptPH6rTihFUV1RC2tqKJVlVFpo93Y3ku59x8JGZW0/qoyU2BG1NGqlpo2opCpqNFBY/8Kub3n2N6qaBmxntv6nRU2xtnP8e9S09LSkvVzMplUMpns9rGvvfaaxowZo7KyMtXW1qq+vl5jx47t9rEdHR0qKyvLuq+8vFxPPfVUj/vINQcAAOQj6IObpJqaGlVWVnbd6uvru93d1KlTtXz5cj322GNasmSJNm/erE984hNqbW3t9vEzZszQTTfdpNdee02+7+vxxx/Xgw8+qK1bt/b4pRbtzAEAAMWkr645aGpqUkVFRdf9uWYNZs6c2fX/xxxzjKZOnapx48bpvvvu00UXXbTX42+55RZ95Stf0YQJE+R5nj760Y/qwgsv1N13393jvjJzAADAAKqoqMi65Roc/F8jRozQoYceqk2bNnXbPmrUKP36179WW1ub/ud//kevvPKKhg0bpoMOOqjHfWRwAABAPt5fPjnMLYRdu3bpL3/5i0aPHu18XFlZmT784Q8rnU7rV7/6lU4//fQe74vBAQAA+eijaw7ydeWVV2rt2rV644039PTTT+vMM89UNBrVeeedJ0m64IILNG/evK7HP/PMM3rwwQf1+uuv6z//8z/12c9+Vr7v66qrrurxS+WaAwAA8jDQ6xy8+eabOu+88/Tuu+9q1KhROvHEE7V+/XqNGjVKkrRlyxZFIv/4G7+9vV1XX321Xn/9dQ0bNkyf+9zn9Itf/EIjRozocV8ZHAAAUIQaGhqc7WvWrMn6+ZOf/KRefvnlPtl30Q4OPC/IWVo5bMllax2DeMTdXhZ1B7oTVrvx/K79W2ssWOscmOWoQy7+ba1z0JFxf+T2ZOLO9nZj+7aU+8KetlTC2W6to9ARsdZBCLdOQsaxzkNgvTVh1ynocG8f6QzZnnI2K+wSHYFnrCHhfmvsdmMdBHN74ySubz2/o93a1mpX3CjJnHT/zlLCWOQi5DoIRWMQ1VYo2sEBAABFJeRphVIaHHBBIgAAyMLMAQAA+eC0AgAAyDKIBgecVgAAAFmYOQAAIA8Dvc5BIRXt4CCiIGe8xYrjWSWTw0YVy2LuTFZ5NFx7wlH7NmnUxS0z8mJxz/3a+zvK2O67o4pWlLElXe5st8phW5+NXV5+a5z3lm+ULXY2W0uvWlFGI2oYbTfa9xjtHc5mRTqNdqvks/XRtFamNdqNauEKouGiklbJ6CDmfv5MwhFldH9t5CeNqGPS2LdZ8tmIyVrltI3PdslEHfchnFYAAABZinbmAACAojKILkhkcAAAQB645gAAAOythP6BD4NrDgAAQBZmDgAAyAfXHAAAgA/imoMc6uvr9eCDD+qVV15ReXm5TjjhBN1444067LDDuh7T3t6ub3zjG2poaFBHR4dmzJih2267TVVVVT3qWCTi58ykW1l2qyRywijZbJVcttYpGBp1B7qHxtyB8CGOQPjwaLtz22TIdQ4SnhU2d8sYdWnbA3cge1emzNluvfdWyWorL23lra32jLGOQdrIe7t37m72UsY6BCHXOYjtce8/arTH2o1y4WZJZyNrbxxaa5mIIBJyHQOr3VqLIGF8thxrEWSM5TkyRsXljG+UEk+424O40R4YHUDR6dE1B2vXrlVdXZ3Wr1+vxx9/XKlUSp/5zGfU1tbW9Zg5c+bokUce0f3336+1a9fq7bff1llnndXnHQcAYEAFfXArET2aOXjssceyfl6+fLkOOOAANTY26v/9v/+nnTt36q677tKKFSt08sknS5KWLVumww8/XOvXr9fxxx/fdz0HAGAADabTCqHSCjt37pQkjRw5UpLU2NioVCql6dOndz1mwoQJGjt2rNatW9ftc3R0dKilpSXrBgAACqfXgwPf93XFFVdo2rRpOuqooyRJzc3NSiQSGjFiRNZjq6qq1Nzc3O3z1NfXq7KysutWU1PT2y4BANB/BtFphV4PDurq6vTSSy+poaEhVAfmzZunnTt3dt2amppCPR8AAP1iEA0OehVlvPTSS/Wb3/xGTz75pA488MCu+6urq9XZ2akdO3ZkzR5s27ZN1dXV3T5XMplUMtm/lfAAAED+ejQ4CIJAl112mVauXKk1a9Zo/PjxWe2TJ09WPB7X6tWrdfbZZ0uSNm7cqC1btqi2trZHHfOUu8KqVXbXao9ZcT4jCmlGGY2oohVHrHRkwoYbebEhEfe+yzyjXLRxbCwZoy5um+8eCFr9s6KKFt+IWqatSJexfWfGnWfzjCuSAlezVZI5bZS9NUomWyWXrahivM392qwoY7TTijpaUcZwf5YFXtiSzEYMNm7EXI2/kVzVyj3js2HlOM1q4Gn3594qRR5Y+3fvvmgMpgsSezQ4qKur04oVK/TQQw9p+PDhXdcRVFZWqry8XJWVlbrooos0d+5cjRw5UhUVFbrssstUW1tLUgEAUNpYIbF7S5YskSSddNJJWfcvW7ZMX/7ylyVJN998syKRiM4+++ysRZAAAChpDA66F+QxbVdWVqbFixdr8eLFve4UAAAoHGorAACQB645AAAA2QbRaYVQKyQCAIB9T9HOHES8IGcFPauyntVuVfZLRtyVCRNGe5lRXm6YEWV0xRVHRHe7t42482ZDPXeezaraaLGijNb+w0YpraiiVVUxbVWVzLhL68WNip8Rqyqjq9mIi0WMqowRoyqjkYI1o4jxPcb3brcRQW53H7toyoixZsJVbbRYVRt9ozJhJulujwxxt3vOmK0RkbUqTloxSzMq6W62oo6lgtMKAAAgG6cVAADAYMXMAQAA+RhEMwcMDgAAyINrWf98ty8VnFYAAABZmDkAACAfnFYAAAAfRJSxCDjXOTCGX2HXQbDarbUAkp57HYQhRu3c4ZHc6yCMiLjXORhhrHMwxFijIW4c26hx0syoqqsy49iFLcmcMc6UpYy6u3syCWd7W8z93u2KuOvuWiWbndsaWXFriQjjrTdLJkc7rHajVPoeYx2DPe4ORjrc7V7K/fyesQ6CJTA+/EHc/dmKlrnXyPB849exY/d+zFjfw/2xVsa9NItkfLaCsCWbrZrRxfKv6iCaOeCaAwAAkKVoZw4AACg6JfTXfxgMDgAAyMNguuaA0woAACALMwcAAORjEF2QyOAAAIA8DKbTCkU7OPADzyyv21t2FNIdyTKjjlZJZ6Ns8VBH7VxXzPG9dncmabhRujVunGmKeO7tU4FRDjtwH5uo3FFMqyRzZ9QdJ+vw3XGy3TF35qsl7Y4qWiWbo2F+OxibWinQsFFHs92IQkY6jGOz2/3Z9TrCtStjHICM+wB6UfdnL0gYUUUjaimVuZ8/lvu7l06GK7lsfnaMqOK+UpIZ/1C0gwMAAIoKpxUAAMAHDabTCqQVAABAFmYOAADIB6cVAABAFgYHAADgg7jmAAAAFNR1110nz/OybhMmTHBus2jRIh122GEqLy9XTU2N5syZo/Z2dwS+OyU5c+C7apdK5voI9vbhxkxRs+yxVRI6d6A8bgSSjbizkp77tZV57o9ExBhP+kb/4sY6B75RGzYVuNdBaA+MdQyi7u1bo+6s+TCjZHMyapQdjliBckeb9VdHyHUQZLRHjJLHZnvaWEfAKrnc7j72as+9PogkBZ3G9tY6CMZ3x0u618jwjHUUrJLPkY7c7ZG0e1vPtz48/btOQVBCfzE7FeC0wpFHHqknnnii6+dYLPfv6BUrVuhb3/qW7r77bp1wwgl69dVX9eUvf1me5+mmm27q0X5LcnAAAMBA84JAXoiRTm+2jcViqq6uzuuxTz/9tKZNm6YvfvGLkqSPfOQjOu+88/TMM8/0eL+cVgAAYAC1tLRk3To6cs96vfbaaxozZowOOuggnX/++dqyZUvOx55wwglqbGzUs88+K0l6/fXX9dvf/laf+9znetxHBgcAAOQj6IObpJqaGlVWVnbd6uvru93d1KlTtXz5cj322GNasmSJNm/erE984hNqbW3t9vFf/OIXdf311+vEE09UPB7XRz/6UZ100kn69re/3eOXymkFAADy0FdphaamJlVUVHTdn0x2X7Nl5syZXf9/zDHHaOrUqRo3bpzuu+8+XXTRRXs9fs2aNVqwYIFuu+02TZ06VZs2bdLll1+u733ve7rmmmt61FcGBwAADKCKioqswUG+RowYoUMPPVSbNm3qtv2aa67RP//zP+viiy+WJB199NFqa2vTV7/6VX3nO99RJJL/yQJOKwAAkI8+Oq3QW7t27dJf/vIXjR49utv23bt37zUAiP69Um3Qw4shi3bmIBN48nJEEtO+e0wTN8oSW9ubUUijPWPEgjLW9iGilO5AkxQ3HhEz2qNGnMvugVuHGdV0x82GOMpd59UedbcnjLrFVns04v6Ceo72fqpg3nesuJzZbmQp0+73Pki5SzYHZtTRKPls/F7x0kaM1Sh3Hil3lwOPpHMfPzOqWOAPj/HSS8ZAL4J05ZVXatasWRo3bpzefvttzZ8/X9FoVOedd54k6YILLtCHP/zhrmsWZs2apZtuukmTJk3qOq1wzTXXaNasWV2DhHwV7eAAAIDB7M0339R5552nd999V6NGjdKJJ56o9evXa9SoUZKkLVu2ZM0UXH311fI8T1dffbXeeustjRo1SrNmzdIPfvCDHu+bwQEAAPkY4EWQGhoanO1r1qzJ+jkWi2n+/PmaP39+Dzu2NwYHAADkYTDVVmBwAABAPgZRVUbSCgAAIAszBwAA5KmUTg2EUbSDgyDwckYGM0YU0WpP++5IR4fvPix2u7syYCpwb+87JnTsGKWbb5XeKzBrKithRB3LPHccLWz7kIi7sl8i6n4H4kZ7xPGbJxN1/1YKjKSSmZA12v2o+7MXxNxPEESNHVgLtBhRQqv0X2BEDYO0EWW0cuJGnDBIuH8vKOXun+eoehn6Hyzr0Brt+0pU0RQE4UpMllB5Sk4rAACALEU7cwAAQDEhrQAAALKRVgAAAIMVMwcAAOTB89+7hdm+VDA4AAAgH5xWAAAAg1XRzhwEgacgR7jWyvqnjUB32HZrnYN2Y52D9qD37e3GGgkdgTurnQrc81pxoySyHzKnmylwzjdhvD7r9ccj7narZHM84j7+kaij3cqiW+sUGN/2TNwoVR5zv3e+sX2QMNZBiBsLNcSNFxBzt3tGydrQn0zju2WWrLa+GyG+O4GxRoS1RoaMUuPWZfheyPZiQVoBAABkG0SLIDE4AAAgD4Np5oBrDgAAQBZmDgAAyMcgSiswOAAAIA+cVgAAAINW0c4cZIKIlKP0shV7sUo2pzLu3E5nxN2eNiJRYUs6t/nJnG3tESsG6S4p3BYYRZ19dxQy7oUbT1pRynbjat6MlecLKWosYRYx5gXjxvZRK8roiozF3NuaJZNjRlTRqCicSRoR4jJ3eyRlRA1TxrFLG8fOKHkso2RzxDeOr7G9Z0QprailYkbU0vH+WuW0rT8DzXLeg6Uks4W0AgAA+CBOKwAAgEGLmQMAAPJBWgEAAHwQpxUAAMCgxcwBAAD58AO7gJa1fYlgcAAAQD645iC3J598Uj/60Y/U2NiorVu3auXKlTrjjDO62oMg0Pz583XnnXdqx44dmjZtmpYsWaJDDjmkzzqdq5Tz+9LGOgfRiFWyOew6Bu723X7C2e5a56DFL3NuG/eMrLc6nK0ps6SxUZbWYG3dbry3u61y2UY57E7jve00SmKHFTFOOjrXOYi6t7W6nkkY+zbWKUin3e1exv29Mj86xntn/WL1jAy55xlhfavkc6d7DRBz+2FDnO2ZIe7fC5my3J9da40Kq1x36HUOBsk6CJ5CXnPQZz3pfz2+5qCtrU0TJ07U4sWLu23/4Q9/qFtvvVVLly7VM888o6FDh2rGjBlqb28P3VkAAND/evxn0syZMzVz5sxu24Ig0KJFi3T11Vfr9NNPlyT9/Oc/V1VVlX7961/rC1/4QrjeAgBQKINohcQ+TSts3rxZzc3Nmj59etd9lZWVmjp1qtatW9ftNh0dHWppacm6AQBQbN6PMoa5lYo+HRw0NzdLkqqqqrLur6qq6mr7v+rr61VZWdl1q6mp6csuAQCAHir4Ogfz5s3Tzp07u25NTU2F7hIAAHsL+uBWIvr00uzq6mpJ0rZt2zR69Oiu+7dt26aPfexj3W6TTCaVTOa+Oh8AgGLgBYGZirG2LxV9OjgYP368qqurtXr16q7BQEtLi5555hldcsklPXou33F+xjPibtbh90NGITv9cFHH3Rl3ZGlXJHdccYfnjkNZJYczxmRRm+eOa5UZUcew2o2ooRVVbDWinrsdMVHJLqdtvbcpMxPmFnGUdPasKGPc3e4n3J/7jDFGt6OIVlDLODZmXM743htli6MJ93sXKTcOQNr47BsllzNl7u99erh7/+khuY9f2iinbUUdg5jx2XJFbGVHdK0UKYpPjwcHu3bt0qZNm7p+3rx5szZs2KCRI0dq7NixuuKKK/T9739fhxxyiMaPH69rrrlGY8aMyVoLAQCAkuPLXqzF2r5E9Hhw8Pzzz+tTn/pU189z586VJM2ePVvLly/XVVddpba2Nn31q1/Vjh07dOKJJ+qxxx5TWZn7LzoAAIoZpxUcTjrpJAWOF+h5nq6//npdf/31oToGAAAKg9oKAADkg9oKAAAgyyBaIZHBAQAAeQi7ymEprZBYtIODIPByVl+0Bl9WVDFjRBWtdivquCfjzg3tjrojTa2Z3BdvxkNGCdsj7n2XRTqd7XG5929FKS0po7SgFWV0VbSU3Mc2n3Yrypg2Yq7WZ9PVGokZMdW4u9034m4Z8xeXkUczk4xWFNH9vfJjVhTT3R4b4n7voh3u74Z5gIy8XiZp/F4ZZlQMHZp7+0y5sW/jpVlVG2VEGT0rpVpK/ypCUhEPDgAAKCqcVgAAAB/k+XksBmZsXyoKXlsBAAAUF2YOAADIB6cVAABAlkG0zgGnFQAAQBZmDgAAyAO1FYqcbxzfXOsj/GN7o2SzUXa30wgFd/rutQDa0u4sfpi1DFJGyePdkQ5ne1nEXbI5YpQVi4bMM2eM9ybsOgi7jHUMdhl1i612q5y39dmLOko2u8o5S5KfMNr9cCWVrXfWKplsHBr58XAlpaPl7h1EO92vL5pyf7Y8c50Dd7O1ToNVdjk1NHd72l3JXX7SKOdtlPuWUdLZWsfAbHfvvXgMomsOOK0AAACyMDgAACAfgSQ/xK2HEwfXXXedPM/Luk2YMCHn40866aS9Hu95nk499dQevtASPa0AAMBAK8Q1B0ceeaSeeOKJrp9jsdz/bD/44IPq7PzHEvjvvvuuJk6cqHPOOafH+2VwAABAPgKFvOag55vEYjFVV1fn9diRI0dm/dzQ0KAhQ4b0anDAaQUAAAZQS0tL1q2jI/eF4q+99prGjBmjgw46SOeff762bNmS937uuusufeELX9DQoUN73EcGBwAA5OP9tEKYm6SamhpVVlZ23err67vd3dSpU7V8+XI99thjWrJkiTZv3qxPfOITam1tNbv67LPP6qWXXtLFF1/cq5e6T55WMEs2hyzpnMq4I1PtnjtOFzOqb0S83HE767XZUUZ37VYrRmm1R0JWFvGNGKn1+tp997G3Si5bUUWrHLdVstmK2boiXzGjZLNvlGzOGFFG650LQpZcNlKoZpQxUmbEQN3VxhVJGcc+7d7eszLURh7POn7GR9cZ5cyU935bSQqMKKMXtUo2uz89RjXr0uErXO7y74epqalJFRUVXXcnk92/QTNnzuz6/2OOOUZTp07VuHHjdN999+miiy5y7uquu+7S0UcfreOOO65XXd0nBwcAABSrioqKrMFBvkaMGKFDDz1UmzZtcj6ura1NDQ0Nuv7663vbRU4rAACQj/fTCmFuYezatUt/+ctfNHr0aOfj7r//fnV0dOhLX/pSr/fF4AAAgHz00TUH+bryyiu1du1avfHGG3r66ad15plnKhqN6rzzzpMkXXDBBZo3b95e2911110644wztP/++/f6pXJaAQCAIvTmm2/qvPPO07vvvqtRo0bpxBNP1Pr16zVq1ChJ0pYtWxSJZP+Nv3HjRj311FP63e9+F2rfDA4AAMjHANdWaGhocLavWbNmr/sOO+wwBX1Qw4HBAQAA+aDwEgAAGKyKduYgCDwzE557W/u5XdLGOgeeZ5SG9d2Z391pd6DZtZaBlaPf47vD2omIu91ex8B9cKNmWj6cjDGetdYxsMptW+sYtBvtVslmiysKH40a5bJjvS/1LUl+xP3eBhH3sc8YZX2DuPG9ci/BIeOjq7TRHskY6xxYh88seW0w10Ewdu9Yi8A3jp1Vstla5yBifPYixmdnn9FH6xyUgqIdHAAAUEwKUXipUBgcAACQD645AAAAgxUzBwAA5MMPJOO6K3P7EsHgAACAfHBaAQAADFYlOXNgRRF9a3RmRhWNWI+xfUcm3GF1lS1OG3knK8oXjxhRRRmv3To2IUs2h2VFPe2oo1ES2ogyWuW+rZLbzpLNRllcP2SU0ZouNaOOaeO1G1FHGVFD32i3trdKLpsfXatic8ioo5XcDhxlkwMrRpowOm+UAzdLNht9t36nlo6QMwfWh6iIlOTgAACAAcdpBQAAMFgxcwAAQD78QKFODZBWAABgHxP4793CbF8iOK0AAACyMHMAAEA+BtEFiUU7OPADT16ObI8VpzOjjsY5IyuO1t+cVRmNKGOnkceyooZ2VDFke8goj2+URLOiglbFzXTgjjJ2Ztzt1vtj9c91/CJGlDEWtWKSzmazsl7GjFK6X3tgRA2DjLG9FRU0ZmzNKq/WRzNsVcaw/y64dm9VRbSiiEaU0fpeW1FFMx5eKlFHrjkAAABZBtHMAdccAACALMwcAACQj0AhZw76rCf9jsEBAAD54LQCAAAYrJg5AAAgH74vMxZjbl8aGBwAAJCPQXRaYZ8cHJh5ZnN7d7u1DoKZ2TWy8q4sfMZYfjNlXPHS33njsOsghGWtI2C1W++ttY6BtX2YVx8xPtZRYx0Ci/nZMDrgGxlu31iHITD+qrLeO7PmscH63pvrLFhvbsj+hWKsg+CFbLfWyEDp2ScHBwAA9DlmDgAAQJZBtEIiaQUAAJCFmQMAAPIQBL6CEGWXw2w70BgcAACQjyAId2qAaw4AANjHBCGvOWBw0L/ClL2V7Khjxpj5iRpXaqSsqKIROcu4ooxhY5QGK85mKXSU0dLfUcewzx9Gf8dUPavsr1kS2t1ufS/7O6Jsb2/1z3qCcK+vP7855vfWijpaJZt73CMUWkkODgAAGHC+L3khrhvgmgMAAPYxg+i0AlFGAACQhZkDAADyEPi+ghCnFYgyAgCwr+G0AgAAGKyYOQAAIB9+IIWJBJfQzEG/DQ4WL16sH/3oR2pubtbEiRP1k5/8RMcdd1x/7S5L2HUQ7Od3t0eMVK9Z1tfRP7+f1yGwhM0rF3qdA4v12bF6b25fxOsgREOuD+IZ63eYywAUsqRxH+zfLPkcdp0Dx/OH7btdrtv93oZdH6VkBIGkMFHG0jlO/XJa4d5779XcuXM1f/58vfDCC5o4caJmzJih7du398fuAABAH+qXwcFNN92kr3zlK7rwwgt1xBFHaOnSpRoyZIjuvvvu/tgdAAD9LvCD0LdS0eeDg87OTjU2Nmr69On/2EkkounTp2vdunV7Pb6jo0MtLS1ZNwAAik7gh7+ViD4fHLzzzjvKZDKqqqrKur+qqkrNzc17Pb6+vl6VlZVdt5qamr7uEgAAoTFzMIDmzZunnTt3dt2ampoK3SUAAAa1Pk8rfOhDH1I0GtW2bduy7t+2bZuqq6v3enwymVQymez6Ofj71ZyZ3R193bV/7MOqIBay3bpw2NzecV132KuCw6YdSCu4hb0i3XccYTvp4Gzug6vxwyUtSCuE279z+7DHzvpeGhU5FTXSDMburd/JLu//WxEMQBIgHXSEOjWQVqoPe9O/+nxwkEgkNHnyZK1evVpnnHGGJMn3fa1evVqXXnqpuX1ra6sk6dWLbunrrgEA9lGtra2qrKzsl+dOJBKqrq7WU82/Df1c1dXVSiQSfdCr/tUv6xzMnTtXs2fP1pQpU3Tcccdp0aJFamtr04UXXmhuO2bMGDU1NWn48OHyPE8tLS2qqalRU1OTKioq+qO7+zSOX+9x7HqPYxcOxy9/QRCotbVVY8aM6bd9lJWVafPmzers7Az9XIlEQmVlZX3Qq/7VL4ODc889V3/961917bXXqrm5WR/72Mf02GOP7XWRYncikYgOPPDAve6vqKjgSxICx6/3OHa9x7ELh+OXn/6aMfigsrKykvhHva/02wqJl156aV6nEQAAQHEpeFoBAAAUl6IfHCSTSc2fPz8r0YD8cfx6j2PXexy7cDh+KDQvGIj8BwAAKBlFP3MAAAAGFoMDAACQhcEBAADIwuAAAABkYXAAAACyFP3gYPHixfrIRz6isrIyTZ06Vc8++2yhu1R0nnzySc2aNUtjxoyR53n69a9/ndUeBIGuvfZajR49WuXl5Zo+fbpee+21wnS2yNTX1+vYY4/V8OHDdcABB+iMM87Qxo0bsx7T3t6uuro67b///ho2bJjOPvvsvQqLDVZLlizRMccc07WSX21trR599NGudo5d/m644QZ5nqcrrrii6z6OHwqlqAcH9957r+bOnav58+frhRde0MSJEzVjxgxt37690F0rKm1tbZo4caIWL17cbfsPf/hD3XrrrVq6dKmeeeYZDR06VDNmzFB7e/sA97T4rF27VnV1dVq/fr0ef/xxpVIpfeYzn1FbW1vXY+bMmaNHHnlE999/v9auXau3335bZ511VgF7XTwOPPBA3XDDDWpsbNTzzz+vk08+Waeffrr+9Kc/SeLY5eu5557T7bffrmOOOSbrfo4fCiYoYscdd1xQV1fX9XMmkwnGjBkT1NfXF7BXxU1SsHLlyq6ffd8Pqqurgx/96Edd9+3YsSNIJpPBPffcU4AeFrft27cHkoK1a9cGQfDesYrH48H999/f9Zg///nPgaRg3bp1hepmUdtvv/2Cn/70pxy7PLW2tgaHHHJI8Pjjjwef/OQng8svvzwIAj57KKyinTno7OxUY2Ojpk+f3nVfJBLR9OnTtW7dugL2rLRs3rxZzc3NWcexsrJSU6dO5Th2Y+fOnZKkkSNHSpIaGxuVSqWyjt+ECRM0duxYjt//kclk1NDQoLa2NtXW1nLs8lRXV6dTTz016zhJfPZQWP1WeCmsd955R5lMZq9KjlVVVXrllVcK1KvS09zcLEndHsf32/Ae3/d1xRVXaNq0aTrqqKMkvXf8EomERowYkfVYjt8/vPjii6qtrVV7e7uGDRumlStX6ogjjtCGDRs4doaGhga98MILeu655/Zq47OHQirawQEw0Orq6vTSSy/pqaeeKnRXSsphhx2mDRs2aOfOnXrggQc0e/ZsrV27ttDdKnpNTU26/PLL9fjjjw+qUsAoDUV7WuFDH/qQotHoXlfmbtu2TdXV1QXqVel5/1hxHN0uvfRS/eY3v9Ef/vAHHXjggV33V1dXq7OzUzt27Mh6PMfvHxKJhA4++GBNnjxZ9fX1mjhxom655RaOnaGxsVHbt2/Xxz/+ccViMcViMa1du1a33nqrYrGYqqqqOH4omKIdHCQSCU2ePFmrV6/uus/3fa1evVq1tbUF7FlpGT9+vKqrq7OOY0tLi5555hmOo96LeV566aVauXKlfv/732v8+PFZ7ZMnT1Y8Hs86fhs3btSWLVs4fjn4vq+Ojg6OneGUU07Riy++qA0bNnTdpkyZovPPP7/r/zl+KJSiPq0wd+5czZ49W1OmTNFxxx2nRYsWqa2tTRdeeGGhu1ZUdu3apU2bNnX9vHnzZm3YsEEjR47U2LFjdcUVV+j73/++DjnkEI0fP17XXHONxowZozPOOKNwnS4SdXV1WrFihR566CENHz6861xuZWWlysvLVVlZqYsuukhz587VyJEjVVFRocsuu0y1tbU6/vjjC9z7wps3b55mzpypsWPHqrW1VStWrNCaNWu0atUqjp1h+PDhXde2vG/o0KHaf//9u+7n+KFgCh2XsPzkJz8Jxo4dGyQSieC4444L1q9fX+guFZ0//OEPgaS9brNnzw6C4L044zXXXBNUVVUFyWQyOOWUU4KNGzcWttNForvjJilYtmxZ12P27NkT/Ou//muw3377BUOGDAnOPPPMYOvWrYXrdBH5l3/5l2DcuHFBIpEIRo0aFZxyyinB7373u652jl3PfDDKGAQcPxSOFwRBUKBxCQAAKEJFe80BAAAoDAYHAAAgC4MDAACQhcEBAADIwuAAAABkYXAAAACyMDgAAABZGBwAAIAsDA4AAEAWBgcAACALgwMAAJDl/wPBdITsqw+j2QAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgcAAAGzCAYAAAC7ErTFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABG9klEQVR4nO3de3hU1b3/8c/MZGYSkIRLIQEJGG9cFJCCImLrjZZa8fHCUUuxUou2tXgB6tPT9CcFaQvaVkXPQRCLUNuD8UrVnipFKlCPIIrlPKIVL2CJyMUbBAJJJrP37w81x5FkfWeycxvyfj3PPDqz9tp7zZ49w8pa+7u+Id/3fQEAAHwq3NoNAAAAbQudAwAAkILOAQAASEHnAAAApKBzAAAAUtA5AAAAKegcAACAFHQOAABACjoHAAAgBZ0DoJWsWrVKoVBIq1atau2mAEAKOgdAFnnttdc0c+ZMvfPOO63dFACHMToHQBZ57bXXdPPNN9M5ANCs6BwAhynf93Xw4MHWbgaALETnAEjTI488olAopNWrVx9Sds899ygUCmnTpk2SpNdff13/9m//pq5duyo3N1fDhw/XE088Eej4S5Ys0SWXXCJJOuussxQKhVLuWTjqqKM0duxYLV++XMOHD1deXp7uuecevfPOOwqFQlqyZMkh+wyFQpo5c2bKa9u3b9f3vvc9FRYWKh6P64QTTtB9990XqO0AsktOazcAyBbnnXeejjjiCD300EM644wzUsoefPBBnXDCCTrxxBP16quvatSoUTryyCP105/+VB07dtRDDz2kCy+8UI8++qguuuiiRh3/q1/9qq6//nrddddd+tnPfqYBAwZIUt1/JWnz5s0aP368fvCDH+jqq69Wv379MjrGrl27dOqppyoUCunaa69V9+7d9dRTT2nSpEmqqKjQlClTGtV2AFnGB5C28ePH+z169PBra2vrXtuxY4cfDof9WbNm+b7v++ecc44/aNAgv6qqqm4bz/P80047zT/uuOPqXnv22Wd9Sf6zzz6b9vEffvjhBuv07dvXl+Q//fTTKa9v3brVl+QvXrz4kDqS/BkzZtQ9nzRpkt+zZ0//gw8+SNnuW9/6ll9QUOAfOHAg7bYCyF5MKwAZuOyyy7R79+6U8MNHHnlEnufpsssu00cffaS//e1vuvTSS7Vv3z598MEH+uCDD/Thhx9qzJgxevPNN7V9+/Zma19JSYnGjBnTqLq+7+vRRx/V+eefL9/369r+wQcfaMyYMdq7d69efvnlJm4xgLaIaQUgA9/4xjdUUFCgBx98UOecc46kT6YUTjrpJB1//PFav369fN/X9OnTNX369Hr3sXv3bh155JHN0r6SkpJG133//fe1Z88eLVy4UAsXLqx3m927dzd6/wCyB50DIAPxeFwXXnihli1bprvvvlu7du3S//zP/2j27NmSJM/zJEk33nhjg3/BH3vssc3Wvry8vENeC4VC9W6bTCZTnn/W9ssvv1wTJ06st87gwYMDthBANqBzAGTosssu0+9//3utXLlS//znP+X7vi677DJJ0tFHHy1JikajGj16dJMfu6F/6F26dOkiSdqzZ0/K6//6179Snnfv3l2dOnVSMplslrYDyB7ccwBkaPTo0eratasefPBBPfjggzrllFPqhvN79OihM888U/fcc4927NhxSN33338/0LE7duwo6dB/6F3y8/P1pS99SWvWrEl5/e677055HolENG7cOD366KN1IZmfF7TtALIHIwdAhqLRqC6++GKVlZWpsrJSv/3tb1PK582bp9NPP12DBg3S1VdfraOPPlq7du3S2rVr9e677+p///d/G33sk046SZFIRLfeeqv27t2reDyus88+Wz169HDWu+qqq3TLLbfoqquu0vDhw7VmzRq98cYbh2x3yy236Nlnn9WIESN09dVXa+DAgfroo4/08ssv65lnntFHH33U6LYDyCKtGywBZKcVK1b4kvxQKOSXl5cfUv7222/7V1xxhV9UVORHo1H/yCOP9MeOHes/8sgjdds0JpTR933/3nvv9Y8++mg/Eomk1O/bt69/3nnn1VvnwIED/qRJk/yCggK/U6dO/qWXXurv3r37kFBG3/f9Xbt2+ZMnT/aLi4v9aDTqFxUV+eecc46/cOHCjNoJIHuFfN/3W7d7AgAA2hLuOQAAACm45wBoAw4ePKi9e/c6t+natatisVgLtQhAe0bnAGgDHnzwQV155ZXObZ599lmdeeaZLdMgAO0a9xwAbcCOHTv06quvOrcZNmxY3ZoFANCc6BwAAIAU3JAIAABStLl7DjzP03vvvadOnTo1aqlYAED74fu+9u3bp169eikcbr6/d6uqqlRTUxN4P7FYTLm5uU3QoubV5joH7733noqLi1u7GQCALFJeXq7evXs3y76rqqpU0vcI7dydtDc2FBUVaevWrW2+g9DmOgedOnWSJPWeeZPCbfzkAQBal1dVpXdn/rLu347mUFNTo527k9q6oa/yOzV+dKJin6eSYf9STU0NnYNMfTaVEM7NpXMAAEhLS0xD53cKB+ocZJM21zkAAKAtSvqekgHi+5K+13SNaWZ0DgAASIMnX54a3zsIUrel0TkAACANnjwF+ds/WO2W1T4mTwAAQNoYOQAAIA1J31cywKLCQeq2NDoHAACkoT3dc8C0AgAASMHIAQAAafDkK9lORg7oHAAAkAamFQAAQLvFyAEAAGkgWgEAAKTwPn0EqZ8tmFYAAAApGDkAACANyYDRCkHqtjQ6BwAApCHpK2BWxqZrS3OjcwAAQBq45wAAALRbjBwAAJAGTyElFQpUP1vQOQAAIA2e/8kjSP1swbQCAABIwcgBAABpSAacVghSt6XROQAAIA3tqXPAtAIAAEjByAEAAGnw/JA8P0C0QoC6LY2RAwAA0vDZtEKQR6a2b9+uyy+/XN26dVNeXp4GDRqkl156qcHtn3vuOY0aNapu+/79++uOO+7I+LiMHAAA0AZ9/PHHGjVqlM466yw99dRT6t69u95880116dKlwTodO3bUtddeq8GDB6tjx4567rnn9IMf/EAdO3bU97///bSPTeegHiFrjUvP3fsLJY39W7GujnKzbca+7WMb7y1gnG7gUTWjAeb+jbGy5q7vW2N14Ybfnx+x9m18OFbbA9Y3Wfu3zn3Qiy/otRfwu2PXD7B/69Qav1myfrOs+lkUvx9EUmElA3wRPjvNFRUVKa/H43HF4/FDtr/11ltVXFysxYsX171WUlLiPMbQoUM1dOjQuudHHXWUHnvsMf3973/PqHPAtAIAAGnwP73noLEP/9MOXnFxsQoKCuoec+bMqfd4TzzxhIYPH65LLrlEPXr00NChQ3Xvvfdm1OZ//OMfev7553XGGWdkVI+RAwAA0tBUoYzl5eXKz8+ve72+UQNJ2rJli+bPn69p06bpZz/7mV588UVdf/31isVimjhxovNYvXv31vvvv6/a2lrNnDlTV111VUZtpXMAAEALys/PT+kcNMTzPA0fPlyzZ8+W9MmUwaZNm7RgwQKzc/D3v/9d+/fv17p16/TTn/5Uxx57rMaPH592G+kcAACQhqQfVtK8cchVP7Pte/bsqYEDB6a8NmDAAD366KNm3c/uTRg0aJB27dqlmTNn0jkAAKCpeQrJC3CrnpfhnZujRo3S5s2bU15744031Ldv38yO63mqrq7OqA6dAwAA2qCpU6fqtNNO0+zZs3XppZdq/fr1WrhwoRYuXFi3TWlpqbZv3677779fkjRv3jz16dNH/fv3lyStWbNGv/3tb3X99ddndGw6B/UxwnbCte7q4YS7PFTb+P2bYZJmSFKw8mYPWTJDBd3nzgz3M8o9q77xjfFyjFDLqFHfFa1mRZMZx1bA8lCO++IIW+Vhd3kkYhzfCGUMW6GSAfnG7j3PffF6xu+KlzTqJxuu79W66/pGuYzfJFm/aY62HU5aOrfCySefrGXLlqm0tFSzZs1SSUmJ5s6dqwkTJtRts2PHDm3btq3uued5Ki0t1datW5WTk6NjjjlGt956q37wgx9kdGw6BwAApCH4PQeZd2DHjh2rsWPHNli+ZMmSlOfXXXedrrvuuoyP80WscwAAAFIwcgAAQBo+uSExQOKlLErZTOcAAIA0eAGXT840WqE1Ma0AAABSMHIAAEAaWuOGxNZC56AeVjigFaoYrnbPK0VqjPqOcuvYkYQRDtbMoZCmgFkLrXLPCjWMGuFkMXf9pFEeirn3n7QyC0Yc9a2Tb4UCxoxQwpj7w4/G3DG88Wiw8ljEffyIEQqZY5SHAw7p1hoXXyLpjoNNGKGO1Qn3xVtT23B5TY372LU17n37smJ8g2V1PFx4CrfoIkitic4BAABpSPohJQPknQ9St6VxzwEAAEjByAEAAGlIBoxWSDKtAADA4cXzw/IC3JDoZdENiUwrAACAFIwcAACQBqYVAABACk/BIg6CRoK3JDoH9bBC0a2Uy9Y6BjkH3eWRqoYbEKk26ta4G2+vgxBwnQRjTs1MuWymRDbWETDWMbDWKUjmusvNeG+j2Gp/yHOcP+sPFmOdg7CxjkE8172IRoe4+8LuFHdfnB2j7vodctzleRF3+2JGLvWI8cW2fvRrjXze1cYiGwdq3fm6KxNxZ/n+RMMX7/6wu67xk6OEkXLZN37zQlmUMwDpoXMAAEAagi+ClD23+dE5AAAgDcGXT86ezkH2tBQAALQIRg4AAEiDp5C8APdXBKnb0ugcAACQhvY0rUDnAACANARf54DOQXYzQhmNiClnymXJHaooSdEDDZfnGHUjB420vNVGWtuEu9wZaidJ1vKgZiijFaro/nJ5cXd5bZ6RdtdIq2tdHL4r5bKksDuaTc5IUiuKMse4rqLuUEYrVLFLrjsgrkvuAWd515i7vFNOlbs84i6PG/nMI8ZnlzROcMIKVTTyfVfUuuNkP6rp6CyPRvKc5S5J47pOJtxhml7ISOmMww6dAwAA0uD5IXlBFkHKopTNdA4AAEiDF3BaIZvWOcielgIAgBbByAEAAGkInrI5e/4ep3MAAEAakgqZN65a9bNF9nRjAABAi2DkAACANDCtkKZbbrlFpaWluuGGGzR37lxJUlVVlX784x+rrKxM1dXVGjNmjO6++24VFhY2RXtbRMhIum2lLTbCrc20y661DKL73QePHDDS1la5y0M1RrmR0jnwOgdRdzx1OGbEY+eZOZ/dxzdTLru/3EbmXDvltaPMN1IOhyPuCzcnx1jnIOq+cAvi7nUOesT3uctj7vIuOZXO8k5h9/E7GguMRIwvtrV6XcJ3X3sVnnsdgo8j7nUM4tYCKg6JpLttB3PcC2xUR9zlnpXHvp1IKtjUgPH1b1Ma3Y158cUXdc8992jw4MEpr0+dOlVPPvmkHn74Ya1evVrvvfeeLr744sANBQAALaNRnYP9+/drwoQJuvfee9WlS5e61/fu3atFixbp9ttv19lnn61hw4Zp8eLFev7557Vu3bomazQAAC3ts2mFII9s0aiWTp48Weedd55Gjx6d8vqGDRuUSCRSXu/fv7/69OmjtWvX1ruv6upqVVRUpDwAAGhrPku8FOSRLTK+56CsrEwvv/yyXnzxxUPKdu7cqVgsps6dO6e8XlhYqJ07d9a7vzlz5ujmm2/OtBkAALQoP2DKZv9wDWUsLy/XDTfcoP/6r/9Sbq47iUi6SktLtXfv3rpHeXl5k+wXAAA0TkYjBxs2bNDu3bv15S9/ue61ZDKpNWvW6D//8z+1fPly1dTUaM+ePSmjB7t27VJRUVG9+4zH44rHjVu8AQBoZUGnBg7baYVzzjlHr7zySsprV155pfr3769///d/V3FxsaLRqFauXKlx48ZJkjZv3qxt27Zp5MiRTdfq5mZF4xnxKJGEkVa5pvFpl61QxZz97nCu0EGjvNrIN10bMBjHCGUMRd2XZCjXnRbXSintB00ZHXPvv7bWGDZsxoiwkPG7E41YoYzuz74g6k6Z3C3qDkXsEXXfT9Qjx13eOexO+dwh7I4RtlI2W6p897W5x+vgLM8NGd8tQ7UjDPdArTsUsSLq/gOs0giDzaLR8GZFVsYGdOrUSSeeeGLKax07dlS3bt3qXp80aZKmTZumrl27Kj8/X9ddd51GjhypU089telaDQAAmk2Tr5B4xx13KBwOa9y4cSmLIAEAkM2SAVM2B6nb0gJ3DlatWpXyPDc3V/PmzdO8efOC7hoAgDajPU0rZE83BgAAtAg6BwAApMFTOPAjU9u3b9fll1+ubt26KS8vT4MGDdJLL73U4PaPPfaYvva1r6l79+7Kz8/XyJEjtXz58oyPS+cAAIA0JP1Q4EcmPv74Y40aNUrRaFRPPfWUXnvtNd12220paQu+aM2aNfra176mv/zlL9qwYYPOOussnX/++frHP/6R0bFJ2VwPKwFZ0KyNZqhjtSOU0cqqaIUqHnCHo/lVRsrIRLBwLDPeLuYOyQp57pMfNkIVw1bWx7iRdbHWKnd/tiEvQKijUTVkXLgxIytjbsSdlTE/x50VsSDHHWrYLbLfWd494g5l7Bp2X7sdw8a14Sy1HfDd5ydqffENVb47THdfsuGF5zrkGG2zzk3Y+tFzF6N53HrrrSouLtbixYvrXispKXHW+SxD8mdmz56txx9/XE8++aSGDh2a9rEZOQAAIA2f3ZAY5CHpkHxC1dX1/1H2xBNPaPjw4brkkkvUo0cPDR06VPfee29mbfY87du3T127ds2oHp0DAADS4AfMyOh/ukJicXGxCgoK6h5z5syp93hbtmzR/Pnzddxxx2n58uW65pprdP311+v3v/992m3+7W9/q/379+vSSy/N6L0yrQAAQBqSCikZYI7ls7rl5eXKz8+ve72hFAKe52n48OGaPXu2JGno0KHatGmTFixYoIkTJ5rHW7p0qW6++WY9/vjj6tGjR0ZtZeQAAIAWlJ+fn/JoqHPQs2dPDRw4MOW1AQMGaNu2beYxysrKdNVVV+mhhx7S6NGjM24jIwcAAKTB84MtZGSkfjnEqFGjtHnz5pTX3njjDfXt29dZ74EHHtD3vvc9lZWV6bzzzsu0mZLoHAAAkJbP7h0IUj8TU6dO1WmnnabZs2fr0ksv1fr167Vw4UItXLiwbpvS0lJt375d999/v6RPphImTpyoO++8UyNGjNDOnTslSXl5eSooKEj72EwrAADQBp188slatmyZHnjgAZ144on6xS9+oblz52rChAl12+zYsSNlmmHhwoWqra3V5MmT1bNnz7rHDTfckNGxGTmojzFsFHSdg1DSPbYUTjR8gFCNsc6BkXLZWsfArzLWQagx1jkwxs1CEXd/NJS0Tp4xpGekfA7XGOW1Rrnx2VnrGFjXTqB1DoxY9ohRHjNSOsfD7muvk7EOQaewe52EzkbK5c5G+zuE3WtYROUu92ScH99ax8D93Uj47nUg9oXznOUdwg3v3/psrHTdYePcmou/tBOeQvIC3JDYmLpjx47V2LFjGyxfsmRJyvMv5jtqLDoHAACkoTGrHH6xfrZgWgEAAKRg5AAAgDS09A2JrYnOAQAAafAUChbKmEVJKrKnGwMAAFoEIwcAAKTBDxit4GfRyAGdg3qYUTvNnNI55AgHtMIgVWvs3Ei5bIUqmqGMvrt9ftJImRw2BrNqjVBO4/2Hat0fjnV+A4UitjLrZynHeHNWuFxu2J022Cw3vhi5RrrvDiF3yuNoyH3tJX3rw3Vf+zUh4/wY5VbK56jj/OeE3XWtz9aKEM6if9Oa1eczKza2fragcwAAQBra0w2J2dNSAADQIhg5AAAgDUwrAACAFK2xfHJrYVoBAACkYOQAAIA0MK0AAABS0DlAswoZawE41wqw6gZlpFwOfHwrltwzyoOcu6ZgrXERcI2MIKxY9bDRuLC5iEMwESMlsjXHGTbeYDjgfG7EWEch3MxhaBHj/EccF0+ElMpoYnQOAABIAyMHAAAgRXvqHBCtAAAAUjByAABAGnwFW6sgm+4MoXMAAEAa2tO0Ap0DAADSQOegnQv8+Rn1fSvmzFUepK4kGeFaoYi73Eq5bIUqhiJGfas8bJQHPj/uYvuzDVY/CCuK0/phsjLGWcOpSWP/SeMWp4Sx/4RxbdVaudCt82NskJB7/1YgaNI8f+7z46pvnfvAy/Zm03g4mgSdAwAA0sDIAQAASNGeOgeEMgIAgBSMHAAAkAbfD8kP8Nd/kLotjc4BAABp8BQKdHNn4BtDWxDTCgAAIAUjB41hdKms5G2+cdb9nIZ7l37UHcoXiho7j0Xd9ZPucK1w2HhzVlZFI1QxFIsZ5e72+8b793OMcL2Iu2dvJuazyoOEOpoZH4OFGtYab67Kc5/7Kt/92R3w4u7yULWzPDeUcJbLq3EWR40wXs+IBa0yQikrPff+rfdf6bnPX7Xj/Nd47uu+1mib52XPX7StqT3dkEjnAACANLSnew6YVgAAACkYOQAAIA1MKwAAgBTtaVqBzgEAAGnwA44cZFPngHsOAABACkYOAABIgy87+6lVP1vQOahPyP0RWimXrVh4z7GOgSQlow3vIBwz1gnINdYJsNYhsFIa19a6y61vjpFy2VrHQHH3+/Nz3fW9uPuS9xznXrI/O3ONi2ZN2WxcV0ase1Wt+9wdTLrL9ybznOV7wh2c5eY6BoaqsPvajAZMuVzlu6/dPcY6Bns89/vf57nP375kboNl1mdjffbmcHcWDYc3J08hhVghEQAAtEeMHAAAkAaiFQAAQArPDynUTtY5YFoBAIA2avv27br88svVrVs35eXladCgQXrppZca3H7Hjh369re/reOPP17hcFhTpkxp1HHpHAAAkAbfD/7IxMcff6xRo0YpGo3qqaee0muvvabbbrtNXbp0abBOdXW1unfvrptuuklDhgxp9HtlWgEAgDS09D0Ht956q4qLi7V48eK610pKSpx1jjrqKN15552SpPvuuy/zRn6KzkE9zKged0STjOypSkbdB/DiDQ/oeHnunYc8d9c0bIUqGimPQ7XucDCza2yFgVopl61QxQ7uUMdkrvvD8+JWmKnx2QUMdQySstlKu5uodb/3g0Yo477ahkPpJOnjREdneW7ICIM1VEXc7evgu1M+R4xgxaQxkGqlrLZCET9MHuEs/7jWff5c5/9Arfu6r0m6P/tk0vrRcxcjMxUVFSnP4/G44vFDQ2GfeOIJjRkzRpdccolWr16tI488Uj/60Y909dVXN3sbmVYAACANn40cBHlIUnFxsQoKCuoec+bMqfd4W7Zs0fz583Xcccdp+fLluuaaa3T99dfr97//fbO/V0YOAABIQ1NFK5SXlys/P7/u9fpGDSTJ8zwNHz5cs2fPliQNHTpUmzZt0oIFCzRx4sRGtyMdjBwAAJCGprohMT8/P+XRUOegZ8+eGjhwYMprAwYM0LZt25r7rdI5AACgLRo1apQ2b96c8tobb7yhvn37NvuxmVYAACANn/z1HyRaIbPtp06dqtNOO02zZ8/WpZdeqvXr12vhwoVauHBh3TalpaXavn277r///rrXNm7cKEnav3+/3n//fW3cuFGxWOyQUQgXOgcAAKShpUMZTz75ZC1btkylpaWaNWuWSkpKNHfuXE2YMKFumx07dhwyzTB06NC6/9+wYYOWLl2qvn376p133kn72HQOAABoo8aOHauxY8c2WL5kyZJDXvOD5JX+VEadg/nz52v+/Pl1vY8TTjhBP//5z3XuuedKkqqqqvTjH/9YZWVlqq6u1pgxY3T33XersLAwcENblJV211rnwIiFT7pDklWb52iAsYiClU46HHU3PlxjrXNgJbY1WOsc5Bgpk42Uy9Y6BrV5RrljjQnJ/uysNS7MlM6OMusuad9a58CIda9MuN/c3hr3OgcxI2Wypcp3n7wKYx2B3FCNszxipGK32OscuM/P3lp3yuYPjXUiPqppuL712VUnjFTlRkpnGddWe+Er2JIP2bRcREY3JPbu3Vu33HKLNmzYoJdeeklnn322LrjgAr366quSPpkfefLJJ/Xwww9r9erVeu+993TxxRc3S8MBAGhJTbXOQTbIaOTg/PPPT3n+q1/9SvPnz9e6devUu3dvLVq0SEuXLtXZZ58tSVq8eLEGDBigdevW6dRTT226VgMAgGbT6FDGZDKpsrIyVVZWauTIkdqwYYMSiYRGjx5dt03//v3Vp08frV27tsH9VFdXq6KiIuUBAECb4zfBI0tk3Dl45ZVXdMQRRygej+uHP/yhli1bpoEDB2rnzp2KxWLq3LlzyvaFhYXauXNng/ubM2dOyjKSxcXFGb8JAACaXdAphSyaVsi4c9CvXz9t3LhRL7zwgq655hpNnDhRr732WqMbUFpaqr1799Y9ysvLG70vAACaS0unbG5NGYcyxmIxHXvssZKkYcOG6cUXX9Sdd96pyy67TDU1NdqzZ0/K6MGuXbtUVFTU4P4aykYFAABaR+B1DjzPU3V1tYYNG6ZoNKqVK1dq3LhxkqTNmzdr27ZtGjlyZOCGtiRr5MezQhmNcLekO+JJCUdYkZlO2kgZHDZC9cK1RihjMmDX1zy3RkrkqBXq6K5vhSrW5hlhqMb+jWg3Mww2SMpm30i7m0i4D36gxt34jyPuUDxLjRHnuT/p/iPhg0jCWZ4bdpeHA074Vhvtt8orat2hmHsTRihkdcP191W7z50ZyljLSvrpaOlFkFpTRp2D0tJSnXvuuerTp4/27dunpUuXatWqVVq+fLkKCgo0adIkTZs2TV27dlV+fr6uu+46jRw5kkgFAED2C3rfwOHaOdi9e7euuOIK7dixQwUFBRo8eLCWL1+ur33ta5KkO+64Q+FwWOPGjUtZBAkAAGSPjDoHixYtcpbn5uZq3rx5mjdvXqBGAQDQ1gS9qfCwviERAIB2qR2tn8xdKAAAIAUjBwAApIFohfbOypxnnDUrc5+d4azhsSfPyFqYjLnHrcJGyFLYCFUMWUkZrWEzKxTTSg5nhGomA2bENEMVjSU5zFDGHOMEBQhllPHZ1hqhjFURo/EGO+uj++Tl5bhDJXNz3KGKOcbFGTbKPePiqzXKa4z3X5V0n98DRmbFA4mG6x+odtdNGNlWfSuUMWAy1sNKFk0NBMG0AgAASMHIAQAAaWBaAQAApGpH0Qp0DgAASEtI5o1TZv3swD0HAAAgBSMHAACkg2kFAACQgs5B+2ambDZi1UMxYwdWrL8jbXHYiLOvrTVSNtda6xhYjXMXh6y0wtaUm7XGRMB1EIysuuY6BWa5sY6ClbLZeX6sHxYjZbNX4z54tbX7pBHnb6T7roy4T040J+ksj4TdwfYR4+ILG+WecXGap9+Ral2Samrd599aJ6LWUd9Kx2199jJ+N8zfBRx26BwAAJAOUjYDAIDPa09ZGYlWAAAAKRg5AAAgHdyQCAAAUrSjew6YVgAAACkYOaiHFS7nG+FsSStkygi3Czv2H3JHe9mhhkZIkpmS2RIwZbMZ5hkw5bNZboUaGmGsVn3PKHfGghpv3ro2/Br3m/eMUMiahLt+IuL+OQmFjVBDozxkxcka5SHr2gk45OsboYxW0h3f+O75jlBS3whFtMJczXJSNkv65BKzLkOrfragcwAAQDq45wAAAKTgngMAANBeMXIAAEA6mFYAAAAp2lHngGkFAACQgpEDAADS0Y5GDugc1MeItzazlzpSLktSyHPvP+koDnyvq3VxtvbFG/QNBss43ezrMFiBzq51GKw1KKxYd3ONDGshgJAx0GiFylvrfzTzZ9/smvm7FXJcXGb8vLVOgbWGRmv/LrQVRCsAAID2ipEDAADSwAqJAAAgVTu654BpBQAA2qjt27fr8ssvV7du3ZSXl6dBgwbppZdectZZtWqVvvzlLysej+vYY4/VkiVLMj4unQMAANqgjz/+WKNGjVI0GtVTTz2l1157Tbfddpu6dOnSYJ2tW7fqvPPO01lnnaWNGzdqypQpuuqqq7R8+fKMjs20AgAAaQgp4D0Hn/63oqIi5fV4PK54PH7I9rfeequKi4u1ePHiutdKSkqcx1iwYIFKSkp02223SZIGDBig5557TnfccYfGjBmTdlvpHNTDTNtrhqsZcUPWeI1j/37QO1oChuo1u6BzcmZW34CxjoHLA6TMNvYdTljHdl945qUVsNy8tJr5s292zfzdcr69gKnMzQ8/6P4PF00UylhcXJzy8owZMzRz5sxDNn/iiSc0ZswYXXLJJVq9erWOPPJI/ehHP9LVV1/d4CHWrl2r0aNHp7w2ZswYTZkyJaOm0jkAAKAFlZeXKz8/v+55faMGkrRlyxbNnz9f06ZN089+9jO9+OKLuv766xWLxTRx4sR66+zcuVOFhYUprxUWFqqiokIHDx5UXl5eWm2kcwAAQDqaKFohPz8/pXPQEM/zNHz4cM2ePVuSNHToUG3atEkLFixosHPQVNrLYBAAAMH4TfDIQM+ePTVw4MCU1wYMGKBt27Y1WKeoqEi7du1KeW3Xrl3Kz89Pe9RAonMAAECbNGrUKG3evDnltTfeeEN9+/ZtsM7IkSO1cuXKlNdWrFihkSNHZnRsOgcAAKThsxUSgzwyMXXqVK1bt06zZ8/WW2+9paVLl2rhwoWaPHly3TalpaW64oor6p7/8Ic/1JYtW/STn/xEr7/+uu6++2499NBDmjp1akbHpnMAAEA6Wnha4eSTT9ayZcv0wAMP6MQTT9QvfvELzZ07VxMmTKjbZseOHSnTDCUlJfrv//5vrVixQkOGDNFtt92m3/3udxmFMUrckFgv38jK6OcYn3AkWHnIsf9wxB0maSbOM7quobCVvs3av7vcDxjq5xvlnpHaz7fKk8YbqDVOsFHfyozoyqxohSqGjKyMYePYVtus/Zv1jUvLLG/uMNOgoYbGpWGHExr1HSHWnvFL7hu/OX6O8b0ywrutTLZovLFjx2rs2LENlte3+uGZZ56pf/zjH4GOS+cAAIB0tKPcCnQOAABIQ3vKysg9BwAAIAUjBwAApKOJlk/OBnQOAABIB/ccAACAz+OeAwAA0G4xclAfq8tkrHMQirkDtsMxd0B4NNpweU6OUTfiLo8Z9SPGOgfWjFnY6Bp7xpxb0ir33B9OotYdkJ1IGuUJd3mtUe7VuMv9GiNtsuvjMQLlrXUMwjXGOgg1Rn1jnQWzvNYqN75X1hIcAddJsFOxG/WN3w3PWEvAN36NXWsZeDGjbtT63hnnPmCm88MG0woAACBFwGmFbOocMK0AAABSMHIAAEA6mFYAAAAp2lHngGkFAACQgpEDAADS0J7WOaBzUA8rZbMVyhgxQhXjue6Yrw7xhmPKOkSNulF3PFpuxF0/ZoRC5hjxZGGj3DPivWqN8qraqLP8oFFemXDHfB2ocdevirjLq52lkmelhHbEjFk/LFbKZCtUMafKqG+8uUi18b0wjh9JuOsHDnUMmLLZDFWMGGmPjd+NZMyo77h0a63rygrjDBthltEs+lcNTYJpBQAAkIKRAwAA0tGObkikcwAAQBq45wAAABwqi/6BD4J7DgAAQApGDgAASAf3HAAAgM/jnoMGzJkzR4899phef/115eXl6bTTTtOtt96qfv361W1TVVWlH//4xyorK1N1dbXGjBmju+++W4WFhU3e+GZjTLaEctwB1dGYOzetax0DSeqSe7DBsoJ4w2WSVBB1B6vn57jrx428ula5xTOCyas8Yx2DpLt8X22us3xvjbv840gHZ7klmXRfPDUJ6+JylFvrHNQaKZmtlMrGOgbRA+4G5Bw0yquNNTKMdRLCCaN+0grmdxdbfHMdA2udAmONj1x3/dpcV33jzRnrGLjSQX+ygVGOw05G9xysXr1akydP1rp167RixQolEgl9/etfV2VlZd02U6dO1ZNPPqmHH35Yq1ev1nvvvaeLL764yRsOAECL8pvgkSUyGjl4+umnU54vWbJEPXr00IYNG/TVr35Ve/fu1aJFi7R06VKdffbZkqTFixdrwIABWrdunU499dSmazkAAC2oPU0rBIpW2Lt3rySpa9eukqQNGzYokUho9OjRddv0799fffr00dq1a+vdR3V1tSoqKlIeAACg9TS6c+B5nqZMmaJRo0bpxBNPlCTt3LlTsVhMnTt3Ttm2sLBQO3furHc/c+bMUUFBQd2juLi4sU0CAKD5tKNphUZ3DiZPnqxNmzaprKwsUANKS0u1d+/eukd5eXmg/QEA0CzaUeegUaGM1157rf785z9rzZo16t27d93rRUVFqqmp0Z49e1JGD3bt2qWioqJ69xWPxxWPxxvTDAAA0Awy6hz4vq/rrrtOy5Yt06pVq1RSUpJSPmzYMEWjUa1cuVLjxo2TJG3evFnbtm3TyJEjm67VrSxshDLGo+5wv05xd8xYl9wDDZb1iO9z1u0WrXSWF+Q0vG9J6hR2h0LmWvFwhqRvhDL67pTKe5N5zvKPEx2d5bGAoZiJZMRZXlPr/kolIsZXznV6gqZsNj46K+WyFaoYrXR/L3IOuhsYqXKXh6vdn12o1oi384P92ebnuAdavZj7s/Xi7msnVGuUO5rvh422GZdd2Pj7LJtupGtO7emGxIw6B5MnT9bSpUv1+OOPq1OnTnX3ERQUFCgvL08FBQWaNGmSpk2bpq5duyo/P1/XXXedRo4cSaQCACC7sUJi/ebPny9JOvPMM1NeX7x4sb773e9Kku644w6Fw2GNGzcuZREkAACyGp2D+vlpDMvl5uZq3rx5mjdvXqMbBQAAWg+5FQAASAP3HAAAgFTtaFoh0AqJAADg8MPIQX3CRna4cLBQxo5Rd1bGrjFHKGPMHcrYI+pefrpbZL+zvFPYnbXRCmWMGOnbkkZ/9IDnjqnaE3ZnTcwNBQtVrDFivioT7vZVRtyhmCHj2vIc447uIFApZETyWVGcEfdlaWZVtEIVcyrd1074gLsBoSp3/VDCeIPWPVMh9xn2c9yhhuGYO2Ool+cuD/nucoUbvjaTUXfbk0YYayhpXF1kZZTEtAIAAPgiphUAAEB7xcgBAADpYOQAAAB8XqgJHpmYOXOmQqFQyqN///4Nbp9IJDRr1iwdc8wxys3N1ZAhQ/T0009neNRPMHIAAEAbdcIJJ+iZZ56pe56T0/A/2zfddJP++Mc/6t5771X//v21fPlyXXTRRXr++ec1dOjQjI5L5wAAgHS0wrRCTk5Og1mNv+gPf/iD/t//+3/65je/KUm65ppr9Mwzz+i2227TH//4x4yOy7QCAABp+CyUMchDkioqKlIe1dUNZ+p988031atXLx199NGaMGGCtm3b1uC21dXVys3NTXktLy9Pzz33XMbvlZGD+hgTQ5GIu/sXi7jjvTvkuOO5O+U0nDa5S447JXOPHPc6B90j7vLOYXc66VwjL7DV20wYJ/dAyDp+sJTRVb77kt+fdK9jkJfjXmchmmOcH3OdA0ehFaZvrnPg3kEkYazvYaR0NlMuG+sYhPe704Wr2l3frzGuDc/IaW2scxByDOdKknLd107YWGfBD7uPH442/O2KJNzfvHCt8d6sU5NFN9I1qyYaOSguLk55ecaMGZo5c+Yhm48YMUJLlixRv379tGPHDt188836yle+ok2bNqlTp06HbD9mzBjdfvvt+upXv6pjjjlGK1eu1GOPPaZk0viA60HnAACAFlReXq78/Py65/F4/R3Lc889t+7/Bw8erBEjRqhv37566KGHNGnSpEO2v/POO3X11Verf//+CoVCOuaYY3TllVfqvvvuy7iNTCsAAJAuP8DjU/n5+SmPhjoHX9S5c2cdf/zxeuutt+ot7969u/70pz+psrJS//rXv/T666/riCOO0NFHH53x26RzAABAGprqnoPG2r9/v95++2317NnTuV1ubq6OPPJI1dbW6tFHH9UFF1yQ8bHoHAAA0AbdeOONWr16td555x09//zzuuiiixSJRDR+/HhJ0hVXXKHS0tK67V944QU99thj2rJli/7+97/rG9/4hjzP009+8pOMj809BwAApKOFQxnfffddjR8/Xh9++KG6d++u008/XevWrVP37t0lSdu2bVM4/H9/41dVVemmm27Sli1bdMQRR+ib3/ym/vCHP6hz584ZN5XOAQAAaWjprIxlZWXO8lWrVqU8P+OMM/Taa69l2Kr60Tmoj/EJhozyiJHSOS/iDrnqFGk4pMtKqdw53HC6Z0nqGnaHi3U22p4bMkKmjHCwhG/tP2CoYsSd9rbCy3OWf2B8Nrk5Rspq4/xZ104Q1q6DhjqGE+4dhKvdKZOtlMtmqOIB97Xv1xg5p61wrrAxy2qFMhqhiiFj/+GYe/8Rx/m3PjvrszcvS99a+JdYx8MNnQMAANLRjhIv0TkAACANLT2t0JqIVgAAACkYOQAAIB1MKwAAgBR0DgAAwOdxzwEAAGi3GDloBCvtbo4R6x4Lu+PB4+GG48E7ht2x3B2MlMsdjbZ1CEfc5aGYszxspGSutXLDeu73V2Wcuw6+lfLZvf9cx7mXpBwjYDxiLjZglTvKzFh0o9yKdbfWQUgasfS1xhoPCfdnZ6VcttYx8Kvdn71vrXNgreERM95fxP3dUcxYByLhbl/IsZZByL1r87O1rg18imkFAADweSHfV8hY7Mqqny2YVgAAACkYOQAAIB1MKwAAgM8jWgEAALRbjBwAAJAOphXaOSs7qSFsXAFWuFvEUT9ihdIZx7aGiqJyh2NFQ0a4lsX4ckSNcLKo3OFeESMmyzr31mcXNs5/2Ez37SwOJmAUZeBQSetObKvcM0INjVBEM1TRbJ+1f/e1HzLSkQc+P45y8y74gP8oZdNweHNiWgEAALRbjBwAAJAOphUAAMDntadpBToHAACkox2NHHDPAQAASMHIAQAAacqmqYEg6By0gqTvjmdLOmIpk37zDvZ4Rihg0gjXihihiJ4xruYZIVlW8rhkMw+Gecb594zPtlnzrhhhkkbTAofwmqw4Tqs8bHy2xrVnhSrax2/mExTg/PjNGiOLOr4f7EtM4iUAAJCtGDkAACANRCsAAIBURCsAAID2ipEDAADSEPI+eQSpny3oHAAAkA6mFQAAQHvFyEF9AmZWrTVi4Ws9d+rXhNfwx5Lw3XWrfPdHesBPOMtjViy4apylYeO9J4yUy1XGOgpV1vv3ooHKqx3nXrI/2zb9h4G1DoK1TEDEWMMhxzg3OUbK4xzj58goD8fc146VctlaxyAUi7nrR41y4/2b59fRPuNrYX62gdfIaCeIVgAAAKna0SJIdA4AAEhDexo54J4DAACQgpEDAADS0Y6iFegcAACQBqYVAABAu8XIQX2MuB3PM8L1jJApK1zugNdwSFSFl+esu8fr4CyPhoKFKtaEao36btYCYZXGud3jxZ3l+4zzs8/LdZZbn02N8dkmjfb7Rrlz2DFoKKJR7hmhdF6OUR6zQg3dYaTKdX+21p3eoYgRKmmEyZqMUMWQ0X4/7q7vG+fPizb8AZphkIFDHbPoT97mRLQCAAD4PKYVAABAu0XnAACAdPhN8MjAzJkzFQqFUh79+/d31pk7d6769eunvLw8FRcXa+rUqaqqqsrswGJaAQCAtLTGtMIJJ5ygZ555pu55jmMZ8aVLl+qnP/2p7rvvPp122ml644039N3vflehUEi33357RselcwAAQBuVk5OjoqKitLZ9/vnnNWrUKH3729+WJB111FEaP368XnjhhYyPy7QCAADp8PzgD0kVFRUpj+rq6gYP+eabb6pXr146+uijNWHCBG3btq3BbU877TRt2LBB69evlyRt2bJFf/nLX/TNb34z47dK5wAAgHQ00T0HxcXFKigoqHvMmTOn3sONGDFCS5Ys0dNPP6358+dr69at+spXvqJ9+/bVu/23v/1tzZo1S6effrqi0aiOOeYYnXnmmfrZz36W8VvNeFphzZo1+s1vfqMNGzZox44dWrZsmS688MK6ct/3NWPGDN17773as2ePRo0apfnz5+u4447LuHGtxpgX8jx3THHCiGU/UOuO966obTgW/+NIR2fd3JB7nQJLwj9g7D/YOgdJI1j/gLGOgbWOw4fJI5zle2vd9Stq3eskVCXdn11NrZGW18p924zrHJjrFOS4L3wv5j6AF3e/dy/Pfe7C1joGYeMNxoxrM2iMuZVy2VrHIM9dnsx1/xwn4w2//2S0edc5IGXzJ0IKeM/Bp/8tLy9Xfn5+3evxeP2/e+eee27d/w8ePFgjRoxQ37599dBDD2nSpEmHbL9q1SrNnj1bd999t0aMGKG33npLN9xwg37xi19o+vTpGbU145GDyspKDRkyRPPmzau3/Ne//rXuuusuLViwQC+88II6duyoMWPGNOpuSQAADjf5+fkpj4Y6B1/UuXNnHX/88XrrrbfqLZ8+fbq+853v6KqrrtKgQYN00UUXafbs2ZozZ448L7NFwDIeOTj33HNTejOf5/u+5s6dq5tuukkXXHCBJOn+++9XYWGh/vSnP+lb3/pWpocDAKBtaOUVEvfv36+3335b3/nOd+otP3DggMJfGGGLfLpyqJ/hsZv0noOtW7dq586dGj16dN1rBQUFGjFihNauXVtvnerq6kNuzgAAoK35LJQxyCMTN954o1avXq133nlHzz//vC666CJFIhGNHz9eknTFFVeotLS0bvvzzz9f8+fPV1lZmbZu3aoVK1Zo+vTpOv/88+s6Celq0lDGnTt3SpIKCwtTXi8sLKwr+6I5c+bo5ptvbspmAACQ9d59912NHz9eH374obp3767TTz9d69atU/fu3SVJ27ZtSxkpuOmmmxQKhXTTTTdp+/bt6t69u84//3z96le/yvjYrb7OQWlpqaZNm1b3vKKiQsXFxa3YIgAA6tGIVQ4PqZ+BsrIyZ/mqVatSnufk5GjGjBmaMWNGhg07VJN2Dj5bqGHXrl3q2bNn3eu7du3SSSedVG+deDye9s0YAAC0lpDvKxTgvoEgdVtak3YOSkpKVFRUpJUrV9Z1BioqKvTCCy/ommuuacpDNS8rlDHpvlWjOuE+rZUJd2foo5qGwxXj4WChhFW+O5xqX9gdymelfI6E3HfEJo2YqUpHumrJTsn8ca071PPDhLt8b8Kd0vlAwt0+K1130KzBTlaoo/FtT8bcO6jNdZeHjDDOkO8OZfTD7v2HjZTGoYSRjtz6YQ5Z4YBGudE+M1Sxg1Ge5whlND67oKGO1rWFw0/GnYP9+/enhFFs3bpVGzduVNeuXdWnTx9NmTJFv/zlL3XccceppKRE06dPV69evVLWQgAAIOt4nz6C1M8SGXcOXnrpJZ111ll1zz+7X2DixIlasmSJfvKTn6iyslLf//73tWfPHp1++ul6+umnlZvr/osMAIC2jGkFhzPPPNMZLxkKhTRr1izNmjUrUMMAAEDraPVoBQAAskILRyu0JjoHAACko5VXSGxJdA4AAEhDY1Y5/GL9bEHnoD5GCjIv6S6vqXWf1v1GOFw04g7Xc6n23Mfel3TfGNoh7M7qGDVCKSPGuJmVlbHac4e7We3f58hoKUkf1bizMu6tdp/7Awl3+2qtrIxGGGzIce1ZvytWOJpxaciIIlVtrtV2d32F3Q0IR937jyTct3qHao0GBA1lNEItPaP9rqyKkjtUUXKHkiaNpWKMr1XgrI04/NA5AAAgHUwrAACAzwt5nzyC1M8WDBYBAIAUjBwAAJAOphUAAECKdrTOAdMKAAAgBSMHAACkgdwK7Z2VsrnWPeBSU+MOGt4fNoKSHayUwAdq3QHNHXISznIrJXRO2ErZbKxzYKwhUWME4x9Mut/fgVp3sH6lscbEvmr3Z3Og2kjZnDDWOag10h67Tp+RNjfwOgbG+h3WF8MPG3H8RtrgSMJdP2ysYxAysplbP8y+mbLZvX8rpbP1/s20y45LM2mk07Y+e+vaIWXzp9rRPQdMKwAAgBSMHAAAkA5fUpC1CrJn4IDOAQAA6eCeAwAAkMpXwHsOmqwlzY57DgAAQApGDgAASEc7ilagc1CPkGeENBmhjLU1Rjiecfyk1/D+D+a4Q/kqou5QvGjYfTdNNOIOVcwJmDnEM2Kiah3vXXKfG0mqMUI9qxPuz8YqTxifrWeEscoKF3ScXittrh9x//B4RiidOeRppSw2fk2S7ihaha0wT+PSMy/NgL/L9vkPVm6FOrrSLluhikmj3M+xwjzd9dsNT8HCOkm8BAAAshUjBwAApIFoBQAAkKod3XPAtAIAAEjByAEAAOloRyMHdA4AAEhHO+ocMK0AAABSMHJQH3eov2TEY/tyBzQnjFj3pCPtb3XEvc5BZcQdSBsOu3uuYWMdBCOrbWCetcaEEXCdNM6tZ6yTYKXjtta4sK4Nc50D1/sz0mH7Oca5sVIWB1zHwMpEHjLee8j43hlvv/ljyK2PzloHoRnXSbA+G2sdA9caCtax25V2tM4BnQMAANJAKCMAAEjFPQcAAKC9YuQAAIB0eH4aN78Y9bMEnQMAANLBtAIAAGivGDmoh5WyWUbqWZkpn42QsVDDcUOeNaRlhdmY5QH3bwnacbZyx5rhbgHfgJk2uPEpmSXj9Fu7NsLNrDBUP2qcvCBtb4r6AfMGBxkNTuvwVqhpwO+mMxQyQISslEaoohEC3X4EHDkI/APYcugcAACQDqYVAABAe0XnAACAdHh+8EcGZs6cqVAolPLo379/g9ufeeaZh2wfCoV03nnnZfxWmVYAACAdvvfJI0j9DJ1wwgl65pln6p7n5DT8z/Zjjz2mmpqauucffvihhgwZoksuuSTj49I5AACgjcrJyVFRUVFa23bt2jXleVlZmTp06NCozgHTCgAApOOzGxKDPCRVVFSkPKqrqxs85JtvvqlevXrp6KOP1oQJE7Rt27a0m7to0SJ961vfUseOHTN+q4wc1MeK5rMy61nZ5QLHA6I9srL6WeFm2XOfdEOy/x0gy3m+Al2Hn95zUFxcnPLyjBkzNHPmzEM2HzFihJYsWaJ+/fppx44duvnmm/WVr3xFmzZtUqdOnZyHWr9+vTZt2qRFixY1qql0DgAASEcThTKWl5crPz+/7uV4vP585+eee27d/w8ePFgjRoxQ37599dBDD2nSpEnOQy1atEiDBg3SKaec0qim0jkAAKAF5efnp3QO0tW5c2cdf/zxeuutt5zbVVZWqqysTLNmzWpsE7nnAACAtPgKeM9BsMPv379fb7/9tnr27Onc7uGHH1Z1dbUuv/zyRh+LzgEAAOloohsS03XjjTdq9erVeuedd/T888/roosuUiQS0fjx4yVJV1xxhUpLSw+pt2jRIl144YXq1q1bo98q0woAALRB7777rsaPH68PP/xQ3bt31+mnn65169ape/fukqRt27YpHE79G3/z5s167rnn9Ne//jXQsekcAACQDs+TmUHMrJ++srIyZ/mqVasOea1fv37ymyCHA50DAADSQeIlAADQXjFyAABAOtrRyAGdAwAA0tFEKyRmA6YVAABACkYOAABIg+978gOkbA5St6XROQAAIB2+H2xqgHsOAAA4zPgB7znIos4B9xwAAIAUjBwAAJAOz5NCAe4b4J4DAAAOM0wrAACA9oqRAwAA0uB7nvwA0wqEMgIAcLhhWgEAALRXjBwAAJAOz5dCjBwEMm/ePB111FHKzc3ViBEjtH79+uY6FAAAzc/3PwlHbPSjnXcOHnzwQU2bNk0zZszQyy+/rCFDhmjMmDHavXt3cxwOAAA0oWbpHNx+++26+uqrdeWVV2rgwIFasGCBOnTooPvuu685DgcAQLPzPT/wI1s0eeegpqZGGzZs0OjRo//vIOGwRo8erbVr1x6yfXV1tSoqKlIeAAC0OYGmFLysWiGxyTsHH3zwgZLJpAoLC1NeLyws1M6dOw/Zfs6cOSooKKh7FBcXN3WTAAAIjJGDFlRaWqq9e/fWPcrLy1u7SQAAtGtNHsr4pS99SZFIRLt27Up5fdeuXSoqKjpk+3g8rng8Xvfc//RuTq+qqqmbBgA4zHz2b4XfApEAtX51oKmBWiWasDXNq8k7B7FYTMOGDdPKlSt14YUXSpI8z9PKlSt17bXXmvX37dsnSXp35i+bumkAgMPUvn37VFBQ0Cz7jsViKioq0nM7/xJ4X0VFRYrFYk3QqubVLIsgTZs2TRMnTtTw4cN1yimnaO7cuaqsrNSVV15p1u3Vq5fKy8vVqVMnhUIhVVRUqLi4WOXl5crPz2+O5h7WOH+Nx7lrPM5dMJy/9Pm+r3379qlXr17Ndozc3Fxt3bpVNTU1gfcVi8WUm5vbBK1qXs3SObjsssv0/vvv6+c//7l27typk046SU8//fQhNynWJxwOq3fv3oe8np+fz5ckAM5f43HuGo9zFwznLz3NNWLwebm5uVnxj3pTabblk6+99tq0phEAAEDb0urRCgAAoG1p852DeDyuGTNmpEQ0IH2cv8bj3DUe5y4Yzh9aW8hvifgPAACQNdr8yAEAAGhZdA4AAEAKOgcAACAFnQMAAJCCzgEAAEjR5jsH8+bN01FHHaXc3FyNGDFC69evb+0mtTlr1qzR+eefr169eikUCulPf/pTSrnv+/r5z3+unj17Ki8vT6NHj9abb77ZOo1tY+bMmaOTTz5ZnTp1Uo8ePXThhRdq8+bNKdtUVVVp8uTJ6tatm4444giNGzfukMRi7dX8+fM1ePDgupX8Ro4cqaeeeqqunHOXvltuuUWhUEhTpkype43zh9bSpjsHDz74oKZNm6YZM2bo5Zdf1pAhQzRmzBjt3r27tZvWplRWVmrIkCGaN29eveW//vWvddddd2nBggV64YUX1LFjR40ZM0ZVZL7U6tWrNXnyZK1bt04rVqxQIpHQ17/+dVVWVtZtM3XqVD355JN6+OGHtXr1ar333nu6+OKLW7HVbUfv3r11yy23aMOGDXrppZd09tln64ILLtCrr74qiXOXrhdffFH33HOPBg8enPI65w+txm/DTjnlFH/y5Ml1z5PJpN+rVy9/zpw5rdiqtk2Sv2zZsrrnnuf5RUVF/m9+85u61/bs2ePH43H/gQceaIUWtm27d+/2JfmrV6/2ff+TcxWNRv2HH364bpt//vOfviR/7dq1rdXMNq1Lly7+7373O85dmvbt2+cfd9xx/ooVK/wzzjjDv+GGG3zf59pD62qzIwc1NTXasGGDRo8eXfdaOBzW6NGjtXbt2lZsWXbZunWrdu7cmXIeCwoKNGLECM5jPfbu3StJ6tq1qyRpw4YNSiQSKeevf//+6tOnD+fvC5LJpMrKylRZWamRI0dy7tI0efJknXfeeSnnSeLaQ+tqtsRLQX3wwQdKJpOHZHIsLCzU66+/3kqtyj47d+6UpHrP42dl+ITneZoyZYpGjRqlE088UdIn5y8Wi6lz584p23L+/s8rr7yikSNHqqqqSkcccYSWLVumgQMHauPGjZw7Q1lZmV5++WW9+OKLh5Rx7aE1tdnOAdDSJk+erE2bNum5555r7aZklX79+mnjxo3au3evHnnkEU2cOFGrV69u7Wa1eeXl5brhhhu0YsWKdpUKGNmhzU4rfOlLX1IkEjnkztxdu3apqKiolVqVfT47V5xHt2uvvVZ//vOf9eyzz6p37951rxcVFammpkZ79uxJ2Z7z939isZiOPfZYDRs2THPmzNGQIUN05513cu4MGzZs0O7du/XlL39ZOTk5ysnJ0erVq3XXXXcpJydHhYWFnD+0mjbbOYjFYho2bJhWrlxZ95rneVq5cqVGjhzZii3LLiUlJSoqKko5jxUVFXrhhRc4j/okzPPaa6/VsmXL9Le//U0lJSUp5cOGDVM0Gk05f5s3b9a2bds4fw3wPE/V1dWcO8M555yjV155RRs3bqx7DB8+XBMmTKj7f84fWkubnlaYNm2aJk6cqOHDh+uUU07R3LlzVVlZqSuvvLK1m9am7N+/X2+99Vbd861bt2rjxo3q2rWr+vTpoylTpuiXv/yljjvuOJWUlGj69Onq1auXLrzwwtZrdBsxefJkLV26VI8//rg6depUN5dbUFCgvLw8FRQUaNKkSZo2bZq6du2q/Px8XXfddRo5cqROPfXUVm596ystLdW5556rPn36aN++fVq6dKlWrVql5cuXc+4MnTp1qru35TMdO3ZUt27d6l7n/KHVtHa4hOU//uM//D59+vixWMw/5ZRT/HXr1rV2k9qcZ5991pd0yGPixIm+738Szjh9+nS/sLDQj8fj/jnnnONv3ry5dRvdRtR33iT5ixcvrtvm4MGD/o9+9CO/S5cufocOHfyLLrrI37FjR+s1ug353ve+5/ft29ePxWJ+9+7d/XPOOcf/61//WlfOucvM50MZfZ/zh9YT8n3fb6V+CQAAaIPa7D0HAACgddA5AAAAKegcAACAFHQOAABACjoHAAAgBZ0DAACQgs4BAABIQecAAACkoHMAAABS0DkAAAAp6BwAAIAU/x+pUUriTDVaGgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# final model\n", "fpath_final_model = './OUTPUT_FILES/final_model.h5'\n", "fpath_true_model = \"./test_model_true.h5\"\n", "\n", "import h5py\n", "\n", "# read vel, xi, eta dataset \n", "with h5py.File(fpath_final_model, 'r') as f:\n", " vel = f['vel'][:]\n", " xi = f['xi'][:]\n", " eta = f['eta'][:]\n", " #zeta = f['zeta'][:]\n", "\n", "# read true model\n", "with h5py.File(fpath_true_model, 'r') as f:\n", " vel_true = f['vel'][:]\n", " xi_true = f['xi'][:]\n", " eta_true = f['eta'][:]\n", " #zeta_true = f['zeta'][:]\n", "# plot\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "# plot vel\n", "plt.figure()\n", "plt.imshow(vel[5,:,:], origin='lower', aspect='auto')\n", "plt.colorbar()\n", "plt.title('vel')\n", "\n", "# plot xi\n", "#plt.figure()\n", "#plt.imshow(xi[:,:,10], origin='lower', aspect='auto')\n", "#plt.colorbar()\n", "#plt.title('xi')\n", "#\n", "## plot eta\n", "#plt.figure()\n", "#plt.imshow(eta[:,:,10], origin='lower', aspect='auto')\n", "#plt.colorbar()\n", "#plt.title('eta')\n", "\n", "# plot vel_true\n", "plt.figure()\n", "plt.imshow(vel_true[5,:,:], origin='lower', aspect='auto')\n", "plt.colorbar()\n", "plt.title('vel_true')\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.1 64-bit ('3.9.1')", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.8" }, "vscode": { "interpreter": { "hash": "02f83e1f4cd9619657a6845405e2dd67c4de23753ba48bca5dce2ebf57b3813a" } } }, "nbformat": 4, "nbformat_minor": 2 }