{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import emipy as ep" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "ep.change_RootPath(r'C:\\Path\\to\\your\\project')\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
FacilityReportIDPollutantReleaseAndTransferReportIDFacilityIDNationalIDParentCompanyNameFacilityNameStreetNameBuildingNumberCityPostalCode...PollutantNamePollutantGroupCodePollutantGroupNamePollutantCASMethodBasisCodeMethodBasisNameTotalQuantityAccidentalQuantityUnitCodeUnitName
3683111007134506-01-01000004838Stadtwerke Flensburg GmbHStadtwerke Flensburg GmbHBatteriestr. 48NaNFlensburg24939...Sulphur oxides (SOx/SO2)OTHGASOther gasesNaNEEstimated582000.00.0KGMkilogram
3684111007134506-01-01000004838Stadtwerke Flensburg GmbHStadtwerke Flensburg GmbHBatteriestr. 48NaNFlensburg24939...Nitrogen oxides (NOx/NO2)OTHGASOther gasesNaNEEstimated378000.00.0KGMkilogram
3685111017134706-01-02000283604Gemeinschaftskraftwerk Kiel GmbH -GKK-Gemeinschaftskraftwerk Kiel GmbH -GKK-Hasselfelde 40NaNKiel24149...Chlorine and inorganic compounds (as HCl)OTHGASOther gasesNaNMMeasured41000.00.0KGMkilogram
3686111017134706-01-02000283604Gemeinschaftskraftwerk Kiel GmbH -GKK-Gemeinschaftskraftwerk Kiel GmbH -GKK-Hasselfelde 40NaNKiel24149...Nitrogen oxides (NOx/NO2)OTHGASOther gasesNaNMMeasured1500000.00.0KGMkilogram
3687111017134706-01-02000283604Gemeinschaftskraftwerk Kiel GmbH -GKK-Gemeinschaftskraftwerk Kiel GmbH -GKK-Hasselfelde 40NaNKiel24149...Sulphur oxides (SOx/SO2)OTHGASOther gasesNaNMMeasured960000.00.0KGMkilogram
\n", "

5 rows × 73 columns

\n", "
" ], "text/plain": [ " FacilityReportID PollutantReleaseAndTransferReportID FacilityID \\\n", "3683 11100 7 1345 \n", "3684 11100 7 1345 \n", "3685 11101 7 1347 \n", "3686 11101 7 1347 \n", "3687 11101 7 1347 \n", "\n", " NationalID ParentCompanyName \\\n", "3683 06-01-01000004838 Stadtwerke Flensburg GmbH \n", "3684 06-01-01000004838 Stadtwerke Flensburg GmbH \n", "3685 06-01-02000283604 Gemeinschaftskraftwerk Kiel GmbH -GKK- \n", "3686 06-01-02000283604 Gemeinschaftskraftwerk Kiel GmbH -GKK- \n", "3687 06-01-02000283604 Gemeinschaftskraftwerk Kiel GmbH -GKK- \n", "\n", " FacilityName StreetName BuildingNumber \\\n", "3683 Stadtwerke Flensburg GmbH Batteriestr. 48 NaN \n", "3684 Stadtwerke Flensburg GmbH Batteriestr. 48 NaN \n", "3685 Gemeinschaftskraftwerk Kiel GmbH -GKK- Hasselfelde 40 NaN \n", "3686 Gemeinschaftskraftwerk Kiel GmbH -GKK- Hasselfelde 40 NaN \n", "3687 Gemeinschaftskraftwerk Kiel GmbH -GKK- Hasselfelde 40 NaN \n", "\n", " City PostalCode ... PollutantName \\\n", "3683 Flensburg 24939 ... Sulphur oxides (SOx/SO2) \n", "3684 Flensburg 24939 ... Nitrogen oxides (NOx/NO2) \n", "3685 Kiel 24149 ... Chlorine and inorganic compounds (as HCl) \n", "3686 Kiel 24149 ... Nitrogen oxides (NOx/NO2) \n", "3687 Kiel 24149 ... Sulphur oxides (SOx/SO2) \n", "\n", " PollutantGroupCode PollutantGroupName PollutantCAS MethodBasisCode \\\n", "3683 OTHGAS Other gases NaN E \n", "3684 OTHGAS Other gases NaN E \n", "3685 OTHGAS Other gases NaN M \n", "3686 OTHGAS Other gases NaN M \n", "3687 OTHGAS Other gases NaN M \n", "\n", " MethodBasisName TotalQuantity AccidentalQuantity UnitCode UnitName \n", "3683 Estimated 582000.0 0.0 KGM kilogram \n", "3684 Estimated 378000.0 0.0 KGM kilogram \n", "3685 Measured 41000.0 0.0 KGM kilogram \n", "3686 Measured 1500000.0 0.0 KGM kilogram \n", "3687 Measured 960000.0 0.0 KGM kilogram \n", "\n", "[5 rows x 73 columns]" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "db = ep.read_db()\n", "data1 = ep.f_db(db,CountryName='Germany')\n", "data1.head()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "ename": "FileNotFoundError", "evalue": "[WinError 3] Das System kann den angegebenen Pfad nicht finden: 'C:\\\\Path\\\\to\\\\your\\\\project\\\\MappingData'", "output_type": "error", "traceback": [ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[1;31mFileNotFoundError\u001b[0m Traceback (most recent call last)", "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[0mexample_res\u001b[0m \u001b[1;33m=\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m10\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 2\u001b[1;33m \u001b[0mep\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdownload_MapData\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34mr'C:\\Path\\to\\your\\project'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mresolution\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mexample_res\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;32mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[1;32m~\\anaconda3\\envs\\emipy\\lib\\site-packages\\emipy\\rawdata.py\u001b[0m in \u001b[0;36mdownload_MapData\u001b[1;34m(path, resolution, clear, chunk_size)\u001b[0m\n\u001b[0;32m 67\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 68\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mos\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0misdir\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mis\u001b[0m \u001b[1;32mFalse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 69\u001b[1;33m \u001b[0mos\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmkdir\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 70\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 71\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mclear\u001b[0m \u001b[1;33m==\u001b[0m \u001b[1;32mTrue\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;31mFileNotFoundError\u001b[0m: [WinError 3] Das System kann den angegebenen Pfad nicht finden: 'C:\\\\Path\\\\to\\\\your\\\\project\\\\MappingData'" ] } ], "source": [ "example_res =[10]\n", "ep.download_MapData(r'C:\\Path\\to\\your\\project', resolution=example_res, clear=True)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "data2 = ep.change_unit(data1, unit='megaton')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
CountryCodeCountryNameLatLongNUTSRegionGeoCodeNACEMainEconomicActivityCodeNACEMainEconomicActivityNameReportingYearPollutantReleaseIDPollutantNameTotalQuantityUnitCode
3683DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017341Sulphur oxides (SOx/SO2)582000.0KGM
3684DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017340Nitrogen oxides (NOx/NO2)378000.0KGM
3685DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140850Chlorine and inorganic compounds (as HCl)41000.0KGM
3686DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140852Nitrogen oxides (NOx/NO2)1500000.0KGM
3687DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140853Sulphur oxides (SOx/SO2)960000.0KGM
\n", "
" ], "text/plain": [ " CountryCode CountryName Lat Long NUTSRegionGeoCode \\\n", "3683 DE Germany 54.805185 9.432377 DEF0 \n", "3684 DE Germany 54.805185 9.432377 DEF0 \n", "3685 DE Germany 54.340605 10.178735 DEF0 \n", "3686 DE Germany 54.340605 10.178735 DEF0 \n", "3687 DE Germany 54.340605 10.178735 DEF0 \n", "\n", " NACEMainEconomicActivityCode NACEMainEconomicActivityName \\\n", "3683 NACE_1.1:40.10 Production and distribution of electricity \n", "3684 NACE_1.1:40.10 Production and distribution of electricity \n", "3685 NACE_1.1:40.30 Steam and hot water supply \n", "3686 NACE_1.1:40.30 Steam and hot water supply \n", "3687 NACE_1.1:40.30 Steam and hot water supply \n", "\n", " ReportingYear PollutantReleaseID \\\n", "3683 2001 7341 \n", "3684 2001 7340 \n", "3685 2001 40850 \n", "3686 2001 40852 \n", "3687 2001 40853 \n", "\n", " PollutantName TotalQuantity UnitCode \n", "3683 Sulphur oxides (SOx/SO2) 582000.0 KGM \n", "3684 Nitrogen oxides (NOx/NO2) 378000.0 KGM \n", "3685 Chlorine and inorganic compounds (as HCl) 41000.0 KGM \n", "3686 Nitrogen oxides (NOx/NO2) 1500000.0 KGM \n", "3687 Sulphur oxides (SOx/SO2) 960000.0 KGM " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data3 = ep.row_reduction(data1)\n", "data3.head()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
CountryCodeCountryNameLatLongNUTSRegionGeoCodeNACEMainEconomicActivityCodeNACEMainEconomicActivityNameReportingYearPollutantReleaseIDPollutantNameTotalQuantityUnitCode
3683DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017341Sulphur oxides (SOx/SO2)582000.0KGM
3684DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017340Nitrogen oxides (NOx/NO2)378000.0KGM
3685DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140850Chlorine and inorganic compounds (as HCl)41000.0KGM
3686DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140852Nitrogen oxides (NOx/NO2)1500000.0KGM
3687DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140853Sulphur oxides (SOx/SO2)960000.0KGM
\n", "
" ], "text/plain": [ " CountryCode CountryName Lat Long NUTSRegionGeoCode \\\n", "3683 DE Germany 54.805185 9.432377 DEF0 \n", "3684 DE Germany 54.805185 9.432377 DEF0 \n", "3685 DE Germany 54.340605 10.178735 DEF0 \n", "3686 DE Germany 54.340605 10.178735 DEF0 \n", "3687 DE Germany 54.340605 10.178735 DEF0 \n", "\n", " NACEMainEconomicActivityCode NACEMainEconomicActivityName \\\n", "3683 NACE_1.1:40.10 Production and distribution of electricity \n", "3684 NACE_1.1:40.10 Production and distribution of electricity \n", "3685 NACE_1.1:40.30 Steam and hot water supply \n", "3686 NACE_1.1:40.30 Steam and hot water supply \n", "3687 NACE_1.1:40.30 Steam and hot water supply \n", "\n", " ReportingYear PollutantReleaseID \\\n", "3683 2001 7341 \n", "3684 2001 7340 \n", "3685 2001 40850 \n", "3686 2001 40852 \n", "3687 2001 40853 \n", "\n", " PollutantName TotalQuantity UnitCode \n", "3683 Sulphur oxides (SOx/SO2) 582000.0 KGM \n", "3684 Nitrogen oxides (NOx/NO2) 378000.0 KGM \n", "3685 Chlorine and inorganic compounds (as HCl) 41000.0 KGM \n", "3686 Nitrogen oxides (NOx/NO2) 1500000.0 KGM \n", "3687 Sulphur oxides (SOx/SO2) 960000.0 KGM " ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "total=['CountryCode', 'CountryName']\n", "ep.change_ColumnsOfInterest(total=total, add=None, sub=None, reset=False)\n", "data4 = ep.row_reduction(data3)\n", "ep.change_ColumnsOfInterest(total=None, add='Lat', sub=None, reset=False)\n", "data5 = ep.row_reduction(data3)\n", "ep.change_ColumnsOfInterest(total=None, add=None, sub='CountryCode', reset=False)\n", "data6 = ep.row_reduction(data3)\n", "ep.change_ColumnsOfInterest(total=None, add=None, sub=None, reset=True)\n", "data7 = ep.row_reduction(data3)\n", "data7.head()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
CountryCodeCountryLatLongNUTSIDNACEIDNACENameYearPollutantReleaseIDPollutantTotalQuantityUnit
3683DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017341Sulphur oxides (SOx/SO2)582000.0KGM
3684DEGermany54.8051859.432377DEF0NACE_1.1:40.10Production and distribution of electricity20017340Nitrogen oxides (NOx/NO2)378000.0KGM
3685DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140850Chlorine and inorganic compounds (as HCl)41000.0KGM
3686DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140852Nitrogen oxides (NOx/NO2)1500000.0KGM
3687DEGermany54.34060510.178735DEF0NACE_1.1:40.30Steam and hot water supply200140853Sulphur oxides (SOx/SO2)960000.0KGM
\n", "
" ], "text/plain": [ " CountryCode Country Lat Long NUTSID NACEID \\\n", "3683 DE Germany 54.805185 9.432377 DEF0 NACE_1.1:40.10 \n", "3684 DE Germany 54.805185 9.432377 DEF0 NACE_1.1:40.10 \n", "3685 DE Germany 54.340605 10.178735 DEF0 NACE_1.1:40.30 \n", "3686 DE Germany 54.340605 10.178735 DEF0 NACE_1.1:40.30 \n", "3687 DE Germany 54.340605 10.178735 DEF0 NACE_1.1:40.30 \n", "\n", " NACEName Year PollutantReleaseID \\\n", "3683 Production and distribution of electricity 2001 7341 \n", "3684 Production and distribution of electricity 2001 7340 \n", "3685 Steam and hot water supply 2001 40850 \n", "3686 Steam and hot water supply 2001 40852 \n", "3687 Steam and hot water supply 2001 40853 \n", "\n", " Pollutant TotalQuantity Unit \n", "3683 Sulphur oxides (SOx/SO2) 582000.0 KGM \n", "3684 Nitrogen oxides (NOx/NO2) 378000.0 KGM \n", "3685 Chlorine and inorganic compounds (as HCl) 41000.0 KGM \n", "3686 Nitrogen oxides (NOx/NO2) 1500000.0 KGM \n", "3687 Sulphur oxides (SOx/SO2) 960000.0 KGM " ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "addition={'Lat': 'Latitude'}\n", "ep.change_RenameDict(total=None, add=addition, sub=None, reset=False)\n", "data8 = ep.rename_columns(data3)\n", "\n", "ep.change_RenameDict(total=None, add=None, sub=None, reset=True)\n", "data9 = ep.rename_columns(data3)\n", "data9.head()" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1gAAAI4CAYAAAB3HEhGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABxq0lEQVR4nO3deZgcdbX4//chCSQBZA2CLAaQHTTESVBAFpEIoiAXUBGBgJDrvYbF5QooXwUURUVABC8/FAQUEUQ2F0SQTeACCSEsYRMQMBAgYY+CkuT8/qjKMJlMku6kZqZ65v16nnqmp7r69Kmqrv70qeVTkZlIkiRJkpbcUr2dgCRJkiT1FRZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkioysLfeeNVVV83hw4f31ttLkmrsrrvumpGZw3o7D9sqSdKCLKit6rUCa/jw4UycOLG33l6SVGMR8WRv5wC2VZKkBVtQW+UpgpIkSZJUEQssSZIkSaqIBZYkSZIkVaTXrsGSpJ725ptvMnXqVN54443eTkWlwYMHs9ZaazFo0KDeTkVSD/M7Wa2i2bbKAktSvzF16lSWX355hg8fTkT0djr9XmbywgsvMHXqVNZdd93eTkdSD/M7Wa1gcdoqTxGU1G+88cYbrLLKKjbkNRERrLLKKu69lvopv5PVChanrbLAktSv2JDXi+tD6t/8DlAraPZzaoElSZIkSRWp1TVYw4/+/SKneeKk3XogE0n9QSPfOc1Y1PfTCy+8wE477QTAs88+y4ABAxg2rLgB/J133snSSy/dPu1pp53GuHHjGDp06EJj7rDDDpx88sm0tbXxyiuvcNhhh3HbbbeRmbzvfe/jjDPOYKWVVlrCOZvX5MmTeeaZZ/jIRz4CwFVXXcUDDzzA0UcfzRVXXMGGG27IpptuWul71oXtlNR9evo7GYojE1/84hf5wQ9+AMDJJ5/MzJkzOe644zjrrLMYOnQoBxxwAOeddx5jxozhHe94R6U59oSJEydywQUXcPrpp8/33Nybqa+66qpL9B6nnXYaK6+8MgcccABjx47l2muv5fHHH2eZZZZhxowZtLW18cQTTwAwZcoUDjvsMJ5++mnmzJnDAQccwLHHHrvIo0TTpk3jwAMP5Oyzz2bdddfl9NNP57DDDgNg/PjxtLW1MXbsWDKTE088kfPPP5+IYM011+SMM85gs80245///Cf77LMPjz32GAMGDOBjH/sYJ510EgBnnHEGQ4cO5eCDD16iZQEewZKkHrPKKqswefJkJk+ezOc+9zm+8IUvtP/fsbiCorH65z//2VT8z372s6y33no8+uijPPbYY7zrXe9i7NixFc5BYfLkyfzhD39o/3/33Xfn6KOPBuCKK67ggQceqPw9Jak7LLPMMlx22WXMmDFjvuc+97nPccABBwBw3nnn8cwzz3QZY/bs2d2a45Jqa2vrsriqyqxZszj33HP59Kc/3T5uwIABnHvuufNN+/rrr7e3GQ8//DD33HMPt912Gz/+8Y8X+T5//OMf+fCHPwzAaqutxg9/+EP+/e9/zzfdmWeeyW233cY999zDI488wjHHHMPuu+/efg3Vl7/8ZR566CHuvvtubr31Vq6++moADj74YH70ox8t1jLozAJLknrRn//8Z7bccku22GILDj74YP71r39x+umn88wzz7Djjjuy4447AvBf//VftLW1sdlmm/GNb3xjvjiPPvood911F//v//2/9nFf//rXueeee3j44Ye58cYb+ehHP9r+3Pjx4znvvPMAOOGEExg1ahSbb74548aNIzOB4ujYUUcdxejRo9lwww35y1/+wr///W++/vWvc/HFFzNixAguvvhizjvvPMaPH89tt93GVVddxf/8z/8wYsQIHnvsMUaOHNn+nn/961/n+V+SetvAgQMZN24cp5566nzPHXfccZx88slceumlTJw4kf32248RI0bw+uuvM3z4cI466ihGjhzJr3/9ay666CK22GILNt98c4466qj2GOeccw4bbrgho0eP5tBDD2X8+PEATJ8+nb322otRo0YxatQobr311vb3PPjgg9lhhx1Yb731FlgYdfV+l19+OTvttBOZybRp09hwww159tln5/n+f+GFFxgzZgybbbYZhxxySPv3PcAvfvELRo8ezYgRI/jP//xPZs+ezezZsxk7diybb745W2yxRZfL6frrr2fkyJEMHPjWiXFHHnkkp556KrNmzZpn2l/+8pdss802jBkzBoChQ4dyxhlntB9FOuKIIzjhhBMAuOaaa9huu+2YM2cOUBRYu+66KwDDhg1jp5124vzzz58vn+9+97vtR6MAxowZw9Zbb82FF17I0KFD29vVpZdempEjRzJ16tT2XIYPH86dd97Z5TJvhgWWJPWSN954g7Fjx3LxxRdz3333MWvWLP73f/+Xww8/nHe84x3ccMMN3HDDDQCceOKJTJw4kXvvvZebbrqJe++9d55YDzzwACNGjGDAgAHt4wYMGMCWW27Jgw8+uNA8xo8fz4QJE7j//vt5/fXX+d3vftf+3KxZs7jzzjs57bTTOP7441l66aU54YQT+OQnP8nkyZP55Cc/2T7t1ltvze677873v/99Jk+ezPrrr88KK6zA5MmTAfjZz37GQQcdtKSLTZIq9fnPf54LL7yQV155pcvn9957b9ra2rjwwguZPHkyQ4YMAYqzEiZNmsR2223HUUcdxfXXX8/kyZOZMGECV1xxBc888wzf/OY3uf3227n11lt56KGH2mMeccQRfOELX2DChAn85je/4ZBDDml/7qGHHuKaa67hzjvv5Pjjj+fNN9+cJ59nnnmmy/fbc889WWONNTjzzDM59NBDOf7441l99dXnee3xxx/Ptttuy5QpU9hzzz156qmnAHjwwQe5+OKLufXWW5k8eTIDBgxon9+nn36a+++/n/vuu6/L7/Bbb72V9773vfOMW2edddh22235+c9/Ps/4KVOmzDft+uuvz8yZM3n11Vf5zne+w8UXX8wNN9zA4Ycfzs9+9jOWWmopZs+ezcMPPzzP6edHHXUUJ5988jxHEF999VX+8Y9/sN56683zHm1tbUyZMmWecS+//DK//e1v20/dnzvdX/7yl/nmsVkWWJLUS2bPns26667LhhtuCMCBBx7IzTff3OW0l1xyCSNHjmTLLbdkypQplZ6Gd8MNN7DVVluxxRZbcP3118/TCP3Hf/wHAO9973vbz59vxiGHHMLPfvYzZs+ezcUXXzzPKSSSVAdve9vbOOCAA5o+jW7uDqYJEyawww47MGzYMAYOHMh+++3HzTffzJ133sn222/PyiuvzKBBg9hnn33aX3vdddcxfvx4RowYwe67786rr77KzJkzAdhtt91YZpllWHXVVVlttdV47rnn5nnfBb0fwI9+9CO+853vsMwyy7DvvvvOl/PNN9/MZz7zmfb3mXuN7p///GfuuusuRo0axYgRI/jzn//M448/znrrrcfjjz/OYYcdxh//+Efe9ra3zRdz2rRp7dcTd3TMMcfw/e9/v/0IVCOGDh3KT37yE3beeWfGjx/P+uuvD8Add9zBVlttNc+06623HltttRW//OUvG44/16xZs9h33305/PDD5ynGVltttQWeCtqMWnVyIUma39/+9jdOPvlkJkyYwEorrcTYsWPnux/HpptuyuTJk5kzZw5LLVXsO5szZw733HMPI0eO5KmnnpqnkZv7+jfeeIP//u//ZuLEiay99tocd9xx88ReZpllgOJoWOdTPRqx1157cfzxx/PBD36Q9773vayyyipNx5Ck7nbkkUcycuTIpo6yL7vssov9fnPmzOH2229n8ODB8z0393sXmv/unTp1KksttRTPPffcPO3BomQmBx54IN/5znfme+6ee+7hmmuu4ayzzuKSSy6Z79qqIUOGdHmPqA022IARI0ZwySWXtI/bdNNN59uR+Pjjj7Pccsu1F2/33Xcfq6yyyjyFztVXX80uu+wy33t89atfZe+992b77bcHimJ52WWXbS8O57rrrrvapwEYN24cG2ywAUceeeQ88d544432I5RLwiNYktRLBgwYwBNPPMGjjz4KwM9//vP2BmD55ZfntddeA4pTHpZddllWWGEFnnvuufYLcjt617vexZZbbsm3vvWt9nHf+ta32GmnnVhnnXV45zvfyQMPPMC//vUvXn75Zf785z8DbxVaq666KjNnzuTSSy9dZN4dc1vUc4MHD+bDH/4w//Vf/+XpgZJqa+WVV+YTn/gE55xzTpfPL+x7b/To0dx0003MmDGD2bNnc9FFF7H99tszatQobrrpJl566SVmzZrFb37zm/bXjBkzZp4OFeaeSt2IBb3frFmzOPjgg7nooovYZJNNOOWUU+Z77Xbbbdd+xOfqq6/mpZdeAmCnnXbi0ksv5fnnnwfgxRdf5Mknn2TGjBnMmTOHvfbai29961tMmjRpvpibbLJJezvW2de+9jVOPvnk9v/3228/brnlFq677jqg6PTi8MMP5ytf+QoATz75JD/4wQ+4++67ufrqq7njjjuA4gjbhz70ofnib7zxxmy66ab89re/bR/3P//zPxx++OG8/vrrQHG08JZbbmk/g+LYY4/llVde4bTTTpsv3iOPPMLmm2/e5bw0wyNYkvqt3u5Oe/DgwfzsZz9jn332YdasWYwaNYrPfe5zQLF3bZdddmm/FmvLLbdk4403Zu2112abbbbpMt65557LYYcdxvrrr8+rr77KqFGj2hudtddem0984hNsvvnmrLvuumy55ZYArLjiihx66KFsvvnmrL766owaNWqRee+4446cdNJJjBgxgmOOOWae5z71qU9x6KGHcvrpp3PppZey/vrrs99++3H55Ze3X9QsSV3p7e/kL33pS5xxxhldPjd27Fg+97nPMWTIEP7v//5vnufWWGMNTjrpJHbccUcyk91224099tgDKI6wjB49mpVXXpmNN96YFVZYAYDTTz+dz3/+87z73e9m1qxZbLfddpx11lkN5bmg9zvhhBP4wAc+wLbbbst73vMeRo0axW67zbtMv/GNb7Dvvvuy2WabsfXWW7POOusAxZGlb33rW4wZM4Y5c+YwaNAgzjzzTIYMGcJBBx3UfgZEV0e4dt11V/bff/8uc91ss80YOXJke2E2ZMgQrrzySg477DA+//nPM3v2bPbff3/Gjx9PZvLZz36Wk08+mXe84x2cc845jB07lj/84Q8MHjyY5Zdfvsv3+NrXvtbepgEcdthhvPTSS2yxxRYMGDCA1VdfnSuvvJIhQ4YwdepUTjzxRDbeeOP2TpfGjx/ffg3crbfeynHHHdfQeliY6Nh7SE9qa2vLiRMnzjPO+4tI6k4PPvggm2yySW+n0SMefvhhdtttN04//fT2+1X1lpNPPplXXnmFb37zm10+39V6iYi7MrOtJ/JbmM5tle2UVJ3+8J08c+ZMlltuOWbNmsWee+7JwQcfzJ577tnbaVVuzz335Hvf+x4bbLBB5bF/8YtfMHXq1PbbgXSXu+++m1NOOWW+jjnmaqat8giWJPVBG2200QJP2ehJe+65J4899hjXX399b6ciST3uuOOO47rrruONN95gzJgxfPzjH+/tlLrFSSedxLRp07qlwJrbKUd3mzFjxgJ3BDbLAkuS1G0uv/zy3k5BknpNx+uP+rKNNtqIjTbaqLfTWCI777xzZbHs5EJSv9Jbp0Wra64PqX/zO0CtoNnPqQWWpH5j8ODBvPDCCzboNZGZvPDCC112Uyyp7/M7Wa1gcdoqTxGU1G+stdZaTJ06lenTp/d2KioNHjyYtdZaq7fTkNQL/E5Wq2i2rbLAktRvDBo0iHXXXbe305Ak4Xey+i5PEZQk9RsRsWJEXBoRD0XEgxHx/t7OSZLUt3gES5LUn/wQ+GNm7h0RSwNDezshSVLfYoElSeoXImIFYDtgLEBm/hv4d2/mJEnqezxFUJLUX6wLTAd+FhF3R8RPI2LZzhNFxLiImBgRE734XpLULAssSVJ/MRAYCfxvZm4J/AM4uvNEmXl2ZrZlZtuwYcN6OkdJUouzwJIk9RdTgamZeUf5/6UUBZckSZVpusCKiHMj4vmIuL/DuJUj4tqI+Gv5d6Vq05Qkaclk5rPA3yNio3LUTsADvZiSJKkPWpwjWOcBu3QadzTw58zcAPgzXZxyIUlSDRwGXBgR9wIjgG/3bjqSpL6m6V4EM/PmiBjeafQewA7l4/OBG4GjliQxSZKqlpmTgbbezkOS1HdVdQ3W2zNzWvn4WeDtXU1kz0ySJEmS+rLKO7nIzARyAc/ZM5MkSZKkPquqAuu5iFgDoPz7fEVxJUmSJKllVFVgXQUcWD4+ELiyoriSJEmS1DIWp5v2i4D/AzaKiKkR8VngJGDniPgr8KHyf0mSJEnqVxanF8F9F/DUTkuYiyRJkiS1tMo7uZAkSZKk/soCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUEQssSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUkYG9nYAkST0lIp4AXgNmA7Mys613M5Ik9TUWWJKk/mbHzJzR20lIkvomTxGUJEmSpIpYYEmS+pME/hQRd0XEuK4miIhxETExIiZOnz69h9OTJLU6CyxJUn+ybWaOBHYFPh8R23WeIDPPzsy2zGwbNmxYz2coSWppFliSpH4jM58u/z4PXA6M7t2MJEl9jQWWJKlfiIhlI2L5uY+BMcD9vZuVJKmvsRdBSVJ/8Xbg8oiAov37ZWb+sXdTkiT1NRZYkqR+ITMfB97T23lIkvo2TxGUJEmSpIpYYEmSJElSRSotsCLiCxExJSLuj4iLImJwlfElSZIkqc4qK7AiYk3gcKAtMzcHBgCfqiq+JEmSJNVd1acIDgSGRMRAYCjwTMXxJUmSJKm2Kiuwyps3ngw8BUwDXsnMP3WcJiLGRcTEiJg4ffr0qt5akiRJkmqhylMEVwL2ANYF3gEsGxGf6ThNZp6dmW2Z2TZs2LCq3lqSJEmSaqHKUwQ/BPwtM6dn5pvAZcDWFcaXJEmSpFqrssB6CnhfRAyNiAB2Ah6sML4kSZIk1VqV12DdAVwKTALuK2OfXVV8SZIkSaq7gVUGy8xvAN+oMqYkSZIktYqqu2mXJEmSpH7LAkuSJEmSKmKBJUmSJEkVscCSJEmSpIpYYEmSJElSRSywJEmSJKkiFliSJEmSVBELLEmSJEmqiAWWJEmSJFXEAkuS1K9ExICIuDsiftfbuUiS+h4LLElSf3ME8GBvJyFJ6psssCRJ/UZErAXsBvy0t3ORJPVNFliSpP7kNOArwJxezkOS1EdZYEmS+oWI+CjwfGbetYjpxkXExIiYOH369B7KTpLUV1hgSZL6i22A3SPiCeBXwAcj4hedJ8rMszOzLTPbhg0b1tM5SpJanAWWJKlfyMxjMnOtzBwOfAq4PjM/08tpSZL6GAssSZIkSarIwN5OQJKknpaZNwI39nIakqQ+yCNYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUqLbAiYsWIuDQiHoqIByPi/VXGlyRJkqQ6q/pGwz8E/piZe0fE0sDQiuNLkiRJUm1VVmBFxArAdsBYgMz8N/DvquJLkiRJUt1VeYrgusB04GcRcXdE/DQilu04QUSMi4iJETFx+vTpFb61JEmSJPW+KgusgcBI4H8zc0vgH8DRHSfIzLMzsy0z24YNG1bhW0uSJElS76uywJoKTM3MO8r/L6UouCRJkiSpX6iswMrMZ4G/R8RG5aidgAeqii9JkiRJdVd1L4KHAReWPQg+DhxUcXxJkiRJqq1KC6zMnAy0VRlTkiRJklpFpTcaliRJkqT+zAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJ/UJEDI6IOyPinoiYEhHH93ZOkqS+p+r7YEmSVFf/Aj6YmTMjYhBwS0RcnZm393ZikqS+wwJLktQvZGYCM8t/B5VD9l5GkqS+yFMEJUn9RkQMiIjJwPPAtZl5RxfTjIuIiRExcfr06T2eoySptVlgSZL6jcycnZkjgLWA0RGxeRfTnJ2ZbZnZNmzYsB7PUZLU2iywJEn9Tma+DNwA7NLLqUiS+hgLLElSvxARwyJixfLxEGBn4KFeTUqS1OfYyYUkqb9YAzg/IgZQ7GC8JDN/18s5SZL6GAssSVK/kJn3Alv2dh6SpL7NUwQlSZIkqSIWWJIkSZJUEQssSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUEQssSZIkSaqIBZYkSZIkVWRglcEiYgAwEXg6Mz9aZexmDD/69wt9/omTduuhTCRJkiT1J1UfwToCeLDimJIkSZLUEiorsCJiLWA34KdVxZQkSZKkVlLlEazTgK8AcxY0QUSMi4iJETFx+vTpFb61JEmSJPW+SgqsiPgo8Hxm3rWw6TLz7Mxsy8y2YcOGVfHWkiRJklQbVR3B2gbYPSKeAH4FfDAiflFRbEmSJElqCZUUWJl5TGaulZnDgU8B12fmZ6qILUmSJEmtwvtgSZIkSVJFKr0PFkBm3gjcWHVcSZIkSao7j2BJkvqFiFg7Im6IiAciYkpEHNHbOUmS+p7Kj2BJklRTs4AvZeakiFgeuCsirs3MB3o7MUlS3+ERLElSv5CZ0zJzUvn4NeBBYM3ezUqS1NdYYEmS+p2IGA5sCdzRxXPjImJiREycPn16j+cmSWptFliSpH4lIpYDfgMcmZmvdn4+M8/OzLbMbBs2bFjPJyhJamkWWJKkfiMiBlEUVxdm5mW9nY8kqe+xwJIk9QsREcA5wIOZeUpv5yNJ6psssCRJ/cU2wP7AByNicjl8pLeTkiT1LXbTLknqFzLzFiB6Ow9JUt/mESxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklSRygqsiFg7Im6IiAciYkpEHFFVbEmSJElqBQMrjDUL+FJmToqI5YG7IuLazHygwveQJEmSpNqq7AhWZk7LzEnl49eAB4E1q4ovSZIkSXVX5RGsdhExHNgSuKPT+HHAOIB11lmnO95akiT1kuFH/36R0zxx0m49kIkk9Z7KO7mIiOWA3wBHZuarHZ/LzLMzsy0z24YNG1b1W0uSJElSr6q0wIqIQRTF1YWZeVmVsSVJWlIRcW5EPB8R9/d2LpKkvqnKXgQDOAd4MDNPqSquJEkVOg/YpbeTkCT1XVUewdoG2B/4YERMLoePVBhfkqQlkpk3Ay/2dh6SpL6rsk4uMvMWIKqKJ0lSb7BDJknSkqi8kwtJklqZHTJJkpaEBZYkSZIkVcQCS5IkSZIq0i03GpYkqY4i4iJgB2DViJgKfCMzz+npPLwhryT1XRZYkqR+IzP37e0cpCVRt+K8bvlIdeApgpIkSZJUEQssSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFvA+WJEktynsQSVL9eARLkiRJkipigSVJkiRJFbHAkiRJkqSKeA2WJEmSWp7XJKouPIIlSZIkSRXxCNYCVLUXxL0pkiQ1Z1Ftp+2mpDqzwGoRFnySpO5i2yBJ1fEUQUmSJEmqiEewtFjc2ylJUuuyHW8drqvWY4GlXuWXhiRJqhN/m2hJWWCpT/CCaEmSpO5nAbpoFlhSyS8MSZIkLalKC6yI2AX4ITAA+GlmnlRlfKkV1K3Hx7rFkXqT7ZR6i9+h/Y/rfOGqOPuorsu4sgIrIgYAZwI7A1OBCRFxVWY+UNV7SGp9VXwZ1vULVfVmOyVJ6glVHsEaDTyamY8DRMSvgD0AGy5JtVS3o3vG6Xa2U5KkblflfbDWBP7e4f+p5ThJkurAdkqS1O0iM6sJFLE3sEtmHlL+vz+wVWaO7zDNOGBc+e9GwMOLCLsqMKOC9IzTGrkYp2fi1CkX4/RMnDrl0micd2bmsAreq10j7VQ5vpXbqjrlYpyeiVOnXIzTM3HqlEt/j9N1W5WZlQzA+4FrOvx/DHDMEsacWFFuxmmBXIzjOjeO67w7h+5op+q2XOqUi3Fc58ZxnffXOFWeIjgB2CAi1o2IpYFPAVdVGF+SpCVhOyVJ6naVdXKRmbMiYjxwDUX3t+dm5pSq4kuStCRspyRJPaHS+2Bl5h+AP1QY8mzjdHucOuVinJ6JU6dcjNMzceqUS5VxmtYN7RTUa7nUKRfj9EycOuVinJ6JU6dcjNOFyjq5kCRJkqT+rsprsCRJkiSpX7PAkiRJkqSKWGBJkiRJUkUssCRJkiSpIpX2IrgkIiKAfYAELgU+COwBPASclZlzmoj1YeDjwJrlqKeBKzPzj03EWAd4PjPfKHMbC4wEHgB+kpmzGohxCvCbzLy10fddSKzVATLz2YgYBnwAeLiZLoYjYivgwcx8NSKGAEfz1jx9OzNf6ck4Zay3AcMy87FO49+dmfc2GGPu/WyeyczrIuLTwNbAg8DZmflmo/l0irstMBq4PzP/1NMxImIFYBfm/Rxfk5kvL04uHeJekJkHNPmaPYGbMvPF8vP3A2BLinX+pcyc2mCcVTNzRof/P0O5fCi2q4Z63alqnUfE4cDlmfn3RqZfSJwl/s4p42xM8b3XMc5VmflgL8ZZE7gjM2d2GL9Lo/NW5fdFb+uL7VQZx7Zq4XH6dDu1pHH6WlvVV9upMpZtVdevr7ydqk0vghHxY2A1YGngVWAZihtA7gY8l5lHNBjnNGBD4AJg7oa0FnAA8Ncm4twPjM7Mf0bEd4H1gSsoGlQy8+AGYkwHngSGARcDF2Xm3Y28f6c4/0mxsgP4LkUjej+wLfC9zDynwThTgPeU94I5G/gnxY+Encrx/9HDcT4BnAY8DwwCxmbmhPK5SZk5ssE4F1LsLBgKvAwsB1xW5hOZeWCDce7MzNHl40OBzwOXA2OA32bmST0Ro3ztAcA3gD9RfOFA8TneGTg+My9oME7nm6gGsCNwPUBm7t5gnAcyc9Py8cXA7cCvgQ8B+2Xmzg3GaV+vEXEsxY+vXwIfBaZm5hcajFPVOn8F+AfwGHAR8OvMnN7IazvEOI1qvnOOAvYFftUpzqeAXzXx2akqzuEUn98HgRHAEZl5ZflcM9tnJd8XddAX26kyjm3VgmP0uXaq4jh9rq3qi+1UGec0bKsWFKP6diozazEA95V/BwEvAEuX/w8E7m0iziMLGB8UH6BG4zzQ4fFdwFId/r+nwRh3l383BP4fMIViT+c3gA2bWTYUG+gqwExg9XL8SsDkJuI82OHxpE7P9UacycAa5ePR5bLZs+OyazDOvR0+K88BAzqs82Y+O3d3eDyBYo8lwLJzP589EaOc/mFgxS7Gr7Sgz/gC4kwCfgHsAGxf/p1WPt6+mXw6PL5rCdb53Z1yW7Z8PKjJ5VPZOqc4VXoMcA4wHfgjcCCwfIMxqvrOeQQY1MX4pXspzn3AcuXj4cBEioar2e2zku+LOgz0wXaq4/rEtqqrGJPpY+1UxXH6XFtFH2ynyji2VQuOUXk7VadrsGYBZHHIdEJm/rv8fxbQ8GkXwBsRMaqL8aOAN5qI8/eI+GD5+AlgbYCIWKWJGAmQmY9k5jczczPgE8BgmrvR5ZuZ+c/MfAF4LDOfLeO+NPc9GnR/RBxUPr4nItoAImJDoJnTE6qKMyAzpwFk5p0Ue6uOLfdGNDNfS5WH4penaNxXKMcvQ/GF2Eyclcp1HFnuIcrMf1B+PnsoBhRfeF0tgznlc41qo/jh9TXglcy8EXg9M2/KzJuaiHNjRJxQHjq/sTwNg4jYEWjm0PmQiNgyIt5Lsf7/Ae3b/ewm4lS1zjMz52TmnzLzs8A7gB9TnO7yeIMxqvrOmVO+f2dr0Nx3YFVxlsryVIvMfILiB8+u5elkzXwGq/q+qIO+2E6BbdXC9MV2qso4fbGt6ovtFNhWLUzl7VRtrsECno2I5TJzZmbuMndkFOdz/7uJOGOB/42I5XnrkOPaFBvW2CbiHAJcEBHHla+dHBGTgRWBLzYYY74Vm8X52vcCxzSRS0bEoHLj3q09eMRgmuuo5BDgh+Uh7xnA/0XE34G/l8/1dJzXImL9LM9rz8xpEbEDxSkumzUR5xyKvYoDKL6cfx0RjwPvozj03KgVKL7gg2KZr1HmtByNb6RVxAA4EZgUEX+iWK4A61CcdvHNRoNkcU3IqRHx6/Lvcyzedj+eYtk+XP7/hYj4B/BbYP8m4kwDTikfv9hh+axCc416Vet8nnVSbmNXAVdFxNAGYxwE/LiC75wjgT9HxF+Zd52/i2L593Sc5yJiRGZOBsjMmRHxUeBcYIsm4lT1fVEHfbGdAtuqhemL7VSVcfpiW9UX2ymwrVqYytup2lyDtSARsSzF4dnnm3zd6nS4aG7unrTFeP9NKE6bGEjxgZyQDV7IPLchXpz37RRnHWBadrogMiLWBDbJzOuajPc2YF3KecrM5xYzryWKExHvAf6ZmX/tNH4Q8InMvLCJWO8AyMxnImJFivOtnyr3OC6R8gvs7Zn5t56MERErAR9m/guHX1qCPHYDtsnMry5BjBWAgeVe6kpExABgmcz8ZxOvWeJ1HhEbZuYjzea7gFhL/J0TEUtRnIbUcZ1PyMxm9ppWEici1gJmdTUfEbFNNtkhQlXfO3XUyu1U+XrbqgW/tt+0U4sbp7+0VX2hnSrj2VYtOFZl7VTtCqwOe786jpunR5dFvH5pitMUsvx/R4qeQKZk872kDCxP/aDcq7Mx8HhmvthknDaKvQSzKc6BfaiZ11ecSzD/h/nOXMIPQkSs3GwuVcYpG/ZXM/PliBhOcbrBQ5l5/5LmVMZfrB8gEfEu4D0U5/c+sIQ57J6ZnS8EXpw4/52ZP17SOB3ibdzsZ3pJt/Ny+qWg2PNZbvebA0808/mJiBVzyXu6argnscWIXdV2tVjrvIplvIj4lfyw72l9sZ0qX29b1Y0x6tpOla+1rZp/+j7TTpVxbKsWL6fF265yMS7c6o6B4rzmqRSH5v4EDO/w3KQm4twDrFQ+/h/gNuBY4FrgO03EGUtxEfMjwK4U57n+meJw4b4Nxtie4mK764CXgN8BtwI3Amv3ZC5lnDHAo8DVwE/L4Y/luDFNxDm2w+NNy7z+RnENwFZNxNmGoteXKcBW5Tp6rJyv9zcR5+jy/R+iOJT7EMWh+SnAFyv6fD7V4HQ3AKuWj/cvl81PKS7CPKyJ9/uPLoZn5z5uIs4XuxhmzH3ck8umnLaq7fzjFBcNT6Po4vWOcpuYCnysiTizyu3zs3RxoXaDMWYDf6U4HWbTJViOVW1XlazzqpZxVZ+dOgwVfn5r006VcWyrFhyjz7VT5bS2VQuers+1U2Uc26pu/NzM97oqPrRVDBS92GxWPt67/BC8r/z/7ibi3N/h8URgSPm42V6e7gNWpThU+Cqwfjn+7Y3Goej5ZW7PPOtS3McAinOT/9STuZTTP9jxi6LD+HXp0INKA3EmdXj8e2DX8vFo4LYm4txJcX7s+8uNatty/Ejg1ibiTAGGUPRc9Rrz9oZ0fxNxutrYvwh8CXhxMT5/E4BVysdDm1xXb1L8yDkX+Fk5vFb+PbeJOK9RdLv8dYoewb5B8QPqG8A3mohz+gKGH1HskW00TlXb+d3A6h22iY3K8e8EJjYR5z6KrncvpPhheCVFF7FDmsxlc4prER6l+PF8dFfb2iLiVLVdVbXOq1rGS7xd1WWo8PNbm3aqw7q2reo6Rp9rp7r4DNpWzRujz7VTHfKxreo6RuXtVJ16EVw6yxsRZualFBXp+RHxcZrrqefViNi8fDyDohckKBquZuZ3dmbOyOI85Jn51gWuzZyPOSDfulfBUxQrm8y8lrdOe+ipXOCt8/M7e5rmerTp6B2ZeXWZz50UDUijBmXmfZn5f8D0zLyljDOpyTizM/N1ivtMvE7xBUSWPf804dsUXcsu32lYjsY/O2+W1xtA0U3x3Bz+RXGha6O2plgGEzLzoMw8CJhRPm7o3jalzShyXxb4fmYeD7yUmceXjxt1EMX9bO7qNEykuYv7q9rOycxny23iqcx8uBz3JM1t529m5u8ycz+K+29cSNF72tSI+GXjqeT9mfm1zHwXcCjFvZJuiYjbmsiloyXZrqpa51Ut4yq2q7roi+0U2FYtTF9sp8C2amH6YjtVvrVt1QJU3k7VqRfBNyNi9XyrW9cpEbETxZ6R9ZuI8zngwoi4h+LGgBMj4maKPVDfbiLOUxHxHYoF/FBE/IDiBnEfojgM2YiJEXEOxY3ydqc43WLuRaTNfIFVkQsUe5gmRMSveKvHlrUp9oQ0dAPI0npR3BQwgLUiYmi+ddFnU93NdnjcuaeqpZuIM6n8klmW4rDw+RHxR4qbbTZzLvkk4IrMvKvzExHRaC8yXwD+FBG/odhjeX1EXENxo82fNZpIZk6IiJ2BwyLiBuAomvxiL+M8BewTEXsA10bEqc3GKE2g2OM535dwFD2YNaqq7ZyIWCqLC/kP7jBuAM19dtp7Zyp//FwCXFJeHP3xZmOUce4E7oyILwHbNZFLJdtVheu8qmVcxXZVF32xnQLbqoXpi+0U2FYtTF9sp+aJU8ayrXpL5e1UbTq5iIgPUewduqfT+BWA8Zl5YhOxBlCcw92xV6VrsomLBKPoSeTzFF8UZ1Dcb2Asxd69b2Z5X4xFxBhEsYdgU4pDsedm5uwo7s+wWlld90guHWJtQnGOascLh6/KJi5qjYjtO42alJmvRcTbgb0z88wG4+wOXJedeuSJiPWBvTLzew3GGQjsQ7F8LqU4TP1piuVzZqN7CCNiI4pDwfPdIT0i3t7oXtjyM/tp5v38XZmLccF4GW9N4FSgLTPXW5wYZZxlgeMozpFu5suUiFgZeKPzulqMHCrZzqO4l8d9mflGp/HDKU7h+UWDcb6cmSc3lPyCY3w6M5vZi7igOJVsV51iLsk6r2oZV7Jd1UFfbKfKOLZVC47RJ9upcnrbqq5j9Ll2qoxjW7XgGJW3U7UpsCRJkiSp1bXE+e8RcXVPx4mIFSLipIh4KCJejIgXIuLBctyKfSmX3sqnL85X3ZZNREyKiGPLva3doje2z6riuJ13f5yemK866KvLtr/n0xfnqaZxWqKtcjuvXz5VxOmOearNNVgRMXJBTwEjejoOxTmu1wM7zD0PN4qbs40tnxvTSrn0UD4HNpNPVXFqNl+1WjYUF22uCNwQEc8CFwEXZ+YzDb6e8r1rtX3WaZ3XKZcq86nbfNVBCy3bsdTrO7Ru+fT4d3rN5qmOcWrTVrXQdt7Sn8GaLefK26nanCIYEbOBm2Dei/BK78vMhnomqTDOw5m5UbPP1TWXmubT5+arhstmUmaOLB9/ANiX4v4kDwIXZebZDcap2/ZZp3Vem1wqzqdW81UHfXXZmk/3xiinrc081TRObdoqt/OWy6eKdV59O5WL0bd7dwwU3WpusIDn/t4Lcf4EfAV4e4dxb6foIee6Vsulpvn0ufmq4bKZ76aIFL2C7QL8rBeWcd3iuJ23yHzVYeiry9Z8+tc81TRObdoqt/OWy6eKdV55O1Wna7COY8HXhB3WC3E+SXFDwJvK8zFfoui6dmWK+w+0Wi51zKcvzlfdls0jnUdk5uzM/GMW9ytp1HHUa/usKo7beffHqWq+6uA4+uayNZ/ujQH1mqc6xqlTW1VFjCrj9NXPYFX5VBGn8naqNqcISpIkSVKrq80RrIhYOSK+HhGHROFrEfG7iPh+RKzUi3E+W8b5arNx6pRLzfPpM/PVD5ZNX43jdl7z+aqDvrpszad/zVPN4/T6tlWnXDrFqdsyrls+Vazzytqp2hRYwC8o7nD+XuAGYHXgu8DrwHm9GKetjLPGYsSpUy51zqcvzVdfXzZ9NY7beffHWdL5qoO+umzNp3tjdIxTh3mqc5w6bFt1yqVjnLot47rlU8U6r66dWpwLt7pjACaXfwN4uqvnWi1OnXIxjuvcOK7zVo5Th6Fuy8Q43R+nTrkYx3VunHqv845DnY5gLVUehlsbWC4ihgNExCrA0i0ap065GKdn4tQpF+P0TJw65dKX49RB3ZaJcbo/Tp1yMU7PxKlTLsbpmTjVt1OLU5V1x0Bxz4PnymEv4DrgWuBpYFwrxqlTLsZxnRvHdd7Kceow1G2ZGKf749QpF+O4zo1T73XecahVL4IRMQCIzJwVEQMp7sD8dGZOa9U4dcrFOD0Tp065GKdn4tQpl74cpw7qtkyM0/1x6pSLcXomTp1yMU7PxKm6nRq4OC/qLpk5OyLaImJtYDbwyOLMWJ3i1CkX4/RMnDrlYpyeiVOnXPpynDqo2zIxTvfHqVMuxumZOHXKxTg9E6fqdqo2R7AiYnvgB8DLFD2B3AqsBLwJ7J+Zf2+1OHXKxTg9E6dOuRinZ+LUKZe+HKcO6rZMjNP9ceqUi3F6Jk6dcjFOz8TplnYqF+O8wu4YgLuBYeXjdYHLy8c7A39qxTh1ysU4rnPjuM5bOU4dhrotE+N0f5w65WIc17lx6r3OOw516kVwQGZOLx8/BbwTIDOvBdZs0Th1ysU4PROnTrkYp2fi1CmXvhynDuq2TIzT/XHqlItxeiZOnXIxTs/EqbydqtM1WBMj4hzgemB34EaAiBgKDGjROHXKxTg9E6dOuRinZ+LUKZe+HKcO6rZMjNP9ceqUi3F6Jk6dcjFOz8SpvJ2q0zVYg4BDgU2Be4Bzs7jgbAiwWmY+2Wpx6pSLcXomTp1yMU7PxKlTLn05Th3UbZkYp/vj1CkX4/RMnDrlYpyeidMd7VRtCixJkiRJanW1uQYrIpaLiBMi4v6IeCUipkfE7RExtlXj1CkX4/RMnDrlYpyeiVOnXPpynDqo2zIxTvfHqVMuxumZOHXKxTg9E6eqXOaJWZcjWBFxJXA5xd2TPwEsC/wKOJbiRl9fbbU4dcrFOD0Tp065GKdn4tQpl74cpw7qtkyM0/1x6pSLcXomTp1yMU7PxOmWdioXo+vB7hiAezr9P6H8uxTwUCvGqVMuxnGdG8d13spx6jDUbZkYp/vj1CkX47jOjVPvdd5xqM0pgsA/ImJbgIjYHXgRIDPnANGiceqUi3F6Jk6dcjFOz8SpUy59OU4d1G2ZGKf749QpF+P0TJw65WKcnolTfTu1OFVZdwzAe4A7gZeAW4CNyvHDgMNbMU6dcjGO69w4rvNWjlOHoW7LxDjdH6dOuRjHdW6ceq/zeWIuzou6awA2AT4ELNdp/C6tGqdOuRjHdW4c13krx6nDULdlYpzuj1OnXIzjOjdOvdd5++sW50XdMQCHAw8BVwBPAHt0eG5SK8apUy7GcZ0bx3XeynHqMNRtmRin++PUKRfjuM6NU+91Pk/MxXlRdwzAfZRVIzAcmAgcUf5/dyvGqVMuxnGdG8d13spx6jDUbZkYp/vj1CkX47jOjVPvdd5xGEh9LJWZMwEy84mI2AG4NCLeSXMXmNUpTp1yMU7PxKlTLsbpmTh1yqUvx6mDui0T43R/nDrlYpyeiVOnXIzTM3Eqb6fq1IvgcxExYu4/5Yx+FFgV2KJF49QpF+P0TJw65WKcnolTp1z6cpw6qNsyMU73x6lTLsbpmTh1ysU4PROn+nYqF+OwV3cMwFrA6gt4bptWjFOnXIzjOjeO67yV49RhqNsyMU73x6lTLsZxnRun3uu84xDliyVJkiRJS6hOpwhKkiRJUkuzwJIkSZKkilhgqV+JiNkRMbnDMLzJ138uIg4oH58XEXuXj38aEZuWj7/aQJzhEfF6mcMDEXFBRAxajFlqJvcREfGRDv+PjYgzysfHRcTTZT5/jYjL5s6PJKln2VbZVqm1WWCpv3k9M0d0GJ5o5sWZeVZmXtDF+EMy84Hy30U2WqXHMnMERQ81awGfaCaXZkTEQGAE8JGFTHZquUw2AC4Gro+IYd2VkyRpgWyrFsy2SrU3sLcTkHpTRCwHXAmsBAwCjs3MK8vnDgC+DCRwb2buHxHHATMz8+ROcW4sp90bGBIRk4EpwGPAi5l5WjndicDz5XsCkJmzI+JOYM1ymvcCpwDLATOAsZk5rXyPe4DtKbbdgzPzzohYGTgXWA/4JzAuM+8tc12/HP8UsE2Z27bAdxa2XDLz4ojYDfg08MMGF6ckqRvYVnXNtkp1ZYGl/mZugwLwN2AfYM/MfDUiVgVuj4irgE2BY4GtM3NG2TAsUmYeHRHjy719lKd1XAacFhFLAZ8CRgPLz31NRAwGtgKOKE+9+BGwR2ZOj4hPAicCB5eTD83MERGxHUVDtTlwPMWdxj8eER8ELqDYA0g5H9tm5usRMRZoy8zx5fuOXcTsTAI2bmS+JUmVsq2yrVILs8BSf/P63AYFoGwkvl02AnMo9sy9Hfgg8OvMnAGQmS8uzptlcUfwFyJiyzLu3Zn5QkQsD6xfNqDrAr8v9+RtTtEQXRsRAAOAaR1CXlTGvTki3hYRKwLbAnuV46+PiFUi4m3l9Fdl5uuLkzuLefdySdISs61qnG2VascCS/3dfsAw4L2Z+WZEPAEMrvg9fgqMBVan2JM312PlHr5VgVsjYneKPZVTMvP9C4jV+cZ1i7qR3T8WI9+5tgQmLsHrJUnVsK1aMNsq1Y6dXKi/WwF4vmywdgTeWY6/HtgnIlYBaPS0i9KbnXpZuhzYBRgFXNN54nLP49HAMcDDwLCIeH/5voMiYrMOk3+yHL8t8EpmvgL8haLxJSJ2AGZk5qtd5PUaHU73WJiI2AsYQ7kXUpLUq2yrumBbpbqywFJ/dyHQFhH3AQcADwFk5hSK88lvioh7KC7kbdTZwL0RcWEZ69/ADcAlmTl7Aa+5AhhKcX773sB3y/edDGzdYbo3IuJu4Czgs+W444D3RsS9wEnAgQt4jxuATcvubT/ZxfNfmNv1LfAZ4IOZOb2B+ZUkdS/bqrfYVqn2InNRR20lLYnyguFJwD6Z+dcliHMj8OXM9FQISVKlbKuk6ngES+pGUdwA8VHgz0vSYEmS1F1sq6RqeQRLkiRJkiriESxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqsjA3nrjVVddNYcPH95bby9JqrG77rprRmYO6+08bKskSQuyoLaq1wqs4cOHM3HixN56e0lSjUXEk72dA9hWSZIWbEFtlacISpIkSVJFLLAkSZIkqSIWWJIkSZJUkV67BktS93rzzTeZOnUqb7zxRm+nIi3Q4MGDWWuttRg0aFBvp9Iwty21glbctqS+wgJL6qOmTp3K8ssvz/Dhw4mI3k5Hmk9m8sILLzB16lTWXXfd3k6nYW5bqrtW3bakvsJTBKU+6o033mCVVVbxB6BqKyJYZZVVWu5IkNuW6q5Vty2pr7DAkvowfwCq7lr1M9qqeav/8DMq9R4LLEmSJEmqiNdgSTU1/OjfL3KaJ07ardJ4zWjkvSOCL37xi/zgBz8A4OSTT2bmzJkcd9xxnHXWWQwdOpQDDjiA8847jzFjxvCOd7yj0hx7wsSJE7ngggs4/fTT53tu7k1qV1111SV6j9NOO42VV16ZAw44gLFjx3Lttdfy+OOPs8wyyzBjxgza2tp44oknAJgyZQqHHXYYTz/9NHPmzOGAAw7g2GOPXeTe7GnTpnHggQdy9tlns+6663L66adz2GGHATB+/Hja2toYO3YsmcmJJ57I+eefT0Sw5pprcsYZZ7DZZpvxz3/+k3322YfHHnuMAQMG8LGPfYyTTjoJgDPOOIOhQ4dy8MEHL9GyqCO3re7htuW2pf6r6u/VRWnm91QjPIIlqdsss8wyXHbZZcyYMWO+5z73uc9xwAEHAHDeeefxzDPPdBlj9uzZ3Zrjkmpra+vyB2BVZs2axbnnnsunP/3p9nEDBgzg3HPPnW/a119/nd13352jjz6ahx9+mHvuuYfbbruNH//4x4t8nz/+8Y98+MMfBmC11Vbjhz/8If/+97/nm+7MM8/ktttu45577uGRRx7hmGOOYffdd2+/1uPLX/4yDz30EHfffTe33norV199NQAHH3wwP/rRjxZrGWh+bltLzm1LUnexwJLUbQYOHMi4ceM49dRT53vuuOOO4+STT+bSSy9l4sSJ7LfffowYMYLXX3+d4cOHc9RRRzFy5Eh+/etfc9FFF7HFFluw+eabc9RRR7XHOOecc9hwww0ZPXo0hx56KOPHjwdg+vTp7LXXXowaNYpRo0Zx6623tr/nwQcfzA477MB66623wB9vXb3f5Zdfzk477URmMm3aNDbccEOeffZZbrzxRj760Y8C8MILLzBmzBg222wzDjnkEDKzPeYvfvELRo8ezYgRI/jP//xPZs+ezezZsxk7diybb745W2yxRZfL6frrr2fkyJEMHPjWCQdHHnkkp556KrNmzZpn2l/+8pdss802jBkzBoChQ4dyxhlntO/pPuKIIzjhhBMAuOaaa9huu+2YM2cOUPwI3HXXXQEYNmwYO+20E+eff/58+Xz3u99t32MOMGbMGLbeemsuvPBChg4dyo477gjA0ksvzciRI5k6dWp7LsOHD+fOO+/scpmrOW5bbltuW1J9WWBJ6laf//znufDCC3nllVe6fH7vvfemra2NCy+8kMmTJzNkyBAAVlllFSZNmsR2223HUUcdxfXXX8/kyZOZMGECV1xxBc888wzf/OY3uf3227n11lt56KGH2mMeccQRfOELX2DChAn85je/4ZBDDml/7qGHHuKaa67hzjvv5Pjjj+fNN9+cJ59nnnmmy/fbc889WWONNTjzzDM59NBDOf7441l99dXnee3xxx/Ptttuy5QpU9hzzz156qmnAHjwwQe5+OKLufXWW5k8eTIDBgxon9+nn36a+++/n/vuu4+DDjpovuVz66238t73vneeceussw7bbrstP//5z+cZP2XKlPmmXX/99Zk5cyavvvoq3/nOd7j44ou54YYbOPzww/nZz37GUkstxezZs3n44YfZdNNN21931FFHcfLJJ89zlOPVV1/lH//4B+utt94879HW1saUKVPmGffyyy/z29/+lp122mme6f7yl7/MN49aPG5bblsdp3PbkurDa7Akdau3ve1tHHDAAZx++untP/Aa8clPfhKACRMmsMMOOzBs2DAA9ttvP26++WYAtt9+e1ZeeWUA9tlnHx555BEArrvuOh544IH2WK+++iozZ84EYLfddmOZZZZhmWWWYbXVVuO5555jrbXWap92Qe/38Y9/nB/96EdsvvnmvO9972PfffedL+ebb76Zyy67rP19VlppJQD+/Oc/c9dddzFq1CigON1otdVW42Mf+xiPP/44hx12GLvttlv73vGOpk2bxiabbDLf+GOOOYY99tiD3XZr/LzxoUOH8pOf/ITtttuOU089lfXXXx+AO+64g6222mqeaddbbz222morfvnLXzYcf65Zs2ax7777cvjhh8/zg3G11Vab58e6lozbltvWXG5bUr0sssCKiMHAzcAy5fSXZuY3Ok0zFvg+8HQ56ozM/Gm1qUpqVUceeSQjR47sci/ygiy77LKL/X5z5szh9ttvZ/DgwfM9t8wyy7Q/HjBgwHynAi3M1KlTWWqppXjuueeYM2cOSy3V2EkAmcmBBx7Id77znfmeu+eee7jmmms466yzuOSSS+a7/mPIkCFd3stmgw02YMSIEVxyySXt4zbddNP2H8hzPf744yy33HK87W1vA+C+++5jlVVWmee6nKuvvppddtllvvf46le/yt577832228PFD/ol112WR5//PF5ftzddddd7dMAjBs3jg022IAjjzxynnhvvPFGU4VAM/prW+W25bYF3bttSWpeI99g/wI+mJnvAUYAu0TE+7qY7uLMHFEOLd1gSarWyiuvzCc+8QnOOeecLp9ffvnlee2117p8bvTo0dx0003MmDGD2bNnc9FFF7H99tszatQobrrpJl566SVmzZrFb37zm/bXjBkzZp6LvidPntxwrgt6v1mzZnHwwQdz0UUXsckmm3DKKafM99rtttuufa/01VdfzUsvvQTATjvtxKWXXsrzzz8PwIsvvsiTTz7JjBkzmDNnDnvttRff+ta3mDRp0nwxN9lkEx599NEuc/3a177GySef3P7/fvvtxy233MJ1110HFHvzDz/8cL7yla8A8OSTT/KDH/yAu+++m6uvvpo77rgDKI4CfOhDH5ov/sYbb8ymm27Kb3/72/Zx//M//8Phhx/O66+/DhRHNG655Zb2jgKOPfZYXnnlFU477bT54j3yyCNsvvnmXc5LBfplW+W25bYF3b5tSWrSIo9gZXEl6czy30HlkAt+haQ6qroL0mZ96Utf4owzzujyubFjx/K5z32OIUOG8H//93/zPLfGGmtw0kknseOOO5KZ7Lbbbuyxxx5AsRd49OjRrLzyymy88cassMIKAJx++ul8/vOf593vfjezZs1iu+2246yzzmoozwW93wknnMAHPvABtt12W97znvcwatSo+U4h+sY3vsG+++7LZpttxtZbb80666wDFHu/v/WtbzFmzBjmzJnDoEGDOPPMMxkyZAgHHXRQ+8XwXe2F33XXXdl///27zHWzzTZj5MiR7T8ehwwZwpVXXslhhx3G5z//eWbPns3+++/P+PHjyUw++9nPcvLJJ/OOd7yDc845h7Fjx/KHP/yBwYMHs/zyy3f5Hl/72tfYcsst2/8/7LDDeOmll9hiiy0YMGAAq6++OldeeSVDhgxh6tSpnHjiiWy88caMHDkSKLqhnnudzq233spxxx3X0HpoVm+2VW5bblt9eduS1Lzo2BPPAieKGADcBbwLODMzj+r0/FjgO8B04BHgC5n59y7ijAPGAayzzjrvffLJJ5c0f6nPWtL7YD344INdXl/Ql8ycOZPllluOWbNmseeee3LwwQez55579nZaldtzzz353ve+xwYbbFB57F/84hdMnTqVo48+uvLYHd19992ccsop83UeAF1/ViPirsxsa+Y9eqqtctvqO/rjtiW1gla5D9aC2qqGTnLOzNmZOQJYCxgdEZ2PQ/8WGJ6Z7wauBebvf7SIc3ZmtmVm29yLXCVpcR133HGMGDGCzTffnHXXXZePf/zjvZ1StzjppJOYNm1at8T+zGc+0+0/AAFmzJjBN7/5zW59D9uq6rhtLbm+tG1Jak5TvQhm5ssRcQOwC3B/h/EvdJjsp8D3qklPkhas4zUSfdlGG23ERhtt1NtpLJGdd965x97LtmrJuW21jp7ctiQ1ZpFHsCJiWESsWD4eAuwMPNRpmjU6/Ls78GCFOUpaTI2cAiz1pqo+oz3dVrltqe78jEq9p5EjWGsA55fnti8FXJKZv4uIE4CJmXkVcHhE7A7MAl4ExnZXwpIaM3jwYF544QVWWWUVIqK305Hmk5m88MILXXb5vRh6rK1y21LdVbxtSWpSI70I3gts2cX4r3d4fAxwTLWpSVoSa621FlOnTmX69Om9nYq0QIMHD57nZrSLqyfbKrcttYKqti1JzWvqGixJrWPQoEGsu+66vZ2G1Oe4bUmSFqaxW6VLkiRJkhbJAkuSJEmSKmKBJUmSJEkVscCSJEmSpIpYYEmSJElSRSywJEmSJKkiFliSJEmSVBELLEmSJEmqiAWWJEmSJFXEAkuSJEmSKmKBJUmSJEkVscCSJEmSpIpYYEmSJElSRSywJEmSJKkiFliSJEmSVBELLEmSJEmqyMDeTkCqi+FH/36R0zxx0m49kIkkSZJa1SKPYEXE4Ii4MyLuiYgpEXF8F9MsExEXR8SjEXFHRAzvlmwlSeqCbZUkqS4aOUXwX8AHM/M9wAhgl4h4X6dpPgu8lJnvAk4FvltplpIkLZxtlSSpFhZZYGVhZvnvoHLITpPtAZxfPr4U2CkiorIsJUlaCNsqSVJdNNTJRUQMiIjJwPPAtZl5R6dJ1gT+DpCZs4BXgFW6iDMuIiZGxMTp06cvUeKSJHVkWyVJqoOGCqzMnJ2ZI4C1gNERsfnivFlmnp2ZbZnZNmzYsMUJIUlSl2yrJEl10FQ37Zn5MnADsEunp54G1gaIiIHACsALFeQnSVJTbKskSb2pkV4Eh0XEiuXjIcDOwEOdJrsKOLB8vDdwfWZ2PvddkqRuYVslSaqLRu6DtQZwfkQMoCjILsnM30XECcDEzLwKOAf4eUQ8CrwIfKrbMpYkaX62VZKkWlhkgZWZ9wJbdjH+6x0evwHsU21qkiQ1xrZKklQXTV2DJUmSJElaMAssSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUEQssSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUEQssSZIkSarIIgusiFg7Im6IiAciYkpEHNHFNDtExCsRMbkcvt496UqSND/bKklSXQxsYJpZwJcyc1JELA/cFRHXZuYDnab7S2Z+tPoUJUlaJNsqSVItLPIIVmZOy8xJ5ePXgAeBNbs7MUmSGmVbJUmqi6auwYqI4cCWwB1dPP3+iLgnIq6OiM0W8PpxETExIiZOnz69+WwlSVoE2ypJUm9quMCKiOWA3wBHZuarnZ6eBLwzM98D/Ai4oqsYmXl2ZrZlZtuwYcMWM2VJkrpmWyVJ6m0NFVgRMYiiwbowMy/r/HxmvpqZM8vHfwAGRcSqlWYqSdJC2FZJkuqgkV4EAzgHeDAzT1nANKuX0xERo8u4L1SZqCRJC2JbJUmqi0Z6EdwG2B+4LyIml+O+CqwDkJlnAXsD/xURs4DXgU9lZlafriRJXbKtkiTVwiILrMy8BYhFTHMGcEZVSUmS1AzbKklSXTTVi6AkSZIkacEssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFRnY2wmoMcOP/v0ip3nipN16IBNJkiRJC+IRLEmSJEmqiAWWJEmSJFVkkQVWRKwdETdExAMRMSUijuhimoiI0yPi0Yi4NyJGdk+6kiTNz7ZKklQXjVyDNQv4UmZOiojlgbsi4trMfKDDNLsCG5TDVsD/ln8lSeoJtlWSpFpYZIGVmdOAaeXj1yLiQWBNoGOjtQdwQWYmcHtErBgRa5SvldQH9MWOVvriPPVXtlWS+otG2q4q2Q42r6lrsCJiOLAlcEenp9YE/t7h/6nluM6vHxcREyNi4vTp05tMVZKkRbOtkiT1poYLrIhYDvgNcGRmvro4b5aZZ2dmW2a2DRs2bHFCSJK0QLZVkqTe1lCBFRGDKBqsCzPzsi4meRpYu8P/a5XjJEnqEbZVkqQ6aKQXwQDOAR7MzFMWMNlVwAFlD03vA17xnHZJUk+xrZIk1UUjvQhuA+wP3BcRk8txXwXWAcjMs4A/AB8BHgX+CRxUeaaSJC2YbZUkqRYa6UXwFiAWMU0Cn68qKUmSmmFbJUmqi6Z6EZQkSZIkLZgFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUauQ+WJKkHDD/694uc5omTduuBTCRJ0uLyCJYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkilhgSZIkSVJFLLAkSZIkqSIWWJIkSZJUEQssSZIkSaqIBZYkSZIkVcQCS5IkSZIqssgCKyLOjYjnI+L+BTy/Q0S8EhGTy+Hr1acpSdKC2VZJkupiYAPTnAecAVywkGn+kpkfrSQjSZKadx62VZKkGljkEazMvBl4sQdykSRpsdhWSZLqoqprsN4fEfdExNURsVlFMSVJqpJtlSSp2zVyiuCiTALemZkzI+IjwBXABl1NGBHjgHEA66yzTgVvLUlSQ2yrJEk9YomPYGXmq5k5s3z8B2BQRKy6gGnPzsy2zGwbNmzYkr61JEkNsa2SJPWUJT6CFRGrA89lZkbEaIqi7YUlzkxSnzT86N8vcponTtqtBzJRf2JbJUnqKYsssCLiImAHYNWImAp8AxgEkJlnAXsD/xURs4DXgU9lZnZbxpIkdWJbJUmqi0UWWJm57yKeP4Oia1xJknqFbZUkqS6q6kVQkiRJkvq9KnoRlCRJHTRyrWGVvG5RkurDI1iSJEmSVBELLEmSJEmqiAWWJEmSJFXEAkuSJEmSKmKBJUmSJEkVscCSJEmSpIpYYEmSJElSRSywJEmSJKkiFliSJEmSVBELLEmSJEmqiAWWJEmSJFXEAkuSJEmSKmKBJUmSJEkVscCSJEmSpIoM7O0EJEn1NPzo3y9ymidO2q0HMpEkqXV4BEuSJEmSKrLII1gRcS7wUeD5zNy8i+cD+CHwEeCfwNjMnFR1olKrcK9/93MZL1x/XD62VZLmauQ7sEp97ftUS66RI1jnAbss5PldgQ3KYRzwv0ueliRJTTkP2ypJUg0sssDKzJuBFxcyyR7ABVm4HVgxItaoKkFJkhbFtkqSVBdVXIO1JvD3Dv9PLcdJklQXtlWSpB7Ro51cRMS4iJgYEROnT5/ek28tSVJDbKskSUuiigLraWDtDv+vVY6bT2aenZltmdk2bNiwCt5akqSG2FZJknpEFffBugoYHxG/ArYCXsnMaRXElSSpKrZVFbKXNklasEa6ab8I2AFYNSKmAt8ABgFk5lnAHyi6vX2Uouvbg7orWUmSumJbJUmqi0UWWJm57yKeT+DzlWUkSVKTbKskSXXRo51cSJIkSVJfVsU1WH1SI+eXe064JEmSpI48giVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFbHAkiRJkqSKWGBJkiRJUkUssCRJkiSpIgN7OwH1rOFH/36R0zxx0m49kIkkSeoNjfwWqJK/K9TfeARLkiRJkipigSVJkiRJFanVKYKeviZJkiSplTVUYEXELsAPgQHATzPzpE7PjwW+DzxdjjojM39aYZ4ty6JRknqGbZUkqQ4WWWBFxADgTGBnYCowISKuyswHOk16cWaO74YcJUlaKNsqSVJdNHIN1mjg0cx8PDP/DfwK2KN705IkqSm2VZKkWmikwFoT+HuH/6eW4zrbKyLujYhLI2LtrgJFxLiImBgRE6dPn74Y6UqS1CXbKklSLVTVi+BvgeGZ+W7gWuD8ribKzLMzsy0z24YNG1bRW0uS1BDbKklSt2ukwHoa6LiXby3eukAYgMx8ITP/Vf77U+C91aQnSVJDbKskSbXQSIE1AdggItaNiKWBTwFXdZwgItbo8O/uwIPVpShJ0iLZVkmSamGRvQhm5qyIGA9cQ9H17bmZOSUiTgAmZuZVwOERsTswC3gRGNuNOUuSNA/bKklSXTR0H6zM/APwh07jvt7h8THAMdWmJklS42yrJEl1UFUnF5IkSZLU71lgSZIkSVJFGjpFUKq74Uf/fqHPP3HSbj2UiSRJkvozj2BJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkipigSVJkiRJFemT98HynkiSJGlxLep3RNX8XSL1LR7BkiRJkqSKWGBJkiRJUkUssCRJkiSpIhZYkiRJklQRCyxJkiRJqogFliRJkiRVxAJLkiRJkirSUIEVEbtExMMR8WhEHN3F88tExMXl83dExPDKM5UkaSFsqyRJdbDIAisiBgBnArsCmwL7RsSmnSb7LPBSZr4LOBX4btWJSpK0ILZVkqS6aOQI1mjg0cx8PDP/DfwK2KPTNHsA55ePLwV2ioioLk1JkhbKtkqSVAuRmQufIGJvYJfMPKT8f39gq8wc32Ga+8tpppb/P1ZOM6NTrHHAuPLfjYCHF5HfqsCMRUzTCOO0Ri7G6Zk4dcrFOD0Tp065NBrnnZk5rNGAvdxWVamqZVxXzl9rc/5aW1+fP+j5eeyyrRrYgwmQmWcDZzc6fURMzMy2JX1f47RGLsbpmTh1ysU4PROnTrlUGae7NNtWVanuy2ZJOX+tzflrbX19/qA+89jIKYJPA2t3+H+tclyX00TEQGAF4IUqEpQkqQG2VZKkWmikwJoAbBAR60bE0sCngKs6TXMVcGD5eG/g+lzUuYeSJFXHtkqSVAuLPEUwM2dFxHjgGmAAcG5mTomIE4CJmXkVcA7w84h4FHiRomGrQlWnaBine2MYp7Xi1CkX4/RMnDrlUmWcdr3cVlWpV05N7EHOX2tz/lpbX58/qMk8LrKTC0mSJElSYxq60bAkSZIkadEssCRJkiSpIhZYkiSpW0XEkRExdAHPjY2IM3o6p+4SEbtHxNG9nUczImLtiLghIh6IiCkRcUQ5/riIeDoiJpfDR8rx20TEvRExMSI2KMetGBF/ioha/raMiCci4r5yPiaW41aOiGsj4q/l35XK8XuVy+EvEbFKOW79iLi4N+dhQSJiow7raHJEvFpuc31m/bUar8GSJEndKiKeANo639S5fG5s+dz4zs+pZ0TEGsAamTkpIpYH7gI+DnwCmJmZJ3ea/jLgcGA4sGdmfikiTgZ+l5k39mTujerqMxgR3wNezMyTyqJ4pcw8KiJuBD4C/Ec57kcRcRHw9cz8ay+k37CIGEBxS4qtgIPoI+uv1fTojYYXJiIC2AdI4FLgg8AewEPAWZk5p4lYH6b4YlizHPU0cGVm/rGJGOsAz2fmG2VuY4GRwAPATzJzVqOxqhARqwNk5rMRMQz4APBwZk5pIsZWwIOZ+WpEDAGO5q15+nZmvtKTccpYbwOGZeZjnca/OzPvbTDG3C6Zn8nM6yLi08DWwIPA2Zn5ZqP5dIq7LTAauD8z/9TTMSJiBWAX5v0cX5OZLy9OLh3iXpCZByxJjN5W1TqPiMOByzPz70uYzxJ/55RxNqb43usY56rMfLAX46wJ3JGZMzuM36XReavy+6LVRMRw4I/A7RSfzwnAz4DjgdWA/YApwI+AzYFBwHGZeWX52p8Dy5bhxmfmbRGxA3AcMKN8zV3AZ+rU3XxELAtcQnEvsgHAr4F3ADdExIzM3DEiDgKOAV4G7gH+1UvpNqXBdbopZcEYEecBrwJtwOrAVzLz0p7PfOEycxowrXz8WkQ8yFvfH115ExhaDm9GxPrA2i3443wPYIfy8fnAjcBRwBxgGd6avw8Az9a9uCrtBDyWmU8WP1+71HLrr8FtD+CHwGDgdeCgzHw4Ir4AbJGZB0fEFsBFwOjM/Ge3JZyZtRiAH1MUVlcBv6D4Qt4f+BXwwybinAb8geLH17bl8KlyXDNx7geGlo+/W+b2GeBciu5/G4lxCrBNBcvmP4G/AU8A/wXcQdHd8MPAZ5uIMwUYWD4+u1xW2wLfAC7rhTifAJ4BJpcxR3V4blITcS4ELgZ+S/GD5PLys3MecH4Tce7s8PjQMq9vALcCR/dUjPK1BwCPAf8LHFsOZ5XjDmgizlWdht8CM+f+30ScPYGVy8fDgAuA+8rlvlYTcVbt9P9ngNOBcZRH1Ht4nb9Sfgb/Avw3RbHf7PZ5GtV85xxVfl6OLpfLZ8rHk5v87FQV53CK75grKL579ujwXDPbZyXfF604UOwdngVsQXFK/l0UbUhQ/LC7Avg2RYEEsCLwCEVRNRQYXI7fgKKreSh+DL5CUbwsBfwfsG1vz2un+d6LYkfk3P9XKD9Dq5b/rwE8VX6XLE3x/XhGb+dd4TodO3d+yu+kX5fTbgo82tvz0OA8PgW8jaKYfwK4t5zPlcppRlD80L2h/Cz+Ctigt3NfxHz9DZhUrrNx5biXOzwfc/8Hdi6n+235+f0TZRtY96FcT+PLx31p/TWy7b2Nt9qbDwG/KR8vBdxM8VtmIhX8Nl9kvr29wDosuPvKv4OAF4Cly/8HAvc2EeeRBYwP4K9NxHmgw+O7gKU6/H9PgzGmlyvySeB7wJaLu2woGttVKH4cr16OXwmY3EScBzs8ntTpud6IM5nilAQojvI8RHGoGuDuJuLc2+Gz8hwwoMM6b+azc3eHxxMof2xT/Ni5r6dilNM/DKzYxfiVFvQZX0CcSRQ7LHYAti//Tisfb99EnI7bw8XAFyi+lMcC1zaTT4fHx1Lcs+hAih8gp/bGOqf44h1DsdNiOsUesgOB5RuMUdV3ziPAoC7GL91Lce4DlisfD6f4Ljti7nJrIk4l3xetOJTL7a8d/r8A2K98vF75HTiRYofe5HJ4CtiE4kfdz8v1MBn4Z/m6HTpucxQ7YT7T2/Paab43pPhR913gA+W4J3irwPo4cEGH6Q+ntQqsRa3TscxbYO3XYfrXenseFjF/y1H85vmP8v+3UxyFXAo4kS52MAPbAaeW6/1iijbn7b09L13kuWb5dzWKo6bb0aHAKp97qYvXHQAcCbyPYmf7Tyh3wNdtKL/nZ8xd/n1s/TWy7a1NscP1/vK786EO069H8Rv6Bz2Rb50uZJsFkMWpPRMy89/l/7MoDtU26o2IGNXF+FHAG03E+XtEfLB8/ATFSmPuxY4NmpqZbRR7Ql4DfhERD0XENyJiwybivJmZ/8zMFygO+z4LkJkvUZxS2aj7y9MyAO6JiDaAMpdmTqOrKs6ALE5LIDPvBHYEji1P22pmvpYqTxlbnqIQXaEcvwxFwd5MnJXKdRyZOb3M7R+Un88eigHFj/OulsGc8rlGtVE0ll8DXsni8P/rmXlTZt7URJwBHR6/KzNPzcypmXkexV7oRnXM/T8oGvHzgU9T7G1qVFXrPDNzTmb+KTM/S3Ea048pTs18vMEYVX3nzCnfv7M1aO47sKo4S2V5WmBmPkHxw37XiDiF5j6DVX1ftKqOp77N6fD/HIodBAHslZkjymGdLE7l/ALFzoP3UGzHSy8g5mxqdLo/QGY+QnEa6H3AtyLi672cUtUWtU4XNn0z206PiohBwG+ACzPzMoDMfC4zZ2dxmcZPKHaGdnxNUOws+ybFUemvlNMd3pO5NyIzny7/Pk/xI3w08Fx5/dnc69Ce7/iaKDpmGQucSXEq2oHALbx1Olrd7EqxI+s56Fvrr7Sobe+bwA2ZuTnwMYpTBefagKLA6qp9rFydvpSfjYjlMnNmZu4yd2R57dG/m4gzFvjf8iLNqeW4tSlOqRjbRJxDgAsi4rjytZMjYjLFKRxfbDBGQntj803gmxHxbmBfitOH3tVonIgYVBafu80dGRGDaa4nyEOAH0bEsRR7OP4vIv4O/L18rqfjvBYR62d5/VVmTiuvL7gC2KyJOOdQHP0aQFFI/DoiHqfY2/SrJuKsQFGMBMUyX6PMaTkabxSriAHFnqZJEfEniuUKsA5Fsf7NRoOUX6qnRsSvy7/PsXjb/Y0RcQLwnfLxnpl5eUTsSLF9NGpIRGxJ8bkdUBaeZOabETG7iThVrfN51km5jV0FXBUL6PGsCwcBP67gO+dI4M8R8VfmXefvApq5+L+qOM9FxIjMnAyQmTMj4qMUp2Rs0UScqr4v+qprgMMi4rDMzIjYMjPvpvgumZqZcyLiQObdyVFrEfEOio4DfhERL1Os59codojMoDjN/YfljqhXKa6/vqeX0hXtP7TPoTjifEqH8WvM3RFKcXrV/Z1eegDwh8x8sfzOnFMOjX5/9ojyusClsri+bFmKsxZOoPi+PxA4qfx7ZaeX/g9wetlGDaH4XVe7+etgX4rri4C+s/6asALFNcfQof0tr2k/neJo3RkRsXd287WQtSmwMnPXBTz1GvDRJuJMArYqC7P2C7znHvVpIs7fgR0jYhOKw6bnUfx4mpCNd7gx3w/qLDpuuJfi4t5G7dnh9VM7jF8F+FKjQbK4mHxsFB1LrEux/qfO3dPR03Eoriebp0Asv/x2obg+q9F8To2y69TMfCYiLqA4GvKT8shYo3GGL+CpOXRYB90do4xzfkRcBXyYtz7HNwLHlEcum1J+bvaJiN0oftA0azxFIfNw+f8XIuIfFOen799EnGkU1yYCvNihAF2FJo7wVbXOgU8u5D0auvg1M++imu+cP5ZHdUYzb+cUEzKz4eKzqjgUDe8866Q8o+CAiPj/msinqu+LvuqbFNel3RtF98h/o2jzfgz8JiIOoDht9R+9lmHztgC+HxFzKI5S/hfwfuCPEfFMFp1cHEdx/djLFKf2qHdtQ/Fdfl+5Mxngq8C+ETGCorB4guKacGCeoztjylGnUOw8/jfFWQl18nbg8qKOZCDwy/K7cgJwSUR8luJyjvbfHuWOgtGZeXw56kcUp/6/THGaa62UhePOdFhHwPf6yPpr1PeA88sder/vMP5U4MzMfKRc1zdExM3l0cxuUbtu2jscqek4btXsomvXBbx+aYpT6rL8f0eKUxWmZPM9eg0sf1BQHoHYGHg8M19s8PXLZYeet5bEkubSIU4w/w+vO3MJPwgRsXKzuVQZJ4peH1/NzJej6GmmjeLc2857axY3r8ValxHxLopTfB7MzAeWMIfdM/OqJYlRxvnvzPzxErx+BYqLSF9Y0lw6xBwALNNoUVO+ZikojtKV2/3mwBPNfH4iYsVc8l4ZG+7xcjFiV7VdLdY6r2IZLyJ+Zd+RkiTVRW2uwYqIHSNiKjAtihudDe/wdDPdW0+gOI2PiPgfilOthgBfiojvNJHPWIpTZB6JiF0pjjp9l+Iagn0biTH3h0NEtEXEnlHcfHDjJualslzKOGOAv1L0KvORcjge+Gv5XKNxju3weNOIeAS4K4qb+G3VRJxtIuLBKG7mt1VEXAtMiIi/R8T7m4hzNHATcHtEHEKxt3dX4OKIaPR0zkVpqDiK4kaNq5aP96fYGzQ3l8MafbOI+I/OA3B2h8eNxvli5wE4ocPjpmXmKx2Lq8X8TM9znVR5ZKXhUxIi4uMUR8Oejog9KHoB/D7FUYCPNZHKjIi4LiI+GxErNvG6ju6O4iaV34yITRczRpXbVSXrvMJlvDBLtNNBkqQ6qs0RrPIw7djMnBIRe1Nc67F/Zt4eEXdn5pYNxrm/vLiNKO7U/YHMfD0iBlJc+PfuBuPcR9HpwvIU54ZvmZmPRcTbKXpwWmSciNge+AHF4eT3UnRFuxLFKRP7Z4P33qkilzLOg8CuWVyw3nH8uhTn4G7SYJxJmTmyfPx7it6Sro6I0cBpmbl1g3HuBD5L0WvRb4GPZ+YtETES+FFmbtNgnCkUR6yGUhwCXy8zp0dxuPyOuZ+HBuIs6MdnAF/LzJUbiNHx8zcB2CUzX4jiUPztTayrNymuzXiet0413ZuiB6PMzIMbjPMaRZE3pUOcIylOSaLDqQ+LLSKeysx1Gpx2R4qe0QZT9HA4bu7nsePnqoE4d1MUrkMotolRWdzr4p0U3bK2NRjnPorTdfel6NjiForz16/MzNebyGX/MsYnKU7lugj4VedtbRFxqtquKlnnFS7jJd6uJElqJbW5BouiW/YpAJl5aVkMXBYRR9Fcj3KvRsTm5alhM3jrZmMDae6I3ezytMQZETEz3+qI4blY8I3bOjsNGFP+2F8XOCUzt4mInSkuJm30qFEVuUB57UMX45+muZ7XOnpHZl5d5nNnFBeBNmpQZt4HEBHTM/OWMs6kJuPMLovof1Os6xfKOP9ocvl8m2IPfVfXAjX62XkzItbMoreimbx13cS/aO4i9a0pLrqdkJn/CxARO2TmQQt/2Xw2oyjylwWOz8x/RsSBzRZWEXH6gp6iPGLcoO8BH+6wI+XaiNg/M2+nyd615l7jVBZ4D5fjnozytLYGvZmZvwN+V37mPkZxD6szI+KazGzkPPQsv2++BnytLIg+BdxS5tZQYdTJkmxXlazz8r2rWMZVbFeSJLWMOhVYb0bE6vlWF+RTImIn4HfA+k3E+RxwYUTcQ7H3f2JE3Exx0e23m4jzVBSnFC4PPBQRPwAuo7iQftpCX/mWAVl2001xb5N3AmTmtRFxWg/nAkXvXxMi4le81bvY2hQ/Bs9pIs56UXTAEMBaETG0w7UzTXWL3uFx504/lqZxkyLilxQ/KP9McYHjH4EP0twpSJOAK7LotGAeUZx62IgvAH+KiN9QHEG4PiKuobix6s8aTSQzJ5SF+GERcQPFzWObPtycmU9RdG6xB0Uxc2qzMUoHUXSo8q8unmv4NFWq25FCRCyVRYczB3cYN4DmPjvtRV15xOoSigueV6Dxi5g790R4J3BnRHyJoseiRlWyXVW4zqtaxlVsV5IktYw6nSL4IWB6Zt7TafwKFHekPrGJWAMojg5tyFtHba5p5mL2KHq8+jzFj74zKE4fGktRKH0z3+r2cmExzi1ffz2wO0XPYl8sTxeblJkNXbtSRS4dYm1Cccfrjp1cXNVMBwzlqY8dTcqi97+3A3tn5pkNxtkduK5zxwYRsT7FfWG+12CcgRTd/CbFKXSjKXrAeYqi15iGet+KiI0ouhae3sVzb88Gez0rP7OfZt7P35WZ+VAjr+8i3poUPeC0ZeZ6ixOjjLMsxfV3W2VmMz/8iYjrgWMz87YunvtbZq7bYJyJwEezQw97EbEW5Y6UzFy+wTijKG7c/Ean8cOBbTPzFw3G+XJmntzItAuJ8enM/OWSxCjjVLJddYq5JOu8qmVcyXYlSVKrqE2B1RdFcSH/ocCmFNcwnJuZs8vTfVbLzCd7NUGpQRGxMvBG52J4MeJUtiNFkiSpjlri/PeIuLqn40TEChFxUkQ8FBEvRsQLUfR4d1I02NtYZr6ZmT/OzPGZ+ZMs70GTma83U1xVkUsD79Gjy6bKOIt4jx6dr7otm4iYFBHHlkcFF1tmvrikxVUZ57rOxVU5/pWqiqve+Cy3Qi7N5lNVnJ6YL0mS6qQ212BF0XNcl08BI3o6DsW1GNcDO3S40Ht1ilPzLqHxDiq6Tibi6lzwzZW7JZceWDYHNpNPVXFqNl+1WjYUvVauSHFTvWcpere7ODOfafD1i9TkZ7mSOHVa53XKpcp86jZfkiS1itqcIhgRsynuZdRVT2Lvy8yGetGqMM7DmblRs891mm5hP1B+l5lr9FQu5bS1WTYVx6nNfNVw2XTs+vsDFB1S/AfwIHBRZp7dYJyqPstVxanTOq9NLhXnU6v5kiSpVdTmCBbFD77/zMy/dn4iIhq6X1TFcZ6MiK8A58+9CDuKi83H8lYPfIsygQX/QFmxh3OBei2bKuPUab7qtmzaZeZfgL9EccPjnSnu2dRQgUV1n+Wq4tRpndcplyrzqdt8SZLUEup0DdZxLDifw3ohzieBVYCbyusGXgJuBFYGPtFgjLk/UHbsPFDco6snc4F6LZsq4xxHfearbsvmkc4jMnN2Zv4xm7unVlWf5ariHEd91nmdcqkyn6riVDVfkiS1hNqcItgXRXEj1fuyvEFnp+c+nplX9HxWUvOq+iy7TUiSpL6uNgVWFN1Ajweeobjp7VeB91Ps8f52Zr7US3GeprhB7zHA1s3GqUJVudRt2fTF+eqry6ZuarrOez2XbsqnFvMlSVKrqNMpgr8AlgXeC9wArA58F3gdOK8X47SVcdZoNk5ErBwRX4+IQ6LwtYj4XUR8PyJW6slcOsXp9WXTTXHqMF99ctlU9Vnuhm2iTuu8Drl0Rz51mS9JklpCnY5gTc7MERERwNTMXLPzc60WJyL+ANwHvA3YpHx8CUUHA+/JzD1abZ6M0zq5VBynqs9yn9sm6pRLX44jSVKrqFMvgkuVe7CXB5aLiOGZ+URErAIs3aJx3pGZH+nww2KHcvxfImJyD+dinJ6JU6dcqoxT1We5L24TdcqlL8eRJKkl1KnA+g7wUPn4YOCnEZHApsDxLRqnqh8WdZon47ROLlXGqduP7Totnzrl0pfjSJLUEmpziiBARAygyGlWRAwERgBPZ+a0VowTEfsCp5X//jfwX0D7D4ts8CavVeRinJ6LU6dcqopT1We5r24TdcqlL8eRJKkV1KrAAoiINmBtYDbwSGY+tIiX1DpOlT8s6jJPxmmtXKqKU8cf2zVbPrXJpS/HkSSp7mpTYEXE9sAPgJcpeq26FVgJeBPYPzP/3opxylhLWqTVap6M0xq5VBmnQ7xa/Niu0/KpUy59OY4kSS0jM2sxAHcDw8rH6wKXl493Bv7UinGA7YGJwHXAS8DvKH5c3Ais3YrzZJzWyaXiOFV9lvvcNlGnXPpyHAcHBwcHh1YZ6nQfrAGZOb18/BTwToDMvBZYc4Gvqnec04BdM/NDwEjgzczcBjiR4sadPZmLcXomTp1yqTLOaVTzWa4qTp2WT51y6ctxJElqCXXqRXBiRJwDXA/sTrFHm4gYCgxo0TgL/GEREaf1cC7G6Zk4dcqlyjhVfZb74jZRp1z6chxJklpCna7BGgQcStGb2D3AuZk5OyKGAKtl5pOtFicizqXoIW3uD4unM/OL5Q+LSZm5cavNk3FaJ5eK41T1We5z20SdcunLcSRJahW1KbD6In9YqK/wx7YkSVJjanMNVkQsFxEnRMT9EfFKREyPiNsjYmyrxsnMNzPzx5k5PjN/kpmzy/GvN/NDsk7zZJzWyaXKOFV9lvviNlGnXPpyHEmSWkVtjmBFxJXA5RS9i30CWBb4FXAsxWlEX221OBGxHPAV4D8ouqT+N/AYcFZmntdIHnWbJ+O0Ti4Vx6nqs9zntok65dKX40iS1DKyBl0ZlkXePZ3+n1D+XQp4qBXjAFcCY4G1gC8C/w/YADgf+HYrzpNxWieXiuNU9Vnuc9tEnXLpy3EcHBwcHBxaZajNKYLAPyJiW4CI2B14ESAz5wDRonGGZ+Z5mTk1M08Bds/MvwIHUezB78lcjNMzceqUS5Vxqvos98Vtok659OU4kiS1ht6u8OYOwHuAOyluPnoLsFE5fhhweCvGAW4Dti0f7w5c0+G5h1txnozTOrlUHKeqz3Kf2ybqlEtfjuPg4ODg4NAqQ68nME8ysAnwIWC5TuN3acU4Vf6wqMs8Gae1cqlwnmr3Y7tmy6c2ufTlOA4ODg4ODq0w9HoC7YnA4cBDwBXAE8AeHZ6b1MJxqvjxVrd5Mk4L5FJlnKzos1xVnDotnzrl0pfjODg4ODg4tMrQ6wm0JwL3zf3BBQwHJgJHlP/f3YpxKvyBUpt5Mk7r5FJxnFr92K7T8qlTLn05joODg4ODQ6sMA6mPpTJzJkBmPhEROwCXRsQ7ae5C6DrFORRoy8yZETG8fP3wzPxhL+RinJ6JU6dcqoxT1We5L24TdcqlL8eRJKkl1KkXweciYsTcf8oG+aPAqsAWLRpnnh8WwA7ArhFxCs39sKjTPBmndXKpMk5Vn+W+uE3UKZe+HEeSpJZQpxsNrwXMysxnu3hum8y8tdXiRMT1wBczc3KHcQOBc4H9MnNAT+VinJ6JU6dcKo5T1We5z20TdcqlL8eRJKlV1KbA6ov8YaG+wh/bkiRJjbHAkiRJkqSK1OkaLEmSJElqaRZYkiRJklQRCyz1KxExOyImdxiGN/n6z0XEAeXj8yJi7/LxTyNi0/LxVxuIMzwiXi9zeCAiLoiIQYsxS83kPiIiPtLh/7ERcUb5+LiIeLrM568Rcdnc+ZEkSVLjLLDU37yemSM6DE808+LMPCszL+hi/CGZ+UD57yILrNJjmTmCoqvqtYBPNJNLM8qe+kYAH1nIZKeWy2QD4GLg+ogY1l05SZIk9UV1utGw1OMiYjngSmAlYBBwbGZeWT53APBlIIF7M3P/iDgOmJmZJ3eKc2M57d7AkIiYDEwBHgNezMzTyulOBJ4v3xOAzJwdEXcCa5bTvBc4BVgOmAGMzcxp5XvcA2xPse0enJl3RsTKFN2crwf8ExiXmfeWua5fjn8K2KbMbVvgOwtbLpl5cUTsBnwa+GGDi1OSJKnfs8BSfzO3+AH4G7APsGdmvhoRqwK3R8RVwKbAscDWmTmjLGIWKTOPjojx5ZEpylMQLwNOi4ilgE8Bo4Hl574mIgYDWwFHlKcJ/gjYIzOnR8QngROBg8vJh2bmiIjYjqKo2hw4Hrg7Mz8eER8ELqA4WkU5H9tm5usRMRZoy8zx5fuOXcTsTAI2bmS+JUmSVLDAUn/z+tziB6AsaL5dFixzKI4ivR34IPDrzJwBkJkvLs6bZeYTEfFCRGxZxr07M1+IiOWB9ctib13g9+VRp80piqZrIwJgADCtQ8iLyrg3R8TbImJFYFtgr3L89RGxSkS8rZz+qsx8fXFyB2IxXydJktRvWWCpv9sPGAa8NzPfjIgngMEVv8dPgbHA6hRHneZ6rDwatSpwa0TsTnFUbUpmvn8BsTrfuG5RN7L7x2LkO9eWwMQleL0kSVK/YycX6u9WAJ4vi6sdgXeW468H9omIVQAaPUWw9GanHgEvB3YBRgHXdJ64PEp2NHAM8DAwLCLeX77voIjYrMPknyzHbwu8kpmvAH+hKBSJiB2AGZn5ahd5vUaHUxMXJiL2AsZQHjGTJElSYyyw1N9dCLRFxH3AAcBDAJk5heLap5si4h6KTicadTZwb0RcWMb6N3ADcElmzl7Aa64AhlJci7U38N3yfScDW3eY7o2IuBs4C/hsOe444L0RcS9wEnDgAt7jBmDTsiv2T3bx/BfmdtMOfAb4YGZOb2B+JUmSVIrMRZ1hJGlJlJ1bTAL2ycy/LkGcG4EvZ6an7UmSJNWUR7CkblTerPdR4M9LUlxJkiSpNXgES5IkSZIq4hEsSZIkSaqIBZYkSZIkVcQCS5IkSZIqYoElSZIkSRWxwJIkSZKkivz/eF1obaNRw8IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Impurity Analysis\n", "\n", "from matplotlib import pyplot as plt\n", "db = ep.read_db()\n", "\n", "data1 = ep.f_db(db, ReportingYear=2015, CountryName='France', NUTSRegionGeoCode='FR30',\n", " PollutantName='Carbon dioxide (CO2)')\n", "data1 = ep.change_unit(data1, unit='megaton')\n", "data2 = ep.f_db(db, ReportingYear=2015, CountryName='France', NUTSRegionGeoCode='FR30')\n", "\n", "testdata=ep.get_PollutantVolume(data1, FirstOrder='FacilityReportID').sort_values('TotalQuantity', ascending=False)\n", "testdata2=ep.get_ImpurityVolume(data2, target='Carbon dioxide (CO2)', impurity='Nitrogen oxides (NOx/NO2)', absolute=True).sort_values('TotalQuantity', ascending=False)\n", "testdata2.loc[:,'Nitrogen oxides (NOx/NO2)'] /= 1000000\n", "\n", "testdata3=ep.get_ImpurityVolume(data2, target='Carbon dioxide (CO2)', impurity='Nitrogen oxides (NOx/NO2)').sort_values('TotalQuantity', ascending=False)\n", "testdata3.loc[:,'Nitrogen oxides (NOx/NO2)'] *= 1000\n", "testdata4=ep.get_ImpurityVolume(data2, target='Carbon dioxide (CO2)', impurity='Nitrogen oxides (NOx/NO2)', statistics=True)\n", "testdata4.loc[:,'Nitrogen oxides (NOx/NO2)'] *= 1000\n", "testdata4 = testdata4.drop('count')\n", "\n", "fig3, ax = plt.subplots(2, 2, figsize=(12, 8))\n", "\n", "testdata.plot(x='FacilityReportID', y='TotalQuantity', ax=ax[0, 0], kind='bar')\n", "testdata2.plot(x='FacilityReportID', y='Nitrogen oxides (NOx/NO2)', ax=ax[0,1], kind='bar')\n", "testdata3.plot(x='FacilityReportID', y='Nitrogen oxides (NOx/NO2)', ax=ax[1,0], kind='bar')\n", "testdata4.plot(y='Nitrogen oxides (NOx/NO2)', ax=ax[1,1], kind='bar', rot=0)\n", "\n", "plt.tight_layout()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 4 }