{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Time interpolation" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import mikeio" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "dims: (time:3, y:31, x:31)\n", "time: 2004-01-01 00:00:00 - 2004-01-03 00:00:00 (3 records)\n", "geometry: Grid2D (ny=31, nx=31)\n", "items:\n", " 0: Sign. Wave Height (meter)\n", " 1: Peak Wave Period (second)\n", " 2: Mean Wave Direction (degree)" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds = mikeio.read(\"../tests/testdata/waves.dfs2\")\n", "ds" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Interpolate to specific timestep\n", "\n", "A common use case is to interpolate to a shorter timestep, in this case 1h." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "dims: (time:49, y:31, x:31)\n", "time: 2004-01-01 00:00:00 - 2004-01-03 00:00:00 (49 records)\n", "geometry: Grid2D (ny=31, nx=31)\n", "items:\n", " 0: Sign. Wave Height (meter)\n", " 1: Peak Wave Period (second)\n", " 2: Mean Wave Direction (degree)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds_h = ds.interp_time(3600)\n", "ds_h" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And to store the interpolated data in a new file." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "ds_h.to_dfs(\"waves_3h.dfs2\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Interpolate to time axis of another dataset" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Read some non-equidistant data typically found in observed data." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "dims: (time:24)\n", "time: 2004-01-01 01:00:00 - 2004-01-03 12:00:10 (24 non-equidistant records)\n", "items:\n", " 0: Sign. Wave Height (undefined)\n", " 1: Peak Wave Period (undefined)\n", " 2: Mean Wave Direction (undefined)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ts = mikeio.read(\"../tests/testdata/waves.dfs0\")\n", "ts" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The observed timeseries is longer than the modelled data. Default is to fill values with NaN." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "dsi = ds.interp_time(ts)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "DatetimeIndex(['2004-01-01 01:00:00', '2004-01-01 02:00:00',\n", " '2004-01-01 03:00:00', '2004-01-01 04:00:00',\n", " '2004-01-01 05:00:00', '2004-01-01 06:00:00',\n", " '2004-01-01 07:00:00', '2004-01-01 08:00:00',\n", " '2004-01-01 23:00:00', '2004-01-02 00:00:00',\n", " '2004-01-02 01:00:00', '2004-01-02 02:00:00',\n", " '2004-01-02 03:00:00', '2004-01-02 04:00:00',\n", " '2004-01-02 05:00:00', '2004-01-02 06:00:00',\n", " '2004-01-02 07:00:00', '2004-01-02 08:00:00',\n", " '2004-01-02 09:00:00', '2004-01-02 20:00:00',\n", " '2004-01-02 21:00:00', '2004-01-02 23:00:00',\n", " '2004-01-03 00:00:00', '2004-01-03 12:00:10'],\n", " dtype='datetime64[ns]', freq=None)" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dsi.time" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(24, 31, 31)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dsi[\"Sign. Wave Height\"].shape" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAGjCAYAAADEoqqHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABz30lEQVR4nO3deVhUZRvH8e+wuwBurAqIu6apgOKO5q6VmpZbi5aW2aZmpWVvZZmVLWq55p6mZpppueFGLuS+JoorIIKIC+DCOs/7x5ERBJVRmGHg/lzXXDLPnBl+o0fm5pzn3I9OKaUQQgghhChGrMwdQAghhBDC1KQAEkIIIUSxIwWQEEIIIYodKYCEEEIIUexIASSEEEKIYkcKICGEEEIUO1IACSGEEKLYkQJICCGEEMWOjbkDFEZ6vZ4LFy7g6OiITqczdxwhhBBC5IFSiqSkJDw9PbGyuv8xHimAcnHhwgW8vLzMHUMIIYQQDyEqKopKlSrddxspgHLh6OgIaH+BTk5OZk4jhBBCiLxITEzEy8vL8Dl+P1IA5SLztJeTk5MUQEIIIYSFycv0FZkELYQQQohiRwogIYQQQhQ7UgAJIYQQotiRAkgIIYQQxY4UQEIIIYQodqQAEkIIIUSxIwWQEEIIIYodKYCEEEIIUexIASSEyLukWNgyXvtTCCEsmBRAQog7HlTgJMVCyFdSAAkhLJ4UQEKIO3IrcJSCW1ch7jhc/M982YQQIh/JWmBCCI1SEL5e+3r9h5CeDNfjtGJIn5Z92x2TIXAI2NiBo7t2E0IIC2L2I0BTp07F19cXBwcH/P392bZt2323DwkJwd/fHwcHB6pUqcL06dOzPZ6WlsbYsWOpWrUqDg4O1K9fn3Xr1hXkWxDCsiXFQuQuWNQLtn6pjUXsgOh9kBCVs/gB+G85zGkPM4Ng9893XkfmBwkhLIRZC6ClS5cybNgwPvroIw4cOEDLli3p3LkzkZGRuW5/9uxZunTpQsuWLTlw4AAffvghb7/9NsuXLzdsM2bMGGbMmMGPP/7IsWPHGDJkCD169ODAgQOmeltCWJadP8KcDnBqY+6PN3wRXg6GV0Og6/famEPZO4/v/0Urgq5FyfwgIYTF0CmllLm+eWBgIH5+fkybNs0wVrt2bbp378748eNzbP/BBx+watUqwsLCDGNDhgzh0KFDhIaGAuDp6clHH33EG2+8Ydime/fulC5dmoULF+YpV2JiIs7OziQkJODk5PSwb0+Iwu/CAVjUG25cBHtH8HsJQn+CpyaDR31tm6ynuC4c1I76vLIRYg/Btu8hMVp7rJQL3LgEgzZBpQCzvB0hRPFmzOe32Y4Apaamsm/fPjp06JBtvEOHDuzcuTPX54SGhubYvmPHjuzdu5e0NO0wfUpKCg4ODtm2KVGiBNu3b8/H9EIUAf/9AXM6a8VPhZraEZ56z2qPedQHzwbaLbf5Pda20GgQDFwHzYdByQpa8QOwpD/snaMVS3I0SAhRSJmtAIqPjycjIwM3N7ds425ubsTG5v5DMzY2Ntft09PTiY+PB7SC6Pvvv+fkyZPo9XqCg4P5888/iYmJuWeWlJQUEhMTs92EKLL0etj6FSwbAOm3oFp7GBQM5as++LmO7hA06k5RdHAR7JgIN+PvbHM9Fv4arh0p2j7xzrjMERJCFCJmnwSt0+my3VdK5Rh70PZZxydNmkT16tWpVasWdnZ2vPnmmwwcOBBra+t7vub48eNxdnY23Ly8vB727QhRuKXehN8Hwtbbp5ibvgn9loKDs3b/7gLnbo7u0Gb0nccDBmpHjl4N0U6bAVRuCdz+f7p3LoR8A2nJ0kNICFGomK0AqlChAtbW1jmO9sTFxeU4ypPJ3d091+1tbGwoX748AC4uLqxcuZIbN24QERHB8ePHKV26NL6+vvfMMnr0aBISEgy3qKioR3x3QhRCCdEwtxMcWwlWtvD0T9BxHFhl+eXg7gLnQRzd75wqy5wz1OELGLINfJpDRjJsGQdTGsM5OQ0thCg8zFYA2dnZ4e/vT3BwcLbx4OBgmjVrlutzmjZtmmP7DRs2EBAQgK2tbbZxBwcHKlasSHp6OsuXL6dbt273zGJvb4+Tk1O2mxBFyvm98HMbiDkEJcvDS6vB74WC+37u9WDA39DlO21+0LUI2PCR9tj+BXB+X845QnKKTAhhQmY9BTZixAhmzZrFnDlzCAsLY/jw4URGRjJkyBBAOzLz4osvGrYfMmQIERERjBgxgrCwMObMmcPs2bMZOXKkYZtdu3axYsUKzpw5w7Zt2+jUqRN6vZ7333/f5O9PiELh8G8wtwtcvwiudWDwFvBpmv/f5+7TZzqdNjE66/wggL2zYdYT2hyh0Kl3xuUUmRDChMzaCbp3795cvnyZsWPHEhMTQ926dVmzZg0+Pj4AxMTEZOsJ5Ovry5o1axg+fDhTpkzB09OTyZMn07NnT8M2ycnJjBkzhjNnzlC6dGm6dOnCL7/8QpkyZUz99oQwL70eNn8O22/37qnZBZ6ZqV3uXhAyT59lFTAQanbWvj6zFTZ+ArYlIe2mNrZrBiRfhcavFUwmIYS4B7P2ASqspA+QsHgp12HFq3Dib+1+i+HwxP/AyowHfTN7CL28AS6fhH+nw8Ujdx4vX10bv1cPIiGEeABjPr9lLTAhipprkbC4L1w8CtZ28PSPUL+PuVPdYWMPDZ+HBv3hzze0S+lBK34AVr99Z9ugUTmPKgkhRD6QAkiIoiTyX60R4c14KOUKfX4Fr0bmTqXJbY5Q2/9B41ch4bx2eX5Gqrb0RqNX7jxHCCEKgJwCy4WcAhMW6cAiWP2Otnipez3ouwScK5k7Vd79/S7smaVdNTbsMNiVMnciIYSFsYilMIQQ+USfAes/gj+HasVP7afh5fWWVfwA1Out/XkzHnZMMm8WIUSRJwWQEJYsOREW99EWMAUI+gCenW+ZR0/KekOd7trXOyZpq8sLIUQBkQJICEt15QzMbg8nN4CNA/SaA20+NO+VXo/C0R2enad1kE5P1i6ZF0KIAmKhPymFKObOboOfn4BLx8HRAwauhbo9H/y8wk6ng07jAR0cXQ4RoeZOJIQooqQAEsLS7J0Lv3SHW1fB00/r7FzRz9yp8o9H/TvLdKwbpTV0FEKIfCYFkBCWIiMd1rwPfw0Dfbp2xGfgGnDyMHey/PfEx2DnCDEH4dCv5k4jhCiCpAASwhLcugqLesHuGdr9J8ZAz9lgW8K8uQpKaVcIek/7etNYSEkybx4hRJEjBZAQhV38KZjVDs5s0dbReu4XaPWeNl+mKAscAuWqaIu4bvpMVooXQuQrKYCEKMxOb9FWTr98Cpwqaf196jxt7lSmYWMPHcZpX++dKyvFCyHylRRAQhRGSsHun2FhT0hOgEqN4dUt4PG4uZOZVs3OUKW1NudJCCHykRRAQhQ2GWnw9whYMxJUBtTvCy+t1ubFFCdJsRBzCPwGALdP9+3+WVtV/sJBORokhHgkshiqEIXJzSvw24twbhugg/afQbO3i/58n9xknvbK6uBC7QayUrwQ4pFIASREYXHpBPzaG66eBbvS0HOWdgqouAoYeOf9XzgIf72jfW3vBN2nQaUAs0UTQlg+KYCEKAxOBsPvL0NKIpTxhr5Lwa2OuVOZl6O7dsuqQk2IPwEbP4VBwWaJJYQoGmQOkBDmpBTs/Al+fU4rfnyaa52di3vxcy+dxoNTRbh8En57SZsvJYQQD0EKICHMJT0FVr0JGz4CpYeGL8ALK6FUBXMnK3wc3bU5P26PQd8lYFsKzobcniiuzJ1OCGGBpAASwhxuxMOCbnBgIeisoNNX8PSPYGNn7mSFk6O7NuHZ0V1rBdBrNqCDffMgdIp2RZg0ShRCGEEKICFM7eJ/MLMNRIZqE3r7LYMmrxfPK70eVs3O0OEL7esNY+DIcmmUKIQwihRAQpjS8TUwuwMkRGrLPAzaCNXbmTuVZWr6BvgPABRsHmvuNEIICyMFkBCmoBRs+x6W9IPU6+DbCgZtApea5k5mua5f1OZNVfSH9GRtbMUg+Oc7iNx952iQnB4TQuRCLoMXoqClJcPqt+HwUu1+wCvQ+WuwtjVvLkuXW6PE+JN3jgZV9NfmVWWkadvV7JzzsnohRLElBZAQBSnponbUJ3ov6Ky1wqfxYHOnKhqyNkqMOaQVmTU6Q/Q+uBGn/TmtGbjU1ra5fjH310mK1YqpgIFSIAlRjEgBJERBiTkEi/tCYjQ4lIHn5msLe4r8kVujxNajwL0eHF0B++dDxE64FKY99utzUMpVu4rMpznU6KgVR0mxcoRIiGJICiAhCsKxP+GPIZB2E8pXh35LoXxVc6cqHqys4fKp2+up3eVGHJzaqN02fQYOzuCS2XRS+gkJUZxIASREflIKQr6BrV9q96u2hV5zoEQZs8Yq8jIbJWYewcnt9FjnCVpxFHsEYg9r7QiSEyAqVNtuz2xoNOjO68nRICGKNCmAhMgvqTfhzzfgvxXa/cDXtV411vLfrMBlNkrMev/uAsarMXg20L7eMl6bI5TVgV+0G8hK80IUA/KTWYj8kHhBm+8TcxCsbKDrd7d71IhCKdtK8wfgr2Ha12514anJ4FzRbNGEEKYhBZAQjyp6HyzuB9djoUQ56P0LVG5h7lQi092nxzLH7j5CZFMSLh6F05sg6H3TZhRCmJw0QhTiURz5HeZ20Yofl9oweLMUP4VN1nXE7qflcO3PrV9B1J6CzyWEMCspgIR4GHo9bP4Clr+idSGu0Qle2QDlfM2dTBgr8whRwxegbi9QGVpH6ZSkBz9XukwLYbHydApsxIgRRr/wmDFjKFeunNHPE6LQS7kOf7wGx//S7jd7G9p9ql1hJCxP1gnUXb+DqF1w9RysHQXdp9z/udJDSAiLlacCaOLEiTRt2hQ7O7s8vej27dt58803pQASRc+1KG2y88UjYG0HT02CBv3MnUrklxJloMcMmNcVDi7UFqr1bnrvTtFKegcJYanyPAn6jz/+wNXVNU/bOjo6PnQgIQqtyF2wtD/cuASlXKD3IvAONHcqkd8qN4eWI2Dbd7B6GDwzK/tRnisREL4WIrbDmRDtOf9Ogyava19LDyEhLEKeCqC5c+fi7Oyc5xedMWMGbm5uDx1KiELn4GKtmV5GKrjVg76/Qhlvc6cSBaX1aDi9WbtEfss4bezkBq0oOrEW9GnZtz+8RLuB9BASwkLolJJjuHdLTEzE2dmZhIQEnJyczB1HmJM+Q1syYcck7X6tJ7VTJPalzZtLFKykWO1KsBWDtEnud3MoB74twckTdk3Txqzt4amJUPUJOQIkhJkY8/ktV4EJcS8pSdpK7pnFT8uR8NwvUvwUB3vnwm/P5178AAQOht4LoH4f7X6lxpCRAhvGaJPkhRCFXp4KoLJly1KuXLk83Yw1depUfH19cXBwwN/fn23bclnAMIuQkBD8/f1xcHCgSpUqTJ8+Pcc2EydOpGbNmpQoUQIvLy+GDx9OcvI9fpAJkZur52B2Bwhfp/1m/8wsaPsxWMnvDMVCwEB4NUS7tXhXG3tq8p2xgJezb99+LHg0gJuXYeEzcD3O5JGFEMbJ81VgmS5fvswXX3xBx44dadq0KQChoaGsX7+ejz/+2KhvvnTpUoYNG8bUqVNp3rw5M2bMoHPnzhw7dgxv75zzK86ePUuXLl0YPHgwCxcuZMeOHQwdOhQXFxd69uwJwKJFixg1ahRz5syhWbNmhIeHM2DAAAB++OEHo/KJYurcDlj6PNy6AqXdoc+vUMnf3KmEKd09kXn7d+BR/85aYlm3Cxql9X/qvwxmt9eK50XPwoC/5WihEIWY0XOAevbsSZs2bXjzzTezjf/0009s3LiRlStX5vm1AgMD8fPzY9q0aYax2rVr0717d8aPH59j+w8++IBVq1YRFhZmGBsyZAiHDh0iNFRb0fnNN98kLCyMTZs2GbZ599132b179wOPLmWSOUDF2L758PcI0Kdrv9H3XazN8xDF14WDMDNIO/JzdwF0t8untSLo5mWo1g76LgFrW1OkFEJQwHOA1q9fT6dOnXKMd+zYkY0bN+b5dVJTU9m3bx8dOnTINt6hQwd27tyZ63NCQ0NzbN+xY0f27t1LWpp2VUaLFi3Yt28fu3fvBuDMmTOsWbOGrl273jNLSkoKiYmJ2W6imMlIh3WjtSu99OnwWA8YuFaKH5H7WmL3Ur4q9PsNbErAqY2w+h3pFSREIWV0AVS+fHn++OOPHOMrV66kfPnyeX6d+Ph4MjIyclwu7+bmRmxs7m3lY2Njc90+PT2d+Ph4APr06cPnn39OixYtsLW1pWrVqrRp04ZRo0bdM8v48eNxdnY23Ly8vPL8PkQRcOsa/Poc/DtVu9/6Q+g1F+xKmjWWKCTyupZYpkoB8Ow80FnBwUWw5csCjSeEeDhGrwb/2Wef8corr7B161bDHKB///2XdevWMWvWLKMD6HS6bPeVUjnGHrR91vGtW7cybtw4pk6dSmBgIKdOneKdd97Bw8PjnnOURo8enW25j8TERCmCiovLp+HX3nD5pPZbe4/p8Fh3c6cSlq5mJ+j6Pfw1DP75Bpw8ck6cFkKYldEF0IABA6hduzaTJ09mxYoVKKWoU6cOO3bsIDAw711xK1SogLW1dY6jPXFxcfdsouju7p7r9jY2NoajTx9//DEvvPACgwYNAqBevXrcuHGDV199lY8++girXK7isbe3x97ePs/ZRRFxZiv89hIkXwOnitpk5wfN8RAirwIGQlIMhHwNf7+rTaiv1cXcqYQQtxldAIE2eXnRokWP9I3t7Ozw9/cnODiYHj16GMaDg4Pp1q1brs9p2rQpq1evzja2YcMGAgICsLXVJhrevHkzR5FjbW2NUgrp+SgMdv8Maz/QVv6uGAB9FknzOpH/Wo+GxGg4sBB+fxleWg1ejcydSgjBQzZCPH36NGPGjKFfv37ExWn9LtatW8d///1n1OuMGDGCWbNmMWfOHMLCwhg+fDiRkZEMGTIE0E5Nvfjii4bthwwZQkREBCNGjCAsLIw5c+Ywe/ZsRo4cadjmqaeeYtq0aSxZsoSzZ88SHBzMxx9/zNNPP421tazWXexlpGm/ja8ZqRU/9Z7TLleW4kcUBJ0OnpwI1dpD+i1trln8KXOnEkIAKCNt3bpVlShRQrVr107Z2dmp06dPK6WU+vrrr1XPnj2NfTk1ZcoU5ePjo+zs7JSfn58KCQkxPPbSSy+poKCgHN+/YcOGys7OTlWuXFlNmzYt2+NpaWnq008/VVWrVlUODg7Ky8tLDR06VF29ejXPmRISEhSgEhISjH4/ohC7cVmpeU8q9YmTUp84K/XPd0rp9eZOJYqD5CSlZgRp+94P9ZRKjDV3IiGKJGM+v43uA9S0aVOeffZZRowYgaOjI4cOHaJKlSrs2bOH7t27Ex0dXTCVmglJH6Ai6FI4LO4NV86AbSno+TPUundrBCHy3fVLMLud1ijRo/7tRomO5k4lRJFSoH2Ajhw5km3OTiYXFxcuX75s7MsJUfBOboRZ7bTix9kbXtkgxY8wvdIu8PwKKFkeYg5pE/Az0h78PCFEgTC6ACpTpgwxMTE5xg8cOEDFihXzJZQQ+UIpCJ0Kvz4LKQng1QQGbwb3uuZOJoqr8lWh3zKwLQmnN8Gqt6VRohBmYnQB1K9fPz744ANiY2PR6XTo9Xp27NjByJEjs01YFsKs0lO1rs7rR4PSQ4Pn4aVV2m/hQphTJX+t0abOCg79ClvGmTuREMWS0QXQuHHj8Pb2pmLFily/fp06derQqlUrmjVrxpgxYwoioxDGuXEZfukO+xdoHzIdxkG3n8BGej2JQqJmJ3jy9uLM/0yAPbPNm0eIYsjoSdCZTp8+zYEDB9Dr9TRs2JDq1avndzazkUnQFuziMW2y87VIsHOEXnOgRocHP08Ic9gyHkK+0gr13gtlbpoQj8iYz++HLoCKMimALNSJdbD8FUi9DmUrQ9+l4FrL3KmEuDelYNVbcOAXbSmWl1aBV2NzpxLCYhnz+W10J+iMjAzmzZvHpk2biIuLQ6/XZ3t88+bNxr6kEI9GKdg5GYI/ARRUbgnPLYCS5cydTIj70+m0U2HXL8LJDdq6dK9sgApF54i6EIWV0QXQO++8w7x58+jatSt169a978KlQhS49BRY/Q4cWqzd9x8AnSeAjZ1ZYwmRZ9a22urx856EC/th4TPwykZwzH1NRCFE/jD6FFiFChVYsGABXboU3UX95BSYhbgeB0v6w/ndoLOGTl9B48Hab9VCWJrrl2B2e7h6VholCvGQCrQRop2dHdWqVXvocELki5jDMLONVvzYO8Pzv0Pgq1L8CMtV2gWeXw4lK9xulPiiNEoUogAZXQC9++67TJo0SVZWF+ZzbBXM6QiJ56F8NRi8Cao+Ye5UQjy68lWh32+3GyVu1iZIy89aIQqE0XOAtm/fzpYtW1i7di2PPfYYtra22R5fsWJFvoUTIhulYNu3sPkL7X6V1trciRJlzZlKiPxVyV/brxf31ea2OVWEth+bO5UQRY7RBVCZMmVyXQtMiAKVdgv+fBOO/q7db/wadPwSrI3ehYUo/Gp01K4OW/22VvQ7eUCjQeZOJUSRYvSnx9y5cwsihxD3lhgDS/ppV8hY2UCXCRDwsrlTCVGw/F+CpBjYOh7WvAeOHtIoUYh8ZPQcICFMKno//NxGK35KlIUX/pDiRxQfQR+A34vaena/vwxRu82dSIgiI09HgPz8/Ni0aRNly5alYcOG9+39s3///nwLJ4q5o8th5VBIT4YKNaHfEihXxdyphDAdnQ66/gBJsbcbJT4HrwRLo0Qh8kGeCqBu3bphb68tJNm9e/eCzCME6PXa+kghX2v3q3eAnrPAwdm8uYQwB2sbaZQoRAHIUwFUtmxZrKy0s2UDBw6kUqVKhvtC5KvUG/DHEAhbpd1v+ia0HwtW1ubNJYQ52ZXSLo/PbJS4qBcMXCONEoV4BHmqYkaMGEFiYiIAvr6+xMfHF2goUUwlnIc5nbTix8oWuk2BjuOk+BECsjdKjD2sNUpMTzV3KiEsVp6OAHl6erJ8+XK6dOmCUorz58+TnJyc67be3t75GlAUE1F7tCu9bsRByfLQexH4NDV3KiEKl/JVof9v2umw05u1y+S7T5MO6EI8hDytBTZz5kzeeust0tPT77mNUgqdTkdGRka+BjQHWQvMxA4tgVVvQ0YKuD4GfRdDWR9zpxKi8ArfAIv7gMqAlu9C2/+ZO5EQhYIxn995Xgw1KSmJiIgIHn/8cTZu3Ej58uVz3a5+/frGJy5kpAAyEb0eNn0GOyZq92t2gWdmyrwGIfJi/wJtqQyALt9qCwELUcwZ8/md50aIjo6O1K1bl7lz59K8eXPDVWFCPJSUJFjxKpxYo91vMRye+B/I5Hoh8sbvRa1J6NYv7zRKrP2kuVMJYTGM/rR56aWXuHXrFrNmzWL06NFcuXIF0Pr/REdH53tAUQRdjYDZHbXix9oeesyEdp9K8SOEsYLeB7+XAAXLX4HIXeZOJITFMHopjMOHD9OuXTucnZ05d+4cgwcPply5cvzxxx9ERESwYMGCgsgpioqIUFjaH25ehlKu0OdX8Gpk7lRCWCadDrp+f7tR4npY3Bte3gAuNcydTIhCz+hfuYcPH86AAQM4efIkDg4OhvHOnTvzzz//5Gs4UcTs/wXmP6UVP+6Pw6tbpPgR4lFZ28Czc6GiP9y6Cgt7agWREOK+jC6A9u7dy2uvvZZjvGLFisTGyn86kQt9Bqz/CFa9Cfo0qP00vLwOnCuZO5kQRUNmo8RyVSAhUmuUmJxo7lRCFGpGF0AODg6GpohZnThxAhcXl3wJJYqQ5AT4tTeE/qTdD/oAnp2v/cAWQuSfUhWyNEo8Ar+9II0ShbgPowugbt26MXbsWNLS0gDQ6XRERkYyatQoevbsme8BhQW7fBpmtYdTwWDjAL3mQJsPZbKzEAWlXBXovwxsS8GZrdpl8nnrdCJEsWP0J9G3337LpUuXcHV15datWwQFBVGtWjUcHR0ZN25cQWQUlujsPzCrLcSf0C7PHbgW6kqBLESBq+gHz80HnTUcXgKbxpo7kRCFUp4bId5t8+bN7N+/H71ej5+fH+3atcvvbGYjjRAf0d45Wl8SfTp4+mlXejl5mDuVEMXL/l+0eXcgjRJFsVEgnaCLk2JZACXFwt65EDAQHN0f7rkNX4Cdk2D3TG28bi/o9hPYlsj/vEKIBwv5BraMA3TQ+xeo/ZS5EwlRoPK9E/TkyZPz/M3ffvvtPG8rCpGkWAj5Cmp2frgCKOQrba5P9D5t7Ikx0HKkLNIohDm1eg8SzsP++bB8ELz4J3g3MXcqIQqFPBVAP/zwQ7b7ly5d4ubNm5QpUwaAa9euUbJkSVxdXaUAslQ3rzz8c+NPan9G7wPbktp6XvKbphDml9ko8fpFCF+nXZH5ygZwqWnuZEKYXZ4mQZ89e9ZwGzduHA0aNCAsLIwrV65w5coVwsLC8PPz4/PPPy/ovOJhJcXClvHZG6QlxcJ/K2FhL1j4jDb25xuwbwFcOJBz26zPT4qFg4thbhdYMUgbcygLT00GZy9pxCZEYWFto12BWTEAkq9p/9/l/6cQxs8Bqlq1Kr///jsNGzbMNr5v3z569erF2bNn8zWgORTJOUAXDsLMIHg1BDwbQMxh7ZB4/Il7P6dOd+0Hp5V1ludvhdQbsHIoXIu493ODRkGb0fn6FoQQj+BGPMxuD1fOgHs9GLAGHIrIzzchbiuQ1eAzxcTEGHoAZZWRkcHFixeNfTlhahf/gy1fausGZarcErwD4Z9vwacZRO3ROjYfWwk/HoKmb2g/MAH+fAsuHrn9RCuo3gEq+WsTLZ+aDB71tYeMnUckhChYmY0SZ3e40yix3zKwsTN3MiHMwugjQE899RSRkZHMnj0bf39/dDode/fuZfDgwXh5ebFq1aqCymoyReYIUFLsnUPdx/+Gf77J8qAOanbRJiu71cl+hMjJU7uSa/dMrZMzaD1FVIb2tZWN9tzmb0OlRjmPLgkhCq/o/TDvSUi7AY/3hh4z5GIFUWQY8/ltdCPEOXPmULFiRRo3boyDgwP29vYEBgbi4eHBrFmzjA47depUfH19cXBwwN/fn23btt13+5CQEPz9/XFwcKBKlSpMnz492+OtW7dGp9PluHXt2tXobBZv71ytMJkZdFfxA6C0ozpudXI+r7SrVhgFvJJl84w7X+vTIWwVnNxYILGFEAUoW6PEpbDpM3MnEsIsjD4F5uLiwpo1awgPD+f48eMopahduzY1atQw+psvXbqUYcOGMXXqVJo3b86MGTPo3Lkzx44dw9vbO8f2Z8+epUuXLgwePJiFCxeyY8cOhg4diouLi2EZjhUrVpCaemf9m8uXL1O/fn2effZZo/NZvBod4cQ6iD14Z6zdp1CljfZ11tNUju7avJ2sY4GvQZ1u2mKmJ9bAtm9zP82V23OFEIVX9fbw9GTtooftP4CjJwS+au5UQpiUWRshBgYG4ufnx7Rp0wxjtWvXpnv37owfPz7H9h988AGrVq0iLCzMMDZkyBAOHTpEaGhort9j4sSJ/O9//yMmJoZSpfK2AKfFngLL2szw3Hb4awSkJGjrAjV9E/75+uFPU8lpLiGKnpAJsOULQAfPLYA6T5s7kRCPpEAnQWdkZDBv3jw2bdpEXFwcer0+2+ObN2/O0+ukpqayb98+Ro0alW28Q4cO7Ny5M9fnhIaG0qFDh2xjHTt2ZPbs2aSlpWFra5vjObNnz6ZPnz55Ln4sWmZDwuh9WlNC0C59fWYmpCRpBZAQQmRqNRISz8O+edpVoaX+BJ+m5k4lhEkYXQC98847zJs3j65du1K3bl10Dzl5Lj4+noyMDNzc3LKNu7m5ERube4+K2NjYXLdPT08nPj4eD4/s603t3r2bo0ePMnv27PtmSUlJISUlxXA/MTHRmLdSeMQc1v48FQw6K60LbKv3wNpWK44e5TSVnOYSoujR6aDLd5B0EcLXwuI+0ihRFBtGF0BLlizht99+o0uXLvkS4O4CSil136Iqt+1zGwft6E/dunVp3LjxfTOMHz+ezz6z0ImAmVd6ndsOG8ZoYyXKQfvPwb0u3LysFS2O7o/Wl+dRny+EKJwyGyXOfwqi98LCnvBKsCxgLIo8o68Cs7Ozo1q1ao/8jStUqIC1tXWOoz1xcXE5jvJkcnd3z3V7Gxsbypcvn2385s2bLFmyhEGDBj0wy+jRo0lISDDcoqKijHw3ZpR5pdeGj4Db07luXYFVb2jje+eaNZ4QwgLYlYR+v0G5qpAQBYuehWQLPRIuRB4ZXQC9++67TJo0iUedO21nZ4e/vz/BwcHZxoODg2nWrFmuz2natGmO7Tds2EBAQECO+T+//fYbKSkpPP/88w/MYm9vj5OTU7abxajRQTviA+BSS/vzqcnaZOVXQ7QJ0UII8SClymuNEku5aM1Olz4P6akPfp4QFsroU2Dbt29ny5YtrF27lsceeyxH4bFixYo8v9aIESN44YUXCAgIoGnTpsycOZPIyEiGDBkCaEdmoqOjWbBgAaBd8fXTTz8xYsQIBg8eTGhoKLNnz2bx4sU5Xnv27Nl07949x5GhIiU5AVa+oR3xcX0MukyAeV20y9TlSi0hhLHK+UL/ZTC3K5wN0S6T7zEDrIz+XVmIQs/oAqhMmTL06NEjX7557969uXz5MmPHjiUmJoa6deuyZs0afHx8AG3ZjcjISMP2vr6+rFmzhuHDhzNlyhQ8PT2ZPHmyoQdQpvDwcLZv386GDRvyJWehlJEGv70El8KgtDv0/01b60cIIR6FZ0Ptkvhfn4Mjv2md4dtb6BxJIe7DrH2ACqtC3wdIKVj9NuxfALYlYeBa7YhP1j5AcrWWEOJRHFgEfw7Vvu78jdYYVYhCrkD7AAkzyixw9Gla8aOz0q7eyDzdJVdqCSHyS8P+kHQBNn8Baz/Qfr7U6WbuVELkG6MLIF9f3/tepn7mzJlHClSsPegITmajw0ydvoKanU2XTwhRvLQcCQnRsG8uLB8MpVylUaIoMowugIYNG5btflpaGgcOHGDdunW89957+ZWreMoscGp2zlkApVyHiCwdsgOHyCFpIUTB0umgy7dw/aK2HuDiPvDyenCtZe5kQjyyh+oEnZspU6awd+/eRw4k0Do6XzgA8eHa1/HhcCPuzuNudeHx3tr6XJlNDoUQoiBY20DP2bDgaTi/Bxb1kkaJokjIt0nQZ86coUGDBpa7jEQWJp0EndnJGeDfqXB4qXHPDxol836EEAXvxmWY3R6unNZ+CRu4BhyczZ1KiGzMMgn6999/p1y5cvn1csXH3rnZ5/XcrUprqNsLyvjA1bPa1V9PTdZ6/YAc/RFCmEZmo8TZHeDiUa1RYv/lYGNn7mRCPBSjC6CGDRtmmwStlCI2NpZLly4xderUfA1XLAQMhJqdYO1oiArVxu4ucDKLnAu3q1lpdCiEMIfMRonzusLZf7TL5HvMlEaJwiIZXQB17949230rKytcXFxo3bo1tWrJxDijObpD1G6t+NHZgEqXAkcIUXh5NoDn5sOvveHIstuNEseaO5UQRjO6APrkk08KIkfxlZKk9dgAaNAXDvxy720d3bU5P3LaSwhhTtXawdM/wsrXYcckcPSEJkPMnUoIo+SpAEpMTDRqMnBSUhKOjo4PHapY2fKl1mysbGWt54ZTxXsXONLoUAhRWDToB4kXYPPnsG6UdlWYNEoUFiRPJ27Lli1LXFzcgze8rWLFitIQMS9iDsGu6drXXb6DcpW1AkeO8AghLEHLdyHgZUBpjRKz9ioTopDL0xEgpRSzZs2idOnSeXrRtLS0RwpVLOgzYPUwUHp4rAdUb2fuREIIYZzMRolJF+HE37cbJW6QRonCIuSpAPL29ubnn3/O84u6u7tja2v70KGKhb1z4MJ+sHeCjuPNnUYIIR6OlTX0nAULusH53bCwJwwK1iZHC1GI5akAOnfuXAHHKEaSYmHnj7Bvnnb/iY+lo6oQwrLZlYR+S7VGiZdPwaJnpVGiKPSkeYOpJcVC6E+Qeh08G0KjV8ydSAghHl3JclqjxFKuWqPEJf0hPcXcqYS4JymATC1q9+0vdPDkRO3wsRBCFAVlK2uNEu1Kw7ltsHIo6PXmTiVErqQAMoWkWG3h0shdEPK1NubbSvvzwsE7a4EJIYSl82wAzy0AKxs4+jtslN5xonCSAsgU9s6FmUEwpwPcjNfGzoZoYzODtMeFEKKoqNYWnv5J+3rnZPh3mnnzCJELoztBR0ZG4uXllW09MNAulY+KisLb2zvfwhUZAQOhZmdAwZ45cGCBLGgqhCjaGvTVmrxuGgvrRoOjBzzW3dyphDAw+giQr68vly5dyjF+5coVfH198yVUkePorh0WzjrpOXO9L88GUgAJIYqmFiMg4BVAwYpX4dwOcycSwsDoAkgplePoD8D169dxcHDIl1BCCCGKAJ0OukyAWk9CRgos6QtxYeZOJQRgxCmwESNGAKDT6fj4448pWbKk4bGMjAx27dpFgwYN8j1gkSMLmgohipOsjRKjdt1ulLhRGiUKs8tzAXTgwAFAOwJ05MgR7OzsDI/Z2dlRv359Ro4cmf8JixpZ0FQIUdzYloC+S2B2B7h8Ehb2gpfXSqNEYVY6pZQy5gkDBw5k0qRJRq0Ob2kSExNxdnYmISGhSL9PIYQwqasRWrfo6xehckutcaKNvblTiSLEmM9vo+cAzZ07V4oCIYQQxivrc1ejxNelUaIwG6Mvg79x4wZfffUVmzZtIi4uDv1dO++ZM2fyLZwQQogixqM+9P5FWy/s6HJtLlCHL8ydShRDRhdAgwYNIiQkhBdeeAEPD49crwgTQggh7qnqE9BtCvzxmrY4tKMnNB364OclxWqNYwMGyoUk4pEZXQCtXbuWv//+m+bNmxdEHiGEEMVB/T6QeAE2fQbrPwQnD/Buev8CJykWQr7SGstKASQekdFzgMqWLUu5cuUKIosQQojipMVwaDQIQ6PE8HVagZPb+oipNyHmkMkjiqLL6CNAn3/+Of/73/+YP39+tl5AQgghhFF0Ouj8jVbwHP8L1n1057HEGIj6F05vgajdEB8OKkN77L8Vd7ZzdJejQeKh5Oky+IYNG2ab63Pq1CmUUlSuXBlbW9ts2+7fvz//U5qYXAYvhBAmkhSrXR7/9wi4eFQbK1EObl3J2/ODRklvNWFgzOd3no4Ade/ePT9yCSGEENntnaud9soqa/HjVg/qdAePulqhtPZ9bbxyS+jwubbIqhAPwehGiMWBHAESQggTSYq9M+fn7DYIHgNN3oDaT4JtyeynuC4chJlBoLMBlQ5dv7+zwLQQFMARICGEEKJA5DaH5/HnwLPBvZ8T+Cr8O1W7esy7KbjVKdCIomh66KvA7r6VL1+eihUrEhQUxNy5cwsiqxBCiOIsczHppm9BtXaQngy/vwxpt8ydTFggowug//3vf1hZWdG1a1c+++wzPv30U7p27YqVlRVvvPEGNWrU4PXXX+fnn38uiLxCCCGKqswC515XdWUuJu3sCd2nQSlXuBQG6z/KfXsh7sPoOUA9e/akffv2DBkyJNv4jBkz2LBhA8uXL+fHH39k5syZHDlyJF/DmorMARJCCAtwejP80kP7uvdCqP2UefMIsyvQxVDXr19Pu3btcoy3bduW9evXA9ClSxdZE0wIIUTBqvoENH9H+/rPNyHhvHnzCItidAFUrlw5Vq9enWN89erVhg7RN27cwNHR8dHTCSGEEPfTZgx4+kHyNVg+GPQZ5k4kLITRBdDHH3/Me++9x9NPP80XX3zBuHHj6NatG++//z6ffPIJAMHBwQQFBeXp9aZOnYqvry8ODg74+/uzbdu2+24fEhKCv78/Dg4OVKlShenTp+fY5tq1a7zxxht4eHjg4OBA7dq1WbNmjbFvVQghRGFnYwe9ZoOdI0TuhH8mmDuRsBAP1Qdox44d/PTTT5w4cQKlFLVq1eKtt96iWbNmRr3O0qVLeeGFF5g6dSrNmzdnxowZzJo1i2PHjuHt7Z1j+7Nnz1K3bl0GDx7Ma6+9xo4dOxg6dCiLFy+mZ8+eAKSmptK8eXNcXV358MMPqVSpElFRUTg6OlK/fv085ZI5QEIIYWEO/wYrBoPOCgb8DT7GfR6JosGYz2+zNkIMDAzEz8+PadOmGcZq165N9+7dGT9+fI7tP/jgA1atWkVYWJhhbMiQIRw6dIjQ0FAApk+fzoQJEzh+/HiOZTrySgogIYSwQH8MgUOLwakSDNkGJWXh7uIm3ydBJyYmZvv6fre8Sk1NZd++fXTo0CHbeIcOHdi5c2euzwkNDc2xfceOHdm7dy9paWkArFq1iqZNm/LGG2/g5uZG3bp1+fLLL8nIuPd54ZSUlId+H0IIIQqJLhOgXBVIPA+r3wZZ6EDcR54KoLJlyxIXFwdAmTJlKFu2bI5b5nhexcfHk5GRgZubW7ZxNzc3YmNjc31ObGxsrtunp6cTHx8PwJkzZ/j999/JyMhgzZo1jBkzhu+++45x48bdM8v48eNxdnY23Ly8vPL8PoQQQhQS9o7Qaw5Y2ULYatgnTXnFveVpKYzNmzcbrvDasmVLvgbIuso8gFIqx9iDts86rtfrcXV1ZebMmVhbW+Pv78+FCxeYMGEC//vf/3J9zdGjRzNixAjD/cTERCmChBDCEnk2hHafwoaPYN1obakM19rmTiUKoTwVQFmv6Mrr1V0PUqFCBaytrXMc7YmLi8txlCeTu7t7rtvb2NhQvnx5ADw8PLC1tcXa2tqwTe3atYmNjSU1NRU7O7scr2tvb4+9vf2jviUhhBCFQZOhcGYLnNqoLZUxeDPYljB3KlHIGH0ZPMC2bdt4/vnnadasGdHR0QD88ssvbN++Pc+vYWdnh7+/P8HBwdnGg4OD73k1WdOmTXNsv2HDBgICAgwTnps3b86pU6fQ6/WGbcLDw/Hw8Mi1+BFCCFHEWFndWSoj7hhsGGPuRKIQMroAWr58OR07dqREiRLs37+flJQUAJKSkvjyyy+Neq0RI0Ywa9Ys5syZQ1hYGMOHDycyMtKwzMbo0aN58cUXDdsPGTKEiIgIRowYQVhYGHPmzGH27NmMHDnSsM3rr7/O5cuXeeeddwgPD+fvv//myy+/5I033jD2rQohhLBUpV2hx+0+cXtmwf4FsGU8JOU+x1QUQ8pIDRo0UPPnz1dKKVW6dGl1+vRppZRSBw4cUG5ubsa+nJoyZYry8fFRdnZ2ys/PT4WEhBgee+mll1RQUFC27bdu3aoaNmyo7OzsVOXKldW0adNyvObOnTtVYGCgsre3V1WqVFHjxo1T6enpec6UkJCgAJWQkGD0+xFCCFGIrB+j1CdOSo3z1P6MPmDuRKIAGfP5bXQfoJIlS3Ls2DEqV66Mo6Mjhw4dokqVKpw5c4Y6deqQnJxcMJWaCUkfICGEKCLSU2FOB7hwQLs/aDNU8jdvJlFgCnQxVA8PD06dOpVjfPv27VSpUsXYlxNCCCEKRlKsNgeo5Uiwvn2hy5J+sGuGVhDJ6bBizegC6LXXXuOdd95h165d6HQ6Lly4wKJFixg5ciRDhw4tiIxCCCGE8fbOhZlBsLQ/ZGjzVbkeC2vfh5mtYetXd7ZNipU5QqZUCP6+83QZfFbvv/8+CQkJtGnThuTkZFq1aoW9vT0jR47kzTffLIiMQgghhPECBkLNztrXMYe07tBVn4Cz/4A+XWuUqE/TVpS/fhFCvtK2d3Q3b+7iICnW7H/fRhdAAOPGjeOjjz7i2LFj6PV66tSpQ+nSpfM7mxBCCPHwHN1zfri2/QQcnGHTZ/DfH3BgIRxdAfWe0x5PSYJb10wetdhJSTJ3gocrgECbDB0QEJCfWYQQQoiCV84Xnp0Hj/eGTWO1eUL752mPzX/SnMmKn5hDd77OrWAtQHkugJ555pk8bbdixYqHDiOEEEIUCEd3CBqV/QP2wkGt+BHms/rtO18HjYI2o032rfNcADk7O2e7/+uvv/LUU0/h6OiY76GEEEKIfOXonvPDNbc5Ql2/B/fHtbHSbjIfKD8lxWpzrQBiD8PfI+CpyeBRXxsz8d91nguguXOzr6r7+++/880338il70IIISxTbqdcKvqDZwOzxCnyynprNwBrbfkqPOqb7e/7odYCE0IIIYSwZFIACSGEELnNERIFpxD8fT/0VWBCCCFEkZHbHCFRcArB33eeC6BVq1Zlu6/X69m0aRNHjx7NNv7000/nTzIhhBBCiAKS58VQrawefLZMp9ORkZHxyKHMTRZDFUIIISyPMZ/feT4CpNfrHzmYEEIIIURhIJOghRBCCFHsSAEkhBBCiGJHCiAhhBBCFDtSAAkhhBCi2JECSAghhBDFzkMVQNeuXWPWrFmMHj2aK1euALB//36io6PzNZwQQgghREEwuhP04cOHadeuHc7Ozpw7d47BgwdTrlw5/vjjDyIiIliwYEFB5BRCCCGEyDdGHwEaMWIEAwYM4OTJkzg4OBjGO3fuzD///JOv4YQQQgghCoLRBdCePXt47bXXcoxXrFiR2NjYfAklhBBCCFGQjC6AHBwcSExMzDF+4sQJXFxc8iWUEEIIIURBMroA6tatG2PHjiUtLQ3Q1v+KjIxk1KhR9OzZM98DCiGEEELkN6MLoG+//ZZLly7h6urKrVu3CAoKolq1ajg6OjJu3LiCyCiEEEIIka+MvgrMycmJ7du3s3nzZvbv349er8fPz4927doVRD4hhBBCiHynU0opY55w7tw5KleuXEBxCofExEScnZ1JSEjAycnJ3HGEEEIIkQfGfH4bfQqsSpUqtGjRghkzZhiaIAohhBBCWBKjC6C9e/fStGlTvvjiCzw9PenWrRvLli0jJSWlIPIJIYQQQuQ7owsgPz8/JkyYQGRkJGvXrsXV1ZXXXnsNV1dXXn755YLIKIQQQgiRr4yeA5Sb/fv388orr3D48GEyMjLyI5dZyRwgIYQQwvIU6BygTFFRUXzzzTc0aNCARo0aUapUKX766aeHfTkhhBBCCJMx+jL4mTNnsmjRInbs2EHNmjXp378/K1euLPJXhgkhhBCi6DC6APr888/p06cPkyZNokGDBgUQSQghhBCiYBldAEVGRqLT6QoiixBCCCGESRhdAGUWPzdv3iQyMpLU1NRsjz/++OP5k0wIIYQQooAYXQBdunSJAQMGsG7dulwfLwpXgQkhhBCiaDP6KrBhw4Zx7do1/v33X0qUKMG6deuYP38+1atXZ9WqVUYHmDp1Kr6+vjg4OODv78+2bdvuu31ISAj+/v44ODhQpUoVpk+fnu3xefPmodPpctySk5ONziaEEEKIosnoI0CbN2/mzz//pFGjRlhZWeHj40P79u1xcnJi/PjxdO3aNc+vtXTpUoYNG8bUqVNp3rw5M2bMoHPnzhw7dgxvb+8c2589e5YuXbowePBgFi5cyI4dOxg6dCguLi707NnTsJ2TkxMnTpzI9lwHBwdj36oQQgghiiijjwDduHEDV1dXAMqVK8elS5cAqFevHvv37zfqtb7//nteeeUVBg0aRO3atZk4cSJeXl5MmzYt1+2nT5+Ot7c3EydOpHbt2gwaNIiXX36Zb7/9Ntt2Op0Od3f3bDchhBBCiExGF0A1a9Y0HF1p0KABM2bMIDo6munTp+Ph4ZHn10lNTWXfvn106NAh23iHDh3YuXNnrs8JDQ3NsX3Hjh3Zu3cvaWlphrHr16/j4+NDpUqVePLJJzlw4ECecwkh7i0uMZkfgsOJS5RTykIIy/ZQc4BiYmIA+OSTT1i3bh3e3t5MnjyZL7/8Ms+vEx8fT0ZGBm5ubtnG3dzciI2NzfU5sbGxuW6fnp5OfHw8ALVq1WLevHmsWrWKxYsX4+DgQPPmzTl58uQ9s6SkpJCYmJjtJoTIKS4phUmbThKXJIsfCyEsW57nAJ06dYpq1arRv39/w1jDhg05d+4cx48fx9vbmwoVKhgd4O6eQkqp+/YZym37rONNmjShSZMmhsebN2+On58fP/74I5MnT871NcePH89nn31mdHYhipOwmET+OHDe3DGEECJf5PkIUI0aNfDy8uLFF19k3rx5nDt3DoCSJUvi5+dndPFToUIFrK2tcxztiYuLy3GUJ5O7u3uu29vY2FC+fPlcn2NlZUWjRo3uewRo9OjRJCQkGG5RUVFGvRchioq7T3Elp2Uwb+dZOk/8h86TtjF7+zkA3v3tIMv3n+dodIJhWzk9JoSwJHk+AhQSEkJISAhbt27ljTfeIDk5GW9vb5544gnatGlDmzZtqFixYp6/sZ2dHf7+/gQHB9OjRw/DeHBwMN26dcv1OU2bNmX16tXZxjZs2EBAQAC2tra5PkcpxcGDB6lXr949s9jb22Nvb5/n7EIUVZmnuB7zdGJvxFWW7Y3i6s20HNuduHidd387BECfRl581fNxw3Pb13HD1UmuuhRCFG55LoBatmxJy5YtGTNmDGlpaYSGhrJ161a2bt3K4sWLSUlJoVq1ajkuP7+fESNG8MILLxAQEEDTpk2ZOXMmkZGRDBkyBNCOzERHR7NgwQIAhgwZwk8//cSIESMYPHgwoaGhzJ49m8WLFxte87PPPqNJkyZUr16dxMREJk+ezMGDB5kyZUqecwlRHKVn6Nl56jIAr/6yzzDu5mhP+zpudHjMjQvXkhm14gh+3mU5EHkVBSzZE8WlpBSefDzvF0EIIYS5Gd0HCMDW1pZWrVrRqFEjmjZtyvr16/n55585deqUUa/Tu3dvLl++zNixY4mJiaFu3bqsWbMGHx8fAGJiYoiMjDRs7+vry5o1axg+fDhTpkzB09OTyZMnZ+sBdO3aNV599VViY2NxdnamYcOG/PPPPzRu3Phh3qoQRd7R6AQW/RvB+mOxXLlx52hPgE8ZOtf1oEs9DzzKlDBsCzC222OUsrfhm/XHWXcklk3H49h0PA6AdUfvnKZ2dbSXo0FCiEJJpzJnEedBcnIyO3fuZMuWLWzdupU9e/bg6+tLUFAQrVq1IigoyKjTYIVVYmIizs7OJCQk4OTkZO44QuQ7vV6x/VQ8i3ZFsOG/i9zvh8A7baszvH0NQCuAnvxxO3+91YK6FZ35ITicSZvuPb8u63PjEpNZtCuS/oHeUhQJIQqEMZ/feT4CFBQUxJ49e6hatSqtWrXirbfeIigo6J4TloUQhc+VG6ks2xvFr7sjibh80zBe19OJzvU8cC5hy5iVR/nqmXrUregMaEdxMrk62vNO2+qGsf6B3rSvo/0M2HI8ju+Cw9GBoaDaH3GVo9EJ1K3oLHOEhBCFSp4LoJ07d+Lh4UGbNm1o3bo1rVq1eqjL3oUQpqWUYl/EVRbtiuTvIzGkpusBcHSwoadfJfoHelPdzRG4c4qrbkVnQwGUlauTg+GITub9rMXMd8HhTH3ejw3/XeTPg9FsOxXPth+3076OG13ryRwhIUThkecC6Nq1a2zbto2tW7fy9ddf07dvX2rUqEFQUBCtW7cmKCgIFxeXgswqhDBCUnIaKw9eYNG/ERyPTTKM16vozPNNvHmqvicl7R5qGuB9eZUtyQ+9G/BGm2p8u/446/67SPAx7Qaw4T+ZIySEMD+j5gBllZSUxPbt2w3zgQ4dOkT16tU5evRofmc0OZkDJCzZsQuJLNwVwZ8HormRmgGAg60VT9f35PkmPjxeqcw9n/so83Rye67MERJCmFKBzAG6W6lSpShXrhzlypWjbNmy2NjYEBYW9rAvJ4R4BMlpGfx9OIaFuyI4EHnNMF7NtTT9A715xq8SziVy75WV1d2nuIyR23OzzhHafPwi3wdnL4YORV0jLCaR2h5OMkdICGFSeS6A9Ho9e/fuZevWrWzZsoUdO3Zw48YNKlasSJs2bZgyZQpt2rQpyKxCiLucjb/Bon8j+H3/ea7dblhoa62j42PuPN/Eh0DfcvddWqag3T1H6Pvgk0zr35C1Ry+y+vAFtoZfYmv4JTrXdaeLzBESQphQngugMmXKcOPGDTw8PGjdujXff/89bdq0oWrVqgWZTwhxl7QMPRuPXWThrgh23G5cCFCxTAn6BXrzXIAXLo6Ft7O5V7lSTO7bkLeeqMY3608QfOwia4/GsvZ2/6DgsIuGbWWOkBCioOS5AJowYQJt2rShRo2HOzwuhHg0F67dYsnuSJbsiTKsxq7TwRM1XXm+iQ+tarhgbWW+oz0Pcvcl9NXdHKnj4WSYHJ1p0saTTNqonSrLnCMk84OEEPntoSdBF2UyCVoUFnq94p+Tl1i0K5JNYRfR3/7fWqG0PX0aedGnsReVypY0b8hHEJeYbCjmNoZdZOLG7HOE2tZy5f1OtUjL0GdrwCiEELkxySRoIUTBuXw9hd/2nufX3RFEXbllGG9apTz9m3jToY47djZWZkyYP+6eIzRx40l+6teQtUdi+ftIDJuOx7H5RBzNq0rPMSFE/pICSIhCQinFnnNXWbQrgrVHYknN0BoWOjnY0Mvfi36B3lRzLW3mlAWvcvlSTOnvR7/T8Uzdcoodpy6z/VQ8AB+uOMKQoCp4ly8l84OEEI9ECiAhzCwxOY2VB6JZ9G8kJy7eaVhY36sM/QO9eepxT0rYWZsxoWncPUdo95kr2SZ5AxyOTmDorwcAeKGJN593rwdIDyEhhPFkDlAuZA6QMIWj0Qks2hXBnwcvcPN2w8ISttZ0a+BJ/0Af6lUq3nNdss4POhqdwKgVR6jr6cTRC4mGbZ6u78nbbauRnCZzhIQQZpoDFBMTQ1paGt7e3vn1kkIUOclpGaw+dIGFuyI5FHXNMF7dtTTPN/Ghh19FnBwe3LCwOLh7fhDAVz0fx0qnY/Kmk6z7L5ZVhy6w+vAFWlWXOUJCCOPkWwH0xBNPEB4eTkZGRn69pBBFxulL11n0byS/74siMTkd0BoWdq7rwfNNfGhUuaxZGxZakjqeTkx/wZ9tJy8xdcspQs9cISRcmyM0ZuURXguqilfZkjJHSAhxX/lWAC1YsICbN2/m18sJYfHSMvRs+O8iC/+NIPTMnbksXuVK0K+xD88GVKJC6cLbsLAwuXt+EMDec1cJPXMl23YHoxJ4feF+AF5s6sPYbnVlfpAQIlcyBygXMgdIPIroa7dYvCuSpXujuHR7DouVDp6o5Ub/Jt4EVXfBqhA3LLQUuc0Rqu3hSFiMNpFcB3RvWJGOj7kxZOF+mR8kRDFgkjlAqampxMXFodfrs43LHCBRHGXoFf+EX2LRrgg2H48zNCx0cbSnbyMvejf2pmKZEuYNWcTkNkdoQq/6gNZPaGPYRf44EM3KA9EAnL96UwogIYSB0QXQyZMnefnll9m5c2e2caUUOp1O5gCJYiX+egpL90SxeHck56/eaVjYvFp5+gf60L6OG7bWlt+w0JLUrejMlz3q0rWeB4t3R7D73FUAXl+4n9Y1XejT2JuGXmUMxZOcIhOieDK6ABowYAA2Njb89ddfeHh4yMRNUewopdh19gqLdkWy7mgMaRna4R7nErb08q9Ev0BvqroU/YaFhcndc4QW7Ypk0qbsy2ooYMuJS2w5cYla7o5Me94f3wqliEtKYdKmk7Sv4yYFkBDFiNEF0MGDB9m3bx+1atUqiDxCFFoJt9L4Y/95Fu2K5GTcdcN4A68yPN/Ehycf98DBtug3LCyMXJ0cGN7+zkLN/QO9aV/HDbgzP+jNNlXZdfYKe85d5XhsEu2+D6F7A22OkBCi+DG6AKpTpw7x8fEFkUWIQunw+Wss+jeSVYcucCtNO8Vb0s6abg0q0j/QW+aVFEK5zQ/qVNeDkR1rcSjqGt+sP86OU5dZvv88Kw6cByAkPO7O8+USeiGKPKMLoK+//pr333+fL7/8knr16mFrm71pm1w1JYqCW6mZDQsjOHw+wTBe082R55t4071hRRylYaFFqu9VhgCfcoZlNjKvg52wPpwJ68MBeKdtdcMRJZkjJETRZHQB1K5dOwDatm2bbVwmQYui4FRcEgv/jWT5/vMk3W5YaGdtRZd67vRv4kOAjzQstDS59RDKeops7dEYpmw5bXjMSqc1roy8fBPv8iVljpAQRZTRBdCWLVsKIocQZpOarmf9f7Es2hXBv1ka63mXK0n/QG96+VeivDQstFh3zw/KHMtazEzZcppvn32c1YdiCAm/xF+HY1h7NJaefhXpUMfd1JGFECZgdAEUFBRUEDmEMLmoKzdZvDuS3/ZGEX89FdB++29X243+TXxoWa2CNCwsRmq5O9HL34v9kVf5et1xdp25wm97z/P7Pm2O0D/hlwzbyhwhISyf0QXQ4cOHcx3X6XQ4ODjg7e2Nvb38tiwKpwy9YuuJOBbtimTLiTjD/A9XR3v6NPamb2MvPJylYWFxcvcpMj/vsjTxLc+u20cDM5tafrP+BN+sPwHcmSMk84OEsFxGF0ANGjS47xwIW1tbevfuzYwZM3BwkB8IonCIS0pm2d7z/LorkuhrdxoWtqhWgeebeNO2tjQsLK5yO0WWdY7Q34cvMC3kjOExax1EXL5B1JWbJNxKk/lBQlgoowugP/74gw8++ID33nuPxo0bo5Riz549fPfdd3zyySekp6czatQoxowZw7ffflsQmYXIE6UUoWcus2hXJOuPxpJ++1f5MiVteda/Ev0CffCtUMrMKUVhdPccoWkhZ/imZz1WH45h28l4Vh68wF+HY2hbW3oICWGpjC6Axo0bx6RJk+jYsaNh7PHHH6dSpUp8/PHH7N69m1KlSvHuu+9KASTMIuFmGsv3n2fRrghOX7phGPfz1hoWdqknDQuF8ep4OvNcI282/BfL1C2nOXj+Guv/iwVg7OpjvBZUBTcnB5kfJISFMLoAOnLkCD4+PjnGfXx8OHLkCKCdJouJiXn0dELkkVKKw+cTWPhvBKsPXyA5TVukt5SdNd0bVqR/oA91PKVHlTDe3XOE/ruQyMHz17Jts/vcFXaf0+YMvdy8Mv976jFTxxRCGMnoAqhWrVp89dVXzJw5Ezs7OwDS0tL46quvDMtjREdH4+Ymh4ZFwbuZms6qg1rDwqPRiYbxWu6OPN/Eh+4NK1La3ujdXAiDvCyzUdWllOFo44LQc6Sk6xnaphoVy8iEeiEKK6M/GaZMmcLTTz9NpUqVePzxx9HpdBw+fJiMjAz++usvAM6cOcPQoUPzPawQmcIvJrHo3whW7I8mKeV2w0IbK56s50H/Jt74eUvDQlEwcltmY1KfhtxISWfSppPsPK3NO/ttbxS9G3kxtHU1PKUQEqLQ0SmVeSFw3l2/fp2FCxcSHh6OUopatWrRr18/HB0dCyKjySUmJuLs7ExCQoIs7VGIpKRnsO5oLIt2RbL77J2GhT7lMxsWelGulJ0ZE4ri5mh0Ak/+uJ2/3mphWBPu3zOXmbTxJKFntKU27KyttEKoTVVpsSBEATPm8/uhCqCiTgqgwiXqyk0W7Ypk2d4oLt/QGhZaW+loV9uV55v40LyqNCwU5nG/PkChpy8zaVO4obu4nbUVfRp78XprKYSEKCj5XgCtWrWKzp07Y2try6pVq+677dNPP21c2kJICiDzy9ArNh+PY9GuCELCLxkaFro7OdCnsRd9Gnnj7ixX2ojCL/T0ZX7YGG44amlnbUXfxl683rqa7MNC5LN8L4CsrKyIjY3F1dUVK6t7N4srKouhSgFkPnGJySzdE8Xi3ZFcSEg2jLesXoHnm/jQtpYrNtKwUFignafjmRh80nC1mBRCQuQ/OQX2iKQAMi2l1O2JoxFs+O+ioWFh2ZK2PBfgRb9Ab3zKS8NCYfmUUoYjQnvOXQW0yfv9GnvzeuuquEn/ICEeiRRAj0gKINO4djOV3/dpy1Ocib/TsDDApyzPN/GhU113aVgoiiQphIQoGMZ8fuf5XMKuXbtYu3ZttrEFCxbg6+uLq6srr776KikpKUaHnTp1Kr6+vjg4OODv78+2bdvuu31ISAj+/v44ODhQpUoVpk+ffs9tlyxZgk6no3v37kbnEgVDKcWByKu8+9shAr/cxBd/h3Em/gal7W14oYkP64a15PfXm9G9YUUpfkSRpdPpaFatAr+91pRFgwIJ8ClLarqeeTvP0eqbLXy2+j/iEpMf/EJCiIeW5z5An376Ka1bt6Zz586A1hH6lVdeYcCAAdSuXZsJEybg6enJp59+mudvvnTpUoYNG8bUqVNp3rw5M2bMoHPnzhw7dgxvb+8c2589e5YuXbowePBgFi5cyI4dOxg6dCguLi707Nkz27YRERGMHDmSli1b5jmPKDg3UtL58+AFFv4bwbGYOw0L63g48XwTH55u4CkNC0Wxo9PpaF6tAs2qlmfHKe2I0L6Iq8zdcY5fd0XSP9CHIUFVZGkNIQpAnk+BeXh4sHr1agICAgD46KOPCAkJYfv27QAsW7aMTz75hGPHjuX5mwcGBuLn58e0adMMY7Vr16Z79+6MHz8+x/YffPABq1atIiwszDA2ZMgQDh06RGhoqGEsIyODoKAgBg4cyLZt27h27RorV67Mcy45BWac+10KfCI2iYX/RvDHgWiuZ21Y+LgHzzfxoaFXGWlYKMRtSim2n4rnh+Bw9kdeA8Dexornm/jwWlAVXB2lEBLifoz5/M7zr9xXr17NtrxFSEgInTp1Mtxv1KgRUVFReQ6ZmprKvn37GDVqVLbxDh06sHPnzlyfExoaSocOHbKNdezYkdmzZ5OWloatrS0AY8eOxcXFhVdeeeWBp9QAUlJSsp2+S0xMvM/W4m5xSSlM2nSS9nXccHVyICU9g7VHYln4bwR7I64atvOtUOp2w8JKlCkpDQuFuJtOp6NldRdaVKvAtpPx/LAxnAOR15i9/SyLdkXwfKAPr0ohJES+yHMB5ObmxtmzZ/Hy8iI1NZX9+/fz2WefGR5PSkoyFCB5ER8fT0ZGRo41w9zc3IiNjc31ObGxsblun56eTnx8PB4eHuzYsYPZs2dz8ODBPGcZP358tvciHk5Mwi1WH7rAsn3nuZKlYWGHOm4838SHZlXLy9EeIfJAp9PRqoYLLatX4J+T2hGhg1HXmLX9LAtvF0KvBVXF5fYCrUII4+W5AOrUqROjRo3i66+/ZuXKlZQsWTLb/JrDhw9TtWpVowPc/YGolLrvh2Ru22eOJyUl8fzzz/Pzzz9ToUKFPGcYPXo0I0aMMNxPTEzEy8srz88vjuISk4lLSiE9Q8/y/dEADF6wz/C4q6M9zzfxoXcjL7miRYiHpNPpCKrhQqvqFQgJv8TEjSezFUIvNNEKoQqlpRASwlh5LoC++OILnnnmGYKCgihdujTz5883rAYPMGfOnBynp+6nQoUKWFtb5zjaExcXd8+V5N3d3XPd3sbGhvLly/Pff/9x7tw5nnrqKcPjer0eABsbG06cOJFrkWZvb4+9vfwAuZ+75/nM+OcMs7efvef2vRt58Xbb6iZMKETRpdPpaF3TlaAaLoSEX+KHjSc5FHWNn7edZeG/kbzQ1IdXW1WRQkgII+S5AHJxcWHbtm0kJCRQunRprK2zX6K8bNkySpcunedvbGdnh7+/P8HBwfTo0cMwHhwcTLdu3XJ9TtOmTVm9enW2sQ0bNhAQEICtrS21atXiyJEj2R4fM2YMSUlJTJo0SY7qPILMeT4VStux49RlNhy7U4iWtLPmZmoG73esSasaLoB2BEgIkb+yFkJbwy8xMTicQ+cTmPnPGX4JjeDFpj4MlkJIiDwx+rpjZ2fnXMfLlStn9DcfMWIEL7zwAgEBATRt2pSZM2cSGRnJkCFDAO3UVHR0NAsWLAC0K75++uknRowYweDBgwkNDWX27NksXrwYAAcHB+rWrZvte5QpUwYgx7jIu6s3Ulmx/zwAH//5n2G8ceVy9G/ijVe5EjwzNZRWNVwMK2ILIQqOTqejTU1XWtdwYeuJS/ywMZzD5xOY8c8ZFoRG8GIzH15tWYXyUggJcU9mbbzSu3dvLl++zNixY4mJiaFu3bqsWbMGHx8fAGJiYoiMjDRs7+vry5o1axg+fDhTpkzB09OTyZMn5+gBJIx39ykupRSbwi6yaFck207Fk56hzbWyt7GifW1XOtXzoHHlcrg6OXA0OsHM6YUonnQ6HW1qudK6pgtbTsQxceNJrRAKyTwiVJlXW1WhXCm56lKIu8lSGLkojn2AjkYn8OSP2/nttaaEX9R69xyPTbrvc95pW53h7Wvctw+QEMJ0lFJsPq4VQkdu/2JS0s6al5pVZnBLKYRE0SdrgT2i4lgArT50gbcWH8DBxorkdG3iuJ21jpY1XOhSz4PUdD2jVxzhq2fqGU5zuTraS8EjRCGkHcGNY+KmcI5Ga33NpBASxUGBNEIUli23ozSRV26w8kA0a47EcDz2OgDJ6XoqlnGgcz0Pegd4Ud3NEcBwmqtuRWeZ5yNEIafT6WhXx422tV2zFULTtp5mwc5zhkKorBRCohiTAqiYyNqt+UZqBr/uimBBaAQpt4/2ZBV9LZlZ285Sys6G4e0dzZBWCJEfshZCG8PimLgxnP8uJDJ162nm7zzHgOaVGdRCCiFRPEkBVEykZ2iFzscrj3Ig6pph3KW0PZ3quuFVriRfrjme4xRXJldHe95pW10ubxfCAul0OtrXcaNdbVeCj11k4saTHItJZMqW08zfGcGAZpUZ1NJXlqgRxYrMAcqFpc4Buvs0V1xiMsdiEll/NJa/jsSQlJxu2DbApwx9GnnTw68S1lY6wyTov95qIae4hCjilFJsuF0IhcVoc4RK29swsHllXmkhhZCwXDIHqJjKPM3VtpYrx2OT+HTVf5yJv5HrtnsjrtG8mgvWVrI2lxDFjU6no+Nj7rSv7Xa7EArneGwSP24+xbwd524XQlVwLpn39R2FsDRSABUhCbfSAHh14T5iE5IN43U9nahbyZklu6PkFJcQwsDKSkenuu50qOPGhmOxTNx4kuOxSUzefIq5UgiJIk5OgeWisJ4Cy+1KrosJtwgJv8Tao7H8czKeDL32z1nSzpq2tVzpH+hNk6oV5BSXEOKB9HrF+v9imbTppKEPmKO9DQNb+PJKC1+cS0ghJAo3OQVWRGW9kquEnTUrD0QzceNJLt9IzbHtzdQMVh+OoYpLaZpUrWCGtEIIS2NlpaNzPQ86PubO+v+0I0InLiYxedNJ5u44y8vNfXlZCiFRREgBZIF+3HyKbScvcTM1AwA7Gx1B1V2o6e7IT1tO53qaS05xCSHyKmshtO6/WCbdLoQmbTrJnB1neaWFLwObSyEkLJucAsuFuU6B5XaKKy4xmairN9l2Mp7l+88TdeWWYftKZUvwXEAlXmqm/SCS01xCiIKg1yvWHo1l0qZwwi9qTVOdHGx4pUUVBraojJODFEKicJBTYBYq6ykuVycHzly6zru/HcrWtyer81dvkaFHfgsTQhQoKysdXR/3oHNdd9YcjWHSxpOcjLvODxvDmb39DINaVmFAcymEhGWRAqgQ2n4qnvFrw9hx6rJhzKW0PX4+ZVj/30W5kksIYRZWVjqefNyTznU9WHMkhkmbTnIq7jrfB4cze/tZBrXwZUDzyjhKISQsgJwCy0VBngLLrVlhXFIKl5JS+OXfc2w+finb9i2qleflFr4E1XAlLCZRTnEJIQqNDL3i7yMxTL5dCIF2RFoKIWEucgqsEMt6mqtCaXvGrQnjz4MX7rm9v085nqjlZsKEQgiRN9ZWOp6u70nXeh78dfgCkzed5PSlG3wXHM6s7WcZ3NKXl5pJISQKJymAzGTZ3vNsPnEx26TmKhVKcib+Jl90r0sDrzKAnOISQhR+1lY6ujWoyJOPe2YrhL7dkFkIVeGlZpUpbS8fOaLwkFNgucjvU2CZp7kiLt9g5j9nOHQ+wfBYKTtrnqzvweCWVUlOy5BTXEIIi5ehV/x1+AKTNp3kzCVtOZ4yJW2lEBIFTk6BFTKLdkUyadPJXB+7kZqBu1MJqrmW5mh0Qq7bCCGEJcl6RGj1Ie2I0Jn4G0xYf4JZ284wuFUVXmwqhZAwLzkClIuCOgKklOL74HC2nLiU40quzAnRd/cBEkIIS5eeoWf14QtM3nSKs7cXaC5b0pZXW1XlxaY+lLpdCMnPQPGojPn8lgIoFwV5FZg0KxRCFFfpGXpWHbrAj5vvFELlStnxaqsqvNDEh7PxN+Tno3gkcgpMCCFEoWNjbcUzfpV4ur4nfx68wI+bT3Lu8k2+Wnucmf+c4en6nuaOKIoROQKUC1P2ARJCiOIqPUPPgtAIZm47Q2xCsmG8c113Xmnhi4OttWGKgBB5IUeACjFXJweGt69h7hhCCGF2NtZWJNxKy1b8AKw9Gsvao7EADG1dlfc71QLkF0iRv6QAEkIIYTb9A71pX0dr9noo6hofrTxKuZJ2XLmZCsCS3ZGULWlH/ybeOdZLFOJRSAEkhBDCbFydHHIUM3MGBBAed52fNp8i8spNxq0JY8Y/p2WOkMhXVuYOIIQQQmRlY23FcwFebHo3iDFda+PmZE/89VTm7DgHwLStp9h77gpHoxOIS0y+/4sJcQ9yBEgIIUShcPdyP7bWViQlp3MxMSXbdn8fieXvI9ocoTdaV+W923OEhDCGXAWWi4K8CkwIIUTeZTaShTtzhMqWtOXqzTRA6yP0Rptq9A/0xsHW2pxRRSEgV4EJIYQoEnKbIzR3QCOOxybx4+ZTRF+7xed/HWN6yGleD6pKPymERB7JHCAhhBAWxcbaij6NvdkysjXjn6lHxTIluJSUwti/jtHqmy3M3XGW5LQMc8cUhZwUQEIIISzC3XOE7Gys6Hu7EPqyh1YIxSWl8NnqnIVQXGIyPwSHy6RpYSBzgHIhc4CEEMLypKbrWbYviimbT3HhdnNFNyd7Xg+qSr1KzvScFirrjBVxMgdICCFEsWNnY0X/QB96+Vdi2d7zTN2iFUKfrj5GuVK2AKSky6kxoZEjQLmQI0BCCGH5oq7eYPGuKH7bG0X8da2ztJODDf0CvenwmBuVypSUjtJFjDGf3zIHSAghRJH0+95opm49bSh+ABKT05kecoZnpobywfLDhiNCMkeo+JECSAghRJHUP9Cbv95qwV9vteCrZ+oB0K2BJ+VL2QGw5cQlWk/Yyi//RnD+2k0mbTpp6Dkkij6ZAySEEKJIyq2H0OCWVfi65+P8tjeKqVtOE5OQzMcrj1KhtFYUpWbIHKHiQgogIYQQxYqDrTUvNq1M65ouLPo3kmX77swRCjkRj5211kjR1dFe5ggVYWY/BTZ16lR8fX1xcHDA39+fbdu23Xf7kJAQ/P39cXBwoEqVKkyfPj3b4ytWrCAgIIAyZcpQqlQpGjRowC+//FKQb0EIIUQhd3cPIYDl+6KZ8c8ZrtxIM4xN2nSSJ3/czpM/bmfRrkhzRBUmYtYjQEuXLmXYsGFMnTqV5s2bM2PGDDp37syxY8fw9vbOsf3Zs2fp0qULgwcPZuHChezYsYOhQ4fi4uJCz549AShXrhwfffQRtWrVws7Ojr/++ouBAwfi6upKx44dTf0WhRBCFAKuTg4Mb18j21j/QG/a13ED4Gh0AqNWHOGrZ+oZ+gRlLZZE0WPWy+ADAwPx8/Nj2rRphrHatWvTvXt3xo8fn2P7Dz74gFWrVhEWFmYYGzJkCIcOHSI0NPSe38fPz4+uXbvy+eef5ymXXAYvhBDFy9HoBJ78cbs0SrRwFnEZfGpqKvv27aNDhw7Zxjt06MDOnTtzfU5oaGiO7Tt27MjevXtJS0vLsb1Sik2bNnHixAlatWp1zywpKSkkJiZmuwkhhBCi6DJbARQfH09GRgZubm7Zxt3c3IiNjc31ObGxsblun56eTnx8vGEsISGB0qVLY2dnR9euXfnxxx9p3779PbOMHz8eZ2dnw83Ly+sR3pkQQghLk9scIVG0mX0StE6ny3ZfKZVj7EHb3z3u6OjIwYMH2bNnD+PGjWPEiBFs3br1nq85evRoEhISDLeoqKiHeCdCCCEsVeYcIbnqq/gw2yToChUqYG1tneNoT1xcXI6jPJnc3d1z3d7Gxoby5csbxqysrKhWrRoADRo0ICwsjPHjx9O6detcX9fe3h57e6n6hRBCiOLCbEeA7Ozs8Pf3Jzg4ONt4cHAwzZo1y/U5TZs2zbH9hg0bCAgIwNbW9p7fSylFSop09xRCCCGExqyXwY8YMYIXXniBgIAAmjZtysyZM4mMjGTIkCGAdmoqOjqaBQsWANoVXz/99BMjRoxg8ODBhIaGMnv2bBYvXmx4zfHjxxMQEEDVqlVJTU1lzZo1LFiwINuVZkIIIYQo3sxaAPXu3ZvLly8zduxYYmJiqFu3LmvWrMHHxweAmJgYIiPvNKLy9fVlzZo1DB8+nClTpuDp6cnkyZMNPYAAbty4wdChQzl//jwlSpSgVq1aLFy4kN69e5v8/QkhhBCicDJrH6DCSvoACSGEEJbHIvoACSGEEEKYixRAQgghhCh2ZDX4XGSeFZSO0EIIIYTlyPzczsvsHimAcpGUlAQgHaGFEEIIC5SUlISz8/3XdJNJ0LnQ6/VcuHABR0fH+3alzk1iYiJeXl5ERUVZ1ARqyW1aktt0LDEzSG5Tk9ymVVC5lVIkJSXh6emJldX9Z/nIEaBcWFlZUalSpUd6DScnJ4vaGTNJbtOS3KZjiZlBcpua5Datgsj9oCM/mWQStBBCCCGKHSmAhBBCCFHsSAGUz+zt7fnkk08sbnFVyW1aktt0LDEzSG5Tk9ymVRhyyyRoIYQQQhQ7cgRICCGEEMWOFEBCCCGEKHakABJCCCFEsSMFkBBCCCGKHSmAiglLnOtuiZmFeVjivmKJmYXpWeJ+YimZpQAyQlpaGtHR0Yb7lvCPnJyczIcffsjMmTPNHcUoGRkZhjXZLElKSgo7d+4kIiLC3FGMYon7Nsj+bWqWuH/Lvm1alrRvSwGUR9999x3Vq1ena9euPPnkk4SGhhq9TpipzZ49G3d3d3bv3o2trS23bt0yd6Q8+eGHH2jevDndu3dn2LBhnD59GtDWaCvMJk6cSOXKlXnttdeoX78+06dPJyMjw9yxHsgS922Q/dvULHH/ln3btCxu31bigb777jtVuXJl9fvvv6s5c+aobt26qQoVKqgtW7aYO9o9nTp1SrVs2VJNnz7d3FHyLDw8XD3xxBOqevXqat68eWrUqFGqSZMmqn379uaO9kBjxoxRNWvWVH///bc6efKk+vjjj1WZMmXUzZs3zR3tvixx31ZK9m9Ts8T9W/Zt07HUfVsKoPvIyMhQaWlpqnPnzur111/P9liLFi1Up06d1IEDB8wT7gG+++471bBhQ6WUUhEREerjjz9Ws2bNUtu2bTNzstxlZGSo7777TnXp0kVFR0cbxpctW6bq16+vwsLCzJju/uLj41VgYKD69ttvDWPh4eGqTp066tKlS0oppfR6vbni5cqS922lZP82JUvbv2XfNi1L3rflFNh9WFlZodfrOXLkCH5+foB2Xha0w8Hh4eGsW7eOlJQUc8bMRt0+v33q1Cnat2/P2rVrCQgIYM+ePUyfPp22bdsybdq0QndIVSlF9erVeeONN/D09DQcMi1dujSxsbGUL1/ezAnvrVSpUhw+fDhbS/ePPvoIDw8Pli1bxoULF8yYLneWuG+D7N/mYGn7t+zbpmXJ+7YcAcpi/vz5atiwYWr+/Pnq4sWLhvFevXqpoKAgw/2MjAyllFJDhgxRAQEB6vz586aOanCvzKNHj1aVK1dWb775ppoxY4ZKT09XSin17rvvqsDAQLVmzRpzRVZKKbV27VrDb42Zf2b+vWYdW7p0qWrYsKFKSkoyfchcZM2tlDL8vX788ceqbNmyqlu3bsrJyUk1aNBADR8+XNWqVUsFBgaqVatWmSuyUsoy922lZP82NUvcv2XfNi1L3bdzIwWQUio2Nla1bdtWVaxYUfXo0UN5e3urmjVrqn///VcppdTChQtVxYoVDTverVu3lFJKRUVFKZ1OZ9jO3Jlr1aqldu7cqZRSat++fcrd3V1ZWVmpHTt2GJ6XkJCgqlevriZNmmTyzEopdezYMRUUFKR0Op36/PPPlVLZ//NkyvxP9M4776gXXnjhntuZSl5yR0REqMGDB6s+ffoYfmjp9Xrl5+enPvzwQ7Pkt8R9+165Zf8uOJa4f8u+bVqWum/fj5wCA7Zt20ZMTAz79+9nxYoVhIeH4+TkxOeff87Ro0dp3749jRo14ssvvwTAwcEBpRT29vZ4eXlx7NixQpHZ0dGRL7/8kgMHDlCvXj169OiBvb09VlbaP7Ner8fJyYly5cqZJfO5c+f4+uuvcXFx4e233+abb74hLi4OKyurHJemZl6psX79elq3bg1oh7bPnTsHmPZS1rzmLleuHHv27OGll17C2tqa5ORkdDod5cqV4+DBg4Z/B1OyxH37Xrll/zZv7sK2f8u+bTqWum8/SLEvgJRSbNu2DRcXFxwdHdHr9djb2zNx4kTi4uKYNWsWLi4uDB48mHPnzjFy5EjS0tLQ6XQcOXIEe3t7wz9yYch88eJF5s+fj5WVFW+88QYVK1bk888/5+TJk1hZWREWFoZer6dv374mzQzg6upKgwYNeO+99xg1ahS+vr4MHz78ntsfPnyYa9eu0alTJ65cucIrr7xClSpVCA8PN+mlrHnN7eDgQGJiIvv27TPcDw8P58aNG/Tu3dtkeTNZ4r79oNyyf5svd2Hav2XfNi1L3bcfyOTHnAqRzEN1H3zwgapdu7ZS6s45b6WU+vDDD1WTJk0Ms++XLFmiSpYsqfz8/NSAAQOUk5OTevXVV9WtW7dMdhVEXjNv3bpVKaXUzp07lZeXl/L29lbPPfecKl++vOrTp49KTEw0Sd67cycnJxvGfvvtt2yHou8+TLpkyRLVuHFjNX78eOXk5KRatWqljh07ZrrQyrjcqampasKECUqn06nevXurt99+W7m6uqqePXuqq1evmiW3Je3bxuSW/Tt/WOL+Lfu27Nv5RQogpZ1ztbW1VcHBwUqpO//QZ8+eVb6+vmrKlCmG52zatEl9++23auDAgeqvv/4qtJmnTp1q2Pbo0aNq0aJF6v3331dr1641eea7ZeZKTExUTz75pGrcuHGu27300ktKp9Op6tWrq+XLl5syYq7ymvv7779Xr776qurRo0eBT1i81w/wwr5vP2puc+3fefnALIz7d37mNtX+Lfu27NsFrcgXQJn/YFkr7bvFx8erZ555RtWrV88wlrl9165dVe/evQs25F0sMbNSect9t507dyp7e3u1YMECw3Mzf8OZM2eOmjlzZv4HvUt+5U5ISCiQfPeSmpqqUlNT77tNYdxPinLuuxWG/Tu/cpty/05OTjb0GFIq9w/nwriPFOXcdysM+/ajKtIF0A8//KBeeumlPG27adMmVbZsWfXFF18YxlJSUlTbtm3Ve++9V0AJc7LEzEoZlzur1NRUNWzYMOXl5aWOHTum+vXrpz766COVkpKS/yFzkZ+5P/74Y3X9+vX8D5mL8ePHq9atW6uuXbuqiRMnGn545Xa1RWHaT4pD7qzMvX/nZ25T7d9ffPGFqlevnmrWrJnq27evOn36tFIq919QCtM+UhxyZ2XufTs/FMkCKCwsTHXr1k2VKlVKubq6qsWLFyul7v8PmpKSoqZOnap0Op1699131ZYtW9SPP/6oPD091aZNmyRzPua+2969e5VOp1M6nU7VqlXLJOeKLTX3gQMHVOPGjVWdOnXUrFmzVN++fZWfn59666237vmcwrCfFKfcd7OU/eRups69c+dO5efnp2rXrq1+++03NWHCBNWsWTP17LPP3vM5hWEfKU6572aOfTs/FckCaN68eeqpp55SS5cuVQMHDlRNmzY1rFnzoEN7EyZMUC1atFC1a9dW3t7eatmyZaaIbJGZlXq03Hq9Xm3atEm5u7srb29vk86pssTcN27cUB988IHq27dvtkmnn3zyierSpYu6fPnyfZ9vrv2kuOa21P3EXLk///xz9eqrr2Y71TZ06FA1cODAbNlyY86fgcUxtzl/duenIlUAZR7WvX79ugoJCVFKKbVq1SrVsGFDw6HGex36vbuTpakqWUvMnPV7P0zuTGlpaWr8+PHq008/LdiwWVhqbqW0RmhTp05V//zzjyGHUkpNmjRJVatWzXD/bubcT5QqfrkzWdp+ksnUuTM/ZGNjY9XJkycN4+np6apdu3bqiy++UMePH8/1uebcR4pj7kzm2rfzm06pQtSV6CHMmzePqKgoGjduTIsWLShVqlS2x69evcr48eNZs2YNf/75J1WrViUjIwNra+tcX08pVeB9Ciwxc37nzsys1+sLvIFaUcjdvHlzSpcuTXp6OjY2NgCGDJ9//jkHDx7k999/v+9+YI79pLjmttT9xFS5H/R/csmSJQwZMgQPDw9Kly7N5cuXefrpp5k4ceJ9cxc0yW3afbvAmaXsygenT59W/v7+ytfXV7Vu3Vq5u7urDh06qPj4eMM2mVX2tm3bVKtWrdSgQYPMFVcpZZmZlZLcppZb7vbt22fLrdfrDb/JdevWzXD0ypyrcktu07LE3Hn5P6mUUjNmzFDLli1TKSkp6sqVK+qvv/5SOp1OHT16VCll+vyS23z7eUGy2AJo4sSJKjAwUN28eVPdvHlThYWFqfLly6tBgwapiIgIpdSdCa1paWlqwoQJqmbNmmrLli1KKaX+/fdfk/+jWmJmyV24ckdGRiql7hRuiYmJysXFJdvEyaioKKWU6X9oSW7J/SiZM/9P5pbpzJkzqnz58mr27Nkmy5qV5C6aLLIASk1NVW3btjX8tp75IbZs2TJVuXJlNW3aNMO2mT8A/vvvP9WrVy/VvHlz9eSTTyqdTqeOHDkimSW3xeZWSpu/5Ovrq/R6vbpw4YLq27evKl++vIqJiZHckrtQ5TY2c1bz589XgYGBKjY21iRZs5LcRZfFncBTSmFra4uzszMXLlwwjAH06tULPz8/VqxYwenTpwEM5yhdXFyIi4tj586d2NjYcObMGerWrSuZJbfF5gY4fvw49evX56uvvqJatWrExcWxb98+3N3dJbfkLjS5HyZzfHw8Fy5cYOzYsYwePZoePXpQoUIFky6mKblNm9vkTFZq5bNffvlFlS9f3nCOMrOd+L59+5S9vb3avn27YdsjR46omjVrqipVqhiu/DEHS8yslOQ2NWNyN2/e3NCDw9zLnEhu07LE3HnNfOzYMfXFF1+oGjVqqDp16qiNGzeaLbNSkruosrgCKLOl+759+1SbNm2ytQ3PPJVRt25d9b///c8wnpycrP744w+T5szKEjMrJblNzdjc169fV2PGjFHz5883fdgsJLdpWWLuvGb++OOPlVJK3bx5U23atEn9+eefpg+bheQu2gpdAbR8+fJcVxbOrUPvjBkzlIeHh1qyZIlhLCYmRtWsWVP9/PPP93xefrPEzEpJbkvPbSqSW3I/iCVmVkpymzp3YVNoCqCQkBBVp04dpdPpsq3ie7elS5cqW1tbtXz5cnXlyhU1fPhwVbJkSTV16lR16NAh9fXXX6tq1aqpw4cPS2bJLbklt+SWzJLbzLkLq0JRAJ04cUL169dPvf7662ro0KHK09NTnT9/3vC4Xq9XCQkJqk+fPsrFxUV9++232RZcGzJkiHrsscdU1apVlbe3t/r7778ls+SW3JJbcktmyW3m3IVZoSiA4uLi1KxZs9TRo0fVjRs3lLe3t3rnnXeybZOSkqKWLFmS7bK8zN4Fer1eXb9+Xe3du1cyP4Dkltx5Ibkld1HMrJTkNnXuwswsBdC2bduyNWFSSmVbn2bBggXKzs5O7du3Tyn14EZdD1q7KT9YYmalJHcmyX1/klsjue/NEjMrJbkzmSq3JTFpAbRp0ybl6+urfHx8lIeHh3rhhRfU/v37lVI5//GaN2+uunTpYpjNbi6WmFkpyW1qktu0JLfpWGJmpSS3eDCTFUBRUVGqadOm6qOPPlIRERFq9erVqkGDBqpt27aGFWmzzmDfuXOnsrKyUitWrFBKadXr3euXSGbJLbklt+SWzJLbvLktlckKoA0bNigHBwcVHh5uGFu/fn2OHgVZDRgwQD3++ONq48aNqlOnTmr06NGGRk6mYImZlZLckjtvJLfkfhBLzKyU5DZ1bktlsgJoyZIlqkGDBur48eOGsYyMDDVlyhRVuXJltX79eqVU9up2165dSqfTKZ1Op9q3b68uX75sqrgWm1lyS27JLbmLc2bJbfrclspkBdCRI0eUvb29WrlyZbbx48ePq+7du6sBAwYYxtLT09WiRYuUvb29CggIUHv27DFVzGwsMbNSktvUJLdpSW7TscTMSklukTcmWwy1bt26PPHEE/zwww9cv37dMF6zZk28vb2JjY3l1q1bAKSmpnL58mUmTpzInj17CAgIMFVMi88MktvUJLdpSW7J/CCSW+SJKautgwcPKhsbGzVt2rRs5yg///xz5e3tne0Sv8LCEjMrJblNTXKbluQ2HUvMrJTkFg9mY8piq379+nzwwQeMHTsWa2tr+vbti16vZ/fu3fTv3x8bG5PGyRNLzAyS29Qkt2lJbtOxxMwguUUemKPqGjp0qHJzc1ONGjVSPj4+qnbt2urIkSPmiJJnlphZKcltapLbtCS36VhiZqUkt7g3nVJKmbroSklJ4dixYxw8eBA7Ozv69+9v6ghGs8TMILlNTXKbluQ2HUvMDJJb3JtZCiAhhBBCCHMy2VVgQgghhBCFhRRAQgghhCh2pAASQgghRLEjBZAQQgghih0pgIQQQghR7EgBJIQQQohiRwogIYQQQhQ7UgAJIYQQotiRAkgIUWRs3boVnU7HtWvXzB1FCFHISSdoIYTFat26NQ0aNGDixIkApKamcuXKFdzc3NDpdOYNJ4Qo1GRZWSFEkWFnZ4e7u7u5YwghLICcAhNCWKQBAwYQEhLCpEmT0Ol06HQ65s2bl+0U2Lx58yhTpgx//fUXNWvWpGTJkvTq1YsbN24wf/58KleuTNmyZXnrrbfIyMgwvHZqairvv/8+FStWpFSpUgQGBrJ161bzvFEhRIGQI0BCCIs0adIkwsPDqVu3LmPHjgXgv//+y7HdzZs3mTx5MkuWLCEpKYlnnnmGZ555hjJlyrBmzRrOnDlDz549adGiBb179wZg4MCBnDt3jiVLluDp6ckff/xBp06dOHLkCNWrVzfp+xRCFAwpgIQQFsnZ2Rk7OztKlixpOO11/PjxHNulpaUxbdo0qlatCkCvXr345ZdfuHjxIqVLl6ZOnTq0adOGLVu20Lt3b06fPs3ixYs5f/48np6eAIwcOZJ169Yxd+5cvvzyS9O9SSFEgZECSAhRpJUsWdJQ/AC4ublRuXJlSpcunW0sLi4OgP3796OUokaNGtleJyUlhfLly5smtBCiwEkBJIQo0mxtbbPd1+l0uY7p9XoA9Ho91tbW7Nu3D2tr62zbZS2ahBCWTQogIYTFsrOzyzZ5OT80bNiQjIwM4uLiaNmyZb6+thCi8JCrwIQQFqty5crs2rWLc+fOER8fbziK8yhq1KhB//79efHFF1mxYgVnz55lz549fP3116xZsyYfUgshCgMpgIQQFmvkyJFYW1tTp04dXFxciIyMzJfXnTt3Li+++CLvvvsuNWvW5Omnn2bXrl14eXnly+sLIcxPOkELIYQQotiRI0BCCCGEKHakABJCCCFEsSMFkBBCCCGKHSmAhBBCCFHsSAEkhBBCiGJHCiAhhBBCFDtSAAkhhBCi2JECSAghhBDFjhRAQgghhCh2pAASQgghRLEjBZAQQgghih0pgIQQQghR7Pwf+43agEcWjnoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = dsi[\"Sign. Wave Height\"].sel(x=250, y=1200).plot(marker='+')\n", "ts[\"Sign. Wave Height\"].plot(ax=ax,marker='+')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Model validation\n", "\n", "A common metric for model validation is mean absolute error (MAE).\n", "\n", "In the example below we calculate this metric using the model data interpolated to the observed times.\n", "\n", "For a more elaborate model validation library which takes care of these things for you as well as calculating a number of relevant metrics, take a look at [fmskill](https://github.com/DHI/fmskill).\n", "\n", "Use `np.nanmean` to skip NaN." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "name: Sign. Wave Height\n", "dims: (time:24)\n", "time: 2004-01-01 01:00:00 - 2004-01-03 12:00:10 (24 non-equidistant records)\n", "values: [0.06521, 0.06771, ..., 0.0576]" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ts[\"Sign. Wave Height\"]" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "name: Sign. Wave Height\n", "dims: (time:24)\n", "time: 2004-01-01 01:00:00 - 2004-01-03 12:00:10 (24 non-equidistant records)\n", "geometry: GeometryPoint2D(x=225.0, y=1175.0)\n", "values: [0.03508, 0.03571, ..., nan]" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dsi[\"Sign. Wave Height\"].sel(x=250, y=1200)" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkkAAAGwCAYAAAC99fF4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAACNEklEQVR4nOzdeViUVfsH8O8M+77vsoOiorKpAaKSilsZamWallYWabmg5Zb5/tyo1DI3XLJy11LLPcUFRMWN1QVBZUeQRdllmzm/P5CpEdAZGHgG5v5c11xvnDnPM/fty8zcnOc85/AYYwyEEEIIIUQMn+sACCGEEELkERVJhBBCCCGNoCKJEEIIIaQRVCQRQgghhDSCiiRCCCGEkEZQkUQIIYQQ0ggqkgghhBBCGkFFEiGEEEJII6hIIoQQQghpBBVJhBBCCCGNoCKJEEIIIaQRypJ0GjNmjNQn3rx5M0xNTaU+jhBCCCFEHvAk2eCWz+fj3XffhYaGhkQn3bt3LxITE+Hg4NDiAAkhhBBCuCBxkZSbmyvxyJCOjg7i4+OpSCKEEEJIuyXRnKQLFy7A0NBQ4pOeOnUKVlZWzQ6KEEIIIYRrEo0kEUIIIYQoGokmbpeUlEh8Ql1d3WYHQwghhBAiLySek8Tj8SQ6oUAgaHFQhBBCCCFck2gk6cKFC6L/TktLw/z58zF58mR4e3sDAKKiorBjxw6EhIS0TpSEEEIIIW1M6jlJgwYNwieffILx48eLte/duxdbt25FeHi4LOMjhBBCCOGE1EWSpqYm4uPj4ezsLNaenJwMNzc3VFRUyDRAQgghhBAuSHS57b+sra2xefNmrFmzRqx9y5YtsLa2lllg8k4oFOLRo0fQ0dGReL4WIYQQQrjFGENpaSksLS3B5798JSSpi6SffvoJY8eOxenTp/Haa68BAK5evYqHDx/i0KFDzYu4HXr06JFCFYWEEEJIR5KZmYlOnTq9tE+z1knKzMxEaGgo7t27B8YYunXrhqCgIIUqGoqLi6Gvr4/MzExa9oAQQghpJ0pKSmBtbY2ioiLo6em9tC8tJtlMJSUl0NPTQ3FxMRVJhBBCSDshzfe3RNuSvCgyMhITJ06Ej48PsrOzAQC7du3CpUuXmnM6QgghhBC5I3WRdOjQIQwdOhQaGhqIiYlBVVUVAKC0tBQrV66UeYCEEEIIIVyQukhavnw5Nm/ejG3btkFFRUXU7uPjg5iYGJkGRwghhBDCFamLpKSkJPTv379Bu66uLoqKimQREyGEEEII56QukiwsLPDgwYMG7ZcuXYKDg4NMgiKEEEII4ZrURdJnn32GmTNn4tq1a+DxeHj06BH27NmDuXPnYtq0aa0RIyGEEEJIm5O6SPr6668RGBgIf39/lJWVoX///vjkk0/w2Wef4YsvvpA6gE2bNsHe3h7q6urw9PREZGTkS/tHRETA09MT6urqcHBwwObNm5vsu3//fvB4PAQGBoq1h4SEoHfv3tDR0YGpqSkCAwORlJQkdeyEEEII6biatQTAihUrUFBQgOvXr+Pq1avIz8/HsmXLpD7PgQMHMGvWLCxatAixsbHw8/PD8OHDkZGR0Wj/1NRUjBgxAn5+foiNjcXChQsxY8aMRlf6Tk9Px9y5c+Hn59fguYiICEyfPh1Xr15FWFgYamtrERAQgPLycqlzIIQQQkjHxOlikn379oWHhwdCQ0NFbV27dkVgYCBCQkIa9J83bx6OHj2KxMREUVtQUBDi4+MRFRUlahMIBBgwYACmTJmCyMhIFBUV4e+//24yjvz8fJiamiIiIqLRSekAUFVVJVruAPh3xU5aTJIQQghpP1p1Mcny8nIsXrwYPj4+cHJygoODg9hDUtXV1YiOjkZAQIBYe0BAAK5cudLoMVFRUQ36Dx06FDdv3kRNTY2obenSpTAxMcHHH38sUSzFxcUAAENDwyb7hISEQE9PT/RQpC1YCCGEEEUk9Qa3n3zyCSIiIjBp0iRYWFiAx+M164ULCgogEAhgZmYm1m5mZobc3NxGj8nNzW20f21tLQoKCmBhYYHLly9j+/btiIuLkygOxhiCg4PRr18/uLq6NtlvwYIFCA4OFv1cP5JECJFfmU8qYKqrBjVlJa5DIYS0Q1IXSadOncKJEyfg6+srkwBeLLIYYy8tvBrrX99eWlqKiRMnYtu2bTA2Npbo9b/44gskJCS8cksVNTU1qKmpSXROQgj3TiTk4Mt9MbA10sLeqX1hoafBdUiEkHZG6iLJwMDgpZelJGVsbAwlJaUGo0Z5eXkNRovqmZubN9pfWVkZRkZGuHPnDtLS0vDmm2+KnhcKhQAAZWVlJCUlwdHRUfTcl19+iaNHj+LixYvo1KlTi3MihMiHlPwyzDuUACEDUgvKMW7LVez79DVY6VOhRAiRnNRzkpYtW4Zvv/0WFRUVLXphVVVVeHp6IiwsTKw9LCwMPj4+jR7j7e3doP+ZM2fg5eUFFRUVuLi44NatW4iLixM9Ro0aBX9/f8TFxYkujzHG8MUXX+Dw4cM4f/487O3tW5QLIUR+PKsWYNqeGJRV1cLT1gA2hprIeFKBcVuikPmkZZ9bhBDFIvVI0po1a/Dw4UOYmZnBzs5ObP82AFLt3xYcHIxJkybBy8sL3t7e2Lp1KzIyMhAUFASgbh5QdnY2du7cCaDuTrYNGzYgODgYU6dORVRUFLZv3459+/YBANTV1RvMK9LX1wcAsfbp06dj7969OHLkCHR0dESjU3p6etDQoL80CWnPlhy9jXu5pTDWVkXo+x4QMIbxW68irbAC7229ir1T+8LWSIvrMAkh7YDURdKLCzO2xLhx41BYWIilS5ciJycHrq6uOHnyJGxtbQEAOTk5Ymsm2dvb4+TJk5g9ezY2btwIS0tLrFu3DmPHjpXqdeuXHBg4cKBY+2+//YbJkye3KCdCCHf+vJmJP25mgc8D1r3nDlNddQDAgc+8MX7rVaQUlD8vlF6DvTEVSoSQl+N0naT2TJp1Fgghre9ebgkCN15GZY0Qc4Z0xpeDnMWezyutxIRt1/AgrwymOmrY9+lrcDTR5ihaQghXWnWdJEIIkTdlVbWYticGlTVC9O9sgun+Tg36mOqoY/+nr6GLmQ7ySqswbstV3H9cykG0hJD2QqIiydDQEAUFBQD+vbutqQchhLQlxhgWHL6FlPxyWOipY+04N/D5jS8jYqxdN4LU1UIXBWVVeG/rVSTlUqFECGmcRHOSfvrpJ+jo6AAA1q5d25rxEEKIVHZfTcex+EdQ5vOwYYI7DLVUX9rfUEsVez/pi4nbr+HOoxK8tzUKez55Dd0s6bI5IUScREVSfHw83n77baipqcHe3h4+Pj5QVpZ6zjchhMhUQlYRlh2v28tx/nAXeNpKNpptoKWKvZ+8hkm/XkNCVjHe2XwFls/XUKpfr5YHntjPAGClrwE/Z2P4dTaBg7FWs3YcKKmsQWxGEbTVlCSOlxDCDYkmbquoqCArKwtmZmZQUlJCTk4OTE1N2yI+uUUTtwnhVnFFDUauj0TW02cI6GaGLZM8pS5aip/VYPJv1xGbUST161vpa6CfkzH8OhvD19EYBo2MYDHGkPX0GW6mP8HNtKeITn+KpMelYAzg84Cw4AE0eZyQNibN97dEw0F2dnZYt24dAgICwBhDVFQUDAwMGu3bv39/6SMmhBApMMYw5894ZD19BhtDTax6p1ezRnX0NFTw52feiM8qQlWt8PnJ//0fJvpvBiEDEnNKEHk/HzdSnyK76BkO3MzEgZuZ4PGAnlZ68HM2gaetAR7mlyE6va4oyiutavC6Kko81AgYLibnU5FEiByTaCTp77//RlBQEPLy8sDj8dDUITweDwKBQOZByiMaSSKEO1svPsTKk/egqszH4c994Gql16av/6xagGuphbh0vwCR9wuQ9JK75FSUeOhuqQcvWwN42RnAw9YAh2Oy8d2pexjSzQzbPvBqw8gJITIfSQoMDERgYCDKysqgq6uLpKQkhb/cRgjhxo20J/j+nyQAwJI3u7V5gQQAGqpKGNjFFAO71H0O5hZX4tKDAkTez8et7GI4GGvBw9YAXraG6NlJD+oqSmLH+zgaAQCuphRCIGRQauJuPEIIt6Safa2trY0LFy7A3t6eJm4TQtpcYVkVvtgbA4GQ4S03S0zoY8N1SAAAcz11vO3ZCW97SrZRdndLPeioK6O0shZ3HhWjZyf91g2QENIsUi8mOWDAAKSnp+Obb77B+PHjkZeXBwD4559/cOfOHZkHSAghACAQMsw6EIfHJVVwNNHCytE9mjUPSR4o8Xl4zaFuNOnyg0KOoyGENEXqIikiIgI9evTAtWvXcPjwYZSVlQEAEhISsGTJEpkHSAghALD+/H1E3i+AhooSQid6QkutfY9m119yu/KwgONICCFNkbpImj9/PpYvX46wsDCoqv57y6u/vz+ioqJkGhwhhADApfsF+PncfQDAitGu6Gymw3FELefjaAygbo5Vdf2ddYQQuSJ1kXTr1i2MHj26QbuJiQkKC2nYmBAiW7nFlZi5PxaMAe/1tsYYD8nm/ci7zmbaMNJSRWWNEHGZRVyHQwhphNRFkr6+PnJychq0x8bGwsrKSiZBEUIIANQKhPhyXwwKy6vRzUIX/xvVneuQZIbH48GbLrkRItekLpImTJiAefPmITc3FzweD0KhEJcvX8bcuXPxwQcftEaMhBAFtepMEm6kPYW2mjI2ve/R4Fb69q7+ktuVhzQKT4g8krpIWrFiBWxsbGBlZYWysjJ069YN/fv3h4+PD7755pvWiJEQooDO3n2MLREpAIBVb/eEnbEWxxHJnq9T3UhSbMZTPKtWjIV4CWlPpL49REVFBXv27MHSpUsRGxsLoVAId3d3ODs7t0Z8hBAFlPmkAnP+jAcATPG1w/AeFhxH1DpsDDVhpa+B7KK6/d38nE24DokQ8h/NvofW0dERjo6OsoyFEEJQVSvA9L0xKH5WAzdrfSwY3pXrkFpN/bykg9FZuPKwkIokQuSMREVScHCwxCf88ccfmx0MIYSsPJGIhKxi6GuqYOP7HlBVlnpWQLviU18kPaDJ24TIG4mKpNjYWLGfo6OjIRAI0KVLFwBAcnIylJSU4OnpKfsICSEK41j8I+yISgcA/PhuL1jpa3AcUeurv8PtVnYxip/VQE9DheOICCH1JCqSLly4IPrvH3/8ETo6OtixYwcMDAwAAE+fPsWUKVPg5+fXOlESQjq8lPwyzD+UAACYNtARr7uYcRxR27DQ04CDsRZSCspxPfUJhnRTjLwJaQ+kHsdes2YNQkJCRAUSABgYGGD58uVYs2aNTIMjhCiGZ9UCTNsTg/JqAfraGyJ4SGeuQ2pTtF4SIfJJ6iKppKQEjx8/btCel5eH0tJSmQRFCFEsS47exr3cUhhrq2H9eHcoK3XseUgvql8vKYrWSyJErkj9STR69GhMmTIFBw8eRFZWFrKysnDw4EF8/PHHGDNmTGvESAjpwP68mYk/bmaBzwPWjXeDqa461yG1udccDAEA93JLUVBWxXE0hJB6UhdJmzdvxsiRIzFx4kTY2trC1tYW77//PoYPH45Nmza1RoyEkA7qXm4JFh+5DQCYPbizaERF0Rhpq8HFvG7T3qspNJpEiLyQukjS1NTEpk2bUFhYiNjYWMTExODJkyfYtGkTtLQ63oq4hJDWUVZVi2m7Y1BZI0T/ziaY7u/EdUic8nWiLUoIkTfNvvCvpaWFnj17olevXlQcEUKkwhjD/EMJSCkoh4WeOtaOcwOfz+M6LE75PJ+8TfOSCJEfUq+4XV5eju+++w7nzp1DXl4ehEKh2PMpKSkyC44Q0jHtvpqO4wk5UObzsGGCBwy1VLkOiXN97A2hxOchtaAcj4qewVIB1ogiRN5JPZL0ySefYPv27fDz88MXX3yBmTNnij2ktWnTJtjb20NdXR2enp6IjIx8af+IiAh4enpCXV0dDg4O2Lx5c5N99+/fDx6Ph8DAwBa/LiFENhKyirDseCIAYP5wF3jaGrziCMWgo66CHlZ6AOiSGyHyQuqRpFOnTuHEiRPw9fVt8YsfOHAAs2bNwqZNm+Dr64stW7Zg+PDhuHv3LmxsbBr0T01NxYgRIzB16lTs3r0bly9fxrRp02BiYoKxY8eK9U1PT8fcuXMbXeBS2tclhMhGcUUNpu2JQbVAiKHdzfBxP3uuQ5IrPo5GiMsswpWHBXjbsxPX4RCi8HiMMSbNAfb29jh58iS6dm35ppN9+/aFh4cHQkNDRW1du3ZFYGAgQkJCGvSfN28ejh49isTERFFbUFAQ4uPjERUVJWoTCAQYMGAApkyZgsjISBQVFeHvv/9u9usCQFVVFaqq/r01t6SkBNbW1iguLoaurm6z8idEkTDGMHXnTZxNzIONoSaOfdmPtuB4waX7BZi4/Ros9NRxZf7r4PEUe54WIa2hpKQEenp6En1/S325bdmyZfj2229RUVHR7AABoLq6GtHR0QgICBBrDwgIwJUrVxo9JioqqkH/oUOH4ubNm6ipqRG1LV26FCYmJvj4449l8roAEBISAj09PdHD2tr6lTkSQv619WIKzibmQVWZj03ve1CB1AhPWwOoKvGRU1yJtMKWfcYSQlpO6stta9aswcOHD2FmZgY7OzuoqIh/0MXExEh0noKCAggEApiZie9TZGZmhtzc3EaPyc3NbbR/bW0tCgoKYGFhgcuXL2P79u2Ii4uT2esCwIIFCxAcHCz6uX4kiRDyajfSnuCH00kAgCVvdoPr87k3RJyGqhLcbfRxLfUJrjwsgL0x3TlMCJekLpIamwTdEi8OJzPGXjrE3Fj/+vbS0lJMnDgR27Ztg7Hxyxelk/Z11dTUoKam9tJzEkIaKiirwhd7YyAQMrzlZokJfWje38v4OBo/L5IK8X5fW67DIUShSV0kLVmyRCYvbGxsDCUlpQajN3l5eQ1GeeqZm5s32l9ZWRlGRka4c+cO0tLS8Oabb4qer1+iQFlZGUlJSbC2tpb6dQkhzSMQMszaH4fHJVVwMtXGytE9aJ7NK/g4GeGns8DVh4UQCpnCrx9FCJc420VSVVUVnp6eCAsLE2sPCwuDj49Po8d4e3s36H/mzBl4eXlBRUUFLi4uuHXrFuLi4kSPUaNGwd/fH3FxcbC2tm7W6xJCmmf9+fu49KAAGipKCH3fA1pqUv9dpnB6ddKHhooSCsurkZxHm4YTwiWpP7H4fP5L/xIUCAQSnys4OBiTJk2Cl5cXvL29sXXrVmRkZCAoKAhA3Tyg7Oxs7Ny5E0DdnWwbNmxAcHAwpk6diqioKGzfvh379u0DAKirq8PV1VXsNfT19QFArP1Vr0sIabnI+/n4+dx9AMCK0a5wNtPhOKL2QVWZjz72hohIzseVB4VwMae7ZwnhitRF0l9//SX2c01NDWJjY7Fjxw783//9n1TnGjduHAoLC7F06VLk5OTA1dUVJ0+ehK1t3XX4nJwcZGRkiPrXLz8we/ZsbNy4EZaWlli3bl2DNZJa+rqEkJbJLa7ErP1xYAwY38caYzxozR9p+Dga1RVJDwvxEa0lRQhnpF4nqSl79+7FgQMHcOTIEVmcTu5Js84CIYqkRiDEhG1XcSPtKbpZ6OLwNB+oqyhxHVa7ciurGG9uuAQdNWXEfjsEykqczYwgpMNp1XWSmtK3b1+cPXtWVqcjhLRTq88k4UbaU2irKWPT+x5UIDVDN0td6Koro7SqFrcflXAdDiEKSyZF0rNnz7B+/Xp06kRD6oQosrC7j7Elom6T61Vv94QdrfPTLEp8Hl5zMAIAXHlYwHE0hCguqeckGRgYiE3cZoyhtLQUmpqa2L17t0yDI4S0H5lPKjDnjzgAwBRfOwzvYcFtQO2cj6MRztx9jAv38jBtoBPX4RCikKQuktauXSv2M5/Ph4mJCfr27QsDA9rNmxBFVFUrwPS9MSiprIWbtT4WDG/53o6KLqC7OZafSMSNtKeIyyyCm7U+1yERonCkLpI+/PDD1oiDENKOrTiRiISsYuhrqmDj+x5QVaaJxi1lqa+Bt9yscCgmC6HhD7BlkhfXIRGicCT6JEtISBCtXC2JO3fuoLa2ttlBEULaj2Pxj7AzKh0A8NO7brDS1+A4oo4jaIADAOD0ncd4QAtLEtLmJCqS3N3dUVhYKPFJvb29xdY3IoR0TA/zyzD/UAIAYNpAR/i7mHIcUcfibKaDgG512yVtfj4hnhDSdiS63MYYw+LFi6GpqSnRSaurq1sUFCFE/j2rFmD6nhiUVwvQ194QwUM6cx1Sh/T5QEecufsYf8dmY/aQzjRSR0gbkqhI6t+/P5KSkiQ+qbe3NzQ06I1MSEf27ZHbuJdbCmNtNawf704LHrYSdxsDeDsYISqlEL9EpmDJm925DokQhSFRkRQeHt7KYRBC2pM/bmbiz+gs8HnAuvFuMNVV5zqkDm2avyOiUgqx/3omvnzdGYZaqlyHRIhCoD/9CCFSuZdbgm+P3AYABA/pDB9HY44j6vj6ORnD1UoXz2oE+P1KWqu+Vl5pJWS0WxUh7R4VSYQQiZVW1mDa7hhU1ggxoLMJLXLYRng8nujfeseVNJRVyf7u4coaAWYfiEOfFedavRAjpL2gIokQIhHGGBYcvoWUgnJY6Knjp3Fu4PN5rz6QyMTQ7uZwMNZC8bMa7L8u27uH80or8d7Wq/grNhsAsPtqOo0mEQIqkgghEtp1NR3HE3KgzOdhwwQPmhfTxpT4PHz2fN2kbZEpqKoVyOS8t7OL8daGy4jLLIKehgpUlfl4mF+OOzLYWDetoByVNbKJkxAuSF0kXbx4sdGFImtra3Hx4kWZBEUIkS/xmUVYdvwuAGD+cBd42tIWRFwIdLeCua46HpdU4e/noz4t8c/tHLyzOQo5xZVwNNHCkem+GNy1bq2ro/GPWnTuyPv5GLg6HF/ui21xnIRwReoiyd/fH0+ePGnQXlxcDH9/f5kERQiRH8UVNZi+NwY1Aoah3c3wcT97rkNSWGrKSvjEr+7ff3NECgTC5l0SY4xhw/n7CNodg2c1AvTvbILD03xhZ6yFUb2sAABH4x5B2MzzA8Bvl9MAAGF3H+NGWsPvDELaA6mLJMYYeLyG8xAKCwuhpaUlk6AIIfKBMYY5f8Yh6+kz2Bhq4oe3ezX6/idtZ3wfG+hpqCC1oByn7+RKfXxljQAz98dh9ZlkAMBkHzv8+qEX9DRUAAADu5hAR10ZuSWVuN7M4ian+BnCk/JEP685I/k6e4TIE4k3uB0zZgyAurssJk+eDDU1NdFzAoEACQkJ8PHxkX2EhBDObL2YgrOJeVBV5mPT+x6iL1LCHS01ZXzoY4d15+5jU/gDDHc1l7hwzSupxNRd0YjPLIIyn4elb7liQl8bsT7qKkoY4WqBAzczcSTuEV5zMJI6xj9uZEHIABdzHaTkl+NqyhNceVAAHydaLoK0LxKPJOnp6UFPTw+MMejo6Ih+1tPTg7m5OT799FPs3r27NWMlhLSh66lP8MPpuhGAJW92g6uVHscRkXqTfeygoaKE29kliLxfINExt7OL8dbGy4jPLIK+pgp2ftynQYFU7y03SwDAyVs5qK6VfHNzABAIGf64mQkACBrgKHqN1WeS6I450u5IPJL022+/AQDs7Owwd+5curRGSAdWUFaFL/fFQCBkCHSzxIQ+jX+ZEm4YaqnivT7W+O1yGkLDH6J/Z5NG+xWWVeHk7Vwci3+EG2lPwBjgaKKF7R/2hp1x05/hfR2MYKqjhrzSKkQk52PI8012JRF5Px/ZRc+gp6GCYa7m8HE0wr7rGYjJKEJ4cj78u9AmyKT9kHpO0pIlS6hAIqQDEwgZZu2Pw+OSKjiZamPF6B40D0kOTfVzgDKfh6iUQsRmPBW1l1TW4GB0Fj749Tr6rDyHxX/fxvXUugJpaHcz/DXd96UFElC33MCbvepGk47ESXcX3f7rdaNIo92toK6iBFNddXzoYwcA+PFMMo0mkXZF4pGkeo8fP8bcuXNx7tw55OXlNfiFFwhoTQxC2rP15+/j0oMCaKgoIfR9D2ipSf0xQdqApb4GAt2tcDA6C+vPP8AYDysci3+EC0n5YpfIeljpYVQvS4zsaQFLfck3Hn/LzRLbL6XibOJjlFXVQluC34P80iqcTXwMoG6Ceb3P+jtg99V03Mouxpm7jzG0u7kUmRLCHak//SZPnoyMjAwsXrwYFhYW9BcmIR1I5P18/HzuPgBgxWhXOJvpcBwReZmgAY44FJOF8/fycP7ev3eTOZlqY1QvS7zZyxL2rxg1akoPKz04GGshpaAcYXdzMdq90yuPORidhVohg7uNPrqY//u7Y6Stho987bHhwgP8eCYZQ7qa0WrtpF2Quki6dOkSIiMj4ebm1grhEEK4kltciVn748AYML6PNcZ4vPpLkXDLyVQbgW5W+Cs2G9aGGnizZ11h5GKu0+I/YHk8Hka5WWLt2fs4EvfolUUSYwwHbtRtlzK+d8M5bFP9HLAjKg1Jj0tx4laO6HIeIfJM6jlJ1tbWdE2ZkA6mRiDEl/tiUFhejW4WuljyZneuQyIS+n5sT1z8yh8Xv/LH18Nc0NVCV2Yj/KOeFzKR9wtQUFb10r5RKYVIK6yAtpoy3uhl0eB5PU0VTPWr21blp7PJqBVId9ccIVyQukhau3Yt5s+fj7S0tFYIhxDChdWnk3Aj7Sl01JSx6X0PqKsocR0SkZCqMh82RpqtMvXBwUQbPTvpQSBkOHkr56V96ydsj3KzhKZq4xcppvjaQV9TBSn55TgS17JtTwhpCxJdbjMwMBB7A5aXl8PR0RGamppQURFfXK6xLUsIIfIr7O5jbLmYAgD44e2er7zziSiWUb0skZBVjCNxj/CBt12jfZ6UV+Of23Wrfzd2qa2ejroKPuvviO//uYefz93HKDdLqCjRPutEfklUJK1du7bVAti0aRNWrVqFnJwcdO/eHWvXroWfn1+T/SMiIhAcHIw7d+7A0tISX3/9NYKCgkTPHz58GCtXrsSDBw9QU1MDZ2dnzJkzB5MmTRL1qa2txf/+9z/s2bMHubm5sLCwwOTJk/HNN9+Az6c3LFEcmU8qMOePOAB1f+UP79HwMglRbKN6WWLFyUREpz9F5pMKWBtqNuhzOCYL1QIhulvqokenly86+qGPLbZfSkHGkwocjM4SuwuOEHkjUZH04YcftsqLHzhwALNmzcKmTZvg6+uLLVu2YPjw4bh79y5sbBq+cVJTUzFixAhMnToVu3fvxuXLlzFt2jSYmJhg7NixAABDQ0MsWrQILi4uUFVVxfHjxzFlyhSYmppi6NChAIDvv/8emzdvxo4dO9C9e3fcvHkTU6ZMgZ6eHmbOnNkquRIib6pqBZi+NwYllbVws9bHguFduQ6JyCFTXXX4OBrh8oNCHI1/hOn+TmLPM8aw/0bdpbb3JCh4NFWV8flAJyw7fhfrz93HGA8rqCnT5V0in3hMylnYJSUljZ+Ix4OamhpUVVUlPlffvn3h4eGB0NBQUVvXrl0RGBiIkJCQBv3nzZuHo0ePIjExUdQWFBSE+Ph4REVFNfk6Hh4eGDlyJJYtWwYAeOONN2BmZobt27eL+owdOxaamprYtWuXRLGXlJRAT08PxcXF0NXVlegYQuTJt0duY2dUOvQ1VXBihh+spFhDhyiWP25k4utDCehipoPTs/uLPXcz7Qne3hwFDRUlXFs0CLrqr97fr7JGgIGrwpFbUon/G9VdtNgkIW1Bmu9vqa8t6evrw8DAoMFDX18fGhoasLW1xZIlSyAUvvzOherqakRHRyMgIECsPSAgAFeuXGn0mKioqAb9hw4dips3b6KmpqZBf8YYzp07h6SkJPTv/+8bu1+/fjh37hySk+t2wY6Pj8elS5cwYsSIJuOtqqpCSUmJ2IOQ9upY/CPsjEoHAPw0zo0KJPJSQ13NoarER9LjUtzLFf/s2/d8wvbInhYSFUhA3Sa601+vG5HacOEBnlXTIsREPkldJP3++++wtLTEwoUL8ffff+Ovv/7CwoULYWVlhdDQUHz66adYt24dvvvuu5eep6CgAAKBAGZm4nsCmZmZITc3t9FjcnNzG+1fW1uLgoJ/N3ksLi6GtrY2VFVVMXLkSKxfvx5DhgwRPT9v3jyMHz8eLi4uUFFRgbu7O2bNmoXx48c3GW9ISIjYpr7W1tYvzY8QefUwvwzzDyUAAKb7O9JeWuSV9DRU4O9Stz/c37H/3pVW/KwGJ27V/Ty+j3SfieO8rGGlr4H80irsvpouu2AJkSGpF5PcsWMH1qxZg3fffVfUNmrUKPTo0QNbtmzBuXPnYGNjgxUrVmDhwoWvPN+Lt60yxl56K2tj/V9s19HRQVxcHMrKynDu3DkEBwfDwcEBAwcOBFA3F2r37t3Yu3cvunfvjri4OMyaNQuWlpZNzr9asGABgoODRT+XlJRQoUTanWfVAkzbHYPyagH62hti9uDOXIdE2olANyucvvMYx+If4euhXcDn83A0LhuVNUI4m2rDw8ZAqvOpKvMxc5Azvj6UgNCIh5jQ14a2wCFyR+rfyKioKGzevLlBu7u7u2heUL9+/ZCRkfHS8xgbG0NJSanBqFFeXl6D0aJ65ubmjfZXVlaGkZGRqI3P58PJqW4o183NDYmJiQgJCREVSV999RXmz5+P9957DwDQo0cPpKenIyQkpMkiSU1NDWpqai/NiRB59+2R20h6XApjbTWsH+8OZbr9mkjI38UUOmrKyC56huiMp/CyNRBdanuvj02z1mka42GFTeEPkFZYgd+vpDWYFE4I16T+hOzUqZPYhOd627dvF42sFBYWwsDg5X9VqKqqwtPTE2FhYWLtYWFh8PHxafQYb2/vBv3PnDkDLy+vBus1/RdjDFVV/64WW1FR0eBWfyUlpVfOoyKkPfvjZib+jM4CnwesG+8GU111rkMi7Yi6ihKGutZtTHskLhu3sotxN6cEqkp8jHG3atY5lZX4mPV8NHPrxRSUVDacW0oIl6QeSVq9ejXeeecdnDp1Cr179waPx8ONGzdw7949HDx4EABw48YNjBs37pXnCg4OxqRJk+Dl5QVvb29s3boVGRkZonWPFixYgOzsbOzcuRNA3Z1sGzZsQHBwMKZOnYqoqChs374d+/btE50zJCQEXl5ecHR0RHV1NU6ePImdO3eK3UH35ptvYsWKFbCxsUH37t0RGxuLH3/8ER999JG0/xyEtAuJOSVY/PdtAEDwkM7wcTTmOCLSHr3lZomD0Vk4kZCDqpq6PyqHuZrDQEvyu5pf9GYvS2y88AD388qwPTIVs4fQJWAiP6QukkaNGoWkpCRs3rwZycnJYIxh+PDh+Pvvv2FnZwcA+PzzzyU617hx41BYWIilS5ciJycHrq6uOHnyJGxtbQEAOTk5Ypft7O3tcfLkScyePRsbN26EpaUl1q1bJ1ojCahbDXzatGnIysqChoYGXFxcsHv3brGibf369Vi8eDGmTZuGvLw8WFpa4rPPPsO3334r7T8HIXKvtLIG0/fEoKpWiAGdTTBtIF3SIM3j7WAEY201FJRV4c/oLADAe1JO2H6REp+H2UM6Y9qeGGy/lIrJPnYtKroIkSWp10kidWidJNIeMMbw5b5YHE/IgYWeOk7M8IMhfQGRFvi/Y3fw2+U0AICdkSYuzB3Y4n3jhEKGkesvITGnBJ8PdMS8YS4yiJSQxknz/S3RSFJCQgJcXV3B5/ORkJDw0r49e/aUPFJCSKvadTUdxxNyoMznYcMEDyqQSIu95WYlKpLG9W7ehO0X8fk8zBnSGZ/svInfL6fhI197mOjQjTKEexIVSW5ubsjNzYWpqSnc3NzA4/HQ2AAUj8eDQECLghEiD+Izi7Ds+F0AwPzhLvC0le4WbUIa06uTHjxs9JHx5Bne8eoks/MO6mqKXtb6iM8swuaIh1j8RjeZnZuQ5pKoSEpNTYWJiYnovwkh8q24ogbT9sSgRsAwtLsZPu5nz3VIpIPg8Xg48Jk3BEIGdRXZ7bnG49WNJn3w63XsupqOqX4OMNejOzAJtyQqkuonUr/434QQ+SMUMsz5Mw7ZRc9gY6iJH97uJZNLIoTUU1HiQ4b1kYifszF62xngRtpTbLzwAMsCXWX/IoRIoVkrye3atQu+vr6wtLREenrdcvJr167FkSNHZBocIUR62yJTcDYxD6rKfGx63wN6GpLtp0UI13g8HuYEdAEA7L+RgaynFRxHRBSd1EVSaGgogoODMWLECBQVFYnmIOnr62Pt2rWyjo8QIoXrqU/ww+kkAMD/3uwOVys9jiMiRDqvORjB18kINQKG9ececB0OUXBSF0nr16/Htm3bsGjRIigp/Tve6uXlhVu3bsk0OEKI5ArKqvDlvhgIhAyBbpZSbzhKiLwIHlI3mnQwJgtpBeUcR0MUmdRFUmpqKtzd3Ru0q6mpobycfpkJ4YJAyDBrfxwel1TByVQbK0b3oHlIpN3ytDWAfxcTCIQMP5+7z3U4RIFJXSTZ29sjLi6uQfupU6fQrRvdskkIF9adu49LDwqgoaKE0Pc9aDd10u7Vjyb9HZeN+49LOY6GKCqpP0m/+uorTJ8+HZWVlWCM4fr169i3bx9CQkLwyy+/tEaMhJCXiLyfj3Xn6/7aXjnGFc5mOhxHREjL9eikh6HdzXD6zmOsPXsfG9/34DokooCkLpKmTJmC2tpafP3116ioqMCECRNgZWWFn3/+Ge+9915rxEgIaUJucSVm7Y8DY8D4PtYY7S67xf0I4drsIZ1x5u5jnLiVg+mPStDNkraAIm2rRXu3FRQUQCgUwtTUVJYxtQu0dxvhWo1AiAnbruJG2lN0s9DF4Wk+Ml3cjxB58OW+WByLf4TBXc3wy4deXIdDOgBpvr+btU5SPWNjY4UskAiRB6tPJ+FG2lPoqClj0/seVCCRDmnWYGfwecDZxMeIzyziOhyiYCS+3Obu7i7R3TIxMTEtCogQ8mphdx9jy8UUAMCqd3rCzliL44gIaR2OJtoY7d4Jh2KysCYsGTs/6sN1SESBSFwkBQYGiv6bMYaQkBAEBQXB0NCwNeIihDQh80kF5vwRBwD4yNcew1wtuA2IkFY2c5AzjsRl42JyPm6kPUFvO/reIW2j2XOSdHR0EB8fDwcHB1nH1C7QnCTChapaAd7ZHIWErGK42+jjwKfeUFVu0VVzQtqFBYdvYd/1DLzmYIj9n3pzHQ5px9psThIhpG2tOJGIhKxi6GuqYMMEDyqQiML48nUnqCrxcTXlCa48KOA6HKIg6BOWkHbiWPwj7Iyq21D6p3FusNLX4DgiQtqOpb4GJvS1AQCsCUtGC27MJkRiVCQR0g48zC/D/EMJAIDp/o7w70J3lRLFM22gI9SU+YhOf4rw5HyuwyEKQOKJ2+vWrRP7uba2Fr///juMjY3F2mfMmCGbyAghAIBn1QJM2x2D8moBXnMwxOzBnbkOiRBOmOqq4wNvW2yLTMWPZ5IxsLMJ7VFIWpXEE7ft7e1ffTIeDykpKS0Oqj2gidukrcz9Mx4Ho7NgrK2GkzP6wVRXneuQCOFMYVkV/H64gIpqAbZM8sTQ7uZch0TaGWm+vyUeSUpNTW1xYIQQ6fxxMxMHo7PA5wHrxrtRgUQUnpG2Gqb42mHjhYf4KSwZQ7qagc+n0STSOmhOEiFyKjGnBIv/vg0ACB7SGT6Oxq84ghDF8KmfI3TUlXEvtxQnbuVwHQ7pwKhIIkQOlVbWYNqeGFTVCjGgswmmDXTiOiRC5Iaepgo+6Ve3Rt9PZ5NRKxByHBHpqKhIIkTOMMYw//AtpBaUw0JPHT+Nc6PLCYS84KN+dtDXVEFKfjmOxD3iOhzSQVGRRIic2XU1HScScqDM52HDBA8YaqlyHRIhckdHXQWf9XcEAPx87j5qaDSJtAIqkgiRI/GZRVh2/C4AYMGIrvC0NeA4IkLk14c+tjDWVkXGkwocjM7iOhzSATWrSHr48CG++eYbjB8/Hnl5eQCAf/75B3fu3JFpcIQokqKKakzbE4MaAcOw7ub4yNeO65AIkWuaqsr4/Pl8vfXn7qOqVsBxRKSjkbpIioiIQI8ePXDt2jUcPnwYZWVlAICEhAQsWbJE6gA2bdoEe3t7qKurw9PTE5GRka98fU9PT6irq8PBwQGbN28We/7w4cPw8vKCvr4+tLS04Obmhl27djU4T3Z2NiZOnAgjIyNoamrCzc0N0dHRUsdPiCwIhQxz/ohHdtEz2Bhq4od3etIieYRI4P2+NjDXVcej4kocuJHJdTikg5G6SJo/fz6WL1+OsLAwqKr+O1fC398fUVFRUp3rwIEDmDVrFhYtWoTY2Fj4+flh+PDhyMjIaLR/amoqRowYAT8/P8TGxmLhwoWYMWMGDh06JOpjaGiIRYsWISoqCgkJCZgyZQqmTJmC06dPi/o8ffoUvr6+UFFRwalTp3D37l2sWbMG+vr60v1jECIjWyNTcO5eHlSV+dj0vgd01VW4DomQdkFdRQnTX68bTdpw/gEqa2g0iciOxCtu19PW1satW7dgb28PHR0dxMfHw8HBAWlpaXBxcUFlZaXE5+rbty88PDwQGhoqauvatSsCAwMREhLSoP+8efNw9OhRJCYmitqCgoIQHx//0gLNw8MDI0eOxLJlywDUFXqXL19+5ajVy9CK20RWrqc+wfhtVyEQMqwc3UO0iSchRDLVtUL4rw5HdtEzfDOyKz7xc+A6JCLHpPn+lnokSV9fHzk5DRfvio2NhZWVlcTnqa6uRnR0NAICAsTaAwICcOXKlUaPiYqKatB/6NChuHnzJmpqahr0Z4zh3LlzSEpKQv/+/UXtR48ehZeXF9555x2YmprC3d0d27Zte2m8VVVVKCkpEXsQ0lIFZVX4cl8MBEKGQDdLjO9jzXVIhLQ7qsp8zBzkDADYFP4Q5VW1HEdEOgqpi6QJEyZg3rx5yM3NBY/Hg1AoxOXLlzF37lx88MEHEp+noKAAAoEAZmZmYu1mZmbIzc1t9Jjc3NxG+9fW1qKgoEDUVlxcDG1tbaiqqmLkyJFYv349hgwZIno+JSUFoaGhcHZ2xunTpxEUFIQZM2Zg586dTcYbEhICPT090cPamr7MSMsIhAyz9sfhcUkVnEy1sWJ0D5qHREgzjfGwgp2RJp6UV+P3K2lch0M6CKmLpBUrVsDGxgZWVlYoKytDt27d0L9/f/j4+OCbb76ROoAXvxQYYy/9omis/4vtOjo6iIuLw40bN7BixQoEBwcjPDxc9LxQKISHhwdWrlwJd3d3fPbZZ5g6darYZb8XLViwAMXFxaJHZiZNECQts+7cfVx6UAANFSWEvu8BLTWJt1IkhLxAWYmPWYM7AwC2XkxBSWXDqwuESEvqT2UVFRXs2bMHS5cuRWxsLIRCIdzd3eHs7CzVeYyNjaGkpNRg1CgvL6/BaFE9c3PzRvsrKyvDyMhI1Mbn8+HkVDeRz83NDYmJiQgJCcHAgQMBABYWFujWrZvYebp27So2AfxFampqUFNTkzg/Ql7mYnI+1p2/DwBYOcYVzmY6HEdESPv3Zi9LbLzwAPfzyrA9MhWzh3TmOiTSzjVrCQAAcHR0xNtvv413331X6gIJAFRVVeHp6YmwsDCx9rCwMPj4+DR6jLe3d4P+Z86cgZeXF1RUmr4biDGGqqoq0c++vr5ISkoS65OcnAxbW1tp0yBEarnFlZh1IA6MAeP72GC0eyeuQyKkQ1Di80SF0fZLqXhaXs1xRKS9k7pIGjJkCGxsbDB//nzcvn27RS8eHByMX375Bb/++isSExMxe/ZsZGRkICgoCEDdJa7/znMKCgpCeno6goODkZiYiF9//RXbt2/H3LlzRX1CQkIQFhaGlJQU3Lt3Dz/++CN27tyJiRMnivrMnj0bV69excqVK/HgwQPs3bsXW7duxfTp01uUDyGvUiMQ4ou9MXhSXo1uFrpY8ma3Vx9ECJHYsO7m6Gqhi7KqWmyNTOE6HNLeMSnl5+ez9evXMx8fH8bj8ViPHj3Y999/zzIzM6U9FWOMsY0bNzJbW1umqqrKPDw8WEREhOi5Dz/8kA0YMECsf3h4OHN3d2eqqqrMzs6OhYaGij2/aNEi5uTkxNTV1ZmBgQHz9vZm+/fvb/C6x44dY66urkxNTY25uLiwrVu3ShV3cXExA8CKi4ulOo4otpUn7jLbeceZ67f/sNT8Mq7DIaRDCruTy2znHWcu35xieSWVXIdD5Iw0399Sr5P0X6mpqdi7dy/27duHe/fuoX///jh//rzsKjg5RuskEWmF3X2MqTtvAgA2T/TAMFcLjiMipGNijCFw0xXEZxbh4372WPwGjdiSf7XqOkn/ZW9vj/nz5+O7775Djx49RPOVCCHiMp9UYM4fcQCAj3ztqUAipBXxeDzMeT43affVdOQWS77IMSH/1ewi6fLly5g2bRosLCwwYcIEdO/eHcePH5dlbIR0CFW1AkzfG4OSylq42+hj/nAXrkMipMPzczZGbzsDVNUKsfHCA67DIe2U1EXSwoULYW9vj9dffx3p6elYu3YtcnNzsXv3bgwfPrw1YiSkXVt+PBEJWcXQ11TBhgkeUFVu0QAuIUQCPB4PwUO6AAD238hA1tMKjiMi7ZHUn9bh4eGYO3cusrOzceLECUyYMAGampqtERsh7d7R+EfYdTUdAPDTODdY6WtwHBEhisPb0Qi+TkaoETCsP0ejSUR6Ui8m2dS+aoQQcQ/zy7DgUAIAYLq/I/y7mHIcESGKJ3hIF1x+cAUHY7Lw+UBH2BlrcR0SaUeavQ/C3bt3kZGRgepq8cW6Ro0a1eKgCGnvnlULMG13DMqrBXjNwRCzB9PKv4RwwdPWAP5dTHAhKR8/n7uPn8a5cR0SaUekLpJSUlIwevRo3Lp1Czwer8HeaQKBQLYREtIOLT5yG0mPS2GsrYZ177lDWYnmIRHCleAhXXAhKR9/x2Vj2kBH2gaISEzqT+6ZM2fC3t4ejx8/hqamJu7cuYOLFy/Cy8tLbBNZQhTVHzczcTA6C3wesH68O0x11bkOiRCF1qOTHoZ2NwNjwNqz97kOh7QjUhdJUVFRWLp0KUxMTMDn88Hn89GvXz+EhIRgxowZrREjIe1GYk4JFv9dt13PnIAu8HY0esURhJC2MHtIZ/B4wIlbObj7qITrcEg7IXWRJBAIoK2tDQAwNjbGo0ePAAC2trYNNo0lRJGUVtZg2p4YVNUKMbCLCT4f4Mh1SISQ51zMdfFGT0sAwI9hyRxHQ9oLqYskV1dXJCTU3bHTt29f/PDDD7h8+TKWLl0KBwcHmQdISHvAGMP8Q7eQWlAOSz11/PSuG/h8HtdhEUL+Y9ZgZ/B5wNnEx4jPLOI6HNIOSF0kffPNNxAKhQCA5cuXIz09HX5+fjh58iTWrVsn8wAJaQ92RqXjxK0cKPN5WD/BAwZaqlyHRAh5gaOJNgLdrQDQaBKRjMR3t509exa+vr4YOnSoqM3BwQF3797FkydPYGBgILrDjRBFEp9ZhOUn7gIAFozoCk9bA44jIoQ0ZeYgZxyNe4SI5HzcTHsCLztDrkMickzikaSAgAAYGBigf//+WLJkCcLDw0VrJBkaGlKBRBRSUUU1pu2JQY2AYVh3c3zka8d1SISQl7A10sI7Xp0AAGvO0GgSeTmJi6TMzExs27YNnTt3xu7du/H6669DX18fgwYNwvLly3H58mXU1ta2ZqyEyBWhkGHOH/HILnoGWyNN/PBOT/pjgZB24IvXnaGqxEdUSiGuPCjgOhwix3isfjVIKWVmZuLChQsIDw9HeHg40tPToampidLSUlnHKJdKSkqgp6eH4uJi6Orqch0O4cDmiIf47tQ9qCrzcfhzH7ha6XEdEiFEQkuO3MaOqHR42hrgYJA3/YGjQKT5/m72MsDW1tbw9fWFt7c3vL29oa2tjWbWW4S0O9dSCrHqdN2SF/97szsVSIS0M9P9naCmzEd0+lOEJ+dzHQ6RU1IVSSkpKfj1118xadIkdOrUCR4eHjh8+DBcXV1x6tQpPH36tLXiJERuFJRV4ct9sRAIGUa7W2F8H2uuQyKESMlUVx0feNsCAH48k0x/5JNGSXx3m62tLUpKStCvXz/0798fX375JTw9PaGkpNSa8REiVwRChpn7Y5FXWgUnU20sD3SlYXpC2qmgAY7Ycy0Dt7KLcebuYwztbs51SETOSDySVFVVBaBuI1slJSUoKSmBz6dNO4liWXfuPi4/KISGihJC3/eAlprUe0QTQuSEkbYapjy/I/WnsGQIhTSaRMRJXOXk5uYiKioKI0aMwLVr1zBy5EgYGBjgjTfewOrVq3Hjxg3RIpOEdEQXk/Ox7nzd5pgrx7jSTuKEdACf+jlCR10Z93JLceJWDtfhEDnT7LvbACAxMVF0h9vp06fB4/FQVFQkw/DkF93dplhyiysxYl0knpRXY3wfG4SM6cF1SIQQGfn57H38dDYZjiZaODN7AJRoS6EOrU3ubnv8+DESEhKQkJCA+Ph4lJaWii7JEdKR1AiE+GJvDJ6UV6ObhS6WvNmN65AIITL0UT876Guq4GF+OY7EZXMdDpEjEhdJeXl5+OOPPzBt2jR07doVlpaW+PDDD3H37l289957OH/+vMKMIhHFsup0Em6mP4WOmjJCJ3pAXYVuViCkI9FRV8Fn/R0BAGvP3keNgKaOkDoSzzo1NzeHiooKvLy8MHbsWAwcOBC+vr7Q0NBozfgI4dSZO7nYejEFALDqnZ6wNdLiOCJCSGv40McW2y+lIONJBQ5FZ+G9PjZch0TkgMRF0qlTp9CvXz9oadGXBFEMmU8qMPfPeADAR772GOZqwXFEhJDWoqmqjM8HOmHZ8btYd+4+RntYQU2ZRo0VncSX24YOHUoFElEYVbUCTNsTg5LKWrjb6GP+cBeuQyKEtLL3+9rATFcNj4orceBGJtfhEDlACx0R0ojlxxNxK7sY+poq2DDBA6rK9FYhpKNTV1HCF687AwA2nH+AyhoBxxERrnH+yb9p0ybY29tDXV0dnp6eiIyMfGn/iIgIeHp6Ql1dHQ4ODti8ebPY84cPH4aXlxf09fWhpaUFNzc37Nq1q8nzhYSEgMfjYdasWbJIh3QAR+MfYdfVdADAT+PcYKVP8+4IURTjvKxhpa+BvNIq7H7+OUAUF6dF0oEDBzBr1iwsWrQIsbGx8PPzw/Dhw5GRkdFo/9TUVIwYMQJ+fn6IjY3FwoULMWPGDBw6dEjUx9DQEIsWLUJUVBQSEhIwZcoUTJkyBadPn25wvhs3bmDr1q3o2bNnq+VI2peH+WVYcCgBAPCFvxP8u5hyHBEhpC2pKvMxc1DdaNKm8Icor6rlOCLCpRYtJtlSffv2hYeHB0JDQ0VtXbt2RWBgIEJCQhr0nzdvHo4ePYrExERRW1BQEOLj4xEVFdXk63h4eGDkyJFYtmyZqK2srAweHh7YtGkTli9fDjc3N6xdu1bi2GkxyY7nWbUAgRsvI+lxKV5zMMTuj/tCWYnzwVZCSBurFQgx+McIpBVW4KuhXTDd34nrkIgMtcliki+6efMmLl68KHH/6upqREdHIyAgQKw9ICAAV65cafSYqKioBv2HDh2KmzdvoqampkF/xhjOnTuHpKQk9O/fX+y56dOnY+TIkRg8eLBE8VZVVaGkpETsQTqWxUduI+lxKYy11bBuvDsVSIQoKGUlPmYN7gwA2HoxBSWVDb9fiGKQ2bfApEmT4O/vL3H/goICCAQCmJmZibWbmZkhNze30WNyc3Mb7V9bW4uCggJRW3FxMbS1taGqqoqRI0di/fr1GDJkiOj5/fv3IyYmptHRqqaEhIRAT09P9LC2tpb4WCL//riRiYPRWeDzgPXj3WGqo851SIQQDr3ZyxLOptooflaD7ZGpXIdDOCKzIuncuXNISUmR+jgeT3yPHMZYg7ZX9X+xXUdHB3Fxcbhx4wZWrFiB4OBghIeHAwAyMzMxc+ZM7N69G+rqkn8RLliwAMXFxaJHZibdHtpRJOaUYPGR2wCAOQFd4O1oxHFEhBCuKfF5otGkXy+l4ml5NccRES5IvJjkq1haWkrV39jYGEpKSg1GjfLy8hqMFtUzNzdvtL+ysjKMjP79YuPz+XByqruG7ObmhsTERISEhGDgwIGIjo5GXl4ePD09Rf0FAgEuXryIDRs2oKqqCkpKDRcQU1NTg5qamlQ5EvlXWlmDaXtiUFUrxMAuJvh8gCPXIRFC5MRwV3N0tdBFYk4JtkamYN4wWi9N0TSrSBIKhXjw4AHy8vIgFIrvcfPi3J+mqKqqwtPTE2FhYRg9erSoPSwsDG+99Vajx3h7e+PYsWNibWfOnIGXlxdUVFSafC3GmGjz3UGDBuHWrVtiz0+ZMgUuLi6YN29eowUS6ZgYY5h/6BZSC8phqaeOn951A592/yaEPMfn8xA8pDOm7ryJ3y+n4eN+9jDWpj+WFYnURdLVq1cxYcIEpKen48Ub43g8HgQCyRffCg4OxqRJk+Dl5QVvb29s3boVGRkZCAoKAlB3iSs7Oxs7d+4EUHcn24YNGxAcHIypU6ciKioK27dvx759+0TnDAkJgZeXFxwdHVFdXY2TJ09i586dojvodHR04OrqKhaHlpYWjIyMGrSTjm1nVDpO3MqBMp+HDe97wEBLleuQCCFyZnBXU/TqpIf4rGKEhj/E4je6cR0SaUNSF0lBQUHw8vLCiRMnYGFh8dL5Q68ybtw4FBYWYunSpcjJyYGrqytOnjwJW1tbAEBOTo7Ymkn29vY4efIkZs+ejY0bN8LS0hLr1q3D2LFjRX3Ky8sxbdo0ZGVlQUNDAy4uLti9ezfGjRvX7DhJxxOXWYTlJ+4CABaM6AoPGwOOIyKEyCMej4fggC748Nfr2H01HVP9HGCuRzd2KAqp10nS0tJCfHy8aM6PoqJ1ktqvoopqjFx3CdlFzzCsuzlCJ3q0qNgnhHRsjDG8uyUKN9KeYtJrtlgWSFcd2rNWXSepb9++ePDgQbODI4RLQiHDnD/ikV30DLZGmvjhnZ5UIBFCXorH4yF4SBcAwP4bGch6WsFxRKStSH257csvv8ScOXOQm5uLHj16NJgwTVt8EHm2NTIF5+7lQVWZj40TPKCr3vSEf0IIqeftaARfJyNcflCI9ece4Pu36btOEUh9uY3Pbzj4xOPxROsbSTNxuz2jy23tz7WUQkz45RoEQoaVo3tgQl8brkMihLQj0elPMTb0CpT4PJwLHgA7Yy2uQyLNIM33t9QjSamptPIoaX8Kyqrw5b5YCIQMo92tML4PrZhOCJGOp60B/LuY4EJSPn4+dx8/jXPjOiTSyqQukurvPCPkvyprBEjIKsaNtCdIyi3F+D42crNytUDIMHN/LPJKq+Bsqo0Vo11pHhIhpFmCh3TBhaR8/B2Xjen+jnAy1eE6JNKKpC6Sjh492mg7j8eDuro6nJycYG9v3+LAiHwrqqjGzbSnuJH+BDfTnuJWVjGqBf8uLHrnUTHOBg+Qi2Lk53P3cflBITRUlLDpfQ9oqspsoXlCiILp0UkPAd3McObuY/x09j42TvDgOiTSiqT+tggMDBTNQfqv/85L6tevH/7++28YGNDaMx3JtZRCHIl/hJtpT5D8uKzB8yY6auhtZ4Dz9/LwML8c8VnFcLPWb/tA/+Nicj7Wn78PAFg5xhXOZvRXHyGkZYIDOiMs8TFOJORg+sASdLOkeakdldRLAISFhaF3794ICwsTbfYaFhaGPn364Pjx47h48SIKCwsxd+7c1oiXcOTAjQyM33YVe69liAokRxMtvNfbGqvf6YWIrwbi+sJB2PS+J4Z2NwcAHI7J4jJk5BQ/w6wDcWAMGN/HBqPdO3EaDyGkY3Ax18XIHhYAgJ/OJnMcDWlNUo8kzZw5E1u3boWPj4+obdCgQVBXV8enn36KO3fuYO3atfjoo49kGijhzraLKVhxMhEAMLKHBd5ys4SnrQGMmtjDaKxHJxyJe4Sj8Y+waGRXqCm3/X54NQIhvtwbiyfl1ehmoYslb9JWAoQQ2Zk1uDNO3spB2N3HiM8sQi+OR81J65B6JOnhw4eN3jKnq6uLlJQUAICzszMKCgpaHh3hFGMMq07fExVInw1wwIYJ7gjobt5kgQQAvk7GMNNVQ1FFDS7cy2urcMWsOp2Em+lPoaOmjNCJHlBXoY2LCSGy42SqjUB3KwDAj2E0mtRRSV0keXp64quvvkJ+fr6oLT8/H19//TV69+4NALh//z46daJLG+2ZUMjw7ZE72HjhIQDg62FdsGB4V4kmYivxeaIPj0Mx2a0aZ2PO3MnF1ot1Bfuqd3rC1ojWMiGEyN7MQc5Q5vMQkZyPm2lPuA6HtAKpi6Tt27cjNTUVnTp1gpOTE5ydndGpUyekpaXhl19+AQCUlZVh8eLFMg+WtI0agRCz/4jDrqvp4PGA5YGumDZQur36xnrUFckX7uWhsKyqNcJsVEZhBeb8GQ8A+LifPYa5WrTZaxNCFIutkRbe8ar7rFtzhkaTOiKp5yR16dIFiYmJOH36NJKTk8EYg4uLC4YMGSJajTswMFDWcZI2UlkjwLQ9MTh/Lw/KfB5+HOeGUb0spT5PZzMd9LDSw63sYhyLf4TJvq2/LERljQDT98agtLIW7jb6mDfMpdVfkxCi2L543RmHorMRlVKIKw8K4ONkzHVIRIaatWAMj8fDsGHDMGzYMFnHQzhUWlmDj3fcxPXUJ1BT5mPzRE/4u5g2+3xjPaxwK7sYh2Ky26RIWnEiEbeyi2GgqYKNEzygqiz1QCkhhEjFSl8D4/tYY0dUOtaEJcPb0Ugu1ocjsiFRkbRu3Tp8+umnUFdXx7p1617ad8aMGTIJjLStwrIqfPjbddzOLoGOmjK2T+6NPvaGLTrnm70ssfx54ZL8uBSdW3GNoqPxj7DrajoA4MdxbrDU12i11yKEkP+a7u+E/TcyEZ3+FBHJ+RjYpfl/XBL5ItEGt/b29rh58yaMjIxeupo2j8cT3eHW0XWkDW4fl1RiwrareJhfDiMtVez4qA9crfRkcu6pO28i7O5jfDbAAQuGd5XJOV/0IK8MozZcQkW1AF/4O2Hu0C6t8jqEENKUFSfuYltkKnp20sOR6b40miTHZL7B7X83taUNbjuekJOJeJhfDks9dez6pC8cTbRldu6xHp0Qdvcx/o7NxtdDXaDEl+0Hx7NqAabviUFFtQCvORhi1mBnmZ6fEEIkETTAEXuuZSAhqxhhdx8j4PmiuqR9a/GkjdraWpSVNdyigrQPeaWVOHErBwAQOtFTpgUSAPi7mEBfUwWPS6pw+YFs185ijOGbv28j6XEpTHTUsG68O5SVaB4SIaTtGWmrYYqvHYC6dZOEwldepCHtgMTfKCdPnsSuXbvE2lasWAFtbW3o6+sjICAAT58+lXmApHXtvZaBGgGDh41+q6wYq6asJLo77pCMtyn582YWDsVkgc8D1r3nDlMddZmenxBCpDHVzwE6asq4l1uKk7dzuA6HyIDERdLq1atRUlIi+vnKlSv49ttvsXjxYvzxxx/IzMzEsmXLWiVI0jqqa4XYcy0DAPChj12rvU79mkmn7+SitLJGJue8+6gEi4/cBgDMCegCb0cjmZyXEEKaS19TFZ/4OQAAfgpLhoBGk9o9iYuk27dvi+3XdvDgQQwZMgSLFi3CmDFjsGbNGhw7dqxVgiSt4587ucgvrYKJjhqGt+Kiiz076cHRRAuVNUKcupXb4vOVVtZg+t4YVNUKMbCLCT4f4CiDKAkhpOU+6mcHfU0VPMwvx5G4tt9xgMiWxEVSaWkpjIz+/Wv90qVLeP3110U/d+/eHY8ePZJtdKRV7biSBgCY0MemVdcU4vF4GOtZN5p0sIWX3BhjmH/oFlIL6iaa//SuG/gyngxOCCHNpaOugs/61/3htvbsfdQIhBxHRFpC4m9GS0tLJCbWbXRaVlaG+Ph4+Pr6ip4vLCyEpqam7CMkreJ2djGi059Cmc/D+31tWv31RrtbgccDrqc+QeaTimafZ2dUOk7cyoEyn4cN73vAQEtVhlESQkjLfehjC2NtVWQ8qcChaNnOxSRtS+Ii6e2338asWbOwa9cuTJ06Febm5njttddEz9+8eRNdutD6NO1F/SjSiB4WMNVt/QnPFnoa8HWsW67/cDM3vY3LLMLyE3cBAAtGdIWHjYHM4iOEEFnRVFXG58/3u1x37j6qagUcR0SaS+IiacmSJfDy8sKMGTMQFxeH3bt3Q0lJSfT8vn378Oabb7ZKkES2npRX40h83aXR1pyw/aKxnlYAgMOxWZBgDVMxRRXVmL4nBjUChuGu5vjo+a22hBAij97vawMzXTU8Kq7EgRuZXIdDmknivds0NTUbLAHwXxcuXJBJQKT17b+RgepaIXpY6cHDRr/NXndod3Noqd5GemEFotOfwstOsm1PhEKGOX/EI7voGWyNNPH92z1pNVtCiFxTV1HCF/5OWHzkDjacf4B3vayhrqL06gOJXKGV9xRMrUCI3VF1e5x96GPXpsWGpqoyhveou4tOmjWTtlxMwbl7eVBV5mPT+x7QVVdprRAJIURm3u1tDSt9DeSVVmH3870lSftCRZKCOZuYh0fFlTDUUsUbPVvvtv+m1K+ZdDw+B5U1r75Ofy2lEKvPJAEA/m9Ud3S3lM2ecoQQ0trUlJUwY1Dd3KTQ8Icor6rlOCIiLc6LpE2bNsHe3h7q6urw9PREZGTkS/tHRETA09MT6urqcHBwwObNm8WeP3z4MLy8vKCvrw8tLS24ubk1uEwYEhKC3r17Q0dHB6ampggMDERSUpLMc5NH9RO23+vNzdBvX3tDWOlroLSqFmF3Hzd4vlYgxO3sYuyKSsOs/bGYuvMmBEKG0e5WeK+3dZvHSwghLTHGoxPsjDRRWF6N359//pL2g9Mi6cCBA5g1axYWLVqE2NhY+Pn5Yfjw4cjIyGi0f2pqKkaMGAE/Pz/ExsZi4cKFmDFjBg4dOiTqY2hoiEWLFiEqKgoJCQmYMmUKpkyZgtOnT4v6REREYPr06bh69SrCwsJQW1uLgIAAlJeXt3rOXErKLUVUSiGU+DxMfM2Wkxj4fB7GeNRN4D4Uk4WiimpcuJeH1aeTMGHbVfT8vzN4Y/0lLD5yB3/HPUJJZS26WehixWhXmodECGl3VJT4mPl84+2tF1NQIqNdB0jb4DFpbzP6j6ysLFhaWoLPb16t1bdvX3h4eCA0NFTU1rVrVwQGBiIkJKRB/3nz5uHo0aOi9ZoAICgoCPHx8YiKimrydTw8PDBy5Mgmt03Jz8+HqakpIiIi0L9/f4liLykpgZ6eHoqLi6GrqyvRMbJUWFaFuzkl8HU0lngxxYV/3cLeaxkY7mqO0ImerRxh01ILyuG/OrzJ53XUlOFmow9PWwN42hqgt50hTXgkhLRbAiHD0LUX8SCvDDMHOWP2kM5ch6TQpPn+lvjutsZ069YNcXFxcHBwkPrY6upqREdHY/78+WLtAQEBuHLlSqPHREVFISAgQKxt6NCh2L59O2pqaqCiIj6hlzGG8+fPIykpCd9//32TsRQXFwOoG4VqSlVVFaqqqkQ//3cfOy7M+TMe4Un58O9igp/GuUFf8+WLKhY/q8Ffz9cnasvb/htjb6wFbwcjRKUUin72sDEQFUVOptpQolW0CSEdhBKfh9mDO2P63hj8eikVk33saCHcdqJFRVILBqFQUFAAgUAAMzMzsXYzMzPk5ja+v1dubm6j/Wtra1FQUAALi7qJyMXFxbCyskJVVRWUlJSwadMmDBkypMkcgoOD0a9fP7i6ujYZb0hICP7v//5PmhRbzbNqAa48qCswLiTl4431l7B5oidcrZqe1PznzUw8qxGgi5kO+tpLdut9awqd6IG7j0rQxVwHRtpqXIdDCCGtarirObpa6CIxpwRbI1Mwb5gL1yERCXA+cfvFeSaMsZfOPWms/4vtOjo6iIuLw40bN7BixQoEBwcjPDy80fN98cUXSEhIwL59+14a54IFC1BcXCx6ZGZytzjYjbQnqBYIYaytChtDTWQ9fYYxoVdw4Ebjc7mEQoZdV7m57b8p+pqq8HEypgKJEKIQ+Hwegp9fZvv9choKyqpecQSRBy0qkhYuXPjSS1QvY2xsDCUlpQajRnl5eQ1Gi+qZm5s32l9ZWVls810+nw8nJye4ublhzpw5ePvttxud4/Tll1/i6NGjuHDhAjp16vTSeNXU1KCrqyv24MrlBwUAAP8upjj2RT8M7mqK6loh5h26ha8Pxje4tT48OQ/phRXQVVdGoLslFyETQojCG9zVFL066eFZjQCh4Q+5DodIoEVF0oIFC6Cvr9+sY1VVVeHp6YmwsDCx9rCwMPj4+DR6jLe3d4P+Z86cgZeXV4P5SP/FGBObT8QYwxdffIHDhw/j/PnzsLe3b1YOXLn8sK5I6udsDD1NFWyd5IWvhnYBnwf8cTMLY0OvIKPw301kd1ypG0Ua19samqotusJKCCGkmXg8HoID6vY43X01HY9LKjmOiLwKp5fbgoOD8csvv+DXX39FYmIiZs+ejYyMDAQFBQGoK8I++OADUf+goCCkp6cjODgYiYmJ+PXXX7F9+3bMnTtX1CckJARhYWFISUnBvXv38OOPP2Lnzp2YOHGiqM/06dOxe/du7N27Fzo6OsjNzUVubi6ePXvWdsk305Pyatx5VDdp3Of5hrF8Pg/T/Z2w86O+MNRSxZ1HJXhjfSTOJT5GSn4ZIpLzweMBk16z4zByQggh/Z2N4WVrgKpaITZeeMB1OOQVOB1WGDduHAoLC7F06VLk5OTA1dUVJ0+ehK1t3Ro+OTk5Ymsm2dvb4+TJk5g9ezY2btwIS0tLrFu3DmPHjhX1KS8vx7Rp05CVlQUNDQ24uLhg9+7dGDdunKhP/ZIDAwcOFIvnt99+w+TJk1svYRmIelgIxgAXcx2Y6IjP5+nnbIwTM/ph2p4YxGYU4eMdN+Fsqg0AGORiChsjTS5CJoQQ8hyPx8OcgC4Yv+0q9l3PwKf9HdDJgD6b5VWL1klSZFytk7Tg8C3su56Bj/vZY/Eb3RrtU10rxMqTiWKru+78qA/6dzZpoygJIYS8zIRtV3HlYSHe622N78b25DochSLN97fUl9syMjIavfWfMdbkStlEduonbfs6GTXZR1WZj/+N6o6f33ODtpoy+tgbop+TcVuFSAgh5BXmBNTd6fZndBbSCjr2bg/tmdRFkr29PfLz8xu0P3nypN1NgG5vMgorkPGkAsp8HvrYN10k1XvLzQpx3w7B7o/7SrwqNyGEkNbnaWuIgV1MIBAyrDt3n+twSBOkLpKaWseorKwM6urqMgmKNK7+rjZ3G31oq0k2nUxZiQ9VZc6XwyKEEPKCOUPq7nT7Ky4bD/JKOY6GNEbiidvBwcEA6iadLV68GJqa/040EwgEuHbtGtzc3GQeIPnXv5fa6NIZIYS0dz066SGgmxnO3H2Mn87ex8YJHlyHRF4gcZEUGxsLoG4k6datW1BV/XffGVVVVfTq1UvsVnwiW0Ihw5WHdVuR0PwiQgjpGIIDOiMs8TFOJORg+sASdLPkbqFi0pDERdKFCxcAAFOmTMHPP//M6YrTiigxtwRPyquhpaqEXtb6XIdDCCFEBlzMdTGyhwWOJ+Tgp7PJ2PaBF9chkf+QerLKb7/9RgUSB+ovtb3mYAQVJZpjRAghHcWswZ3B5wFhdx8jPrOI63DIf0j9bVteXo7FixfDx8cHTk5OcHBwEHuQ1nHpQd2lNh+61EYIIR2Kk6k2At2tAAA/hiVzHA35L6lX3P7kk08QERGBSZMmwcLCQi52lO/oqmoFuJ5K85EIIaSjmjnIGUfiHiEiOR83057Ay655m8cT2ZK6SDp16hROnDgBX1/f1oiHNCImvQiVNUIYa6uhs5k21+EQQgiRMVsjLbzr1Qn7rmdizZlk7Pv0Na5DImjG5TYDAwMYGlKF25auPF8fqZ+TEY3cEUJIB/XF685QVeIjKqUQV57PQyXckrpIWrZsGb799ltUVFS0RjykEZdofSRCCOnwrPQ1ML6PNQBgTVhyo1uAkbYl0eU2d3d3sRGMBw8ewMzMDHZ2dlBRURHrGxMTI9sIFVxJZY3obgcqkgghpGOb7u+E/TcyEZ3+FBHJ+RjYxZTrkBSaREVSYGBgK4dBmnL1YSGEDHAw0YKlvgbX4RBCCGlFprrq+MDbFtsiU/FjWDIGdDahaRYckqhIWrJkSWvHQZog2orEkUaRCCFEEQQNcMSeaxlIyCpG2N3HCOhuznVICotWJZRzNB+JEEIUi5G2Gqb42gGoWzdJKKS5SVxp9t1tLz6MjIxgZWWFAQMG4LfffmuNWBVObnElHuaXg88DvB2MuA6HEEJIG5nq5wAdNWXcyy3Fyds5XIejsKQukr799lvw+XyMHDkS//d//4f//e9/GDlyJPh8PqZPn47OnTvj888/x7Zt21ojXoVSf6mtRyd96GmqvKI3IYSQjkJfUxWf+NXtYvFTWDIENJrECakXk7x06RKWL1+OoKAgsfYtW7bgzJkzOHToEHr27Il169Zh6tSpMgtUEdUXSf2caBSJEEIUzUf97PDblVQ8zC/HkbhsjPHoxHVICkfqkaTTp09j8ODBDdoHDRqE06dPAwBGjBiBlJSUlkenwBhjNB+JEEIUmI66Cj7tXzea9PO5+6gRCDmOSPFIXSQZGhri2LFjDdqPHTsmWom7vLwcOjo6LY9OgT3IK0NeaRXUVfjwsDHgOhxCCCEcmOxjB2NtVaQXVuBQdBbX4SgcqS+3LV68GJ9//jkuXLiAPn36gMfj4fr16zh58iQ2b94MAAgLC8OAAQNkHqwiqR9F6m1nCHUVJY6jIYQQwgVNVWUEDXDE8hOJWH/+AUZ7WEFNmb4T2orUI0lTp05FREQEtLS0cPjwYRw8eBCampqIiIjAxx9/DACYM2cODhw4IPNgFcllutRGCCEEwMTXbGGmq4bsomc4cCOT63AUitQjSQDg6+sLX19fWcdCnqsVCHE15QkAoB8VSYQQotDUVZTwhb8TFh+5gw3nH+BdL2u6wtBGJBpJKikpEfvvlz1Iy8VnFaOsqhb6miroZqHLdTiEEEI49m5va1jpayCvtAq7r6ZzHY7CkKhIMjAwQF5eHgBAX18fBgYGDR717aTl/rsVCZ9Pe/YQQoiiU1NWwoxBTgCA0PCHKK+q5TgixSDR5bbz58+L7ly7cOFCqwZEaCsSQgghDY3x6ITQ8IdIK6zA71fSMN3fieuQOjyJiqT/3qlGd621rvKqWsRmPAVA85EIIYT8S0WJj5mDnTH7QDy2XkzBJG9b6KrTbgytqVkb3EZGRmLixInw8fFBdnY2AGDXrl24dOmS1OfatGkT7O3toa6uDk9PT0RGRr60f0REBDw9PaGurg4HBwfRsgP1Dh8+DC8vL+jr60NLSwtubm7YtWtXi1+3rVxPe4IaAUMnAw3YGGlyHQ4hhBA5MqqXFZxMtVH8rAbbI1O5DqfDk7pIOnToEIYOHQoNDQ3ExMSgqqoKAFBaWoqVK1dKda4DBw5g1qxZWLRoEWJjY+Hn54fhw4cjIyOj0f6pqakYMWIE/Pz8EBsbi4ULF2LGjBk4dOiQqI+hoSEWLVqEqKgoJCQkYMqUKZgyZYpoNfDmvG5bKq2shYWeOo0iEUIIaUCJz8PswZ0BAL9eSkVRRTXHEXVsPMaYVLvmubu7Y/bs2fjggw+go6OD+Ph4ODg4IC4uDsOGDUNubq7E5+rbty88PDwQGhoqauvatSsCAwMREhLSoP+8efNw9OhRJCYmitqCgoIQHx+PqKioJl/Hw8MDI0eOxLJly5r1ugBQVVUlKgiBurv8rK2tUVxcDF1d2d6BxhhDZY0QGqp0iychhBBxQiHDiHWRuJdbimkDHfH1MBeuQ2pXSkpKoKenJ9H3t9QjSUlJSejfv3+Ddl1dXRQVFUl8nurqakRHRyMgIECsPSAgAFeuXGn0mKioqAb9hw4dips3b6KmpqZBf8YYzp07JxZzc14XAEJCQqCnpyd6WFtbS5Rnc/B4PCqQCCGENIrP52FOQBcAwG+X01BQVvWKI0hzSV0kWVhY4MGDBw3aL126BAcHB4nPU1BQAIFAADMzM7F2MzOzJkejcnNzG+1fW1uLgoICUVtxcTG0tbWhqqqKkSNHYv369RgyZEizXxcAFixYgOLiYtEjM5NWPSWEEMKNwV1N0auTHp7VCLA5/CHX4XRYUhdJn332GWbOnIlr166Bx+Ph0aNH2LNnD+bOnYtp06ZJHQCPJ74OEGOsQdur+r/YrqOjg7i4ONy4cQMrVqxAcHAwwsPDW/S6ampq0NXVFXsQQgghXODxeAh+Ppq062o6HpdUchxRxyT1tiRff/01iouL4e/vj8rKSvTv3x9qamqYO3cuvvjiC4nPY2xsDCUlpQajN3l5eQ1GeeqZm5s32l9ZWRlGRkaiNj6fDyenuvUj3NzckJiYiJCQEAwcOLBZr0sIIYTIm/7OxvCyNcDN9KfYeOEBlr7lynVIHU6zlgBYsWIFCgoKcP36dVy9ehX5+fmiSdGSUlVVhaenJ8LCwsTaw8LC4OPj0+gx3t7eDfqfOXMGXl5eUFFpeq0Ixpho0nVzXpcQQgiRNzzev3OT9l3PQNbTCo4j6niatcEtAGhqasLLy6tFLx4cHIxJkybBy8sL3t7e2Lp1KzIyMhAUFASgbh5QdnY2du7cCaDuTrYNGzYgODgYU6dORVRUFLZv3459+/aJzhkSEgIvLy84OjqiuroaJ0+exM6dO8XuZHvV6xJCCCHtgbejEXwcjXDlYSE2nH+A78b25DqkDkXiImnMmDES9Tt8+LDELz5u3DgUFhZi6dKlyMnJgaurK06ePAlbW1sAQE5OjtjaRfb29jh58iRmz56NjRs3wtLSEuvWrcPYsWNFfcrLyzFt2jRkZWVBQ0MDLi4u2L17N8aNGyfx6xJCCCHtxZyAzrgSGoU/o7MQNMARdsZaXIfUYUi8TtKUKVPEft67dy/efPNN6OjoiLX/9ttvsotOjkmzzgIhhBDSmib/dh3hSfkY426FH8e5cR2OXJPm+1vqxSTr/XchSUVERRIhhBB5kZBVhFEbLoPHA8Jm94eTqc6rD1JQrbqYJCGEEELkS89O+gjoZgbGgJ/O3uc6nA6DiiRCCCGkA5g9pG5PtxMJOUjMKeE4mo6BiiRCCCGkA+hqoYs3eloAAH4MS5bZeQVCho0XHuBCUp7MztleSHx329GjR8V+FgqFOHfuHG7fvi3WPmrUKNlERgghhBCpzBrcGSdv5SDs7mMkZBWhZyf9Fp/zaHw2Vp1Ogo6aMm58MxjqKoqzt6jEE7f5/FcPOvF4PAgEghYH1R7QxG1CCCHyKPiPOByOycaAzibY8VGfFp2LMYY3N1zC7ey6y3frx7vjzV6WsgiTM60ycVsoFL7yoSgFEiGEECKvZg5yhhKfh4jkfNxMe9Kic91MfyoqkADgr9jslobXrtCcJEIIIaQDsTXSwrtenQAAa860bG7Sr5dSAQDeDnX7o0Yk5yO/tKplAbYjVCQRQgghHcwXrztDVYmPqJRCXHlQ0KxzZD6pwOk7dZvB/99b3dHLWh8CIcPR+EeyDFWuUZFECCGEdDBW+hoY38caALAmLBnNWTd6Z1QahAzwczZGZzMdjPWwAgD8FZsl01jlGRVJhBBCSAc03d8Jasp8RKc/RURyvlTHllXVYv+NTADAR772AIA3elpCmc/D7ewSJOWWyjxeeURFEiGEENIBmeqqY9JrdRu3/yjlaNKh6CyUVtbCwVgLAzqbAAAMtVTh72IKADisIKNJzSqSioqK8Msvv2DBggV48qRu5nxMTAyysxVr1jshhBAiz4IGOkJTVQkJWcUIu/tYomOEQobfLtdN2J7iawc+nyd6rv6S25HYRxAIm7X1a7sidZGUkJCAzp074/vvv8fq1atRVFQEAPjrr7+wYMECWcdHCCGEkGYy1lbDZB87AHWjSUIJCpsLSXlIK6yArroyxnh0EnvO38UUehoqyC2pRNTDwtYIWa5IXSQFBwdj8uTJuH//PtTV1UXtw4cPx8WLF2UaHCGEEEJa5tP+DtBRU8a93FKcvJ3zyv7bn9/2P76PDbTUxDfmUFNWEm19cjim419yk7pIunHjBj777LMG7VZWVsjNzZVJUIQQQgiRDX1NVXzsVzf5+qew5JdeJkvMKcGVh4VQ4vPwwfMRqBfVjy6dup2L8qpamccrT6QuktTV1VFS0nB34aSkJJiYmMgkKEIIIYTIzkf97KGvqYKH+eU4Etf0/OH6uUjDupvDSl+j0T4eNvqwM9LEsxqBaB2ljkrqIumtt97C0qVLUVNTA6Buv7aMjAzMnz8fY8eOlXmAhBBCCGkZXXUVfNrfAQDw87n7qBEIG/QpKKvC33F1C0V+1M+uyXPxeDyMdq8bTToc07Fv2JK6SFq9ejXy8/NhamqKZ8+eYcCAAXBycoKOjg5WrFjRGjESQgghpIUm+9jBWFsV6YUVOBTdcD7R3msZqK4VolcnPXjYGLz0XKPd6+5yu/ywADnFz1olXnkgdZGkq6uLS5cu4dChQ/juu+/wxRdf4OTJk4iIiICWllZrxEgIIYSQFtJUVUbQAEcAwPrzD1BV+++m9FW1Auy6mg6g7tIcj8dr9Bz1bIw00cfOEIwBR+I67jYlUhdJaWlpAIDXX38dc+fOxddff43BgwfLOi5CCCGEyNjE12xhpquG7KJnOPB8RW0AOJGQg/zSKpjpqmG4q4VE5xrzfM2kQ9FZzdr2pD2QukhycHBAv379sGXLFtFCkoQQQgiRf+oqSvjC3wkAsOH8A1TWCMAYE932/4G3HVSVJSsNhvewgKoyH/fzynDnUcMbujoCqYukmzdvwtvbG8uXL4elpSXeeust/Pnnn6iqqmqN+AghhBAiQ+/2toaVvgbySquw+2o6bqQ9xZ1HJVBT5mNCHxuJz6OnoYIh3cwAdNwJ3FIXSR4eHli1ahUyMjJw6tQpmJqa4rPPPoOpqSk++uij1oiREEIIITKipqyEGYPqRpNCwx9i44UHAOrWPzLQUpXqXPXblByNz270jrn2rtkb3PJ4PPj7+2Pbtm04e/YsHBwcsGPHDlnGRgghhJBWMMajE2yNNFFYXo2I5HwAwEe+dlKfx8/ZBEZaqigoq0bk/XwZR8m9ZhdJmZmZ+OGHH+Dm5obevXtDS0sLGzZskGVshBBCCGkFKkp8zBrsLPrZz9kYzmY6zTrPKDdLAMChDnjJTeoiaevWrRgwYADs7e2xY8cOvPvuu3j48CEuXbqEzz//vDViJIQQQoiMjeplBRfzusJoqp9Ds88z9vk2JWF3H6P4WY1MYpMXPCblfXvW1tZ477338P7778PNza2VwpJ/JSUl0NPTQ3FxMXR1dbkOhxBCCJFaQVkVMp5UvHLxyJdhjGHo2otIflyG78b0wHtSTP7mgjTf31KPJGVkZGDVqlUyK5A2bdoEe3t7qKurw9PTE5GRkS/tHxERAU9PT6irq8PBwQGbN28We37btm3w8/ODgYEBDAwMMHjwYFy/fl2sT21tLb755hvY29tDQ0MDDg4OWLp0KYTCjjfpjBBCCGmKsbZaiwokoGNvUyJ1kVS/CmdFRQXu3buHhIQEsYc0Dhw4gFmzZmHRokWIjY2Fn58fhg8fjoyMjEb7p6amYsSIEfDz80NsbCwWLlyIGTNm4NChQ6I+4eHhGD9+PC5cuICoqCjY2NggICAA2dn//h/3/fffY/PmzdiwYQMSExPxww8/YNWqVVi/fr20/xyEEEKIwgt0twSPB1xPe4LMJxVchyMzUl9uy8/Px+TJk/HPP/80+rxAIGi0vTF9+/aFh4cHQkNDRW1du3ZFYGAgQkJCGvSfN28ejh49isTERFFbUFAQ4uPjERUV1WQ8BgYG2LBhAz744AMAwBtvvAEzMzNs375d1G/s2LHQ1NTErl27Gj1PVVWV2FpQJSUlsLa2pstthBBCCICJv1zDpQcFCB7SGTMGOb/6AI606uW2WbNmoaioCFevXoWGhgb++ecf7NixA87Ozjh69KjE56murkZ0dDQCAgLE2gMCAnDlypVGj4mKimrQf+jQobh58yZqahqfLFZRUYGamhoYGhqK2vr164dz584hOTkZABAfH49Lly5hxIgRTcYbEhICPT090cPa2lqiPAkhhBBFUL/p7eGYjrNNidRF0vnz5/HTTz+hd+/e4PP5sLW1xcSJE/HDDz80OvrTlIKCAggEApiZmYm1m5mZITc3t9FjcnNzG+1fW1uLgoKCRo+ZP38+rKysxPaXmzdvHsaPHw8XFxeoqKjA3d0ds2bNwvjx45uMd8GCBSguLhY9MjMzm+xLCCGEKJphrubQUFFCWmEFYjKKuA5HJqQuksrLy2FqagoAMDQ0RH5+3eJRPXr0QExMjNQBvLjTMGPspbsPN9a/sXYA+OGHH7Bv3z4cPnwY6urqovYDBw5g9+7d2Lt3L2JiYrBjxw6sXr36pYthqqmpQVdXV+xBCCGEkDpaasoY7moOAPgrNovjaGRD6iKpS5cuSEpKAgC4ublhy5YtyM7OxubNm2FhIdnOwQBgbGwMJSWlBqNGeXl5DUaL6pmbmzfaX1lZGUZGRmLtq1evxsqVK3HmzBn07NlT7LmvvvoK8+fPx3vvvYcePXpg0qRJmD17tlQjYYQQQggRN+b5mknH4nNQVSv5HGV51aw5STk5OQCAJUuW4J9//oGNjQ3WrVuHlStXSnweVVVVeHp6IiwsTKw9LCwMPj4+jR7j7e3doP+ZM2fg5eUFFRUVUduqVauwbNky/PPPP/Dy8mpwnoqKCvD54qkrKSnREgCEEEJIC3g7GsFMVw3Fz2pw4V4e1+G0HJPQ/fv3G20vLy9n0dHRLD8/X9JTiezfv5+pqKiw7du3s7t377JZs2YxLS0tlpaWxhhjbP78+WzSpEmi/ikpKUxTU5PNnj2b3b17l23fvp2pqKiwgwcPivp8//33TFVVlR08eJDl5OSIHqWlpaI+H374IbOysmLHjx9nqamp7PDhw8zY2Jh9/fXXEsdeXFzMALDi4mKp8yaEEEI6qpUn7zLbecfZ1B03uA6lUdJ8f0tcJPF4PNapUyc2adIk9ttvv7HU1NSWxCiyceNGZmtry1RVVZmHhweLiIgQPffhhx+yAQMGiPUPDw9n7u7uTFVVldnZ2bHQ0FCx521tbRmABo8lS5aI+pSUlLCZM2cyGxsbpq6uzhwcHNiiRYtYVVWVxHFTkUQIIYQ0dC+nhNnOO86cFp5ghWWSf6+2FWm+vyVeJykyMhIREREIDw9HVFQUKisrYWNjg9dffx3+/v7w9/eHlZVVK413yR/aloQQQghp3Mh1kbjzqARL3+qOD7ztuA5HjDTf31IvJgkANTU1iIqKQnh4OMLDw3H16lVUVVXByclJNKm7o6MiiRBCCGnc9kupWHb8Ltys9fH3dF+uwxHTqotJAoCKigr69++Pr776CgsWLMC0adOgra2NBw8eNCtgQgghhHQco3pZQonPQ1xmER7ml3EdTrNJVSRVVlbi/PnzWLx4sWgT2RkzZqCsrAyhoaFN7rlGCCGEEMVhoqOGAZ1NAAB/teNNb5Ul7ThgwADcuHEDjo6O6N+/P7788ksMGDCgyTWNCCGEEKK4Rrtb4fy9PPwVm43gIZ3B5ze9ULS8kngk6cqVKzA2Noa/vz8GDRqE119/nQokQgghhDRqSDcz6KgpI7voGa6nPeE6nGaRuEgqKirC1q1boampie+//x5WVlbo0aMHvvjiCxw8eFC0PQkhhBBCiLqKEkb2rNuJ43BM+9ympFl3twFAaWkpLl26hAsXLiA8PBzx8fFwdnbG7du3ZR2jXKK72wghhJCXu576BO9uiYK2mjJuLBoMDVUlrkNq/bvbAEBLSwuGhoYwNDSEgYEBlJWVkZiY2NzTEUIIIaSD8bI1QCcDDZRV1SIs8THX4UhN4iJJKBTi+vXr+OGHHzB8+HDo6+vDx8cHmzZtgrm5OTZu3IiUlJTWjJUQQggh7Qifz8MY97qFptvjJTeJ727T19dHeXk5LCwsMHDgQPz444/w9/eHo6Nja8ZHCCGEkHZstEcnrDv/ABeT85FXWglTHXWuQ5KYxEXSqlWr4O/vj86dO7dmPIQQQgjpQOyNteBuo4/YjCIcjXuET/wcuA5JYhJfbvvss8+oQCKEEEKI1MZ4dAIAHG5nC0s2e+I2IYQQQogk3uxpARUlHu7mlCAxp4TrcCRGRRIhhBBCWpW+pioGudQtQP1XbPsZTaIiiRBCCCGtbrRH3V1uf8dmQyBs1hKNbY6KJEIIIYS0Ov8uptDXVEFeaRUuPyjgOhyJyKxIysnJQUZGhqxORwghhJAORFWZj1G9LAEAB25molYg5DiiV5N4CYBXef3115GcnAyBQCCrUxJCCCGkAxntboWdUek4kZCDsLuP4WyqjS7mOnAx14GLuS5czHVgoqMGHo/HdagAZFgk7dy5ExUVFbI6HSGEEEI6GDdrfYz16IRTt3NQUS3AnUcluPNI/G43A02V54WTLl5zMMQwVwuOom3BBreKjja4JYQQQppHKGTIfFqBe7mlSHr+SMwtQVpBOf47p3tED3Nset9Tpq8tzfd3s0eSqqurkZeXB6FQ/JqijY1Nc09JCCGEEAXA5/Nga6QFWyMtDO1uLmqvrBHgQV4Z7uWW4l5OCVyt9DiMshlF0v379/HRRx/hypUrYu2MMfB4PJqTRAghhJBmUVdRgquVHufFUT2pi6TJkydDWVkZx48fh4WFhdxMriKEEEIIkSWpi6S4uDhER0fDxcWlNeIhhBBCCJELUq+T1K1bNxQUtI9FoAghhBBCmkvqIun777/H119/jfDwcBQWFqKkpETsQQghhBDSEUi9BACfX1dXvTgXSdEmbtMSAIQQQkj7I833t9QjSRcuXMCFCxdw/vx5sUd9m7Q2bdoEe3t7qKurw9PTE5GRkS/tHxERAU9PT6irq8PBwQGbN28We37btm3w8/ODgYEBDAwMMHjwYFy/fr3BebKzszFx4kQYGRlBU1MTbm5uiI6Oljp+QgghhHRMUk/cHjBggMxe/MCBA5g1axY2bdoEX19fbNmyBcOHD8fdu3cbXW8pNTUVI0aMwNSpU7F7925cvnwZ06ZNg4mJCcaOHQsACA8Px/jx4+Hj4wN1dXX88MMPCAgIwJ07d2BlVbcD8dOnT+Hr6wt/f3+cOnUKpqamePjwIfT19WWWGyGEEELaN6kvtyUkJDR+Ih4P6urqsLGxgZqamkTn6tu3Lzw8PBAaGipq69q1KwIDAxESEtKg/7x583D06FEkJiaK2oKCghAfH4+oqKhGX0MgEMDAwAAbNmzABx98AACYP38+Ll++/MpRq5ehy22EEEJI+9OqK267ubm9dG0kFRUVjBs3Dlu2bIG6unqT/aqrqxEdHY358+eLtQcEBDRYqLJeVFQUAgICxNqGDh2K7du3o6amBioqKg2OqaioQE1NDQwNDUVtR48exdChQ/HOO+8gIiICVlZWmDZtGqZOndpkvFVVVaiqqhL9TJPUCSGEkI5N6jlJf/31F5ydnbF161bExcUhNjYWW7duRZcuXbB3715s374d58+fxzfffPPS8xQUFEAgEMDMzEys3czMDLm5uY0ek5ub22j/2traJpclmD9/PqysrDB48GBRW0pKCkJDQ+Hs7IzTp08jKCgIM2bMwM6dO5uMNyQkBHp6eqKHtbX1S/MjhBBCSPsm9UjSihUr8PPPP2Po0KGitp49e6JTp05YvHgxrl+/Di0tLcyZMwerV69+5fmauktOmv6NtQPADz/8gH379iE8PFxsVEsoFMLLywsrV64EALi7u+POnTsIDQ0VXZJ70YIFCxAcHCz6uaSkhAolQgghpAOTuki6desWbG1tG7Tb2tri1q1bAOouyeXk5Lz0PMbGxlBSUmowapSXl9dgtKieubl5o/2VlZVhZGQk1r569WqsXLkSZ8+eRc+ePcWes7CwQLdu3cTaunbtikOHDjUZr5qamthcq/rijC67EUIIIe1H/fe2JFOypS6SXFxc8N1332Hr1q1QVVUFANTU1OC7774TbVWSnZ3dZKFTT1VVFZ6enggLC8Po0aNF7WFhYXjrrbcaPcbb2xvHjh0Taztz5gy8vLzE5iOtWrUKy5cvx+nTp+Hl5dXgPL6+vkhKShJrS05ObrT4a0ppaSkA0GgSIYQQ0g6VlpZCT+/lG+lKXSRt3LgRo0aNQqdOndCzZ0/weDwkJCRAIBDg+PHjAOrm/EybNu2V5woODsakSZPg5eUFb29vbN26FRkZGQgKCgJQd4krOztbNFcoKCgIGzZsQHBwMKZOnYqoqChs374d+/btE53zhx9+wOLFi7F3717Y2dmJRp60tbWhra0NAJg9ezZ8fHywcuVKvPvuu7h+/Tq2bt2KrVu3SvzvYGlpiczMTOjo6Ei8yW/9JbrMzMwOfUecIuSpCDkCipGnIuQIKEaeipAjQHm2FGMMpaWlsLS0lKiz1EpLS1loaCibPXs2mzVrFtu8eTMrKSlpzqnYxo0bma2tLVNVVWUeHh4sIiJC9NyHH37IBgwYINY/PDycubu7M1VVVWZnZ8dCQ0PFnre1tWUAGjyWLFki1u/YsWPM1dWVqampMRcXF7Z169ZmxS+N4uJiBoAVFxe3+mtxSRHyVIQcGVOMPBUhR8YUI09FyJExyrMtSb1OEmk+RVlbSRHyVIQcAcXIUxFyBBQjT0XIEaA825JEl9uOHj2K4cOHQ0VFBUePHn1p31GjRskkMEIIIYQQLklUJAUGBiI3NxempqYIDAxssp8ibXDbHGpqaliyZInEK5K3V4qQpyLkCChGnoqQI6AYeSpCjgDl2ZbochshhBBCSCOkXnGbEEIIIUQRSFwkXbt2DadOnRJr27lzJ+zt7WFqaopPP/1UbG8zQgghhJD2TOIi6X//+x8SEhJEP9+6dQsff/wxBg8ejPnz5+PYsWMICQlplSAJIYQQQtqaxHOSLCwscOzYMdEK1osWLUJERAQuXboEAPjzzz+xZMkS3L17t/WiJYQQQghpIxKPJD19+lRsq5GIiAgMGzZM9HPv3r2RmZkp2+gIIYQQQjgicZFkZmaG1NRUAEB1dTViYmLg7e0ter60tFRs/zTScSnCDZGKkCPpWBThd1YRciTyReIiadiwYZg/fz4iIyOxYMECaGpqws/PT/R8QkICHB0dWyVIeVdTU4Ps7GzRzx31jVxZWYmFCxdKtcddeyQQCEQbGHdUVVVVuHLlCtLT07kOpVUpwnuT3pcdB70v5Y/ERdLy5cuhpKSEAQMGYNu2bdi2bRtUVVVFz//6668ICAholSDl2Zo1a+Ds7IyRI0fijTfeQFRUlMQb3rYn27dvh7m5Oa5fvw4VFRU8e/aM65BaxU8//QRfX18EBgZi1qxZePjwIQBAKBRyHJnsrF27FnZ2dvjss8/Qq1cvbN68uUMuAqsI7016X9L7sr1pd+9LaTd7KyoqYrW1tQ3aCwsLWVVVVXP3kGuX1qxZw+zs7NjBgwfZr7/+yt566y1mbGzMLly4wHVoMvXgwQPm5+fHNm/ezHUorSY5OZm9/vrrzNnZmf3+++9s/vz57LXXXmNDhgzhOjSZ+uabb1iXLl3YiRMn2P3799nixYuZvr4+q6io4Do0mVKE9ya9LzsOel9e4Dq0JkldJBHGBAIBq6mpYcOHD2eff/652HP9+vVjw4YNY7GxsdwE1wrWrFnD3N3dGWOMpaens8WLF7NffvmFRUZGchyZbAgEArZmzRo2YsQIlp2dLWr/888/Wa9evVhiYiKH0clOQUEB69u3L1u9erWoLTk5mXXr1o3l5+czxhgTCoVchScTivTepPclvS/bi/b8vqQVt5uBz+dDKBTi1q1b8PDwAFA3LwCoGzJNTk7GP//80+4X12TPrxM/ePAAQ4YMwalTp+Dl5YUbN25g8+bNGDRoEEJDQ9v9ED9jDM7Ozpg+fTosLS1FQ/ja2trIzc2FkZERxxHKhpaWFhISEsT2QVq0aBEsLCzw559/4tGjRxxGJxuK8N6k9yW9L9ub9vy+pCJJAjt37sTs2bOxc+dO5OXlAQBUVVXx2muvYffu3QAAdXV1CIVCeHp6IiAgAIcOHUJBQQGXYUvtxTzrrxPr6+vjjz/+wMmTJ7F8+XIcP34cN27cwJdffokdO3YgPDyc28Cl9M8//4i+aBhjUFJSwsiRIzFixAgAEOVdUlICS0tLaGhocBZrc/03R6Bu0qu6ujrmzp2Lb7/9FoGBgdDT08P9+/fRs2dPrFu3DmPGjMHx48c5jFp6ivDepPclvS/pfckhroaw2oPc3Fw2aNAgZmVlxUaPHs1sbGxYly5d2NWrVxljjO3evZtZWVmxkydPMsYYe/bsGWOMsczMTMbj8UT95F1jebq4uLArV64wxhiLjo5m5ubmjM/ns8uXL4uOKy4uZs7Ozuznn3/mKnSp3L17lw0YMIDxeDy2bNkyxljdMPCL6oe2Z86cySZNmtRkP3kkSY7p6els6tSp7L333hPNLxQKhczDw4MtXLiwXeSqCO9Nel+Ko/el/OfaEd+XNJL0EpGRkcjJyUFMTAwOHz6M5ORk6OrqYtmyZbh9+zaGDBmC3r17Y+XKlQDqKmPGGNTU1GBtbd1uVh9vLE8dHR2sXLkSsbGx6NGjB0aPHg01NTXw+XW/MkKhELq6ujA0NGwXeaalpeH777+HiYkJZsyYgR9++AF5eXng8/kNbj+t/4v19OnTGDhwIIC64eK0tDQA8nu7qqQ5Ghoa4saNG/jwww+hpKSEyspK8Hg8GBoaIi4uTvT/sTxThPcmvS/pfUnvS+7J/786RxhjiIyMhImJCXR0dCAUCqGmpoa1a9ciLy8Pv/zyC0xMTDB16lSkpaVh7ty5qKmpAY/Hw61bt6CmpiZ6I8uzl+X5+PFj7NixA3w+H9OnT4eVlRWWLVuG+/fvg8/nIzExEUKhEOPHj+c6jVcyNTWFm5sbvvrqK8yfPx/29vaYPXt2k/0TEhJQVFSEYcOG4cmTJ/j444/h4OCA5ORkub1dVdIc1dXVUVJSgujoaNHPycnJKC8vx7hx49o6bKkpwnuT3peNo/el/Oqw70suhq/kXf2w7rx581jXrl0ZY0xs2YOFCxey1157TXQXyf79+5mmpibz8PBgkydPZrq6uuzTTz9lz549k+u7EiTNMzw8nDHG2JUrV5i1tTWzsbFh7777LjMyMmLvvfceKykpafvgpVCfZ2Vlpajtjz/+EBvefXEoe//+/axPnz4sJCSE6erqsv79+7O7d++2XdBSkibH6upqtmrVKsbj8di4cePYjBkzmKmpKRs7dix7+vRpm8cuDUV4b9L7kt6X9L6UH1QkNaL+/6To6GimoqLCwsLCGGP//qKnpqYye3t7tnHjRtEx586dY6tXr2ZTpkxhx48fb/ugm0HSPDdt2iTqe/v2bbZnzx729ddfs1OnTnETeAvU51FSUsLeeOMN1qdPn0b7ffjhh4zH4zFnZ2d26NChtgyxxSTN8ccff2SffvopGz16tGiOgLxo6oOyI703W5pje3lfSvKl197fl7LMkd6X8kVhi6T6/1MbWxizXkFBARszZgzr0aOHqK2+/8iRI9m4ceNaN0gZoDybduXKFaampsZ27twpOrb+r+9ff/2Vbd26VfaBtoCsciwuLm6V+GSlurqaVVdXv7RPe/+dVYQcGZMszxe1t/elrHKU9/dlZWWlaN0mxhovmDrC7+yLFHJO0tq1azFlyhQAgJKSUpP9jIyMMH36dGRlZWHFihWi/tXV1aisrISNjU2bxNtclOfLeXl54fPPP8eiRYuQmJiIDz74AN9//z2qq6sxZcoUTJ06tbVClposc1y9ejXKy8tbK9QW+e677xAQEIDRo0fj559/Ft0S/OL2E+35d1YRcgQkz/NF7el9Kcsc5fl9uWLFCvTu3RtvvfUWJkyYgJSUFPB4vAbbprT339lGcV2ltaXExET21ltvMS0tLWZqasr27dvHGHv5X+ZVVVVs06ZNjMfjsTlz5rALFy6w9evXM0tLS3bu3Lm2Cl0qlKfkIy03b95kPB6P8Xg85uLiInfzGxQhR8YYi42NZX369GHdunVjv/zyCxs/fjzz8PBgX375ZZPHtLffWUXIkbHm5fkief+dVYQcGasb8fLw8GBdu3Zlf/zxB1u1ahXz8fFh77zzTpPHtMff2ZdRqCLp999/Z2+++SY7cOAAmzJlCvP29hbtjfOqa8qrVq1i/fr1Y127dmU2Njbszz//bIuQm4XyfHWeQqGQnTt3jpmbmzMbGxu5vSauCDmWl5ezefPmsfHjx4tNUF2yZAkbMWIEKywsfOnx7eF3VhFyZKzlebaH31lFyLHesmXL2Keffip2KXDatGlsypQpop+b+hxqL7+zr6IQRVL93QNlZWUsIiKCMcbY0aNHmbu7O1u+fLlYn6aOZazul0Eeq/16lOer86xXU1PDQkJC2P/+97/WDbaZFCHHesXFxWzTpk3s4sWLjLG6uBlj7Oeff2ZOTk6in1/Unn5nFSFHxpqfZ7328DurCDnWFz65ubns/v37ovba2lo2ePBgtnz5cnbv3r1Gj21vv7OvwmNMTlfgaqHff/8dmZmZ6NOnD/r16wctLS2x558+fYqQkBCcPHkSR44cgaOjIwQCQZPzPRhjcrkOB+VZR5o863MUCoVytUCbIuQIiOfp6+sLbW1t1NbWQllZGQBEMS9btgxxcXE4ePDgS38n5fF3VhFyBGSbp7z+zipCjsCrP3/279+PoKAgWFhYQFtbG4WFhRg1ahTWrl3b6Pnk9XdWahwVZ63m4cOHzNPTk9nb27OBAwcyc3NzFhAQwAoKCkR96ivdyMhI1r9/f/bJJ59wFW6zUZ4dJ09FyJGxxvMcMmSIWJ5CoVD0V+xbb70lGjWTt7VTmqIIOTKmGHkqQo6MSfb5wxhjW7ZsYX/++SerqqpiT548YcePH2c8Ho/dvn2bMda+cpZGhyuS1q5dy/r27csqKipYRUUFS0xMZEZGRuyTTz5h6enpjLF/J73W1NSwVatWsS5durALFy4wxhi7evVqu/g/m/LsOHkqQo6MvTzPjIwMxti/xWBJSQkzMTERm+iZmZnJGJPvD2NFyJExxchTEXJkTLLPH8Ya5pGSksKMjIzY9u3b2zrkNtWhiqTq6mo2aNAg0V/Z9V8sf/75J7Ozs2OhoaGivvW/3Hfu3GFvv/028/X1ZW+88Qbj8Xjs1q1bbR+8FCjPjpOnIuTImHR5MlY3/8re3p4JhUL26NEjNn78eGZkZMRycnLaPHZJKUKOjClGnoqQI2PS5/lfO3bsYH379mW5ubltEitX5OuiaAswxqCiogI9PT08evRI1AYAb7/9Njw8PHD48GE8fPgQAETXg01MTJCXl4crV65AWVkZKSkpcHV15SYJCVCeHSdPRcgRkD5PALh37x569eqF7777Dk5OTsjLy0N0dDTMzc05yeFVFCFHQDHyVIQcgeblWVBQgEePHmHp0qVYsGABRo8eDWNjY7ndXFgm2rgoa3W7du1iRkZGouuk9cuiR0dHMzU1NXbp0iVR31u3brEuXbowBwcH0R1E7QXl2XHyVIQcGZMuT19fX9H6MfK2zcbLKEKOjClGnoqQI2OS53n37l22fPly1rlzZ9atWzd29uxZzmJuSx2mSKpfFj46Opr5+/uLLX9ef5nC1dWVffvtt6L2yspK9tdff7VpnC1FeXacPBUhR8akz7OsrIx98803bMeOHW0fbDMpQo6MKUaeipAjY5LnuXjxYsYYYxUVFezcuXPsyJEjbR8sh9pNkXTo0KFGd0JubOXhLVu2MAsLC7Z//35RW05ODuvSpQvbtm1bk8fJA8qz4+SpCDkyJvs85ZEi5MiYYuSpCDkypjh5tja5L5IiIiJYt27dGI/HE9tB+EUHDhxgKioq7NChQ+zJkyds9uzZTFNTk23atInFx8ez77//njk5ObGEhIQ2jF5ylKe49pynIuTImGLkqQg5MqYYeSpCjowpTp5tRa6LpKSkJDZhwgT2+eefs2nTpjFLS0uWlZUlel4oFLLi4mL23nvvMRMTE7Z69WpWVVUlej4oKIh1796dOTo6MhsbG3bixAku0nglyrNOR8hTEXJkTDHyVIQcGVOMPBUhR8YUJ8+2JNdFUl5eHvvll1/Y7du3WXl5ObOxsWEzZ84U61NVVcX2798vdhti/XoOQqGQlZWVsZs3b7Zl2FKjPP/V3vNUhBwZU4w8FSFHxhQjT0XIkTHFybMtyVWRFBkZKbZ4FWNMbB+cnTt3MlVVVRYdHc0Ye/UiXa/a24orlGedjpCnIuTImGLkqQg5MqYYeSpCjowpTp5ckosi6dy5c8ze3p7Z2toyCwsLNmnSJBYTE8MYa/h/qq+vLxsxYoRoZn57Qnl2nDwVIUfGFCNPRciRMcXIUxFyZExx8pQHnBdJmZmZzNvbmy1atIilp6ezY8eOMTc3NzZo0CDR7sP/nY1/5coVxufz2eHDhxljdZXvi3vMyCPKs+PkqQg5MqYYeSpCjowpRp6KkCNjipOnvOC8SDpz5gxTV1dnycnJorbTp083WLfhvyZPnsx69uzJzp49y4YNG8YWLFggWgBLXlGeHSdPRciRMcXIUxFyZEwx8lSEHBlTnDzlBedF0v79+5mbmxu7d++eqE0gELCNGzcyOzs7dvr0acaYeGV87do1xuPxGI/HY0OGDGGFhYVtHre0KM+Ok6ci5MiYYuSpCDkyphh5KkKOjClOnvKC8yLp1q1bTE1Njf39999i7ffu3WOBgYFs8uTJorba2lq2Z88epqamxry8vNiNGzf+v727B2lrDcA4/hyOZoiCEZGUlGpAalAELVRcWmg7OTmoECQgdS10UnRXEJz8WDImxUEnuxTpVjuqKIIoIviBOKi0IlQrRsx7pysNPcu9VU963v9vM4TD+98ezjHJQx/3f6MzOJ02NBpjR6cNjcbY0WlDozH2dBYL33/gtqmpSW/evNH4+LjOz89vX08kEqqpqdHR0ZEuLy8lSblcTt+/f9fExISWl5f1/Plzv479n9EZnE4bGiU7Om1olOzotKFRsqezaPi90owxZm1tzZSUlJh0Ol3wnHRkZMTU1NQUfKTxb0ZncDptaDTGjk4bGo2xo9OGRmPs6SwGJX6PNElqbm7W0NCQhoeH5bquenp6lM/ntbS0pFQqpZKSojjmH6MzOJ02NEp2dNrQKNnRaUOjZE9nUfB7pf3q3bt3JhqNmtbWVlNbW2saGhrM+vq638e6c3QGhw2NxtjRaUOjMXZ02tBojD2dfnKMMcbvofavq6srbW5uam1tTaFQSKlUyu8j3Qs6g8OGRsmOThsaJTs6bWiU7On0U1GNJAAAgGLh+6fbAAAAihEjCQAAwAMjCQAAwAMjCQAAwAMjCQAAwAMjCQAAwAMjCQAAwAMjCQAAwAMjCYBVFhYW5DiOzs7O/D4KgCLHN24DCLRXr16ppaVFExMTkqRcLqfT01NFo1E5juPv4QAUNX4qGIBVQqGQHj165PcxAPwFeNwGILDevn2rr1+/anJyUo7jyHEcZbPZgsdt2WxWkUhEnz59UiKRUDgcVnd3ty4uLvThwwfF43FVVlbq/fv3urm5ub12LpfT4OCgHj9+rLKyMrW1tWlhYcGfUAD3gjtJAAJrcnJS29vbampq0vDwsCRpY2Pjt/f9/PlTU1NTmp2d1Y8fP9TZ2anOzk5FIhHNz89rd3dXXV1devHihZLJpCSpr69P+/v7mp2dVSwW08ePH9Xe3q719XU9ffr0QTsB3A9GEoDAqqioUCgUUjgcvn3EtrW19dv7rq+vlU6nVVdXJ0nq7u7W9PS0jo+PVV5ersbGRr1+/VpfvnxRMpnUzs6OZmZmdHh4qFgsJkkaGBjQ58+flclkNDo6+nCRAO4NIwmA9cLh8O1AkqRoNKp4PK7y8vKC105OTiRJq6urMsaovr6+4DpXV1eqqqp6mEMDuHeMJADWKy0tLfjbcRzP1/L5vCQpn8/LdV2trKzIdd2C9/06rAD83RhJAAItFAoV/MP1XXj27Jlubm50cnKily9f3um1ARQPPt0GINDi8bgWFxe1v7+vb9++3d4N+hP19fVKpVLq7e3V3Nyc9vb2tLy8rLGxMc3Pz9/BqQEUA0YSgEAbGBiQ67pqbGxUdXW1Dg4O7uS6mUxGvb296u/vVyKRUEdHhxYXF/XkyZM7uT4A//GN2wAAAB64kwQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAOCBkQQAAODhHwatDfFEjjRrAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "diff = (ts[\"Sign. Wave Height\"] - dsi[\"Sign. Wave Height\"].sel(x=250, y=1200))\n", "diff.plot()" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.03509921674697619" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "mae = np.abs(diff).nanmean().to_numpy()\n", "mae" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Clean up" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "import os\n", "os.remove(\"waves_3h.dfs2\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.6" } }, "nbformat": 4, "nbformat_minor": 4 }