{ "cells": [ { "cell_type": "markdown", "id": "c4d5ef32-f24e-4ce3-b4df-9edb52de1bfb", "metadata": {}, "source": [ "# ECMWF API plugin for EODAG\n", "\n", "This tutorial will show you how use eodag to download data from [ECWMF](https://www.ecmwf.int/) using [EcmwfApi](https://eodag.readthedocs.io/en/latest/plugins_reference/generated/eodag.plugins.apis.ecmwf.EcmwfApi.html) `eodag` plugin. The API plugin and this tutorial have been developed in the context of DOMINO-X." ] }, { "cell_type": "code", "execution_count": 1, "id": "1298613c-5383-48c4-bc0a-f4136130daaf", "metadata": {}, "outputs": [], "source": [ "from eodag import EODataAccessGateway, setup_logging\n", "\n", "setup_logging(1) # 0: nothing, 1: only progress bars, 2: INFO, 3: DEBUG\n", "dag = EODataAccessGateway()\n", "dag.set_preferred_provider(\"ecmwf\")" ] }, { "cell_type": "markdown", "id": "fd6acf76-19ac-4e27-be87-644d2d978ba0", "metadata": {}, "source": [ "### Search (build download request) from an existing product type:\n", "\n", "In this example, we will request data from [TIGGE public dataset](https://www.ecmwf.int/en/research/projects/tigge), using a product type already configured in eodag (`TIGGE_CF_SFC`) or using a custom request to access the same data.\n", "\n", "Retrieval request can be customized from ecmwf web ui on https://apps.ecmwf.int/datasets/data/tigge. \n", "\n", "For performance purpose, we add a `param` (total cloud cover, `tcc`) to the request as `TIGGE_CF_SFC` is configured to ask for all available `param` values. See [ecmwf availables parameters](https://apps.ecmwf.int/codes/grib/param-db/). You can add more [MARS keywords](https://confluence.ecmwf.int/display/UDOC/Keywords+in+MARS+and+Dissemination+requests), that will overwrite default values configured for `TIGGE_CF_SFC` product type." ] }, { "cell_type": "code", "execution_count": 2, "id": "2181ef4d-cc88-458c-a633-f009cedc43ef", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "1 product built EOProduct(id=TIGGE_CF_SFC_20211101_31613e05afb70ccea479a4e2ce057227b15de05b, provider=ecmwf), \n", " having param=59/134/136/146/147/151/165/166/167/168/172/176/177/179/189/235/228001/228002/228039/228139/228141/228144/228164/228228\n", "\n", "1 product built EOProduct(id=TIGGE_CF_SFC_20211101_28b0b1a0da723aba1682661621bd7b21fe12f1ea, provider=ecmwf), \n", " having param=tcc\n" ] } ], "source": [ "# Request for all parameters\n", "products_from_product_type, total_count = dag.search(\n", " geom=[-60, 30, -130, 70],\n", " start=\"2021-11-01\",\n", " end=\"2021-11-02\",\n", " productType=\"TIGGE_CF_SFC\",\n", ")\n", "print(\n", " \"%s product built %s, \\n having param=%s\\n\"\n", " % (\n", " total_count,\n", " products_from_product_type[0],\n", " products_from_product_type[0].properties[\"param\"],\n", " )\n", ")\n", "# Request for total cloud cover\n", "products_from_product_type, total_count = dag.search(\n", " geom=[-60, 30, -130, 70],\n", " start=\"2021-11-01\",\n", " end=\"2021-11-02\",\n", " productType=\"TIGGE_CF_SFC\",\n", " param=\"tcc\",\n", ")\n", "print(\n", " \"%s product built %s, \\n having param=%s\"\n", " % (\n", " total_count,\n", " products_from_product_type[0],\n", " products_from_product_type[0].properties[\"param\"],\n", " )\n", ")" ] }, { "cell_type": "markdown", "id": "59577515-7446-4e1f-ba33-3ef87fb7c80d", "metadata": {}, "source": [ "### Send product retrieval request and download when available\n", "\n", "- download performed using [ECMWF credentials](https://github.com/ecmwf/ecmwf-api-client#registered-access-recommended) set in `~/.config/eodag/eodag.yml` as for other EO providers:\n", "```yaml\n", " ecmwf:\n", " priority:\n", " api:\n", " outputs_prefix: /data/eodag_data\n", " credentials:\n", " username: john.doe@csgroup.eu\n", " password: my-ecmwf-api-access-key\n", "```\n", "- you can check your request status from https://apps.ecmwf.int/webmars/joblist/\n", "- See [ECMWF guidelines](https://confluence.ecmwf.int/display/UDOC/My+request+is+queued+for+a+long+time+-+Web+API+FAQ) if request is queued for a long time" ] }, { "cell_type": "code", "execution_count": 4, "id": "c5844159-d67f-405b-aaf8-2017ff0010e7", "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "d19e81147f054136982b51ffd1e75d8a", "version_major": 2, "version_minor": 0 }, "text/plain": [ "0.00B [00:00, ?B/s]" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "'/data/eodag_data/TIGGE_CF_SFC_20211101_28b0b1a0da723aba1682661621bd7b21fe12f1ea.grib'" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "product_path = dag.download(products_from_product_type[0])\n", "product_path" ] }, { "cell_type": "markdown", "id": "db727f7e-68b7-4570-a15b-28c7cd6d1c59", "metadata": {}, "source": [ "### Open dataset with xarray and [cfgrib](https://github.com/ecmwf/cfgrib), then plot using matplotlib" ] }, { "cell_type": "code", "execution_count": 5, "id": "a733a270-b984-45cb-80a4-a35a1b2112f4", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 6, "id": "1d1e5b4c-0919-493a-af0c-1b5604909978", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:           (latitude: 81, longitude: 141)\n",
       "Coordinates:\n",
       "    number            int64 0\n",
       "    time              datetime64[ns] 2021-11-01\n",
       "    step              timedelta64[ns] 00:00:00\n",
       "    entireAtmosphere  float64 0.0\n",
       "  * latitude          (latitude) float64 70.0 69.5 69.0 68.5 ... 31.0 30.5 30.0\n",
       "  * longitude         (longitude) float64 230.0 230.5 231.0 ... 299.5 300.0\n",
       "    valid_time        datetime64[ns] 2021-11-01\n",
       "Data variables:\n",
       "    tcc               (latitude, longitude) float32 100.0 100.0 ... 8.406 10.62\n",
       "Attributes:\n",
       "    GRIB_edition:            2\n",
       "    GRIB_centre:             ecmf\n",
       "    GRIB_centreDescription:  European Centre for Medium-Range Weather Forecasts\n",
       "    GRIB_subCentre:          0\n",
       "    Conventions:             CF-1.7\n",
       "    institution:             European Centre for Medium-Range Weather Forecasts\n",
       "    history:                 2022-05-11T09:27 GRIB to CDM+CF via cfgrib-0.9.1...
" ], "text/plain": [ "\n", "Dimensions: (latitude: 81, longitude: 141)\n", "Coordinates:\n", " number int64 0\n", " time datetime64[ns] 2021-11-01\n", " step timedelta64[ns] 00:00:00\n", " entireAtmosphere float64 0.0\n", " * latitude (latitude) float64 70.0 69.5 69.0 68.5 ... 31.0 30.5 30.0\n", " * longitude (longitude) float64 230.0 230.5 231.0 ... 299.5 300.0\n", " valid_time datetime64[ns] 2021-11-01\n", "Data variables:\n", " tcc (latitude, longitude) float32 100.0 100.0 ... 8.406 10.62\n", "Attributes:\n", " GRIB_edition: 2\n", " GRIB_centre: ecmf\n", " GRIB_centreDescription: European Centre for Medium-Range Weather Forecasts\n", " GRIB_subCentre: 0\n", " Conventions: CF-1.7\n", " institution: European Centre for Medium-Range Weather Forecasts\n", " history: 2022-05-11T09:27 GRIB to CDM+CF via cfgrib-0.9.1..." ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds = xr.load_dataset(product_path, engine=\"cfgrib\")\n", "ds" ] }, { "cell_type": "code", "execution_count": 7, "id": "92531254-85c0-49d0-a676-f2972efb42b5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEXCAYAAACkpJNEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACatElEQVR4nO2dd3gc1dXGf2dX3ZIs994LBgM2vfdOaAm9BQgJIfQQCJCQBAIkJF/oEHrvvYbee7HBFXdb7t2WbdlW273fH+8debVWWVVLeN7n2Wd3p9y5c2d25573nPMec84RIkSIECE2P0Q2dQdChAgRIsSmQfgACBEiRIjNFOEDIESIECE2U4QPgBAhQoTYTBE+AEKECBFiM0X4AAgRIkSIzRQ/iQeAmRWa2YGbuh+NgT+H9Wb2eIrbH2hmxWYWb+i5m9leZjalIfuGCNEaYGZnmtnnm7ofbRU/iQdAa4WZjTSz0Wa2zr+PrGOXI51zpyfsf52ZjTezCjO7JnFD59z7zrlcYE49+uPMbHBCG58557ZIdf+Whpl1NbOnzWyBma0ysy/MbJekbU4xs9lmttbMXjGzjn55ppk96NetMbMxZnZYwn4ZZvaCf/A6M9u3jr7Uur2Z7WdmH/l+FqZwbmZm/zKz5f71LzOzhPX3mdkU/4A/s+7Rqhtmdo2ZPdEUbTUlzKy/H7t1Zja5rU/m2hLCB0ACzCytCdvKAF4FngA6AI8Cr/rlqWI68Efgf03VrzaGXOA7YAegIxrD/5lZLoCZDQfuBU4HugHrgP/6fdOAucA+QHvgauA5M+uf0P7nwGnAohT7U9v2a4GHgMtTbOsc4BhgBLAtcCTw24T1Y4HzgO9TbK8t42ngB6AT8GfgBTPrsmm7tJnAOdcsL6AQuAwYB6wCngWy/Lozgc+TtnfAYP/5EfRDfgsoBr4AugO3AiuBycB2Sce6CvjRr384OJZffwQwBigCvgS2Tdr3Ct/PUiCtic7/YGA+YAnL5gCH1jJeB9aw7gngmvrul7Tdp36M1/oxPRHYF5iX1NblfizWAg+iP9a3gDXA+0CHhO139eNZhP6w9m2u+ynhmKuBHfznfwBPJawbBJQBeTXsOw44tprl8+rT99q2Bw4EClNo40vgnITvZwNfV7Pd58CZ9RyjK/y9twaYAhwAHOrHptxf/7F+2/b+Oi/0+1wPRP26M/1v7070G54MHNDE13Oo/93lJSz7DDi3hu07Aa/5++Bb4DoS/kuA29CDfzUwGtjLL++OJgidErbdHlgKpDf3fdtaX81tAZyAbrwBaJZzZj33vRrojG6Qr9BsqDPwAnBz0vanAoegP4Ghfl/MbDs0M/stunnuBV4zs8yEfU8GfgYUOOcqkjtiZuPMrKiG13+Tt/cYDoxz/k7zGOeXtzicc3v7jyOcc7nOuWdr2PRY4CA0hkeiP/8/AV2QxXgRgJn1QpbJ9Wh2fhnwYk0zNzN7o5YxfCOVc/AUWgayjEBjOTbhHGegP7mh1ezbzS+fmMqxWgBV+u4/N/reMLMtgAuAnZxzeeg3Ueicexs9MJ/113+E3+URoAIYDGyHJi6/TmhyF2AG+t39DXgpoNmqOXZDrvFwYKZzbk3CstrG4i6gBOgB/Mq/EvEdMBLdk08Bz5tZlnNuEfAx+l8JcDrwjHOuvIZj/eTR3A+A251zC5xzK4DX0YVJFS8750Y750qAl4ES59xjzrkYsia2S9r+TufcXH+sG9CfOsjUvtc5941zLuacexQ9UHZN6udc59z66jrinNvWOVdQw+u8Gvqfi2ZNiVgF5KU6AJsIdzjnFjvn5qOZ2DfOuR8SrkMw7qcBbzrn3nTOxZ1z7wGjgMOra9Q5d0QtY3hEXZ0ys3zgceBa51wwrimNsZmlA08CjzrnJqcyCC2A5L6vAnIT/QANRAzIBLYys3TnXKF/MG4E/1A8HLjEObfWObcEuAU4KWGzJcCtzrlyP2mYgiZLG6GB1zjl34mZRdEE5a++vxMQLZjYhyecc8udcxXOuZv8WAR+rkfRfRu0dTK6pzZbNPcDIJErXYcudqpYnPB5fTXfk9uam/B5NtDTf+4H/CFxNgL0SVifvG9ToRjIT1qWj8zy1oxUx70fcHzSuO6JZmZNCjPLRhOIr51z/0xYVecYm1kE/cjL0Mw4leP1NUVYFZtZcaM6r/b+lNDePTX0PR8oTrIY6w3n3HTgEuAaYImZPWNmPWvYvB+QDixMuIb3Al0Ttpmf1KfE31ZToD6/ky5s8O0k9qcSZnaZmU3yzvgiRHF19qtfRQ/GAcjKXeWc+7bxp9B2samcwGuBnOCLmXVvgjb7JHzuCyzwn+cCNyTNRnKcc08nbF/rj87MJib+ISS97qlht4nAtkkzum1pPRREYzEXeDxpXNs5526sbmMze6uWMXyrpoN4qu4VxLv/Nmn1ROREDbYdiGZ8U/13Y4Mf49hUTX3n3BxPk+Q6RVo1Cs65fyS0d251ffefm+TecM495ZzbE/3BO+BfwaqkTecia7hzwjXMd84l0i+9ku7hxN9WFTTwGk8EBppZ4oy/prFYiuiq5N96cPy9UNDECchXVYCsCQPwVuxzyAo4nc189g+b7gEwFhhuCpPMQrOVxuJ8M+vt+ck/I5oI4H7gXDPbxYR2ZvazpBuuVjjnhif+ISS9zq1ht4+ROX6RKSQxmH1+mOpxzSzdj08ESDOzLG+61rT9mVZ7COJiYGCqx68DTwBHmtkhZhb1fdvXzHpXt7Fz7rBaxvCw6vbx1M0LyPI4wzkXT9rkSd+HvcysHfB34KUEPvluYEsUXrsRveevS5b/muHPoUYKprbtzSzi16Xrq2VZ7RFfjwGXmlkvP0P/A+Ljg2Nl+PYMSPftRfy6fc2s2kmLmW1hZvv7B2cJGrtg3BYD/YN2nHMLgXeBm8ws35/DIDPbJ6HJrugeTjez49F4vlndsRtyjZ1zU1GAxt/8Of4cTZRerGbbGPAScI2Z5ZjZVsAZCZvkoQfEUvR7+SsbWxePIV/kUYQPgE3zAPAX/e8oqmQainRoLJ5CN/NM5LS63h9rFPAbFMmwEjkQz2yC49UK51wZCvP7JYqS+RVwjF+eKu5HP+CT0UNtPZq51IQ+KGqjJlwDPOrN/RNq2a5OOOfmAkcjB/FSNJu8nKa9p3ZHEVwHA0UJs8m9fB8mAueiB8ES9AdwHoCZ9UMWw0hgUcK+pya0PwWNaS/gHf+5Xy39qW37vf33N9GsdD26H2vCvYjWGg9MQA71exPWv+vb2B24z38OHPl9UBRRdcgEbgSWIQq2K4qQA3jevy83syC89JfIsR5E0L1AVRrvG2CIb+8G4Djn3PJazqshOAnY0R//Rn+MpTVsewGiIRehB+bDCeveAd5GFuBs9ACsQu86575AD8TvnXOzARLvKT+ZqKT9TPRdjRZqW4c1knIM0UQwZeT2QM7vM1LY/gA0S8oEDnfOfWRm7wIXO+cmNW9vQ2xKmNkDwPPOuXea+ThnAr/2dNJPBmb2IQoffmBT92VTI3wAhAgRolr8FB8AZrYT8B7Qx1UNPd0sEWYChwgRYrOAmT2KaOdLWvLP38weMrMlZjYhYVlHM3vPzKb59w5+uZnZ7WY23ZR/tH2z9i20AEKECBGi+WBme6Nw18ecc1v7Zf8GVjjnbjSzK1HU0hVmdjhwIcrP2AW4zTm3S01tNxbNagH4iIQxCa/VZnZJTU+/ECFChPipwTn3KbAiafHRbEhiexQFjATLH3PC10CBmTV5bk2AJhM/qw7OuSn47F8fvjgfZZNeCXyQ8PS7EumX1IjOHaOuf5/05uxuFTgfMl3mLaRy/6xcWq7o0eISKUlYROtdmdZHMmMA9MouIt30OctH7JnCkSvfl8TURqeoIhSj/hippIK6pE+W0l6pIe7bDMYgiCEs9wctirUDIDtSRjtTUFMQPVnht1kVzwZg+XJF4ZkX2EjvpO17ZhQB0K7Ria8b+rkmrvGbv6Az0RVrddwcRW1WZCl6Nq1YqQBlHXQv9emyDICo6SyzkqIro0Qo9yOwKqaozjKnn023tHUALIvpGEvW6lyzFvsRK/NpB3F9T7S2K6NNM32kaLIl7vdNttAtU/dMWUEaGat1f7n1JclDUi0sovGp6KhrM6iH8jQzao4sbhYEZ9R0d+zGGD2udJlzrsGCcofs184tXxFL9VgTUcRRgPucc/fVsVs3H4YLimjq5j/3omrk0jy/bCHNgGZ9ACThAGCGc262mR2NhMhAT7+PqeMB0L9POt++07e2TZoU5V4SaE5M13VRhfLW7lu8LwCfTZWqclqmtnNztT5r8GoA/rHNK3SPKsN9i3Rtk67wazJNfz63rhwEwFntlfPSPqI/kkgKhlnc/ynFfGh8etMJmVLq86XKnX4A6/xYLIjpGC+v2gGArbPnsUvWPAAy/B/acv8n/MaabQF4/MmDAMjygYM9TpkFwF/7vQ7ATpmNf6gH1+rjEo3fX//6a/Kf+gaAyLBhABQNbw9Ax8+VwzTvF70A+M/59wNQENFDOLhWEf/3lBvJYklMUYFvrlUKxZyyTgBc0lGRlA+v2hKA2749AIBhN/sowvlKonZr1Xa8dMN/RCTTX+v+Pm0i5h8a/g8/Pk+/d1fhHwQxXYu0frpn5h3RjV4frNS2E6dW2aYmRLJ1jy77+UgAXvzrvwHom9ay6iTBvZvKfd5QRHtMm133VjVj2YoY37xTbUrLRkjvMaPEObdjQ4/lnHM15XU0N1ryAXASkn2Fmp9+VWBm5yAtH/r2asmuhggRYvOGq5xcNRMWm1kP59xCT/Es8cvnUzXTubdf1ixoESewz4hcAAx3zi02syKfph2sX+mcq9UPsP2ITPfF2z0rZ8/NhaacnRTHNeNbFi/3bQoFET3MciOZTXasVFHqykn3Jn9xvBSAb0o1A8zyPE2Wqb/jSzUD+lW+ZrLHztBsvqhUs9eTeo3iuDzpjC3wk89r5h4FwM86jwNg/xytb4lZ5iI/U3+4aHsefkOz8SGP+N9Vseia5fsrb8v5Ie80WtTsghu04A/D3gfgkJxCALpGN1aCCO6R4fedD8CAF9XG5HNlZUw65i4A3lpXAMBdZx8PQNp3k4mvr5qQHC3QNtZe41M6WDI8na7VBHZesdpc9r3mSHmFfr8yR6dROq6boZpAyW0HiGSL8nHbyGo9/XEl8p6al0xL/3QQ7TFtdGNm5cH/TSrI6VlY57FMdSjeSHAC/x+wPIEG7+ic+6OZ/QwluwVO4Nudczs39DzqQktNqw9DmXeBsFhNT78QIUKEaBWI0zQWgJk9jSjvzmY2D8lq34gKFJ2NspaDzPw30Z//dCSgeVaTdKIGtNQD4GQ20D+ggg5noEE4A6n0hQgRIkSrgMMRayJ2xDl3cg2rDqhmWwec3yQHTgHNTgF5ka45wEDnddzNrBNS5euLf/o56fjXiB1HZLmWdAIHWB5XNMmlc6VldWuftwHI8jRKiXeUron7d6flw9NzaCpMKy/mzeKtAdglR7VQds2s/dkdOC4vniNKZs6aAgBWFOdQWqzIk9v2fAaAmJPD8/mlOwFQVCbK4JhuYwD4TfuqAQjBzKjUlTPBRz9dPv04HWeyhF1zZ2t5wcHa99OtX0ntZJsAceKVDuyA7kqm2Q748UgA0g6sWlJ51em7AXD2nzQnST73RATO5wvn7wXAu6O2AeC3e0vv7+T8HwDonaaoqbNm78v3L0tos+e/q0r5BDTN8pNGArBCPnSsm2jEl3a/u8r2R398AYMf0PHtC9FtpMhZz7hF5zj9xJqEbJsHyTPq4BoFaEp6t7EU0MgRGe6Dt7rWvSHQudf8Rh1rU6LZLQDn3FpUiStx2XKqefqFCBEiRGtBvHaV+J8EwtCaECFChEiCgyajgFoz2owUREtRQIGZGpinK32kzAfrdOxRxQMA2D5XURpHtBOFEFBCq3zET5lzNUb7lCeVHa4rhr/UlW9kHq+MK6plhw9UZqDj5zrGmgNFWb2xq0oV907TfsH+k8pL6BOtmoIzpVzHD2Ly1zud869nH6Jzi+qcftFpNAC7+qD+RTHjRZ8TMGmNkhWfHPB+lXMN8N56HWMbnwDWvZromppQ3/GqDYUVkoA5YpRqy/Q6VsKptoOomegiH1mzRhTawtOHc9mFKi2RHDUT9Cvoz+RyjX3PqM49PyJaJ7inZpSv45r5qoxYdKwvSe0TvmKrlD8SxP3P+/PuAPzi+M/UdrGotbGfDwGg+zcx2r0zXu0HiWB1UEBpXVQYa8Z/Fd0yajflQORGsmrcpyE4b76UC97/QDI22+wh2nLJOl3z14Y/CUCHSNPRpMloLAU0YkSGe+vNznVvCPTqvTCkgEKECBHipwKHI7YZUEBtxgLYYUSm+/LtXlSgmXm5n+0Es6zmRpAdu9BnBi/16f9dovpe4gPLZ5TL3fHE4t35bY+PABiYptldXkRWQo6fjacnpeDXlA+QGLsfbHPBAoUGf/i6ZuAFUzUeBWM0Sy08QbOXbQ5Sluj9/d8AINvSN5pBBzPZwgqdS8xnwQ5Lb1dlfVNmG29q/GuFYuKfmKZx7PlzZWPP+pccpB23kURE8add6PNWEQDr+2gGG71EEgpvb/kyAAtjir9PJdfhD4s0K/72Bjncc15KKknr72uL6novukCz6cPPUp2fL5fKAo3d3Y28D6foc1FRncdNbDPaTc7NaReprbuOv58DsiXR0Rw5KcH9E3DqQZZ1cD81R2ZwYy2AbbZNd6+maAEM6rMotABChAgR4qcDq5wI/ZQRPgBChAgRIgkOiLcNcqRRaDMPgDiO9a6c78vkONo3q2WvTuBE7ZsmM3qi78f9S1U/O+LN21nFHQGYs6ID/7z8lwCU52qYF+6h+Pv4MDkLb97hOQD2ypKoV3YlNVT1slQXH31d9491vP1EOc3fXpIBR/9Zsef/9THoi2KiyALBsjcWbc3Ng1QadpuM7CrHG5Je1TEb0F61xWcnm+/NKVLXlLiioxyT5++kGh23jhsJQNkVXviuQkKSgw8tZFJfSbO07y1xv+v7vwdAsS/vXB+Zi2WlGuP0Yh8DX4PjNhB2637XdwB8NVVU1fwDNJ4d84GIn6F6kcG6nMBBmxULlNsw4Eq933r/Eax9W+d0TLu1KZ9Lqki+B4IclcvnHQ7AtCKN9TPDHwGgRzSbqD+n+tJCTZW9C4QWQIgQIUJsjnCED4AQIUKE2GwRdz/9B0CbiQIKvPJ9fVp9SypowgbTMlDQnFYhKuiTtVsA8NQsRXWsWKyiIFvcvhY3QVEazpOJ0a0UebJkD0UXrB6k5fvup1T++3pXlQdIpT9f+RDw/60eCcA5HRUt0r8GWmLIx2fS9TXFoC/eRTf4W7/4DwCDvHxF8tgGOQftk+LFE7cL+jOxTOOzJCaqI1AW3S1r432aE8n5HDVJQtSExJoII5+/BIAuo7Qu/QxFAV0xULIgwzOWAtAjqpP8/QLVUH/765EAZHTX+A26eGkl/VJfpG2he6dopO6drBUVZI6W0mpsxcoGtVmlfR8ZNP8kHefcc18D4Ox85bukSsnEibMspvOdWaHxGJGha5BtmdXu88pa/aavn/Izf6w4AwsUzXZcV1FgLyzR72vUl0MBGHfKbQAsiOk6HfKp8mG6va5jfP3MZY2KzNly20z32BupFeLaud/sMAooRIgQIX5K2BwsgDbzAIia0T4S3cjZ2NKWQGaQ9VkqffYPlqriVMV7csZm+4l3pHg9Lk9fgjjt+BRVw4ruqG27jJAK9tXd3wGg3Mkpm8psK1i3faZm3Heu1czw6ajyAs7vMBbYOE9i2r6PMLDoHB3/G93gh3a+EIAjtpRDdOxKZYpeOfAtALr7XIeRGXX3J3AsQ6xyDSTOqOU4zYlkVLt/bajPNd9wn+i4S2Iap1QzkAPHd6al02GoMp/zn/eZqwfrOl59oZR6c34mlfPPtn0BgN93/QCAj3ooa7d8jo4587e59Pu7rnld1bsCBLH7QT2DtPUagyXbZ5DVT479rh+qylnFzMKU2kxu29LSwX/u8WkRAC/76mZPd9GMeuCfJgNweEdZq8fnrqq2zQiRSjfsqZ//BoBdB2m8HumnLPFkp3DgeD5mewVFlLsKVvts9F2f+QMAOfN1r3ZaJqv5Z++dB0BGkbYbMl79i5eUpnTudSH0AYQIESLEZgsj5lp2crkpED4AQoQIESIJDoi3MLuwKdBmHgBpRKqIRzU39VMT3fBtqaiBV5YqpX/mx/0B6DFJ1Eb2GDnNKpYs3ajNQOir49MSVVseEV1z43kHAnBdD+nId7J2G/VldKkog/eLJVp2VaepVbZ5rP+7QKJ5Xb1ExuiyMm47+AkAtjpCdMR3JYpzPyJHDsq9H/o9ABfueQoAb+9550b9gfpdg0RKpbGo6/jlrmIjGm1Gue6d8WUanxzT9do5U20F2z9XLHrukByVYf31rKPpeo7i1ivmVx3zrnfIaZ/2v/4AHHjvMQDMXabqpr2f8DRXTLkFRUPSie+mmgGRryU9EdwTyYh6+rB0FwUZzNpfbUWHSswuI72ClbOU+xEtlbMyr49yUMrzdY7F3fW+3sval3XwBefT9V45wc2Ic+j2onZmFWth2q91vOxy0Vffvap6FPGjzI+PaMvqpFhOmaz7pucrauPHM0SXLuydmmRGuqWxx5ei1wY/KaopPkYUT7SdrmPJ7qJep56vY5w+UuPyxCdywHPhc7UeIxWEFFCIECFCbIZwzih30bo3bOMIHwAhQoQIkQQ5gUMKqNEwswLgAWBrNK6/Ag4BfgMEPMmfnHNvNndf6oOa6IWPihV5MecBxUt3XSHzPuMDSS9UJEZ3JKXom9fmj+0iGqfPrxTHvT4mM3ZBhba/tUjm9oAsDc/JeXPYKVPRGDtlVqUhgtjqRT69fpo32XfNrKrguesYlWwsKU8j71HlKmxzlSKFbur5adBhAEZdfRcA583bA4Dj/3U5AE/9UfkCQ9NbRoE1GcnXpCapiiXx9fSKimZI91b8rQsOBmDdsVpw+Vfv+62r0kk/bydabNtHFH0y4K/f1B2xk6HjzxmnSKz8mTrGnIO9HMYaXweiXwknnPMVAK/ctD8AHZ8do16sX1+1TR+Vs8/Nopmu6DRuo3Mt3VHnHzlexwtorC9K9L48Jipx/2wpm37nZSjWxBWfXxBVZFH3aDG9/fEqo9w+UNuPrNA9cKqve3F63jJ/9Kr3QFBDYrsvfk3/k0VvzbuhFwB/20IyE0FpzOTrllhiFGB6eZyKMvVn+qnaZ+Qdoua2b69+3P+9qLHB9+vajCrsB8AWxaKKCmksQidwU+E24G3n3HFmlgHkoAfALc65/7TA8UOECBGiXgidwE0AM2sP7A2cCeCcKwPKzFqXcyWx4tSSuGZiT68aCcB9EzQD6vimnE+d3pqm927KVIwF2b5BA34WZtHoRjPHwOGXPlNx42N+GAjA5GM14/79QgnLffiqnMPtZ6nV6/aOMfbw24ENTrfiuGLzS328fVFcY7pFuiyBqFWttvTiNg+rrYUHs+oPauOXnZU1XFOG5j29NVstvkp1Dcb6YvFLIkFmcHrlLK4l8zFW+2sUVFpLnlFmWYQdr/kdAH1Ok5VVdqzGqXh3jflZ7/8agFlHqCrW6+s0Xrtkesf4GTcDcPJ/f0HF/AVVjp/WQ9W5Tv9EWv7XPKfqXV+e8H8AdPW5BsGsOBjfL0rinHvP+WpDRhjxEcpste+UgxFprxWLT5SleXDeB1XOLcCiWHFllbp1cbU/r8w7gT13/fP2o/3xte+eWT5tHL1vyI5uVzmmwf3fL01/DTd1/57aEJzjfmNPBaBsaTaHjdNv4+IOKmS/PK44//+tk+jbWN/v1V6ocKtsOdxnlWr9E2N3ptdL/louU39WvK0Z/ufzCgDYYo4qucXX6V4MfsGVeRNNgNhmkAjW3L/aAYjmedjMfjCzB8wqQ1wuMLNxZvaQmXVo5n6ECBEiRMpwGOUuLaVXW0ZzPwDSgO2Bu51z2wFrgSuBu4FBwEhgIXBTdTub2TlmNsrMRi1dnlrmZIgQIUI0FoETOJVXW0ZzP77mAfOcc9/47y8AVzrnFgcbmNn9wBvV7eycuw+4D1QUvqGd2Gv8z9WZ2XLSRdt551lUzqctusvsn7uqgLVTCoAN8ccD2nkz8EvRIZVk0dLAGVYVFkkwG2vQZ4/7fYddK/P5mOsP1QpvvvZf/6PaypWx1H5yR/60veihO3uKdoh6Gu2UKScAsEPHuQB0TJO5fVz+GLXlY65jfvT+0/OjShmGSIqXf5WnWMaWyPG9KKaxaRcpZVCazP2glkBTFnCvCSVe4C27suZAVbM/3zLZ9RzRH7lRjfEl34vuWuvHYVCa+jutXJTZPXvrGty5kVjbgkpaL22wSiguvkXntMaXBTU/Nwmon6Ao+lePKVekw1RPUX06gZ7rqxf8C27u2EqNbfdXZgLw22zJdGQcoHtm+WwZy5H2ZewzWHTkJd3k0B7efn6VNiOIRkmmou4oEgXz6/baf3ZFBd2iup8C53mA5HKgY8pEH53+w68AKHhC2y85QGfwj4Oe56TcquJ0j65SUMMbC5QDsW83HffGbrpGe449HoCV3yphoffoCoIIzJXDRM0t21n96NTTS6B09fd/ZznHCyKiin7xsSQiOONZGgOHhRRQY+GcWwTMNbMt/KIDgB/NLFFm7+fAhObsR4gQIULUF3EiKb3aMlqCwLoQeNJHAM0EzgJuN7ORaOJTCPy2BfoRIkSIECnBOcIw0KaAc24MkKyVfXpD22uIDMGHW6sEYvHwsirL/7JoX61/TVE3vT9aT9oP43WcQFWwhlT9umAZ6TWm+cfL1I/4cqlMpvX0BlFE5+R8ZIMrVR+m/rk3d3WVTEQcmcSBOf/OsP8BcMbsvQF4/nNFpPy3434A7DNcUUoP9v1Eh6hl3JKjaQLqoNzzEz3TZdqv8Fr/C8o7sLRCUSsRKwQ2UCtNgZquddckRc8NtRpET0TNuKqbome6RgKFVY1bLImW6xH1P4F832bVgB+iBQUU7ysD9or/PAbAz3J0nCsXjwCg866qDxBE0oxZ1lvLx2m7yCff+37WjCrKnEB8uai13q+IZprZXgqtBdtr+VfbP5EQGZRdZRymlqsfa53W75BRNcrrwgLF0pf6P7iCiOP9dWr/X1MPAWCLjqJFB7dTLsoTH0tiod+b4rv+fMeLAJT/Q/2+a8a+AExa35P17TQeQf8u7TCrynuAA348Wv18T9RPj8NEYb129nPk+toThRWSeOgcUVtB1FsmVWm/WRW6Sccd+F8A2tNYGPFQCiJEiBAhNj84oKyNR/ikgjZTEWzYtpnuodd70T2q2U1NFa+qQzCzvW/VIABu+eYgrYjo3Ps9o5lQzqhZVCzTDKuuAtv1gnciznpaTrDJez+iw9dhxczxs59Xi4dz58uqlvSfEx8FYJ8sWQ/veaWvh/aQ47F4VzkqF+6um/ffx2rWGmR9bpNRTKeInMvJM+wgm7hr1MfE/3AiAPcMfxKAHTIyNtovuQj3Kj8LTxTuaygCC6SmXIPAQRkUZx9bpvPaJbOkcpvp5epfl6hmjp38zLLU7xs4uDt6x/j76+VkfWyRLKnr+rzKoDTNoJMF5nb8Xg74vLtkBa0Ypn7mz9Ex8z+Ws5MOmo8uvjWd6POK1e8wSWNd2lH9WbaN9t3teGVn756vovUn5sm5X1OuRnMjsKqOn/YLAEYUaJZ+ZHtlvu+RpbEIrtXYsih/nqGgi9/0/RyAg3LmARtm7XHv8g6c+ckWXXUIqtKtiWufHlFZPYFzOjn4IKvnrEZV6eq7db774wup7X7hlh+FFcFChAgR4qeEth7imQrCB0CIECFCJMEB8dAJ3HpQHM/ik7XDeGH2dgB8NOJxgEpnkbaRuRoIWgUme0AhbJ9VCEDHL2Xur5VWFbEsmZWuc0ciq2Wax0s30AiNRVo/6e3fsP2rAKyLl23U9+oQ6KafmPcjx/9SuQHpPv5/97suA2DIYZI6WL1nfwDyfxSFlTFcafU3TjsMgG7tRCeNm9KXdjM1PnHvQ+z7TxXejr6vff43VKUg3x7xCACXzlUbO7YvBOD8gsLKPgZ0SGCCJxeODyiiwPlaQe0FwhORRtVC7ks8RdXZU1SB+b+4QvTDHfNVV2GXgW9UXvNunvra9RWJuw27WhIC1lU5IZe/8zIA+/o0kyNztP2hA972x2hXSW+MK9M5rPHSC6N8CUMe1Nvul50LwEF/k7je1XdMrNJ/QGmR9ULLUj/BdZwfE9Va5G+Sg7pq3AIqcauM4PdRNShh10z4YKvXgQ20zTp/7TtFq9a5yK+hD8vja0lPmn0HlOKKuK/N4O+juKewm6LORFVYWA8gRIgQITZHhBZAiBAhQmymCAvCtDIsX5bPkw8fxLo+MidPbKdIg1eGvF5JBcyOad0zK0cC8O0KKQjO+Uyp7xGfBtBrsszbrt94M3a2D/6ORjCv7Y6P1W9MNFClMmG5Ik36pyudvy7qJxldo7kM+eQMAAbeJpPX9tK6yUsUBXTudSoJeXSe8hjaR3TsC+YcDsAlPaTJvsOQ6EbRNIt8ycOnVitKachTojK6f6VjLZAKBXsepKiWch+9kWkb2gnotiDiKjhGYJqXe5M9FkSdJVnXo8vKKqOMAiTLSNQULTIsXdTC84Pe8cfeQJvk+fjxvXYUhTHxBNViyFqp6xpICEDVaxIcO068kt7YNl30yLhyUUKrfURKoND66D+lbh7IYjS/1FZV1JQ3EahxPrVaCqOvLPD5C1miaA7rPJ6CqLbZK1vx/8H9E0TqnNVedNaG6K6qUV4B3XPzsp2ZsFo5Bfnp+n090Fc5LPNjoiEDuYkgb2JKuca6Z5q27xFtV2OEXFPmmdSFpkwEM7PfA79GxsV4lBDbA3gG6ASMBk73iskthp++jRMiRIgQ9YTqAVhKr7pgZr2Ai4AdnXNbA1HgJOBfqC7KYGAlcHbznVENfWsreQBZg3q53v/8HbHFmm3dcJgccMnCUwB/XOwdxXfuCkCnCV4jf5aEvuJrvKM3uQoTbFTFa6PvSYgWFOiD1/6PDZfe/Kxj2pGrhEs6TC3zyzSr+tfBEqo6PndVlbZqqnAFGxzc+1yvgu2d76leVCySnV2lP8VHaCwW7KMbddDzJSzdTrO38v10/Pd3ug/YEP+fapZ1Yg7A4Pekr//8PvcAG3IGUhWHWx5fSwc/k25obYHq8hMC53NyDH9jkFyPILCIZlXoOudFdMxkYbWmxjnzlKvw4Zey3DKX6dwOPGYUANd1lzM6sFAakkUfIJjhZ3nhvSBf4h1fve6LBco/eXjbR9kyPagulppj9osS9Wtuhap+dYquoX9aEQDdvBVdFK96H2X4YIj5vpresyuVB9MrQ/8Hl271XqNi83sO7+DOfmbflLa9fttXaj2WfwB8DYwAVgOvAHcATwLdnXMVZrYbcI1z7pCG9rkhaDMUUIgQIUK0FOQETjkKqLOZjUr4fp9XMlZbzs03s/8Ac4D1wLuI8ilyrnKGNA/o1eiO1xPhAyBEiBAhkqCCMCk7gZfVYQF0AI5GBbKKgOeBQxvbx6ZAm3kAdMtezeUj3uOR2bsB0N2biNPKi7ms8FgAZhfJLO1wjxxFmbne5C0WfVKxZCl1Ion6ieYnmfFBOcsKPbjX7q2Sfjnz5ERLW6R+9f9fBHeNnL5n95EW/dX/k7RCn3TJOJQnsW8BlRDxvGIiXRLo/w89XeJus9ZqHNo//lWVNgIRuzQf575od53HD8feAsD4IzO5YMLJAGR5KYyP1ytPIaDTAmfdpbOPAWB4vpzkf+08qcqxiuOl7Pj4pQB0lHw9Wx+ofgZ01jr/npwfECDYDmBWhWiGbhGdd32d5YnURvA5YvWjO2qj4QIaJNdEOwSia0H5xAFpQZ2FmmeOyfr6DcXHJUbnDFGZH/5Czue+G8mjVC3c3hj6K1naI8iXODJH9SnwdSqSnemwQSZinhdHfHrVTgD8skD7LIrJadwuUur76Sqpn2Asu0SDsU0WBtR2O3T7ocryS1M9sVrQhFLPBwKznHNLAczsJWAPoMDM0rwV0BuYX0sbzYLQCRwiRIgQSZActKX0SgFzgF3NLMdUEP0A4EfgI+A4v80ZwKvNcjK1IHwAhAgRIkQ1iDtL6VUXfEXEF4DvUQhoBFU6vAK41Mymo1DQB5vvbKpHm6GAOkfLObv9fMZ2lrLgxXcoVr3LmFIylsoUtn+I0lh3QREAeTcr2dzWyFwNInZiRUVVvrvKmH+3oQ6Ap4Jiq3ykTlAWcIByCuYfKbO1x+erAajIV0TIjN8qJn34FnP5cz9Vutw1U8N8won/9WejtmKVJQ2rSlcsiam/QVROqSvn0dWKLpr0vLTpuz2uKKBonsz+oHaA89E/Sw+X8mmsi85tWoVM5e5p6zhviKJDXvN69gH1M7Fcx71ylkr07dlZipQn58u8XhlXG1Fvlt+6Ynu6jPHyCH3U95fXKi9hpyypWCbHbQd0QIAgxj7donSIBKOz6eYl1VE/QZ8D6ieQSRjq8w+C/iZHHgX7JeZENFWJzH2zHNtnqNJqfkT3wItr9b5Tpii73mnV96+lxze4xoPSNbaXdxoDwIKY7qO+vqzo1hnq3zpXwQJfZvO6+VLBXbRO51ZSoTZyM/w1Sdf9/ePibgBkpgc+1esa1WeHNWkmsHPub8DfkhbPBHZusoM0AG3mARAiRIgQLYnNXgvIzF5LoY0Vzrkzm6Y7tfQFzVxu9LHN2/WQqlb6JyW4qao01PUozQaCqkqzrpWzKWN7zdp7fK2ZW6TU68j31Aw76r2x2V9NIzpQDtGgOheF8svEt1Kcc3mGZsG9npAzNtB6X7S/rI2h98uBOuPKzowYomnM6njgWEzz717fvgYHZaBZHzhFD/rwIgY95GPLp2hWHsxzYmt0vIiPu1/2K8VDFw3z5zRVs6+//lV67hVz5rP0POVHvH7lv3VueKslXePx0hBZLsFsdb3LqNLvwFl9TofveDuiSmS5C3S8f085GIDnRzzot63q9Fzlx2JMmXTxd8+UhRXEqrcWJOY4JAvX1VSLInAgl3sLIMgTaK4Zd9B+gGO94B9U37/m6EcwTjO89fhtSV/eWDYS2GCtru2leyPmgzJuO+gJAHbJVNZxTz/zD+6vbMukk+/qQ/2UwR4EEwQicUFQxE3LZcV+P1lZ/+2+bRpROIdREQ+lILZE6cs1wYC7amvAzAqAB4CtUXjtr4ApwLNAf1QT+ATn3MYZXSFChAixiRCWhIQ/O+c+qW0DM7u2jjZuA952zh3nC8PnAH8CPnDO3WhmVwJXIodIiBAhQmxyBFFAP3XU+gBwzj1XVwO1bWNm7YG9gTP9tmVAmZkdDezrN3sU+JgUHwBBCb94DzmBZh7XnoETtC6aK4djrFhO4UE3ee1373zCx2u71TKVc75XfkBQvD0GsEKGSECp2JD+2sabnOmT5dwMnK74EpI9P5PZPf8AvZcXOrb30h7f7yY6pKbU+IAmCZzAgZNxfGl3QCUrMybNBqBi6bJq2wgKzQeTFpcms7v/o9KjqJgvGQy35wiKdtK27X1/Dpx4DABzxsuxffGhbwLwizwJgF23ULTO1wtlZl/j9d63yVzMJ/+RAbjcS1V094JtD62WE/pfY5XZfuZWXwNwVaepAByarT6UVlN3tali5RuDhtAlAVWU3UL/Gy3pzA2kSCaVixYZX6rC9zv5Ghs903TSXdNWc38/3T+5f/ywjlZ1rwTX+4cyvXePllbmyCyNiQ6dUyEn76T1SpZdF9fvc9sc/R7/s6/kVV7ZWtTwD5V5uA1HKAftYWZDgcuBfon7OOf2r2PXAcBS4GEzG4HSny8GujnnFvptFgHdajjuOcA5AH17hf7qECFCtAwUBbSZWwAJeB64B7gfP1GuR/vbAxc6574xs9sQ3VMJ55wzs2oV6byexn0AO47IahuqdSFChPhJIPQBbECFc+7uBrQ/D5jnEyFAyRBXAovNrIdzbqGZ9QCWpNpgoLBoy2Rud5jswKeDx9ZWjfc3T6lUqn+uW5dyxysplYlTqywPom8CiijSTWUU05YomqX326KTFu1dwKsnaMgyrWopvGSs97RWNqJk5lb4mGcffePSjHjx2pT63fl70VsZa3TMwjNF20RL9F68ZTlZM3Wcn09REmLhTBlgUf9ov+tFxV5vf2ohAO9NHqZ+rNPtctWYXwKw+6HjWFYiM/6RQSqtGOQSHN5OpSp/tddjOkZF9REq1dFim5L62ZwQRPCsjK/n1WLlmXyycosq2zzaT1F3gSzHTj7waKfMxX6LqkqjB2VXUJ0cRG0IrvcOGUG/0irprY6+XsO2GaJaf96uqoJuMg12bDu5LJ+qVw82hoMwCsjMOvqPr5vZecDLQGUmj3NuRW37O+cWmdlcM9vCOTeFDSnQP6LU5xvZRCnQIUKECFEjUszy3ZQws1Qkj9Y65+6taWVdU63R6GEYjMTlCescMDCFDlwIPOkjgGaiSjgR4DkzOxuYDZyQQjshQoQI0SIICsK0clwO3M1G9fWq4FygYQ8A59wAADPLcs6VJK4zs5TsPOfcGKA6qdQDUtm/Jkz1sgrlJ8TYetsLAej7jqgUt95H1YxWsla1hV8aiKDMY7xcx4h4asZ5CYmILynZ872l7LvDxQCMPvBOYGNFxQDZngaZURGUGhSt8ue39FwcNn4+FSmegxulkpB5o2UaF3g1U/MRUEUHDeGuG6UMuiim5LWLRv0KgE7j1EbuPF3qX8fPB+D8496ucoyHnpCS7SfTh1QuW9Dfv1eozcFZ6m9ADfSKtq5Er1Rx4swDeGSAolqSy1w2BVoi4ik4xrvrdV/dME0UX16mrvNtg57juNxCAM7MV2xGfSOMmjIiKbGtTZkg2NotAOBx59zfa9vArHb+OdWrVl35qepLUoUIESJEG0dQEKYpxOCarY/O/bGx29TlA+iOqtRkm9l2bDA18kmuCt3CCGYJ6QbxdAUIZY0pBDbEygfJ/MGsPRBKawgCeQkbISeZ+VKarlDCW9ZJtQgokxVS1rcjXT6SV2uvfCVTf7Lz/QB0ilR9KI8pU1unPX0ZAO22kWul3RydY8XsufXvb0SXKhCzW/+upCwWTzC29Q7skcjRPfa026rsG8g1/H7OUQDc+fmBAPxxb82Ey6R+Qd/H0vjDXUrrP/67cwCIjJbFce9vlB+wm7cTo/XU5W8teHbgB0Bmnds1FM0585/jHe9/nn8EAEd5EbZPtn02qQ+1BylsjpAURNu6Z81sV+Aa5IW/zTn3cl371HX3HYKSuHoDN7HhAbAaZfOGCBEixE8Srd0HYGbdnXOLEhZdCvwc/U9/g4J2akVdPoBHzexx4GTn3JON6WyIECFCtBm4NuEDuMfMvgf+7X20RajATBxN0utEnfancy5uZr9HFexbHSJEmHG8nNzx40X6DHn+d3r/w2hgA/XTGCookIuIzJhXdbmXhIitTIpPnllIpy1VLjKWpfKMDwyVcuGlHSYDG2iRz9Zpuwwvh5f9ZAEA+dN0DTfoUtajv3HRSksu2F3H2EqO39yts0h2/SSrXV65eE8Aru4tVdCeSiHgwaJtAOg0Xm2nrynnb/93FgDdj5Nq6pzhuqVGZIhe+t+6jv674sY3LlvY9Ch3Fa0ql6CldPiD4wRqrd2i4t8e75cs59V6xqa1op5F4TcJnHPHmNmRwBtm9hhwCXAKouePSaWNVO/I983sMjPrY2Ydg1dDOh0iRIgQbQGt3QkM4Jx7HVH17RHlM9U5d3tQf7gupPoAOBE4H/gU5QaMBkbVv7shQoQI0foRaAG15geAmR1lZh8BbwMT0P/00Wb2jJkNSqWNlGzBIB+gtWPfCSp6MvTKsQA4X5gjKO/ovFc/iOgJaJ36ICgnmQriUyWHkDuoAIBHJ6sQS5etFZ1xYI6KuzzwjOLqe/uCNaUd1L/46An17l8l/Dl3vVPRurMv1/ctI/GNqIhkiuLGHp9rnwpRPWPLvNzDVBWb6VIkimHWUTnkSaSUwvmSxHh87/urtN0/bTkAM31+QO807RsUT0mmnxqDQLEykC3Y1EiO8S+OlzRL35KvX2YbjbhqbYi1fjXQ61FJyWzgHefczsAfzGwIcANwUl0NpKoGmg78Dkk7g+Sb73XO1f8fNESIECFaOVzbcAKvAn6BOP9KPTXn3DRS+POH1L1BdwPpQFDV/HS/rLZqYS2GYKZVVqHTyctVXHNshRyRgTgcvSV6FhSJr5gjx2UQMx/t3ImKJV5v3zXE9VoVgbM558tpAHSJytl7/VFHAvD8QOUQdP/Wl7tbUKR+ZXcCIJLthbYakckc9Y7oGxcUAPBA3/c3miEGM8dgHINZee+oLIBCX4g78oUSAAqP9jPb5bBqD/Vtn8GydjJM53zA2NMAuGDwxwAc0U51Cc6YfRgAf+6pnIJh6U1nAbSWmX+AZEd01IzVcY1XU2a4tnSR980FrvU/AH4OnAyUI+dvvZHqA2An59yIhO8fmtnYhhwwRIgQIVo/Nr2Dty4455YBdzSmjVSnDrFEp4KZDaR+dQFChAgRok3BOUvptangcwAatU2qFsDlwEdmNhNlmfVDqp4tjmSH1/iy9fzyJqmidrv7OwBigdPXUzCVjlv/XnGgtOmW/aIPAJ1+FAWzdNt0KtoNBiBLvkt6vlSofeYvaHCf3VqZ/e0+U22BLh2krz+lQuXthn3xo47h6xlkV4hiiTdCuiIoj7lkb+UgPNnnUQDSq9HwC8Y0oCwCB+2+Xvd/zVrt03usKLVeH6qfU8/M59N9NAHJMeVYbP/2RQD8cqevACh3Wh4I4T3U9wMAZgSFFVoILRWLn3is4D6s8HOlmHObVNwsROpoC3kAwJZmNq6W9YbCQ2tEqlFAH3jPclAtYopzrrS2fUKECBGizaJtFIUflsI2tc4i65MSuAPQ3+8z0sxwzj1Wj/1DhAgRok3A0fqdwM652Y1tI9Uw0MeBQcAYNjxRHNDiD4DAhD9vvmLSZx/fja6FXpnaSz1Ypo8s8cqcyfH+GV9MBKD3vN4ArBskJc++D00htmp1lX2agqmIlyo+Pa2/jlcwRRRKziKpcsaKi6ts3xD1z2QEbXa6V2PT4W8bi7cGUT8rffx856i2uXi+on3TIqIw/rujVEBuePBMAEq7KsrK5VewLu5VR330T/pSRQw98YnkJK455EUAvivVeG6Zru2GpauNIHY/09KaVb6hJaifoOxlfy93EfHRVunBz6wJ/09am9zFTw+t3wncFEj1DtoR2Mo5FxZmDxEixGaBzeHfLtUHwASgO7CwGftSKxyOclfBF6WaYc78napRucLxG7bxAmiujuLvQVx9fLLi8zMm+/23HIpbvrxJ+10FphlF2gpVEbNvlenbnPdZ2sD+ACyPyxLoFGlXOfNf7mffnXz8/N+WbAvA3HUFAHw24hkAlvjY9WvveRCAMz49G4CMnDJmVEgSaqnP9E0fKguqx6OyJq7vdjgA5ct0jO+OlijdfSuVXL6qQtv9rcv4Sufpc8XKgzgpd2XjTr6F0b8FhO4CVBDbYFmEaBa0dgoIwMyiwPvOuf0asn+qd1Bn4Ecz+5aqReGPSqGDhcAaRB1VOOd2NLNrgN8AgWDRn5xzb9aj3yFChAjRbHAOYm2gIIxzLmZmcTNr75xbVfceVZHqA+Ca+jachP180kIibnHO/aeR7YYIESJEs6ANUUDFwHgzew9YGyx0zl1U146phoEmC4pXgZl95ZzbLZW2GoqJy7uy9aMXMvBa5TW40vEbb+TjriO+5GG8rKxex4hNmtq4TtaBCk85tSRiHUVLnDbteABeG/papfMwiN2/fpmSvF98aS8AbjtDgm7Bdv9cfAAAb07cGoCOneXszEyr4I2VIwGY72mjWEyzpsfvvBmA/V9Smcs7jlAeQr6XmfhhdV9t57XqJ5aXcPo4pZa8OfJB33vlMgTO1W5RXdemFJBrq2jKwvQhqkdboIA8XvKveqOpSMTaRFgc8K6ZOSQgd59ffoGZ/RLJSv/BObcR4Wtm5wDnAKQVdGiiroYIESJE7XBs2izf+sBXbswG+jrnptRn36YiuWozlvZ0zm0PHAacb2Z7IyG5QcBI5Fi+qdpGnbvPObejc27HaLuwcHWIECFaDi7F16aGrwo2BtUFwMxGmtlrqezb7GEEzrn5/n2Jmb0M7Oyc+zRYb2b3A2/U1U7m8goGPb4M80qfQZx+deUd60v9/JQRXSH6hFNFj23xj3MYd6BEXQNZggPzlBdx3K9U4yfLx/Qv8kN7VmfVB7jzoG8B+Pn0g/V9wIvkeZro5OmqxZD2g2ib77brCcBlh+jS/mnCzwF4t7dqIBzbWcd6e71onUOz4fkRon4eLtoegAs7KMv9saKdAbiik74nlz5sSjqkJSUjGoPW3r82D9e0FJCZFQAPAFurdX4FTAGeRQm2hcAJ1TEhKeAaVBfgYwDn3Biv11YnmuouqnakzKydmeUFn4GDgQlm1iNhs5+jMNMQIUKEaDVwcUvplSJuA952zg0DRgCTgCuBD5xzQ4AP/PeGoLyaCKCU9OxTzQRuB6z3BeKHIg2KtxIKwpxew67dgJdN8e9pwFPOubfN7HEzG4mehIXAb+vsRGkZbuYcYj6rNhGRTLkg4tWs29yxbsuuACz5lXIjCjLWUhokc/s8gL2yNH7L47pnSnz4Q7rPZB1bItG87TKUBvLi4LcBiLlsdvn+VAAynpePZt2+sr7ivppS97QiAD7b8SEARrx9IQA9d9VE59IOSsI4tfBATu+mXIXD8jTTf2Od5gknt5fIXxzN9M+YfSAA+wX75s0DNlgCcTauepYqwpl1iABNFQVkZu1RMa0z1a4rA8rM7GhgX7/Zo2gGf0UDDjHRzE4Bol6z7SLgy1R2TJUC+hTYy8w6AO8C36H6k6cCOOeqncE752aip13y8poeGCFChAixyVFPLaDOZpZYI/2+hGAXgAEo5+lhMxuBaqpfDHRzzgXJtYvQhLkhuBD4M8rRegp4B5WLrBOpPgDMObfOzM4G/uuc+7eZjWlIT0OECBGi1cMBqT8AljnndqxlfRqwPXChc+4bM7uNJLrHOed8pGRDMMw592f0EKgXUn4AmNluaMZ/tl8Wre/BGgPnXI0UT6OpH091RLIyiZf4ROcmKAnZGpDxPzlu8zopTeOjf97GCk/1fFcm2qZXVPThF+tVPvLsfIkMBjr2Z+ZrknJXUX8AzspXvkQcx3vbidp5e6hooucW7QTAsV7GYVaFqKfZFbpd0op0y71w60EAPH64HLwjui/g71NVKvPPQ/4HwDBPOS2Pi6Lq789pjwI5knumV+8vawiNM61c4nlD0nPrvW+InyaaMBFsHjDPOfeN//4CegAsNrMezrmF3i+6pMYWasdNZtbdt/tsTYxMdUj1l3IJcBXwsnNuovcwf1TvboYIESJEm0BqDuBUnMDOuUXAXDML6qkcAPwIvAac4ZedAbzakJ56HaD9EM10r5mNN7OrU9m3PpnAn5hZjv8+EzkaQoQIEeKniaYN8r8QeNLMMoCZqKJiBHjOU+uzgRMa2rh/yNxuZh8BfwT+Sgp+gFSjgHYDHkS5+X29I+O3zrnzGtrhVgVP9wQqoQBpXVRKsWJpVQmjSE6O37akyr6tHR3fUMTML8afyQGPixa6tKOkKc6dJ6mHK7q/C8Cv5vwMgM8nSHH1+8NuB2CbLNUpmOdzL7IsXqmA2TVNKqBZUSmNTi3XWC6PabzWONE41xz5HAB/6aC8gfxo3O9XTrkX3/rbZGkMHtZbpTKnr1Uk0596iRr6Ra6SHefGlEOQaV76owEx/DMqRP28ukaxCltmzQfgZzk104or46K1gjKXIX6CaOI8AOfcGCSrn4wDGtu2mW2JgnKOA5ah3II/pLJvqj6AW4FDkMmCc26sz+gNESJEiJ8mWkOab2p4CHgGONg5V6/i5SlnAjvn5vp4/gANr1jeBlCxbEW1y+NBrYFdNWOMTpwJQGzNmirbWVo60W5d1NaCRVq4Ca2F2Ao5TNOiUT5YKiqyX4asm4MKlAlcENH17Ziuc3z4AGXmzq7QjLpPms4xywcr3LZsb0Zfrqzd+XtrFn7NSaoh0MWHCBRWaOb//DI5h/doLwfuzMMe2LiT8iNXFqXf7R9iGTOLdLzTu6jfO5ysPIHLvcUCOnZDnL+Ty2TpXdxBVkUqVbaCmX+QEX1FHymZ75oZ6vP/tNBmtIB289TSUDPriGq2l9e1H6T+AJhrZrsDzszSUQzrpIZ1N0SIECHaANqIBWBm+6DyvIXoqdXHzM5IlNypCak+AM5Fqcy9gPkoGez8BvU2RIgQIVo7HJC6zMOmxs2I/pkC4NUangZ2qGvHVKOAluGzftsSLE3SANHeEiarLLaeChVT1zZfjwU25sEsuiE9omL+gqr96KDyiZX00iaghCqWLlPAGPCXv58CwKdn/R8Aty/fBYACTwGNXd8PgHML5CyeVaGzfbJIdE6EOOlFypso8/TN7TP2B+C+NFmgwwoWA/CLTqMBuGfevgAcl/sysEGQLhGlvmTl61f+G4CsyuLqes/0NE26bbxvXVgUk9P3jGknAvDi0Bd9W7UpmldF4AReXaa6BNcWymn91hZhUbufEtpQQZj0RBlo59xUz9TUiZRIUzMbamYfmNkE/33bVONMQ4QIEaJNoq3oQcMoM3vAzPb1rwdQnZU6karX7H6UCFYO4JwbB5zUoK6GCBEiRFuAs9Remx6/Q4llF/nXBL+sTqTqA8hxzn2bFAVUUZ8ebgoENQMqCmdXuz6gZoLtmuSYlfUJNpBDlf3w1E8kS9RBYt4BbKCPqqtx0Bzo91cJBh4zSSHDiw7WJf3jrlL7PLf9fL+lxunyWZJq2KmjxnP8qp6s7af6DHsMU55BXpri5//bS1nvr69TxExA/dwz6Hnfpm69+TFFFvWK5lXG8QcUT0APlXtKaFy51m+brvswlaifgPJ5drXKWa6LK2Jo1jcqSZm5RVqVY9QWBbTeie66dbkosFlTpFaatlrXbcjCMwHI+U79Hnv53XX2L0TrRYOVeVoIZtYF6OKc+xH5AW72y4cD+SgzuFakagEsM7NBeIPHzI5DlbxChAgR4qeHVOmfTfuQuAPoXM3yjihop06YS8HT4bV/7gN2B1YCs4BTnXPVT62bAfnW0e1ijU6aq4K0QQMAcCuLAB8r7x2OzeGgratuQUOL2TcVolsrzn7W8Z3Uj2FrAZi01yMATCzTDPiFVUpoPLngW469X0Xf8+boPur0hXIe4u01859+uayHx3dVTsF4X1vgN+01fwhm/eviZeREMqr0p6Ha/Mvjazln1jHqT6bO4foeyhl4eNVIAM5qPwbYUPMg6mO+A6vj5pW6N56841CO+J2i6S73juzDfxT7eWF/yWHlRWXJvVO0DQAdvBP9nQVbAvD+Nk8DYTH7lkS0x7TRdSh01orMfn1cj6suTmnb2b+7vFHHaijMbFRNxzWzCc65retqo04KyMyiwHnOuQN9YZiIc25NXfuFCBEiRJtGK6eAgLxa1jVNFJBzLgbs6T+vDf/8Q4QIsVmg9VNA083s8OSFZnYYEpyrE6k6gX/wVeafB9YGC51zL6W4f6NhkQiR7JwNUgxNgIoZszZe2ITUT6VTN667ZCPqJ6hDEFA/m7ikZWyCQon7TfaThx1EYQyfpJy/jCIt/s+F9wPw0Io9Kd1S9Ee0RNRJBy97EUnvDUD356Wvf9pitTH4Kd0+2z+nOgJbpGm8F8Yq6Ogd55fNOxSATyarPsGb+0mM7pwpSkU5ttcYAE7Ol+x5lmmcV8TlyO0cSefRAa8D8EKxnL3flEpQLihBOcbH8KdbVYf7A8tUn+Cbu5VDs8u53/Pcy5K9Wn6EHN4Lx3UH4A6fUDF/keoqpM0VxXf1cXJ037rFuCptN6ZUZYgWRv0KwmwqXAL8z8xOQFXGQIJzuwFHpNJAqg+ALGA5sH/CMgfU+QAws0JgDQqLqXDO7ej1Kp5FNT4KgROcc9VX9wgRIkSITYDWHgXknJtmZtsApwAB3/8JUmpOaTaZaibwWQ3rYiX289nEAa4EPnDO3WhmV/rvDSmGHCJEiBDNg1b+AABwzpUCDzd0/1TrAdxezeJVwCjnXEOq2BwN7Os/Pwp8TB0PgMFbF/O/d75gtz+pBEHBo766Wn0om2aM8KkVNR3PL2+puP9UEeQtRKfMAaB/iaQ0Iov0DL/1dVmXA55dwJh97gHggILTtO9nqiHgApXU7RRRNPRe5UAsPFBRa8e9K0po/OF3ALAklsM9yxXQ8NWc/gCMOfBOAKZViJK6YfArAFw26TgAPs4TRfTbXh8DsE3GcgAOmXAKC+bouLOOEF31danoocIKTYxO+vxCAJ7ZU7W7B6frnN8cJZVXdtA1+XZxP9p5gd3p5+l4A6OiIUs/7uLHR9dxhRgzTs/TOE0rV9RUhEz/HtI/bQmt3QJoCqR6R2YBI4Fp/rUt0Bs428xurWNfB7xrZqPN7By/rJtzLsgjWAR0q25HMzvHzEaZ2ahly1vXn2SIECF+4mg7mcANRqo+gG2BPXxEEGZ2N/AZig4aX8e+ezrn5ptZV+A9M5ucuNI558yqf9Y65+5D+QfsMCJrM3gehwgRolVg00f4tAhSfQB0QOUgV/nv7YCOzrmYmZXWtqNzbr5/X2JmLwM7A4vNrIdzbqGZ9QCW1NWBCSu6sMWT5/HjDaIFDlouqYucwtUwV5EnsaKi2hsJqJgWooJSpXbqI0URyVa0TaSjIk/iPomtKaOjKvsVtJmm8ZpxwWAAdjtI0Tc75hbyu7mHAPDE1o8A8H93qEjKh18qKarTGM2Qluwt6sf5obeYlk/3giLbZJSyR/fvAbjJv7+yVvsMSvfuo4gS5G7bSkVnOkWq+rnyvIzDwsJOPH2IZBjOnrsPAP/upUSwTH/L37SrSlN+vX4QAPcsUYJah7Fav3Y/SUj0yF3NlVfqeOePl3pqryvU6dLOonbW9BVFtU7KEJUFbfqnKSrotpVKsHv21oP59NpbgTAprE2glT8AzGw8tfTSObdtXW2k+gD4NzDGzD5GBQf2Bv7hE8Per6WDlYlj/vPBwN9RackzgBv9e0P8CCFChAjRbGgDPoAg1DOozfK4f09Zuj8lKQgAP1Pf2X/9LpXak15C4mX/NQ14yjl3g5l1Ap4D+gKzURho9TUYPbbZNsO99GZnBqRJYiBwqJW7Cna7Vg69bp/IkKiY6nMgWtDZG8T8R/J8cl5FBbHi4ka1mSgNkdZDsefxFUU63hBp9Qex+82CQCZhuBy700+VY/XBEzW73ipjLSV+jMeXad32GbqM08p1naKm9UVxfZ+wXjPtp+6WpdDvxBkALFiTT8kHcqpe8JtXAOiYpvEbmq7rWuI0XwkctqNK830/FEHcI6o4/ReLO1SewkPz91DfF6ntq7aT0N3u2bpHDntX6f6dvtUsvttrKlnpyr1VVlZObFtZPmXttU3O2HnapkQWSMmOsiJW99P1OvA8iexd3kXvz6+R8/j2p48iR+URuOqyJwA4tl3D8iqfL27PgnKd58UdChvUxk8ZjZaC6NPH9b749yltO/PyP2wSKYgAZvaDc267pGXfO+e2r2vfVKOADFWvH+ic+7uZ9TWznZ1z39a2n3NuJjCimuXLfXshQoQI0epgrk1YAAHMzPZwzn3hv+xOigE+qVJA/wXiKBHs7yix60Vgp/r3NUSIECHaANpOhM/ZwENm1h5R9CuBX6WyY6pqoN8757ZPNDXMbKxzbqPZfXNhhxFZ7pt3eldSP595/98eWRvS6/+4WFbQ+D286mYzOEY3BaJbb4GV+JoC01OS+Gia47ZvD2ygQyLtROMQ1IXIy2X4CxKEvdrTHadO/zkAf+0nKYaevj5AYbnomYAKuvAT5Q3cu+8jABRE1vPE8t0BWFCi4w7JFfXz8iuicQbeNVX9KfP0TD/lJ6zaUtuvGqj7oNPECub+XNRT2hLRNnmF6nPRbupPnx6iqnIuEW0z5bcq1znsrsR8RWB9CWTrfoq1lwO+op32Ke2g+VNaqX5DsQwdY9k2ogNtxGoAPt1FuQY7v3QpHxzzHwB+fYaop3effACof47AjIpiTrz+cgC+veauBrVRHwS1EJKd1xPL9Ru7a4lkMe7s9VWryHdoLAWU1buP633hpSltO+PKSzcpBRTAPwBwzq2qa9sAqVoA5V4VNKgH0AVo+YK2IUKECNFCaCsUkJn9Nek7AM65v9e1b6qP6tuRM7ermd0AfA78o37dDBEiRIg2hNavBhpgbcIrBhyGdNbqRKpaQE+a2WjkuDXgGOfcpAZ1tcFwxFycGRUyOX/zhMzftHXw3YW3AvDvbj8A8LPuRwMQn1lYa4tBgZYgDr+1STIEaNZIn9qOu1oRKmk9FYHkCqTsGftRkTLR9SWM31V0yEnpBwIw6WZF+aAgJX4sU6RKl6giegb5YiljDpEExAFjfwnAjcNe4pae31Q5/tb3KLotx2eJHP+Zcg5vfO5Y9cPnBnb/Ru/9HvH0WGYmg1YrKintS4kk2jBF6qwcrsih4s9EH628oQiAvI80F3KLfBW9PjrnhQf3pet3UjBdM0AUkItohpWxRvfL4h30M+r9oc6tm48oKp2m8SraSf07a99P+HCdIopmnKrjrYsrtyE3kkV9UFjeno4/6njvrte+h2Y3XyGhgPoJSmeW+vfh6aL0buqpojmjSyOsjuv8D8hunb+nlODA2gjH4Zy7KfG7mf0HeCeVfWt9AHjVzgBLgKcT19UVuhkiRIgQbRatY3bfEOQgqZ46UasT2MxmoWEwFLO/0n8uAOY45wY0tqepYscRWe7bd/pWu27L+5UV3ONLzeQzlsvRF13gHXrlVYuy15QfEMnMIjJAM9j4HBVEr9ORbEksWksLzTUA0Vw/k/d5CpXZxbly1LoSOfxia5Ji1P25mp8BJ1pMkRyfn9FeM+wgXyHSXnkRrruyehfvIYvguav+D4CsBKK1vWnmGBSFXxnXdcyLaPmZs1T7YsZjiqvv+pWuZ0kPnc+qgbJG1nWH/i/r+PGxP6rPaelV+m6+v6U7akY+8yQt33FYIQA/FOr3kzElm7KO6mPWUm0T837Qvm/LMlg1JKfKMO37+68A+HpZfwAWjJa18bODvqVPlvr8yTKdw3OD3/DnXHsBp6B0ZpBlvCZeziH/JyvYJ0gz+q8tV4R+UUz3zu7vy5k9+qA7fP/idI/mtlg/akKjncC9+rh+56bmBJ76103rBE7KCI4CXYC/O+furGvfWi2A4A/ezO4HXnbOvem/HwYc04g+hwgRIkSIpkFi8ZcKYLFznqOrA6k6gXcN/vwBnHNvoQLxIUKECPHTRBtxAjvnZiNW5kjg58BWqe6bahjoAjO7GnjCfz8VqFMKoqVwy2kPAnDRFicBMPAWmerx5ZIIiO8wDAC3hcx6+3xMte3ES0uIT55W67GC2PjYcLFfafOlQY+n0tbs0AuAvM9n1Ek5NQQBXRNfr1KMiXIRtSGtq6QQ4qvWYJ1Ew0SDlRmiH4L+Rj2NE9AmkSAO3lNCrhrfXkCVJVNm8aW+X0vkXO2+Wt7hQwZfpjbzNVGxighbbjkXgOO7y3F7dG4hABPKNE+5urfokomXiVKZWiJH7WcX7QpA7kJ1LH1dhKW7FADQrUQUT3yubldX6jX6czSOZe31E+jxgbq5qI8oq3iFjhnLgoJJup9W7asxz8vT+/zSjr7v2jd7qe6Bl9/U3OjmE1Wn45ncXdTPe3fm4T/dDMBxeSoXmWm11fXeQLXkelqs1JfNbB/JpP9xktEoO0fXa49jfgHAF9s2f6XWgObp+LXuv/aH6B6pTw7AYVNE6b21xZt1bLkJ0IacwGZ2MfAbNlRofNLM7nPO3VHXvqk+AE4G/oZCQR3wqV8WIkSIED9NtILZfYo4G9jFObcWwMz+BXwFNM0DwEf7XNyYHoYIESJEW4HRdhLBUHcT7fKYX1Yn6goDvcY5d01jt2kKOBzlroJ027jL+2UpGmP8XqKCMvcWdbHVl5IbaP+KzP3V/TUmv7xH23+ymyiEgNqwaLTuXIC47MLIelEbZQNErUQ/V4x69hwpRVawIdom7qNq6qP7X+PhPfUTqI+W7ivJ7+yxKt8483eiPMoHa7usCTr3Pu9JliC+dDnx2XOrtBn1tQUC1U9b42mcYo1TzL+nDRLttWJXjVvmqhjZH4yv0q+N4OmvIOcCT1UN/buicyq2UZvTf53GNf1e07b+l7fOc00FEd+Gv6XXxUU7bJklWueTqxXttep1UXxrBsUZ+ILXClkmGtCyFLpTsvdwADKK1I/ydmq00zeiqKYdLUomd4K2r9h1DQXv63jd3tFxlh6g45T+TGO6brGip7JW6JoMfK4IgEvSzwQg1kkcUbsOcP4lmkfNO0rn9t3BqrYa8xRi16QImoBq+W+RjvnRCtGZ38/qQ+73GtOuXTX2BefruCs/0fX7cL0K7TVUcTQVFO2tYz+yWsUQFpR14PyOqufQIZJT434AF/T5UO8LJDJ8e8+vgVZUOrPtPAAeBr7x9VZAAToPprJjXRbAr81sdS3rDTgJuCaVg4UIESJEm0AzqIF6OZ1RwHzn3BFmNgB4BugEjAZOd87VO5vPOXezr9Wyp190lnPuh5T6VEcewN9SaKM4OROtOVBbHsDquGYhQfz4Mh8/3iuq2dywB5Qn0OcDLU9fIsearda7W+WfcdEocR8bn2wJVMbK5/lZfS/N/ANhtPj3EwFI66QM1NignkTWaOYfnypnXeBUDWbFlY7belQoC+oCrN69PwAlHbRvl2/kwF2xnWbzl1z9LADXjvsZAAP+or7YuvVUJFkANSGwMoJ3F9e9Elgyke2HE1koJ3jFoiUpnwNsGIugLYtGK9ufc42cup13VaW3G4e8CMAapxnv0go5PQ/KKQRgYpkc8xc+qpLTOQshf47azZ4oK6Fi3vwq5xItKFBHvGVQeQ8M0Ew7Pslfs4hV9qsy/8F/D5zj8S37A1DWQd8zl2oGHllSBMDa7RQYECmNs3SErIlOk2UVDLpaltCDfb4ANsT5TyjXMTpFdN1O/72c5vmf+2zn/DzwfV5+qKy+Tp/qXKefo+ONO0PWRbr5WhXNMLO+o0i/yQsLZIEOfPUc3jlcju4h6bXnAwS5DfMqZGHOrZAltX1meaMrpjU2DyC7Rx834Fep5QFM+kdqeQBmdimwI5DvHwDPAS85554xs3uAsc65lJM5khJ1N0Iqibp15QFcm2pnQoQIEeKnhKaMAjKz3sDPgBuAS32Nlf2BU/wmjyImpT7ZfKPZkKgLG0gr858H1tVAqlFAIUKECLF5IXUKqLOZjUr4fp9z7r6kbW4F/ggEcb+dgKKEhK15QK96da8JlBh+Eg+AbJ9GP7FcY3nc5yoR2e8Rmb4/PHorAH87UtTCh//Ve3k7yRN0nKb9sr+Yim2rAt7uhx+rHCOgDpYfInM7d57om7RPxwAb4vFjAxSjPuuYPAbfKgpj6p2yDofdKYckEV9q0TuUbYm31AJKYaloleocqzEfT5/3vtZF91Z/l96otrpcVQRA9zS9T9z9MQBufl4O3t1ypnHOD6cDsF8/5Tz8sEy0x+LlolYGnTZWYxDQYJ6iqqR+vEPXChfg/JgnO3vjpVULticjkfoBsMxMon11//e/UfRl4Dy/ocP+AEy9Wg7QKw5TCen9vzkXgIJ2GotfnaDC7y/M2Y7iF3Rtd/2XnPI/7JJdpV/x3nKQxsdNrtJ/Jkz1p+ypvbINVGBlrv2WknFYN6gAgJw5crJmLRR9WDTc54psq/Esy1dbqwfF6f6VjhMIyi08W+e89xDRV/P21/JeQ0WpzV8kSq+Tp/ryO6jtyX/OY8YBCvse+qh33sd17w2+c5baOk3j1zfqY/STZUuaAGfkS6jw9wtFP3f+LsLxUyRR0euYQgD+N/StKvsElG1+RNekm+9fj6jGJr2R9E+ToH5JXstqo4DM7AhgiXNutJnt2+i+bWj3ECDPOfdC0vJjgdXOuffqaqNFHgDVOD8eAfYBgsIFZzrnxrREX0KECBEiFTShE3gP4CgzOxzIAvKB24ACM0vzVkBvYH492/0r1UvyfAK8DtT5AEhpSmBmQ83sAzOb4L9v6zODU8XFQLJ89OXOuZH+NaYebYUIESJE86OJpCCcc1c553o75/qjqMkPnXOnAh8Bx/nNzgBerWcPM51zS6s53jKgXSoNpGoB3A9cDtzrDzDOzJ4Crq9rx2TnR4rHqxeC3IDBaaIjPt9HInjn91N5wr3+fgkA0RJdrbgfmpculSLl0phM0Ot2PBjzkRvmY/grtpN5PXtvmavdv6kapbXsbKX5d3tbURCzD/JSAhmO/Jdl0g66Svu4DPVz0V4FAKzZRaZwp3cVOZQ3R9vFtlJM9ZxTRD8cu80Y3npatFWf/0pCwDJFOWUuk5mfeYQ0+mOeyvjnQOUH/P0QlW0eep2ilE5rP5YLt/oYgLPyCwFY3FW0SO80DcyyOUEcuZRRHz9IJRkrqalA7iFBfiKI6gkim4IoGxds46PNXJJkRbSTD2TIz6vcprLNIOpmvfrX8zOd221LjwEgx6twLNpF12/AYNEmh/eeyFdP6Pq9F1HfC/xYZ0zUdWKh/92k65rEk6is6uQuAppv1Qj1OZ7mlUWdv+ZeW6P8NFF6y5eKAhoxUFFXjw98jQO+vwSA1X009yoaoGu/7YmiHI/KE2V1fgd9z95GdMjFw3Ud3+shpsGtrWBauSinqWfIb/jeep3LrR+qNsPTRdr26s6TNz6ZJkLE+x8n76nr3zlzEu3e0OfnB75f7T53F20DwPkFEwDIifgSmz4CanJZOU+u1P3+/PjtAZhx4MPN0f1a0QKJYFcAz5jZ9cAPpBi7n4D8BAuiEmaWDmSn0kCqpGCOc+7bpGUpqc2xwfmR7FO/wczGmdktZtWTfmZ2jpmNMrNRS5e34eISIUKEaFtw6B8rlVd9mnXuY+fcEf7zTOfczs65wc65453zhZdTx0vA/WZWOds3s1zgHjboAtWKVC2AZWY2iA01gY8DFta1Uy3Oj6uARUAGcB96Em5Uv9J70u8D5QHUdbxgJpHpdFpPDJR42N0XyXm4pEwzsrdmSyxvTKlm2kGm5IJTtuCoX6uy0fcrNfudPFvPyKf2lqzGuWvlYC6Rj5HOY/wdkKlnWLuF6ubTv72Fc6+8BICYz0AuLfACZOoGvbvKKXzrtYrZ/+Wdvwegx+eKix7sA8LGlQ+nT4mmu2W7yOkbT1O/0t/5DoC0njoX0jQNXXmvZsXvbiNrKOrzFbItt3Lm/+AqBRH8Ml/WQxAn3jmqDM6j28mJvcOnzwAwsUyO04KoLIBFFQXcN3svAFb8T87Mno+J6Yut8A7vGhDkVbhSbx21z4FJinFPdn47/z3vG83es1/WrZfWRRehpJOcsn/rcCQAa1fkMAw5srs8I4tp6rWadaZfJEfpwPMXA3I+A1WsGa3YODcjsG4KxuharB0iB+3is9S/0vWa+Q64WRbBoTepwtlb83S/ZVoav71KyZo3Pybhtut/+ygAF314KgD77a3ZenIc/O+6fAzAN7N3AKDrGHhhD82O55aoH3PX6n3+8bquXx2r6zjnA90jfdNqF55rCIL7KqgHUbFiJaOnqV+FffW76uzrOQRVz35XoOzxbV++BIBJv9A9Gpzzz1/5La6Txvqo4eOavM+pwEhRS2HT4mrEwsw2s9l+WV9kSfwllQZSfQCcj/6Ih5nZfGAWcFoK+23k/DCzJ5xzwb6lZvYwcFmK/QgRIkSIlkErl4Lw1M+VZnYtMNgvnu6cq0GXZWOkKgY3EzjQmxoR51xK4iLOuavQbB9vAVzmnDvNzHo45xb6ZIhjgAmpdjhEiBAhWgJtRQzO/+GPb8i+dUlB1Oq0dc7dnPKBNjwAjjCzD1HZMgPGAOc654pr2782KYiaEKSar/Rxx0U+hX+d99adN0VJeGf3Uxr+XjkzifrHfreoLy/o/SuvF8usPjlPkVr/Wi4n66Mf7Q1ARpGP057pU/hf+bGytGLE6+9P/X1/AKac+l8AfjNX9MlOnpK5ZdwBAJSvlCkcyZVTLLYunbwpelbveqKojfnHiEea9A9RGq/vJzN6RUymeHdfhL23d87WllofjFOyVEBQADxwsgfbFcd1XrmRTD4t0VjOL9c57potmuaEfyoWvPPdX9Z43KZCkEsQIJKbW0kbBbINCy6REzV7mZdx8B6slUNl6Pe/9ptq27LcdsRWrqqyLK2rqKfiXfoDkHupHLdT5okiG3q+qCznS5Gufkn03PLV7RjeQ7Ta6d01LiMz9b1/HfTMTj8cD0DHG3R902YtYua5SvQs7afrsccWkq/Yrb3eV/l74dkHdF+NuaLuJNOa7oW6cPjBJwBgi5fjynyOR4Hu0de/eq1Km8F9NT+ma3TAKyIAxv7iNgAOmXAKy1eL1t6mp+i+h/uLzg1opLrQWCmInG593JCTUotZGXf7pi0J2RjUZQEEd+UWwE7Aa/77kUCyU7hWOOc+Bj72n/evz74hQoQI0aJoQwVhGoOUtIDM7FNg+4D6MbNrgP81e+9ChAgRYlOhlVNAZrZ9beudc9/X1UaqTuBuQGKYRJlf1qoRmJydIjInP14vg+aFJaID3t/6OWADnfOhi/LvN44GIN5dZvUX+yr6556ZonrW9Zdm+fkdJf0x3peAXPJvUUSRhGhVGyzKavUWokcGPy31xsHZkjCY8nNRQQf9eCwA5Utl3l51wOsA/PMblczr/7Ijni6z+f0JWwJw26dPA3BkTlCCUeb+f9dK1qFfrqJxUlFVrMncT669EGwXpPAXx0so8Lvu4emqLZ+ROT/oB0UyrTxrNwC6vC15grhXsEwuHdkg+EidytKVnX18/tJlG5XI7HGTrlsQ1RPdWtFU7Z/xdI2XvUhWgY2WZRAJcgV8m2VbaowrcnT8ea/0B2DA99735tuI9BE9F71Xsf6RraPctZuigDp5KqOwQhTUyrjGI1lD/5HVXQHofLxopkgP/ezKhvSo3KbgW13j2a/onKa30z2yfGsfxyIBWZb48pLJNQc+WB/ljaKRANzYXeOU7netiwr6ulT3ZWyCJCEiOTmVZTdL91R/RpdqPHbK9JFr/r4KaK8Zx93rW9OYfLjNs3xVKgr2gvv1W8m9MDXqpynRBnwAtakwOyQ2VytSfQA8BnybVHDg0RT3DREiRIi2h1b+AHDO7dfYNlKNArrBzN4C9vKLUi440JowNF2Zonnpyiw9zM+8o77i1KwxvdlhNwmkfTdDxcuP/1HCaUtmaXb5n0LNyh/64cgqbefGNRPKWOF193PbVd4/+Z/JKUfHAgCGPKFh33G6cgq2PkGz9bh3SD76l6PU3zmaRa/tnUOpn2rnj9O+F0dUkvnIg6smD55XoJniepdR4zgsj6vdXPN5E15Mr77IjWQxMukwLx97CwCLjtHs7ta5BwGwbqHyKrLnyqGaVuydtLk+YXH2AkjTucWKilLrQFBXIRB4m7+gzm0rMVvbRnLbVTlmcu2DeElpZUZyIHQXWSNLILdMbZbmadYe8d9dIPLngwByp6ntinNitPcx8ftPkFP3s200pwqyYJMd72fm65795XQ5iwPt/L5peWw/+kQA1pXoInS+XX1ftq3GdJCvTFbSTefY9YyqM//AIvixdCtu6TEagOJ4VQvorqJBAFzaYRbVYXmsapvxkUOxb5V1njNDuSBzKvTb2Skztcpk6RZl3yDtx7+N9tbXVula0Nh6AamgDVgAlTCzrYGtCMwowDn3WF37pfQAMLO+wDJUFL5ymXNuTv27GiJEiBCtHPVTA92k8IW79kUPgDeBw4DPEXNTK1KlgP7HhuHIBgYAU4Dh9exriBAhQrR6GG0qCug4YATwg3PuLDPrBjyRyo6pUkDbJH733ufz6tvLTY1h6TK/j+w4BoBtuksO4Bdjfw1A/kxjwgrJCvT+UaZwlhf04tcyzbNmyPTMW+BrCHwu59e6veTwihbL7K/o24W0lXLsrdpXSXrZi73UR1BG0luxX43WvrmFonkivizEmt6iUbIOWsqqYtEMBW/p/aSRSu8f8Kb6/u5BtwIbyvDVZiIHTvHmwDYZoh+28VJRB3kt+OIHRNP8ZrYotB8+1Dm386xNfmF7csb7LwEF5J28gRPWvNxAfI2ohGSHbX0QW1OVjojtJ4mF6EejqxwbF8fFvbPZpwhYieiImScFUdKaG3UZJRpn2SkjAcgo1j9I3vS1frv1fFMqCz2gfv6yZGsALu+s6zm+TNct+O/ZMVP3zDxfP6Gzp6jKXQWxmPrV40Htk7ZC9NqAB0QXLTxGgQnDT69a2yLAcZOUkD9vVhd+4cs4vrRGc7pD22mfU/KVX/T2etE4h2ZXda6viet6p23hE1GXrK6kOmPTCwG4+ikd59jfpFbsKtHxnLWnZDdmlqkEa6+o8nCyoy1QM6CNWADAeudc3MwqzCwfWAL0SWXHBtUDcM59b2a7NGTfECFChGgLsFqSZFsZRplZAVJtHg0UA1+lsmOqPoDElLgIsD1Qi8ctRIgQIdow2pAPwDkXsDH3mNnbqOh8Sip6qVoAiXnqFcgn8GLqXWwdCCIrDs1e7ZfIfN2vlyJ/PinbhXQvSLG6r0ztVf1FO2z5R21TsaRq/YV4EIPuI0BsjS/VaIabpzT29l5KYN2O/XXUmUXa2YkOSFsrkzfTL67wDM1Z52wopXfnW4cCsL6z6KNnxiiX4cLdPgDgmvlHAPBk/4+rPfeGpvg3FYIU/qcHfKgFZ39YZf3q+Hp2fELzjAGvSmqhrKPM/Iz/+aTzZMXOpoCneqIfVw1qi+brli/bfjCZk0U7xHw9hEiJ6JgeX2hMl+yge6XwmAIAnI+hz/NKslc+9yQA+2Y5nilWTshZc0T9XNbjHQDyffz/sAzRRRtoOo3BmFLlAyytUL9+kTeZ9DRRYPP31T3Y8wsdP7petGX3z1SX4NteUiPlrE+ADfdC0bvKJRj6WTH/3kEh4yd30sTxT3OPATZIMIxaK9mJPTPHABuu511XSQIib4Eif1xZOZFuun5B2c+B9yqCaNqZ+nEFNGVNWB1fT6b/rT66zSMAHP25/uPG7HNPrfs2JdpKFJCZfeCcOwDAOVeYvKw2pPoA+NE593zSQY8Hnq9h+xAhQoRo02jtTmAzy0IZoJ3NrAMbFKzzSbHAfKoPgKvY+M++umVtAoElEMRBvzFNs7Hy7WNESjWGsRw/o8/ULGvFUjm5Cp73zsOgCHqHAgCyvvEWwpb9AYhOKsTa+axOr32f/ani/c3PLvvcKyf0qkOVuZmxyh9rmGZ0D99/GADbnPgjcR//vLaP3jPmKfZ7sq8e9ki/oPpS1UuaHFdeHTaldRAUCM+2dKb+Uk7CiSfLeX7RNMW5X3CbxnZg+jIAblygcZl7gxz20RL1P+tHzdRdSUllxnGdjuLk/ACP2CpZbdGPRuPydL0sQ9fFrSwCIBKTYzKWqWvS61Mda/FOGuvup2vmOyQ9sDjzOCJHVuHTC3cG4NIZygc4p4/qUByeU31NkONz1Z84iq0f/OqllfdqcNXS1un4aT5PgTk6VsFUWR1B1u7vblb+SZ9XFcVdMXcBU3ZRW6+MkjP8Z53FIARiiEvLNAa7fvsrLV+oWXzHzjp6rheAi5eWsG6kMqCz5ignpWKB+nHRUb/RMd58zJ+Lxj45YCE/ks1/i5Rt/dx89efcEZ8BcPOK7QC4uvNEf+7NeM+2fgvgt8AlQE8gUfZhNXBnKg3U+gAws8OAw4FeZnZ7wqp8Uq8IFiJEiBBtC671U0DOuduA28zsQufcHQ1poy4LYAEwCjgKeZcDrAF+35ADhggRIkSbQCt/ACTgXjO7CNjbf/8YuNc5n15eC2qtB1C5UTWFh1saDakHkCoCmuTiBXswMFs0w2sLlPrwt0ESZrt9ngptd8wUPfHJZBWL3+J2xbdHZiv2umL58sp2o8NFUViRqKb4UrVdKWI2UOezbLdOVfpTWiBzvMNUb36PTKOsQNcpvb/a6vS0nITz99fy9n1EEXy4vYpnv1rcH4D9cyRD0RzlABuDYMzj/lcWJ14pSVFTXYLCCo31cbep1kDvF1VsPb5E4xrpKKqjYv6CDZIOjcgVAAnNBc7MymW+7Uh7OVvnPCDK44ddRG18XOJLH76qHI09dlNM/aP9Pq2kPWKeegqouf+skOTCZR1nVDnW+DJRZEF+RYBSV87WT4vKGfiStklb6ulJL10RL6wqILfNy6oa+NZjuwPQ8x45vl1ZGbat8jIq8tX3WceIltllF5WoHLNQlHL2uzrnoqG6blvcrDYrFuj+x8Vh95EArBqkPnd4tmpuxZzLJWL51bnKPQjEBYPr/Nq6DvRNkwO7yOcZ3LdwXwB+mC1qaPp+j1AbGlsPILdTH7f1YanNcb958g+btB6AmT0ApLNBn+10IOac+3Vd+9ZFAT3nnDsB+MFsY4PIObdtA/obIkSIEK0frTwPIGFivpNzbkTCqg/NbGwqbdRFAV3s349oSAdDhAgRok2ibRSE+RblZMXMbJBzbgaAmQ0EUjJ96yoIs9B/PM85d0XiOjP7F3DFxnu1PQRm+PldP2J4uiJ3TmuvB+h3Pv46wOEdFR0xpqPM/sXXannRfJnwW9wrczuelc7k36rdV/cVNRDxRtSX67RtQVTm9RHtpPgY8zOOVZ66e3ilYv1fvns/+r0samnWCYqx9pn5tCvUe9oPWrDDCtECp45UicOD2olSWO8UXZJp6ZssFyARqUQlBTryZU6US6eoqI2uPxP1Mydd2e5r++laXLWf6LpVsRweeepgAPrc5IMjAuVQn0sQ1BAI6J2aKKNk+idxm+WP61r8MFLXd/hnipCJL9Q9NGInjf1+HTyNUlbCyAxRLKtciW9M/bnQbxPQXIFWfkD9BNFS6Z5GOWrysbSfGkT9+aZ8noLztRbi5aJU3FxRQeMPV2GAnkWifiJ9faRg8ToqMjUeS3ZQ/3pvrTzPPtmKOlp8p/IAlitgji1u8/RbEBG17TAACv67mDleZmPlaI1Dh61Fl0aWKxqq/2OKPrrrRE1aP1+u30Oa/8c9q+fn/OHS8wH423+kdpufrvPvUBDIagjJEWzvrW+QuEG1aAMPgOAGuAz4yMxm+u/9gbNSaSDVf4KDqll2WIr7YmZRM/vBzN7w3weY2TdmNt3MnjWzmrWLQ4QIEWJTwKX42nTo4lUaRgL3Ah/61/3Adqk0UOsDwMx+Z2bjgS3MbFzCaxaQUqqxx8XApITv/wJucc4NBlYCZ9ejrRAhQoRodphL7bUJEQVykVJDGl7E1H9OKeqj1iggM2sPdAD+CVyZsGqNc25FSgcw64280zcAl6KC8kuB7s65CjPbDbjGOXdIbe00ZxRQgDhx1sVlkq+Iy/RfFZdJ2c9blu+skxl9WI6SuEo91bbTq4oY6PeGTyC7dAm/66fU+0NyFCGRHO0QDdQu6zDESl05h/76dwBkrFT/IuvVv7KuMrNX95MRteYQRQlN3uNxAK5cLDP7/bsV+THqmtQUGZsLjUk6C4qmLIyJJhlfJrrtrZWK2Hp7lM61xydGeY6s45KOeu84xau5vqmIlEZFB/nrds7U6QAc2676QicBbRNc90RMLheVMbFM99NryzRhe7SfEsImlovGCSjJAIt88uKez19WSQBscZfuRbzkSKCWWkl3BVFLuUrecr5QTRAtFOnRjXiBosoW79Ze2/jLkz9b47ZqgH4Aq7fU92F3FAFQeKxoMP8zocsuiyip0JfMR0RLphdrrI/5j5IV3z1IdNHsM6RWmraraKYBHfSXsuba3szfS1FIncdp35VDdQ5fXqDIoZHviOoM/r6O2U60Vtx3/M4dnmpcFFCHPm7EARfXvSHw5YuXb5IoIDP73jlXa13gulCXD2AVsAo42R+wK6o4k2tmuSkWhLkV+CMbnkidgKKEsNJ51JC2bGbnAOcA9O3VdNxeiBAhQtSF1p4IxgYfQMMbSDEP4EjgZpRyvAToB0xyztVaEMbMjgAOd86dZ2b7ImfFmcDXnv7BzPoAbznntq6trea0AIJZ6Vcl8OlazU4u7CCGKxC9ClLTg5KLwQwu5q/BNK9X3j9djrhfjT+D3w+VUNtLi/WQ7palmdltvTTLi/h9A0sgQHVx8Ev8bPKwW/8IQIZXF8hZqhnS6n6aIXUZo9nxmn6aQS3ZT7PmnPba/7OdH6gsPN6SEhDB7H1WhWal762V/MUt3ym/YodBcxjzleQ2Bj2vWW5knbaN58i6Wd9dM+myPO+w9d1OW697OH+yZpK2rhTncwMs21fI8zUYgrKTwdSxYuGilM8hmEmv/5nkCeJpavOhWzQrDUTOAod7uXc8J1oAgVXw6lo5sKeulwVwdZcf/D66nqPKtM/SCsXd758jp+yB38vRvHZSB7p/pfbnHq5zGXbBmKr9zfCuNS9hYZleciGigQti99OGDoQKHbeiq46XtkCz8RV76r7PnymLxGI65uJdNZ9bt4d+B+ljZUGkrYOI95t3mKLrF/XO/OhqjUuk1G/gLRZ8nY4V++r3XTB+FVMu1fl3fV/nkLtAbc0+RN+POlACgTkRLb+u6wQAvvNt79q/sFGz8rwOfdzIfVOzAD5/ZZNZAB1TZWJqQqrT6uuBXYH3nXPbmdl+wGkp7LcHcJSZHY4sh3zgNqAgIYa1NzC//l0PESJEiGaCc60+D6Cxf/6QehRQuXNuORAxs4hz7iOgzieec+4q51xv51x/4CTgQ+fcqcBHqIwZwBnAq/XveogQIUI0H9qAE7jRSNUCKDKzXOBT4EkzWwKsrWOf2nAF8IyZXQ/8ADzYiLYajYAC2SML9sia6pdmVdlm/3Yq/VgcFw2Q44OEZ3t99qcX7wpASUxDumO3uTw4Zw8A5syXoyw9W+bpW55eur1Qct3n9/8IgDsKpcl+SHcFTF3RSbHhcRyPFYl2SF9X2WkAVg1Qfzr/TNTUogyZ7D0/1+WxmOgeO60IgAeLtuHuHyQZsvNApfHf2U/x881ZKjKQeQjKcg4rkPsoa2fVPLjt8WPo/6Xoq+gyX/IxW+Z+1JfWzC0UZeG88mQwQ7M879ws1jnH162rdIBSXNx0J+GpunkH6v2un0l2o0dU1/yOIlEYX64UlfVQf53b/auk2Pra4hHMeUWOz3ZLvATGaaKqLu0sSuNXM48B4O4BL2k703ksqNAx18yQ3MWZh33MF3/RvTfsnao5DNFOPknE016VTl8/bmXD1U8boByXigrHul6iXCLluq9zStR2h3GeVivzOQVp6kfPV4vU5mv6Huuq/VYNySV7mbZdvJOuX9fvvVJohq97Mc9fG0+pxTvJ8Zy9RNutHZBH1y6a3Hb6Um3FZuv+HlAq1nnJXurfF76c6nW/EAUU5Iw0Cdr4n3sqSPUBcDRQggTgTgXaA3+vz4Gccx8jkSKcczOBneuzf4gQIUK0JNr67D4VpFoUPnG2/2iNG4YIESLETwEOiP30nwB1icGtoXpDyADnnMtvll61AqyMi3ZI91zLIysUR98hXc/CKzoqBvy/yxV7/rueKnF41vsS4BvwoqMi30sYnOnN2WsVhXFT71MBKO6t9dcfeDgA923zBACLYjKJvyjRsXfMLK+kg976ucr7LZio6JHuX+ryzP9S1E/BYn0v2kLUT6fviwBwf9axP1oxnA6HiN5a0qMq5fPiWpnVHxTpGDf2+BioPo49GUFEUYBUI4t+015qI7+54G7GnCMK6KmVu6mvC0SlrPtK8f65c0VZZKzRsaJlOtftrpXcwzdL+mv5/Z1o9+oooO54/5RUQz31E8hCDH1UtNJ5GWcAkD9ZP6NVW/tSkb1EmzyxWv35z0vHADD44UXkb+vj2k9TG9u0173xxlrRMsd1U7+7R3P9u7pwR5HkEobepyizL//SiXiZ2ogWFAAw/3zlEvR5pWpkU1DAJuZpsbSvRbHZUNFRhcd2ovNY9St3nM8p8NIQldFT/t18URl8hJFLUwfn76N7J2sFLN9K+8Y8i7q6r753mKrr++PlorG2+pf/a/G3SmlHjeO6rhFWL1N7nWK+DKcvrrSmlxr9YpyUdv/vkKernOv40j7+U1VV1YZgs7cAnHOtS0M4RIgQIVoKTRQF5EPdHwO6oQn1fc6528ysI/As0u4pBE5wzq1skoOm2rdU8gBaA1oiD2BSeQkZ/vP1CzQrD4p250U08+kY0TOzxMdr/1iuZ+TjS2Uh/Lf3xwDct2oQT/xLbawa7PXZ/eO2p3dszTlWx7p2NwVB7ZOtEoLrfIB7lnc0F8XT6RLV7PKz9RqDa549CQBfYZAe7/mZ22rvQC31WaAdC7TcO0hdaRnmy1j+6cNXABiRoW0/KpHzcJ8szbpeX6v8vL98qICt0UfeClCZRzCmrISOPuj72L8rP+Hta/4DNK1DOcic/rpUA3jF1ecCkLVMxy7P1Sw0f6IXQ8tII57tY9/HqpxkvLSkyfoTWASRdPUntosck/P3kaWUN1u/qUAH342Uo7KsQxad/lIIwMT3JJBW0kX30et+bIPM3+eLZQVe8Y6u87D7inQekzSzTRSpi+T4vA4vAmc7KzPaxihwwXxJUqK+DoW//qtHestq1hqswmewrwpqV/ixLJUDuTKbeIDuv5U7KO8lZ6GPz/+N9r9o5Ifc/LmS+tt11j1XPknnMnB3BR2sul9tdPhcgnLFvoRkRbb6t2KrCAOe93U15nrH/3rlTyw7XcGHBado30XvaMb/1UW3APDyWlnCZw79qnF5AO17ux13vTClbT9+98paj2VmPYAezrnvzSwPFdc6BuVErXDO3WhmVwIdkkU3mxubXhYyRIgQIVobUhWCS2H+7Jxb6Jz73n9eg3TReqHgmsCn+ih6KLQoQn2FECFChEiCAZY6O9LZzEYlfL/POXdfte2a9UdKnd8A3RIk9xchiqhFET4A2FCeb0Ush8sni+74ZIScS5lWvQM08H53jeom6dVdVNEvZx0NwKDcpSzb3t9AgTfJv+1+nbT6b8gfA8Bu3lkWSRLwm1Yuc/z4Fy6h3VzRSAUzRIf0iIkCmH2EjLiC6XKspa1TG9Hvlc8Qm6V4+7TOonfc+vVYV+UlXPRvaa6v3E5tTj/yXgCWeYfoQTna96ktZW4/vVqOt3N9DP/aeAan3SMRvKevlhxCU1I/wfkf8oFS8qPZ6mfFvupfu5lyRGYv1fZZ3fy5fzmhkiJpFoLT3y8Vu0m9JH2l6Ik+NyoWPTh2ZIgct7MOUb+u++UTdIrqnG74SFTF/P08bWQ6p12vEL2VP0ttDhuvnJDYKi+bkCAbEg3E3fz1iu8tJ/CcQ9Tm4HkqNRrzshgBSneRcz1rqeideEYaafP1P+Q6F6jvWaJlAomI9YPU1vLhGvMOh0uaYulb2q5HZ+3/xL8PJ6e37tXOr+teyFosWnLtKJ1zu5U6bqyr7tnMFWV+uSis7CV5WLEoOxc4oYerz8t30j2wU74u+oz+Ov6kclFUE9dVKyvWIFjqUUDLUqGbfC7Vi8AlzrnVFpwbiqiprupicyOkgEKECBEiGU1IAQGYWTr683/SOfeSX7zY+wcCP8GSpjuB1BA+AEKECBFiI7gNekB1veqAaar/IBLQvDlh1WtICgc2kSROGAWUhK99nHNRXJEV+3gFz1fXKu7+zj+fAEDxaTLJjxswBoBZ62Uin9L5awD2zCqpVPsMyh+mqr4ZbPfyWkVP/PPfp9H1+R+rbDPrEsXql3aWid7ta9+mv5ztvHpi2qdjquxXqRAJTPm/bbXtXJnPMb/qoKNFZx7X8TsAfvmhchtePvBOAL5eL2rj5teO5OHj/wvAHlkbK5gmnnsqCBRDj5t2pPr3TX8A2k/TOGYv07laXCcZUAf2ha9N5Jqhhp9FNmo3bbDKI1Ku/lbMnrthWyDixzjiabfiHUR9ZC0rY+lI3Vc9XlHE15qddU/nfaMImdgSURuNqVcQRAXZILVd3knf09Z4CYYJymGppKpycrBAPsJH27guomdWjtTyoL5CwQyvdJqr65q2ztNPXhm1LD9K+3HKbSjvKooqY+qCKm2Wd8j2/VFbq4eKIitrp/Hr+sVSWKo2nB/7Dm+oe6PmehmLKaKX0rfR73DUzvKlBpIj0R7TGhUFlJ/Xy+28/fkpbfvBp3+uKwpoT+AzYDxUJsz8CfkBngP6ArNRGGijBd7qg9AHECJEiBDVoYkmx865z6lZu/+AJjlIA7FZPACCGfW1SxUf/bcu44HqZ+K7ZgZDotnSiTN9LP8f5Gyyv4ume21LFQJ/09cP2CJLTrDzRinLd2i3JTw46AUAOviaAq+t0wxop0zNiHpENRNKniUHTuk9fJsrto/R5Vsd301UXPvAx9TG/KO0/Pir3gbgneFyT0e3lMMWPxuM+SpRlFdUzmiHXCBn9II/KofhqF98CcD/ChXXPslbGVstlBPx2ExVJdt1kGavW+46a6OZf4D6zPwD3LxCx111m2Z58X30A+z6ZqHvu59xL1220b71RUoZwKCx8jP7tO5ePK2zd9Z/PRbYMOM3P9ZuvRyY8RXK6clapv3sm/F0+dIfb2B/ADJW+cxaL9RWY38C529tVo7fJsgHSFuiyeSSAzWL73mvsskrHeT+WLHitaR1U1x/xSDdT2lLdb+083H+Ea/p73wVsdzvfC2oHN3DgXBfTnms0nGcMb9I2/gsXlunGX/GerXpMnQNCkbpN7VusKzoiikbsngXXSKRxWf6ijnZ55lLATjvd68A8K/vlXNQ7NTmhPImYjRcmygK32hsFg+AECFChKg34m2DHm8MwgdAiBAhQlSDeuQBtFlsFg+AgOop91rh9SmB+OxAlXXcs89vAVj9rhxbb/YV9XNu+6rFzI7Z8yEASl0FuZFc/1kmd0FEpvmCWJBbIIdb1LODORbUGkiv8v3onUcz5c+iHeLefC/rI7O+19tyGt7b9VAABmbIgbtucAEA2TMKq/QvUUIgQJ+3RVU8t5VqDmTNlHBccR/ZwO2/E3XQ/X/KU/n26H4A/H7EBxu1VV8EpTZvf+lIsjyz03Gd+jjoEjmhg+LRWOrXLYiRj9VQDyCZaol4uYTASV4Zd08C9bNIVEVkRREA8cDp66UVYsu8fEIwcwzomq/kpF51yi4sl9+dJ0+4A4DTv1OJx0GX++Lvy5dXf0L1cHCndRWdE8iBdL/lyyr9NU/jBDRY+f7bkTZJdGPE0zMVXXS/ZU7XOQeO7rStJWtRspXi7dN9mcfoEtUoLd62WyV1Ejjp0+b7c/LlOGOddG2sRFd27WBRo7nfyhHu8vOw9jr+jidp7H435zA/Dnp78B/Kt9ny7EJgQ/5Jpw0xDo1H+AAIESJEiM0QDgh9ACFChAix+cFwIQXUWJhZFiojmemP9YJz7m9m9giwDxDY2Wc658Y0Z18AhmZJWfDEmYq8uqyXImfSfRp+O6tgSHputft+frtkEmqK5Q9kC4L9E6Ng0j2Vs1OmtnlitZQgt8iXXMOPZYoSGpaxvkqbQR7BN0v60yFXkSVp3mxnhtQ/V+0mNcSMLTWUkQLlDuR85PMGAiXIoERiNXHt8R8VFz7sX9KHj02SimWlAqSnFJaO1Dn3eUQU1XkPz0s5tyFAQIe9sU4U1k2jDgLA9Sqj/yuiE9yo8dXvXA8axHKCEoeiGepSA437+PcgDj6aJwoitmYNFQu9vn4QZRO0tftIACq+qj0PYfa1imTJWA2Drh0DwF//osirfmXat6LaPWuGpaVXXp/kc6vwuQSV2wYRT56aCuivtC6SBImOnQ1eNTYy3+8bqMd6SmzF2epv/hzdR1lTNCYVvRS5s34L0WSZK8tJn+UTWrNEJQYRTiVDRCGmF+k6R4o11pFS0TfxrronIiVlFB6v9m7q/hwAH64V9TQ+c0sA7vi7KLSdgroFzYH4T98EaG4LoBTY3zlX7FOhPzezt/y6y51zLzTz8UOECBGi/ggpoMbDKc048MKl+1eL2VXBbDPIDszy2vU9s4sAGJyu74G+fSoo93UA0r3jdllMjt2aLAfYENe/2M+8Spz6E/dDkR/RjOhfS/cE4IVxEvV6dZ+7AMhOL2P9liosnrlUxxt0n2btHdIUj31sgZy/5x56CQAFT6rIeNTP3oM6AdH8PFx/OfAiPi67Et5pmIxgBjzkbvV75q/6VK5LdeYfWAoHXHwBAMWny2Lp+5Rmpxlvf9dkN4ZFo5Wx+A1FZd5EAirj/beQpVSernOPBk7VzCRdfj/z7vfXLyvbaKr/lGjXzjhfTD06Rw7cRMd1IoIZf9D/VcfJ2Z+2zheAX7ge+0GOfsvVbHzG31Tprt08nds63X50GuuttAIflODHILgvI4tWQL7/LazyY+jzALJmyRm8bpivJbBGFkDWUr1PuUjH7jCqgPHn6d6/eomqw/2ps+7nhad9BjTzzN9jc6CAml0LyMyiZjYGCR2955z7xq+6wczGmdktZpbZ3P0IESJEiHqhibSAWjOa/QHgnIs550YCvYGdzWxr4CpgGLAT0BGotgqOmZ1jZqPMbNTS5Q3XRgkRIkSI+qHpxOBaM1osCsg5V2RmHwGHOuf+4xeXmtnDwGU17HMfcB9IDK6+xwyonwAn5MoEfWeFdNyfWK1Y/gu9vn1D2uzqi3evdzKNs70xU+4qKh3B6z0V1dmX5PtFnhx/OSaTPN3kAnzlHTkLh7wqB9yJM5X2joOsC5TWv2699vln508BOHPsmQD8uEY2etYKn7Lvzf6Avgliw2efPYS+ryjgfvlucrR1+lwUwpoRctJlzyys0kYlPKXxt5OfqWF0asapsw4EIPdlxfa3/1riehXzJtW7rbrgYrFqKZyGwtI8ZeedrZHpilcvel45DF1n672icLbfIVLZj7oQ0DLxsurpt5pQsWAh5nX+q8vtAEjzchNuufI85p4rqY1yz9D0ec8fc/QklvtSi1G/KHO5F3+brra7fav7e10f0TnR9aKPMlZ5qYhlCfRTseigSirI696X9BFllf2Fr1UxVFTi9Iv0O9ny+iIA5hzZhS9KNIZ754mayvVyKtd1nVDtuTY5HG3+zz0VNKsFYGZdzKzAf84GDgImJ2hgGyqD1kJXNUSIECFSg8VcSq+2jOa2AHoAj5pZFD1snnPOvWFmH5pZF6SQNwY4t5n7ESJEiBD1w2ZgATR3FNA4VP8yefn+TXmcQHs+oFryI1XLOCbHqj/c97Mq6x9aLerjV/mLG9yH7CQ/dmIewNK4qIAuEVEo7f17gH4+RX7f/aUuOfda9bfP114v/ZRdWOvj5rtOV1tnf6RSjGk+rGT9k16Tfj/faJJ6ZBAb3ufmNcQ8ldFpjaKBAtXKHJ/uX9NtP+NMKUWelLuyhi02RhCJtfoYq9Kvinnza9ql1SGZYgmifDoeoUgsBkoaI1pQAECsqCjltutL/VTbLz+m0a1UNrG0pxRhF/vyjdnLfAiPv1f6veXVQmeJHizdZ1vK83R9Onwn+rHgR7W9YL8CALqO0m+s3WwF9VmZvsfyfD3TmG88IwNKdd8WD1OOQO6PoqrSV2l5xTaKolqwp36n8WK1NeNU0ZQ3nfwwe2SpvZhTX1/xtTGOaleV2tsg81L/+hO1whGKwYUIESLE5om27+BNBT+JB8DKuGa0D68aCcAVHadXWV9XrHqQCRwnXi+huFQxKE3OsHtWacZ9eDs5tt5cL+fhcbmyPK7o/i4AF/U7U/2ZLO3//Ce/or13RAazvuRS9dbXF9weJYd2RQ1ZqYlZo6nOwtO20MwypwEG0n7jTgIgb+n0OrZsOwicwoGoWmz2PKBmZ2xz9yPqxepWbKcZ96qB6lfUp3ksPlD96vOyz+xeK6sj3lOZwGv6ZNBpvDaecbycvANf1n3SeZwXdPOz98jaqvkVaeXe0V3hc5nT0ykZrnsxp1Cz9YrOuv8jJepH+gpZGRU+Sz2jo9qMrVEeQNTi7DP+WAAWLCkA4O87q1piuSsCNg7GCGb+xfHG5X9UQfgACBEiRIjNFOEDIESIECE2QzgHjajL3Fbwk3gAtI/IAXtG+zF+Sc2yDNVhWIbi4K9fti1/7dx0cemB8zmQj7j96aMA+OIwUSqj5isOOmuEzNuh6RLRWrKnnGEdPQUU33s7MmZ7sXwvcRA4dQO5gVW7yOxu/5UcuUFh8MBh2RisHi5q4bY/3J3yPoHzN++6+l2LtoDKkoqbsA+RnBycF1kLqLyOH3pph8WiGitydG90G6WeWkzbx/LkHC7PE43ScdxqyjqKjsmb5WmtLC/RscxLPCz1jv8gtt+XfQxi/Eu38qKEy9aSuUQUT/AHOv0kOaXzZori6TzBy0n4WIjyderPHT9/BIBXlu/A/Bn6DVx3gOTCTs1THkycqgEUyVjn6iurVwtCCyBEiBAhNkOEUUAhQoQIsRkjtADaBoKIgBKXWgRAQM0Ux2WK7pQp83dkxnjiPgqoKaKBgjb2HnMCAAPuEaWz/BXFZffsKurqj0eeAkDWEm3f/61CAFYdJ2mIrf44nnkH6xyTY8wDuYF2z38NQMVuUnHE0wJlh+4EQM53M7W+ppKD1SC+t1I43r9DyozJkRe1Ybd/XARAt1FSKd0MlHUbLOvQELitBhGds6jqQk/PlHTSzzp3tlQ21/b20Tarda9k/TADgPQeih6adHkB/bwwe948bRPL1L0Ya6dzKhkhijFrkY/gyZcaaPoS5QWkr9Jvb+2AfHLHLACgbIDa7+7FUPOnq2xkSVf159iffQHA2R21wTvFkmb54OMR/OawDwE4Oc9Tnyn+LgNpliZB+AAIESJEiM0QoRO4dcHhqgisBYgT5weflZjvHVJ1IZhFTClXWzv5JN50S6uMIw7Ep+rCn5aoyve5Hb/w/RH6p+VVblP8pRxaBUtkAeAduJm+6lT2ztsA0GmibrhpF/QHoMsYtTb/pC7EigpT6g9fKZu47DA/85/mZ/wFcsSxYmXKlbVefEpO3wg1z/yTRfAGvX8WAEMfHgO0zGx4UyEoJB9UEwvOta6C9IlI6ylrsGLBwvodu7ScKX8apC+ddQ2G3CInb7t5uoetPO6/V7UE2HWojpmj38HAJypIW6d9Vw2WozZnsb5bhWbB2fOqZuCmr1Sbpb11XwVZvtGSOPHOytqtyNXvq+ALn2Gep7ZPuVmz+39+cgQAz+VtD0CHAlkXP556Z+XvfLW30uvK7q9vZbqUEFoAIUKECLGZInwAhAgRIsTmCBdGAbUmOKCCGFFv4s2qUHzytPKODEn3WvlxXyQ7RXNw64yqVMh6V5oy9RPgL12+85/kLJtevjG9krZTkfqT5CQMNOtz5wVx2npP98xB7nMqg1eTrENtyHhL/WpIVPSiS1QAvBw55/KTaLeVcY39OhdjZrlorMJyyQr0eslr5zdB/kFrR2Uh+SSkRP0MEX1T0q8DAJnr1Fbg5A+oIedLeQbx9pP/o6Lol+31FmOK+wLw4/+JQpxygY+RN91Hgx7R/b90pO7ptT21vNu32m7lMK1vP8OIp/vcgHY6TqTC0zdT5Wgu7ycas6Sztls6Qm0EchN93tb1zly2nshyOXtzfIH5FQf6c+2ktv/9+tHaqYtopqO3VH2Mm7p/D0Cpcwm/Ye1TnhTfn0wFB3knyaKMDYYD14DfXVtDm3kAhAgRIkSLIrQAQoQIEWIzRBgF1LoQwaqYdz+WKcb491+fQL/uinT5YKvXgQ3mYKbXSf/Z1MMA+N/Qt6q0mWwuNsR8TN5nm4yNt9m+u9Qil7RXxIRbqfJ5gaRAp7ellDnpLzKVt7xNZnfMq026Fr4PS1V6gE6RdtWuP+CffwDgpSv+zbk/nApA+XRRQUPGKZolmXqKZIqGSFQj/akilXN1mfrpZf0wC4CKVaJNFv5B9Fvx1qIJ88Z7Tf8DRafs1WEKAM9edTjufEmHZJyjMR8S9dILExSzXyIFD4r7aSYbTxelseAIbbdFP+03vX1fYu28jMRcUTsdJ+n4c08QzRSogq7rrPWTzlGE2F+WqLzqV9/uDEDmmMJK+s8GSh4iUlk1S/fzByf9HwB9EyLlEpGYb1IXJRtQQ/XJUUkZoRM4RIgQITZPuHjoAwgRIkSIzRBhQZhGw8yygE+BTH+sF5xzfzOzAcAzQCdgNHC6cy6ljKHJ5UoW+f3rKiPsoo4393vJr5UZmGwOJlM/AZqijNyTa8SXBGqF1SE/TVTA2ONlLnf7qgiAom2UMNP+sa8AyFg1RDus8VEkQVlHmpYDShs80LfvlR9n+SIy+yip7e0z/+23rN5Ev+73DwMwtyKPSEQ/kiH3eurHt7XwMlEZkb2kItnrYtECcV92MtVktLaIgPpJThRL69SpUoojPlnSHJahe3X2tZL9GHv27UDCPXyI3oL7/vJZKpSScdECFr4vimX4z6YC8O0klVq0fFGL8w7V/XPh7h8AcOfHBwFw3W4vA5AV0XY3rjuUjjdIPXZNP1Euxb1EbfZ5uhCAin6iXF2aKKntr/8dAPlzfanIGaoWFF+3DrYSlblwD68CeoTujSv6fwLUTP00BKX+N7wqrnPZ5/nL/Zo/NK7hzUQMrunLX1VFKbC/c24EMBI41Mx2Bf4F3OKcGwysBM5u5n6ECBEiRMpwSGcrlVdbRnMXhXdAEBSd7l8O2B84xS9/FLgGqFVsPo6j1JVz+oQzASiYrNnrqL/dA7VIFdSGpiggnTzzH/LJGQAM67GYSd/3B6Dv25qd9JgqgbaiXRXj3ektSUNU+Jl+v6sVd9+EiubVomK6Zp+BZMFvJsmxePmnI4GqMhbV4YIPTwfg5J2/oXSqd2yv0LmtPmUXALY7diIA869W7YNz338NgJtmaRaaeVBho8+jIaiUaVjr8xQaYIkETl43QhZbdLZmv863GeQBJOcJJArxVdYU8H8gg56UQ3efHU8G4OuRUmebUaG2hqWr38X/p1n/8mFplO6g9pf/RUXpd75ODuVnB2rG/8paOfEv/UA/taEXK85+j5mzATjwBc2We34WJ//f/p44wktAFMg6XTdSDuXZx6ufQ+6RJTLvQLWd/U3VcVx81nakrdfMee9fSgjwrM6fAfDMSt0bQzI0XjtkVI2YCCzyqEVSlnTYc/SZAPS6UPkSW2RqHGeltHctcK5JrVQzOxS4DYgCDzjnbmyyxhuB5rYAMLOomY0BlgDvATOAIucqMzvmAb2aux8hQoQIUR+4uEvpVRfMLArcBRwGbAWcbGZbNXP3U0KzPwCcczHn3EigN7AzMCzVfc3sHDMbZWajli1v26ZWiBAh2hhcPLVX3dgZmO6cm+l9nc8ARzdr31OEuRb0dJvZX4H1wBVAd+dchZntBlzjnDukjn2XAmuBZbVtt4nQmbBf9UHYr/oh7Ff90Blo55zr0tAGzOxt304qyAISkz7uc87dl9DWccChzrlf+++nA7s45y5oaP+aCs0dBdQFKHfOFZlZNnAQcgB/BByHnoRnAK/W1ZZzrouZjXLO7dicfW4Iwn7VD2G/6oewX/WD71f/xrThnDu0ibrTqtHceQA9gEc9BxYBnnPOvWFmPwLPmNn1wA/Ag83cjxAhQoTYVJgP9En43tsv2+Ro7iigccB21SyfiXixECFChPip4ztgiM9/mg+cxIYoyE2KtpYJfF/dm2wShP2qH8J+1Q9hv+qHVtUv7+u8AHgHhYE+5JybuIm7BbSwEzhEiBAhQrQeNHsYaIgQIUKEaJ0IHwAhQoQIsZmi1TwAzKyPmX1kZj+a2UQzu9gvv87MxpnZGDN718x6+uVmZreb2XS/fvuW7FfC+j+YmTOzzq2hX2Z2jZnN9+M1xswOT9jnKt+vKWZWa95FU/fLr7vQzCb75f9OWL7J+mVmzyaMVaHPWm8N/RppZl/7fo0ys5398k19f40ws6/MbLyZvW5m+Qn7tMR4ZZnZt2Y21vfrWr98gJl944//rJll+OWZ/vt0v75/c/SrzcI51ypeKGR0e/85D5iK0qbzE7a5CLjHfz4ceAtVmdgV+KYl++W/90GOndlA59bQL6SrdFk1228FjEXKrAOQJEe0Bfu1H/A+kOnXdW0N/Ura5ibgr62hX8C7wGEJ99THreT++g7Yxy//FXBdC4+XAbn+czrwjR+H54CT/PJ7gN/5z+ex4T/jJODZ5hivtvpqNRaAc26hc+57/3kNMAno5ZxbnbBZOyQmB0qlfswJXwMFZtajpfrlV98C/DGhT62lX9XhaOAZ51ypc24WMJ1mCMWtpV+/A250zpX6dUtaSb8AzayBE4CnW0m/HBDMrtsDCxL6tSnvr6FI4h2k7XVsQr9aYrycc64mgckX/PJHgWMS+vWo//wCcIC/1iFoRRRQIryZth16umNmN5jZXOBU4K9+s17A3ITdml1ULrFfZnY0MN85NzZps03aL7/oAk8PPGRmHVpJv4YCe3kz/BMz26mV9CvAXsBi59y0VtKvS4D/8/f9f4CrWkm/JrJBx+Z4NiQ4tVi/rH4Ck5X98utXoTokIWiFDwAzywVeBC4JZv/OuT875/oATwKbRD8jsV9IsflPbHgYbTJUM153A4NQ/YWFiNZoDf1KAzoic/1y4LlNMROr7v7yOJkNs/8WRzX9+h3we3/f/55NlC1fTb9+BZxnZqMRNZRSIaemhGuEwGSIqmhVDwAzS0c325POuZeq2eRJNpicLZZeXU2/BiGec6yZFfpjf29m3Tdxv3DOLfY/kDhwPxvM8E3aLzQre8mb8N8CcSS2tan7hZmlAb8Ank3YfFP36wwg+Pw8reQ6OucmO+cOds7tgB6YM1q6XwGcc0VIV2w3RIUFia2Jx67sl1/fHlhOCKAVPQD8bPBBYJJz7uaE5UMSNjsamOw/vwb80kdF7Aqscs4tbIl+OefGO+e6Ouf6O4lOzUMOs0Wbsl9+eSIf/HNggv/8GnCSj4oYAAwBvm2pfgGvIEcwZjYUyEBKkpu6XwAHApOdc/MSlm3qfi0A9vGf9wcCampT319d/XsEuBo5XIN+tcR4dTGzAv85EJicxAaBSagqMPma/45f/6FzLsx+DdDcXuZUX8CeyJkzDhjjX4ejGcgEv/x15BgGRQPchWYg44EdW7JfSdsUsiEKaJP2C3jcH3ccuvl7JOzzZ9+vKfgIkxbsVwbwhL+W36NSoZu8X37dI8C51eyzKcdrT1Qveyzi3ndoJffXxSgiaCpwI15NoAXHa1skIDnO30tB1NZA9MCZjiymINosy3+f7tcPbI5+tdVXKAURIkSIEJspWg0FFCJEiBAhWhbhAyBEiBAhNlOED4AQIUKE2EwRPgBChAgRYjNF+AAIESJEiM0U4QMgRIgQITZThA+AnwjMrLjurerd5lFmdqX/fIyZbdWANj42sx3ruf0UMzuqmnX9zWxCdfv9FGFmf0r4nG2Shi4zLz0eIkRjET4AQtQI59xrzrkb/ddjkORvS+BU59xrzXkAM4s2Z/tNhMoHgHNuvZP+zYKaNw8Ron4IHwA/MXiJgP8zswmmoh0n+uX7+tn1C6aiLE8GYmxmdrhfNtpUbOQNv/xMM7vTzHYHjkLqlGPMbFDizN7MOntNpGCm+oyZTTKzl4HshL4dbCom8r2ZPe+Fxuo6nx1MxT/GAucnLI/68/zOpHz6W788Ymb/9efznpm9aWbH+XWFZvYvM/seOL6m/vhjfuLH451AXsPMLjIVSBlnZs/U0ud2JiXWb83sB5NybGDBfOaP970fV8ysh5l96sd2gpntZWY3AsGs/8mULn6IEPXFpk5FDl9N8wKK/fuxSCI3CnQD5qDiHvsiKdze6MH/FUr3z0JyuQP8/k8Db/jPZwJ3+s+PAMclHO9jvAwBEnUr9J8vBR7yn7dFyqk7+m0+Bdr5dVfg0/iTzqOyXf99HLC3//x/wAT/+Rzgav85ExiFBPqOA97059gdWBn0G0l2/DGhzxv1B+nLfwl08ctPTDifBWyQGCio5Vr8Azgt2A7JJrQDcoAsv3wIMMp//gPwZ/85CuQlXtOktgvxsiPhK3w19hWo54X46WBP4GnnXAxYbGafADsBq4FvnRc8M+mp9weKgZlORTxAD4BzGnH8vYHbAZxz48xsnF++K6KQvvCGRwZ6CNUIL/pV4JwLCpA8DhzmPx8MbBvM7pHK4xB0/s87qaEuMrOPkpoNFD9r6s8WwNbAe355FMlqgx5GT5rZK0jcriYcDBxlZpf571lAX/QAudPMRgIxVCMBVGXrIZP65ivOuTG1tB0iRJMhfABsXihN+Byjcde/gg0UYlYK2xvwnnPu5EYcM7m9C51z71RZmFADuQasra0/ZrYNMNE5t1s1+/4MPeCOBP5sZtu4DUVIkvt2rHNuSlLb1wCLgRFo7EoAnHOfmtnevv1HzOxm59xjdZxHiBCNRugD+OnhM+BEz5F3QX9YtcnyTgEG2oZi2SfWsN0aVAAkQCGwg/98XMLyT4FTAMxsa0QDAXwN7GFmg/26diZZ6BrhpPdeZGZ7+kWnJqx+B/idnzVjZkPNrB3wBXCs9wV0Q9RXdaipP1OALma2m1+ebmbDTfLHfZxzHyG6qD1Qkw/jHeDCBB/Ldn55e2Cht05OR9YFZtYPVSO7H3gACAq9lwfnFyJEcyB8APz08DKiKsYCHyLOe1FNGzvn1qPC2W+bqjytQb6CZDwDXO6dmoNQmcLfmdkPiE8PcDeQa2aTgL8jSWOcc0uRT+FpTwt9RWqVnM4C7vKUVWIFsQeAH1EhngnAvciieRHVZ/gRyU9/X9351NQf51wZeqD9yzuexwC7oz/rJ8xsPJIjvt0/oKrDdciXMM7MJvrvAP8FzvDtDmODNbIvKi70A3oA3+aX3+fbCJ3AIZoFoRx0CMws1zlX7GesdwHTnHO3bKK+fAxc5pwb1Yg2gvPphKyfPWp7CLYl+GirHZ1zyzZ1X0K0fYQWQAiA3/gZ9kREU9y7CfuyAvHgGyWC1QNv+PP5DLjup/Dnbz4RDFkW8U3cnRA/EYQWQIgQDYSZnYUqZCXiC+fc+dVtHyJEa0P4AAgRIkSIzRQhBRQiRIgQmynCB0CIECFCbKYIHwAhQoQIsZkifACECBEixGaK/wd5AWKROnXL7wAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ds.tcc.plot()" ] }, { "cell_type": "markdown", "id": "2ffc20a5-9b96-4db5-a69a-15a4de78d1c7", "metadata": {}, "source": [ "### To go further, search using a custom MARS request:\n", "\n", "Here we use a set of custom parameters corresponding to `TIGGE_CF_SFC`, which should result to the same request sent to ecmwf.\n", "\n", "If you want to eodag from a MARS request obtained from https://apps.ecmwf.int/datasets/data/tigge, you can put in this `ecmwf_req_params` the content of your MARS request without `area` and `date` settings, which must be passed to eodag as `geom`, `start` and `end` parameters." ] }, { "cell_type": "code", "execution_count": null, "id": "c060d3b3-3121-46ea-89ae-2a7b5f99ebe2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Request using productType or directly MARS parameters result to the \n", " same product title TIGGE_CF_SFC_20211101_28b0b1a0da723aba1682661621bd7b21fe12f1ea \n", " and downloadLink https://apps.ecmwf.int/datasets/data/tigge?date=2021-11-01/to/2021-11-01&area=70.0/-130.0/30.0/-60.0&class=ti&dataset=tigge&expver=prod&grid=0.5/0.5&levtype=sfc&origin=ecmf¶m=tcc&step=0&time=00:00&type=cf\n" ] } ], "source": [ "ecmwf_req_params = {\n", " \"class\": \"ti\",\n", " \"dataset\": \"tigge\",\n", " \"expver\": \"prod\",\n", " \"type\": \"cf\",\n", " \"levtype\": \"sfc\",\n", " \"origin\": \"ecmf\",\n", " \"grid\": \"0.5/0.5\",\n", " \"step\": 0,\n", " \"time\": \"00:00\",\n", " \"param\": \"tcc\",\n", "}\n", "\n", "products_from_mars_req, total_count = dag.search(\n", " geom={\n", " \"latmax\": 70,\n", " \"lonmax\": -130,\n", " \"latmin\": 30,\n", " \"lonmin\": -60,\n", " },\n", " start=\"2021-11-01\",\n", " end=\"2021-11-02\",\n", " **ecmwf_req_params\n", ")\n", "# title and downloadLink properties must be the same with the two request methods,\n", "# as they are built from ECMWF request arguments\n", "if (\n", " products_from_mars_req[0].properties[\"title\"]\n", " == products_from_product_type[0].properties[\"title\"]\n", " and products_from_mars_req[0].properties[\"downloadLink\"]\n", " == products_from_product_type[0].properties[\"downloadLink\"]\n", "):\n", " print(\n", " \"Request using productType or directly MARS parameters result to the \\n\",\n", " \"same product title %s \\n and downloadLink %s\"\n", " % (\n", " products_from_mars_req[0].properties[\"title\"],\n", " products_from_mars_req[0].properties[\"downloadLink\"],\n", " )\n", " )" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }