{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
""
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import numpy as np\n",
"import os\n",
"import pandas as pd\n",
"\n",
"from bestiapop import bestiapop\n",
"\n",
"#Extend screen\n",
"from IPython.core.display import display, HTML\n",
"display(HTML(\"\"))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Working with lists of lat/lon combinations\n",
"\n",
"In this example we will download climate data based off a list of lat/lon datapoints\n",
"\n",
"### Define your climate extraction variables\n",
"\n",
"We need to define the scope of what type of data and output we expect from BestiaPop.\n",
"\n",
"> **NOTE**: the `coordinate_list` variable contains lists of lat/lon pairs where the first value is `latitutde` and the second value is `longitude`"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example: getting climate data for the entire potato production areas of Tasmania, Australia \n",
"\n",
"In this example we retrieved 1724 grids of 5 km resolution from SILO. This lat/lon combinations were masked for the entire potato production areas in Tasmania, Australia. Then we put it in a dataframe and work with the data to create statistical summary and maps of each climate variable."
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"# Define variables to pass in to Bestiapop\n",
"\n",
"#You need to host your lat/lon list in some web server like GitHub and then provide the url. This example is using a csv hosted in our GitHub repo (sample-data)\n",
"from csv import reader\n",
"url = \"https://raw.githubusercontent.com/JJguri/bestiapop/master/sample-data/bpop_tassie_example.csv\" #specify url\n",
"df = pd.read_csv(url)\n",
"coordinate_list = df.values.tolist()\n",
"action='generate-climate-file' #specify action\n",
"data_source='nasapower' #specify dta source\n",
"year_range='2019-2020' #specify range of years\n",
"climate_variables=['radiation', 'max_temp', 'min_temp', 'daily_rain'] #specify climate variables\n",
"output_type='dataframe'\n",
"output_path=None"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"#Get MET filed directly in your computer\n",
"\n",
"# #You need to host your lat/lon list in some web server like GitHub and then provide the url. This example is using a csv hosted in our GitHub repo (sample-data)\n",
"# from csv import reader\n",
"# url = \"https://raw.githubusercontent.com/JJguri/bestiapop/master/sample-data/purdue.csv\" #specify url\n",
"# df = pd.read_csv(url)\n",
"# coordinate_list = df.values.tolist()\n",
"# action='generate-climate-file' #specify action\n",
"# data_source='nasapower' #specify dta source\n",
"# year_range='2019-2020' #specify range of years\n",
"# climate_variables=['radiation', 'max_temp', 'min_temp', 'daily_rain'] #specify climate variables\n",
"# output_type='met'\n",
"# output_path=r'C:\\Users\\uqjojeda\\Dropbox\\github\\ClimateDataAutomation\\sample-data'"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Instantiate and run BestiaPop\n",
"\n",
"In this instance we will iterate over the lat/lon combinations"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST\u001b[0m - Extracting data and converting to met format\n",
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mGenerating DataFrames\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Fetching data for 40.297382/-86.901067\n",
"Total Progress: 0%| | 0/1 [00:00, ?it/s]"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable radiation\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 40.297382 - Lon -86.901067 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:31:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 40.297382 - Lon -86.901067 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:31:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:56\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:56\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 40.297382 - Lon -86.901067 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:56\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:31:56\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 40.297382 - Lon -86.901067 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:31:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:31:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:03\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:03\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 40.297382 - Lon -86.901067 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:03\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:03\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:07\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable max_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:07\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:07\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 40.297382 - Lon -86.901067 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:07\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:07\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:10\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:10\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 40.297382 - Lon -86.901067 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:10\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:11\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:14\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:14\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 40.297382 - Lon -86.901067 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:14\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:14\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 40.297382 - Lon -86.901067 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 40.297382 - Lon -86.901067 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable min_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 40.297382 - Lon -86.901067 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 40.297382 - Lon -86.901067 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 40.297382 - Lon -86.901067 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:39\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:39\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 40.297382 - Lon -86.901067 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:39\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:39\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:42\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:42\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 40.297382 - Lon -86.901067 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:42\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:42\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:46\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable daily_rain\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:46\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:46\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 40.297382 - Lon -86.901067 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:46\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:46\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 40.297382 - Lon -86.901067 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:49\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 40.297382 - Lon -86.901067 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 40.297382 - Lon -86.901067 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:32:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:32:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 40.297382 - Lon -86.901067 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Total Progress: 100%|####################################################################| 1/1 [01:16<00:00, 76.20s/it]"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST\u001b[0m - Extracting data and converting to met format\n",
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mGenerating DataFrames\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Fetching data for 39.032738/-85.528156\n",
"Total Progress: 0%| | 0/1 [00:00, ?it/s]"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable radiation\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 39.032738 - Lon -85.528156 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:06\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 39.032738 - Lon -85.528156 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:13\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:13\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 39.032738 - Lon -85.528156 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:13\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:13\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:18\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:18\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 39.032738 - Lon -85.528156 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:18\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:18\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:22\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:22\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 39.032738 - Lon -85.528156 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:22\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:22\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable max_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 39.032738 - Lon -85.528156 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:29\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:29\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 39.032738 - Lon -85.528156 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:29\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:29\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 39.032738 - Lon -85.528156 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 39.032738 - Lon -85.528156 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 39.032738 - Lon -85.528156 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:44\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable min_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:44\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:44\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 39.032738 - Lon -85.528156 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:44\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:44\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 39.032738 - Lon -85.528156 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:51\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:51\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 39.032738 - Lon -85.528156 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:51\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:51\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 39.032738 - Lon -85.528156 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:58\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:58\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 39.032738 - Lon -85.528156 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:33:58\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:33:58\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:01\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable daily_rain\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:01\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:01\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 39.032738 - Lon -85.528156 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:01\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:01\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 39.032738 - Lon -85.528156 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:08\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:08\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 39.032738 - Lon -85.528156 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:08\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:08\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 39.032738 - Lon -85.528156 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 39.032738 - Lon -85.528156 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Total Progress: 100%|####################################################################| 1/1 [01:13<00:00, 73.99s/it]"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST\u001b[0m - Extracting data and converting to met format\n",
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mGenerating DataFrames\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Fetching data for 41.105156/-85.397344\n",
"Total Progress: 0%| | 0/1 [00:00, ?it/s]"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable radiation\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 41.105156 - Lon -85.397344 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:20\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 41.105156 - Lon -85.397344 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 41.105156 - Lon -85.397344 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:27\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 41.105156 - Lon -85.397344 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:30\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable radiation - Lat 41.105156 - Lon -85.397344 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:34\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable max_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 41.105156 - Lon -85.397344 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:37\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 41.105156 - Lon -85.397344 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:41\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:45\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:45\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 41.105156 - Lon -85.397344 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:45\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:45\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 41.105156 - Lon -85.397344 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:48\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable max_temp - Lat 41.105156 - Lon -85.397344 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:52\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable min_temp\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 41.105156 - Lon -85.397344 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:55\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 41.105156 - Lon -85.397344 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:34:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:34:59\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 41.105156 - Lon -85.397344 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:02\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 41.105156 - Lon -85.397344 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:05\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable min_temp - Lat 41.105156 - Lon -85.397344 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:09\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for climate variable daily_rain\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 41.105156 - Lon -85.397344 for Year 2011\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:12\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 41.105156 - Lon -85.397344 for Year 2012\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:16\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 41.105156 - Lon -85.397344 for Year 2013\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:19\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 41.105156 - Lon -85.397344 for Year 2014\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:23\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing data for year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mProcessing Variable daily_rain - Lat 41.105156 - Lon -85.397344 for Year 2015\u001b[0m\n",
"\u001b[32m2021-04-12 15:35:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - Extracting data from NASA POWER Climate DataBase\n",
"\u001b[32m2021-04-12 15:35:26\u001b[0m - \u001b[34mPOPBEAST.NASAPOWER_CONNECTOR\u001b[0m - \u001b[32mNeed to get data from the NASA Power Cloud\u001b[0m\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Total Progress: 100%|####################################################################| 1/1 [01:10<00:00, 70.12s/it]\n"
]
},
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
\n",
" \n",
" \n",
" \n",
"
\n",
"
"
],
"text/plain": [
"Empty DataFrame\n",
"Columns: []\n",
"Index: []"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"final_df = pd.DataFrame()\n",
"\n",
"for coord in coordinate_list:\n",
" print('Fetching data for {}/{}'.format(coord[0], coord[1]))\n",
" climate_data = bestiapop.CLIMATEBEAST(\n",
" action='generate_met',\n",
" data_source=data_source,\n",
" output_path=output_path,\n",
" output_type=output_type,\n",
" input_path=None,\n",
" climate_variables=climate_variables,\n",
" year_range=year_range,\n",
" lat_range=[coord[0]],\n",
" lon_range=[coord[1]],\n",
" multiprocessing=None\n",
" )\n",
" data = climate_data.process_records(action)\n",
" final_df = final_df.append(data, ignore_index=True)\n",
"\n",
"# Show the final dataframe\n",
"final_df"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Calculating yearly statistics\n",
"\n",
"Let's calculate some statistics base on the climate data we retrieved from SILO\n",
"\n",
"> NOTE: here we presente the sum of rain and radiation and the mean of daily minimum and maximum temperature across years. You can change the parameters you want to get lookin at the [GroupBy](https://pandas.pydata.org/pandas-docs/stable/reference/groupby.html) documentation."
]
},
{
"cell_type": "code",
"execution_count": 96,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
\n",
"
\n",
"
radn
\n",
"
maxt
\n",
"
mint
\n",
"
rain
\n",
"
\n",
"
\n",
"
lon
\n",
"
lat
\n",
"
year
\n",
"
\n",
"
\n",
"
\n",
"
\n",
"
\n",
" \n",
" \n",
"
\n",
"
143.85
\n",
"
-40.05
\n",
"
1991
\n",
"
4876.8
\n",
"
16.6
\n",
"
10.2
\n",
"
1020.9
\n",
"
\n",
"
\n",
"
1992
\n",
"
4663.2
\n",
"
16.2
\n",
"
10.0
\n",
"
1399.2
\n",
"
\n",
"
\n",
"
1993
\n",
"
4901.0
\n",
"
16.9
\n",
"
10.5
\n",
"
922.1
\n",
"
\n",
"
\n",
"
1994
\n",
"
5016.2
\n",
"
16.6
\n",
"
10.0
\n",
"
870.2
\n",
"
\n",
"
\n",
"
1995
\n",
"
4913.9
\n",
"
16.3
\n",
"
9.5
\n",
"
765.7
\n",
"
\n",
"
\n",
"
...
\n",
"
...
\n",
"
...
\n",
"
...
\n",
"
...
\n",
"
...
\n",
"
...
\n",
"
\n",
"
\n",
"
148.45
\n",
"
-40.40
\n",
"
2016
\n",
"
5105.7
\n",
"
18.1
\n",
"
11.9
\n",
"
961.0
\n",
"
\n",
"
\n",
"
2017
\n",
"
5356.2
\n",
"
17.9
\n",
"
11.5
\n",
"
633.9
\n",
"
\n",
"
\n",
"
2018
\n",
"
4942.9
\n",
"
18.0
\n",
"
11.6
\n",
"
709.5
\n",
"
\n",
"
\n",
"
2019
\n",
"
5088.6
\n",
"
17.9
\n",
"
11.2
\n",
"
585.0
\n",
"
\n",
"
\n",
"
2020
\n",
"
4664.0
\n",
"
17.5
\n",
"
11.1
\n",
"
741.3
\n",
"
\n",
" \n",
"
\n",
"
51660 rows × 4 columns
\n",
"
"
],
"text/plain": [
" radn maxt mint rain\n",
"lon lat year \n",
"143.85 -40.05 1991 4876.8 16.6 10.2 1020.9\n",
" 1992 4663.2 16.2 10.0 1399.2\n",
" 1993 4901.0 16.9 10.5 922.1\n",
" 1994 5016.2 16.6 10.0 870.2\n",
" 1995 4913.9 16.3 9.5 765.7\n",
"... ... ... ... ...\n",
"148.45 -40.40 2016 5105.7 18.1 11.9 961.0\n",
" 2017 5356.2 17.9 11.5 633.9\n",
" 2018 4942.9 18.0 11.6 709.5\n",
" 2019 5088.6 17.9 11.2 585.0\n",
" 2020 4664.0 17.5 11.1 741.3\n",
"\n",
"[51660 rows x 4 columns]"
]
},
"execution_count": 96,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"#calculate yearly mean of mint and maxt (degree C) and cummulative radn (MJ m-2) and rain (mm) during across years and lat/lon combinations\n",
"mini = final_df.groupby(['lon','lat','year'], as_index=True)['maxt'].mean().round(decimals=1)\n",
"maxi = final_df.groupby(['lon','lat','year'], as_index=True)['mint'].mean().round(decimals=1)\n",
"radn = final_df.groupby(['lon','lat','year'], as_index=True)['radn'].sum().round(decimals=1)\n",
"rain = final_df.groupby(['lon','lat','year'], as_index=True)['rain'].sum().round(decimals=1)\n",
"st0 = pd.merge(radn, mini, on=(['lon','lat','year']), how='outer')\n",
"st00 = pd.merge(maxi, rain, on=(['lon','lat','year']), how='outer')\n",
"st_final = pd.merge(st0, st00, on=(['lon','lat','year']), how='outer')\n",
"st_final\n",
"\n",
"#Save the statistical summary in the current dir\n",
"#st_final.to_csv(os.getcwd())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## This code re-structure the datasets before mapping"
]
},
{
"cell_type": "code",
"execution_count": 84,
"metadata": {},
"outputs": [],
"source": [
"#Preparinf dataset to plot\n",
"radnrain = final_df.groupby(['lat', 'lon','year'],as_index=True).sum()\n",
"temp = final_df.groupby(['lat', 'lon','year'],as_index=True).mean()\n",
"#Reseting index\n",
"Pre_nc = radnrain.reset_index()\n",
"#Set indexes\n",
"Pre_nc.set_index(['lat', 'lon', 'year'], inplace=True)\n",
"df1 = Pre_nc[~Pre_nc.index.duplicated(keep='first')]\n",
"#From 2D dataframe to xarray\n",
"nc_r = df1.to_xarray()\n",
"#Reseting index\n",
"Pre_nc2 = temp.reset_index()\n",
"#Set indexes\n",
"Pre_nc2.set_index(['lat', 'lon', 'year'], inplace=True)\n",
"df2 = Pre_nc2[~Pre_nc2.index.duplicated(keep='first')]\n",
"#From 2D dataframe to xarray\n",
"nc_t = df2.to_xarray()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Mapping rainfall\n",
"\n",
"> NOTE: Each 5 km grid represents the mean (across years) of cumulative rainfall during the year"
]
},
{
"cell_type": "code",
"execution_count": 93,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAi4AAAHYCAYAAABqaWNGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAACyw0lEQVR4nOzdd3zbx3n48c9hcu8pLk1rb1mWLFny1LDkEY84w0mcNGnSJnFG21/2bJs2bZrRxGmaOKmdxFlO4j1kyZZly5KsPa0tapMU9ybm/f4ASOIgboIkQD7v1wsv4vC97xdHCQQe3D13p7TWCCGEEELEAstIN0AIIYQQoq8kcBFCCCFEzJDARQghhBAxQwIXIYQQQsQMCVyEEEIIETMkcBFCCCFEzLCNdAOGwpo1a/TLL7880s0QQggR3dRIN0D036jscamqqhrpJgghhBBiCIzKwEUIIYQQo5MELkIIIYSIGRK4CCGEECJmSOAihBBCiJghgUuE+P3+kW6CEEIIMepJ4DIIbrcbj8fDkSNH2LlzJx6PZ6SbJIQQQoxqo3Idl+HgcrnYsmULTqeTb/3gOErBJz94hTvWr8Vut49084QQQohRSXpcBkhrjdfrZeXKlTS32WlqtfPT35Tx0ksbcbvdI908IYQQYlSSwGWAbDYbTqfTyG1pbHHw349d4Ge/+APnzp3D7XZL7osQQggRQRK4DJDNZmPSpElcvnzZeLyxxcETz3nxer289NJLXL58WXpghBBCiAiRwGUQCgoKKC0t7eKIYsKRf+cO9TxVr36PQwcPSPAihBBCRIAk5w6C1prW1la2fXCz8bhKcnbcn5PezMYDz2OZN3+4myeEEEKMOtLjMgh+vx+tda/1ZiRUcubMGXw+3zC0SgghhBi9JHAZJKfT2WudgkQPp06dksBFCCGEGCQJXAbI4/GwceNGli5d2qf6K1asYPfu3TLLSAghhBgECVz6ye/343a7eenpP7G49hlsT/9tn85L2P11GmvKZXVdIYQQYhAkObeP/H4/Ho+Hmpoa3n7p99zoOEyyMzD0oxLCVsqNCxs+cjjQGuzxKdhs8k8uhBBCDJR8ivaBx+Nhz549VFdXU1JSwvqEAyjVv2s0uRUZuVlYrdahaaQQQggxBkjg0gdHjhxh8uTJLFmyBADfkf5fo95lITMzM8ItE0KIgfF6vdIDLGKS5Lj0QXp6Om1tbYO6RoNLkZKSEqEWCSHEwHk8HioqKiTnTsQkCbf7IDc3l8OHD1NYWAiA9YHHjeP+1z5rlFVionkBiwV7s8Lr9Q5lM4UQoldut5vXX3mZlu2buPGfvk1KSgoWi3yH7U1K2gzt9TRH5FqtLec3aK3XRORiY5AELn3Q1taGw+EY1DUSHZqmpiYZLhJCDBuPx4PFYunIrXO73ezds5vCFx8lz93Aq798hDs+/U+Dfn8bC7yeZqbN+kJErrVv5yezInKhMUrC7F54PB42b97MjBkzBnWdjAQoKyuTdVyEEMNGa8358+dxu9243W5OnjyJ9fc/ZJy7AQtw/fHX2bBhg+yl1gdKgbKoiNzE4EiPSy+UUvh8vkF3pcY7oPr8CdzuecTFxUWodUII0T2LxcKzf/0LkydMwNtUz4Q9LzKlparjeLzfw4IFCzhw4ABz586VnpdeKKsEHdFAApde2Gw2Vq5cydGjR5k5cyYA+u2vDOhaCe5zksUvhBhWa9bfgeunX2a8uxYFhG88Yv3MvTRlj6fqa98nNycHq7xHiSgnr9A+yMrK4p133un2uMrMMB/o5lvLhAI/ZWVlFBUVRbJ5QgjRJZvNxpQpU3hp1UeYs27dVcfL7wss8TCt8izbf/YD1n7p28RL4NI1BUgSc1SQ/4U+UEpFJDelMENx5swZyXMRQgwLt9vNoUOHGDduXI/1FHDt0bd44de/kinS3VIoa2RufXo2pZKUUt9USj2nlCpTSmml1GM91B+nlPq5UuqiUsoV/PkXpVRKWD2bUuprSqlSpVSbUuqYUupTSl29rGp/6g4nCVz6QGuN1nrQ17FaFW1NtfLGIIQYcl6vl2PHjmG325k/f36v9e3az+wtf+HNza9Jsm50yAK+ASwEdvdUUSk1DdgH3Az8L/B3wE8AJ5AQVv1/gG8DG4FPAQeBHwNf6+LS/ak7bKRPsA8iEbQMxbWEEKI7NpuNSZMmcfTo0T6fk+puofb3j3BuwkTGjx+P3W7v/aSxIjiraBiVAYVa60tKKRvQ5TfeYO/Hb4GLwEqtdVN3F1RKzQU+CvxAa/354MOPKqWeBL6slPqF1rqsv3WHm/S49IHWGkvZZvw7voR/x5cGfB23V5OSkSuZ+0KIIae1prm5ud/vN+Mbr3DiZ/9Ja2srWoa1DcM5VKS1dmmtL/Wh6s0EemW+obVuUkrFK6W6izgfCP78UdjjPyLQO3P3AOsOKwlc+qCxsZGsxJA/4KwM85aSbN6Swm4pqZCSyoWWeCZMmMAIDw8KIUYxl8uF1+vl7Nmz7Nq1i9mzZ3db12qzGDeLRWGxKBZeOMyGnz+CR1b7jgWrgz+blVI7gBagTSn1mlJqZljdRUCF1vpc2OM7AT+BAGggdYeVBC698Pl8XLx4kbzkwQ/xlJYrsrOzI9AqIYQwud1uamtrefXVV9myZQt2u53bb799QF+ULGiu3fMyrzzzlOS7tGufVRSJW2RdE/z5JwLDRfcDnwfmAG8opQpC6o4DrurF0Vq7gWpgoHWHleS49KK+vp6KigpmJA4ucNFa47EW4vP5ZC0XIUREeTweSktLqa6uZs2aNRHZeyjB56b4+cc4NHEyc+bMGfP5Loq+D/P0QZZSKjTh9uda658P8FpJwZ/7tNb3tT8YvP5W4B8IBDIA8UBDN9dpCx5nAHWHlfS49MDlcrFjxw5uvvnmQV+rut7H+InTcTqdEWiZEEIE+P1+GhsbuXLlCtdff31EN0zMa63D5/NRU1ODzxe+dJ0YhCqt9aKQ20CDFoDW4M8nQh/UWr8FnAVWhtXt7kMoLuRa/a07rCRw6YHWGqvVGpGclEtXPJSUlESgVUIIEeD1emltbWXr1q0sX758SJ5j/Pc/y7b//aEs4xC9exVdDv4s7+JYBZAeVveqRX2UUg4gM+Ra/a07rCRw6YHVau349qKSkztuJKaat+TMsFuGcVPJWRCXLEm5QoiI8vv9VFZWkpiY2OWxS5cuceLECZqaup4he9WMl64+aIHrjm7lxSceH/PBy3DOKuqHXcGfhV0cKwQqQ8p7gDylVHFYvWsJxAN7Blh3WEng0gObzRaxfJTcLCculysi1xJCCACHw0FhYSELFizgueee63iPqays5Nlnn8Xj8ZCRkcFrr702qOexax+ztvyZbW++Icm60ecZAsM2f6OUsrY/qJS6nUAC7YaQun8K/nw47BoPA27g6QHWHVaSJdoDj8dDQkL4ooMDk5Xh5KU3NrH29rsBxnyimxAiMmw2G6mpqaxevZqNGzdisVjIyMhg3bp12O12tNYUFhbS0NBASkpK7xfsRrq7mco//g81M2aSl5cXwd8gRgz/AnQopT4FpNHZyTBHKfXV4P1ntdYHtdaVSqmvAd8DXgsuEDcO+AxQCvyg/Xpa631KqV8Bn1dKJROY2rwKeDfwLa315YHUHW4jHri0R4ha66jL/HK5XBELXOx2C8vmudj08uNk5s5h7ty5kqgrhIgIi8WC0+nktttuw2634/V6O74cKaWYMmUKZ86cYe7cuYN6nin1l3l71y5Wr14tC2kOj38EQpMj5wdvEJj6fBBAa/1fSqlq4HMEAphG4EngS1rr2rBrfgI4D3wYeIhAAu9nCCzlH64/dYfNgAIXpdQNwBvBYpHW+mLIsRTgO8B9QCpwAPia1npjF9f5LoFuJ59S6oda66+GHc8AvgCsB8YDLgLzyt8Cfqi1PjaQ9veV3+/HarX2XrGPMtOdrF3pZPO2PVitCyJ2XSGEADq+DIUHFQkJCVz43c8o+I+DxGcOfBarAvJffoyKefMoKCiI6AymqKcUWIf399Vaj+9H3ceAx/pQzwN8K3iLWN3h1O/AJbhnwiNAM5AYdkwBzwGLge/TGam9qJS6VWu9JaTue4D3EYjoAL6jlDqgtX4yeDyNwMZS2cDjwKHg880gsMDODmBIA5f4+HguXrzIjBkzUDO/3Hmg4Umzoi1s2Ccu2Sxrc9nsxOQGWc9FCDFsPB4P/mnz4eJB43FbvPne5ff23PGtLBYmtlax+cnfcuen/xEVfGysGO6hItG1gbziPg3kAY92ceweYAXwd1rrr2it/xe4ETgH/FdY3eXA97TWj2utHwf+M3huu48CE4B7tdaf0lr/r9b6+1rrjwL5wPMDaHu/WCwWJk2axKuvvkpra+SmrNc3EdGeHCGE6ElcXBxTltxAtSOp98q9UMCio6/z3M9/zOkzZ/B4PLLGixhW/QpclFL5wDeBLwN1XVR5d/Dx37Y/oLVuA34JLFRKTQqpexr4gFJqvlJqPvAgcCLk+OTgzzcIo7V2a62r+tP2gbDZbBQXF7N8+XI2btxITU1NRK7r99nlD10IMawmT57MqazxEblWiq+Nm/Y8g+u/Ps+GDRs4fPgw9fX1+Hw+XC4XbW1teEfZPkdKRe106DGnv2MV3wNOAr8Cvt7F8QUElh0Of8XuDDl+Onj/ZwR6aPYGyxuB0NUDzwZ/fgT4aT/bGTHtCW5r167lxRdf5M4772SwLztNINNfCCGGi8/nQ89aDHtO9165j8a56pm5fj01NTVcuHCB6upq2tra8Pl8OBwOxo0bx7Rp00ZNLowMFUWHPgcuSqmVwHuB5VprfzeLqeXTGaSEKgv+7FiFT2vdGrzmTEADR7T5af4o8FngEaXUPwBbCOy78KLWuqsVAoeU3W5n2rRpVFVVkT3oyUCSjS+EGF42m424vMjvi9f8udtwEhjXn/WDzjkYWmteeOEFJk6cSFxcXMSfV4xdfQpcQhJyn9Bab+uhajyBmT/h2kKOd9Ba+wkk3V5Fa12llFoEfIlAz8yHgzetlHoceFhr3diX9kdKQUEBp0+fJquwzXhcOa9etTKUz9J53O32YrElSY6LEGJYeb1e3D4/if+9qeOx5odvDatlM463/L81xlF/i7lyrrJ335OilCIuLm70TEJQMswTLfraf/cZAnPJv9BLve42ZYoLOd5nWuuLWutPaq3zg8//QQLLDD/ECMwjT0xMHFSei8vlZcNL5axYsVIWoBNCDCur1Up8/PBu6KuUGjXDRABYLJG5iUHp9V9QKZUKfINAXotDKTVeKTWewGp+AIVKqfY9EsroYlMmAkNIMIhNmbTW57XWvwGWAaeA9wZ7gobVYHJTTp6o5KabVku3qRBiWPl8Pg4ePMi8efNGuilCDFpfQr90IJnAQnGlIbfPBI9vJ5B7AoFE23ldBBTXBX/uG1RrCcwoIrConQPIGuz1+kMpNaiNEqurW0hOTu69ohBCRJDf76esrIz09PTeK4tuRenu0GNOX3osrgDv6uLx9wAPEFhv5ULwsT8TmBL9IMEV/JRScQRmBu3TWp/qa8OUUksIJOw2hj2eDiwFajB3vRwWAw1c3G4vtTUKr9crS2ULIYaV3W6nubmZDRs2kPvc/zK+LTDkbU+4+iOg9ZvrB/w8TZ+6xSire74y4GtFneB0aDHyeg1ctNYtdLELpFJqXvDuhpAl//9CoPflf5RSkwkENA8RWK7/tn627UHgg0qpZwhs290EFAMfIjAc9cnh3t+oY7Vbf9jThq2Mi7szlcfv1+w6cJn6hkzWrbtv9CSqCSFihsfj4a677sJms3Fp2jQ2//5XLDu9FWdy2JeosGRblWAj/pvdr/XZ8g+rOu77XbI2lRgeEf0UDU6TXg/8G/AxIIXArKF1WuvN/bzczwgEK7cAqwkMWdURSM59WGv9XKTa3Vfl5eUUFBQQSLHpndaa51+pYsWN7yIhIUF6WoQQI8Jut3dMCCgsLCT3s19iy7c+w4qOznLRG4UM80SLAQcuWutvElhFN/zxeuDvg7cB01ofBr44mGtEUn19Pdu3b+fuu++G8L02u7F7XxXLbriT1NTUQeXGCCFEpFitVrxeL75Ft6D3PjboBTW7oxncZIaooxj2TRZF1+R/oY/q6upYtGhRn4Z6PB4/Z842UNecT3JysgQtQoio4nQ6yZ84hZbIdrobqqwJY28HaTEsJOGijzIzMzm14xG0dkBS15n51TVtbNvVRFrGJMaPn89NN+VLTosQIiolJCRQoRIp0E04iWx+ig/Fvoyp3DFpUu+VY4gMFUUH+VTto/j4eKobgt8caq4YxzyVFWyvWorDkcWa2xditVrlW4YQIqoVFBRw7r1f5WhDA2V7djBny1Mkel3kP/V23y9i71wBXPkDw0JNysHWSTex6oMfj3STR5bMKooaErj0kcfjoU3lAtXG46fLfbxTXcItqxficDhkRVwhRExISkpi5syZ+P1+vAsX8gyK61/9/aCuecGaQuny93Pnujs63gv9fr98kRMRJa+mPnI4HIy/ZiEXKwMbXze1al46EAcF97Nm/X0kJiZK0CKEiDkWiwWHw8HKd7+Ps6l5A75Og3JwefXHueX29R3vhT6fz/gZ2yKz+JwMNw2eBC59ZLFYmDJlCrvP5fH6IS+7quZz692fYMKECTLNWQgR89LS0qhctHLA5+8oWsoNN97U8X7o9/vxer387ne/w+/393J2DAgOFUXiJgZHhor6weFwcNe970drjc/nk4BFCDFqWK1W0mctoP6Lt2MnkK/imGXuqqKbPV2dihsLOYtv6phB6XK5cLvdvPbaa8ydO5fq6mry8gbemyNEKAlc+ql9rNZqtfZSUwghYofVaqW4pIQaFU+ubunTOZbgyrtxgOOFn/FS6XHyp82h/PRxEo6+xfqv/RilFBs3biQjIyOmv+wpQEmuTlSQwEUIIQQQWK+qQHfdq9Kbmd5KZhz8E20H/8p0vIH12kK+4MV8gq4Ci+SnRIUYfyUJIYSIlPr6ehIZWOACgV6J+GDQAlBRUcGlS5coKCigqakpIm0UQgIXIYQQAMyYMYN9joKIXMsP7N69m8rKSpqamjh+/HiMJ+kqLJbI3MTgSOAihBACgPT0dFLue5jj1sxBX6tKJTB16lTmzJlDcnIyly9fxuMZeG/OSFMyqyhqSI6LEEIIILCL9LRp0zhy3z9wPj6eadOmGcd9Tz5klHVLZyCi27zGsSqVxDU5OR1LSZSUlMR+nouICvIqEkII0cHhcDBz5kzq6+uprKzstX6z38arlum8PuN+XstbwR6VSxN2CnUjb7zxBh6PB4fDQVJSUswv0ilDRdFBelyEEEIYlFLYbLY+9ZC8kbKUtR/6+46l/f1+P7t37qT5rWdYunQpL7zwAtdffz1paWkxPR1aRA8JXIQQQhhaWlpoaWkhM7PnXJcGv52Sa2+8Ksi5bulSPIsWsWHDBm6//XZefvll1q1bN9TNHloyHTpqyFCREEKIDn6/n/Lycsaf/DW+332wx7pVfifjxo27qmfGbrcTHx/H2rVrefXVV8nLy4vxGUWBqd4yVBQdpMdFCCFEB4vFwqRJk3h2+3xWeHYTuui/9f7HjLqX/+tTTExL6/I6SlmIi4vjlltuQSklw0QiYiRwEUIIYbDZbNzx/o/x5uuTcGzdyrx580hKSuo4fvr0ad555x0mrP5Er9dyOp1D2dRhpGRWVJSQwEUIIcRV7HY7y1bciNaavXv3UldXR3Z2NuXl5cycOZPVq1f3OYF3NFAKLLIGS1SQwEUIIUSX2ntLFi9ejNYaj8fD3LlzsVqtstGsGDESuAghhOiRzRb4qIj1dVgGSxJro4MELkIIIURvFCgJXKLC2BicFEIIIcSoID0uQgghhpTb7aa+vh6Hw0FCQkJMDjkpZA2WaCGBixBCiCHh9/vx+Xy8ufk1bC//liZHEku/8B3S0tJibzaSrJwbNWLslSOEECIWuN1u2traeOaR7zF9w0+Y23iR66uP8dqv/gev19v7BYTohvS4CCGEGBS3241SCpfLRUtLCzU1NZw7dQL/kR3cdmUvduXHZw98T76h9A1eeXYKq+68O6ZW01XIOi7RQgIXIYQQXXK5XCilaGhowOPx4HK5AEhPTychIQGr1YrFYsFisfDkd7/O5PqzxGsPqf42luOiq4/5BL+HiW/8nqNTpjJjxozYyXeRoaKoIYGLEEIIg9vtxuVyseWV50it309ughuHRWO9Uo8GznrtnI6/hpsf+BhJSUlcuXKFmXWnmOSt7bhGT1OHizz1bH/8h+R94d/JyMiIneBFRAUJXIQQQgDg9Xrx+/28teV11Gt/ZNUNTmyJncd9tW0A5NvbmKb3sPlHF0ld9R6q3niOm0KClr6YX3GCA//+RVj9bpYsvwGHwxHlQ0eyV1G0kMBFDAmfz9exjb3P58Nms3WsvimEiD5ut5uTJ09S+vTjLL5yiETtwWa5ptv6FgU31L1D65PfZo7uf7KtAuZdOYH7t//GK68+w+LPfJX8/Pyo3UpAIQvQRQv5JBER5/f7sVqtPProo0yYMAGv10t2djbz5s2TLmEholR9fT2Nf/oet/ivQIoFcOJ665xRx5qbaJSVRZGAL7ADod0MOHz1LrNuWGJrexDgxM/K8sNs+cWPWfflb0dt4CKih/R7iYizWCz4/X6Ki4tZtWoVt99+O+PHj+evf/2rTIMUIgr5fD4OHtjPHP+VEWvD/CNvsfGZp3C73SPWhh4Fk3MjcRODI4GLGDCPx9Ptsfr6enJzczvK2dnZxMXFobUejqYJIfrB6/Xirq3uchbQcInzech97jeUlpZG7Rcci1VF5CYGRwIXEVE+nw+3283rr7/OvHnzOh5vamqioKBAuoGFiFKWlHQ2Z1zLluRZ1DEySbIp7mZeeOEFmpub8fl8I9IGEf0kx0UMSFVVFTabjbS0NCAQsFitVjZv3ozL5WLVqlVGBv6lS5cYP368ZOULEYWcTic33ngj3HgjFouFl576M8sP/YF4eu750MBlazIAubRio7NHtdqawP4Ji0meuxS/24WnuZGsnc8zo7W82+vtXnYnn3/48/zpT3/ivvvui8SvFjFK1nGJGhK4iAG5cuUK48aN6yhbrVbcbjctLS2sWbPmqmmNjY2NFBUVDXczhRB95HQ6O+6vufteni67xOST28nxNNJkddJsdVJfZaElPgntTKCw4gznpyxm3vs+hsViYcuOt5i550mStJvdqdPIvf89rJszB4vFgtYav9/PpsRUkm6/HQDLF9d2PJ8GTpVZmX/rGrTW3H///fh8vij7oiP5KdFCAhfRbz6fj7KyMgoLC43HHQ4Ht912G5s3b2bVqlUA7Nu3j/LychwOR5S9CQkhuuNwOLjnk5+luvoDlJeXk5aWRm5KComJidhsNvx+PxcuXKA4JYX09HQsFgu5uffwDFayMjK49YYbsFgsVy2BUFJSQlVVFVlZWR2PubHw+vjbmPvxuygsLESpQHAgyyeI7sgrQ/Sbx+OhsrKyy0DEZrOxYMECnn/+ec6cOcN73vMeZs6cic1mk/wWIWKIzWYjNzfXSLJvZ7FYmDBhAtrvRwXfBxwOB/fffz9ut7vbheQmT57Mnj17jMBlW+pM1jz0CeLi4obmF4kQGSqKHhK4iH7xer1oramsrOz4ZhTKbreTnp7ObbfdhtVqRSkl35yEGKVUF19eelv9tqGhAb/fT6WKp17FkbriLlnfSfSLfKKIXrlcLs6ePUtTUxMul4umpiYWL17c7dCPrJIrhOiK1ppLpaU8+/vfMvMj3yUjPp6ixMSY6Y3tKlAbsudSKgn4R2AhsAjIAx7XWj/Uy3k3A68Gi1O01qfCjtuALwEfAfKBs8BPgEd02HoV/ak7nOTTRfTKZrOxc+dO7r33XiwWC3FxcezevbvLHhchxNjk8/nweDy9DvmsWr+e3Nzcjp6ZqF1wLpxiuNdgyQK+AZQBu4H1vZ2glHIAjwDNQGI31f4H+CjwC2AnsAr4MZABfHsQdYeNZEuKXvl8Pm644QaOHz/e8abU2toaM9+ShBBDw+1209bWRmlpKa+/vpm33nqrx/VXAkm8ucZwUnRvrDiiyoBCrfU44F19POcfCAQVv+jqoFJqLoFA5Ada67/VWj+qtX438Gfgy0qp/IHUHW4SuIheORwO8vPzKSsrA+D06dP4/X5ZBVeIMcbn9dLW1obP56OqqopXX3yeNz//YVq+8Qn4/f8yf/78Hr/QRP8O0N1TRGa5/74m+GqtXVrrS31un1IlwFeBLwL13VR7IPjzR2GP/whwAncPsO6wkqEi0SftSba1tbVUVFSwdOnSmH0DEkL0n9vt5sqVK+x59s/YTx0mp6mS2Q3VOPyBHpZTM64nKSlphFs5hKJ/VtF/AweBxwgMMXVlEVChtT4X9vhOwE8gn2YgdYeVBC6iT5RSKKUoLy9nxowZErQIMYb4fD601hz553/i2vLTHY9rv6Z9xzJLWhZ+v39kGjjGKaXWE8iBWay11j3kH44DrurF0Vq7lVLVQMEA6w4rCVxEv/j9fknKFWKMsVqtHDlyhOLasi6PV8cn05CcNryNGgER7HHJUkrtDin/XGv984FcSCkVT6C35VGt9Z5eqscDDd0cawseH0jdYSWBi+iT9mBFay2BixBj0PTp03l26rUsPbjFeNyP4sw9H+OB++4f1e8NCrBE7ver0lovitC1vgKkBX/2ppVAfkpX4oLHB1J3WEngIvpMKYXf75el+4UYg5RSXPfxz/DO0dtoOnOSjB2vMqnqAlWJqcxbslTWbhoBSqlxBNZ6+SGQppRKCx7KCP4sVkr5tNalwfJlYHYX13EAmcHjDKDusJJXmugTi8VCfHw8LS0to/pblRCiazabjcLCQvLz8/GuWEnpyluo/ManqI1L5prs7JFu3tCLzuTcHAK9Il8I3sK9SmCGUVqwvAe4TSlVrLU+H1LvWgKzjEOHmvpTd1jJV2fRJxaLhQULFnDs2DFZnluIMcxqteJ0Opk8eTJH56ygfNykqN9nKBIUCqslMrcIKgXu7+L2ZPD4p4EPhtT/U/Dnw2HXeRhwA08PsO6wkh4X0WdWq5UHH3xQelyEECiluO//fYnGxkZZ02mIKKU+RaC3pL2TYY5S6qvB+89qrQ8SWBAu/LxZwbsvhy75r7Xep5T6FfB5pVQynavhvhv4ltb68kDqDjcJXESfyRRoIUQ7q9WK1WolIyOj98qjxAgMFf0jUBJSnh+8AVwksG5Lf30COA98GHiIwP5DnyGwlP9g6g4bCVyEEEKIXigV0VlFfaK1Hj/A874JfLObYx7gW8Fbb9fpc93hJDkuYsi53W48Hk/vFYUQQoheSOAihpTb7eaJ3/6W48ePy6qaQoiYNpx7FYnuSeAieuX3+wfUY+J2u9nw0jPcc201Z45ux+fzSfAihIhZErhEBwlcRK/cbnePO752d85bW7dwbdE5UhIsrJxWx8t//TknT57E4/Hg8/mGqLVCCCFGMwlcRI88Hg+///3vKSvreo+SrrhcLg4ePEhJ3H5y0wMBT0qChfWL2khpeIaX//wIjY2NeL3eoWq2EEJElFJgsUTmJgZH/glFj2pqarjOvYNdb7yA2+3u0zlaa86fPsCE/KsnreVlWLl9nos9r/2Kt7Zupby8HK/XK0NIQogop7CqyNzE4Mh0aNEtr9fLkUMHWZbmIjW3jMOHDzN79uw+rZybkFqAmvA3AOgrvzKOWbxubl4CLvc7HD6yl+0VDm6940MkJibKPkhCCCF6JJ8Solt+vx+vqwWbBcalKxpPv0xbW9ugrun1arYfauPF3clsO30Ncbmr0MrKS08/LkNHQoiopZDk3GghPS6iZyF/Y8sme9j4xqvcfOtanM7udjvvmtaaY2daOX0ll5tueQCbzYbNZsNqtTJlyhSeffopjh07xvTp02UvJCFE9InOTRbHJAlcRLdsNhuOhFR8TYGuOatFQUvfknQ9Hg+1tbWcPHkSa0sTF6sTWbTkXaxemH3V1gFOp5M77ro78r+AEEKIUUcCF9Eti8VCbm4u9ZdsZAUfi9fVWK1WPB4PlZWV1NbWUlhYSGJiIhAMdhwOVqxYwdGjR1mwYAEu13TmxMVhs9m63aBxsLvLhiYOy55KQohIUwz/kv+iaxK4iB5lZmZSlz6OrKRAYDJnqp8NLz1DUkomCf5DTCx0cuGImwtVcWTmzWH+/MD+XykpKSxcuBC73T4sgYTD4eDo0aO0trYya9YsCV6EEGKUksBF9Cg1NZWzzZ3l9CQLq6aWQlIdNlsgmMlItTFrVjJ19SfZsnE/9sRrWLJkaZfX83q9eL3ejhyXSMrJyWHXcz9kv8fN/AULJVdGiF643W5aWlpISkrC6/WilOp3/tpYIjku0UECF9Ejm81Gg8ccxrFZFdiu/gNOS3Vw63IHLZ4a3t76G5pdWWRmjaOoqIiEhATKysq4cOEcfm8LJROmM3HixIj2jCQnJ2NRcPLwDhZduzhi1xViNHK73Wx74w2qfv8rLOlZWLKysU6cwpp33StBf1eUzAiKFhK4iB55PB489gygude67RISbKy4PgeA1taLlF05RU2Fj7yceKYud0JrC8fPbOGdd9qYMWNGxIIXv99P/MRV3H/ddfh8PlkTRogeNDc34330RyysON/xWE1iKqfnzGPatGkde4tJECOijQQuolfKmUl/ApdQ8fE2Jo5PMR7TwNSJSex9ZwfHjlmYNm0adru928TdvoqLi+O6666T/BYh+uDixYvk1VUaj6U317PnzdfJy8vjzJkz1NfXs3z5cgleCCTnWqXHJSpI4CJ65HA4cCRm4E9sMLpJVVKGUc8fn9lxv8VnHnNazKDHbg287BZel8mFi8d56fm9zJ63gqKiogG/QXo8nmFLBBYi1vn9fqqqqrgmXmOJ7xwK9nv9zH3pCQ69+SJZLfXE+TQvXv4Ya+9/QP62kFlF0UL60kWPLBYLObm5NDQPzW7ORYXJ3LHMh6XhZV5+8aku90Pyer14PB7a2tpoa2vD5XIZq+y2l/uzEaQQY5nH48Htdnf5AZDoaWNibRkprhZSXC2UPPUY+/bu6fNeZUIMNelxEb3Kysqi9rKXtOShebkopSgZF0ebu5KKigoKCwuNYSOLxUJNTQ1//vOfKSlIRVnjsdnjSUhMJjExkcsXz2Bzn6VFFbP29jvkm6EQvbBarWit+1Q3q6WeI7/8ES3Tfjam/7aUrJwbNSRwEb1KTk6mfIh6XEJdUxLHS9teRKl15OXlAYFZTRaLBbvdzsKFC1lUsAeoB8Dj1bQ2u5k9S6GUjfMV53hjy2ZW3ngzVqtVknOF6ILWmjNnzjB79uw+nzPt8hm2btzAbXfePaanS0vgEh3knV30SmvNcPy5KqVYe70N15Vn2fTCo+zduxe32015eTmbX3uOBXOqjfp2myIl0dLRO1Oca2dRwVFefP4JqqqqpGtbjAp97RnpK7fbzYkTJygoKOjzOVbtx/bK07IRqogK0uMieqW1xpKQjEruTLptcxQZdRraxnXc92kzwbbN0mqUbZYso5yW1LnjtAImjm9m4nioqX+H/XuOkpkZz/q1WVitFsidZDZOmbF3mruVOydo9ux/kbKEeUyfHrnp1kIMl/apyBUVFVy+fJl58+ZF9HWcnp4OQEpRsvG4p8UMTPzezp7W+d7LvPT733DnBz/S0RM6lihZxyVqSOAiehVY6Xb436QyUu1kFmX3+zylFIvmZ7H/0AFOn3YyadIklFIyfCSinsfjAeDo0aOcfP4pCvZuJdHVyovv/hi3P/DeQQcvXq+XPXv2MHfu3H6fa9d+Fr7+R55vbiR30fUUFRXhcrlobm5m8uTJJCQkjPq/L5lVFB0kcBG9amxsJDkp9tZxmDc7k1Nn9rJp435s9jRsNidxcXEUFRWRl5fXbSDjcrlwuVw4HI5Bb/4oRF+53W5OnTrFhQsXyH/kX7m2tbHj2Li//prj8xZ0rHk0UFpr6urqSEpKGtD5SV4XN2x/Cvfbz1LrSMLqcZPk9/Hi3Ju4+/Nfkt5NMSxGd3gsIqKuro6U5NgLXAAmT0zl9luSWbXCx43zy1k0pRRfwyu88OyvuXLlSkceTPv00MuXL/PKK6/w3HPPyaJbYlhpramsrGT16tWkhQQtALkN1Zz/yX9SVVWF1+vF7XZTU1PTr1wuj8fDyy+/zNKlXe8j1h8Ov4/ctnrSXc0kedq4dt+rbHzu2VGdV6YIJOdG4iYGR3pcRI98Ph/V1dXMmhT7H+JKKex2RUlBAsXjNIePP8P+/QVkZeVRW1uFy1XLlCkLWL9+PYcOHaKlpYXk5OTeLyzEIPn9fs6ePcu0adO6rTPz1EF2/etXSXzX+6jbt4vswzvw2B3Y3/93LLl+Wbe9HX6fD6/Px/O//T+W7HkK51s/pSXC7Y/3ucl57jHOzJzF5MmTI76BalRQsnJutBiFry4RST6fj9amC7ic5ph4Zes1RrmqtTOJL/6qfBhzuCXJYc6SiIsrMcq2OFfH/XqPOfOhwW3mvDgs5jc8pfxGudXb+dxZCReNY7OXXGKGz4/LVUd8fBxKjUMxFYC8vDzKy8slcBHDwufzcezYMe68804AtN/8G7E6A2/VC84cwPf9g9hDNjmt/+1/sLnp77hp1Zqrghe328358+c59PTvWH5mK/HaC3ROZ7YVmq9vq8f8+9Etnh7b7Q85PgMXOy5dorCwcMBDUUL0hQwViS65XC7cbjfHjx9ncom/9xNilNVqISHh6n2SsrOzeeedd2hpacHnG/o1bMTY5fV6OXToEEuWLOnTfl3WsOnROa5Gxj/1U1544tfU19d3DCWVlpby/C9/hvu//4lbT78eDFqG1sqVK9m0aVNHkvFooggk50biJgZHelzEVVpaWtj86gYsqpHpk2B8UQqtvZ82qiilWLduHa+88gqLFi0iPT1dcl7EkPD5fFRVVbFgwYIBXyPL3cSy1x5n/97NNOUVE3/lPHmuBm50N2B3WiPY2p7ZbDaWL1/Ojh07RuGGp5KfEi0kcBEGt9vNa5teZPUNYLeP7e5em83G2rVr2bZtGykpKUyfPh2/3z/K3ozFSHK73bz55pssX7580NeyoplVdw5L88XeKw+hrKwsEhMTqa6uJicnB6t1+AInMTaMeOCilLICaK2lP34E+f1+PB4PVVVVlGSWYbf3fVXN0eRSc4NRLkhMYdmyZVy+fJlXX/sj6HTmzJlPdnZ2xyqiNpttdCYjiiHX0tJCSkoKVXctG+mmRNT8+fN54YUXWLVq1egJXBSM8mVqYsaA3m2VUjcAbwSLRVrri8HH84GHgWuBRUAq8GGt9WPdXOe7wfo+pdQPtdZfDTueAXwBWA+MB1zAJeAt4Ida62MDab+4mt/v59lnn2XK+GbmLhhPk32icby6bfJV55QkdSb5FffQObP1cr1RbnSbf/02ZSbn6pDUq11l5lh5Rc0VozxtXIpRdljNaze4OudPnCaDVSUZIUfN3yk8aAk1btw4cvIm4fP5KS3dzaG9l/DpRFDx2BxJJCenkZmZSVpaGmlpaVgsltHzhi2GhNvt5o033mDdunVcCjvmTDcT2pPHmX9gluQeev3ChjOuSsAN/RsIX6eorc0o+i7WGmV/vcsoq4Suh0+VUtx8881s3bqV5cuXj4peyvYcFzHy+h24KKVswCNAM5AYdngq8EXgDLAfWNnDdd4DvA/4RPCh7yilDmitnwweTwN2A9nA48Ch4PPNAO4HdgASuESIUoq5c+di9e3uU4LgWGW1Wpg8OZupOU3BRzxALa66MzQ0+6k742PnRRtLbnwPmZmZI9lUMQQCq0h3vm26XC6UUv3+YPb7/Vy6dIkZM2aM2gA3ISEBh8OBx+MZFYGLiB4D6fj6NJAHPNrFsT1AttZ6EvDVLo6HWg58T2v9uNb6ceA/gRUhxz8KTADu1Vp/Smv9v1rr72utPwrkA88PoO2iG1arlfHjx3P0VBxe7+idRTRUnA4L2ek2Jhc5WX2dhf379sqGdKOMx+PhwIEDeDwePB4Ply9f5uWn/0pjY2PvJ4fx+XwcPnyYyZOv7skcTRYvXswbb7yBy+XqvXIMsKjI3MTg9KvHJTgU9E3gH4DC8ONa60agr3/Fp4EPKKXah5weJNCz0q79L/oNwmit3UBVH59H9JHD4WDa9DnU1O4jIX+kWxO7rBZFW9M5/H4JAEcDv9+Py+Xixcd/TubxXWxe/SCt7+yh5MweFrsbKZ00pV+9a263m82bN7NixYreK8c4h8NBdnY2ra2tOJ3O3k+IYjJUFD36O1T0PeAk8Cvg64N87p8B9wB7g+WNwM9Djp8N/vwI8NNBPpfoA601tbW1FEyxM1ozpX998LJR/uCccd3UHJyijEaam5uli3yE+f3+QW385/F4aGxsZNu2bazY8wx27Yc//3vHcQ1cOnqIOXPmYLVaex32cbvdHD16lIzHf0DTT6pp6rH26DBv3jxeffVVbrrpJvl7EBHR58BFKbUSeC+wXGvtH2wehNa6NXjNmQT+/o9obays9CjwWeARpdQ/AFuArcCLWuvyQT256JLH46Gs/DSLrs0HqlAhKUqDWQhz+bhUo/xCabVRrmyBdRM6v7E+d6azM622yVwZt7HFLB+6UGeU453mS9ph7/wg8YetRhquIDGlx+NlrbM77uelhCUYc9ooTyvy8vLrG1iz7h55sx4BXq8XrTVHjx5l5syZA8ojcbvdXLp0ibNnz7Ju3Tratl/9/Sn+X19kcVkZGzZsICkhgWumTSMnJ6fLNX98Ph81NTW0tbVR2Gj+DSTkJJjXDSvbJ5h/Q5akzteUJS2sJyMsUFP5eWZ55sKO+7qpxjy3wgzsrWHJu9YaM1nXV2Mm83alfdbdaNg5WnpcokOfXkkhCblPaK23RerJtdZ+rfUhrfXhsKAFrXUVgZlJPwUSgA8DvwQuK6X+Tykla7FHiM/nw+12c/jwIWbMiO3u3GhhtynmF5dz6NChUbmKaDRrX/H55e9/i+Yn/4s9e/bgcrnwer20tbX1uhGgDi4NcGD/fip/8Q1uuummHhPW8/PzueOOO5j/p39lx+9+2W2bKisrOXDgAIsXLx7U7xeLJkyYQGVl5Ug3Q4wSfe1x+QxQAqwawrZcJTjN+pPAJ5VSxQRmKT0MPERgyPGh4WzPaOXz+Th//jxNzUeYv6Ck9xNEnxTm2Ck/toWLFzMoLCyUlXeHgdvtZtfOnaQ/9UNu9tRhSXZQ9cx32f1cCo2JWfiyirHlFZOVX0h+fj4ZGRnYbDa8Xm9Hj0BlZSU7n/kD80q3kutv7vNzW9Bcc34/Z8+epbi4GKvVis1mw+12s3f3blp/89+s+unvx9ysPZ/PR21tLR6Ph/z8GE6ek8TaqNFr4KKUSgW+QSCvxaGUGh88lBb8WaiUag8yhozW+jzwG6XUH4EjwHuVUh/Vehg24BjlHA4HxcXFvPPOAdraPMTFyQdspCya5mT74aex2+8nNzdXgpch4na7aWxs5M0332TWlsfI93Vmj2TRSpa/FRor8F/eBwehRdmptiZQ6kikISkbnTceW24xrrpqig+/wipvDQP5jCpx13Lo5//CielLcIwbT35hESd27WDqq39kfFvjmAta6urqeP3111m2bBlpaWkj3ZxBUYB1jP3/Rau+9LikA8kEejoe7uL4duAcgQXihpzW2q2UOkBg1lEWIPkuEeBwOFi9ej0bN/6R29cVX/UGe7S2xSiXNV09vfHmonQAfrbz3FXHPrG4556c72w60XF/9sSMHmoOnYef2GOUF8wycwNumZw2oOsuneVg664ncM25j6KiIsl5iYD26bVWq5Xq6moO7tlFwptPckvbZSy+nnOZABK0hwRvPcWqEWrLofYQHO3muf7rnm6v0/rN9Vc9NrvhPLx9nubqVpptcSzxto2Z3Ww9j7634/6Vgw3sveEDrLv7XgnYRUT1JXC5Aryri8ffAzxAYL2VC5FsFIBSagmBhN3GsMfTgaVADSCDphFks9lYuPA2du/ZzvwFhdiH6N02NBEXzKAFYM/xzv/WxDjzJZoQb74Bhq85c7oybJXehs7kQbuzi5d7P2YVldZ15qo0uWcZx3ITcoxyaqqZ5HhDejnHTr3Jq5syuX7ZCuLi4ga0cNlY5vP58Pv9VFRUcGjX2/jOHiP+8inyaGKhtxFncC6cJaMzodSWa66RGZ5MakkO+0ANSyDVrs4O3dCEWAC/y+zsdRSnXFXO/8ozXf4uBc/t6vLxdt5fvc9s1jXmFhwqt/P1pqvNBFuVaibyEjYtX5d1JpKrnGLzWL65yoVKiDePh/37WHydifShq+r6NOyYdxd3jrKgRYaKokOvgYvWugV4OvxxpdS84N0NocNESqn2hefa/yLuUEq1/zX8WGttfrJ070Hgg0qpZ4BdQFPwmh8CxgGflP2NIstqtZKZmUlpaS5VVU3k5/R+jui7aZNTmFTi4vC+/6O8PhWLI5P4xDTi4uI6bu09MqN1NdWBcrvdHDp0iAuvPs+kiwe43lXbMZRjTZB9oqLNibZkbn7w/lEVtCiUzCqKEkPxF//PYeV7gjeA3wJ9DVx+RiBYuQVYTWDIqo7A6rwPa62fG3RLxVUcDgfXXruYjZueJH/tSLdm9LHbLcybmkBgq4By/P4yPNU1uL3gqtU8vzOVux/4WwlcwrS2ttLyf//BiqbL2LrqORNRJcvm4s2Xn2f+8hvJzc2VnkURUQN+B9Baf5PAKrrhj0ckJNVaHyaw75EYZkopMtIn0NDQQEpKz2ubiMGxWBROu8Jph+R4xYycOk6dOsXEiRNpbGwkNTUVv98f86uO9oXP50MpddV6Hx6Ph11v72BWc8UItUz0V7bdzU1vPcaJXc/T8DdfZ/r0aVgssR+MS49LdJCvLuIqNpuNBQsWsumNl7nuhvnBR3tf0vzH20oBsNt6To753J/2G+XsjISuKw6D9/zn6x33c8aN/NJAM4qsnK98igONt5KWlsaRI0dwu1q48aZbR1W3e6j2JfV37NhOUlIyc+fOxWKxdGxmqJSi9a0NOGRkOKYoIN/bSKXXixrQHK3oomQ6dNSQwEV0ye/3oz2df6WtYUmwzR7zQ+RseWcOdXNr/xZcK68wt7dKS+tMCPSEzRA5UWqu3Olwmt/iLGHvLI31ncmYbpfZZl/Y71RZbi7Avrm21Cg3L+pMXEwMSxIuTDd7popSco1yfmJn/npCgplMaUu6ZJRLCqEkObBI2TXXXEPd8Z/w/NNlrLr9AeLi4kbVMJLb7aalpYXN//ff3Fy1ExdWNm26hlnv/zTjxo3DZrOhlCLBXYczPZB0a0k1e59shWbAac0KCYTDElMtWWayqQr7f6MtbLacrfPfWteba7roNvP1dFWi7yDYPvI787mOfMeskNC5lLUld7xZ1xu2wF5rg1l2dP4btCTOMC9rPWley2r+Tsod9rcdspGoajG3j7Mk2Cm1ZTN/yhTUKFg1V0QPCVxEt/zyDTcqpCZZWX9dG69t/CXFU29jypQp+P3+mMwb8Hg8uFwuzp8/j8vlouzcaZIObGSV+wIWwIGfm5qO8NoTj5D7+W9js9nw+Xy0ZRbCpYhPXhRDrLZgYNstRCsZKooOEriILtntdqw2+SONFjarYtW1dsqrN7L5mZdIyL2WxdctibnhI5fLxZZHvsGCplLs+JiKt8tBhMW177B7104WX7cEgMT5N9BweQ8puufl+kX08GAhafqiURa4jHQLBPRxryIx9rhcLpLTBrGzohgSeZl2bp2nmZq4g+efeozm5uaOxdiincvl4s1XN7KSUtJwkdhN0AKQiIe6N5/vWOtm2Yob2THnbk7bMnDL21ZM2O8cx4xZszpylYSIFHlFiS4ppWis7/s+Lb35wI/e7LiflT/ySbDRzheyl2lXH9OZKRbWz6tn+0s/ojluGlm5RRQXF5ORkYFSKuo+LPx+P83NzaQc3YDdovH3fgrjGi7Q1NREWloadrudde/7ENXV6zl85gwNZefxV5SSUXWKhbqVaO3B97/yqc5CyP+Jys4y6l21iFyiuXAeYQvBRbsrbgep936K9PT0kW5KxChkqChaRNe7m4gaNpuN8QWTuXS+jILifPzm5t3kJzlZlBMSgEzsfCP+6M93GHVdYcm6R3aZ21pl5HTfs9PSbJ5bc8VMoK09by4L1HK2zigre2c39bil5qqgDVXmNgY128122cOSQJsbOocpcgrM4Ot4mvnBsmKuuZmcy9e58mmq09z+INlhbi2Q6jBX3fXlzyeczbIEO7DimkC5paWFsrIy9u3bR3p6OgsWLAj8DsM8lNTW1obT4TCSMT0eD83NzWx98meszWtCqUQcq83fyX/stFHWSnG5bgKTQq7jcDjIz88nNzcXj2cRAE1NTbzy/C9ZXXAZaxcRnkoOC5LrzdeLmtDLpqKhybqOK+a5YcGhSk8zf4eyCO5GklvY7aEr9hVG2eo0h9N88WYu1LGaztf11tfN3+kjy8xr5SVsM8o6zqyPI+T1Ze38UD/vTmD+9OkxN5TZM1mALlpIn6vokt1u55prruHcqUu9VxYjLiEhgUmTJrFmzRqmTp3K5s2b2blzJ263G4+nf7O8+svn8+HxeDh79iwbN27E6+tM6na73Zw4cYL9T/wLa5zv9LlnpMFrpWDxXSSG9zwAFosFp9OJ0+kkLS2NFe/6BK9fSkf3vkWRGCb1lrRRldsioosELqJbXq8XpeXNJ9akpaWxatUq5syZw5t//k+2vbUVl8uFy+VCR+jTXWs/bW1teL1eDh06xIYNG2j90eeJ2/pUxweW2+3m7W1vwaNfZ3lSeb8SG1NsPs4f3Y3P1/XMNrfbjcvlCkyVTkhg9uqPs/FijgQvUcKflIfX6+29YgxpX8clEjcxOBK4iG7ZbDacjviIfdiJ4ZWcnMxNE5tYYHmFbX/9DzZvfIm33nqLY8eOUVFRQUtLC36/H7fbTVtbGy6Xq9tAAQKBbPu6KwcOHOTVP/+ON998k6KiItavX0+Rp54MTzP1weEYq9VKxcl3GOdp6Paa3VEKprXtYdfOnTQ3N6O17mhfU1MTr214kdeef5ojR46gtSYnJ4fcOWupbJG3tKgwClbJ7YpFqYjc+kIplaSU+qZS6jmlVJlSSiulHuui3kKl1A+VUgeVUo1KqXKl1KtKqVu7ua5NKfU1pVSpUqpNKXVMKfUppa5uWH/qDifJcRHdstlsJCYk4vXKei7RpqotMISXGXfKPFB9/qq6iQ5YOaEFeBtveRX1l200eKxcclupb7Xidqbjj89CxWdgT84gMyub+Pj4jmEmt9sd6OFobaHt+H4yLhxmoqeWKX4X9slpcAjasypKPLXs2fk2K266GYfDwR0f/Civ/Wcpt9H/fI9JSa3kXvgVh1700mBNoO2axZTMWcj5V/7ITQ0HsaIp35/AK9OWo+OyyW49RHZ+X9J+hYgJWcA3gDJgN7C+m3r/RGBPv78APwGSgA8DG5VSf6+1/p+w+v8DfBT4BbATWAX8GMgAvj2IusNGAhfRrfZvuRarE6+/770uF06bMySazpir3eqwVXfbysyE27J3Kjvu+1rM/Ay/y+x+Dj+3x3a9aq6E628zr+XMMfMpXFfMWVVnnjnWcb8uLNF34ixzpdwDZ6qN8qXUzuTd9BQz6TczPs0op8WZM07awlb4nZxmrjTcEzVlBirzAx1l244vkUnnBg7K7gC8QDm69ihaQ/NlhcunsLe2YLNobEpjs2iobQ28JU4HSMXfZCaB2vKTSALydj7BNruDOfPmcebMGVLi2rBkmv+2atIso2yJizMbHlzxNgVYnB4Ixvx6N7Wl+5iRV4fKD8xWKQKK1R5wEeg/rgDGdf776dZW83lzw7Y8bwqbOVc4wSy7Oo9bJpptDrej/kajvGjmO0bZ1tKZdK3bzNetKppkXixtnHm81UwqbrV3JhXvumC+jZ+6ZP5OqclmL9Sf/7K/437TJfO1NHfSTUZ57XhzVpCKM7fnUM6Q/7fQVa6tozOAHOZhnjKgUGt9SSllI7Aza1f+G3hIa92xTLhS6n+A/cC/KqV+obX2Bh+fSyAQ+YHW+vPB6o8qpZ4EvhysW9bfusNNAhfRLZ/P17HxnRgblIIkhyYJDWHDRn39KJpsaaBw7/9wcVcixZZWUi0eIHvQbbMoyIzz4e97rCpGiGqpvWqzzFg33NOhtdYuoNfZEVrrbV081qqUeh74PJAHtE+ZfCD480dhp/wIuA+4m0AvS3/rDisJXES3ArNFZKVS0X/xyscUW/9zW8To4GytkllFI28cge7U0C7vRUCF1vpcWN2dBL6bLBxg3WElgYvolsViwe2XwEUI0T+FzmbKy8spKCjovXIMscZI57NSagZwD/Cs1jp07HAcXfTiaK3dSqlqoGCAdbtrRzZwAzCLzm7XSuAI8KbW+kp35/ZEAhfRLaUUFtmvKPZ5Xh7pFogxpjiulU1vv0X2+rtjcjPQrgSmQ0fs/TBLKbU7pPxzrfXPI3FhpVQK8CTQAnwu7HA80F1XaFvw+EDqhj6/DXgP8DfAiq7qhNR9k0Di7x/b83D6QgIX0S2/309WSg6+OsWSvJQ+n1e721z5NTyhNpyn/uq9dnZt/lifn2+oLLjup90eC19lt+mkmZBcs6zYKIevtGscC1s5OD3ZTN7NyTATIh3W1I77mmuMYylpZhKs019llFVJWPKps7O+CjtEmdlDbA1LZLWEJb7S0JnoqVvNnjqVb64OrN3muQs/ePXqyXvf/vvAnbgfmM+bZL6+dFmlUaatI0cRNW+BeSxvinnuthfNdtaZXwBd42/uuL/xXJFxLHw16T8/f8AoT5lq5vXcPGc63Xn7iPk7rJxuJnuXpIQl84YkbM/ONl8fP/yxmfKgwroJxoe0a8m7ZhrH5uWYK91avOb/ufabmU46ZCViS3LnuRagePdTnJ01j4kTJ0bdFhRRoEprvSjSF1VKxQPPAROBNVrr8GmGrYDzqhMD4oLHB1K3/fnfC/wz0P5u0lukd0Pw9i2l1Ne01r/vpT4g67iIHjgcDhYuXMiOHTtkLRchRL9M8tdx/g8/Dq7DMzpmGQ3nOi79pZRyAE8BS4H7tdZbuqh2mcAQUFfnZgaPD6RuuycIBE0KaAZeBL4FvB9YC6wDHiQwlfqlYB0VPOc3ffg1AelxEb1wOBwsXbqUw4cPM3v27JFujhAihixpOcUrv/4Z6z/+uVEzZBSNgsMzfwJuA96vtX6+m6p7gNuUUsVhvTHXEujI2DPAuqE2AP8LvBScGdVTu+MIBDSfCLa9T6THRfQqPT2ds2fPjrolvIUQQ8sC3HDpLd7esT3m3z8C06Gjb8l/pZQF+DVwF/AJrfUfeqj+p+DPh8Mef5jAOpJPD7Buu7la67Va66d7C1oAtNZtWuuntNargbm91W8nPS6iVw6HgxUrVrB7926WLFky0s0RY0SD5wQA3WcHiViQgIfGygp8Pl+M57oM/+7QSqlPAWl0djLMUUp9NXj/Wa31QeB7wHuBLUCrUurBsMts1FpXAGit9ymlfgV8XimVTOdquO8GvqW17hj+6U/dkHMODfR37c+5sfwqEsMoLi6O+vp62traiAtf5TTMzlc/Okytih7uqhajHLrKLsClkFV5HVlmMuWpLDM5PzPP/Kh+z90zjPL+i3Ud9+uyzGTc/CSz7uTUw0bZnmOureFRnYm+9W5zdmNmwRtGmVZzlVXlCftCVd25WvBVb+/hKy831xnFTW+sM8o2S0hyb1hCKC1hq+GGJZ/6SjvbYal/0zhmmROWrDxjXti1zOGMd2o6k3ndPvP/2Os1f6fUDPP/8Z1D5jYH+3d3JnTXXTQnazQcNpOCJz52r1Fu8ZjXfj3k2jMnmQm1OUWpRjk+0fydJk7sXA135SRzleY4m5kkrBvNFaC5YrZTX6nruO85Za6Q7a1tw4uFrUWLuH7lzTid3eV5ih78I1ASUp4fvEFgUbmDQHsG+srgLdxNBNaUbvcJ4DyBbQEeAs4CnyGwlH+4/tQdNhK4iD5xOp2sXLmSHTt2cOONN450c4QQUa7B6uTtObdz+4c+hsPhwO/34/V6sVgssdnzMgI7O2utx/ehzo39vKaHQMLstyJZtztKqeXAJ4EpQDpXf6/RWutJV53Ygxh89YiRYrFYsFgsNDY2kpwsHfhCiK6V6wSOrn4Pd96+Hrvd3rGr96FDh1i0aFFMBi7DveT/aKCU+giBdVq6rQL0e8qqJOeKPnM4HCxZsoRt267aGkMIIQBo1VYOTruTW9au6whaGhoaeOyxx1iwYIEMGY0tXyIQnHR3G5DYC3vFiLJYLGRlZVFZWUl29uA3zhNCjC47Uuew5l3343A48Hq91NXV8eYj/87MZatjfP8ijerzVqMiqJBAj8rTwHeAqh5r95EELqJfbDYbc+bM4ZVXXmHdunW9nyAAcF1p7vI+QGNY3fC/7GfDkivbWjt3t58911yR9uaZ+Ua5rMVceycj7oxRdqrOtuiwDlh3srnKrMP7jtkwFdZhm5nZed/VZh4LGxpQ2SVGOcFmJnbaVEjib0qmcUwVeYxyeLKuNakzodSz39xqxVKzz6w7v94oM9lc3XZBSmeC8oIMMwm22jvNKF+uSjfKE4rTjHKcs/Pf4GKFuRLukTxz5eBHfrLdKPs95gdm6/nOdr8xyXxef5s57XjFXebv1Obq3PX7xcNlxrEPXRv2b+sKWzk3bNNVFRcSiFiDr4ec8R35LNXV1ez+4bcpqC1j4rz5Mb+Wi1ISuPTTMWAO8KjWurt1X/pNhopEv1ksFtLT02lpaem9shBiTPHUXiEuLo6Kigr2/dfXue7iAcrzJ5OWljbSTRPD7+vBn58M7qEUEdLjIvrNarUyd+5c3nnnHRYtivh2G0KIGJZbe4pNmzahX/od1145gR9wzL8ef/i09hij0Fjw9V5RdNBaP6eU+iaBWUmXlVLHgfqrq+lb+nNdCVzEgNhsNioqKnqvKIQYU+aqSjw7foL7SmDQsyIulelzY3+YCGSoqL+UUquBrxHIc0kA5oVXQWYVieGilIrJKY1CiKFnV52fRefTCygoKOihthjF/p1AB0n7DCKZVSRGlpI1DYbFgd+ZK2E7cztXy01MNqeWjssykzyzJ2QY5VavWY53dvbaZjhLjWN+HZZQm5hmlHVLWI9v6Eq6NXXmscJCs2w32x0XnpLs6UwS1Q1hq7c2mOnMujLs3JAhCUuq+TyeEzVmudRM1nXMvGCULVM6P3BVqrkirWvCtUb5gYUwLiGto/zbsMTXlRM6h/gT5xuHeHXaPKP8i0d3GuVPf3qpefzxvR3373nXTOPYB+deNMqVbblGOdnR2a44i9lrai0z/z30yZNmQ8Py2rSvM0DRHnMYJWFSGvHffB79wgt4vd5R8EVHZhUNwFQCPSqvAz8FamDw/4ix/koSI0Qphc8n471CiO5VV1dTUlLS6zYhsUKGivrtBDAb+A+t9YZIXVSGikS/+P1+mpub2bRpE+PGjRvp5gghotiJEyeYMmVK7xXFaPUVAj0u71UR7KKXHhfRLz6fj5deeom77roLu93e+wlCiDGrpqam90oxQmYVDci9QCnwAeBmpdQuup5V9Df9uagELqJftNYUFRVRVVVFfn5+7ycIIcakVqykp6ePni84SoaKBuAhOmcNFQLdZWlL4CKGjsPhYMGCBWzatImmpiYmT548apN097799yPdhC4tuO6nHfdP7jQTMadPzTLKB6+YiayZ8QlGOcU5t+P+NWknjGNxyvy2rJvrem6YDnlTnzrLPJYWNqzYZCbc6lozkZUrVzqPXTB/R89pc5VdX4WZMGor7kyCVdaw16bdXHK++az55a/2qNmutAnlHfedi81APT/pcaPcVnAjrd5zHeX3jztsPnd9yP+Fy1zt9/7J5r+X98Pm+kiZ8eZU4ue/2Hm+rnrKfJ59ZrJy1u4fGmVLcue11Djz9aLd5sq4/rI6s1xjtttX1Vn2VJnHTsbnMHPmTCwWyUgY41Q399v1ezq0BC6i3+x2O7feeiuXLl1iy5Yt3HjjjSPdJCFElLlSMJ0FCQm9V4whMquo3z48FBeVwEUMiN1up6ioiCtXrtDQ0EBKSsRWcxZCjAKOa+ahdb+/TEcxLUNF/aS1frz3Wv0nfXhiwKxWK/PmzWPfvn29VxZCjCnKYhm1w8hiZEmPixgUrTWtra2jZIEpIUSkeJsb8fl8oyc5F1Ayq6jflFLLgU8CU4B0rs5z0VrrSf25pnzSiEFxOp0sX76c7du3c8MNN4x0c4QQUcJ/uXRU9bgoNBYZKuoXpdRHgF/0VAVJzhUjweFwYLFYqKmpISMjo/cTRMT4WjxG+cjRSqN8+qw5A2dCcZpRXjY1p+N+jWtC2NXNcnZGslG2Vew3q9tCvlnXXzGPnTOXjvc3N5vHw3cOvtI5M6ZtZ7lxqOVCg1EOXXYeQIfMFHKkmEv+h7PGmbOM2i6as7DKQmbKZLV4jWPxl5uMsnORuWWCDptN4z/X+W9inWV+wbyQc71RPld+ySi3usznXubo/DfQh4+Yz9vcZpbDzm07Ejrr6LJxzJJhrnBryze3kLBNSDOv7e/8t/ecN/9f/GWXsVrNf18x5nyJQexJ1B3JcRGD5nA4WLx4MVu3bh1lyXhCiIFSzY2jLnBR+CNyG0MKCfSoPAVcS+DbUPhtYn8vKj0uIiKsVisLFizg6NGjzJgxY6SbI4QYYRbtR2s9qoaLRL8dA+YAj2qt90TqotLjIiLCYrGQk5NDaWkpLper9xOEEKOaJ3vcKNuINTAdOhK3MeTrwZ+fVEpFbM0M6XEREeNwOLj55pvZunUrt9xyy0g3RwgxgmwTrhlVs4oUMquov7TWzymlvgl8C7islDpO13sV9esDQwIXEVE2m420tDQqKyvJzs4e6eaMSqFbESy9w1zfqarcTC71+82co7xcM8G2tq0zubcg2UzMTHWYCaLhY/MqLc8o69aQ524MS6A9d4Ee2cxcCN3W+QERukQ9gC3e/CD0tpoJyl5X57mNF8MSRj3m72B1mm+BFrvZCe1xdV67vrTOOGZPNtthT4g3n+uUuVWBpaQzEdo3f61xLJ9D2CxLOspfvPkaeqKrQ4ZfkswEWsu1S4yy9/+eMcqNIUnFniZzif+sa81tDaxzzQRtVWKWHSWdyb0pqceNY/6mWlnuf4xTSq0GvkYgzyUBmBdehQHMKpJXlYgou93O3LlzOXTo0Eg3RQgxgiae2sHly5fxh88Yi1kyVDQA/06gg6Q90lZhtwGRHhcRcT6fj7a2NknME2IMK3bX8erzf6bg458ZNT0vlrE1IygSphLoUXkd+ClQA4P/R5TARUSc0+lk6tSpXL58mYKC7nYxF0KMZgqYduItzp69i0mTJo2a4EX0ywlgNvAfWusNkbqovJLEkCgpKeHYsWMj3QwhYobLt9u4jQaFnnqObHgGr9fbe+UYIOu49NtXCPS4vFdFsPtdelzEkPD5fLhcLhkuGmLbn/tQj8dXPPgHo7zfYf7J7999iT9/+eYuzz1eV2SUi5LMae62ODPRl+a6zvuWsP/zLHNFZWU3E24B1LX/3Fmo/2Pn41s2hZ0blkB7Pnwl3c7Vbn0ucxaIu8n8HbzV5iqz4Z+voZ0EqslMAq47VWc+7293mucmmMm7ob+x9eQbZrsmd/1/0J2NTes67t9qOWi26wv/Z5QbzpmTOFrCfudQtsNVRjnVZW6g6pxnJmxbbriu89z7phrHEpJfBiCv5tyoWJhSKdkdegDuBUqBDwA3K6V20fWsor/pz0UlcBFDwul0Mn36dC5dukRhYeFIN0cIMUKS7V5aWlpwOnvefkGMSg/ROWuoEOgud6BfgYsMFYkhU1RUxPHjx3uvKIQYtWxKj5qF6BS+iNzGmNAZROGzigbUHS89LmLIyHCREKLeY6MkbK2Z2CS7Qw/Ah4fiohK4iCHjdDqZMWMGFy5coLi4eKSbI0RM0bu+1nFfzb4u7Oj1xIo6nTxqVs8V/aO1frz3Wv0ngYsYUoWFhWzevFkClxHyxm/fM+Bzp6bFhz0y2yiVt5qJnGnpaR334+xh+QyZ4666vsp4f7fPfZI7O+5PWF9iHLM7fm2UdfhquCGrwfqre+6Wd5sLx+Ltobrfbx70X2oyyu6wVWizZ2aZxw9Wdty3nHnTOGbN3WuUbQumdN8Q4NajX+2479px3jhWc6LGKJdfMJOK01I77+cvMFe3Tv/cbeYTjQv7u7WaAYhO7zyuKk+ax4L/L97sWXg8nlGxU/QYmxEUtSTHRQwpn883aqZCCiH6p9ptp3jaXOLi4nqvHOWUrJzbK6XUn5RSiwZw3rVKqSf7Wl96XMSQslgs2GzyMhNiLNrrKeHWqVN7ryhGi/uAe5VS+4E/AhuAw1pro7tSKWUl0IW7GngAmNufJ5FPFDGktNajaK8SIURfNXhtZC9cN6oS88fgjKD+OgFcQ2AzxXnAvwEupdQ5oDZYJwMoBkLHkxXQ5ymoEriIIeX3+0fF4lNCiP6pcDkYP378KOpxlQXo+mAG8BHg/wGTg4/FEQhmQoVGs6UEApxf9fVJRssrSkQpm82Gw+HA6/WOojcwAZAXnxX2SGe5McHMaah1jb/q/J7StaekhiYGLzGO6dnmaq62CnPlXFvISrrabwbN7lbzG3ObuZDuVcm6oS/Z8GP+RvPacc3mxdz1l42yJWTFXxW2snB8lpkInV5qLi7qOH7WKHuOVnfcP/HUGeNYXdi6pEWTzFWKC9/Vmfhrmx22Hlh6ulmuuWIU1bhJZtnfuQqvrjB/34YGPxNGxTRo0Vdaaz/wKPCoUuomAkNHK4FpdObUagK9K28CT2qtN/b3eeSTRAwpm83GzJkz2bZtGxMmTGDcuHGjYnaBEKJn6RY3tbW15Ofnj3RTIkZp6XHpK631ZmAzgFLKAmQGD1UHA5wBk1lFYsilp6dzzTWBnsIXXniBqqqqXs4QQsS6AlsLJ04cHzWr5gKg/ZG5jTFaa7/WujJ4G/Q/gAQuYsg5HA7y8vIoKipi7dq1bNu2baSbJIQYYnalaTp/DI/H03tlIfpBAhcxrKxWKzNnzuTEiRMj3RQhxBCLqzg6inLbtPS4RInR8ooSMcJisVBUVMT+/fs5ceJExxCSGF2S7dPDypG7tkoPy5mwmomufl9n0qwvbFVdd9iX//CEW5fLTLi1hCTRtrb2PDuutdUs19aGX6uzLSWTzYTZunONPZbbnio1yi0tnffD44KpS9OMcs4nFhtlf3ltx333NvO6/lfMLxRxt88wL55bZJbP7O+4qy+Zybm2a3IYP+VhqqqqyMvLI+ZpQGZIRgXpcRHDzuFwMG/ePEpLS3GHf3IIIUaNcePGcfr0aVnLSUSUBC5iRDgcDpYvX86ePXtGuilCiCGgNbz55pvMmzcPi2WUfNTIUFFUGCWvJhGL7HY79fX1krwnxCj0dmUKs2bNwul09l5ZiH6QwEWMGIfDwbJly9i7d2/vlYUQMeNsg530az9Ednb26ErO9fsjcxODMlpeUSJGORwOqqurZWVd0We6qabH435v9+uGhG9UHv4ZEn48NME2vG5vox9XXTukWQ17zdwuhyOsbi/tDP1Tyc01k5PLjtQZ5Qt/90q358aFdYaMv8NcGZe2NqOod7xplFV2duexmqaO+4fc17J20iTs9ghmZUcDGebpN6XUg8AngSlAehdVtNa6X2/+0uMiRpTT6WTZsmXs379/pJsihIiA8jYHM5atkS8iAqXUl4HHgcUENldU3dz6RQIXMeLi4+Pxer3s379fNmQUIsZ5tcLhcIyehNx2WtZxGYC/Y4DBSU9G2StLxCKHw8GCBQvIzMzkwIEDI90cIcQg5Me5OHn8ON7w8a7RQAKX/soisALOH4FCwK61toTd+r15nQQuIio4HA5ycnIoKysb6aYIIQbBqqDt4gHpPRUAR4M/f621vqy1jsjGVTIIKaKKUhHtURSjkMr5iFG2PWiWEyvu6bjfXN5sHEtKNJNiw5Ne/f6wVXj9nR++4avuel09t7OnFI/UVPN50tLCjmeb2bquBrPdNZ2L316VyBv+vDnFcUbZmW6WjecpM/+9fM8fN8r2gmTzBEd5Z92KzuV8J1h3U1lZybhx47p9rpg0jDOClFJJwD8CC4FFQB7wuNb6oS7q2oAvAR8B8oGzwE+AR3RYBDlUdbvxbwR6W+4HXupD/T6RwEVEDbvdTmJiItXV1WRmZvZ+ghAiKhX7Gti+fy/Z2dmjaGaRHu4l/7OAbwBlwG5gfQ91/wf4KPALYCewCvgxgYTYbw9T3a7MAE4DH1JKLQA2A7XhlbTWfblWBwlcRNSwWCwsWbKEN998k4kTJzJ+/PiRbpIQYgCsaFwHtsNtq0e6KbGsDCjUWl8K9nx0uVKnUmougeDiB1rrzwcfflQp9STwZaXUL7TWZUNZtwffIJDjAjA7eOtKvwIXyXERUcVut7N8+XIuX77MqVOnRro5QogBSq44O7r2KNIMa3Ku1tqltb7Uh6oPBH/+KOzxHwFO4O5hqNuT7qZAD3i2kfS4iKjjcDhYtGgRBw8e5OjRo0yfPr33k4QQUaWkrYZDhw4xd+7c0TNcFJ0zghYBFVrrc2GP7wT8BHJkhrpud77Vhzr9JoGLiEoOh4M5c+Zw9OhRDh48yJw5c0a6SSJGOP/hrx33s+rWGcdaXjhjlCcsTDPK1SfM4feq6s4vhImJ5vNUXDHzHXqb/RuaNNvYaJ4bF2d+8Uz2mB+Q6ZPNduandSbYJi7JN59n2VyzXW+ZSwzUvnK2476rwcww9jSZScC2eDPgyJxktkOHJC+3XGo0jiV6/Jz8zY+pH/9fZKSnY7H2e9ar6JtxwFU9M1prt1KqGigYhrpd0loPSeAiQ0UiajkcDqZPn47D4ZBdpIWIQXPLTrDlZz/E64vILNgRptHaH5EbkKWU2h1y+9tBNCwe6G6OW1vw+FDXHVbS4yKimsPhYOLEiVy4cIHt27ezdOnSkW6SEKKPFDDj+NucO3eOKVOmjHRzBkcTyenQVVrrRRG6ViuBnJOuxAWPD3VdAJRSXw/e/ZXW+mJIuUcxN6tIKWUFiNTCNGL0cTgcFBcXo5Ti8OHDzJo1a6SbJIToo3RXExu3b2fixIlYZbhoKFymi9k6SikHkBk8PtR1232TQIi3CbgYUu7N0M8qUkrdoJTSwVthyOM3K6UeVUodU0q1KKXOKaX+qJS6ppvrfBdoAuqVUv/SxfEMpdR3lVJHlFLNSqkapdQhpdTPlFLTBtJ2EZvsdjuFhYVUV1dTVVU10s0RQvSRBbBfKh0FM4yidq+iPUCeUqo47PFrCfzz7xmGuj0Z+VlFwfnkjwDNQFi6Gv9BIBL7C3CMQHLPJ4F9SqnlWut9Idd5D/A+4BPBh76jlDqgtX4yeDyNwKI72QR2lzwUfL4ZBFbh2xF8DjFGOBwOrr/+ep5//nnWrVuHw+Ho/SQxpsX/8wtGOfngYqOcep2Z2JqQk2BeYGtnbmL4564rbPS/rc0sN5mL0BrJueGr8LrNnNir+Fxmh3RoEq2/xnxi3/aDRtlbWm+UW6s7e/jdTWZD/GFJwXHpZrlhd7lRvmBN4fzkeWSfP05JQ9e/hE8pMudGalRkhEXnNgZ/Ar4IPExgpd12DwNu4OlhqNuuPRn3fFg5ogYyVPRpAksPPwp8JuzY54GtWneGlEqp3xEIOr4G3BNSdznwPa3148F66cAK4Mng8Y8CE4DVWutXQp9EKfX3QMoA2i5inN1uZ82aNWzatIm1a9fKFgFCjJCL8RmoT/0zd8yZw45t22j+/pdI9LRdVa8iKYMZc+eNninRw0gp9Skgjc7RkTlKqa8G7z+rtT6otd6nlPoV8HmlVDKdK9y+G/iW1rpjSGeo6oac862eypHSr8BFKZVPYMzqHwjs9GjQWr/RxWOnlFIHCfSUhDoNfEAp1X7OgwR6VtpNDv7s6ppuQMYLxiiHw8G8efPYu3cvCxf2ZSkBIUSktVntTJ8wAYvFwvSZM3knMZXEuqsDlytJGSzIyxuBFg6B4V/H5R+BkpDy/OANAjkk7d1rnyDQy/Fh4CECewp9hsDy/OGGqu6w6W+Py/eAk8CvgD5lC6vAV+I8Av/IoX5GoAdmb7C8Efh5yPGzwZ8fAX7az3aKUcxqtZKVlcXp06dpaGggJUU634QYbtmuBkpLS5k9ezbJyclUp2Uzvq7iqnqeksl4vd5RkJirhz1w0VqP72M9D4FhmV57OIaqbneC6SXvAhYD6VydW6u11n/Tn2v2OTlXKbUSeC/wcOhQUB+8HygG/hD6oNa6FVgJzCGQubxaax06avwoUAE8opQ6rZT6lVLqI0qpURK6i8FwOBxcd911sr6LECMk1dNK6bN/xO/3o7Um5c73ciKr6KopJPYJ19C3jYTFaKOUSiWQq/oHAqkkHwY+FHJ7KHjrlz71uIQk5D6htd7W14srpWYS6C3ZTaCHxRAMgA51da7WukoptYjAltr3EPiFPwxopdTjBAKoxq7OFWOHz+fD7XZLoq7ok3HP7OzxePiiFY7lnbNBq8rMRNbwlXItvXwNDK0fnujr85kf7H6vWcHbYj53QnZnErE/7JguM89tOd9glF0hSbT1dWbd8Hal+MMSfz1+/ErxzvT5eCZNo+78Werq6sjOzub6lTdSO2cur//6/5i/4UnsPi8FG3aRtmvX6Pj7jOw6LmPFVwh0TnRlwNFsX4eKPkNgnG1VXy+slCoCXgLqgLuDXU79orW+SGBW0ieDU7JWEshofojANKqH+ntNMXo4HA6WLFnCvn37uO6660a6OUKMCQdW3s5t//gl4uLiaGhoID09HQj8Pebm5rLsvQ9y8a1XyK2v5sqVKxQVFWHpLaqLCcM/VDQKrCMQoLwMrA3e/08CCcd/A2wFHuvvRXt9NQW7er5BIK/FoZQar5QaH3xigMLQtVyC52QDrxBYEnhVH3e47JHW+rzW+jfAMuAU8N5gT5AYw+x2O1VVVaNgjQghYkP6yltJSkrC4XCQlpaGL2w5f601luAHfE1NDRkZGSPRTBEdxgd/ho64PKu1/gTwHQIzifv95t2XMDgdSCbQ01EacmufCr2dQNQEdAQ6GwhswLRWax3RtVaCM4oOAA4gK5LXFrHH6XSyaNEijh2TJX2EGA46mNMCYLPZjGEgj8fD2dOnSWlpAsDv94+S3pag6FyALpq1Z2TXAe2DpanBn1sJjJx8ob8X7csr6gqBjODw2x+Dxz8K/C2AUioBeBGYDtyptd7d3wa1U0otCc4dD388HVgK1ACVA72+GD3S09M5ffr0SDdDiDGh4ehhPJ6uR/611lx57i84vYHjVd/75uhaa0nryNzGjprgzzg6lzD5e6XUDAI5qxBYr61feh1q0Vq30MUKeUqpecG7G4K5KABPANcDvycwhPRg2LV+24+2PQh8UCn1DLCLwNYAxQQykccBn5T9jQQEhoumTZvGxYsXKSy8ankhIQYsbUJqx/1zp8zvSeEr44Yn6wK8Tx/vuP87NbXbY392msfi481v5VkWM1CwJ3S+dbeeNVfGjS8ylwcIXx23uanz2i0tZns9nrAkYb8ZdNia3bRs343jU2ayrdvtxuPx8OKvH2PmobcBsNos+JUaXYGL6K9zQC6B1JJdwB3A7cEbBHJezvT3opHOEWlfGOe9wVu4/gQuPyMQrNwCrCYwZFVHYH+Eh7XWzw28mWI0UUoxceJEXn/9dQlchBhCNTnZtNx8Y8fwj8fjobGxkVef+C3JO15jbsUF7L7OCM6v1OgZKtJaZhX13x7gOgKjMD8kkKwb/oL49/5edMCBi9b6mwRW0Q19bPxAr9fF9Q8T2CdBiF55PJ5uu6+FEIPXnJxE0nf/jQ8tXozf78fr9XLx4kXe/s6/MHvndhIzr57yHOdx09TURHpaGmq0BDCiPz5PYEq0S2vdppS6E/gsUESgN+anWutn+3tReSWJUUMWuRJiCClFWno6Vqu1oxclPz+fyR//O4596G84Nm4ibXYzeBlffZmNv/gpnq7G0WKRJOf2i9barbWu11q3Bcsvaq1Xaa2na63XDCRoAQlcxChhs9mw2+0SvAgxRBIbGtn/s5/jcrnw+/04HA7i4+NZsGABd3zq06z8xRNc+ty3uZDZubi5RWvmv/4cL/7pD6OjR1QCl4hSSt2ilHq1v+fJOihiVLDZbGRkZNDS0kJiYuJIN0eMEnHpcR33wzsNekvG7Q+32yw3NoYH4OY8hKRLTd1ey++pM8qeJvPiPaVphB9rbNSsqjqOy+Xi0KFD5NTW8uijj3LvvfeSnZ2NUgq73d6xD9GKW27h9PjxnPzqPzHLeyXQVjQTn/81O8dP4NrF142OVXRFr5RSSwlMqLmotX4r7Ni7CKyKP6BdciVwEaNGdnY2tbW1ErgIESEaeOutt2htbWXp0qUkJCRw66234vf7r9o0UWuNw+Fg0qRJNPzjV6j472+Q21IHQFZrA3W//D6XCr5PUVERNlsMfvRIcm6fKKWSCKyUuzTksX0EVt5PIjDreEn7IQaw9L8MFYlRIy0tjZqamt4rCiH6pKEwl3HjxrFy5UoSExNRwenNXe303D7t2eFwMH/hQi7f9zEaHfEdxyfXXKK0tJSWlpbYXena74/MbXT7JwLLoqiQ23zg+wQWp20PaNrnyW/o7xPEYNgrRNcSEhKora0d6WYIMWq471lDQUEBdru9X+c5HA5W3Xk32zKzOLpzK0prnJfPceuKFfzud7/j/e9//xC1WESBO0Pu1xMIUFKAD9AZrPiBJ4Hvaq0P9PcJJHARo4bf78flco10M8QYFr7IXKzLWjBnwOuw2O12rl+2HN+SpdhsNvbu3csrr7zC3XffHaNJ9DJU1EeTCAz//Exr/UkV6Ir7KfDx4OO7gA9qrU8M9AkkcBGjhsfjITW1c6XT2tpaLl68SFFREWlpaSPXMBGzEv97U+f9Z6aZxxJh2bnOPbIGE7SEJ/W+nGFeKzxZ9+LhztVy45zmtdIKE4xyQo5ZTqpu67jf1tb9Srleu43U7OxB5aPY7Xbsdjtut5vS0lKuvfZaEhISYjTHBfDHYsA17NqTDP8CoLXWSqk/EwhcAN6nte73armhJMdFjAper5e9e/cyZ84cAI4cOcI777zD+PHjOXLkCEeOHBnhFgoRWzxxTlJSUnqv2Au3282lS5dITk6msLAwNoMW0R/t0W9ryGMd0fJggxaQHhcxSmitaWhoID4+ntOnT6OU4tprr8XhcLB06VI2btzIzJkzR7qZQsQMpfWgkmg9Hg9aa7Zt28ahQ4f4+Mc/HvtToWWoqD/+rJRqH7vv6BdUSoUHLlprPak/F5bARYwKFRUVTJ06Fb/fz6FDh1i3bl1HQqHH4yE9PX2EWyhEbBls4ALQ3NxMdXU1c+bM6XImUkyR6dD9lRdWbh9nKwl5bEDToSVwETHP5XJx8OBBbrrpJp555hluvvlmI6HQarVK0q4YEk8nRCYZNzw/JiMK4mzl1/h8vt4rdsNut5OcnMz1119PRUUFPp8v9oMX0VdDuiW4BC4i5imlSEpKYs+ePdx+++04nWa2os/no6mp+5VGheiLvCJzSvChveaKtAnxcHfLwFbODdfWZpbDU01CJ/okZpjDL7Z4s51+j9lLYHV2Bg8pyebyv3WdOb/4rYPvcbHZbGRmZrJ3714uXLjAokWLyMzMxGazxeau0dLj0hcfHuonkMBFxDSfz8fhw4eZPXs2e/bs6fIbndVqJTk5mdbWVuLj47u4ihAiXEN2JtPzwnv7+89ms5Gens7ChQs5e/Ysu3btYtWqVcTFxfV+crSRwKVXWuvHh/o5YjDkFWOdy+UyNmybOHEib7/9NklJSV1u5Gaz2Zg3b57MLBKijxoy00n88j+QnZ096Gt5PB5qampwOp1MnTqVuLg4mVkkBkUCFxFztNacO3cOt9vd0ZuyePFi3nnnHUpLS7s8x263U1lZOcwtFSL2NGSmY/+Xr3Ljqtv6vWJuOJ/PR1VVFRMnTux4TGsdmwvQtSfnypL/I04CFxFzHA4H9fX1nDx5Eq/Xi9VqJTU1lfvvv5/Jkyd3eY5SKjbH1IUYRh6nA77xBW5cdVtEpi77/X727t3L9OnTjccGk/Q7ovw6MjcxKNJfJ2KOxWJh5syZ/PjHP2bKlClAZx5LdxwOB1lZWdTW1srU6CjU+pXbjXLZmxc77l8qNZNgw1eRDY9Hc3M7JzQkp5lvcSlF5msksSDJKFcfrjLKzZWda2g5knrvfQidZdTSah5LCEmvitbPrqqicay55eaIBC0ej4etW7dy3XXXdWzA6HK5SEtLi/31XMSIksBFxCSlFJMnT+btt99m6dKlfRoznz17Ni+//DJz586lpKSk1/pChLp8zsyfSk01Z3zW1w88Ggmf/hwafAEkJpnRmcUWMt3fbh6LzzQT0N0N5lIArubOmUQ1YXuStkxMishquT6fj1//+tfccccd5OTkdDx+8uRJJk+eHJu9nxoZ5okSMfjqESIwTh4fH09jY2Ofp2w6HA7WrFnDmTNnqK+v7/0EIcaQmoI8WHdblwnu/eXz+bj11lvZs2cPmzZtwuVyceXKFa5cuWLsJxZbJMclWkiPi4hpGRkZ/arvcDhYvnw5mzZtYu3atUPUKiFij/3/fZqH7rsXr9fbe+VeOBwOSkpKyMvLw2q18vzzz6OU4vbbbx90wq8Q0uMiYtqUKVPYuXMnbre7zzMVLBYLEydO5NSpU0PcOiFig18pLA4Hbrc7olOVnU4nSikWLlzIrbfeGvNBi9b+iNzE4EiPi4hpmZmZlJSU8PLLL5Ofn8/EiRNJTU1FKdXt8uJWq5UJEyawefNmcnJyIjKmPxaUt5qJq/HWOqOc6jrQcV83XDFPvnzZKPrPnI9o28TgVE0oZOWS64YkadZqtZKXlxebU6BDyV5FUUMCFxHzioqKKCoqorm5mQsXLvD2228za9Ysxo0b1+23R4fDwc0338zzzz/PnXfeKXuojDB/i5lX4WnqLId/liYmmomrzc3mB2JlZWe5ocEc9rAnmFN9Ko+Z2amXLpvXCl16vy1su6vwz7DeOircIb9iWliah93e89Yutnjz4qH/Pn6v2ZCaEzVGOTx5N7Qdra0arRT6wXeTlpbWYxsGI9Z7WkR0kaEiMWokJiYybdo0br/9dg4fPtzrNzy73c4tt9zCli1bhqmFQkSfqukTuOW+eyLS2+L1eqmrq6O+vh6v1xuRRN+oIsm5UUECFxGT2teF6M7ixYs5frz3De/i4+MpLCzsqHv58mVOnToVkQRFIaKd32LB+amHIra2kc/n4+WXX+b8+fO8+uqrvPnmmzQ2NuJ2u3s/WYg+ksBFxCSbzYbdbu+2VyUrK4uzZ8/2+oZpt9sZP34858+f5+TJk5w/f56UlBSee+650fdtUYgwFcvn8+73PIDFYuHMmTODDtidTidr1qzh5MmTpKamMn36dHbv3s2xY8dGR/AiK+dGBQlcREyyWq1kZ2fT0NDQbZ0JEyZw6NAh3G43Lper23oOh4Mbb7yR/fv3M2/ePHJycrjpppv61GMjRKyqzs5m0T9/CYA//OEP7HjgY+zds2fQwUtSUhLr169nwYIFXLp0iaSkJOrq6kZPcq4MFY04Sc4VMSs3N5eqqqpuF7SaOXMmNTU1vP3227S2tnLjjTd2O45vt9u55557OsoJCQmcPn2akydPkpKSgs1mY+XKlUPye8SKM3VmguXUDKdZwR5SbjAX+POfvWCUXXsqjHLdqTqjXFXd+SEX/kXd5zM/AMOPNzR23k9KNOueP9FGT1pawsqtXdcDiIszy+GfR+GLwyYldn8sIcEs28JyxVvrzF/SfK7wgMAMPBqrzIbV1gbqVz64jltnz+bAgQPY//m/Sa6o5uKe/cyZO3dQU6JDz507dy5lZWUcPnyYo0ePMmfOnNhcNVdEFXkFiZiVmZnJpUuXeqyTkZHBDTfcwJIlS3jhhRdobW3tdoM3q9XaMbtIKUV+fj7r16/nlltuobW1h08wIWKQ2raLX3/tWxx/6NMkV1TTlJ1O9vSpEQ0s7HY7+fn5LFy4kJMnT8bu5ortpMclKkiPi4hZXq+X5ubmPtVNSUlh/fr1vPLKK1x33XWkpqb2OEVTa83ixYsBKC8vp7i4OCJtFiJa5Ow5jPftwwC0pKcQ/29f5Lpl10d8LRe73U52djbvete7YnvZAdmrKGpIj4uIWVarFZvN1ue9iux2O7fffjtHjx6ltLS0x2TB9jdvt9vN22+/zfTp0yPSZiGiUdM9q7l13boh27XZYrFgs9l6nQ0oRF9I4CJilt1uZ9q0aVy4cKH3ykFKKW644Qa8Xi8HDx7sMXjx+/2cPXuW+fPnyxuuGNXiCvJjfxhnyElybrSQoSIR0/Lz83n99dcpKSnp13kzZsxg27ZteL3ebr9lejweTp48ybp16yLR1Jh3fb6ZBF3aaCbnZlo7k0J1rbkirTcs+bbmsLl9wLkjTUY5dAXbwUxyCU+gTUk2y5VVYSvlhk0+C033CM9XTU8zg9nwpOHwVJHQ4PeqZNywa4ef64g3h1jiM+PpjqfVnMZfXWX+UqGfm+3/ttbQzGHRPQk6ooL0uIiY5vP5+rXBYqgpU6Zw8eLFHutkZmYOtGlCxAyb0yG9iiJmSOAiYprT6WTKlClUVFT0XjlMfHw8LeHzX0NUVVWRn58/mOYJERNayytjf52VoSbruEQNCVxEzBs/fvyAFotzOp1XTXMOzXkpKysjOzt70O0TItq1ll6Q9VX6QlbOjQryShUxz+fz9XladCibzUZra6vxTbO2thaPx4PH46GyspKE8GQEIUYhz5nzg1p0TojhJK9UMSoMZHxeKdURpLQn6J44cYKmpkCiqAwTDUJYlqu3zEy+rSw1y2XlVye23ufquhftd2pqj08d+vmbmGi+LtILzKTW06XmUGFPvfgZYfsQ5oR1xjU0ms8Vnggcn9aZBO5MNRObvS392xfLltD5S9qcZuKuq8H8tw/PPW8MWVm4/d8qrqGB1tZWEhMlSbdbso5L1JDARcQ8q9U64G5uv99/1Towq1evlm5zMaY4m1poamqSwKVHWgKXKCHvziLmaa1JSUkZ0LkWiyW2V/MUIhK0luRcETOkx0XEvIaGBjIyMgZ07kADHiFGFaVkOnRfSI9LVJAeFxHzUlNTuXLlyoDOLSgooLy8PMItEiK2aAlceqdB+3REbmJwpMdFxDy/39+RUNtfBQUFbNmyhfz8fJlV0U9bzlQb5aSp13Xcz0rfb1a2mh+KbW3m4R52XrjK+7SZtBuerBueFBtq71tmMm54O8KFviRSU83foSlsIlt4WlT4l/PUks7evfgCs5E6LDm3rrTevJbHXI7f2+INOWY+kbL0/G/dVZuVRQIXETvknVrEvPYN3Px+f7+Taq1WK21tbfh8PpRSuFwueQMXY44n3klc+P4I4mqyBktUkKEiEfPaN1ssLS0d0PlpaWkdiYler1cCFzFqdfex25aZJjOKeqM1+CJ06wOlVLFS6pdKqVKlVKtS6oxS6n+VUkVh9WxKqa8F67UppY4ppT6lungj60/daCaBixgVcnNzOXz48IBmRixYsIBNmzYNqMdGiFhxpbiAuh/8G2XXzzcCGJ/NSsLNgR3TRXRQSmUCbwPvAp4APg08C3wQ2K6UCt3x9H+AbwMbgU8BB4EfA1/r4tL9qRu1ZKhIjAoOh4OFCxdy4sQJpk7teYGycHFxcSxZsoRDhw5x4403Dk0DRb+F5q6E57WI/mssKuCD97yLsqVL2PSej5J75BR1hbnEf/Wz3LV2bbe7pIsADejhGyp6AMgD7tJaP9v+oFKqFPghsAp4Uik1F/go8AOt9eeD1R5VSj0JfFkp9QutdVnw3D7XjXYSuIhRIycnh5deeomJEydit9v7dW5WVhZZWVlD1LLR6aG5BUb5V/s6d9p+35y7jGNx1581ykXHa41yfX2jUa4xD/coPKjZlNMZ8FRUmB80LebWVITnY4d3uM1b0Pk6SipIMo6dfLvOKLe2hn+omb3v1cdqOu7nOM0n9rnM3o7wBFtXvZm97GntrG+1m422xZvXttqgLjsLy7UL8Hg8ZGRkkPSR9+D+wnexfu7j3LZ+vQQtfaHp8zBPBLRncl8Oe7y93J4a/kDw54/C6v0IuA+4m0AvS3/rRjXpFxejhsPhYNWqVbzwwgvS7S1E0PmpU5j828d44FOfxGq14nQ6ufHO9Zz64kdZfe+7JGiJTq8Ff/5YKXW9UqpAKXUb8K/ADuCV4PFFQIXW+lzY+TsBP7Aw5LH+1I1qEriIUcXpdLJq1So2bNgw0k0RIiok1Tfg8wWmU7fncCUmJvLwww/j6GfP5JimAZ8/MrfenkrrncAngWnAW8BFAsHKceA2rXX7N7NxwKUuzncD1UBot2h/6kY1CVzEqGKxWLDb7VgsFlnCXAggs7yCvd/+V5qbmzv25XI4HIFlBCRwGSlZSqndIbe/7aLOZQK9K58H7iKQQLsCeE4p1b5baDzg6uJcgLbgcQZQN6pJjosYVfx+P1VVVeTl5cm05lGktx2hQ4Xv2ixgwtu7eP4HP+LdX/mSDA0NmI5kcm6V1npRdweVUvcAfwTmaa2PBB9+Vim1F3gB+ATwA6AVcHZ9FeKCx9v1p25Uk8BFjCper5ddu3Zxxx13jHRTxpyPzC8MKRUaxy6v+6FRzrn4N0Z5uuu0UT66r/M9tD+JugDWkFV6bTbzgyYt1awbHx+WBOsy65863rmi7WTM1ZnDE3vDV/+trzevlZLcOUTgbujui2/X/F5zeCE0hcvmNH+HhJwEo5zhCgwTJb72DC9dM5m199/f5+DF7XZjsVjw+/1YrdaxvSHp8CbnfgY4GRK0tHsJaCHQ8/IDAr0ys8NPVko5gEzM5N7+1I1qEriIUcPtdvPaa69x8803S2+LEGHi3C5yHv8Zh6ZOZc7cub3OvPN6vVRWVrJjxw6SkpKw2+3MmTOHlJQU6bUZenndPK4IpHi0/+ftAW5TShVrrc+H1Ls2WG9PyGP9qRvVJMdFjBpWq5X4+Hi2bt3Kpk2baOttIxohxpis+hpOPvdiR65LV7xeL263m127dnHy5EnuueceVq9ezcqVK9m2bdvYXqTRryNz690xYIpS6rqwx99NYFhnd7D8p+DPh8PqPQy4gadDHutP3agmPS5i1LBardxwww14PIHu/ddee421a9eOcKuEiC7efQe6HfLxeDxUV1ezd+9eli1bRmpq59ia1WqluLiYxsZG0tPTh6u50SO4O/Qw+S6wFtiolPopcAaYA/wtUAb8FEBrvU8p9Svg80qpZAJTm1cRCHC+pbXuGP7pT91oJ4GLGFUsFgtOpxOfz0daWhrNzc2yB4sQIVRwQ9Hu7Nu3j7Vr13ZZJzMzc+wGLsNIa71NKbUI+DrwXiCfwJTl3wNf01pfCan+CeA88GHgIeAsgRyZH3dx6f7UjVoSuIhRyWq1MmPGDE6fPs2cOXNGujlCRI9ugha/38/+/fu5/vrruw1s8vPz2bx5M/n5+f1enTr2aehhiC3iz6b1QQIr2vZWzwN8K3iLWN1oJoGLGLXi4+M5cOAAFRUVFBQUMGPGjJFu0pg1LiHNKOsHzVlfOVV/MsquhrOdhdM+41j4LKPwlIu0kKe6ZlGKceziO+bWAuGfQz7zqair77x/7IjHOBaen9oWNlGoqdksx4eskmE5VWccS8wzewVtCeZbc/qkNKNsjesc6rHYzGEfZ645q8jnMn+pnNY6rly5Qn5+vvG41+ultrbWGB4KZ7PZiIuLo7GxkZSUFGzhU6tGs+GdVSR6MIazrMRoZ7PZePDBB7nllls6ZkgIMdbllpez+Te/7nJbDK11j4m7AMuXL+eJJ57A5XJ1rMgrxHAaQ+GyGGvaZz8opZg6dSr79+8nO1tWJxNjm0KT/MYWXJ/8lNFj4nA4WLBgAS+88AJ5eXkUFxezb9++wDnBoSOtNVprVq9ejd1uH3Prugzj7tCiBxK4iDFBKUVtbT9XMhNilLL4/V1uiZGamsqaNWu4ePEir732Gg888AAej6cjcPEHz7PZbGMvx0WGiqKGBC5izFBKobWWxelGwHc2nTDKX5o7Qg0RAGiluvw7cDgcHdtmzJ07l4sXL1JcXDwCLRSiexK4iDHB4XAwYcIEKisrycnJGenmjHkqxfw/iPu4ud5Ofs1THfdbm83lJfxh3fXhqRqhS+9nLMw1jmXdWGSUzz5x1CjX1prXDr1W+JL+/RV67bg481hCjplXklyUbJQtSWYmsK0kJOk4rBdAe8xrpTWaDfe2ePDnpREX3oggt9tNZWUlixYt4uWXXyY3Nxens7stbsYSLT0uUUKSc8WYUVxczMmTJ0e6GUKMOL30xi6TcwHsdjt2ux2llCztH0oHclwicRODI4GLGDMsFguJiYkdCYdCjFX21O4XkLty5Qrjx48HYNmyZTz33HO4B9vdJEQEyVCRGDMcDgczZ85k06ZNI90UIUaUu6a6y+Rct9vN9u3bGT9+PK+88gper5fc3FwcDgc+n2/MzSK6im/4FqAT3ZPARYw5Ml4/9F45VzPSTRDdKM0uYP6a23GEzQryer1cunSJxMREpk6dypw5c/B6vR3DRWM9aNFapkNHCwlcxJjh8XjYtWsXCxcuHOmmjDnXTzOTcV+vzOHGwrSOsgqbuBK/bGfH/cl2c0Q744QZFF083mKUm5s7P1ze/ImZ0zRhkvmW19RsfhAlJpozberqu0+oTQtbYLa8gh6FrqRbUWE+r8Vi/g72sGTclBJzBeCekkT99eYSvu6mwIq/LmXlYOZExn31s+Tn52MJC0S8Xi8Wi4WbbrqpI1gZUyvjipghr0oxZni9XlpaWnpc0lyI0ehEYi4N6z7MshtW4HA4ulyDRSlFSUnJCLQuVsisomghgYsYM1paWmhoaKC6uprMzMyRbo4Qw+JkQjaZn/kO84uLe5wlJEOovdCADBVFBZlVJMaMlJQU7rjjDvbt20dVVdVIN0fEGA20ZaXiSk/utW40KZ+8iOJeghYhYon0uIgxo717fOXKlbzyyiusW7duhFs0tl1qbui4nxl32jgWbW9M7vQkPJ97N0vefTf19fW88+EvkFZ9ufcTR5hXK5JnL+rYt0sMjpahoqgQbe8PQgw5i8VCQUGBDBkNoY17LhrlnPR4o/y+xeaKtuG87/08cdb5AFiv/Q/jWNz+I0Y59bXzRrn5Smei65Xj9caxpgZzN+OMcebwyPRr883n+stZ3A4Hdf/8ZdbcdSd2ux2Px8OJe2/E/n+/N+rmZJsfald62Iy8odEsuz3muXHOOqNstZmBR0LIkIUl1fwdmo51Ji/X2+PJv69QkmwjQYaKooaE4WLMsVqtzJw5kwMHDox0U0QMUGiSMzoXbLNaLCz54PupzY/+rSM0Xe9JJEQsk8BFjEk+n4/W1taRboaIAXa3hys//inV1dV4PR4sVivFxcW0vO9eov37twWNx+MZ6WaMEjqwAF0kbmJQJHARY1JcXBxFRUXU1taOdFNEDCg+eJhNn/siVyorOxZlu+X976N8yoRBX9uTmkDVLQvxffezNH32fbTkZUSgxQEaeOfAPlwuV691hYgVEriIMeuaa67hxIkTI90MEQMUUPTGDjZ96Ru0tLSgtSY1NRXn33wA3wBWlNUKmq4poO2f3sPMTb/gob8+xh1/+2He++0vU/Sr79CQPPiZS9XOJMo+9g3mXXvdoK8lkE0Wo4hkbIkxSylFdXX1SDdjVPrPe+b084z53R5RU/+f+cD+DxnF5M+uNsuAmv5FALIefa9xzDZ7onmtoiKjqE+dMspzkzoXajv6+5PM3LWV5370I979xS/icDi49e67+N3jT5P55gHiwpZBSUq8+ne5s+k4VVVVbN26lZsWLyY7OxuLxWIsp3/98uX8aeENjH/mxY7HnOcajOu4Gjo3PUzIMROfW6sDQ6CHl67kvhUr8Hq9kpwbKTKrKCpIj4sYs6xWK8nJybS1tY10U0SMsGhNyTNP8/a2bbjdbux2OxM+/1G8CX1bvK2yspI9e/awbt068vPzsdvtV+0B5PP5SJw7a1DtrItPYvI9D+Dz+TqCFrfbTW1trez0LGKeBC5izLLZbCxYsICDBw+OdFNEDElsbubir36JUgqHw8GKm2+k9q4V6D7M3tmzZw8333wzdru929k+NpuNhMICyiaNpyEzfUAJwGeKpzJr1qyOoMjj8bBp0ya2b99ObW1tlztDi17IUFHUkMBFjGk2m40rV67IG7nol/zjxzl9OrBonsVi4T0/+FeqPvchGgpzeg00fD5fj8dtNhs333ILq57+E1P+9BvO3rO+38GL9dpleL3ejnJdXR15eXmsXbuWvXv3dvS6yOu+f7RPR+QmBkcCFzGmOZ1Opk+fzsWLF3uvLIbM3sqmq24dGv9i3qJASn09B197k7a2Nmw2G6mpqbz/y//E0hd+S80XP4YnIa7L8xqPn+lTsGCz2UhJSWHixIms+coXOVU8uX8NDOnM8fv9pKam4vV6efp7/0L9icM4nU7q6+vZtm2bEeAIEQskY0uMecXFxbzyyisUhSVqiuhkWb3mqsdU2nu7qAn2j5qr2/76oLlM//rMsHyPTMhwju8oNi2+1HF/RuKnjapHNrwCH/lAR9lms1FSUsIDn3+YvyQlkvi1/8YSsmZHmwv8Pl+/FoSzWCykpKRgf+jjnHjkRzRl5+C49jpsycl433iVqUf24WoITHXWQIszHndeLm3nz3Qs82+xWHA4HCxcuJC5c+ey4eufZsM3Pk1Saz1Wn5fmWbNkx/Q+0FqGeaKFBC5izLNYLJSUlFBeXk5eXt5IN0fEiHGHj7H55Q3ccvtaYwNDh8PBsjvXse3HvyWpzNzMs/bVHVg//Yl+PY/D4WDlrbfStGQJycnJeL1elFL477ufDU/+Cc/+XTgmTiF+/ETyiksozMnh/LZt/Pa3v2XNmjXk5ubi9/tRSnH48GFKPLVMag5sC1BtjefSpUskJibKzKM+8MswT1SQV6oY86xWK1OnTpWNF0W/2N0eKn/6KE0rbiAjw1w0Licnh9acNJLKqmjOy6RpyWxISSX+pW2cOXOGqVOn9u+57HbS0wPbDoQGGGsfeA/W970fr9eL1Wrt2Ej0jjvu4Nmf/JCL/++j7MsrxB+fhMXVwuTKUgoSO4eGMn2tnC4vJzU1lfz8fNmMUcQECVyEINDrMmnSJM6fP09xcfFIN0fEiPxjp3jtyT+z/kMfJC6uM6/FarXiXTSLimsmsPz/fZopU6Zw+vRp3n5mC0feepuSkhKj/kA5nYFp2OG9JT6fDzt+ChurmKjNjSZJNJ/3+uuvZ8uWLWRnZxs9RyJMcFaRGHkjHrgopawAWuueU+2FGEJWq5WJEyfy0ksvMW7cOOk2jwK66tcAKGcXK7mF193zDaOsFn4r8PiFR8JqvqvXa/n0tpBSSY91ldZ4Ss9flbeitea2z/wd+fn52Gw2tNZcuHABR1UDlX/dgHr/A722YzBsNhtcONOnujt++01WPvhNCVp6oQHtl32GosGA3p2VUjcAbwSLRVrri8HHlwH/RGAZzBygCTgCfE9r/XwX1/ku8DDgU0r9UGv91bDjGcAXgPXAeMAFXALeAn6otT42kPYL0RWHw8Ett9zCtm3bWLFixUg3Z0zJTzSTZPNaOt8udM3VM75UyWc6j4cFLT15sOZ75nWqrzHK/nxzxd/Mxo0d991u80OrcLITt82O5+YVVwUuDoeD4uLijqGXtrY2Kvcf5oHWozQ3N7P9rbdYdsMNHUM7kXbp0iUmNV8iLt1J4uxs45h9SmZnwe+nzTVRdpAWMaXfgYtSygY8AjQD4V+FriEwEe8XQDmQBNwDPKeU+qTW+qch13kP8D6gPVPtO0qpA1rrJ4PH04DdQDbwOHAo+HwzgPuBHYAELiKinE4n8fHxVFRUkJubO9LNEVHuZMlUbrnhhi57K0LzRex2O81v7wUgMTGR0l/8huuXL+/z8/h8Po4cOYJSioKCApKSknrsIUlNTWX/4vWUH97GUt2EtYe4xNFyGZfLhcViGbJAalTQsgZLtBhIJtangTzg0fADWuv/01rfpbX+F631o1rrHwI3AQcI9MSEWk6gJ+ZxrfXjwH8CoV9zPwpMAO7VWn9Ka/2/Wuvva60/CuQDV/XgCDFYdrudefPmsXXrVs6cOcOFCxcoKyvjypUr1NXVyYJdooPXYsG6/r4+JbR6vV6sFZ0zjLK37mX71rfweDx9ei6tNSdPnqS4uJhXX331qm0CwqWkpLDmwQ8z5XP/wjlfzxs23pB2ga2/+zfZ+qIPZOXc6NCvwEUplQ98E/gyUNeXc4K5KxeBtLBDp4EPKKXmK6XmAw8CoVv1tq+49EbYeWit3VrrqvDHhYgEu93OHXfcQUJCAhaLBbfbTWNjIxcvXuTZZ5+lpaVlpJsoosDl9FwWLVvekSDbE6vVinvG5I7F3hytbZz5xa/7HAh7vV4SEhJITU29agZTV9qnP1dUVJCgel5gzqqgJK5RFqITMaO/PS7fA04Cv+qpklIqWSmVpZSaopT6R2ANsDGs2s+AVmBv8FYD/Dzk+Nngz4/0s41CDJrD4SAvL4+CggJKSkqYNGkSM2bMYO3atWzcuFF6XgRVk6Z3TFHui/f+x7+wffv2jnL21r3seKtvvS5WqxWr1Yq/j8mhNpuNP//5zzh++c/kWVt7rHux2cERzwSSkpL6dO0xS8uS/9GizzkuSqmVwHuB5Vprfy/JXP8H3Bu87wOepjOXBQCtdWvwmjMJJGwf0eanwaPAZ4FHlFL/AGwBtgIvaq3L+9puISLFYrFgs9mYO3cuZ8+eZcKECSPdpJj13v/aYpS//3ez+3W+/68f7rivSq6evu5/7bOBY1Om9L9xfRQ3d6Gx+3JPbDZbIM+lubnjMXtrG5ePHMN//fW9nm+xWHA6nfj9frxeL36/v8fhora2NlIsmsK2WgJpglfTGrZU55N//fu5Z8oUWcOlD2SYJzr0KXAJSch9QmtjrmB3vkWgR2UccFfweRII9Kp00Fr7CSTdXkVrXaWUWgR8iUCC74eDN62Uehx4WGvd2Jf2CxEpFouFcePG8fLLLzN+/HiZjREhhyrNcl5JyDTk0pPGMX3BnGXk3XfCKFsndiZVe54xAyRrboJRtsw3e0wsh18zyp4tnTuH1+wsM681LanbD3uPx4PW2kigtdlsHcFGfLyisqSQFffc1eNQk8/nQ2tNZWUlbW1taK1RSvXa4xcXF8e85SspO/kyM2aaOS4qNweA8noonP4+Jk6cKNP/RUzpa4j9GQILGnyhL5W11oe01pu01r/WWt8LtACvKaV6Hww2r3NRa/1JrXV+8Pk/COwBHgJ+3J9rCREpdrudOXPmUFpaOtJNEcOoxeYwdmm2OBxdBq4ej4fHH3+cI0eO4A0ZBqqpqaGgoAAAv8WC98H7ycrK6vK5tN+Px+PhyJEj/PUbX+PYh9/PbbfdxuHDh7nuuuv6tOZKTk4Op8bfiNvXdXB9vMIuQUs/aA1+v47ITQxOr4GLUioV+AaBvBaHUmq8Umo8ncm2hUqpwl4u8wdgCv+/vTuPj+o6Dz7+O7PcmZGQEEgIBIhFLAIhdrE7IIQFEthgu47tN3Zct87HaVM3aZM3bdokTZrmbdLEbZq0SZo0SdPESZO4SbywY1YjjNnBgNnFIiSBJIR2zZ3lvH+MNGikQQtopBnxfD8ffaSZOffqSLqaeeac5zwndNVQj2itr2qtfw4sAS4A/6dlJEiIPqWUYuTIkZw8eVJyXR4Q75tJVHzmmxQtWYfHEhg18TY1hs058fv9OK+Wcexb/4Gv5XHTNNnz1gamTJmCz+ejuCCXNR997q4BiNfn4zevfIOmP3+JmZteZ+TNMiwWS3DUxdeNRFrDMCj48LNsbVzInpJ49l6P43r9nafMJsc4Scjtkd7Jb5Ecl/vXnRGXIUACgUJxxW0+WitAvUsg96Qzrjbnui9aa5PA8moDCP92RYgIs9vtzJw5kzNnpJTQQFfudeJZ9AI5OTk8+vKneW/hGrzKgre+PmzgYrPZmP4Ha3EdPs0bP3sVj8fD3t17qP+v1zh06BDr16/nyX/+RqdTROXl5YzdvpmkutBy/QsWLGD//v3U1NZ2K+hwOp0UrHuaZS9+g8V/+HWOsCj4mMXqkKlOEZO6M2Jxk/B1sp8BniZQb+UagFIqVWt9s20jpZQFeIlAAu6h7nZMKbWQQMJuXbv7hwCLCOTLVIQ7VohIU0qRlpbGiRMnOHnyJNnZ2f3dJREhp/wjKFi4CKUUVouFR17+NG+hoCl83ROr1UpWVhYnnn8M5ys/4pUrVxgzJp3sf/0KWVlZzJw5s8sl1Ha7HV+75NtbL3wIgHnAnu0zWfzZv2Po0KFdTvUYhoHH46GkpIThjceB1r2K/DJi2BOyV1HU6DJw0Vo3ElgVFEIpNavlyy2tJf+BzUqpMgJVbUsJFKp7BsgG/kVr3b3NMwKeA55XSr0BHCSwfcAY4A8JJP3+mexvJPqTYRjMmDGD0tJSNmzYwLBhw5g7d26XxcEE3LhUHXLbaQsd/G3y3qlV4pqUFfKYmjAl5LatJuS9Evra9eDXlqTQAMGSFVrin8TU0HNnp4JafudcW59glLeCyspKhg0bhrJYcDocPPap/8vBgwfD/WhAIPB4/M8+zr7Z0/nLlu0jDMPo9sqdlJQUDmZmMMkZqBnkHHLn51BAvrrAvv/4KtM+9jlSUlI6zXkxTZPi4mJuX9nEghVtkp5LZbSlpyRwiQ69nSPyI+ApAtV1hxAINo4BH9Fa/08Pz/UfLcevAFa1nO82geTcT2qt3+qdLgtx7wzDID09nbS0NJqamigqKpJ9jgaYCdZatu/aTsG6JzAMA2WxYBgGOTk5nY52uFwucnNz72nzQp/Ph+VDq/BfOYWFji+WClhcf47j3/sCF/OfZ+HiJfh8PgzDCE5fud1umpqaOLB/LxlDLzN/WrtRHnnfJ2LUPQcuWusvE6ii2/a+7wHfC9f+Hs5/Evhcb5xLiEhqLQ4mdTAGJpvSzCp+k13b4snNXxUMRLqa7mm9Lu6Fw+FgYd7DvL/xV0y+dT1sGwXMMkupP/Ej3j65E506iYQRY/F6TNxnd5CqbhNn8ZG/agY2W2hfPV6N1yfXa0/olgJ0ov/JqhwheoFpmhQVFTF37tz+7oqIgOHWZqp/9q/stdlZkru8W2X+71dCQgKNa5+j+lffJY27V9cdYvPwMGeg8gzeCoUl0cA2+M7jFlvHKaHtRzQrCh7F6XR2eEyIaCeBixC9wGaz4fF4SExM7O+uiPtU67lT0K5teJLSWEv5b3+Kf+mykPamad7X6MrdGIZB/iNrOZg6guItrzGx/BypZuc1N21K052BPy9JUr+lxzS6m1suiMiSK1eIXtDc3BxSTKx1tYYsNw1v7LTQMvSXq0I3rkx03NlOYcqQ0GW/TkL3V9W+dqMRzReDXyp7u1fxoaEVFHT5eVTaS8HbbYMWAPuEOxUcrO+V445PwDRNXK5AhQePx8O5c+dITU0lNTU00bc32O125s2fj1qwgIsXL7L33d0kn9/HzIl+nJbAi6i1MC/0II875KZKTg+5rZ0JZM++zblz55g8eTJ2u73X+z0gyVRR1JDARYheYLfbuXTpEsnJyRQXF9Pc3IzD4WD58uVdHyyiXrNP8c6ix1n5xy8FgxYIFJs7/b1/o2TVavIKC+8pEbcrrYFFZmYm48aNw+P5CIcPHqDh7DsssBfT9V7RHY0bl8SO7UcYP368BC4i5kh2lhC9wGKxsHbtWiwWCwsXLqSgoIC6OtlKqz+YPthaNoqdPML6prncNu/v/VmzT7GlYQ5rP/FJBg8eHBKcnDhxgomH36PxrddpamrqVkXb++FwOBg0aBBLPrSUgpe+yEE9ueuD7uJDS0exr2hPt3anFgHar3vlQ9wfCVyE6AVWqzW4NNrpdOL1ekPemYu+4Tb9bKzKZtlTf8GyvIdZ89wnOFl3f3lH795KZfXTL3ZIyDVNk5s3b+Jwu8k8tJ+3/upzVN261Sdl9LXWnDlzhnGey/d8DrvdirepGp9PlkV3h+xVFD0kcBEiAvx+f9hy8CKytr3nYfWTL+ByubDZbPh8PhpdaZ0eow9/KfjRntunsE5YETZXyWq18tBDD/H+Mx/F7XSSsXUn2/7tuxH/u3s8Hk6cOEHtrh8wMa6x07Zer2ZXUSXv7K/qUCW3ocHEmTBcCiaKmCM5LkJEgN/vx+v1YppmRPIeYl1TQ+j0xJbtF0JuD3t8WvDrEfGhyaXD646G3NYXzgJwq1ExbtKzHX/fien4K48Evt4Zuq2aJXN8aFNv6N5TV1QikzKnhP0bWq1WEhMTeeZzf8trfg/Zv/kfrhafjVg9H4/Hw61bt3h3w/8yf/R1RsyJA+JQcUmhDW2Bvvp8fjZsqeKxJ16gvLyc4pJ3yBg/GKUDgdX+d8tYuvQZyXHpAUnOjQ4y4iJEBDidTvLy8ti0aRM1NTVdHyDu29kKG5MnTw4JHJRSWB0dp+xqPVbKmrsOKGvdFhISEkLu8/v9mKaJz+fj1q1bnD9/Ht0yPRRfXx/RqaLDhw9T2LCbEYldv4AeOlbF8hVPoLUmOTmZmxWhozONDRbq6upobm7GNM1IdXng0L2T3yI5LvdPRlyEiBCXy8WaNWvYs2cPeXl5XR8g7kuDP75box1aw0ZLLmMyMohrfJshcXdv67DpDi/qpmmy4cf/jv3gQZLqb5PUUEtOcyNuIK6pkfr6+ogUdrPb7QwZMgQPFroaI9FaU1UzFJfLhVIKi8VCXX1oLkt+fipXru7gwoUxWCwW5syZI9NGIibIiIsQEeTz+eTdbB/x2wZ1K9G0wjTImb+Q6dOnc72687aJhp/q6tBGVqsVrTXjy6+QXlVOQnMjrRkwKbcqOHvmTMRWF2VlZfGObzxdberc0OhlxIj04BSXx+PB9CSE5Lm4XHamZKaQk5NDfX19RPo70Gif7pWPnlBKjVRK/VApVaKUcrd8/q1SKrFNG5tS6otKqWKlVLNS6oxS6mUVJjmrJ22jlQQuQkSQ1hqfz0dTU1N/d2XA09bwQyeaO8/HWsMRM52MjAwSEhK42dD52EWyy0fJtashAZHNZmXdx/+C+i99i9Mf/2uOr36Go1nzKE8dgeHxcH39G/i7iizuUXx8PPOe+ww7iwd1Grw0NHoZPPhO3X/DMJg1azEn3q/s0Pb8+fNMmDBBRlu6ovt+ObRSagpwFMgDfgD8KfDvBIo6t73gvw98BdgGvAycAP4N+GKY0/akbVSSqSIhIsjpdJKfn8+mTZtYsmRJSHXdcHw+HxaLZcBX3P3VZ3M7ffzrO+5UsL3dGLqc+YnMGSG3jYpyAFSp2eH3ZrNZycyeyf4zbzPfUsp7VUOY+8cfx2Kx4PV6cbsm428+TPCwhtuhxw9PoelSGR6PJ/jCrpQFuz1Qr8ecMwetNXa7neLiYqqqqpiTnExVVRUjRozo5m+j+2w2G0lJSUxf8TH27H+VZdMVvvjQ72NtuklSgoVLH9wkIyMjOFU0bNgwKiqmcfL902RnBa7D0tJ6KivtzJ8/v9f7Ku5PywjIq0AJsExrHXZYTCk1E/gY8C2t9adb7v6RUuo14G+VUv+ptS7radtoJiMuQkSYYRisXr2a9957j8rKju94W5mmyfr163nzzTe5detWH/ZwYLD6ajvsv2OxWBk1ahSTP/oFfpWQT8bzn2fo0KFYrVYcDgfDRo3H3cWsjtUsD7uvj1IKh8OB0+nEarUyduxYUlNTaW5u5ty5cxFL0rXZbAwdOpQp8z/C3lMd3727TR/vn67m4sWLIcXlDMMgKysLLDM43jLycuJkPQsWLJCVb92g6fMRlzxgLvAlrXW9UsqllAo3RPh0y+dvt7v/2wRGZh67x7ZRSwIXIfqA3W5n5cqVHDhwgNra2rBtrly5wvLly1m3bh1FRUWdBjmio9QED9XVHQuq2e12UlKS+cgLL5KWlhZ8kTZNk7LLH+DoYtx5kLWhW99fKYXP5yMhIYG6urqI1nMJ/EwpZMx+mnf3XQ3eb5o+Nm6vZdrsZ3n66ac7BCSGYTB58mQcrjkcPnYTZYnvUN9F3IXu8xyXVS2fG5RS+4FGoFkptUMpNa1Nuxzghtb6SrvjDwB+AsHPvbSNWhK4CNFH7HY7+fn5HD9+POzjSinc7sAGeQUFBZw8eZLTp0/3ZRdj2vQxive2/Yy6uroOCdFWqw273R6c7vF6vZw7d46Z8R/Q1aycy06XZfE9Hg/Nzc1cunSJsWPHMnLkyODfMlLsdjvDhw8nffQijh0rBWDXzisUrH4Gl8t11/oshmGQkZFB4uAFXL9+jWvXrknZ/+jUup/DbwhMF30Y+DQwA9ijlBrV8vhI4Hr7g7XWJlAFjGpzd0/aRi0JXIToQ63VXMOx2+3BFxC73c7ixYvx+/1cudL+zZEIRylFQXYTp3/2Bfbv29fpi7HWmsvbfkN6UtejIhYLnY6emKbJ2bNn2fX3nw0uex8/fnyf/N3sdjtpaWlUVcRTWlpD6rCp3SooZxgG48aNo3D10xw6dChiRfMGlt4p999S8j9FKXWozcdLYb7hoJbPR7XWT2qt/1dr/W1gHTAU+EzL4y7gblFyc8vj3EPbqCXJuUL0IZ/Pd9d8gsGDB1NSUsLIkSOBwItLdnY2u3btYuzYsX3ZzX73ubw7mwf+/kJFyGPvV2WG3M4ZXxL82gosrK6momEjG35bzKp1z2AYdqzW0Ke6hoYGquKGs/mCBosNGqoY777ElISrIe1odqPNzpOlbTYbJSUlPPLN7wdHWQYPHszBgwfJzMyMeGVawzBYunQ5v//973nyyYfD5uOE0zpis27dOpku6gYN9OLsX6XWOqeLNq1LEX8R0g+ti5RSl4FlbdqFbqR1h7PNeXraNmpJ4CJEH1JKdZjG8Pv9FBcXU1JSQlxc6JJerbW8qNyDYfF+Cpyn2PHf/8SstS+RkpISEjDGxcXx9B++GLxttVrZt6+IceY7OI3QIMWv6TRw8Xg8uN1uNm3aRGJiIn6/n9raWpKTk/vsb2e1Wnnqqad6vKTZYrFEpFie6BWlLZ/Lwzx2A0ht0256+wZKKQNIbnOenraNWhK4CNGHvF4vDQ0NaK2DOS0bN25kwYIFLFmypMOUhN/vl91775Fhhfz4ixx742tcmfkMOfPmB0c/2o96+f1+XK44woUZvi4CF6UUEydOZOLEidhstmBwarVa+2wfINlvqA/oXh1x6Y6DwEvA6DCPjeZOrsphIF8pNUZr3XbIcB6BdJDDbe7rSduoJRObQvQhh8NBdnY2V68GnjNOnz5Nfn4+I0eOxGazdXhB9Xq9Ed37ZqBTCmYn3Cb9g5+w+Y3X7pr3Ypomly8X420XI2oNZZ5RHXJA2v5Ngit1HA6sVmtwFEOCiYHH7++dj256g8DUzYtKqeBQmlJqNYEk2i0td/2m5fMn2x3/ScAEXm9zX0/aRi0JXIToY6NGjeKDDz4AoLKykvj4+OBjbrc77AhLuORQv99PfX192OkI0zQpKipi3759MmIDjHCYLK7bxsbXfhk2eDEMg8cee5zT/ny2nI7H7dE0ezQbzyfxoYKnQvJGmpubqaurCzleghTR27TWFQSq2S4CdrSU5f9H4DWgGPhWS7ujwE+ATyulfqCUelEp9WvgKeBrWuvSNufsdttoJlNFQvQxrTVxcXE0NTUxevToYPJteXk5FRUV5OXlBXckdjgczJgxg/Pnz5OZeScptaGhgW3btgVXrxQUFISM1mzZsoUVK1bg8/lYv349jz766AO/cmSw3cuc8+s5cngys+fMDfl9WSwWDMMgJycH76xZ7Nm1HcPhYsWTizAMI/i783g8rP/lD0kak0Vubm63E2FF7NME8p369Htq/c9KqSrgL4FXgDoCgcvfaK3bbqL1J8BV4I+AF4DLwKcIlPJvrydto5L81wnRx+x2Ozk5OezevZv8/HwyMjK4fft28L7200UpKSkcOXKEsWPHBhMp9+/fz+rVqzEMg/Hjx3PmzBlmzAiUwi8tLWXSpEk4nU6UUixbtoy9e/eydOnSPv9Ze0NFdehCh/xxN0MbNLUbUTLujH4oR2iy6siseE7X1N41idVms2Gz2Vi2/GEsFkvISIppmryzeyf5k6soKr+M1+uVwOVB0vc5LoFvq/VPgZ920cYD/H3LR1fn63bbaPVgvwUTop8YhsHixYvZsWMHO3fu5OjRo2zdupVp06Z1eFE1DIP8/Hy2bNkSnPYxTTM4fRQfH09NTU2w/bFjx5gwYUJwz6P4+HjS09PZsWNHRKu5RiOvnw6bESqlulzt43A4QoIWv99HRUUFKY3vkuiCeF+pBC1C9BP5zxOiH9hsNhITE1m2bFlIMOFwOMKuYLHb7eTl5bFp0yYefvhhnE5n8IVTa01zczNAcHO/ttNCdrud0aNHk5yczJtvvklBQcGAXgLb7IEDJS4aHBkkjlyN3+/n9uWT5BqHqXQbJI1N6nEA5/X6OLr7f1k9ORDwpCf5uHnzZrDmjhj4ermOi7gPErgI0Y+6u7mdxWIhLi6OvLw8NmzYwOLFi4OBi8/no7m5maqqKvbv38/KlSs7jNrY7XZsNhuFhYXs2bOH/Pz8Xv9ZosHFy3Wcu72A5Y8XAIFAEMCzaBG/+ImXqTPnMWfWrB5tKmiaJhtf+zkPj78V3B5g9BDYf/48w4cP73HtFCHE/ZHARYgYYbVaiYuLY82aNSEvlg6Hg4KCAt555x3mz59/1+NbR3IG6iqj2joPF8tGsmJVbofAxGq18tyLn0Br3eOgZcNrr5Ln3UO8Y1zwfrsV6iov4/EskMDlQdFPOS6iIwlchIgx4aZ57HY7Dz30EDabrdPVQw6Hg8TERDweT8ws4f3Vjw+F3F7yj6tCbk+zBTai3HekjvzCD+P3+2lubkYphdVqxev1Bmus9JTFYsHGLRJnTaA+68Mhjzlv7YuZ36HoHRK4RAcJXISIcl6vN1hpV2uNz+dDKYXdbg8JUro7kpCcnExNTQ0pKSmR6nK/SE22sHnT63jczThtJko3oyrLcQ8ax5KV60hISOhxoKG1ZlDabJrMvR0eGz0yjvLyctLS0h74peZC9CUJXISIclprzp49y4kTJ0hJSQkGLk6nk2HDhjF8+HASExO7HbgMGTJkQAYuc6YlAn707VvB+/SREm403eDw4dHdWg7u8Xior68nPj4ewzCw2+3MnjOX4n17GdOu7eRJSby5cRuFqz/So+knEZskOTd6SOAiRBTyeDwopfD7/cGS8sXFxeTl5QWTcls38ysrK6OoqIjFixczdOjQkCkRt9uNw+EIqTlis9kemG0EbrltHI1fSX5ubpejLV6vl0OHDnHtnd+x8qXPYxgGfr+fyspKkuI7tldKkbdsKLt2bWfZsrxgIrAYoCTHJWrI+KYQUchisfDDH/6QvXv3YpomNpuNJUuWcObMmZA2SUlJTJ06lbVr13Lo0KEO5ewdDge7d++mtLQ0GKy0TjsNdFpr9jZnk1/4aLemiGw2G2PGjGGao5rdG1/HNE18Ph+nDm4jbXD4YxIGGcQ5ymWqSIg+JCMuQkQhn8/HQw89xK1bt9i1axdAp1M7Siny8vLYvn07K1asCL77P3nyJFVVVbjdboYNG4bNZsPtdsdUUunSdVNDblc0hgZnTSMmh9x2qUAQ8cG5GhauWt2jn3XEiBGccSUzQx/k+PFpDB8+nKkjbqNcLq7WhU4WDXMNDXy/YSlUV1eTmpra7e8jYo9MFUUPeZsgRBQyDIMpU6aQmZmJz+cjNzeXjIwMZs+efddjHA4Hw4cPp6kpUCLf4/Fw8eJFCgsLyc3NZfPmzQAcOXKEMWPaZ2wMPBevGQwdOrRHx3g8HsxBoxmbYJJ87r8pP/ELJqZ1/jTp8/llxOVBoPt8d2hxF/LfJkSUMgyD1NRUli1bxoYNG1BKdTnFM2vWLHbv3h2cMrLb7bhcruC2AXv27GHu3LkDvvZIaXkTk6fOu7ey/PZBeP1gUZq5Eyxd/s4rym+SlJR0bx0VQvSYBC5CRLHWonOPPPIIe/fu5fbt2122X7JkCZs3b8ZqtYaMBAwaNIilS5eSlpYW4V73L4/Hz4ETbjIyMsKOhJimidfrxTTN8Cew2tl0awos/izrj8Xh8d59XyOPxwMeNWCL+olQWute+RD3R3JchIgBdrudlStX8vbbb1NYWNhp25SUFHJycti8eTMLFizoox72j7LTx1n/7g1sdoXDYUEDNhXHysI1dx1V8vv9vPrqq2RmZjJv3ryQpcxKKRq1i6WPPk9ycjKjRv0Jm373Yx6dVRv2XMfePUbe8hWyougBIDku0UMCFyFihFKKkSNHcuPGDYYPH95p27S0tAEzsnLsYEnI7QVT7yTBmjVVfPjxJ1BKBafHlFKdBhKtv8fBgwdz+PBh5s6dGwxelFK4XC6Sk5OxWCyBERvHUPxOL0OcoUnBtXVekgYNkV2ihehjMlUkRIyw2WxkZWVx8ODB/u5K1BiSPY9XX/15SFG+rkY/HA4HeXl5XLhwgbi4OI4dOxacNvL7/bhcrpBl5Y21pVgsHfNcLp+/yOzZs2NqhZa4D5KcGzUkcBEihlitVtLS0qivr+/vrkQFuyuOtU9ks33n79mx423Onj2L2+3G7XZ3mndiGAarV6+muLiYhIQETpw4gWmawZyY1jwEn8/H2EkLaGzu+GpTf7tGRlseMBK4RAcJXISIIRaLhVGjRlFRUdHfXQnr2rVrvP3221RVVfXZ93TFGSzLy2RyVhxXrh3hzV9/i9d/90v8XbxCtK60qqysZPDgwezcuZN9+4pCqgxbrVbmzZvHpeueDsfblVUSLYXoBxK4CBFjhgwZck+Bi9aapqYmGhoaItCrwPmPHz9OXl4ep06d4ty5cxH5Pu01NrjZ9NZZLP5xLFrwCAuXP8vEyTOorq7uUEm4PavVSl1dHRMnTmT58uWk1Wxl/PjxISMpNpuNyprQLRKam5oZlpIiexQ9QFqTc2XEpf/JOKcQMcZms1FTU9OjY27cuMH+/fsZPXo0SinKysoYMWIEc+fO7bV+VVdXM378eCwWCwsXLmTr1q1MmjTpvrcXeP3v8kNub7ocOprzwdlqClY9Qnx8YEOhhIQE0tPTu7Ufk9PpZOLEiZSUlJCens4lYyn548aFtDl9+jR1agpNZVVMmDABgPcvvt8rP5uIIbJXUdSQwEWIGNSdSq2texJprTlw4ACFhYXBEQK32015eTlnzpxhypQpIccA9/SCbLHcKdZmGAazZ8/m0qVLwRf7SKmubAibINvd/JNx48axceNGzp07R2ZmZofpn4kTJ5KZmcnly5fZvn07ubm5XLt2jczMzF7pvxCiZyRwESLGKKXCTlForSktLeXUqVP4/X5sNhuJiYmMGTOGKVOmhBzjcDhIS0tj69atwcClqamJbdu2AZCfn4/L5epRv5xOZ8g01IgRI9i7d2/EAxfU/W0aaRgGjzzyCF6vF6vV2iEIav29jR07ltTUVDZt2sSkSZOkzP8DRuq4RA8JXISIMaZpdggqqqqqKCoqIisri7y8PPx+P0opiouLOXr0KPPnz+9wHqvVSnx8fHD1zZYtWygsLEQpxYYNG1i9enW3C6tprTl37hw3btzANE0Mw8Dr9Qb3TYqk7NkT2LlzB7m5y+8558Rms3U5QmOz2UhISCA/Px+llKwoetDIVFHUkP88IWJMuMDl4MGDFBYWdhgtGDduHEVFRWGnUqxWK5mZmezZs4f6+npWrFiBYRgopVi9ejU7duzA4XCQlJREYmIiCQkJJCYm4nQ6Q0Y4/H4/27dvZ/r06WRnZwdX80S6vLnX3YzVbjBkaCKmvtBn+SZSJVeI/iWBixAxprm5GZfLFZwaqq6uxjRNfD5f2GmOlStXBhNX20tJScHpdDJo0CBsNlvwxb+1SJvdbsftdlNfX09tbS1Xr16lpqYGr9cbkg+zaNEi4uPjQ/Jc7HY7hmHg9/f+7sl1N0pJM29TXnmLovM+vE0WWZosIs4vl1hUkMBFiBjT1NTEsWPHOHv2LJMnTyY9PZ2SkpK7VnAdMWLEXc9lGAZDhw4N+1jryILL5cLlcjFs2DAmTJiAz+frELi0HYVoDVxaR3SKi4t7Pc+lvvQqi9cWBGuueL1eWZosIkpyXKKHBC5CxJi0tDTWrl2Lz+cLBgx5eXn4fL6wGwvebbPBe2W1Wrt9zuHDh7N+/XrGjh3bqzkh8cqHx+PB6XQG+ySEeDBIWrwQMcYwDGw2W8goh2EYUTniYLVaWbVqFZs2beLMmTO9Np3jsltlVY/oW7JXUdSQ/3whRMRYLBbi4uJYtWoV8fHxvPnmm1RWVt73eWscg7l06RJerze4QaIQkSSVc6OHBC5CiIgzDINRo0YFNzbcsWNHtyrb3k3KxKkcq2ri529u5Njx4xK8CPEAkRwXIUSfsFgsWCwWZs2ahWma7Ny5k/z8/K4PbKdwXHLgi3HJwGwOHz5MY2NjVE6ViYFFRkuigwQuQog+1br6KSUlhfLy8g6rnnw+H5WVlcTFxTFo0KAu67M0NTXddUWVEL2lGPeWZzmX0kunu//50geYBC5CiD5nt9vJzs5m8+bNPProo2ituXr1KqdOncLhcDB27FjKysq4fv06GRkZTJ06tcM5vF4vBw4cIDExUQIXEXFa64L+7oMIUAOxaFNOTo4+dOhQf3dDCNEJv9/PsWPHuH37NvX19WRnZ5Oeng7cGZUxTZOysjIuXLjAihUrgscdP36c0tJSli1bFrUrqkRMkO29Y5CMuAgh+oXFYmH69OnU1taSkJCA3W7vMC1kGAajR4+mvr6emzdvUl1dzdmzZ1myZAlZWVlSfl+IB5CsKhJC9Bu73U5ycnJwj6RwrFYrEyZMYO/evcTHx1NQUEBycrIELUI8oGTERQgR9axWK2vWrJFgRQghgYsQIvpJ8q0QopVMFQkhhBAiZkjgIoQQQoiYIYGLEEIIIWKGBC5CCCGEiBkSuAghhBAiZkjgIoQQQoiYIYGLEEIIIWKGBC5CCCGEiBkSuAghhBAiZkjgIoQQQoiYIYGLEEIIIWKGBC5CCCGEiBkSuAghhBAiZkjgIoQQQoiYIYGLEEIIIWKG0lr3dx96nVKqArjS3/0QQggR1Sq11gX93QnRMwMycBFCCCHEwCRTRUIIIYSIGRK4CCGEECJmSOAihBBCiJghgUs/UUoZSinV3/0QA49cWyKS5PoS/e2BC1yUUoOUUl9WSm1USlUopbRS6svdOO6jLW293Wj7xZa2F8I8ZiilfgrUA7eUUp9o93huy7F3+/hq939a0ZcieW0ppdYqpXYrpWqVUnVKqSNKqefatZFrawCL1PWllJqnlFqvlCpTStUrpU4ppf5WKeVq106uLxEVbP3dgX6QAnwJuA4cAVZ2dYBSKhH4BtAAOLtoOw74m5a24XwayAGea+nLPyqlDmut32vX7vvAvjDHv99Vf0W/ici1pZT6K+CfgNeBzwF+IBMY066pXFsDW69fX0qphcAe4DzwCoGgZBnw/4D5wGNtmsv1JaLCgxi4lAGjtNalSqnRwLVuHPMPwC1gO/BMF22/DewFDGB0mMcfAj6vtX4DQCk1ueW+9v/8+7TWr3ajbyJ69Pq1pZSaC3wN+Cut9Te7OJdcWwNbJJ67WkdNlmqtq1q+/oFSygZ8WCk1RGtd3XK/XF8iKjxwU0Vaa7fWurS77ZVSM4A/Az4FdDpNpJRaDaxuaXs3F4E/UUpNVUotAdYC57rbHxG9InRtfRq4CfxLyzEJnZxSrq0BLELXVyLQDFS3u7+MwMieu819cn2JqPDABS490ZKA9l3gLa312120dQLfAf5da/1BJ02/RmCI/zSBkZl3tNZvhWk3SCmVEubDfm8/jYgmPbi28oGDwJ8qpW4CtUqpSqXUV5RS7f9/5doSQI+ur11AAvAjpVSWUmqMUup54EXgX7TWjW3ayvUlosKDOFXUE88TmNPN6kbbvybwBPDlzhpprcuVUrOA6UCt1rpDAm+L77d8tJcPdBpEiZjQ5bWllEoChgGLgBUEhv3PA08CXwTigc+0tpdrS7TR3eeu7wKTgY8Bf9Rynwa+pLX+h7YN5foS0UICl7tQSg0mkNT2ita6uIu24wkkTb6sta7p6txaaw+B5LrOfJ3AvHR7XR0nolwPrq1BLZ9TgGe11r9suf3blnP8uVLqa1rrytYD5NoSPXnuAnxAMfAWgeTvZmAd8BWlVI3W+jttG8v1JaKBBC5391XAJDA82pXvACeB/+rF73+qq+kpEbO6e201tXz2Ar9u99gvgEICKz829vD7y7U1sPXkuesVYA0wXWttttz3WmCmiW8opX6ttb7Rw+8v15eIKAlcwlBKTQL+lMCy5lR1p9bSoJbHxwHNLUOny4FHgI8AY9q0dQK2lrb1bd8ViwdXT64tAgmTTUCN1trX7lStLyZDIt1nETt6+NxlB14mkJdntjvV74BngXnA+j7ouhDdJoFLeKMAK4Hh1m+EebwY2A3kAmNb7vtlmHatbf8beKFXeyhiVbevLa21Xyl1DJivlDLavbi0LrWviGRnRczpyXNXCmBvad+erd1nIaKGXJThnQQeD3P/JwkUZ/oDoHUEZcdd2n4VSCbw7udqBPooYlNPri2A/yGQnPsiLQmPLauJ/hioA96NZGdFzOnJ9XWTwKjeY0qpz2ut69u0f5ZAku7RCPZViHvyQAYuSqmXgSQCNQwAliqlvtDy9c+11lcIJKq1P+4xQGutg49pra8SJjBRSv0F4GzbtocWq/DbgVzRWr9zj+cUEdab11aL/yQQtHynpeDXBQIvTB8CPqW1rruHbsq1FaN6+bnLp5T6OoGqzIeUUj8mMDX5GIFVbD9pOV9PyfUlIkpprfu7D31OKXWZO1M87S3XWu+6y3E/BZ7TWncZ8CmldgGjtdYTe9i3XGBnJ01+q7V+sifnFH0nEteWUiqZQKLlY8Bg4CzwTa31z3vYt1zk2oppEbq+ngD+EphCICC6SGB6+5+11l3uzdbmPLnI9SX6wAMZuAghhBAiNknlXCGEEELEDAlchBBCCBEzJHARQgghRMyQwEUIIYQQMUMCFyGEEELEDAlchBBCCBEzJHARQgghRMyQwEUIIYQQMUMCFyGEEELEDAlchBBCCBEz/j/A9n5EyUs2NgAAAABJRU5ErkJggg==\n",
"text/plain": [
"