Skip to content

Instantly share code, notes, and snippets.

@alexlib
Created March 28, 2023 20:20
Show Gist options
  • Save alexlib/e6efee07680a56c15a602ca1761c1a29 to your computer and use it in GitHub Desktop.
Save alexlib/e6efee07680a56c15a602ca1761c1a29 to your computer and use it in GitHub Desktop.
pivpy_plot_vorticity.ipynb
Display the source blob
Display the rendered blob
Raw
{
"nbformat": 4,
"nbformat_minor": 0,
"metadata": {
"colab": {
"provenance": [],
"authorship_tag": "ABX9TyMuOF9IYECgNSds60tqXJhN",
"include_colab_link": true
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3"
},
"language_info": {
"name": "python"
}
},
"cells": [
{
"cell_type": "markdown",
"metadata": {
"id": "view-in-github",
"colab_type": "text"
},
"source": [
"<a href=\"https://colab.research.google.com/gist/alexlib/e6efee07680a56c15a602ca1761c1a29/pivpy_plot_vorticity.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {
"id": "NQ946M7r1Ecl"
},
"outputs": [],
"source": [
"# PIVPY example"
]
},
{
"cell_type": "code",
"source": [
"%pip install git+https://github.com/alexlib/pivpy"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/"
},
"id": "KUILWy5r1Hxg",
"outputId": "d8e37634-abfa-4c3a-cb80-2b119b60f27b"
},
"execution_count": 48,
"outputs": [
{
"output_type": "stream",
"name": "stdout",
"text": [
"Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n",
"Collecting git+https://github.com/alexlib/pivpy\n",
" Cloning https://github.com/alexlib/pivpy to /tmp/pip-req-build-1ijqkqy_\n",
" Running command git clone --filter=blob:none --quiet https://github.com/alexlib/pivpy /tmp/pip-req-build-1ijqkqy_\n",
" Resolved https://github.com/alexlib/pivpy to commit b8e1d52beb8fe4b03d37999bed1c6889fe5cadb5\n",
" Installing build dependencies ... \u001b[?25l\u001b[?25hdone\n",
" Getting requirements to build wheel ... \u001b[?25l\u001b[?25hdone\n",
" Preparing metadata (pyproject.toml) ... \u001b[?25l\u001b[?25hdone\n",
"Requirement already satisfied: numpy in /usr/local/lib/python3.9/dist-packages (from pivpy==0.0.18b0) (1.22.4)\n",
"Requirement already satisfied: xarray in /usr/local/lib/python3.9/dist-packages (from pivpy==0.0.18b0) (2022.12.0)\n",
"Requirement already satisfied: scipy in /usr/local/lib/python3.9/dist-packages (from pivpy==0.0.18b0) (1.10.1)\n",
"Requirement already satisfied: matplotlib in /usr/local/lib/python3.9/dist-packages (from pivpy==0.0.18b0) (3.7.1)\n",
"Requirement already satisfied: pytest in /usr/local/lib/python3.9/dist-packages (from pivpy==0.0.18b0) (7.2.2)\n",
"Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (4.39.2)\n",
"Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (0.11.0)\n",
"Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (23.0)\n",
"Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (3.0.9)\n",
"Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (2.8.2)\n",
"Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (1.0.7)\n",
"Requirement already satisfied: importlib-resources>=3.2.0 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (5.12.0)\n",
"Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (8.4.0)\n",
"Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.9/dist-packages (from matplotlib->pivpy==0.0.18b0) (1.4.4)\n",
"Requirement already satisfied: exceptiongroup>=1.0.0rc8 in /usr/local/lib/python3.9/dist-packages (from pytest->pivpy==0.0.18b0) (1.1.1)\n",
"Requirement already satisfied: pluggy<2.0,>=0.12 in /usr/local/lib/python3.9/dist-packages (from pytest->pivpy==0.0.18b0) (1.0.0)\n",
"Requirement already satisfied: tomli>=1.0.0 in /usr/local/lib/python3.9/dist-packages (from pytest->pivpy==0.0.18b0) (2.0.1)\n",
"Requirement already satisfied: attrs>=19.2.0 in /usr/local/lib/python3.9/dist-packages (from pytest->pivpy==0.0.18b0) (22.2.0)\n",
"Requirement already satisfied: iniconfig in /usr/local/lib/python3.9/dist-packages (from pytest->pivpy==0.0.18b0) (2.0.0)\n",
"Requirement already satisfied: pandas>=1.3 in /usr/local/lib/python3.9/dist-packages (from xarray->pivpy==0.0.18b0) (1.4.4)\n",
"Requirement already satisfied: zipp>=3.1.0 in /usr/local/lib/python3.9/dist-packages (from importlib-resources>=3.2.0->matplotlib->pivpy==0.0.18b0) (3.15.0)\n",
"Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.9/dist-packages (from pandas>=1.3->xarray->pivpy==0.0.18b0) (2022.7.1)\n",
"Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.9/dist-packages (from python-dateutil>=2.7->matplotlib->pivpy==0.0.18b0) (1.16.0)\n"
]
}
]
},
{
"cell_type": "code",
"source": [
"from pivpy import io, pivpy, graphics"
],
"metadata": {
"id": "oUMss0hu1KM-"
},
"execution_count": 49,
"outputs": []
},
{
"cell_type": "code",
"source": [
"d = io.load_openpiv_txt('test.txt')"
],
"metadata": {
"id": "ysFZWvHD1mf6"
},
"execution_count": 50,
"outputs": []
},
{
"cell_type": "code",
"source": [
"d = d.fillna(0.)"
],
"metadata": {
"id": "WKs-Pbeq1tga"
},
"execution_count": 51,
"outputs": []
},
{
"cell_type": "code",
"source": [
"d = d.piv.vorticity()"
],
"metadata": {
"id": "hwpZJMa911Il"
},
"execution_count": 52,
"outputs": []
},
{
"cell_type": "code",
"source": [
"graphics.quiver(d,arrScale=5)"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 387
},
"id": "8NcBy-7o13Yt",
"outputId": "4f466469-7936-447f-934b-bcf8696d0ac2"
},
"execution_count": 53,
"outputs": [
{
"output_type": "stream",
"name": "stderr",
"text": [
"/usr/local/lib/python3.9/dist-packages/pivpy/graphics.py:327: UserWarning: Warning: function for a single frame, using the first frame, supply data.isel(t=N)\n",
" os.chdir(p)\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": [
"(<Figure size 432x288 with 1 Axes>,\n",
" <Axes: title={'center': 't = 0 [frame]'}, xlabel='x (pix)', ylabel='y (pix)'>)"
]
},
"metadata": {},
"execution_count": 53
},
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
],
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUEAAAEWCAYAAAD4hSV+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAABhCUlEQVR4nO2dd5gcV5X239Pdk4NGOVmy5JxtbAE2JhjbxCV9uywLyy6wGLyEXXJeWFhYWGCNwbAmGBsDtuQs25KTbCtYOYeZ0eScZ3o65+6qOt8f3bJnps5tq0cTe+7veUaaOdX31q2u6tO36pz7HmJmaDQazVzFMd0D0Gg0mulEO0GNRjOn0U5Qo9HMabQT1Gg0cxrtBDUazZxGO0GNRjOn0U5Qo9HMOojoy0R0kohqiegBIioes72IiB4iohYiOkBEa1R9aSeomVUQ0RoiYiIKE9GtI+yfJaLBjH3hNI/x5sw4LCK6eTrHko8Q0UoAXwCwjpkvA+AE8OExL7sFgI+ZzwPwSwA/U/WnnaDmtCGijsn6UBPRPxJRJxFFiOgJIlrwKk2qmPmuTNsCALcDeDszlzOzZzLGeLow84vMXA6gazrHkee4AJQQkQtAKYC+MdvfD+Avmd8fBXATEZHUkXaCmmmHiC4F8AcA/wxgKYAogN/m0MVSAMUATir6d53pGDUzB2buBXAb0l8y/QACzPz8mJetBNCdeb0BIABAvEPQTlBzWhDRfQBWA9icudX7xgR2/1EAm5l5JzOHAXwPwN8SUcVpjOsCAI2ZP/1EtC1jZyL6PBE1A2jO2O4gom4iChLRESJ604h+fkBEjxDR/UQUIqIaIrqAiL5NREOZdm8f8fp5RHQPEfUTUS8R/TcROSfuLdGoIKL5SM/01gJYAaCMiP5pvP1pJ6g5LZj5n5H+5n1v5pbz52NfQ0Sricif5ecfFd1fCuDEiH21AkgCuOA0xtWUaQ+kb5FvHLH5AwBeD+CSzN+HAFwFYAGADQAeGfNA/b0A7gMwH8AxAFuQ/oysBPBDpGerp/gzAAPAeQBeA+DtAD71auPVTAg3A2hnZjczpwBsBPCGMa/pBbAKePlOYB4A8TGJdoKaCYOZu5i5KsvPBkXTcqRvV0YSAPCqM8FX4X+Y2cvMscz47mdmDzMbzPwLAEUALhzx+l3MvCVz+/QIgMUAfpr5oD0IYA0RVRHRUgDvBvAlZo4w8xDSD9/HPpzXTA5dAK4lotLMc76bANSPec0mAB/P/P5BANtYoRajn5VoZgJhAJVjbJUAQmfYb/fIP4joa0hHDVcA4Mw+Fo14yeCI32MAhpnZHPE3kHbYKwAUAOgf8azdMXZ/msmBmQ8Q0aMAjiI9Gz8G4C4i+iGAw8y8CcA9AO4johYAXmT5gtJOUJMLWXXXiGg1gLosL/lXZl4v2E8CuHJEP+cgPUtrGs8gR/DyeDPP/76B9KzhJDNbROQDIEYMX4VuAAkAizKzRs0Uw8zfB/D9Meb/HLE9DuDvT6cv7QQ1uTAI4BzVRmbuQnqmlCvrAezLOKqjSD9/28jMZzoTHEkF0rMGNwAXEX0L9tnnacHM/UT0PIBfENH3kJ7JrgVwFjO/NFEDnmsQ0ZUA/MzcOZX71c8ENbnwPwC+mwlyfG2iOmXmkwA+g7QzHELaYX1uovrPsAXAc0jPLjsBxHFmt68fA1CI9MzXh3Qu2vIzHOOchplPTLUDBADSytKa2QQRnY10SkwcwNeZ+Y/TPCQbRHQTgMeQvqV/NzNvn+YhzQmI6G4AtzNztkcy9nbaCWo0mrmMvh3WaDSzisz68QYiWk9E9UT0aCZdZgcRrSOis4momYgWEZGDiHaNTHQfi3aCGo1mNnIhgN8y88UAghjxDDnzXPFnAH4H4KsA6oRldS8zK26HFy1axGvWrJnuYWhmAX1DfiyoKkNxYcF0D0UzSSSTSTQ0NOCKK64AAIRCIQwODsI0TaxatQqlpaUAgObmZsTjcVx66aU4duzYMDMvlvqbFSkya9asweHDh6d7GJpZwPfv3Ih/es/1OP/spdM9FM0k0dHRgbe85S0v+4Rt27bhN7/5DXw+H2677TasW7cO0WgUr33taxGPx/H0009jxYoVyqizvh3W5BUz/75GMxF0dXVh3759AIANGzbgjW9846jt3/zmN/HRj34UP/zhD/HpT386a1/aCWryDlk1TpNPXHjhhbjzzjtx8cUXw+fz4bOf/ezL21566SUcOnToZUdYWFgIKGS0gFlyO6zRnC6z4Rm35sxxuVy4//77R9l27Njx8u/79+9/+feNGzeCiJRCu5M6E8wobjyaCWfXE9F1RLSAiF7IhLBfyGiDaTQThkJAWKMRmezb4TsAPMfMFyG9QL4ewLcAbGXm8wFszfyt0UwIeiKY/6xZswa1tbUT1t+kOUEimgfgzUhL2oCZk8zsx2jt/78gLXyp0UwYeiKoyYXJnAmuRVqx414iOkZEdxNRGYClzNyfec0A0vUhbBDRrUR0mIgOu93uSRymJp/QzwQ1uTKZTtAF4GoAv2Pm1wCIYMytb0bpVbxqmfkuZl7HzOsWLxZzHJVYppmTnZmVH55c+wKA0IDstMfTl2qbmUqJ9mzH4mlqE+1GMqncv2nIcnkx31gh6Ff6Um3zNLeL9kQ4Iu87lUKge2wRsTSq9yUZjiifCaqORXleLEu0Z22T5RpT9TeR10W2Y5moa3wmH8t4mEwn2AOgh5kPZP5+FGmnOEhEywEg8//QRO7Usiw07zwgbmvasV+0ezp7MNwmV0dUtWnZfUg8SaH+IbQ8L0vKNW7fJ9p7axoQ9vhsdmZW7r/vuCyU4evug7ulw2YfbmzFyUeelvd/TCzShngwBF9Hjziu4xseF9sc+sN6OFz2ekPhwWE0PLFFbFP/1FbR3rXnsNJxNiveS09rJyJur7it/8RYBfY0TYq++k82ITg0LLdRnBeV3d83iCHFl5CqTdu+o0glEjZ7MhZH+4HjYhvVNTbQ0IrgoPzlnOuxBAfcGGhozalN+4HjSMbiNnsqHkfb/mM59TXY2IZA/8S5jUldNkdEuwB8ipkbiegHAMoymzzM/NOMsOUCZs5auWzdunU8W1aMmKkUkuEoSubPm+6hjCLq8cGIxVF51sRI3pmGAafLnmEVD4RQPM9eGiTmCyAyNIxFF5572vsI9Q/BWeBC6aJXK0H8Ct/46X3413+4CeeuXXHabTT5DxEdYeZ10rbJzhP8dwDriagQQBuAf0F69vkwEd2CtLjlhyZ5DFOKs6BgxjlAAChdOLGZSJIDBCA6QAAomT8v5/elYvmSnMflKimGs6gw53aaucukOkFmPg5A8r43TeZ+NXMXHRbR5IpeNqfJO3SGjCYXtBPUTBmq58+RoWFEFA/tVVHIpCKirFNkNLmSd07QsizUvbBL3HZyixy1HW7vxqAiCqlq07B1jxgdjgVDaFdEu3qPyVnuXUdqEB62RzQT4Qh8Xb1im77qevED7+nowWCjHLlTRdSatu0VU268bV0YUkQB+2saRHv30VqE3PIyza79R0X73h/8AsULqmz20IAb7vpmsY0q0g3Iy+ZMw8BgfYv4+nrFOe45UaeMqNY9v1Mel6IvX08/+uoUx6Jo07hjnzI63LzrYE599dY0wN83KG5TtVHZ/X2D6K1tzKlN884Dyuhw00tyFFjVV9/JJvh6B8Rt42FWiKrOpuhwNhLhCIrKy179hTOERDiCwrJSpVNRBUckLNNEsHcQVatHR22HTtRhy6e/jn8+aHdq4cFhsGXZAiRsWdh72x9w/Tc+a2vzrV8+jM99+CasXj5aNMQyTRiJJApLS057zJr8YTqjw5oRzCYHCGQfby4OEAAcTqfNAQLAossuxAcev0dsU750kWgnh0N0gC9vF5y2w+nUDlAjkne3w5rZhcPpRMXKZRPW32y4s9HMLLQT1OQdWkBBkwvaCWqmDCMhr1E2s61dzrJNgrOs99VoJPLOCTIz3Ip1wEPCmloAiPgCiHj9ObVxt3WJ6RvMDK9i/9KaXgBIxeIwhCggAEQV40rF5ddbpomYP5jT/oN9g2J0mJmVYhDKvgbcSESi4ramJ+W1w82bXhDtMa8fXYooaESIpgOAs12O8gOAp7VDtrfIbQIDQ4grUnHcrXLdHtX1EguGlFFzVZvhjm4xA8FMpeDptK/pztZXcGgY8VA4pzYqezwUVq6pVrXxdPaIAhamYYzrWGLBkLhtPOSdE8yGSl0k2+2Tug2J22K+ANwNciqGakeeti44C+1LvYxkEkFFWouhcIL+rj71bEix/55Dx+EssJeoHDzZBGeBIgCS5U1TvWfH7lov2lue3Y5Qnz3l4cQ9D4AFJ2CZJvwKdRlXayvCPfZtZiqFrr1HxDZtCtEB1TnORrbXZ7uWcuory7iyjnei9p+FqeprItXDdYqMZkpgZnRs3Y01N14Pcoz+7u0/fAILLjgXRZXlo+zu2kaUr1iKEiGHUMU3fnQv/v1T78Wq5XJkWTM30SkymmmHiLD25jeJ25avu1K0L77swpz3wyUlcDrtUl4ajYo5dTusyX9YSyhockQ7Qc20M9GPZHS1OU0uaCc4CSSjsQnrayIdhKqvbPtQbVNJ8pvJJAKK9c5JRdS4/+Bx5f5TwnpTQF0qgOPy6wF1SYD4BEYaNbOPvHOC0y2gEA+G0KxYXK9aqB/oGxSdDTPD3Sh/2FV9RX0BcaE6M6O/WhY9cDfJxx4aGoav0+7QLMvCiQc2iW2O/P5+MaIbHvKgYfOLYpsT92wQ7f0n6sRIu5lK4fifHhLbOBsakBK+hBKhMBqekvdf89Bm0a4FFOaGgMLLhXlm8s8111zDswUjmeSIxzfdw5h0jGRStB/5w/3itlQ8zjF/UOznyY98TuwrEY5wKh632UP9Q1x932Nimy9/7w/c0+u22ZPRGA83t4tthltkuyZ/AHCYFf5Fp8ho8oqv/HwDvnnL32DpwplX4kAzfWRLkcm722HN3GY2fKlrZhbaCWryC2YdHdbkhHaCmqxkm1nlOuuayFmaaq11ebushA0A4Q71Ns3cJe+coGVZaNi2V9xW/+Ju0T7c0a1cEK9q07hjnxgdtixLKT2v6qvneJ0or8/MaNi6Rx6zanF7MISwYqF+33G5yHp/TYMooJAIhdG52x6FNFMpHLlrPVJCyountRMJxUL9zj3yc11V2QEzlUIsIItBNCkizRXNTYgK65Bjw14MvyBHtEOKCHBvTYNSKEB1LlV2X+8ABhRrylVtmncdVEaHWxXroFWZEX0nm5TlFXI9lkD/kDLSrWrTsuewMjrcsvtQTn3117coI93jQhUxmUk/uUaHTdPMyW5ZFluWNSF9MTMbqdSE9aXa5u3oEe3hYa8YUbUsi3uPnxTbeFo7RXvHroPs6+y12dt37OPGzS+KbfprG0R7Mhbn5ud3itu6Dh4X33/TNMVIs2kY/OJ3fir29bMbP8p7f/Zbm33fT+7g4bt/ILaJhyOifTznZaqusfHsf6L6ytbfdB+LCujosGau8I2v/RIfu+5CXPZ37x5ljw0Owtj/OMrfe6tNwEGT/+josGbOkFi0GMtutAs1FC9ZAmfFfO0ANTYmVUWGiDoAhACYAAxmXkdECwA8BGANgA4AH2Jm32SOQzN3UN7ZmAag1WU0AlPxtfhWZr5qxFT0WwC2MvP5ALZm/tbMMiL96gfTcZ8/p77G80hGCuQAAKVSYooMp5KggqKc96PJf6bj3uD9AP6S+f0vAD4w0TuI+uWF8ip7KpFASrHwPte+LMtSSn+rRAfi4YgoPZ5tP/GAvA/LspTR2YETdaI9GY2phQoUUejG9RtF+9DRGnhPNonbVKUC+hQCCpZlIeqRbxJan9km2svq62Rl7VRC6QQjCmGFRCSqdLa5XhdGMqkU1si1L2ZWRs1VbSbyWMxUKudjiQWCyi+7XI8lGY0pP0vjYbKdIAN4noiOENGtGdtSZu7P/D4AYKnUkIhuJaLDRHTY7ZZTGCQsy4K7RU53GVQIBUQ8foSH5Q+bqs1wm1z/IRmNwdcjL+4OqOx9g0gItSyYGe7mDrFNt8JxeFs7ERU+1P7OHjQ8Idf4qHv8ObiK7Q6iesPjiAlOyEwm0Xj/Y2JftXfLYgjMjPZdB8RtRxWy+31Ha8VaIkYiiZP3yfsvHB7G0MGj9v0bSaDAXsIAAIYU4hmB/iHEQ3KNEdV1obJH/UFluo2qjaejR/ywG4kEPIKwBQAMKgQ3goPDiCm+OHM9llgghOBgjsfS2SvmdqYSCXg65Bojqr6yHcu4UIWNJ+IHwMrM/0sAnADwZgD+Ma/xvVo/s0lAYaZiJBIccXtyamMahpiOEPP6+Ml3fZStMSkMlmXxC5/6KvtyECSwLIt7Dx3PaVxRj4+P3Plncdt3P/I1rt74rM2eGujgeO3enPajyR8wE1JkiOgHAMIAPg3gBmbuJ6LlAHYwc1YddZ0iozlTUj1N4GgIhRdcM91D0UwD05IiQ0RlRFRx6ncAbwdQC2ATgI9nXvZxAE9O1hg0mpfRgRGNgslMkVkK4PFMpM4FYAMzP0dEhwA8TES3AOgE8KFJHINmBmGmUmJpz6xtkkmxHGmucCoBKq084340+cekzQSZuY2Zr8z8XMrMP87YPcx8EzOfz8w3M7NcRVsDYHwiBao2qgi0ZZrKOsYqIor1yUY8gfCAvEa1QxHRVY0rGQihTxFMcZ+UlY1VEv6cVEeHNXObvEufn255/VgwhPb9x8Q2dSq58t4BsS/TMNC5T14or5LXD/UPYVg4Fss0Uff4c2KbI3c/AIfLnkgcD4XRplioX62IzjY8+jQSfnvKQ3RoGB1PbxXb1D7ytGjv3roLYaHIOjOj+k8Pim1ObnwWlpAiw6kEwj45dUh1XrS8vpbXnzE/czU6bFmWUsZehZFMigIKRirF/i67GAIzc9uO3KKmlmXxvl/+Udy2679+yQFhPx3Pv8Sb//YWsY1KwGHvd37CJ+6812b3tXXyE//wGbFNx66DnIzZjz928Dk2g16xjSb/wUyIDp8JOjo8+0mGwoj0DWL+heeedht/SzuKF1SheMH8M9r3yxc7kRZcnaNkiw5P6tphjeYUhRXlKLywPKc2VeetHde+2EyBnCMCMEYCZKRAJRXj6k+T3+TdM0HN3IYTEXBwTNDGSIGDcrBGo9FOcA6hevQhrrXNkBSW873aPkzFus6kYk1zTLGm+FR/ElIgCQCs4DA4OmbZoGWAw1qoSCOTd07Qsix0HqkRt7Ur1tsG+ofg6+kXt6nadB2tlYuvhyPor5dl1DsVfQ01tyMqRFSZGV2Hq8U2MeH1QDoXL6FIE+lSyNu3KqLpyUgULYqIrkqqf+BoDYLCe8mWhfqH5SLndRseF+1Rjw/uejmi2q4oOxBsa0/LZo2ktApUUIQBRURTJezgbu1ERLGt49AJeVyKcxwcGoanU14jq2rTc6JOXDucisfRqyjhoOrL3dYllnDI1kZlDw97MdzenVObnup6UaQklUigp7o+p76GO7oRUqRojYe8c4IOhwOLzlktbltyvvyMqXR+JcoWyg/fVW0Wrl0Fh6BPV1hSjKqVoiYEFiv6qly2GMUVZTY7EWHRuWeLbVQqGql4Aq5Ce0JyPBCCu05Wd2nc9IJob3nqRUQUi/6HFEoxNRuegJm0q5X07DkEj1Bjg5lRc+9D4mzU09SOyJB8sTduel60BwZ9MAtKR9mICAzAo0iDUqmrVCxdJJ4XAFh83hrRrrzGqipRsXhhTm3mr14Jl5Ao7ioqwvxVy3Pqq3LpIhRXys9kVW1U9uLKclQsye1YFqxeAVeRPU/TVViIBatX5NRXxeKFKJk3cc93dXRY83LkdCwxjw+WYaJs6aLT7isZjsBVUmz7gkgEQoh5fahaO/oLyjJN9Ow+iJVvWJfzahIJc7AdVFgCx/xlo/fTVQNadZmODs9RdHRYkxWVYyhRzI6zUVguz5yK5lWgSPj2djidWP2W63LejxLLBCQJfacLsAzAeeaOVpNf5N3tsGb6sRS3l1OzcxNwCDL6BcVAKrelgZq5gXaCmqyoHpckg3KkFwD6npKfMapQrVu2UilEFcvWlBHtrE5QVs/WzG3y0glKa0ez2U+tKJiIvsa7/1z7UooORGNK6XNJ2TebnZmV0vfNjz4l2t279iNQK0f7+g/LEdXWzXKQw3uyEcPH5Eh/9055vSlbJkhwglRQBE7KTnAqzjGQJd1nhl5js/FYxkPeOUHLstC0Q/6ANCjSKjydPXC3ypL8qjbNL+2XBRQCQXQpUnRUffUcPymmLzAzGrftFdv0V8spEt6WDvQckAUcOhQpMh0KkYS+g8cxKKTCMDPq7n1IbNN53yOouvxim93X2onW53aIbRoffBIpwXEPHT6hrFdS/Sd5/+OZCTYpzktfbaNSEl91XlTn2NfTj8HG1pzatO45rBRQaNtnLyGQra/+umYE+uWEcdWxqOyB/iH0K8QgVPtv23dUKaCgEulQHkt9i1IMYlyoFhXPpJ+5KqAwHlKJBCejsQnpq3vvYe7Zd8Rmjw57ecO6d7KZSo2yW6bJ+z/2eVGSf6imgbt3H7TZLcviR9/+Yfa1dNi21f3pAa75rV1G30gkeNvXfyTuJ9VymK2kXUDBMk02u2ptds3cAFpAQTPRWIYBh2t0ckGksxtD23Zj7b98ZEL2EXN7AGaULDn9FB2j+SCc51wNctoTH6yuGjhWXz4hY9PMLnSKjGbCGesAAaDs7FUT5gABoESRXJwVy5JvhzUaBXn3TFAz/agit3HFUqekok5ztr6y7F2dEE004Q/VNbMf7QTnEErp+XE4hlCfrOybDITgPioHhlrXy2rU7YpIs5lKoW/PIXk/CmGHVFxdlNu0CDCEQEOOIhGa/CLvnOB45PXdbV0TJ68fCKJtvxy5U0nidx+pUS5udzfJxbQHFNG58NAw/IIkPQCcfFh2Nt2KaLJlmuhWCAXU/OVR0d700BNICc420NyGwd1yvZDOJ59Dwue32Xu270VQUZj7pEKSv2nbYVhC+pART6Bld7WYMK3qq/v4SaW8fq6S9N7uvomV198pv5cnFRH4CZfXVwg4aHl9HR2ecBKRqGiXJOSZmaP+gHLb7p/9VoyoDrd0iHZmZn9Pv81mGgb/5Y3/z9bGsix+/OYPce+u/bY2vS/u5P1f+p64j6du/DsOtNqjw9v+7T+45q77bfZUPM57bvu92Nfx23/B4aFhmz004OaWTU+x6R+wbXvpv38t9qXJH6Cjw5rJgAXhhe6tu7DyLdfZAie9z78ENg2c9a6bbP0c/59f44qvf87WZvBINeatXZWTvL7RuA+uC+W1yBwNgKN+OBbJyjya/EVHhzWTghSAWHXTm8TXrnz7W5T9XPXtL4j2pddcMb6BqSgoBhSrRjRzl7x7JqiZ2zgvuFa90VWY/tFoRjDpTpCInER0jIieyvy9logOEFELET1ERPqqnCIm8tFHVLGcjJmV0daw4sF8IosYg1TDGEgLOEhiqFYqJQZGgPTM1bFkfMWbNPnLVMwEvwhg5Ir6nwH4JTOfB8AH4JaJ3qFK+tvd1iXao/6AUkZd1Wa4o1t0KqlEQhm5UgoYJJNKQYSUQmFFJQlvplIIKpxNtyLdZEihOA1A2Vedovj64NEa+IV12GxZaLhfbtP4iCy7H+odQNdL+8Rtx+68VxRhbXt2O1JR+y0vM8PT2iH2FVfkKQYH3cpSBarrQmWPBUPqDABFG09nj5iBYBoGvIoMAFVfwaFhxBVfTrkeSzwcUcrbq9p4u3rFa9wyTXi7enPqK+T2IK6oVzMeJtUJEtFZAP4GwN2ZvwnAjQBO5Vf8BcAHJnKfzKxMsFUV52HOvaAPm5bcJsv+44oZT8Tjh1NYgQEA8YA8E/Ip6lV0Hzgu9+MPou7RZ8RtrS/uFnMIE6EwPMKFmAxH0PjgE2Jf1X96UHzPencdhK9JFhBoe3orIgP2xf119z0KM2HP+2t69CkE2mTBi9q/ym1CfYNKef1gnywscCp6KG7L8RoDs1otRXWNWYqZO7O60FSW/SPHazxbX7kefzZFmJzflyznZVyowsYT8YO0s7sGwA0AngKwCEDLiO2rANS+Wj9zOUVmorAsi+OB4IT0FR4Y4pe+8SNxW/3Dmzk86LbZe17aJ6a7MDPX3PsQh/sH7X098AQPVdfb7K2bX+DunfY0HMuyuOnJLZwIhV/tEDRzDExHigwRvQfAu5n5c0R0A4CvAfgEgP2cvhUGEa0C8CwzXya0vxXArQCwevXqazo75W9+jUajeTWypchM5u3w9QDeR0QdAB5E+jb4DgBVRHTq3u8sAOIDAWa+i5nXMfO6xYsXT+IwNRrNXGbSnCAzf5uZz2LmNQA+DGAbM38UwHYAH8y87OMAnpysMWhmFgHFM0FgPEIJucGWpZTx18xtpiNP8JsAvkJELQAWArhnGsYwo+BXnpHaUD0c9isCI8ysrNUbUkR6VdFRQF3fWCXGEOzpF4/FTCbR9oBcZL1v72GxZknU7UW4Xx7zYHWdaA8Pyqk7MbcHfbsPits0c5spcYLMvIOZ35P5vY2ZX8fM5zHz3zPzhH49z2QBhQZFX/2KxeiWZaFZITF+/J4HRXv79r3wKiKnJ/4qix5Ub3hCtMf8QbRtl1NUahRiDH2HT2CottFmD9Q3Y1ioMcLMOPjT/0OhUI7T09CCXoWM/Il7ZHn96vWPi3mC4d4BdD0vv/8+RbqJFlDQAgoz5mcuR4clYQPLsvjoH9eLr+89UsNRr1/cVv/4c6J9qL5ZtJuGoZTq93X2ivbBmgYeECK64Z4+3vkvX7TZvU2tfP9r3yX21fnSfq57eLPNnorFecsX/lNs07ptD6cEAYn+/Uf4ifd8TGyjyX+gBRQ0040RjaHj0c0472Mfsm2zTBMO5+SqQTMzQl29qDz7rEndj2ZmMl3RYY3mZVylJaIDBDDpDhBIL5nTDlAjoZ2gZtrJFhmeDXcqmtmNdoKzFJVzYGalgEBCEQXO5mhMRXRYEi8A1FFjAHALNYwBoG+X/JDfMgy4j1SL2+IKYYWgImoOALFOdYqOZu6Sd07Qsiw0bs+tMPZwR3fOxdebVMXXgyF0HpY/uKq+QooIJDOjV1Gvo32HHLX1d3TDXd8ibqt7WBYq6FWMNxWNofV5OdLe9NSLor1710GEB+zHY8TiaHzgCbFN4/2Pic62f+9h5XrjOkVfTRufQVJYXB/3B+Hd+Yzo8HuP1Yp99dU0IKRQy1GdS2Xx9d4BDORYfL1l96EsxdflguX1L+4W7X11zQgI67Oz7V9lD/QPKSPdykLye48oo8Mte+Tn/VNVfD3vnKDD4cB5b3qduO38t7xetC88+ywsXLsqpzbnXr9OfJZVUlmBVa+5NKe+CsvLRHsqGhNLWzIzYr6A2CY8MIywKq3j4c1yDl8qJarVDJyog5mSZ4L9R2Tn7GvrxHC9/QPia2zBwAE53cV9rAaBlg6bvePpFxH3+Gz2RDCM5s2yEw509KBzm/3DE+zuRVFiGJywz1StlCG+L8suOR8ViprHqnOpsletWIol58syXqo2a699DQqKimz2wpJinP3aK8U2F9wg6ykuu+hcVC6VV17leiyVyxZj2UXn5tRmzeuuRGFJsc1eUFyMta+/Kqe+ll54DuYtXyJuGw86OjyHMFMpUX5qKkiFI2h/eisu+If32bYd/9VduPTjH0TB/AWj7N0v7kLl2lWYd+6aUXZmRiocQWFFua0vT30zihdUoWzMB94yUohtexClb/0QqMDuWDT5jZbX1wDAtDlAACgoLxMdIACcc8M1SAz225zgqptlqX4iEh0gACy8+HzR7nAVwFFUrB2gxkbe3Q5rZh/BwweQGpafvQWq5WCKRjNRaCeYh+QqRjCRj0Sy9RVVPJg3YzEUzLdXlLMMA13rH86pL0Beb81GEnDJM2EtrDC3yUsnGFOoMavsqUQCqbhchSzXvizLUkp/Z1PwVTmPiEKS3dsuS48nAiEMKCLK/UflKGivYq0zWxYGFWtEhxVF4QeP1yEesKfisGWhcf1GsQ0AlJxjv40denEHkl57YASAsi9/ayf6hWLyHI+CikrENi2bnhftiUhUWfYg1+vCSCbF6Oh4+mJmxBTpTqo2yWhMmdaU6/7NVArJaCy3vrKIdOR8LLG48ljGQ945QcuyMNgkiyH0K1JHIh4/Qm7Z2ajaDDV3iE4tGY3B2yUvyI8IkU4ACPQNiuUrLdNEqxDpBICj98ozpOanXhBVZIx4AicUdUEan3weUaFmRLB3AK2KRexH/yQLGHTtPghvs91BemobMKTI+bNiMThL7A4q5Q9g5Qffb7Mng2G0PPa0zR7zePHMJ76Epa+xafSCEzFQUam4/9p75WMJ9A8hJjh0ABhokNNdVNdL1BdAUEgdytbG3dYlqvgYiQQ8ijo6qr7Gcywqe9QfRCDHY/G0d4sTjVQigWFFLRFVX8EBNyKK7IhxoVpUPJN+5rKAgoqkIBLAzNx36Dj72rtsdiOZZLdCKKH30HFRkj4Vj7O/SxZKiHp8orhD/9EaDvXZpfL79x3hHf/2HbGvtp/JUv0nvvofnIpEbfbIoJs3vfsfORWNjLKH+gb4xN0bxL5SvS2caDpst0djvOPrP2LTMMR2mvwAWkBBM5Np//l/Y+03vjvKxpaFE1/8Jq76zf+KbaItjYDTidK1553WPlLtNYDThYLVF5/xeDWzj3ELKBDRdUR0JxFVE5GbiLqI6Bki+jwRzZuc4WrmEswMCI8C2LJw3pc+p2zn270dVkx+LiXuJxFV3g5r5jZKJ0hEzwL4FIAtAN4JYDmASwB8F0AxgCeJSE780mgEpDXNVjQKR2mp7fmqw+VC+blrEVM8e4q1t8JRbH+OGBsalqPD8SjYJecIGoqAhWZukG0m+M/MfAszb2LmPmY2mDnMzEeZ+RfMfAMAeZGuRkT16IFZXXtVFQVTvR4AYorC7CnFh92viDQDaqECKQIMpGdwqv13PLvdZnOWlcFx0TqYwtiMWBxtD8hRYNe8+Sg9x34r3Ld1F+LCel8rFkbfQcU67CefE+2auYHSCTLzMAAQ0SVjt2VKaL78mpnETJbX71Gst217ab+Yq2Ykk2h4aqvYpu6JLaK973A13HVN4raTG58V7fWPPI1kOGKzx7x+NAuOCwBqH9okj+uxZ8SC8eGefnQ8Ix9L9Z33wllsn6U13b0+q7OX6H1xJyI99ui8t6kTg0fkxOvGvz4insuBumZljRUtrz+H5PUB1CJdHIkAlAD4DYB9r9ZuIn/yJTqsikAG+4fEbfFQmD0tHWKbQE+/aG9+bgd7W+U2rdv2iPZt3/oJe5vbbXZvexe3bd8rtjmxfqNoP/jbv4ry9jV33c87vmSXxI8Ne/nPa14n9rX701/l4WM1NnsqFOSu3/9KbPPi33+KA81t9vF+91v8wqe+KrZ5+Jq3c6jbHgW3LEuMgGtmHziT6DARlQH4GYBrAFQAWA/gZ8w8uTUSR6Cjw7OfSP8gyOFA6VhhA9OEt64Jiy63R21P/vIPuPTL/2qzR1saEa6rwZL3fdC2reb23+Pyr3zGZo/3dcOwXCg/a/koeyrgx8DBapx10xtBjrxLm9VkOFMBhRSAGNKzwGIA7VPpADX5QdnypaLd4XSKDhCA6AABINHXi6LlK8VtkgMEgOIVslTa8PPPYeX/+6B2gHOY0znzh5B2gq8F8CYAHyGiRyZ1VJo5h6UIAEl3KomBXrjmLzzt12fDDIfg3SZrE2rmBqfjBG9h5v9k5hQz9zPz+wHIT8U1MwKVI0gIBc4BIKpYnwzIaS1AeumWRKi7DynFulJPTb04Nm91HcIKWfwBoWB6ydrz0L/PXsMYAAb3H5FTZJhF6f3ySy5FpNleJ1kzd8iWJ1iZ+bWNiBaM/AFgX7g5Q7AsC10KoQCV7H2gf0gZbVK16T5+UvywxcMRDDTIax5Vwgpht0fpuAYUUTi3QsAg5gvAqyoV8Igsr9/0yFMwEvYoqJlMok0lvb5RTivp33sYfiHSngyFcfL3fxXXSNcrosPuI9UY2HvIZq+4ah3aH3xS3H/P1l0InrRH4SMt9WjbZB/zvNddh5RP/hJIhCPK86K6LlT2kNsDb1dvTm16quvF6HQqHkevUOA+W1/DHd0IK9au53osYY8Pwx3y2mVVm97aRjHSnUoklJFmVV+ezh6Es3xx50q2meCGzP9HABzO/H9kxN9ZIaJiIjpIRCeI6CQR/VfGvpaIDhBRCxE9RESFZ3gMo3A4HFhwtvy8SCWhX1JVibL58gIYVZv5q1aI8vqFJcVK6W/LkFVkEqGI6BxCA26ljL5PkdvnbmgV011SsTjqNjwhthk8WoPh6nqbvXvfUfiFhfqJYAj1Cofat/cQfEItjepf342iBVU2uxGNwb3/iCigcOJ/78Q8QZI+1N4FIyoXdBrcfwSx4/Y6G+7du2B57DNBIoKztAxGWKH8o5gJq64Llb24shxlC+1yYdnazD9rOVyF9o+Hq6gIVSuX5dRX+aIFKK6QyzjkfCwVZShftEDcpmpTtXKZfCyFhTkfS9nC+ShSHMu4UIWNz/QH6ZSa8szvBQAOALgWwMMAPpyx/x7AZ1+tr3xJkZlOUvEE9x08Lm7reHEXh3oHbHbLsthMpUR7eHBY7Kt+/UZ2n6iz2at/fTfHhr02e9zr4/1f+p6YirLv6z9kI5Gw2b21DVx922/F/e//1n/z0B3fYMscnXI08PA93HOf3Mb94hb27tklbtPkBzhTAQUi+lsAbwTAAHYx8xO5OFoiKgWwG8Bnkb6VXsbMBhFdB+AHzPyObO11ikz+wswwIlEUKIpNSRiRKIxYHMXCbIRTSQSevAfz/t+toBEz9URrDTiZQPHF9iyJpGcY/Q+ux9mf/+L4DkIz4xm3gEKm8W8BfAZADdKJ058hojtPc8dOIjoOYAjACwBaAfiZ+dQ9Rg8A8d6ViG4losNEdNjtltePamY/RJTVAUpf0q6yUtEBAoAVj8K5cOkoBwgAVFgEmPKtbeHCRcrngpr853SiwzcCeAcz38vM9wJ4d8b2qjCzycxXATgLwOsAXHS6A2Pmu5h5HTOvW7xYLhWoyQ2VsvXp3A2MZTyiA5ZhiNL/zIy4oCBtGQa8x+UgV0jxYN70DsK1wJ6TSK4CsCGn4YS6erHkPR8AK94fTX5zOk6wBcDqEX+vythOG2b2A9gO4DoAVUR0Kkn7LAByyOwMUH2ox/NhH09fubaxLEu5RtXfKb89qhQVyzSVwZROxZrqAZXis2nCo1hT3fy4vA45FYsjqIi0N6zfKCYld2zagviw3Ql6jtbAV2MP2ABAw70PinbDMwin5ASdaifY/NdHUHnVa2yzx6jXr1y7PNuusYne/0T2NZH7Hw+n4wQrANQT0Q4i2g6gDkAlEW0iImW+IBEtJqKqzO8lAN4GoB5pZ3hqvdPHAci5DuPEsixl5fp6hRMYbu/GkFD8O1ubxu37lAIKHQeP59RX85aXxOiwmUyiesPjYpvWl+QF9D0HT2BYiM6yZaH2z7KMfMNjzyAuOM7IkAdNm18Q29Sv3yiKPhz8zb1i1Lzl8WfhPWkXdkhFojj4/f9FxZqzbNsGdu5DWEgrCTS3o3+HfI5VM0G4CtLFlsa+PplE0/rHxA9c155D8CtmnKpzqbL7evqVcvHK62LnAaWAQuse+Rl53fM7RXtvbSMC/XJxqlyPxd83qEzRUbVp2X1IKaDQssueC5qtr+kQUHhLtp8s7a4AcAxANdLPEv8zYz8HwEGkZ5OPACh6tTHM1ehwKh7nwdrGnNokIlFORmM2e8zj42c/+RUxClv9l0c44vbY7PFAkDt3HRD389i7PsrumnqbvefAMbZMc5TNsize9m//wb1CX8PVdbzlw/8q7qP5rw9z19Mv2OwD+w7zsa99TWzjf/JPYkTbCHo5vGuTzZ4MhvmZ9/yTGIWOh8KcEuya2QfGEx0mImLVxhxeMxHo6HD+YhkGjGgMhZUVtm2Dew6idPlSVJxz9ih7pG8QnvvvwKqv/9g2g/Y+8gcUXf4mlF00WgHOioYRPfg8ym/4W9t+3IdPYPG6KyfgaDQzlfFGh7cT0b8T0cjngSCiQiK6kYj+gvTtrEYzbhwul+gAAWDp9a9D+drVNnvJoipwLAgrNjop3ErEEHd7kBiw6wmqAiMJ97B2gHOcbE7wnQBMAA8QUR8R1RFRO4BmAB8B8Ctm/vMUjFEzR2Fm+E/YhVAdhUUoPvt8OEvLR9lN7xCSvrCc7uIqAAQn6N13EJH2jokasmYWkk1ZOs7Mv2Xm6wGcDeAmAK9h5rOZ+dPMbK9wrRkX/MpzVBuqtBZVUXBAHTkOD8r5ltIyu5e3ReTlaSp5fSMWh0ehbK0qFZBQrKkONjTDc+j0LzPDO4iCpctQfpl9ZkcOB/w99sBA4aKFaLn9/057H5r847RE1PgVBRn/JI/njJnJ8vr1ir56Dp1QOpsWxbE0Pb1NtHtbOzBYI0fuahXF1+see0ZcI5sMR9CikGs/fs8Dor3x4c1ICPL6iWAIJx95Sh7XA3KCwODWXYjlEAU0PYNgi+CqrLRtC3X3ItRpX29dft45MEJhMUcwNOhWftloef38kdfPOyVJh8OBS972JnHbpe94i2hffM5qLBUW6mdrc9FN14upICXzKnHOtVeLbS5W9OUqKoTDJevbOgsU9sICMbcw6vHDWVhgsyeCYXQrUhGSoTCiQp5e1OMTZbFCfQNICrM3ZkbNH9fDUWDf/0B1A5ZeYStXAwDwKtKTChdUYcHVV4jbxNeffQHMZAquCrsTHNq5AwUuexpS0ZLFqLziUkQ7BAe5ZJF4jgH1daGyL1i1AisuOT+nNhfecB0KiuzqOoUlxTj/za+X+3rnDaJ95eUXoWqFLGyb67FUrViKlZfL6x5Ubc5/8+tRWFJssxcUF+OCt1ybU18rLr0A8xWiC+NBF1/XTBjMjOGaBlSdezYKys68xm/zb+/FouvWYf5rLrdtG/zzHVj6Cfta3/bbfow1X/2OLWoca61DsuEQ5v2NPZYXqK6Fd/8hrL31X854zJqZyZmuHf53IpI1gDSaERARFl9x8YQ4QABY9o63ovz8c2x2Nk0gixy+lHheuGgJChfJs4fKyy5BuNmeYK6ZG5zO7fBSAIeI6GEieidJV5hGMwlUnLtGFFdgIwVy2W+5s8IAFJcuORy47Oc/GscINfnAqzpBZv4ugPMB3APgEwCaiegnRHTuJI9NM8FI6tFAeumYiol8XJJQFHL3K1SymRlBYdkap5KI+OWIduXr3iDvnFnME7R8g4i6PeLsUTM3ON3oMAMYyPwYAOYDeJSIfj6JYxs3HkW9CpUkeCwQREQhOqBq4+nsER2EkUwq12gOCyrNAJAS1uCewt9lT/wFoJTQT8UTCPT0i9uan3xetLc8s11Ud7FME12KYEqPIgIeHhwWRR+YGXXrN4pt6v/6sFx7pLYRQ5K4AznQdUSO5vtClngsAIOD9vNitB1HvSJqHvH4lF8CqutCZY+Hwkp5e1Ubb3efmIFgGgZ8inOsusbCw14kFBkIOR9LOKKUt1e18fX0i8dimSa83fI1ruorPOxFPEtaV66czjPBLxLREQA/B7AHwOXM/Fmk6xD/3YSNZIJgZhhJOR9NNROyTEspo6ScPaUM8QPClqXMh1MpxUS9fsW4THQftOfJMTNatskCAn1HapT9Hb97gzjmpsefRUxoU7P+cbE4UzIcQY1Cqn/fr+5GyXx7dLZr+174WmTH1bVlByJC+kbnlh2I9NlTIfxtPfB2yl80rc/tQEjoiw0DVsCeJ2m6uzG4eyciQg5l1ONDRPFhN3O8xtiyshStyu0aAzPMlKIvxTVmGoY651RxLCo7W5YydSjrsShk1JTHovy8mhMre6ZaVDwigfe/AJyt2Hbxq7WfiJ+5KqAw0TRtfoGNZNJm79p5gBPhiM0eD4ZEYQFm5lQsLtol8QJm5s7te7npiedsdss0+YVPfJHj/oBtW+1d93PHs9ts9qFjNXzox3eI+xk4fpKTkah9vL1tHHzof9myRos7JI5s4SP/8wv2NraI/WnyA5ypvP50o1NkNGeKMdyPxP7NKHvPraPt7dWg4nI4l9uj0Jr84YxSZDSa6cRU1DDOFTYNOKsWgHn0LRlVLIAV1tL6cxntBDUTjuruQhJhPUVKeGjPpok+oVYwAER65cCAEU8g4fMLG5JwlpUDqdFjcFQsgOVX17CZDXdKmjNDO8EZgOphNgB4FIXcjURCGZhRydu7BWVnAPAoUlQA9ZK2gKKQeCoaw+BRe/FzAGh9TF473LN9D6ID9kCH9+BRhBRja7nvUbmv53cgNuSx2dlIph/MO8YsgyssgeWVHWqgoxsxYTmhJr/IOyc4kwUUVIvbW7ftFVVZkpEo6h6Si5wfu2+jmNvWufOA6LiYGdX3yvL6tRueQEpY3D5YXY/uPYdsdjOZxOE7/yL2deRXd4vroL0nG9ErLPq3TBP7v/dzlK1cbtvmP1EHR5E9KTrm9qDzKTndp/XhzTCkVBDDACwGuUYXACciWBGfmEPYuW0PAp1ymoYWUMgfAQV5df4sZrwCCiqyCShIZBNQUC1uX3P9OhSUltjssWEvypYuEttc9sG/Ee3kdGDeqhU2e9zrR2FFudACKKoshxlPoGDMAvdQ3wCWXnGxvS9fEKve+Fqb3Uyl0Lv3EK7+d/sa3M4tO1Ap1BEJtnbCUeCCq9guFEAELL3Z/v6ngiGUr7JXamVmpEIhlEjvmdOR1hQU2oAZVjQIZ+XCUdvi/iAsRZrIeAQUVGQTUJAYr4BCrvvPJqCQqxiDarzjFVCYSHR0WDOhMLM4Q40Ne+EqKbatKzaiMcS9fpSfZZ8J1nznx7jgK59B0aLRzinm9sB3sgErbrB/EXVs2oI173uHzW70t4H9fSi4+I32sT37RxRd+z445ssfbM3sJ1t0OO9mgprpRbX8rERRLN1VWoJyYRYMAGs/+Y8oXGhvV7J4IUoEBwgAZ7/37aLdSsTgIMXTH1chrJBXO8E5St49E9TkD+XnrVU61YjiWZ33wBG5s6BbrTzjKgCHdJrMXEU7wVmKSomamZWpKKpHH6rlXADE5XRAdkn+uEIoQRJoBQB/a4dSXMF9TI40d9wvR4e77n9EXlIVD6nTXSwLVtgeBfa3dmRN69HkB3nnBC3LQuOOfeK2hm17RbunswfuNruycLY2zTsPyNHhYAidR+QPrqqvkKL2BzOj77i90BCQLpgu4W3pwGB1vbitY5ccUWx+drvoIOL+IFqfkWX8axVrh7t2HoBXiLQb8QRq/io7ruN33ANXqV11ePDgMQwKa6dTwRB6nnhW7Cva0wePNBu0TJDqSyASEAuz9x84pqyX0rhdPpeqc+zrHcBAo6xZqGrTsuewMjqsykBQ9dVf34KAkIaUrY3KHugfUhaSV7Vp23dEGR1u3SvP3lV9DTS0KEVKxkPeOUGHw4Hzrheff+L8N9kjmgAwf9UKLBQil9nanHPd1bK8fmUFzrrSHlHN1peruEh0QslwBKSQd2/dskOcjQ3VNCAoKIwws2gHgOGmdrGv7t0H4RdUSSJDHrQLAg7MjJ3fvw0Lzltj23bsrvtRunihzR51e9Dy+DNwFo5OXWFmHP7Jr+UosGWhdNUKcUH+ive+Q0x3YYcLKLQ7WgDgRASctOddDh2tgftEndjmPCE6DqjP8bzlS7BEeF+ytVn7uiuV8vpnX2NX287W19IL1qJy6eKc2qjslcsWY+kFcjkKVZuz112hlNdf81q5hIKqryXnr0XlMvlYxoVqUfFM+tECCtNDKp7geCAkbpOEGJiZYz67EAIzs2VZbBqGvZ9Eguvve1Rs0/zIZk4Kwg7MzAdu/apo9x45zk0//rHNnqx9iY2OGrFNeOMvOXFyL1vJ0aIQ3Tv3s6e+WWyjmV0gi4DCpEWHiWgVgL8irUzNAO5i5juIaAGAhwCsAdAB4EPMrNPyZyCuokK4igrFbU6hmBIAFFfZZbSAdNRYmtU6Cwtx0T/JimznffA9yrGd8y8fEe2lq1eAA/YZr6OsElDMqgFC4SX2nLyz3iTntmnyi8m8HTYAfJWZLwFwLYDPE9ElAL4FYCsznw9ga+ZvjSYnFl17jXg7nOxpR7Fwq8RwACyJrQKYBbmymslj0pwgp+sUH838HgJQD2AlgPcDOLXm6i8APjBZY9CcOTwOB6ES75QVn9VR42zbPEerkZQiykSoXCdI7FNall9Glt4PtnaIwg6a/GJKAiNEtAbAawAcALCUmU/drwwgfbsstbmViA4T0WG3W63yIREXauJmsxvJpBiFG09flmUpZcxVbVTpLgAQGpCPPahYB2oahlJAQWUfqpWLtVumiYHjcmBgQBG1DvUNYOCwIIkPoO15eS1o00Y50p0MhtH21As2OzOj9qe/QeH8efY2A70oWGy/pIgIxoBckgBgsJAiM3ToOHyK6HCu14WZSiEVt0dHx9MXAKW8vPIai8WVatA5H4thTOyx5NgmFVcfy3iYdCdIROUAHgPwJWYe9dWdeWApTjWY+S5mXsfM6xYvPv1IkGVZ6FcsVO9TfNhDbi9Cg8M5tRloaBVnPMlIVFnnQdnXySYYghO2TBO1jz0ttql5cBMsYWbVsuUl0XEyM6ofeFLsq3HT84i47corga5etL1oF6OI+QI49idZjGHvf/8axYJz6tl3BP0HT8jj+uMGsa9DP78TrhL7ahLv8VrEBofEROpkfw84aD+XVjwG0yNHx8Es5gkOH6mG76R8zlTnUmWPeP0I9MtfaKo2Q80dosJQKh6HW6Huo+rL3zuAqKKOTq7HEvUF4O+Vv4SVx9LSITrOVCKBoVyPpW8IEUX+6rhQRUwm4gdAAYAtAL4ywtYIYHnm9+UAGl+tHx0dtmOapmg3kkkxCputjb+rl1Nxu1x+PBjiobomez+pFPs7e2x2yzT50B13ixHdvsMnuP9ItX0f/iA/8/EvsWVZtr62ffF7oux92wOPc9eTz4rH0v2bn7D73p/Z7Mnalzj8uCzJH370F2yG/Tb7idv/wD3bdottNLMLTFN0mJAu01nPzLeP2LQJwMcB/DTzvzw90WTFoVgCporaZmsjqc4AQFFFORZffL69H5cL81bb8/fI4cC6L9wi9rX8GjkXrGheBd7151+Kfb31Vz8U25SuWIqF11wJIP2ckUYeF1uwhKRcKikDFdtrGAMAFZeChBzCcz/0XhQolHc0+cNk3g5fD+CfAdxIRMczP+9G2vm9jYiaAdyc+VujEWHTtAVnlr75OrjKSuHZvXvUEjkzEoaztBRF510mdARl8XUwA077l0fZyuUorKw4k+FrZgGTNhNk5t0AVBWtb5qs/WpmH0Zm5uYSVhQMPL8dy991s80e6+nB0NPPYOEbX5HGSg70wFU5D64FS2yvp3lLQCoVGfDo2eSIcUlj0uQXebdsLt8YOws6HXuubaSgzClUkWt/u7zWmi0LYUVE210rqxF3PPEcjJg9FSbh8aLnYflpSe8DD6LyqitH2chVgIL5VbITLCxVzwQV39U9W3ciKdRd1uQXeecExyOvP9zePTXy+oq++msbxYLtlmWhQ5C3B4CWF3eL0eFg7wDcjXJdDlVf9U/IUvWJYAgtijEfv+dB0X7ywU1IBOz5e6lIFPUPbhLbNNz7gJgj0PvYUygRZPfNRALJ4WGs+NCHRtmLzz4XDgfBOd/uBNkyAcVMkJmRcvfZ7PGWOgzskhfx1z2/U7SrzrGvp39i5fV3HcyprwmX11dEbmejvP60rws+nZ98jw6bpmmLjjKn19vGQ2GxTdTrF9skIlGOBYJim3hQXgfcte+IGDkO9Q9xz8HjNnvM6+fNn/iy2NfWb/+PWJi9a+d+fuZT9vW+qUiUn3jT+zkhjPnYl77DCa/PZu/buJGHd+wQ9+99/B62TJOtMcdjePo4+uJ9ttdblsXBh2/n8OGdtm2N//MDrrv9l+J+NLMLTEd0WHP6qKK2RISicjmiWSLk4gFAoUKlGUhHeyVWKWqilC9bjHJhCVrx/Hl4z723Cy2AG38ir4Jc9abX46w3vs5mdxQW4KJPfsQWgDDCEXDKQOH8qlF2Ngz49u7FxT/7mbgfsAVYJiL1NSi/fMRxWSaYCMx2+X8jnoIrYp+9Ljx3OVBpn1Vq8ou8ux3WzFyk5GaHy4ULP/EP9hc7CBd89fM28/C27Vj41reKgQy2LIAc8Dy/2bZEjy0L4R6Pbf0wEYEKClGwzC6lVnLOBai87EqbXZNfaCeomXYkZWtXaSnK1q5GaszysPJLL8Hit70NABDr6UV88BVxTTPohaOiCv7tz6Fw0eilc2ykMHy0CTS27jAAy0jBtVjIlSQHWF7QpMkjtBPUjAuVJL+ZTCqj0KaiyHzvC3KQgZnR+ejousslK1fCUVCAcHMrmn7+KxQteeV23fQOgUrKUXrR5ShcPtqpxTq6EJPEZi0TViIBZ0WVOAatMJP/5J0TtCwL3YrF/V1Ha0V7YGBIGTlTtemprhejw4lIFINNcnRW1Ze7tROxoL34OqAWUAgLa32BdHRaldbiaZUFBPoV6smWZWFAUWS7VREd7dl/DFG3vWgRM6Np84tim/ZHN9tmfADQ/8JLiHT1in2d/M4PsOxv3jHqFtvwDiHe3Y3ytatsMz5nSRGqrrwYVmxMyouRApsGSKhJDCKlE1SdS5U9POyFt9segc7Wpre2UV47nEgoI82qvjydPcr1tjkfi8cHT2dPTm366prFSLeRTKLvpCxSoerL29WLsGfiJEjzzgk6HA7MF2rYAsCC1fLysJJ5lShViIGq2lStWCrK6xcUF6FckJEHgPmr5HGVLZwvSo8DQFRx4SbCsqNzFRYqHadPceH2H5MvtlDfIMKK+ifV6x8X7b37j2C4wV5/wtvUhu699trRlmHAc7wWccGpD+zYgwpBkr57w8NY+aG/xZKb3zrK7lq8Asm+bjhd9mePCfcwSuYXw4qNdrZWKg5W5UlmcYKq60JlLyovRakimKVqM2/5EnEZpKuwEPMU8vLKa2xBFQrL5KDZeI6lbEFVTm0qly6Cq9Au0OssKFBK5av6Kp0/D0XlpeK2caEKG8+kn3xPkZmNSGkwzOn0GWmbmUpxzGsXKWBmfumTdgEFZuY9n/oyG4nEKFuwsYlrvvE9tiyLk37fqG1Jj5sH/vhzDr/4oK2vgUcf5KHffZ+TPa2j7ImuRh68/SviuKLHd3Oi0y4goZl9IEuKTN7NBDVTg6vYXgAISKfPSNscLpcosQUAb/jNT8TI8SVf+ldbAabhHbtw4Xe+Bk6l4Htp9O118OBulF50BVzL7UWAkpEEHJUL4Fw0enZBpVWg4lLFc0xSq1Fr8gbtBDXTjkuR21h16YWwDAPmiCV1a2/9JArmVcK3cyvM6OhHAtGGWhQuWQIqtt8qpXx+OEtL4Sga89jBTMFRVCoXeVcus9PkE9oJaqadYHuXMqLc+ecHQC57Tn+ssx2l513w8t8pvxcFi5aAk3HRCa75/GfFZ7hmyA9HobxmIDrkQSKLIrImP9BOMA9R1fiIB+QItEp2HwDCQ3IUOqqIzg3XNysFGVSS/K0PPQlTWFcaamjGwLNb4RCCA4bfi8qrX1mBUlC1AEv/8RZwPApHkd0Jmu5ukNN+uZsBL5wF8sdguKYZwVZZJVyTP+SdE7QsC/Uv7ha3TZWAQvuBYzn11XW0FuFhOa2kQSEGEegfEu2RYa8yOlz9kCxgcPReWSq/71gteg8dF7cdv3+jaO/ZfxTDdfbo8HBDC1q37LDZmRltj2wWCxoNbNmOpW9/q81uGSmQ02lbNUIOJzgRFWeC5kAHyGm/vSUCHA55Fuo52YKAIq1oKgQUml7aP6ECCqprRnUsKvu4BBR2HVQKKDTvPJBTX1pAQTNuktGYaB9SFBiPBYKcUBQ/D7s9or19+14ebmy12fuP1XL/0Vqb3YgneMvf3cLJ4GihCMuy+ODHP8+WUCog3FjHffffzVbKXgA+tGWDaI+89BiHnxcEFOJRDj16uxidPv7TX3Pv1l02u2b2AS2goAGAAkUu4uKLzhPtxVlUlcsWLRDta26wFzEHgGVXXSranUWFuPwLn0JBxWihCP/Rasxfd5VYsD1SXwuYSRhBPwoWjM4xi/b1oygSQcG8qlF2hgOOeYtsfVFRCQCHGBhZfsP1KD/bvqZYk19oJ6iZdpa/+VqbrerqKzDv8kvE18faW8ARHwy/z+YE4wNDiHd1oODyq0bZrWgIzsoqeQAk1CoBsOT1srqOJr/Iu2eCmvyAiOAoLECsfxBm4pWlY2yZsEwTrnlVKFxmL/aU9AWQ8tqDORyLykvjAICcYmH2SGfPqBommvxEO8EZDk/gAn5VX5JC9Zm0UQklGHE5amzEE4gIy/PYNNHwy9/DWTQ6YXrlxz4NchXAWWrXWiyoLB+VOnMK4mRaXXrsPhJxwOEAp+xjGz50DIHGVnHMmvwh75zgbJTX7zxSI0aHLdNE2zZZ3r3vuCx6kAhHEBPk7QGgfcc+0a7ah5FIoPdojdxm6x7R3rn7IMID9uLnqVgcJx95Smyz5/u3oUh4/tjyh7/aZmLkcKJw8VLbawGA2QKbJhwF9jWqxCYQs6cIJZqPAUTiTDB48iQCh+T3pj5HSXotr6/l9XV0eBxEPD4earAXH2dmDvYPiXbLssRIp2kY3CMUP2dmbt95QLRHfQEO9A6I21q37RHtx+97TIxCd+05xO079tnsfQeO8WPv/bjY14FPf4VDrR02u2UY3PuH2zgV8I22J+Lc/qOvcdLnHW23TI489yeOvvSwra/wjsc4vPkuNiP20gMtt/8v9973R3FsmtkFdHR4dlK6oAqlCrWOCoXyhrj8C4DD6cTKqy8Xt615k132HgBKqipRolDXOeetbxDtV/7T34r2VW9YJ9qXvfZKvGfDnTZ7sKEFJcuXovycs23bjFAAKChEtL4ala9/88t2KiwCVS6xzQSJHHCtuhCW1z57oMJiMABHqb30QFGpEy6HfiaY72gnqJlWiAiFQh2VigvOwcXf/HexjeH3wvR7YSXst1fx/j6YiTicZaP7pIIisGkXgi1YdQESHllibP6114mrTDT5hT7DmhkJORxKYQXD7wUbKRQssj8bTHq8iHXbl7pRQTFg2suaOioXgBNxsKAW4ywsgKNYXbhKkx9MmhMkoj8R0RAR1Y6wLSCiF4ioOfP//Mnav2Z2kVQIFQQ77A7NDPgAhwOlF15m2+YoKU1HfMdSWAQIJQEc5fPSkWFhlgjLAMXl9daa/GEyZ4J/BvDOMbZvAdjKzOcD2Jr5e8JRyZh7Bal2AIgFQ8qIqqqNt7tPTB8xkkkEBuQ1mqq+ktGYaAeAlCqtRCFSYJkmTEX9j7jK0SjWlAJqGX9vm7ymNh4MKQUZ+hUCCi1PbkGg3e7sAm2d6NlqXwdedvnVcBSXiKtJKq98DUpW258jUlEJ2LK/L+RwgtgEJ+znwEomYPnsdUkAYFARHVWd43g4opS3V7Xx9Q6IqUiWaSojvaq+wh6f8jpTXuMKeyISVcrbq9r4+wbFbArLNJWRXlVfEa8fCUUJifEwaU6QmXcCGJv38X4Af8n8/hcAH5iE/SpPdiIi282UASNpv1XK1iYVi4tO0DItGAk5T07VV9jtFS92I5FQSjlFfbLTtkxTWWPE3yV/OfQerhbtiXAEvg778zLTMNCiSEM6fNcGlC2xlxfoPXQCQ0Itib59h3Ho53di0eUX2bbV/v6vIJfd0bkq59t1ATOwYYrSWygoAYQ8QQAgWLACguhE1A/E7e+/ZZqoefBJsS/VObYMQ6wXkq1NKha3lQ49tf+UkG6SrS8zmYKZkq9x1edFZbcMA+Z4Pi/CsbBlKfNHVX0ZiaTyWMaFKmw8ET8A1gCoHfG3f8TvNPJvoe2tAA4DOLx69erJiZtrJhwpPeeU3TJNm33gSDW3bn5BbLPnm//N/pZ2mz3R38N9f/4/cV/td95hS5FhZrbiEQ488htxP6EHf8rJXnsqUmTP0xzb+6S9L0Uakmbmgpkor58ZmHI5BDPfxczrmHnd4sVyOohm5qFK0SEisWD60qsvxznvuVls84af/gfmnbvGZk8M9MII+GGG7bNhK2XAEGYwTE6EuuXbLiqtABzC7NEygEJ7YISIlMepmX1MtRMcJKLlAJD5X/0wSqMBkBj2wBpz65Uc6EGipxuGz77KJuXxINbRYe+ntwvxAfk5GhWWgMMBm52NFOCyrz7R5BdT7QQ3Afh45vePA5AfrGjmHKpi7h333g9H4Wjhg0R/DwqXLkeBsHzOsghJr/2hPRUWwVWqKNPocMAK2h0qzBSQionPfiWbZnYymSkyDwDYB+BCIuoholsA/BTA24ioGcDNmb814yBbRFn1AVWJHvArz2FtqB5Aq6LTnCUwpQrYBNu70L/nkN1e1wDvPru9aMkKOIqK4SyxO7WCefPhFJydo6AQrnL7qhAAoFQCZkBQnjFT6fXGhj2gMahYU62ZfUxmdPgjzLycmQuY+SxmvoeZPcx8EzOfz8w3M7Pw9XtmjEdAwd3WNTUCCs/tEO0qeX0A8CvSB04+9oz8+u4+UV7fMk00KeTSG57dLtoDvQPo3CcfS83DshjCiQ1PjJK+OkUyGkPdE1ts9oQ/gBc+9gUsFpb0hRqasezdb7PZF/7NB6F6nOwoKkTlZXYBVyuZAIRIMwBwKgGW8gFNIz0THKMwYySSOLleLi+gul683X0TK6CgkqRXXGMTLqBQ05BTGy2goAUUJhxve3fObbJFbiVSiQSn4nKRdZW8vr+7T7QnozFxW929D/Kur/xAbHP8y9/mVChss5vJJPf87jaxTcMP/stWsJ2ZOdraxN0//7rYJvLErzl2/CWbPXboOY7v32R7fwJdvbzlc98W+9LMTKAFFPKP+Wtyl33PFrmVcBWqgwIqef15Zy0X7QUlxeK28z/0Plz4UbvoghEOg4jgEtYVGwEfXPPkxUZmMoFEV4dNU9CMx8WKdgAAhxPFV77ZZqbCYlDStL0/latW4O13/kTuSzPr0GuHNdOKq7RELKmZ8gex6sN/BzMehxEenbBs+L1wlMnqNhwNI9bZYe/POyhqBnIiohwbJ2NAgZyUrckftBPUzEhKzlqB+a+9Gt33PwjHmBmpEfAh0CQv2ytetggpj/2ZqMPlQum559rslrcfIMXHIBYGlcjBFE3+oJ2gZtphyxKFEhLuYfQ/udnmBMnhxNDOQ2JEu+Kii1F2kb1Ak6uAULRAuIVPRkGFRfLAklGgWHaCcZ89r1AzO8k7J2hZljJy1qiQl/d29WJYWMCfrU3L7kNydDgYQtfRWqEF0Lhdlmrvra5HRFiQzsxo3iFHznyd8uJyI5FAPCgrnwQUETVp30B67XCXYl3xcEuHch/SulbLsuBRFDLf/4NfwCGs9w03t2D+a+1irMXnXozEsBcsrF81QkG4yoTUmaoFcopMYUl6AaeEaYjPS5kZJx+QU1xV14u/bxCDTW05tWnbd0QZHW4/cEzuS3GNDTS0ICjUccm2f5U9MDCEgYaWnNq0HzimjA6rsilUfQ00tioLyY8LVcRkJv3kGh1OCdHBbHbTMNgUinyPpy9mZiNpL/6drY3q9dnaeDt75NfH48o2AYUkfyISzXn/Ktn94IC8j2QszuFh+5retqe38qM3/r3Ypvm2OzjS3mEfUyTKJ772PbFN5x0/5fAxu/R/rOkEh3ZstNkN3wBHn79X7Cv2wr1s+O3H4+/o5o0f+ozYRnmNmSYbqVRObbJdY7m2MVIpNoW12+PpayYfiwrMteiwKqqpsjsESabx9gUATuFBf7Y2qtdnazN/tb3cJAC4ihS3dgAqFZL8hQrx0mz7r1whFzuqWCrvo6C4CAXF9rEtvPQCvPX/fmyzMzOinV0oXWOXxYq0tqNcWFMMABz1wxKWwFmJGEh4n8lVoL4dBoEK7YGRylUr8PZf/0hsobzGHA5AWDudrU22ayzXNk5JWWecfc3kYxkPeekENbOHyrPlVJ9QfQMqLrHLawFArLcP5eedI24jZtGpcSIOkqrQFZbAUaT4EigqEYMm5HCgbOkiuY1m1qGdoGZGUnLWWVj5Dx8Uty17581gZsQH3SgeM/NkIhTMt+sZmvEYCkqEmYVLPXMGkVKDUJM/5F1gRDM74TGR3oLKChTOrxKjxgAw8NxWJIbsD/qLli4Dol5bf1Y8BjLsD+bJ4VCumyZygJPqNdKa/CAvnWA8LCfAquxmKqVU/c21L84iIKBqk4zFleIGCUUbKTINpFWfVR/QXPvKdiwqYQVmVpYECCrWrg4eq0VIKIkQ6R1A1zNb7ftOJNHw41+iZKV9BUrp2nOAoNsmeuAsLgUi6QJNY7EEXUI20++jOSCvKe/Ze1i0K68xwxAjvdnaqOwAlPLyqjapeFxddiHXz4thIBWXV9+o2mSTw8/5WBIJ5bGMh7xzgpZloU+xuLvnRL1oDw4OK0Puqjb9J5tE55EIR+BWpIKo+vJ09CAWsKe1MDP6quVjCQoiCUA6506VCtOr2H9ocFi0J6MxDDXJTsAryO4DQMeeQ3AIQgWB7j60Cekb4YEh7Pzuz1GxaoVtW82d98JRZL+FTfn8KF6xDAXzq2zbOBkDhGBG2VXXgmNhsCCXb3jt554TUVipFEy//X1OxeKo3fCEzQ6oz3F42At/j3xeVG0GG1qVKTKDjXK6jaovb3e/ssaJqo3KHvH64cv1WBrbRMeZSiQw0NCaU1++7n5lWtd4oNkwrV+3bh0fPix/82pmD5ZlpSOLIwj1DWDwSA3Oe69dLWbbJ7+Ma77zBcw7b+0oe8+jm+AsLcFyQWHGv/k+lK5egcIrbxq970gAyUNPo+gtHxmV+2eZBoJ//iEqPvINOEcUYLe8fYgfeRFwFqL0xg+P6isZjqB3/1GsvflNp3/wmmmFiI4wsz3pFHk4E9TMXMY6QACoWLFMdIAA8Ibb/tPmAAFgePd+LHnrG5ESlKWtVCJdY3gMpxyfLfmZCM4Fi0EFo2OEbKTgqFgAR6V9lUlheZl2gHmEdoKaGUvxAlkp5orbfghHcTEGnnjcto0ckNf7KpRyHA4nyFUIGvOs0LnkbFCpLNKgyS+0E9RMO5ZpKoM5kR57wMThcsGzfRuSbkEowQFQsV1+iy1LvTyOHLBiUmnTmf+oSHPmaCeomTJU9WW7t+4GBCcYHRhC3wuyGnaothqLbrJXqSMikCR6kErKFeUAwOGE2SsUU+eX/7GhiqhrZh955wRno7x+55EaUV6fmVGvkMT3CSklQDrdwNPWJW4bUkQUm7btFVMOzFQKvi5ZqCGkWIzv6+oVZ3UxXwB1G+WSACfv2SDauzc/j7hHeO5nGEh5vZh3jf05N6lmgkZaVFXSFAQzrEH7e8bMsBRpLc05XmNaXn/myuvr6HAewsxK9ZPpqpfrbe0AACwYs+Y37vNj/ZU345MdB211iWt+8TtUrFmFNX/3nlF23/59SLrdWPre99n2k6rZAdelb7b1ZfY1I7RzMyre9Qk4540OdkSeugvOokIUv+0To+17NiPlHkDVBz5tP562Liw4Z3WWI9bMJHR0eI6Rq4z+VLDg3DU2BwgARVXz8L5NfxYDFwVlZVh87TU2u2fHNiy84a3yjpjFIu9WMoGUZ1h+9udwilHg1PAQzJCsG6gdYP6g1w5rphUiwuKrLhO3XfSZj9lsqWAA5HDAVVEBtiyFw0vahFg5FYezYh5ci+1J2bBMwGVXmLFSKTjK5p3mkWhmK3omqJlVxDo7seht70C0pw/hltHPcdkyYaYMuHfstrUjZjhKK0BO4XvfMuBce5XN7CwoQMFZ503U0DUzFO0ENVOG6vlzKhJFMmi/TU1FojYBhcrLr8C811yN9j/cC2fJmKToeATDx5qQ9AiF1OMRwKnQbbRMOCrsyjNWKqXUGmTFWm/N7CPvnKBlWeipltccdh8/KdqDg27l2mFVm96aBnntcCSqjDR3HZNl991tXYgJkvjMjJ4TdWKbnuOyPeT2IKCIAobdducAQCnHn16jKq/rVC56V/QFAN17j4j2E7/7KxyFdgc1fPwkhk/Y338jEsXAcy+ioLJilJ3jEfhqWuAssesDWmE/nMJzv3RDS9QaRCoJKpa1BttetM82AfX1Evb44Ovpz6lNX12zKFSRSiTQXy/L26uuMW9Xr3LtsGr/KnvE64dXkTWgatNf3yJGuo1kUhk1V/Xl7e5DeALXDk+LEySidxJRIxG1ENG3JrJvh8OhVFCep1BDLq4oR3GlXFBH1aZi6SJRkdpVVIjyhfJKh6qVy0R72YJ5KBw7q0H6edm85UsU45LtRWWlKKqwp4gAQNxvV0sBACMhK+i4CgtQrqgvrCIVi4szPk9LB3zt9jSUqNuDuvsehUtQnR46Ug1PrT1/z4zFsPTtN8FVMfqccTwMBmHpu+zL8KxIGI4q+2wvG5aRgkMQY0hFY6h//Dmxjep6KSorQcm8CnGb8hpbvECsveIsKEDFYvm8qK6xkqpKFJbJDl21f5W9sLQYJVXyappsxyIpRTsLClC5RD4vqr5K5lWgSHEs42HKAyNE5ARwJ4C3AegBcIiINjGzPLUZB5VLZNVflT2bvHyufTldLpQtqMqpTWmV+uF7haKNyp7tWBadb1+HCwDli+WL0OF0okzh0IuEYkaAWl5/4XlrsPC8NeK2q79wi2iff9F5sISZkGfvQSy6/vW2oAjHI2CL4RRKDFiJGFzCTJCZ4Zwvf9ic85eABNVpV0kx3v2bH4ptVOe4oLhYWcJY1aZCdV4cDuWXk6qvkkrZAWdrk/1Y5INRtVGNl4gm9FjGw3REh18HoIWZ2wCAiB4E8H4AE+YENbOH0sULcdktHxG3rX77W0S779BRnP/lz9nslmGopfIrFsNZaXfoHA/DUSk7GyquEGeCRASa4DoXmuljOs7kSgAjn3b3AHj92BcR0a0AbgWA1at1TpbmFS7+/jfF20TXOVfh/K/L14qVSIjPBDkSACnSYNhIisWZNPnFjA2MMPNdzLyOmdctXizfYmnmJpIDPGUvXio/K7ViEVCx/Raeo0GlWgynknLARJNXTIcT7AWwasTfZ2VsGs2kIi4ljATVM0HtBOcE0+EEDwE4n4jWElEhgA8D2DQN49DMIZT5ftEAHKWyEyy59BqQsJJEk19M+TNBZjaI6N8AbAHgBPAnZpYTgjSaCWLeuz4q2jkRS9cXFii58MrJHJJmhjAtIS5mfgaArKuk0UwhBRdfO63CEprpZ8YGRjSaqcAxX04u1swdZoWeIBG5Ach1LF9hEQC5duTUoscxGj2OmTUGYG6O42xmFtNMZoUTPB2I6LBKNFGPQ49jJoxjJoxBj8OOvh3WaDRzGu0ENRrNnCafnOBd0z2ADHoco9HjeIWZMAZAj2MUefNMUKPRaMZDPs0ENRqNJme0E9RoNHOavHCCk6lU/Sr7/RMRDRFR7QjbAiJ6gYiaM//LqqQTN4ZVRLSdiOqI6CQRfXGaxlFMRAeJ6ERmHP+Vsa8logOZc/NQZr34pENETiI6RkRPTdc4iKiDiGqI6DgRHc7YpvS8ZPZZRUSPElEDEdUT0XXTcH1cmHkfTv0EiehL0/F+jGXWO8ERStXvAnAJgI8Q0SVTtPs/A3jnGNu3AGxl5vMBbM38PZkYAL7KzJcAuBbA5zPHP9XjSAC4kZmvBHAVgHcS0bUAfgbgl8x8HgAfAFlGeuL5IoCRxWamaxxvZearRuTDTfV5AYA7ADzHzBcBuBLp92VKx8HMjZn34SoA1wCIAnh8qsehGtys/gFwHYAtI/7+NoBvT+H+1wCoHfF3I4Dlmd+XA2ic4vfjSaRLF0zbOACUAjiKtFjuMACXdK4mcf9nIf2BuhHAUwBomsbRAWDRGNuUnhcA8wC0IxMEna5xjNn32wHsme5xnPqZ9TNByErVK6dpLACwlJlPlRUbACAXsJgEiGgNgNcAODAd48jcgh4HMATgBQCtAPzMbGReMlXn5lcAvgHgVF3MhdM0DgbwPBEdySilA1N/XtYCcAO4N/N44G4iKpuGcYzkwwAeyPw+neMAkAe3wzMZTn+9TUkOEhGVA3gMwJeYeVRZuakaBzObnL7dOQvpWjIXTfY+x0JE7wEwxMxyfc+p5Y3MfDXSj2o+T0RvHrlxis6LC8DVAH7HzK8BEMGYW84pvk4LAbwPwCNjt03lOEaSD05wpilVDxLRcgDI/C8XNJ5AiKgAaQe4npk3Ttc4TsHMfgDbkb7trCKiU5JtU3FurgfwPiLqAPAg0rfEd0zDOMDMvZn/h5B+/vU6TP156QHQw8wHMn8/irRTnK7r410AjjLzqeLY03adniIfnOBMU6reBODjmd8/jvQzukmD0mJ49wCoZ+bbp3Eci4moKvN7CdLPJeuRdoYfnKpxMPO3mfksZl6D9LWwjZk/OtXjIKIyIqo49TvSz8FqMcXnhZkHAHQT0YUZ001IV3ac0nGM4CN45VYY0ziOV5jqh5CT9KD13QCakH4G9R9TuN8HAPQDSCH9jXsL0s+ftgJoBvAigAWTPIY3In0LUQ3geObn3dMwjisAHMuMoxbAf2bs5wA4CKAF6Vugoik8PzcAeGo6xpHZ34nMz8lT1+VUn5fMPq8CcDhzbp4AMH+axlEGwANg3gjblI9j7I9eNqfRaOY0+XA7rNFoNONGO0GNRjOn0U5Qo9HMabQT1Gg0cxrtBDUazZxGO0HNrIKISojopYxwRrbXPXMqb1GxfTERPTfhA9TMOrQT1Mw2PglgIzOb2V7EzO/m9MoV1XY3gH4iun6Cx6eZZWgnqJkRENFriag6o0tYltEkvEx46UeRWVVARDcQ0U4ieprSepK/JyJHZlsHES16lX6fyPSnmcPoZGnNjIGI/htAMYASpNe7/s+Y7YUAuph5WebvGwA8h7SOZGfm9z8w86OZtcPrmHlY1S8RrURaZ+/yKTg8zQzF9eov0WimjB8ivRY8DuALwvZFAPxjbAeZuQ0AiOgBpJcRPnqa/Q4BWHHGo9bMavTtsGYmsRBAOYAKpGduY4kJ9rG3MtKtjarf4kyfmjmMdoKamcQfAHwPwHqk5fBHwcw+AE4iGunIXpdREHIA+AcAu3Po9wKkxR40cxjtBDUzAiL6GIAUM28A8FMAryWiG4WXPo/0Le8pDgH4P6Rlu9qR1u073X7fCuDpCT0QzaxDB0Y0swoiuhrAl5n5nzOBka8x83vG2ddOAO/PzDA1cxQ9E9TMKpj5KIDtr5Ys/WoQ0WIAt2sHqNEzQY1GM6fRM0GNRjOn0U5Qo9HMabQT1Gg0cxrtBDUazZxGO0GNRjOn+f8DeYdYIKXjtAAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"d.isel(t=0).w.T.plot.pcolormesh()"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 313
},
"id": "T9VcQT_r2NPV",
"outputId": "4569037c-6d22-426d-db0c-0e90d0fc0781"
},
"execution_count": 54,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x7f1c7f001970>"
]
},
"metadata": {},
"execution_count": 54
},
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
],
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqlUlEQVR4nO3deZhlVXnv8e/vVFUP0C1Td9qWbiZBiBgFaXBWBDTojeM1BM1D0BDbMRevJkY0DlET1IjE5yaS4JAQRRFwgPAYEBHHKMokU0MQAwI2tM08dnXVee8fe5duirNW7ao6Y53f53n20+fsfdbe65yqPm/ttd/9LkUEZmY2vBq97oCZmfWWA4GZ2ZBzIDAzG3IOBGZmQ86BwMxsyDkQmJkNOQcCM7Mh50BgfU/SbpJC0v2S1lfWv0nS7eX6nXrcx8PKfjQlHdbLvpjNlgOBZUm6sVNfbJJeI+kmSQ9I+rqkHWdosn1EnFy2HQM+AbwwIpZFxB2d6GNdEfGtiFgG/LKX/TCbCwcC6wlJ+wL/AhwFrAIeBD41i12sApYAVyf2PzrfPpoNCwcCS5L0eWAX4D/KYY93tnH3fwz8R0R8LyLuB94LvFLS8hr9egJwXfn0bknfLteHpLdIuh64vlz3SUk3S7pX0iWSnlPZzwcknSHpC5Luk3SlpCdIOk7SprLdCyuv307SZyVtlHSrpA9LGmnfR2LWGw4ElhQRR1EMdbykHH752PTXSNpF0t2Z5TWJ3e8L/KxyrBuAceAJNfr132V7KIaLDqlsfjnwNOCJ5fOfAvsBOwJfBM6QtKTy+pcAnwd2AC4DzqP4f7Ez8EGKs5Yp/wZMAHsC+wMvBP5spv6a9TsHApuXiPhlRGyfWb6YaLoMuGfaunuAGc8IZnB8RNwZEQ+V/ftCRNwRERMRcQKwGNi78vrvR8R5ETEBnAGsBD4SEVuB04DdJG0vaRXwYuBtEfFARGwCTgSOnGd/zXrO46jWK/cDj5m27jHAffPc783VJ5L+AjgGeBwQ5TFWVF5ye+XxQ8DmiJisPIciaD0OGAM2Spp6fWP68cwGkQOBzSRbp1zSLsA1mZe8ISJObbH+auAplf3sQfHX+n/PpZMVv+lveT3gncChwNUR0ZR0F6BU44ybgS3AivLswWzBcCCwmdwO7JHaGBG/pPiLebZOBX5UfllfSjEe/9WImO8ZQdVyijH9XwOjkt7Fo89CaomIjZK+CZwg6b0UZzS7A2si4rvt6rBZL/gagc3keOCvywu/f9GunUbE1cAbKQLCJoov7Te3a/+l84BzKc4ybgIeZn5DOX8CLKI4A7oLOBNYPc8+mvWcPEOZ9TtJu1Kkiz4M/GVEfLrHXXoUSYcCX6EY3npxRFzY4y6Z1eZAYGY25Do6NFSm3Z0p6VpJGyQ9Q9KOks6XdH357w6d7IOZmeV1+hrBJ4FzI2IfigyRDcC7gAsiYi/ggvK5mZn1SMeGhiRtB1wO7BGVg0i6Dji4zMJYDXwnIvZO7AaAFStWxK677NKRfj5a6vNIZxxGZptS+4vmLPpU3WHr2B3KZERmfsbJ/pVbZyv/WbR+z7k2ZN6XEu8ru7/c+53Dsaz/XHrZZZsjYuVc26/V0ng48bs63WbGz4uIw+d6rH7RyfTR3SnS9v5V0lOAS4BjgVURsbF8zW0UxcMepSw3vB5g7dq1/PCHP+xgVyvHndjScn2MLk62aWa+eBqJL57G+APJNpH4sgdgZKx1HxLrASab6S+xseZ4+lipoNNI/9rkPouRiYdbrt86kv5sRxqZz3Zya8v1k430Z9HMfKHnjuVAMDiWbrPNTfNpv4UmR6heMtin4qYVM7+q/3VyaGgUeCpwUkTsDzzAtGGg8kyh5f+wiDg5ItZFxLqVKxbEZ21mA0DAiFRrWSg6GQhuAW6JiIvK52dSBIbbyyEhyn83dbAPZmazNqJ6y0LRsUAQEbcBN0uaGv8/lOJGnLOBo8t1RwNndaoPZmazNYxnBJ0uMfHnwKmSFgG/AF5HEXxOl3QMxd2eR3S4D2ZmtUmwKHO9aCHqaCCIiMuBdS02HdrJ45qZzVVxRtDrXnSXi85NMzm6pOX6kWbrDBWA3BRVqWye5tjSZBtlUktTGUXZzKBMscxcBtBEYuQw90uTygwCaI4lPtvM/v7v0n2S2z7+wIaW63d8Zrpk0Z0//Kfktlxi0JB9Lwy5hTXsU4cDgZlZhRi+apwOBGZm0/iMwMxsiGmBpYbW4UBgZlYhnDVkZjbUpu4jGCYOBGZm03hoaMiNTCaKzs0hzRJgNJGTmGvTaKQTKhsx2fo4ueKjyvQ9kzKZKsLW2HJ/ulFGKj1zMpO3+Q8Ptk4Rzbn7vz6V3LY1k2YbmcqkQzZSMNSKawTD9QN3IDAzm8ZnBGZmQ6yB2naxuKy19uXKqj2A9wHbA6+nKNUP8O6I+EZbDjoHDgRmZtO064wgIq4D9gOQNALcCnyNou7aiRHx8fYcaX4cCMzMKjp4jeBQ4IaIuEl9dg1i2O6kNjPLmio6V3M+ghWSLq4s6zO7PhL4UuX5WyVdIelzknbo4FuakQOBmdk0s5iPYPPUTIrlcnKr/ZWl+F8KnFGuOgl4PMWw0UbghM6/qzQPDXVYalL5XATOzavbSM1nnKlYmpp7GPJpp42H7225vrlo2/ShMv1Izd+cmyY8Nwdy6nMai3Sl2K2ZX/lMZqnLjw6RDpWhfhFwaUTcDjD1L4CkTwPntP2Is+BAYGZWIcFYo+2DJa+mMiwkaXVEbCyfvgK4qt0HnA0HAjOzRxBq4ymBpG2BFwBvqKz+mKT9gABunLat6xwIzMyqBI02BoKIeADYadq6o9p2gDZwIDAzqxCgkeHKo3EgMDOrEm0dGhoEDgRmZlVSW4eGBoEDwTRbG4tars+lWY5Njie3hVrvL59mmRZqXZlUmVZqpievz5lY/JjW+8tVOs1MRT+ZyM8czSSQNsYfSm5LpbE+1Ez3IZeau2Q0+8lnttlCIsHIWPp3aCFyIDAzm8ZDQ2Zmw0zyxWIzs2Em2ps+OggcCMzMqgQasinpHAjMzKokRhb5YnHbSLoRuA+YBCYiYp2kHSlm7NmN4tbqIyLirk72YzZS2UG5zJsYaZ0ZBJmic5OZTJ5MkThi9mOXzZGx5LbGROs5miFdJO7hzETHo5m/pMZi9tlLubmiR6/5dsv143senGzzk1+l51s+bLfWWVI2XDSE9xF044rI8yNiv4hYVz5/F3BBROwFXFA+NzPrG42RRq1loejFO3kZcEr5+BTg5T3og5lZayqKztVZFopOB4IAvinpksrMPasq5VdvA1a1aihp/dSsP7/evLnD3TQzKwhoNFRrWSg6fbH42RFxq6TfAc6XdG11Y0SEpJYDzuVMPycDHPDUp/q2TjPrDrnoXFtFxK3lv5skfQ04CLh9alIGSauBTZ3sg5nZrEiMLBquQNCxdytpW0nLpx4DL6SYheds4OjyZUcDZ3WqD2Zms6XyjKDOslB08oxgFfA1FemTo8AXI+JcST8FTpd0DHATcEQH+zBrqTlyG5k0RmWKzpFol0s5zUmlsebSLMcn0yNri0YXJ7c1tj7ccv3SRqawnNKpqs1G622NyfQcww830v1jn0Narh7NTD588K7pFNFUqi+AMsXqbOHxncVtEhG/AJ7SYv0dwKGdOq6Z2bz4zmIzs+EmtKDuEajDgcDMrGoI7yx2IDAzq5JojA3XV+NwvVszsxlItHVoaBBqrg3XQJiZ2YzUifTRvq655jOCaVIVN3NyqZvN1BzDmSFINSdn3YdcddTFmf6RSYtsji1pub4x/mCyTa5477ha90OJ9QBjuZTOxKYHM+mjuZ9uak5lyM9ZbQtMd+4sfhlwcPn4FOA7wF91+qApPiMwM3sEoUaj1gKsmKqJVi7rW+xwzjXXusVnBGZmFZJoLErfGDnN5spwT8qca651iwOBmVmVoNFo32DJINRc89CQmdk07bpYPCg113xGYGZWJbXzYvFA1FxzIDAzqxBMXQiet0GpueZAUFO2MmVmhK0RqVTQdBtlJpSPsaUt12/NXGpalEktJZrpfiTSWCfHtkm2mcykoy5KTF6fS79tXPz15LbmA/e2XH/eqhcn29w3nk7NPerJucQNVx8dGu09IxgIDgRmZlWCkUXD9dU4XO/WzGwGkto2NDQoHAjMzKbx0JCZ2TDzNQIzM/PQkJnZEJNEYyRXPnHhcSCoaS4VQYt2iZTJzOT1qRRRSKeWjo2ka6NszVX3JJNaSuv3nE6Jzd+qnupHbl74sXUvS25r0jql91Xj9yfbbBnbNrltfDKdSrt4tPU786T2C5Cg4awhM7NhNjhZQ+UENzNpRsTduRc4EJiZVag78xG0y6/KJTdjxgiwS24nDgRmZlWDlTW0ISL2z71A0mUz7cSBwMxsmkEZGgKe0Y7XDMy7NTPrCgmNLqq19FpEPAwg6fPTt02tm3pNjs8IOiyZHZQp9paVm+w4oZFrM5nJhkpkIk3kiuxljjWayDZKzesMkPuUmomMnQdH0kXxckmBow1PTGxQzkzT607M1r7VJ5JGgAPqNnYgMDOrEmhA7iOQdBzwbmCppKlyvALGgZPr7qfjYU/SiKTLJJ1TPt9d0kWSfi7py5J6f35lZvYbgsZIvaXHIuL4iFgO/H1EPKZclkfEThFx3NTrJO2b2U1XrhEcC2yoPP8ocGJE7AncBRzThT6YmdUjBiYQTKl+6Sc86hpCVUcDgaQ1wP8CPlM+F3AIcGb5klOAl3eyD2Zms6HyhrI6ywDJXgDr9DWCfwDeCSwvn+8E3B3xm6mqbgF2btVQ0npgPcDatWs720szsykS9EFGUJtla6F0LKRJ+gNgU0RcMpf2EXFyRKyLiHUrV6xoc+/MzNIW4BlBVifPCJ4FvFTSi4ElwGOATwLbSxotzwrWALd2sA9doUwqaChRrGxya3qHmbHHyZHFLdenUikhH+1jtPX+csa2PpTeX6Zg3pZo3ZOxzEnr1mb6fS1OvLHJTArrwxPp/W07tnD+Y9s8SH01/t8m47mNHfvNj4jjImJNROwGHAl8OyL+GLgQeFX5sqOBszrVBzOz2Wtf1pCktZIulHSNpKslHVuu/4CkWyVdXi4vnnevpR0kHSTpuVPL1LaIeHqubS/uI/gr4DRJHwYuAz7bgz6YmbXW3vsIJoB3RMSlkpYDl0g6v9x2YkR8vB0HkfRnFBmaa4DLgacDP6JIzplRVwJBRHwH+E75+BfAQd04rpnZ7LXvzuKI2AhsLB/fJ2kDiQSZeToWOBD4cUQ8X9I+wN/VbexBUTOzqtnVGloh6eLKsj69W+0G7A9cVK56q6QrJH1O0g7z7PXDlbpDiyPiWmDvuo1dYsLMbLr6ZwSbI2LdTC+StAz4CvC2iLhX0knAhyjSOj8EnAD86Rx7C3CLpO2BrwPnS7oLuKluYwcCM7MqCbUxa0jSGEUQODUivgoQEbdXtn8aOGc+x4iIV5QPPyDpQmA74D/rtncgqCuRBgr5OzXScxan5xjOSc0XnPvFbfu0upnPQpPpLLXFyUqs6QqoizIX7VLzNy+aQ0oswGTmgxpJ3JjpeqULUfvSR8tqCp+lmEDmE5X1q8vrBwCvAK6a53E+HxFHAUTEd6fWAUfVae9AYGZWJdpZhvpZFF/GV0q6vFz3buDVkvaj+DvyRuAN8zyOy1CbmbWLJDTWnhITEfEDWp84fqMd+29XGWoHAjOzRxicO4sj4njgeEnH16hAmuRAYGY2zaDUEZL01PLhGZXHvxERl9bZjwOBmVnVYNUaOiGzLeinO4vNzAZKJjOun0TE89uxHweCuuY62XxCrmLpnCQqe0I+xbGZ2dpIJcZm+j5yz8bktuaylS3X/zq2TbbZdix9rGXX/6jl+tj7Ock23//l/cltz1yzPLltNFUitd2pudYHNDCBYIqkbYC3A7tExHpJewF7R0St+xMG692amXWaIBqjtZY+8q8UmULPLJ/fCny4bmMHAjOzR1BxnaDO0j8eHxEfA7YCRMSDzOJ+x74KaWZmfWFAsoYqxiUtpRyslPR4oPWt9y04EJiZVQTpmQX72PuBc4G1kk6luKP5tXUbOxCYmVVp8C4WR8T5ki6lmJBGwLERsblu+2QgkLRjjfbNiLi77sEGWTbLJzdnceqCUp/8omWHOZuJ95XJsW4uX5Xclpof+Xfuuy3ZZnLxY5PbblzzrJbrV2UufR2welly20RmfuTJxLbRvhomtvYQ9NeF4KQWN5FNpe3tImmXdtxQ9qtyyf2qjwC71DmQmdmgGKChoakbypYA64CfUXxnPxm4GHhGnZ3kAsGGiNg/11jSZXUOYmY2UAYkEEzdUCbpq8BTI+LK8vmTgA/U3U8uENSJJLWijZnZwOi/1NA69p4KAgARcZWk363bOBn2KvNfHjZ9m6Sjq68xM1tQ1Ki39I8rJH1G0sHl8mngirqN67yT90k6SdK2klZJ+g/gJXPurplZnws1ai195HXA1cCx5XJNua6WOpfGnwe8A7i8fP6+iPjS7PpoZjYgJBgZmKyhkynmJv5WRJwInDiX/dR5tzsABwE3AGuAXSUpou0z4fa1XPTPjSam0k5z+4vM+GRjcmvmaK1NNtLzIzczKZPp3OJM3zNzMafe89btdk62yQ3VrknUqms8eEeyTXPpTsltG+9Pf7bLFqVSZofqv8GQGKj7CD4LvAh4u6Rx4JvAuRHxs9nspM67/XG548OBA4HHAT+cZWfNzAbHgFwjiIiLIuIDEfEc4Ajgl8A7JF0m6XOSjqiznzpnBIdFxC/Lgz4E/B9Jz51zz83M+lyfjf/XEhF3AF8qFyQdABxep23uzuJ9IuJaYIWkFdM2p4u6m5kNsgEsMZHw5Ij42zovzJ0RvB1YT+up0GacAk3SEuB7wOLyOGdGxPsl7Q6cBuwEXAIcFRHjdTprZtYVbbyPQNLhwCcpKjF8JiI+0rad5/0NxTwFM0oGgohYX/4716nQtgCHRMT9ksaAH0j6T4oAc2JEnCbpn4FjgJPmeAwzszZT2yadkTQC/BPwAuAW4KeSzo6Ia9q0/9S9AgLShb+mmfHdln/Zvxl4NsWZwPeBf57pZrIyq2hqCGmsXKbOJF5Trj+F4jZoBwIz6x/tGxo6CPh5RPwCQNJpwMso8vzbYRXw+8Bd09YL+K+6O6kT9v4duA/4f+Xz1wCfB/5wpoZlNLwE2JMiKt4A3B0RE+VLbgHSuYP9JPeLkatM2pxsvTqT0pnNSEzsL1cRNGeUTN8T7zmbStucSG9LfU5K933RxvT/F219qOX6q5c/KdnmvvsfTG7bc4clyW1bE2m2ixfEULJVhZRN4Z5mhaSLK89PjoiTK893Bm6uPL8FeNo8u1h1DrAsIi6fvkHSd+rupE4geFJEPLHy/EJJtaJZREwC+0naHvgasE/djklaT3GNgrVr19ZtZmY2PwGzuEtqc0Ss62BvsiLimMy216S2TVfn75lLJT196omkp1GUN62tnLPgQooiddtLmgpAaygmWW7V5uSIWBcR61aumJ60ZGbWKUEz6i013ApU/5JNfuf1Up1AcADwX5JulHQj8CPgQElXZi5UIGlleSZAOZfmC4ANFAHhVeXLjgbOmnv3zczaK4DJqLfU8FNgL0m7S1oEHAmc3a6+lrOSzfs1dYaGat2Q0MJq4JTyOkEDOD0izimHlU6T9GHgMopbpM3M+ka7KuhExISktwLnUaSPfi4irm7Lzgu/m/uDnOKi8XYz7WTGQBARN82mV5V2VwCPmtimvHp+0Fz2aWbWaQFkSnDNfn8R3wC+0b49PkKd666JDJPfyt1ZfGlETJ8Pc9avMTMbNINSSnCuf6hPlzsjaMspx0LRzNQYHcmljyY0Mr9quRTMVJpoM1P1M1sdNTVBPcDElparJ8e2STbZmvmVWjLZOt1zye3Xpbuw/eOS2655qHU/dlue/iwml6W33T+e/iyWLUp9ioPylWG1RXvPCAZBLhC05ZTDzGzQDFmV/WyJibaccpiZDZKprKFBIunPgS9ExPQ7jGvxfZFmZtM0o97SR1ZR1DE6XdLh0uyq5jkQmJlVRBRDQ3WWfhERfw3sRZGO/1rgekl/J+nxddrPGAgk/bmkHebVSzOzAdKsufSTstDnbeUyQTHN8JmSPjZT2zpnBPM65TAzGzQR9ZZ+IelYSZcAH6OYSvj3IuJNFJUh/vdM7evcUPbXkt4LvBB4HfCPkk4HPhsRN8yr9wMkl+6Zq12eSgXNpojmJFJVc5Pa51JLc31PTUQ/kkgrBRjJ/e9IpL6Or9432WQis7vf29S65NXEdk9Otvn+bek5kPZ/7LbJbam6Mg3/WbTgFDeU9dG3fD07Aq+cnuQTEU1JfzBT41rXCOZzymFmNmjaWGuoW/aYHgQkfR4gIjbM1LjOxDTHAn8CbAY+A/xlRGyV1ACuB945l16bmfWrwTsh4BGn1WWNtwPqNq5TdG5epxxmZoMkCJoDcse4pOOAdwNLJd07tRoYB05ONpymzjWC92e2zXjKYWY2UPrsQnBORBwPHC/p+Ig4bq77ac8MzWZmC0if3SyWJGmfiLgWOEPSowqARsSMcxGAA0FbZOfpTWTzNBelM1RymTJjtD5Wbh7hnNyxRudyepzJLk59FjmTWpTe39jilutHrvtBss2yFc9IbvveTfcktz1v19b1FUdHnTa00BQlJgYkEsDbKab0PaHFtgAOqbMTBwIzs2kGJQ5ExPry3+fPZz8uMWFmVjF1H0Gb5izuCklvmZoauHy+g6Q3123vQGBmVhUw2ay39JHXR8TdU0/KKqSvr9vYQ0NmZhUDemfxiCSVN/9O3UeQvsA2jQOBmdkjxCBdLJ5yLvBlSf9SPn9Dua4WBwIzs4oI2Npn9SNq+CuKL/83lc/Pp6gEUYsDQTtkUjdThds0mS5+NjqSPqObVGJ/mSxGNdMzio4kCsFBejbeGG2dtgkzFL8bW9pyfS6tdOnEfcltW//n6pbrR/fcP9nm/vG5za46eN8LNlfdGhqS9PfASyjuAr4BeF1E3C1pN2ADMDWZ948j4o25fUVEEzipXGbNgcDMbJouDQ2dDxwXEROSPgocR/GXPcANEbHfTDuQdHpEHCHpSlr87RYR6VK8FQ4EZmYVxRlBF44T8c3K0x8Dr5rDbo4t/51X3Tenj5qZVQVMNqPWAqyQdHFlWT/Ho/4p8J+V57tLukzSdyU9J9nViI3lwzdHxE3VBah9H4HPCMzMKoJZ3Sy2OSLWpTZK+hbw2Bab3hMRZ5WveQ/FPC+nlts2ArtExB2SDgC+LmnfiLi3xX6mvIDfDitNeVGLdS05EJiZVQSwtU1jQxFxWG67pNdSDOscOnUPQERsAbaUjy+RdAPwBOBRU/JJehPFX/6Pl3RFZdNyiikra+lYIJC0Fvh3ijmPAzg5Ij4paUfgy8BuwI3AEeVdcGZmvVcODXWapMMpJvZ6XkQ8WFm/ErgzIiYl7QHsBfwisZsvUgwpHQ+8q7L+voi4s25fOnlGMAG8IyIulbQcuETS+cBrgQsi4iOS3kXR+VqnLwtJZFJEI1fBM7W+A1kOc9nnZCM9P3IjEqmbiXmYAUY3XZ/cNj7ROu10ctmKZBvm+CdHJD8LVx9daLp4Z/E/AouB81X8n59KE30u8EFJW4Em8MbUl3pE3CPpfmD/6ZOHzUbHAkF5EWNj+fg+SRuAnYGXAQeXLzsF+A5DGAjMrH91476RiNgzsf4rwFdmsZ9JSddJ2iUifjmXvnTlGkF5g8T+wEXAqsqV7tsoho7MzPrCgNYa2gG4WtJPgAemVkbES+s07nggkLSMIrq9LSLuVWXYIyJCUstPvEzDWg+wdu3aTnfTzAwohgEHsMTEe+fTuKP3EUgaowgCp0bEV8vVt0taXW5fDWxq1TYiTo6IdRGxbuWKzJivmVmbDdp8BBHxXeBaimyh5cCGcl0tHQsEKv70/2zZoU9UNp0NHF0+Pho4q1N9MDObrampKuss/ULSEcBPgD8EjgAuklT7TuVODg09CzgKuFLS5eW6dwMfAU6XdAxwE0Wnzcz6Q0BzUGav/633AAdGxCb4TQrqt4Az6zTuZNbQD0jn1h3aqeP2nTlMKp+rFtpUolpoJuW00UynZyozzVIzUTk194fQSDMzQX2iyqiaE+k2meqoi/bar/VhMvtbsc2y5LY7H8r03YZGcUbQ617MWmMqCJTuYBYjPr6z2Mxsmn4a/6/pXEnnAV8qn/8R8I26jR0IzMwqIoLxPpuQeCYR8ZeSXgk8u1x1ckR8rW57BwIzs4qgOyUm2knS24EvV7IzZ8WBwMysIrpUa6jNlgPflHQnRS23MyLi9rqNPR+Bmdk0s5iPoC9ExN9ExL7AW4DVwHfLEti1+IzAzKwi6K8v+VnaRFG65w7gd+o2ciBog8ikiCpVWTNTcTPZhtwpXHrLnPpHOo01UimskEwRLXbYuh8xujjZZGL7NcltsWR5y/WNh+5Jttl3UXpuj/uWbZvcdtWvH2y5/pk7p9NRbTBFwPjEYF0slvRminuyVgJnAK+PiGvqtncgMDOrGNBrBGsp6rldPpfGDgRmZtMMWiCIiOPm096BwMysYsCvEcyJA4GZWUUETDgQmJkNN58R2OxlsnKSM91msnWikf6xaHI80Sj9i5vLymlmirqlsoayheUytLV15k1jywMt1wPZ7KrJRNbQA0vT81dcf+eW5LatmWJ1T1y5TXKbLSwRDFyJiflyIDAzq/A1AjOzITeg6aPz4hITZmbTdKPEhKQPSLpV0uXl8uLKtuMk/VzSdZJ+f95vaAY+IzAzqyiqj3btGsGJEfHx6gpJTwSOBPYFHgd8S9ITIiI9Y9U8+YzAzKwq6p0NdHD46GXAaRGxJSL+B/g5cFCnDgYOBGZmj9AM2DLRrLUAKyRdXFnWz/Jwb5V0haTPSdqhXLczcHPlNbeU6zrGQ0Mdlp2PN6GZnOoZGiOLZn2cyMxnnJNLcU0eK5OqqomHE40y8yYvSheC03jrtNNtMqf1k5FOA71/PH3mncrOzf6sksnD1s9mOTHN5ohYl9pYloJ+bItN7wFOAj5UHvJDwAnAn86qs23iQGBmVtXGrKGIOKzO6yR9GjinfHorRRG5KWvKdR3joSEzs4qp+wi6kDW0uvL0FcBV5eOzgSMlLZa0O7AX8JN5HWwGPiMwM5umS/cRfEzSfhRDQzcCbwCIiKslnQ5cA0wAb+lkxhA4EJiZPUIETHRhYpqIOCqz7W+Bv+14J0oOBGZmFRHQHLI7ix0IzMweIYhMEceFyIGgDbKpm4lKos3MvL/ZbM853PGYqiJabJzbXMfJ3eWqqo4uad1m8o70/prpdFQmW/9nzf08nrRyp+S2DZkP/urEnMUHPi49Z/Fip2IMrPAZgZnZEBvCoaGO/c1S3im3SdJVlXU7Sjpf0vXlvzvk9mFm1m1BcY9jnWWh6OTJ678Bh09b9y7ggojYC7igfG5m1j8CJiebtZaFomOBICK+B9w5bfXLgFPKx6cAL+/U8c3M5iaIZr1loej2NYJVEbGxfHwbsCr1wrJ403qAtWvXpl5mZtZWxdDQwvmSr6NneQ1R5GclP+2IODki1kXEupUr0nPQmpm1VUAzotayUHT7jOB2SasjYmNZZ2NTl4/fEbnJ5lNXlBqZO8Yjk1oaicnmg0w6au4XNnfFK5E+mqtmqol06qbGW6dg5lJim9um0z1Tn/uiX12ZbLM01QcA7ZLcdN+W1u8rW4qgMbeqr9Z7PiPorLOBo8vHRwNndfn4ZmYz8jWCNpH0JeBgiokbbgHeD3wEOF3SMcBNwBGdOr6Z2VxExILKCKqjY4EgIl6d2HRop45pZtYOC+kegTp8Z7GZWYWLzpmZ2YIa/6/DgcDMrCocCGwO8hOYt5avWJpJH01Noj7HnOaRzLFSaafZdNRMpdPUe57YIX3DYGqCeoDG+EOtu7B0u2Sb5pLlyW233z6e3LZlyC4eDrPAF4vNzIabzwjMzMwXi83MhpxnKDMzG2IRC+uu4TocCMzMpunG0JCkLwN7l0+3B+6OiP0k7QZsAK4rt/04It7Yyb44ELRBI11ENS0zH3AuKycxTS8jmQJn+aJzbf6Fz7yvZAG5zG2cI7dcm9w2eU/ruY6165OSbbYsSU+Kd9dD6bmTU1xXbgGKoDmRziBr32Hij6YeSzoBuKey+YaI2K/jnSg5EJiZVQRBZNKg202SKOquHdK1g07Ts/kIzMz6UkBMTtZaKIpqXlxZ1s/hiM8Bbo+I6yvrdpd0maTvSnpOW95Xhs8IzMweYVZnBJsjYl1qo6RvAY9tsek9ETFVhv/VwJcq2zYCu0TEHZIOAL4uad+IuLdup2bLgcDMrCraNzQUEYfltksaBV4JHFBpswXYUj6+RNINwBOAi9vSqRYcCMzMpuniNYLDgGsj4papFZJWAndGxKSkPYC9gF90shMOBGZmFdGlrKHSkTxyWAjgucAHJW0FmsAbI+LOTnbCgWDAjCbSFbt6+0sm3TNbTC/RThMPp9s8lC46N/LY3Vuuf3D79NzDiya3JLc9fselyW13PtT6feXSdm1QBc0unRFExGtbrPsK8JWudKDkQGBmVtXGawSDwoHAzKwi6Oo1gr7gQGBmVhUxdY/A0HAgMDOr6u7F4r7gQGBm9gi+RmBmNtSKawSeqtLaKVGNM5fumZsDeSSR/hiji+e0P2WyH5X4qyg/33LmVyqRPhpj2ySbbN330OS20Xs2tly/7a9+lmwzsWKP5LYHt6Y/jB2Xtn5fk5lyxWMjTi0dSM4aMjMzBwIzs2EW3buhrF84EJiZVUQEza3OGjIzG2K+RmBmNvQcCMzMhtkQZg0p2j15eQdI+jVwU6/7UVoBbO51Jyrcn7R+6gv0V3/6qS/Q3v7sGhEr59pY0rllf+rYHBGHz/VY/WIgAkE/kXRxbmq6bnN/0vqpL9Bf/emnvkD/9WfYePJ6M7Mh50BgZjbkHAhm7+Red2Aa9yetn/oC/dWffuoL9F9/hoqvEZiZDTmfEZiZDTkHAjOzIedAkCHpc5I2Sbqqsm5HSedLur78d4cu9WWtpAslXSPpaknH9rg/SyT9RNLPyv78Tbl+d0kXSfq5pC9LWtSN/pTHHpF0maRz+qAvN0q6UtLlki4u1/XkZ1Uee3tJZ0q6VtIGSc/oRX8k7V1+JlPLvZLe1svPxhwIZvJvwPSbRd4FXBARewEXlM+7YQJ4R0Q8EXg68BZJT+xhf7YAh0TEU4D9gMMlPR34KHBiROwJ3AUc06X+ABwLbKg872VfAJ4fEftV8uN79bMC+CRwbkTsAzyF4nPqen8i4rryM9kPOAB4EPhaL/piFRHhJbMAuwFXVZ5fB6wuH68GrutRv84CXtAP/QG2AS4FnkZxd+houf4ZwHld6sMaii+QQ4BzAPWqL+XxbgRWTFvXk58VsB3wP5TJIb3uT+X4LwR+2A99GfbFZwSztyoipqbGug1Y1e0OSNoN2B+4qJf9KYdiLgc2AecDNwB3R8TUFGa3ADt3qTv/ALwTmJoGbace9gWKSei+KekSSevLdb36We0O/Br413Lo7DOStu1hf6YcCXypfNzrvgw1B4J5iOLPl67m30paBnwFeFtE3NvL/kTEZBSn+GuAg4B9unXsKkl/AGyKiEt6cfyEZ0fEU4EXUQzjPbe6scs/q1HgqcBJEbE/8ADThl66/btTXq95KXDG9G29+H817BwIZu92SasByn83devAksYogsCpEfHVXvdnSkTcDVxIMfyyvaSpqrZrgFu70IVnAS+VdCNwGsXw0Cd71BcAIuLW8t9NFGPgB9G7n9UtwC0RcVH5/EyKwNDL350XAZdGxO3l857/Hg8zB4LZOxs4unx8NMVYfcdJEvBZYENEfKIP+rNS0vbl46UU1ys2UASEV3WzPxFxXESsiYjdKIYbvh0Rf9yLvgBI2lbS8qnHFGPhV9Gjn1VE3AbcLGnvctWhwDW96k/p1fx2WIge98V6fZGinxeKX9SNwFaKv6qOoRh7vgC4HvgWsGOX+vJsitPlK4DLy+XFPezPk4HLyv5cBbyvXL8H8BPg5xSn/Yu7/DM7GDinl30pj/uzcrkaeE+5vic/q/LY+wEXlz+vrwM79PB3Z1vgDmC7yrqefTZewiUmzMyGnYeGzMyGnAOBmdmQcyAwMxtyDgRmZkPOgcDMbMg5EFhfk7SbpIfKUhYzvfYzZSG+1PalZcXLcUkr2tpRswHm9FHra2VdpXMi4klt3OeNwLqI2NyufZoNMp8RWM9IOlDSFeXcBtuW8xpkv/DLM4RrJZ1a1tU/U9I25bbvSFonadeyrv0KSQ1J35f0wu68K7PB40BgPRMRP6UoLfBh4GPAFyLiqnwrAPYGPhURvwvcC7x52n5vopiL4CTgHcA1EfHNdvbdbCFxILBe+yBFnaJ1FMGgjpsj4ofl4y9QlN94hIj4DPAY4I3AX7Shn2YLlgOB9dpOwDJgObCkZpvpF7YedaGrHC5aUz5dNufemQ0BBwLrtX8B3gucSjGcU8cukp5RPn4N8IMWr/louc/3AZ+ebyfNFjIHAusZSX8CbI2ILwIfAQ6UdEiNptdRTPaygaKK5knT9vs84EDgoxFxKjAu6XXt7b3ZwuH0Uetr09NH25FO6vRRs0fyGYH1u0lguzo3lM1k6oYyYIzfzm1sNvR8RmBmNuR8RmBmNuQcCMzMhpwDgZnZkHMgMDMbcg4EZmZD7v8DXOdtk4RLfFEAAAAASUVORK5CYII=\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"d.piv.showscal(colorbar=True)"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 291
},
"id": "wZBrnSoU2xc7",
"outputId": "5a5160fe-ba2d-4935-d99a-75b93399ab04"
},
"execution_count": 55,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
],
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAESCAYAAAA48DgcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAAr+klEQVR4nO3dd7xcVb3+8c8DBEI1QACRInCpoUOkCGioopcLXEUviAj8UC6KCgoqiIANRax4vZZIVQNIlaICkSJyqQECafSaUEKkBEJL4Pn9sdfAcDhlzpnMmXOS5/16zevMXnutPd8J4Xyz9lp7LdkmIiKiGQu0O4CIiBj8kkwiIqJpSSYREdG0JJOIiGhakklERDQtySQiIpqWZBIREU1LMokBT9JqkizpRUkH15V/TtJTpXzZNse4U4njDUk7tTOWiHZIMoluSXq4Vb8cJX1S0iOSZkn6s6RlemgyzPbo0nYI8FNgF9tL2P5XK2JslO2/214CeLSdcUS0S5JJtIWk9YHfAvsBKwAvAb/qxSVWAIYCk7q4/kLNxhgRjUsyiS5J+gOwKnBpuYXztbl4+X2BS21fZ/tF4Fjgo5KWbCCutYF7yuFzkq4u5ZZ0qKT7gPtK2cmSHpM0U9Jtkraru863JJ0n6Y+SXpA0QdLako6WNL2026Wu/rsknSrpCUnTJH1P0oJz748kYvBKMoku2d6P6rbNf5RbSSd1rCNpVUnPdfP6ZBeXXx+4s+6zHgBeA9ZuIK57S3uobn3tUHd6T2BLYEQ5vhXYBFgGOAs4T9LQuvr/AfwBWBq4A7iC6v+LlYDvUPWeas4A5gBrApsCuwCf6SneiPlBkkk0xfajtod18zqri6ZLAM93KHse6LFn0oMf2H7G9sslvj/a/pftObZ/AiwCrFNX/5+2r7A9BzgPWA440fZs4BxgNUnDJK0AfAQ43PYs29OBnwF7NxlvxDwh95WjXV4ElupQthTwQpPXfaz+QNKRwEHAewCXzxheV+WpuvcvAzNsv153DFXiew8wBHhCUq3+Ah0/L2J+lWQSPel2jwJJqwKTu6ny37bHdFI+Cdi47jprUPUa7u1LkHXejLeMj3wN2BGYZPsNSc8C6qpxNx4DXgWGl15MRNRJMomePAWs0dVJ249S/cu9t8YAN5Zf+LdTjU9caLvZnkm9JanGOJ4GFpJ0FO/sDTXE9hOSrgR+IulYqp7V6sDKtv8xtwKOGKwyZhI9+QHwzTKYfuTcuqjtScAhVEllOtUv/s/PresXVwCXU/V2HgFeobnbUp8GFqbqiT0LnA+s2GSMMQ+TdFqZGTixruzjkiaVB1xH1pUvLOn0MqvwTkmj6s5dK+keSePLa/lSfoCkp+vKP1PKt68rGy/pFUl7tvS7ZqfFGOgkvZdqKvArwFdt/67NIb2DpB2BC6hu1X3E9jVtDikGAEkfoOrF/t72BqVsPeANqpmCR9oeV8oPBUbaPrAki78B7yu3Z6+tr1t3/QNKmy90E8MywP1UveiX5vZ3rMltrhjwbD9C9YDigGX7KmBYu+OIgcX2dZJW61A2BaBuIkfNCODqUme6pOeAkcAtTYaxF/C3ViYSyG2uiIiB4k5gd0kLSVod2BxYpe786eWW1bF6eyb6mKS7JJ0vaRXeaW/g7BbGDaRnEhHxNksNG+E5s2c1VPfllx6dRHX7tWZ0bf24PjgNWA8YRzXGdwNQm6a+r+1pZYWIC6iWIfo9cClwtu1XJf03cCbw5kO8klYENqQaP2ypJJOIiDpzZs9i3Q2+3lDdO2459BXbI3uu2bMy5fzLtWNJN1CmytueVn6+IOksYAuqcZj6BU5PATquUvEJ4KLyEG5LDYpksqQW9HIMaXcYMUAtvkR1t3bWi2+0OZIYCB7i1Rm2l2t3HL0laTGqSVGzJO0MzLE9uSxaOsz2jLJa9m7A30ubFW0/US6xOzClw2X3AY7uj/gHRTJZjiF8j/e2O4wYgLYeNZSVdl8TgGmX3M+N177SQ4uY1+3LvY+0O4YaSWcDo4DhkqYCxwPPAP9DtXTPXySNt/0hYHngCklvANOobmVBNUPwipJIFqRKJLUZjV+StDvV81TPAAfUffZqVGMu/fIc1KBIJhGdqSWSIbtsDutvw0r8gq1JQomBw/Y+XZy6qJO6D/P2deNq5bOoBuM7u/7RdNHzKNdbqcFQm5bZXDEo1SeSCe8+jAn/GsmQL3+JlXZfk61HDehZxBHzpPRMYtDpmEiumvxkdWLESDb88pfSQ4logySTGFQ6SyT/uObBtyokoUS0RZJJDCr1YyRX/bNKJI/+5T6gbpRxxEg23GVzVoIklIh+kjGTGFSmXXL/m+93HPHuTutsuOy4TssjonXSM4lB5cZrX4Ev/InVfgkbrg/fOfT9HFfOfXD7NTh8u6kw6f+YfeVtmSoc0Y+STGLQ6SyhXDX5yaqnMulPSSQRbZBkEoNSx4Sy4XYkkUS0UZJJDFr1CQVIIolooySTGNRqCWWl3ddMIolooySTGPRuvPYVuHZizxUjomUyNTgiIpqWZBIREU1LMomIiKYlmURERNOSTCIiWkTSaZKmS5pYV7aMpLGS7is/l647N0rSeEmTJP2jlK0i6RpJk0v5YXX1/1Tqj5f0sKTxpXyLuvI7Jf1nq79rkklEROucAezaoewo4CrbawFXlWMkDQN+Bexue33g46X+HOAI2yOArYBDJY0AsP1ftjexvQlwAXBhaTMRGFnKdwV+W7b/bZkkk4iIFrF9HdV2uvX2AM4s788E9izvPwlcaPvR0nZ6+fmE7dvL+xeo9nl/2w6KkgR8Aji71HvJ9pxyeijgufetOpdkEhHRv1aw/UR5/ySwQnm/NrC0pGsl3Sbp0x0bln3dNwVu7nBqO+Ap2/fV1d1S0iRgAnBIXXJpiTy0GBFRZ+Fhi7Dqv6/VUN07bmG4pPo9D0bbHt3oZ9m2pFqvYSGqvd53BBYFbpR0k+17ASQtQXUr63DbMztcah9Kr6Tu2jcD60taDzhT0t9st2yJiCSTiIi+m2F7ZC/bPCVpRdtPSFoRmF7KpwL/sj0LmCXpOmBj4F5JQ6gSyRjbF9ZfrIyFfJQqEb2D7SmSXgQ2AFq22U9uc0VE9K9LgP3L+/2Bi8v7i4FtJS0kaTFgS2BKGQ85FZhi+6edXG8n4G7bU2sFklavDbhLei+wLvBwK75MTXomEREtIulsYBQwXNJU4HjgROBcSQcBj1ANnNd6EJcDdwFvAKfYnihpW2A/YEJt6i/wDdt/Le/3psMtLmBb4ChJs8u1Pm97Rou+JpBkEhHRMrb36eLUjl3U/xHwow5l1wPq5jMO6KTsD8AfGg50LshtroiIaFqSSURENC3JJCIimtbSMRNJDwMvAK8Dc2yPlLQM8CdgNarZBZ+w/Wwr44iIiNbqj57J9mXtmNpc7E7XpYmIiMGrHbe5ulqXJiIiBqlWJxMDV5Z1Zg4uZV2tSxMREYNUq58z2db2NEnLA2Ml3V1/ssO6NG9Tks/BAMPzOExExIDW0p6J7Wnl53TgImALyro0AB3WpenYdrTtkbZHLsmCrQwzIiKa1LJkImlxSUvW3gO7UG3Y0tW6NBERMUi18v7RCsBF1RplLAScZftySbfSybo0ERExeLUsmdh+kGr55I7l/6KLdWkiImJwyhPwERHRtCSTiIhoWpJJREQ0LckkIqKFJB0maaKkSZIOL2U/knS3pLskXSRpWClfWNLpkiZIulPSqLrr7FPK75J0uaThpXxjSTeWc5dKWqr/v2WSSUREy0jaAPgs1TN2GwO7SVoTGAtsYHsj4F7g6NLkswC2NwR2Bn4iaYGyBe/JVGsdbkS1G+MXSptTgKNKm4uAr/bLl+sgySQionXWA262/ZLtOcA/gI/avrIcA9wErFzejwCuhjcf9n4OGEm106KAxcue8EsBj5c2awPXlfdjgY+19Bt1IeuURES0zkTgBEnLAi8DHwHGdajz/6i25QC4E9i97B2/CrA5sIrtWyR9DpgAzALuAw4tbSZRLaD7Z+DjpV3DJG3WQLXZtid0VyHJJCKizpJLLsIHt1+joboXH89wSfXJYbTt0bUD21Mk/RC4kioJjKfa3wkASccAc4Axpeg0qt7MOKqHum8AXpc0BPgcsCnwIPA/VLfGvkeVjH4h6ViqFUZe6+VX/gdwK93sMw+sTrUHVZeSTCIi+m5G3V5NnbJ9KnAqgKTvA1PL+wOA3YAdbbvUnQN8udZW0g1UYyqblPMPlPJzKXtB2b6barkqJK0N/Hsvv8OttnforoKkq3u6SMZMIiJaqKyajqRVgY8CZ0naFfgasLvtl+rqLlbWMkTSzlQ71E4GpgEjJC1Xqu4MTOlw/QWAbwK/6U18PSWSRuukZxIR0VoXlDGT2cChtp+T9EtgEaqtOQBusn0IsDxwhaQ3qBLIfgC2H5f0beA6SbOpboEdUK6/j6Ta+MmFwOl9CVLSENuzO5QNtz2jkfZJJhERLWR7u07K1uyi7sPAOl2c+w2d9Dpsn0w1bbhPJG0P/AEYKul24OASB1RjPY0M0Oc2V0TEfO4k4EO2hwOjqXpLW5Vz3Q3Kv016JhER87eFbU8CsH2+pCnAhZK+TrX1ekOSTGKesPWood2ev/HaV/opkohBZ7akd9t+EsD2JEk7ApcB/9boRXKbKwa9nhJJo3Ui5lNHUW1m+CbbU4EPAic2epH0TCIi5mO2/95F+fPACY1eJ8kkBr36W1hd9UBymyuie5J2tj22r+2TTGKekqQR0XuSPg3sQ7VQZJ8kmUREzMfK+mA70ftlWN4mySQiYv52HLBe/bIufZHZXBER87cDgfNrOzf2VXomERHzMdtnSXqSaj+Ubft6nSSTiIj5nO2rJU1v5hq5zRUREdie2Ez79EwiIgJJI4HtgPdQbTE8ERhr+9lG2qdnEhExH5N0YFl6/mhgUeAeYDrV+MnfJZ1ZNvbqVnomEREtImkd4E91RWtQTcW9hmpvkiWAh4F9bc+UtC/w1br6GwGb2R4vaXPgDKpf+H8FDrNtSd8CPgs8Xdp8w/ZfexHmYsA2tl/u4jtsAqwFPNrdRVreM5G0oKQ7JF1WjleXdLOk+yX9SdLCrY4hIqIdbN9jexPbmwCbAy8BFwGnAEfZ3rAcf7XUH1NXfz/gIdvjy+V+TZU01iqvXes+6me1dr1MJNj+364SSTk/3vZVPV2nP25zHUbZq7j4IdUXXxN4FjioH2KIiGi3HYEHbD8CrA1cV8rHAh/rpP4+wDkAklYElrJ9k20Dvwf2nJvBSTpJ0lKShki6StLTkj7VaPuWJhNJK1M9on9KORawA3B+qXImc/kPJCJigNobOLu8nwTsUd5/HFilk/r/VVd/JWBq3bmppazmC5LuknSapKX7GN8utmcCu1HdeluTt99y61arx0x+DnwNWLIcLws8Z3tOOe74B/ImSQcDBwMMz9BORPSTpYYOYccR7260+nBJ4+qOR9se3bFSuZ2/O9UgN8D/A34h6VjgEuC1DvW3BF5qcLrur4HvUu2K+F3gJ+X6vVX7RfvvwHm2n6/+/d+7xnOdpN2A6bZvkzSqt+3Lf5DRAGtoaMNbR0ZE9KMZtkc2UO/DwO22nwKwfTewC4CktXnnIov1vRiAacDKdccrlzJq1yzX+h3VDol9cZmku6mmBX9O0nJAw8twt/I21zbA7pIeprrvtwNwMjBMUi2JvfkHEhExD9uHuuQgafnycwHgm1Qzu6gr+wRlvATA9hPATElbleGCTwMXl/or1n3Of1I9H9IwSe8pn3EU8H5gpO3ZVJMF9uiubb2WJRPbR9te2fZqVFn2atv7Uk2J26tU25/yBxIRMS+StDiwM3BhXfE+ku4F7gYeB06vO/cB4DHbD3a41Oepxp/vBx4A/lbKT5I0QdJdwPbAl3sZ4imSbpJ0ItVUZAHYnlXbF74R7RiM+DpwjqTvAXcAp7YhhoiIfmF7FtV4cX3ZyVR3ajqrfy2wVSfl44ANOinfr8n4PiJpKDCKqmfzY0mPApcDl9vu9vmSmn5JJuUP59ry/kFgi/743IiI6JntVyjJA6rnAanGeX4p6d22e/ydnWlSERHxJklLAc9TjdmcA7zYSLskk4iIQNJ/A9+mmsFVm0Fr22s00j7JJCIiAI4ENrA9oy+Ns2pwRERANUOsz/vAp2cSERFQPZ1/g6SbgVdrhba/1EjjJJOIiAD4LXA1MAF4o7eNk0wiIgJgiO2v9LVxxkwiIgLgb5IOlrSipGVqr0Ybp2cSERFQrR8Gb61sDNUU4UwNjoiIxthevZn2uc0VETEfk7RtD+eXkvSONcE6Ss8kImL+9jFJJ1Gty3Ub8DQwlGqnxe2B9wJH9HSRJJOIiPmY7S+XgfaPUW0hvCLVBllTgN/avr6R6ySZRETM52w/A/yuvPokYyYREfOxssV603WSTCIiWkjSMEnnS7pb0hRJW5fyL5aySWXMor7NqpJelHRkOV5F0jWSJpf6h9XV/Xgpe0NSI/vRd/QjSZtK2qyrF/D9ni7S5W2uBh9WecP2c70IOiJifnMy1Y6Fe0laGFhM0vZU+6tvbPvV2p7wdX7KW9vyAswBjrB9u6QlgdskjbU9mWrP949SLYfSF0+Vz+vOfT1dpLsxk8fLS93UWRBYtacPiYiYH0l6F9We7gcA2H4NeE3S54ATbb9ayqfXtdkTeAiYVSuz/QTwRHn/gqQpwErAZNtTSrs+xWh7VJ8adtBdMplie9PuGku6Y24EERExUCy60Cw2XHZco9WHS6qvPNr26Lrj1amm2p4uaWOqqbeHAWsD20k6gWozqiNt3yppCeDrwM5U+4u8g6TVgE2Bmxv/Vq3XXTLZuoH2jdSJiJhXzbDd3TjFQsBmwBdt3yzpZOCoUr4MsBXwPuBcSWsA3wJ+ZvvFznoaJdlcABxue+Zc/SZN6nIAvmwwj6SdOp6TtH99nYiI6NRUYKrtWi/ifKrkMhW40JVbqJZ8Hw5sCZwk6WHgcOAbkr4AIGkIVSIZY/vCfv0WDWhkNtdxkn4taXFJK0i6FPiPVgcWETHY2X4SeEzSOqVoR2Ay8Geqp8uRtDawMFUvZzvbq9leDfg58H3bv1TVTTmVavihp8HyPpF0m6RDJS3dl/aNJJMPUm3nOB64HjjL9l59+bCIiPnQF4Exku4CNqGaZnsasIakicA5wP623c01tgH2A3aQNL68PgIg6T8lTaUadviLpCv6GOd/Ae8BbpV0jqQPqRej+o08Ab80sAVVQlkZeK8k9fDFIyICsD0e6Gxc5VM9tPtW3fvr6WJmre2LgIv6HuGb17kfOEbSscBuVAnvdUmnAyeXp+S71EjP5CaqOdK7Ug0UvQf4v+bCjoiIgUbSRsBPgB9Rjc98HJhJtZ1vtxrpmexk+1EA2y8DX5L0gb6HGxERA42k24DnqMZmjqo9AwPcLGmbntp39wT8urbvpppHPbzD6Rf7GG9ERAxMH7f9YGcnbH+0p8bd9Uy+AhxM1eV5x7WBHbq7sKShwHXAIuVzzrd9vKTVqQaclqV6gGe/8lRoRES0z56djLc/D9xWxn261WUysX1w+bl9HwN7FdihPHwzBLhe0t+oktTPbJ8j6TfAQcCv+/gZERExd4wsr0vL8W7AXcAhks6zfVKXLWlgAF7SUElfkXShpAskHV56Hd0qD+PUbocNKa9aj+b8Un4msGdP14qIiJZbGdjM9hG2jwA2B5anbm2x7jQym+v3wPrA/wC/LO//0EhkkhaUNB6YDoylml78nO05pcpUqsXKIiKivZanuqNUMxtYoUy8erXzJm9pZDbXBrZH1B1fI2lyI5HZfh3YRNIwqnnQ6zbSDkDSwVRjNgzPhpAREa02hmrm1sXl+D+AsyQtTvXUfrca+S19u6StbN8EIGlLoOElNQFsPyfpGqonNIdJWqj0TlYGpnXRZjQwGmANDc0DkhERLVKedD+Dag+V2jTgQ2zXftfv29M1GkkmmwM3SHq0HK8K3CNpAtXQyEZdBLccMLskkkWpllT+IXANsBdlCQHg4s7aR0RE/7BtSX+1vSG97CzUNJJMdu3LhYEVgTMlLUg1NnOu7cvKLbJzJH0PuIPqAZmIiGiv2yW9z/atfWncYzKx/UhfLmz7LqoNXDqWP0i11ldERAwcWwKfKsvfz6JaC6zLu08ddfcE/O22N+uucSN1IiJiUPhQM42765msV5ZM7oqAdzXz4RERMTDYfkTStsBatk8v495LNNq+u2TSyDTe1xv9oIiIGLgkHU/1BPw6wOlUD5r/kbdmd3Wru+VU+jRWEhERbyljEC9Q/eN7ju2Rkr4L7EG1Xe904ADbj0tal+oX+WbAMbZ/XK7R6VqH5dw/gSXLxy0P3GJ7zz6E+p9U49y3A5R4luy+yVvyNGBEROttb3tG3fGPbB8LIOlLwHHAIcAzwJd45zJTna51aPsm29vVKkm6gL4/bvFamSLscq3Fe9O4keVUIiJiLrI9s+5wcap1C7E9vUzNnd2hfldrHb5J0lJUax/+uY9hnSvpt1QPln8W+Dvwu0Yb99gzkfRF4I+2n+1jgBERg8fLL8KkhjeTHS6p/iG/0WX1jnoGriz/4v9t7bykE4BPUy3z3uPq7OWZvduANYH/tX1zhyp7Ald1SFQNs/1jSTtT7ay4DnCc7bGNtm/kNtcKVBvM3061J/AV2f89IgKAGbY729+93ra2p0laHhgr6W7b19k+hmrP9aOBLwDHd3eRjmsdStrA9sS6KvsAp/T9q0BJHg0nkHo93uay/U1gLaon1Q8A7pP0fUn/1pcPjIiYn9ieVn5Op1rwtuND22OAj/Xies9RLUv15uokZTfcLYC/9DVOSR+VdJ+k5yXNlPSCpIZ7OQ2NmZSeyJPlNQdYGjhfUrebpUREzM8kLV6bEVUGtHcBJkpaq67aHsDdPVxnudIjoW6tw/o2ewGX2X6liXBPAna3/S7bS9le0vZSjTZuZMzkMKr7ejOoulBftT1b0gLAfcDX+hh4RMS8bgWqW1JQ/b49y/blZaPBdaimBj9CNZMLSe+mWmhxKeANSYcDI+hircO6z9kbOLHJWJ+yPaWvjRsZM1kG+GjH505svyFpt75+cETEvK6sRbhxJ+Wd3tay/STV1hwddbrWYV27UX0Msd44SX+img325mZYti9spHEjCz12OSjUTBaLiIgBZSngJapbcTUG5k4yiYiIeZ/tA5tpn4cWIyICSWtLukrSxHK8kaRvNto+ySQiIqB62v1oytP3ZU+qvRttnGQSEREAi9m+pUPZnEYbJ5lERATAjPIwem2hx72AJxptnAH4iIgAOBQYDawraRrwELBvo42TTCIiovZMzE7lSf0FbL9Qf17S/rbP7Kp9bnNFRMSbbM/qmEiKw7prl55JzNO2HjX0zfc3XtvMskUR8z11dzI9k5hv1CeWiOi1brceSTKJeVZnySMJJaLP0jOJqJeEEtEn3W4/mWQSERGUTQ+H1R0vLel7tWPbX+iufZJJREQAfLjs4giA7WeBjzTaOMkkIqJFJA2VdIukOyVNkvTtUr6DpNslTZR0pqSFSvm7JF1aV//Aumu9Lml8eV1SVz5G0j3lWqdJGtLHcBeUtEjddRcFFumm/tu0LJlIWkXSNZImlz+Uw0r5MpLGlr2Gx0paulUxRHQmU4SjH70K7GB7Y2ATYFdJ7wfOBPa2vQHVTov7l/qHApNL/VHATyQtXM69bHuT8tq97jPGAOsCGwKLAp/pY6xjgKskHSTpIGBsibMhreyZzAGOsD0C2Ao4VNII4CjgKttrAVeV44iIeY4rL5bDIeX1OvCa7XtL+VigtvOigSVV7fO7BPAMPSy2aPuv5XMM3ELnOzU2EusPgROA9crru7ZParR9y5KJ7Sds317evwBMAVYC9uCtbHcmsGerYoiIaDdJC0oaD0ynShy3AAtJGlmq7AWsUt7/kuoX+ePABOAw22+Uc0MljZN0k6Q9O/mcIcB+wOV9jdX232wfWV5X9KZtvzwBL2k1qv2LbwZWsF1bifJJYIUu2hwMHAwwPA/qx1ySW1zRE898idlX3tZo9eGSxtUdj7Y9+m3Xs18HNikzpS4C1qfaJ+RnZYziSqreCsCHgPHADsC/AWMl/dP2TOC9tqdJWgO4WtIE2w/UfdSvgOts/7MXXxdJ19veVtILvP3BRFXhe6lGrtPy39KSlgAuAA63PbPqvVVsW1KnT1WW/yCjAdbQ0G6fvIyIaJMZtkf2XA1sPyfpGmBX2z8GtgOQtAuwdql2IHBiuWV1v6SHqMZDbrE9rVznQUnXUv0D/YFyjeOB5YD/7u0XsL1t+blkb9vWa+lsrtLtugAYY7u2Kf1TklYs51ek6vpFzHUdeyHplUR/k7Rc7dmNMjtqZ+BuScuXskWArwO/KU0eBXYs51YA1gEeLM98LFLKhwPbAJPL8WeoejT71N0S60usf2ikrCutnM0l4FRgiu2f1p26hLdmLuwPXNyqGCJqkkiiTVYErpF0F3ArMNb2ZcBXJU0B7gIutX11qf9d4P2SJlBNUPq67RlU4yjjJN0JXEPVe5lc2vyGarjgxjJt+Lg+xrp+/UGZrrx5o41beZtrG6rBoAll8AngG8CJwLll6tkjwCdaGEPM55JEop3KPuqbdlL+VeCrnZQ/DuzSSfkNVFN/O/uMpn6PSzqa6nfzopJm8tYaXK9Rhhoa0bJkYvt6ul4YbMdWfW5ERDTO9g+AH0j6ge2j+3qdPAEfEREAx0j6lKRj4c0Hz7dotHGSSUREAPwvsDXwyXL8YilrSB7giIgIgC1tbybpDqgWeqxbyqVH6ZlERATAbEkLUh5clLQc0PBU4ySTiIgA+AXVE/rLSzoBuB74fqONc5srIiKwPUbSbVSzbQXsaXtKo+2TTCIiouY+YCYlN0ha1fajjTRMMomICCR9ETgeeIpq4UlRjZ9s1Ej7JJOIiAA4DFjH9r/60jgD8BERAfAY8HxfG6dnEhERAA8C10r6C9V2wwB0WKi3S0kmEREB1fL3jwILl1evJJlERAS2v91M+ySTiIj5mKSf2z5c0qW8fdteAGzv3sh1kkwiIlpE0irA76k2rzLVHvEnl3NfBA6lmob7F9tfk7QaMAW4p1ziJtuHlPoLA78ERlEtc3KM7QskrQqcCQwDFgSOsv3XXoRZ203xx338mkCSSUREK80BjrB9u6QlgdskjaVKLnsAG9t+tbaNb/GA7U06udYxwHTba0taAFimlH8TONf2ryWNAP4KrNZogLZvK29frHsPgKTdGr1OpgZHRLSI7Sds317ev0DV61gJ+BzV1ruvlnPTG7jc/wN+UOq/UbbzharHs1R5/y7g8T6G+ztJG9QOJO0DHNto4ySTiIh+UG5hbQrcDKwNbCfpZkn/kPS+uqqrS7qjlG9X2g4r574r6XZJ50laoZR9C/iUpKlUvZIv9jHEvYDfS1pX0meBz9PJFsJdyW2uiIg6rz33KtMuub/R6sMljas7Hm37HfumS1oCuAA43PZMSQtR3abaCngfcK6kNYAngFVt/0vS5sCfJa1P9bt6ZeAG21+R9BWqMY79gH2AM2z/RNLWwB8kbWC74eXjAWw/KGlv4M9UU4R3sf1yo+2TTCIi+m6G7ZHdVZA0hCqRjLF9YSmeClxo28Atkt4Ahtt+mvLAoO3bJD1A1Yu5DXgJqLU/DziovD8I2LW0uVHSUGA40MitMyRN4O2zuJahGsi/WRK2szZXREQ7SRJwKjClw5Pkfwa2B66RtDbVQ4IzyoZUz9h+vfRU1gIetO0ydXcUcDXVMvGTy7UeLcdnSFoPGAo83YswGx5k706SSURE62xDdStqgqTxpewbwGnAaZImAq8B+5eE8QHgO5JmU03/PcT2M6Xd16luYf2cKlkcWMqPoBo8/zJVD+OA0uNpiO1HmvmCNUkmEREtYvt6qqXcO/OpTupfQHVLrLNrPQJ8oJPyyVRJq60ymysiIpqWZBIREU1LMomIiKYlmURERNNalkwknSZpepmtUCtbRtJYSfeVn0u36vMjIqL/tLJncgblQZo6RwFX2V4LuKocR0TEINeyZGL7OuCZDsV7UC2VTPm5Z6s+PyIi+k9/j5msYPuJ8v5JqmWYIyJikGvbAHx5QrPLpzQlHSxpnKRxL/B6P0YWERG91d/J5ClJKwKUn10uRGZ7tO2RtkcuyYL9FmBERPRefyeTS4D9y/v9gYv7+fMjIqIFWjk1+GzgRmAdSVMlHQScCOws6T5gp3IcERGDXMsWerS9TxendmzVZ0ZERHvkCfiIiGhakklERDQtySQiokW6WFbqu5LukjRe0pWS3lPK9y3lEyTdIGnjujZfljRJ0kRJZ5etees/5xeSXuy/b/ZOSSYREa1zBu9cVupHtjeyvQlwGXBcKX8I+KDtDYHvAqMBJK0EfAkYaXsDqv3Z965dTNJIoO3rHCaZRES0SGfLStmeWXe4OOXhbds32H62lN8ErFxXbyFgUUkLAYsBjwNIWhD4EfC1lnyBXsi2vRER/UzSCcCngeeB7TupchDwNwDb0yT9GHgUeBm40vaVpd4XgEtsPyF1tTtw/0gyiejE1qOGcuO1r7Q7jGiDWS++0Zv/9sMljas7Hm17dE+NbB8DHCPpaKqEcHztnKTtqZLJtuV4aapFclcHngPOk/Qp4Grg48CoRoNtpSSTiDpbjxrac6WIt8ywPbKJ9mOAv1KSiaSNgFOAD9v+V6mzE/CQ7adLnQuB9wPPAmsC95deyWKS7re9ZhPx9FmSScz3kkCiP0lay/Z95XAP4O5SvipwIbCf7XvrmjwKbCVpMarbXDsC42z/BXh33XVfbFcigSSTmI8liUSrlWWlRlHdDptK1QP5iKR1gDeAR4BDSvXjgGWBX5Wexpyy2O3Nks4HbgfmAHdQZnoNJEkmMd9JEon+0sWyUqd2UfczwGe6OHc8deMqXdRZotcBzkVJJjFfSAKJaK0kk5inJYlE9I88tBjzrGYTSRJRROOSTCIiomlJJhFdyEOLEY3LmEnMs2689pVe3apK8ojouySTmK8lgUTMHUkmMd9I4ohonSSTmKclgUT0jySTGPQ6joskgUT0v8zmioiIpqVnEoNeeiIR7ZeeSURENC3JJCIimpZkEhERTUsyiYhoIUm7SrpH0v2Sjmp3PK2SZBIR0SKSFgT+F/gwMALYR9KI9kbVGoNiNtdDvDpjX+59pM1hDAdmtDmGegMtHhh4MQ20eGDgxTTQ4oHmY3rv3ApkLtgCuN/2gwCSzqHa931yW6NqgUGRTGwv1+4YJI2zPbLdcdQMtHhg4MU00OKBgRfTQIsHBmZMTVgJeKzueCqwZZtiaalBkUwiIvrLQ7x6xb7cO7zB6kMljas7Hm17dCviGuiSTCIi6tjedS5ebhqwSt3xyqVsnpMB+MYNtH9tDLR4YODFNNDigYEX00CLBwZmTH11K7CWpNUlLQzsDVzS5phaQrbbHUNExDxL0keAnwMLAqfZPqG9EbVGkklERDQtt7k6Iek0SdMlTawrW0bSWEn3lZ9L92M8q0i6RtJkSZMkHdbOmCQNlXSLpDtLPN8u5atLurk8nPWn0q3vV5IWlHSHpMvaHZOkhyVNkDS+Nkjbzr9H5fOHSTpf0t2Spkjauo1/j9Ypfza110xJh7f7zyj6Jsmkc2cAHQfhjgKusr0WcFU57i9zgCNsjwC2Ag4tDz61K6ZXgR1sbwxsAuwqaSvgh8DPbK8JPAsc1E/x1DsMmFJ33O6Ytre9Sd1U13b+PQI4Gbjc9rrAxlR/Vm2JyfY95c9mE2Bz4CXgonbFE02ynVcnL2A1YGLd8T3AiuX9isA9bYztYmDngRATsBhwO9Xc+RnAQqV8a+CKfo5lZapfPjsAlwFqZ0zAw8DwDmVt+28GvAt4iHJ7eyDEVBfDLsD/DZR48ur9Kz2Txq1g+4ny/klghXYEIWk1YFPg5nbGVG4njQemA2OBB4DnbM8pVaZSPbDVn34OfA14oxwv2+aYDFwp6TZJB5eydv49Wh14Gji93Ao8RdLibY6pZm/g7PJ+IMQTvZRk0geu/snU7zMXJC0BXAAcbntmO2Oy/bqr2xMrUy0ZsW5/fXZnJO0GTLd9Wzvj6GBb25tRrct0qKQP1J9sw9+jhYDNgF/b3hSYRYdbSO34u13GsXYHzut4rl3/r0XvJZk07ilJKwKUn9P788MlDaFKJGNsXzgQYgKw/RxwDdUtpGGSag/C9vfDWdsAu0t6GDiH6lbXye2Myfa08nM61VjAFrT3v9lUYKrtm8vx+VTJpd1/jz4M3G77qXLc7niiD5JMGncJsH95vz/VuEW/kCTgVGCK7Z+2OyZJy0kaVt4vSjV+M4UqqezV3/EA2D7a9sq2V6O6ZXK17X3bFZOkxSUtWXtPNSYwkTb+PbL9JPCYpHVK0Y5UCw62LaZiH966xcUAiCf6IM+ZdELS2cAoqtVLnwKOB/4MnAusCjwCfML2M/0Uz7bAP4EJvDUe8A2qcZN+j0nSRsCZVA9hLQCca/s7ktag6hUsA9wBfMr2q62Op5P4RgFH2t6tXTGVz72oHC4EnGX7BEnL0qa/RyWuTYBTgIWBB4EDKf8N2xFTSbSPAmvYfr6UtfXPKPomySQiIpqW21wREdG0JJOIiGhakklERDQtySQiIpqWZBIREU1LMokBTdJqkl4uS7f0VPeUsgBmV+cXLavTviap0W1ZI6IBmRocA1pZi+wy2xvMxWs+DIy0PWNuXTNifpeeSbSNpPdJuqvsj7J42Rul26RReip3SxpT9uM4X9Ji5dy1kkZKem/ZC2O4pAUk/VPSLv3zrSLmT0km0Ta2b6VaOuN7wEnAH21P7L4VAOsAv7K9HjAT+HyH6z5CtY/Jr4EjgMm2r5ybsUfE2yWZRLt9h2ptr5FUCaURj9n+v/L+j8C2HSvYPgVYCjgEOHIuxBkR3UgyiXZbFlgCWBIY2mCbjgN97xj4K7e+Vi6HS/Q5uohoSJJJtNtvgWOBMVS3phqxqqSty/tPAtd3UueH5ZrHAb9rNsiI6F6SSbSNpE8Ds22fBZwIvE/SDg00vYdqs6kpwNJUYyP11/0g8D7gh7bHAK9JOnDuRh8R9TI1OAa0jlOD58ZU4UwNjpj70jOJge514F2NPLTYk9pDi8AQ3toXJiLmgvRMIiKiaemZRERE05JMIiKiaUkmERHRtCSTiIhoWpJJREQ0LckkIiKa9v8B5WkZC/nFLEAAAAAASUVORK5CYII=\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"d.isel(t=0).w.T.plot.contourf()"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 313
},
"id": "Gvz90wj62-24",
"outputId": "91cbb234-ae65-4f85-b7fe-56950728b798"
},
"execution_count": 56,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"<matplotlib.contour.QuadContourSet at 0x7f1c7ee133d0>"
]
},
"metadata": {},
"execution_count": 56
},
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
],
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEWCAYAAACjYXoKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAAktUlEQVR4nO3de5xdVX338c+XcAm3ECCYxgRMLBFFXoIQAR9piwRCoEioRR/wQqQ8UitarFqFeuFRwQJaEZ8qLyMgAcGAEQQpGkOAWmy5JECFEGlS5JI0JIaEBLkEgt/nj70GDsPMmZM5OXPOMN/367Vfc/baa+3zyzCc31l7rb22bBMREdGMzdodQEREDH5JJhER0bQkk4iIaFqSSURENC3JJCIimpZkEhERTUsyiYiIpiWZRMeTNF6SJf1e0sk15X8jaUUp37nNMR5a4viDpEPbGUtEOySZRF2SHmrVh6Ok90l6WNJTkn4iaac+moy0PaO03QL4BjDF9na2H29FjI2yfaPt7YBH2hlHRLskmURbSHoz8F3gg8Bo4GngOxtxitHAcGBhL+ffvNkYI6JxSSbRK0mXAbsBPy2XcD6zCU//fuCntn9p+/fAF4B3S9q+gbjeADxQdp+QdFMpt6RTJC0GFpey8yU9KmmdpAWS/qTmPP9X0o8k/UDSk5LulfQGSadLWlnaTampv4OkiyQtl7RM0pmShm26X0nE4JVkEr2y/UGqyzbvKpeSzu1eR9Jukp6os72vl9O/GfjPmvf6b+A54A0NxPVfpT1Ul74OqTl8DHAAsGfZvxPYB9gJuAL4kaThNfXfBVwG7AjcDcyh+v9iLPBlqt5Tl0uADcDuwFuBKcD/6SveiKEgySSaYvsR2yPrbFf00nQ7YG23srVAnz2TPvyj7dW2nynx/cD247Y32P4nYCtgj5r6/2Z7ju0NwI+AXYCzbT8PzALGSxopaTRwJPAJ20/ZXgmcBxzXZLwRrwq5rhzt8ntgRLeyEcCTTZ730dodSZ8GTgJeC7i8x6iaKitqXj8DrLL9Qs0+VInvtcAWwHJJXfU36/5+EUNVkkn0pe4zCiTtBtxfp8pf2768h/KFwN4153k9Va/hv/oTZI0X4y3jI58BJgMLbf9B0hpAvTWu41FgPTCq9GIiokaSSfRlBfD63g7afoTqm/vGuhz4j/KBfxfV+MTVtpvtmdTanmqM43fA5pJO45W9oYbYXi7pF8A/SfoCVc9qAjDO9r9uqoAjBquMmURf/hH4fBlM//SmOqnthcBHqJLKSqoP/o9uqvMXc4CfU/V2HgaepbnLUicAW1L1xNYAs4ExTcYY8aqgPGkxOp2k11FNBX4W+Hvb32tzSK8gaTLwY6pLdUfavrnNIUUMqPRMouPZftj28DI7rOMSCYDteSW+rZNIoouki8s9S/fVlH1N0m8k/VrSNZJG1hw7XdISSQ9IOrymfGopW1Iu13aVT5B0eym/UtKWpXyrsr+kHB/f6n9rkklEROtcAkztVjYX2Mv2W6guwZ4OIGlPqqnmby5tviNpWLkx9tvAEVT3Tx1f6gKcA5xne3eqS68nlfKTgDWl/LxSr6WSTCIiWsT2L4HV3cp+UTMj8DZgXHk9DZhle73t3wJLgP3LtsT2g7afo7r/aZqqOeqHUI3dAcykumm361wzy+vZwGTVzGlvhUExm2tLbeXhbNvuMCJiEHiSNats79Lf9oe/c1s/vvqFvisCC369fiHVWF6XGV2LkTbor4Ary+uxVMmly9JSBi+fOLKUapWHnYEnahJTbf2xXW1sb5C0ttRftRGxbZRBkUyGsy0HaHK7w4gO9cLk/QAYNm9BmyOJTnCjZz/cTPvHV7/AHXN2a6jusDGLn7U9qT/vI+lzVFPXe7oPa9AZFMkkojcvTN6PFZO2AmA0+yWhxKAg6UPAUcBkvzSldhmwa021caWMXsofB0ZK2rz0Tmrrd51raVlBe4dSv2UyZhKDVm0ieXafp1kxaasXeykRnUrSVKqVGY62/XTNoeuA48pMrAnAROAOqsVKJ5aZW1tSDdJfV5LQzcCxpf104Nqac00vr48FbqpJWi2RnkkMSt0TSdfPFWyTHkp0DEk/BA4GRklaCpxBNXtrK2BuGRO/zfZHbC+UdBXVTbEbgFO61omT9DGqm3CHAReXm34BPgvMknQm1arXF5Xyi4DLJC2hmgDQ8gVJk0xi0OmeSI5700uJYxb7JaFEx7B9fA/FF/VQ1lX/LOCsHspvAG7oofxBqtle3cufBd6zUcE2KZe5YlDpKZF85TX38ZXXVPeEHfemBbnkFdEGSSYxKHXvkQAvJpSIGHhJJjGoDJu3gNHz1zP8nm2YtejlPY8vrNyLWYv2Y/g927QpuoihK2MmMegMm7eA0WVsZBYvJZTa5DJ6/vqMmUQMoCSTGJR6SyjD79kmiSSiDZJMYtCqTShdkkgi2iNjJjHojZ6//sWfSSQR7ZFkEoNaV/JIIolor1zmikEvSSSi/dIziYiIpiWZRERE05JMIiKiaUkmERHRtCSTiIhoWpJJREQ0LckkIiKalmQSERFNSzKJiIimJZlERETTkkwiIqJpSSYREdG0JJOIiGhakklERDQtySQiIpqWZBIREU1LMomIiKa19EmLkh4CngReADbYniRpJ+BKYDzwEPBe22taGUdERLTWQPRM3ml7H9uTyv5pwDzbE4F5ZT8i4lVH0sWSVkq6r6ZsJ0lzJS0uP3cs5ZL0LUlLJP1a0r41baaX+oslTa8p30/SvaXNtySp3nu0Ujsuc00DZpbXM4Fj2hBDRMRAuASY2q2sty/URwATy3YycAFUiQE4AzgA2B84oyY5XAB8uKbd1D7eo2VanUwM/ELSAkknl7LRtpeX148Bo3tqKOlkSfMlzX+e9S0OMyJi07P9S2B1t+LevlBPAy515TZgpKQxwOHAXNury5DAXGBqOTbC9m22DVza7VwD+qW9pWMmwEG2l0l6DTBX0m9qD9q2JPfU0PYMYAbACO3UY52IiEGoty/UY4FHa+otLWX1ypf2UF7vPVqmpcnE9rLyc6Wka6i6aCskjbG9vGTWla2MISJiYyzbsDVfWLlXg7UXj5I0v6ZgRvki3JB6X6g3lYF4D2jhZS5J20ravus1MAW4D7gO6BpAmg5c26oYIiJabJXtSTVbI4lkRfkiTbcv1MuAXWvqjStl9crH9VBe7z1appVjJqOBWyX9J3AH8C+2fw6cDRwmaTFwaNmPiBgqevtCfR1wQpnVdSCwtlyqmgNMkbRjGXifAswpx9ZJOrDM4jqh27kG9Et7yy5z2X4Q2LuH8seBya1634iITiHph8DBwChJS6lmZZ0NXCXpJOBh4L2l+g3AkcAS4GngRADbqyV9Bbiz1Puy7a5B/Y9SzRjbGvhZ2ajzHi3T6gH4iIghy/bxvRx6xRfqMiPrlF7OczFwcQ/l84FXDPC040t7llOJiIimJZlERETTkkwiIqJpSSYREdG0JJOIiGhakklERDQtU4MjIoaw2qXu63je9r31KiSZREQMbf9KdUOk6tSZQPVAw14lmUREDG132j6kXgVJN/V1koyZREQMYX0lkkbrJJlERASStuihbFSj7ZNMIiKGMEnvLItQLpf0C0njaw7/otHzJJlERAxt5wKH2x5F9XTbuWUJfKg/KP8yGYCPiBjatrS9EMD2bEmLgKslfRZo+AmNSSYREUPb85L+yPZjALYXSpoMXA/8caMnyWWuiIih7TSqJ+O+yPZS4M/YiCfhpmcSETGE2b6xl/K1wFmNnic9k4iIQNJhzbRPMomIGOIknQB8splz5DJXRMQQJulzwKHAnzdzniSTiIih7YvAm2w/3cxJcpkrImJoOxGYvTFLp/QkPZOIiCHM9hWSHgN+AhzU3/MkmUREDHG2b5K0splz5DJXRERg+75m2qdnEhERSJoE/AnwWuAZ4D5gru01jbRPzyQiYgiTdKKku4DTga2BB4CVVOMnN0qaKWm3vs7T8p6JpGHAfGCZ7aMkTQBmATsDC4AP2n6u1XFERESPtgHeYfuZng5K2geYCDxS7yQD0TM5FVhUs38OcJ7t3YE1wEkDEENERFtI+jtJCyXdJ+mHkoZLmiDpdklLJF0pactSd6uyv6QcH19zntNL+QOSDq8pn1rKlkg6bWPjs/3t3hJJOX6P7Xl9naelyUTSOKq7Ki8s+wIOAWaXKjOBY1oZQ0REu0gaC/wtMMn2XsAw4Dh6/1J9ErCmlJ9X6iFpz9LuzcBU4DuShpUrP98GjgD2BI4vdfsT67mSRkjaQtI8Sb+T9IFG27e6Z/JN4DPAH8r+zsATtjeU/aXA2J4aSjpZ0nxJ859nfYvDjIhomc2BrSVtTnVJaTm9f6meVvYpxyeXL+HTgFm219v+LbAE2L9sS2w/WIYLZpW6/THF9jrgKOAhYHfg7zfmH9kSko4CVtpeIOngjW1vewbVIyQZoZ0aftpXREQzVj+7DbMW7ddg7WtGSZpfUzCjfHYBYHuZpK9TjTc8Q/VM9QX0/qV6LPBoabtB0lqqL+Fjgdtq3qe2zaPdyg9oMPjuuvLBnwM/sr22ymMb17gV3gEcLelIYDgwAjgfGClp8/KLHAcsa2EMERGttMr2pN4OStqRqqcwAXgC+BHVZapOdL2k31Alvb+RtAvwbKONW3aZy/bptsfZHk91re8m2+8HbgaOLdWmA9e2KoaIiDY7FPit7d/Zfh64muqL9shy2Qte/qV6GbArQDm+A/B4bXm3Nr2VN0zSawFsnwb8L6rxneeBp9mIS2btuM/ks8AnJS2h6r5d1IYYIiIGwiPAgZK2KWMfk4H76f1L9XVln3L8Jtsu5ceV2V4TqKbq3gHcCUwss8O2pPrift1GxnihpNsknQ28BRCA7ae6ngvfiAG5A972LcAt5fWDVINGERGvarZvlzQbuAvYANxNNRb8L8AsSWeWsq4v1RcBl5Uv26upkgO2F0q6iioRbQBOsf0CgKSPAXOoZopdbHvhRsZ4pKThwMHAXwBfl/QI8HPg57br3l/SJcupRES0kO0zgDO6Fff4pdr2s8B7ejnPWfTwTHbbNwA3NBnjs5TkAVB6P0cA/yzpj2z32QFIMomIiBdJGgGspZpmPAv4fSPtkkwiIgJJfw18iWoGV9ftGLb9+kbaJ5lERATAp4G9bK/qT+OsGhwREQD/TTUduF/SM4mICKiWoP93SbfDS2tY2f7bRhonmUREBMB3gZuAe3lpPcWGJZlERATAFrY/2d/GGTOJiAiAn5XV2sdI2qlra7RxeiYREQFwfPl5ek2ZgUwNjoiIxtie0Ez7XOaKiBjCJB3Ux/ERkvbq6zzpmUREDG1/KelcqnW5FgC/o3oG1e7AO4HXAZ/q6yRJJhERQ5jtvysD7X9JtcjkGKoHZC0Cvmv71kbOk2QSETHE2V4NfK9s/ZIxk4iIIUzSUZuiTnomERFD29ckLaM8YbEXXwWur3eSXpNJgzer/MH2Ew3Ui4iIzrQC+EYfdRb3dZJ6PZP/KVu9bDUM2K2vN4mIiM5k++BNcZ56yWSR7bfWayzp7k0RREREDG71BuDf3kD7RupERMSrXK/JpDxgHkmHdj8maXptnYiIGNoamRr8RUkXSNpW0mhJPwXe1erAIiJi4EhaIOkUSTv2p30jyeTPqB7neA9wK3CF7WP782YREdGx/jfwWuBOSbMkHS6p3gSsl2kkmewI7E+VUNYDr9uYN4iIiM5ne4ntzwFvAK4ALgYelvSlRm4VaSSZ3Ab83PZU4G1UmetXTcQcEREdSNJbgH8Cvgb8mGqtrnVUj/Otq5E74A+1/QiA7WeAv5X0p/0PNyIiOo2kBcATwEXAabbXl0O3S3pHX+3r3QH/Rtu/AUZJGtXt8O/7GW9ERHSm99h+sKcDtt/dV+N6PZNPAidTdXlecW7gkHonljQc+CWwVXmf2bbPkDQBmAXsTLV2/gdtP9dXoBERg5GkkcCFwF5Un51/BTwAXAmMBx4C3mt7TRmPPh84Enga+JDtu8p5pgOfL6c90/bMUr4fcAmwNXADcKpt9yPUY3oYDl8LLLB9T1+N691ncnL5+c4etrqJpFgPHGJ7b2AfYKqkA4FzgPNs7w6sAU5q4FwREYPV+VTjzm8E9qZ6TshpwDzbE4F5ZR/gCGBi2U4GLoAX10o8AziAakLUGTVTeC8APlzTbmo/45wEfAQYW7a/Luf6nqTP9NW4zwF4ScMlfVLS1ZJ+LOkTpddRlytdl8O2KFtXj2Z2KZ8JHNPXuSIiBiNJOwB/SjUOge3nyuK406g+/+Dln4PTgEvL5+dtwEhJY4DDgbm2V9teA8yl+oI+Bhhh+7bSG7mU/n+mjgP2tf0p258C9gNeU+L/UF+NGxmAvxR4Evh/Zf99wGVUo/x1SRpGdSlrd+DbVNOLn7C9oVRZSpUBIyI6wmZPb8bwe7ZptPooSfNr9mfYnlGzP4HqMbjfl7Q31efhqcBo28tLnceA0eX1WODRmvZdn5H1ypf2UN4fr6G6otTl+RLnM5LW99LmRY0kk71s71mzf7Ok+xuJzPYLwD7lmuE1wBsbaQcg6WSqbh7Dafg/bETEQFple1Kd45sD+wIft327pPN56ZIWUF3FkdSfMY5N7XKqmVvXlv13AVdI2hbo8zO/kftM7ipjHQBIOgCYX6f+K5Ru3c1UC0OOlNSVxMYBy3ppM8P2JNuTtmCrjXm7iIhOsRRYavv2sj+bKrmsKJeoKD9XluPLgF1r2nd9RtYrH9dD+UYpA/+XUH2Bf6JsH7H9ZdtP2X5/X+doJJnsB/y7pIckPQT8B/A2SfdK+nWd4HYpPRIkbQ0cRjXwdDPQtRzLdODaHk8QETHI2X4MeFTSHqVoMtW3/OuoPv/g5Z+D1wEnqHIgsLZcDpsDTJG0Yxl4nwLMKcfWSTqwJIQT6MdnahlvucH2fNvnl22jOg2NXObq78yAMcDMMm6yGXCV7evLJbJZks4E7qYMTEVEvEp9HLhc0pbAg8CJlM9ESScBDwPvLXVvoJoWvIRqavCJALZXS/oKcGep92Xbq8vrj/LS1OCfla0/7pL0Ntt39l31lfpMJrYf7s+Jbf8aeMXDtcpNMfv355wREYNNuUejp3GVyT3UNXBKL+e5mGq9rO7l86nuYWnWAcAHyhWop6iesmvbb2mkcb074O+yvW+9xo3UiYiIQeHwZhrX65m8qd6YCFXW2qGZN4+IiM5g+2FJBwETbX9f0i7Ado22r5dMGpnG+0KjbxQREZ1L0hlUl+P2AL5PdaP5D4A+F3mEOsmkv2MlERExKP0F1Tj3XQC2/0fS9o02bmRqcEREvPo9VyYAGKDcrNiwJJOIiIBqqvJ3qW4s/zBwI/C9Rhv3OTVY0seBH5TFxSIi4lXI9tclHUb1ZMU9gC/antto+0ZuWhxN9YD5u6jmOM/p51r5ERHRwUryaDiB1OrzMpftz1OtkX8R1TLEiyV9VdIf9+cNIyKi80h6t6TFktZKWifpSUnrGm3f0JhJ6Yk8VrYNwI7AbEnn9ivqiIjoNOcCR9vewfYI29vbHtFo40bGTE6lWjxsFdWjJ//e9vOSNgMWA30+gSsiIjreCtuL+tu4kTGTnYB3d7/vxPYfJB3V3zeOiIiOMl/SlcBPqHlIlu2rG2ncyEKPZ9Q51u8sFhERHWUE1UrFU2rKDGyaZBIREa9+tk9spn1uWoyICCS9QdI8SfeV/bdI+nyj7ZNMIiICqrvdTweehxefSXVco42TTCIiAmAb23d0K9vQaOMkk4iIAFhVbkbvWujxWGB5o40zAB8REVA9LngG8EZJy4DfAu9vtHGSSUREYPtB4NCy9Pxmtp+sPS5puu2ZvbXPZa6IiHiR7ae6J5Li1HrtkkwiIqIRqncwySQiIhpR99EjSSYREdGI9EwiIqJpv6p3MMkkIiIoDz0cWbO/o6Qzu/Ztf6xe+ySTiIgAOML2E107ttcARzbaOMkkIqLFJA2TdLek68v+BEm3S1oi6UpJW5byrcr+knJ8fM05Ti/lD0g6vKZ8ailbIum0JsIcJmmrmvNuDWxVp/7LJJlERLTeqUDt85/OAc6zvTuwBjiplJ8ErCnl55V6SNqTatHFNwNTge+UBDUM+DZwBLAncHyp2x+XA/MknSTpJGAu0OtNit21LJlI2lXSzZLul7SwPP4XSTtJmlseXD9X0o6tiiEiot0kjQP+nOqx50gScAgwu1SZCRxTXk/jpQ/w2cDkUn8aMMv2etu/BZYA+5dtie0HbT8HzCp1N5rtc4CzgDeV7Su2z220fSuXU9kAfMr2XZK2BxZImgt8CJhn++zSJTsN+GwL44iIaNgWT5nR89f3XRG4H0ZJml9TNMP2jG7Vvgl8Bti+7O8MPGG7a0XepcDY8nos8CiA7Q2S1pb6Y4Hbas5Z2+bRbuUHNBR8D2z/DPhZf9q2LJnYXk5ZcdL2k5IWUf3jpwEHl2ozgVtIMomIwWmV7Um9HZR0FLDS9gJJBw9YVBtB0q22D5L0JC+/MVGAbY9o5DwDstBjGUR6K3A7MLokGoDHgNG9tDkZOBlgONsMQJQREZvcO4CjJR0JDKd6zvr5wEhJm5feyThgWam/DNgVWCppc2AH4PGa8i61bXorb4jtg8rP7fuqW0/LB+AlbQf8GPiE7XW1x2ybXm7Rtz3D9iTbk7ZofEJBRETHsH267XG2x1MNoN9k+/3AzcCxpdp04Nry+rqyTzl+U/mcvA44rsz2mgBMBO4A7gQmltlhW5b3uK4/sUq6rJGy3rS0ZyJpC6pEcrntq0vxCkljbC+XNAZY2coYIiI60GeBWeWmwLuBi0r5RcBlkpYAqymPzbW9UNJVwP1U49Gn2H4BQNLHgDnAMOBi2wv7GdOba3dKz2i/Rhu3LJmUGQgXAYtsf6PmUFfmPZuXZ+SIiFct27dQjRF3PTtk/x7qPAu8p5f2Z1HNtupefgNwQ3/jknQ68A/A1pLW8dIaXM9RPSyrIa28zPUO4IPAIZLuKduRVEnkMEmLgUPLfkREtIHtfyzjJV+zPcL29mXb2fbpjZ6nlbO5bqX3VSYnt+p9IyKiXz4n6QPABNtfkbQrMMb2HY00zh3wEREB1Z30bwfeV/Z/X8oakmfAR0QEwAG295V0N1QLPXatGdaI9EwiIgLg+bLWlwEk7QL8odHGSSYREQHwLeAa4DWSzgJuBb7aaONc5oqICGxfLmkB1QQpAcfYXtRHsxclmURERJfFwDpKbpC0m+1HGmmYZBIREUj6OHAGsAJ4gbLQI/CWRtonmUREBFQP8NrD9uP9aZwB+IiIgOq5KGv72zg9k4iIAHgQuEXSvwAvPh2s29qKvUoyiYgIgEfKtmXZNkqSSUREYPtLzbRPMomIGMIkfdP2JyT9lB4eVmj76EbOk2QSETG0dT1N8evNnCTJJCJiCLO9oLz8fc1rACQd1eh5MjU4IiIAvidpr64dSccDX2i0cXomEREBcCwwW9L7gD8BTgCmNNo4ySQiIrD9oKTjgJ9QTRGeYvuZRtsnmUREDGGS7uXls7h2AoYBt0vCdtbmioiIPjU8yF5PkklExBBm++FNcZ7M5oqIiKYlmURERNOSTCIiomlJJhER0bQkk4iIFpG0q6SbJd0vaaGkU0v5TpLmSlpcfu5YyiXpW5KWSPq1pH1rzjW91F8saXpN+X6S7i1tviVJA/8vTTKJiGilDcCnbO8JHAicImlP4DRgnu2JwLyyD3AEMLFsJwMXQJV8qJ7PfgCwP3BGVwIqdT5c027qAPy7XiHJJCKiRWwvt31Xef0ksAgYC0wDZpZqM4FjyutpwKWu3AaMlDQGOByYa3u17TXAXGBqOTbC9m22DVxac64B1bJkIuliSSsl3VdT1mPXLiLi1U7SeOCtwO3AaNvLy6HHgNHl9ViqZ7F3WVrK6pUv7aF8wLXypsVLgH+mypRdurp2Z0s6rex/toUxRERsnCefZti8BX3Xq4ySNL9mf4btGd0rSdoO+DHwCdvraoc1bFvSKx5KNdi0LJnY/mXJxLWmAQeX1zOBW0gyiYjBa5XtSfUqSNqCKpFcbvvqUrxC0hjby8ulqpWlfBmwa03zcaVsGS99dnaV31LKx/VQf8AN9JhJb127V5B0sqT5kuY/z/qBiS4iYhMqM6suAhbZ/kbNoeuArhlZ04Fra8pPKLO6DgTWls/MOcAUSTuW4YEpwJxybJ2kA8t7nVBzrgHVtrW5+urala7iDIAR2mnQdwEjYkh6B/BB4F5J95SyfwDOBq6SdBLwMPDecuwG4EhgCfA0cCKA7dWSvgLcWep92fbq8vqjVMMKWwM/K9uAG+hk0lvXLiLiVcf2rUBv931M7qG+gVN6OdfFwMU9lM8H9npli4E10Je5euvaRUTEINbKqcE/BP4D2EPS0tKdOxs4TNJi4NCyHxERg1wrZ3Md38uhV3TtIiJicMsd8BER0bQkk4iIaFqSSURENC3JJCIimpZkEhERTUsyiYiIpiWZRERE05JMIiKiaUkmERHRtCSTiIhoWpJJREQ0LckkIiKalmQSERFNSzKJiIimJZlERETTkkwiIqJpSSYREdG0JJOIiGhakklERDQtySQiIpqWZBIREU1LMomIiKYlmURERNOSTCIiomlJJhER0bQkk4iIFpI0VdIDkpZIOq3d8bRKkklERItIGgZ8GzgC2BM4XtKe7Y2qNZJMIiJaZ39gie0HbT8HzAKmtTmmlti83QE04knWrLrRsx9ucxijgFVtjqFWp8UDnRdTp8UDnRdTp8UDzcf0umbe/EnWzLnRs0c1WH24pPk1+zNsz6jZHws8WrO/FDigmfg61aBIJrZ3aXcMkubbntTuOLp0WjzQeTF1WjzQeTF1WjzQ/phsT23Xew9mucwVEdE6y4Bda/bHlbJXnSSTiIjWuROYKGmCpC2B44Dr2hxTSwyKy1wdYkbfVQZUp8UDnRdTp8UDnRdTp8UDnRlTv9jeIOljwBxgGHCx7YVtDqslZLvdMURExCCXy1wREdG0JJOIiGhakkkPJF0saaWk+2rKdpI0V9Li8nPHAYxnV0k3S7pf0kJJp7YzJknDJd0h6T9LPF8q5RMk3V6WjbiyDDgOKEnDJN0t6fp2xyTpIUn3Srqn616Edv4dlfcfKWm2pN9IWiTp7W38O9qj/G66tnWSPtHu31H0T5JJzy4Bus81Pw2YZ3siMK/sD5QNwKds7wkcCJxSlmRoV0zrgUNs7w3sA0yVdCBwDnCe7d2BNcBJAxRPrVOBRTX77Y7pnbb3qblvop1/RwDnAz+3/UZgb6rfVVtisv1A+d3sA+wHPA1c0654okm2s/WwAeOB+2r2HwDGlNdjgAfaGNu1wGGdEBOwDXAX1V29q4DNS/nbgTkDHMs4qg+fQ4DrAbUzJuAhYFS3srb9NwN2AH5LmXjTCTHVxDAF+FWnxJNt47f0TBo32vby8voxYHQ7gpA0HngrcHs7YyqXk+4BVgJzgf8GnrC9oVRZSrWUxED6JvAZ4A9lf+c2x2TgF5IWSDq5lLXz72gC8Dvg++VS4IWStm1zTF2OA35YXndCPLGRkkz6wdVXpgGfUy1pO+DHwCdsr2tnTLZfcHV5YhzVYnZvHKj37omko4CVthe0M45uDrK9L9WKsadI+tPag234O9oc2Be4wPZbgafodgmpHX/bZRzraOBH3Y+16/+12HhJJo1bIWkMQPm5ciDfXNIWVInkcttXd0JMALafAG6muoQ0UlLXjbADvWzEO4CjJT1EtTLrIVTjA22Lyfay8nMl1VjA/rT3v9lSYKnt28v+bKrk0u6/oyOAu2yvKPvtjif6IcmkcdcB08vr6VTjFgNCkoCLgEW2v9HumCTtImlkeb011fjNIqqkcuxAxwNg+3Tb42yPp7pkcpPt97crJknbStq+6zXVmMB9tPHvyPZjwKOS9ihFk4H72xlTcTwvXeKiA+KJfsgd8D2Q9EPgYKqlsFcAZwA/Aa4CdgMeBt5re/UAxXMQ8G/Avbw0HvAPVOMmAx6TpLcAM6mWh9gMuMr2lyW9nqpXsBNwN/AB2+tbHU8P8R0MfNr2Ue2KqbzvNWV3c+AK22dJ2pk2/R2VuPYBLgS2BB4ETqT8N2xHTCXRPgK83vbaUtbW31H0T5JJREQ0LZe5IiKiaUkmERHRtCSTiIhoWpJJREQ0LckkIiKalmQSHU3SeEnPlKVb+qp7YVkAs7fjW5fVaZ+TNGqTBhoxxGVqcHS0shbZ9bb32oTnfAiYZHvVpjpnxFCXnkm0jaS3Sfp1eT7KtuXZKHWTRump/EbS5eV5HLMlbVOO3SJpkqTXlWdhjJK0maR/kzRlYP5VEUNTkkm0je07qZbOOBM4F/iB7fvqtwJgD+A7tt8ErAM+2u28D1M9x+QC4FPA/bZ/sSljj4iXSzKJdvsy1dpek6gSSiMetf2r8voHwEHdK9i+EBgBfAT49CaIMyLqSDKJdtsZ2A7YHhjeYJvuA32vGPgrl77Gld3t+h1dRDQkySTa7bvAF4DLqS5NNWI3SW8vr98H3NpDnXPKOb8IfK/ZICOiviSTaBtJJwDP274COBt4m6RDGmj6ANXDphYBO1KNjdSe98+AtwHn2L4ceE7SiZs2+oiolanB0dG6Tw3eFFOFMzU4YtNLzyQ63QvADo3ctNiXrpsWgS146bkwEbEJpGcSERFNS88kIiKalmQSERFNSzKJiIimJZlERETTkkwiIqJp/x8G2y6bTfj4SAAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [],
"metadata": {
"id": "M8nbP6d-3ub5"
},
"execution_count": 56,
"outputs": []
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment