{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Plot SWE From Pits \n", "\n", "**Goal**: Calculate SWE from the Density Profiles and plot them\n", "\n", "**Approach**:\n", "\n", "1. Grab All the site ids and dates associated to all the density data\n", "2. Calculate SWE for each layer and sum them up \n", "3. plot it!\n", "\n", "\n", "## Process \n", "\n", "\n", "### Step 1. Grab all Density profiles and Calculate SWE" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from snowexsql.data import LayerData\n", "from snowexsql.db import get_db\n", "from snowexsql.conversions import query_to_geopandas\n", "import geopandas as gpd\n", "import matplotlib.pyplot as plt\n", "from geoalchemy2.shape import to_shape\n", "import pandas as pd\n", "\n", "# Connect to the database we made.\n", "db_name = 'db.snowexdata.org/snowex'\n", "engine, session = get_db(db_name, credentials='./credentials.json')\n", "\n", "q = session.query(LayerData).filter(LayerData.type == 'density')\n", "df = query_to_geopandas(q, engine)\n", "\n", "# Convert density to float\n", "df['value'] = df['value'].astype(float)\n", "\n", "# Calculate SWE\n", "swe_lambda = lambda row: row['value'] * (row['depth'] - row['bottom_depth']) / 100\n", "df['swe'] = df.apply(swe_lambda, axis=1)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Step 2. Prepare SWE as point data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Prepare the data to be a single point by summing the SWE by site and date\n", "point_swe = gpd.GeoDataFrame(columns=['date', 'swe', 'geometry'])\n", "\n", "sites = df['site_id'].unique().tolist()\n", "\n", "# Loop over data by site and date\n", "for site in sites:\n", " ind1 = df['site_id'] == site\n", " dates = df['date'][ind1].unique().tolist()\n", "\n", " for date in dates:\n", " # Grab all density at this site and date\n", " ind2= df['date'] == date\n", " \n", " profile = df[ind1 & ind2]\n", " \n", " # Check if there is data on this date/site\n", " if len(profile.index) > 0:\n", "\n", " # sum the swe column and assign data to a dictionary\n", " data = {'swe': profile['swe'].sum(), 'geometry': profile['geom'].iloc[0], 'date': date}\n", "\n", " # Add the data to a dataframe\n", " point_swe = point_swe.append(data, ignore_index=True)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Step 3. Plot the SWE as points coloring by mass" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "tags": [ "nbsphinx-gallery", "nbsphinx-thumbnail" ] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhkAAAE4CAYAAAD2AuGFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3RUxdvA8e+TRhJCQoeQAKH3Jh0lSBGQXhRQUFRUigVBAfG1/lQs2EARAZGqoAhSRJp0kE7ovXcIHRLSn/ePexOXZVOQJLT5nHMPu3Nn5s69Wd3ZqaKqGIZhGIZhpDe3210AwzAMwzDuTaaSYRiGYRhGhjCVDMMwDMMwMoSpZBiGYRiGkSFMJcMwDMMwjAxhKhmGYRiGYWQIU8kwjAwmIrVFZLKIHBORGBG5LCLrRORDEQm83eVzJCLvi0iq89pFZImIqIj8k8z5Mfb5Y+lfylvjUPbE46SIzBWRmg5xxorIIYf3IfazKZrGa2QRkT4isllErth/810iMk5ESthxBtrXL+SUtqYdvtZFvoNFJF5EAhzKqckc0//jIzKMdONxuwtgGPcyEXkdGAwsBt4GDgB+QB3gRaAa8OhtK+CtuQLUFpHiqrovMVBEfIHH7PN3qi1Ad/t1CNbfZqmIVFHVncCHwBCH+CHAe8AKrL9haiYBjYHPgdWAO1AGeBwoC+wFltlxQ4GJDmlDgUigioj4qepVp3ObVfWSQ1g40MpFGc6noZyGkaFMJcMwMoiI1MeqYAxR1T5Op/8SkU+wvnRSysMTiNM7c9W8LUA+oAvwvkN4O/vfeUDtTC5TWl1R1dX269Uisho4CPQEXlXV/f81Y7u1oy3wmqo6VlTmAF+JSGIL8jrgGq4rGROBp7Aqo/PtfLMCDwDDnC4Z43AvhnFHMd0lhpFxBgBn7X9voKoRqjo28b3dJK8i0ktEPheRE0A0kF1E8ojICBHZIyKRInJURH4RkSDHPBO7O0SkhIjMFpGrInJYRN51+HJLjFtFRJaLSJSIHBeRdwC5yXucgFXJcPQ0MA2IcI4sIh52N8EuEYkWkRMi8qWIeDvF+VBE9ttlOysiK0TkIYc4nURkkYiE2/cYJiJdb7LsSVT1EFaLQHE7/6TuEhF5GKslCmCBQ3fEw8lkl9P+91Qy10qw/43BauUIdbgvN+BBYCGw1vEcVoXDA1h6UzdnGLeRqWQYRgYQEQ+gHrDA/jK5Gf8HlMTqTmkLRGF9cUUBA4GmQD+gBLDS8QvawR/AIqANMB34AEj6EhaR3Pb53Hb4S3a+z91kWScARUWkjp1vAaAhMD6Z+BOxuiZ+AZoDnwDdgJ8d4gwA+gBDgSbAs1hfujkd4hQFfgc62/c4C/hRRHrcZPmxyx1g53/RxemNWM8H4FWs1pnadrgru4DLwKci0kVE8qVw6WVAKRHJa7+vAOQAltuHYyUjFFA73Ln8Hi6Om60wGkb6U1VzmMMc6XxgdSMo8ImLcx6Oh0N4iJ1mIyCp5O8OFLTjt3UIf98Oe9Yp/lZgvsP7j4EYoKBDWFaslhdNw/0tAVbYr5cBP9iv+wNHsH7AjAWOOaSpa5ftaae8Otvhle33fwLTbuJZu9nPchTWeIU0ld3hb1AMqyKmQBs7zljgkEOah+3zjdJYppZYLSNqH/uB74DSTvHq2+cfs9+/AuyzXzfGqlhmsd8vBbY7pR/rcA3n443b/d+BOcxhWjIMIxOJSH4g1vGwWz0cTVfVG8ZgiEhPe7bCVSAO68scoJSLS812er8NcJzFUBtYrapHEwNUNQKrReBmjQc6iEgWrK6Sn9XuEnDSFKti87vjL27sMQf8+6t9HdBMRD4WkYdExMs5I7s7aJKIHOffZ/k8rp+FKw86pNuH1RXRQ1XTZUaGqs7CqjS2A77FaiHpBYSJSCOHqKuxnknivYfyb0vFKqxKUE372dbAdVfJGaC6i2NCetyLYdwKM/DTMDLGOaxfoYWcws9ifQGA1R3ygou0J50DROQVrO6Dr7C6Si5g/YJfDbjqLnGeWRDtFC8Qq+Lh7LSLsNRMwfoifRcoR/KDWfMCXrgYq2HLZf87COvZdQHeAq6KyO9AP1U9KyJ+wAKsGRhvYrUSxGAN2kxrd89mrEqJYt3zcVcVu1thV9r+sA9EpBbwN/Ap1qwiVPWaiKzn30pGXax7RlWviMhm/u0m8ebfGSmOYlV1fXqW3TDSi6lkGEYGUNU4EVkGPCIiXmqPy1DVOGA9gIi0SC65i7BOwEJVfT0xQESK3EIRT2J16ThLafyAS6p6SURmYH3hr1drCqgriRWvusmcP2HnFwt8Bnxmt/y0wKpc+QIdsVphCgN1VXVFYmIXLUIpuZrZX8yqulpE5mO16DhaCgwQa52OfFw/5iJxXEbiZ8JVJcMw7limu8QwMs7nWAMrP0uHvHyxmvYdPXsL+a0CaolIwcQAe4pky/+Y33dYXS2fpxBnLtav8QBVXe/iOOGcQFVPqeqPWC0A5e1gX/vfpOchIjmA1v+x7GkRbf/rk1pEEclmP0vncHeswbrOLVXLsP5f/BZwWlX3OpxbgVWpaoA1VuOGZ2QYdzLTkmEYGURVF4rIm1izDCpijV04iPVFWxKrdSIC1y0XzuZi/dp9C2tqYwOsBa/+q6+xxgjMF5H3sb5E+2Gt23DT7BaFFanEWSIik7DGZHyFdR8JWGMXmgEDVHWP3SqyGWsA7AWgCtav/xF2Vv9gzd4YJiLvYQ1YfRurKyrgv5Q/DfZgjYN5TkTOYz2v3arqasGxUsBc+16XYI2ZCMTqnimP9dwdrQTisSp4U53OLcdavK0+MCaZsnnZXTHOIlV1Syr3ZRgZylQyDCMDqernIrIS6I011iAPVpfBbuBXrFkZ8WnI6n9Adqypnd5YTexNSNvqk67KdVZEGmKtajkOqyvjB6z/J7z7X/JMoy5YMyiew5qqGw0cwlq4K3E8yDKscR0vYbVaHMFqIfnYLnu4iLQFvsSaxnrCvo+cWKtypjtVPSciL2NNr12KNbunPlYlwtk+rPEzjbAqgrmBq8Am4HFV/d0p7ysiEoY1TmO507nTIrIXqwUkua6SPFgtU86282/rj2HcFpLOY50MwzAMwzAAMybDMAzDMIwMYioZhmEYhmFkCFPJMAzDMAwjQ5hKhmEYhmEYGcJUMgzDMAzDyBCmkmEYDkTkMRGZam+Pfk1EdovIJyKSzUXcHCLyo70VeYSI/C0iFZziVBORkfbW5pEickREfna1WqeIuNnboB+ytzjfLCLtb6Ls3iIyWERO2mVfJSKhLuL1FZFZdjy118lI6zVKisgQEdlib7F+UkRmikglF3F9ReQDsbanvybW9vTjRSQkrddLQ3ncReQNsbZ9Py0iV0Rko4h0E6et7e345URkvl32cyIyRkRyOsVJ189AKuUvKCK/i8glEbksItNEpJBTnIYiMlFE9tvl2S8iw+XfnVtTu0awiHxrfx4i7b95iIt4af6sGkZamUqGYVzvDayFkd7CWgBqONaeGAscv7RERLBWuGyKte5De8ATWCwiwQ75dcLaz2Mo8CjW0tsPAOvFYbVN24dYu6h+Z8ddDUwRkWZpLPtorL1Q3sVaivskME9EKjvFewFrH5H/shlYY6z1IcZhLR7VC2udhtUiUtUp7o9YC3yNwlps622sJbIXirX/SHrwsfPdhrUXTBtgsX3N61ZaFWsb+iV2msew1uFoBPzpVCFJ78+ASyLiCywCSgNdgaew1sNYLNevGNoDa1+Xj+xrfQK0wnrmaXmOxYEOWAub3bBNvIOb+awaRtrc7m1gzWGOO+kA8rgIexprVc4GDmGt7bD6DmEBWBuTDU0lv8JYK13+zyEsL9bCVB84xV0IbElDuSvhtMU71sJau4GZTnHdHM4r8P5NPJ/cOG1Db9/3BWC8Q5gv1gqZg5ziNrWv2SSd/l7uQE4X4T9hLXrm4xD2NdZuqNkdwhL3BWmXUZ+BFMreG6syU9whrIj93PqmUp7Ecj+Xhuu4ObxO3BQuJI2f/Rs+q+Ywx80cpiXDMByoariL4HX2v0EOYa2AE6q62CHtJaxftq0dwm7IT1UPA+FO+TXB2qF0olP0iUCFNDRZt8Lay+NXh+vEAZOBJmJtFZ4Y7mob9jRR1bOqqk5hl7CW3Xa8H3f7uOyUxUX733T5f4+qxquq846zYP3NsmBVihK1AmaramIZUNVlWCuKpvg34xY+AyloBaxW1X0O6Q9iLTP+X8rjUlr/3jfxWTWMNDOVDMNIXT37X8fdRcvheqv07UChlJqxRaQMVsuFc37RWEtSO+cHUDaVMpYDDqpqpIv0XlhN5hnCHtNQHof7UWtPjwnAqyJSX0T8RKQcMBhrX5KFGVUeWz2sCs1Ju4w+WK0Eyf3NUnu+6foZSEP6/1KedJfMZ9Uw0sxUMgwjBSIShLVvyN96/dbgObG6CJwl/qrOkUx+Hlh7hIRjjaFwzO+icyuBQ345SVlq5Ukt/a34FhDgG6fwZ4E/sMYdXMH6QvUEHlHVmIwqjIg0wRqD8IXdmgPW30NI/hkl+3zS+zOQxvTJprUHoH6D9cX/X8bVpEkKn1XDSDOzQZphJMP+JToDq4/8VrZVd/QdUAdorqquvmBSK5M71pdlooRb6f5I4TpuXP8jxOV1RGQg8CTQzbHZ3/YR1oZob2A17xfC2sBsjojUU9WIZK4tWF0tiVTTtokcIlIWmIQ1+POzVKKnJb9b/gzYX9aO4l1UJm8mr0lY3RcPOlSi0vU6tlv6rBoGmJYMw3DJbl6fBRTFGqR4zCnKBVz/2szpcN45z0+xZkA8p6rzXeSX3f6CdZVf4q/jhVhjLxKPxB1TUyuPq3ELKfnJ6To/OUcQkR5YO8u+rao/OZ0rhzU7oa+qfqmqy1R1ItYsk6pYAxCT09Xp2vvTUmARKQosAA4CbR2/gLG6TpTkn9ENzycdPwOxTkc9h/PJpXf1+XHDmtXTCGijN27jntx1bloqn1XDSDPTkmEYTkTEE2sL8WpYTftbXUTbjjWd01lZ4IiqXnXK8/+wtgl/RVUnJJNfFqAY14/LSOyb32H/2x1wXK/hhEP6tiLi6zQuoywQw41jPVLzPtYv2URnHU+KyFPA98CXqvqxi/SJa0WscwxU1b0ichEok8K1ZwHVHd5Hp1ZYe8roQqyBpk1U9boBp6oaKSKHsMZBOCuLtX27Y37p+Rmo7nR+t0P65Mqzw0X4D0BH4DFVdTWmJbnr3JQ0fFYNI+1u9/QWc5jjTjqwWvd+A64BDVOI1wbrl3E9hzB/4BzwrVPcV+24b6WQX16sysB7TuF/A1vTUO4q9jW6OoR5YPXbz0omzU1PYbXTtcXqPhiZQpyH7byfdwovaYe/k45/szzALqwWjOAU4n2D1UIQ4BD2kF2e9hn5GUgm/Wv2cyzqEBaC1QrxulPcL7Gmkj51i88q2Smsaf2smsMcN3OI6q102RnGvUVEhmMtfvQx8KfT6WNqN5nbTdcrgIJYC05dAAYCFYFKqnrUjtcJ+AWYB3zglN9lVU36xWo3Ub+GtQjURqxfrd2BVqrqXBZXZZ+MNRW2H9YXbk+sRbnqqOpGh3jVsL7M3LCmvE7B+lIF+EtvnKHieI1QYD7Wr/BXsL74EkWrapgdzx3YgDWj4yNgPdaYjLexKgUVVfVIaveUGrtLYyXWr//ngANOUXao3aphD+DcbJf9E6w1LT7HmoFSR+0xJ+n9GUih7Fnt8lzDei6KtSBbNqznc9WONwD4FKvLapRTNuGqmmp3kog8Zr9saN9bL6wBneGqutSOk+bPqmGk2e2u5ZjDHHfSARzC+p+9q+N9p7g5sf7Hfx6IxGqur+QUZ2wK+S1xiuuO9WVzGKuLYAtW03hay+4DfAWcwlqIag3wsIt4KZUpJJVrvJ9C2kNOcXNh/QLfi/VFehSrUlMqHf9eISmUR53vH6sbZwEQgVUpGAvkysjPQCrlLwRMxermuYI1WyTEKc6SFMozNo3XSfUzeDOfVXOYI62HackwDMMwDCNDmNklhmEYhmFkCFPJMAzDMAwjQ5hKhmEYhmEYGcJUMgzDMAzDyBBmMa5U5M6dW0NCQm53MQzDMAwj02zYsOGsqua51XxMJSMVISEhrF+/PvWIKYiIiGDIkKGMGDWOS5cv4evjS8cOrXn3nXfIkSO1PZQMwzAMI3OJyOH0yMd0l2SwsLAw8gUWZPi4v6nQ5Aua91hIjTYjmLnwIPkLFGLGjBm3u4iGYRiGkSFMS0YGOnfuHA/WrU/t1l9RtvZzSeE5A8tSsHQjDm+fw+MdH2fdmpVUqlTpNpbUMAzDMNKfqWRkoAED3iSweIPrKhiqyon9y9mx+ieunD9M1hyFaNuuPVs2b8LPz+82ltYwDMMw0pfpLslAU6bOpHL9vknvo69dYsawR1j0a3eyF6lCxdbvUK7J61yK9SNfYBDjx4+/jaU1DMMwjPRlKhkZKCLiMjnyWztaa0ICf/3Ylqz5i9P2w20Uqd6BPStGs3bK62TJWYDsJUJ5vuer5MhbgKeffpp58+YRFxd3m+/AMAzDMP47012SgdzdPYm6Eo5P1lwc3bOQyIizPNJlAVGXTvPnZw9RsHYnWg45iJefNcMkPi6GI/9MZtJPvZj8+3T8svrw9ltv0ee1VxGR23w3hmEYhnFzTEtGBipSJISda8cCsGPVaEo/3BM3N3dWT36NQg92ptITg5IqGADuHl4UCX2ah/r+gZuHNxFRCXwyZBS9+7xxW8pvGIZhGLfCVDIy0Ef/e5ety4Zx/uQOLp87QK6QB4i8eJITO/+mTMt+yabLX/ERsgWWJLD0w2QrWJWff53KypUrM7HkhmEYhnHrTCUjAz322GM0aFCP37+uQ2xMBLGRlzi1dzl5yzyMp69/imkL1nwMTx9/jm+cSaHQ7nw1ZFgmldowDMMw0oepZGSwObNn0f2FZ7h6/jB7VowmIS4Gd0/vVNO5e3ojbu5kzVWIgEKVWLZsWSaU1jAMwzDSj6lkZDARYeiQbzh0YA+nd84DhLN7V5GQEJ9iuvA9K8leoCzxcdF4ePkQGRnBvHnzSEhIyJyCG4ZhGMYtMpWMTFKgQAEmjh/Lpql9cXP35MSGmcnGjTx/nJMb/yRP0ZrERF4i8sIJfD3z8vJzAylauCQLFixAVc0UV8MwDOOOZioZmaht2zb8NWsaxYJysPr7rpzdu/qGONcunGTZ5y0o1/QNdvw9lKL1nuPA3O/oVO4j3q2/gTp5etK6dUe8sniTJUsW8gcW5J133+f06dO34Y4MwzAMI3lmnYxM9uCDDxK2YQ3jxo2jx0tNyBZckcJ1n8LNw4vwXcs5vnYaJR/uzpWzh7h8ei9u7p74R3tTtWAbpm/7mIUHR1Gn6UeUeaAzXt7+nD21jTmLv2fkyKosWjiPcuXK3e5bNAzDMAwARFUz9gIi7sB64LiqthCR0UA1QIA9wDOqelVE+gLPA3FAOPCcqh4WkfrA1w5ZlgY6qep0ESkCTAZyARuAp1Q1RkSyAOOBqsA5oKOqHrLLMxDoBsQDr6rqvJTKX61aNb3Vrd6Tc+3aNfr07cvPv04jLkHw9suLl19Ozh/eiH/+kkSFH6FUnrr0rDmWbScXMHHbADq8soKs2fLdkNfODRPZtOxt9u/bhbd36gNLDcMwDCM5IrJBVavdaj6Z0V3SG9jp8L6PqlZS1YrAEeBlOzwMqGaH/w58DqCqi1W1sqpWBhoAkcB8O81nwNeqWhy4gFV5wP73gh3+tR0PESkLdALKAU2B7+1K0G3h4+PDD8OHc+nsSebN/I1Xu7Ulp1s4ef2K8GC2przXYBl9H5qKj2c2/tz9BXVbfXFdBSPyyhnWLfmCBb/34OSx9YhHLn777bfbdTuGYRiGcZ0MrWSISDDQHPgxMUxVL9vnBPAB1A5frKqRdrTVQLCLLB8D5qhqpJ2+AVaFBGAc0MZ+3dp+j32+oR2/NTBZVaNV9SCwD6iRHvd6K9zc3AgNDeWDDz7g2We6UiX/I3So9BEFAkoDcDbiCKev7KNomeaAtQ/Ksr8G8tPnpTh3Zie5C1Qka7ZAoqJj6PVyH8LCwm7n7RiGYRgGkPFjMr4B+gPZHANFZAzQDNgBvO4iXTdgjovwTsBX9utcwEVVTZxicQwIsl8HAUcBVDVORC7Z8YOwKjC4SONYvheBFwEKFSqU4g2mtxo1ajD8q5dRHZy0X8nV6HP4+QXi5m79uZb8+Qanjqzh2Tf34OuXJylt9Qb92b3pVxo0bMKqf5ZRunTpTC27YRiGYTjKsJYMEWkBnFHVDc7nVPVZoABWN0pHp3RdsMZsDHYKDwQqACmOoUgPqjpSVauparU8efKkniAd1atXD8+sCYQd/zMpzD9LHq5cOU58fCwXz+5n54YJtOn253UVDLDW5ChdpRMVHuxHj169M7XchmEYhuEsI7tLHgRaicghrMGZDURkYuJJVY23w9snholII+D/gFaqGu2UXwfgD1WNtd+fA7KLSGJrTDBw3H59HCho5+kBBNjxk8JdpLkjiAg/jvmBsRu7sfHoTFSVnFmDCQoox76tf7B59UjKVX8Wb98cyeZRsdaLrFyxnI8+GgRAREQEY8aMYcCAgbz77nusWLGCjB7waxiGYRgZVslQ1YGqGqyqIVjdHIuAp0SkOCSNyWgF7LLfVwFGYFUwzrjI8glgkkP+CizGGqcB0BWYYb+eab/HPr/Ijj8T6CQiWeyZKSWAtelzx+knNDSU6bN+58+jb/LOgnJM2vQ6AVlysmj6K5w5toHgYqEpps/iE0CeoCp8/sU39OjRiwL5C/L9oOns+zsbYTPj6NTuOSqUrcK2bdsy6Y4MwzCM+1Fmr5MhwDgR8bdfbwZ62ucGA37AFHsswhFVbQUgIiFYLRBLnfIbAEwWkY+wZqeMtsNHAxNEZB9wHquSg6puF5HfsMaCxAEv2S0qd5zQ0FB279vO0qVLWbduHZCfw4eDGTNuMvFxMammT4iPpWSNF/j55x95u8k68voXSzrXpsqHrNo/gfqhjVi5ehklS5bMwDsxDMMw7lcZvk7G3S4j18n4L157rQ/zlp2kWZfJyca5fOEIE76qTJf/28WE94vxxeNHiIg6S0z8NXJkDcYvS04A5m4bTFS+NcyY9XuyeRmGYRj3n/RaJ8Os+HmXeeedtxldpDjnTu8kV74yLuOsWTiIUtW64JstL55Z/Hjvz6rEucXh5R3A1YtHKR/clFbl3qJeye4MnP4JJ06coECBApl8J4ZhGMa9zuxdcpfJlSsX3w4dwm/DQjm85+/rBnBeizzP4umvcfzgcmo2+4CY6KtERV2kSo8faDX0MI8O3krrIYehRi0+W9CQQ+c2UCB7OZ599gX69XuTDRtumAhkGIZhGP+Z6S5JxZ3WXZKoe/ceTJz0B96+uckb/ACx0Vc5tn8JIeVaULf9N/hkzcXWFcPZv28+9fr8cUP6UzuXsOqbDuTxDSF7qZr4+uZgz+bxFCtWmD+mTiYwMPA23JVhGIZxJ0iv7hJTyUjFnVrJuHLlCsVLlCGkSjf8cxXF3dOb4BL18c2WF4Dzp3cxbWg9Ql/7gzwl67jMY/mXbTmxbR4vvH0E36y5SUiIZ93ijzm25xc2blhN9uzZM/OWDMMwjDvE3bR3iZEBsmXLxt8L5nAobDRnDi7BP1cRPDx9uBS+n7Vz3uW3L2pQ5YnByVYwAAqHdsE3exC+WXMD4ObmTs2G7+KXqxrfffd9Zt2KYRiGcY8ylYy7WIUKFdi+bROdWpRm3R9PM+69QGYPf4hAn13kKhBC0bpPp5jeyzcHPnbLh6PKdfoy7PsRZsEuwzAM45aYSsZdLnfu3AwcOICD+3dxLfIqZ8NP8t3Qr7ly7jhxMddSTHvhyGZy5Cp+Q3i+4Ae4dOkCly9fzqhiG4ZhGPcBU8m4BwUFBVGjRk0Or05+LY2EhHj2LRpJxerdbjinCQnEx8fh7u6ekcU0DMMw7nGmknGP+uiDt9k+9U3OHbxxWmpCQjzrxr5EgH8wwUXq3nD+0J75FC9eGj8/v8woqmEYhnGPMotx3aPq1KnD+J9G8lTXxgRVaUmB6k/g6RvAhUMb2fHXFwRkK0i7Z2YmbSefKD4+lo3LB/F//V+6TSU3DMMw7hVmCmsq7tQprGl19uxZRv80hqkzZnPt2jWKFyuKxMexat1uAos04ODe+URePkkWn+yElGxK5IV9FC2Uhb9mT8fDw9RBDcMw7kdmnYxMcrdXMlzZsGED9Ro2JnvFxhRp9RpZA4sTffEUh+eN5MSicUwc+yOPPfZY6hkZhmEY96S7Zp0MEXEXkTAR+dN+P1pENovIFhH5XUT87PC+IrLDDl8oIoUd8igkIvNFZKcdJ8QOLyIia0Rkn4j8KiJedngW+/0++3yIQ14D7fDdItIko+//TnP+/HmaNGtJ6Z4jeaDfJHKUqomXfy6yFSpH+ReGUOPDhTzXvRdbtmy53UU1DMMw7nKZMfCzN7DT4X0fVa2kqhWBI8DLdngYUM0O/x343CHNeGCwqpYBagBn7PDPgK9VtThwAUicKtENuGCHf23HQ0TKYm37Xg5oCnwvIvfVFIofR/9EQIWGFKjT3uX57MWrUrBFXz778ptMLplhGIZxr8nQSoaIBAPNgR8Tw1T1sn1OAB9A7fDFqhppR1sNBNvxygIeqrrAjndVVSPt9A2wKiQA44A29uvW9nvs8w3t+K2ByaoaraoHgX1YlZb7xuhxE8nf6MUU4wQ/0o1pv/9GbGxsJpXKMAzDuBdldEvGN0B/IMExUETGAKeA0sC3LtJ1A+bYr0sCF0Vkmt3tMthufcgFXFTVODveMSDIfh0EHAWwz1+y4yeFu0hzXwg/cxrf/EVTjJMlIA+KG7NmzSI+Pj6TSmYYhmHcazKskiEiLYAzqnrDQg2q+ixQAKsbpaNTui5ANWCwHeQB1AXeAKoDRYFnMqrcdhleFJH1IrI+PDw8Iy+V6QKy5yD6/IkU48RGXCImJoqX3xpEgUJF+G7YcLPEuGEYhnHTkp2jKCLt0pA+SlX/Subcg0ArEWkGeAP+IjJRVbsAqGq8iEzGaukYY1+zEfB/QD1VjbbzOQZsUtUDdpzpQC3gJyC7iHjYrRXBwHE7zXGgIE5oRfQAACAASURBVHBMRDyAAOCcQ3gixzRJVHUkMBKs2SVpeA53jaee6MjkRWPIXqJ6snGOLZ5AgRptqfn6r1w8EMaHX3XlyLFjfP7Jx5lYUsMwDONul1JLxiigBdAyhcNVVwcAqjpQVYNVNQRrsOUi4CkRKQ5JYzJaAbvs91WAEUArVT3jkNU6rMpEHvt9A2CHWj+tFwOJcy27AjPs1zPt99jnF9nxZwKd7NknRYASwNoUnsE9p0f3Fzi9agrnti93eT7i1AH2ThlE8ea9AchetArV31nEyJ/GsXHjxswsqmEYhnGXS2m1pTmq+lxKiUVk4k1eT4BxIuJvv94M9LTPDQb8gCn2KpRHVLWV3eLxBrDQrphswKoAAQwAJovIR1izU0bb4aOBCSKyDziPVclBVbeLyG/ADiAOeElV76tBBwUKFGDqr5No16EdQY17ENykOz65g4m5eoFjC8exb9rnlOnwPrlKW1vEx0Rc5OiyX3D3zkObNk/w+OOt6dHjBUqUKHGb78QwDMO405nFuFJxLy7GBbB3716+/Hoo4yeMJ+paJOLmToGa7SjevDc5S9YE4NjKKWz+oTvFijamXKl2eHh4c/zESjZvGUOnTh0YNmyI2UTNMAzjHpRpK37aMzmaAyE4tHyo6le3evG7wb1ayUh05swZipQoTZPRZxG3f3vPTofNI2zoM3TtNJf8+StdlyYq+jK/T21F00cfYOjQ++JjYBiGcV/JzBU/Z2HN5sgFZHM4jHtArly58PHx4cqxndeF7/r5bRo89D7Hjq9hzbphHDy0JGmGiXcWf9q3m86YMWM5fvyGcbOGYRiGAaStkhGsqu1U9T1V/SDxyPCSGZnC3d2dni++wOG/vk4KO7luFtdOHmTp4nc5v3sV1/ZvY8GfrzLsu9Ls3j0LAB/v7FQo/wSjR4+5XUU3DMO4b128eJEhQ4ZSr04jqlWszWNtOrFgwQISEhJST5yJ0rLN5hwRaayq8zO8NMZt8VrvVxg3sRb7pn9G/lqPs+nbbrQq25/GpXvj6Z4FAFVlx6lF/DDrKWJjr1G+fAfy5a3Ozp2Lb3PpDcMw7i8LFiyg02NPUjF3I0Jzv4ZfQE4OH9pKr6f6kaeQP3/OnU7OnDlvdzGBtLVkrAb+EJFrInJZRK6IyOWMLpiReXLlysXKpYtw2zmdf/rX4tESr9K8XP+kCgaAiFAusCGv15/N7L96EhMTQXTMFXx8vG9jyQ3DMO4vmzZtotNjnXmjwjT6lJ9E9fwtKJOzDk0Ld2dwjY3kuVSFFk1a3zEtGmmpZHwF1AZ8VdVfVbOpqn8Gl8vIZAULFmTWH1MgIZYmpXsnGy8kZxWK56nD1u2T2bvvV1q2bJqJpTQMw7i//e/dQbQr/DblctW94ZybuPFsqa85e/QqCxYsuA2lu1FaKhlHgW1q5rre89atW0fZ4Lr4eKY8rrdqgZbs3PY7EVeP0LJly0wqnWEYxv3t/Pnz/L1wPg2Dn0k2jpu48Uj+nvzw3Y/JxslMaRmTcQBYIiJzgMSlvu+bKaz3k4SEBNwk9XUv3Nw8OHlqDQsW/YWHR1o+QoZhGMatOn78OHmyBePrmXJnQoh/JVYeGJVinMySlm+Ig/bhZR/GPapSpUrsPtWL2Pjo68ZjONt84k96vvwCtWrVysTSGYZh3N98fX25GnURVcVeGduliNiLZM2aNRNLlrxUKxlmuur9o3jx4lSqVIkVB8ZRv8SLLuOEXz3E7rOLmT9wtMvzhmEYRsYoWrQoATn92X5uGeVz10s23orwn2nb/c7oyk52TIaIvJ9a4rTEMe4uXw/9nBk732bd4ak3bO9+8vIevl7ehA/+9z45cuS4TSU0DMO4P4kIr73+Cr8cHEhMfJTLOLsvrGHDmdk8+9wzmVu4ZCS7rLiIHMOaWZJsWuAFVS2dEQW7U9zry4q7sm7dOp7o8BQJUV5UztceTzdvDl5eyZ4zK/n44w956ZVeSXETEhJYtGgRa9asAaB69eo0atQIN7e0jCk2DMMwbkZCQgKdOz7N1pUH6Vz0E8rmrIuIcC3uKkuOTeS3g+8yftIYmjdvfkvXyfC9S0TkvTSkv6qqX95qIe5k92MlA6zFtxYuXMiihYuJiYmlXPkydOzYEV9f36Q4S5Ysoevz3YkSbwIqP4oA5zbMJvLsceKiIhCEChUeoPdrPenYsSNeXmZIj2EYxq1KSEhg+PAf+OaLb7lyMYJs3jk4c/kIoaEP8+7/3qJ69eq3fI1M2yDtli9gbbC2Hjiuqi1EZDRQDaslZA/wjKpeFZG+wPNYW7CHA8+p6mE7j3hgq53lEVVtZYcXASZj7auyAXhKVWNEJAswHqgKnAM6quohO81AoBsQD7yqqvNSKv/9WslIzZIlS2jVvgNlXh5L3qqPJg1CUlXObvqbTV90oVPDUYgI/2z+gtz53Zg7bxZ+fn63ueSGYRj3BlVl7969REREEBQURN68edMt77upktEXq1Lhb1cy/FX1sn3uK+CMqn4qIvWBNaoaKSI9gYdVtaMd76qq3vDtJCK/AdNUdbKI/ABsVtXhItILqKiqPUSkE9BWVTuKSFlgElADKAD8DZRU1fjkym8qGTdSVYqWKkveToPJV72FyzjhYQvYNaQnA5/dQ2xMJGNmNsfb7yLPPvs0zZo1o0yZMplcasMwDCOtMnMX1v9MRIKxtolPWhXEoYIhgA+gdvhiVY20o60GglPJW4AGwO920Digjf26tf0e+3xDO35rYLKqRqvqQWAfVoXDuAlLliwhIt6dvNWS7/PLXbkR4peNP/7uxSc/FsHN3Yec/o8y7bcDPFinAaGhjTh69GgmltowDMPIbBm9ktI3QH+ctoYXkTFAM2AH8LqLdN2AOQ7vvUVkPVZXyqeqOh2ri+SiqsbZcY4BQfbrIKyVSlHVOBG5ZMcPwqrA4CKNY/leBF4EKFSoUFrv9b6xfv16sldsnOI8bRHBzS+AfceW8sIz/5A7Z4mkc40bfM2qdV9Ru3Yo69evIn/+/JlRbMMwDCOTpVrJEJGhLoIvAetVdUYK6VpgdYVsEJGHHc+p6rP2WI1vgY7AGId0XbC6VxwnARdW1eMiUhRYJCJb7TJkCFUdCYwEq7sko65zN0u+emGJPH2IK/vD6NNjN9n8rq9EeLh7UbfWm0RGhvPOOx8watTwjCuoYRiGcdukpbvEG6gM7LWPilhdGd1E5JsU0j0ItBKRQ1iDMxuIyMTEk/Y4iMlA+8QwEWkE/B/QSlUdlzA/bv97AFgCVMEa0JldRBIrSsHAcfv1caCgnacHEGDHTwp3kcZIoxo1anBh09wb1tFwdGTOCCqX73JDBcNR7eqv8+vkyVy5ciUjimkYhmHcZmmpZFQE6qvqt6r6LdAIKA20BRonl0hVB6pqsKqGAJ2ARcBTIlIcksZUtAJ22e+rACOwKhhnEvMRkRz2bBFEJDdW5WWHvWHbYuAxO2pXILFlZab9Hvv8Ijv+TKCTiGSxZ6aUANam4RkYDkJDQ8mWRTi9dlaycc6HLaRMiVYp5uOfrQC5cxdl165d6V1EwzAM4w6QljEZOQA//u2eyArkVNV4EYlOPplLAowTEX/79Wagp31usH2dKXZff+JU1TLACBFJwKoUfaqqO+w0A4DJIvIREAYkrnU9GpggIvuA81iVHFR1uz0jZQfW+I6XUppZYrgmIowZOZwWbdoDo8lXo+V1U1jPbJxL5KkDacorIcH0RhmGYdyrUp3CKiLdgLexuikECAUGYU0FfV9V+2VwGW+r9JzCevjwYWbPnk1ERASFCxemVatWeHt7ExMTw+7du4mLi7PWpg8ISDGfixcvsmHDBuLi4ihXrhzBwSlOxMkwy5Yt45kXenA1VshRuQlxcQmcWjsL9+g4iuarTVb//DRrnHyP2uUrJ/juxzKcOnWMbNlS3l7eMAzDyDyZuk6GiATy71TPdap64lYvfLdIj0pGeHg4Tz/fneXLl5H7wdaQNTuxh7Zydd8mHqhShU1hm8iSNSfuHl5cCj9Cu3bt+d/7bxMSEnJDPm/0G8i0aVMJDKpEXGw0J0/swMvLG79sflSuVIlXX+nOI4884nJZ77i4OGbPns2uXbvw8vKiUaNGVKhQAYAdO3bw/cgf2b5rDz7e3rR+tDFdunROdSc/VWXZsmXMnTuXr74awoOVXmHjzl+4EnEKdw8v+vban+y4jNnze5MjzyH+/HM6ixYt4qtvhrFsySLi4mIpXrIsvV/pTufOnfHx8flvD94wDMP4TzK7khEEFMahe0VVl93qxe8Gt1rJuHDhAlVrP4hWbUlQ1/dw9/53We7Io3vY8VZrgovXp85T3wNw7fIZ9iwaxpFVo1i+ZGHSolWnT5+mVq265C/Ukup1B7Bp3Sg2/PMtFWu9SIkK7XD39Ob4gWXsWPsdZUsXYsb0Kdd9OY8fP5E3+r9J1hwh5CxYi/i4axzeMp0SxYuSL39elq78h7yPPo9vyerER13l8vLfuLxtJb/+PIGmTZumep/vvfsBS+eeo2XoUJas+4IzJ7eQJ6AYGw5PpdNjU6+bwhoXH8OKVYNZ+s8g2rdvR2xsDMv/2Uj5um9QtFI7PDy8OXFwJTtXfkMWwlm8aB65c+f+z38DwzAM4+ZkWiVDRD7Dmma6HUiwgzVxae973a1WMt7o/ya/7QynSD/XW6PHXAxnY9fyNO+/jIDAUknh+1aM4fTKr9i9YwsiQtt2HQm/VIy6jwxi8/oxrFn+OY/3WISff+B1+cXHx/L3lKcpXcSNKb/9DMCIEaP4v/cG8fAzU8hbuNp1cbcvHcba2e9Redgq/ELKXpfXpW3/sPOdlsyZMY169ZLfVnjbtm20b9WR48dOEhN/jRy+Bbh47RSvtlzIoTOrmRP2MUH5q5E/sDKRkWfZtnMK2bIVoPoDPYiOvsia9d+TPV9pHn1+Kj5ZcyXlq6qs/rM/WWLCWL7077Q/dMMwDOOWZGYlYzfWEt03O8jznnArlYzo6GjyBRWk5NCV+ASXSDbewRFv4Xf2GjWf+DopTFWZ90EFfhnzLaVKlaJU6fK80PcQXl5+/PBlcZp3nkRg4Zou84uLjWLMZ4XYsO4fcufOTcFCRWj5+lqy53Vdhg1zB3H4fBhl/zflhnMn547l0A+vs3TBHGrUuHFx1DE/jaV/7wG08+tB6xzPksMjDwejd/HL2W9ZePF3Oj08gipF27F2z0Smrx1IgaAaNG34Bfny/FuhSUiIZ+6i/uw/voz2ry3Hw9P7unOTPinGvL+mUrVq1WSfoWEYhpF+MnNZ8QOA561e6H504MABPLLlSLGCAZCzdjPOHFxzXZiIUKDak/wxYxaLFy+maImGZPH25+ih5WTxDiB/oeRXQ/fw9KZMlaf5acw4xo4dR+FyjyZbwQCoUO9lzof9TfS5kzecy9ugE3FxCTSq2Yig7EHULlebET+M4OrVq6xcuZIBvQcyKngZPfN/QIEsIfi4Z6Wsb1U+KjSW9wr+yOQlL/LljLqs2jOekJD6dO3013UVDAA3N3cebfgFWb2ys2fj5BvOlaz2HOPG/5ziM8xIkZGRbNmyha1btxIVFXXbymEYhnG3SUslIxLYJCIjRGRo4pHRBbsXqCpxcSnPkI2PvkbUqYPERUcQF3v9F5iXb3auXo0kJiYGD09rAOaVS8fIkbd0ikt6AwTkLsORI8dZtnIN+Uo9mmJcLx9/8hSpyZW9G2845+7lTUBQGQYxiJWXVvLujncZ32c8hQML8/QTXWnu+zSFs5R0mW/THJ2omrUeV6PPc/ryHurW7p9sGUSEOtV6s27O/4iNjrjunF+OEA4fPsrHH39CSKESeHp6EZAtB12efIaM3Lzu9OnT9HylN/mCCtKo/ZM0aNuJfEEF6d33Dc6dO5du14mLi2PGjBkMGjSIwYMHs27dunTL2zAM43ZKyzoZM+3DuEkxMTFcCT9B1MlDeAeGXH/uYjjHJn7K6Xnj8PDOhqgy5fWCFH/wGSo0fxNvv1xcOb6ZEo8UpWTJkpw6NhhVxcvbn2tXw1O99rWIM2TP50/kySgk1UXAQST5+qZnVDyVqMQudtHLozd4+VKxUGdiNZq/Dv3OP/sX8n6BkZTxfeCGtE/m6c37p3sQE3OV7AEhN2YOnDqzlY2bfuL8ub1EXz3HpE8r0+GNNXhnzYkmJHD57D7+WfE3F/dlo8sDkwjKWYGI6POs2T+Rpo+04u133+S1Pq8mW/5jx44x4odRzJ+9kJiYGEqVKUHZiqXw9PTE19eX8uXLEx8fj5eXF9WrVydr1qwcO3aMmg+F4lm9JeVGhOGdz9rD5tqJA/zx6+fMqP0ga1csu+WtlSf9MpnX+ryBv38IgfkfIiEhmsGDO1CgQB4mTBidNPvHMAzjbpThW73f7W5lTEb3ni/zx6IwPMtWoXif75LCo8+dZMvLoeSv8AglW/cja74iAFw9tZ+90wdzdutiGr46gwWf1Gbv7h3kz5+fEiXLUy30G4IK1+H7zwrxVJ8w/HO43rxNVfnlm7L8Nmkka9euY+y0TdTtMj7ZcsZGXWX8O4Wo9tNWsuS5fr+4a8f3s6NbLcbEjOB5z1482WwipYIbJrWkJGgCG/ZMYsbSPvwQMo/SvlWuS38oajc9jregSIkmREZfoEObf7s9YmOvMW3mMxw5uoLQEs9TKEclImMusnzPaI5d2k6B4qEc27uI+PgYBDdqFHuChuV7UyhX5aQ8zl09wlfz6zLu5xEuZ8GMHPkj/fv2J7TAk9TO2Z4s7r5sv7CcGQeH4J3Fn8J5q7Hx0HQ8/bLhE5iPyBOHeOqpp9gQtolzZZoQ3Hmgy2d2dPT/Uez8Nub/mez2PS5dvnyZ8eMn8NNPEzl4cB+xscKTnf4kOOjf7q+EhHg2bxnPkmX96devD+7u7uTKlYs2bdqkOstm69athIWF4ebmRu3atSlWrNhNlc8wDAMyYeCniPymqh3szchuiKSqFW/14neDW6lkBOTIzUN95rH0+w7kadGN4Cf74ebhybY3mpK3aG3KdnjPZbodk97h6KKxPN2pHd8PGwLArFmz6PpMD9o+NZdtG8dz/sJ+Wj71G27uNzZGbVr5Laf2jmHrlg1cuHCBwiHFaN1vI/65i7i8XtiCwRw4vYpyH027LlxV2fdeZ9quCuRXJvNk618pVuAhl3ms2TmOsLUjmFjsn+vCV1/5m8+vvE3nJ2bz+fAivP7SAbL65kZVmfx7e3xiPXmh7ng83bNcl27z0dmMXNqZ3g1nUip/KJejwlmxdwzztn9Jp9pDqV6s47/X3j+J3dE/smzlwuvy+OOP6XR/5mUa5H2OPVfWEx0fSaBPUZoEd6O4f1UGb+mC+njzdL3RzNk2mEWHfqLEoCmcHP8h5zcvp9a0U7h5uB6OFB8VSVinQmzduI4iRVw/V2fbtm2j0SOPkit3ZQoVbcHieW/y/LMryJu33PV5x8cyf+GbbAwbRXDxugTkKsLVCwc4dmAFHR/vxPffD8HHx4eoqCiWL1/O5cuXuXjxIiOGj+HgwcOUDKpHfEIcu44spFq1anz/wxBKlnTdpWUYhuFKZlQyAlX1pIgUdnVeVQ/f6sXvBv+1kqGquLu78+SYWK5dPMmKUV25dHInAdUe4cKq2TQbdRx3zywu08bHRDH3+fxs3rCGUqX+ndb6yy+T6N6jF4WLNSD8zEG8vP2p0+QDgouGIiJcCN/L5n+GcPLAnyxbupDjx4/z/bBRrFy5gouRcTR7aQ65CpRPyi8hIZ6d/4xm5ZTXKNz1HYIf64O7lzWz49rx/Rz74W3ybDjMG1Ev81HuH+jVKfmlUeIT4vjfmBC+LzSbUj6VksL7HX0C78q1eLB6b0b+XI/qD/SgYvknOHJsFdP+eIpBbXbg4e7lMs85Wwdz9Pwmuof+2/px7MI2Bs9rSN9mCwnKad1LbHw0A34LZM++HUnbxqsqIQWLcfbMBUoXakzlMk/g4xXAsTNhrNoynGCfEvStMI6Xlpfnufo/c+jsOjYcmcbFnO74Va4L0dco0fs7l+VKdOiLbvRr+gAvvfRSivHAWi+lSNHieHsU4PLlo2TxyEpE9AXKlmpL7Yf6EZjfap1J0ASmTO3INY2k0sN92P7PSI7sno9fQBBRkReABALz56Bl61aM/ukn/INKkYAbZ3dt5Il6I3igREfc3ayKZ0zcNVbuGMHSHZ8xYeIYZv81j2l/zCAi4gpBQYXo2f05unZ92qy2ahjGDdKrkpHsmAxVPWn/e19UJtKbiBCQIzcR4YfIlq8YTQYs5OKx7WyY3J/gWo8lW8EAa7Bl4Yc6Mnfu3OsqGU8++QTNmzdj/PgJLFjgxoEDe1n4+5NERV3Dy9MbNI5u3Z6j59hF9OrxKtu37KN2iV50rNmLDfunMGPwgwQElqZgmUdIiLvG7jU/ox6e5K7UiMO/fMrxX7/Cv3AF4q5d4dLRbbyU0ItP40bzrrxP8eIpDx51d/OgTMEmbItYk1TJWHJpJusil9C7grWVu5dXNhYt/4C8ecqxfuNI6pfskWwFAyC05PP0n1KEq9Hn8cuSE4DgHOVpUOZlFm3/lhrFn+TgGWtWjo9XNsLDw5MqGT///DOnz4TTo808ihSok5RniYL1Ca3yKr8ueJ7BmzuTzScfw5a0p3iNJwms2xHfM3vZP30k3oFFiI24TOT+LVzZbQ3EzFamJv7laid1FcV7+xMRcf0g1eS0b9cBr7gAnij/AQ8EtcTDzZOrMRdYdmAsEyc0pnWbsZQs0Yw9e/4k/MJearX4iHnjn6B2w3do1u5Hsnj7o6ocObCERbNeY8SPE6j90Ur8AouzuEdJnm08mQohLa67ppeHD/UrvsbRM2G0aduBSg+/zMNP/4l31lycO7mNYWO/59PPv2Tp4gWmW8UwjAyR6sBPEWkHfAbkxdq7RLAW4/LP4LLd9Tp37syK5T9S8bFPAMgeXI48xWsSJbGppnXzzeHyCywgIIBXXnmZV155GbA3JDtzhri4OPLmzYunpyctm7fl0vHs9G+zFXd3q7m/aP7atKn1CUu2DGXeov/xROcO7F4bS7wKmjc3Bar0I3zSUCr4P0Sx8g1Yce1jqp2ohC++Vl9ZKrNZsOMoCRyI2slv539g7uVf6dxhJj7e2UlIiOd0+DbKlmzLuMmPIgkJNK73XIrZZc2Sg9x+IZy7ejipkgEQWvw5BkwrzoETy6ma51FEQaLcaNKgGXVCaxEYGMTMmbNpVfeL6yoYidzdPGjfYBjvjSlIwYot6NJlDR5e/66OWqfDNyyb8ALrHw/B2y83xUo+iqIcmDIcsvpS9I3hBJSrzeUdq/Fv2PWG/J1NnTqVrev2MqjJxuvuw88rB81K96FE7tp8OrUp5co/wcEjS6j48Kss+PkZ2nWdQXDIgw6PVyhcrD5dXlrNxGG1OLt1MZGnD5DVPYDyhZu7vPa+E8vZfnw+7V9fQ878ZZLCff3zUbBUQ7Ys/ZbKVWry1+w/qFu3bqr3YhiGcTPSMrvkc6Clqu78LxcQEXdgPXBcVVuIyGigGlZlZQ/wjKpeFZG+wPNYu6OGA885tqLYO7fuAKar6st2WFVgLOAD/AX0VlUVkZzAr0AIcAjooKoX7O3lhwDNsKbmPqOqN87bTCd9er/MhOq1yF+pJXlLWF92WXMV5kzYtFRSQtSxLRRu3znVeCJCvnz5kt5v3LiRNas38M7j+5MqGIk83bPwSJV+5MxWmB1bvyMmLo6CPT4ksF0vALKVqUHYjx+yetuP+Prkph/9aUxjamhV5u8fBVXfTLYcCQnxbD4yi5XXRpPdJ5AqFZ+hZ9X1BPhbm7ft3DsTv6z5aPHI1zza4HNGjqtDdFzqrQDRcVfxcBqvEeAbSFx8NMMe2oq7mzsA3cp8wfrwOXwxqwtlirbg1KmTVGvZJdl8w/b8in/+0jz8zFjEaZ8XDy8f6j87numf1aJM6fbUqj8AAE1IYM+2acx9sxVFXhtKxKEdnD6T+kyfLz4dQufKX1xXwXBUInctqhZoxZUTh7hy+RhRkecoUrLJdRUMR56ePjRo+Q3Tx7Uja1ApKuapy5HwDVyKOIG3lz9F89dJaiGav+lzarT88LoKhqOK9V7hQNhUHm3SghEjh9O5y5Op3o9hGEZapWWdjNP/tYJh6w04pu+jqpXsgaNHgJft8DCgmh3+O1blxtGHgPOggOHAC0AJ+0icXvAmsFBVSwAL7fcAjzrEfdFOn2GKFSvGb5Mmsuq71mz+9Q0uHd9JYMWmnN25nMgzyfdCRZw6wPl9a2nbtu1NX3PkiJ+oVfLFGyoYjioXacvOnbvIkjcoqYIBkL1mY8qMWk6pkUvJ+9bXaK2aVJLKXOEKJ89t48jp5MemhO2bgmoC//fqGQb0PkHj+oOSKhhHT65j+vyeNH54EADu7p6ULd2etYdvXGHU0eGzG4lPiCXQv9R14eFXDuCfJU9SBQOsylb1vM3oX/kX9h9dTP6cZcnimfzmbos2f80Dzd++oYKRlJ+bGw80f5eNq4ddF1aq4mM82mYEB758iVxVHuHIsZT3Cjx37hxbt22ianDrFOM1KPYCl6PDcRcPDu+aT9kqKVcwCxV9GC8PH3JFeLNy7xhGLOzAgkOj+G3DAN6eWJi/1n/IpcjT7D+xjJLVUq44lK/3MnlzVuSVl15j9+7dKcY1DMO4Gcm2ZNjdJADrReRXYDqQtLS4qqb6c1xEgoHmwMdAXzvdZfucYLVAqB2+2CHpaqCLQz5VgXzAXKxWkMSdYf1VdbX9fjzQBpgDtAYetpOPw9qmfoAdPl6t0a6rRSR74gDX1O7lv2rSpAmbN67j2++GM3FIYy6cD0fl/9k7z7Coji4Av7P03sQCCKhgF0XsvceKvUWNMXaTWGJsMWqKGrvGFnuN3Rh77x1FUSygYgFpAtJ72fl+sCKEovmixuh9n2cf7p45M3fOXcGzM2fO0eHCtNY0nHoOXZOc325TYiO4saAraruarAAAIABJREFU48eNw9DQMJ9R8+fxQ39KWLQoUEdLSwdLk5Ik18z7m7JBcWcMijtj5taEK62s+F49ibikKJbsbcGQdgcpUbRWlq6UkluP97Ll5ED0tAyZtdSRKuV7UsqxGenpyXjd+Z3gZ150brMGJ8emWf3cKn/BgstlCSw3EjvL3Lkg1OoM/vSaRKMyQ1BlcyYAzt5bQWPbPnnOvar1JxhoGZGQHFngMwh/7oNN2SYF6tiWbUpCbDApybHo6b/cHTQxK46QkuTghzw1iichISHfarVxcXEYG5ijrSo4aa6JXiGS0+MpbOpMbHw4BgZ5r3q8QAiBjq4xoYl+tBiyC5uyTbJiRaKC73J56wgenbqIoVFhdHQL/ndkVqgUSamx1CozkEW/LmXx0l8L1FdQUFB4XQpayWineZmSubXQIpusbQH9srMAGMvLwmoACCHWAqFAWWBRHv36k+ksIDKzRM0Fvv2Lji0QmO19oEYGUCSb4xBKpoPyos/TfPq8NRwdHZk7ZybPQp6SmpJMalIs7o2rc3JEGXw2TyTC5zwRd8/hu2UiZ0dXoE+HFowf+1dzXw8jI0MSU6JeqZeYHIVuYbsCdYQQ6FnbsWLzciJjI1m4ZA6bznZl8f66HPCYxJ6L45mz0YXjx8cxRftnFmotYbXuRho9sCPm+B8cOT6G0FAvenbYRvnSOb/JGxsVpt0nS5h9tDmeT3aSoU7PaguNuc+ikx1Jz0ilZYXROfr5hpzmvN862jrkfaJDCEEd6w4kp8bw9Nm1gqxDZqQX0A5qdTpCqIiNDgAgMuIB2+bV4viKbgxO68+Ip53QuqCNfWF75syYS14ntaytrYlPjiY+pWCnJyjWB0sDWzIyUklJiiE89FaB+ikpccTGBdJ8xEFsy73MWwJgYVOeFsMPEE0McXHBZGQUHAOUEBuCoZ4FNZz6snPnHwXqKigoKPwd8nUypJT9pJT9gFUvrrPJ8i4pmg0hRFsgTEqZ6y+9ZgwbMrdRuv+lX28yVytma0TDgINSykD+DzSrFn8r45gQYpAQwlMI4Rke/uo997+LEILfN6zD89JZGhaJI3r3WKL3jKNR0QSuXDjD7Fm/vDJteH506tqOGwEF1/kIjrxDXHIoKu38T3ZAZgyCVkIkpUuXxsjIiP79vyAw6DEz5o+hYsNkznv/yuyMGVwz8GGo/kja63ahnW4nvtP7gc26O/hWaxwmWhbsPjiQ51EPc41fulQrdPSM2XDlK0ZuLcLCsy2YcbwGs47VwS/iDEXMHXka5U1ccgT+z73Y6DGEhSfbM67KVooavsxNoZZq0jJe1u8TQkXxwm7svzgh3/9gDQwt8b+1v0D7A7z3oaNvipaWLjGRT/hjcSNGRH7GfZUfM3TmMEFnErvU+zmXcZUN0zbw0+Sfc41hZGREu7btOfOk4F+Z4w+WUt2uExGJARha2nL1/DykWp2v/p1rGzCzdqKQXd7parS0dXF1n4K2riGPbu4u8N73Lq6leolumBgUJiExrkBdBQUFhb/D68Rk5LXSkJfsr9QF3IUQT4CtQBMhxO8vGqWUGRp55xcyIUQzYCLgnq3qa23gK804c4DPhBAzgCAg+1dxO40M4JlmO+XFtkqYRh4EFM+nTxZSyhVSympSymrW1tavYer/R7ly5fhtyUJuXr3IzasXWbJoAWXLlv1HY3bp0oVn0be57X8wz/YMdTr7PcfSvkM74g+vLXCsaI8j2BYtkuMYrba2Nh06dGD27Nk0b9iCYHVQvg7RMJ2RJGXEk5aezJLVruw+MBBfvwM8eHSUQyfHMmupI8b2FVDpShYsnskvi79h5e9zCA0P4sFDXxp0Lc5G7158f6A06250x7ZaMtYmRXGxagyAZ9ghJnq1ocNhPTofNab3meJsevAjHuH7qOfyFdpa+qzY2ybHikZKWgLnvX8jOTEar4NTyUjLu7hweloyN47MQiW0MLcqxeWDkxmSMYgh2sPQEjm3b0qoSrIn/Qi/zllAYGBuX/jbcSPZdftn7odfzNUGsOfuDB5HXcfBojIm5na0GneG1LQEju8dnqejEex/mbOHJlC5Vd7ZSF9gV645aanxXNj9LYlxYXnqPL69n9CH56lRujfPonwpWvitL+wpKCh8RBQUk1EbqANYa05+vMAU0Mq710uklBOACZqxGpG53dFHCOEkpfTTxGS4A74aHVdgOdBSShmWbZysCDghxOdkBoeO17yPFULUAjyAz3jp/OwF+gIzND/3ZJN/JYTYCtQEYt5mPMa/gb6+Prv3/kGb1u1pGPUtdcoOwEjfEikl/mFXOeT1PTYldVi9eiVVatQidPNsin46Jtc4qc9DCVk4nOWzp+XrRHw/fSJtG7elUnoV3LRzVoVVSzWj1CNIrVyOpHFfozd6Ip631uJ9bwcqlQ46JlZYFncmLugyG9etxN3dPUf/YsWK8cOPk/nhx8lZMiklZUtV4OqzA9yIOc256H2U6TSBLrV3oqWjT/TTW1w99CvBTx6ToU7lizZ/cNprPkt2NcVI3xIjfSsiYh5iqG+JtW1lIsN8OPJbexr0WYWxxUt/NT4qkLMbB5KaGE21eqNISYrmoc8+hmn75fvci4qidFf1YvmSFfz8y0852nx9fbG0dGDO+fZUtXWnoeNnmOhZExzrw1G/33ieEoy1ZTkeRFwiMT4MLR093H/w4tj81qz9tQpudYdTuJgLyYlR+FzfxCPfg+hqG2HtUHCeHKFSoWNgik297mybVZWarX/E2a0nOrqGxD5/wu1zy7h3aS3DWu7FQNeUiw+W8cXAVx/JVVBQUHhdCjrCqgsYa3SypwSMBbr8n/cTwHrNcVQB3ASGatpma+63Q/OfWoCU0j3PUV4yjJdHWA9pXpDpXGwXQvQH/IFuGvlBMo+v+pEZZ9Lv/7TjvaZ27dpcuHiGn36Yzk/bSlLIrDhJqfHo6gm+Hj6MkaNGoKOjw7ED+6jftDn+Ph6Yd/wK4zJupCfGEnl8K893zOfbr7+kW7du+d6nRo0arN2+lu7dW9NUtqJL+qeYCXNuZXixULWEiHLWJB3ZDUlJpPfuguGyZXzZvz/hz2NQaaloVL8OXbp0QU8v/8Rk2RFC8P2PE/ji866YFClF0+ke6BqZZ7Vb2LtQY/BqzEtVZ+OWvvRqvBa7wlUx0DFlaLNdpKtTeR7nz+ZLX2KUrmLE10NZsmwNO36sQFGn+hhb2hMfGcAzvwvomxbG0NCaGvW/IfipB0465bHCqsD5NUlrzrrzK3LJf13wG81qTcXRpg5Xb61h/a0xpKTGYWZiRzW3L6jk3JkHASc4eWYKhQztCbp1mOKV29B24iUC7xzB84+JpEdHYG1SAldbdwZ2XsDy83149vgyZkWc851PbPhDpFRTqec0zB1cuLBmOKe3DkVLSxdtLV1qlfmc8R09sDJ15KrfFh48O8bAgXNe67NQUFBQeB0KLJCmyXGxXUrZOV+lD5x/UrvkfSA6OpqAgAD09fVxcnJC9Zcjm3Fxcaxbt56Fy1fy9PFDdPX0admqNaOHf0nNmjVf6x6RkZGsXbOOfVv2kZCQgNBR8Sg6iIS0NLSNjMiIjKRr9+5MGTuWkiVL/iN7Jk2awpz5S6k7cgdFyjfKV+/wlJrEBNxGD30GN99BUfOynPVdwWmfpTiUbYaxzlMuXzrLgQMHGDthMk+fPsXAtDBSnUZSbDgmpqZEhEVgZuFAEbvqiFveXNX2KnBuezN283utNew7k7Nosb6+IZMGh6Gna5xv3/SMVCYuMGJYzQ1sfvQzrSaeR88484RJUmwYe76vyLD6Wylnk3kixst/DzvuT6Pdd5dynb55wYVtI0g30sOlzyySo59xe/N4As9to7JTW+qUHoylsT0RcY/x8FvB0ygPjhw9gIvLR1GSSEFB4RW89dol2W50SUpZ+5/e6L/Kf93J+LeQUvL06VNSUlKwsbHJ94jn38WlclUe+AfTcUlIgcGxfqdWcefYfBKD/NDTMkBKiW3J+sTHBmFqlMrFC6fJHm/j7e2Nv78/RkZG1KlTB319fTIyMti/fz8XL15i6bwleGs/oKgomu89h2kNpNQYR76fPDGH3NDAmO8GBqOvl3+S3LS0JL5faEp95x5cCdiHtpEpVTpMpUT1bmjr6PPw8hY8NgylUenBNCn3JRaGtvxypBG6js7U67MiR6E8KSX3zq/k0o5vMLYohK6+EfGRwXTp0pURXw3lxPGTbNywlaioSAoXLkK//r357LM+Sg0TBQWFLN567ZJs3BBC7AV2AFkpGl8nT4bCx4sQAnv7vEvR/79s2bqNB/6BGBVxfuXpG0NLO7SNLdA1NEdbrY1anU56QgwmuiaM/uZT/hrQ6+LikutbvJaWFu3bt6d9+/Y8D3rOwj/mMj3r0FNO/NVP2C134jMod966atVq4fPoAK7leuY7X59HByhRwpneYxqypsmP+Pj4MHveIjZvGIiOjh46Ojp069EZ1OH8vMuF1NQUMtQZFBKB7JjkSOm6gzAtXIakuGcEeK5HVyTidc0DLS0t0tLScHR0zHIiXF1d+XbM6HznoqCgoPCmeB0nQx94DmTPXCQBxclQeGdIKZk4+UfKfjaLe5snI9XqfLN1AiSEP0anUDGS/X357stQhBBIKdmytz2BgYEEBgZiZ1dwnpDsTJ3zM3WO18E0xoyRfIu+0M9qu6H2ord2V6b89ENWgbbsjBw1jNEjp1LJuRPa2rnjT9IzUrnoPZtpM6fQvXvmie6SJUvSpk0bUlJSSEpKwtTUNGura9361SQmJqKvr4+WlhZeXl6sXrOegKd/YGNuyveLpvLJJ5/k2hpTUFBQeNe8crvkY0fZLnk/8PDwoF2Pz6k5/w7nR1enSuep2FTJuzKslJKDk6pi3rInQZtmUda+BSWK1ubStUUkp8RgbWXH82h/XF2r8uNPE2nYsOFrzSEoKIgBPQdyzdOT1lrumKhNuK7jib/qMT/N+pkvBuQdR6xWq+nUqTuPfBPp0HgVpsbFstpiE0LZe2owNo6wb/8utLReeXBLQUFB4a3zzrZLNKnBF5GZ9wLgHJmFyP6v5FgKCv8PAQEBmBYvj0qlwqnLeDzXf01zh3MYWBTLpXt791QytARxvp5Y1m/Pg3P7CAu5QbeWa3GwrYMQgvT0FG7e207HDt1ZsmQBPT/t8co52NracujsQfz8/Dh69CjJyck0d2pMq1at0NHJP224SqVix47NjBs3kV9XVaCkfQPMjByITQjAz/80fft+ztx5MxUHQ0FB4YPjdQI/jwGbgY0aUW+gl5Sy+Vue23uBspLxfnD48GEGTZiB64+nAXjwx0we71tImU+G41i7JzoGJkQ+vobPwXk897+OTa9v8V81GWsjBzLSU0lMi8ZAxxz3pr9S3ullVvzQ8Nus+qMhd328sbV9+4moYmNj2bdvH2FhYRQqVAh3d3fMzMze+n0VFBQU/g7v8nTJDSlllVfJPlQUJ+P9IDExkaJ29tSc6YlhEUcAoh948vjgEsKuH0Kdmoxh0VIkR4eiZVmYtMCHGGiZYGJRnFKVO6OjY0jQo7M8uLuXymW7073l2qzg0b2nvqRxS2t++vmHf89ABQUFhfeId+lknADWAls0op5APyll0/x7fTgoTsb7wzffjmPX5Qe4jN6ByGNrwf/Yam6v+YZihqWITHjKJz3WUqp8zlp+cTFBbF/WhKrO3fmkXmZmzieBFzjrPYJbt5XPWUFBQQHenJPxOuHnX5CZMTMUCCEz2+cHmSlT4f3ml2k/4aCfgNe01kT6XMyqepoU/pQ768fhs3kSqgyJgZ4Z9dpMz+VgAJiY2dJj6BnOXV9AXMIzAPT1zEhKSnyntigoKCh8DLzSyZBS+ksp3aWU1lLKwlLKDlLKgHcxOQWF7Ojp6XHs4D5GftqSJ8v7cmaALeeGlODEsDIkpMfgOGwWRawqEBx1h/JV++Q7jpFpUZwrduTqrTUABIfd+MeZSBUUFBQUclNQgbTJ+bWRWUE9d11rBYW3jK6uLqO/GcWokSMIDAwkNTWV0eO+40aGRCKJiwvG2qYK2jr6BY7j4NycIO9DSCnxvLuEOQvGvyMLFBQUFD4eClrJSMjjBdAfGPeW56XwHvE+5lJRqVTY29vj5OTE0oXzSbp0iCczh2BiVIT09ORX9s/ISEUltDh8fgymFpI2bdpktQUGBjJx8hScK7ti61SaWk2asXnzZlJS8i4Lr6CgoKCQN/k6GVLKuS9ewAoyK532A7YCytryB87Tp08ZPXos5maWqFQq9PUN6dqlJzdv3sxT/+7du0ycOIkBA4by/feT8fHJnV77bZCUlESH9l1xsmnGpEGBDOh4gIjQWyQlPC+w372b2/EPOklCxmWOHN2Ptnbmot6OHTsp61KZNQ+fkzF2KYYL9hHUbggjF62iolt1goOD34VZCgoKCh8EBcZkCCEshRBTAW8yt1aqSinHSSnDXvcGQggtIYSXEGK/5v1qIcRNIYS3EGKnEMJYI/9GCHFXIz8hhHDQyB2EENeFEDeEEHeEEEOyje0mhLglhPATQiwUmjOJmnkfE0I80Py00MiFRs9Pc5+qf/N5fRRcvHiRcuUqseS3legXrkDZZl9h69qRPQcOUqtuYxYuXJylGxUVRfNmralevQ4bNhzj9MmH7Nl5mzq1G9Pyk3ZERUW91bkuXrSE1IRCdGqyGkN9C4wNrSlfsh1Xz+Rfsjw08BohARdYv3EZlz3OUahQIQAuXbpEv2FfUmTlCay+W4xOETtiD28n6dRB1PbleO7oQuOWrUlPT3+rNikoKCh8KOR7hFUIMRvoROYqxhIpZfz/dQMhvgGqAaZSyrZCCFMpZaymbR4QJqWcIYRoDHhIKROFEEOBRlLK7kIIXc08UzQOyW2gjpQyWAhxBRgOeAAHgYVSykNCiFlApGbc8YCFlHKcEKI18DXQGqgJ/CqlLLCe+X/5CKtarebo0aMsWLSMGze8UKlU1K1Th47t23DB4wrnL18FoH7tmgwfNoTSpUsTERGBQwlndIyKUHfkDiyKV8oaLy05nls7p/Dw9Er27dlJ/fr1KVeuEsFh4ZiVrYlxzWbIjHRiz+4n7VkgDoVd0TUM5eKlM+jrFxwj8f/a5+jgTIf6m3CwqZUlj00IZdHWOpRx60X1RmPQ08+sfirVah75HuTUnoGsXrmETp065Rivebv23K3WFvMO/QibPpLYA5up4PIpNkVcSU6O5ubdzcRFPOKXqZMZNWrUG7dHQUFB4X3hrefJEEKogRQgncyCaFlNZAZ+5l+3+uUYdsB6YBrwjZSybbY2ASwFnkgpZ/6lnyuwWEpZ9y9yK8ALqKWZ0ykpZVlNW08yHZPBQoh7musQIUQx4LSUsowQYrnmeoumT5Zefjb8V52MlJQUOnftyTXvBzg1HIFt2aZkqNO5uOVLnj3xwKbDECxqZsYhxF09TPjB1QwdOADS01iweDntZvvmmbIb4NKyz9EKu0q9erXZtGsf5ecfxdipcg6d6Gsn8ZvUgyJm9kz6YSj9+/d/4zY+e/aM0s4VmDQwPFdV1tj4EP48PZwH/icoXrIRUkhCAjzQUmUwauRQpkyZkqNPREQE9k7OOB4L5NnUrzG4/ZB27quwtCyVVYhNSsmdW1s5sncw3t7XcHZ2fuM2KSgoKLwPvPXaJVLKN1HCcQEwFjDJLhRCrCVzNeEukFfN6f7AoWz6xYEDgBMwRrOKUQ3IXj8lEHiRF7pINschFCiiubYFnubRJ4eTIYQYBAwC3ni58nfF0C9HcD9ETatxnmjpZFb+9Do0g7jUSKpvfYS2sTlSnYFKRxezKg0o3HUUa8d+QlKQP04N++frYABU7PA9+ydU4nGAPxWXns/lYACYuzWhxISVBM/7lgXzl762kxEdHc26devZunkX8fFx2NraMnhoP9zd3bPiJl6gVqtRCRVCCKLjAomOe4qutiFFC1XE1LgYfdvuICYuCC/frZy8NpuWA3YQG+nPstXTCI+IZvGi+QghePLkCT/98DMkSfzqFkFbapGAmg1La6Kvb06VmkNwq/0lurpGVHTpSfTzh3z99ShsitlwYP8BklOScChegsHD+tO372cYGxv/jU9KQUFB4cPlrdWCFkK0JXMr5Npf26SU/QAbwAfo/pd+vcncXpmdTf+plNKFTCejrxCiCK+JzFyq+VvHI6SUK6SU1aSU1aytrf9O1/eCkJAQtm/fRu0+67McjLSUBLyPzMK6VV/uTejEhaZ6XGxhxPXuzgRumYNKRw/HKTtIU0ts3NoVOL5JUSe0dQ0xKFUJY+f8s8tb1m1Hho7Ax9c7V1t4eDi7d+9m+/bteHtntp8+fZpSJUqzZfVlXItNoKXLCgqruvPd6HlUdqmWI+hSrVZz7do1klOSWLqlAfPXV+Hg8VFs2tuNGaucOH11LhnqdMxMbNHW1sPOuRE2TvUpW6M37sMvs/vAGVauWs3BgwepWrk6YVdM+KnxeZa1DmBSoxPUK94TXaFLH8fJJNy+zJaVjUlOjgFApdLmzKlzRN4tzvBm55jS8QFNS8xi3cJjVK7kRkBAAM+fP+fevXuEh4f/7c9PQUFB4UPhlVVY/wF1AXdNHIQ+YCqE+F1K2RtASpkhhNhK5krHWgAhRDNgItBQSpnrvKBmBeM2UB+4ANhla7YDgjTXz4QQxbJtl7wIVA0CiufT54Nhy5YtlKjaGV3Dl4W3Hl/dgY6RGbF/rqN+wwk4fXoMlZYOwQGXuXJhPrf2r6XCwhOo9AwQ4tW+p0SNZb32BeoIlQrTqo1JPvIkSxYWFsaIr0Zz4OB+ytjURUfLgIfPPChUxJKngU/5rPkOShdvkqVvX6Qa1cr04ti1qTRt0pKb3p4EBwfTspU7z58nIKQWTR36UatBD3S1DQB4FHGVrTfG4x98kR4t13PuxiIa9F6RNaaegRk13Ofz89T+xEfGMLrmfkpbvYzpMNazpJRlNa4Ft2PF1YEsqneTTQ9/5OjuYZSr0pMbl5cxubM3ViYOWX3K2TalnG1TjnnPxqWCG2kZKZibFiE2LhzXqtX4buK3tGzZssDnFRoaym/LV7Bq/QYiQkMwNjOnW5cujPxqGGXKlHnlZ/JPUKvVeHp6Eh4eTqFChahevToq1Vv7DqKgoPCR8Nb+ikgpJ0gp7aSUjkAP4CTQRwjhBFkxGe6Ar+a9K7AccM9+ekUIYSeEMNBcWwD1gHua7ZBYIUQtzVifAXs03fYCfTXXff8i/0xzyqQWEFNQPMZ/lcCgEAyscsYL+F35HXMjO/oMuUAZly5oaesihMDWoTYdP91OxbJduT+pG9oWhQnxPpzVLyboLsE3DxN2/wIZ6akAJD4PRJ2ahMx49SmLjLRkSpbMnEtYWBi1atQj5n5RprV/yLB6+xlYewfT3B8hYwpRv+KoHA7GC4QQNHf7HnWyOevWraNhw2YUse9GYlw0E5qeoEGpflkOBkDJQtUZ2/gQqXHP+XVzDawc3bBxapBjTFunhkQ8i6SJ/bAcDkZ23Gza4WrTlqOBq/mizEwePjjE+ZM/0a3WvBwORnaau4zBTN+Bro1WMbbHA6Z8HoqDUT/69h7C7Nnz8n1Onp6elHOpwhrvICx+2E6lA6HYLT7DgUQTqtWpx9Zt21/5rP8fpJQsX7ESe6cytO75OUN+WULb3v0pXqo0S35b9l7mSFFQUPjv8DZXMvJCAOuFEKaa65vAUE3bbMAY2KEJyAuQUroD5YC5Qgip6TNHSnlL02cYsI7MHB6HeBnHMQPYLoToD/iTWXsFMk+gtAb8gEQ+0BosFuZmpMY/y3qfkZbC8yfX6P2VR76ZMOs2nYT3zDWkx0fzMHQlFg5V8Dk4j6SoYMyLlSMl4TmJ0cE4NxtGUkwolV0q43t8M/Z9v893Hur0NJ5f3MfE6T8C8M3IsZQ1d6dr1dk59NLVqTx+fpXebbbkNQyQ6WjULT+ckcP741yxNSotParYtqWElVue+tpaunSrMp355zrSqffaXIGhAOlpibQoNTSP3i9pVnIQv3n0pbvTd9Qo3JZLoX9S2cG9wD4Nyw/jjNcCLvqu5mHAKaQ6AyvL0vz4w8+4uVWhSZOcjtTz589p0aYd1mOWYdGwQ5Zcy8iEYoOnYtakGwO/ak5pZyeqVn2zp65HjRnH7/uPYjNuAyYVaiGEQEpJ/N0rTJk/lLu+91i8YF7W85NS8vDhQxISErC1tc06/qugoKCQF+/EyZBSngZOa97WzUenWT7yY4BLPm2eQMU85M+BXFViNfEZX77OnP/LdOnSmTkLGlO14wy0tHUJvncKyyLlsbQunad+XEwQAQ9PoZehooPjaG5HX8Bz3VfU77ua4i5tUKkyK55GBd3h4tbhRD29zoAv+nJ73e9EXT2GRfXmeY777MAaUKvx8fHF19eXvXv3MNX9QS69+OQI9HVMMTYsOP6lmFVF1EDlmiM5d2A8XUsVnHjWuVBtdLT0iYvyx6JwTttTk+OQUmJhkH+AK0Bh45JEp2Q6bOY6hdHXNUNLlf+vjVqque6/i+fqMFxaDMOt2i60tPUIf3CBu/tn0anrp9z38aZw4cJZfVavWYt+5QYYOLmQHh+DtrFZjjENnV2w+HQMX33zLWeOHUFHR6fAOQMkJiaye/du/P39MTIyom3btrnqs5w8eZL12/+g9PKr6JhaZsmFEJhUqInT/JNsGVKD9m1a0axZM5YvW8G8uQuJjo7D2NCCiKinuLq6UrtOdezt7XFzc6NmzZp5OnQKCgofJ8qm6wdIuXLlcHN15ea+yUgpSUmIxNS8eC69+LhQ9vzelbXzKnHn5haSk6Kpb9MDv7hrtBt3AYcq7lkOBoCFbQVajTiMlV1Fjh8/iUv9r3jwQy8izuxCqtVZeur0NEJ2L+Pp0gk0cx3L5VNRVHWtgYl+YYz1rHLNQ0/biOS0WDLUBW+/JCQ/JyUljqJ21UhPS8JAx6xAfSEE+rrxcQpEAAAgAElEQVSmZKQm5WoLenAy8xmkRBY4RmRSEBkyg433JvMw1ovk1BjUUp2v/mHvWYRrPafNtJuUrNsbHT0jVFraFCnbkEaj92NTszftO7+MdT569CjTps8h/OxBHg1sgneb4jz6piNxNy/kGNe67Rd4XLyIo4NTjqyrXl5ezJ8/n9mzZ3Pw4EHS09P5ZfpMbIvZM+eHTZzbEcOfq27jVqUmrT5x59mzlytccxYuwarbtzkcjOxom5hj1X0MsxcupvennzN/1gZauv7GxN4BdKq3BgtjJ27euM+JA6FsWnOTTh36ULFCVc6fP1/gM1VQUPh4eNfbJQrviG1bNlC/YVPOrX6MVakGREc+ztGeEPeMzUvr4ezWg7591uDntQP7BEtOBm+kVM1emBXNe9VDpa1DFfdpnFnRmdq1y9O6/x+cXTmCgCVjMa3SkPSUBOKvn6GweWm+7nqWYoUyF5pi4oNZtr0xB2/NpI1LzmJkxvpWFDUty64zIzAztsVQ3xKXkh0wNSqaQ+/SndVoqbRJS43HzKokTyKv42xdO99nkJgaQ1R8IEYWdjnkcVFPubR9FC5m9Tn1eC3tyuZ1ijqTk49XY1PIhT3+C0jPSEFb2wCfwGNUKP5JLt30jFRO3PmVJhNPoqNnlKtdCEHV7jM5NN4JT09PTp06y+zZC2naeA7ly3dBW0ef1JR4vG9t4uSYLhQbOYtCrTOryWqbWaKta0TDst/TpHELJk2ewIoV6wgLi6SMkztaKj2WLd1BRGRftNGjU9X52FpWws6yMkIIulT7lUPeU6lbuyFXPC9iaWnJ6RPHqDRoRa55ZseqSTdOL/gax6KuDGp3El1tAwKeebJ8TyvaNplPlXI9sxxRtVRz5/6ftG3bkb17/6BBgwYFjq2goPDho6xkfKBYW1tz5fJ5+rZ35fHpWUSG3iEs+OU34LOHv6NE5Q7UajcVXX0T1Bmp6KkM8Iw6ikPNbgWMDEWd6yMl+Htvw8apPt3HXKNlr03EXzhC6chCfNX+OMM7n8lyMADMjG0Y3OU4B2/NICHlZarxdHUam6+MIjTGh+cJ/iSmx/Ao9CLTfy/HpmP9SEnLrMsXGOaF1/2t6OgY4XNzK5VqDuCY32+o1Rn5zvPswzXoGZjz7IkHibHPiA57wJXDP7FlekU+LTSCoY4z2X9vNsGx9/Lsfz/iEhcCttC3xSYGtt6NrrYRanU6a09/Tmzis9z6IacxsLTDzLZcvnNSqbSwr9uPmbPnMnPmfHr3Po9L5d5ZsTK6esZUqzaYL3qfImjeaJKe+GY+p/gY0lMS8A06SVx8POPGTcLZYQjDBz2idbPFfNJkLl/0ukL3DrtJSktmz90ZLD7dlUm7ynPOdzU6Wvq0rzoNO8P6TJ+emfsuLTUFlb5hvnMFUOkbkp6WRotq09HVNkBKyfaTg3FvupCqFXrnWOlSCRWVynSmY9O19P1sAGp1/is+CgoKHweKk/EBY2pqyncTxhP09DGTJ33Hyd39SU6KJjkxige3d+HW7GVMg3mRsvjGXCZNnYJ2Ht/CsyOEwMDIlECfY8SE+SGEQAiBrrYBHZospGihCnn2MzctTmmHFlx6+DuQ+c13+dleBCXcZ9wAP/p33E+bBjPp2XoD3w14TAYZLNvTmjM3F7FkTwv0DcxpXG0sF47/iHWxyuhbFGX11SF5Ohp3Q0+xy2c61q36cuHET2yaWYldS5sSkhaIltCmlnlLyppUY6jDTKacrMcen5nEpWQWVYtMDGL77R+YecGdz5pvxNLEnjJ2TXAp2QFbswokpkUzfntJVhzvmeUEhUT5cPjGDHQtXp3CxcCiOB5XrlOj5jhMzezy1ClkXZZqroMI37E0c/z1s9DW0sPAqRI2JerRrOFUargNyRUf4lC8Lv17nSQhPpSuE7yo22c5hx8sZOfVzM+6WYWxrF29lpSUFBycyhB/x6PAucbf8UDb0ARnu0YA+D+7QmJKNC7luufbp2ypNqgzDDh58uQrn4WCgsKHjeJkfASoVComTpyAe5t6bF3qxvmjkyhk44KByctAS1unhiSrUjHVtiL80ZUCx0uMCSU5PpLp06bx57x6PHvsQVSoD/bFaqF6RY6NEsUbcDfkBADegQcJibtH3/a7MDEsnEPPQN+c7i3XkiFTuXBvDVaFy9G42lia1ZyAjtBj07LauDX4hkB1AGMOlOeQzzzuhJzgiv9O5p5pxdwzbXH5fjule/9A7QWXaboljEYbAqg0YgV2bYawKXQuAK2K9KWorgN/+Exl6H47eu3UZ9SRcoSJMEZ2Pkd5h5e5LWqX+4KUjHj0jawwsS/P3fDTjNpgxbBVuszd14gnYVdJjnx1ldak548JCfLHxaVPgXquLp8Rc+pPkvzvE7F9MY06LaRM1Z6EB16nuuugfPsVLeKCvV1d/Ly2YevUgPajTuP5dBd3A49RxMwZA11zHj9+zIghg4jc9WuBc4j4YwHG2oZZwZxPgi9StlSbAj9nIQROxdty4cLFAsdWUFD48FGcjI8EIQSLFs1nx7Y16Etv0jNScrXX7DQL/5ib3Du+CHUBOTDunf2Nrl27MWLE11SuWJpDv7njefBnUpKjXzmPtPQk/CLOMO9UI3Zdn0T9at+gra2Xp65KpUWLOj+SlPAcdXIitSoNBKBZjYnYFjHi+rlviU94iJ6FEft8Z7DMoxdnI2fy+ehWaOvrYGSXdwIrh87fcCn1DGsCfiRdpvM06R662obMG5zIrIFRzB4YQ/dGSylqUTZHPzMjG1LSE0mMDcF1+hGklkBH15CWFccyv2coE9peICbsPlFPb+V5XwB1Rjr+F9Yh1Rno6xccuGpoWIi0+Ch8+9XC2Lgo5ap/xrOAqzjYN0BHx6DAvmVKtuLZ48xVCn1DC6q0GMcJ3yVZ7RkZGRQrVpTkOxfxGdOGhEe3c/SXUhKy8Re0/b1JS0sgOSU2U44aldDiVaiEtrJdoqCgoDgZHxNCCBo2bMi2bZtIiHxAWkpCjvYSLu7U6b6ImJinnF3TF3V6Wq4xnlzbyaMLy5n43VgAfv11Lvq62tSpPoJHT8+SnBpX4BzuBfzB2vUr+X7GUCIS/SjrWHAWTGeHZsTFBzGgwwH0dDK3cQz0zXBwKIXfgzscObSTNStncf7icSJjQ7hx+yojRgyne/ceBOyZn+eYuqZWVPrhD7YHL6Dz1eJkyAxS0xNISolGV9sg3yOY4bEPMdazQkffBF1TS2xbD8KyaAVuBO1DCEFxSxfaVZrI5RVfkJaU+zlIKbm1czwVy5fFxMScyMiHBdoeHuGDvrYJrsXbUbZ6n5fzeo0jon/N2upctTs+AYcJi31IdHwoTVu0ZtTkuRSt8SnmekW4O7QBXj3LEfzHYkK2zedev0oYXNnFxdMnada0BVd91wNgZ10VP/8Tr0zS5R96nKpVXV85TwUFhQ8bxcn4CClevDh16tTF5/KaXG1lqn9Kr4k+RAXeZvMYG67v+5En1//E9+xKjs1vwO19ozlx7BClSpUCoHr16gwfMZQbNxZTrKgrp6/OzjXmC3weHSQpNZSOHTvSvXt39PXyXsH4K0KoMNS3yHofHHGd0qVLIoSgSpUqNGvWDBcXlxzOwdQfJxN3eTtPds3J5Swlhj7Ga3oXkmUS8TIBbV0jSlhU45JP7ueRnUu3V2JhYEvR2pkl4gtVb0lKWjzhcY+yYjPaVBpPecNqHJlUnQenVpIS95y0lASCvY9wYUErZNAZ/ty5hX79Puf69aUF3u/a1WU0q/wterqm6Opl1hgsbFeVJ/5nSEtPLrDv/cdHsC7+MnGXtp4RGeo0Np4fAHpGuAzdTr3Jl3DpPZ8ag9fQYdkzyjYeQvCK73EJuczOlYu5de0KdnZ2TP5hAie8fuZR8AWc7BqRkZHK/cdH8r33k6CLxCU+pU2bNgXOUUFB4cNHcTI+UubOno73iZ/xu/FH7m+lAvT1tOnZrQPVbSJQPVlPkZQzzPlxOAGP/XB1zfkNdcqU75kzZwoq7VDOXZ/PgbPjSEzOdoIkPYUrt9bwx8m+bNu+CS2tzOX2qlWrc9//aIHz9As4SRHL8mhr6QKQmpaIp88aBg0uuKqrjY0NHhfOYuy7n4sDHLm3ciT3f5/MtUmtuTS8GrXL9mfK+CQG97uEvWN9opKCOHZtBoHhN/Ic7+q9zTwKucD9qEvYdvwq27MS6OjosOpCd1LTkxBC8FmtpXzuupCEU/vZPcqR7UMsePjnSL4b1gWPi2exsrJi5MivuXvnd+7fP5Dn/by81vLo4TFqlvmMwqZOhD+9DoC5tRN6+qZcv5G/QxQe4csT/9OUrt7rpSzgGoYGhfALv0yDiSexcqqZo49KW4fSrUZQtt040tIyaNSoUZbTVrVqVbZu28iGo+3ZerIX1cv2YduBPjwMOJ3r3v7Bl9l2sAu//bYwV9VcBQWFjw+h1CYomGrVqklPT89/expvBU9PTzp37YlaZYJdha7o6BoRFXSFJ3cO8O3o0UyePPFvZW+UUnLixAmmTZ3J1atXKWVfGy0tXR4HXsbFxYV582fg5vYyDfjevXsZ8eUUhnS+nGdchlqqWfVHSyo6daJu5SGkpSWx5VgPyruas3nz+tee161btzhw4ACXPTy4cvk+X3zmgZ7ey3LsGep0tm5pT0ZsFEGxPtQs25e65QdibFiYsOj7nL/1Gz4BR1BrCUr0n4pdm8ygS78NP5DueQlb6zTsbW05eewM9UoNpkShWqRlJHMjaCc3Avaw9LfF9Or9aa55eXh40Lx5K2xtG1Kt+jBMTGyIjnrM9esrCQq6irGBNUlJzzHSsyIs9gH9Jj3GyKQIq6bYok5LoXXzX3Gp2DNHEGZwqBeb/uiIW6tJlK/z0hE7tKozMUEXMCtdn+pf78j3WaUlxnJopAP3fW5ja2uboy06Opq1a9exc8deQkKCCAsLx9qqFM4O7ggE/iFHiYp9yLLli+nYseNrfz4KCgrvH0KIa1LKav94HMXJKJgP2cmAzOqbx44d49DhoyQnp1K+XGn69OmNhYXFqzsXQFhYGB4eHmRkZFChQgWcnZ1z6ajVatq7dyHgYTqdm6zFyOBlNtDklFh2nfyS4Gc3aVP/F0IivLhydzlNmjZgw8Y16L3mVkt2kpOTqepak2JFO9Kw3pQcDlRaejJbt7TnScAZrIztiEmOIC0jBW0tPdLSk7Cq0BCHPhOxdGkIQGpcJJcGVsRAywhtVTxeXp6EhYWxdMkK7t7yRVdPl5atm9Lns96EhYWRlpZGqVKlsLTMmV0zJCSEmjXrEBcn0dLSw9ioKFVc+lCpQk8SEyPY8WdHKrnYYW1tzekLPrTqd4DfZ1Ska6s1HDo1lrSMJMqVbo+Wlh6P/E8SGfOY2u1nUbZm5skVKSXepxfgd/lXatSuRYh1S0o2/LzA53RlXmsWTBpKu3btCtTLyMjgwIEDXLrkgZSS6tXdcHd3f6205woKCu83ipPxjvjQnYx/m9TUVEYMH83vv2+kfKk2mBqWJCE5iFsP/qRMmbIYGZmQkZ5B+QplGPblIKpUqfKP7hcSEkLLT9oRGyupVGEQVlaliY8P5a7POiKj77Jhwxru3buHv78/q1avxaDKJ5QZNA+9bPkvkkKf4D21G1rxyZgZFcPRsR4pyac4f/5lXoiEhASmTf2FlStWYahnha62PqGRD2ndqg0/TZ2co3R7XFwcgwd/xb59eynj3BoDgyJER98jIPAyI0eMYMoPmUXoRn0zhnXrN6KlbU6retOpWKYz/kEXeRxwmsSkSG4/2IW2vgnl6w3CxMqR+KinPLq2GgOdFI4c2seYCZN4at6MEg37UhBX57dl3sRBuLsXXAhOQUHhw+U/42QIIbQATyBIStlWCLEaqEZmRdX7wOdSynghxDfAACAdCAe+kFL6CyGqAL8BpkAGME1KuU0zdglgK2AFXAP6SClThRB6wAbADXgOdJdSPtH0mQD014w1XEqZfwQbipPxroiIiGDnzp08e/YMS0tLOnXqlGu5/k2hVqs5cuQIq1ZtICgoBHNzM3r16kLXrl3R139ZpTYsLAz3zt3x9PSkSL3O6JkXIe7hTaJ9LiDQonSpVrRvvxqVlg5Llzhy+vRhKlWqRHx8PI0aNEM7yZGWLlMoZpGZATQhOZLzvis4c38+x44fylVRNTQ0lH379hEbG4uNjQ3u7u4YGeVMjObn58eIEaPwvhHKsD5XcqzGqKUavyfH8bq7iachlyla1JhZM3+mZcuWqFQq5s2bz7K913Ab8nu+zyY9OYFDI+256+2Fvb39m3jcCgoK/0H+S07GN2Q6FaYaJ8NUShmraZsHhEkpZwghGgMeUspEIcRQoJGUsrsQojSZBVQfCCFsyHQmykkpo4UQ24FdUsqtQohlwE0p5W9CiGGAi5RyiBCiB9BRM1Z5YAtQA7ABjgOlpZT55qZWnIyPm9WrVzN/wR5MTUoTFHwVlUobO7tauLp+gaVlqSy948e+pUULMyZPnsSQwV9x+0I8vevmXWL+2qOdHL47lsf+D7KCYP8Oqamp1KndEEOdmrRsODdHam8pJVduLOOK9wy8blzNUe01MjKS4o4lafLzdYyLlMxraO4dWkCh8JMcPbj3b89LQUHhw+FNORlvNfxbCGEHtAGmAd8AZHMwBGAASI38VLaul4HeGvn9F0IpZbAQIgywFkLEAE2AFxF164EfyFz1aK+5BtgJLNbcrz2wVUqZAjwWQviR6XBcepN2K3w4xMXFYWrqSIsWcwrU09cvxKXL53F1q8Wdmzf5uefDfINm3Up24bRvZtXUV8U95IWuri5Hjx2gQ4duLFzrTJXyA7CycCY2PoTb99eho5fCmbMncjgYAJaWlkyfNpWfZ35CjRG7MSv+Mv27VKt5cm4jjw/8wqZzp//2nBQUFBTy4m2fMVsAjAVMsguFEGuB1sBdIK8SmP2BQ38VCiFqALrAQzK3SKKllC9SUwYCL9bXbYGnAFLKdI1DYqWRX842ZPY+Cgq5sLe3JzJy3yv1QkKuERb9EIcq3XAIM8TcyKZA/SrFP2XP7gP/l5MBmQ7DmTPHuHr1KmvWbCDw6TUsi5rz2+jptGjRApUq79PpI77+Cn1dPcZ/1xhzhyoYl6yDOi2JZ9f+oLClKWdPHadcufwLvCkoKCj8Hd6akyGEaEvmVsg1IUSj7G1Syn6aWI1FQHdgbbZ+vcncXmn4l/GKARuBvlJK9d85Wvl/zH0QMAhQ9qU/ctq0acPAgUOJCPelkHXZPHUSEsJ54HeYT6c+Jvj+KdLu33nluAa6ZiTEJ/6juQkhqFGjBjVq1Phb/QYPHkjfvn34888/uXvXBz09E5qO30itWrX+1pFlBQUFhVfxNlcy6gLuQojWgD5gKoT4XUr5YhskQwixlcyVjrUAQohmwESgoWZLA43cFDgATJRSvliJeA6YCyG0NasZdkCQpi0IKA4ECiG0ATON/gv5C7L3yUJKuQJYAZkxGf/4SSj8Z9HT02PChPEsXNiTHj2OYWhUKEd7ako8W7e4U77BEAxMCmFSqATXI26gluoCi4gFR9+gdr284yLeBfr6+vTs2fNfu7+CgsLHwVvL+CmlnCCltJNSOgI9gJNAHyGEE2TFZLgDvpr3rsBywF1KGfZiHCGELvAnsEFKuTPb+BI4BXTRiPoCezTXezXv0bSf1OjvBXoIIfQ0J1OcgYJLjip89IwePZIePVqxalUlzpyeTFDgFYKDr3HhwkwWL3EmQ0+Hmh1nAGBt74bKwBifpzkzmaZlpGRlVk1OjeOK30b6D+j3zm1RUFBQeJe867y/AlivWZkQwE1gqKZtNmAM7NAs2QZIKd2BbkADwEoI8blG93Mp5Q1gHLBVCDEV8AJWa9pXAxs1gZ2RZDo5SCnvaE6k3CXzqOyXBZ0sUVCAzG2JWbOm07dvLxYt+o1z54YipRo3N1ccHG1xbDQ564SHEAK39j+zftMA+jfcyK3AQ1y8t5bklGhUQosKjm1JSA6ka9euODg4/MuWKSgoKLxdlGRcr0A5wqpQEK3bdiLFuhOla/bOIb+ydzLex+dSru4AKjX6GlPrUqQmxXDv8nq8jkzjh0njGDMmr5hnBQUFhX+fN3WEVSmQpqDwD+jTqyv+19blkKUkRuN7YSXNB2yjXrdfMSvshBACPUNzXJqMoPO4a8yYs5Bjx47lGi8qKop5c+dTr2ZjXCvWoJN7N44cOYJarX5HFikoKCi8ORQnQ0HhH9CpUycSI/3wu7olS3b/8gZsSjfCsVLbPPsYWxanSsuf+WXmvBzyw4cPU8rRmf3Lr9NQdxzdCi+icFBzvv58AnVqNiAyMvKt2qKgoKDwplG2S16Bsl2i8Cq8vb1p2qwlxSt3p2z9Lzm57jOqt/2R4uWa59snPS2Z3ycUJijQHwsLC65du0bzJq0Y7baHMla1c+iqpZqNd74h3NiTix5n882BoaCgoPCmULZLFBTeE1xcXLjhdZXGlXU5uKA2USF3MDIvOMebto4+RiaWREVFAfDT5Ol0KjUll4MBoBIq+lSYx/PA+Dy3WBQUFBTeVxQnQ0HhDWBra8vcOTOJjHhGpUqViIt4XKB+WnI8CXGRWFlZERERwanTJ2ho/1m++iqhoonNUJYvWZ2vjoKCgsL7huJkKCi8QVQqFYMG9OXhlZUF6j3w3EzjJs0wMzMjODgYa9PiGOiYFNjH3syFJ4/93+R0FRQUFN4qipOhoPCG6dXrU2JCvLjvsTHP9qhQX24ensJ34zOPsBoZGRGfHMWr4qMSUqNylX5XUFBQeJ9518m4FBQ+eIyMjDh+9CBNm7Xkmd9xnGsPwaJoOZLjI/C7ugHfC8tYuGAudevWBaBkyZKYWZhyN+IsFawb5jvuxWeb6DTU/V2ZoaCgoPCPUVYyFBTeAhUqVODunZv0aV+Jm7v7s3WKPUeXNqSqfRSXL56hb98+WbpCCEZ9O5xt98eTkpGU53j3Ii7iHXaEz/v1zbNdQUFB4X1EOcL6CpQjrArvArVaTZ9PP+fGOT+6OU2nvHVDhBAkpsVyJmADux/+xKZtG2jZsuW/PVUFBYWPgDd1hFXZLlFQeA9QqVRs3LyO5ctWsGDul0R7x2CkZ05EXBBNmzTj2KrDVK1a9d+epoKCgsLfQlnJeAXKSobCu0ZKyaNHj0hISMDW1hYrK6t/e0oKCgofGW9qJUNxMl6BECIceJ1zg4WAiLc8nXfNh2gTKHb91/gQ7foQbQLFrv8aBdnlIKW0/qc3UJyMN4QQwvNNeH3vEx+iTaDY9V/jQ7TrQ7QJFLv+a7wLu5TTJQoKCgoKCgpvBcXJUFBQUFBQUHgrKE7Gm2PFvz2Bt8CHaBModv3X+BDt+hBtAsWu/xpv3S4lJkNBQUFBQUHhraCsZCgoKCgoKCi8FRQnQ0FBQUFBQeGt8NE7GUKIMkKIG9lesUKIkdnaRwshpBCikOa9EEIsFEL4CSG8hRBVs+n2FUI80Lz6ZpO7CSFuafosFEIIjdxSCHFMo39MCGHxL9nUS2PLLSHERSFE5Wy6TzTyG0IIz2zyPOde0PP5F+xqJISIyaY/OZtuSyHEPc08x2eTlxBCeGjk24QQuhq5nua9n6bd8V+0a0w23dtCiAwhhKWm7b3/vIQQPwghgrLJW2frM0Ezl3tCiE+yyd+Lz+vv2iSEaC6EuKb5TK4JIZpkG+u0xqYXfQq/au75PZ9/wS5HIURSNvmybGP9rb9379O/QZH5tzC7vloIUUXT9l58XvnZpGn7WgjhK4S4I4SY9ap5vJPfKyml8tK8AC0glMwkJADFgSNkJuMqpJG1Bg4BAqgFeGjklsAjzU8LzbWFpu2KRldo+rbSyGcB4zXX44GZ/5JNdbLNtdULmzTvn7zQ+8u4ec49v+fzL9nVCNifT9+HQElAF7gJlNe0bQd6aK6XAUM118OAZZrrHsC2f8uuv+i3A/7X3p1HyVWWeRz//kwgQAIJAUUkeBKURQgYQsQEAQmgRkbB0YzKOIPOcA6iuICDyAhicGYOosygiIDAYMSFVWCEMwMiiwsYQhKydAhLgAwhQEIg7E5Y8swf71vpm6Kru7rTt6ty8/ucU6dvvXd7n7q3qp9631v3vXVDOl7ANODELpbZPR+LIcCYfIwGtevxajKmvYG35emxwLLCvNuBCV2s02XdG70+LYprNNDRYP1efd610zlYt/yewEPtfLzqYpoM/A4Ykue9pR3eVxt9S0adQ0gnVe0On2cDJwHFq2OPAC6NZAYwQtL2wIeAmyPimYhYBdwMTMnztoqIGZGOzKXAxwrb+lme/lmhfEBjiog7c50BZgCjmthuo7o3en36WzPHqpF9gcUR8XBEvAJcDhyRv3EdDFydl6uPqxbv1cAhtW9o/ay3cR0JXNbEdtvteDWq4+URsToiHgEWk45Vux6vHmOKiHsi4vH8dCGwuaQhPWy3Ud0bvT79rZlj1aU+ft610zlYdCTpXOtJK49XMaYvAN+NiNUAEbGiUL+Wva+cZKzr0+QPbElHkL51zKtbZgdgaeH5Y7msu/LHuigH2C4insjTTwLb9UMM9ZqJqeho0reKmgB+q9TUe0yhvFHdG70O/a3ZuCZJmifpfyTt0UMdtwGejYjXuqj72nXy/Ofy8v2t6eMlaQtgCvDrQnHbH6/sS7lp/BJ1dhP29r3V6uPVTExFnwDm1P4JZD/NTd7fKnxYN6p7Ox0rgDGS7pH0e0kHFOre28+7dour5lO8MYFvt+NVjGkX4IDcjfF7Se+pr19dPQbkfeUkI8t9TocDV+UP728Cp3W/Vv/JWX+//p64tzFJmkxKMr5RKN4/IsaTulGOk3TgQNS9O72Iaw6p2+HdwI+A6waqjn3Rh3Pwo8AdEfFMoaytj1cuOh94BzAOeAL494GqS3/pbUw5wT0T+Hyh+DMRsSdwQH78fcnV7lEv4noCeHtE7A18DfiVpK2a3U+7n4OS3pkIK5sAABHMSURBVAu8HBEdheK2Ol5dxDSY1F0/Efg6cGVJra294iSj04dJ3zKWk06+McA8SUtI3QdzJL0VWEbqJ68Zlcu6Kx/VRTnA8lqzYP67gv7VbExI2gu4GDgiIp6ubSAiluW/K4Br6Wzua1T3Rq/DgMcVEc9HxIu5/v8NbKJ08WSjOj5Naqod3EXd166T5w/Pyw94XIXl67+ZbQjHi4hYHhGvR8Qa4KJCHXv73mrl8Wo2JiSNIh2LoyLioVp54Vi9APyKLl6Hurq3zbHKTe9P5+nZpL79Xejb513bxFXQ3XurXY7XOjGRWhyuyd1OM4E1pAHQWvu+ihIusNkQH6T+qH9oMG8JnRcT/hXrXqQ0M5ePBB4hXfS5dZ4emefVXwh1WC7/PuteCPW9FsX0dlI/3X51ywwFtixM3wlM6a7ujV6fFsX1VjpvOLcv8Giu12DShblj6LzgaY+83FWse8HTF/P0cax7wdOVrYorPx8OPAMM3dCOF7B9YfoEUn8xwB6se4Haw6SL09ruePUiphG5vh+vW39w4TzdhNS/fWx3dW/0+rQorjfX9k26cHAZffy8a6dzMD9/U45np3Y+Xl3EdCzwnTy9C6lbQ43qwQC9r/r1QG6oD9IH8tPA8AbzlxROMAE/JmXuCyhcbQz8I+mf9eK6gz8B6MjrnEvnP75tgFuAB0lXBY9sUUwXA6uAufkxK5fvlE+8eaSL1k4prN9l3bt7fVoQ15dyveeRLmjdr7DcYcADuZ7FuHYifUguzm+02pXam+Xni/P8nVoVV37+OQofjBvS8QJ+nvc1H/gN637gn5Lrcj/5Vwntdrx6ExNwKvBS4b01F3hL3sbsvPxC4Id0/tNuWPdGr08L4vpErvdcUrfkRwvr9Orzrg3PwYOAGV1so22OV4OYNgV+kV/7OcDB7fC+8m3FzczMrBS+JsPMzMxK4STDzMzMSuEkw8zMzErhJMPMzMxK4STDzMzMSuEkw8yapjTia3EEyJN7XusN2zhI0n6F58dKOqqf6lcbhXZCL9Z5R47lxf6og5l18k9Yzaxpkl6MiGHruY1pwIsRcVb/1GqdbS8h3WdhZR/WXe/YzGxdbskws/Um6TRJd0vqkHRhbcwESV+RdG8eiOpySaNJdyY8IbceHCBpmqQT8/K3SzpT0kxJD9QG3pK0haQr87auzYNA9dhakVs2zsj7miVpvKSbJD0k6djyXhEzg3RbUTOzZm0uaW7h+RkRcQVwbkR8B0DSz4GPANeTbh89JiJWSxoREc9KuoBCS4akQ+r2MTgi9pV0GPBt4FDgi8CqiNhd0ljSnSab9WhEjJN0NjAdeB/pzoUdpFsmm1lJnGSYWW/8JSLGdVE+WdJJwBakcXwWkpKM+cAvJV1H86PgXpP/zgZG5+n9SbdyJiI6JM3vRZ1/k/8uAIZFGuDqBUlrE59ebMvMesHdJWa2XiRtBpwHTI00FPZFpJYCSINf/RgYD9xdGNmxO6vz39fpny9Cte2tKUzXnvuLllmJnGSY2fqqJRQrJQ0DpgJIehOwY0TcBnyDNHLsMOAFYMte7uMO4JN5u7sDe/ZDvc2sZM7izaw36q/JuDEiTpZ0EekahyeBu/O8QcAvJA0njbR5Tr4m43rgaklHAF9ucr/nAT+TdC9wH6k75rl+iMfMSuSfsJpZ25M0CNgkIv5P0jtIQ4XvGhGv1C23BP+E1axtuLvEzDYEWwB/kjQPuBb4Yn2CkT0F3NKXm3EBy/unqmZW45YMMzMzK4VbMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFE4yzMzMrBROMszMzKwUTjLMzMysFINbXYF2pylTgpUru54ZoO5WjryNvi4TdXXpYZmG++lpmTbdTzvVZX33M1DnyYbyemxUdenFNvrjPKnS67Fe+6lwXfp8DjS7DWA2s2+KiCk9LNYjJxk9WbkS3T0LANUOTuFvV2VVmNcOdWinee1Qh3aa1w51aKd57VCHDWVeO9Shnea1Qx3q5wEIbUs/cHeJmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVwkmGmZmZlcJJhpmZmZXCSYaZmZmVQhHR6jq0NUk3AtvWFW8LrGxBdQaa46yWjSVO2HhidZzV0k5xroyIKeu7EScZfSBpVkRMaHU9yuY4q2VjiRM2nlgdZ7VUMU53l5iZmVkpnGSYmZlZKZxk9M2Fra7AAHGc1bKxxAkbT6yOs1oqF6evyTAzM7NSuCXDzMzMSuEko46kSyStkNRRKBsp6WZJD+a/W+dySTpH0mJJ8yWNb13Ne0fSjpJuk3SvpIWSvprLKxWrpM0kzZQ0L8d5ei4fI+muHM8VkjbN5UPy88V5/uhW1r+3JA2SdI+kG/Lzqsa5RNICSXMlzcpllTp3ASSNkHS1pPskLZI0qWpxSto1H8fa43lJx1ctTgBJJ+TPoQ5Jl+XPp0q+R2ucZLzRdKD+t8EnA7dExM7ALfk5wIeBnfPjGOD8Aapjf3gN+KeI2B2YCBwnaXeqF+tq4OCIeDcwDpgiaSJwJnB2RLwTWAUcnZc/GliVy8/Oy21IvgosKjyvapwAkyNiXOEnf1U7dwF+CNwYEbsB7yYd20rFGRH35+M4DtgHeBm4lorFKWkH4CvAhIgYCwwCPk2136MQEX7UPYDRQEfh+f3A9nl6e+D+PP0T4MiultvQHsB/AR+ocqzAFsAc4L2kG94MzuWTgJvy9E3ApDw9OC+nVte9yfhGkT6MDwZuAFTFOHOdlwDb1pVV6twFhgOP1B+XqsVZF9sHgTuqGCewA7AUGJnfczcAH6rqe7T2cEtGc7aLiCfy9JPAdnm6dtLUPJbLNii5GW5v4C4qGGvuQpgLrABuBh4Cno2I1/IixVjWxpnnPwdsM7A17rMfACcBa/LzbahmnAAB/FbSbEnH5LKqnbtjgKeAn+YusIslDaV6cRZ9GrgsT1cqzohYBpwFPAo8QXrPzaa671HA3SW9FimtrMxPciQNA34NHB8RzxfnVSXWiHg9UlPsKGBfYLcWV6nfSfoIsCIiZre6LgNk/4gYT2o6P07SgcWZFTl3BwPjgfMjYm/gJTq7DIDKxAlAvhbhcOCq+nlViDNfU3IEKXl8GzCUN3bNV46TjOYsl7Q9QP67IpcvA3YsLDcql20QJG1CSjB+GRHX5OJKxgoQEc8Ct5GaJEdIGpxnFWNZG2eePxx4eoCr2hfvAw6XtAS4nNRl8kOqFyew9lshEbGC1H+/L9U7dx8DHouIu/Lzq0lJR9XirPkwMCcilufnVYvzUOCRiHgqIl4FriG9byv5Hq1xktGc3wCfzdOfJV2/UCs/Kl/tPBF4rtC819YkCfhPYFFE/EdhVqVilfRmSSPy9Oak604WkZKNqXmx+jhr8U8Fbs3fotpaRPxzRIyKiNGkJudbI+IzVCxOAElDJW1Zmyb143dQsXM3Ip4ElkraNRcdAtxLxeIsOJLOrhKoXpyPAhMlbZE/f2vHs3Lv0XW0+qKQdnuQTvIngFdJ3ySOJvWD3QI8CPwOGJmXFfBjUh//AtJVwy2Pock49yc1P84H5ubHYVWLFdgLuCfH2QGclst3AmYCi0nNs0Ny+Wb5+eI8f6dWx9CHmA8CbqhqnDmmefmxEDgll1fq3M11HwfMyufvdcDWFY1zKOlb+vBCWRXjPB24L38W/RwYUsX3aPHhO36amZlZKdxdYmZmZqVwkmFmZmalcJJhZmZmpXCSYWZmZqVwkmFmZmalcJJhNoAknZJHYZyfR5x87wDs81pJHys8v1/SqYXnv5b08W7WP17SFutZh2mSltWNtjmij9s6VtJRfVx3uqSpPS/Zv/J+H5F0bAnbPkBpNOWOnpc2G1hOMswGiKRJwEeA8RGxF+kOgEu7X6tf3AHsl+uwDen21JMK8ycBd3az/vGkweWaJmlQF8VnRx5tMz+e7c02ayLigoi4tC/rttjXI+KC/t5oRPyRdI8bs7bjJMNs4GwPrIyI1QARsTIiHgeQtETS6ZLmSFogabdcPlLSdbnlY4akvXL5Akkj8l0Pn659s5d0qaQP1O33TnKSkf9eD7w5rzsG+EtEPCnpfEmzckvL6Xl7XyGNs3CbpNty2Qcl/TnX9ao8/k0thjMlzQH+ppkXRNLmki6XtCi3uNwlaUKe92JhuamSpufpaZJOlLSbpJmFZUZLWpCnT5N0t6QOSRfmOyzW73sfSb9XGmTtJnXewvr2HMdMSQ9IOiCXD5J0Vt7mfElf7m47PcQ9Pb/eMyQ9LOkgSZfk12F6YbkXJX0/H5PfSdo31+9hSYc38xqbtZKTDLOB81tgx/yP6zxJ76+bvzLSoF/nAyfmstOBe3LLxzeB2jf4O0jjHuwBPAwckMu7apWYDYxVGoBqP+DPpOGx35Wf15Y/JSImkO6S+n5Je0XEOcDjwOSImCxpW+BU4NBc11nA1wr7ejoixkfE5V3Ef0Khq+S2XPYF4OWIeBfwbWCfRi9evYi4D9g0J0oAnwKuyNPnRsR7ImIssDmpBWktpXF7fgRMjYh9gEuAfyssMjgi9iW14nw7lx0DjAbG5ePxyya2052tScfrBNItpM8mHc89JY3Lywwl3U56D+AF4F9Jt8b/a+A7Te7HrGUG97yImfWHiHhR0j6khGAycIWkkyNiel6kNkjdbKB2jcT+wCfy+rdK2kbSVsAfgQOB/yUlJcdI2gFYFREv1e13taSFpMG1JgLfI93KeD9gb1LCAvBJpWHTB5NaXXYn3c66aGIuvyM3DmxKSlpqrqCxsyPirLqyA4Fzcj3nS6rfX0+uJCUX381/P5XLJ0s6idTNM5J0+/HrC+vtCowFbs5xDCINJ1BTPBaj8/ShwAWRh+WOiGckje1hO925PiIit74sj4haK8zCvM+5wCvAjXn5BcDqiHg1rzP6jZs0ay9OMswGUES8DtwO3J7/UXwWmJ5nr85/X6fn9+YfgOOAtwOnkL7ZTiUlH125g/QPfcuIWCVpBvAlUpLxk9wacCLwnjx/OmnshHoCbo6IIxvs56UG5X1RHPOgq7pASmquknQNaUTwByVtBpxHGtNiqaRpXawvYGFETKJrzR6LnrbTndo+1hSma89r+3w1Osd+WLtcRKxR58idZm3L3SVmA0TSrpJ2LhSNI7VEdOePwGfy+geRulSej4ilwLbAzhHxMPAnUpLwhwbbuRP4PGlQMUgtFBNJSUoHsBUpQXhO0nakYbdrXgC2zNMzgPdJemeu01BJu/QQQ3f+APxt3tZYUldNzXJJ75L0JlIS9QYR8RApEfgWna0otYRiZb5epKtfk9xPui5lUt73JpL26KGuNwOfr/1zlzSyj9sx22g4EzYbOMOAHyn9dPM10uiKx/SwzjTgktyN8DKdQz8D3EVqnoeUjJxBSja6ciepi+QMgIh4TdIKYGlErAHmSbqHNELkUjq7UAAuBG6U9Hi+LuNzwGWShuT5pwIP9BAHpGsy/q7w/GOkrp6fSloELCJ1T9ScDNwAPEW69mNYg+1eAXwfGJNje1bSRaTk6Ung7voVIuIVpZ+yniNpOOmz8AekbpVGLgZ2AeZLehW4KCLO7cN2zDYaHoXVzNqGpNuBEyNiVqvr0p9y99MNEXF1Sdsfnbc/toztm/WVu0vMzMr3HPAvKulmXKSLWlf297bN1pdbMszMzKwUbskwMzOzUjjJMDMzs1I4yTAzM7NSOMkwMzOzUjjJMDMzs1I4yTAzM7NS/D+BOklWw0M+nQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Build a singl plot\n", "fig, ax = plt.subplots(figsize=(8,8))\n", "\n", "# Colorbar keyword arguments\n", "kwds = {'label': \"Snow Water Equivalence [mm]\", 'orientation': \"horizontal\"}\n", "\n", "# Plot it all up as a scatter using SWE as the color\n", "ax = point_swe.plot(ax=ax, column='swe', cmap='cool', markersize=100, legend=True, edgecolor='k',legend_kwds=kwds)\n", "\n", "# Set some style/ labeling choices\n", "ax.ticklabel_format(style='plain', useOffset=False)\n", "ax.set_xlabel('Easting [m]')\n", "ax.set_ylabel('Northing [m]')\n", "\n", "title = 'Grand Mesa Pit SWE\\n{} - {}'.format(point_swe['date'].min(), point_swe['date'].max())\n", "ax.set_title(title, fontsize=16)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Close the session to avoid hanging transactions\n", "session.close()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 4 }