Skip to content

Instantly share code, notes, and snippets.

@ccarouge
Created December 2, 2019 02:55
Show Gist options
  • Save ccarouge/6c0eef4c4e360ef1820a003b606caddb to your computer and use it in GitHub Desktop.
Save ccarouge/6c0eef4c4e360ef1820a003b606caddb to your computer and use it in GitHub Desktop.
Upsample data and fill with 0.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import os\n",
"import xarray as xr\n",
"import numpy as np\n",
"import sys\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"year = '1995'\n",
"fpath='/g/data/w35/ccc561/scratch/martin_763'"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"fname = f\"{fpath}/AWAP.Rainf.3hr.{year}.nc\"\n",
"ds = xr.open_dataset(fname)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'Rainf' (time: 8, lat: 681, lon: 841)>\n",
"[4581768 values with dtype=float32]\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1995-01-01 ... 1995-01-01T21:00:00\n",
" * lat (lat) float32 -44.0 -43.95 -43.9 ... -10.099998 -10.049999 -10.0\n",
" * lon (lon) float32 112.0 112.05 112.1 112.15 ... 153.9 153.95 154.0\n",
"Attributes:\n",
" long_name: Rainfall rate\n",
" standard_name: rainfall_flux\n",
" alma_name: Rainf\n",
" units: kg m-2 s-1\n",
" _fillvalue: -999.0"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"rain = ds.Rainf.sel(time='1995-01-01') # Limit time for test"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x7ff660de28d0>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEXCAYAAAByAUkhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzsvXu8NFlZ3/t9VlVf9uW9zp2ZgWFgiDqIRkHQeIkBBfmoGA0eAwfP8RLEy9HEaMRLvH1QI+ecxKgIYoJIDmiiiBAlgiSSoOMoF8fR0RmYGWCYG3N5573s3Xt3d9V6zh9rrapV1dW9q/fu/b57v2//Pp/au+u2anV19/rV8zy/51miqiyxxBJLLLHE+YS50B1YYoklllji0sOSfJZYYoklljjvWJLPEkssscQS5x1L8lliiSWWWOK8Y0k+SyyxxBJLnHcsyWeJJZZYYonzjiX5HECIyJNFZENEkgvdlyWWWGKJ/cCSfA4AROQTIvKCsK6q96nquqrmF7JfsyAi3yMiHxKRoYi8uWH/t4vI3Z5E/1BEnhTtOy4ivyEij/jlJ2vnfkJEtvy5GyLy3h36coOI/LGIDETkzvheisg1IvIuEXlQRFREbmjx3l4mIp8UkU0R+T0RORnt+0YRucVf6/3R9i+J+rvpr7URLU+edm7Uhvpzwzn/YYd+/gsReVhEzojIm0SkF+07KSLv8O19UkRetkNbz/f3buDv5VOifSIiPy8ij/vltSIiO93HJZaYhSX5LLFbPAi8BnhTfYeIfBnws8BLgJPAx4HfjA75d8AqcAPwBcArRORbas18jSfgdVX9yh368pvAXwKXAT8K/I6IXOH3WeAPgW9o86ZE5GbgV4FXAFcBA+BXokNOAb8A/Jv4PFX9QOgvcLPffDx6D/dNO7eGz4nO+fYZ/Xwh8Grg+bj7eCPwU9EhrwNG/j28HHi9f29NbV0O/C7wr3Gf14eA/xwd8krg64DPAZ4FfDXwHTPewxJL7AxVXS4XcAH+E26A3AI2gH+FG0wUSP0x78cN9Lf4Y/4rbqB9K3AW+CBwQ9TmZwB/hBvs7gK+cR/7/xrgzbVt/w/wumj9Sf79PM2vPwY8J9r/I8AHovVPAC9oef1nAEPgSLTtA8Craselvg837NDezwJvi9afhhvEj9SO+3bg/VPaqHx+Dfsbz/XnPL3l+34b8LPR+vOBh/3rNd/nZ9S+Z/9mSluvBG6J1tf89/Ez/PotwCuj/d8G3Hohfi/L5eJZlpbPBYaqvgK4j/JJ/7VTDv0m3NP4tbgB8c+AX8c9qf4d8BMAIrKGI563AVcC/xT4lRlPvb8iIqenLLfv8m2JX+J1gGc2bAuv430AbxWRR0XkvSLyOTOudTNwr6qei7b9FaX1MS9u9ucDoKr34AfyXbY3L/6Xd6X9buwiFJEvFpHT0/rpX18lIpfh+pqr6kdr+4t74j/fL25qS1U3gXui45uutdv7u8QSwNLtdpjw66p6j6qeAf4bcI+qvk9VM+C3gb/vj/tq4BOq+uuqmqnqR4C3A/+kqVFV/S5VPT5ledYu+/pu4BtF5FkisgL8OO6pftXv/0Pg1SJyRESeDnxrtA+cm+gG4CnAHwPvEZHjU661DpypbTsDHNll3xfd3jz4Mtz7/gycW/P3RSQFUNU/UdX4HtT7GV4fadgX9hfvwX++fzKlrfrxTddaX8Z9ltgLluRzePDp6PVWw/q6f/0U4LmxBYMbzK8+P90EVf3vOEvs7cAncW60c8D9/pDv9X3+GPBOXMzm/uj8P1XVLVUdqOrPAaeBLwEQkTuigPyX4NyQR2tdOOqvNxM1kcAdfvOu29srVPV/qepIVU8D3wc8FfjMKYfX+xlen2vYF/ZPew87Hd90rQ1VXVYlXmLXWJLPwcAif8SfAv5nzYJZV9XvbDpYRN5QU2VtNAzIc0NVX6eqN6nqlTgSSoG/8ftOqerLVfVqVb0Z9z38i1nN4d10qnqzlgH5DwB3ADeKSGyZfI7fvlMfPxC1FdxId/jzARCRG4Ee8NGmNvYZxftuQKWf/vWnVfVxXF9TEbmptn/aPam/5zWca/eOpv07tLXEEq2wJJ+DgU/j1EqLwO8DzxCRV4hIxy/PEZHGJ2hVfVU0ANeXqX59EUlFpA8kQCIi/eAi8q+f6SW6TwbeCPx7VX3C73+aiFwmIomIfBUu4P0av+/JIvIPRKTr2/lB4HLgT6f0/6PAbcBP+OP/MU6R9faor30cgQD0/Po0vBX4Gm8VrQE/DfxuiCn5PvdxZGr8NTsz2ovv2dRzReRmEflcf8w68P8CD+DieU14C/BtIvJZInIC+DHgzf6ebOLUaz8tImsi8g9wysP/NKWtdwDPFJFv8P37ceB2Vb0zutb3i8i14iTz/zJca4kldo0LrXhYLgpuYLgP5176AZrVbt8eHV9RmAEvAO6O1v8e8AfAo8DjwP8APnfBff5J38d4+Um/7zhwO7AJPAz8HJBE534jLqYxwBHHC6N9N0fnPg78d+DZO/TlBn+PtnDqvhfU9tf7qTu09zL/eWzi3IIno33/Z0N7b27oz4Tabda5wD/yfd8EHgF+D7gpOvdLcK6uuL3vxz24nMWJT3rRvpO+jU3/Xl5WO3cD+JLad+hOfw/fT1U9KcBrcerJU/61XOjfzXI53IuoLt22SyyxxBJLnF8s3W5LLLHEEkucd1ww8hGRl3rlkhWRZ9f2/bC40ix3+UzuJZZYYoklLiKkF/DafwN8Pa6USQER+SxcQuXNuMz494nIM/QA1zlbYoklllhiPlwwy0dV/05V72rY9RLgt1R1qKofB+7G1f9aYokllljiIsGFtHym4Vrg1mj9fr9tJrrS0z5r+9apAwuRaiKIFH9Aov9he3htDCoCBtSIWxLcYsq2NGqO+mtAm7JQ2uS9Tzumhf6l96nNFhcoMbxuDYnbrV8jWpf6/oZjO480Xz+7cm12/zV62zP6UPTDIxlk5Ur4rIwgVhGrkFmwblFV8CIiKT53sOsrU/qh1fcfBEhBJKhU15v2hXaItqsWr7Wy/2DgHE88pqpX7HzkdLzwy9f08VPtHDIfvn34HlV90V6ud7FhX8lHRN5Hc2b9j6rqO6ed1rCt8ZsrIq/E5YjQZ5XnyvN31c8DDTHRS5nYV2wTA8aRixjjSCZJkCSBxECSuKXbgU6K7XfRfkreT8lWE8brCaMjwvCoYLslCalxrxFHNPFrjJbk46u5xcTl/lePUXHbALciWv1fnuaPqa4/7V/8RTt7XW3x8t5//oWIBawg1rcV/mvDunUL0Tq4/1f/u1umEuej3/hFDf1w7brrh9davU79lhR9UUTh2F+dcjtTg13tMjrWRTuG9NyY9NwQ88QmDAboYAu7PUSzMQCm23WfOcBI2f7yz/aE5UhGcgVPYJIpYq3fZpHMOsLI1RFS7tfDa5uD1WhbXuzTPIcsc/+tuv9qURveZPnZXCi8T3/nk3tt4/FTOX/xnie3Oja55mOX7/V6Fxv2lXxU9QU7HzWB+4Hro/XrcDkhTe2/EZfAyFE5ebAerfaKHUinsl2MIxhwxJMkjoQ6KRgDaQpp4v53UrSbVolnLWG0LozX3ZOyWEciYTAmrFN9rSqIZ4eCNxpy8mNecW1KdLwU/2MeCqTz9H8+q/BBAxoGtht/8M8A+MTPflEz8TT8J14nIp4ZuOJ1t/DodzUQkH+bIuX7Krb7+1xYmBr/FxR1ZGsMttdheKLHeD2hs7nzE7cdjRwBSY2t434E60hm/Hxiq0V18h7XrRrjHnywjsHFCGoNYjwBiTkQBLRXuK/S4X8fFwoH0e32LuBtIvJvcYKDm5hdeuXiQm2gmMvaCaSTJG57GsgngU4H7aTQSbD9lHwlJVtxxDNeF7LVQAIUT+XxgFgQSNhnqPrkgquHaHytWwhabq+PdYGDbvq/bp0cLNugxWB2w4848rjvp76oSjjeyolJL+7zk37hz1sPmlf8yi089p01App4r+KsiQoJ1I/x/fGuUzWGfDVleCLBjP0pnkjdAYpaO9E/tUplPlyR6S6wJnYsGtLJ85q2RW1Jkjg3oDVAjSwvAgJSlPFSB7VrXDDy8WVQfgm4AvgDEblNVV+oqneIyH8B/hbIgO++JJRuswbcJkvHvxZPOHjCkTR168HaCaSTuidn7SbkgXTWDOM1gw0z3dho/AnjQsRJAi5GRDTwSXAVudGysIqKk7w1hU61jlRqpBMGpR2IeLdunCf/hCOh+3/si6puvToZWbj2lz/iXVd5eb0dcPnrb+Hx7/iiqrUX+lyQCsV9iFstNnuyP3H7E/7zFfKVFJsKZhzHVYDEoNtD5/Kq9VHz3N23JKluD5+L+s8lkeLzqsSDpqHpXkxYQMYREDlCguaU1g9cFAS0tHx2jwtGPqr6DlxNqaZ9PwP8zPnt0QVEA/HERNNEOoWFE0gnSZyVEyygwsWWoKlzs9luQt5PyPuG0ZGE3Fc7q7ua1PpBUMqYBOLHCe+lUevHDuIH+PKpWiYGVpkkIH/KTd/7F80DUURCExZg/f7tYhC77jW3gBge+OHnTRDP9W+8wwXwE+Me2kX9wNnuWpf9qiOgWajcXyIXn793YuH0Z5/g+O0+5lMb3EPshtyiWYYdjZr7FmJ/8XXr3CF+owFsHHurudxCnCeGtTtbQOSI1txvcKgJSFHyAyakOEw4iG63SwvTiKe+PQweswQFwcXmrR7tddBOgu0m2F5C3kvIVgzZqrd2oEI6xQCMG0vERgOk4AcmSuvGW0qxtVNYQF6AsBCo9YNW9b5UBrHdwLd13Ws/WFgHBaGLuCf3og6VGzgn3Ecz+z1lve5uCwZcTEbBzQmEwH8ytIUYIQgSUIXxGElTJMnQrDqQS3gw6XZY/atPsfWs66eHd+LpeeLXe7nHcXtiCgK/WAjItpFnLtGIZXmdCwUxjW6leIAtBg4xxRNkMTh2Uvc6EhLQ60K3i3Y72H4Xu9IhX+2QrXUYr6eM1xOyVVMRE4iX2oYBLSiw4vXCJRe5o+oxk1KhBUE8EM4JSq+wPR78nv59Ppy3m8GnbhUuGlIjJKY8GEzBZW+8pUrsNFgcE9es/Qd3H/OcZDujd8aRX5BZy9gpyxqFAIB0u8hKH1lfg2Mz5sPzT/AaHmym7K+sz4r51PthQnzShA1Va3Y/Pr99hgI52mpZYhJLy+dCoIW1U6xPExTE8unElG62TuKtnRTbd9ZO3jfkXZfHU7mmRvGFYN3Eg6WtjoFF/DseIL3wQE1pITm5tXezER2r/oUfgedSsgXrJ5nc5dRUu/iBB99hEwLxW3UCAfCyYZjH+qkpyKccRIXMXd/K88VasGAGI3qPCaPjXSeNHucwztDxGM1zJ2muv40j63BkHXt0hexoH4qHDS1capPiD5lutKqd3xIKIgefV+biP3lkCR3ewXlp+eweS/I535gWRI9FBeGYIJ9uiu3UrB7tJGiSoD3vZuskZGuJIx0jZf4NVAkmDHzx9vAAnZTxn7LDFBs0xH8kctHFQgIjbtxJtJQTQ6XB1sQRixH8U3OhQCO44BbkvhFBOp3S2sky1IjLf1HrCKitq0iZHMinueMqfaB0vWXe2sktSWbpD8doNy2tnnGGjrPa+Yb08pPoZSfIj68wOtZleDJl+4Rw4q6Rt0SpWi7h9TQjZDckYYyXXMfbxLsxPfnFn+Uhcr8pMF7GfHaNJfmcL+wgKijWd0oWTaPcncSg3Q4U1o5b8p7B9gw2LdmglEdrhSFC7EC9OiDEHOJYT0FA9TiFif5LtARSo7SI1JTbUOGmf/kXpWLP54E0Djwz3DHh/pUklO9uAIutKhGk24GVFUfu4klne+gO3VIkb698O/kfb+HUt1WFB9IwYE3IvOPXWUksMh4jW/47YhVGYzQkdIZj0g7JZSfQK08yvGqd0dGE4XHD8KgwPgb3PaXDDe8eEpJYUY2SXb37TcRVSFj0k32wgpIEIfdxNDiMOUC6dKntCUvyOR+YJRkW00w6RdDby6eD9VMkjHrC8fJp2zXkvQTb8VZOIJga2UAcm6GQRAsKVgqhgPEJ6/iuFa+9i42gdgvkElxwQXwQl+hRVwqGyvY6qTTkgYR7NY2A/D2TBJepb8S5n+axgLQMgGPVEXsM499Yv+cG5sLKaE90Yp35U+T3zDpWG9ZH40lFmVeeaZ776gJlP8z6GnryONvXrHP6pi7jdchWIe8rtgs2Ve75ppSnv3VcEo/6Uj2qziqqayuKsju1//XXlc6X6sfwvVZrqRCQ9/lOJKHCwSchhXzJPbvGknz2G3sgnsLa6XQc4XihgabGVSnopo5wegk2EZer0dSFyOUVJ4vG7raShCjyd+rWTVDBFcQTBik/PgcLJ1g5hUoO54ILT9qNXiYjFNlcUlW2TZVZx8Ssio7H/nJzEhDeespzpxzrdvy5PsAUl5fZLQLBx+t1a6fpNMGRD6Dh/RQkEUmsI+SnT5MeXSdfOcHGk5Xs6hHdlTGaG+wwgdygufCxV3R4xptHjriDeCDItyudCIRja+uuDzMTV2soCKi68VCq4GIP9RLzY0k++4kZ8Z3GKgVNsZ00LcknTdBuiu2maMc4+XR3Up1UIZgZwe5K1QJwg5oASBG/KUrt5DUXml8PllIgnFCc1PrXMZGJH9+e/uoPTfRlwo0V36M4QdJqdUBKEldGCEDEZ9Qr80miI1+jKjrOkCx39z63LuYyHhcurnkHxKJu26z9Wh4Xiw2O/+GdUTKtFvGTws3WIDIAoNslWxF4yoDveOYtvGD9b/nU+ATveOzz+cjD17F5ZgUGCR97eY9nvGWrJB5vYUk+aTHPheAWhOa4TyCscJyUplYlDnigCUjIF5ZPcOnh8OkbDwt2EBYUr2vEI0ninuYD6XRcIVDtd7ArYXF12Wxa+/jiATxWMoXXUSmWeKCLpdNiQXLF5IrkjnRM5pexW5JR9Hrs9oX/klORWMfybTe4l/ejYtFMSygV4wbYYiDzFtE+SHML6ycPgfxxuYxG6GhUHexbDorH3+xqy+0os57ZOS1db6rOQoaJygXgCopqLyXvCMYoVg3HzYivXTvH665/L1/31NtZPbrlPvNM+NjLVpFcC+LBUlTJbt23aWjIHZJG12bkYj0kMmzFfS3bLEtMYmn57Ad2UrSFAaPI44lEBcG9FivZ4mTRblJYM85V4gTNEqyV4NYKr+sPZrE15A2dSiJpvFtKkiraNuW6GimsHjVgfWWfIMcOVlHRtnWtT7jRmvJTgpjAxzOcGMBbIjWJ9I5WQFtoWZGZ8djXXxPIcnQ0QodDH1OafzSpFGuduG6z1XPijz7mVo34OEzVggguLEmSquCg20XTFNsBMcq27fBo3qfDJn0Rvu3kn3HnxtV8+NE1zJbxyjctK1gHgsu1JBYb/d8pbiXGuQibrB9v8RTut0JhJ1CL+xx0F5wCo+Xz+66xJJ9Fom2ZnIa6bBMVqL2rTbsptp+iaYImwbcecmhkgkyK69aIp15ZOhBWcX4UBwrHuDZjC6q8lgbSScAmuJpj3tUmFsj9eBJccnEcqYYml1vYHsjHrK46It7YxI4zZyGatHxSh8bg+7xQqzDOQIbQcU//mufoaFz0ZddP4uFt1j6X+v5im0S5RsbfwIqWXZ0irUY+pE56n3eFxFg28h6fGl9Grk9w1Awxonz5yTv5kLkBLCQjccQTW1aBePLIApoQPLS8zzsVM433yeEiILtjEtcS07Akn31EI/EEMUFdWBAk1IF0wtQHvQRNTDVPJ8DHxCtEY7UqPKhbOrHs2T/ZS/AfEMUfCleZFu4zCFaOeEvHE6GAyRUrfr6c3LWv0etAbDf+xGS8Z+qAohYxgllZR66+Era2sWc3vIXi76+17j7meTGXTFOyZWuEtjNx8R3Ajmtxnl0QkHhxxkRtu6n90NKlZi2CjwGGeXM8MasXI1RgnPgk64OqsJn1uHd0BSNNuCzdoC9jbuw94ry+Y8GM4L6vOc5Tfu/xSeIp6rnZVlZPsHYK66eyb4r1c0gJSGEZ89kDljbjIhH9KKbGeCr5O8ZZPIWSrYzv6EoH7aeu3AlRvGYXiCcmq5TD8fEdybUS10lG6pahJd2yJONoQAoxpNBOrT2TgfGkIzkFGRWkNOOeAd6iKS3DwkV55iyECs1RkqkOh+WSjfdGPFGf7NipyKYW65xjADz3si/050QW5ER77l8coyvq9XU60O/B6gqs9KHXc5MCGimIUNKOb8AgaeqqVoOTMgOPjI7y0PgEj2ZHOZ2vcjpfQ3MpY3kjnPVTt3iCm7OttLqOqBCuWzc7Hxuai9YPYgxIEXJMq2UniMiLROQuEblbRF7dsF9E5Bf9/ttF5PN2OldEXioid4iIFZFnR9tfLiK3RYsVkc/1+97v2wr7rtzzjZqCpeWzSDRMfTAhpfYDivgk0SJ/p9ctJ3rrJG7wUPVGg/+bA0koWuPiLa4kfji2FvcJ26CI7QAFcdQHPJPj64X5nA/xeT8+h8a5t4rWiv8GcV41cceFuI/Jq1LUckbT8um14nIL9ysqogq42Tk3B+VhcYVmq2Czybb2ikU+Xau6gquGymcRq9xiN2fxOk2Lidnseg+70sUMM2Tkl2wFWVlx92ZrG83G7v6pIpll5XHlkUfWuK1/LdcfPc2p0RqP9da5vLPBqWyN9ExCuukEJJLhy/ho6VILFk88i2lAk8utEr+JrJ9C0RYUbqUFFMQHRfwnVr95yxcS1OqBtIAW4XYTkQR4HfAVuMk0Pygi71LVv40O+yrc3GY3Ac8FXg88d4dz/wb4euBX4+up6luBt/prfzbwTlW9LTrk5ara4KJYLJbksyg0zbkT/sc5PMZ40jFF3g5pinY7aDedmqtToEIyUWUCKOXVoUs+thMnnBYTlNWsIZN799rYNgsVvDsvTlBd/Z1b2Xzp81DjBteqRUUZogivm8aK2FqMKkoXr/McZVxze0l5f+skdj4GpTmv4Vyd/qZGcbcCtc+s2OYrLOhKl/HxFcZHEpLtlGRoSbYzZDvDGIOMx+jQV2Cwim5tkzyxwdrDPfoPdHlkfBmPX77G1SfPMljrstHv8ZFHr6dzVki3ndUjOTz2hVdw+QceLl1t8VTZRV919hQK9ffeJD6YfnAjAVUOOUAybEUY1Qsm7g5fANytqvcCiMhvAS/BzWkW8BLgLaqqwK0iclxErgFumHauqv6d3zbr2v8U+M1FvIl5sSSfRWAK8UjkcihiPIkppdRBYOBL5LiETh97AdQKmIg8YoSSAxbElNaPKwfjs/19Q3EB0SLWE36/1pNO5uTVwKRQwbfT+/3JQqBrv31rZf2Jb/1CJJdCbCDhzeTw5Nd+sHQVwUQ5FanLzjtpMRnaxLCl0VPw+RyA4oeKeQgo+gwah+DKAO//+fwl208ZH0nYPpGQDg3JtpJsJ6SDhNSA2Yx+xmqxW1uYJ07Te7DDZXec4Oy5hO1zqzywlTK6KmGQdXj4wRMcPQfpAMxYScaQbtlJ4mnq39Q3GcVuIrKpxH8arJ82bcclmA4KASlg20cuLheR2Jp4o6q+0b++FvhUtO9+nHUTo+mYa1ueOwv/G46sYvy6iOTA24HXeMJbOJbks1fU/c8NddoqxFOUyakWBQXKsiZQcbkVFk5D1qioOpIKA30Y3IKrLMQQhNLiCe2HeE2mZfmbwjVWvU73v32w1e048aY/c7XMggVUK9OioV/125h2qtMX1P3/aaeI51Tca/sx8DSVd6lVXQglYdqV16GQWQciL/bFbjctjwecFSyCpgm2YxivCXlXSLtK2nGiApMpplP9GWueY8+cJRHDejfFZOucG6Wc63Q4d7TPMy9/iDtGhmQb0i11cb2Rkm7l5JcfIXliAFu+nlwkva6+qYg04hJOUwiost5AOBPig2nWjxwsAppDcPCYqj57yr6mRuoD/rRj2pzbfFGR5wIDVf2baPPLVfUBETmCI59XAG9p0968WJLPXtA0BUL8upLHY0pVWzzTaGqKH29BJLG1YwETEZF34RQTFoTAsieoUEcMmcz/iX09Tg7t54QJY0gD8XTe0450Ypz8j7fw+D/7IpegCk6ODYTgeJPFIp3UzT0T4hwhcBTFd5pEHJqz88CzqFphdfFI2yoKEeEX6sJ4XyTggIh8QnDe112zXXcvMf4zzYWkl5B6OXblknmOPXcOc+YI3aN90ssTxMKRlW3+4fG7OHXzGp+442mk20pn05IOMpKNIebcNmwOsGfOgi9Ymlw1JeY8jYBitFG/7YRK/IcDYwGpCnmjDHVu3A9cH61fBzzY8phui3On4ZuoudxU9QH//5yIvA3nElySz4HBjHptzcSTVImn61Vtcab3tB9vHT54XRCQRsRjiuwdRzzWlZspXF9xn0NiIRR11+LrN4kD5sFlv3ZL7YLGj8ENFQKCCrDbcfcICpWXqKIDCjdbURkhyKtJgFqeS/1Jea+oZdyHQW9mNe5wbbXFxHrFNAmUr+uT9sXKRPXxP7FKOrCkmwnjNaADkokrIpu4wV2zScm1HY0wVr00Xsj7livXNnhS+gQ/cN0f8qq178aMLenmmPTsENnchsEWOtgqiEetkj34EADpNVc3vMeG72xc1BZqyaYNsZ8JqXXN+pl1by8wAdnFSK0/CNwkIk8FHsCRwstqx7wL+B4f03kucEZVHxKRR1ucOwERMcBLgS+NtqXAcVV9TEQ6wFcD79vzu5uCJfnMi53qtUUCA8JUzEEyW0yFkFR/tNEPrxLrCeuuNHRV+RYIyFJWKAhGuMFfX0vSafKtx2RjymB4QO+Pbmtnv8+LqYN1ZOnEcmOoDDCysoJ0O6669GDgtnXSosJBODbEi4ptCxicSqm3d7vtZP0UVhfeuvGDahRzq+yDQsK++ZQjrH38HDLOSTfHrD1i2LgmcVOgi1tE1ZX+aZCYS+pyxrK1hGwV9EjGNatnsQgf2b6B8ZqP95wbIRtbsDlAt4eumkOo5BDdr+yhhwFPQm0TTHeDtiGGC0xACq1k1Du2o5qJyPcA78E9Tb1JVe8QkVf5/W8A3g28GLgbGADfMutcABH5x8AvAVcAfyAit6nqC/1lvxS4PwgVPHrAezzxJDji+bU9v8EpWJJPW8xbnbouqe50yuoFUJKJ9e0E0e4vAAAgAElEQVREopmK+w1cUU+Rwv2GUhT7LMmmONlV+5eSTDSal0UTt11NOX1CcOPEWP2Tj6JJMldsYyZmna/WPWmPxm4StzRxxTwHA5dEmSTOPWgE6fWQtVXodqufQZxxHxIxYTF5P1B8/hOS3x2O3/q651RdoaYkmanJvCEep7B1/TqrnzhDmlnMdk73iZRsNcX2DCZT0nMj7OkzzV1IXLIyAu4jFB4aHOXnzryY+x8/Tncb0s0Mc2YTNjcd8QyHTjFXL+AaIXvo4aoVFBNRkwUDzhNgG1xvbVEXplB9oLgQBKQIY13MEKqq78YRTLztDdFrBb677bl++zuAd0w55/3A82rbNoHPn7Pru8aSfNpglutmBvEUkuqgbkuidqyW0xQ0uR0C6uo3tJi+ACKrR8rjiMnG4PplpEo8fntRlTq65PrfPubcgiMD+RyxjT1ArXuCl26n0ZUjnbSwJMmyMn7S7ZQuw9xNMSCqBelUBqV51GnBpUaDW7XpuCmuvsIyhcn6+4XLLZpTR8vtYbIYGedgM5LhGLNp6HghAplFNgdkgwGNMIJkOZ2zGb0zCauf6HDn6RvonBOOPgwn7tym++lzjuS3tp3FE5U0ao3YfTyHDLu8D4o2lfBpm7d1AQkoX5bX2TWW5DMNUwin0eIp1qu5PBVlG0wnmIBG9xsVAgLcdAbB2xc7xgzlAGdAk9gkciQX/1ZCfba4L8lWjvY7yKDjSTQD3SG2sQBoniN54qyf4JaMVVaerHU8LkvqhLhapwPdrnPBjSwT88XA/H0XP7PpRDM6c33iexMl+hbJpvG+wvrxBOTdcLEFtH39cfp3urhLXPiTPEeH1bl8JrC1TedUhxNjy/qDHfKekA4s3VNbmE8/gW4OnLXjp+FuaylmDz1M+qRrKu+lvAcNLt4mImmaIG/a8VOKz9bJ5nwSUKhwsMTucMHIR0ReCvwk8JnAF4SMWhG5Afg74C5/6K2q+qrz17HpX6bpeTyRpDpWtgWJNfgfQVL7gVFaPzWiAapuORMNxGFGzHpbuPYwxseGIMSeQrmV4PrRRDj75DQElTBj6GwpPatoN3WWRif1E5ntv+UDOBffaFSQSlEotKhaDYhB/cArxqCpOtdmljn3U1EPLRqA5um/uM+umCOo5n6aKvOe9rASqRMLZVu0r+Juq1lAhXvO4ioYxJPJhRp20RTbE8hzdGMTsZb0iXPux56mMBrBcIQdDIp2diVfr0+XEPJ3oNEC0rg+nN8/YfHU+7EDgexIQPsMuxi12yWJC2n5NJZ+8LhHVT/3vPamDelExzVNBlco2yQSF0BZ7LMINmv5w43db00EVLiU4h5FxBMRkpNXR0QTIzKCNq7tYDt4F5xvLoG8A3lHsN0U03GWjxYS8v0loMJKBBd7qLy3cgByhOQHHTHVQU6MIyofq2qtSis6URKP9Hrl9jDQ57aM9+jOcZ/Ri11aR1FFfIrVUycdFFdjjWDpgN3aKs4DClK04+nko3nuSGs0cuQ8zjD9XnFuk6hgHmQPPEh63bXNVk8ogur7OpV4au9n8k20yKOaRUD7aP0sSnBwqeKCkU/L0g8XFBPJkA15PMCkpDrx5BMTgypgXeQ3EFCAxZeu0UYCanxtKeu6eddZUUyySEoNx0thZQ2uKj/yeMprC0ju5bsdH6sqJufZX0jacYOiJ2sdjZzQQK3L/VnpuwOzrLCE3GCTI2opErB9HTDA1UPbHpZPz236kSRIt4NZW4W1NXcPVGE4QkIwfpzRmoijShPOwAyEQiWuU1g56snGRlZPkF3Xq1er3ZE4KhaNP64gsUUOyk0JphXpfvxdrxFP9B5mkfnUqdRr+xsJZ58ISJFlzGcPOKgxn6eKyF8CZ4EfU9UP7NuVWsZ2Ktsj4pmUVAeLxtYCsfjqBYGAKM+H4sdaEBAUardwvoQfdu5/SCbk9IgfpCLSKTrs1vOeIVtLiviCShAf+Mv7osh51wkTNDWVAp77hTDgy+qKI5/4Papzv8m6J4IsR0YjdGu7kBerNQh5daoKkWrQXK2vkEDzIBQsnm4Hs74GJ44xvO4YeT9BMqV7ekjy+AZy5iyab0LWQu1Wh0bRt2DZeAJy8R4tP+Mo3oMq8qe3TynJ02JAnTFdxUKQRA9Z4EQScQXryPoJKsSCeBrIsRFRYnIbEpo6FcNO15kTqixM7XYpYl/vnIi8D2jITONHVfWdU057CHiyqj4uIp8P/J6I3KyqZxvafyXwSoA+q3N2bhekE7bFhUKDAqtuIQT3QlMJea0REFRK68TCgkp7obpwaM9bTI3wIabOw2exR1dQ08N2BdMx5B2fZZ9BkpUEZDKwHdz03MGKMqZ54FsExCDdLmZlxU0ZYBLIMhdvyVxdN9IU+j10pevUXVvDwlrW7aEPkCdIKs5KCg8AeY5RLQLpRXXkOgHFxLO2CieOMbr6KBtP6jJec9fpP5FwBEgGzmpoLd8OcZzw3Qk30laJxxcJK4inYvVYFjpgLhrZJz9F+hSfYK/qyCjqv/t9lG5jjaZpqLgvpzwUuH+l9TqtPFP1tCkEFNpc2P2URSWZXpLYV/JR1Rfs4pwhMPSvPywi9wDPACZKfPvCfG8EOConZ4+RbcimdtwEEYXYDlSJp+iPRaZZNaruqbCQW08hIHdCtU+xiwKqku1wXuKvc+cnQASTJEi/B70eJjGkiaBJz2W6dwSGik2liP2E31AyCutVebEuOuQTBv2isKo3vay/x2laBv4T4+rfddMiPibqZ94ceeup30P6fad6AyRzs53qOCtiNpKHuFH5vvAJqtLvwcnjbF9/nMHVHYbHBU18xWcfj5lXQtz9ww8yetFzGpRbVeKRXCuEU5LOvtH+YhHK+7jAZam6jMoDFQQU8q/axM0Kd5ypEFAFTdtqMR8xNbeemImf2G7gvKbLmM9uceBsRhG5AjilqrmI3Iibv+LeHU6b0WBLIUF0bCPpuM75TbW8hhDfsaAmEJCU1k84VyLRAUxaQBU/OeWAF6oMhwoA9eMAvfdTlfM0z2E0dpySGEwipInBpkLSF0SlTFT1cBOLhSfx6tNh66D9HCjK5PgcHTdA1ay7cM/Cew7S9U7HJaTmtqwU3nVzImGkkLcLuG3jDB2NkKA6jK4vaYr0+9i1HtlaQt6lmBgvHSidgcVsjcuYU/s3SPc9H2b0wjJvL9R1i2M8sXKxLLXjP4db/2pP9/h8ILv3E6Q33uBWcvxDGaULTsTVdis8BHEMs53rsJGApv22I4tpaiLqgrAUHOweF1JqPa30w5cCPy0iIbL7KlU9tbuL7MG1BhOVlQPUWkdAYUCsxXcaLaAGibXb3mABuUbKYGwYmGFSMacKDz1SSpOh9KVnGToeY6yLBCXGkHYMtitkK8YVTbagxhUkTYaQDJVkmCOZnW+g3S2ClDorqygXEutwj8cZMkomC1F6l04Rm8oyn1CbFtaTO98F68UYNElc/k5wmwbrK01QY5Bc6QyUdNtZJOmWpXNmhAyc2KD14BV/r+JTaoq2IsYTjjss1k4d4bsYau6FB4egfPPWhxiD5raxFt1MRA8NE6j/TuNSTPtIQIosZDK5SxUXUu3WWPpBVd+OK+W9e8zpYmuc8noHv3KBOgFZd+4EAU0VIPgVW5g//l/N1Tbt2qfOFMo0UUXD1M+hpdEIu7VNkmUI0AkDtqaYzGBTihpx6baSbOWYrQzGWUEEO5Vb2RWCyyXPYTwurcnRuLR+8hwdCzIYIFlWStezrKxyAKWVF6uuwM9JY51SLM8JJfvLfCwfH/JuTDMc0zuVkG65mWTN2E3aZs76gpvTptSuo/ZA033fRxg9/+9HKjctBAeFxRPVdTtMVk9Ads/HSZ/21CoBGQFMOW9PIc7Z3YBdljYyk7/lGPE+20xAE9UmdtMfloKDveCSuHNzkU5A7B4LiNxBhfVTh3e/Vc6N3WQTFlCwimq/hiAw8NdrxGDLxUZE0CxDrS2Ip9qWJT9zhiTPEVU6gMoqSWqxiTjRQq6YTN1guzmE7WE5aM9bbqUlNM+RTGAUTbNQXM+vW2+BDUdlrCskWWZZISgoYiWx0MMr44Jsu3DxxbL48DmNM+TMBunWiDR15CPWwnAMwyG6seksyZZPzZW5f4Du//grsi97liMdr24rpk5ocKMeRlQISELgJ3K9JQmqWsq958E0F1vNJV4glt83ENBiIPPM57NEDRcX+TR8QWcmiE45p4CNvsANqLjfoGrdWC3jP+B8KvH+KF/HXaN2TbsD8YCLcWRZafXs4LLJNzZIcE/ZXSBf7aKpj2NZkFGO2RoiZ86h/kl/4RZPDJ/oyGgUVTHIJwd405xIOaGWCu66bsfFGEYjR1ChzluQZJsojmcVVWc5sYWzCsX4CV+dy1PH4yJJs1Wiao14AtL33+b6+tzPRjJbSOdV3ANAJTn4EFk9FRQya7w7OTyMedfbbkk2yteauL9NbcbbfF9aFYSdp0ssKxzsBRcX+UTYMa5T2z4XoifVCQsoEIwPlKuxruJBE0HFiJ+Cw9L0NCc4gur5uEZiykC4CKbbda6yKfGafGODxJNVutlzwXhwbrbhyFkLw5IMikoC+5UlbhW1GRJPh1CTQu+otFOL5t4tYxXJ8zKkHd0LMaHihBsE1fh4WlHBwPeh1r82CZ2hr+WUC1KeX8ef/zWqFvOZz3DqPePVfMY4Evrz23d4wwcY9UFfvcDAWz/56dO7b7r+EFhXn9ag8QNchYAWRxhLy2f3uCjJp8md1ugjbjvpWJMF1ERA9fhP2D8t/tNUuaCeH1R5Y06goOsrvj/GPT0bA2Lcb90TCKPx1KCuHQwwxkBwf4xL66BIdvVJsyqCEYNm48VNTxCjKJsze3/btjRz1lScIBsIQa0rBaniknPJQT1xF2VqdkOybUinQVac/91H57/WAUf2sXtIb3qaWwnuN6OOhBbhWozFQDOIp75dra0Q0CKgKkvLZw+46MhnavHPndCmDEs9DtSGgEKSvakRUFAE1X84jRWZg+sC9IhPpjWuLdtJsP0OeS8h2cpIzmxhzm3C1jZ2Y9PFgGoDquY5dmOj8Wle0k4xo6gkic/xyZFxBx2NG9s7cPAkVH34iHJGsoxQnHRP9c3iWnBhkkDwZXiGk23uQ5b9gUQI6AeLA58HNIUkWqPBpSnRQ97sU905KnahZQqXeT67x0VFPrsmnoA2llATAUH1B9DGAgoy1KaJuOpthR9tr4fkvjQOLoCeHemxfXmXwRWGzmaX/ukevcdWSE9tYozAxmZZ4yx6b9MGXM1zjOk5y6HbcYF5cHLlJMGEgpTnQ4a9V0QPFHFmvFoFG1t7uyMe4wuQSr/nkmSNceKE0YjGxNyLnXQ8srvudtZPHP/xZXjSq64k+/Qjc7cpqU9CbqowArOJLajt8C7OBcnZ3WRy+1+C6mLFRUU+wO6Jp44wUExTwsGEECH4mCsWEFQUcI0EtGNffFu5kw6Lj/VokjA+knL2yYZsFYYnhO3LUvrHE1Y/3aFvjHOviZSKr6KvM64r4pI0Q+WBxEDWdU/3AFvbro3DQkBQEm782BvIaUqW/CxIkjjiObIOqytoJ0XGWSEb1/NQCfzAoZ6qYGv7xJZTXszddCCZGvE0la+KoTWLy9qF1Sx0goNlzGe3uHjIR2RxxBOjbVXkJhl2bNk05QBRs6JmuQ+sdYqs4keYQGrIewYzdmK6vA+269VT0iEdrNIZjZHxGKI6Z7NIo6g64CsGaC91iZmqzr0U3meo8H8YCAhmkonEsuvCDScza45Jkrhq20fWsKu+pE8Rf7sEB6SdfiNRPlV6zdVkDz08/zVqxWMniKeS31PzIhTS68W6yZYVDnaPi/LOhaTI+rKHBsslhtVymTglBHtqogJ/bGVisLiuV9F2VGZme+gUbeOxU6RtD2EwQLbHpJs53XNKuunrkOWgCT55VCau0YosfHka7SZoN0V7CbabYlc6sLaKrK64Wmph0N6tavAAQIw4Ul1bwRxZxxw5glnpF/GumfklaYp2EtTLtzXkEC1B9tF7S1FARQHnSCO9pqnecDNMr99MPEHtaWQyHSJsq1efX+DnEyoctFl2goi8SETuEpG7ReTVDftFRH7R779dRD5vp3NF5KUicoeIWBF5drT9BhHZEpHb/PKGaN/ni8hf+7Z+UfZxzpuLx/I5X5jmjquXvYFJCyiqAQdOhAA4KXZTUmu45GhUZuOHsjl5jhGhe/oo2hHEJiRDwaauOGhn0yWM4nNUNM/blzSpk6mIz9cQtJcito8MR27QrbuyDhNCVe3VVTiyjq73QRU5s4meesJPATDD/VYUAbW0cp9eaigkzpQWR563toDM6mrFo1DMDly07eM4Ux4Qige8+CEsnLsg2AU8v4tIArwO+ArgfuCDIvIuVf3b6LCvwtW5vAl4LvB64Lk7nLubCTtfj5sp4Fbg3cCLgP+25zfZgEuKfPYyWVVDY/HJ7n9MQFC44NwhpuqCg6oSLv4B1eJF6isTuMtap6byarnkzBY9INnuoEZc0miupFs5ydltN5PlOGtNPKEmHNtDF+OxIKOyioBkFsZedr2fCajnAaaTetfZOvayI4yP9lwFCGOQjU0/gdsUYvX3Sbb9fTXGkdBhcUPOg8p3c8rn3cY9XbjGbWG9pNdcXc7zU3yntGr5h3MLq6dGOk3uNqqkpGp3XdZnGlyh+oW0+QXA3ap6L4CI/BbwEiAmn5cAb1GnlrhVRI6LyDXADdPOnXfCTt/eUVX9M7/+FuDruBDkIyK/2KKNs6r6YwvqzwXDXGVTJk8OO6tWQ+U3G5FQXWZta1LsojaZJX/8iTLREypKMzsYkDz0COmpLqm1ThywuuLEAtsjGAywm4O5izja7SGS58hgq7C4RAyq1iVg5raYS2dfKyDsF4qYzQpyZJ388qMMr1phtG5IRooZ9UlXV5DBYOokbpqNsRubrmhrr1sm61p7eOJgbSDGzTKL+17sWJ4m2pfd+THSz7iptH7Aq948CYVK7V54I4UCtDbtdhTTaSQcKI9JaiKeUDOvTkQLgCJktrV44XIRiaeFeaOfEgbgWiAqS8/9OOsmRtMx17Y8twlNE3Ze68+vX6OCRfHCTpbPS4Af3+GYVwOHnnzaYhpJiZFJa6hORE1ihGn14fDH+8nSwhN4fWDTwQAdDh1hxDknndQRxGg8X7zLVwsAyqd4iSaUi6Y8PsyDrBhBuh10dYXx8R7DYwnjFaEzsNiu8dUGpg8soSyQtYqMRki36xSIoS7dYSPkWWg74+isfcZZ5EDpjiNxCrgwlQZ4WbQiKiBaEIU0TG0yFbHIoO5yo0pEe8UcFQ4eU9VnT9nX1Ej9RzvtmDbn1tE4YeccbS2EF3Yin3+nqr8x6wARObFDG5cE6gN8JYu6iOmEY6OSPFEyamUuIHD5EOIGtGluwaIsTMjuz8aEHJS9WCbOwgrdr0vKYxHF4RtkC0Vfr4eudhmvJ4zWBNsFk4mbyTXZWbUW7q/kOTIal4VPm+rTHVbolGK1LZHd+THSv/d0txKTQlTyxuWsRUKCYh4ri2gtbWEn4qm76uICowELm1JhYVLr+4Hro/XrgAdbHtNtcW4FMybsvN+fv1NbC+GFmY8AqvoLOzXQ5pjzgkUHfGOFW5tl7uYbgtR1d0PoB6WCb6JaQZOSTy02lMzZDTnU3ls8rULlCfgQEg/gc06cWs32UvK+kPeCQhD3/Nf26Ti44IbbWF9Be7+qgB9WZHfdXa7EirMQvxHjUgeCHDuqPF4Ugo3PM6a6xKgr2ppUbguL/bjyOm2WHfBB4CYReaqIdIFvAt5VO+ZdwDd71dvzgDOq+lDLc6u9FrnCCxWIJ+z07Z0Tked5lds3A++sn78oXmglOPCzi/4zXHCrOEdVv7XN+ecN037wu0gi3NW1Y79ylFFf7Ksp4goLqFYJASB79LH5+hr74vejEOhFPJiGh9diFlHYNWlrdkDu0z4Wg90NsrvudhZQrDiLLSAAEtfnJHJZx5MmTkNTiap6fk98zQXCLqCwqKpmIvI9wHtwPvY3qeodIvIqv/8NOOXZi4G7gQHwLbPOBdjlhJ3fCbwZWMEJDeYSG4jIj6vqT7c5tq3a7Z3AB4D3cRh1tRfoR9iGgNw+rbjfstOP7+6Chzjf5nzCqQdTGA4xWyM6m12yvvssuhtKspVBVsa7DtIg3gZF7pXKgYrLBQsovelpVVKJSSGo4cJXeZ5bHyeSxhWt94mAVGHcXnCwQ1v6bhzBxNveEL1W4Lvbnuu3zz1hp6p+CHjmPH2v4duBhZLPqqr+0O77c4mgbTWEiIDq8Z/83Ln5rwnldQ/ZQHlBoNbJz7e2MWcG9Hsdku0OopCeG5OeHrhE3sOYu+OVfKEwbPiOFbleB4CMso/dA0D69BubLaAQBwK/LRIkxKh/PpUYzw4EtABcqtNoi8jZabtwFlMrtCWf3xeRF3uGvSixU4B4tyV7QruFGq4hJ0itxQ4Gu2o/utDezr/EoHmODodw+gxpnpM+3nWD1XAEwyF2MChnST0sCMVOu12k13Puq07H5YgNR+47NhodCAICyO6+F/Ak1GQBqa2SEzSTyLSHhFkEtCAswu12CHEaeI6qfrq+Q0Q+1XB8I3bK8zkHhZzvR0RkCIz9uqrq0bm6fIAwrxqpQiKzD5xq/exFNbTEguGtH/JNdGu7HJyCuOKwigZCzbmjR9B+xyUFZ04R6eZnMhw0z3lBQk97aqSC0yoB1d1w01xpTbXc2swKvAssUO122PAW4CnABPkAb2vbyEzyUdUjc3bqwGPREtjmp8jJbYfQgXPRo8ihsm567H1JmD0fYpc60hRd7ZEf6SOjnOT05vm79h6Q3fNx0htvcCuhlFQ95lZsjwb9nUQJMQEtGJfiZHKzkkfnCc+0unMi8t/bbDtoWGhx0dBmmIDsgPjPL3WkN95AesNTSG94yq7OLz7LfSIeMVIUYF1UKf+pCDJ9r/bKewkacpZsKZk/yMju/USlokGBNlLpusQ6xoKrWQPQsqjopWAdichPznvOzE9ERPoichmuNMQJETnplxuAJ+2ql/uIRRNNrfED/dR4KaJ4St4t4s90n4inIKCw7HMVcBfLGiGDIelgjBmOYXvoXYwH3JXo51bK7v3k5HZo5zY7jwSkQKam1XIJ4GvnPWGnu/IdwIeAzwA+HC3vxFVS3TVE5P8WkTt9efB3iMjxaN8P+5Led4nIC/dynT1jSToXP/br842nfY6snplTNewRhZDizFmSx85hzm65yf+yzD2QHVQ5fq1f2cfva3ee6mxS2kcCCjGfpeUDML/yYqcKB/8eeDrwGlW9UVWf6pfPUdVf3m0vPf4IeKaqPgv4KPDDACLyWbgs3Ztx5bx/JWTjnncsSedgI3bPLLhi8cIRCVYkSTCdFNPtFstCSSHP0e0hbGy6quZZVhS2PUzIPn5fu891loptn+XyS/Ip8PnznrDjN15Vc1xm7UKhqu9V1aBlvZWyptBLgN9S1aGqfhyX0fsFi77+7M4trZ3DiN3GffYTQTlXxAeDFZSmbhrulRUnjd4Pd5wv4VTEnLzr78BaP1B1UVZiPjXX27xFRmMsyPpZ5GRyhxEi8loROSoiHeCPROQxEfnf257fNs/nvSLyDcDv+kzbReNbgf/sX1+LI6OAxrLeACLyStzER/RZnTygac6dnXCRkU56rQ/N7fCjzT7ZWp5/sDCteOSFhq8OPjHYG0HSFOn3XKWFJEG3LeKnZNhreZ5GconrpB3Uyf/iFIWo//l9D5BcN2d4eR+rGtRxieb5BHylqv4rX8bnfuClwB8D/1+bk9uSz/cDa0AuIlu0zPMRkfcBTfPl/qiqvtMf86NABrw1nNZwfPO0Km4+jDcCHJWT7pjdTnZ1kZEOuJkiK3WuoMyTqN2K9MnO8Mzuu59Dh30cXPaEMD2FhDJLbipoOqmzfApSyPZlEjrNfcF/6/KXDtTkf3NU5Mjvf3CSgJqmnj/f0Es2zyeg4/+/GPhNVT01z6zbrchnt/k+qvqCWftF5P8Avhp4fmRRtSkvPuWCO3yRpxHQQflBLhDplVdUcxtC4dIQoI3nVglQJb32SWQPtLvdFxrZPR93iYmw8ATChUItag0Tkcs0ddOch5k5w6SBe/0++krkgs9ry3M028dcpkUhuk87JnPH+T56YYhIgcweXBfmecB/FZE7gS3gu3wB6u22J7eeRltEvhZXDRXg/ar6+3N1c7K9FwE/BHyZqsa1Zd4FvE1E/i1Ozn0T8Bd7uVYFB/WHt2iEOVGiWVGBcr6UGaVJ0iddU9mXPfTwfvd2Icju+fiF7sJ0+IFf/EydjKLZZcNTfJJg8HGiOWefLS+jztIK06+DI56DlJfmVYDalESqFlecuYZCEu/v1QHIV7pUa7sFqOqrReTncbOW5iIywMXsW6HtlAr/BngOpWvs+0Tki1X11XP3uMQvAz1coArgVlV9lS8l/l9w85dnwHd70cMSc0Dz3M/YWKuiXbeGdoK1pFddGTWsldeVqR8uAA404dQQJgXU0aiUXqu6HJyQA9TrOeLY2nLlmHZRMDaQjCSJV7odkHJBoehpWE18X1sQkH3o05hrrppss621s0+uWb2EyQdAVZ+IXm8CmzMOr6Ct5fNi4HPVz2krIr8B/CVuqtRdQVWfPmPfzwA/s9u2l4D81BMkJ91kghJcA4Zma2gn1AgnRnrZZWSP73IKiEsN3h3mispmhFlqNRtjVlddTTZfEFSOH8WEyekGg/ZVGLxruVDYXWjCgTLRttdz67mtEGRwt4W+qlXnopyHeHcilqXg4MChtdsNOA6ECYeO7UNfllgw8lNPkBw/jpK7J+vCP90wp9AsNMy4GmNJQO3hnvS1ooSTtFNMg8BKH7veA2OQ7RVkY8vd/60t1Lassn0hCGcHopAkQfo9RAx2MEBUovmukuKeNJ7fJg52AUQnuhQc7AltyefngL8UkfnfDNAAACAASURBVD/GqdG+FJ8UusTBRn76NMmxY2ApSAiYJKJZqJGNNgwEyckT5KeemNi+RAOKJ3xXYbpS880IWLD9BEPXk9AQHY2QPK9YCPuGuihn1vUmZN0N/QsW3/awEL7I+hoyzlw8ChDxYtUZDvb8gYdJnnRVs9KtiXj2nYyE/BIWHIhIR1XHtW2Xq2orX3yrO6eqvwk8D/hdv3yhqv7WvJ29qODrUBXLAUZ+5oyv66VuUq6Q+KjWLcFFY7VxKY7zC1Duj14nx48Xrr4lWsAnM2ueO/fa1rab/sCAGYxREbTjLSKRdhbAPPBWV6iyULzupFVCnPVdVztRNLWpgKrmuZsjyVrMsaMuubbfcwm2aepcjaEIazz7b7hH/ruWP/DQ4t7/AqAqrZaLCSLy5SJyP/CgiLzX1/oMeG/bduZxuxngMX/OM0TkGar6v+Y4/+KA//GFH0hRwPSgzyQapuqGUmhQ5Pz49zKr79PURQ3bk5MnnNjBSOHWyx55dDe9vjSgFjvcRvIcA5itbVhdQTup+9WNx169uNjvlumkftK5xM33E7thh0PwrrEYYqRqnQQRQRBMGJmqrEvW15Hjx2BtxVlBY3EVt0VC4iBCjuZELrnoN+ZJKHvgIdJrrykbnmbh7LMr7hKez+e1wAu9OOyf4ERjr1DVW5mjxlvbKRV+HvhT4EeBH/TLD8zf50MO/0MznbR4Cqw8qfljDiyC/NraovQKWrN05mmnvkDRhoby/n57etWVLul1ianQzIsLzm3A2XPImXPImU10NF58Eqov8+Osj06R+EonrRRBbTy1ZuVIt4P0ephjRzAnjmOOHyM55hbT6xe/CQ3fuywvrx8Iq6a8rEx14U6uEGHWxgLab7fblJ/BtJ/GLIjIi3wh5btFZELIJQ6/6PffLiKft9O5IvJSEblDRKyIPDva/hUi8mER+Wv//x9F+97v27rNL5HUtUBXVe8AUNXfAb4O+A1f6aD1TW9r+Xwd8PdUddi24YsOscUjBjGlMiesH+S5UvKNDaeoivN8rJ3M+wnvISTxhf8xpj1RFu3W4kjRsemTrgGRQ5PIer6heQ7DYVF6B2PQsXeriwFZkAWktixA6i0fDa7ZGdOSVLYFq6eTIv2+s9b6neK7JFsjZGsbnjiNHQywW9uIMUiWOaKDCcFLcd2oOsQ0ZA88WJaPqjRy/n6Hi1C7+cLJrwO+Apdk/0EReZeq/m102Ffhch5vAp4LvB547g7n/g3w9cCv1i75GPA1qvqgiDwTeA/VEmYvV9UPzejyWESuVtWHAbwF9Hzg94GntX3fbcnnXlwphUuXfIjqZhkpymMVyXKHBKG6sYSKBzEBwSQJTXtvdQKKXXqegNRYJ2wIRBSuBYeqksL5RHjalzxHw+CeeCsBP/NqcIftkYQC0cUu48Z2w3oxrXV0fJ5jkgRWV7BH+mRH+2TrKckgJ9nukWx0MaqQZdjRCLs9xIBz96W+3lyWeWu86YEmunbD9uyBB8ukaDivxKOwqHjOFwB3q+q9ACLyW7hkzZh8XgK8xVeCuVVEjovINcAN085V1b/z26r9Vv3LaPUOoC8ivTmMi1cDVwFF9rmq3i8iXwZ8T8s2WpPPALjNz15adFBVv7fthQ4d4rlYoh+j5vnsuowHNOZjVlYmCaZuAcXb2qBS1NOU7rxGAvLnxNdaohlxOR4jzjXW6bjZSEdjl6Q6zqZLk9tepqi07T7Dwn1cDPS1PKGGifc0G7sk2U7K+PgKG9d32T4hdM8mdLY6rH8S0u1V5NwGjLPyunmOJMYRj1e9zcxNmvE+swedC65CQucFQm5bk8/lIhJbE2/0tSnBWR1xZd/7cdZNjKZjrm157ix8A/CXNeL5dRHJgbfjptOpMLqqvq/eiIh8nqp+hDnyM9uSz7v8ckmgKEHf7UKSOLdHmIxryo/gwFs/sStMFbV20gUHzVbQLNQtoBBTikQNSwLaHcSIU4OdPEF2+RFs15CeGWIeP4OeOYvd2t47AfmETjHiip1C8Z2QTLCeMGZdw25sklhLvpqwdZkwOg5bVwq9J4Rku8/6ue2SVNQLEkYj/5uyToG5gIe27MGHZscV92FOozksn8dU9dlT9rUppjztmNaFmCcuKnIz8PPAV0abX66qD4jIERz5vAJ4S4vm/gPweTseFaFtYdHfmLVfRN6uqt8wz4UPKiRJMCsryJF1WF2BNEE2t1A/MZcdHZCs8Tlgen33ohaLmSAgmLSCoL0lFNqNlVPGQK6oKKggYpYJqS0QHn5kdYX8snXOPGOVUzdDdqzD8duPc83/eBS57wF0a+9CBDGCrK4i62suFhOIYWsbAyUBzThfjWF4zAkRki0wKaQDSMaKbGyTD1z5xsKtGCT8TQ9t0yTlLfKPQh3CCgnVS0otiIScmGAhbrc2xZSnHdNtce4EROQ64B3AN6vqPWG7qj7g/58TkbfhXIJtyGfuG7Gox88bF9TOhUXIY+h2nA/76Aq233U/yMMwW2Yb1BNGww+xKX5THDTjQWrKPg1qulhZZ3VJPG0QXL5+Dh41gljQBMidZWFXOo35NLu+XrfjJNBGXCym23WWP0wlArO66pZjR9GVDjYBM4bVR5T+Y5BuKcl27sgshlrnriusoYalnlPUpCKdkWOXPfRw+d2LsWDrZ0GTyX0QuElEnioiXdxMznVP07uAb/aqt+cBZ1T1oZbnViAix4E/AH5YVf802p6KyOX+dQc348DftLwVP9XyuALz5PnMwgH3ObVHSHRzriiQ8AMJc6IcNqun260GiWdhmgUU9kHVCmogrIJ0iCwrvz8/d263b+PSRFCNWaWzYVm7PyFbTVh7QJFxvrjgunrX19kNt97tQDJ7wDSd1Flm3Q4cWUeNYe2hMUfus6gRzj61SzqEZDNzcZ3d9KmBXKRiuc9+/9mnH6kWxd0HLOIjUNVMRL4HpzpLgDd5Bdmr/P43AO/G1di8GxeD/5ZZ5wJ46fMvAVcAfyAit6nqC3GigKcD/1pE/rXvxlfiioK+xxNPArwP+LVZfReRZ+FED0ZEvt736XfbvO9Fkc9FBc1zZDzGbBuXk7A9PFgl6VvChCfXJkwTFswioLC/jibrya8HAloSzxyIpw+wFjMY0TudcuRTQtYXVh8ZIxvb2EU8xXvXl45GTgTQSYucIh2NJr7zRTWDNHXHrq6iK10kz+nf+xj2oU8z+uKb2T4pdDbgqCr23EbUQDSdQpsiqf4c928K8cyo/LCfBKQIdkHldVT13TiCibe9IXqtwHe3PddvfwfOtVbf/hrgNVO68vlt+ywibwKehVPMhQ9AcVVwdsSiyOci8EdFyC1sbTt1Tp6jw1FRmv4wYIJ0atM411GxUOrYSRwQEY9G7rX4WroPgd6LHXacYdIxMh4jgyGpMaz4KbbTM9u+AoFOVhzYBTQbY7dwZBJEAOBeR9Wni1I4ndRVsEhTdKUHuSJ55uYo+swbeeyZXbauUszYff52a8u1J4bkmJv8WLe2scPW847t8AZmf7/2l4AuaTxPVT9rtyfPTT4icgK4XlVvjzb/0G47cFAgaaf0s3tJq8s/cE+FukPQ9aBA0k5twxTiiayeqcQT0ERAMaFU3G9aywXRcvBZoj3UusF/e4ikKSa3mG0fO9keuoE+1Hubqf2nkpcz9XLZuDqBXS3eIt0uZnUVel0XEwq7zm26bYAeWePR5xzj7M1jOusjtodrDE/0WPVVz9MrLoMj64UbTkLcJ1xvB5FB48NfSzf4vhDQ4gQHhxV/JiKfVUuGbY22k8m9H/haf/xtwKMi8j9V9fsBVLV1MbkDhVAep9txwdUwOOcu+5scimKcIRB6gFEQ6H6gIp+eJJ441rPEYmDHGUYMbGwiozFsJV49WLqAFzb1dh2Ry0uSxFk8vS6s9F0JHnAu6cSgxjC+Yo3sSIfNa2Dl+DaXHdnk/qu6bFzXof/I1aQrffTkUWw3JdnYRgZb7jcXPdTNLRuf8z1njzxGeuXlc52zcx8W29whw2/gCOhhXP6nL9Gnz2pzclvL55iqnhWRbwd+XVV/QkRu3/GsgwwxroLv+hqc8NMTqbqnyu0hGkqcBKXbIcAE8YSnxontc1g9Maao4CrutrbihiV2hleFhUoC4tVvbp+6mUpVd55qoe0gHX9u0TnF96rbJb9s3W0b50juhg/bT8lXUobHEjSF8XbKo3ad5FzqJNXGoEfWsF0/3Ix8Yqp/TzMnvZv2fZqXbH0b2SOPkV69OAvoErd83oTLA/pryphPa7Qln9SXcvhGXHHRQ48wuRVH1p0bIVg43qWgeY7d2sJ0u4ci1lOJ8zS52mrigrlIp47IwlnGc/YXcbmdIL1GpJBZl4SU7yn246ybjlOvgU9g9cIDqy7n5/QZzFoPu94HEWw/QdOEbC1l85qUwRWCZLD24RXWHrYcvXeT5O4HHHEeO4KsdDCbQ/TUaXQwIJSqEhWXPzf1Jiz2O5YvsML6JW7s36equy4+0JZ8fhon5ftTVf2giNwIfGy3Fz0IKJ7msgwyXz9rPC5KfYQ4hQ05Cgf4ab4gnnofa8SzJ8IJaCKe2OqpYRnvWQB8kc2QnBksnuKBIsR9dut+i2rISd/Pojoq4zEhFpSfPk2S5ySXnYR+l3y9j+0YxusJ41XB9mD103Dk/ozVe09j7/o4WTZ2SdtyDNtN0U5CwuXIxsBXuHbJrIyzA+/WrsPNcnFwx4XzgDt9Iup/pVp2bXFqN1X9beC3o/V7cTWBDi2K4o2DARIIJrdOXFDPSzjAxAOOIIsqBlNcbBXiacrVmXM67R2J55C4Kg8VQqVnQFQiF9yUgqBzNa1IbrGbg9IV5tuOkZ87h8kyzBWXYVJDvpqCgslg9SFYfTTj/2/vzcNlqcp7/8+3qnvvs/cZmZRJBPE4gBpELqgxxiRGwUQwUX/RmGiMebzkys3gz5tIyODVeEOuN8nzmKiEJEZJIETjTSSGRIWfQzQQBQUEFEWcDiIIHM64h+6u9/fHWtVdXbu6u3rau3v3+jxP7V29qlbVqh7qW+9a73rfub2raKXuxqMAzc9jC/McOGWB2taIeGUr27+9QnXvEtp3EBoNotWKi6QwZgEa9ZjojFs+CzjRyYbnGa2rtaQn4EJ4P9rMnuInFp3v/cWnEms0SJZXmt0ZzXGL/FyeSU8S50lWll3w0JROYpIvLxs6p2xXm3+/gsUzJizB6omLeJ3JdTOS4zZAljSDfLYdNx9gd+8+JBFvnad6sMG2hqEE5h5yD8C143ZSqVaI73sAVStYJeLASREHT62zsKdCVJ9jey0hrjVQzQUXVWJDx6rrdn3NiBCjfJicYfExs9cOU79st9tf4BLI/bk/6W3e3Jpa8ck+RebdqNsm1k246GRJb/jR4mKzrKPFU+S51olOwtPB6hnZ/I1AZzLf3/T1UMdKXNdd85h+oL8Z8SOOm5+11epYrU5cb6B6QuVQjcphES2tEu09ROPRO3nkcVtY3FZh8fAS9tBedGCJxe8ZjfkKaoAM1EhcgC+plVguf12D0ivj8EjYfCmy15Oy4rNoZp/L5YWYjokv3ej0RZ8iwSki8UEc423bWoWdLJx+RChLh/6GIDzryCi/pwXHiqoVF3C0WnEPaN4iUkXNhxrVGkRJghqGlladG3bDWHygTlQzbHEBHn4EDh9m5z3LzO+bo74YMXcgcZNTDxzGDh123qX12micezLpUFLxUey7Fkedf2uGLZ9hKSs+D0o6Ff9W+7zdJfLYdkbSO4AXA6vA14HXmtkjkk4Gvgzc5Xe90cwuHOZcXZlyoelG46ALa9ImQlAcsaCTCHXKWAqtJ2H/o07yASRnjPjJT6Dx5a9udDNGg6JWpOu5Kjp4yAU49VELqFacu/fBJTf/J4Jk+xZqR2yBBCqH6s4jbscW4h3bXZiglTrVgzHxcsLCdw+i7+/FDi/R2H+gbV7RwI4Tay5BrYyp4CKtg0/KNwLCJNOhKCs+bwAuB54k6V7gG8DPDXnuj+OiqtYl/SFwMa1ICV83szOGPH7A0xSh7dvXbswmloOuIec7jfWMbWLrlFB54uNJ8xrFpz+Bxh3TL0BRGjQ0jpzQLGyBTIw2yafiWK0hCZurUt82x+qOClHNUD0iXnHWDQtbIElYftQC9YUINfz8sJUVdMQu4iQhOXzYhe9Z2AK1unP1zkZcGIK26N9mozVWZlh8JL2xoHgfcLOZ3dKrfllvt3uA50vaCkRmNnSUyFxUhBuBlw17zEB30uCe8c6drcLUZTefFK4Xua6LWbV6Ko9/nM+TRDMS+tQLkJ9/Y0vLzlpYrblYbAcPARBtXXTejEnixoJqNWRGdW9EtJpgFVHbVoUIGlurWGU7yXzMI6dWaMxD9SDADuKTtiODRvVYEFgkGnOw7TvLVL96L/XvP9SfBZQPxSNDZq0JummG3VEGCJ7tbrez/PLP/vVP4FI8XCjpg2b2v7tVLuvt9mjgfwHHm9l5kk4DnmVmfzV4u9v4ReDvM69PkfRFYD/w22b27x3a9Xrg9QBbWCzaJVBAY9++dgFKSQWoW8KtGfctLcR8Qsl+UpBPKuk4SeJu2nbIJ4Hz836iasVZQ+l1Li176yhBElEUkSxUqSw1SKoRVhVWEaob1UO4nERAY0tEYwvU52F1p9h/xgrVhTq172/h0TcscOSeRRQ93J/zQS4KduocoTh2Lt/j+Hxm++dwFHCmmR0EkPR7wD8AzwVuBrqKT1nXj/fhJpke719/Ffi1XpUkXSfp9oLlgsw+l+CcF670RfcBJ5nZ04E3AldJ2lF0fDO73MzOMrOzqsyXvJTAyLA0h88M/wKz41/4G9y0DiNmBurTwKbU60540i6wSqXl/RZFmCUkBw85kWokqF5H9QQ1EqKaW1RLiFYabNtTY24fVA8bFkFSgWRO1BfhnN3f5KSjHyaqi8rykJ5p2e+lD01E3XnoWb0+OsvHcN1uZZbNyUm4MfuUGvBYM1siM+m0E2XHfI42sw9IuhhIExj1/ATN7Pndtkt6DS5b3o/5fBWY2Qq+4WZ2s6SvA08AbirZ1kAJ2qyffiaaZrvnxhHQctpwuZSb6cKJI8QUOh/kbvSWmUOUrjfdlhsNJLlJ2UvLJEtLaGWFeK6KKjFRHEMsrG5Eq3Ws6sydqGEsPNTAYqGGYbGIakZSifjPWx9P5UDE0XfAtm8ewPbvL5c/q9N3sClALRfyprCOkBnvCLgKuFHSh/3rFwN/54dneka6Lis+hyQdRcvb7Zm4gaWBkXQuzsHgh83scKb8GOBhM2v4MD67gXuGOVd/DctNKg03WKDHxFLcAPVMjvs0rb7Ejye4/1PvhJERnVaRoXodWwZjBRpJ8zO3RgPzkaolESX+/agnLp5bJSLaUiGOXTdcZamBau47tfVeOPpWQ3Uj3reEHnyY+t595X53fSSkc9cz4oSQyZR/zkNgZm+TdC3wHFxE6wvNLDUSXtWrflnxeSMuL/ipkj6LS8s6rIPAnwHzwMf9/KHUpfq5wFsl1XHflAvN7OEhz1WOdGJaHOOyyKbE7T/GTSJGa8Z+ujzGKYraBciPC7hts/sDbH4Xkqi9E7sxZY/EflKpc5tP2ssz60mt7qNor51cbLU6WlmBJEErcUuMzVAcEW/bQn3BuVpXDtaIHzmMff8h7OAhlz6iWsEqFdc1Nq4UESNGU/YxjxJJr/Pj/jdnyi41szeXqd9TfCRFwBbgh4En4hTuLjMbyg/SzB7fofxDwIeGOfYwRNVKq187M/fF6u5Hl9TqwRrKknkvBrV+KsccTf37D466ZeuDZSyfBk1PsakkI0DN1wX7dHMCsOWVthw9aSpwqhWipe1Utrj5QdFSDfbuI9m3v3m+xIfYmZqHGWPWHQ5eJmnZzK4EkPRuKD/43rMD1MwS4I/MrG5md5jZ7cMKz0SiyM3onp8n2rYV7dyBjtiFdu1A27e5kDWVyvT8MNaDZtTs1tdoTQrvLsRHHkHlqKMAqDzqmJE2bd1oJO2LH6Sv3/PNjW7ZYKRJEwd4uDKf5M5WVlyA3pUVGvv20ThwgMbefWj/YSr7l6kcWEEHl0gOHGzFkEvPl/YwTHgwX0dJZ4MSDgeSzpV0l6S7Ja2xHOR4p99+m6Qze9WV9HJJd0hKJJ2VO97Ffv+7JL0wU/4MSV/y294pdR0I/mngFyS9UtIVwKqZva7nxXrKfsIfk/TSHg2ZfuK4NaFu26LLT79l3pXF0/BjGJD8eE6vUDttMeKKBagZZbuAeOdO4l27Mudzj48jT3O8HqQOB2lg2jQv1IyiyGc+jWMXksdHLFAkqNWI9h4ievgAHDhYmJrevIv31PQsWMmlC5Ji4F3AecBpwCv9dJYs5+HGv3fjppe8p0Td23EC8enc+U4DXgGcDpwLvNsfB3/c12fOdW5Be4+UdCQuqvUvAb+BmxbzVl9ein7GfLYCdUnL0EyXWugCPY2kPxhVq1CtYnMVTIJ63U1Uq9fd7OghQ9dPHKnQpP8HyfmTjv/knDWi+S1tOYWUTvJrpuQWZklrtjxOgOr3PzDEBa0z2XGyGXd90lwVLS6AhNVqqFFFDZ+mfsd299mvrDTdtzVXbVk640pTP4osqN0YzUd+NnC3n8yPpKuBC2j3GLsAuMJ7Bd8oaZdP8Hlyp7pm9mVflj/fBcDV3rP4G5LuBs6W9E1gh5nd4OtdAbwE+Ndc/ZtxV67M/5/wiwGPK3PRZSMcFMRl2YRI7qbos5kqSVywxKVlbHnFPaltJuEpIh+1uiwZB4ROXSaWJC0BAh8ZYK0xPU0CVL/ve+0Wmxn1EWbKnCYUx7Btq3vQ8AnizKwZBw4JVmsu+roPUtrzYW6Y8dVOXXeKRiMaxqi83U4AvpN5vQc4p8Q+J5SsW3S+GwuOVfPr+fI2zOyUHscvRdkIB2cWFO8DvmVm0x/dOiXtPlmtuXDvtTocOkRy6HD7RLtZICc8a1yti8Lx9Dse5idkWpSgrLeYROW4Y6nf973+jrdRZN6HWRUe8PN/0lBDZiRLy96JofW9SD3bygiKKlUXhdqPJZVGkRt79FbVmrojfIDsw9vtaEnZuYqXm9nl6WEK9s8fudM+ZermGepYks40sy90PUGJfcp2u70bOBP4kn/9VOBW4ChJF+bitE0n6UzoWg0dOuy62NKYVqurm9/iGYQy8eCKhglTke/WxZckUyNAsyw4WWy1hu3f79YzgUFdtmC/XvIBTpUq0Y5tqFJxjgl9JCdUpDYHhvZGjtqFu/SeD5rZWR227QEek3l9IvDdkvvMlahb9nx7/HqvY/21pOdRLFYpfwU8vVsjyorPN4HXmdkd0Byw+h/A23ApU6defCwxZ+kkS1i0DKnnziiDEE4iJeJddc9cmouKnS+H9oR2uW2F7Un3mSIBCjhhabpO5+cA1WstR4QSXW3R1kV0xC5scR7V61BWfLyDw5RNeP48sFvSKcC9OGeAn83tcw1wkR/TOQfYZ2b3Sfp+ibp5rsGFLftjXMi03cDn/MT+Az6IwH8Crwb+tKD+Tty4T7cbR88nsrLi86RUeADM7E5JTzezezaNA5yfv9Apd/1mpJnnp4P10iuqQRsdvgc9M6mCC2LZ7HrLiJnfv3L8cSBRv7fXA11goyl8WJP3gKtW3Dhfo4HVe3y3Gg1YXkErq677rnQDRpQJtSSjmGTqw5VdhIufGQPvNbM7JF3ot18GXAu8CLgbOAy8tltdAEk/hROPY4B/kXSLmb3QH/sDOIeGOvAGs+a79su4WJ4LOEeDvLMBZnby8FddXnzukvQe4Gr/+meAr0qaxw1SbQ5mQHDK0lF48pZSVrjKPoikwpO1jtLB4XTMICXtnjNrihAS9T33ljtXYENJx22oVFwqBtxE1ORw0rlXwRIaBw+iZRebsu+x1vX8HY8oaKiZXYsTmGzZZZl1w+VVK1XXl/8j8I8d6rwdeHtB+U3AU/pp+6CUFZ9fAP4bLpK1gM8Ab8IJz4+MpWWBsZJaPU0PtAxrhKeb4OT3S+kmRAXjPeZvGFLU8oJrRozOuYGbUTnh+GAJTSp+wF/VSut7EMcwPw+REKBGA1ZWunZrT7yDjzG9EcwngLKu1ks+dMJHzOyu3OaDRXUCk01WYLpaOUXrfZ5HGdEoFDGpKThNEUq74HJjQEBThCrHH0f9u0Nlcw+MGMUx0eKiy0g6Nwe1mgtP5cWGSqXdJV8bMG9uVK7WzHZst2EpNaNQ0vnALcC/+ddnSLpmnA0LbB6saO5QkeAl1rwxmSWtm1Q+gkCSuDEBMyrHHTvm1gf6IVpcREcdSeOU46idfAzs2gGLC0gRyb79JA89THL4sAtDtBHd3KMO2zOCCAfTig/583OSfte/PknS2WXrl/0kfg83C/cRAJ+f++Q+2xqYEKLFHllf05v8MHSbJ9RNgKApQpa6v6eiZLa2fj9OEYHxkXa1bdtKcuQ2Dp24yPIx89hcBbbMY40GyeqqW5ZXSGr19YsW4j3gxhIvbobFBzcF51nAK/3rA7hQP6UoO+ZTN7N9m8azLTCcuJTJXJrvMis6f0FXWtu+hV1xufq4oKRhrs0EELfSkMztq1Pdv0KyOEd8eAVbySS23ACLR3EuLcoojmkz3+12jpmdKemLAGa2V1LpyMJlxed2ST8LxJJ2A78C/Ef/bQ1MAsnhw0QLC31UGOAXlg2d40Uoa/3Ie7C1Qu10iC2XNZgil7ag6Rm3iWO9Th1pN2m9TrTvMFsOrTjrdb4CBw6SLPfMqjw2lIqiXFbTkTLDyeSAmg9IakCaCLT0k0XZn+9/x0VAXQH+DhfB9Nf6a2dg6siMwXQkG4I//0Sb1s92m/nFksQtfuymFRG60epOy4ztNKNF+3ApaZdc2i2XpmYIbBxWdxFB2H8A9h+Eg4fQ/Q+R7DuwodMYbIwTxlPrp9eySXknzpX7UZLejvOC/oOylct6ux0GLvFLYBboJjrdbiRF29rm8OQiIkjFFlFbW7IZwiQbswAAIABJREFUVNNjZU7pM2/GRx5B4+G9ndsWGCuWmMvjkyRNa8NWVjY+2sA4hW/zCktPzOxKSTcDP4abgvOSNJJ2GbqKj6R/psvba2bnlz1RYEoYVHS6kdZL5/BkifK7rj1H2zykIiFq9A4RFFgHsg4ikYs3nBTk7Nk0bG6rpieS/sbMfh74SkFZT3pZPv/H//9p4Fjgb/3rV+LivQWmlX7GcbqITtkB3LZgj63CgnZkrKPmObzDQdYqyjsrzHgunUnBfR8aYOr+3cjlfppaZvtrd3r2hR//eUbZyl3Fx8w+5Q/6NjN7bmbTP0v6dIdqgWkhm4OnEx1uDv16DWX3b6Yi79RFlz12xmtuTfcctIlQ48CBvtoUGAP+M7Wk83Cy4rgpPj1jvE06Myg+ki4GfgtYkLSfVoDRVeDyjhVzlPV2O0bS4zLZ8k7BBasLTDMDdLF1FJ0yT7DpDSd3DGUtnWwXXbaNqQj57rVUiJLDh3ufN7D+WNJKBJdJCKdKtT3sDlMQRqcLs9jtZmZ/APyBpD8ws4sHPU5Z8fl14JOS7vGvT8bl+Q5MKcnKsktzXUSBkHS1dMp2nWT3y0z4s8TaBSjdN7NP0kc+l8CEUJBTR3HsQu9A28ThngI0TEbTcTKD4pNiZhdLOgKXkmFLprxUr1hZb7d/8/N7nuSLvuLzfwc2G/3+wIdxQsgJELRbQclKH6H0AxuPIvf5pRZuPj2JIohbFq1lveB6iIvznou9a/6EiFBwOPgl4FdxSeduAZ4J3AD8aJn6Xef5ZNNnm9mKmd3ql5Wiffps+Nsk3SbpFkkfk3R8ZtvFku6WdJekFw5y/EBvkpXltUsaAiW3WL0gjXiPm4Al1lw67LDmGNlzBqaHNFGcKlU0VyXatpV421aXzjrjXGCrNWx5xcd28/O6SohJc57OpAhPymyH1/lV4L8A3zKzH8FlLi0daqSX5TOSdKkdeIeZ/Q6ApF8Bfhe40GdJfQXOk+J44DpJT8gkOwpsIFkBUiacSts+BWJTZNlkNm7+jLGbGR83TZEg9onj5qpg7vYSWeJcrtPPeSXBvKXT/K70EpVJE52UzSssZVg2s2VJSJo3s69IemLZyr3EZyTpUosws/2Zl1tpfYwXAFd76+obku7GBTW9YZDzBMbHIIKRekFlRSgIzyYhtXCqFRfFIo1OkaWZaXT6P3Mx291uwB5Ju4B/Aj4uaS9QOslWL1frk4drW3d8SIZXA/toJaU7Abgxs9seX1ZU//V4x4ct9IjUHJgMUlfc6b/3BFIyqasVR9BIsEaC1WpufYTBPCcKA02oQbYemNlP+dW3SPoEzlj5t7L1xxqaUdJ1km4vWC4AMLNLzOwxwJXARWm1gkMVfnvN7HIzO8vMzqoyP56LCAQCvbHEjQuu1rClpea4zrhiqk0MMzrmIymSdHv62sw+ZWbXmFnpwdqxio+ZPd/MnlKwfDi361XAS/36HuAxmW0n0ocpFwgENg5rNFyunnrN/c97vK0H48rdU8SIxEfSud7B6m5Jby7YLknv9Ntvyzp6daor6UhJH5f0Nf//CF/+Ku/olS6JpDP8tk/6Y6XbHlV42S7Pya2STir7VuXZsKD03nU75Xxa8YGuAV4had5PZt0NfG692xcIBAanmSiuKNr5+E++buccRVRrH5bmXcB5wGnAK73jVZbzcPfC3bihhveUqPtm4Hoz2w1c719jZlea2Rlmdgbw88A3fYLQlFel283sgS5NPw64Q9L1kq5Jl+5X26LUPB+5LHKvAh5nZm/1anesmQ0jCpd6z4gE+BZwIYCZ3SHpA8CdQB14Q/B0CwQmmCIrY1K900bNaLrUzgbuzkSQuRrneHVnZp8LgCvMzIAbJe2SdBxuwn+nuhcAz/P13w98EvjN3LlfiUuTMwj/c8B6QPkIB+/GicSPAm/FpUv9EM7HeyDM7KVdtr0dePugxw4EAmMgLzJ+orAiObf7SsU7GzScA8JmF6DRORycAHwn83oPcE6JfU7oUffRZnYfgJnd16EL7WdwIpXlryU1cPf43/eCt4Y09ueglO12O8fM3gAs+5PuBUqnSw0EAlNMOo8nnUgax22RDNL/iiI3vydXb1NTfsznaEk3ZZZseLIyTlad9intoJVH0jnAYTO7PVP8KjN7KvBDfimVHmEQylo+Q6VLDQQCU0omZI4iOesGsNUaWK0ZFseSCGp1iOQsoagCcQxpFtF0DGhMbXSNWv9bUh/zfB40s7M6bCvjZNVpn7kude+XdJy3eo4D8uM3ryDX5WZm9/r/ByRdhesSvKLbhQ1K2ceSonSp/2scDQoEAhOGt3qoVFAcEy0stCycNEpBGr2g0YA4RgsLaGELmp/vGAljUzAab7fPA7slnSJpDicK+YH7a4BXe6+3ZwL7fJdat7rXAK/x668Bml7GkiLg5cDVmbKKpKP9ehX4SSBrFY2UsoFFh0qXGggEphRLgNj/90n7fAidZjDQdJKpDEwonw5jXdq4AYxoDo+Z1SVdBHwUiIH3eser1AnrMuBa4EXA3cBh4LXd6vpDXwp8QNLrgG/jxCblucCe1FHBMw981AtPDFwH/EW+vZK+1OHK5ZpkTytz3b3SaB+ZefkAGRNN0pFm9nCZkwQCgenF6jU0N+dC5iQJrKxi9YL02F6ELDHXpRJ3zjCbWkNj7Y4bM6J73LF+MLNrcQKTLbsss27AG8rW9eUP4QyGojqfxEWhzpYdolwm0p8ssU9Pelk+N9Ma1DoJ2OvXd+GU9JRRNCIQCEw2lhiiAYqw+lLP/E7WaHS8MadOC+6FTXU201kMr2Nm3xrFcXrFdjsFQNJlwDVeYZF0HvD8UTQgEAhMPi5IqEgDgva0WNLuukg4zVIz/pslmZmXmUynLhr2HEguTM8kplDIswlD5/RC0gG6d7vtKHOcst5u/8XMLkxfmNm/SnpbybqBQGDaST3a/Hrv3Q3qteLU6ZlApGuy20qtLr5piAk3g+JjZttHcZyy4vOgpN8G/hb3dv8c8NAoGhAIFBHNtaaRhcRyE0I/VkjWCaHksazRQHW5iNhTYvXMeEoFAPzk1Wwa7W+XqVfW1fqVwDE4d+t/Ah7lywKBkdMUHj9JMZrf0r1CYDLJRD8oTCKYx5LJzVjaiRmNag0g6XxJXwO+AXwK+Cbwr2Xrl3W1fhiXMjUQGCttwtMsFNHCAsnS0sY0KjAwiuPWZNMynm1esCyJvOv2ZIvQLDocZHgbzmPuOjN7uqQfoQ+jpGxg0U9QoN9m9qNlTxQIdCPbzZYXnsDs4MZ54ukY72Hmu91qZvaQz+0TmdknJP1h2cplx3zelFnfgsu9U+DoHwj0jyrV4g1Z4ZGIFhdJDh9en0YFhscSklodJRnPthJ1Cp0RJpFN3KVWkkckbQM+DVwp6QH60IWy3W4354o+K2moiKaB2aaj4BQFolRLhKJFly49iNAEooio6mO6JeYcRSzpfx7PpItOlhkUH0nzZraCi4a9DPw6LuXOTlzWg1KU7XbLRjqIcLNgjy3d2sDM0y2+V3Mwuqi7Tcrt6/aJt20DoHHw4OgaGRgKxTFaXETzzlWag4dIlpanS0z6QMxst9sNwJnAZWaWRr1+f78HKdvtlo10UMd5N7yu35MFZoeuwSQzIlPoBVVQlopOnnjbtiBAk0S93jIGGptTdNqYTfGZk/Qa4NmSfjq/0cz+b5mDlBWfJ5vZcrZA0nzJuoEZYajoxV262zoJT0qwgiYDazRIlleAlVaU682M0RrPmi0uxHWz7QJenNtmwEjF5z9wZlaWGwrKAjNEc9ymn26VXsnFct1ta4RHBZaSD14ZrKANpsz4zgbm3hkHs9jtZmafAT4j6SYz+6tBj9MrqvWxuDStC5KeTiuI6w5gcdCTBqaLNjdocqFSoBVav0/WjPV0GOfpWJbdFgRo4lGliuJ4c6XZnkHxSTGzv5L0bOBkMlpiZqWSz/WyfF4I/AIuO94fZ8oPAL/VT0MD08WaqAKjuFEUjfV0EJ42i6eb8GT3CQI00Vij0frcczHd0myp6xLd4Fk/AP/xDyM51CxaPimS/gY4FbiFNOKsk+PhxcfM3g+8X9JLzexDwzQ0MD1ECwtupUN/tpuB3uevzotMGWun4xhPUXmSvYkFAZpo/LyfjuJSIvpB6X3XixkWH+As4DSfa6hvenW7/ZyZ/S1wsqQ35reb2R8XVAtsRnp1rRVtz43vDCQ8vayedP9kbduCAE0gRd+h7MTSQehHlJ71A0OcKEcILHo7bsrNfYNU7tXtttX/31awbbbf9k1K0+oBJw79WD8dnAna3KnLdLNlyosPmHbdFLQtY/0E1pkBx/5K0+PhphRlunDLHoqZj+12NHCnpM8BK2mhmZ1fpnKvbrc/96vXmdlns9sk/WCfDW3D5wO6AEhwKbp/wcy+K+lk4MvAXX7XG7O5hAIbSJ83lzVzeMo4FnQrTzFz+2SFJoqC9bORKPKu9vH6RaVOE9FlX3fj2We43QbQrM5tmOkHnbcMU7nsx/CnJcv64R1m9jQzOwP4CPC7mW1fN7Mz/BKEZ51os3pKoEhdl9zO/iQ9JpD2+2RaYv80JE9gfCgSqlbQXLVc+oQyx6xUO4dhSulD5CwasfDgut3KLD2PI50r6S5Jd0t6c8F2SXqn336bpDN71ZV0pKSPS/qa/3+ELz9Z0pKkW/xyWabOMyR9yR/rnVLnH5iZfapoKfve9RrzeRbwbOCY3JjPDmCIGYVgZvszL7cSuvE2lGhxsfgprkvXWyk6dY0M6tWWklo/nY4920+k6441GrC66jzWRjHxMvWAA6yXxV0ms+pzzhi+TWsOykjuWpJi4F3AjwN7gM9LusbM7szsdh6w2y/nAO8BzulR983A9WZ2qRelNwO/6Y/3df/gn+c9wOuBG4FrgXPJ5eiR9Bkze05BOu2+0mj3eg6Yw433VIDtmWU/8LIyJ+iGpLdL+g5utmzW8jlF0hclfUrSD3Wp/3pJN0m6qdbqcpwcBumT3gD6sgx8gre2pdu+bScaYB5PFLWWPH0ITLB+xo81Gli9NpouN+8Z19U7rh/Sbtp0GRFKyi09OBu428zuMbNV4GrckESWC4ArzHEjsEvScT3qXkAr5tr7gZd0vRZ3vB1mdoP3YLuiqI6ZPcf/325mOzLL9rLCA73HfD4FfErS+8zsW2UPmiLpOooDkF5iZh82s0uASyRdDFwE/B7Oc+IknyfiGcA/STo9Zyml7bscuBxgh46crEfd9MY77kHYIWnelLvdyHtZPwOKbNPqyd8MeoTTGZSQkmHKGNHvJnnu09u6vmyEKaJG5HBwAvCdzOs9OOum1z4n9Kj7aDO7D8DM7vPprlNOkfRFnCHx22b27/5YewrO0ZVB02iXDa9zWNI7gNNzJ+maTM7Mnl/y+FcB/wL8ng/VveLr3yzp68ATgJtKHmsySAdDJ1R4+rYEhu1+6+Xd1tyvT+EJ3WuBlALng+S5T3cvx5GT0Ojn+3e0pOw97HL/8AxQ1Lr8gTvtU6ZunsIH/H6PJel84I+A43FOY4/FOYud3uP8QHnxuRL4e+AncUHlXgN8v2TdQiTtNrOv+ZfnA1/x5ccAD5tZQ9LjcH2c9wxzrg1jAoUnWlgYvNshFZB+RajbAHTaljFZO2uaMsnWzyaLe7au5K1vdbCqu40VDnLa8j+FB83srA7b9gCPybw+EfhuyX3mutS9X9Jx3uo5DicQdHnA3+Prd2tHlqHSaJf9xR/lA8jVvEfDL/qTDsOlkm6XdBvwAuBXfflzgdsk3Qr8A3ChmT085Lla9Bqn2GRECwttCzC8tTCoN1Pe6hnmJpAkha7VZQjjP5uMDr/n5LkF4+kjFB6g5XTQa+nO54Hdkk6RNAe8Argmt881wKu919szgX2+S61b3WtwhgL+/4fBPeB7RwWyD/j+eAckPdN7ub06rdOBmpk9BDTTaAOlPTvKWj41//8+ST+BU8MTu+zfEzN7aYfyDwHjCeWT/ZKuV5dY/ocxwDnXxFnrunPJH1f2CXAQMSrbDVemPcNYPUMI0MRZQMHiGYz8fB8g+eGCe2CnTqoBGVUyOTOrS7oI+CjOi/i9ZnaHpAv99stwnmcvAu4GDgOv7VbXH/pS4AOSXgd8G3i5L38u8FZJdVxMtuwD/i8D7wMWcF5ubZ5uOcafRhv4fUk7gf8XN79nB/BrZU8ysYyzm6PoaazP8zSjSRf8uIor9PnLGrYLopcADdLdlm3PoBZapp51EaeJFKDAYGR+W8nzOmR6GfW4j9nIxhzN7FqcwGTLLsusG/CGsnV9+UPAjxWUd3zAN7ObgKeUbPYFwBLjTKNtZh/xq/uAHwGQNF3iU8YleJQiNMSx8ikMxtpNOIouuDIWUJnEcHkhHGS+Tp/7BwEKDMMsh9cxs0N+NcEFoI5x3X5Xlqk/zF1tTaDRqWcCxoPWCA9MZ3dMv1bYoKF2sgwopGEMaPOQt3pMav7PLqNiVBEOpglJOyRdLOnPJL3Aj0NdhHMM+3/KHqdst1thG4aou770Kygb4HVUKDopk+4gkbd+egnPODzchrTgggW0eRml2LQfmOGmH0wvfwPsxWWz/iXgf+C87i4ws1vKHmQY8ZmOd32YG3cu+dk4ctJ3FZ2UsmM+iQ3uiTYsJcZ3Ona59bo5DND91m2spxNBgKabrNXTUXAiXCfRqJiOu+CoeZyZPRVA0l8CD+LmDR3o5yC9YrvlY/c0N+G8ISabEVkMaYwpF7XXMYwQ9QyWOAsM627dIYr1sAQB2hzIbK0ARbn/IznP6I41RaTez/j5mN/oV3igd3id7YO0bBbIClFKkSAFockxri6QDINYPVmCAE0pI3alLnfOmVSfH5CUhjsTsOBf9xVYdJhut5mgn/DwTUHqw+KyxEYWgn7q6RZBYVSRqvPHWAcxDIyf5Ic7uFePE5tNbzczGyqjQcrmFZ8RdLlNlCiUHfeZNIrGe4q63KLc+jgGcvuMgq0oConopoV0zto6Wj9ukulMWj4jYQrvZhPONArELDDgTSLeVpRBPjANdBKGkXq/JSWXwBo2p+Wz2ayelGm1fnpR9F6Py/opSdZSSwUoWEGTif3Q0/1KD+snwT1ujzLETrB8BmZzis+UMdJxn410ty7C/zhTJwClXmrrFMV6lMTbW/43jQN9O/cExogslzahgwCN1OoZUSbTWWXzic9mtAyyTJP102a5OFG0RgMk98QoeREarVgqiob2eFt70PY2ZoWoiKkWJ0VE1crosoiOG/9dkoGRDZgLIuNy7b9qo8vtY2g2J5mOhCm5i60/I8lDv1Hn28gfRGKtxZL2pVluThzMnBXUWI/o4h3uOGPqNuklTuvGACGjFMdQqRBVK5P/oPNsH706+znmPtM1XWMjtX6s3BJYw+ayfEb8QxmoO2wIy2Tq3a5TwcmQiqoitbYlEUS0rKB053iYaBQlXLGL9hmjq3W8fftUWkDWaMBK0nogmtQkd8/6ATBDibAIt46cZVMQsd3FdRvh+WfU1XpUTPhjzcYysBAM8SMtZQGVOf56Wj9ZSwd3DemS0vY6bwU1GhtvAY26jmddLaBOaTz6/T5a4j6TQequJx2+4s2oA2atfTK9u6MVoGD5DMrmsnwmiXGPzUzC2E9WTMgJZ9FNS1HGElprBQ1lAfXzAx+VtTNpThOjEIpuSRaz37cNFiWd8zSnK/5P1voptni81TOOTKaBgZiwX89kMfQ4zIA/0PUebxqIzJjOGqum03Vntq2xgvIW0KS/B30IT7xz5xgbMkKKHmZyZZPQLayznuq62PJjO5mvXVvMtdH7tGTOY6WWwFqC5TNuNtJCGYfbda57rUmHsZ4szRuXf0+aY1yWtFtAZhDHkxv6ZgCLJ965k8a+fWNozAjJPzT4760iNT9PKxjXWy+iM08Hqc3YaAYQzaeFT9ezX6Gm9TOiBhnQCMIyKMHymWbW8yZQMK7TbEOmHfmxnixFFlKzLGsBJQk0GpPTVz5MO6TmEu/aNbo2rSM9u1PXgejppzXXZeYmjKafS6Z5eQeAti63ESPKWT3B8ikmWD4Tysg834axfkqO6fTbTdh2bW1WkLeAMMzMPaA25wJNEEWTZEtaafGuXTQeeWQMjRozGyQ68VOf1LJ20pt40XudDgA1M5fivvd5R4O0fFQEYRmYCftVb1I22mOo3/GT3DydTmM63ayctvk9hZtzxyTTpeOtIGs0WuNAk0yf3YPxrl2TbwVtsDNL/JQnEj/lia2Cfr4C3tIcq6NByoi83SSdK+kuSXdLenPBdkl6p99+m6Qze9WVdKSkj0v6mv9/hC//cUk3S/qS//+jmTqf9Me6xS+PGvo96kAQnwlmZG7XUE6AhhWdToLTRYgKBSh9nXVEyOdKGqMgdcy4OorzZ566J1KAiialrvPDU3z6E3xbWiLiLJhu2XLdP4uyFk6rfCwYIwksKikG3gWcB5wGvFLSabndzgN2++X1wHtK1H0zcL2Z7Qau96/BZR59sc9G+hpcWuwsrzKzM/zyQJm3YhA2XHwkvUmSSTo6U3axV/G7JL1wI9s3FQxzc0gFJzOvY2DRKdvWTs4JmfM3X+fHgSbRCur0ZBtp7ZIrj488gvjII0bbnlRABlkmnZwYZWO1WSzyVk9zW6bLbZTjPyMa8zkbuNvM7jGzVeBq4ILcPhcAV5jjRmCXpON61L0AeL9ffz/wEgAz+6KZfdeX3wFskTQ/+LswGBs65iPpMcCPA9/OlJ0GvAI4HTgeuE7SE8xs8LzVg7ZvAtxKR0o6/pOzNtxqf55rRfv1Rc4LMD8OZEnqaZXxhEuSVmDSMVHK6imiGatusO9MKkCNh/cOdv4SwpG+v9moE+n73tXK7jb3ZxyoXUSaVk+R8GS71CQsUsvqgfF2uWH9fBePlnRT5vXlZna5Xz8B+E5m2x7gnFz9on1O6FH30WZ2H4CZ3dehC+2lwBfNbCVT9teSGsCHgN83G8/A1kY7HPwJ8BvAhzNlFwBX+zfjG5Luxqn7DRvQvg2ntONBWZfuTPfVmhvOeohO0XF8u8s6IiA5kVgvV+x+zjOke3vWCiolRJ3m4XT5LjQn+CpCafLdTG7KonTw6ylAjdvvIn7ak9uFJoo6io5FeNFx5ZbdZ5wY/TgcPGhmZ3XYVtTS/IE77VOmbvFJpdOBPwRekCl+lZndK2k7Tnx+HriizPH6ZcPsbEnnA/ea2a25TZ0UvugYr5d0k6SbaqyM9MexaayeXHDPNaFvcuMxRaFx1hxrHG1srvZ2RGgLTDpJpO0Z0QTZ+Mgjuo8LZQRGkVxA0Ez3mSI1l171svulrwvPN86uuUz3Xy/hsViZJWoTnvxdLdvNNnKX6xGM+eDucY/JvD4R+G7JfbrVvd93zeH/N8dvJJ0I/CPwajP7elpuZvf6/weAq3AP/mNhrJaPpOuAYws2XQL8Fu2K26xWUFb4a/Zm6+UAO3RkOgsuc6QJ6sNej4CjvSyXbgP+JY85Ngq64aD1lN7WDdfA3aCyVtAIGNVxgDYLSF0+d+v1/kZr5wc19u1vritSm9hkNjTrg7dsvEVjiXmrJ3bb5fupzJCs2e3ZtIay43DjoOj9iYDIi1DcLjqtyaKZsrSpmcmlbV1uo3axTps+mgegzwO7JZ0C3IsbdvjZ3D7XABdJuhrXrbbPd6V9v0vda3AOBZf6/x8GkLQL+BfgYjP7bPNapAqwy8welFQFfhK4bhQXWMRYxcfMnl9ULumpwCnArXJfjhOBL0g6m3JPASUbsIHRBYoYowB1tFTK7lui3tgpeH/aoiB0ioqQTVQ3ajp1uZXyiHMCZJZ0FKC0vKMIFXxW8c4dbtPBQ+4Y+a625utc2yvu554tbXZfptZkBBAjckIVjUmAOrwvyS1fQc843X2+fizNmpZQpnstLu5e62rhjFKDRiA+ZlaXdBHwUSAG3mtmd0i60G+/DLgWeBFwN3AYeG23uv7QlwIfkPQ63Lj6y335RcDjgd+R9Du+7AXAIeCjXnhinPD8xdAX2IENGfMxsy8BzcEvSd8EzvKKew1wlaQ/xjkc7AY+N/jJ2scVylebzPQGXW8Ag4pNl7rrTsHn1SY45KwgPxZEpJZDwgD0VW+Ac3QToO7nUsduvGjbVuzwYfdeFQjOmnGxtN15UWo0aEu+N4zIFLlpD/oA2HQ0oKvwWHOcMNfuMUU2aGI2skjsZnYtTmCyZZdl1g14Q9m6vvwh4McKyn8f+P0OTXlG+VYPx0Y7HKzBK/4HgDuBOvCGkXi69dkdN9HC00MoSj+hTorgFNGxGy5vBflryDkkjIRRHCfT/dbLAiq0fnp8llpcxFZWmm1tE5y2dblzZ0XKDFutQRQhMo4Gkdz7Oaoo2QPXTbvbsl1tLeFpCZLf39QeUHQ9mLRxxyliIsTHzE7OvX478PbxnbD309hEWD4lfvx9dYVMstgU0a8VBM2uOFc4hs9vEMsqJ0DQfRyo/XydLZ8m8s4GafuKBKdSyYx/+P+rq2h+DqvX212Gmy75re62fr5nqcNC9mFpkC47y1g96fhNyxGBliD5uTtqmHP/Sk+bfvxjGu9xJwniMygTIT4bQocIvu27WG6XEXyBBxCAgfvap01sOlHSCnJlaVccvWfFe0p3uQ0zppRzwc5aOV2FqMRnr7m51ous4KRWQxRBtQJR3F7Rx6hL02arkWCrq06M/E3eatayiPr5PqVu3M05WwN0Xvznl+BZT/PH82M7fgzImg4Fmf+xmgJUxEgjWoNzg5r01B8TzOyKT54S3XLdRGBYYdqo3EFTQy9nBCjuiispQGvI1xljcNN+nA26oWrGukktoVR8ohjmq64sHadI/M+/EoNVoJG493OZZkijZkbT0heTsKaTfNjvZjNYaMZzLWPdpeM+yk6iTt+6fp0R+sI2/+9ujATxKaKEVbS2yjo9Ac3yl71DN5wr6u2QMPZJqb2OP+QE1J54b7ZshIDmdccRVkktwsjcycWOAAASAUlEQVS7CM9BveEyx1Zity53E7elJbfvIN83S7B6B0efbl3eRVMBml2o/nXWicBbPVkriNT5rcvbPNJxodDtNjBBfMrQ6wc4CnfuWRaVfulLhBqurIMIjcw9u4ywjXkM0ZaW0Y5tmbGdjOu15CY7ViKQuXumKiiOYHnFWUPzVWAOGg1kCaxGREBSqw/+/Syq1+exmmITFZRl94vkRDU/TpZ730dm+Rgj83abRYL4jIIgHBtDGREaZVfcFGAHD6FdOwuvT0mCZe7gLilb5iZdq3vrLELz8+59rVSIGw1sZWU4ERoHWc894XrBUusnK0A5izNYPpNBEJ/A9FNiPKiwKy4egaPBOEWsOX6RucF1nPRaMJ8nTyNBtQYmoSRxT+3psVN37WoFtsxDvQ5xjFZWnEUBRIqweq049tt60e3tXvfnCQviMwRBfAKbg0GsoAbFIXo20CpaM98nEjSyEccL5jHlu/P27oMjd7WXp0KzWkN+jk8zVXlibWnLbcdWksXt1BerzH3vANp3AGPZWUqp2/R6WUBlP4tOITYH3a/sOScxxceUEMQnsLkYJEKCD9EzlvA8Jcm6W68JuZOKRT/HyYxtmZzlQiNxFk0jad000+Nm/mv/IWrHbmPf4+aoHj/Hzrvmib/3MNY46LzoavVhLnU4suKRGMQqFuQClHbLjbQ9wfIZlAkKfBYIjJCCIKsdo2VnI2WPiy7OBr3ivpX11Gs7zsOPNKM+E+EiQ8eRE+UkgVrNLWmWWMl1p9XrUG84j7cIlo4R9Z1zsLAFzVW9EA74Pg2S1K7Dzb1w3GYjdGBEabRnkWD5BDYvvaygonxBsDYmWidG0D3XK9JBx7A7KUWi1jYQn86H8ZfciKGm5rhNdkKqqlUnTBLxoTqVw3NsvS8hWvVddlu3ouUV9372K0D9eoQW7Z8AcbsFIz+lJy1ritJ6TH0w29jxrykniE9g81OQtM69LHDJzo4DxXHR0TLHLRlHroPVUyg8+X0TcwIUJe7m28kSKIhorQceJjnhmOZrGk5Uhfd8W1nx6bwj151WrUClgs1VkRm7vnqY+rYqFolkxwJaqqGdO4hWaySrq+UEqCjdQwc6zZVTajzIeehZ2g3pFcgJkCEElhGg7ETTcc2xChEOBiaIT2B26JS6u8gKisAajfUdB+pixSjJCFB+W6dUCj66QTMeWoQbI5FQvdEad5+rOtHZUsW2zFHfNsfSMXMcOjZCCSw8lLBtpUGyfZ5qNUYHDxHPVUmWltuf/HMin801VIamI0iWBCwN1mBOYmgKkRespqt1xrU6qwl5J4Pgaj0RBPEJzBadUndnnBHWOCJAK1xNv5S1eno9lUfKCFCJG54XoOi7D1J/7KPdzTcGswhFEZEZqtfdONCWeZLFORpb56htq7K6M+bQsRGrO6F6AFZ2RWz7jljdWaFyMCaaq2JLDTRXRXU1x86ajhxpJtT0GssGyO0mVKmVWdD1Bt7asdb/tCzvYDBSpwOz4O02BEF8ArNJxgpqi2CeZvHMjwNB05Loi4LunlIRrfPdeX7mfkcBKkog58euKt9+gNqpx7bGRCKRMEe8sAUOHfZBRQ01jKhuRHWI6hAvt8K/1RcrbP3aI876qHtvt0bi6qbnW3HjQYpj55wAPjZcwd0+7TrzwlXULZec+eTWtStr8bR3vWXHfdL3RUUCNI6ut2D5DEwQn8DskhMgaM0JagoQNMeBFMfNSNCtY/QY9ylzs2tLf91h//RmmxUgOghZ/hgS1XvuZ+UJx3mHChfnwBbm0OoqHDiI6nXixgJqJCCoLVZJqmL+ESNehajm34sDh6Dhx4sAxVWoVrDlleYYmaoVtGUeq9Xb8gQpl8zOGg3XFZh0sZD8YE/LoaDV9SYAyzw8+O62jp5wIx/yCQ4HwxDEJzDbdBkHys8HMhrNAftBo1x37W4rEI1COmUoKMpamtk2//X7WX7isVg6UL9QJWosoMNL2CP7Ua1O7AOCzi9EJNWY1W0irsFCNROWKI7WdENGWxfdjdiLuKVzgXyiumabsnH12lJ3g8vcXHBZ6bhPruvNncirUmIoyZRlrnss83uAkFJhOIL4BAJ9OiI0+/rTG/wosqeWEZ7sjTuO2zOPdjpm1t0aSGJ3M06IUCNC81W0cwd2732+iogqEZWleaqHvFUosXxUheqBLcT7DqFKBRa2uMmmcYy2bHHzgOp1WFr2yems1YZmoruWg4T5uUWdrtWecmrLykmdCZRaPM7CSbc1x3ms1eXWFJvsZ5PpihsZkxTrbsoIk0wDAXA3kcyNpNuEVEvz3Aw72NxNNJr7RO1Lur2TK3i32G8Si1/+HlZxk0+TakyypUqydQvRo47GVldddtPVBtFKg3jFd7nVDYtFbcdc0zNO8/Noy3xr4mkaOSGN8hy59jWXaqW53pZIL4p6exRmRKVJ6jdirX2U6XJTrvttHPN/DJqTl3stvZB0rqS7JN0t6c0F2yXpnX77bZLO7FVX0pGSPi7pa/7/EZltF/v975L0wkz5MyR9yW97pzS+WFNBfAKBLJ0EKBvTLI2IkEYH6EQ3R4NO3W3dus7SsnReThR1cM/uFgg1YvGO+0nmRDIX0ZiPSbbOkRy1g+jE4901LS0TrSZUlhIqS0ZcA8yobY1ZfcyR2M7tToSqVdeOet2lZajVYK7qhGl+3u2TLpVKK99QNlhqRkjThSc/bs2cprYuNXNddkrX0y43yy1k6owDH+uu1NIFSTHwLuA84DTglZJOy+12HrDbL68H3lOi7puB681sN3C9f43f/grgdOBc4N3+OPjjvj5zrnMHem9KEMQnEMiTE6CsCBWG5EmS0dzgOglPtvusLVFc1F1oOhGJrXc8gFVEUnUilGypkmzfAscdA40G8eFV4pUG8UpCvGJENUBQX4xJdmzBti24rrdUhKoZMapWWkscu8VbR82uwmz7M4vtPqn9LbHMnB1v3TRJvAgl1hSjpig16/f/9vTDiCyfs4G7zeweM1sFrgYuyO1zAXCFOW4Edkk6rkfdC4D3+/X3Ay/JlF9tZitm9g3gbuBsf7wdZnaDmRlwRabOyNk0Yz4H2PvgdfYP39rodoyIo4EHN7oRG8jGX3/+frG+Xfvrc/1f77Ltu2M/ezE3r9tn/9hhD3CAvR+9LvnA0SV33yLppszry83scr9+AvCdzLY9wDm5+kX7nNCj7qPN7D4AM7tP0qMyx7qx4Fg1v54vHwubRnzM7JiNbsOokHSTmZ210e3YKML1z+71T9O1m9mouqSKxlXyjz+d9ilTt+z5BjnWwIRut0AgENhY9gCPybw+kbW2Z6d9utW933el4f8/UOJYJ/Zox8gI4hMIBAIby+eB3ZJOkTSHcwa4JrfPNcCrvdfbM4F9vkutW91rgNf49dcAH86Uv0LSvKRTcI4Fn/PHOyDpmd7L7dWZOiNn03S7bTIu773LpiZc/+wyc9duZnVJFwEfxU2ffa+Z3SHpQr/9MuBa4EU454DDwGu71fWHvhT4gKTXAd8GXu7r3CHpA8CdQB14g5mlbpu/DLwPWAD+1S9jQRZiEwUCgUBgnQndboFAIBBYd4L4BAKBQGDdCeKzAUh6r6QHJN2eKXuHpK/40Bn/KGlXZlthKIxppcP1v81f+y2SPibp+My2TXP9Rdee2fYmSSbp6EzZprl26PjZv0XSvf6zv0XSizLbNtX1B1oE8dkY3sfasBUfB55iZk8DvgpcDD1DYUwr72Pt9b/DzJ5mZmcAHwF+Fzbl9b+PgpAlkh4D/DhuYDgt22zXDh2uH/gTMzvDL9fCpr3+gCeIzwZgZp8GHs6VfczMfBx6bqTlb18YCmPdGjsGOlz//szLrbQmt22q6y+6ds+fAL9B+6S+TXXt0PX6i9h01x9oEcRnMvlFWi6OncJqbDokvV3Sd4BX4S0fZuD6JZ0P3Gtmt+Y2bfprz3CR73Z9byb68ixd/8wRxGfCkHQJzvf+yrSoYLdN6R9vZpeY2WNw136RL97U1y9pEbiElti2bS4o2zTXnuE9wKnAGcB9wB/58lm5/pkkiM8EIek1wE8Cr7LWBKwyoTc2G1cBL/Xrm/36TwVOAW6V9E3c9X1B0rFs/msHwMzuN7OGuUQ/f0Gra20mrn9WCeIzIUg6F/hN4HwzO5zZVBgKYyPaOE4k7c68PB/4il/f1NdvZl8ys0eZ2clmdjLuhnummX2PTX7tKWn8Mc9PAakn3Exc/6wSwutsAJL+DngecLSkPcDv4bzb5oGP++SBN5rZhT1CYUwlHa7/RZKeiEte8C0gDS2yqa6/6NrN7K+K9t1s1w4dP/vnSToD16X2TeC/wua8/kCLEF4nEAgEAutO6HYLBAKBwLoTxCcQCAQC604Qn0AgEAisO0F8AoFAILDuBPEJBAKBwLoTxCcQCAQC604QnxlA0sExH/8vfQRiJP3WAPVPLkox0GP/JUm3dNj+Fklv6rcd04ikXZL+W+b1qT4twVg/80BgWIL4BIbGzH7JzO70L/sWnwH5uk+/MDamJHz/LqApPmY29vclEBgFQXxmFEmPlXS9jyR8vaSTfPn7JL1T0n9IukfSy3x5JOndku6Q9BFJ12a2fVLSWZIuBRb8k/eVeYvGJ0t7i19/hqRbJd0AvCGzTyyXWO/zvm3/teT1XOITjl0HPDFTfqqkf5N0s6R/l/SkTPmN/jxvTS0FSc+T9AlJVwFf8mU/J+lz/rr+PBUlSS+QdIOkL0j6oKRtvvxSSXf69v+fLm0+RtKHfBs+L+kHffnZ/v3/ov//RF9+eqYdt/mQRJcCqbXzjjLvVSAwEZhZWDb5AhwsKPtn4DV+/ReBf/Lr7wM+iHswOQ2425e/DLjWlx8L7AVe5rd9Ejgrfy7gZOD2zOs3AW/x67cBP+zX35HuB7we+G2/Pg/cBJySa3v+uM/ACcUisAOX9+VNftv1wG6/fg7w//n1jwCv9OsXpu3GhX45lJ4TeLJ/r6r+9buBVwNHA58Gtvry38RFpj4SuItW9JBdXT6Xq4Dn+PWTgC/79R1Axa8/H/iQX/9TXNBZgDlgIf9edPvMwxKWSVpCbLfZ5VnAT/v1vwH+d2bbP5mLMHynpEf7sucAH/Tl35P0iUFPLGkn7qb8qcz5z/PrLwCellpVwE5cQMlvdDnkDwH/aD4gq6Rr/P9twLOBD/p4eeAEDdz1v8SvXwVkLZTPmUteBvBjOHH7vD/GAvAA8EycOH/Wl88BNwD7gWXgLyX9C07kOvF84LRM23ZI2u6v+f3esjGg6rffAFwi6UTg/5rZ1zJ1A4GpIohPICUb5G8ls67c/36o0961uyVzrE5BBQX8dzP7aJ/nKjpeBDxi/Y+BHMq15/1mdnF2B0kvBj5uZq/MV5Z0Nk60XoHLS/SjHc4TAc8ys6Vc/T8FPmFmPyXpZJxliZldJek/gZ8APirpl4B7+ry2QGAiCGM+s8t/4G6O4DKHfqbH/p8BXurHfh6N654qoiYpfVK/H3iUpKMkzeNyFWFmjwD7JD0nc/6UjwK/nB5D0hMkbe3Rtk8DPyVpwVsOL/bn2Q98Q9LL/bEk6Qd8nRtp5Qx6Rf6AGa4HXibpUf4YR0p6rK//g5Ie78sXfVu3ATvN7Frg13AJ0jrxMVpJ85CL7AzO8rnXr/9CZvvjgHvM7J24dANPAw4A27ucIxCYSIL4zAaLkvZkljcCvwK8VtJtwM8Dv9rjGB/C5Zq5Hfhz4D+BfQX7XQ7cJulKM6sBb/X7foRWjh6A1wLv8g4H2Sf/v8SF0P+Cd1b4c3pY6Gb2BeDvgVt8O/89s/lVwOsk3QrcAVzgy38NeKOkzwHHdbgWzHnx/TbwMf9efRw4zsy+jxOGv/PlNwJPwgnBR3zZp4Bf79L0XwHO8s4Dd+LTSOC6QP9A0meBrMfdzwC3y7mYPwm4wswewnX93R4cDgLTREipECiNpG1mdlDSUbikXj9oLunZerfjZOAjZvaUIY6xCCyZmUl6Bc754IJe9aYFSQfNbNtGtyMQ6EQY8wn0w0ck7cINrr9tI4TH0wB2SrplgPGclGcAfyY3Yv8IzuNv6pF0Ks76u3+j2xIIdCNYPoHAmJF0CfDyXPEHzeztG9GeQGASCOITCAQCgXUnOBwEAoFAYN0J4hMIBAKBdSeITyAQCATWnSA+gUAgEFh3/n97ZLh1lhx5sgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# First date with some rain\n",
"rain.sel(time='1995-01-01T15').plot()"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray (time: 48)>\n",
"array(['1995-01-01T00:00:00.000000000', '1995-01-01T00:30:00.000000000',\n",
" '1995-01-01T01:00:00.000000000', '1995-01-01T01:30:00.000000000',\n",
" '1995-01-01T02:00:00.000000000', '1995-01-01T02:30:00.000000000',\n",
" '1995-01-01T03:00:00.000000000', '1995-01-01T03:30:00.000000000',\n",
" '1995-01-01T04:00:00.000000000', '1995-01-01T04:30:00.000000000',\n",
" '1995-01-01T05:00:00.000000000', '1995-01-01T05:30:00.000000000',\n",
" '1995-01-01T06:00:00.000000000', '1995-01-01T06:30:00.000000000',\n",
" '1995-01-01T07:00:00.000000000', '1995-01-01T07:30:00.000000000',\n",
" '1995-01-01T08:00:00.000000000', '1995-01-01T08:30:00.000000000',\n",
" '1995-01-01T09:00:00.000000000', '1995-01-01T09:30:00.000000000',\n",
" '1995-01-01T10:00:00.000000000', '1995-01-01T10:30:00.000000000',\n",
" '1995-01-01T11:00:00.000000000', '1995-01-01T11:30:00.000000000',\n",
" '1995-01-01T12:00:00.000000000', '1995-01-01T12:30:00.000000000',\n",
" '1995-01-01T13:00:00.000000000', '1995-01-01T13:30:00.000000000',\n",
" '1995-01-01T14:00:00.000000000', '1995-01-01T14:30:00.000000000',\n",
" '1995-01-01T15:00:00.000000000', '1995-01-01T15:30:00.000000000',\n",
" '1995-01-01T16:00:00.000000000', '1995-01-01T16:30:00.000000000',\n",
" '1995-01-01T17:00:00.000000000', '1995-01-01T17:30:00.000000000',\n",
" '1995-01-01T18:00:00.000000000', '1995-01-01T18:30:00.000000000',\n",
" '1995-01-01T19:00:00.000000000', '1995-01-01T19:30:00.000000000',\n",
" '1995-01-01T20:00:00.000000000', '1995-01-01T20:30:00.000000000',\n",
" '1995-01-01T21:00:00.000000000', '1995-01-01T21:30:00.000000000',\n",
" '1995-01-01T22:00:00.000000000', '1995-01-01T22:30:00.000000000',\n",
" '1995-01-01T23:00:00.000000000', '1995-01-01T23:30:00.000000000'],\n",
" dtype='datetime64[ns]')\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1995-01-01 ... 1995-01-01T23:30:00"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Create a DataArray with the new times\n",
"ntime = rain.sizes['time']*6 # 3hourly to 30 min.\n",
"dates = xr.DataArray(pd.date_range(start=f'1/1/{year} 00:00:00',\n",
" periods=ntime,\n",
" freq=\"30min\"),dims=['time'])\n",
"dates"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'Rainf' (time: 48, lat: 681, lon: 841)>\n",
"array([[[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" ...,\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]],\n",
"\n",
" [[0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" ...,\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.],\n",
" [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1995-01-01 ... 1995-01-01T23:30:00\n",
" * lat (lat) float32 -44.0 -43.95 -43.9 ... -10.099998 -10.049999 -10.0\n",
" * lon (lon) float32 112.0 112.05 112.1 112.15 ... 153.9 153.95 154.0\n",
"Attributes:\n",
" long_name: Rainfall rate\n",
" standard_name: rainfall_flux\n",
" alma_name: Rainf\n",
" units: kg m-2 s-1\n",
" _fillvalue: -999.0"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Broadcast old times to new times and fill with 0.\n",
"rain_test = rain.broadcast_like(dates)\n",
"rain_test = rain_test.fillna(0)\n",
"rain_test"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x7ff62a20ba20>"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEXCAYAAAByAUkhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzsvXu8NFlZ3/t9VlVf9uW9zp2ZgWFgiDqIRkHQeIkBBfmoGA0eAwfP8RLEy9HEaMRLvH1QI+ecxKgIYoJIDmiiiBAlgiSSoOMoF8fR0RmYGWCYG3N5573s3Xt3d9V6zh9rrapV1dW9q/fu/b57v2//Pp/au+u2anV19/rV8zy/51miqiyxxBJLLLHE+YS50B1YYoklllji0sOSfJZYYoklljjvWJLPEkssscQS5x1L8lliiSWWWOK8Y0k+SyyxxBJLnHcsyWeJJZZYYonzjiX5HECIyJNFZENEkgvdlyWWWGKJ/cCSfA4AROQTIvKCsK6q96nquqrmF7JfsyAi3yMiHxKRoYi8uWH/t4vI3Z5E/1BEnhTtOy4ivyEij/jlJ2vnfkJEtvy5GyLy3h36coOI/LGIDETkzvheisg1IvIuEXlQRFREbmjx3l4mIp8UkU0R+T0RORnt+0YRucVf6/3R9i+J+rvpr7URLU+edm7Uhvpzwzn/YYd+/gsReVhEzojIm0SkF+07KSLv8O19UkRetkNbz/f3buDv5VOifSIiPy8ij/vltSIiO93HJZaYhSX5LLFbPAi8BnhTfYeIfBnws8BLgJPAx4HfjA75d8AqcAPwBcArRORbas18jSfgdVX9yh368pvAXwKXAT8K/I6IXOH3WeAPgW9o86ZE5GbgV4FXAFcBA+BXokNOAb8A/Jv4PFX9QOgvcLPffDx6D/dNO7eGz4nO+fYZ/Xwh8Grg+bj7eCPwU9EhrwNG/j28HHi9f29NbV0O/C7wr3Gf14eA/xwd8krg64DPAZ4FfDXwHTPewxJL7AxVXS4XcAH+E26A3AI2gH+FG0wUSP0x78cN9Lf4Y/4rbqB9K3AW+CBwQ9TmZwB/hBvs7gK+cR/7/xrgzbVt/w/wumj9Sf79PM2vPwY8J9r/I8AHovVPAC9oef1nAEPgSLTtA8Craselvg837NDezwJvi9afhhvEj9SO+3bg/VPaqHx+Dfsbz/XnPL3l+34b8LPR+vOBh/3rNd/nZ9S+Z/9mSluvBG6J1tf89/Ez/PotwCuj/d8G3Hohfi/L5eJZlpbPBYaqvgK4j/JJ/7VTDv0m3NP4tbgB8c+AX8c9qf4d8BMAIrKGI563AVcC/xT4lRlPvb8iIqenLLfv8m2JX+J1gGc2bAuv430AbxWRR0XkvSLyOTOudTNwr6qei7b9FaX1MS9u9ucDoKr34AfyXbY3L/6Xd6X9buwiFJEvFpHT0/rpX18lIpfh+pqr6kdr+4t74j/fL25qS1U3gXui45uutdv7u8QSwNLtdpjw66p6j6qeAf4bcI+qvk9VM+C3gb/vj/tq4BOq+uuqmqnqR4C3A/+kqVFV/S5VPT5ledYu+/pu4BtF5FkisgL8OO6pftXv/0Pg1SJyRESeDnxrtA+cm+gG4CnAHwPvEZHjU661DpypbTsDHNll3xfd3jz4Mtz7/gycW/P3RSQFUNU/UdX4HtT7GV4fadgX9hfvwX++fzKlrfrxTddaX8Z9ltgLluRzePDp6PVWw/q6f/0U4LmxBYMbzK8+P90EVf3vOEvs7cAncW60c8D9/pDv9X3+GPBOXMzm/uj8P1XVLVUdqOrPAaeBLwEQkTuigPyX4NyQR2tdOOqvNxM1kcAdfvOu29srVPV/qepIVU8D3wc8FfjMKYfX+xlen2vYF/ZPew87Hd90rQ1VXVYlXmLXWJLPwcAif8SfAv5nzYJZV9XvbDpYRN5QU2VtNAzIc0NVX6eqN6nqlTgSSoG/8ftOqerLVfVqVb0Z9z38i1nN4d10qnqzlgH5DwB3ADeKSGyZfI7fvlMfPxC1FdxId/jzARCRG4Ee8NGmNvYZxftuQKWf/vWnVfVxXF9TEbmptn/aPam/5zWca/eOpv07tLXEEq2wJJ+DgU/j1EqLwO8DzxCRV4hIxy/PEZHGJ2hVfVU0ANeXqX59EUlFpA8kQCIi/eAi8q+f6SW6TwbeCPx7VX3C73+aiFwmIomIfBUu4P0av+/JIvIPRKTr2/lB4HLgT6f0/6PAbcBP+OP/MU6R9faor30cgQD0/Po0vBX4Gm8VrQE/DfxuiCn5PvdxZGr8NTsz2ovv2dRzReRmEflcf8w68P8CD+DieU14C/BtIvJZInIC+DHgzf6ebOLUaz8tImsi8g9wysP/NKWtdwDPFJFv8P37ceB2Vb0zutb3i8i14iTz/zJca4kldo0LrXhYLgpuYLgP5176AZrVbt8eHV9RmAEvAO6O1v8e8AfAo8DjwP8APnfBff5J38d4+Um/7zhwO7AJPAz8HJBE534jLqYxwBHHC6N9N0fnPg78d+DZO/TlBn+PtnDqvhfU9tf7qTu09zL/eWzi3IIno33/Z0N7b27oz4Tabda5wD/yfd8EHgF+D7gpOvdLcK6uuL3vxz24nMWJT3rRvpO+jU3/Xl5WO3cD+JLad+hOfw/fT1U9KcBrcerJU/61XOjfzXI53IuoLt22SyyxxBJLnF8s3W5LLLHEEkucd1ww8hGRl3rlkhWRZ9f2/bC40ix3+UzuJZZYYoklLiKkF/DafwN8Pa6USQER+SxcQuXNuMz494nIM/QA1zlbYoklllhiPlwwy0dV/05V72rY9RLgt1R1qKofB+7G1f9aYokllljiIsGFtHym4Vrg1mj9fr9tJrrS0z5r+9apAwuRaiKIFH9Aov9he3htDCoCBtSIWxLcYsq2NGqO+mtAm7JQ2uS9Tzumhf6l96nNFhcoMbxuDYnbrV8jWpf6/oZjO480Xz+7cm12/zV62zP6UPTDIxlk5Ur4rIwgVhGrkFmwblFV8CIiKT53sOsrU/qh1fcfBEhBJKhU15v2hXaItqsWr7Wy/2DgHE88pqpX7HzkdLzwy9f08VPtHDIfvn34HlV90V6ud7FhX8lHRN5Hc2b9j6rqO6ed1rCt8ZsrIq/E5YjQZ5XnyvN31c8DDTHRS5nYV2wTA8aRixjjSCZJkCSBxECSuKXbgU6K7XfRfkreT8lWE8brCaMjwvCoYLslCalxrxFHNPFrjJbk46u5xcTl/lePUXHbALciWv1fnuaPqa4/7V/8RTt7XW3x8t5//oWIBawg1rcV/mvDunUL0Tq4/1f/u1umEuej3/hFDf1w7brrh9davU79lhR9UUTh2F+dcjtTg13tMjrWRTuG9NyY9NwQ88QmDAboYAu7PUSzMQCm23WfOcBI2f7yz/aE5UhGcgVPYJIpYq3fZpHMOsLI1RFS7tfDa5uD1WhbXuzTPIcsc/+tuv9qURveZPnZXCi8T3/nk3tt4/FTOX/xnie3Oja55mOX7/V6Fxv2lXxU9QU7HzWB+4Hro/XrcDkhTe2/EZfAyFE5ebAerfaKHUinsl2MIxhwxJMkjoQ6KRgDaQpp4v53UrSbVolnLWG0LozX3ZOyWEciYTAmrFN9rSqIZ4eCNxpy8mNecW1KdLwU/2MeCqTz9H8+q/BBAxoGtht/8M8A+MTPflEz8TT8J14nIp4ZuOJ1t/DodzUQkH+bIuX7Krb7+1xYmBr/FxR1ZGsMttdheKLHeD2hs7nzE7cdjRwBSY2t434E60hm/Hxiq0V18h7XrRrjHnywjsHFCGoNYjwBiTkQBLRXuK/S4X8fFwoH0e32LuBtIvJvcYKDm5hdeuXiQm2gmMvaCaSTJG57GsgngU4H7aTQSbD9lHwlJVtxxDNeF7LVQAIUT+XxgFgQSNhnqPrkgquHaHytWwhabq+PdYGDbvq/bp0cLNugxWB2w4848rjvp76oSjjeyolJL+7zk37hz1sPmlf8yi089p01App4r+KsiQoJ1I/x/fGuUzWGfDVleCLBjP0pnkjdAYpaO9E/tUplPlyR6S6wJnYsGtLJ85q2RW1Jkjg3oDVAjSwvAgJSlPFSB7VrXDDy8WVQfgm4AvgDEblNVV+oqneIyH8B/hbIgO++JJRuswbcJkvHvxZPOHjCkTR168HaCaSTuidn7SbkgXTWDOM1gw0z3dho/AnjQsRJAi5GRDTwSXAVudGysIqKk7w1hU61jlRqpBMGpR2IeLdunCf/hCOh+3/si6puvToZWbj2lz/iXVd5eb0dcPnrb+Hx7/iiqrUX+lyQCsV9iFstNnuyP3H7E/7zFfKVFJsKZhzHVYDEoNtD5/Kq9VHz3N23JKluD5+L+s8lkeLzqsSDpqHpXkxYQMYREDlCguaU1g9cFAS0tHx2jwtGPqr6DlxNqaZ9PwP8zPnt0QVEA/HERNNEOoWFE0gnSZyVEyygwsWWoKlzs9luQt5PyPuG0ZGE3Fc7q7ua1PpBUMqYBOLHCe+lUevHDuIH+PKpWiYGVpkkIH/KTd/7F80DUURCExZg/f7tYhC77jW3gBge+OHnTRDP9W+8wwXwE+Me2kX9wNnuWpf9qiOgWajcXyIXn793YuH0Z5/g+O0+5lMb3EPshtyiWYYdjZr7FmJ/8XXr3CF+owFsHHurudxCnCeGtTtbQOSI1txvcKgJSFHyAyakOEw4iG63SwvTiKe+PQweswQFwcXmrR7tddBOgu0m2F5C3kvIVgzZqrd2oEI6xQCMG0vERgOk4AcmSuvGW0qxtVNYQF6AsBCo9YNW9b5UBrHdwLd13Ws/WFgHBaGLuCf3og6VGzgn3Ecz+z1lve5uCwZcTEbBzQmEwH8ytIUYIQgSUIXxGElTJMnQrDqQS3gw6XZY/atPsfWs66eHd+LpeeLXe7nHcXtiCgK/WAjItpFnLtGIZXmdCwUxjW6leIAtBg4xxRNkMTh2Uvc6EhLQ60K3i3Y72H4Xu9IhX+2QrXUYr6eM1xOyVVMRE4iX2oYBLSiw4vXCJRe5o+oxk1KhBUE8EM4JSq+wPR78nv59Ppy3m8GnbhUuGlIjJKY8GEzBZW+8pUrsNFgcE9es/Qd3H/OcZDujd8aRX5BZy9gpyxqFAIB0u8hKH1lfg2Mz5sPzT/AaHmym7K+sz4r51PthQnzShA1Va3Y/Pr99hgI52mpZYhJLy+dCoIW1U6xPExTE8unElG62TuKtnRTbd9ZO3jfkXZfHU7mmRvGFYN3Eg6WtjoFF/DseIL3wQE1pITm5tXezER2r/oUfgedSsgXrJ5nc5dRUu/iBB99hEwLxW3UCAfCyYZjH+qkpyKccRIXMXd/K88VasGAGI3qPCaPjXSeNHucwztDxGM1zJ2muv40j63BkHXt0hexoH4qHDS1capPiD5lutKqd3xIKIgefV+biP3lkCR3ewXlp+eweS/I535gWRI9FBeGYIJ9uiu3UrB7tJGiSoD3vZuskZGuJIx0jZf4NVAkmDHzx9vAAnZTxn7LDFBs0xH8kctHFQgIjbtxJtJQTQ6XB1sQRixH8U3OhQCO44BbkvhFBOp3S2sky1IjLf1HrCKitq0iZHMinueMqfaB0vWXe2sktSWbpD8doNy2tnnGGjrPa+Yb08pPoZSfIj68wOtZleDJl+4Rw4q6Rt0SpWi7h9TQjZDckYYyXXMfbxLsxPfnFn+Uhcr8pMF7GfHaNJfmcL+wgKijWd0oWTaPcncSg3Q4U1o5b8p7B9gw2LdmglEdrhSFC7EC9OiDEHOJYT0FA9TiFif5LtARSo7SI1JTbUOGmf/kXpWLP54E0Djwz3DHh/pUklO9uAIutKhGk24GVFUfu4klne+gO3VIkb698O/kfb+HUt1WFB9IwYE3IvOPXWUksMh4jW/47YhVGYzQkdIZj0g7JZSfQK08yvGqd0dGE4XHD8KgwPgb3PaXDDe8eEpJYUY2SXb37TcRVSFj0k32wgpIEIfdxNDiMOUC6dKntCUvyOR+YJRkW00w6RdDby6eD9VMkjHrC8fJp2zXkvQTb8VZOIJga2UAcm6GQRAsKVgqhgPEJ6/iuFa+9i42gdgvkElxwQXwQl+hRVwqGyvY6qTTkgYR7NY2A/D2TBJepb8S5n+axgLQMgGPVEXsM499Yv+cG5sLKaE90Yp35U+T3zDpWG9ZH40lFmVeeaZ776gJlP8z6GnryONvXrHP6pi7jdchWIe8rtgs2Ve75ppSnv3VcEo/6Uj2qziqqayuKsju1//XXlc6X6sfwvVZrqRCQ9/lOJKHCwSchhXzJPbvGknz2G3sgnsLa6XQc4XihgabGVSnopo5wegk2EZer0dSFyOUVJ4vG7raShCjyd+rWTVDBFcQTBik/PgcLJ1g5hUoO54ILT9qNXiYjFNlcUlW2TZVZx8Ssio7H/nJzEhDeespzpxzrdvy5PsAUl5fZLQLBx+t1a6fpNMGRD6Dh/RQkEUmsI+SnT5MeXSdfOcHGk5Xs6hHdlTGaG+wwgdygufCxV3R4xptHjriDeCDItyudCIRja+uuDzMTV2soCKi68VCq4GIP9RLzY0k++4kZ8Z3GKgVNsZ00LcknTdBuiu2maMc4+XR3Up1UIZgZwe5K1QJwg5oASBG/KUrt5DUXml8PllIgnFCc1PrXMZGJH9+e/uoPTfRlwo0V36M4QdJqdUBKEldGCEDEZ9Qr80miI1+jKjrOkCx39z63LuYyHhcurnkHxKJu26z9Wh4Xiw2O/+GdUTKtFvGTws3WIDIAoNslWxF4yoDveOYtvGD9b/nU+ATveOzz+cjD17F5ZgUGCR97eY9nvGWrJB5vYUk+aTHPheAWhOa4TyCscJyUplYlDnigCUjIF5ZPcOnh8OkbDwt2EBYUr2vEI0ninuYD6XRcIVDtd7ArYXF12Wxa+/jiATxWMoXXUSmWeKCLpdNiQXLF5IrkjnRM5pexW5JR9Hrs9oX/klORWMfybTe4l/ejYtFMSygV4wbYYiDzFtE+SHML6ycPgfxxuYxG6GhUHexbDorH3+xqy+0os57ZOS1db6rOQoaJygXgCopqLyXvCMYoVg3HzYivXTvH665/L1/31NtZPbrlPvNM+NjLVpFcC+LBUlTJbt23aWjIHZJG12bkYj0kMmzFfS3bLEtMYmn57Ad2UrSFAaPI44lEBcG9FivZ4mTRblJYM85V4gTNEqyV4NYKr+sPZrE15A2dSiJpvFtKkiraNuW6GimsHjVgfWWfIMcOVlHRtnWtT7jRmvJTgpjAxzOcGMBbIjWJ9I5WQFtoWZGZ8djXXxPIcnQ0QodDH1OafzSpFGuduG6z1XPijz7mVo34OEzVggguLEmSquCg20XTFNsBMcq27fBo3qfDJn0Rvu3kn3HnxtV8+NE1zJbxyjctK1gHgsu1JBYb/d8pbiXGuQibrB9v8RTut0JhJ1CL+xx0F5wCo+Xz+66xJJ9Fom2ZnIa6bBMVqL2rTbsptp+iaYImwbcecmhkgkyK69aIp15ZOhBWcX4UBwrHuDZjC6q8lgbSScAmuJpj3tUmFsj9eBJccnEcqYYml1vYHsjHrK46It7YxI4zZyGatHxSh8bg+7xQqzDOQIbQcU//mufoaFz0ZddP4uFt1j6X+v5im0S5RsbfwIqWXZ0irUY+pE56n3eFxFg28h6fGl9Grk9w1Awxonz5yTv5kLkBLCQjccQTW1aBePLIApoQPLS8zzsVM433yeEiILtjEtcS07Akn31EI/EEMUFdWBAk1IF0wtQHvQRNTDVPJ8DHxCtEY7UqPKhbOrHs2T/ZS/AfEMUfCleZFu4zCFaOeEvHE6GAyRUrfr6c3LWv0etAbDf+xGS8Z+qAohYxgllZR66+Era2sWc3vIXi76+17j7meTGXTFOyZWuEtjNx8R3Ajmtxnl0QkHhxxkRtu6n90NKlZi2CjwGGeXM8MasXI1RgnPgk64OqsJn1uHd0BSNNuCzdoC9jbuw94ry+Y8GM4L6vOc5Tfu/xSeIp6rnZVlZPsHYK66eyb4r1c0gJSGEZ89kDljbjIhH9KKbGeCr5O8ZZPIWSrYzv6EoH7aeu3AlRvGYXiCcmq5TD8fEdybUS10lG6pahJd2yJONoQAoxpNBOrT2TgfGkIzkFGRWkNOOeAd6iKS3DwkV55iyECs1RkqkOh+WSjfdGPFGf7NipyKYW65xjADz3si/050QW5ER77l8coyvq9XU60O/B6gqs9KHXc5MCGimIUNKOb8AgaeqqVoOTMgOPjI7y0PgEj2ZHOZ2vcjpfQ3MpY3kjnPVTt3iCm7OttLqOqBCuWzc7Hxuai9YPYgxIEXJMq2UniMiLROQuEblbRF7dsF9E5Bf9/ttF5PN2OldEXioid4iIFZFnR9tfLiK3RYsVkc/1+97v2wr7rtzzjZqCpeWzSDRMfTAhpfYDivgk0SJ/p9ctJ3rrJG7wUPVGg/+bA0koWuPiLa4kfji2FvcJ26CI7QAFcdQHPJPj64X5nA/xeT8+h8a5t4rWiv8GcV41cceFuI/Jq1LUckbT8um14nIL9ysqogq42Tk3B+VhcYVmq2Czybb2ikU+Xau6gquGymcRq9xiN2fxOk2Lidnseg+70sUMM2Tkl2wFWVlx92ZrG83G7v6pIpll5XHlkUfWuK1/LdcfPc2p0RqP9da5vLPBqWyN9ExCuukEJJLhy/ho6VILFk88i2lAk8utEr+JrJ9C0RYUbqUFFMQHRfwnVr95yxcS1OqBtIAW4XYTkQR4HfAVuMk0Pygi71LVv40O+yrc3GY3Ac8FXg88d4dz/wb4euBX4+up6luBt/prfzbwTlW9LTrk5ara4KJYLJbksyg0zbkT/sc5PMZ40jFF3g5pinY7aDedmqtToEIyUWUCKOXVoUs+thMnnBYTlNWsIZN799rYNgsVvDsvTlBd/Z1b2Xzp81DjBteqRUUZogivm8aK2FqMKkoXr/McZVxze0l5f+skdj4GpTmv4Vyd/qZGcbcCtc+s2OYrLOhKl/HxFcZHEpLtlGRoSbYzZDvDGIOMx+jQV2Cwim5tkzyxwdrDPfoPdHlkfBmPX77G1SfPMljrstHv8ZFHr6dzVki3ndUjOTz2hVdw+QceLl1t8VTZRV919hQK9ffeJD6YfnAjAVUOOUAybEUY1Qsm7g5fANytqvcCiMhvAS/BzWkW8BLgLaqqwK0iclxErgFumHauqv6d3zbr2v8U+M1FvIl5sSSfRWAK8UjkcihiPIkppdRBYOBL5LiETh97AdQKmIg8YoSSAxbElNaPKwfjs/19Q3EB0SLWE36/1pNO5uTVwKRQwbfT+/3JQqBrv31rZf2Jb/1CJJdCbCDhzeTw5Nd+sHQVwUQ5FanLzjtpMRnaxLCl0VPw+RyA4oeKeQgo+gwah+DKAO//+fwl208ZH0nYPpGQDg3JtpJsJ6SDhNSA2Yx+xmqxW1uYJ07Te7DDZXec4Oy5hO1zqzywlTK6KmGQdXj4wRMcPQfpAMxYScaQbtlJ4mnq39Q3GcVuIrKpxH8arJ82bcclmA4KASlg20cuLheR2Jp4o6q+0b++FvhUtO9+nHUTo+mYa1ueOwv/G46sYvy6iOTA24HXeMJbOJbks1fU/c8NddoqxFOUyakWBQXKsiZQcbkVFk5D1qioOpIKA30Y3IKrLMQQhNLiCe2HeE2mZfmbwjVWvU73v32w1e048aY/c7XMggVUK9OioV/125h2qtMX1P3/aaeI51Tca/sx8DSVd6lVXQglYdqV16GQWQciL/bFbjctjwecFSyCpgm2YxivCXlXSLtK2nGiApMpplP9GWueY8+cJRHDejfFZOucG6Wc63Q4d7TPMy9/iDtGhmQb0i11cb2Rkm7l5JcfIXliAFu+nlwkva6+qYg04hJOUwiost5AOBPig2nWjxwsAppDcPCYqj57yr6mRuoD/rRj2pzbfFGR5wIDVf2baPPLVfUBETmCI59XAG9p0968WJLPXtA0BUL8upLHY0pVWzzTaGqKH29BJLG1YwETEZF34RQTFoTAsieoUEcMmcz/iX09Tg7t54QJY0gD8XTe0450Ypz8j7fw+D/7IpegCk6ODYTgeJPFIp3UzT0T4hwhcBTFd5pEHJqz88CzqFphdfFI2yoKEeEX6sJ4XyTggIh8QnDe112zXXcvMf4zzYWkl5B6OXblknmOPXcOc+YI3aN90ssTxMKRlW3+4fG7OHXzGp+442mk20pn05IOMpKNIebcNmwOsGfOgi9Ymlw1JeY8jYBitFG/7YRK/IcDYwGpCnmjDHVu3A9cH61fBzzY8phui3On4ZuoudxU9QH//5yIvA3nElySz4HBjHptzcSTVImn61Vtcab3tB9vHT54XRCQRsRjiuwdRzzWlZspXF9xn0NiIRR11+LrN4kD5sFlv3ZL7YLGj8ENFQKCCrDbcfcICpWXqKIDCjdbURkhyKtJgFqeS/1Jea+oZdyHQW9mNe5wbbXFxHrFNAmUr+uT9sXKRPXxP7FKOrCkmwnjNaADkokrIpu4wV2zScm1HY0wVr00Xsj7livXNnhS+gQ/cN0f8qq178aMLenmmPTsENnchsEWOtgqiEetkj34EADpNVc3vMeG72xc1BZqyaYNsZ8JqXXN+pl1by8wAdnFSK0/CNwkIk8FHsCRwstqx7wL+B4f03kucEZVHxKRR1ucOwERMcBLgS+NtqXAcVV9TEQ6wFcD79vzu5uCJfnMi53qtUUCA8JUzEEyW0yFkFR/tNEPrxLrCeuuNHRV+RYIyFJWKAhGuMFfX0vSafKtx2RjymB4QO+Pbmtnv8+LqYN1ZOnEcmOoDDCysoJ0O6669GDgtnXSosJBODbEi4ptCxicSqm3d7vtZP0UVhfeuvGDahRzq+yDQsK++ZQjrH38HDLOSTfHrD1i2LgmcVOgi1tE1ZX+aZCYS+pyxrK1hGwV9EjGNatnsQgf2b6B8ZqP95wbIRtbsDlAt4eumkOo5BDdr+yhhwFPQm0TTHeDtiGGC0xACq1k1Du2o5qJyPcA78E9Tb1JVe8QkVf5/W8A3g28GLgbGADfMutcABH5x8AvAVcAfyAit6nqC/1lvxS4PwgVPHrAezzxJDji+bU9v8EpWJJPW8xbnbouqe50yuoFUJKJ9e0E0e4vAAAgAElEQVREopmK+w1cUU+Rwv2GUhT7LMmmONlV+5eSTDSal0UTt11NOX1CcOPEWP2Tj6JJMldsYyZmna/WPWmPxm4StzRxxTwHA5dEmSTOPWgE6fWQtVXodqufQZxxHxIxYTF5P1B8/hOS3x2O3/q651RdoaYkmanJvCEep7B1/TqrnzhDmlnMdk73iZRsNcX2DCZT0nMj7OkzzV1IXLIyAu4jFB4aHOXnzryY+x8/Tncb0s0Mc2YTNjcd8QyHTjFXL+AaIXvo4aoVFBNRkwUDzhNgG1xvbVEXplB9oLgQBKQIY13MEKqq78YRTLztDdFrBb677bl++zuAd0w55/3A82rbNoHPn7Pru8aSfNpglutmBvEUkuqgbkuidqyW0xQ0uR0C6uo3tJi+ACKrR8rjiMnG4PplpEo8fntRlTq65PrfPubcgiMD+RyxjT1ArXuCl26n0ZUjnbSwJMmyMn7S7ZQuw9xNMSCqBelUBqV51GnBpUaDW7XpuCmuvsIyhcn6+4XLLZpTR8vtYbIYGedgM5LhGLNp6HghAplFNgdkgwGNMIJkOZ2zGb0zCauf6HDn6RvonBOOPgwn7tym++lzjuS3tp3FE5U0ao3YfTyHDLu8D4o2lfBpm7d1AQkoX5bX2TWW5DMNUwin0eIp1qu5PBVlG0wnmIBG9xsVAgLcdAbB2xc7xgzlAGdAk9gkciQX/1ZCfba4L8lWjvY7yKDjSTQD3SG2sQBoniN54qyf4JaMVVaerHU8LkvqhLhapwPdrnPBjSwT88XA/H0XP7PpRDM6c33iexMl+hbJpvG+wvrxBOTdcLEFtH39cfp3urhLXPiTPEeH1bl8JrC1TedUhxNjy/qDHfKekA4s3VNbmE8/gW4OnLXjp+FuaylmDz1M+qRrKu+lvAcNLt4mImmaIG/a8VOKz9bJ5nwSUKhwsMTucMHIR0ReCvwk8JnAF4SMWhG5Afg74C5/6K2q+qrz17HpX6bpeTyRpDpWtgWJNfgfQVL7gVFaPzWiAapuORMNxGFGzHpbuPYwxseGIMSeQrmV4PrRRDj75DQElTBj6GwpPatoN3WWRif1E5ntv+UDOBffaFSQSlEotKhaDYhB/cArxqCpOtdmljn3U1EPLRqA5um/uM+umCOo5n6aKvOe9rASqRMLZVu0r+Juq1lAhXvO4ioYxJPJhRp20RTbE8hzdGMTsZb0iXPux56mMBrBcIQdDIp2diVfr0+XEPJ3oNEC0rg+nN8/YfHU+7EDgexIQPsMuxi12yWJC2n5NJZ+8LhHVT/3vPamDelExzVNBlco2yQSF0BZ7LMINmv5w43db00EVLiU4h5FxBMRkpNXR0QTIzKCNq7tYDt4F5xvLoG8A3lHsN0U03GWjxYS8v0loMJKBBd7qLy3cgByhOQHHTHVQU6MIyofq2qtSis6URKP9Hrl9jDQ57aM9+jOcZ/Ri11aR1FFfIrVUycdFFdjjWDpgN3aKs4DClK04+nko3nuSGs0cuQ8zjD9XnFuk6hgHmQPPEh63bXNVk8ogur7OpV4au9n8k20yKOaRUD7aP0sSnBwqeKCkU/L0g8XFBPJkA15PMCkpDrx5BMTgypgXeQ3EFCAxZeu0UYCanxtKeu6eddZUUyySEoNx0thZQ2uKj/yeMprC0ju5bsdH6sqJufZX0jacYOiJ2sdjZzQQK3L/VnpuwOzrLCE3GCTI2opErB9HTDA1UPbHpZPz236kSRIt4NZW4W1NXcPVGE4QkIwfpzRmoijShPOwAyEQiWuU1g56snGRlZPkF3Xq1er3ZE4KhaNP64gsUUOyk0JphXpfvxdrxFP9B5mkfnUqdRr+xsJZ58ISJFlzGcPOKgxn6eKyF8CZ4EfU9UP7NuVWsZ2Ktsj4pmUVAeLxtYCsfjqBYGAKM+H4sdaEBAUardwvoQfdu5/SCbk9IgfpCLSKTrs1vOeIVtLiviCShAf+Mv7osh51wkTNDWVAp77hTDgy+qKI5/4Papzv8m6J4IsR0YjdGu7kBerNQh5daoKkWrQXK2vkEDzIBQsnm4Hs74GJ44xvO4YeT9BMqV7ekjy+AZy5iyab0LWQu1Wh0bRt2DZeAJy8R4tP+Mo3oMq8qe3TynJ02JAnTFdxUKQRA9Z4EQScQXryPoJKsSCeBrIsRFRYnIbEpo6FcNO15kTqixM7XYpYl/vnIi8D2jITONHVfWdU057CHiyqj4uIp8P/J6I3KyqZxvafyXwSoA+q3N2bhekE7bFhUKDAqtuIQT3QlMJea0REFRK68TCgkp7obpwaM9bTI3wIabOw2exR1dQ08N2BdMx5B2fZZ9BkpUEZDKwHdz03MGKMqZ54FsExCDdLmZlxU0ZYBLIMhdvyVxdN9IU+j10pevUXVvDwlrW7aEPkCdIKs5KCg8AeY5RLQLpRXXkOgHFxLO2CieOMbr6KBtP6jJec9fpP5FwBEgGzmpoLd8OcZzw3Qk30laJxxcJK4inYvVYFjpgLhrZJz9F+hSfYK/qyCjqv/t9lG5jjaZpqLgvpzwUuH+l9TqtPFP1tCkEFNpc2P2URSWZXpLYV/JR1Rfs4pwhMPSvPywi9wDPACZKfPvCfG8EOConZ4+RbcimdtwEEYXYDlSJp+iPRaZZNaruqbCQW08hIHdCtU+xiwKqku1wXuKvc+cnQASTJEi/B70eJjGkiaBJz2W6dwSGik2liP2E31AyCutVebEuOuQTBv2isKo3vay/x2laBv4T4+rfddMiPibqZ94ceeup30P6fad6AyRzs53qOCtiNpKHuFH5vvAJqtLvwcnjbF9/nMHVHYbHBU18xWcfj5lXQtz9ww8yetFzGpRbVeKRXCuEU5LOvtH+YhHK+7jAZam6jMoDFQQU8q/axM0Kd5ypEFAFTdtqMR8xNbeemImf2G7gvKbLmM9uceBsRhG5AjilqrmI3Iibv+LeHU6b0WBLIUF0bCPpuM75TbW8hhDfsaAmEJCU1k84VyLRAUxaQBU/OeWAF6oMhwoA9eMAvfdTlfM0z2E0dpySGEwipInBpkLSF0SlTFT1cBOLhSfx6tNh66D9HCjK5PgcHTdA1ay7cM/Cew7S9U7HJaTmtqwU3nVzImGkkLcLuG3jDB2NkKA6jK4vaYr0+9i1HtlaQt6lmBgvHSidgcVsjcuYU/s3SPc9H2b0wjJvL9R1i2M8sXKxLLXjP4db/2pP9/h8ILv3E6Q33uBWcvxDGaULTsTVdis8BHEMs53rsJGApv22I4tpaiLqgrAUHOweF1JqPa30w5cCPy0iIbL7KlU9tbuL7MG1BhOVlQPUWkdAYUCsxXcaLaAGibXb3mABuUbKYGwYmGFSMacKDz1SSpOh9KVnGToeY6yLBCXGkHYMtitkK8YVTbagxhUkTYaQDJVkmCOZnW+g3S2ClDorqygXEutwj8cZMkomC1F6l04Rm8oyn1CbFtaTO98F68UYNElc/k5wmwbrK01QY5Bc6QyUdNtZJOmWpXNmhAyc2KD14BV/r+JTaoq2IsYTjjss1k4d4bsYau6FB4egfPPWhxiD5raxFt1MRA8NE6j/TuNSTPtIQIosZDK5SxUXUu3WWPpBVd+OK+W9e8zpYmuc8noHv3KBOgFZd+4EAU0VIPgVW5g//l/N1Tbt2qfOFMo0UUXD1M+hpdEIu7VNkmUI0AkDtqaYzGBTihpx6baSbOWYrQzGWUEEO5Vb2RWCyyXPYTwurcnRuLR+8hwdCzIYIFlWStezrKxyAKWVF6uuwM9JY51SLM8JJfvLfCwfH/JuTDMc0zuVkG65mWTN2E3aZs76gpvTptSuo/ZA033fRxg9/+9HKjctBAeFxRPVdTtMVk9Ads/HSZ/21CoBGQFMOW9PIc7Z3YBdljYyk7/lGPE+20xAE9UmdtMfloKDveCSuHNzkU5A7B4LiNxBhfVTh3e/Vc6N3WQTFlCwimq/hiAw8NdrxGDLxUZE0CxDrS2Ip9qWJT9zhiTPEVU6gMoqSWqxiTjRQq6YTN1guzmE7WE5aM9bbqUlNM+RTGAUTbNQXM+vW2+BDUdlrCskWWZZISgoYiWx0MMr44Jsu3DxxbL48DmNM+TMBunWiDR15CPWwnAMwyG6seksyZZPzZW5f4Du//grsi97liMdr24rpk5ocKMeRlQISELgJ3K9JQmqWsq958E0F1vNJV4glt83ENBiIPPM57NEDRcX+TR8QWcmiE45p4CNvsANqLjfoGrdWC3jP+B8KvH+KF/HXaN2TbsD8YCLcWRZafXs4LLJNzZIcE/ZXSBf7aKpj2NZkFGO2RoiZ86h/kl/4RZPDJ/oyGgUVTHIJwd405xIOaGWCu66bsfFGEYjR1ChzluQZJsojmcVVWc5sYWzCsX4CV+dy1PH4yJJs1Wiao14AtL33+b6+tzPRjJbSOdV3ANAJTn4EFk9FRQya7w7OTyMedfbbkk2yteauL9NbcbbfF9aFYSdp0ssKxzsBRcX+UTYMa5T2z4XoifVCQsoEIwPlKuxruJBE0HFiJ+Cw9L0NCc4gur5uEZiykC4CKbbda6yKfGafGODxJNVutlzwXhwbrbhyFkLw5IMikoC+5UlbhW1GRJPh1CTQu+otFOL5t4tYxXJ8zKkHd0LMaHihBsE1fh4WlHBwPeh1r82CZ2hr+WUC1KeX8ef/zWqFvOZz3DqPePVfMY4Evrz23d4wwcY9UFfvcDAWz/56dO7b7r+EFhXn9ag8QNchYAWRxhLy2f3uCjJp8md1ugjbjvpWJMF1ERA9fhP2D8t/tNUuaCeH1R5Y06goOsrvj/GPT0bA2Lcb90TCKPx1KCuHQwwxkBwf4xL66BIdvVJsyqCEYNm48VNTxCjKJsze3/btjRz1lScIBsIQa0rBaniknPJQT1xF2VqdkOybUinQVac/91H57/WAUf2sXtIb3qaWwnuN6OOhBbhWozFQDOIp75dra0Q0CKgKkvLZw+46MhnavHPndCmDEs9DtSGgEKSvakRUFAE1X84jRWZg+sC9IhPpjWuLdtJsP0OeS8h2cpIzmxhzm3C1jZ2Y9PFgGoDquY5dmOj8Wle0k4xo6gkic/xyZFxBx2NG9s7cPAkVH34iHJGsoxQnHRP9c3iWnBhkkDwZXiGk23uQ5b9gUQI6AeLA58HNIUkWqPBpSnRQ97sU905KnahZQqXeT67x0VFPrsmnoA2llATAUH1B9DGAgoy1KaJuOpthR9tr4fkvjQOLoCeHemxfXmXwRWGzmaX/ukevcdWSE9tYozAxmZZ4yx6b9MGXM1zjOk5y6HbcYF5cHLlJMGEgpTnQ4a9V0QPFHFmvFoFG1t7uyMe4wuQSr/nkmSNceKE0YjGxNyLnXQ8srvudtZPHP/xZXjSq64k+/Qjc7cpqU9CbqowArOJLajt8C7OBcnZ3WRy+1+C6mLFRUU+wO6Jp44wUExTwsGEECH4mCsWEFQUcI0EtGNffFu5kw6Lj/VokjA+knL2yYZsFYYnhO3LUvrHE1Y/3aFvjHOviZSKr6KvM64r4pI0Q+WBxEDWdU/3AFvbro3DQkBQEm782BvIaUqW/CxIkjjiObIOqytoJ0XGWSEb1/NQCfzAoZ6qYGv7xJZTXszddCCZGvE0la+KoTWLy9qF1Sx0goNlzGe3uHjIR2RxxBOjbVXkJhl2bNk05QBRs6JmuQ+sdYqs4keYQGrIewYzdmK6vA+269VT0iEdrNIZjZHxGKI6Z7NIo6g64CsGaC91iZmqzr0U3meo8H8YCAhmkonEsuvCDScza45Jkrhq20fWsKu+pE8Rf7sEB6SdfiNRPlV6zdVkDz08/zVqxWMniKeS31PzIhTS68W6yZYVDnaPi/LOhaTI+rKHBsslhtVymTglBHtqogJ/bGVisLiuV9F2VGZme+gUbeOxU6RtD2EwQLbHpJs53XNKuunrkOWgCT55VCau0YosfHka7SZoN0V7CbabYlc6sLaKrK64Wmph0N6tavAAQIw4Ul1bwRxZxxw5glnpF/GumfklaYp2EtTLtzXkEC1B9tF7S1FARQHnSCO9pqnecDNMr99MPEHtaWQyHSJsq1efX+DnEyoctFl2goi8SETuEpG7ReTVDftFRH7R779dRD5vp3NF5KUicoeIWBF5drT9BhHZEpHb/PKGaN/ni8hf+7Z+UfZxzpuLx/I5X5jmjquXvYFJCyiqAQdOhAA4KXZTUmu45GhUZuOHsjl5jhGhe/oo2hHEJiRDwaauOGhn0yWM4nNUNM/blzSpk6mIz9cQtJcito8MR27QrbuyDhNCVe3VVTiyjq73QRU5s4meesJPATDD/VYUAbW0cp9eaigkzpQWR563toDM6mrFo1DMDly07eM4Ux4Qige8+CEsnLsg2AU8v4tIArwO+ArgfuCDIvIuVf3b6LCvwtW5vAl4LvB64Lk7nLubCTtfj5sp4Fbg3cCLgP+25zfZgEuKfPYyWVVDY/HJ7n9MQFC44NwhpuqCg6oSLv4B1eJF6isTuMtap6byarnkzBY9INnuoEZc0miupFs5ydltN5PlOGtNPKEmHNtDF+OxIKOyioBkFsZedr2fCajnAaaTetfZOvayI4yP9lwFCGOQjU0/gdsUYvX3Sbb9fTXGkdBhcUPOg8p3c8rn3cY9XbjGbWG9pNdcXc7zU3yntGr5h3MLq6dGOk3uNqqkpGp3XdZnGlyh+oW0+QXA3ap6L4CI/BbwEiAmn5cAb1GnlrhVRI6LyDXADdPOnXfCTt/eUVX9M7/+FuDruBDkIyK/2KKNs6r6YwvqzwXDXGVTJk8OO6tWQ+U3G5FQXWZta1LsojaZJX/8iTLREypKMzsYkDz0COmpLqm1ThywuuLEAtsjGAywm4O5izja7SGS58hgq7C4RAyq1iVg5raYS2dfKyDsF4qYzQpyZJ388qMMr1phtG5IRooZ9UlXV5DBYOokbpqNsRubrmhrr1sm61p7eOJgbSDGzTKL+17sWJ4m2pfd+THSz7iptH7Aq948CYVK7V54I4UCtDbtdhTTaSQcKI9JaiKeUDOvTkQLgCJktrV44XIRiaeFeaOfEgbgWiAqS8/9OOsmRtMx17Y8twlNE3Ze68+vX6OCRfHCTpbPS4Af3+GYVwOHnnzaYhpJiZFJa6hORE1ihGn14fDH+8nSwhN4fWDTwQAdDh1hxDknndQRxGg8X7zLVwsAyqd4iSaUi6Y8PsyDrBhBuh10dYXx8R7DYwnjFaEzsNiu8dUGpg8soSyQtYqMRki36xSIoS7dYSPkWWg74+isfcZZ5EDpjiNxCrgwlQZ4WbQiKiBaEIU0TG0yFbHIoO5yo0pEe8UcFQ4eU9VnT9nX1Ej9RzvtmDbn1tE4YeccbS2EF3Yin3+nqr8x6wARObFDG5cE6gN8JYu6iOmEY6OSPFEyamUuIHD5EOIGtGluwaIsTMjuz8aEHJS9WCbOwgrdr0vKYxHF4RtkC0Vfr4eudhmvJ4zWBNsFk4mbyTXZWbUW7q/kOTIal4VPm+rTHVbolGK1LZHd+THSv/d0txKTQlTyxuWsRUKCYh4ri2gtbWEn4qm76uICowELm1JhYVLr+4Hro/XrgAdbHtNtcW4FMybsvN+fv1NbC+GFmY8AqvoLOzXQ5pjzgkUHfGOFW5tl7uYbgtR1d0PoB6WCb6JaQZOSTy02lMzZDTnU3ls8rULlCfgQEg/gc06cWs32UvK+kPeCQhD3/Nf26Ti44IbbWF9Be7+qgB9WZHfdXa7EirMQvxHjUgeCHDuqPF4Ugo3PM6a6xKgr2ppUbguL/bjyOm2WHfBB4CYReaqIdIFvAt5VO+ZdwDd71dvzgDOq+lDLc6u9FrnCCxWIJ+z07Z0Tked5lds3A++sn78oXmglOPCzi/4zXHCrOEdVv7XN+ecN037wu0gi3NW1Y79ylFFf7Ksp4goLqFYJASB79LH5+hr74vejEOhFPJiGh9diFlHYNWlrdkDu0z4Wg90NsrvudhZQrDiLLSAAEtfnJHJZx5MmTkNTiap6fk98zQXCLqCwqKpmIvI9wHtwPvY3qeodIvIqv/8NOOXZi4G7gQHwLbPOBdjlhJ3fCbwZWMEJDeYSG4jIj6vqT7c5tq3a7Z3AB4D3cRh1tRfoR9iGgNw+rbjfstOP7+6Chzjf5nzCqQdTGA4xWyM6m12yvvssuhtKspVBVsa7DtIg3gZF7pXKgYrLBQsovelpVVKJSSGo4cJXeZ5bHyeSxhWt94mAVGHcXnCwQ1v6bhzBxNveEL1W4Lvbnuu3zz1hp6p+CHjmPH2v4duBhZLPqqr+0O77c4mgbTWEiIDq8Z/83Ln5rwnldQ/ZQHlBoNbJz7e2MWcG9Hsdku0OopCeG5OeHrhE3sOYu+OVfKEwbPiOFbleB4CMso/dA0D69BubLaAQBwK/LRIkxKh/PpUYzw4EtABcqtNoi8jZabtwFlMrtCWf3xeRF3uGvSixU4B4tyV7QruFGq4hJ0itxQ4Gu2o/utDezr/EoHmODodw+gxpnpM+3nWD1XAEwyF2MChnST0sCMVOu12k13Puq07H5YgNR+47NhodCAICyO6+F/Ak1GQBqa2SEzSTyLSHhFkEtCAswu12CHEaeI6qfrq+Q0Q+1XB8I3bK8zkHhZzvR0RkCIz9uqrq0bm6fIAwrxqpQiKzD5xq/exFNbTEguGtH/JNdGu7HJyCuOKwigZCzbmjR9B+xyUFZ04R6eZnMhw0z3lBQk97aqSC0yoB1d1w01xpTbXc2swKvAssUO122PAW4CnABPkAb2vbyEzyUdUjc3bqwGPREtjmp8jJbYfQgXPRo8ihsm567H1JmD0fYpc60hRd7ZEf6SOjnOT05vm79h6Q3fNx0htvcCuhlFQ95lZsjwb9nUQJMQEtGJfiZHKzkkfnCc+0unMi8t/bbDtoWGhx0dBmmIDsgPjPL3WkN95AesNTSG94yq7OLz7LfSIeMVIUYF1UKf+pCDJ9r/bKewkacpZsKZk/yMju/USlokGBNlLpusQ6xoKrWQPQsqjopWAdichPznvOzE9ERPoichmuNMQJETnplxuAJ+2ql/uIRRNNrfED/dR4KaJ4St4t4s90n4inIKCw7HMVcBfLGiGDIelgjBmOYXvoXYwH3JXo51bK7v3k5HZo5zY7jwSkQKam1XIJ4GvnPWGnu/IdwIeAzwA+HC3vxFVS3TVE5P8WkTt9efB3iMjxaN8P+5Led4nIC/dynT1jSToXP/br842nfY6snplTNewRhZDizFmSx85hzm65yf+yzD2QHVQ5fq1f2cfva3ee6mxS2kcCCjGfpeUDML/yYqcKB/8eeDrwGlW9UVWf6pfPUdVf3m0vPf4IeKaqPgv4KPDDACLyWbgs3Ztx5bx/JWTjnncsSedgI3bPLLhi8cIRCVYkSTCdFNPtFstCSSHP0e0hbGy6quZZVhS2PUzIPn5fu891loptn+XyS/Ip8PnznrDjN15Vc1xm7UKhqu9V1aBlvZWyptBLgN9S1aGqfhyX0fsFi77+7M4trZ3DiN3GffYTQTlXxAeDFZSmbhrulRUnjd4Pd5wv4VTEnLzr78BaP1B1UVZiPjXX27xFRmMsyPpZ5GRyhxEi8loROSoiHeCPROQxEfnf257fNs/nvSLyDcDv+kzbReNbgf/sX1+LI6OAxrLeACLyStzER/RZnTygac6dnXCRkU56rQ/N7fCjzT7ZWp5/sDCteOSFhq8OPjHYG0HSFOn3XKWFJEG3LeKnZNhreZ5GconrpB3Uyf/iFIWo//l9D5BcN2d4eR+rGtRxieb5BHylqv4rX8bnfuClwB8D/1+bk9uSz/cDa0AuIlu0zPMRkfcBTfPl/qiqvtMf86NABrw1nNZwfPO0Km4+jDcCHJWT7pjdTnZ1kZEOuJkiK3WuoMyTqN2K9MnO8Mzuu59Dh30cXPaEMD2FhDJLbipoOqmzfApSyPZlEjrNfcF/6/KXDtTkf3NU5Mjvf3CSgJqmnj/f0Es2zyeg4/+/GPhNVT01z6zbrchnt/k+qvqCWftF5P8Avhp4fmRRtSkvPuWCO3yRpxHQQflBLhDplVdUcxtC4dIQoI3nVglQJb32SWQPtLvdFxrZPR93iYmw8ATChUItag0Tkcs0ddOch5k5w6SBe/0++krkgs9ry3M028dcpkUhuk87JnPH+T56YYhIgcweXBfmecB/FZE7gS3gu3wB6u22J7eeRltEvhZXDRXg/ar6+3N1c7K9FwE/BHyZqsa1Zd4FvE1E/i1Ozn0T8Bd7uVYFB/WHt2iEOVGiWVGBcr6UGaVJ0iddU9mXPfTwfvd2Icju+fiF7sJ0+IFf/EydjKLZZcNTfJJg8HGiOWefLS+jztIK06+DI56DlJfmVYDalESqFlecuYZCEu/v1QHIV7pUa7sFqOqrReTncbOW5iIywMXsW6HtlAr/BngOpWvs+0Tki1X11XP3uMQvAz1coArgVlV9lS8l/l9w85dnwHd70cMSc0Dz3M/YWKuiXbeGdoK1pFddGTWsldeVqR8uAA404dQQJgXU0aiUXqu6HJyQA9TrOeLY2nLlmHZRMDaQjCSJV7odkHJBoehpWE18X1sQkH3o05hrrppss621s0+uWb2EyQdAVZ+IXm8CmzMOr6Ct5fNi4HPVz2krIr8B/CVuqtRdQVWfPmPfzwA/s9u2l4D81BMkJ91kghJcA4Zma2gn1AgnRnrZZWSP73IKiEsN3h3mispmhFlqNRtjVlddTTZfEFSOH8WEyekGg/ZVGLxruVDYXWjCgTLRttdz67mtEGRwt4W+qlXnopyHeHcilqXg4MChtdsNOA6ECYeO7UNfllgw8lNPkBw/jpK7J+vCP90wp9AsNMy4GmNJQO3hnvS1ooSTtFNMg8BKH7veA2OQ7RVkY8vd/60t1Lassn0hCGcHopAkQfo9RAx2MEBUovmukuKeNJ7fJg52AUQnuhQc7AltyefngL8UkfnfDNAAACAASURBVD/GqdG+FJ8UusTBRn76NMmxY2ApSAiYJKJZqJGNNgwEyckT5KeemNi+RAOKJ3xXYbpS880IWLD9BEPXk9AQHY2QPK9YCPuGuihn1vUmZN0N/QsW3/awEL7I+hoyzlw8ChDxYtUZDvb8gYdJnnRVs9KtiXj2nYyE/BIWHIhIR1XHtW2Xq2orX3yrO6eqvwk8D/hdv3yhqv7WvJ29qODrUBXLAUZ+5oyv66VuUq6Q+KjWLcFFY7VxKY7zC1Duj14nx48Xrr4lWsAnM2ueO/fa1rab/sCAGYxREbTjLSKRdhbAPPBWV6iyULzupFVCnPVdVztRNLWpgKrmuZsjyVrMsaMuubbfcwm2aepcjaEIazz7b7hH/ruWP/DQ4t7/AqAqrZaLCSLy5SJyP/CgiLzX1/oMeG/bduZxuxngMX/OM0TkGar6v+Y4/+KA//GFH0hRwPSgzyQapuqGUmhQ5Pz49zKr79PURQ3bk5MnnNjBSOHWyx55dDe9vjSgFjvcRvIcA5itbVhdQTup+9WNx169uNjvlumkftK5xM33E7thh0PwrrEYYqRqnQQRQRBMGJmqrEvW15Hjx2BtxVlBY3EVt0VC4iBCjuZELrnoN+ZJKHvgIdJrrykbnmbh7LMr7hKez+e1wAu9OOyf4ERjr1DVW5mjxlvbKRV+HvhT4EeBH/TLD8zf50MO/0MznbR4Cqw8qfljDiyC/NraovQKWrN05mmnvkDRhoby/n57etWVLul1ianQzIsLzm3A2XPImXPImU10NF58Eqov8+Osj06R+EonrRRBbTy1ZuVIt4P0ephjRzAnjmOOHyM55hbT6xe/CQ3fuywvrx8Iq6a8rEx14U6uEGHWxgLab7fblJ/BtJ/GLIjIi3wh5btFZELIJQ6/6PffLiKft9O5IvJSEblDRKyIPDva/hUi8mER+Wv//x9F+97v27rNL5HUtUBXVe8AUNXfAb4O+A1f6aD1TW9r+Xwd8PdUddi24YsOscUjBjGlMiesH+S5UvKNDaeoivN8rJ3M+wnvISTxhf8xpj1RFu3W4kjRsemTrgGRQ5PIer6heQ7DYVF6B2PQsXeriwFZkAWktixA6i0fDa7ZGdOSVLYFq6eTIv2+s9b6neK7JFsjZGsbnjiNHQywW9uIMUiWOaKDCcFLcd2oOsQ0ZA88WJaPqjRy/n6Hi1C7+cLJrwO+Apdk/0EReZeq/m102Ffhch5vAp4LvB547g7n/g3w9cCv1i75GPA1qvqgiDwTeA/VEmYvV9UPzejyWESuVtWHAbwF9Hzg94GntX3fbcnnXlwphUuXfIjqZhkpymMVyXKHBKG6sYSKBzEBwSQJTXtvdQKKXXqegNRYJ2wIRBSuBYeqksL5RHjalzxHw+CeeCsBP/NqcIftkYQC0cUu48Z2w3oxrXV0fJ5jkgRWV7BH+mRH+2TrKckgJ9nukWx0MaqQZdjRCLs9xIBz96W+3lyWeWu86YEmunbD9uyBB8ukaDivxKOwqHjOFwB3q+q9ACLyW7hkzZh8XgK8xVeCuVVEjovINcAN085V1b/z26r9Vv3LaPUOoC8ivTmMi1cDVwFF9rmq3i8iXwZ8T8s2WpPPALjNz15adFBVv7fthQ4d4rlYoh+j5vnsuowHNOZjVlYmCaZuAcXb2qBS1NOU7rxGAvLnxNdaohlxOR4jzjXW6bjZSEdjl6Q6zqZLk9tepqi07T7Dwn1cDPS1PKGGifc0G7sk2U7K+PgKG9d32T4hdM8mdLY6rH8S0u1V5NwGjLPyunmOJMYRj1e9zcxNmvE+swedC65CQucFQm5bk8/lIhJbE2/0tSnBWR1xZd/7cdZNjKZjrm157ix8A/CXNeL5dRHJgbfjptOpMLqqvq/eiIh8nqp+hDnyM9uSz7v8ckmgKEHf7UKSOLdHmIxryo/gwFs/sStMFbV20gUHzVbQLNQtoBBTikQNSwLaHcSIU4OdPEF2+RFs15CeGWIeP4OeOYvd2t47AfmETjHiip1C8Z2QTLCeMGZdw25sklhLvpqwdZkwOg5bVwq9J4Rku8/6ue2SVNQLEkYj/5uyToG5gIe27MGHZscV92FOozksn8dU9dlT9rUppjztmNaFmCcuKnIz8PPAV0abX66qD4jIERz5vAJ4S4vm/gPweTseFaFtYdHfmLVfRN6uqt8wz4UPKiRJMCsryJF1WF2BNEE2t1A/MZcdHZCs8Tlgen33ohaLmSAgmLSCoL0lFNqNlVPGQK6oKKggYpYJqS0QHn5kdYX8snXOPGOVUzdDdqzD8duPc83/eBS57wF0a+9CBDGCrK4i62suFhOIYWsbAyUBzThfjWF4zAkRki0wKaQDSMaKbGyTD1z5xsKtGCT8TQ9t0yTlLfKPQh3CCgnVS0otiIScmGAhbrc2xZSnHdNtce4EROQ64B3AN6vqPWG7qj7g/58TkbfhXIJtyGfuG7Gox88bF9TOhUXIY+h2nA/76Aq233U/yMMwW2Yb1BNGww+xKX5THDTjQWrKPg1qulhZZ3VJPG0QXL5+Dh41gljQBMidZWFXOo35NLu+XrfjJNBGXCym23WWP0wlArO66pZjR9GVDjYBM4bVR5T+Y5BuKcl27sgshlrnriusoYalnlPUpCKdkWOXPfRw+d2LsWDrZ0GTyX0QuElEnioiXdxMznVP07uAb/aqt+cBZ1T1oZbnViAix4E/AH5YVf802p6KyOX+dQc348DftLwVP9XyuALz5PnMwgH3ObVHSHRzriiQ8AMJc6IcNqun260GiWdhmgUU9kHVCmogrIJ0iCwrvz8/d263b+PSRFCNWaWzYVm7PyFbTVh7QJFxvrjgunrX19kNt97tQDJ7wDSd1Flm3Q4cWUeNYe2hMUfus6gRzj61SzqEZDNzcZ3d9KmBXKRiuc9+/9mnH6kWxd0HLOIjUNVMRL4HpzpLgDd5Bdmr/P43AO/G1di8GxeD/5ZZ5wJ46fMvAVcAfyAit6nqC3GigKcD/1pE/rXvxlfiioK+xxNPArwP+LVZfReRZ+FED0ZEvt736XfbvO9Fkc9FBc1zZDzGbBuXk7A9PFgl6VvChCfXJkwTFswioLC/jibrya8HAloSzxyIpw+wFjMY0TudcuRTQtYXVh8ZIxvb2EU8xXvXl45GTgTQSYucIh2NJr7zRTWDNHXHrq6iK10kz+nf+xj2oU8z+uKb2T4pdDbgqCr23EbUQDSdQpsiqf4c928K8cyo/LCfBKQIdkHldVT13TiCibe9IXqtwHe3PddvfwfOtVbf/hrgNVO68vlt+ywibwKehVPMhQ9AcVVwdsSiyOci8EdFyC1sbTt1Tp6jw1FRmv4wYIJ0atM411GxUOrYSRwQEY9G7rX4WroPgd6LHXacYdIxMh4jgyGpMaz4KbbTM9u+AoFOVhzYBTQbY7dwZBJEAOBeR9Wni1I4ndRVsEhTdKUHuSJ55uYo+swbeeyZXbauUszYff52a8u1J4bkmJv8WLe2scPW847t8AZmf7/2l4AuaTxPVT9rtyfPTT4icgK4XlVvjzb/0G47cFAgaaf0s3tJq8s/cE+FukPQ9aBA0k5twxTiiayeqcQT0ERAMaFU3G9aywXRcvBZoj3UusF/e4ikKSa3mG0fO9keuoE+1Hubqf2nkpcz9XLZuDqBXS3eIt0uZnUVel0XEwq7zm26bYAeWePR5xzj7M1jOusjtodrDE/0WPVVz9MrLoMj64UbTkLcJ1xvB5FB48NfSzf4vhDQ4gQHhxV/JiKfVUuGbY22k8m9H/haf/xtwKMi8j9V9fsBVLV1MbkDhVAep9txwdUwOOcu+5scimKcIRB6gFEQ6H6gIp+eJJ441rPEYmDHGUYMbGwiozFsJV49WLqAFzb1dh2Ry0uSxFk8vS6s9F0JHnAu6cSgxjC+Yo3sSIfNa2Dl+DaXHdnk/qu6bFzXof/I1aQrffTkUWw3JdnYRgZb7jcXPdTNLRuf8z1njzxGeuXlc52zcx8W29whw2/gCOhhXP6nL9Gnz2pzclvL55iqnhWRbwd+XVV/QkRu3/GsgwwxroLv+hqc8NMTqbqnyu0hGkqcBKXbIcAE8YSnxontc1g9Maao4CrutrbihiV2hleFhUoC4tVvbp+6mUpVd55qoe0gHX9u0TnF96rbJb9s3W0b50juhg/bT8lXUobHEjSF8XbKo3ad5FzqJNXGoEfWsF0/3Ix8Yqp/TzMnvZv2fZqXbH0b2SOPkV69OAvoErd83oTLA/pryphPa7Qln9SXcvhGXHHRQ48wuRVH1p0bIVg43qWgeY7d2sJ0u4ci1lOJ8zS52mrigrlIp47IwlnGc/YXcbmdIL1GpJBZl4SU7yn246ybjlOvgU9g9cIDqy7n5/QZzFoPu94HEWw/QdOEbC1l85qUwRWCZLD24RXWHrYcvXeT5O4HHHEeO4KsdDCbQ/TUaXQwIJSqEhWXPzf1Jiz2O5YvsML6JW7s36equy4+0JZ8fhon5ftTVf2giNwIfGy3Fz0IKJ7msgwyXz9rPC5KfYQ4hQ05Cgf4ab4gnnofa8SzJ8IJaCKe2OqpYRnvWQB8kc2QnBksnuKBIsR9dut+i2rISd/Pojoq4zEhFpSfPk2S5ySXnYR+l3y9j+0YxusJ41XB9mD103Dk/ozVe09j7/o4WTZ2SdtyDNtN0U5CwuXIxsBXuHbJrIyzA+/WrsPNcnFwx4XzgDt9Iup/pVp2bXFqN1X9beC3o/V7cTWBDi2K4o2DARIIJrdOXFDPSzjAxAOOIIsqBlNcbBXiacrVmXM67R2J55C4Kg8VQqVnQFQiF9yUgqBzNa1IbrGbg9IV5tuOkZ87h8kyzBWXYVJDvpqCgslg9SFYfTTj/2/vzcNlqcp7/8+3qnvvs/cZmZRJBPE4gBpELqgxxiRGwUQwUX/RmGiMebzkys3gz5tIyODVeEOuN8nzmKiEJEZJIETjTSSGRIWfQzQQBQUEFEWcDiIIHM64h+6u9/fHWtVdXbu6u3rau3v3+jxP7V29qlbVqh7qW+9a73rfub2raKXuxqMAzc9jC/McOGWB2taIeGUr27+9QnXvEtp3EBoNotWKi6QwZgEa9ZjojFs+CzjRyYbnGa2rtaQn4EJ4P9rMnuInFp3v/cWnEms0SJZXmt0ZzXGL/FyeSU8S50lWll3w0JROYpIvLxs6p2xXm3+/gsUzJizB6omLeJ3JdTOS4zZAljSDfLYdNx9gd+8+JBFvnad6sMG2hqEE5h5yD8C143ZSqVaI73sAVStYJeLASREHT62zsKdCVJ9jey0hrjVQzQUXVWJDx6rrdn3NiBCjfJicYfExs9cOU79st9tf4BLI/bk/6W3e3Jpa8ck+RebdqNsm1k246GRJb/jR4mKzrKPFU+S51olOwtPB6hnZ/I1AZzLf3/T1UMdKXNdd85h+oL8Z8SOOm5+11epYrU5cb6B6QuVQjcphES2tEu09ROPRO3nkcVtY3FZh8fAS9tBedGCJxe8ZjfkKaoAM1EhcgC+plVguf12D0ivj8EjYfCmy15Oy4rNoZp/L5YWYjokv3ej0RZ8iwSki8UEc423bWoWdLJx+RChLh/6GIDzryCi/pwXHiqoVF3C0WnEPaN4iUkXNhxrVGkRJghqGlladG3bDWHygTlQzbHEBHn4EDh9m5z3LzO+bo74YMXcgcZNTDxzGDh123qX12micezLpUFLxUey7Fkedf2uGLZ9hKSs+D0o6Ff9W+7zdJfLYdkbSO4AXA6vA14HXmtkjkk4Gvgzc5Xe90cwuHOZcXZlyoelG46ALa9ImQlAcsaCTCHXKWAqtJ2H/o07yASRnjPjJT6Dx5a9udDNGg6JWpOu5Kjp4yAU49VELqFacu/fBJTf/J4Jk+xZqR2yBBCqH6s4jbscW4h3bXZiglTrVgzHxcsLCdw+i7+/FDi/R2H+gbV7RwI4Tay5BrYyp4CKtg0/KNwLCJNOhKCs+bwAuB54k6V7gG8DPDXnuj+OiqtYl/SFwMa1ICV83szOGPH7A0xSh7dvXbswmloOuIec7jfWMbWLrlFB54uNJ8xrFpz+Bxh3TL0BRGjQ0jpzQLGyBTIw2yafiWK0hCZurUt82x+qOClHNUD0iXnHWDQtbIElYftQC9YUINfz8sJUVdMQu4iQhOXzYhe9Z2AK1unP1zkZcGIK26N9mozVWZlh8JL2xoHgfcLOZ3dKrfllvt3uA50vaCkRmNnSUyFxUhBuBlw17zEB30uCe8c6drcLUZTefFK4Xua6LWbV6Ko9/nM+TRDMS+tQLkJ9/Y0vLzlpYrblYbAcPARBtXXTejEnixoJqNWRGdW9EtJpgFVHbVoUIGlurWGU7yXzMI6dWaMxD9SDADuKTtiODRvVYEFgkGnOw7TvLVL96L/XvP9SfBZQPxSNDZq0JummG3VEGCJ7tbrez/PLP/vVP4FI8XCjpg2b2v7tVLuvt9mjgfwHHm9l5kk4DnmVmfzV4u9v4ReDvM69PkfRFYD/w22b27x3a9Xrg9QBbWCzaJVBAY9++dgFKSQWoW8KtGfctLcR8Qsl+UpBPKuk4SeJu2nbIJ4Hz836iasVZQ+l1Li176yhBElEUkSxUqSw1SKoRVhVWEaob1UO4nERAY0tEYwvU52F1p9h/xgrVhTq172/h0TcscOSeRRQ93J/zQS4KduocoTh2Lt/j+Hxm++dwFHCmmR0EkPR7wD8AzwVuBrqKT1nXj/fhJpke719/Ffi1XpUkXSfp9oLlgsw+l+CcF670RfcBJ5nZ04E3AldJ2lF0fDO73MzOMrOzqsyXvJTAyLA0h88M/wKz41/4G9y0DiNmBurTwKbU60540i6wSqXl/RZFmCUkBw85kWokqF5H9QQ1EqKaW1RLiFYabNtTY24fVA8bFkFSgWRO1BfhnN3f5KSjHyaqi8rykJ5p2e+lD01E3XnoWb0+OsvHcN1uZZbNyUm4MfuUGvBYM1siM+m0E2XHfI42sw9IuhhIExj1/ATN7Pndtkt6DS5b3o/5fBWY2Qq+4WZ2s6SvA08AbirZ1kAJ2qyffiaaZrvnxhHQctpwuZSb6cKJI8QUOh/kbvSWmUOUrjfdlhsNJLlJ2UvLJEtLaGWFeK6KKjFRHEMsrG5Eq3Ws6sydqGEsPNTAYqGGYbGIakZSifjPWx9P5UDE0XfAtm8ewPbvL5c/q9N3sClALRfyprCOkBnvCLgKuFHSh/3rFwN/54dneka6Lis+hyQdRcvb7Zm4gaWBkXQuzsHgh83scKb8GOBhM2v4MD67gXuGOVd/DctNKg03WKDHxFLcAPVMjvs0rb7Ejye4/1PvhJERnVaRoXodWwZjBRpJ8zO3RgPzkaolESX+/agnLp5bJSLaUiGOXTdcZamBau47tfVeOPpWQ3Uj3reEHnyY+t595X53fSSkc9cz4oSQyZR/zkNgZm+TdC3wHFxE6wvNLDUSXtWrflnxeSMuL/ipkj6LS8s6rIPAnwHzwMf9/KHUpfq5wFsl1XHflAvN7OEhz1WOdGJaHOOyyKbE7T/GTSJGa8Z+ujzGKYraBciPC7hts/sDbH4Xkqi9E7sxZY/EflKpc5tP2ssz60mt7qNor51cbLU6WlmBJEErcUuMzVAcEW/bQn3BuVpXDtaIHzmMff8h7OAhlz6iWsEqFdc1Nq4UESNGU/YxjxJJr/Pj/jdnyi41szeXqd9TfCRFwBbgh4En4hTuLjMbyg/SzB7fofxDwIeGOfYwRNVKq187M/fF6u5Hl9TqwRrKknkvBrV+KsccTf37D466ZeuDZSyfBk1PsakkI0DN1wX7dHMCsOWVthw9aSpwqhWipe1Utrj5QdFSDfbuI9m3v3m+xIfYmZqHGWPWHQ5eJmnZzK4EkPRuKD/43rMD1MwS4I/MrG5md5jZ7cMKz0SiyM3onp8n2rYV7dyBjtiFdu1A27e5kDWVyvT8MNaDZtTs1tdoTQrvLsRHHkHlqKMAqDzqmJE2bd1oJO2LH6Sv3/PNjW7ZYKRJEwd4uDKf5M5WVlyA3pUVGvv20ThwgMbefWj/YSr7l6kcWEEHl0gOHGzFkEvPl/YwTHgwX0dJZ4MSDgeSzpV0l6S7Ja2xHOR4p99+m6Qze9WV9HJJd0hKJJ2VO97Ffv+7JL0wU/4MSV/y294pdR0I/mngFyS9UtIVwKqZva7nxXrKfsIfk/TSHg2ZfuK4NaFu26LLT79l3pXF0/BjGJD8eE6vUDttMeKKBagZZbuAeOdO4l27Mudzj48jT3O8HqQOB2lg2jQv1IyiyGc+jWMXksdHLFAkqNWI9h4ievgAHDhYmJrevIv31PQsWMmlC5Ji4F3AecBpwCv9dJYs5+HGv3fjppe8p0Td23EC8enc+U4DXgGcDpwLvNsfB3/c12fOdW5Be4+UdCQuqvUvAb+BmxbzVl9ein7GfLYCdUnL0EyXWugCPY2kPxhVq1CtYnMVTIJ63U1Uq9fd7OghQ9dPHKnQpP8HyfmTjv/knDWi+S1tOYWUTvJrpuQWZklrtjxOgOr3PzDEBa0z2XGyGXd90lwVLS6AhNVqqFFFDZ+mfsd299mvrDTdtzVXbVk640pTP4osqN0YzUd+NnC3n8yPpKuBC2j3GLsAuMJ7Bd8oaZdP8Hlyp7pm9mVflj/fBcDV3rP4G5LuBs6W9E1gh5nd4OtdAbwE+Ndc/ZtxV67M/5/wiwGPK3PRZSMcFMRl2YRI7qbos5kqSVywxKVlbHnFPaltJuEpIh+1uiwZB4ROXSaWJC0BAh8ZYK0xPU0CVL/ve+0Wmxn1EWbKnCYUx7Btq3vQ8AnizKwZBw4JVmsu+roPUtrzYW6Y8dVOXXeKRiMaxqi83U4AvpN5vQc4p8Q+J5SsW3S+GwuOVfPr+fI2zOyUHscvRdkIB2cWFO8DvmVm0x/dOiXtPlmtuXDvtTocOkRy6HD7RLtZICc8a1yti8Lx9Dse5idkWpSgrLeYROW4Y6nf973+jrdRZN6HWRUe8PN/0lBDZiRLy96JofW9SD3bygiKKlUXhdqPJZVGkRt79FbVmrojfIDsw9vtaEnZuYqXm9nl6WEK9s8fudM+ZermGepYks40sy90PUGJfcp2u70bOBP4kn/9VOBW4ChJF+bitE0n6UzoWg0dOuy62NKYVqurm9/iGYQy8eCKhglTke/WxZckUyNAsyw4WWy1hu3f79YzgUFdtmC/XvIBTpUq0Y5tqFJxjgl9JCdUpDYHhvZGjtqFu/SeD5rZWR227QEek3l9IvDdkvvMlahb9nx7/HqvY/21pOdRLFYpfwU8vVsjyorPN4HXmdkd0Byw+h/A23ApU6defCwxZ+kkS1i0DKnnziiDEE4iJeJddc9cmouKnS+H9oR2uW2F7Un3mSIBCjhhabpO5+cA1WstR4QSXW3R1kV0xC5scR7V61BWfLyDw5RNeP48sFvSKcC9OGeAn83tcw1wkR/TOQfYZ2b3Sfp+ibp5rsGFLftjXMi03cDn/MT+Az6IwH8Crwb+tKD+Tty4T7cbR88nsrLi86RUeADM7E5JTzezezaNA5yfv9Apd/1mpJnnp4P10iuqQRsdvgc9M6mCC2LZ7HrLiJnfv3L8cSBRv7fXA11goyl8WJP3gKtW3Dhfo4HVe3y3Gg1YXkErq677rnQDRpQJtSSjmGTqw5VdhIufGQPvNbM7JF3ot18GXAu8CLgbOAy8tltdAEk/hROPY4B/kXSLmb3QH/sDOIeGOvAGs+a79su4WJ4LOEeDvLMBZnby8FddXnzukvQe4Gr/+meAr0qaxw1SbQ5mQHDK0lF48pZSVrjKPoikwpO1jtLB4XTMICXtnjNrihAS9T33ljtXYENJx22oVFwqBtxE1ORw0rlXwRIaBw+iZRebsu+x1vX8HY8oaKiZXYsTmGzZZZl1w+VVK1XXl/8j8I8d6rwdeHtB+U3AU/pp+6CUFZ9fAP4bLpK1gM8Ab8IJz4+MpWWBsZJaPU0PtAxrhKeb4OT3S+kmRAXjPeZvGFLU8oJrRozOuYGbUTnh+GAJTSp+wF/VSut7EMcwPw+REKBGA1ZWunZrT7yDjzG9EcwngLKu1ks+dMJHzOyu3OaDRXUCk01WYLpaOUXrfZ5HGdEoFDGpKThNEUq74HJjQEBThCrHH0f9u0Nlcw+MGMUx0eKiy0g6Nwe1mgtP5cWGSqXdJV8bMG9uVK7WzHZst2EpNaNQ0vnALcC/+ddnSLpmnA0LbB6saO5QkeAl1rwxmSWtm1Q+gkCSuDEBMyrHHTvm1gf6IVpcREcdSeOU46idfAzs2gGLC0gRyb79JA89THL4sAtDtBHd3KMO2zOCCAfTig/583OSfte/PknS2WXrl/0kfg83C/cRAJ+f++Q+2xqYEKLFHllf05v8MHSbJ9RNgKApQpa6v6eiZLa2fj9OEYHxkXa1bdtKcuQ2Dp24yPIx89hcBbbMY40GyeqqW5ZXSGr19YsW4j3gxhIvbobFBzcF51nAK/3rA7hQP6UoO+ZTN7N9m8azLTCcuJTJXJrvMis6f0FXWtu+hV1xufq4oKRhrs0EELfSkMztq1Pdv0KyOEd8eAVbySS23ACLR3EuLcoojmkz3+12jpmdKemLAGa2V1LpyMJlxed2ST8LxJJ2A78C/Ef/bQ1MAsnhw0QLC31UGOAXlg2d40Uoa/3Ie7C1Qu10iC2XNZgil7ag6Rm3iWO9Th1pN2m9TrTvMFsOrTjrdb4CBw6SLPfMqjw2lIqiXFbTkTLDyeSAmg9IakCaCLT0k0XZn+9/x0VAXQH+DhfB9Nf6a2dg6siMwXQkG4I//0Sb1s92m/nFksQtfuymFRG60epOy4ztNKNF+3ApaZdc2i2XpmYIbBxWdxFB2H8A9h+Eg4fQ/Q+R7DuwodMYbIwTxlPrp9eySXknzpX7UZLejvOC/oOylct6ux0GLvFLYBboJjrdbiRF29rm8OQiIkjFFlFbW7IZwiQbswAAIABJREFUVNNjZU7pM2/GRx5B4+G9ndsWGCuWmMvjkyRNa8NWVjY+2sA4hW/zCktPzOxKSTcDP4abgvOSNJJ2GbqKj6R/psvba2bnlz1RYEoYVHS6kdZL5/BkifK7rj1H2zykIiFq9A4RFFgHsg4ikYs3nBTk7Nk0bG6rpieS/sbMfh74SkFZT3pZPv/H//9p4Fjgb/3rV+LivQWmlX7GcbqITtkB3LZgj63CgnZkrKPmObzDQdYqyjsrzHgunUnBfR8aYOr+3cjlfppaZvtrd3r2hR//eUbZyl3Fx8w+5Q/6NjN7bmbTP0v6dIdqgWkhm4OnEx1uDv16DWX3b6Yi79RFlz12xmtuTfcctIlQ48CBvtoUGAP+M7Wk83Cy4rgpPj1jvE06Myg+ki4GfgtYkLSfVoDRVeDyjhVzlPV2O0bS4zLZ8k7BBasLTDMDdLF1FJ0yT7DpDSd3DGUtnWwXXbaNqQj57rVUiJLDh3ufN7D+WNJKBJdJCKdKtT3sDlMQRqcLs9jtZmZ/APyBpD8ws4sHPU5Z8fl14JOS7vGvT8bl+Q5MKcnKsktzXUSBkHS1dMp2nWT3y0z4s8TaBSjdN7NP0kc+l8CEUJBTR3HsQu9A28ThngI0TEbTcTKD4pNiZhdLOgKXkmFLprxUr1hZb7d/8/N7nuSLvuLzfwc2G/3+wIdxQsgJELRbQclKH6H0AxuPIvf5pRZuPj2JIohbFq1lveB6iIvznou9a/6EiFBwOPgl4FdxSeduAZ4J3AD8aJn6Xef5ZNNnm9mKmd3ql5Wiffps+Nsk3SbpFkkfk3R8ZtvFku6WdJekFw5y/EBvkpXltUsaAiW3WL0gjXiPm4Al1lw67LDmGNlzBqaHNFGcKlU0VyXatpV421aXzjrjXGCrNWx5xcd28/O6SohJc57OpAhPymyH1/lV4L8A3zKzH8FlLi0daqSX5TOSdKkdeIeZ/Q6ApF8Bfhe40GdJfQXOk+J44DpJT8gkOwpsIFkBUiacSts+BWJTZNlkNm7+jLGbGR83TZEg9onj5qpg7vYSWeJcrtPPeSXBvKXT/K70EpVJE52UzSssZVg2s2VJSJo3s69IemLZyr3EZyTpUosws/2Zl1tpfYwXAFd76+obku7GBTW9YZDzBMbHIIKRekFlRSgIzyYhtXCqFRfFIo1OkaWZaXT6P3Mx291uwB5Ju4B/Aj4uaS9QOslWL1frk4drW3d8SIZXA/toJaU7Abgxs9seX1ZU//V4x4ct9IjUHJgMUlfc6b/3BFIyqasVR9BIsEaC1WpufYTBPCcKA02oQbYemNlP+dW3SPoEzlj5t7L1xxqaUdJ1km4vWC4AMLNLzOwxwJXARWm1gkMVfnvN7HIzO8vMzqoyP56LCAQCvbHEjQuu1rClpea4zrhiqk0MMzrmIymSdHv62sw+ZWbXmFnpwdqxio+ZPd/MnlKwfDi361XAS/36HuAxmW0n0ocpFwgENg5rNFyunnrN/c97vK0H48rdU8SIxEfSud7B6m5Jby7YLknv9Ntvyzp6daor6UhJH5f0Nf//CF/+Ku/olS6JpDP8tk/6Y6XbHlV42S7Pya2STir7VuXZsKD03nU75Xxa8YGuAV4had5PZt0NfG692xcIBAanmSiuKNr5+E++buccRVRrH5bmXcB5wGnAK73jVZbzcPfC3bihhveUqPtm4Hoz2w1c719jZlea2Rlmdgbw88A3fYLQlFel283sgS5NPw64Q9L1kq5Jl+5X26LUPB+5LHKvAh5nZm/1anesmQ0jCpd6z4gE+BZwIYCZ3SHpA8CdQB14Q/B0CwQmmCIrY1K900bNaLrUzgbuzkSQuRrneHVnZp8LgCvMzIAbJe2SdBxuwn+nuhcAz/P13w98EvjN3LlfiUuTMwj/c8B6QPkIB+/GicSPAm/FpUv9EM7HeyDM7KVdtr0dePugxw4EAmMgLzJ+orAiObf7SsU7GzScA8JmF6DRORycAHwn83oPcE6JfU7oUffRZnYfgJnd16EL7WdwIpXlryU1cPf43/eCt4Y09ueglO12O8fM3gAs+5PuBUqnSw0EAlNMOo8nnUgax22RDNL/iiI3vydXb1NTfsznaEk3ZZZseLIyTlad9intoJVH0jnAYTO7PVP8KjN7KvBDfimVHmEQylo+Q6VLDQQCU0omZI4iOesGsNUaWK0ZFseSCGp1iOQsoagCcQxpFtF0DGhMbXSNWv9bUh/zfB40s7M6bCvjZNVpn7kude+XdJy3eo4D8uM3ryDX5WZm9/r/ByRdhesSvKLbhQ1K2ceSonSp/2scDQoEAhOGt3qoVFAcEy0stCycNEpBGr2g0YA4RgsLaGELmp/vGAljUzAab7fPA7slnSJpDicK+YH7a4BXe6+3ZwL7fJdat7rXAK/x668Bml7GkiLg5cDVmbKKpKP9ehX4SSBrFY2UsoFFh0qXGggEphRLgNj/90n7fAidZjDQdJKpDEwonw5jXdq4AYxoDo+Z1SVdBHwUiIH3eser1AnrMuBa4EXA3cBh4LXd6vpDXwp8QNLrgG/jxCblucCe1FHBMw981AtPDFwH/EW+vZK+1OHK5ZpkTytz3b3SaB+ZefkAGRNN0pFm9nCZkwQCgenF6jU0N+dC5iQJrKxi9YL02F6ELDHXpRJ3zjCbWkNj7Y4bM6J73LF+MLNrcQKTLbsss27AG8rW9eUP4QyGojqfxEWhzpYdolwm0p8ssU9Pelk+N9Ma1DoJ2OvXd+GU9JRRNCIQCEw2lhiiAYqw+lLP/E7WaHS8MadOC+6FTXU201kMr2Nm3xrFcXrFdjsFQNJlwDVeYZF0HvD8UTQgEAhMPi5IqEgDgva0WNLuukg4zVIz/pslmZmXmUynLhr2HEguTM8kplDIswlD5/RC0gG6d7vtKHOcst5u/8XMLkxfmNm/SnpbybqBQGDaST3a/Hrv3Q3qteLU6ZlApGuy20qtLr5piAk3g+JjZttHcZyy4vOgpN8G/hb3dv8c8NAoGhAIFBHNtaaRhcRyE0I/VkjWCaHksazRQHW5iNhTYvXMeEoFAPzk1Wwa7W+XqVfW1fqVwDE4d+t/Ah7lywKBkdMUHj9JMZrf0r1CYDLJRD8oTCKYx5LJzVjaiRmNag0g6XxJXwO+AXwK+Cbwr2Xrl3W1fhiXMjUQGCttwtMsFNHCAsnS0sY0KjAwiuPWZNMynm1esCyJvOv2ZIvQLDocZHgbzmPuOjN7uqQfoQ+jpGxg0U9QoN9m9qNlTxQIdCPbzZYXnsDs4MZ54ukY72Hmu91qZvaQz+0TmdknJP1h2cplx3zelFnfgsu9U+DoHwj0jyrV4g1Z4ZGIFhdJDh9en0YFhscSklodJRnPthJ1Cp0RJpFN3KVWkkckbQM+DVwp6QH60IWy3W4354o+K2moiKaB2aaj4BQFolRLhKJFly49iNAEooio6mO6JeYcRSzpfx7PpItOlhkUH0nzZraCi4a9DPw6LuXOTlzWg1KU7XbLRjqIcLNgjy3d2sDM0y2+V3Mwuqi7Tcrt6/aJt20DoHHw4OgaGRgKxTFaXETzzlWag4dIlpanS0z6QMxst9sNwJnAZWaWRr1+f78HKdvtlo10UMd5N7yu35MFZoeuwSQzIlPoBVVQlopOnnjbtiBAk0S93jIGGptTdNqYTfGZk/Qa4NmSfjq/0cz+b5mDlBWfJ5vZcrZA0nzJuoEZYajoxV262zoJT0qwgiYDazRIlleAlVaU682M0RrPmi0uxHWz7QJenNtmwEjF5z9wZlaWGwrKAjNEc9ymn26VXsnFct1ta4RHBZaSD14ZrKANpsz4zgbm3hkHs9jtZmafAT4j6SYz+6tBj9MrqvWxuDStC5KeTiuI6w5gcdCTBqaLNjdocqFSoBVav0/WjPV0GOfpWJbdFgRo4lGliuJ4c6XZnkHxSTGzv5L0bOBkMlpiZqWSz/WyfF4I/AIuO94fZ8oPAL/VT0MD08WaqAKjuFEUjfV0EJ42i6eb8GT3CQI00Vij0frcczHd0myp6xLd4Fk/AP/xDyM51CxaPimS/gY4FbiFNOKsk+PhxcfM3g+8X9JLzexDwzQ0MD1ECwtupUN/tpuB3uevzotMGWun4xhPUXmSvYkFAZpo/LyfjuJSIvpB6X3XixkWH+As4DSfa6hvenW7/ZyZ/S1wsqQ35reb2R8XVAtsRnp1rRVtz43vDCQ8vayedP9kbduCAE0gRd+h7MTSQehHlJ71A0OcKEcILHo7bsrNfYNU7tXtttX/31awbbbf9k1K0+oBJw79WD8dnAna3KnLdLNlyosPmHbdFLQtY/0E1pkBx/5K0+PhphRlunDLHoqZj+12NHCnpM8BK2mhmZ1fpnKvbrc/96vXmdlns9sk/WCfDW3D5wO6AEhwKbp/wcy+K+lk4MvAXX7XG7O5hAIbSJ83lzVzeMo4FnQrTzFz+2SFJoqC9bORKPKu9vH6RaVOE9FlX3fj2We43QbQrM5tmOkHnbcMU7nsx/CnJcv64R1m9jQzOwP4CPC7mW1fN7Mz/BKEZ51os3pKoEhdl9zO/iQ9JpD2+2RaYv80JE9gfCgSqlbQXLVc+oQyx6xUO4dhSulD5CwasfDgut3KLD2PI50r6S5Jd0t6c8F2SXqn336bpDN71ZV0pKSPS/qa/3+ELz9Z0pKkW/xyWabOMyR9yR/rnVLnH5iZfapoKfve9RrzeRbwbOCY3JjPDmCIGYVgZvszL7cSuvE2lGhxsfgprkvXWyk6dY0M6tWWklo/nY4920+k6441GrC66jzWRjHxMvWAA6yXxV0ms+pzzhi+TWsOykjuWpJi4F3AjwN7gM9LusbM7szsdh6w2y/nAO8BzulR983A9WZ2qRelNwO/6Y/3df/gn+c9wOuBG4FrgXPJ5eiR9Bkze05BOu2+0mj3eg6Yw433VIDtmWU/8LIyJ+iGpLdL+g5utmzW8jlF0hclfUrSD3Wp/3pJN0m6qdbqcpwcBumT3gD6sgx8gre2pdu+bScaYB5PFLWWPH0ITLB+xo81Gli9NpouN+8Z19U7rh/Sbtp0GRFKyi09OBu428zuMbNV4GrckESWC4ArzHEjsEvScT3qXkAr5tr7gZd0vRZ3vB1mdoP3YLuiqI6ZPcf/325mOzLL9rLCA73HfD4FfErS+8zsW2UPmiLpOooDkF5iZh82s0uASyRdDFwE/B7Oc+IknyfiGcA/STo9Zyml7bscuBxgh46crEfd9MY77kHYIWnelLvdyHtZPwOKbNPqyd8MeoTTGZSQkmHKGNHvJnnu09u6vmyEKaJG5HBwAvCdzOs9OOum1z4n9Kj7aDO7D8DM7vPprlNOkfRFnCHx22b27/5YewrO0ZVB02iXDa9zWNI7gNNzJ+maTM7Mnl/y+FcB/wL8ng/VveLr3yzp68ATgJtKHmsySAdDJ1R4+rYEhu1+6+Xd1tyvT+EJ3WuBlALng+S5T3cvx5GT0Ojn+3e0pOw97HL/8AxQ1Lr8gTvtU6ZunsIH/H6PJel84I+A43FOY4/FOYud3uP8QHnxuRL4e+AncUHlXgN8v2TdQiTtNrOv+ZfnA1/x5ccAD5tZQ9LjcH2c9wxzrg1jAoUnWlgYvNshFZB+RajbAHTaljFZO2uaMsnWzyaLe7au5K1vdbCqu40VDnLa8j+FB83srA7b9gCPybw+EfhuyX3mutS9X9Jx3uo5DicQdHnA3+Prd2tHlqHSaJf9xR/lA8jVvEfDL/qTDsOlkm6XdBvwAuBXfflzgdsk3Qr8A3ChmT085Lla9Bqn2GRECwttCzC8tTCoN1Pe6hnmJpAkha7VZQjjP5uMDr/n5LkF4+kjFB6g5XTQa+nO54Hdkk6RNAe8Argmt881wKu919szgX2+S61b3WtwhgL+/4fBPeB7RwWyD/j+eAckPdN7ub06rdOBmpk9BDTTaAOlPTvKWj41//8+ST+BU8MTu+zfEzN7aYfyDwHjCeWT/ZKuV5dY/ocxwDnXxFnrunPJH1f2CXAQMSrbDVemPcNYPUMI0MRZQMHiGYz8fB8g+eGCe2CnTqoBGVUyOTOrS7oI+CjOi/i9ZnaHpAv99stwnmcvAu4GDgOv7VbXH/pS4AOSXgd8G3i5L38u8FZJdVxMtuwD/i8D7wMWcF5ubZ5uOcafRhv4fUk7gf8XN79nB/BrZU8ysYyzm6PoaazP8zSjSRf8uIor9PnLGrYLopcADdLdlm3PoBZapp51EaeJFKDAYGR+W8nzOmR6GfW4j9nIxhzN7FqcwGTLLsusG/CGsnV9+UPAjxWUd3zAN7ObgKeUbPYFwBLjTKNtZh/xq/uAHwGQNF3iU8YleJQiNMSx8ikMxtpNOIouuDIWUJnEcHkhHGS+Tp/7BwEKDMMsh9cxs0N+NcEFoI5x3X5Xlqk/zF1tTaDRqWcCxoPWCA9MZ3dMv1bYoKF2sgwopGEMaPOQt3pMav7PLqNiVBEOpglJOyRdLOnPJL3Aj0NdhHMM+3/KHqdst1thG4aou770Kygb4HVUKDopk+4gkbd+egnPODzchrTgggW0eRml2LQfmOGmH0wvfwPsxWWz/iXgf+C87i4ws1vKHmQY8ZmOd32YG3cu+dk4ctJ3FZ2UsmM+iQ3uiTYsJcZ3Ona59bo5DND91m2spxNBgKabrNXTUXAiXCfRqJiOu+CoeZyZPRVA0l8CD+LmDR3o5yC9YrvlY/c0N+G8ISabEVkMaYwpF7XXMYwQ9QyWOAsM627dIYr1sAQB2hzIbK0ARbn/IznP6I41RaTez/j5mN/oV3igd3id7YO0bBbIClFKkSAFockxri6QDINYPVmCAE0pI3alLnfOmVSfH5CUhjsTsOBf9xVYdJhut5mgn/DwTUHqw+KyxEYWgn7q6RZBYVSRqvPHWAcxDIyf5Ic7uFePE5tNbzczGyqjQcrmFZ8RdLlNlCiUHfeZNIrGe4q63KLc+jgGcvuMgq0oConopoV0zto6Wj9ukulMWj4jYQrvZhPONArELDDgTSLeVpRBPjANdBKGkXq/JSWXwBo2p+Wz2ayelGm1fnpR9F6Py/opSdZSSwUoWEGTif3Q0/1KD+snwT1ujzLETrB8BmZzis+UMdJxn410ty7C/zhTJwClXmrrFMV6lMTbW/43jQN9O/cExogslzahgwCN1OoZUSbTWWXzic9mtAyyTJP102a5OFG0RgMk98QoeREarVgqiob2eFt70PY2ZoWoiKkWJ0VE1crosoiOG/9dkoGRDZgLIuNy7b9qo8vtY2g2J5mOhCm5i60/I8lDv1Hn28gfRGKtxZL2pVluThzMnBXUWI/o4h3uOGPqNuklTuvGACGjFMdQqRBVK5P/oPNsH706+znmPtM1XWMjtX6s3BJYw+ayfEb8QxmoO2wIy2Tq3a5TwcmQiqoitbYlEUS0rKB053iYaBQlXLGL9hmjq3W8fftUWkDWaMBK0nogmtQkd8/6ATBDibAIt46cZVMQsd3FdRvh+WfU1XpUTPhjzcYysBAM8SMtZQGVOf56Wj9ZSwd3DemS0vY6bwU1GhtvAY26jmddLaBOaTz6/T5a4j6TQequJx2+4s2oA2atfTK9u6MVoGD5DMrmsnwmiXGPzUzC2E9WTMgJZ9FNS1HGElprBQ1lAfXzAx+VtTNpThOjEIpuSRaz37cNFiWd8zSnK/5P1voptni81TOOTKaBgZiwX89kMfQ4zIA/0PUebxqIzJjOGqum03Vntq2xgvIW0KS/B30IT7xz5xgbMkKKHmZyZZPQLayznuq62PJjO5mvXVvMtdH7tGTOY6WWwFqC5TNuNtJCGYfbda57rUmHsZ4szRuXf0+aY1yWtFtAZhDHkxv6ZgCLJ965k8a+fWNozAjJPzT4760iNT9PKxjXWy+iM08Hqc3YaAYQzaeFT9ezX6Gm9TOiBhnQCMIyKMHymWbW8yZQMK7TbEOmHfmxnixFFlKzLGsBJQk0GpPTVz5MO6TmEu/aNbo2rSM9u1PXgejppzXXZeYmjKafS6Z5eQeAti63ESPKWT3B8ikmWD4Tysg834axfkqO6fTbTdh2bW1WkLeAMMzMPaA25wJNEEWTZEtaafGuXTQeeWQMjRozGyQ68VOf1LJ20pt40XudDgA1M5fivvd5R4O0fFQEYRmYCftVb1I22mOo3/GT3DydTmM63ayctvk9hZtzxyTTpeOtIGs0WuNAk0yf3YPxrl2TbwVtsDNL/JQnEj/lia2Cfr4C3tIcq6NByoi83SSdK+kuSXdLenPBdkl6p99+m6Qze9WVdKSkj0v6mv9/hC//cUk3S/qS//+jmTqf9Me6xS+PGvo96kAQnwlmZG7XUE6AhhWdToLTRYgKBSh9nXVEyOdKGqMgdcy4OorzZ566J1KAiialrvPDU3z6E3xbWiLiLJhu2XLdP4uyFk6rfCwYIwksKikG3gWcB5wGvFLSabndzgN2++X1wHtK1H0zcL2Z7Qau96/BZR59sc9G+hpcWuwsrzKzM/zyQJm3YhA2XHwkvUmSSTo6U3axV/G7JL1wI9s3FQxzc0gFJzOvY2DRKdvWTs4JmfM3X+fHgSbRCur0ZBtp7ZIrj488gvjII0bbnlRABlkmnZwYZWO1WSzyVk9zW6bLbZTjPyMa8zkbuNvM7jGzVeBq4ILcPhcAV5jjRmCXpON61L0AeL9ffz/wEgAz+6KZfdeX3wFskTQ/+LswGBs65iPpMcCPA9/OlJ0GvAI4HTgeuE7SE8xs8LzVg7ZvAtxKR0o6/pOzNtxqf55rRfv1Rc4LMD8OZEnqaZXxhEuSVmDSMVHK6imiGatusO9MKkCNh/cOdv4SwpG+v9moE+n73tXK7jb3ZxyoXUSaVk+R8GS71CQsUsvqgfF2uWH9fBePlnRT5vXlZna5Xz8B+E5m2x7gnFz9on1O6FH30WZ2H4CZ3dehC+2lwBfNbCVT9teSGsCHgN83G8/A1kY7HPwJ8BvAhzNlFwBX+zfjG5Luxqn7DRvQvg2ntONBWZfuTPfVmhvOeohO0XF8u8s6IiA5kVgvV+x+zjOke3vWCiolRJ3m4XT5LjQn+CpCafLdTG7KonTw6ylAjdvvIn7ak9uFJoo6io5FeNFx5ZbdZ5wY/TgcPGhmZ3XYVtTS/IE77VOmbvFJpdOBPwRekCl+lZndK2k7Tnx+HriizPH6ZcPsbEnnA/ea2a25TZ0UvugYr5d0k6SbaqyM9MexaayeXHDPNaFvcuMxRaFx1hxrHG1srvZ2RGgLTDpJpO0Z0QTZ+Mgjuo8LZQRGkVxA0Ez3mSI1l171svulrwvPN86uuUz3Xy/hsViZJWoTnvxdLdvNNnKX6xGM+eDucY/JvD4R+G7JfbrVvd93zeH/N8dvJJ0I/CPwajP7elpuZvf6/weAq3AP/mNhrJaPpOuAYws2XQL8Fu2K26xWUFb4a/Zm6+UAO3RkOgsuc6QJ6sNej4CjvSyXbgP+JY85Ngq64aD1lN7WDdfA3aCyVtAIGNVxgDYLSF0+d+v1/kZr5wc19u1vritSm9hkNjTrg7dsvEVjiXmrJ3bb5fupzJCs2e3ZtIay43DjoOj9iYDIi1DcLjqtyaKZsrSpmcmlbV1uo3axTps+mgegzwO7JZ0C3IsbdvjZ3D7XABdJuhrXrbbPd6V9v0vda3AOBZf6/x8GkLQL+BfgYjP7bPNapAqwy8welFQFfhK4bhQXWMRYxcfMnl9ULumpwCnArXJfjhOBL0g6m3JPASUbsIHRBYoYowB1tFTK7lui3tgpeH/aoiB0ioqQTVQ3ajp1uZXyiHMCZJZ0FKC0vKMIFXxW8c4dbtPBQ+4Y+a625utc2yvu554tbXZfptZkBBAjckIVjUmAOrwvyS1fQc843X2+fizNmpZQpnstLu5e62rhjFKDRiA+ZlaXdBHwUSAG3mtmd0i60G+/DLgWeBFwN3AYeG23uv7QlwIfkPQ63Lj6y335RcDjgd+R9Du+7AXAIeCjXnhinPD8xdAX2IENGfMxsy8BzcEvSd8EzvKKew1wlaQ/xjkc7AY+N/jJ2scVylebzPQGXW8Ag4pNl7rrTsHn1SY45KwgPxZEpJZDwgD0VW+Ac3QToO7nUsduvGjbVuzwYfdeFQjOmnGxtN15UWo0aEu+N4zIFLlpD/oA2HQ0oKvwWHOcMNfuMUU2aGI2skjsZnYtTmCyZZdl1g14Q9m6vvwh4McKyn8f+P0OTXlG+VYPx0Y7HKzBK/4HgDuBOvCGkXi69dkdN9HC00MoSj+hTorgFNGxGy5vBflryDkkjIRRHCfT/dbLAiq0fnp8llpcxFZWmm1tE5y2dblzZ0XKDFutQRQhMo4Gkdz7Oaoo2QPXTbvbsl1tLeFpCZLf39QeUHQ9mLRxxyliIsTHzE7OvX478PbxnbD309hEWD4lfvx9dYVMstgU0a8VBM2uOFc4hs9vEMsqJ0DQfRyo/XydLZ8m8s4GafuKBKdSyYx/+P+rq2h+DqvX212Gmy75re62fr5nqcNC9mFpkC47y1g96fhNyxGBliD5uTtqmHP/Sk+bfvxjGu9xJwniMygTIT4bQocIvu27WG6XEXyBBxCAgfvap01sOlHSCnJlaVccvWfFe0p3uQ0zppRzwc5aOV2FqMRnr7m51ous4KRWQxRBtQJR3F7Rx6hL02arkWCrq06M/E3eatayiPr5PqVu3M05WwN0Xvznl+BZT/PH82M7fgzImg4Fmf+xmgJUxEgjWoNzg5r01B8TzOyKT54S3XLdRGBYYdqo3EFTQy9nBCjuiispQGvI1xljcNN+nA26oWrGukktoVR8ohjmq64sHadI/M+/EoNVoJG493OZZkijZkbT0heTsKaTfNjvZjNYaMZzLWPdpeM+yk6iTt+6fp0R+sI2/+9ujATxKaKEVbS2yjo9Ac3yl71DN5wr6u2QMPZJqb2OP+QE1J54b7ZshIDmdccRVkktwsjcycWOAAASAUlEQVS7CM9BveEyx1Zity53E7elJbfvIN83S7B6B0efbl3eRVMBml2o/nXWicBbPVkriNT5rcvbPNJxodDtNjBBfMrQ6wc4CnfuWRaVfulLhBqurIMIjcw9u4ywjXkM0ZaW0Y5tmbGdjOu15CY7ViKQuXumKiiOYHnFWUPzVWAOGg1kCaxGREBSqw/+/Syq1+exmmITFZRl94vkRDU/TpZ730dm+Rgj83abRYL4jIIgHBtDGREaZVfcFGAHD6FdOwuvT0mCZe7gLilb5iZdq3vrLELz8+59rVSIGw1sZWU4ERoHWc894XrBUusnK0A5izNYPpNBEJ/A9FNiPKiwKy4egaPBOEWsOX6RucF1nPRaMJ8nTyNBtQYmoSRxT+3psVN37WoFtsxDvQ5xjFZWnEUBRIqweq049tt60e3tXvfnCQviMwRBfAKbg0GsoAbFIXo20CpaM98nEjSyEccL5jHlu/P27oMjd7WXp0KzWkN+jk8zVXlibWnLbcdWksXt1BerzH3vANp3AGPZWUqp2/R6WUBlP4tOITYH3a/sOScxxceUEMQnsLkYJEKCD9EzlvA8Jcm6W68JuZOKRT/HyYxtmZzlQiNxFk0jad000+Nm/mv/IWrHbmPf4+aoHj/Hzrvmib/3MNY46LzoavVhLnU4suKRGMQqFuQClHbLjbQ9wfIZlAkKfBYIjJCCIKsdo2VnI2WPiy7OBr3ivpX11Gs7zsOPNKM+E+EiQ8eRE+UkgVrNLWmWWMl1p9XrUG84j7cIlo4R9Z1zsLAFzVW9EA74Pg2S1K7Dzb1w3GYjdGBEabRnkWD5BDYvvaygonxBsDYmWidG0D3XK9JBx7A7KUWi1jYQn86H8ZfciKGm5rhNdkKqqlUnTBLxoTqVw3NsvS8hWvVddlu3ouUV9372K0D9eoQW7Z8AcbsFIz+lJy1ritJ6TH0w29jxrykniE9g81OQtM69LHDJzo4DxXHR0TLHLRlHroPVUyg8+X0TcwIUJe7m28kSKIhorQceJjnhmOZrGk5Uhfd8W1nx6bwj151WrUClgs1VkRm7vnqY+rYqFolkxwJaqqGdO4hWaySrq+UEqCjdQwc6zZVTajzIeehZ2g3pFcgJkCEElhGg7ETTcc2xChEOBiaIT2B26JS6u8gKisAajfUdB+pixSjJCFB+W6dUCj66QTMeWoQbI5FQvdEad5+rOtHZUsW2zFHfNsfSMXMcOjZCCSw8lLBtpUGyfZ5qNUYHDxHPVUmWltuf/HMin801VIamI0iWBCwN1mBOYmgKkRespqt1xrU6qwl5J4Pgaj0RBPEJzBadUndnnBHWOCJAK1xNv5S1eno9lUfKCFCJG54XoOi7D1J/7KPdzTcGswhFEZEZqtfdONCWeZLFORpb56htq7K6M+bQsRGrO6F6AFZ2RWz7jljdWaFyMCaaq2JLDTRXRXU1x86ajhxpJtT0GssGyO0mVKmVWdD1Bt7asdb/tCzvYDBSpwOz4O02BEF8ArNJxgpqi2CeZvHMjwNB05Loi4LunlIRrfPdeX7mfkcBKkog58euKt9+gNqpx7bGRCKRMEe8sAUOHfZBRQ01jKhuRHWI6hAvt8K/1RcrbP3aI876qHtvt0bi6qbnW3HjQYpj55wAPjZcwd0+7TrzwlXULZec+eTWtStr8bR3vWXHfdL3RUUCNI6ut2D5DEwQn8DskhMgaM0JagoQNMeBFMfNSNCtY/QY9ylzs2tLf91h//RmmxUgOghZ/hgS1XvuZ+UJx3mHChfnwBbm0OoqHDiI6nXixgJqJCCoLVZJqmL+ESNehajm34sDh6Dhx4sAxVWoVrDlleYYmaoVtGUeq9Xb8gQpl8zOGg3XFZh0sZD8YE/LoaDV9SYAyzw8+O62jp5wIx/yCQ4HwxDEJzDbdBkHys8HMhrNAftBo1x37W4rEI1COmUoKMpamtk2//X7WX7isVg6UL9QJWosoMNL2CP7Ua1O7AOCzi9EJNWY1W0irsFCNROWKI7WdENGWxfdjdiLuKVzgXyiumabsnH12lJ3g8vcXHBZ6bhPruvNncirUmIoyZRlrnss83uAkFJhOIL4BAJ9OiI0+/rTG/wosqeWEZ7sjTuO2zOPdjpm1t0aSGJ3M06IUCNC81W0cwd2732+iogqEZWleaqHvFUosXxUheqBLcT7DqFKBRa2uMmmcYy2bHHzgOp1WFr2yems1YZmoruWg4T5uUWdrtWecmrLykmdCZRaPM7CSbc1x3ms1eXWFJvsZ5PpihsZkxTrbsoIk0wDAXA3kcyNpNuEVEvz3Aw72NxNNJr7RO1Lur2TK3i32G8Si1/+HlZxk0+TakyypUqydQvRo47GVldddtPVBtFKg3jFd7nVDYtFbcdc0zNO8/Noy3xr4mkaOSGN8hy59jWXaqW53pZIL4p6exRmRKVJ6jdirX2U6XJTrvttHPN/DJqTl3stvZB0rqS7JN0t6c0F2yXpnX77bZLO7FVX0pGSPi7pa/7/EZltF/v975L0wkz5MyR9yW97pzS+WFNBfAKBLJ0EKBvTLI2IkEYH6EQ3R4NO3W3dus7SsnReThR1cM/uFgg1YvGO+0nmRDIX0ZiPSbbOkRy1g+jE4901LS0TrSZUlhIqS0ZcA8yobY1ZfcyR2M7tToSqVdeOet2lZajVYK7qhGl+3u2TLpVKK99QNlhqRkjThSc/bs2cprYuNXNddkrX0y43yy1k6owDH+uu1NIFSTHwLuA84DTglZJOy+12HrDbL68H3lOi7puB681sN3C9f43f/grgdOBc4N3+OPjjvj5zrnMHem9KEMQnEMiTE6CsCBWG5EmS0dzgOglPtvusLVFc1F1oOhGJrXc8gFVEUnUilGypkmzfAscdA40G8eFV4pUG8UpCvGJENUBQX4xJdmzBti24rrdUhKoZMapWWkscu8VbR82uwmz7M4vtPqn9LbHMnB1v3TRJvAgl1hSjpig16/f/9vTDiCyfs4G7zeweM1sFrgYuyO1zAXCFOW4Edkk6rkfdC4D3+/X3Ay/JlF9tZitm9g3gbuBsf7wdZnaDmRlwRabOyNk0Yz4H2PvgdfYP39rodoyIo4EHN7oRG8jGX3/+frG+Xfvrc/1f77Ltu2M/ezE3r9tn/9hhD3CAvR+9LvnA0SV33yLppszry83scr9+AvCdzLY9wDm5+kX7nNCj7qPN7D4AM7tP0qMyx7qx4Fg1v54vHwubRnzM7JiNbsOokHSTmZ210e3YKML1z+71T9O1m9mouqSKxlXyjz+d9ilTt+z5BjnWwIRut0AgENhY9gCPybw+kbW2Z6d9utW933el4f8/UOJYJ/Zox8gI4hMIBAIby+eB3ZJOkTSHcwa4JrfPNcCrvdfbM4F9vkutW91rgNf49dcAH86Uv0LSvKRTcI4Fn/PHOyDpmd7L7dWZOiNn03S7bTIu773LpiZc/+wyc9duZnVJFwEfxU2ffa+Z3SHpQr/9MuBa4EU454DDwGu71fWHvhT4gKTXAd8GXu7r3CHpA8CdQB14g5mlbpu/DLwPWAD+1S9jQRZiEwUCgUBgnQndboFAIBBYd4L4BAKBQGDdCeKzAUh6r6QHJN2eKXuHpK/40Bn/KGlXZlthKIxppcP1v81f+y2SPibp+My2TXP9Rdee2fYmSSbp6EzZprl26PjZv0XSvf6zv0XSizLbNtX1B1oE8dkY3sfasBUfB55iZk8DvgpcDD1DYUwr72Pt9b/DzJ5mZmcAHwF+Fzbl9b+PgpAlkh4D/DhuYDgt22zXDh2uH/gTMzvDL9fCpr3+gCeIzwZgZp8GHs6VfczMfBx6bqTlb18YCmPdGjsGOlz//szLrbQmt22q6y+6ds+fAL9B+6S+TXXt0PX6i9h01x9oEcRnMvlFWi6OncJqbDokvV3Sd4BX4S0fZuD6JZ0P3Gtmt+Y2bfprz3CR73Z9byb68ixd/8wRxGfCkHQJzvf+yrSoYLdN6R9vZpeY2WNw136RL97U1y9pEbiElti2bS4o2zTXnuE9wKnAGcB9wB/58lm5/pkkiM8EIek1wE8Cr7LWBKwyoTc2G1cBL/Xrm/36TwVOAW6V9E3c9X1B0rFs/msHwMzuN7OGuUQ/f0Gra20mrn9WCeIzIUg6F/hN4HwzO5zZVBgKYyPaOE4k7c68PB/4il/f1NdvZl8ys0eZ2clmdjLuhnummX2PTX7tKWn8Mc9PAakn3Exc/6wSwutsAJL+DngecLSkPcDv4bzb5oGP++SBN5rZhT1CYUwlHa7/RZKeiEte8C0gDS2yqa6/6NrN7K+K9t1s1w4dP/vnSToD16X2TeC/wua8/kCLEF4nEAgEAutO6HYLBAKBwLoTxCcQCAQC604Qn0AgEAisO0F8AoFAILDuBPEJBAKBwLoTxCcQCAQC604QnxlA0sExH/8vfQRiJP3WAPVPLkox0GP/JUm3dNj+Fklv6rcd04ikXZL+W+b1qT4twVg/80BgWIL4BIbGzH7JzO70L/sWnwH5uk+/MDamJHz/LqApPmY29vclEBgFQXxmFEmPlXS9jyR8vaSTfPn7JL1T0n9IukfSy3x5JOndku6Q9BFJ12a2fVLSWZIuBRb8k/eVeYvGJ0t7i19/hqRbJd0AvCGzTyyXWO/zvm3/teT1XOITjl0HPDFTfqqkf5N0s6R/l/SkTPmN/jxvTS0FSc+T9AlJVwFf8mU/J+lz/rr+PBUlSS+QdIOkL0j6oKRtvvxSSXf69v+fLm0+RtKHfBs+L+kHffnZ/v3/ov//RF9+eqYdt/mQRJcCqbXzjjLvVSAwEZhZWDb5AhwsKPtn4DV+/ReBf/Lr7wM+iHswOQ2425e/DLjWlx8L7AVe5rd9Ejgrfy7gZOD2zOs3AW/x67cBP+zX35HuB7we+G2/Pg/cBJySa3v+uM/ACcUisAOX9+VNftv1wG6/fg7w//n1jwCv9OsXpu3GhX45lJ4TeLJ/r6r+9buBVwNHA58Gtvry38RFpj4SuItW9JBdXT6Xq4Dn+PWTgC/79R1Axa8/H/iQX/9TXNBZgDlgIf9edPvMwxKWSVpCbLfZ5VnAT/v1vwH+d2bbP5mLMHynpEf7sucAH/Tl35P0iUFPLGkn7qb8qcz5z/PrLwCellpVwE5cQMlvdDnkDwH/aD4gq6Rr/P9twLOBD/p4eeAEDdz1v8SvXwVkLZTPmUteBvBjOHH7vD/GAvAA8EycOH/Wl88BNwD7gWXgLyX9C07kOvF84LRM23ZI2u6v+f3esjGg6rffAFwi6UTg/5rZ1zJ1A4GpIohPICUb5G8ls67c/36o0961uyVzrE5BBQX8dzP7aJ/nKjpeBDxi/Y+BHMq15/1mdnF2B0kvBj5uZq/MV5Z0Nk60XoHLS/SjHc4TAc8ys6Vc/T8FPmFmPyXpZJxliZldJek/gZ8APirpl4B7+ry2QGAiCGM+s8t/4G6O4DKHfqbH/p8BXurHfh6N654qoiYpfVK/H3iUpKMkzeNyFWFmjwD7JD0nc/6UjwK/nB5D0hMkbe3Rtk8DPyVpwVsOL/bn2Q98Q9LL/bEk6Qd8nRtp5Qx6Rf6AGa4HXibpUf4YR0p6rK//g5Ie78sXfVu3ATvN7Frg13AJ0jrxMVpJ85CL7AzO8rnXr/9CZvvjgHvM7J24dANPAw4A27ucIxCYSIL4zAaLkvZkljcCvwK8VtJtwM8Dv9rjGB/C5Zq5Hfhz4D+BfQX7XQ7cJulKM6sBb/X7foRWjh6A1wLv8g4H2Sf/v8SF0P+Cd1b4c3pY6Gb2BeDvgVt8O/89s/lVwOsk3QrcAVzgy38NeKOkzwHHdbgWzHnx/TbwMf9efRw4zsy+jxOGv/PlNwJPwgnBR3zZp4Bf79L0XwHO8s4Dd+LTSOC6QP9A0meBrMfdzwC3y7mYPwm4wswewnX93R4cDgLTREipECiNpG1mdlDSUbikXj9oLunZerfjZOAjZvaUIY6xCCyZmUl6Bc754IJe9aYFSQfNbNtGtyMQ6EQY8wn0w0ck7cINrr9tI4TH0wB2SrplgPGclGcAfyY3Yv8IzuNv6pF0Ks76u3+j2xIIdCNYPoHAmJF0CfDyXPEHzeztG9GeQGASCOITCAQCgXUnOBwEAoFAYN0J4hMIBAKBdSeITyAQCATWnSA+gUAgEFh3/n97ZLh1lhx5sgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"rain_test.sel(time='1995-01-01T15:00:00').plot()"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"<xarray.DataArray 'Rainf' ()>\n",
"array(0.)\n",
"<xarray.DataArray 'Rainf' ()>\n",
"array(0.)\n"
]
}
],
"source": [
"# Check data is the same:\n",
"diff_test = rain_test - rain\n",
"print(diff_test.max())\n",
"print(diff_test.min())"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"# Remove _FillValue for netcdf file\n",
"enc = {'lat':{'_FillValue':False}, \n",
" 'lon':{'_FillValue':False}, \n",
" 'Rainf':{'_FillValue':False}}\n",
"\n",
"rain_test.to_netcdf('test.nc',encoding=enc)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python [conda env:analysis3-19.07] *",
"language": "python",
"name": "conda-env-analysis3-19.07-py"
},
"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.6.7"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
@ccarouge
Copy link
Author

ccarouge commented Dec 2, 2019

Do not try the whole file at once. The 30min dataset is about 80GB... You might do better working on months and concatenate afterwards? You might want to set the unit of the time coordinate to a common start.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment