{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Big Data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Dask\n", "* A [package](https://dask.org/) which facilitates chunking of data to fit in memory and parallelism to your code\n", "* Setup is as simple as enable it and forget\n", "* It integrates with popular packages like Numpy, Pandas and Xarray, so no additional syntaxes to learn.\n", "* Can be [deployed](https://docs.dask.org/en/latest/setup.html) on simple laptops, HPCs and even on cloud computers." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Exercise\n", "* Open OSCAR surface ocean currents (large dataset with 10 GB size) and chunk it\n", "* compute monthly mean currents using all cores of your cpu\n", "* use 3D ocean data (EN4) to draw temperature and salinity vertical profile\n", "\n", "The data is too large to be distributed through Github. You can dowload the data from link below:\n", "* OSCAR: https://podaac.jpl.nasa.gov/dataset/OSCAR_L4_OC_third-deg\n", "* EN4: https://www.metoffice.gov.uk/hadobs/en4/" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import numpy as np\n", "import cmocean\n", "import cartopy as cr\n", "import cartopy.crs as ccrs\n", "import matplotlib.pyplot as plt\n", "import cartopy.feature as cfeature\n", "from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "\n", "
\n", "

Client

\n", "\n", "
\n", "

Cluster

\n", "
    \n", "
  • Workers: 3
  • \n", "
  • Cores: 6
  • \n", "
  • Memory: 25.25 GB
  • \n", "
\n", "
" ], "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from dask.distributed import Client\n", "from dask.diagnostics import ProgressBar\n", "client = Client() # your local cluster\n", "client # click the dashboard link to visualize the task and worker status" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "ds = xr.open_dataset('/home/lijo/WORK/oceanCurrValidation/data/oscar_data/oscar_vel2003_2011.nc',chunks={'latitude':100,'latitude':1000,'time':10})" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "11.97881186" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds.nbytes/1e9" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Dimensions: (depth: 1, latitude: 481, longitude: 1201, time: 648)\n", "Coordinates:\n", " * longitude (longitude) float64 20.0 20.33 20.67 21.0 ... 419.3 419.7 420.0\n", " * latitude (latitude) float64 80.0 79.67 79.33 79.0 ... -79.33 -79.67 -80.0\n", " * depth (depth) float32 15.0\n", " * time (time) datetime64[ns] 2003-01-01 2003-01-06 ... 2011-12-26\n", "Data variables:\n", " u (time, depth, latitude, longitude) float64 dask.array\n", " v (time, depth, latitude, longitude) float64 dask.array\n", " um (time, depth, latitude, longitude) float64 dask.array\n", " vm (time, depth, latitude, longitude) float64 dask.array\n", "Attributes:\n", " CDI: Climate Data Interface version ?? (http://mpimet.mpg.de/cdi)\n", " Conventions: CF-1.4\n", " history: Wed Oct 18 12:18:41 2017: cdo mergetime oscar_vel2003.nc ...\n", " source: Kathleen Dohan, ESR (kdohan@esr.org)\n", " VARIABLE: Ocean Surface Currents\n", " DATATYPE: 1/72 YEAR Interval\n", " DATASUBTYPE: unfiltered\n", " GEORANGE: 20 to 420 -80 to 80\n", " PERIOD: Jan.01,2003 to Dec.26,2003\n", " year: 2003\n", " description: OSCAR Third Degree Sea Surface Velocity\n", " CREATION_DATE: 16:06 04-Apr-2017\n", " version: 2017.0\n", " contact: Kathleen Dohan (kdohan@esr.org)\n", " company: Earth & Space Research, Seattle, WA\n", " reference: Bonjean F. and G.S.E. Lagerloef, 2002 ,Diagnostic model a...\n", " note1: Maximum Mask velocity is the geostrophic component at all...\n", " note2: Longitude extends from 20 E to 420 E to avoid a break in ...\n", " CDO: Climate Data Operators version 1.7.2 (http://mpimet.mpg.d..." ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "array(['2003-01-01T00:00:00.000000000', '2003-01-06T00:00:00.000000000',\n", " '2003-01-11T00:00:00.000000000', ..., '2011-12-16T00:00:00.000000000',\n", " '2011-12-21T00:00:00.000000000', '2011-12-26T00:00:00.000000000'],\n", " dtype='datetime64[ns]')\n", "Coordinates:\n", " * time (time) datetime64[ns] 2003-01-01 2003-01-06 ... 2011-12-26\n", "Attributes:\n", " standard_name: time\n", " long_name: Day since 1992-10-05 00:00:00\n", " axis: T" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds.time" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "ds_mean = ds.groupby('time.month').mean().compute()" ] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 37, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeUAAAGoCAYAAAB49v3QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9e3hdxXnv/321JVkySLYRlmRkc+TYOMEtdoAESJ0LSYASfgFSQtoSt01SWpoLJ/RycnLhnFxLDy1tWlLScNwkNDl1mqQJTdyUNEBaIHExIS63YEhiYzWxkC0QwhK2ZG9J8/tjrVmaPXtm1sy67aWt+TyPH2uvvdas2XvPmnfed94LMcbg8Xg8Ho+n8bQ0ugMej8fj8XgCvFD2eDwej6ckeKHs8Xg8Hk9J8ELZ4/F4PJ6S4IWyx+PxeDwlwQtlj8fj8XhKghfKixgiWk5E7xZen09E38rpXt1ENExEt+TRvsfDKWpcE9GpRHQnET1BRHuIaDDre3gWH14oL26WA3h37FnZ8AkA9xZ0L8/ipqhx/UUANzHGTgdwDoDRAu7paXK8UF4AENEgET1JRJ8loh8R0XYiuoCIdhLRT4nonPC8k4joG0T0KBHtIqJN4fGPEtHniegeInqKiN4bNn0jgHVE9DAR3RQeO5GIvhbebzsRUQb9PxtAH4A707blaR4W8rgmoo0AWhljdwEAY+wFxtjRNG16PAAAxpj/V/J/AAYBzAA4A8FCajeAzwMgAJcD+EZ43l8D+Ej49+sAPBz+/VEA/wFgCYCTAYwBaAvb/ZFwn/MBHAawOrzP/QBeqejP+wA8rPj3KcW5LQDuAbAGwNsB3NLo79P/K8e/BT6u3wTgWwBuB/AQgJsAVBr9nfp/C/9fq0Fee8rFfsbYYwBARI8D+C5jjBHRYwgmIQB4JYA3AwBj7N+IqIeIloXv/Qtj7BiAY0Q0ikBzVfEDxtiB8D4Ph21/XzyBMXYTgknIhncDuIMx9vMMlG5P87FQx3UrgFcBOBPAzwB8BcGi83OW13s8SrxQXjgcE/6eE17PYf53VEk9ntxcvH4W+t8+9jwieh+ArYpr72OMvVc69goArwodb04E0E5ELzDGPqC5v2dxsVDH9QEADzHGngqv/QaA8+CFsiclXig3F/chmFQ+QUTnA3iWMTZh0FAnAXS53sRFo2CMRZMcEb0dwMu8QPY4UrpxDeBBACuIaCVj7BkEZvUfut7T45HxQrm5+CiA24joUQBHAbzNdDJjbCx0qvkRgG8D+Jf8u+jxOPNRlGxcM8Zmieh/APhu6DS2G8DfZn0fz+KDGPOlGz0ej8fjKQM+JMrj8Xg8npLghbLH4/F4PCWhoUKZiDqI6AdE9AgRPU5EHwuPn0REd4UJBO4iohWN7KfH4/F4PEXQ0D3l0EHiBMbYC0TUhiBu8DoAVwB4jjF2IxF9AMAKxtj7G9ZRj8fj8XgKoKHe1yxYEbwQvmwL/zEE2XzOD49/AUFGKKNQPvnkk9ng4GAe3fR4rNm9e/ezjLGVWbbpx7an0eQxruO4tnMNu796GP85O9nCFpFHcsNDooiogiCcYD2ATzPGHiCiPsbYCAAwxkaIqFdz7TUArgGAU089FT/8oQ8T9DQWIvqvjNrxY9tTGrIa1y7cVX0Ox9kcfrdjYA7qBDJNScMdvRhjs4yxlyLIS3sOEf2iw7XbGGMvY4y9bOXKQhdxHk+u+LHtWcy8u3MN60AL3rKkD989/hzf6lwUNFwocxhjzyMwU18M4BARrQKA8H9fEs3j8XgWCd+tPocL2k/CypZ2LG9p5dryoqDR3tcriWh5+HcngAsAPAlgB+az9rwNwDcb00OPx+PxFAnXkvtblgAAXt920qLSlhu9p7wKwBfCfeUWAF9ljH2LiO4H8FUiuhpBBZa3NLKTHo/H4ymG71afwxvbT45eR9py++LYW2609/WjCEqfycfHALy++B55PB6Pp1HIWjLn9W0n4evHRkFE1Oye2KXZU/Z4PB7P4obvJcsspr1lL5Q9Ho/H03B0WjJnsewte6Hs8Xg8noaj05I5i0Vb9kLZ48mI6ampRnfB48mcIsZ1nJbMWQzashfKHk9GtA8/0ugueDyZU8S4jtOSOYtBW/ZC2ePJiOqBfY3ugseTOXmPa1stmdPs2nLTCGU2cxzV0SFUR4eM501PTdWZY6qjQ9Fxb4L0JOHYPdtzbPwI5vbuMp6iGr/iMT+uPUmIm0+zwFZL5jS7ttzo5CGZQbPV6G8+kGa7+rTn10xShvM8njjm9u5C2+p1uWkUrNKG4wObgakpVCYPAagf2/w4AFQnw2Pha9Nz4PHoqI4OoTJxEJXV63K7h6uWzGnmuOWm0ZRlxIlInLBkKpOHon8ejyvTU1OY7e4HALTlNXm1zK+d+bh2Ha+VyUOFaD2e5mB4/Ej0d57ma1ctmdPM2nLTaMqs0lZ3TJy4RA1DN6Hx49VJoK13MPtOepqGSMB19c0LyomDud5THrfacTxxMFooyOdVR4f82PYYefTpw+jpbA3myhzHNBH90qbKic5aModry81G0whlrk3wCSlO8OrgQttPXh4VfMLqkRZ3s4Jwzpy5GWsP2Nnu/hqBrMKPbY8KPrZPbx0HqgCqwXhqW53bLc/cUFma+OKVLe2YRtMpyk0klEPiJiSgXpPQwfedOzo7U/fLs/DhJr2ejpZIIC/UbQ8vmD2caKHZ2Yre6jO5W3xE0rpPN6P7dVMJZZNpuuY8g0CuTB6q24/25mzP9NQUejpa0FsdB7gjVYGTlzxmVfe2WWiKeMHsmZ6awumt46iM1Y4nvo+cm5+ER0vzCOW5GQD2gpljcpyR9+I4fiJbXOjCiWa7+yPh6CoQXRAjC8R7Z4Ef14uX6ugQKphf4FUP7FMK4eMDmwvu2eKmeYRyyNj0XN1+nysL1STpSQ43TY9NzQSmvI4gMKEyeUgZMscXc7Ze/mlQOTHaotuqsd3C8Sxs+LgeWHECHn36MDasaI/8E6KQuXBxKQrkttXrAv+EojvsaSKh3NKKsWn3TX/ZXK2Dn7MY9peHx4+gt/pM9HoxaVB1Ahm1nvt8HCjHWtvKYjqJWouQbMrmwpYfVwlg/noxjeu5vbsw0nNG9Psuhs/MGR4/gg0r2usWjtE46O6vszL6+PbG0DxC2QHVYNOZveVzp6emmvJhFk20okBuZvhn5gI2cuCqtxYDmF/AiQJ5bGoGQ88H7Qwuz2dcqMzXxtj7UAhbOT2Gn6lZxzUQ/M5j03MYW3o6MDWDsXCob2riz8upWWBW9X4QKsuP+Fx4iqN5hPLczLwwCecwcQDKmoGIq9lRN4GV3Vs7mpwkEy0QJlHR7S0tAE1Zt+8r/ha6c2TNmGOrKQwu74wEcx6w48cwt2dn9JvUxCBn4GwWF59f9nGtS0HKtT+Zns6FNe2J45b/VuLvJCb6kJ9pQD0PxlGz8ExggfQkp6Gjk4jWAPgigH4AcwC2McZuJqKTAHwFwCCAIQC/yhgbN7Y1W8Xcnp0A6j0G0whkXVpDmeHxIwtqRTk2NQOgNdIOxQd3IXtc1v2enYMAgol7rOcMAOqJC6g185p+b3792PRcNMH3dHal7rsOVj2G6vB8ViU5bjSrveHZrj7MSgvOhZAzWxS+dQJkeg49HS3hv3YA8wvQuTGgZf15RXfXCfn755+zTfF+3baLZltDRrlwWUBzWbNhFMpENBFzPQEYYYxtSHj/GQB/xBj7TyLqArCbiO4C8HYA32WM3UhEHwDwAQDvNzXEjssT17qaQWgSyDYetHJGMFGrqI4OoRfAbEdo2pw+goEVJ8R89OKpTB5CT1cfejra6yYvnTPQQkIWsABqtHyTMAbMizcVRU9c1eF9aBtYF1kysnbUkjXmub27UBG+k7KbuPmYDhacATqtmIf8LCm5UBapTB7CaOi3MCD8DCYLISdurJj2kr3ja7HEacr7GGNnmk4gooeS3pwxNgJgJPx7koieADAA4HIA54enfQHAPYgRypy2gXU1E5ZugtWtJuMQB+jc3l018XyzXX3o6WgppblnemqqxotYJVCUC5eJg5jbu6uUGoXKrMeJfv8wFlfO8ma7tQHUal+N0iDaBuo9Y12wnaj5eXOKVJ1ldfyR+yUKYtGqwV/PdvWhsnFLace2GKZWQa2/C59fhsfnF/6mOUw1tmXrX9Lsh558iBPKb7Zow+acWIhoEMCZAB4A0BcKbDDGRoioN/b69iVYeu5FANw0HtXEZJv0f7a7Hy0bw3sJx8umJVdHh4yVsHTfE/c2rkweioRbmeCCwiQwdL9hnTey5MBVJqgtyA2cdFtBN2mbvLKj18L3WkYtmW8b2S6cIuEspkUtmQVAFMKiA5YsmI1tGOZA2zA+15wPnmwwCmXG2FNxDdicEwcRnQjg6wB+nzE2YVu7moiuAXANAKxZs8YY5G7jWW163+QEJAoHANE+ZiOpjg5Fe+wAUNm4Jfrbdm9d9sicLYnDT3V0KOiTwmNUhi8k2noHazSQsu+jyWObnXsFjgvvOyXISWHmHgv3ZMuCGE8OIHKwO3tVsKcv95V/T3xVrxsrZTHNc+uPKUJktqsvWvh3dHZiDrDeqtMhbuF4QdxYrBy9iOgKAH+KYGxT+I8xxrrTdoCI2hAI5O2MsdvDw4eIaFWoJa8CoCwFwhjbBmAbAJx91ll1NTVtB5fJ9d8krOQqVKp9zEYwt3cX5g7si/Yg21avAwStqBkeOpV2rPJM5cQJZPF6FUWbbuPGNsfkzJPkc8r0Vp+ZT8HYYDOvyelsfh85dLzTLCTKbIa3oa13EGIqmeroEJBCIMf50/gkM8Vj6339ZwAuZYw9keXNKVCJPwfgCcbYJ4W3dgB4G4Abw/+/adtmEoFjK4xN74nhRDiwD0vO3+rcjyyY27sLRx+4E4f3DgMAlmHeQQioT6OX5IErg1aR1BnFZP0w7buqvPCThlC5MsNqNVYbfwhTX4zlSw1tNnL/lW/BzD+rfOoKwtGGnp/C4PLOKNzPhE4wl2Fcq0izkEjsO7PAnDxNENHFAG5GsEX/WcbYjdL7ywD8PYBTEQysP2eM3VZ4R0NshfKhrAVyyBYAvwngMSJ6ODz2IQTC+KtEdDWAnwF4i01jSXJep21HppHhRLJAFuGe6fMCOuircl8xnLRNHplc8yzLPrNovRD7Nrd3V03IjGlB5YKqHdvyikmwEcgujlnaZDkWhS+KZm7vriglJKcn/Gw9He11QlgVd26r8ZVtXIuIUR98XIvoFoxxv2GsZ3YJxkBSiKgC4NMALgRwAMCDRLSDMbZHOO09APYwxi4lopUAfkxE2xljxxVN5k5cSNQV4Z8/JKKvAPgGgGP8fcHcnAjG2Pehr771eqfG5mbizwlxMeslrTxV9C7csXu2R4J32foB47lyaI2I7YpcjpcsEm7C4xNojZNW28oo9ro6OoRZwc/AxnEvjkaY/lspxrSuMz1mEObS6HEtIgoH8e9e+fNLCdAWklDp6OysGdfy/5XJQ0qBDMQLYvl5jxvvTWK2PgfAXu77RERfRhDdIwplBqArtNyeCOA5BOG6DSFOU75U+PsogIuE1wxAKqGcKS12Sr/LnovK87GM8IxGctgMMB+PyQW2KbRmoe21idqP6JE6Nj0XCWebVIEuk49tQoZMUSw4nUOiDBYCl7aKLvfIKxnVHJPGNIAa64+OiuRXYfJMbsRi0xbT/q8I/55ska+P0rWWfF4QnSJDtoU+GUAQYvtz4b0DAM6VmrgFwZbp0wC6APwaY6xhoRhx3tfvAAAi2sIY2ym+R0Rb1Fc1DlV4TFxGrri9Y1VOWJc28oZryHI42PHwu2jZOF8BRnxIa1bMFiEStpnNGgH/nWoybWUcLx5X3g7IV0C7OtzELR5krbOMWpHKbF0NHRj5Fs3E0Ai6B1cBe4exbP1A3faMqTiHbEmoGfsliKBw2Uuu8WlRoBqzNvHNedPf0YoXL21PfH3HMQKbnXeKVKCyxMqOk78M4GEArwOwDsBdRPQ9xlhc8qxcsN1T/msAZ1kcaxxzM8pBXIQwVd2jyJW2KJAjU25YwlKEC2iOzYOnM3uWQTi3rD8PLQi1KeG3l7Vim5AeWwcw2fyniinNGl6QwlZ4qibbGmEU835ZkHN8ywJnYmik5n/OMsxrzmJxDlPVrDJZwsRMge3Dj9SEeqrGqU2KXPG3N40D8T5l+k5ScADAGuH1agQascg7ANzIGGMA9hLRfgAvAfCDYrpYS9ye8isA/BKAlUT0h8Jb3UDdItYD4eEuYKX9wvaPR5NPjUCWEB9k1cOWlYMc59GnDwMANp2yzLrdNHDHF47Lvng0sTl4I7um4ywaUQCJx1R/l5Fj92xHy8Z6Q5xs7RHpHlwFIPCnELdoRFTfS905BvN10YU52noHcWzPTlRiTMhxn0v328fuKYfPR1mT6ljyIIDTiGgtgGEAvw7grdI5P0Pgw/Q9IuoD8GIAqfNvJCVOU25HsPHdisDWzpkAcGVenUpDnMmnqNVfXntv3MP653f/EGsueFlgrhNWx6JWqHqYeJpBIPl3YdIMuUDmYSqXbSxWAET7pDXhM/F7iPwcE6r9dx6fPtvdX8rJKysBLGpYeYzt6uiQegHEt6QEwdM2sA7LoHZoVO0tu1gBdONE3NvWVdPKG9Nz52L1MoX3qZIijY03xAk5ExhjM0R0LYDvIFAkP88Ye5yI3hm+fyuATwD4OyJ6DIG5+/2MsWcb1ee4PeV7iej7AM5gjH2soD6lJssEAWJburjUos08bQPrsOaC+b+Dfc76fqn2VW0zNNnswcvVhCqTh3B6KzDatjLXqkkquCn72J6d0XeRxe+jS8Mqmw5nu/rQk/guZlilzVm42k7Srl7YeTpC8Xu0CAu5mn3f7v7ot5UdGcVjNvcwIWvCPG98IxIE8XwHbQDmhg9a5XCPXiue4THFtpbqXPEeG1Yk3/MtA4yxOwDcIR27Vfj7adQ6MTeU2D1lxthsWEqx6XF1BitaSIsFMAC9x6m46lUl5E+K/Lm4xsQFMhBMEj0DffMxuyuKTzYhL1LqTLkx8btAvAOgzbFGkqSGbqPgIUA6E6zJ5CoK8FnUksS7nPeHa8LcZD02PYfYBPwFo/OW1pHEklO2cb0YsHX0eoiIdgD4RwBRRe20ccpFkSZ9oinW0/VeSeACT3Z0WXruRXV7SfJErBLOWSOHrMx29+eaRCOOJedvxbF7tgNQT1oq4kzYpsVa2QUeoBdOScZqlj4TXODpQp44bavVWwfG/gHKXPi2oWHRQiF83dPVB1QFbV7bg8Yi57AWvyubOcAL4cZjK5RPAjCGwGWcU6445RBbbTUrAWuKcc7SxCea7vhr07m8D0AyRy7XxYiuD42YvJacv7XG8UtG1ihcFm1yCcCiWUiJMGyoC2G0dMRKeq+6YzHadJmc+VrWn4dZMYe7hTe9VeKjhClZPflgJZR5vHKpaWmN9TIW3ysC0ayb1DGko7MT0+jTamU2+0hxuFgLTNdHmnnJtcckoWBArfmvjE5dRZKFJzIf2zKq/Oxpn9msBHIZkouIST1080Jl4mDqLQwezaHag/bkh22VqNUI4pK3INCQvw/gOsbYgRz7lhoXz1oR2z3kIonb7856oaFa1NjuuZdhdd2y/rxIW7adkHSfL04A8wpFcYUQssAmxlRHGWORZUz1za2uV3gPaxfncTmfpetr2mpQcpG23kHMicLWco7TCWhTXoeaDHmewrCdRW4D8CXMF4b4jfDYhXl0Kg/kmNQiPKezalsWvC4m1DRmfHlCU05OCe7ZKNL85qr9uEZOVlYxt1LyDfF/1TkuZCWcdNskWS02045F1fgvKzULCIvxwVEtqr0gbhy2QnmlVMrq74jo9/PoUFG4CivXBzLLB7girYx5qJPKo9o6xMWwV2p7vckpKq4f3BkLQG5lLuW9dR4O5mIFANTfb63zXL4aMs1WnYVp3baGwyRdJDoTctF7uXGZsmwXL2LN57ySjFhbfhSLsjJnMvME2M4mzxLRbwD4h/D1VQgcv8qDJs1mVpTBJCsiF2DglZFUmL6TpHHMvF2dpsP/ntXsPbJzr4juba5plRyVKU71ni06q0Fenu1xuApa+XzXkn1K03lG8brywjNVWwU8p1XB4Ur0F6m5d04mbt1cZHRYdLSINGpMe+yF8m8jqKTxlwj2lP8jPFZ6sspa5aKBqkxvsxkVUBcXHnJyEFcNUJfxS3V9R2dnjRagu6dqYTQ2PYcB6aPfuf8w9o8fBQBct+VF2j6mocYxa2oGPZ2tNceyDBHJc0GovadFXK/Le7bncuE+292P2QyzeyXZ885CALf1DkaLR+19hL6ZFuhze3bOO6rllGAkT+cr1RgukzKyGLD1vv4ZgMty7ksmODl0WUxqjTRbA+GK3OCYYSxJ6PjgJl0dywsX2/jItSuWAsjGS11FJIin1KVRXawCIjqv3LwmL3b8WF11qiQCuWxMT03BNldUnB9I2u9DuRWTwuRfPbAPbQhipSOL0IoTErWlQ7ZSJXVs5dd7yoGt9/VKAL8LYFC8hjG2ILRlV7LyUs1ae7Jpy+V+cQJJfLBVoSvyPrRuH17VpzecPD3/YnI6l/SF3Bta9IrmgjopKpO47vNnjU1KyaI8rMUwpaLCg/Ja9IimaBnRKmDTj5aNW4CwVGpev0Nv9RmgmkvTSrzALhbb2embAL4H4G6gLptdqUmy0s3yYUojmE2TRV64THy2NYsrk4dq9temp6a0JcZ4GFPL+mzSc8oFOlwFcly8slzDOW94veD6rFfxeZ/LhGgd4RjN7gktEUlTbdr2S2R4/EhdKs7K5CH0Iuj/dFZbWIq9/jTzjBe65cJ2hlrKGHt/rj0pAUniWW3Cq5KGj2QxoVjdJ4UGEmc+5ylCuQYsp+WMzpVzGifuUcD01FSd8xufHNPQiFARVj2G6vC8EOZ/60oUFkWasKhoXBRE1s+SapHQ09FSo8HKKXBnu/tTO38du2d7XQlLnoZUnIvSbqfYhj+m5YTepejpW5H4+spk81UQthXK3yKiS8JqG5lCRJ8H8EYAo4yxXwyPnQTgKwjM5UMAfpUxNm5sqEX9UdIkjnBNUZd17HNc7t4sVrgujmHi+Tb5h6MJors/ErK25lVuJXDdY66ODtXGpBuS9if5/nSLkIEVJygd4fJAFsZ1BUliwozK4LhjM37irql5FhxqCifplziebCpmqTzVKxMHI5mdxHfihe0fB1C7KOO/fZoFjkrLFh08q5OJm/YkwFYoXwfgQ0R0DMFakAAwxlh3Bn34OwSe3V8Ujn0AwHcZYzcS0QfC12ZNfa7WmcdUnL7mWA6mG93kJ5uj4x5MV0/hJOi8xU3nu7YPzE9cYnpAG+KcwMSsXbNdfYApGYomRjNNeBS/tjo6ZLx3WnRasc5srav9nBRdG6KgESfyOEQP4hrvfcOizWSVykrrNoX4WUdgKKpaySFfo+NHIgfETacsq2tjePwIeqvPRNfJ2xXa/jv6FETWLOm7rRnXnkKx9b42Fsclol9gjD2epAOMsfuIaFA6fDmA88O/vwDgHsQIZVWCBRl5wMoZsjguDlByW6pjSSbDR58+PN+fMJRHl/giCaJQyWph4mJZiPscKguESSNRXV+ZtKvi5LofJ/abC6DqZPEaqEogmwRaFshjOok5eHg8KDQnjmnepk3qy7pjlgtwbZu6eGzNuKi7n2CS1jo7Sv0J5phW7B6ZxNDzU7hso+L9cLUjl2vlfgWiM5nSVB6D/CyqxjWv4ewpjqxSEf0/AGdl1BYA9DHGRgCAMTZCRMpSpkR0DYBrAGBNX09d2EhaoondIruXSevSTfgmr1XZIUlcKGQhRBvt3KHTlADUpfmL9oal/bi5vbtqwsVkXLKy2QpmnfUja69neWyLuAhjIJ/FgilxjInq6FDN3r68D8qP5YWNRm36LeMSqMRZt+TvaXB5vVVB1pI5Uf10g0BOivzdV0eHtKlZPfmSlVCmjNpxgjG2DcA2ADj7xS9iwPyqMm7i0k0iuuQZtvBJJq3pedXYYwCAkZ4zonZ5f2xSYuockvL0EFbFK0coHFxUHsu673u2q7YUphi/rUpWAqT/rC777VlrE6qxraqe5EKavNdZhfcdu2t+l6r1ivfVtR2XGEM2dWvPk8OYChAslclDQNvKKOwubnHR09GCno52/GT8eHRsemoqCHlSIP7+XnttXrISyiyjdjiHiGhVqCWvAjBqe6HNHptMHoJKDJVRodL8gED7qx7YFzlzDGwNQoOqo/UPahJPYFlgqQSjTbs2Xtf8b2DeKqDKDKaK/dXdh1/LfS6z0vhtterZrr7A+UWsa5uz1YHal0Tj2joZRg5CyMbSMDY9B0wfqUuUMbd3F44+cCd+fvcPo2Onbw2ciCLHNM12ElA/DvIWsmI9YlmTF5H9HLhZuqejBbMdemuCeHzTKUEb3H/CxmGt5hlLUVRE/kxtvYORluwFf2PIv9ZcMnYAeBuAG8P/vxl3gThxydhoyDpUGm/dvpqhTZXwSmKe4w9LVsi1gV0XJuI1qolHnmDFNJsdnZ3A9JG6vsglEHV9kicTU/+TeqvHnStubWSlmRvvl3dYXEbVo3TfAbdgdQ+uqjkulyKsyekuJXqRPep1zk9JPJJN3tuq8TY2E4TxbJI+S2XvLvR292O2Qy3I4xY2uj5wc7XJkmR7HyuHTouUop58yEooH48/RQ0R/QMCp66TiegAgI8gEMZfJaKrAfwM8yUjY4lLQGCLarVuyuAkC4a4gTzb1VeXD5ubZFs29mPJxi01pj5gPkevKmlHmpq+8kQqJ9xwRc45DQADilDEmvvECGSugcnlKvm9VAIyj8mkMnkoym7G76fKq50lrJIuZ5ZpS0dHpCk6hBmJi1QxUUZ1dCi698nhvuiz994XXdey/jzMjg5pTeQ1z5UkkLnTUxbInzeptihq/jZ7zLOC9cfWLyZusal1OIsJ+eLhT/NWj/zybHvU2KbZ3ALgYcbYkbBa1FkAbmaM/RcAMMYSp19ijF2leev1Tu1U2qw14qyyL6kqLAHxGpPYt8gcq3hoTtz6YeU9Ve0nTR0Z19e6jFia87We7Jp+iXG9PR0tGJvS90cn7Gq+72rt9Vzr0EYNSggAACAASURBVE2q/LOYJk2lGVX6nHwxkTZ9Z57YOIapBHBSE7E8wfP2WzbOa3oDUrlO0QN4Xgi0RosrlTCOwyYtqYyNEOZ7wSpa1p8XhenxzyIuHFWIFgB5xHEhLXqEq663sezoULUpW9I8xWE7i3wGwGYi2gzgfwL4HIK44tfk1bEk2GjCYrlDQC+UdCkk4/aKxb64aGoq7U+urGTTvzxxFcic3uozqI4+U7f/Jgp5sXCE7PijurepFKNttiiTNSTumqjfHe2FmK9d0Wm5LhWkXJNxxN1vtG0lYDG2+XPTCwDVoPLSHNSCOC6Rig7Rc1mV39qE+Fwfe2AngNqa4IEgnb8HH9e9itm2ffiR6Jq5muvXRZ8ni7rIytAugyMj9/f3Arl4bIXyDGOMEdHlCDTkzxHR2/LsmDOajF42k6xJMKfBJcRGJdiGx+sdZlwyR4nVkZJqcS6xx7r9wGifWaqtPDY1U6dxcO9bLphtNXl+X1UoidgHEzYTkGqfM4rXK7BIgA157UO7xgCL589bFg7XJM3gmaPE307Uig/vHY6OL1sfVOCOM1vbOn2aajnLkQSRPwkPF1q9DkcfuBO4ZzuAWuHM4QtOeaFemThY49QpfiZRQzbFhidNH6p7ruXa7N50XTy2M/UkEX0QwG8CeBURVVBccZhMyUqbUWmqadoWhVAW5vU4gZy0/TiNIuq/sDcsm4pVWhI3CYrfKf9OTPtndeZnQTCa0mnaWjz4udEktSLcRy6oAEWZSOzlGwqOnoE+2Ew5oqAShXEcNlqyymwf/a3ZmzU6TXX316S7FIk+d1cfxqZqF5wqgWxCZ/1J40Oj+2w1/cwoFM5jj61Q/jUAbwXw24yxg0R0KoCb8uuWOzNMrfVmNXGqTD22bZv2LVUDnpdmU703PTWl9U4VifNgBpJ5gZsSFig15um5QHsdO4je7n78zU/nojrKYhYjHg6yKhT6cR7uxpAsg5eqjM6pTdU+d7Qbm5rB0PPcWtGBs1d1FS6g48zzSWOYVSkibfogml/lggkcOrAPq2rMsfNWoOmpKbQLgurw3mFMDI3guZ+O46TTAi9B2XNbdPJyCYXUheCNhfHCPZ2txudU1Jbn9uzEs/feh2Xr9wVx11NT6BC05RYAc/dsx+kbt2DukZ1g514xHya4Wv9dqfaRTdnrksAFLrdcjY4fwe6RSewfP4q1K5ZicPmK4LtI1LonKbZpNg8S0dcBnBYeehbAP+XWq4S4Toy26TJd3o8jLi5Y9i4WvVh1AtkYDmRpUk06yesEYE2fhD783lmrapIlcDo6O4HOQRzbsxNtq+c1ctvvO2nok67Pcnt8P7yjsxMYPx4JYv5ZitScbbNSAW6C1fY6+RpZqMgOVuLr4O/ANDs3cTByjGoXrhMF8ti+oA4NF8z8fW7C5vewQWdlERdZg8s7QytT67yVRxFlMd/WvqgvfH8YUtlRbtJ+4v3vxfqNWwKPa6HNlo3q/f8as3qKIip17SrSaQLB1tjukclQIM8fX+j7ykR0MYCbETi4f5YxdqPinPMB/BUCC/CzjLGG+UvZel//LoKUfycBWAdgAMCtcPSQLhMmT0ibEANbVDG4Ks9f+fwa86oY0yu0o4oTdsHFu1acrMW/R9tWYmz8eNQfec8smkzC6yoTB7FBUf1K5OgDd2LpuRfVJZWo61PM51ZtLyQRnnKBhU2nLMPw+JHAu7ZzZdR2GUnsQa3ZppC9oNtWz5tuTXHDNYJ5eN5s2yYdVwlkzsTQCLoHV2n3lF1NubVbLcFvzPd/d49MYnB5J4aen4oElGyVmu3qQzv/7MP7os+4RFMLvHtwFeb27ERl45aaPok5v5M6KKo+WxyqAi+XbezHjj3Nk1Yz3Gr9NIALARwA8CAR7WCM7RHOWQ7gbwBczBj7mS6tc1HYmq/fA+AcAA8AAGPsp43uuEyrZaJPLnRd4/gAu0xDqvPF/d35v1vr3gPmH3xVrK+MrTDOMvuR7BTTA0ThTDLcC1UsQVmZOKjVKIBAq6iGJeq0fYhJjJDGuc+WgRUnKLOs5U3RNYh1qSpV2qnOFCtS57QlvCcKZADoWRdoxyedtiIyXS9bP6A0WafNTyAucsem1DmpOeKiumdgM7irIl9oLIm5l/i9mnLAqz6HGBefF5dt7I8K4pR1senAOQD2MsaeAgAi+jKCgkd7hHPeCuB2xtjPAIAxZp1BMg9sf9ljjLHjRIHkI6JWZJ9asxBsM+KkRSWQa97XHNehLFBhaZqOSwJhc72JmtJzQurQY3t2BpN0KJRlz1UdPD5bzGCmy5Ptara2zfxlU++2rXcQAwhSSOaZccumAlpeuN7XRjBzJoZGMDE0AiDQIvnfoqmav6cSxnHfuct4ELd6NqwIruvo7NSWVeTCcWwK2DCwGZXufiDGaWvgI/8XQH05RHFs6zyrRYEMmJMFmSxotmN70ynLgn5OZj83Zo1YvCVkW5g7Hgisuj8X3jsA4FypiQ0A2ojoHgBdCCKMvogGYSuU7yWiDwHoJKILAbwbwD/n163yY7uCjDRfjRC2WfGK58gPXJ45amXPVFcz+ZLzt+JYGCoSmeW7+jC3J4zt1Jj55HuKCSVki4bsEJYE+VqbesBA6Jw0/Eim+30LCV1yDlEwuwhojujQZRLGut/eNgWriHLxochNDwRWEl5+EgB+Mn4cPZ0ro2xlcWSpGMhac9yWlu3Y5gsHG8tTGk7oO6FuEeZC65OtYC/MF29RoLKhygplK4CzEWzHdgK4n4h2McZ+krhjKbAVyh8AcDWAxwD8HoA7AHw2r07lRdokDy7X1WWLUghnk6OWbjJRVclJK5hlTVoXJiLme7Zlyflbo7hqVVyyiZqFh4MDXpLEKvzzt8QsFOT+ic5ojXaI0cXfZlniLwuWrR+oMWHLXtU1TlxCmcLZ8JhpLOg8q3UoM2Tx78mgUQ6sOAEDK4Brv/4oTus/EWtXLMUbXdKYCnHAun7amt1114qfzUY7Vt270WM6Aw4AWCO8Xg3gacU5zzLGjgA4QkT3AdgMoLxCmTE2R0R/D+A+xtiPc+5TZsSlZwTy3TNRt21e1YphRRzbcnXyOVaFzh334pJqg8vu+Avgkj8KNMvQ0SUO0amoIgi/uGxfHJMznXgN9zGI9rot4LWceZYq1+2IrFGG/uRoQUkq4NsG1qE6vK9G8AL6hCC6Qgzib18GC8XH/uFhXPShwGG3w+J8VWUs0bkzLoNdXO4BeuB2VBwri3Gqo0MYawscGGvD/xYkDwI4jYjWAhgG8OsI9pBFvgnglnBbth2BefsvC+2lgK339WUI4pLbAawlopcC+Dhj7LI8O5cUl9VdlglAALX3tmhylduXK99wVNl/8kRONWg70bUPP4I5BB7TgDpfd0dnJ14AsHtkEr88ug8tG/vx7Wc7gGcP1sQqq+DtLpE8VkVM8elcMItauq5Yha2WPD01Fe2Tc4HcyInLdtJN61vgci/T3rIsmEVBrKoZnTSPvWphK6djBdTPVvXAPkAq5KEaH7e8eRNu3vkUhnbvxtj0q4KUspNmzbQyeWj+uxGcIG3yTWvzEhjmLRcteXpqCujqi2K2h56fwv7xo9bXlw3G2AwRXQvgOwhCoj7PGHuciN4Zvn8rY+wJIvpXAI8CmEMQNvWjRvXZ1nz9EQRebPcAAGPsYSIazKdL6cjL3BLn8RuXbF7XhipOlNdxVb0vV/zh6fji4hxtSKJZ8WuWnnuR8bwTt34Y2DNvbgcGsX/8KG7e+RQA4LotL9Jee3jvMHovrD8uJymRj8970wZCk3vYyhqx6+e+c//hGu9cPnH99OALTu2kJU/nsrhUjrPd/Whb7damLp7Yxota1iQ5VvvFMc+gipZwEWjzDF235UX4OMLtmRWBxWX4M/8HwLxzV01/Jg7imFTdSlUlLelcVpk8BGzcAiRYgP1EEMYAGjKus4YxdgeCLVfx2K3S65tQkoRYLrmvD3Pv6zLCM3qlQafFxmFT/SUOUdjGJWeQj9mWe3PBdu+YL0h42JPJdHfZxn5g41bM7d2Fs1d1xWbCOj6wGUvPBdoGgs+qMlerJi+VqW/o+UBTv2htvTet637bnU+M4rT+E2uO/fTgC3hyZMKpnTRkKZBjvZgT3ss0JnWCWGcJEf8H7C1aVrnbFQUpIv8J4ZjpjmN3XI9Hnz6Mn4wfx+nd/XXmeRdkK48NNee1rURv9Rkn/wiOLIyfHJnA+KGFLZQXGrZC+UdE9FYAFSI6DcB7AfxHft0qnqQCWYfN9SZh6xJaomsrraDWmeLlY657ei3rzwPGj2DZHX+BsdDhR6VRAIFgjvaT4+KPJW1D1K7OXtWFr+05BAhCOcmkBQD3P/w0nuwLhPJLVnUDQKGTV1YCOW07tpWi4u5n2i9WkUQgxzlD6tJxuswDm05ZhkefPozRtpXoueJ9xnO5A9uItP2xYUX9olGHLl45KGJzQv0FFtz5RBCiyxeY+x49iIkDC8aNqCmwFcr/HcD1AI4B+BIC+/wf59WprEiyPyw/hK5J4OU2eOJ5EatasJaJ6mVEk5hJqHOBXbeXbVESTpdoRc4+BqCuylXN8a0fBtc3eQiGqLXy8I2z/+RJ3PauFcb9ND6RB5NbJ/7w9scAAJeeeQqu3NiHkx/5Bi76wlfQcdsOZdISF3Z/4mL0XHIDAOChVC25kUSIZl3EQHufmPzMuraDMRPvkGm8p8aHI3odE6VgmzFOjC/WWVdUcc28vjJfBLasPw9L1p+Hm3c+hdcuDRaS4naIeF/unS06En5nw8ujv9/6vtfi5Ne8GoC6QpUrt7x5UzS2PY0hVmqFaco+xhi7njH28vDf/2KMTRfQv0QEFYfSCWSeElJ5nuE9+TwRMb2gC4f3Dkf/4rBtXyewnVJvhjHEWfAH35vAH3xPb/799/1jAGo91sX7y85Wn7zijEz6tRDQjUfT+M2auDAlOVSJ/8sa06LaJJDFPo1Nz1nPBTZUD+zTPm98zPZ0tqKnszWqBiXfb+j5qYXuBe2xJFZTZozNEtHZRXQmC1yqIumSD9g+gKbQIzlHsAwXsLq9J50AlpPxq6hKTiSuJI17FjXWLEPNdn/iYlz79UfRc9Yqg3f7fAICnq/47mtfOX/S+VtxegaaxGLC5resyUaleJ5c44aTYOP/ECeQtW3nGPVw3ZYXYceeg7VVqRRZ+riWvH/8KPaPH8X7XmhI+KynIGzN1w8R0Q4A/wggsk8yxm7PpVcJsUmhGHdcGx4hIO/VikJMVaQd0GuwNkJWdQ0wL9BVBeBFdEI6CwcxMbGHPIlXJg9hGn2RR6fKtMe55c2bjPd5cmSirqasuCjizlwc072yYOyO6wGg4aa+OG1Yl3XLZeGlCzsSHbDkuG/duVlRt5DWPOOqMDAbE71pW6cycRDHwqx0HJPp2PQeN1uLY1uMqBD7WpQX9Ngd1zd8XC9mbIXySQDGALxOOMYA5CaUbcptiegKUuhifW2cVHRarhyWFHetLJBdCrebcGnHVXuOQrNivjPt9dzzfPIQfvEobyPZXm7PJTfgqqsvjfbn5VJ3lYmDQM/KqKpPESzESSvNIkwWvGlyAaQhiQXG5N1t0zc5dn+2qy8qQnE8puKZiQtu+X60zVKZOIijD9ypfUb3jx8tzLt/IY7tZsI2o9c78u6IiE25LR1arVhTN1aFrdezbpIrQiDrSKJ5p4WbD1UT5rG7vojDe4fR8z//KigAb5l7l3Pt1x8FAFx0ei+qB76LZ++9D6skz/S21esC03kB8njtNV8rjTdqnGnVxXvfBpMAi9Omy4CtKX1sagY88acumU5kGZs8hP63h/ndQ+uJDY8+fRhf2Hpm9FqeM7IOcbTh5p1P4Z8fkjNQeorGatlJRJ9S/PsEEV2eU7+icluMseMAeLktI7YCWYfJISM6J8aRqnpgX+TQZeN0VYQAbRtYl2iP2WUhY2Lorj2oTB7C//3PkShZiC3/8Lmg7sng8k5Uh/dhYmgk+o45/Dfr6WzF4PJODC7vjErPZc3vv+0sXHX1pTjzktdg8Oyz8eHrr8KZl7wGZ15SbE10F8cu1XeluyYrx728cfUSd93bnu3ur9GwRcc08frZrj5cdfWlGDz7bJz9v//Vuj89na3Ap/4w2i8WS1qqntWfHnwB44deiELw8uDKjX146I57c2vfY4et+boDwEsQ7CkDwJsBPA7gaiJ6LWPs9zPul025rZqSXaeuPkXZkG6P15WaiY0XaB9YV6Mtx7Wft5ZcRtoG1mFs31cBBA+9Kzsm/hWb33UhDne2Rt8f//5Fk/x8Jan86swCwGvX9kR/33/oBVy5sQ9rVywFAFx27StB9L8yuU/N2D5l3l/BZgtBDMMzLQyzKmgCILFpO+l1SUlyDzFGXrye7/ne//DTWNF3ou5yJavGHsPR9QMYC1+3DawD9g4rk6qIWzJ5mrBXjT2GwbMDn96h3buj42N3XJ/ZuPbEY7tBsx7A6xhjf80Y+2sAFwA4HcCvADDnV0yGTbktMMa2McZexhh72ck9PVEogfiPk5VAVr4vaNguIU95aMl5a95JtOWz37kF3362IwpVm9u7K4rdjOM3Nl6Hz/ze30eT93M/Ha8JD0saZpaU01vHAQCn9Z+IFX0nYtkdf4E3tg/hje1Dmd5HHNsrT1oeHY8Lz7HRlsuCKl+5SzijTRiWGI5lK5B7Olvrr52aURYdaR9+BK946Sl4yapuPHTHvdb7sT9aejrOuHdj9Lpt9Tqc/JpX4/jAZhwf2BwJZLHPzx16IdcENSPbb8NvXbwBv3XxBnSvfnFu9/GYsVUrBhCkiOE2wRMAnBKGSx3LoV825bZqoFlFLAHqhXGWk5PsPMXbjgt3inuv7LhoVsFkE/w99me/j5Nf82qMbH4TgGBQ2UIP3I6JoRGM7RvH5ncN4Od3/xATQyNYc8HLoj4BiOouD7z5k8E9Hfb5bJjt6sPaFcFj8JJV3fi3a7+KN34vj3VpPCqHw5r3hbHu5OQn+AjYCLIkzldx19jeO0k5UUBdQlV7rqIftfebwKc2H8UHt7sbDLfd/1/4y1fNm6TrwsnCe3MN+SNXvRQ79sQXcknCxNAIrtzYh7GpGRz/0ofw4i7uzpbtM+QxYyuU/wzAw0R0DwIt9tUA/oSITgBwdw79sim3VYdJG26EpiA6XZlM1yrBLtecLTvyxCjmDm5bDVy2vh/HXvNqtGzcEk2IAxa1zfdvuxLYdiWO3bMd3YOr0LNuXHlesI0w//qqqy9N9Xl00AO3Y3DDGwEE2vL+r+0AMJrLvUxklftct7iSnfdcTdLaKkc5lEq1Ecg1FZgsymy6fN5b3rwJ1dGhoDwpACgqpclsOmUZ9m+7Ejv2HAQwHT0nYs1onsmLV2k6qe9EZe72LLjpxA34pZf2YW1HS1BN7ScPYuhsuxKraejsWVZXT9uFSntbhr0pB7be158jojsQOGARgA8xxrjmak7ymgBduS3jNcf1CnvWAlkUtlxblrVk1bkcWy3ZNbFII7Rv27SKs6NDaNm4JTgWxi270LZ6HZat34fN4WfsHlyFiaGROg2Qa8x/fsnmyLHMVIHKluroEI7d9UW0DaxDT2crxqZmor3kpDm084aPS52W7LKPrPOsdhGyedYuzxNRiIv5pqPxHabf5M+ly+7y4PJOjLZ1RXWh+WKIC2RViB9PZ6tLYevCzle+CoMXbsQvvbQ28UtP50psGmlY9cJFjW09ZQLwegAvYox9nIhOJaJzGGM/yKtjqnJbxvOrx6yFr0qouQg0+VzTffm5qvZVkyU/ZtOmyUxuG5esSy4RhyqbmcqMODY9BwgF0zF+PNZcqGJk+23R5zwxzJk98Y7Lon7LnvNj22/DlQBWbX0H5vaOphacf/C9CXzwkj8K9jwB9HS0J/oceWCsxmSoVSwjx6TrMnNxTEJW1LBtkvrYmp9VGmwaZ7FY07VGq5arMgEIwv4c2XTKMgyPHxES7AxidPwIdo9M1pz3klXdwKpu/I87fooPvm4dVo09hrmx9AvCV/zdTZjt7seqrYFzYLQNNJAsB7onPbZL178B8AoAV4WvJxHEEXtywkagLls/kKl2nEVcqzjRyibDNIk9dFWkuKCRBc7E0Eim382TIxMYm5qJJs+yECdo21avi/6ZiAt/E3ONu4RNpRHItpq1ay7tqGaxJJBtHM1MQjxpPm+VxstTanLTNRBsl1zziv+GsakZjPScgZGe9PndR3rOCBzauvujGtJ51un2xGO71D+XMXYWET0EAIyxcSJqj7vI01iS5MCOSymaBp5ScF77OBL9bZMW80Rhr274Y78X7UXxSWTJ+efh2D3bUR3eh+7BVanyf8uMH3ohyqnNaaQ51kYYm0g78SZxrioTtlYOvlUhni+XCBUJ3jtSMzZsEubw8T88fqTuPb79wvNk7x6ZxL/vP1rzXlKC/rYCbSsX7PZCs2ErlKthli0GAES0EkB5UvVIqMy6ollYNP+atCkbU7KMrXZmIzDEveqk/dHmvZba0PUnzolI5Ymt8t7lkxtnw4r2ILZ4RV8i7UL1PcthVrzfs939GBl3nyhFeBzq0PNT2LBi4axH0wjfuNrIC10wiyTRjlXCmY/ruUeE3NgOhVBU5vLAGSxY1PJ+cn8Gnoxn7Yql0YLRJe87L3Ua3BvYsKI5fs+FjK1Q/hSAfwLQS0Q3ALgSQKmjyVUC2XSOCRvh6CIw02pwcn/ijsu4Or7ZCGZALQBq9haFyasyeQhze3aiffU69CTIH6z6DoP7B4sQcVriubFrhP90sG930dplaB9+JMypHZjvVLVyLz3zlBotuUhYpU1Z8MRE3ibILBKO2CBXHnPeMw6vSXKtVfuSwObjOmmlNtmJkDO4vDO6l5zjnReq4Kbuf98/FtUQf/be+wCgxh9DZP/4UeW4Hpuew0Bjhvuix9b7ejsR7Ubg7EUA3sQYeyLXnqXAJJC1GmGM5mhyylIds9VE49BdJwrgGmceS8Eso2qr5n3HsBtZi9KF1lQP7AuS+zs6rLRs3KIUUG1CXmyxz2LlKs4bTgYqw0PavfQXtn8cT732OgDz2chswmnyJC9BqKqmJKMsxFKQtsx/P5cQLfkal9zcpvvIVh+R2a4+tGv8HGzg7eoWgJXJQ+jpXAmu3QLzWrK4QOitPoM5BHNh28C6aMEJAMfu2Y7v9L4eQG2GOqBcucoXK0ahTEQnCS9HAfyD+B5j7Lm8OpYFKs9T3QScWnuV2s9yP1N7TwuBDdSHbamuVaWulFF9d0YNOhSAKuE829GHSjhRjLStBMbt95enp6ZQEe+jEs7CZ9BVBOPHo88QHquODmG2qw9Lz72optatStuvTByM9nFyC41qSeblrasXrj1f3obQZQeTFmiugtk1v7ZKmOqEpuu+aNJ9VJ1gHpueQ293P1o29mMWwGjbSow9fTjSbOOSfpj2uns6WjDbMZ9QpO7eUrhWZeMWLOHOW7zC1dQUKhu3YHCm0/g5gPn97SxCrzz2xD3tuxHsIxOAUwGMh38vB/AzAGtz7Z0D1LZEL6RW1wtnFTbexy4ONrpwI5XGIR6TJ70kXtEmTTru/KzRTdr82PxaPZnwkb8v8TvUaZZR9R/pfR4rWpk8hOMDm9ED1CXS4J+pSHRCT7XwyfJeomA21hFPoTGnubbRzkmqzGCjgtNUD4K92rNXdVm15/J5RBM2d0Kc70/rfKiWvIBpW4keIedGWUL7PAHGX4MxthYAiOhWADvC2GEQ0RsQ5L8uLTrhaTL/iRmhbON3TXWZddepjpsWDiYBbXzP0nRv+74Kk1lbt++om4RdJiSel7gyeQhtq+sXOqJAljON2bRtc04kEMP7NEI86OoDi9+l2FeXtgA47WWrFgi2ixf5vLhFiHyeS8yz7b1cTbmipipe6yL02ocfUT4zvI+mkLxaIS2awFvr+hf1Lae9dk9ybEfLyxlj7+QvGGPfJqJP5NSnVNgIUNVDyFNCzrejvtYG07lJyx9yjIkipDrDrm2kjVM2FkqQJkKThjc9NVUzUcjmM37taNvKKNc1UNt3VVKTJJMPz7SkQix6DwB5JfybYfFJMwCzENGO+ZiiDvLCI7rW4bdOQlJhzo/ZCvW487NGDHlSmYWPPnAnlp5bm0td5+nOhS4XxmJMs/hadBoThXVPZ6txzDTaCrFYsRXKz1JQu+vvEZizfwOIqo6VAmpfYjRXcuImWKBeQNugm/TqzjNoHjWCJKHwThNXbGsyd838pfpcKrOn+H31IvhO79x/GF/bcwj//NDTuPvaV0bvz3b1YWz8eBTGwfenoQrRCisr9QxsVk5CZuetMNGEJmsUXxTw1IdFYppQbVJgJrEa6IR0GnRjPVU4l2HLxFbgm5y9tJm+pOMqLZlrpzv2HIwEJ483Prx3GG0D9VY6+bNECVA6WqLYZS58948fjTyygcA7+7R+VfJPvXs1D+0CgLnhdIqExw1boXwVgI8gCItiAO7DfHavUiCGjQDp9thMAtalXZtzXUNLbDWWtNp6neNaAmFvu8gQJ0nxvMrEQbzh5H7Mrl2F167twaNPHw4mo7aV6JHaECfbmr1Osb2o5rKMLgmE2uSnoiitQicMlIksBMGchYnSZDVIqh3blqFMIqDTmNNFsjTxip7gF61dhrFwr3l4/Ah6q8+g910fxNye+RhnXmiFfxend/fXLiyqwQL2orXz+bKBQDuWhTNn//hRrF2xtC4ZDiDsTU/PRc9Ks8SiLxRsQ6KeA3Bdzn1JxQy11jhYuGBT0SbrgSmbPnXv27YD2E9CNhOcal9cW1FImlhNzlX83LiFgdhG+/AjOJ1fezSYqJacvxWb5AQgnYOYDj/7nft5ldHB6O2LugKvbv5791afQWXiIFYhMBvyZDJtA+uwSkg8gmrgLc4p2qw3M8fqElXIwlmvuc3/ndahxyiYLJ49m7GoGxc2i1fTtTpsHdVkS4n8Xaq+f9X3Le4JKbQi/gAAIABJREFU93S2orf6TE0fjz5wZ10sMefYPdsB7MPM8D6QFFlSmTiI3u5+9Kzow4t/8q3gue1D8A/AbHd3TV1pk3WIJxQZm4JScC80iOhiADcjKG70WcbYjZrzXg5gF4BfY4x9rcAu1hAXEvVRxthH057TCFQOLzrK6uzgWsvW5OHsqiFktYcun2MTD6u6ll9zfGAzxnrOAMbVVXLmM3Ud1r430FmfzlAsUhJ5rUthP3X9KliDSC1UDekhsyDumbPZ4jE5TsYJ5iRjS+6LzW+qsj7YfKcqQSiWOAWA1iv0RfeWnL81FMy18FC8FgR+DWK9vJHttwGojVcONOD6tKFAbXw036tOk7O+0YSZKD8N4EIABwA8SEQ7GGN7FOf9KYLKhA0lbiT9DhFNGN4nBLWOP5pZj1JS86AkKCtXU0QhpwLvMmkWBLaf13bPuwhsJk/dBNw+/AhWhW1UR59RZt8C4uNBA4F+AhBev2z9ncasbCpTaqO+v7QJTGz2PVO177AgtjGHi1sRaR0la+7tEBmgIk11KmA+73SPsHCuTB5CdVKdVQ4IBPOSmHaXCGk9J4b+BN2DqwIrEAKhXQkFM9/e6OloV9ablh3JFijnANjLGHsKAIjoywAuB7BHOu+/A/g6gJcX27164p7GvwUQF2D3txn1pRS4eLRmQZbtuy5I4sJJrJNOODjPRNdoBLO8Dy3HE1cmDkb7bDzRh24CSwI79wosDdNu8vsBtXvrvF+m9KJ5kFc2Mdt9c+d2LZzNdGNWuS3jmGpUR6wZPIE3to2A1ubQ5osYoZ7y7NSUc352G6rD+7A0tAD1AqiMBWNZ9LWQhbTJIawMENE1AK4RDm1jjG0L/x4A8HPhvQMAzpWuHwDwKwBeh7ILZcbYx4rqSFnRmbbLHi7gUoDetK9tI2yTaoyxnvKa9+dDv0KNVhDK3LwnFqMA9JoHh+/j7dhzEBet3YxKd3/kyV2ZOFjXnk3/FwJ5J47QjUOlR7NBQKs05zzRCWabGOekGrToXKXihe0fr6sFrytpyjn9th3R33N7d0V1x1skS5IuLW7wd34FWDr7VqQqsVpZ0oZQAG/TnEKKY0x6/VcA3s8YmyVSnV4sTZXKpa7EmoNgMiGbtm3atL13bonyM1o0JNGC09xLROcYxIVk9cA+VIf3obr949o2q8OBWfrYnp01Zj0dn/q3vfgUgG//6uqoT7NdfWgffqSmfi130FH1Oy9EM2IWzjdlzuQkP2tJfSNkXCxBSe9lDFVTWCX4saiMIoLPXR0dmu/LxEEsPfciLLmw1lowt3eXdXrX4wObgYHNoAdurzkuf7dZfdcl4QCANcLr1QCels55GYAvhwL5ZACXENEMY+wbxXSxlvI+lQnIe5JxKbou/h93nen9pAkvsiTLhzRNSJlsxmxbXZswRfX34b3DAC/lGeb1blu9LpigANADt+PZzW+KJsZPXnEGejpbMQoAPSujPMMD688DBOcwUfOenVqYe25ZmK2ThFw5a48GZ8Y0aUbT5utWXR8nkMX/Te+NTSEsPBHSszJ4L7y0pyfwdq9MHkI1rCEOBA6Lq7a+AyPbb8Oqre8AEGy9PHvvfTj5Na8Oxv65V2A2bLatdxCQxi83m1cntR9lIfEggNOIaC2AYQQ+UG8VT+CZKwGAiP4OwLcaJZCBBgplInoLAgex0wGcwxj7ofDeBwFcDWAWwHsZY7Eeca0aq0NW2rItOrNc4sT3GcaZpkWnZeSxslYVGpAnQa7Fc1PcLFD3d69QTUqsHsUNclUAJz/yDSwLJ7WThVATdu4VUfjI2NThmkIZw0IBjaLqK/OSfVmGqKT1yLYZl2nHrslXwkWwJnUCdb2PDdzqYfotReEtO1vxzFw9nSvRc+4VWBI+e70XBmNfFMhA4H3NBXcljHXmCW+iKIbOweCasCBLe4ZOdY2CMTZDRNci8KquAPg8Y+xxInpn+P6tDe2gAqsnkYg2APgMgD7G2C8S0SYAlzHG/jjFvX8E4AoANZsiRLQRwWrmFwCcAuBuItrAGJutb6IWnSk4a8HcCAHpIpxdzOxJUZm1s5q44j5jzftS0n0VUb+Eus2zYd+5EF8ilYKc7Q6q/IjeqcPjR6QygMWvaeMEcs32TYPLTCbFtFBQLdYSeT7n9IzY9kf3O8rC1/R7Dz0/haHn+XkrpDYGg2s3nB58j5uDvjEg0pL5Z5+emqpJvIMm8JUQCWs23CEdUwpjxtjbi+iTCdtZ5W8BvA+hAGWMPUpEXwKQWCjzesyKjfXLAXyZMXYMwH4i2ovArf3+pPfSYbNyzioTkqntJKFb4vV5EbeYyVJ7cBLECd6XqUmPGR1cWXuO+HdoKkS1vo08vGR1qISuTstVleXTHdPerwEx/C5ae5q+uQrkIr4Hbg3hf4vHAXV40v7xo3V5rwHgzidGa17zVJuvXdtTk6ZThG/rjE3PBVs2nsKxHf1LGWM/kARoXsvwAQRZVTgHwmNW6EKasp7UXdAJ37Sr9Dw0Yfl7sMl2lvYeRVPrWTovdMamZiTzcO3E1ch0g8ocyjbpP4UMYNFE3CAHr7zjo/Mk6Zi1+czyOSbtWJUWUyw4wdk/fhSn9Z9YlwMb4AUqOmsLUoQWp7Ew49hC+m2aDZeCFOsQupIT0ZUARuIuIqK7AahsIdczxr6pu0xxTHZh5+1H8Wlr1qxRnVKqPdk8vLezoAjzX9LvP+vfz9RfOYMR37cLnG7UHrNAkCEsa8SxPbB6fmzrvg9VXV+RPPeMrdopoUDO0oRt4zRnsmiIbZjym6sEuMnkzXNgA/NCeV6rrhfw/P21K5YqM+Z58sf2yXgPgjiwlxDRMID9CCpFGWGMJam5bOPCztuP4tPOPussBrjVVdURtxeXZEIRTY22E53KvJ5Gc81jYeIS+pWGrISzqBWbkNMN8r078b15T9kgpWeWgkYe27ptC5Mnb9akGfdJr7dpN23bmWxRxQjc2D4YPLPlc+IybA09P6U0Z4uVorjJW1WwgnPzzqcMFaY8eWFbkOIpABcQ0QkAWhhjeTrL7wDwJSL6JAJHr9MA/MD2Yl3RdxfinGXSTHpJJxMX83ujspJlpdnn5UDG0f1+Yp1ZkZ7OrpoJUVddJ8/MR3GCuAhcQ6jSCGSXz1gGrTtL0qa1VAlkHbKJW8QksD35EVeQ4g81xwEAjLFPJr0xEf0KgL8GsBLAvxDRw4yxXw7d1b+KIDfpDID32Hhec7IWQCZzUhbkkeKwUab6vDKfidmGstCYXb7v+f7zazrrBHPeiftnws2bMnhT56X5FkGZtrJkxHlGtaesG1uqAhKiyZq/FpGFtk4wew25McQ9VTzv9YsR5ATlOdsuRVBTOTGMsX9CUJ9Z9d4NAG5waW+GQVuSrE7zcXw4ixDM/D6Npkjvbh2quGc5SUQWJTpdvd759fKkqdrXyxJeujHrjF5pcRm3rotPl2fOth+ZL9iz9ncw9j/+9xbHocr5i6MS0qIANl3ryR+r3NdEdCeAs7jZmog+CuAfc++dA2LNWRndhOAS7pHUbOdCo4Wzyju8bFpFkgQlKkEuhzfFXcvR7x125bpwK3OlHlfhKWIMx3IQzEU/M7IHf1JMPiI2pSFNntu62ttxlGHBt5ixHcmnAjguvD4OsXp8iWiW5All0JoB84RjEthZmK1Ne8kuwtl0bmylKsl0btq7zytx/7FZtYZuQuUtXkYySfVZ0LOSRzSEqU2bCls6ZaP29fy4VG4xdcYVAvQUie1o/n8AfkBE/4QgPOlXAHwxt15lRFyptDKTlWC23QNMMuE0qoJWIxPlR0lEAKDa2NhlHaaQqDIKZ5WFqGz9jBvTNmMyba5tuZ2kz6yn3Nh6X99ARN8G8Krw0DsYYw/l1y13Wlso9d5xHLZZv7KaALMwZxdRpCNNTmFXshbINmUAjZp2AQuEJZWWOs/wpGOrzIJPRDZdZ7lATebolx6bqII8S6W6UsYF52LANvf1qQCeheCYRUSnMsZ+llfHsiSJcE6S+lInnF1IGoKVdtJKE9aUSfKFDE3hpvJ83Fy9kGohKxecJdnesMEmi5Xr9Vn1Q3tugjHnKiwbafFRbdtw+LOhqyntyRfbUfovmM+q1QlgLYAfIygaUQpaKb1jUhkSXCT19HZ1otG245gPvIzpNk2TSZ7C2E9g9rikCC2apGPaph54WuSiKUnO0wlk+bwixnNleS/ahMpsrlDbkgx7Uw5szddniK+J6CwAv5dLj1Kgy9tsI6yzNFOphLNLwYu8Q7BskbN1ZVmBy9rrPeHvYqqzm3Vd6IWoUZhSOdqcZ8tC0ejz2HaxGRN5acuRtjtx0FoIe8pBoieGMfafRPTyrDuTFm1e4BwFsin21Tb9ZB6COY/JMMsQqTydxOLqO+smS9tC9knOSYPOCpTGMmM7PtLsYWdmdm6wBSxPXM3eroJ0tru/Tiib2shrEeuxx3ZPWczs1QLgLADP5NKjhMwoS1aYSfqwmiaIJPvWJsEMZGea1t0/72vStCVrvC4C0FVYll3b1X1fRWSq0o3FrBd/cZ9xIRJnMTP5QGSBqxAW8QK5eGyfKDGQbQbBHvPXs+9OOlwnJZ3ptShv4iwSl9jEMia5rpHklfvaZNKOu66RAtu04Myr3rd23CRwlipy6ygLko4THbafTzXObPuSVHjqxnYWhX08ybB9wvYwxmoyeBHRW1CirF48o5dz4vuc6iznkQ2r0ZptUmw8u11NzUWT9UTtStGx9WnGb1zdcNX2UhnGKSdOc7UVkFnEJRfibCX11wvkxmIrwT6IegGsOtYwWltUZZjdyTKbVtZ7sHlQZIyx6+Sb9YTgJxg1ZUimIy/csq48lpSF6NCX1MHMm6rLQVyVqDcAuATAABF9SnirG4EZu1SoEiPY5oyV38sq7lfGxVRouiYP8pwI08QaF4WtRq9zzinLhN3IakhJt1N075dFOAP5aMjidUnJ04rjBXXxxEmepwH8EMBlAHYLxycB/EFencoCk1AtOtwobkKJ06iTmPdsPZtVE3jWiwSb84uaUEQTne09ZSGsm6jynMDiFpm2v63umqxM1TILrZZ3UtLs6cqvbcdlmmuTJDvxFENclahHADxCRNsZY6XTjJMSF+aRpUep7WQRN4G6CGbXNpI4yAGGTGYl0GpsBKdLOci4MKs8mZkzhxbYLPrisBXQWf+2aXIIyAvJRoy7tPuvqv3rpDH2SaMNymz1WYzEma+/yhj7VQAPEVHdzMAY25RbzwqgiLzQSa/RJepI2q7YRlbmzSwWEHlRROKGoieuhVINyVWoJ/W9kK8rctzlFSFguoeOtPeWhbMXyI0l7im/Lvz/jXl3ZKFQ1EOv01zSJCRRtZeULEyfSchyVZ/Ey7YRE5Yq97WISRgVvUBKGoJnGrOmxag8DoswZ6fJcV1m50WxrSIdQD21xJmvR8I/380Ye7/4HhH9KYD3119VfhZabG/WGnfW/Wj0Q5t2hZ9UYy6T2T5rn4QsyXILx9SeeH1ewjntolAnoJOMQa/RNie29rALUS+A36A41nDiYiJd9tfKHk9pg6m/qs9TtgWJjCnpf5YhK7p2GlE/Oi22v7Hq3KT3yrLqmGtbRQhiF3ROVeJi0tVnIU0MtM214nfoHcKKJW5P+V0A3g3gRUT0qPBWF4CdaW5MRDcBuBTAcQD7ENRofj5874MArgYwC+C9jLHvJL1Pmgc0jwQgZcLGHNiIhB8u2kgSb2oTccK4EQKY575W9cdEkrScWS08i07HmidxnspZ7sXaCGfTwrTR2e/KCBFdDOBmABUAn2WM3Si9vxXzCuYLAN4VOjk3hDhN+UsAvg3g/wD4gHB8kjH2XMp73wXgg4yxmdAU/kEA7yeijQB+HUFZyFMA3E1EGxhjszaNZj15xrWTxpxVlnAOHbp+xa2cbSaJuDaKzFvtqjlwinD20eES7qSiMnkIvViYk25Sa0WaZ9UmhMg07rPWhE3RBbp726ASxqNtK62vLxtEVAHwaQTW3gMAHiSiHYyxPcJp+wG8hjE2Hubm2Abg3OJ7GxC3p3wYwGEAVwEAEfUC6ABwIhGdyBj7WdIbM8buFF7uAnBl+PflAL7MGDsGYD8R7QVwDoD7XdovwtSc1V6kqd2iJk1X7dT2s+s0WBvto6wCg4+tRoZJcVwc7uJMqFnRyMWKqR+6c/LIK20a965tqhbvWRSx0Anwsek5YAFu1Wk4B8BexthTAEBEX0YgYyKhzBj7D+H8XQBWF9pDCdsqUZcC+CQCzXUUwH8D8AQCbTYLfhvAV8K/BxB8MZwD4TFVv64BcA0ArFmzptSDKKlTUBHCyWXVzSlCIDXahKZbqMi/XVvvYPB/hvcWx/apq0+xGgdJt1uyNnmqjif17FUtUJNqxa7klWlLlfRDdZ4tsoC23VbS3Y+Po47OTgDAQGeibhWG+KyEbGOMbQv/HgDwc+G9AzBrwVcjsA43DFtHrz8GcB6AuxljZxLRaxFqzyaI6G4Aqrph1zPGvhmecz2ClJ3b+WWK85XZE8IvfhsAnH3WWQmKN6bD5MBhesBc92kbLZxMuJr1kuCyd12khs0nrTyQx7Zrak+VpWihp2G0eQ7yNA8n3YpJa1VKqmzwOsquNZgBYGDFCYnu6UplxUq0rV6X+HpqX1LzrKhOURxTyopQrl0N4JWJO5QBtkK5yhgbI6IWImphjP17uA9shDF2gel9Inobghjo1zPG+Bd1AMAa4bTVCNJ9mpmbydTs6xq3GvfgZbG6Typw4tJNqo67kMaUFqdFqM6XMx/JZLGIkdvOQxt2xTU1YtY+C0kFsPw7cE0ssbDJMXyo0dmtslhUcmEMxAvkPBeWJcFKnhDRJgCfBfAGxthYQX1TYiuUnyeiEwHcB2A7EY0iZUGK0CPu/Qg22I8Kb+0A8CUi4uby0wD8wLV91cRsK7TTmHPTkJXW6ZrNKq8EBGkmT9NiwVY4uewXmhZXXCCXAfGzm7SqPHOZx/XPhriQG95O2m0SVX9czOdls04lfW5V5y4CgQwADwI4jYjWAhhG4ET8VvEEIjoVwO0AfpMx9pPiu1iLrVC+HMA0giIUWwEsA/DxlPe+BcASAHcREQDsYoy9kzH2OBF9FcFG/AyA99h6XgP2+ytZmd2yznyTxV6tPGmb2sxi0tF9B3GfxXTvuH5lJZjjFhFlEsgcV2/xrGLP477zrKw44nFxDGfhsd/MIY6A/b7xIhHICKN7rgXwHQQhUZ8PZcw7w/dvBfBhAD0A/iaURTOMsZc1qs9WQpkxdkR4+YUsbswYW2947wYANzg12NKqHJBZ7nnFPfRZenzn4Ujl+n0kMW2rBLTqvllbFlzai/udyiiIbdB9x0mzu6m0sqS/Wxbj2DR+k2jorvewuSYPzd7l/bjzFoswFmGM3QHgDunYrcLfvwPgd4rul4645CGTUG+KEwDGGOvOpVcZ4eL0kVQDyDvBiM1DXpk4mMiZw0QSBzYZcU8zLzOgi6k+znt3oQpjwO73snUwlLcKsuhLzfsOe54yquch67EVZz0x+ZbIr3ULpawd5Ezfx2IUxAuZuDjlrqI6kiVJzFxpHpa4JBt570slFcg672VX8siz7Upc+IluMi/SWzsvjAu2mM9XlDk6el8QyHHEbTuIDn95L/rEz+VyPxsflLjfz9QX0zULfVwvVoqpBVcgac1baa9xuT5PT08Xh7SswkjyKPnn6rWu2083XRdd3zlo29ViMUQW2O6ri2RhBubtpt0yUll50jhn5imcTZ8rq317WwENxFscSj+uPUqaRyhLE5eMbnWbdHJK6wXt4uGdhYOL7f2TtJ/V4kKV39nGkU41ibv+PlmbE7MmK6uG7rqk410l8HVtVSYPOVl16oRQggUI75NtO0mdD12SAcltJdnHFgVy9cA+APuieN/Z7v4aK1Ajw/g87jSPULYgiZNTmskwC2GVpeYf16ZLuzYmtLRaVxLN21bQqjS80k5eLfWPaVaLiKSmbde2XO6VxHLjomE2Aus5wzLhx2x3f3SumHyDC2TAm68XKotGKCf1xJaFRdoHPKnWm5UwT4vKhKabQGwFs067TeJwo3o/D8/vRpOVw5D4G6VpK6utnjRey1mMedcxy+8tH7O9XmtVsHDeVL3fDGN7sbMohLJJ0CU1iRnPT5jezsarE4iPjU4zObmaz7Pw+tb9Dllrg7Em2LLuvc3p8/RkIYyKFGZZYSU4E0YlxArMBGb1uHtxbPpr+uxiuF/px7VHSdMIZZqtat8zaVlZrbBr2k4hjG0cVpLUOHYhSRtlWaHHOQa57teXAoX5WsRli0V2DjJ5Q7uO47TbNa6aZpyTWdZhgmWm2ZOiLCaaRiiziv2OoG1MoUia5O7yxCe3YdISXYWdjak9awentAI5jcOS6bdz8cReaMRppjaLzSwFsvYejt7TebStI+45yHrryqY/tpY7sV8qy9mCWGx6lDSNUHZBN6GZJniuVeSRqMOVNOEfaUNYklLUHrfpHN1vzt+vjg6VM4GIJrIgC8GsvC7F+DaF6aTZL467XhbQrr4bunsUJdzSLBhVljN5XAMLOznOYmJRCmUgfqIWz4mOW05WrmY0kyYdNxFlof25TGJJhEPaiS3LcKC4e5TWA1uD6xjgY8s2q5brIjSrBWvcAspEUmdK8b1GCWNXS5dMmgW7pxwsWqHMMU1qaUxX1hNIjGk7a3QPrWmv1SamsogQnSyv136GBegUk6c2XLRVKKlHfdp2Xe+T1UIxbv5xMdEnCSvzlI+mE8pJHHlsQmZcTIEuD0TaSS+L/VgVps/rGjJiukcS0oQ4ReFXEwfDpAtAy8YtTT+JFWFtsCGJfwWQzA8kC5ysEAn7ZKsY6K6teS2N67LDOpelmgNdfIkWCs0jlBVVooBs4xeTsCC9fSWSCMG4713nqAKoa+3GTVo2mb8qk4ciocAnrrbV68Drgvo9t2JJkt0rjeCL7pujaTfP51v3HdRn9/LjeiHTPEI5Q/I0xea599ps6Mxxqu8wriiIjJySEABmp6aasqJO3ILG1SKSZpw6h1lJe9p5xFSnsQjJ7ZraySKe2eQxXyOQw+9stquvacd1s9I8Qjkm93VWuOwzZ706z0LrFicIp7hQYXK0deyyXYD0JDTd6dozIU7wi8EZJm68znb1Kc/JwxnKdvzKgidNOKKOJKFHrtfL52YyB0ix5bPd/dH30rIx+P+44nuuTnqNeaHQPEK5IJJ6PDfShJ5FGyphlmY/27ZPtt+1yzaFUTgsQEcvE0m9lbNc4CYVRo0OPVSRVtvNQsjrvpdm2CrzeKGciqz2q4uKFU6Li5XA1YlMF6LjYm2wmZTEPWjZ5L2QTX1JJuTIbJ/hOJZxtcgA5RTGgHk86d5P2m7SdkzPSGnj8D01eKGcAJPjU1KTXhYPN28nby0nywnc1Qzo6g2r0spnu/oWfeajPEz3Sa1IOsqwWHWx4tj013Y/G9B7qpvuo1vslnWh46mnYUKZiD4B4HIAcwBGAbydMfZ0+N4HAVwNYBbAexlj38ninkWEUbjeI8uJLGsBH3cf44RiqeHyCcZ20shKMJvani6xtpxFnLjJATHteMla406yKMx7/KvazitkMsn1Rmcwry2XnkZqyjcxxv43ABDRewF8GMA7iWgjgF8H8AsATgFwNxFtYIzN6puKp2inHteHNC5Bh43DTtp+uGI7kes+i43Xq6sAivN4b7TmlZS48eGCbVxsmRavru2Kr4v4zXNZ/GZQZc6z8GiYUGaMTQgvTwDAwr8vB/BlxtgxAPuJaC+AcwDcn+Q+LrHLSSYk171THbaCxkUgJ+mHqV2X/V2T2cxVq8vSxC22pzIPRvcqk8NXTJWoJNhYHPIQaCqhyfujel98z6Y93TllXYxZ1U12nJe8w9fCpqF7ykR0A4DfAnAYwGvDwwMAdgmnHQiPqa6/BsA14csX2vvWjgF4Np/eNoyT0XyfCWjez/XiLBqRxvaxzqVLf5RFuyWjGcdAM34mIKNx7YknV6FMRHcDUC0Dr2eMfZMxdj2A68M95GsBfAQAKc5nimNgjG0DsE243w8ZYy9L3/Py0IyfCWjuz5VFO+LYbubvqtk+VzN+JiC7ce2JJ1ehzBi7wPLULwH4FwRC+QCANcJ7qwE8nXHXPB6Px+MpHerchAVARKcJLy8D8GT49w4Av05ES4hoLYDTAPyg6P55PB6Px1M0DRPKAG4koh8R0aMALgJwHQAwxh4H8FUAewD8K4D3OHheb4s/ZcHRjJ8J8J+r0W2WgWb8XM34mYAF/LmI6GIi+jER7SWiDyjeJyL6VPj+o0R0ViP6yWmk9/WbDe/dAOCGBG0u2IGjoxk/E+A/V6PbLAPN+Lma8TMBC/dzEVEFwKcBXIhga/RBItrBGNsjnPYGBBbZ0wCcC+Az4f8NoZGassfj8Xg8eXIOgL2MsacYY8cBfBlB2K3I5QC+yAJ2AVhORKuK7ijHp9n0eDweT0OY6+hOl/SkpVUOHwSAbYJmPwDg58J7B1CvBavOGQAwkrxjyfFC2ePxeDwLFjk0VsImxNY6DLcIvPna4/F4PM2KTYhtqcJwvVD2eDweT7PyIIDTiGgtEbUjqKuwQzpnB4DfCr2wzwNwmDHWENM14M3XHo/H42lSGGMzRHQtgO8AqAD4PGPscSJ6Z/j+rQDuAHAJgL0AjgJ4R6P6C3ih7PF4PJ4mhjF2BwLBKx67VfibAXhP0f3S4c3XHo/H4/GUBC+UPR6Px+MpCd587fF4PB5PhhDRh1XHGWMfj7vWC2WPx+PxeLLliPB3B4A3AnjC5kIvlD0ej8fjyRDG2F+Ir4noz1EfiqXE7ykvYohoORG9W3h9PhF9K4f7/GlYEexHRPRrWbfv8RQ4lv+ViJ6X2w7jYB8gop8S0VfCmFiPh7MUwItsTvRCeXGzHMC7Y89KARH9fwDOAvBSBDln30dE3Xne07MoyX0sh9wE4DcVx/8UwF8yxk7ido1wAAAgAElEQVQDMA7g6gL64ikpRPRYWAbyUSJ6HMCPAdxsc60XygsAIhokoieJ6LOhtrmdiC4gop3hyvyc8LyTiOgb4UDYRUSbwuMfJaLPE9E9RPQUEb03bPpGAOuI6GEiuik8diIRfS2833YiUuWFdWEjgHsZYzOMsSMAHgFwcco2PQuUBT6WwRj7LoBJ6TMRgNcB+Fp46AsA3pT2Xp4FzRsBXBr+uwjAKYyxW2wu9HvKC4f1AN6CoBrKgwDeCuCVAC4D8CEEk8DHADzEGHsTEb0OwBcRaKgA8BIArwXQBeDHRPQZAB8A8IuMsZcCgckPwJkAfgFB7tedALYA+L7YESJ6H4Ctij7exxh7r3TsEQAfIaL/n71vD9OrKu/9vXPLTMgkDENmMplgJ+SCpJgAQaGNF1BQ5Ci0iG0tnqKlh9rair14FDmnKJYeeujxFGtr5WgRa7zVao1KFbAiNQWUyD0RSchUM0lmwjBJhmQmc1vnj73XnrXXt6779l1m/54nT+bb395rrb2/d7/veu8fRWDCuQjAzoorSywk1Cst69AN4DBjbCb8zLsMlVigYIz9Z9JrS6FcP9jLGHsSAEJzyHcZY4yIngQwEJ7zSgBvAQDG2L8RUTcRLQu/+xZj7ASAE0Q0AkDXL+2HjLF94TyPhWPHGBlj7DYEZjwrGGP3ENHLAfwHgEMAHgQwY76qRIOjLmnZgJrqMlSivlEK5frBCeHvOeHzHOZ/RxNzEK+fhf63t57nq10wxm4BcEt47ecBPKuZu8TCQN3SsgbPAziZiFpCbbmqXYZK1DdKodxYeAABg/lIaL57njF21OBKG0dgAvSCj3ZBRM0ATmaMjYZ+wY0A7vGds8SCQ83RsmEMRkTfA3AVgC8CuAbA19OMWWLhohTKjYUPAbiTiJ5A0O3kGtPJoaDcTkRPAfhXAN/KYU2tAP49ZKZHAbxd8L2VKKHDh1B7tAwi+ncEPu0lRLQPwLWMse8AeD+ALxLRnwN4FMCn85i/0fDiNMPo5Fzi62ca0ElAQYOMEiVKlChRojgQ0bv/6mN/9/HfePtvJR7jjRe9Ek889mjqqPpaQpkSVaJEiRIlStQIqiqUiaidiH5IRI8T0dNE9OHw+ClEdG+Yt3gvEXVVc50lSpQoUaJEEai2pnwCwGsZY5sQ5CBeSkQXIMg5/G5YHee74ecSJUqUKFGioVFVocwCvBh+bA3/MQBXIKiKA5TVcUqUKFGixAJB1aOvw5SZHQiq/PwtY+xhIupljB0AAMbYASLq0Vx7HYKqQDjppJM2v/SlLy1q2SVKKLFjx47nGWPL045T0naJWkJWdF3CjqoLZcbYLICziehkAF8jorM8rr0DwB0AcN5557FHHnkkp1WWKOEGIkpcXk9ESdslaglZ0XUJO6rtU47AGDsM4H4EzQqGiagPAML/R6q4tBIlSpQoUaIQVDv6enmoIYOIOgBcDOAnCJpB82IBZXWcEiVKlCixIFBt83UfgLtCv3ITgC8zxr5JRA8C+DIRXQvgZwg6ypQoUaJEiRINjaoKZcbYEwjaq8nHRwG8rvgVlShRokSJEtVDzfiUS5QoUaJEiYWOUiiXKFGiRIkSNYJSKJcoUaJEiQUJl5LORHQaEX2PiHaF5aCvz3NNpVAuUaJEiRILFS4lnWcA/Alj7EwAFwB4NxFtyGtBpVAuUSIjzO1+qNpLKFEiczQ4XVtLOjPGDjDGfhz+PQ5gF4D+vBZU7ZSoEiUaBtP79lR7CSVKZI486frF6VmMTswkvn5mjsVK0oa4I6yI5wKnks4cRDSAIGPo4QTLdULDCGU2M4XpkUEAQGvPgPa8yYkJAEB7R0fFMQ7xuxIlXDA9MoimDVtyGZvNzWFyYsJIlyIN8/OmRwYx29lbcbxECVdMjwyiddWaai/DCLEkrQpEdB+AFYqvbvSZh4iWAPhnAO9ljB31WqQHGkYoo2n+VrhwBhBjSiI4E2seH0az/GXHQLZrK9HQ4PSmo7WsINKsaa7oPIdzS5TQgW/qRlrruw8FY+xi3XdENExEfaGWrC3pTEStCATyVsbYV3NaKoAG8ymrmE/z+HDFP/G7EiXSQNRQ86Sn0cm56O+kQla2CJUoocPQ2DGMtC6P0V2DwlrSmYgIwKcB7GKMfTTvBTWOpuwBG/OcHhk0msBLlBgaOwYA6G5vwmxnb7DhO3owl7lmmPl7Gz2L30+Pm907JUpMTkxEft7ujhZ0tzeU7ibjVihKOhPRSgCfYoxdBmALgP8K4Ekieiy87oOMsbvzWFDDCWWdFsEZk8g4Z5eq3AwBSsFcwoTu9iaMTs5hdHIuYlomekqDFoKSMSbVzEvaLmFDd0dLJJgbWVvWlXRmjO0HcFn49w8AUFFrahyhPDdTwaREAc21mRIl0mByYiJiUoVpEHMzaBt6PBL6SjfN0YPWTYF4XWu2KyzRABgaO4bRiRlBO24c8VBPaJinTrPTFVqwSQi7MLBZRaR2iYWLobFjkSDubm8qbpM38SKm9+1B6yo9XavoWRbe/LrZzl5rNHeJhQNxoymCW4NKFIuGEcps6kSUT9e6ak2Ffy/SMhxNjGLUasnASkxOTMQEsggXLTUrcLr20Yp1KOm6hCiQI3P1xAyAQFtucH9yTaJhhDIHz6kTE965kLZpE7IGIgfIROOV/rgFBTnCWhZ4RQhkOVc0yUaAr1vUmMX0wZKuFxZ4sCIH9yMPHp5Ad0dndLx0+xWLxhHKHUvAzr8SB3jgTf+miohYFSMTmazMtESYvms0TI8MYtdMUJd948plVV5NMRAZFNcOouBA6Vw5aDB3oXzSyZjq3xSb2wZdfrIqJXAhxVvw3NuFYiHg9+tqhu7uaEF3R2dMQy7z3ItF4wjlsHiIzFBN2nHEXD0Y0mxnb0MGyYja4GjrcmAm2DEPHp7A5RuKMc0WDX7PIsPqmT4ETNuvFQVy3kxLTIkyCVBb5oHp/Eala0DymbYuRzca23QfE8StywFJILuYpEWTdndH44iJekBDPW2bcE3DPF18zKoSnrUEuewix+jkXKz+bLBbbsHg4SJXlw66aldyOdUYjYTnyRs5GTq64cfzDoZpYZWZBap1uMBWCUxFu7whQdPaC5znKRJ8fWJ0OqcD8bepVwHDXQyi1U/8LU7cvzVyb/BNokyTOkEsZxLwFD/OD4KNeTb3UcIN9UehOszNZG5OlLWSGGOUSnFOjwzOmzlrtEznbGcvfjo2hcHDE9g7dhwAsLprMQBg4ORKZvz61fVnuraVVTXFENjGiDG6hBpIUrjEROhgO4cLsObx4YriItMjg5gNTeftfksuDFP9mwK6fp7T9QEAwOquIxFd840mj5pvPnoQc6jdjQaHuJGe7exF29DjFd9zt10E4W+bMFZ9Hp2ci/zLKr5QIl8YhTIR2YpuE4ADjLH1SSYnotMAfBZBsfA5BN09bieiUwB8CcAAgEEAv8YYGzOONTtvc8xKOJu0E7EAw9zuhzAnBJad2Lkdiy68OtXceUAWyM8efBGruxZXvHiiRtE8PlyzmwwftA097uSXFTdirhpo3hGqnLZFwZyVyVzlwuG0zTW00ck57DgwDuBITboyRLp+9uCLAIB1K5Yozx2dnEO38OxqPbZYF1TI1z0SmuNF8FQmUftNAr6JKdOiioVNU97DGDvHdAIRPZpift48+sdE1AlgBxHdC+AdCBpP30pEH0DQePr9KeYBEDdt6gqN+PiXI7MSgujYWm7d98T+Ixg8HGiLXDvmAlln0hudnEMParcClKmWM/8dfcpKygF/MtIyubTIOqDMtukEwmcxOYfNfZ2hYK49cLoGAmEsW39k+ua+1u72JszWoG/ZpcGJ/E7KwjOtMBU3mmVaVLGwCeW3OIzhco4SYR9L3stynIh48+grAFwYnnYXgPthEcps6kT0t4l5ZR1ZPdvZG2lXrauCY7VmEnti/xEAwOa+TsuZ8+AvdSSoapB56ZBnbAFQqYkUAVXWAOBmmk48p+SbrEUtmdP2wMkdMYuPaqMpCyr+G9Za0JfpN9XxKvm+Yt9J/Yp1z6ZEbcAolBljz9kGcDnHBVLzaKfG02Jz65es9CsOIiKpMK44JpmWagFcQ460BsXLpxMw8jGXftVFQa5CpErhaB4fjtbatPYCjAoVuYDK+/bZsPUAmG23+6eTIkbbq1Zq15NW6HLohHwtB0dt2xm4qWRrj0nAyL8/p4FasQbJGQHyvfD18rX2d50UXhjPORah+/1Mz2mhpMjVIpzeNiK6EsBfIuBFFP5jjLGlWSxCbh4ddMqyQ2xufd7LzmS69JRYsFbGwWDi2LzvaH8mI6cD1yC4QDYxVlP0cSzCuEbuT2W2FtfO16xisqoAFx/mJNJPnoxLpO3NZ29kWaddVfgqdSb7GusSJLphODht69apc1WJTUXQurzqdK0Cp8/oHlqXV2j1/F23bZ5cn0+RGD0+XfF7+mBqtvH83a5b4P8N4M2MsV1ZL0DTPNqp8bQI1tzqZZZMWxqxlneSYiGMiiAuhRCq5XsxQec3U5kip0cGo02FeL0P5EI0haFp/jXN4rdyfU+ax4fRA6B5NLzXKrtlOF0PnNwRY+Rcm3exfPDjsmCupY2HCZFmHEJ813VWDZ/3vciSsSXUcKXE4ZwEsq55tLXxdAWa4sSYhXnPF5z45fJ1ReKJ/UciH9L6rrbonwhRkLk+H7HusynAqmj41Ofl56quqYeNiWqNPpsDHv/gMo+q1SnPB64GeAcjjs19ndE/ANo2g0o3k3BMPL8W6Fo0W4t06vLbySb8JLWrRYHsSi8lsoUtJerK8M9HiOhLAP4FQBRRJWi2SaFsHg1N42kXaPOKVedmvCtsHh+OpVsUDW7GAoIXVPQTZqEJiFpFLfiXZZM1X5vY3SsSJIrf2dTq03SeDE5DRWnPSefRld40YXbpith7Ug19UqRrIG6u5jm1/LMIMZreFhhXS3QNxH8rUVDzoLTpkUGMtC6vCOLi55ngEksgzqsM6CmRG2zm6zcLfx8H8HrhMwOQSihbmkdXNJ42QtFPWW7lKP7Pv1MxuEzM2gXn9g6NHYsxJ7nud1LIjFzcuQPV6cvLzdPTUoZOhTl78hjQ/bLoO1H7s/3GIiOMKkRprilam0hLn6rGFLHxFZ+rZUm4fftzUYoTEAZ1JUjXsaW5cYyEqVLVKjna33USJicmYlXJOH46NgUAGNx7BEEpl+AF4DEjqpQ9ObBNPl4R2KcoiCMGNJbIH7bo63cCABFtYYxtF78joi15LiwL8F0+UMmIxe9k1JNfRTa5yQLZ5T5sVa58i2lUA5yBqfzMFRqUZKKTx1Bdo0LhwtihGpkLVN3QXKwEvMd4UelDPLoaiEdYu8ZFJPl9asW3LLrfurkVKKyyJcIUxCnfi8u9yeU2SxQP10CvvwFwrsOxqoFmp5WCyCSYGwHy7lfVGUsHXdlGVSQzoI5wria4iVFsPwjotQMR8vp9YxBUpsXcnslcpYky665lOsFc7Q5Sq7sWR35jlZBQWccia5iDuV4sMRpDFavYtXd0AB1BRbXmowfRfPQgenhN6wlU5GO7pDOqUA+b7YUIm0/5lwD8MoDlRPTHwldLUdnRrqpgUycwvW+PsXcyMP/S+gqrWsTkxEQsnYlHyfLnIEN33z6CqEJ4KZhX0Q0MeElIFyYsCjMXgSNv9FRNKIrQKLIQmkl95hXn5SywJicm8NOxqZhZtnl8OOre5bMhsdFEbNMpvhsKn3I16PrEzsBAyQsTdXcsj0VZ+wpk38DO5vHhKNWzRDGwacptAJaE54nloI4CuCqvRaWBKJh1yCogx8Qc8ipGIGqFzYjfi1jmkz8HeYOh2nDoTN1JdtC8Mcf0vj3Avj2F1wBvG3o8JkB1fWFVwT+mnGQROoHMfX55IS9t1mUjo6qNnSVUHczWd7VFwljWgLVr1QT0ZaEN8trpqFIFsOajB9Hd3wvOtl2aTaiC37SFYhRWt9nOXozmTNcl4rD5lL9PRD8A8DLG2IcLWlNVkURLLtLEJ75IbUOPR4J4emgPWvvntWOVQObIypcudxNqPnoQzQieYVPBJRn5WgLNYg+aNmzx0pL455iJ3vIsRG1itrO3IvUsK3DXjFiwpOhgrLzHH2ldHmP+PHtAdsWYNpUq6DZUKmHFn6kscKPGHEK088aChLK4qZ3b/ZB3hoePb1iMxxB/67zouoQaVp8yY2w27NpUQoLI9LP28akgB3WJAlmEynStg6hdy6Z/271wjWlyYiLqTQzEayYXXSWJb05UDNwlGMhF+KgKVBSxKVP9Nq7m+iyQl2CWc5AjP+m04SK4WbySxJGItbCHxo5FzSswOVd1ARVZbwy/aybNKMJnXw/5+40G10CvR4loG4B/AhBVxsggT7mmkIUvOWvGJQpikRE3jw+jOexMJWvIgPledL5lF598NIbQqCLKDxb8rZzJ9ne53GV+sPmO05g2a4Fh1cIakoLXL5cFMqDZQAk06eOCaj560DmqPhJmilrSsXKXVQoEC97RyuMyDbt0NFNt7MqI6+rDVSifAmAUwGuFY6nzlLMEtS3y0hA5suq5nLemrCqKMrt0hfIFzWQ+g2BWMQC5EpE7aWWLJVf/GYB49SlrMJenYK6mIPTR/FzTnHyQR366HLSkqiaWBtzs73Ov1W7TqULT2gsATXtYecOu24T45OqXqA6cOCfPV65lsGZ/NpF1Na9qIFFXLIHZJdnIyFC1w+uZPhQeOUl9Uc7QpXupwMtKJtWYi04pSVvwRo6ir3ZKjEkYazeGjlpz7DzHCm71DNff06W9bREuuWojdM1+CcAAgEEAv8YYG9Oc2wzgEQBDjLE35bUm1y5RqxDkJW9BoCH/AMD1jLF9eS0sCWypTkWCF1rgSBqt2d7REVX4SYKYhu1TJ9ny0tYTQ3NNY6pV5sOaW71pO+/fZ1IoZZoU7R0dFWZiZRU+TU65DWk23SoNeXRyLgqyqlbFLxVcyoiarhUhP1vxnhsUHwDwXcbYrUT0gfDz+zXnXg9gF4KU4NzgamO8E8DnMV+D+u3hsUvyWFQRqCfTjU0IqhiIj28oaQ/qWhbMPHeZI40Z0lZARRfRmxdc/f31CFVanryZ0sUIZPFOmyrC1YK/ddGFV0euGTnFj8NkETDRhkzbtWS6zxFXALgw/PsuAPdDIZRDxfS/ALgFwB/L32cJV6G8nDF2p/D5M0T03jwWlBRiRS8VihbCWVcIEoWgzKRc288lqd2s00pE01aP4rgNWWhaLlAxl+72JrQNPa6/RtOXm48nI9aXNwfoqtXp1uMCXx+6VqNKSdc904ei2srKvFmHIi06jdoE7rKIxnAIdnLdfFWrqYWcQhYdD/3H4rvAg+tixVmksRpcQ+boZYwdAICwVbCu/8ZfA/jviNfryAWuQvl5Ino7gC+En9+GIPCrLlCkQM6qPrELbIX5ZR9dFgzd9xpeFEIWvlGhjbEpbFy5zHsdLhhpXY5uzAvNnulD0bPQVT3jSBIYVKRmkbU2bAr+0hUsAdKbccWxdc/PtQ55mutdII7DI8eBeI9jMT0wLxN3VMQE8zEcygJGYVCYWEypu7M3evd0NbOByoCxvHDo2BT2jh1PfP2JmTkQ0XUArhMO38EYu4N/IKL7AKiEwI0ucxDRmwCMMMZ2ENGFiRfrCFeh/NsAPg7g/yLwKf9HeKxmwKZOxD77CuIkqUw+OaJpCvnr0h1MUN2LqUmFSgsQx4jykR3mFQsx8LH7iy+AFJQjbA/ySrvbm9A8erCi6pkIWUi7MiX+nPLULHRNNDiKNqFzTBZc3col11yGyookXtPe0RGL/3BBNdP+0m7+5L7LgF4I17qbCgBCAXyH4fuLdd8R0TAR9YVach+AEcVpWwBcTkSXIWjPtZSIPscYe3vatavgGn39MwCX57GAouEqSGsl6Gdo7FgiRmvMzRWqQ+k0jLT3z5mWaTfOv+O1jgHkojXze5GFsAid1lwLTIlnFti0yqwEclH3PDRWmQssQ47MllHkb2N6X8RzgGBdojYNxDXqNFC9V5FQFQICVTRtMktXm86rhG0ArgFwa/j/1+UTGGM3ALgBAEJN+U/zEsiAe/T1cgD/DUHYeHQNY6xmtGVqW5R6DFNAhFLz1Ozas2ZqfKdva0mog9xDWobLOLymdUUpSoPwj1rtTcxUaBNRXWMAs+29RuGdFHyO5qMHjQKZQ1uW1CGvN2+GphPIWQniIkvFclg3bQUJZG4BskUwu84nvmd5WU/4s5Nru5t+Q937v0CFMcetAL5MRNcC+BnCYGYiWgngU4yxy4pekCsn/DqAfwdwH4DZ/JaTHXS5jkkFpquvWC67mIbgdVpyUvO1i0nfpdCGy/li8Bk3IXMMjR2L8pj578SjKyYngnvLwiTKg7mOP3xP7LhYAQ1Il6tdiCBrir+meZmobbSaBU1zDI0dy3wjpiv4YXteKpp2KWMpr39yYiKq/S7zi+ajBzGH7DpMiQK5Ito61JYryswKLWz5eTpk+VvXMhhjowBepzi+H0CFQGaM3Y8gQjs3uL4VixljutytBQfZxG0S2LHPnpGqPdOHYu3q8mDGeQgVztT4etuGHgdCZjQ5MSEUFonnlosNF9JG9Z64fyuAoBb2z+97BACwdKAPy9bOV+MWhbGOgbloz6ZjjQJVRbkICX6rPCtm6bIRfCqZuWzeVe+jysKWZe2EJ/YfMZYijdYhRKRnXdGuRL5wFcrfJKLLGGN3Z70AIvoHADy67azwmHOVFQ5dgYW8I69Vkap5CjrdZ926dIwozRpVz7hCM9D5sy3mNf57JW0PKArjI7uHcHTwAF54NiCdpQN9Fef7dh3Srb+9owPT497LdcIM8+v2kwdUtJaKhqTgOF+TfFKBqapkpntnbe+YvIFUvaOAFKSXou3l5MREvCTptN4iKLviZH+z9l4k8OJFJYqDq1C+HsAHiegEAt2NADDGWBaVTT6DILL7s8IxnyorMRgbMaTYDape0CL8iKadto0p2oR2Eqbqu+vna2iSPtuC0GY7exNHrKsE8inrurBsbT9a+9dYm3Yk3cjlTQ/VFswy0qRFcfrj/a97oNYyTfO6QhbAprFdx+pub4ryq1ul701IQ9eqCHvXcrKiYHaB+NsukCIiNQPX6GtjwjQR/SJj7OkkC2CMPUBEA9JhpyorrsiCWbpqCqa5xApTgL24gFjirvnoQe8ymTpTetGmqpHW5cYWjrod/OjkHEbHjkSfVZHZfBcv9pYWccq6rshsvfj810fzZY2FrE3we1fl7eoQD8DTdydz8X+a5hD/lgWyThPWvR+2eBGfje6kpQTvifu3onXVGkz1b8oksC/pJkRuq1kif2QVafGPAM7NaCzAvcqKEkUIHaM/JwMzsUu6CIc6sniPViPMyrzuUklJt8uu8MsrBDPvNmViCnJ0Ne8tzX3HonacVgNWPbdoYxUWSVnIcBUYPCBKxHyhC/U1ctxBEnO6SxyAb1UwADGfuq8Fqnl8GNPjlRt0kV5VmRe+c/gi2ihMZh+QV8KMrJ42ZTSO36RCJZeXrFoJoPoCWWYYpvWYzH4qBqeqQiULo4o5HPsju8K3Kpgq2CaW69m6PBb4pRxDwxQmJybQJglyMbL61PDvpg1bopQBmxYkQueuEM/l1o+sG6GItH3aaaflZrY2WU98AqNcwXN351OFzKlqYlUqIJ5jL64r7yA71zTA0ck5oHV5RZUt8R2Mab4KF8DQ2DH0SO+ZqqiPS4ZJkufCTezyvCWKQVZCmWU0DodLlZVYJZfNZ2/0XoNqx13tCNrpkUE0Hz0I2rcH00N70HLl+6LvVKk7ravWBMK6f02FYOafxTJ7LnD1y+tqEystBYI2wQsYcOYk+/tEKDcn4TMSE62i5yBoGHy8EYMv1ua3NAkm+ZlmKZhF2t50zrmZvF8uqXwu39s2NirND5gXyIHlowU9iP9uQNzq07pqTQXNm6w+Li0fTRB/c1v8iEnoqd4D8W8dLYrWsZFQsOvmEK1DWtO/Iy9Tva/8GVSbHy5E1KpdwlplRQWnSEnH42ngYl5z0Txa+9dgUUdHINB6BjC3+yFjBCXXEiNh3K/Pv5U1EBfwNbuWdOQakXjfnFk3734I3WENX1Veqc6sOLf7oQrTJ4CYQJbXGVuP487ftGGTO1DVcsexIunbBNnnPjoxg/6VA0DPAFox775Q0aMxn1YXBOmZhaH6rUWhJG9GuUtl48qB6Hh7Rwe6Id2nQ860yqdsSreih7+KE0N7Yu93kqYcKrT2DHiXHC2RLbISylNJLySiLyAI6jqViPYBuAmaKitJUK2dnsu8ct1gsfhA66pKX3HT2gswKwWLqSALY5XZW6eZALCau3X5nyphyM15qkAvMS3GViS/behxzEG9mdBV65LXaWJ0rnTSfPQg5o4eBJauUKaC5QXfyOss1mIKslIJK1GIiSltskDu7ojHCajo2iZERaHrssFM4sKRA8VGWpdHnwcPB/ekKwsrblzEjkwioo5NYTCjrUEEH1cXC23Nvdb8njF3zDiAzt7IxSRvqkvkD9cym1sAPMYYOxZ2izoXwO2Msf8EAMZY4jI1jLG3ab6qqLLiC19iShPp6YIoLUF4ecSqXTFNeukKLLqw8rHOdvYqGZIJJmalrfkspCa5wkeDEhlx8/gwujsChmcTPHLBD7HOr63QR5JUmuhvhWk05hsHAIFp1yt0rf9EuLo3VBoX/337u+LCLE36nqu1x1Uwa10Xwm+9uU+dkMLzem2uGLmFogxxExbb+IWVwWxwrQegOiZuQOL+/xJFwFVT/gSATUS0CUFPyU8jyCt+TV4LS4K0O7qshLGK6HV+U1uRAxNUAjlNyUhgXtC5Pgse/SyeL0eL9kwfwvTIIaWfsfnowaAVHW9g0d6mvHedRs8/y4EuurVG9+mz2ZAEsvj7FNWy0aVUZNZBjiZXCZDMVKqKMVCNJ2upLmvzeX/ld0znR0nZS3QAACAASURBVBbPESuF8d/c1jvZlkEgd2xypSVuDRPfdxdBrCq1qfqNywCv6sFVKM8wxhgRXYFAQ/40EV2T58LqGS4Ryk5CV1H9p72jQ7lTNpWMBMxBSBXrtQg4XQBL9KJ3xkuCzrYHzI2bMUWT/fS+PWheuiLSMn86NoX1XULlpXBM8Z7lwKCKlouSEFVZQHyFmOyzUwXu5CWcWxxyG3wL22gD3Gy+14w2rjJty5WjXN4ZXyGcBmLREEAI5tLQti2uRKd9qoSh2DNZ1wfctpGO6FdoTiMfA+Y3IryoCz9Woji4CuVxIroBwH8F8CoiakZ+PbwXBMRdeJKoXVlzNL1AvJKPTbPwibbUBXup7mVW8FEBUPZWFjWK0ck59AhjTYeR6KKvXLYSqNZu0/REs7YpCl9+TnIxB5lh1wKy0pzTCOEklit53SKd+aSzAXZhYqN3OdhLPpdbeWI0JlgAXFo8mtbsu8nzcb+ZKoFxF5mcflaiGLgK5V8H8JsAfpsxdpCIXgLgtvyWVV1knQrg69O07XpNAtlU/ce1JJ8Ouqhm0/jiyy1rB609A0AYydw29Dh6lq5Ad1dwD5Ffj5uoDZHkJk0/y5KZzePDUeQtD/TZ3NeZv89tbsZLyMpCpBqaTpJUpVhwlCJyHnAP1tPdc0U0tEXQa10hYUpSG4IuZBF99gwE2nLHABCmOImd0lw3ArE1Cvec1j3FIbsKWnsGsGv/fPW87u6XFbLJHD4yiWcPvpj4+hMzjVcC1LXM5kEi+mcA68JDzwP4Wm6rqgFkIZhl4eo6XmCiCj9I5usXt94c8ydpg44c5lQxLl0aCGCvzqWab3bpCszt3A5gD9pWrcFTi89Et1SOsbVnIIhoDq/vFnbswX0GY4mbEa45t65ag1nVvCrNOWWkdGvPQMRoOXYcGMfAyR2FVD1yEbLKnFaPe/XJ9XXOg1X4MgHEaHtu90NAguhonStFB9lv65KOp4KYIXFk9xCWAWi58n2xjApO33L0uatgVhUMUfVG9wleU0F0IwwensDeseO4aHU3ajdjtrHhGn393xBUFzoFwBoA/QD+HhlESFcTOnNuFsJY9XdMQBvMquJueG73Q1EP1he33mycV/RtujIXFUzCWGxMzxlcLOVDKHoAzD/L6aHABH3W+cBsC38O8zWSm9ZegCaEzFlYg6gx80ITQFxzjvnHHJm06TsVoxYZ7ejEEQycHHzmGrP4XPJG0oA167iKYii2+ARfS4SqYxL/zQG1QKqY3xKgJZ8DVMYAqDaYpvx4+fye6UPhvezBqa95NQCgSRPExmnH5DOvMI1b3l3fxhIcsiVNjO/YuHJZRM+DhycweDg4zmm9RDFw5SLvBvAKAA8DAGPsWd961LUGXeqHC7G7nqcqEeiKWPnMh++JjstmXNVaYkIxZZUjGWJzdS6Qv7d3FACwumtx+AIL7eUg5EQP7akw9akwt3M72oRCIDrBLMNmXnc1x8kCVi7usHHlMjyx/wjWd7WFzGuiaozLxaydtNKYS8Cgc5pRqC3L108DFRqyy2ZSTuGxni9tZOTgvFhgomB1GlXkz4sdopo3bIlM2dG9aejatJnSac/y802zCZPfC1XRkss3rMDt25/D3rHjiecpkQ6uQvkEY2yKKAgDJaIWZF9aszAog5EEpqHNd9VowCrETK2SEHH1C8llM1VpEDbGpEtxykrDctEOZ5euQNOGFVFkoOn+m9ZeAAjPSxQ6qqpPuiAsOQWlu6PFWoCDz7W+q80aZLNx5bKY5lPNov1yQFJm44qCL+Pa3oA66tc2lypi2BexNCeFS0P126ssIWIOr66RhgixEhz/zWw0KT8jE12axvEJGrt+y+m4fftzAIKNdoli4cpJvk9EHwTQQUSXAPh9AN/Ib1n+oNlpJ8YxG1Zi0n0HuDEgl6IdsfOF6GHdtRVVt4SymSo/sivSasiiWU+M8OQ9cK/fcrryusmJCYy0Lo+06dU9r4sVXdC1c1x04dUAgCmNuc+0oehub4qqg6lgYoLiuC7tB0cn56L7ybO9nY62nXKzUwpTl+tdrTGmaHiVQFZV6kpUmcun7ntnL7rDv9s7OtDfVXmOuBkLaK0L3d1BYxXbXK09A/HNnKShm6DacIobBdcWjy60fdWG3mgebtIuUQxchfIHAFwL4EkAvwvgbgCfymtRecHVNG1DEoEs/q+LJBZzEFX+06wiL1UwBRGJDJOvMZYuoSgEAYRMrQO4ZusP8OZzVuKeXSMxM6+K4enAhanORG0qvOCixfJ7tPW4rlhXFRmXTUCpaDlJzXPTNa4umiw0XI4k9y2vQ1VEwxVRfEEHMDpxJPztO7Dse58FACxZay5wGOTityldD7pa2TGBq6Bnm2vGRRDH1lEK5KrBNfp6jog+B+ABxtgzOa+pasjDRKeCts1i/xq9YJZ6I5vMlS6MxlT1ShWlbCy/6CDIvvHofjx69/exbsXb8Lvn9lnPr1hj6/KKgDNgnkG5CGQuuH86NlXhG3dpISlicmKignFVww+XJE6B/61rA2q6XifMkwh7F+S1EbXlsPvgDSPfRWvbGhzYPQQAWOJwTVQJz1E7NkEOrOT3NRU2fPGFHC1eoli4Rl9fjiAvuQ3AaiI6G8DNjLHL81xc1nCJJq0muJlaFswqyHmgaX1/FZHoikpdKqY9vW8PWjHPAFTBI/f9wStx8cd/gKWrzsAf9o7g+Fc/h0WX/Ja2xZ88L783pV9PEM7i94OHJyLTskqLUBX9aO3Sr0XE3O6HAInh7R07nirf0gWqUqOu5yf53nV8m2tG21rQk1Z1rQptJSNdIAeOze3cjjnpGTcpNOCNK5fhxE+Dv48OHsAZN34wljGhQndHC5pH4+v08fnatGJebGfx+fp1a68dGQQ6e0sNuYpwNV/fhCD6+n4AYIw9RkQD+SwpX1TkDmte4rQMKy/o1u/SSIDDxhx116kin5s2bAGOHjRW6wICwbz6uq8AAI7sHkLPJcE8J3ZuBzDvR9aBC19tAX/h++6OliidI3ZfUQ50S+yYb9Dbge6XYXRsKvCTVykQplboU3bNAPouZWmEp0rLd6VbFWy+7aYNWzArfacThYsuvDqW1nVg650A7gQA9N/0SeU1vs9Al3In1+LmKTGmYjsm7JrpAoSYjDIKu3j41L4+wqOv6x1pNWNfs58PVNqyjKKCd1SMQ9fvdt70q/dd7b3jqvCvOzG3czuaHbW91lUAFndFu3fRLy23AZQhB9I0jw8HjGt6/hyVdm/CV3YGjPsbj+7Hm89Z6XVtGqjabcpQxS3o3CWxsR2ZuAutq7IGRHpWBXMBZp+16biPidt1oxr7jkdYG8bl2ihf1+C9OwEA/Tepzx9pXQ50L486T3EaVkX+yxtOGZymow16uHn20ZA57njwPwEA61YERvi8rT8lKuEqlJ8iot8E0ExE6wC8B8B/5Les6sOUD6s6F8hOOIuC2TRfljB1YeLQlen0TcXpv+mTOHH/1qgqkgmtq9Zgqn8TBvceMZ8ogQtuU0Sqb1AXx1/f9WOc0rsELwy/iGf7lkbHf3LgaKLxfCHTpkroughi1fkq4az7TpzjSOhPBYBla/srzjPRc3ReQr90Hv5sTtc+VpQz79yGE/dvxamveXUklHXYcWA8qpwlC1pbg5OK2Iqw7Wkf5t9Xm+VJh299J7DFr9k4/66PDZeCuUi4CuU/BHAjgBMAPg/gOwD+PK9FZQXXClpZQSXIfZmjeJ0pStu0Bt25tiAd+bMcXCZCV8ZyJCxDuePAOICgGIEKiy68OsrZBIDtr3wVAGDLD/49Ota09gLM7X4I9ygEciyAS6Etr+8Kag2L6VF9o7sSaQ8qDO7YAQD41neAU3oDreKFKjAvF8EoHhe/V51f0fhjqFL4y3Qpjy9/Xra2X2ne1uXhK2MXNOZxlT/b5x3Q1iSQal/vmunCeswLRF0086ILr8aJ+7fGjt22ZD0A4H0v/jS6dseBcXzj0f0AgvQjIPQ1C1kQuowCrjFzmp9/17Kh7aP7nsGj+57B0lVnRJ9LFAerUA47Qn2YMfY+BIK5ZiGayFRNCkwBSxxawSUwJJNpWT4/DWIMx+Czi10jtTWUv/OaX7pPXXtAlTahawIvgjekAIBNv3eJ8pymtRfgw//z27jpbWc7rjoAF8ZntowBCDZlxx++x5qu4oK9d1yF7stuARAwrKP7Ug+ZCDqBrPosQ3W+rOG6zOsCPpco+HXjmjR102e5g5iL1uyTBsU3eS5V4RZdeDW2/GBeUz2js63inMs3rMA7//TTAF4Tacvd7U2RW0XU0lUCePBwPv5embYbHUR0CoAvARgAMAjg1xhjY4rzTkaQBnwWgsJZv80YezCPNVmFMmNslog25zF5NeCT9sE1X1lbsPnf0ghkFwZmy3fOA7qWcjKTStpZ5rmLrgcAbDScw03SYjqTvJ5oHRZfcz2BTZ0wfi8K2KODBwAASwf6Ko6JEL8Xx+FabWv/GqXQdxXctnXK0GnqKshrU13H32PVxjSJmbvCRaPJzVfh9c/9SPsdNw3zbADZotfNq36F9DxwckcZGZ0tPgDgu4yxW4noA+Hn9yvOux3AtxljVxFRG4DcIjxdzdePEtE2AP8EIGqTwxj7ai6rSgiTubXINZjykAE3oa3SXHy0mbwgC0BrQwAP5qXDjo9ciif2H4mZrGUzH+/FDMyb9zauXAZgWXBBz0AmWnK14CJMVMJXPPbCs2M4ZV2X9tylA30xwcwhClNZONs0cleItG0SzEYrkaXntgs4LbcNPR58FixsWfOV0btvjDIS2oYex/S+PRDfLh5vwQUzh84lVCIRrgBwYfj3XQgyjGJCmYiWAng1gHcAAGNsCoC+bGBKuKo1pwAYBfBaAG8O/70pr0UBABFdSkTPENHucAdTUxBfeJeXVWQWrf1rYv8SryEDE7kN4r3Zcimbx4ejf0DQWGJu90N4Yv8RPLHfHKi1ceWyUIjqEfSlDf7xOuXifKLQ9q1g5IvRu2/E6N3FenOivGDBeuMqFF94diz2vwoqQc3HPzp4QPm9K/j1rmPYUq208RYZvxOq/P3mowcxPTKIyYkJDI0dw5DUzlNGe0eHNsL/vdeci8HDE5jetwfPf/+BqJsaULmp6O5oKaTGetF0XWX0MsYOAED4v6rR0ukADgG4k4geJaJPEVFuDMa1otc781qACqEf+28BXAJgH4AfEdE2xpg5pFGC3PxcjPb1iWyWzWWqnbxYyEDWiLW7/tC8Jp7vq4HIaxHXb4sgd/FTq4o+iGbhHoGCKpoK8FzP8PTbtz+nrZOtw+TEBN706R246+pz8NOxKZx1fFew9iioJ1zHUv5758+0uM+tGhBLr5oEkKwhixA1Zvkzv47T4bK1/fj5fY/ExhU1ald6XTrQl0qoq2B7BjrorA7GVp/CO/DNqQEMzHRgPeZrrfuUjAWCd4Gn1A2c3IHnbvsSnv32c9j8ri3hMw6ebyuAtlVr0BN2TPMpMuKLatD1wSOTaEqRtTA5PQsiug5Ba2GOOxhjd/APRHQfAJXm5Lr7aAFwLoA/ZIw9TES3IzBz/8+Ey7ZOZgURfUxx+AiARxhjX892SQCCQiW7GWPPhfN/EYGZQSuUWXOQRZjExGQyC6oEmyjAVPO1rlqD40K7RdM8co1rzmRUZmr5WFqfstV3p2l+YfPXjrQuR3d/fEM0eHgC33h0P77x6H7c9wev9FrnXVefg9GJGZzZMjbvK5Q2I2LJwv/1b8Ez/PhbTB7qZJAZ18Dmzfils4Nc5S98urgeLSphJAq+tEIwjVlaF/mt8mPzKHCdW0amUfEd4hvavC1G/B3ghTV4gNX6rja0DT2O7u6XYXJiwivf/bPf/ikGd+zAS/veHERfD/ThmfGfoPvenThlXdf85ig8v3VV2CVqYsZ7A5AEf3bj23DR6qA1x0W/8/H8J0yBUADfYfj+Yt13RDRMRH2MsQNE1AdgRHHaPgD7GGMPh5+/gkAo5wJXtaIdwEsR+JQB4C0AngZwLRFdxBh7b8br6gfwc+HzPgDnyyeJO6SXrFqp7TuqSt8xCeIs2tWZmIk8h06j1aWmJBXGPvnULsEw3R0tUV/Z2HG5800YCX3X1ee4LjVC/1s+iqWrzsDt730lzmo7GPncWvvXxJ5b04Yt0TW88EEeWLrqjFhU6rY/flVUTGT07htB9D8ymUek7dN6uy1nC+tTCD4Zopas+pwESWIdRG3cBRVd1DxqCahgamxRUcKzsxfdYf/w+b7hAYINagsweczZbfJLZ6/E4I4duO6XfgEAYnnNsmXj1PB9bx4fRneYk5wHRNq++ZYv4ObweJZ0XYPYBuAaALeG/1comYyxg0T0cyI6I+z98DoYFMS0cBXKawG8ljE2AwBE9AkA9yAwLz+Zw7pUpcMq+jeLO6TNZ2/Mtb+zTchyuJTPi20ehB7BoinbNr8rbNq57VwRYt6kzbclBoCt72pDd0cL+kafxP85+7fwy2cHx8WcZBPWbFwRMMHj89p704YVmOrsBSnuI8/Sl//lDevxkwMroqjZr+wcjjSKLBGj7TNOT0TbeZiMXcHplWt6qtxl3TUuEdhFQRbaAd13xOj/rOO7gOPAU4vPxOjEEWtsBAA8+Nh+LF11Bs46vguzLSuwY8+8ID5lXRdeeHYMLzw7hqUDffNusFWILFB54Ad/+UZcszWgbZ6HvwBwK4AvE9G1AH4G4K0AQEQrAXyKMXZZeN4fAtgaRl4/ByA3l66rUO5HUD+RR+ucBGBlmC5lztdIhn0AThM+rwKwP4uBfWrOygLTZ47WVW5FOEQk2fmrOkklgeu1phaJgLnC1/GH78F7vnlTInPjmS1jeOaWvwAAnHbxeQD0G5WBkzsik1vWQSsf23Qc3zxzbdSAYuoNl+I74Xcbw+IQecP3+XEmz/8G5jVqX6HNr3M1N3OYhLDLcfkd4u9wEppPW/mLC+Tm8eGIRxx/+B58byCIEXIRyhxPLT4TZ2IM6y49Hfj2cxVWC9mM3bx0hbWRS1Isu/v/4D2vvQ57x45j6v++R/imcQO/GGOjCDRf+fh+AJcJnx8DcF4Ra3IVyv8bwGNEdD8CLfbVAP4ijEC7L4d1/QjAOiJaDWAIwG8A+E3jFU0tziXxTOZpneBMIvyMJnIhCE3cKJiu0QWo6OpRmzYfaYS5TTCr0N91Ek70r0HThi1Y7DEfF6q73jnfkEzVX5o3D+hub8Jojmmcxx++B5sv+xNs7uvEjq7FwA8fwDOveHV+E3pA5cvl2rLI7EUTtyycZaF7ZPdQNIZJINsErK56lwuS0HeeiNwzYaGPprUXYNHSFfjr/xEExLkEM+74yKUAgKGxYzhx92eNUfFF4LYl6/Gmt74Umy/7kyAXuoboeqHBNfr600R0N4IALALwwXAnAQDvy3pRjLEZIvoDBOU8mwH8A2Ps6TRjypHYgFtQmElgGovYS0zDNpcLk3EVnmJDeZtg1l1rg6lIiLw5ah4fDnb3q9ZgVpjDp8yIHDFsexYDmzfjheEXsfq6rwiNMJJj+ytfFc2/9pJDmO3sjaqWXV6QhqyCTmMVy2jKQlX+np9jgkogm7IKOMQgRtN5InSbRVOVPhNMG1C5ZaNYpU7MhefgdC/ykcmJCYy2Lo/Vi3bFqY//Cx75xL3B2Gv0vv0ju4ewDMCiDVuCgLP9bmZyG25bsj6qOPbCs2NYO30I3V296O5oqSpdL2S4Rl8TAhX/dMbYzUT0EiJ6BWPsh3ktjDF2N4C7sx43SQMFlWZdwSDE/sNCCpZuDar1iPP4aPPGdelq+3psGIBkLQ6j8cN7lJuum7ruyBi4ZEP0N2+VV/EceH/kjpaoLOfAyR14IgMGJpp9o2fXml/QTRaQBTM/Jn7PoaubLWrMvgJZ/FyRwWDJePCBag7dOnygq2KnwujETOKOYbLJOnDPPKI8d27ndpw+tActV77PO+JbhTM627D5XVuwLty01QttNzJczdd/B2AOQfGQmwGMA/hnAC/PaV1VhUrAAm5C1nVc43mKQDCbVu6LoqueqXIsdb2Xdej+738d/d209gJMTUxYN1hZFVvgfmyO5qMHcw26McElD96UWqeK6tedJ39vMjtnKWhFuFqk8uiexqGzDIl0nSTgb9GFVwN3fSnaNPXf9Em8uPVmfRnU8G9e231jSqG87tLTceprXo1la+c3NNWk7RLuQvl8xti5RPQoADDGxsIotLqESVvWNV0AFPVvDeNnibQCNJF1QGHuV43heq+jYTpJUiHZf9MnMTlR6Sy2/T5Ja3HLWHz+6zMZJwmobZHyuFFAWhqZ6KL9K4rROGrEaZCavqvgX1bRfVLaPvPObRj68O9WHFcVZjmyewinvubV+N7eUQB+QWUqnH7Nr3vHeZTIF65UNB1W2WIAQETLAeRXWqYAJBGcSr90gnFchJvrpsEkbGUzueqarDcQUcqUpmGE2Hauu30iWo9LNKloqpseGUQz4ppKf88AZicm0N3eFBXu520fL9+QTqOI/PSdvWgberxwS4NLGpst110+X2n21VSH84Xq+ajcMlk9xyTjqCxQSd00vEnK6ORc4GMWaN4ld7n/pk8CAE7cvxWt/WuwTPHdi1tvxpHdQ2hdtQa/vzSbzWbThi2Y7eytqPldonpwFcofA/A1AD1EdAuAqwA0bDa5DT4vrUpIJR1PdV6azUUWEJlYTPhOzmnveXRiBuu72kAPB/1M5oCohaMPZJP40NixaM6gTnBnVHls286DeO2j80V/llz9Z8nmmpxDT5UantiEpKk4TCxSPYO0orQuE9frVRthl+981yPCVzCLdN43+iSm9+1BX6wann+ZZHFzJPdn5pjt7MWu/Ufwvb2juGpDb7QOHz/zSOtyQHiHqkXbJebhGn29lYh2IAj2IgC/whjblevKUsBFELqea2pTKENmDnKgiOp6VTBJd3uTlino1iMf90lZ8o2ktq1LO4/U5SmqzpUiCEd1r83jw+gJz+lun/eyPHfR9ZFGw5sI9EwfwtzO7ZHZlp1/ZcTU5nY/hL8ZDnJPY/5CXlWJa+jplHArWHOrc7cinlevgikI0DfwD3B37diEp0kARu+EqS61hwBNKsBdMdW/CTNhiV1elz3JhlP3G57aH5Tw5e6U9V29WN/VF9Nw5xD8hifu/Wx0bNElvxWMG67lxa0349/OuS5WmQwA0LE8eqfypusSahiFctgAmmMEwBfE7xhjL+S1sKRw1tYkRm46NwlUjEKeQxfdqVqLSujpBKFvDrHPxkP2T5uK5FuLjISlMac6ewEhaMtmyp4Uzp1v3ahvK6n6reV1LArNeHx8ICjUcNFicwnK7vam6Py0kbA26MrIAvaIf5+xrecmtNj4bDRdadhnA14EmseH0Xb+6zG7dAVmEdaBF8zZNlO2Sx1+lZ9f3rQ1Hz0Yi4OYRfD8Z8N357mLrseAZh7Of4qi6xJx2DTlHQj8yATgJQDGwr9PRlCSbHWuq/PADDMLLt3Lzs2b3R0tSmGoetnTMgKZObkIXLEBRJKAEtcgqyzN7bZ5oj6xk3PomT7kfK0cfCZqzfL6eIBZ9Fn4vfn3PP2jG/ENEa93rL2HsOZx7gIhLIyjypvlUG2WgOw1Q5ff3nWTl+VG03VO16DHpL5lMc8ZCGiKW24CuJmydelduuqAol9YnB9QbE5bl6Pbko8Y84mXMrlQGLk0Y2w1ABDR3wPYFuYOg4jeCEDbeaMWYRJy/DMXzBwmgcz/TsqQXTQBXScmWbDYrosEkIdg163PV5OxIRKA7e4R4nM7twPQFD8xCGSOwcMTGDw8/3nehCcIagui36BADc0mJFw1WNfgQPFcFy3X93sdfYv0mmgDqphTZ6HiSLJ5kedRbQr5Z+fMjVDYRqZvxIupqKwltt/AZPo3gf8+RXSlKjEPV4p/OWPsXfwDY+xfiegjOa0pd5jaDsbOsxG7xJCVaRKSEHM2yzmsURayNiHuMn5SYa09X7MR0PnXRyfn0M3Nxxpz9vPffwAA0Hd12D1HiqKNnrV0j4OH51OqePs9/vfqrsUxQQ3MC2tx3fL9xJhtR3yd1YLN2qESEjqhm1QYV5zv+M7J5yal04oxHTfQWWcjyOPOCml9KrMwp2VTmpf2N9E+Y/2z0sWf+PxeJbKFK2U/T0Hvrs8hMGe/HcBobqvKACqi0gmuJLtxHw3Jh4ElfRmyeol0WrZpDleNXdR8dH5EVRnDobFjGJ2YiXIyeclIXS1weV3ybzV4eAKruxYrBTP/m2Pg5A7jZqdHuPXpkUHtWtJA55oB9MF+NuQlfGTkxdxdx43Rr0GDTpJRIY9hA6dpLngn0TtfBMSQbxyLnFe4JbgbSGUR459l3qd6D8VrxefGgyJLFANXafQ2ADchSItiAB4Ij9UsVALYZvb1QRrzrcoslyXzEjVCAJURlh4waZt8bNPa+fl8DSbGoDIxNo8Ph8dbMDR2DH2jT+KMGz+IA1vvjM3DGV3z+HDMXNc8Pozm0UDjWC+V+dzcN58yJd6XKLAHD0/Enp98PyOty6saZJSVGyHpXFltcm2Q6Y7DRtumTSZgfn4295V4zDUzg2u//P/1XYFAHRo7FsRVKALuYlHywjjiplSkeZ5xULF5EHiii2ugCLo+dnQyaoOaBDPTdV0uQwnXlKgXAFyf81oyRxJfat7I2yykYl6yYHG9XiWMXMbWrUEeUxl0BbVw7m5vQvP4MI6H6Sa8oALH3O6Hor+bx4cjpje9bw/mEBRJ4Axtfdd8oYQ+BExv/ep4lLbteYn3XaRg9EFSjc42TnTcYEEwuVNUMG0kdbTHv/PZdCo3hboYDEuQaGxcQyyKqlRv89GD0T+eG8xTmKIc+p6BGF3z64D5jnFiMw0OHnPRpwgK6+7SZ0sAldalrLNSSthhS4n6EGPsQ2nPqTZcBXGWmrTvWvLUnAE35iUzPxMzTHKeeK4snGUzo8w8e6YPRWkeXDCLu3JPRAAAIABJREFUaFp7AZowb0KWC2jM7dyOpjAFi2M+gCb4Xiz4IDMvvh75GcYqlOVENzNzTHnc253gqRWmha9gFpGEpgA/q5BrIGXSFEOTaZz79UXNWUXbnK6BgLbFQDBTGhsvryrXR29btQY9S1f4Nduo0U1no8LGRX6HiI4avicEvY4/lNmKagBpGGxWgl23i7fBpGGoGJ14Pi9NWQRkRireZ+UaOtDdNc9EFl3yW5geGVTmM/NjJ0JtQfVdK1BRR/v57z8QFeXn+Z3dQhqJbtMUaRZVtsCkhS44zEWrMsElJiFLuGrOtncB8MxWcPDXyugO05h4T/XZzl60XLlJe35rz4AybkF+D0Ta53WzlyFsMLJvD5o2rFCmEirXWGrJhcP2Vv0/AJ0O5zQETIES8jkqyFGjWTDqpMIZUGsNNqGbRDCLQVJJYJpPNBF3C9rF7NIVWsEMhN13DGjv6ABi58z7qAPTYFzDqcxfV/vt8oLPZi+NnzfP+zEJaFdztQtMmnNa37TqOSaNXZGrlXHankOgIasQ0XvPgLYf+aILr0bUwuT7D0SNLbjWvFhh9i4q8K+EHbY85Q8XtZC0mJlj1qAOV8jXiYIxzW4/aaSsak0yVIxAtVYXLYKfY2Ng/HubQDaNp/pOt0YumGP5mWPHnAr+uyLqrrRvD9oAdIfBYbo806Q5oElRhIulCJhM20njH1ziKVTn+Zq8XbV/V34kCucoJcqw4UwKuRKYaD43oRTYxaK+32wLdNpq0gIYabRW3RrSIBZcorrPlH5qG7NKE9UtX6saS6mVKAJqpkcGYwEwHGLbQlsDCh44Nj0yiJHW5egbfRJAZUctWTjn2XGL48TsnNHcb6PLrAW4L+1q07l8c43F6nsTM040pEIauhXXkea6io1/+Eyn+jdFwpKboUU6Bub7XHMha7IKiQGRUYQ3L8cpVR/jMNF4ifzR0EJZhyJMjiaBkpWfxjVlQ15TtQsDuDJFm9tA1DAixiKkPbUgbHyxczte3HozHv/EvdF3m37vEgDAv/3pl7H5XVtw6mtejee//wD6rn4nTt23HQhNfHzc9o6OWOEHUYvhfr68mNei5qaYhidrf7bfM2/tOmlJSkBfvEJ5rlh8pmBfddZQBZPJgWFNGwIT8yIpQBEITNDT+/ZgemgPnnvn5Xjh2TGcsq4LLzw7hk2/dwke/8S9GLhkA059zasBIKJtIBD8HK09AzG65seA/PLuS5jRMEK5pYmUx7OMinV5+V3rS+cVQOHC5NJExaZBWg0FUJsEVeU0o9+hdTm6z78Si8aHcZ6it/AbvhBoHc9//wEcHTyAo7f8BZYO9GHZ2j1RwNdsZy8mJyYCP3RYtYszrNnOXjQDTmbAtHD199uKSMjHbOPI8KVfHU1GbTY9gstqFUk3vNqCOpqGNvMnLEd3/yYsGh/GGee/PrIQnY6ArrnYfV7wKR/YeieWre3HoktWCK1IK90/vB90D0pUA07SiojWA/gEgF7G2FlEtBHA5YyxP086MRG9FUHU9pkAXsEYe0T47gYA1yJobvIexth3ks4DZKMpZG0azLO7jYvvulqCOSl8SojKx0cnAKAL3Zt+pfLE7pcBm4Jn1iNpulMSYxTN5M2IaxxFwbax8YkSzgouWnKWQhzQu26qQdOyBp9mDfPX2qvkRXTdvTygYxGXvSxG0z2XVJaj5c9zUipr24xwUzCd+DZKpICrRPl/AN4H4JMAwBh7gog+DyCxUAbwFIAr+ZgcRLQBQZrVLwJYCeA+IlrPGJs1DcZzOU0vaBqtOS/Tn0v5xDQwac4mBlIN5u6DJBG6cm3rSrRHfw2c3AGMTUX9lzlUuaGtPQOJeubmgaJ+l7yLSiQdOy8XTZ6bWNW4trlE+ud0XZnSBQBCN4nFXcBMcC9Cd3BtHfSsA81KuMFV0ixmjP2QKGYiTkWhjLFdACCNCQBXAPgiY+wEgL1EtBvAKwA8aBvTFDVdb8iS6RnbT2awSfF9vi4C1TWVJQl0Y9+zayQ6tm7FkjCqvAP8NeF5y5HlhffIzbm1XUsTOaXiFI1a62UsIivtOWkKWt4wWUxEvsdpm5eNFWla9mFHm/eC6LqEGj4NKdYgqHsNIroKwIGc1tQPQKwtty88VgEiug7AdQDQv+o0AGoTVzWFs8+curq8eZu45fm8x3LUIpJWacoaqqpiMp49OF+Pd/Awj3gOro0zOHvj+iSQaVvld60VF0S1SzHa3DVp3/+0z9i3FGgaVAbDxeMQ5putBPQsPxvx3ShbNlYHrkL53QDuAPBSIhoCsBdBpygjiOg+AKpacDcyxr6uu0xxTFlnkDF2R7gubDrn3IpzVKbbPMshykgS9AHohTNHnkFiqvmcrnUUELYcaBvEDk5p4TIWP4f/Lwpr/vflG/TlDpNCpO3N51bSdhHweVdcNpBporRl6ObxsQpZm1WkrEvA4ZN3nUZ4B+9UZREWkXbXrViCvWPHte1K+fl7x44bO1c1CojoFABfAjAAYBDArzHGxhTn/RGA30Egi54E8E7G2GQea3JtSPEcgIuJ6CQATYyxccfrLk6wpn0AThM+rwKw3/ViW0oQf9FqufiCbX2m4vdZIG0eN4epalPSymG+EIWoK9atWKIc4ycHgoqzL+1bGvv+9u3PRQyvaOg2Q1nVUvetDmYSztUqQqENHKsSD7A12HCBrgKaaMmRwelT1aJUhW07D6au1lcH+ACA7zLGbiWiD4Sf3y+eQET9AN4DYANjbIKIvowg7ukzeSzI1pDijzXHAQCMsY/msKZtAD5PRB9FEOi1DsAPfQex5rhmnLuZVTUxcTzX9CoVshDWujF880pNlZt8BLPcA1mFJEKYQxSqIuNat2IJnj34Il7atxQ/OXAUPzlwNCaY08zpArmfcpK620noMVUUcY35m303mnm4BbJ0yajGMr0bKsFqE85ZWqVqGFcAuDD8+y4A90MSyiFaAHQQ0TSAxfBQFH1he1N53eszALwcgcAEgDcj6KmcGET0qwD+BsByAN8ioscYY29gjD0d7kR2Iggme7ct8hpQl9nkMAm4NMLU9uJWM1XDlSnqCl6YtBqfgg82+JqzZeYiMw5XbZULUt35LuPyY3ysvISzjrZr3eIDFO9vVtGzroa5a7ONWvDbJ4X4vuhq4auEs7gBrnXhLMZfhLgjdP+4oJcxdgAAGGMHiKgiPZsxNkREfwXgZwAmANzDGKtsVZcRnGpfE9E9AM7lZmsi+hCAf0ozMWPsawC+pvnuFgC3+IynKkUowvXFctZQPV/UPFszmjYWeWksKs1DVwTChcHZtGYbQ3GByFxUWrELkmgcaaGj7bwFsutG1nadLefYBFehrkvrMcF1c5n0ObgivR85gEiHNmEsfqerYV+E2frY0RNoGk6+mZ2dmY3FX6hgim1ymYOIuhBo1KsBHAbwT0T0dsbY5xIs2QrXt/olAKaEz1MIHOM1iaRE7sPkbIFNWbywaQW/C0zMy6WDjG/bt+72NmXwHRD/3dxMsrYGZnEkFWa2GuoRVi/De71Gdkce0bu+1bx8Nrb2c9x/BxfBnKbMpziP6Xv5HNeOcRwDJ3dkX4e8o7PCVbe5z+29UL134jjR3zXigkgKU2wTEQ0TUV+oJfcBGFGcdjGAvYyxQ+E1XwXwywCqKpT/EcAPiehrCKLPfhXAZ/NYUFr4Mq8szdXVhOk+krxUouaRdXs3mREmzpX20LaA5Jrl+q62RNdlhUXNGt9+CvO1a73sPOG6/qKEgm0eH5eN7b5sPm7TWiqyMXIqiFTvwtgR2wBcA+DW8H9VVtDPAFxARIsRmK9fB+ARxXmZwDX6+hYi+lcArwoPvZMx9mhei0oCXrR/oSFrYcxh6oiUZRRtVi++yXyZdo4kfve8Ueu+5DSoV2GQZN15XmP0mXvM2+Bdom4F8GUiuhaB8H0rABDRSgCfYoxdxhh7mIi+AuDHCOKcHoXBXJ4WrrWvXwLgeQg+YCJ6CWPsZ3ktzBe86lHS3GAdfMbM2/fkOh+QDWNTCeZa7K2qYhoujMR2L1EtYKHZhKrMZt5IStsqZF7DPetyljUskNP0Q88bqhaiOm2cH7e9I6r3v9HAGBtFoPnKx/cDuEz4fBOAm4pYk+ub+C3MF/DoQODwfgZBfeqaQEtYciTrKke+QS2u5+RtHswq6jVvIewq9PNgDqZeyCqBLKIamxNfH2aScXyRxaZVRF7R2ln3vU6Sy+8iBH1h2pDKwtl1HeJ5tbgJb3S4mq9jLUiI6FwAv5vLilKAR1PmXX4wC6aWJBq7kU2WWcHUsN208xcZUOyaUDOWO+zUCooIBix6zGqX7fSBzzptmqevELQJV9V4vpvbUjAXj0RvE2Psx0T08qwXkyXyEMx5CsU8xrb5lIoq32kCZ1S2F18lMFXf246bxjExLC6Qi9TkOVrIXJUqS021VpB1Kp8oEKspaIo2CWcxV6ObsGsNrj5lsbJXE4BzARzKZUUZIknqy0KBisFXqwqTD4NM6udOw1hcfc95w6VcpG9ZzIWEtII4bw2+GhuFamwyS5jh+raKiW8zCHzM/5z9cmoXtczY0lTYqhczoYjSnKaHWNo0Dc1mWbVNNz6QrrBI0chrXTp6tj0bm/XIBXlkVJRIB9e3didjLFbBi4jeipRVvbKEXB/YF76FFExIytC8826F89PWqa5lFMU4kviNi9AoktB2WoEs/g9kR0cyndaqAM4TNjr2fdZZCedSMNcGXN/cG1ApgFXHqgZeHzjrwvuJxnNt2JCQIflcZ2KseUe5AtkI0qzG0DGeSOhLAtnHbN1IDE1lpjVVttKhkQRuraQG+nbhytNtUyIf2LpEvRFBrlY/EX1M+GopAjN2zaClSdWG2Q2qYJlcolRriElVK+3E1jAgT9iCwVxSpBYKbP7TWqLlLOGTLiTDpLGqaE8nYF2qebn6t8V5XaK1S1QfNsmzH0E5scsB7BCOjwP4o7wWlRQ+JfuyKlXnilpgYnmvwUW4qoRdrWghaYJeGpGhufp7a8liYIyg96RPV+EqRnWL5+VRqjZtlT4RpZCuTdi6RD0O4HEi2soYqynN2AfqxgjuxF2Pftlq53r67ORV5TxlBlJr/XlFuKZ15QWXZ+JTcjFNQYxqPIe8imEkieoXr5EFvDxetWm5YuNQo7n4Cw028/WXGWO/BuBRImLy94yxjbmtLAHyInLVuEkKyRcFXxNXHvOnKaqgE8hJxk4DH02h2hqiDUnjEGzwzbstwiriY7JNg1qnD++cfkEgZxE8ViIZbDbb68P/35T3QuoNaQRD3oKl2huDvHo35w0X02c1Te26oL1qW0Vcn0HepVJVaxF/tyx/q7yFlS0oEcg2AFI+luUcJfxgM18fCP/8fcbY+8XviOgvAby/8qrqoCV5nFehqLYWWzSSBHYVmTvpy1yrwaR4RS8R9Ug7qo1NUmFZq7+b7X5M3+u6smUxlgnVFLyTLxwE2/dM4utnp05kuJragGt00yWoFMBvVBwrYUE1mKmPdpeVRmGr8SuvxaVVZNro6LxNcrVm6qvV3FPZvOxLn0UhqVVEt6nUlfnMoviH7zVZ1tgukS1sPuXfA/D7AE4noieErzoBbM9zYSWSwTVS2LVWdJbzAvNdl6JGD8JadP5Jl044Sdcor8d2rc7cV0uMSxX1CyT/bUVTedabSpdAP9dnm6Qcah7CSSeMaynNrhbXVCKATVP+PIB/BfC/AHxAOD7OGHshzcREdBuANwOYArAHwDsZY4fD724AcC2AWQDvYYx9J81c1YQvQ0sTZZwV8/KdY7az10m70bVBlMfyDRxKvJFwWI9qPhVqRSN1sVD4rlWsUJeH20W2YPiuL03tc5s/NQszcJL7SaU9e0RRl0K49mDzKR8BcATA2wCAiHoAtANYQkRLGGM/SzH3vQBuYIzNhP7pGwC8n4g2APgNBL2aVwK4j4jWM8ZmU8zljCyZrhyhnQVDM41RK4JBB1dGkcTMnISRzS5dkSoNxBQlXi24PIck2nMRbpesrTccaTZ5rrSY1bsnWo3keZNqtUlN5bXOTxoVrl2i3gzgowiE5AiAXwCwC4HgTATG2D3Cx4cAXBX+fQWALzLGTgDYS0S7AbwCwINJ5wLcNCqblpFEwxCRtjRhLQSJuQo/X1NhFjv2RObLjARyrpib8aI/n01NkUF1ecBmWUrqcwXcn4krXWdtktdFZ6tSm5LOUaJ4uAZ6/TmACwDcxxg7h4guQqg9Z4TfBvCl8O9+BEKaY194rAJEdB2A6wDgtNNO0worXYCFDFskY9rAmbSCNGnBk7zMjSqk9dvZ5tT5l5PM4Rq5ajOvtnd0AAD6O7yXoIVI2y9ZtTLRGHkHE1WTmbvQd9qNHqcrl1Q5l7mLMBX7CGHbejhdlygWrkJ5mjE2SkRNRNTEGPteaHI2gojuA6BSRW5kjH09POdGBHW0t/LLFOdXFC4BAMbYHQDuAIBN55zLim6tpkLWeYRJUaQmnYeWaxMEtuAgGUkiaW3n5Mm0RNrefO65zCRgVeuUhZZJuBQR7Zz1e+FL30kj77MSyC6w5QynhU+AYimQqwdXoXyYiJYAeADAViIagUNDCsbYxabviegaBIVJXscY44J3H4DThNNWIajBbUQLm9G++NXKd02bUmFCEkaqY+xpzPpAZcCUbBK2aZ6qdcgMw8XXZ3Q/GKK+ZaiOt3d0AB0DAIBW7SzFwIWOXLsI5bEGH5+2bSzXMVVjpApoLLiNpy4DIakP2NUV19ozUHV6LhGHq1C+AsAkgiYUVwNYBuDmNBMT0aUI8pxfwxg7Lny1DcDniYj7sNcB+KHruEkFX5IXOMmu21dDS/Kiu/h90zKsivE8GJhrJKzv9z6Qn7HuedSTxpCllsWfiSlNzcW64XKOiw83by0yWotA2ybBnNay4ztO1qgnul5ocBLKjLFjwse7Mpr74wAWAbiXiADgIcbYuxhjTxPRlwHsRKCNv9s18jqPVB+X89LOY5ozqbk2K4ZlCo5S5fjmwWTSpogA5jXqfMT1Dl8hpqIxn/NdYIsNkM9TXZsUNouNLUXOGnMivQ8my0/WGwrX8Vp7BjKdt0T2sBUPGYfan0sAGGNsadKJGWNrDd/dAuAWrwGb3Fsvyi+HSXNL6otKglqLgrRuIByiPLWalkPRjjyeh2mjVc8Mq97yTbPYaCWdU7Upk6OWVdeqUCHoC6Znl3Hrma4XImx5yp1FLSRL+Ppkk0b2ZuEL85nLV1u2QRfBmrWZUGl6DHOEs5rDNr8OdSHM5vTxEkDyGIQi7t2XpnXIKpvAZS1JaTGJC8d1059EOYj8xs5XlKgFuKuXtQ6JcSWBKfUl6c7e97okglleqyuKjNBWMRWRiZnuO42ZFDBvPuqxh6zNgpOHqdQW52DaKCQRzOI1aQSyz4alGoFwvtHdWSsCJWoPjSOUM4RPsJhr4JaL8Ix8UwmZmG18n/NMUBVBcWWWvlYJXwYjX68r2BI9hzoTyBw2wZjlBscEmxXJlaZV6+PXmGoP6ObzQVKLk6tFyWVNojvHtBbTBqUUxv4gorcC+BCAMwG8gjH2iOa8SwHcDqAZwKcYY7fmtaYFK5Sz0iSSRFSbznFhYhVrEF5ok9BNmmohgzODNHW6TQFs4ve+48lwYVz1buYrOmjIZz7xd07y27qaeSvoSHSNZLDx0s5rsbRkLSjFOuSqefh6ai0+pYbxFIArAXxSdwIRNQP4WwTdEvcB+BERbWOM7cxjQQ0rlF27/wB2QZB0982hGzdpuUxT8JRLHq7vXNH40pzdBsFvM3dWXCcxN77+LO7B10RYD9D91i4RzKpo4yItBkloxMsFJEVRi8/K1xJjEsaqv33pPrrO0ZUDlBpylmCM7QKAMANIh1cA2M0Yey4894sI0oRLoWyEEH3tkmuoYkgu5lPVeK7R3DJsL5erGVvHUH3MZtG4GeQbu84tr0ElmIHkgrMRmZWOtk0uDR2KEMRZ0Yh1noQbDP4eZbWJTepuSfosVAGUdbTRHJk9/nzii+emJwHGYiVpQ9wRVsTLCv0Afi583gfg/AzHj6FxhHJGyIKRZeGfUh1XpXKYrrXO6WBNcNEU0jJ3117G0fkp/M51xLAqoMqjTVvIJSsk/Q0y93E70JLJhy2vT0yj0gUpatfiQKc6i1AS6512vrASXQ3iX6YP7cbc6lehqW2x98Un/vNBLDrtPBx7eltUklYFl3LPFjiXfs4CDSWUXfINK66pUqBPGpO4MVjMQ1imufesn5uPhSEpI69ngQwIFh1LkQvT967unCKeVZYCWb4vWzS4y7rkqOy07i3XzUgjWnhUYIxNL9l4JSYHt2Px+ku8rp2bnsTUwacxO36w3WEeY7lnByQq/ZwUDSWUAfcXx2bSVp3r8lLZfHpZvHBGH6ljWcCs5rOhaEEop0LpXASqdU1OTNRmNa+mlnS/n0d8RXRNhtYJ29hJkSTgzGcNundZDrTq9tC8XeaVAz5Va1BdoxtncmICQG1Wqjv25Nfamjv7puYGtnhpy4GW/HIce3rbiRyXx/EjAOuIaDWAIQC/AeA385qscYSyZ89ZV3gxMklLTZPOk4Rx+eZE+4zL4Tu+/JvIgjKJluECsfc0EPjvdQwzre+61pFn9LF17pQuliwh0l6SgC9+TZJcaZPQ1JnNTetqlHzlJNqyj5ZsAxH9KoC/AbAcwLeI6DHG2BuIaCWC1KfLGGMzRPQHAL4DoBnAPzDGnk47t3ZN882Z6hubz97IHrpnW27aqQvS+lnTCD8XC4FqfJ0GYEt7Sro+eTybRpYUqntxWXPakoREtIMxdl6qQSRw2s4CTtaeMIDMxV9to/kkNO0ryF3zelVIskG0ze0rSEWhnSbN0DZ/GtrOg66FsVubO/umOl/xTidteWL390Cti3Hs6W3GkOl6ReNoyiHyNI25BIG5+vuyzm10ihxP8WxkDdeUgmMKaJGZjazdu5jrKsyHCgZmal2oSgkCqhdfUBRcBWTkrlHQs42+ZYh06WvJcXVFuZp4OXQCz3V9cl0AVQCmaRzdu8TXNjo5l7iKWfP4cCJ3RTXhoy1nqSXXKhpOKGeBNNq2S3CVjrEVbW6S79Mk3HSMQmeWA/TCTz5Pm3etMftxxuVS6Uk3Z+xY+JvVauF+mp3OpBxoIk3QoBmbzkszb5ZuqKxLyapS0rJ2m8lWqtj8mrlMz7hW6VqEq2+5YF9yVdBwQjmPWsbV9M243I+rf8w0lq8ZOclOXt7Bu2hRutQPk0A2mV1NaS3TI4N1wcBMSBpX4GoyzsQ/7WgtSrI5dg2kSgrT++P6DE3vi3itXJxHDv7SbX7rka5dtOWFoCUDDSiUfeAjCFTX6s43MUZf81+W0dTKiHPZFCjs/sV5dOX9VNBFurv0q3XZUJhM5uIc0/v2oHXVGuucqjXXElhzq7MwTOsXzStYEEgXc+GtbRsEv02zTUoTSTYMKj4iFtMB4u+PqriOy5j1IJht2vJC0JKBBhTKvi97klQJm0BMkrLkC9ML6ewbUwhk8X9ZOLvu8HVIW3BBV4VNnoOfZxPIqusbAb5CJY+shdj4CX33Pq4I2yZYV4XP9KxcvtONV7EO3fpUY0suLp8Khap87XqhcZO2vFC0ZKABhbIP0jKjtMSexwtj3BBY5uPMwEWjdh3TF86CRFpnLLVJ7tVsKCxRLwzLFUnpmb8LWQRiiWMCZrpKA5fNrzXw0nET7rJB8AkSi12rilrXrF9rOfN4Z2sZOm15oWjJQIMJ5STMVhbMLi9XUmJ3CYrKE1lps75j5nG9LJDF/5vH9VXdbJGytWzmc9IMLRqbCVkKZH6OKJirgbSWM/G4Kksgi/fX1rbRxeVlyoPmqGXa5lBpywtJSwaqKJSJ6CMIOm3MARgB8A7G2P7wuxsAXAtgFsB7GGPf8R2/3rQhXQ5xra3dR+PxNeUZx/LwrbukZsnn19pzluErZE1+9yyQxMpUzRoCqnWY5vbd1HMaUj332NwOGrypyEnFeFJ+MwBtuhVQH+1JZW15IWnJAJBtvoAfbmOMbWSMnQ3gmwD+DACIaAOCMma/COBSAH8X9rM0I+wS5eL/qRb4S+vjt3YJSMkLJqYin6cSdK7jNh89WGmmU4ypgul5uozBmZl43vTIoHXeQtFUuXf2+e1dac4XPuO6rrcaNK36Lovna3rfZ5euiP7xzxyqXH7dPzlV0QZecrOWwRibbh+4AJOD2yMt+fjObywILRmooqbMGDsqfDwJ8103rgDwRcbYCQB7iWg3gn6WDxoHDMtsAmotwaQ5V2O3Lq9RtSadYBbvU6d9ZKGVJNUW5GuT+iqtbgRNQJrL2BxZ57HmAoG2OUwbjqzo2fl3SBBUqJtL/lt1vbwuH59vWu3YJTDUNcrcxptkLVhXR8C1el2tKCcu4Noym51eUFoyUF1NGUR0CxH9HMDVCDVlqHtX9muuv46IHiGiR54fHXWas9bNlByuPrsSjYkktN2oqJd3tijkFfNSS+Da8vTBpxaUlgzkXPvatY9l6ENuZ4zdRER/C+BBxtjnwu8+DeBuxtg/W+Y6BOAYgORds2sTp6Lx7glo3Ps6gzHWmeWARDQO4Jksx6wRNCINNOI9ATnQtQ2h23IxY2y8yHmrjVzN1x59LD8P4FsAbkLC3pWMseVE9EheRdOrhUa8J6Cx7yuHYZ9p1GfVaPfViPcE5EbXRjDGZgEsKIEMVNF8TUTrhI+XA/hJ+Pc2AL9BRIvC/pXrAPyw6PWVKFGiRIkSRaOaecq3EtEZCFKi/hPAuwCAMfY0EX0ZwE4AMwDeHe6YSpQoUaJEiYZGNaOv32L47hYAtyQY9o7kK6pZNOI9AeV9VXvMWkAj3lcj3hPQuPdVc8g10KtEiRIlSpQo4Y46SNIsUaJEiRIlFgZKoVyiRIkSJUrUCOpWKBPRyUT0FSL6CRHtIqJfIqJTiOheIno2/L+r2uv0BRH9ERE9TURPEdEXiKhFpXfmAAAIr0lEQVS9Hu+LiP6BiEaI6CnhmPY+iOgGItpNRM8Q0Ruqs2ozNPd0W0iDTxDR14joZOE773sq6bq20Yh0DRRD2yXcULdCGcDtAL7NGHspgE0AdgH4AIDvMsbWAfhu+LluQET9AN4D4DzG2FkAmhHUAa/H+/oMgtrlIpT3kbjeefH4DCrv6V4AZzHGNgL4KYAbgFT3VNJ1beMzaDy6Boqh7RIOqEuhTERLAbwawKcBgDE2xRg7jKBu9l3haXcB+JXqrDAVWgB0EFELgMUICqfU3X0xxh4A8IJ0WHcfUb1zxtheALzeeU1BdU+MsXsYYzPhx4cQFLsBEtxTSde1f1+NSNdA/rRdwh11KZQBnA7gEIA7iehRIvoUEZ0EoJcxdgAAwv97qrlIXzDGhgD8FYCfATgA4Ahj7B7U+X0J0N2Hc73zGsdvA/jX8O8k91TSdX2i0ekaSE/bJRxRr0K5BcC5AD7BGDsHQc3rejB9GRH6oq4AsBrASgAnEdHbq7uqQkCKY3WVq0dENyIodrOVH1KcZrunkq4bC3VP10BmtF3CEfUqlPcB2McYezj8/BUEzGyYiPoAIPx/pErrS4qLAexljB1ijE0D+CqAX0b93xeH7j4S1TuvFRDRNQDeBOBqNp/4n+SeSrquTzQkXQOZ0nYJR9SlUGaMHQTwcwrKdALA6xCU5dwG4Jrw2DUAvq64vJbxMwAXENFiIiIE97UL9X9fHLr7qNt650R0KYD3A7icMXZc+Mr7nkq6rrv74mg4ugaype0SHmCM1eU/AGcDeATAEwD+BUAXgG4E0Y/Phv+fUu11JrivDyNozvEUgH8EsKge7wvAFxD4D6cR7KyvNd0HgBsB7EHQovCN1V6/xz3tRuBfeyz89/dp7qmk6+qvNQEN1DVdG+4rU9ou/7n9K8tslihRokSJEjWCujRflyhRokSJEo2IUiiXKFGiRIkSNYJSKJcoUaJEiRI1glIolyhRokSJEjWCUiiXKFGiRIkSNYJSKJcoUaJEiRI1glIoa0BEL+Yw5uVExDvI/ErYbcV3jPuJ6DzP858hossV3w2IrdoaHUT0DiJaKXzeSkQvENFV1VxXkSjpuvFQ0nVjoRTKBYIxto0xdmv48VcAeDOvhLiaMbYtzwnqpHXbOxDUXgYAMMauRlCdqEQKlHRddbwDJV03DEqhbAEFuI2C5uxPEtGvh8cvDHfrvCH91rCEIIjosvDYD4joY0T0zfD4O4jo40T0ywAuB3AbET1GRGtETYGITiWiwfDvDiL6IgWNxr8EoENY2+uJ6EEi+jER/RMRLXG4n81E9DgRPQjg3cLx5vA+fxTO9bvh8SYi+jsKGtR/k4ju5jtwIhokoj8joh8AeGt4H9/+/+2dXYhVVRTHf//JaswpIzULKTRLhigTMkxmcHwIo3qI0JLwpQ+SipKKNIIKwoeKoqCSMiIGwqbQqSAtx8hGzbRvZ8wygvKhByNo1OxDyFYPa93mNN17BzGZM7p+cDjr7r3v3mufu2af/TV7SfpM0iZJzZFunKTOyPsTSS0R3hb13yb3inRyHb0XF3R7uBD+ZpS3Q9LCQl3aC7/Z3aHzdGBFlDeyVlnHAmnXaddJSRnqI8XKegH74z4Xd/Z9HDAeP8f3TGA2sBc/jL0B2AK0Ao340XST4vsdwOqQbwCeDbkdmFcorxt3Ag8wFtgV8j3ASyFPxb21TI80G4FREXcf8FCVevyTb3zuBdpCfhz4MuSFwAMhn4gf9TgJmAe8HXU8A+ir6A3sApYU8n4POC/kGcD6kF8BWkM+G/g65LeAlpCbgBE1fos5wAu4d5oGYDUwK+JOi/tI/AjHMcDFwLuF759a7VlU+x2O9ivtOu06r3JfI0gGoxXoMLODuDeYDcAlwD7gYzP7AUDSNmAisB/4ztz5N3jjtfAwyp8FPA1gZr2SeiP8UnyacHMMZE7AG9CaSBqN/yFviKCXgStCngNMVf861Gj8oPlWYKWZ/QXslvT+gGxfi7ybcM8/K0Mf8EYQ3EvQ+YXwU2L0sBl4UtIK4PXKs6zCnLi+iM9NodtGYJGkayL8rAj/BjhH0jPAGmBdvedyjJJ2nXadlJB8KQ9ONd+hFQ4U5IP486yXvh5/0r+c0DggrtoB5cJ7zdcfQhmqkVcl7k4z6/pXoHTVIHn+GvcGYI+ZTauSpgGYaWa/Dwh/VNIa4Epgq6TLzGxnDd0eMbPlA3SbjTeMM83sN0ndQKOZ9Um6CLgcn8q8DnfSnvSTdl2ftOtkSMg15cHZCMyP9ZxxeA+/npuynXhvdmJ8nl8j3S9Aca1pFz49BT61Vix/AYCkC/CpPoCtQIukcyPuJElT6lXEzPYAeyW1RtCCQnQXcJuk4yO/KZJGAR8Ac2MNbjw+vVkt733A95Kuje8rGhDwHv0dlbSSpsV9spltN7PH8GnF5hqqdwE3VdYWJU2QdDo+6umLhqsZH2UhaSzQYGadwIO4T2L47zM/lkm7TrtOSki+lAfnDXy9qgdYj6817a6VOHrNtwNrY6PIj/ga3UBeBRbHRpDJwBN44/Ehvq5W4TmgKab3lhANp5n9hK/ldUTcVmr/8Re5EVgm3xBT7OG/iPvu/Vz+7yTL8RFSJ+7KrRL2UY36gDeGN0vqAXYAV0f4ImB6bGb5Crg1wu+KTSs9ocs71TI1s3X4+t0WSduBVXgjtBYYEfVfGs8AYALQHVOv7cD9Ed4OPJ8bYoC067TrpJSk68YjgKQmM9svX2xaBnxrZk8NkS7dwL1m9ulh5FGpzxi88Wyp14APJyS14xuWVg21LmUn7Xr4kHY9fMmR8pHhlujN7sCnopYPkv5I8jPQriqHLBwCq6M+m4ClR1HDtQJoA/4Yal2GCWnXw4C06+FNjpST0iDpQnznbJEDZjZjKPRJkv+DtOvkUMiXcpIkSZKUhJy+TpIkSZKSkC/lJEmSJCkJ+VJOkiRJkpKQL+UkSZIkKQl/A+KlcTIyoUxZAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ds_mean.sel(latitude=slice(30,-30),longitude=slice(60,120),month=[4,6,9,10]).u.plot(col='month',col_wrap=2,levels=np.arange(-1,1,0.1))" ] }, { "cell_type": "code", "execution_count": 56, "metadata": {}, "outputs": [], "source": [ "ds_en4 = xr.open_dataset('/home/lijo/WORK/oceanCurrValidation/data/en42003-2010.nc',decode_times=False)" ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Dimensions: (DEPTH: 42, LAT: 173, LON: 360, TMON: 96, bnds: 2)\n", "Coordinates:\n", " * LON (LON) float32 1.0 2.0 3.0 4.0 5.0 ... 357.0 358.0 359.0 360.0\n", " * LAT (LAT) float32 -83.0 -82.0 -81.0 -80.0 ... 86.0 87.0 88.0 89.0\n", " * DEPTH (DEPTH) float32 5.0215898 15.07854 ... 5050.9897 5350.272\n", " * TMON (TMON) float64 2.436e+03 2.437e+03 ... 2.53e+03 2.531e+03\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " DEPTH_bnds (DEPTH, bnds) float32 ...\n", " TEMP (TMON, DEPTH, LAT, LON) float64 ...\n", " SAL (TMON, DEPTH, LAT, LON) float64 ...\n", "Attributes:\n", " history: FERRET V7.4 (optimized) 5-Jul-18\n", " Conventions: CF-1.6" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_en4" ] }, { "cell_type": "code", "execution_count": 58, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "4.017810764" ] }, "execution_count": 58, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_en4.nbytes/1e9" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "ds_temp = ds_en4.sel(LAT=slice(10,12),LON=slice(80,82)).isel(TMON=5).TEMP\n", "ds_sal = ds_en4.sel(LAT=slice(10,12),LON=slice(80,82)).isel(TMON=5).SAL\n", "t = xr.where(ds_temp<1000,ds_temp,np.nan) - 273.15\n", "s = xr.where((ds_sal<1000) & (ds_sal>10),ds_sal,np.nan)" ] }, { "cell_type": "code", "execution_count": 109, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAU0AAAGgCAYAAAAuKjWYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de5yUZf3/8deH81EBRUVAQQVUUFFJS81MJQ95yjKx8lCW1VdTs/LYNy3T1F+pfStPaYamIJ7xLCimZamAohxEOayygicE5Ay7+/n9cd0jwzK7O7s7M/d9z7yfj8c8Zue+577va2bv/ex1vszdERGR/LSJOwEiImmioCki0gwKmiIizaCgKSLSDAqaIiLNoKApItIMCpoiIs2goJkAZuaNPP5iZs9FP99X77gHo+3PZW3b3symm1lttG+Nmd1W77iaaH/vrG1/N7OlBf5cp5vZYjOrix4LzGyPHO/ramZrzaymgfOcHX2WF7K27WFm72d9zgPyTFOuc5mZPR19L25m75hZ36z9r5jZumjfWjO7ud456//O3szad4uZrYy215jZS2bWMWv/0nrHrsvad0COc0/M2v9TM1uSOXeOz7ok+t7dzFab2RX19t9rZuuj/SvN7MdZ+9bkuPb7WftnZX0nt+bz3ZcLBc0EcHfLPIBa4PysbWdGb1sPHFHv0MOi7UAIPsBMYCvgIKALcDlwmpk9XO9YA/5R8A+zsW2BvwP9gd7AamBijvc9GO3bhJl1Bq4GVtTbVQM8B/wo38Q0cq6bgS8B+wHbAB2Af2btXwkcC7QDfgz8wMzOqHeOQ7J+Zztnbe8GnB897wsMAx6qd+zorGM75Eh6p6z9h2ZtXwLcQ8O/x+8DnaP76kzg4sw/LTP7HvAN4DRCHHgY+IuZtQdw905Z92QbwvedfQ+9DlwMrGrg2uXL3fVI0INwc/6i3rbngHnAWuCYaNsx0et5wHPRttuBOqB3veOvBxzom3WNJ6L3bh9t+zuwtMif7Vvhltto2xejz/G/QE2OY54AXgLmAC/k2N8x+mwH5HH9nOcCqoFHs16fEZ1ziwbO8z4wPuu1Awfn+R08DLyf9Xop8PcG3ntAdO6OTZzzF7m+u3rvOS0612lZ98SKrP29o/175Dj27Ghf7xz7lgG3FvO+SdpDOc10eY7w3x3gEmBSvf0HA/Pc/aN62y+Ink/N2vYssIg8c5tNVCE8kWf6T2DTXN59wHXA8hzX3A84hPAPolWaOJdFj4zM38XBOc7TkxBgXqq3a0JUVfBeE1UF+wDz6207JSpGf2pm5+Y4ZlVUtH/bzAY3cu5NmNkHZuaEf6gfA3dGu/5f2G3fi3KXNxBy+6/nOM1ZwNwc91VFUtBMl0uBEVExc2/gV/X2dwM2ubHdfS0hV9m33q4zgQPMbOf6x+Q4hzXyqF9tsAkz+zqhiHth1rargDbufmEDh90P3ODuHzR1/jw0dq4JwMioDrE/cFm0vUeO974CfAJcmbXtJ8BmwA7AYmBidr1lRlS33JuN/3mdR/i9bA6MBa4zsy9H+94GTiFUs+wBdAZebPKTZnH3raPjfwX8291ro10LgX8BtwHrgK8DP/Yo+5iV5i2AQcCfm3PdshZ3VlePjR80UjyPfv4EeBn4JHqdXTx/B5iT45yZIuzF9a9ByPVMpojFc0KOrRa4KWtbb8If66HR63PJKmICv8l8xuh1i4vnTZ0LaBt9xzXR46HonPvUO8/LhPrNvo1cq3107PH1tv8u+g6Ob+jY6H0fAfc2sG83sqpZsrY3WTzPOvdvo5/viL7/r0Sf/5IofXvWO+aG6DuxBs6p4rkk3gPA5wjF2vomATtkt4pHro6e72RTpxNyrds1dtEmiudPNXLcfsDTwFh3z260OZgQYJ4ys1rgWqBtVMQ9ADga6Bm9rgV2JOSKFzWWzgY0ei53r3X3g9y9nbu3A/5LCCBTsj7Hc8DuhDq/9/K45mfFfTO7hFBFcrq7P9CcY+upi55b+nfbhtAQBTAcmO7uT0ef/wpCEP12vWNOJORQNR1aRtxRW4+NHzSd0+xJyFn0jF5n5zS7E1ozPyQ0InQm1IHWsnFDx0bXAN4k/EEWNKdJCMbrs6+dta8jIeeUeVwdpXM3QjDtU2//u8BUYIesc2zOhgaMw4HNG0hHo+ciFKu/TAhWRwNrgDuzjn+SEFByNZIcQwgs7YGtgVej93aO9p8Xfbdn5Th2++j3s3n0ffwl81mi/d+LPldbYKco3Uuyjm8bHfu/0e90c6BrtO+IaHvP6D64ITr3t6L9t0Tp/FL0uS+I9h+Rdf4R0bYv50h71+h6nxJyrZsDbeP++ynJ32jcCdCj3i+kiaCZ4/2fBc3o9Q6Ebke10Q2/hnqts/WvQWiccAofNJ+NzrvRo4H3blQ8z7F/k+J5Y+cGZgAz8jkXoYi6LjpHDfBwU9cBnor2nZd1bB2hce3QrGOX5Dj2o2jfzoTifmb7cuCCrGP/RPin49Hvcy6wW73vrP65l0b7vkpodMtsXwlclXWsEbpV1UT715JVfRK950lgWQPf4dIc1z437r+fUjws+gJERCQPqtMUEWkGBU0RkWZQ0BQRaQYFTRGRZlDQFBFpBgVNKRtmdomZzTCz183sNTPbN9reO5oC7Yf13l9lZlvGk1pJq3ZxJ0CkEMzsC8BRwF7uvjYKhplp1k4gjPI5iTANnEiLKacp5aIP8LGHyUlw94/dfWG07yTgZ0C/7MmFRVpCQVPKxdNAfzN7y8xuMLMvAUSzFm3j7i8D4whDHkVaTEFTyoK7ryCMdT+DMJvPPWZ2GjCKECwhTL12UiwJlLKhYZRSlszsG4R5K/sSJtLILAuyLTDU3d82sypghLt/HE8qJY2U05SyYGZDzGxQ1qbhhIbOru7e190HuPsAwryWo+JIo5QHBU0pF92A0WY208xeB3YlzAr0YL333c/GRfTXzaw6elxborRKiql4LiLSDMppiog0g4KmiEgzKGiKiDRDqoKmGZ3MeNmMaWbMMOPX0fbLzXjdjNfMeNqMbRs5R1szXjXj0VJf34weZtxnxptmzDLjCyW89k+j46abMcaMToX47Fn7f26Gm5FzLLcZh5sx24w5ZjS0ZG/Br21GfzMmRd/3DDPOac61W3v9rPe06L4rwPfe4nuuQNdv8X3XyD1/mRnvRff8a2Yc2cDxLb7nGhX3ehvNeYAbeLfo5/bgL4F/HnyzrPecDX5TI+c4D/xu8E0W+yr29cFHg38/+rkDeI9SXBu8L/h88GixLx8HflohPnv0uj/4U+DvgG+Z49i24HPBd4g+9zTwXUt07T7ge0U/dwd/qznXbu31W3vftfbarbnnCvDdt+q+a+Sevwz8500c26p7rrFHqnKaUZpXRC/bRw9359Ost3UlLPK0CTP6ERacurXU1zdjM+BA4LboXOvcWVqKa0faAZ3NaAd0ARY28L5mXT96fR1wfiPX3geY4848d9YRRuYcW4pru7PInanRz8uBWYQO73lr5Wdv1X3Xmmu39p5r7fUjLb7vmrh2U1p1zzUmVUETPivmvEZYpnaCOy9F268wYwFh3eZfNXD49YRfcl0D+4t5/R0Iw/tuj4ppt5rRtRTXduc94PeEJWAXAcvcebo5127o+mYcA7znzrRGDu0LLMh6XU0zA1crrp19jgHAnhC+txJev1X3XSuu3ep7rjXXL8R919A9D5wVVUv9zYyeOQ5t9T3XkNQFTXdq3RkO9AP2MWNYtP0Sd/oDdwFn1T/OjKOAD92ZEsf1Cf9x9wJudGdPwpKqzapnacVn70n4LzuQMIywqxnfac61G7j+7sAlNPxP6rMk5Dpdia4dEmB0I3RsP7de7ryo1y/EfdeKz97qe6411y/EfdfAPX8jsCNh1Nci4A+5Lp/rdM25dkNSFzQzomLGc8Dh9XbdDXw9xyH7A8eYUUXIqh9sxj9KeP1qoDrrP+V9hBu6FNc+FJjvzkfurAceAPZrybXrXT/zBzEt+l77AVPN2KbeIdVA/6zX/Whm9UArro0Z7QkB8y53HmjJdVtx/YLddy383gtyz7Xw+gW777LveXc+iIJpHfBXQlG8voLdc/WlKmia0duMHtHPnQm/lDfNyB5zfAzwZv1j3bnInX7uDCCMPX7WvXn/9Vp5/feBBWYMiTYdAswsxbUJxaPPm9HFDIuuPSvfazdy/Vfd2cqdAdH3Wg3sFX3WbK8Ag8wYaEYHwvc/vhTXjj7vbcAsd1o0TLI112/tfdfKa7fqnmvt9WnlfdfIPd8n621fA6bnOLxV91xj0jZzex9gtBltCQF/nDuPmnF/dGPUAe8APwKw0P3mVvfcXRJiuP5PgLuiX+I84LuluHZUB3UfMBWoAV4FbinEZ2/ozfWuX2PGWcBTQFvgb+7MKMW1CTm9k4E3oroxgIvdebxE12+t1l67Nfdcq65fgPuuoXv+TjOGE4rbVcAPc1y7tfdcgzT2XESkGVJVPBcRiZuCpohIMyhoiog0g4KmiEgzJC5omtnhZjbbzOaYWeEG2YuIFECigqaZtQX+AhxBWK7gJDPbtYXnOqOQaUvLteO+vj57fPTZSyNRQZPPBtn7PHdv7SD7OH+Jsd5AMV9fn70yr18xnz1pQbNog+xFRAohaSOCmhxkH2XDM/9V9u7SpUvOE3Xo0IGuXbvG0nM/zmvHfX19dn32crj+qlWr3N1zZiqTFjSbHGTv7rcQDcXq2rWrr1y5Mu+T19VBm6TlrUWk8NzDo4V/8Ga2uqF9SQsh0SB7G2hmBR1kD/C978Exx8D48VBTU6izikiifPghjBgBd95ZlNMnKmi6ew18Nsh+FjDO3QsyyB5gxx1h8mQ49ljYbju4+GKYO7dQZxeRROjdG9auheuuC7nNAkv1hB3NLZ5DyGE+/jjceis89lgosn/5y/D978Pxx0OnZi03JiKJdNtt4Y/6mWfg4IObfbiZrXL3nLPcV1zQzPbeezB6dAig8+dDz55w8snhu95ttwImVERKa82aUJzcd1945JFmH95Y0ExU8bzU+vYNRfQ5c2DiRDjsMLjpJth99/Bd//WvsHx53KkUkWbr1An+53/g0Udh9uyCnrqig2ZGmzZwyCEwZgwsXAjXXw8rV8IZZ0CfPnD66fCf/xSlekREiuXHP4aOHeGPfyzoaSu6eN4Yd3j55VB0HzMmBNGhQ0PR/dRTQ1FeRBLu9NPDH/CCBbDFFnkfpuJ5C5htKKIvWhSeu3WDn/4Utt8eLroo9GwQkQQ7+WRYvRomTSrYKRU089C9e8hh/ve/8OqrcOSRcPXVMGAAnHtuaFASkQR6/fXwvFeLF+HchIJmMw0fDmPHwqxZcOKJ8Oc/ww47wA9/CPPmxZ06EdnIk0/C4MHhj7RAFDRbaMgQuP320PJ++unw97+H380pp4SAKiIxW70annsODj+8oKdV0GylAQPghhtCP8+zz4b77w8NRt/8Jrz2WpOHi0ixPP98CJwKmsm07bZw7bVQVRUaiZ56CvbcE44+OtSFikiJPflk6HL0pS8V9LTqclQkS5eG+s7rr4fFi+GLX4QDDoBhw8JjyJDw+xSRIli+PDRADBoUgmczaRhljFasgJtvDnWeb765YXaltm1DHWgmiGYeO+4Y9olIC82ZE2blmT0b7rsPjjuu2adQ0EyIdevgrbdg+vSNH/PmbRht1KkT7LLLpsG0f//Qd1REGvHUUzBqVBjmN25cGOrXAgqaCbdyZWhxrx9Ms/t/du++aSAdNgy22iq+dIskhjv8v/8XGhSGDYOHHoKBA1t8OgXNlFqyBGbM2BBEZ8yAN94IdaQZW265aSAdOhR69Igv3SIltWpV6Pc3dmzotvK3v0HXnPEubwqaZcQ9DN+snyudPj3Un2b067dpMN1lF2hgSSWRdKqqgq99DaZNg9/9Ds4/vyD1WAqaFcAd3n1300A6a1aYxBrCvbTjjrkbnzT5sqTOpEkhZ7l+fZiU44gjCnZqBc0KVlMTlvTIDqQzZoQGqdraDe/bZpvQUX/77Td93n77Vpd2RArnlVfgmmvCSJKdd4aHHw5diwpIQVM2sXZt6JExfXoIqu+8E0o6VVUhx7p+/cbv7907d0DNPG+2Wck/glQS99Df8pprwtDIHj3CfJkXXliUm09BU5qlri5Mh5cJpLme16zZ+JiePRsOqAMGhHtcXaak2davDw0811wT/sP36wfnnRemHevevWiXVdCUgso0RjUUUKuqQjeqbN27h+DZUGDdcksFVcmyfHmYAfy668IEwsOGhUaeUaOgffuiX15BU0rKHT75JHcwzTw+/XTjY7p02RBEcwXWrbdWUK0I778Pf/pTmAVn6VI46KAQLA8/vKQ3gIKmJM7SpY0X/7P7okIYp59plMoVWPv00fDTVHvrLfjDH8LysOvWwde/Dr/4BeyzTyzJUdCU1Fm+PATPhgJr/aVG2rcPQ02zg+nQoWHOhoEDw6g6SaDJk0P/ygcfhA4d4LvfhZ/9DHbaKdZkKWhK2Vm1KrTyZxf7swPrwoUb3tu9O+yxRwigmcfQoeqbGqvXXoNLL4Xx40Mr4plnwllnhXqYBFDQlIqzejXMnBn+NjOPadM2rGPftm0YIVU/mG65ZbzpLnszZoRgef/9oUvFz34WZu9OWJ81BU0RQleq+fM3DqSvvQbV1Rve07fvxkF0+PCwvIyK9600ezZcdhncc8+GZV1/+tPETpKgoCnSiI8/DrnQ7EA6a9aGEVPduoUcaXaudNgw6Nw53nSnwpw5cPnl8I9/hPqQs8+Gn/+8WWuQx0FBU6SZ1qzZtHj/2msbivdt2oQRfPVzpb17x5vuxKiqgt/+Nsy+3b59qLM8//zUzGWooClSAHV1IRbUD6QLFmx4z7bbbgigmZzpTjtVUPF+wQK48kq47bbQr/JHPwpDHfv0iTtlzaKgKVJEixdvWryfOXND8b5rV/jc50KVXoHX+EqONWvg178Oqwu6h2GOF18chj2mkIKmSInVL94/9FDIhJ1wQphgfPvt405hAU2ZAqecEj7wKaeE4DlgQNypahUFTZGYrVoFv/89XHVVyIj94hdwwQUpn3Jv3bpQb3nllaF/5a23FnROyzgpaIokxIIFIViOGRNKrtdcE+agSN24+mnT4LTTQjb6lFPCWtU9e8adqoJpLGhWSvW0SCL07w933w0vvBBa2r/1LfjiF0MJNxVqakLu8nOfC/MHPvRQGC9eRgGzKQqaIjE44IAwAfmtt8Lbb4cY9P3vwwcfxJ2yRsycCV/4Avzv/4YJNaZPD+uLVxgFTZGYtG0bFlF8660wr+7o0TB4cJjsZ926uFOXpbY2tF7ttVfoczVuXKhfqNAxp6rTFEmI2bPDUOzHHgtF9yOPhKOPhpEjYxyaPXVqWFbi5ZfDqo833piYSTWKSQ1BIiny9NMh1/nEE7BkSRhQc9BBcNRRIYgOHFiCRCxbBr/8ZZgMuHfvMIN6KlusWkZBUySFamrgxRfh0UfhkUfgzTfD9l13DcHzqKNCFWNBJ192D0Xv884Lk5b+z/+Ehp+ETqxRLAqaImVgzpwQQB99FP75zxBUe/UKxfijjgorQmy+eSsu8OabYYz4s8/CiBFw002w994FS3+aKGiKlJlly0Ix/pFH4PHHw1DOdu1C96VMLjTvpcBXrYIrrgiNPV26hJnUzzijotcPUdAUKWO1tfDf/24oxs+YEbYPGbKhHnS//RpYxPHRR+EnPwmt4iefHAJnBTT0NEVBU6SCzJ+/oRg/aVJYOrxHj1B8P/roUJzvsewdOOccePjhMIX9jTeW8WwizaegKVKhli+HCRNCDvSxx+Cjj6BD2xqO8Yc5ud1YDv/VPnT4xTlhUTP5jIKmSKVbtoy6P1zHS79/gbGrj2VMp9P4aM1mbLFF6El0yilhVFKF9Chqksaei1Sq1atDPeUOO9Dm8l/zha/24o8zR/Lep5vxyCNwyCFhKOe++4aZ6H/721C9KQ1TTlOkHK1fH2ZPv/zysJ7xYYeFFvIcXYiWLoX77oM774Tnnw/bDjwwtAudcEIruzGllIrnIpWitjZ0Tr/0Upg3D/bfPwTLPBt5qqrgrrvgjjvCmPiOHcOcHCefHOJuzhb4MqSgKVLu3GH8+DD0cfr0sEDRlVeGSYFbUFHpHmZhuvPOEIMXLw6jKU86KQTQvfcu7/pPBU2Rcvbss3DRRWFSjcGDQ5H8G98o2Gpu69bBk0+GADp+fHi9884heJ5ySmqXAWqUgqZIOVq5Es49N7Tk9O8fiuSnnhqGBhXJkiVw770hgP7rX2EA0f/9H3zve+WV81TQFCk306aFvkKzZ4f1My69FDp1KmkS5syBH/4wZHRPPBFuvrl8Go3U5UikXLjDn/8c+ggtWxZ6rv/udyUPmBDWc3/66dDOdN99sOeeoYag3CloiqTFxx+Hpuyf/CR0sJw2LTzHqG3bsLz588+Hhvv99w+LxdXVxZqsolLQFEmDSZNCi/hTT4WVHx99NDRnJ8R++4WFKY89NtQWHHFEwtc7aoWiBU0z+5uZfWhm07O29TKzCWb2dvTcM2vfRWY2x8xmm9lhxUqXSKrU1ISFzA45BLp1C9MZnXNOIltdevYMjUQ33RRynnvsEWoPyk0xc5p/Bw6vt+1C4Bl3HwQ8E73GzHYFRgFDo2NuMLPKncxPBEJP8wMPDGMbv/vdsM7vnnvGnapGmYXGoVdegS22gK98BS68MAxQKhdFC5ru/jzwSb3NxwKjo59HA8dlbR/r7mvdfT4wB9inWGkTSbyxY2H48DA55t13hyGR3brFnaq8DRsWAucZZ8DVV4fnclG8Dl25be3uiwDcfZGZbRVt7wv8N+t91dE2kcry6adw1lmhI+TnPx/GNO6wQ9ypapEuXUI3pKVL4Zln4k5N4SSlIShXBU3ODqRmdoaZTTazyTU1NUVOlkgJ/fvfoSLwrrvgssvghRdSGzCz7b47LFgAK1bEnZLCKHXQ/MDM+gBEzx9G26uB/lnv6wcszHUCd7/F3Ue4+4h2RRz5IFIyNTWhc/qBB4ZKwRdeCK/L5P7eZZfwPHt2vOkolFIHzfHAqdHPpwIPZ20fZWYdzWwgMAiogG6yUvHmzg2rof3mN/Cd74R+O/vtF3eqCmrnncPzrFnxpqNQivavzMzGAAcBW5pZNXApcBUwzsxOB94FTgBw9xlmNg6YCdQAZ7p7bbHSJhI79zD/2llnhR7iY8eGsYhlaKedwkfMrNuedhp7LlJqS5aEfjn33hvmubzjDthuu7hTVVQ77wxDh8L998edkvxo7LlIUrzwQmgZefDBMGb8mWfKPmBCCJrlUjxX0BQpldtvDyN7OneG//wn9PpuWxljOL7whRA0H3ss7pS0nornIsVWVxdmtbj6ahg5EsaNCwuRV5C1a8PETO+9B6+/Dn36xJ2ixql4LhKXVavC6mRXXx3qMR97rOICJoS1hsaMCfMmn3JKumdBUtAUKZaFC0PfywcfhGuvhRtvrJyVyXLYZRf44x9h4kT4wx/iTk3LqXguUgzTpsFRR4WW8jFj4Oij405RIriH5YvGjw/VuiNGxJ2i3LTchUgpPfJIWLaxZ8/w8/DhcacoUT75JIwW7dQJpk6F7t3jTtGmVKcpUgrucN11YSbeXXYJaz8oYG6iV68wvH7evDAJfdooaIoUQm0tnHkmnHcefO1r8M9/Jr+JOEYHHgiXXAKjR8M998SdmuZR8VyktWprwxq2d9wB558fOq0XaM3xclZTE+ZU7tYt1G8miYrnIsVSWxtmVb/jjjDpxtVXK2DmqV270D42eXK6po3Tb1ekpTIB88474fLLw1o+0iwHHRRynC++GHdK8qegKdIStbVw2mkhYP72t/DLX8adolTab7+Q43zuubhTkr/ymOVUpJRqa+HUU0MT8BVXhCGS0iLdusHnPpeuoKmcpkhz1NSEcYB33QVXXqmAWQAHHRQWYUtLvaaCpki+MgHz7rtDC/lFF8WdorKQtnpNBU2RfNTUwMknhyGRV10VpnWTgkhbvabqNEWaUlcXGn3Gjg1dis4/P+4UlZVu3cIY9Oefjzsl+VFOU6Qx7nD22RsafRQwi2LPPWHGjPB1J52CpkhjLr0U/vIX+PnPVYdZREOGwNKl8NFHcaekaQqaIg257rrQaf300+Gaa8Ka5FIUmWV+07A2uoKmSC633x4m3/jGN+DmmxUwi2zIkPCchmV+FTRF6nvgAfj+9+ErX4F//KNiFj+L03bbhfk1ldMUSZuJE8MEwvvuG4Jnx45xp6gitGkDgwYpaIqky3//C8cdF8qKjz0GXXPODCZFsvPOKp6LpMf06XDkkbDNNvDUU2GpCimpIUNg/nxYty7ulDROQVNk/vxQf9mpE0yYoBnXYzJkSJgLZc6cuFPSOAVNqWxLloQc5urVIWAOHBh3iirW5z8fOimMGRN3ShqnoCmVa+1aOP54mDsXHnoIhg6NO0UVbaedwvJKf/oTLFsWd2oapqAplckdfvCDMEvE7bfDl74Ud4qEMNPesmVwww1xp6RhCppSmX796zDr+m9+A9/+dtypkcjee8Phh4fBWKtWxZ2a3BQ0pfKMHh2C5mmnaZmKBLr44jAG/dZb405JblrCVyrLs8/CYYeF4vjjj0OHDnGnSHI48MDQqWHu3Hh+RVrCVwRg5szQ8DN4MNx3nwJmgl1yCVRXh5WRk0Y5TakM778f+rSsXRtG/my/fdwpkka4hwXXli4No4TalXi6dOU0pbKtXAlHHx0qyh55RAEzBcxC3ebcuXDvvXGnZmPKaUp5q62Fr38dxo8PfTGPOSbuFEme6upgt93CZB7TpoXnUlFOUyrXJZfAww/DH/+ogJkybdqEyfKnTw8FhKRQTlPK15NPwhFHhE7st9wSd2qkBWpqwpj0LbaAl14q3VzQymlK5Vm4MKxRPmxYyGVKKrVrF1ZLfuUVePrpuFMTKKcp5ae2FkaODK3kkyfDrrvGnSJphXXrYMcdQ/vdCy+UJrepnKZUliuvhEmT4M9/VsAsAx06wAUXwL//nYy10ZXTlPLy/PPw5S/DqFFhfR8tiFYWVq8Os/bttluYwa/YlNOUyvDxx/Ctb8EOO8BNNylglpHOncPS8xMnhlqXOCloSnlwD6v14fIAABpRSURBVBNwfPQR3HMPdO8ed4qkwH70I+jVC664It50KGhKebjuurAY2u9/D3vtFXdqpAi6dYOf/hQefRRefTW+dKhOU9LvlVdg//3hq18Ny+6qWF62li4NregjR4Y5V4pFdZpSvpYtgxNPDIuh3XabAmaZ69EDzj4b7r8fZsyIJw0KmpJe7nDGGfDuu2E1rl694k6RlMA554Ql6X/3u3iur6Ap6XXbbTBuHFx+Oey3X9ypkRLZckv48Y/D/8k4lvtVnaak09KlYZjIbruF2dhLOQWOxG7RIujXL0zo8dvfFv78qtOU8nP11fDJJ6HVXAGz4vTpE8Yw3HNPqKUpJd1tkj4LFsD118N3vgN77hl3aiQmJ54Yiuel7n6koCnp86tfhRlqi1Euk9Q4/vgwC9I995T2ugqaki5vvBGW4P3JT7RsRYXbYgs49NDQFljKIrqCpqTLBRfA5puHBWSk4p14IlRVwcsvl+6aCpqSHs8+C088EZawUJ9MAY47LkwdV8oieoNdjszs+DyOX+PujzdwfH/gDmAboA64xd3/aGa9gHuAAUAV8E13XxIdcxFwOlALnO3uTzV2cXU5qiB1dbDPPmFCjtmzoVOnuFMkCXHMMTB1ahjjUKiOFI11OWpsNeG/Ag8DjY1LOxDIGTSBGuBn7j7VzLoDU8xsAnAa8Iy7X2VmFwIXAheY2a7AKGAosC0w0cwGu3ttI9eXSnHPPTBlCtxxhwKmbOTEE8PCay++CAccUPzrNRY0n3D37zV2sJn9o6F97r4IWBT9vNzMZgF9gWOBg6K3jQaeAy6Ito9197XAfDObA+wD/CevTyLla+3aUCTfYw/49rfjTo0kzDHHhP+j99xTmqDZYGbW3b/T1MH5vAfAzAYAewIvAVtHATUTWLeK3tYXWJB1WHW0rf65zjCzyWY2uaamJp/LS9rdeCPMnw/XXKOO7LKJ7t3hyCPDrEe1JSiXNnkHmtkJUfEaM/ulmT1gZnlPWGhm3YD7gXPd/dPG3ppj2yYVru5+i7uPcPcR7do1llGWsrB0aeiPOXIkfOUrcadGEmrUKHj//dKsIZTPv+3/jYrXBwCHEYrUN+ZzcjNrTwiYd7n7A9HmD8ysT7S/D/BhtL0a6J91eD9gYT7XkTKWGS559dVxp0QS7KtfDTMflaIVPZ+gmcnwfhW40d0fBjo0dZCZGXAbMMvdr83aNR44Nfr5VEJjU2b7KDPraGYDgUFACXtfSeJkhkt++9saLimN6tIFjj46zLNZ7Fq7fILme2Z2M/BN4HEz65jncfsDJwMHm9lr0eNI4CpgpJm9DYyMXuPuM4BxwEzgSeBMtZxXuN/8RsMlJW8nnhjW1ps0qbjXaXJqODPrAhwOvOHub0dF6t3c/eniJq1p6qdZxlauhK22CrnMW26JOzWSAsuXw2abhWXvL7qodedq1dRw7r6KUO+YacyvAd5uXZJEmjB+PKxaFWYyEslD9+6w7bbw1lvFvU4+reeXEvpRZmJ3e6DB/pkiBXH33WGW2VJ0vJOyMWgQvF3kLF0+dZNfA44BVgK4+0JAi0pL8SxeDE8+GfqRqF+mNMPgwQnIaQLrPFR8OoCZ5SznixRMpgn0W9+KOyWSMoMHh+kJliwp3jXyCZrjotbzHmb2A2AiYVy6SHHcfTfsvDMMHx53SiRlBg8Oz8UsoufTEPR74D5CJ/UhwK/c/U/FS5JUtOrqMKzjpJO0hrk0WyZoFrOI3ug4RDNrCzzl7ocCE4qXDJFIZqWsk06KOyWSQjvsEKrBixk0G81pRp3LV5nZ5sVLgkiWMWNgxIjQDCrSTB06wMCBMeY0I2uAN6K5MD/rSe7uZxctVVKZZs8Oc2Zee23T7xVpQLFb0PMJmo9Fj2wlXmlYKsKYMaEe88QT406JpNjgwaFa3L041eL5BM0e7v7H7A1mdk7hkyIVzT0EzYMOCsM6RFpo0KAwCvf996FPn8KfP58uR6fm2HZagdMhlW7q1FCmUt9MaaVit6A3mNM0s5OAbwEDzWx81q7uwOLiJEcq1t13Q/v28PWvx50SSbnsoPmlLxX+/I0Vz18krPGzJfCHrO3LgdcLnxSpWLW1MHYsHHEE9OwZd2ok5fr3h44dY8hpuvs7wDvAF8xse2CQu080s85AZ0LwFGm9F16AhQtVNJeCaNMm1GsWK2jmM8vRDwgjgm6ONvUDHipOcqQijRkT1io4+ui4UyJlopjdjvJpCDqTMAv7pwDu/jYbVpAUaZ116+Dee+G448KaBSIFMHgwzJ1bnKUv8gmaa919XeaFmbVD/TSlUP797zAlzQknxJ0SKSPbbQfr14cZjwotn6D5TzO7GOhsZiOBe4FHCp8UqUgTJ0LbtnDwwXGnRMpI12gCy1WrCn/ufILmhcBHwBvAD4HHgV8WPilSkSZMgM9/PqxVIFIgmaBZjCXEmhwR5O51hPkzNYemFNaSJTB5MvzqV3GnRMpMrDlNMzvKzF41s0/M7FMzW25mnxY+KVJxJk0KwydHjow7JVJmMm2KseQ0geuB4wlL+KoBSApn4kTo1g322SfulEiZibtOcwEwXQFTCm7ChDBBR/v2cadEykzcOc3zgcfN7J/A2sxGd9ekh9JyVVUwZw6cdVbcKZEyFGtDEHAFsALoBHQofBKkIj3zTHhWfaYUQSanWYzieT5Bs5e7f6Xwl5aKNnFimOxwl13iTomUoWLmNPOp05xoZgqaUjh1dSGneeihWnFSiqJTp3BrxdUQdCbwpJmtVpcjKYjXXw/j2w49NO6USJkyC0X0uDq3a6iGFNbEieFZQVOKqGvXEhfPzWybpg7O5z0im5g4EXbdVWsBSVF16VL64vnjeRyfz3tENli7NiwVqFymFFmxcpqNFc/3aKLu0ojm2BTJ24svwurVCppSdMXKaTa23EXbwl9OKl5mKrhirHglkqXkdZoiRTFxYpgKbrPN4k6JlDkFTUm/zFRwKppLCcTRECRSWJMmhY7tCppSAnE0BH3GzNoCW2e/393fLXxypKxlpoLbd9+4UyIVoOQNQRlm9hPgUuADoC7a7MDuhU+OlLWJEzUVnJRMnDnNc4Ah7r648JeXivHuu/D223DmmXGnRCpEly6wZg3U1oYOG4WS7yTEywp3SalI//53eFZXIymRzPRwq1cX9rwN5jTN7Lzox3nAc2b2GJqEWFpqyhTo2BGGDo07JVIhOnUKz2vXhqr0QmmseJ6ZqOPd6NGBDZMQa+kLaZ4pU2CPPVSfKSWTCZpr1hT2vI2NCPo1gJmd4O73Zu8zsxMKmwwpa3V1MHUqfPvbcadEKkixgmY+dZoX5blNJLe5c+HTT2HvveNOiVSQ7OJ5ITVWp3kEcCTQ18z+L2vXZkBNYZMhZW3KlPCsoCklVPLiObAQmAwcA0zJ2r4c+GlhkyFlTY1AEoM46jSnAdPM7G7CNHA7ExqAZrv7usImQ8ralCmw++5qBJKSirNOcyQwF/g/4M/AnKjoLtI099AIpKK5lFjHjuG5lMXzjGuBL7v7HAAz2xF4DHiisEmRsjR3LixbpqApJRdnTvPDTMCMzAM+LGwypGypEUhiEkdDUMYMM3scGEeo0zwBeMXMjgdw9wcKmyQpK1OmQIcOagSSkoszaHYizHCUGTT8EdALOJoQRBU0pWFTpsBuu4XAKVJCsQVNd/9uYS8pFSPTCPTNb8adEqlAxerc3mSdppkNNrNnzGx69Hp3M/tlYZMhZWnePFi6VPWZEos4G4L+Shg2uR7A3V8HRhU2GVKW1AgkMcrUCMURNLu4+8v1tjU5jNLMOpnZy2Y2zcxmmFlmApBeZjbBzN6OnntmHXORmc0xs9lmdljzPookzpQpoUP7sGFxp0QqkFnIbcYRND+O+mZ6SIh9A1iUx3FrgYPdfQ9gOHC4mX0euBB4xt0HAc9ErzGzXQk52KHA4cAN0dpEklZTp4ZGoEwvY5ESiytongncDOxsZu8B5wI/buogD1ZEL9tHDweOBUZH20cDx0U/HwuMdfe17j4fmAPsk+8HkYRxDzlNFc0lRh07xtN6Pg841My6Am3cfXm+J49yilOAnYC/uPtLZra1uy+Kzr3IzLaK3t4X+G/W4dXRtvrnPAM4A6CDurEkV1VVWOdcQVNiVIycZj7LXdTfDuS33IW71wLDzawH8KCZNVa5ZblOkeOctwC3AHTt2lUzyCeVGoEkAUoaNNmw3MUQ4HPA+Oj10cDzzbmIuy81s+cIdZUfmFmfKJfZhw1DMquB/lmH9SNMTydplGkE2m23uFMiFaykdZru/utoyYstgb3c/Wfu/jNgb0JAa5SZ9Y5ymJhZZ+BQ4E1C8D01etupwMPRz+OBUWbW0cwGAoOA+q32khZTpoRWczUCSYw6dSrhzO1ZtgOy589cBwzI47g+wOioXrMNMM7dHzWz/wDjzOx0woJtJwC4+wwzGwfMJHRpOjMq3kvaZBqBjj8+7pRIhSt18TzjTuBlM3uQUMf4NTa0fjco6gS/Z47ti4FDGjjmCuCKPNIkSfbOO/DJJ6rPlNh16hQGpRVSPq3nV5jZE8AXo03fdfdXC5sMKStqBJKEiCunibtPBaYW9tJStqZMgXbt1AgksYurc7tI82QagTIzJojEREFTkk8jgSRBijEiSEFTCuv992HxYhg+PO6UiChoSgrMnx+ed9wx3nSIEKrWawvccVFBUwqrqio8DxgQZypEAGjbVkFTki4TNLffPtZkiICCpqRBVRVsvTV06RJ3SkQUNCUF5s9X0VwSQ0FTkq+qSkFTEqNttPZDXV3hzqmgKYVTVxfGnStoSkJkgmYhc5sKmlI4ixbB+vUKmpIY7aKB4jVNLgWZPwVNKZxMy/nAgbEmQyRDOU1JtkzHduU0JSEUNCXZMjnN7baLNRkiGQqakmxVVbDNNtC5c9wpEQEUNCXp1N1IEkZBU5JNQVMSRkFTkqu2Ft59Vy3nkiiZLkcKmpI8Cxeqj6YkTianqX6akjyaEk4SSMVzSS4FTUkgBU1JLvXRlARS0JTkqqqCbbfVCpSSKAqaklyaR1MSSEFTkkt9NCWB1OVIkqmmBhYsUNCUxFGXI0mmhQvDXamgKQmj4rkkk7obSUIpaEoyZebR1BBKSRgFTUmmqiowg/79406JyEYUNCWZMn00O3aMOyUiG1HQlGRSdyNJKHU5kmRS0JSEUpcjSZ5MH001AkkCqXguyVNdHe5I5TQlgRQ0JXnUR1MSTEFTkkdBUxJMQVOSR300JcEUNCV5qqqgXz/o0CHulIhsQl2OJHnU3UgSTF2OJHnefVdFc0msTA4zEzwLQUFTWq6uDt57T0FTEmvt2vBcyBG+CprSch99BOvWhTpNkQRS0JRkqa4Oz8ppSkIpaEqyZIKmcpqSUAqakiwKmpJwCpqSLNXV0L499O4dd0pEclLQlGRZsAD69oU2uo0kmRQ0JVmqq9UIJImmoCnJUl2t+kxJNAVNSQ53BU1JPAVNSY7Fi8MdqaApCaagKcmxYEF4VtCUBFu3LjwraEr8NBpIUiCVOU0za2tmr5rZo9HrXmY2wczejp57Zr33IjObY2azzeywYqdNWkEd2yUFMkGzffvCnbMUOc1zgFlZry8EnnH3QcAz0WvMbFdgFDAUOBy4wcwKOKGTFFR1dZjhdaut4k6JSIPWrg25TLPCnbOoQdPM+gFfBW7N2nwsMDr6eTRwXNb2se6+1t3nA3OAfYqZPmmF6mrYdtvCTlQoUmBr1xZ+UYFi5zSvB84H6rK2be3uiwCi50xWpS+wIOt91dE2SaIFC1Q0l8TL5DQLqWhB08yOAj509yn5HpJjm+c47xlmNtnMJtcUcg57aR6NBpIUSFXQBPYHjjGzKmAscLCZ/QP4wMz6AETPH0bvrway/wr7AQvrn9Tdb3H3Ee4+ol1m1SQpLXVsl5RIVdB094vcvZ+7DyA08Dzr7t8BxgOnRm87FXg4+nk8MMrMOprZQGAQ8HKx0ietsGQJrF6toCmJV4ygGUdW7SpgnJmdDrwLnADg7jPMbBwwE6gBznT3Ai68KQWjju2SEqkNmu7+HPBc9PNi4JAG3ncFcEUp0iStoD6akhKpKp5LGdNoIEkJBU1Jhurq0D9zm23iTolIoxQ0JRmqq6FPH3Vsl8RT0JRkUMd2SQkFTUkG9dGUlFDQlPhlOrarEUhSQEFT4rdsGaxcqZympIKCpsRPfTQlRRQ0JX4aDSQpoqAp8VNOU1LCXUFTkqC6Gtq0Cf00RRJs/frwrKAp8aquDiOBCrnoikgRFGNRNVDQlOZSH01JCQVNSQaNBpKUKMaa56CgKc2lnKakhHKaEr9PP4XlyzUaSFJBQVPip+5GkiIKmhI/BU1JEQVNiZ9GA0mKKGhK/KqrwQy23TbulIg0SUFT4lddDVtvDR06xJ0SkSYpaEr81N1IUkRBU+Knju2SIgqaEj/lNCVFFDQlXsuXh1nbFTQlJRQ0JV4LF4bnvn3jTYdInlavDs+dOxf2vAqakh8FTUmZ5cvDc7duhT2vgqbkJxM0NfmwpMTy5dCpE7RrV9jzKmhKfjJBUx3bJSWWL4fu3Qt/XgVNyc+iRdC1a3HuQpEiWLFCQVPitHBhyGWaxZ0SkbwopynxygRNkZRYvrzwjUCgoCn5UtCUlFFOU+LjrqApqaM6TYnPsmWhp7C6G0mKqHgu8Vm0KDwrpykpouK5xEd9NCVl3FU8lzgpaErKrFkDtbUKmhIXDaGUlCnWuHNQ0JR8LFwIm21WnDtQpAgyQVM5TYnHwoXKZUqqrFgRnhU0JR6LFqk+U1JFOU2Jlzq2S8qoTlPio9FAkkLKaUp8liwJi60oaEqKqE5T4qM+mpJCymlKfBQ0JYVUpynxUcd2SaHly8PSve3bF/7cCprSuMxkHQqakiLFGncOCprSlIULoUcP6NIl7pSI5K1Y08KBgqY0Rd2NJIWKNS0cKGhKUxQ0JYUUNCU+CpqSQqrTlHi4h4YgNQJJyqhOU+KxeDGsX6+cpqSOiucSD3Vsl5RKbdA0syoze8PMXjOzydG2XmY2wczejp57Zr3/IjObY2azzeywYqZN8qCgKSlUzPWBoDQ5zS+7+3B3HxG9vhB4xt0HAc9ErzGzXYFRwFDgcOAGM2tbgvRJQxQ0JYXWroWamvKq0zwWGB39PBo4Lmv7WHdf6+7zgTnAPjGkTzI0hFJSqJiTdUDxg6YDT5vZFDM7I9q2tbsvAoiet4q29wUWZB1bHW2TuCxcCL16hUG8IilRzGnhANoV57Sf2d/dF5rZVsAEM3uzkfdajm2+yZtC8M0EYDez1dHP7YCaVqW2tNKTXvvsV5OeNG+QtjSnLb2Q0DSfdlp45JBPejs3tKOoQdPdF0bPH5rZg4Ti9gdm1sfdF5lZH+DD6O3VQP+sw/sBC3Oc8xbglvrbzWxyVr1p4qUtvaA0l0La0gvpS3Nr01u04rmZdTWz7pmfga8A04HxwKnR204FHo5+Hg+MMrOOZjYQGAS8XKz0iYi0RDFzmlsDD1oo2rUD7nb3J83sFWCcmZ0OvAucAODuM8xsHDCTkHU+091ri5g+EZFmK1rQdPd5wB45ti8GDmngmCuAK1p4yU2K7AmXtvSC0lwKaUsvpC/NrUqvuW/S1iIiIg3QMEoRkWZIf9A0Oxyz2ZjNwezCuJOTF7MqzN7A7DWi4aWJY/Y3zD7EbHrWtl6YTcDs7ei5ZyNnKK3c6b0Ms/ei7/k1zI6MMYWbMuuP2STMZmE2A7Nzou3J/J4bTm9yv2ezTpi9jNm0KM2/jra3+DtOd/E8DLN8CxhJ6LL0CnAS7jNjTVdTzKqAEbh/HHdSGmR2ILACuAP3YdG2a4BPcL8q+gfVE/cLYkzlBrnTexmwAvffx5iyhoUud31wn0roaTKFMELuNJL4PTec3m+S1O85tER3xX0FZu2BfwHnAMfTwu847TnNfYA5uM/DfR0wljAcU1rL/Xngk3pbGxoCG7/c6U0290W4T41+Xg7MIoyCS+b33HB6k8vdcY/GCNE+ejit+I7THjTTOvTSgacxm8KG4aVpsDXREFg2HgKbZGdh9npUfE9GMTcXswHAnsBLpOF73ji9kOTv2awtZq8RBtJMwL1V33Hag2ZeQy8TaH/c9wKOAM6MipZSeDcCOwLDgUXAH+JNTgPMugH3A+fi/mncyWnSpulN9vfsXov7cMIow30wG9aa06U9aOY19DJxouGluH8IZIaXpsEHUb1Wpn7rw8bfHjP3D6I/mDrgryTxew71bPcDd+H+QLQ1ud9zrvSm4XsGcF8KPEeYerLF33Hag+YrwCDMBmLWgTAf5/iY09Q4s65RJXr4ecPw0jRoaAhsMmX+KIKvkbTvOTRS3AbMwv3arD3J/J4bSm+Sv2ez3pj1iH7uDBwKvEkrvuN0t54DUfeG64G2wN8Io4qSy2wHQu4SouGliUyz2RjgIGBL4APgUuAhYBywHZkhsO7JaHzJnd6DCEVGB6qAH35Wj5UEZgcALwBvAHXR1osJ9YTJ+54bTu9JJPV7Ntud0NDTlpBJHIf7bzDbghZ+x+kPmiIiJZT24rmISEkpaIqINIOCpohIMyhoiog0g4KmiEgzKGhK6pjZFmb2WvR438zey3rtWT+/ZtHMV2b2nJm9a7ZhlTgze8jMVkQ/DzCz1dExM83sJjPT34dsotirUYoUXDT7/3AAi2Yy8miGHTNb4WHIXC5Lgf2Bf1no8Fx/Qfe57j7czNoBzxImcXgAkSz6TyqVZCxh1BiEqcFyBkR3rwFeBHYqUbokRRQ0pdx0rlc8PzFr3zPAgRbmYR0F3JPrBGbWhbCO1RvFT66kjYrnUm5WN1I8ryVMQnsi0Nndq7KqOAF2tDCFmAMPu/sTxU2qpJGCplSasYSx/5fl2De3kYArAqh4LpXnBeB3wJi4EyLppJymlJvOURE740l3/2zBPQ8z1CRvLRtJDc1yJCLSDCqei4g0g4KmiEgzKGiKiDSDgqaISDMoaIqINIOCpohIMyhoiog0g4KmiEgz/H/WHvGCJEs8RwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig,ax = plt.subplots(figsize=(5, 6))\n", "t.mean(dim=['LAT','LON']).plot(ax=ax,y='DEPTH',yincrease=False,ylim=(500,0),color='r')\n", "ax2 = ax.twiny()\n", "s.mean(dim=['LAT','LON']).plot(ax=ax2,y='DEPTH',yincrease=False,ylim=(500,0),color='b')\n", "for ticklabel in ax.get_xticklabels():\n", " ticklabel.set_color('r')\n", "for ticklabel in ax2.get_xticklabels():\n", " ticklabel.set_color('b')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Further resources\n", "* Documentation [page](https://docs.dask.org/en/latest/)\n", "* Dask+Xarray [tutorial](https://github.com/pangeo-data/pangeo-tutorial-sea-2018/tree/master/notebooks)\n", "* Dask [tutorial](https://github.com/dask/dask-tutorial)\n", "* [Example](https://pangeo.io/use_cases/physical-oceanography/sea-surface-height.html) of cloud computing using dask\n", "* Follow Matthew Rocklin's [blog](http://matthewrocklin.com/blog/)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }