{
"cells": [
{
"cell_type": "markdown",
"id": "aca93ab3-f452-4725-b864-b41cd0327954",
"metadata": {},
"source": [
"# Spatial partitioning in Dask-GeoPandas\n",
"\n",
"Dask-GeoPandas has implemented `spatial_shuffle` method to repartition `Dask.GeoDataFrames` geographically.\n",
"For those who are not familiar with Dask, a Dask DataFrame is internally split into many partitions, where each partition is a Pandas DataFrame.\n",
"Partitions are split vertically along the index; \n",
"\n",
"\n",
"\n",
"There are numerous strategies that can be used to partition Dask DataFrames, which determine how the elements of a DataFrame are separated into each resulting partition.\n",
"Common strategies to partition Dask DataFrames include using a fixed number of partitions or partitions based on column values.\n",
"Dask-GeoPandas can partition GeoDataFrames by taking advantage of their spatial structure.\n",
"\n",
"In the rest of this notebook, we use an example dataset to describe how Dask-GeoPandas takes advantage of this spatial structure and some best practices."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "4a5fbccd-acca-4988-83db-948bf1d2bf6a",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/tstatham/anaconda3/envs/gsoc/lib/python3.9/site-packages/geopandas/_compat.py:111: UserWarning: The Shapely GEOS version (3.8.0-CAPI-1.13.1 ) is incompatible with the GEOS version PyGEOS was compiled with (3.9.1-CAPI-1.14.2). Conversions between both will be slow.\n",
" warnings.warn(\n"
]
}
],
"source": [
"import geopandas\n",
"import dask_geopandas\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"id": "b6139e44-48d1-4906-b33a-eaa0c9b2bc3b",
"metadata": {},
"source": [
"We first begin by reading USA boundaries from [gadm.org](https://gadm.org/)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "a4d93451-bd17-4ca1-99c1-aed7f4d828b2",
"metadata": {},
"outputs": [],
"source": [
"path = \"https://geodata.ucdavis.edu/gadm/gadm4.0/gpkg/gadm40_USA.gpkg\""
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "dd57d7c7-91dd-4e09-937b-1707b7b2ce2e",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
ID_0
\n",
"
COUNTRY
\n",
"
ID_1
\n",
"
NAME_1
\n",
"
VARNAME_1
\n",
"
NL_NAME_1
\n",
"
TYPE_1
\n",
"
ENGTYPE_1
\n",
"
CC_1
\n",
"
HASC_1
\n",
"
ISO_1
\n",
"
geometry
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
USA
\n",
"
United States
\n",
"
USA.1_1
\n",
"
Alabama
\n",
"
AL|Ala.
\n",
"
\n",
"
State
\n",
"
State
\n",
"
\n",
"
US.AL
\n",
"
\n",
"
MULTIPOLYGON (((-88.12063 30.23321, -88.12122 ...
\n",
"
\n",
"
\n",
"
1
\n",
"
USA
\n",
"
United States
\n",
"
USA.2_1
\n",
"
Alaska
\n",
"
AK|Alaska
\n",
"
\n",
"
State
\n",
"
State
\n",
"
\n",
"
US.AK
\n",
"
\n",
"
MULTIPOLYGON (((179.20071 51.35992, 179.20010 ...
\n",
"
\n",
"
\n",
"
2
\n",
"
USA
\n",
"
United States
\n",
"
USA.3_1
\n",
"
Arizona
\n",
"
AZ|Ariz.
\n",
"
\n",
"
State
\n",
"
State
\n",
"
\n",
"
US.AZ
\n",
"
\n",
"
MULTIPOLYGON (((-110.87869 31.33084, -110.9720...
\n",
"
\n",
"
\n",
"
3
\n",
"
USA
\n",
"
United States
\n",
"
USA.4_1
\n",
"
Arkansas
\n",
"
AR|Ark.
\n",
"
\n",
"
State
\n",
"
State
\n",
"
\n",
"
US.AR
\n",
"
\n",
"
MULTIPOLYGON (((-91.90713 33.00801, -92.01710 ...
\n",
"
\n",
"
\n",
"
4
\n",
"
USA
\n",
"
United States
\n",
"
USA.5_1
\n",
"
California
\n",
"
CA|Calif.
\n",
"
\n",
"
State
\n",
"
State
\n",
"
\n",
"
US.CA
\n",
"
\n",
"
MULTIPOLYGON (((-117.23300 32.77416, -117.2330...
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" ID_0 COUNTRY ID_1 NAME_1 VARNAME_1 NL_NAME_1 TYPE_1 \\\n",
"0 USA United States USA.1_1 Alabama AL|Ala. State \n",
"1 USA United States USA.2_1 Alaska AK|Alaska State \n",
"2 USA United States USA.3_1 Arizona AZ|Ariz. State \n",
"3 USA United States USA.4_1 Arkansas AR|Ark. State \n",
"4 USA United States USA.5_1 California CA|Calif. State \n",
"\n",
" ENGTYPE_1 CC_1 HASC_1 ISO_1 \\\n",
"0 State US.AL \n",
"1 State US.AK \n",
"2 State US.AZ \n",
"3 State US.AR \n",
"4 State US.CA \n",
"\n",
" geometry \n",
"0 MULTIPOLYGON (((-88.12063 30.23321, -88.12122 ... \n",
"1 MULTIPOLYGON (((179.20071 51.35992, 179.20010 ... \n",
"2 MULTIPOLYGON (((-110.87869 31.33084, -110.9720... \n",
"3 MULTIPOLYGON (((-91.90713 33.00801, -92.01710 ... \n",
"4 MULTIPOLYGON (((-117.23300 32.77416, -117.2330... "
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"usa = geopandas.read_file(path, layer=\"ADM_1\")\n",
"usa.head()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "247440da-ada6-4c25-bcd1-e1f0eb221d0c",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"RangeIndex: 51 entries, 0 to 50\n",
"Data columns (total 2 columns):\n",
" # Column Non-Null Count Dtype \n",
"--- ------ -------------- ----- \n",
" 0 State 51 non-null object \n",
" 1 geometry 51 non-null geometry\n",
"dtypes: geometry(1), object(1)\n",
"memory usage: 944.0+ bytes\n"
]
}
],
"source": [
"# Filter and rename cols\n",
"usa = usa[[\"NAME_1\", \"geometry\"]].rename(columns={\"NAME_1\": \"State\"})\n",
"usa.info()"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "20cdb8ab-f58c-436b-af5e-919f9f8b41a0",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAABnCAYAAAAZvRmoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAWWUlEQVR4nO3de3hcZZ3A8e9v7pPMJJN70rRp0ltKW6iFkJZLAUGEBRSvKKLyKMiz6np5dt0VRF12cd11L674WEUUWBAUcUXh0Yf1goAK0hsUCi20TZPQpmmTNrfJZWYyM+/+MSftpE2atJ1LTvr7PM88mXnPZX7n5MxvzrznPe8rxhiUUkrZjyPfASillDo5msCVUsqmNIErpZRNaQJXSimb0gSulFI2pQlcKaVsypXLNysvLzf19fW5fEullLK9zZs3HzTGVBxdntMEXl9fz6ZNm3L5lkopZXsi0j5RuVahKKWUTZ1WCVzvOlVKzSanVQIXEWLxJG3dYeLxeL7DUUqpU5LTOvB8MMYQjsZ55C8ttBwc5qW9YVoODFLlhWXzSygtcHP9ubX09/QSKi9nWW0JHo8n32ErpdSUZnUCb+ns4eb7X6B9wJA8atq+KOzb0QvAz17uIlW70oYDqC/1EnInCHidxB1eLl9ew7Vn11FaqIldKTVzzMoEvrm9l56hGHc9+QqtA1PXe6dXjSeB3T1R61UciPJ82wA+r4frm+uyEa5SSp2UWZfAB0eifOt3r/OnXT0ZXa9kdG1KKXXqZt1FzIDfy49uPo8nP3shIZ8z3+EopVTWzLoEDvDirg7W7z7Ic198K/d9eCVr6otPeZ3aAFEpNdPMmiqURCLJzgP9fPaB59kzCCMJ+NH6PZhEPK1OWymlZo9ZcwbudDporAnx5WtX4nSmqk52dw8xEjcsn1OU5+iUUirzbH8GnkwmcThS30MiwkXL5vLcbZUMhcMUBIL43U7++YmtvLZvIM+RKqVUZtk+gW/Y3c2O/f0sri7hvEWpzrqK/R5ebBnCcXCE/31xH1s7h/IcpVJKZd60EriIhIAfAitIXc/7OPAG8FOgHmgDrjPG9GYjyONZs6iKNYuqxpW93tnP5x7aip5zK6Vms+nWgd8F/J8xZimwEtgO3Ao8ZYxZDDxlvZ4R/nrdnzV5K6VmvSkTuIgUAxcB9wIYY2LGmD7gWuABa7YHgHdlJ8TJDY1EJiy/9zNrmW7DQZcDbmiq4dIlpficcE2Dg5ubynBnLkyllMqK6VShNADdwP0ishLYDHwOqDLGdFrz7AeqJlk+K5LJJA7n+PCNMTz0QjvzSgr41OUNPNvSR9Dr5DfbD064DgE+sbqawaSLZ3Z0csWSYi5cVErT4mquWzXAxv2jdPZHSIiLJ19q49lte/n8lStYWl3EcCyOMVDotf1lBKWUTclUfWSLSBPwAnCBMWa9iNwFDACfMcaE0ubrNcaUTLD8LcAtAHV1dee0t084sETW/H7bAW597BUisQRxY0gkDaOJ1DZ/pHkOX7lmBT/+03aaFlaysKoYv98/rmVLuu17D/HQxg7Orivhvj+10FBWwDevb8LjmjWtMZVSM5CIbDbGNB1TPo0EXg28YIypt16vJVXfvQi4xBjTKSI1wDPGmMbjraupqcnkeki17oERDOBzOdjZPcQTW/bxwF9SXyJfvaKBrVta+fonL8fvm15Pg8YYHt3Qyn/8YjsNRRCsLuX7N67G7dQkrpTKjskS+JRZxxizH9gjImPJ+TJgG/AEcKNVdiPweIZiPWGj8cSk0yqK/FQW+Skq8NJQHuBPuw5yVqWbxz+1hqdf78JZWjjt5A2ptuYfWL2AxQ4oLha6dvSw7g+7MrEZSil1QqZbgfsZ4GER8QC7gY+RSv6PishNQDtwXXZCnJoA8UQS13HOgp9vOch3/rCT82u93HLxMkpCAd63qoorm5ac1Ht+//bLeMedT9EOxJ7eSdP8EBcsrkBE+y1USuXGlFUomZSPKpTndnazqbWb9S3dXNdUyzWr6olEYxT4/Tgcp5Zsw0MjnHnnHwBYEYILzlrAbVedkYGolVLqiMmqUGZtE4pt+wbY1r6fHW17ePu5S1lb5+PspQsACLgys9nBQj+PfqiB637cSmEfSHziZo1KKZUNsyqBG2PoHYpyqKePaNLNhYsrOLNMaFxUm7X3bD5rGS3LG3E4HFp9opTKqVmRwI0xREYTbGk7SO/gCJefVYfbleqRsLr8mJaNGTfW+6FSSuWSbRP4wEgMj8uJ2yGEIzHio3HOnBMgEKjOd2hKqdNUa1c/ZQVuigIFJBKJrJ/c2SaBd/f00TM8yp6OfSydP4f+4SiN8+fgdDoIFfryHZ5SStFQeaQTD6fTSV94mEg0SnEwgN+b+Q46bJPAiwoLqSh10zg31WXs3DzHo5RSUwkFC0gU+LJ2fcw2CdybhW8vpZTKNmcW79K2TQK3s3giSTgSJ/1LOGlSNyBN54s5aQyOCWY0JrX8vr4Iy3TYOKVOO5rAs6y1a4DNmzZy2x8jjKaVn1lbjN/tZENbz5TrWD6naMIh4c6tL2FjWy/lAQ9vX17NmbXFXH1WDUGvi4GROA4HBH36y0Wp2UoTeJate3onT74SobbEQ1tvbNy0gcjoJEsdbeK7ZcOROAAHB2P8eP2bAHz5l68yt8TPocEoXpeTL1zRyHVN83CexF2nw7E4XpfzpJZVSmWfdqGXZdesqGI4wTHJG+D1/WGaG0qnsZaJE2jv8Ph1hgrcvGVeiKoiL4PRBIeGYtz22Fbete65E/iySNnc3kPT137P/c+1ntBySqnc0TPwLKsMelhUBG1DgsvhoKG8kEKvE7/biVNgV9fgNNZiOLsuNL7Eqv/uG4oRtfo3b6wKsr61h4qgl+aGUtoPDnEgHGVrRz/v+e7zPHzzaqqKptfk8qEX3mQ4lmAkNnlPj0qp/NIEnmUOt4+dA+B2GEZGE2zrTNVl+9wOvC4HAyMxmutL2NDWi8/loLbEf/jC5lg/Y0U+F3/Zfex40R6Xg2K/h+7BKAAjo6lk2x2O0h2OsqK2iPllhWxo62FX1yDvXvccz/z9WyccgKKle5DBSJyEMbzROcDjWzoAeOzFvbx3RTFzqiozvWuUUqdIE3iWjSXX0eT48saqIC/v7QdgQ1svTfNL8LgcPN9y6Jh1+D0T380Viyfxe44k40QySYHHybB11vxqxwDL01qn7OuP0H5oiMVVwXHreeD5Nv7xiddYWh1kT88wK+YUk7S+PDoODVNeVnZiG62UyglN4FnW2nVs6xGAeHL8hclN7b0srCiccN7j9fhbV1pAZfBItUhzfSnP7Og+/Hpv7/C4+bvDURZXBQlHRln3dAuR0QT/83zb4Zgaq4McsurW/7MZ3vvuq7STLqVmqGkncBFxApuADmPMNSLSADwClJEa6Pgjxphjr9Sd5nZ2DU1YfiI5UYAlVQEk7WJmgceBx+Vkx4FBusLRw+UN5YU0lBcQTxoEiMUN/SNxaxknq+pKMMbw4pt93PPHFtK/R4wxBL0uvJLgZ7eupSSkbcuVmslO5Az8c8B2YOxT/Q3gv40xj4jI3cBNwPcyHJ/tLa8J0lxfSqopoOBxOYjFk7icWOXgdTuIxpP43Q5auidO+DsOTHyx85z5JeMSeOvBocPtw4NeJwsrg5QFPERiCS5bWs7D69t58tVORIRz60sPN1Dc0zNMaaGHodEEg/1hdneHOUcTuFIz2rQSuIjMBa4G/gX4W0n9pr4U+JA1ywPAHWgCP0ZZ0DvuZp1FFYXsOipJL6sJsq0zTE2xl1XzQjgcHDkzNhDwuTi7LsT+gQj7+sYPGnH0iEoLKwrZ3J664DmvtJAte/qAVD36PX9ODea8cm4xL+/tG7fcRYvLMdFBzqjw8Z6rz6WqeOLqHKXUzDHdM/BvAf8AjF39KgP6jDFx6/VeIHujJtjYgfD4WqWJqrPHyjr7o3T2R1ndUMr61mPv0Dx/YSnJoy6GDkUTLKooZDRpSCQNVUW+tLP49PqRI08LJrgoOhgZZZUjwm3vu1TrvJWyiSlv5BGRa4AuY8zmk3kDEblFRDaJyKbu7u6pF5hl3t80b8p5pjssaTxh2D8QGfd440CYsoCX9kPD7O0doW/4yA07E/WfUh7wMBCJc96CMi5prGBFbZDLGsto3dPPJz9wgSZvpWxkOndiXgC8U0TaSF20vBS4CwiJyNgZ/FygY6KFjTH3GGOajDFNFRUVGQjZXrwuJ28740gb6pICz0mva7I8Pxg9krS3dQ5wTl1qFKKx9t6r6kLMCflori+hPODltX0D9A7HeOnNXkr8HlyJKP/6/iWUl4ZOOjalVO5NmcCNMbcZY+YaY+qBDwJ/MMbcADwNvM+a7Ubg8axFaXP/9M7lzAmkdvWm9l7OrT+5Yd4mPzseXz4Ui+MUKPK7ObsuxLZ9A7R0D7G5vZeByChBr4u9vcP0j8TZ1HaIKucwC0LaolQpuzmVT+0XgUdE5GvAS8C9mQlp9qktKeDiZXP5yYZUh1Mb23ppbihhQ2vqYmN1sY9i/5FeA4M+1zF9pLgdgsspE/adUuAe/z38+v4wFy8p57ldB3GKEIunKs6XzSmidyiGOCA8krrZZzQJJcEClixsyNwGK6Vy4oQSuDHmGeAZ6/luoDnzIc1O7zir5nACB9jQ2ktzQyl7eoZp6RpkT+/I4WmrG0rZcNRFzOb6EloOjhw+1+7sP9IaZUlVgJriIzfzzC/zs7NrEKdDiFi3gDodwr6+EeaXFdJ3IEyowI2x+hlvmaSJolJqZtPeCHPk/EXlXN9cN65sQ2sPS6oC45L35IT9/RE6+yO4jureNRyJ02lN6+yPkEjCocEYQqoe3O0UltUUEfC5iYwmcDiEvuFR+kfiDEdHWV5fnrkNVUrljCbwHLphdR3nLxzfr8jLe/vxuk6s5Ud18fgeBdOrxgXYPxAhGk9SG/ITiydZNS/E1o5+qot8bOsM43Y4DtfDv2eJl3evWXJS26OUyi+9cpVDjZWFrLt+FU9v38dXf7WDwWicoM81runfZExaG5SNbb0U+VwsqAjgdAhelzAwMspgNMHyOUW8ao3eMzbgw9iSm9t7WFwZYGfXIAORUb5+5Tw+dMlZGd9OpVRuaALPIbfbRbHTydUr53JgKMEjG/ewoLyQPT3TqUIZbyASP3yXJcDaxWUMx5IEPEf+pX0jqZuIhmOpRB5PwmgiVSf+sVUhTd5K2Zwm8BxzOASvx80nL1nEJy9ZxEgswWMv7eW+P7eO6welMuilqsiHz2phcnTvhUc7NBhjW2eYlfOKAfC6Uv2rAIc7swJoOzTMVSXwpfedn+lNU0rlmCbwPPN7nNywej7Xn1vHszu6ebNnGL/byaVnVOJ3O7nhh+vHnWlPJuBL/Su9Lie1JX66w1EKPU7KAh5KCz3sTbtQ+pWPr8nW5iilckgT+AzhcAhvXXrsqDcP3tTM2m88Tf/I8evJx7qa3dDaQ1XQy8pQksa6cnz+An7+YgcNpV4KfV7cDqip0AEalJoNNIHPcEU+N+8/Zy6/3XaAO965jDue2MabPeMHaVgYgOtX13FeQ4hfv7CbneEoh8Lwsy+sBuCaM6sp9rn48N3P0xGBNw8NUlcWyMfmKKUySBO4Ddx+9RnctLaBmmI/iSR84sFNh6e9pQy+e/PFzCkJALVcXBljYCBMUSh0ZJ76VB80F3rhpxG444fPct8Xr87xViilMk0TuA2ICDXFfgAuaaxg+ZwiXts3gA+4/xNrKQkdOZte9ZaVk67nbZfN56ePtRMJFGQ7ZKVUDuiNPDbjdjq45aIFAHzm/JITGvasqXE+JX4XG/cME40nshWiUipHNIHb0KVLK7lyeTWhyjkntFxJcZCPnlfHKBxuYqiUsi+tQrGhoM/N3R85h8QUbcMncv2qCn69fjdBr/7rlbI7PQO3MafjxEfPGQ6HefTT5+nIO0rNArP2NMwYo0lqAgsWaL/fSs0Ws/YMfCx5Hz1qu1JKzRbTGdR4nog8LSLbROQ1EfmcVV4qIr8TkZ3W35MbJyzLxhL5rq6wJnOl1KwynTPwOPB3xphlwBrg0yKyDLgVeMoYsxh4yno9IxljWFgRoG94lJYDfSQS2gJDKWV/0xnUuNMY86L1PAxsB2qBa4EHrNkeAN6VpRhPmYggIpQUeqgtDfJqR1++Q1JKqVN2QnXgIlIPrALWA1XGmE5r0n6gapJlbhGRTSKyqbu7+1RizQif20l/JMEP/rg736EopdQpmXYCF5EA8HPg88aYgfRpJlW5PGEFszHmHmNMkzGmqaKi4pSCzZTFVQHetao232EopdQpmVYzQhFxk0reDxtjHrOKD4hIjTGmU0RqgK5sBZlpY/2KKKWUnU2nFYoA9wLbjTHfTJv0BHCj9fxG4PHMh6eUUmoy0zkDvwD4CLBVRLZYZV8C/g14VERuAtqB67ISoVJKqQlJLttGi0g3qWQ/U5UDB/MdxDTYJU6wT6waZ+bZJVY7xDnfGHPMRcScJvCZTkQ2GWOa8h3HVOwSJ9gnVo0z8+wSq13inMisvZVeKaVmO03gSillU5rAx7sn3wFMk13iBPvEqnFmnl1itUucx9A6cKWUsik9A1dKKZs6LRO4iLzf6ho3KSJNaeX1IjIiIlusx91p084Rka0isktEvi05Gi1islitabdZ8bwhIleklV9ple0SkZz3Eikid4hIR9p+vGqqmPMl3/tqKiLSZh13W0Rkk1WW966cReQ+EekSkVfTyiaMS1K+be3jV0Tk7BkQq22O0eMyxpx2D+AMoBF4BmhKK68HXp1kmQ2kutMV4Engr/Ic6zLgZcALNAAtgNN6tAALAI81z7Ic7987gC9MUD5hzHk8DvK+r6YRYxtQflTZvwO3Ws9vBb6Rh7guAs5O/7xMFhdwlfWZEesztH4GxGqLY3Sqx2l5Bm6M2W6MeWO681t9vRQZY14wqf/yg+So+9zjxHot8IgxJmqMaQV2Ac3WY5cxZrcxJgY8Ys07E0wWc77M5H11PHnvytkY80eg56jiyeK6FnjQpLwAhKzPVE5MEutkZtoxelynZQKfQoOIvCQiz4rIWqusFtibNs9eqyyfaoE9aa/HYpqsPNf+xvq5fF/aT/yZEtuYmRbPRAzwWxHZLCK3WGXT6so5DyaLa6buZzsco8c1awc1FpHfA9UTTLrdGDNZx1udQJ0x5pCInAP8UkSWZy1Iy0nGmlfHixn4HnAnqeRzJ/BfwMdzF92scqExpkNEKoHficjr6RONMUZEZlxTspkaV5pZcYzO2gRujHnbSSwTBaLW880i0gIsATqAuWmzzrXKMuJkYrXef17a6/SYJivPmOnGLCI/AH5lvTxezPkw0+I5hjGmw/rbJSK/IPVzfqZ25TxZXDNuPxtjDow9n+HH6HFpFUoaEakQEaf1fAGwGNht/SwcEJE1VuuTj5L/7nOfAD4oIl4RaSAV6wZgI7BYRBpExAN80Jo3Z46q33w3MHb1f7KY8yXv++p4RKRQRIJjz4G3k9qXM7Ur58niegL4qNUaZQ3Qn1bVkhc2OkaPL99XUfPxIPUP20vqbPsA8Bur/L3Aa8AW4EXgHWnLNJH6J7cA38G6CSpfsVrTbrfieYO0VjGkrvrvsKbdnof9+yNgK/AKqQ9EzVQx5/FYyOu+miK2BaRaRLxsHZe3W+VlpAYS3wn8HijNQ2w/IVXlOGodnzdNFhep1ifrrH28lbTWVHmM1TbH6PEeeiemUkrZlFahKKWUTWkCV0opm9IErpRSNqUJXCmlbEoTuFJK2ZQmcKWUsilN4EopZVOawJVSyqb+H/+CtLm+B5yLAAAAAElFTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"usa.plot()"
]
},
{
"cell_type": "markdown",
"id": "c313f7bd-9748-4d8c-9e61-8bc112229f3e",
"metadata": {},
"source": [
"We can quickly filter for Contiguous United States, which represents the 48 joining states and the District of Columbia (DC) by the `cx` method.\n",
"This filters out any geometries using a bounding box."
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "ea92949a-7767-4d01-8541-e41ea2f71648",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Int64Index: 49 entries, 0 to 50\n",
"Data columns (total 2 columns):\n",
" # Column Non-Null Count Dtype \n",
"--- ------ -------------- ----- \n",
" 0 State 49 non-null object \n",
" 1 geometry 49 non-null geometry\n",
"dtypes: geometry(1), object(1)\n",
"memory usage: 1.1+ KB\n"
]
}
],
"source": [
"us_cont = usa.cx[-150:-50, 20:50]\n",
"us_cont.info()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "061cd916-390a-451d-ac46-5e273569c972",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAADTCAYAAACY5kjLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABF70lEQVR4nO2dd3hU1dbG3x2SEEKA0HsHEaUJQaQjTUAFvEjxWrCBglxQqoiCnwUEsSB6UUSsXAVFEMECogh6BWkRQUR6h1ClhNTZ3x/vnDslM5l22iT79zznOTN7Ttln5sw6a6+9ipBSQqFQKBTRR4zVHVAoFApFeCgBrlAoFFGKEuAKhUIRpSgBrlAoFFGKEuAKhUIRpSgBrlAoFFFKbDAbCSEOALgIIBdAjpQyRQhRBsBCALUAHAAwQEp5Lr/jlCtXTtaqVSuC7ioUCkXhY/PmzaellOW924MS4E5ulFKednv/OIDVUsoXhBCPO99PyO8AtWrVwqZNm0I4pUKhUCiEEAd9tUdiQukD4H3n6/cB9I3gWAqFQqEIkWAFuASwUgixWQgx1NlWUUp53Pn6BICKuvdOoVAoFH4J1oTSTkp5VAhRAcAqIcSf7h9KKaUQwmdMvlPgDwWAGjVqRNRZhUKhULgISgOXUh51rtMALAFwPYCTQojKAOBcp/nZd66UMkVKmVK+fB4bvEKhUCjCJKAAF0IUF0KU0F4D6A5gO4BlAAY7NxsM4AujOqlQKBSKvARjQqkIYIkQQtv+P1LKb4QQGwEsEkI8AOAggAHGdVOhMJHz54GPPgIaNwY6drS6NwqFXwIKcCnlPgBNfbSfAdDFiE4pFJayfz+wdy+wZAnw8cdAhQpW90ih8ImKxFQofNGxI/Dii8D331vdE4XCL6EE8tib338Hvv0WGDcOGD8eSEgAhACk5Nqb8+eB5OS87f6294e2/TffANdeC1Sv7vosNRVITwfatAl83LS0/DU998IbofQvWLT+ea/9bReoLT0d+PproF+//LcN9rq07Xxtk5UFxMXl36/t23mP3HFH4L5fvAgMGQI0aAAsXgwMGuS/XwqFhRQcAd64MfDBBxSWM2YAO3cCV19t3vn37QMefhho2dLV9s47wJEjwJQp5vXDLpw9C2zbBjz9tNU9IQ4HEB8fen8yM4HsbCA21pgHp0IRAQVDgOfm8g/mztixwPLl5vUhKwsoVsy889mdmBhP7dpqYpzWwkuXgKSk4PcbOhQYNoz3V7lyXN97L1+//z5HcomJwMCBQJUqRvRcofBLwRDgRYrkbVuxwtw+5ORwGO9NYdXa7CbAAT7oQxHeAHDVVcC8ea73584BTzxBbf7OO2k2O3cOmDaNk58ffACUKaNvvxUKP0SPABcCSEkB3nsPuOaavILR4eDw+MMP+UeygpiYvO8dDmv6YjV2FOB6ULo0MGeOZ1vx4sBLL9HmP3Ik78169fQxnV2+DGzcCFx3HVCqVOTHUxQookOAZ2YCrVpxkvLdd2keadWKf6YyZYCSJYG+fYFNm4BVq1wC/MoV88wavoSVEuBW98I84uOBPn2Ao0eBhg2p7Q8bBrRuDdx9N003s2YBP/5IYV+zJlCiBJCRQROMN3v2UPMvUQJo0QJ49VWONB98EKhUyfTLU9iT6BDg8fFcSpUCHn2UWs6uXfQW2LsXWLaMbSkpnvslJgInTgAVTciz5XDktcMLQdNKYcSOArxVK7oFdu5s3DmGD3e97toVmDgRuP9+Ct0RI4AxYyjIDx2iOScpiXb0ihWB336jfX3tWnokTZniUkB69ADOnAHefpuKydChQNWqxl2HIjqQUpq2tGjRQoZNr15Svvde/ttcvixl06ZSUnRI2bWrlLm54Z8zFPr2lfLAAc+2gwelbNBAysOHzemDncjNlbJjR6t74cnq1VK2bWt1L/Jy8qSUaWlSVq8u5bJlUv79d/7bnz8v5YsvSvnUU1IeOmROHxWWAmCT9CFToyeQZ/lyYOlS4KDPvOZk505qMRorVuS1S5tJjRrAG28AHToUPlOKld+7Pzp35v2Rnm51TzypUIGmwcGDgVtvpUkwP0qVohlx3Djg00+Bp54CDh82p68KW2HDf5kfhOBM/9Spvj93OPKaUPbuNb5fgejShd4pW7ZY3RMFQAE5ZIjVvfDE4WDI/v/9X2j7lSgBjB7NwLWFC4HJk2mDVxQaokeAAwzMOXLE92cDnLm03ANpihc3vk/BcOoUtaXChF1HHEWLcrLbTkyaBNx1V/ijlhIlqJGPGQMsWEBvrJMnde2iwo2sLAbtPfCAq01KBq+ZTHRMYrpTty4167p1Pdu1IIqNG83vk4Y/oXXkCNCrF9C2LfCf/9ADoaCTk2NPM0qLFvRUsoqffwYmTOCEtzZbc/iwPq6vpUpRGz93DnjzTXrCDB+u/NL1REpOLk+YwIjvEyeY9GzdOgZ3xcVxJBVqvEGYRJ8AHz+eduVp0zzbmzXzfL9jB23QZpFfwE5iIrBmDf2Eu3WjBhjIzhntOBz2FODt2tHrw2xOnKCHSUwM53LKlWP7qFHUvvWkdGl6v6SlAbNn8/57+GFq6orImDcPuOkmoHZtauBz5tB99OGHKQOOHwceewyYOdMcv31fM5tGLRF5oWg4HC4vk88/d7VnZUlZq5aUCxdKee5c5OcJldtuk3LPnsDbffONlHXqSHnxovF9spIrV6Ts0sXqXuTlyhUpY2LMO19mppQPPCBlq1ZSrl+f97PmzY3vw+HDUk6aJOWsWVKmpxt/voLK2rVSzpkTeLvffpNy3jxdT42o90LR0DTd5s1pltCIi+MwdMAA31kGzehXMHbfm24CuncHbrnFvnbigkxCgvHnyMqi9jtmDAN5OnUC1q+nH7o7o0YxfsFoqlUDnnuO/5dQJ0oV5PBhesI99FDgbUuXNr4/TqJPgAPAfffRq8MOXiYawQpwgCYgh4NmlG7dGIhUELFbII+GkaadY8cYfNOgAYfZmzf7NpFkZHC+ZvDgvJ8ZRb16jOYsrMFlkTBrFh9+NsttFJ0CfM8erq+5xtp+uBOKAI+JYbTd+vXcLzXV0K5Zgh3t3xq5ucYJsX37gP79OZE4YoT/7X74AfjjD6aHMJPBg5mM6/33aatdvNjc80cjUnIewYzRW4jY+F+WD5rHyY4d9tHyLl0K3Y2oUSPu4x5+rTCewYOB224z7vi+cpt407MnzSz/+hewerV55rSrrmIY/tVX05yn/MYD88gj/N5CwSS5FJ0CXPNxbdSIml5cHKPsvvjCuj49/TTw5JOh7ZOezmvRPBIKGnZ5uHrz7rueEbtWkZTE7JkLFrCEW7t2wO23A599ZqxAr1eP9vh27VTJuEDs3ElzWCieQiYmsYs+N0KAw0/NFvXdd3Qp7NyZQ9dy5WiH1Itt21w+ukLwxylShMKpSBE+POLj6Zq2b19ox963j+lCFeZTubL+x8zJof91KMnTWrXynNzcto2ualOnso/t2jEVQ9u2+vc3KYkJt+69l4mxqlZVo0Fv1q2zdaGW6BTgDgeHgDt38n2XLlyPHRvZH/PSJeZPmTmT9i7N9tWoET/PzWVbbi4FuWZL1eypUvIP2KRJcOe7+moK8JycvJkMox0h+F106hT+MdzrVfqqpen9mXeb1g9fxwz1YRuIL75gWHv79nljFEKhSRPgtdd4j7/1FpWVlStp6rj9dv+pJMJFs9NnZuZvsy+s/P03FbRQyMpi9LUJRKfUiIkB/vyTX5T7lxsfnzdCE6CN8cAB/hjt2zPcPiuLrlWXL1Nj+vxzCtGWLYFFi+hBECoHDzLt5w03cJgeiOeeoztkTAw9aooVKzhluYRgndIffrC6J3lZvZoJoPTi7Fm6DO7apd+DOCaG+cSHDeP7n36i6+mRI6z6ozeHD/Oc+/YBderof/xo5PJl4JdfaOIKhWPHTPNWiU4Bnh9t2wLXX0+TxvXXU0uvUIGJruLiOPN+4gSjpP75T36Wk0O7Y7VqkZ27Zk2er3FjniNQ4v3Ll6lp/fe/NP9kZHBitiAIcTvmAwf4+zz8sL4pFz75hA98Iz1vtAjSChWA3buB+vX1PX69eiwaMWcOR6IjR1oTT2EX9u4Fpk/naCpUE0rDhvz+zpzhg93I/7Ov6B5fC4AiALYCWO58/x6A/QBSnUuzQMfQJRLTFZrEvMi+yM5mvuQnnuBrM/n7bykbNw5u2z/+kLJlS0Zmbt0qZe/eUj76qKHdM5UOHazugW9atJDy00/1O97ff0s5cKCU48bpd0x/VK4s5b59rvf33Sdl9+5SDh/OKFM9OHlSyiefZNRhVpY+x4wmsrOlvOEGKTdujOw4Q4ZIuWGDLl2CDpGYowDs9GobJ6Vs5lxSI3mQhISWz9lfroHYWNrDn3/efNvyAw8AL7wQ3LYNGwK//sqnfbNm1IC++cbI3pmLzYIe/sfy5TSZ6RUFWbIk7zMzKuRkZQFt2nDkdtddHPUtXMiIz3bteE2R+rhXqAA8+yyPN2aMPv22O1lZ/N/eeSdNoNWqRZ7NtE4d34XOdSQoAS6EqAbgZgDzAm1rCnadFXY4KIzdQ/xDoXZtuixNmKBvv6zCjiYUgKatl1/WN5J31y6GxhvNF19QsPTsycn1KVNo6rjrLmbHW74cOH1an3M1akSvruxsfY5nZ5YupWK4ezcTgf3731SwIqFrV048h+peHALBqqevAhgPwDud2fNCiMkAVgN4XEqZqWPf/LN7N3Mo242HHuIkaSQsWkQNKNTRQ7du1CK88eedoX2mtfny6HDfRsPf/t5t2dmFy0XywgVzztO2LevC7t4NPP64q33bNmrL27cHF0gULNdcw4jRpk31O6ZenD5Nm/2RI8w/okVKXrnier1lCx0bNHu+P6Vi6lTud+IEty1fPvL+paTwwR4Xx/+0AXlvAkoIIcQtANKklJuFEJ3cPpoI4ASAeABzAUwA8IyP/YcCGAoANfRK7zptmj3d7lJTme87EuLjmQb37rupCQTrkpidzYrndkHTZOyMniYeM1O1+nL3GzKE6Yr1FN4A86f/8IO9BHhODicYMzMZJVmhgn+lo2NHKlb5adOXLtEMFRvL0Zme0amaEmNQ4rJgTChtAfQWQhwA8AmAzkKIj6SUx5329UwA7wK43tfOUsq5UsoUKWVKeT2eagCHdY89ps+x9KJXL6BvX328A1JT+Yfp0yf/iK6ffmIwxvPPR35ORWRIybJoVvDLL7TZBvJ6CodatfKvQ2s2Fy4A999Pl8pnnqELsCawhcj7UI6JCZyXe9YseqwBfBhEGggoJR8y8fEuzX/o0MiO6YeAAlxKOVFKWU1KWQvAIADfSynvEkJUBgAhhADQF8B2Q3roi5Ur7VMuDeAw6dQp/cw6MTEcHleuTJvmvHm0ffbsSY28bl2uu3YF7riDN+CZM/qcW0/sagM3gg0bGIBzww30MzeT557TP8BHw12ztZrz51ma8KWXgh8RxMYGFsj169MVFADKlg3NHLZ6NU0vn3zC72r4cFd6D23u4PXXGZRlAJHYIRYIIcoDEKAb4cO69CgYihc3t9pOIJ5+mkmJ9GbWLKYIaNKE3gbVqnGol5DA8P4GDbjdk08CN9+s//kjwc7ZCI0gNpaa8NmzQO/eFKgLFhijFXtz7pzrXjACOwjwQ4foGTN1amj26dtuY8kz74LnAB8IL73EfDQrVrCtTJngKzbdfDPw1VeebXPmcJ2URBm1Y0fwfQ2DkAS4lHINgDXO150N6E9wXH21vVzUunalfWvfPgpzvWjZkoEnycl5a+y5/2Fr1iy4CbGijTJlaNpatgy45x6OFo1kxw7jHxIVKzLpWig5XvTC4WDOoaFDgcmTQ59c3LKFLpa+6NqVGRn/+MM1d1CxYvAC/Ndf87YNGcKHjEn/x+hUk/Ryk9KLBx5gqP7bb+t/7GrVTCuQqjt20NysondvamBdu+qbXM2buXMZUWwkjRrRu8UKXniBuYnefBOoXj30/bds8awe787+/dTq3Sd+S5d2xZm4IyVTcYwYQfv7wIF55VCXLvw9TFSmbOjKEQR2HJ5fuGBtoWK7CUu7htKbybx5HGL37MncK7ffrv85tmwBXnlF/+O607gxzQxa0jiz+P572qf79w//GEWK+G4/e5afecuS8uXpleKNP5lToQKVrJUraT8/d06VVMsXLROh3Xj2WWDQIOvObyeTEmDPh6w7Doc5D5hevZh3xUhPIaO/67JlzZ8kz8ykGSrSh15uLnMMefPoo76dDhITXRW/gLxZMAGaL4cM4cR11arAqlX8jgBThTcQjRr4jz8Gb6Myk9hYhh4rogezHjLx8f41wUjIyjIvHkKPh925c9Tkg6lcdfgwS79Fqpjs2MEHwYABrrYjR4CtW4H33vO9T82artf33st1QgK9Sg4dckWnFi/OEZCFRJ8Af+UVTjwoPLGjucKOfdLo2hW4eJH5yn1pWYEIdZ+//uKfXs/cGHv3ujQ/I5GSGnG4XLpEV8cSJVjOLtKsn6Hw0ktc3AX41Kn0GvP3AD93juvMTKbuHTaMofXuhJoj3CCiS4CfPg18+SWXRx6xujf2wm4mFDv5D/vi/Hm6/blrW0YyciQjAufP1++YCQk0ERjNzz9HliIiJoYPLyvSX4wYAXz6qWfbX38B993nf58zZzi66daN772Ft42wuaHSiyNHOKyyo2Cwuk9Wnz/aqFSJ2R/N4l//0r8OZ82a5timV6wIP0EbQLtyqVLWmT69lZtz54DrrvO/fZkyTJi3bl1kIw8TiC4NfO1aTh7YFSu1YLtp4FICmzZxyOwd6iwlJxGFcA1FNY8AX8NabcLRPa2A9lorIKs9wJo2Da40WN++LKRhFitW6O+FYoYN//BhfreR2tql5GRfJB4l4eCdWtfhoMdYftejuUwOGWLPnEtu2Lt37uTkAOvXG5IQpkBgNw28SBGmx3UvkpuT47Idx8ZSu3E4uGj1Rb3Rto+J4T6+siNqwt/hYCm7li05RM7IoNDIzOSxx4xhfhmAfsVm5l6/9Vb2ye4JvrxZs0YfoTt6NN1sL1/WP+FWfly44HnPLF/uu+yixrx5jLu4dInunzb3pooeAV6kCN17Ll+2Vx4Ud+wmRK2mbFnPiutm0KkT/7T/+AcFxt1301xSty5HA6NHMxDj1KngMz3qQd26tKu+9BIFeZky5p07Eg4dYjh6uOTkMMnX++/TjHL2rLkC/MwZptrt2JHv09MZ3ay9j4ujPClaFEhLY9BV69Z88NtceAPRJMCFoK/1/PnG5B3RA7uZMawmv0yKRlKyJPDdd3nb16zhukcPCpapUzmnYhZr1lCAt2gBvPEGhYfD4fJMiY/nxGR8PNsSE1053jVh4j4K0UYlaWkMKNGb336jG144kcBZWXQZ3LWL+x86xDkss6OKFy1iCL6/Yhvp6XSOuHiR76+9ltWxAJrZbK6URY8AB/hl2jXnh1XCyq7YWXv55hu6lb34orkCPCGBnhj79zOL5KBBHFlqpqPsbC6aOSk729N/3H0eQBMsv//OEG+9BfiePRTA774b/D4ZGbQf//QTCyMMGAAsXsyH0fr11qSE+Pbb/E1liYmeifGSk/l9zpwZFeXkokuAX7miBKVCHxYtYpJ/o7TX/Jg3z1W4w72qTji0b69/hXohaOZ5/XWaFnzhcHCSetMmCmuA2zZqBPTrx9FN//4svuLug20mDgcfhKGYbP7+m+uxY6NCgNtYTfJB9+5MsWpX7Kx1KvLSooW+ftmhMHMm/au1SD+7MH481/Pn+xd827YxFP3oUdrHn3mGy8SJTKdaowZtz//4h3XCG6D2nZzM/h46RDfGjAwu/hRBKbktQM09N5eJqzTBbjOiSwMvXpyTVMnJtFPlZ5/yDiTxFVjiXtNRSg5t69TxrCPp7xze25w4wfBcbYLV37n9HfPAAVY/ya+/+V3rqVPBbatw8frr1L4HDqTHjJkkJjIgrWXL8I+xY4e+oweHA5g921VFxpu0NH5nVaqwKLS7i93OnbQZP/ggBd+WLazwvnFj3qhV79Jnvt7nV5cVoO3aVzpXd06dYjWhadNc21etysnU3Fwe/8gRHr9sWX72+eeurIc9e9KMtWgRFxvaw6NLgAMcovXsqX/5qpwcRl798IO+xw2WTp1ck2zh7q8IjeRkToj/+9+0h5vNwoXhe6N8/DEjDIsWBd55x9Wu2ck1Ieidg6V2bRYJ8cbh4PxSv37AjTdykvWhhyikMzKoWV+8yFKGvhI27d3LFBfjxvF98+YU6kYxfjwF8alT/sPaL1yg55E2qvDFnDnUzPv3pynqmWc46amtd+5kitihQ+mO6s+kZBHRJ8CvuYbeKHpjdfpTGz7dCwXly1NLtIJp0wJrkb44fZrmioEDGTGoFeF195nXNNjcXM/R4qRJrmIg3lquELTPJyRQu3/8cU48XrlCAZaf/7TZzJhBLf/HH10h794cPx64UpF2rdpDacoUml60IC93V9O0tPBykhtI9AnwlBRXzl49Z7Wttl8rF0RraNGCQu3ECeMr2+TkAJs308f4449pHw4nKVJODt3dpk8Pfd/vvw9upHfttbTT25kuXajxr1zp25R06lT+Oczr12eiqv37Pff3jtBdt4753G0mvIFom8TU6N/fd9L1aEZp4NbQrRuH/k89Zfy52rWjH7jDQR/rKVPCO47DEb7CYdR9ZoUCNHEiI1yrVMmbsAqgANfmNnJzgdRUvj55kmaiPXvoaZKUxCynEyZwG3eNvkkTavqRmDcNJPo0cIA17PTWlqx2T1QauHU0bOiqSm4kp04xN7UeE492fOBb0ad9+5iUbMkSTrDOm0fN2uFgVOWPP7L95Zd971+mDCdeN25k+TaAvuvz5zNtx9GjttS8NaJTgGsVpPUkEq1GDyK9+c+d46QY4HoYZGe7jutw8H1ODu1+GkK4gkfc8X4fE+OKEszOznv+mBhOdmk5r8+e5cTWvHmBH47Jyda6mz3/PAWr0UyfzohAPSbg7fbAt6I/L7/MlMALFjCB2c6dtNUnJwNt2nB+4ddfgauuohbdrBlHPu6ZCM+epca9eLGrrUQJV+RmtWr2fFg6iU4BbvMMYZbw4IMc/rlH6sXGenohxMYCS5dyYueqq/Ju5z6ppb3XtnE4KKBzclzh3N43dlwcl9xcbt+7tyvIw98fXEprgz0AzweakWzZ4ukqGgnhChW7Cf5IeO01z4o4DRvSi6xDB0/Pk7/+oldMr16sUepOWhoF/o03mtJlvYk+SVi5skso6ImWzS5aCTY/zPHjnEMwu0BtfvjKW2Imly4ZU/LMm2+/5SSmlRilTZqtpe7eTe8Yb0cG9wfUAw/QDPLkk9zeXXinp9ODJ8qJPgFeogTwxRf6H9dqN0KFdSQm8g9dowYnreLiOMqIj3fVs9TS3mZnM1GTFgiSleXKYaIpAO45z5s35zJ7NiMT9cBqjylfPPywubEIWpIs99H4hQue27zzjqeP/M8/07RSgAhagAshigDYBOColPIWIURtAJ8AKAtgM4C7pZRZxnTTybZtwP33A9dfb+hpCjRaIQWFi5gYBvLcdx+DNjIzaS66coVLbi630UxExYpRsMfE8HVSEs0wvkx7nTrRbDV5MjVCPYhkpGjUb5+b6z+C0wh69eLcxccfM+ITYBItjV276G5scpV4swlFAx8FYCeAks730wG8IqX8RAjxJoAHAMzRuX+eNGzofzY52jFTqNpRg7OagQNZfOK//9W3cs5339HrxF+wSThkZoZfHNmoUaavYhxGEhND18FevVwpYHv1AjZsYA56bY6ngBPUP1kIUQ3AzQDmOd8LAJ0BfObc5H0AfQ3onydxcZxBnjfP8FMVWJSZyD87d9L7IEvHgWRODs0y7ilLIyUjI/yJfCMUhdRUa9I8V6tG+/ZNN7na0tPN74eFBHsXvApgPIASzvdlAZyXUmqP3SMAqurbNT+MGsUw2g8+AO65R7/jFibTgt008CtXaGrwTnoE+G6LBPfj+UqWdOYMU7T+8os+31NCgjFBZ+F+J0Y8wN96i54fZv9/7rmHvt7uqRCi1JskXAIKcCHELQDSpJSbhRCdQj2BEGIogKEAUEMvLWT8eObrPXiQlbn1oDAJcLtRrJi9It0mTqTGnJhIF7RIyc1luLaeGQ/tNJJaupRh9xs2GH+ukSP5cHU4+H368mB66CHj+2ETgtHA2wLoLYToBSABtIHPApAshIh1auHVABz1tbOUci6AuQCQkpKi3103fDjwn//oVyQ2K8scVzJ/mPWHtNMf365Mm0aPhg8+0Od4kyYx54ZeRZQj8ZjSW0lJT+ckbrNmwHPP5fVE0R5cWsrmgwf5cPzzTyblqlWLphBttOM9OnLn1ClmIz1/nvlPvB+IJ09yrddkcRQQUIBLKScCmAgATg18rJTyTiHEpwBuBz1RBgMwwLcvH2rXduU20AOrIzEL+ySm1d+/Ow4H8z9rdRIj5fbbGQF4yy2sih4pkXxPej7A09IYTzBjBicSfaWPbd+eo6vUVOCRR1gsIS2Ndvz165kf5uTJ/G3XJ04w/gOgK+amTXmFd26uK71GJDnWo4xI/jETAIwWQuwBbeLvBNheX9at03em2exZdG/M0oztGKwkhL36FRPDSjla5XI9ePZZRsDedFPkk6S7d1t7v86eDTRtSuH7/vv5V62PiWGukubNKYSlZDBZQgLw5puc/8jIYE4TXwjhEt4AcPgwtfo+fVxtWnQwwOMVIkIS4FLKNVLKW5yv90kpr5dS1pNS9pdSZhrTRT80b86bRy/spAEajd1s/XbrD0C/8LQ0BoLoJSxfeokBL+7CJxwef5zavBWMGsVJyzfeoBmkefP8t//zT+boyckBPnM6rR0/Tpfghg1pcqlRg+USvYtqeN8XkyfThRJgkiohWEFLc6ncu9e8tAg2IfoiMTW0ENqsrPByKnsTG2svLbCwYcfvftUqmggSElzBIpFy220MQImEffuYFz9U9u+nKadhQ2rCK1cG796YkcEsoJ984rI1B8Nvv+XNHDp8OHP3DB8O3HADJ0Fr1mR1nCZNPN0C3RkzhuaoRYv4YGja1PUw3LOHwryQEb0CHKAdTC8BnpREzeDnn4G2bSM/nl2x40jDjho4wMm1m2/Wt3AIQGWhQwdqk088EbpGPm0ac9+sWxfafkWLAiVL0la9di1w111cB0P//iwGEepv5S68s7JoRmrShMIb4PfQoQNf33kn/3v79+ctqgBQ687IoB28enVXJaKSJe1VLchEbPZPDpHZs/kn08PuFRND4f3YY/qaZuxGdrbt6vrZGq0kmZ4sWMD77IsvXGlLQ6F7d3pwBENODvNcL1nC+pAaHTpQ8LVtyxQCDofLX333bk/vrrQ0av0rVlDzDYchQ3i+UaP85zIqVw74/Xd+576KdPfu7arzuWSJyxvHphXjzSC6NfAmTaiFJCbq48ddsiRnxm+5hUMyI2pv+sMsLTQ3157peO1oQgEoAMMNW/dH3bpc1q71zE0dLI88QhNEMOzaRXfbOnWYfrVhQ9dnmkfMhAnUaGvVouAsUoT3Sf/+tHH//jtHuaEmq1qyhB4hlSpRe05Pz3/0d/68K3eJlu/o3nuB995zbbN2Lb2Ddu/m+wsX+L8tpNjwnxwijRtTE69enZMYkWqXMTFMOzl0KCO99PIFtgsOh7X+7r6wqwkFoIZnlMnp2WdZLSYUfvqJSd2CScHrcNBE0749Jx39MX26q77miRPM3bJrF00sXbuy2MHXX4fWTwAYPZr/x1q1OHkb6Ht092bRlAx34V2pEvs3eDBt6wAjMQuxAI9uE4rGiBHAo4/qOwM9dy61lW7drHcxVFiHw2HMiCUjg1WUAlVNd2fNGk7wff554G3T05m5s2zZ/IW3N5UqUdA2bEivkXfeYXmx+vWDPwbgcjVcuZLzCKNHB97HPRrXV2ZDrQ7A4MEuxapixdD6VcAoGAIc4Aw1QPvimTN8f/Ysb0CHw1XvLhQmTuSfoHXrglNEuVw5BlPYDbuaUIxKsRAfD1y+nDeHdX6MHAlcfTX9r32RmkpzYlYWlZpSpaiIhMuqVUCLFq5JxmC5cIFFgufOpZdLsMVGpPSdM93bzNSnj8tHv4Cniw1EwRHg2p9s+HDa7F5+mdrHgw/SZBBuyP0dd/BmbNMGOHJEv/56Y1YgT1KSfhGGemG3QB53jDKhxMRQs23Vim5xwbBjB33JU1N5L2rf2bvvsrTYwIFA1ap0zfvuO2DWrMhGD717A9u3h668TJrESctwijf7Mllu3er5vlYt2uivvTb04xcwCo4AB6hxjxzpykhWt64rB0MkArJdO9oBe/TQN3zfCnJzlQ08FIy0gbduzQRQd95JjTK/CE2tIlDv3tRIu3Shdty8OYX1DTdQq9+9G/jwQ31C9itVohY/blzw+8yezbmocBWm4sVdoxJ/Zf8OHuR6+/bwzlGAKFgC/OJFTgx99BHf792rX27n+vVpo7v3Xg4toxm7+YED9tXAAWO/r5Il6drZowe1cX/88guF6Ycf0vNq1y4W8G3VimbDY8dcmuq119JDSw9SU6nhb9oUeFuHA5g6NfKHR4kSNPOtXu1qc890WL++ve8XE7HhPzkCtKAB98rf//mPfscvV44BBk8+yVJO0YgdsxHaWQM3S1BMnMhz+dMqt27l5F2fPi4beHIyoxcB3ptGTLbGx9Nvu2VLavz5sX0755/0eOD9/rvne3dzSbdu9r5nTKRgCXAtIrN6daBvX76+7z6aDfQiMZEBP/Pm0R6pFyobodU98I9Zv82QIcDbb3u2ZWXRm+Ojj8KzKetB27Y026Smun6njAyu09OpETdtylGE9kAJF83cqY2iNQYOdL1+7bXIzlGAiH4/cF/UqMGCpkZpm7GxHN7168coNc2HNhpQGnhomPl9LV8OdO5MITllCrXto0epfeo5kgyF55+niebNN1n5pmFD9mv+fE6YxsezFuW0aVRuIiG/ILMVK7iuX99+czgWYkNVTCfKlaO90EgWL+bE6dChxp5Hb+wsMO2GkZOY3kydynu2e3cK8SefZCrWXbvMOb83bdvSA2X7droSjhnDe75iRSou69ZRqM+aFbnwBiiYpfTt+jhsGNfffx/5eQoQQpqoYaSkpMhNwUyG6MX48XRLeuwxhgsbeZ5Dh5ipLVzq1mVeF8BTwPqq2+j+m3m3ay55/oTOX3/Rd7Z8ed/Hcz+mlK7ITX/n9NXH/Lbdsydv4qFjxygUYmP998X7Gr1f+/rcG3+fae3u16G1/fUXA1ratPF9TL1p1Qr49Vdm/mvYkImnHA7zXeZuvx0oUyYyP/JIGTaMmj9A54GXXy60ft9CiM1SyjwpKAumCUVjxgwu/fsbm4VvxgyWk7r5ZuDLL8M7T7VqzBdhNI89RjfL3r2NP5cvOnWyV/3LQAwerH8ulPzYsIFCqqqzRrh77hIzOHuWLomNG1srvAFPm/9779FJYdo0y7pjRwq2AAeYnOrIEePNBk8+Cfz738wd8d13oQtxZdawL2ZNsKalcRKweHHfJomrr+bDpEoVuh46HJ4jlkuXGKhVrBgfAqVK8TgXL3IyNDmZD4YSJXiOhAQer3RpnvfHH2nTvvVWPuQ/+4zn0JwAjh0DDhxw5eyuVEmfVM7+mDIFePppRlE/8QTXSoB7UPAF+Jo1jMg0Q0AOH06/3vbt6aNr5M0drRw6lDernbsJIz/zCMCsevPnG97N/5GezkRMpUrl7ZO36QXIa4pxx9c1uV9PlSoUluXKMYjF+1yxsfSAysmhoI6N9VQUtFHm2bP8ns+cYVBM6dIU6idPUpk5f57CODOTx7p4kUmnhg7leTIyaMbR+qtNGhYtyvNv2gQsW8ZUAFq6Xe/f0Jdp1peZyvt70sx22vUCwafOLYQUfAF+553M4GYWd91FDadtW07yFLISTwGpXj0yE0qoKU0jpWlTjqoeesiY49eqRWGanExBdvQoBXk4aMK8TBkuerN6Nb1Sjh/nQ8ZIcnMZc9G+Pc19O3awfds2/YKUCgAF1wtFo1o1aisvv2zeOfv0oUeBNovvj8WLOSxMT6c2owiM2W6QrVsHF4UYDj//zFFauXKsc5mUFL7wNpIdO+hKOHo0tXqjhTdAP/AHH2TpuC+/dLU3bWr8uaOIgq+BC0HvkEceoQbhL7+C3nTrxiFnu3ZMQu+ds3jMGGqideuyNJxZXg5WB8xEm62/dWs+jPUmJ4fugqNH01Vwy5bgs/aZQU4ORx5aWHulSsBTT5mXe7tKFXoAafTt66rpmZ+nUSGj4GvgGq++Si8RPaMyA9GhAyt4t2tHu6Q7DRqwkOuiRTTxfPqpef2yYyRmsJj9x01MNCYffEwM8MwznJ/RonvDqc5jFJUqAddcw0nOjRuBb7+la6FZ1KrFeQctMd3Spa4iDkaZs6KQKP4nh0hcHIMAGjc297ytWnGSqmNHehlo1KvnEkb169uzzJmCGGG2iYnhKOzRR4GFCzmZGWpxYyN45hmgUSNgwAB6VX3+uTWT8fXrc7L1++9d6Y8XL2Zf3nnH/P7YlMIjwAGaKa67ztigHl+kpDD5VZcurqoiKSkcEioKN3v3ukZgVo2MvvqKD4+UFJpytm2j8LaC8+eBsWM56tGSVmmKTr9+QLNm1psBbUTAO0YIkSCE+FUI8ZsQYocQ4v+c7e8JIfYLIVKdSzPDe6sHCxbQE+Kee8w9b6NGNJd060avg6Sk6DZlFCYSEz1HT3rhcHAkdu21zOttBRMmMBDt//6PGu7Spdbel3FxTBIXF+eq+zl2LNc33si4DsX/CGbcngmgs5TykhAiDsBPQgitwuk4KeVnxnXPIEaM4FN93jxzh4dancGbbqJNsWhRYyNEFfpQsiSrtRuVDdCq/B5jxtB8E6jSlJmV34sX53yRuxvkm28y0+HXXzO8Xk1g/o+AAlwyWYrmCxfnXGyY0i5E5s2j3fG33ziRZBYNGriE+LlzDFKwKk2oIjgyM+mNoTfa5KgeiaBC5cgRemV99ZXvz/fvZ2Da9u2cA/DOXwNQkJ4+zUyBNWvq06+9ezkqee895j8BeP5581yCO9xqPwWQoGbOhBBFAGwGUA/AG1LKDUKIYQCeF0JMBrAawONSykzjuqozDzzAmy872/xzN2hA18Y2bZTwjgYyMozRQKtV4+9/yy2+P3ePVNy7l6Y/LS+Lr2hHX5GQ3tqsxt9/8/4fOZLbaqX2pKTfd8uWDKL5+uu8+7pz6620TevlK795M9ea8AZcgnv/fs9iLYrgBLiUMhdAMyFEMoAlQohGACYCOAEgHsBcABMAPOO9rxBiKIChAFCjRg19eq0XWVnAn3+6KvmYybXXAlddxUmb5GTzz68InqwszlnoSWwsR36bN7uyUObHokUMX/cudBCIDh1CS5JWqRJtzsFquV9+yeRXPXrwARAby4dMXByva/r04E2Ely8z+G33bnqhaMTH8zPlqZWHkL4RKeV5IcQPAHpIKWc6mzOFEO8CGOtnn7mggEdKSoq9TC+VKjFtp9nh2RqvvQb07Ml6hwr74nDoKzy2bKHGu3t38Jq9lK4MhUZx9iyvU8u9HSxr1zLiOCuLy5UrNDt9+SVjLwJp8Rr+HpLnzyvh7YdgvFDKOzVvCCGKAegG4E8hRGVnmwDQF0D0lYg+fZqz7j/+SE3ICE+D/GjdmiacQ4fMO6cdK/KEwpUrwP338/WFC9HpUpadTZtyKGaZbdvCywmekEABGAxlyjDXSfPmrC4fLPHx3LdSJVbDatCA+UomTaLXVTD+7evWcR0b60qktXQp/dKLFQu+L4WMYMY2lQH8IITYBmAjgFVSyuUAFgghfgfwO4ByAJ4zrpsGMXEia2Zu2wbs2wcMGsSSVl268EY2qw/dutHGqQjMhg0ckmsFMFq3trpHoZGWxshcdxtvMOzfH14Sp2bN6PEULIMHU3B+9FHwmnN+jB5Ne/qgQf63GT+epp5//pMPt5Yt2d6nD8P3Ff6RUpq2tGjRQkYN/fpJOWCAOefq1k3KmTPNOde//iXll1+acy5fdOyo7/Fat5ayXj0pJ06U8o8/pLxyRcqTJ6U8c0a/c3TooN+xACl79Ah9v86deW2hsmKFlMOGhb7fuXNStmkT+n7+GD9eygcfzNu+bh2/k7lz9TtXAQTAJulDpioH5LVr+fSfPZtDv6uvpnbUvTsnjho2ZBTlu+8y1afetG3LAIUxY/Q/tj+s9DvX24Szdi01xv37qeUlJ3MEVakS7bJ2MrFs2cJ1OKHgOTnhpSbu1Cm8mprJyfTJ/uKL0Pf1xfTp9CYZP96zvX175jwZMkSf8xQy1MxAYiJ9skeOpDDdto2z/XXrUthcvMi8xOXKMTvbzp36nXvxYrpzKb/W8ImNpW3444/5XguMGjmSJpb4eG4zcCDXc+awrNyZMwwXz8pi24kTzPvRsiV/fymBgwf1zaudksIUqWamjI0kGdfnn9MeXrcuI4kjZe5c5lipUIETplpiOZVSImyUAE9J4TJ2LCMjvSlRggIe4J/cV+HecElK8nSXUkSONrp49VUWSh4wgD7KM2dSA77xRtrRd+0CWrQAtm6lf/VVV1Ggf/89R2M7dtDNr39/fX7rnTt5HCsy6YU7CklKYjj7zTczIZtmmw4XrRBEcjIdB666ism8CmmhYj1QAlzDl/D2Zv16mlgAV/X3SMJ6d+/mMFWhPzEx9IIA+JC84w6aA26+2eWS9uKLwAcfeHp3uOeL79aNpjSAGmjFipzYCyf9wogRXKfkKSxuPLGxHGmE0+8aNTgi6dKF+evDjeVwOGjO2brVM/I0FG8XRR6UDTw/pGSdQIA21UGDOByVEhg3jkIi3EjO8+c5pHzkEd26GxVYmceiTx9Pf+Jx4/J3zevWjb/1N9/QfXHDBpe7W6hs2hT+taenR+YHXa9eZGUFa9dm+ocOHcI/xowZ9GyxIm1AAUYJcF98+CFtdEJQyAIcWqemsko2wBtyzBhOcobDihWclIo2N7hIiUY/9JtuopvnP/8ZvhCeODH8a//jj8iihVu2DP/Bo9G8ObNovvhi6PsKwWrykZpgFHlQAtwXW7e6PAY0rrmGw8Bn3LIFzJzJ/OLhpKa96SaVVa0wcd114WvRf/wRWQ6QLl14T0fKjh1UYLaHGLOXmMigK4XuKBu4L3wVQL71VmZEq1jRs/3TTymIP/ggtHM4HIzAHDmSIfUK+3PpEj1Ywsld89df4ee82bOHE37hUrdu3pJ+4VC/PtM+tGzJEWTnzvlv/+efHKGWKZN/cW9F2CgBHiwDBnAYe+FC3tSizZszY10ofroVKnAS84Yb9O1nIKw2YUTzqKN4cU66hWMLnjSJvtDhcOAAFYhIOHMmsv01mjRh8YfHHnPVqPTmwAFOGq9fT4Fv9j1eiFACPFjuvpuTOb6E9JYtzNfgSzhevEj3qRMngJMnuT57ltr3l1/61vaNJpqFqNWE65LXqxfwxhvhFfA4diy8PCjulCnD3N+9ekV2HIDBOFOn5m3PyPDMWzJsmBLeBqMEeCi0a+e7/ehRZoqbMiXvZ0lJQOXKXBo2pCtV2bL0cDhyxPyCsXaKTIw2InnwVavGoK0tW0J3JczOjtx744kn6BuvhwAHOAp95hlg8mT60Rct6pnPZ/XqwCYWRcQoAa4HVarQxpeYGPyffNw4JuqpXJnDTW/bupGoEm7hE+4DsGZNTmKOHs3wf7Pp2VPfIKKDB+kTPmqUZ0m2pUuDyz6o0AUlwPUi1ICcHj24HDtG/9iTJ5kXomdPY3MfW20Dj2YiefCdPcuH9J9/6tefULlyhTUwBw6M/FgxMZwTeuUVV1vJkkp4m4wS4FZTpQr9yaWkr+7kyXwYDBzIAAwjUBp4+IT7AFy+nKY2Pdz5wuHYMYas9++v3zE7d+aEJgAsWQL07avfsRVBoQS4XRCC3g0dOtBW+sknDCJq0gT4xz/0i2ALRwBVruxyY/OXC8bddOT+mff2zZqFfn67kJLCB6x7uH2wHD1KO3io15+To8+k85dfUuDq+fBu394V2KOEtyUoAW5HSpVy2Su3baP7mcPBP0nz5pH/ob3/xDNmcHjt77gXLlB4TZ9euEtbPfggg7eWLAFuuy34/ebM4STilSuhnzMnh7UmI+WTTxjHoBfp6UDv3nydlaXfcRUhUYj/jVFCkyZcMjM5QbRwIVCnDic+S5XS5xyXLrlSBLhXQtcqnE+axKouc+YA//qXPueMVu69l6lrgxXgly4x381TT4WXzzshIfx0sO6kpDA97LhxkR8LcOUAWruWBYwVlqCModFC0aK0i8+YQRfE2bOp1a1fH/nEpBDUymNiqO0VKeJ6ra179aLnweHD+lxPtNK9OzXwYDl8mN+v9oC0invv1df7RYsq9edaqzAFpYFHI3XrAk8+SS1o2TIGViQkeJpGfJlDpGTB2IsXXX9mIYLX8J5+mr7ExYszEq8w0rw5ff4/+gi4667A2z/zDDXUSGzPetjA165lPh89qVZNBYVZjBLg0UxcHNCvHxfN3BEIzWsgHJKS+OB46y1mbLz77vCPFc288gpw//3BCfA1a0LT2I3iwAEGkOnBvn1c//Of+hxPETZCmugXnJKSIjdt2mTa+RQGISWwYAFHAPfdZ3VvrCE2NriRS5UqTBlcp45nu8PhWgDmxtEKUHjTsSMr2GisXcv83NqEshA0scXH5/UGys7mJOPlyyxGoUf5svr1mWBLxRSYhhBis5QyTwiv0sAVoSMEcOed9F+fNYsTm4XNtzxYt87Fi+mDreE+UoqNdXmYDBjgX4B7M3YsC42UL+96CGRk+PYGiYujeS0ujqOBY8fCr8m5di3Qpg1T2+7ZE94xFLqiBLgiPISgD/Ds2UDjxgxUqV3b6l6Zw88/B1/LNJiCHe++G5x9Wgh+zw89xMr2a9eG5taZmkrvmebNmdyqQgUu1aszaKxCBf8P4u+/Zzm6okUZp6CwBUqAK8KnSBEWpR04EBgyhDUn9fBZtjsDB9KlUy+KFw88gtHMFe+8Q3fA+fMpyOfP57zG5MlMT1y7Nn8DX/MhU6YAt9/O3CWnTjF9w7ZtwMqVQFoaXR59BWT9/TdjEG69FVi0iGkBhg3T5dIVkRHQBi6ESACwFkBRUOB/JqWcIoSoDeATAGUBbAZwt5QyX49+ZQMvwMyfTy2yMKQPTUrSv0BB48YUzMeP015evTrby5alwHzpJWDjRtrCt21jdff80NM+vWsXo0/btWNAUHY20KpV3qpVCsPwZwMPxnCZCaCzlLIpgGYAegghbgAwHcArUsp6AM4BeEDH/iqijd69w6uXGI2UKaP/MR99lHUzO3ak8O7Vi6aOixcpjEePphnj2LG8wjs5mYFdw4bRRg1Qk9aWSGnQgJr76NF8HxfHfiosJ6AJRVJF19SNOOciAXQGoPkRvQ/gaQBz9O+iIiooV47RfidPmpsa1wouXgyvMEN+FC/OCc9y5ajhNm4MjBjB5FdZWRSa5crR/PHgg0xPmx9ff+3K/b1smSvsXWPbNqBpU74OpmrOq696vk9M1P87UIRMUDZwIUQR0ExSD8AbAPYCOC+l1PyojgCoakgPFdHDlSsUbgVZgB85wmIGegqulSupQQPA6dNcjxsHtGhB/+1w6NmT/tp16jDFq7tJ5dQpl/DW3BxnzAgtzL5SJRY3btIkvP4pdCGou1BKmSulbAagGoDrAVwd7AmEEEOFEJuEEJtOnToVXi8V0cHttwP//ndw254/z2LOlStzaH7ffUy2NGMG7cBWIwTNFS++COTmMvISYPDKW2/R/vvhh8DmzfwccAnJt94Chg9nxGww19Kpk6dJonp1pkg4dy6ya6hd29W3tDTaz6+/nkI9Lo6TnkePcpslS/g+WHr1Ar75JrL+KSImJDVCSnkewA8AWgNIFkJoGnw1AEf97DNXSpkipUwpX758JH1V2J0mTagxBhI8Fy4wejQuDti5k6aC4cPpopaeDnTtysATq6lalWkKcnMZdTp9OvD770wpoAm/Awc4uffCC3TDA+gf/9VX1ICrVHFp1f6Ij/cM1Pnvf2lGefPNyK9BGylUrMgHxcaN9FrJynLlm4+JYaKyBg2CLzhRpw41cIW1SCnzXQCUB5DsfF0MwDoAtwD4FMAgZ/ubAIYHOlaLFi2kooBz5IiUt90m5Vdf+d/mhRek/O9//X9+991Sjh4tZWam/v0LlqFDpaxXT0rq1VKWLCnlvn1S5uRImZWVd/tffuF2bdpI2a+flKVL8/2tt3J9/fV595kyxXV8bblwQcqlS6V88019r2fXLim7d5fy1Cn/2/TpI+Xw4cEfE5AyNTXirikCA2CT9CWffTV6bAA0AbAVwDYA2wFMdrbXAfArgD1OYV400LGUAC9EDB0q5erVFHje9O0beP9Vq6QcP17K3Fz9+xYM5865hOqqVeEdY+tWKR0OKe+4g8dZsIDtu3ZJWbcuvyNvAa4tVpCZyXNXqRLc9lWqWNfXQkbYAlzPRQnwQkRurpQLF0rZvLmU337rEuRbtkh5yy3BHePnn6UcMcK4PvrDXTN+7bXIj+f+MAhlsYI1a3juFi2k3L8//22zs5UANwl/Alz5ACmMISaG+T1+/pm24u7dge++o0154cLgjtGmDUu5bdhgbF+90TI2nj6tTwGL5GSK5DVrPNsTE1kI+O+/+bkdJvm1ydTNmwOnRtDC+AtD8JZNUQJcYSwJCfQweeUVhtq//HJo9T379AE++MC4/nmzbBnX2dn6pV/V6NiRk7Pp6RTYly9TeJcsyc/LlmWCqptvtjbTn5SulLGBKth//LH5D1jF/1ACXGEOTZq4El+FQo0aLCcXTj1Jfxw/ztFBq1ZAw4YUpC1a0HUwNZVpYo2q/ZmYCBQr5vszIejut3y5MecOBU37DuQ5NmgQ16qwgyWoZFYK+zN4MP3Lx4zR53hCAJ9+6nqflMT1qlV0YVS4CKVgcVYWXSIVpqE0cIX9ad+ewlUPs8LFiwweGjXK1fb66zy2Et55uf/+4Lf96Sfj+qHwiRLgiuhg4EBAj0yWf/zB9axZXF+8yKrxCk80zTuYCcpXXuH6xhuN64/CJ0qAK6KDQYOYC3vz5vD2z8mh6UQTSKVKMRmTZj5ReOKd/Co/tNGMFravMA0lwBXRQbFitIPPnk0vjlBx196bN2cuFjXx5p/Dh7kWInDOE+17DKbIs0JXlABXRA8xMZzQXLQotP2E8CxtFq4WX5jQTE0AMxsG4sMPg/fvV+iGEuCK6KJ9eyaK+u674PfZtcv1+uxZ/ftUENm40fU6mGLLmvb99dfG9EfhE+VGqIguYmOB999n/uz09OBstW+/7XpdurRxfStIpLhV75o8Ofj9evWyNgipkKE0cEX0UawYvUh27gy8rRDAzJl8feaMsf0qqATr2+2eEldhCkqAK6KTjAzgt99C28eIWpYFmZMnQ9u+USOuU1N174rCN0qAK6KTBg0CV2YH1HA+EipU4Gjn7ruZziAQ2gNSC69XGI4S4IroxeEA/vrL6l4UbDZuZDm5hITg3C5PnfKcNFYYihLgiujlhRdYCkxhHNde63qtVbnPj2BGRQrdUAJcEb0kJbnyaeeHVmg5GDOAIi/t2nHdqZOl3VDkRQlwRXRz660s/pufrVtLRdu0qTl9KmisXMn1+PHAZ59Z2xeFB0qAK6Kb5s1ZOadhQ2DGDN/5OIYP51rZZsOjWDFg7lygaFGgf//A2z/wgPF9UgBQAlwR7cTEAC++COzYwQCf2FgKcncSEqzpW0FiyBDgwoXgtg1k0lLohhLgioJBkSIU4jk5zEtdvDj9koUA1q8HJkywuofRjxbQ8+qr/re5csVVjk1hOEKa6CebkpIiN+mR01mhCJZz5zwDeJRfeGS4uxI6HHldC2Ni+B2r71lXhBCbpZQp3u1KA1cUbEqXZnGCAweUUNGDK1eAbt34OsaH+FDfsakoAa4o+MTFATVrWt2LgkFCAr1Sxo7N+5k2SaxyophGQAEuhKguhPhBCPGHEGKHEGKUs/1pIcRRIUSqcwnCy1+hUBQIXnyRa/c6onPncn3qlPn9KaQEk042B8AYKeUWIUQJAJuFEKucn70ipZxpXPcUCoWtWb2advAePVyVe/r1s7ZPhYiAAlxKeRzAcefri0KInQCqGt0xhUIRRWjCu0sXa/tRyAjJBi6EqAXgOgAbnE0jhBDbhBDzhRA+M+ULIYYKITYJITadUkMrhaLg4GvCMpRKSYqICVqACyGSACwG8KiU8gKAOQDqAmgGaugv+dpPSjlXSpkipUwpX7585D1WKBTWUqsWzSaVKnm2Kw8U0wlKgAsh4kDhvUBK+TkASClPSilzpZQOAG8DuN64bioUClNZs8bz/Q8/AEeOUHAfPMg2reDDN98o4W0RwXihCADvANgppXzZrb2y22a3Adiuf/cUCoUlNGni+b5zZ6B7d9f7Z591BezcdJO5fVP8j2C8UNoCuBvA70KIVGfbEwDuEEI0AyABHADwkAH9UygUVuBdfk5p2LYkGC+UnwD4KsXxlf7dUSgUCkWwqEhMhUKhiFKUAFcoFIooRQlwhUKhiFKUAFcoFIooRQlwhUKhiFKUAFcoFIooxdSKPEKIUwAOmnbC0CkH4LTVndAZdU32p6BdD1Dwrsnq66kppcyTi8RUAW53hBCbfJUtimbUNdmfgnY9QMG7JrtejzKhKBQKRZSiBLhCoVBEKUqAezLX6g4YgLom+1PQrgcoeNdky+tRNnCFQqGIUpQGrlAoFFFKoRXgQoj+QogdQgiHECLFrb2bEGKzEOJ357qz22ctnO17hBCvOXOl24J8rqesEOIHIcQlIcTrXvvY9noA/9fk/Gyis9+7hBA3ubX3cLbtEUI8bn6vg0cI0VQI8YvzN/hSCFHS7TOf12dnhBDNhBDrhRCpzjKK1zvbhfP+2uMswdjc6r4GixBiofN6UoUQB9xSatvjN5JSFsoFQEMADQCsAZDi1n4dgCrO140AHHX77FcAN4Dpdb8G0NPq6wjieooDaAfgYQCve+1j2+sJcE3XAPgNQFEAtQHsBVDEuewFUAdAvHOba6y+jnyubyOAjs7X9wN4Nr/rs7q/QVzPSu0eAtALwBq3118777MbAGywuq9hXt9LACbb6TcqtBq4lHKnlHKXj/atUspjzrc7ABQTQhR1ViAqKaVcL/kLfgCgr3k9zp98rueyZE73DPd2u18P4P+aAPQB8ImUMlNKuR/AHrCk3/UA9kgp90kpswB84tzWrlwFYK3z9SoA/Zyv/V2f3ZEAtFFEKQDa/6gPgA8kWQ8g2auil+1xjk4HAPjY2WSL36jQCvAg6Qdgi5QyE0BVAEfcPjvibItWovl6qgI47PZe67u/druyA64HTH8A1Z2vo+06NB4F8KIQ4jCAmQAmOtuj9XrcaQ/gpJRyt/O9La4pmJJqUYsQ4jsAlXx8NElK+UWAfa8FMB1A9/y2M5NIrseuFMRrcie/6wPNJq8JIZ4CsAxAlpl9C4cA19MFwGNSysVCiAFgLd2uZvYvHIK8B++AS/u2DQVagEspw7p5hBDVACwBcI+Ucq+z+SiAam6bVXO2mUa41+MHy68HCPuajsKlrQKefffXbglBXF93ABBCXAXgZmdbftdnKfldjxDiAwCjnG8/BTDP+dq21wME/o2EELEA/gGghVuzLa5JmVC8EEIkA1gB4HEp5c9au5TyOIALQogbnPawewBErYYY5dezDMAg59xEbQD1wQnZjQDqCyFqCyHiAQxybmtLhBAVnOsYAE8CeNP5kb/rszvHAHR0vu4MQDM3LANwj9Mb5QYAfzvvv2ihK4A/pZTuJkd7/EZWz+xatQC4DbRbZQI4CeBbZ/uTAC4DSHVbKjg/SwGwHZxxfh3OQCg7LP6ux/nZAQBnAVxybnON3a8niGua5Oz3Lrh5z4AeD385P5tk9TUEuL5Rzr7+BeAF9+/f3/XZeQG9nTaD3hkbALRwtgsAbziv53e4eRRFwwLgPQAP+2i3/DdSkZgKhUIRpSgTikKhUEQpSoArFApFlKIEuEKhUEQpSoArFApFlKIEuEKhUEQpSoArFApFlKIEuEKhUEQpSoArFApFlPL/23yEw93sjJQAAAAASUVORK5CYII=\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"us_cont.plot(facecolor=\"none\", linewidth=0.5, edgecolor=\"red\")"
]
},
{
"cell_type": "markdown",
"id": "e4509aa4-4036-4621-ad0c-2e87eef2340a",
"metadata": {},
"source": [
"We can then transform the `geopandas.GeoDataFrame` to a `dask_geopandas.GeoDataFrame`, with a fixed number of arbitrary partitions - this does not yet spatially partition a Dask GeoDataFrame.\n",
"\n",
"```{note} \n",
"We can also read (and export) parquet files in Dask-GeoPandas containing spatially partitioned data.\n",
"```"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "7dd696d3-6492-4890-96c5-8f62b88e3d7e",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
Dask-GeoPandas GeoDataFrame Structure:
\n",
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
State
\n",
"
geometry
\n",
"
\n",
"
\n",
"
npartitions=4
\n",
"
\n",
"
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
object
\n",
"
geometry
\n",
"
\n",
"
\n",
"
15
\n",
"
...
\n",
"
...
\n",
"
\n",
"
\n",
"
28
\n",
"
...
\n",
"
...
\n",
"
\n",
"
\n",
"
41
\n",
"
...
\n",
"
...
\n",
"
\n",
"
\n",
"
50
\n",
"
...
\n",
"
...
\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Dask Name: from_pandas, 4 tasks
"
],
"text/plain": [
"Dask GeoDataFrame Structure:\n",
" State geometry\n",
"npartitions=4 \n",
"0 object geometry\n",
"15 ... ...\n",
"28 ... ...\n",
"41 ... ...\n",
"50 ... ...\n",
"Dask Name: from_pandas, 4 tasks"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d_gdf = dask_geopandas.from_geopandas(us_cont, npartitions=4)\n",
"d_gdf"
]
},
{
"cell_type": "markdown",
"id": "2ff72689-2f40-4e43-87c4-a847c9f1bf4a",
"metadata": {},
"source": [
"By visualising the convex hull of each partition, we can get a feel for how the Dask-GeoDataFrame has been partitioned using the fixed number.\n",
"A useful spatial partitioning scheme is one that minimises the degree of spatial overlap between partitions. \n",
"By default, the fixed number of partitions does a poor job of spatially partitioning our example data - there is a high degree of overlap between partitions."
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "a5346f4a-cd57-41b2-80a6-f70d11ed5d49",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0 POLYGON ((-81.96430 24.52042, -82.87569 24.610...\n",
"1 POLYGON ((-89.42000 28.92833, -90.92028 29.048...\n",
"2 POLYGON ((-109.04997 31.33190, -114.63287 35.0...\n",
"3 POLYGON ((-97.40186 25.83808, -97.40523 25.838...\n",
"dtype: geometry"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d_gdf.calculate_spatial_partitions() # convex hull\n",
"d_gdf.spatial_partitions"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "c7b09d92-c3cb-4a5d-a453-656336d14a1b",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmwAAAFUCAYAAACZVGB4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACeqUlEQVR4nOz9Z7Bk6XnfCf6f97i015Q3bap9F4BGNwByQIAACRpRxBQpStqhpB3NrgIjzW7E7sZs7JiN3tmI/TK7ETWz1ITMSDMkxemVRqJEASAAAgUQlgAaluhutEVVd5m+5a436U7mce/77Idzzr15fea9aU7mfX8Rp25l5jHvyTzmfx5LzAyNRqPRaDQaTXYRwx6ARqPRaDQajWZvtGDTaDQajUajyThasGk0Go1Go9FkHC3YNBqNRqPRaDKOFmwajUaj0Wg0GUcLNo1Go9FoNJqMowWbRqPRaDQaTcbRgk2j0Wg0Go0m42jBptFoNBqNRpNxtGDTaDQajUajyThasGk0Go1Go9FkHC3YNBqNRqPRaDKOFmwajUaj0Wg0GUcLNo1Go9FoNJqMowWbRqPRaDQaTcbRgk2j0Wg0Go0m42jBptFoNBqNRpNxtGDTaDQajUajyThasGk0Go1Go9FkHC3YNBqNRqPRaDKOOewBaDQaTTtXn75IiB8mjX2mTubZNHmGZX3+sY9/sGnmWh9YeueHH1i6MQugAcC9eO1qOLi91Gg0mu7Qgk2jOQJcffriXuKma+HTw+V3W7YvELMZkfH/tGWIa9MPexdq8/982m/UAeDq0xc9AC5iAbfX1Lx47ars1xg1Go1mJ7Rg02gOQGIFGoRw6dXy1J9vYnSRJHL3SicfmPYbV5O3csl0fL9lrz59sYm9RV0q/JoXr13lPgxfo9EcMbRg02SCLW6wYQmfbpbV8Z+jCNEm4bpQOPboMyvvXt1t9j0oJNOpfeZTV5++uJvVbuv7vhZ3Go1mN8ZGsP29P3mZAOATHzyLb78yh/T/nfCpJ86N5UUycYMdVMAMQzRpNH1l64m+lJ96RpL4isFK9WmTAkA5mfYjuvr0xb0E3fqk4+00mqMHMY+2Vnnh+mwOwM/V3OA3DEGFYt5q3JlvPPXubP1v8w5uIAIYBIXt1+7dEDvOyx0vv2WxTWPabR1bx80AwAwYBlVafnQsjBTiiRFGsu3/iicba2sfmXn5GPwAHEbxguvr4fg76GK83cy/0/Lx7nS33YNua6dtbIxhR7cgbyyy6zwAmDsZP4No9+1s/u+m9RFRIZdDMZ9HIZ+nYj6HQj6Psm15ISO3aU00PO/mYY+HjrfDvPElSulVWy3HbbXgNlvstlpoeh6UUpvOp85+HxCDNj0cfGz2jX9xoT5/v7d70Hd87CPq0s91vJ1GMx6MrIXtheuzEwA+DODnADgTRVshvkChmDfXiChgZmfrcgwQeDStOUSAILIcyxCOtfsulGWuOHX2A+suOxVG4DCECkJwEILDKP5/GEIF8WcchFDJX442X997c3em/t/lcdCxdiKAejN+IkLeSYVZ8jeXQz7nbPXWpSMToOy4X2nTnz5up+27KBgUwpooTE5MbMzAQMv34bZaaDSbaLY8NFoteL7f9bYWCtPnR1CwOcnUSbxdC/snUriI4+36ZWnUaDSHZOQE2wvXZ08C+CiA92OLG00pJqnYOD6ZWzxzPP+t5ar3oTBSx5hHbz97ibBMwDJhFPIdzc/MbQKuTeylQq9d3CXvoW8epdEl7zjroiwVZoVcDiS60TuMo54vIEHbz18C8jkH+ZyDE9NT628rpZBY4uC2vPj/rRbCKNp1/U0rN9mPcWeIfDKd3Gc+7iLeztPxdhrNYBkJl+gL12cJwIMAfhHAU1s/bzTDciFnNpYq3rmqG5w7PuHcrjej4yennDs/e7fyO64XPTHwQfcJ0xDVSKo9bzBl6bUeWbvfmTrrESzlZmtdasELo13EXzgAx9pgcGw7EWU5FHLJ33weoithtjM22AtAuf3nHG9soBXEouNABGHYJuKa2HCrxq76U821V457tXtN0ykf92pzF9du3+jl+McQif1FXRpvFwxrkBrNOJFpwZYItacQC7UHd5uv0QzLq3X/bBCqYqUePFXMm/fOHMtfdWyjtbjWujAz1/jdgQ26z5iGqEVSTew1zzAE20HYLOb2cNUm72111Q4a27I2uTJTgWYY/fNYasEW2xctwAuA3n4PDDQ9D03PQ6PZTKxxHnzfl/nIXwSAnAzWTjdXb3xo8Z2fpovNlM+cvzF1/llFwsxFQR1gPucu33qsOnvnaNtCdyVAB4kUiOPtdjeFajRHnEwKtheuz5qIXZ4fBXBiv/nDSFk/fGPxvz4xlXvRMChwLFEHAENQWG2EDy1XvY/1ecgDY5wEW7fs7aoNwEHUE1etZZpbhFnszjTNwYc+HtayNA6Y4CAC2YPanpQKTS92qzZaLchmo3rh/rWvtMLIaFq5iZuT5/8K71DWpRS27h7z6u9+bPb1bwtk8MI6GrTQefFiHYehOVJkSrClGZ8AfgFAqZtl6244UWsGJ5Yr/tPnThZecVvh8YXV1s8zIDxf7mqdGzU6E2x+65G1e0f6Jp/S7qrdSLqIXbVCSuSEgZwwUDBN5E0TRcuGbWYn5FELtjj7k2LpPdRkoSAI12PiNiYP7dfQE1719XON5bePe7Xl8+7y4hCHO+4wgL2KF7eLvpaOt9OMA5kRbP/sjTulm/drf5q3Te+hM8U3bMs4cJ0hZsa9pebj1XrwsBdEJxutaFvc26hiGlSLJGvB1gVEBNs0YVkWLNOEZZmwTBOGsfP9nyIJkhIk1cb/o/i12PRagiIFKNW3tAAt2AABSAbETmV6hg0zo+X5m0Vcs4XJysLLf+XuS1+6Xzxxci1XPl4MvcYjtbl7wx7vEUViZ6vdTu8FWtxpskpmTAmFnFlwLCO8t+T+9vxq8xMnpnI/PHMs//ZkyV7dqdzBXijFxuyS+/FCzrr74OnSi/cW3bDeDN/Xp6EPmLim2V4c1asNEcWCzDRhWyYs04Jl7S7MdoNNA9yF+5OY1wXdZjEXC7pY/G2IPookKCMPSqOAAgwT8KO4jEWmICIU8jkU8jmcxPT6+0q+5wOm9/7aiqs+sOwFU81WU55+53v/uFBbrg9jnAxgvnDsxL3SqYerTvHMI7W5Nx6rzt4ZxliGgAFgIpn2I+yweLGrixdrBk1mLGx/eO3eQ6++s/q1RivcZA2zLbH46Lnyl05N5+8LQYqZt9WrYmYohjAEKQBwW2Hpnbu1Xz05lXtrdrn5Uccyltfq/n8wwN3pG51Y2ErKbz26Ot4WtliUbbaYmRlyZe4HKbXJSrdN0CUiz4lUK4xknqTMnnlpQCSViJUasXZgBhAqwEwtgwJQxdC7b7rVxfO1havR2lIjXFl0o6W5JkehqtmFwmsnHv/5ueLx5z60+PbneiGoJAnx2onH339z6vwv+8Kcav+sFLbuPlm5++LF1dvXj+qxdUg8dB5vp4sXaw5NZgTb7//s7snXrq9+f7cSHAQoIchjhmWZYsU0qFHImfctUzSbXnQ6jNTEEw9NftnzZXGl5j124UzppfmV1hMgMCsY95bc3x70PvWDoybYNixmFszEcmYaZgadY/3BACKZWMKz5KodNOPkGhaAOiPkrZMkFyxCxMw863oTbzXlM7WWZ7mtFqJGfe2R+Xe+5hl2oRD5DVNFQUSGFRhWrhi2ak9U79/eaxsKRG8ef+TiO9MP/apnWHsW181H/sJj1dnvPbN88y1DJ0v0i73i7doteU3tktXsRmYE2wvXZ803b659a7nqffyg67BNsWxbYkUqOCennJ+2fHms0QofOXO88EMplX1vsXlJMQ8s26wfmAbVI8l79iUsqaD16Ordkbq5maYB22y3mFmwTGOobZiygAmEEWAdZNkdXbVbrXoj4qotEFebTGNX4NYGWgowImDbdUlKtT3JodlCJCUerc1+/am1u28u5adOmioKJwO3Num7DUtF0bXphx+/euzhX2uazpluxpKT4cqj1fsvPrt8840+9lbV7I3CHm7YLa99Le6OFpkRbADw33zr2j++v+T+571Y14mp3IsPnS6+XG0Ep5Yq3vscy6gsVbxf6sW6h0knFjZLRfI9q7eZMxSjmGIaRiLKNrszu41TPCocRrAdhE5dtZvi8wYwLhPwJWBnMfFg0KTZqo1WC81WC41mC0Z97bZrOGdtGay1TOf0YdbvqHDtkerci88tXX/N1MIty0TovHixjrcbAzIl2P7g6r1fvzpT+X+v1vwPH2Y9RIhytjGbd8y58ycLrxiGCF99Z6UnQnDYdCLYAOCR+nw44buCh1QKwTBEbDFLBFmaAKCFWXeYhDDiwQm2g7Cbq1bsmIhxcFdtDtzwQF2V+zkKCEARc1T3fLvdEue2WvCCgzcZsFVUebQ6++KzSzdes1jqGKzRxkfnxYv1b51RMiXYXrg++wCAf/DmrbVfXa4czDVqGqJWLlhvE0Gu1f2fPz2d/1YrkCeqjeADPR7uUOhUsJXCFh6v3Q8VU19v9oYQbS7M2GJmmyZIjFR8eGYZBcHWLQd11WrBtjN7uYullJt6qqZCLupCf9kqqj5anf2uFm5Hhhb2T6RwoYsXD5ysCTYB4L9Uikvff33hv5KKC92uY6pkvzJddm7cnm/8zmTJfu3hM8Uf37xf/7V6M3xPH4Y8cDoVbADw1Nod5FXgK6ZDl0MQQmwul5GIM6GFWV8ZR8F2EFJXrSlVoCJp7+iq3Sr6hj3oASAASQDLLsMf/CCM3altIq7pbS4CvBVbRdVHqrMvPrd041Ut3DSIq8Xsaa1r+7+n4+0OT6YEGwC8cH32rwN47tV3Vn6r0gg+1M2ytiUWCznzbt0NL1qmWM3ZxuKxSeftSj149LBu1qzQjWA77lVxvrHMBPBOrXR2QhCtW8xsKy6VYZsWRB/7ZWp2Rwu2zVgEP+TO6rHt5KoVkdxu3RvhrNpeJmO0FwFuNFvrgs7f4lZNLG7f0TFumi6Q6Dze7uB+/DEnc0HpAN4B8NxE0brTrWBzLGO5Uk+WiRQpZuehfOlHCyutPdPaR4yO7ykVp4xz7qoEq21xbGmR2a0Ws26LzGo0A6bj4//gBZC3JlxkM6vWACKPe+ci3lQE+NhGEeB1t2qzBTfusTp53XL+2ruTZz/+WGX2u88u33hdZ5Vq9sEAMJlMe3L16YsBOo+3i/o14CySRcF2E4DK58xKNwsV8+Y7jVb4pCHIlYqLUnFeKXZuzzc+cnzSecP1osf7M9yB0/EdQpJANVfGGdmMDMO02hMAhtHIXKPJMkwEtkzA6vyy2FFWbZ9ctQ6xO4hSJ4ZhYKJUxESpuOl9Pwin3dbTv/PT5i/86unK7I8fvn/tNbU850LKbLltNKOGDeBYMu3J1acvtpCIuhuT53PfeOhDHz3nrnzzr936/g8ArAGojJOoy5xg+9QT5/wXrs/OTJedvBDkKcW5TpazDFG3TWNZKs7lbHHfENRSzLZpkFfMWyv9HvewEQByBiFvCOQFIW8QCgYhP3HWPOlV1ShWitfg6PYZ2wXF2TqGWQiwLdBp5RViBlJ3bCLqNhc97sxVawJBi6mj0Ih+4dgWHHsSmJws89mzvz538YMfPU3Ru8da1duqXmlE1VU3Wll0w+UFV1ZXPF2TV9MH8gDynmGdff3EY/9pKWidIeaPAPg3yed89emLNcTibfVfXvzNZ6tO6c6xVvVf/3++8E9GLg4zczFsAPDC9dlfAPCbt+caF9+dq/+tTpcjgioX7DcMgWB6wrmxWvWfMgxqVRrBs1LyWGSXWQbVTGAiLwTyRiLMBMExCLTLs3s5cFGIfBaAYoC0cBsd2jsdaAADkHJIpWoGiQCUQQgBEJiZIski7mjBpCSLUKrQC0siCCD8EBRGmYm/y4EbZw357jFSq+mYOIoi2ai6srbmRpWVRrSy6IaLs67ymmNj/dAMBwWiP7/w4d9ddcoXAWAiaN76a+9+/39NP5ckxIvn3v+risiYLx7/kAJZkTB+xER/87//0/9hbngj756sCrZjAP5zZqbvvbbwX0vF+1btF0ShYVBjsmS/KYgi2xJ1z5fHlqvexwYw5J5DAJzEUpY3BArJ/4um8FWHQdcpBiscb1UBcHpRZ12AdDTQgm0zYswtxQTAJPgR714kWIAlg4z2KzcpBeGH8ZSIOOEHEEEIUsO5xufBtXOGfHeaVGW3eZTX8mSj6kbVVVeuLbvh8kIjXJ5rareqplNePPf+X75dPv2J9veOe7U3j3u1O3dLJz/kmc6pnc4lJcSKb9iX/ofP/Pc/HthgD0kmBRsAvHB99v8C4MRf/mzpf9f0okf3m18QhTnHuGMZohFKVSrmzXvLFe8XmbN/s1sXZkLEbkyDkBMEsUORWYPgyy4FGwCUgiYKkQcAIIAJ4HG+8Y0LWrBtZt1iM9RR9AeLECiG2O/37iZTFgBEGMXibYugG5RVrgiunBfRrQnB9U7mZ2ZWbr0p65VGVFlxo9UlN1yaa8jamq/dqpp23jp24amfnnzi7xx0eVtFd20Z/vb/48t/+Fovx9UvsizYfgPAR9++U/1wEKryStX7xf2WIYIiUKiYD113rB/YgpAXcWxZrs1qtpMw242DCjbBCicSKxsQizbEr7SlLcOMu0WpW8ZRsBlARAQV8fZ+ojvMKxVg9GL/B22VKxOvnBfRuyVi9yDLcxRGslFzo+paQ1ZX3HBlwY0W57Rb9YgyVzh24tsPfPA/k0SH6g9OwMzfffvrj45CnbgsP7m/A+CjxZy58vj5/Evvzhkri6utjwaROrHbAswQjOGLNYtoI76sLQnAGGJbJkUCTdNZt7IxQCJuNKxFm2ZkYMSWqLADcZN10ji1iOF0+twsCJHk3sTwsRCQeQcyv/2S2Q+rXJ3p+DVpHZ8kXjwvopkCcaub5cm0THPq+KQ5dXwS2Ej6V17Lk/WKG9XW3Gh1KY6P027VsaZpOs73z73/7xxWrAFAIfIsAKcBzB9+ZP0ly4LtLgBvsmQvvvLOyt92W9ETwx7QVsx2YSaSDE2DYGa0X2bTyiEf+aDEPqEAcVSCuDVjxUjfiNuEmt2Na9MAougA1vWDoCwTyjKBLalavbDKVZlO1aR16nEjenWSVPWwYxW5fE7k8jnr5NnjeCx+j5VSqllvynrVjd2qi41wad7VbtXRR4HoLx744N/wDKuX9VWfhBZsB+dTT5yTL1yfvVEuWDnPl+cRX6SHooQMwqbA/7R0hiWyKcx2Q5FA03JQDL2N9wBDizbNKDFqt9u2qwSbhKBboZYiCFIOOSZ3T6tcEMZTm4jbzSrHAKospnsh2HaChBBGabJklCZL9tmHNrbb7latLLvh6qJ2q44Y3zv//k+sOaWnerzaJwF8t8fr7DmZFWwJ7wB43/sem/4jP1CFd+fqn/QDea5fGzMIyLUF/ucTq5k9YsJsL1pmDoVww8rGiEVb4h5lFX8NOotUo+kRBiGIwzUgDiLUgNi6dtBlB4WyLSjb6tgqVw+CaQAzgxzjJrfqwzu4VaurjWht2Y2W591wZUG7VTPGzclzD90pnfqlPqz6/NWnL5YuXrva6MO6e0bWBdsNANz0oikvkOVeiTVBQF5sDvzPGwLOGAmz3YitbDkUw43wEY4nQQCnQe462F2TVUbpLE0taoddDxEkRiDjfSd2s8o1gYmI5OOlMKw7vh86vh86nh/mPD+0glAO8nfe5FZNx526VWsVN6quxm7VxTlX1ivarTok5gvHHtp/ri7hde/d4wBe7fn6e0imLwCfeuJc84Xrs3fDSL33+IRz9+5Cd8lFhLj6f6GtyGws1AijddnvLU0zh0LorVvZUhggBkiX/NBkGc5Yt4N96NWFhsYl2SKFAayycU4WzHmrVAgixRZA7EkuFwWv5YPAt4MgcrxgXcw5nhcaajBqaZNb9dzDG+OOwkjWq42oVnFlZdkNVxYa4eKcy35r5Crnjxp1u3i65yvdqDj/JLRgOzTvKMale0vN53abIS0yW9jSmim3R/X/owwTwbVyKIU7J2kJQOmYNk1WkYB11CzAqZWu2xpsWYcBqgV82hTsBQqFpLg3F01R8fI5x8vnnK3tws0wjBw/CG0/GIpVjkzLNKdPTJnTJ6Y2u1WbnqxVGkm2qhutLGi3ag+JSIiGleu9YNvg8atPXzQuXruaWeE9EoLNbYVna274FADktgT+p1YzymhmZlZpmTkUIg9ih4dVfXXRZBkGYBAijrsBHCkUj1+MqQJELNY22CvDLLIsM7Is0y0VN3XAEUope13IBdGgrXIiV8iJXCFnnTq3XnqKlVJJEeC4Jdfakqvdqt2z5pRKL5579m94hn2yD6tPDzUbwMMAbvVhGz1hFATb0uMl+207Z56AZJ7MmYtuEB1nhggklyM1/Lprg4R79JsxEZpmHqWwuel9irdxZCwXmtEkYthHsQtEur8GEAmCZACKYaoxsogzQJLZEkRBN8spIcS6VW4Lw7LKkRDCKE+WjPJkyT738Lp1SIVBqBo1N6qtubKyot2q+/DS6ad/pWYX9u14dBAYm6w9T0ILtoPzqSfO8fVbb/7slbz1UcsxZkPJuZwp6mut6CFLUAuCWBAiIpJBpCbG6clzJxTDMAihZFiHXVcrKaQrWLW/PbTyKRpNNxBBHVVzsATM9hIfAiwNogiIBdyohzR4EZdNm1Z6dSHq0CoXOr4fDcIqJyzbEutu1Y0So9vcqsvzjXBloQU1pIawGSAiIVZyE08PaHNPXX364lez2vUg84INAB5UHL2x0Fien3A+zKD1C5EvuU20MExBLUGQgeTSTusZF3olSpkITSuHUrBhZSOwUm3fsUaTVUh779dRIEO1dUBIBVzEcEbxS/IUSlbEnqC4WLBiGI6gpmXA7+V2smaV28utGtXWGjLOVnXDpfkj41b1Ddu2ZVSLTKOw/9yHZhrAcQDLA9hW12ResHlfeUMAKJ+W6s7tpeZH7ZPFXeeNFOcN6u0J3QsIYEHYsTDjHqfbrtcEQRQp5p78di3DQYHarWzauKbRHBYBSIOGVzstFXACUCYhRJwBDjBIAtYo3OYbER9HOm4AnuSJCdCCbZC313K9IitWuXa3Ks5fWH8/dqtW3ahWcaO15Ua0uuiOo1v1funE2abpnOnjJrbe9J6EFmwHpgRAXGuE/4EMFUJTwJrO7zqzZDg5U1S9SE3uOlOfoNhoJQEwAQoUZwwrZrMXLswUZu5ZjFmcMZpHOYhLpujOB9mCAGVoS1I7TEkqYRZps2wdqJtBr1GAUFvGQQBbsYjjLFvgdvIkEEHtNO8g2csqZwVhZAeJmPPi2nL9ssrFbtWTU+b0yU1uVdlyW6pejWvHrS01ouUFd5Tdqucay3P54/5Cy3T6kiEaCSMvQWRsmCufAvCDfmzrsIyCYJsEAEdQHQCiNQ/CFDDKu18LB3lUEqCSGDpmZpLtdZK4P6NhQBgEX8VZcoe+DnimjULYgsEKAhjNs3pMYUAcpfIVHcOAAWTKkhCfO3TgbgaDggEKGTYBECOUuFEwac0S3SUiDJrQtszQ7sAqF4u5qB9WOSNfzBv5Yn6zW1Uq5TZit2plxU2yVRuyUQ2y7FZlAG+dePRZ37R72Td0E5aKGsbmL+HBq09fzF+8dnXnuldDZBRO1EkAiJjXhVCw3IRjCoj8zkarXrkL2zEF+e0ZqQYhAADJsCXDHvRBLxkOAcoUFBw2U5YRW9kmAjctnqsFgibTJHXYMmUFNghh1sXaVoggLUBKhpX1BwNfctEkBINyifaSXa1yDLbCULZb5XKJZa6XVjkSxt5u1epaI6qsxLXjFmddDvxMPAzdK506c33y/G/2cxu5KKhseUsg7nrwRj+3exBGRrA9UrDeWgn8JxkwwIC/4CJ3rgyyt1+zmfe/kOdMUWUAzGxIxbZpiBYAKMVWqLiQxpzZhnAjxXakOO8YVJMMJ3FxDr3iOAMiUuwYhOCw4/FMG8XQg81yW9FcA5AKMLL7HKY5SmTxeDTjLgQjJdYYQJSM2QAiAkSWvtOtRAy7GvJpW3KraNKaKSgc9pgODYG6ssp5Xuj4ftQrq9wmt2obsuW2ZL3iyuqaG60tNcKleTdaXRyoW/XWxNkH3jzx6K/1ezuBYZXulE6deaixON/29pPIoGAjzrA5FAC8r7zxSQAf/upS85euu+GvtH9GpoBzvgwytj8YWoKahqBQMRuKYUaKc+l7oVSF3WLKDEJgCApGKdO0F4INAHKRj8nA3VSUM41nE4Dk+IKe0eih8eSoVfTfj6x+H6PegYCS2MAsfrc7kXS3qRctqogMxLUNlNQq5wehHbRZ5Tw/tML+ZbCykkq6dVfWKnHtuNXFRrQ07/bDrfrSqac+8Pb0Q789qPuNwco701x9JRcFjeeWrr+cl0EVwP/34rWrmTq2RsbC9nDefPdmM/yYahNaHCkE8w0458rAlk4HoeJC2PYwIAhy63s7IRm2lDx061k39MqF6Zk2LBWhGHmSQYIATq1takO0ZcoNpTlaKEBkMSkmjEMUmDLoqu2EpHvCSIg1ILYOeorLvs+lgkmVvEF12toceVxpt8phi1VOKrXuXu2xVY6EIczyVNksT5W3ulVlvdqQtTU3WluOy44cwq1aswuFd6Yf/A8HaRyQJHL3iyc+arJ0H6/cu5qXQQDgQQC3BzWGThgZwfaOG75X7WAVU75EsOjCPlXcJto2zdeBm3QU6ZV1LYZQtwvEIKMYeduetrVpTZMFkocHxRlz4TFABkGO4rVmr/g7kxASNqxYDFCUgZAQIB6LG/G0F3G5YNJazqTm/kuNL8rYJ1YuEXK23zurnLBsSxw7OW0dOzmNC0/Gm2Nm1Wp6slFxZWW1EVWW3U7cqr5hWUu5qeMnWtU3F/NTzx1iWAeiEPoLx71aNXn5JLRg65pJZkZTqondZpBuiHC1Bev4IOrqZYveW7wIDTsPJqLilubwo+Iu0Yw3HE+CABYZc+ONcjFfK65hSWGbGBOAUgxDYfPDshknXVHUw3JFh0ECZj3ik55kv2yJZUPsXPfyyNJulSvvZZXzQ8cLQsfzwpznR+KAMVNEREahmDcKxTxOnd+UrRq7VddcWVl1w9XFRrg45wbNRnS3dOrclF+v/vDse/9Tk9VQMjTrduGRtVy5fNyr1RALtq8PYxy7kWnB5n3lDQtAgYhg7tNXLqr6IFPAnMwNaHTDp1ctqrYTZ40ygFKbaGOABFjqTgiaLBAfj9mJXxKAGtU4tpDhpFbL9D0ByCQsYtt9IrWwmYQADIqQDeEWMhxfcaEgqDbssYwKnVrl0m4Ph7HKpW5VozxVrp19bGKJxdlpUks/vXn72DsrlTNPrd39MgBEJHYvttpH8pG/cCwWawBw8urTF6cvXru6Noyx7ESmBRuAdataJ8IkXGmBTAGjmAlrfd8hItXPciJNKw+ANjWIT2PbdPKBRjNeJJ0ZYiEGhBIw1T7neZKRnymrYqTY1penHtBHqxwBuK2MiwEoV2FxJl8sLvJK1bg2/dBv921/OqBpOmd+duzCU5aSwUTg1s40V58E8ONhjqmdrAu2SQAIFRu1SJ3vZIFg0YVzTkA4Wd+1w0GAlIr7/lTbtHJgovVOCG1WDcrUVVqjGT4jfUqEDNsi+Iz4hsr79JMw4/hZK2sPb1JlI75unNnHKhfZfhDtZZWrKpoIQOvusGK+kBmvzasnn/hbACBYhR+fff2di1qwdcwkANQilfcVd9ZqigF/voHc+QmQmZnQlp4jCHJQteBapgMGMLHRvqqvmXqJIExvfrvdDHa4ObYvkt52eko3K+zm5r11vX0ZvKa/jMMP1qlL10QcjpE1sQYAErBU7AoYaQHdLxRDBJLzIm6jCAIxiLk9/pJBgpmp7felVMAnxdUJ6YWWQclcLADlCRMomIxCIV0jAQApxaYfKNMPVVhplFBx18eUz+czc7NOwwIkCefW5Lnfufr0xT+4eO1qJnqUj4pg273j+05Ihp+U+yCRuevJSOKZDpgIk34DQH9NCe3lRLpecsf/jy6jHMR+1MiieOkHcb02FrxXWv4QYQCRYmcUOyL0E6lgtCRPtiSXN3+y0yVm58uOCYQHj1ckwHIAy0HRi3wDG4LNNE3DMk2EUXZyRWwVVRnkfPPBD33kIvDtYY8HyFB21S5MAsCNZvh4twtyIBEsNJDlPmkHRcTJBgM3+/uGjapTQnpfSupOaTQaAGN5sdkBkxBkPfEoK9mrWYAZVA/42FqgHtgu1rqDelCk2AAkB+Gm34cJVMwPJc9gVywVuQ0rd4KBjw17LClZF2wTAFAN1cmDLKxaEYLl8SrLk/QwHZpO8g0bFaeE9ILN6E0Tbm1J0ow6dAQMbARwVmqw7YVUo1X8vF9IBbMSqLOe4nIvLrC9sCILQoQg2qQ9FEMU8tmq8OCauXORMAsMenXYY0kZCZeoSThwzzhZDxCaAtZ0ttT7QSCAVQbiRgLDwppTwpTfIAJDAcZhMkfT3pBZbTuk0XSCBIwsdmHoJaPSMzULvZ6HTSDZqYd8qpfX1F49WFO4+ZbOIFHIDV+wEcBTfuMdAAgMs/zR2Tf+25Ne9cqwx5WSWcHmfeUNQiLYDmvejtY8CMuAURrtc9gQFESKM3GxDA0LlVwZU14dACOp4bTjjSpNAdjts7SRdxrFqtGMKqloU2Pad3dU2lcddZdoIDlfC/lkFo9BxWyIYHOsGmNwLlGTVau9zltOBstSGHnBHD61evvr+ch3y2GrcbxVnbVY/ruL165mxvuTWcEGII8kuNFTvGuXg04Jllw4poDIZXmXd4cApZgzNfhQmFjLlTHl12G0he+kjeKxcbFgAXBihKO0NEhbz0UGQAoQInavZuYE0Wi6JbWwJQf/2NQsFIAcFSHEAEUKpnkEOx74kvP1kE/14yIaMezD1uHMNZoBbbmXMXrvEjVY+QTIiEQBAEyWTYCoFDTvNuzCg4rIFMzR+1be/fOGlZ80WEUPNhbvTvmNerJzX7x47Wqjp4M6JJkSAFtYL+MhuQf1xtbLfZRB1uh5LJJEg0xY19qJhImKM4FpvyYMZknYfMNKmmEn1oaNc3yL22j9A5VYKAY1fo2mX6RW43GxuAlCNEp9UiWzbYKOlGBTDOGGfKxfT7wMkAFEO3W/2A8C2CQEvFzb1kNSAYZlmuhVpuhJr/r6Y5V7L//s2IVfrtnFRx0VVU43V18PhZE731h+p2YXbtwvnXzujLvyxoXa3B1bRpFr5XKlsNVKTtK3Ll67+uahB9JjsmzengAAZoYnebona1QMf64BlpnpZtMNmb3YR8LAqjNBEQmj/abE0P1HNUcbBih5OKFRfxBRnOkH/G1EajSsgb1AMagZcbniq7MHEVOdslursr0wgMgi+AYQST90jMb2RMA0nKZXblGhVBgYds5gDgmQp5prbx7zanPLuan3/PTkE7/79Nqdt967MvOND89f/bEjw1CAubwh1lwAmYlbayfLJ+AkALiSnZDRXR22PeBIIUhqtGW0jNCOKIZhCvKZWfSnf+jhkMKI3aNeXRg8koJYo+kb6cPLXvGcWSYZ98hY14DYwpbh59yeIRWMtUCdH5QFt1uXKAMiTMS+U6nvuCAjLpheyOdEpV4/1PgEOHykNvv6Ma++drK19ufXph9+zpGBOxG41XzkLb5ndeY7E0GzORE0391lFV+8eO1qJstLZNn6MQkAXh/aLylfIlh0958xQzBgRIodybBMQZmourwVSQbWchOQ2S7RpNEMhSR2cyStbBTHpY4U8ghY2BohT1cCdW5QYs0kBN1uK52fmGGt7S7GCFC9yBQ91Vx79ZHq3N254vGzppKSQXRx9fabq87EiV+599PPXKjN39tj8dcuXrt67dCD6BOZF2xddznoEOmGCFcyKaL3hTm79WoVCazlyoiEFm0azVZkUr5m1OoOZvaCswcSsMa5lHEj5OmW5IlBhZ0Q4qSDbpdLWg3CaDQhwt3j0whQB3WJCnCY/i2FrVUpDGM5P/kAg/Bwff5twcxPrd15m0FkstrNBVQH8OcHGsCAGAXBdqjKzHsRVX1EtUwaq0YaRQJrThmhyLLHXaMZDskNlgxApjezrDOKCRMMIOLxLKDrSS605OGrJ3RLtw8aFiFIE8ys1dp+61adZIoarPycDFdMVi1bRTUBjk60qm85MlwtB83bj1Rn37ZUFL1ndeaVicB1JQlRiDwvL4NgImzuZaX5wsVrV1vd7N+gyfIddQIA6n0UbAAQLjdBpoBRGB3rucpgDNtWmAQqTlzyw1JHKlFLo9kXRmxtSzKpM136I4lfy/LD/a5IBcsSCIY9jl6iGOSG6tigD5m0q02XFj0GABFGMHdINmhHgFUnmaImy6ZgFZ5o1d+qOOVHbBmuPbd0/bsrucnjx7zqyslWtcIgiv8Cj9bmtrlA0wzuNl6+eO3qjS72ayhkUrB5X3lDYEOw9SZDdA+ChTgJQTiZ/Do2YRD8LJb32AkmQsUpY9JvwFYHblah0YwtSbB11gXbyHYgSVx4oxWwvA/NiKeG1cuVCKpTG5sAZNoVw6rUQfstR8TgOFN0p8QDS8n6ca96jQAOhFnIRUHtbLT86kP1hesnW5W1icBtWDKKCHHMUN3K54uh59EOPX63nGwVAF/rbK+GS1YVShnJd1qN1Km+b40RZ46enwCZ2b4uEdFI9ZhmIlSdEiaCBhypRZtGs5UkezSzom3U4u3aicYsUzRUbHtDcIUeBKOtbt9+7lAAII5PgUI+t02wOTJcfW75xhfBzAYrdXPy/LOP1mavnmxV1iwlJQAIZjbabo5tNdX24wsXr10didiorAq2SSCuwVYL1QOD2CBLjgvrnisnyfeaXsFEqDklTPguHDlW3gmN5tAkLkd90ekDUsFWDBI0uqIzhQG4IR8f5o50WjjZaus5a+6TbJBCBAYDWzNFCeBHanPfP9msLJcD1xVgPu8uz9syDNtPGltt9qN2eEL9+OK1q7uV98gcmRZsy4Eq97IG235wIOEvNOCcKY1UjbZRgEGoOkVM+IScHImHGY1mIKQ1qLLqdhxlpaMA0Yx4UgBKEOTGRFG/RRwz4EsuepInetXSa5jCngAIihNlwn2yRdu/2E6sa8n6FbC9eC4DtJKbeOjngrdfSd9zZNgLd80qgG/2YD0DI6uCbQIArrnB44PesGpFCJabsE8OTCceIQg1pwAOgHykRZvmcKR9Z0e9g8AIMLJPrwYQJYYbIRltRccZRKQENok4KejwWbuKITzJ5VakysOKNesHDCBkOBbtncQhAJWW/xCRhFnvLoSwPVOUAP747Gt/eL6+NH+AIe8FA/jcxWtXR8rlk1XBNsnMWA3VGUugESqUBrlxWQ8QWQbMqd42o9UAAKFuF8AgFCJv2IPRjDBJ39qsXsO6QovO/mAKCnazETKzkNgm4lgA0WFEXDVQZ2KL2sjq3P3Y0zKZWMoEAJidJBtsrJUAoD1T1GDlTfjNutH7wO0fXLx29W6P19l3snqxm/xRxX9m1oveP6x08nC1FZf7KI1lGZ8hQ2jYeTARimGmy95oMsw4xX0lXRAy6xYdRQhgg9CV64yZSQLWDiJOEiEyNoTcriIub1ClHvHJw41+dElj0YDO3aHxghvLFfN5YHnunafW7vxgKmg0ejzEJQB/0eN1DoTMCraZVvheyRiqiStYcuGYAiKX1a9plCG4Vh4MoKRFm+YAEMCj2ptzKwoQGS6iO5JfsUnY1brWDYmIM8Ew26LaWRCkQZBEkEYcE6cAIGdS05Mc7BfnNa6ksWim24IRdK6Xm0yTBeJqk2mykM/h/fNvffmEV6v2eHgKsSt0JIuDZlWJTLoyA08ojDhz9HwZZI1NKEKmaFp5AIRSOJptwjTDQwKGAUg5Yk3JR41RNWPu5Q49PEyKYSpO76EMgLgtsSHCERVs6TfelXUtIc1Cnc7ZjUnf7bVlDQBevHjt6mwf1jsQMifYvK+8Yc/70RlP8rFhjwUAoOJyH865MsjQ3op+0LRyYCKUg7Gqb6nRdIx2h/aWJNlgwBbLDREnOXv31kFBACiSMGvdXc8JgEPctIHWybxx22Ipd+hIcBjmAXy3d6sbPFk8qCZmmtGDwx5EOxyquLDuubIu99EnWqYDBjChRdsmkkzIgR506UUyfVLu9cYVYPTK7qHGxLqW5dZPo+gPNYbcjsokBAdplD4OMENY1QaoyzwBBlBlcWqa1NwDBWdtBT299kjErtCRTu7JomCbvNkMnxv2ILaifIlg0YV9eqAJq0cKz3QAEKaDxijeI/pCIm7G6ilBAJIAOqxFKcvdAbol47XYRuo7JgAGUTRMqSkI0hHk+oqPVH0oi+CHDKd4AHcoEF8bHhDyXeV2aZ7bn29fvHZ1ocfrHDiZuzh8eq7x5FqoHhv2OHZCuiHCVR0g308800Y1V1Yjdo/QdIECjDTIXoAlEIsvAUgjmUQ8qXTaaT0ZDtI/KJl8UMmy9W8nCJA79Y8cJAxAMVtjeIzui9H0YPidGzhNICiAayYQ5IjrNnEYzt9d7eGQ7gP4fg/XNzQydyL+7tnSqyWT5oY9jt2IKh6imi762k8C0+aKU4IWbeONAoSKI0MlAawAQyZTKurSKRVyaV/LdP5h70Mv4aSn6LDHsZXERZ65ce3O8IbKACKG40lMSIY5LhbgdoigTIJvEQKLEJjxFJpACICKazXPBAIDCAWg9voCLMA3CGETNBEBdsBUDBmmf+dmrwRbhNgVOhbCOXOCDcDiLx/Lf0YQMpt2Gy43IZu6kXk/CQwLFaeU1lLUjDGpONtrnlTIMUBpkVkxZsVm05t7amHM0pE/SoItdocOFgYhTIRaoJBL3dvjKNgUw4gYTsiwQ4YdxZMVAVYUSpsrjVwE2BKwku+g3YIemuDAAjwbaCnA8NvaT0rAvOPLs98WUx/o0QH3zYvXri73ZlXDJ3OCLffJZ+QjBeu1kpFdKxsABAsNKD+zmnIsCAwLlVxZizbNJiRgJDdEY5SERCcwQKmFMWOCdCS+ZwGouJzHYGAQhQqOJ3kiVMglBZDHttQMASz3SKawV2s7JhtwHLNqSMCKQHYI5AIgv7VTCQE8V6kVjnu1Xtz/7wD4cQ/WkxmymHSAUPFsU/KJYY9jTxhx5uj5CZCZOd07NoTCxFqujCm/DjHcsBRNhmgTNKwAcxyPjNQVPOxxJGT+KzYJoS3QHMRQGUSRYididtKtpfFq4+aqb8eIs1+dnT4jZlirh6tze07I64U7r3/LW7uzeKgVxe7Zz4+LKzQlk4Lt9XpQixj5YY9jP1jGNdpy58px9ImmL0TCRMWZSETbWJ1/mgPCiC1ttPFy7E5ARtyGkTPgWst68LwAlCXQ9+rbDBKJULPbZeFRKeC810FoVhoQ0cGeLwSgzgh56zRF8yu331k72Og28bWL1672MnEhE2RSsC0H8jYBahSykziQ8BcacM6UdI22PhIJA2tOGdNatGmOEKnbd5gtuJJtZ/JabBJCAqQpKKAef0OKSSiwoRgGM8y4xM7mbSS/jToKYg0ANnqsbsdeqRx4vTa4OUGqEtQqrmo2DhsgfgvAS4dcRybJ5En4V08W7gtCpmPY2lGtCMGKLvfRb6QwsJYrQ1ImD1vNEEjrlw17HP2E40fBoe1j8vCcOVIXqCXgH7aMh2ISEcMKFfK+RKklMekpLgcKhYjhyC3FnglgI85ApHF2gbazV51A023B8A4eOuiDSnPKeOgHa957bkyee+jAKwJ8AF+4eO1qFg/ZQ5PJO1/uk8+wZHxz2OPoBlnzEVW8YQ9j7JFkYC03AUlH4hqp0QCILVzGkDJHs5rYQXTY+D5CxLBbEhOpOAsZ9lZxtjF3LFoMQCYuasrkF9Mn9joO7OXKodaddjm4WfcuvHzqqf8oFMZBL/B/fvHa1V43jM8MmRRsANCS6l8Dw20v0i3haguyMVJDHkkUCazlyoiEFm2a9azRrATn94U0Zm/c97NTCIcTkopheJJLgUK+k/jAxPUpGRBHxf25E6l7vh3hhzDqhw8fjCKJmusiFGb5lVNP/Vz7Zzcnzj14t3Ty9D4/+DsAXj30QDIMcYYz7/7HP/j+b0SK/9ggOnCs3faYg86eUdPlCAQGb1puy+v1gGdOKilZZ8s25fZOHe10HDsO7ZAByFv3Z6fPU9L939pY8hDj7wSOFOf2m0mwwpTfgKnGt7zKOLVf6idtGXqZfQjtBcMIbs9iyywBlgZRZAl05daIy3BwrtM+n8lxtWOR5qQobHZvoG304pjZeuylF6X87JJ0Vms7xpy3fT/sgFvAeuVrIgY3QRPpvEura7h6610QwL8w99YLj9Vm76afNU3HqdnF0szEmSd+Yf5nP9pheC0A//zitav1w+5nlslk0kHK792qvHrWMUqXThWcnOjV9aLb84u3/N38vkHw5ZY052C2gdz5Msja6xwZ9nne2fYJvFkw8Lb/9BxTUEetJBQJrDlxyQ9rjEWbZn8yVgJjrMhawoEAJIOEZFgWOhNsDCBScCLmXDdXLgJ4N7GTNRG7F70Q+rzFVsAASEoYlbrRybolEPqgQrpwgbjafhtZq9XSj+gnZy7+ndVc+Vs/v/j2ywBQiHzfUFK+fuKxU5KEMFhtjen88riLNSD7B9wH5nzp/Om8i4YcobhiFZf74FEa84jCRKg4ZQRi1+QlzRGBj0Ce9qAf80TGEg6S8RhJIVYRKBR4n589YtiexETI6EqsjRNJ3cJDWgS3f8/Wag3U4RGy9XdKM07TN1erGw3jIxKFqlM63T6/o6Jo2qvdX8xPTW9Z9c8AvNnRIEacTFvYAEwCuL0aqoc/M9fA75wuYnpPq9VQ2PFo5VAhWHDhnNXlPvoNE6HqlDARNOBI3TLsqJJYgsayJhsQW5YGkZG49cbeL8slY3uRWdr8cfreJjdH+wU3YlhSsmUJtEyioP1TyTBDhdxRyeLci7SeXy9d6sQMe7W2/4xti7S/MAmBA262mErVlmcHYSgAwFZRtRB6C3WrcPbb55/79QcaS9c9087fLZ16b8PKnz/vLt9uW40L4Mq4ZoVuJeuCrQ7g2wD+XkMyPjvv4rdPFXDayfqwY5QXIVh0YZ8uDXsoYw8ToeaUMOG7cKRO/DiijG28X9qKaxDb2ssN2EsMIKJNMcCbxFh7kPCevykDCBTyEdgRSeYoM2hr2yPN4dj6O5hVFyLsPBSlfXkC+BjJeQYJAYrerVbPAIDJqvVE5d63IhLW29MP/VbTyp2bLx7/gCODSjH0Fk43V1/PR357yMyXLl676h5230aFrB/QJQC3ATQBFDzF+PyCi0+eLOCh/Gi4wKQbIlxtwTqW+cYNIw+DUHWKmPCBnBZtR5JhFpjtJ6MS3N4l1EuBrQChOPNhPqPMpt/qMIVyASBPaK0oOumDirVq1S9G3vIn7v30T6b8Rv3Fc8/+SjHyZh+r3P/+xdXbbxOYTVaqZhUKE2EzTUl9/eK1q1cPNYgRYxQEGwO4B+BJAAgZ+NJiE79+Io8nix0l+gydqOKBLAGzvGMLNk1PIdScIjgg5KOOchc0Y0LSeFshbuU08iR1v2TaEH7Y49GMPmkP3oNYaze14mp6MFr7X1/bs9zbiz8zQB5TTgKmGYX187ff/MFJd3luym/U3504+8Cj1ftvfXj+rRcdFW0y4bWJtTqAr3S7D6POKAg2AJt7xCkAX1tuoaUYz46ICAqXmiBTwBgRy+BoQ6jbBTAIhUgXMz5KqB7H6QyatChrmpk5qvuhySaJ23k9AaHT42treaFOC+VuFmwbqyiAa8eFXD4JLDavvnRjeun6Pd+wrOX85JSlZPTi+Wc/ddZdefm5pRs/mgoajXR9EkRG3NXizy5eu3rk2guNimA7vtOHL656aErGR6b2LdmVCYIFF865MoStr8H9h9Cw82ACiqEWbUeNUXSNiiPUk1IzXBgQBEgjsbbtd660Cy8RhDBre4eNJa27QgMcCYI0gMgBN0+SunuM5LJJkAJg2Ww03Vd/cB8AHBmGlpTRO1OnHv/A0vXPGKwkgVmCqGU6DgG4Uz71wMW1O//+4rWr13vyRYwYmRVsF56/YgNwED9pntltvperPjyp8Inj+X4Xc92NzjeqGMF8A865MmjvurqankBwrQIAQjE8cg9jRxZOSkOPGmMap6bJIKLNCk0AC7BS2L3fnwJEWkDZXqmmnSZQJF5FHIogBVhybAGLCFAGWAKABBl5YveMkPMAEDEMkRzr7ms/ugUp14/7QuR5H51783tbz15FQrx+4rEPLRSmZd0ufv1ib7+OkSGzgg1AMfn7GIA9/YhvNUK0FOM3ThRgZryEBkcK/nwDuXPl+BlH03dcKw8GoRQevn2KJvswSIxihwgFGKPsztWMBrSlc0Psgqe07RltLQicxlIqwCCl4KzVVY7YLYCrqSjbINZeEchSoLgjBViWSK0Xta2zmJgmtRatLVW8628uty9tb4lZS5kIm82Pzb3xzwB88eK1q0e2dlOWzTypO3RrkbwdudWM8MUFF0GmyjzuDAcS/qIb31I0A6Fp5VC3C8MehmYAjJpQS0n7hQ57HJrxJbXi7nTnUYCRxoC2W3tTaxwBarJaq5zk6G4ZanW7WIvPvSbERADKEVhFIIsA1WJav/hOC7XGzFz/yXdvdHgPVAC+CuBzR1msAdm2sKWC7W0An+xkgfu+xJ/ON/DXThdRMLKsRQHVDBGutGCd0CJiK8wsOm1P1QEKyYNJaOfhCkETXkNgdG7qPSsEq+LsyVHZ70NBYMV7uHiyyijG3mlGgzaxtuc1IBFnEGCZukkdgYZBFE0tr9UBhgKl1xJKkmQoAtkBKJcHN1qMUgvGhEPs2sTecaE2WdKC2ZmFcO5OY6ftb6EJ4NMXr11992B7PV6MgmB7qJuFlkOFz867+GunC5g0s329jmo+yBQwRyRpYlBIhtUr66MpyI8Ur6cSR8KGtOOuCEeNpK3PERFso7mbQrtENX0iLXvTyTUgfkokCEDZBjUls0nVBrUCOQUIKIax3hc0QQCqQFxtMk0SWJVJLdtg/xTJOdH2HMJSysZffqcTATYP4N9dvHa10uWuji2jINgMJIVzO12wGil8ds7Fb58u4mT/MzIPpSzC1RbIEjBGpKbcOOCZNphKmPSPnmjTZBedJarJCkk9QzIF/EByQQKGuVhBi3fWeiJ2oyoDCCdINQRYAYBD8GzCprg07+bP7sra2n4elDcQl+440i7QrWRdsJUQu0O7VjNNxfjcfAOXThVxPpfl3QSCRRfOWQGRzXGOZW9G37BRccqY8hvQTqjxY2vg9Cigs0Q1/aTTgrlx1iikBMxQIc8AjJYPo+mtf06AssC+SfAFoBxwMxVp7RTbkg0AQPme33j5xbt7bJ4BfA3Aj45Kf9BuyKRCSJgA8HdxALGWEjDwZwsufuNkAY8VMlywlgF/Ic4cpew1tx9bAsNCxYktbaTvlWNFUlV93/pSWSApkSCPirt6VEh+jFE4hDpl3+OLADYIAZIWX8l5JJyVinSAwCQObLBnAgElcSu7PWgIgiwQb0rNb7718rsc+NuSFRJaiOPVbnW1V0eILAu2X8Qe9dc6RQL486UmPnE8j/eWMux2lAw/rdGW8YSJcSIwLFRyZUx5dS3axoikREbqihmIEOL1f9t7XG/6PG3Rw2mgdhJTJAbV2F3TOYlYGRsRnVjG9ix3IwDJDMMUCAQQESBFJHG8VmuYSQHcTjVsAdxoj12T9Wqj+dZLC7vMPg/gTy5eu7rW1U4dMbIs2D7WqxUxgL9YaaEpFX5+srcB/r10Y3Co4m4IZ0tAdurJjb2KCYWJtVwZ034dpEutjA0cZ8UOuK/o/udtW89TGvz4NJ0ybi5quU+dPwJgCgQGUdheerq8Umk5zF1n7ZdIbQoSbrzyvRtQO9bdehPAF3S82v5k2ZTT85PlxxUf31ltgTN8HiovQrCUnQKvmZGNfSYSJtacCSjK8imh6QYVl6beFleTBZKxyayOTzOeMNYTCnb8LFAotCRPSI5DkUgxF1YqQbfbSZIN1gVYuDi77M+8U9lhk18D8Fkt1jojyxa2FwG8v9crfaMewJOMXz+Rh5EdK9YmZCNAaAlY0/lhD+VIFYaKhIE1J7a0ie3xs5oRJHEDZfIQ5tjSkcWhacYU1db4fSfXaOLQJ6VgwkCQr9QCo611VKe0JxvERXK/szUurQXgMxevXb3Z9U4cYbJsTvg+gL6o7uvNEF9abCLMsPsrWvMg672qHavpFJmINqktbWOBjG9QuwU5DxU6QnXxRpFx/W1Sy/NuO5dkEjAAFJf2Lb+xja3JBv7t6/ej5fl2t9ECgD/QYq17snxX+gyAfwqg2o+V3/UifG7eRUtl15ISLDWhWtpSPGikMFDJTUCOXqF8zQ6kVoVhj2MHsvvEeMRJSmBk+f54KHifBxmDKHLqzdDyg67Pm/ZkA47CqPGT79xu+/hNAH+kkwsORmYPyJnLl0IA30Ys2pb3nvtgLAYSn51zUY8Ofi3v91OYv+BCBZk0EIw1kgTWcmVEQou2USdx82RRHI2lBWccGOcfJi1yu0fyAQtCVFzu3roGbE42aL392oxqNkLE59/XEcerdR0Tp4nJcgwbALwK4EOIrW2fAuDsOfcBqEQKn0n6jx7PYg00xQjmG3DO63Ifg0aRWC+ua6po/wUyDgFKdNzzK4386uWta6f19W4b7WtpX2u6wwYQyQxd8zh2Te1ZZkEzHLKo7g9L3B90/9ZnJiG0PD9y6m7XF732ZAPZclvuqz+ahY5X6xmZuXjtwj1sWNe+AuCv92MjrmT86byL3zpVwFmn66+k7xdbjlQs2s6W40u8ZmAoElhzypjy67BGXLTF9b66Cc7r9bG20/oGezyLuECtyEIbKAZgZGAcms0ktfHG7jfptE8tJda1g5yZ7ckGzdd/fIujcA5xfbXVA6xOs4VMm2xmLl9ixFY2AHgdwGv92pavGF9YcDGT0Zgx5Uv4iy561RRd0zlMhIpTRiAy3C1D0xEq7n5AWUlG0Gdz9qAxuMoagBSANJKJAO5IrAFsShkV1mpduy3bkw2iykq19fbr30Ycr6bFWo/ItGBLeA3xdY0BXO/nhiIGvrzYxNVGV8fqwEwEqhkiXGkNanMxo37l6hFMhKpTQmBo0TYOZCUJQQFGRuPrjixbfwwBlvtP2OW99mm397cuu9Pne61n87KJ699Iu2fI5CGlk30ngAur1eAgBcTbkg3Yn3nnBTB/Rser9ZbMC7aZy5fqAG4kL3+GjXYzfUEB+OZKC6/UOou3HPSVNqr5iKregLeqATZEm29kuMWZZl/SO1eW2kFp0ZYNksjN9ftiLHxo3wkAtr8Pg+NJqGQ9HIckJO1j42xNtW0b2GfaOt/mZdN1E8DdHuMmwS8vVw5Uub1EqgEhImNy+tMP/ckLn9fN23tP5gVbwqvJXwZQG8QGf7Dm4ftrA7ZmdUi40oJ09YPLMGAQqk4RnhZtI0vWSja09xgd9liOMonVtadJIIkVTKW/LQFKxQIu7R/LvbD2EsBG3PuTN7/d+fIWwStV664RRV0fhzbBcxynknv8vT/OPfL0v+l2eU1nZOaitQ9vI840KQKYHNRGf1oL8PXlJlQGr6PBogvlj3YQ/OhCqDlFtMyeJy1rBkBSODRTJ3WbaNMMgSRuTfRCrKXZmGnfztQlSVh3ga/DST/Zw8RTEmLLWNKwXba9v2Pjzp0wgEgQVGl59UDum3yxdCf/5PtfMQqlb+Y++UxfaqdqRkSwzVy+FAF4A8BJDNiN8bYb4st7d0UYzjWWAX++AQ77HDet7yC7QKjbBS3aRpDk5tkTy0YvScalsjaucSUJxkcicqgXCj5dn0qEWvtnvDHPpt+YNxbteltJUkEoCJFiJG7YeI0KEDtY3XZaDwuBKOe2PLvld3dDIeL65Kml009c/BGZZhVxhyJNnxgJwZbwKgB3GBueaUX4woILv+PnlQEhORZtUl/fh0Ms2ppWbtgD0XRJFq1sQDyuw1pcDkNiBRql+8KBSK1fSOK8enAgUNK3ds8yLW2WPGFs+Y2TzM6OLuZJEJw0BTxTkBcoFCOGE1sIN5pkE4FzBjV2O54EIC0BzyAKy8vdxQApYail4w8uOqfOXjMEKQDfyn3yGd1PsY+M0ok5h7itxVDMrfO+xGfnG3AzJo44VAh0uY8hQmhYBbhWftgD0XRJVjM0Gdkd26iTWMBUGozfq3i15Lei/cRuUuMNjDihIRVtadmNTsQ6AWwSAlMgEIQoVJznXSyEEcMCgJxBDRMIE4vbOgYhMgmB6QeRU2t0HGMTWk4wf/qR2WZhovVQ3pxD3B/01U6X1xyMkRFsbTXZfjKsMayGCp+Za2Ct327ILlGtCMHSgRJ7OkHfNDrAtfJoWIVhD0PTBemBnbTqydxxPgzXaFJYeFCBEOtu6UFtkJKEk17/2LJDKx1vOczSeoCp0EvE3H73ZZLMhkEUJV3aN5VT3/n3Y9gGNR0Dbs6gWuq6NQU8AOimUG4rX3bnTj06F1q5qGRQZdoymgC+mvvkM9myZowhIyPYEl5PpqFRl4zPzrtYSAL+s3Khl40AYUazWo8KTSuHhq1F2yjBsVUkM+U9hs2gXMW04YpMLUoswH1/Eh5mSG4szGjbsZYWc05fJxmf+zxEEBSzIABWLLp2zQ5VnN7n41mYWQBgS5BHAIRUqrBW3b/sABFXJ0+uLZ54aImNWJ+ez5lzAK7nPvnMrX2X1xyakRJsM5cvuRiSS7QdTzE+v+DiTiuESZQZlRSteZB1HUIwTJpmDnW7OOxhaMaIQSZIMECD2Faa+JG6f2PhTMbWuK5eYGyU1hgaRmy57OjBQAFCJi5xAag0OWILBACh4lziDl2fZavKCxQX2jVcxLBNQmAQRwCQX60GYp/47DRerTJ5upqGyBmE8MG8uQDga53sl+bwjJRgS/gx4hIfQyVk4EuLTbxZ9ycdQUMXkSnBchMqo+21jgot00FNizbNIVGJ9SmxRgkxRhmkiXzYphJ6nfCQJlGk3+WgramJ6JKdukzbSRNQGOuZrADiL00QIkGIiOIab1uPi60uUsm83iTbEuQlVjmAwaXltT2f8kPLCeZPxfFq7e+ftI0Fk+gnuU8+s9TlrmkOyCgKtq8B+J8BDL3cvwLw1eUW/rLqTdoCNWTBPcqAv+CCg549qOrCHgfAMx1UndKwh6EZYVJ3bSpiVI/qhHWw3UHAO+1IYuGTSRbkoYfSHpN3ENF0mO2mlr3DiMRErG1bR/z9UOxa5ti1nL7ftug6gUJRMpkMArVlqOWr9dAIdy+Uux6vZue2JSQ8mDNnAHz7QDumORAjJ9iS5IMXAPwbDKnMx1ZeXPXw3VVvwiI0BfW3dVZHKF3uIwv4ho2KU4bWvJpeMaCY2Z4esG3iZb3qPwO0W/mLtMxGu3hLBrVt39N193K8PYLSormHWkn8z7rCSr8LQRQmianJv+vfazrfpu+KAfiKi6rN0gbEyQY7b5i4OrE5Xq2dkkHV47bxldwnn8nEPfioMHKCLeEtxIV0Pz/kcazzUtXH15abRQCRQRh6IBlHCsF848iX+9ialTVoAsNCxSmBsyHajvbBMB709TdMyk709GAlpD05Y2thN+uPy5uwTGLAKI3pSv+mBWqzkvwF9K71mQBLk+BFjG198BSzAQCCODQJYSrY9lqfAUhBtC5u7aYX2U1vm9iN49UeWKxMbcSrbeW4bVxDHJ6kGSAjKdhmLl9SAH6AuCn8d4Y8nHXeaoT40qKbC5kNi4Zv/VO+RLCga7QNm8CwUMmVsyLaNCNKIqb6GoOVWr96uc7DXn0UaN1SlQbkp3/TeZLuAdsK0fa7CPAuFr6efH+moDBibGulQuueUgAgihQcrPetXy+qy1uWiev7MdYtbMWl7da1jXi1yV3jxA1CeNI2/jj3yWeG7006YoykYEt4FbFL9NsA/hhA3wqRdcOtZoQvzLumqzhvC9SHPR7ZDBGuHipHQ6u9HhAKE2u5MniXJ1aNZj+oz49eSaxUT+8Jicjs+30mcaMam4PzuWeFcfcirqUWlyTplXUtiU0TSLJp0/fTBATJMBmAZLYkYKUWzGQMUfJXpVPSV1QosAAAI4hUrlrflJ3Wypebu8WrtSOAN0/YxlDLax1VRlawzVy+FAL4UfLyOoCvDHE4m7jvS/zpXENUQlXOQgZpVPURVYeeo3HkiYSJNacMRSN72mmGxGGD1zuh19a1VLwMqhBvEsuVFkLGTjXPeklitUoSQsjopQWUABYEaRICS6CVtL1iJNsUBBkq5AUooo3iw5wcJ2ZqhWyf4g4JFABAcaWyUSh3I15tcad4ta0sB/KF3Cef0Q/yQ2DU7xw/AdbjxW4idpFmguVQ4bPzLhYDmYmyH+FKC9Ldvzaipr9o0aY5CP20rqWFWg8bIN+2vsQlN/h7eppZO4gtb+0K0csOCgyISMGOGLZUsC0BL20rxYCQDJMAJYiVLdAUiMWdkSS97ZaBC4BJKS6sVHwA4A7i1dqZ96O3/t7f/6i2rg2Jkb5rzFy+5AF4KXnZAvBVZMQ1CgDVSOGzcw3c86JMlP0IFl0o/2iFHVAG48akMLDmlCG1aNN0SD+TZ8Q+Dcv3YqezKxUy/bZw7TKegVxj+13Tra2gMCRgAmDLoJYt0DTjpDZSSZ9QQSRtQU1TUECxYNsm7g1A2gJNAqOwWguEUohMO5w79cjcXvFq7UTM6jsrrX/by/3UdMc43DF+hI24hQqwLW6sjjirND2G72GjhluQTH2jqRifm2/g3WY0YQk0iYaYgs6Iy31EutzHsJHCQCU3ATn4e5pmBOlXF4DDEFvSeFvQfb8D/YdNauUakBVPCUCGCgVfclFx7A61BJoEKMUwAAYRmBlg3jlj1hTwBZECwMXlNb+VLzdnTz82G9q5jqusv14Lbv2w4uvM0CEy8ifWzOVLdcQJCADwMNYfTtb5LoDPIM4mnUPsNr2JOO7t7iDGGDDwZwsurjXCogEEgvorEvdEJjXa9mlFouk/kgTWcmVEQou2DJA9U+wWJNBz0ZaIjwPtOyW9Mbkt/5k2/Rk8agAlPmiLK7SfJG5Wg5IYxpCRCxQKFFvc1r1JBGYiwBQUbDWwmYQgLufBKldzQzc3tdJpvFpKS6rgyqL7Z0nsuGZIjLxgS/gB4pN0AcAfAHi57bMpAH8FwC8BOAvgEwDeC+AJAI8B22vc9AMJ4M+XmnitHuQBwKDhdWrgQCJY0DXasoAigYpTRiTM/WfWHHnijMSetqjidotYGn+2i6tTbRGMBGwE9yfB/kN9FEw7A/Rr/bSR2zAw0u8z3S8GECgUJLNNbXXVFEMoZkOAlEkI0uPETA0EJNiPcvc7jVdr5zur3sxKqH60/5yafjIWgm3m8qUVAFcBNBAfz18B8GkAy8nrt4HhdyBgAH+x0sKPK56tAMsSw6vVploRguXMhPsdaRQJrDklhFq0aTogzfjrhbUtvakbST0xSnqXxp/x+v/TmmYqsfIlhWzX7x9pzFUW3KFp39VerzfN1BzGPibb5DbRRoFC3pdc5kR9GYRIEEV2bHmjTS3NhBE27fKKL/JdJ8CtBtL9xnLzxzOXLy32cp803TP0k6uHfK/t/wzgZwD+EBsXtcyokx9XfHxnpWWECkOt1SbrAcK1Q9Vo0/QITixtWrRpOiFte9Re7uEgJMkM3NZGidJsx8RyRkZbgH0SBG8MKhPzIHCSVNtrMxiBe16nrhsUYAhCZALrbknRJrCTRlUcW9ricdoCTbadZqM81RCuXDvId/KVpebNiNeT+zRDZGwE28zlS7MA3t3ydgDgIcSxbZ9FRnqPAsDr9QBfXWqKluSSLTC0sh/RmgfZ0OU+sgAToeKUERjWsIeiGRE4Lqx6IBcpJf+2ZztuFSSpQNsqzrIq1lIUyKCkmG0viOubDT9DKGLYROC0QLAQkOmvkXYVDRXnFWBYBC/MF91mcdKFQmg0gq7vf3da4epPqv4sYgOIZsiMjWBL+MYO77mIr033EHdFyAzXmyG+tNgkV/JQa7UFSy5US8eSZgEmQtUpwTcGElqpGQMOWl6CwH0tTTFsGCR6maSRlayUiGGbAoFJ8GnL/kXMlgJMIlLN4kTDyxU9ADDrQb3b8TMzf37evQngNZ1skA3GSrDNXL50H3FT+HbOAvg4gN9BnCWaKe56ET4376ISDbHALgP+ggsOM1U14MjCIFSdIjwt2jSd07XRi5Pq/P0YTBZggNLM2sPuZ7/bgnU4BgDxDx0q5HYq4REp5CISqBUmPWU7XrKAMmt+16E3b9SDuVutyMXmJD7NEBkrwZbwTWxOMPg64sKDzwH4VcTdETLFYiDx2TkXy4EcXoFdxfDnGvEztyYDEGpOES1zW+9njWYb3GU5i7S+Vz/HlBWStM6ur6lJqyfZq/6gvcBI7m1JRwezvTm8YhK+YVGzNFUl02yKJIPUcIMGKe7qwh4ojr6w4M4AuKuTDbJDJg7CXjJz+VIFGz1GgbjeWiP5/6MAPjjoMXVCJVL4zHwDc76csARcGsKTL0cKwbwu95EdCHW7oEWbZl86LWeRCBClACEzktXZT9KEiW5cv2KjGDCpjQSLoXtE00zctEWVSfAVQ6RZok0nL2u5cmiahmsZtF42yqp6XVvXfrTm3VkJVQBtXcsU43qyfg8bWaEKQK3ts8zGbLiS8afzLu60opJB8JI2IwNF+RLBYmZyMzSJaGuauWEPRJNhumkOP+4iLaW9Pyonr9stbbTx3ro1LS5XQsZOiRaDJi23kpRXAbCeZMJJ7Bo5BrkgUm5x0q2ZedM0RIsInNZnE62oKULVVfxZI1LelxbdtCPQW73dK81hGMsTN+kx+u22t0amVoKvGF9YcHGjGRYEEA2jK4J0Q4QrmamCogGhYefhWvlhD0STUTiuzbbJwtZWW23Y2mNgtIkw3tofNX2dCrP0veR7I5WUNsnKl5VaQZNYw8gi+CYQJvXvSBBFLIyoUZqqR5YTArHIY94oi2vVvNruW9iZbyw3bwUMBZ1skDnGUrAlvIy4cC4A/BkwvM4C3RIx8OXFJt6sBzkCYMTNfgc7hqqPqDbwzWp2heBaeTSswrAHoskgW+O0UutSKlJSEZPEuo0lSUweyz1cmLyDMEvrzmWVtAF8yHAUIATFoq1BprNamJTKMBUA2IZwhaDAEBQCAIUyEK2oq/vegh/VvrPqpTFrr/R2TzSHZWwF28zlSxJxwgEA3Efc+WBkUAC+udLCS1XfZkCYhIFXuA2Xm5DNMMvXsiNH08qhYWvRptmOWm9btblcR5otCaz3wRxT+MB9UUcFBRhEUFGuUFtzSoYvUUw/MwX5tkFe+gWYNb/W7ZfxZwvNG8nxcXfm8qWF3oxa0ytGxlV4QN4BMAPgAoBbAP4YwM8jtry9D0B5WAPrlO+veWhKZX3sWF6YhGbEGOjdOlxoSPtMmUTe3Dj3t1zxe3AD6Oi6QnSgRAwleh8LuKmK+g7738n+7FmJnfdYj2flACKa8N2tw9hvu9TPMg79ulludekdZjvt69qynpG/0afHDIEg2uK32j4nznAM72HISmHbfsMgLFnFAuxctSiwBJCSzIZBJJPKuQAAUizNLgvlXneDxbcaQepC1ckGGWSsBdvM5Ut84fkrXwXwf0B8Qb6eTEBcSPd3hzW2bvhpLUBLsvFrJ/KOLagRKpQGtnEGWnN10z5dglEYbgV+5u4twoIgFY/fcd40HCgbmEiuyWnm35CH1Rd6KQR3W9fIq7WERJRRLwvGZpkkWUBtFafjiCQDVaeESBgwpMpJRTbAFCnOGQTfEBQVbWMZAIy6Xyfu/FlaMqvPzbu3kpc62SCjjOUFvp2Zy5fmALy+w0dXMUIXtWtuiCuLTaMluTDw/qMMBAsNSFe3sMoSnumg6gxOu2tGA8LRyAQVsXuXjoJYCwwLq7kyIhHvqmLYkeK8VOwAgGQ4UfJ/MNiq+V0lG7xS9e/N+jKNd3tdJxtkk7E/qRO+BWxzixUwYif6TCvCFxZcUY+4lBTYHQSx8YGBYMHVfUczhm/YqDjlo5MGqNmXNOtx2OPoFxvlN47G/cu1cqg4JTBt7G5qKW63GDNDMIOMZtAg2XmhXE+q8M8W3Dttb2l3aEY5Egf8zOVLVQA/aHvrHIAPDGk4h2LOl/jsfIPWQjUxjKbxwaIWbVkjMCxUc2Xm8b1Hazok7V7QqYBP65ClImivSSQ1wcT6/+Np/3l5y2csk04LKqkpxkkttF2P4HRbSaar2KuYLe0ypUWDR6XMSdyirgTXKqAT/W0K8iSzZVW7s659d9V7ty45NWjoZIMMM3axPXvwfQAfAlAE8ACAXxvucA7OaqjwmbkGfud0cfK0Y1Z9xZOD3H6w6MJihlnWFfizQmBYXHHKmPLr7bHHmiNCm1BL64p1RBLzhs6TEbp5KKAtf7tbRzIXJ/90XB9tt/m4rQ4b4rp1vHWRrDzySGFwLVeSUpid3qSVKciXVU+KQHbszlwLpfvVpWZ7j21tXcswR0awzVy+5F94/spfAPgtAKvDHs9hqUvGZ+dd/NapwuRDeWvgoi1cagKKYU7qCvzDRhBCZhiRYapqrowpr057iLb2HtIH3+aYBXr30o3Yq6zWnSxBOxWETd9vc5OJPTbYs4zbfrNXpvQh1yvS5Ixer7sXBIaFml2EgjChOjtNOZTs1Pz71ApDWJ2fll9ZbN6UG8eETjbIOEfCJdrGKwCWANzZb8ZRoKUYn19wcd0NJh1BA3ePhistRJWRqUc8thDWb0IiEKZYzZVJkqD0prRlQjLt9nmn09jRg++EksKson3q1XqSydilaTslbkLsM23bzsC+4AwgAJXlfW4m8WqKOrw1MyNca8G7V6Nq3T9rCWqFiu1OFr3XitZ+UvXbjRc62SDjHCnBNnP5kgLwNQABgNtDHk5PCBn40mITb9b9yWHEtIWrLYRrA6/pq2ljq3iKhImKU+78oq/Z0ZqVRXYaZLt1TbM7Im4En8mTgkGo2aWkk8nuP6WgDbGuWiG8ezWg4jXN+PndrEXqTMT7CzZm5s8tNG5sOZ60OzTjZPLg7TM3EBfR/TTG5ABVAL663MJLVX+y19mjncTDRGsewlUt2obFTpf3SBhY06Jt7DBi0bHu80oC6tVRyZg8KFu/tywhSWAtV4Zn7m8YE0QBS4VgyUUw34gQKhUxCgVDLBkEX4Aig2hfK9lbjXD+ZjNqL6x7TycbZJ8jd5LPXL7EiK1sDQBfAvDScEfUO7676uG7q97EQUSbSWjZAlVboGYLVCl5khOEjkzkUcXTDeMzhhQG1nITkONfAP5IkPYHbX9vW9S8BkBsdWzPSM2qZS0wLKzlJhCJDsPJG0Hg3a1B1gOYBN8xaK1k0n1boPlgzvxpyaQVT3J5L7doqFh+br7x7pa3x8J4Me5k8iDuNzOXL80DeDV5+S1gYDXN+s5LVR9fW2pNWLFo6/g6TgQZKEwGChOBwqRFcE1CM1Cdt++Kqj7CZS3askT69B5p0bYnWXcp7pbkkcQv6h8XG6U/gI3kDAUy0jjCYY9vK02z83g1DiX82TrCJTcqEBYEITpmGbfOO8ZbZxzzxpRlLBARCoaoFU2quFJN77auH1e8Oyuhaq/N5EMnG4wER1KwJXwLQAigBeCfo489FgfNm40AX1hwJwxCHR2Ktq2uz0BhIulb2tWFLqr5CJbcOKpGkwkUCVTaqqRrtpPFG3o7SeYnb5k6ClkYd+L2VBsN77McqxaTxKvZHVxeN5IKoLwIkjk3bYm7Jyzj7aZU0w3J03db0XN3W9H7A8W5eqSOhYodi8iTzNtOeDdS/pcWm3e3vP3azOVLurjmCJDhg7q/zFy+VEdcmw0YwyeMW80In5lvTACdibZexsDIeoBgqalFW4ZQJLDmlBF26no5YjBAWVZsu2R4IhFuR1q0tTd+V4CR1Vg1ILZ4r3YYr5YmFURr3voVPFQo3vfkzy0G8r31iB9Y8OUzvuLJkLnYlGpyMZAXFwP5ZMjseJK39a37xkrrlqe2dUF4pRf7puk/R1awJfwAWO/L+WcAZoY3lN5z35P4k7n6RMTc2O+i3uuLvmwECBa1pS1LMAlUtGjblVETPrE7lI9E39CdSLolZNyatkEgOotXY6kQLLrw5xrgsLNDUgARA4IZZqAwUTbFqinIZ+b155AlX9a/vdLamlhwLwkR0owAI3Gg94vEDPwXycsIwI+GOJy+sBQo/NvZRtmVyt3nhtRzA4N0Q/gLWrRlCSZCxSkjMKxhDyWLjMyBaiSxWqll6aiR9hLdq0VVlmiaOVRy+8erRXUf3t1a1+3/GDACxQUAMAgeM8gR5IWM9XY0X1x0t5bxAHSywUhxpAVbwqsA0qeOdwD8KwCVYQ2mH1QjhX832yivRdKlnYtuIopP7J7fsFQzhD/f0KItQzARqnYJvhZtm8j8Xb+NxIU7UhbBXpD2IJVdtKoaLoSaXdw3Xo2DNKkg7iDTLZJhu5JPGwRfECJTxKU90iSMG2649Ho92Fqnc+xCgcadIy/Y2orpArFgeRfAleGNqD80FePfzzbK837UJNou2pJMs77cs1Qrgj/XONCFSNMfmAg1pwTf6KgouiZDpA3MR8UVeBgEoJJm8irpUiBGZb834tX26LnMjHC1Be9+nFSwEwbBb3+dM2i1aNA8Nj9gEwA65Rg/eyBnvra+eoAks/r8QuPWDqt+XScbjBYjceD3m5nLl24iLqibcgNxYd17wxlRfwgY+NN5t3y7GbrbRBuDdrO+9QLlRfDn6mAt2jIDg1B1invfUI4Qo3BkpmJtnPq47oZI3J5xeY71Nl8jQSfxarIVwrtbi9v77bFjDIiiQQsFgxYFISSAQ+a8IGxSeBOmuFs0RM0gWn9fEOSrteD+PU/uVNlcu0NHDC3YNvgaNp82PwPw7wHs/NgzokgAX1xsTlxrBJtEGwPCFtTo57aVLxHM1eMEfE1GINTsAlpatAEj4hUdhZitXsAjJNDa2S9ejaVCsNBAMNcAR/tfCxXDciWfDhSXj1ni5pQp7p60jRunbeOtokHzloBbMun+cUts65Hdkqy+sODu1Ibxvk42GD20YEuYuXxpEdvTm09gDL8jBvD15dbEK1W/sVm0cd9vBMqXcfaTFm0ZglC3C2iauWEPRLM/Y+8KJWS/SfvO7B+vFtWSpAK3+x7rESO/EqgnFwP59FqoHiiaonraMa4/nLdeOuOYN2hLYRpfce4Ha96tWqR2MjqMTYefo8RYn/gH4C8QN4ZPWQOwBHTWnmnU+P6aN/m91VY9FW2KMZAo9DTAtpOnS82gIDTsPFwrP+yBaPZAAYYY42SDjezXg9+bEsGXtqVKJ5VObUWHe4baJ15NBRL+/VrcCeYQYSEMCIuoVTbFIgAIom3HgmIWkWLruhu4X19ubrO6QScbjCxasLUxc/lSA8D32t6qAPgDAP8IgDeEIfWdn9aCqa8vNasMlgbRwAJQOVRxTJsWbRmC4Fp5NI6uaMu8B87oY5xpVjjoj5AmKABpAV0y0kK6aQxc2qYqKUSMRLhxLOZYHkTEBcLC6m7xaooRrjTh36tB+Qf76SzCepN2QYgk2LIFbesBGCq2kx6i6ocV7+0/ult/MeQdv06dbDCiaMG2nR9ic29RBaCJMYtla+eaGx67stisBDzY2hscqtjSFo79PWikaFr5xK2jyRJGUiS2Vy7RHVpdrVueem2B6pRu3KCpJc0AJMWFg0S35T7SDhLxd0qE5HvodPm94tVkM+lUUPV3WLJzIsCZssS7ydjkMct4VwDKjdRkoNgBAF9x3iCEiuH+m9nGa5+dd+/ssRM62WBEoQHfo0eCC89feQ7AX297qwzgvxjKYAbIKduY++TJwtmyOVgdTwbBOVcGWb1PfDMF+VFyURtXDEIo++DOzkc+yoG742dJhfmxylSM+1FmtxCtAUQKOFSbivRqn4qdTrJN2xRUVzeLnYRPD+LSKB2I0eH4u6Wz4yBO1tnJBcqRQrjSPFCc2m44gtbyBlUcQY1apM4Gikt5g1YMonDKFLMGIVoJVe1/uVt7a86Xe3mD7s9cvvSHPRuYZqDoHjU78xqADwM4m7z2ACwjTkIYWxYDGXx6roFPnirgrDO4Q4Mlw5utwzlbhrAze788crRMBwxgYhfRNm5QRmPcBVgCtN4/tBcwYrdhYpnac8fbttnVF5QU9913/QehX2ItZu+2sooEKk5puwuUGVHNR7jm9bzmpK942lc8jTbNHSieOGmLd0xCdNUN5/7Vvfo7O/QJ3Yq2ro0w2iW6AzOXLzE2iukCcdLBUfDbnW4q9j8/7+JnXbZGOTSS4c/Wofyx9TyPJJ7poOps6yGtGRCxm5KQuEJ7KlB4fRP9ox/dGOJm7/2z7jIg2pMU2idpmLKSnwjYMAODEJiCfABQfhR3Klhp9VSsbf3+bIGGJeCahOZJ23jbElT/6nLr7T+8U7vWgVjzAbzZs8FpBo52ie7Bheev/McAnkxe/i0AF4c4nEGxbkl8tmzjY8dyg7U8CIJzpgSR642FT7tEe4MjA0z6LtLb/Di6RPtrtTkYIo7PQj/GFTeP7y/9+E6H9Ts1zRwadh6bNK5ioOo1Wmtez55qBCGyiFxf8aQtUCUQFwxatoh8BTYEUZSMYPHfztZf+1kjrO29xnVemrl86Uu9Gqdm8GiX6N58DcDjiC2R54Y8lkFxAsAqgGOv1QOshgq/ebIARwxItCmGP5e4R3sk2jSHxzdsVBzClN/ACCRTjg2Ew5W42GfdikfQyyIBI0nAGFA/0Z3j1aQbIFxpgSNVcgRVfcWTvdkaFBEik9Cctow7eYNqZlv3AgC43Qq9f3Wv/spKqLpxhejaayPOyJ2sg2Tm8qVlbPj8v4Yxawq/B+tBS3e9CJ+ea2BtkJmcDPhzdajWWJa/G1kCw8KaUwKDwEk9sG4y6jTdk2Yx9nH1faVfY5dx/J3q92PkTvXVOFII5hsIFtz1skS+4klH0Nbm6gdCMqycoOr5nPla2RSr7WKNmfkHa97MP3m3+uMuxZrubDAGaBPG/nwbwPsRt6pqAngYwAriuLYJAA8idpUayXvj0E27hHgfjwNAJVL49FwDf/VkAQ/nB1JbNxZt8w3Yp0swCgPapmZfQsNCJVfGlFcHgYWIYzvXLR1pg26t4npDPwXJIGzm/bTgKUAkNdQ4Scjo6S6FwkTVaSvZkSYVrLZ2lLqpaDuMpc0WqB+3jRtFQ2xzc7qR8j8937j6ai2oHGDVOtlgDNAxbB1w4fkrHwPw63vMYgAoAGgAOAbgrwB4agBD6ychgEUA59M3CMBHp3P4wMQAQ8IIsE8VYRQPpoN1DFt/MFWEKb8OwYy48j4zg0RqcRtF0ZalGLa0mGta6LUf2xjU/vYrU7SdXu/L1ng15UcIlprgYH9Pw0FFW96glVO28Y4lthcwv9UMl//VvfrblUgdxO3gA/iHulju6KNdop3xYwB7mbslgDrim9UK4qbxf47RbiFjIbawre8DA/j+modvLDchByX0GQgWXMhBZ61q9iQSJirOBBSJJIORjLQAaWL5OApZ1X0jzYQcvX6a2+lHpuhWkmOuB2zpB6oY4XIT/v16R2INiC1tttjzfrENW6B22jGubRVrEbP62lLznX86U33zgGINAN7QYm080IKtA2YuXwoBfLOLRRRikXerPyMaCE3EHR+2HSPX3BCfW3DRHGAD92DRhawfrmK4prdEwsCaU96xyjsAMtqq0A96bKNMkoHbdwvlzl2Les+AXK902IeErfFq0g3g3a0iqnV/3QkUJm2BTrI3GQAft4xbWxML1kLp/k+3qy9/eak5e8hfSrtDxwQt2DrnDQCzXS7zZ4ib7I6ateFu8vfUbjPM+xL/fq6BpQ6fOntBsNQ80MVT0z+kMLCWK0NuKQyvkjZBiZuKkvIUnAi49b6PR6E3ZjfQjv/t29ZG3nq3hQPvUCjM9X6gHCn4aVKBPLhUChTKlkBjr3kMQlAyabZoik0Wuddr/ux/d7Py8s1mdNiq1bMzly/NHXIdmoygY9i64MLzVx4F8HsAbgA4iThp4zaAv434xvMmYhH8s+R1+uUWAbwXwKMAlhDHFLwPwOkBDr9TFOIs0XInM5sE/NrxPJ44YIzZQbCO52FO5jqa94jEsAWSh5vsIlhhyqvD5P31V3v9r8T61rMK/ochCzFsaezauMSWAevu3YEYBw7yG67HqzEQVX2EazsnFRwEApRJaIWM4k6fC0J02jbeTAWbrzj64oL79vfWvKXejABfnLl8SVvYxgQt2Lqkrc/ofwLgDID/FfHF6DSAdxG7ETv5UssA/q/ISJBzG7cRZ8J2xc9NOviFqc5EVC+wjuVhdrA9LdgGh2CFKb8OU3VvNEtv6qnFbacaWwKQ3GUQ/lZBaGy4G7etY9iCLWm+3tcq/lu2NxCn6KCEYbotArgzgbhRX015EYLlzpIKuh4TQRpAEDHyQCzSEGfpGAWDFk45xi2DKJrzouq/ul+/uk8v0G4IAPyejl8bH7RgOwAXnr/yiwB+E/FF6DC9lH4XwHt6MqjecQfAQwdZ8JG8id84WYA1IE+LOZ2DNZ3fex4t2AYKscKU34ClujstUqGW1NdaR2wIOLE5eJ3BoE3CKxU87UJExfXiZPpe2j9zJ2E07O4Ng7Q2DlJEDXp7nTS2VyRQdUoIYCBabfU91MISaBQNsVg0aDUnqElELJmNluRy0aC1H1X825+Za9yWvdXQurPBmKFj2A7AzOVL3wfwIg4n1oA4EzNCnFmaFeV84HG824rwmbkGatFgkhGiNS+uiaTJDEwCFaeMcGtj7H1oi3dLi8Wuvw9sZE1uTHHQXCr00iK+KllPOh8QizTZ9jqp2WUIcKYSIhKxOZCYvkHv9yC3x9johrDT52m8mtdUB04q6BQC1HFbvPNQznzlhG3czRvCJYo1uUEkCVj4l/frr/3JXGOmx2IN0MkGY4cunHtwvo64wOz7D7GOP0Xc8uo2gCcQx8KNNCthXGT3kycLODeA1lJRxQMUwzpR6Pu2NJ3BRKg4ZUz6DdgHrkSw2UqGHW5mDFB6Y+7WlRmLQjLSpt6DirHajVg8Ut/aUO3A2Aq2lMRau8my1zId1CiHYN6Fah32eXtvCJBnc8ZrBUPUd/r83Wa4/C8PXlttP3SywRiiLWwHZObyJQbwBQDXD7EaD3HpDwngGoB/jbhY7TA5tGutpRifX3Dx1oBqp0U1H8FycyDbyiiZsRKlMBGqTgm+cfh6vvu1ZzpM3FmazZrGPgHxhgYtMBhk9LNI7lbGLT10N2i9fRqhZhWx1jTg3av1XawBQMmkuZ3EWsSsvr7UvP5PDldbbT+0dW0M0YLtEMxcviQBfBrAvR6t8iaAPwLwdo/W1y0eOswO3Q8F4C9WWvjOagtqAPc+WfMRLLnx87QmEzARak4JvpGJ8Lo9SYSbGXdtiA/YQZYcSRIqxlZHDeusVIAhSdAa5cLKwu5tpXqNLVA/aRvvbn2/Esrm79+uvXJlqXm/j8MIEFcs0IwZWrAdkiQD548BLO8x21Yzw2kAH0Dci3QrAYDPA3ipF+PrggBxyZFSL1f6Rj3Any248FT/49pkPUCwqEVblmAQqk4RnjEaeR+plWvQ4onBAxWIg6afPUX3ImSBhbpAY94DB4coqtYFBKhTtnlNEG266L1ZD+b+u5uVl683wz1rs/WAN2YuX9IFK8cQLdh6wMzlS03E5T12q2y91ey9gPg57+8j7kG6FQ/AFQDf79UYO2AecSP7nnPPk/j0nIvVsP/3I+mGCBa0aMsWhJpTQMscDdE2DLZmvPZ/e4M9P3gIVXq9psLcMiNoyHT7fd9pg+BPWeLdnEHrMRq+4uhP5xtv/Yu7tbdbqoNChYdHu0PHFC3YesTM5UtVxKLtSXRmpXoNwO9j74vINwC8fvjRdURfj4VqFCcjzLT6FbKxgWyG8I+SaBsJRxqhbhfQNAdXq2+USAXF4H7KwXc5GERPUQAQQYTaoo+FmoFUHjHDtAQO2zVgX4qGWMgbtP7gPu9H1X/0buWl7672rBDufszOXL7UbUcezYigBVsPmbl8aQnA/wvAVezv3mDE/Tr3q0vxFcSN5fvNVL83EDJwZbGJVwbQXko1Q/jzDbA6KqptFCA07Dxca+/aeUeV9nImA9jWEGR+f89FYoaxXMP8ksSa2v5gECiUu23K3i31SD1gE7WYmX+05t3+h7cqr/awEG4naOvaGKML5/aJC89fKQP4KICfw/YYtm55DsDvHHZM+zCLuMTIQHiyaOFXj+dh9vlBP1ewanSqOAExEmaoAzGKxYELYQulMNs19AbZ+SBJdhhYWY9BdTlop5/fp+m2IJcauGcfQ7RHDUBBiEyCGyhMoMei1SD4J2zjbUGY/8xc4+pPa0Gll+vvgADAP9Txa+OLFmx95sLzV4oAfgHAfwAgvameQRwz1ik5AP8l+ls3LwJwH7G4HIhwO2UbuHSqgKLRv3uUI6gaWmLSPlOKq26NIaMo2ACgEHoohdktx2LE50RHN3Xe4X87LLrvugYh2IYh1oD+CDaKJJyFFdRaAvOlE+AONZgl4EYK+V4kQ+QMWjUJftkQC4uBvNnH2mr78fLM5UtfHMJ2NQNCC7YBceH5KznEou0jiAXb4+jOfP2fYTBCagXA8QFsBwBQNAiXThVwyu6PFnUEVX3Fk2QbcM6Op2gbVcEGAPnIRznoe2jRgRikha29NVe/tzWsQsG93q61Voc9v4K5/HGs5XZKuN9nPIRQ8aG9H3ggZ7xkCar/xUrr5pXFvpbr2I8/0PFr440WbAPmwvNXbMRu0v8WcWzaqx0u+g8AnO/TsNqpYADxbO2YBPzq8TyeLPa+Xlcq2ABgXEXbKAs2AMhFPiYyKNoG3Qx+UNsbVs/U1PZ12DuO8EPkZpfArRC3J86gdcBElvZrw0EwCc0JU9wnws0/vt94awDlOvZibubypd8f4vY1A0C3phowSd22H1x4/spvIq7Fdh5xW6r9eistob+CLY15HmSALAAgYuBryy0sBwofmXZAfYqH5kDCn63DOVsGmeMl2kYZz3QAECaCYd7vhs+gRNSwojkZ6wWCD7SfxAx7aQ32UgUtM4fbUw/sGa+2H5K5qydEApQtqCqZcwVDLOUNqt5uRW/9m/v1GwMq17EXOtngCKAF25CYuXwpBPCXAHDh+SuE2N35RDKdw/br6k3ErsqD1kprAnABNJLJ3TL9HIDHELtrh8IrNR9rocRvnCzA6lMyAocqFm3ntGjLEp5pg6mESd9FBjttDQQ6hJgZd0y3Bef+EowgxEpuEnNdxKvthuoyGSxv0HLIXDCIfALq31hufu87q96wWwkCcbLBG8MehKb/aJdoBkkSFR5HLN4eA5DWQfgQgN/aMnsTsWu1lvxNp0bb3wb2roH0EIBP9Wj4h+aYJXDpVBGTPRBUu7k9yBSxe9Qa/fvjqLtE27FliCm/gSyItmG4RBn9t7QNer/aSYrXUqe/LkUSufkVmJU6AML90skDxavtMhYlCKFk7HnuWAINk8gjgAUhbEm+8em5xg9nB1uuYy90ssERQVvYMsjM5Usu4sK6r114/orAhtvUAeADqGJDmPXCFN+bK2CPWA3jIru/ebKAB3L9OUQ5UhvuUXv0Rdu4EBgW1pwSpvwGKAOibQj03WM5zG+VAepUMFprNTjzKxBSIRImbk+c6WnhZQaEZDi2QDVQWH+oswRck6jFDCEIkUEIAoXSKdt4+9W6//Zn5hoz4aBbReyNdoceEbSFbcRI6ru1W996ZVl5EHG9t/cByES3bgHg48dyeKZ88F3cN7DYIDhnyxAjLNrGycKWYqkIU159qKJt0JYoASgCuN/bHFbSQdv298wWNfwAzuwSTDc2YDWtPG6XTx8qXm0viBAxbxgvigbN+4onDaLAEnBP2sa7nuTmZ+fdq6/U/LW+DOLgzCHODtU38iOAFmwjzIXnrxiIhVYa+3aqB6u1ALwHsXi70IP1HZr3lWz80vEcxAGMDx1lggmCc7YE4YymwXkcBRsAmCrClF+HGNI1ahiuQxFbzKmfZTeGVdZjv+2TUrCXKrCXK6DkN1/NTWK2B/Fq+2EL1JJiulw0aIEB8hVPFgxarkXqzT+ebbyyFg6lttp+fGnm8qWXhj0IzWDQgm2MuPD8lUlsWN8exeEtZdOIM1mfA1A+5LoOxXnHwG+eKiAvurvPdJy6LwjOmRJEn1yw/WRcBRsAmEomom0gbSg3MSTBpjiO8eqbQhlW4dx2tn63ZqOJ3OwyRBBrIu5xvNp+tF8nTELrmGXcCpjNv6x4b359ufXmsL+vXdCdDY4YWrCNKReev2IiTiZIrW8nDrE6QiwEPwDgKQzp6XzCJFw6VcTxLhIFuqq1RIhFW/7QtTQHyjgLNgAwlMT0EETbIAUbYaM5ej+tX1kQa2lpDwYgIglnbhlmtbGuUPsRr7YfbXFsXDJpdiVQ5ndWW//6HXeotdX245WZy5f+bNiD0AwOLdiOCBeevzKNDfF2AQfvb1oE8CyAD2KAHRFSLAL+yokCHi10Nvyui2MSYJ8pwRgh0Tbugg0ADJaY8uowBijahpQl2ld3KMUd2IfaWDd1idqrNTgLKyC58Zs2zRxuT5zpW7zabhgEP54omGmG17602PxeS3HWLVd/OHP50v1hD0IzOLRgO4JceP6KBeBhxOLtScSuz4PwEGKr23tx+Ab3XfHhKQc/P7n/E/iBqpkTYJ8uwdhHFBIgBW0rl8Lo7Ia434nX8U2ViNS4CzYAEKww7dVhDKhG6Th2Ohh2/BoAmH7A9v0lMpubq2J0Eq8mmSEZiNr+Rrz5/fR1+/th8ldu+bx9PQBuB4qXlwL1fQBZSy7YyjyA39fJBkcLLdiOOEnR3mPYbH3r9qbhIBZtH8Rg2mcBAJ4oWPi1E3mYexTZPUz7Gft0EcYe7bIMAV+qnmXpajpAsMKUX4ep+i/aUjflIATOoDI3e7UdZt5hSt7HDu8xg5UC3Z0LgnsLtpIKihWUUogYeO34I7hVOrOjkGoXWn28W0nEbfm+BeBt9KZcUj/RyQZHEC3YNJtIep0+gjhm7UkA3Yqdk4itbs9i/3Zbh+aUbeA/PFVAaZf+oIftF2ifKsIo7SzaDEIgORslUI4SghUm/QYsFfV9W73qf7kXB84MbRdFiUhSbQIKiVhSiYgCMwgIpWJrXURhB2HV9pna4bODECyvoXn1Jky3UQ2EuX4+tkwH3z3/HJYLUwdabw95HcBVALcQB/NnmRDA7+lkg6OHFmyaXUmsbyewYX17GJ3fVAzEgu+DiMVf3ygahE+eLODMDmU59hNseYOWSwYtbX2f2/7SyWKZS/aO/ldmGIFUJcmDdQkfdYgZU359IKJtP1flrpamvQQRUlHFADiSCub2z3cWUen7B0EQ+4ppYFZh5Qdovn0LwfwyAMBW0bpgW8pP4cXzz6JlDS65YBdqABYAfANAFlpN7YdONjiiaMGm6ZgLz19xEJcLSQVcp6U+JhGXBvkg+tRVwQDwK8fzeHqLNWw/wTZpitsnHWNmr3UzgOBE4bgs2TvuLwOoe9GZ/VrcaHoLMWPSb8BWYcfWptRqhD1E1Nb3CRxIZpsZm9bZi2unAAI1wELVg9oeM8O/N4/W9RlwtOFdtFVUC4Q5cWPqAfzk9NNQIhMFq10AXwbws2EPpEN0ssERZfSKTmmGRmKCvwrgamJ9O42Num8PYnfrWxXAdwB8N5n/g4jLg/QsW00C+MZKCyuhxEenc6AeJsIRAHu5uRIyc1R2tglOApCzjKobyF4ULh5pFDMUYsORanutkvij3V4zxymSsaBC2zo2XqstnysGwCYeqK1gInD7tk8GgWUfnmsFcUigbDU56gFR3YX7sxuQ1fq2zzzDLr98+in/nemHsvRwU8SQ60x2wTyA2WEPQjMctIVN0xMuPH8ljw3r2+MASvssUkIc6/ZBAFO9HMvDeRN/9UQBtqCeWNhSGEB4LD8VTThTO33W8OXpSPHQ/TvroqYTkYQkzmmHzzeE0x6vE1GVCrRhQGA8WF9Mmsb3Y/1QBES9tkwJIGDA5AFmbQriUDH1xX3PUqJ18w6827OJq3cznungu+efxVJheg0Hz0zvNQGANwB8FXFsWNa5MnP50k+GPQjNcNCCTdNzEuvbWWy4Ts9jd2saIRZ6HwLw9B7zdcW0JfBbpwo4ZZs9E2xAItqmc5PRZG7bDUcxG26gjkeK88nrbdYmRmpJ2t0atel1F9am9PVRPKMJjAfqS5j2t1t1ekE/XIkGsS8HGE8GACIWnj33rARLq2hevQnl7RwHv5yfwnfjeDWJIfYx3YGfAmgB+PqwB9IBIeLOBt6+c2rGEi3YNH3nwvNXCogb1afWt92yR1OrW5qkwG0TsF3M7XTwrif25QTRb54sKEdQOZmZty5WNMTdsilmImYVMVSoWEYMFTGrSLEK4/dUyJv/f/J8+dTJM+WHFMDrE8fRUqFka6Xu/zxzpm5MR4LzjSUc92p9WTfFurhn1jCD4A867jH27lLPjkvlB3Cv3US4sLLrPFvi1SrosUX9kLwK4GuIRVvW+enM5UtfGPYgNMNDCzbNQLnw/BUB4Bw2rG/nhjuig/O+x6bPP/3w5BM7fRZFnK80/A8y6zjRQXPWXcHJVqXn6+2lwBLgACAMMuEA6J3oZGb4d+fQvH47zqHdaR4i/OT0RVyfegDYqJU4i+yc8y0A7yK2rlWGO5SO0MkGRxwt2DRD5cLzV0rYSFx4DMDQY8C64T0Xps6+55GpJ3fKcvB8eaLRCi8yD7ey/FHkTHMVp5q9LVYvgCgx9x5ahA86O7SXRPUG3LduQNZ2jxlsi1fb+tEigKwk5/wMcULU15H9SALd2UCjBZsmOyTWtwewYX07M9wRdcZTD02efuax6ad3Em1RxPmaG7xPKu57EWHNZk4113CmudrTdfZKaPUz+L9fcCTRunkb3p3ZPeXNclJfrbl7fbX7GGBHlD14FcCXkP2uBoBONtBAl/XQZIiZy5cUgDvJ9M0Lz1+ZwIb17VEgm3XO3r5TXZCK1bNPHHsPbRFtpkmtQs6cqTfD9wxrfEeVxcI0FAmcc5eHPZRtMJMhiH3wHn3VNi8R/6E2qcQggMQBkggYeyT3CCDYqsiixRXZ/NmNiD0fFkAg0E6a7ebk+eiHZ9/nyzhejXeYgOzEi81jNMRaiDiTVXPE0RY2zUhw4fkrBuJm82niQlbcKus8eq584gNPHn8Pie0u0KYXnXVb0RPoYe05TWcc82o431jqyRefNVfmQeLq9opjoziJeT0pQXk+3Ks3ES7tbalkIrx0+mm8M/Vge7zabniItzFsC+NPAYxCxwCdbKABoAWbZkS58PyVKWy2vg374g8AePhM6diHnj7+PiFoR9HW9KJHdSLC4Jn263igvnho0dbrTNHDcpCsz72WSV21zAz/9iyaN+/smlSQ4pk2Xjz/LBYLx7oZxqDdoh7ih6V2cXsbwP9vgGM4KP9i5vKle8MehGb4aMGmGXkuPH/FRNznNI19Oz7M8Tx4qjj18+858YwQ22+KzBBNLzrf9KJHhzG2o8yk38BD9YUDi7asWdcOOp6ttdi2idBavVV960Ze1vfvHrGSn8R3zz+3V7zabtxGfM4OgghxFmgBwBLizgaTAD4D4J0BjeGgLAD4n3WygQbQgk0zhlx4/soxbLhOH8EQYjXPnShMfvi9J58xDNpx234gj9eb0VPMnAnL4FFhInDxcG1hUyhYpwy6cfp+HLTMiCCEimOLtAAiS3DDVzTFUYTmjdvw787BlhtN2nfj5tR5/OXpi4fpB7qGOD5LJpPC5uA5wobr1EJsHcuhewtnFbFAa0chdol+qetRD5Yvz1y+9JfDHoQmG2jBphlrLjx/xQJwARvWt4G1xDlzLF/+yDOnnt1NtHm+PFFvhu8d1Hg0MeWwiYdr8xBdXPuymNV50K4F7fviCK44guvL91cfdK/dAvsBAMBk6UZkFAlgAksArBAv02W8Wj8IAPjJlAq+drEnsCH0TiTvudjeLu9VAFmODdOdDTSb0IJNc2RIWmYdx4Z4exh9bpNzajpX+ugzp541TbHjzX6tFjwXSbWnJUPTe4phC4/U5iG4sw6ow+hK0AkHaW/V7krNB61ZvPr64v1K+Fz7PLaKLWwEqHLQnJFCOK6ZO78lXq0KYAJxfJiNbLWcSlkFcAxAHdsbvL+F2C2aVXSygWYTOvhZc2RI4kCWk+mHF56/YiNOWEiTF3ounBbXvMb3Xlt49Rfff/pZyxLb4o0Mg9xI9n67mr1xrTxuTZ7DI9VZGB2KtiwSizXmnUxdBMg0wYABI41VI4JixfBu30fujdfuGmEQiuKJQBGtH5+COTJZtQisCpG3okhY84XjxW8++EHPtQsKcWzVW4jFmkBsvfpdZC8Luol4bFvFGhCf9wayW9rj5WEPQJMttIVNo8G69e0kNqxvD6GH2YDTZTv/8efOPGdbYpM1JJKcr9T8n8tS5uFRIh/5eKQ2B1Ptfc/OWsJBO+3WP5u4JgHHJjQsAZfArEBmqFAQQBQyirlaZS56+bVKxY1OH/drt4iZPcMuLRamnwMAS0UNRwZrERn5qaBxW5Jh3y6fvvPFR3/xWmhYEWJXHRC7HIsAHgTwJICL2G5lG3az9ypi0XZ2l88vI3atZg2dbKDZhhZsGs0OXHj+Sg6x9S1NXtjpCb0rJkt27peeO/2cYxubUuqSGm1PHnb9moORkwEeqc7C2kO0CXCgQJkWbILYLxm8uG5Ji+O6AAAMGAijCG9dXeV3btYZIAYJRSQkCZuJRNPMHVdEVilozVkq8gxWkSLiH515743vPPDcfd5eqWYrZQC/gLhG4ruIRdwZZLTgdcLvIy6gmzV0soFmG1qwaTT7kFjfTmPD+vYgDuj6mShYzi994MxzOcfIr7/JoOWK93HOnjvpyGDLEI9WZ2GraMfPs2xhI7C0BWoWoWUQhzvNw/fnXH7ltVXV8pgYTNj/wu8ZdvDlR37hZ9eOXah0OaSzAN6PWLB1VZxtCHwJ2XM96mQDzY5owabRdMmF56/kETeqT61vxW6WL+VN+5c+cObZQs6Ml2PQUsX7OLRgGyq2jPBIbRaO3K55slbSI4UAVTLVHMWuR2IGMSAMStyWzVaoXnltFffnumoHtZqbqH/miU+8uZyfOqi70ALwFGIL20fQmctfIq6LlroxLyC2cveTl5G90h6vzly+9PlhD0KTPbRg02gOQWJ9O4cN69s5dCC8CjnT+uUPnHm2mDdLzBArFe9j2sI2fCwV4ZHqHHIy2PR+Fst62MQ1R6AmiKWnaFIQIovQZGYhwBFfv1lVr1+tQUnu5kJ/c/Lc/Ocf+6V3fNPuRTbG+xG7SXeLIWvnLQCfxUZ5DgPAJwB8GP3rZHIfwL/o07oPyh/NXL50d9iD0GQPLdg0mh5y4fkrRWy2vuV3mzdnG+YnPnjm2VLBKq9W/Z+XiguDGqdmd0wl8UhtFvkoFm1Zta4J4tAm1BzBjWokHrIF1x1CndZWXX7pp8tYq4YAIEmYBu/i622jy3i1bngAwH+M3c8Fibi5uQvgW2iLvUuYBPC3ED8M9RoJ4PcQZ7tmAZ1soNkVLdg0mj5x4fkrAnG/xNT6ts3K4NjC/OUPnH3GscSFmhu+v9djKOWta4rZavnRBTBEYsVjHGFrHhEk896ZiwYrPFKdRTHyJOK+CIfOdDQIPoPRa/EngIABw5Bh07j69n3rnXcW0ws7A2AiIXjv2iW+YYVXHvnIWweIV+uEIoD3AvjkDp/dA/AVALP7rGMKwP8J/bG0/Vtkp0WVTjbQ7IoWbBrNgLjw/JUyNmq+PYYke84yhfGJD555Rin+lUjyxGG3YxpUt01j2TColbONZcS1HdYFh2K2glBNESFiQLS86CGpuKs4vFHGNERVKlUEwyjmrbdbfvTQTtZNgxUeqc1FhbBlCEJgEgKRlLSgWAspCeRCRcW9yrLYgqsGIbCJmy1FxxTDUAyzV1mnJrEXzC5Gxhtv3rGbbsOW4Xq8GhPRfhf5tVy5/uknfuUw8Wr74SDuOPAPktcSwHcBtBBb1tJuBfvx1wE824fxfQPA9/uw3m6JAPyeTjbQ7IYWbBrNELjw/BUDsavoCQBPmAad/dDTJ34nZxuHsrIRUXh80vkR0Ta30q74gZxuNKOnFXMmsyD7QSFn3gLAxbx5r1IP3hdG6nj754IosG2xYIH9C2v37amgseeF0ldU9hVN7mSJM4lbRYMXt8w/4Sk6dJs02WwBb16bK8zeucsgcmTY7Gb5W5PnFj732C+93aN4tb04i7jcxwSAHyK2lBFi0dbpTYgAfAzAr/Z4bFmJY9PJBpo90YJNo8kAF56/Mpl3jKfe++j0Py7mzF846HoKjvlusWDe6Xa5lhedabSipw663VGjlDffzufMeaXYXK35vyCIfKm4IAR5ece4U3DM+fYO8VYUGsVmPVfw3Fzeb+asKNjmmmOAQkWFiJGPmBxOOslYxA1LoOlJmhLEKid41SCETUknQqauLJvrnQuUgjdzH/LGu5Xj7trbloq8rT5uBtFu5TsUEf/4zHtufvuBD9zrcbzafpiIy6MEiC1K3VIE8F/1dEQx/whxduow0ckGmj3Rgk2jyRB/709epmLe/L9Jyf+FF8jzXS7OJyZz3yfRfasdpdhyW9FDyTbHPr7NtsTSRMF+WzFbDFAUqZIQ5NumqLcLtV2XDwOz2IoFXMFzc4aMtrX5U0wGwBAU/x6SyfYUphyBikkcMGAEivIRoxgx5bZvZRtcMtRca6Waq7x5syDdllBE9kTg3poMmrOCNyr/7hWk6BtW+OULH3nr6vG+xKv1GwHgPwHwSI/X+y0AL/Z4nd2wCOB/0skGmr3Qgk2jySD/5y++/jcbzfBP0EW/X8OgxrEJ51BFQMNQlSuN4IOHWceoUCpYV21TVMNITeRsY+mgMpUYsEPPKrYasQXOc3OGkh0nKUgm21eY2M/apsIQrXdmZO7dW++WwtZiKMyCa+VORWTki5E3Vwi9tf2yQddy5fpnHv/EW0uF6VGOkzIB/BriciG9YgXA/9jD9XXLV2YuX/rxELevGQF083eNJoP8s99+/5/+b//1T37FNOhPhKCOyhnkbeP+YbdrWaJeylvXXC98Yr9MypGHIYJQTudj9+fBV0OAb+dC386Fq5Mn6sSAE7SsYquRL3iNXN5v5YSSe/kdOWkvZW3tpiBit6EKZhdQf2fGNj2vZbAKBbOyVNSyZNTIs7/iGfZUpXTiiVwULJfC1kI+8re59wYYr9ZvIsSJAu9BHBPXC6YQx9gt7jNfP4gAvD6E7WpGDG1h02gyzD9+dSbX9KK/t1YPfndhtfVre807VbJftizR6MV2g1BN1tzwvcycqWKxvaKYM68zYORsY9EwqK/Nv4kZOb9lF1v1fBwD13KI1bqAi5gcX2ECgIjikh8ExGKt6DXu8iuvLkYLK7xYmHq/JCN3zl3+EYMMS0UeACgShiQy54onfj4ng5VQGKVpr3G9EHmV+HPivzz9npt/8eDA49X6zSkAfx+9axn2QwBf69G6uuG1mcuXPjeE7WpGDC3YNJqM88L1WQLwwet3a/9NtRG8x22FT24tI3HuZOHLj54r/7TeDCcbzfBY05PHI6U6iYvalUYzvNDy5cOHGnzGEIK8Ut68QURhp/FqvYaYUfBcpxC7UPP5oOWAmRTIVAzDV5iEYhY3b87bb755j6RiBlB1SudrduGRU83KqzkZbBPmvmEVGUQLhekP5SN/4ZhfvylJNEc4Xq0TfhVxxul5xD1+D0MNwD/FwZIhDoNONtB0hBZsGs2I8ML12VMAfvf63eon7y81L6XvH5twfvzMY9N/TrTh12NmtHxZqLnBcdeLjvmBnNyrVtiOMFBpBM+Ekcp6A++OsS2xVMpbtwyDMhPDJZSiguc6BS+OgbNmZyl66dWaUa15xEoyxWYxSYaliExHhu5O62EQtUxnsmYXHnJkWAHh6mce/8QbIx6vthfpAc+I6xv+3R6s8wqAl3qwnk7RyQaajtGCTaMZIV64PmsB+ORqzf+NOwuND3u+PPvMY9P/rpi39nSFRlIZNTecbrTCY00vOi4Vd1Rtv9GMHm750YVejD0LTJXsVyxL1Ic9jp2IIhVenaneuD6zujgZuPaU37AnAtcuB027HLacYtiyi6Hn5CPfzke+nYsCm8DrKl2RMHzDKhpK+ncmTs9+4dGPvzUG8Wrd8H9E3Gz+MFQB/D7i+nCDQCcbaDpGCzaNZgR54frsMwB+C0m3hG5gZjS9qFhrhsebreiYH8rJ3RrP+4Gc7kfLrGFQKlhX844xjKDyfZlbbs698vbKrZYvw06XEUrSZOBaU37DmQhcOx/5xdXcRPVO+XTDM52uS7uMAY8B+JsADtuT998BePvww9mXCMA/nLl8aVDiUDPiaMGm0YwoL1yfPQbgP8Ihm2JHUpnVRjDdaEXHW350TKrNHQ+q9eC9QaROHGYbw8Q2xXIxb86YptjRlThMml7kvnZ99Z37S81hF20dBz4M4Ndx+OoHP0GchRocekR7o5MNNF2hy3poNCPKp544t/rC9dk/QnyT+shB12MaIjo+mVs6PoklZobrRaW6G8axb6GcyOfMe0EjyJRgI4CThAEFQHBc1l8yJ90FTLEmiHzbEis5x1ge7mi3oxSrW7P1mddvrN1VSj8194hZ7NAa7ADUAJQArPZgXXtxqJqJmqOHFmwazQjzqSfOSQBffeH67LsA/gaA/GHWR0Qo5a1GKY6Jux1GylqpemcazfDJnRqkDwvDoIZtGUteIM9ahqhZllhtNMOnDYMatmksO5ZYMQ3ROkjXh36zVgtWX357+XqlHmhXWG+ZRRx7dtjjtIy40fxXAMwdcl27sQRAZ4ZqukK7RDWaMeGF67MTAP43AHpeiuP+kvvI9bu1/32v13tQCFClgvUzQ5BvmcKVim1mNtJQPNOkrpqgD4IwVMHPZio3rt+tZTKObgwgAP93AIcqZwPgXQAuYivbvzzsoHbhz2cuX/pRn9atGVPGqoqiRnOU+dQT52qIbzDfAXpbX+zMsfydnG1kwiJAhIgBUW+G7wulmkDc4VwQkTJNamZRrM0uNWe/+uP7f6nFWl9x0Jsiug8ACAGsAehHSZsIwGt9WK9mzNEuUY1mjPjUE+cUgL944frsDGJrW6kX6zUMIYt5644XyMMWJz0UQpCfs4zZVhA9BAZFkotRxHnTIG8YRXD3o9mK3P9/e/cW2+Z93nH8+/IkUtTBkmzHpg+hHTPJsjVd0rTBtizrOgQpyg7DLjqs2LCBF8VWDOsuigG8zO4IDOnNgm1dsbFDVhTD1g3owjaHBV3SJI0TH5LYbh3Tlhhbog+SZZ0okSL5vrv4U5acyBJFvjzY+n0AIgbJ988/nCB8+Pzf53lOZq9/eHlqaa7Te9kG9uFOEsIPDAIngVYE/2dUGSqNUIZN5C6UiEXGgH8Azru1ZmRnqOPzDm3b6fH5rPlQjy/XHw6c9nut+W4M1mzbqWYvzV148ejEMQVrbfNLLq7lBWaAVjQdVrGBNEQBm8hdKhGLFIDvAa9gqimbMtTfMxkO+rJNb6xBHssqWeCAhd/nmQfHEwr6rnRbsDY9W7r+6rHL776fnVYFaPtYmD5sbtlJa4I1FRtIw1R0ILINpLP5/ZiebTuaWae4XA2euzj75PRcqeE2Ilvl9ViFgN87aTtOwGNRcRw8vUHfRDeNlwJTVHBmbCZ7fnxustN72YYOAgmX13wVeMPlNVVsIA1Thk1kG0jEIuPAPwI/b2adYMBbfPjI8MsDYf8pd3a2uZ6A94plYfeF/GO2TU9v0Hepq4I1B2f8WmHixaPj7yhY65i6J0RswUHcvc9bxQbSFAVsIttEIhYpAv8BvID58mjYob39r/m8VlvuzSqWqgfKFXuwXKkO9AZ9H3m9Vqs70NetsFRZeOP9qyfePj2ZLS3bTf2dSlNGgAngPRfXPAQMuLjez1VsIM1QwCayjSRiEScRixwDvgM0PAFgaKDn+shg8F33dnYrC2zLwvZY1rLjOP5wyDfmsayyz2d1xReebTvVcxdnz7/49vjxK9NLXTlMfpu5D6gCr2OqOydcWNMHPObCOitUbCBNUVsPkW0oEYtcTWfz/wR8CfjVRtbYPRT88Or00u+4urGaYI93/J7h0NG9I70XJqYW78fB2rszlJ0rlEcKS+Xh4nJ1h9OhH5zXZ0tTx38xlZ1bLJc68fmyringYcyItguYxrf7XFj3PkxrnIUm15kELja/HdnOVHQgss2ls/mHgS/TQNPR8+Nzj45fK/yu+7uCgN9z9bEHd/5LqWz3VKu2f0d/z83ZjtWq7Z0rlHfML5VHloqV4YrtNNvdflPLZbt0ZvRG9sLEfNfNJhW8wO8BZzDB29cwjXTd8BNM5q4ZKjaQpilgExHS2fwIpop071avPXdx9rNXbyw9Wa06rjTpXWtHX+DErqHg6X27wmO3e4/jOCyWquH5wvJwoVgZLpnsm+XaJhycS9cKEyfPXR9bLttdN5tUbhrEZNa+jrsTCqaAv6fx6SEV4Fu5VLzrJnDInUUBm4gAkM7mfcBTwONbvbZUrgZOfHj9T0vL1Yibe4rs6v3RYDiQv2c4VPc9SRWTfRtaWCyPLJYqw1XbaTjTsrBYnj95bvrcVd2ndifwYgoF/sjldQuYez5nG7z+g1wq/l8u7ke2Kd3DJiIAJGKRCvDjdDY/hjleCtV7bY/fu9wX8uXcCthqxQa+gM+zuHsouKUbyH1eT3V4oGdqeKBnynEcFouVvrlCebhQrIyUytUB6si+VatONTs+N3pm9EbecbqrMa/c1j3Ak8Ac7lZ3hoE9NB6wqdhAXKGATURukYhFzqaz+cuYWaQH673O7/M0e2M2fp9nulyxh3sCniuOg9fCsi2r8dNNy7IIh/wL4ZB/AbhYrti+ucLy8MJSZXipWBmuOs4n7tubmilOHj97/fy8igruNJPAKUwhjdt2Ax82cN0UKjYQlyhgE5FPSMQis+ls/rvAbwNPUEdWaqlU3dXs51oW1cG+wHvhoO8yQH/Y7+oN/n6fpzIyGLw2Msg1x3EoLFX65xbLw4vFysh8oew/NTp9fiy/cN3Nz5S2KQNngc8DvS6v/cvAu2x9XNXxXCquDK24Qn3YRGRdiVjETsQirwLPU0dbg8P7+l+3LJq6Kb9qOyGvxypVbSfg9Vqlof5AyyoyLcuir9c/H9nZO3Zk/8BzY5fnvzqWX/gOJkvTFf3eZMsWaU1GK4zJsm2FJhuIq1R0ICKbSmfzfcDvs8mA7Tc/uPqNcsUecuMzY/sHnt+3OzzqxlobmAD+JxGLXFn7ZDSZ8WD6eB0BYoCrxRTSUvcDX3V5zQqQA763hWtUbCCu0pGoiGwqEYsspLP5fwN+A/gC62TnHcexqrbjylGU3+eZ7gl4C26sdRsl4H+B44lYxP74i7lU3AYu1R4/iSYzfawGb/cBLe/7Jg0bxRyP+l1c04f59x8B8nVeo2IDcZUCNhGpSyIWcYA30tn8R5iebYNrX7csyxkI+0/NzC83Pc4n4Pe0sjntGeDFRCxSd6uOXCq+gJlT+V4t+3YAE7zFMNWJ0j0qmEa3v4U5Il35ceFGn8D7qC9gU7GBuE5HoiKyZelsPoRp/fHg2ucXlsp9752b/lqlajfVVmEwHHh/R38geyjSf6aZdT5mBsgkYpGsi2sSTWYGuDX7tuWJEeI6P/AXwBVgDPiiS+vmgH+t430v5VLxn7n0mSKAMmwi0oBELLKUzub/Hfgs8DSmaSl9If9CdG/fD8+Pz/1xM+v7fJ4Fr8eqOI5jWZbV7K9KG3gLeC0Ri5SbXOsTcqn4HHACOBFNZryYVigr2bemK2elIWUgA3wauObiurswGbtPHKOvUUXFBtICyrCJSFPS2fwe4CvACJg5nz87fe2vKlWnv9E1/T7Pjc88uPPbS8VK39BATzNtNi4BLyRikatNrNGwaDKzg9Xg7RDu3lclGwthgqeHMNlgNxSBZzHHrrdzKpeK/8ClzxO5SRk2EWlKIha5ks7mvw3EgU97vZ5qX8ifnVlYfrTRNSsVe4fXY9kej7VRJmMjRVaLCjr2qzSXis9g+ne9G01mfMC9rAZwI53a1zbhYGbj7nd53c3+m1SxgbSEAjYRaVoiFlkG/judzY8C8YG+QK6ZgM0Ba3KmGAn4PFttVAqmj9pLiVik6ckLbsql4hXgQu3xYjSZGWY1eIui/x+7rYgpEHjAxTUdNg7YpoCPXPw8kZt0JCoirkpn8ztt2/mDY2envrlYrBxudJ09I6GXe3t80/t2h7Pe+jJt05iigguNfmanRJMZP+bIdCWA29HRDd09dtUeB4HP0HxQPA383Qavq9hAWkYBm4i4Lp3N+147eTkVDvqeLhQrDzrO1r8ovV5r4YGDg9/3eqzKyGBwoxvHq8CbwE9bUVTQbtFkxsIcl64Eb/dSK+qQhkSAy8CnMM2fm5Hj9lWiVeDZXCq+2ORniKxLAZuItEw6m+8vV+zfHJ2Y//PL1xe/TB0zSVdYFvbjD+3628nZ0v4Du8Pnb/O2jzBFBZOubLgLRZOZHm7NvjXVMmUb2wP8WZNrnAH+8zavqdhAWkr3TIhIy9Sa0/7o2cWxn84sLP9lOOT7/OzC8qPlir3pDfeOg+fK9NLh4YGeiXVeXgJeBt7rZFFBO+RS8RJmqPnZWvZtN6t93w6imdD1uoa5B63uHw3rmNngNRUbSEspwyYibZPO5ntn5ktfOD8+/0yhWH7EcTYONnYOBt84tK//7XDQt3ZM1fvAy4lYpJWjq+4I0WQmCBxmNfvmRjf/u9k3ae7vKAMcW+f568BzuVRcX6jSMgrYRKTtnnn9XPjcpdlnBvsCv+bxWMMei0ql6vQtlSoH1t7v5vN6ZveMhF47sn/gJOZL8YVELDLWuZ13r1r2bQ+rwdt+mssm3Y2+gunL1qjvA+fWef7lXCr+VhPrimxKR6Ii0nbPPHl/AfjrdDZvYY71fgV4yLad/lOjN56+MVf6HEDVtsOYJqX/B7yRiEU2ali6rdWyO5drj9ejyUwIMyprJYDr7eD2usVRTG+2oQavL63zXBUzZ1akpZRhE5GukM7mPcChpVLlkbO52W8UipVPPXBw4G92DYWeT8QizUw72PZq2bcIq8FbhO2ZfQtiMmyHMT8UtjKNYwn4Z0ymd63TuVT8doUIIq5RwCYiXSedzfteeWfisac+t+/o3V5U0AnRZCbMrQPrQ53dUVt5MG1S7gUexrT7qMcU8APMQPm1vptLxXOu7U7kNhSwiYhsY9FkxgPsYzX7trezO2qrRzAj1erpc/dDTMHL2ibOKjaQtlHAJiIiN0WTmX5uzb71dHZHLfeHbD6+agxTcPDxxswqNpC2UdGBiIjclEvF54GTwMloMuMFDmCCtyPAPZ3cW4u8w+YB22k+GaxVMRk3kbZQhk1EROoSTWYGWc2+HQYCnd2RK4LA19l4gsRzwAK3Vomq2EDaShk2ERGpSy4Vn8V09D8eTWZ8mErLlXvfdnZyb00oAh8AT2zwnl2YKtG1NNlA2koBm4iIbFkuFa8Ao7XHS9FkZojVo9NDgL+D29uqEmZG6BOY5sNrXcGMtFo71H0aMwhepG10JCoiIq6KJjN+TNuMlezbcGd3tKk/wQRkC8DjtecWgQnMKKobwOSa97+SS8XfbOsOZdtTwCYiIi0VTWZGWA3eotTXRqOd9gFPAXngLDACXMAEcA9hig5WVIFv5VLxbT/LVtpLAZuIiLRNNJkJYI5MjwD3A4Od3dG6fJhJEAFMgFZc85qKDaQjFLCJiEhH1EZm7WQ1+3YvZhJBJ3gw96o5mGAtXPvnfO31AnAKOJpLxW90ZIeyrSlgExGRrhBNZnow7UJWAritzPp0UxAz0aCAOSL9ALiQS8XtDa8SaSEFbCIi0nVq2bd7WO37dgB3s299mMICL6tNcb2Y49BfYJrins2l4qX1LxdpLwVsIiLS9aLJTIjV7NsRTMDVCC9m2P1I7c+jmIzaEvAW8HZt2oNIV1HAJiIid5Ra9m0vq0ene6ivr6gF9AK/jhm/dRF4Dfgwl4rPtGSzIi5RwCYiIne0WgDXh6k43VF7DAJDwH5MoDZXe24eyABjuVRcX4Byx1DAJiIiItLlOlU+LSIiIiJ1UsAmIiIi0uUUsImIiIh0OQVsIiIiIl1OAZuIiIhIl1PAJiIiItLlFLCJiIiIdDkFbCIiIiJdTgGbiIiISJdTwCYiIiLS5RSwiYiIiHQ5BWwiIiIiXU4Bm4iIiEiXU8AmIiIi0uUUsImIiIh0OQVsIiIiIl1OAZuIiIhIl1PAJiIiItLlFLCJiIiIdDkFbCIiIiJdTgGbiIiISJf7f5vrbak5B/btAAAAAElFTkSuQmCC\n",
"text/plain": [
"