{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# _Dataviz Code Session_: Dados Geoespaciais" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Objetivos da DCS\n", "\n", "- Aplicar técnicas de dataviz para plotagem e manipulação de representações visuais geoespaciais relacionadas a campos e poços de petróleo em bacias sedimentares brasileiras." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Ferramentas utilizadas\n", "\n", "- Módulos Python\n", " - `numpy`\n", " - `pandas`\n", " - `matplotlib`\n", " - `seaborn`\n", " - `geopandas`" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Aplicação do modelo referencial\n", "\n", "- Vide {numref}`Capítulo %s `." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "import numpy as np, pandas as pd, matplotlib.pyplot as plt, seaborn as sb, geopandas as geo" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "tags": [ "remove-cell" ] }, "outputs": [], "source": [ "from warnings import filterwarnings; filterwarnings('ignore');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Dados de entrada" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Shapefiles dos campos de produção (ANP) - [[Fonte]](https://www.gov.br/anp/pt-br/assuntos/exploracao-e-producao-de-oleo-e-gas/dados-tecnicos/shapefile-de-dados) " ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "geocp = geo.read_file('../data/ANP_campos_shp/CAMPOS_PRODUCAO_SIRGASPolygon.shp')" ] }, { "cell_type": "code", "execution_count": 12, "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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
NUM_RODADANOM_CAMPOAREAOPERADOR_CNUM_CONTRADAT_ASSINADAT_TERMINNOM_BACIACOD_CAMPOSIG_CAMPODAT_DESCOBDAT_INICIOETAPAMED_LAMINAFLUIDO_PRIIDgeometry
0Rodada 11IRARA6.0BGM Petróleo e Gás Ltda.48610005457201381NoneNoneEspírito Santo4841.0IRAR03/01/202213/07/2022Produção0.0ÓLEO2460.0POLYGON ((-39.81288 -19.25048, -39.81288 -19.2...
1Rodada 0PESCADA55.0Petróleo Brasileiro S.A.48000003912978406-08-1998NonePotiguar213.0PE24-10-198001-05-1999Produção19.0GÁS1100.0POLYGON ((-36.76596 -4.73997, -36.76856 -4.739...
2Rodada 0SALGO11.0Petróleo Brasileiro S.A.48000003841973806-08-1998NoneSergipe275.0SG01-12-198431-12-1986Em Devolução14.0ÓLEO1102.0POLYGON ((-36.57586 -10.63324, -36.57846 -10.6...
3Rodada 0CAXAREU34.0Petróleo Brasileiro S.A.48000003560974906/08/1998NoneCampos1164.0CXR07-08-2002NoneDesenvolvimento1141.0ÓLEO1125.0POLYGON ((-40.07812 -21.39323, -40.07812 -21.3...
4Rodada 0CEXIS54.03R Candeias S.A.48000003641974906-08-1998NoneRecôncavo88.0CX07-03-196630-06-1966Produção0.0GÁS1126.0POLYGON ((-38.4639 -12.62544, -38.4613 -12.625...
......................................................
422AM2Porto do Mangue6.0Máxima 07 Exploração e Produção de Petróleo Ltda48610007488200648NoneNonePotiguar4535.0420110/10/2011NonePlano de Avaliação (Teste de longa duração)0.0None2623.0POLYGON ((-36.79461 -5.09674, -36.79461 -5.099...
423Rodada 0MIRANGA LESTE1.0Petróleo Brasileiro S.A.480000036759761NoneNoneRecôncavo182.0MGLNoneNoneEm Devolução0.0OUTROS2661.0POLYGON ((-38.14358 -12.32336, -38.14619 -12.3...
424Rodada 13 (AM)SÃO JOÃO6.0Petro-Victory Energia Ltda.4861000200520169112/5/2016NoneBarreirinhas4695.0SJ01/08/196601/07/2020Produção0.0Óleo2662.0POLYGON ((-43.32852 -2.59932, -43.32852 -2.596...
425Rodada 9GAVIÃO CARIJÓ163.0Eneva S.A.48610001413200815NoneNoneParnaíba4629.0GVCANoneNoneDesenvolvimento0.0GÁS2700.0POLYGON ((-44.33073 -5.66927, -44.33073 -5.684...
426Rodada 0LESTE DE POÇO XAVIER22.0Potiguar E&P S.A.4861000400098NoneNonePotiguar167.0LPXNoneNoneProdução0.0ÓLEO2701.0POLYGON ((-37.67708 -5.52122, -37.67448 -5.521...
\n", "

427 rows × 17 columns

\n", "
" ], "text/plain": [ " NUM_RODADA NOM_CAMPO AREA \\\n", "0 Rodada 11 IRARA 6.0 \n", "1 Rodada 0 PESCADA 55.0 \n", "2 Rodada 0 SALGO 11.0 \n", "3 Rodada 0 CAXAREU 34.0 \n", "4 Rodada 0 CEXIS 54.0 \n", ".. ... ... ... \n", "422 AM2 Porto do Mangue 6.0 \n", "423 Rodada 0 MIRANGA LESTE 1.0 \n", "424 Rodada 13 (AM) SÃO JOÃO 6.0 \n", "425 Rodada 9 GAVIÃO CARIJÓ 163.0 \n", "426 Rodada 0 LESTE DE POÇO XAVIER 22.0 \n", "\n", " OPERADOR_C NUM_CONTRA \\\n", "0 BGM Petróleo e Gás Ltda. 48610005457201381 \n", "1 Petróleo Brasileiro S.A. 480000039129784 \n", "2 Petróleo Brasileiro S.A. 480000038419738 \n", "3 Petróleo Brasileiro S.A. 480000035609749 \n", "4 3R Candeias S.A. 480000036419749 \n", ".. ... ... \n", "422 Máxima 07 Exploração e Produção de Petróleo Ltda 48610007488200648 \n", "423 Petróleo Brasileiro S.A. 480000036759761 \n", "424 Petro-Victory Energia Ltda. 48610002005201691 \n", "425 Eneva S.A. 48610001413200815 \n", "426 Potiguar E&P S.A. 4861000400098 \n", "\n", " DAT_ASSINA DAT_TERMIN NOM_BACIA COD_CAMPO SIG_CAMPO DAT_DESCOB \\\n", "0 None None Espírito Santo 4841.0 IRAR 03/01/2022 \n", "1 06-08-1998 None Potiguar 213.0 PE 24-10-1980 \n", "2 06-08-1998 None Sergipe 275.0 SG 01-12-1984 \n", "3 06/08/1998 None Campos 1164.0 CXR 07-08-2002 \n", "4 06-08-1998 None Recôncavo 88.0 CX 07-03-1966 \n", ".. ... ... ... ... ... ... \n", "422 None None Potiguar 4535.0 4201 10/10/2011 \n", "423 None None Recôncavo 182.0 MGL None \n", "424 12/5/2016 None Barreirinhas 4695.0 SJ 01/08/1966 \n", "425 None None Parnaíba 4629.0 GVCA None \n", "426 None None Potiguar 167.0 LPX None \n", "\n", " DAT_INICIO ETAPA MED_LAMINA \\\n", "0 13/07/2022 Produção 0.0 \n", "1 01-05-1999 Produção 19.0 \n", "2 31-12-1986 Em Devolução 14.0 \n", "3 None Desenvolvimento 1141.0 \n", "4 30-06-1966 Produção 0.0 \n", ".. ... ... ... \n", "422 None Plano de Avaliação (Teste de longa duração) 0.0 \n", "423 None Em Devolução 0.0 \n", "424 01/07/2020 Produção 0.0 \n", "425 None Desenvolvimento 0.0 \n", "426 None Produção 0.0 \n", "\n", " FLUIDO_PRI ID geometry \n", "0 ÓLEO 2460.0 POLYGON ((-39.81288 -19.25048, -39.81288 -19.2... \n", "1 GÁS 1100.0 POLYGON ((-36.76596 -4.73997, -36.76856 -4.739... \n", "2 ÓLEO 1102.0 POLYGON ((-36.57586 -10.63324, -36.57846 -10.6... \n", "3 ÓLEO 1125.0 POLYGON ((-40.07812 -21.39323, -40.07812 -21.3... \n", "4 GÁS 1126.0 POLYGON ((-38.4639 -12.62544, -38.4613 -12.625... \n", ".. ... ... ... \n", "422 None 2623.0 POLYGON ((-36.79461 -5.09674, -36.79461 -5.099... \n", "423 OUTROS 2661.0 POLYGON ((-38.14358 -12.32336, -38.14619 -12.3... \n", "424 Óleo 2662.0 POLYGON ((-43.32852 -2.59932, -43.32852 -2.596... \n", "425 GÁS 2700.0 POLYGON ((-44.33073 -5.66927, -44.33073 -5.684... \n", "426 ÓLEO 2701.0 POLYGON ((-37.67708 -5.52122, -37.67448 -5.521... \n", "\n", "[427 rows x 17 columns]" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "geocp" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Shapefiles do território nacional" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "geobr = geo.read_file('../data/gadm36_BRA_shp/gadm36_BRA_1.shp')" ] }, { "cell_type": "code", "execution_count": 14, "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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
GID_0NAME_0GID_1NAME_1VARNAME_1NL_NAME_1TYPE_1ENGTYPE_1CC_1HASC_1geometry
0BRABrazilBRA.1_1AcreNoneNoneEstadoStateNoneBR.ACPOLYGON ((-68.10553 -10.72192, -68.10547 -10.7...
1BRABrazilBRA.2_1AlagoasNoneNoneEstadoStateNoneBR.ALMULTIPOLYGON (((-35.88986 -9.84431, -35.88986 ...
2BRABrazilBRA.3_1AmapáNoneNoneEstadoStateNoneBR.APMULTIPOLYGON (((-50.8257 2.52208, -50.8257 2.5...
3BRABrazilBRA.4_1AmazonasAmazoneNoneEstadoStateNoneBR.AMPOLYGON ((-58.13699 -7.35614, -58.13682 -7.356...
4BRABrazilBRA.5_1BahiaBa¡aNoneEstadoStateNoneBR.BAMULTIPOLYGON (((-38.69208 -17.95958, -38.69208...
5BRABrazilBRA.6_1CearáNoneNoneEstadoStateNoneBR.CEMULTIPOLYGON (((-40.83181 -2.88125, -40.83181 ...
6BRABrazilBRA.7_1Distrito FederalNoneNoneDistrito FederalFederal DistrictNoneBR.DFPOLYGON ((-47.93915 -16.05135, -48.0211 -16.05...
7BRABrazilBRA.8_1Espírito SantoEspiritu SantoNoneEstadoStateNoneBR.ESMULTIPOLYGON (((-40.88403 -21.16125, -40.88403...
8BRABrazilBRA.9_1GoiásGoiáz|GoyázNoneEstadoStateNoneBR.GOPOLYGON ((-49.0048 -18.53718, -49.01248 -18.53...
9BRABrazilBRA.10_1MaranhãoSão Luíz de MaranhãoNoneEstadoStateNoneBR.MAMULTIPOLYGON (((-44.54986 -1.88458, -44.54986 ...
10BRABrazilBRA.12_1Mato GrossoMatto GrossoNoneEstadoStateNoneBR.MTPOLYGON ((-52.47143 -16.12672, -52.48842 -16.1...
11BRABrazilBRA.11_1Mato Grosso do SulNoneNoneEstadoStateNoneBR.MSPOLYGON ((-54.16836 -23.9992, -54.17063 -24.00...
12BRABrazilBRA.13_1Minas GeraisMinas|Minas GeraesNoneEstadoStateNoneBR.MGPOLYGON ((-50.84402 -19.96758, -50.85332 -19.9...
13BRABrazilBRA.14_1ParáNoneNoneEstadoStateNoneBR.PAMULTIPOLYGON (((-47.29792 -0.63375, -47.29792 ...
14BRABrazilBRA.15_1ParaíbaParahybaNoneEstadoStateNoneBR.PBMULTIPOLYGON (((-43.01208 -9.40801, -43.01453 ...
15BRABrazilBRA.16_1ParanáNoneNoneEstadoStateNoneBR.PRMULTIPOLYGON (((-48.36514 -25.73403, -48.36514...
16BRABrazilBRA.17_1PernambucoPernamboucNoneEstadoStateNoneBR.PEMULTIPOLYGON (((-35.13597 -8.83792, -35.13597 ...
17BRABrazilBRA.18_1PiauíPiauhyNoneEstadoStateNoneBR.PIPOLYGON ((-42.69078 -9.54547, -42.69122 -9.545...
18BRABrazilBRA.19_1Rio de JaneiroNoneNoneEstadoStateNoneBR.RJMULTIPOLYGON (((-44.67125 -23.35458, -44.67125...
19BRABrazilBRA.20_1Rio Grande do NorteNoneNoneEstadoStateNoneBR.RNMULTIPOLYGON (((-35.10958 -6.19375, -35.10958 ...
20BRABrazilBRA.21_1Rio Grande do SulNoneNoneEstadoStateNoneBR.RSMULTIPOLYGON (((-52.07069 -32.02847, -52.07069...
21BRABrazilBRA.22_1RondôniaGuaporéNoneEstadoStateNoneBR.ROPOLYGON ((-62.89219 -12.86014, -62.89212 -12.8...
22BRABrazilBRA.23_1RoraimaRio BrancoNoneEstadoStateNoneBR.RRPOLYGON ((-60.0379 0.26349, -60.04676 0.23211,...
23BRABrazilBRA.24_1Santa CatarinaSanta CatharinaNoneEstadoStateNoneBR.SCMULTIPOLYGON (((-48.54236 -27.9318, -48.54236 ...
24BRABrazilBRA.25_1São PauloNoneNoneEstadoStateNoneBR.SPMULTIPOLYGON (((-48.08236 -25.28431, -48.08236...
25BRABrazilBRA.26_1SergipeNoneNoneEstadoStateNoneBR.SEMULTIPOLYGON (((-37.38458 -11.45986, -37.38458...
26BRABrazilBRA.27_1TocantinsNoneNoneEstadoStateNoneBR.TOPOLYGON ((-49.23736 -12.88397, -49.23532 -12.8...
\n", "
" ], "text/plain": [ " GID_0 NAME_0 GID_1 NAME_1 VARNAME_1 \\\n", "0 BRA Brazil BRA.1_1 Acre None \n", "1 BRA Brazil BRA.2_1 Alagoas None \n", "2 BRA Brazil BRA.3_1 Amapá None \n", "3 BRA Brazil BRA.4_1 Amazonas Amazone \n", "4 BRA Brazil BRA.5_1 Bahia Ba¡a \n", "5 BRA Brazil BRA.6_1 Ceará None \n", "6 BRA Brazil BRA.7_1 Distrito Federal None \n", "7 BRA Brazil BRA.8_1 Espírito Santo Espiritu Santo \n", "8 BRA Brazil BRA.9_1 Goiás Goiáz|Goyáz \n", "9 BRA Brazil BRA.10_1 Maranhão São Luíz de Maranhão \n", "10 BRA Brazil BRA.12_1 Mato Grosso Matto Grosso \n", "11 BRA Brazil BRA.11_1 Mato Grosso do Sul None \n", "12 BRA Brazil BRA.13_1 Minas Gerais Minas|Minas Geraes \n", "13 BRA Brazil BRA.14_1 Pará None \n", "14 BRA Brazil BRA.15_1 Paraíba Parahyba \n", "15 BRA Brazil BRA.16_1 Paraná None \n", "16 BRA Brazil BRA.17_1 Pernambuco Pernambouc \n", "17 BRA Brazil BRA.18_1 Piauí Piauhy \n", "18 BRA Brazil BRA.19_1 Rio de Janeiro None \n", "19 BRA Brazil BRA.20_1 Rio Grande do Norte None \n", "20 BRA Brazil BRA.21_1 Rio Grande do Sul None \n", "21 BRA Brazil BRA.22_1 Rondônia Guaporé \n", "22 BRA Brazil BRA.23_1 Roraima Rio Branco \n", "23 BRA Brazil BRA.24_1 Santa Catarina Santa Catharina \n", "24 BRA Brazil BRA.25_1 São Paulo None \n", "25 BRA Brazil BRA.26_1 Sergipe None \n", "26 BRA Brazil BRA.27_1 Tocantins None \n", "\n", " NL_NAME_1 TYPE_1 ENGTYPE_1 CC_1 HASC_1 \\\n", "0 None Estado State None BR.AC \n", "1 None Estado State None BR.AL \n", "2 None Estado State None BR.AP \n", "3 None Estado State None BR.AM \n", "4 None Estado State None BR.BA \n", "5 None Estado State None BR.CE \n", "6 None Distrito Federal Federal District None BR.DF \n", "7 None Estado State None BR.ES \n", "8 None Estado State None BR.GO \n", "9 None Estado State None BR.MA \n", "10 None Estado State None BR.MT \n", "11 None Estado State None BR.MS \n", "12 None Estado State None BR.MG \n", "13 None Estado State None BR.PA \n", "14 None Estado State None BR.PB \n", "15 None Estado State None BR.PR \n", "16 None Estado State None BR.PE \n", "17 None Estado State None BR.PI \n", "18 None Estado State None BR.RJ \n", "19 None Estado State None BR.RN \n", "20 None Estado State None BR.RS \n", "21 None Estado State None BR.RO \n", "22 None Estado State None BR.RR \n", "23 None Estado State None BR.SC \n", "24 None Estado State None BR.SP \n", "25 None Estado State None BR.SE \n", "26 None Estado State None BR.TO \n", "\n", " geometry \n", "0 POLYGON ((-68.10553 -10.72192, -68.10547 -10.7... \n", "1 MULTIPOLYGON (((-35.88986 -9.84431, -35.88986 ... \n", "2 MULTIPOLYGON (((-50.8257 2.52208, -50.8257 2.5... \n", "3 POLYGON ((-58.13699 -7.35614, -58.13682 -7.356... \n", "4 MULTIPOLYGON (((-38.69208 -17.95958, -38.69208... \n", "5 MULTIPOLYGON (((-40.83181 -2.88125, -40.83181 ... \n", "6 POLYGON ((-47.93915 -16.05135, -48.0211 -16.05... \n", "7 MULTIPOLYGON (((-40.88403 -21.16125, -40.88403... \n", "8 POLYGON ((-49.0048 -18.53718, -49.01248 -18.53... \n", "9 MULTIPOLYGON (((-44.54986 -1.88458, -44.54986 ... \n", "10 POLYGON ((-52.47143 -16.12672, -52.48842 -16.1... \n", "11 POLYGON ((-54.16836 -23.9992, -54.17063 -24.00... \n", "12 POLYGON ((-50.84402 -19.96758, -50.85332 -19.9... \n", "13 MULTIPOLYGON (((-47.29792 -0.63375, -47.29792 ... \n", "14 MULTIPOLYGON (((-43.01208 -9.40801, -43.01453 ... \n", "15 MULTIPOLYGON (((-48.36514 -25.73403, -48.36514... \n", "16 MULTIPOLYGON (((-35.13597 -8.83792, -35.13597 ... \n", "17 POLYGON ((-42.69078 -9.54547, -42.69122 -9.545... \n", "18 MULTIPOLYGON (((-44.67125 -23.35458, -44.67125... \n", "19 MULTIPOLYGON (((-35.10958 -6.19375, -35.10958 ... \n", "20 MULTIPOLYGON (((-52.07069 -32.02847, -52.07069... \n", "21 POLYGON ((-62.89219 -12.86014, -62.89212 -12.8... \n", "22 POLYGON ((-60.0379 0.26349, -60.04676 0.23211,... \n", "23 MULTIPOLYGON (((-48.54236 -27.9318, -48.54236 ... \n", "24 MULTIPOLYGON (((-48.08236 -25.28431, -48.08236... \n", "25 MULTIPOLYGON (((-37.38458 -11.45986, -37.38458... \n", "26 POLYGON ((-49.23736 -12.88397, -49.23532 -12.8... " ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "geobr" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Planilhas de poços públicos (2022 e 2023)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "p2022 = pd.read_excel('../data/ANP_pocos_publicos/planilha-pocos-publicos-2022.xlsx')\n", "p2023 = pd.read_excel('../data/ANP_pocos_publicos/planilha-pocos-publicos-2023.xlsx')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Processamento" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "def filter_by_state(state,fluid,phase):\n", " \"\"\"Seleciona poços e campos para estado, bacia e propriedades especificados.\"\"\"\n", " \n", " a = p2022[p2022.ESTADO == state[0]]\n", " b = p2023[p2023.ESTADO == state[0]]\n", " c = geobr[geobr.NAME_1 == state[1]]\n", " d = geocp[(geocp.NOM_BACIA == basin) & (geocp.FLUIDO_PRI == fluid) & (geocp.ETAPA == phase)]\n", " \n", " return a, b, c, d" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Escolha de parâmetros e aplicação de filtragem" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "state = ('BA','Bahia')\n", "basin = 'Recôncavo'\n", "fluid = 'ÓLEO'\n", "stage = 'Produção'\n", "p2022_s, p2023_s, geo_s, geo_cp = filter_by_state(state,fluid,stage)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "# casting de coordenadas\n", "p2023_s['LATITUDE_BASE_DD'] = p2023_s['LATITUDE_BASE_DD'].str.replace(',','.').astype(float)\n", "p2023_s['LONGITUDE_BASE_DD'] = p2023_s['LONGITUDE_BASE_DD'].str.replace(',','.').astype(float)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Visualizações" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Mapa geográfico de território e campos" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAGdCAYAAADqhP93AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABjTUlEQVR4nO3dd3Tc9Z3v/+f0GZUZ9d4sybJky0XGuILBGDDBtGwCmEAS4ixskuXuZrO7914295CQc9ncc7K7v9xt9xxvnACBhFBDMcEFF+FeZavYktVlFVt9VEZTv78/xhY2WLZGntF3Rn4/ztHBmvZ9zxd7XvP5fppGURQFIYQQIgBatQsQQggReSQ8hBBCBEzCQwghRMAkPIQQQgRMwkMIIUTAJDyEEEIETMJDCCFEwCQ8hBBCBEyvdgE3yufz0dHRQWxsLBqNRu1yhBAioimKwtDQEBkZGWi1E7cvIj48Ojo6yM7OVrsMIYSYUdra2sjKyprw/ogPj9jYWMD/Rq1Wq8rVCCFEZLPb7WRnZ49/tk4k4sPj0qUqq9Uq4SGEEEFyvW4A6TAXQggRMAkPIYQQAZPwEEIIETAJDyGEEAGT8BBCCBEwCQ8hhBABk/AQQggRMAkPIYQQAZPwEEIIETAJDyGEEAGT8BBCCBEwCQ8hhBABk/AQQkQ0t9tNT08PXq9X7VJuKhG/qq4Q4ubk8/k4e/Ysez/6CG9/P/ElJZQtW0ZMTAypqamyOVyISXgIISKO1+tl+9attO7dS9bICKXp6ew9cICdx45BVBSLHnqIpUuXql3mjCbhIYSIGENDQ5w4doz6o0fxdHRwa1QUC0tK0Gq1PBEfj9vr5WBjI3UnTrBo0SKMRqPaJc9YEh5CiLDn8Xioq6tj70cfoWttZaHVSkF2NvHR0eOP0Wg0GPV65mVkUF9ZSfmuXdy9bp2KVc9sEh5CiLDU0tLCvm3bGBsawjM2Bn19FCoKt82di1E/8UdXYkwM2Xo9QwMD01fsTUjCQ4gI5HA4OHr0KLNnzyYtLW389q6uLjweD16vF71ej8PhYNasWeh0OhWrnZrOzk6GjxxheVoaZoOBpMxMEi5raVyLUa9ndGSEzs5OaqqrKS4pITMzM8QV31wkPISIQJ999hnNR49yescOcleuJCUlhdOVlQw3NoLDATodeL0QF8dXnnqK7OxstUsOWFpaGqfi4shNTMRmsQT03NkpKdRWVPBhfT309dGQkkJUejpphYWkpKaSk5NDXFxcaAq/SUh4CBGB9DodusZGlmVmUrdlC0c1GsweD4usVmanpTHkdFLT0UFrWhopKSlqlzslWVlZEBdHx8BAwOGRERfHE2YzY243tqwsztvtdLa20nXmDEeAI9nZ3Pbgg8yZMyc0xd8EJDyEiDB2u536I0dYmp5OaWYmpZmZuDwe9FotWq1/3m98dDQGnY7WpiZeeeUVsrKyyMvLY+7cuSpXf31ut5vu7m7/SCmtFq/PN6XXiTWbiTWbAchOSCA7IQEAj9fL3oYG9rz+OqNf/So5OTn09PTQ3d2Nz+cjOjqa/Px84uPjg/aeZiIJDyEiyPDwMG/85jfQ2krxvHnjt1+tAzndZuPOsTEqT56k58gRzuXkYLPZiI6ODqtLNv39/XR2dtLf349er+dsRQWjjY1gNKK320kIcl+FXqfjzqIifLW1HHn7bY7o9TA4iM3rxaDR0AyczMvj3sce87d+vmBsbIze3l6Sk5Nv6qHAGkVRFLWLuBF2ux2bzcbg4CBWq1XtcoQImc7OTnZ8/DHew4dZk5dHbmLipJ/bMTDAR0NDcPHDLjo3l/yCAlasWhWqcieltbWVT377W7Td3Vh9PsaANJ2OxdnZuL1ebBYLUSH8gG7r68Pr85GdkIDuYqvN4/WyvbaW9rQ0lt5/P3PnzkV/MZztdjvvvvIKrvZ2NDYbRStWsHLVKgwGQ8hqnG6T/UyVlocQEWB0dJSt772H7+hRHpk/f9Kjji5Jiokhv7OTVKuVvpER6j79lMq6OlLT08nPzw9R1dfX3t6OqaODJ0tL0aswIuzSpazL6XU61hUXs7+xkUMvv0zHXXdx3wMPAFC+cyem2lrWz5pFl93O4XffxWgysWLFiukuXXWyMKIQEeDkyZO46ut5uLQ04OAA/2Wtu0tKmJ+ZyR1FRTyzahWJAwOcqakJQbWTl5+fjzMpiebeXlXr+CKtVstthYXcHhdHa0UFTqeT8+fP03H0KKuyskiOjaU4LY0orxeXy6V2uaqQ8BAizI2OjlJ3/Dj5ijKl4JhIaUoK5w4c4MyZM0F7zUClpqaSUlpKS5iFxyWJ0dHgcDA0NITD4QCHg5TYWFweD+/X1DBWWEhJSYnaZapCLlsJEaYaGhqoPnWKrq4uOH2aZUVFQV0pdk5aGqdPnuTovn0UFRWNj9SaTj6fj1G7nfgwncQYazajGR2lo6MDs9kMPh86rZb2/n76EhN54PHHI3Yo9I2S8BAizJw7d46amhqajx2D06fBbKbIZhsfdhpMc1JS+Ky9nbGxMaKiooL++tez97PPGK2spDQnZ9qPPRlmg4G5BgOHP/gAfVISqV4vBp2ONJsNfV0dDQ0NZGRkjD++s7OTpqYmcnJyrjpSayaR8BAiTLS2tnJs/366Ozrg5ElK09JYumIFOq02ZHtTZMXHo21ooLq6mltvvTUkx5jI2bNnObN9O6sTEkiMiZnWYwdieX4+ts5Ous+epezi4AKzwcCKhAQ++/BDxkZHWbFqFW1tbex95x2UlhbqFizgkW9+M6yGRAebhIcQYaLq5EmGt21jQUoKcxcvxhrgrOqpiDWbydfpqK+qClp4dHZ24vP5yMjIuCL03G43o6Oj45MA93/0EUUOB8Vh2uq4RKfVUnqVuSYl6emYe3r47KOP+N3hwzA6ylyvlyWLFvFmZSXHjx3jrrVrVah4ekh4CBEmPC4XWfHxLJ/mobMZcXHUd3QwMjJC9A12yB87doxj778PHg+Zy5ez9t57MZvNdHR0sPXtt3H39IDHg2ZsjFzg9qKi4LwJlcxKSiIrPp6W3l4scXFkxsdzsLGRsYQE4mb4DHUJDyHCgN1uZ6CrixQVJpul22zoGhp4ffNmUnNzuffee7FModUzPDzMsU8+YbHHQ5rNxs7t23m7t5f07Gwa9u4lta+PJZmZmAwGYkwmzDNkYp1Bp6Pwsk7zup4eKC2luLhYxapCT8JDCJU1NTWx6733iG5qYr4K38RtFgv3Z2ezt7qa83V1/LahgUW33UZubi6pqamTfp2Ghgb0Fy6wsLgYg07H/QYD+yoq6D1+nNsSEymZN++m2Fc83mSiNzraPzprBpPwEEJFHo+HT999l8yWFtaWlFxzk6NQSrfZePSWW+gdHqaivZ2Kl1/m9Lx5PLZx46RbIX19fSQoCoaLw26TYmJ4+LL1t24WaTYbXT09dHZ2zug9RGSSoBBB0t7eTkdHR0DPOX78OLS3s2zWLNWC43KJMTGsnTOHJQkJOLu78Xq9k35udHQ09hDWFiluyckhvbOT/Tt3ql1KSEl4CBEkWz76iI8++mjSjx8ZGeHM4cOUGo1BnTkeDB0DA1jz84kJYAjtkN2OLE3qX9okPTaW/s5OtUsJKfW/6ghxDYqi0NPTg9frJSUlRZVZ0JPh8/lAUYi6LAQURaGuro7BwUHAP1Q1NjaWgoICDh06RP3+/dDfT2EYXtro8HgwBriPxvnGRvJVmGgYjmwWC4yO0tPTQ1JSktrlhISEhwhbbrebTz78kM7jx8HrJXHBAh782tfCcg+Fo0ePwsXdDSoqKlAUhbbmZroOHMA8NsaYx+N/YFQUBzMz0TQ2siA6mpSEBJLCaIKcoigMOhzgdBJvs036eb29vQx3dJAm2yIAkJ+UxNHqao4fPcq9992HoigzbrCAhIcIO263m6qqKs5WVzN6+DD3ZWRg0OnYeuAAe5OTuevuu9UuEQCv10tNTQ0dHR20HDoE588zajZz8sgRAOIUhfszM8mKj8fhcqHRaBh1uahqb6cwO5uMMJt9rCgKH1dV0e7zwYIF3Lps2aSfe/zoUaw9PWSXloawwsih1WpZlJxM+ZEjbGppoWzx4mmfwR9qEh4irHi9Xj756CO6y8vJAtbm5Y0vXbFybIzd5eXMKSlRfRTL2NgYu3ftonXHDtLcbm61WplfWorL673q5kWWi7eZDQZWh+nEOK/PR5fHAwsX8sif/dl1F/xzOBzjmyA1nzjBisTE8Q2VhL/1UV5TA8nJqv99DQUJDxE2fD4f5Xv2cKG8nAcyM0n9wiWQotRUzlRV8en771O8ZAmLFy8e3+EN/MNem5qa6O/vp7i4OKQ7Sx4+fJiuP/2Ju1NTyU9OHr9djQ2NgkWv05Gu1zMSE3Pd4Ghra2Pr739PVEICRYsWoYyMYA2zTn+1GfV68vR6ms+d8y/nPsNIeIiwMDw8zKeffMKFfftYk5j4peC45N45c9hbX8+pujrONTWxeu1aBgcHGRwcpK+vj8YdO8DnoyIlhYw5c1hUVobdbicqKoq8vLwbrnN0dJT9e/fSePAgpRbLFcExE4x5vURPsKyGoihcuHABk8nE0QMHiGtuRn/uHMerq8nXaMiagd+ub9S9c+fyTkUFx/fvp6CgQO1ygkrCQ6jG4XBQX19Pe1sbHTU1mJubeSgvb8LgAP9ln7uLi+kaHOSTPXt4t7kZ+vrA6QSrlSK3m4y4OJrPnaO5pYWOQ4dgdBRiY7nrm9+ksLDwhmr++L33cBw4wO3JycyZNeuGXivcNHZ30xMdTekE4fHptm00lpeDwQCjo9yfm0tWfDwujycs5qiEK71Ox3m7naqqKqxWK1qtlpSUlLAc+BEI+T8uVNHX18cHr7+Ot6GBVEWhLDaWeXPnTvpDKM1m457cXJp7eylKTcXt9TLsdFI4ezZarZai1FSae3qwGI1Y09J4/8wZOjo6phweZ8+epaenh766Otamp1Mww1ocALXnz4PVStWePThGRlhYVkZiYiKjFzdDajx0iKVaLTF6PRqbjayLISPBcW23Zmdz/MwZ9re1gVYLHg/a7GzW3H8/ubm5V1x6jSSRWbWIaKOjo+z85BNi6up4oLh4ygvkZcbHk3mNlUvzLhtf79Nqxzt3J0NRFAYGBsb3p9719tvQ0kJhdDSzZs+eUr3hLjchAU1fH5k6Hac++ICG8nIsGRk47XZ8PT1kezzMnzNHOsUDlBEXR0ZcHN1DQyiKgk9R2H3mDJ+2txNfWsrDjz0Wka0QCQ8x7fbs3Ilj/37uz8+ftpVVM41GGmtqyMvLIzU19UuTDX0XJ8Rdur18zx5qd+0CtxuMRjQ9PXy1qCis5mQE29yMDOZe3BVvbno65+122lpbMen15GdlTcv+IjNZcmzs+J8fX7CA9oEBPj52jKZbb2XOnDkqVjY1Eh5iWp07d462gwe5JyNjWnePm52SQsPJk3zY1ET+6tXk5eeTmJhI/MWWyydbtnC+uZmChQs539FBf2UlC30+8hITcXo8pJSUzJglxCdDp9WOf2MWwafRaOgcHISkJNLS0tQuZ0okPMS0amhoIKavj1nTvHtcRlwcT0ZH09jTw2fbttGo1WItKWHDd75Df38/5+rr4cQJOs6eJVZRKLDZWDzDOsRF+Bh1uTjR3o5u1SpsAczkDycSHmJapaSkUGsyqTJCx2QwUJKeTqzZTFNPD6c7OhgYGOCPr71GTH09jyxdetUJfkIEW8/wMOTlsTiCZ51Lz5eYVkajEcxm/z8elWTFxzM/MxO6u3nzP/8T99mz3FtQIMEhpk2MyYRmYIDmhgaGhobULmdKJDzEtEpISIDYWP+wUBXFRUWxLjWVOzwevj5r1ozuCBfhJyE6mmKjke59+/jwj39Uu5wpkfAQ0yo+Pp4VDz7IWY2GbpW/ceUmJjInLS3s9tIQM5/d4eB0fz/mhQtZvWaN2uVMiYSHmHbz5s3DWlpKzQzfLEeIiXQMDkJBAQ889BBZWVlqlzMlEh5i2mm1WqJiY3G43WqXIoQqzHo9uN2YzWa1S5kyCQ+hirzZs2k1m9lbX8/Q2Jja5QgxrRKio2FoiLq6OrVLmTIJD6GK0tJSbt2wgbMZGfyhro49dXU4Li4FIsRMZ7VYyHO7ObJlC0ePHsXpdKpdUsA0inJx78wIZbfbsdlsDA4OhnT/BhEaTqfTv5/3iRNo6+r4+pw5xERwU16IyfL5fOxvbKRGo2HpE0+waNEitUsCJv+ZKi0PoSqTycTq1atZdd99uGw2+kZH1S5JiGmh1WqZnZICGs34MjmRJKTh8dJLL7Fy5UqioqKIm2CNnNbWVh588EGio6NJSkrir/7qr8ZXMhU3h/Pnz3O4vJy4oSFyEhLULkeIadM5OIg+QrepDWl4uFwuHn30Ub7//e9f9X6v18v69esZGRlh7969vPHGG7zzzjv87d/+bSjLEmHm4Gef4ThwgLvz89UuRYhpdXZgAHNqakTu6RHSil988UUAXn755avev23bNmpqamhrayPj4lLQ//zP/8zTTz/NSy+9JH0YN4mE1FTOx8ejl30ixE2m0GbjaHMzHR0dJCYmYjKZ1C5p0lSNuwMHDlBaWjoeHADr1q3D6XRy7Ngx1kTozEsRGI1GAw4H9rEx2TNC3FQWZGVRU1XFR//1X+hjY0nOzyc9K4v58+eHfZCoGh5dXV2kpqZecVt8fDxGo5Gurq6rPsfpdF4xrM1ut4e0RhF6ixcvpqO1lU/27+c7NpvsVCduGjqtloeLiugYHGSwp4fe5mZOaTScLinh0aefDutJhAH/K/3pT3+KRqO55s/Ro0cn/XoajeZLtymKctXbAX7+859js9nGf7KzswN9CyLMmM1mYqxWjBoN2gn+vwsxU8WYzRSlpnJrXh73lZby1bw8HPX1E36BDhcBtzyee+45NmzYcM3H5OXlTeq10tLSOHTo0BW39ff343a7v9QiueT555/nRz/60fjvdrtdAiRC9ff309DQQEdHB13797MmJWXCLw1C3CxsFgsatzvsl2oPODySkpJISkoKysFXrFjBSy+9RGdnJ+np6YC/E91kMnHLLbdc9TkmkynsrwUGyuv10tXVRUpKCoabaKvT7R9+yNDBg+i0WqK8XvIn+aVDiJlMq9WSq9Vy5uRJ5s+fr3Y5Ewppn0drayt9fX20trbi9XqpqKgAoLCwkJiYGO69917mzp3LN7/5TX7xi1/Q19fH3/3d3/HMM8/cNCOtqqqqOHXoEMP19STMn88jjz0WkcP2psLtcDA3NpYVBQVqlyJEWJmfkcGHlZVUV1czb948tcu5qpB+Sr3wwgu88sor47+XlZUBsGvXLu688050Oh1btmzhBz/4AatWrcJisfCNb3yDf/qnfwplWapzuVzs2LqVc21t0NJCgcvF8qQkdhw7RuOSJRQVFald4rS41EcmhLhSus1GYWcnJz777OYMj5dffnnCOR6X5OTk8NFHH4WyjLB0rqEBqqu5r6CAnIQEPF4vmp4efD6f2qVNi97eXobb28mYYOUBIQRow3jouoyJVIHRaKR06VIMBsMVy3FofD56e3uJ8LUqr8vj8fDZzp1Ye3vJlPAQ4qpSYmMZrq+nvb1d7VKuSsJDJfHx8bgv69vQ63Qstdmofu899uzcqWJloaUoCts/+YS+zz7jrrw8mdMhxARKMzMx2u1cuHBB7VKuSv7lquTS5Ea31zt+24KsLEp1OtrPnlWrrJAbHByk7fBhbktOJuUmGRQhxFQlAT3d3WqXcVUSHirwer1U79/PPLMZg053xX3JMTGMtLSwfevWmbu6sFaL+SYakizEVCVZLPS0tqpdxlVJeKigqakJT2cnJRfntlxudmoqd0VH07FlC59u2zajOtD7+vr407vvYuzpISU2Vu1yhAhrHq+XMbeboc5ORsNwn5ubY0JBGLhw4QIHysvpbWnBMzpKrtvt38f4KgpTUjDq9WzduZPtWi3r7r9/mqsNPpfLxZY33yS6upr7Z8+WlocQ13Cmq4uK8+exJyYSl5uLNgz7BiU8psnh/ftx7NjBrampGHU6Cq4zlyMnIYHlDgcHjh3DuXZtRM6qHxkZwWKxoNVqOXzwIK7qar46e7ZsMyvENbx7/Dg9NhuJt9/OV1atIisrKyznQ0l4TBOD0Qh6PfMD2DHMpNeDVhuWf3Gup66ujt1vvUVUUhJzFi+m7vRpZhkMEhxCXIPH66XH54NZs/izxx8P63/7Eh4hVlNTw9Hduxmz28kKcNmRVKsV7dmzfPDWW3ztG98I679IX1RfW0tCWxsp/f2crq5GB2QGaU00IWaqEZcLoqMpKi4O+3/vEh4hVFVVxf633qJoeJjk2FhyA9xm1WIwYPN6Ge7ru+Yy9eFIbzSiMZlYXVTEarWLESJC2CwWinU6Gk+eZHTpUqKiotQuaUISHiHg8XjY9qc/cW7vXkp8Pm6fM2dKr7Ojro6hoiLu+9rXwrLDbCJDQ0Ocq6kh8wvDkIUQ17ds1izqams5c+YMixcvVrucCUl4BJHH4+Hw4cM01tbiqKjgvvT0K5YfCZTd7Sa9uPiKbXrVNjAwgNlsvuYOZw0NDdDUxF0lJdNYmRAzg09R8LlcYT/PS8IjSJxOJ3v27KF5xw5yXS7mZWSQFR9/Q6/p1WqxhNHCaKdOneLg+++ji4ri3scem3ATrrbGRjK12i9NgBRCXN/+hgbIziYnJ0ftUq4pcq6FhLHh4WH27dtH844d5Hu9rJs794aDAyBeq2VE5T3aHQ4HAwMDnDlzhoPvvkvp8DDpra18tm0b3suWVgEYGxtj68cf03nyJDlBeP9C3Iw8Ph/4fGHd3wHS8piypqYmKk+eRFEUehsa8HR1Mcvn4+7i4qAdI0qvp9/hCNrrXY/X66WyshKfz8fs2bPRarW89etf4+ruBrebOW43K4uKGHQ4eOvUKT764APco6PoDQaW3X47dbW1dHz8McsTEpgzwTbCQohruyU3l5aeHrlsNRO1tray/be/xdDcTEJMDAtjYynMzsYa5DkM6TYbdRUVvDE6yp1f+QppaWlBff0vOnz4MJVvvYXR4+FUfj5aoxFdbS3rMjOxWK3jCxnaLBbuTU/nxO7dRAFDXi8fVlaCy8XKxERKA5jLIoS4ks1iAYeDwcFBUlJS1C5nQhIeUzA2NgYOB6sLCihITg7ZceakpWHq6aHi2DF2azRsePrpkBzn+LFj1J08ib2+nhWxscxJTaWirQ2fojB39mysV+l3yUlIGB8MoCgKzb29WKxW0my2kNQoxM3CoNNh8XoZGBhQu5RrkvCYgtmzZ3O0qIjO06dDGh4AeUlJaLVaPmluZmhoiNgQLCh4au9eYk+f5u7sbPIvvp+ls2ZN+vkajYZZMgFQiC/x+nzsrq+n0+UiVatFq9XS5XJh8fmYk5LCvMtGUro8HurOn6emrw+HzRb287okPAI0PDzMzq1bGa6pIXWa9qNIiY2F8+epqqpi+fLlQf9LFZeejq21dTw4hBDB0TM8TIPZTPodd9DT04NzdJSk3Fx8Xi/79u3D1t+PSa/ndGcn9W43vvR08h58kNtKS8NqiP7VSHgEaNuf/kTPrl2sTktj9jR1CpsNBpbFxHDonXc4e+wYOXPnsvK22zAajVN+zeHhYY4cPkxXYyMjAwMky7BaIYLOpyjg87Fs2bIv9V/8YWCAj/fvh/h4YoqLWbR4McXFxWE/yuoSCY8AzS4upqeqiphpXuV2YXY2aXY7bR0dnKyvxxoXR1lZ2ZRaIe3t7Wx55x3MtbXkmM2Y9HrKcnNDULUQN7eDLS2Y584l/ipD12+7+2465s4lNTWV7OzssL9M9UUSHgE4dPAgp3buJLq3l2QVRkGkWq2kWq3Ya2s5+sknnDx0iLIVK7BarRgMBgwGAxaLBds1Oq0HBgb4+Le/xdjSwl15eUGZjyKE+LJRl4tuo5HZJSUYrrJ/TWZmJpkRPDJRwmOS7HY7J7dtY8HoKIvnzsUY4Aq5wbQyPx9LaystbW1UnDmDS68HnQ6MRjRmM8WrV7Ni5Ur0V6mxoqKCqLY2NixciC6C1ssSItKY9XoYHb1qcMwEEh6T1NvbC0NDFGdmqhoc4O8DWVFQwIqLvzvdbhwXf9r7+6l47z3a29rIKyjg1ltvRavV0tDQwJFDhxhqbyffaJTgECLEfIoCZjPWaRpYM90kPCbJ4/GAzUbn4CBxYdahZTIYMBkMxOGfWJja10fVoUNUHT7MhfZ2XG43fadOQW8vWWYzdwRxFrwQ4ur0Oh3JHg+HDx2ioKCA6Am2nY5UGkVRFLWLuBF2ux2bzcbg4GBIEt7n81FZWcmhP/6RvIEB7igsxBQhzdDzdjvbm5vRA7emp5MVFxcxtQsxE5zu7OSz/n6YPZvb1qxh7ty5apd0XZP9TJWWxzV4vV527thB0/btzNfpWBFhS4ynWq1smDcPjUYjl6mEUEFJejrDTicn2tro6elRu5ygkvC4hj998AFdu3dzb2oqeRE6g1ov8zeEUNUtOTk0VFZy5uBB+vr6WLlyZVivWTVZ8nV0Ah0dHXScPMltiYkRGxxCCPVptVq+OncuOe3tXNi2jU8++EDtkoJCwmMCB8vLSevro2gGfEMQQqjLZDBw37x5LEtMZKytjTNnzqhd0g2T8LiKzs5OeqqqWJyZGVF7hwshwtuCrCxobqalqUntUm6YfDJeRVdXFya7XWZfCyGCatTlgrg4enp6iPCBrhIeV+NwOIiKsHVmhBDhL9pkolCnY6S52b8vUAST8PgCn8/Hufp64mSUkhAiBOamp0NfHx9/9BEVFRXY7Xa1S5oSGap7mZaWFk4cPcrgiRPclZ2tdjlCiBkozWbjbpeL4zt2cPjIEQ6npfH1p58m4eLOnJHipg8Pr9dLa2srJw4doufkSVJHRliXkUFSTIzapQkhZqj85GTyk5PZVF4OCQmYzWa1SwrYTR8eW//0J8599hlpDgcPZGWREcD2q0IIcUO0WrBaI2YDqMvd1OHh8/noaGrC2t/PQ0uWqF2OEOIm4vX5wGwmJkJX3b2pO8wVRcGg0ZAil6iEENNMA6DRYIvQKQE3dXjodDpy5s2jx+nE7fWqXY4Q4iai1WrB5aK9oSEi53zc1OEBUFJSgj0vj9dqamieYateCiHC29KsLOjsZP/+/WqXErCbPjzS0tJ4/Hvfw1dczOmuroj8BiCEiEzZ8fEwOEj1kSNs+uUvOXz4MKOjo2qXNSk3fXgAxMbGsuS222iLiaGtv1/tcoQQN4nEmBjuTE+Hqio4cYKK7dt57dVX/dteT5HL5eLgwYNBrPLqburRVpcrLi7mkMUifR9CiGlVlJpKQXIy5+12Pjp1CvLyqKqq4o477gj4tTZt2gRuN3i9pKamMiuEUw+k5XGR0WgkNjOTYx0dNPX0MBzh684IISKHTqslIy6Op5YvJ667m9pduzh44EBAr7Fp0yb6e3sZPX0ajMaQBgdIy2OcRqPhvq9+lR0WC9tPnoTubkr0eqKMRmLNZiwGA8mxsZhlD3AhRIhEGY18vayMXx08yKk9e6g8cYKypUvJysoiOTkZ3WVr7nm9XjZv3gwuFygKjs5OOvbvpz89ndumYSsJjRLhPcST3ax9sjwej38XwY4OTn36KUarFVdvLygK+oEB0oC8+HhK0tPRyMq7QogQcLhcHGluprG3F1dyMsTHQ2IiTz75JNHR0TQ3N7Pt7bfB4YCzZ0GnA7cbe0oK2qwsfvSjH0352JP9TJWWxxfo9XpycnLIyclhwYIFWCwWBgcHAWhsbKSro4O9R44w0NhIcVoasWYzWo0GnWwaJYQIEovRyOqiIlYDTT09bD9zBmw2du/cyfoHH6SzsxOam2Fg4LInWbAWFfHs9743LTVKeFzDpfVm4uLi6OvrIyoqiszsbIxmM1W7d1N1/jyMjWHw+SiNieHWvDx1CxZCzDizkpLI7+6mKT2duaWlDA8PU3v4MEVmM3euXs2oy8VrBw9CTo5/raxpIuExSX987TU8TU3gcqHPzydp8WIA+vv7cbtcnDh9moUeD0a9nFIhRHAlxcTQ6HIRGxvL9k8+wVVVxcqyMsB/iQuLBZKTefbZZ6etJvmkm6SU7Gz6amp4YuFCjrW0UN/SwihgVhRy9XqSbDYJDiFESGg1GvB4eHfzZhgdhfh4Pqqu5sHSUt45fhyKiqa11QESHpMWZbPRaTSi1WhYnp/Pcvyr8mo0Guk4F0KEVEl6OmOtrXgdDnyKQrXXS09GBu+cOAFmM6SkTGurAyQ8rmtkZISGhgbqy8u51Wq9omNcK53kQohpYNDpWHrZvI3Eri721tZi9/lg9uxpb3WATBK8rg/efpuDr78OTU0szMpSuxwhhKA4LQ2fzwcxMZCaOu2tDpDwuK7Fy5fDrFlgNuPx+dQuRwhxM7Ba/Z3g15OWpkqrA+Sy1XWlp6djMhiI1uulQ1wIEXqzZ8OaNf4///a3/omAXzDqcvknBk7zCKvLScvjOt589VWcVVUUJSaqXYoQYqaz2T4PDoAJLpW/dvAgpKeDisslSXhch29sDLq7SbNaqe3qom9kBJ/PR3VHB619fWqXJ4SIdFYrJCb6/xsdfeV9E+3todFARgbf+ta3Ql/fBEJ6Heall15iy5YtVFRUYDQaGbh8Kj1w8uRJ/s//+T/s3buXnp4e8vLy+N73vsdf//Vfh7KsgJStXk2V0cgf29v9/3M7O9E6HPhSUtB0dvJNWSxRCBGoy9eM2rDhS3d7fV50Wh0MDX3pvk3l5ZCcDGYzZrM5lFVeU0jDw+Vy8eijj7JixQr/6o9fcOzYMZKTk3nttdfIzs5m//79PPvss+h0Op577rlQljZpty5dStnixTQ0NADQ2tpKU0MDDAwQ5/NhuGyVSyGEuC6r9fPA2Lr1S3ePB8fWrWC3X/01MjNZsWJFCIu8vpCGx4svvgjAyy+/fNX7N27ceMXv+fn5HDhwgHfffTdswkNRFKqqqjjypz+h9PeD2028x8Oc+Hhmz5kjCyIKIQJzcQ08f0hc9vmxbRvce68/OACusqvppvJyf/hYrcyfP38aip1Y2A0fGhwcJCEhYcL7nU4nTqdz/Hf7RMkcJIcOHeLUe+8xH5iflYXZYEAvrQ0hxFRdXKUbnw+6uuCjj8DjgQsX4I03/J3gbvfErY6cnOmr9RrCKjwOHDjAm2++yZYtWyZ8zM9//vPxFk2oNTY2cmrHDhbrdCyRFXOFEMHQ0gJvvIHO7fYHxQMP+G9/442JA4OLrQ6bDRISVBuee7mAr7n89Kc/HV/PaaKfo0ePBlxIdXU1Dz/8MC+88AL33HPPhI97/vnnGRwcHP9pa2sL+FiTdba2FpqamJeREbJjCCFuQnb7FfM3vD7v5CYF5ueHsKjABNzyeO6559hwldEBl8sL8Ft6TU0Nd911F8888wz/63/9r2s+1mQyYTKZAnr9qUpKSaElPp4RlwuL0TgtxxRCzECXRld9sWVht8PWrf6+j/PnJ3z6pvJySEiA2NiwaHXAFMIjKSmJpKSkoBVQXV3NXXfdxbe//W1eeumloL1uMCQkJEBSEu8eP86zq1erXY4QItwlJsLXvgbvvAO9vf7bMjI+vzT13nvQ3X1lmLS0TO618/Ku2R883ULa59Ha2kpfXx+tra14vV4qKioAKCwsJCYmhurqatasWcO9997Lj370I7q6ugDQ6XQkJyeHsrRJycvLA7OZgjCoRQgRAS5dFbn0X6v18+AA+OpXYft2uHRp/jr9HHCx1ZGSAjExfP3rXw9B0VMT0nGmL7zwAmVlZfzkJz9heHiYsrIyysrKxvtE3nrrLbq7u3n99ddJT08f/7n11ltDWdaktba2Qm8vOWGU9kKIMHJx2Oy44WHYsgU6Ovy/XzaBuLyl3P+Hy/t0rzPJT1EU/xpW+fk8+OCDwao6KEIaHi+//DKKonzp58477wT8ne9Xu7+5uTmUZU1ae3s7jI3hVRS1SxFChJtLk/02bPg8RDZsgPXrPw+Uiy0Qr8/L6tyrXPr2eq95iP/67DPIzQWjkfT09GC/gxsiM9yuYcGCBZCfT/2la5dCCHHJ5csSxceP/9HrcfvnaYB/DSrwT/w7derLr3HpcVexqbzcv9ZVRkbYdJJfLqzmeYQbrVYLo6OMXGVJZCHETe7yD/516/yT/d54A53Z7A+WhAR/K2TXLhgZ+bzvY9cu/x4cPt91+zsoLFRtv47rCc+qwoTFYkEfFcXg6CgXQjyTXQgRYex22Lnz898feACKi+GRR/yXry7t/9PXd+XS6iYT3HGHf+l1m+2qL72pvNy/jInNFpatDpDwuCaNRkNuQQGUlPDH2lo817k+KYS4SVgs/n6NLy6ZvmjR53/WauFPf4KCgitvv9SHeuTI50uVXGZ8TkdOTtgGB8hlq+sqKS2l4dAhGBvD7fXKulZC3OwsFvjmN6//uKuslOH1etA1NsK5c1cNjs7BQf8lr6Iili5dGoxqQ0ZaHteRlpbGkvXrYf58fltZyfDYmNolCSHUdLX9e76wV9EX1Vyo8f/hUqvjKsEB8OHJk/6WitHIostbK2FIwuM6tFotixcvZsHq1RATQ9Wl8dtCiJvTF/s6AA4d8i+pDhxuP/ylp8xNmQuATm+YcOvYTeXl/hnqKSlhfbnqErlsNUkWiwWSkzk1MEDShQsUpqSoXZIQQi1f3Gvj0mgrYGnmVS43bd3qH501PHzVEVavHDjgnzBYVERZWVkoKg46aXlM0sKFC/nO3/4tUQUFdNntON1uzvX343C51C5NCDHdrjY/o7//yzsD7trlX4KkpcU/63yCUZtOtxvmzAGDIWxW2LgeCY8AGAwGRvv6qO3q4nd1dXw8MkL52bNqlyWEmG4XV8O9gsHgb1lc5PV6/R3jk1m7KicnrIflXo1ctgpQ8bJlnLFY0NtsMDqKo69P7ZKEENPt8vWsLndZi0Sn0/lX1G1omPBlGrq7/ZerwnxY7tVIeATA5XLR39MDvb24GhowKQpFQVyeXggR5mw2/4ipL+5ptG+f/792u39hxPXr/b+vXetfgn2C1senp09DUVHYziK/lsirWEXV1dWc37+f1TodG0tK+HZpKXNll0EhZqxN5eX8+44dDIyOwuzZ8Pjj/tYE8OsTv/Y/aPduWLXq8wUS29uhshLwz+uYaP2qTeXl/jkjETK66ouk5RGA2upqYux2iouL1S5FCBFiPp+PEaeT5u5utrW28tgPf+i/4+L+PhvLNvp/dzo/f9KlYbgHDkBrK7q+viu2m71kU3m5fzHFuXMjstUBEh6TNjY2hv3cOcxeL4qioLm4WqYQYmbYVO7fb+PSrqHlZ88yqNdztqeHWy/fEK6k5PM/b9165bDdy1sZ7e1XPY7X5/P/ITeX+YsWsWLFiqDUP90kPCbJZDIRl5/PQH8/F4aGSJ2ow0wIEXE2lZf7O64v+/CvO38eT3Q0eStX8tTzz1/9iRqNvz/jjTf8v19nZJXX52Pz3r3+fg6rNWKDA6TPY9KGhoYYaGhA29VFUkyM2uUIIYLE7fX6g8NsHt9/w+XxgMHA7GXL+M9Nm8YfW/Xpp1c++d57/f0WdvukhuRu3rvXv8x6WlpE9nNcTsJjkkwmE7q4OHxmM3LBSoiZ4zf79vk7uk0m8HgAeHn/fkhNJeYLW1C/ezFIvD4v/3743/03RkVN7kAaDcyaFbabOwVKwmOSTCYTmTk5aLm4SZQQYuZIT/e3IC53eUc48K2nnuL3n37Kb/7jP9BpdSxMXQjAQbd7vL9kIuPrVmVlzYjgAAmPgNgHB0mWJdmFmHm6u+l3Oj9fNVujgfR03Jf1gRiB1WvW8J2//EsAbs+9HQC3Xg8WC8dbW6/60uPbyRYVsf7SboIzgHSYB8DncuGSDaGEmFn0ekYHB6msrQWX6/MP+7g4NlycDLh7925KU1OJ/8Kk4P/6l3+htbaWbIeDxTk5X3rpl/fv9w/fnT8f9HoyMzOn5S1NBwmPAEQlJNBlNOJ0uzFNsKyyECJyeH0+hrVajldXU5Caii0/3z/iatasKx53yy23UHf4ME78g2daWlr4zW9+w/naWlbHxVGclval1+4bGcGlKDBvHhiNM+Zy1SUSHpPU399PV1UVaT6fBIcQM4TD5cI9NARuN1GKQoxGAzEx/n3H09L49a9/zcaNG4mNjeWJ738fn89HbGwspaWlPPHEExzftQtOnmR1UdEVr7upvNw/+a+4OOIWPJwsCY9J0mq14PNd8Q1jaGyMc/39JMXEkBwbq2J1QoipiDGbyYiNZcMjj7C1uhq6usDnA6OR4eZmOuLjxx+7bds2vF4vjz32GADHjx/3h8wXjLpcoNPBggUQGzsjgwMkPCatvb0d+vupczoZdjo5OzjIoMkECQlom5r4ZkmJtEiEiEDfXL58/M9Rej1PLV/OpvJy3BkZV/6brq9nIC4OgE8++cS//PrIyBWv5fP5eO3gQf8kwBkcHCDhMWmDg4Ngs9GfkcEFu520ZcuYm5dHY2Mj/V1d6GUUlhAR7dKyJABotcRnZPDst789fpP7wgXwemlubqa/v9/fQunrY+OqVf7LVDqd/5JXQcGMmAR4PRIek7Ro0SIKCgrQarW0trYyMjzMoY8/RtvRwYrERHQy90OIGWF7TQ1YLPxp1y42//73/MXXvkZUdzf28+chNRWDwcDQ0JB/qXXg1/v2+edwzJ4NRiPAjA8OkPCYNIvFQmtrK3vefht6etC63cyPjaW0sJBok0nt8oQQQbL34uWp48eOYTEa6T58mHizGfLzIS7OP9zW64XRUf/EwsREmDWLO9esoegLHeczmYTHJHk8Hg7u2UN0UxOPLFyIQafDqJfTJ8RM4XS7+Y/du9lWXc3Y2BgLMzLQe73+7WQzMiAzc7xFYY6OZmz+/PHn/vmf//lNt/KEfPpNUnV1Na7aWu4vKZGWhhAz0I4zZ/DpdPgMBnR2O0k+H/G5uSTdfjvExPCNb3xj/LHf+ta3VKw0PEh4TILdbufQn/5EhtMpQ3KFmKHWz59Pe3k5hRkZ2BITmVNSgmXRIjCbb4o+jEDdXO2sKYqOjiY2NZU+rxffpY1chBAzkiYhAW10NJYFCyQ4rkHCYxJ0Oh2G6GjGdDrGLi7ZLISYWcY3hAIG0tIgKkqC4xrkstUk9XV3w/AwUReH4gkhZoaW3l7/7PKYGFiwgNV6vYTGJEh4TJIOSL84u1QIMTOMtzby8yEzEzQaCY5JkvCYhObmZrznzqFcYzn2nuFhXB4PGRIwQkQOrda/6m10tIRGgCQ8JsFkMkFuLj1VVV+6r7mnhxMdHXTHxICicNvoKHMzMlSoUggRiOaeHoiPl+CYIgmPSUhPTyfKYGDUbsfl8dA3MkKX3U7r4CBdMTHELl9O2Zw5nNixA+ULW1cKIcLT4eZmSEtDo9GoXUpEkvCYpPlLlnCorY0/HD2KIz4eUlPJuOUW7lu8mLaWFk5s2UJiby/5c+aoXaoQ4hou7Te+sqCA/Xo9Tz75pMoVRSYJj0lauHAh5zs6uJCezi1LllBWVsbAwACVp05Ru3UrZRoNt8ybd9MtUSBEJPH5fP5tYbVaPm1rw2218tprrwE3x2KGwSThEYB7v/IVANra2njjN79huLkZ7eAgy+PiWJCVpXJ1Qojr+dXevZCcDElJnDh5Ep/DwS1WK5aoKDo7O0lPT1e7xIghX5MDdO7cOf70u9+hP3yYB2Jj+eacORIcQkSSpCQ4fZrFikIZYGlshLExdu3apXZlEUVaHgHq7e0Fp5OUmBgZlitEpDEYGPL56Dx7lp88/jixMTHsPHOG+r4+hoeH1a4uokjLI0AlJSVgs1E3OKh2KUKIANU3N3PiwAEAPG63/7YLF0C+CAZMwiNALpcL3G5yoqLULkUIEYBN5eVkzp1LWl8febm5xMfH++/Q6yE5WTrMAySXrQJ09uxZaGtjRUGB2qUIIQJk8fkoKiwc36+8obvb3wcie/QETFoeAYqKigKzmY/PnmV7bS3uayxZIoQIM19YFfvTM2ckOKZIwiNAhYWFLHn8cWLvu48mm43ari61SxJCTNbY2JW/Kwp0dsLgIEeOHFGnpggl4REgnU5HVFQUF+rqiO/upjAlRe2ShBCT8Ozq1SzNyxu/ZHXpNlwuGBjgxIkTKlYXeaTPI0C1tbWU//73zBkZYUVJCUa9nEIhIsWi7OzxP+9vaOB4czMagwFDTw/k5qpYWeSRlkcAXC4XB7ZtI/rcOe4oKpLgECJC7Kmr448VFQyNjY2vbXW8uZma8+c57HDgkuAImHz6BUCr1eLS6XAZDGqXIoSYJI/Xy65TpxgYG+OC3Q5aLXaHA0wmhuPimFtUxAPf+AZ5eXlqlxpRJDwCcO7cObhwARwOvD4fOlkEUYiw9+t9++iy23FYrYy4XBj0et44cgTjrFnckpTE008/jc1mU7vMiCOffgGIjo6GjAzIzaWpp0ftcoQQk+ByuYjRaBjt6+O1vXuprK7G5XZDdDT/7b/9NwmOKZLwCMDO7duhooKikRGyLs1OFUKEtR/cdRd3zJ9Pnl5PvsPB2KXNn1wufvWrX6lbXAST8AiALSEBYmIw6fWYpd9DiIig1Wp5cMkSfrFxIwuXLuWW1asxxsXBhQs0NzezadMmFEVRu8yII+ERgNWrV5O6ahWVo6N4ZGa5EBHF4/WSkpyMeWgIbDbIy6O1tZX29na1S4tIIe0wf+mll9iyZQsVFRUYjUYGBgYmfGxvby8LFy6kvb2d/v5+4sJwlUuj0chQVxeFRiN6nU7tcoQQk7SpvByfXs/ugwfxajRoAcOyZaxetQocDt544w2eeOIJtcuMKCFtebhcLh599FG+//3vX/ex3/3ud1mwYEEoy7lhPT09jDY3Mz8jQ+1ShBABUpxOxvR66hSF9thYjEND+I4fh/Z2hoaG1C4v4oQ0PF588UX+5m/+hvnz51/zcf/v//0/BgYG+Lu/+7tQlnPDhoaGwOuld2RE7VKEEAF4dvVqdHo9t956K2aLhfScHBSjEa3H41/fSgRM9T6Pmpoafvazn/Hqq6+iDfN5E7m5ueStWUP52BjHW1vVLkcIEQhFwRwVxcq4OHIcDuYUF/v3M4+JUbuyiKTqp7XT6eSJJ57gF7/4BTk5OZN+jt1uv+JnuhiNRu79yldIWrSIDtlJUIiIE+vzUZCejsXh4OSxY4y1tIDVKhtBTUHA4fHTn/4UjUZzzZ+jR49O6rWef/55SkpKeOqppyZ9/J///OfYbLbxn+zLFjqbDqOjo/R2d6NIh7kQEWXjxc5xo9HIkNPJhXPnqDtzBsxmtUuLSAGPtnruuefYsGHDNR8z2TVidu7cSWVlJW+//TbA+FjrpKQkfvzjH/Piiy9+6TnPP/88P/rRj8Z/t9vt0xog5Xv2oLS0UJSQMG3HFELcOL1ON74c+78MD6Pp62NWbi4YjSpXFpkCDo+kpCSSkpKCcvB33nkHh8Mx/vuRI0fYuHEjn332GQUTbPNqMpkwqbTz1+DgIK1nzoDTSVJamio1CCFu3I/uv583jx5lwGJRu5SIFdJ5Hq2trfT19dHa2orX66WiogLw78YXExPzpYDoubheVElJSVjO89Dr9f5vKfHxfNjQwCOzZxMXFaV2WUKIKRgYHYXsbDSXlisRAQlpeLzwwgu88sor47+XlZUBsGvXLu68885QHjokoqOjefYv/5Kuri4++O1v2VJbS7TBQGZMDLfKcs5CRJ74eJ555hm1q4hIGiXCF3Wx2+3YbDYGBwexWq3TdtzOzk5qa2vpOHeO4WPHmG80kpuYSEYYtpiEEF+26fBhWLpURlp9wWQ/U2U/jylKT08nPT0du93OeyMjVHZ2UtnWxtMxMTfVDoM+nw8g7OfoCHG53uFhkC96N+Tm+ZQLEavVyteefJKjR49S98EH7GtoYF56OinT2AqaLl6fj121tQy5XOiArsFBMJuJAu4qKsJmseDyeIiPjla7VCGu6eyFCyDD7W+IhEcQxMTEsGrVKsbGxjh7/DhnW1q4xWrllhmyL/LA6CjHW1up7+2F3FxITweHw7+sQ1QUo0NDfHRxh0W0WhgaAo8H9HrQ60k2GHho4UK0F+cBCaG2KKPR/3dUTJmER5AYDAbuu+8+zpWWsmfHDo5VVZE3PExiBC99MDQ2RlV7OzWjo3hnzYLCQkrmz+f2228H/LP93W43v3v9df8lALcbDAYYG4PhYf+y1xoN3Z2dbD52DHw+ohQFjUbDnUVFZMqGWkIluYmJHKyqAhnoMmXSYR4CXq+X3/x//x/5HR2sKijAFEEbR3UODrKtpganXg8mE6Smkl1Wxt3r1mGY4H243e4r7lMUBbfbjV6vR6vV8uabb/qX4/f5/K0ShwM6OlieksKCrKxpemdCfG5TeTlkZkJBgXSYf4F0mKtIp9Mx//bbObllC+01NTy1YEHYX67x+Xx8euYMTW43FBcTk51NUlISixcvvu6k0C+GikajwXjZrN3HHnts/M+jo6M0Nzezd9cuKiorJTyEehISePLJJ9WuImJJeASZ2+1Gq9WydOlSTh46hLGnJ+yDo3NwkA9PnoTiYmwFBaz7yldCNkkzKirKP9ny3DlWTfO6ZEKM0+nAZiNaBndMmYRHkL312msMO53kz5oFFy5wS5hvHLXj9Gka+/th6VIWLF3K8uXLQ37Mhro6MkdGKJg1K+THEuKLvD6ffxl2GV5+QyQ8gmx4bAzq62mvqmJJXBz5QVoHLNg8Xi976+tp9Hhg+XJW3nYbpaWloT+ux0PX2bMsCpP+KXHz2VpdLXt4BIGER5AZXC6sIyN8ddGisJw41zcywpmuLqrOn/ePNMnI4KmnniJqmtboqqqqwtvaymxpdQiVnOvvh+Rk1q5dq3YpEU3CI8gS8vI439CAy+vFHGbhUdXezv7z5yEtDVas4Lbbb2fu3LnTdnyn00lFeTklBgMxsoeCUFNc3IQrd4vJkfAIspSUFM5HR1Pb1cXCMOsQ3t/QACUlFCxfrsq3rsbGRjzNzZTl50/7sYW4QgQNnw9X4fXVeAaw2WxgsVBx7pzapVzB7nCA1Qrx8dx2222q1NDZ2YlvaMg/u1cIFSiKAhaLLE0SBNLyCLJLk2pWhNG369OdnXxWXw9FRVhiY1XbTKu+rs4/C10IlfzXZ59BmI+AjBQSHkE2NjYGDgemMBhNVNHWxuGmJkhNhWXLuHXVKhYuXKheQVotWK2Mud2Y5bKBUIvNxooVK9SuIuJJeASRz+ejra0NDAYcKn/DfvvYMfpiY2HJEmbNm8ddd92FLhya6sPDeLzeG77m7PX5GHW5aOvro290lGGXi9b+fpbn5pIVH098VFTYT84UKomJYf78+WpXEfEkPIKora2Ns7t2QWcn2QsWqFJD/YUL7Gxp8bc2MjP59saNql2m+qLFixdzvK+P39XXk6XRkGCxkBkXR0J0NNGTrNHn89ExOMjetjbsBgMkJWEtLkZrMMDwMAf7+6G9Hfr6SDAaSYiKIichgTSrVUZ4Cb9w+BI1A0h4BFFMTAzExJAfGzvpD8Ng29ncDIsWsf6rXyUjIyOsvn0vWbKE3Nxczp8/T0tjIzVNTZw6fx5aW1lgtXJrXh66CYY3n7fbOdHaynnAabUSc8stlJWWUlhYSPxlq/M6nU6qq6s5+tln9Gk09Hk81Hd3w8AAtrEx1hUWyr7zQgSBhEcQJSYmkpiXR9fZsyE/lqIoHGhspKqzE3Q6FqWl4fH5IDOTB/7sz8gI007B5ORkkpOTKS0txefzUVNTw/6PP+bU0BCVhw6xKjcXs8FAtMlE6sV+o67BQT5obITcXIpXrqSkpITExMSrTsI0mUwsXryY+fPnYzAYGBsbo6enh8HBQfZt28abR4/y9MqVN9Vuj+ILLu5+KW6M/AsKMo2iYA3xUFSvz8fHlZV0XuwIN0dHU+9yMdLfz6I1a8I2OL5Iq9VSWlrqnxtz/jx11dXsb2zENzaGtrOTlUlJ5CQk8EFtLcyfzxMbNxIbGzup17600q/ZbCYrK4usrCy8Xi8HHQ5ePnmSeTYbqwoLQ/n2RLhyOtm0aRMbN270L9IppkTOXBApisLQ4CDZIbxk5fJ4ePngQSgsZMkDD7B48eLx+3w+X1guiXI9KSkppKSkMH/+fIaHh7Hb7ZyqqGDvnj1w4gRkZJBZWDjp4JjIggULsFgs7Nqxg+qzZ0nq6mJOWlqQ3oWIGPX1MGcO586dI082g5oyCY8gcblcvPeHP+A8c4Y5qakhOca5/n4+rq6GefNYvHbtFcEBRGRwfFFMTAwxMTEkJCTw6pkzkJ3NratWBW10zOzZsykoKGDnp5+yZ/t2ehsaKEpNJUkWyrspPLt6tX8jqMFBOjo6JDxugIRHkHg8HuwXLhBtt5NZXBz0169oa+NwWxuUlbFw5UqWLFkS9GOEE7PZzNPPPYdOpwv6EGOtVsvd99xDZVoah7dto6qzk1XR0cyLkMt94sZkxMXRoSjX3eRMXFvkf1UNE1FRUeQtWMCIVkvP8HBQX/tkWxuHh4dh2TLWP/YYy5YtC+rrhyuj0RjSuSnz58/n23/91+StXcu+oSE2HTwY9P93Ivw43G7w+Th+/LjapUQ0aXkE0fLly+nv7mbXZ5/xaFnZDb9e7/AwFW1tNHq9MGcOT3/3u1ds7ypunF6vZ+0993A0Pp6T+/bxbnMzeVotxWlpmPV6UsJgpQARXOk2G/0uF3a7Xe1SIpqERxDFxsaSlZdHVUUFFW1tLApwVd1Pz5yh1eUiGRgaG2NIq4XiYmYvXMjqO+8MjxniM5BOp2PZsmXMmzeP3/3udzS3tdHc1wceD4t6e1mSmzsj+pOE3/JZs6iprFS7jIgn4RFkK1euxOl0cvSPf2RBAKOfXjlwAGdODlGzZ9PR2wtjY5gSE1n/yCNybXaaxMTE8Oyzz9LU1MTw8DCd7e1U7N1LR1UVd82ejdViUbtEEQR6nQ7sdnC51C4lokl4hIDi8/nnKkwiOJxuN8dbW3GmpmKbO5fHH38c8C+waJblNFQx6+Iuh/PmzeNMTg5Hdu3i7dpa1qSlkWGzYZJFHWeGwUFGR0enbRfNmUba4iHgUxSSJvEB4/P5eOXoUSqB5CVLePTRR8fvk+BQn1arZe7cuTz5F39Byl13sd3h4HenTqldlggWr9e/CraYEgmPIDt79ixNhw6RMYl5A4MOByQkYCkt5cFHHpHr6mFKr9ez/uGHybv9dtwaDRVtbdRfuIDT7cYnS11ENEVR1C4hYsllqyAaGRlhzx//SGFfH8uKiq752D8cOcKgVgslJWzYsEGWSQhzGo2GxYsX4/N4OHziBIyOQns72QYDdxQV4fF6aenro31gAIfPx9rCwkn1kThcLtr6+8lLTJT1tqZRWU4OJxRFWvg3QP62BomiKGx5/318LS0sKy7G7fXSNzJCjMlEjNnMR6dO0eHxgEbj300vPR2ys3nooYfG12ES4S0pKYn77r8fx5o11NbWUnXqFJ29vbxWV+dfbC89HTIzQatlX0MDd5eUYJhghFxbXx+ftbczrNdDXByJp0/zlTlzZIveaTLqckGYbFUQqSQ8gsTpdDLQ2MjcqChGnE62NDXhio9HOzyMaXAQR0oK5OdTUFCATqfDarVSVlYWVkumi8mxWCwsWrSIRYsW4XQ6+eyzzzAZjSxdtgyj0cjp06fZ9/77bD9zhvvnzWPU5eKC3Y5BpyPNZqOqvZ2TfX34ysq44/bbiY6OZveWLXxQWcn6oiJi5dtwyCmKAoqC0+kkOjpa7XIikoTHDfJ4PNTV1XHq6FEYGiIlNpat1dW4ysr4+pNP0tLSQn11Nelpaaxdu1bCYoYxmUzcfffdV9w2d+5cjEYjO3/9azYdPgxxcWC1wvAwtLRAUhI569dzx5o1WC5e2nr4G99gyzvv8P6xY9yfn0+CfKCFlO9iX4fMnZo6CY8bVFdXx94334TmZuJNJnaPjUFmJqW33EJCQgIJCQmUBWG2uYgshYWF1K9ahc/rpWT+fJKTkxkcHOT8+fMkJSWRk5NzxeNjY2N5+PHH+dhk4oP9+/lKbu74fiYi+OovXIC4OLxer9qlRCwJjxuUn5/P3oQEYs+do9/jwbpsGctWrCA3N1ft0oTK7lu//orfY2JiyMzMnPDxFouFB7/2NbaazWzZvZt70tPJTkgIdZk3NRmqO3UyNvQGmc1mFixdylBKCsyaxS233sqsWbNk2K2YEqPRyFcefJDMdevYeuGC/xuyCDqDTgdOJw0NDWqXErGk5REEy5cvJyoqiqioKApldzpxg/R6Pffcdx/lFgs7t2zB6fHIcvFBdraxkf7z53GkpHD77berXU5EkvAIkgULFqhdgphBtFotd9x5J2aLhX3vv4+juZklsnFR0LhcLoiKwul0ql1KxJJrK0KEKY1Gw/Lly1n6+OMcNxgor6uTGe1B8oOHHqIsJ4e4uDi1S4lYEh5ChLlFixZx55NPUhsfz/YzZ/DICKEblp+c7F9V1+1Wu5SIJeEhRAQoKirivm99i/acHD4+fRqXx6N2SZHP64XhYRwOh9qVRCQJDyEiRHZ2Ng889RT9s2fz4enTOGQ/ihs3OMhvf/tbtauISBIeQkSQlJQUHnzySRylpXxQW8uwzFO4MaOjalcQsSQ8hIgwCQkJPPTEE/jKyni/rk4uYU3RXcXF/h0FZRDClEh4CBGBrFYrax94gJH4eHpHRtQuJyIVpqT4O82l9TYlEh5CRKi4uDgwmfzLi4uADY2NgU4nLY8pkvAQIkIZjUb0UVESHlMUazb7R1zJyrpTIuEhRASLjotjRGZJT51G4w8QETAJDyEiWFRcXEAtD5/Px6GmJvbU1dEzPAxAe38/22tq6Lus7+SmmcmuKDA2dvO83yCS8BAigkXZbIwG8M259vx5ThqN1CYk8G59Pc09PWyprqYpP58/NjRQ1d7Om5WVvFJZSe/FcJnxnE5+9atfqV1FxJHwECKCRUdHc72xVkNjY3QMDHDq3DkOdXWRsWgRG557Du3s2WwbHkZTVsaqtWvxZGSwv6GBgTlz0C9ezHs1NVS2tzM205fwkMt+UyKr6goRwcxmM4MeD4eamogxmegdHibKaGRBVhYen4/ari5ODAzgiYlBm5hI4rx53LVuHVFRUWx45hlGR0dJTExEp9MxZLfTlp1N2YoVxMfHszMqigMnTnDq1Cm+On8+UUbj+HF9Pt/M2bNmYADsdoaHh4mJiVG7moihUZSLm/lGKLvdjs1mY3BwEKts2yluMufPn+f93/8ek8GAs6+PmIwMRvv60J47h2dgAE1pKbmlpcxbsIDU1FT0+sl/X3Q4HBw+fJjagwehspJvr1hBQ3c3zb29dLhc3J6ZyZy0tBC+u+mxqbwccnIgL49nn31W7XJUN9nPVGl5CBHBUlNTefaHPwT8H/Zms5nR0VFqa2upr6lh6W23kTfFfUAsFgt33HEHY2NjtIyM8MqBA5CWBrNmQXQ0e2pqSIiOJjk2FgBFURh0OBhxOrEYjSRERwfpXYbWY0uW8GZFBWRlqV1KRJGWhxDiuvbs2kX9yZOkzprF3evWodFoeOU//gOam8Hj+Xy+hMkEMTHQ1cVTCxdecakrnG0qL4e8PMjJuelbH5P9TJXwEEJMiqIoaDSa8d/tdju7du1ibGwMi8VCWloaYxeHvdbt3w81NTy1fHlEBEjn4CBvVlVxxOvll7/8JUlJSWqXpBq5bCWECKrLgwP862s9/PDDV32sY3SUNuCNU6f4VlkZXp8PnVaLPkxncydEReFzuUCrla1pJymkwyVeeuklVq5cSVRU1DW3e3z55ZdZsGABZrOZtLQ0nnvuuVCWJYQIsa/cfz9ppaV4YmP5dWUlr5w6xa8PHAjb2fCvHDhAdGEhS269lczMTLXLiQghbXm4XC4effRRVqxYwebNm6/6mH/5l3/hn//5n/nFL37BsmXLGBsbo7GxMZRlCSGmwcKFC+nq6iI+Pp7+3l44e5aT586xsqBA7dK+TKeD9HR++IMfqF1JxJiWPo+XX36ZH/7whwwMDFxxe39/P5mZmXz44YesXbt2Sq8tfR5ChL/Ozk4+fPVVOHOGtYWFFCQnq13SuE3l5f5RZEVFN31nOURIn8f27dvx+Xy0t7dTUlLC0NAQK1eu5J//+Z/Jzs6+6nOcTucV1yTtdvt0lSuEmKL09HQWr13LcZOJT+vr+bS2Fm1UFElGI3cXFhJjNqtbYFLSlL/A3qxUnSLa2NiIz+fjH//xH/nlL3/J22+/TV9fH/fccw+uCRZ7+/nPf47NZhv/mShkhBDhZcmSJTz29NNQWgpLluBbtIgLOh0DDoe6hRkMEBdHQTheTgtjAYfHT3/6UzQazTV/jh49OqnX8vl8uN1u/vVf/5V169axfPlyfv/733P27Fl27dp11ec8//zzDA4Ojv+0tbUF+haEECqJi4vj2e99j2f/6q/4yvr1kJvLnsZG1bbS3VReDunpMFOWWplGAV+2eu6559iwYcM1HzPZGa3p6ekAzJ07d/y25ORkkpKSaG1tvepzTCYTJpNpcsUKIcJWdnY29/3Zn7HtD39gd10d9172OTCtkpOlr2MKAg6PpKSkoE2gWbVqFQC1tbVkXVwaoK+vj56eHnJzc4NyDCFE+MrJySF/8WLq29uxOxxYLZbpLcBshqio6T3mDBHSDvPW1lb6+vpobW3F6/VSUVEBQGFhITExMRQVFfHwww/z13/912zatAmr1crzzz9PcXExa9asCWVpQogw0N3dTf3p0+DzocpSFyaTfzdBEbCQhscLL7zAK6+8Mv57WVkZALt27eLOO+8E4NVXX+Vv/uZvWL9+PVqtljvuuINPPvkEg8EQytKEEGFAq9VCfz95Wi22aW51jLndIPu/T5msbSWEUM2ZM2cof/ttVrrdlE7zzO5N5eVQWAgZGdLncZnJfqbKEAMhhCp6enoof+89rK2tFKakTOuxN5WXQ1wcpKRIcEyRLIwohFBFVVUVNDXxyMKFmKfxMvXOM2f8neTz5jF/0aJpO+5MIy0PIcS083q9tLS0QHw8vcPD03rs+gsXIDERdDpWrFgxrceeSaTlIYSYVhUVFRzeswfa2pir1ZIyzX2Vs1NTORumS8NHEml5CCGmTUtLC4e3boXjx5mn03FbYSGGaf4gL0hOhvPnweeb1uPONNLyEEJMm8oTJ0hobubh5cunPTQu+aSqCoxGWZLkBkl4CCFCbnR0FEVR6DpzhqXJyaoFh9fnA70epKP8hkl4CCFCqre3l3d+9zv/h/a5c8wuLlaljlGXi9cOHoSCAjCbZYjuDZLwEEKEVHl5OdTXQ1wcc8xmLEajKnW8dvCgf8dAm02CIwjkop8QImTq6uroPnKEPI2GpT4ftxcWqlbLs6tXg9fr7+8QN0xaHkKIkKivr2f3228zx+XiDrWWW7/MmNsNMTESHkEi4SGECKqzZ8+i1+tpbGyEjg5uKSpSuyQAXj1wACa515C4PgkPIURQjI6O8oeXX8bd2QkeD1itEB3NoMOh/h7llxiNPPDAA2pXMSNIeAghgqK+vh53XR3LYmNJjIujta+PUZeLdJtN7dI+Z7GQkZGhdhUzgoSHECIoLpw/T5LHw8LsbACy4uNVrugqZFmSoJHRVkKIG+Z0Omk/fZqs2Fi1S5mYXg/Tvc3tDCbhIYS4Yd3d3Tjb28mMi1O7lKvy+Xz+Zdil5RE0Eh5CiBtmMBjAbMbt9apdylX9au9eSEhQu4wZRcJDCHHDzp49i76/n7Rw6hy/6ILdDmYzZGXJzPIgkvAQQtyQ/v5+asrLMQ4PU3/hAsNjY2qXdIU/VlRAcrKsohtkcjaFEDdEo9GA1cpoXh77tVo+qq1Vu6QvM5lYJCvpBpUM1RVC3JC4uDj+/Ac/wOFwUFFRQfV77zE0NkZsuEwMBFAUBgcH1a5iRpHwEEJMms/n4/jx47jdbnJycvD5fJytq2NeaSlJSUnUnDwJIyMoiqJ2qVdyu2lqalK7ihlFwkMIMWnt7e0cf/ddGBykMi3NvwxJVxf1FRX+YbB9fSxPTsYaRvMpnlq+nNeOHoWUFJqbm8mT9a2CQsJDCDFpDfX1xAwN8eCcOTR0dxNjMpFYXMyhlhbaenuZl5FB6axZapd5hdcOHoT4eLBYyL44+13cOAkPIcSkjI2N0VRRwUKbjVizmUWXfRDfFwZLrl/TrFksXLQInUwSDBoZbSWEuC6n08nbr76KvrGRopQUtcsJjFYLJhPLli1Tu5IZRcJDCHFd586dY7Sujq/MmhU+y6tPwt76ekhJAYNB7VJmHAkPIcR1tZ87R4zTSWJ0tNqlBKSmo0N2DgwRCQ8hxHW11tRg4+KEwEgjrY6QkPAQQlyXNSMDeyQGB0B0tKxpFQISHkKIa1IUBZ/LhTXCwmNTeTnExkJMjNqlzEgSHkKIaxoYGOBCZSVFSUlqlxK4vDz/JlAi6CQ8hBDXFB0djTY+niNdXf5NlSKJySSXrEJEwkMIcU1Go5H0oiLcRGCHuSzDHjJyZoUQ1zXQ1YXR45HwEOPkYqAQ4rpGvV4Urxe314shUpb40OlkjkcISSwLIa7pxIkTKN3dxEfaB7F0lIeUhIcQYkJ2u50j27bB6dN8bcGCyGl1ALjdEG77iswgEh5CiAnFxMSAwUCK1Yo20voPfD4YHVW7ihkrwv42CCGmU1dXF4yNYTEYwm93wMkYGaG1tVXtKmYkCQ8hxIS6urrA4aDFbGbb6dMMOhxqlxQYu51PPvlE7SpmJOlREkJMaPHixaSlpdHb28uBrVtpqa3l8TlzsIXRNrMTeXb1ajYdPYrJZFK7lBlJWh5CiGvKyMhg/vz5PPn972MqLeW96moGIqUvQVFwOp1qVzEjSXgIISYlOjqatQ89hCs7my67Xe1yhMokPIQQk5aUlITBZqOlr0/tUiYnEjv5I4SEhxBi0sxmM8vuvpsWvZ5NFRVsq6lRu6Rri7SFHCOIhIcQIiDFxcWUfvWrUFpKc18fp86do2twUO2yri7S5qZEEDmzQoiAaLVaVq5axfxFi6CwkINOJx80NdEajpeyxsbAbmfTpk1qVzLjyFBdIcSUrFi5kty8PI4fPkxHRwfbDh3iiSVLiFZxaOym8nIA9DodG1et8t/Y0AALF6IoSuStChzGpOUhhJiyjIwMHnjkER545BFYuJADjY1qlwR6PR6vF5/Px7OrV8PQELS28l//9V9qVzajSMtDCHHDMjIyyJk7l8amJrbV1LAoK4sUq3Xajr+pvNzfv6HTMWQ0Yhgd5Vd79/rvNBggMXHaarlZSHgIIYLijjVrcIyM0NzcTHNLC1+fNYuE6OhprcHlctHW1ITbbGZhVpZ/P4+SErBYZDvaIJPLVkKIoDCZTDz8ta+x4ZlnIDmZvpGRaTv2d1atAp8Po8lEbmYmc+Li/JtB5eVJcISIhIcQIqgMBgMYjTT39k7bMfWXhuRGRRFtNmM2GiE+HuLjJThCRMJDCBFUFouFWYsX0+h20zEwMC3H7Bwc9Lc0HA7/JlAJCVBQwDMSHCEjfR5CiKArKyujqbKSxq4uMuLiQnac8Y5yAIsFvF5/eKSmkpqWJkNzQyikLY+XXnqJlStXEhUVRdwEf4GOHDnC2rVriYuLIz4+nnvvvZeKiopQliWECKELFy6wZ8cOopqaWJSVFdJjZSck+NevMhrBZgOTCYqLISmJhx9+OKTHvtmFNDxcLhePPvoo3//+9696/9DQEOvWrSMnJ4dDhw6xd+9erFYr69atw+12h7I0IUQIVFZW8sf//E98+/Zx36xZxJjNIT1eW18fmM2Qnu7vHF+4EFJSePYv/iKkxxUhvmz14osvAvDyyy9f9f7a2lr6+/v52c9+RnZ2NgA/+clPWLBgAa2trRQUFISyPCFEkI2MjIDdTr7VSlJMTEiPtam8HFJTobAQdLqQdYz39PTgdDrJzMwMyetHKlU7zOfMmUNSUhKbN2/G5XLhcDjYvHkz8+bNIzc3V83ShBBTsHTpUgrWr+eYw+FvFYTIpvJySEuDoiK+vXHjVYOjp6eHTZs2UVtbO+Xj9PT08MMf/pA//OEPDExT53+kUDU8YmNj2b17N6+99hoWi4WYmBi2bt3Kxx9/jF5/9UaR0+nEbrdf8SOECA9arZbVq1eTtWYNfxoYYHdtLc4gX4LeVF4OVisUFrLxu9+dcJvZ0Yu7He7Zs2fKH/xOp5PFixezevXqCfttb1YBh8dPf/pTNBrNNX+OHj06qddyOBxs3LiRVatWcfDgQfbt28e8efO4//77cTgcV33Oz3/+c2w22/jPpctdQojwYDAYuP/BB7njO9+hOTeXt86cobmnJ3gH0GqhqAi02gm/ZALk5ORw3333sX79+il/8GdmZrJx40aWLFkyxWJnLo2iBLbVVk9PDz3X+YuQl5eH+bKOspdffpkf/vCHX0r/zZs38w//8A90dnaivTjczuVyER8fz+bNm9mwYcOXXtvpdF6xJ7Hdbic7O5vBwUGs07iWjhDi+kZHR/ls925a9u4l3+FgVUEBFqNxyq+3qbwcsrNh1iyZ/Bcidrsdm8123c/UgDvMk5KSSEpKuqHiLhkdHUWr1V4xFvvS774JdgAzmUwTNlOFEOElKiqKdfffT31REfu3buWtM2dYmZpKYUpKwK/l9fn8EwGzsiQ4wkBI+zxaW1upqKigtbUVr9dLRUUFFRUVDA8PA3DPPffQ39/PX/7lX3L69Gmqq6v5zne+g16vZ82aNaEsTQgxjQoLC3l040YyHnyQnWNjbK2pYeSyKwiTsXnvXpg1y79KrlBdwJetAvH000/zyiuvfOn2Xbt2ceeddwKwfft2XnzxRaqqqtBqtZSVlfHSSy+xfPnySR1jsk0sIUR4aG5uZu/WrXiqq1melERxWtp1n7OpvNy/5EhpqbQ6Qmyyn6khDY/pIOEhRORxOp0c2L+ful27yBwcZHVBAbHXmFC4ae9eWLKERcuXs3Tp0mms9OYz2c9UWRhRCDHtTCYTd65Zw/3PPsvgggW81dBAVXs7V/suu6m8HDIzwWyW4AgjEh5CCNVkZWXx6Le+xZzHHmO/TscH1dUMXJyfMU6jgawsnnnmGXWKFFcl4SGEUJXBYGDVqlU89OyzOJcs4e3mZk60tuLz+TjX3++fEGgwyAq5YUaWZBdChIW0tDS+9uSTHCsu5uj27TRWV9Pb3+9f8FCEHQkPIUTY0Ol0LF26lPz8fPZs3w4nToDNxne+8x21SxNfIJethBBhJykpia8+/rh/XkdsrH9rWxFWpOUhhAhLWq2Wb33nO+h0OrVLEVch4SGECFvmEG8mJaZOLlsJIYQImISHEEKIgEl4CCGECJiEhxBCiIBJeAghhAiYhIcQQoiASXgIIYQImISHEEKIgEl4CCGECJiEhxBCiIBJeAghhAiYhIcQQoiASXgIIYQIWMSvqqsoCgB2u13lSoQQIvJd+iy99Nk6kYgPj6GhIQCys7NVrkQIIWaOoaEhbDbbhPdrlOvFS5jz+Xx0dHQQGxuLRqNRuxzsdjvZ2dm0tbVhtVrVLkdVci4+J+fic3IuPheO50JRFIaGhsjIyECrnbhnI+JbHlqtlqysLLXL+BKr1Ro2fxnUJufic3IuPifn4nPhdi6u1eK4RDrMhRBCBEzCQwghRMAkPILMZDLxk5/8BJPJpHYpqpNz8Tk5F5+Tc/G5SD4XEd9hLoQQYvpJy0MIIUTAJDyEEEIETMJDCCFEwCQ8hBBCBEzCIwScTieLFi1Co9FQUVHxpftffvllFixYgNlsJi0tjeeee276i5wm1zsXAL29vWRlZaHRaBgYGJjW+qbTROfi5MmTPPHEE2RnZ2OxWCgpKeH//t//q16h0+Bafy9aW1t58MEHiY6OJikpib/6q7/C5XKpU2gIPfTQQ+Tk5GA2m0lPT+eb3/wmHR0dVzzmyJEjrF27lri4OOLj47n33nsn/Hc03SQ8QuC///f/TkZGxlXv+5d/+Rd+/OMf8z//5/+kurqaTz/9lHXr1k1zhdPnWufiku9+97ssWLBgmipSz0Tn4tixYyQnJ/Paa69RXV3Nj3/8Y55//nn+/d//XYUqp8dE58Lr9bJ+/XpGRkbYu3cvb7zxBu+88w5/+7d/q0KVobVmzRrefPNNamtreeedd2hoaODrX//6+P1DQ0OsW7eOnJwcDh06xN69e7Faraxbtw63261i5RcpIqg+/vhjpbi4WKmurlYA5cSJE+P39fX1KRaLRdmxY4d6BU6ja52LS/7zP/9TueOOO5RPP/1UAZT+/v5pr3M6TOZcXO4HP/iBsmbNmukpbppd61x8/PHHilarVdrb28dv+/3vf6+YTCZlcHBQhWqnz/vvv69oNBrF5XIpiqIoR44cUQCltbV1/DGnTp1SAKW+vl6tMsdJyyOIzp8/zzPPPMNvf/tboqKivnT/9u3b8fl8tLe3U1JSQlZWFo899hhtbW0qVBta1zsXADU1NfzsZz/j1VdfveYCbJFuMufiiwYHB0lISAhxZdPveufiwIEDlJaWXtEqWbduHU6nk2PHjk1nqdOqr6+P119/nZUrV2IwGACYM2cOSUlJbN68GZfLhcPhYPPmzcybN4/c3FyVK5bLVkGjKApPP/003/ve91iyZMlVH9PY2IjP5+Mf//Ef+eUvf8nbb79NX18f99xzz4y6pjuZc+F0OnniiSf4xS9+QU5OzjRXOH0mcy6+6MCBA7z55pv8xV/8RYirm16TORddXV2kpqZecVt8fDxGo5Gurq7pKHNa/Y//8T+Ijo4mMTGR1tZW3n///fH7YmNj2b17N6+99hoWi4WYmBi2bt3Kxx9/jF6v/pq2Eh7X8dOf/hSNRnPNn6NHj/Jv//Zv2O12nn/++Qlfy+fz4Xa7+dd//VfWrVvH8uXL+f3vf8/Zs2fZtWvXNL6rqQnmuXj++ecpKSnhqaeemsZ3EDzBPBeXq66u5uGHH+aFF17gnnvuCfG7CI5gn4urba2gKEpYbLlwPZM9F5f8/d//PSdOnGDbtm3odDq+9a1vjW/C5HA42LhxI6tWreLgwYPs27ePefPmcf/99+NwONR6i+NkeZLr6Onpoaen55qPycvLY8OGDXz44YdX/AX3er3odDqefPJJXnnlFX7zm9+wceNG2trarlhGPjU1lf/9v/83zzzzTMjeRzAE81wsWrSIysrK8ccoioLP50On0/HjH/+YF198MaTv5UYF81xcUlNTw5o1a/jzP/9zXnrppZDVHmzBPBcvvPAC77//PidPnhx/TH9/PwkJCezcuZM1a9aE7H0Ew2TPhdls/tLt586dIzs7m/3797NixQo2b97MP/zDP9DZ2Tl+WdflchEfH8/mzZvZsGFDSN7DpKnW2zLDtLS0KJWVleM/W7duVQDl7bffVtra2hRFUZTa2loFuKLDvLe3V9FqtcrWrVvVKj3oJnMu6uvrr3jMr3/9awVQ9u/fr5w/f17ldxA8kzkXiqIoVVVVSkpKivL3f//3KlYbWpM5F5c6zDs6Osaf98Ybb9wUHeatra0KoOzatUtRFEX513/9VyUtLU3x+Xzjj3G73Up0dLTy+uuvq1Tl5yQ8QqSpqemqo2oefvhhZd68ecq+ffuUyspK5YEHHlDmzp07PsJiJproXFxu165dM3q01SVXOxdVVVVKcnKy8uSTTyqdnZ3jPxcuXFCv0GlwtXPh8XiU0tJSZe3atcrx48eVHTt2KFlZWcpzzz2nXqEhcOjQIeXf/u3flBMnTijNzc3Kzp07ldtuu00pKChQxsbGFEVRlNOnTysmk0n5/ve/r9TU1ChVVVXKU089pdhstivCVS3S5zHNXn31VZYtW8b69eu54447MBgMfPLJJ+MjLMTN56233qK7u5vXX3+d9PT08Z9bb71V7dKmnU6nY8uWLZjNZlatWsVjjz3GI488wj/90z+pXVpQWSwW3n33XdauXcucOXPYuHEjpaWl7NmzZ3x59uLiYj788ENOnTrFihUruP322+no6OCTTz4hPT1d5XcgfR5CCCGmQFoeQgghAibhIYQQImASHkIIIQIm4SGEECJgEh5CCCECJuEhhBAiYBIeQgghAibhIYQQImASHkIIIQIm4SGEECJgEh5CCCECJuEhhBAiYP8/N3fhWourAiMAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f, a = plt.subplots()\n", "geo_s.plot(ax=a,edgecolor='k',facecolor='r',alpha=0.4)\n", "geo_cp.plot(ax=a,edgecolor='w',facecolor='g',alpha=1.0);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Selecionando poços (onshore ou offshore) que se tornaram públicos em 2022 e em 2023" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [], "source": [ "# para onshore, 'T' \n", "p2022_m = p2022_s[p2022_s['TERRA_MAR'] == 'T']\n", "p2022_tcm = p2022_m['CAMPO'].unique()\n", "\n", "p2023_m = p2023_s[p2023_s['TERRA_MAR'] == 'T']\n", "p2023_tcm = p2023_m['CAMPO'].unique()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Selecionando campos relacionados aos poços" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [], "source": [ "geo_cp_m_2022 = geo_cp.apply(lambda row: row[geo_cp['NOM_CAMPO'].isin(p2022_tcm)])\n", "geo_cp_m_2023 = geo_cp.apply(lambda row: row[geo_cp['NOM_CAMPO'].isin(p2023_tcm)])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Plota campos de interesse" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAN8AAAGdCAYAAACSHqb/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAZMklEQVR4nO3df0xd9f3H8dcttrewwm3ZxVW/vS21M211/lFN1OK2FBz9oXEkcxjQbSxuSDXoXGOWYQm07ovdHyZro1vqMoIuZX9QVzczF51b0cSUWKle1Nsfrl0rFGRfavVeKi0g/Xz/cNysQoFb7uF9W56P5CZez7kfPvfmPnt+3HvA55xzAjDlZlhPAJiuiA8wQnyAEeIDjBAfYIT4ACPEBxghPsDIZdYTsHD27Fl1dXUpMzNTPp/Pejq4hDjn1NvbqyuvvFIzZoy9bZuW8XV1dSkUCllPA5ewjo4OLViwYMx1pmV8mZmZkj5/gbKysoxng0tJLBZTKBSKv8fGMi3jG97VzMrKIj54YiKHM5xwAYwQH2CE+AAjxAcYIT7ACPEBRogPMEJ8gBHiA4wQH2CE+AAjxAcYIT7AyLS8qgG2Ntc9rt7TZ0x+dmb6bNVufNTkZ3+Rp/HV1dXpxRdfVDgc1qxZs/TJJ5+cs7ytrU2//OUv9frrr+vEiRPKzc3V+vXr9ZOf/GTMcX/729/qD3/4g9566y319vbq448/1ty5c717Ikiq3tNn9J31D5v87F3bt5r83NF4uts5MDCg4uJi3X///aMu37dvn3JycrRjxw5FIhFt3LhRVVVVeuqpp8Yct6+vT2vXrtWjj6bGv2DAhfB0y7d582ZJ0jPPPDPq8nvvvfec+1dddZVaWlq0a9cuVVZWnnfchx9+WJL06quvJmOagImUO+aLRqPKzs5O6pj9/f3q7++P34/FYkkdH7gQKXW2s6WlRU1NTaqoqEjquFu2bFEgEIjf+OVJSAUJx7dp0yb5fL4xb62trQlPJBKJqKioSDU1NSosLEz48WOpqqpSNBqN3zo6OpI6Pi4eQ740PVJdo0eqa7S57nHTuSS821lZWamSkpIx18nNzU1ozP3796ugoEDl5eWqrq5OdErj8vv98vv9SR8XF5/iigfj/2195jPh+ILBoILBYNImEIlEVFBQoLKyMtXV1SVtXCDVeXrM197ernA4rPb2dg0NDSkcDiscDuvUqVOSPg8vPz9fhYWF2rBhg7q7u9Xd3a2enp74GJ2dnVq2bJn27t0b/3/d3d0Kh8M6fPiwJOndd99VOBzWyZMnvXw6QFJ5erazpqZGzz77bPz+ihUrJEnNzc1atWqVdu7cqZ6eHjU2NqqxsTG+3qJFi3Ts2DFJ0uDgoA4dOqS+vr748u3bt8c/xpCkb37zm5KkhoYG/fCHP/TwGeFSMnz8N2yqv/3ic865KftpKSIWiykQCCgajfJLcw08Ul1j9g2XsezavlVP/O9jkxojkfdWSn3UAEwnxAcYSblvuABWpvoYkPiA//jvzwAl7z8HZLcTMEJ8gBHiA4wQH2CEEy7AeZz57KynZz+JDziPuys3nHM/2Wc/2e0EjBAfYIT4ACPEBxghPsAI8QFG+KgBUy4zffaop+2HfGkjvtycSpJ91QPxYcqd7w3732/sVJTsqx7Y7QSMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIARPmQHzmPn008qzQ2dd3lm+uxJjU98wHmkuaFJ/+2GsbDbCRghPsAI8QFGiA8wQnyAEc52Av/xxY8WJvtRwniID/gPrz9a+CJ2OwEjxAcYYbcT09ZUH+N9EfFh2prqY7wvIj6Y21z3uHpPn9GQL816KlOK+GCu9/QZfWf9w9bTmHKccAGMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUa4pAjmMtNna9f2rSY/15LPOedMZ2AgFospEAgoGo0qKyvLejq4hCTy3mK3EzBCfIAR4gOMEB9ghPgAI8QHGCE+wIin8dXV1SkvL08ZGRmaO3fuiOVtbW0qLS1VKBRSenq6li9frm3bto055smTJ/Xggw9q6dKlysjI0MKFC/XQQw8pGo169CwAb3j6DZeBgQEVFxdr5cqVqq+vH7F83759ysnJ0Y4dOxQKhbRnzx7dd999SktLU2Vl5ahjdnV1qaurS0888YSuueYaffDBB1q/fr26urr03HPPefl0gORyU6ChocEFAoEJrfvAAw+4/Pz8hMZvampys2bNcoODgxNaPxqNOkkuGo0m9HOA8STy3kq573ZGo1FlZ2cn/JisrCxddtnoT6e/v1/9/f3x+7FYbFJzBJIhpU64tLS0qKmpSRUVFRN+zEcffaRf/OIXYz5my5YtCgQC8VsoFErGdIFJSTi+TZs2yefzjXlrbW1NeCKRSERFRUWqqalRYWHhhB4Ti8V0++2365prrlFtbe1516uqqlI0Go3fOjo6Ep4fkGwJ73ZWVlaqpKRkzHVyc3MTGnP//v0qKChQeXm5qqurJ/SY3t5erV27VnPmzNHzzz+vmTNnnnddv98vv9+f0JwAryUcXzAYVDAYTNoEIpGICgoKVFZWprq6ugk9JhaLac2aNfL7/XrhhRc0e7btdVnAhfD0mK+9vV3hcFjt7e0aGhpSOBxWOBzWqVOnJH0eXn5+vgoLC7VhwwZ1d3eru7tbPT098TE6Ozu1bNky7d27V9LnW7zVq1fr008/VX19vWKxWPxxQ0NDo84DSElennYtKytzkkbcmpubnXPO1dbWjrp80aJF8TGOHj16zmOam5tHfYwkd/To0QnNi48a4JVE3ltcyc6V7EgirmQHLgLEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYMTT+Orq6pSXl6eMjAzNnTt3xPK2tjaVlpYqFAopPT1dy5cv17Zt28Ydt6KiQkuWLFF6erpycnJUVFSkgwcPevAMAO94Gt/AwICKi4t1//33j7p83759ysnJ0Y4dOxSJRLRx40ZVVVXpqaeeGnPcG264QQ0NDTpw4IBefvllOee0evVqDQ0NefE0AG+4KdDQ0OACgcCE1n3ggQdcfn5+QuO3tbU5Se7w4cMTWj8ajTpJLhqNJvRzgPEk8t66zLj9EaLRqLKzsye8/qeffqqGhgYtXrxYoVBo1HX6+/vV398fvx+LxSY9T2CyUuqES0tLi5qamlRRUTHuur/5zW80Z84czZkzRy+99JJeeeUVzZo1a9R1t2zZokAgEL+dL1JgKiUc36ZNm+Tz+ca8tba2JjyRSCSioqIi1dTUqLCwcNz177nnHr399tt67bXXdPXVV+uuu+7SmTNnRl23qqpK0Wg0fuvo6Eh4fkCyJbzbWVlZqZKSkjHXyc3NTWjM/fv3q6CgQOXl5aqurp7QY4a3YldffbVuvvlmzZs3T88//7xKS0tHrOv3++X3+xOaE+C1hOMLBoMKBoNJm0AkElFBQYHKyspUV1d3weM45845rgNSnafHfO3t7QqHw2pvb9fQ0JDC4bDC4bBOnTol6fPw8vPzVVhYqA0bNqi7u1vd3d3q6emJj9HZ2ally5Zp7969kqR//etf2rJli/bt26f29na1tLTorrvuUnp6um677TYvnw6QXF6edi0rK3OSRtyam5udc87V1taOunzRokXxMY4ePXrOYzo7O926devc5Zdf7mbOnOkWLFjg7r77bnfw4MEJz4uPGuCVRN5bPuecM6neUCwWUyAQUDQaVVZWlvV0cAlJ5L2VUh81ANMJ8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxhJub/PN12szFurjz4emNC6X543Sy17XvJ4RphqxGfko48HdPmC8f8OoST93/GnPZ4NLLDbCRghPsAI8QFGiA8wQnyAEeIDjBAfYIT4ACPEBxghPsAI8QFGiA8wwherp9jw1Qx9fbOspwJjxDfFErmaAZc2djsBI8QHGGG3M0Ude79BGRmfX+meNcd4MvAE8aWojIwBvX9gt/U04CHimyITPcs5vMX78jzOhl7qiG+KTPQsJ1u86YMTLoAR4gOMsNuZIjjWm36Ib4p8ed6sMX//Zuh/ZqllD8d60wnxTRF+4zS+iGM+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMeBpfXV2d8vLylJGRoblz545Y3tbWptLSUoVCIaWnp2v58uXatm3bhMd3zmndunXy+Xz605/+lLyJA1PA00uKBgYGVFxcrJUrV6q+vn7E8n379iknJ0c7duxQKBTSnj17dN999yktLU2VlZXjjr9161b5fD4vpg54z02BhoYGFwgEJrTuAw884PLz88ddLxwOuwULFrgPP/zQSXLPP//8hOcTjUadJBeNRif8GGAiEnlvpdzFtNFoVNnZ2WOu09fXp9LSUj311FOaP3/+uGP29/erv78/fj8Wi016nsBkpdQJl5aWFjU1NamiYuxfsffTn/5UeXl5KioqmtC4W7ZsUSAQiN9CoVAypgtMSsLxbdq0ST6fb8xba2trwhOJRCIqKipSTU2NCgsLz7veCy+8oN27d2vr1q0THruqqkrRaDR+6+joSHh+QLIlvNtZWVmpkpKSMdfJzc1NaMz9+/eroKBA5eXlqq6uHnPd3bt368iRIyPOnt555536xje+oVdffXXEY/x+v/x+f0JzAryWcHzBYFDBYDBpE4hEIiooKFBZWZnq6urGXf/nP/+5fvzjH5/z/6677jr96le/0h133JG0eQFe8/SES3t7u06ePKn29nYNDQ0pHA5Lkr761a9qzpw5ikQiys/P1+rVq7VhwwZ1d3dLktLS0pSTkyNJ6uzs1K233qrf//73uvHGGzV//vxRT7IsXLhQixcv9vLpAEnlaXw1NTV69tln4/dXrFghSWpubtaqVau0c+dO9fT0qLGxUY2NjfH1Fi1apGPHjkmSBgcHdejQIfX19Xk5VWDK+ZxzznoSUy0WiykQCCgajSorK8t6OriEJPLeSqmPGoDphPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUYus57ApWBz3ePqPX1mQutmps9W7cZHPZ4RLgbElwS9p8/oO+sfntC6u7Zv9XQuuHiw2wkYIT7ACPEBRjjmm4ThEy1DvjTrqeAiRHyTkMiJFuCL2O0EjBAfYIT4ACPEBxjxNL66ujrl5eUpIyNDc+fOHbG8ra1NpaWlCoVCSk9P1/Lly7Vt27Zxx121apV8Pt85t5KSEg+eQfIN+dL0SHVN/La57nHrKcGIp2c7BwYGVFxcrJUrV6q+vn7E8n379iknJ0c7duxQKBTSnj17dN999yktLU2VlZVjjl1eXq7HHnssfj89PT3p8/dCccWD59zn62bTl6fxbd68WZL0zDPPjLr83nvvPef+VVddpZaWFu3atWvc+DIyMjR//vykzNPS8JZwGF+8nj5S7nO+aDSq7OzscddrbGzUjh079JWvfEXr1q1TbW2tMjMzR123v79f/f398fuxWCxp850stoTTV0rF19LSoqamJr344otjrnfPPfdo8eLFmj9/vt577z1VVVWpra1Nr7zyyqjrb9myJb4VBlJFwidcNm3aNOJkxxdvra2tCU8kEomoqKhINTU1KiwsHHPd8vJyfetb39LXvvY1lZSU6LnnntPf//53vfXWW6OuX1VVpWg0Gr91dHQkPD8g2RLe8lVWVo57ZjE3NzehMffv36+CggKVl5eruro60Snp+uuv18yZM/XPf/5T119//Yjlfr9ffr8/4XHHk5k+W7u2b9WQL23E7uOFOvPZWT1SXcOx3zSQcHzBYFDBYDBpE4hEIiooKFBZWZnq6uoueIzBwUFdccUVSZvXRAzH8d8nTCbr7soNkjj2mw48/Zyvvb1d4XBY7e3tGhoaUjgcVjgc1qlTpyR9Hk1+fr4KCwu1YcMGdXd3q7u7Wz09PfExOjs7tWzZMu3du1eSdOTIET322GNqbW3VsWPH9Ne//lXFxcVasWKFbrnlFi+fDpBUnp5wqamp0bPPPhu/v2LFCklSc3OzVq1apZ07d6qnp0eNjY1qbGyMr7do0SIdO3ZMkjQ4OKhDhw6pr69PkjRr1iz94x//0LZt23Tq1CmFQiHdfvvtqq2tVVoal/bg4uFzzjnrSUy1WCymQCCgaDSqrKysEcsT+YVIkpJ6zDds1/ateuJ/Hxt/RaSU8d5b/y2lPmpIFVynh6nAF6tHMXwWc+fTT1pPBZcwtnyj8OIsJvBFxJdidj79pNLckDLTZ1tPBR4jvhST5oY40TJNcMw3Bo794CW2fGPg2A9eYssHGCE+wAjxTQDHfvACx3wTwLEfvMCWDzBCfIAR4ksAx35IJo75EsCxH5KJLd8FYAuIZGDLdwHYAiIZ2PIBRogPMMJu5yQk4/d2Dl+/999jYnogvklIxrEf1+9NX+x2JgFnP3Eh2PIlAWc/cSGIL4ky02dr59NPjnn8xzEehvFLc8f5xaaJGu8X7vIHUC5tiby3iC/J8WF6S+S9xQkXwAjxAUaIDzBCfIAR4gOMEB9ghPgAI8QHGCE+wAjxAUaIDzBCfIAR4gOMEB9gZFpeTDt8FVUsFjOeCS41w++piVypNy3j6+3tlSSFQiHjmeBS1dvbq0AgMOY60/Ji2rNnz6qrq0uZmZny+XxJGzcWiykUCqmjo4OLdMdwKb9Ozjn19vbqyiuv1IwZYx/VTcst34wZM7RgwQLPxs/Kyrrk3lReuFRfp/G2eMM44QIYIT7ACPElkd/vV21trfx+v/VUUhqv0+em5QkXIBWw5QOMEB9ghPgAI8QHGCG+BH3729/WwoULNXv2bF1xxRX6/ve/r66urnPWefPNN3Xrrbdq7ty5mjdvnlavXq1wODzmuBUVFVqyZInS09OVk5OjoqIiHTx40MNn4h0vXqOTJ0/qwQcf1NKlS5WRkaGFCxfqoYceUjQa9fjZeIf4EpSfn6+mpiYdOnRIf/zjH3XkyBF997vfjS/v7e3VmjVrtHDhQr3xxht6/fXXlZWVpTVr1mhwcPC8495www1qaGjQgQMH9PLLL8s5p9WrV2toaOi8j0lVXrxGXV1d6urq0hNPPKF3331XzzzzjF566SX96Ec/mqqnlXwOk/LnP//Z+Xw+NzAw4Jxz7s0333SSXHt7e3ydd955x0lyhw8fnvC4bW1tCT8mVXn1GjU1NblZs2a5wcHBpM95KrDlm4STJ0+qsbFReXl5mjlzpiRp6dKlCgaDqq+v18DAgE6fPq36+npde+21WrRo0YTG/fTTT9XQ0KDFixdf9FdeePUaSYr/JaDLLrtIv6JsXf/F6Gc/+5nLyMhwktzNN9/sTpw4cc7y9957zy1ZssTNmDHDzZgxwy1btsx98MEH447761//2n3pS19yktyyZcsu6q2eV6/RsBMnTriFCxe6jRs3JnvqU4b4nHO1tbVO0pi3N998M75+T0+PO3TokPvb3/7mbrnlFnfbbbe5s2fPOuec6+vrczfeeKP7wQ9+4Pbu3etaWlrcnXfe6a699lrX19c35jw++eQT9/7777vXXnvN3XHHHe766693p0+f9vS5T1SqvEbOOReNRt1NN93k1q5dG9+VvRjx9TJJJ06c0IkTJ8ZcJzc3V7Nnj/wTzsePH1coFNKePXu0cuVK1dfX69FHH9WHH34Yv55rYGBA8+bNU319vUpKSiY0p+HH/O53v1NpaWniTyrJUuU1Gj5Zk5GRob/85S+j/ryLxUW6s5xcwWBQwWDwgh47/G9Xf3+/JKmvr08zZsw45yLd4ftnz55NeOzhca2lwmsUi8W0Zs0a+f1+vfDCCxd1eJI45kvEG2+84Z588kn39ttvu2PHjrndu3e7r3/9627JkiXuzJkzzjnnDhw44Px+v7v//vvd/v373Xvvvee+973vuUAg4Lq6upxzzh0/ftwtXbrUvfHGG845544cOeIef/xx19ra6j744AO3Z88eV1RU5LKzs92///1vs+d7Ibx6jWKxmLvpppvcdddd5w4fPuw+/PDD+O2zzz4ze76TQXwJeOedd1x+fr7Lzs52fr/f5ebmuvXr17vjx4+fs97wcU4gEHDz5s1zBQUFrqWlJb786NGjTpJrbm52zjnX2dnp1q1b5y6//HI3c+ZMt2DBAnf33Xe7gwcPTuXTSwqvXqPm5ubzHmsePXp0Cp9h8nDMBxjhcz7ACPEBRogPMEJ8gBHiA4wQH2CE+AAjxAcYIT7ACPEBRogPMEJ8gJH/B+hNibWqBDJwAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f, a = plt.subplots()\n", "\n", "#geo_s.plot(ax=a,edgecolor='k',facecolor='gray',alpha=.2)\n", "geo_cp_m_2022.plot(ax=a,edgecolor='k',lw=.5,facecolor='lightblue',alpha=.7)\n", "geo_cp_m_2023.plot(ax=a,edgecolor='k',lw=.5,facecolor='darkblue',alpha=.7);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- Lista de campos participantes da RV" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Atividade de divulgação em 2022 nos campos de: \n", "TIÊ, ARAÇÁS, MATA DE SÃO JOÃO.\n", "\n", "Atividade de divulgação em 2023 nos campos de: \n", "TIÊ.\n" ] } ], "source": [ "cp_list_2022 = geo_cp_m_2022.NOM_CAMPO.tolist()\n", "cp_list_2023 = geo_cp_m_2023.NOM_CAMPO.tolist()\n", "\n", "print('Atividade de divulgação em 2022 nos campos de: ')\n", "print(*[x for x in cp_list_2022], sep=', ', end='.\\n\\n')\n", "\n", "print('Atividade de divulgação em 2023 nos campos de: ')\n", "print(*[x for x in cp_list_2023], sep=', ', end='.\\n')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### RV finalística" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABk4AAAIcCAYAAACwxlKLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACm3klEQVR4nOzdeVRV1f//8ddlBhFwwCkRVJwNsbTSHHHAITWzTLSUnNIy8+PPBs0xw8ysNP1kE0IOOVCapmmKijmWZVcTUUtRHD/ixCAJCOf3h3C/3hjExAB7Pta6a3HO3mfv9znnnrX0vO/e22QYhiEAAAAAAAAAAADIpqgDAAAAAAAAAAAAKC5InAAAAAAAAAAAAGQhcQIAAAAAAAAAAJCFxAkAAAAAAAAAAEAWEicAAAAAAAAAAABZSJwAAAAAAAAAAABkIXECAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAADgrklLS9O0adP0xx9/FHUoQLETExOjkJAQpaenF3UoAAAAAICbkDgBAAB3zSuvvKJDhw7J19fXar/JZNI333yT53HHjx+XyWSS2WyWJEVFRclkMunKlSt3L1iUKOHh4fLw8CjqMP62jIwMBQcHa9u2bXrrrbeKOhwAAAAAwE1InAAAUISCg4NlMplkMplkb2+vGjVqaMyYMbp69WpRh3bHVq9erf379+vzzz+/47aaN2+us2fPyt3dvRAiK76ioqLUo0cPVa5cWaVKlZK/v78WL16co97WrVv14IMPysnJSTVq1NDHH39sVf7ZZ5+pZcuWKlOmjMqUKaP27dvrp59+sqrz9ttvq2nTpipdurQqVKigxx9/XIcPH77jc8hOemV/HBwc5Ovrq7feekuGYdxx+9mefvppHTlypNDay9amTRur2GvWrKmxY8cqNTW1UPt59913FRgYqNWrVysyMtKSJAQAAAAAFD27og4AAIB/u06dOiksLEzp6enatm2bBg8erKtXr2revHlFHdod6d69u7p3714obTk4OKhSpUqF0lZxtnPnTvn5+em1115TxYoVtXbtWvXv319ubm7q1q2bJCk2NlZdunTRkCFDtGjRIu3YsUMvvPCCPD091atXL0k3EjBBQUFq3ry5nJycNGPGDHXs2FHR0dG67777JN1Ivrz44otq2rSprl+/rjfeeEMdO3bUwYMHVapUqTs+l8jISDVo0ECpqanavn27Bg8erMqVK2vQoEF33LYkOTs7y9nZuVDa+qshQ4bozTffVFpamvbs2aPnnntO0o1kU2F5/fXXLX/v2LGj0NoFAAAAANw5RpwAAFDEHB0dValSJXl5ealv377q16+fZRqr1NRUjRw5UhUqVJCTk5NatGihPXv2WB0fHR2trl27ys3NTaVLl1bLli119OhRSVJmZqbefPNNVa1aVY6OjvL399f69estx6alpWnEiBGqXLmynJyc5OPjk+/L4eDgYD3++OOaMmWKKlSoIDc3Nz3//PNKS0uz1PHx8dGsWbOsjvP399fkyZOt9p09e1adO3eWs7OzqlevroiIiDz7zW2qrh07dqh169ZycXFRmTJlFBgYqMuXLxfoul2+fFn9+vWTp6ennJ2dVatWLYWFheXZv2EYmjFjhmrUqCFnZ2c1atRIX331VY74vv/+ezVu3FjOzs4KCAjQ+fPntW7dOtWrV09ubm4KCgpSSkpKnv2MGzdOU6dOVfPmzVWzZk2NHDlSnTp10sqVKy11Pv74Y1WrVk2zZs1SvXr1NHjwYA0cOFAzZ8601Fm8eLFeeOEF+fv7q27duvrss8+UmZmpTZs2WeqsX79ewcHBatCggRo1aqSwsDDFxcXpl19+yTO+21GuXDlVqlRJ3t7e6tevn5o3b669e/dayvfs2aMOHTqofPnycnd3V+vWra3KJenKlSsaOnSoKlasKCcnJzVs2FBr1qyRlHOqrqNHj6pHjx6qWLGiXF1d1bRpU0VGRv6t2F1cXFSpUiVVq1ZNvXr1UocOHbRhwwZL+a2+D9KdPZfZo3ZWrFihtm3bysXFRY0aNdKuXbus+sjvGVi/fr1atGghDw8PlStXTo899pilf0lq1qyZVfJGkuLj42Vvb68tW7ZIuvGc9O/fX2XKlJGLi4s6d+6s33///W9dUwAAAAAoSUicAABQzDg7O1sWi3711Vf19ddf64svvtDevXvl6+urwMBAXbp0SZJ0+vRptWrVSk5OTtq8ebN++eUXDRw4UNevX5ckzZ49W++9955mzpyp/fv3KzAwUN27d7e8/Pzwww+1evVqLV++XIcPH9aiRYvk4+OTb3ybNm1STEyMtmzZoiVLlmjlypWaMmXKbZ/nhAkT1KtXL+3bt0/PPPOMgoKCFBMTU6BjzWaz2rVrpwYNGmjXrl3avn27unXrpoyMjAJdtwkTJujgwYNat26dYmJiNG/ePJUvXz7P/saPH6+wsDDNmzdP0dHR+s9//qNnnnlGW7dutao3efJkzZ07Vzt37tTJkyfVu3dvzZo1S19++aXWrl2rjRs3as6cObd1nRISElS2bFnL9q5du9SxY0erOoGBgfr555/zXGQ8JSVF6enpVu3k1o+kfOv8XT///LP27t2rhx9+2LIvKSlJAwYM0LZt27R7927VqlVLXbp0UVJSkqQbyYXOnTtr586dWrRokQ4ePKjp06fL1tY21z6Sk5PVpUsXRUZG6tdff1VgYKC6deumuLg4S53Jkyff8vv9V/v27dOOHTtkb29v2Xer78OdPpfZ3njjDY0ZM0Zms1m1a9dWUFCQpY1bPQNXr17V6NGjtWfPHm3atEk2Njbq2bOnMjMzJUn9+vXTkiVLrKZPW7ZsmSpWrKjWrVtLupEo/fnnn7V69Wrt2rVLhmGoS5cuLGYPAAAA4N5nAACAIjNgwACjR48elu0ff/zRKFeunNG7d28jOTnZsLe3NxYvXmwpT0tLM6pUqWLMmDHDMAzDGDt2rFG9enUjLS0t1/arVKlihISEWO1r2rSp8cILLxiGYRgvvfSSERAQYGRmZhY43rJlyxpXr1617Js3b57h6upqZGRkGIZhGN7e3sYHH3xgdVyjRo2MSZMmWbYlGcOGDbOq8/DDDxvDhw83DMMwYmNjDUnGr7/+ahiGYWzZssWQZFy+fNkwDMMICgoyHn300VxjLMh169atm/Hcc88V6JyTk5MNJycnY+fOnVb7Bw0aZAQFBVnFFxkZaSl/++23DUnG0aNHLfuef/55IzAwsED9GoZhREREGA4ODsaBAwcs+2rVqpXjnu7YscOQZJw5cybXdl544QWjZs2axp9//plreWZmptGtWzejRYsWBY4tL9n3ztnZ2ShVqpRhb29vSDKGDh2a73HXr183SpcubXz77beGYRjG999/b9jY2BiHDx/OtX5YWJjh7u6eb5v169c35syZY9meM2eOERAQkO8xrVu3Nuzt7Y1SpUoZDg4OhiTDxsbG+OqrrwzDKNj34U6fy+xr+Pnnn1vKo6OjDUlGTEyMYRj5PwO5OX/+vCHJ+O233yzbdnZ2xg8//GCp06xZM+OVV14xDMMwjhw5YkgyduzYYSm/cOGC4ezsbCxfvrzA/QIAAABAScSIEwAAitiaNWvk6uoqJycnNWvWTK1atdKcOXN09OhRpaen69FHH7XUtbe310MPPWQZmWE2m9WyZUurX8NnS0xM1JkzZ6yOl6RHH33UcnxwcLDMZrPq1KmjkSNHWk1HlJdGjRrJxcXFst2sWTMlJyfr5MmTt3XezZo1y7F9uyNOclOQ6zZ8+HAtXbpU/v7+evXVV7Vz5848+zp48KCuXbumDh06yNXV1fJZsGCB1dRHkuTn52f5u2LFinJxcVGNGjWs9p0/f75A5xgVFaXg4GB99tlnatCggVWZyWSy2jayRg38db8kzZgxQ0uWLNGKFSvk5OSUa18jRozQ/v37tWTJkjzjiYuLszr/adOm5Rv/smXLZDabtW/fPi1btkyrVq2ymhrq/PnzGjZsmGrXri13d3e5u7srOTnZMkLEbDaratWqql27dr79ZLt69apeffVV1a9fXx4eHnJ1ddWhQ4esRpyMGDHCarqyvPTr109ms1m7du1S7969NXDgQMv6MQX5Ptzpc5nt5u9T5cqVJcny/cnvGZBuPAd9+/ZVjRo15ObmpurVq0uS5Xp4enqqQ4cOWrx4saQba+fs2rVL/fr1kyTFxMTIzs7OapRQuXLlVKdOnQI/pwAAAABQUrE4PAAARaxt27aaN2+e7O3tVaVKFcvL1rNnz0rK/SV59r6CLI6d3/EPPPCAYmNjtW7dOkVGRqp3795q3759jvUaCiK7TRsbG6vpfyQVeGqf3F785ya/884riXDzeXfu3FknTpzQ2rVrFRkZqXbt2unFF1+0WickW/bURmvXrrUsrJ7N0dHRavvmF+UmkynHi3OTyWRpLz9bt25Vt27d9P7776t///5WZZUqVdK5c+es9p0/f152dnYqV66c1f6ZM2dq2rRpioyMtHoJf7OXXnpJq1ev1g8//KCqVavmGVOVKlVkNpst27ea0svLy0u+vr6SpHr16unYsWOaMGGCJk+eLCcnJwUHBys+Pl6zZs2St7e3HB0d1axZM8t6Obe78Psrr7yi77//XjNnzpSvr6+cnZ315JNPWq2/U1Du7u6W2BctWqQGDRooNDRUgwYNKtD34U6fy2x//T5J//d9vFUf3bp1k5eXlz777DNVqVJFmZmZatiwodX16Nevn15++WXNmTNHX375pWW9m+x4cpNbnAAAAABwr2HECQAARaxUqVLy9fWVt7e31YtSX19fOTg4aPv27ZZ96enp+vnnn1WvXj1JN36Rvm3btlwTE25ubqpSpYrV8ZK0c+dOy/HZ9Z5++ml99tlnWrZsmb7++mvLWiC52bdvn/7880/L9u7du+Xq6mp56e7p6WlJ+kg3fmEfGxubo53du3fn2K5bt26e/d7Mz88vz5EDBblu2XEGBwdr0aJFmjVrlj799NNc26tfv74cHR0VFxcnX19fq4+Xl1eB4r0dUVFR6tq1q6ZPn66hQ4fmKG/WrJk2btxotW/Dhg1q0qSJ1ffn3Xff1dSpU7V+/Xo1adIkRzuGYWjEiBFasWKFNm/ebBmRkBc7Ozurc7/dtVBsbW11/fp1y4v7bdu2aeTIkerSpYsaNGggR0dHXbhwwVLfz89Pp06d0pEjRwrU/rZt2xQcHKyePXvq/vvvV6VKlXT8+PHbijE39vb2GjdunMaPH6+UlJQCfR8K47m8lfyegYsXLyomJkbjx49Xu3btVK9ePcui8Td7/PHHde3aNa1fv15ffvmlnnnmGUtZ/fr1df36df34449W7R45cuS24gQAAACAkojECQAAxVSpUqU0fPhwvfLKK1q/fr0OHjyoIUOGKCUlRYMGDZJ0Y+qhxMRE9enTRz///LN+//13LVy4UIcPH5Z041f477zzjpYtW6bDhw/r9ddfl9ls1ssvvyxJ+uCDD7R06VIdOnRIR44cUUREhCpVqiQPD48840pLS9OgQYMsi6tPmjRJI0aMkI3NjX9WBAQEaOHChdq2bZsOHDigAQMG5Lqgd0REhObPn68jR45o0qRJ+umnnzRixIgCXZuxY8dqz549euGFF7R//34dOnRI8+bN04ULFwp03SZOnKhVq1bpjz/+UHR0tNasWZPny+DSpUtrzJgx+s9//qMvvvhCR48e1a+//qr//ve/+uKLLwoUb0FlJ01GjhypXr166dy5czp37pxVImvYsGE6ceKERo8erZiYGM2fP1+hoaEaM2aMpc6MGTM0fvx4zZ8/Xz4+PpZ2kpOTLXVefPFFLVq0SF9++aVKly5tqXNzUuxOXLx4UefOndOpU6e0bt06zZ49W23btpWbm5ukGwmuhQsXKiYmRj/++KP69etnNYqidevWatWqlXr16qWNGzdaRkatX78+1/58fX21YsUKy/Rgffv2zTG6Z+7cuflOb5WXvn37ymQy6aOPPirQ9+FOn8uCyO8ZKFOmjMqVK6dPP/1Uf/zxhzZv3qzRo0fnaKNUqVLq0aOHJkyYoJiYGPXt29dSVqtWLfXo0UNDhgzR9u3btW/fPj3zzDO677771KNHj9u+hgAAAABQohTV4ioAACDn4vB/9eeffxovvfSSUb58ecPR0dF49NFHjZ9++smqzr59+4yOHTsakgxJRsuWLS0LkmdkZBhTpkwx7rvvPsPe3t5o1KiRsW7dOsuxn376qeHv72+UKlXKcHNzM9q1a2fs3bv3lvFOnDjRKFeunOHq6moMHjzYuHbtmqVOQkKC0bt3b8PNzc3w8vIywsPDc10c/r///a/RoUMHw9HR0fD29jaWLFliKb/V4vCGYRhRUVFG8+bNDUdHR8PDw8MIDAy0lN/quk2dOtWoV6+e4ezsbJQtW9bo0aOHcezYsTzPOzMz05g9e7ZRp04dw97e3vD09DQCAwONrVu35hlfbouXT5o0yWjUqFG+1zf7Pt78ad26tVW9qKgoo3HjxoaDg4Ph4+NjzJs3z6rc29s713b+eg9y+4SFheUZX0Fk37vsj62trVG1alVjyJAhxvnz5y319u7dazRp0sRwdHQ0atWqZURERBje3t7GBx98YKlz8eJF47nnnjPKlStnODk5GQ0bNjTWrFljGEbO6xsbG2u0bdvWcHZ2Nry8vIy5c+carVu3Nl5++WVLnUmTJhne3t75xv/XY7KFhIQYnp6eRlJS0i2/D4bxf8+li4uLUbp06dt6Lv/6/TcMw7h8+bIhydiyZYtlX37PwMaNG4169eoZjo6Ohp+fnxEVFWVIMlauXGl1XmvXrjUkGa1atcpxzpcuXTKeffZZw93d3XB2djYCAwONI0eO5Hv9AAAAAOBeYDKMPCYwBgAAJUpycrLatm2rnTt35roodWEIDg7WlStX9M0339yV9gEAAAAAAIoaU3UBAHAPOH78uM6dO6f4+HgdOHCgqMMBAAAAAAAosUicAABwDwgPD1f9+vVVsWJF1a5du6jDAQAAAAAAKLGYqgsAAAAAAAAAACALI04AAAAAAAAAAACykDgBAAAAAAAAAADIQuIEAAAAAAAAAAAgi11RBwAAAFASZGZm6syZMypdurRMJlNRh3PPMAxDSUlJqlKlimxs+E0PAAAAAKDokTgBAAAogDNnzsjLy6uow7hnnTx5UlWrVi3qMAAAAAAAIHECAABQEKVLl5Z04wW/m5tbEUdz70hMTJSXl5fl+gIAAAAAUNRInAAAABRA9vRcbm5uJE7uAqY/AwAAAAAUF0wkDQAAAAAAAAAAkIXECQAAAAAAAAAAQBam6gIAAABuISMjQ+np6UUdBoA7ZG9vL1tb26IOAwAAAMUciRMAAAAgH8nJyTp16pQMwyjqUADcIZPJpKpVq8rV1bWoQwEAAEAxRuIEAAAAyENGRoZOnTolFxcXeXp6sog9UIIZhqH4+HidOnVKtWrVYuQJAAAA8kTiBAAAFBshISFau3atzGazHBwcdOXKFavyffv2afr06dq+fbsuXLggHx8fDRs2TC+//HK+7X766af68ssvtXfvXiUlJeny5cvy8PC4eyeCe0Z6eroMw5Cnp6ecnZ2LOhwAd8jT01PHjx9Xeno6iRMAAADkicXhAQBAsZGWlqannnpKw4cPz7X8l19+kaenpxYtWqTo6Gi98cYbGjt2rObOnZtvuykpKerUqZPGjRt3N8LGvwAjTYB7A88yAAAACoIRJwAAoNiYMmWKJCk8PDzX8oEDB1pt16hRQ7t27dKKFSs0YsSIPNsdNWqUJCkqKqowwrznLFyyVCmpaXfUhoujg54N6lNIEQEAAAAAUHRInAAAgBItISFBZcuWLfR2U1NTlZqaatlOTEws9D6Ki5TUNLV8rOcdtbFtzcpCiubekJycrCmTQ/Tdtxt0Pc2QnYPUpVugJk1+o1AWpb5y5YratGlj2U5JSdGxY8d0/vx5y/PQpUsXnTlzRtKNae7uv/9+2djYqHTp0tq2bZtMJpOSkpLk6uqqNm3aKC4uTm5ubpY2Z86cqfbt2+foe+DAgdqxY4ecnZ3l5uamDz/8UP7+/pKk4OBgRUZGytPTU1evXlXFihU1dOhQPfvss1ZtTJw4USEhITp27Ji8vb1zPcfw8HCNGjVKPj4+km6sT/Hmm2+qR48eljqGYah69ery9fVVZGTkbV/HwpacnKwPQkK0a8MGORqZumYyqXnHQP3nDe77zSZPnqyPPvpIVapUkSTZ29tr9uzZat68uaVOUlKSKleurD59+ujzzz+/o+sGAAAA3C4SJwAAoMTatWuXli9frrVr1xZ622+//bZlBAxwO5KTk9W2VaAuHS0tZ5uGcjCZZBiGVn7+q6I2B2rLD9/f8Ut0Dw8Pmc1my/bMmTO1detWqyTid999Z/nbZDJp586d+fb74Ycf6rHHHrtl348//rg+/fRT2dnZac2aNerdu7eOHDliKX/99dctI8D279+vp59+WvHx8Ro9erQkKTMzU+Hh4WrVqpXCw8M1adKkPPtq3769vvrqK0nSTz/9pK5du1olTiIjI+Xu7i6z2azY2FhVr179lvHfLcnJyXo6MFBDS5fSyPvry5R136PMe/V0YKCWff/vvu9/1b9/f82cOVOStHz5cr388svas2ePpXzp0qV64IEH9PXXX2vWrFmFkngCAAAACoo1TgAAwF01efJkmUymfD8///zzbbcbHR2tHj16aOLEierQoUOhxz127FglJCRYPidPniz0Pu4lGSYbfRK+QJ+EL9DCJUuLOpwiNWVyiC4dLS0X20qW9RRMJpOcbSrp8rHSmjI5pND7DAsL06BBgwq93dx0795ddnY3fn/1yCOP6MSJE8rMzMy1rp+fn2bPnq133nlHhmFIkjZs2KCKFSvqvffeU1hYWJ7H/tWlS5dUpkwZq32hoaEaOnSo+vbtq7CwsDs4qzv3QUiIhpYupbaVre9720oVNdTNVR+E/Lvve34uX76c67197bXX1LJlSy1fvvzOTwAAAAC4DYw4AQAAd9WIESPUp0/+a19kT8VTUAcPHlRAQICGDBmi8ePH30F0eXN0dJSjo+Ndafte1Kbr/40C+LdP2/XdtxvkbNMw1zInU0V99+0GvTvz7ULrb9euXbp48WKBRg3kZ+TIkVbP0zfffHPLZ3P27Nnq0qWLbGzy/j1W06ZNdf78ecXHx6tChQoKDQ3VwIED9cADD6hMmTLatGlTnsnPyMhI+fv7KyUlRadPn9ayZcssZZcuXdL69es1b948xcXF6bHHHtPkyZPzjeVu2rVhg0beXz/XsjYVK+iTDRukt/+99/2vFixYoMjISCUkJCgxMVHff/+9pSw6OlonT55Up06ddP36dc2YMSPHGlcAAADA3UTiBAAA3FXly5dX+fLlC6296OhoBQQEaMCAAQq5C7/gBu7U9TRDDlkjDv7KZDIpPe3Wv8C/HfPnz1f//v0towH+roJO2ZRt0aJFWr58ubZt25ZvvZtHHFy4cEEbN27UZ599JkkaNGiQQkND80yc3DxV14EDB9S+fXvt3btXVapU0aJFi9S5c2eVKVNGZcqUUcWKFbVhwwZ16tSpwOdQmByNTMtIk78ymUxyNAo2sqagStJ9z83NU3Vt2rRJTzzxhA4fPixnZ2eFhoaqf//+srW1VdeuXTVs2DDFxMSoXr16BY4TAAAAuBNM1QUAAIqNuLg4mc1mxcXFKSMjQ2azWWazWcnJyZJuJE3atm2rDh06aPTo0Tp37pzOnTun+Ph4SxunT59W3bp19dNPP1n2nTt3TmazWX/88Yck6bfffpPZbNalS5f+2RP8l8ietuvfOmWXnYMpz5fGhmHIziH3l+t/x9WrV7Vs2TLLr/GzR2j4+/sXSmIxr/aWLVumKVOmaOPGjbmOJrjZnj17VKFCBXl6emrhwoW6fv26/P395ePjo3feeUerV6/WxYsXNX36dEtfN48+yNawYUNVq1ZNO3bskHQjcbBlyxb5+PjIx8dHsbGxCg0NveNz/rtSTTb53vdUU+H916uk3feRI0da2vvtt99y1G3Xrp2uXbumAwcOKD09XYsWLdKCBQvk4+MjX19fpaSkaP78+Xd8XgAAAEBBMeIEAAAUGxMnTtQXX3xh2W7cuLEkacuWLWrTpo0iIiIUHx+vxYsXa/HixZZ63t7eOn78uCQpPT1dhw8fVkpKiqX8448/tlrovVWrVpJurA8QHBx8F8/o3yl72q5/65RdXbp11MrPf5WzqVKOsmvG/9SzW8dC6ysiIkJ+fn6qW7eupBsjNG5ePPxO5dbe8uXLNX78eEVGRqpatWr5Hr9//36NGjVKr732mkwmk+bPn6+vvvrKalTIk08+qcWLF+v111/X66+/btkfHh5u1dapU6f0+++/q3bt2vr5558VHx+vM2fOWKaLunTpkry9vRUfHy9PT887O/G/oVnHjooy71XbShVzlEX977yadfz33vcPP/zQqvzrr7+22t63b5+Sk5Pl4+OjVatWqUaNGtq9e7el/MCBA2rXrp2mTZsme3v7Ozs5AAAAoABMRkFW6wMAAPiXS0xMlLu7uxISEuTm5lbU4RSqT8IXqOVjPQu93ai1q2RrZMrF0UHPBuW+zk1xv67Xrl1TbGysqlevLicnpwIdk5ycrLatAnX5WGk5mSrKZLoxAuWa8T+VrZmkzVu/l6ura6HE17JlSw0cOFDPPfdcvvVMJpOSkpKs+r15X5s2bRQXF2d1D8aOHaunn346R1v29vaqVKmSypUrZ9m3adMmlStXTsHBwYqMjFT58uWVkpKiChUqaOjQoerfv79+/PFHde/eXadPn7aaXmrVqlWaOHGi9u3bZ9VPeHi4Ro0aJR8fHxmGoevXr2vEiBEaPny4hg8fLhcXF7333ntWx/Ts2VMtW7bU6NGjC3YBC1FycrKeDgzUUDdXtalYwXLfo/53Xp8mJmvZ9//O+56byZMn66OPPlKVKlVkGIZMJpMmTJigXr16qXPnzurSpYteeuklq2MaN26sCRMm6IknnijQNcrL33mmAQAA8O9D4gQAAKAAivsL/jtxtxIn2batWanng3N/gVrcr+vffcmanJysKZND9N23G3Q97cb0XF26ddSkyW8U2stzFD/Jycn6ICREuzZskKORqVSTjZp17Kj/vMF9Ly5InAAAAKAgmKoLAAAAKGSurq56d+bbenfm20UdCv5Brq6umvD229Lb3HcAAACgJGNxeAAAANxVqRmZ/+rF4gEAAAAAJQsjTgAAAHBXte/RS9K/d7F4AAAAAEDJwogTAAAAAAAAAACALCROAAAAAAAAAAAAspA4AQAAAAAAAAAAyELiBAAA4F/OxdFB29as1LY1KxW1dtVd6yfDZPOvWSQ+OTlZ4155Qw/Ve0QP+TZT03qPaNwrbyg5OblQ2r9y5Yr8/f0tn9q1a8vOzk6XLl2y1OnSpYul3GQyyc/PT/7+/mrZsqUkyWQyWeJp06aNatSoYdVmZGRkrn0PHDhQderUkb+/v1q1aiWz2WwpCw4OVtWqVdW4cWPVrl1bLVu21MKFCy3lx48fl52dnaWPevXq6a233srzPH18fFS3bl1L3b59++rq1atWdebPny+TyaTt27ff9nUsbMnJyZo6cZwea/Ownmz3sLq2fkhTJ477R+/7zY4ePSobG5sc1zg8PFweHh7y9/dXw4YN1bZtWx05ciTXNmJjY/Xggw/K399f999/v5566ildvnxZknTw4EGreHx8fFS2bFnLsevXr1eTJk3k5+enRx55RPv27bOU3fydq1Wrljp06KC1a9fmee7Z3y1/f3/5+fmpVatWOnToUIHOFwAAALhdJsMwjKIOAgAAoLhLTEyUu7u7EhIS5ObmVtTh3DWfhC9Qy8d63tU+tq1ZqeeD+0sq/tf12rVrio2NVfXq1eXk5FSgY5KTk9W5TVdVPl1dVR2qy2QyyTAMnUo/rrNVjmld1Fq5uroWapwzZ87U1q1b9e233+ZabjKZlJSUZNXvzfvatGmjMWPG6LHHHrtlX6tXr1aXLl1kZ2enNWvWaPTo0ZaX7sHBwWrSpIlGjBghSdq/f7+efvppDRkyRKNHj9bx48fVpEkTXbhwQZKUlJSkWrVqadOmTWrQoEGOvnx8fLRmzRo1bNhQhmGoW7du6ty5s1588UVLnUcffVQODg7y8fFRWFhYwS9aIUtOTlbQ4100slVVtfOrZrnvm/bH6cMfTmnJN9/94/f9jTfe0M6dO3XixAkdPXpUJpNJ0o3EyZo1a/TVV19Jkl599VUdOHBA3333XY42UlNTlZmZKWdnZ0nSqFGjZGNjo/fffz9H3REjRshkMmnOnDm6fPmyatWqpW3btqlevXraunWrXnzxRR04cECScnznfvjhBz399NOaO3euevXqlaPtv363ZsyYoT179igiIuKW53uzv/NMAwAA4N+HEScAAABAIZo25W1VPl1dXo41LC9uTSaTvByqq/KZGpo25e1C7zMsLEyDBg0q9HZz0717d9nZ2UmSHnnkEZ04cUKZmZm51vXz89Ps2bP1zjvvKLffayUlJckwjAIlzVJTU3X16lWVKVPGsu/QoUOKjY3VokWLtHLlSiUlJf3Ns7pzH8yYppGtqqp9I2+r+96+kbdGtqqqD2ZMK/Q+87vvGRkZ+uKLLzR37ly5urpqy5YtebbTtm1bnThxItcyR0dHS9IkIyNDycnJsrHJ+d/I1NRUffnll5Z4jh49qgoVKqhevXqSpNatW+vEiRPau3dvrv20atVKU6ZM0fTp0/M+4SyGYejKlStW34XbOV8AAADgVkicAAAAAIUocs0mVXWonmtZVXsfRa7ZVKj97dq1SxcvXizQaJH8jBw50mrapePHj9/ymNmzZ6tLly65vkjP1rRpU50/f17x8fGS/m+6qfvvv1/Vq1fX0KFD5eXllefxTz75pPz9/VWxYkWZTCb17t3bUhYaGqpnn31W9913n9q2baulS4tuGrgff9ikdn7Vci1r51dNP/7wz9739evX67777lODBg00aNAghYaG5lovMzNTK1euVJ8+ffLsKy0tTf7+/ipfvrz++OMPTZw4MUedFStWqHr16vL395ck1apVS/Hx8dq9e7ckaeXKlUpOTs73e9W0aVNFR0fnWT59+nT5+/uratWqWrRokcaNG3fb5wsAAAAUBIkTAAAAFNiZUye1clG4Vi4K15mTcUUdTvGUbsp1iiDpxggEpede9nfNnz9f/fv3t4wC+bs+/PBDmc1my8fHxyff+osWLdLy5cv1ySef5FvvryNNPDw8ZDab9dtvv+ns2bNas2aNVq9enefxX331lcxmsy5evKjq1avrtddekySlp6dr4cKFGjhwoCQV+ctyJ1vle9+dbAu3v1vd99DQUMu1eeaZZ7R27VrL2iSSFBkZaUmGbN68WS+88EKefTk4OMhsNut///uf6tSpo48//jjXeG4e/eLu7q6vv/5ar7/+uh588EFFRUWpfv36sre3z7OfW80i/frrr8tsNuv06dOaPHmynnzyyQKfLwAAAHA7SJwAAADglq79+aemjhiohW8MUcP4rbo/fqsWT3heU14I1p8pKfkeG7V2lWXx+W1rVsrF0eEfirqI2Bt5vgA2DEOyL7wlBq9evaply5ZZXhhnvwz39/dXSEjIHbefV3vLli3TlClTtHHjRlWoUCHfNvbs2aMKFSrI09MzR1nZsmXVoUMHff/991YLn/fsmXOdHTs7O/Xq1Uvr16+XJK1Zs0ZXrlxRYGCgfHx89OKLL+qXX36xrKHxT7uWkfeLf8MwdC2j8Pq61X0/f/68vvvuO7311lvy8fHRgw8+qPT0dH355ZeWNtq3by+z2axTp06pdu3alsTJzSOPfvvtN6t+HRwc9Nxzz2nhwoVW+0+cOKGdO3eqb9++VvtbtWqlqKgo/fLLL5oxY4bOnDljmborN3v27FHDhg0l/d9II39/f128eDFH3T59+uiXX35RfHx8gc4XAAAAuB139rM0AAAA/Cu8PWqoXnjQVW0atLLs6/ZQLW2LOa1pLw/W1M/yfkFpa2RaFoP/N2j/WDv9sfi4vHKZrutU+nG1f6xdofUVEREhPz8/1a1b90bfWS/DC0tu7S1fvlzjx49XZGSkqlXLfWqqbPv379eoUaP02muv5ToaIzU1VTt27NDTTz9tGYmSn82bN6tOnTqSbowwmDVrloYNG2YpHzNmjObPn5/rwuV328Ot2mnT/uNq38g7R9mm/XF6uNU/d99nzpypxx9/3GrqsjVr1mjixIl68cUXrdpycXHR559/rjp16mjv3r368MMPrcrj4uJUrlw5lSpVSpmZmVq+fLn8/Pys6oSFhalnz57y8PCw2n/27FlVrlxZkjR16lQFBATI19c313Patm2bJk+erLlz50qSZeH6vGzatEnlypVTuXLl9P777xf4fAEAAICCYMQJAAAA8nXs9yPyzLysNg1yrkPRst59qmyTpD8OxRRBZMXTuEljdbbKMZ1Mi7WMQDAMQyfTYnW2yjGNmzS20PoKDQ0ttEXh/7rGybJly3Kt169fP127dk09evTIdURA9joU2aMYxo4dq9GjR1vKbx5Z4u/vr0aNGmn48OF5xpU98qBBgwaKiYnR7NmzdebMGW3evFlPPfWUVd1nn31WCxcuVFpa2h1ejdv3n1fH6cMfTily3wmr+x6574Q+/OGU/vPquFu0UHC3uu/z589Xv379rPZ16tRJp06dynVx9ipVqmjMmDGaNGlSjrIDBw6oWbNm8vPzk5+fny5cuGCVXDEMQ+Hh4bnGM2HCBNWtW1e+vr46ceJEjqnUsr9ztWrV0pQpU/TZZ5+pV69eeZ5X9nerUaNGmjp1qr7++mvZ2Njc9vkCAAAAt2IybjWRLAAAAJSYmCh3d3clJCTIzc2tqMO5az4JX6CWj1lPk7Tgw3fVxemYWuaSOJGkXYdOa0VCVQ0cnXtCYNualXmOOCnu1/XatWuKjY1V9erV5eTkVODjkpOTNW3K2zcWgk83SfaG2j/WTuMmjZWrq+tdjBhFKTk5WR/MmKYff9gkJ9sb03c93Kqd/vPqOO57MfF3n2kAAAD8uzBVFwAAAPJ1PT1NjqXzXtna0d5WGen//C/8izNXV1dNezdE096983VGUHK4urpqwpvTijoMAAAAAHeIqboAAACQrwdatte3+87kWb7afFoPtGqfY3/2ovD3/GLwAAAAAIB7CiNOAAAAkK9GTR9S+MxrOnbuimpU8rAqO/6/K9oZd1W9H2me47h/26LwAAAAAIB7A4kTAAAA5MtkMumNOWF6/oX+6ljbQ0829ZZJ0tc/x2ldzEVNmLdAJpPJUj9q7SrZGpmMNAEAAAAAlEgkTgAAACBJWrhkqTJMuc/k6lmxomZFrNPOLZsUsmG1DMPQwx2f0+y3OsjGxvoYRpoAAAAAAEoyEicAAACQJKWkpqnNYz3zLLexsVGLdh3Uol2HfzAqAAAAAAD+WSwODwAAAAAAAAAAkIXECQAAAFDIkpOTNenVcWp7/8NqX+9htWn4sCa9Ok7JycmF0v6VK1fk7+9v+dSuXVt2dna6dOlSjrpt2rSRg4ODzp8/b9l37Ngx2djY6Mknn5QkHT9+XOXLl8+1r8zMTI0ZM0YNGzZU3bp1NWjQIKWlpVmOs7Ozs4rl6NGjlmNNJpP8/PzUqFEj1a5dW0FBQTp48GCe5+Xj46O6devK399f9erVU9++fXX16lWrOvPnz5fJZNL27dsLfsEAAAAA4DaQOAEAAAAKUXJysnoGdJHT2hN61SVAY9zb6bVSAXJae0I9A7oUSvLEw8NDZrPZ8hk6dKg6d+6ssmXL5lrfz89PCxcutGzPnz9fDz74YIH6Cg0N1f79+7V3717FxMRIkmbPnp1nLDVr1rQ6fufOndq3b58OHTqkgIAAtWjRQrGxsXn299VXX8lsNuvgwYNKTExUeHh4jnjatGmj0NDQAsUPAAAAALeLxAkAAABQiN59c5oCrnqpcWlvmUwmSTdGXjQu7a2Aq156981phd5nWFiYBg0alGf5c889Z0lAZGZmatmyZerbt2+B2t63b5/at28vBwcHmUwmdenSxSoJU1A2NjYaMmSIOnXqpI8++uiW9VNTU3X16lWVKVPGsu/QoUOKjY3VokWLtHLlSiUlJd12HAAAAABwKyROAAAAgEL0w7pN8netlmuZv2s1/bBuU6H2t2vXLl28eFGPPfZYnnW8vb1VsWJF/fjjj9qwYYOaNGlilZDIT9OmTbVq1SolJSUpLS1NS5cu1fHjxy3liYmJatq0qR544AG9+eabysjIuGV70dHReZY/+eST8vf3V8WKFWUymdS7d29LWWhoqJ599lndd999atu2rZYuXVqgcwAAAACA20HiBAAAAChEttdlGWnyVyaTSXbXC7e/+fPnq3///rKzs8u33sCBAxUaGqrQ0FANHDiwwO33799fgYGBatWqlQICAtSgQQPZ29tLkipXrqxTp05pz549ioyM1LZt2/Tee+/l255hGPmWZ0/VdfHiRVWvXl2vvfaaJCk9PV0LFy60xD5o0CCm6wIAAABwV5A4AQAAAApRhl3eyQHDMHQ9//zGbbl69aqWLVtmSSZERkZaFmkPCQmxqvvEE0/o+++/1759+9SuXbs828we8eHv76+LFy/KZDJp4sSJ+vXXX7V9+3bVrVtX9evXlyQ5OjqqQoUKkqSyZctq4MCB2rZtW74x79mzRw0bNrRa4L5nz5456tnZ2alXr15av369JGnNmjW6cuWKAgMD5ePjoxdffFG//PKLDhw4UPALBgAAAAAFUIj/bQMAAEBJ5uLooG1rVhZKO/9mrTq3k3ntCTUu7Z2jzJwcp1Zd805a3K6IiAj5+fmpbt26kqT27dvLbDbnWtfJyUkffPCBXFxcZGOT9++nvvrqK6vta9eu6dq1a/Lw8NCFCxc0ffp0TZ06VZJ0/vx5lSlTRvb29kpNTdWKFSvUuHHjXNvNzMzU/PnztX79eu3du9eyqHx+Nm/erDp16ki6MU3XrFmzNGzYMEv5mDFjNH/+fL3//vv5tgMAAAAAt4PECQAAACRJzwb1KeoQ7gmvTBynnlFdpKQba5qYTCYZhiFzcpw2u57UyomfFFpfoaGh+S4K/1dPPPFErvuvX78uJyenXMsSEhLUunVr2draKiMjQ6NGjVK3bt0kSdu3b9fEiRNla2ur69evKyAgQG+88YbV8c2bN5d0IwHzwAMPaMeOHapevXqeMT755JNycnJSenq6fHx89PHHH+vMmTPavHmzvvjiC6u6zz77rNq3b6/p06fLweHfnbADAAAAUHhMxq0mGQYAAIASExPl7u6uhIQEubm5FXU494zifl2vXbum2NhYVa9ePc/EQm6Sk5P17pvT9MO6TbK7Ll23uzES5ZWJ4+Tq6noXI/57IiIi9NFHH2nLli1FHQpwV/3dZxoAAAD/LqxxAgAAio2QkBA1b95cLi4u8vDwyFG+b98+BQUFycvLS87OzqpXr55mz56db5uXLl3SSy+9pDp16sjFxUXVqlXTyJEjlZCQcJfOApBcXV01ZcY0bfntR22M+VFbfvtRU2ZMK5ZJk6CgIL311lt6++23izoUAAAAACgWmKoLAAAUG2lpaXrqqafUrFkzhYaG5ij/5Zdf5OnpqUWLFsnLy0s7d+7U0KFDZWtrqxEjRuTa5pkzZ3TmzBnNnDlT9evX14kTJzRs2DCdOXMmx1oOwL/RkiVLijoEAAAAAChWmKoLAAAUO+Hh4Ro1apSuXLlyy7ovvviiYmJitHnz5gK3HxERoWeeeUZXr16VnV3BfkdS3KeUKqmK+3VlWh/g3sIzDQAAgIJgxAkAACjREhISVLZs2ds+xs3NLd+kSWpqqlJTUy3biYmJfztGlHz81gi4N/AsAwAAoCBInAAAgBJr165dWr58udauXVvgYy5evKipU6fq+eefz7fe22+/rSlTptxpiCjh7O3tZTKZFB8fL09PT5lMpqIOCcDfZBiG4uPjZTKZZG9vX9ThAAAAoBgjcQIAAO6qyZMn3zIBsWfPHjVp0uS22o2OjlaPHj00ceJEdejQoUDHJCYmqmvXrqpfv74mTZqUb92xY8dq9OjRVsd6eXndVowo+WxtbVW1alWdOnVKx48fL+pwANwhk8mkqlWrytbWtqhDAQAAQDFG4gQAANxVI0aMUJ8+ffKt4+Pjc1ttHjx4UAEBARoyZIjGjx9foGOSkpLUqVMnubq6auXKlbf8tbGjo6McHR1vKy7cm1xdXVWrVi2lp6cXdSgA7pC9vT1JEwAAANwSiRMAAHBXlS9fXuXLly+09qKjoxUQEKABAwYoJCSkQMckJiYqMDBQjo6OWr16NQsC47bZ2tryshUAAAAA/iVsijoAAACAbHFxcTKbzYqLi1NGRobMZrPMZrOSk5Ml3UiatG3bVh06dNDo0aN17tw5nTt3TvHx8ZY2Tp8+rbp16+qnn36SdGOkSceOHXX16lWFhoYqMTHRclxGRkaRnCcAAAAAACi+GHECAACKjYkTJ+qLL76wbDdu3FiStGXLFrVp00YRERGKj4/X4sWLtXjxYks9b29vy/oT6enpOnz4sFJSUiRJv/zyi3788UdJkq+vr1V/sbGxtz1NGAAAAAAAuLeZDMMwijoIAACA4i4xMVHu7u5KSEiQm5tbUYdzz+C6AgAAAACKG6bqAgAAAAAAAAAAyELiBAAAAAAAAAAAIAuJEwAAAAAAAAAAgCwkTgAAAAAAAAAAALKQOAEAAAAAAAAAAMhC4gQAAAAAAAAAACALiRMAAAAAAAAAAIAsJE4AAAAAAAAAAACykDgBAAAAAAAAAADIQuIEAAAAAAAAAAAgC4kTAAAAAAAAAACALCROAAAAAAAAAAAAspA4AQAA+bp+/boiIyP1ySefKCkpSZJ05swZJScnF3FkAAAAAAAAhc+uqAMAAADF14kTJ9SpUyfFxcUpNTVVHTp0UOnSpTVjxgxdu3ZNH3/8cVGHCAAAAAAAUKgYcQIAAPL08ssvq0mTJrp8+bKcnZ0t+3v27KlNmzYVYWQAAAAAAAB3ByNOAABAnrZv364dO3bIwcHBar+3t7dOnz5dRFEBAAAAAADcPYw4AQAAecrMzFRGRkaO/adOnVLp0qWLICIAAAAAAIC7i8QJAADIU4cOHTRr1izLtslkUnJysiZNmqQuXboUXWAAAAAAAAB3ickwDKOogwAAAMXTmTNn1LZtW9na2ur3339XkyZN9Pvvv6t8+fL64YcfVKFChaIO8R+TmJgod3d3JSQkyM3NrajDuWdwXQEAAAAAxQ1rnAAAgDxVqVJFZrNZS5Ys0d69e5WZmalBgwapX79+VovFAwAAAAAA3CsYcQIAAFAAjIy4O7iuAAAAAIDihhEnAADAyurVqwtct3v37ncxEgAAAAAAgH8eiRMAAGDl8ccft9o2mUz66wBVk8kkScrIyPinwgIAAAAAAPhH2BR1AAAAoHjJzMy0fDZs2CB/f3+tW7dOV65cUUJCgtatW6cHHnhA69evL+pQAQAAAAAACh0jTgAAQJ5GjRqljz/+WC1atLDsCwwMlIuLi4YOHaqYmJgijA4AAAAAAKDwMeIEAADk6ejRo3J3d8+x393dXcePH//nAwIAAAAAALjLSJwAAIA8NW3aVKNGjdLZs2ct+86dO6f/9//+nx566KEijAwAAAAAAODuIHECAADyNH/+fJ0/f17e3t7y9fWVr6+vqlWrprNnzyo0NLSowwMAAAAAACh0rHECAADy5Ovrq/3792vjxo06dOiQDMNQ/fr11b59e5lMpqIODwAAAAAAoNCZDMMwijoIAACA4i4xMVHu7u5KSEiQm5tbUYdzz+C6AgAAAACKG0acAACAPL355pv5lk+cOPEfigQAAAAAAOCfQeIEAADkaeXKlVbb6enpio2NlZ2dnWrWrEniBAAAAAAA3HNYHB4AAOTp119/tfocOHBAZ8+eVbt27fSf//yn0PsLCQlR8+bN5eLiIg8Pjxzl+/btU1BQkLy8vOTs7Kx69epp9uzZt2z3+eefV82aNeXs7CxPT0/16NFDhw4dKvT4AQAAAABAyUfiBAAA3BY3Nze9+eabmjBhQqG3nZaWpqeeekrDhw/PtfyXX36Rp6enFi1apOjoaL3xxhsaO3as5s6dm2+7Dz74oMLCwhQTE6Pvv/9ehmGoY8eOysjIKPRzAAAAAAAAJRuLwwMAgNu2fft2devWTZcvX74r7YeHh2vUqFG6cuXKLeu++OKLiomJ0ebNmwvc/v79+9WoUSP98ccfqlmzZoGOYRHzu4PrCgAAAAAobljjBAAA5OnDDz+02jYMQ2fPntXChQvVqVOnIorKWkJCgsqWLVvg+levXlVYWJiqV68uLy+vPOulpqYqNTXVsp2YmHhHcQIAAAAAgJKBxAkAAMjTBx98YLVtY2MjT09PDRgwQGPHji2iqP7Prl27tHz5cq1du/aWdT/66CO9+uqrunr1qurWrauNGzfKwcEhz/pvv/22pkyZUpjhAgAAAACAEoCpugAAwF01efLkWyYg9uzZoyZNmli2CzJVV3R0tNq2bauRI0dq/Pjxt4wjISFB58+f19mzZzVz5kydPn1aO3bskJOTU671cxtx4uXlxZRShYypugAAAAAAxQ0jTgAAQJ4GDhyo2bNnq3Tp0lb7r169qpdeeknz58+/ZRsjRoxQnz598q3j4+NzW3EdPHhQAQEBGjJkSIGSJpLk7u4ud3d31apVS4888ojKlCmjlStXKigoKNf6jo6OcnR0vK24AAAAAABAyUfiBAAA5OmLL77Q9OnTcyRO/vzzTy1YsKBAiZPy5curfPnyhRZTdHS0AgICNGDAAIWEhPztdgzDsBpRAgAAAAAAIJE4AQAAuUhMTJRhGDIMQ0lJSVbTWWVkZOi7775ThQoVCr3fuLg4Xbp0SXFxccrIyJDZbJYk+fr6ytXV1TI9V8eOHTV69GidO3dOkmRraytPT09J0unTp9WuXTstWLBADz30kI4dO6Zly5apY8eO8vT01OnTp/XOO+/I2dlZXbp0KfRzAAAAAAAAJRuJEwAAkIOHh4dMJpNMJpNq166do9xkMt2VhdMnTpyoL774wrLduHFjSdKWLVvUpk0bRUREKD4+XosXL9bixYst9by9vXX8+HFJUnp6ug4fPqyUlBRJkpOTk7Zt26ZZs2bp8uXLqlixolq1aqWdO3feleQPAAAAAAAo2VgcHgAA5LB161YZhqGAgAB9/fXXKlu2rKXMwcFB3t7eqlKlShFG+M9jEfO7g+sKAAAAAChuGHECAAByaN26tSQpNjZW1apVk8lkKuKIAAAAAAAA/hkkTgAAgJX9+/erYcOGsrGxUUJCgn777bc86/r5+f2DkQEAAAAAANx9JE4AAIAVf39/nTt3ThUqVJC/v79MJpNym9nTZDIpIyOjCCIEAAAAAAC4e0icAAAAK7GxsfL09LT8DQAAAAAA8G9C4gQAAFjx9vbO9W8AAAAAAIB/AxInAAAgT6tXr851v8lkkpOTk3x9fVW9evV/OCoAAAAAAIC7h8QJAADI0+OPP57rGifZ+0wmk1q0aKFvvvlGZcqUKaIoAQAAAAAACo9NUQcAAACKr40bN6pp06bauHGjEhISlJCQoI0bN+qhhx7SmjVr9MMPP+jixYsaM2ZMUYcKAAAAAABQKBhxAgAA8vTyyy/r008/VfPmzS372rVrJycnJw0dOlTR0dGaNWuWBg4cWIRRAgAAAAAAFB5GnAAAgDwdPXpUbm5uOfa7ubnp2LFjkqRatWrpwoUL/3RoAAAAAAAAdwUjTgAAQJ4efPBBvfLKK1qwYIE8PT0lSfHx8Xr11VfVtGlTSdLvv/+uqlWrFmWYJc7CJUuVkpqWa5mLo4OeDerzD0cEAAAAAACykTgBAAB5Cg0NVY8ePVS1alV5eXnJZDIpLi5ONWrU0KpVqyRJycnJmjBhQhFHWrKkpKap5WM9cy3btmblPxxN8XXs2DFt2bJFgwYNKupQAAAAAAD/IiROAABAnurUqaOYmBh9//33OnLkiAzDUN26ddWhQwfZ2NyY8fPxxx8v2iBxz5o9e7a8vLy0Z88eywgnAAAAAADuNhInAAAgXyaTSZ06dVKnTp2KOhTcQ9LS0tSxY0dlZGRIkrp27aqTJ0+qcePGWrlypYKDgxUdHa1SpUppzZo1+vzzz+Xr62s5/tq1a3r++efVq1cvvf/++7py5YouXryomjVrqlmzZnJ0dFSTJk302GOPqXTp0nrwwQclSb6+vvr88881YMAAffTRRypVqlSRnD8AAAAAoPgicQIAAPK1adMmbdq0SefPn1dmZqZV2fz584soqpIvIyNDc6eH6PdDMap3v5+G/b9XZWtrW9Rh/WMcHBwUFRWlxMREubu7a+3atbp06ZLGjBljqRMWFqaGDRvmevyyZcvUpUsXde/eXd27d1dUVJTWrFmjmTNnSpImT55sqVunTh1FRUVZHf/EE09o0aJFev755wv93AAAAAAAJZtNUQcAAACKrylTpqhjx47atGmTLly4oMuXL1t98Pd8veRL9aztI4cN6xWUfEUZ367S47V9NCdkalGHVmKsWrVKLVq0+NvHBwQEaPXq1YUYEQAAAADgXsGIEwAAkKePP/5Y4eHhevbZZ4s6lHvGe++9p8RfzdrevYtc7G78U+xpXykp3U/dF4Ur6fJjej64fxFHWTw899xzlqm01qxZI1dXV0vZsWPHVLly5QK1c/jwYbVp00aS1KJFC7311lsqXbq04uPjCz1mAAAAAEDJR+IEAADkKS0tTc2bNy/qMO4p/506TesC21qSJtlK29srrPWj6r7yGymcKdCk/KfqcnBwkI1NwQZP5zZVFwAAAAAAeWGqLgAAkKfBgwfryy+/LOow7hkZGRnysLNVZReXXMt9SpeWu92/Z52TO1G1alWdP3/+bx+fnJwsT0/PQowIAAAAAHCvYMQJAADI07Vr1/Tpp58qMjJSfn5+sre3typ///33iyiykikjI0MOt1gA3qGAoyj+DW6eqmvWrFny9/e3lPXo0UORkZHq27fvLdu5eaquSpUqaenSpdq0aZMee+yxuxE2AAAAAKCEMxmGYRR1EAAAoHhq27ZtnmUmk0mbN2/+B6MpWomJiXJ3d1dCQoLc3Nz+djs13Ty0+/GucswlgXI1PV3NV69TbOKVO4i0ZPm71zU1NVWDBw/WwoUL/1a/AwYM0H//+1+rdVMAAAAAAJBInAAAABRIYSVOGtW7Xx3sbTTxwUY5ysb+uFf/q1VbX61acSehliiFdV0BAAAAACgszAUBAAAK5NSpUzp9+nRRh1GiLVyyVINfGaNV5y8peMt2HbpyRemZmYq+fFl9N/2grVf//FclTQAAAAAAKI5InAAAgDxlZmbqzTfflLu7u7y9vVWtWjV5eHho6tSpyszMLOrwSpyU1DS16/6EVh04qHrDX1Tf7XvU5Jvv1PeHH9Vz0gQdOHG0qEMEAAAAAOBfj8XhAQBAnt544w2FhoZq+vTpevTRR2UYhnbs2KHJkyfr2rVrCgkJKeoQS6x+g4aq36Chilq7SqUd7PRsUJ+iDgkAAAAAAIg1TgAAQD6qVKmijz/+WN27d7fav2rVKr3wwgv/qqm7CmMtjoVLliolNc1qn4ujw786acIaJwAAAACA4oYRJwAAIE+XLl1S3bp1c+yvW7euLl26VAQRlWz/5gQJAAAAAAAlBWucAACAPDVq1Ehz587NsX/u3Lny8/MrgogAAAAAAADuLkacAACAPM2YMUNdu3ZVZGSkmjVrJpPJpJ07d+rkyZP67rvvijo8AAAAAACAQseIEwAAkKfWrVvryJEj6tmzp65cuaJLly7piSeeUHR0tMLCwoo6PAAAAAAAgELH4vAAAOC27du3Tw888IAyMjKKOpR/DIuY3x1cVwAAAABAccOIEwAAUGyEhISoefPmcnFxkYeHR47yffv2KSgoSF5eXnJ2dla9evU0e/bsArdvGIY6d+4sk8mkb775pvACBwAAAAAA9wzWOAEAAMVGWlqannrqKTVr1kyhoaE5yn/55Rd5enpq0aJF8vLy0s6dOzV06FDZ2tpqxIgRt2x/1qxZMplMdyN0AAAAAABwjyBxAgAAio0pU6ZIksLDw3MtHzhwoNV2jRo1tGvXLq1YseKWiZN9+/bp/fff1549e1S5cuVCiRcAAAAAANx7SJwAAIAcnnjiiXzLr1y58s8EUgAJCQkqW7ZsvnVSUlIUFBSkuXPnqlKlSgVqNzU1VampqZbtxMTEO4oTAAAAAACUDCROAABADu7u7rcs79+//z8UTd527dql5cuXa+3atfnW+89//qPmzZurR48eBW777bfftoyAAQAAAAAA/x4kTgAAQA5hYWGF1tbkyZNvmYDYs2ePmjRpclvtRkdHq0ePHpo4caI6dOiQZ73Vq1dr8+bN+vXXX2+r/bFjx2r06NGW7cTERHl5ed1WGwAAAAAAoOQhcQIAAO6qESNGqE+fPvnW8fHxua02Dx48qICAAA0ZMkTjx4/Pt+7mzZt19OhReXh4WO3v1auXWrZsqaioqFyPc3R0lKOj423FBQAAAAAASj4SJwAA4K4qX768ypcvX2jtRUdHKyAgQAMGDFBISMgt67/++usaPHiw1b77779fH3zwgbp161ZocQEAAAAAgHsDiRMAAFBsxMXF6dKlS4qLi1NGRobMZrMkydfXV66uroqOjlbbtm3VsWNHjR49WufOnZMk2draytPTU5J0+vRptWvXTgsWLNBDDz2kSpUq5bogfLVq1VS9evV/7NwAAAAAAEDJQOIEAAAUGxMnTtQXX3xh2W7cuLEkacuWLWrTpo0iIiIUHx+vxYsXa/HixZZ63t7eOn78uCQpPT1dhw8fVkpKyj8aOwAAAAAAuDeYDMMwijoIAACA4i4xMVHu7u5KSEiQm5tbUYdzz+C6AgAAAACKG5uiDgAAAAAAAAAAAKC4IHECAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAAAAAAAAAAAAWUicAAAAAAAAAAAAZCFxAgAAAAAAAAAAkIXECQAAAAAAAAAAQBYSJwAAAAAAAAAAAFlInAAAAAAAAAAAAGQhcQIAAAAAAAAAAJCFxAkAAAAAAAAAAEAWEicAAAAAAAAAAABZSJwAAAAAAAAAAABkIXECAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAAAAAAAAAAAAWUicAAAAAAAAAAAAZCFxAgAAAAAAAAAAkIXECQAAAAAAAAAAQBYSJwAAAAAAAAAAAFlInAAAAAAAAAAAAGQhcQIAAAAAAAAAAJCFxAkAAAAAAAAAAEAWEicAAAAAAAAAAABZSJwAAAAAAAAAAABkIXECAAAAAAAAAACQhcQJAAAAAAAAAABAFruiDgAAAOBetHDJUqWkpuXY7+LooGeD+hRBRAAAAAAAoCBInAAAANwFKalpavlYzxz7t61ZWQTRAAAAAACAgmKqLgAAAAAAAAAAgCwkTgAAAAAAAAAAALKQOAEAAChkC5csVYaJf2YBAAAAAFAS8T96AABQbISEhKh58+ZycXGRh4dHjvJ9+/YpKChIXl5ecnZ2Vr169TR79uxbttumTRuZTCarT58+d2+B9pTUNLXp2uOutQ8AAAAAAO4eFocHAADFRlpamp566ik1a9ZMoaGhOcp/+eUXeXp6atGiRfLy8tLOnTs1dOhQ2draasSIEfm2PWTIEL355puWbWdn50KPP9ufKSk6f+6sypb3lJ0d/9wCAAAAAKAk4X/yAACg2JgyZYokKTw8PNfygQMHWm3XqFFDu3bt0ooVK26ZOHFxcVGlSpUKJc687N65SyGvT1LiifPaOnep/peZrAc7Bmjw2Ndlb28vScow2eiT8AVycXTQs0F3b9QLAAAAAAD4e0icAACAEi0hIUFly5a9Zb3Fixdr0aJFqlixojp37qxJkyapdOnSedZPTU1VamqqZTsxMTHf9n+I2qqpQ1/R824t5FrZSZJkGIZ2bT6msQf7a8aXi2VjY2OZwitq7Sp9Er5AkkiiAAAAAABQjJA4AQAAJdauXbu0fPlyrV27Nt96/fr1U/Xq1VWpUiUdOHBAY8eO1b59+7Rx48Y8j3n77bctI2BuxTAMTRk9TiM9WsvRxt6y32QyqblrTZ05tVe7orbo0YB2lrKb10DZtmZlgfoBAAAAAAB3H4vDAwCAu2ry5Mk5Fmb/6+fnn3++7Xajo6PVo0cPTZw4UR06dMi37pAhQ9S+fXs1bNhQffr00VdffaXIyEjt3bs3z2PGjh2rhIQEy+fkyZN51o2JiVGlVGerpMnNAl3radWnOddsyZaakamFS5bmew74Z23atEkvvvii3nnnHRmGUejtL1++XC+99JI+/vjjQm8bAAAAAHBnSJwAAIC7asSIEYqJicn307Bhw9tq8+DBgwoICNCQIUM0fvz4247pgQcekL29vX7//fc86zg6OsrNzc3qk5f4+HiVNZzyLC9t56w/k5LzLG/fo5dSUtMKFnwJc/z4cZlMJm3ZskWSlJaWpjJlymju3LlW9Xr27ClJevPNN9WmTRv5+Pjo/vvvV5s2bbRs2TIlJibK2dlZO3fuzLWf0qVLKyAgwPK9uHTpkqQbibvsdtq0aaPt27dbHbdkyRI98sgjatmypR577DHL/oCAAK1atUpeXl4ymUySpGvXrmn48OFq06aNmjVrpoiIiBxxhIeHW87t5MmT6tatm1q3bq2AgAD99ttvlnqdOnXS559/Lj8/P8u+y5cvKygoSK1bt1aLFi0UFRWV67kOHjxYiYmJioqKko2NjWJjYyVJ586dk52dndasWSNJatKkiSTpnXfe0Z49e3JtCwAAAACQE1N1AQCAu6p8+fIqX758obUXHR2tgIAADRgwQCEhIX+7jfT0dFWuXPmOYlm4ZKlSUtN0+dIlHf7zf+qaR27lzLVLquBb9Y76KsmaNGmiFStWqG3btoqMjFStWrVy1HF2dtbFixc1ceJETZw4UZMnT1aTJk0syYxFixbpueeeU0REhJo3b57j+Dp16mjz5s2SpM8//1zDhw/XsmXLJN2Ydu3mpMjN3nnnHf30009ycHDQ5cuXLft37NihLl26aOXKlerbt68kaerUqfLz89O8efOUlJSktm3bys/PT3Xq1Mm17QEDBuiDDz5Qo0aNdOjQIT311FPau3ev7O3ttXr16hzn89JLL6l///7q3Lmzzp07p3bt2mnbtm1Wa/jExsZaknqS9OCDD+rrr7/WmDFjtGLFCj3wwAM54hg8eLBGjRqlhQsX5honAAAAAMAaI04AAECxERcXJ7PZrLi4OGVkZMhsNstsNis5+cZojejoaLVt21YdOnTQ6NGjde7cOZ07d07x8fGWNk6fPq26devqp59+kiQdPXpUb775pn7++WcdP35c3333nZ566ik1btxYjz766B3Fm5KappaP9VT3/oMkrzI6n5qQa72VKQf09MgX76ivkszb21txcXEyDEMrV660jC6RpH379km6MeJk5cq813pZtWqVQkJCtH///ltOnTV48GD98ssvysjIuGVs165d048//qiMjAyVKVPGsj8iIkIvv/yyUlNTlZSUJEnasGGDnn/+eUk3RrgMGzZMy5cvz7XduLg4ubu7q1GjRpKkunXrqlGjRtq9e3eu55ORkaEjR46oc+fOkqRKlSqpZ8+eOdbv+fbbb60SRy1btrSMoomMjFT79u1zxFKuXDmdPn26QNcDAAAAAEDiBAAAFCMTJ05U48aNNWnSJCUnJ6tx48Zq3LixZQ2UiIgIxcfHa/HixapcubLl07RpU0sb6enpOnz4sFJSUiRJDg4O2rRpkwIDA1WnTh2NHDlSHTt2VGRkpGxtbe8oXhdHB0WtXSVJen3uh/rvtd3anxRnebF/JT1FoZd2qlb3Nqpdv8Ed9VXSNWvWTD/88IPi4+OtRvp88803kqQuXbrkSBJkS0pKkmEYKlOmjFq0aGFJPuSnfPnyunDhgqQb69VkT9X1xx9/WNVbvHixZs2aJV9fX02dOlWSZBiGDh8+rAYNGqhnz56Wqa9MJpNsbP7vn89eXl46c+ZMrv2fPXtWVatajzLKrp/b+Vy4cEEVK1bMtf7NDh06JB8fH8u2jY2NKleurF9//VVubm5ycHDINZ4KFSooLi4ur8sFAAAAALgJU3UBAIBiIzw8XOHh4XmWT548WZMnT863DR8fH6sRCV5eXtq6dWshRWjt2aA++iR8gSSpUpX79MHqFVr20Tx9uzlKdoaNnMq56+nxr+ihlq3ybSdq7SqVdsz9hfe9olevXnr66afVv39/q/0bN26UJD377LM6cuSILl68qHLlylnVWb16tY4cOaJOnTrp6tWrSkpKUrNmzfLt78KFC5Yp4vKbqit7qqu0tDR17dpVhw8f1oULF3T06FF16tRJaWlpKl++vIKCgiRJmZmZluTJqVOnVKVKlVzbrVKlik6fPm2179SpU+rSpUuu5/Puu+/q/PnzOerXrl07R9tOTtbr6TzxxBMaOHCgJk2aJLPZnGs8d2OBewAAAAC4VzHiBAAA4A7cPOrEo0xZPf/GG/pk0wb9d/N6vRex7JZJE0myNTL1bFCfux1qkapVq5ZatGihJ5980rLPbDZbpktbsWKFQkJCLCNQbvb1119ry5YtWr9+vbZt26a9e/fmmwgICwtTkyZNCjSi6Pfff5d0Y2SSh4eHMjMzFRERoeXLl2v9+vXavHmzUlJSlJycrI4dO+rzzz+XJCUnJ+uTTz5R7969c23Xy8tLCQkJlgXhjxw5on379umRRx7J9XxsbGxUq1YtbdiwQZL0v//9TytXrlTXrl2t2q1Tp46OHj1qtS97rZVOnTrleZ7x8fHy8vK65fUAAAAAAJA4AQAAuCPPBvWRrZFZ1GGUCB9++KHV9FURERFq2bKlZbtdu3aKiIiwOiYpKUlnzpyxGoVSu3btHNN1HT58WAEBAQoICNCOHTv03//+11J281Rd69evtzpuzJgxat68uVq0aCFfX1/VrVtXmzdvtqxNIknNmzfXt99+q/Hjx2vv3r1q06aNOnbsqDFjxuRYGN4wDEvCJiwsTOPGjVPr1q01bNgwLV68WNeuXcvzfObMmaOwsDC1adNGvXr10ty5c60Whpek7t27a/PmzVb77Ozs9MUXX1iNRLk5josXL6pKlSqys2OwOQAAAAAUhMlg3D4AAMAtJSYmyt3dXQkJCXJzc7MqW7hkqZLSrqtN1x5/q+1ta1bq+eD+t654D8rvupZEU6ZMUaNGjfT444/ftT4GDx6s9957T+7u7nnWOXHihF5++WV98803mjFjhtq0aaOHHnrorsUEAAAAAPcSRpwAAADcIUadQLqxlsr27dsVGBh4V/v5/PPP802amM1mPfHEExo1apQk6dVXXyVpAgAAAAC3gREnAAAABXCrkRELlyxVSmqaMkw2tz3yhBEn986IEwAAAABAyceIEwAAgELwbFAfPR/cX6Ud7CyLxQMAAAAAgJKHxAkAAEAhYtouAAAAAABKNhInAAAAAAAAAAAAWeyKOgAAAIB7jYujg6LWrsp3rZOotassI1NcHB3+qdAAAAAAAMAtkDgBAAAoZM8G9dHCJUvzTZ7YGpn/2gXhAQAAAAAozpiqCwAA4C54NqhPrgvFR61dpW1rVjLKBAAAAACAYspkGIZR1EEAAAAUd4mJiXJ3d1dCQoLc3NwKfNzCJUuVkppm2XZxdNCzQX3uRogl0t+9rgAAAAAA3C1M1QUAAHAXkSQBAAAAAKBkYaouAAAAAAAAAACALCROAAAAAAAAAAAAspA4AQAAAAAAAAAAyELiBAAAAAAAAAAAIAuJEwAAAAAAAAAAgCwkTgAAAAAAAAAAALLYFXUAAAAAJYFhGJKkxMTEIo7k3pJ9PbOvLwAAAAAARY3ECQAAQAEkJSVJkry8vIo4kntTUlKS3N3dizoMAAAAAABkMvh5HwAAwC1lZmbqzJkzKl26tEwmU4GOSUxMlJeXl06ePCk3N7e7HOHfU9QxGoahpKQkValSRTY2zCILAAAAACh6jDgBAAAoABsbG1WtWvVvHevm5lZsEyfZijJGRpoAAAAAAIoTftYHAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAAAAd4mjo6MmTZokR0fHog4lTyUhRgAAAAAA/kksDg8AAAAAAAAAAJCFEScAAAAAAAAAAABZSJwAAAAAAAAAAABkIXECAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAAAAAAAAAAAAWUicAAAAAAAAAAAAZCFxAgAAAAAAAAAAkIXECQAAAAAAAAAAQBYSJwAAAAAAAAAAAFlInAAAAAAAAAAAAGQhcQIAAAAAAAAAAJCFxAkAAAAAAAAAAEAWEicAAAAAAAAAAABZSJwAAAAAAAAAAABkIXECAAAAAAAAAACQhcQJAAAAAAAAAABAFhInAAAAAAAAAAAAWUicAAAAAAAAAAAAZCFxAgAAAAAAAAAAkIXECQAAAAAAAAAAQBYSJwAA4K5JS0vTtGnT9McffxR1KECxExMTo5CQEKWnpxd1KAAAAACAm5A4AQAAd80rr7yiQ4cOydfX12q/yWTSN998k+dxx48fl8lkktlsliRFRUXJZDLpypUrdy9YlCjh4eHy8PAo6jD+toyMDAUHB2vbtm166623ijocAAAAAMBNSJwAAFCEgoODZTKZZDKZZG9vrxo1amjMmDG6evVqUYd2x1avXq39+/fr888/v+O2mjdvrrNnz8rd3b0QIiu+oqKi1KNHD1WuXFmlSpWSv7+/Fi9enKPe1q1b9eCDD8rJyUk1atTQxx9/bFX+2WefqWXLlipTpozKlCmj9u3b66effrKqM2/ePPn5+cnNzU1ubm5q1qyZ1q1bd8fnkJ30yv44ODjI19dXb731lgzDuOP2sz399NM6cuRIobWXrU2bNlax16xZU2PHjlVqamqh9vPuu+8qMDBQq1evVmRkpCVJCAAAAAAoenZFHQAAAP92nTp1UlhYmNLT07Vt2zYNHjxYV69e1bx584o6tDvSvXt3de/evVDacnBwUKVKlQqlreJs586d8vPz02uvvaaKFStq7dq16t+/v9zc3NStWzdJUmxsrLp06aIhQ4Zo0aJF2rFjh1544QV5enqqV69ekm4kYIKCgtS8eXM5OTlpxowZ6tixo6Kjo3XfffdJkqpWrarp06dbRgN98cUX6tGjh3799Vc1aNDgjs8lMjJSDRo0UGpqqrZv367BgwercuXKGjRo0B23LUnOzs5ydnYulLb+asiQIXrzzTeVlpamPXv26LnnnpMkvf3224XWx+uvv275e8eOHYXWLgAAAADgzjHiBACAIubo6KhKlSrJy8tLffv2Vb9+/SzTWKWmpmrkyJGqUKGCnJyc1KJFC+3Zs8fq+OjoaHXt2lVubm4qXbq0WrZsqaNHj0qSMjMz9eabb6pq1apydHSUv7+/1q9fbzk2LS1NI0aMUOXKleXk5CQfH598Xw4HBwfr8ccf15QpU1ShQgW5ubnp+eefV1pamqWOj4+PZs2aZXWcv7+/Jk+ebLXv7Nmz6ty5s5ydnVW9enVFRETk2W9uU3Xt2LFDrVu3louLi8qUKaPAwEBdvny5QNft8uXL6tevnzw9PeXs7KxatWopLCwsz/4Nw9CMGTNUo0YNOTs7q1GjRvrqq69yxPf999+rcePGcnZ2VkBAgM6fP69169apXr16cnNzU1BQkFJSUvLsZ9y4cZo6daqaN2+umjVrauTIkerUqZNWrlxpqfPxxx+rWrVqmjVrlurVq6fBgwdr4MCBmjlzpqXO4sWL9cILL8jf319169bVZ599pszMTG3atMlSp1u3burSpYtq166t2rVrKyQkRK6urtq9e3ee8d2OcuXKqVKlSvL29la/fv3UvHlz7d2711K+Z88edejQQeXLl5e7u7tat25tVS5JV65c0dChQ1WxYkU5OTmpYcOGWrNmjaScU3UdPXpUPXr0UMWKFeXq6qqmTZsqMjLyb8Xu4uKiSpUqqVq1aurVq5c6dOigDRs2WMpv9X2Q7uy5zB61s2LFCrVt21YuLi5q1KiRdu3aZdVHfs/A+vXr1aJFC3l4eKhcuXJ67LHHLP1LUrNmzaySN5IUHx8ve3t7bdmyRdKN56R///4qU6aMXFxc1LlzZ/3+++9/65oCAAAAQElC4gQAgGLG2dnZslj0q6++qq+//lpffPGF9u7dK19fXwUGBurSpUuSpNOnT6tVq1ZycnLS5s2b9csvv2jgwIG6fv26JGn27Nl67733NHPmTO3fv1+BgYHq3r275eXnhx9+qNWrV2v58uU6fPiwFi1aJB8fn3zj27Rpk2JiYrRlyxYtWbJEK1eu1JQpU277PCdMmKBevXpp3759euaZZxQUFKSYmJgCHWs2m9WuXTs1aNBAu3bt0vbt29WtWzdlZGQU6LpNmDBBBw8e1Lp16xQTE6N58+apfPnyefY3fvx4hYWFad68eYqOjtZ//vMfPfPMM9q6datVvcmTJ2vu3LnauXOnTp48qd69e2vWrFn68ssvtXbtWm3cuFFz5sy5reuUkJCgsmXLWrZ37dqljh07WtUJDAzUzz//nOci4ykpKUpPT7dq52YZGRlaunSprl69qmbNmt1WfAXx888/a+/evXr44Yct+5KSkjRgwABt27ZNu3fvVq1atdSlSxclJSVJupFc6Ny5s3bu3KlFixbp4MGDmj59umxtbXPtIzk5WV26dFFkZKR+/fVXBQYGqlu3boqLi7PUmTx58i2/33+1b98+7dixQ/b29pZ9t/o+3Olzme2NN97QmDFjZDabVbt2bQUFBVnauNUzcPXqVY0ePVp79uzRpk2bZGNjo549eyozM1OS1K9fPy1ZssRq+rRly5apYsWKat26taQbidKff/5Zq1ev1q5du2QYhrp06cJi9gAAAADufQYAACgyAwYMMHr06GHZ/vHHH41y5coZvXv3NpKTkw17e3tj8eLFlvK0tDSjSpUqxowZMwzDMIyxY8ca1atXN9LS0nJtv0qVKkZISIjVvqZNmxovvPCCYRiG8dJLLxkBAQFGZmZmgeMtW7ascfXqVcu+efPmGa6urkZGRoZhGIbh7e1tfPDBB1bHNWrUyJg0aZJlW5IxbNgwqzoPP/ywMXz4cMMwDCM2NtaQZPz666+GYRjGli1bDEnG5cuXDcMwjKCgIOPRRx/NNcaCXLdu3boZzz33XIHOOTk52XBycjJ27txptX/QoEFGUFCQVXyRkZGW8rffftuQZBw9etSy7/nnnzcCAwML1K9hGEZERITh4OBgHDhwwLKvVq1aOe7pjh07DEnGmTNncm3nhRdeMGrWrGn8+eefVvv3799vlCpVyrC1tTXc3d2NtWvXFji2vGTfO2dnZ6NUqVKGvb29IckYOnRovsddv37dKF26tPHtt98ahmEY33//vWFjY2McPnw41/phYWGGu7t7vm3Wr1/fmDNnjmV7zpw5RkBAQL7HtG7d2rC3tzdKlSplODg4GJIMGxsb46uvvjIMo2Dfhzt9LrOv4eeff24pj46ONiQZMTExhmHk/wzk5vz584Yk47fffrNs29nZGT/88IOlTrNmzYxXXnnFMAzDOHLkiCHJ2LFjh6X8woULhrOzs7F8+fIC9wsAAAAAJREjTgAAKGJr1qyRq6urnJyc1KxZM7Vq1Upz5szR0aNHlZ6erkcffdRS197eXg899JBlZIbZbFbLli2tfg2fLTExUWfOnLE6XpIeffRRy/HBwcEym82qU6eORo4caTUdUV4aNWokFxcXy3azZs2UnJyskydP3tZ5/3VkQ7NmzW57xEluCnLdhg8frqVLl8rf31+vvvqqdu7cmWdfBw8e1LVr19ShQwe5urpaPgsWLLCa+kiS/Pz8LH9XrFhRLi4uqlGjhtW+8+fPF+gco6KiFBwcrM8++yzHmiMmk8lq28gaNfDX/ZI0Y8YMLVmyRCtWrJCTk5NVWZ06dWQ2m7V7924NHz5cAwYM0MGDB3ONJy4uzur8p02blm/8y5Ytk9ls1r59+7Rs2TKtWrXKamqo8+fPa9iwYapdu7bc3d3l7u6u5ORkywgRs9msqlWrqnbt2vn2k+3q1at69dVXVb9+fXl4eMjV1VWHDh2yGnEyYsQIq+nK8tKvXz+ZzWbt2rVLvXv31sCBAy3rxxTk+3Cnz2W2m79PlStXliTL9ye/Z0C68Rz07dtXNWrUkJubm6pXry5Jluvh6empDh06aPHixZJurJ2za9cu9evXT5IUExMjOzs7q1FC5cqVU506dQr8nAIAAABAScXi8AAAFLG2bdtq3rx5sre3V5UqVSwvW8+ePSsp95fk2fsKsjh2fsc/8MADio2N1bp16xQZGanevXurffv2OdZrKIjsNm1sbKym/5FU4Kl9cnvxn5v8zjuvJMLN5925c2edOHFCa9euVWRkpNq1a6cXX3zRap2QbNlTG61du9aysHo2R0dHq+2bX5SbTKYcL85NJpOlvfxs3bpV3bp10/vvv6/+/ftblVWqVEnnzp2z2nf+/HnZ2dmpXLlyVvtnzpypadOmKTIy0uolfDYHBwfL4vBNmjTRnj17NHv2bH3yySc56lapUkVms9mynde0X9m8vLwsbderV0/Hjh3ThAkTNHnyZDk5OSk4OFjx8fGaNWuWvL295ejoqGbNmlnWy7ndhd9feeUVff/995o5c6Z8fX3l7OysJ5980mr9nYJyd3e3xL5o0SI1aNBAoaGhGjRoUIG+D3f6XGb76/dJ+r/v46366Natm7y8vPTZZ5+pSpUqyszMVMOGDa2uR79+/fTyyy9rzpw5+vLLL9WgQQM1atTIEk9ucosTAAAAAO41jDgBAKCIlSpVSr6+vvL29rZ6Uerr6ysHBwdt377dsi89PV0///yz6tWrJ+nGL9K3bduWa2LCzc1NVapUsTpeknbu3Gk5Prve008/rc8++0zLli3T119/bVkLJDf79u3Tn3/+adnevXu3XF1dVbVqVUk3fsmenfSRbvzCPjY2Nkc7f12EfPfu3apbt26e/d7Mz88vz5EDBblu2XEGBwdr0aJFmjVrlj799NNc26tfv74cHR0VFxcnX19fq4+Xl1eB4r0dUVFR6tq1q6ZPn66hQ4fmKG/WrJk2btxotW/Dhg1q0qSJ1ffn3Xff1dSpU7V+/Xo1adKkQH0bhqHU1NRcy+zs7KzO/VaJk7+ytbXV9evXLS/ut23bppEjR6pLly5q0KCBHB0ddeHCBUt9Pz8/nTp1SkeOHClQ+9u2bVNwcLB69uyp+++/X5UqVdLx48dvK8bc2Nvba9y4cRo/frxSUlIK9H0ojOfyVvJ7Bi5evKiYmBiNHz9e7dq1U7169SyLxt/s8ccf17Vr17R+/Xp9+eWXeuaZZyxl9evX1/Xr1/Xjjz9atXvkyJHbihMAAAAASiJGnAAAUEyVKlVKw4cP1yuvvKKyZcuqWrVqmjFjhlJSUjRo0CBJN6YemjNnjvr06aOxY8fK3d1du3fv1kMPPaQ6derolVde0aRJk1SzZk35+/srLCxMZrPZMj3PBx98oMqVK8vf3182NjaKiIhQpUqV5OHhkWdcaWlpGjRokMaPH68TJ05o0qRJGjFihGxsbvweIyAgQOHh4erWrZvKlCmjCRMm5Lqgd0REhJo0aaIWLVpo8eLF+umnnxQaGlqgazN27Fjdf//9euGFFzRs2DA5ODhoy5Yteuqpp1S+fPlbXreJEyfqwQcfVIMGDZSamqo1a9bk+TK4dOnSGjNmjP7zn/8oMzNTLVq0UGJionbu3ClXV1cNGDCgQDEXRHbS5OWXX1avXr0sI0scHBwsiYphw4Zp7ty5Gj16tIYMGaJdu3YpNDRUS5YssbQzY8YMTZgwQV9++aV8fHws7WRPKyVJ48aNU+fOneXl5aWkpCQtXbpUUVFRWr9+faGcy8WLF3Xu3Dldv35dv/32m2bPnq22bdvKzc1N0o0E18KFC9WkSRMlJibqlVdesRpF0bp1a7Vq1Uq9evXS+++/L19fXx06dEgmk0mdOnXK0Z+vr69WrFihbt26yWQyacKECTlG98ydO1crV64s0HRdN+vbt6/GjRunjz76SGPGjLnl9+FOn8uCyO8ZKFu2rMqVK6dPP/1UlStXVlxcnNU0adlKlSqlHj16aMKECYqJiVHfvn0tZbVq1VKPHj00ZMgQffLJJypdurRef/113XffferRo8dtXT8AAAAAKHGKanEVAACQc3H4v/rzzz+Nl156yShfvrzh6OhoPProo8ZPP/1kVWffvn1Gx44dDUmGJKNly5aWBckzMjKMKVOmGPfdd59hb29vNGrUyFi3bp3l2E8//dTw9/c3SpUqZbi5uRnt2rUz9u7de8t4J06caJQrV85wdXU1Bg8ebFy7ds1SJyEhwejdu7fh5uZmeHl5GeHh4bkuDv/f//7X6NChg+Ho6Gh4e3sbS5YssZTfanF4wzCMqKgoo3nz5oajo6Ph4eFhBAYGWspvdd2mTp1q1KtXz3B2djbKli1r9OjRwzh27Fie552ZmWnMnj3bqFOnjmFvb294enoagYGBxtatW/OML7fFyydNmmQ0atQo3+ubfR9v/rRu3dqqXlRUlNG4cWPDwcHB8PHxMebNm2dV7u3tnWs7N9+DgQMHGt7e3oaDg4Ph6elptGvXztiwYUOesRVU9r3L/tja2hpVq1Y1hgwZYpw/f95Sb+/evUaTJk0MR0dHo1atWkZERITh7e1tfPDBB5Y6Fy9eNJ577jmjXLlyhpOTk9GwYUNjzZo1hmHkvL6xsbFG27ZtDWdnZ8PLy8uYO3eu0bp1a+Pll1+21Jk0aZLh7e2db/x/PSZbSEiI4enpaSQlJd3y+2AY//dcuri4GKVLl76t5/Kv33/DMIzLly8bkowtW7ZY9uX3DGzcuNGoV6+e4ejoaPj5+RlRUVGGJGPlypVW57V27VpDktGqVasc53zp0iXj2WefNdzd3Q1nZ2cjMDDQOHLkSL7XDwAAAADuBSbDyGMCYwAAUKIkJyerbdu22rlzZ66LUheG4OBgXblyRd98881daR8AAAAAAKCoscYJAAD3gOPHj+vcuXOKj4/XgQMHijocAAAAAACAEovECQAA94Dw8HDVr19fFStWVO3atYs6HAAAAAAAgBKLqboAAAAAAAAAAACyMOIEAAAAAAAAAAAgS4lJnISEhKh58+ZycXGRh4dHjvJ9+/YpKChIXl5ecnZ2Vr169TR79uxbtvv888+rZs2acnZ2lqenp3r06KFDhw7lWjc1NVX+/v4ymUwym813eEYAAAAAAAAAAKC4KTGJk7S0ND311FMaPnx4ruW//PKLPD09tWjRIkVHR+uNN97Q2LFjNXfu3HzbffDBBxUWFqaYmBh9//33MgxDHTt2VEZGRo66r776qqpUqVIo5wMAAAAAAAAAAIqfErfGSXh4uEaNGqUrV67csu6LL76omJgYbd68ucDt79+/X40aNdIff/yhmjVrWvavW7dOo0eP1tdff60GDRro119/lb+/f4HbzczM1JkzZ1S6dGmZTKYCHwcAAAAAAFCcGYahpKQkValSRTY2JeY3ugAA5MmuqAO4mxISElS2bNkC17969arCwsJUvXp1eXl5Wfb/73//05AhQ/TNN9/IxcWlQG2lpqYqNTXVsn369GnVr1+/4MEDAAAAAACUICdPnlTVqlWLOgwAAO7YPZs42bVrl5YvX661a9fesu5HH32kV199VVevXlXdunW1ceNGOTg4SLrxq4ng4GANGzZMTZo00fHjxwvU/9tvv60pU6bk2H/y5Em5ubnd1rkAAAAAAAAUV4mJifLy8lLp0qWLOhQAAApFkSZOJk+enGty4WZ79uxRkyZNbqvd6Oho9ejRQxMnTlSHDh1uWb9fv37q0KGDzp49q5kzZ6p3797asWOHnJycNGfOHCUmJmrs2LG3FcPYsWM1evRoy3b2PyLc3NxInAAAAAAAgHsOU5MDAO4VRZo4GTFihPr06ZNvHR8fn9tq8+DBgwoICNCQIUM0fvz4Ah3j7u4ud3d31apVS4888ojKlCmjlStXKigoSJs3b9bu3bvl6OhodUyTJk3Ur18/ffHFF7m26ejomOMYAAAAAAAAAABQvBVp4qR8+fIqX758obUXHR2tgIAADRgwQCEhIX+7HcMwLOuTfPjhh3rrrbcsZWfOnFFgYKCWLVumhx9++I5jBgAAAAAAAAAAxUeJWeMkLi5Oly5dUlxcnDIyMmQ2myVJvr6+cnV1VXR0tNq2bauOHTtq9OjROnfunCTJ1tZWnp6ekm4s0N6uXTstWLBADz30kI4dO6Zly5apY8eO8vT01OnTp/XOO+/I2dlZXbp0kSRVq1bNKg5XV1dJUs2aNVnwDAAAAAAAAEAOGRkZSk9PL+owANzE3t5etra2BapbYhInEydOtJoWq3HjxpKkLVu2qE2bNoqIiFB8fLwWL16sxYsXW+p5e3tbFnRPT0/X4cOHlZKSIklycnLStm3bNGvWLF2+fFkVK1ZUq1attHPnTlWoUOGfOzkAAAAAAAAA94Tk5GSdOnVKhmEUdSgAbmIymVS1alXL4Ih86xo8wf+IxMREubu7KyEhgcXhAQAAAADAPYN3HsD/ycjI0O+//y4XFxd5enrKZDIVdUgAdGN5jvj4eKWkpKhWrVq3HHlSYkacAAAAAAAAAEBxlp6eLsMw5OnpKWdn56IOB8BNPD09dfz4caWnp5M4AYC/+u9/v9C5c4lFHUaJV6mSm158cUBRhwEAAAAAQLHDSBOg+Lmd55LECYB/nXPnEmVr+0BRh1HinTu3t6hDAAAAAADgnpCcnKypb4Vo/fpIZWQasrExqXOn9pow/o0CrceAkik5OVnvzgjRtq2RsrMzlJ5uUqs27fXKq9z3omZT1AEAAAAAAAAAwL9VcnKy2nfopMhth1S+2sOqVL2ZPKs9rE3bD6l9h05KTk6+o/avXLkif39/y6d27dqys7PTpUuXLHW6dOliKTeZTPLz85O/v79atmwp6cYv9bPjaNOmjWrUqGHVZmRkZK59P/fcc5a2mjZtqk2bNlnKgoODVbVqVTVu3Fi1a9dWy5YttXDhwhxtTJw4Uba2tjpx4kSe5xgeHi4PDw9LPI0aNdKqVaus6hiGIR8fH7Vv377gF+8uSk5OVs/HO8nV9pDGvfywXhvRTG+MeliutofU8/E7v+9Syb/32SZPnqwKFSpY+mzatKl27txpVScpKUmurq4aPHjw375eN2PECQAAAAAAAAAUkalvhSjdxlNlyle17DOZTPIoV1VXLt4of2f623+7fQ8PD5nNZsv2zJkztXXrVpUtW9ay77vvvrPqe+fOnfmOePjwww/12GOP3bLvDz74QB4eHpIks9ms9u3bKz4+3jJl0uuvv64RI0ZIkvbv36+nn35a8fHxGj16tCQpMzNT4eHhatWqlcLDwzVp0qQ8+2rfvr2++uorSdJPP/2krl27qkePHpbyyMhIubu7y2w2KzY2VtWrV79l/HfTuzNC1LGlpx7ws77v2dvvzgjRlDf//n2XSva9/6v+/ftr5syZkqTly5fr5Zdf1p49eyzlS5cu1QMPPKCvv/5as2bNuuMRO4w4AQAAAAAAAIAisn59pDzK3ZdrmXvZ+7R+fe6/6P+7wsLCNGjQoEJtMy/ZL86lG6Mf8ltjws/PT7Nnz9Y777wjwzAkSRs2bFDFihX13nvvKSwsTJmZmQXq99KlSypTpozVvtDQUA0dOlR9+/ZVWFjY7Z9MIdu2NVKN78/9vje+/z5t21q4910qWfc+P5cvX871/r722mtq2bKlli9f/rfjzkbiBAAAAAAAAACKSEamkedLZZPJpIyC5QoKZNeuXbp48WKBRgzkZ+TIkVbTNR0/fjzPuq+//rpq1qypJ554QhEREfm+QG/atKnOnz+v+Ph4STdehg8cOFAPPPCAypQpYzXd019FRkZapqPq1auX3n//fUvZpUuXtH79evXt21eDBg26rSTM3WJnl/99tyvkuaJK2r3/qwULFsjf31/Vq1fXuHHjNG3aNEtZdHS0Tp48qU6dOmnQoEEKDQ392+eXjcQJAAAAAAAAABQRWxtTnr+yNwxDtoX4Bnf+/Pnq37+/7O7wrfyHH34os9ls+fj4+ORZd/r06Tp69KiWL1+uV155RWlpaXnWvfk6XLhwQRs3blRQUJAk3fKFePv27WU2m3XkyBH9+OOPGjx4sM6cOSNJWrRokTp37qwyZcqoUaNGqlixojZs2HCbZ124rl/P/75fv164/ZWke5+b/v37W6ZZW758uZ544gn9+eefkm4k2Pr37y9bW1t17dpVx44dU0xMzN86v2wkTgAAAAAAAACgiHTq1F4Jl07nWpZw6bQ6dSqcxcyvXr2qZcuWaeDAgZL+b4SGv7+/QkJC7rj9W7XXvn17JSUl6bfffsuzjT179qhChQry9PTUwoULdf36dfn7+8vHx0fvvPOOVq9erYsXL2r69OmWvr7//vsc7TRs2FDVqlXTjh07JN1IGmzZskU+Pj7y8fFRbGxsoYxKuBMtW7fXr7/lft9//e20WrYuvEXsS9q9v3lUS27HtGvXTteuXdOBAweUnp6uRYsWacGCBfLx8ZGvr69SUlI0f/78OzonFocHAAAAAAAAgCIyYfwb2tqhk65cvLGmicl0YyRCwqXTssuI14TxCwuln4iICPn5+alu3bqS/m+ERmH5a3vXr19XbGysatWqJenGgu3nz59XjRo1cj1+//79GjVqlF577TWZTCbNnz9fX331lTp16mSp8+STT2rx4sV6/fXX9frrr1v2h4eHW7V16tQp/f7776pdu7Z+/vlnxcfH68yZM7KxuTGO4NKlS/L29lZ8fLw8PT0L6QrcnldefUM9H79xbo3v/7/7/utvp7VhW7xWflM4910qeff+ww8/tCr/+uuvrbb37dun5ORk+fj4aNWqVapRo4Z2795tKT9w4IDatWunadOmyd7e/m+dE4kTAAAAAAAAACgirq6uity4XlPfCtH69ZHKyJRsbW6MRJkwfqFcXV0LpZ/Q0NBCWxh85MiRGj9+vGV77Nixevrpp63qZGRkKDg4WAkJCbK1tVWpUqX01VdfWS3qPX36dH3++edKSUlRhQoVNHbsWPXv318//vijzp8/r/btrUddPPvss5o4caJGjhyZI6bsUQ83prm6rmnTpqlRo0YaPny4+vTpY0maSFLZsmXVvn17LVy4UKNHjy6Ua3K7XF1dtfKb9Xp3RoimzY6UnZ10/fqNkSgrvym8+y6VrHuflwULFigyMlKGcWNtmIULF8rT01OhoaHq16+fVd2GDRuqSpUq+vbbb/XEE0/8rfM0GQVZph53LDExUe7u7kpISJCbm1tRhwP8q02YMEe2tg8UdRglXkbGXk2d+lJRhwEAAACgiPHOA/g/165dU2xsrKpXry4nJ6eiDgfATW7n+WSNEwAAAAAAAAAAgCwkTgAAAAAAAAAAALKQOAEAAAAAAAAAAMhC4gQAAAAAAAAAACALiRMAAAAAAAAAAIAsJE4AAAAAAAAAAACykDgBAAAAAAAAgCKUnJysN94Yp0ceeViPPPKwHn74Ib3xxjglJyffcdtXrlzR/2/vzqO0qu51Ub8fRSslBUjZhYJCQDEaBVGjRCOgoJfkgJpobLaoGMVsTeJ2qwk7NihBY2yCaTyeXEvE5iRqErdu49YogkfFREUJil2iIIpwRdEqQGks6v5hWSeERigKPpDnGWONwVpzrvn9Zo2h4xv11pyzd+/eDdeuu+6a5s2bZ8GCBQ19hgwZ0tBeKBSy1157pXfv3jn44IOTJIVCoaGW/v37Z5dddllpzIcffni1n33qqac2jLXffvtl4sSJDW2nnHJKOnfunD59+mTXXXfNwQcfnFtvvbWhfdasWWnevHnDZ+y+++758Y9/vMZ5VlZWplevXg19TzjhhCxevHilPjfddFMKhUIef/zx9f9BslVpXuwCAAAAAAC2VosWLcqQIUNS2bVzBh02IIVCIXV1dXnjjVkZMmRI7r///pSWljZ6/Pbt22fatGkN91dffXUeffTRdOzYseHZ/fff3/DvQqGQKVOmrPUzf/7zn+frX//6Z372z372s7Rv3z5JMm3atBx22GGZP39+CoVCkuSHP/xhzj777CTJ9OnT861vfSvz58/Pueeeu0rtCxcuTM+ePXPUUUdljz32WO3n/e53v8uee+6Zurq6/I//8T9y880356yzzmpor6qqSv/+/VNVVZWDDjroM+tn62XFCQAAAABAkVxxxeWp7No5lZVdGwKFQqGQysquqezaOVdccXmTft748eNz2mmnNemYa/JpaJJ8svLl0/mtzl577ZXrrrsuV155Zerq6lZpX7hwYerq6tKuXbvP/NylS5dm8eLF6dChQ8Ozl19+OTNnzsxtt92Wu+++OwsXLly/ybBVEZwAAAAAABTJxIkT07Vrl9W2de3aZaXtrTbUk08+mffee2+dVouszfe+972VtuqaNWvWGvv+8Ic/TPfu3XP00UfnrrvuWmt4st9+++Wdd97J/Pnzk/zfbca+9KUvpVu3bjnjjDNSUVGxxve/+c1vpnfv3tlhhx1SKBRy7LHHNrRVVVXlpJNOyhe+8IUMGDAgv/3tb9d/4mw1BCcAAAAAAEW0pjBhbSFDY9x0000ZPnx4mjffsBMcfv7zn2fatGkNV2Vl5Rr7/uQnP8lrr72WO++8M+eff36WLVu2xr7/vNLk0626nn/++cydOzf33Xdf7r333jW+/7vf/S7Tpk3Le++9l27duuUHP/hBkmT58uW59dZbM2LEiCTJaaedlqqqqvWYMVsbwQkAAAAAQBGtbmuqtT1vjMWLF+eOO+5oCA8efvjhhhUjY8eO3eDxP2u8ww47LAsXLszzzz+/xjGefvrpbL/99ikvL1+lrWPHjhk0aFAefPDBlQ68P+qoo1bp27x583zjG9/IAw88kCS577778sEHH+Twww9PZWVlzjrrrEydOjUvvPDCBsyYzzOHwwMAAAAAFMmhhx6aN96YlcrKrqu0vfHG7Bx66KFN8jl33XVX9tprr/Tq1SvJJ0HGPx4av6H+ebyPP/44M2fOTM+ePZMkTz31VN55553ssssuq31/+vTpOeecc/KDH/xgtSttli5dmieeeCLf+ta3VjnwfnUeeeSR7Lbbbkk+2aZr3LhxOfPMMxvazzvvvNx000259tpr13OmbA0EJwAAAAAARTJq1H9kyJAhST4506RQKKSuri5vvDE7b7zxVq7/n79uks+pqqpqskPhv/e97+XCCy9suB81alS+9a1vrdSntrY2p5xySqqrq1NSUpK2bdvmd7/73UoHtv/kJz/JjTfemA8//DDbb799Ro0aleHDhze0f7qyJPkkOBkwYEC+853vrLGub37zm2ndunWWL1+eysrK3HDDDXn77bfzyCOPZMKECSv1Pemkk3LYYYflJz/5SVq2bLkhPw4+hwQnAAAAAABFUlpamvvvvz9XXHF5Hnr4/x4Ef+ihh+b6//nrlJaWNsnnPPbYY+vUb3Xbg/3js8mTJ6/TOK1atcoTTzyxxvabb755re9XVlbm448/XqfPSrLWA+o//PDDVZ7tvffeDYfQwz8TnAAAAAAAFFFpaWnGjr08TXDUCNAEHA4PAAAAAABQT3ACAAAAAABQT3ACAAAAANCEVndOCFBc6/PfpTNOAAAAAACaQIsWLVIoFDJ//vyUl5enUCgUuyQgn4Qm8+fPT6FQSIsWLT6zv+AEAAAAAKAJlJSUpHPnznnrrbcya9asYpcD/INCoZDOnTunpKTkM/sKTgAAAAAAmkhpaWl69uyZ5cuXF7sU4B+0aNFinUKTRHACAAAAANCkSkpK1vkXtMDmx+HwAAAAAAAA9QQnAAAAAAAA9QQnAAAAAAAA9QQnAAAAAAAA9QQnAAAAAAAA9QQnAAAAAAAA9QQnAAAAAAAA9baY4GTs2LHp169fttlmm7Rv336V9r/+9a85/vjjU1FRkTZt2mT33XfPdddd95njjhw5Mt27d0+bNm1SXl6eYcOG5eWXX16pT2VlZQqFwkrXD3/4w6aaGgAAAAAAsJloXuwC1tWyZctyzDHH5MADD0xVVdUq7VOnTk15eXluu+22VFRUZMqUKTnjjDNSUlKSs88+e43j9u3bNyeeeGK6dOmSBQsWZPTo0Rk8eHBmzpyZkpKShn6XXXZZTj/99Ib70tLSpp0gAAAAAABQdFtMcHLppZcmSW6++ebVto8YMWKl+1122SVPPvlk/vCHP6w1ODnjjDMa/l1ZWZkf//jH2XvvvTNr1qx07969oW3bbbfNjjvuuAEzAAAAAAAANndbzFZdjVFdXZ2OHTuuc//Fixdn/Pjx6datWyoqKlZqu/LKK7Pddtuld+/eGTt2bJYtW9bU5QIAAAAAAEW2xaw4WV9PPvlk7rzzzvzxj3/8zL7XX399LrjggixevDi9evXKQw89lJYtWza0f//7388+++yTDh065KmnnsqoUaMyc+bM3HjjjWscc+nSpVm6dGnDfU1NzYZNCAAAAAAA2OiKuuJk9OjRqxy6/s/XM888s97jzpgxI8OGDcvFF1+cQYMGfWb/E088Mc8991weffTR9OzZM8cee2yWLFnS0P5v//ZvOeSQQ7LXXnvl29/+dm644YZUVVXlvffeW+OYV1xxRcrKyhquf17BAgAAAAAAbH6KuuLk7LPPznHHHbfWPpWVles15osvvpiBAwfm9NNPz4UXXrhO73wabvTs2TMHHHBAOnTokLvvvjvHH3/8avsfcMABSZK///3v2W677VbbZ9SoUTn33HMb7mtqaoQnAAAAAACwmStqcNKpU6d06tSpycabMWNGBg4cmJNPPjljx45t9Dh1dXUrbbP1z5577rkkyU477bTGPq1atUqrVq0aXQMAAAAAALDpbTFnnMyePTsLFizI7NmzU1tbm2nTpiVJevTokdLS0syYMSMDBgzI4MGDc+6552bevHlJkpKSkpSXlydJ5syZk0MPPTS33HJL9t9//7z++uu54447Mnjw4JSXl2fOnDm58sor06ZNmwwZMiTJJ2el/PnPf86AAQNSVlaWp59+Ov/2b/+WoUOHpkuXLkX5WQAAAAAAABvHFhOcXHzxxZkwYULDfZ8+fZIkkyZNSv/+/XPXXXdl/vz5uf3223P77bc39OvatWtmzZqVJFm+fHleeeWVfPjhh0mS1q1b57HHHsu4cePy/vvvZ4cddshXv/rVTJkyJdtvv32ST1aO3HHHHbn00kuzdOnSdO3aNaeffnouuOCCTTRzAAAAAABgUynU1dXVFbuIrUFNTU3KyspSXV2ddu3aFbsc2KpddNEvUlKyT7HL2OLV1j6bMWO+W+wyAACAIvM7DwA+b5oVuwAAAAAAAIDNheAEAAAAAACgnuAEAAAAAACgnuAEAAAAAACgnuAEAAAAAACgnuAEAAAAAACgnuAEAAAAAACgnuAEAAAAAACgXqODk48//jgPP/xw/tf/+l9ZuHBhkuTtt9/OokWLmqw4AAAAAACATal5Y1564403csQRR2T27NlZunRpBg0alG233TY//elPs2TJktxwww1NXScAAAAAAMBG16gVJ9///vez77775v3330+bNm0anh911FGZOHFikxUHAAAAAACwKTVqxcnjjz+eJ554Ii1btlzpedeuXTNnzpwmKQwAAAAAAGBTa9SKkxUrVqS2tnaV52+99Va23XbbDS4KAAAAAACgGBoVnAwaNCjjxo1ruC8UClm0aFEuueSSDBkypKlqAwAAAAAA2KQatVXXz372swwYMCBf/OIXs2TJkpxwwgn529/+lk6dOuU3v/lNU9cIAAAAAACwSTQqONl5550zbdq0/OY3v8mzzz6bFStW5LTTTsuJJ5640mHxAAAAAAAAW5JGBSdJ0qZNm4wYMSIjRoxoynoAAAAAAACKZp2Dk3vvvXedBx06dGijigEAAAAAACimdQ5OjjzyyJXuC4VC6urqVnmWJLW1tRteGQAAAAAAwCbWbF07rlixouH605/+lN69e+e///u/88EHH6S6ujr//d//nX322ScPPPDAxqwXAAAAAABgo2nUGSfnnHNObrjhhhx00EENzw4//PBss802OeOMM/LSSy81WYEAAAAAAACbyjqvOPlHr732WsrKylZ5XlZWllmzZm1oTQAAAAAAAEXRqOBkv/32yznnnJO5c+c2PJs3b17+/d//Pfvvv3+TFQcAAAAAALApNSo4uemmm/LOO++ka9eu6dGjR3r06JEuXbpk7ty5qaqqauoaAQAAAAAANolGnXHSo0ePTJ8+PQ899FBefvnl1NXV5Ytf/GIOO+ywFAqFpq4RAAAAAABgk2hUcJIkhUIhgwcPzuDBg5uyHgAAAAAAgKJpVHBy2WWXrbX94osvblQxAAAAAAAAxdSo4OTuu+9e6X758uWZOXNmmjdvnu7duwtOAAAAAACALVKjgpPnnntulWc1NTU55ZRTctRRR21wUQAAAAAAAMXQrKkGateuXS677LJcdNFFTTUkAAAAAADAJtVkwUmSfPDBB6murm7KIQEAAABYT6+//nqqqqqKXQYAbJEatVXXz3/+85Xu6+rqMnfu3Nx666054ogjmqQwAAAAABrnuuuuS0VFRZ5++unst99+xS4HALYojQpOfvazn61036xZs5SXl+fkk0/OqFGjmqQwAAAAAFa1bNmyDB48OEkyderU9O3bN7NmzUqfPn1y991355RTTsmMGTPStm3b3HfffbnxxhvTo0ePhveXLFmSkSNH5hvf+EauvfbafPDBB3nvvffSvXv3HHjggWnVqlX23XfffP3rX8+2226bvn37Jkl69OiRG2+8MSeffHKuv/76tG3btijzB4CNrVHBycyZM5u6DgAAAADWQcuWLTN58uQkyb777pvJkydn1qxZOe+88xr6jB8/Pnvuuedq37/jjjsyZMiQDB06NEOHDs3kyZNz33335eqrr06SjB49uqHvbrvt1vBZnzr66KNz2223ZeTIkU06LwDYXDTqjJMRI0Zk4cKFqzxfvHhxRowYscFFAQAAALBx3HPPPTnooIMa/f7AgQNz7733NmFFALB5aVRwMmHChHz00UerPP/oo49yyy23bHBRAAAAADTeqaeemv79+6d///5ZtGjRSm2vv/56dtppp3Ua55VXXmkY58ILL0ySbLvttpk/f36T1wwAm4v12qqrpqYmdXV1qaury8KFC9O6deuGttra2tx///3Zfvvtm7xIAAAAANbd2rbqatmyZZo1W7e/pV3dVl0A8Hm3XsFJ+/btUygUUigUsuuuu67SXigUcumllzZZcQAAAAA0rc6dO+edd95p9B+/Llq0KOXl5U1cFQBsPtYrOJk0aVLq6uoycODA/P73v0/Hjh0b2lq2bJmuXbtm5513bvIiAQAAAFh3p556atq2bZskGTduXHr37t3QNmzYsDz88MM54YQTPnOcT7fqSpIdd9wxv/3tbzNx4sR8/etf3xhlA8BmoVBXV1e3vi+98cYb6dKlSwqFwsao6XOppqYmZWVlqa6uTrt27YpdDmzVLrroFykp2afYZWzxamufzZgx3y12GQAAwHpaunRpvv3tb+fWW29t1Psnn3xyfvWrX6W0tDSJ33kA8PmzzitOpk+fnj333DPNmjVLdXV1nn/++TX23WuvvZqkOAAAAACaVqtWrRodmiTJhAkTmrAaANj8rHNw0rt378ybNy/bb799evfunUKhkNUtVikUCqmtrW3SIgEAAAAAADaFdQ5OZs6c2XDw18yZMzdaQQAAAAAAAMXSbF07du3ateFMk65du6712hjGjh2bfv36ZZtttkn79u1Xaf/rX/+a448/PhUVFWnTpk123333XHfddZ857siRI9O9e/e0adMm5eXlGTZsWF5++eVV+v3xj3/Ml7/85bRp0yadOnXK0Ucf3RTTAgAAAAAANiPrvOLkH917772rfV4oFNK6dev06NEj3bp126DC/tmyZctyzDHH5MADD0xVVdUq7VOnTk15eXluu+22VFRUZMqUKTnjjDNSUlKSs88+e43j9u3bNyeeeGK6dOmSBQsWZPTo0Rk8eHBmzpyZkpKSJMnvf//7nH766bn88sszcODA1NXVrfWMFwAAAAAAYMtUqFvdQSWfoVmzZqs94+TTZ4VCIQcddFD+8z//Mx06dGiyYpPk5ptvzjnnnJMPPvjgM/ueddZZeemll/LII4+s8/jTp0/P3nvvnb///e/p3r17Pv7441RWVubSSy/Naaed1ui6a2pqUlZWlurq6rRr167R4wAb7qKLfpGSkn2KXcYWr7b22YwZ891ilwEAABSZ33kA8HnTqBUnDz30UH70ox9l7Nix2X///ZMkTz31VC688MJcdNFFKSsry8iRI3PeeeetdnXIplJdXZ2OHTuuc//Fixdn/Pjx6datWyoqKpIkzz77bObMmZNmzZqlT58+mTdvXnr37p2rr746e+yxxxrHWrp0aZYuXdpwX1NT0/iJAKv1q19NyLx56//f1oIFy1J/ZBMAAAAAwEoaFZx8//vfz69//ev069ev4dmhhx6a1q1b54wzzsiMGTMybty4jBgxoskKXV9PPvlk7rzzzvzxj3/8zL7XX399LrjggixevDi9evXKQw89lJYtWyZJXn/99STJ6NGjc+2116aysjLXXHNNDjnkkLz66qtrDGauuOKKXHrppU03IWAV8+bVNGrliNAEAAAAAFiTdT4c/h+99tprq1162a5du4agoWfPnnn33XfXOs7o0aNTKBTWej3zzDPrXd+MGTMybNiwXHzxxRk0aNBn9j/xxBPz3HPP5dFHH03Pnj1z7LHHZsmSJUmSFStWJEl+9KMf5Rvf+Eb69u2b8ePHp1Ao5K677lrjmKNGjUp1dXXD9eabb673PAAAAAAAgE2rUStO+vbtm/PPPz+33HJLyuv/dHv+/Pm54IILst9++yVJ/va3v6Vz585rHefss8/Occcdt9Y+lZWV61Xbiy++mIEDB+b000/PhRdeuE7vlJWVpaysLD179swBBxyQDh065O67787xxx+fnXbaKUnyxS9+saF/q1atsssuu2T27NlrHLNVq1Zp1arVetUOAAAAAAAUV6OCk6qqqgwbNiydO3dORUVFCoVCZs+enV122SX33HNPkmTRokW56KKL1jpOp06d0qlTp8aUsFozZszIwIEDc/LJJ2fs2LGNHqeurq7hfJK+ffumVatWeeWVV3LQQQclSZYvX55Zs2ala9euTVI3AAAAAACweWhUcLLbbrvlpZdeyoMPPphXX301dXV16dWrVwYNGpRmzT7Z/evII49syjoze/bsLFiwILNnz05tbW2mTZuWJOnRo0dKS0szY8aMDBgwIIMHD865556befPmJUlKSkoaVsXMmTMnhx56aG655Zbsv//+ef3113PHHXdk8ODBKS8vz5w5c3LllVemTZs2GTJkSJJPth8788wzc8kll6SioiJdu3bNVVddlSQ55phjmnSOAAAAAABAcTUqOEmSQqGQI444IkcccURT1rNGF198cSZMmNBw36dPnyTJpEmT0r9//9x1112ZP39+br/99tx+++0N/bp27ZpZs2Yl+WSlyCuvvJIPP/wwSdK6des89thjGTduXN5///3ssMMO+epXv5opU6Zk++23bxjjqquuSvPmzXPSSSflo48+ype//OU88sgj6dChwyaYOQAAAAAAsKkU6urq6hrz4sSJEzNx4sS88847DQeof+qmm25qkuI+T2pqalJWVpbq6uq0a9eu2OXA58JFF/0iJSX7FLuMrVZt7bMZM+a7xS4DAAAoMr/zAODzplErTi699NJcdtll2XfffbPTTjulUCg0dV0AW53a2trcf/9Vmf//PZe6QrPsscdROeCAY4tdFgAAAABsVRoVnNxwww25+eabc9JJJzV1PQBbpZdeejSPPXpVDjtk1xx80sFZtqw29z04MTf8z1tz/An/b8rKdix2iQAAAACwVWhUcLJs2bL069evqWsB2CotW7Y4j/2fq3LtmCNT3qm04Xmvnjtk2gtzcs2vRuaMM+8pYoUAAAAAsPVo1piXvv3tb+d//+//3dS1AGyV7rnnxxn2/3xppdDkU733/EK6VpTl1VefKEJlAAAAALD1adSKkyVLluTXv/51Hn744ey1115p0aLFSu3XXnttkxQHsDX44P2X0/8rX1tj++ABu+a23/02u+76lU1YFQAAAABsnRoVnEyfPj29e/dOkrzwwgsrtTkoHmD9FJLUrVhz+4oVdf7fCgAAAACbSKOCk0mTJjV1HQBbrY6d9srE//Nqjjt6n9W2P/DIKzmw3/c2cVUAAAAAsHVq1Bkn/+itt97KnDlzmqIWgK3S17/+w/zxoRfz9rzqVdr+MvWNzHl7Ubp3378IlQEAAADA1qdRwcmKFSty2WWXpaysLF27dk2XLl3Svn37jBkzJitWrGW/GQBW0bJlmww6fHR+cOl/5VdVj+fFV+bmueffyqVXPZBf3DglJ/zLTcUuEQAAAAC2Go3aqutHP/pRqqqq8pOf/CRf+cpXUldXlyeeeCKjR4/OkiVLMnbs2KauE+BzrVu3fXPa6fdk0qRf5+rrpyRpnn33OyEDBw8udmkAAAAAsFVpVHAyYcKE3HjjjRk6dGjDs7333jtf+MIX8q//+q+CE4BGKCkpyWGHfSfJd4pdCgAAAABstRq1VdeCBQvSq1evVZ736tUrCxYs2OCiAAAAAAAAiqFRwcnee++dX/7yl6s8/+Uvf5m99tprg4sCAAAAAAAohkZt1fXTn/40X/va1/Lwww/nwAMPTKFQyJQpU/Lmm2/m/vvvb+oaAQAAAAAANolGrTg55JBD8uqrr+aoo47KBx98kAULFuToo4/OjBkzMn78+KauEQAAAAAAYJNo1IqTJNl5551XOQT+r3/9ayZMmJCbbrppgwsDAAAAAADY1Bq14gQAAAAAAODzSHACAAAAAABQT3ACAAAAAABQb73OODn66KPX2v7BBx9sSC0AAAAAAABFtV7BSVlZ2We2Dx8+fIMKAgAAAAAAKJb1Ck7Gjx+/seoAoAnMn//ndOzYcpN8Vnn5NpvkcwAAAABgU1qv4ASAzVvHji0zZsx3i10GAAAAAGyxHA4PAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQT3ACAAAAAABQr3mxCwBIkl/9akLmzatZr3cWLFiW8vKNVFCRzZ//53Ts2HK939txx3YboRoAAAAA2HoIToDNwrx5NSkp2We93vm8hiZJ0rFjy4wZ891ilwEAAAAAWx1bdQEAAAAAANQTnAAAAAAAANQTnAAAAAAAANQTnAAAAAAAANQTnAAAAAAAANQTnAAAAAAAANTbYoKTsWPHpl+/ftlmm23Svn37Vdr/+te/5vjjj09FRUXatGmT3XffPdddd91njjty5Mh07949bdq0SXl5eYYNG5aXX365oX3y5MkpFAqrvZ5++ummnCIAAAAAAFBkW0xwsmzZshxzzDH5zne+s9r2qVOnpry8PLfddltmzJiRH/3oRxk1alR++ctfrnXcvn37Zvz48XnppZfy4IMPpq6uLoMHD05tbW2SpF+/fpk7d+5K17e//e1UVlZm3333bfJ5AgAAAAAAxdO82AWsq0svvTRJcvPNN6+2fcSIESvd77LLLnnyySfzhz/8IWefffYaxz3jjDMa/l1ZWZkf//jH2XvvvTNr1qx07949LVu2zI477tjQZ/ny5bn33ntz9tlnp1AobMCMAAAAAACAzc0WE5w0RnV1dTp27LjO/RcvXpzx48enW7duqaioWG2fe++9N++++25OOeWUtY61dOnSLF26tOG+pqZmnesAAAAAAACKY4vZqmt9Pfnkk7nzzjszcuTIz+x7/fXXp7S0NKWlpXnggQfy0EMPpWXLlqvtW1VVlcMPP3yNwcqnrrjiipSVlTVcn9UfAAAAAAAovqIGJ6NHj17jweufXs8888x6jztjxowMGzYsF198cQYNGvSZ/U888cQ899xzefTRR9OzZ88ce+yxWbJkySr93nrrrTz44IM57bTTPnPMUaNGpbq6uuF6880313seAAAAAADAplXUrbrOPvvsHHfccWvtU1lZuV5jvvjiixk4cGBOP/30XHjhhev0zqerQnr27JkDDjggHTp0yN13353jjz9+pX7jx4/Pdtttl6FDh37mmK1atUqrVq3Wq3YAAAAAAKC4ihqcdOrUKZ06dWqy8WbMmJGBAwfm5JNPztixYxs9Tl1d3Urnk3z6bPz48Rk+fHhatGixoaUCAAAAAACboS3mjJPZs2dn2rRpmT17dmprazNt2rRMmzYtixYtSvJJaDJgwIAMGjQo5557bubNm5d58+Zl/vz5DWPMmTMnvXr1ylNPPZUkef3113PFFVdk6tSpmT17dp588skce+yxadOmTYYMGbLS5z/yyCOZOXPmOm3TBQAAAAAAbJmKuuJkfVx88cWZMGFCw32fPn2SJJMmTUr//v1z1113Zf78+bn99ttz++23N/Tr2rVrZs2alSRZvnx5XnnllXz44YdJktatW+exxx7LuHHj8v7772eHHXbIV7/61UyZMiXbb7/9Sp9fVVWVfv36Zffdd9/IMwUAAAAAAIqlUFdXV1fsIrYGNTU1KSsrS3V1ddq1a1fscmCzc9FFv0hJyT7FLmOzUVv7bMaM+W6xywAAAPhMfucBwOfNFrNVFwAAAAAAwMYmOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKgnOAEAAAAAAKjXvNgFAHzezZ//53Ts2HK93tlxx3YbqRoAAAAAYG0EJwAbWceOLTNmzHeLXQYAAAAAsA5s1QUAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFBPcAIAAAAAAFCvebELAEiSHXdsl3nzni12GRvFjju2K3YJAAAAAMA6EpwAm4Wzzjq52CUAAAAAANiqCwAAAAAA4FOCEwAAAAAAgHqCEwAAAAAAgHqCEwAAAAAAgHqCEwAAAAAAgHqCEwAAAAAAgHrNi13A1qKuri5JUlNTU+RKAAAAAJrOp7/r+PR3HwCwpROcbCILFy5MklRUVBS5EgAAAICmt3DhwpSVlRW7DADYYIU6fw6wSaxYsSJvv/12tt122xQKhWKXw1aipqYmFRUVefPNN9OuXbtilwMAbMZ8bwCgserq6rJw4cLsvPPOadbMrvAAbPmsONlEmjVrls6dOxe7DLZS7dq18wsQAGCd+N4AQGNYaQLA54k/AwAAAAAAAKgnOAEAAAAAAKgnOIHPsVatWuWSSy5Jq1atil0KALCZ870BAADgEw6HBwAAAAAAqGfFCQAAAAAAQD3BCQAAAAAAQD3BCQAAAAAAQD3BCQAAAAAAQD3BCWzBhg4dmi5duqR169bZaaedctJJJ+Xtt99eqc/TTz+dQw89NO3bt0+HDh0yePDgTJs2ba3jjhw5Mt27d0+bNm1SXl6eYcOG5eWXX96IMwEANraN8b1hwYIF+e53v5vddtst22yzTbp06ZLvfe97qa6u3sizAQAA2HgEJ7AFGzBgQO6888688sor+f3vf5/XXnst3/zmNxvaFy5cmMMPPzxdunTJX/7ylzz++ONp165dDj/88CxfvnyN4/bt2zfjx4/PSy+9lAcffDB1dXUZPHhwamtrN8W0AICNYGN8b3j77bfz9ttv5+qrr87zzz+fm2++OQ888EBOO+20TTUtAACAJleoq6urK3YRQNO49957c+SRR2bp0qVp0aJFnnnmmey3336ZPXt2KioqkiTPP/989tprr/z9739P9+7d12nc6dOnZ++9916vdwCAzdvG+t5w11135V/+5V+yePHiNG/efGNOAQAAYKOw4gQ+JxYsWJDbb789/fr1S4sWLZIku+22Wzp16pSqqqosW7YsH330UaqqqrLHHnuka9eu6zTu4sWLM378+HTr1q3hlygAwJZtY31vSJLq6uq0a9dOaAIAAGyxBCewhfvBD36Qtm3bZrvttsvs2bNzzz33NLRtu+22mTx5cm677ba0adMmpaWlefDBB3P//fd/5i8zrr/++pSWlqa0tDQPPPBAHnroobRs2XJjTwcA2Ig21veGT7333nsZM2ZMRo4cubGmAAAAsNEJTmAzM3r06BQKhbVezzzzTEP/888/P88991z+9Kc/paSkJMOHD8+nO/B99NFHGTFiRL7yla/kz3/+c5544onsscceGTJkSD766KO11nHiiSfmueeey6OPPpqePXvm2GOPzZIlSzbq3AGA9bO5fG9Ikpqamnzta1/LF7/4xVxyySUbbc4AAAAbmzNOYDPz7rvv5t13311rn8rKyrRu3XqV52+99VYqKioyZcqUHHjggamqqsp//Md/ZO7cuWnW7JOcdNmyZenQoUOqqqpy3HHHrVNNn75z44035vjjj1//SQEAG8Xm8r3h04Plt9lmm9x3332r/TwAAIAthY2HYTPTqVOndOrUqVHvfpqDLl26NEny4YcfplmzZikUCg19Pr1fsWLFeo/96bgAwOZhc/jeUFNTk8MPPzytWrXKvffeKzQBAAC2eLbqgi3UU089lV/+8peZNm1a3njjjUyaNCknnHBCunfvngMPPDBJMmjQoLz//vs566yz8tJLL2XGjBk59dRT07x58wwYMCBJMmfOnPTq1StPPfVUkuT111/PFVdckalTp2b27Nl58sknc+yxx6ZNmzYZMmRI0eYLADTexvresHDhwgwePDiLFy9OVVVVampqMm/evMybNy+1tbVFmy8AAMCGEJzAFqpNmzb5wx/+kEMPPTS77bZbRowYkT333DOPPvpoWrVqlSTp1atX/uu//ivTp0/PgQcemIMPPjhvv/12Hnjggey0005JkuXLl+eVV17Jhx9+mCRp3bp1HnvssQwZMiQ9evTIsccem7Zt22bKlCnZfvvtizZfAKDxNtb3hqlTp+Yvf/lLnn/++fTo0SM77bRTw/Xmm28Wbb4AAAAbwhknAAAAAAAA9aw4AQAAAAAAqCc4AQAAAAAAqCc4AQAAAAAAqCc4AQAAAAAAqCc4AQAAAAAAqCc4AQAAAAAAqCc4AQAAAAAAqCc4AQDYglVWVmbcuHEbZez+/fvnnHPO2ShjAwAAwOZKcAIAsImccsopOfLIIxv17s0335z27duv8vzpp5/OGWec0XBfKBTyn//5n40rEAAAAEjzYhcAAEDjlZeXF7sEAAAA+Fyx4gQAYDNw7bXX5ktf+lLatm2bioqK/Ou//msWLVqUJJk8eXJOPfXUVFdXp1AopFAoZPTo0UlW3qqrsrIySXLUUUelUCg03K9upcs555yT/v37N9wvXrw4w4cPT2lpaXbaaadcc801q9S4bNmyXHDBBfnCF76Qtm3b5stf/nImT57chD8FAAAAKD7BCQDAZqBZs2b5+c9/nhdeeCETJkzII488kgsuuCBJ0q9fv4wbNy7t2rXL3LlzM3fu3Jx33nmrjPH0008nScaPH5+5c+c23K+L888/P5MmTcrdd9+dP/3pT5k8eXKmTp26Up9TTz01TzzxRH77299m+vTpOeaYY3LEEUfkb3/72wbMHAAAADYvtuoCANgM/OMh7N26dcuYMWPyne98J9dff31atmyZsrKyFAqF7Ljjjmsc49Ntu9q3b7/Wfv9s0aJFqaqqyi233JJBgwYlSSZMmJDOnTs39Hnttdfym9/8Jm+99VZ23nnnJMl5552XBx54IOPHj8/ll1++PtMFAACAzZbgBABgMzBp0qRcfvnlefHFF1NTU5OPP/44S5YsyeLFi9O2bduN+tmvvfZali1blgMPPLDhWceOHbPbbrs13D/77LOpq6vLrrvuutK7S5cuzXbbbbdR6wMAAIBNSXACAFBkb7zxRoYMGZIzzzwzY8aMSceOHfP444/ntNNOy/Llyzd4/GbNmqWurm6lZ/847j+3rc6KFStSUlKSqVOnpqSkZKW20tLSDa4RAAAANheCEwCAInvmmWfy8ccf55prrkmzZp8cQXfnnXeu1Kdly5apra39zLFatGixSr/y8vK88MILKz2bNm1aWrRokSTp0aNHWrRokT//+c/p0qVLkuT999/Pq6++mkMOOSRJ0qdPn9TW1uadd97JwQcf3LiJAgAAwBbA4fAAAJtQdXV1pk2bttJVXl6ejz/+OL/4xS/y+uuv59Zbb80NN9yw0nuVlZVZtGhRJk6cmHfffTcffvjhasevrKzMxIkTM2/evLz//vtJkoEDB+aZZ57JLbfckr/97W+55JJLVgpSSktLc9ppp+X888/PxIkT88ILL+SUU05pCHGSZNddd82JJ56Y4cOH5w9/+ENmzpyZp59+OldeeWXuv//+jfCTAgAAgOIQnAAAbEKTJ09Onz59VrpuuummXHvttbnyyiuz55575vbbb88VV1yx0nv9+vXLmWeemW9961spLy/PT3/609WOf8011+Shhx5KRUVF+vTpkyQ5/PDDc9FFF+WCCy7Ifvvtl4ULF2b48OErvXfVVVflq1/9aoYOHZrDDjssBx10UPr27btSn/Hjx2f48OH593//9+y2224ZOnRo/vKXv6SioqIJf0IAAABQXIW6ddnUGgAAAAAAYCtgxQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEA9wQkAAAAAAEC9/x/tQNfa71OOXwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f, a = plt.subplots(2,1,figsize=(14,6),sharex=True)\n", "\n", "# --- 2022\n", "\n", "geo_cp_m_2022.plot(ax=a[0],edgecolor='k',lw=.5,facecolor='lightblue',alpha=.5)\n", "p2022_m_1 = p2022_m[p2022_m['CAMPO'].isin(cp_list_2022)]\n", "p2022_m_1f = sb.scatterplot(data=p2022_m_1,\n", " x='LONGITUDE1',\n", " y='LATITUDE_1',\n", " hue='POCO',\n", " marker='o',edgecolor='k',\n", " palette='magma',\n", " ax=a[0])\n", "p2022_m_1f.legend(bbox_to_anchor=(1.4,1),ncol=2,fontsize=8)\n", "\n", "a[0].set_ylabel('Longitude')\n", "a[0].set_title(f'Poços públicos em 2022 - Bacia: {basin}',fontsize=10)\n", "\n", "geo_cp_m_2022.apply(lambda x: a[0].annotate(\n", " text=f'{x.NOM_CAMPO} ({x.SIG_CAMPO})',\n", " xy=np.array(x.geometry.centroid.coords[0]) + np.array([0.15,0]), \n", " ha='left', \n", " va='center',\n", " color = \"black\",\n", " fontsize=6, \n", " weight='light'), \n", " axis=1)\n", "\n", "\n", "# --- 2023\n", "\n", "geo_cp_m_2023.plot(ax=a[1],edgecolor='k',lw=.5,facecolor='darkblue',alpha=.5)\n", "p2023_m_1 = p2023_m[p2023_m['CAMPO'].isin(cp_list_2023)]\n", "p2023_m_1f = sb.scatterplot(data=p2023_m_1,\n", " x='LONGITUDE_BASE_DD',\n", " y='LATITUDE_BASE_DD',\n", " hue='POCO',\n", " marker='o',edgecolor='k',\n", " palette='cividis',\n", " ax=a[1])\n", "p2023_m_1f.legend(bbox_to_anchor=(1.4,1),ncol=2,fontsize=8)\n", "\n", "a[1].set_xlabel('Latitude')\n", "a[1].set_ylabel('Longitude')\n", "a[1].set_title(f'Poços públicos em 2023 - Bacia: {basin}',fontsize=10)\n", "\n", "geo_cp_m_2023.apply(lambda x: a[1].annotate(\n", " text=f'{x.NOM_CAMPO} ({x.SIG_CAMPO})',\n", " xy=np.array(x.geometry.centroid.coords[0]) + np.array([0.15,0]), \n", " ha='left', \n", " va='center',\n", " color = \"black\",\n", " fontsize=6, \n", " weight='light'), \n", " axis=1);\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "dataviz", "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.13.0" } }, "nbformat": 4, "nbformat_minor": 2 }