{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Calculate Geospatial Weighted Averages from Monthly Time Series\n", "===============================================================\n", "\n", "Authors:\n", "\n", "* [Tom Vo](https://github.com/tomvothecoder/)\n", "* [Stephen Po-Chedley](https://github.com/pochedls/)\n", "\n", "Date: 05/27/22\n", "\n", "Related APIs:\n", "\n", "* [xarray.Dataset.spatial.average()](../generated/xarray.Dataset.spatial.average.rst)\n", "\n", "The data used in this example can be found through the [Earth System Grid Federation (ESGF) search portal](https://aims2.llnl.gov/metagrid/search)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Overview\n", "\n", "A common data reduction in geophysical sciences is to produce spatial averages. Spatial averaging functionality in xcdat allows users to quickly produce area-weighted spatial averages for selected regions (or full dataset domains).\n", "\n", "In the example below, we demonstrate the opening of a (remote) dataset and spatial averaging over the global, tropical, and Niño 3.4 domains." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2018-11-28T20:51:35.958210Z", "start_time": "2018-11-28T20:51:35.936966Z" } }, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import matplotlib.pyplot as plt\n", "import xcdat\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1. Open the `Dataset`\n", "\n", "We are using xarray's OPeNDAP support to read a netCDF4 dataset file directly from its source. The data is not loaded over the network until we perform operations on it (e.g., temperature unit adjustment).\n", "\n", "More information on the xarray's OPeNDAP support can be found [here](https://docs.xarray.dev/en/stable/user-guide/io.html#opendap)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "ExecuteTime": { "end_time": "2018-11-28T20:51:36.072316Z", "start_time": "2018-11-28T20:51:36.016594Z" } }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 1980, bnds: 2, lat: 145, lon: 192)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n",
       "  * lat        (lat) float64 -90.0 -88.75 -87.5 -86.25 ... 86.25 87.5 88.75 90.0\n",
       "  * lon        (lon) float64 0.0 1.875 3.75 5.625 ... 352.5 354.4 356.2 358.1\n",
       "    height     float64 2.0\n",
       "Dimensions without coordinates: bnds\n",
       "Data variables:\n",
       "    time_bnds  (time, bnds) datetime64[ns] 1850-01-01 1850-02-01 ... 2015-01-01\n",
       "    lat_bnds   (lat, bnds) float64 -90.0 -89.38 -89.38 ... 89.38 89.38 90.0\n",
       "    lon_bnds   (lon, bnds) float64 -0.9375 0.9375 0.9375 ... 357.2 357.2 359.1\n",
       "    tas        (time, lat, lon) float32 -27.19 -27.19 -27.19 ... -25.29 -25.29\n",
       "Attributes: (12/49)\n",
       "    Conventions:                     CF-1.7 CMIP-6.2\n",
       "    activity_id:                     CMIP\n",
       "    branch_method:                   standard\n",
       "    branch_time_in_child:            0.0\n",
       "    branch_time_in_parent:           87658.0\n",
       "    creation_date:                   2020-06-05T04:06:11Z\n",
       "    ...                              ...\n",
       "    version:                         v20200605\n",
       "    license:                         CMIP6 model data produced by CSIRO is li...\n",
       "    cmor_version:                    3.4.0\n",
       "    _NCProperties:                   version=2,netcdf=4.6.2,hdf5=1.10.5\n",
       "    tracking_id:                     hdl:21.14100/af78ae5e-f3a6-4e99-8cfe-5f2...\n",
       "    DODS_EXTRA.Unlimited_Dimension:  time
" ], "text/plain": [ "\n", "Dimensions: (time: 1980, bnds: 2, lat: 145, lon: 192)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n", " * lat (lat) float64 -90.0 -88.75 -87.5 -86.25 ... 86.25 87.5 88.75 90.0\n", " * lon (lon) float64 0.0 1.875 3.75 5.625 ... 352.5 354.4 356.2 358.1\n", " height float64 2.0\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " time_bnds (time, bnds) datetime64[ns] ...\n", " lat_bnds (lat, bnds) float64 ...\n", " lon_bnds (lon, bnds) float64 ...\n", " tas (time, lat, lon) float32 -27.19 -27.19 -27.19 ... -25.29 -25.29\n", "Attributes: (12/49)\n", " Conventions: CF-1.7 CMIP-6.2\n", " activity_id: CMIP\n", " branch_method: standard\n", " branch_time_in_child: 0.0\n", " branch_time_in_parent: 87658.0\n", " creation_date: 2020-06-05T04:06:11Z\n", " ... ...\n", " version: v20200605\n", " license: CMIP6 model data produced by CSIRO is li...\n", " cmor_version: 3.4.0\n", " _NCProperties: version=2,netcdf=4.6.2,hdf5=1.10.5\n", " tracking_id: hdl:21.14100/af78ae5e-f3a6-4e99-8cfe-5f2...\n", " DODS_EXTRA.Unlimited_Dimension: time" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "filepath = \"http://esgf.nci.org.au/thredds/dodsC/master/CMIP6/CMIP/CSIRO/ACCESS-ESM1-5/historical/r10i1p1f1/Amon/tas/gn/v20200605/tas_Amon_ACCESS-ESM1-5_historical_r10i1p1f1_gn_185001-201412.nc\"\n", "ds = xcdat.open_dataset(filepath)\n", "\n", "# Unit adjust (-273.15, K to C)\n", "ds[\"tas\"] = ds.tas - 273.15\n", "\n", "ds" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2. Global" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ds_global_avg = ds.spatial.average(\"tas\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'tas' (time: 1980)>\n",
       "array([12.52127071, 13.09115223, 13.60703132, ..., 15.5767848 ,\n",
       "       14.65664621, 13.84951678])\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n",
       "    height   float64 2.0
" ], "text/plain": [ "\n", "array([12.52127071, 13.09115223, 13.60703132, ..., 15.5767848 ,\n", " 14.65664621, 13.84951678])\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n", " height float64 2.0" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_global_avg.tas" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEiCAYAAAD+lrlwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABbAUlEQVR4nO29eZhj113n/Tlaq0q1L7253W63t3i3E9s4wR6yOBsQCFtmAoEwk2F9mZl3huUF8rIMkMBL4BnCMAkECMkDSRiGbRggkJ0sju20E+/x7m67eqtNUkmlXTrvH/ceSa4uVUn3nnOuqut8nqef7lKpSr+W7v3d3/3+NiGlxOFwOBx7h1jUBjgcDofDLs7xOxwOxx7DOX6Hw+HYYzjH73A4HHsM5/gdDodjj+Ecv8PhcOwxnON3RIYQ4oQQ4q6AP/uoEOKVpl9nWBFCSCHEhhDiXQF//i4hRFEI0brQ3hvHzjjH79iVSCmvlVJ+LuzvEUK8UgixqMGkrX73bwkhnhJCFIQQjwshfmCH53+vEOKk79D/Vggxu8NL3CilfGcQ26SUn5JSjgPPB/l5x+7GOX6HwxwbwJuAKeDtwHuFEK/Y6olCiGuBPwC+H9gPlID3WbLTscdwjt8RNTcJIR4SQuSFEP9TCDGiviGE+FYhxANCiJwQ4m4hxA1d32vLN0KIUSHEh4UQWSHE14UQP7NFFH/e6wghMsDHgUO+7FEUQhzS9R+TUv6SlPJxKWVLSnkv8AXg5T2e/n3A/5FSfl5KWQR+AfhOIcREP68lhDjqyz//Vgjxgv9e/KgQ4lb//50TQvyenv+ZY7fjHL8jat4CvAG4FLgB+EEAIcRLgQ8CPwLM4UXDfyeESG/xO34JOAocA14LvK2f15FSbgBvBE5LKcf9P6c3/6AQ4md9x7nln37+k0KIUeBW4NEeT7kWeFB9IaV8BqgBV/bz+7v4BuAK4F8DvwO8E7jL//1vEUJ804C/z3EB4hy/I2p+V0p5Wkq5Bvwf4Cb/8R8C/kBKea+Usiml/DBQBW7f4ne8BXi3lDIrpVwEfneA19kRKeVvSCmne/3p89f8Pp5j/+ce3x8H8pseywN9Rfxd/KqUsiKl/ASe1PQxKeWSlPIU3h3HzQP+PscFiHP8jqg52/XvEp4DBLgE+MlNkfXFwFZSzCHgha6vX9jiOb1exzhCiPcA1wFvkb2nIhaByU2PTQKFAV/uXNe/y1t8be3/7RhenON3DCsvAO/aFF2PSSk/tsVzzwCHu76+eIDX2XE8rRDi57tyAOf92eFn/yuenPQ6KeX6Nk99FLix6+eOAWngyf7+Gw5H/zjH7xhW/hD4USHENwiPjBDiW3okO/8C+DkhxIwQ4iLgJwZ4nXPAnBBiqtcTpJTv7soBnPen188JIX4O+F7gtVLK1R3s+AjwJiHEnX7S+VeAv5ZSDhrxOxw74hy/YyiRUh7H0/l/D8gCT+MnfrfgV4BF4DngU8Bf4uUD+nmdx4GPAc/6kpK2qh7g3cAR4KmuO4SfV9/0v77Tt+NR4EfxLgBLeNr+j2u0xeFoI9wiFseFhhDix4B/I6W8YCtYhBAVvIvb70opfyHAz78G+Cs8OembpZSf1WyiY4hxjt+x6xFCHMQr5fwyXinjPwC/J6X8nSjtcjiGlUTUBjgcGkjh1flfCuSAP8d1vTocPXERv8PhcOwxXHLX4XA49hjO8TscDsceY1do/PPz8/Lo0aNRm+FwOBy7ivvvv39FSrmw+fFd4fiPHj3K8ePHozbD4XA4dhVCiJNbPe6kHofD4dhjOMfvcDgcewzn+B0Oh2OP4Ry/w+Fw7DGc43c4HI49hnP8DofDscdwjt/hcDj2GM7xOxwOR5/85F88yB994dmozQiNc/wOh8PRB4VKnb/52iJ/9IXnaLV293BL5/gd2rjn2VW+431folxrRm2Kw6Gdrz2foyXh7HqFBxZzUZsTCuf4Hdr4yL3P87Xnczy0y08Kh2Mrjp9YIyYgGRd8/OEzUZsTCuf4HVqoNpp89vElAB5azEdsjcOhn+Mns1xzaJI7Lp/n44+cZTfvMnGO36GFe55do1htIAS7/jbYNPlyndO5ctRmOAag3mzxtedz3HLJLG+87iCL2TKPnl6P2qzAOMfv0MInHj3LWCrOq6/a56SeHfjlv3uUt/zBl3d1xLjXeOz0OuV6k1uPzvLaa/YTjwn+cRfLPc7x90GzJfnCU8vuRO1BqyX55GPn+KYrF7j10lleWCuztlGL2qyhRErJF59eYTFb5pSL+ncNXzmxBsAtR2eYyaR4+bE5/mkXyz3O8ffBXxx/ge//4/v4+CNnozZlKHlwMcdSocrrrt3PDYenAFzU34MTqyWWC1UA7j+ZjdgaR78cP5HlyOwY+ydHAHjj9Qd4dmWDJ88VI7YsGM7x98FH7vV2GXzo7hPRGjKkfOKxc8RjgldftZ/rL5pCCJfg7cV9z60CEBPwVef4dwVSSo6fXOOWozPtx153zQGEYNfKPc7x78DDi3keObXOVfsnuO+5Nb5+ZjgSOrVGi4eHxLl+8rFz3H5slqmxJBMjSS5bGHcRfw/ufW6NuUyK24/NcXyIHP9KscpSoRK1GUPJidUSK8Uatx6dbT+2MJHmpUdm+OLTKxFaFhzn+Hfgo/edZDQZ54/efgsjyRgfHpKo/8/uOcm3/Y8vcmJlI1I7nlku8vRSkdddc6D92A2Hp3jghfyu1T9Nct9za9x26Sy3XDLD18+ss1FtRG0SAP/+w8f5oQ8P13rTQqVOcQjeH6Xv39oV8QMcmR3j3PruvFg6x78NhUqd//3Aad5040Eunh3jzTddxN8+cIpcKfrE5eeeXEZK+ELEEYeq3b/rmv3tx248PM1KscqZ/O48KUxxKldmMVvm1qOzvPSSGVoSHnwhF7VZvLBW4oEXcjy4mB+aMtNKvcl3vf9u/t2ffCVqUzh+Yo2ZMe9OtpuFiTQrxequDHCc49+Gv3vwNKVak7fedgSAt7/iKJV6i//5lRcitatSb7a14rsjdvyL2TIT6QQXTY+2H3MJ3q35ynNe5HjbpbPcfMSLHochwdutU3/66+citKTDez/9FE+eK3LfiTWeWY42gXr8RJaXXTKLEOJFj8+Pp6jUW0NxVzIozvH3QErJR+99nqsPTnLTxdMAXH1wktsuneVP7zlJM8IhTV89maVSb3FoaoS7n1mN1JbVjRpz46kXPXb1wUkSMcGDQ5CD+MqJNX7t7x8biqjs3ufWmEgnuPrgJFOjSa7cP879zw+H47/+oikunc/wiceid/wPLeb4wOef5a6r9xMT8Ff3L0ZmS63R4rnVDa45NHne9xYm0gCsFKNXAAbFOf4ePHwqz6On1/nebzjyoiv9D77iKIvZMl+KMNL+wtMrxGOCn3j1FeTLdR49HZ2DXS1WmRtPv+ixkWScqw9ODkXE/6dfPskfffE57vOj7Sj5ygmvMiQe846nl10yw1dPZiOd9LiYLfHgYp5vvv4gr71mP/c8u0qhUo/Mnlqjxc/85UPMj6f47bfcyDdducDffO1UZMHN2XwFKeHwzOh535v3j3tVnrubcI6/B6rU7vXX7n/R4yqz/1yESdUvPb3CzRdP8zrftigrC1aKVeY3RfzgyT0PLeYjH1+rpJQ/vedkpHasFKs8vVTktkvn2o+99MgM65VGpFLGxx/2elO+5fqD3HX1fupNyeefjO54et/nnubxswXe/R3XMzWa5Ltedpgz+QpffmY1EnsWsyUADk+f7/g7Eb9z/BcMZ9YrpOIxFjZFs3OZFPGYiKz0LbtR4+FTee64Yp758TQvOTAR6d3HarF2XsQPXoK3UGnw3Gp0F8gzea87djaT4p8fPRtpueLxEx19X/GyS6LX+f/h4TNcd9EkR+bGeNklM8yMJfnkY9E1Kv7N107xTVcu8JqrvaDmrqv3MzmS4C/vjyavtugnuy9yEf/e4EyuwoGpkfMSOrGYYH48xdJ6NB/23c+sIiXcecU8AHdcPs9XTmSp1O3PwG+2JGulGvOZ8yP+YwsZwEv+RsVXT+YA+MVvvYZ6U/IXESbl731ujZFkjOsvmmo/dul8htlMKjLHfypX5oEXcrzxuoMAXhPeS/bzmceXqDdbkdi0XKhy+b5O9cxIMs6bbjzEPz16NhIJ6lS2jBBwcOp8xz8zliImXMT/IoQQHxRCLAkhHtn0+H8QQjwhhHhUCPGbpl4/LGfznuPfin0TIyxH9GF/8ellxtMJbjg8DcA3XjFPrdHi+An7zmNto4aUMD9xfsSv7gJWIzwp7j+ZZSQZ41tuOMgdl8/z0XufpxGRQ3v09DrXHpoileicckIIXnpkJrIEr5op/y3XH2w/9tpr9rFeaURyPG1UG5RqzbaEoviulx2mUm+1ZSmbnMqV2TeRftHnpojHBHPjaRfxb+JDwBu6HxBCvAr4duAGKeW1wG8ZfP1QnFkvc6in409HFvF/8ekVbj82RzLufXS3HZ0lGReR6PyrG957MJc53/Er3T/KaOj+k2vceHiaZDzG226/hNP5Cp/x+w5sczpX5uIt5IIbDk/x7PJGJFvLPvHYOa4+OMnR+Uz7sTuvWCCViPHJCKp71LEyv0k6vPniaY7NZ/irr9qv7jmVLb+oVHkz8+NpF/F3I6X8PLC5lOLHgN+QUlb950RzFu5AqyU5l69yYIvbO4B9k2mWIrjKn1zd4IW1clvmAcikE9x8ZCYSnX+l4JWxbZXcHU8nSCVirEZU6lauNXn09HpbR7/r6n0cmBzhz+593rotzZbkbL7CoS0ciHIqZyPoAH1uZYMbuqQn8I6nV1w2x+eesH9qqsh58/EkhODOK+Z5LIJxKadyZS6aGev5/YWJNMuunHNHrgTuFELcK4T4FyHErZZfvy9WN2rUmi0O9oj4FyZGWN2oWpcNVFR/R5fjB0/nf+R0nqzlUcjtiH+L5K4QgvlMKrIa54cWczRasj1YKxGP8a9vvZjPP7lsPUJbKlRotOSWCUJ1jJ3J282FNFuS1WKVfZPnf3ZX7Z9gMVu23vugPpfNUg/AvskRCpWG1VxWsyU5ky9vWcqpmB9PseKknh1JADPA7cBPA38hNmdPfYQQPyyEOC6EOL68vGzTRs76owZ6a/xppPQuEDZ56lyRTCrOsa5bc/BKTKWEh07ZredXTn2riB887T+q22Clm998cWe+ynV+dGt7LIF6va0i/oP+Y2dydiP+1WKVlvSO5c0sTKSpNVusl+12pKrIeXMlnbIJ7FbQLBUq1JtyW6nHi/h339gG245/Efhr6XEf0ALmt3qilPIDUspbpJS3LCwsWDVSRV+Hekk9/kFoW+f3bjtHz6s0utS/ELywVrJqz0qxSiImmBpNbvn9uUyqfVdgm/tPZLlsIcNMV8VRFM4D4JTv1LdyIAf8+e62pR4lVS5MnB/ctN+nol2blgtVhIDZLarE2uecxZLcU9nepZyKhfE0tUaL9cruGttg2/H/LfBqACHElUAKGLq5puok7BXxL0RwEIJXGnl4C71RVR3Ydvxe127qvAuRYm483c4D2ERKyf3PZ9v6vqLjPOw6/u0i/tFUnOmxpHWpRx27W0k9CxEFNivFKrNjKRLx893SPv8CZdMmtSFtq+YtxW5t4jJZzvkx4MvAVUKIRSHEO4APAsf8Es8/B94uh/Ae6Uy+QjIumNsi8gBPbwT7DuRUtrRl1BiLCS6eGeXkqm3HXzuvAqOb+fE0qxv2b4OfXdkgV6pzyyWzL3pc2Wr9Ti1bZmo0yXg6seX3D06NWpd61HuwldSjnKztkuXlQnVLfR86Fyib59xiHxH/bm3i2vpI1ICU8q09vvU2U6+pizO5MvsnR4jFto5kFyJwIOuVOuuVRs9E05HZMZ6PQOrZKrGrmB9PUW9K1iuNnnKQCe73a9BfuiniTyVizIwlrUsYp3PlLaN9xcGpEesjrDtSz3Do6aDGf2x9PM2O2e+YP5UrMzOWZCzV200qe13EfwFwJl/pWdEDHQcyTHrjkdkxXlgrWY2uV4pbd+0q5iKq5X9gMcfkSOK8JDh40WwkuZnp3seT5/jtSz3TY0nSifh535sc8Upxbd/RbhfxR9Exfypb3jbah+gukmFxjn8Lzq5XtmzR7mbfxEgkt51bafwAR+YyFKoNciU7be1SSlY3qlt27Srm2927dnX+F9ZKXLowvuUdWxQ9GP1E/NlS3Wqp4tJ6dUuZB7xS3AXLHalSyp4D/xS2O+YXe0ir3UyPJonHhIv4dztSyh0jfrDvQE75UwJ7HYhHZr0Lgi25Z6PWpFJv9cyDQKej1/bYBi8JvvX7ZNuhFXyJbjsHohoFz1qUe5YKVfZP9j7G903afZ/U8dQr4ge7HfNSSk7lti6m6EbdiURRxBAG5/g3sbZRo9Zo9azoUSxMpFm2WIK3mC2TTsR6RkS2Hb9y5ttq/BP2pZ5WS3IqW+5ZibHgOzRbkthpP2m7XcSvRoOctij3bCergHeBtClldrp2t3H8FoOttY0alXprx4gfPJujmt0VFOf4N6GSbDtG/P5tpy0H4kUf59fwKy6e9Q5QW46/M1eld8Q/O6Ycv71oaKVYpdZsbRvx15ot8mU7kth2pZwKFWTYivillCwVKu3qna1YmLAb8fea09ONzY75U9uMY97MQoSNikFxjn8TZ9uOfyeNP029Kcla0tQXs9vPDBlLJZgfT/O8pZLOTtdu7xM1EfeS4DabuF7YIReiSnFtObW2A9lW4/e7dy05/mypTr0pe2r84AU22VKdWsPOWJLlbaqMOjbZ65hvF1P0G/G75O7u5sx6nxF/u67Yzsm6mC1tOzME4JI5eyWdq304fvV9m/qn2pjUK1Kz3cR1KlcmERPbOjTbTVzbNW8pbDcm9Rfx2yujbjdvDRDxD2FLUk+c49/E2bx3om6nXUNXk4sFB7JRbZAt1XeMPmzW8qsTdav2+m7mxu2ObdgpwrbddX06V+bg9Eh7z24vDkyOWJN6Os1b20s9YO/OaLlQJdZjXIPC5tiGxWyZTCreV//J/Lh3929LPtSBc/ybOJOrsH9y5xPV5ryefqOPi2fHOJMvW7k9Xy1W2/Xe2zE3nraq8S9mvVWLmR5dsvssO7TTuXLPmU/dHJoebSeCTaPudraXeuy+TyvFKrOZ9LbnnU2ZbjG79VysrdiNYxuc499EP6Wc0B05WnD82f4c/5HZMVrSzvTJlY3atjX8igXLiyq2K+UEb0/AaDJurSzwdK7Sl058YGrE2qC2QaQeW5LYTlVG0NUxb+Ocy22/gKUbVeAQxY6OoDjHv4mz671XLnaTSSfIpOKWbjs9+WanmmJV0nnSgtyzUqgyv8Xmrc3MZVIUKg2qDTvNSb3mGSmE8PR2Gydpo9ni7PrWC1g2c2hqxC8hNP8+La1XGU8n+hpFYE3qKda2rRADux3zp7KlHc83xUJ7bMPuqeV3jr8Lr3mr3FfED96tpw0Hspgrk4rHtpxT3o3NWv7VjVq7Tn875ix270opd4z4wZMxbDi0c4UqzR4LWDZjs4lrudC7a1dhe67RSh8RP3h3Iqbv1sq1JuuVRl8BoLIJdtfYBuf4u8iX61TqrZ4rFzfjNXHZ0RsPTfceGqfYN5EmbWk882qxuuWu3c2oKM6G418p1qg2WjtGal4jkHmH1k8Nv6Kzicu8XUuFytA4WfAu2MvF6o6BDdgZlZIre8fq9Fh/gwWnRpMk47trbINz/F20uyz7jfgn7DiQfoZFgT+eeXbMeC1/vdkiW6q3h7Bth4r4VyxU9vRTMw/2xjZ0asF3Pp6U4z+7bj4/s1SothOl26G2S5mmUG1Qa2w/rkFh425NbR7rd6KsEMIvW3aOf1eiTrp+b/FsDWpbzJY5PN2f3mijpFPt9t2phh+69E8r75OfC5ndQeqZHGHdwv7WUwNF/N5zTFf2SCm3HdDWzb6JESsXyH7GNShsjNxQZZmDjBLfbWMbnOPv4kyfXbuKfZNpSrUmxaq5tWuVepOVYrWvRhKwM555uY9xDQp1V2Cj23Kxz27LBUuJy9N9zHNXqCYu0xp/sdqgXG/25fhVEtx0Y9LKAI5/38SI8ZEbQRz/bhvb4Bx/F6rDtB8JA7pr+c2drIPMDAGvlr/gN3yZQun1OzW5AYyl4owkY1YmdKpNVxMj25+wC5a2Oe00jnkzBybNL2Rp1/BvU8qpsLVPVgUS/eYdwOxnF8Txz2V214RO5/i7yJfrTKQTJLfY+bkV7T2gBg/Cnebwb+YSC5U9qhO3nwitrX9aSO72M9YCupuTzDrZs+vV9jL1fjg0PWp8bEM/XbsKdXEwfWfUifh3DrhsNE4GcvzjadY2artmbINz/F3kyjWm+szkg509oDtt3tqMijBNNnENemc0Z6mJq59STrDXnLRerg/kPA5MmR/b0G7e6rP5rvtnTLFSrBGPCWbGBnD8Bm1Sjn+nO8du5jIpas0WBYOyr06c4+8iX6r3XcIFdrTixWyJeEywv48TFTolaOsGNdCVjSqpRIyJHmMRNjOfSRmP+NXijIv6SILPZdLEhPlItlCpMzmA4z84OcKq4SauQSJ+W/Xpy4Uqc5nUjuXK0BnbYPKive7f+e80tqUbFQSt7ZImLuf4u8gNGKFNjSaJCciVzH3Yq8Uac5kUiT7lJ2W/yeTXSsHbtdvPHBPwJCHTGn+2VKdUa/YV8cdjnvxkUi5otSSFaoPJkf4ujtC5gzSZCF8qVEgnYkyO7myXrUGEK8X+mrcAv+M4btSm9fJgF2zoDJezOZAwDM7xd5Er1Zge7U++AK9ufnosxZrBE3V9wKhxLBUnHhOsVwwmd3fYtbuZuXHvPWq1zOmfnbEW/TffmZQLirUGUjLQZzftSx1Zo46/yr7JdF8X7cnRBKl4zHzEX6z2lS9SmB65kR8wAITo9ksHxTn+LvLl+kAaP8DMWJKswYh/0INQCMHUaNJsxF+sbrtrdzNz42kaLbNja08NmATfZ7g5ab2tE/cf8U9buFvzavj7SziruUY2kruDOH5v965ZjX9Qx9+J+J3j31VIKcmV6u2Tr19mM2Yj/iAH4eRIot19aILVYm2gE7U9tsHgbfDigElw0+MI1Ps/OUCCcMZ3HiYDCW/l4mDRtckLpJTSm/vUZ6EAmG8sC+P4TfoCnTjH77NRa9JoyYGSuwAzYymyG2abSQY9CE1G/FJKL+8wkOM3P71wMVtiYiTR93u1b2KElaI3RM0EBV9qG0jq8Z9rsgdjZcCLtumIv1RrUm20dlzos9mmYZN6RpJxxtOJXdPE5Ry/j3KUg2j84F3pTUZo6+X6QAlC8JyNKY1/vdKg1mwNFKGpigeTJ8WZfKWvhSeKfZNpWtJchKaangaJ+JXGnzd0PDVbkvVKvX1n0Q+mnax6/wexaX48RbFqbtR3EMkXOrms3YBz/D6qMmfQD3x6zHP8Jho3VGXIwFKPwYi/n92om5lVSUuDkWyuPFyluErj76d6RpFKxMik4sbep3y5jpReXqpf5v3GpEbTzFY3FTTN9lHDr1AXyJyB96lS9+5ABj3nwAsCXXJ3l5EvqYh/UG0vSb0pjczrKVQGrwwBL8o0pfH3u2S9GxXNmaxxHrQET13g1Qhe7fYoqWeAiB88p2bCoUFXdD2Ak50dM5twDhLxm9TTO5/bYHfZ4PWHuOTuLiOn2rQDaPyAEZ0/SOu4ev56uW7kLkRF/P127QIk4zEmRhJGJbFCpTGQk1XvqakLpPq94wM6kOmxpLG+EPV7B3GynYSzWcc/SJXYjMGy186dWgCpJ5OyMpNKB87x++RKwTV+MFOJsR4gQeg9P0Gt2aJqYOn6agCpB7yTwmi/Q7k+kKyiLhKmciGFSp2xVLzvuU+KaYPlwcp5DyL1tJ2sIZvCRPwmLkZBgy3oaPy7YV6Pc/w+g27dUbRlDAMnRpiIv/vndbJcrCHEYM4DvPfJlONvtrtk+7dpsh3xm0qC1weWecCXegzZpJz3QFJPxlx0rWyKx8RA0oo69obpnAPvvWq0pNFSal04x++TL9VJJ2KMJOMD/ZzJ286gB2E7mjXgQFaLVWbH+h8hoZg12OFcrAy2MQlgIp1ACIOOv9wY6A5EMT2aNKbxq2N0kOBGPddkxD8z1v/4D8+m4TvnoKt7dxeMbXCO3yc34IA2hapGMOHUhjHiXylWB9L3FSYb3YJIYrGYYCKdMDZrPmjEPzOWIlcyM94iW6qTjAvG+xyuB92JVHMa/yD6PtAeEGjknCuFi/hhd3TvOsfvky/XB9b3wWvJj8eEGY0/YKKpLWMY0K8H7dpVzGZSrBkqe1UXuED9DgalnkHGNSimx5K0JEbG++ZKNaYHjK5Hk3FSiZixhHN2o85MJoAkljGTBM+rjuuAGj/sjnk9zvH7DDqLXxGLCX9ej5lEUzwmyKQGk5/MR/yDO/6ZTIpao0Wppr/pJnASfMRcv0Oh0gjkPDpNXPrtypZqA+dmhBBGZbq1Um2grl3F7FiKNQPv0XrApDx45ZzgpJ5dRa40eJu2whvbYEbqmRpNDhShQSfyNZFk8iL+YFIPGKq9DjAXB/yyV1MdzuWgUo85TT1bqrcvLIMwk0kZLecM4vhnMmbPuSDMWuhX0YUxxy+E+KAQYkkI8UjXY78shDglhHjA//PNpl5/UDypJ7jjN9NMMnjXLnQiX93RbKXepFBtBJN6DOZCOhH/oFKPmWF2UkrWKwGTuyYd/0ZtoA5ZhakJtM2WJFcKZtPsmJlRKWEcfyrh9avsdY3/Q8Abtnj8v0kpb/L//KPB1x+IoMldgJmMmRMjH2BOD3gNU2OpuHb9Wh3QgSL+cXNlr4FzISNmIv5SrUmzJQOXc4IZmS5bCqanm4qu18t1WnKwGn7FtMG77CASnWJ+fHd07xpz/FLKzwNrpn6/Tir1JuV6M9BtMKhBbWY0/qAHoYkJnZ2l2CEifgO3wevlOkLAeGqwi6SpKabqYjLIzlaFqfJgb+x4LZjUYyjiV0FAII0/k2Sj1tS+pnLQPcmbmd0l3btRaPw/IYR4yJeCZiJ4/fNYD1G7Cx2NX3fFSpiD0EQ0q5JWQZO7YKrDucFEOtHXztZuJkeTlGpN6poHkBXUZM4AUo+6w9PdxFWsNmi05MDJXfAu2rlyXfsIayX7BdX4Qf+gtjBSD5jvUNeFbcf/fuAy4CbgDPDbvZ4ohPhhIcRxIcTx5eVlo0apkyyo1KM69nSX4IVx/GYi/uBSz+RIgkRMGLkNDrIjVdkEHUet0x7v9w9uUyIeY3Ikod2hqVlSg3TtKmYyKaTU3+wWZGicYtbQKInQjn88ZXTvhC6sOn4p5TkpZVNK2QL+ELhtm+d+QEp5i5TyloWFBaN2BZ3FrzDRSSilDCX1mEhcrmwEl3qEEOa04oDNUqaS4EHLSxXTBhKXQcY1KNTP6M7PZENE/CbOuXrTKzcOK/VkDTXg6cSq4xdCHOz68juAR3o91ya5EN164OmNoLdipeRvBAss9RiK+MfTiYHHWihM1YMHHY/QmdCp2fG3y0sHtwk8TV17xN+ezBksuQv68w6rIRz/rIEZWZ07tWCfG3i1/E3D+6V1EPx/uANCiI8BrwTmhRCLwC8BrxRC3ARI4ATwI6ZefxBUB2Dgqh4DiyHCzAwBcxp/kHENClNjG9YrdS6Z62/JejemOpzDJHcBpvyxDTppT58NmNwF/dMwsxs1xlLxQIGEuoDpvBi1z7mAfgC6unc3aoGqlWxhzPFLKd+6xcN/bOr1whD2AzfRnKScRxiNv1Bp0GxJ4gMmPXuxUqwGknkUs5kUXz+7rsWWboI2S3WG2emVxFTOIMjIBvAc7cnVDZ0mhdLTTVUarZVqgeyB7nHRwxNsQVf3brHK5fvGtdhlAte5ixcNxf2hXUGYNpBoUi37QRwadKLZosbE5Wpx8IFa3cwa0/iDjUdQ8pB2jb8cbNKrwkQneK7kjdMOM3xMe94hYNcu+Mt9NA9q0+H4TXao68Q5fvw5PQFGIygm/UFtw3QQmpjXs1KsMj8RPOKfyegvC2w0WxQHnMWvmDIo9YRpAvJGSTS07rnN+iNJgtz9jaXipOIx7cndtZByyExGbxJch+NXFW8rzvEPP7lS8HEN4FesjOlt4gqv8fvzejQ5tWZLsrZRYz5MxD+WREq06tdq13GQ5O5oMk4iJowkd8MkCJWmrnNkdDaErOJVZCWNSD1h7iB1L/cJs3ZRYWO/tA6c48ev3Q2R0AGvskfnibEeYLlIN7oj/mypRksSKuKf9fMDJqK0IBG/EMIbzTxkEb8J6TDMSBJAe2AD/kjmgBcj8AIJnQUVYc858CSoqdHk0E/odI6f8E0b4A9qM+DQBl3WrdC9WrC9ZD0TwvGP6Z9Xvh5ifjqoRjfNDVyVRuCKHuhUl+l0amsBB7QpdOcdqo0mxWqjXQod1Cbd8upIMkY6ESw3o5jLpIZ+Xo9z/ISXekD/ibFerreXvARBd8SvnHWQrl1FuwRPZ+21apYKeoEcSWiXegoBh+spOuXBOiP+YHN6FLOa9fR2J/EwafylYNVhm5kbH/55Pc7xE/6kAP0HYdhhUbpr1NsRf4hyzs6iiuHRZYdT6tFfN58t1QPN6VFMa142pCL1sFViJY2D2nTc+YPZNaO62POOv9mSgefedzOb8U4MXYPawh6EmVSceExoi/jV/JGFEI7fRNNN2PEIkyN61y9KKf3kbniNX1fEr6bPhomuZzN6dwGHGSGh0N04qcvxz42nh3794p53/OqkD5P4Au8gVBcRHeQDNiUphBC+jKHHnpVilWRcBKqeUaQTccbTCa2Lu8OOR5jUrPFXGy1qzVao92kinSAm9Dk0HU52eixFS+q7gwwzrkGh7mB0RdfaHP8umNez5x1/2MmcilnN80x0HIQ6J3SuFqvMZdKBex0UM5kkaxorHtYrdWICMgPO4ldMjia0Sj0qkAiT3I3FBNNjKXJlPcdSZzJnmMYkvU5WnSdhNX7QlzNar+hx/DP+RXKY5/U4x6/m9ASczKnQPcFQx0GoU79e26hrmT2ie0m2Gsk86Cx+xeRIklqjpU0nVnd8YZK7ANOj+jT1ziyq8LKKLpvWNrxO4jBFFbo7isNu31J0z+sJw+lcmV/+u0d56lwhtE2b2fOOX8dgJtA/wVBHb4HOiD9froWKGBW6xzasV8Lp6brLXsPmHBTTY0ltGr9y1kEmcyp0z+vJlrxu+UQ8uAvSaVOzJSloyPVB90UynF3Pr5X40N0nWCrorxByjl9DmzZ0L4YI70CqjSaVeit01KgzcZkN2QCkMNFtGUZP1z22IcwSlm5mxlLaNP72isOQ5ZygL7peDdlXAB15VkfOqKDpgg2d9ypsgrczUTW8TZvZ846//eaGrePXWLGi62KkM3GZK9WZCimHgf6Z/EGXsCjUxVXX+9Tp/gx30Z7S2JWa2wgv9XRKTPVp/GGlw2Q8xsRIQotNus450HeR1CHR9cI5/pBLWBTjaW+1oA6NP2w3qkJX4tLbBqZJ6hlPUa43Kdf01V5rkXo0R/xhkrugIn59Uk8mFSeVCH66j6cTJONCW0XWWojJnN3oaixTfkCXnAnhE+FZjTZtxjn+srdVKozWCJ25LwUNDkRX9DE1qidx6S0kl1puOWd1J8EDbt9S6N7C1V60HtLxT48m2ag1qTXCT+jU0aCoBhHquxiFl3pA39gGnRH/SDLOWCoe2q5cqUYqEWM04Hjv7XCOX5N2DWirm9cxJdCzR49Ty2qqfAL9Za/hpR79yd1kXDCSDHdqTWf0NXFlS3qia11OVkpJdqPOdIhkc8empJ6IX1NZt0JH927O77YOW0K9Fc7xhxhXuxld5ZM6I/7u3xeUthym8TZYx9gGtRw7zAVyoj2+WpPG70tPofsdNI5tWNOWmNeTdyhWG9SarVDjGjo2pdp9CmHI+xcPHXks0OP4s6WalmBrK5zjL+uM+PVU0YRdu9i2R5N+rS4cOi6Q7XnlGpq4Chpq5keScdKJmMaIP9g2sM2oE15HxK8ruNE1gTbM/t/NzI7p1fh1SD2gL+I3UdEDzvH71Sp63ly1OSksYdcudtsD+iJ+LVGjxvkquiQxnf0O6yEncyrao5k12JXd0JOYn9HUg6EcohaNX9Ogtny57m0aC5EA70ZH9VqY5Tk74Ry/VqlHz4jffLnOaDL8QdjewhUy76BGB4Qtee22SYej7YxkHp4O50KlHrqiB7pn8odzHo1mi/VKQ1t0nSvXQ8+gac8O0iD16NpdkCuHH83ejZaIv1wP1XS3HXva8bdakrxmqUeXQ9NxF6J+R1jnoU4qHRJGIh5jPJ3Q8z7pKnvVOMzOk3rCR/y6RiSoO1Adx/j0WLLd4RqGztA4jXeQIeca5Up1pjRG1zOZcGXLUkpypZq2nMNmdnT8QojvEUJM+P/+f4UQfy2EeKkRayxTqDRoSX0NEpOjSaoayid1TQlUvyOs88iVaowm44xoKivTJa3oyoXolHpy/iiCsKgF52Ej2bzGahVdXc4qGauj0mh6VE/E7+2/CH/BVqjEddCcyIZfQm2ihh/6i/h/QUpZEELcAbwe+DDwfiPWWEanhAEdGSNsRKTL8SfiMSZHwkfXqqxMF5OjmpLgbY0/5GgLTVKPWkgfZmeBQgjhd+9q6v7UEDnqyhllSzViIrxEB51Ks9DvU1lvBU3Y2V259l1RdBq/Cl+/BXi/lPJ/A2assUxW8ywMfVU0eipDwLubCd06XtZ7Gzw1qkfqCbNovRtd1VhrG+EX0nejo0ZdJYd1lOLqdPzTY6nAE1W70VUsoLuCZi5k2bLJOT3Qn+M/JYT4A+AtwD8KIdJ9/tzQo3sWhq5moLBrF7uZ0bAyL1eqaU186ZR64jHBWCqcBOWNtmiE3p6m1lPOa4j4wTsuQ0s9mmZRQSewCe34N/Q5WV3VTzmN5xyEj/izmn3TZvpx4G8B/hl4g5QyB8wCP23EGsvovqoqySFsSacuqQeU8xiuemJtjr/cYHIkEbpZamrUS1puhJwftFzQ7PhHwzdM6QxudEb8Oko5Aa/6LR4LdWdUqXujMXTcFSnCRvwm5/RAH45fSlmSUv41kBdCHAGSwONGrLFMR//UV9UD4SL+RrNFsapnLjjokwuG0vGHXGqu0HWnpiL+BW1Sjx6ZDsIvhgF9jn9tI/zsIIUQgumxZPvOJgidCb36ouvJkSTxmAgc8eejjviFEN8mhHgKeA74F//vjxuxxjJt/VOX49eg8esa66sIKxdIKcmX6loPwKnRJJV6i2ojXIStxiOERVdupiP16HmvpjNJcuV6KAkqX64zoWEIIXiVRomY0BPxa6xPnw4Z3LSLPDQGN7GYYGYsGTrij1Lj/1XgduBJKeWlwF3Al4xYY5lcqc7EiJ6TAronPQaXenRtBFPMjKUoVBo0msGmPJZqTWrNlnaNH8JHjvmQS1gUnYg/nES3UqyRTnh9CjqYHk1Ra7QohygPzpfCb3JTCCG87vQQn5uUkmyprrVaZXo0XHCje1yDYmYseKdztlTzR2GbSaf281vrUspVICaEiEkpPwvcZMQay+js2gVIJ2Kk4rFQkaPO8bAQPvmle2oh6Ft3uLpRYy4TXlbRtWRkuVBlfjz8QnqFjkFtwybTlfxR0zq6dhXTY+Fs0n3OKcJ075qc0wP9Of6cEGIc+DzwESHEewE93S4Ro2udoMKbyR+uVNGY4w/o1NTP6ewg1BXxrxZr7cXWYdA1KnqlWNVWygkdfTdMcj6necLjZEjHry6uupK7EF7qyRuSVWYzwYfa6Q5KN9OP438QKAH/Gfgn4BkulOSu5hIuCF8Trtvxh239N3FS6HD8lXqTYrWhpYJG16jo5UJVS/OWQsccGq8HQ+9nF+b4Vl27Oo+nsHkspfEPU8SvOyjdTD+O/1VSypaUsiGl/LCU8neBW41ZZJG8gavqRMgJne2mJN2OP+gtp8aRzAodjl9nIlXXxqSVYo2FCY3dn+2LdrhoVmd+RlvEr1nqCTMqJV/2+kF05WYUsxmvlLoZYKidjq1p29HT8QshfkwI8TDwEiHEQ11/ngMeMmaRRUxcVb2BXyGqeoZM4+80khiI+ENEaatFzy4dGj+ES8SBGtdQ1VbD79kULuL3diXr1vjDSZkmGpOUlBX0ApnzL466N13NZlK0ZLAAJ1fWOyZlM9td4j6KV7b568DPdj1ekFKuGbPIEs2WZL2iNxoCLyI6lSsH/vl8uc5IMkY6oWcgWniNX3/iq9MBGvzOqB3xa9LU58ZToaSe9rgGjY4/7ByajVqTRktq1fjVzgkpZSBHqS6uOiP+7gvkwanRgX/ehOQLL166Psj/t9meGmwu4u/p+KWUeSAPvNXYq0dIoVJHapzMqfA0/hBSj8bFMIC3SD4mgmv8/oVI12ROgGQ8RiYVDxU5diL+4ViVp7t5CyCd8CSooJ9dJzGv925NdTkHkUbWSnWE0GxTyDsjnSWv3XQ7/kFYL3u+yWTEb2zmjhDig0KIJSHEI1t876eEEFIIMW/q9XfCVIOEN/clXHJX50nhdTYGH9uguypEEbYscGVD73iE2bFU+2ISyB7Nc3oUMyESlzp3JSvC5mfU2Oq4hgFtirBrKvOal7AoVI5mUMdvQl7djMlhax8C3rD5QSHExcBrgecNvvaOmBp7OjWapBYy0aS/kSQZeCG1qeqCsEnClUKNTCrOaMgBbYqwEX9nTo/+4ymMQwN9I0mUPRA8P7O2ob+gQm2pCt6rYiaRqkqNBz2uOr0z0ZZzBkJK+Xlgq1zAfwN+Bgg3CjEkJqIh6OoCDRj1m3H8wWe+5A05/rBlgasbVeZ0JlJDbkzSnXNQzGSC16jrXGquCDuh09sjq/d40pHcNaHxB63KMj2LHyyPVxZCfBtwSkr5YB/P/WEhxHEhxPHl5WXttuhewqII25XqjSHQa9NUiM5G3QsqFGGlHl3NW4qwG5NWijVSiRgTmksCp0dToSJZGK4ejOxGXWtiF2AkGSOViAW6C1GrJE04/pFknEwqPrCEqO7Od6XGvxkhxBjwTuAX+3m+lPIDUspbpJS3LCwsaLdHRUO6r6qdZeLBErw6Z/ErwkzoNNU6HlrjL2ounQzZvbviN2/pLgmcHgs+mtnEKIKwk0yzBurThfAGogV5n9bLZnJ9ipnM4HfbbY3f0L5dsBvxXwZcCjwohDgBHAa+KoQ4YNGGNiq5qzu6DjPpsdFsUdA4klnhST2DT3mUUmrv/FSEd/w1rXp62Pnpy5rHNShm/MR8K0ATUL6kvyJLHQthpB7dET94TjJIcKN7Qu9m5jKDlwnny3ViAiY0jNLuhTXHL6V8WEq5T0p5VEp5FFgEXiqlPGvLhm7ypRqTIwmt1QUQLiLqjGTWrIGOBZvyWPYXVJjQGqdGk+3fPygmmqXCzuvxxjUYcGhjSVoSCtXB7yBNaNfjqQQxEczxl2tNKvWWmTvIsWQgSSxnuIJmJjN4Y6DO1ZS9MFnO+THgy8BVQohFIcQ7TL1WELKlutYJgYowW7hMTQkMOuWxs6DCzIkKAbsaS16zlK4afgg/r8e7A9Ef8YcZ1GYiPxOLicDL6U0MaFPMBFzG0jnnzMgqQarFsprHbGyFsXsJKeW2jV9+1B8ZOUO1u2EiflOOv9t5XDTdf2ejyYXP3UnCQZuelHPWWdWjNiat+f0Bg6DuQHQ2bym6u1IvmRvsZ3OGGpOCynTKAZooU/SkntzAP5c3rPHPjg3u+E1V0nVzQSxND0K+VGPKwAE4koyTSgSbyW/O8QfrbOxMLTQj9UCwiH9F825b6GxMWgvQ75At6R/XoAizK8BUY1JQx6+OPyMa/1iwbWWmlrAogpQJZw2PZIY97Pi9LUBmPuygYxvMST1B64nNR/xB7oxW/AhKd7OUd1s+eMSve8l6N527tWCBxDBVZK2169P126TyWJX6YDkjU+ecop07GuDcy2ledboVe9bxe6MIzHzYUwHHNgyrxm8quQvBIv5Vv1lKp9QD3kkapMNZ967dbmbCaPyGHEjQrmuV5DSRWwt6Z5Qr1Y2uOJwJYJeX3HVSj3YazRbrlYaxq+pkwK5U9TO6S0zbUePArePmKh5CST3FKvGY0H7hns2kWA0Q8ZsY0KZQfSGDXrQr9SbletNIJBu067pTn24g4h8NLmeaivZh8Du2aqNJqdY02rwFe9Txq4obU1fVoFu48uU66YTeumuAVMKbhjmo88iXzNgD4Vr/V4teLbjucreg83pWCr70ZMDxJ+IxJkcSA0f8uvc6dDM54kX8g+rp2Q2vhDphILoOWv2kexruZgaVWTsb75zUox3TtbuTAbdwmTwIp8dS7Qi+X0zecoYZzbxSrGot5VTMjnnjEQbdmLRcrBoZ16CYyQw+tiFnsFplajRJvSkH1tOzJf3jGhRBFw6ZyoMoBpVZswbl1W72pOPPGr6qBt3CZWJAmyJI67+3mcjcARg0SeitONQfXc9mUkg5eNRoalyDYno0GaIHY3jyMybGNSiCV66ZdfyDyqw2RjLDHnX8eUMD2hSqwWXQW2GTjj/IhE5T4xoUQZOEqxuGIn4/WTzo+7RcrBpJ7CqC7FMwWZ8exvGbivjDVK6ZlHoGlVnbCXAX8evHxLjabiZHgt0KD1vEn92oGXGwisARf6GmvaIHOh2lg05TNHUHoggygMzE9i1FYMe/YS66HknGSSdiA9nk7SSuGevaVQxy4V7yS4P3TZo7nmCPOv6OjmYq4ldjGwbXG4cp4s+WakZK7xRBqkNKtQbletNIzXzQVXm6J4VuZjrAZze0Eb/BSHZmwDujcr1JvSmNyyqD7FRYKlRIxITR9wn2qOPPl2oIARMj5qp6YPDmpHUDs/gV0/5M/n6nPLZa0kvGGTwAg0T87V27BqSV2QAz+aWUAy/THpTpsSSFSoNGs/87yFypTjwmAu3F3Ykgjr9S98sUDb9Pg+RCTHftKtR03H5YWveCCJMD2mCPOv6sr+vpnsypCDKaudmSRkYyK6bHvMRlvzatV7zqFtMR/6COf1nVzBuIsNUKv7UBpJ5itUGzZTZqVHrvIO+VN6AtaSThHMTxm+wCV0yNDjaozWRfQTeD3G0vFarGZR7Yo45/qVAxrl0DA41tMFl3DYOXlSm5YzZj9kQt1ZrUB4hkTUb86USc8XRioIjfdMs/dHelDuZoTdk07jeVDXJHqxyfyaTlzIAly7b0dG/n9QCO32C+SLEnHf9itszFs2PGfr/qthwk4jftQAaterBxogYZzWxqXINi0CYuO45/8OakvMGKrHhMMDGSGOhzs1GmOKjUczZfAeDAVP8Ta4MwPZZivU+pbrlQYWFixKg9sEcd/wtrJS6eMej4A9wKm3YgnTrn/pyHmlI5lzEXfQSRDNR4BFN3bEEdv6ncDHSqjQZxaqYmcyoGTcznLTQmTfkz+fstoz6TryAExiPsmT4DnEazxeqG2QoxxZ5z/PlynfVKg8Mz5q7yEwFuhduO31BEpKLGfoeQdQZqma3jh0Edf42JdMLIGAkY3PG3d7YaLAls5x4GmCNkesLjoPmZrAWNf2YsRa3ZotTnCOSz+TIL42ljA9radmX6u3CvFGtIaf5CBHvQ8S9mSwBGpZ50Is5IMjbQ2AbzUs9gLe1K5zZZrRI04jeh7ysGdfztyhCDDk3ddQ2yKyBXMjt8bHDHb146VHc4/cqZZ9erHJwyL6v0K9UtFTzpyTl+A7ywVgYwKvWAd2IMqsmqnzPBoBumshs10okYo4Yiawg2k3/V0IpDhXL8/coFNjT+0ZQXSPT72TVb0p8+OzyOP1/Wv/h9MypJq5K2O3E2X+aABcffb2HF0rpKNjuNXzsq4jcp9cBgtbtg3oHEYmKgaFbVppuaPwPBI36TdyGzmRTVRv9yQb7s1ctnUuYcGqgVfv29TzYuRgNH/Bv69/9u5sCkd06rpO1OnMlXOGDByfZbWNGuMnIRv34Ws2XG0wnz3XoDdhGul+ukDI1AVsxlUqz0WaNuY/2bckyDSCsrRTO7bRWDdu+qbmuTF0iA2fH+a8HVnYGpyifwk7sDVK1lLeyRVbLNmT4cf7HaoFBpGK/ogf4LK5TUY/KOVrHnHP8LayUOz4yaP1EDVIeY7iCcH0+3yyF3wnQ3KnijmefHU5xb78+merNFtlQ3K/WMDe74TTcAgRdIrPZpU7vXweDnNzmapFJvUW30e2dkPpCYHkuSTsQ4my/v+Fx1V2BD4x9PJ0jExM5ST8G7m00lzLvlvef4syUOG9b3YfCaYhuOf268f+eRLdWNdu0q9k2MsLTe3625csYmFp4oZscHd/wmSzkVc5lU301A6jM2mQQftCLLRsQvhODg1EhfEf+5dVXDb97xCyGY6ePzW1q307wFe8zxSyn95i3zt3ezGU/q6Xc2jhXHn0n3PXlytVhl1vCJCt6Jd67Qn+NXS80XTFb1qAmdA0o9ppkZ4A5ydcN8RVY7Ydln3iFncBZ/NwemRtpOfTvOWIz4we/e3UHq8Zq3nOPXztpGjVKtabyiB7xb89YAs3FWi2ZHIIMXARarDSr17W/P6/5OYhsR//7JNGfz/Uk9naXmBkcgZwbrkrXl+Ocy3mfXj7Si5DyTA/bUrKSVPqRDKSW5Ut34HlmAg1OjfUX8Sg7abyG5C2rC6s5Szz4LXbuwxxz/C1m/lNNgDb+i03TTb5RWNZqMA9rLQnaKZlVtummNHzypZ3Wj2te8HpWYNun4J0cSxGOi/x2phjc4KdpNQH1E2Gsb3spME7ttFUpu68fxF6sNGoYH2Sn2T3oR/0532mfyFWbGkkaLKbrxdir0PqZaLcmypQFtsMccv61STugu4dr5RG22vNG+Jrc4QacRaKcEb9ZC85biwNQIUvbnQNRzTN4OCyH6XnXYaknWLUX8gySdbdw9qs9guY+aedOLj7o5ODVCvSl3DG7OrVesVPQodirvzpZqNFrSafwmaDdvWYj4Z9sR2s4naq5UoyXNl3GpZN9OOn97MqeFE3W/H+H0U3u9XKgymoyTMbTUXNFPIg6gUG3QkubnuUPX8dTHnYjX3Wz2WJpIJ0glYn05fhtduwqVrN3peDqTr1jT96GzhatXY2Cnht9JPdp5IVtiZixpZDnFZtRB3s+I3xWDo4a7me9Tl+3M6bEj9QB9lXSuFKvMT5i3qZ9EHHQ6jm1U9SjH30/Sec3wykzw7owWxtPt/QjbYWMWv6JTy799SefZfMVKRY9iZsxbx7rRozHQ1ohoxZ5y/IvZspVSThgsSdgeNWxwEiZ0Rfw7OA8bc3oUKrm21Edlj+kVhwovOttZ6rHRIasY5A5ydaNmPIgAT+4ZLOK3UyUGcHabyp5KvcnqRs1K166iLf32+PxUSbOTegywuFayUsoJkEnFScVjfbXZr6j6dMMn61gqwWgyvrPGv2F+drpiLpMiERN9ST0rBbNzehTeeISdHWx7r60Fxz89lkKInS/azZb0dtsaDiLAu4PspxPcpsY/n0nveDypmTg2I/5O9+7W/sBJPYZotfwafksRvxDCa+LqKxlnvsVeMTee6kPjrzOeTpBOmK94iMUE+ybS/Us9NiL+TJJcH3PdTY/S7iYe85POOxxP2ZI32td0EAH9R/y2dtuCdzztnxzZ1vGruwGbGv/MDjma5UKViXSCUcMznxR7xvEvF6vUmi0OW0jsKmYzqb40/tVirX1im2ZuPN2+w+hFtlQzOod/M/smR3aUehrNFmulmtHmLUW/c91tOjTor4lrzULzlmJhPMXaRpXmDqWT2ZK3Q8H03HvFTt27Sv+36vh3GNS2VKiwYEnfhz3k+F9Ys1fKqeh3UNvqhjejI2Zo+Xs385nUjlLP2kbNSkWP4sAOERp4eQcpzZZyKvotnbSp8UN/EtSKpXwReJ9FS+78PuVKNaYtBhIHpka21fhtrVzsptPp3EvjtzeuAfaS41cLWCxJPeA1cfWjFS8XzFdhKPqTempWKnoU+yfTO7bZK0nBTnJ3ez1WkS/XScaF0Z0F3cxmdp7QqY43G8ld9VnsJPfkynXjI5m7OTA5wpl8uadUdyZfYTydsFLdp1DBQa9afptdu7CHHP+iX8NvO+LvpxFodcOOdg2e1LO6Ud1Wv7Yd8e+bHGG90qC8jbTS7tq1EBXtpMcqvHENZncWdDOb2XnIno3JnIqFPrt3bQxo6+bA1AiVeqvnADmvecuekwVIxGNMjiS2VACklCwVKi7iN8FitszCRNpaizb0P6httWin/A48h1Bvym3XQnpVIXYjNGDbqH/FYsTf7+IMr2vXXtQ46zeWbXfRXi1WEcJOBU2/3bs5C7sdujnoSzi9dH7bzVuKmczWgWCh2qBSb1mr4QeDjl8I8UEhxJIQ4pGux35VCPGQEOIBIcQnhBCHTL3+Zk7nyxyathftg3fy9TOobbVYtaLJQsdx9tL5K/UmpVrTstTTh+NvD2iz08AFO9fM58pm99puZjaTotHa/qK96t+txW3ki/psCMxFEPFD71r+s5Y2b23GG9R2/jHVXrl4gUg9HwLesOmx90gpb5BS3gT8PfCLBl//RZzKlblo2u6HPdvHoLZyrclGrWkv4t+hicvmnB5Fe2zDDo4/nYhZ0WV30mMVtiZzKvpp4lot2rtby6QTjKXi20b83v7fupU7EMXBbcY2NJotlgoRRfxjyS3zRucsN2+BQccvpfw8sLbpsfWuLzNAf8Pqw9vCmVylfQtoi34Gta1uqBnzljT+HQa1qYuUzVvz/f5JuLRNLf9KscbCRNqKnp6Ix5ga3Xlsg23Hr+7CtisRXrPUtauY32FsQ75cR0o7XbuKhYk0MbG11LNcrNKSdit6FDM9Iv7PPr5EMi645tCkNVvsCZQ+Qoh3AT8A5IFX2XjNXKlOud60LvXs1KYN9ub0KJRU0qvjUo39tRnxT6S9juKdpB5bCXBQ83p2iPhLdiPZdpnpNlVZKxtVrj5gz4EsTKS3lXpUMtOm1JOMx1iYSG+5gvHkqlfdd5HFIg/F9BYRf7Ml+bsHT/NNV+6zeixZT+5KKd8ppbwY+AjwE72eJ4T4YSHEcSHE8eXl5VCveSrnHQD2pZ6dIzSbXbvQiRp7lXR25vTYO1GFEN5Clm0c/3LBsuPPbN+D0WxJCtWGlQFtin6OJ/sRf2pbqSdrcVxDNwd6LGR58lwBgKv2T1i1B7xAsFhtUGt0dk/c8+wqS4Uqb77ZWroTiLaq56PAd/X6ppTyA1LKW6SUtywsLIR6oc6aNcsRfx+D2myW34EXDU2PJdsS02ayEUg94CV4t5d6qixYmMypmNmhWapQ8SSMKDT+XnbVmy1ypbrVuzUv4u/9PuVK0RxPBybTW2r8T5wtMDmSaOeVbKICl8fOdBTvv/3aKcbTCe66er9VW6w6fiHEFV1ffhvwuI3XPe1H/LalnkwqTjIuth3UtrKhIn57J8ZcpncT19pGDSHsOjTwNyf1GNvQWVRj72Td6ra8G9tduwBjqTjpRKyndKj0Y1t3jwAL4yOsbdR6blBrD2izfDwdnBrd0vE/ea7AVQcmrPVedPPN1x9g30San/2rh6g2mlTqTf7pkbO8/toDVsvMwWw558eALwNXCSEWhRDvAH5DCPGIEOIh4HXAfzL1+t2czpdJxWPWomqFEMJL6OxQhTGWijOWspdumRvvrctmS16Josm1fVux34/QtqpRX9uws6imm9keiTiFzcmcCiHEtk1c6mI+b/E4V/sRegUSNpewdHNgaoRCtUGx2il9lVLyxNkCV0Yg84And/36d17P42cL/PdPP81nH1+iUG1Yl3nAYHJXSvnWLR7+Y1Ovtx2ncxUOTo9YmYWzmZ3a7FctJy3B02WfPFfc2h7LXbuK/ZMjVBst1suN86Zd2liyvpmZTIpSzYvKtorGbE7m7EY1cW2Fcr5WpZ6uWv6tumFzpToxARMjdutIOiWdZS7f5zn6c+tV1isNrjoQjeMHeM3V+/nulx3m/f/yDFfun2BhIs0rLpu3bsee6Nw9nStzKILyLfAkg+0c/4rFrl3FXCbds5wza3lOj6LdxLWF3GNj1+5mdprXE4XUA9tPfF3dsFsoAJ0RGr0SvOoO0nbQdfm+cQAeWsy3H3vCT+xGFfErfuFbr2HfRJqvn1nnTTccstJst5k94fjP5MoctFzRo5jdYZTuisWuXcXcuNc63thCl31+rWQ9FwIdx7+VLmuza1cxu8PYBtsjmRXbHU+2CwWgE/H3quXPlevWZR6Aqw9MMj+e4gtPrbQfe/LscDj+qdEkv/ndNzA5kuAttx6OxIYL3vE3mi3Orle4KAJnBmo083YNXDWrDg06EeHmyHGj2mAxW+ZKP1qyiaqy2KqWvz2Z02rEv30PRlQR/3bVRmsb3l4HmzbtNK8nV6pZreFXxGKCOy6f5wtPrbRnZT1xrsDCRNqqFNaLO69Y4MFfeh0vsdhz0c0F7/jPFbxOPdulnArVrbfVoLaWX61iW+qZ71HL/9SSp/tfGYEG2tm9e74DWSnWSCViTFgco6sW0fRq4lov10knYtarMWYzKQqVF9eCK1Y3qsyM2dnroBhJxplIJ3oXC2zYbXLr5s4rFlgpVnncj/SfOleIpH6/F1FUFikueMd/pl3KGY3UM5PpPagtX67TbMkIpB41tuHFjv/JCDXQkWSchYk0Ty+dn3ReKVRZGLczrkGxk9Rje1yDYnab3pDVov27R9h+BWO+bHdAWzd3XuElTb/w1DKtluTJc8XIZZ5h4YJ3/J2u3Wgi/tltIsfVCGr4u19vcxPXk2cLpBMxjlhcT9nNbUdnuefZ1fNKOpeLVesOrR+pJ0rHv1VJ5+qG3XHaivnx3o4/a3kkczf7Jkd4yYEJvvDUCovZMuV6k6sO2Jcxh5EL3vG3u3Yj1Phh627L5YL3mK0BbYr5jCrB2xTxLxW5fN94JFUGALdfNseZfIXn/TWZinPrFeslr6lEjEwq3lPqyZWidfxbXZA82dB+R2qveT3Vhjfi23bzVjd3XjHPfSfWeGAxB0Sf2B0WLnjHfzpXZnLE7pq1brYb1BZF+R3A5GiCREycd7I+dS665haAlx+bBeDLz6y2H1tar/DkuSI3H5m2bs9283qiivjV6N7F7PkDyLwKsSgi/q3n9SgpcTrCZOqdVyxQa7T4s3tOAnCFc/zAHnH8UZQnKrYbrLVqeTKnQgjB5fvGeeD5XPuxfLnOmXwlUsd/2cI4CxNpvvxsx/F/5vElwGt8sc3MWO+a+Xy5br15C+DoXIb58RR3P7Pyoscr9SaFSiMSqWdhIs16pUG18eLVmXf7F/CbDk9bt0lx26WzpBIx7ntujcMzo5EFgMPGHnD8lUgd/3aD2tSavCg00FdetY+vnFij4Cedn15Sid3oNFAhBLcfm3uRzv+pry9x0fQoL4mg0mh6m9HMuVLN6gJxRSwm+MbL5/ni0y/OhXz1ZBYgkq7UziauFx/jn/76OfZPprnuomhKFsErGviGS707yWGq6ImaC9/x58uRVfTA9oPaViyuydvMq65aoNGSfOlpLyp74qxfyhnxyXH7sVnOrVd5bmWDSr3Jl55e4TVX74uk9G22h9STL9fZqDUj2eIE8I2Xz7NSrLY7UQE++8QSqXiMOy633/6vdsWeWNloP1apN/mXJ5d5zdX7Iy1bhE51TxRlysPKBe34S7UGuVI9shp+2H5Q20qhal3mUbz0khkm0gk+94QnpTx5rsBYKh5Z9ZPi5cfmALjn2TW+/Mwq5XozEpkHejdLRTXtVaEc2Re7ulI/8/gS33BslkwEUsY3XDrH5EiCj977fPuxe55dpVRr8tqIPrtuXnXVPoSAGw9PRW3K0HBBO/7TOa+iJ2pnNtdjPd2pCPMPyXiMO6+c53NPLCOl5KmlAlfsG49kkF03l85n2Ofr/J/6+jnGUp1bddtMjyUpVBrnjRw+5SdWo9jiBF4z4mULmfY4gudXSzyzvMGrrtoXiT2ZdIK33naEjz9yhsWsV5H16a8vMZqM8/LL5iKxqZsr9k/wuZ96Ja+/9kDUpgwNF7jj907QqG7JFUfnxniu6zYYvBGxJ1dLHJ3LRGSVp/OfXa/w+NkCT5wdjuYWIQQvv2yOLz+zymceX+LOK+atd8cqOs1SL5bpou4NAa9a5d7nVqk2mnzm8XMAvPol0Th+gLe/4ihCCD589wmklHzq6+ci/ew2c8lcJnLJaZi4oB3/mXy0t+SKYwsZnl8rvShyXN2oUaw2uGQummYpgFde6W02+5uvnWKlWB0Kxw+e3LNSrHImX4lM5oFOE9dmnf90rkwqYX+/QzffePk8lXqLr57M8Zknljk2n+HofHRBxKHpUd543QH+/L4XuO+5Nc7kK9x1TfQyj2NrLmjHfypXQQi2nBNuk2Pz4zRb8kWNSSdXvTuAKCP+fZMjXHtosq3NXhFhRU83t/s6vxDRRrGdsQ0vjvgXc2Uumh6NVBa7/dgs8ZjgE4+d5Z5nV3lVhO+T4h13XEqh2uCn//KhyD87x/Zc0I4/u1Fj30SapOVtUps5tuA592eXO3LPiRXvIhBlxA/wyqsW2luKolxQ0c0lc2McmhrhpounrXfsdqNmzGxO8J7KRlspBjAxkuTmi6f5yL3PU2u0ItP3u7n5yAwvPTLN82slbo74s3NszwXt+H/1zdfx+Z95VdRmcGzei6SfXe4MIDu5ukFMwOGZaB2/chgT6QQHJqN1ZgohBL///S/jPd99Y6R29OrBOO1H/FFzxxXz1BotMqk4t0WUAN/MO+44BkTTcOfonwva8QOkE9Enl6bGksxlUi9K8J5YLXHRzCipRLQfwU0XTzM5kuCK/eNDlfy64fB0e4tSVMyPez0WKpkL3vyZpUKVi6ajvWAD7Zr9O66Yj/w4UrzhugP8+ndez/e//JKoTXFsg+tftsSxhcyLpJ6TqxuR6vuKRDzGu7/z+kjmzgw76UScyxYyPHp6vf2Y2hAWtdQDcOPF07z6Jft42+3D42TjMcFbbzsStRmOHXCO3xLH5sf5tF92B17E/6YbD0ZoUYdvveFQ1CYMLdcdmuJLXXNxoq7h7yYZj/HBH7w1ajMcu5DhuD/cA1y6kGGlWCNfrpMreX8PQ8Tv2J5rDk1ybr3anj65OAQ1/A5HWJzjt8Qxv8b6uZUNTqyqih7n+Iedaw95bf6Pns4DXmJXiOhWeTocOnCO3xLHFjqVPZ0a/ugThI7tueaQN1lS6fynsmX2TaSHJpnqcATBafyWODI7RjwmeHZ5g2Q8hhBwcUQrDh39MzWa5OLZUR5Tjj/i/Q4Ohw6c47dEKhHj4plRnl0pMpKIc3ByZGjmmDi257pDUy+Seq67yE15dOxu3P2qRY4tjPPs8gYnVjecvr+LuPbQJCdWS+TLdU7nKkNR0eNwhME5foscm8+0k7tH553Ms1tQCd4vPLVMrdlyFT2OXY+Teixy6UKGaqNFtVFzEf8u4lo/wfvJx7w+DOf4HbsdF/FbRM3sAVfRs5vYNznCwkS6vfjdJXcdux3n+C1y2UInyncR/+7i2kOTFCreFFOn8Tt2O87xW2RhIs24vxM16nHMjsFQcs/ESILJETfXyLG7cY7fIkKI9k7ZsZRLr+wmVILX6fuOCwHnfSzzPbccZqVw/uJ1x3BznXP8jgsI5/gt8wMvPxq1CY4AXDw7ysJEOvIdAQ6HDpzjdzj6QAjB3/+HO5gYcaeMY/djTOMXQnxQCLEkhHik67H3CCEeF0I8JIT4GyHEtKnXdzh0s39yxOVmHBcEJpO7HwLesOmxTwLXSSlvAJ4Efs7g6zscDodjC4w5finl54G1TY99QkrZ8L+8Bzhs6vUdDofDsTVRlnP+O+DjEb6+w+Fw7EkicfxCiHcCDeAj2zznh4UQx4UQx5eXl+0Z53A4HBc41h2/EOLtwLcC3yellL2eJ6X8gJTyFinlLQsLC/YMdDgcjgscqyUKQog3AP8P8E1SypLN13Y4HA6Hh8lyzo8BXwauEkIsCiHeAfweMAF8UgjxgBDi9029vsPhcDi2RmyjtgwNQohl4GSAH50HVjSbE4ZhswecTf3ibOoPZ1N/2LLpEinleVr5rnD8QRFCHJdS3hK1HYphswecTf3ibOoPZ1N/RG2Tm87pcDgcewzn+B0Oh2OPcaE7/g9EbcAmhs0ecDb1i7OpP5xN/RGpTRe0xu9wOByO87nQI36Hw+FwbMI5fofD4dhjXBCOXwghorZh2BnG92hIbbogzgmTCCGGbimBEGLC/3vojqlhZNce5EKIa4UQrwTYbuaPTYQQB/2/41HbAiCEuFoI8XIYqvfoOiHE64UQiSGy6XohxE8CSClbUdsDIIS4TQjx7mG6EAkhXi6E+EPg1qhtUQghXiqE+EvgHTAcx7kQ4iYhxA8JIQ5EbUsvhu7KvRP+ifB7wKuB54UQrwH+t5TyuBAiFsWJK4QYB94PfJ8Q4kYp5cNCiLiUsmnbFt+eKeC3gNuAZSHEvcCfSCmfjsIe36YZ4F3AK4BngLuEEL8vpXwmKpu6eBfweiHE/VLKz0X82U0Cv47nXD8kpWwJIUTUDk0I8UPAfwTeB3wtyvfIt2cO+GW892kWb78HEX92STzfdAvwdeB2IcQHpJT3RmHPdgxNNDEAM3jzfq4Gvg9YBX5SCDEeYbT2JuAF4HfwLgBEeVIAP4NXsXUj8CPAHHA0QnvAs6kqpbwJ+PfAtUCkt+VdksXngfcCvwbeZxdhpP1O4HbgdVLK9/n2RB7FAkeAd0op3y+lrER8fIMX2Egp5e140f73Q+Tn3fXAlJTyZVLKt+H512EbFQHsEscvhHitEOK1/peTwMuBMSnlMvBXeJu+/i//uVaciW/T6/wv/wH4b1LK/wIcEUL8G/851u6oNtnzfuAXAfyIehrvoLSKb9Pr/S9/WUr5n/1/vw4vSrtWabOWbboLQErZ8I+X1wN/CCwJIf69/72W5WNJvU8fBJaAfUKI7xZC/JYQ4t8IIY7YsGWTTXf5/57Cu1DfJ4R4tRDin4UQPy+E+E7/+1Gccz8qpfyP/r+XgceEEFfZsGMLm5RvagJvEUJM+e/N7cBrhBA3+88dmvzDUDt+X8f/c+DngSyAlPI54EvA/+0/7Qzw18DNQohDpqOjTTat+Tat07my/xfgN/3HG1v+EjP2vLPLnkUp5emuC08ZT16xwqb3aNW3qep/718BPwV8GPgO4BeFEMZXcG6yKec/pvIMD+Ddsf0a8NNCiP8lhDhs+VhS79MTeLLFx4EfB54Avse3K5L3SUqZB/LAnwFvxpN7zuB9djdafJ+6j/Fq112ZBA4BJf/5xh1sD9/0IPAbeO/P7wPvBi4GfkUIceWQ3Ll5SCmH8g9eRLgKvG+L770K+EfgUv/r6/CWu18WlU3+91VD3N3Af/X/PRKhPQn/708AN/v/jkX5Hm167jV4F4BXRXgsjeE5tJcAv43n7P7F/17cpk3qswFGgR/c9D79CXCn7fepy6YjeJH1z3Z9793ALw7D8QR8BvgJ/98iwvdJAL8A/Cv/6zngd4F/a9KmQf8MXcQvhDjoJ2jW8E7EtP/4D/rVIJdIKT8LfA14D4CU8hHgEqAagU2vE0Ic85+qqnneDPxHIcQvA+8VQuyPwh7pSRmXA2tSyq8JIX4M+AUhxLROewaxSfj49j2GN572hG57+rTpctlZCPQVYByvaOCIEOIGaUAv3s4mvIT3FVLKMt4FEWi/TweA53Xb06dNl0spn8cLrr6760f34QU5Vm3adDwpH/a/8OSxuPQ9rk2b6Hx2Em+t7FsApJSrwEXAYyZsCkzUV56uK+ZrgC/gafZ/4D82CjxIR855L/AQcBmQwlv08t+BR/AOygk0Xu37tOl3fJsu7/q5a/H0vs8B10dkzxX+91+LJ/N8Fk8+uCqCz639HuFFRAng24FP42naUXxu78WTeF4CfAtwZdfPfz9wJMr3qevnvs1/n/44wuP7EeCo//2/wpMz7vG/fzDC9+mKrp/7NbzKNW22BLDpYf8Yv9Y/597j/9yfAwsmbAv8f4raAP9NvBK4Fy+a2Oc7qNf433sj8Pau534Q+P/8f+/HKw/8toht+mPg1/x/H8bT9/51xPb8uv/vt+HponcN0Xt0F3Af8OaIbfoT4Fe6vo5hQAoL8T69Arh/CN6nP+465ybxLpavGwKb3tX19aXAG4fgeFKf3c14FXXfodsmHX8iq+NXt2jSK8G8CbhPSvmXwqtj3gCeEEIkpZQfV8/3n/sPeI4DKeU54NwQ2PSPeNl7IaVcBH50COy5y/81H5VS/pkOezTY9Br/Zz8FfGoIbPp7ZVPXz0dtU/exdDfwsmGxyf/ZdWAdeHxYbPJ//jnguYht+ns6vulreHL0cBLF1Qb4t8Bp/Cs2XgRxN15J3Um8uuo/A/7U/75Kmr4d7xbv2y90mzTY8+YhfI+G0aZhPJbc++RsMvrH/gt6CbS/Bf4T8FV8zRlYAH4a+DH/6xG8KoLb/X//Np5mfuuFbtOw2eNscjY5m6KxydSfaF7UT5zhJYk+5v87hqfb3dn1vN8D3uB/75K9ZNOw2eNscjY5m6Kxycj/M9IX90rU7sNPygA/hzfj4iq8xogHgGN72aZhs8fZ5GxyNkVjk84/kW/gEkL8CPA2KeWd/te/BRzEu5L+jJTyhb1u07DZ42xyNjmborFJF5E6fpUNF95Y1SW8luu/AB6WXhPLnrdp2OxxNjmbnE3R2KSTSDt3/Td2DK8+9i3A81LK+6J8Y4fNpmGzx9nkbHI27X6GYR7/j+Nl0F8r/UFeQ8Cw2TRs9oCzqV+cTf3hbLLIMGj8kSxP2Y5hs2nY7AFnU784m/rD2WSXyB2/w+FwOOwydNM5HQ6Hw2EW5/gdDodjj+Ecv8PhcOwxnON3OByOPYZz/A7HJoQQ00KIH/f/fchv4nE4LhhcVY/DsQkhxFHg76WU10Vti8NhgmFo4HI4ho3fAC4TQjwAPAVcLaW8zt+t+ma83crX4Y3jTeGtaqwC3yylXBNCXAb8D7xxviXgh6SUWhaXOBw6cFKPw3E+Pws8I6W8CW8OezfXAd8L3Aa8CyhJKW/G2//8A/5zPgD8Bynly4CfAt5nw2iHo19cxO9wDMZnpZQFoCCEyAP/x3/8YeAGIcQ43q7c/yWEUD+Ttm+mw9Eb5/gdjsHontnS6vq6hXc+xYCcf7fgcAwlTupxOM6nAEwE+UHpLSN/TgjxPQDC40adxjkcYXGO3+HYhJRyFfiSEOIR4D0BfsX3Ae8QQjwIPAp8u077HI6wuHJOh8Ph2GO4iN/hcDj2GM7xOxwOxx7DOX6Hw+HYYzjH73A4HHsM5/gdDodjj+Ecv8PhcOwxnON3OByOPYZz/A6Hw7HH+P8B61hCWSW5C3cAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot the first 100 time steps\n", "ds_global_avg.tas.isel(time=slice(0, 100)).plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3. Tropical Region" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "ds_trop_avg = ds.spatial.average(\"tas\", lat_bounds=(-25, 25))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'tas' (time: 1980)>\n",
       "array([25.24722608, 25.61795924, 25.96516235, ..., 26.79536823,\n",
       "       26.67771602, 26.27182383])\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n",
       "    height   float64 2.0
" ], "text/plain": [ "\n", "array([25.24722608, 25.61795924, 25.96516235, ..., 26.79536823,\n", " 26.67771602, 26.27182383])\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n", " height float64 2.0" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_trop_avg.tas" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEiCAYAAADu2jXpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABnsElEQVR4nO29eZxkZ13v//7W3l3Ve/esmSUrySRAQoYQViEkkUQviwqikhv9KbiAghf1x9WrFxeUH7IoV0VBULyCCBIFBYSAEVkTkpB1hpDJZDJbT0/v3VXVtT+/P855qqqrT9U5p7qqTk3P83695jXdVXWqn6o6db7Pd/t8RSmFwWAwGAxeCQW9AIPBYDCcWxjDYTAYDAZfGMNhMBgMBl8Yw2EwGAwGXxjDYTAYDAZfGMNhMBgMBl8Yw2E4ZxGRYyJyY5vHPioiL+723+lXRESJSEZE3tHm8TeKSFpEKlvtvTG4YwyH4bxEKXWlUuo/N/s8IvJiETnZgSU5Pfe7ReRxEVkVke+JyH93efxPishTtkH4FxEZd/kTz1RK/VY7a1NKfVkplQKOt3O84dzGGA6DoX/JAP8NGAFuB/5URJ7n9EARuRL4K+A2YDuQBf6iR+s0nGcYw2E417laRB4SkWUR+UcRSeg7ROSHReQBEVkSkW+KyDPq7quGn0RkQEQ+KiKLInJYRH7DwYvY8HdEJAl8Adhlh23SIrKrUy9MKfW/lVLfU0pVlFJ3A18Dntvk4T8F/KtS6r+UUmngt4EfEZEhL39LRPbb4aufEZET9nvxCyLybPt1L4nIn3XmlRnOdYzhMJzrvAZ4GXAh8AzgpwFE5FnAR4CfByawduOfFZG4w3P8b2A/cBFwE/A6L39HKZUBbgFOK6VS9r/TjQeKyNvsC6/jPy8vUkQGgGcDjzZ5yJXAg/oXpdQTQAG4zMvz1/Ec4FLgx4E/AX4LuNF+/teIyA/4fD7DFsQYDsO5zvuVUqeVUgvAvwJX27e/HvgrpdTdSqmyUuqjQB643uE5XgP8oVJqUSl1Eni/j7/jilLqnUqp0Wb/PD7NX2IZhi82uT8FLDfctgx48jjq+H2lVE4p9SWsUNk/KKXOKqVOYXk81/h8PsMWxBgOw7nOmbqfs1gXUIB9wFsbdvZ7AKdQ0i7gRN3vJxwe0+zvdB0R+WPgKuA1qrkqaRoYbrhtGFj1+edm6n5ec/i9Z6/b0L8Yw2HYqpwA3tGwux9USv2Dw2OngQvqft/j4++4ykuLyG/W5UA2/HM59nexwmE3K6VWWjz0UeCZdcddBMSB73t7GQaDd4zhMGxVPgT8gog8RyySIvJDTZLFnwT+p4iMichu4E0+/s4MMCEiI80eoJT6w7ocyIZ/zY4Tkf8J/CRwk1Jq3mUdHwP+m4i80E7a/x5wh1LKr8dhMLhiDIdhS6KUuhcrz/FnwCJwBDtx7sDvASeBJ4EvA/+ElQ/x8ne+B/wDcNQOiXWsqgr4Q2Av8Hidh/Kb+k779xfa63gU+AUsA3IWK7fxSx1ci8FQRcwgJ4NhPSLyi8BrlVJbtoJIRHJYxvH9SqnfbuP4lwKfxgqH3aqUuqvDSzT0McZwGM57RGQnVinut7BKUT8H/JlS6k+CXJfB0K9Egl6AwdAHxLD6PC4EloBPYLquDYamGI/DYDAYDL4wyXGDwWAw+MIYDoPBYDD44rzIcUxOTqr9+/cHvQyDwWA4p7jvvvvmlFJTjbefF4Zj//793HvvvUEvw2AwGM4pROQpp9tNqMpgMBgMvjCGw2AwGAy+MIbDYDAYDL4whsNgMBgMvjCGw2AwGAy+MIbDYDAYDL4whsNgMBgMvuia4RCRPSJyl4gcFpFHReTNdff9sog8Zt/+Lp/Hvl1ETonIA/a/W7v1Grzwbw+d5vaP3EO5YjS/DAbD+UE3GwBLwFuVUvfbU9fuE5E7ge3AK4BnKKXyIrLN67FKqUP2/e9TSr27i2v3zDeOzPHV78/ytcdnefHTnF6KwWAwbC265nEopaaVUvfbP68Ch4HdwC8C71RK5e37zvo4tu+YSxcA+MQ9JwJeicFgMPSGnuQ4RGQ/cA1wN3AZ8EIRuVtEvioiz/ZxrOZNIvKQiHxERMaaHPcGEblXRO6dnZ3tyOtwYi5tTRj98uEZZlc9TRs1GAyGc5quGw4RSWGNmHyLUmoFKzw2BlwP/DrwSRERj8cCfAC4GLgamAbe43SsUuqDSqmDSqmDU1MbNLo6xny6wNN3j1CqKD59/8mu/R2DwWDoF7pqOEQkinXh/5hS6g775pPAHcriHqACTHo8FqXUjFKqrJSqAB8Cruvma3BjPp3n2fvHuW7/OP/4nROYwVgGg2Gr082qKgE+DBxWSr237q5/AW6wH3MZ1tjOOY/H6vnQmlcBj3R88R5ZK5TJFMpMpGK89ro9PDmX4dtHF4JajsFgMPSEbnoczwduA25oKJ39CHCRiDyCNdv5dqWUEpFdIvJ5l2MB3iUiD4vIQ8BLgF/t4mtoic5vTKXi3Pr0nQwnInziO8eDWo7BYDD0hK6V4yqlvg445i6A1zk8/jRwq9uxSqnbOrXGzTKfsSqqJlIxEtEwr7xmN5/4zgneVSoTj4QDXp3BYDB0B9M5vgnmbY9jIhUH4NLtQxRKFZbXikEuy2AwGLqKMRybQIeqJlMxAIbilgOXzpUCW5PBYDB0G2M4NoFu/ptIWh5HShuOvDEcBoNh62IMxyaYTxdIxsIMxKx8RiphPA6DwbD1MYZjE8yl80wOxau/a49j1XgcBoNhC2MMxyaYz+SZSMaqvw8nogCsGo/DYDBsYYzh2ATz6UK1ogrqQ1WmqspgMGxdjOHYBHPpPJN1hiMZt3IdJjluMBi2MsZwtEm5oljIFKqluADxSJhYJGRyHAaDYUtjDEebLGULVBTrchxg9XKYqiqD4fzlkVPLvP7v7uXeY1tXt66bEwC3NLqHo76qCqw8hwlVGQznH/lSmf/zlSN84KtPUK4onrZ9iIP7x4NeVlcwhqNNqnIjyQbDYTwOg+G8I18q8yN/8U0ePb3Cjz7rAu48dIbFbCHoZXUNE6pqkzlb4LA+xwGW4TA5DoPh/OKJsxkePb3C7/zwAd7zmmcyORRnaQtr1hnD0SaNAoeaoYTxOAyG842Ti1kADu63JlmPDcZYMh6HoZG5dJ5wSBgdiK67fSgRZTW/dXcaBoNhIycX1wC4YGwQgNGBKEvZrXsdMIajTebTBcaTMUKh9WNDTI7DYDj/OLm4xmAszNigtZEcHYwFbjgen1nllj/9Gvc82fnqLmM42mQuXdhQigu1qioze9xgOH84uZjlgrEBrKnXMDoYDTxUdWppjcPTK4S7cJU3hqNN5tJ5phpKccHyOIplRb5UCWBVBoMhCE4urlXDVGCFqjKFMoUArwNnV6w87LahRMef2xiONmkUONQMJcxMDoPhfEN7HJpR+9qwtBac1zGzkgNg2/DGDe5m6ZrhEJE9InKXiBwWkUdF5M119/2yiDxm3/6uJse/zH7MERF5W93t4yJyp4g8bv8/1q3X0IpGgUNNykwBNBjOK5bXiqzkSusNh100sxxgnuPsap7RwSjxSLjjz91Nj6MEvFUpdQVwPfBGETkgIi8BXgE8Qyl1JfDuxgNFJAz8OXALcAD4CRE5YN/9NuArSqlLga/Yv/eUbKFEtlBeJ3CoMVMADYbzi1MNFVVgleMCLAZoOGZWcmzvQpgKumg4lFLTSqn77Z9XgcPAbuAXgXcqpfL2fWcdDr8OOKKUOqqUKgCfwDI22P9/1P75o8Aru/UamjGvR8amnJPjYGZyGDbPqaU1zq7mgl6GwQXdw7HO47Crq4JMkM+s5rsSpoIe5ThEZD9wDXA3cBnwQhG5W0S+KiLPdjhkN3Ci7veT9m0A25VS02AZJ2Bb1xbehDm7+a+xaxxgKG6dMMbjWE+uWGbFzCnxjFKK2z58N//7M48GvRSDC409HAAjdqgqyO7xsys5tg+fYx6HRkRSwKeBtyilVrD0scawwle/DnxSdA1b3WEOT+WrvlVE3iAi94rIvbOzs22svDlVgUOnUFXV4zAXyXp+/98Ocdtf3x30Ms4ZHptZ5ehshtnVfNBLMbjQ2MMBMKaT4wF5HJWKYnY1z/Zz0eMQkSiW0fiYUuoO++aTwB3K4h6gAkw2HHoS2FP3+wXAafvnGRHZaT//TsAp1IVS6oNKqYNKqYNTU1OdeUE2zeRGwFRVNePobIZD0yuUK6a/xQtffGQGMCHPc4HGHg6AZCxMJCSBNQEuZAuUKurc8zhsL+LDwGGl1Hvr7voX4Ab7MZcBMWCu4fDvAJeKyIUiEgNeC3zWvu+zwO32z7cDn+nKC2jBvC1w6NgAGDc5DifmM3mKZcXppbWgl3JO8MVHzwDGcz0XaOzhABARRgdjgSXHq6W4Dr1mnaCbHsfzgduAG0TkAfvfrcBHgItE5BGspPftSiklIrtE5PMASqkS8Cbgi1hJ9U8qpXSw953ATSLyOHCT/XtPmV3NMxSPkIhuLHOLR0JEw2I8jgZ0QcFT89mAV9L/nFjIcmh6hXgkZDYg5wCNPRya0cEoywH1cVSb/7rkcXRtHodS6us45yoAXufw+NPArXW/fx74vMPj5oGXdmiZbbGQKVRjmI2IiNGraqBcUSzYsd5j8xlecGljZNJQz5cOWWGqH7xyB//60GkqFbVBE83QHzj1cGhGB6IsZoL1OM65UNVWZjHb3HCAmQLYyGK2gJbuemo+E+xizgG++OgZnrZ9iKfvHkEpyBTMudSvOPVwaEYHY4FVVZ21iyqmHPKwncAYjjZYyhbXVVA0kopHTYihDh2mAjhmQlUtmU/nuffYAj945fZqoYU5l/oXpx4OzehglOWAqqpmVnJMJGPEIt25xBvD0QaL2UK1M9SJoXiEtJnJUUX3vYwMRI3H4cKXD89QUXDzlTsYSlibE2M4+henHg7N2GA0wOS4swhrpzCGow2WssVqZ6gTJlS1Hm04nrV3lKfms1RMSW5TvvjoDLtHB7hy17DpCToHcOrh0IwOxlgrlskVyz1f19nV7jX/gTEcvimUKqTzpZYeRyoeMbvEOnSo6tp9Y+RLFWaMjIYjlYrinicXeMnlU4hILVRlNiF9i1MPh0Z3jy/7zHM8dHKJ//utY5ta18xKrmvNf2AMh2+0THLLHIeZO76O+Yw1ZvcZF4wCcGzO5DmcOLGYJZ0v8fTdIwAMmxxH3+PUw6HRm0u/TYB/841jvP1fD1EqtzfLo1xRzKULxuPoJ/RJMNoqx5GImF1iHXrM7oWTScBUVjXj0OkVAK7YOQxQl+MILlSllOI3/ulBvnGksUfXAM17OKB9ocOjcxnKFcX0cnue+XwmT7miutbDAcZw+GYxoz2O1snxQqlCvtT72GY/osfs7hodIBoWU1nVhEPTK4RDwmXbhwD6oqpqNV/ik/ee5HMPTwe2hn6lVQ8H1AyHnwS5Uoqjs2nAagRth9rkPxOq6hsWqx5Hq3Jc6wufyRvDAdYOaGooTjgk7BkfNB5HEw6dXuHiqWRVkWAgGiYckkA9jhl719vuRWwrM71sVVTtGm1mOKzNpZ/u8bl0obpRON7me97t5j8whsM32u1s3QBoS6ub2DRgT0u036/9E0njcTTh0PQKB+wwFfSHCoEOl7R7EdvKrKxZn8vogPO1QE8B9ONxPDlX21SdWGzT47Cb/0xyvI/Q0hmtGwB1NUywZZTLa0Xe9PH7ue3DwcqZz6fzVSXhfROWx6GUKcmtZzFTYHo5x4Fdw+tuH0oEW6F3xjYcpxbX2k7WblUydh5zMO48mnUwFiYWDvlKjuswVTwS4vhCe4KgMys5RJzHPnSKrmlVbVWWskVikRADDgKHmqq0eoBf+AdPLPHGj9/PycU1omFBKeVYMtht1gplMoVydVri/okk2UKZ2XSebV0aa3kucnjaSowf2Dmy7vahRJSVPvA4Snayds+4cwXR+Yju1dIbxUZEhJHBqK/k+NG5DLFIiGv2jrYdHpxZyTORjBMNd88vMB6HTxYzBcYGoy0vwkHPHf/3R87wY3/5TZSCH3r6ToplRb4UzG6xOi0xWfM4wKjkNnJoWldUDa27fSgeCTTHcWaltus14ar1ZG0NsWQTwwFWZMKvx7F/YpD9E8lNJMdzXU2MgzEcvlnMFltWVEHwc8fvuP8kU6k4n/uVF3D9ReOBrqU6u6TO4wA4NmcS5PUcOr3CjuHEhuFg/RCqGrfzU8bYrydtF7+kYs0Nx+hArNr75YWjcxkumkyxZ3yQ+UyhGg7zw8xqd5v/wBgO3yy56FSBtUuE4Dp+F7MF9k0kGR2MVXsBgvJ+Gqcl7h4bIBwScxFq4ND0yob8BliGI0j5munlHFfvGSUaFuNxNOCW4wDsUJU3j6NYrnB8PstFU8lqSLCdBPnZlXxXK6rAGA7fWJLqzRPjUGvcCirHsZApVHeJtYmEwYQ7tNyIrqqKhkPsHh3gmCnJrZIrljlyNr2uokozlIgGHKrKsWs0wZ6xQY4vmM+snky+RCwSaplL8BOqOrGQpVRRXDiZZK82HD4T5KVyhbl0vqvNf2AMh28sgcPWHkciGiIcksAUchezxapxSwWcqJ/L2DmOuhCMVVlldq+aI2fTlCqqqcexmisFUoWWK5ZZyhbZOTLAnvFB43E0kCmUmibGNdb4WG+hKl2Ke9FUij12U6Hf93w+U6CiuluKC8Zw+EIpxdJa61kcEGz9fbmiWMoWGLeNW9BCefPpAslYmIFYzZ3fN2EuQvVoqREnjyOViFCqKHLF3hc36FLcHcMJ9o4PcrxPjL1Sijd+/H7ueuxsoOvI5MskW4SpwGoUzpcqnhRyj85ahuPiqSTjyRjJWNh3grw2a9x4HH3DSq5EuaJccxxgK+QGcLFeWStSUbWu1aF4sDMd6ns4NFOpBMtrRQoBVXr1G4emV0jGwtXwRD1B6lXpUtwdIwn2TQyykiv51l3qBstrRT730DR/841jga4jnS+RbJEYh1pzoBev4+hcmvFkjNHBGCKWyoJfwzGrJ/+Zqqr+QX9p3EJVYCc1A7hY6xO0muOohqoCynFkCtWKKs3kUMy+Lx/EkvqOQ6dXuGLnsONc8eEAPUZdirtjJFFN1vaDp3jG3lV/64m5QPM/mXypZSku1BqFveQ5npjNVIVAActw+EyO6w2ilnTvFl0zHCKyR0TuEpHDIvKoiLzZvv3tInJKRB6w/93qcOzT6u5/QERWROQtXo/vFlo6wC1UBZbHEUQ1zGKDJEotOR5QjiNdYCK5fvej8x31I2XPZx6bWeXyhv4NTZBCh9MNoSroD8Oh11UsK772eHCqvV4Mx4gPw/HkXIaL6g3H2CAnFtZ85bdWXZoSO0U3PY4S8Fal1BXA9cAbReSAfd/7lFJX2/8+33igUuoxfT9wLZAF/rnuIS2P7xaLPjyOoKYALmSsE1TnOGKREPFIKNBy3MlGj8P+fTZtPI5yRbG8VmQq5RyTDjJUNbOcYzgRIRmPVA1HPxQ1aOHFWDjElw/NBLaOTKFMyi3HMaBncrTeJK3misyu5rloKlW9be/4AGvFMnM+Nlg6yqE3HN2ia4ZDKTWtlLrf/nkVOAzsbuOpXgo8oZR6qpPra4clDzpVmqCmAFZl3+tKhocS0UBCHZWKYsEpVGV7HHOrxnBkqt3HzhegID3G6eUcO0Ysg5aMR5hMxfpCJXd62dJietlVO/iPx84GpqGV8ZDj0N/DJZcpgDoxXh+q2jvh38vL5EtEQkI80t0sRE9yHCKyH7gG0Gp7bxKRh0TkIyIy5nL4a4F/aLjNz/EdYzGjQ1XechxBfNlrIoy1NQa1luW1IqWKah6qyphQldsOcSjAueNnVnLsGKlJhvdLSe7MSo7JVJxbrtrBUrbIvU8tBrKOtIdQldfkuC7FvXhqfagKrGFRftfUbV26rhsOEUkBnwbeopRaAT4AXAxcDUwD72lxbAx4OfCpups9HS8ibxCRe0Xk3tnZ2c2/ECyPQwSGPSSerBxH77/si5kCsUiIwbryV6s0uPdr0cnvRo9jMBYmEQ0Zj4NaR3+zC1AtVBWMx7GzrpFs73h/9N9ML+fYMZzghZdNBRauUkrZOY7WoaqBWJh4xF0h9+hsmpDUvAygOpLWTxn0as69t6QTdNVwiEgUy2h8TCl1B4BSakYpVVZKVYAPAde1eIpbgPuVUtUzw+vxSqkPKqUOKqUOTk1NdeT1LGaLjAxECTtUvzQylIiSK1Yo9tiNXrR7OOp3HEHJVujYbKO8s4gwmYpXBRDPZ9wUVoMKVRXtDmQdqgLYNz7I9PJa4GXUMytWCC0Vj/Dciye48/BMzxskc8UKFdVa4FAznoyx4OJdH5vPsmt0gHikZogGYmGmhuK+KqvS+WLX8xvQ3aoqAT4MHFZKvbfu9p11D3sV8EiLp/kJGsJUPo/vKIsedKo0tSmAvf3CL2SKG4ZMBZVvqcqNpDa+Z5OpuAlV4R6qCoeEZCzc88/v7GoepVhnOPaMD1JRcHqpvTkRnUJ7HAA3HtjOU/NZnrDnWPQKnZvysrv3YjhmV531pfb6DA+m8+e+x/F84DbghobS2XeJyMMi8hDwEuBXAURkl4hUK6REZBC4Cbij4Xkdj+8FltyIt/rooBRyF7MFxhu0tFIB5TjmHeRGNJOpWLVZ6XzGLVQFwehVnVmu9XBoqpVVAeY51gpllteK1XXdeMU2AO481Nsucr0hdEuOg2U43DZJcw7VhwB7xgZ86VWlc6XqtaebdO0vKKW+DjjFdBzLZ5VSp4Fb637PAhMOj7utU2v0y2K24Fl1ctiOTa/0+Au/mCls0DwaDkgoby5t5YScvLTJVJwHTiz3fE3NuO+pRTL5Ei+6rDNhTa+4haogmFCj7pXYWR+qsiXxg0yQ6+Y/7XHsHBlg9+gAj8+s9nQdNYPfOscBlsCnm6jnXDrPdReOb7h9z/ggn33wNKVyhYiHwUzpfKknw7ZM57gP/HgcwwPWhUDPJe4VCw7hNN2M2Os48Hw6z/hgzDEnNJmKs5DJU6kEP0J2JVfkDX93L3/wuUM9/9vVUFW8+XkVRFWc1qnaOVyrqto2FCcWCXE8QGXjMw4GbXIo3vOeoIw9i8NbjiPOQotejGK5wmK26OiZjw7GqKja33MjnS+d2zmOrYifHEcQHoduJmvMcQwlIlQUrHkQWusk8+mNPRyaiZT1hfCqHNpN/vyuI8xnCoFIxHjZuQYTqsqRiIaqGyCAUEiYSsUD7fjXMijb6wzHVCre87BnxsP0P81EKkamUG4qdKjzH5MO+lK12T7ePv90zr23pBMYw+GRXLFMtlD21PwHNa2YZZfGn06yvFZEKRgf3JjjgN7nWxayhaZd9v3Sy3F8PsvffP0YIsHoQWXyJRLRUMswRBA5qumVHDtHBjb0AwwlIoHOQD+zbBmIHXUh46mh3lfoZXxIe2jduGbnul77lMMmq6o15+HcLFeU1c1uPI7+QddhN+7mm1H1OHpoOBaqXeMbQ1XQe8ORyZeqO6ZG+qV7/I++cJhwSPjRZ11AJoBw3mq+RKpFmAosocNeX6zP1FUurV9LwIOlltcYsmVQNFMpK/ncyw7yjIeiBo02HM3CVc3K1qH23fXiDfup9NosxnB4ZNGhI7sV2ur38gvfqIyrGQ5I76iVCFw/6FV9++g8X3jkDL/04ou5aCpJRdHzuRfpXMlV72goEe15M+mZOrmR9WsJeAb6Sm5dfgMsj0OpmmpCL9DzxpMxb8lxaK4GrTdPjeMHoC5a4MHj6JVOFRjD4ZmawKG3UFU4JJZbH4TH0Zgc9+HudpJMofmgm35QyP3zu46wcyTB6190UdUz6vl7lHcvnxyKR3raTFqpqGqT3Ya1JCKe4+3d4MxybkNlo5490cs8hx+PQxuEZr0cOlTlVI475MPjqFXodVdSHYzh8Ew1VOXR4wBrp9/L5Phik1CV3oH0OvnbSgRuZCBKJCSBdo8fnc1w/UUTJKLh6gWg14Zj1YNQntccVbFc4dDplU2H25ZsjbEphx2wlajvP48Demw4Cu7zxjXVUFULwxGPhBxDTH42fVXDYTyO/sFvqAosTateluPqeSHjDuW40NscR6WiyBbKTXdkoZAwnowFZjjK9q56Z536K/S+0z+dcy+f9Cqt/s/3n+LW93+NV/75N/jq92fbNiCtEr9BzkAvlSvMruY35F6q+bIeeq8ZHx3aw4kI0bC0SI4XmEzFHYUJ/eQ49GPcQp+dwBgOj2iPw2uoCnRSs4ceR7ZAIhpaN98b6sbH9vCi6CVRNxlgaedcOk+potg5avUppIIKVRXcL0BehzkdX8gSEutCdPtH7uHHP/jttka96rLtxvMIrM1Q2d4U9JrZdJ6KYp1iL9QMR29DVe7zxjUiwthgrEVyPO9Yigu1znRfHocJVfUPi5kCA9Ewiah3a255HL3NcTR6G1A/PrZ3F0V9YRls8eWaDKCMUnPK1lvaZXscQWmLpXPu0txeDcecPd/9rl97Mb//iiv57vFFfv2fHvLtHVQ/OwfD0RcTCUfWX2ST8QjJWLinhsPLvPF6WsmOzKULjqW4YHnmyVjYm+HImVBV37G05r1rXDOc8Gc4/uW7pzY1KGcxU3AsFw6HhMFYuKdVVV6kNCaTsbbCC50Ik0wvWRehXbbHEWSOw+2L7rUqbi6dZypldXff9tz9/L8vu5w7D83wd9/yNwMta3uLTh5HkBMJz1RH2Q5suK/Xm5Bswd3g1zORijWvqkrnHUtxNalExNOmr1djY8EYDs9YciPe8xtgyY54Lcd9fGaVt/zjA/z93e0POnSSG9H0ega6FxE4/WX3YwiWs0We8btf4qvf39yMlell7XEEF6oqlCoUSpWmvS4arzmq2dX1IY+ffcGF3HD5Nt7xucM8etq7Ltha1eNwznFAb8vMNVXD4VDt1evu8XS+ef7OifFk3DE5rqdktjQcHr+7tRyHMRx9w/JagVEPA5zqGRmIks6XPDUm/eN3TgAwu9L+yb+U3Sg3oul1/b0X1deJZIx8qeLrYv3kfIbVXInHzqxsan2nl3IMxsJVSQ0dr+5lqMprSafXKYBWkrX2+YsI7371MxlLRvnlj3+3ahDc0KGqAYew7HDAEwljkZCjesNUj/WqrOS497D1RNI5x7GYLVCuKMdSXE3K4+jnTKHEYCzsaV7QZjGGwyN+BA41OsTgdmEslCrc8d1TwOYa4qwch/MavZ58nSJbFYFrkeNooxpGy30vZDZ34Tq9tMbOkUS1kqWWhOxd0tdLOA+8TQFUSjGbzldLUzXjyRi/+/IrOTqX4f7j3kas6uS4c44juImEupvdqfqo14PBMvmSo0fWjPFkjNV8iXxp/fk1V51Z09zjGPI4wXPVQ76sUxjD4ZGltWJVf8oresSsW0nulw/PsJApMDoYbdvdLpUrjgKH1bUk2hsf++Rcpi0pBy8icDqsMu/jC68TpIub1LiaXl6r5jegloTspcehDYdbOW4sEiIeCbXcgKzkShRKFcfeCy2H7tVL0J6Jc44juOR4s252sDyOpWxxw4XZjc89NM0L3/Ufvo/zOzBJi30uNmx4as1/HQhVtZD46TTGcHhAKUt1dsS3x6Hjwa2/sP/4nRPsHEnwsit3tG04luwkfKPciKadKYBPzWe48b1f5V8fOu17PV5201qKwc9OURuOzcpLnF7OVfMbmmTcWxKyU3gJ52mGEtGWeYVWFyD9GXjNS7SuqgowOb7irJ8FtSZAv+Xdn33wFCcW1nwVpSile5T8hapgo+xIVeBwqFWoytt5mc4Ve1JRBcZweCJXtJKYowN+k+PuQoenl9b4r8dnefW1F7B9OMFCttCWtMRiE7kRTTvJ8S8fPku5oji16H9UqA5VOV18NNWOXx9f9k54HPlSmdnVPDtH11+EUvEI6UIPDYePZKaVo2p+Hmm9o8ZQlT62/u+5sWa/B4nIxs8uGQsTkt4PKFNKOXaNa6aqYU/vm5ByRfHtowsAPDXv3XDkSxXKFeU7OQ4bu8dbCRxqUvGIN62qHo2NBWM4PLG05k+nSqNzHK2k1T9170mUglcf3FMTa2vjothMp0rTjlTEVw7P2M/t/yKR9lBVVZWb9vFl1zmOzczxmLGluZ08jiBCVd4NR/O1zXrwOLxuHLKFMgPRMCGHJKuItC07kiuW+a7HPEsjS9kihVKFbU08jsk2ZEcOT69Uv5t+53qDv+qlZrIjc+k80bC0DIMPJSKelJtXc8Zw9BXVrnHfOY7WoapKRfGp+07wgksm2TM+uCnNncWq7Huz5LjlcXiduLeSK3LPkwv2c/u/SFvJQ+eLjyYaDjE6GG0rVKVfbzuc1qW4o+sNRyoow+EhvODV43CqzomEQwxEvffxZIvllp5iuxV6f/21o/zYX36rLW9F/71mF9h2vjvffGIOgGhYfBkOfY74SY7XwrINhmM1z0TSWW5Ek4p7G8SWKfRm3jgYw+EJbTh85zhckuOnl9c4ubjGy67aAWxOrK2ZpLpGJ80yHkMxX/v+HKWKIh4JteUBZTw2SE2m4sytent+rdoqAkt2GWM76B6OxlBVso080GbwE6qaSLYuN51LFwiHpHmo0sfc8rVC2TExrml3IuHXj8xRrqi2zm+tyNvsvdIG088m5BtH5rlkW4qLp1Ic9xGq0mNc/ZTjjgxECYeEBYccx2SL/AZ4V35I57ZAclxE9ojIXSJyWEQeFZE327e/XUROicgD9r9bmxx/TEQeth9zb93t4yJyp4g8bv8/1q3XoFnWoSqfOY5ULNIyHqx3z3vt4fJTm9DccQ9V+auG+cr3ZhgdjHLdheNtehxlT7MKJlt01DYyl8lTLCsunLBmZ7Qr53Jad42PNHocYc+GtRN4CedpdowkmFlp3iw5l84zkYw19fD8TO7L2v0AzWhnCmCuWOb+40tAe1L6brMm4pEww4mI5+9OoVThO8cWeN7FE+wdH/TncfgYG6sJhYSxwahjjqNVfgPqGkBbGH6llCWDcq4bDqAEvFUpdQVwPfBGETlg3/c+pdTV9r/Pt3iOl9iPOVh329uAryilLgW+Yv/eVdoROATrZBlqITtyWusl2Ttf7XGcXc35XuNipsBgrLmWlt8RlP/52Cwvedo2JlPOHa9utBriVM9EKu65j0N3Dl+xaxhov7Lq9NIao4PRDbtqK8fR2z6OpEs4T7N9OEGhVKmei43MrraWrRjyUTGmcxzNGG4jVPXAiSUKJavow09OS+Mlr+CnCfChk0tkC2Wed/EE+yYsw+E1jOunGq6e8WRsg9G0DL43w9Hq88uXKhTL6twPVSmlppVS99s/rwKHgd0deOpXAB+1f/4o8MoOPGdLdALNbx8HtJYd0TvfnfbONxENM+Rj11RPK7kR8Ne49cCJJRYyBW64fBtjg7G2Kpi87n6mfDRu6ffrwE7LcLSj/AqWp9fobYC/cE4nSOe8x6R1GeqZFedNRSuFVfD32nLFzoeqvvXEfPXnuTbPJ2idD5oa8h72/OYT84jA9RdZHke+VPFsdPzMG69nPBlbtwlTSjGfLriHqjwUN1R7graAx1FFRPYD1wB32ze9SUQeEpGPtAg1KeBLInKfiLyh7vbtSqlpsIwTsK1b69YsrRWJhqWl+96MVkKHeue7bn5ym9IJbp3tNb0j9y/8f3xvhnBIeNFlU4wno2QKZXIuiblGvMiFg5U0XM1t7Kh1QldUHdAeR5vd46eX1qpeXj2pWKSqH9UL0j6E8rQibHPDUXBs/tMMxb1f7LOFcsvEbzsex7ePznOFbfDb8Tj032t1YZxMef/ufOPIHFfuGmZ0MMZeu0HSa0luLTnu73ow0eC9r+RKFMrOTZv1eBnk1UtlXOiB4RCRFPBp4C1KqRXgA8DFwNXANPCeJoc+Xyn1LOAWrDDXi3z+3TeIyL0icu/s7OYE8ZayRUYGYi0rH5oxnIg2Lce1ZC/W73zbFWtbdlHvHfIRqvrK4bMc3DfGyEC0Wn/uN8+RzbeuzNHoTvdmIZh6pldyxMIhLplKWWtqs5djejm34X2H3g9z8pPM1ONSZ5Y3Gg6llC1wuPkmMvCeHPcqTpkrlvnuiSVecMkEY4PR9nIcHj0OL9+dtUKZ7x5f4nkXTwK1HKPXPEctOe7vIj3RIK3upWscavN0vHgcvZjFAV02HCISxTIaH1NK3QGglJpRSpWVUhXgQ8B1TscqpU7b/58F/rnucTMistN+/p3A2SbHf1ApdVApdXBqampTr2N5reA7v6GxQlVNDMdyjt0NO99tw4n2qk5yxeoJ5oTX5PiZ5RzfO7PKDZdbjty4Xd7rN8/htRnJbaxm49p2jCSq8g3t5Dgy+RLLa8UNFVXQe4VcL/PGNduGbMPhIILpZefqtYkMbI+jRY5jKGGVh2Y8iibef3yRQqnCcy+eYCIV91wMUU86VyIkzsKLmqmhOOl8yVXM8b6nFimUrfUA7B4dICRwfD7jaS1+5o3XM56MsbxWrDb41kqovXkcrSSD2ukt2QzdrKoS4MPAYaXUe+tu31n3sFcBjzgcmxSRIf0zcHPd4z4L3G7/fDvwmc6vfj1L2aLvHg6NFapqluPonMexstZ6BKnXEZTH7C/PVbtHgFqVVqPGjhtek+O153c3AtNLluEYiIaJRUJteRy6FHf3qHOOA7yXLGvOruZ4xZ993fcsFT/DgGKREBPJmGOoqiZb0fwCNOyjj8e9qsqf7Mi3jy4QEji4f5yJNmew6I1IK69/0mP3+D3HFgiHhOv2jwPWe7tzZMCzx5H2MW+8Ht3Lob33ate4S46jqtzcwiD2UlIduutxPB+4DbihofT2XXaZ7UPAS4BfBRCRXSKiK6y2A18XkQeBe4DPKaX+3b7vncBNIvI4cJP9e1dpRxlXMzIQdfQ49M63sQltaihOplD2HS5ZzRWrX2gnkrEIIu7jY+erap3WyVz1CHzs7isVRabFvPF6xqtfJi+hqpqi7fhgrK0y4caChHraDVV958lFHjy5zHeOLfg6btVHchyscNWMg+GY9bBzTSUiKGU197mxViwz0MKg+S3t/vbRea7aPcJwImqPC24vx9Hq/Ib6qsTWz39meY3JVGzd+blvYpCnPIeqSp5KzRvRYV/9HfMaqopHrI1SyxyHj2bSTtC1v6KU+jrgtD1wLL+1Q1O32j8fBZ7Z5HHzwEs7tExPLK8Vq4k9vwwPRMkWyhTLlXU7lOogoYaQiT7559J5z65w2b5Qt/I4QiEhFWvdfaz/LtROZp2D8LO71x2uXr5ceraCm2GqVBQzy/nqBX8sGWsrOV5t/nPQPNINXX4Tv0/OpQE46VPTy2sBgWbHSKJaklyPN4VVO07uIktRLFtlnW59HODN48gVyzxwfImfef5+QE/Ca8fjKLq+V171qhYyhQ0lsHvHB/myLbHjRtbnECdNY1h2Lp0nJM17r+oZikdI55u/372c/gcePA4ReXVd2Oh/icgdIvKs7i+tf1huQ1JdM9xkd3aqYXSppp3ucbfmKI2XBGnjyaxDdH5yHH5iwKMeQ1UL2QKFcqV6wR9PRtvyOE4tWZ3nTvLcNY/DXwXZ0VkrvOdHDFIp5XoRb6SZx9FKbkTj9WLfShm39ly2IoIHA3u/nU+4/iIrnzCRtOTP/Qp5pj3kg7Z5/O7MpQtVj1qzd2KQuXTBs3x5Oxdo/Tfn6wzHeDLuafCS23fX6zWgU3gJVf22UmpVRF4A/CBW78QHurus/qFYtibUtZ8cd1bIna42/23McYA/w6FDYcMurrwXjaG5dIHxZKx6MkfCIUYG/F2k/STqYpEQQ/GIq2HSM8L1Bb/d/pLppTW2DcUd49OpNkNVT8xZhuPkkj+F1VLFX8PW9uE485nChnLh2XS+pdwI1JV0ury2VrM4NM02Q05859giInBwv1V1X5tL4bPYwoORHU/GEHH/7lgeR4PhsCurvOSpvMrpOK0PYCGdZ61Q5vD0aktjX4+bunU6XyQSEuKR3qhIefkrevv1Q8AHlFKfAfxpb5zD6FLatg1HE4Xc00trhAS2NyQ0a1Lj3g2H/gJrUcVmeJFWn0tv7EBubFxyQ+9avX65xpIx12a+xhDT2GCsraqqZqW40F5VlVKKJ2etUJUfj6OdKhjdBNioLDC3WmgpNwK1/gc3jzNrFwa07Bwf8J4cn1nNMZGMVb2UmqaUv89u1YPHEQmHGB+MuYaq5u2dfj37xr33cvidN64ZG7QM2zeemOeH3v81HjixxI9de4GnY9101DL5MqlE6+KBTuLFcJwSkb8CXgN8XkTiHo/bElQFDtsNVWmPo+FLdmopx/bhBJGGne94MkbIw66pHv0FdkseehkfO5/Ob3Djxwbb8zi8JhDHkjEWXJLjupqoPsexvFb0LXTYaqZDO8nx+UyBFXs3fGppzbtsRRtVMNtHdElug+FwMPaNeFUOaDU2tvZctuqzy2RLsDyL0TpPSI9I9VuS67XnZdKlKjFXLJMplDeGqvx4HG0mx7VXeOehGfKlCh9//XP4uRde5OnYIZdNXy8l1cGbAXgN8EXgZUqpJWAc+PVuLqqfqAocekhgOVGVVm/4kk0vrzlewMIhYcJnSe6qx/immzQ3OIuujftMRPutcx8bjLqGLk4v5YiGpRpiGB+MolTrWSdOrOaa56ui4RAxlxGtjej8xnMvnqBYVq4VPZrNeBxnljcqrLYqxYV6rbLW71ctVNV8XQPRMOGQePI4FjIFxusNR3UGiz+PI+OjL6jVJkfnFxpDVSODUUYGojy14N7Lkd2EmOCPXLOb112/l39/ywurDYhecJOM8VI80ElcDYdSKms37y2LyF4gCnyv6yvrE9qdxaEZaeJxWLIXziETv70cWnLazePwInQ37yC65jef4FcEbnzQPRR2ZnmN7cOJajhmzEfjYD1WR3vrfhd/hsMKU73oUusicMpjnmNThqPB43ATOKz/O24eh5fkuDXMyZvsyGK2sG5GzITHyqd6dNWgl3yQW1hVlwJPOLxflkque7hxM5P2/tcPH+APXvl01+9qIymX724vp/+Bt6qql9s9E08CX7X//0K3F9YvtKuMq9E5jvrkuFLKmnndzHD41Kvy43G0uiiuFSw3vrEhaTxp5RO8SkzUchzeQ1XuOY71woTjDc1UXlBK2YnN5uvyO8zpybkMsUiI6y60qoa8luS2oy00OhglFglxts5wKKUsL9GjUJ5Xw9EqxwHevFew+nPqZ8QMJyJEw+KrJFc3ZHr1OFoajkzzuTV7JwZdu8etc8jfvPFOkEq07vz3I5jZCbyEqn4fSxb9+0qpC4EbgW90dVV9xGaUccHauYVDss7j0JUxu5rNT/aouaPxajhS8VpPiRPVfoBGjyMZo1CqVC8qbvgNVY0nY65CimdWcutKaHUFkR+PI1esUFGtJ7cl4xHSPspxn5jNsH9ikD3jllHzbDja8DhEhO3D8XUeh1ehvHBISMbCrt7UWtGbgJ8lmtj6uZRSG3IcIsJE0l8ToJ9S07FkjKUWua/56oxvB8MxPsjJxbWWebN25o13gqF4awHO1X7zOICi3XQXEpGQUuouLIHC84KltSIi7mGgZogIw4nIuhzH6SaluJqpIUtq3M+Y11gkRDzS+ss+nmp9sa0ajkaPw+dF2s+AIqh5c82EDpVSdjVUneGoiiP637m29jjCrnmAeo7OpbloMsVgLMJEMtZVwwFWuKq+CVBvMNxyHOCtj6cWqvKSL2v9XKv5EqWKWpfjALsJ0EeOw4+A30QyhlLNzws9ga9ZqKpUUdUKvlZr8Xpudwq3UvFMvrXkUKfxYjiWbIXb/wI+JiJ/CrQ/8PkcYzlbYDgR9dSk04zhBtmR002a/zRTqTjFsvKc+F3Nlaq19a1w6xGpyo04eBzgPSyUyZeqCVQvuBmmBdtDq/c4asd4PxWzed3R3trj8NoAWCxXOD6f5aIpq5Rz99gAp5Z8Gg6fX/bGJkCvshVgq9p6To633oQ0ntNO6LzYWLLRcMR9zeRY9RHWcztX59MFYpGQY1WUziE5CUlqqudQjz2OVKK1Qq7fZtLN4sVwPAhksTSl/h14gvMpOe4iV+6FRml1Lx4HeO/lWFlrrVNVe97Wc5lrHsfGqirw7nH4jQG7eQ96xG69xzEQC5OIhtorE+5QjuPk4hqliuIiW+b9grEBTi56TI57UHt1YsdwgjMruWq+yY/hSHmYqe4lOQ7ePA6tPzaeXH9uTiZj/kJVPrwzt6qtuXSByaTziISaakPzCZy6d6jd0HW7tMpRlX1ow3UKL4bjJUqpilKqpJT6qFLq/cCzu72wfmEzyria4YHIuuT49PIa8UioqtPUiF/ZEUsAzluNOzRvvmpWqug3Ee1VGbfx+Zs19OmGt23D63NCYx6qserRzW2dqqrSFVUXTtoex+gApxbXPBUR6AmJfhu2tg8nyBUrVbkPP6Eqt+IIsAxHJCSuyq/DCe8eR2Mpu9MI1Vb4ynEMtj5XFzL5asi2kW3D7iKJtc980HUtnaTVPB0/xQOdounZISK/KCIPA5fb0/r0vyeBh3q2woBZWisy0mYPh8ZSyK3PceTYPTrQ9KLhVXNHYynjejcczZ53djXPUDyyYW6537BQxodcOLhLq5+1QwfbHQyHvxyHe7VX0pfhsCpwLrZDVReMWSNIvXRFp/PehzjV09gEeO9TiwwlIp42N168BLexsfXP5SbTro36xhxHnLViuWrI3dA5J08eR4MeVCPzDgKH1WOTcUJSO9+ceHIuQzgk7LU7zXuF3og55d96rVMFrdVxP45VdvtHwNvqbl9VSvnTjz6HWc4W2De+ud1F4/jYUy16OKA9j0MP+mlFMh5hMBZuGqqazxQc51YPJSKEQ+K5l8NvTbkOBTYzTDrm3Fg55FsKxUO1VzIeIVsoU6molhIeAEfnMownY9Ud9QVjurIq6+oBpHPtNZHVmgBzDETDfOHhaV7/wotc1wruvQDgPotDM2TLtGcKzeXO9a5/Y46jFk4aHHd/D/zkOPS51OxcnU8XqhMkGwmHhKmh+AZJl3qOzmbYMzZArEeaUJpWoapeT/+DFh6HUmpZKXVMKfUTSqmn6v6dN0YDNqeMq2lMJDbrGtek4hES0VDLE7ger6EqsLyOpjmO1fyGMBVYkuxjg1HP2lDZQplBHzmOaDjEUCLSNLwws5pjPBnb8GUdS8Y8zfHQeKmI0V6Al2FOR2fTXDRZ23nutg2HlwR5ptBe3X19E+BHvvEkIRF+2pYsdyMVj3oKVblVVIE3CZPFbIFwSDYUbtT0qvx12XvxYuORMEPxiKPHoZRiPrNRUqeebUOJlqGqJ2bT1ZxWL2kVqvJjWDvFeaM51Q6VinKd5e2F4USEXLFCvlSmUKpwdjXf0uOoDSryWlVVrGpiuTGZai4CN59p3oE8NhhjwWNcOu0zxwGtpSLOruSr4bt1xwxG2xJfbLWj9iOtfnQuU81vQG2qoJeS3Ha1hXQc/vGZVf7xOyd4+TN3NRVtbESHl1r1KawVyp4S9l6GOS1kiowNRjeEZCcaBhq5kc6Vqv1QXhhLOisdZAtlcsXKBoHDerYNxZuGqioVxbH5zLrNQq9oNcFTb0r7IsdhsOrQK2rzFRT6on52Jc9fffUJlHIeXVrPUCLqqTPXyxCnelqJwDnNKdCMJb2r0WbyJVI+69xbJbrPruY2JMbBSrourxUpeZztUOvjaBWqsi6abjvz1VyR2dX8ut3nUMLSO/KiktuuREQiGmZ0MMrff/s42ULZs0ietT53b8ryOLwYDneF3MVMwVHqvZaH8O5x+HmvxpPOw6L0+dXS4xiON/U4Ti+vkStWuHCq94ZjMBZGxLmP4ylb1l83ofYCYzhasFyVG9lcclzLjtz6/q/xnju/z4sum+KWp+9oecxQQ9NgM2qJMW/GzWou3PilKpUrLGY3Chxqxn3oVbUzIc3N42iUn9fHgHehw2y+TNhlZoFXafUn7S/rRQ0XEa8luV5F+5zYMZxgrVjmBZdMcmCX98mUXubOZz0mx73M5FjIFjbkN6DmcXiVVvciqV5Ps3OpVrrc/Ps8NZRgPpN33IzoYoiLJnsfqhIRq5za4bw8MptmOBFxVQ/oJMZwtGBJK+Nu0uPYYyfXr9gxzCd//rn83f9znbsgYSLi2qwFNTfVj8exmC1s+GJYWlTNv1RjLqqjGi96UE6MDkZZdEiOlyuK2XS+GqJpXBP4KBO2E7+tSmC9DnM6YYvh7W0onNg96q0JsN3kONSqy37uhRf6Os5LXmLNc3LcWbyznqVsYUNFFVg9OMlY2Feoyk8F2njSOaxarfJyCVUp5VyVpUtxLw7A44DmIqVHzqa5ZFuqZ7M4oIszx7cCmxU41Fy7b4x7fvOlTA3FPX+4wwNRjs65SzxXhzh5NRz2F2MhU1gX/plb1Ro+TTyOZJTFbNG12mitWKai/HfWWjkd5y97uaI2lOLqY6zHdK5MOOnR49BhlsbqqQvGBvn6kTmUUk0/60pFkS546/Z34jkXjVMsV/iBy6Z8HedFWt1rclyvvdX42IVMkWv3OX93JlJxf6Eqnx6HU1i1pozQKjlu93Ks5Decc0/OZUjFI556ZrpBM2n1I2cz3HC5v3Nhs3TN4xCRPSJyl4gcFpFHReTN9u1vF5FTIvKA/e9Wr8d6Pb5TLG1y+l8924YTvnYEXmWrvQ5x0kzZHkVjHHe+hYYPWDmIckW5rsmvpHr1+ZMxO3m5Pimt+xWckuNarttPR7tbtZdXj2Nudf1sds0FYwNkC+WmulsA6UIJpdrXP/ulF1/Cx19/ve8dpheF3LWC1z6O1jkOpZQlqd4kzOtHr8qvnMZ4MkauWNnQJzLvKcfhPGURrGKIi6aSPd3Z1+PUnLqcLTKXznPJtt6Gz7rpcZSAtyql7heRIeA+EbnTvu99Sql3+z1WKXXI4/EdQcfOvVYsdZIhu/ej1c4VvCvjavRuqbGyyi3+W9/dPdLCkNb0oPyFquq70+urhHQi3yk57tYlvHFt7hcgrx7HXGb9bHbN7rFaZZVTfB/8f2adwkteIuuxqioRDREJSdPnWslZ1VtO8uVg5Tk8y7PkS756FOq1z+q9p/l0noFouKVHVfU4HBLkR2czPNuenR4EqYZ+MIAjs6sAPTccXfM4lFLTSqn77Z9XgcPA7m4f20mWA9KlAeuiUqoocsXWFUNehzhpmsmOVOWmm7jhXgcnte1xDDp7D9rjcApV+TUcGQ8VQ16T41bPy8b3qr4JsBl+vcROkWrRCwCWl7BW9FZVVRvm5Oxx6I7+Zh7HZMq58skJr8oImmbn6rxt7Fuhvx+NJblrhTKnlta4MIDEuMZpfOyRs1be5ZKpoZ6upSfJcRHZD1wD3G3f9CZbvuQjItLShDsc6+v4zbCULTIYC7vKlXeDYQ/ljtb9/navk00msM2m88TCoaZJyHEXWRBNptrF6r8cFzZKqzfrGgdr1xsNi6fqM702txxHIhoiHBLXUJXVZb/xIqS9pcYpffUE5XHUQlXO55TepHgJVYEuGXd+n/RFeyzZLMdhlV+7jQ5QSrVVjlu/Bs18ptCyogogFgkxnoxtCFU1q6LrJcl4eENy/MjZNPFIqOrp9oquGw5bkv3TwFuUUivAB4CLsWZ6TAPv8XEsXo8XkTeIyL0icu/s7Gxba98xkvA1F7iTDHlIPkJtsqDXi1AyHmEgGt7QyzGftr5UzcJibkKEGi+9Ei2fv+HLfrZJ1zjoXa+3fhfQHe2t1yViDTxyawB0GrELlucUDUtLae52pv91gmQsgkjzctyqCKRHxd6xFg2Yiy4ex0QyTrniPjpAF1v4TY6Dg+FI5109DrCbABu+H0fnrJ19kIbDqfP/yFmrk30zYx/aoauGQ0SiWBf+j9lzy1FKzSilykqpCvAh4Dqvx/o5Xin1QaXUQaXUwamp9ioOfu6FF/HXtx9s69jNMuyh3BGs3auXIU716EFR9cyl800T4+BdWj3TZo6jWWntTJOucY3XIgJrbSVSHsqEvSjkNmuWFBFLtqKFx6E/03arqtolFBJSseYjSL0OcdJMDSWaNpPq0upmF2qvG5F2Bl41O1cXMoWW57hmyslw2D0cFwbQNa7RVVX1XtqR2XTP8xvQ3aoqAT4MHFZKvbfu9p11D3sV8IjXY70evxXwIukAlkfi9wLkJDuiPY5mWCG7kOdQlV+PQ/fKNH7ZZ5t0jWu8zr4G76WmSRcxwFyxTDpfalq6vH04zkwLnbFVn02bnSTVwtCuFb0NcdK0GnHcTOBQk/RYvdaO8utwIkIkJOvOJaUU8y2UEerZNpRgdmVjqGrXSMKzUe0GOoyctT+nXLHMycW1pqKN3aSbHsfzgduAGxpKZ98lIg+LyEPAS7AGRCEiu0Tk8y7H0uz4rYYXSQd9/7DPC9BkKl7t29C4eRzWrOiYa7dvu8nxSDjEcCLimONo6XF4mH0NdY2JHi6KqUSkpSyHWwWaNaWveagqqByH/pvNQ1Xehjhptg3Fmc8UHGfYL2QKRELSNGfmVdql3dnsjQ2rq3lrNnurHg7NtuE4s+n8urkqRwMSN6xH59QeO2NF7Z+YTaNU7yuqoIvluEqprwNOgbfPO9yGUuo0cKvLsSilbuvUGvuZ4QFvHocfZVzN5FCc+55arP6ud2NuU+S8NG3pi4/fUBVslEmv2F3j2x26xjVDiQjHF9zLOteKZZTCNccB7qGqagVaU48jwdePzDU9fjVXJBwS39P/OkGr16ZzHF49Dt3NP58urBvrC5bHMTrYPGemixTccknVfFAbDaX159JCk7HITmwbskY3L2aLjCdjKKU4Opvhldf0vLBzHTcf2MFw4lH++mtPcu2+8VpF1VYKVRk2R1XSwSV5aJUq+vM4plJxFupkR1Zy1m7MreLES9NWJl8iHgkRcZkg50TjLnHe7hpvNWukVWXP+nV5N2jJWOvxsdrjaOahbR9OsJorNR1UpI19EI1kqRbFBGt+cxy6dNUhLLeYKW4YGVuP9kjdhjnpfIzfQoLGTYje8DSb/lePPt/065pN51nNlwJNjIP1nr3u+n38+6NnODaX4YmzaUIC+3s8jRCM4ehbkrEwIemex6FlR4Dq/GdXjyMZ99TH0a54X+MusdbD0drjcCsggNoFyksILdUinAP1HkezUJVzL4DGb19CJ7E00FqHqrx6Qjr35JTnWGjRNQ7e+2WqOQ6fQ4o2GA79mXnxOBo+v289MQ/AM/eM+lpDN/jp5+0nGgrx4a8/yZHZNPsmkoG0CxjD0adU1TA9VFX5vQhp2ZFZ22DoChK3xOGEnVRvNVPb77zxekYHY+tyHK26xjXDHkaYQu0C5WU3PZ6MMZcpNH2dc1qepclFSDcrzjSprFrNlXxfCDtFM6E8qCXHveY4plp0WTeTVNfoHIdrcrxjHoctcOjJ41j/ur50aIapoThXXzDqaw3dYNtwgldes4tP3XeCB44vcXFAeRdjOPqY4QH3MEw7oarG7vFvHpkjJHBgZ2uJ7olkjHypUp3d7USm4F9SXTOejDp6HK3LcaPVEaatyHqYN67ZNhSnUKo0bSycWy2QjIWb5gK0xzHTpOJoNe/f2HeKVuXLOlTluaqqSZc1WDmOZhVVUMtxpN1yHNViC//l3UtrxerQquosDk99HLVQVb5U5quPzXLjFds8jeftBa9/4UXkihVOL+cCyW+AMRx9zVAi2rIB0O8QJ03jTPMvPHKG6y4cd61x1/fPtxj5aXVnt+c6jyVjrBXL1QuY3vG1UiP1Wrac8eFxaA+nWUntfCbfVJql/vhmvRzteImdIhWPslYsO86b8FtVFYuEGBuMMpte/zotgcPWOY5QSBiMhV09jnb6lMAyEErVpE/m0nmSsTAJD2G4gZg1fvbsSp5vH10gnS9x4xXbff39bnLp9iFuuHwbEJzEuzEcfYxb/N7vECdNvezIkbNpHj+b5mVXth4sBbVQVquSXL8S2PXsn7C+BHc/acWUZ1ZyjA1GW140vMyYgNpF0Uv+ZftQ6xzFXNp5Nnt1TXZ3fvNQlX8vsVPoz8apmmnN9toSPptJG98nLXDYKlQFlhF3S46n80Vfszg0jQ2l9x5bZL+P5r2p4ThnV3PceegMA9Ewz78kGAWJZrzxJZcwFI9w7b5gRBeN4ehjhl26ov0OcdJo2ZG51TxffPQMAC+7aqfLUTU3v1WCfClbbFsU8sYrtjM1FOdvvnEMsDwOJ3HDemoeR+tcUM3j8BCqaiGtDbiWLouI1QTYNDkebKgKnBUJtDKun5DMtqHEhhyHbhJ1MxypeNg9VJVrbyOiz9X5dIEHTyzx8KllfvzZezwfv23I+vy+fOgsL7ps0pOn0kuu3TfGQ2+/ObDeEmM4+phhFx0mv0Oc6pkcshLdX3hkmmv2jm6ow3fCS6hqea3Y9sTEWCTEbdfv46vfn+XI2TRnV3KuQ3P8hqq85F90TqXZhd+tWRIs4+PkcWjRvsAMR4tqpqxHZdx6tjl0j2sZETddqGS8ddmzXmc7VXr1ysl//+2nGIyFeZWPPoxtQwkeOrnEmZUcNx1w98aDIKi5IGAMR1/jpsO0sgl57qlUnAdOLPHIqRVPYSqo28U18TgqFcVKrn2PA+Ann7OXWDjE337zSY8ehzdNr4yP+H0yHiEVjzh6HOWKYsGDyur24Y07cbB29eWKCjxU5WQ4vA5xqkfLjtRXoC25yI1okh40wVZ9DnHS6LDq0bkM//rQaV5x9W5f77luAgwJ1XyCoYYxHH2MVn5tVha6GemKyVScY/NWx/XLrvJmOBLRMKl4ZIPOVXU9eWuy3WYGX02m4rz86l18+r5TtuFovbP3MpwIrD6OSEiIO6jsOrHNIXYP1kWxotx7XrYPxZlZyW347IKUG7H+bnMpm6zHeeP1TA3FKZQr61Ru9Shfp3nj9aQ8ehztvFd6auffffMpcsUKP/Wcvb6O170cB/ePe1LUPd8whqOPGUpEqCialr+uVlVW/V+odVXQFTuH2TfhPWnYWB9fz7Ldg7HZwVc/8/z9rBXLrl3j4D05nslbYRiv7v02OznaiC4McOt52T6cIFsob9hRBzXESTNSFZNskuPwKeLXWKEHtRzHaIuqKsBTVVW7oap4xKqMOrOS4+o9o1y1e8TX8fq8u/lA/1RT9RPGcPQxeufeLM+xWY8D4BaP3oamlezIcnVG++Z2aFfuGuE5F44DrbvGoX6EqXty3E9/iVPSF7x32esda2OeRHdtB+Vx7LRzWdNLaxvuWyuUGYj6uyTUeh5qr3Mh21rgUJOKR1r2BEH7yXGohcped/0+38deu2+Mq/eM8sPP2NXW397qGMPRx7glfjeze71g1JoY5ttwJDfO8tBow9GJUbs//wMXIYJr1UhthKl7Oa6fMMy2JqGmWRdlXM32Jr0cVWPfZpPkZklEw0ymYpxyMhxFb7Lz9VTlOeq8s7lVa2CSm3fnJTm+ugklgvFkjJGBKD/8DPeKwUb2jA/yL298vqeikfOR4MTlDa64CR2u5ixBQafpeG684ppdXLFzmEu3+5tVPJmK8eDJJcf7Omk4brh8O/f9r5s8xZe9TAHMFPyFPLYPJ8gVK6zmS+tCgW7KuPXHw8YmwqBDVQC7RwecDUebyXFYH6p6+NQyl7uoEIBlOLKFMpWKciwBzpfKFEqVto3sW268lFJZ9V0p7VbAGI4+xs3jWMmV2r4AxSNhnn6Bv7gvWKGqRXtWdOOXfWnNuqh2wnCAezmnxpPHkfe3m64J3eXWGY65dJ5ISFzzSs1KeoNOjgPsGh3g+zOrG27PFsqex8ZqhuIREtFQtZAgky/x/ZlVbvZQqaenMWYKzuexblJsVzTzxU8z1VDdwoSq+hi38bHWEKfeXoDGk3FKdtltI7UcR293014MRzpf8qV3VI3dN1z459MFxpMx1ya5ZDzCUDyyoZdjtc2mzU6yy/Y4GsNw7VRVVUfl2h7HQyeXqSi4Zu+o67G1KYDOeY7abPbgvDODM8Zw9DFupabLa8VNlb62w2QL2ZHltSKxSKjnoQFL06t1qMq6KLbhcTQkyOczedcwVf1zbDQcJURqIn9BsHt0gFyxwmLDtMW1ov+qKlg/QvaBE0sAnpRkq8OcmsiOtKuMYOg+xnD0MW7NbScWslwwNtDLJVWlxJ26x1fWNtf81y5ePA6/qr21UNP6C/+sx7nV4DxCVje0Bam0ussujDhdl+colisUy8q3xwF2z4udy/nu8UX2Twy6Nv+B+9zxfgjrGZwxhqOPqZWabvxilcoVTi6usXe8t9O/9EXTqZdjMzpVm8FNmgUg61O1NxWPMBgLb/Q40vmqnLgb2x1kR1Zzpbb6bjqJ3mycXKwZDh1mbOciPTUU56zdPf7AiSWu2etNeM9t7ri+PajZJYbmGMPRx9RKTTdeFKeXc5Qqin0TPTYc9k5yzsFwLAfocaTzpaYd9pWKsnSYfHgcVux+fahJKWXrVHnzOLYNW93n9esKcvqfxsnjeMKeX92OaN62oTiruRLH5rOcXc1ztcdJeSmXHEc/5IMMznTNcIjIHhG5S0QOi8ijIvJm+/a3i8gpEXnA/ndrk+NfJiKPicgREXlb3e3jInKniDxu/x+MrnCPaDbM6fiCJReyp8cex1hVdXRjqGozAoebwa3Dfq1YRilv88br2dagN5UtlMkVK55zHNuHEhTKlXVTDYNUxtWMDUZJREPrSnIftw3HpW0MBtKFBF+ylZa9JMbBhKrOZbrpcZSAtyqlrgCuB94oIgfs+96nlLra/vf5xgNFJAz8OXALcAD4ibpj3wZ8RSl1KfAV+/cty1Ai4tjH8ZStM+VHLqQTRMMhRgejjt3jwXkcrTvsdfLVj8cBG5VfdeOjmzKuxqmXYzVfbLu8tFOICLtHB9Z5HEfOpknGwtXOcj/oXo4vHZohFglx+Q73Hg6onwLYOlTVbue4oXt0zXAopaaVUvfbP68ChwGvusbXAUeUUkeVUgXgE8Ar7PteAXzU/vmjwCs7tug+ZCje3OOIhUPscFGP7QYTTfSqlrO9r/IC936XbLUfwK9k+PochVedKs12B9mR1U303nSSXQ6G45JtqbakurXhuP/4IlftGvbckKpzHM2GOa3kim1N/zN0n57kOERkP3ANcLd905tE5CER+UiTUNNu4ETd7yepGZ3tSqlpsIwTsKW7fJpVDB1fyHDB2ADhAKpznGRHyhXFar7U3x6Hz1LT7cPxdUKFOjznJzkO6yuz+iFUBRu7xx8/u8ol2/ypCGh06bJSeE6Mg/vccauQIPj3yrCRrhsOEUkBnwbeopRaAT4AXAxcDUwD73E6zOE258xn87/7BhG5V0TunZ2d9bfoPsLKcTiHqvb2ODGumUjFNszkWAmo+Q/qp9o571x18tVv70RNqNC68Os8gNuMEM324QSRkPDUfAawkutBjo2tZ/foAHPpArlimeW1IjMreS7d3t40uYlkHL1/8ZoYB/e54+k2Z3EYuk9XDYeIRLGMxseUUncAKKVmlFJlpVQF+BBWWKqRk0D9nMcLgNP2zzMistN+/p3AWae/rZT6oFLqoFLq4NTUVGdeUABYc8fXf7GUUhyfz7Kvx4lxjaWQu97jWOqgTpVf3BolazkOfyGP7XXd40op/uW7p7h235jrVEJNLBLikm0pDk9b8h75ktUr0Q8eh66sml7OccQ2iJe0OYY0HJJq3seP4YDWQof9YmQNG+lmVZUAHwYOK6XeW3d7vVTlq4BHHA7/DnCpiFwoIjHgtcBn7fs+C9xu/3w78JlOr72fGEpESedLlCv1E9aKrOZLPa+o0kwk4yxmi5TKleptnRQ49ItbqCrbpuZRvfLrw6eWefxsmh991gW+nuPAzmEOnV6x19f+qN9Oow3HqcW1ailuux4HWOG7yVTcd0NqqsUUwH4J6xk20s1P5fnAbcDDIvKAfdtvYlVIXY0VejoG/DyAiOwC/lopdatSqiQibwK+CISBjyilHrWf453AJ0XkZ4HjwKu7+BoCR19k0nX5g6cWgqmo0ujk8GK2WN19B2s4PHocvpVfax7H/U8tEouE+CGfEt0Hdg1zx3dPMZ/O94UyrkZf4E8vrfH42VXikRAXjLW/EXnF1bsoliu+k+uDsTDZJmXU6XyJvclgNkeG1nTNcCilvo5zrmJD+a39+NPArXW/f97psUqpeeClHVpm3zNct5vWF+XjVcMRnMcBlm5To+EIIscxEA0TbjHMSYdC/OY4hhORar/DZx88zc0Htvs2jFfY8uKHp1erBq4fdtHbhxOIwMmlNR4/m+biqdSmCi1+/gcubuu4VnPH+6UCzbAR0zne51QTv2u1L9dxO9m6ZxM7xM2gPY76Xo7lrPVzEOW4bsOc9I7Wb45DK79+5oFTLGaL/Oi1/sJUUDMch6aX6xragr8YxiIhtg8lOL20Vi3FDYJWc8dX+qDL3uCMMRx9jlP8/qn5LNuG4r6H7nSKmkJuLUEeZKgKWgsdZvIlomFpqx9g+3C8GpJ74SWTvo8fT8bYOZLg0OmV6mfYL5VCu0YTHDmb5uTiWlsd452gWXJcKUU6b3Ic/YoxHH2OU/z++EK25+KG9ehQVX0T4PJakUQ0uGYtq1GySXK84H8kqkbLabzy6l1Ewu19XQ7sHObQ9ErfSWjsGh3gIXua42YS45shFQ879nFkCpZMTL+8V4b1GMPR5+jQz2q+dlE8vhBcDwdYXkU4JBs8jtEBbx3V3cCpbFmT8amMW4+urGonTKU5sGuYJ2Yz1XnlQavjanaPDaCL9dpt/tssg7GIY+d4PxUSGDZizHmf0+hx5Iplzqzk2DceTEUVWI1bO0cSHF+odR4HJamuGUpEHedog1VV5VenSvOag3vYNTLgWX/JiQM7hylXFPc/tQj0j/bSbrskNxKSwAotms0dr07/65OwnmE95lPpc2rJcWsHdnJxDaVg70RvBzg1cvmOYQ5Pr1R/D0rgUDOciPC9plVV/oY41XPFzuFqgrtdDuyyjr/n2ALJWDgQmRgndo1Y59CFk0mibYbhNkuzueMrfRbWM6zHhKr6nHgkTCwSqnocxxesiqq9AXocAAd2DnF0Nk2uaMWngxhjW0/rqqr2Q1WdYM/YIKl4pO/KS3fbvRxB5Teg+dxxE6rqb4zhOAcYTkSrO7CanHqwjVFX7BymouDxGavreGWtGEgPh0Z32DsNc8rk20+Od4JQSLhip5VD6Kcd9O6xAUICl20PJr8BtVBUYy9HdfpfH71fhhrmUzkHGB6I8KVHz1CpKI4vZBmMhauT+ILi8mpj2wpPv2CEpYBDVUOJCOWKIuswW3x5rcjlO4K7OIJlaL9zbLGvLoTDiSj/92efw1W7RgJbg27KbCzJ7bcKNMN6jMdxDvDrNz+NZ1wwwr8/eoZvHZ3n0u1Dbc1N6CT7xgcZjIU5NL1CsVwhWygHnhwHZ9mRhUyB8YAN7QHb0Kb6LPTy/EsmGQnQUxysy3HUY0JV/Y0x5+cAtzx9J7c8fSeViuLoXDrQXIImFBKetmOI751ZCbz5D2o703S+CNRkz9cKZdaKZcY9Dl/qFjpBbnbQ62k2dzydKyECg1EzxKkfMR7HOUQoJFyybajalBY0V+wc5vD0anWmdrA5DueZHIu2FMr4YLCG47LtQ4RD0hfKuP1Es7njK/YsjlCfVKAZ1mMMh6FtrtgxxPJakcfOWPMmgq2qcg5V6e72sYBDVYlomN/54QO85uAe9wefRzRLjlvT/4L3rA3OmO2PoW10f8O3j84DwYaqasOc1vdyaMMRdDEBwO3P2x/0EvqOZh5HOl80zX99jPE4DG1zeR8ZjlSTmRw6VBW0x2FwRucwMoXGPg4jcNjPGMNhaJtUPMLe8cHqLO7RvghVrfc4tPR7P3gcho00mztuDEd/YwyHYVPoxjYINseRjIUJibPHEZL+ERY0bMRJWn01V+y70mVDDWM4DJtCi/8lY+HA9I7AGro0Ohhjvk7qHawcx9hgzFTn9DFOc8fNLI7+xhgOw6bQCfIg8xuaHcMJZpZz627rh+Y/Q2uS8Y2hqhUTquprjOEwbArdEd0PTYm7RhOcdjAcJjHe3yRjkXUNgPlSmUKpYsKLfUzXDIeI7BGRu0TksIg8KiJvbrj/10REiciGeZwi8jQReaDu34qIvMW+7+0icqruvlu79RoM7lwwNkAqHgm0+U+zYyTBmeX1MzkWs4XAm/8MrUnGI+skR8wsjv6nm59MCXirUup+ERkC7hORO5VSh0RkD3ATcNzpQKXUY8DVACISBk4B/1z3kPcppd7dxbUbPBIKCbc+fQc7RoKdDwKwc2SAxWyRXLFMwi7zXMgUOLjfGI5+pjE5bgQO+5+ufTJKqWlg2v55VUQOA7uBQ8D7gN8APuPhqV4KPKGUeqpbazVsjnf92DODXgJg5TgAppdzXDiZpFJRLGaLxuPocxrnjtcMR/BerMGZnuQ4RGQ/cA1wt4i8HDillHrQ4+GvBf6h4bY3ichDIvIRERlr8jffICL3isi9s7Ozba/dcO6wc1QbDitctZorUa4ok+Poc6wcR53Hkbd6cUyoqn/puuEQkRTwaeAtWOGr3wJ+x+OxMeDlwKfqbv4AcDFWKGsaeI/TsUqpDyqlDiqlDk5NTbW7fMM5xE47XHbGTpDPZ/KAaf7rd5LxCGvFMuWKNYTLhKr6n64aDhGJYhmNjyml7sC64F8IPCgix4ALgPtFZEeTp7gFuF8pNaNvUErNKKXKSqkK8CHgum6+BsO5Q32oCozcyLlC0p7JkbUT5NpwmKqq/qVrJl2sSUMfBg4rpd4LoJR6GNhW95hjwEGl1FyTp/kJGsJUIrLTzp8AvAp4pMNLN5yjDMTCjA1Gq6GqhYwV8jA5jv6mfu74UCJK2paNSRmPo2/ppsfxfOA24AYvpbMisktEPl/3+yBW5dUdDQ99l4g8LCIPAS8BfrULazeco+wYGWB6yfI4FuxQVdBDnAytaZRWN6Gq/qebVVVfB1rqPCil9tf9fBq4te73LDDhcMxtnVulYauxcyRRDVUZj+PcQM8d1wKVq/kSiWgoUAkbQ2vMJ2PYUuwcSXBmpZbjSERDDMTM+NF+5nJbKPM7xxYAy+NIxU1+o58xhsOwpdg5kmAhUyBXLDOfLjCRjAe9JIMLF4wNcsXOYe48ZNXArOaKZsRun2MMh2FLsaOuJHcxW2AsaXau5wI3HdjOfU8tMp/Om1kc5wDGcBi2FLtGaiW5WlLd0P/cfGA7FQX/8b2zpPMlU1HV5xjDYdhS7BipdY8vZAqm+e8c4cpdw+wcSXDnoRlWc0WGTI6jrzGGw7Cl0N3j08s5Fo2k+jmDiHDjFdv52uNzzKcLJlTV5xjDYdhSDMTCjA5GOT6fZTVfMqW45xA3HdjOWrHMfKZgQlV9jjEchi3HjuEEh6ZXANP8dy5x/UUT1WZAo4zb3xjDYdhy7BxJ8NjMKmCa/84lYpEQP/A0S5DUlOP2N8ZwGLYcO0cHKJQqAGbe+DnGzQe2A0ZSvd8xhsOw5dhpq+SCMRznGi+9Yju3XLWD6y4cD3ophhYYs27YcuiSXDCS6ucaqXiED7zu2qCXYXDBeByGLceuUaskVwRGB0yS1WDoNMZwGLYc2uMYGYgSMQqrBkPHMd8qw5ZDTwI0FVUGQ3cwhsOw5UjGIwwnIiYxbjB0CWM4DFuSS7al2DM+GPQyDIYtiamqMmxJPvTfDxKNmH2RwdANjOEwbEkmUmaAk8HQLbq2JRORPSJyl4gcFpFHReTNDff/mogoEZlscvwxEXlYRB4QkXvrbh8XkTtF5HH7/7FuvQaDwWAwbKSbvnwJeKtS6grgeuCNInIALKMC3AQcd3mOlyilrlZKHay77W3AV5RSlwJfsX83GAwGQ4/omuFQSk0rpe63f14FDgO77bvfB/wGoNp46lcAH7V//ijwys2t1GAwGAx+6En2UET2A9cAd4vIy4FTSqkHXQ5TwJdE5D4ReUPd7duVUtNgGSdgWzfWbDAYDAZnup4cF5EU8GngLVjhq98CbvZw6POVUqdFZBtwp4h8Tyn1Xz7+7huANwDs3bvX97oNBoPB4ExXPQ4RiWIZjY8ppe4ALgYuBB4UkWPABcD9IrKj8Vil1Gn7/7PAPwPX2XfNiMhO+/l3Amed/rZS6oNKqYNKqYNTU1OdfWEGg8FwHtPNqioBPgwcVkq9F0Ap9bBSaptSar9Saj9wEniWUupMw7FJERnSP2N5KI/Yd38WuN3++XbgM916DQaDwWDYiCjVTn7awxOLvAD4GvAwULFv/k2l1OfrHnMMOKiUmhORXcBfK6VuFZGLsLwMsMJpH1dKvcM+ZgL4JLAXqyrr1UqpBZe1zAJPtfEyJoG5No7rFv22HjBr8opZkzfMmrzRqzXtU0ptCNl0zXBsBUTk3oZS4EDpt/WAWZNXzJq8YdbkjaDXZDQZDAaDweALYzgMBoPB4AtjOFrzwaAX0EC/rQfMmrxi1uQNsyZvBLomk+MwGAwGgy+Mx2EwGAwGXxjDYTAYDAZfGMNBtVnR0IJ+fI/6dE3mO+WCiPTdHKC6huO+O6f6kfP2JBeRK0XkxQCqTxI9dVIq4aDXAiAiV4jIc6Gv3qOrROQHRSTSR2t6uoi8FUApVXF7fC8QketE5A/7yZCJyHNF5EPAs4Nei0ZEniUi/wT8LPTHeS4iV4vI652kmPqFvrP83cb+Iv0ZcANwXEReCnxGKXWviISC+OLbQpAfAH5KRJ6plHpYRMJKqXKv12KvZwR4N5Y+2KyI3A38jVLqSBDrsdc0BrwDeB7wBHCjiPylUuqJoNZUxzuAHxSR+5RS/xnwZzcM/BHWxflvlVIVEZGgL4gi8nrgV4C/AL4b5Htkr2cCeDvW+zQOfNu+PcjPLop1bTqINYbiehH5oFLq7iDW04q+2Y30kDFgCLgC+ClgHniriKQC3C3+N+AE8CdYBoQgv1RYs1JEKfVM4OeBCWB/gOsBa015pdTVwM8BVwKBhhXqQi7/Bfwp8AdgfXYB7vR/C2tw2s1Kqb+w1xP4LhpLIui3lFIfUErlAj6/wdoYKaXU9Vjexm0Q+Pfu6cCIUupapdTrsK7P/SZ1ApwnhkNEbhKRm+xfh4HnAoNKqVks9d4F4I32Y3tyMbLXpOXlPwe8Tyn1P4C9IvJa+zE98wgb1vMB4HcA7B39KNZJ3VPsNf2g/evblVK/av98M9Yu8Uodm+7xmm4EUEqV7PPlB4EPAWdF5Ofs+yo9Ppf0+/QRLMXobSLyYyLybhF5rYj0dLZA/ftke7BXAveIyA0i8kUR+U0R+RH7/iC+c7+glPoV++dZ4JCIPK0X63BYk742lYHXiMiI/d5cD7xURK6xH9s3+ZctbTjsPMYngN8EFgGUUk8C38CaDwIwDdwBXCMiu7q9O2tY04K9phVqO4v/AbzLvr3UzbU0rOe36tZz0p6Fog3XGlZ4qCc0vEfz9pry9n0vAn4Na/rjq4DfEZELerymJfs2nWd5AMtj/APg10XkUyJyQY/PJf0+PYYVdvkC8EvAY8Cr7XUF8j4ppZaBZeDvsSZ2/gXW9+537NBsr96n+nM8X+cVKmAXkLUf3/ULdJNr04PAO7Hen78E/hDYA/yeiFzWJ56jhVJqS/7D2pHOA3/hcN9LgM8DF9q/XwX8LXBxUGuy79cNmd8Eftf+ORHgeiL2/18CrrF/DgX5HjU89gCWAXlJgOfSINYF8XLgPVgXy6/a94V7uSb92QADwE83vE9/A7yw1+9T3Zr2Yu3s31Z33x8Cv9MP5xPwH8Cb7J8lwPdJgN8GXmT/PgG8H/iZbq7J778t53GIyE47wbWA9UWO27f/tF2Ns08pdRfwXeCPAZRSjwD7gHwAa7pZLBl5AF1N9UrgV0Tk7cCfisj2INajrFDMJcCCUuq7IvKLwG+LyGgn1+NnTWJjr+8Qlrz0sU6vx+OaLlFKZe2HfwdIYRVd7BWRZ6guxMtbrQmrYOBSpdQalkEFqu/TDqwxBB3Hw5ouUUodx9qc/VjdoduwNkk9XVPD+aSvgZ/CCu+FlX3F7uWaqH12CmtS6msAlFLzwG7gUDfW1DZBW64OWvGXYs3/+DTwV/ZtA8CD1MJRfwo8hDWJMAZ8C/g/WEOi/hYrad6x3YbHNf2JvaZL6o67Eive+Z/A0wNaz6X2/Tdhhanuwgp/PC2Az636HmHtyCLAK4CvYMX0g/jc/hQrRHU58EPAZXXH3wbsDfJ9qjvu5fb79OEAz+9HgP32/Z/GCsd8275/Z4Dv06V1x/0BVuVgx9bSxpoets/xK+3v3B/bx30CmOrG2tp+TUEvoEMfzGXA3Vi7mW32Be6l9n23ALfXPfYjwP9n/7wdq7zz5QGv6cPAH9g/X4AV3/zxgNfzR/bPr8OKC9/YR+/RjcA9wCsDXtPfAL9X93uILoTyNvE+PQ+4rw/epw/XfeeGsYztzX2wpnfU/X4hcEsfnE/6s7sGq6LxVZ1eUyf+nbN9HNrFVFYJ7dXAPUqpfxKrjj0DPCYiUaXUF/Tj7cd+DuvCg1JqBpjpgzV9Hqt6QpRSJ4Ff6IP13Gg/zceVUn/fifV0YE0vtY/9MvDlPljTv+k11R0f9Jrqz6VvAtf2y5rsY1eAFeB7/bIm+/gngScDXtO/Ubs2fRcrnN6fBG252vkH/AxwGnvHgLWD+SZWSeRTWHX1fw/8X/t+nXS+HctFfcVWX1MH1vPKPnyP+nFN/XgumffJrKmr/wJfQBsfTAr4F+DNwP3YMXdgCvh14Bft3xNYVRzX2z+/Bytn8OytvqZ+W49Zk1mTWVMwa+rWv8AX0OYHtNf+/53AP9g/h7Dili+se9yfAS+z79t3Pq2p39Zj1mTWZNYUzJq68jqDXsAmP6QdWEnSW+zf/yeWxsvTsBprHgAuOp/X1G/rMWsyazJrCmZNnfx3zk8AFJGfB16nlHqh/fu7gZ1Ylvw3lFInzvc19dt6zJrMmsyagllTpzinDYeuRhBLFvkslmTAJ4GHldUEdd6vqd/WY9Zk1mTWFMyaOsk53TlufzCDWPXRrwGOK6XuCfKD6bc19dt6zJrMmsyazn3O2T6OOn4Jq4LhJmUL4fUB/bamflsPmDV5xazJG2ZNPeScDlXBuuaZvqHf1tRv6wGzJq+YNXnDrKm3nPOGw2AwGAy95ZzOcRgMBoOh9xjDYTAYDAZfGMNhMBgMBl8Yw2EwGAwGXxjDYTB0GBEZFZFfsn/eZTeBGQxbBlNVZTB0GBHZD/ybUuqqoNdiMHSDrdAAaDD0G+8ELhaRB4DHgSuUUlfZs6VfiTVb/iosOe0Y1qjZPHCrUmpBRC4G/hxLjjsLvF4p1ZHBRwZDJzChKoOh87wNeEIpdTXWHIZ6rgJ+ErgOeAeQVUpdA3wL+O/2Yz4I/LJS6lrg14C/6MWiDQavGI/DYOgtdymlVoFVEVkG/tW+/WHgGSKSwpoV/ikR0cfEe79Mg6E5xnAYDL2lXrOoUvd7Bev7GAKWbG/FYOhLTKjKYOg8q8BQOwcqpVaAJ0Xk1QBi8cxOLs5g2CzGcBgMHUYpNQ98Q0QeAf64jaf4KeBnReRB4FHgFZ1cn8GwWUw5rsFgMBh8YTwOg8FgMPjCGA6DwWAw+MIYDoPBYDD4whgOg8FgMPjCGA6DwWAw+MIYDoPBYDD4whgOg8FgMPjCGA6DwWAw+OL/B5EoT2vLKLUTAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot the first 100 time steps\n", "ds_trop_avg.tas.isel(time=slice(0, 100)).plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4. Nino 3.4 Region\n", "\n", "> Niño 3.4 (5N-5S, 170W-120W): The Niño 3.4 anomalies may be thought of as representing the average > equatorial SSTs across the Pacific from about the dateline to the South American coast. The Niño > 3.4 index typically uses a 5-month running mean, and El Niño or La Niña events are defined when > the Niño 3.4 SSTs exceed +/- 0.4C for a period of six months or more.\"\n", ">\n", "> — https://climatedataguide.ucar.edu/climate-data/nino-sst-indices-nino-12-3-34-4-oni-and-tni\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "ds_nino_avg = ds.spatial.average(\"tas\", lat_bounds=(-5, 5), lon_bounds=(190, 240))" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'tas' (time: 1980)>\n",
       "array([27.00284678, 27.06796429, 26.18095324, ..., 27.17515272,\n",
       "       27.30917002, 27.38399379])\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n",
       "    height   float64 2.0
" ], "text/plain": [ "\n", "array([27.00284678, 27.06796429, 26.18095324, ..., 27.17515272,\n", " 27.30917002, 27.38399379])\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00\n", " height float64 2.0" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_nino_avg.tas" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEiCAYAAADjxEWuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABXkElEQVR4nO29eZhbd3no/3klzUgzI82+j8e7HW+x48TZWZKQhCTs3JalwKW9bWkpULilUNr+ulwolNtSKBRKC4XSS1kamtCWEKDQhiVAFmexE8d2vNtjz75qmZFG0vf3xzlHI89ol46Wme/nefSMpLPonaOj8553F6UUGo1Go9Esx1FpATQajUZTnWgFodFoNJqUaAWh0Wg0mpRoBaHRaDSalGgFodFoNJqUaAWh0Wg0mpRoBaGpekTkrIjcXuC2R0TkFrs/p1oRESUiQRH5cIHb3y4iARGJr7Zjo8mOVhCaVY1SardS6ofF7kdEbhGRoRKIlGrfHxOREyLiF5FjIvI/s6z/SyJyzrzw/5uItGf5iH1KqT8sRDal1A+UUl7gfCHba2obrSA0msoTBF4BtABvBT4pIjelWlFEdgN/D7wF6AFCwN+WSU7NGkMrCE2tcJWIHBaRWRH5FxHxWAtE5OUi8rSIzIjIz0Rkb9KyhNtIRBpE5J9EZFpEjorI+1NYBSs+R0SagO8A/aa7JSAi/aX6x5RSf6KUOqaUiiulHgV+AtyYZvU3Ad9SSv1YKRUA/gh4rYj4cvksEdloup1+RUQumMfiN0XkWvP/nhGRT5fmP9PUOlpBaGqF1wF3AZuAvcAvA4jI1cAXgd8AOjDurv9DRNwp9vEnwEZgM3AH8OZcPkcpFQTuBi4ppbzm49LyDUXkA+YFNuUjl39SRBqAa4EjaVbZDRyyXiilTgERYHsu+0/iemAb8Hrgr4E/BG439/86EXlxnvvTrEK0gtDUCp9SSl1SSk0B3wKuMt//deDvlVKPKqViSql/AsLADSn28TrgI0qpaaXUEPCpPD4nK0qpjyqlWtM9ctzN32EogO+lWe4FZpe9NwvkZEEk8SGl1IJS6j8xXFxfU0qNKaUuYlgw+/Pcn2YVohWEplYYSXoewrhQAmwA3rvsTn0QSOUC6gcuJL2+kGKddJ9jOyLyl8Ae4HUqfRfNANC87L1mwJ/nx40mPZ9P8bps/7emetEKQlPrXAA+vOxuvVEp9bUU6w4D65JeD+bxOVnbHovIHyTFKFY8smz7fzDcWHcqpeYyrHoE2Je03WbADTyf27+h0eSOVhCaWufzwG+KyPVi0CQiL0sTtL0X+H0RaRORAeCdeXzOKNAhIi3pVlBKfSQpRrHikW47Efl94JeAO5RSk1nk+ArwChF5oRk8/yBwv1IqXwtCo8mKVhCamkYpdRAjDvFpYBo4iRnATsEHgSHgDPAD4F8x4hW5fM4x4GvAadOVVbIsJuAjwHrgRJLF8QfWQvP1C005jgC/iaEoxjBiD79VQlk0mgSiBwZp1ioi8nbgDUqpVZuxIyILGErwU0qpPypg+5cA92G4se5RSj1UYhE1VYxWEJo1g4j0YaS4/hwjxfPbwKeVUn9dSbk0mmrFVWkBNJoyUo9RJ7EJmAG+jq5C1mjSoi0IjUaj0aREB6k1Go1GkxKtIDQajUaTklUVg+js7FQbN26stBgajUZTMzzxxBMTSqmuVMtWlYLYuHEjBw8erLQYGo1GUzOIyLl0y7SLSaPRaDQp0QpCo9FoNCnRCkKj0Wg0KdEKQqPRaDQp0QpCo9FoNCnRCkKj0Wg0KdEKQqPRaDQp0QoiibmFRV739z/nwWeGKy2KRqOpMVZjXzvbFISIDIrIQyJyVESOiMi7zff/RUSeNh9nReTpNNvfJSLHReSkiHzALjktlFK8/xuHeezMFI+fnbL74zQazSri0IUZ9n/o+5wYXV2D/ey0IKLAe5VSO4EbgHeIyC6l1OuVUlcppa7CGERy//INRcQJfAZjRu8u4I0isstGWfnSz87y3SPGvPrAQtTOj9JoNKuM54bnmAkt8umHTlZalJJim4JQSg0rpZ40n/uBo8CAtVxEBHgdxhjH5VwHnFRKnVZKRTD69r/KLlmfOj/NRx48yu07u9na7SUQ1gqiEOJxxe/f/wxPnZ+utCgaTVmZDBiTa7916BLnJoMVlqZ0lCUGISIbgf3Ao0lvvxAYVUqdSLHJAHAh6fUQScqllMzOL/LOrz5Ft8/Dx35xH80el1YQBTIRCPO1x87z5UfStnbRaFYlE4EIbpcDl9PB3/3oVKXFKRm2KwgR8WK4kt6jlJpLWvRGUlsPAJLivZQRIBF5m4gcFJGD4+Pjecvndbt4/bWD/O2brqa1sR6vp4457WIqiOHZBQB+dnJyVQbsNJp0TAYj9LV4eN2BdfzrE0MMz85XWqSSYKuCEJE6DOXwFaXU/Unvu4DXAv+SZtMhYDDp9TrgUqoVlVKfU0odUEod6OpK2bE2I06H8Nsv2ca+wVYAfG4XgYXFvPejWVIQI3MLnJlYPWa2RpONCX+YTq+b33jRFuIKPv/jM5UWqSTYmcUkwBeAo0qpjy9bfDtwTCk1lGbzx4FtIrJJROqBNwD/YZesyfi0i6lgRpLumn52arKCkmg05WUyGKbDW89geyOvvmqArz52jm8dulTz1xI7LYibgbcAtyWltd5jLnsDy9xLItIvIg8CKKWiwDuB72EEt+9VSh2xUdYEXrdLZzEVyPDsAvUuB30tHn52aqLS4mg0ZWMyEKHD6wbgXbdtpbWhnnd97Smu/tD3+dUvPc7QdKjCEhaGbQODlFIPkzqWgFLql1O8dwm4J+n1g8CDdsmXDq/HRTASIxZXOB0pxdekYXh2gb4WDwc2tPPfx0aJxxUOfQw1q5xoLM5UKEKnqSA2djbx0w/cxsGzU3zn2RG+9LOz3Lilg1974eYKS5o/upJ6GV63oTNr3TSsBCOzC/Q2e7hpSwfToUWOjsxl30ijqXGmQ4soBZ3e+sR7Todw/eYO/uQVu/C6XQxN12bQWiuIZfg8WkEUyvDcPH0tHm7e2gnAz3UcQrMGmAwaNRAdTe4Vy0SEgdYGrSBWCz5PHaCrqfMlHleMzobpbWmgt8XD5q4mfnpSxyE0q5/JQASAjiQLIpl1bQ1cnNEKYlWw5GLSqa75MBmMEInF6W/1AHDTlg4eOzPFYixeYck0GnuZMKuorRjEcgbaGrhYo0FqrSCW4TVdTH5tQeTFiFkD0dtsKIibt3QSjMQ4PDRTQak0GvuZMC2IzjQWxEBrA3MLUeZqsL5KK4hl+NxaQRSCVTna19IAwA2bOxCBh0/oOIRmdTMZCONyCM2me3o569oaAbhYg3EIrSCW4dVB6oIYmTMtiBbDgmhrqufq9W1851k9W0OzupkMRGhvqk+b0j3QZtw01WKgWiuIZeggdWFcmlmg3umgo2nJzH7F3j6OjfhXXY98jSaZiUA4bfwBDBcTUJNxCK0gltFY50QE/NqCyIuR2Xl6WtyX3UXds7cPh8C3DmsrQrN6mQhG0mYwgRGbcLscNZnJpBXEMhwOwVvvwl+DAaVKMjy7QF9zw2Xvdfs8XL+pgwcOXdLdXTWrlsksFoSIMNBWm7UQWkGkwOvR/ZjyZWRuIRF/SOYV+/o5PRHkuWFdVa2pTaaCEZ69OJt2+WQgcplrNRUDrbVZC6EVRAp0R9f8UEol+jAt5649vbgcwrcOaTeTpjb58weP8tYvPpZyWTAcZX4xRqcvvQUBRiaTzmJaJXjdWkHkw1QwQiQaT6kg2pvquXlrJ9/SbiZNDaKU4scnxpkKRVKev4kq6iwWxLq2BiaDEUKR2rquaAWRAq+nTtdB5IE1KKi3pSHl8lfs6+fizDxPXZgpo1QaTfGcGg8wOhdGKQhFYiuWTwQzV1FbWJlMl2rMzaQVRAp8bh2kzgerijqVBQFw5+4e6p0Ovq2zmTQ1xsMnlvqJpfIqTPhzUxDrarQWQiuIFGgXU34Mz2VWEM2eOq7f3K6b92lqjodPLnUCSOVVmAxmbtRnUavFclpBpMCns5jyYmR2HpdDMt5FXb+pnWMjfqbMH5RGU+1EY3EeOT2ZuPFJddM4aTbqa88Sg+j2eXA5pOYymbSCSEHyVDlNdoZnFuhp9mScHnfD5g4AHjszVS6xNJqiODQ0QyAc5aW7ewEjY2k5E4EIPrcLT50z476cDqG/taHmMpm0gkiBniqXH+lSXJPZu64VT52DR07r5n2a2uDhE5OIwJ27eoDULqaJQDhriquFMTiottptaAWRAj1VLj/SFcklU+9ycM2GNh7VFoSmRvjpyQn29LckurGmdjFlL5KzGKjBwUFaQaTA69YN+3LFKJKbp781dYprMtdv6uDYyBwzIR2HyIVjI3Mp3Roa+wmGozx5fpqbt3YudXhOkdk4GQxnDVBbrGtrYMwfJhxdmS5brdimIERkUEQeEpGjInJERN6dtOxdInLcfP8v0mx/VkSeEZGnReSgXXKmYsmC0Kmu2ZidX2RhMU5Pc2YLAow4hFI6DpEL08EIr/ibh/nHn56ptChVzWIszge/9RznJ0vrunnszBTRuOIFWztpchvxhWCKOojJQISOLCmuFgOtDShlxOxqBZeN+44C71VKPSkiPuAJEfk+0AO8CtirlAqLSHeGfdyqlCp7bqSeKpc7VppfumlayewbbMHtcvDI6SnuNAN/mtT8+MQ4izHFuRJf+FYbz16c5Ys/PcOof4HP/NLVJdvvI2cmqXc6OLCxDbfLSb3LseJ6EI3FmQpFstZAWFiprhdn5tnY2VQyWe3ENgtCKTWslHrSfO4HjgIDwNuBjyqlwuayMbtkKBQ9VS53rGOUbppWMm6Xk6vXt/HoGR2ozsYPj48DS1XqmtQcHzFmjXznmWHOTARLtt+TowE2dzUlspOM2qjLPQrToUWUyu3mCGDQjGXUUqC6LDEIEdkI7AceBbYDLxSRR0XkRyJybZrNFPCfIvKEiLytHHJa6KlyuTM3b/xoLLdcNm7Y3MFzw3PMhrT7Lh3xuOJHz1sKoraCmuXm2IgfT52DOqeDv//RqZLt98xEkE1Jd/le98raqEmzzUZHU24WRG+LBxG4WEMuJtsVhIh4gfuA9yil5jDcWm3ADcD7gHtFJFUC/c1KqauBu4F3iMiL0uz/bSJyUEQOjo+Pl0RmPVUudxIWREN2CwLg+s3tKAWPn9VxiHQcvjjLVDBCt8/N8OyCbnKYgeMjfq7obeZ1Bwa578mhkijUxVic81MhNnctUxDhy2MQVqO+XC2IOnPi4rhfKwgARKQOQzl8RSl1v/n2EHC/MngMiAOdy7dVSl0y/44B3wSuS/UZSqnPKaUOKKUOdHV1lURuPVUud+YW8rMgrhpspd6l6yEy8dCxMUTgf1yzjlAkps/DNCilOD7qZ0ePj7e9aDNxBf/wk+KD+uenQkTjis2d3sR7Xs9KF9OEWUWda5AaoMvnYWwuXLSM5cLOLCYBvgAcVUp9PGnRvwG3metsB+qBiWXbNpmBbUSkCbgTeNYuWZdjTZXTFkR2rKaGucQgADx1Tnb2NXNcz6lOyw+fH2ffulZ29TUDS80QNZczHggzFYxwRa+PwfZGXrmvn689dp7pItu5nBk3YhmbVlgQy1xMObb6Tqbb52bMrxUEwM3AW4DbzFTVp0XkHuCLwGYReRb4OvBWpZQSkX4RedDctgd4WEQOAY8B31ZKfddGWVfg9eiOrrkwNx/F6RAa6zO3GkimpUG3U0/HZCDM4aEZbr2iO1GdrgPVqTk2bNxk7Oj1AfD2W7YQisT48INHi3LLnZ4IALA5SwxiJhRBJHf3KlgKona+T9vSXJVSDwPpmvO8OcX6l4B7zOengX12yZYLuqNrbvgXFvG6XaQOI6XG53bVXF/8cvHjE+MoBbdc0ZUowBrWxyolVgbTFaaC2N7j4123beVv/vskXT43v3fXjoL2e2YiSHtTPa2NS5aB4WK6PAYxHVqkpaEOZ4YeZMvpafYw7g8Ti6u8tqsUdtZB1DR67GhuzC1EaW7I7zRKdTemMXjo2DgdTfVcOdBCTClEtAWRjmMjfrp87stiAL9zx3YmgxE++8NTtDfW8+sv2pz3fk+NBy+zHsC4qVmZ5hqhrTF39xJAd7ObuDIyoLp92YtLK41WEGnweuoSKZya9PgXFvG5czexwbob0wpiObG4Md7ytiu6cTgEB0KX161jEGk4PjqXcC9ZiAgfetUeZkOLfPjBo5ybCvLS3b1ct6kdtys3N+iZiSC3XnF5wkuT28XCYpzFWJw6p+GZnwkt0tqY37nfbTb2G5vTCqKm8bldXKyhgpZKMTdfoAURjhKPq4wtwtcaZyeDzIQWuWFLR+K9vhZPYiCTZolYXHFiNMBbbtiwYpnTIXz89ftwuxzce3CIf37kPI31Tv6/l+3il65fn3G//oVFxv1hNiVlMMFSh+dgOJpwPU2HIokLfq50mUphvEYC1bpZXxp0DCI35hYWE3UjuWKlxAZrbIC73ZwaM4Kj27qXLk69LR5GdLHcCs5OBglH44n4w3LcLicff/1VHPrjO/niLx9gc1cTn3noZNbg9Wkzgym5BgJSF8/OhBbzdjH1NJsWRI0EqrWCSIOeKpcb/oVozimuFnreRmpOT1gXpyUF0dfSoGMQKbAC1Dt6mzOu11Dv5LYdPbz5+g1cnJnPml5ttetIFYOAy8/Z6VDkskB2LnSZFsdojdRCaAWRBj1VLjcMCyI/F1OT9WPTCvgyTo0F6PK5aUlKm+xt8eBfiGpluoxjI34cAtt6vNlXBm7bYfQE/a+jmVu/nR4P4BBY39F42fvLz9lwNEYoEqMtzxiE2+WktbFOWxC1TsLnqN0gaYnHFYFwNK88cEjqlqsvepdxajyw4s7VqoXQbqbLOT4yx8aOpqyjPi26mz3sXdfCfx0dzbje6Ykg69oaVwS0l5+zVi+x1jyK5BKy+Nw1U02tFUQafLrld1YCkShKQXOeFoRPWxArUEpxajzIlu7L74h7m3WxXCqMHkyp4w/peMmOHp66MJNokZGK0+PBFfEHWDpnrQFO06aCyNeCAKMWolaqqbWCSIOeKpedfFp9J2PdjelpaUtMBSPMzi+ypetyBWFN6tMKYolQJMq5qVD+CmJnN0ottVJfTjyuODMRvKwHk8XSVDlLQRhtNvINUoMRhxirkcw0rSDSoKfKZSffVt8WlvtOu5iWOGVmz2xZdvfabWa96FqIJU6MBlCKFTUQ2djd30xPszutm2nUv8D8YuyyHkwWTcuC1NbY3HzrIAC6fR7GA+Ga6NKrFUQa9FS57FjHJu80V22dreD0uJHiutyCcLucdHrrtQWRRK4ZTMsREW7b0cOPnx9PORfaSnHdkmLaW1P95deDJRdTYTGIxZhK7KOa0QoiDXqqXHYsCyLfQjlrxq/OzFni1HgAt8uRcCkl09vi0YODkjg6MkdDnZP17Y3ZV17G7Tu7CUZiKeeiW2nGqSwIp0NoqncmxSAKdzF111AthFYQadBT5bLjD1supvwsCJfTQUOdUx/bJE6NGxPMUjVw621u0C6mJI4NGwHqQqrwb97aiafOkTLd9clz07Q01CUSA5aT3CJmJrSI2+WgIY8uxhY95v5roRZCK4g0eHWmTVbm5q0gdf4dW4x26vrYWpweD6xwL1n0tXi0i8lEKcWxkTl29uUXf7Dw1Dl50bYuvv3MMIuxeOL9YDjKd58d4Z4r+9J2Jm5yuxJxs+lg/o36LJb6MVX/d6oVRBqa6l3GVDk9EyIt/oXCLAiwumNqBQFG0dX5qdCKALVFX6uH2flFQromhzF/mOnQYt7xh2Ref+0g4/7wZVbE946MML8Y47VXD6TdzpfUhXi6gEZ9FlaTvlpIddUKIg0Oh+Bzu5jTd7lpmVuI4qlzUO/K/zTyelwEtPIF4NxkiLhiRQ2ExVKxXPXfcdrN0eE5IP8MpmRevL2LvhYPX3vsfOK9bz51kcH2Bg5saEu7ndfjSsQgZgpo9W3RUO/E53bVRMM+rSAy4PPUJWYu1yrvvfcQH7jvsC379hfQqM+iqX5tWRCzoUW+8PCZlK1b0mUwWfQ2G4FrrSCMFhuQfwZTMi6ng188MMiPT4xzYSrE6NwCPz05wWuuGsg4+Cq5ged0KEJbU2HnPhiBah2krnGaG+oSfvZa5ccnxrn/yYvM2jDbYm4+WlD8AaovBvF/vnWEf/jJadty0//6v57nQw88x9MXplcss2ogNqVIr4QlC+KSVhAcG56jv8VDS4HuHYvXXzsIwL0HL/DvT18kruDV+9O7l8CMQSwsBanzbdSXTLfPUxNBaj0PIgPNHldNWxBzZm97gO8/N8ovXLOu5Psv1IKophhELK74yiPnicTijM4t8Af37MxrhGo2poIRvv7YBQCeHw1wzYb2y5afGgvQ1+JJFGMtp1f3Y0pwbMTPjr7CrQeLgdYGbtnexb88foG2xnr2DbZe1kU3FdY5q5RiZn6xoDYbFt3Nbp48v/JmodrQFkQGDAuidhWEVfgD8MDhSyXfvzFutLAfSTVNlRuenScSi7O128vnf3KG37vvMNGkDJdi+dJPzzC/GKPOKTyfot30qYlgWvcSGJk3bY11a96CiETjnBwLFBV/SOaN161nzB/m+Kif12axHmApBjG3ECUWVwXHIGCpYV+1V1NrBZGBZk9dVjfIsZE5jo3MlUmi/LB827fv7ObhExNMByMl3b+/gFbfFl638WOrhh/IuUljcuAHX7mb337JNu49OMSHHniuJPv2LyzypZ+d5aW7e9jR28yJ0cBly5VSnB4LpGwQl8yWLi/Pj2SeZbDaOTUeIBpXJbEgwGgB3u1z43IIr9jXn3X9JreLaFwlYkEtBd4cgVELEY7Gqz4JRiuIDDQ3uDJaEAuLMd76xcf41S8drMq5EafHgzgdwm/dupVoXPG9IyMl3X+xMYjFmCIcLd2deqFYCmJDZxO/c8d2Xrt/gH99YihlO4Z8+eqj55lbiPJbt2xlW493hQUx7g/jD0czWhAAewZaeG54rirPs3Jh3YjtLJEF4XI6+ONX7OL37tpBew5tu63uCkPmKOJiLIiuGqmFsE1BiMigiDwkIkdF5IiIvDtp2btE5Lj5/l+k2f4uc52TIvIBu+TMRLOnjkDEmJ2cim88McToXJiLM/P88HjmQSSV4PREgPXtjewfbGVDRyMPHB4u6f79C4t5d3K1SDWhq1KcmwxS73LQZ1a4vnxfH8FIjEdOr2zHkA8LizH+4eEzvGBrJ/sGW9ne42PMH07MEgB49tIsQNbOpHsGWghFYpyZCGRcbzVzbNhPvdORNphfCC/f28+vv2hzTuta3RWGpo1YUFFZTDVSC2GnBREF3quU2gncALxDRHaJyK3Aq4C9SqndwMeWbygiTuAzwN3ALuCNIrLLRllT0txQh1Kpu45GonH+7oen2DfYSrfPzZcfOVdu8bJyejzI5s4mRISX7+3jZ6cmMvbCz4dwNEY4Gi/cxeSpnkr1s5NB1rc3Jlo33LSlk4Y6J99/rjiL6ycnJhj3hxMXoO3m9LPnx5asiINnp3E5hH3rWjPu68qBFgCeuThblEy1zNERP9t6vLiclXF8WCMALkwZFkRRWUw10o/JtiOtlBpWSj1pPvcDR4EB4O3AR5VSYXNZqlvv64CTSqnTSqkI8HUMpVJWLPdJKjfTN58a4uLMPO95yTbecN16fvT8eOLEqQZiVm9707f98r39xBV859nSuJkSsyAKDVJbHV2rwoIIsTFpxKSnzsmLtnfyg+fGioqRXJox7jR39xs+823dhpWQ7GZ64tw0u/ubs/b02dLVhKfOwTND1RnvKgfHhueKqn8oFqvJ5IUSuJisfkzVPlmuLKpYRDYC+4FHge3AC0XkURH5kYhcm2KTAeBC0ush871U+36biBwUkYPj46kHgRSKlcK5PNU1GovzmYdOceVAC7dc0cUbrxvEIcJXHj2fajcV4dLMPOFoPJG6t6PXx5auJr5domymQmdBWHirpFuuUopzkyHWt1/utrhjVy8jcwtF3bGPzi3gcgjt5oVkoLWBxnpnIlC9GItzaGiGqzNU71q4nA529jUnXFJrjclAmDF/uOAeTKXAalM/ND2PSHFBaq/bRWO9s+prIWxXECLiBe4D3qOUmsOovWjDcDu9D7hXViadp0pCT3krp5T6nFLqgFLqQFdXVwklX2pjvbxY7j8OXeL8VIh33bYVEaGvpYGX7Ojm3oMXShLYLAWnllXnigi37+zhiXPTLCwWL2Oh0+QsfFXSLXfcH2Z+McbGzstbR9+2oxuHGPUjhTLmD9PpdSdcVw6HsK3bywnTxfTcpTkWFuMcWFYXkY4rB1p47tJc2pjYaua4aXXlO0WulFhu0QtTIZo9dSk77+ZDT7OH0bXqYgIQkToM5fAVpdT95ttDwP3K4DEgDnQu23QIGEx6vQ4ofSJ/FprTWBD/7+fn2NHr445dPYn33nLjBqaCEb7zTGkzhQrFqoFITp+8dmM7izHFoQszRe9/rohGfZA8oauydSZnrQymjsstiPameg5sbC9aQfSYvmaLbT0+njctiCfOGYVSV29ozWl/ewZaCISjnJkMZl95lTFqZvsMpJiXUS4sq3duIVpUkZxFT3P1jx61M4tJgC8AR5VSH09a9G/AbeY624F6YGLZ5o8D20Rkk4jUA28A/sMuWdNhmZDLYxBD0/PsX992WbXtzVs62djRyNcfrw430+mJAM0eFx1J6XvXmK6Mg+eKr+BcikEU52KqdJD6nHmx3ZBi+Mydu3o4NuIvOLY0NrdAl+/y2QLbe7yM+8PMhCI8cX6agdYG+lpyu+jt6TcC1c+uwUD1VND4DeaSjmoX3qRK92IC1BY9zdXfbsNOC+Jm4C3AbSLytPm4B/gisFlEnsUIPr9VKaVEpF9EHgRQSkWBdwLfwwhu36uUOmKjrClZsiCWLmLxuGI6FLnswguG++DGLR0rCqEqxenxIJu7vJcpsbamerZ1ezl4trj0TUiOQRTnYqr0XOpzkyGcDmGgbeVF2rIQ/7NAK2LMH05kq1hs67EC1QGePDedU/xhaVsv9S7HGlUQYZwOKdilWQo8dY6EW6k0FoSHkbmFqigWTYdtvZiUUg+TOpYA8OYU618C7kl6/SDwoD3S5YY3RRbT7Pwisbiiw7vyDqKn2cNkMEI4GsPtyn/SVCk5PR7kpq0dK94/sLGdBw5fIh5XBU3ksliKQRR2CrldDlwOqbgFcXYyyLq2BupSpE5u6Ghie4+XHzw3yq++YFNe+41E40wFI/SssCAMBfGj58cYnl3I2F56OXVmoHotprpOmQN6ijlni0VE8LpdzM4vFpXBZNHT7CESjTM7X1zjPzvRldQZcCZmQiwpiMmgYRKmMnWtrpuVTl0LhKOMzC2krM69dmMb/oXoZbn4heBfWERkaZh7vohIVfRjOjcZWhF/SOamLZ08fWEm7wrmcbPeZLkF0d/ioaneyb0Hh4Alt1+uXDnQzJGLay9QPRlYabVXAsvNVBoXk3FujFRxHEIriCw0N1zej2kyYPQz6vS6V6y7NGu2sl/4GTNAnWpC2bUbjYyZx88WF4eYW4jidbuKuqPzVrijq1KKs5PBlPEHiz0DLcwvxhJ9rXLFCj5a4yUtRIStPT7G/WEa6px5N57b09+CPxzlfBXV3JSDqWCkovEHC0tBlMLF1FsDs6m1gsiCz3N5P6ZJs+FdKhdToi1zhRXEabMdQ6r2xevaGuj2uYuOQ8wV0WbDwps0wrESzIQW8S9E2dCRXkEUWsFstVCwbhqS2W5OjrtqsDXvquA9a7SieioYoT3Fb67cWG7n1hIoq2q5ocyEVhBZaF42VW4ykN7FZN0RVHry16nxIA4h5YVPRLh2YzsHi7Ug5qMFF8lZ+CrsYjprZjBtzOBi2tLVREOdM38FkcaCgKU4xIGN+bmXrG3rnWsvUD0ViiQKDitJKS0Iq2HfaBW3cdcKIgtGR9ckF5NpQaQ6WVsa6nC7HBW/Izg9HmBdW2PaQPmBjW1cnJnn4kzhA2j8C4sFt9mwqLSLyXLTZLIgXE4Hu/qb874gj/nDOAQ6UrgirdYb123KrUAumXqXg+29Xo6uodbf0VicmdBidbiYPJaCKF4Wa85HNRfLaQWRhZUWRIS2xrqUrgERobfFw0iFfYpGimv6u2IrDlGMm2luofBW3xZeT11FXUxnJ0KIwGCGGAQYbqYjl/JrtT02Z1RRp6q2vXFLB/e9/SZesHV5fWhu9Lc0VPVdZ6mZNrvfpnLrlhtvvRWkLk26bU+zh5FZHYOoWZZPlcsWLOtp9lT0xxuLK06NBzLOF9jR66Op3lmUm6mYVt8WXreronUQ5yaD9DV78NRlTkkupNX2qH9hRQaThYhwzYa2gseaWvnza4Upy2pfZRYEGN9lNXd01QoiC80e4yJmpRVOBMIp3QYWvRX+8Q5NhwhH44nW0qlwOR1cvaGNx4uxIOYLnyZn4XU7K2tBTAYzprhaFBKoHpsLJ3r+l5reFg+z84sl6amVC6Ucv1oImVLLy411U1Q6BeGuuEs6E1pBZMGaCRGIGBeyyWDmfGzDxVS56kirktuq2E3HjVs6ODbiLygOEY8rAuHC51FbeN11zC/GKnIBsrq4Zoo/WBTSajtVH6ZSYQW+y3FhuTAVYs+ffo+HjlVuIJZlQXQ02XM88+H11w7yqTfuz9qePVd6mj2M+8MVV8Lp0AoiC4l2G6abaSoYyegLtaojZ0KVaUJ3YsxQEFu7M4+wfNmVfQA8WMCUuWAkSlwV3urbwjLXg+Hyd8A9MRZgMhhht2kdZMLldLCrL/dAdTQWZzIYXtGHqVRY6dTlyJ8/ORZgYTHOn3/naMXGnU5XkYupt8XDK3OYX50rPc0e4mop+aXa0AoiC8ktv6OxuNmHKbOLCSpXC3FizE9vsydrfGBDRxNXDrTwQAHzIYpt9W1hjR31V6Cj63eeGUEEXprUkTcTRqB6NqcK5olABKVSp7iWgp4ynmOWf/z50QD//vRF2z8vFdbFs1SB4WqipwSp8Q8cvsRf/efxUol0GVpBZCG55fd0aBGlMmdT9LZUtnz+5FiAbRniD8m8fG8fh4ZmOT+ZX1XuuL80PmFvBWdCfOfZYa5Z30Z3ikK2VOwZaCEYiXF6InurbeuimqpIrhQsTSMrg4IwrZQdvT4+8YPniUTL7wqZCkZoaahL2S+r1uktQbHcQ8fGuf9Je5T36jviJcbys/sXojn5QhPVkRXIZIrHFSfHAlndSxYv22u4mR54Jj8r4pxZP7A+B/99JirV8vvsRJBjI37u2tOb8zZXrsu91bbl+rHLgmj2uPDUOcpSkDnmD9PSUMfv3b2DC1Pz/EsF2tlni/vVMlacqhgFYSTO2HN8tILIgi+po2umKmoLK3OlEhbEpdl5QpFYYvZxNta1NXLVYCsPHMovDmHNR1ifpX4gG5WyIKy53PkoiK1dXiNQnYOCsCyIdGmuxSIi9DZ7GPXbH4MY8y/Q7XNzy/YurtvYzqf++yShSHm/r6lAdfRhsoMOs1ammHjSZDBsmwLVCiILyS4myxfamUFb17scdHrrK5K6tpTBlJsFAYab6bnhubya0Z2fDNHpddNYYCdXC5+7Mgriu88Os3ddC+vacldwrjxabY/NhRFJ3dCxVHSXqd7GmmkhIrz3zu2M+8Nln5pYLY367MDpELq8xaW6TgYitp1rWkFkYcmCiCYsiEx1EGBVR1ZAQZgtvLfl6GKCJTfTt/PIZjo3FcwpPTQbCQuijC6mizPzHBqazct6sNi3rpVnhmaz1h+M+RfoaKq31WfeW6Z5xsn1HNdsaMPpkEQzyHIxmSVzsNbpaXYX7HFQShmt0LWCqAwup4OmemfCgnAItGbJ/zeK5cpfPn9iNECXz51Xr/q+lgau3djGA3koiAtT80W7lyApBlFGC+K7pnvp7j19eW/7ou2dzC/GshYYjs3Zl+Jq0dPsZmTW3nobpRTj/nAiluJyOhhobeD8VOE9vAqRYToUKVlhWjXS0+wpeIbM3EKUSCye0atRDFpB5IDVbmPSNHWzzUDoafFUxsU0FsjLerC4c1cvx0f9ieykTISjMS7NzmftX5QL1rAhfxktiO88M8yOXh+bOrNXUC/nxs2d1Lsc/Oj4eMb17CySs+hp9hA2p5HZxUxokUgsnug6CkbcqZyzKObmo8TiatW6mKC41ikTpldDu5gqiNWwbzIQzulE7W32MBWMlK0VAhh3WicLVBBWd9FjI9krhS9Oz6MUGYfs5IrDITTVO8tmQQTCUZ44P82dOdY+LKeh3sn1m9r54fOZFcTo3IJtGUwWPWUYNmPNtEhOBR5sb0wkKZQDq83GancxFdo6JdMAs1KgFUQOWC2/jbGH2b+I3kSeevncTCNzCwTCUbZmabGRip19hoI4OpxdQZwvUYqrhddTvqFBE/4wSsH6HPovpePF27s4ORZgaDr1RTIWV0wE7OvDZFGO4VSJeo5lFsRUMIJ/oTzFjUuN+irfZsMueoq4XkwE7FWgWkHkgGVBZGuzYdFTgclyiQymAiyItqZ6eps9HB3OPmMgMUOhBBYElHcmxHTIutgUXpF7yxXdAPwwjZtpMhAmrrDfxeSzfxqZdcFKtiCs2NOFMsUhEhMcV7mLCQq7XkzWqotJRAZF5CEROSoiR0Tk3eb7fyoiF0XkafNxT5rtz4rIM+Y6B+2SMxeaGwwFMRHILd+4Eu02nh/NP4MpmZ19vtwsiMkQbpfjMr90MXg9dWVr+W31xypm4PyWribWtTWkVRCWW8buILVVY2FnqmvCxZT0XVvZa+WKQ1RTq2+7WOqtlf93OR6IIFKaCXepKC6RPTNR4L1KqSdFxAc8ISLfN5d9Qin1sRz2catSasI+EXPD53ExFYgQjMRySidLfOFlTHU9ORagvam+4HS3nX3N/OTEBOFoLO0kOjCqqNe3NxY8y2A5PreLQJndFcVkxIgIt1zRxf1PXiQSjVPvuvwe63tHjCyp/lZ7FYSnzkmrzdPIxvwLNNU7aXIvXSYGExaEVhClohhrcDIQpr2xPu/Z5rlimwWhlBpWSj1pPvcDR4EBuz7PTpo9dQQjRgAplxO12eOioc5ZXhdTHi02UrGjr5mo2aojExemcmuRnStet6tsWUwJF1ORKZMv3t5NKBJbMZHvCw+f4W/++ySvvqo/MUPCTnptnkZmFMldruhaGupoaagrqwXRWO/MOtSplmlucBU8qtjONhuQg4IQkV80LQBE5P8TkftF5Op8PkRENgL7gUfNt94pIodF5Isikm5yuwL+U0SeEJG35fN5pcbq6AqZq6gtlkaPlteCKEZB7OozgtvHMsQhlFKcnwqVJMXVoq2pLnHhtpvpkFHHUmyb8pu2dFDvdFyWzfS1x87zoQee4+49vXzsF/eVzMLKRLfN08jG58IpXYnlTHVdzVXUFiJCf2tDQXEdO6uoITcX0x8ppb4hIi8AXgp8DPgscH0uHyAiXuA+4D1KqTkR+SzwIQwF8CHgr4D/lWLTm5VSl0SkG/i+iBxTSv04xf7fBrwNYP369bmIlDfJba1zdeH0NLvL5mIKRaLMzi+yrq2h4H1s7GjC7XJkjENMBCKEIrGSBajB8NVPBiNEY3HbzGSL6dAibY3Z61iy0eR2ce2mNu5/coiTYwEmAmGeuTjLrVd08ck37Lf9/7DobXZzLIe4UaGM+RfYk8ISWt/emFO8qhSs5kZ9yezqb+bp8zN5bzcRCHPlutaSy2ORy5lsJee+DPisUurfgZy+MRGpw1AOX1FK3Q+glBpVSsWUUnHg88B1qbZVSl0y/44B38yw3ueUUgeUUge6urpyEStvkien5Xo3U87Ro1Zbj94i2ku7nA6u6PVxNEMtRKlTXMEIgCplKB+7mQ5GSjZT4PXXrsftcjI6t0BbYz2/9oJNfPbN16yISdhJT7OHiYA908iUUoymGZs62N7I0PR8WQYITQVzqz2qdfata+HizHwiKylXDAvCvuOTiwVxUUT+Hrgd+L8i4iY315QAXwCOKqU+nvR+n1LK6uvwGuDZFNs2AQ6llN98fifwwRxktYVkC6Izx3xsYxh5eeogEgqipbjA6M7eZr5/dBSlVEoXyfkpYxZCKdpsWCRywP0LRcufjelQ6dwVr9zXX9LJYoVgTSObCERKfuwC4Sjzi7GUHWnXtzcSicUZnVugv7VwqzUXpgIRruhptvUzqoF9phVweGiWW3d057TNwmIMfzhqq4spl9ud1wHfA+5SSs0A7cD7ctjuZuAtwG3LUlr/wkxfPQzcCvxvABHpF5EHzW17gIdF5BDwGPBtpdR38/nHSokVg3A55LJ4ROZt6ohE44Sj9ldTW5ZKMRYEwI4+H1PBSFrFdn7S8JHm0wU1G1YKZTmKCqeDi0WluFYbdk6Ws86BVPUc1g1COeIQU6HV3ajPYs9ACw6Bpy/M5LxNOWpEsl7tlFIh4H4R6RYRy8l/LIftHgZSOXsfTPGe5VK6x3x+GtiX7TPKhWVBtDfV5xx8TB6G4/bam4ExXCoLIqmiOtU0tPNTIXqbPSXNKLHuUMthbU2HIlw12Gr755SLUkwjS0eiSC6FiylZQdywuaPkn20RikRZWIyvCRdTk9vF1m4vh4dmct5mwm9vkRzk5ip6pYicAM4APzL/fsc2iaoQKwaRT41BOTuVjs4t0NJQV/R8hp29loJIncl0fipY0vgDGCe3iL0VwbDUFbS1iCrqaqMU08jSkRh6lCKLqa/Vg9MhttdCWH2Gik1LrhX2rWvl0NBszh16y9GnKhcX04eAG4DnlVKbMGIRP7VNoirESovMx5SztilHjv/w7ELR7iWAlsY6Blob0maonDeL5EpJndNBe2O97RZEMBJjMaZW1cVmaRpZ6RXEuD+9BVHndNDf6rHdxbQWiuSS2TvYylQwwtB0bumuE357G/VBbgpiUSk1CThExKGUegi4yjaJqpA6p4OGOmdemtpbRgUxOle6AO+O3stbbsTiioXFGJOBMKNz4ZIrCIAun5txm4ffTJegirracDqEbp/blmK5MX+YepcjbcytHLUQCQWxBmIQAFeZgepDObqZJoL2u5hy8UnMmLUMPwa+IiJjQHl6I1QRL9vbx81bc/e3+tyGK6McLqbh2YWEe6hYdvY189DxMe755E8YmVtI/EgtNhYwRyEb3WXI+LKK8dpW2d2oXcVyY2bL8nQxt/XtjXz/udGSf24y41YjulXcyTWZK3p91DsdHB6a5eV7s2fITfgjNNU7aai3L8aZi4I4BIQwso3eBLQAhZfs1igf+8X8YuaJcZphe3XpYizORCBcMgvirj29/Pz0JK0NdVy1vpVOrxu3y0GdU2hyuwqepZCJbp+b50eyd5IthmmzUZ9dTc0qxaaORh4+OZk2NblQxpImyaVisL2RiUCEYDh6Wa+mUnJ6PEidU2zva1Ut1Lsc7Opv5lCOmUyTwbBto0YtcvlmbzWL2uLAPwGYKaqaDJQrBjFmzjgolYLYM9DCfW+/qST7ypVun5uJQJh4XBVd5ZyOhItplVkQN23p5N+evsTxUT87SmRFgnFebe1Kfx+YaPs9HSrp5yZzcszP5k5v2SrTq4F961r4xhNDxOIKZ5bfgt1FcpAhBiEibxeRZ4AdZt8k63EG0AoiC1YWk90KolRFcpWk2+cmGldM2diTKeFiWkUxCIAXbOsE4OETpW16PDq3kLJIziKR6jppXxzi5FiArT1ry1mxb7CVUCTGqfHMTTPBatRnrwWRSTV/FXgF8O/mX+txjVLqzbZKtQqw3DJ2xyBK0Waj0nSXYQLfdNDom9/SsLpcTP2tDWzpauInJVQQC4sx/AvRlLUwFnYXyy0sxjg/FSp4vkmtstcMVOdSMDdhc6M+yKAglFKzSqmzSqk3KqXOJT2m0m2jWUJEjGlpdlsQZopjX41bEICtnUmnQ4u0NNRlNdtrkRdu6+LRM5Mlq9q3FHWmoVAtDXV0et0cuWRP077T40HiiqI6FNcimzub8LldWQvmYnHFVDBcOReTpnh8njrbZ/eOzM7jdjlq+s7YyrW3M5NpKhRZVTUQybxgaycLi3GeODddkv1lKpKzEBGu39zOI6cncy7syocTY9aExPxnrNcyDoewd7CFp7J0dp0JRYgre1NcQSsIWynHvOWRuTB9LZ6yzB+wC8vXPW6jgpgJla6Ta7Vxw5YOXA4pWRzC6qyb7eJzw+YOhmcXbHEznRwL4HQIGztLX3dT7Vy9vo1jI35CkfTXDus7srtPlVYQNuL12D8tbWR2PqOvuBbw1DnxeVyM2dhuYyq4uGorcr1uF/vXt/LwydIoCCugn+3ic+PmdgAeOT1Zks9N5sRogA0djRnH365Wrl7fRiyuOHRhNu06VltwbUHUML6yWBALNR1/sOj2uW11MRkWxOpUEAAv2NrFMxdnE+m8xZDr7O4tXV46vW4eOV36sOSJMf+aC1Bb7F/fCsCT59O7DBNFhNqCqF18NlsQSilGZ8P0rAIFYff8jKlgZNUVySXzgm2dKAU/O1X83fx0jnOgRYQbbIhDRKJxzk2G1lz8waK1sZ4tXU08mSGmZDUy7LC5ylwrCBvxeuy1IKaCESKxOH017mICy4Kwx8U0H4kRjsZXXZFcMvvWteDzuHj45Hj2lbNgKNPcjpUVhzhXwnqIc5NBonG15jKYkrl6fRtPXZhJq3gnAmFcDrE9OUUrCBvxuutsTXMt1RyIaqC72cPYXNiWjJipVVokl4zL6eDGzR0lsSCm8pi8Z82DKGUc4sSYUSS2phXEhjamghHOplG8kwHjO7Kr84CFVhA24vO4iMTsmypntXnubbF37GM56Pa5CUfjzNmgUFdjJ9dUXDnQwrnJEMEirdbpYCRna2tLV5MZhyihghgNIGLEONYq12xoA0iZurwYi/PImUk2lHg2Syq0grARu/sxlWrUaDVgFWXZ0fZ7ZpU26lvO9l7DZ2/dgReKUTOS27FaikNMlcz6OzHmZ7Ct0dYupdXO1i4vPo8rZaD63oMXODcZ4jdfvMV2ObSCsJHksaN2MDK7gNMhGStea4VEsZwN7TYsF9NqTXO1uKLHUBDFdsadDi7mFa+5YXMHI3Oli0OcHAus2QwmC4dDuGqwdUWgej4S41P/dYJrNrRx245u++Ww/RPWMHaPHR2ZXaDLnCpW61jFcqO2WBCGgljNaa5gtOD21Dk4Plq4gghHYwTC0byqzksZh4jG4pyeCK65Jn2puGZDG8dH/Zd1Y/inn59ldC7M7921oyzFsVpB2IjdU+VGSjhJrtIk+jHZYUEELQWxul1MToewrdvH80UoiIQ7Lg8LwohD1PPY2eLrIS5MzxOJxtdsimsyV69vQykSBXOz84t89oenuOWKLq7b1F4WGWxTECIyKCIPichRETkiIu823/9TEbkoIk+bj3vSbH+XiBwXkZMi8gG75LSTZo+9U+VGSjSLuhrwul001DltqYWYCS3i87ioWwNzBbb1eItSEJYyzWf+uohw5UALz5Wgcd8JU/a1nMFkcdX6VkSMgrmjw3O8996nmZ1f5H0vvaJsMtj5i4kC71VK7QRuAN4hIrvMZZ9QSl1lPh5cvqGIOIHPAHcDu4A3Jm1bMyzNhLCnYd/I7OqxIESE7mZ7qqmngrmnbdY6V/T4GJ0LJ9xq+VLoYKU9Ay2cGAuwsFhcxt7xEa0gLJo9dWzv9vHZH57i7k/+hIdPTvDeO7azu7+lbDLYMysQUEoNA8Pmc7+IHAUGctz8OuCkUuo0gIh8HXgV8JwdstrF0tjR0lsQgXAUfzi6ahQEmMVyNvRjml7lbTaSsTKZnh8NFOSGKDSgv7u/mVhccWzEz1WDrXl/rsWhoVk2dzUlbq7WOi/f28c3n7rI668d5PXXDpb9PC6LzS0iG4H9wKPmW+80p9N9UUTaUmwyAFxIej1E7sqlarBzqtylmXnAGBizWuj2eWzp6DqdR9pmrWNlMhUaqC60ZsS6q332YvoGc7lweGiGfebQHA286yXb+O/fvYXfePGWitzk2K4gRMQL3Ae8Ryk1B3wW2AJchWFh/FWqzVK8lzLJWkTeJiIHReTg+HjxbQZKiZ1T5S6aCmJgFQ10X9fewND0fMHukXRMBxdXfZGcRV+LB5/bVXCq61TQcIfmG9Bf19ZAS0MdRy4VriBGZhcY84fZt658LhRNZmxVECJSh6EcvqKUuh9AKTWqlIoppeLA5zHcScsZAgaTXq8DLqX6DKXU55RSB5RSB7q6ukr7DxSJiNg2NGh4xnDFrCYL4jX7B4jE4vzL4xeyr5wH06HcK4NrHRFhe6+vYAtiKhimuYCAvoiwZ6C5qAlzh8wpanuLcFFpSoudWUwCfAE4qpT6eNL7fUmrvQZ4NsXmjwPbRGSTiNQDbwD+wy5Z7cSusaOXZuZxOiRRYLYa2NHbzPWb2vnyI+eIxUtTlbuwGCMUia36KupktvcYqa6FVDZPhQqfm7G7v4Vjw34WY/GCtj90YQaXQ9jV11zQ9prSY6cFcTPwFuC2ZSmtfyEiz4jIYeBW4H8DiEi/iDwIoJSKAu8EvgccBe5VSh2xUVbbsGuq3KWZeXqbPauiSC6Zt960kaHpef772FhJ9ldIXn+tc0WPl5nQYkHxnHz6MC1nd38zkVicE6OFtfo4PDTLFb2+rG3GNeXDziymh0kdS1iR1mqufwm4J+n1g+nWrSXsmip3cWaegVXkXrK4c1cPfS0e/t/Pz3LHrp6i9zcZNC6Sq3UedSq2JwWqu/Osk5kKRgoeQLVnwAxUX5plV39+VkA8rjg8NMPL9/UX9Nkae1j9lUMVptkmBXFpdp6+VRSgtnA5Hbzp+vX85MQEJ4tsOgdLc66tVh5rgeRU13wpJl6zqaOJxnpnQQVzZyeDzC1EdYC6ytAKwmbscDHF44qR2YVVFaBO5g3Xrafe6eDLPz9b9L4sBdHlXX3KNB2dXjcdTfWXZTLlEo9QShVVVOgw4weFpLoeHjK22atTXKsKrSBsxo6pchOBMIsxtWoVRKfXzcv39vGvTwwVHPC0SMzu9a0dFxMYbqbHz07xoQee4/aP/4irP/T9xKD7dMwvGpP3iqk63zPQwnPDc3knGRwamsFT51jzXVyrDa0gbMaOqXKrsQZiOS++ootgJMap8eLcTOP+MF63i8b6tVWZu7u/mdMTQb7883O0N9UzHVrkgcPDGbex+jAVE6/Z3d9MKBLjzEQwr+0OXZjhyoEWXGugX1Ytob8Nm7GmyhXboyaZS6uwBmI5O81Ux6PDxTWAG/eHV8W8jHx5123b+OqvXc/Tf3IH9/7Gjezqa+b+py5m3GY6WHzGl1VRnU/B3GIszpFLc9q9VIVoBWEzPhv6MVltNvpWwajRdGzubKLe5eDocHHDb8b9Ybq8a09BtDTWcdPWzoTl9Jr9Axy6MJPRIlvqw1R4zci2Hi/1Tgc/Oj7OfCS3m6LnR/2Eo3H26gB11aEVhM3YMVXu4sw8XreLZs/qdZu4nA6293iLtyACa9OCWM4rr+rHIfDvGayIUszurnM6uHVHF/c/dZGrP/R93v7PT2QNWlvzDnQPpupDKwibsWOq3PDsPP2tnrJMlKokO3ubi1cQc1pBAPQ0e7h5ayfffPpi2oymRAyiyKLCz/zS1Xz1167nF65Zx89OTfK+fz2ccf0TY34a6pxs6Ggs6nM1pUcrCJvxmUOD5krYj+nSzOpNcU1mZ18zE4EIYwWOIZ2PxPCHo1pBmLxm/wAXpuZ5YtmcY4vpUASHLA26KhSX08FNWzv50Kv38JYbNvD8qD9jDO78ZIgNHY2r/oanFtEKwmYSMYgSupguzcyvGQUBFByHmDDTOrWCMHjp7l4a6pxpg9WTwQhtjfU4Sti+Zc+AMSfieIbusuemQqxv19ZDNaIVhM2U2sW0sBhjMhhZlW02lrOzz6gILtTNZE2n0wrCoMnt4qW7e3jg0KWU9SXF9GFKR2JORJqspnhccX4qpN1LVYpWEDZT6qlySxlMq7cGwqK1sZ6+Fk/BCmKpilorCIubtnQytxBlZHal224qGCl5zyprTsSzF1N/h6P+BSLROOs7mkr6uZrSoBWEzVguplL1YxqeXf01EMns7Cs8UG1VUXdrCyKBZU2Np6iqng6Vfnb30pyI1BbEuckQABu0i6kq0QrCZtwuJ/VOR8kUxFIV9VpRED5OjQcLKjQc94cRKT4rZzWRUBApWoFPBRdtaYu+J8OciPOmgtioLYiqRCuIMmD0YypNFtOlmXlEjLTFtcDOPiPIWUhn13F/mI6met2+IYl0CkIpZVoQpR+stHugJe2ciHNTQVwOoX8Vt42pZfQvpwyUcqrcpZl5un1u6l1r46uzMpmeK8DNNO4P06njD5fR3lSPyFKGl8XcQpRYXNkyu3uPORsiVaD63GSIgbYGrcSrFP2tlAFfCWdCXJpZWNUtNpazsaMJT52joDiErqJeSZ3TQVtj/QoLYrpERXKp2NjRRFO9kyMpKqrP6xTXqkYriDLgdbvwlzCLaa3EHwCcDuGKAiuqJ9Zoo75sdHndKxSE1YfJjhiEwyHs7m/hSIpBQucmdYprNaMVRBnweUrjYlJKcclss7GW2NXn4+iwP6ehNxZKqTXbyTUbXT73iiym6RK0+s7E7oHmFXMiZkIRZucX2dCuA9TVilYQZaBUU+WmghEWFuNrysUEcEWPj9n5xZSpmemYm48SicV1DUQKunzuFTEIy6KwK+Nrd3/LijkRVorrem1BVC1aQZQBn6cOfwl6MZ0aN35cm7rW1h2XlbE14Y/kvM14wKgX0RbESjq9Rgwi2SK7MB3C6RDbCjD3DBiB6uR6iHNTZg2EVhBVi20KQkQGReQhETkqIkdE5N3Llv+uiCgR6Uyz/VkReUZEnhaRg3bJWQ6ssaP5uEhS8fyo0c9me4+vFGLVDB2mFTAZzN2C0G020tPlc7OwGL/Mqr0wZbgu7com2trlxe1yXNb6+/ykccOjg9TVi50DBaLAe5VST4qID3hCRL6vlHpORAaBO4DzWfZxq1JqwkYZy4LX7WIxpghH43jqnAXv5+RYgKZ6J/1roM1GMpbbw2pHnQuWy0RXUa8kuRbC6jZ8YTrEYJt9F2qX08GOvmYODyVZEJMhunzuNTcOtpawzYJQSg0rpZ40n/uBo8CAufgTwPuB4m6pa4ROr3GBW+73zZfnR/1s7fGtubbIS8cvfwXR5V1byjQXrGOSfDwvTM3bfid/85YOHj87xckxwxI+NxXSLTaqnLLEIERkI7AfeFREXglcVEodyrKZAv5TRJ4QkbfZLaOdWD70VA3S8uH50QDbu72lEKmmaPbU4XIIk3ko2HF/mHqng+YGfXe6nE6foXAtJRqKRJkIhBm0+WL9qy/YhKfOySd+cAIw2mzoAHV1Y7uCEBEvcB/wHgy30x8Cf5zDpjcrpa4G7gbeISIvSrP/t4nIQRE5OD4+XiKpS4uVdTQyV7iCmA5GmAiE11z8AYw8+ram+rxdTF0+95qztnLByuwaNwcxDU0b/b3WtdmbHdfhdfO/bt7Etw8P89T5aUbmFnSKa5Vjq4IQkToM5fAVpdT9wBZgE3BIRM4C64AnRaR3+bZKqUvm3zHgm8B1qT5DKfU5pdQBpdSBrq4ue/6RIuktgQVxwuxFtK1n7VkQAB1N9fm5mHQVdVraGutxOiSRNnzBzCYqR7D411+4GZ/Hxe/dZ4wh1RlM1Y2dWUwCfAE4qpT6OIBS6hmlVLdSaqNSaiMwBFytlBpZtm2TGdhGRJqAO4Fn7ZLVbpobXDTUOYtSEFYG07Y1aEEAdHrdTOWRxaSL5NLjcAid3vpE2vB5U0HY7WICaGms420v3MzzZuM+7WKqbuy0IG4G3gLcZqaqPi0i96RbWUT6ReRB82UP8LCIHAIeA76tlPqujbLaiojQ2+JhuAgX08mxAF63a81lMFm0N9UzWYCLSZOaTq87yYKYp6HOSUeZ2qL/ygs20dZoZE/pIHV1Y1sETyn1MJDRAWxaEdbzS8A95vPTwD67ZKsEvc2eoi2Ird3eNetT7/DWM5mji2kxFmcqFNFV1Bno8i31Y7owbTTMK9e55XW7+P27d3L/U0N6VkeVoyupy0RfS7EKIsC2NZjBZNHpdRMIR3MaHDQVjKCULpLLRHLDvgtTIQbby9u+5XXXDvL1t924Zm94agWtIMpET4uH0bkF4vH8Sz/WcgaTRT7FctaFT8+CSE+Xz81kMEw8rrgwFWKdjUVymtpFK4gy0dfiIRpXefnRLdZ6BhOQ8I/n4mayfOvagkhPp9fNYkxxZjJIMBLT7S40KdEKokwUUyy31jOYIL9+TBMJC0L7t9NhKc8nz00D5clg0tQeWkGUCatLZiHFcms9gwnysyAsK027mNKTUBDnLQWxtlrIa3JDK4gysVQsN5/3tms9gwmMLCbIzYKYDITx1DlorC+8MeJqx1IQT52fAbC1UZ+mdtEKokx0eN24HFKQBfH8aIDtazj+AEZqZL3LkVMMZyIQoaNJt9nIhKUgjo/66Wiqp8mte1ZpVqIVRJlwOoRun5vhPGMQs6FFJgJhtq7hFFcwig07mnKrhZgIhOnUAeqM+EyFqxSs0/EHTRq0gigjvWaqaz4MzxkuqYFW/SM2iuVycTFF6NQFWBkRkUQh4aDNTfo0tYtWEGWkr6UhbwtibM64IPY06zvijiZ3TnUQE4FwImahSY/lZtIprpp0aAVRRnrMdhv5jB61LI5u39rNYLLIpaNrPK6YCkZ0BlMOWApCp7hq0qEVRBnpa/EQisTwJ80CzoY1W7lbWxB0eLPPhJhbWCQaV4m6CU16OhMuJq0gNKnRCqKM9LTkXyw3NrdAS0NdUbOsVwsdXjfzizFCkfQK1hrrqovksqNdTJpsaAVRRvoKUBCjc2G6dUYOsNSPKVMmk+WC0i6m7Nx6RRev2NfPgA5Sa9Kgk5/LSCGT5cb8C4k2HWudzkSxXCSt39xSHjpInZ3969v4m/VtlRZDU8VoC6KMWHGEfIrltAWxREeT2Y8pQ6qr5WKy1tVoNIWjFUQZcbucdHrrc051VUox7g/TrS0IIDcX02QgjAh6EI1GUwK0gigzPc25F8vNhBaJxOLagjDpSHIxpWM8EKG9sR6nQ7fZ0GiKRSuIMtPX4snZghj1G+vpGIRBY72LxnpnRhfTpC6S02hKhlYQZSYfC8KqotY1EEu0N9VntCAmdZGcRlMytIIoM30tHqaCkZxmK1uKpEdXUSfo8LozKgijzYZWEBpNKdAKosz0txo552cmglnX1VXUK+lsytywbzIQSQwX0mg0xWGbghCRQRF5SESOisgREXn3suW/KyJKRDrTbH+XiBwXkZMi8gG75Cw3L9jaiQh8/7nRrOuOzS3Q7HHpKuok2jO0/F5YjBEIR/Usao2mRNhpQUSB9yqldgI3AO8QkV1gKA/gDuB8qg1FxAl8Brgb2AW80dq21ulu9nDN+ja+8+xI1nVH58I6QL2MDq/R0TVVw8OlGghtQWg0pcA2BaGUGlZKPWk+9wNHgQFz8SeA9wPp2ppeB5xUSp1WSkWArwOvskvWcnPXnl6ODs9xNoubacy/oN1Ly+j01hOJxVM2PJzUbTY0mpJSlhiEiGwE9gOPisgrgYtKqUMZNhkALiS9HmJJuSzf99tE5KCIHBwfHy+VyLZy155egKxWxOhcWAeol2EVwFkZXslY86p1mqtGUxpsVxAi4gXuA96D4Xb6Q+CPs22W4r2U1oZS6nNKqQNKqQNdXV3FiFo21rU1snddC999djjtOlYVdZe2IC7j6vVtiMA3nxpasWzCry0IjaaU2KogRKQOQzl8RSl1P7AF2AQcEpGzwDrgSRHpXbbpEDCY9HodcMlOWcvNXXt6OTQ0y8WZ+ZTLrSpqbUFczsbOJu7e08v/+/k5/AuLly2b0BaERlNS7MxiEuALwFGl1McBlFLPKKW6lVIblVIbMRTB1Uqp5b6Wx4FtIrJJROqBNwD/YZesleDuPX0AfDeNm8mqotYxiJW8/cVb8S9E+edHLs9xmPBHaKx30livmxRrNKXATgviZuAtwG0i8rT5uCfdyiLSLyIPAiilosA7ge9hBLfvVUodsVHWsrOps4kdvb60bqalWdTagljOletaeOG2Tr7w8JnLCg4ng7rNhkZTSuzMYnpYKSVKqb1KqavMx4PL1tmolJown19SSt2TtOxBpdR2pdQWpdSH7ZKzkty9p4+D56a5MBVasWxpFrW2IFLx9lu2MBEI840nlmIRkwHdZkOjKSW6krqCvPbqARrqnLz7608RicYvW5aootYxiJTcuLmDqwZb+dyPTxGNGcduIhDWcyA0mhKiFUQFGWxv5C9+YS9Pnp/hIw8evWzZ2NwCPo+LhnpdRZ0KEeG3btnChal5PvPQKcAYN6pnUWs0pUNH8yrMy/f289T5Gb7w8BmuGmzl1fuNco8xv66izsYdu3p47f4BPvGD51nX1sBUMKxdTBpNCdEKogr4wN07eGZolg/cf5jeFg83bO5gdG6BHp3BlBER4aP/Yy/Dswu8/77DxJVOcdVoSol2MVUBdU4Hn37TfvpbG3jzPzzKl39+1pxFrS2IbNS7HPzdm69hU2cTgG71rdGUEK0gqoRun4d/e8fNvGh7F3/070e4ODOvayBypKWxjn/85Wt5+d4+rtvYXmlxNJpVg1YQVUSzp47P/88DvPPWrQBs6fRWWKLaYbC9kU//0tX0tmirS6MpFToGUWU4HcLvvvQK3nDdIH0tDZUWR6PRrGG0gqhS1rU1VloEjUazxtEuJo1Go9GkRCsIjUaj0aREKwiNRqPRpEQrCI1Go9GkRCsIjUaj0aREKwiNRqPRpEQrCI1Go9GkRJRSlZahZIjIOHAuz806gQkbxCkGLVNuaJlyQ8uUG2tVpg1Kqa5UC1aVgigEETmolDpQaTmS0TLlhpYpN7RMuaFlWol2MWk0Go0mJVpBaDQajSYlWkHA5yotQAq0TLmhZcoNLVNuaJmWseZjEBqNRqNJjbYgNBqNRpMSrSA0Go1Gk5I1oyBERCotQy1QjcepSmVaM7+dQhGRqps3IyI+82/VnVPVyKo+yUVkt4jcAqCqJNgiIn3mX2elZbEQkZ0iciNU1XHaIyIvFRFXFcl0pYi8F0ApFa+0PAAicp2IfKSaFJaI3CginweurbQsFiJytYj8K/CrUB3nuYhcJSK/LiK9lZYlHVWn4UuB+WP5NHAbcF5EXgL8u1LqoIg4KvHjFhEv8FngTSKyTyn1jIg4lVKxcsuSJFML8DHgOmBcRB4F/lEpdbKCMrUBHwZuAk4Bt4vI3ymlTlVKpiQ+DLxURJ5QSv2wkt+fiDQDf45xEf6SUiouIlLpC5+I/Drw28DfAk9VwTneAfwpxnFqBx4x36/kd1eHcX06ABwFbhCRzymlHq2EPJmomruOEtMG+ICdwJuASeC9IuKt4J3fK4ALwF9jKAoq+cMxeT9GJts+4DeADmBjRSUyZAorpa4Cfg3YDVTUHZDkKvkx8Engz8D4/ip45/6HwA3AnUqpvzXlqfhdMbAe+EOl1GeVUgtVcI5/DOPQ3IBhPbwFKv7buxJoUUpdo5R6M8Z1uNpafACrSEGIyB0icof5shm4EWhUSo0D9wFTwDvMdctywTFlutN8+W3gE0qp3wHWi8gbzHXKasUtk+mzwB8DmHforRgnb1kxZXqp+fJPlVL/23x+J8Zd327Ld1xmmW4HUEpFzXPmpcDngTER+TVzWbzM55N1nL4IjAHdIvILIvIxEXmDiKwvhyzLZLrdfN6CodAfE5HbROR7IvIHIvJac3klfne/qZT6bfP5OPCciFxRDjlSyGRdn2LA60SkxTw2NwAvEZH95rpVEx+peQVhxhm+DvwBMA2glDoD/BR4j7naMHA/sF9E+u2+01om05Qp0xxLdwm/A/yF+X7UTllSyPSHSTINKaUuJSmpeQy3TllYdpwmTZnC5rIXAb8L/BPwGuCPRWRdmWWaMd+z4iBPY1iBfwa8T0S+ISLrynw+WcfpOIa75DvAbwHHgV805arIcVJKzQKzwD8Dr8ZwMw1jfHf7ynicks/xcJKVp4B+IGSub/uFOM316RDwUYzj83fAR4BB4IMisr1KLEEDpVTNPjDuLieBv02x7FbgQWCT+XoP8CVgS6VkMpdbxYk/A/6P+dxTYZlc5t//BPabzx2VlGnZurswFMWtFTyfGjEufDuAv8K4KP7IXOYsp0zWdwM0AL+87Dj9I/DCch+nJJnWY9ypfyBp2UeAP66G8wn4b+Cd5nOp4HES4I+AF5mvO4BPAb9ip0z5PmrSghCRPjPINIXxY3Wb7/+ymfmyQSn1EPAU8JcASqlngQ1AuAIy3Skim81VreylVwO/LSJ/CnxSRHoqJZMyXChbgSml1FMi8nbgj0SktVIyiYkp33MYbY/PllqeHGXaqpQKmas/DngxEiDWi8heZYM/O5NMGIH7bUqpeQzFCSSOUy9wvtTy5CjTVqXUeYwbsV9I2rQb44aorDItO5+sa903MNxyTmVemcspE0vfnQKiwOsAlFKTwADwnB0yFUylNVSeGvklwE8wYgp/b77XABxiyY30SeAwsAWoB34O/A3wLMaJ66OEdw45yvTXpkxbk7bbjeGL/CFwZQWOkyXTNnP5HRjupYcw3BZXVFCmrRh3WC7gVcB/YfjcK/HdfRLDtbQDeBmwPWn7twDrK3mckrZ7pXmcvlDBc/xZYKO5/D4MN8oj5vK+Sp/j5jp/hpGpVzJZCpDpGfMc323+5v7S3O7rQJcdshX8P1VagDwO/nbgUYw7k27zIvYSc9ndwFuT1v0i8H/N5z0YKZOvrLBMXwD+zHy+DsP3+PoqkOnPzedvxvDb3l4FMlnH6XbgMeDVFZbpH4EPJr12YIMLrojjdBPwRBUcpy8k/e6aMZTqnVUg04eTXm8C7q6C88n67vZjZBC+ptQyleJR1XUQllmojNTUq4DHlFL/KkYOeBA4LiJ1SqnvWOub634b4+KCUmoUGK0CmR7EyFQQpdQQ8JtVItPt5m6+qpT65yqR6SXmtj8AflAFMj1gyZS0faVlSj6ffgZcUy0ymdvOAXPAsWqRydz+DHCmwjI9wNL16SkMV3h1UmkNlUEj/wpwCVP7Y9yN/AwjzfAcRk76PwNfNpdbwd+3YpiVr9Iy5STTq7VM+rtb5cdpTchkx6PiAqQ5+F7g34B3A09i+sOBLuB9wNvN1x6MjIkbzOd/heHTv1bLpGXSMmmZ1opMdj0qLkCGL2G9+fejwNfM5w4Mn+ILk9b7NHCXuWyDlknLpGXSMq1FmWz5PystQA5fRC9GoPJu8/XvY/QvuQKj+ORpYLOWScukZdIyaZlK+6iJiXIi8hvAm5VSLzRffwzow9DK71dKXdAyaZm0TFomLVNpqXoFYUX+xWjVO4ZRJn8v8IwyCoW0TFomLZOWSctkA1VfSW0e/EaM3OLXAeeVUo9V8uBrmbRMWiYtU7XKVEqqug4iid/CyBa4Q5nN3KoALVNuaJlyQ8uUG1qmMlL1Lia4rMCkatAy5YaWKTe0TLmhZSovNaEgNBqNRlN+qj4GodFoNJrKoBWERqPRaFKiFYRGo9FoUqIVhEaj0WhSohWERlMgItIqIr9lPu83i6U0mlWDzmLSaApERDYCDyil9lRaFo3GDmqlUE6jqUY+CmwRkaeBE8BOpdQec/bwqzHmj+/BaPNcjzGiNAzco5SaEpEtwGcw2kSHgF9XSpVkwI5GUwq0i0mjKZwPAKeUUldhzAFIZg/wS8B1wIeBkFJqP8aM9P9prvM54F1KqWuA3wX+thxCazS5oi0IjcYeHlJK+QG/iMwC3zLffwbYKyJejFnS3xARaxt3+cXUaNKjFYRGYw/JPXniSa/jGL87BzBjWh8aTVWiXUwaTeH4AV8hGyql5oAzIvKLAGKwr5TCaTTFohWERlMgSqlJ4Kci8izwlwXs4k3Ar4rIIeAI8KpSyqfRFItOc9VoNBpNSrQFodFoNJqUaAWh0Wg0mpRoBaHRaDSalGgFodFoNJqUaAWh0Wg0mpRoBaHRaDSalGgFodFoNJqUaAWh0Wg0mpT8/9YPeqAkeRTfAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot the first 100 time steps\n", "ds_nino_avg.tas.isel(time=slice(0, 100)).plot()" ] } ], "metadata": { "interpreter": { "hash": "937205ea97caa5d37934716e0a0c6b9e4ffcdaf6e0f20ca1ff82361fb532cef2" }, "kernelspec": { "display_name": "Python 3.9.12 ('xcdat_dev')", "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.9.13" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }