{ "cells": [ { "cell_type": "markdown", "id": "cb353fdd", "metadata": { "papermill": { "duration": 0.0122, "end_time": "2022-10-20T19:14:09.862731", "exception": false, "start_time": "2022-10-20T19:14:09.850531", "status": "completed" }, "tags": [] }, "source": [ "# Breakdown the PIV into separate functions and details" ] }, { "cell_type": "code", "execution_count": 1, "id": "bcf33c44", "metadata": { "execution": { "iopub.execute_input": "2022-10-20T19:14:09.884467Z", "iopub.status.busy": "2022-10-20T19:14:09.883874Z", "iopub.status.idle": "2022-10-20T19:14:09.893190Z", "shell.execute_reply": "2022-10-20T19:14:09.891809Z" }, "papermill": { "duration": 0.024058, "end_time": "2022-10-20T19:14:09.896292", "exception": false, "start_time": "2022-10-20T19:14:09.872234", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "# test the idea of vectorized cross correlation for \n", "# strided images, rectangular windows and extended search area\n", "# in one function" ] }, { "cell_type": "code", "execution_count": 2, "id": "c3b039ba", "metadata": { "execution": { "iopub.execute_input": "2022-10-20T19:14:09.905785Z", "iopub.status.busy": "2022-10-20T19:14:09.905490Z", "iopub.status.idle": "2022-10-20T19:14:10.644236Z", "shell.execute_reply": "2022-10-20T19:14:10.643225Z" }, "papermill": { "duration": 0.746244, "end_time": "2022-10-20T19:14:10.646814", "exception": false, "start_time": "2022-10-20T19:14:09.900570", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "from scipy.fft import rfft2, irfft2, fftshift" ] }, { "cell_type": "code", "execution_count": 3, "id": "54beb307", "metadata": { "execution": { "iopub.execute_input": "2022-10-20T19:14:10.659964Z", "iopub.status.busy": "2022-10-20T19:14:10.659390Z", "iopub.status.idle": "2022-10-20T19:14:10.941134Z", "shell.execute_reply": "2022-10-20T19:14:10.940500Z" }, "papermill": { "duration": 0.290236, "end_time": "2022-10-20T19:14:10.942644", "exception": false, "start_time": "2022-10-20T19:14:10.652408", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "from openpiv.pyprocess import moving_window_array, \\\n", " normalize_intensity,\\\n", " find_subpixel_peak_position, \\\n", " get_field_shape, \\\n", " get_coordinates,\\\n", " correlation_to_displacement,\\\n", " fft_correlate_images\n", "from openpiv.tools import imread, transform_coordinates" ] }, { "cell_type": "code", "execution_count": 4, "id": "9a21fa0a", "metadata": { "execution": { "iopub.execute_input": "2022-10-20T19:14:10.950468Z", "iopub.status.busy": "2022-10-20T19:14:10.950187Z", "iopub.status.idle": "2022-10-20T19:14:11.642135Z", "shell.execute_reply": "2022-10-20T19:14:11.641331Z" }, "papermill": { "duration": 0.697919, "end_time": "2022-10-20T19:14:11.643048", "exception": true, "start_time": "2022-10-20T19:14:10.945129", "status": "failed" }, "tags": [] }, "outputs": [], "source": [ "frame_a = imread('../test1/exp1_001_a.bmp')\n", "frame_b = imread('../test1/exp1_001_b.bmp')\n", "\n", "# frame_a = frame_a[:128,:128]\n", "# frame_b = frame_b[:128,:128]\n", "\n", "# frame_a = normalize_intensity(frame_a)\n", "# frame_b = normalize_intensity(frame_b)" ] }, { "cell_type": "code", "execution_count": 5, "id": "7cd2d11c", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "# for debugging purposes \n", "# frame_a = frame_a[:64,:64]\n", "# frame_b = frame_b[:64,:64]" ] }, { "cell_type": "code", "execution_count": 6, "id": "0679abbf", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "# parameters for the test\n", "window_size = 48\n", "overlap = 8\n", "search_size = window_size #not extended search for a while" ] }, { "cell_type": "code", "execution_count": 7, "id": "248481be", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "# for the regular square windows case:\n", "aa = moving_window_array(frame_a, window_size, overlap)\n", "bb = moving_window_array(frame_b, window_size, overlap)" ] }, { "cell_type": "code", "execution_count": 8, "id": "578ea6a4", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "c = fft_correlate_images(aa,bb)\n", "n_rows, n_cols = get_field_shape(frame_a.shape, search_size, overlap)\n", "u,v = correlation_to_displacement(c, n_rows,n_cols)\n", "x,y = get_coordinates(frame_a.shape,search_size,overlap)" ] }, { "cell_type": "code", "execution_count": 9, "id": "02b9e299", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "# let's assume we want the extended search type of PIV analysis\n", "# with search_area_size in image B > window_size in image A\n", "window_size = 32\n", "overlap = 8\n", "search_size = 48" ] }, { "cell_type": "code", "execution_count": 10, "id": "d9c353ce", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPsAAAD6CAYAAABnLjEDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAbNUlEQVR4nO3de3DW1ZkH8O8DhIRwRxEDAYFCocgAOiBodYsiXXQVUJTZrrXoOEMvroWplqJrHa22Re2otGNFXMqyDAW8zjKM1lUWCsoKIiKGe6ACwUDAJAIikMCzf+Sly7mE95df3ls4388MA8/hvL/f8U0e3/wO5zxHVBVEdP5rlu0BEFFmMNmJAsFkJwoEk50oEEx2okAw2YkC0ahkF5ExIrJNREpFZHqqBkVEqSdx/51dRJoD2A5gNIAyAB8C+J6qbq7vNXl5eZqfn2+0ffXVVw2+96BBg5y2jRs3Nvg6ceXl5Rlxjx49nD4HDhxw2o4ePZq2MdkuuOACI27evLnTxx5P//79nT4lJSVO28mTJxs5Oj8Rcdr69u1rxDt27HD6xPkebtbM/ZzztdXW1ia9lj3uKOPxfT1OnTqV9HVRqKr7RgJo0YhrXgGgVFV3AYCILAIwDkC9yZ6fn4/BgwcbbatXr056I/uL8NZbbzl9unXrFmHIqXHRRRcZ8VNPPeX0+cMf/uC0rVixIl1DcowdO9aI27Vr5/RZtWqVEa9Zs8bp069fP6dt165djRydX0FBgdM2a9YsI77xxhudPsePH2/wvQoLCyPd/9ChQ0mvZf/PP8r/DH1fj6qqqqSva4zG/BjfDcDes+KyRBsR5aDGfLJHIiKTAUwGgJYtW6b7dkRUj8Z8su8D0P2suDjRZlDV2ao6VFWH2j/uEFHmNGaCrgXqJuhGoS7JPwTwL6q6qb7X5Ofna1FRkdHWp08fI162bFms8WRS586djdg3GTdy5EinbeXKlUY8f/58p89DDz1kxKdPn3b6fPnll0YcZeLPfoYH3DkE33U6derktF1zzTVGfPvttzt9pk6dasQVFRVJxxiFbzyVlZVG7Jtos99H3wdPTU1NrDEVFxcb8eHDh50+vrZ0SfkEnarWisi/AngbQHMAfzpXohNRdjXqmV1V3wTwZorGQkRpxBV0RIGI/cwe62Yizs3sGXrfeLp27WrEu3fvTnov37+Z2v+2marnyCFDhjhtvkU+vudv24wZM4z4+eefd/rs3bvXabPZz632ewgAZWVlSa/Tvn17p+3hhx824gULFjh9NmzYkPTa2WSvlQBS9/3gY7+PrVq1cvrs378/Jfeq75mdn+xEgWCyEwWCyU4UCCY7USCyPkFnT5SkapLk8ccfd9p+/etfG3GcDRTpZk/2ZXKiyzfR2LFjR6etQ4cORvzGG2/Eup+982vixIlOn4ULFzb4ut27d3fa7E0mbdu2dfrYOwUB4NixY0acrk1AqcQJOqLAMdmJAsFkJwpE1p/Z4/jGN77htO3cuTMVl8661q1bG3GcSj5RTZgwwYg/+OADp8+wYcOcNnvBkG+x0Geffdbg8bz88stO21//+lcj9i0yyqQrrrjCaVu7dm0WRlI/PrMTBY7JThQIJjtRIJjsRIHI6ARdr1699NFHHzXa7rrrrgZfZ+jQoU6bvfgkSglgSs73XtsTUr4+69evT3rtMWPGGLGveoxdOnrr1q1Jr+sriR3ldecLTtARBY7JThQIJjtRIDL6zJ6Xl6f2xoqDBw9m7P65ZuDAgU7bZZddZsS+Wvtz58414igVcNLJrrYLuHMmvtNO7I0wt9xyi9Pn1VdfbeTowsNndqLAMdmJAsFkJwoEk50oEBmdoGvWrJlzPrtvQYatTZs2RuybtNq2bds5Y8A9aurEiRNOnyhlmjPJV3XlyiuvNOK//OUvTp9MHjeUKlGObbK/hgBQWlqatjHlmuHDhxtxSUmJEX/99dc4deoUJ+iIQsZkJwoEk50oEE2yUo19RC/gVq/xLeKwn/VPnjwZ6/72YpBTp05Fep29QSOkzRmUGnZl31GjRhnxsmXLUFVVxWd2opAx2YkCwWQnCgSTnSgQTXKC7uKLL3baCgsLjdh3TI9dCeWZZ55x+syfP9+IU1kmOO7EXrr06NHDiO2yzQBwySWXOG32YiTf1+PLL7804ttvv93p89prrxlxtnfvRdG7d2+nLdeOhOKuN6LAMdmJApE02UXkTyJSISIlZ7V1EpF3RGRH4nf3qE8iyilJn9lF5B8AHAXwn6o6MNH2FIBKVZ0hItMBdFTVXyS9WYqe2dMpk8cv5Rr7+GzAnWcAgMGDBxuxbyNOFE3xve7Zs6fTFueoqz//+c9Omz2vcfXVVzt91qxZk/TasZ/ZVXUlAHvJ2jgA8xJ/ngdgfNIREFFWxX1m76Kq5Yk/7wfQJUXjIaI0adHYC6iqnuvHcxGZDGByY+9DRI0T95P9gIgUAUDi94r6OqrqbFUdqqrJq1QQUdpEWlQjIj0BLD1rgu5pAF+cNUHXSVWnRbhO2iboBgwYYMSbN2+OdR27dLGv4s1TTz1lxPYCEsA/kdO1a1cjXr16dYwRAsXFxUZcVlaW9DW+xTFHjhwxYt9uwmzr16+fEe/bt8/pc/To0QZf1y5pDgAtWrg/6Nq7KX197Pdt3bp1Tp9vfvObRvz55587fexFRfaiJwDYs2ePERcUFBjxiRMncPr06XgTdCKyEMD/AugnImUicg+AGQBGi8gOANcnYiLKYUmf2VX1e/X81ah62okoB3EFHVEgmuRGmFRZtWqV03b33XcbsV3NEwAWLFiQtjHF4Tvq2K7I66vcM22aOc3y6aefOn1mzHCf0L773e8a8euvv+70GT9+vBFn8pivmpoap23YsGFG/OGHHzp9Zs6c6bQ98MADSe9nVwCeOHGi02fWrFlGnM4FRNwIQxQ4JjtRIJjsRIFgshMFIqgJOnvhzaFDh5w+FRXmYkDfYpTdu3endmCN9OCDDzptL7zwghFXV1cnvY49qQcAI0aMcNp27NhhxL4jqvbv32/EviOaxo0bZ8RLly51+kSp5mMvdPGde29/zXwTlqk6WspeCOTz7rvvOm0TJkww4ihVkuydipWVlaipqeEEHVHImOxEgWCyEwWCyU4UiPN2gs63O8nedWaXVwLcnUa/+93vnD6tWrVq5Ojq2DuWAOD48eMpuXY6zZkzx4h/+tOfOn3irBCzy1TFvU4UvrPg27Zt67T5djTGYZf32rhxo9PnuuuuM+IDBw7EuhdX0BEFjslOFAgmO1Egzptndvvc6muuucbpY1ev2blzp9PHfpbr1auX0+dvf/ubEfvmB6Kc/X7VVVc5bXGr12STb9edb+dZLrEXWAFASUmJ03bHHXcYse/rk6pFVnEqEPnwmZ0ocEx2okAw2YkCwWQnCkSjD4nIhk6dOjlthw8fNmJfqaSf/OQnRuyboLPL+fr62GWifaWMfTvqbFOnTnXa7MkeX+nkXJOqyThfCbAoZ5vFsX37dqfNt8Nu7969RpzOHY9RJuTshVgNWYTFT3aiQDDZiQLBZCcKRJN4ZrefU/r37+/0sRc7+BZ62M/69vM5ADzxxBNGvGXLFqfPsmXLjDjK87nPz3/+c6fN3jDRVDfLRPma2Rs9fBtR7PcjSuUaH3ux1MUXX+z06d27t9N2+eWXG3G7du2cPiLmGpYom2cKCwudtkceecSIn3nmGafPY489ZsS//OUvjfhcFYn4yU4UCCY7USCY7ESBYLITBaJJ7nqzJ20Ad+Jm0KBBTh+7OohvcY591nZTKCWdi+zJJru6DeCWrt62bVtK7u37/igqKjJie+ciAHznO99x2uyJX1+57Tjnw0fh201ZW1ub9HXc9UYUOCY7USCY7ESBaJLP7OezCy+80Ih9z4P2GfL2UU8+9jMr4M5ZbNq0KcoQHb4FIva4X3rpJafPlClTjPjZZ591+jz33HNG/PXXXzt97AovH3zwgdPHXuRz7NixpNcB3PfIVxU21/CZnShwTHaiQDDZiQKRNNlFpLuILBeRzSKySUSmJNo7icg7IrIj8XvH9A+XiOJKOkEnIkUAilR1vYi0BfARgPEA7gJQqaozRGQ6gI6q+osk1+IEXQ67/vrrnTZf6WTf5FYynTt3dtoOHjyY9HVRKrPYk4++8aXqGKemIPYEnaqWq+r6xJ+PANgCoBuAcQDmJbrNQ93/AIgoRzVoP7uI9ARwGYA1ALqoannir/YD6FLPayYDmNyIMRJRCkSeoBORNgBeAzBVVY3qjlr3LOD9EV1VZ6vqUFUd2qiRElGjRPpkF5E81CX6AlU9U7b1gIgUqWp54rm+Il2DjKNly5ZOmz0/ccMNNzh9lixZkrYxpcqLL75oxL///e+dPnEWyLz77ruxx2S79NJLjXjXrl2xrhOlKk/fvn2N+P333491r/NdlNl4ATAHwBZVPbtOzhIAkxJ/ngTgv1I/PCJKlSif7N8GcCeAT0VkQ6LtIQAzALwsIvcA2A1gYlpGSEQpkTTZVfU9AN6pfACjUjscIkoXrqAjCgR3vaWJrwT0pEmTjNieaIvKLmdsH30VV5QKQHG1b9/eaYuy0OXaa6814qqqKqfPhg0bYo/rfMRdb0SBY7ITBYLJThQIPrMnEeX4pR49ejh9ysvLnTa7Co2vT1Ngb2oZOtRdHGlXiv3ss8+cPr7jt2ytW7c2Yl8VGrsqrF0h2Mc+1gnwb6DZunVr0mulypAhQ4zYNxcxffp0I54xY4bTh8/sRIFjshMFgslOFAgmO1EgOEEXCPt8ciDaBFlIokzGplOrVq2M2Fc2O9lrjh8/jtOnT3OCjihkTHaiQDDZiQLBZCcKRIMKTmaCvcsJAJYvX56FkeSukydPGrFvNZi9M85XEppMmZyMs8+QA6Kt/LM1ZFciP9mJAsFkJwoEk50oEBl9Zs/Pz3fOwN65c6cRnzhxIpNDSom8vDynraamJm3385XJjqNDhw5GXF1dnZLrUnJxns997Pmbc+EnO1EgmOxEgWCyEwWCyU4UCO56y6LbbrvNaXv11Vczdn97YjGdk4pk6tixo9M2fvx4I/7Wt77l9JkzZ44R2+W/AJalIgoek50oEEx2okBk9Jm9oKBAu3fvbrSVlpZm7P72EURRjh+68847nbbFixcbcUMWNlBu6Nevn9Nml/oGUnfW+8SJ5iHHl1xyidNn9uzZRnzrrbc6fQoLC43YLnW9du1aHD58mM/sRCFjshMFgslOFAgmO1EgMrrr7dSpU0nPErdL4wLRSura7HOzgHjneM+fP99p842RmhZ7t2Vj2Gf9Pfzww06fyZMnG7GvUo1ddWbHjh1On7179xrxvn37jLi2trbecfKTnSgQTHaiQCRNdhEpEJG1IvKJiGwSkccS7b1EZI2IlIrIYhFJTUUFIkqLpItqREQAtFbVoyKSB+A9AFMA/AzA66q6SERmAfhEVV9Icq2kK3jsqqgAkj7nx2VXzQGAsrKytNzLxz57HAC++uqrjN0/rhYtzKme0aNHO33eeuutTA0nkrlz5xrxzJkznT7Dhg1z2l566SUjvummm5w+9957rxGPGjXK6dO/f38jrqiocPocPXrUiO0FaID7zH7RRRcZcWVlJWpqauItqtE6Z0aRl/ilAK4DcGaL1jwA45Ndi4iyJ9Izu4g0F5ENACoAvANgJ4BqVT0z9VcGoFtaRkhEKREp2VX1lKoOAVAM4AoA/c/9iv8nIpNFZJ2IrIs3RCJKhQbNxqtqNYDlAK4E0EFEzjy8FQPYV89rZqvqUFUd2piBElHjRJmg6wygRlWrRaQVgP8G8CSASQBeO2uCbqOq/jHJtZyb2RU7ioqKnNdt3rzZiJs3b+70acgxOOebgQMHOm0lJSUZu3+mS2lnkv3f5ivjPWbMGCNesmSJ08fOM9/iF/vaUXZTdu7c2YirqqrqnaCLsoKuCMA8EWmOup8EXlbVpSKyGcAiEXkCwMcA5pzrIkSUXUmTXVU3ArjM074Ldc/vRNQEcAUdUSByrrqsb4OAvbigbp2PKc4ml6agTZs2Tpu9+MLHrlS6Zs0ap49dqefYsWNOn6a68CdkrC5LFDgmO1EgmOxEgWCyEwUi5yboKLf4FjmVl5cnfZ1dzSdKtaEok5H2jjvg3NVZmrI+ffo4bdu3bzdie4fdunXrWEqaKHRMdqJAMNmJApHR6rK5xt5EAAAHDx7MwkhSz94sFGWjkK9yj2+RU5Rn9jgVgX2VYgYNGmTEvgoz5wu7mk3fvn2dPva8hm8hVH34yU4UCCY7USCY7ESBYLITBSLnFtX4zkNfuHChEUdZRDF8+HCnza7eYu8MA4AFCxYkvXa23X333UbsK79tv0fPP/980uvaVYMA4IsvvnDaRowYYcTbtm1z+lRXVxvxfffd5/SxzyP3TSJ27drViO1SynH5Jh4rKytjXcue6M3kJO/3v/99I37zzTfxxRdfcFENUciY7ESBYLITBSKjz+x5eXl64YUXGm32MTinT5/O2HgyzT6qx3dE0pEjR4z4lltucfpUVVUZsb0YAwC2bt0aZ4gO33zAD3/4QyN+5ZVXnD4dOnQwYl8lId+1bfbRX7/97W+dPg8++GDS64SElWqIAsdkJwoEk50oEEx2okBkfVFNjx49jHjPnj0ZG08Uvgkye2Lt+PHjka5lL1rxLVixy2RfffXVTp/3338/0v0a6umnn3bannzySafNPtrJV1o6TvUY37Febdu2NWLfe22/Z3F23DUVURbwcIKOKHBMdqJAMNmJAsFkJwpE1ifoQuY719wuwWyvIMs0Xznjzz//3IgbUhqpoX7zm98Y8UMPPZS2e0Xhm7B94403Mnb/9u3bG7F9Xh/ACTqi4DHZiQLBZCcKREaf2Vu3bq0DBgww2j755BMjthdsUPrYZ6/7FqP4diG2bNnSiE+ePJnagWVJs2buZ5/93+9b+GNXvfEtdLFLQNvHWqUSn9mJAsdkJwpE5GQXkeYi8rGILE3EvURkjYiUishiEWmZ7BpElD0N+WSfAmDLWfGTAJ5V1T4AqgDck8qBEVFqRZqgE5FiAPMA/BrAzwDcDOAggItVtVZErgTwqKr+Y5LrODezzxcrKytLOp5LL73Uadu0aVPS11FqxJls8i3O2bVrlxGnqiTZ22+/7bTdeuutRuzbqZcqvkm8nj17GvHOnTuTXscu4Qa45+zZ+TtixAh89NFHjZqgew7ANABnvhoXAKhW1TP7GMsAdIt4LSLKgqTJLiI3AahQ1Y/i3EBEJovIOhFZF+f1RJQaUY5s/jaAsSJyI4ACAO0AzATQQURaJD7diwHs871YVWcDmA1wbTxRNjVoUY2IjATwgKreJCKvAHhNVReJyCwAG1X1j+d6fUFBgdrPLvbxPsuXL488nkzo3bu302Y/azYFl19+udM2duxYI3788cedPlHOdfexN/mcOHHC6WOXm4676cde1HLo0CGnz4QJE4w46uaVKOfc24uT7EVHgFv+O66CggIjHj16tBGvXLkS1dXVKV9U8wsAPxORUtQ9w89pxLWIKM2i/Bj/d6q6AsCKxJ93Abgi9UMionTgCjqiQDDZiQLRoB/jU8Ge4LDPNksVXxWYKDvq7rnHXAi4aNGilI0pm9avXx+pLVXs99q3o6xLly5GHHeCzj5X3XevQYMGGXG/fv2cPr5z5qNMUNoLdKIs2Il6f9uLL75oxPn5+UZ8rsVl/GQnCgSTnSgQTHaiQGS9uuy1115rxL4qH4WFhUb88ccfO32yWeHGXlQBAPPmzXPafvzjHxux77+V0sPePOXbcOWr1BqFPR8wfvx4p8+vfvWrBl/X3iQGuM/oEydONOK5c+eivLyclWqIQsZkJwoEk50oEEx2okBkfYLOnljbsmWL3QU/+tGPjNi362z//v2NHV7a2Ys9UlWZJdPs3WrV1dVZGUeu8pXWts+5t99DALjvvvuSXtvOl3bt2hnx8ePHcfr0aU7QEYWMyU4UCCY7USAyvhHG5tuwYps2bZoR51qlmBYt3LextrbWacu1Z/S4cwhRntGjVHg5X11//fVOW/fu3Y141apVsa5t54v9zO6rCHQGP9mJAsFkJwoEk50oEEx2okBkfVHN4MGDjfgHP/iB87r777/fiO1dcABw7NgxIx4xYoTTxy4xXFpaeo7R1u/mm282Yt8xPXPnzo11bTJle6KvY8eORpyqktBRquls2LAh1rV5PjtR4JjsRIFgshMFIuvP7Dbfs8xVV11lxPaxtQCwZ88eI05V5ZohQ4Y4bXGfpWwjR4502lasWNHg6/gqldpVVzK9UcheaORbZJQuvjmU/v37G/F7773n9LGfzwG36oyvAuzq1asbOEJg+PDhTtvAgQONePHixU4f+7gn3zFWfGYnChyTnSgQTHaiQDDZiQKR6Qm6gwB2A7gQgHuIdm5rimMGmua4Oeb4LlHVzr6/yGiy//2mIutUdWjGb9wITXHMQNMcN8ecHvwxnigQTHaiQGQr2Wdn6b6N0RTHDDTNcXPMaZCVZ3Yiyjz+GE8UiIwnu4iMEZFtIlIqItMzff8oRORPIlIhIiVntXUSkXdEZEfid3chdRaJSHcRWS4im0Vkk4hMSbTn7LhFpEBE1orIJ4kxP5Zo7yUiaxLfI4tFpGW2x2oTkeYi8rGILE3EOT/mjCa7iDQH8DyAGwAMAPA9ERmQyTFE9B8Axlht0wEsU9W+AJYl4lxSC+B+VR0AYASAexPvbS6P+wSA61R1MIAhAMaIyAgATwJ4VlX7AKgCcE/2hlivKQDOPr4o58ec6U/2KwCUquouVT0JYBGAcRkeQ1KquhJApdU8DsCZQ9fnARifyTElo6rlqro+8ecjqPtG7IYcHrfWOZoI8xK/FMB1AF5NtOfUmAFARIoB/BOAf0/EghwfM5D5ZO8GYO9ZcVmirSnooqpn9tbuB9Alm4M5FxHpCeAyAGuQ4+NO/Di8AUAFgHcA7ARQrapn9sTm4vfIcwCmAThTaP8C5P6YOUEXh9b9E0ZO/jOGiLQB8BqAqap6+Oy/y8Vxq+opVR0CoBh1P/n1P/crsktEbgJQoaofZXssDZXpE2H2ATj7aIziRFtTcEBEilS1XESKUPdJlFNEJA91ib5AVV9PNOf8uAFAVatFZDmAKwF0EJEWiU/KXPse+TaAsSJyI4ACAO0AzERujxlA5j/ZPwTQNzFz2RLAPwNYkuExxLUEwKTEnycB+K8sjsWReG6cA2CLqj5z1l/l7LhFpLOIdEj8uRWA0aiba1gO4LZEt5was6o+qKrFqtoTdd+//6OqdyCHx/x3qprRXwBuBLAddc9m/5bp+0cc40IA5QBqUPf8dQ/qnsuWAdgB4F0AnbI9TmvMV6PuR/SNADYkft2Yy+MGMAjAx4kxlwB4JNHeG8BaAKUAXgGQn+2x1jP+kQCWNpUxcwUdUSA4QUcUCCY7USCY7ESBYLITBYLJThQIJjtRIJjsRIFgshMF4v8A73K2sd4BBFgAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPsAAAD6CAYAAABnLjEDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAASgklEQVR4nO3da4xVVZoG4Pe1KO4Ul+JiSRWDoqEDMqIiaaPR0RGC2gGNnRbTP5iEBGNmEjs9Ru0xGe3oD5kfdrfJXCSjDqPY2NOQYDTK2DQGNSNyK6BAaS7SEVNQXCyugly++XE2ndprLTi7zq3OqfU+yUnVt2qdsz/L+thnrbP22jQziEjvd0VPJyAilaFiF4mEil0kEip2kUio2EUioWIXiURRxU5yFskdJHeRfLpUSYlI6bHQz9lJ1gH4E4AZAPYBWAfgETPbfpnn6EN9kTIzM4baizmzTwewy8z2mNn3AJYCmFPE64lIGRVT7GMBfN0l3pe0iUgV6lPuA5BcAGBBuY8jIpdXTLF/A6ClS9yctKWY2SIAiwCN2UV6UjFv49cBuI7k1ST7ApgL4J3SpCUipVbwmd3MzpH8BwArAdQBeM3MtpUsMxEpqYI/eivoYHobL1J25fjoTURqiIpdJBIqdpFIqNhFIqFiF4mEil0kEip2kUio2EUioWIXiYSKXSQSKnaRSKjYRSKhYheJhIpdJBIqdpFIqNhFIqFiF4mEil0kEip2kUio2EUioWIXiYSKXSQSKnaRSKjYRSKhYheJhIpdJBIqdpFIqNhFIlHM/dmlAl5//XWvrb29PRU/+uijXp+XX345Fffr18/r09HRkYrXrFnj9Tl27JjX9t133+V9bfd5R44c8fo0NDSk4tOnT3t9mpubU/Htt9/u9Vm6dGne1xGd2UWioWIXiYSKXSQSKnaRSNDMKncwsnIHK6PGxsZUfPjw4bIda9iwYV7bkCFDUvGJEye8PmfOnEnFgwcP9vq4bXv27CkgQ6B///5e2xVXpM8jp06dKui1n3zyyVQcmuh74403UrH73x4bM2OoXWd2kUio2EUikbfYSb5GsoNkW5e2ESQ/JLkz+Tq8vGmKSLHyjtlJ3gHgBID/NrPrk7Z/AXDEzF4k+TSA4Wb2VN6DVdmY3R37AsDo0aNT8e7duyuVTpA79gWACxcudPt1+vbt67V9//33BeVUiNDCmyxj64kTJ6Zid0EPAOzfvz8VZ/nvamlp8dq+/vrrvM+rBQWP2c1sDQB3VmQOgMXJ94sBPFBMciJSfoWO2ceY2cU1m/sBjClRPiJSJkWvjTczu9zbc5ILACwo9jgiUpxCz+wHSDYBQPK141IdzWyRmU0zs2kFHktESqDQM/s7AOYBeDH5uqJkGVVQ6Aqq1atXp+JBgwZ5fU6ePFmS4w8dOtRrO3r0aCoeNWqU1+fAgQOpmPTnY+6///5U/P7773t93Mm/0MRf6LVdoUled0Kw0IUuO3bsKOh5+fSWybjuyPLR228B/B+AiST3kZyPXJHPILkTwD1JLCJVLO+Z3cweucSP/rbEuYhIGWkFnUgkdCFMlRs4cKDX5l5UElow415A4+5KU0rTp0/32ubOnZuK3Z1zAGDv3r3lSsnjXqzTm3ez0YUwIpFTsYtEQsUuEgkVu0gkNEGXR2hRSTl/Zy+99FIqXrVqldens7MzFW/atMnrU+jOMIUYMGCA1zZhwoRU3NbW5vVx1dXVeW0zZ85MxVu3bvX6uFfCnT171uvj7uYTWiwV2vGnkvVRKpqgE4mcil0kEip2kUio2EUioQm6PK666iqvzb3q7Pz582U7/uTJk722bdu2peJSbV0VusLu4MGDeZ8X2t7r+PHj3T5+iLudVejqOXeCMNTHXUFXyQnMStMEnUjkVOwikVCxi0RCY3aHO24NLapxr5gK3cO8FmXdlWfGjBmpuLW11evjjuMLvbWUq6e3xK4FGrOLRE7FLhIJFbtIJFTsIpHQBJ2UhbtNtrtFdqVl2Ta7t9AEnUjkVOwikVCxi0RCY3apiEIvsnF3r5k1a5bX54MPPkjFoQuTGhoaUnHWhVDu8UO76VTboh6N2UUip2IXiYSKXSQSKnaRSGiCLg93hxOgd98nLItrrrkmFZfqirbGxkavzd3euampKW+fQ4cOFXT80aNHe21TpkxJxevWrfP6VNtVj5qgE4mcil0kEip2kUhozF7lhg8f7rW5cwbu7Y8qLbS7rTvXUV9f7/Up1cUxV155ZSoOjdnPnTuX93VCu+Bs2LAhFYfuRd/Tv3+XxuwikVOxi0RCxS4SibzFTrKF5GqS20luI/l40j6C5IckdyZf/cGliFSNvBN0JJsANJnZRpJDAGwA8ACAvwNwxMxeJPk0gOFm9lSe19IEndQUd/JvzJgxXp/NmzdXKp1MCp6gM7N2M9uYfH8cwBcAxgKYA2Bx0m0xcv8AiEiV6tOdziTHA7gRwFoAY8ysPfnRfgD+P3m55ywAsKCIHEWkBDJP0JEcDGAZgJ+ZWWoxsOXGAsG36Ga2yMymmdm0ojIVkaJkOrOTrEeu0JeY2fKk+QDJJjNrT8b1HeVKUqTUsl7gtH///svGtSTLbDwBvArgCzN7qcuP3gEwL/l+HoAVpU9PREoly2z87QA+BrAVwMXNtv8JuXH77wCMA/BnAD8xsyN5Xkuz8VIVevOly5eajdfaeIlSjMWuFXQikejWR28SNm7cOK8tdAWVu6NL6Iq2w4cPlyQndzFIaBeWLVu2lORY1eauu+7K22fbtm1eW285s1+KzuwikVCxi0RCxS4SCc3G5xG63c/kyZNT8c6dO70+/fr189o6OztLllc+AwcOTMWDBg3y+rg7xRR6G6Pnn38+b59ly5Z5ba2trQUdL5+VK1d6bc8++2wq/vbbb70+7e3tXlu17RybhWbjRSKnYheJhIpdJBIqdpFIaIKuBEILaCp5z253whAILxopl9B/f0tLSyqePXu218edSDtw4IDXx11kFLrPu9vnwoULXh83x2q7p3opaYJOJHIqdpFIqNhFIqFiF4mEJugiEbofW2giqxB33nmn17Zx48ZUPGHCBK9PW1tbKg7dj83dujk0iSdpmqATiZyKXSQSKnaRSGjMXuWGDBnitR0/fjwVDxgwwOuT5Z7h7m42tbxNcj4jR45MxadOnfL6hNpqkcbsIpFTsYtEQsUuEgkVu0gkot5KOnSjgLNnz6bi8+fPVyqdoHvvvddrW7p0aSpevny51+eJJ55IxaHJt4ceeigVr1jh38Hr4MGDXtuZM2fCyRapVAt/QvdQdxfs9JbJuO7QmV0kEip2kUio2EUioUU1Vc5d+AL42yAPHTrU63PzzTen4o8++sjr4z7vnnvu8fq8+eabWdKsObqxo4j0Wip2kUio2EUioWIXiYQm6KqMe0+20NVrhSw0CW337C40Cd0PLrTlsnslXuie8u4CmT59/PVbZHoeKbRYx/37DN1nvqOjw2srxB133OG1HTlyJBVff/31Xh93kVNP0wSdSORU7CKRyFvsJPuT/JzkZpLbSP4yab+a5FqSu0i+TdJ/nygiVSPvmJ25gdUgMztBsh7AJwAeB/BzAMvNbCnJ/wCw2cz+Pc9raczeRejCjxdeeCEVv/XWW14fd4waGrMuWbIkFTc3N3t9nnvuuVQcuuhl8ODBXtv69etTcWhX2CxmzpyZim+66Savz+bNm1Oxu2stUNiOs/369fPaynWBT6UVPGa3nBNJWJ88DMDdAH6ftC8G8EDxaYpIuWQas5OsI9kKoAPAhwB2A+g0s4v/pO8DMLYsGYpISWQqdjM7b2ZTATQDmA7gB1kPQHIByfUk1+fvLSLl0q3ZeDPrBLAawK0AhpG8+OFpM4BvLvGcRWY2zcymFZOoiBQnywTdKABnzayT5AAA/wtgIYB5AJZ1maDbYmb/lue1NEFXIfPmzUvFR48e9fq4C1ReeeUVr8/DDz/stQ0fPjwVf/nll14fd9Lsq6++8vrMmTMnFW/fvt3r416d5k4Oiu9SE3RZtqVqArCYZB1y7wR+Z2bvktwOYCnJFwBsAvBqybIVkZLLW+xmtgXAjYH2PciN30WkBmgFnUgkdCFMmbgXeQBAfX19Kg5dZJKFu3vNwIEDvT4jRoxIxZMmTfL6tLe3p+LQuPrEiRNeW5bbRLmLVkI5uhfehOYVTp48mYpDF9QU+nt0hRYe7du3LxXXwmIcXQgjEjkVu0gkVOwikVCxi0RCE3QFqKurS8VZbxFV6PNc7mRX6HXcSaMs93kPKdUtmcaPH5/3tQ8dOuT1OXbsWLePFTtN0IlETsUuEgkVu0gkNGavQe7CjtCiDvfWTqEFK+7inCyLZUpp1KhRqTi0U450n8bsIpFTsYtEQsUuEgkVu0gksmxeIY6pU6em4tbWVq9P6Kq3Uk2GtrS0pOIpU6Z4fbZs2ZKKH3vsMa/P2rVrU3FoS+pCFtBkVYsTcqGr7grdSrvSdGYXiYSKXSQSKnaRSGhRTQ1yby28Zs2avM+54YYbvLZbbrklFX/88cdeH3fhCwB88skneY/XW0yblt4BPXSxzt69eyuUTTZaVCMSORW7SCRU7CKRULGLREITdI4HH3wwFYd2c3EXVrS1teXtAxQ2kRPaYcZ97dDCl9BVbq4JEyak4t27d3czu9IaN26c19bU1JSKQ79rd7vpEPdWV/Pnz/f6LFy40GsbOXJkKg5dYZjld11JmqATiZyKXSQSKnaRSKjYRSKhCboSCK0yC13R5W4DNXHiRK+PO9m0detWr0+13VusnNz7wWWZjAtxt7J27x8PAI2NjV7b2LFjU3Fo4vXTTz8tKKdy0QSdSORU7CKRULGLREJj9jzcK8yAbFeZubd6AoCGhoZUHBp/uvdwDzl9+nQqDi2qqeT/10I988wzqTg09+HupvPee+95fbLcIirL2P/aa6/12nbs2JGKP/vsM6/Pbbfdlvf4laQxu0jkVOwikchc7CTrSG4i+W4SX01yLcldJN8m2bd8aYpIsbpzZn8cwBdd4oUAfmVm1wL4FoB/ZYGIVA8zy/sA0AxgFYC7AbwLgAAOAeiT/PxWACszvI7poYce5X1cqv6yntl/DeBJABenfRsBdJrZxQ2z9wEYG3ieiFSJvMVO8kcAOsxsQyEHILmA5HqS6wt5voiURpY7wtwGYDbJ+wD0B9AA4DcAhpHsk5zdmwF8E3qymS0CsAiozc/ZRXqLvGd2M/uFmTWb2XgAcwH80cx+CmA1gB8n3eYBWFG2LEWkaMV8zv4UgJ+T3IXcGP7V0qQkIuWg5bIivYyWy4pETsUuEgkVu0gkVOwikVCxi0RCxS4SCRW7SCRU7CKRULGLRELFLhIJFbtIJFTsIpFQsYtEQsUuEgkVu0gkVOwikVCxi0RCxS4SCRW7SCRU7CKRULGLRELFLhIJFbtIJFTsIpFQsYtEQsUuEgkVu0gkVOwikVCxi0SiT4WPdwjAnwGMTL6vJbWYM1CbeSvnwv3VpX5Q0Vs2/+Wg5Hozm1bxAxehFnMGajNv5VweehsvEgkVu0gkeqrYF/XQcYtRizkDtZm3ci6DHhmzi0jl6W28SCQqXuwkZ5HcQXIXyacrffwsSL5GsoNkW5e2ESQ/JLkz+Tq8J3N0kWwhuZrkdpLbSD6etFdt3iT7k/yc5OYk518m7VeTXJv8jbxNsm9P5+oiWUdyE8l3k7jqc65osZOsA/CvAO4FMAnAIyQnVTKHjP4LwCyn7WkAq8zsOgCrkrianAPwj2Y2CcAPAfx98rut5rzPALjbzG4AMBXALJI/BLAQwK/M7FoA3wKY33MpXtLjAL7oEld9zpU+s08HsMvM9pjZ9wCWAphT4RzyMrM1AI44zXMALE6+XwzggUrmlI+ZtZvZxuT748j9IY5FFedtOSeSsD55GIC7Afw+aa+qnAGAZDOA+wH8ZxITVZ4zUPliHwvg6y7xvqStFowxs/bk+/0AxvRkMpdDcjyAGwGsRZXnnbwdbgXQAeBDALsBdJrZuaRLNf6N/BrAkwAuJHEjqj9nTdAVwnIfYVTlxxgkBwNYBuBnZnas68+qMW8zO29mUwE0I/fO7wc9m9HlkfwRgA4z29DTuXRXpdfGfwOgpUvcnLTVggMkm8ysnWQTcmeiqkKyHrlCX2Jmy5Pmqs8bAMysk+RqALcCGEayT3KmrLa/kdsAzCZ5H4D+ABoA/AbVnTOAyp/Z1wG4Lpm57AtgLoB3KpxDod4BMC/5fh6AFT2YiycZN74K4Asze6nLj6o2b5KjSA5Lvh8AYAZycw2rAfw46VZVOZvZL8ys2czGI/f3+0cz+ymqOOe/MLOKPgDcB+BPyI3Nnqn08TPm+FsA7QDOIjf+mo/cuGwVgJ0A/gBgRE/n6eR8O3Jv0bcAaE0e91Vz3gD+GsCmJOc2AP+ctF8D4HMAuwD8D4B+PZ3rJfL/GwDv1krOWkEnEglN0IlEQsUuEgkVu0gkVOwikVCxi0RCxS4SCRW7SCRU7CKR+H9erbKtev4VZAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPsAAAD6CAYAAABnLjEDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAa70lEQVR4nO2de4zV1bXHv8vhIW8BESggSBAUKWKkRIsaApfiW2OMisZgtcWYa0R7jXI12kv1GnttVFJvVXJtRWscH9BKKSmKjoBW3iAXhvdEYCgwIIyAFuSx7h9zuGGvvZnzm9+c17C/n4Qw3806v9+eM2dxzl6zHqKqIISc+pxW7A0QQgoDnZ2QSKCzExIJdHZCIoHOTkgk0NkJiYRGObuIXCki60Rko4hMzNWmCCG5R9L+nl1EygCsBzAaQDWAxQDGqmplPY/xblZWVubo008/3Xvct99+m2qPhJQqzZo189asLxw6dCjVtVVVgvdMdbU6hgHYqKpVACAi5QBuAHBSZw/RoUMHR/fv39+zWbBgQepNElKKnHnmmd6a9YV169bl9J6N+RjfA8DWE3R1Zo0QUoI05p09ESIyHsD4fN+HEFI/jXH2bQB6naB7ZtYcVHUKgClA+MxOCCkMjQnQNUNdgG4U6px8MYDbVXV1PY8pKWd/4IEHvLUVK1Y4uqamxrNZu3ZtvraE5s2bO/qJJ57wbGbNmuXogwcPejarVq1y9JEjR3KwO5Ir7rvvPm+tXbt2jp4+fbpnc+zYMUdXVVV5NjkP0KnqERG5H8BsAGUAfl+foxNCikujzuyqOgvArKyGhJCiwww6QiIh9Zk91c1K7Mwewib1tGrVyrPZu3dvobYTxCZkJDmPh36vu3v37pztqdSxP9fDhw8nepyIe/zNVeyjY8eO3ppNqkn78znZmZ3v7IREAp2dkEigsxMSCXR2QiKh6AE6G5Q4evRowfaThMsuu8xbW7x4saPTVieFsAHBFi1aeDbffPNNzu5HGka3bt28tR07duTk2kOGDHF0KBi4adMmR//zn//0bBigIyRy6OyERAKdnZBIKPqZvdTp0cMv0d+2zSvu8xg1apS3dumllzr6mWee8Wx+8pOfOLq2ttazsWfE0Llt586dWfeYFluwsX///rzdqyny/fffe2v33nuvo8vLyz2bzp07O7q6ujrV/XlmJyRy6OyERAKdnZBIoLMTEgkM0OUAmwyR1O4vf/mLZ/P11187evTo0Z6NDRCGupXY7jWnneb/v26vPXv2bM/mvPPO89Y2bNjg6JYtW3o23333nbdW6tguQQDQqVMnRycJfPbp08db27dvn6P37NnTsM01AAboCIkcOjshkUBnJyQSCn5mt2dH2y2zKdC3b19Hh87Mv/nNb7y1iRPdcXihIhfbwSR09rVdVg4cOHDyzWa45ZZbvLVdu3Y5uqKiIut1AOCLL75w9EsvveTZ2A68S5cuTXTtYnLXXXd5azNmzHB06GeWq0KYXMEzOyGRQ2cnJBLo7IREAp2dkEgoaICuVatWaoNbtvNGLru+FBNbwQT4CTNDhw71bGyCig2iAcD69esdHUqYsRV2oQ5AlZXudO1QtVYoYea2225z9MKFCz0bey17r7Tceuut3pqtDgsFA0MjsoqJbW0N5G6PDNAREjl0dkIigc5OSCSwEKaAtG/f3tGhEUC2iGLu3LlZr2uTdQDgjDPOcPT777+f9Tp2zDOQ7BwZ+j5sh59QR1w7RitJclCIQYMGOXrLli2ezbnnnuvoppDkkxae2QmJHDo7IZFAZyckEujshERC0QN01113naNDc8QvuOACR4c6qnzyySeOLuQYqVCnmi5dunhrtn3w66+/7tk8/fTTjg7NgrcdVUKzxm1Sj03oAfzAWqglda4SPeyYL8D/GYU6xdgEnmnTpnk2NhHrueee82xeffVVR69bt+7km62HUOegAQMGOHrFihWezWeffZbqfmlggI6QyKGzExIJWZ1dRH4vIjUisuqEtU4i8pGIbMj87f+ilRBSUmQ9s4vIFQAOAHhDVQdl1v4LwB5VfVZEJgLoqKqPZr1ZyqSawYMHO3rr1q2ejU3ayFUHnNC97LlxwYIFnk2ooOfxxx939M033+zZ2ISQtN9HPkdh9+vXz9GhOItN0AklzNgCnlDM4Pzzz3d069atPRsbj9i+fbtnk09sItRXX33l2diCopEjR3o29nUUitckIfWZXVXnAbB9b28AMDXz9VQAN6baFSGkYKQ9s3dV1eP/fe4A0DVH+yGE5Ilmjb2Aqmp9H89FZDyA8Y29DyGkcaR9Z98pIt0BIPN3zckMVXWKqg5VVb9TAyGkYCRKqhGRPgBmnhCgew7A1ycE6Dqp6iMJrlPUqrdmzdwPMqHZ65s3b3Z0KGFm7Nixjn7qqac8m1Al2MCBAx0dmmtuEzLs+CHA78wSsrHJMKGElVAlmiWUDGMr6kIJOzYglaQDUf/+/b0125UnCaHnI5/jlmzQ8KyzzvJsbAvqUItw25L6yJEjqfaTOkAnIm8D+ALAABGpFpF7ADwLYLSIbADwLxlNCClhsp7ZVXXsSf5pVI73QgjJI8ygIyQSil4Ikwbb8QXwEzJCxSFNgTfeeMPRoS6kV111laMvv/xyz8aeEUNnVnvWbtOmjWdj4wOFplWrVo4OJd4UkgkTJnhrjz7q5pPNnDnTs5k3b56j//jHP3o2bdu2dXTazj0shCEkcujshEQCnZ2QSKCzExIJJRegCyWjdOjQwdGhqqJc8eCDDzr6xRdf9GyGDRvm6EWLFuVtP6HAmt1TqCuO3ZNNKAKAP/zhD43bXIbQiKrHHnvM0aHA2nvvvefoUOKNTZCxFXdAfp//JPzsZz9z9JtvvunZ2O8tSeeetDBAR0jk0NkJiQQ6OyGRQGcnJBJKLkBHXEIz3G3rqpqak1YYF4RQsOmRR9wiyBdeeMGzsZV5d9xxh2fz29/+1tF/+9vfPJuXXnrJ0TfddJNn8/bbbzs67ay3JHPVQ226bIVjkirAtDBAR0jk0NkJiQQ6OyGRwDN7AbHn79Ac8WKfv9PQt29fb62qqirr4+xZ3475AvyOLuecc45nY1swh7rA2JbUoVn0uSLUfvyBBx5w9Pz58z2b3bt3O9p2+wH8ysRQ0hXP7IREDp2dkEigsxMSCXR2QiIh6gBdkgSJfBKqRLPtg0OzzUIBKIude79w4ULPpikEA3v16uXo2tpazybUkrvUsK3U9u3bl+o6tk2X5eDBgzh27BgDdITEDJ2dkEigsxMSCY0e7NiUKeT5PESS8T4/+MEPvLWNGzc6OlQsM2DAAEeHEj1ydWa/9dZbvTWb6GLHaiXFzijv1q2bZ2NHW+Vz1FMSQufqiy++2NEVFRWprt2YVtp8ZyckEujshEQCnZ2QSKCzExIJUSfVlCIrV650dCiw9dZbbzm6vLw8r3uy2Ln2oRbINmEo7cy46dOnO3ry5MmeTbt27Ry9du3arNe1Qc5CkyShKi2seiMkcujshEQCnZ2QSOCZ3WBHIoXOiNu3b3f0zp07c3Z/25m0T58+ns2SJUtydr802OIce2YGcvec2HFPoVFTdo556DV97NixnOwnCaEkp7lz5zraxiIAf/RYqFjGfm+hcz7P7IREDp2dkEigsxMSCVmdXUR6iUiFiFSKyGoRmZBZ7yQiH4nIhszf/qxlQkjJkDVAJyLdAXRX1WUi0g7AUgA3ArgLwB5VfVZEJgLoqKqP1net1q1ba//+/Z01m2wxaNAg73Gff/65o22VE+CP81m9erVn8/DDDzs6lLByySWXOPrw4cNZ7xUb9vkPBb9yNWu81Ojdu7e3tm3bNkd37drVs3nooYccHQqs/epXv3J0ko5EIVIH6FR1u6ouy3y9H8AaAD0A3ABgasZsKur+AyCElCgNqmcXkT4ALgKwEEBXVT3+O6gdAPz/zuoeMx7AeCD8jkwIKQyJA3Qi0hbANAAPqqrzC0CtOwsEzwOqOkVVh6rq0FA+MCGkMCRKqhGR5gBmApitqs9n1tYBGKGq2zPn+k9VdUCW6+QtqcaOIDr77LM9m+XLlzv6m2++afB1gWSjjU5l7r33XkdPmzbNs7GjjApJPotM0mK711x44YWeje3uk5bUZ3YREQCvAVhz3NEzzAAwLvP1OAAfNHaThJD8keRz9XAAdwL4XxFZkVl7DMCzAN4VkXsAbAZwS152SAjJCVmdXVU/AxD8WABgVG63QwjJF8ygIyQSmmTVW9u2bb01W/lESotcjdqyM90BP4Hnvvvu82xefvnlBt+rqcKqN0Iih85OSCTQ2QmJhKKf2Vu2bOnoH/7wh97jbGeWJGf2Dh06eDZjxoxx9IwZMzwbW8QQKpax589Q0UeogCZmTjvNf1+xBTRDhgzxbGwn24kTJ3o248aNc7TtJAT4Y5NCMYS6lBKXQ4cOOTqU8m2vlSRZK5/wzE5I5NDZCYkEOjshkUBnJyQSChqga9++vQ4bNsxZs3PDQ9Vq3bt3d/Ts2bM9m1zNGk/CwIEDHV1ZWVmwe5cCtoLroosu8mxsYNN2c0mKDeCGAmRJOrrYwGsoiBfCVqeF5qPbmfHz5s1LdO18wQAdIZFDZyckEujshEQCnZ2QSChogK5FixZqZ5nZzDcb/AGAr7/+2tGhIE2aCipCTsS2EQf812NorrsNMhcbBugIiRw6OyGRQGcnJBKKXvVGSofBgwd7aytXrizCTkqHjh3dEYahpJpSixfxzE5I5NDZCYkEOjshkUBnJyQSSm7SYr9+/bw1m7RgWwUVG1uVBwA//vGPvbVOnTo5+o033vBs8vW9XXPNNd7anDlzHN2iRYu83LsUGDRokKNDiTChdlJ79+7N257ScMEFFzh69erViR/Ld3ZCIoHOTkgk0NkJiYRTNqkmVNRgZ3SvW7fOs9m/f3/Wa9tCHHuOAsKtitesWePofLYctgkyoeQY271l586dqe7Vu3dvb8221w49r7n6/m2nmOHDh3s2FRUVjt6zZ09O7l2KMKmGkMihsxMSCXR2QiKBzk5IJDSJAJ0NiIXmuO3evdvRoY43tmKpWTM/p8gG8fJJkvlnabn++usdvWXLFs9mx44d9eqkhJKKbLCvf//+no3tUlRdXe3Z2Ll+obl69ueaZBbgqUK7du0c/e233+Lo0aMM0BESM3R2QiIhq7OLyOkiskhEvhSR1SIyKbN+jogsFJGNIvKOiJy6idWEnAJkPbNL3dDqNqp6QESaA/gMwAQAvwAwXVXLReQVAF+q6stZrsVONQWiMQUT2Zg/f76jL7vsMs9m4cKFjg7FUF555RVHv/xyvS8fkpDUSTVax/HoRvPMHwUwEsD7mfWpAG5s/DYJIfki0ZldRMpEZAWAGgAfAdgEoFZVj4euqwH0yMsOCSE5IZGzq+pRVR0CoCeAYQDOS3oDERkvIktEZEm6LRJCckGDovGqWgugAsClAM4QkeO/qO4JIDiTV1WnqOpQVR3amI0SQhpHkgBdFwCHVbVWRFoB+BDArwGMAzDthADdSlX9XZZrFSxAZ+d6h9b27dvn2ZSVlWW9diixI19cfPHF3prtZrNq1SrPxlbdhZ6PtFVn9tqffvqpZzN+/HhHh56zf/zjH44OjfWqqalp8P6KnSxVbE4WoEvSlqo7gKkiUoa6TwLvqupMEakEUC4iTwNYDuC1nO2WEJJzsjq7qq4EcFFgvQp153dCSBOAGXSERELJdZdNy5133unoDz/80LOxBRtdunTxbGwySOhcX0hC3U1Hjx7t6F27dnk2dsy17a4K+B1wFy1a5NkcPnzYW1u6dKmj7777bs+msrLSWysUv/zlL721d955x9GhLkWh79V2vfn8888bubvkPPnkk97asmXLHD1v3jxH11fww3d2QiKBzk5IJNDZCYkEOjshkdAkOtXkCpts0VQTLewop1Cb6BEjRjh68eLFns1ZZ53l6FA3m1DQynahKbXnMUmXoqTYNtWhANh3333n6PPO87PJ7Zz3H/3oR57NlClT6r0u4Cc0hWbDs5U0IZFDZyckEujshERCVGf2UwU7bmnbNr/gMMk52sYwQmfN0BnVdsWtqqrKeq8khApYrr32Wke3adPGs7EJMkuW5K6aetSoUY5esWKFZ9OrVy9HT58+3bOxCVzXXXedZ7N8+XJHpy1U4pmdkMihsxMSCXR2QiKBzk5IJJwyVW9psKNzAL8zTCjRJFcBqRBjxoxxdGgkk+3e0rp1a8/GBo1CwR47JinU8SZXCSqhLjQ2YadHD79n6eWXX+7oUHKQHf2VhD//+c/emg3GAf5z1KlTJ8/G/oxGjhzp2fz0pz919ODBgz2bUMefXMJ3dkIigc5OSCTQ2QmJBCbVGDp37uxo2/GlFAmNfr799tsdXV5e7tkMHep29w5195k8ebK39swzzzg6dIa3cYRQUYcdG7Vp0ybPpra2Nuu9kmDvFRoPbbvSAH7Rz5w5c1Ldv5AwqYaQyKGzExIJdHZCIoHOTkgkRJ1UEyJJQM7OPg+1ck4ztigtoaSagQMHOjrUkvr555939BNPPOHZ/PWvf/XWkgTJbBecUMLMmWee6ej169d7Nt9//33WeyXh73//u6OPHTvm2Xz11Vfemk2qmTRpkmfz3HPPObpFixaezZ49e5Js06Ffv37e2tixYx1tu9vU9/rlOzshkUBnJyQS6OyERAKdnZBIYAZdFkLBL/uchQIyaVsK5QobIOvbt69nYyvI8jl3/sorr/TWbDZaqbWkLkVsNeOQIUMcPW/ePNTW1jKDjpCYobMTEgl0dkIigWd2UhCSdKoJtZI+Vc/xV1xxhbdmZ62HsGOk9u/f7+gjR46w6o2Q2KGzExIJiZ1dRMpEZLmIzMzoc0RkoYhsFJF3RMT//RMhpGRoyDv7BABrTtC/BvCCqvYDsBfAPbncGCEktyQK0IlITwBTAfwngF8AuA7ALgDdVPWIiFwK4D9UdUw9lym5AF2SoFGIYreu+vnPf+5oG6QBgDVr1jj6yy+/zOue0mBbXk2YMKFIO0mO/dkDwFNPPeVo28oKAN58801HhxKxbNvs0M/1T3/6k6MPHTrk2TQ2QPcigEcAHK8L7AygVlWPh0qrAfg1jISQkiGrs4vItQBqVHVpmhuIyHgRWSIiuRutSQhpMEmaVwwHcL2IXA3gdADtAUwGcIaINMu8u/cE4M8NBqCqUwBMAUrvYzwhMdGgpBoRGQHgYVW9VkTeAzBNVctF5BUAK1X1d1keT2dvILbQAfDHBNnuJQCwYMECR4fiE7YTSmj2+MGDB721UJeXU5Vu3bo5evbs2Z7N+eef7+hQQZEdLRXq3LN161ZHhzrnhMaBWfKRVPMogF+IyEbUneFfa8S1CCF5pkE96FT1UwCfZr6uAjAs91sihOQDZtAREgl0dkIioUlUvdkZ4WnnfaUh1KmmS5cuju7evbtnYwNkQO6+D9uCOZRYEUrIsLRs2TLrdWIiNIvezmPfts3/pZMNdIaee/uaWbVqVZotJoJVb4REDp2dkEigsxMSCU3izF5M2rVr563ZTq2hTrKhc5vtMmJHCwF+YktZWZlnk88usMXkoYce8tZCiT6WioqKBt/rlltu8dbefffdBl8H8F8jSeIlIdq3b+/offv2pboOz+yERA6dnZBIoLMTEgl0dkIioaDz2cvKyrygVLHHJGUjlDBj9xyqTgphAzm7d+/2bGxA7lQNxoUItZK2c+bTzDkH/ISZUDDOjlYCgJEjRzo69LOeO3du1vt36NDB0Xfffbdn8+qrr2a9TmPgOzshkUBnJyQS6OyEREJBz+xHjx7FgQMHCnnLRrN+/fqcXWvLli2ODiXV9OzZ09GhDjOWzZs3e2tJuuQWG9tNdv78+Z7NypUrHW2Ld0Kcdpr/HpYkeSzUOdZ2gV22bJlnY5OlRowY4dnY131lZaVnE4pZ5BK+sxMSCXR2QiKBzk5IJNDZCYkEVr1lIRQQsmtpq5NCTJ8+3dGhJJLy8nJHz5kzJ9W9bIAw1JJ65syZ3lrv3r0d/cEHH3g2w4cPd/TGjRvTbDEVd9xxh7f28ccfOzr0cxXxi8Wqq6sdnWRevK2KBICqqipH9+nTx7NJmpyVDVa9ERI5dHZCIoHOTkgknLJn9kmTJnlrdtxtqAuK7dwaKliYNWuWo3PZKdR2M73ppps8m7Vr1zo61Mk2ydnSJuzYMUYAMG7cOG/t/vvvd3QoZmHPrWm7tyQhSZdcm7ASSlZq06aNt2aLlUIdaM8++2xHb9iwwbOxBU6hBJ4ko52SwDM7IZFDZyckEujshEQCnZ2QSCh0gG4XgM0AzgTgt2kpbZrinoGmuW/uOT29VbVL6B8K6uz/f1ORJao6tOA3bgRNcc9A09w395wf+DGekEigsxMSCcVy9ilFum9jaIp7BprmvrnnPFCUMzshpPDwYzwhkVBwZxeRK0VknYhsFJGJhb5/EkTk9yJSIyKrTljrJCIficiGzN8d67tGoRGRXiJSISKVIrJaRCZk1kt23yJyuogsEpEvM3uelFk/R0QWZl4j74hIi2zXKjQiUiYiy0VkZkaX/J4L6uwiUgbgvwFcBWAggLEiMrD+RxWF1wFcadYmAvhYVc8F8HFGlxJHAPybqg4EcAmAf808t6W870MARqrqhQCGALhSRC4B8GsAL6hqPwB7AdxTvC2elAkA1pygS37PhX5nHwZgo6pWqer3AMoB3FDgPWRFVecBsC1ibgAwNfP1VAA3FnJP2VDV7aq6LPP1ftS9EHughPetdRzvsdw880cBjATwfma9pPYMACLSE8A1AP4nowUlvmeg8M7eA8DWE3R1Zq0p0FVVt2e+3gGgazE3Ux8i0gfARQAWosT3nfk4vAJADYCPAGwCUKuqx2t0S/E18iKARwAcy+jOKP09M0CXBq37FUZJ/hpDRNoCmAbgQVV1Cs1Lcd+qelRVhwDoibpPfucVd0f1IyLXAqhR1aXF3ktDKehEGADbAJw4KrNnZq0psFNEuqvqdhHpjrp3opJCRJqjztHfUtXjnStLft8AoKq1IlIB4FIAZ4hIs8w7Zam9RoYDuF5ErgZwOoD2ACajtPcMoPDv7IsBnJuJXLYAcBuAGQXeQ1pmADjetmUcAL+lahHJnBtfA7BGVZ8/4Z9Kdt8i0kVEzsh83QrAaNTFGioA3JwxK6k9q+q/q2pPVe2DutfvJ6p6B0p4z/+Pqhb0D4CrAaxH3dns8ULfP+Ee3wawHcBh1J2/7kHduexjABsAzAHQqdj7NHu+DHUf0VcCWJH5c3Up7xvAYADLM3teBeDJzHpfAIsAbATwHoCWxd7rSfY/AsDMprJnZtAREgkM0BESCXR2QiKBzk5IJNDZCYkEOjshkUBnJyQS6OyERAKdnZBI+D/Tzi+kuuv7KQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# for the regular square windows case:\n", "aa = moving_window_array(frame_a, search_size, overlap)\n", "bb = moving_window_array(frame_b, search_size, overlap)\n", "\n", "aa = normalize_intensity(aa)\n", "bb = normalize_intensity(bb)\n", "\n", "# make it use only a small window inside a larger window\n", "plt.figure()\n", "plt.imshow(aa[-1,:,:],cmap=plt.cm.gray)\n", "\n", "\n", "mask = np.zeros((search_size,search_size))\n", "pad = int((search_size - window_size) / 2)\n", "\n", "mask[slice(pad,search_size-pad),slice(pad,search_size-pad)] = 1\n", "\n", "mask = np.broadcast_to(mask, aa.shape)\n", "aa *= mask.astype(aa.dtype)\n", "\n", "\n", "plt.figure()\n", "plt.imshow(aa[0,:,:],cmap=plt.cm.gray)\n", "plt.figure()\n", "plt.imshow(bb[0,:,:],cmap=plt.cm.gray)" ] }, { "cell_type": "code", "execution_count": 11, "id": "0d1e46f5", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "c1 = fft_correlate_images(aa, bb, correlation_method='linear')\n" ] }, { "cell_type": "code", "execution_count": 12, "id": "3c0ef682", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAaAElEQVR4nO3da8wdxXkH8P+DsY2lEDvEQJDNrSSSg6o2SIjm0g+IKApNUEilBDVNIyoh+UsrETVRgHxJW6US+ZKL1FxkBRRXimJoLgUhJQhRUJovtBDSJmCqYFToi2zcqOEmRXYdnn44e8z6ePfsbWbneWb+P8nye8573rOzs7P/M2d2d1ZUFURE5M8ZqQtARETjMMCJiJxigBMROcUAJyJyigFOROTUmXMubMsZ23TbprPnXCQRkXsvn/ifX6nquavPzxrg2zadjXfv/OiciyQicu9HR772bNPzHEIhInKKAU5E5BQDnIjIKQY4EZFTDHAiIqdmPQuFhjm+Z/fa3295amOmkhCRRQxwY7pCu+21DHOi8jDAjRgS3Ov+nkFOVA6OgRsRIngZ3kRlYQ/ckGUAD+2NM7iJysQAN6gryBnYRAQwwE1jUBPROhwDJyJyigFOROQUA5yIyCkGOBGRUwxwIiKnGOBERE7xNMKB+lxkw9P/iGgO7IFHMHVeEyKiPhjgA4ydKZCIKAYGeEQMcSKKiQHe09gwZoh3O75n98l/RNQfD2L2EGKubh7YbLZat011zbojasYA7xCqVxgjxL3fkadv3c69nvwQIS8Y4GuE/kofKsSbyrX6nPXAmTokFWv92srlrX6pDBwDb8DxWPtibB9uc/KGAb4i5k48V6/NQ+9wShm3PLWRdB091C+VoXeAi8gmEXlcRO6rHl8qIo+IyNMicpeIbIlXzPhi9rpDB8669/IULkPLGjO42fsmj4aMgd8M4CCAN1aPvwDgS6p6QES+AeAmAF8PXL7ocuhxp1peCFue2li7DSytk6WyEAE9e+AishvABwF8s3osAK4B8N3qJfsBfDhC+aKKEd7LXiLDu7962ev1Z2mdLJWFaKnvEMqXAXwGwGvV4zcDeFFVT1SPNwDsavpDEdkrIo+KyKPHX/vNlLIGF3qnnHMnXw0976wFdp3VchF1DqGIyHUAjqrqYyJy9dAFqOo+APsAYPvm83To38fW9RW+z9+nwmAJh3VJHvUZA38PgA+JyAcAnIXFGPhXAOwQkTOrXvhuAM/HK2ZcY0KcO7xtsc8XJ7KgM8BV9TYAtwFA1QP/tKp+XET+EcBHABwAcCOAe+IVM74+Ic4wsKXPhy6nMaCcTbkS8xYAB0Tk8wAeB3BHmCKls9zR2Xuzbei3JYY45WpQgKvqwwAern5+BsBV4YuUnvedPcTZNdbqwOtxCqKYOBdKBua+rHyuQGRoE63HAHcs1dWDMSd2GrNODGsqFQPcKUuXfk8J9KnrwfCmknEyKwrO0ocLUc4Y4E5Z73n2DfGxV2BavnKTaC4McIpmSE98SCAzuIkWGOCO5RhkXUGe4zoTjcWDmIV66bKtpzzefuhYlOWMvYhm9aKq+nNEtMAAd2zKwcJYgR0aQ5uoHQPcKQ9nejB8ieJigDvE8KYUYl7ARePwIKYzDG9KoandxbyPLPXDAHfEw87C8CaaD4dQnPAQ3oCPr9mcJGu4UtfbOga4A17Cu4mVQA9VhzytkSxhgBvnObybzH2zjJj1N8eNIviBQeswwA3LLbzrYgf5XHUXaz3aDhoyxPsp5YOPBzEpK7mcGdEUOjkHEY3DHrhhfW60TAsp6ylWsDKwxyul7tgDN66UhjhFjuFN1Ad74A6wJ96MwU2lY4CTSwzvONrqNed19oxDKE7kuANNWaeU9ZHzt6HVeuWdj2xjD9wR70MpOQVBzqf05bpeOWIP3BmPPSKPZe7D84cp5YEB7pT1QFyGdsxyWq8Dotg4hOKYpSGVVGGaug7mnhqAqI49cOcYHDYsrwC18oFKZWCAZ4AhbguDnObCIRSjhp6Pm3oogU6X85kqZAMD3Ii+4ds2yxrDO4yXLtt6yuPth45Nej+OkVNMDHDHrIS2x5s0rAb10NcNDXb2xikGBrgRnnrTFoJoTB31De0h7zUkyNkbt8/KHaT6YoAbtGw01oLcSmMeUi99Q/uVi6X1d2c/q2vfn71x/9ralPUPXQa4YVaC3FLj7VsX64J7XVive31bkI/tjVuq11JNPfaUGk8jdCBlg7HUWGOH97GLjp/yr+lvl/+GLrcJTzdMa2zdW9pm7IE7kKLBWApuIF54NwX1Rbt/BQB4DjsBAFuf29Jr2WNxCtc0hn7Dtbg9OgNcRM4C8GMAW6vXf1dVPycilwI4AODNAB4D8AlVPX1vIFcsNtIQH2D18G4KbeD14K4/fm5jJ45ddPy0EH/lYlk7Nh5C03pb3D7erdap9XHvuj498GMArlHVV0VkM4CfiMgPAfwVgC+p6gER+QaAmwB8PWJZizRX79tiYx267m29767wXgb3+99y8LTf3Y+3A2jujc8R4qusjsVaMvVMEk/12hngqqoAXq0ebq7+KYBrAPxp9fx+AH8NBnhQc4S3xcYacr3bwrve264H9x+/8fGTP//g5Svw/rccxP1H3t7aG08R4ksM8wVLY9Jz6zUGLiKbsBgmeSuArwI4BOBFVT1RvWQDwK6Wv90LYC8AnHXGG6aWlwIpeYcnykWvAFfV3wJ4h4jsAPADAHv6LkBV9wHYBwDbN5+XpqviUKxeRc7Bvf3QsdOGUZa941culpM952MXHcdzGztPvmY5TLLq/iPV8MlG88HM1Z731Mvux8p5m/ZR8voPOgtFVV8UkYcAvAvADhE5s+qF7wbwfIwCUjheGvqUibmaQhxoDnJgfZi3BXf9/erLnZOXbUlx9TkL5VwA/1eF9zYA7wPwBQAPAfgIFmei3AjgnpgFLUno3ndpO3tbiAOL4K2Pi6+G8/JgZdPvln/ftLw5lLYdU8ntLJQLAOyvxsHPAHC3qt4nIk8COCAinwfwOIA7IpaTBvLQ+NaZOj1uV4g3We2dr3ttfTmxjN2GXfXmvW3E0FRnHg4S9zkL5T8AXNHw/DMAropRKJrGamMbKkSI13VdKTnkbJIYwT1luw2pJ28TNs2hra1ZrxteiWmQ5ZszzD2HR8j5YNpCt88l8LGHSWL1tse8h/XQisXjejPAjRob4n0b4dAgrpclxRhhzIm9Up09spQyvNe9r8dAK02xk1lZ7eHWxdiB6hModd2It8/v57blqY2sgsVaeM+9DJqmyB64p4Y5tCe+rme97n2mzMyWIlAtDzPFNud6eziQV7IiA9xbQ5wa4rF3+BRfuUOsU5/yxqy7MfWV8kPL86x9uSoywEuwDPG5e2tz7LxD1mlqeWKNvXsL7yE4hj6fYsfAvfGyw8e+ScGc4b36XqHez8u2nMpjmb1hD5yiiNELSxXe6953aEiVEt5L7I3HxQB3xOOBu9zvJhR7Wd62dxveAzQODqE4k9tpdDHkEHq8Xyb1wQB3ikGer1yDO9f1SolDKM7FvELRsyFf2S1N/pT7duRQSlgM8EwwyLtNuVipLtYVsqVgiIfDAM9MnwOdTTtPbgHiKSByq3uaDwM8Q2Omxpzzys0YPAV2nce6DoG98DAY4IUYsrN4Go7xGgIe6jY2hvh0DPBMLXvhU3YQ6+edp5pEayrLddqkbb70ENPwMsSnYYBnLMSOYTXEQ51hMidLZemj60YXL122Nflc6qVjgFMni0MqXT03C5NP1Vmqu1AY3ukxwKk1XFZDy1qQzzFtbolDJkC/28xRegxwp8begHXMzW/bgnzo+8UwZfmcxyQ9jn9PwwB3ZshddUKNX6+bUc5SmA8RMzg81UOT0nvfnmZQZIA7MjQYQgdJV6/f6gHPJjHOfvCy7l3qY9sxz0CxGJCrN++2WMa6YgPc06esZZ5DK8QO6nn9+4hxoNLiPrfuxt0Wy7tUXICvbqixY8lzyz0oUhl7015uj+G87leWQ7yYAA9x8C4VhsU82rY76386K/tS3ZQbhVtRRIB73gE9l90r1nlYFoNvzDa2GOK8oYNhDBKbOz/5ldudjhjgZFL9jkMMcd+sBKaVcoRUxBDKGKlDI8fG1nXOeOo6p3hSDz/kuD8BTgO8z8aYu7GMPZuh671y0HUBEM2rre5TtDve7WgaEwFusYKnlCl1b8OaEPXh6SIhr7xeVdsmxjpYO0NNVHW2hW3ffJ6+e+dHTz6O3UhS7/RTNnIOO9AYXXVWar2ENte3xKHz84QKxjnayZwh/qMjX3tMVa9cfT5ZD5w7IjXp6q2n/lBObd08JUOumgxxs49QQtx8pG6u9mGhNz57gJe0801plCUHFYegXjdkYqmm164LdUv17C28V5eZqh5nPY1Qz9oy5+JMBKCFMuTGSujEFmJWwK73yKl9plyXVMvmeeBG5bRjxVBKiFM3CxfnmO2Bi8iFIvKQiDwpIk+IyM3V8+eIyAMi8svq/zfFL65PQ+dcSN0Yvcg5xEPNyZ37bc9y31e68qDPGPgJAJ9S1Z+KyNkAHhORBwD8OYAHVfV2EbkVwK0AbglQ5mLl3hhpHrmHNmBrX4nVkeizjp0BrqqHARyufn5FRA4C2AXgegBXVy/bD+BhMMBbtR3osNQQPSr5YO9SCYG9ZG1bpwxvYOBZKCJyCYArADwC4Pwq3AHgCIDzW/5mL4C9ALB1644hi8uWhUY45x3daZg+wyclhbZVqcMbGBDgIvIGAN8D8ElVfVlETv5OVVVEGq8IUtV9APYBiwt5Su8xeVh3a3ef75Jbm2q7pRlDO705jrsMac+9AlxENmMR3t9W1e9XT78gIheo6mERuQDA0SEFBPwERCje1jfVdhqzk+QW4kuxQ9vbgeCQ29nyuvddz84Al0VX+w4AB1X1i7Vf3QvgRgC3V//fM6aQgL9gK42X7ZRriHu1bltMvQq0hO3cZz379MDfA+ATAH4uIj+rnvssFsF9t4jcBOBZADdMKehSCRvGKw9BXsrObc3QOp96GXop27mrfvqchfITANLy6/eOKNNavB+hbdweeZkyL3uItjAlyKeEuOXhk1XH9+xenCbSwMR0sut46PWVYM76nzK3BNvJdJ7mExnS4fMU2n2ZD/ClUr4yzWFoQ/ayQ7N9+BZq3vicdLVpV3Oh5LZxrEt9Wf/QKQiISuMqwAGGeAhtYbcM7NTBXdenHFbKStNZanupBTmN0CIOp0znqf4sTJxP8+I278dlgAP9D25auoErTWPpBgQ0jxi3W/MgylwoFvE8Uh9CXRLOnlm5uO1P5z7Ap2CIx9c0MdPqc2MCnb3xcuXeK48ymRXZktPUtC9dtpUTNdEoufXKh+7DxQa417Bb8lL+7YeODb7ZLtFQOXwjG7NPFxngXsIvFzHC2vvOWpIht4ebenzEa7sYm0nFBTjD2z+vO2nuQtzHc/keY4Pc45DKlExydyHPFAxvysmWpzZMBNVLl20NdhPmUO9Zyr5eRA+8lI1ZCo/bc8qUyV0hnXrCt7bjHCHel9bLOsA97uie8WyS0zWFb1eYW+hVD7Ua4n3bQYhTSpt4rMMxsgxwBnd8bT2uqWOYORgSHjkFzZhtzgPc02QT4B5Ce0rDsrJ+fb8ql3Yz3tShYaV9pJZ6O4wx5YJC1wHuqdGGmOc49fqOHefMdWjFY1jkzPP2GLt/uwzw1EE2RIhG5Wl9iVLwHN5LYybecxXgnoIsh+ESWs/zhSPkx7o25iLAPQXamB3aw/rFOE0sBwxxSsl0gHsItjruyGWKGeK5z7xH05gOcC+83UF96Lm3lnvfVg6Qhr6Eu6lteOvQzKnUb0KmA9zCmRcxWAnupufrYWg5uIHXy2fp3POpF+bk2N7nUmKImw5wYNolyHPo22BSln1IEHu5JLrtRhEWQnxVn9v+WWzbHpUW4uYDvC7m19NYvAS3J+vWy2qIr+MtvK1/4JQU4q4CPJRQvfp1jSR1A7cS3qHDtM96eQxxD+rt3fo3Y4/Tyo5RZIDXTZnJbd0nfepeSqwZ4oaWIcZ7dq0XwzusPrMhWgxxYHiQx16P0B8oxQf4VNZDvC6X4FsX4l7WwYs+gWM1vOu6gnyudQg9vMMAD8ByiNflFG5jpy+l/nIJ77rVIE9R/pDDO0XdkWedqZW5riHkPg6XyjK0Gd7h5Rjedcf37E5e/hDLZw88oK6e+OpraTqGd3jscMxnam+cAV4zdrhjTOUP/RsGPsU0pQ33Oc+9Cdv068aOjTPAVzRVooXbXqW+7+FqOYD0ZaHpQrTjse+R68RvY40JcQZ4D5a+UqYM8tV6YJj7ZKk9D9W3g+XV0CEVBrhTcwZ5n8bEMLfLc2D3YeXbaUh9e+MMcOe6xhe7GkGfRm/lHFrqJ/fAbmPplN0QTlmXI82v6QxwEbkTwHUAjqrq71bPnQPgLgCXAPgvADeo6q+nFZfGWtdoV383paec086Rm1JDe1Xoeujb5lN9A+3TA/8WgL8H8A+1524F8KCq3i4it1aPbwlfPFpn7OX/lA8Gd1xTzzCLvb91Briq/lhELll5+noAV1c/7wfwMBjgRLNhcPsQO8zHjoGfr6qHq5+PADg/UHmoB/aiy8Xg9ivGGP3kS+lVVQFo2+9FZK+IPCoijx5/7TdTF1e80sL7pcu2Jp9VkSiULU9tBP0QHhvgL4jIBQBQ/X+07YWquk9Vr1TVK7ecsW3k4sjC3A0pMcQXSm4DOQkV5GOHUO4FcCOA26v/75lckoJxp2y2Gtq8UcNCSXecyd3Uc9j7nEb4HSwOWO4UkQ0An8MiuO8WkZsAPAvghlFLLwhDeph1N2BmiDPEczM2yPuchfKxll+9d9CSCsbwphgY4vkZGuScD3wG3MmG237o2Gk97abnSsfOQZ76jpHzUvqZWLknnzcM7G71NsPOQl66ticD3BjLd7on+2LOZ29hWuWldVNElIQB7kiOs65Rem09+CETndWFCtMh7XzINKxD9x/LHw4cA3co9MUAMS3L6qW8JbO2jWK1myHvaa1OVrEH7ljM6V2tN1wKZ919XMdOmBay/fS5icPQ5eXSvhngGRoT7HM06Nzma/au7406LIT4qlwCeCoGeEEsNHqGuA0W2gJNxzFwIlqLYW8Xe+A0O/bC05l6B3lLpxISA5wS4SmR8wsRtAxrWxjglJS13nifMx68KS10S7oylQFOyaW6IWxbGbp+5yXQcw8voHtbDLnAxyMGOJkyd1hOvWltG4sfRJQfBjiZZTG8h7x3n/JPHbJhYJeNAU5FmuvCJWD4mRtN30IY1NSE54GTSTF733OHYX15Y+b3YHi3K32aZvbAjSvpiPocUtUhtx3FwAA3ztqOz6/zZE2fs5hybbMMcGrVtjPM8a3A2vnh5EOuQd2GAU6niTWRPqXHO9nkhQcx6aTje3aP7vVO+ds2DJdw2rYPv+X4xh44mcahlGHGzt29xA9NXxjgVIwcwyn0hxuHxHzhEAqdNHWnjXlA08J7WBPzm0mMITEKjz1wmizHcKSFOU4b5YHV8dgDp1NYvUpwynJyDIQ5e8exe/pNz9X/UTsGOJ2m7+XeKS9JHyKnEGCoUR2HUKhVTj1X7wfnUob2HHeX9759UmGAkytTTyv0EhRWetkehshKxgCnIqUOcisBvQ5D1T4GOLkT8uKeOSfnYmhTaAxwcsnbFZrWysqgzgMDnNxad8eboe8RC++NSTExwMm9qfeVDC3XM0bIHgY4ZSlFkDG4aW4McKIApoZ3n7vKtL2eysUAJ5ooRM976Htw/hACJl5KLyLXish/isjTInJrqEIReWHl7BJeYl+m0QEuIpsAfBXAHwG4HMDHROTyUAUjouEY5GWZMoRyFYCnVfUZABCRAwCuB/BkiIIRWWc5KFevNG0qK4dd/JsS4LsA/Hft8QaAP1h9kYjsBbC3enjsR0e+9osJy8zZTgC/Sl0Io2zWzZHUBTipu37slHVuNtvOcBc3PRn9IKaq7gOwDwBE5FFVvTL2Mj1i3bRj3azH+mmXe91MOYj5PIALa493V88REdEMpgT4vwF4m4hcKiJbAPwJgHvDFIuIiLqMHkJR1RMi8pcA7gewCcCdqvpEx5/tG7u8ArBu2rFu1mP9tMu6bkRVU5eBiIhG4D0xiYicYoATETk1S4DzkvtTicidInJURH5Re+4cEXlARH5Z/f+mlGVMRUQuFJGHRORJEXlCRG6uni++fkTkLBH5VxH596pu/qZ6/lIReaTav+6qTiookohsEpHHReS+6nHWdRM9wHnJfaNvAbh25blbATyoqm8D8GD1uEQnAHxKVS8H8E4Af1G1F9YPcAzANar6+wDeAeBaEXkngC8A+JKqvhXArwHclK6Iyd0M4GDtcdZ1M0cP/OQl96p6HMDykvtiqeqPAfzvytPXA9hf/bwfwIfnLJMVqnpYVX9a/fwKFjvjLrB+oAuvVg83V/8UwDUAvls9X2TdAICI7AbwQQDfrB4LMq+bOQK86ZL7XTMs15vzVfVw9fMRAOenLIwFInIJgCsAPALWD4CTQwQ/A3AUwAMADgF4UVVPVC8pef/6MoDPAHitevxmZF43PIhpkC7O7Sz6/E4ReQOA7wH4pKq+XP9dyfWjqr9V1XdgceXzVQD2pC2RDSJyHYCjqvpY6rLMaY4bOvCS+35eEJELVPWwiFyARQ+rSCKyGYvw/raqfr96mvVTo6ovishDAN4FYIeInFn1NEvdv94D4EMi8gEAZwF4I4CvIPO6maMHzkvu+7kXwI3VzzcCuCdhWZKpxi3vAHBQVb9Y+1Xx9SMi54rIjurnbQDeh8UxgocAfKR6WZF1o6q3qepuVb0Ei4z5Z1X9ODKvm1muxKw+Fb+M1y+5/7voCzVMRL4D4Gosprp8AcDnAPwTgLsBXATgWQA3qOrqgc7sicgfAvgXAD/H62OZn8ViHLzo+hGR38PiQNwmLDpfd6vq34rI72BxcsA5AB4H8GeqeixdSdMSkasBfFpVr8u9bngpPRGRUzyISUTkFAOciMgpBjgRkVMMcCIipxjgREROMcCJiJxigBMROfX/jFRh7ZxO0HwAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.contourf(c1[2,:,:])" ] }, { "cell_type": "code", "execution_count": 13, "id": "59dd5db4", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "n_rows, n_cols = get_field_shape(frame_a.shape, search_size, overlap)" ] }, { "cell_type": "code", "execution_count": 14, "id": "385c1e95", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "u1,v1 = correlation_to_displacement(c1,n_rows,n_cols)\n", "x1,y1 = get_coordinates(frame_a.shape,search_size,overlap)" ] }, { "cell_type": "code", "execution_count": 15, "id": "7ba344b8", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [ "x,y,u,v = transform_coordinates(x, y, u, v)\n", "x1,y1,u1,v1 = transform_coordinates(x1, y1, u1, v1)" ] }, { "cell_type": "code", "execution_count": 16, "id": "dda78bcb", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAKrCAYAAADyAksxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAADcZUlEQVR4nOz9eawk2Z7Y931PbBmREZH73au6q7tfv50zb4ZvZkiJEiXSgijK9pCwJNN/SIIggDZAAZIhGLb8jyXAAmwDEgH9IQIUKIsyBMuESIhDeixxNJwRNZx5M69nea/37lq69qq770sucfxHZN6bkcut6qpblXlO/T5Az1TevNXvREfEOb/4nd85obTWCCGEEEII8aZzZt0AIYQQQggh5oEExkIIIYQQQiCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAGAN+sGALRaLX3jxo1ZN0MIIYQQQlju93//9ze11guTvpuLwPjGjRt88MEHs26GEEIIIYSwnFLq7rTvpJRCCCGEEEIIJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACeI7AWCkVKqV+Tyn1E6XUx0qp/6D/8/9CKXVHKfVH/X9+0P+5Ukr9J0qpm0qpnyqlfv4VH4MQQgghhBAvzXuO3zkD/ozW+lAp5QO/pZT6//W/+z9orf+bkd//F4D3+//8EvDX+v9fCCGEEEKIufXMjLHOHfY/+v1/9CV/5ZeB/7L/934E1JRSKy/fVCGEEEIIIV6d56oxVkq5Sqk/AtaBX9Na/27/q/+wXy7xV5VSpf7P1oD7Q3/9Qf9no//Ov6yU+kAp9cHGxsaLH4EQQgghhBBX4LkCY611T2v9A+Aa8ItKqe8D/x7wbeAXgAbwf/w6/8Na67+utf6h1vqHCwsLX6/VQgghhBBCXLGvtSuF1noX+A3gz2mtH/fLJc6A/yfwi/1fewhcH/pr1/o/E0IIIYQQYm49z64UC0qpWv/PEfDPAZ8N6oaVUgr4C8BH/b/yK8C/1t+d4k8Ae1rrx6+g7UIIIYQQQlyZ59mVYgX4m0oplzyQ/lta67+vlPqHSqkFQAF/BPzv+r//q8CfB24Cx8C/ceWtFkIIIYQQ4oo9MzDWWv8U+LkJP/8zU35fA3/l5ZsmhBBCCCHE6yNvvhNCCCGEEAIJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhLhSuttD97JZN0MIIcQLeJ7t2qyz//CAJz/6CgCloPpei8WfXZlto67QZ3/rpzgOOJ6Dn5R4+3/2/qybdGUe/e59Dh/sohyFUrD6p96l3CrPullX4mjzhMc/uovudMnaXdLrNVb/iRuzbtaVufvrN+kcntE764LWfOt//YNZN+nKbN/aYffLDbqHp/ROO6z8E+9Qe7cx62Zdif37ezz57dvgeSjfo/puk8XvL866WVfm1q98TK+ToRyFE3i89y9+e9ZNujKPP3jIwd1tAJSjWPmTN0iWkxm36mocPDrg0W/dRrkOjueQvt1g6Qf2jONf/u2fojU4vosb+rzzz39z1k26Mk9+8pT925uQZThkLP3SDdJr1Vk369wbGRj3uhm9k/bF57PuDFtztbKeRrc79IDerBvzCnRPu4VzpzM9w9Zcvc7u0fmfbbouAdoHZ3QPT88/616Gcu2YtOqedDjbPDj/3Dk8m2FrrpbudOl1Mui04aRN93R+BrCr0D1uk3XzDL9q29Vr9s6K/eXgOG2QdTOydt5H9rCvv+yedkBD77Rj1XkD6LV79I7ysaBHPss2T+wYlV6SUrNuwdXJesVA0fHsOsWjgbBN584Jis+pumNXR+8GbuFz99Se4wsqYeFz5+B0ym+aZxB8DIxep8bTF32KcizqUCaxqMMcLVey6dxl3QyGhjrHd6f/soG0KsYlSs9X4G9X1PS87EoyFvQ6I52FJRm5gdHA2HHt6QxHA47srDOjlrwabql4fN2z+coSvAybA+PRBzTl2xUYF/oUiwLHSWw6vNEkkE1jXdYp9o3KszswJpPAeOb0SGBsU2cxGhi7vmWneOTk2ZQlcD2Fdi+CDusyxqOBsU0Z4yQodPaDaUIbWJ0xtnkwgPHjs8hoxtimJElvpKTHtplfnJGMcTZfSRLL/ms/H5un40drkWwKHGH8wdKm43McioFxL5u7J+mXMZYRb9sTGHseqLB0/rl9mkHHjoy/1Rlj2wPjURYdn2SMzSUZ4zlkc19oe8bY6lKKkcC4lwFde4JHmzPGAF5yUU6hNYWFhiazPWNc6FFsGgwsZ3ON8WjGeHR9hvFGM8ZSYzx7Fs8ujWWMrZuCsbiUAkAF/vmfdYY1WUcYD4xtW0U+HBgDtPftCIxHV4xbtRDI5iwJjPeXFh3eWJLEorHO9oxxNhp6SsZ49mzuC20PjMemzywLjIezcZllGWMvtDswHl2A1z6wY8u24VIKrZTdgbFjV385lgSyaLCzeawbDYytuudAMsbzyObAuNe19yka7K4Ph2L9ptZ2LcCzucYYwItLhc+2lFIUrkHXsyt21Lq4S5FtHYrFxsYCi5Iko4GxbaUUUmM8h2wOjG1+igaKJ0/Z1RlCsZQC7NqyzfaMcalq4ZZtvV5hlkbbGBgPs2kwALtLKUZrjGXxnTG0UzweyRjPIZs6C9sD4+EsgQK7Th52Z1W9UrEzHF1gYrogctHexYNN9+jM/AUN3W7hECQwNpxFxzc21lm0EHv02CRj/HrZ1MU9N5v7wrEbyrPo4CjeP0ph18ljfCssqwJjyzPGQQCZf1FO0ekAZ4bVGWtdvMm63cJH2wNj22agTH8uu8xYKYXNGWOLjg2QGuN5ZHNnMbpdm20Z48LrWy0MjG3OGDuuAvci82HTsUF/L+Poopyi3QZODSunUAr92eewt5d/7vUoxB+uezGmZebv1TwaXGns6k9sHuyszhi/aTXGvfmaPbQsano+b1LG2LbAuFBK4WDXyQOckr01xlAM/HudnnUDtw17GW/teWz+6CbZFzfh6Ag9mjE+O4F79+j+wU/Z3pyvTM/XZvNgMIFVGXFt0ULzkcDQ+l0plEIP3WvzljG2aKf2S2RZfhP1s1WX9oVa5xepZ+Z/mrHXZJrcWUyg9UVOx6Iu/txoB2hbVtUteWQn/fICnZdTuKF/+V8ySFAJGQ6FOwenxnWyYSvh/q19dnb2WFjYY7hL8XY36f70CdvbsN1Jqa6Upv+LTGB5KYXNbHrLa299C/fsGNbWwPfH9g4/zxhnGezvQ632+ht5hbQmz2zp/nFKjfHrd3Ss2Pvdz+DxY+j1JieptIbtbfZ/9InRWbrR7dpsevNdloEaLqUwuCOcxnFV8bXQFm3XBuOlIrbVGY+95MPAvYzLCzFK5VtoP34M3aHuUPW63LoFOzvQrrRoNGbXzqswmkiwLWM8WipiE5te8LHbTbj/ky3af/ARPHky/rZJ3YOHD/NZmrsHM2rlFRuqM5aM8QyUQsXtrZjk0SNay+s4WVL43jk5Yu93PmHr4SkEAZU4mlFLX95oR29TYKw1xRpjCwNj182nq1Uv7xizbpY/EViy4mm0Vq572iWozqgxr0CQBGilzh/gTCylcNKYKILj4+m/ox0Xf6FGZG5XmZNSCmPZNDuaLkY8PXI5vNmjsfWQ3si91/vpxzw+gL1dKP98A8OfR4G8zvj8atQ6/2dO7j9zr6SvwfOgtFjl9BQefNVl//5u4fuHXxzx6PYpZ2eQrJk9ShcW3ym7VrO+CYGx41DY8iuz/LXQ3bP5WnTxskqhQg/tTNE+6Zr39kLXJW6Gl/5Kt9qg0bKgb3nDAmObjO1jbPB4EJQUbjWBfOJ6bDObe/dgdwd6fkj9WjybRl4hrRlP9sxROYUFPdvzSVYrF8XeI33h8PmovFV7bW16FYbrrlzLFqeNjWEGd4TTOA6FUorh10KfL+YyeMHaaGBsWylFEEBWuggqux3Od6bIMjg6mlHDvqbywuWDb7fWol5/TY15lWwPjA3uK0Zt3z/i7PTieIZfPDOaBDo9np8g63mVF5Nn/o5qNkjT19CY10APl1IoJDCehUrdJYsuv/CC0CFsPfvinGfDT9HKwZopeOjfN4U335k/iHW78OSr0/Mpd9cF7RVfC023y+52xhe/8TDfKcCg4x7t66bVGHe7eVak3X5dLbs6w8dYKlHIGA/2Mj48hE8+gZNjbUQGOVpIpnYdWSkiapQJL08qz6deD31yUd4yVoM7em/N2TZSz3Lw6IDj7UvKd/rHpzPNzp3d19OoK5IdHHHrVz7m4Yfb+S00dOM5Kt9qb+dpmy9/4wG7H96fXUNf0PMExpUbZhZRjD6fnS++GzZHgfEbUWMMEMfkKzmPpxeuxysV4wLJzccd6os+rtvvw7ORjPHQ8cxRCc9z29toU10IgMHNdXGHjZ6qbjvDC8w6f54Hu4ceu7/9UyqrCe5SC4Zel3l6Crc/2OZk8x4BbYK3a7Nr7AvIMvjsM2g2839GX/LRPe2ysQGPHoE6O+WtBQ2YVbh6/34e6zYaUK2O72X86PYpjx1Q7TOiszuQvD33u96oJKZchsPD8e+6tRZLpmaLHYeHv/EFpcUqze+vTN6VQmu669tsf/IEvvENFq+bs1WW68LdX/2Y8vUmCz+zMva9zjQbn22x/fFjwkaZ+ju119/IF5SsVtj4g/vs/+QOe188hbPueY1qphVf/IOv6G1uo7Sm8s++O9O2vohkKeapcqYuROtFCQ1Dd4HZ3ISTE1hZAb9fKaiH9rSnnzHe2cm3UL/x9myDlfnuna+QUlBeqXL2aPqTpIllFCdP9tj8vS1aP7hGZSVmeNPRQcb46AjWb+6TZnu0fu767Br7ArY/ecLmWYfFH6zipdHYrhRaw+7TMzY/esLimk/1O6szbO2LaSx6PI1r7D3YgQcHDG9edngIHG7jAOUK/Sc8c3iuJt57zMO9Jo8elSj3il3O5pMuZ1kPf/Mx9c46/Mx3Z9TSF7fQ0nzxoy12N6u4JR/88Lxj1Rq2Hp3iJpsET+4Tvp/laeV5F4aUE4fDw+IgrZWiW2mYuxuFUoQLKRufbbJ7a4vqSrnwtd4/4Omvb7HztE0viHjnlww4V0PKiwlu4HJ8f4uvHmzjjJSb3fu1z893ukmvm9VXhrUQL/LpnnRQJyOr07Qm29hCAX7Jobxi3lqhUqhwKzHZ3uTknbfUpFye+NXcS1N48PkRO/c1zbeTPOM/lDE+2IcnH2WcOLCW7MHNdXj//Zm1940JjAEqCyWeBiFOe3yqyfMx8maKVyocfnyXjf/pMzabddTQ3kq9Lnz5413ajzZxTo+59k+b1RFCfnxbv3eLB/9gF2+hmKY6O8n44r+/Q7a5jVKQ/tK3Z9TKl1Ovw8NaC+9g59LfK8cYFxijFPXyKccffkQvrnBaTgmGvz46JLr9MarbIVp2MXF+vhwrKuqQsy/v0osSsrhYBOgd7Ob/eOR7NpswK6VUXmf8pDhI95IacdUjCKb8PQOkaxU2Pt+m29Fs3SsWfff2Dtnuv/TPadbNC0SUorxc4eDeDkprdK+YER8ExVopKtcqs2jhS4mWKhx8tXXp7yRrVTPusQmihYSjCYGxVoraDVOnafJuPYh9nI8/ZPdJTLexhDuUEX76FDr6kNLhQ2pLB7A62ydvM6+eF1StQi+tTfwuWYznfnpzkqQRkPVrGvXWDk7nokiz04Hu7Xs4p8fEMfiL5t1Y6Wp6vmiyu1EMHDsnXbLNbQCSuo+TGhY09vl+/gCQ+Zdnp+IymDdSQ3qjlQeFR/sEGw8L3zmds/OHubBp5vkDaH53CQD35BB/8/HE3ymVMCNb3Be14uE3eAN5GYWx2eK+0kLlfDr3Mun12itvy6uQXnt2gidsJmNlTSaIl5+98qzytnnj3EC8NLnOuJdUaSyYU9IzSaUV0Iti3JMjSg9v4+0VH3CC9QeUOgd5bmTGKwzfqMDY98FvTe40TO0EgwC8+nN0FkuRkdm4KHFx0mcvSkgMPX8DjQZ0a82p33t+/la1sUjFAKpaodJ6RopRPd/ik3lVXY5wm7VLf8e0wJgkKTyHac8ni1Pzd6PwfeLW5XXsWRDmfaaBkrXqM18LmlyrvZa2XLVk9fIst6llFAPJclJ4VfJAuNo0quuYpFqFXuXyp+rzob4y29mMNyowBkiX48J2WJDHGvGquTdTtHT5RaQUpG+bm+aJV559fJW3zB6t63Xo1VoTO0XoZ4tNK6MYUn13etAPEJYwNuM/UP/O8qXfByXMejiN48Il1621SCvqubKt8+5ZARa12qyTVi/MLXmEjcvvpcp1M8c7v+zjV6Y/sKTXzC2jAChFDu7IC8a061F728zzNSxJIKtePk6nKXm2b8a1WuZeQS+oWlP04mKnmDQCVNnM7AA8e3opScBtmRs4PuuhJal5OBVzs42QP5xVWz69ZPKxlg0PjMNrrUvfkhZGGH18AM234kunAEPTMsaeR7l+0d5utWl+trgvWa1cuug9vlY3bgefYZeVU/iViLBq0HU4orw0/R4zuYxiIFosHl+vUqfeMPhi7HMcSBs+vfLk8+c4/XFuDp5I37jAOElAV2vFnxn+tru07pGVpkcdleWyWQPyiHQxKrwNblQ+dWh+x5GXU7Qmfmd6YEwQULk2vcMr180sExnmOFD91vSscRBg3H1YasR4HvlgVipZExirNCFOJw9/mR/kfabBLguM47Xa62vIKzBtBjEIHaJls8dyGC8pK19rmLj8aaJqNX9r5iRp2h/GZ1xGAW9gYKwUlJcv3oLnOObWFw8EAbi1yUGH65pdRgH92vDmlJtFmX/+BqpVIK2MPQR4fv7K0EtTrgaovtua+vxi8sK7Ya13K+hoPKjy/H7cb1hgTJzvZ9ytNklTI9cnT+Y4U2fasrSW34sGK9XLeNHkZMLzLM6bZ8lyMjERYvJuFMOS5YvAOPNL1K+bPRs6rFrNd7aZVDJ4niiWjPFsVBsXb8FLKg6qMvsT8bLCxcmBY1oB1TA/zTNtQUWSOrj12T9hXgXHgXpDjWWNz3ejMDwr7jZrVOrjWWHXhVLTjs4/CCD+xviLFUoBeVRpWlY8SYgrLr1K3fjdKEala5P7jXClbtxpmmTSTKhT8okXzX4IdQOX0oQaalPfCjcqTDyccr4WIas3qNVm256rFAQQpR7ZSDkrqj8hGobMwyKGNzIwrlQ4r+VMVivGBxwwfTVrZTmeeSH7VUhXJz+8mL7YYtSkcgrjyygGHIfqO+ODV2TL8fW1vlEjC4qL7EqmLbwbiCLK15so17FqgIZ8Fmr0lGjPp7pqx7U4KTAur9ZsGO4oLxcDqyB0nrkI3SRhv844eatp0/AG5PFXd2R3iiTuD+NzkC2GNzQwLpXAa9WsmoZPay5ZWJzC9XyIr9vxFJ3UPHR5ZMCy6PwNpCm4UVBYIBqb+GKPKeK3W2MJgajsmBk0ThHHELxVrDU2bqu2AaUIbqyyuGhesvuZwpC4XkwadNMatboFkSOQrKQ4bvFYTC+jGEhWigGUbQmSeDGhF8XGvgL6MoP3SeihN98lc1RGAW9oYAz5W/DiRohTt6OjKJXAqRafmKsVsGW1jONAuFQ8V0micBt2nL8BpfJTNsgae35/ZsnAF3tMVC6PvYY3bMZWzNoMa32zQeZfBF3GBsYArsvKeHWIFUa3bQsWajZMsAGgPLewkEu5jpFvu5skXoxR3sWTmi1lFAPxUgKN5jysQ7tySQKu7xR2YEoGl6kExrNVqUDyrTWLVpOMb2NTXUvmol7nqoyuRk5WKxamsfJyil5aQ7teXl/smvmq5GkKexor8lcPW6ZWV6ilpfPPRgfGWHmbAfnahcGxacelsjYfA/NVGd7qMlyqjmWQTaUcRdjKo6lS5ExdY2OqqBrQ/GbTtnwB0H/vQAV6/d0pwgh8j3xx+ZzEY29sYJymUHmrNutmXKl4MT6fnghDKK3Y9RRdWS5fvJzFwjKKgSTJd6HoVpv21BcPCZYbJGne44cljN+DehKloPmtFtr18Pz+LK9FDzfWSNPza7FnURnFwHDphC1lFAODOuPkWs2qMoqBpRX7jmmgWoVeXEE77sXb7uYkWwxvcGDsOHPzcHJl0qpDVs6vskoFbFstE5UVql8uEpf7deKWGizCs6m++JznnT+U2vBij2laiw7ZwlL+Yg8wOmNsLdc936VB1WvWVCwN+GlIqVJCK2Xs2+6mGZTB2PBSj0ksjPXPVSqA49BLa3O1TduAxf/p3zxhyPnWc5W11KoyioFBOUW6ktj3ZDOk0cgHNb+e2FNfPKTybgvXhXKtZO15dF2ofmOBUtk1c6u2N0SyWgHHnvrbUfFqlVIjwY/sus/KjZByI7SujOJN4Pv5sOYtNi7yBcn8zBzadacIykspzh74S3aVUQwkqxWOPob0LTuzBANRBIuLgNeyMqOqKimVVmDNiz2mWVxxOT5cAO9g1k0RU7j1CuHSCbWGnXmi9FoVVbJkReGIpR9etzu1arFqFXSawpaXbyk7RwmS+WmJuBLxQhl/07OujGIgbfjsN8tWl1EMLC4CNKzbsQEApah/o0kptrsLKpUg+OYiPO7OuilimnKZynsL8zSTe6XKSyleYmd9u2SLzZW/XVJBVp+7hxu7R6U3UKWq8L9j124bw3wfWn9sxYqXljxLHg9bGBT3ldZa0LU/YFSBj7X7ndlAKRpvp1Y+fwL5XtSJ/f2lMMv5y1x1A3q9WTenwM7o6Q0WhkCp+czfM1lyrTbrJoirEARvxAMO8OYcp6GsDYqFmFPn91ySQJbNtC2j5it/La6G9PJCCCGEMMGclVLMV2uEEEIIIYSYEQmMhRBCCCGEQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYztpPesWCCGEEEIYx5t1A2Zh7842G39wH6cc4iUhlRsNatfTWTframQZt//OH0EY4sQRfj1h7QcLs27Vldn46ClHj/bQGnSmWf6F65Rb5Vk360qcHbRZ//ApvaNTuoenRAsJa3/qnVk368o8+oMntHeP6R2dos/afON/9bOg1KybdSWe/N49jp4c4IQBbuiz+INVwkow62Zdib27u6z/+C5uHOInJdK3G9Teqsy6WVfm9t/9EDwPJwzwkpBrv7Q26yZdmfWfPObo0R44Dsp1WPy5NcrNaNbNuhJH60c8/eA+ynVQjiK5VqP1HXvGutu/+hkAju/iRT7X/tSN2TboCj3+8QOOHu/jhR5e5LPwMyuUquGsm3XujQyMOwendE66cHLI2dYhpUYM2BEY9w5PODvVcHoCuyd0T7tgUWB8tnfKyfrB+edeuzfD1lwtreHw9vr5Zy88nWFrrt7xo106u0fnn3snbdxyaYYtujrt/VPa+6ewn5+z7I+tzLhFV6dzcEr3tEv3tN9f1iLAjsBYd7qcHbSBNnCMl0aAPYFx++CM062he669PMPWXK3uaZez7YtjC1I7+pKB4WNzIzsesgfa+6d09k/o7Oefm99dmm2DRryRpRTdw2LA4SXz86TysnqHJ4XPKrYjmzowmmC0qWrEjwO0urgle0d2BcZuXLzPOgf2HF922i58dkN/Ri25ep3Ds8JnP7WovzztFD6rkl0BiO5lhc+Oa8cMDYwfm3LtCWdGEz6O786oJa9G97Rb+OyX56u/tOdK+hpGAw4/sedJMzs8Lnx2YjumzQaUU+zYdWZPZOy6oKKLoKNzlkGnc8nfMMvoA6hVgfHZxXnSrocX2NO1jp4nmzJzvZPiA41Tmq8B+mXZHDxmvWLf73j2HJv1gfHJ0LimwAvnq3jhmVeSUipUSv2eUuonSqmPlVL/Qf/n7yilflcpdVMp9f9RSgX9n5f6n2/2v7/xio/h69G6EBhrx527p5WX0R3JGFsXGFucMQYKpQVa25U1Hg2Mu0dnU37TMN0uve5FAKI9H8eeMXrsPNkUGBcGaOybsh5NHFgVGHdHsuEWB8ZuYE9grDX0hhIJXsmbu7Umz3MlnQF/Rmv9s8APgD+nlPoTwP8d+Kta628AO8C/2f/9fxPY6f/8r/Z/b3602/S6F51FVgrx7YmL0UcXgbFWCqdsz7TnJNqyyHg0eGzvWxwYH1pybJ0O2dA4pgJ/3vr5F5dlheBRByX8wJaDGy+BsT5j7Fh07rr2HpvNGePOaQ81NG7PY9nZMwNjnTvsf/T7/2jgzwD/Tf/nfxP4C/0//3L/M/3v/6xSczRMnJ7SHS5vKYV485XFf3HtNt2hG0qXIjx/fv7TXwWbOr9JRus3rQkeGS9ZsubY2m2yoTHaCe3JOmbHp/SGxmgvtidbDBMCY4vOHUDhwsSurOpoNtymY7M6MD4uztL40fwFYM91JSmlXKXUHwHrwK8Bt4BdrfUgxHzAxVLeNeA+QP/7PaB5hW1+KfqkGBiPLggy2vFxIXOVlSJce+6nyexKGFtdh+sHisy/CKy6px0KUZehRhdw2ZR1PNsvllHYFhiPnrt5zF69jNE6XCmlMEPWsTgwPikuvJvHe+65riStdU9r/QPgGvCLwLdf9n9YKfWXlVIfKKU+2NjYeNl/3XMbzVJZVWpwclKIM7JSZE82vM/mxXcwnjG2qcY4CECXLo6v2wVOzT++sQVcFmUdbd6RAibsJmJZjfFoxtimGTebd9wYzRh7JXsC49G6fi8yNDAe0FrvAr8B/EmgppQahF3XgIf9Pz8ErgP0v68CWxP+XX9da/1DrfUPFxZe3z67Nm/VxskJw31FFpatzxhbVmJMEDpk/sXg3D1ujw1upvI8yIKhwLgDnJm/AG94RwqYzwzIixqdsbBpBx8Y2U1EOVYFIFAMHrVSVgWPYzXGFmXDbS6lGAuM52xHCni+XSkWlFK1/p8j4J8DPiUPkP+l/q/968Df7f/5V/qf6X//D/UcrZAaDoy1UnZNDR4f0xuapbAyYzzar8/PpXUlfB/0ULlBx5LgEfJzN7zrRreXlzaZrntsb8Z4dEcKqzLGvV4hANF+YF0ioTCjphx7FoXyZm3XZtOuFGM1xnO4K9jzhE0rwN9USrnkgfTf0lr/faXUJ8B/rZT6vwJ/CPyN/u//DeD/pZS6CWwDf+kVtPvFdLt0zy4iR+1btMI6y+Ds7Dy5qD0fPM+6jn6slMKuuJggyHdKcY/zt/udlxtEdmy75yUh53kenT+ozl+3+PWMZoxtqjEeLqXQSlFK7Qn66XQKpWfa863rLwvlBkpZFRjbXEphc43x6Ms95rGU4pmBsdb6p8DPTfj5bfJ649GfnwL/8pW07qqdnY1kVC3aqu0k36Zt0FdkpTyQsj1jbFtg7LqgwqGXfHSwog53wEtChvOrVgTGQ3WqWjm4JYNvul6P8+iw16MzvFWbFxCU7Ak+RncT0X5gXX9Z3C7FnowqvFmL72zKGI+WUgTl+bvp7LmSnsfIVm06sDAw7h9fFuavgrYtA2L74jso1r1nWb5lli38so92Li5K41/ykWWFQcz0rOP+H9wk293PP5yd5XXgA6FF/SWMZ4xtK6XQupA4sClwBLu3a7M5MB7eCUap+VyTYc+V9DxGAuPMpsD4OH8VdDaUMXacuXuhzJWzLWMM41sIWrVlm5/P1Ax02xm025f8jTnX6RQWvJqedTw6cbjzD77k+PP7ZEcnxa0ty6W8P8ky2N2ls38y9d8zl7pdOu2hDqPdtruUIsuK63YtGwysyhiPrCOxNTDWGnrti07F85jLaW2Dr6TnsL/fn4vuG8sYl+wJjAcZ4/79ZOPCO7B/8R1AkARodXFr2hQYB0FxcaHxW7a124XyLNODq1I1pN2Guz9e58nv3i18F6gOZ5/dYf1/+Ak3f/0ue6eGLVw+O2PrH3/GV58c58NCp1MMHIPgvH/pdmFv1/C+RevCNu827doAdtUYnzzdh88/v0hwTQuMz87gNW5ve9U6HVBD01BeOH+vg4bnW3xnrnabg5/eofzdG7iN6sSMsTXB49AexlopdCk0eoCeStm9+A7y4PEoKKHO8oedbjvLexTfPy+tMHWbwbGM8SAwrlRm16iXMRJcmT4dX6pcBLujgeHZ0x1uP8n/nC0v02gZFmhFEbE6ZvsPP+Pz+0sspWeFjLEb+hxtnbL95Rb7jw5Z/FPfnF1br0KWFV6AZF1gnNnz8pK9XsKTD++xtP4p5etNsqGsKgqckyO4t8HJ412Oa6s0X98Ot1eqfTL/r4MG2wPjUonToy5Pfu0mi99doNo5Ow+Mtevhlrz8gtvdpdOB4/oa1epsm/xC2m3o9S52pChFoJQ9Qf8QGzPGWhePy/fzhzanHxh3+sHjwanPgz/c4BuVdfjB92bT2Jfk+3lt/4CJiwv39vLa/SQhzxgPB8YjGePRczvvwtolD1z9W00rh9r7C+at5XIcwmoJ9eAMd+MJmyOJN3X3K+590Z95W1yj3jDoxE0wmlG16eUeUCylUAqjFxcmrZDttsvduz0qO1v0hiorXAd2f/wlOzt5V7nwzzZm19CXNLpV2zzuSAFvQGAchvmLBB79ZIP9ZCiO0pr49occ7LbZ3YWNhe/wdjrLxr6E/vRLd6iMAuxbeAdYmTHWGr74yQl174DaO3WCwM/fEJfv2Eb7DB58dsjuw6f4x3v4PzT1Qr0I+gdGSymGN0WYV66bz3o2GrCWdQqvYR9kjLtd2NqCw0N4773ZtfXrcuMQz6e46G5Er95icdXMocOvlvH9s0KF3cDgQVQrh/jtlvGJhdF9fk0OHCcZDvwdB7OeQEeUY4WOE9TBHvv7xe96PXj8OP+zjhPqy4aVMA0x4eUeYHtgHASEkWKQ6jg8vPhKZT26xz0eHEMvruDEZer12TTzpfXri7M3IDA2uO+bynGg3AjZ+J1bbP3RffxajNO9uDX394H9R7hAWAZKZnWMw1lTzwMVjtQY9/fffvw4/91rtcN+OnYO7e4SOy5KpWxvw/HDNt5Qwlt7Po8ewe5uPpN97drMWvpigoAwVBx2Jj9xaqVI3lsyd21GFFEu77C3N/1XetUGrWUzh8b19fyBzfMmZIyHSg2yLH8eLZdfdwuvju5lDIYD5SijBwfHgaCR0D245MIE4reaRo/rnZP538MYbF98B/hx6ZlP/p3WCs2mQQ/UR0fFVOnowrv+Vm2mZzwmGdvH2JLt2hYWFZ3WMgCd3SPcw8kdZKmEcYHx48dw9+7Fg6lfunjtda8Hh9ttPvko48mDDuVHN2Fzc4atfQbXRX35BdXdu9DroUdqjDM/YHsb9Fkbb3eThoGznkFlejlFL62zdN3gl3yUy0TPCAZ1a8HYkvfup1/y8Y8OePy4vzZhiHLUefbx49/e42z7aEatfDGPP3jI0VcboHV+zw2NgcM7UmSnbY7ub8+ghS+nvHh5MkArRf1dU7N3udGM8Ty+9Q5szxhDXk4RnXJ4MPnrXpSQlRMWTCpmPzyEe/foXH+X9b0SywfHuIy/3GP0yTLLDAr++7KzDnt7UFvwUeryN99tPzqlVHaJa/N5s10mDCG+3qSz8aiwanfS7xGatfCu0YDPfmeH3c9O8Vp1el6IG4TQybdpu38feqdPibbXid7qQml1xi2+RP8NhOnZJmfre8UXKADu4T7e3hbu8QHhjWUjM6ulSgkeTt6KrfT2stFZRqKI6JLbpxclNK6VjU0+JlWXvY++YGujwlaaMjwEnJ4pPv5HW7ibT/HbJ1R/4Wdm1s4XEkXc+607pJ+vU//+WuHeU46is3/C1sdP2Ptqh8U/fp14hk19EfFizJ5ShcVpw5x6jbRmVrp4ezsvn0v71X9SSjEv+nXG0wLjTmuZNDUs1ggCOD7Gv/kpu/vXOX50xuJi/nKPwaugAbzOCb3tNlvdKpub8P775gXGjuew81sfsl5p0vr+8thb0nSm2X1yysZPH9PZ3ONb/8rPzqSdV2FxSfFVY4lg/cHU3ymFGJcxDkOIl1M6f/AVauMRXhCOFYcHG49AQVhivm9Gz4MgICq3UVvjDzClx1+d/7n2vklP2xdK1cn//XtJleW3DH81eRAQJh6O0x19pgGgW1+g2Xz9zboq8XIKH+/gHu3DUbFY1dnexCGfjSkvp8a9ujxZSdn9QzjYOOXgN24x/OzSPu5y8+99Av2yrcpbtVk184UlqSILY9yTw4nfpzeaxj2wOQ588QW0WrC2Nv466EHG+Pg4z/ctLs6ilePsD4zDcOo4m4VlsqRqVrYYOE9D9XpUdr7isAuPHuUdgvYc/Cf3cA/3OPyqzaOFb9CNoVrFyOwVrkuyFLN1Z52nTzcKr0sG2L+1wcEn9wEIl2q4nmE9x5BKBdzlBfTWE1SvO/4Lg8DRsMAYYGHF416thb+zjtOevAtFKeiXysz78UURUXj5S0l0tUZ1wcySg2mBsVpeolZ7vW15JaKIcvmgsOYE8qRCsFQfTAoYya1XiMLz6rqp0rfMm5KPaz46jFCn4wc3/JwdL6dzuw3YZTwP/HpCNiEw1q5H44Z59T2D2aXNzXzdhXfcuajfVbC557P1OL9e52k9hmH5wxfQzxhP0mku43mY19kHFwNueagT1xqczhn+zgZOp80RMd0433/O5CxIer0GkE8xjfT46uwiyCqvmLjX3gWlYGHJodOY/NhcCkAFvnlpf/IHM2dlCX1JyuP8Pp33wLhcxvMub2by7qKJpwkAJyoNdzEA9KKY1jupcRmricrlicFvt9aitWD4AZZKxPVnPJApqLxtXmDsOBAuPHtHHhOD/oFoYXKdsbfYICqbd20GvRM8N39q6XYpvnBNw4OnHicn+bmdpxjF0K77ayiV8L3xhWhZENJLa7RaBi5mHUr9Xlbv12nltZqua2DwPyRarj5Xtrty3ezAGPLOIWsuFt58N2DiwrsBpaC5EtCrTF+NFobkN+q8L7vuR1XTFnFlQUjD2L0fgTAcu8yyhWVardk058pF0VhgrJWiW18wcrHkqHjl8sxieTHFLZk5WfysY1OGBv0D8dLkwLj67hxFjV/H3h7VJ5+jTvMtZRmdCXXz67Be03hn87MY1P7AOMjnZ0ezxp3mMihlZmev1HlwHIbAhMC+FyVkSd6JNBoGBv/DgoBk8fIVPyouU66aN302yvOgseDSrY/X94QG1hcPa7Wg2995YxJjjq//NFqeMuXuLC/O7W5zz8XzCKKLh5MsCKm9XZ3755XnNsgYD/WJvaRGteVbsZNPtJheOluRvm1u9J8sJ5fOOsXLKU5g7klMqu754vmBLAhpXDN0xWu5TFkfEd35FP/pg8LCcu35qPYZ/tMHLD39KWMbOM+Q/YGxUv29OS9+lPkBvWqDatWMcXii/lynUkxcZd1ZuFjZb2TwP+JZiylML6MYtrgIneZ42YHJGWPIn+VqyyHddHJGx5jjK5XAcSZOx2vHpfauuYHHwHCdcae5xNKyyU/WI8IQxy0mS7qNxbmayn0ZqlqhPG1LBqWovl17nc25UuXURcXT95uovmP2vVcqgVMpPlWXrzfNXB8EEMfn95m//bSw44bqdohuf0x68pQo6M1VoGJ/YAz9LdsuPnYbS+ZmiweG64xHHiZ75ZQszqdyo8jsTdwHyqu1S7M5yZo9gXEUQVLz6daKF6iJW7WNWljgfL/mYb7fr6Aw5fjKZYJgvESrW23SWDA/tTp4NbT2fJK3mkY8rzw3pfIFeP0xIStFOJWEqi1diOcRNydPZ5SXzC2jgPzUTaszdlx1vh7FZNFi8fiMDvZdl7B6eedRrzF3uwO8OYFxf7zVrke31iIIMLsjHAqMRzNXw9liW7IgRBFJc/INpl2PyrIF0f+QxcX8AW6QNfYGdfKGRyhJAmG9TC8u1goas/BuYFBnPHLvhdcXjDmEywRpCRR06ot2ZYsHhuqMO41F88vNRkyrxU2um1t/OzD12FZSlG9u0D9QeNFHmlJbNHN3m4GgVp5a2uM4+W5M87Y12BsTGA8W4HUaS+A4Zi66Gzb0dDU8OPfiClk5v7GUworFJAPTssKlhQqeb/LJHFergReX6PXLDs6DLQuirrxUpJg1Ni4wHtQZDz2P9copjVVDMt7PoKKQIHTxVxfON+e3Sv8NeNpx6VUa9iQQ+koLlfEEnOFlFAPJUjxxcbLJi+6GJXWfzM/7weTtprG725yL46ndeq0GThgwb6+aNP0/+fPpn5VSOV/UpNRclbO8mKGMseddxMnD2eJKZa5mJ15acq02sZOwqb54QKl+2UE/gDRmx4bn0GiAqqT0oovMSMnQwHh4Z4qstUjdjrEZwhB/dYHFFfOvt4miCN8DZ7FFFDtWlJsVJAlxUkwWRIvp3L5p7OsoxwpVKT6tuZ4iuVabTYOuWBSBShO0cqi/U5t1c15euTy1Qq5eZ+6yxfCGBcbBtUVw3XkrZ3kxIxuNlsv5m6my6GJhgm1ZEJUmJLXxjr1ybb6eNq9KqwVEEd20Zs7CtOcw2LNyuNY4DMmDflO2BQjDfLebUv4Qk/kB6TWLdm4olSjfWLRqxqmgP80WXjf7TXdTOQ7xYnGRmsn7+44qLxUDY1vKKCDvT0rNBLdZI6la0KFMCYzjGIKSmstA5c0JjB2H8K38xQlz+IDy9Y0ExlEZ2kPZYtc1vIZ6EqVIrxUPSqUJ5YodHeIo38+fqLvNZXO2MntOCwuQJVWyUpTPeJhWP+04UCqhVP5Q2q0t0GxZVM7jODSXffOncafxPFhYIGmWrA3+45XK+ZZ02pIyioHycjEZUrlh10ksLyZU3pm/gPGFuO75Yt5h9Tp5LcUcZilt7faKHAdWV4mrHqXS3JWzvBjPKxRJR8s1dHgxH9hoGPmCtGdK1qqF2vDyStXsWvFnWFwEXY4JmqlZgeMzhCGkaV4qcn5YpuxIMXBeTqFQCy3ranFtvq8AWFuj0ZjLcflKDF4PDRAt2FFGMZAuRuj+yyE8354yioGkZfhLgkaU6uXCvuGeny/EnteaVgtDpykWF/F9WF199q8aYeglHwDhOyuFadw5nJ24Ek69Wqidi1dtS4sXxXHed6jVFasCY8iD/l6lTqnSPy7Tjm8QGK81aCx69geStnFdK5MH58pl4ko+KJj8Uo9JwjDfrxkgWa1YsfZiWJpCKbSnQ1FJTDjUvddqoML5zVLa3C0U9Uctq6bNBuUU9TqUy+dv2wrDPKCykuPkHSH5HquVpSmvH7PI2hp5T2lZbUy1mteYld7u1xqbFhj361TLby/Ma+JDvMmUIl5OrSujGBjUGduyG4XVhuuMVR4Yz3On+eYExjYaZIxXVoCLYHiOr7crkV6vgQK/VbV2GnTYeTLElIVpz2mwO0x0rZlfy6YFxuUyxDFOGhtXBSLeDNFiSrxcwQ/tyqhCnin2fEW8Vpt1U8SzDAXGlRR8f763BpPA2GRBkKfA+5mrQcbYqqz4BF6zSjmy6213b6qFhf6U4dKSeYGx51lUmyVspKoVlr9jZ0Y1bQYk7y5aV0ZhJcehVMvjlPNFd3Oc6JnflolnC4LCU1ccz+0iz6vl+yTLCaW1+axPEs/vvG9cWDBzteic1sgJAUAYUlo2+81p0wQBtP7YyqybIZ5T2CgThif50ow53xrMwJFInGs2Cyv5Hadfj/oGqHz3GnEql681TAyKhTCBxfeWjSUitnKSMouL5DODc76Fj2SMTTZhCulNqXX0a7auLhRCCCEsE8f5Oqg5zxaDZIyFEEIIIcSrFEUXrz2dcxIYCyGEEEKIV8dx8h205njR3YAExkIIIYQQ4tVaWpp1C56LBMZCCCGEEOLVMuT1oBIYCyGEEEIIgQTGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBADerBswCyfbJ2x9ug5aA1BertD4RmPGrbo6D/7xXdAajcIruaz8wrVZN+nK7D865OjxPt2jM7rHbRZ+sEaynMy6WVeic9xh58tN2odtuicdwmbM8s+tzLpZV+bR794n62agNcp1WPuTb826SVdm64stTtYPzj+3vr9MWAtn2KKrc/T0kK2PHuOUfJzAI16pUL1emXWzrsy9X/8S5bk4JR+vHLD0M0uzbtKVefLj+5xuH+OlEUElpP5eAz+yY9g/enrIxk8e4UU+ThiQrFaorKWzbtaVefiPv0JnGhwHJ/BY/YW1WTfpymx8usnJ+gGO5+D6Do1vLVCqzk9/accd8jV1jtoc3Nk8/6w8FywKjA/vbuU3FOCGvlWB8fHTA3Y/fXz+ub13ApYExr2zLlsfPjr/PDiHtjj4aous0wPA8V2wKDA+2Tjk4Kut88/VdxrWBMadwzOOHu+ff1auY01grHtZ4djcOLQqMD7ZOOR0+xg2DgGovl2bbYOu0NneaeFhVDnKqsD48P5OnkgAlG9XYHy6dcTR/e3zz+n12lwFxlJKgX0BiLbrcAq8clD43Dlqz6glVy9IisfWPenMqCWviFLnf7TtGh09nqFDNd7odeiW7Mmn9E7tPTa0pnt4evFROWN9jMm6x8W+f3RsMF6Wnf9RuXaFar1OVvjs+vN1fPPVmtdkdNCybZAusGmEBry4VPjcPT6bUUuunuO7hQ6ie9y26tosXIo2HZjlsna38Nkp+TNqydXrndl7bNlpm277IgDx0siq4cDqwFjrYhfp2BWqDTLhAxIYzwHlFHsH+8Zo6w7o3GjGo3dsT8YYwIuGBuZej85ZNv2XTWNzxhiLIo4R2Vkxq+oE9mRVRzPGKrAnMG7vnRQ+e2k0o5a8Gr2TYt/vxxYFxllWGMVtyxhLYDyHxjLGs2nG62FTigDwyz566Ji6R/ZkjGG8c28f2hP4D1+K2rbIeIRNt91oxtgN7QkeRwNjL7Qn6B8NjP10fmo4r8LoubMtMB4OTLSyK1QbLaVwvPk6vvlqzWsyOmgpywdpm3i+Av+iA+yddQu1WKYb7dzPDi2rM7aUzT3IaMbYqsDY4lKK9v5p4bNfsTdjrB0Xv2RPOKN7xTHN5oyx48zf8c1Xa14Tq2uMrTqYcUqBG10Ej90u6DN7sqp+uTgw27S4cPjGU6M1dJaxNWOslcIN3Bm25mqNlYlYFBh3D4oZ46BqUWDc7Rbqp1UpsKsMdyTZM1r+abrhwN9xmbsa6vlqzWtifY1xoTjJrhsKwC0XF+B1Du0ppxitoe4c25MxtvBSvDDah1jUqRRKKVwPz55qA3t3pdCazsHQjhSOSymxJ+jPTtuF2HE4WWKD0YzxvAWOLyPLQA+dPNdVczc42PNf+yVYNIa9EWzesu1NyRiD3ffdnPXzL67bJetdnCjt+bj2JIztrZ8+O6PTvjhvWSmiVLrk9w0z2i/aHhjPW6nBy+j1QA0FxvNWXwxvaGA8Vkph0z7GNm+o2ucnI1u2WbQAz+a9jK2+72zV6dDrDX32PJuSV2OlFNYsvjs5oT0cO4ahVZn+0Zk0q7Zqw+6Mca9HoVREAuM5YXWN8QgL42KrM8ZOyS8MYFbtZWzjxdg3ul2bNYfa7RYCY5u2M4NiKYV23LnbNupFdQ9PC/2GbQvvrN7DmDcgY6yHa4znr7O057/212BbIXuBNVHUdKMZ455FL/nA9/OdNwY6HTr2JI0LJGM8/3S7U+hSbNrDmCwrrI63qUxkbKs2ywLj0T2MC/u/W2C0b7Rpu7Z536oN3tTA2OK4eIyFBzu2l7FlL/kYrjNW3Q7ts+ENLc0NJm1+IDX4tFzK5u3M6HQYTsxpixYWdvZHdqSo2LWHsWSMzSWB8ZyyudbR1gF6mO+D9i+yxt3TLufzvVkG6+szatnVmLqX8dER2e2v6HYn/CUD2HzfjVKW7Gw8umuD6Rnj7ue34KQfNHa79IbuJWsyxlrTPriYRdOuZ9WOFGD5W+9gfLs2CYxfq/lr0SsyHExcVmOs9UW/aSRbI+Ozs/POIghA+8W9jGm34fAQPvmEzvbBjBr5gvb24PFjBqtlRqcFOwencP8+R7//GV9+afD7TCzelcKaY9nfZ+vzzfPjGd21wfSM8eZ6xuN/+Cndrx7A2VnhXlKBnz/Q7O/T+/I2B48M60cGTk+LO1IEoVU7UqB1sTbc8/EDu2ajrF5815nv10EDmP34/zXcv5+Py8vLEEwZoLe34dEjWFqCyJKSrNFFQabSnS7qk0+gXsep11GlAI7y7zod4P59sr0D1teht1Bnbaat/Xp6UcKj//E2afyIdDWlpIodR/fzWzx2MnZ3oLtcJjA0OTJW1WNNNDnO2AqmUomDj75k5+YWS7/wFmp0n9/QB605fbrH1q1dku+9TbVmzsGGzZjHd/fZ/+2nLCyqQilF0D7k7Pc/Yudpm92TEo1/OiWdXVNf3OlpYV1CFkbG9hkTdbt0O0OBvx/gm/28BmdnPN1waK3ksxaXlVK029DtaMqxOffdsOG6fpjPjPGbERhrTc074vZ6wtYW1Ec6iW4XPvmknynudqh4HaA8i5a+kKOHu8Sr1Xw0Hgk2hgdo3e1xuntK1Ipfcwtf3pODmO5mRG1/i2Rzi/AQBpNpWsP+wwOerkO3A4vfNOupxg1cnIUmj29t8PjJAdFIOeDB3kVHUqqZdWzDRh/SbIqLbXkApVSinLpsrR9y/9c+pRQWB63uo3W++mKPk/0O2eIyKxWzjjteyPv1LIOnT4oXYO/gmNuf5X/uLLZoNF53617QyUleXzYokD45oT0UGGvT9zDWOv9nkDVtt4vlZJ5vfm2449D+g4/4NFyi9b0l0u74m+9OTzTrX+xycHuDG3/2PcDMuh8TMsbz16JXQSkqBw8J732Bc3TA3l7x615Xc3LQwX9yn9q9D43rRHa/2uX2r3zE3q3NifsY626P7Y8fc+u//ZDDbTMXqjWbsOUt8eA+3LxJcW9V4OHDPCgGCGvmLTRpfWch/8MzSnmCmjkPbKPeqBpjs+LFgqiZX2NKa9onxRut/WiTk/38RkvfaRk3w+vX4mdmF7VSlK83zQm2ul2OfvwJHOSlH/r45LwvBPCS0OjrEeDwx5+i9/YB0GdtukOXpRUv9/B9kkaAu/6Y7X/0Efc+2i98ffxol1u/8jFHH97GDT3iiplBMZhRY2zKrf/S3EaVinrI4b0DsqAYOLknh0Q3P0LpjHjRh7JZwUdyrcr+nS0e/c5dtkYWWWRHJ9z6bz+kc9pDK0X1reqMWvlyggDS6zVO1wO6ncuD+6huXmBcqkWUl1KOn06va9SeT5Sae8vatCtFe/eYoHucP7FdFnV0u3QfreOtLWHKyq58Runy+tpeXKGxYlgGAcD3CVOfzvb0PRB7SY3WikFz80HA1pMO2+tfsPyzy2STtmrTGo6O4OgIvbhkVqCsFLv7Dhv/w5esfL/VL+e5+NqNgvz49vboPNmivbBG3DRvDIiXU/jqFNXrwkmxtl/tbJ/PSSXX6q+/cS/h9KhH6XQP1cynYEZLKcYyxr3ezPvK+QvVX5VajaRfMOa0TwtfqW7nfMPpeNW8wDFdrZwHHec7GPRlnR6d0/zx2m9UCEJzT/nCoqJbX7z0d/ykNJdPoM+j8e3Ljy0Ly9bUvoPZpRT77ZA7/+M9Dn70MWxujmW/Va9L+85Dnvzahzz6bH/mHf3X4VXjZ86aOYstkuT1tOeqDTLiU7VaVE0aBoKAIIDDA7jzO0/YeTK0I4VSlLcf0PvDn7L7u59z+zZjM6YmCOsRpydw54NNnv7R48J3wdEO+7/1U+7/5i2++KRLzzcvKAYImukza8G1coxLbp20XW79T484+PFncHg4OWOcZZw+2GTjNz9mHjbuNzf99HWFIWmzxJPH018GoZSZgbHju0QLyaXZRsgzriarVMBbbqE3H6Oy3sTfMbkGN71WJUgC2oeTM+JZKSI0s88H7CqlqNYdnrgVTm/tET66ix7ZQu/p791l/wDQsPgDw/qUcpmonG8EM4l2Perv1F5rk65SeSGGLydHh5lfovZWxbiMapAEsNOm14OdnaGvtObkwRZf9EA7Htn3FnjHsMsRIGr0+3V9XjFy7mTziIeb+Z/VO4tUKq+3bVcmTUmSfBOAaVS9RlIxK/GTJPDAS3nw+Sbxg8/JsmL0f/J0nycfbXK836W+GjEPg5xZ/4VfkteqXVolUY4VTtXIdcjPFdBXrhvYI45YWHbp1lpTvzexvvicUtS/uTD969isHSkOD/N/BqaVUrTbeY24SXwfgoX8fjo9oVDTCbC/z/l0b+Wt2mtt20sLAsqV6aUE3VqLRtOkyLEobJSn1kZ3ay2azdfbnqtQSqd3DN0uoKHTXGZx2TEr6O97ZpaffHeK5nu1V9+YV8XzzheHThNfqxt3/nwf3Hr+tHJ0NL4H9eNbxxzv55mF6jvzseL1mYGxUuq6Uuo3lFKfKKU+Vkr92/2f//tKqYdKqT/q//Pnh/7Ov6eUuqmU+lwp9c+/ygP4Wmo10kvi3ng5NWrKc1jl2uWPyU4aU64aVDc3RbMJWWt6yUFYNzdjDFB/v1V8JfSQoGrWsXkefP45fPHFeJYHoNPW3L8PH30EZ6fmZY+TtWc/aIaVwMjX8V4WiITXWsYtUB6m4sklSVop/JWWaUtMAChVLj8h2vXQrQVa03MKc82Jn73lXK++QLNlWNQ4oryUTg18teMaV0YxEC89u+6qVIJodT7qp58nY9wF/l2t9XeBPwH8FaXUd/vf/VWt9Q/6//wqQP+7vwR8D/hzwH+qlJqPaDOOSevTq0dMLKMYCKrRpVmD8mrt9TXmFXJdaCwHdNPJN5CJC++GKd+jemP82LRyCKtmRSPh0RbJ3kMOdrp88QXsHxR7/Fu3YP2pxtnZpL5ze0atfHHVhYAsvDyKKq+Y2acE9Xjirgy9cmrmorthvk9UHe8re0mN5pKZ1YVeObg0p9NpLLGw5Bi3i8g57/K392nlkL7TMmcnkSmcako8ZTdVXa2RGrY94kBS98lKlycIKisx8/LE/czbRGv9WGv9B/0/HwCfwqXvT/hl4L/WWp9pre8AN4FfvIrGvjSl8FvViSUsvg/hkpmD2MBlgb0NZRQDCwvQbS6N/dzz862JTNf49uJY1iALI6KyYZ1imlI/e0J080P8pw/QI3vseXtbhLc+pvT4LknLvPNWLgO12qW/Y+x9V56cVc2aC9TnI6nzUiZlxLv1BXP2Lh7VX4A3iXY9eo0FFi9f2zv3LpsN7FUbLKwYHhUDpClxMrmfL1+rG/tgkyT5Q/VlJiWEZuVr/WdWSt0Afg743f6P/i2l1E+VUv+5UmpwVGvA/aG/9oAJgbRS6i8rpT5QSn2wsbHx9Vv+oqrVieUUcTOcm6eVFzV1ardUIl00bzp3mijKF9D0ouKjdSktWfHqTK8aky4XB+4sLM/DmoSvJwhIFssoneFvP8U9OSp87W8+xumcEYbgt8wLIJ+1WNcvOUSLZq5ZII7HSgq065Fer5labVZQbpUZfidL5pdIV1Nz36BWKk0dvjqNJRoLrrnH1ndZYOytLhq7S0qB45AsjaeMtetRfUa55Dwrlbh0/VYcg79oYGCslEqAvw38O1rrfeCvAe8BPwAeA//R1/kf1lr/da31D7XWP1xYmL7g6MpVKhOnI0yd8hxWXkonvkUmWqnZEC8WLC5Ct1HMGpu8I8Wo5sjWbboUGblVW3m19szpzbjiMnX+cM5Vlstob3LEES+n5j6oed5YuUG32jS+hnPASWPCoUCyW18wctHduSCYGBhrx6VbX2BpfILNOOc7U4zolVNa1w3sHKcImuMPaFmlZtSr1ycpL00PjNOVhHlaWf5cvbZSyicPiv8rrfXfAdBaP9Va97TWGfCfcVEu8RC4PvTXr/V/Nh9cl1JrZL9A9XwLaeadcp2JF19yrfb6G/OK1WrgNGpk/sWJNHpHihHhSp1y5SKiVHHZzIzP0P7h05y/ztxAlQp009rE79JrZvcpYTMuxPVqoXXp4mWjDJWKaKXQjaZZexePmlJK0WkuUam7Rj5Uj3LTybvyZK1Fc0tgJqlUxrLf0VrD+JmapOqOzfJC3vXPUxkFPN+uFAr4G8CnWuv/eOjnK0O/9heBj/p//hXgLymlSkqpd4D3gd+7uiZfgZHdKcqJi1u1YR5mPMDXrkd11cxs3GWUgtZC8YUfpu9IUeA4NL55sYTc2Gx4FJE2pmcCHMfs2RrPg2jC2gTXNXsxLxR3b+iVU+or5r9a+JznEdXy67KX1mksesYm9wFQamzxdZ4tXrQiWwxAGBJGxQsw8wOqb1WNDxoL4pg4vbgYtedTXTM/PklTyCbUGVcq4DQNC4yBfxL4V4E/M7I12/9DKfWhUuqnwD8L/O8BtNYfA38L+AT474C/orWe/DaGWalWSYfKdcrLFWMzVqPSa9VC7VxpsYof2HFsoxYWoFdvoR0Xx4FS1Z6MMUD67gJBAFkQEsXmjtrx2vRSnjgGVTW3dg4gXU3RqniA8UIZVZqfqcEXEsdE/TpjU/f3vUzUzBMGthybHweFYazbWCRKXHNfeDFKqbE+vltbYGHRsvFNqXx7s/5hddMatbr5xxiGMGnKqbKWMm/Toc9cxqm1/i0Koda5X73k7/yHwH/4Eu16tYKAqFHGe3BMt2NHGcWAGwVEtZCTnfy11zaWUQz4PtSaLofrTeKz9bl4Y86VCgLq79R4+EgZfWiqXiNJ1vOXXoyIl2JM32OpWnfYSCp4B7vnP7OiTymXKUf5rFOwVLdiOn5YUCvjxicEzelbZBmlX2d8eppnizuNJdZsyRb3hfUI7p0A+RZtpTUz951+FqeaEpf3OTqCcKVhehcJ5LnH8mJC+65C6Xzfes+D+Pr81cGYm4Z6WbUaadKfCl22YBAbMpjC1Uo988UfpltYyDMjQSU0d6HTJWrfXMRJJm+dZYwkIa1NnutMDC83gHyXFKdeO/+slPn1xQC4LlE9pFezZ9FdQRxTWrMjWwxAqXReg9ttLOKHrl21txQX4FmzRdsklQpxnJdR1K6ZX0YxkFQcsvDiKbRSVaihvnNe2BdJPK/+tm3lBfMzVqMGe6d69YrRU/DPI00hrJYo3Vh59i8byKmm1N+pGZ0xRqn+Arvij03dpm2S4QxxueLhpDakIPPdG0prLesCLCDPiF9v2nNs/YyxVg6d+iKLi9ZUCJ5zk+h81j1rLVqxp/ZEUURSdelWGs/aKt0oSQK95CJZV71emcv4a/5a9LqUy5RrATqyY2AeVmomBJFLuFabdVNei4UFiBJbRrdxi2+F89h3fC1us0YSb3N4ePGzuOZjyzxotelxECW4J4d2lFEMLC6y0gznrQTwarguzSWLJpr6exl3G4s4gcfr3AX1tYkiwhBO/ZTGWmTPuRulFKVWSlKuz9MuZi8tjkHHKWzkexuHK/P5ZGPrZfVcVK1q1yA2oBTl5Yq5b936mhoN+8qLh5keFANQrZKkxfRVvGJPmU+aQpbm95sVZRQD5bI9W7RNYFVgFQQEoUOnsUSrhV07NQwEAVHi0m0s2hn4D2u1uPYtO2aeBpSCqBWjlZPvyzyn6XCbuoWvb2HBmozVqOq3V0jqNqZ5xrmufVOG1nEcktWLQNh1zd6mbZTr5i8zicoKv2lPwC8MEgQE1/JssTVbtE1QWqpRXqlanQwBoFq1Kls8kKSKrJzkibs5fXqzIRf14oxe0XS5ctPeYxNm8hdqRNEeJycQJ8r4bdpGpQshzpmtqTox95RCrSyz7M7VS8SuXPSNNRZOJBNiqjSFs6UUf2nKO8znwJsdGAshXp/+gteTk3zbHtsCyGoViFdn3QzxJnNdlpdn3YhXyw196rZniy0Wx9B9uwrV+Q2M3+xSCiHE6+P7pMt5zZxVdbh9YQhhIrkGMVtSVibmmetCYy2a6wL/+W2ZEMI6wWKNahW8hl1lFEIIIewggbEQ4vWp1Vi8Flhd3y+EEMJcEhgLIV6fMMSz7T21QgghrCGBsRDi9bJ+A1IhhBCmksBYCPF6yeogIYQQc0oCYyGEEEIIIZDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAIAb9YNmIWDRwfsfrmBn4YEaYl4KaFUKc26WVfm6U+egNYAOJ7DwveXZtyiq7P34ICznWN0t4fuZdTfb1GqhrNu1pU43T1l9/Y2ylEoBWEzpnKtMutmXZlHP7qH8ly80MONAhrv1WfdpCuzd3+f0+3j/L7Tmtp7TWv6lM5xh92vdsk6PbJ2l7AeUf9Gc9bNujL3fvM2XuTjJyVK1dCqe27jix3Odo4JQge/5JC+VccrB7Nu1pU43jph+/MNvNDDCz3KCzHlhXjWzboyT//oMTrTKNfJx/HvLc66SVdm7+4uJ1vH+VjnKKo36nPVX76RgfHJ5hGH93fOP7d+cM2q4HHnk8foXgaAW/KsOraDezsc3N44/1xeSq0JjM/2z9j55PH55+r7i9YM0rqXsXfz4rx5SWhVYHxwf7dwXUYL9jxstw/bbP7BvfPPyfW6NYFx97TL0YOLscCvRFSufXeGLbpahw/3OH24df45bJStCYxPd04K91ztW0tWBca7nz8l6/QAUL5nVWB8+HCP/dub55/DRnmu+ss3spSifXBW+OyndgRWA/1ksZWU5xZ/0OvNpiGvwsiJU2pG7XgFemfdwmcnsOyZfPRkWXQTjp6rwWBtg7P94ljgxfMzOF+F9mlW+ByE9gz53dNin+JF/oxa8moMklsAyrFoMAB0NjLWzdnx2XOXfA2dg9PC57BqV2c4PEjbFFwB4Ixcslk2+fcMNBZLWXTybA+Mxzp2mwJjv/gwqjvdKb9pns5hMTAOUrvGgm77on90XXB9e4b8zknxOgzKdvUpergPGR33TDeaBHLn6/jmqzWvyXBnqJWaqxT+VShccxYN0DCeMdZde7JXo0/RNumddgqfndCu7I7VGeORwNimjPHY7GFiz1jQ6RSzjr6PVQFW96TYp1iVMc4yGI6LXXuSJFC8LkEyxjOXdXqFQdqLgrk7KWK60SfL0RvMJjZdl6MZY9eyjLHNgbHrKbS6uO9sehgdDYxtyhi32xSuQ+sC41N7M8ZjSRKLzhtMKKWQjPFsjdaUBRW76osBu0sp3DcnY2zTubO9lEJTPFkKewJjxwHtXZwvmzLGo6UUNs0enp2B0pIxNtG8Z1RfWjbfx2fPXfKczvaK9cU2TZ1NZFHmCiY8WdpUYzznCxJeRnY2UkpRsmcQA8aeYmwri3GGSph0pq2577pHxbK6ILFjxwboZ4yzkcDYkqdtrYsP256X79xgi3nPqL6ssf5xzq5Lu/5rP4c3bUeKObveXprVNcZv0OI7L7RnEAOsLqWA4n2XZUDX/AV4upcVso5uGFhVy3l2BgxnjAOsyRh3u+RF1H1e6FnVX2bdkQdPi44N5j/wn6/WvAajO1LYNHV2zrKbqMCxuMZYtmsz1ti5siwwHl6Ap7UdD6Q2L7yDPGOshjPGgT3Dffs0K5SJuCW7+pN5DxxfmpRSzJfRztCWl0MMWDYej7F5H+OxbP+cdRYvY7SUwpVdKYwyOk1tQ53x6HoTGwPjwuK7kj3Dfed4pL7YshmosRpjywLjeS8btOu/9nMYzhhrpSgllg3QI2xaBAR270phW13qsKw9uiuFO+U3DWV7YOzZt2Wb7Rnj4VIK1813F7HF6B7Gftmucdz2xXfznhGfr9a8Yr2zLr32RYfuxyXrLjjruW9OxtimWorhwFi7Hp5vz7EB9i++G93LePhBR5uxGO90+7hQ0jO88A7s2qqt08mvQdXvVGzbkcL6jPGcB44vbai/0ErN3VBn2X/tIq3h8PDi89jUmWUL76Bf/zdvV9kVsjljbG2NcZbRG1pMol1v7PnGeLZljDsdODk5/zgaGJ8nGA4P4bPPjAiMO8cdbv/9Tzh4sAe8AXsYD50Tz7LA+E16HTRg0WCQKwb+8xcY2/WYNUIpuH0bmk1YXYX2fnHhnU0d4TTWlVIo0I6LyvKB2YZFQNNYM5vR6ZANnSbt+RIYzzvPg08+gSiChYXx3WCOT+DOHdjeRgclcD3m/WoNayHdkw4PfvMm1fdaExdin20fsXtnh7PDDm/96Xdm1NKvr9PJd2qIovxzXkZxcQ0GtgXGI3sY+5FdoUzWsytjfLjXI0rci35/5HXXw93n2RmUZhya2XU1jTo8JDo+40mnyd4epBO2auv1YG8PdnZgbQ1C+5LIVlGKvIMfRFoGZKqe17wvSHhh3W7xNHmeTWN0zrLA+PhEkQUNkoN1ODggeKI4Hvpeb2xCAgcHsN5LeM+AS9WP8+3Ysp5m79bm2Pf3//tPzmcVF//ke6+7eS/FdeHDD2FlBZaX+ztSTHm5h9b5Pybfg7bVGGtd7EJsqzFu33nIvXWH5veWWVjxisfXL6U4PNA8/XwXX7d56xeWZtdYbA+MfZ9k83NOe3ucLL9F78lp4YC3Dks8+GkeW1WiDmHJg7nPe1wuH4/NPobLnGeMyTMG5y8bGPTyoz2MQax98123y3Bi37atlQDrAuMggA+3m9Qer1Ovj888nZzA1hYcH0P4rWRGrfyalMJPQ852TyZ+PQiKnXJI/Z3aa2zYy3Oc/Jxt/OQRO+sLeaA48nKPo2PF1r38YeY738ry6NmQTNDRXpe4etFv9E4756OcUhb0KZ0OX37eY/GtkGqV8bK6fsY4y2BzE8phRlIx58kmWU5QH91ha3OTreVlnE7x+G79aIOze09xOme0/qnZP5QafjU9QxAQlBTe0x3ck8Oxr4/bLu7uOv7BLis3APXN197EV87w4CrLYOPOIa1rIW7Ju8gY92md/1I3c1j/fIcqe8TfuzGr5r4UWxffZWcdhuMq6956B9YFxp4H6VKZ44chp49Px77fHEq4xkuGBMZAUJkeGA9Uv7lkZDY1DOGkfQYffcRJcwlVTs+/e/QI2nsO7bMOq94Gzkfr8N57xgTGW7d22djdZ/UXr+GVA7J2l8GsvO8P/k+ufdQhiM3qY5TvEd39mHt36wQ3VqmOzIRmODx+pNm6uYOz8ZRv/YvfwKQlYkGrQhBAu92DRw8Lj9mq26Fz6x4O+QNAulaZVTPP2R0YK4UfB8AZqtsZ+zq68ymQ9w3Je994zY17jQzOojoO7O/02P29n1J7r4m3uoh2im/hevLFPjufr6OOj1j+0yszbO1LsnTxne0v9wCsC4wB6nV4UG0SbDyc+jva9UgXzAiuIN+3/uCS75Xv0fpm47W15ypFERzGKd7+Nv7mYzz19Py7LAP3+JDyzQ9pvafzQDIx54EmbMZsfXKX239vj8q3VwvjuefBWddl+/Yu+7c2WFxxCf74uzNs7QtQivJCzMGXm/R+ss1WKSyEvfv391C3nuJ22pQXE/P2gfc84oUy7YfHl/5aebWG480+4J99C16x53mBR2MtIp+/MJ/WFAZpG2KryrUKPVy2v9hk/Tc/wRnK/u/twc4f3IGjI8plcCrmdPZnZ8U36172go9ez9xyagmMzVSrQVZvXvo7Kk3OF3yZ4FnjQfLuIp6hb4iLIsiGssTDNcYAqtelkuo8uVqrGfXkHdVDtHLQvYy9jx/gtC9mMU5O4Nbf+4Td37+FOtin+l5rhi19cclqnilVOsM5LQaQ7tE+TqcNQPpW/bW37SrEK8/OBFfeno9jM7MH+BqCyuUdoe9D5ZvLr6k1r8bj2yecnT57ID48hN2HR6+hRVer3lB0Kxc3jJoSdCSpMioLok9O+eR39nj8uL8d84SMcbcLDx/CF79/gJN1J/+L5pz1b70Do4KM5+W6UG359OLpA1q0YM79BpcHxlo5LHx34TW25mqFIeigROZPX9LfaIz+wQzlWJGF5Ynfac15Brm6EOBU04m/N+8G5QaXUvMTPH5d8Url8m7SmY8yCngDAmMVlobLj8Y0lgNUw8wLbSDbP+TW3/uEe3+4NTFAPtjX3Pxgl7v/3af4e+OrseddEECw9OyOPF6MjVpqHVYCoke32Prtz/jkR/uFPbcB1jcUH//2Hjs/+pza7lf5nKEBsp7mdKgs1caM8fGR5smdk/Ms/tiq8f5DTrsN9748o7dv3gMp5PFTtzo9axwvxq+xNS8vqIRTB+fyWy1KsbnX5qBcOIsnB4ZR1N/OzbAyCsgf0rzK5MB4WP39lrkPqVFUWGA48VdaibF7NjuVhHIyfXyOlqu4/nyM3/PRilcpDKfuied5UPv2srk3Ul96rYo6O+Xo46+4+6sfF+qvTk7g3j/4jM5nt/A7x5RXa7Nr6EuoriVk/vTHac+HUtOszh7HIV5OcU+O8O58Se/xeuHrs5v3Ce7dxD05JLlWm00bX4Czuc6tH21w/56m1xsPjAcryLMMNjf08HskjBGVFdt/eJfPf/UWT28f0cuKfUinrbl3s80Xv3aX7MOPjX0db7UK1GqFuv4BrRTJklmBsXKd/rqTca3vLr7m1lwtx8n3f+2VJwfG50niupmJoFLj8mstSSBYNbOMYqC8fHnG1NQyCgCUunSh7jxlwu0PjEvTM8b1BQ9n4fIaOhOkzQAd5U/TTufs/OUXA85ZHnkkqUJVzJxmqtehV5meNU5iIDXv2OLVi9r20fOmenlA6TgYFRgThlT27rH/e5/xyY+P2N+6eFDTStFTHg++POHTf3Cf9Q/umbIwvkApiNdqsLvL9m9/xv5Pvyp8v37rgMPf/Qh3Z5Ok5kH52dmueaQU1JsOvcr4oOUkMVHZvIB/UnldsFQnaZr/wqcoYmLpi+cPdY+GlVEMhM3LA+PajRqXTg8bIFlJL10YVL0xP8Hji5haZ+w4VK7Pzzov+wPjIKAUjR+mUlD/1qJRU+/TKMVzZYKTldTY4w1D8BYvCYwNqy8eSFYrz1whmVRdVGrQsZXLRCE4p8d4Nz/jaLdYY7z/4885+N1PcDbXqbxdN3bCJrl+MUgNHj4HVK97Xgs/WFRjqnp9cjmFafXFA5MC48Z3ZvtCgasSRYDnkZWKKyIb9f7EaBBAbFaWf6BcL02cuYA8Hk7fMTtbDODUKpSnLGYNDS6jGAgXK/nrySf9fI4Wvc5PS16hSa9+rjcd3BWzp86GJWvPeNpSz/E7c666HI11+AAoKLfKRgb9bhwSVS/PVKXXa2aV+/g+5drk6Wql9fme4pkf0LhhbtBYXSxNvh6HeH7e6ZssTcGtJmOLusqLZgbGowvw3GpCbc3MYHHUYPZluJxCqXwTCsDYMgro77oRTT5P9aUAKmbfZwAEAXFz8hTa8IO4saKIpDYeGc9biYh5kcQL8JORwENB41sLXLy423yV5fKlNbhJDE7d7MC4XofuhHKKKAS3Zl4ZxcBwOcUoZegDTdSKn5kJ95dbplYYAPkahdJS7dLfSRKMH7CV6t97tWLW2LT64oHRwLj+nWWjnjsvM9g6r5dcXHO12tBQZ2gZBeTJbhWPdxhKQe0bBi+6GxEvTx7Lau/MV/D4okbLKbRSVN+arzHujQiMR1ciV6sKf82ebDHknYbfnH5xJa2QqasQDVEug7sw3rEnCUaWUQxcNtUeJ8rIBxoV5+UUl6m+a0F9/zMyHclibMxuIpdpNIo1/ioKiVIzjyushecPbSoKab5r3v01TRjm8WEWJej+oHceC5dKxta6D0xagFetgrtofl8yEC1VxnJ2pWaSv+bbAvFKsXywtFDFD+crSflGBMaEYWF/wOb7DZ69YaB5LlugZWLWcZLqQkAvKgbBseGBcbSYTq2vSlYrRpaIEMeXjsG9pEpj2fx7sLocTd03Vqnn29TeBHGcl6QN7j1TyygAVOAT9AdiU1//PI1S/fyH65KFZZJ0aKgzuIxiYNICvPqNql3jeZqOlYFbUUbR59bSQtJk3soo4E0JjEul8/smSaH0ttkv9JimspZOXJwQRuC37AiM63XoVS8yV64LUTM2uyzGcSgvTZg+U/36YhOVy0SXBMbR9ZbpC8iBPEPnNmsTv4vjfDGNLRqNi3IKkwNjAD8NwfdpfcueTOPAoJwiiys0hyfYLAiMy1UfPbR6K4ogvGb+orsC16U8sj+46btRFPg+cSu/SOexjALelMB4aC/jxjs1jNwf6jnEiUJVxwfitOoanVEdFsdAvX4+TRjHWHFsk8opyhF4rdrrb8xVcF2i+uSFadrz8yyPJaZlPJKqa+wOAJPU69BL62jlkCybfc8FlZDk3UX8wI661GGDwLjUSi9mbcLQ+DIK6G9HN7QAr77o9zfbtksyNNMUNBJKiQVZhCGD8c5vVCiV5y+p9WYExp6HH7pEEcTv2ZkthunbtqVrFWsWJigFtZZHL8k7wyTByP2LRyVr1bHFaslyYnR9qlcpT5zh1I0mtbod1yNAdTUuZLEG8q347DnOKIIocVGLC4TP2Ell3pXqZVrfMff1z5cZ5H0W3kkuyrAsyBZD/9j6U1GeB5V37Vl0N8xvpOfJvHksNXhZgzrqeT22NyMwJq+Pa95IrcrgTJJeq55nUyHf37G0aNcT9fDLPmzJGLvlEuHIanljyygGptQZJzdaVtV1xjFD+2HlwhD8pj1lFAONBpRurMy6GS+t+l6LKJm/TNVViKK83683hvZ2tyQwVupiAV6tBmrBsjKKgSQhTvNOcp7eCHdlkoS44lJ9uzbrlkxk0fB0ubAWkr5vb7Z4oFJ3yYYWp6Up1k01DY6pFDl4lbLZ9cVDhhdqhREEi7XZNeYqxPH5tO5AL67QXDU72zhKKSivFQcvG7Zpm6TRgLRm/v3mePYOfaUSLC31E6lp2k/1X77ftklK9TIoqF6v2LXobphSlJdS/HpCmNpVRgGA41D7zsrcPpza2zuMcJp1KweqUZ4H4XLt/HOyZMd2UcOUgmrdyYMRC8ooBoZ3DklakfHb6xFFlOPiNKez2LIhwT+mupag3Yv7LG6GVg7aQQBN+9arWUUpWBzsRlqpWJMtHihXPJJGiWDNzlKYgXg5tTNb3Be/M79b5r4xgfHoVKfNBtu2uS6UV+zKFg/U65C81bCijGKgvJTi+vktaXwZBYBSBLXyeUJfux71d2ozbdKrUqkqsjS/12y+78DM3QPfNOfVdFFk3ZNMudzfos2ymdBRTq1C6317A+N5rg2XLs5C1YWArBSRJKBqdnYelUr/DUEWZYxRivJyhSAoZv2NFsfn27Z1q02arfntDF+G60K0mg9itpZRCAMpZd3MRbkMyTdX5zqwuhJRhFOysIzCABIYWyiKwGnUSOq+FVv0TOI4oBxlTX3xQLJaseu8xTHlfnljdL1l2xhdkK6m+VZmFceqmQwh5olSWNfvi/liV/GpOJder5H4nVk3Q3xN6bUqpexk1s24OuUy5TL0ooTFVTv3Dx+oNRzW0yrxcib1BkIIYSgJjC21/G4Z58zSrWws5kYB0Y2lWTfj6oQhYexCq2V9mX8QQP3dOm5dHkiFEMJUEhhbyvMAz+49m61l+m4UI1S1wvJa/Y2Y/Vz+ZgV0d9bNEEII8YIkMBZCvFpra7RKb0ZpgRu4wBvwBCCEEJZ65millLqulPoNpdQnSqmPlVL/dv/nDaXUrymlvuz//3r/50op9Z8opW4qpX6qlPr5V30QQog5ZlkGXAghhL2eJ43TBf5drfV3gT8B/BWl1HeB/xPw61rr94Ff738G+BeA9/v//GXgr115q4UQQgghhLhizwyMtdaPtdZ/0P/zAfApsAb8MvA3+7/2N4G/0P/zLwP/pc79CKgppVauuuFCCCGEEEJcpa9V+KeUugH8HPC7wJLW+nH/qyfAYCn9GnB/6K896P9s9N/1l5VSHyilPtjY2Pi67RZCCCGEEOJKPXdgrJRKgL8N/Dta6/3h77TWGtBf539Ya/3XtdY/1Fr/cGHB7neeCyGEEEKI+fdcgbFSyicPiv8rrfXf6f/46aBEov//1/s/fwhcH/rr1/o/E0IIIYQQYm49z64UCvgbwKda6/946KtfAf71/p//deDvDv38X+vvTvEngL2hkgshhBBCCCHm0vPsY/xPAv8q8KFS6o/6P/s/A/834G8ppf5N4C7wr/S/+1XgzwM3gWPg37jKBgshhBBCCPEqPDMw1lr/FqCmfP1nJ/y+Bv7KS7ZLCCGEEEKI1+rNeB2VEEIIIYQQzyCBsRBCCCGEEEhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEAN6sGzALB48O2Lu9RbSQUF6ICesRSs26VVdn4+N1HM/BDVycwKOyls66SVfmaOOY9v4pOtOgNem1Kn7Zn3WzrkT7sM3Bw32Uo1CuQ6kaUm5Gs27Wldm9s4PuZWgNynWov1ufdZOuzPHmMWf7Z+hell+XaxX8OJh1s65E57jDwaMD0BqlIEhLxEvJrJt1ZR7+9l2CSkipFhE1y/iRPcPi5mebdI/bhPWIsBZSqoYox47B7mT7hN3b2wRpiSAtETUivNCec7dzaxud6Xw8cBS1dyzrL/f64ziQrKQEyfz0l/ZcRV/D4cM9Dr7a4uCrLQBaP/8WC99dmHGrrobONJt/eP/8sxeXqPzF78+wRVdr54sNDu5snn/24/etCYxPto5Z//Hd88/V9xcpN6/PsEVX6+nv3SXr9ABQvmdVYLzz5Sb7tzbOP3t/+hvWBMYnO6c8/dGd88+Vd1vWBMadozb7t4f6k2qZb/wvvjPDFl2t/TtbnG0dnn9+53/+PcJaOMMWXZ3jjSN2P3ty/rn+3RWWf351hi26Wuu/f5+s3QXA8RyrAuPdW1vsfbl+/tn9p96bq8D4jSylONs+KnwuL8QzasnVy7pZ4bNyLTvFo6l9rWfTjldg8PR8zrHs3Nls5FyNnUsxl062jgufSzV7ZmgAOvsnFx+UolQpza4xV6x9cFb4HKT2HBtQHNtsmtKG8VmLORvH37yRV2vOdoY6Q8chatjTGQ4ycgOO786oJa/ISAdhUwCie5Y/1Aybs47wZY129EpnU37TAhadu7Pdk8LnUr08o5ZcvfZRpzAeBBV7yigAOod2B8YW3WZjNCPj+Jwdq8Uj72Sn28eFrGpQK+O49nQWvU5xQLYtMJ73J82XMRYYWzSIgXVJjyKbH9iw98QVkiTYlTE+2R4J+i06NhgPjG3KhgOFsc36sWDOxvE3LjA+3iiWUYQtO2rlBkYzxtZlHd+gUgrrzt0Qi04bMH6ubAqMR9n0gNPeLQbGYcOejPFoNjyo2FFbPNA9ugiMletYs9ZkwLY+smDOEwn2jrxTnG4VA+OoZU99Mbx5pRQ29R6jGWPbaoxtCqjGjHb0o+fSZJbec7rbo33YvvhBEFAq29Nfnu2dFj7blDHunnYLY50XW5YtBrTNGeM5Px67Rt7nMLxCF+xaeAcTSik8y06xpYM0AJndpRTDbDptYHfG2LZzNdDePS4cW1CzJ1sM4xljm9bSnO4Vyyj8xL7A2GqSMZ4fut3h9OAiQ6DCEqXEsumXrmSMTWV7KcWblDEefcixicKOe+50x94aXK2LO1I4nmPV4rT2fjEbbtOxQX7+lM27UkiN8fw43jhiuE8P6rFt1xu99psVGM/bk+bLeKN2pWDu+sKXMnaurDo4yzrJvrGFdxbtSHF20C48nNlWX2z7Vm2j3Yd1t+Ccj+N2j7wjTjeLZRS21ReD/cGVzbtS2F5KYV3nPsziGmObbrFho6UGNi28OxnJhgdVe7LhYP9WbaOsGwvm/Hjsipqe4WTT7h0pYHzxnRvYnTG2adS2/QUf87535cuwucZ4jA0nTmvaexfBo3Zcwsr8vHnraxt5qD6zuEwEJmzVVrUrI251/8GEsWDOjteukXfU3t7Fn7UuvPFOK4e4ZVdnAeNvvpNSCnPIPsbmsnomY+TE2XAe9ckp7fbQqv+kjG/wcpPtO3vsf3z//Loby4ab/hro02JNcSEwVsr+rdpsuOmGSI3xDD3+4CHdO3ln0Ts4pn128R/fTcsEJbsuNpiwXZttu1KMZlHn7IZ6GbaXwYyy6NSNBcZSSjFnzs548NnheWL1dPu4sN6kZHipgSoFPPzDdR7/xmd50D+yOK2wI0W7jWke/vgRp1/cgyyje9Yja3fPv/OTknVJhNGbzvi4OMsuf8X18He9HnS7zJI30//1V+yk4/HV761z7eCEblgsmyg1hz5rnZ841/zsquxjbLDRztCywNj4zv0SWtn7wDbGxGNTirMPv+DL9TXe/sUl6xbeBUleBrL76JiTX/2UzunFZLXrO9Bus3t7g6OHuziVhJU/8fbsGvsCen7InR8/ZuHRPtGNpcJ351u1nZ5y9nibrBQRrdZn0MqXtLkJtRp43vjiu+HAv9PJM+hp+lqb9zI6HdCffEJwfQmazckPMlkGGxuc3HmC9zPfxZ9hdGp1YOyWPM468NWHB5TLB4XvotiBJ0842z7keOOY6Oe/Q5iaF0R2z3p4pYt2j5ZS2F5jLKUUhtLayPhqGpsyxr0enJ1BeRArXlJK0enkcXIw7+W5vo/vadr3HnBr95DI6xS+LtUizs5g78kJh48PWPvBAqXQnPsviH20clA64+wkK1Rw9roZX/z9z0GDVop3f+kbM2vniwrSEkcaNh6c4T25V/jO6Zyx+Y8+4WDjhJNTh3d/+Y/NqJUv5+GtU2qdD4nfXUJXWuO/0OnAkycc3Nmkd/0GNXPiYlzf4bP7JWoP7tJcfoSTjTyI7u1x8uMNNh93OMzKfPuXZlsaY3VgPAgYtYaj4ro7Tm494svPodsB5/oa3zR0P+PDe9ts3ztk+Y+vUa4FU0sp2sdddr7aY+m7zVk088qMBiCje6pqbW5mUmeWvZzl9BTCi9rGS8/LyQlE5k5n27Rdm+vCrY9PKTdCVlcn/06nA0+fwsaTHt//ZhuCOT93SuGFHhx1cfZ3ORv5ev2jdZ7+zgFOp035xqJRQTGQ10cHAZydjn85dClW3m4QVs3bwWF4cd3oLPvh5imD/abKb7cIE0PDmmqVex88Jd54TK3xtPCVap+x/48/ZHNDc9pW3Pi5yowa+WIcB6KVGnuf77G318H39wrf7zw8YqP/rBq9U5n5GG7oFfR83NL0wzvpr03I/IDVn1mc+Yl4UfFSwuPfvcdX/99d0m8s0Tsr9hon+x2e/uQxR3e3qH5jYUatfHHHhxlbt3ZZ+nY9rwmfUkpxeKB5+skWN74d4lYN3W3Ethd8bG/nC2BXV6FaHftaa/LvHz1CpxXUtbXX38YrYlPGGKB+eJ+dRw6fPlkhjJxCBnJvV/PwJxnu1jq1k6f433wHmPPAGPDjALYm1y6q7S0UeUnM4s8sv96GXRGvHNCbFBj3aaVY+JmV19iiq/M8iwe1UrS+t/TM35tXldWEPeVwdJRxdFTsP06OMh72k3tuLaVs4Ox25a0ap5/fBfIH62Hdoc/J2vhY8bq9sYHxQOmdNap1cwMQvxIRRC7tkx6HXz4e+/7Jb3wKgAKq7zRec+teXjlxuP/VOjc/f0T6/jKZW8zsHx9mbP14m9Pbj4i9M9w//rMzaunLs62UYrNTRd97TOPgJm5aRunh49Pw2WfosyN2d2F36TrvXJtVS6+AZYtCk7UqBw/v4x3s0otihofh7sEJ4f6HqF6XdMmBxIwHUS969qxg/M4iUcXM2UOvHNDbmf595e0GpYp52WKAUtlFez6q25n6O+FynaQx7zU90yWpyhME+7uX/9612mtpz1WrNj2elBOc48Opv6Mdl+rq7N8vYW5E+Bz88PKnql4Us/I984LFUdHis4uN3DgkWTRzgUl5rY5qn3H48V2OPrxd+G7v9hbtz+/gdM6IF2PwzH3W00PBlFbK+MA4WYrZ2PG5eRM27h4XAn+lNbuPjrh5Cx5veITN2XeGL+NZGWPT4uT0eo1Bmtg9KdahqW4H1cszr+lqasx+24MFatNox2Xxj5mbcTxfhDaFqdliyC8xN748a9z4jrnnDvLJ0PLKs7Ollbdqr74xr4DnQbhy+aJIv1nBD2Y/7pnRo72gZ2WM429eIzZ7PAYgWnh2xiZ929wHgOGOoJh1LEpWzaq7ardHAqZCMKUKVSOmBVaQlxd7zSpZli+4HtmKlM2NfAqtF1eo1mbfGX4dvR7cvXtxTNP2MT49hdu34ejQrBPoRgFR/fLyCN+H0uLspz2f17P2uo3fXSJKzX2w9srTA//U4GzxwGWvtfbqKbVVMxM/w9K1Z4xhcUxSN3NGA/oP3JeYlzHc6sDYC6d3cr1KndVvmjEF+CzJyrMzxrV3zQ2M01YJHV3e6Xk+REvzcVM9r+Nj+PRT2N3NPxeyjI6DUvkONuvrcOvWTJr40p5n2s9p1C52QDCE60Lv0VM++f0Tbt+G03axK22fab76Cj7+aY/jm4+ID59O/hfNsXitdun3acrE2vF55ccBTHn+0q7H4vcXX2+DrpgfTwmMldnZ4oHLAuPGd82sCx9VaQVkpekPpMm1mrHroQCqCwFZOL2zT6/NR39id2BcmlxKoZWi+t01SmY/QJ8b1BlP/b6eENXMPVjHgfJq7dLfSSoupqX/0/2HtDf2uHULPvsMOkNv4tKOw/o6fPgh3L/TIe7uXfJvml/Va+n4Hr9DtFLEK2Y90AxUKhDd+YTjn3zJ4y+LdXOHez32vnhKdOsjaiePUYlZ1yY8R3ZnMTJgn7YLKvCnVlrF7y0TJeYtaBo2LSOcvtWw4pXJ047PiSPqb5vZh4wKAvBb04PDZ92T865UAn+hNvnLKCKuzUc23O7AeErGWLcWWblhbqA4RqlLazRNLqMYSN+6vDYpXq0Yt0+bG/rUNm9S+upzTjYOC29mBMXmg1Oc+3eJbn5ILZntm4BeVJw6ONXpMxpZlFBtmBmQVN+uAeAe7RPe+6Lwnep2CNYfoHpdkoo5C9SGhY3y1Ol514Xy8nxkd55bEEx8aYD2fBa/Z96OPaPyvYxH+kBLssXA1OC+9u1lU8rcn0u8OuW+CkMqi+Y/4EybRYxXq3MzhFt0OY1Tvjf2MjvtejS+t2LyGq2JykuTgw+tFPV3DXwL0IjqUkgWTO4UlILExKxjmpIk4J4cEt79vPCV6raJbn+Mt7tJ4Od7QJrqsin5LK1SMfDUQb7Y6XnemBYvp8Y9tA1MG6STBFTNsMDY9/P9fkeU313OX/hkON8H7RUfZNLrdSuyxQBhJRgP/IOA1vvmj2/DKisx2hlPFkTLVSseAGorEZk/npicp5lDC/4zX0Kp/HWYQ/TyCkurZmaoLjOtzjhcqBCUzX8K8DwIl2sTv4tjcGrzc1M9t0umjtTQart4pWL068or16cHUOFyzeRDI7l++aBcKoHfMOgVVSOmZXeSqnmlS7guflgcDzI/YOn75meLob9zQzQUGFuULQbwA4UKi0F+9f1FXM/Mh85pBtu2jTK9jGKgXM7XlQzTjkt1bX5m1ewOjCnuTJH5JZa+v2DFU9cov1rGL40fWOWG+WUUA9OCkLgZGlXrOCxaTJ8ZGFbeNjsjUmn6ExdPZkFIZcHskqZnnZskAWNT4uTbsSlv5AJV/U34DcyC+yOlIfF7K0Rl845jGi++uJ/S63VKtfl/8crX4adDgbHr0vr2hFcnG27Stm3a9eZif9+rMhrkl1opnj8/96GFIWLRcGDsvHWNZmt+/uNfKaWIWiM3juNQu1GbSXNehdpqmcwfD4DnZYuXF1KpXJp483z1zN0B5p1Sk8speknVpE0NJgqrJbzK9HKKpO6b/aprR43t9pLE4NTNPHHDW7Zlfoml7zZn2Jqrd74zhWXZ4oHhnSniGwuXLjo32ejuDOFyba4Cx5dVXUvQ3sW9OLWuekbemMC4V05Z/a7ZW508y2idcbRaxwvsOcWlEnitYoYuCCBoGRwYp+mlgXGyanYZxcCkbXjcRtXkmPHctJkMxzFvC8FJRsspTM6CD29pFr+/alW2GC6Oz8ZsMVwExlopFr5n9vZ6l6k0PHrRxcBgSxnFQJKArlyMCc/cv/k1sydqmsIN8qAiePea8dmpZ0mWizU6NpVRDIx2EEmq+huqGioI8lKQKZK37DiH1ZUyeijbrx137Ho1VfXG5MC4HIOqGHxt9lWuVwuLnpIlc98wOcgYZ0HI0rfNLlGaxE9K1maL4WJnimitaeyru59HEEDQzINFrRyq18zvR4YpBfG1/P5T5Yi4Pl+lkNYHxl7o0a02WXvfsDcIvAC/Fl+8TtH3qV2362aCCVMwy+a8knYav5FO3FPb9VRey2kB181XVQ/0kqpxb7ubplwv4Sbj2TmTM6vD/NCl1MgfYsLo8n1W550TBriundliyDPG6bWaldliyEuXAFrfM/v1z89jUF4QtCoEodlj3CTVaynacedqN4oB+/5rj3CjgPj9VeMWUL8QxzmvM46v1XFc+zr+OAZVrwF5PDyPN9XXVqlM3OY2Xq2OL3wy2HA5hcnbtE0yyH4UftaKmLg/mIEGNeJJglFvuxvj+3iVspXZYoBS4rPws6uzbsYrUyq7eCsLpAt2bEF3mcpKjHa953p7qIkqVUWWVucy+WN9YOytLLD2znyl6V+lQZ1x9V27FpUMGwQhcQyqakF0NaXO+FkvNTFNZS1/C55WivJyxfREf8Ho7hSl0sVUqA0GJUxp3ce493cPCwLq31u1orZ9kqisrM0WQ54MWf2FtVk347UY1OFW35q/wPEquC6UrzWorM5fSZ1FQ9NkSdW15tXPzyNeSlBRSGXZ4MHrGaprSf4kbfiK/3OuS3khLiwMdVw1N++NvypB6BC0KmRhTKVhZo3qNEkrRJUvrsU4xuza9xFxPaBUiwiXDL8mo4j6DcOP4Q0XV+yZRbuMUlD9zipRaldfOWzl29W53Ifa+sD4TRM0EprfbFq9+0ZaUVCtmr1N2whVLW7bVl6xq4xiILlWo5fWjJ6Nn2R4MQlYsCh0gpXv1MwuowAj914Wb66Vt+2e7Z7XvJYExrZxHFrftXcbG8in01rv1/Gahg/Sw0bKKUx/qcc01beqeM0qoYUlgoPSF6WgvJgYvyh0VLRat2IxoRCmsGCnTiPZ1XMLAJRr/2ldeK9i1yCdJCSV/Lwp17GujGIgSj1a1yyMioHKYpi/5juxpPZ9VBRZF+wLIcQo6eWEkZSj7HqcVoqgkeD75AvTfIuObcSipRMajgPRSo0kxq6HNiGEeINIYCzEvKhUiBNI37bjpR7T2FzmWXm7TlLz5rd4TgghxKXsXe4ohGnSlLTiEF23s4ziTVBfjVBqwe7oXwghLCaBsRDzolwmebsJvkzkmEopYHl51s0QQgjxgmQEFmKerNr71qo3hixQE0IIY0kPLsQ88WQSRwghhJgVCYyFEEIIIYRAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEIIQAJjIYQQQgghAAmMhRBCCCGEACQwFkIIIYQQApDAWAghhBBCCEACYyGEEEIIIQAJjIUQQgghhAAkMBZCCCGEEAKQwFgIIYQQQghAAmMhhBBCCCEACYyFEEIIIYQAJDAWQgghhBACkMBYCCGEEEII4DkCY6XUf66UWldKfTT0s39fKfVQKfVH/X/+/NB3/55S6qZS6nOl1D//qhouhBBCCCHEVXqejPF/Afy5CT//q1rrH/T/+VUApdR3gb8EfK//d/5TpZR7VY0VQgghhBDiVXlmYKy1/kfA9nP++34Z+K+11mda6zvATeAXX6J9QgghhBBCvBYvU2P8bymlftovtaj3f7YG3B/6nQf9nwkhhBBCCDHXXjQw/mvAe8APgMfAf/R1/wVKqb+slPpAKfXBxsbGCzZDCCGEEEKIq/FCgbHW+qnWuqe1zoD/jItyiYfA9aFfvdb/2aR/x1/XWv9Qa/3DhYWFF2mGEEIIIYQQV+aFAmOl1MrQx78IDHas+BXgLymlSkqpd4D3gd97uSYKIYQQQgjx6nnP+gWl1P8b+GeAllLqAfB/Af4ZpdQPAA18BfxvAbTWHyul/hbwCdAF/orWuvdKWi6EEEIIIcQVUlrrWbeBH/7wh/qDDz6YdTOEEEIIIYTllFK/r7X+4aTv5M13QgghhBBCIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBPAc27XZSGton2m6Byd0t/fxSi7xDXteMrL36AhHaRwHPA9KzWTWTboye3d30b0MPylRqoZ4JXfWTboynaM2xxtHOIGHF/n4Zd+q49t/dIhS4DjgeoqwGc+6SVemc5bRPuqgz9roszbRYopXDmbdrCvRa/c42T45/+yXfUqV0gxbdLWOtk7R3R5kGUpnxKvVWTfpyhw9PaR72sWLfNzQJ4h9HFfNullXon1wxsHDfYJKSJCWCJIAZcehAbB3dwfHc3ECD7fkEVbs6E8ATndO6Bx3UK6Dch3Cagk3mJ+x7o0MjA93Otz/7z9B9boAVJYiqwLjh7/xBUpnAHihx/v/0s/OuEVXZ+ezp5xsHJ5/vvEvfo+oHs6wRVfn6Okhj3/7zvnn6reWWf2FtRm26Go9/M0vIcuvS7fk8c1/2Z7rcv0nj9n/4sn557U/9Q6VG40ZtujqHG8c8eA3vjz/nL7T4to/+fYMW3S17v/6F+h2BwDHVXzrf/PzM27R1dn+bJ3D+zvnn6//c98mWbLjgfTo6SHrH9w7/2xbf/n4H99BZ/l2uk7J51v/8s/MuEVXZ+vTdfZvb55/Xv3T71O9Xplhi4reyFIKv+wz/GiZHZ1ApzPDFl2dXo/CsTmWneHu4en5n7VSVmWuuqfdwmcvtOe5Vev8nwHbrksV+IXP2Wl7Ri15DWxKywG4Q5mq0QvVcL3T4rjmRf6U3zRPe/+08DlI7RkLdLd3HhQDOIE9YwFwniAZUO58DQjz1ZrXxPOgV07PP3e7wMHB7Bp0hUauN1zHnk5ed3t0hoJHJyxZMy0IkJ0VBzGnZM8gNhpr2BZbqdLINGfb4sDYMlqNDIM9e17WOtqnuKE9fUr36Kzw2U/tmDkEyNrFJIny7QqMdU8C47njutCLL9L23R6wvz+7Bl0hmzPGnYPT/CXkfX5iT4YAoHdW7Azdkj2d4egDm03XJYwHxlZljC3KoE7kjtQ2jl6sBhsOjLXj4vn2PJEOzx4CVs0ejmb6bcsYS2A8h5QCp3qRMe51sSswHuIoewa1s71iR2hThgAmZHcsDoyVY88ADeOBsT6zKDC2nLUZ416PrHtx46nAt2qmpnN4kTHWjksQ25MNH80Y2zR7CEgpxbxyo4AsyAMrrftByenpM/7W/LM5YzyaIbAtMB7NGNtUD5j3gxcPaTYN0JBndIYDrMFiLivZdvJszRh3OoUY36bgSp+16bQvzpMTlawa60bHAimleL3mqzWvkedBLx7KGltSTmFzjfHoYgvbAuPRLIHNGWObBjHIj0d7F4GH7nTtCbAsZ2vGeGw63qLAuHN4Viyrs3wssK2UQjLGc8rzIButM7ZgAV6vl+/WMGBTANI5GKkpq1rWGVpcDzhWSmHPoQH9wNi/KKfQGnsW4NleYzzaSVryQDO6y41NC+9GkyRebE99MUiN8azNV2teo9GdKXqDnSkMHwSsrTHWurAKWbseQdmiziLL6A3VA2rPH5vhNdlYxtii3URgPDDOMuwJjEdZ9lSjnYsbTSkkY2yA4fpiAC+xK0miO/Y+1ACFAUErZ+66lDc6MMZ16UX5Zue9Xv//HB9f/JLW/b3czJFfb0MZY1uCq3abTvsiyM+CkMCeFwFBt0tv6FKzLUOQZaCGHjptmsmAPKAayxhbsje69SzNGI8GxjYFV6OzhzbtYQz2b9dWuMccCYznhte/zgblFOdJgkGdca8HX34JZ2fjf3mOjSY7XAfjs+AA2fFp4RnFi0tzdzO9DN3uFF+AYVF2B+zfx3i0xtjqjLFl9Gj2wJLA2OZ90Uf3MA4qdmWMbV5vAsVSCi2B8fwYBMaDcorzgPLgIB/QPv88/7Nhqa3RXSlG15WYyvaFd6OrkG3MGA8z7LZ6JseBbLiUQmqM51eWFWcCRy/GwWCgNayvw8nJ62vbFbI2Y6x1cas216NUtmVqNGd1YJxlFHoUKaWYE/v7eN080MqiGK2cfPEdwOEhO7/zGUeb/c7QsBHc1oyx7YHx2EIZmzpC7N/HOM8YD9X2mBwY93rw9OnUzKkeKtVib8+4wLHdUfDpp3DvHpyejmeMe738uD75hM5XD6Bk5jS9ta+DHiurK9lVVsfIQmzXMy0MuVyWFXYUmcdSCrtG30toDbdvQ60G9V6b4Msvie65ZFEMjkOvm3FyCk+eaE5POnzjG/2/OG9n7BmyDCv3MbZuD2OtC+fJ5mlPePN2pZhYStFuY8IIrh2X7pMt/KdPYXkZsgnZuN1dePwYfXxC93s/i0lX6/qGortTZWlvg2hjA399pC/Z3OTk7jqbmxA0EpYM7UQLwZVSuIGhWVWt6WXqfDGyPj1jeG2aE4XWjHNA/mKW3kXkaNtCbLKskKuTxXczpBSEIXz1Ffz0UYutswSV9XCP9lG9LkdH8NUdOD3Jyyz8QU8/53fc7i4cHl58HtuVYihjbPJi6+HFFlop81//eXSUZ6z6J2V8a6WRZ1bDs/7WllL0z4tSgOue74mbDS++63bzjufp05k08etSCu5sV3lyv0Pnzn2c+3cL37ubT9E3b7Hz6JgvHqe0e2aN2tUq7DhNvvoK7t8Hzi76liyDe19lfPVV3q8mK+m0f83cOby/w/FRv5/QulCepT3/vHzQOO02D37zJttP8/upc3Ba3MM4MXwsANr3nlwM0N1u3n/0adczPjA+PBwaA0YC49GM8Tws67JleHouzWb+/3s9WA/fKuz3Ozxwh9HQX5rzEdxx8nLo27fzC2rSm+9OT+Hu3XxsNsb29sWfez3aJ0OdvF8iKM3ZI+bXFUWwsQEffwzb2+Ovgw68PKB6+jQ/wYYbzoDA3N9Wz+/JE7hzB+dgD7Q+zxrrjPxm3NzMz/HWVn7ODRGvVtnZgZs3YWuzeO5OTzQ3b8GTx9CJa5TLM2rkC0oSUElM5gd5UmE4CNH5MyvkL8QrL5kTGLsHu9z51U+5/+kh3ZMOw1vFDmcd99dPuffJoTnP2kFAcLrPk3/4CXf+cJfjnWLkNJg9bJ9p1r/cY2/bvAzQ8fohd/7+xxw92MkD46FDOA+Mu13O7j7h6Rd7M2vnizo51nz4YV6yn3VHMsb9wPjwMN/v4MmT2bVzwNRnyBcShnmneHgIOozo1hfxt8ezONHwzNq85fhHJAl4+9vsHbrs7lZRCoYnax8+hE00ePBO6wDu78L167Nq7nPbu7lBOd7E/+Y70G7THYobdVAytezvnHZcVBjmTy137pDtFK+zYPNRPkJrDUtLc38djjo+zuPCtB9XXFZK0evln00LltttCJpN+PBDnK1topv+eQb5/Hjv3qXdzk9zxaAIsrISs+d657Npw4YzOsla1bRLE6WgUoHDSgNna/oonKQKlcSvsWUvJ2yUcdvbHP7+53z5ZQVvuPTbcdm9tcWDrzbp7R1S+cVvm3PelMJPSqidU04/vsUj5QxXudM5POPWr3/F2dNdiCK+9b/81qxa+sLKCzGPP9vj3m/eJm2VCjs9ql6Xk49v8/TeLgeHisqf+pnZNfQFVTdu8mB3kfvdKhvdDHd4glQ5fP55Hpc5x4e8+54CZnvfvVGBMUCrdVF60FlYxT3YwekUawHPEztKzX1A4jgQVgKyTz8nC8t0msuFTPjxMbgnhwS769SyA1hbmWFrn19blXj60y2Wtz+hcq1SLNcMQ3ynBzv7sLNDu7pA0DQnswN5oLS+GbMSnBIE0BvJqLqnh9B/QOtUWzg9jJpO83347KdtKq2Aa9cmv+Cj3c4T4vv78L1vdsAxqVI1z2ycnQWseDWS7i5q6OlNkyf8t7Zgeycv1SU0py4+SfvR48721N/JwjLVhfmvmZ6kVoPdahP/ksA4XoyNelpTSUwU9bfi399nOPZwj/bZ+0m+FWlWimhcNyfgh/5bTu/nJS9KFzuT7t2H+c+B9JvLRvWTA0E9xvfz6quDzWJG3D0+4NHH+Z97lTqNBfMOMCh7VJ7e5GQrpZtUcUZKKU4e7VDafop/dkT6C39sZu0ceOMC43o9ryvr9QDHobN0ndKDW4XfOR+/5jwoHogXY/a/8HBOjyk9vD32fenBLSqVfh9frb7+Br6AUqVErwcP73bZ294ulrqc7ND+8Tp7e5q9s4j4F1LMCPcvRBHs92KOb2/Rao6/R8Zz8wFuq52y+1XI978/m3a+KN+HZP8RJ08zPnm6ghuHhZmM9XXY/QjU0SHL+jFspf3o0Rz1OnzxBRweLVJb381H5n6H3+3CrVsXDwRJMzQryFJ5OcXxJYFxL61RqbzGRl2hSgV0KSQrRThnE3bVUGbVFwMQRReB8SX81QWTqnoACKvPniLMShEL3zBjfBtTLlMu55uhXMZbbhlXugRAmhLH27R3DnCPDwpfebubeLubAMQrCU44+4dtc3rqK+I40GhcfO6lNXrJxc0UBEOZOUMGsrSi6KW1S3+nWiVfVWjIXRXWLrJrh8X7iM5xm1s3NZsbcFxZZmHhNTfuikStGHRehrq7W/zuyZO8LnxLtajXjXlGKyivVPEOdojufIJzt/jAdnrYJfjqC8K7n1Nz9vMUnmGSJH8AyOKUYx0VtyDSF0Gx74NfMSwSAdJr1cLs06jyStXI7BzkXWGSQLfanPh9XAa3mrzmVr0k1yWqXz4roZVD/b3Gpb8zj/z02TtPlN9ZMrfEzvMo1y9vfOaXqF037GFtIE1JnuN2Sq7XX31bnoMZkd8Va7WKn9tL189Xk0cGLbwbSBLoXhIYuy7EMXmaxJAIK0gvebNdPwDJ/BK1d+rGrraOW9HUwKPbzRdd9NLa+aJR0ySrlfPj8w52C9+pXhf3+ADPJx/MDSozGFDqIp7vNBan/l6cYMwD6bBqw0VHk6fcMz+gsmzeMQ2rVvOp6UnSiuK5RvI5EzUvPye9WpN6y8CnmdLl60oyP2DhW+YF/MOi1uXlLd1aq5DUM0qpRLkWXB5+KKi8LYHxzJTLxQBYByU6rXwyvhAYGxJEui6ErfQ8uB91Hg8bUkYBPLMjBOi0lllcMuMcTZKkiiycPpB1Kw3KiWPctOdAWnPJosuDi2oFs67LEYOBqldpjL8ooi8uY9SOFAOuC+HS5HPTS6omnzYgv+y0H9CbcI2aVl884FXKl26VnbyzYGYiIQwvHQ/8taW8Lt5gpUY89dxopQjXmuZmxAGnVrk0PxAvpXOzf795d/4VGc0ad5tLZEFo1FZtw5KKQy+ZXPB3PoCZVBDougTl6T249nzSt5smJhrPlctcHhjXF8auU5O4LgQLl0dPaQUjyygGBuUUOA7d2uSTVTY0MIZ814lJ/FbN6EEa8lMSBNCtFadkSqX8xR5GKpenXmq9KKZ5zczrEN+nFE0ej7Xr0fq2wR3lQBxPDRx7SZXG0nwEjS/sGeUUyVvzkw43J/K7Ys3mSEJYKbqrbxGWij8zRZJAL6mN/TwI+mNyHGNaqiBIp4+8ncaS0dli6GfkGpOnz3pRAmFIfT5mll5YvDo9MPZ9iBKvX+djrsE56k4opwhD8EquEW+8m6S2EpH5xbZrx6WyZmit44haLe83h0ua0pSLfQZNUy4TTQmu1MKCsYcFEFSmZEEWF6k3LQhloohyPHlM69UXjB8L8gV4U75TispbtdfZmktZcDW9GNdl7EILF1JUayh7YFDGOE0pLCIcqAx+ZOC857S322nXI1gzu5MfCJuTe4purUWtZtyzzJh0ISTzJ5/HarX/fwx6AJ1k0I9oPxir9Tc5Wwx59tRvFvuOXlyhWjP7nA0MFiVn8cVsWpKaWV8M5AvwJuzgoF2P2jt1o2+1SeOBVg6Nby0YfVznHIdyc7yvyPyAdDU1fizA9ylVQ7wJie9oqYI3+rbXGTIn8nsFRqep4xi4du1iWwqDAmPPgzDx6JWL0WJ10N9bFBh3GossrZhzbi6TtMKx2lTtuPQqdWMX3Q1LEsjSyddemmLkdTkqSS4Swt16MWtcjjFy4d2w0ay/qtdMT/KfS9O8mx/sTuF5/QVsBvX9o8JGeaz53WqT5oK5xwT5zhSju6BkjRat5fkJqF5WqRGPHWO31qLRtCHyJy+nmNB3zMtuFANm3ykvKU0p1MnFMXnPuLaW/8Cwx9A0pbBtW9ivocP3jRycg0o4dgq0clALFkwr9SXJeJ1xt9rALzlGlYRP4zgQLY8Hv0EAYaTMqnu/xOB6zOKUrNTP+ihzF94Nq167WNirlSJeqZjWNU41eAteL6milZMnik2fippQZ1xaaxm9HgOAUqlQkaSVovL+kvmZ1GFxPF4K02jakD/ITagz1kpRfbs2k+ZM80YHxlDMGp9nQRYW8g+GZQ3yOuOLO6hm4qK7ISoc35miW19gcdWzZmAOAnDS4iN0t7ZAo2Hcc9lU8fL4jimVKnkAYupGuCOGH9QGW7dFg3d6GB4Yx6mDquTBYhYl1Fo2RSL9tZ+OQ69StyMwHgmueuWU5prpUTH5zhRDh9GrNFi8Zmbt/lRxnD9M9/WSKrXFwLRQZLo0zXN0Q2NbuFjFD+drHLDlP/cLG0xX+/7I+pi33jIyMNZBKc9YqaH+3dTHzQkZgl5rydgXekxTbl30hL0oRoeRFWUUA2lFje2YUkkxejeKUXF8Mfs02Lrt/EHb8MBYqfxlHgBZxdy33U0zOJ6s1iBOlPGLQSmXiYYCyKxpyQxbqURpaDwov7NkfhZ8VBgSJRdBotF7F0/iebhp8fqcp90oBsyK/F4B38/jxrG+sFyGpaWZtOlFBUE+OPfSGslgEwpl8HS15xFExU6iteLbkmQ8N7yxe7fWIo6Nj6UK4hj0UJ1xqdQPIk19YJviPPjob90Wx+QHatgD9iSV6/m5ipbNfdvdNL6fX6PxcopTq5g/i+G6RLX8KU17Pun1mvGHBIDnUYrz2YpuWmPxbYs6ySGDGnHteriNqvETGGOGyim0cqi+NX/jgPk99hVotaYkCQyMTgZvwTuPOZLE6I5+sABPK0W3scTi9BeMGStpBGjPRyuHXqVhVbYYihlH6D+nlcvGbmE2zXBWLmsu5N2HgbX9k1RaAb1qg8qC4ZsXT1Gtku+0cf36rJtyJdy0TKnUTyYsWFKTxcV44K0u2Rcw9qkk38/YqkV3w9I0fxso4C/Upu5PPUvz16IZqFbtSV6lKThxmaTRDzpMzRb3DTrCXlqntvTst+GZqFwGHZXpVRso17Fj2nNE2vDPFxlWDH+pxzTl8kU5RdL8/7d3ZzGSbOlh3/8nloyIjIxca6/q7tt9t7lDc7iIkihro0RJJkFB0oNgSbYhwhDAB+tBBmwYsl8M+0l+sWzBgADZFkQZhixBokSatASORqS4DMnhcDhzl75L77fX6qrqqq69conjh8isysileqnqzjynvx/Q7IzMvM0TkxEnvvjOd04EqFrVyJvrUVwXgrdXbPzZgOxwrFbBmg6mu56xMzdjVQDplUJUUmLuiqHL6T2POCaKspVErCqj6CmViCKF60JycTovdhIYk2W0LLl+USplHbxTr2ZvGB7xF5IAFLQaC1ZmiyE7/oJ6TKs2a8XaxaP01tkOw26i2PDjcpzehSxJgLk5ezoWYH7FtyZuHBRF3ScY2qJYJJyrUF+wa1SGICB4a8HOgLEnjinOJ4TV0Kbu44TrQrFIXHansowCJDC2ThB0S6Or1b7H3plLhQFeo0qxERm75v7ziFYa6LBoXRlFT7EIVLsTtwoFa0oMBvWy/cdPT7MoXWf44NObpVikeGnWvv6kWmX5g4o1K/aM5PtEVxatm2SekyTULleJitP5Q1qYmxJRBOgSVtxWhyHeSsSMWfMgX1hcK+Dv2B18lOaKlHuzXS0VRdmozXHcb3B9vzCY5xHMWXieRRGW5cBHUuWEGT3pVrxC5TLFZHp3UAJjWykFCwuTbsXZxTHlRWVtbWNPqZQtHWhzJmRhAXy3YmV9cb+LFyfdAiGE6Wy+FhDHU72DEhjbzIZslVL2DQeO4HnGrQ74wopFsloDi8oLRjG8ekkIIV6tKV/CUgJjIaaEjZPuhthcKyKEEMJ40x22CyGEEEII8ZpIYCyEEEIIIQQSGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAYA36QZMUvuwzcH6HspRlJbKk27Oueo0OzR3m7QP2yRLyaSbc650qmnuNmk+PSCaifEif9JNOlftwzYHG/s4riJesOe36zQ7pO0UnWp0qgnKwaSbdG6au03SdopbcHF8F9e3J+fQPmjRfHqAFwd4xQKOqybdpHO18+UmXhxQKIdW/W4Ah5sH6FTjRT5e5KMs+uk6zQ5HO028wMULsvPOJs29FmiNUqAcZdV1rn3QIm2nKNcBpfACF+VMz8H5RgbGB1tHfPmNa6QHRwAU5xOrAuPrP/8JrZ3DbEMpvvLXfmiqDrqzWP1wlY2PHqB0CsDiH3mL6tuNCbfqfBw+PeLOL39OetQCoLhQtiowvv4vPyJtdQBwfJf3/8oPTrZB52j1O/fZ/fLJ8faFP/M+pYXSBFt0fnYfbPPwt24DoJWi+sESSz+8MNlGnRPd7nDv124eb3u1hHd/6r0Jtuh8rX34kN27m0D22136c18hni1OuFXnY+fhLg9//frxduWdWZZ+9OIEW3S+bv3SVdJmGwDlOnzlr/3QhFt0fla/+5DtG2vH20t/8l0qF6YnBrPr9vg5FUoFOofN4+3m5h5oPcEWna/cnbPWtA9ak2vMOXMK3nFQDNB8ejDB1pwvv+jT6XaEAM2t/Qm25vwp96S7SdvpKd80Ty/g77Epe9XZPzp+rbTGLdizb8cJhC7bMsbt/ZPrnNIaN7AnF2bzOQf5kMSmTD+A7uT7/2lL3NnVCzwn11N4ycldc7uZku4fnvJfmMUrFnLbzd3mmG+ap1CJctutbXsCY9d38OLweLt92D7OHluhv3e36EYUhgN9m4LH/sAYwIstKoEZuLF2S9GYb5opPcz3/X7RouH4o3xg7AX2nHODbAuMSQcCY3e6QtHpas1rVKjFue2D9b0JteT8+XE+MG7tWRQYl8PcdnvXnhsaAL+cvzAfPrEnazzY+enUnuBYt+3NXg0Gxm4cjvmmeQYzxn5iz76hNZ3Dvhtr38f17ImwOs38OWfTzajtBvt+4wJjpdQ/VEo9Vkp93PdeXSn1daXUte7fte77Sin195RS15VSHyqlfvhVNv4s/Fq+/u9oY3dCLTl/gxljmwLjIHJI/ZOMVWuvCZ3OKf+FWYJavv7vcNOejPjgcJlN5RSDgbFNQ/K9uRiQ1alaNQlo4MbaK9kTGOtmK9c1ulFh/JcNZHtgnBtUsyxlbEMpxT8CfmLgvb8NfENr/S7wje42wE8C73b//Azw98+nmecvbOQzxs1NizLGpfxQZ2v3aMw3zeM44MQnWdVWGziwJ3gMqvmMsU11xjZnjPvrHbVy7Fm5odOhc3RS9669Ap5vyb4B7Z183zE4ImWyzkE+IeKEdgfGUkphDoXhGWOt9a8BTwbe/ovAz3Zf/yzwl/re/8c689tAVSm1eE5tPVdRpYD2TjIfze1DazKPhVK+A+yfgGGD3HCntmsCXljPZ4ytCowHsgKDWQOT9WeMlWfRBfroiPZJXIwuBHi2zN9KU9p9wWPqFyhE9vx2rf38/ATbM8a2BcaWTcPIGcoYmxYYjzGvtX7Yff0ImO++Xgbu9n3vXve9IUqpn1FKfVsp9e21tbVRX3mlwhA60UnW+OgI2LMjazwYGNtUSgHgJQMT8HbsqTMOS17+hm3naGiigqkGO7+0Y0nPrzWdvn2xqb6Yo6NcvsCqwPjggHZf7KiDCN+eKpGhjLEbWrRz2B0Ya00uMrYtY2xDKcWptM5+wpf47/6B1vpHtNY/Mjs7e9ZmvDDHAb9yEhi329B+akdg7AUuuCedxGAHabrBCWo2ZYxdF5zSSda41YJ0z479szVjrNudXA9oXca4v061aM+KFBwe0soFxqFVgfHgSKHNGWOlLLshHWBbYDyYDrclY7zaK5Ho/v24+/594ELf91a6702lwZUpjp7YERhDfgJep5UeLxRug0I5RPf1FIN1gqYbXJLOlpUpbK0xHpoEZNHEu87+US7odyJ7AuN0/zCXDbdp4h0MJ0QGJ2WbTOv8eed65JJBprO5jALszRj/AvDT3dc/Dfx83/t/vbs6xY8CT/tKLqbO4AQ8m1emsGkt4yBU6P6VKQ7a5AohDVeo2llnbOuqFIMPGrApYzy0VJvhGeP+DPHgGuiDJVqmyy3Vhl2BcacDpCfnnecigbFJ+soDtZq+RMLzLNf2T4DfAt5XSt1TSv0N4O8Af1YpdQ34M91tgP8PuAlcB/534L94Ja0+J2HRIQ1PgpCjgw4c2lGvanNgXChAGpxcxJpN8itTtFpGT6QMa/kL9NHWyb7t7cHulpk3AUOlFJZmjE0e0tWpzl2U23sWBcYPH3L9e3tsbGSbg0u1WbWGMflSCq0UfmRLcXiWB1F9fbzJ59zzmLaM6lnl+n7HmbpSkWeeKVrrvzbmox8f8V0N/M2zNup1iaJsAp5zmGXkmr0JeKH5HaTND/kIgqwekJ1s+3jJtiShufaUvU9uU/vj/8FE23gWYSVAK+f40dfNpwfoVLP6WLH6yTrvLu5C9a3JNvIlDJVSWFJjbFXG+O5dPl9vsPRuTDnRueF47fl4hZPfUOvs/tOkyXjxvc+5t7HIk0tzVLf61mf2fPzQJU1haws21lKWVhziePy/NdW0ptP31Ezt+XYts9cG1ZcxtnoNYxv1TyhXyrzA2GZhCDqKYTNbFeN4ZYpGY7INOwe2rWV8sNshbG6jKmWU62aPbl3PPmu3QO/ts/XRPR5/tEpt1jNuWO3OjTZB7FGtQhgp0jDCPchq3g8PUq5/dMDhnVUK208ovD1/+j82pYZKKSxZlcKmx0ErBdGdz7j9eIHkUp1SX7eR+tmKFEdHsH7vkCdrHd77odicwDgICANNYe0Bze1N1vv3LYjYXj1g8+oW6ukW/kKD+P25ybX1rFqt/DJ7XsHoiYVpCmh9vD744OOgTT7njqVptioAw4HxUODY910T7G13iEpu1mStcxljrfIZ41aLiR+rpnRpr4RS4FdjeJBtdzrQ2trDv5ht9w7OabubeR6DGWPT1zLeP3K5/ZtrVN1bJEsJYSeif4/ufGfjuJoiSMyrpSu7ezz491+yFpdx65WshrobGLdboD/5Ai/toBTGPHlM6/y5MzTzuK/3bzbBc7URD8bodLL753I5235Wxnjwf4epFoZEEexvPOJga53D/l1zHLY+e8TmxhOcowOir31AYFJlRaFw3F7nKF9f7OzvwNWruGTZ1cb7M6+/fWe0tQVJ0s0JNJt0coGxb84NzBjXv36L+Q9qVC7VTg2MW4cdaDaHVi+adg9+9z7VKhSvLAD+6OXaWi12rj0i1YrK961Mopkv5ejal3z5xGfm+xeZmVX5wL9bSrG/p3n8+SZKp1z6A5M9/ww/Vc4urIYcOu7xsEzz6QF+mtJsO9y+De+9qwFTrmonTnvIh9ZZKXVkUL9Rr8PDuSW2rn/O1uY2sJ37vL/E2MTAuHqpwtqHPp2tddhaHzoxe8dnodD7P9Pv4MkBq9d3mP1ghlLZGblcW7sNqzd22f38Pu//1DvA9Gd+3NYhdz86pDBXZWVlODDurUqxswOrq7C0qCnGhvQh3cAYQHXytezu3jbsbaOATlxm4UJx+L+fZkEwNpBXfVdqPTdPrWFONg6Adpvml+t8oudZuaCo08xljJ2wcBxc7e9n/X+9PpmmvgzHAQoFHvz6TZ5+UcKZywdObsFlZ+2QJ58/5vDeBu/8R28DBl3gAJ2UufO71yl9tkb57YElbLXm6cd3efL5GocHmos//u5kGvmSSgslVj/9krW1dTYW5lD9XYvW3P7WYw7urOK0miz8yfcn1s6eNy4w7nS6ax52+70ogv0ozjp94OhQ03qwz+31Eh5tuHYT3ntvgi1+fv2jK4XYRyt13OH3AuODA7hzW1N3nxK9V51QS1+cUtC4VGLjQYK7v3PKF80MjJWCxteWefwbX5z6Pd9n8uNMz6lYD2nf/YIvbz7EX57DU/mSg837+zz4ZBO1/ZTSfGxO+UuhQPXJVTYfJ3z6eIXSwOS7vUOXTz9JOXq0SXFvjeLldzEh4AcgDJ9rioWem6dafeWtOV++TxAqUHrsyvva9ai9N2vSKHXG8yhuP8J9+ITbu2+xqVq5rJwKCqyvw9pjzdHqFhc/iAGz+slotsTug1X2VndhNb+C1PqtHdJrnwBQbvi41WQSTTyT8nLC09932N1J2f3u41w67mC3w8F3s1Vx3YJLvGDW/hVmygQBHB110A8e5k4/p3nI0fW72UoQvk9luTShVp4w7fQ/Fx9/DDduwMZGNnEk7XsC3sYG3Plkl3Rvn9LdT/OpyCl3/1aTGzdgczPbdsOTjq912Obh/ZTPv7NH5+NPiVrbY/6V6TU7C525058wXvBBhSaN756oX0pwKqd3eL6PMRljlKK0UkV12rS/fMDhndXcx4dfPkZtPwUgXixPooUvx3EozhRx97aJbl2lvbqR+/jg/hP0hx8RPLxNda5gTsAPUCjgB86pw+5pEFG9WDanPKSPCgOCU06fdmOe2XkzL4tRo4hzdEB4+zMO7+afJts+aPHgd75Ef/d7RE/uZ8elYYpz4wOm/jr/6pW6QbVLJ0plB116dsAbL5XN278goFh99jEXLdWm4qZ0CprwerkuzNQ67H1xny8/fsqXtzp0opMTrt0Gb2ud8PbnBKpp1JTr2XCHw29/xN3fvs9Hv3tIx+07EDWs/94dglufZZ1nzaxhJsh+iuqFJPd7DSoEmBM4DlAKGt+/dOp3TMoYA5Qv1Z7re6UlgwJj8hdp1c6vF+vuPj0uQ6hcNmi8uicMCU/pHlr1eWbMK8HN9NUZD9KOS+nyrEmnV45byh5prbTGaeUnW/sbj/A311Bph+TtuakIPl5UXPFIC6cPZxQKEF80c/K8UlBcrDzze+WVZ39nGj1PH59cfL7rxatm4OlxdnOLLq5uE969TnTtexQe38t97rSOUDrF8zEqMA6X6pSiFH/jEYVrn8BefrjJ234CZDcHpk1M6Jmfh9bs+OAxNDgwhqxcxKmO70D8om9UtiBeSPDC08+hQuQOPYVy2gUzyTO7hrDkEcwZeBELQ6Ix8Yf2fMKlulHzE3KCYGypSKsxz9yiQdn9QVH0zN9FK4f6u2YGjmEIlE4fZi8vFM2aPDOgtHx6f+E4z/7OtIoXktMvXVNSRgFvaGDseVD7vqXuWrF6aIZyj+9hVGCMUtT6lhhSYxZDDAKM7TzCEJKl8VnjQgHMmiqfpxTUv3957OeF2LB0llLEy9VTv1KcKxkV7ANQKlF8xtyz5GLNvP2C3AS8Qa36PDOzBu5Tz5gJeNpxKSzPPSvumm7F4jO79WCpQRibG/yHM6f8QKpbRmGwymwh9/CqQfF8CeUbFJP0UeXk1HXBg4Xa1IRbb2RgDDC37JPOLZz6nUIBswJjoHxlBr9w+oUrSAyrexwwPw+tmdG1xkHRNXrfABoXiqhadeRng8vwmaB8sXrq50bVF/e4LsXG6VFI5YqZmblxE/C045LWZ4xazWBIoUAwYt/a9Tmzs8XwXBMnq+8ZvDYz2QS8cUol8OdNPjizmMOfrY79vLQy/rOp5/un9pnTUkYBb3Bg7HlQ+8o82hufgfMNK6UAwPefeddcKJv9ZL8kgWC2PDJrPLhMnYkcJxvRGOS64IaGZYzJasvGLsCvDA2MgWhu/ESZYi3I1kg3URjiusMDL+3qDPVZ18j61GNBgD/w/B+tHPTsHLXpuS6/HMchrIZjVxfVSZnaotl9f1wPxl6zqxcrRs2/GCdeGlMqoSAxtL64Z1ydceoXqK5Mz3CNyV3cmc0vOrTnxtSrqm5MbFpgDNTeP31yhYkT7waNyhr7frZepw1mL0ZDC40atSJFv+7qFKNEiY8bm3mxDmdKY7uH5JLBmatu2rE/+6iVol2fM3fSXU/3/OkP+tu1WWYXPbMD/i4njsauulF+x8xJd/3imJEJEc+D0kWDz7k+5YUi2h3uWOJ6aGxf2RPMJCPvXYL52lTd0xh+mpyN50HjvcbImh7P65YHGhgYu0mR8tL4uy8bAuN6HdxamU7fUntBiNH1xf0cB2pfHQ78p6r3eAHjVqcoLpiZLQbG1xkrg8soIOv4giBXr9op1wnLhVNrBI3gZeniXtCvlZMt0TZ7+n9mjDF1xqkf0Lhs8LnW5XngVYevbZW6O7b8zDSlRKHLw5lhUyfd5SQJcWl4SGOayijgDQ+MAeYXFO2F4UcrFnrxh4GBMUD9g/mR7/t+llUwnVIwNwetmZOMf2DQU+Gex+yFEF0/CbCMzRiTlUuMKqcwtYwCAN8nqg7fiJXmYtyi4TdoA0u2Gb1E26C+Jdva1Rmqs76p95vDomjkUnuF5VmiosGTJvuMqjOuvFXD+HR417hl2541V8MIjpM9zKlP6heoLk/XHbcdR9IZeB7MXCnTifMXaN/wwDiYqxDXhoOoIITnerSVAWZmgPJJ1jgwfKm2Qa6bZY11d2WDQgFjM8bKGV6dwnEMzxgDxfnhOuPyZYOzxT1hSBhkF+lOXIYoMnvSXb8gIAiz8pDWzAJzZs9Hy4uioe5dK4f6e7bc1UCxEaGdk5vsKIJg0ZaDM1NaKh/3+wBRycWrTFfw+LLixXKuDt6frU3dQO8bHxhDVq86mDU2PTAeXLqtJyiH1txZex40Gie1xoUAa0opeuYuBOh6dlEzOWMMw+UUxUZk7NJDPYN1xo6rpm5Y8KWEIUpl99CtxgK1mrld4ZAge/pduzpDseybvUTboEKBsOjmVgnU9Qa1GcNX3OgTlxRp8eRHq84VshnZFqnUXdLiyT7FSxUzl34cwamWKfaNajzvQ6BeJzsipDPyPJi5ENGuntxVH098NXjpr+TyDIUw/xMHVfPLKPrNzUFaqtCJYutKKaCbNf5gAa0UfugafVNTWirj+Cfnk9FlFD2lUq6ms7RUxilYEEH2JuDVi6RxYk8ZBUChgOMqvGXLssVdKi7mlqQrvz1rcrcxpFjkODB2HCi/ZVe2GLrLts2clFNYUUbRUyxSqmTXgdQvUFmavky4RafL2czPQ3sue+gHmF9jDIDrUn0nf0ULKnaUUfSEIVQq4Cwt4niO2b/XGHMrBXRj1vil6JSTX53CisA4CChWT34XK8oo4CQwvjRPEFiWkAsCaDRIGgXzl2gbJYqOn1zYKSY0VuxKhih18qCPcgWcWUvOuQG9ZduCUBHMWtBX9ihFvJB1KG6jNpXPGpPAuMv3YXbJp92YP962IcjqX7pNKfsyxpDd1BRmK9h5lcsOw7mvLaBC88tEemUGfkERNOwYw47ns/1wC64dM8chO+iKRYrLNbuyxZAFxgsLLC4aPQAzXt+joQsrc898QqOJijNFtFJUFmNr5swMqswFpIUwCyINHrkeJZxN8LzpW42ix8Zu4aXNz0NndgHt+dYExk4UUL7QvfMMQBXtC4yTJKs1ZmV4dRFbzCz6sDz+UdGmSJazcoriXMmaqCRolHDdrJNXrh37BMDFixQCZV9gHIYQBDZ076MVi4RhNkxdv2zJjdqAOHEoVGOKK/aVUfSUSqDLFeMf6jFSuUxcm77VKHos6sXPzvdhdt5BLy4Zu4bxKI0PskK6IFTWTU7rqdWw5vcaRSmsyIwoR1FarthRRtGTJBSLULls2UW6u2ixxaeVncKQIABnbpZ6w44JW4PiGKoXkqGHINlEKaheqlBcqk66KecvDKm8Nz+166JLlzdgfh729xrQeWzNFaEwUyaeiSjEypqZrcJc5Uu1XF2u8cKQ6lJx5NJtQrx2jgPFIgsXZ2wZlBlSKEDtK/Pg2VViMGjxvQQs3cX4rel9qo4dkd858n1YXlHQWYHNzUk359zUvzKH3tmddDOEyIYGLbtBK3314qSbIMSJ5WVqZbsv76MeGGQby0qL86b4GmD3mfOSsvR+eap/uBdVulgnXZ90K4TAqvPq2LSOCYo3U9miUiUhXjNLB1rOiU1rFDmOtcvaCCGEEEKcBwmM3yQ2ZuqEEEIIIc6JBMZCCCGEEEIggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgBveGCsU83ew212721NuimvRNrqcLS5P+lmvDKdZged6kk34/xonf15A7whuykM0D5s29WP9OkctUlbnUk345VI2yntw7a1nUnatHff+qWt6buOe5NuwCS0Djs8+NY9Dh5uoVttwmpIaaU66Wadm1vfuEnzyS7pUQuU4oP/5IdAqUk361w8+s4Dnl57fNzZL//Yu5RXyhNu1fnYvrfN/V+7DihQisq7cyz9weVJN+vcXP8X36N91CHV4HguX/krPzDpJp2b+799l+0ba90tzcqPvUuybMlxeX+He796HYUGNJUrMyz9kUuTbta5uf5zH6JTjeO7eKWQt3/qK5Nu0rl5+K277Nx5Akrhhj7Lf+Jt4tnipJt1LnYe7PDg164D4BQ8qu/NMf+DixNu1fm59nMfkbZTlOfiFgPe/QsfTLpJ52b1e4/Y/GwV3WoDsPhHr1C9XJtwq068kRljr+Cw//Dp8Y9y+PSQzu7BhFt1ftKjVhYUA2hNa/dosg06RzrVuQxI56g9wdacL52mZLGHhjSdurvos0o7Gp1qlNbW7ZtO9UnG37LEv9agdHqybxb9dmn75DxLWx10267satrs9o9a0zloolx7Lvmdw9bx67RpWdZfa9JOmr1sd+zL+mt9HH8BU7d/9pwlL0A5imip7+5Ew969zck16Jx5pTC33dw+nFBLzp9bcHPbneZ0nVBnoQajKceu07N/0MKmwBGyPiXHoh3U5PctyxzboX2Yv7F2CnYNoh4nSLq8yJ9QS87fYFLEDe3ZN91q03+aKd+u41J5030dt+vK+wKKy/m0/f69JxNqyfkbDIxbOxYFxkG+g5i2O80zSdPcplZ2nZ6Dcb5FseNwqZJNOzewb5ZUZQHDwZVtgXH//mmlhhILJhu8qfFCe36740x/l22B8eBxOG3XcbuuvC8gWSyR+oXj7f3NI9i3Y6Kan1gcGA9mjC0qpRgMjK2KQBjOqg7urtEsDozfpIzx4I230bQmbfZljF0P1564eCgbblPGeDAwtu2GzRk4z6athOmNDYyjCJz6Sdb46Ahaj+0opxgMjNt79tQYO/50D8GcyRtUSoHWdgXGg2wKjAd2xabAeGg43qbAuNMh7eseVcGewBGGfzurMsYDGdTB657pnMFSiilLcNl15X1BuTpj7CmnKJQK6L4opL1rT8Z4sPOzKjAejBRtC4zfpFIKiwxmjG1ic2Csm61cl2Jb1rF/8h3YVT8tpRSTZdeV9wUlC3GunGJvswl7exNs0fkoBArtB8fb7cM2tKfrjuxlDZ1QTTv2C7C/lGJg26qMsZRSGMnmwHiwb7Sp1ACG66e9wJ6sav+KDWB/YCylFFMkSaCTnGSN9/aATfPLKXwfdHBSTtFqAYd2ZI2n/U7zTAaCKesm37n5AMui2NHqwHiQ1YGxRcPx7YN8RtUJ7AqM+wN/5XlWDbDZPil0sDRk2q7jFh1KL65QAG+ufrzdbsPRI/MDY8ivTNHpQLpvZ2A8bbVJZ2J7KUUudrS7xtimNVWHMsYWDWQMZVUtCkCsDq7abTrtk3PMtmz4YKA4uLyZ6aY9wWXXlfcllOaKpH1lB7aUU9i6lrFyVO5uc9qGYM7E9lIKm5OqFu/cm1RKYdUEriN7M8a2rz89WEph2/65Lmhneq/jb3xgXC5Dp9xXTrEPPDF/Et7Qkm0WPf2u/+45bXXsiUEsX5VicHeMzhgPLdVwyk2M1tmwjaGGzi9rTrg3KzC2qX7a6mw49tcYKwX9awemrU6uX5l0F2PXlfclJAm0yyflFPt7oJ+YX07hl4Lctk0rU9j6kA+lLc8YW7SO8f52m+3f/RyePgXGPPlOazqr66z/6sekHXODSZsyxq2n+bXqB/sOmx+AYVXG+MCuNYw7B81cMDhU4mPBTc3u5vhVUjodIE3RGtZu7/Hod+++/gb2eeMDY9eFYiMiLWQZ1jSFw50W7O5OuGVnE8Qe2j058Np7R8drx967Z9gjsAduHwcL948vAO02O188fF2tOrN0/QlPPnt8Uo867sl3WrN3b9O458+kj9dz59HYaoPdXfZvPjIqqRqUfO5+qbn5y9fZ/p1PUbs7uc/TrW3W//0n3Ph3d3i67xuV0erkkzfPrDGedHbnRWxce8Ldf/s5zacHwEDm0Tt5AEZ62GT77tMJtPDltVr5gQmrMsYDfaNtq4l0Hq1x7d/c4OlaE8iXFmilcLzsWtBpdli9vpNNqDfMweo2n//Lqzz8aJ3mYZob+dUatlcPuP5vb7P+G58RhJNNCpl9NJ2TchnWyzWc9Syo2t2D6MkTKJWyL2idXeCTZIKtPF2awuPHUKtBEGQTC9NCiHuQBSatpmZ344g7dx30rdssfp8D1E7/R6fE9p1NWmtb1N6fwymXhh/ycdTm4M5THn/3AWmpTPLehBr6gpw44umvXmXj08c0vraM18l3/grN/s1HrH3ymG2d8JWfNOP36ml2XB7868+pX0qofGUxF1AprUl39mg+eMDj69s8LV/g+65Mrq0vynWzddCPbuxx/9o+qPxdy6Nr2/QSq5Wl8gRa+PLabbj1u+vMv1+lNusNPxK6u2Pb2/D42lPeWmnjzTcm0dQXFs3GbH62yt4vfUr9vRk67ZNzTvkeB6vbbH6xxs7dLao/+BYm/XL64SOu3otZfC+h0RifdUzTrFqwWgXPkAigvd/k4e98yfwfvEihHJ6eMU5T0pTjYNIEhXoJZ/sR9395m/VLCwS7fb+d69HaPuDB9TW2bz9BLS0y/870xiLjVN+qsv6dL9n63h22Prk/9Pm9X/kC1b3LjudLr7t5OYacFq9WksCjpIbfDYz3dmF2awsuXAClSG/dyeojpzgwdhw4XH3KF9/eIlyoUlwoZ0u2dQPjnR3Y+O1V3O1NCm4Htzg34RY/v2ixyv1v3mH92ibl+Qia+c/Xfus6+ztt0FC7PNkT6oVEEXHVZ+PREau/dXOoBnf/O5+x271wB+/V8A0bLQznyqAUD7/YYe3WDoNlctvf+ozHR4CG4vuVibTxLEorVY5u3Ms2hupwT14mK2btWxCAOjrk4S9/xNrCDM78bO7zvV3N6rd3aN25T6L28N59d0ItfXHxXNY/6FSz8dla7jN9cMjtr1/LXjsu9bfNuhEtRC7h3S94uF5n7a0VitsnwaN2PTqp4t7tNls3Ngj8lJk/sTjB1r4YLw44erLHzV+8Sv2DeQZTpl7o0dk9YPPaOk/vbnPxJ75q1nB4HBPH0NpKad5+kLvEqXaLrW9eBbLsce2t+uh/Y8r5oUu4UOXwwZORz1XoBcWeB0FDAuOJK5VAFSPSIMI5OuDgMBuGcnZ3OXzwhPXPNlj5/unvJGffLrP32Ze0P1vn6RcuTl8pRZqCt7UOQFgku/oZwg8coqUaB/c22Hp4MPT5/vbJSRbNxK+zaWeWLCVsPMomew7W3KbdoFgrh2TZpNxVl+tSWkw4urVNuwXtgeG/3tLaaSEkmTHneOypzgesdfuMcdyCS3Gm+BpbdT7i5So7D1fpPHxM+9Farphie/UAVr/ABZJFZ6oTBoO8yMcv+rT2Tx+LLi7X8AOjQiuIY4oRtLefoD9+yv7AkoGrv3UTd3cLV2uq/6Ehw2o9ShHUY1oPd9j45NHQ04J2rt5l43eO0BrCS/P4BcPmZ3gexVrA1tbpk+Q7pQr1ecMyJH0qVxpZYHyKsBpOfCjDsDP/1VAqC47bvdUpdLY6xcEnN/nyO+vZzY0BM4XikqKwMg+ASjs4rdEnWRAAYTjys2lVufLsoVqtHEpzZgUh0Xz5mX1AJ6lSqZl5qpYvPDtb2ilVKBsY9xcK4M9WT/1OvFg2chJleTE+nqOgTikiLi0mxu1fOPPsbFTl7ZnX0JJzFkUU4+y3UOnAOridNt7OJkprdBBSu2jOzUxP1Ojr2wcOyaPto+Na99q7Bv52QHHu2cdlsNgwKac1pHoh4VlDn8/zv8OrZubV9hx0Ot2ZrVevwhdfUH16B6d1MoCxsQ53brTpdLr9vgGBMUDj/ZncpLtRggCjMsaQnVAqKJz6nUItxvXMukiTJMel7OM4jRpFs+L9Y+FC9ZmHmj9ToXD6Tzu1kgvVUz8vLZtVRtFTShRpuXrqd8IICrPm7d+zRpVUMaK6bNbIEwBK5YPHMYoXZyadkHspzzUaWE6ozJuV9Onxq/Gp/aB2PSoXzTvf+jmuIr54epJLAuMJUgo+ve5zdXORe5/u0H60flxqAHBwcDLb2nEwZup1reHA/Pyp3wlChWmRiFKQvHX6CWVaGQUAhQKlmfEduXZcM8soegoFkvnxF2vtuBOfaHEW5YUiqT/mXFJQMmziXY/jQHGpeup3khJQMe9CXZw9vZ9I3mqYlgQ/FtSKpy5/rpWi+rYZEyUHFWfjoRKKQeXLM8b+dsQxxVMOzXalQb1h6s6dqFweXyPtuhDNTv568MYGxo6Tza07CGtsOg22tsZ/16SMsVJQe28291SZ/BcgSArGDX8Cz+zQp+GEehnxQjL25+iUKsaWUfScNvmsE5cpV8w7FnviGJz66PkHUb1o9PqqpaXy+H4ESGZD426wIRuSd9zRx5xWito7ZgaOAKoUE0WnfF6rUWkYmC4GVFAgjMe3XTsujberr69B5y2KiEvj+/pwZca4CdijlOcjVDw6WRKV/akYzTb7intG1Wr2pzl/IfdY6EEmBcYAswsunfrsyM8KPqhw8gfey4jrAV51fPBraubRqZaJx2QKdLVm0tymkU4rp+iUKsbvX7xcHfm+qdninkpV0YlH70OhAMGcedliAByHoDo6eizMVimWzQwcASgWiU6ppii9ZXBGFQgb41OqwVKDsGhwSKPU2NGMThRTWzSzRGTQaaO/05LcMvgoOh8XLmQPjGguX0aP6TEcB6MCY8+D8rvzJw+I6GNifXG/8uXRJ5QTRwRFQ59aNabOWDsu8WLZ/CdDF4sk9dGpjnC+0v9kUCNVluKRdf2mLdM2KAjAm6mO/CxJwMgZk13jyq6eZ5LvVAtDivHoDiP1A+qXzL4LPS0wNnLC5ACvEo+8PHeqM9Smf2Gs51a9Uh8Zb01DfTGcMTBWSt1WSn2klPquUurb3ffqSqmvK6Wudf+e6p+zUIDFRUijmNbMmHUdDcsYA8wtebSrwx1FEGJ0YFx/uzYy4D+tw5x6rks8N9x+k1ejGDRqkloniinXDc7OdSXl4YlqhcglqBt8THbFS5WRF7CkYtYybYNGZaa0X6B2ydxgHziZgDcix+MuzI4dmTLFuKUPdVyivnxKDYkpSqWhidZaORSXa0ZOmBwnrnh49XziwHHsyhj/Ka31D2qtf6S7/beBb2it3wW+0d2eavPzEEXQbizQiYZ/GEdhzOS7nijK1nMcvKiFBi7V1s8PXaLF6tD703JCvazCTHnoZ2mXaybObRopnK8M3Y/ZUEYBWYc+WE5RXDBzmbZBlbpLWsz/SJ4P0Zx5y7T1G3UjGl9s4Pnm7lOPk8RDfYnJk+76BfV45AhTcnnG/JE1yB70MXAp6yRV6rOGD6uNUB54UEkUO6jidNzcvIpD6S8CP9t9/bPAX3oF/z/OlVJw8WL2orn01tCEE9NKKXrmVgp0yvmDz/RSChg93GlqffGxgXIK7bhEc2UrJlsAkCRZlrGPTirPXKrOFOWV/EQ1o1cS6ZMk0CnnB/1MXY2in1sMKBTzKTgbAkcAikUG44tOUqMxb0HK0fOIqvnrl4lPKRzL8ygO7F+nNkO1OpnmvEq1y9Vcn5mtOjIdN6ZnDYw18MtKqd9TSv1M9715rfXD7utHwMi1w5RSP6OU+rZS6ttra2ujvvJalUrQaIAuBDQXLuY/7GWMDcsaVyrgLC0cb6veKm0GziTvN7Smse9TrJkd7FMqUSqfnI6dpEqlOh2dxLlwHMorJ8Gi9nxKc8Vp6QfPrFJVdEpZsOg43Qd7WMBxoLiYD4JNry/u6S+/UuWE8qzhfUhPHA8Nx0cX7FjRACCs53fOm6sTJzakizNu+STjn/oByVJi/DyMUYJQUVg4SdxNS30xnD0w/mNa6x8GfhL4m0qpP9H/odZaM/SMmuPP/oHW+ke01j8yOzt6BYXXbWUlm7jWqdRp92Vand7F27CssVIwsxLSTrK76SAkC4oNH3NyXEXp4snvE9Rj8wMspYhmS8d1ZO1yzbosQTB/sjqFLWUUPZ6Xrb4BUGxEz3wYjUnKDf+4xMx1oVgPjR91gnxgXL5i9moNOUFAGJ9EUjZMuus3OHGy8s50xA/nplQ6rgVvVxvUxy/7a7zeZHqlpqsc8kwRktb6fvfvx8C/BP4QsKqUWgTo/v34rI18XTwPlpez182Fi8cL9x93mIZljAFmZqAzm2WNQwvKKHr6hz2n6YQ6k3KZUikbGvRq5VPXIzVSpULSTTSa+hjo0yTLZbRSxi/TNqhSyUYwIMsWq4od+9erM9aOS/1KdbKNOWd+NT7OEOuZWatusvsD4zSKaaxY1lH2Zfx1rWF61dKp6hdidCGgGCtUaXpmhr50YKyUipVSSe818OeAj4FfAH66+7WfBn7+rI18nWZmsrIKXJfm0mUAjhdBMCxjDFmGp75SpBOXragv7inNhLjl7ESapiGYM+nWGVtXRtHjeZQXS2ilUBX7Av9qwyWNy8Yv0zYoDMFtVAEoJRhfX9wTNmIcB4LFutnr345SLFIsZpPukrcapg8S5rhJ8fgyZurjrU8VRRRLDmmpTGW2YNVvN8jzIFppENUjpqle5Cz/k88Dv6GU+h7wLeCXtNb/Bvg7wJ9VSl0D/kx32ygXL2ZZ4rRYojWzaGwpRc/cHLQaC1YFxpBNwtPKIZ49ZUV7kxSLxBWPjoVlFD3BfJVCPaFcta+3D8NsaPC0h9CYKpkJTp7MZcuMSdclqoVWrH87JI6JitApVe2YdNfPcQhrEVo51N+xsM5AKZwkprBodxlFT+VKY+qSWy99xmitbwI/MOL9DeDHz9KoSYuiLJhcXSVb27izDewZGxiHISRLCcFRyeil2gbV3q6zc3fLiiWWepxGjXK3pMJKlQrJBQfPjtH4IYtfrY1cQ9Z0lQqwVMWp7Bs/R6FfeGGW2rIlN9b9uitTuKVZK/uSsF5EtWMrb7ABqFQoJ1Xrys1Gqc4VULWRazRMjGW3kudnaQk2N6HZVHD5Mty/amSNcc/SEnitBeNXpOhXiFwaX1uedDPO19ISK6myKfbIC0OqV3yUpT2PLTP/ByUJ6IvV7JnyFpn5YBZnekZwz0+hQNAoMVOzZ9Jdv2gmplqxYNL1OLOzzCnH3v3r4zjAlE1WtvTydHaOkz0u+sYNUGGQbRiaMQa6xfwVo4P7UazL9nge09VFnD9jH939BnPd7rnWtuvodFyLI4+LF5m1Z4AwJ1ys0dAWhy+OY+PAkzEsPrLOrlrtrgXskM3Ka7cn3aSzexNuQYUQr4Z1M50sFkXWBlfKtz+BICZHerlnuHixL1EsFwUhhBBCCGtJpPcMFpXkCiGEEEKIU9g6xUcIIYQQQogXIoGxEEIIIYQQSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAUhgLIQQQgghBCCBsRBCCCGEEIAExkIIIYQQQgASGAshhBBCCAFIYCyEEEIIIQQggbEQQgghhBCABMZoDe32pFshhLCG1tkfIYQQxvEm3YBJONrc5/Hv3eVor03roE28XOXiH7806Wadm9u/8CHtjiLFQYUB7/7kO5Nu0rlZ+3iV7dtPSFsdOs0OS3/sCuXlZNLNOhe7j3Z59Nu3j2Oq5FKdhR9emmyjztHNf/05abONTjXKdXjnL3x10k06N+tfPOHpjXXSgyPSwyYLf/QKlUu1STfrXOzc3+bhb91GuQ7KdShfbjD3/fOTbta5+eyffg/HVTgFj0IScPFPvT3pJp2b+9+8w979LZyCh+O7LPzhSxQb0aSbdS52Huzw8Ju3To7Lt+rMfW1h0s06N9d/8TN0q41S4PguV37qg0k36dw8/vAR27c20KlGp5qFP3yJ8kp50s069kYGxo6TBSE9R3sWpYy1prnXotPpbrftyly1D1o0t/aPtztH9vx2aTultXt0vN0+tGffAFo7h6TNbJ8cz67BqvZBi+bGzvF2Z+/olG+bpdPs0DlsnWz3vTZd2tHoVptOq7tfaTrpJp2rzlH7+A+ATu25HqSt/HHZPrDnuARo7x2hW9nv5vp29ZedozatncPjbd2ZrvPOrv+1n5MfeTh9e97ca9kz8pmmqL59Sy37ib3AzW13mp0x3zSPUpNuwSvWt4M2XaAB3GKQ2+7s2xMYp+38RUu59vQpgzefbmBXrqh3I9pj0/7ZfFwOsfziMG3XA4uPpFP4PkH/dazd5siW61ia4vSdQ9pxrEqCOL69gfEga27WulTfgakt27nBwDg9sKVDGQ5ABs9Bk7UO8/2HW7Bn3+DNCoxtG4XqDxZti4v7rwUggfF0UIpCdNIBqk6bg4MJtuc8pWkuG46yKzD2wnzH3j6yJzC2rfMbMriDFgXHXmxvxli38+eYTZm54YyxZYFx6+S300pZNSRvc2A8OH/XtmvDUGA8ZZcCe46kFxTEJwGWSjscHkzZL/Oy0jR/ElmeMe7v+G0zbZ3FmfUfmBqrdtD1HbTnH2+nh01r9s/mAGRwxMmzKKMKHNeoAigvX0JousG6VJtu2Ia6Dtsi4wGSMZ4SQeLntg+2LSncH8gYa9syxoM1xhZNvhvs+xTT1Vmc1bQPn52F50Hqn2SNOx2g2Zxcg86RzaUUgyNOVpVSpCmdvsnXyrNo3xg+Lm3KhuuBxIFtcfG0XwvsOZJeUFDMZwYOdy0JsAYm39mWMR4MjCVjbBCbSyk80IWBwNiSiQuDpRRWZYyP7K3B1a127hRzCvbsGwxnjG06Lge7Rpsy/SClFFPLjzzcvhirud+euh/npVg++W4wo2PT5DvbsgJDBnr3acsSnIXrWhwYWxyADJdS2JNVHdw32wJjm1elGIxFBgNJ4ynJGE8nz6NQ6Ntutey4jlk++c7x3dw5ZXNgbMWNWp9Th88s2Nn+lSnaFgXGVpdSWLxc2+CKFDb9bmD3DVvWHZrfJ44zlASasv7fniPpRfk+QXiyaezKFO02HJ4slG375DtcN5fpl1KKKba5CevrJzsyKvLXGjY24ObN19++c+ZEfStTaHtWprC6lMLmjPFAmYjtGWObjsuhjLFtCePBDPiUXezsOZJelOcR9GWMjQ2MXRc++wy+/BJarfGT73Z24Pr1fBBtqFzmo9M5ecqf4WwbLntyELHzyR346CNYXR2aTKifbMLVq+x+fJv7G+GYf8UcucAYewLjTsveAMTqwLhp8cRC7A6MB5NZtl0bpr3G2K5byBfhebmHfKh2y8yYUSkolWBtLcu8RVFu8p1zuIf67FMI9iGK0EGITs0u5ncLLhxknb5KO3TaGte1q+OA6essXlShHHLncURx7YDGxj2c/fzn+9fu8+QJHBxA9Q+VJ9PIM9A6n8nxQo8j5aB0dlXrBcatFjx8CBcvTqKVL+7pUwgCCLv3KqcNWbfb2cRDU/VnVZWyK6s6WEqhfHv2DewOjLUGZfOqFIMJY6kxnhIDT78zNmMMHBUS1te7mZ29vdzkO29nC72XRSRP3TpXrw7fjZpmsFau3TR8h7ps6/ziGKhWOTyE+/fhcOD8un8/C4q1ckgWSxNp41ns7MCNG7C9nW0PrkzRbqWs3mvx8cfQ3DnKRnVM8PAhX/zmGl/e0bTbowPj/X24dSPl7ifbE2rky1n9fIvmXrY0p9b5UizPg1ydluFsL6UYOi4tqqF+00oppi0wtutMeRGeh+dl/WCnA3Syx0IPZoFMEMwk3N3Kyjmr1eFs8MEB3L4NT5ZqRFWzMzwwPEGmfdSBxJ5OsWewc9SpNmpITSlILtbYW314+heThLhkzn71JAnc+eyA648dCkmA60JaCHCOsjuA+/fgwDnCOdqinN6DdyoTbvHzKc1G+N++wc7GY64uLFM47H96msONz9sc3l3D21xj6QdmAXOy/e3DNjd/4WPq37dI9f15VOckeHQ9Ba5L6yhl49oTCvqI+vcvT7C1LyYdGAkcnH9hWymF/ZPvTpjU7z+XwVUppisufoMzxt3osJc1Vu0WWg+X4E7bDzZSsUi55qJ1Nt9pYyP/8fY27BKjCwFJMpkmnqfBzEB/ZqS9tWtOzbHW3PvV62zeybJu427ImjtH3PvNOzy5s/MaG3c+KgtR7sEXoxQXysbdjEL2e9VnHKIbH6M++5Tm3VX665i0hsKDWxQefUkUpCe1CVPOrSaEITjNQ/wvb+SfnqZT2t/9CH/9IU7apvxWfYItfXHRTIzupGx8eJ9bv/gJ9AXGqVbc+c17XP+5D9n87h3i2eIEW/ri2mubXP+NR+ztZhetoVUp+jLGu087tA1bun/wWjxUSmFwOZ3W3adl9m33G+ofDfvx0laHrbsn169xGWOtYfPWFnvrkx2+f3MDY5VlBwrHgXF2oPUC43Y7GyY1peygvHx6xNtJagBGBsbNjR12bjymF/EOPf2u2UEfNVn/1k2u/8pdc4IspfBDl0e/fo1r/+9nPL2XD3zTwyb3v3mHG7/wCdu3n1BaMK/coFwGXa2d+p1k2ZyM46CZ5YB2UsM53Kfw+B7e9pPc504ru9hFIcYExrguUWN8UNirfSzNFXGKhuxTVzwToZ2s/0gPjnJ1nM3DlP1bq9DpUExcgjkzMvw9hUqEvn+f27/0CXc+2uZof3gd4427+1z7t3e4/+s3jRs5fPAAHj06uSb3Z4yVo46jx1YrSxCZRCm48RsP2fjWDfT+wVBpQe+a1tzc4+GvX2f39voEWvnyHN9l4/duc/0XP2Pz9tPhL2jN1p2n3PjFT3n0zZuEldOTKa+aYafG2fVOKscBfJ8w6E7i0imkKQcHDp4Ht25l3zWl5CycK1MobI19Cm27Uj+ep2cav1Lk9tevsfr7D6hcaTCYEt67ucr2d3Y4Okjxl+eMmliYXKjy5PoT2k/3WP9oL/dZa2ObVjf77zfKBJFBO9alFMTLVQ7XHo38PPULlOfMCq76BQEEF+fpfDL+SlwodAeoTAmMgeJ8wub9/VO/Y1q2GLLfSyWlbIbhKSqX6+bNUA5DiiWX9uYR+9+7hlaK/hzB6jdvoA6y37T6A5cm08YziJ8+4PajkLW1OsvLA4Gx63BwAKur8GStw4X6HtTMuuEuzJR5/PFNntzaIlqo5j5rH3W49yvX2HmwjaNg4Q+YU+LTU1qp8uSzxzz6jetDE0G3rz9m+4vsGhHMlif+eG/DzvzzcfVqtnLZ5k7+6Xeq02Z9Hb74AlpNnX8AyLRLEspjEhydqASeTxybE+j3U55LcqFK67DD+tXHrF/fyn2+//ApRwdZJxnOmBX5FxfKeP6zU9zxklnZq36VpZjUH30yFRplfP81N+icNS7G2Tk2RhR1XxgUGMfzJTjlsHScrH7cRGEjPvVzx4Hy5cZras35imZO9k0NjMf3gmKtHKqXzfvtSvUCwYNbOJ9/yp2Pd0j7AuOOdvj0e022P7lL+MWHVEPzlpgqLZXRStFuwc7drdxnu1ttdu5vg4a4HqCK0eh/ZIqVL1aPX/eXZwG5ofl4cfI3NG9WYKw1jgNzc1nC4N6qx/37Jx+rdotWC9ThPoX7N3OrVky9MKRcHx1hdCpZZsfEMoqe6tvPd6GK580KjHHd5yqRSFYMDowrkCbVkZ+Vls3dr55aDfTc/NjPo4gsbWxQBrJXZzxOaaGECkzKHJyIZk8/38pzIU5yevA8rZ61bwCFuSpRybwMiVdLCAJwDvcJv/wiH/h32oQ3PsbffEwcpfhz5gX+5ZpLesoNdo+p14JotoQfPbtIYRr2z5ye+jzs7MAXXzDjbuJ7Gu3lA0nVPMJ/dJfo1qcAZmWMyVanGLyYaaVoG1xf3BPNlwlLp59UqR8QV81LPyYXqqd+rqMipbphB2Mf14VoafhCpZUiWTL4oOxynCwbMm6SYRhhVLYYANeleEqdsYllFD3F2Rh9ykSE570Jn0bRTPzM+6/ylZnX05jzFgQUq6P7QaX1caCcLCeYOAzl++DPnB4UKvXs68XUUuqZI59OWKBYn3xf+WYFxuUyNJs4t2+yvPEh7l5+slPw4Bb+5mMA0rBoVsYYIEmGgt+0mIDnZbWeZiZBMko9c3jTrZSMu5mBrEzitItZcbFizoTCMarL8fCNaBwTl83LXI0yMwPtxnDWWCkIA8wLjMnqjEdxPUXpgnkZuZ5iyUEHo4eiwwiiFXMDY+KY4imLaaR+gdpFc29Gxx2TxxSUL5l70/aswLFUcXHKho2K9ilfOr3fiBYmny2GNy0wBmhknV4taeN3xtchpWFsXpCVJJQHynPa5exALJWMGskdqXqlfmrdo2n1xT0qDIgb4wMnG8oNKlVFe6CcwtRl2kaJY/DmG2g3P6oRRd0Z5QYGxvFCMvJ8S5bLRj9FzfPAq47uKyorZSOzjcdcl6g+vv40Wm4YvXvPqj8tlRTeTPX1NOYVqMyHpy5vWVqumPeghT7xQnLqxLppKKOANzEwnpkBpXAcqJ9yY5mGRfMC4yCgUCpkQ7dkQ9UmL9M2yE2KJLPjO/3inJmBMWQzdkfRrkd50eRUf8b3GZppXVqa/CSL8zQz59Cqz+XeM3HiXY9TLmXLzA1IDM7I9YyagKcU2ao3hjutHzS5TASeXfueLJeNfoJVHIMeN4ve5DKKHschHpPomabSujcvMPb9bDYQWWA8apWG1A/A88wrpQAol6l04400Lh/v4GAm2VTjLlza9U7Nuk67ZKUyMjvnz1TMu0Ebo7ycHK8hq12P8oJZD1B4lkYDOrVZdN9DPo4D48i8WeS4LsWZ/G/kFZyxN3EmGTVJrVxzcRvV19+YcxbNxCOTirqUUJ038aLWx/eJaqP7eaXMrn2HbFR3XDlBqaRwDVuCbpTymOC+UE/wCtMRkk5HK163mWzywbiscRrFuK6ZS5uRJCTdc6ddznbOcTi17swkpYt1vBFDgTouGb2PbjkmLg9nOuLl6utvzCtSrSk63XKKQiPBL5g7JDiK50Gl4dGunty8hRFZR2JoFiteyGdwSssV82uygLhWGFpC0Mi1i0dQSWnkfVjpUsOG3RtbZ1wqOzi16RiKP4tkKcndXPeUFhNDg5K80nJl5FMKp2kE0YLT5CWUT+rIarXhvjCNYjOzxQBJgu9BFDvHQUiSGF2WlKMKPuURT0oLGiWz91GpoVpirdQzn2hokiDgeBklk9dlPs3MDLTr2SQ8zwffw8xscVdxPsmdV5XLZmfkeopF0NFJOUWhAPFFs8sMjgUBUZK/EdOOS+2KuRMm+42rfS9fqFgROJYrik5p4Bqn8usAm0x57tANN0zXXJo3MzBW6jhr7LpQH+gP08jAiXc9vg9hmA3NdyN+G+qL+42qk3ue9Tun3VD9WCmxZtWGnl45hcmPgT5NuQxeHNBOqhQNri/uUUmJYpxFIYXIpbg4PRevs1Aqu5nuqcyHhi/bk1ecze+Lqtcole243I+qM3YcjF4ppV8YgjuQ+Y6L4NbtOPdg+FqngsJUlULacaa8jJmTtRzrfVljrRRpEJkbGAMkCZXL9eNMj22BcTCfX6BeK4d41uA6ii6/USYqnqRCbFimbVC17pDOLRi53vTzUCqrNW43FoyeeHfMdYm66xmXLtTsGXqi72ZamT8pbVBxLv/kwsqVhj0/nedRHFh5I6nYUUbRMziilswXzXuwwimSC9Xc8RhN2bXuzQ2MC4XjGWmue1JrrIMIHMfcUgqAahWvXibpliSZXHs7kuNQeeskO5CGReLSFJ1VL8txsjqyrmlZuuY8FYtQfX9+qjrB8zYzk406HS8faHJgzEmdsS1lFD3FRoRWDkkJvDm79s1J4uPDLi2E1C6YP6LWb3AoPrlQtaI+vKfc8EnDkwu38atRDHACn3juZFSjNAWPge6n9MDz1CfSCKXWgDuTbod4pWaA9Uk3Qkw9OU7E85DjRDwvOVbEKJe01rOjPpiKwFjYTyn1ba31j0y6HWK6yXEinoccJ+J5ybEiXpQ9Yw9CCCGEEEKcgQTGQgghhBBCIIGxeH3+waQbIIwgx4l4HnKciOclx4p4IVJjLIQQQgghBJIxFkIIIYQQApDAWAghhBBCCEACY3EOlFL/UCn1WCn1cd97daXU15VS17p/17rvK6XU31NKXVdKfaiU+uHJtVy8TkqpC0qpX1FKXVVKfaKU+lvd9+VYETlKqVAp9S2l1Pe6x8r/0H3/slLqd7rHxD9VShW67wfd7evdz9+a6A6I10op5Sqlfl8p9YvdbTlOxEuTwFich38E/MTAe38b+IbW+l3gG91tgJ8E3u3++Rng77+mNorJawP/ldb6q8CPAn9TKfVV5FgRw46AP621/gHgB4GfUEr9KPA/AX9Xa/0OsAn8je73/waw2X3/73a/J94cfwv4tG9bjhPx0iQwFmemtf414MnA238R+Nnu658F/lLf+/9YZ34bqCqlFl9LQ8VEaa0faq2/0329Q3YhW0aOFTGg+5vvdjf97h8N/Gngn3ffHzxWesfQPwd+XCmbHzwuepRSK8BPAf9Hd1shx4k4AwmMxasyr7V+2H39CJjvvl4G7vZ97173PfEG6Q5h/hDwO8ixIkboDo9/F3gMfB24AWxprdvdr/QfD8fHSvfzp0DjtTZYTMr/Avw3QNrdbiDHiTgDCYzFK6ezNQFlXUABgFKqBPwL4L/UWm/3fybHiujRWne01j8IrAB/CPjKZFskpo1S6s8Dj7XWvzfptgh7SGAsXpXV3rB39+/H3ffvAxf6vrfSfU+8AZRSPllQ/H9rrX+u+7YcK2IsrfUW8CvAHyErp/G6H/UfD8fHSvfzCrDxelsqJuCPAn9BKXUb+H/ISij+V+Q4EWcggbF4VX4B+Onu658Gfr7v/b/eXXHgR4GnfcPowmLdWr7/E/hUa/0/930kx4rIUUrNKqWq3dcR8GfJatJ/BfjL3a8NHiu9Y+gvA/9Oy9OrrKe1/m+11ita67eAv0r2u/+nyHEizkCefCfOTCn1T4AfA2aAVeC/B/4V8M+Ai8Ad4D/WWj/pBkf/G9kqFvvAf661/vYEmi1eM6XUHwN+HfiIk3rA/46szliOFXFMKfU1sklSLlkC559prf9HpdQVssxgHfh94D/TWh8ppULg/yKrW38C/FWt9c3JtF5MglLqx4D/Wmv95+U4EWchgbEQQgghhBBIKYUQQgghhBCABMZCCCGEEEIAEhgLIYQQQggBSGAshBBCCCEEIIGxEEIIIYQQgATGQgghhBBCABIYCyGEEEIIAcD/D+d3KXpLvdNgAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(12,12))\n", "plt.quiver(x,y,u,v,scale=100,color='b',alpha=0.2)\n", "plt.quiver(x1,y1,u1,v1,scale=100,color='r',alpha=0.2)" ] }, { "cell_type": "code", "execution_count": null, "id": "ba6fdd7d", "metadata": { "papermill": { "duration": null, "end_time": null, "exception": null, "start_time": null, "status": "pending" }, "tags": [] }, "outputs": [], "source": [] } ], "metadata": { "jupytext": { "formats": "ipynb,py:percent" }, "kernelspec": { "display_name": "Python 3.8.12 ('echopiv')", "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.8.12" }, "papermill": { "default_parameters": {}, "duration": 3.115691, "end_time": "2022-10-20T19:14:11.961668", "environment_variables": {}, "exception": true, "input_path": "./notebooks/extended_search_area_vectorized.ipynb", "output_path": "./notebooks/extended_search_area_vectorized.ipynb", "parameters": {}, "start_time": "2022-10-20T19:14:08.845977", "version": "2.4.0" }, "vscode": { "interpreter": { "hash": "f83b0c3a4910470a1212112b1707d582432916ed4ba8aec962241a050aa18fae" } } }, "nbformat": 4, "nbformat_minor": 5 }