{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Impact Data functionalities\n", "\n", "\n", "Import data from EM-DAT CSV file and populate Impact()-object with the data.\n", "\n", "\n", "The core functionality of the module is to read disaster impact data as downloaded from the International Disaster Database EM-DAT (www.emdat.be) and produce a CLIMADA Impact()-instance from it.\n", "The purpose is to make impact data easily available for comparison with simulated impact inside CLIMADA, e.g. for calibration purposes.\n", "\n", "\n", "## Data Source\n", "The International Disaster Database EM-DAT www.emdat.be\n", "\n", "Download: https://public.emdat.be/ (register for free and download data to continue)\n", "\n", "\n", "## Most important functions\n", "- clean_emdat_df: read CSV from EM-DAT into a DataFrame and clean up.\n", "- emdat_to_impact: create Impact-instance populated with impact data from EM-DAT data (CSV).\n", "- emdat_countries_by_hazard: get list of countries affected by a certain haazrd (disaster (sub-)type) in EM-DAT.\n", "- emdat_impact_yearlysum: create DataFrame with impact from EM-DAT summed per country and year.\n", "\n", "\n", "\n", "## Demo data\n", "\n", "The demo data used here (demo_emdat_impact_data_2020.csv) contains entries for the disaster subtype \"Tropical cyclone\" from 2000 to 2020.\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "\"\"\"Load required packages and set path to CSV-file from EM-DAT\"\"\"\n", "\n", "import numpy as np\n", "import pandas as pd\n", "from matplotlib import pyplot as plt\n", "\n", "from climada.util.constants import DEMO_DIR\n", "from climada.engine.impact_data import (\n", " emdat_countries_by_hazard,\n", " emdat_impact_yearlysum,\n", " emdat_to_impact,\n", " clean_emdat_df,\n", ")\n", "\n", "# set path to CSV file downloaded from https://public.emdat.be :\n", "emdat_file_path = DEMO_DIR.joinpath(\"demo_emdat_impact_data_2020.csv\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### clean_emdat_df()\n", "\n", "read CSV from EM-DAT into a DataFrame and clean up.\n", "\n", "Use the parameters countries, hazard, and year_range to filter. These parameters are the same for most functions shown here." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " Dis No Year Seq Disaster Group Disaster Subgroup Disaster Type \\\n", "0 2005-0540-VNM 2005 540 Natural Meteorological Storm \n", "1 2005-0540-THA 2005 540 Natural Meteorological Storm \n", "2 2005-0536-VNM 2005 536 Natural Meteorological Storm \n", "3 2005-0611-VNM 2005 611 Natural Meteorological Storm \n", "4 2006-0362-VNM 2006 362 Natural Meteorological Storm \n", "5 2006-0648-VNM 2006 648 Natural Meteorological Storm \n", "6 2006-0251-VNM 2006 251 Natural Meteorological Storm \n", "7 2006-0517-VNM 2006 517 Natural Meteorological Storm \n", "\n", " Disaster Subtype Disaster Subsubtype Event Name Entry Criteria \\\n", "0 Tropical cyclone NaN Damrey Kill \n", "1 Tropical cyclone NaN Damrey Kill \n", "2 Tropical cyclone NaN Vicente Kill \n", "3 Tropical cyclone NaN Kai Tak (21) Kill \n", "4 Tropical cyclone NaN Bilis Kill \n", "5 Tropical cyclone NaN Durian (Reming) Kill \n", "6 Tropical cyclone NaN Chanchu (Caloy) Kill \n", "7 Tropical cyclone NaN Xangsane (Milenyo) Kill \n", "\n", " ... End Day Total Deaths No Injured No Affected No Homeless Total Affected \\\n", "0 ... 30.0 75.0 28.0 337632.0 NaN 337660.0 \n", "1 ... 30.0 10.0 NaN 2000.0 NaN 2000.0 \n", "2 ... 19.0 8.0 NaN 8500.0 NaN 8500.0 \n", "3 ... 4.0 20.0 NaN 15000.0 NaN 15000.0 \n", "4 ... 19.0 17.0 NaN NaN 2000.0 2000.0 \n", "5 ... 8.0 95.0 1360.0 975000.0 250000.0 1226360.0 \n", "6 ... 17.0 204.0 NaN 600000.0 NaN 600000.0 \n", "7 ... 6.0 71.0 525.0 1368720.0 98680.0 1467925.0 \n", "\n", " Reconstruction Costs ('000 US$) Insured Damages ('000 US$) \\\n", "0 NaN NaN \n", "1 NaN NaN \n", "2 NaN NaN \n", "3 NaN NaN \n", "4 NaN NaN \n", "5 NaN NaN \n", "6 NaN NaN \n", "7 NaN NaN \n", "\n", " Total Damages ('000 US$) CPI \n", "0 219250.0 76.388027 \n", "1 20000.0 76.388027 \n", "2 20000.0 76.388027 \n", "3 11000.0 76.388027 \n", "4 NaN 78.852256 \n", "5 456000.0 78.852256 \n", "6 NaN 78.852256 \n", "7 624000.0 78.852256 \n", "\n", "[8 rows x 43 columns]\n" ] } ], "source": [ "\"\"\"Create DataFrame df with EM-DAT entries of tropical cyclones in Thailand and Viet Nam in the years 2005 and 2006\"\"\"\n", "\n", "df = clean_emdat_df(\n", " emdat_file_path,\n", " countries=[\"THA\", \"Viet Nam\"],\n", " hazard=[\"TC\"],\n", " year_range=[2005, 2006],\n", ")\n", "print(df)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### emdat_countries_by_hazard()\n", "\n", "Pick a hazard and a year range to get a list of countries affected from the EM-DAT data." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['China', 'Dominican Republic', 'Antigua and Barbuda', 'Fiji', 'Australia', 'Bangladesh', 'Belize', 'Barbados', 'Cook Islands', 'Canada', 'Bahamas', 'Guatemala', 'Jamaica', 'Saint Lucia', 'Madagascar', 'Mexico', \"Korea, Democratic People's Republic of\", 'El Salvador', 'Myanmar', 'French Polynesia', 'Solomon Islands', 'Taiwan, Province of China', 'India', 'United States of America', 'Honduras', 'Haiti', 'Pakistan', 'Philippines', 'Hong Kong', 'Korea, Republic of', 'Nicaragua', 'Oman', 'Japan', 'Puerto Rico', 'Thailand', 'Martinique', 'Papua New Guinea', 'Tonga', 'Venezuela, Bolivarian Republic of', 'Viet Nam', 'Saint Vincent and the Grenadines', 'Vanuatu', 'Dominica', 'Cuba', 'Comoros', 'Mozambique', 'Malawi', 'Samoa', 'South Africa', 'Sri Lanka', 'Palau', 'Wallis and Futuna', 'Somalia', 'Seychelles', 'Réunion', 'Kiribati', 'Cabo Verde', 'Micronesia, Federated States of', 'Panama', 'Costa Rica', 'Yemen', 'Tuvalu', 'Northern Mariana Islands', 'Colombia', 'Anguilla', 'Djibouti', 'Cambodia', 'Macao', 'Indonesia', 'Guadeloupe', 'Turks and Caicos Islands', 'Saint Kitts and Nevis', \"Lao People's Democratic Republic\", 'Mauritius', 'Marshall Islands', 'Portugal', 'Virgin Islands, U.S.', 'Zimbabwe', 'Saint Barthélemy', 'Virgin Islands, British', 'Saint Martin (French part)', 'Sint Maarten (Dutch part)', 'Tanzania, United Republic of']\n", "['CHN', 'DOM', 'ATG', 'FJI', 'AUS', 'BGD', 'BLZ', 'BRB', 'COK', 'CAN', 'BHS', 'GTM', 'JAM', 'LCA', 'MDG', 'MEX', 'PRK', 'SLV', 'MMR', 'PYF', 'SLB', 'TWN', 'IND', 'USA', 'HND', 'HTI', 'PAK', 'PHL', 'HKG', 'KOR', 'NIC', 'OMN', 'JPN', 'PRI', 'THA', 'MTQ', 'PNG', 'TON', 'VEN', 'VNM', 'VCT', 'VUT', 'DMA', 'CUB', 'COM', 'MOZ', 'MWI', 'WSM', 'ZAF', 'LKA', 'PLW', 'WLF', 'SOM', 'SYC', 'REU', 'KIR', 'CPV', 'FSM', 'PAN', 'CRI', 'YEM', 'TUV', 'MNP', 'COL', 'AIA', 'DJI', 'KHM', 'MAC', 'IDN', 'GLP', 'TCA', 'KNA', 'LAO', 'MUS', 'MHL', 'PRT', 'VIR', 'ZWE', 'BLM', 'VGB', 'MAF', 'SXM', 'TZA']\n" ] } ], "source": [ "\"\"\"emdat_countries_by_hazard: get lists of countries impacted by tropical cyclones from 2010 to 2019\"\"\"\n", "\n", "iso3_codes, country_names = emdat_countries_by_hazard(\n", " emdat_file_path, hazard=\"TC\", year_range=(2010, 2019)\n", ")\n", "\n", "print(country_names)\n", "\n", "print(iso3_codes)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### emdat_to_impact()\n", "function to load EM-DAT impact data and return impact set with impact per event\n", "\n", "#### Parameters:\n", "- emdat_file_csv (str): Full path to EMDAT-file (CSV)\n", "- hazard_type_climada (str): Hazard type abbreviation used in CLIMADA, e.g. 'TC'\n", "\n", "#### Optional parameters:\n", "\n", "- hazard_type_emdat (list or str): List of Disaster (sub-)type according EMDAT terminology or CLIMADA hazard type abbreviations. e.g. ['Wildfire', 'Forest fire'] or ['BF']\n", "- year_range (list with 2 integers): start and end year e.g. [1980, 2017]\n", "- countries (list of str): country ISO3-codes or names, e.g. ['JAM', 'CUB']. Set to None or ['all'] for all countries \n", "- reference_year (int): reference year of exposures for normalization. Impact is scaled proportional to GDP to the value of the reference year. No scaling for reference_year=0 (default)\n", "- imp_str (str): Column name of impact metric in EMDAT CSV, e.g. 'Total Affected'; default = \"Total Damages\"\n", "\n", "#### Returns:\n", "- impact_instance (instance of climada.engine.Impact):\n", " Impact() instance (same format as output from CLIMADA impact computations).\n", " Values are scaled with GDP to reference_year if reference_year not equal 0.\n", " impact_instance.eai_exp holds expected annual impact for each country.\n", " impact_instance.coord_exp holds rough central coordinates for each country.\n", "- countries (list): ISO3-codes of countries imn same order as in impact_instance.eai_exp\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2021-10-19 16:44:58,210 - climada.engine.impact_data - WARNING - ISO3alpha code not found in iso_country: SPI\n", "2021-10-19 16:44:59,007 - climada.engine.impact_data - WARNING - Country not found in iso_country: SPI\n", "Number of TC events in EM-DAT 2000 to 2009 globally: 533\n", "Global annual average monetary damage (AAI) from TCs as reported in EM-DAT 2000 to 2009: USD billion 38.07\n" ] } ], "source": [ "\"\"\"Global TC damages 2000 to 2009\"\"\"\n", "\n", "impact_emdat, countries = emdat_to_impact(\n", " emdat_file_path, \"TC\", year_range=(2000, 2009)\n", ")\n", "\n", "print(\n", " \"Number of TC events in EM-DAT 2000 to 2009 globally: %i\"\n", " % (impact_emdat.event_id.size)\n", ")\n", "print(\n", " \"Global annual average monetary damage (AAI) from TCs as reported in EM-DAT 2000 to 2009: USD billion %2.2f\"\n", " % (impact_emdat.aai_agg / 1e9)\n", ")" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2025-11-26 10:27:26,085 - climada.engine.impact - WARNING - The Impact.tot_value attribute is deprecated.Use Exposures.affected_total_value to calculate the affected total exposure value based on a specific hazard intensity threshold\n", "Number of TC events in EM-DAT in the Philipppines, 2013: 8\n", "\n", "People affected by TC events in the Philippines in 2013 (per event):\n", "[7.269600e+04 1.059700e+04 8.717550e+05 2.204430e+05 1.610687e+07\n", " 3.596000e+03 3.957300e+05 2.628840e+05]\n", "\n", "People affected by TC events in the Philippines in 2013 (total):\n", "17944571\n", "2025-11-26 10:27:26,187 - climada.engine.impact - WARNING - The Impact.tot_value attribute is deprecated.Use Exposures.affected_total_value to calculate the affected total exposure value based on a specific hazard intensity threshold\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHJCAYAAACMppPqAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAASHNJREFUeJzt3XlcVGX///H3oLIYi+KCYgguuSCmgli4pOQWJpbpN9Nyy7pTqTSycvmaud95W1mKlllamablkpVZaJqWlpTaN/clDFPIHZTCBc7vj37M7QgoAzMsh9fz8ZjHw7nmmnM+MwecN9d1rjMWwzAMAQAAmJBLcRcAAADgLAQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQd2M1iseTrtmnTJoftc9GiRbJYLPrpp58cts2i0qFDB3Xo0KG4y3CIuXPnatGiRQ7b3tGjR2WxWBy6zSVLlmjWrFl57mvmzJkO21derv9d8PHxUYcOHfTFF1/k6Pfkk0/edHvZP/9Hjx61tg0aNEhBQUE2/YKCgjRo0CAHvIIbK6r9OEpycrL+93//VxEREapataq8vb0VFham+fPnKzMzM0f/ixcvauTIkfL395e7u7uaN2+ujz76KEe/7777To899pjCwsLk5uaW4xhlS09P10MPPaSGDRvKy8tLt9xyi5o0aaIpU6YoPT3dGS8Z1yhf3AWg9Nm2bZvN/cmTJ2vjxo365ptvbNqDg4OLsqwSa+7cucVdgsPMnTtXVatWddiHXM2aNbVt2zbVq1fPIduT/gk6u3fv1siRIx22zYLo3bu3nn32WWVlZem3337TlClTFB0drc8++0z33nuvXdu69957tW3bNtWsWfOG/VatWiVvb+/ClJ0vRbUfR/n555/1/vvva8CAARo/frwqVKigL7/8UsOGDdMPP/ygd99916b/Aw88oISEBP373/9WgwYNtGTJEvXt21dZWVnq16+ftd+GDRu0fv16tWjRQt7e3nn+cXflyhUZhqHY2FjVqVNHLi4u2rx5syZNmqRNmzZp/fr1znz5MIBCGjhwoHHLLbc4dR8LFy40JBkJCQlO3Q9urEmTJkb79u2Lu4wbuvfee43AwMAc7YmJiYYk4z//+Y/Ta5BkxMTE2LQdPnzYkGR06tTphv3ya+DAgbm+TuR09uxZ4/LlyznaY2JiDElGUlKSte2LL74wJBlLliyx6du5c2fD39/fuHr1qrUtMzPT+u///Oc/hiQjMTEx33U9//zzhiTjyJEjdrwa2IupKzjckCFD5Ovrq7/++ivHY3fffbeaNGlivZ89dP/WW2+pQYMGcnNzU3BwcK7DxJJ04cIFDRs2TFWrVlWVKlX0wAMP6MSJEzZ9srKyNGPGDDVq1Ehubm6qXr26BgwYoD/++CPH9t599101a9ZM7u7u8vX1Vc+ePbVv3z6bPoMGDZKnp6cOHz6sbt26ydPTUwEBAXr22Wd16dKlm74f109dZU+h/Oc//9HLL7+soKAgeXh4qEOHDjp48KCuXLmi0aNHy9/fXz4+PurZs6dOnjxps82goCB1795dq1at0u233y53d3fVrVtXb7zxhk2/jIwMPfvss2revLl8fHzk6+uriIgIffrppznqzMrK0uzZs9W8eXN5eHioUqVKuvPOO7VmzRrrPvfs2aNvv/3WOiWTPXWSlZWlKVOmqGHDhtbn3n777Xr99ddv+N7kNnX10ksvyWKxaM+ePerbt698fHzk5+enRx99VKmpqTd9r7/44gv9/vvvNlNH13v11VdVp04deXp6KiIiQj/88EOOPj/99JN69OghX19fubu7q0WLFlq+fPkN938j9erVU7Vq1fT777/neOyDDz5Q48aNVbFiRTVr1kyff/65zeO5TV3l5voppU2bNslisWjx4sWKjY1VjRo15OHhofbt22vnzp02z83+Od+zZ486duyoW265RdWqVdOTTz6Z43c5r/0sXbpU48aNk7+/v7y9vdWpUycdOHAgR53r169Xx44d5e3trYoVK6pNmzbasGGDTZ9Tp07pX//6lwICAuTm5qZq1aqpTZs2BRr9qFy5sipUqJCjvVWrVpJk83/DqlWr5Onpqf/5n/+x6Tt48GCdOHFCP/74o7XNxaVwH6HVqlWTJJUvz+SKMxF04HAjRozQuXPntGTJEpv2vXv3auPGjYqJibFpX7Nmjd544w1NmjRJn3zyiQIDA9W3b1998sknObb92GOPqUKFClqyZIlmzJihTZs26ZFHHrHpM2zYML3wwgvq3Lmz1qxZo8mTJ2vdunVq3bq1Tp8+be03ffp0DRkyRE2aNNHKlSv1+uuv6//+7/8UERGhQ4cO2WzzypUr6tGjhzp27KhPP/1Ujz76qF577TW9/PLLBX6f4uLi9P333ysuLk4LFizQ/v37FR0drSFDhujUqVN69913NWPGDK1fv16PPfZYjufv2rVLI0eO1DPPPKNVq1apdevWGjFihM05KJcuXdLZs2c1atQorV69WkuXLlXbtm31wAMP6P3337fZ3qBBgzRixAiFh4dr2bJl+uijj9SjRw/rh+uqVatUt25dtWjRQtu2bdO2bdu0atUqSdKMGTP00ksvqW/fvvriiy+0bNkyDRkyROfPny/w+9OrVy81aNBAK1as0OjRo7VkyRI988wzN3zO3Llz1aZNG9WoUcNa4/VTrXFxcYqPj9esWbP04YcfKj09Xd26dbMJURs3blSbNm10/vx5vfnmm/r000/VvHlz9enTp8DnE507d05nzpyxfrhl++KLLzRnzhxNmjRJK1assAbu3377rUD7yc3YsWP122+/acGCBVqwYIFOnDihDh065NjHlStX1K1bN3Xs2FGrV6+2/hHSp0+ffO/n999/14IFCzR//nwdOnRI0dHRNufBLF68WF26dJG3t7fee+89LV++XL6+vuratatN2Onfv79Wr16tF198UV9//bUWLFigTp066cyZM455UyR98803Kl++vBo0aGBt2717txo3bpwjfNx+++3WxwvKMAxdvXpVaWlpWrdunV555RX17dtXtWvXLvA2kQ/FPaSE0i+3qav27dsbzZs3t2kbNmyY4e3tbVy4cMHaJsnw8PAwUlJSrG1Xr141GjVqZNSvX9/alj11NXz4cJttzpgxw5BkJCcnG4ZhGPv27cu1348//mhIMsaOHWsYhmGcO3fO8PDwMLp162bTLykpyXBzczP69etn8/okGcuXL7fp261bN6Nhw4Y3fnP+/3tx7XRP9hRKs2bNbIa+Z82aZUgyevToYfP8kSNHGpKM1NRUa1tgYKBhsViMXbt22fTt3Lmz4e3tbaSnp+day9WrV40rV64YQ4YMMVq0aGFt37x5syHJGDdu3A1fS15TV927d89xvPMj+71YuHChtW3ChAmGJGPGjBk2fYcPH264u7sbWVlZN9zmzaaumjZtajP9sH37dkOSsXTpUmtbo0aNjBYtWhhXrlyx2Ub37t2NmjVr2hy33GT/DF65csW4fPmysW/fPiMqKsqQZMTFxdn08/PzM9LS0qxtKSkphouLizF9+nRrW/bP/7XTIrlNXQUGBhoDBw603t+4caMhyQgNDbV5344ePWpUqFDBeOyxx2y2J8l4/fXXbbY5depUQ5Lx3Xff3XQ/1/8+LV++3JBkbNu2zTAMw0hPTzd8fX2N6Ohom36ZmZlGs2bNjFatWlnbPD09jZEjRxrO8tVXXxkuLi7GM888Y9N+2223GV27ds3R/8SJE4YkY9q0abluLz9TV0uXLjUkWW+DBw/O8TMGx2NEB04xYsQI7dq1S99//70kKS0tTR988IEGDhwoT09Pm74dO3aUn5+f9X65cuXUp08fHT58OMd0U48ePWzuZ/+VlT0dsHHjRknKcbJsq1at1LhxY+tfjNu2bdPff/+do19AQIDuvvvuHMPoFotF0dHROfad2zREfnXr1s1m6Ltx48aSlONE1ez2pKQkm/YmTZqoWbNmNm39+vVTWlqaduzYYW37+OOP1aZNG3l6eqp8+fKqUKGC3nnnHZspui+//FKScoy25VerVq30yy+/aPjw4frqq6+UlpZWoO1cK7djnZGRkWMaz1733nuvypUrZ7Nd6b8/Q4cPH9b+/fv18MMPS5KuXr1qvXXr1k3Jycm5Tsdcb+7cuapQoYJcXV3VuHFjbd26VZMmTdLw4cNt+kVGRsrLy8t638/PT9WrVy/Uz9b1+vXrZzOFFxgYqNatW1t/X66V/bqvfa6kXPte72a/n1u3btXZs2c1cOBAm/c1KytL99xzjxISEqyrkFq1aqVFixZpypQp+uGHH3TlyhU7XvGN7dixQw8++KDuvPNOTZ8+PcfjuU135uexm+natasSEhL0zTffaOrUqVqxYoV69eqlrKysAm8TN0fQgVPcd999CgoKUlxcnKR/zjFIT0/P9YO0Ro0aebZdP0xdpUoVm/tubm6SpL///tumf26rU/z9/a2P57dftooVK8rd3T3HvjMyMnI8P798fX1t7ru6ut6w/fp95ed9W7lypR588EHVqlVLixcv1rZt25SQkKBHH33UZnunTp1SuXLlct1mfowZM0YzZ87UDz/8oKioKFWpUkUdO3Ys1OUAbnasnbXdP//8U5I0atQoVahQweaWHVKunQLNy4MPPqiEhAT99NNPOnDggM6cOaPx48fftJ7smgr7Oq+V18/K9T/n5cuXz1FPXr+Lucnve9u7d+8c7+3LL78swzB09uxZSdKyZcs0cOBALViwQBEREfL19dWAAQOUkpKSn5ecp507d6pz58667bbbtHbtWmuN176G3F5rdl3X/37ao3LlymrZsqUiIyM1duxYzZ8/X2vWrMn1nDk4DmdAwSlcXFwUExOjsWPH6pVXXtHcuXPVsWNHNWzYMEff3P7jym7L7UPgRrL7Jycn69Zbb7V57MSJE6patWqOfte7tl9Jlp/3bfHixapTp46WLVtm85fo9SdRV6tWTZmZmUpJSbnpEubclC9fXrGxsYqNjdX58+e1fv16jR07Vl27dtWxY8dUsWJFu7dZXLKP/ZgxY/TAAw/k2ie3n+PrVatWTS1btnRobQWV18/K9b9fV69e1ZkzZ2zaC/q7mJvs93b27Nm68847c+2TPbpbtWpVzZo1S7NmzVJSUpLWrFmj0aNH6+TJk1q3bl2B9r9z50516tRJgYGB+vrrr+Xj45OjT9OmTbV06VJdvXrV5jydX3/9VZIUEhJSoH3nJvtk6IMHDzpsm8iJER04zWOPPSZXV1c9/PDDOnDgQJ4XRtuwYYP1Lz1JyszM1LJly1SvXr0cYeVm7r77bkn/fMBfKyEhQfv27VPHjh0lSREREfLw8MjR748//tA333xj7VeS7dmzR7/88otN25IlS+Tl5aXQ0FBJ/wyzu7q62oSclJSUHH9BRkVFSZLmzZt3w33mZ6ShUqVK6t27t2JiYnT27NmbrhRytMKOhjRs2FC33XabfvnlF7Vs2TLX27VTTaXB0qVLZRiG9f7vv/+urVu35nohyw8//NDmfvaiAkdc9LJNmzaqVKmS9u7dm+d7mz2Cea3atWvrySefVOfOnW2mZe2xa9cuderUSbfeeqvi4+NVuXLlXPv17NlTFy9e1IoVK2za33vvPfn7++uOO+4o0P5zkz0dWL9+fYdtEzkxogOnqVSpkgYMGKB58+YpMDAwxzku2apWraq7775b48eP1y233KK5c+dq//79eS4xv5GGDRvqX//6l2bPni0XFxdFRUXp6NGjGj9+vAICAqyrdipVqqTx48dr7NixGjBggPr27aszZ85o4sSJcnd314QJEwr12ouCv7+/evTooZdeekk1a9bU4sWLFR8fr5dfftk6gtK9e3etXLlSw4cPV+/evXXs2DFNnjxZNWvWtFlZ1q5dO/Xv319TpkzRn3/+qe7du8vNzU07d+5UxYoV9dRTT0n656/djz76SMuWLVPdunXl7u6upk2bKjo6WiEhIWrZsqV1CfWsWbMUGBio2267rUjfl6ZNm2rlypWaN2+ewsLC5OLiYvfIyltvvaWoqCh17dpVgwYNUq1atXT27Fnt27dPO3bs0Mcff+yk6p3j5MmT6tmzpx5//HGlpqZqwoQJcnd315gxY2z6ubq66pVXXtHFixcVHh6urVu3asqUKYqKilLbtm0LXYenp6dmz56tgQMH6uzZs+rdu7eqV6+uU6dO6ZdfftGpU6c0b948paamKjIyUv369VOjRo3k5eWlhIQErVu3zmaUbdOmTYqMjNSECRP00ksv5bnfAwcOqFOnTpKkqVOn6tChQzY//9lL/6V/Qn/nzp01bNgwpaWlqX79+lq6dKnWrVunxYsX25zfderUKX377beS/jvi8+WXX6patWqqVq2a2rdvL+mfn6ctW7aoS5cuCggIUHp6urZs2aLZs2erdevWuu+++wr93iJvBB04VZ8+fTRv3jwNGzYsz2tO9OjRQ02aNNH//u//KikpSfXq1dOHH36Y7yWt15s3b57q1aund955R3FxcfLx8dE999yj6dOn2wy/jxkzRtWrV9cbb7yhZcuWWa9lM23atCL/cC6I5s2ba/DgwZowYYIOHTokf39/vfrqqzZLsAcPHqyTJ0/qzTff1Lvvvqu6detq9OjR+uOPPzRx4kSb7S1atEihoaF65513tGjRInl4eCg4OFhjx4619pk4caKSk5P1+OOP68KFCwoMDNTRo0cVGRmpFStWaMGCBUpLS1ONGjXUuXNn61Voi9KIESO0Z88ejR07VqmpqTIMw2Y0Iz8iIyO1fft2TZ06VSNHjtS5c+dUpUoVBQcH68EHH3RS5c4zbdo0JSQkaPDgwUpLS1OrVq300Ucf5bgidYUKFfT555/r6aef1pQpU+Th4aHHH39c//nPfxxWyyOPPKLatWtrxowZeuKJJ3ThwgVVr15dzZs3ty4OcHd31x133KEPPvhAR48e1ZUrV1S7dm298MILev75563bunjxoqTcz7W71rZt26zn3eT2B9fChQttFiasXLlS48aN04svvqizZ8+qUaNGWrp0qR566CGb5+3ZsyfH9Xayz+Nq37699UrJTZs21eeff64xY8bo9OnTKl++vG677TaNHTtWsbGxXEfHySyGvf8DAHZ49tlnNW/ePB07dizXOX6LxaKYmBjNmTOnGKorvYKCghQSEpLjwnLAtbJHPD7++GP17t37hn0HDRqkTz75xBoeSoPnn39eS5cu1aFDh3IsFgCyESPhFD/88IMOHjyouXPn6oknnnDIiYwAcK2NGzdq/PjxhBzcEEEHThEREaGKFSuqe/fumjJlSnGXA8CEEhISirsElAJMXQEAANNieTkAADAtgg4AADAtgg4AADCtMn8yclZWlk6cOCEvL69CfVkbAAAoOoZh6MKFC/L398/zOm0SQUcnTpxQQEBAcZcBAAAK4NixYzf8uqAyH3Syv7Pm2LFj8vb2LuZqAABAfqSlpSkgIOCm3z1X5oNO9nSVt7c3QQcAgFLmZqedcDIyAAAwLYIOAAAwLYIOAAAwLYIOAAAwLYIOAAAwrTIbdOLi4hQcHKzw8PDiLgUAADhJmf/28rS0NPn4+Cg1NZXl5QAAlBL5/fwusyM6AADA/Mr8BQMBAIDjZWYZ2p54VicvZKi6l7ta1fFVOZei/05Jgg4AAHCodbuTNfGzvUpOzbC21fRx14ToYN0TUrNIa2HqCgAAOMy63ckatniHTciRpJTUDA1bvEPrdicXaT0EHQAA4BCZWYYmfrZXua1yym6b+NleZWYV3Toogg4AAHCI7Ylnc4zkXMuQlJyaoe2JZ4usJoIOAABwiJMX8g45BennCAQdAADgENW93B3azxEIOgAAwCFa1fFVTR935bWI3KJ/Vl+1quNbZDURdAAAgEOUc7FoQnSwJOUIO9n3J0QHF+n1dAg6AADAYe4Jqal5j4Sqho/t9FQNH3fNeyS0yK+jwwUDAQCAQ90TUlOdg2twZWQAAGBO5VwsiqhXpbjLYOoKAACYF0EHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYVqkPOgcOHFDz5s2tNw8PD61evbq4ywIAACVAqb9gYMOGDbVr1y5J0sWLFxUUFKTOnTsXb1EAAKBEKPUjOtdas2aNOnbsqFtuuaW4SwEAACVAsQedzZs3Kzo6Wv7+/rJYLLlOO82dO1d16tSRu7u7wsLCtGXLlly3tXz5cvXp08fJFQMAgNKi2INOenq6mjVrpjlz5uT6+LJlyzRy5EiNGzdOO3fuVLt27RQVFaWkpCSbfmlpafr+++/VrVu3G+7v0qVLSktLs7kBAABzshiGYRR3EdksFotWrVql+++/39p2xx13KDQ0VPPmzbO2NW7cWPfff7+mT59ubfvggw/01VdfafHixTfcx0svvaSJEyfmaE9NTZW3t3fhXwQAAHC6tLQ0+fj43PTzu9hHdG7k8uXL+vnnn9WlSxeb9i5dumjr1q02bfmdthozZoxSU1Ott2PHjjm0ZgAAUHKU6FVXp0+fVmZmpvz8/Gza/fz8lJKSYr2fmpqq7du3a8WKFTfdppubm9zc3BxeKwAAKHlKdNDJZrFYbO4bhmHT5uPjoz///LOoywIAACVciZ66qlq1qsqVK2czeiNJJ0+ezDHKAwAAcL0SHXRcXV0VFham+Ph4m/b4+Hi1bt26UNuOi4tTcHCwwsPDC7UdAABQchX71NXFixd1+PBh6/3ExETt2rVLvr6+ql27tmJjY9W/f3+1bNlSERERmj9/vpKSkjR06NBC7TcmJkYxMTHWs7YBAID5FHvQ+emnnxQZGWm9HxsbK0kaOHCgFi1apD59+ujMmTOaNGmSkpOTFRISorVr1yowMLC4SgYAAKVEibqOTnHI7zp8AABQcpjiOjrOxDk6AACYHyM6jOgAAFDqMKIDAADKPIIOAAAwLYIOAAAwLYIOAAAwrTIbdFh1BQCA+bHqilVXAACUOqy6AgAAZR5BBwAAmBZBBwAAmBZBBwAAmFaZDTqsugIAwPxYdcWqKwAASh1WXQEAgDKPoAMAAEyLoAMAAEyLoAMAAEyLoAMAAEyrzAYdlpcDAGB+LC9neTkAAKUOy8sBAECZR9ABAACmRdABAACmRdABAACmRdABAACmRdABAACmRdABAACmVWaDDhcMBADA/LhgIBcMBACg1OGCgQAAoMwj6AAAANMi6AAAANMi6AAAANMi6AAAANMi6AAAANMi6AAAANMi6AAAANMi6AAAANMqs0GHr4AAAMD8+AoIvgICAIBSh6+AAAAAZR5BBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmBZBBwAAmFaZDTpxcXEKDg5WeHh4cZcCAACcxGIYhlHcRRSntLQ0+fj4KDU1Vd7e3sVdDgAAyIf8fn6X2REdAABgfgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWgQdAABgWqYIOomJiYqMjFRwcLCaNm2q9PT04i4JAACUAOWLuwBHGDRokKZMmaJ27drp7NmzcnNzK+6SAABACVDqg86ePXtUoUIFtWvXTpLk6+tbzBUBAICSotinrjZv3qzo6Gj5+/vLYrFo9erVOfrMnTtXderUkbu7u8LCwrRlyxbrY4cOHZKnp6d69Oih0NBQTZs2rQirBwAAJVmxB5309HQ1a9ZMc+bMyfXxZcuWaeTIkRo3bpx27typdu3aKSoqSklJSZKkK1euaMuWLYqLi9O2bdsUHx+v+Pj4PPd36dIlpaWl2dwAAIA5FXvQiYqK0pQpU/TAAw/k+virr76qIUOG6LHHHlPjxo01a9YsBQQEaN68eZKkW2+9VeHh4QoICJCbm5u6deumXbt25bm/6dOny8fHx3oLCAhwxssCAAAlQLEHnRu5fPmyfv75Z3Xp0sWmvUuXLtq6daskKTw8XH/++afOnTunrKwsbd68WY0bN85zm2PGjFFqaqr1duzYMae+BgBFKzPL0LYjZ/TpruPaduSMMrOM4i4JQDEq0Scjnz59WpmZmfLz87Np9/PzU0pKiiSpfPnymjZtmu666y4ZhqEuXbqoe/fueW7Tzc2NVVmASa3bnayJn+1VcmqGta2mj7smRAfrnpCaxVgZgOJSooNONovFYnPfMAybtqioKEVFRRV1WQBKkHW7kzVs8Q5dP36TkpqhYYt3aN4joYQdoAwq0VNXVatWVbly5ayjN9lOnjyZY5QHQNmVmWVo4md7c4QcSda2iZ/tZRoLKIPyNaKzZs2afG+wR48eBS7meq6urgoLC1N8fLx69uxpbY+Pj9d9991XqG3HxcUpLi5OmZmZhS0TQDHbnnjWZrrqeoak5NQMbU88q4h6VYquMADFLl9B5/7777e5b7FYZBiGzf1s9gaHixcv6vDhw9b7iYmJ2rVrl3x9fVW7dm3Fxsaqf//+atmypSIiIjR//nwlJSVp6NChdu3nejExMYqJiVFaWpp8fHwKtS0AxevkhbxDTkH6ATCPfE1dZWVlWW9ff/21mjdvri+//FLnz59Xamqq1q5dq9DQUK1bt87uAn766Se1aNFCLVq0kCTFxsaqRYsWevHFFyVJffr00axZszRp0iQ1b95cmzdv1tq1axUYGGj3vgCYU3Uvd4f2A2AeFuPaoZl8CAkJ0Ztvvqm2bdvatG/ZskX/+te/tG/fPocW6GzZIzqpqany9vYu7nIAFEBmlqG2L3+jlNSMXM/TsUiq4eOu7164W+VcLLn0AFDa5Pfz2+6TkY8cOZLrVI+Pj4+OHj1q7+aKTVxcnIKDgxUeHl7cpQAopHIuFk2IDpb0T6i5Vvb9CdHBhBygDLJ7ROeuu+5ShQoVtHjxYtWs+c9SzZSUFPXv31+XL1/Wt99+65RCnYURHcA8uI4OUHbk9/Pb7uvovPvuu+rZs6cCAwNVu3ZtSVJSUpIaNGiQ6xdyAkBRuSekpjoH19D2xLM6eSFD1b3c1aqOLyM5QBlm94iO9M8F++Lj47V//34ZhqHg4GB16tQpx4X9SgNGdAAAKH3y+/ldoKCTLSMjQ25ubqUy4GQj6AAAUPo47WTkrKwsTZ48WbVq1ZKnp6cSExMlSePHj9c777xT8IoBAAAczO6gM2XKFC1atEgzZsyQq6urtb1p06ZasGCBQ4tzJlZdAQBgfnZPXdWvX19vvfWWOnbsKC8vL/3yyy+qW7eu9u/fr4iICJ07d85ZtToFU1cAAJQ+Tpu6On78uOrXr5+jPSsrS1euXLF3cwAAAE5jd9Bp0qSJtmzZkqP9448/tn6NAwAAQElg93V0JkyYoP79++v48ePKysrSypUrdeDAAb3//vv6/PPPnVEjAABAgdg9ohMdHa1ly5Zp7dq1slgsevHFF7Vv3z599tln6ty5szNqBAAAKJBCXUenNIuLi1NcXJwyMzN18OBBTkYGAKAUcdrJyHXr1tWZM2dytJ8/f15169a1d3PFJiYmRnv37lVCQkJxlwIAAJzE7qBz9OhRZWZm5mi/dOmSjh8/7pCiAAAAHCHfJyOvWbPG+u+vvvpKPj4+1vuZmZnasGGDgoKCHFocAABAYeQ76Nx///2SJIvFooEDB9o8VqFCBQUFBemVV15xaHEAAACFke+gk5WVJUmqU6eOEhISVLVqVacVBQAA4Ah2X0cn+0s8AQAASjq7T0Z++umn9cYbb+RonzNnjkaOHOmImooEX+oJAID52X0dnVq1amnNmjUKCwuzad+xY4d69OihP/74w6EFOhtf6gkAQOnjtOvonDlzxmbFVTZvb2+dPn3a3s0BAAA4jd1Bp379+lq3bl2O9i+//LJUXTAQAACYn90nI8fGxurJJ5/UqVOndPfdd0uSNmzYoFdeeUWzZs1ydH0AnCgzy9D2xLM6eSFD1b3c1aqOr8q5WIq7LABwGLuDzqOPPqpLly5p6tSpmjx5siQpKChI8+bN04ABAxxeIADnWLc7WRM/26vk1AxrW00fd02IDtY9ITWLsTIAcJxCfannqVOn5OHhIU9PT0fWVKQ4GRll0brdyRq2eIeu/+XPHsuZ90goYQdAiea0k5El6erVq1q/fr1Wrlyp7Jx04sQJXbx4sWDVAigymVmGJn62N0fIkWRtm/jZXmVmFfhvIAAoMeyeuvr99991zz33KCkpSZcuXVLnzp3l5eWlGTNmKCMjQ2+++aYz6gTgINsTz9pMV13PkJScmqHtiWcVUa9K0RUGAE5g94jOiBEj1LJlS507d04eHh7W9p49e2rDhg0OLc6ZuGAgyqqTF/IOOQXpBwAlmd0jOt99952+//57ubq62rQHBgbq+PHjDivM2WJiYhQTE2Od4wPKiupe7g7tBwAlmd0jOllZWcrMzMzR/scff8jLy8shRQFwnlZ1fFXTx115LSK36J/VV63q+BZlWQDgFHYHnc6dO9tcL8disejixYuaMGGCunXr5sjaADhBOReLJkQHS1KOsJN9f0J0MNfTAWAKdi8vP3HihCIjI1WuXDkdOnRILVu21KFDh1S1alVt3rxZ1atXd1atTsHycpRVXEcHQGmW38/vAl1H5++//9ZHH32kn3/+WVlZWQoNDdXDDz9sc3JyaUHQQVnGlZEBlFYODTqhoaHasGGDKleurEmTJmnUqFGqWLGiQwsuLgQdAABKH4deMHDfvn1KT0+XJE2cOJELAwIAgFIhX8vLmzdvrsGDB6tt27YyDEMzZ87M82sfXnzxRYcWCAAAUFD5mro6cOCAJkyYoCNHjmjHjh0KDg5W+fI5M5LFYtGOHTucUqizMHUFAEDp47STkV1cXJSSklLqVlflhaADAEDp49BzdEJDQ3Xu3DlJ0oQJE0r1t5UDAICyw+6TkSdNmmSKk5HN8F1XmVmGth05o093Hde2I2f4tmkAAK6Tr6mriIgIeXp6qm3btpo4caJGjRplmpORS+vUFRd7AwCUZQ49R4eTkUuWdbuTNWzxDl1/4LIv8zbvkVDCDgDA1DgZOZ9KW9DJzDLU9uVvbEZyrmWRVMPHXd+9cDdXuAUAmJZDT0a+VlZWVo6Qc+7cOc2ePVvNmze3u1DYZ3vi2TxDjiQZkpJTM7Q98WzRFQUAQAmVrwsG5mX9+vV65513tHr1alWtWlUPPPCAo+pCHk5eyDvkFKQfAABmZnfQSUpK0sKFC7Vw4UJdvHhR586d0/Lly9WrVy9n1IfrVPdyd2g/AADMLN9TV8uXL1eXLl3UuHFj7d69W6+//rpOnDghFxcXNW7c2Jk14hqt6viqpo+78jr7xqJ/Vl+1quNblGUVCsvkAQDOku8RnX79+un555/XihUr5OXl5cyacAPlXCyaEB2sYYt3yCLZrLzKDj8TooNLzYnILJMHADhTvkd0Hn30Uc2dO1f33HOP3nzzTeuVklH07gmpqXmPhKqGj+30VA0f91K1tDx7mfz1J1enpGZo2OIdWrc7uZgqAwCYhV3Ly//++28tX75c7777rn788Ud17dpVX3zxhXbt2qWQkBBn1uk0pW15+bUyswxtTzyrkxcyVN3rn+mq0jKSwzJ5AEBhOGV5uYeHhwYOHKhvv/1Wv/76q4KDg+Xn56c2bdqoX79+WrlyZaELR/6Vc7Eool4V3de8liLqVSlVgYBl8gCAomD3dXSy3XbbbZo+fbqOHTumxYsX66+//lLfvn0dWRtMjGXyAICiUKjr6Ej/XCk5Ojpa0dHROnnypCNqQhnAMnkAQFEo8IhObszytRBwPjMukwcAlDwODTpAfmUvk5eUI+yUxmXyAICSqcwGnbi4OAUHBys8PLy4SymzzLJMHgBQctn97eVmU5qXl5tFaV4mDwAoHvn9/C7Qycjnz5/XJ598oiNHjui5556Tr6+vduzYIT8/P9WqVavARaNsyl4mDwCAo9kddP7v//5PnTp1ko+Pj44eParHH39cvr6+WrVqlX7//Xe9//77zqgTAADAbnafoxMbG6tBgwbp0KFDcnf/77kVUVFR2rx5s0OLAwAAKAy7g05CQoKeeOKJHO21atVSSkqKQ4oCAABwBLuDjru7u9LS0nK0HzhwQNWqVXNIUQAAAI5gd9C57777NGnSJF25ckWSZLFYlJSUpNGjR6tXr14OLxAAAKCg7A46M2fO1KlTp1S9enX9/fffat++verXry8vLy9NnTrVGTUCAAAUiN2rrry9vfXdd9/pm2++0Y4dO5SVlaXQ0FB16tTJGfUBAAAUGBcM5IKBAACUOg69YOAbb7yR7x0//fTT+e4LAADgTPka0alTp07+Nmax6Lfffit0UUWJER0AAEofh47oJCYmOqwwAACAolKoby83DENl/BQfAABQghUo6LzzzjsKCQmRu7u73N3dFRISogULFji6NgAAgEKxe3n5+PHj9dprr+mpp55SRESEJGnbtm165plndPToUU2ZMsXhRQIAABSE3cvLq1atqtmzZ6tv37427UuXLtVTTz2l06dPO7RAZ+NkZAAASp/8fn7bPXWVmZmpli1b5mgPCwvT1atX7d0cAACA09gddB555BHNmzcvR/v8+fP18MMPO6QoAAAAR7D7HB3pn5ORv/76a915552SpB9++EHHjh3TgAEDFBsba+336quvOqZKAACAArA76OzevVuhoaGSpCNHjkiSqlWrpmrVqmn37t3WfhaLxUElAgAAFIzdQWfjxo3OqKNQypcvr5CQEElSy5YtWeoOAAAkFXDqKtsff/whi8WiWrVqOaqeAqlUqZJ27dpVrDUAAICSx+6TkbOysjRp0iT5+PgoMDBQtWvXVqVKlTR58mRlZWU5o0YAAIACsTvojBs3TnPmzNG///1v7dy5Uzt27NC0adM0e/ZsjR8/3u4CNm/erOjoaPn7+8tisWj16tU5+sydO1d16tSRu7u7wsLCtGXLFpvH09LSFBYWprZt2+rbb7+1uwYAAGBOdk9dvffee1qwYIF69OhhbWvWrJlq1aql4cOHa+rUqXZtLz09Xc2aNdPgwYPVq1evHI8vW7ZMI0eO1Ny5c9WmTRu99dZbioqK0t69e1W7dm1J0tGjR+Xv76/du3fr3nvv1a+//prnxYMuXbqkS5cuWe+npaXZVS8AACg97B7ROXv2rBo1apSjvVGjRjp79qzdBURFRWnKlCl64IEHcn381Vdf1ZAhQ/TYY4+pcePGmjVrlgICAmyu5ePv7y9JCgkJUXBwsA4ePJjn/qZPny4fHx/rLSAgwO6aAQBA6WB30GnWrJnmzJmTo33OnDlq1qyZQ4rKdvnyZf3888/q0qWLTXuXLl20detWSdK5c+esIzR//PGH9u7dq7p16+a5zTFjxig1NdV6O3bsmENrBgAAJYfdU1czZszQvffeq/Xr1ysiIkIWi0Vbt27VsWPHtHbtWocWd/r0aWVmZsrPz8+m3c/PTykpKZKkffv26YknnpCLi4ssFotef/11+fr65rlNNzc3ubm5ObROAABQMtkddNq3b6+DBw8qLi5O+/fvl2EYeuCBBzR8+HDrFJKjXX/xQcMwrG2tW7fWr7/+6pT9AgCA0q1A19Hx9/e3+6TjgqhatarKlStnHb3JdvLkyRyjPAAAANez+xwdSdqyZYseeeQRtW7dWsePH5ckffDBB/ruu+8cWpyrq6vCwsIUHx9v0x4fH6/WrVsXattxcXEKDg5WeHh4obYDAABKLruDzooVK9S1a1d5eHhox44d1hOBL1y4oGnTptldwMWLF7Vr1y7rlY0TExO1a9cuJSUlSZJiY2O1YMECvfvuu9q3b5+eeeYZJSUlaejQoXbv61oxMTHau3evEhISCrUdAABQclkMwzDseUKLFi30zDPPaMCAAfLy8tIvv/yiunXrateuXbrnnntyTDPdzKZNmxQZGZmjfeDAgVq0aJGkfy4YOGPGDCUnJyskJESvvfaa7rrrLrv2k5e0tDT5+PgoNTU1z2vvAACAkiW/n992B52KFStq7969CgoKsgk6v/32m4KDg5WRkVHo4osSQQcAgNInv5/fdk9d1axZU4cPH87R/t13393w+jUlDefoAABgfnYHnSeeeEIjRozQjz/+KIvFohMnTujDDz/UqFGjNHz4cGfU6BScowMAgPnZvbz8+eefV2pqqiIjI5WRkaG77rpLbm5uGjVqlJ588kln1AgAAFAgdp+jk+2vv/7S3r17lZWVpeDgYHl6ejq6tiLBOToAAJQ+Dj9H56+//lJMTIxq1aql6tWr67HHHlNQUJBatWpVakMOAAAwt3wHnQkTJmjRokW699579dBDDyk+Pl7Dhg1zZm0AAACFku9zdFauXKl33nlHDz30kCTpkUceUZs2bZSZmaly5co5rUBniYuLU1xcnDIzM4u7FAAA4CT5PkfH1dVViYmJqlWrlrXNw8NDBw8eVEBAgNMKdDbO0QEAoPRx+Dk6mZmZcnV1tWkrX768rl69WvAqAQAAnCjfU1eGYWjQoEFyc3OztmVkZGjo0KG65ZZbrG0rV650bIUAAAAFlO+gM3DgwBxtjzzyiEOLAQAAcKR8B52FCxc6sw4AAACHs/srIMyC77oCAMD8CnxlZLNg1RUAAKWP0769HAAAoLQg6AAAANMi6AAAANMi6AAAANMi6AAAANMqs0GH5eUAAJgfy8tZXg4AQKnD8nIAAFDmEXQAAIBpEXQAAIBpEXQAAIBpEXQAAIBpEXQAAIBpEXQAAIBpldmgwwUDAQAwPy4YyAUDAQAodbhgIAAAKPMIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLQIOgAAwLTKbNDhu64AADA/vuuK77oCAKDU4buuAABAmUfQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAApkXQAQAAplVmg05cXJyCg4MVHh5e3KUAAAAnsRiGYRR3EcUpLS1NPj4+Sk1Nlbe3d3GXAwAA8iG/n99ldkQHAACYH0EHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYlmmCzl9//aXAwECNGjWquEsBAAAlhGmCztSpU3XHHXcUdxkAAKAEMUXQOXTokPbv369u3boVdykAAKAEKfags3nzZkVHR8vf318Wi0WrV6/O0Wfu3LmqU6eO3N3dFRYWpi1bttg8PmrUKE2fPr2IKgYAAKVFsQed9PR0NWvWTHPmzMn18WXLlmnkyJEaN26cdu7cqXbt2ikqKkpJSUmSpE8//VQNGjRQgwYN8rW/S5cuKS0tzeYGAADMyWIYhlHcRWSzWCxatWqV7r//fmvbHXfcodDQUM2bN8/a1rhxY91///2aPn26xowZo8WLF6tcuXK6ePGirly5omeffVYvvvhirvt46aWXNHHixBztqamp8vb2dvhrAgAAjpeWliYfH5+bfn6X6KBz+fJlVaxYUR9//LF69uxp7TdixAjt2rVL3377rc3zFy1apN27d2vmzJl57uPSpUu6dOmS9X5aWpoCAgIIOgAAlCL5DTrli7Amu50+fVqZmZny8/Ozaffz81NKSkqBtunm5iY3NzdHlAcAAEq4Eh10slksFpv7hmHkaJOkQYMGFVFFAACgNCj2k5FvpGrVqipXrlyO0ZuTJ0/mGOUBAAC4XokOOq6urgoLC1N8fLxNe3x8vFq3bl2obcfFxSk4OFjh4eGF2g4AACi5in3q6uLFizp8+LD1fmJionbt2iVfX1/Vrl1bsbGx6t+/v1q2bKmIiAjNnz9fSUlJGjp0aKH2GxMTo5iYGOvJTAAAwHyKPej89NNPioyMtN6PjY2VJA0cOFCLFi1Snz59dObMGU2aNEnJyckKCQnR2rVrFRgYWFwlAwCAUqJELS8vDvldngYAAEqO/H5+l+hzdJyJc3QAADA/RnQY0QEAoNRhRAcAAJR5BB0AAGBaBB0AAGBaBB0AAGBaZTbosOoKAADzY9UVq64AACh1WHUFAADKPIIOAAAwLYIOAAAwLYIOAAAwrTIbdFh1BQCA+bHqilVXAACUOqy6AgAAZR5BBwAAmBZBBwAAmBZBBwAAmBZBBwAAmFaZDTosLwcAwPxYXs7ycgAASh2WlwMAgDKPoAMAAEyLoAMAAEyLoAMAAEyLoAMAAEyLoAMAAEyrfHEXYEaZWYa2J57VyQsZqu7lrlZ1fFXOxVLcZQEAUOaU2aATFxenuLg4ZWZmOnS763Yna+Jne5WcmmFtq+njrgnRwbonpKZD9wUAAG6MCwY68IKB63Yna9jiHbr+Dc0ey5n3SChhBwAAB+CCgUUsM8vQxM/25gg5kqxtEz/bq8ysMp0rAQAoUgQdB9meeNZmuup6hqTk1AxtTzxbdEUBAFDGEXQc5OSFvENOQfoBAIDCI+g4SHUvd4f2AwAAhUfQcZBWdXxV08ddeS0it+if1Vet6vgWZVkAAJRpBB0HKedi0YToYEnKEXay70+IDuZ6OgAAFCGCjgPdE1JT8x4JVQ0f2+mpGj7uLC0HAKAYlNkLBjrLPSE11Tm4BldGBgCgBCDoOEE5F4si6lUp7jIAACjzmLoCAACmVWaDTlxcnIKDgxUeHl7cpQAAACfhu64c+F1XAACgaPBdVwAAoMwj6AAAANMi6AAAANMi6AAAANMi6AAAANMi6AAAANMq81dGzl5dn5aWVsyVAACA/Mr+3L7ZVXLKfNC5cOGCJCkgIKCYKwEAAPa6cOGCfHx88ny8zF8wMCsrSydOnJCXl5cslv9+8WZ4eLgSEhJy9M9ve1pamgICAnTs2LFivxBhXjUX9fbseV5++t6oT0Ee4xg69nlFfQxzazPrMSzMtvL73MIevxs9zu9g6fgdvNHjJeEYGoahCxcuyN/fXy4ueZ+JU+ZHdFxcXHTrrbfmaC9XrlyuB8Xedm9v72L/Bc2rtqLenj3Py0/fG/UpyGMcQ8c+r6iP4Y36m+0YFmZb+X1uYY/fjR7nd7B0/A7e6PGScgxvNJKTjZOR8xATE+OQ9pLA0bUVdHv2PC8/fW/UpyCPcQwd+7yiPoYl+fhJjq2vMNvK73MLe/xu9Di/g6Xjd/BGj5emY1jmp66che/QKv04hqUfx7B04/iVfiXhGDKi4yRubm6aMGGC3NzcirsUFBDHsPTjGJZuHL/SryQcQ0Z0AACAaTGiAwAATIugAwAATIugAwAATIugAwAATIugAwAATIugU0wSExMVGRmp4OBgNW3aVOnp6cVdEvLpwIEDat68ufXm4eGh1atXF3dZsNNrr72mJk2aKDg4WE8//fRNvxgQJc/MmTPVpEkThYSEaPHixcVdDvKhZ8+eqly5snr37m3T/vnnn6thw4a67bbbtGDBAofuk+XlxaR9+/aaMmWK2rVrp7Nnz8rb21vly5f5b+QodS5evKigoCD9/vvvuuWWW4q7HOTTqVOndOedd2rPnj2qUKGC7rrrLs2cOVMRERHFXRry6ddff9XAgQO1detWSVLHjh31xRdfqFKlSsVbGG5o48aNunjxot577z198sknkqSrV68qODhYGzdulLe3t0JDQ/Xjjz/K19fXIftkRKcYZP/n2q5dO0mSr68vIaeUWrNmjTp27EjIKYWuXr2qjIwMXblyRVeuXFH16tWLuyTYYd++fWrdurXc3d3l7u6u5s2ba926dcVdFm4iMjJSXl5eNm3bt29XkyZNVKtWLXl5ealbt2766quvHLZPgk4BbN68WdHR0fL395fFYsl12mLu3LmqU6eO3N3dFRYWpi1btlgfO3TokDw9PdWjRw+FhoZq2rRpRVg9Cnv8rrV8+XL16dPHyRXjeoU9htWqVdOoUaNUu3Zt+fv7q1OnTqpXr14RvgIU9hiGhIRo48aNOn/+vM6fP69vvvlGx48fL8JXUPY48v/Oa504cUK1atWy3r/11lsdeiwJOgWQnp6uZs2aac6cObk+vmzZMo0cOVLjxo3Tzp071a5dO0VFRSkpKUmSdOXKFW3ZskVxcXHatm2b4uPjFR8fX5QvoUwr7PHLlpaWpu+//17dunUrirJxjcIew3Pnzunzzz/X0aNHdfz4cW3dulWbN28uypdQ5hX2GGafW3X33XerZ8+eCg8PZ2TcyRz1f+f1cjuDxmKxOKTm7B2gECQZq1atsmlr1aqVMXToUJu2Ro0aGaNHjzYMwzC2bt1qdO3a1frYjBkzjBkzZji9VuRUkOOX7f333zcefvhhZ5eImyjIMVy+fLkxfPhw62MzZswwXn75ZafXitwV5vcw25AhQ4zPP//cWSXiOoU5Zhs3bjR69eplvf/9998b999/v/X+008/bXz44YcOq5URHQe7fPmyfv75Z3Xp0sWmvUuXLtaT5sLDw/Xnn3/q3LlzysrK0ubNm9W4cePiKBfXyc/xy8a0VcmUn2MYEBCgrVu3KiMjQ5mZmdq0aZMaNmxYHOUiF/n9PTx58qSkf1ZCbt++XV27di3SOvFf9vzfeb1WrVpp9+7dOn78uC5cuKC1a9c69Fgyzudgp0+fVmZmpvz8/Gza/fz8lJKSIkkqX768pk2bprvuukuGYahLly7q3r17cZSL6+Tn+ElSamqqtm/frhUrVhR1ibiJ/BzDO++8U926dVOLFi3k4uKijh07qkePHsVRLnKR39/D+++/X+fPn9ctt9yihQsXMnVVjPJ7zLp27aodO3YoPT1dt956q1atWqXw8HC98sorioyMVFZWlp5//nlVqVLFYbXxU+Ek188vGoZh0xYVFaWoqKiiLgv5dLPj5+Pjoz///LOoy4IdbnYMp06dqqlTpxZ1WbDDzY7hzUYKUPRudszyWk3Vo0cPp/2xwdSVg1WtWlXlypWzSbDSP0Os1yddlDwcv9KPY1j6cQxLn5J8zAg6Dubq6qqwsLAcq6ji4+PVunXrYqoK+cXxK/04hqUfx7D0KcnHjKmrArh48aIOHz5svZ+YmKhdu3bJ19dXtWvXVmxsrPr376+WLVsqIiJC8+fPV1JSkoYOHVqMVSMbx6/04xiWfhzD0qfUHjOHrd8qQzZu3GhIynEbOHCgtU9cXJwRGBhouLq6GqGhoca3335bfAXDBsev9OMYln4cw9KntB4zvusKAACYFufoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoAAAA0yLoACiwoKAgzZo1q7jLKJEWLVoki8Uii8WikSNHFnc51loqVapU3KUARYqgA5hA9odYXrdBgwbd9PmrV692eF0vvfSStYby5curatWquuuuuzRr1ixdunTJ4fsraby9vZWcnKzJkydb2/IKh7NmzVJQUJD1fnp6ul544QXVrVtX7u7uqlatmjp06KDPP//c2qdDhw7W99fNzU21atVSdHS0Vq5cmWP7ycnJhFKUSXx7OWACycnJ1n8vW7ZML774og4cOGBt8/DwKI6yJElNmjTR+vXrlZWVpTNnzmjTpk2aMmWKPvjgA23atEleXl7FVpuzWSwW1ahRo0DPHTp0qLZv3645c+YoODhYZ86c0datW3XmzBmbfo8//rgmTZqkK1eu6Pjx41q1apUeeughDRo0SPPnz7f2q1Gjhnx8fAr1eoDSiBEdwARq1Khhvfn4+Fg/YLNvS5YsUb169eTq6qqGDRvqgw8+sD43exShZ8+eslgs1vtHjhzRfffdJz8/P3l6eio8PFzr16+3u7by5curRo0a8vf3V9OmTfXUU0/p22+/1e7du/Xyyy9b+y1evFgtW7aUl5eXatSooX79+unkyZPWxzdt2iSLxaKvvvpKLVq0kIeHh+6++26dPHlSX375pRo3bixvb2/17dtXf/31l/V569atU9u2bVWpUiVVqVJF3bt315EjR2xq3Lp1q5o3by53d3e1bNlSq1evlsVi0a5du6x99u7dq27dusnT01N+fn7q37+/Tp8+bff7kV+fffaZxo4dq27duikoKEhhYWF66qmnNHDgQJt+FStWVI0aNRQQEKA777xTL7/8st566y29/fbbBTpegNkQdACTW7VqlUaMGKFnn31Wu3fv1hNPPKHBgwdr48aNkqSEhARJ0sKFC5WcnGy9f/HiRXXr1k3r16/Xzp071bVrV0VHRyspKanQNTVq1EhRUVE2UyyXL1/W5MmT9csvv2j16tVKTEzMdcrtpZde0pw5c7R161YdO3ZMDz74oGbNmqUlS5boiy++UHx8vGbPnm3tn56ertjYWCUkJGjDhg1ycXFRz549lZWVJUm6cOGCoqOj1bRpU+3YsUOTJ0/WCy+8YLPP5ORktW/fXs2bN9dPP/2kdevW6c8//9SDDz5Y6PciLzVq1NDatWt14cIFu587cOBAVa5cOdcpLKCsYeoKMLmZM2dq0KBBGj58uCQpNjZWP/zwg2bOnKnIyEhVq1ZNklSpUiWbaZZmzZqpWbNm1vtTpkzRqlWrtGbNGj355JOFrqtRo0b6+uuvrfcfffRR67/r1q2rN954Q61atdLFixfl6elpU0ebNm0kSUOGDNGYMWN05MgR1a1bV5LUu3dvbdy40RpWevXqZbPfd955R9WrV9fevXsVEhKiDz/8UBaLRW+//bbc3d0VHBys48eP6/HHH7c+Z968eQoNDdW0adOsbe+++64CAgJ08OBBNWjQoNDvx/Xmz5+vhx9+WFWqVFGzZs3Utm1b9e7d2/rab8TFxUUNGjTQ0aNHHV4XUNowogOY3L59+3J8OLZp00b79u274fPS09P1/PPPKzg4WJUqVZKnp6f279/vkBEdSTIMQxaLxXp/586duu+++xQYGCgvLy916NBBknLs7/bbb7f+28/PTxUrVrSGnOy2a6e8jhw5on79+qlu3bry9vZWnTp1bLZ74MAB3X777XJ3d7c+p1WrVjb7/Pnnn7Vx40Z5enpab40aNbJu3xnuuusu/fbbb9qwYYN69eqlPXv2qF27djYnNt/I9e8vUFYxogOUAdd/4OXnQ/C5557TV199pZkzZ6p+/fry8PBQ7969dfnyZYfUtG/fPmvoSE9PV5cuXdSlSxctXrxY1apVU1JSkrp27ZpjfxUqVLB5Xdfez27LnpaSpOjoaAUEBOjtt9+Wv7+/srKyFBISYt1ubu+FYRg297OyshQdHW1zTlG2mjVr2vW6vb29lZqamqP9/PnzOU4WrlChgtq1a6d27dpp9OjRmjJliiZNmqQXXnhBrq6uee4jMzNThw4dUnh4uF21AWZE0AFMrnHjxvruu+80YMAAa9vWrVvVuHFj6/0KFSooMzPT5nlbtmzRoEGD1LNnT0n/nLPjqKmQ/fv3a926dRozZoz1/unTp/Xvf/9bAQEBkqSffvqp0Ps5c+aM9u3bp7feekvt2rWTJH333Xc2fRo1aqQPP/xQly5dkpubW677Dg0N1YoVKxQUFKTy5Qv332ajRo2s50FdKyEhQQ0bNrzhc4ODg3X16lVlZGTcMOi89957OnfuXI5pO6AsIugAJvfcc8/pwQcfVGhoqDp27KjPPvtMK1eutFmRExQUpA0bNqhNmzZyc3NT5cqVVb9+fa1cuVLR0dGyWCwaP368zUhJfl29elUpKSk5lpc3b95czz33nCSpdu3acnV11ezZszV06FDt3r0731M0N1K5cmVVqVJF8+fPV82aNZWUlKTRo0fb9OnXr5/GjRunf/3rXxo9erSSkpI0c+ZMSf8dCYuJidHbb7+tvn376rnnnlPVqlV1+PBhffTRR3r77bdVrly5fNcUGxurNm3aaNKkSerdu7ckacWKFVq3bp22bt1q7dehQwf17dtXLVu2VJUqVbR3716NHTtWkZGR8vb2tvb766+/lJKSoqtXr+r48eNauXKlXnvtNQ0bNkyRkZEFfu8A0zAAmMrChQsNHx8fm7a5c+cadevWNSpUqGA0aNDAeP/9920eX7NmjVG/fn2jfPnyRmBgoGEYhpGYmGhERkYaHh4eRkBAgDFnzhyjffv2xogRI6zPCwwMNF577bU8a5kwYYIhyZBklCtXzvD19TXatm1rvPbaa0ZGRoZN3yVLlhhBQUGGm5ubERERYaxZs8aQZOzcudMwDMPYuHGjIck4d+7cDV/rhAkTjGbNmlnvx8fHG40bNzbc3NyM22+/3di0aZMhyVi1apW1z/fff2/cfvvthqurqxEWFmYsWbLEkGTs37/f2ufgwYNGz549jUqVKhkeHh5Go0aNjJEjRxpZWVm5vvbcaru2pnbt2hmVK1c2KleubLRt29aIj4+36TNt2jQjIiLC8PX1Ndzd3Y26desaTz/9tHH69Glrn/bt21vfX1dXV6NmzZpG9+7djZUrV9pdE2BWFsO4bjIaAMq4Dz/8UIMHD1ZqamqBL7a4aNEijRw5UufPn3dscYVQEmsCnI2pKwBl3vvvv6+6deuqVq1a+uWXX/TCCy/owQcfLPQVpVNTU+Xp6amYmJhcT2QuSp6enrp69arN6jKgLCDoACjzUlJS9OKLLyolJUU1a9bU//zP/2jq1KmF2mavXr3Utm1bSSoRX6SZfZVne84nAsyAqSsAAGBaXDAQAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACYFkEHAACY1v8DnsWK5yfpn98AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\"\"\"Total people affected by TCs in the Philippines in 2013:\"\"\"\n", "\n", "# People affected\n", "impact_emdat_PHL, countries = emdat_to_impact(\n", " emdat_file_path,\n", " \"TC\",\n", " countries=\"PHL\",\n", " year_range=(2013, 2013),\n", " imp_str=\"Total Affected\",\n", ")\n", "\n", "print(\n", " \"Number of TC events in EM-DAT in the Philipppines, 2013: %i\"\n", " % (impact_emdat_PHL.event_id.size)\n", ")\n", "print(\"\\nPeople affected by TC events in the Philippines in 2013 (per event):\")\n", "print(impact_emdat_PHL.at_event)\n", "print(\"\\nPeople affected by TC events in the Philippines in 2013 (total):\")\n", "print(int(impact_emdat_PHL.aai_agg))\n", "\n", "# Comparison to monetary damages:\n", "impact_emdat_PHL_USD, _ = emdat_to_impact(\n", " emdat_file_path, \"TC\", countries=\"PHL\", year_range=(2013, 2013)\n", ")\n", "\n", "ax = plt.scatter(impact_emdat_PHL_USD.at_event, impact_emdat_PHL.at_event)\n", "plt.title(\"Typhoon impacts in the Philippines, 2013\")\n", "plt.xlabel(\"Total Damage [USD]\")\n", "plt.ylabel(\"People Affected\")\n", "plt.xscale(\"log\")\n", "plt.yscale(\"log\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### emdat_impact_yearlysum()\n", "\n", "function to load EM-DAT impact data and return DataFrame with impact summed per year and country\n", "\n", "#### Parameters:\n", "- emdat_file_csv (str): Full path to EMDAT-file (CSV)\n", "\n", "#### Optional parameters:\n", "\n", "- hazard (list or str): List of Disaster (sub-)type according EMDAT terminology or CLIMADA hazard type abbreviations. e.g. ['Wildfire', 'Forest fire'] or ['BF']\n", "- year_range (list with 2 integers): start and end year e.g. [1980, 2017]\n", "- countries (list of str): country ISO3-codes or names, e.g. ['JAM', 'CUB']. Set to None or ['all'] for all countries \n", "- reference_year (int): reference year of exposures for normalization. Impact is scaled proportional to GDP to the value of the reference year. No scaling for reference_year=0 (default)\n", "- imp_str (str): Column name of impact metric in EMDAT CSV, e.g. 'Total Affected'; default = \"Total Damages\"\n", "- version (int): given EM-DAT data format version (i.e. year of download), changes naming of columns/variables (default: 2020)\n", "\n", "#### Returns:\n", "- pandas.DataFrame with impact per year and country" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2014\n", " 2015 2016 2017 2018 2019 2020]\n" ] }, { "data": { "text/plain": [ "Text(0, 0.5, 'Total Damage [USD]')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOydd3hUVdrAf29CpPcgIi2hE0KCgFSlWABRURAQRUXUxc7adVkVy64FXT51dRVURJBmF7siIgKC0gWRIjX0XkNIeb8/zp0wSWYmk2RakvN7nvtM5p5zz31n5ua+95y3iapisVgsltJLVLgFsFgsFkt4sYrAYrFYSjlWEVgsFkspxyoCi8ViKeVYRWCxWCylHKsILBaLpZRjFUExQ0RURJqEW47SjIjcKCLzvLQ1EJFjIhIdarkiCREZKiLfBXC8OSJyS6DGs+TEKoIA4twAXFuWiKS6vR/q9GkmIh+IyD4ROSwiK0XkvtJ+4wglwbypqOpWVa2kqpmFkKuHc90cy7V1dtrnOA8CybmO+9TZ38PLuHFOu2u83SLyhYhc7KX/HBE5KCJl3fZ97XZ8uoiccnv/hofvYYqq9irod+Cc6wkRea8wx/o5fp6HqdznFJFRIrLJ+XwpIjLDwzgTRSRDRM4OlqyhwiqCAOLcACqpaiVgK3C5274pItIYWARsA1qralVgENAeqBw+ycOHiJQJ4blERCL9mt/hfh052y9u7euAG1xvRKQm0AnY68fY1ZxrMxn4HvhERG507yAiccD5gAL9XPtV9RK3a3sKMMZNvtsK80EjFREZBlwPXOR83vbAD7n6VASuAg4DQ0MuZICJ9H+KksaTwAJVvU9VdwKo6lpVvVZVD3k6QEQeFJGdIrJDRG7K1XapiCwTkSMisk1EnnBrcz0FDnfaDorIbSJyrjMLOSQir7r1bywis0VkvzNbmSIi1dza2zrnOurMaGaIyL/c2i8TkeXOuAtEJMnbl+DIdaeIrAfWO/taiMj3InJARNaKyGC3/hNF5A2n/aiI/CQiDd3au4jIb84M6zcR6eLWNkdE/i0i84ETwGTMje5V52nvVT/OX1NEZjrf869AYx+fzfW9l3E7/9MiMt+R/TsRifV2vB9MAa6W0zPIa4BPgFP+DqCqu1T1ZeAJ4PlcyvEGYCEwERhWWCEl1/KZ853cJiLrnWvxNRERD8f1AUZhPuMxEVnh1tzQ2/coIp2c6+6QiKwQL7MjPzkX+FZV/4Ls72t8rj5XAYeApyjC9xQxqKrdgrABmzFPFO77dgHDCzBGH2A3kAhUBKZintSaOO09gNYYhZ7k9L3SaYtz+r4BlAN6ASeBT4EzgbrAHqC7078JcDFQFqgFzAVectrOALYAfwdigAGYG8+/nPa2zlgdgWjMP8ZmoKyXz6WYJ9IaQHnns20DhgNlnPH2Aa2c/hOBo0A3R76XgXlOWw3gIOYJrgzmxngQqOm0z8HMzlo57THOvlvc5Mnv/NOB951+icB21/k9fDbX917G7fx/Ac2czzoHeM7LsT2AFB/XwxzgFuA74BJn369AZyAF6OGPTG77Gzn7W7rt2wDcAbQD0oHaHsab6Prtfch6o/t35JznC6Aa0AAzg+nj5dgngPc8fHaP3yPmWt4P9MX8L1zsvK/l4/pr4u2cwHXAAeBBzGwg2sMYPwBjgNpABtA23PecomzFckYgIhNEZI+IrPKjbzcRWeqs5Q3M1faN8wTxRfCkzUFNYGcB+g8G3lHVVap6HHOxZqOqc1T1d1XNUtWVwDSge64xnlbVk6r6HXAcmKaqe1R1O/AzcI4z1gZV/V5V01R1LzDWbaxOmBvkK6qarqofY25ALv4GjFPVRaqaqarvAmnOcd54VlUPqGoqcBmwWVXfUdUMVV0KfAS4/15fqupcVU0D/gl0FpH6wKXAelWd7Bw7DfgTuNzt2ImqutppT/cgi9fzO0/eVwGPq+pxVV0FvOvjc3niHVVd53zW94E2Pvqe7VyT7lvFXH0mATeISHPMcs8vHsbxhx3Oaw0AETkPaAi8r6pLMDfeaws5tieeU9VDqroV+BHf34MnvH2P1wFfqepXzv/C98BijGIoMKr6HnA30Bv4CdgjIo+42kWkAdATmKqquzFKoVjPCoqlIsA8kfTxs+9WzNPJVA9tL2CeJEPFfqBOAfqfjXlSdbHFvVFEOorIjyKyV0QOA7cBuZcddrv9nerhfSVnrDNFZLqIbBeRI8B7bmOdDWxX51HIwV2uhsD97jcvoL5znDdyH98x1/FDgbM89VfVY5gntrOdLcf34ryv6+VcnvB1/loYJej1d/CDXW5/n8D5zr2wQ1Wr5dqO5+rzMXAB5mY1OfcAktPQ3MDHuVzf0QHndRjwnaruc95PJbA3uIJ8DwU5viEwKNfvdx7e/9cyMTNDd2IwMyAg29h9EWYGcxvwlIj0dpqvB9ao6nLn/RTgWhHJPWaxoVgqAlWdy+mLF8he4/5GRJaIyM8i0sLpu9l5Ws7yMM4PmCWHUDEL83TpLzsxN1QXuf+ppwIzgfpqDM9vAHnWXf3kWcyUOUlVq2Ceslxj7QTq5lrTdZdrG/DvXDevCs7TuTdyK5Wfch1fSVVv93Q+EamEeYrd4WwNyUkDzPKNp3N5eu/r/HsxU39fv0NIUdUTwNfA7XhQBJrT0LzVx1D9MUt6a0WkPGYG2l1EdonILuBeIFlyeSmFgIKmRN4GTM71+1VU1ee89N+KWS5zJx4PCt6ZAX8ArMQsC4KxozRy+57GYh6aLimg3BFDsVQEXhgP3K2q7YAHgP+FWR5PjAa6iMgLInIWgIg0EZH3xM0w68b7wI0ikiAiFZzj3akMHFDVkyLSgaJN4ysDx4BDIlIXsz7q4hfMU9RdIlJGRK4AOri1vwnc5sxQREQqijFk++sJ9QXQTESuF5EYZztXRFq69ekrIueJyBnA08AiVd0GfOUce60j29VAgjOmN3Zj1sfzPb8aN9CPgSdEpIKIJBAZywCjMPadzQU9UERqi8hdmOvpH6qaBVyJ+Y0TMEsubYCWmOXDG7yNFSR2A3Hiv4fXe8DlItJbRKJFpJwYV9x6XvrPAB4VkXoiEiUiF2GWEj+EbEP3pSJS2Wm/BGNjWiTGlbcx5vp3fU+JBH72FFJKhCJwnhC7AB+IyHJgHAVbggkJarwQOmOeRlY7yzkfYdYz88xMVPVr4CVgNsaINztXlzswU9ajwOMYxVFYnsQYSQ8DX2Jufi45TmEMxDdjPCWuw9w805z2xRg7wasYQ+0GzHKcX6jqUYwxewjmCX8X8DzGMOxiKubGdQBjyBzqHLsfs8Z/P2bp7SHgMrflDU+8jFn/Pygir/hx/rswyxC7MMuS7/j72QrB2ZI3jiDPLFJVd6iqx6A2HxwSkePA75j180GqOsFpG4ZZg9+qxktml6ruwvymQyWEbr7AB87rfhFZml9n54HgCoxy3IuZITyI9/vbU8ACYB7meh0DDHXsPwBHnLG2Yq73McDtzvc9DPjMsc25f08vA5eJSI0Cf9oIQHIu+xYfxPg7f6GqiSJSBVirql5v/iIy0en/Ya79PYAHVPWy4Elb8hCRRcAbqhrMm6LrXBMx3jSPBvtcFktppETMCFT1CLBJRAZBduBQqNc1SzQi0l1EznKWX4Zh3FW/CbdcFoul6BRLRSAi0zDr1s3FhH/fjFkquFlMAMpqzFQRZ603BRPBO05EVruN8zNmGnqhM07v3OeyZNMcWIFZOrofGKhOUJzFYineFNulIYvFYrEEhmI5I7BYLBZL4AilJ0BAiI2N1bi4uHCLYbFYLMWKJUuW7FPVWp7aip0iiIuLY/HixeEWw2KxWIoVIuI1It4uDVksFkspxyoCi8ViKeVYRWCxWCylnGJnI7BYihPp6emkpKRw8uTJcItiKSWUK1eOevXqERPjfzJUqwgsliCSkpJC5cqViYuLQ/IW5LJYAoqqsn//flJSUoiPj/f7OLs0ZLEEkZMnT1KzZk2rBCwhQUSoWbNmgWegVhFYLEHGKgFLKCnM9WYVgcVisZRyrCKweGb/fli3DjIywi2JpQhs3ryZxMTEHPueeOIJXnzxRQAWLlxIx44dadOmDS1btuSJJ54AYOLEidSqVYtzzjmHpk2b0rt3bxYsWBBq8f3mmWeeCdhYcXFx7Nvnq5xFycMqAotnPvwQmjeHnTbBaElm2LBhjB8/nuXLl7Nq1SoGDx6c3Xb11VezbNky1q9fzyOPPMKAAQNYs2ZNQM+fketBI/d7fwmkIiiNWEVg8cyKFeb1gw9897MUa/bs2UOdOqaeU3R0NAkJCR779ezZkxEjRjB+/Pg8bbt376Z///4kJyeTnJzMggUL8sxEXnzxxezZRo8ePRg1ahTdu3fn5Zdf5sYbb+S+++6jZ8+ePPzwwxw/fpybbrqJc889l3POOYfPPvsMMLOUAQMG0KdPH5o2bcpDDz0EwCOPPEJqaipt2rRh6NChOWR7/fXXs/u5xrj77rsBuPLKK2nXrh2tWrXy+Ll8fYa//vqLPn360K5dO84//3z+/PNPn99zpGPdRy2eOXTIvB48GF45ShD33APLlwd2zDZt4KWXCn/8vffeS/PmzenRowd9+vRh2LBhlCtXzmPftm3bMm7cuDz7R44cSffu3fnkk0/IzMzk2LFjHMznujl06BA//fQTADfeeCPr1q1j1qxZREdHM2rUKC644AImTJjAoUOH6NChAxdddBEAy5cvZ9myZZQtW5bmzZtz991389xzz/Hqq6+y3MOXO3DgQDp37syYMWMAmDFjBv/85z8BmDBhAjVq1CA1NZVzzz2Xq666ipo1a/r1vY0YMYI33niDpk2bsmjRIu644w5mz85dSbb4ELQZgYjUF5EfRWSNiKwWkb976NNDRA6LyHJnezxY8lgKyLZt5nX//vDKYSkS3jxIXPsff/xxFi9eTK9evZg6dSp9+vTxOpa32iWzZ8/m9ttvB8ysomrVqvnKdfXVV+d4P2jQIKKjowH47rvveO6552jTpg09evTg5MmTbN26FYALL7yQqlWrUq5cORISEtiyxWseNQBq1apFo0aNWLhwIfv372ft2rV07doVgFdeeYXk5GQ6derEtm3bWL9+fb5yAxw7dowFCxYwaNAg2rRpw6233srOYr6EGswZQQZwv6ouFZHKwBIR+V5V/8jV72dbLzgCcf7xKGVGs2BSlCf3wlKzZs08T+cHDhzIEWzUuHFjbr/9dv72t79Rq1Yt9ntR/suWLaNly5Z+nbdMmTJkZWVlv8/t116xYkWv71WVjz76iObNm+fos2jRIsqWLZv9Pjo62i+bwtVXX837779PixYt6N+/PyLCnDlzmDVrFr/88gsVKlTIVjj+fIasrCyqVavmcQZSXAnajEBVd6rqUufvo8AaoG6wzmcJMK+8AiJWERRzKlWqRJ06dfjhhx8AowS++eYbzjvvPAC+/PLL7Cf99evXEx0dTbVq1fKM89NPPzF+/Hj+9re/5Wm78MILef311wHIzMzkyJEj1K5dmz179rB//37S0tL44osv/Ja5d+/e/Pe//82Wa9myZfkeExMTQ3p6use2AQMG8OmnnzJt2rTsmcjhw4epXr06FSpU4M8//2ThwoV5jvP2GapUqUJ8fDwfOPYzVWWFy6ZWTAmJsVhE4oBzgEUemjuLyAoR+VpEWnk5foSILBaRxXv37g2ipJZsrrjCKIMrrgi3JJYiMmnSJP71r3/Rpk0bLrjgAkaPHk3jxo0BmDx5Ms2bN6dNmzZcf/31TJkyJXuJZsaMGbRp04ZmzZrxzDPP8NFHH3mcEbz88sv8+OOPtG7dmnbt2rF69WpiYmJ4/PHH6dixI5dddhktWrTwW97HHnuM9PR0kpKSSExM5LHHHsv3mBEjRpCUlJTHWAxQvXr17GWkDh06ANCnTx8yMjJISkriscceo1OnTnmO8/UZpkyZwttvv01ycjKtWrXKNmgXV4Jes1hEKgE/Af9W1Y9ztVUBslT1mIj0BV5W1aa+xmvfvr3awjRBZudOWL0aOneGXFN4S8FYs2aN38spFkug8HTdicgSVW3vqX9QZwQiEgN8BEzJrQQAVPWIqh5z/v4KiBGR2GDKZPGDWbPg4ouNwfjQIQjyw4LFYgkvwfQaEuBtYI2qjvXS5yynHyLSwZHHuqmEG5cnxvvvQ/XqcPx4eOWxWCxBJZheQ12B64HfRcRlXh8FNABQ1TeAgcDtIpIBpAJDNNhrVZb82boVatWC+vXN+337oFKl8MpksViCRtAUgarOA3ymwVPVV4FXgyWDpZBs2QING4IruGbfPoiLC6tIFosleNgUE5a8uBRBrGOusS6kFkuJxqaYsORlyhQoUwbKlzfvbXSxxVKisTMCS17atYPkZKhTBx5/HFp5DO+wFANKQxrqOXPmcNllJjnBzJkzee6554o8Zo8ePfDkpv7SSy9x4sSJAo01dOhQmjdvTmJiIjfddFN24JuqMnLkSJo0aUJSUhJLly7NPqZPnz5Uq1Yt+3O5mD17Nm3btiUxMZFhw4YVOltrbqwisOQkJQUmToS9e6FyZXjySZPZzFIiCXca6szMzICO169fPx555JGAjulOYRXBn3/+ye+//05qaipvvfUWAF9//TXr169n/fr1jB8/PjtfE8CDDz7I5MmTc4yTlZXFsGHDmD59OqtWraJhw4a8++67Rf9QWEVgyc2iRTB8uFEIAHv22JoEJZhApKG+8cYbGTlyJF26dKFRo0Z8+OGHgHniffDBB0lMTKR169bMmDEDME/wPXv25Nprr6V169Zs3ryZFi1acMstt5CYmMjQoUOZNWsWXbt2pWnTpvz6668A/Prrr3Tp0oVzzjmHLl26sHbt2jyyTJw4kbvuuguANm3aZG/ly5fnp59+8priOjU1lSFDhpCUlMTVV19NampqnrFfeeUVduzYQc+ePenZsycA06ZNo3Xr1iQmJvLwww97/O769u2LiCAidOjQgRTnf+uzzz7jhhtuQETo1KkThw4dyk5ed+GFF1K5cuUc4+zfv5+yZcvSrFkzAC6++GI++ugjj+csKNZGYMmJK4agYUPz2r07JCbaugSBokePvPsGD4Y77oATJ6Bv37ztN95otn37YODAnG1z5hRJnECkoQbYuXMn8+bN488//6Rfv34MHDiQjz/+mOXLl7NixQr27dvHueeeS7du3QBzU1+1ahXx8fFs3ryZDRs28MEHHzB+/HjOPfdcpk6dyrx585g5cybPPPMMn376KS1atGDu3LmUKVOGWbNmMWrUKJ83QldSuM8//5wxY8bQpUsXRo8e7THF9bhx46hQoQIrV65k5cqVtG3bNs94I0eOZOzYsfz444/ExsayY8cOHn74YZYsWUL16tXp1asXn376KVdeeaVHedLT05k8eTIvv/wyANu3b6e+y0UbqFevHtu3b89WzLmJjY0lPT2dxYsX0759ez788EO2ubIEFxE7I7DkZMsWEzNQvbp5HxtrvYaKMaFIQw2myEtUVBQJCQns3r0bgHnz5nHNNdcQHR1N7dq16d69O7/99hsAHTp0yJEBNT4+ntatWxMVFUWrVq248MILEZHsGQOYRHGDBg0iMTGRe++9l9WrV+f7+devX8+DDz7IjBkziImJ8Zrieu7cuVx33XUAJCUlkZSUlO/Yv/32Gz169KBWrVqUKVOGoUOHMnfuXK/977jjDrp168b5558PeP4+fRWeFxGmT5/OvffeS4cOHahcuTJlygTmWd7OCCw52brVzAZcF2RsLPz1V3hlKkn4eoKvUMF3e2xsgWcAoUpD7Z4e2nWD86U4cqehdj8+Kioq+31UVFS2QfSxxx6jZ8+efPLJJ2zevJkenmZXbhw/fpzBgwfz5ptvcvbZZ2fL5CnFNfi+CXuiILGvTz75JHv37s0xo6pXr16OJ/qUlJRsOb3RuXNnfv75Z8DUbVi3bl2BZPaGnRFYcuKKIXBRs6adERRjQpGG2hvdunVjxowZZGZmsnfvXubOnZud/bMwHD58mLp1TSb7iRMn5tt/+PDhDB8+PPsJHLynuO7WrRtTpkwBYNWqVaxcudLjmJUrV+bo0aMAdOzYkZ9++ol9+/aRmZnJtGnT6N69e55j3nrrLb799lumTZtGVNTpW26/fv2YNGkSqsrChQupWrWq12UhF3v27AEgLS2N559/nttuuy3f78Ef7IzAkpOvvwZ3Q5lraUj19CzBUqyYNGkSd955J/fffz9AnjTU9957LxUqVKBMmTJ50lDPmzePEydOEB8f7zUNtTf69+/PL7/8QnJyMiLCmDFjOOusswpd3/ehhx5i2LBhjB07lgsuuMBn3y1btvDhhx+ybt06JkyYAJgb8mOPPcY999xDUlISqkpcXBxffPEFt99+O8OHDycpKYk2bdp4VVgjRozgkksuoU6dOvz44488++yz9OzZE1Wlb9++XOEhbfttt91Gw4YN6dy5M2DqIzz++OP07duXr776iiZNmlChQgXeeeed7GNcdZCPHTtGvXr1ePvtt+nduzcvvPACX3zxBVlZWdx+++35fg/+EvQ01IHGpqEOMQsXwsqVcPPN4NwgLP5j01BbwkFEpaG2FDO2b4ennoJNm07v69QJRoywSsBiKcFYRWA5zcqVMHo07Nhxet+JE7BsGRw+HD65LBZLULGKwHKa3DEEACtWQNu28Msv4ZGpBFDcll8txZvCXG9WEVhOs3WrSTbn7rngnoraUmDKlSvH/v37rTKwhARVZf/+/V6DAr1hvYYsp9myxRSjcbcH2FTURaJevXqkpKSwd+/ecItiKSWUK1eOevXqFegYqwgsp0lJybksBFCtGkRF2VTUhSQmJiZH8JbFEolYRWA5zezZcOxYzn1RUTaozGIp4VhFYDlNdDRUrZp3/xtvQIMGoZfHYrGEBGssthh27jTxAk7GxhwMGADtPcahWCyWEoBVBBbDunXw5puel4DWroWffgq9TBaLJSRYRWAxeIohcPHCC3DNNaGVx2KxhAyrCCwGlyJwK5SRjXviOYvFUuKwisBi2LIFatcGT4EosbGQnp7Xo8hisZQIrCKwGNLSoEkTz202uthiKdFY91GLYfJk70s/7tHFNjjKYilx2BmB5TTeCs906gTffgvNmoVWHovFEhKsIrCYJ/1LLoEff/TcXqsW9OrlOdjMYrEUe6wisJhCNN98A04t1jxkZsJHH8Hvv4dWLovFEhKCpghEpL6I/Cgia0RktYj83UMfEZFXRGSDiKwUkbbBksfiA18xBGCWjAYPhvffD51MFoslZATTWJwB3K+qS0WkMrBERL5X1T/c+lwCNHW2jsDrzqsllOSnCKKioEYN6zVksZRQgjYjUNWdqrrU+fsosAaom6vbFcAkNSwEqolIHSyhZcsWqFLFpJz2hiuozGKxlDhCYiMQkTjgHGBRrqa6wDa39ynkVRaIyAgRWSwii22BjyBQqRJ07uy7T2ysrUlgsZRQgh5HICKVgI+Ae1T1SO5mD4fkcWZX1fHAeID27dvbPAeB5pln8u9TsyZs3Bh8WSwWS8gJqiIQkRiMEpiiqh976JICuCe3qQfsCKZMlkIyZozNNWSxlFCC6TUkwNvAGlUd66XbTOAGx3uoE3BYVXcGSyaLB44ehcRE+OAD3/2aNYPmzUMjk8ViCSnBnBF0Ba4HfhcRV7WTUUADAFV9A/gK6AtsAE4Aw4Moj8UTW7bA6tUmVsAXa9bArFlw881QoUJoZLNYLCHBqyIQEU9LObk5oKq3eGpQ1Xl4tgG491HgTj/OYwkW+bmOuli0CEaOhEsvhUaNgi+XxWIJGb5mBK2B23y0C/ByYMWxhJytW81rforAlXhu/36rCCyWEoYvRTBaVX/wdbCI/DvA8lhCzZYtEBMDZ53lu597BlKLxVKi8GosVtWp+R3sTx9LhFO/Plx5pYke9oWtSWCxlFh8/veLyFAR+VVEDjvbQhG5NlTCWULAnXf6l0PIfWnIYrGUKHwZi68DHgLuB5ZibAJtgTEiIqo6JTQiWiKCatVg/XqoYzOAWCwlDV8zgjuB/qr6varuV9V9qvodMAC4OzTiWYLKqVMmmdxrr+XfV8SUsqxYMfhyWSyWkOJLEVRV1b9y71TVjYCtUFISSEmBgwehfHn/+k+cCFPsRNBiKWn48hpK9dF2ItCCWMKAvzEELiZMgOhoGDo0eDJZLJaQ40sRtBSRpR72C2CL15YECqoIataEDRuCJ4/FYgkL+QWUWUoyrmCyevX86x8bCwsXBk8ei8USFrwqgtz2ARGpDpwHbFXVFcEWzBICWrWCESOgXDn/+rtqEqga47HFYikReDUWi8inIpLo/H0WsBq4A5ghItZrqCRw1VUwbpz//WvWhPR070XuLRZLscSX11BTVV3l/D0c+EFVLwE6AH8LumSW4HP8eMH633GHOaZKleDIY7FYwoIvRZDu9veFmJTROFXGsoIplCUEZGWZJ/zHHvP/mAoVbApqi6UE4ksRbBeR20XkcqAd8A2AiJQDzgiFcJYgsmcPpKXln2zOne3b4b77YIU1EVksJQlfiuBmjAK4DbhWVQ86+7sA7wZbMEuQcbmONmjg/zEnTsD//R+sXBkcmSwWS1jw5TW0C8hTdEZVZwOzgymUJQQUNIYAbOI5i6WE4ivp3CeAe7VyBfYBP6rq9GALZgkyhVEEVauadNU2FbXFUqLwFVD2qod9NYChIpKkqqOCJJMlFHTsCKNGmZu7v0RFGQOzVQQWS4nC19KQx+pkzkxhMaYQvaW40q2b2QpKrVrGVmCxFFdUYe9ecy3bwEjA94zAI6qaIfbLK/789ZfxGCpoWumVK03iOYuluNK1K/zyC+zaBbVrh1uaiMBXZHEVD1tDEXkMWBNCGS3BoF07ePjhgh9nlYClOKNqlADAunXhlSWC8OU+uhpY5by6/v4EqI5xKbUUVw4fNltBDMUupk6FW28NvEwWSyg4ePD03+vXh0+OCMOXjaB+KAWxhBBX1tHCKIKVK+Gdd+CNN+z6qqX4sWnT6b/tjCAbX0tDSfkd7E8fSwTiUgQFCSZzERtrE89Zii/uisDOCLLxZSyeLCLnYQrReONd4JzAimQJOn7GEGRmeh0qpj4AACAASURBVDAJuAeV2eRzluKGSxG8/HLhZsQlFF+KoCbGNuBLEewJrDiWkNC9O7zyik+PidWroW1bWLIEEhPdGlyKYN8+iI8PrpwWS6C56y64/HJo0SLckkQUvmwEfpatshQ7WrUymw/mzIFTpzwogjPPNFuqr5LWFkuEUr68UQKpqbBmDTRtCpUrh1uqsFPgOAJLCWDhQqhTx+fUeNky8+q+pApAhw6we3fwZLNYgsmzz0KbNqYq3wUXwKxZcOGF4ZYq7PhyHy0SIjJBRPaIyCov7T1E5LCILHe2x4MliyUXAwbAk0/67OJSBJs3B18ciyUkqMJTT8EPP5iZAFjPIYegKQJgItAnnz4/q2obZ3sqiLJYXKSlwc6dPj2G0tNhlaO+PSqCQYPgzTeDIp7FEjR27YKTJ/l8dTzdhpxtiixZzyHAT0UgIkNE5J/O3/VFpF1+x6jqXOBAEeWzBJqUFPPqY1nojz+MfaBiRQ9LQwBz58LixcGRz2IJFs7FPD8ljvm/RJHVuImdETjkqwhE5FWgJ3Cds+s48EaAzt9ZRFaIyNci4tV6KSIjRGSxiCzeu3dvgE5dSvHDdXT5cvPat6/RG+npuTrExtqaBJbihzO9XbQnnqwsOFK7mZ0ROPgzI+iiqrcCJwFU9QCBKVW5FGioqsnAf4FPvXVU1fGq2l5V29eqVSsApy7F+BFVvGyZmTVffLEpbbxtW64ONhW1pTiyYwcAi/bEAfBr9wfg9dfDKFDk4I8iSBeRKJwiNSJSkwAUr1fVI6p6zPn7KyBGRGKLOq4lH3r1go8/hvreM4gsWwZJSdC4sXmfx04QG2sVgaX48cAD/LHoKKlUAGDeqY7Gc8jilyJ4DfgIqCUiTwLzgOeLemIROUucfNYi0sGRxa43BJuzz4b+/eEMz5O6rCyzNHTOOafjxfIogubNzTgWSzFj3Y5KgImY37DyBHz2GWzcGGapwk++cQSqOklElgAXYaKMB6mqR5dQd0RkGtADiBWRFGA0EOOM+QYwELhdRDKAVGCIqqqX4SyB4quvzBN9hw4emzdtgiNHjCKoV88UJctjMH722eDLabEEmjvvpMyh84BrOO882PLHcbjySnjpJfj738MtXVjJVxGISBVgG/CO275oVc30dZyqXpNP+6t4LodpCSYjR0L79jDdc9lpV/zAOedATIxZQbKxBJZiT0YGjBtHdOuq1KhxDV27wpjnY9GqVRHrOeTX0tAfwEFgq7MdBLaJyK8iYhPOFSdclt98DMXR0afTSsTFeZgRzJplZhQuw7PFEumkpEBmJmtOxtOkCbRsCRmZwskG1nMI/FMEnwH9VLWaqlYDLgemAPcSODdSSyjYvdsECOTjOpqQYCLwwSiCPDOCkyfht99MgI7FUhxwLuIlB4wiSEgwu/dWbWpjCfBPEXRU1S9dbxwPn56qOh8oFzTJLIHHFUPgI6p42TKzLOQiPt543aWluXVyz0BqsRQHnGntoj1GEbRoYeoqbYppZma2J0+GWcDw4o8iOCQi94tIXWe7z9kXDfi0E1gijHyCyXbvNtkn2rQ5vS8uzqRoybEK5F6TwGIpDqSmklE9li00oEkTEycTFwczKt9iQum9eNGVFvxRBNcATYBvgG+BpsBQjKHZp0HYEmFceqlJDdGsmcdmd0OxC48upDVrmlc7I7AUF+64g28m7SWDGJo0MbsSEmDBlrpmehAVzLRrkY8/7qN7gdu9NK8NrDiWoFKpErTznibKpQhyzwggl8G4alXo0uW0QrBYigEbNphXlyJo2RJmfa9kvfIaUa1alup01P64j8YC9wOtcLMJqGqvIMplCQbvvWeKcFxxhcfmZcvMDKBatdP76taFMmVyzQiiomD+/KCKarEElH79OOtIX6pUuS17ZTMhAdJOCVlPPU1U/36lWhH4Mx96D9gMNMNEFO8ClgdRJkuwGDMG3n7ba3NuQzEYV9IGDWwsgaUYk5YGX3xBespumjQxRmI47Tl0+EzrOeSPIqilquOAU6r6AzAM8ByWaolstmzxaig+etRMnXMrAvASSzB8OAweHHARLZaAs2ULqLLiaHz2shCYpSGA7RVsLIFfSeec110i0htIBLxnLLNEJocOmdwRXhTBihXm1ZMiiI/3MCM4fNjUfLVYIh3nKWbxvpyKoEoVs/S5NqupcZc7ejRMAoYffxTBMyJSFXgAeBSzVPRgUKWyBJ580k97MhS7iIszsWM56tXbDKSW4oKjCDZkxWdXqHSRkAC/HWmWo19pJF9FoKozVfWwqq5U1fNVNVlVPw6FcJYA4ociqFXLc1JRl+eQKwwBOF2TwOYJtEQ65ctzuPE57ODsHDMCMIpgws6+ZB06YnKvl1L88RpqANwFxLn3V9UBwRPLEnAuvdREjLm7BLnhMhS7DGnuuMcStGjh7IyNNYm8jh41c2yLJVIZNozJR4ehd+NREew/UZ5th6Bh1fCIFwnkqwiAmcAk4HsCUJDGEiZE4MwzPTadOgWrV0Pv3p4P9RhLkJgIV11lDrZYIpwNG0wN7tq1c+53eQ6lPjkG2pQz2XlLIf4oglOqOjbokliCy+uvm+yjd96Zp2n1alOX2JOhGKBOHROBn8Ng3Lu3d81hsUQSyckknbqOJk0ezDPjdXkOVfj5W1hzotQqAn+Mxf8VkUdF5FwRSXJtQZfMEljefRc++cRjk6tYvTdFEBVlTAs2lsBS7Dh2DFau5MC+rDzLQmBMXWeeCRvLlO5YAn8UQTPgTuAlTNnK17AFZYofPmIIli0z02ZP/ygu8sQSbNli7ARTpgRUTIsloDgX7dKD8V6v74QEWJnaDA4cKLWJFP1RBIOBOFXt6ngNna+q3YItmCWAnDxp/D99KILkZN95t/LEElStav5p9uwJqKgWS0BxuY5mxvlUBAv2On6lpTSwzB9FsBKoHGxBLEEkJcW8elAE7sXqfREXB3v3mpk2YBRBdLSNJbBENo4i2IT3GUHLlrDsRDOyqlYrtTMCf4zFNYE/RWQRkF2exLqPFiN27zbWXg+K4K+/zM09P0XgciHdsgVatcJ4IdWsWWr/cSzFhHr12JR8JftWxPqcEayjGT9+eIALL/LgP10K8EcR/DvoUliCS9euJizYQ/CXpxoEnnC5kG7e7CgCsNHFlsjnqqt47ZerKLfWc7AkuFxIhT/WwIUXhVK4yMGfegQ/hEIQS5DxYgBYtsykmc6+uXvBYyzB1Vd7DVCzWCKCzEw2bIimcWPvNrDataF6dag79QX47XeYNCm0MkYA/kQWnwv8F2gJlAUESFNVG05aXPjPf4xR9/nn8zQtX26UQNmyvoeoXdsUtM9hMH788YCKabEEFFWIjaV3zEjo8qTXbiJmVpC1KQVWfWKO8xRiX4Lxx1j8P0zq6Y0Yo/FdGFdSS3Hhyy9h3jyPTcuWeU40lxsRL+mo09M9dbdYws/+/XDoEBsPVvfpGg3GYLz4SFNjMNu1KzTyRRD+KIIoVV0LlFHVdFV9EyilK2nFFC8xBDt3GjtyfvYBF3lcSEePNlXAbeI5SyTiXKzrMrx7DLlISIClx5wspKXQhdQfRXBcRM4AVojIMyJyN1ApyHJZAkVWFmzbZsqM5cJfQ7GLuDgPsQQZGabOgcUSafjhOuoiIQHW48QSlMIIY38UwY1Ov7uATKApMDCIMlkCyc6dZvnGw4zAVw0CT8THm+DL7Pu+q3i99RyyRCKOItiM92AyFwkJsJUG7K3f1rhalzL8qUewESgPiKo+pqojVbX0qcziyqFD0Lix2XKxbJnZ7W8WaXcXUoDsKuBWEVgikTZt+LndPZyMqUL9fGoq1qsHFSpF8/SVS+CGG0IjXwThVRGI4VER2Y0pXr9VRHaJyKiQSWcpOq1amRy8vXrlafJUrN4XeVxIXYrABpVZIpFevXip4f/RqJEJgveFiDEYl9bqq75mBCOBC4DzVLWqqlYGugEXiEi+uVpFZIKI7BGRVV7aRUReEZENIrJSRNoW6hNYCsXhw7BxY8EUgXuBGsBohvvvJ9/HLYslHOzcycb1mfkuC7lISIBzf/ufudCzSlfpFV+KYBhwtapmm9CdJaFrMXaD/JgI9PHRfgnG3tAUGAG87seYloLy9NMweHCe3a5i9f7aB8CYBCpWdFMEtWvDiy9C69ZFFtNiCShZWWhcHDes+UeBFMHewzHmAt+2LajiRRq+FMEZqro3905V3YMJLPOJqs4FDvjocgUwSQ0LgWoiUie/cS0FZNEij+5wBfUYAjN9jo/PFUtw7JixQ1gskcTOncipU6zPyN9Q7KI0ew75UgS+ahCm+Wjzl7qAu9pNcfZZAomXGIJly8wDfZ0Cqt48LqSNG8PDDxdFQosl8BTAddRFy5Ym+RxQ6mIJfCmCZBE54GE7CARiLcBTDLfHyCQRGSEii0Vk8d69eSYpFm+o+lQEBZkNuHBFF2fHkNnEc5ZIxHlaKYgiiIuDg2XrkBZT0c4I3DgDqOVhiwXKBeDcKYC7lbEesMNTR1Udr6rtVbV9rVq1AnDqUsKhQ3D0aB5FkJYGf/xROEUQH2/iCLJXg2JjrdeQJfJwZgTbouK81WPKQ3Q0tGgpfFfnxtNV7UsJXhWBqmb62gJw7pnADY73UCfgsKruDMC4FhfHjkH37nku6lWrTEBwYWcEkCuWwM4ILJHGxRfzbvJY6sSXIybG/8MSEuBueRVGjAiebBGIP/UICoWITAN6ALEikgKMBmIAVPUN4CugL7ABOAEMD5YspZb69WHOnDy78ytW7wuXC+mmTc7xNWtaRWCJPDp14pXoTn4vC7lISICpU+HY/jQqVY/xXb+1BBE0RaCq1+TTrsCdwTq/Ba/pdJctg8qVoVGjgg+ZZ0YwcKB1H7VEHLroVw6sa0TnzrEFOq5lS7iKD6l45tUmuqxZsyBJGFmUDnVXWrn5ZuiTN5TDn2L13qhWzaSkyHYh7dUL7r67aHJaLIEkPR26dObmYy8VakawnbpIVlapMhj7SjFx0JvXkIj4ig+wRAo//WTSRLuRmWmCyQqzLASnYwmyZwSpqbB2rXm1WCKBlBQkK6tAHkMuGjeGzWWcWIJS5ELq65kwFu9eQ9Z1J9LZtcvkkOjaNcfuDRvg+PHCKwLIFUsweza0aAG//174AS2WQFKIGAIXMTFQs1lNjsZUtzMCyOs1BFQFarttlkhm/nzzmksRFCaiODeu6GJVbAZSS+SRnX46Ptu5oSC0TBA2RjW1MwJ3RORSEVmH8ftf5LzODrZgliIyf74pMtw2Zy6/ZcvMU09R3KTj4sysYv9+bE0CS+SxaROZEo00qJ9vLW5PJCTAf0/dSvpVQwIvW4Tij9fQv4GuwHeqeo6IXAxcFVyxLEWmSxfjGpSryMby5ZCYWLTaG+7pqGOb2hmBJcIYOpRHZ7QhvmHhnCITEuApvYmRXSEpwKJFKv74jWQ4yeeiRERU9XvApoyOdAYOhCefzLFL1f9i9b7IkY66alUTkmmjiy2RQsuWjD8wsMD2ARcJCRBFJlt+2mwi80sB/iiCwyJSEZgHTBKR/wClK1l3cWPPHo9pdHfsgL17i2YfgFyxBCIwbhz071+0QS2WAHH03Y+pfmBDoRVBs2bQVpZz+ch4+P77wAoXofijCK4ETgL3AHOA7cBlQZTJUlTeftsUqz+Q08s3EIZiMHEENWq4xRLcfDO0b1+0QS2WQJCaSuUbr+IaptG0aeGGKFsWMuJLVzpqfxTBPxzPoXRVfVtVxwL3BVswSxFYsMC4dNaokWP3smXmAT45ueinyOFC+tdfsHRp0Qe1WIrKli1A4VxH3WnYugr7omuXGs8hfxSBpypjlwZaEEuAyMoyiiCX2ygYRdCkibEhFxVXOmoAHnqoVBb8tkQgbjEEhUmh4iIhAdZkNiNrbSmfEYjIrSKyDGguIkvdtvXAH6ET0VIg1q41S0JeFEFRl4VcuKKLVbGJ5yyRg6MI0urEU7584YdJSIB1NCXrz9IxI/DlX/U+8APwLPCI2/6jTrlKSyTiJZDs0CFz4w5Udt24ODh50tila7tqEnhJcmexhIxNm0iTslRpdlaRhmnZEu7mZhJv6EXHUnBd+4osPqiqG1R1EFAeuNjZbHqJSObSS00e3VyWsqKknvaEezpqYmNNgYMjRwIzuMVSWO6/n/7V5tC4adHyabZoAb/Qhe+qX13ilQD4F1l8J2Z20MDZ3heRO4ItmKWQ1KkD11yT5+INlMeQixwupDbNhCVCOFLhLL4+WPA6BLmpWBGaNEwn66efTc6uEo4/avNWoIOqjlLVUUBH4LbgimUpFPv3w//+BzvzFnpbtszoiNoByhLlKv+3aRPQowd88gmceWZgBrdYCsmRp1+mLUuKrAgAEltkMPqHbvDee0UfLMLxRxEIkO72Ph3Phect4ebnn+HOOz0+wQTSUAxQqRLUquXMCBo0gCuvDIw7ksVSWI4cod6L93AhPwREETRpXZ6tNCgVnkO+vIZchuTJwEIReVREHgUWAO+GQjhLAZk/3yQRatcux+7UVFNsKZCKANxiCdLS4OuvS8UU2hLBuLmONm5c9OFatjSeQ6dWlXzPIV8zgl8BVHUMMAJTVzgVuE1VXwyBbJaCMn++ifAtVy7H7lWrTEGaQCsCVzpqTp6Evn3h008DewKLf/zwg4njyCrlmV8cRXC0ZjyVKhV9OONC2oyov9Y5ftIlF1+KIHv5R1V/U9WxqvofVf0tBHJZCsrJk7Bkicf4AZfHUFGTzeUmLs4EcmZVqgJlylhjcThQhe3bYfLk0z90acVRBNFNClGEwAMtW8J6mnLG8UMlPqmirziCWiLiNZWEk2rCEimsXm1qtXoJJKtShUIV6fBFXBycOgU7dwl1bVBZeHjoIZgzx/z91Vd56k+UKjZt4phU4swWNfLv6wdVq8K82gN5+pyOPFalSkDGjFR8zQiigUpAZS+bJZJo185EFPfqlafJlXq6MMXqfZEjHbUrqMwSOrKyTMxI3bpw7rlGEZRijj/5Igm6miZNA+fLUj2pPjP3dSlaAY9igK8ZwU5VfSpkkliKTrVqeXZlZsLKlfC3vwX+dO6xBF1jY+2MINTMm2dyiw8ZAn/+CU89ZX4DV1xHKWNjyhlso0FAPIZctGwJB8Z9QNasmkRddEHgBo4w/LIRWCIcVRg8GD7/PE/TunVw4kTgDcWQK5Zg7Fh49dXAn8TinenToUIFuPxyY6xv08YohtKIKuX/eR/dmRNQRZCQAI+njSL1pTcCN2gE4mtGcGHIpLAUjXXr4IMP4OKL8zQFOqLYnfLl4ayznKWh0rw2HQ4yMsxvfvnlJgy2Q4fSnQp83z6afP5/JNOQxo17BGxYl+dQ7TUlO5bAV66hA97aLBGGl0RzYBTBGWeYKW4wyE5H/ccfpiBOCXezixhOnYJHHoHbcgX5nzxZOt1IHY+h/ZXjPa2QFpqEBOM5VD5lfYm+tgNsPrSEhQULoHp1kykrF65i9TExwTm1Kx0133wDt9wChw8H50SWnFSoAPffb9J7uJgzxxQj+q0Ueni7imME2DWuZk3YWakZMadOlOhlN6sISgLz50OXLnncglzF6oOxLOQiLg62boWsGhGSeG72bBgzJrwyBJu0NJgyJW9h9datTVtp9B5yFEH5lnEBH7o0lK20iqC4c+qUySR3YV6TTkqK8egMpiKIjzfL1ftwFEG4XUgfeQQeftitjmYJ5Jtv4LrrzEzQnZo1oVOnUqkI0nfvZw+1qNcy8J7tGR270rrKFrRb94CPHSlYRVDcOeMMsyRw7715moJpKHbhciHdfrKm+SOcM4L1608vi0ydGj45gs306cZF9AIP7ox9+8LixbB7d+jlCiPr//YCddkeUI8hF02TK7DqSAN27Sm5t8ugfjIR6SMia0Vkg4g84qH9RhHZKyLLne2WYMpTIvFhGJw929gGkpKCd3qXIth0NAKWhk6cgEsugfPOK7mGvePHYeZMGDgw2/AzdSq8+abT3revef322/DIFyY2bIAMYoKiCBIS4Abe5cgzJdg9WlWDsmEik/8CGgFnACuAhFx9bgReLci47dq1U4sbPXuq3nRTnt0nTqhWr646eHBwT3/ypKqI6lOPpqkuX656+HBwT1jamT5dFVTnzFFV83VXqaJaoYLz1WdlqT7zjOqff4ZXzlCSmalrk67SS/lc9+0L/PA7dqh+zJW6v3bLwA8eQoDF6uW+GswZQQdgg6puVNVTwHTgiiCer/SRlmbWiatWzdP04Ydw8GDgahR7o2xZOPtsE9VJcrJJahQOtm83mwtV2LYtPLIEk/nzzRd+3nmA8dg9csRMht5/H1OZ7h//gObNwytnKNmxg2YrP6J5+W3UCEyaoRycdRZsLduMynv/MqH6JZBgKoK6gPt/YoqzLzdXichKEflQROp7GkhERojIYhFZvHfv3mDIWjxZutQoAw/xA+PGQZMm0LNn8MXITkc9aRJ88UXwT+iJMWNMneYTJ8z7O+80+XcyMsIjT7B45RVj/ImOJiMDXnoJzj/fxIlMmOD0OXUKvvyyRHu55MDxGEqvFx+U8sIikFqvKTFZp4yLXAkkmIrA00+Se+H2cyBOVZOAWXgpeKOq41W1vaq2r1WrVoDFLMa4vEZyKYLVq82D44gRgU8054nsAjVjxsA77wT/hLnJyDAG1L59jX89mCjr3bth1qzQyxNsnJKgH35o7ksPPgjDh8Mvv5iUQ6SmmopxEyeGVcyQ4XiIlWka4PS6bkS1aGb+KKHKNZi3iRTA/Qm/HpAjIkNV96tqmvP2TSBnaS2Lb+bPh0aNzNzVjfHjjTPRjTeGRoy4OOOqqjXClIr6hx9gzx4YOvT0vr59TZDd5MmhlydYDBliXGMxK18vvmhWgC69FK6/HqKjHT1ctap5OCglbqQZG8yMoGpSw6Cdo0r7ZmQQzZF1eeuBlwSCqQh+A5qKSLyInAEMAWa6dxCROm5v+wFrgihPyaNvX7jnnhy7UlPNCs2AAaamcCiIjzdLpycqhCkD6ZQp5ubn8pgBY7wYPBg++SRv4FVxZN8++OgjXGsfc+eaOkT33WdmfWedZRTCpEnOaljfvrBiRU67SQnl0EFYTjJxLcrl37mQxHWsTQVOsDz5xqCdI5wETRGoagZwF/At5gb/vqquFpGnRKSf022kiKwWkRXASIwXkcVfbrkF7r47x67334dDh4JvJHbH5UJ6KCYMNQlOnTrtTlm2bM626683mrEklND8+GNzhx8yBDCzgVq1zEd0cdNNsGuXiTfLVorffBN6WUPMr30e5xyWB8V11EVCKyGdM1hTUh9VvbkTRepm3UcdtmxR3bUrz+7OnVWbNTNehKFi40bj0bi87z9Uo6NDe3JV49+3eXPe/VlZqjNnqqamhlaeYNCzp2rz5qpZWbpmjfm+n3giZ5dTp1TPPFO1f381n71ePdVhw8IhbUh5+WXzfXj4dwgYWVmqt5Z9Rxe2HBa8kwQZwuQ+agkm//qXWSB2Cyj7/XdjMBwxgqB4T3ijXj2zPPFlwoPmkTTU1KlzujiCOyImTXO54C0ZhISdO030+JAhIMLYseYj3XFHzm4xMWaG8PnnsGevmIsh25WohJKezhX/PpdhZae7bOhBQQTa19xI+zWTzSy0hGEVQXFl/nzo3DmHW5DLSDxsWGhFiYkxyuDP3dVN6oNQaaH9+00ksVu2zWPHcoUPZGXB008Xf6Px/ffDtdeyZ4+xAwwb5tkGNHy4WUF67z1Oa+iSzLZtNNyzmIZnpgb9ssuMb0o0WbBxY3BPFAZK+FVSQjl40OT/d3MbPXHC3OsGDgxPpcL4eEhdsxkefTR0/ygffGDWwN1ybN9wg0nCmW2qiIoyj8hjx4ZGpmBQpw688AI0a8Zrr5kHUg+ppQBo1Qo6djQTAVXg738v2dlYg5R+2hPlk40L6fHl64N+rlBjFUFx5JdfzKubIpgxw5QCCKWR2J24ODi5ZTf8+9+EzKI2ZYpJBJOcDMCvvxonocOH4fnn3fpdf70pzLBqVWjkCiQ7dxr32MxMTpyA114zq12+AoeHDzexJIsXY36Ldz2G55QIMh3X0Qqtgq8Izuxq0lHvW1DyYgmsIiiOzJ9vnMY7dMjeNX68qUvTrVt4RIqPhzV7Q5iKessWU7x96NDspajHHjOZmAcMgP/+181z8uqrzfc1ZUrw5Qo0kyfDRRfB1q1MmmS+2gce8H3IkCHGhjBhAmbp7I8/Smxa7iMrNpFBNLXaeEpaEFiadqzB7ySWyMSuVhEUR0aMMD7lFSsCsHIlLFwYeiOxO3FxsJcQZiB1pZm+5hrA+NV/950pR/DCC2ad/F//cvqeeSb07m0UQXEr4zh9OnTsSFbDeMaONVkznDRDXqla1SwRTpsGJy9w3Ei//jr4soaB3Vm1+Iq+NGnhq/x6YIiLgw7lfuf9+vcH/VyhxiqC4kjDhnDF6fx948YZF/pQG4ndiY+HI1QhK7pMaBRB/fpw880QH4+qMU3UqWM8aRo1Mkrxrbfgr7+c/sOHm8IMBw8GX7ZAsXatySs0ZAiff27KLTzwgH/K/qabzBLZx6uamS+khEYZz066hyuYGdQYAhfR0WZJ7o8/gn+uUGMVQXFjwwaTfP7QIcCkp3/vPRg0iKBkXvQXE1QmnKwYCwcOBP+E111n7vTA99/Dzz/DP/95OtXQo48aG/Lo0U7/gQPhs8/M2lFxYcYMc9cfNIgXXzT6f8AA/w7t3t38JhPeEfOE0KhRUEUNFxs2QPny5iEgFAyv9AEvzWoFJS35pbcAg0jdSn1A2dixJnomJUVVVd96y7ydOze8YmVkqJYpozr6gWPBDyhbulT16FFVNadq3161YUNTG8Gdhx82tRJWrnTbuWlT8Qkw69ZNtXt3XbjQ/MYvvVSww5980hy3aVNQpAs/x4/rgbK19dG6E0J2yvF3LNNTlNHU3v1CHzhZRLABZSWIBQvMo2Fdjjp0DwAAFvRJREFUYxwbP96kIM5v3TjYREdDgwawbnvF4BoqsrKgXz/jJ4rJLrF4MTz+eN4MEw89ZMojPPqos+O338wa1syZFAt++AGmTuU//4Fq1cxyT0EYNsz8FO++i/El3VnCEqZt3kz1tN3UqFM2/74BouudbXgk5j+U+3Ymmc+/GLLzBhurCIoTqsZjyHEbXb7cuEzeemv4jMTuxMVBi98mmzWaYDF3rkl1OmgQWVnGU6hp02y9kIMaNUyK5pkzjTGdtm1NUZf33guefIGkTBk2pZ3NRx+Z37hyAeuyN2xoHI7eeQd06HVmvagE4XIdjWkWfNdRFwkJ0H7i3bzPIBj1D3M9lgCsIihObN5snuocRTBunHET9HQTDAfx8RC3Y3722n1QmDLFeEv168f775u0Gk8+CWW8OI38/e/GaWjUKNCoaONu+vXX4a2tnB+q5qY9YQIvvWRmW7lyC/rN8OHG03Zdzc7G2ry+5ARDHV6xGYCqbUKnCACuuVZYefdb/KWNWPbiDyE9d7CwiqA4sXy5ee3alWPHzD1x8GCTdj8SiIuDbSdqovv3B8dNMy3NVGPp35+MshUZPRoSE02YgDcqVTITlB9/NCstXHed8S2dMSPw8gWKX3+FuXM5fjKat982HrJ1C+kmf+WVZlnp9a0lz4302O+bOEF56ratHfJzPzG2Cg/2WEynb59k0aKQnz7gWEVQnOjf33grJCYybZpJsx+uSGJPxMXBPmKRzEzjuxho5swx3lJDhzJ5sikW9dRT+afTufVWY78YNQq0dRIkJRkn+0hl+nQ44wzG77mS48dNmqHCUr48XHstjPuuEZlNm5coN9KNFRKZwE00aRr6ddEyZWDCh1WoWxdGX7aEI0+/HHIZAoo3K3KkbqXea8ihfXvVVq0iy3Fh3jzV65hkXFXWrw/OSVas0LTj6dqwoWq7dv5//gkTjFgff6yqS5ao7tkTHPmKSkaGap06mtnvSq1TR/Xii4s+5OLFTprwC+9VLVtW9dixog8aATz4oPk4mZnhk2HFCtXxZW5TBT31yRfhE8QPsF5DJYDDh00Q2fz5LF1qPGUixUjsIi4O9lOTU2UrwZEjwTlJUhJvTSzDli0mctjfz3/99SYFx6OPQmZy29CVbyso8+bBzp38XHcIO3fmn07CH9q2NYn4nt19s0nU55akrzizeW0ajRqFN8FqUhLUeGcsy0nm1JDrjUGmOOJNQ0TqVmpnBN98Yx7rvv9eR4xQLV9e9eDBcAuVk8xM1TNisvThh4Mw+LRpqtdfryf2HNU6dVTPO6/gs6EPPjBf4bvvqurXX6tee21kTalUVZcv16xhw7RDq2PaunXgxPu//zOf/fffAzNe2Dl4UBX0f61fC7ckqqo6ZsR6PUQV3R3XIW9AS4SAnRGUAObPh6gojiZ0ZOpUYyCtVi3cQuUkKgoaxkl2ZuCA8tZbsGABr0+qyM6dBZsNuBgwwDwdjx4NGSm7TL6ihQuDIGwRSE7m+2sn8uvqitx/f+BmfEOHmonAF/9Za1JyqwZm4DChG81FVr5h6A3Fnrjvf014uc1Eztz8K5sefiPc4hQYqwiKC/PnQ3IyUz+vzLFjkWUkdqdpgzRumnWtSYoXKHbsgNmzSbvqWp59TrjoosK5xEdFwTPPGC/cCYcGGEtqJMUUrFsHq1fz4osm3MHJpxcQatUycXj7PpxjrM9//hm4wcPAwaWhSz/tD9HRcPfs/tx81pecP/1OduwIt0QFwyqC4kBGBixaBF27Mn68We/t1CncQnmmXqMzuODAB04y/AAxfTqoMiFtKPv2uWUVLQS9ehkl8viLVci49ArjRhoppQeffZbMTl346fs07r7bVJsLJMOHw/vHLjFvirn30MFlRhHEnhsZigCMG/e93/fl0LEy3HL5btL+2BBukfzGKoLiwN69kJzMhvo9WLo08ozE7sTFC/uIJX1XAGsSTJlCRpt2/GNicy6/3FTgKiwiZlawezd8UvE6k+D/m28CJ2thSUuDjz9mQe0BxFQsy623Bv4UvXtDRp0GbK6cWOwVQdqfmzhMFeLaRNb6aGIiTHpX+ffSPuzreiV67Hi4RfILqwj8JS0tfOuqderA/Pk8v/4qKlQwMVGRSny88RxKTQlQ5G5GBvTowcx6d3L4sCk/XFS6dIHLLoM7Pu1Fes9e3sOSQ8k338CRIzy7aQg33xycIMEyZUz+oQ+O9UV//tkEohRTllS/iDFRj9CgYeQ9EQ24Slh+7QvUOfQH6y68rXjYY7xZkSN1C7nX0NGjqiNHqkZHqyYkqL74ouru3aGVITNTDx9WrVhR9aabQnvqgvLLL6o/0l33tTo/YGPu2WM+++DBARtSly83XjSjRgVuzCIxYIAeLR+rZ8gp3bgxeKdZu1a1G3P0VJlyJvCjuHHkiOpvv+mgQapNm4ZbGO9kZKhObmbSv657YFy4xVFV6zVUNGJi4NtvTXhm1arGsbtuXZMEP1Q0a8a6a0Zz/HjkGoldxMXBXzTmKAXMkOYJVZg7l+efySQ11eQUChTJycYY+9JLsHvdYWOoDRcZGWQtXsJ7mddyxcCYoNZhb9YM6NKV9vEH0C5d8+0fMRw/Ds8/j8bHc+LiK5j3Y3pIitEUluhouOyXR5lboTcNXrybXV8tDbdIvvGmISJ1C8mMYP581X79Tkdgnjjx/+3deZhU1ZnH8e/b0A000LILsoiIig6jsigwalzDLo6jxBU1PASZjAuKUTTGLY8QyBijKIkKRsyMxieJC4qKu6PYKK4ooLLJ2q0IsjQ0Nt31mz/ORUuEbrq7tqbez/PUU1X33jp13q7b961z6t5zvl+3cGEY6H7z5vD8oYfC808/TU5dVqyQQBPaT9FRR2Xeae+7isWkhg2lceMSUFhhoQS6JPevuvjiBJS3i88/Dw29Va2Okk48MfFvUJmionBp7NatkqSHr/tEDSjV3LnJf+vp00Nr6K23kv9etbZ1q3THHapo1VoCvZQ3SL15R927h90j0302Z51m5I7UqUd/nfZpMKikRZD2A3t1b0lNBGvXSiNGhD9L+/ah/6Aq48aFowlIxx0X/st2JolEeOQRCdSD9zR1auKKTaZu3aSzzkpAQZddpm/rNVTzepu0dGkCytuN0aOlm3NuC5/fihXJeZN4JSVhxpjGjaX69VU0Y7ZuvFFq1ixcJJcKmzdL/Rq+r+Wteu/dPp5Gq6aFCylfzPmp+vKWBgyQXngh878QxXvyybB7jRxRplh5+sbD8ERQlfJyadIkqUkTKS8vdBxHM2BJ4fjw8MNhB9ztzwNFRdLkydJhh4U/6cCB1Xv/HTukNWvCGDizZn33LVFPPy11767S+o1VkL9DmzbVPMRUurnX03q3yYnSxo01L6SsTOUtW+vvNlxjxiSsaj+yapXULW9p+NwmTkzeG1VUSPffL7VtK4FW9z1LI0/4XGZhFrUhQ0JjM1XGnrNWAn17y4TUvene2L5dsXunasmY32vYMMmIqW/9eRo5sm5fFT1h/Ca9RV+9OWxS2urgiaAqsVjoGjj9dGnJEpWVSa+9Jl17bRjYLXRWf39r104aNEgaP17629+kRYtCLlEsFrqV5swJ5RYXhwImTgxDJNx5Z+hG+uKLsP6xx6Q2bcKRIP4Ndu7x06ap/NBumlD/Nxo1KvFhJ8tDp8xQrQeee/ZZCXRW/Se1alXi6rY748ZJb3CcSg8+InlfNWMxlfY+Xl+076ehLeZ81+i8+ebUNER29cYb0rv0VPGhKWqGVKWsTDv+/IBKWnaSQM/TX61axnTTTeHfqK6rKI/pjQOGawf1NKfjz/T6iPu18vUknhWwG54IdmfxYmn48PBNXFLx0hI9+KB09tlSQUH4y+TmSqeeKt1xR2hBv/xymDL4ooukI48Mc/TuPHY3aiQde2zoarj33pALSuZ+LJ1wwg8P8rm50iuvhDoUFkqjR6v81zdp86SpKpr6uD59qFCvP7dVTz0VWiGjRoWXvfNOYsJOhb//fJYEKnmp5p2435xzqTbQTOMuS/64LevWSVc2+FP4Qy9YkLiC581TxaDBeuaBterfX2rBeuVYTEOHSjNnhoZgusRi0j0tblQ59aQNG9JXEUlbZr6ib1ocJIEK6aORHWbrvj/HfvDT3L5g06pNer3bL7Q254Dvjgd9Oq3V6NHSrClLtf7zr5P6/mlLBMBA4DNgCTB+N+sbAI9F698GOldVZq0TwZYt0vjxiuXlqTy/iR4+d5Z69fr+OH3AAeHg+8QTVXf1b98uffBB+L147Fjp5JOl5s1/eNzv2lUaM3C5rh/2iS4c9LVOPimmnj3D8jZtwjC6u7Y4dr0dd1zd6hN9ccLbEmjZlJoPy3v+z3aoZ8MFKTtTd+L4jTqcBZo3LwGFLV+uLaefJ4HWWWudxCvq0EG65RZp5coElJ8gMy6dI4GK7nosNW8Yi4XuwkWLVPrcq1ry+Ee64grp2Ebz9TbHaPyRs/TM07G0DiudCrGKmJY+s1AvXPAXDRsmNW0q/ZMzVYFpQX4vvdLnOr036UVtW5/YTFhZIrCwPvHMrB7wOfBTYDUwDzhP0sK4bX4JHClpjJmdC5wpqZL5pqB37956t4bDF5Q88Cg5111D/jdreazBRYz99nd8ldOOfv1g8GAYMiQMK1ubq3alMKXuhx/CRx+F+/nzwygGBQXhDNT4+6qWFRRAq1bhdLS6Yv6TyzjyzIP54MqH6PHHi6v/+vnh9M7rrw9XAafC5s3QpQuc0/EtrurxGtY4H8tvRE6TfOo1aURZ/6HkNm1Iow1ryNu6gdyCfPKahXXWOB9ycykrg+XDr+WgZ+6mQjncydXMH3QdI/6rgIEDM+8zXLOygjkHnk/JBZcy8n9OqXE5EpQsKWbjoiK2Lili+/IiKtYUU5Tbidc6jmDdOrj18e603baMhrHS7143nZH8Z+50zjsPrhorju6ReReHpUJ5OSycMY8NjzxPs/de4ohNheSxg7nWlxtPKeS002DIoYs5YmgX6uXVfCcys/ck9d7tuiQmgn7ALZIGRM+vB5A0MW6b2dE2hWZWHygGWquSStU0ETz6KJSdfwn/wif8Zr8ptDq9H0OGhLFnWrSodnGuEuuXb2ZNl+OZ0vTXvNrmHGat+lfyY1sRhqIsO6vJudzVOgwa9PzywwCQGcJoW7aSCQ1uY9zacSn9bKZMgeIrbud2bvzRupZ8zQZacjs3cAMTf7S+TdNStlY05N5tl5DfOIe1Y25j+FUdajzFZKoMGRKm8ZyRN4pepW9iEiByiLEs9zB+0W4WEkwt/g+6f/sehE8RQ3yc24tLmj/F+vXwSdkhHMIPx9Z5nDO5sNHjtG4NE7ePpUF+fb5t2Y5Ym3bUa9+WnEMO5vgLDsz4v1GqlRSXsOj+N/jogwqmLBvKZ/O3s4EWLK/XldmT53P11TUrt7JEkMxr69sDq+KerwZ2HSXmu20klZvZJqAl8IPxCcxsNDAaoFOnTjWqTJ8+8Mh199DtjHyeOTYn476d7UtadC7g7pvmU7oU+gjW0I+8itK4S+1F/Tad6dElPFtb2hMj9IQZ4it1Z/A5nVOeoC+/HD456QbmrLuGii3bKN9SSqxkGxUlpdzdphnbd0DjFRfyzMqesL2UnNJtUFpKzvZtnN+jATHB/gMepP/AurN/3XorNG0K5QsOonjTFsCQ5YAZW/I70qdHaCFv+bgHy0oKwHYmc2NzQVcG94CWLWFB8WRW7wcNOrejcdd2NDtsfwZ0aMi2xjvf6Y/pC7KOadK2CcfcNIhjgFHAlyvgw99P47OPSpOWNJPZIhgODJA0Kno+AjhW0uVx2yyItlkdPV8abbPHEctq0zXknHPZqrIWQTKHmFgNdIx73gHYdZTu77aJuob2AzYksU7OOed2kcxEMA84xMwOMrM84Fxg5i7bzAR2/pp4NvBKZb8POOecS7yk/UYQ9flfBswG6gEPSlpgZrcRTmOaCUwH/mpmSwgtgXOTVR/nnHO7l9SB2CU9Czy7y7Kb4h5vB4Ynsw7OOecq58NQO+dclvNE4JxzWc4TgXPOZTlPBM45l+WSdkFZspjZOmBFDV/eil2uWs4CHnN28JhdVQ6U1Hp3K+pcIqgNM3t3T1fW7as85uzgMbva8K4h55zLcp4InHMuy2VbIrg/3RVIA485O3jMrsay6jcC55xzP5ZtLQLnnHO78ETgnHNZrs4nAjPraGavmtkiM1tgZldGy1uY2Ytmtji6bx4tNzO728yWmNl8M+sZV9bF0faLzaz6k+2mQCLjjdYXmNkaM7snHfHsjQR/xpOjMhZF22TkRLk1iLmbmRWa2bdmdk1V5WSiRMUcrWtmZv8ws0+j8vqlI6Y6Y0+z2teVG9AO6Bk9bgp8DhwBTAbGR8vHA5Oix4OB5wAD+gJvR8tbAMui++bR4+bpji9Z8caVdxfwCHBPumNLwWf8b8AcwrDo9YBC4KR0x5egmNsAxwC3A9dUVU6640tmzNG6GcCo6HEe0Czd8WXyrc63CCQVSXo/erwFWESYC/kMws5AdP/v0eMzgIcVzAWamVk7YADwoqQNkr4BXgQGpjCUvZLAeDGzXsD+wAspDKHaEhizgIaEA0MDIBf4MmWBVEN1Y5b0laR5wI69LCfjJCpmMysAfkKY7wRJZZI2piSIOqrOJ4J4ZtYZ6AG8DewvqQjCDkb49gBhx1oV97LV0bI9Lc9YtYnXzHKAO4Bfpaq+iVCbmCUVAq8CRdFttqRFqal5ze1lzNUtJ6PVMuYuwDrgL2b2gZlNM7PGSaxunbfPJAIzawL8ExgraXNlm+5mmSpZnpESEO8vgWclrdrN+oxU25jNrCtwOGH+7PbAKWb2k8TXNHGqEXNKykmFBNS1PtAT+JOkHsBWQpeS24N9IhGYWS5hx/lfSY9Hi7+M6wJpB3wVLV8NdIx7eQdgbSXLM06C4u0HXGZmXwD/DVxkZr9LQfVrJEExnwnMlVQiqYTwO0LfVNS/JqoZc3XLyUgJink1sFrSzpbPPwiJwe1BnU8E0Vkf04FFkv4Qt2omsPPMn4uBp+KWXxSdWdIX2BQ1N2cD/c2seXRWQv9oWUZJVLySLpDUSVJn4BpCn3pGfmtK4Ge8EjjRzOpHB5wTCf3QGacGMVe3nIyTqJglFQOrzOywaNGpwMIEV3ffku5fq2t7A44ndHXMBz6MboOBlsDLwOLovkW0vQH3AkuBj4HecWWNBJZEt5+nO7ZkxxtX5iVk9llDCYmZcKbQfYSD/0LgD+mOLYExtyV8E94MbIweF+ypnHTHl8yYo3VHA+9GZT1JBp4BmEk3H2LCOeeyXJ3vGnLOOVc7ngiccy7LeSJwzrks54nAOeeynCcC55zLcp4InHMuy3kicC4NzKxeuuvg3E6eCJyrgpn9Nn4cfzO73cyuMLNfmdm8aM6DW+PWP2lm70Vj6o+OW15iZreZ2duEIT6cywieCJyr2nSiIQ6iUVvPJQxffQhwLOEq1l5xA9iNlNQL6A1cYWYto+WNgU8k9ZH0ZioDcK4y9dNdAecynaQvzGy9mfUgzN/wAWFClP7RY4AmhMTwf4SD/5nR8o7R8vVABWFANecyiicC5/bONMKYTG2BBwkDmU2UdF/8RmZ2EnAa0E/SNjN7jTAZDsB2SRWpqrBze8u7hpzbO08QZqw7hjAq7WxgZDR2PmbW3szaAPsB30RJoBsZPMy1czt5i8C5vSCpzMxeBTZG3+pfMLPDgcIwejIlwIXA88AYM5sPfAbMTVedndtbPvqoc3sh+pH4fWC4pMXpro9zieRdQ85VwcyOIMxR8bInAbcv8haBc85lOW8ROOdclvNE4JxzWc4TgXPOZTlPBM45l+U8ETjnXJb7f4Ql2mfaqqRSAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "\"\"\"Yearly TC damages in the USA, normalized and current\"\"\"\n", "\n", "yearly_damage_normalized_to_2019 = emdat_impact_yearlysum(\n", " emdat_file_path,\n", " countries=\"USA\",\n", " hazard=\"Tropical cyclone\",\n", " year_range=None,\n", " reference_year=2019,\n", ")\n", "\n", "yearly_damage_current = emdat_impact_yearlysum(\n", " emdat_file_path,\n", " countries=[\"USA\"],\n", " hazard=\"TC\",\n", ")\n", "\n", "import matplotlib.pyplot as plt\n", "\n", "fig, axis = plt.subplots(1, 1)\n", "axis.plot(\n", " yearly_damage_current.year,\n", " yearly_damage_current.impact,\n", " \"b\",\n", " label=\"USD current value\",\n", ")\n", "axis.plot(\n", " yearly_damage_normalized_to_2019.year,\n", " yearly_damage_normalized_to_2019.impact_scaled,\n", " \"r--\",\n", " label=\"USD normalized to 2019\",\n", ")\n", "plt.legend()\n", "axis.set_title(\"TC damage reported in EM-DAT in the USA\")\n", "axis.set_xticks([2000, 2004, 2008, 2012, 2016])\n", "axis.set_xlabel(\"year\")\n", "axis.set_ylabel(\"Total Damage [USD]\");" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.6" } }, "nbformat": 4, "nbformat_minor": 4 }