{ "cells": [ { "cell_type": "code", "execution_count": null, "id": "1c9cdeda-2cea-40ef-9f80-1a265993e124", "metadata": {}, "outputs": [], "source": [ "# processes historical precipitation data from CHIRPS" ] }, { "cell_type": "code", "execution_count": 10, "id": "cc2aacca-23d1-4774-9fc6-c4d10aec81bf", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import pandas as pd\n", "import geopandas as gpd\n", "import numpy as np\n", "import datetime as dt\n", "import matplotlib.pyplot as plt\n", "import seaborn as sns\n", "from pathlib import Path\n", "import glob\n", "import os\n", "import math" ] }, { "cell_type": "code", "execution_count": 2, "id": "694afb9c-96b5-4c94-83f1-7b66627bbfbd", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\Owner\\AppData\\Local\\Temp\\ipykernel_16880\\3510852079.py:2: UserWarning: Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n", "\n", " centroid = gpd.read_file('AOI/luanda_catchment_level4.shp').to_crs(epsg = 4326).centroid\n" ] } ], "source": [ "# get the centroid of the catchment in order to determine the appropriate data source\n", "centroid = gpd.read_file('AOI/luanda_catchment_level4.shp').to_crs(epsg = 4326).centroid" ] }, { "cell_type": "code", "execution_count": 21, "id": "1d618768-dbbf-4540-b68f-18c9cdab3f7a", "metadata": {}, "outputs": [], "source": [ "# select data source\n", "data_folder = Path('D:/World Bank/CRP/data/CRU/tmp') # change file path as needed" ] }, { "cell_type": "code", "execution_count": 4, "id": "0cd8103d-9c89-40e0-929f-6e10a273f52d", "metadata": {}, "outputs": [], "source": [ "output_folder = Path('output')" ] }, { "cell_type": "code", "execution_count": 17, "id": "1eff7aff-fd6f-42fc-90f4-dbd00339999d", "metadata": {}, "outputs": [], "source": [ "first_year = 1993\n", "last_year = 2022" ] }, { "cell_type": "code", "execution_count": 23, "id": "a553478f-d024-4244-abc6-e904ecf5283f", "metadata": {}, "outputs": [], "source": [ "for i in range(math.floor(first_year / 10) * 10 + 1, math.ceil(last_year / 10) * 10, 10):\n", " if i == 2021:\n", " data_file = 'cru_ts4.07.2021.2022.tmp.dat.nc'\n", " else:\n", " data_file = f'cru_ts4.07.{i}.{i+9}.tmp.dat.nc'\n", " \n", " # open data source\n", " nc = xr.open_dataset(data_folder / data_file)\n", " \n", " # select data based on coordinate and time period (last 30 years)\n", " df = nc.sel(time = slice(f'{first_year}-01-01', f'{last_year}-12-01')).sel(lon = centroid.x[0], lat = centroid.y[0], method = 'nearest')\n", " \n", " # write data to csv\n", " df = df.to_dataframe().drop(['lat', 'lon'], axis = 1)\n", " df.to_csv(output_folder / 'temperature.csv')" ] }, { "cell_type": "code", "execution_count": 24, "id": "68d6e58f-1fe8-41f9-818c-48abba16d0d2", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzYAAAIOCAYAAABqEZg1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACU+ElEQVR4nOzdd3hUdfY/8Pe0TNqk997oJNTQQSxgV7B3EFFXkS36c9eyq2Jdddev7uqKqyxWLFgARREQpUiAEFrokF5IT2ZSp97fHzeZzJgAKZPcKe/X8/g4M7mZnNwJM3Pmc875yARBEEBEREREROTC5FIHQERERERE1F9MbIiIiIiIyOUxsSEiIiIiIpfHxIaIiIiIiFweExsiIiIiInJ5TGyIiIiIiMjlMbEhIiIiIiKXx8SGiIiIiIhcHhMbIiIiIiJyeUxsiKjX3n//fchkMshkMvzyyy9dvi4IAtLS0iCTyTB79uwBjWXnzp145pln0NDQ0OVrSUlJuOqqq857H7/88stZfxdHx+Tpjh49imeeeQaFhYVSh9Kt119/Hddddx2Sk5PP+/f7448/Yvr06fDx8UFgYCCuvvpqHDlypMtxBoMBTz31FJKTk+Hl5YXExEQ8/vjjaG1ttTuupKQE8+fPR0pKCvz8/BAYGIhx48bhzTffhMlk6lH8X3/9NW699VakpaXBx8cHSUlJuP3223Hq1Kluj9+8eTOmTp0KX19fhIWFYeHChaiqqrI7JicnB0uWLEF6ejo0Gg0iIyNxySWXYMuWLV3u77333sO8efOQlJQEHx8fpKWl4YEHHsCZM2d6FD8AGI1GDB8+HH//+9+7/G7JyckIDAzEQw89BKPR2OV7dTodXnjhBUycOBEBAQFQq9VISkrCokWLsG/fPutxK1asQGxsLJqbm3scFxE5PyY2RNRnGo0GK1as6HL71q1bkZeXB41GM+Ax7Ny5E8uWLXOqJMIZY3IWR48exbJly5w2sVm+fDmKiopw0UUXITw8/KzHrV27FpdffjkiIiLw1VdfYfny5Th16hRmzpyJvLw8u2NvvfVWvPrqq7jvvvvw/fffY/HixXjttddw88032x3X3NyMgIAA/O1vf8O6devw2WefYcaMGVi6dCl+97vf9Sj+l19+GS0tLXjyySexYcMGPP/889i/fz/Gjx/fJenaunUrLr/8ckRGRmLt2rV44403sHnzZlx88cXQ6/XW4z799FPs2bMHixYtwtq1a/Hee+9BrVbj4osvxocffmh3n08//TT8/f3x4osvYsOGDfjzn/+M7777DhMmTEBlZWWPfof//Oc/qK+vx9KlS623FRcXY9GiRXjiiSfw5ZdfIisrC2+88Ybd9+Xl5WHcuHH4+9//jgsvvBCffvopNm7ciGXLlqGyshITJkyAVqsFACxYsAB+fn545ZVXehQTEbkIgYiol1auXCkAEBYvXiz4+PgIWq3W7ut33HGHMHXqVGHUqFHCBRdcMKCxvPrqqwIAoaCgoMvXEhMThSuvvPK89/Hzzz8LAISff/55wGNyN83Nzb06fvXq1Q491x0MBoNgNBr7fT9ms9l6+Vx/v8OGDRMyMjIEi8Viva2wsFDw8vISbrvtNuttWVlZAgDhn//8p933v/jiiwIAYePGjeeN6aabbhKUSqXQ1tZ23mMrKyu73FZWViaoVCrhnnvusbs9MzNTGDlypN15+/XXXwUAwn/+859z3qfJZBIyMjKE1NTU8/787OxsAYDw3HPPnTd+o9EoxMbGCo899pjd7atXrxbmz59vvb527VrhqquusosnPT1dCAgIEHJzc7u97++//97u7/Uf//iHEBgY2Ou/YSJyXlyxIaI+u/XWWwGIn+h20Gq1+Oqrr7Bo0aJuv6eurg4PPvggYmNj4eXlhZSUFDz55JN2nxADgEwmw0MPPYSPPvoII0aMgK+vL8aMGYPvvvvOeswzzzyDRx99FACspUPdlZRt2LAB48ePh4+PD4YPH47//e9/5/y9PvroI8hkMmRlZXX52rPPPguVSoXy8vJuv7cnMX3++eeYOnUq/Pz84O/vj0svvRT79++3u5+FCxfC398fx48fx6WXXgo/Pz9ER0dby3N27dqFGTNmwM/PD0OHDsUHH3xg9/0d5YKbNm3C3XffjZCQEPj5+eHqq69Gfn5+l7g7PqkPCAiAr68vpk+fjp9++qnL7yaTybBv3z7ccMMNCA4ORmpqKgBg7969uOWWW6wlSElJSbj11ltRVFRkF9ONN94IALjwwgut5+b9998HIJYOLly4sEtss2fPtisJ6ygd/Oijj/DII48gNjYWarUap0+f7vHvcjZy+flfFmtra3HixAlcfvnlkMlk1tsTExMxevRorFmzBmazGQDw66+/AgCuuOIKu/voKJH86quvzvvzwsPDIZfLoVAozntsREREl9tiYmIQFxeHkpIS621lZWXIzs7GnXfeCaVSab192rRpGDp0KL755ptz3qdCocCECRPs7vNsx06YMAEKhaLLsd1Zt24dysrKcOedd9rdnpqaih07diA7OxvV1dX44IMPMGzYMOvX16xZg9zcXDz++OMYPXp0t/d9+eWXw9fX13r99ttvh06nw2effXbeuIjINTCxIaI+CwgIwA033GCXKHz66aeQy+VdymwAoK2tDRdeeCE+/PBDPPzww1i/fj3uuOMOvPLKK7juuuu6HL9+/Xq8+eabePbZZ/HVV18hJCQE8+fPt74xX7x4sbVc5euvv0ZWVhaysrIwfvx4630cPHgQjzzyCP70pz9h7dq1yMjIwD333INt27ad9fe6+eabERUVhbfeesvudpPJhHfeeQfz589HTExMt997vphefPFF3HrrrRg5ciS++OILfPTRR2hsbMTMmTNx9OhRu/syGo247rrrcOWVV1pLnx5//HE88cQTWLBgARYtWoRvvvkGw4YNw8KFC5GTk9MlnnvuuQdyuRyrVq3C66+/jj179mD27Nl2ZXIff/wx5s6di4CAAHzwwQf44osvEBISgksvvbTbhOC6665DWloaVq9ejeXLlwMACgsLMWzYMLz++uv48ccf8fLLL+PMmTPIzMxETU0NAODKK6/Eiy++CAB46623rOfmyiuvPOtjcS6PP/44iouLsXz5cnz77beIiIjo9e/SFwaDAQCgVqu7fE2tVqOlpcVajna2YzuuHzp0qMt9CIIAk8mE+vp6fP7553j//ffxyCOP2CUgvZGfn4+ioiKMGjXKetvhw4cBABkZGV2Oz8jIsH79bEwmE7Zv3253n2ezdetWmM3mHh27fv16REREYOTIkXa3jxs3DosXL8bkyZMRERGB8vJyPPnkk9avb9y4EQAwb9688/6MDlFRURg+fDjWr1/f4+8hIicn9ZIREbmejlK07OxsaxnX4cOHBUEQy1sWLlwoCELXUp7ly5cLAIQvvvjC7v5efvnlLmU5AITIyEhBp9NZb6uoqBDkcrnw0ksvWW87Xymat7e3UFRUZL2ttbVVCAkJEe6//37rbd2Voj399NOCl5eXXWnN559/LgAQtm7des7zc7aYiouLBaVSKSxdutTu9sbGRiEqKkq46aabrLctWLBAACB89dVX1tuMRqMQHh4uABD27dtnvb22tlZQKBTCww8/bL2t4zGyLd8RhM5So+eff14QBLGULCQkRLj66qvtjjObzcKYMWOESZMm2Z0TAMJTTz11zt9fEMTSoKamJsHPz0944403rLefqxQtMTFRWLBgQZfbL7jgAru/o47Ha9asWXbH9eZ36YmzlaKZzWYhJCREuPjii+1ur6+vFzQajQBA2LlzpyAIgrBmzRoBgPDRRx/ZHbtixQoBgDB06NAu9//SSy8JAAQAgkwmE5588slexW3LaDQKs2fPFgICAoTi4mLr7Z988okAQMjKyuryPffdd5/g5eV1zvt98sknBQDCmjVrznmcTqcTRowYIcTHxwuNjY3njXfEiBHCZZdddtavV1dXC6dPn7YrARQEQbjssssEAD0q17N1++23C5GRkb36HiJyXlyxIaJ+ueCCC5Camor//e9/yM3NRXZ29lnL0LZs2QI/Pz/ccMMNdrd3lB/99hP1Cy+80G4AQWRkJCIiIuzKm85n7NixSEhIsF739vbG0KFDz3sfDzzwAADg3Xfftd725ptvIj09HbNmzerxz7f1448/wmQy4a677oLJZLL+5+3tjQsuuKBLCZ1MJrMrYVIqlUhLS0N0dDTGjRtnvT0kJOSs5+X222+3uz5t2jQkJibi559/BiAOOqirq8OCBQvsYrJYLLjsssuQnZ3dZXLU9ddf3+XnNDU14S9/+QvS0tKgVCqhVCrh7++P5uZmHDt2rNfnqid+G0dffpe+kMvlWLJkCX766Sc899xzqKqqwunTp3HHHXegpaXFegwglj+lpaXhL3/5CzZt2oSGhgZs2LABTzzxBBQKRbelbwsXLkR2djZ+/PFH/PnPf8arr75q10gvtK/o2P7XHUEQcM8992D79u348MMPER8f3+UY21K6ntwOiJPPXnjhBTzyyCO49tprz3pcW1sbrrvuOhQVFWH16tXw9/c/67EdysvLuy1n6xAWFobU1NRzxtcbERERqKqq6vHUOSJybn1b1yYiaieTyXD33XfjX//6F9ra2jB06FDMnDmz22Nra2sRFRXV5U1JREQElEolamtr7W4PDQ3tch9qtbrLmNxz6et9REZG4uabb8Y777yDxx57DEeOHMH27dvxzjvv9Phn/1bHVKjMzMxuv/7bN7m+vr7w9va2u83LywshISFdvtfLywttbW1dbo+Kiur2to5z3RHTb5NNW3V1dfDz87Nej46O7nLMbbfdhp9++gl/+9vfkJmZiYCAAGti1pvHqzd+G0dffpe+euqpp9DU1ITnn38eTz31FACx1O7uu+/Ge++9h9jYWADi4/LDDz/gzjvvxNy5cwEAfn5+ePHFF/Hcc89Zj7MVFRVlfdzmzp2L4OBgPPbYY1i0aBHGjRuHDz74AHfffbfd9wiC0OX64sWL8fHHH+ODDz7okoB0/Lv47b85QDxH3f2NAcDKlStx//3347777sOrr7561vOj1+sxf/587NixA9999x0mT5581mNttba2dvmb74mODy8KCgowfPjwHn+ft7c3BEFAW1tbjxIvInJuTGyIqN8WLlyIp556CsuXL8cLL7xw1uNCQ0Oxe/duCIJgl9x0fGIaFhY2GOH22B/+8Ad89NFHWLt2LTZs2ICgoKAuKyC90fH7ffnll0hMTHRUmOdUUVHR7W1paWl2Mf373//GlClTur2PyMhIu+u/TUy1Wi2+++47PP3003jsscest+v1etTV1fU4Vm9v7y5DJACgpqam27+N38bRl9+lr5RKJV577TU8++yzKCgoQFhYGKKjo3HppZciOTkZcXFx1mPT0tKQlZWFsrIy1NXVITU1FVqtFn/4wx96tPo3adIkAMDJkycxbtw4XH311cjOzj7r8R1JzcqVK7FixQrccccdXY7paLDPzc3tMtggNze32wb8lStXYvHixViwYAGWL19+1lUTvV6PefPm4eeff8batWtx8cUXn/d37BAWFtarv5kOl156Kf773/9izZo1dn+D51NXVwe1Ws2khshNMLEhon6LjY3Fo48+iuPHj2PBggVnPe7iiy/GF198gTVr1mD+/PnW2zv2wujNG6AOHU3YA7EqMGHCBEybNg0vv/wyDh8+jPvuu69Hn/afLaZLL70USqUSeXl53ZZzDYRPPvnE7mft3LkTRUVFWLx4MQBg+vTpCAoKwtGjR/HQQw/16WfIZDIIgtClQf69996zTgfrcK7HKykpqUsz/cmTJ3HixIkeJb2O+F16y9/fH+np6QCAffv24aeffsI///nPbo+NjY21rtD89a9/hZ+fH+65557z/oyOssGOZDQ0NLTblUhATGruvfderFy5Eu+8806XlR3bWCZNmoSPP/4Y/+///T/rxLVdu3bhxIkT+OMf/2h3/Pvvv4/FixfjjjvuwHvvvXfOpGb+/PnYsmULvv76a1x66aXn/f1sDR8+vMs+QD1x7bXXIj09HS+99BKuuuqqbhOzH3/8ETNnzrSbjJafn99lUAERuS4mNkTkEL/dJbw7d911F9566y0sWLAAhYWFSE9Px44dO/Diiy/iiiuuwCWXXNLrn9vxpvKNN97AggULoFKpMGzYMIdtDvqHP/wBN998M2QyGR588MF+xZSUlIRnn30WTz75JPLz83HZZZchODgYlZWV2LNnD/z8/LBs2TKHxN1h7969WLx4MW688UaUlJTgySefRGxsrPV38ff3x7///W8sWLAAdXV1uOGGGxAREYHq6mocPHgQ1dXVePvtt8/5MwICAjBr1iy8+uqrCAsLQ1JSErZu3YoVK1YgKCjI7tiON5z//e9/odFo4O3tjeTkZISGhuLOO+/EHXfcgQcffBDXX389ioqK8Morr5xzo0xbjvhd9u7da908VKfTQRAEfPnllwDEEsKOlbZffvkF2dnZyMjIgCAI2LNnD15++WVcdtllXZKqV155BVFRUUhISEBlZaU1uf/oo4/sStGefvppVFZWYtasWYiNjbX247z77ru48cYbMWHChPOeg9///vdYsWIFFi1ahPT0dOzatcv6NbVabdeb9fLLL2POnDm48cYb8eCDD6KqqgqPPfYYRo8ebZcQrV69Gvfccw/Gjh2L+++/H3v27LH7mePGjbMmrDfccAN++OEHPPnkkwgNDbX7+QEBAedNImbPno1nn30WLS0tdgnI+SgUCnzzzTeYO3cupk6digceeAAXXngh/Pz8UFRUhC+//BLffvst6uvrrd9jsViwZ8+eHiWXROQiJBpaQEQuzHYq2rl0N1WqtrZW+N3vfidER0cLSqVSSExMFB5//PEu04wACEuWLOlyn91Nznr88ceFmJgYQS6X203cOtsGnWebstXdpC69Xi+o1epzTmrqztliEgRxUtaFF14oBAQECGq1WkhMTBRuuOEGYfPmzdZjFixYIPj5+XUb+6hRo7rc/tvfteMx2rhxo3DnnXcKQUFBgo+Pj3DFFVcIp06d6vL9W7duFa688kohJCREUKlUQmxsrHDllVcKq1evth7TMRWturq6y/eXlpYK119/vRAcHCxoNBrhsssuEw4fPtzt4/X6668LycnJgkKhEAAIK1euFARBECwWi/DKK68IKSkpgre3tzBx4kRhy5YtZ328bGPr7e9yNh3T6Lr7ryNOQRCny02ePNn6GI4ePVr4xz/+IRgMhi73uWzZMiE1NVVQq9VCUFCQcNlllwnbtm3rcty6deuESy65RIiMjBSUSqXg7+8vTJo0SfjXv/7V481HExMTzxp/YmJil+M3btwoTJkyRfD29hZCQkKEu+66q8smm+c6J/jN9L9zHdeTzXpPnz4tyGSyLpMTe6qhoUF47rnnhPHjxwv+/v6CSqUSEhIShDvuuEP49ddf7Y796aefBABCTk5On34WETkfmSD8puOQiIisvv32W1xzzTVYv359l14EZ/b+++/j7rvvRnZ2NiZOnCh1OEQ9dvXVV8NkMuGHH34Y0J9z5513Ij8/37qJKhG5PpaiERF14+jRoygqKsIjjzyCsWPH4vLLL5c6JCKP8NJLL2HcuHHIzs4+6wTB/srLy8Pnn3+OLVu2DMj9E5E0uI8NEVE3HnzwQVxzzTUIDg7Gp59+6rB9M4jo3EaPHo2VK1d2O9HPUYqLi/Hmm29ixowZA/YziGjwsRSNiIiIiIhcHldsiIiIiIjI5TGxISIiIiIil8fEhoiIiIiIXJ7TTUWzWCwoLy+HRqNhsy4RERERkQcTBAGNjY2IiYmBXH7uNRmnS2zKy8sRHx8vdRhEREREROQkSkpKEBcXd85jnC6x0Wg0AMTgAwIC+n1/RqMRGzduxNy5c6FSqfp9f9RzPPfS4vmXFs+/9PgYSIfnXlo8/9LjY+A4Op0O8fHx1hzhXJwusekoPwsICHBYYuPr64uAgAD+YQ0ynntp8fxLi+dfenwMpMNzLy2ef+nxMXC8nrSocHgAERERERG5PCY2RERERETk8pjYEBERERGRy2NiQ0RERERELo+JDRERERERuTwmNkRERERE5PKY2BARERERkctjYkNERERERC6PiQ0REREREbk8JjZEREREROTymNgQEREREZHLY2JDREREREQuj4kNERERERG5PCY2RERERETk8pjYEBERERGRy2NiQ0RERERELo+JjRMTBAEtBpPUYRAREREROT0mNk7KZLbglv/uQsYzG/H8d0dhMFmkDomIiIiIyGkxsXFSPxyuwO6COpgsAt7bUYCb3slCaX2L1GERERERETklJjZOSBAEvLMtz+62AyUNuOKN7dh0tFKiqIiIiIiInBcTGyeUlVeLw2U6AEBquB8SQnwBALo2E+79cC9L04iIiIiIfqNXic1LL72EzMxMaDQaREREYN68eThx4oTdMTKZrNv/Xn31VYcG7s6Wb8u3Xv7jJUPx3e9n4PLRUdbbWJpGRERERGSvV4nN1q1bsWTJEuzatQubNm2CyWTC3Llz0dzcbD3mzJkzdv/973//g0wmw/XXX+/w4N3R0XIdtp2sBgDEh/jg8tFRCPBW4T+3j8eya0bBSyE+ZCxNIyIiIiLqpOzNwRs2bLC7vnLlSkRERCAnJwezZs0CAERFRdkds3btWlx44YVISUnpZ6ie4d3tnas1985MgbI9kZHJZFgwLQnjE4KxZNU+FNe1WEvTFs9Ixp8vGw4vJSsLiYiIiMgz9eudsFarBQCEhIR0+/XKykqsX78e99xzT39+jMcorW/BuoPlAIBgXxVunBDf5Zj0uECWphER0YCobtTju0Pl+GJvCXs5icjl9GrFxpYgCHj44YcxY8YMjB49uttjPvjgA2g0Glx33XVnvR+9Xg+9Xm+9rtOJTfNGoxFGo7Gv4Vl13Icj7mugvbctD2aLAAC4Y3I8lDILjMauLyw+CuCNm9KRmRiElzacgNEs4EBJA67813a8PH80Lh4RMdihd8uVzr074vmXFs+/9PgYnF9Nkx57Cuqxu7AOuwvqkVfdWVr+4+EzeOvWsVDIZb2+X557afH8S4+PgeP05hzKBEEQ+vJDlixZgvXr12PHjh2Ii4vr9pjhw4djzpw5+Pe//33W+3nmmWewbNmyLrevWrUKvr6+fQnNJbWYgKdzFDBYZFDJBTwz3gx/1fm/r6QJWHlSgVp95wvP7GgLrk6wgJVpRERkS2cA8nQynNLJcFonQ2XruZOW2dEWzE/iyg0RSaelpQW33XYbtFotAgICznlsnxKbpUuXYs2aNdi2bRuSk5O7PWb79u2YNWsWDhw4gDFjxpz1vrpbsYmPj0dNTc15g+8Jo9GITZs2Yc6cOVCpepApSOTtrfl4bfNpAMDtk+LxzNUjevy9jW1GPP7NEfx4tMp625i4QLxxcwZig3wcHmtPucq5d1c8/9Li+ZceH4Nzr8j8llIuQ3psAEbFBOCz7FKY2isInrl6BG6f1LU0+lx47qXF8y89PgaOo9PpEBYW1qPEplelaIIgYOnSpfjmm2/wyy+/nDWpAYAVK1ZgwoQJ50xqAECtVkOtVne5XaVSOfQPwdH350htRjM+3FUCAJDLgPsuSO1VrCEqFZbfOREfZhXhhfXHYDBbcLBUi2v/swv/uHEM5oyMHKjQe8SZz70n4PmXFs+/9DzpMahu1GN3QS125ddiV34dTlc1nfVYpVyGjLhATEkJxZSUUExIDIafWnxbMDo2CI99nQsAeG79cSSF+WP2sN6XOXvSuXdGPP/S42PQf705f71KbJYsWYJVq1Zh7dq10Gg0qKioAAAEBgbCx6dzZUCn02H16tX45z//2Zu791jf7C9DTZO4anV5ejQSQ/16fR/dTU3Ttho5NY2IyI05KpH5rVsmJaCgthnvbM2H2SLgoVX78eUDUzE8qv+VFEREA6VXic3bb78NAJg9e7bd7StXrsTChQut1z/77DMIgoBbb7213wG6O7NFwLs2G3LeP6t/Y7E7pqb95ctD+OGwmHi+t6MAe4vq8eZt4xAX7Dl9S0RE7magEpnu/OXS4SiqacGGIxVo0ptwz/t78c2SaYjQeDviVyEicrhel6L1xH333Yf77ruvTwF5mk1HK5FfI9Y8T00JRUZcUL/vs2NDT9vSNHFq2g6nKE0jIqKeGcxE5rfkchn+7+axOPPfLBws1aKsoRX3frAXn903FT5eij7fLxHRQOn7Mx71myAIeGdbnvX6/Rc4bhPTjtK0cQlBWLJqH0rqWlmaRkTkIg6VNuCxr3Jx9IzurMc4OpHpjo+XAu8umIj5b+1EWUMrDpZq8afPD+A/t4+HvA9joImIBhITGwntLarH/uIGAMDwKA0uGBru8J+REReE75bOxGNfsTSNiMgV1DcbsPiDvahq1NvdPhiJTHciNN7438JMXP/2TjTpTdhwpAIvbziOx6/o+fROIqLBwMRGQu9s7VytuW9WCmSygfn0K9Dn7KVp/7xxDC5haRoRkVMQBAF/XXvYmtSkhPnhstFRg5rIdGdYlAZv3T4ei97Phtki4J1t+UgK88OtkxIkiYeIqDusRZLI6apGbD4m7jsTE+iNq8fEDOjP6yhN+/KBqYgPESfYaVuNWPzhXjz/3VEYTNyAjYhIausOlmP9oTMAgCBfFT69bwr+fNlwzBoaLllS0+GCoeFYds0o6/W/rjmM7aeqJYyIiMgeExuJ/NdmEtqiGclQKQbnoegoTbt8dJT1tvd2FOCmd7JQWt8yKDEQEVFXZ7St+Nuaw9brz88bjcgA55pAdseURCyeIe5hZ7YIePDjfThZ2ShxVEREIiY2EqjUteGb/WUAAI23ErcM8lJ+R2nasmtGwas9oeooTdt8tHJQYyEiIrEE7c9fHoKuzQQAuGZMDK7KGNiV/L56/IoR1umajXoT7l6Zjerf9AMREUmBiY0EVv5aCKNZHJ1955RE+EtQXsDSNCIi5/HRriJsP1UDAIgMUOO5a0dLHNHZKeQyvHHLWKTHBgKAOAb6w71oM5oljoyIPB0Tm0HW2GbEJ7uKAABeCjkWTkuSNJ5zlaYV1TZLGBkRkWfIr27Ci98fs15/9YYxCPRVSRjR+fl6KfHegomIDhRL5Q6UNOCRLw7CYunZfndERAOBic0g+3RPMRr1YqnBdeNjEeEE9dNnK0274NVfcOE/fsGjqw/ii+wS5Fc39XiTViIiOj+T2YI/fXEQbUZxlfyuqYmYNQCj/wdCZIA3VizIhF/7Zp3rc8/gHxtPSBwVEXkyjnseRAaTBf/bUQgAkMmAe2c5bkPO/upuQ08AKKhpRkFNM1bnlAIAQv28MDEpGJlJIZiYFIJRMQGDNviAiMjdvP1LHg6WNAAAksP88Njlw6UNqJdGxgTgzdvG454PsmERgP/8koekMD/MHxN1/m8mInIwJjaDaN3BclTo2gAAl4yIRGq4v8QRddVRmrZiRwF2nq7BoVItDObOfpvaZgN+PFKJH4+IQwa8VXKMjQ+yJjrjE4Kg8XbuEgoiImeQW6rFGz+dAgDIZcBrN42Br5frvSxfODwCT189Ck+vOwIAeOLrXERp+DpARIPP9Z5BXZQgCPjvts4NOX93gfOs1vxWoI8KD88ZiofnDEWb0YzcMi2yC+uwt7AeewvrrFN7AKDNaMGu/Drsyq8DIL44D48KQGZSMMbFB0LLQTlERF20Gc340xcHYGrvSVlyYRrGJQRLHFXfLZiWhIKaZry/sxAmi4CHPj2Ih1xr8YmI3AATm0Hyy4lqnKxsAgBMTAzGhMQQiSPqGW+VAplJIchMEuO1WAScqmpqT3TqsLeoHqX1rdbjLQJw9IwOR8/o8EEWACjx3/xt1hWdzKQQDInwh1wuk+YXIiJyAq/+eAKnq8TXhNGxAVh60RCJI+q/v101EiV1LfjpeBV0bSa8c0yBa5oNiAri6g0RDQ4mNoNk+dbO1Zr7nKi3prfkchmGRWkwLEqDO6YkAhA3letYzckurMexCh1sZwyUNbSh7EA51hwoBwAEeCsxMSnE2quTHhsIb5VCil+HiGjQ7cyrwYodBQAAL6Uc/3fTWHgpXb9XUSGX4V+3jsONy7Nw9IwOtXoZHlx1AKvuncLneCIaFExsBsH+4nrsLhBLtVLC/XDJiEiJI3Ks6EAfXD3GB1ePETeT07UZsb+4AbvzarBp/2mUtCqtE3/Er5uw5XgVthyvAiCOvU6PCxQTncQQTEgMRrCflyS/CxHRQNK1GfHo6kPW63++dBiGRGokjMix/NRKrFg4EfPe/BWVjXrsK27Ao18ewhs3j+VKPRENOCY2g+C/2/Ktl++fleL2T+4B3ipcMDQc05KDMMxwEnMuvQgnq1vbV3TEXp3aZoP1eIPZgpyieuQU1eMdiOdqSIQ/MpNDsGBqEoZFuc+LPhF5tme/PYqyBrF8d0pKCBZNT5Y4IseLDvTBO3eMw03vZMFgkeHbg+VIDvXFw3OHSR0aEbk5JjYDrLCmGRuOVAAAwjVqzBsXK3FEg0+lECenjY0PwuKZKRAEAQU1zdhbWC8mOkX1KKix3wz0VFUTTlU14cu9pXj8iuFYOC0JMpl7J4RE5N5+PFKBL9tH52vUSvzjxjFu+0HXqJgALBhiwYqTClgE4F9bTiMx1A/XT4iTOjQicmNMbAbYu9vzrf0md09PglrJOmOZTIaUcH+khPvjpsx4AEB1ox45Re19OkX1OFKmhckiwGC2YNm3R7H9VA1evSEDof5qiaMnIuq9miY9nvg613r96WtGIS7YV8KIBt7oEAGPXz4ML3wvbtr52NeHEBvsgykpoRJHRkTuyvW7FZ1YTZPeurGln5cCt09OlDgi5xWuUeOy0VH461UjsXbJdBx8ei4Wz+gs0dhyvAqXv7Edv56ukTBKIqLeEwQBj32Vay3BnTsyEteP94zV+wVTEnBn+6AZo1nA/R/lIL+6SeKoiMhdMbEZQB/uLITBJDbN3zopAYE+HHnZU35qJf561Ui8f3cmwvzFQQJVjXrcsWI3/v7DcRhtNg0lInJmq3NKsfmYuKlxmL8XXrou3WNKa2UyGZ6+eiRmDwsHAGhbjVj0fjbqbPosiYgchYnNAGnWm/BBVhEAQCmXYdEM92sQHQyzh0Xg+z/MxMwhYQAAQRBHZ9+wPAtFtc3n+W4iImmV1LXg2W+PWq+/dJ3nldQqFXL8+9ZxGN4+CKawtgX3f7QXepNZ4siIyN0wsRkgX+wtgbbVCAC4ZmwMYoJ8JI7IdUVovPHB3ZPwxBXDoVKIn3IeLGnAlf/agTX7yySOjoioexaLgEdWH0ST3gQAuGliHOaMdK9x/z2l8VZhxcJMhGvEpC67sB6PfZULwXbTMyKifmJiMwBMZgve215gve7KG3I6C7lchvtmpeKrB6YhKVRsuG3Sm/DHzw/g4S8OWN84EBE5ixU7CrCnfQ+zuGAf/O2qkRJHJK3YIB+sWDAR3irxrcc3+8vwxk+nJI6KiNwJE5sBsD73jHWfgtnDwjE8KkDiiNxHRlwQvvv9TFxn03j79b4yXPWv7ThU2iBdYERENk5UNOLVH8VpYDIZ8I8bx0DjzT7LjLggvHHLOHS0GL2++RRX3onIYZjYOJggCHhnq+2GnKkSRuOe/NVKvHbTWLx+81j4q8WJ5YW1Lbj+7Z3477Y8WCwsbSAi6RhMFvzp8wMwtA85WTwjmSOObVw6KgpPXD7Cev3PXx5CdmGdhBERkbtgYuNgv56uxdEzOgDAmLhATEkJkTgi9zVvXCzW/34GxsQHARBHib74/XEsWLkHVY1t0gZHRB7rjZ9OWl8Hhkb645G5wySOyPksnpmM2yYnAAAMZgvu+3AvCms4EIaI+oeJjYO9sy3Pevm+WakeM9JTKomhfvjyd1PxwOxUa2nD9lM1uOKN7fj5RJW0wRGRx8kpqsfbv4ivAyqFDK/dNBbeKm7M/FsymQzLrhllnXhZ3yKOgW5o4RhoIuo7JjYOdLhMi+2nxA0kE0J8cdnoKIkj8gwqhRx/uWw4Plo02Tpxp6bJgLtXZuP5745ypCgRDYoWgwmPfHEAHdWwf7xkKEbHBkoblBNTKeR46/bxGBrpDwDIr2nG7z7Ose7/RkTUW0xsHOjd7Z29NffOSoFCztWawTRjSBg2/GEmLhoeYb3tvR0FuO4/O7nTNRENuBe/P4bC2hYAwLiEINzPiZjnFeCtwooFnRsx78qvw+Nfcww0EfUNExsHKalrwXeHzgAAQvy8cOOEOIkj8kyh/mqsWDART189El4K8c/7SLkOV/17B1bvLeGLJRENiF9OVOHjXcUAAB+VAq/dNBZKBV9ieyI+xBfv3jURaqV4vr7aV4q3fj4tcVRE5Ir4rOsgK3YUwNxef7BgahJrqiUkk8lw9/RkfLNkGlLD/QAALQYzHv3yEH7/2QHo2owSR0hE7qShxYA/f3nIev2JK0cgOcxPwohcz7iEYPzfzWOt1/+x8SS+PVguXUBE5JKY2DhAfbMBn2eXABA/qbtraqLEEREAjIoJxLdLZ+CWzHjrbd8eLMcVb2zHvuJ6CSMjInfy1zWHUdWoBwDMGhqOO9qnfVHvXJEejb9cNtx6/ZHVB1HASWlE1AtMbBzg411FaDWKDeo3TYxDsJ+XxBFRB18vJf5+fQbeum08NN7injel9a24cXkW3vr5tHWVjYioL9YdLLeWIQf6qPDqDRmchtkPv7sgBTdPFD+MMpgs+L9NJyWOiIhcCRObfmozmvH+zkIAgFwGLJ7JZlFndGVGNH74w0xMSAwGAJgtAl798QTueG83KrTc84aIeq9C24a/fpNrvf78vNGIDPCWMCLXJ5PJ8LerRyKk/QPCbw+V43iFTuKoiMhVMLHpp6/2laK2WZy7f2VGDOJDfCWOiM4mLtgXn983Bb+/eAg6BtZl5dfi8je2YfPRSmmDIyKXIggCHv3yIHRtJgDANWNicPWYGImjcg/+aiUeuCAVACAIwGsbuWpDRD3DxKYfzBYB727rHPHM0Z7OT6mQ4+E5Q7Hq3imIDhQ/Wa1vMWLxh3vx9NrDaDNyzxsiOr+PdxVZ9y2LDFDj2WtHSRyRe7lzaiIi2vcl23i0EodKG6QNiIhcAhObfth4pMK6Z8H0tFBuxOZCpqSE4oc/zMSloyKtt32QVYR5b/2KU5WNEkZGRM6uoKYZL3x/zHr91RvGIMiXvZWO5K1SYOlFadbr/+SqDRH1ABObPhIEAcvtVmtSJYyG+iLI1wvL75iA5+eNtu6fcLyiEVe/uQOrdhdzzxsi6sJktuBPnx9Am9ECALhzSiJmDQ2XOCr3dHNmAmKDfAAAW09WI7uwTuKIiMjZMbHpoz0FdThY0gAAGBEdgJlDwqQNiPpEJpPhjimJWPfQDAyL1AAA2owWPPFNLpas2gej2SJxhETkTJZvzcOB9uf+5DA/PH7F8HN/A/WZl1KOP1wyxHr91R9P8AMnIjonJjZ99M5vems43tO1DYvSYO1D03HnlM49iL7PrcAXe0skjIqInMnhMi1e33wKgDgF8583jYGvl1LiqNzbdeNikdK+2emegjrsOF0jcURE5MyY2PTBycpGbDleBQCIDfLBlRnREkdEjuCtUuC5eaPxfzePsd62O5+lD0Qkjvb/0+cHYGrf+2rJhWkYnxAscVTuT6mQ449zhlqv/2PjSa7aENFZMbHpg//arNYsmpEMlYKn0Z1clREDr/aem9wyrcTREJEz+MePJ3CqqgkAMComAEsvGnKe7yBHuSo9GsOjxFLhgyUN+OlYlcQREZGz4jvyXqrQtmHtgTIA4i7Tt2TGSxwROZpKIcfI6AAA4vQjbatR4oiISEpZebVY8WsBALHv4/9uHmv98IMGnlwuw8N2qzYnYLFw1YaIuuIzcy+t/LUARrP4hHrnlET4qVlf7Y4y4jpHdx/hqg2Rx2psM+L/rT6IjuqnP186DEPbB43Q4JkzMtL6vHy8ohHfHz4jcURE5IyY2PSCrs2IT3YXAxA/tVswLUnagGjApNvsSXSIiQ2Rx3r226Moa2gFAExODsGi6ckSR+SZZDIZHpk7zHr9tU0nYeLUSiL6DSY2vbBqdzGa9CYAwPXj4xDevisyuZ+MuCDr5dxSJjZEnmjjkQqszikFAPirlfjnTWMgl3MCplRmDQnDpKQQAEB+dTPWHCiXOCIicjZMbHpIbzLjfzvEGmuZDLh3Jj+1c2ep4X7wUSkAAAdLG6QNhogGVavBjA92FuLPXx2y3vb01SMRF+wrYVQkrtp09tq88dNJGExctSGiTkxsemjtgXJUNeoBAHNHRiIl3F/iiGggKRVyjIoRBwiU1reirtkgcURENNAaWgz410+nMP3lLXh63RE0tIiDQ+aMjMQNE+Ikjo4AYHJKqHVD7JK6Vu41RkR2mNj0gMUi2I14vv+CVAmjocGSbjNAgGOfidzXGW0rnv/uKKb9fQte23TS7oOMS0ZE4tUbMrgJsxOx7bX595ZTaDOaJYyGiJwJR3r1wJbjVTjdvn/BpKQQbsrmIWwno+WWNuCCoeESRkNEjna6qgnvbM3DmgNl1mmXAKCQy3B1RjTuvyAVI9pHv5PzGBsfhEtGRGLzsUpU6vT4ZHcx7pnB8nAiYmLTI7arNffNSpEwEhpMtgMEDnGAAJHb2F9cj+Vb87DxaCVsN7FXK+W4OTMe985MQXwI+2mc2SNzh2LzsUoAwH9+Po1bMuO5/QIRMbE5n5yieuwprAMApEX446LhERJHRIMlOdQP/molmvQmlqIRuThBELDtVA3e/uU0duXX2X0twFuJu6YmYeH0JIT5c9qlKxgRHYCrMqLx3aEzqG024P2dhVhyYZrUYRGRxJjYnMd/t+VZL983K4WjPj2IXC7D6NgA7MqvwxltG6oa2xCh8ZY6LCLqBbMArM+twH+3F+LoGZ3d1yID1Fg8IwW3Tk6APz/tdzl/vGQovs89A4sAvLM1D3dMSUSgj0rqsIhIQhwecA751U3YeFRc6o7QqHHt2BiJI6LBZluOdpirNkQuo81oxqo9JXhhvwJ//OKQXVKTEuaHl69Px7Y/X4h7Z6UwqXFRaRH+uG68OK1O12bCiu355/kOInJ3fDY/h3e3F1jrrxfNSIZaqZA2IBp06bGdAwQOlmhx0fBICaMhovPRtRnx8a4i/G9HIWqa9AA6V9nHxAXigdmpmDMyCgquvruFP1w8BGv2l8FkEbBiRwEWTk9GiJ+X1GERkUSY2JyF3mS2Nib6q5W4bXKCxBGRFDI48pnIJVTp2rDi1wJ8sqsYTXqT3ddmpIXiwdlpmJoayrHNbiY+xBc3Z8bjk93FaDaYsXxrHp64YoTUYRGRRJjYnIVaqcDP/282PttTDLNFQIA363Y9UUKILwK8ldC1mXCoVAtBEPjGiMiJFNY0451t+fgqpxQGc+cu9HIZcNmoSIyUleH+myZApeJzuLtaetEQrM4phcFkwQc7C7F4RjIiAtgPSeSJmNicg79aicUzOd7Zk8lkMmTEBWHH6RrUNOlRoWtDdKCP1GERebzcUi2Wb83DD4fF5vEOXgo5rp8Qh/tmpSAu0Avff18mXZA0KKICvXHnlESs2FEAvcmCN38+jWevHS11WEQkASY2ROeREReIHadrAIj72TCxIZKGIAjYmVeL5VvzsP1Ujd3X/NVK3DElEYumJ1k/rTcajVKESRJ4YHYqPt1TjBaDGZ/uKRYT22DuRUTkaZjYEJ2HXZ9NqRaXjoqSMBoiz2O2CNh4pAJvb83rsllumL8ai2Yk4fbJHPXrycL81bh7ehLe+jkPRrOAf/90Gi/fkCF1WEQ0yJjYEJ1Hus3I50McIEA0aPQmM9bsL8M7W/ORX9Ns97WEEF/cf0EKrh8fB28VJ1YScN/MVHyYVYTGNhO+3FeK381ORXKYn9RhEdEgYmJDdB4xgd4I9fNCbbMBuaUNHCBANIBqm/TILqzD7oI6fJ97BpU6vd3XR8UE4HcXpOLy0VFQKrgVG3UK9FXhvpkp+OemkzBbBLy++STeuGWc1GER0SBiYkN0HjKZDOlxgfjlRDXqW4worW9FfAhrt4kcobyh1ZrI7Cmow+mqpm6Pm5oSigdmp2LmkDB+sEBndfeMZKzcWYi6ZgPWHSzHg7PTMCxKI3VYRDRImNgQ9UBGrJjYAOIAASY2RL0nCAIKa1uwp6AWuwvqkF1Yh5K61rMeL5MBc0ZE4oHZqRiXEDyIkZKr8lcr8cAFqXjh+2MQBOC1TSfwzp0TpQ6LiAYJExuiHrDvs2nAlRnR0gVD5CIsFgEnqxqxp6BzRaa6UX/W4xVyGdJjAzE5OQSTkkMwMTEEgb4cCEC9c8eURLy7PR9VjXr8eKQSuaVapNsMgSEi98XEhqgHfjsZjYi6MpotOFKuw56CWuwpqEN2YT20rWcfuaxWyjEuIQiTkkIwKTkU4xKC4KfmyxL1j4+XAg9dlIan1h4BAPxj4wl8sGiSxFER0WDo1SvISy+9hK+//hrHjx+Hj48Ppk2bhpdffhnDhg2zO+7YsWP4y1/+gq1bt8JisWDUqFH44osvkJCQ4NDgiQZLZIA3IgPUqNTpkVumhcUiQC5nnT95tjajGQdLGrCnoA57CuuQU1SPFoP5rMf7q5WYkBiMSckhmJwcgvS4QKiVnGhGjndzZjze2ZqPsoZWbD1ZjezCOmQmhUgdFhENsF4lNlu3bsWSJUuQmZkJk8mEJ598EnPnzsXRo0fh5yeOVMzLy8OMGTNwzz33YNmyZQgMDMSxY8fg7e09IL8A0WBJjw1Cpa4SjW0mFNW1cIwoeZwmvQk5RfXYU1CL7IJ6HChpgMFsOevxwb4qTEoWV2MmJYVgRLSGk8xoUKiVCvzh4iH481eHAAD/+PEEPrtvCgdPELm5XiU2GzZssLu+cuVKREREICcnB7NmzQIAPPnkk7jiiivwyiuvWI9LSUlxQKhE0sqIC8TmY5UAgEOlDUxsyO1pW4zY3V5WtqewDkfKdTBbhLMeHxXg3Z7IiCsyqeH+XNkkyVw3PhZvb81DQU0zdhfU4dfTtZgxJEzqsIhoAPWrmFmrFXsNQkLE5V2LxYL169fjz3/+My699FLs378fycnJePzxxzFv3rxu70Ov10Ov72wm1el0AACj0Qij8ey12T3VcR+OuC/qHXc79yOjOhOZg8X1uGJUhITRnJ+7nX9X4+rnf3dBHRZ/tA9txrOvyCSG+CIzKRiZSUGYmBiM+GAfu0/EzWYTzGevTBtwrv4YuDJnOfdLL0zBw6tzAQCv/ngckxInecSqjbOcf0/Gx8BxenMOZYIgnP3jt3MQBAHXXnst6uvrsX37dgBARUUFoqOj4evri+effx4XXnghNmzYgCeeeAI///wzLrjggi7388wzz2DZsmVdbl+1ahV8fTlSl5xHoxH4617xs4BUjYDfj5bwHRvRAPvfCTkO1tmXjUX7CEgN6Pwv0Eui4Ih6yCIArxxU4EyrmMzcO8yM0SF9ettDRBJpaWnBbbfdBq1Wi4CAgHMe2+fEZsmSJVi/fj127NiBuLg4AEB5eTliY2Nx6623YtWqVdZjr7nmGvj5+eHTTz/tcj/drdjEx8ejpqbmvMH3hNFoxKZNmzBnzhyoVBwbOpjc8dxf8I9tKNe2wc9LgZwnL4LCicts3PH8uxJXPv96oxmT/v4LWgxmBPmo8Pf5ozA+MQjBvq6VybjyY+DqnOncbzpahQc/PQAAGB6lwdoHprh9iaQznX9PxcfAcXQ6HcLCwnqU2PSpFG3p0qVYt24dtm3bZk1qACAsLAxKpRIjR460O37EiBHYsWNHt/elVquhVqu73K5SqRz6h+Do+6Oec6dznxEXhHJtBZoNZpQ06DEk0vl3tHan8++KXPH8/5rfOd3s4hGRuCwjVuKI+scVHwN34Qzn/vKMGKRvK0BumRbHKxqx6UQNrsqIkTSmweIM59/T8THov96cv16NpxEEAQ899BC+/vprbNmyBcnJyXZf9/LyQmZmJk6cOGF3+8mTJ5GYmNibH0XklGw3eTvE/WzITXUMyQCAOSOdu5eM6HxkMhkemTvUev21TSdhOsc0PyJyXb1KbJYsWYKPP/4Yq1atgkajQUVFBSoqKtDa2mo95tFHH8Xnn3+Od999F6dPn8abb76Jb7/9Fg8++KDDgycabGPigqyXc8uY2JD7EQQBm49WAQC8lHLMHBIucURE/XfB0HBkJgUDAPKrm7H2QLnEERHRQOhVYvP2229Dq9Vi9uzZiI6Otv73+eefW4+ZP38+li9fjldeeQXp6el477338NVXX2HGjBkOD55osKXH2q7YNEgXCNEAOVymQ4WuDQAwPTUUfup+Dc8kcgriqk3nZuKv/3QSBhNXbYjcTa9esXo6Z2DRokVYtGhRnwIicmaBviokhvqiqLYFR8p1MJkt3HCQ3MommzK0S0ZGShgJkWNNSQnFzCFh2H6qBiV1rVidU4LbJ7NMnsid8B0ZUS91rNroTRacqmqSOBoix9p8tDOxuXg4ExtyLw/P6ey1+fdPp9Fm5Nh+InfCxIaolzJsBgjkcoAAuZGyhlYcPSNukpwRF4ioQG+JIyJyrHEJwbhkhDgQo0LXhk92F0scERE5EhMbol5Kjw2yXj5U1iBZHESO9pNtGdoIrtaQe3p4Tmevzdu/nEaz3iRhNETkSExsiHppdGzn5lAc+UzuZNNR2zHPTGzIPY2MCcCVGdEAgJomA97fWShtQETkMExsiHpJ461CSrgfAODYGR30JtZok+vTtRmxK78WABAb5IPhUc6/+SxRX/3pkqGQy8TL72zNg7bVKG1AROQQTGyI+iCjfYCA0SzgZAUHCJDr23ayGkazOPlyzshIyGQyiSMiGjhpEf6YPy4OAKBrM2HF9nyJIyIiR2BiQ9QHGTYbdbLPhtyB7TQ09teQJ/jDxUOgbF+2WbGjAHXNBokjIqL+YmJD1AecjEbuxGi2YMvxKgCARq3EpOQQiSMiGngJob64OTMeANBsMOOdrXkSR0RE/cXEhqgPRsYEWOuzOUCAXN3ewnro2sTJULOHR8BLyZcG8gwPXZRm/Xv/IKsQVbo2iSMiov7gqxdRH/h6KTEkQmyuPlnZyE3eyKVtthvzHCFhJESDKzrQB3dMTgQAtBkteOvn0xJHRET9wcSGqI/S28vRTBYBx9o3NSRyNYIgWMc8K+UyzB7KxIY8y4MXpsJHpQAArNpTjNL6FokjIqK+YmJD1Ed2fTZlLEcj13SqqgnFdeIbuUnJIQj0VUkcEdHgCvNX4+7pSQDESZf//omrNkSuiokNUR+lx3YmNgdLmNiQa9rEaWhEuH9WKjTeSgDAl/tKUVDTLHFERNQXTGyI+mhEdIB1VGguRz6Ti7Ltr5kzkokNeaZAXxXunZkCADBbBLy++aTEERFRXzCxIeojb5UCw9p3Zz9d1YRmvUniiIh6p6qxDQdKGgAAw6M0iA/xlTYgIgndPT0Jwe2lmOsOluNERaPEERFRbzGxIeqHjj4biwAc5QABcjFbjlVBEMTLLEMjT6fxVuGB2akAAEEA/vXTKYkjIqLeYmJD1A/psUHWy9zPhlyN3ZhnlqER4c4pSQjx8wIA/HKiCiazReKIiKg3mNgQ9YPdZLTSBukCIeqlVoMZ20/VAADCNWpk2AzDIPJUPl4KTE0JBQA0G8w4UclyNCJXwsSGqB+GRmrgpRD/GR3iyGdyITtO10BvEj+NvmREBOTtgzCIPN34xGDr5X3FDdIFQkS9xsSGqB+8lHKMiBYHCORXN6OxzShxREQ9s5ljnom6NT4hyHp5X1G9dIEQUa8xsSHqp3Ru1EkuxmIR8NNxMbHxUSkwPS1M4oiInMeomEB4KcW3RzlMbIhcChMbon7KsBkgkMsBAuQC9pc0oKbJAACYOSQM3iqFxBEROQ8vpRxj2j+wKq5rQXWjXuKIiAaf2SJIHUKfMLEh6qeM+M4VG/bZkCvgNDSicxufYNtnw1Ub8jzPfnsEv/90P6p0bVKH0itMbIj6KS3cH94q8Z8SV2zIFXT018hkwEXDIySOhsj52A8QYGJDniW3VIuPdhVh3cFyXP3mDuhNZqlD6jEmNkT9pFTIMSqms2yhocUgcUREZ1dY04xTVU0AxE+lw/zVEkdE5HzsVmzYZ0MexGwR8Nc1ueioRFs0PRlqpeuUKzOxIXKA9FgOECDXYFuGNodlaETdCteokRDiCwA4WKqFwcSNOskzfLqnGAfbq0+GRvpj0YxkiSPqHSY2RA5gu1HnIZajkROz66/hmGeis5rQXo5mMFlw9IxO4miIBl51ox6vbDhuvf78vHSoFK6VKrhWtEROyjaxYZ8NOauGFgOyC8WymuQwP6SG+0kcEZHzst3PhmOfyRO89MMx6NpMAIDrx8dhUnKIxBH1HhMbIgdIDvOHn5dYg8pSNHJWP5+oso7wvGREBGQymcQRETkvDhAgT7IrvxZf7ysDAAT6qPD4FcMljqhvmNgQOYBCLsOo9j6bsoZW1DRx3wNyPpuPVlkvswyN6NyGRWrg2/6BFQcIkDszmCz425rD1ut/vmyYyw6WYWJD5CAZsSxHI+elN5mx9WQ1ACDIV2XtHyCi7ikVcoyNDwIAnNG2obyhVdqAiAbI/34tsE7LHBMfhFszEySOqO+Y2BA5SEb7CyDAAQLkfHbn16FJL9ZOXzQ8AkoXawglksIElqORmytraMUbm08BAOQy4IV5oyGXu26ZMl/ZiBzEbsWmrEG6QIi6YTfmmWVoRD1iu58NBwiQO1q27ghajeIGnHdNTcJom/cyroiJDZGDJIb6QuOtBMAVG3IugiBg81ExsfFSyDFzaLjEERG5hnE2k9HYZ0Pu5qdjldjY/toQrlHj4blDJY6o/5jYEDmITCazjn2uatSjUtcmcUREoiPlOpRrxb/Hqamh8FcrJY6IyDUE+XpZx6IfKdehrf2TbSJX12ow4+l1R6zX/3rlCAR4qySMyDGY2BA5UHpskPUyV23IWdhtyjmSZWhEvdHRZ2OyCHxeJ7fx1s+nUVovDsSYlhqKa8bESByRYzCxIXIg+406G6QLhMiGXWIzIkLCSIhcDwcIkLvJq27CO9vyAAAqhQzPXjvabfY1Y2JD5EDpNk13h7hRJzmBM9pWHC7TARD/PqMDfSSOiMi1cIAAuRNBEPC3NYdhNIubNd8/KxVpEf4SR+U4TGyIHCgu2AfBvmKN6qFSLQRBkDgi8nSbj3FTTqL+SA33R0D7YJh9RfV8XieXtu5gOXbm1QIQ37MsuTBN4ogci4kNkQOJAwSCAAB1zQaUcUM3ktimo7b9NSxDI+otuVyGce2rNrXNBhTXtUgcEVHf6NqMeH79Mev1ZdeMgo+XQsKIHI+JDZGD2ffZsByNpNPYZkRWXg0AICbQGyOjAySOiMg12fbZsByNXNVrG0+iulEPAJg7MhIXu+EqPhMbIgdjnw05i+2naqx11JeMjHSb5lCiwcYBAuTqDpdp8WFWIQDAR6XAU1ePlDagAcLEhsjBOkrRAK7YkLQ225ShzeGYZ6I+GxMfBHn75wI5RQ2SxkLUWxaLgCfXHIalvT3s9xcPQVywr7RBDRAmNkQOFhmgRrhGDQA4VNrARlOShMlswZYT4uAAf7USk5NDJY6IyHX5q5UYFiWWcp6o0KFJb5I4IqKe+zS7GAdLGgAAQyL8cc+MZGkDGkBMbIgcTCaTIaO9HE3XZmKjKUlib1E9GlqMAIALhoXDS8mne6L+GJ8QBACwCLC+SSRydjVNeryy4YT1+nPzRrv164H7/mZEEkq3GSDAnapJCnZlaG7YIEo02DhAgFzR3384Dm2r+CHXdeNiMSXFvVfvmdgQDYAMu8SmQbpAyCMJgoBNx8TERiGXYfawcIkjInJ9HCBArmZPQR2+zCkFAAR4K/H4FSMkjmjgMbEhGgDpsUHWy1yxocGWV92EolqxBHJSUgiCfL0kjojI9SWE+CLUT/y3tK+oHhYL+yfJeRnNFvx1Ta71+qOXDbf2/7ozJjZEAyBco0ZMoDcAccQiXwBpMG06WmW9fAmnoRE5hEwmw/j2VRtdmwl51U0SR0R0dit/LcDJSvFvNCMuELdNSpA4osHBxIZogHT02TQbzMivaZY4GvIkm4919tdcMiJCwkiI3Mv4BJajkfMrb2jF65tPAQBkMuCFeelQyD1jHzMmNkQDxG4/m7IGyeIgz1LdqLe+4Roa6Y/EUD+JIyJyHxwgQK7g2W+PosVgBgDcOSXRbqCRu2NiQzRA0mM5GY0G38/Hq9CxddIlnIZG5FAZcYFQtn/yva+4QdpgiLrx8/EqbDhSAQAI81fjkbnDJI5ocDGxIRogtolNLhMbGiSbbMvQ2F9D5FDeKgVGxYgbdZ6uakJDi0HiiIg6tRnNeHrdEev1v145AoE+KgkjGnxMbIgGSLCfF+JDfAAAR8p1MJktEkdE7q7NaMb2U9UAxE/qxtqUQxKRY4y3KUfbz1UbciL/+fm0dVPwqSmhuHZsjMQRDT4mNkQDKKN97HOr0Yy8ag4QoIH16+katBnFBPqSERGQe0izKNFg4gABckb51U1YvjUfAKBSyPDcvFGQyTzvNYCJDdEAst2o8yA36qQBZj8NjWVoRAOBAwTI2QiCgKfWHoGhvTLk3pkpSIvQSByVNJjYEA0g20kk7LOhgWSxCNh8TNy/xlslx/S0MIkjInJPMUE+iG7fp+xgSQPLjEly3x06gx2nawAAsUE+WHrREIkjkg4TG6IBNNp2MloZExsaOAdLG1DdqAcAzEgLh4+XQuKIiNxXRzlas8GME5WNEkdDnqyxzYjnvjtqvb7smlEe/fzPxIZoAAV4q5ASJu4jcuyMDgYTP9mjgWFbhjZnJDflJBpItgME9rEcjST02qaTqGr/UOuSEZEePw2TiQ3RAOsoRzOYLDjJT/ZogGw+KpahyWTARcM9+4WNaKDZ9tlwPxuSypFyLT7YWQhALEF++uqR0gbkBJjYEA0wu/1sWI5GA6C4tsVaDjMuPgjhGrXEERG5t5HRAVArxbdQHCBAUrBYBPx1zWFY2jdk/v3FQxAf4ittUE6AiQ3RAMuw2UvkEAcI0ADYzE05iQaVl1JunXpZXNdi7W8jGiyf7y2x7qOUGu6HxTNSpA3ISTCxIRpgo2IC0DFKPresQdJYyD1tOmrTX8Mxz0SDgvvZkFRqm/T4+w/HrdefmzcaXkq+pQeY2BANOD+1Emnh/gCAExWNaDOaJY6I3Im2xYg9hXUAgMRQX6RF+EscEZFn4AABksrffzgObasRADB/XCympXK8fwcmNkSDoKMczWgWcLyCAwTIcX45WQVze5H1nBGRHrnTNJEUuGJDUsgurMPqnFIAgMZbicevGC5xRM6lV4nNSy+9hMzMTGg0GkRERGDevHk4ceKE3TELFy6ETCaz+2/KlCkODZrI1WTYbdTZIF0g5HZsy9DYX0M0eMI1aiS0N2sfLNVynD8NOKPZgr9+c9h6/dFLhyFC4y1hRM6nV4nN1q1bsWTJEuzatQubNm2CyWTC3Llz0dzcbHfcZZddhjNnzlj/+/777x0aNJGrSbdJbDhAgBzFYLJg64lqAECgjwoTbUpjiGjgdYx9NpgsOHpGJ3E05O7e/7XQOgEzPTYQt09OlDgi56PszcEbNmywu75y5UpEREQgJycHs2bNst6uVqsRFRXlmAiJ3MDI6AAo5DKYLQJHPpPD7CmoQ6PeBAC4aHgElApWFxMNpvEJQfhmfxkAcezz2PggaQMit3VG24r/23wSgLhf2fPzRkMhZ+nxb/XrVVCrFd+ghYSE2N3+yy+/ICIiAkOHDsW9996Lqqqq/vwYIpfnrVJgaKQGAHCyshGtBg4QoP7bdLTCevkSTkMjGnQcIECD5bnvjqKl/b3D7ZMTMIZJdLd6tWJjSxAEPPzww5gxYwZGjx5tvf3yyy/HjTfeiMTERBQUFOBvf/sbLrroIuTk5ECt7rppnF6vh17fOf9dpxOXco1GI4xGY1/Ds+q4D0fcF/UOz7290TEaHDujg0UADpXUYXxC0ID+PJ5/aQ30+RcEwdpfo1LIMDU5iI/1b/DfgHQ85dynhHjDz0uBZoMZOUV1TvP7esr5d2aOfAy2narB97niB1mhfl7440WpHvXY9uZ3lQmCIPTlhyxZsgTr16/Hjh07EBcXd9bjzpw5g8TERHz22We47rrrunz9mWeewbJly7rcvmrVKvj6cgdVch87KmRYXaAAAFyXZMYF0X36p0cEAChrBl45JH42NTzQggdGsnGZSApvHpHjlE4sgHlmvAnBXT/DJeozgxn4+0EFavVi2dkdaWZkhnvW+4eWlhbcdttt0Gq1CAgIOOexfVqxWbp0KdatW4dt27adM6kBgOjoaCQmJuLUqVPdfv3xxx/Hww8/bL2u0+kQHx+PuXPnnjf4njAajdi0aRPmzJkDlUrV7/ujnuO5txdfpsXq5bsBAJbAOFxxRfqA/jyef2kN9Pl/8+c8AHkAgFtmjsQVkxMc/jNcHf8NSMeTzv0Jr9M4tTUfABCcNh5XpEvfY+xJ599ZOeox+NeW06jVi39fk5KC8dRdEz1urH9HNVdP9CqxEQQBS5cuxTfffINffvkFycnJ5/2e2tpalJSUIDo6utuvq9XqbkvUVCqVQ/8xOvr+qOd47kWj4oLhpZDDYLYgt1w3aOeE519aA3X+fz5ZY708d3QMH+Nz4L8B6XjCuc9MDgXaE5sDZTpcOz5e4og6ecL5d3b9eQwKaprxzrZCAIBSLsML89Ph5eXlwOhcQ2/OX6+GByxZsgQff/wxVq1aBY1Gg4qKClRUVKC1tRUA0NTUhP/3//4fsrKyUFhYiF9++QVXX301wsLCMH/+/N79FkRuRq1UYHi0OEAgv6YZjW2eUx9LjlWhbbOODR8VE4CYIB+JIyLyXONs+iU5QIAcRRAEPLX2MAxmscx48cwUDGkfQkRn16vE5u2334ZWq8Xs2bMRHR1t/e/zzz8HACgUCuTm5uLaa6/F0KFDsWDBAgwdOhRZWVnQaPhgEKXHivvZCAJwpJx7HlDfbD5msyknp6ERSSrI1wtpEf4AxOf1NiOnXlL/fZ9bge2nxJX52CAf/P7iNIkjcg29LkU7Fx8fH/z444/9CojInWXEBeITsc0GuaVaTEkJlTYgckm2ic2ckUxsiKQ2PiEIp6uaYLIIOFSqxaTkkPN/E9FZNLYZ8ex3R6zXn756JHy9+jzI2KNwNzeiQZQeG2S9fIgbdVIfNOtN2Hm6FgAQHeiNUTH9H7JCRP0zwXY/m2KWo1H/vPnzaVTqxK1QLh4ewQ+weoGJDdEgGhLpD7VS/GeXW9ogbTDkkrafqrbWXF8yItLjpuMQOaPxCZ2JTQ77bKgfBEHAmv1lAAAvpRzPXDOKz/O9wMSGaBCpFHKMbP+EvbC2BdoWDhCg3tl0tMp6+RJ+ikfkFFLD/RHgLZYK7SuqP2/pPtHZ5FU3W1drpqaEIj6Eezr2BhMbokGW0T5AAAAOl7McjXrOZLZgy3Gxv8bPS4EpKazjJ3IGcrkM49vL0WqbDSiua5E4InJVWXmdo/ynpbIPt7eY2BANsoy4IOvljpG9RD2xr7gB9e2rfBcMC4daqZA4IiLqwHI0coSdebXWy9PTwiSMxDUxsSEaZBlxnSs2h9hnQ73AaWhEzosDBKi/LBYBWfliYhPoo8KIaA6H6S0mNkSDLCXcH75e4iftXLGh3th8VExsFHIZLhwWIXE0RGRrTHwQ5O093jlFDZLGQq7pWIUODe2r8lNSQqCQc2hAbzGxIRpkCrkMo2PEVZuyhlbUNukljohcQV51E/JrmgEAExODEeTrJXFERGTLX63EsCjxE/YTFTo06U0SR0SupmOUPwBMS2UZWl8wsSGSQLpNOVou97OhHuhYrQFYhkbkrCYkBgEALAJwsKRB0ljI9ezk4IB+Y2JDJAHbPptclqNRD2yySWwuGcHEhsgZcYAA9ZXRbMGegjoAQLhGjbQIf4kjck1MbIgkkG4z8vkQV2zoPGqb9Mhpb0YeEuGPpDA/iSMiou5wgAD11aFSLZoNZgDiag035ewbJjZEEkgK9YNGLW7mxhUbOp8tx6vQsd8fN+Ukcl4JIb4I9RP73/YV1cNi4Uad1DPcv8YxmNgQSUAul1n7bCp0bajStUkcETkz2zHPLEMjcl4yWedGnbo2E/KqmySOiFzFrxwc4BBMbIgkwgEC1BNtRjO2nRQ/yQvz98LY+CBpAyKic2I5GvVWm9FsLTeOC/ZBfIivxBG5LiY2RBLJiA2yXj7IcjQ6i6y8WrQaxbrri4ZHcF8DIifHAQLUW/uK6mEwWQCwDK2/mNgQScR+MlqDdIGQU9vIaWhELiUjLhDK9g8gmNhQT+zM6yxDm57GMrT+YGJDJJG4YB8E+aoAiKVogsAmU7JnsQj4qb2/Rq2UY+aQcIkjIqLz8VYpMCpG3Kgzr7oZDS0GiSMiZ2e7f83UFK7Y9AcTGyKJyGQy69jnmiYDzmg5QIDs5ZZpUdWoBwDMHBIGHy+FxBERUU+Mt+mz2V/cIF0g5PQa24zWcvS0CH9EBHhLHJFrY2JDJCHbcrRD7LOh3+A0NCLXxAEC1FPZhXUwt48FZ39N/zGxIZJQus0AgdyyBsniIOe0yaa/5qIRERJGQkS9wQEC1FM77cY8M7HpLyY2RBLiig2dTUldC45XNAIAxsYHIULD8gQiVxET5IPoQPHf7IGSBpjMFokjImfVMThAJgOmsL+m35jYEEkoOtAbYf5qABwgQPZsy9DmjGQZGpGr6Vi1aTGYcaKyUeJoyBnVNxtw9IwOADAqJgBBvl4SR+T6mNgQSUgmk1lXbRpajCipa5U4InIWTGyIXJvtAIF9LEejbmTl25ahccyzIzCxIZJYx2Q0ADjEPhsCoG01Ynd+HQAgIcQXQyL8JY6IiHrLfoBAg3SBkNOyG/PM/hqHYGJDJDH7jTrZZ0PA+kNnYGqfkjNnZCRkMpnEERFRb42MDoBaKb7N4gAB6k5Hf41SLsOkpBCJo3EPTGyIJGa3YsPEhgB8mVNivTx/XKyEkRBRX3kp5dYProrrWlDdvicVEQBUaNuQX90MQBwQ46dWShyRe2BiQySxiABvRLVvyHW4TAuLhQMEPNnpqiZr2crwKI11B3Micj22Y5+5nw3ZysrvLEPjmGfHYWJD5ATS2z/Va9SbUFjbLHE0JKWv9pVaL98wIY5laEQujAME6Gx+tdm/ZioHBzgMExsiJ5BhU46WW8ZyNE9ltgj4uj2xUcplmMcyNCKXxhUb6o4gCMhq769RK+UYlxAkbUBuhIkNkRPIiA+yXmafjefacboGlTqxDv/C4RHWPY6IyDWFa9RICPEFABws1cJg4kadJPZclTWI2ztMTAqGt0ohcUTug4kNkROwHSDAyWiea/XezqEBN0yIkzASInKUjrHPBpMFR8r5/E6d09AA7l/jaExsiJxAiJ8X4oJ9AACHy7Uwc4CAx9G2GLHxqLgpZ6ifFy4aHiFxRETkCONtyoy4nw0Bv01sODjAkZjYEDmJjrGgLQYz8qqbJI6GBtu3h8qtZSrXjo2FSsGnZyJ3wAECZEvsrxEnovmrlXYVG9R/fOUkchLpsUHWy+yz8Tyrc+ynoRGRexgWqYGfl9hDwQECdLKyCTVNBgDA5OQQKPkhlkPxbBI5iY4VGwDILW2QLhAadKcqG3GwpAGAuFv5SO5dQ+Q2lAo5xrQPiDmjbUN5e9M4eaadeZ3710xlGZrDMbEhchKjYzoTm0Mc+exRvrTZu+bGiVytIXI3E2zK0XJYjubRODhgYDGxIXISgb4qJIWKY0GPlutgNHMsqCcwmS34el8ZAEClkOHasdy7hsjdcD8bAsS9ynbli4lNiJ8XhkdpJI7I/TCxIXIiGXFBAAC9yYJTlRwg4Am2n6pBdaO4d81FwyMQ4uclcURE5Gi2GzBygIDnOlymRWObCQAwNSUUcrlM4ojcDxMbIidi12dT1iBdIDRovrQZGnDjhHgJIyGigRLk64W0CH8AwJFyHdqMZokjIinYlqGxv2ZgMLEhciK2Yx85Gc39NbQYsKl975owfy9cMCxc4oiIaKB07Gdjsgh8fvdQtoMDuH/NwGBiQ+RERsUGQta+Ms0XPve37mA5DO29VPO4dw2RW+MAAc9mMFmQXVgHAIgK8EZymJ/EEbknvooSORF/tRKp4WK5wvEKHfQmliu4M9sytBs4DY3IrXGAgGc7UNKANqP4Qda0tFDIZOyvGQhMbIicTEZ7OZrRLOBERaPE0dBAOVHRaF2VS48NxPAo7l1D5M5Sw/0R4K0EIA4QEARB4ohoMP162rYMjWOeBwoTGyInkx7HPhtP8GVOifXyDRO4WkPk7uRyGca3l6PVNhtQXNcicUQ0mLI4OGBQMLEhcjJ2k9GY2Lglo9mCb/aXAwC8FHJcMyZG4oiIaDDYlqOxz8ZztBhM2F8iPt5Job6IDfKROCL3xcSGyMmMjA6Eon22/aEyJjbuaNvJatQ0iXvXXDIyAsHcu4bII3CAgGfaW1gPo1ksPZzKMrQBxcSGyMn4eCkwpH2/g5OVjWg1cICAu7EbGsAyNCKPMSY+CB17Mu4rbpA0Fho8tvvXTE9jGdpAYmJD5ITGxAUBAMwWAQdKGiSNhRyrrtmAzcc69q5RY9YQ7l1D5Cn81UoMax8UcqJChya9SeKIaDDY7l8zJYWJzUBiYkPkhCYlh1gv7ymokzAScrR1B8qsJQnXjY+FknvXEHmUCYlBAACLABzkB1duT9tqxOH2svLhURqE+asljsi98RWVyAnZJTaFtec4klzNl/tYhkbkyThAwLNkF9bD0j7Zm9PQBh4TGyInFB/SOTUlp6geBpNF4ojIEY6d0eFwmQ4AMCYuEEMjNRJHRESDjQMEPEtWfmfVBfevGXhMbIicVMeqTZvRglxOR3MLHBpARAkhvghtn4S4v7geFgs36nRnu9oTG7nMvhqDBgYTGyInxT4b92I0W7BmfxmAjr1rYiWOiIikIJN1btSpazMhr7pJ4ohooOgMwMkq8fFNjwtCoI9K4ojcHxMbIic12Sax2V3APhtX9/PxKtQ2GwAAc0ZFItCXL3BEnsq2HG1fMcvR3NVpncx6eRr7awYFExsiJ5Uc5mednrK3sB5mliu4NJahEVEHDhDwDCe1TGwGGxMbIiclk8msqzZNehOOndFJHBH1VW2THluOVwEAIgO4dw2Rp8uIC4SyfadOJjbu61R7YqNSyDAxkf01g4GJDZETm2RXjsY+G1e19kA5TO0rbvPHxUEhl53nO4jInXmrFBgVI27UmVfdjIYWg8QRkaOVNbSiRi8+149LCIaPl0LiiDwDExsiJzY5xSaxyWefjatabVeGxqEBRATrAAEA2F/cIF0gNCBsxzxP55jnQcPEhsiJDY3QWKeoZBfWcSyoCzp6RmctIxwbH4S0CO5dQ0QcIODudtnuX5PG/prBwsSGyInJ5TJkJomrNvUtRpzmWFCX8/X+cuvlGydyaAARiThAwH0JgmBNbHxUcoyJC5I2IA/CxIbIyU1hOZrLMlmAdQfPAAC8lHJclREjcURE5CxignwQHegNADhQ0gCT2SJxROQo+TXNqGzUAwAmJgbDS8m324OFZ5rIyXGAgOs6Ui9DfYsRAHDpqChuzkZEdjpWbVoMZpyobJQ4GnKUnXmdH0LafjhJA4+JDZGTGxkdAH+1EgCwp6AOgsA+G1exp7pz+tmN3LuGiH7DdoDAPpajuY2dp2usl6elsL9mMDGxIXJySoXc2mRa1ahHYW2LxBFRT9Q06XG0XkxsogK8MT2NU3GIyJ79AIEG6QIhh7FYBGS1l437KASMiObAmMHExIbIBdiWo+0pYJ+NK1h38AwsEBOb68bHcu8aIupiZHQA1O39Fxwg4B6OVejQ0F6CnBYg8Ll/kDGxIXIBk9ln41IEQcBX+zqnod3AMjQi6oaXUo6MuEAAQHFdC6rbG87JdWXZ9NcMDWTp+GDrVWLz0ksvITMzExqNBhEREZg3bx5OnDhx1uPvv/9+yGQyvP766/2Nk8ijpccFWj/V28PExukdKdfhZJU4mnt8QhBSwv0ljoiInJXt2GfuZ+P6bAcHDGFiM+h6ldhs3boVS5Yswa5du7Bp0yaYTCbMnTsXzc3NXY5ds2YNdu/ejZgYjjcl6i+1UmF98Sutb0VZQ6vEEdG5rN5bYr183Tg+BxLR2XGAgPswmi3WbRnC/L0Q5SNxQB6oV4nNhg0bsHDhQowaNQpjxozBypUrUVxcjJycHLvjysrK8NBDD+GTTz6BSsXxpkSOwD4b16A3mbH2oFiGppILuGJ0pMQREZEz44qN+zhUqkWzwQwAmJIcAhnbawadsj/frNVqAQAhIZ1vuCwWC+688048+uijGDVq1HnvQ6/XQ6/vrCnV6XQAAKPRCKPR2J/wrPdj+38aPDz3jjUhIcB6eVdeDa46zxtmnn9pbDxSaW0czQgR4K3gYyAV/huQDs99zwV5yxEf7IOS+lYcLNWiuVXf7w0def6lseNklfXypMRAoLaUj4ED9OYc9jmxEQQBDz/8MGbMmIHRo0dbb3/55ZehVCrx+9//vkf389JLL2HZsmVdbt+4cSN8fX37Gl4XmzZtcth9Ue/w3DuGwQwoZAqYBRl+PlyK71VFPfo+nv/B9c4xOToWwyeHCzz/ToCPgXR47nsmUiFHCeQwmCx476sNSHLQhGCe/8H13ZHO539j2RHAm4+BI7S09Hybiz4nNg899BAOHTqEHTt2WG/LycnBG2+8gX379kHWw/W3xx9/HA8//LD1uk6nQ3x8PObOnYuAgIBzfGfPGI1GbNq0CXPmzGFZ3CDjuXe8Tyv2YF9xA6raZMiceTHCNeqzHsvzP/iqGvV4ePc2AAKiAtQYEtjM8y8h/huQDs9979TvLsbe744DAHziR+GKaYn9uj+e/8GnN5rxaPbPACyIDfLGLVdNwebNm/kYOEBHNVdP9CmxWbp0KdatW4dt27YhLq5zjOn27dtRVVWFhIQE621msxmPPPIIXn/9dRQWFna5L7VaDbW665szlUrl0D8ER98f9RzPveNMTgm1buK2v7QRV2acf9oWz//gWX+4GGaLOAVn/rgYyA2neP6dAB8D6fDc90xmSucGvgdLdQ47Zzz/gye7SAuDyQIAmJYaBi8vLwB8DByhN+evV0WcgiDgoYcewtdff40tW7YgOTnZ7ut33nknDh06hAMHDlj/i4mJwaOPPooff/yxNz+KiLphv58NBwg4E0EQ8GVOqfX69eNiJYyGiFzJsEgN/LwUALhRp6uyHfM8PS3sHEfSQOrVis2SJUuwatUqrF27FhqNBhUVFQCAwMBA+Pj4IDQ0FKGhoXbfo1KpEBUVhWHDhjkuaiIPNSExGHIZYBG4n42zOVSqxclKce+azKRgJIb64ojEMRGRa1Aq5BgTH4SdebWo0LWhvKEVMUGcFexKdubVWC9PTQ09x5E0kHq1YvP2229Dq9Vi9uzZiI6Otv73+eefD1R8RGRD463CqBhxl+rjFY1oaDFIHBF1sF2tuWFC3DmOJCLqaoLNfjZctXEtTXoTDpaKk4JTw/0QGeAtcUSeq1crNoLQ+x1Uu+urIaK+m5wcgtwy8Ql0T0Ed5o6KkjgiajOasa597xoflQJXZsQA4I7TRNRzv93P5uox3NzXVWQX1Fn7K6elsgxNSv0blE5Eg85+o06WozmDzccqoW0V5+xfPjoK/up+bRFGRB5oXEKQ9fI+rti4lF9Pd5ahTWMZmqSY2BC5mMwkm8SmkImNM2AZGhH1V5CvF9IixEmXR8p1aDOaJY6IeqpjcIBMBkxJYWIjJSY2RC4m2M8LwyLF3dsOl2nRpDdJHJFnq9S1YdvJagBAbJAPX9SIqM/Gt6/amCwCDrX3bJBzq2824OgZcZ+VkdEBCPbzkjgiz8bEhsgFTU4RV20sArCXqzaS+npfGdpLq3H9hDjI5T3bnJiI6Lc4QMD17MrvHPPMMjTpMbEhckHss3EO4t41Jdbr14/n3jVE1He/HSBAzs92/xoODpAeExsiF8TExjkcKGlAXnUzAPExSQz1kzgiInJlqeH+CPAWh4/sK6rv0zRaGly/tu9fo5TLkGnz2kzSYGJD5IIiNN5ICRPfRB8sbUCrgU2mUrAdGnAjhwYQUT/J5TKMby9Hq202oKi2ReKI6FwqtG3Ib/9wa0x8ECdiOgEmNkQuqmPVxmgWsL+EJQuDzXbvGl8vBa5Ij5Y4IiJyB7blaOyzcW5Z+Rzz7GyY2BC5KJajSWvj0Uo0tokT6S4fHQ0/flJHRA5gO0BgLxMbp7bzdGd/zVQmNk6BiQ2Ri5psM1Z4dz4Tm8HGvWuIaCCMSwiCsn264u6C2vMcTVIRBME6OMBLKbdbaSPpMLEhclGxQT6IDfIBIE7PMZgsEkfkOc5oW7H9lLh3TXyIDyazYZSIHMTXS4n0uEAAQH51M6ob9RJHRN0prmtBWUMrAGBiYjC8VQqJIyKAiQ2RS+t4Q603WZBb1iBtMB7k631l6BhWdP147l1DRI5lW2qczb3KnJLtmOfpaRzz7CyY2BC5sI6NOgFgN/tsBoUgCPjKpgzt+vEsQyMix5qS3FlqzB5K52Sb2LC/xnkwsSFyYZOS2Wcz2PYVNyC/RhzvOTUlFPEhvhJHRETuZkJSMGTtC8G2O9uTcxAEAVnt+9f4q5XIiA2UOCLqwMSGyIUlhfoiXKMGII4FNZnZZzPQvswpsV7m0AAiGggB3iqMjA4AAJyobERDi0HiiMjWqaom1DSJj8mk5BAoFXw77Sz4SBC5MJlMZq3FbtKbcOxMo8QRubdWgxnfHTwDAPDzUuDy9CiJIyIidzW5fUVeEIC9hRz77Ex+Pc39a5wVExsiFzcl2bbPhiULA2nj0Qo06sW9a65Ij4avF/euIaKBMYnP7U7Ltr9mWioHBzgTJjZELs6uz4ZNpgNq9d7OoQE3ToyXMBIicnfchNk5mS2Cte8p2FeF4VEaiSMiW0xsiFzckAh/BPmqAIhjQS0WQeKI3FNZQyt+bW8WTQjxRWYSN2MjooET4ueFoZH+AIDD5To0ta8Wk7SOlGvR2CY+FlNTQznu38kwsSFycXK5DJOSxE/2GlqMOFnFPpuB8M2+UuveNTdMiINMxhczIhpYHas2ZouAnCL22TiDX0/bjnlmGZqzYWJD5AZYsjCwBEHAl+1718hkwPWchkZEg2CS3X427LNxBjvzODjAmTGxIXIDk9lnM6D2FtWjsLYFgPhCFhvkI3FEROQJJvNDK6diMFmQXSg+DpEBaqSE+UkcEf0WExsiNzAyJgD+anFC1+78OggC+2wc6UuboQHcu4aIBktkgDeSQsVNgA+WaNFmNEsckWc7UNKANqO4X9z01DCWJDshJjZEbkAhl2FiezN7TZMeBTXNEkfkPloMJqzPFfeu8VcrcdmoaIkjIiJP0lFqbDBbsL+4QdpgPJxtGdpUlqE5JSY2RG6CfTYDY8PhCus0oqsyouHjpZA4IiLyJPZ9Nnxul9JOu8EBTGycERMbIjfBWuyB0TE0AGAZGhENPrvn9kIOEJBKi8GE/SXiZLrEUF/EBftKHBF1h4kNkZtIjw2Ct0r8J80BAo5RWt9i3WE6OcwPExK5dw0RDa64YB/EBHoDAHKK6mEwWSSOyDPtLayH0Sz2r3IamvNiYkPkJryUcoxPEN94lzW0oqyhVeKIXN9XOWXWy9y7hoikIJPJrKXGbUYLcsu0EkfkmTo+5AKAady/xmkxsSFyI7Z9NtmF3MytPywWAV/uKwEg7l0zf1ysxBERkaeanMI+G6ll2QwOmJLCFRtnxcSGyI3Y7mezh4lNv2QX1qGkTlz1mpEWhhjuXUNEErH90Go3N+ocdNpWo3WlbFikBuEatcQR0dkwsSFyI+MSgqBSiOVSXLHpn9UcGkBETiIlzA9h/l4AxF4Ps4V7lQ2m3fm16DjlnIbm3JjYELkRb5UCY+KCAACFtS3QGqSNx1XVNxvw7cFyAIDGW4lLR0VJHBEReTLbPpsmvQnHzugkjsiz2PfXMLFxZkxsiNzM5JTOkoU8HZvd++Kz7BLo2ycP3TghHt4q7l1DRNKyLTXm5MvBldWe2Mhl9v1O5HyY2BC5GdvN3JjY9J7JbMFHWYUAxKEBd01NlDYgIiL8ps8mn302g6W6UY8TlY0AgPTYQAT6qCSOiM6FiQ2Rm5mQGAyFXExomNj03uZjlSjXtgEALhwWgaQwP4kjIiISm9Y73lRnF9bBwj6bQZFlk0RO5Zhnp8fEhsjN+KuVGB0TAAA40ypDfQsbbXpj5a+F1ssLpyVJFgcRkS25XIbMJHHVpr7FiNPVTRJH5Blsxzyzv8b5MbEhckO2JQt7CxukC8TFHDujs9aup4T7YUYaP50jIucxmeVog65jcIBKIcPEpGCJo6HzYWJD5IZs+2yyizj2uac+bO+tAcTVGrmcpXxE5Dzs97PhAIGBVlrfgqLaFgDAuIRg+HopJY6IzoeJDZEbykwKhqz9PTn3s+mZhhYDvtlfBkAs57tuPPeuISLnMiomAH5e4pTGPQV1EAT22QykLI55djlMbIjcUJCvF4ZF+AMAjp7RQddmlDgi5/d5dgnajOKI5xsmxMFfzU/miMi5KBVyTGjvs6lq1KOwfTWBBob9/jUsTXYFTGyI3FRmey2wRQByWI52TmaLgA+ziqzXF3BoABE5Kds+mz0F7LMZKIIgYGf74ABvlRxj44OkDYh6hIkNkZvKtGly3MNa7HPafKwSZQ2tAIDZw8KRzBHPROSkJrPPZlDk1zSjUqcHAGQmhcBLybfMroCPEpGbsk1sOD3n3N7niGcichHpcYFQt7/J5odWA4dlaK6JiQ2RmwrzVyPCW2wsPVSqRavBLHFEzulERaN1A7bkMD/MGhIucURERGenViowLiEIAFBa32pdbSbHst2/ZnoaBwe4CiY2RG4sNUBMbEwWAfuL2WfTnfd3FlovL5iayBHPROT0bEf6s8/G8SwWwToRTeOtxKiYQIkjop5iYkPkxtICOkeBsha7K22LEWvaRzz7eSlw/QSOeCYi5zfFboAAn9sd7ViFDvUt4jTRKSmhUPADL5fBxIbIjaXaJTb8VO+3vthbglajWKJ3w4Q4aLxVEkdERHR+4xKCoWx/s707n4mNo9nuXzM1hWVoroSJDZEbC1YDcUHeAID9xQ3Qm9hn08FsEfBBVqH1+l0cGkBELsLHS4GMOLE8Kr+mGVWNbRJH5F522Qzcmcb+GpfCxIbIzWW2lyzoTRbklmoljsZ5bDlehdJ6sel21tBwpIb7SxwREVHPTbZZScguYA+lo5gtgrV0O9hXhaERGokjot5gYkPk5jITbcY+sxbb6v2dBdbLd3O1hohczCS7/WxYauwox87o0NhmAiD213CgjGthYkPk5iYlMbH5rVOVjfj1tPhGICnUFxcM5YhnInItExOD0fGemwMEHMe2v2YK+2tcDhMbIjeXEOKDCI0aAJBTWAeT2SJxRNKzHfF819QkfiJHRC5H462yjiE+XtGIhhaDxBG5B9v+mqmpTGxcDRMbIjcnk8mstdjNBjOOlOskjkha2lYjvt4njnj29VLghokc8UxErsm2HC27kH02/WUyW6yrXyF+XhgSwd5LV8PEhsgDTOKeB1arbUY8Xz8+DgEc8UxELsquzyaffTb9dfSMDo36jv6aEMhkXM13NUxsiDzAZLsmU89NbMwWAR9mFVmvL5iWKGE0RET9k5lk86FVoec+tzuKXRka+2tcEhMbIg8wJMIfIX5eAIDswjpYLMJ5vsM9/XKiCsV1LQCAmUPCkMYxnkTkwkL8vDAsUnweO1ymRVP7agP1DQcHuD4mNkQeQCaTIbN9Opq21YgTlY0SRyQN26EBCznimYjcQEc5mkUA9nLVps9MZou1TynM3wtp7K9xSUxsiDzEpOTOT588sc/mdFUjtp+qAQAkhPhi9rAIiSMiIuo/9lA6xuFynXXFa3JKKPtrXBQTGyIPMdnDX/w+2NnZW3PX1EQoOOKZiNyApz+3O4ptfw3L0FwXExsiDzEiOgAatRKAuEu1IHhOn42uzYiv9pUCAHxUCtw4MV7iiIiIHCMiwBvJYX4AgIOlDWg1mCWOyDVxcIB7YGJD5CEUchkmtvfZ1DQZkF/TLHFEg2f13lK0tL/YXz8hFoE+HPFMRO5jUvt0NKNZwP4S7mfTW0azBdntq11h/mqkhvtJHBH1FRMbIg/iiX02FouAj7IKrdcXTE2SLBYiooEwOYXlaP1xuEyL5vYPv7h/jWtjYkPkQWxf/DxlM7etJ6tRWCuOeJ6eFoohkRzxTETuhQME+mdXfuc5m5rKMjRXxsSGyIOMjgmEj0oBQNyo0xP6bFbajXhOli4QIqIBEhfsi9ggHwDAvuJ6GEwWiSNyLVkcHOA2mNgQeRAvpRzjE4MAAGe0bSitb5U2oAGWV92EbSerAQBxwT64aDhHPBORe+qYjtZmtOBwuU7iaFyH0Wyx7v8ToVEjJYz9Na6sV4nNSy+9hMzMTGg0GkRERGDevHk4ceKE3THPPPMMhg8fDj8/PwQHB+OSSy7B7t27HRo0EfXdZJs+m91uXrLwoc1qzYKpSRzxTERuy7YcrWOjSTq/Q6Va63CZKdy/xuX1KrHZunUrlixZgl27dmHTpk0wmUyYO3cumps7pysNHToUb775JnJzc7Fjxw4kJSVh7ty5qK6udnjwRNR79rXY7ttn09hmxJc5nSOeb+KIZyJyY3bP7YXu/aGVI3H/Gvei7M3BGzZssLu+cuVKREREICcnB7NmzQIA3HbbbXbHvPbaa1ixYgUOHTqEiy++uJ/hElF/jY0PgpdCDoPZ4tZNpl/mlFqn3MwfH4tAX454JiL3lRzmhzB/NWqa9MgpbsA8vkfvEbv9azg4wOX1KrH5La1WCwAICQnp9usGgwH//e9/ERgYiDFjxnR7jF6vh16vt17X6cS6UKPRCKPR2J/wrPdj+38aPDz30jrb+VcAyIgLwN6iBhTWtqC0thGRAd4SRDhwLBYBH9iUod2eGTvof4f8+5ceHwPp8NxLY1JSEL4/XIlmvRnlzTz/52MwdfbXRGrUiA1QOeyc8d+A4/TmHMqEPo5FEgQB1157Lerr67F9+3a7r3333Xe45ZZb0NLSgujoaKxZswaZmZnd3s8zzzyDZcuWdbl91apV8PX17UtoRHQe64vl2FgmVqIuGGLG+DD3mo52tF6Gd46L09+GBFjw0ChOCCIi97e9QoYvC8TnvnmJZlwY417P7Y5W0Ai8flj8jH9CmAV3DeFrhTNqaWnBbbfdBq1Wi4CAgHMe2+fEZsmSJVi/fj127NiBuLg4u681NzfjzJkzqKmpwbvvvostW7Zg9+7diIjoOpGouxWb+Ph41NTUnDf4njAajdi0aRPmzJkDlYqlKIOJ515a5zr/20/XYNEH+wAAt02Kw7KrR0oR4oBZ/OE+bD1VAwD4z61jMWfk4E9D49+/9PgYSIfnXhonKhpx1VtZAID0YAs+X3oxz/85vL01H69tPg0AeOHakbhpYtx5vqPn+G/AcXQ6HcLCwnqU2PSpFG3p0qVYt24dtm3b1iWpAQA/Pz+kpaUhLS0NU6ZMwZAhQ7BixQo8/vjjXY5Vq9VQq9VdblepVA79Q3D0/VHP8dxLq7vzPzklHAq5DGaLgL1FDW71+BTUNFuTmtggH1yaHiPpNDT+/UuPj4F0eO4H18jYYAT5qtDQYkReowwKhZLn/xz2FDZYL08fEjEg54r/BvqvN+evV1PRBEHAQw89hK+//hpbtmxBcnLPNrsTBMFuVYaIpOWnVmJ0bCAA4GRlE+qaDRJH5Di2vTV3TU3kiGci8hhyuQyZSWLfc4tJhtPVTRJH5LwMJgv2Fon9NdGB3kgMZfuDO+hVYrNkyRJ8/PHHWLVqFTQaDSoqKlBRUYHWVnGTv+bmZjzxxBPYtWsXioqKsG/fPixevBilpaW48cYbB+QXIKK+mWw39tk9pqM16U3WEc/eKjluzuSIZyLyLHbP7dzP5qwOljagzSj21HD/GvfRq8Tm7bffhlarxezZsxEdHW397/PPPwcAKBQKHD9+HNdffz2GDh2Kq666CtXV1di+fTtGjRo1IL8AEfXNpCT3S2y+yilFk94EAJg/LhZBvl4SR0RENLhsN2HmRp1ntyvPdv+a7qf7kuvpVY/N+eYMeHt74+uvv+5XQEQ0ODKTQiCTAYIA7HaDjTotFgEfZBVary+YliRZLEREUhkRrYGfWoFmvRnZhfUQBIGrEd3YZfO6NzUlTMJIyJF6tWJDRO4j0FeF4VHidJGjZ3TQtbn2rP3tp2uQX90MQCzF6PjdiIg8iVIhx4SEIABAdZMBBTXN0gbkhPQmM/a2r2bFBHojPsRH4ojIUZjYEHmwjlpsQQByXLxkwXZowN3TkySLg4hIapmJwdbL7lJq7EgHS7TQm9r7a1LZX+NOmNgQeTDbJtNdLlyOVljTjJ9PVAEQP327ZESkxBEREUlnkhsOh3GkXfm2/TWh5ziSXA0TGyIPlukmL34fZhWhowXwzqlJUCr41EZEnmt0TABUMvFJcbcLP7cPlKw82/4aJjbuhK/+RB4szF+N1HA/AEBuqRYtBpPEEfVes96E1XtLAABqpRy3cMQzEXk4L6UcSRoxsSlraEVpfYvEETmPNqMZ+4rF0uvYIB/Eh3D/GnfCxIbIw01u/7TKZBGwv7hB2mD64Ot9pWhsH/E8b2wsgv044pmIKDWgc5KtK6/IO9qBkobO/hqu1rgdJjZEHs62z2Z3vmv12QiCgPdthgZwxDMRkSjNZjDk7nwmNh1s+2umpjKxcTdMbIg8nG2TqavVYu84XYO89hHPk5JCMDKGI56JiAAg0V+ASiFO+9pT6FrP7QPJNrGx/WCP3AMTGyIPFx3og4T2GuP9JQ3Qm8wSR9RztiOeF3LEMxGRlZcCyIgNBAAU1DSjStcmcUTSE/trGgAAccHsr3FHTGyIyLpqYzBZcLBEK3E0PVNc24KfjosjnqMDvTF3JEc8ExHZykyy2c+GqzbYX9wAQ3t/DaehuScmNkT0mz0PXKPP5sOsQuuI5zumJHLEMxHRb9gmNuyzAbK4f43b4zsBIrIfIOACfTbNehM+bx/x7KWU49ZJCRJHRETkfMbFB0EuttlwMhp+szEnBwe4JSY2RISEEF9EBXgDAHKK6mE0WySO6Ny+2V+GxjZxxPO1Y2IQwhHPRERdaLyVGN3eZ3OishH1zQaJI5JOm9GMA+39NQkhvogN8pE2IBoQTGyICDKZzFqO1mIw40i5TuKIzk4QBLuhARzxTER0dpOSbEqNPbjPZl9RPQzmjv1rOA3NXTGxISIArtNnszOvFqeqmgAAExODrZ9GEhFRV/bP7Z6b2HD/Gs/AxIaIANh/guXML37vc8QzEVGPMbER7bIZnjA5mYmNu2JiQ0QAgNRwf2uvyp6COpgtgsQRdVVS14LNxyoBAFEB3rh0VJTEERERObcgXy8Mj9IAAI6Ua6FrM0oc0eBrNZixv6QeAJAY6osY9te4LSY2RASgvc+mvRZb12bCiYpGiSPq6qNdRTYjnhOg4ohnIqLz6li1sQjigBhPs6+4Hkaz+OLB/WvcG98VEJHV5BTn7bNpMZjw2Z5iAICXQo5bOOKZiKhHbEuvPLEcLSuP+9d4CiY2RGQ1yYn3s1mzvxy69hHPV4+JQZi/WuKIiIhcQ2Zy50adnpjY7OLGnB6DiQ0RWQ2PCoDGWwlAfPETBOfos/ntiOeFHPFMRNRjERpvpIT5AQAOlTag1WCWOKLB02Iw4WBpAwAgOcwPUYHe0gZEA4qJDRFZKeQyZLb32dQ2G5BX3SxxRKKs/FqcqBR7fsYnBCE9jiOeiYh6o6PU2GgWsL/Yc/psxE2nxQ/puH+N+2NiQ0R2JtuUo60/dAZtRuk/2bNbrZmeLF0gREQuyplLjQcSy9A8i1LqAIjIudi++P3f5pN46+fTSI8LxMSkYGQmhmBCYjCC28dCD4bS+hZsOiqOeI7QqHH5aI54JiLqrUk2AwR2O9lwmIFku38NExv3x8SGiOyMjg1EUqgvCmtbAAAGswU5RfXIKarHO8gHAAyJ8MfEpBBkJgVjYmII4kN8IJPJBiSej3YVwWId8ZzIEc9ERH0QG+SD2CAflDW0Yn9xA/QmM9RKhdRhDahmvQkHSxoAAClhfogMYH+Nu2NiQ0R2VAo51i6ZgR+PVGBvUR32FtYjv8a+1+ZUVRNOVTXh0/bxyxEaNTKTQsRVnaQQDI/SQOmABKTVYMZne0oAiCOeb+WIZyKiPpucEoKv95VBb7Igt1SLiUnu3XOSU1QPU/snY1NSuVrjCZjYEFEXgb4q3JQZj5sy4wEA1Y165BTVY29hHbKL6nGkTGt9sQCAqkY91ueewfrcMwAAPy8FxieKqzkTk4IxNj4IfureP92sPVAGbau4S/ZVGdEI13DEMxFRX01OFhMbQOyzcffEJov9NR6HiQ0RnVe4Ro3LRkfhsvb+lhaDCQdKGrC3sB57i+qxr6geTXqT9fhmgxnbT9Vg+6kaAOK0tVExAZiYKJavTUgKRoTm3CUBgiDgfZuhAQs44pmIqF/s+2zqsORCCYMZBPaDA9w7iSMRExsi6jVfLyWmpYZhWmoYAMBsEXC8Qoe9hfXILqxDdmEdKnV66/Fmi4BDpVocKtXif78WAACSQn2tfToTEkOQGu5n16ezu6AOxyvEEc/jEoIwJj5o8H5BIiI3lBTqiwiNGlWNeuQU1sFktjikbNgZNetNOFSqBQCkhvud98M0cg9MbIio38QVmUCMignEgmlJEAQBpfWt1h6dvYX11n1oOhTWtqCwtgVf5pQCAEL8vDAhMVgcSJAUgv/tKLAeyw05iYj6TyaTYVJyCL47dAbNBjOOntEhIy5I6rAGRHZhHcwd/TUsQ/MYTGyIyOFkMhniQ3wRH+KL+ePiAAANLQbsK65HdqHYq3OwRAuD2WL9nrpmAzYdrbSOdu4QrlHj8tHRgxo/EZG7mtye2ADA7vw6t01sbMc8T+XgAI/BxIaIBkWQrxcuGh6Ji4ZHAgDajGYcLtMiu7AeOUV1yC6stw4KsHX75AR4Kd2zVIKIaLBNTrHvs7l3VoqE0Qwc2/6ayclMbDwFExsikoS3SoGJSSHtU3lSYbEIyKtusq7oHChpQGywDxbNSJY6VCIit5EW7o9gXxXqW4zILqyDxSJALh+Yfcik0qQ3IbdM7K8ZEuHPiZoehIkNETkFuVyGIZEaDInU4LbJ3K+GiGggyOUyZCaFYOPRSmhbjThR2YgR0QFSh+VQ7K/xXKzvICIiIvIgk5I7Rx/vKag7x5GuaVce96/xVExsiIiIiDyI7Zt9t0xsbPtruH+NR2FiQ0RERORBRkQHwF8tdiPsLqiFIAgSR+Q4jW1Ga3/N0Eh/hPmzv8aTMLEhIiIi8iAKuQwTk4IBADVNBuTXNEsckeNkF9ahvb2GZWgeiIkNERERkYexHYHsTuVodvvXMLHxOExsiIiIiDyMuw4QyMqz7a9hYuNpmNgQEREReZj02EB4q8S3gbvz3aPPRttqxJFysb9meJQGIX5eEkdEg42JDREREZGH8VLKMSFR7LMp17ahtL5V4oj6by/7azweExsiIiIiDzQpyb36bLLs9q/hmGdPxMSGiIiIyAPZ9tnsLqg9x5GuYZfN72A7HIE8BxMbIiIiIg80LiEIXgrxraCrr9iI/TU6AGJ/TTD7azwSExsiIiIiD+StUmBMfCAAoLC2BZW6Nokj6rs9BXXomH8wNZWrNZ6KiQ0RERGRh7IvR3PdVZtd+bb9NUxsPBUTGyIiIiIPNcluo07X7bPpGBwgkwGTkzk4wFMxsSEiIiLyUBMSg6GQywC4bp9NQ4sBxyrE/poRUQEI8mV/jadiYkNERETkofzVSoyOCQAAnKxsQl2zQeKIes+2v4ZlaJ6NiQ0RERGRB7Pts3HFVZssm/4aDg7wbExsiIiIiDzY5GTX3qhzV74Ys0wGTEpif40nY2JDRERE5MEyk0IgE9tssKfQtQYINLQYcLy9v2ZkdAACfVUSR0RSYmJDRERE5MECfVUYFqkBABwt10HXZpQ4op7blW+zfw37azweExsiIiIiD9fRdG8RgJzCeomj6TnuX0O2mNgQEREReThX3aizI7GRy4BM7l/j8ZjYEBEREXm4zCTbxMY1+mzqmg04XtEIABgVE4hAH/bXeDomNkREREQeLlyjRkq4HwAgt1SLFoNJ4ojOb0+BbRkaV2uIiQ0RERERoXPss8kiYH9xg7TB9EBWHvevIXtMbIiIiIgIk12sz6Zj/xq5DJjI/WsITGyIiIiICL8ZIJDv3H02tU16nKgU+2tGxwYiwJv9NcTEhoiIiIgAxAT5ID7EBwCwv6QBJXUtEkd0drYrSty/hjowsSEiIiIiAMDsoREAAIPJgns+yHbazTq5fw11h4kNEREREQEAHp4zFMlh4nS0k5VNWPLJPhjNFomj6qpjcIBCLsPEpGCJoyFnwcSGiIiIiAAAwX5e+N/CTAT5ij0r20/V4Ol1RyAIgsSRdapp0uNUVRMAsb9Gw/4aaterxOall15CZmYmNBoNIiIiMG/ePJw4ccL6daPRiL/85S9IT0+Hn58fYmJicNddd6G8vNzhgRMRERGR4yWH+eGdOyZApZABAFbtLsaKHQUSR9Vpd35nfw33ryFbvUpstm7diiVLlmDXrl3YtGkTTCYT5s6di+bmZgBAS0sL9u3bh7/97W/Yt28fvv76a5w8eRLXXHPNgARPRERERI43OSUUr9yQYb3+wvfH8OORCgkj6pSVX2O9zMEBZEvZm4M3bNhgd33lypWIiIhATk4OZs2ahcDAQGzatMnumH//+9+YNGkSiouLkZCQ0P+IiYiIiGjAzR8Xh4KaFvzrp1MQBOCPnx3AF/dPRXpcoKRxdexfI/bXcMWGOvWrx0ar1QIAQkLO/kel1Wohk8kQFBTUnx9FRERERIPsT5cMwbVjYwAArUYz7vkgG+UNrZLFU9XYhtPt/TXpsYHwV/fqM3pyc33+axAEAQ8//DBmzJiB0aNHd3tMW1sbHnvsMdx2220ICAjo9hi9Xg+9Xm+9rtPpAIj9OkZj/0cMdtyHI+6LeofnXlo8/9Li+ZceHwPp8NxLy9Hn/4VrRqC0rgU5xQ2oatRj0co9+PTeSZIkFTtPVVsvT04Kdtq/Mf4bcJzenEOZ0McxF0uWLMH69euxY8cOxMXFdRvEjTfeiOLiYvzyyy9nTWyeeeYZLFu2rMvtq1atgq+vb19CIyIiIiIHajIC/5erQI1eHCgwIsiCe4db0D5fYNB8kS/Hr5ViwdHvRpgxIsh5prXRwGhpacFtt90GrVZ71nyiQ58Sm6VLl2LNmjXYtm0bkpOTu3zdaDTipptuQn5+PrZs2YLQ0LM3dnW3YhMfH4+amprzBt8TRqMRmzZtwpw5c6BScRzgYOK5lxbPv7R4/qXHx0A6PPfSGqjzn1/djJve3Q1tqwkAcMfkeDx15XDIZIOX3Vz6xg7k17RAKZdh7xMXws9JS9H4b8BxdDodwsLCepTY9OqvQRAELF26FN988w1++eWXcyY1p06dws8//3zOpAYA1Go11Gp1l9tVKpVD/xAcfX/Uczz30uL5lxbPv/T4GEiH515ajj7/w2KCsPyOibjrf7thNAv4eHcJUsI1WDSj6/vBgVCla0N+TQsAICMuEEH+PoPyc/uD/wb6rzfnr1fDA5YsWYKPP/4Yq1atgkajQUVFBSoqKtDaKjaRmUwm3HDDDdi7dy8++eQTmM1m6zEGg6F3vwUREREROZWpqaF46brOMdDPrT+KzUcrB+Vn7yqw3b+GY56pq14lNm+//Ta0Wi1mz56N6Oho63+ff/45AKC0tBTr1q1DaWkpxo4da3fMzp07B+QXICIiIqLBc8OEODx0YRoAQBCA33+2H4fLtAP+c7Pyaq2Xp6YysaGuel2Kdi5JSUnnPYaIiIiIXNvDc4aisLYZ3x06gxaDOAZ67ZIZiAr0HrCfuTtfTGxUChkmJAYP2M8h19WvfWyIiIiIyPPI5TL848YxGJcQBACo1OlxzwfZaNabBuTnVerakF/TDAAYExcEXy/nHBpA0mJiQ0RERES95q1S4N27JiI+RGziP1Kuwx8+2w+zxfHVO7vyO8vQ2F9DZ8PEhoiIiIj6JMxfjZULM6HxFldQNh+rwvPrjzr85zCxoZ5gYkNEREREfZYWocHyOyZAKRf3s1n5ayE+zCp06M/oGBzA/ho6FyY2RERERNQv09PC8ML80dbrz6w7gp+PVznkvs9oW1FYK+5fMzY+CD5eCofcL7kfJjZERERE1G83ZybgdxekAgAsAvDQqn04Wq7r9/3uzuf+NdQzTGyIiIiIyCH+fOkwXJEeBQBobh8DXaVr69d92u1fw8SGzoGJDRERERE5hFwuw2s3jcWY+CAAwBltG+75YC9aDH0fA72rQExsvBRyjGd/DZ0DExsiIiIichhvlQLv3TURsUHiGOjcMi3++NmBPo2BLm9oRVFHf01CELxV7K+hs2NiQ0REREQOFa5RY+XdmdCoxTHQG49W4u8/HOv1/XDMM/UGExsiIiIicrihkRq8dft4KNrHQL+7vQCf7C7q1X3YJzYhDo2P3A8TGyIiIiIaELOGhuPZa0dZrz+19gi2nqzu8fdntSc2Xko5xiewv4bOjYkNEREREQ2Y2ycn4t6ZyQAAs0XAkk/24URF43m/r7S+BSV1rQCAcfHsr6HzY2JDRERERAPqsctHYO7ISABAk96ERe9no6rx3GOgbfevmZrK/ho6PyY2RERERDSgFHIZXr9lLNJjAwEAZQ2tuPeDvWg1mM/6PVkcHEC9xMSGiIiIiAacr5cSKxZMREygNwDgYKkWD39xAJazjIHeZdNfM7Z9Xxyic2FiQ0RERESDIiLAGysWZsK/fQz0D4cr8MqPJ7ocV1LXgtJ6sb9mQkIw+2uoR5jYEBEREdGgGREdgH/fNg7tU6CxfGsePttTbHcM96+hvmBiQ0RERESD6sJhEVh2TecY6L+uOYxfT9dYr++yGRzA/Wuop5jYEBEREdGgu3NqEhZNF8dAmywCfvdxDk5VNkIQBOuKjVopx9iEIAmjJFfCxIaIiIiIJPHklSNwyQhxDHRjmwl3v5+NAyUNKGto769JDIZayf4a6hkmNkREREQkCYVchjduGYtRMQEAgNL6Vty1Yo/161PZX0O9wMSGiIiIiCTjp1ZixYJMRAWIY6Ab9Sbr16ZwY07qBSY2RERERCSpqEBvrFg4Eb5enWVn3io5MuICJYyKXA0TGyIiIiKS3KiYQPz71s4x0NNTw9hfQ72ilDoAIiIiIiIAuHhEJN6/exK2HK/CPTOSpQ6HXAwTGyIiIiJyGrOGhmPW0HCpwyAXxFI0IiIiIiJyeUxsiIiIiIjI5TGxISIiIiIil8fEhoiIiIiIXB4TGyIiIiIicnlMbIiIiIiIyOUxsSEiIiIiIpfHxIaIiIiIiFweExsiIiIiInJ5TGyIiIiIiMjlMbEhIiIiIiKXx8SGiIiIiIhcHhMbIiIiIiJyeUxsiIiIiIjI5TGxISIiIiIil8fEhoiIiIiIXB4TGyIiIiIicnlKqQP4LUEQAAA6nc4h92c0GtHS0gKdTgeVSuWQ+6Se4bmXFs+/tHj+pcfHQDo899Li+ZceHwPH6cgJOnKEc3G6xKaxsREAEB8fL3EkRERERETkDBobGxEYGHjOY2RCT9KfQWSxWFBeXg6NRgOZTNbv+9PpdIiPj0dJSQkCAgIcECH1FM+9tHj+pcXzLz0+BtLhuZcWz7/0+Bg4jiAIaGxsRExMDOTyc3fRON2KjVwuR1xcnMPvNyAggH9YEuG5lxbPv7R4/qXHx0A6PPfS4vmXHh8DxzjfSk0HDg8gIiIiIiKXx8SGiIiIiIhcntsnNmq1Gk8//TTUarXUoXgcnntp8fxLi+dfenwMpMNzLy2ef+nxMZCG0w0PICIiIiIi6i23X7EhIiIiIiL3x8SGiIiIiIhcHhMbIiIiIiJyeUxsiIiIiIjI5UmS2Lz00kvIzMyERqNBREQE5s2bhxMnTtgdIwgCnnnmGcTExMDHxwezZ8/GkSNHrF+vq6vD0qVLMWzYMPj6+iIhIQG///3vodVq7e7nhRdewLRp0+Dr64ugoKAex5ibm4sLLrgAPj4+iI2NxbPPPgvbOQtnzpzBbbfdhmHDhkEul+OPf/xjn87FYHOHc2/r119/hVKpxNixY3t8/1Jxh3O/cOFCyGSyLv+NGjWqbydlkDn7Y9DW1oaFCxciPT0dSqUS8+bN6/a4rVu3YsKECfD29kZKSgqWL1/eq/MgpcF6DAoLC3HPPfcgOTkZPj4+SE1NxdNPPw2DwXDeGPn877zn3pYrPf8D7nH++RowsI+BJ7wGDDhBApdeeqmwcuVK4fDhw8KBAweEK6+8UkhISBCampqsx/z9738XNBqN8NVXXwm5ubnCzTffLERHRws6nU4QBEHIzc0VrrvuOmHdunXC6f/f3v2ERPW9YQB/BR0xJkcMmkkCK7L/UdFCg2pKQo2IFhGpFLYIVLIkg2gRTdIii2gTVBg2tYggSqlFUhBai9wIAw66KE0jU5GC0alQ1Hm+C5n53XHU35gzd+65PR+YRfce7zk8b5wz58q9dnfj3bt3yMnJwdGjR8P6unLlCm7fvo2amhrYbLaoxjcyMgK73Y7i4mJ4vV68ePECS5cuxa1bt0Jtent7ce7cOTx+/Bjbt29HdXX1onPRgxmyD/L5fFizZg0KCgqwbdu2v85EL2bI3ufzYXBwMPT59u0bMjMz4XK5Fp2PHoxeg1+/fqGiogL19fUoLCzEkSNHItp8+fIFS5YsQXV1Nbq6uvDgwQOkpKTg+fPnf52LnvSqQXNzM06dOoU3b96gp6cHL1++xPLly3HhwoV5x8f539jZB6k2/wPmyJ9rQHxr8C+sAfGWkI3NTMPDwxARvH//HgAQCATgcDhQV1cXajM2NgabzYb79+/PeZ1nz57BYrFgYmIi4pzb7Y76y8Xdu3dhs9kwNjYWOnb9+nVkZWUhEAhEtHc6ncosbDOpnP3x48dx+fJluFwuZRY2LZWzD2pqakJSUhL6+vqi6sNojFYDrbKyslkXtYsXL2LDhg1hx8rLy5GXl7fgPoxAjxoE3bx5E6tXr553PJz/1che9fkfUDv/IK4B02JVA61/ZQ2INUM8YxP89V1mZqaIiPT29srQ0JAUFBSE2qSmporT6ZSPHz/Oe5309HRJTk5e1Hja2trE6XSG/VGlwsJCGRgYkL6+vkVd22hUzd7tdktPT4+4XK5F9ZdIqmav1dDQIAcOHJDs7OxF9Z0oRqtBNNra2sLGJzJdp/b2dpmYmIh7/7GmZw1GRkZC/cyF87/xszfD/C+ibv5aXAP+d51Y1CAaZlsDYi3hGxsAUlNTI7t375YtW7aIiMjQ0JCIiNjt9rC2drs9dG6mnz9/yrVr16S8vHzRYxoaGpq1b+3YzEDV7D9//iyXLl2SJ0+e6PJFMh5UzV5rcHBQmpub5fTp04vuOxGMWINozFWnyclJ+fHjhy5jiBU9a9DT0yN37tyRioqKecfE+d/Y2Zth/hdRN38trgHTYlmDaJhpDYiHhG9sqqqqpKOjQ54+fRpxLikpKezfACKOiYiMjo7KoUOHZNOmTQu+g7N582axWq1itVrl4MGD8/Y923GVqZj91NSUlJaWSm1traxbt25B/RmJitnP9OjRI8nIyJjz4UajM2oNomGW+UmvGgwMDEhRUZEcO3Ys7EsY53+1sjfL/C+iZv4zcQ2ITw2i8S/MUX8robc7zp49K69evZIPHz7IypUrQ8cdDoeITO9KV6xYETo+PDwcsUv1+/1SVFQkVqtVmpqaJCUlZUFjeP36dehXd2lpaaH+Z+7Mh4eHRSRyF68qVbP3+/3S3t4uHo9HqqqqREQkEAgIAElOTpa3b99Kfn7+gsahN1Wz1wIgDx8+lJMnT4rFYllQ30Zg1BpEY646JScny7JlyxY0hkTSqwYDAwOyf/9+2bVrl9TX14ed4/yvVvZmmP9F1M1fi2tAfGoQDbOsAXGjz6M84QKBAM6cOYOsrCx8+vRp1vMOhwM3btwIHRsfH494eGtkZAR5eXlwOp34/fv3vH0u9CHqjIwMjI+Ph47V1dWZ4uFR1bOfmpqC1+sN+1RWVmL9+vXwer1hbzYxGtWz12ppaYGIwOv1RnVtozB6DbTme3B048aNYccqKiqUeXBUzxr09/cjJycHxcXFmJycjGp8nP+Nm73K8z+gfv5aXAPiUwMts64B8ZaQjU1lZSVsNhtaW1vDXhv458+fUJu6ujrYbDY0NjbC6/WipKQk7HV7o6OjyM3NxdatW9Hd3R12He1/oK9fv8Lj8aC2thZWqxUejwcejwd+v3/O8fl8PtjtdpSUlMDr9aKxsRHp6ekRr5wMXmvnzp0oLS2Fx+NBZ2dnjNOKLbNkr6XKW3HMlP2JEyeQm5sbw3T0YfQaAEBnZyc8Hg8OHz6Mffv2hX4uKPiqz/Pnz6OrqwsNDQ1KvepTrxp8//4da9euRX5+Pvr7+8PazIfzv/Gz11Jl/gfMlT/XgPjUADD/GhBvCdnYiMisH7fbHWoTCATgcrngcDiQmpqKvXv3ht0ZCN4tmO3T29sbaldWVjZrm5aWlnnH2NHRgT179iA1NRUOhwNXr16NuGMx23Wzs7NjkFD8mCV7LVUWNrNk7/P5kJaWhvr6+ljEoisVapCdnT3rz2m1trZix44dsFgsWLVqFe7duxeLeHShVw3cbvecbf4fzv/Gzl5LlfkfME/+XAPiWwOzrwHxlgTM8Sd9iYiIiIiIFJHwt6IREREREREtFjc2RERERESkPG5siIiIiIhIedzYEBERERGR8rixISIiIiIi5XFjQ0REREREyuPGhoiIiIiIlMeNDRERERERKY8bGyIiIiIiUh43NkREREREpDxubIiIiIiISHnc2BARERERkfL+A5TfCub9KXDJAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot\n", "\n", "# Create the line chart using Seaborn with connected lines\n", "plt.figure(figsize=(10, 6))\n", "\n", "sns.lineplot(data=df, x='time', y='tmp', linewidth=2)\n", "\n", "# Add labels and title\n", "plt.xlabel('')\n", "plt.ylabel('')\n", "plt.title('Monthly temperature 1993-2022 (°C)')\n", "\n", "# plt.xticks(rotation=45)\n", "plt.grid(True)\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "crp", "language": "python", "name": "crp" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" } }, "nbformat": 4, "nbformat_minor": 5 }