{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Metocean gridded spatial skill\n", "\n", "Demonstrate plotting functionality etc by the SpatialSkill class" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import modelskill as ms" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\jem\\Source\\modelskill\\modelskill\\timeseries\\_track.py:134: UserWarning: Removed 22 duplicate timestamps with keep=first\n", " warnings.warn(\n" ] } ], "source": [ "fn = '../tests/testdata/NorthSeaHD_and_windspeed.dfsu'\n", "mr = ms.model_result(fn, name='HD', item=0)\n", "fn = '../tests/testdata/altimetry_NorthSea_20171027.csv'\n", "df = pd.read_csv(fn, index_col=0, parse_dates=True)\n", "o1 = ms.TrackObservation(df, item=2, name='alti')\n", "cc = ms.match(o1, mr)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ss = cc.gridded_skill()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
                            "Dimensions:      (y: 5, x: 5)\n",
                            "Coordinates:\n",
                            "  * y            (y) float64 50.6 51.66 52.7 53.75 54.8\n",
                            "  * x            (x) float64 -0.436 1.543 3.517 5.492 7.466\n",
                            "    observation  <U4 'alti'\n",
                            "Data variables:\n",
                            "    n            (y, x) int32 3 17 0 0 0 0 49 9 0 ... 70 33 15 27 37 79 0 20 75\n",
                            "    bias         (y, x) float64 -0.02626 -0.188 nan nan ... nan -0.2475 -0.1156\n",
                            "    rmse         (y, x) float64 0.03291 0.2105 nan nan ... nan 0.2635 0.1475\n",
                            "    urmse        (y, x) float64 0.01984 0.09468 nan nan ... nan 0.09066 0.09174\n",
                            "    mae          (y, x) float64 0.02626 0.1923 nan nan ... nan 0.2475 0.1217\n",
                            "    cc           (y, x) float64 0.9973 0.996 nan nan ... nan 0.6293 0.9842\n",
                            "    si           (y, x) float64 0.04669 0.08355 nan nan ... nan 0.0844 0.1475\n",
                            "    r2           (y, x) float64 -1.061 0.9003 nan nan ... nan -18.78 0.9078
" ], "text/plain": [ "\n", "Dimensions: (y: 5, x: 5)\n", "Coordinates:\n", " * y (y) float64 50.6 51.66 52.7 53.75 54.8\n", " * x (x) float64 -0.436 1.543 3.517 5.492 7.466\n", " observation " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ss.rmse.plot();" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUMAAAE6CAYAAACMHjAOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABUjklEQVR4nO3dd1QU19sH8O/Qlo5YaAoIKBoUe0SxYDf2ktiVYoktKqJiSwQVwYqoURONEUhiNIlgokYRC9ijWKIoIioCKsiLDUX63vcPwvxYdmF3VpYiz+ecOYeduXvnYYhP7sy9cy/HGGMghJBaTq2qAyCEkOqAkiEhhICSISGEAKBkSAghACgZEkIIAEqGhBACgJIhIYQAoGRICCEAKBkSQggASoYqExwcDI7jEBMTU9WhqMy+ffsQFBQk8xjHcfD19a3UeKpS6d/37t278PX1xePHj6XKuru7o3HjxpUWG1EMJUOitPKS4aVLlzB16tTKDagauXv3LlauXCkzGX7zzTcIDw+v/KBIuTSqOgCieu/fv4eurm6lnrNTp06Ver6axM7OrqpDIDJQy1AJ58+fR+/evWFgYABdXV04Ozvj6NGjMsu+evUKHh4eqFu3LvT09DBkyBA8evRIosyNGzcwePBgmJiYQCQSwcLCAoMGDcKTJ0/4Mowx7NixA23atIGOjg6MjY3xxRdfSNXVo0cPtGzZEmfPnoWzszN0dXUxefJkDB8+HNbW1hCLxVIxOjk5oV27dvzn7du3o3v37jAxMYGenh4cHR2xfv165OfnS5zn6NGjSEpKAsdx/FZM1m1ybGwshg0bBmNjY2hra6NNmzYICQmRKBMVFQWO4/Drr79i+fLlsLCwgKGhIfr06YP4+Pgy/iKq8X//93+YNWsWHBwcoK+vDxMTE/Tq1Qvnzp0r93vBwcEYNWoUAKBnz578tQkODgZAt8nVFSVDgaKjo9GrVy+8efMGe/bswa+//goDAwMMGTIEBw4ckCo/ZcoUqKmp8beUV65cQY8ePfD69WsAQFZWFvr27Yvnz59j+/btiIyMRFBQEKysrPD27Vu+nunTp8PT0xN9+vTBoUOHsGPHDty5cwfOzs54/vy5xDlTU1MxceJEjB8/Hn///TdmzZqFyZMnIzk5GadPn5Yoe+/ePVy5cgUeHh78vocPH2L8+PH46aefcOTIEUyZMgUbNmzA9OnT+TI7duxAly5dYGZmhkuXLvFbWeLj4+Hs7Iw7d+5g69atCAsLg4ODA9zd3bF+/Xqp8suWLUNSUhJ++OEH7Nq1CwkJCRgyZAgKCwvL/fswxlBQUKDQJs/Lly8BAD4+Pjh69Cj27t0LW1tb9OjRA1FRUWV+b9CgQfD39wdQ9D+W4mszaNAgueckVYgRQTp16sRMTEzY27dv+X0FBQWsZcuWrFGjRkwsFjPGGNu7dy8DwEaMGCHx/QsXLjAAzM/PjzHGWExMDAPADh06VOY5L126xACwTZs2SexPSUlhOjo6zNvbm9/n4uLCALBTp05JlM3Pz2empqZs/PjxEvu9vb2ZlpYWy8jIkHnuwsJClp+fz0JDQ5m6ujp7+fIlf2zQoEHM2tpa5vcAMB8fH/7z2LFjmUgkYsnJyRLlBgwYwHR1ddnr168ZY4ydOXOGAWADBw6UKPfbb78xAOzSpUsyz1es+PuKbImJieXWVVpBQQHLz89nvXv3lvq7lv59f//9dwaAnTlzRqoeNze3Mq8bqTrUMhQgKysL//zzD7744gvo6+vz+9XV1TFp0iQ8efJE6lZuwoQJEp+dnZ1hbW2NM2fOAACaNGkCY2NjLF68GN999x3u3r0rdd4jR46A4zhMnDhRomVjZmaG1q1bS7VSjI2N0atXL4l9GhoamDhxIsLCwvDmzRsAQGFhIX766ScMGzYM9erV48veuHEDQ4cORb169aCurg5NTU24urqisLAQ9+/fF37hAJw+fRq9e/eGpaWlxH53d3e8f/9eqlU5dOhQic+tWrUCACQlJZV7nvbt2+Pq1asKbRYWFnLj/u6779CuXTtoa2tDQ0MDmpqaOHXqFOLi4hT5tUkNQh0oArx69QqMMZibm0sdK/6H9eLFC4n9ZmZmUmXNzMz4ckZGRoiOjsaaNWuwbNkyvHr1Cubm5pg2bRq+/vpraGpq4vnz52CMwdTUVGZctra2Ep9lxQcAkydPxqZNm7B//35Mnz4dERERSE1NlbhFTk5ORrdu3dCsWTNs2bIFjRs3hra2Nq5cuYLZs2cjOzu7nCtUthcvXgi6biWTMwCIRCIAkHt+fX19tGnTRqGYNDTK/88/MDAQCxYswIwZM7B69WrUr18f6urq+OabbygZfoQoGQpgbGwMNTU1pKamSh179uwZAKB+/foS+9PS0qTKpqWloUmTJvxnR0dH7N+/H4wx3Lp1C8HBwVi1ahV0dHSwZMkS1K9fHxzH4dy5c3xSKKn0vpIdGSU5ODigY8eO2Lt3L6ZPn469e/fCwsIC/fr148scOnQIWVlZCAsLg7W1Nb//5s2bMutUVL169QRdN2VFR0ejZ8+eCpVNTEwstyPj559/Ro8ePbBz506J/SWf5ZKPByVDAfT09ODk5ISwsDBs3LgROjo6AACxWIyff/4ZjRo1gr29vcR3fvnlF3z++ef854sXLyIpKUnmGDyO49C6dWts3rwZwcHBuH79OgBg8ODBWLt2LZ4+fYrRo0d/0O/g4eGBmTNn4vz58zh8+DC8vLygrq4uEQMgmWAZY9i9e7dUXSKRSOGWYu/evREeHo5nz55J3J6GhoZCV1e3wobiFN8mK0LebTLHcVL/o7l16xYuXbokdbtfmqItWVJ9UDIUKCAgAH379kXPnj2xcOFCaGlpYceOHYiNjcWvv/4q1SqLiYnB1KlTMWrUKKSkpGD58uVo2LAhZs2aBaDoeeCOHTswfPhw2NragjGGsLAwvH79Gn379gUAdOnSBV9++SU8PDwQExOD7t27Q09PD6mpqTh//jwcHR0xc+ZMheIfN24cvLy8MG7cOOTm5sLd3V3ieN++faGlpYVx48bB29sbOTk52LlzJ169eiVVl6OjI8LCwrBz5060b98eampq6NChg8zz+vj44MiRI+jZsydWrFiBunXr4pdffsHRo0exfv16GBkZKRS/PAYGBmXGINTgwYOxevVq+Pj4wMXFBfHx8Vi1ahVsbGzk9ka3bNkSALBr1y4YGBhAW1sbNjY2Urf/pBqp0u6bGurcuXOsV69eTE9Pj+no6LBOnTqxw4cPS5Qp7k0+ceIEmzRpEqtTpw7T0dFhAwcOZAkJCXy5e/fusXHjxjE7Ozumo6PDjIyMWMeOHVlwcLDUeX/88Ufm5OTEn9fOzo65urqymJgYvoyLiwtr0aJFufGPHz+eAWBdunSRefzw4cOsdevWTFtbmzVs2JAtWrSIHTt2TKp39OXLl+yLL75gderUYRzHsZL/OaFU7ypjjN2+fZsNGTKEGRkZMS0tLda6dWu2d+9eiTLFvcG///67xP7ExEQGQKq8KuXm5rKFCxeyhg0bMm1tbdauXTt26NAhmb3Bsn7foKAgZmNjw9TV1SVip97k6oljjFbHI4QQGlpDCCGgZEgIIQAoGRJCCABKhoQQAoCSISGEAKBkSAghAGrBoGuxWIxnz57BwMCgzNfUCKnpGGN4+/YtLCwsoKamWBsnJycHeXl5CpXV0tKCtrb2h4RY7X30yfDZs2dyX50i5GORkpKCRo0ayS2Xk5MDG2t9pKWXPz9kMTMzMyQmJn7UCfGjT4YGBgYAAOslK6BWjf6Q18f/UNUhyDTC3rGqQ5CSPaRiXq+rSLvXbavqECS8eyeGc8cM/r93efLy8pCWXojEa9YwNCi/JZn5Vgyb9knIy8ujZFiTFd8aq2lrV6tkaGigLr9QFdDgNKs6BCkamtXn71bMQE4CqSpCHwXp6DPo6Jf/Elp+LXlJ7aNPhoSQsokhhvSqONJlagNKhoTUYoWMoVBOy0/e8Y8FJUNCajExGMQoP9nJO/6xoGRISC0mBkMhJUMAlAwJqdWoZfg/lAwJqcXomeH/UDIkpBYT/7fJK1MbVOlgKV9fX3AcJ7GVXFrT3d1d6nhFLRxECAEK/3tmKG+rDaq8ZdiiRQucPHmS/1xypTYA+Oyzz7B3717+s5aWVqXFRsjHLp8VbfLK1AZVngw1NDRkLrReTCQSlXucEKI8MTgUovy3VsRyjn8sqvydooSEBFhYWMDGxgZjx47Fo0ePJI5HRUXBxMQE9vb2mDZtGtLT08utLzc3F5mZmRIbIUQ2MVNsqw2qNBk6OTkhNDQUERER2L17N9LS0uDs7IwXL14AAAYMGIBffvkFp0+fxqZNm3D16lX06tULubm5ZdYZEBAAIyMjfqMZawgpW+F/LUN5W21QpbfJAwYM4H92dHRE586dYWdnh5CQEHh5eWHMmDH88ZYtW6JDhw6wtrbG0aNHMXLkSJl1Ll26FF5eXvznzMxMSoiElEGRZEfJsAro6enB0dERCQkJMo+bm5vD2tq6zONA0TNGkUikqhAJ+aiIGQcxk/PMUM7xj0WVPzMsKTc3F3FxcTA3N5d5/MWLF0hJSSnzOCFEGLpN/p8qTYYLFy5EdHQ0EhMT8c8//+CLL75AZmYm3Nzc8O7dOyxcuBCXLl3C48ePERUVhSFDhqB+/foYMWJEVYZNyEejEGoKbbVBld4mP3nyBOPGjUNGRgYaNGiATp064fLly7C2tkZ2djZu376N0NBQvH79Gubm5ujZsycOHDig8Gy+hJDyMQVuk1ktuU2u0mS4f//+Mo/p6OggIiKiEqMhpPbJY+rQZOW3/PIoGRJCPnZicBDLuQ2uLbPW1I6HAYQQmVTRgdK4cWOpOQU4jsPs2bMBFC1r6uvrCwsLC+jo6KBHjx64c+eOKn49QSgZElKLFTI1hTYhrl69itTUVH6LjIwEAIwaNQoAsH79egQGBuLbb7/F1atXYWZmhr59++Lt27cV/vsJQbfJhNRiRbfJFftucoMGDSQ+r127FnZ2dnBxcQFjDEFBQVi+fDn/4kRISAhMTU2xb98+TJ8+XWadW7duFRQDAHh4eAjqbKVkSEgtJlZg6EzxM8PS7/kr8oJDXl4efv75Z3h5eYHjODx69AhpaWno16+fRD0uLi64ePFimcnQ09MTjRo1kprVqiwpKSkYPHgwJUNCiGIUuQ0unum69GutPj4+8PX1Lfe7hw4dwuvXr+Hu7g4ASEtLAwCYmppKlDM1NUVSUlK5dcXExMDExKTcMsWUGX5HyZCQWkwMNYV7k1NSUmBoaMjvV+S11z179mDAgAGwsLCQ2F96sXvGmNS+knx8fKCvry/3fMWWLVuGunXrKlweoGRISK1WyDgUyhlHWHzc0NBQIhnKk5SUhJMnTyIsLIzfVzw3aVpamsRrtenp6VKtxZJ8fHwUPi9QNGGLUNSbTEgtls80FNqUsXfvXpiYmGDQoEH8PhsbG5iZmfE9zEDRc8Xo6Gg4OzsLPsfatWvx+vVrpeIrjZIhIbWYqt5NFovF2Lt3L9zc3KCh8b9kynEcPD094e/vj/DwcMTGxsLd3R26uroYP3684PP4+/vj5cuXgr8nC90mE1KLiQG5t8nKrI538uRJJCcnY/LkyVLHvL29kZ2djVmzZuHVq1dwcnLCiRMnlOr0YBW4jKlCyVCZqfOFPFsghFQNxTpQhLcM+/XrV2ai4jgOvr6+cnuiK5tCybBOnTrl9vSUxnEc7t+/D1tbW6UDI4SonmJDa6rv07S7d+9K9VQrS+Hb5D/++EOhrmrGGAYOHPhBQRFCKocq3kCpTBW5pIdCydDa2hrdu3dHvXr1FKrU1tYWmpqaHxQYIUT1akrLsG7durh//z7q16+vUHkrKyucO3cO1tbWCp9DoWSYmJiocIUAEBsbK6g8IaRqKNJbXB1mun79+jWOHTsGIyMjhcq/ePEChYWFgs5BvcmE1GI1aUEoNzc3ldavVDI8deoUTp06hfT0dIjFkh3vP/74Y4UERghRPcUmaqj6lmHpPKMKgpPhypUrsWrVKnTo0AHm5uaCepkJIdVLPlOHOit/Jpj8ChzLV50JTobfffcdgoODMWnSJFXEQwipRGKmBrGcDhJ5xz8WgpNhXl6eUu8QEkKqn0JA7rT+wrohai7BKX/q1KnYt2+fKmIhhFSy4pahvK02UKhl6OXlxf8sFouxa9cunDx5Eq1atZIaTxgYGFixERJCVKamjDOsDAolwxs3bkh8btOmDQAaT0hITccUeAOFVbM3UK5fvw5NTU04OjoCAP7880/s3bsXDg4O8PX1hZaWllL1KpQMz5w5o1TlhJDqrSa2DKdPn44lS5bA0dERjx49wtixYzFixAj8/vvveP/+PYKCgpSqV3AHyuTJk7Flyxap6XaysrIwZ86cajvO0PpoFjQ0qs+j4KOjtKs6BJmeLK9+nWMW53OqOgQpIwK8qzoECYV5OQCWCf5eTRp0Xez+/fv83envv/+O7t27Y9++fbhw4QLGjh2rdDIUnPJDQkKQnZ0ttT87OxuhoaFKBUEIqRqqmtxVlRhj/CDskydP8hPDWFpaIiMjQ+l6FW4ZZmZmgjEGxhjevn0Lbe3/tWwKCwvx999/K7xyFSGkeihQYNB1AVP92x9CdOjQAX5+fujTpw+io6Oxc+dOAEVzKJS3joo8CifD4jkNOY6Dvb291HGO47By5UqlAyGEVD4hC0JVF0FBQZgwYQIOHTqE5cuXo0mTJgCKphn8kDHQCifDM2fOgDGGXr164eDBgxJzG2ppacHa2rrCJlkkhFSOmvjMsFWrVrh9+7bU/g0bNii8yLwsCidDFxcXFBQUwNXVFR06dKjQSRUJIVWDKTCominRm/z06VMsXrwYx44dQ3Z2Nuzt7bFnzx60b9/+vzoZVq5ciV27dvHroGzfvh0tWrRQ+Bx5eXkyJ4uxsrISHC8gsANFQ0MDBw8eFDxPGCGkeioEp9AmxKtXr9ClSxdoamri2LFjuHv3LjZt2oQ6derwZdavX4/AwEB8++23uHr1KszMzNC3b1+8fftWbv33799Ht27doKOjA2tra9jY2MDGxgaNGzeGjY2N0EvAEzy0pnfv3oiKioK7u7vSJyWEVA9iJv82WCxw0pp169bB0tISe/fu5fc1btyY/5kxhqCgICxfvhwjR44EUDRKxdTUFPv27cP06dPLrd/DwwMaGho4cuRIhc6cJTgZDhgwAEuXLkVsbCzat28PPT09ieNDhw5VuC5fX1+pThdTU1OkpaVJlZ0+fTp27dqFzZs3w9PTU2jYhBAZhMxaU3qVTJFIBJFIJFX+r7/+Qv/+/TFq1ChER0ejYcOGmDVrFqZNmwagqNc3LS0N/fr1k6jLxcUFFy9elJsMb968iWvXrqF58+YK/Y6KEpwMZ86cCUD2O8gcxwm+hW7RogVOnjzJf5b1APTQoUP4559/qIOGkAomZEGo0v0EPj4+Mpf7fPToEXbu3AkvLy8sW7YMV65cwdy5cyESieDq6so3dkoPgzE1NUVSUpLcmB0cHD5oPGFZBCfDip5xVkNDA2ZmZmUef/r0Kb766itERERg0KBBFXpuQmo7IUNrUlJSJNZDl9UqBIpyRIcOHeDv7w8AaNu2Le7cuYOdO3fC1dWVL1f69pYxptAt77p16+Dt7Q1/f384OjpKTRaj7JrtVb4GSkJCAiwsLCASieDk5AR/f39+vWWxWIxJkyZh0aJFCvcy5ebmIjc3l/9cumlPCPmfAqYONbG8QddFxw0NDRVKNObm5nBwcJDY98knn+DgwYMAwDd+0tLSYG5uzpdJT09XaNB0nz59ABT1X5RUnEyV7eBVKhlGR0dj48aNiIuLA8dx+OSTT7Bo0SJ069ZNUD1OTk4IDQ2Fvb09nj9/Dj8/Pzg7O+POnTuoV68e1q1bBw0NDcydO1fhOgMCAmjwNyEKUsWsNV26dEF8fLzEvvv37/PLdtrY2MDMzAyRkZFo27YtgKJhMtHR0Vi3bp3c+lU1cYzgZPjzzz/Dw8MDI0eOxNy5c8EYw8WLF9G7d28EBwdj/PjxCtc1YMAA/mdHR0d07twZdnZ2CAkJgYuLC7Zs2YLr168L6i1aunSpxPyLmZmZNCaSkDKoYtD1/Pnz4ezsDH9/f4wePRpXrlzBrl27sGvXLgBFt8eenp7w9/dH06ZN0bRpU/j7+0NXV1eh/OHi4iIoHkUJToZr1qzB+vXrMX/+fH7fvHnzEBgYiNWrVwtKhqXp6enB0dERCQkJUFNTQ3p6usQAysLCQixYsABBQUF4/PixzDrK6uEihEhTxRoon376KcLDw7F06VKsWrUKNjY2/Ct0xby9vZGdnY1Zs2bxg65PnDghNRtWWV6/fo09e/bwd6cODg6YPHmywusqyyJ4aPmjR48wZMgQqf1Dhw4VvNh8abm5uYiLi4O5uTkmTZqEW7du4ebNm/xmYWGBRYsWISIi4oPOQwgpUtwylLcJNXjwYNy+fRs5OTmIi4vjh9UU4zgOvr6+SE1NRU5ODqKjo9GyZUuF6o6JiYGdnR02b96Mly9fIiMjA4GBgbCzs8P169cFx1pMcMvQ0tISp06d4l+OLnbq1CnBt6MLFy7EkCFDYGVlhfT0dPj5+SEzMxNubm6oV68e6tWrJ1FeU1MTZmZmaNasmdCwCSEyCBlaU13Mnz8fQ4cOxe7du6GhUZTCCgoKMHXqVHh6euLs2bNK1Ss4GS5YsABz587FzZs34ezsDI7jcP78eQQHB2PLli2C6nry5AnGjRuHjIwMNGjQAJ06dcLly5f5B62EENWqiRM1xMTESCRCoGiInre3Nzp06KB0vUoNujYzM8OmTZvw22+/ASjqNj9w4ACGDRsmqK79+/cLKl/Wc0JCiHJqYjI0NDREcnKy1BsoKSkpCj9zlEWpoTUjRozAiBEjlD4pIaR6qInJcMyYMZgyZQo2btwocXe6aNEijBs3Tul6lR50XdHT5xBCKl8h48DJXRCqeiXDjRs3guM4uLq6oqCgAEBRf8LMmTOxdu1apesVnAwTEhIwefJkXLx4UWL/h47+JoRUvprYMtTS0sKWLVsQEBCAhw8fgjGGJk2aQFdX94PqFZwM3d3dVTJ9DiGk8tXEZFhMV1eXXzu5IghOhqqaPocQUvlqSjIcOXIkgoODYWhoyM+BWJawsDClziE4Gapq+hxCSOWrKcnQyMiIvws1NDRUyR2p4GSoqulzCCGVjzEOTE6yk3e8MpScNTs4OFgl5xD8Ol6fPn1w+fJl9O7dGyYmJjA2NoaxsTHq1KkDY2NjVcRICFGR4jdQ5G3VSa9evfD69Wup/ZmZmejVq5fS9QpuGapq+hxCSOWrKbfJJUVFRSEvL09qf05ODs6dO6d0vYKToaqmzyGEVL6acpsMALdu3eJ/vnv3rsRaSYWFhTh+/DgaNmyodP1VPtM1IaTq1KSWYZs2bcBxHDiOk3k7rKOjg23btildPyVDQmoxsVgNhWI58xnKOV5ZEhMTwRiDra0trly5ggYNGvDHtLS0YGJiInNBOUVRMiSkFmMAmJx1kQUum6wyxbNZVfSidMUoGRJSi4nBgath8xkWu3v3LpKTk6U6U4Ss3V6S4GSYnZ0Nxhj/HmBSUhLCw8Ph4OAgsSg0IaT6q0kdKMUePXqEESNG4Pbt2+A4Duy/pm3xQGxl50cQ/DBg2LBhCA0NBVC0DoGTkxM2bdqEYcOGYefOnUoFQQipGqqa9l+V5s2bBxsbGzx//hy6urq4c+cOzp49iw4dOiAqKkrpegUnw+vXr/NLgv7xxx8wNTVFUlISQkNDsXXrVqUDIYRUPsYU24Tw9fXle32Lt+K1kovOyeDr6wsLCwvo6OigR48euHPnjsL1X7p0CatWrUKDBg2gpqYGNTU1dO3aFQEBAYKWFS5NcDJ8//49P5vsiRMnMHLkSKipqaFTp05ISkpSOhBCSOUrvk2WtwnVokULpKam8tvt27f5Y+vXr0dgYCC+/fZbXL16FWZmZujbty/evn2rUN2FhYXQ19cHANSvXx/Pnj0DUNTBUnq9ZiEEJ8MmTZrg0KFDSElJQUREBP+cMD09nd5LJqSGUVUy1NDQgJmZGb8VD4NhjCEoKAjLly/HyJEj0bJlS4SEhOD9+/fYt2+fQnW3bNmSH4Dt5OSE9evX48KFC1i1ahVsbW0Fx1pMcDJcsWIFFi5ciMaNG6Njx47o3LkzgKJWYtu2bZUOhBBS+YQ8M8zMzJTYcnNzy6w3ISEBFhYWsLGxwdixY/Ho0SMARWMF09LSJDpbRSIRXFxcpCaMLsvXX3/ND6/x8/NDUlISunXrhr///vuDHtUJ7k3+4osv0LVrV6SmpqJ169b8/t69e9O6KITUMGIxwInlDK35b1hf6aWAfXx84OvrK1XeyckJoaGhsLe3x/Pnz+Hn5wdnZ2fcuXOHf4XO1NRU4jvFfQ+K6N+/P/+zra0t7t69i5cvX8LY2PiDpvZSapyhmZkZ3r17h8jISHTv3h06Ojr49NNPadZrQmoYIUNrUlJSJB6FiUQimeUHDBjA/+zo6IjOnTvDzs4OISEh6NSpEwBI5YriZUMUERISgi+++AJ6enr8vrp16yr03fIIvk1+8eIFevfuDXt7ewwcOBCpqakAgKlTp2LBggUfHBAhpPIwBTegaK7SkltZybA0PT09ODo6IiEhge9VLjnJAlDU51C6tViWhQsXwsTEBGPHjsWRI0f4RaE+lOBkOH/+fGhqaiI5OVliAZYxY8bg+PHjFRIUIaRyqKoDpaTc3FzExcXB3NwcNjY2MDMzQ2RkJH88Ly8P0dHRcHZ2Vqi+1NRUHDhwAOrq6hg7dizMzc0xa9YshZ85lkXwbfKJEycQERGBRo0aSexv2rQpDa0hpKYp2fQrr4wACxcuxJAhQ2BlZYX09HT4+fkhMzMTbm5u4DgOnp6e8Pf3R9OmTdG0aVP4+/tDV1cX48ePV6h+DQ0NDB48GIMHD8b79+8RHh6Offv2oWfPnmjUqBEePnwoLODieoV+ISsrS+aSfBkZGQo3mwkh1YQiLT+BLcMnT55g3LhxyMjIQIMGDdCpUydcvnyZn2jB29sb2dnZmDVrFl69egUnJyecOHGCH78shK6uLvr3749Xr14hKSkJcXFxgusoJjgZdu/eHaGhoVi9ejWAogehYrEYGzZsQM+ePZUOhBBS+RR5w0ToGyj79+8v9zjHcfD19ZXZE62o4hbhL7/8gpMnT8LS0hLjxo3D77//rnSdgpPhhg0b0KNHD8TExCAvLw/e3t64c+cOXr58iQsXLigdiMpdjQU4TfnlKsminzyqOgSZRE4vqzoEKWpr7lV1CFLqR1d1BJIKWL5S36uJEzWMGzcOhw8fhq6uLkaNGoWoqCiFnzeWR6mlQm/duoUdO3ZAXV0dWVlZGDlyJGbPng1zc/MPDogQUokYJ/82uJolQ47jcODAAfTv3x8aGhU3C6HS4wxXrVpVYUEQQqoGExdt8spUF/n5+UhNTUXTpk0rNBECSgytAYBz585h4sSJcHZ2xtOnTwEAP/30E86fP1+hwRFCVKsyhtZUJE1NTcTGxqrkBQ/ByfDgwYPo378/dHR0cP36df79xLdv38Lf37/CAySEqJgiI66rEVdXV+zZs6fC6xXczvTz88N3330HV1dXiV4jZ2dnunUmpIapiR0oeXl5+OGHHxAZGYkOHTpIvJYHAIGBgUrVK7hlGB8fj+7du0vtNzQ0lLnKfXnkTQLp6+uL5s2bQ09PD8bGxujTpw/++ecfoSETQsoi5H28aiI2Nhbt2rWDoaEh7t+/jxs3bvDbzZs3la5XcMvQ3NwcDx48QOPGjSX2nz9/Xqm5xFq0aIGTJ0/yn0su9Wdvb49vv/0Wtra2yM7OxubNm9GvXz88ePBAYplAQoiyuP82eWWqjzNnzqikXsHJcPr06Zg3bx5+/PFHcByHZ8+e4dKlS1i4cCFWrFghPID/JoGUpfTrOYGBgdizZw9u3bqF3r17Cz4XIaQUFbyOV1kePHiAhw8f8jNnCZn5RhbBydDb2xtv3rxBz549kZOTg+7du0MkEmHhwoX46quvBAdQPAmkSCSCk5MT/P39ZbYw8/LysGvXLhgZGUnMo1habm6uxKSTmZmZgmMipNaogcnwxYsXGD16NM6cOQOO45CQkABbW1tMnToVderUwaZNm5SqV9Azw8LCQkRHR2PBggXIyMjAlStXcPnyZfzf//0f/3qeEMWTQEZERGD37t1IS0uDs7MzXrx4wZc5cuQI9PX1oa2tjc2bNyMyMhL169cvs86AgAAYGRnxW+kJKQkhJRQPupa3VSOqmjlLUDJUV1dH//798ebNG+jq6qJDhw7o2LEjvziLUAMGDMDnn38OR0dH9OnTB0ePHgVQNHljsZ49e+LmzZu4ePEiPvvsM4wePRrp6ell1rl06VK8efOG31JSUpSKjZDaoHjQtbytOjlx4gTWrVtX4TNnCe5NdnR05NczqGglJ4Esua9Jkybo1KkT9uzZAw0NjXLHGIlEIqlJKAkhZaiBLUNVzZwlOBmuWbMGCxcuxJEjR5Camiq1SMyHKDkJZFkYY+UuREMIURzHFNuqk+KZs4pV1MxZgjtQPvvsMwDA0KFDJXpuintyCgsLFa6rvEkgs7KysGbNGgwdOhTm5uZ48eIFduzYgSdPnmDUqFFCwyaEyFIDO1BUNXOW4GRYkWN8ypsEMicnB/fu3UNISAgyMjJQr149fPrppzh37hxatGhRYTEQUqvVwFlrimfO2rlzZ4XOnCU4Gbq4uCh9stLKmwRSW1sbYWFhFXYuQogMNbBlCBTNnLVy5coKrVNwMixeyb40juOgra0NKysrmv6fkJqiBibD48ePQ19fH127dgUAbN++Hbt374aDgwO2b98OY2NjpeoV3IHSpk0btG3bVmpr06YNmjdvDiMjI7i5uSEnJ0epgAghlagS3k0OCAjgF4LiT8sYfH19YWFhAR0dHfTo0QN37txRqL5FixbxnbW3b9+Gl5cXBg4ciEePHsHLy0vpOAUnw/DwcDRt2hS7du3CzZs3cePGDezatQvNmjXDvn37sGfPHpw+fRpff/210kERQiqJiofWXL16Fbt27UKrVq0k9q9fvx6BgYH49ttvcfXqVZiZmaFv3754+/at3DoTExPh4OAAoGhKwSFDhsDf3x87duzAsWPHlI5V8G3ymjVrsGXLFvTv35/f16pVKzRq1AjffPMNrly5Aj09PSxYsAAbN25UOjBCiOopMnRG2aE17969w4QJE7B79274+fnx+xljCAoKwvLlyzFy5EgARS9amJqaYt++fZg+fXq59WppaeH9+/cAgJMnT8LV1RUAULdu3Q8a3ie4ZXj79m1+yb+SrK2tcfv2bQBFt9KpqalKB0UIqSQCbpNLjymWN9539uzZGDRoEPr06SOxPzExEWlpaejXrx+/TyQSwcXFRaGF4Lt27QovLy+sXr0aV65cwaBBgwAA9+/fl3orRQjBybB58+ZYu3Yt8vLy+H35+flYu3YtmjdvDgB4+vQpTE1NlQ6KEFI5OCgw6Pq/spaWlhLv/QcEBJRZ7/79+3H9+nWZZdLS0gBAKkeYmpryx8rz7bffQkNDA3/88Qd27tyJhg0bAgCOHTvGj4NWhuDb5O3bt2Po0KFo1KgRWrVqBY7jcOvWLRQWFuLIkSMAgEePHmHWrFlKB0UIqSQCxhmmpKRIvN5a1qiRlJQUzJs3DydOnIC2tnaZ1ZaebkvRKbisrKz4XFPS5s2b5X63PIKTobOzMx4/foyff/4Z9+/fB2MMX3zxBcaPHw8DAwMAwKRJkz4oKEJIJREwtEbRd/2vXbuG9PR0tG/fnt9XWFiIs2fP4ttvv0V8fDyAohZiyUHS6enpCt9RFhYWIjw8HHFxceA4Ds2bN8fw4cM/aMU8pb6pr6+PGTNmKH1SQkg1oYJxhr179+b7D4p5eHigefPmWLx4MWxtbWFmZobIyEi0bdsWQNF8pdHR0Vi3bp3c+mNjYzF06FA8f/4czZo1A1D0vLBBgwb466+/4OjoKCzg/yi1VOhPP/2Erl27wsLCgp8yZ/Pmzfjzzz+VCoIQUjVUMVGDgYEBWrZsKbHp6emhXr16aNmyJT/m0N/fH+Hh4YiNjYW7uzt0dXWlZreXZerUqWjZsiWePHmC69ev4/r160hJSUGrVq3w5ZdfKnkllEiGO3fuhJeXFwYMGIBXr17xEzMYGxsjKChI6UAIIVWgihaE8vb2hqenJ2bNmoUOHTrg6dOnOHHiBP+orTz//vsvAgICJN40MTY2xpo1az5oQSjByXDbtm3YvXs3li9fLnF/3qFDB6mmMSGkmqukZBgVFSXRWOI4Dr6+vkhNTUVOTg6io6PRsmVLhepq1qwZnj9/LrU/PT0dTZo0UTpGwc8MExMT+fv8kkQiEbKyspQOhBBS+VQ56LoilRxM7e/vj7lz58LX1xedOnUCAFy+fBmrVq1S6JljWQQnQxsbG9y8eVNq4PWxY8f4V2QIITWEmCva5JWpYnXq1JGaP3X06NH8PsaKMvaQIUMEzalakuBkuGjRIsyePRs5OTlgjOHKlSv49ddfERAQgB9++EGpIAghVaOmtAxVtVZySYKToYeHBwoKCuDt7Y33799j/PjxaNiwIbZs2YKxY8eqIkZCiKrUkCm8KnIe1bIoNc5w2rRpmDZtGjIyMiAWi2FiYlLRcRFCKoMiQ2eqQTIs7fXr19izZw8/6NrBwQGTJ0+GkZGR0nUqNc6wWP369SkRElKTVdHQmg8RExMDOzs7bN68GS9fvkRGRgYCAwNhZ2eH69evK12vQi3Dtm3bKvTOIIAPCoYQUslqyG1ySfPnz8fQoUOxe/dufnhfQUEBpk6dCk9PT5w9e1apehVKhsOHD+d/zsnJwY4dO+Dg4IDOnTsDKOrWvnPnDk3OQEgNU1M6UEqKiYmRSIQAoKGhAW9vb3To0EHpehVKhj4+PvzPU6dOxdy5c7F69WqpMikpKUoHQgghijA0NERycjI/ZWCxlJQUhd5gKYvgZ4a///47P7NsSRMnTsTBgweVDoQQUgVq4DPDMWPGYMqUKThw4ABSUlLw5MkT7N+/H1OnTsW4ceOUrldwb7KOjg7Onz+Ppk2bSuw/f/58uXOXEUKqH44BnFh+mepk48aN4DgOrq6uKCgoAABoampi5syZWLt2rdL1Ck6Gnp6emDlzJq5duybxKsyPP/6IFStWKB0IIaQK1MAOFC0tLWzZsgUBAQF4+PAhGGNo0qQJdHV1P6hewclwyZIlsLW1xZYtW7Bv3z4AwCeffILg4GCMHj36g4IhhFSumtiBUkxXV1fpuQtlUWrQ9ejRoynxEfIxqIEtQ1VRfo5sQkiNV5NbhhVNod7kunXrIiMjQ+FKrays+BmwCSHVWA3sTVYVhVqGr1+/xrFjxxR+7+/FixdKT6NDCKlENeQ2uV27djh16hSMjY2xatUqLFy48IM7TEpT+DbZzc2tQk9c2zHHt1UdgkwW3tXvf2IFVR3AR6ym3CbHxcUhKysLxsbGWLlyJWbMmFE1yVAsljMQiRBSM6mgZbhz507s3LkTjx8/BgC0aNECK1aswIABA4qqYwwrV67Erl278OrVKzg5OWH79u1o0aJFmXW2adMGHh4e6Nq1Kxhj2LhxI/T19WWWVXaIH3WgEFKLcWIFBl0LbAs1atQIa9eu5dcjCQkJwbBhw3Djxg20aNEC69evR2BgIIKDg2Fvbw8/Pz/07dsX8fHxZb5OFxwcDB8fHxw5cgQcx+HYsWMy10jmOI6SISFECSpoGQ4ZMkTi85o1a7Bz505cvnwZDg4OCAoKwvLlyzFy5EgARcnS1NQU+/btw/Tp02XW2axZM+zfvx8AoKamhlOnTlX49IEfNJ/hh/L19QXHcRKbmZkZACA/Px+LFy+Go6Mj9PT0YGFhAVdXVzx79qwqQybkoyJk3eTMzEyJLTc3V279hYWF2L9/P7KystC5c2ckJiYiLS0N/fr148uIRCK4uLjg4sWLCsWsqgmlqzQZAkXPE1JTU/mteLnR9+/f4/r16/jmm29w/fp1hIWF4f79+xg6dGgVR0zIR0TA0BpLS0sYGRnxW0BAQJnV3r59G/r6+hCJRJgxYwbCw8Ph4OCAtLQ0AICpqalEeVNTU/6YIh4+fIg5c+agT58+6Nu3L+bOnYuHDx8q/H1Zqvw2WUNDg28NlmRkZITIyEiJfdu2bUPHjh2RnJwMKyurygqRkI+XgNvklJQUGBoa8rtFIlGZX2nWrBlu3ryJ169f4+DBg3Bzc0N0dDR/vPRk0YwxhSeQjoiIwNChQ9GmTRt06dIFjDFcvHgRLVq0wOHDh9G3b1+F6ilNqWT48OFD7N27Fw8fPsSWLVtgYmKC48ePw9LSstweIVkSEhJgYWEBkUgEJycn+Pv7w9bWVmbZN2/egOM41KlTp8z6cnNzJZrvJddbJYRI4v7b5JUBiuYRLJkMy6OlpcV3oHTo0AFXr17Fli1bsHjxYgBAWloazM3N+fLp6elSrcWyLFmyBPPnz5eaoWbJkiVYvHix0slQ8G1ydHQ0HB0d8c8//yAsLAzv3r0DANy6dUtiElhFODk5ITQ0FBEREdi9ezfS0tLg7OyMFy9eSJXNycnBkiVLMH78+HL/IAEBARJNeUtLS2G/ICG1SSW9gcIYQ25uLmxsbGBmZiZx15eXl4fo6Gg4OzsrVFdcXBymTJkitX/y5Mm4e/eu0jEKToZLliyBn58fIiMjoaWlxe/v2bMnLl26JKiuAQMG4PPPP4ejoyP69OmDo0ePAijqXSopPz8fY8eOhVgsxo4dO8qtc+nSpXjz5g2/0ezbhJRNSAeKopYtW4Zz587h8ePHuH37NpYvX46oqChMmDABHMfB09MT/v7+CA8PR2xsLNzd3aGrq4vx48crVH+DBg1w8+ZNqf03b978oI4VwbfJt2/f5qfuKqlBgwYyW3RC6OnpwdHREQkJCfy+/Px8jB49GomJiTh9+rTcZrpIJCr3WQYhpAQVDK15/vw5Jk2ahNTUVBgZGaFVq1Y4fvw4f/vq7e2N7OxszJo1ix90feLECYWn7J82bRq+/PJLPHr0CM7OzuA4DufPn8e6deuwYMECYcGWIDgZ1qlTB6mpqbCxsZHYf+PGDTRs2FDpQICi531xcXHo1q0bgP8lwoSEBJw5cwb16tX7oPoJITJU8Ot2e/bsKfc4x3Hw9fWFr6+vUvV/8803MDAwwKZNm7B06VIAgIWFBXx9fTF37lyl6gSUSIbjx4/H4sWL8fvvv4PjOIjFYly4cAELFy6UuTZKeRYuXIghQ4bAysoK6enp8PPzQ2ZmJtzc3FBQUIAvvvgC169fx5EjR1BYWMh3vdetW1fiFp0QohxVvIGiahzHYf78+Zg/fz7evi16x/9DFoIqJjgZrlmzBu7u7mjYsCEYY3BwcEBhYSHGjx+Pr7/+WlBdT548wbhx45CRkYEGDRqgU6dOuHz5MqytrfH48WP89ddfAIreSyzpzJkz6NGjh9DQCSGl1JSJGspSEUmwmOBkqKmpiV9++QWrVq3CjRs3IBaL0bZtW6kFohRR/HqNLI0bNwZj1fivQMjHoIZM4VUZlB50bWdnBzs7u4qMhRBSyWp6y7AiKZQMvby8FK4wMDBQ6WAIIZWMWoY8hZLhjRs3JD5fu3YNhYWFaNasGQDg/v37UFdXR/v27Ss+QkKI6tSwZJifn49+/frh+++/h729fYXWrVAyPHPmDP9zYGAgDAwMEBISAmNjYwDAq1ev4OHhwQ+JIYTUDDXtNllTUxOxsbEKv8cshOA3UDZt2oSAgAA+EQKAsbEx/Pz8sGnTpgoNjhCiYjVwQShXV1e5YxmVIbgDJTMzE8+fP5eakCE9PZ0f80MIqRk4xsDJGbUh73hly8vLww8//IDIyEh06NABenp6EseV7bcQnAxHjBgBDw8PbNq0CZ06dQIAXL58GYsWLeJnriWE1Aw1cdB1bGws2rVrB6Cov6KkD7l9FpwMv/vuOyxcuBATJ05Efn5+USUaGpgyZQo2bNigdCCEkCpQwzpQAMk+jIokOBnq6upix44d2LBhAx4+fAjGGJo0aSLVVCWEVH81rQOlpAcPHuDhw4fo3r07dHR0BE0QK4vSg6719PTQqlUrpU9MCKkGamDL8MWLFxg9ejTOnDkDjuOQkJAAW1tbTJ06FXXq1FG6I1dwMuzZs2e52ff06dNKBUIIqXw1sWU4f/58aGpqIjk5GZ988gm/f8yYMZg/f37lJcPSkybk5+fj5s2biI2NhZubm1JBEEKqSA1sGZ44cQIRERFo1KiRxP6mTZsiKSlJ6XoFJ8PNmzfL3O/r68svAUAIqTmqW8tPnqysLOjq6krtz8jI+KCJnStsqdCJEyfixx9/rKjqCCGVgTHFtmqke/fuCA0N5T8Xz6u6YcMG9OzZU+l6K2yp0EuXLkFbW7uiqiOEVIKa+Mxww4YN6NGjB2JiYpCXlwdvb2/cuXMHL1++xIULF5SuV3AyLD2wmjGG1NRUxMTE4JtvvlE6EEJI5eMKAU7O/SFXKKzOgIAAhIWF4d69e9DR0YGzszPWrVvHT+wCFOWNlStXYteuXfw6KNu3b1doqWEHBwfcunULO3fuhLq6OrKysjBy5EjMnj1bYvlRoQQnQ0NDQ4neZDU1NTRr1gyrVq1Cv379lA6EEFIFVNCBEh0djdmzZ+PTTz9FQUEBli9fjn79+uHu3bv8eOT169cjMDAQwcHBsLe3h5+fH/r27Yv4+HiFZq82MzPDypUrhQUmh+BkGBwcXKEBEEKqjipuk48fPy7xee/evTAxMcG1a9fQvXt3MMYQFBSE5cuX83eaISEhMDU1xb59+zB9+nS553j16hX27NmDuLg4cByHTz75BB4eHqhbt66wYEsQ3IFia2src0nQ169fw9bWVulACCFVQEAHSmZmpsSWm5ur0CnevHkDAHyiSkxMRFpamsSdpEgkgouLCy5evCi3vujoaNjY2GDr1q149eoVXr58ia1bt8LGxgbR0dFCrwBPcDJ8/PgxCgulHyLk5ubi6dOnSgdCCKl8QhaRt7S0hJGREb8FBATIrZ8xBi8vL3Tt2hUtW7YEAH6VS1NTU4mypqam/LHyzJ49m19LPSwsDGFhYXj06BHGjh2L2bNnC7wC/6PwbXLxSnUAEBERASMjI/5zYWEhTp06hcaNGysdCCGkCgh4ZpiSkgJDQ0N+tyJj+r766ivcunUL58+flzpW+k02Rd8tfvjwIQ4ePAh1dXV+n7q6Ory8vCSG3AilcDIcPnw4gKJfoPSbJpqammjcuDFN7kpIDSPkmaGhoaFEMpRnzpw5+Ouvv3D27FmJt0XMzMwAFLUQS/b+pqenS7UWZWnXrh3i4uIkeqcBIC4uTuoNOSEUToZicdGkZjY2Nrh69Srq16+v9EkJIdWEIoOqBQ66Zoxhzpw5CA8PR1RUFGxsbCSO29jYwMzMDJGRkWjbti2Aoglbo6OjsW7dOpl13rp1i/957ty5mDdvHh48eCAxp+r27duxdu1aQbGWJLg3OTExUemTEUKqF1X0Js+ePRv79u3Dn3/+CQMDA/45oJGREXR0dMBxHDw9PeHv74+mTZuiadOm8Pf3h66uLsaPHy+zzjZt2oDjOIm11L29vaXKjR8/HmPGjBEW8H8USoZbt27Fl19+CW1tbWzdurXcsnPnzlUqEEJI5VPFTNc7d+4EAPTo0UNi/969e+Hu7g6gKJFlZ2dj1qxZ/KDrEydOlDnGsDIaYRxj8tvANjY2iImJQb169aSavBKVcRwePXpUoQF+qMzMTBgZGaEHhkGD06zqcHjJfzhWdQgy2X6dU9UhSCm4l1DVIVR7BSwfUfgTb968Uei5XvG/C+e+K6GhWf5rtAX5ObgY6aNw3TWVQi3DklmZbpMJ+YjUwCm8AODp06e4cOEC0tPT+f6MYsrenQp+Zrhq1SosXLhQagqd7OxsbNiwAStWrFAqEEJI5eOgwDPDSolEcXv37sWMGTOgpaWFevXqSQzH4ThO6WQoeND1ypUrZc5b+P79+wp/V5AQomI1cAqvFStWYMWKFXjz5g0eP36MxMREfvuQx3SCW4ZlDYz8999/P+i9QEJI5auJU3i9f/8eY8eOhZpahU3HCkBAMjQ2NgbHceA4Dvb29hIJsbCwEO/evcOMGTMqNDhCiIrVwGeGU6ZMwe+//44lS5ZUaL0KJ8OgoCAwxjB58mSsXLlS4nU8LS0tNG7cGJ07d67Q4CqSuoE+1Dmtqg6DZ767ek6EW3DvdlWHQCoRxxg4ObfB8o5XtoCAAAwePBjHjx+Ho6MjNDUlR4kEBgYqVa/CybD4FTwbGxs4OztLBUAIqYHE/23yylQj/v7+iIiI4F/HK92BoizBzwxdXFz4n7Ozs5Gfny9xXMg4JF9fX6lOl5IzV4SFheH777/HtWvX8OLFC9y4ceOD3j0khEiqiS3DwMBA/Pjjj/wA7ooi+Ank+/fv8dVXX8HExAT6+vowNjaW2IRq0aIFUlNT+e327f/dpmVlZaFLly4f9L4hIaQcYqbYVo2IRCJ06dKlwusV3DJctGgRzpw5gx07dsDV1RXbt2/H06dP8f333yuVtDQ0NPhZLEqbNGkSgKI5FAkhFa8m9ibPmzcP27Ztk/tqsFCCk+Hhw4cRGhqKHj16YPLkyejWrRuaNGkCa2tr/PLLL5gwYYKg+hISEmBhYQGRSAQnJyf4+/t/0IzZubm5EjPwZmZmKl0XIR89Fcxao2pXrlzB6dOnceTIEbRo0UKq/yIsLEypegUnw5cvX/LvJxsaGuLly5cAgK5du2LmzJmC6nJyckJoaCjs7e3x/Plz+Pn5wdnZGXfu3EG9evWEhgagqKeJBn8TohhVTNSganXq1JFapbMiCE6Gtra2ePz4MaytreHg4IDffvsNHTt2xOHDh1GnTh1BdQ0YMID/2dHREZ07d4adnR1CQkLg5eUlNDQAwNKlSyW+m5mZCUtLS6XqIuSjVwNbhnv37lVJvYKToYeHB/7991+4uLhg6dKlGDRoELZt24aCggKlx/cU09PTg6OjIxISlJ+lRCQSKTQdOSEENXLQtaoITobz58/nf+7Zsyfu3buHmJgY2NnZoXXr1h8UTG5uLuLi4tCtW7cPqocQopiaOLTGxsam3PGEyr6fLDgZlmZlZQUrKyukpKRg8uTJ+PHHHxX+7sKFCzFkyBBYWVkhPT0dfn5+yMzM5Ad4v3z5EsnJyXj27BkAID4+HkDRGgpl9UATQgSogbfJnp6eEp/z8/Nx48YNHD9+HIsWLVK63g9OhsVevnyJkJAQQcnwyZMnGDduHDIyMtCgQQN06tQJly9fhrW1NYCiFfk8PDz48mPHjgUA+Pj4wNfXt6JCJ6T2YpD/hkn1yoWYN2+ezP3bt29HTEyM0vVWWDJUxv79+8s97u7uXuGjzAkh/8OJGTg53cVcNRt0XZYBAwZg6dKlSnewVGkyJIRUsRp4m1yWP/7444OmEazYCcEIITWLWMFNgLNnz2LIkCGwsLAAx3E4dOiQxHHGGHx9fWFhYQEdHR306NEDd+7cUbj+tm3bol27dvzWtm1bmJubY9myZVi2bJmwYEtQuGUob5Dj69evlQ6CEFI1VNGbnJWVhdatW8PDwwOff/651PH169cjMDAQwcHBsLe3h5+fH/r27Yv4+PgyV8crafjw4RKf1dTU0KBBA/To0QPNmzcXFGtJCifDkvMXlnXc1dVV6UAIIVVAwG1y6VdbyxrTO2DAAIkXKiSrYggKCsLy5cv5BlZISAhMTU2xb98+TJ8+XW7IPj4+cssoQ+FkqKpR34SQKiQgGZZ+k0uZUR2JiYlIS0tDv379+H0ikQguLi64ePGiQslQVagDhZDaTEAyTElJkZivVJk3vYrnKjU1NZXYb2pqiqSkpHK/q6amJnfyVo7jUFBQIDgugJIhIbWbGPLXAv2vA8XQ0LDCFpEvndTKWmiupPDw8DKPXbx4Edu2bQP7gJ5vSoaE1GKV/Tpe8ZtjaWlpMDc35/enp6dLtRZLGzZsmNS+e/fuYenSpTh8+DAmTJiA1atXKx0bDa0hpDYrFCu2VRAbGxuYmZkhMjKS35eXl4fo6Gg4OzsrXM+zZ88wbdo0tGrVCgUFBbh58yZCQkJgZWWldGzUMiSkNlPBoOt3797hwYMH/OfExETcvHkTdevWhZWVFTw9PeHv74+mTZuiadOm8Pf3h66uLsaPHy+37jdv3sDf3x/btm1DmzZtcOrUqQqb2IWSISG1mgLJUODLyTExMejZsyf/uXh+UTc3NwQHB8Pb2xvZ2dmYNWsWXr16BScnJ5w4cULuGMP169dj3bp1MDMzw6+//irztvlDcOxDnjjWAJmZmTAyMkJvgwnQqEbrJuc4Kz84VJU0I65WdQhECQUsH1H4E2/evFGok6P430UfmznQUCu/V7hAnIuTidsUrltV1NTUoKOjgz59+kBdXb3McpU27T8h5CMiVmB212oyUYOrq+sHrYssDyVDQmozJi7a5JWpBoKDg1VaPyVDQmqzj2jWmg9FyZCQ2qwG3SarGiVDQmozahnyKBkSUpuJFZj3n1qGhJCPnliB2VvF1aMDRdUoGRJSm9FtMo+SISG1GSVDHiVDQmoz6k3mUTIkpBZjTAwmZ1C1vOMfC0qGhNRmjMlv+dFtMiHko8cUuE2mZEgI+eiJxQBXM95NVjVKhoTUZtQy5FEyJKQWY4WFYFxh+WVY+cc/FpQMCanNxAzgqGUIUDIkpHZjCrybXEuSIa2OR0gtxsRMoU0ZO3bsgI2NDbS1tdG+fXucO3eugqOvWJQMCanNime6lrcJdODAAXh6emL58uW4ceMGunXrhgEDBiA5OVkFv0TFoGRISC2mqpZhYGAgpkyZgqlTp+KTTz5BUFAQLC0tsXPnThX8FhXjo39mWLz4XwHLr+JIJBUU5FR1CDJx1ew6EcUUoOjvJnSxywKWK7flV1x3ZmamxH6RSASRSHplvby8PFy7dg1LliyR2N+vXz9cvHhRUHyV6aNPhm/fvgUARL/7rYojKeVkVQdAPkZv376FkZGR3HJaWlowMzPD+bS/FapXX18flpaWEvt8fHzg6+srVTYjIwOFhYUwNTWV2G9qaoq0tDSFzlcVPvpkaGFhgZSUFBgYGKh0mUGg6P+clpaWSElJqdL1ZZVBsVeNioqdMYa3b9/CwsJCofLa2tpITExEXl6ewvWX/vcjq1VYUunysuqoTj76ZKimpoZGjRpV6jkNDQ1r3D/KYhR71aiI2BVpEZakra0NbW3tDzqnLPXr14e6urpUKzA9PV2qtVidUAcKIaRCaWlpoX379oiMjJTYHxkZCWdn5yqKSr6PvmVICKl8Xl5emDRpEjp06IDOnTtj165dSE5OxowZM6o6tDJRMqxAIpEIPj4+cp+lVEcUe9WoybGXZ8yYMXjx4gVWrVqF1NRUtGzZEn///Tesra2rOrQycUxoXzwhhHyE6JkhIYSAkiEhhACgZEgIIQAoGRJCCABKhoIJnZYoOjoa7du3h7a2NmxtbfHdd99VUqT/ExAQgE8//RQGBgYwMTHB8OHDER8fX+53oqKiwHGc1Hbv3r1KirqIr6+vVAxmZmblfqc6XHMAaNy4scxrOHv2bJnlq8s1r60oGQogdFqixMREDBw4EN26dcONGzewbNkyzJ07FwcPHqzUuKOjozF79mxcvnwZkZGRKCgoQL9+/ZCVlSX3u/Hx8UhNTeW3pk2bVkLEklq0aCERw+3bt8ssW12uOQBcvXpVIu7iQcijRo0q93vV4ZrXSoworGPHjmzGjBkS+5o3b86WLFkis7y3tzdr3ry5xL7p06ezTp06qSxGRaSnpzMALDo6uswyZ86cYQDYq1evKi8wGXx8fFjr1q0VLl9drzljjM2bN4/Z2dkxsVgs83h1uea1FbUMFVQ8LVG/fv0k9pc3LdGlS5ekyvfv3x8xMTHIz6+6qbLevHkDAKhbt67csm3btoW5uTl69+6NM2fOqDo0mRISEmBhYQEbGxuMHTsWjx49KrNsdb3meXl5+PnnnzF58mS5kxVUh2teG1EyVJAy0xKlpaXJLF9QUICMjAyVxVoexhi8vLzQtWtXtGzZssxy5ubm2LVrFw4ePIiwsDA0a9YMvXv3xtmzZysxWsDJyQmhoaGIiIjA7t27kZaWBmdnZ7x48UJm+ep4zQHg0KFDeP36Ndzd3cssU12ueW1Fr+MJJHRaIlnlZe2vLF999RVu3bqF8+fPl1uuWbNmaNasGf+5c+fOSElJwcaNG9G9e3dVh8kbMGAA/7OjoyM6d+4MOzs7hISEwMvLS+Z3qts1B4A9e/ZgwIAB5U6xVV2ueW1FLUMFKTMtkZmZmczyGhoaqFevnspiLcucOXPw119/4cyZM0pNa9apUyckJCSoIDLF6enpwdHRscw4qts1B4CkpCScPHkSU6dOFfzd6nDNawtKhgpSZlqizp07S5U/ceIEOnToAE1NTZXFWhpjDF999RXCwsJw+vRp2NjYKFXPjRs3YG5uXsHRCZObm4u4uLgy46gu17ykvXv3wsTEBIMGDRL83epwzWuNKu2+qWH279/PNDU12Z49e9jdu3eZp6cn09PTY48fP2aMMbZkyRI2adIkvvyjR4+Yrq4umz9/Prt79y7bs2cP09TUZH/88Uelxj1z5kxmZGTEoqKiWGpqKr+9f/+eL1M69s2bN7Pw8HB2//59Fhsby5YsWcIAsIMHD1Zq7AsWLGBRUVHs0aNH7PLly2zw4MHMwMCg2l/zYoWFhczKyootXrxY6lh1vea1FSVDgbZv386sra2ZlpYWa9euncTwFDc3N+bi4iJRPioqirVt25ZpaWmxxo0bs507d1ZyxIwBkLnt3buXL1M69nXr1jE7Ozumra3NjI2NWdeuXdnRo0crPfYxY8Ywc3NzpqmpySwsLNjIkSPZnTt3yoybsepxzYtFREQwACw+Pl7qWHW95rUVTeFFCCGgZ4aEEAKAkiEhhACgZEgIIQAoGRJCCABKhoQQAoCSISGEAKBkSAghACgZEkIIAEqG1cLjx4/BcRxu3rypkvo5jsOhQ4eU/n7J6eiHDx9ebtkePXrA09NT6XMRae7u7vz1/5C/IylfrU+G7u7ucv+Bq5qlpSVSU1P5+QWLk8/r16+rNK7S4uPjERwcXNVhfLTK+rtv2bIFqampVRNULULzGVYD6urqchc5qg5MTExQp06dqg4D+fn5VTYDTVUwMjKCkZFRVYfx0av1LUN5oqOj0bFjR4hEIpibm2PJkiUoKCjgj/fo0QNz586Ft7c36tatCzMzM/j6+krUce/ePXTt2hXa2tpwcHDAyZMnJW55St4mP378GD179gQAGBsbg+M4fnbkxo0bIygoSKLuNm3aSJwvISEB3bt3589VejorAHj69CnGjBkDY2Nj1KtXD8OGDcPjx48FX5usrCy4urpCX18f5ubm2LRpk1SZvLw8eHt7o2HDhtDT04OTkxOioqIkyuzevRuWlpbQ1dXFiBEjEBgYKJF0fX190aZNG/z444+wtbWFSCQCYwxv3rzBl19+CRMTExgaGqJXr174999/Jeo+fPiwxEp5K1eulPj7+fr6wsrKCiKRCBYWFpg7d265v7O8+gIDA+Ho6Ag9PT1YWlpi1qxZePfuHX88KSkJQ4YMgbGxMfT09NCiRQv8/fff5f7dSeWglmE5nj59ioEDB8Ld3R2hoaG4d+8epk2bBm1tbYkEVDzr8j///INLly7B3d0dXbp0Qd++fSEWizF8+HBYWVnhn3/+wdu3b7FgwYIyz2lpaYmDBw/i888/R3x8PAwNDaGjo6NQvGKxGCNHjkT9+vVx+fJlZGZmSj2/e//+PXr27Ilu3brh7Nmz0NDQgJ+fHz777DPcunULWlpaCl+fRYsW4cyZMwgPD4eZmRmWLVuGa9euoU2bNnwZDw8PPH78GPv374eFhQXCw8Px2Wef4fbt22jatCkuXLiAGTNmYN26dRg6dChOnjyJb775RupcDx48wG+//YaDBw9CXV0dADBo0CDUrVsXf//9N4yMjPD999+jd+/euH//PurWrYuIiAhMnDgRW7duRbdu3fDw4UN8+eWXAAAfHx/88ccf2Lx5M/bv348WLVogLS1NKpmWJK8+AFBTU8PWrVvRuHFjJCYmYtasWfD29saOHTsAALNnz0ZeXh7Onj0LPT093L17F/r6+h/0dycVpIpnzalybm5ubNiwYTKPLVu2jDVr1kxiNbPt27czfX19VlhYyBhjzMXFhXXt2lXie59++ik/f92xY8eYhoYGS01N5Y9HRkYyACw8PJwxxlhiYiIDwG7cuMEYK3uVNGtra7Z582aJfa1bt2Y+Pj6MsaLpotTV1VlKSgp//NixYxLn2rNnj9TvlJuby3R0dFhERITM6yArnrdv3zItLS22f/9+ft+LFy+Yjo4OmzdvHmOMsQcPHjCO49jTp08l6uvduzdbunQpY6xoiq5BgwZJHJ8wYQIzMjLiP/v4+DBNTU2Wnp7O7zt16hQzNDRkOTk5Et+1s7Nj33//PWOMsW7dujF/f3+J4z/99BMzNzdnjDG2adMmZm9vz/Ly8mT+3qXJq0+W3377jdWrV4//7OjoyHx9fWWWlbc6Xsm/I6l41DIsR1xcHDp37iyxdkaXLl3w7t07PHnyBFZWVgCAVq1aSXzP3Nwc6enpAIo6HSwtLSWeCXbs2FFl8VpZWUlM6d+5c2eJMteuXcODBw9gYGAgsT8nJwcPHz5U+FwPHz5EXl6eRP1169aVWMPj+vXrYIzB3t5e4ru5ubn8FPzx8fEYMWKExPGOHTviyJEjEvusra3RoEEDid/j3bt3UlP5Z2dn87/HtWvXcPXqVaxZs4Y/XlhYiJycHLx//x6jRo1CUFAQbG1t8dlnn2HgwIEYMmQINDRk/7OQV5+uri7OnDkDf39/3L17F5mZmSgoKEBOTg6ysrKgp6eHuXPnYubMmThx4gT69OmDzz//XOq/H1I1KBmWg8lY7InJWFyo9MN8juMgFovLrENZampq/PmLlVz+svSx0nECRbfS7du3xy+//CJVtmSykUfWuUoTi8VQV1fHtWvX+FvbYvr6+nw9ZV3jkvT09KTqNjc3l3r+CIB/3igWi7Fy5UqMHDlSqoy2tjYsLS0RHx+PyMhInDx5ErNmzcKGDRsQHR0ts4NGXn1JSUkYOHAgZsyYgdWrV6Nu3bo4f/48pkyZwv+dpk6div79++Po0aM4ceIEAgICsGnTJsyZM0eqTlK5KBmWw8HBAQcPHpT4B3vx4kUYGBigYcOGCtXRvHlzJCcn4/nz5/zCUVevXi33O8XP7QoLCyX2N2jQQGKIRWZmJhITEyXiTU5OxrNnz/hV2C5duiRRR7t27XDgwAG+00FZTZo0gaamJi5fvsy3kF+9eoX79+/DxcUFQNH6v4WFhUhPT0e3bt1k1tO8eXNcuXJFYl9MTIzc87dr1w5paWnQ0NBA48aNyywTHx+PJk2alFmPjo4Ohg4diqFDh2L27Nlo3rw5bt++jXbt2gmuLyYmBgUFBdi0aRPU1Ir6Jn/77TepcpaWlpgxYwZmzJiBpUuXYvfu3ZgzZ06Zf3dSOag3GUWLqt+8eVNiS05OxqxZs5CSkoI5c+bg3r17+PPPP+Hj4wMvLy/+P3Z5+vbtCzs7O7i5ueHWrVu4cOECli9fDqDspSutra3BcRyOHDmC//u//+N7I3v16oWffvoJ586dQ2xsLNzc3CRaXH369EGzZs3g6uqKf//9F+fOnePPVWzChAmoX78+hg0bhnPnziExMRHR0dGYN28enjx5ovA109fXx5QpU7Bo0SKcOnUKsbGxcHd3l7gu9vb2mDBhAlxdXREWFobExERcvXoV69atw99//w2gaMW+v//+G4GBgUhISMD333+PY8eOyW1N9+nTB507d8bw4cMRERGBx48f4+LFi/j666/5ZLpixQqEhobC19cXd+7cQVxcHA4cOICvv/4aABAcHIw9e/YgNjYWjx49wk8//QQdHR1YW1vLPKe8+uzs7FBQUIBt27bx9X333XcSdXh6eiIiIgKJiYm4fv06Tp8+jU8++QRA2X93Ukmq6mFldeHm5iZzfRA3NzfGWNF6Gp9++inT0tJiZmZmbPHixSw/P5//vouLC99hUGzYsGH89xljLC4ujnXp0oVpaWmx5s2bs8OHDzMA7Pjx44wx6Q4UxhhbtWoVMzMzYxzH8XW9efOGjR49mhkaGjJLS0sWHBws0YHCGGPx8fGsa9euTEtLi9nb27Pjx49LPXhPTU1lrq6urH79+kwkEjFbW1s2bdo09ubNG5nXqKwH+2/fvmUTJ05kurq6zNTUlK1fv17qeuTl5bEVK1awxo0bM01NTWZmZsZGjBjBbt26xZfZtWsXa9iwIdPR0WHDhw9nfn5+zMzMjD/u4+PDWrduLRVXZmYmmzNnDrOwsGCamprM0tKSTZgwgSUnJ/Nljh8/zpydnZmOjg4zNDRkHTt2ZLt27WKMMRYeHs6cnJyYoaEh09PTY506dWInT56UeQ0UqY8xxgIDA5m5uTnT0dFh/fv3Z6GhoRLX7quvvmJ2dnZMJBKxBg0asEmTJrGMjAz++7L+7sVK/x1JxaI1UKrAhQsX0LVrVzx48AB2dnZVHY5cUVFR6NmzJ169elUpg66nTZuGe/fu4dy5cyo/V03CcRzCw8Or/I2pjxXdJleC8PBwREZG4vHjxzh58iS+/PJLdOnSpUYkwpIaNWqEcePGVXi9GzduxL///osHDx5g27ZtCAkJgZubW4Wfp6aaMWMG3+FEVIdahpUgNDQUq1evRkpKCurXr48+ffpg06ZNUsNCqqvs7Gw8ffoUQNGzwop+dXD06NGIiorC27dvYWtrizlz5mDGjBkVeo6aLD09HZmZmQCKhm2V7lknFYOSISGEgG6TCSEEACVDQggBQMmQEEIAUDIkhBAAlAwJIQQAJUNCCAFAyZAQQgBQMiSEEADA/wO+OxZ6RKZ9VwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ss['n'].plot(figsize=(3,3));" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Multi models" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "fn = \"../tests/testdata/SW/HKZN_local_2017_DutchCoast.dfsu\"\n", "mr1 = ms.model_result(fn, name=\"SW_1\", item=0)\n", "fn = \"../tests/testdata/SW/HKZN_local_2017_DutchCoast_v2.dfsu\"\n", "mr2 = ms.model_result(fn, name=\"SW_2\", item=0)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "o1 = ms.PointObservation(\"../tests/testdata/SW/HKNA_Hm0.dfs0\", item=0, x=4.2420, y=52.6887, name=\"HKNA\")\n", "o2 = ms.PointObservation(\"../tests/testdata/SW/eur_Hm0.dfs0\", item=0, x=3.2760, y=51.9990, name=\"EPL\")\n", "o3 = ms.TrackObservation(\"../tests/testdata/SW/Alti_c2_Dutch.dfs0\", item=3, name=\"c2\")\n", "cc = ms.match(obs=[o1, o2, o3], mod=[mr1, mr2])" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## gridded_skill\n", "Group by model" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "ss = cc.gridded_skill(by='model', bins=6)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['SW_1', 'SW_2']" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ss.mod_names" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAp0AAAEiCAYAAAC7oWkDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAyR0lEQVR4nO3deXQUZb7/8U91Z2MgiSSQhGgSArIv4g9ZgnIBRWQRlTuggBNQ0VEBBUFU1DmChwvXuYq4gYAMi8MA1wWvOgzbZXUgyKocfgxGBQlLwDCasGbprt8f+dFjS7ornaRTHfJ+nVPnkNr6W+3xySfPU/WUYZqmKQAAACCIHHYXAAAAgKsfoRMAAABBR+gEAABA0BE6AQAAEHSETgAAAAQdoRMAAABBR+gEAABA0BE6AQAAEHSETgAAAAQdoRPVpmfPnho/fny591+0aJGuueaaoNUDoHrRBgC1G6ETtdrevXt15513KiEhQVFRUWrcuLHuu+8+5eXlSZIaNWqkV155xeuYZ599VoZh6H//93+91t92220aPny45WeePHlSw4cPV4sWLeRwOAL6JQygatnRBnz88ce6/fbb1bBhQ8XExCgjI0Nr1qypuosCQhShE7XW6dOn1bt3bzVo0EBr1qzRwYMH9ac//UmNGjXShQsXJJX2zGzcuNHruE2bNiklJcVrfVFRkbZv365evXpZfm5hYaEaNmyoF154QTfccEPVXhSAcrOrDdiyZYtuv/12rVq1Srt371avXr00cOBA7d27t2ovEAg1Jmq1Hj16mGPHjjXHjRtnXnPNNWZCQoI5d+5c89y5c+YDDzxg1qtXz2zSpIm5atUqr+M2bdpkdurUyYyIiDCTkpLMZ5991iwuLvZsP3funJmZmWnWrVvXTEpKMl999VWzR48e5rhx4zz7FBYWmpMmTTKTk5PN3/zmN2bnzp3NjRs3erYvXLjQjI2NDdq1r1y50gwLC/Oq+9fmzp1r1qtXz7NPQUGBGR4ebr7zzjvmzTff7Nlvy5YtpiQzOzs7oBp+/Z0A1Y02wN424LLWrVubU6dOrdCxQE1BTye0ePFiNWjQQF9++aWeeOIJPf744xoyZIi6deumPXv26I477lBmZqbnL//jx4+rf//+6tSpk7766ivNmTNHCxYs0LRp0zznnDRpkjZu3KiVK1dq7dq12rRpk3bv3u31uQ8++KD+/ve/a/ny5fr66681ZMgQ9e3bV9nZ2eWuvU2bNqpXr57PpU2bNj6PTUpKUklJiVauXCnTNMvcp1evXjp37px27twpSdq6dauaN2+uwYMHa+fOnZ7vZOPGjbruuut0/fXXl7t2IFTQBtjbBrjdbp09e1ZxcXEBHwvUKHanXtirR48e5i233OL5uaSkxKxbt66ZmZnpWXfy5ElTkrl9+3bTNE3z+eefN1u0aGG63W7PPu+8845Zr1490+VymWfPnjUjIiLM5cuXe7afOXPGrFOnjqeX49tvvzUNwzCPHz/uVc9tt91mTp482TTN8vVyHDlyxMzOzva5HDlyxO/xzz//vBkWFmbGxcWZffv2Nf/4xz+aubm5Xvtce+215vTp003TNM1JkyaZo0ePNk3TNFu2bGmuXbvWNE3T7NWrl9d3Vl70dMJutAH2tgGmaZp//OMfzbi4OPPUqVMVOh6oKejphNq3b+/5t9PpVHx8vNq1a+dZl5iYKKn0/idJOnjwoDIyMmQYhmefm2++WefOndOxY8f03XffqaioSBkZGZ7tcXFxatGihefnPXv2yDRNNW/e3KtXYvPmzfruu+/KXXtaWpquv/56n0taWprf4//jP/5Dubm5evfdd9W6dWu9++67atmypfbv3+/Zp2fPntq0aZOk0nu5evbsKUnq0aOHNm3apMLCQmVlZenWW28td91AKKENsK8NWLZsmaZMmaIVK1YoISEh4OOBmiTM7gJgv/DwcK+fDcPwWnf5F4vb7ZYkmabp9cvm8rrL+5o+hql+ye12y+l0avfu3XI6nV7b6tWrV+7a27Rpox9++MHn9rS0NB04cMDvOeLj4zVkyBANGTJEM2bM0I033qhXX31VixcvllQ6vDZu3DidOXNGe/fu1b/9279JKv2F89Zbb6lPnz66ePFiuR4gAEIRbYA9bcCKFSs0atQoffDBB+rdu3dAxwI1EaETAWvdurU++ugjr18827ZtU3R0tK699lrVr19f4eHhysrKUmpqqiTpp59+0jfffKMePXpIkm688Ua5XC6dPn1a3bt3r3Atq1atUnFxsc/tv/5laiUiIkJNmzbV+fPnPet69eql8+fPa+bMmWrWrJmn16dHjx4aOXKk/vrXvyo9Pd2yRwW4WtAGVL4NWLZsmR566CEtW7ZMAwYMCKhGoKYidCJgo0eP1qxZs/TEE09o7NixOnTokF566SVNmDBBDodD9erV06hRozRp0iTFx8crMTFRL7zwghyOf93N0bx5c91///0aMWKEXnvtNd14443Ky8vThg0b1K5dO/Xv379ctVQm6H3++edavny5hg4dqubNm8s0TX322WdatWqVFi5c6NmvSZMmSk1N1VtvvaX777/fsz45OVlpaWl69913NWTIkIA+e9++fZKkc+fO6ccff9S+ffsUERGh1q1bV/h6gOpCG1Cqom3AsmXLNGLECL3xxhvq2rWrcnNzJUl16tRRbGxsha8HCHWETgTs2muv1apVqzRp0iTdcMMNiouL06hRo/Tiiy969vmv//ovnTt3TnfddZeio6M1ceJE5efne51n4cKFmjZtmiZOnKjjx48rPj5eGRkZ5f5lU1mtW7fWb37zG02cOFE5OTmKjIxUs2bN9N577ykzM9Nr3169emnx4sWee7ku69GjhxYsWBDwsNqNN97o+ffu3bv1l7/8RWlpaTpy5EhFLweoNrQB/1KRNmDu3LkqKSnRmDFjNGbMGM/6kSNHatGiRZW5JCCkGWZ5br4BAAAAKoGn1wEAABB0hE6givmbrHrp0qV2lwcgyGgDgLIxvA5UsR9++MHn07SJiYmKjo6u5ooAVCfaAKBshE4AAAAEHcPrAAAACDpCJwAAAILuqp+n0+1268SJE4qOjr7itW0AQp9pmjp79qySk5O9JhcvL9oAoGarbBtQFS5duqSioiLL/SIiIhQVFVUNFdVMV33oPHHihFJSUuwuA0Al5eTk6Lrrrgv4ONoA4OpQ0Tagsi5duqT4OvV0QS7LfZOSknT48GGCpw9Xfei8/JTgnskPKDoqwuZqylZ/+BjrnWwWntDY7hJQSxUUFCglJaXCT/xePm5SZGNFGqF5R9HE7fPtLsGSo8lNdpeAWqqybUBlFRUV6YJcGqFrFeHnrsQiubUk97iKiooInT5c9aHz8nBadFREyIbOmBowfUZ4TIzdJaCWq+jQ+OXjIg2HokI0dMZE17O7BEsO2gDYzO7bY+oYTkX4aUOcpiExH5BfV33oBAAAqCyHITn95F6HROi0QOgEAACwEOEwFOGnt9U0DcldjQXVQIROAAAAC07DkNNP6HSK2TGsEDoBAAAsOC2G153VV0qNRegEAACwQE9n5RE6AQAALIQZhsL9hE4XodMSoRMAAMACw+uVR+gEAACwUBo6/Q2vw4qtMyVPmTJFhmF4LUlJSZ7tDzzwwBXbu3btamPFAACgNrrc0+lvgX+293S2adNG69ev9/zsdHr/rdC3b18tXLjQ83NERGi+VQgAAFy9wi3m6XSZpE4rtofOsLAwr97NX4uMjPS7HQAAINh4er3ybH8RcXZ2tpKTk5Wenq6hQ4fq+++/99q+adMmJSQkqHnz5nrkkUd0+vRpv+crLCxUQUGB1wKg9qANABAMDK9Xnq2hs0uXLlqyZInWrFmj+fPnKzc3V926ddOZM2ckSf369dPSpUu1YcMGvfbaa9q5c6duvfVWFRYW+jznjBkzFBsb61lSUlKq63IAhADaAADBQOisPMM0zZB5Pf358+fVtGlTPfPMM5owYcIV20+ePKm0tDQtX75c//7v/17mOQoLC71CaUFBgVJSUpQ99feKjgrN+0HjRl55raEmPDHd7hJQSxUUFCg2Nlb5+fmKiYmx3N9XG/BiVBNFGbYP7pRp8v6ldpdgydG0s90loJYKtA0I1ufPrt9MdQzfz6hfNF0a/VO2bXXWBLbf0/lLdevWVbt27ZSdnV3m9kaNGiktLc3ndqn0HtDIyMhglQggxNEGAAiGcMNQhMN3d2aJm65OKyH1Z39hYaEOHjyoRo0albn9zJkzysnJ8bkdAAAgGBherzxbQ+fTTz+tzZs36/Dhw9qxY4cGDx6sgoICjRw5UufOndPTTz+t7du368iRI9q0aZMGDhyoBg0aaNCgQXaWDQAAapnLT6/7W+CfrcPrx44d07Bhw5SXl6eGDRuqa9euysrKUlpami5evKj9+/dryZIl+vnnn9WoUSP16tVLK1asUHR0tJ1lAwCAWsbyNZhkTku2hs7ly5f73FanTh2tWbOmGqsBAAAoW7jDoXCH7wHicIXMc9khK6QeJAIAAAhFhtOQ4edBIoPhdUuETgAAAAsOpyGHn9DpIHRaInQCAABYcTpk+Blel8HwuhVCJwAAgAVnuENOp+/Q6XQROq0QOgEAACwYDv89nUbovOAxZIXU5PAAAAChyOE0LJeKmD17ttLT0xUVFaWOHTtq69at5Tru73//u8LCwtShQ4cKfa4dCJ0AAAAWDKdhuQRqxYoVGj9+vF544QXt3btX3bt3V79+/XT06FG/x+Xn52vEiBG67bbbKno5tiB0AgAAWHBGOOSMcPpZAo9UM2fO1KhRo/Twww+rVatWmjVrllJSUjRnzhy/xz366KMaPny4MjIyKno5tiB0AgAAWDCM0nk6fS7/f8qkgoICr6WwsLDM8xUVFWn37t3q06eP1/o+ffpo27ZtPutYuHChvvvuO7300ktVd3HVhNAJAABgweF0WC6SlJKSotjYWM8yY8aMMs+Xl5cnl8ulxMREr/WJiYnKzc0t85js7Gw999xzWrp0qcLCat6z4DWvYgAAgGpmdd+mYZZuy8nJUUxMjGd9ZGSk//P+alJ50zTLfLuRy+XS8OHDNXXqVDVv3jyQ0kMGoRMAAMBCeUNnTEyMV+j0pUGDBnI6nVf0ap4+ffqK3k9JOnv2rHbt2qW9e/dq7NixkiS32y3TNBUWFqa1a9fq1ltvDeSSqh2hEwAAwIIjwilnmNP3dkdg83RGRESoY8eOWrdunQYNGuRZv27dOt19991X7B8TE6P9+/d7rZs9e7Y2bNigDz/8UOnp6QF9vh1qTejM//6EXBHhdpdRpvjiC3aXAFz1fi52KTJEX1NnFpf9oAGA0OEwqv7d6xMmTFBmZqZuuukmZWRkaN68eTp69Kgee+wxSdLkyZN1/PhxLVmyRA6HQ23btvU6PiEhQVFRUVesD1W1JnQCAABUlOF0yPDzGkzDHfiz2ffdd5/OnDmjl19+WSdPnlTbtm21atUqpaWlSZJOnjxpOWdnTULoBAAAsGD11iGHu2JvJBo9erRGjx5d5rZFixb5PXbKlCmaMmVKhT7XDoROAAAAC44Ipxzhfu7pDNHbd0IJoRMAAMCCwymLns5qLKaGInQCAABYuPzmIX/b4R+hEwAAwILD8a+3DpW53cVLHq0QOgEAACw4IpxyRPi5p1OMr1shdAIAAFgwHA4ZDj9TJvnZhlKETgAAAAsOp8Xwup9tKEXoBAAAsGIxObwInZYInQAAABYMh8UbiRhet0ToBAAAsOAID5MjPNz3djeTw1shdAIAAFiwfPc6w+uWCJ0AAAAWHA6HHH6G0P1tQylCJwAAgAV6OiuP0AkAAGDBERYmR7jv2ORwMTm8FUInAACABXo6K4/QCQAAYIHQWXmETgAAAAu8kajyCJ0AAAAWHOFOi3s6XdVYTc1E6AQAALDA8HrlEToBAAAsGIbD76suDYPQaYXQCQAAYMFwOuVwOv1uh3+ETgAAAAuOiDA5IrinszIInQAAABYMh8XwOq/BtEToBAAAsMCDRJVn6zc0ZcoUGYbhtSQlJZW576OPPirDMDRr1qzqLRIAANR6hsPwBM8yF4dhd4khz/ZY3qZNG508edKz7N+//4p9PvnkE+3YsUPJyck2VAgAAGq7y8Pr/paKmD17ttLT0xUVFaWOHTtq69atPvf9+OOPdfvtt6thw4aKiYlRRkaG1qxZU9FLqna2h86wsDAlJSV5loYNG3ptP378uMaOHaulS5cqPDzcpioBAEBtZoRFWC6BWrFihcaPH68XXnhBe/fuVffu3dWvXz8dPXq0zP23bNmi22+/XatWrdLu3bvVq1cvDRw4UHv37q3s5VUL20Nndna2kpOTlZ6erqFDh+r777/3bHO73crMzNSkSZPUpk0bG6sEAAC1msNhvQRo5syZGjVqlB5++GG1atVKs2bNUkpKiubMmVPm/rNmzdIzzzyjTp06qVmzZpo+fbqaNWumzz77rLJXVy1sfZCoS5cuWrJkiZo3b65Tp05p2rRp6tatmw4cOKD4+Hi98sorCgsL05NPPlnucxYWFqqwsNDzc0FBQTBKBxCiaAMABIPhdPqdizPQeTqLioq0e/duPffcc17r+/Tpo23btpXrHG63W2fPnlVcXFxAn20XW0Nnv379PP9u166dMjIy1LRpUy1evFg9evTQG2+8oT179sgwyn9z7owZMzR16tRglAugBqANABAUDmfp4m+7rvxDNzIyUpGRkVfsnpeXJ5fLpcTERK/1iYmJys3NLVdJr732ms6fP6977723XPvbzfbh9V+qW7eu2rVrp+zsbG3dulWnT59WamqqwsLCFBYWph9++EETJ05U48aNfZ5j8uTJys/P9yw5OTnVdwEAbEcbACAYjLAwGWHhfpbSfryUlBTFxsZ6lhkzZvg/76861kzTLFdn27JlyzRlyhStWLFCCQkJFb+wahRS83QWFhbq4MGD6t69uzIzM9W7d2+v7XfccYcyMzP14IMP+jyHr78ozvwjT5fCQupyPQrcof9U/q6sI3aX4NfjXRvbXQJChK824LzLrWIb6imP7Po32F2CpZ17jtldgl+Z/+c6u0vA1c6w6Ok0Srfl5OQoJibGs7qs9kiSGjRoIKfTeUWv5unTp6/o/fy1FStWaNSoUfrggw+uyEqhzNYU9vTTT2vgwIFKTU3V6dOnNW3aNBUUFGjkyJGKj49XfHy81/7h4eFKSkpSixYtbKoYAADUSuUcXo+JifEKnb5ERESoY8eOWrdunQYNGuRZv27dOt19990+j1u2bJkeeughLVu2TAMGDCh//SHA1tB57NgxDRs2THl5eWrYsKG6du2qrKwspaWl2VkWAACAl2C8BnPChAnKzMzUTTfdpIyMDM2bN09Hjx7VY489Jqn0dqHjx49ryZIlkkoD54gRI/TGG2+oa9eunl7SOnXqKDY2tgJXVb1sDZ3Lly8PaP8jR44EpxAAAAB/wsIlf3NxhpUEfMr77rtPZ86c0csvv6yTJ0+qbdu2WrVqlafz7eTJk15zds6dO1clJSUaM2aMxowZ41k/cuRILVq0KODPr26heZMjAABACKnqKZMuGz16tEaPHl3mtl8HyU2bNlXoM0IFoRMAAMCK1QTwFXwNZm1C6AQAALBSzgeJ4BuhEwAAwILhLJ2P0992+EfoBAAAsGA4nDL89Gb624ZShE4AAAArDofF8Dr3dFohdAIAAFjhQaJKI3QCAABYCNaUSbUJoRMAAMBKWITF5PDF1VdLDUXoBAAAsBCM12DWNoROAAAAK4bFPJ0Gw+tWCJ0AAABWDEMy/PRmGkb11VJDEToBAAAsmI4wmQ7fscnfNpTiGwIAALBiOCx6Ormn0wqhEwAAwIph+B9CZ3jdEqETAADACpPDVxqhEwAAwAL3dFYe3xAAAIAV7umsNEInAACAFUJnpRE6AQAALJiGIdNPsDR5kMgSsRwAAMCKw2m9XGVWr16tL774wvPzO++8ow4dOmj48OH66aefAj4foRMAAMDK5eF1f8tVZtKkSSooKJAk7d+/XxMnTlT//v31/fffa8KECQGfj+F1AAAAC6bhsBhev/pC5+HDh9W6dWtJ0kcffaQ777xT06dP1549e9S/f/+Az3f1fUMAAABVzXD8a67OsparMHRGRETowoULkqT169erT58+kqS4uDhPD2gg6OkEAACwUgufXr/llls0YcIE3Xzzzfryyy+1YsUKSdI333yj6667LuDzXX3fEAAAQFVzhFkvV5m3335bYWFh+vDDDzVnzhxde+21kqS//e1v6tu3b8Dnu/q+IQAAgCpWG6dMSk1N1eeff37F+tdff71C56OnEwAAwEqQnl6fPXu20tPTFRUVpY4dO2rr1q1+99+8ebM6duyoqKgoNWnSRO+++26FPjdQFy9eVEFBgdcSqIC/oQceeEBbtmwJ+IMAAABqLMOwXgK0YsUKjR8/Xi+88IL27t2r7t27q1+/fjp69GiZ+x8+fFj9+/dX9+7dtXfvXj3//PN68skn9dFHH1X26sp0/vx5jR07VgkJCapXr57q16/vtQQq4NB59uxZ9enTR82aNdP06dN1/PjxgD8UAACgJjEdYZZLoGbOnKlRo0bp4YcfVqtWrTRr1iylpKRozpw5Ze7/7rvvKjU1VbNmzVKrVq308MMP66GHHtKrr75a2csr0zPPPKMNGzZo9uzZioyM1HvvvaepU6cqOTlZS5YsCfh8AYfOjz76SMePH9fYsWP1wQcfqHHjxurXr58+/PBDFRcXB1wAAABAyCvn8Pqvh6ALCwvLPF1RUZF2797tmYbosj59+mjbtm1lHrN9+/Yr9r/jjju0a9euoGSwzz77TLNnz9bgwYMVFham7t2768UXX9T06dO1dOnSgM9XoQeJ4uPjNW7cOI0bN0579+7Vn/70J2VmZqpevXr63e9+p9GjR6tZs2YVOXXQhNUJU1hYaD431Sk18C7q6jZhxT67S/Dr8a6N7S6hRvvL3mN2l+DThXNnq+Q8bWMiVccIzdfUtUyMsbsES48v22d3CX5l/p/Ap2/Bv/zf3MDvz6su586GRm2lDxL5HkK/vC0lJcVr/UsvvaQpU6ZcsX9eXp5cLpcSExO91icmJio3N7fMz8jNzS1z/5KSEuXl5alRo0bluZRy++c//6n09HRJUkxMjP75z39KKp1K6fHHHw/4fJV6kOjkyZNau3at1q5dK6fTqf79++vAgQNq3bp1hZ9sAgAACDWmab1IUk5OjvLz8z3L5MmT/Z7X+FWQNU3zinVW+5e1vio0adJER44ckSS1bt1a//3f/y2ptAf0mmuuCfh8AXf9FRcX69NPP9XChQu1du1atW/fXk899ZTuv/9+RUdHS5KWL1+uxx9/XE899VTABQEAAIQal2nKdTlZ+tgulfYIxsRYj140aNBATqfzil7N06dPX9GbeVlSUlKZ+4eFhSk+Pt7yMwP14IMP6quvvlKPHj00efJkDRgwQG+99ZZKSko0c+bMgM8XcOhs1KiR3G63hg0bpi+//FIdOnS4Yp877rijQgkYAAAgFLnN0sXf9kBERESoY8eOWrdunQYNGuRZv27dOt19991lHpORkaHPPvvMa93atWt10003KTw8PLACyuGXnYe9evXSP/7xD+3atUtNmzbVDTfcEPD5Ag6dr7/+uoYMGaKoqCif+9SvX1+HDx8OuBgAAIBQZJqmZyjb1/ZATZgwQZmZmbrpppuUkZGhefPm6ejRo3rsscckSZMnT9bx48c9T4o/9thjevvttzVhwgQ98sgj2r59uxYsWKBly5ZV7KIClJqaqtTU1AofH3DozMzMrPCHAQAA1ERV3dMpSffdd5/OnDmjl19+WSdPnlTbtm21atUqpaWlSSp9duaXc3amp6dr1apVeuqpp/TOO+8oOTlZb775pn77298G/uE+vPnmm/r973+vqKgovfnmm373ffLJJwM6d2g+zg0AABBiKpArLY0ePVqjR48uc9uiRYuuWNejRw/t2bMnCJWUev3113X//fcrKirK70PhhmEQOgEAAKqay23K5ac709+2muSXt0f+8t9V8ZQ8714HAACw4C7HcjVasGCB2rZtq6ioKEVFRalt27Z67733KnQuW0PnlClTZBiG15KUlOS1vWXLlqpbt67q16+v3r17a8eOHTZWDAAAaqPyztN5NfnDH/6gcePGaeDAgfrggw/0wQcfaODAgXrqqaf04osvBnw+24fX27Rpo/Xr13t+djr/9caQ5s2b6+2331aTJk108eJFvf766+rTp4++/fZbNWzY0I5yAQBALRSMB4lC3Zw5czR//nwNGzbMs+6uu+5S+/bt9cQTT2jatGkBnc/20BkWFubVu/lLw4cP9/p55syZWrBggb7++mvddttt1VEeAABAuSeHv5q4XC7ddNNNV6zv2LGjSkpKAj6f7fd0ZmdnKzk5Wenp6Ro6dKi+//77MvcrKirSvHnzFBsbW6EJSQEAACrKlMXwut0FBsHvfvc7zZkz54r18+bN0/333x/w+Wzt6ezSpYuWLFmi5s2b69SpU5o2bZq6deumAwcOeF7n9Pnnn2vo0KG6cOGCGjVqpHXr1qlBgwY+z1lYWKjCwkLPzwUFBUG/DgChgzYAQDC4TVNuP72Z/rbVJBMmTPD82zAMvffee1q7dq26du0qScrKylJOTo5GjBgR8LltDZ39+vXz/Ltdu3bKyMhQ06ZNtXjxYs9F9+rVS/v27VNeXp7mz5+ve++9Vzt27FBCQkKZ55wxY4amTp1aLfUDCD20AQCCwZT/3syrI3JKe/fu9fq5Y8eOkqTvvvtOktSwYUM1bNhQBw4cCPjctt/T+Ut169ZVu3btlJ2d7bXu+uuv1/XXX6+uXbuqWbNmWrBggSZPnlzmOSZPnuyV0gsKCpSSkhL02gGEBtoAAMHgdksuP/Miua+SOZM2btwYtHOHVOgsLCzUwYMH1b17d5/7mKbpNXT2a5GRkYqMjAxGeQBqANoAAMHglim3n/5Mf9tQytYHiZ5++mlt3rxZhw8f1o4dOzR48GAVFBRo5MiROn/+vJ5//nllZWXphx9+0J49e/Twww/r2LFjGjJkiJ1lAwCAWqY2ztNZ1Wzt6Tx27JiGDRumvLw8NWzYUF27dlVWVpbS0tJ06dIl/eMf/9DixYuVl5en+Ph4derUSVu3blWbNm3sLBsAANQytXGezqpma+hcvny5z21RUVH6+OOPq7EaAACAstXGeTqrWkjd0wkAABCKrIbQyZzWCJ0AAAAWass8ncFE6AQAALDgspgyyd82lCJ0AgAAWKCns/IInQAAABZK3KaK/cwAX8Lj65YInQAAABYYXq88QicAAIAFhtcrj9AJAABggXk6K4/QCQAAYKHYbarY5TtYFnNPpyVCJwAAgAXTYnjdpKfTEqETAADAgsssXfxth3+ETgAAAAs8SFR5DrsLAAAACHXFLrflEiw//fSTMjMzFRsbq9jYWGVmZurnn3/2XWtxsZ599lm1a9dOdevWVXJyskaMGKETJ04ErcbyIHQCAABYuDy87m8JluHDh2vfvn1avXq1Vq9erX379ikzM9Pn/hcuXNCePXv0hz/8QXv27NHHH3+sb775RnfddVfwiiyHWjO8ftP//FUxMTF2l1Fj5R752e4SarSj/zxndwl+Db/xOrtL8KmgoECPVMF5Hs7ZTRtQCceyz9hdQo128dIlu0vwq3VS6P6/UfAbuysoZdfw+sGDB7V69WplZWWpS5cukqT58+crIyNDhw4dUosWLa44JjY2VuvWrfNa99Zbb6lz5846evSoUlNTg1KrlVoTOgEAACrK7Tbl9jMt0uVtBQUFXusjIyMVGRlZ4c/dvn27YmNjPYFTkrp27arY2Fht27atzNBZlvz8fBmGoWuuuabCtVQWw+sAAAAWSt+97nu5/O71lJQUz72XsbGxmjFjRqU+Nzc3VwkJCVesT0hIUG5ubrnOcenSJT333HMaPny4rSM+9HQCAABYKO8biXJycryCna9ezilTpmjq1Kl+P3Pnzp2SJMMwrthmmmaZ63+tuLhYQ4cOldvt1uzZsy33DyZCJwAAgIXyDq/HxMSUqzdx7NixGjp0qN99GjdurK+//lqnTp26YtuPP/6oxMREv8cXFxfr3nvv1eHDh7Vhwwbb72sndAIAAFhwyWJy+ADP16BBAzVo0MByv4yMDOXn5+vLL79U586dJUk7duxQfn6+unXr5vO4y4EzOztbGzduVHx8fIAVVj3u6QQAALBw+el1f0swtGrVSn379tUjjzyirKwsZWVl6ZFHHtGdd97p9RBRy5YttXLlSklSSUmJBg8erF27dmnp0qVyuVzKzc1Vbm6uioqKglJnedDTCQAAYKHY5ZbTzwTwwZwcfunSpXryySfVp08fSdJdd92lt99+22ufQ4cOKT8/X5J07Ngxffrpp5KkDh06eO23ceNG9ezZM2i1+kPoBAAAsOBym3L5uafT37bKiouL05///Ge/+5i/6Glt3Lix18+hgtAJAABgwc7QebUgdAIAAFhwuf0HyyCOrl81CJ0AAAAWikrccpT4TpZFfrahFKETAADAgttieN3fHJ4oRegEAACw4DIt7ukMwQd3Qg2hEwAAwAIPElUeoRMAAMBCYYlb8nPfZiH3dFoidAIAAFigp7PyCJ0AAAAWeJCo8gidAAAAFlym6fdhIR4kskboBAAAsFBkcU8n83RaI3QCAABY4J7OyiN0AgAAWHCZbrncvnszXSY9nVYInQAAABZ4kKjyCJ0AAAAWXG5TDobXK4XQCQAAYKGwxJTbz8NCxSWETiuETgAAAAv0dFYeoRMAAMACobPyHHZ++JQpU2QYhteSlJQkSSouLtazzz6rdu3aqW7dukpOTtaIESN04sQJO0sGAAC10OUHiXwtPEhkzfaezjZt2mj9+vWen51OpyTpwoUL2rNnj/7whz/ohhtu0E8//aTx48frrrvu0q5du+wqFwAA1EIlLv+Tw5e4mDLJiu2hMywszNO7+UuxsbFat26d17q33npLnTt31tGjR5WamlpdJQIAgFrObdGbSU+nNdtDZ3Z2tpKTkxUZGakuXbpo+vTpatKkSZn75ufnyzAMXXPNNT7PV1hYqMLCQs/PBQUFVV0ygBBGGwAgGEzTlOnn/er+tqGUrfd0dunSRUuWLNGaNWs0f/585ebmqlu3bjpz5swV+166dEnPPfechg8frpiYGJ/nnDFjhmJjYz1LSkpKMC8BQIihDQAQDKbbtFzgn62hs1+/fvrtb3+rdu3aqXfv3vrrX/8qSVq8eLHXfsXFxRo6dKjcbrdmz57t95yTJ09Wfn6+Z8nJyQla/QBCD20AgGBwlZiWC/yzfXj9l+rWrat27dopOzvbs664uFj33nuvDh8+rA0bNvjt5ZSkyMhIRUZGBrtUACGKNgBAMDC8Xnm29nT+WmFhoQ4ePKhGjRpJ+lfgzM7O1vr16xUfH29zhQAAoDa6/CCRvyVYfvrpJ2VmZnpuG8rMzNTPP/9c7uMfffRRGYahWbNmBa3G8rC1p/Ppp5/WwIEDlZqaqtOnT2vatGkqKCjQyJEjVVJSosGDB2vPnj36/PPP5XK5lJubK0mKi4tTRESEnaVXqYufv2N3CZbuu/NWu0uo0VLj6tldAkLYhY9etbsES4P797e7hBqtTlSU3SWgkqzu2wzmPZ3Dhw/XsWPHtHr1aknS73//e2VmZuqzzz6zPPaTTz7Rjh07lJycHLT6ysvW0Hns2DENGzZMeXl5atiwobp27aqsrCylpaXpyJEj+vTTTyVJHTp08Dpu48aN6tmzZ/UXDAAAaierh4WCFDoPHjyo1atXKysrS126dJEkzZ8/XxkZGTp06JBatGjh89jjx49r7NixWrNmjQYMGBCU+gJha+hcvny5z22NGzfm/ggAABASXG635GcCeJc7OJPDb9++XbGxsZ7AKUldu3ZVbGystm3b5jN0ut1uZWZmatKkSWrTpk1QagtUSD1IBAAAEIrKO7z+67mBK/twY25urhISEq5Yn5CQ4LntsCyvvPKKwsLC9OSTT1b4s6taSD1IBAAAEIrcbquHiUr3S0lJ8ZoreMaMGWWeb8qUKTIMw+9y+bXfhmFccbxpmmWul6Tdu3frjTfe0KJFi3zuYwd6OgEAACyUd8qknJwcr+kdffVyjh07VkOHDvX7mY0bN9bXX3+tU6dOXbHtxx9/VGJiYpnHbd26VadPn/Z6ZbjL5dLEiRM1a9YsHTlyxO/nBguhEwAAwIKrxJScvkPn5cnhY2JiLOcUl6QGDRqoQYMGlvtlZGQoPz9fX375pTp37ixJ2rFjh/Lz89WtW7cyj8nMzFTv3r291t1xxx3KzMzUgw8+aPmZwULoBAAAsGDXlEmtWrVS37599cgjj2ju3LmSSqdMuvPOO70eImrZsqVmzJihQYMGKT4+/oq5zcPDw5WUlOT3afdg455OAAAAC3a+e33p0qVq166d+vTpoz59+qh9+/Z6//33vfY5dOiQ8vPzg1ZDVaCnEwAAwILbNGX4uafTHcRpHuPi4vTnP//Z7z5W00zadR/nLxE6AQAALLhdbhklvufidPuZwxOlCJ0AAAAWTIv3qwdzeP1qQegEAACwUN4pk+AboRMAAMCCXU+vX00InQAAABZcJSUyHSU+t7tLfG9DKUInAACABdPtkul2+d0O/widAAAAFky32yJ08vS6FUInAACABdPlkunyEzr9bEMpQicAAIAF07QYXjcJnVYInQAAABbcJUWS4fS/HX4ROgEAACzwIFHlEToBAAAs8CBR5RE6AQAALLjdLslP6HTT02mJ0AkAAGCh9J5Oh//t8IvQCQAAYMXlkunw05vJlEmWCJ0AAAAWTNP/8DpTJlkjdAIAAFgw3W7/oZMHiSwROgEAACy4S4plyPc9nWZJcTVWUzMROgEAACyYFk+vM0+nNUInAACABUJn5RE6AQAALLjdLhmEzkohdAIAAFhwlxTLMA2f200X93RaIXQCAABYMN0uyaCnszIInQAAABYInZV31YdO0zQlSQUFBTZX4tvFCxftLsFSofOc3SX4Fcr/fVE5l//bXv5/OVA1oQ24cOGS3SVYKnTTBsAelW0DqopZfMl/sGR43ZJh2v1fMciOHTumlJQUu8sAUEk5OTm67rrrAj6ONgC4OlS0DaisS5cuKT09Xbm5uZb7JiUl6fDhw4qKiqqGymqeqz50ut1unThxQtHR0TIM3zcAFxQUKCUlRTk5OYqJianGCkMf341/fD++VcV3Y5qmzp49q+TkZDkcvidm9oU2oPL4bvzj+/EtFNqAqnDp0iUVFRVZ7hcREUHg9OOqH153OBwB/WUUExNDo+ED341/fD++Vfa7iY2NrfCxtAFVh+/GP74f3+xsA6pCVFQUYbIK2PMnAwAAAGoVQicAAACCjtD5/0VGRuqll15SZGSk3aWEHL4b//h+fKtJ301NqrW68d34x/fjG98Nfumqf5AIAAAA9qOnEwAAAEFH6AQAAEDQEToBAAAQdLU+dG7ZskUDBw5UcnKyDMPQJ598YndJIWPGjBnq1KmToqOjlZCQoHvuuUeHDh2yu6yQMGfOHLVv394z91xGRob+9re/2V1WSJoxY4YMw9D48ePtLqVMtAG+0Qb4RhtQfqHeBqD61PrQef78ed1www16++237S4l5GzevFljxoxRVlaW1q1bp5KSEvXp00fnz5+3uzTbXXfddfrP//xP7dq1S7t27dKtt96qu+++WwcOHLC7tJCyc+dOzZs3T+3bt7e7FJ9oA3yjDfCNNqB8akIbgGpkwkOSuXLlSrvLCFmnT582JZmbN2+2u5SQVL9+ffO9996zu4yQcfbsWbNZs2bmunXrzB49epjjxo2zuyRLtAH+0Qb4RxvgrSa2AQiuWt/TifLLz8+XJMXFxdlcSWhxuVxavny5zp8/r4yMDLvLCRljxozRgAED1Lt3b7tLQRWhDSgbbUDZaAPwa1f9u9dRNUzT1IQJE3TLLbeobdu2dpcTEvbv36+MjAxdunRJ9erV08qVK9W6dWu7ywoJy5cv1549e7Rz5067S0EVoQ24Em2Ab7QBKAuhE+UyduxYff311/riiy/sLiVktGjRQvv27dPPP/+sjz76SCNHjtTmzZtr/S+dnJwcjRs3TmvXrlVUVJTd5aCK0AZciTagbLQB8IU3Ev2CYRhauXKl7rnnHrtLCSlPPPGEPvnkE23ZskXp6el2lxOyevfuraZNm2ru3Ll2l2KrTz75RIMGDZLT6fSsc7lcMgxDDodDhYWFXttCCW1A2WgDyoc2oFRNbgMQXPR0wifTNPXEE09o5cqV2rRpE79sLJimqcLCQrvLsN1tt92m/fv3e6178MEH1bJlSz377LP8sqlBaAMCQxtQijYAvtT60Hnu3Dl9++23np8PHz6sffv2KS4uTqmpqTZWZr8xY8boL3/5i/7nf/5H0dHRys3NlSTFxsaqTp06Nldnr+eff179+vVTSkqKzp49q+XLl2vTpk1avXq13aXZLjo6+op7/urWrav4+PiQvBeQNsA32gDfaAN8q2ltAKpPrQ+du3btUq9evTw/T5gwQZI0cuRILVq0yKaqQsOcOXMkST179vRav3DhQj3wwAPVX1AIOXXqlDIzM3Xy5EnFxsaqffv2Wr16tW6//Xa7S0OAaAN8ow3wjTYACBz3dAIAACDomKcTAAAAQUfoBAAAQNAROgEAABB0hE4AAAAEHaETAAAAQUfoBAAAQNAROgEAABB0hE4AAAAEHaETAAAAQUfoBAAAQNAROgEAABB0hE7USD/++KOSkpI0ffp0z7odO3YoIiJCa9eutbEyANWBNgCoeQzTNE27iwAqYtWqVbrnnnu0bds2tWzZUjfeeKMGDBigWbNm2V0agGpAGwDULIRO1GhjxozR+vXr1alTJ3311VfauXOnoqKi7C4LQDWhDQBqDkInarSLFy+qbdu2ysnJ0a5du9S+fXu7SwJQjWgDgJqDezpRo33//fc6ceKE3G63fvjhB7vLAVDNaAOAmoOeTtRYRUVF6ty5szp06KCWLVtq5syZ2r9/vxITE+0uDUA1oA0AahZCJ2qsSZMm6cMPP9RXX32levXqqVevXoqOjtbnn39ud2kAqgFtAFCzMLyOGmnTpk2aNWuW3n//fcXExMjhcOj999/XF198oTlz5thdHoAgow0Aah56OgEAABB09HQCAAAg6AidAAAACDpCJwAAAIKO0AkAAICgI3QCAAAg6AidAAAACDpCJwAAAIKO0AkAAICgI3QCAAAg6AidAAAACDpCJwAAAIKO0AkAAICg+3+B8xgtilcaQAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ss.bias.plot();" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHFCAYAAAAT5Oa6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABjS0lEQVR4nO3de1hU1foH8O8GBBSFBJGLIgKKiKIiXrikZihm5TUTtUDLS1bmLTMpTbySp1TUwtIspF95OalpZQewQDTREkGt1DAxUCGCkhEVEFi/PzjMcZwBZpgZhoHv53n2c5i1917zrr2rec/aa68lCSEEiIiIiJoxE0MHQERERGRoTIiIiIio2WNCRERERM0eEyIiIiJq9pgQERERUbPHhIiIiIiaPSZERERE1OwxISIiIqJmjwkRERERNXtMiIiM2NWrVyFJEmJjYzU+Nzk5GZIkITk5Wedx1SYnJwcvvfQSPD090bJlS9ja2sLHxwczZ85ETk4OAODdd9+FJEk4deqUwrmVlZWwtbWFJEm4dOmSwr6ysjK0atUK48ePVzuWX375BS+99BICAgJgZWVlkOtBRI0DEyIiajDXrl1D3759kZiYiIULF+Lw4cP4+OOPMXnyZPz000+4cuUKAGDo0KEAgKSkJIXzz549i3/++QdWVlZK+06dOoW7d+/Kz1XH6dOn8eWXX8LW1hbBwcFato6IjBkTIiJqMNu3b0dBQQESEhIwe/ZsDB06FGPHjsUbb7yBjIwMDBo0CADg6+uLhx56SKm3Jjk5Gc7Ozhg9erRSQlR9rCYJUVhYGG7cuIFvvvkG4eHhWrWNiIwbEyIiLURGRkKSJJw7dw5PP/00bGxsYGtri4ULF6K8vByXLl3CY489hjZt2qBz587417/+pVRHdnY2nn32WbRv3x4WFhbo3r071q9fj8rKSoXjbty4gYkTJ6JNmzawsbFBaGgo8vLyVMZ1+vRpjB49Gra2trC0tISvry/27t2rl2ugicLCQpiYmKB9+/Yq95uYmMj/d/Dgwfjhhx9QXl4u35+cnIxHHnkEQ4YMUZks2dvbo0ePHmrHU/19RET8rwGRDkycOBG9e/fGvn37MHPmTGzcuBELFizA2LFj8cQTT+DAgQN49NFH8frrr2P//v3y8/766y8EBgYiISEBq1atwqFDhzBs2DAsWrQIc+bMkR939+5dDBs2DAkJCYiKisK///1vODo6IjQ0VCmWpKQkBAUF4ebNm/jggw9w8OBB9OnTB6GhofUaawQA5eXlam1CiFrrCQgIQGVlJcaPH4/4+HjIZLIajx06dCiKi4vx008/AagaP5SSkoIhQ4ZgyJAhyM/Px6+//gqgavxQamoqHnnkEUiSVK82ElEzJ4io3pYvXy4AiPXr1yuU9+nTRwAQ+/fvl5fdu3dP2Nvbi/Hjx8vLlixZIgCIU6dOKZz/4osvCkmSxKVLl4QQQmzdulUAEAcPHlQ4bubMmQKA+OSTT+RlXl5ewtfXV9y7d0/h2CeffFI4OTmJiooKIYQQSUlJAoBISkqqtY1ZWVkCgFpbXXVVVlaKF154QZiYmAgAQpIk0b17d7FgwQKRlZWlcGxGRoYAINauXSuEECItLU0AEBcvXhRCCOHg4CDee+89IYQQR48eFQBETExMrd9fm3//+99qtYGImib2EBHpwJNPPqnwuXv37pAkCSNHjpSXmZmZoUuXLvjjjz/kZd9//z28vb0xYMAAhfOnTZsGIQS+//57AFW9Pm3atMHo0aMVjpsyZYrC58uXL+PixYt45plnACj27Dz++OPIzc1VejurLs7Ozvjpp5/U2vz8/GqtS5IkfPDBB7hy5QpiYmLw3HPP4d69e9i4cSN69OiBo0ePyo/t1asX7Ozs5I/GkpOT4ejoiG7dugEABg8eLB9HVJ/xQ0RE9zMzdABETYGtra3CZ3Nzc7Rq1QqWlpZK5fc/JiosLETnzp2V6nN2dpbvr/5fBwcHpeMcHR0VPv/5558AgEWLFmHRokUqYy0oKKijNYrMzc3Rp08ftY41NTVV6zhXV1e8+OKL8s979+7F5MmT8dprr+HHH38EUJU8DRkyBPHx8bh37x6SkpIwZMgQ+TlDhgxBZGQkhBBISkqCo6MjvLy81G8YEdF92ENEZEB2dnbIzc1VKr9x4wYAoF27dvLjqpOd+z04qLr6+IiIiBp7cdRNbqpdvXoVLVq0UGu7v4dHExMnTkSvXr3w888/K5QPHToUt2/fxqlTp3Ds2DGlhKigoABpaWk4efIke4eISCvsISIyoODgYERFReHMmTPo27evvDwuLg6SJMl/5IcOHYq9e/fi0KFDCo/NPv/8c4X6unXrhq5du+Ls2bNYu3atTmKsfmSmjurHWTXJzc2Fk5OTUnlxcTFycnLkPWPVqtu/ceNGFBUV4ZFHHpHv69GjB+zs7BAVFYWSkhImRESkFSZERAa0YMECxMXF4YknnsDKlSvh6uqKb775BjExMXjxxRfh6ekJAAgPD8fGjRsRHh6ONWvWoGvXrjh8+DDi4+OV6vzwww8xcuRIjBgxAtOmTUOHDh3w999/48KFCzhz5gz+/e9/axSjubk5+vXrp5P2rlmzBj/88ANCQ0PRp08ftGzZEllZWXjvvfdQWFiId955R+H4Hj16oH379jhw4ADs7e3RvXt3+T5JkjB48GAcOHAAQP3GD925cweHDx8GAJw8eRIAcPToURQUFMDKykphDBgRNW1MiIgMyN7eHidOnEBERAQiIiIgk8ng7u6Of/3rX1i4cKH8uFatWuH777/HvHnzsGTJEkiShJCQEOzevRuBgYEKdQ4dOhQ//vgj1qxZg/nz5+Off/6BnZ0dvL29MXHixIZuooKwsDAAwO7du/HOO++gqKgItra28PPzw+HDh1UmII888gj27t2r8Lis2pAhQ3DgwAF06NABXbp00Tie/Px8PP300wplkZGRAKrGOV29elXjOonIOElC1DFxCBEREVETx0HVRERE1OzxkRkRNTkVFRW1zpotSZLaUwQQUfPAHiIianKCg4NrnR7Aw8PD0CESUSPDMURE1ORcunQJt27dqnG/hYUFfHx8GjAiImrsmBARERFRs8dHZkRERNTscVC1CpWVlbhx4wbatGkDSZIMHQ4RETViQgjcunULzs7OMDHRXz9DSUkJysrKtK7H3NxcaZ1FYkKk0o0bN+Di4mLoMIiIyIjk5OSgY8eOeqm7pKQE9i1boxgVWtfl6OiIrKwsJkUPYEKkQps2bQAAvw92RxszvpprcWC7oUNoJHwNHQARNUIymQwuLi7y3w59KCsrQzEqsABusNBitEspKrExLwtlZWVMiB7AhEiF6sdkbcxMYc2ECBbWrQ0dQiNhbegAiKgRa4ghFi1hAkvU/3eJA4drxoSIiIjISJhAu6SGCVHNmBAREREZCSZE+sNrQ0RERM0ee4iIiIiMBHuI9IcJERERkZFgQqQ/vDZERETU7LGHiIiIyEhI/920OZ9UY0JERERkJPjITH94bYiIiKjZYw8RERGRkWAPkf4wISIiIjISTIj0h9eGiIiImj32EBERERkJCdr1ZPAts5oxISIiIjISfGSmP0yIiIiIjAQTIv3htSEiIqJmjz1ERERERoI9RPrDhIiIiMhIMCHSH14bIiIiavaYEBERERkJEx1s9RETEwM3NzdYWlrCz88Px44dU+u8H374AWZmZujTp49CeWxsLCRJUtpKSkrqGaH2DJoQRUZGKl0MR0dH+f5p06Yp7ff396+z3n379sHb2xsWFhbw9vbGgQMH9NkMIiKiBmGIhGjPnj2YP38+3nzzTaSnp2PQoEEYOXIksrOzaz2vqKgI4eHhCA4OVrnf2toaubm5CpulpWU9ItQNg/cQ9ejRQ+FinD9/XmH/Y489prD/8OHDtdaXmpqK0NBQhIWF4ezZswgLC8PEiRNx6tQpfTaDiIioSdqwYQOmT5+OGTNmoHv37oiOjoaLiwu2bt1a63kvvPACpkyZgoCAAJX7qztB7t8MyeAJkZmZmcLFsLe3V9hvYWGhsN/W1rbW+qKjozF8+HBERETAy8sLERERCA4ORnR0tB5bQUREpH+66iGSyWQKW2lpqcrvKysrQ1paGkJCQhTKQ0JCcOLEiRrj/OSTT/D7779j+fLlNR5TXFwMV1dXdOzYEU8++STS09PrbL8+GTwhyszMhLOzM9zc3DBp0iRcuXJFYX9ycjLat28PT09PzJw5E/n5+bXWl5qaqnTjRowYUeuNKy0tVfqHg4iIqLGRdLABgIuLC2xsbORbVFSUyu8rKChARUUFHBwcFModHByQl5en8pzMzEwsWbIEn332GczMVL/M7uXlhdjYWBw6dAi7du2CpaUlgoKCkJmZqfa10DWDvnY/cOBAxMXFwdPTE3/++SdWr16NwMBA/PLLL7Czs8PIkSPx9NNPw9XVFVlZWVi2bBkeffRRpKWlwcLCQmWdeXl5Gt04AIiKisKKFSt02jYiIqLGKicnB9bW1vLPNf2mVpMkxVXQhBBKZQBQUVGBKVOmYMWKFfD09KyxPn9/f4UxwUFBQejbty+2bNmCzZs3q9sMnTJoQjRy5Ej53z4+PggICICHhwd27tyJhQsXIjQ0VL6/Z8+e6NevH1xdXfHNN99g/PjxNdar7o2rFhERgYULF8o/y2QyuLi41KdJREREeqOrxV2tra0VEqKatGvXDqampkqdCvn5+UqdDwBw69YtnD59Gunp6ZgzZw4AoLKyEkIImJmZISEhAY8++qjSeSYmJujfv3/z7SF6kJWVFXx8fGq8IE5OTnB1da31gjk6Oqp946pZWFjUmR0TEREZWkNPzGhubg4/Pz8kJiZi3Lhx8vLExESMGTNG6Xhra2ull6NiYmLw/fff44svvoCbm5vK7xFCICMjAz4+PhpGqDuNKiEqLS3FhQsXMGjQIJX7CwsLkZOTAycnpxrrCAgIQGJiIhYsWCAvS0hIQGBgoM7jJSIiakiGmKl64cKFCAsLQ79+/RAQEIBt27YhOzsbs2fPBlD1lOX69euIi4uDiYkJevbsqXB++/btYWlpqVC+YsUK+Pv7o2vXrpDJZNi8eTMyMjLw/vvva9E67Rg0IVq0aBFGjRqFTp06IT8/H6tXr4ZMJsPUqVNRXFyMyMhIPPXUU3BycsLVq1fxxhtvoF27dgpZanh4ODp06CAfEDZv3jwMHjwY69atw5gxY3Dw4EEcOXIEx48fN1QziYiIjFZoaCgKCwuxcuVK5ObmomfPnjh8+DBcXV0BALm5uXXOSfSgmzdvYtasWcjLy4ONjQ18fX2RkpKCAQMG6KMJapGEEMJQXz5p0iSkpKSgoKAA9vb28Pf3x6pVq+Dt7Y27d+9i7NixSE9Px82bN+Hk5IShQ4di1apVCuN7HnnkEXTu3BmxsbHysi+++AJLly7FlStX4OHhgTVr1tQ65uhBMpkMNjY2yH+0K6zNTHXZZKNkEf9/hg6hkfAzdABE1AhV/2YUFRWpNS5Hm++IhQdaof6/S3dQgWn4Xa+xGiuDJkSNFRMiRUyIqjEhIiJlDZkQxekgIQpnQqSSwechIiIiIjK0RjWomoiIiGpmiEHVzQUTIiIiIiPBhEh/eG2IiIio2WMPERERkZFgD5H+MCEiIiIyEkyI9IfXhoiIiJo99hAREREZCfYQ6Q8TIiIiIiMi1X0I1QMTIiIiIiPBHiL94bUhIiKiZo89REREREaCPUT6w4SIiIjISDAh0h9eGyIiImr22ENERERkJCQAkhavmUlCZ6E0OUyIiIiIjISJJGCiRVZjAgEwKVKJj8yIiIio2WMPERERkZGQJC0fmQHsIaoBEyIiIiIjIUG7mao5y3XN+MiMiIiImj32EBERERmJqkdm9X/mxR6imjEhIiIiMhI6GUNEKjEhIiIiMhJMiPSHY4iIiIio2WMPERERkZHQycSMpBITIiIiIiPB1+71h4/MiIiIqNljDxEREZGx0HJQNdWMCVEtTidXwIr/4GFw+R1Dh9A48N8WIjIwvmWmP3xkRkRERM0e/z8vERGRkZAkoeVM1XzLrCZMiIiIiIyEiVS11ft83YXS5PDaEBERUbPHHiIiIiIjwUHV+mPQHqLIyEhIkqSwOTo6qjz2hRdegCRJiI6OrrXO2NhYpTolSUJJSYkeWkBERNRwJAitN1LN4D1EPXr0wJEjR+SfTU1NlY758ssvcerUKTg7O6tVp7W1NS5duqRQZmlpqV2gREREBsYeIv0xeEJkZmZWY68QAFy/fh1z5sxBfHw8nnjiCbXqrK2niYiIiOhBBh9UnZmZCWdnZ7i5uWHSpEm4cuWKfF9lZSXCwsLw2muvoUePHmrXWVxcDFdXV3Ts2BFPPvkk0tPT9RE6ERFRg6ruIdJmI9UMmhANHDgQcXFxiI+Px/bt25GXl4fAwEAUFhYCANatWwczMzPMnTtX7Tq9vLwQGxuLQ4cOYdeuXbC0tERQUBAyMzNrPKe0tBQymUxhIyIiamyqV7vXZiPVDPrIbOTIkfK/fXx8EBAQAA8PD+zcuRNDhgzBpk2bcObMGUgapLT+/v7w9/eXfw4KCkLfvn2xZcsWbN68WeU5UVFRWLFiRf0bQkREREbN4I/M7mdlZQUfHx9kZmbi2LFjyM/PR6dOnWBmZgYzMzP88ccfePXVV9G5c2e16zQxMUH//v1r7SGKiIhAUVGRfMvJydFBa4iIiHSLj8z0x+CDqu9XWlqKCxcuYNCgQQgLC8OwYcMU9o8YMQJhYWF47rnn1K5TCIGMjAz4+PjUeIyFhQUsLCzqHTcREVFDkKDdm2LMh2pm0IRo0aJFGDVqFDp16oT8/HysXr0aMpkMU6dOhZ2dHezs7BSOb9GiBRwdHdGtWzd5WXh4ODp06ICoqCgAwIoVK+Dv74+uXbtCJpNh8+bNyMjIwPvvv9+gbSMiIiLjYdCE6Nq1a5g8eTIKCgpgb28Pf39/nDx5Eq6urmrXkZ2dDROT/z35u3nzJmbNmoW8vDzY2NjA19cXKSkpGDBggD6aQERE1GC4uKv+SEIIXp0HyGQy2NjY4LCJO6wk5Ykim5vBJTsMHULjYDbI0BEQUSNU/ZtRVFQEa2trvX7HT21d0Fqq//DfYlGJ/v/kaBxrTEwM3nnnHeTm5qJHjx6Ijo7GoEF1/zfxhx9+wJAhQ9CzZ09kZGQo7Nu3bx+WLVuG33//HR4eHlizZg3GjRunaZN0plENqiYiIqLGZc+ePZg/fz7efPNNpKenY9CgQRg5ciSys7NrPa+oqAjh4eEIDg5W2peamorQ0FCEhYXh7NmzCAsLw8SJE3Hq1Cl9NaNO7CFSgT1EithD9F/sISIiFRqyhyitrQtam2jRQ1RZCT8Ne4gGDhyIvn37YuvWrfKy7t27Y+zYsfLxu6pMmjQJXbt2hampKb788kuFHqLQ0FDIZDJ8++238rLHHnsMbdu2xa5duzRvmA6wh4iIiMhIVI8h0mYDoDQZcWlpqcrvKysrQ1paGkJCQhTKQ0JCcOLEiRrj/OSTT/D7779j+fLlKvenpqYq1TlixIha69Q3JkRERERGRNJiq+bi4gIbGxv5VlNPT0FBASoqKuDg4KBQ7uDggLy8PJXnZGZmYsmSJfjss89gZqb63a28vDyN6mwIjWoeIiIiItK/nJwchUdmdc3F9+CKEUIIlatIVFRUYMqUKVixYgU8PT11UmdDYUJERERkJLSdbbr6XGtra7XGELVr1w6mpqZKPTf5+flKPTwAcOvWLZw+fRrp6emYM2cOgKqF2oUQMDMzQ0JCAh599FE4OjqqXWdD4SMzIiIiI6GrMUTqMjc3h5+fHxITExXKExMTERgYqHS8tbU1zp8/j4yMDPk2e/ZsdOvWDRkZGRg4cCAAICAgQKnOhIQElXU2FPYQERERUY0WLlyIsLAw9OvXDwEBAdi2bRuys7Mxe/ZsAFXrgV6/fh1xcXEwMTFBz549Fc5v3749LC0tFcrnzZuHwYMHY926dRgzZgwOHjyII0eO4Pjx4w3atvsxISIiIjISJlLVps35mgoNDUVhYSFWrlyJ3Nxc9OzZE4cPH5avKpGbm1vnnEQPCgwMxO7du7F06VIsW7YMHh4e2LNnj7wHyRA4D5EKnIdIEech+i/OQ0REKjTkPETn23dAGy3mIbpVWQmf/Ot6jdVYcQwRERERNXt8ZEZERGQkdPWWGSljQkRERGQkJAitVqznavc14yMzIiIiavbYQ0RERGQk+MhMf5gQERERGQnJRIKkxXv3EpgR1YQJERERkZGQTKq2ep+vu1CaHI4hIiIiomaPPURERERGgmOI9IcJERERkbHQdu0OPjSrER+ZERERUbPHHiIiIiIjwUHV+qNWQiSTyTSumIvGERER6ZYkSZC0GAikzblNnVoJ0UMPPaTRRZQkCb/99hvc3d3rHRgRERFRQ1H7kdkXX3wBW1vbOo8TQuDxxx/XKigiIiJSJklaPjLjUmY1UishcnV1xeDBg2FnZ6dWpe7u7mjRooVWgREREdED+N693qiVEGVlZWlU6c8//1yvYIiIiIgMgW+ZERERGQmt3zLjI7Ma1Ssh+u677/Ddd98hPz8flZWVCvs+/vhjnQRGREREirRe3FXwkVlNNE6IVqxYgZUrV6Jfv35wcnLiK3xEREQNhEOI9EfjhOiDDz5AbGwswsLC9BEPERERUYPTOCEqKytDYGCgPmIhIiKi2mg5hggcQ1QjjS/rjBkz8Pnnn+sjFiIiIqpN9eKu2mykklo9RAsXLpT/XVlZiW3btuHIkSPo1auX0nxDGzZs0G2ERERERHqmVkKUnp6u8LlPnz4AON8QERFRQ+Kgav1RKyFKSkrSy5dHRkZixYoVCmUODg7Iy8tTOvaFF17Atm3bsHHjRsyfP7/Wevft24dly5bh999/h4eHB9asWYNx48bpMnQiIqIGx9fu9UfjMUTPP/88bt26pVR++/ZtPP/88xoH0KNHD+Tm5sq38+fPKx3z5Zdf4tSpU3B2dq6zvtTUVISGhiIsLAxnz55FWFgYJk6ciFOnTmkcGxERETUPGidEO3fuxN27d5XK7969i7i4OI0DMDMzg6Ojo3yzt7dX2H/9+nXMmTMHn332mVrro0VHR2P48OGIiIiAl5cXIiIiEBwcjOjoaI1jIyIiakyqZ6rWZiPV1L40MpkMRUVFEELg1q1bkMlk8u2ff/7B4cOH0b59e40DyMzMhLOzM9zc3DBp0iRcuXJFvq+yshJhYWF47bXX0KNHD7XqS01NRUhIiELZiBEjcOLEiRrPKS0tVWiPTCbTuB1ERET6JkmS1huppvY8RA899JD8Ynp6eirtlyRJaTxQXQYOHIi4uDh4enrizz//xOrVqxEYGIhffvkFdnZ2WLduHczMzDB37ly168zLy4ODg4NCWU3jkqpFRUVpHDsRERE1HWonRElJSRBC4NFHH8W+fftga2sr32dubg5XV1e1xvjcb+TIkfK/fXx8EBAQAA8PD+zcuRNDhgzBpk2bcObMGY0z2gePF0LUWkdERITC1AIymQwuLi4afScREZHeSajHYJf7VNZ9SHOldkI0ZMgQlJeXIzw8HP369dNLwmBlZQUfHx9kZmbCxMQE+fn56NSpk3x/RUUFXn31VURHR+Pq1asq63B0dFTqDcrPz1fqNbqfhYUFLCwsdNIGIiIifeFr9/qjUZ5pZmaGffv2oaKiQi/BlJaW4sKFC3ByckJYWBjOnTuHjIwM+ebs7IzXXnsN8fHxNdYREBCAxMREhbKEhAQuN0JEREavamC0pMVm6BY0XhqvZRYcHIzk5GRMmzZN6y9ftGgRRo0ahU6dOiE/Px+rV6+GTCbD1KlTYWdnBzs7O4XjW7RoAUdHR3Tr1k1eFh4ejg4dOiAqKgoAMG/ePAwePBjr1q3DmDFjcPDgQRw5cgTHjx/XOl4iIiJqmjROiEaOHImIiAj8/PPP8PPzg5WVlcL+0aNHq13XtWvXMHnyZBQUFMDe3h7+/v44efIkXF1d1a4jOzsbJib/S3kDAwOxe/duLF26FMuWLYOHhwf27NmDgQMHql1ntZRKCZZg/+Kgy4cMHUKjIFKWGzqERsFk1veGDoGo2dL21Xn2ENVMEkJotPbt/cmHUmWSpLfHaQ1JJpPBxsYGS+ABS5gaOhyDe+uC+kluUyZS0gwdQqPAhIhIUfVvRlFREaytrfX6HTcCPWBtVv/fJVl5BZxP/K7XWI2Vxj1ElZUcok5ERERNi8YJERERERkGH5npT70uzdGjRzFq1Ch06dIFXbt2xejRo3Hs2DFdx0ZERET30e4NM+0Whm3qNE6I/u///g/Dhg1Dq1atMHfuXMyZMwctW7ZEcHAwPv/8c33ESERERKRXGj8yW7NmDf71r39hwYIF8rJ58+Zhw4YNWLVqFaZMmaLTAImIiKgKJ2bUH417iK5cuYJRo0YplY8ePRpZWVk6CYqIiIiU8ZGZ/micELm4uOC7775TKv/uu++4/hcREREZJY0fmb366quYO3cuMjIyEBgYCEmScPz4ccTGxmLTpk36iJGIiIiAqsVdtenkYQdRjTROiF588UU4Ojpi/fr12Lt3LwCge/fu2LNnD8aMGaPzAImIiKgKX7vXn3pdmnHjxuH48eMoLCxEYWEhjh8/zmSIiIhIzyRJyzFE9RxVHRMTAzc3N1haWsLPz6/WqXaOHz+OoKAg2NnZoWXLlvDy8sLGjRsVjomNja1qywNbSUlJveLThXpPzFhWVob8/Hylmas7deqkdVBERETUOOzZswfz589HTEwMgoKC8OGHH2LkyJH49ddfVf7mW1lZYc6cOejVqxesrKxw/PhxvPDCC7CyssKsWbPkx1lbW+PSpUsK51paWuq9PTXROCHKzMzE888/jxMnTiiUCyGazFpmREREjZEhXrvfsGEDpk+fjhkzZgAAoqOjER8fj61btyIqKkrpeF9fX/j6+so/d+7cGfv378exY8cUEiJJkuDo6Kh5QHqicUI0bdo0mJmZ4euvv4aTk1O9u9+IiIhIM9q+Ol99rkwmUyi3sLCAhYWF0vFlZWVIS0vDkiVLFMpDQkKUOkZqkp6ejhMnTmD16tUK5cXFxXB1dUVFRQX69OmDVatWKSRSDU3jhCgjIwNpaWnw8vLSRzxERESkZw9Ok7N8+XJERkYqHVdQUICKigo4ODgolDs4OCAvL6/W7+jYsSP++usvlJeXIzIyUt7DBABeXl6IjY2Fj48PZDIZNm3ahKCgIJw9exZdu3atf8O0oHFC5O3tjYKCAn3EQkRERLUxQT1fh7rvfAA5OTmwtraWF6vqHbrfg0+DqofJ1ObYsWMoLi7GyZMnsWTJEnTp0gWTJ08GAPj7+8Pf319+bFBQEPr27YstW7Zg8+bNmrRIZzROiNatW4fFixdj7dq18PHxQYsWLRT233+BiYiISIdMpKpNm/NR9Vutzu91u3btYGpqqtQblJ+fr9Rr9CA3NzcAgI+PD/78809ERkbKEyKlsExM0L9/f2RmZqrTCr3QOCEaNmwYACA4OFihnIOqiYiImhZzc3P4+fkhMTER48aNk5cnJiZqNN2OEAKlpaW17s/IyICPj49W8WpD44QoKSlJH3EQERFRXXT0yEwTCxcuRFhYGPr164eAgABs27YN2dnZmD17NgAgIiIC169fR1xcHADg/fffR6dOneRjjY8fP453330Xr7zyirzOFStWwN/fH127doVMJsPmzZuRkZGB999/X4vGaUfjhGjIkCH6iIOIiIjqoqNHZpoIDQ1FYWEhVq5cidzcXPTs2ROHDx+Gq6srACA3NxfZ2dny4ysrKxEREYGsrCyYmZnBw8MDb7/9Nl544QX5MTdv3sSsWbOQl5cHGxsb+Pr6IiUlBQMGDKh/27QkCSGEwb69kZLJZLCxscESeMASpoYOx+DeujDa0CE0CiIlzdAhNAoms743dAhEjUr1b0ZRUZHextFWf8ff47vDukX9f5dk9ypgu/+CXmM1VlzVhIiIiJq9ei/dQURERA3MAI/MmgsmRERERMaCCZHeaPzI7O7du7hz54788x9//IHo6GgkJCToNDAiIiKihqJxQjRmzBj5q3U3b97EwIEDsX79eowZMwZbt27VeYBERET0XxL+9+p9fTZ2ENVI44TozJkzGDRoEADgiy++gIODA/744w/ExcUZbLptIiKiZqH6kZk2G6mkcUJ0584dtGnTBgCQkJCA8ePHw8TEBP7+/vjjjz90HiARERGRvmmcEHXp0gVffvklcnJyEB8fj5CQEABV65pwTgMiIiI90uZxmbazXDdxGl+at956C4sWLULnzp0xYMAABAQEAKjqLfL19dV5gERERPRffGSmNxq/dj9hwgQ8/PDDyM3NRe/eveXlwcHBCgu/ERERERmLenWeOTo6ok2bNkhMTMTdu3cBAP3795cv5EZERER6IOlgI5U0TogKCwsRHBwMT09PPP7448jNzQUAzJgxA6+++qrOAyQiIqL/4iMzvdE4IVqwYAFatGiB7OxstGrVSl4eGhqK//znPzoNjoiIiO7DhEhvNB5DlJCQgPj4eHTs2FGhvGvXrnztnoiIiIySxgnR7du3FXqGqhUUFMDCwkInQREREZEK2r46z9fua6TxpRk8eLB86Q4AkCQJlZWVeOeddzB06FCN6oqMjIQkSQqbo6Ojwn4vLy9YWVmhbdu2GDZsGE6dOlVrnbGxsUp1SpKEkpISzRpKRETU2PCRmd5o3EP0zjvv4JFHHsHp06dRVlaGxYsX45dffsHff/+NH374QeMAevTogSNHjsg/m5qayv/29PTEe++9B3d3d9y9excbN25ESEgILl++DHt7+xrrtLa2xqVLlxTKLC0tNY6NiIiImgeNEyJvb2+cO3cOMTExMDU1xe3btzF+/Hi8/PLLcHJy0jwAMzOFXqH7TZkyReHzhg0bsGPHDpw7dw7BwcE11vlgTxMREVFTIEmApMVjL4kdRDXSOCECquYhWrlypU4CyMzMhLOzMywsLDBw4ECsXbsW7u7uSseVlZVh27ZtsLGxUZgQUpXi4mK4urqioqICffr0wapVq2qdRbu0tBSlpaXyzzKZrP4NIiIi0hdtH3vxkVmN6pVnHjt2DM8++ywCAwNx/fp1AMCnn36K48ePa1TPwIEDERcXh/j4eGzfvh15eXkIDAxEYWGh/Jivv/4arVu3hqWlJTZu3IjExES0a9euxjq9vLwQGxuLQ4cOYdeuXbC0tERQUBAyMzNrPCcqKgo2NjbyzcXFRaN2EBERkXHTOCHat28fRowYgZYtW+LMmTPynpVbt25h7dq1GtU1cuRIPPXUU/Dx8cGwYcPwzTffAAB27twpP2bo0KHIyMjAiRMn8Nhjj2HixInIz8+vsU5/f388++yz6N27NwYNGoS9e/fC09MTW7ZsqfGciIgIFBUVybecnByN2kFERNQguLir3mh8aVavXo0PPvgA27dvR4sWLeTlgYGBOHPmjFbBWFlZwcfHR6E3x8rKCl26dIG/vz927NgBMzMz7NixQ+06TUxM0L9//1p7iCwsLGBtba2wERERNTp8y0xvNE6ILl26hMGDByuVW1tb4+bNm1oFU1paigsXLtQ6OFsIoTDepy5CCGRkZNRrwDcRERE1DxonRE5OTrh8+bJS+fHjx1UOhq7NokWLcPToUWRlZeHUqVOYMGECZDIZpk6ditu3b+ONN97AyZMn8ccff+DMmTOYMWMGrl27hqefflpeR3h4OCIiIuSfV6xYgfj4eFy5cgUZGRmYPn06MjIyMHv2bE2bSkRE1Liwh0hvNH7L7IUXXsC8efPw8ccfQ5Ik3LhxA6mpqVi0aBHeeustjeq6du0aJk+ejIKCAtjb28Pf3x8nT56Eq6srSkpKcPHiRezcuRMFBQWws7ND//79cezYMfTo0UNeR3Z2NkxM/pfX3bx5E7NmzUJeXh5sbGzg6+uLlJQUDBgwQNOmEhERNS6cqVpvNE6IFi9ejKKiIgwdOhQlJSUYPHgwLCwssGjRIsyZM0ejunbv3l3jPktLS+zfv7/OOpKTkxU+b9y4ERs3btQoDiIiIqNgAi1fu9dZJE2ORpemoqICR48exauvvoqCggL8+OOPOHnyJP766y+sWrVKXzESERERKfn0008RFBQEZ2dn+QLz0dHROHjwoMZ1aZQQmZqaYsSIESgqKkKrVq3Qr18/DBgwAK1bt9b4i4mIiEhDfO1ebuvWrVi4cCEef/xx3Lx5ExUVFQCAhx56CNHR0RrXp/Gl8fHxwZUrVzT+IiIiItISB1XLbdmyBdu3b8ebb76psA5qv379cP78eY3r0zghWrNmDRYtWoSvv/4aubm5kMlkChsRERGRvmVlZalclsvCwgK3b9/WuD6NB1U/9thjAIDRo0dDum+VOCEEJEmSd1kRERGRjknQ7rFX0+kggpubGzIyMuDq6qpQ/u2338Lb21vj+jROiJKSkjT+EiIiItIBLu4q99prr+Hll19GSUkJhBD48ccfsWvXLkRFReGjjz7SuD6NE6IhQ4Zo/CVEREREuvTcc8+hvLwcixcvxp07dzBlyhR06NABmzZtwqRJkzSuT+OE6Ny5cyrLJUmCpaUlOnXqBAsLC40DISIiojpwYkYFM2fOxMyZM1FQUIDKykq0b9++3nVpnBD16dNHYezQg1q0aIHQ0FB8+OGHsLS0rHdgRERE9AA+MpO7e/cuhBBo1aoV2rVrhz/++APR0dHw9vZGSEiIxvVpnCseOHAAXbt2xbZt25CRkYH09HRs27YN3bp1w+eff44dO3bg+++/x9KlSzUOhoiIiEgdY8aMQVxcHICqZbsGDBiA9evXY8yYMdi6davG9WncQ7RmzRps2rQJI0aMkJf16tULHTt2xLJly/Djjz/CysoKr776Kt59912NAyIiIqIasIdI7syZM/Klur744gs4OjoiPT0d+/btw1tvvYUXX3xRo/o0TojOnz+v9IobALi6usonQurTpw9yc3M1rZqIiIhqwzFEcnfu3EGbNm0AAAkJCRg/fjxMTEzg7+8vX8ZDExpfGi8vL7z99tsoKyuTl927dw9vv/02vLy8AADXr1+Hg4ODxsEQERFRLThTtVyXLl3w5ZdfIicnB/Hx8fJxQ/n5+bC2tta4Po17iN5//32MHj0aHTt2RK9evSBJEs6dO4eKigp8/fXXAIArV67gpZde0jgYIiIiInW89dZbmDJlChYsWIDg4GAEBAQAqOotUjWDdV0kIYTQ9KTi4mL83//9H3777TcIIeDl5YUpU6bIu66MnUwmg42NDc6374A2Jk2of7GeXHOzDR1Co3DBtrOhQ2gUuv991dAhUCNy7wnNZwRuamT3KtAu8TcUFRXVq2dCre/47+/SP2v9YW2pcV/G/+opKUfbN07qNdaGlJeXh9zcXPTu3Rsm//29/vHHH2FtbS1/aqWuel3V1q1bY/bs2fU5lYiIiOpLkqo2bc5vQhwdHeHo6KhQNmDAgHrVVa+E6NNPP8WHH36IK1euIDU1Fa6urti4cSPc3d0xZsyYegVCREREpK6SkhJs2bIFSUlJyM/PR2VlpcL+M2fOaFSfxgnR1q1b8dZbb2H+/PlYvXq1fDHXtm3bIjo6mgkRERGRvkjQboHWJtRB9PzzzyMxMRETJkzAgAEDap00Wh0aD5DZsmULtm/fjjfffBNmZv/Lp/r16yd/7Z6IiIj0oPqRmTZbPcTExMDNzQ2Wlpbw8/PDsWPHajz2+PHjCAoKgp2dHVq2bAkvLy/5fEH327dvH7y9vWFhYQFvb28cOHBAo5i++eYbfPnll9i6dSsiIyOxfPlyhU1TGidEWVlZKkdvW1hY4Pbt2xoHQERERI3Xnj17MH/+fLz55ptIT0/HoEGDMHLkSGRnq37hxsrKCnPmzEFKSgouXLiApUuXYunSpdi2bZv8mNTUVISGhiIsLAxnz55FWFgYJk6ciFOnTqkdV4cOHXT6MpfGCZGbmxsyMjKUyr/99lt4e/NtAyIiIr2StNjqYcOGDZg+fTpmzJiB7t27Izo6Gi4uLjUuj+Hr64vJkyejR48e6Ny5M5599lmMGDFCoVcpOjoaw4cPR0REBLy8vBAREYHg4GBER0erHdf69evx+uuv12sSRlU0Tohee+01vPzyy9izZw+EEPjxxx+xZs0avPHGG3jttdd0EhQRERGp0MCPzMrKypCWlqa0WGpISAhOnDihVh3p6ek4ceIEhgwZIi9LTU1VqnPEiBFq1wlUDdUpKSmBu7s72rRpA1tbW4VNUxoPqn7uuedQXl6OxYsX486dO5gyZQo6dOiATZs2YdKkSRoHQERERA1LJpMpfLawsICFhYXScQUFBaioqFBafcLBwQF5eXm1fkfHjh3x119/oby8HJGRkZgxY4Z8X15eXr3qvN/kyZNx/fp1rF27Fg4ODloPqq7Xa/czZ87EzJkzUVBQgMrKSrRv316rIIiIiEgNOlrLzMXFRaF4+fLliIyMrPG0B5MNIUSdCcixY8dQXFyMkydPYsmSJejSpQsmT56sVZ33O3HiBFJTU9G7d2+1z6lN/ae7BNCuXTudBEFERERq0NHEjDk5OQozVavqHQKqfudNTU2Vem7y8/PrXLPUzc0NAODj44M///wTkZGR8oTI0dGxXnXez8vLC3fv3lX7+LqolRD5+vqqnbVpOhESERERqUlH8xBZW1urtXSHubk5/Pz8kJiYiHHjxsnLExMTNZp3UAiB0tJS+eeAgAAkJiZiwYIF8rKEhAQEBgaqXefbb7+NV199FWvWrIGPjw9atGihsF/TpUnUSojGjh0r/7ukpAQxMTHw9vaWL6R28uRJ/PLLL1zQlYiIqIlZuHAhwsLC0K9fPwQEBGDbtm3Izs6WL+EVERGB69evIy4uDkDVIvCdOnWSryV2/PhxvPvuu3jllVfkdc6bNw+DBw/GunXrMGbMGBw8eBBHjhzB8ePH1Y7rscceAwAEBwcrlFc/equeOFpdaiVE909wNGPGDMydOxerVq1SOiYnJ0ejLyciIiINGGAts9DQUBQWFmLlypXIzc1Fz549cfjwYbi6ugIAcnNzFeYkqqysREREBLKysmBmZgYPDw+8/fbbeOGFF+THBAYGYvfu3Vi6dCmWLVsGDw8P7NmzBwMHDlQrpnv37mHw4MGYPHmyxou41kTj1e5tbGxw+vRpdO3aVaE8MzMT/fr1Q1FRkU4CMySudq+Iq91X4Wr3VbjaPd2Pq9038Gr30YNg3VKL1e7vlqPt/GNNYrV7e3t7nDhxQikfqS+Nf+1btmypskvr+PHjsLS01ElQRERERLUJDw/Hjh07dFafxmnm/Pnz8eKLLyItLQ3+/v4AqsYQffzxx3jrrbd0FhgRERE9wACPzBqrsrIyfPTRR0hMTES/fv1gZWWlsH/Dhg0a1adxQrRkyRK4u7tj06ZN+PzzzwEA3bt3R2xsLCZOnKhpdURERKQuHc1D1BT8/PPP6Nu3LwDgt99+U9hXn0ka6/UgcuLEiUx+iIiIyGCSkpJ0Wp9WEzMSERFRA+IjM71Rq/PM1tYWBQUFalfaqVMntVafjYyMhCRJCpujo6PCfi8vL1hZWaFt27YYNmwYTp06VWe9+/btg7e3NywsLODt7Y0DBw6oHTsREVGjpc1K99pO6tjEqdVDdPPmTXz77bewsbFRq9LCwkK1J0Tq0aMHjhw5Iv9samoq/9vT0xPvvfce3N3dcffuXWzcuBEhISG4fPky7O3tVdaXmpqK0NBQrFq1CuPGjcOBAwcwceJEHD9+XO35DYiIiKh5UfuR2dSpU/UTgJmZQq/Q/aZMmaLwecOGDdixYwfOnTunNDNltejoaAwfPhwREREAqmbQPHr0KKKjo7Fr1y7dBk9ERNSQ+MhMb9R6ZFZZWanx5u7urlYAmZmZcHZ2hpubGyZNmoQrV66oPK6srAzbtm2DjY1NrSvbpqamIiQkRKFsxIgROHHihFrxEBERNVbV+ZA2G6lm0EHVAwcORFxcHDw9PfHnn39i9erVCAwMxC+//AI7OzsAwNdff41Jkybhzp07cHJyQmJiItq1a1djnXl5eUqr5To4OCitqnu/0tJShUXnZDKZli0jIiLSA/YQ6Y1BZyQYOXIknnrqKfj4+GDYsGH45ptvAAA7d+6UHzN06FBkZGTgxIkTeOyxxzBx4kTk5+fXWu+D8w9UL/RWk6ioKNjY2Mg3FxcXLVpFRERExqZRTdFkZWUFHx8fZGZmKpR16dIF/v7+2LFjB8zMzGqdqtvR0VGpNyg/P1+p1+h+ERERKCoqkm9cpJaIiBolvmWmN40qISotLcWFCxfg5ORU4zFCCIXHWw8KCAhAYmKiQllCQgICAwNrPMfCwgLW1tYKGxERUaMjSYCJFhsfmdXIoGOIFi1ahFGjRqFTp07Iz8/H6tWrIZPJMHXqVNy+fRtr1qzB6NGj4eTkhMLCQsTExODatWt4+umn5XWEh4ejQ4cOiIqKAgDMmzcPgwcPxrp16zBmzBgcPHgQR44cUbkgLRERERFQzx6i33//HUuXLsXkyZPl43n+85//4JdfftGonmvXrmHy5Mno1q0bxo8fD3Nzc5w8eRKurq4wNTXFxYsX8dRTT8HT0xNPPvkk/vrrLxw7dgw9evSQ15GdnY3c3Fz558DAQOzevRuffPIJevXqhdjYWOzZs4dzEBERkfHjIzO90biH6OjRoxg5ciSCgoKQkpKCNWvWoH379jh37hw++ugjfPHFF2rXtXv37hr3WVpaYv/+/XXWkZycrFQ2YcIETJgwQe04iIiIjALfMtMbjXuIlixZgtWrVyMxMRHm5uby8qFDhyI1NVWnwRERERE1BI0TovPnz2PcuHFK5fb29igsLNRJUERERKQCH5npjcYJ0UMPPaQwZqdaeno6OnTooJOgiIiISAVOVa03GidEU6ZMweuvv468vDxIkoTKykr88MMPWLRoEcLDw/URIxEREZFeaZwQrVmzBp06dUKHDh1QXFwMb29vDB48GIGBgVi6dKk+YiQiIiKAj8z0SOO3zFq0aIHPPvsMK1euRHp6OiorK+Hr64uuXbvqIz4iIiKqZoKqCRa1OZ9UqvfEjB4eHvDw8NBlLERERFQbbXt52ENUI7USooULF6pd4YYNG+odDBEREZEhqJUQpaenK3xOS0tDRUUFunXrBgD47bffYGpqCj8/P91HSERERFU4MaPeqJUQJSUlyf/esGED2rRpg507d6Jt27YAgH/++QfPPfccBg0apJ8oiYiIiI/M9Ejj4VXr169HVFSUPBkCgLZt22L16tVYv369ToMjIiIiaggaJ0QymQx//vmnUnl+fj5u3bqlk6CIiIhIBU7MqDcaJ0Tjxo3Dc889hy+++ALXrl3DtWvX8MUXX2D69OkYP368PmIkIiIigPMQ6ZHGr91/8MEHWLRoEZ599lncu3evqhIzM0yfPh3vvPOOzgMkIiIi0jeNE6JWrVohJiYG77zzDn7//XcIIdClSxdYWVnpIz4iIiKqZiJpOTEju4hqUu+JGa2srNCrVy9dxkJERES14Wv3eqNxQjR06FBItVzQ77//XquAiIiIiBqaxglRnz59FD7fu3cPGRkZ+PnnnzF16lRdxUVEREQPYg+R3micEG3cuFFleWRkJIqLi7UOiIiIiGrAhEhvdLbu7bPPPouPP/5YV9URERHRgyQT7TdSSWdXJjU1FZaWlrqqjoiIiKjBaPzI7MHJF4UQyM3NxenTp7Fs2TKdBUZEREQPkLR87Z6PzGqkcUJkbW2t8JaZiYkJunXrhpUrVyIkJESnwREREdF9OIZIbzROiGJjY/UQBhEREZHhaDyGyN3dHYWFhUrlN2/ehLu7u06CIiIiIhU4qFpvNO4hunr1KioqKpTKS0tLcf36dZ0ERURERCrwkZneqJ0QHTp0SP53fHw8bGxs5J8rKirw3XffoXPnzjoNjoiIiKghqJ0QjR07FgAgSZLSjNQtWrRA586dsX79ep0GR0RERPfh4q56o3ZCVFlZCQBwc3PDTz/9hHbt2uktKCIiIlJB23FAHENUI42vTFZWFpMhIiKiZiQmJgZubm6wtLSEn58fjh07VuOx+/fvx/Dhw2Fvbw9ra2sEBAQgPj5e4ZjY2FhIkqS0lZSU6LspNVKrh2jz5s2YNWsWLC0tsXnz5lqPnTt3rk4Caww6Zf4Ka2trQ4dBjUT8PxaGDqFR6G7oABoJkRlh6BAahRbf/GroEAyuhUwG3DeuVq8MMKh6z549mD9/PmJiYhAUFIQPP/wQI0eOxK+//opOnTopHZ+SkoLhw4dj7dq1eOihh/DJJ59g1KhROHXqFHx9feXHWVtb49KlSwrnGnLFC7USoo0bN+KZZ56BpaVljYu7AlXji5pSQkRERNSoGCAh2rBhA6ZPn44ZM2YAAKKjoxEfH4+tW7ciKipK6fjo6GiFz2vXrsXBgwfx1VdfKSREkiTB0dFR43j0Ra2EKCsrS+XfRERE1IAkScsxRJolRGVlZUhLS8OSJUsUykNCQnDixAm16qisrMStW7dga2urUF5cXAxXV1dUVFSgT58+WLVqlULC1NA0vqorV67EnTt3lMrv3r2LlStX6iQoIiIi0h+ZTKawlZaWqjyuoKAAFRUVcHBwUCh3cHBAXl6eWt+1fv163L59GxMnTpSXeXl5ITY2FocOHcKuXbtgaWmJoKAgZGZm1r9RWtI4IVqxYgWKi4uVyu/cuYMVK1boJCgiIiJSofq1e202AC4uLrCxsZFvqh593U96oGdJCKFUpsquXbsQGRmJPXv2oH379vJyf39/PPvss+jduzcGDRqEvXv3wtPTE1u2bKnHRdENjROimi7C2bNnlbrD6hIZGak0wrz6eeK9e/fw+uuvw8fHB1ZWVnB2dkZ4eDhu3LhRa52NceQ6ERGRTlSPIdJmA5CTk4OioiL5FhGh+iWBdu3awdTUVKk3KD8/X6nX6EF79uzB9OnTsXfvXgwbNqzWY01MTNC/f3+D9hCpPQ9R27Zt5cmFp6enQlJUUVGB4uJizJ49W+MAevTogSNHjsg/m5qaAqjqcTpz5gyWLVuG3r17459//sH8+fMxevRonD59utY6G9vIdSIiosbE2tparbeozc3N4efnh8TERIwbN05enpiYiDFjxtR43q5du/D8889j165deOKJJ+r8HiEEMjIy4OPjo14D9EDthCg6OhpCCDz//PNYsWKFwtId5ubm6Ny5MwICAjQPwMxM5ShzGxsbJCYmKpRt2bIFAwYMQHZ2tspX/ao1tpHrREREOmGAiRkXLlyIsLAw9OvXDwEBAdi2bRuys7PlnSARERG4fv064uLiAFQlQ+Hh4di0aRP8/f3lvUstW7aU5w4rVqyAv78/unbtCplMhs2bNyMjIwPvv/9+/dumJbUTourlOtzc3BAYGIgWLVroJIDMzEw4OzvDwsICAwcOxNq1a+Hu7q7y2KKiIkiShIceeqjWOjUduV5aWqowoEwmk9WrLURERHplgNfuQ0NDUVhYiJUrVyI3Nxc9e/bE4cOH4erqCgDIzc1Fdna2/PgPP/wQ5eXlePnll/Hyyy/Ly6dOnYrY2FgAwM2bNzFr1izk5eXBxsYGvr6+SElJwYABA+rfNi1JQghR35Pv3r2Le/fuKZRpMpHht99+izt37sDT0xN//vknVq9ejYsXL+KXX36BnZ2dwrElJSV4+OGH4eXlhf/7v/+rsc6TJ0/i8uXL8PHxgUwmw6ZNm3D48GGcPXsWXbt2VXlOZGSkygHhRUVFnJiR5KKlboYOoVGYLy7VfVAzwIkZq0hdax+M2xzIZDLY2Njo9Tej+jtuJj4Layvz+tdzuwwPDf8//r6poHFCdOfOHSxevBh79+5FYWGh0v6Kiop6B3P79m14eHhg8eLFWLhwobz83r17ePrpp5GdnY3k5GSNbmJlZSX69u2LwYMH1zjLtqoeIhcXF/4DQwqYEFVhQlSFCVEVJkQNnBB9F6Z9QhT8KX/fVND4YeJrr72G77//HjExMbCwsMBHH32EFStWwNnZWf78sL6srKzg4+OjMMr83r17mDhxIrKyspCYmKjxDVRn5LqFhYV8gJm6A82IiIgansn/xhHVZ9P8Z7/Z0PjKfPXVV4iJicGECRNgZmaGQYMGYenSpVi7di0+++wzrYIpLS3FhQsX4OTkBOB/yVBmZiaOHDmi9BhNHdUj16vrJCIiInqQxgnR33//DTc3NwBV44X+/vtvAMDDDz+MlJQUjepatGgRjh49iqysLJw6dQoTJkyATCbD1KlTUV5ejgkTJuD06dP47LPPUFFRgby8POTl5aGsrExeR3h4uML8CStWrEB8fDyuXLmCjIwMTJ8+HRkZGfWaEoCIiKhR0dE8RKRM7bfMqrm7u+Pq1atwdXWFt7c39u7diwEDBuCrr76q8+2vB127dg2TJ09GQUEB7O3t4e/vj5MnT8LV1RVXr17FoUOHAAB9+vRROC8pKQmPPPIIACA7OxsmJv/L6xrjyHUiIiKdMMBbZs2FxgnRc889h7Nnz2LIkCGIiIjAE088gS1btqC8vBwbNmzQqK7du3fXuK9z585QZ7x3cnKywueNGzdi48aNGsVBRERkFJgQ6Y3GCdGCBQvkfw8dOhQXL17E6dOn4eHhgd69e+s0OCIiIqKGoPVw806dOmH8+PGwtbXF888/r4uYiIiISBUTE+03UklnV+bvv//Gzp07dVUdERERPYiDqvWGqSIRERE1exqPISIiIiID4aBqvWFCREREZCwkScvV7pkQ1UTthGj8+PG17r9586a2sRAREREZhNoJkY2NTZ37w8PDtQ6IiIiIamAiVW3anE8qqZ0QffLJJ/qMg4iIiOrCMUR6w7fMiIiIqNnjoGoiIiJjIZloOaia/SA1YUJERERkLPjITG+YEBERERkLJkR6w74zIiIiavbYQ0RERGQstF2glYu71ogJERERkdGQ/rtpcz6pwlSRiIiImj32EBERERkLDqrWGyZERERERkPLeYj4YKhGvDJERETU7LGHiIiIyGhwULW+MCEiIiIyFhxDpDd8ZEZERETNHnuIiIiIjAUXd9UbJkRERERGg2OI9IUJERERkbGQoOUYIp1F0uSw74yIiIiaPfYQERERGQ0TaNeXwX6QmjAhIiIiMhZ87V5vmCoSERFRs8ceIiIiImPBHiK9YUJERERkNPjavb7wkRkRERE1e+whIiIiMhacqVpvmBAREREZC44h0huDpoqRkZGQJElhc3R0BADcu3cPr7/+Onx8fGBlZQVnZ2eEh4fjxo0bdda7b98+eHt7w8LCAt7e3jhw4IC+m0JERERGzOB9Zz169EBubq58O3/+PADgzp07OHPmDJYtW4YzZ85g//79+O233zB69Oha60tNTUVoaCjCwsJw9uxZhIWFYeLEiTh16lRDNIeIiEiPJB1spIrBH5mZmZnJe4XuZ2Njg8TERIWyLVu2YMCAAcjOzkanTp1U1hcdHY3hw4cjIiICABAREYGjR48iOjoau3bt0n0DiIiIGgrHEOmNwa9MZmYmnJ2d4ebmhkmTJuHKlSs1HltUVARJkvDQQw/VeExqaipCQkIUykaMGIETJ07UeE5paSlkMpnCRkRE1Ng8OMykPhupZtCEaODAgYiLi0N8fDy2b9+OvLw8BAYGorCwUOnYkpISLFmyBFOmTIG1tXWNdebl5cHBwUGhzMHBAXl5eTWeExUVBRsbG/nm4uJS/0YRERE1MTExMXBzc4OlpSX8/Pxw7NixGo/dv38/hg8fDnt7e1hbWyMgIADx8fFKxzW28b4GTYhGjhyJp556Cj4+Phg2bBi++eYbAMDOnTsVjrt37x4mTZqEyspKxMTE1FnvgxmwEKLWrDgiIgJFRUXyLScnpx6tISIi0reGH0O0Z88ezJ8/H2+++SbS09MxaNAgjBw5EtnZ2SqPT0lJwfDhw3H48GGkpaVh6NChGDVqFNLT0+XHNMbxvgYfQ3Q/Kysr+Pj4IDMzU1527949TJw4EVlZWfj+++9r7R0CAEdHR6XeoPz8fKVeo/tZWFjAwsJCu+CJiIj0zQBjiDZs2IDp06djxowZAKrG6sbHx2Pr1q2IiopSOj46Olrh89q1a3Hw4EF89dVX8PX1lR/T2Mb7GnwM0f1KS0tx4cIFODk5AfhfMpSZmYkjR47Azs6uzjoCAgKUBmMnJCQgMDBQLzETEREZmwfHzZaWlqo8rqysDGlpaUpjc0NCQmodm3u/yspK3Lp1C7a2tvKy+oz31TeDJkSLFi3C0aNHkZWVhVOnTmHChAmQyWSYOnUqysvLMWHCBJw+fRqfffYZKioqkJeXh7y8PJSVlcnrCA8Pl2eYADBv3jwkJCRg3bp1uHjxItatW4cjR45g/vz5BmghERGRLunmkZmLi4vC2FlVPT0AUFBQgIqKCo3H5t5v/fr1uH37NiZOnCgvq894X30z6COza9euYfLkySgoKIC9vT38/f1x8uRJuLq64urVqzh06BAAoE+fPgrnJSUl4ZFHHgEAZGdnw8Tkf3ldYGAgdu/ejaVLl2LZsmXw8PDAnj17MHDgwIZqFhERkX7oaKbqnJwchSEodQ0b0XRsbrVdu3YhMjISBw8eRPv27XVSp74YNCHavXt3jfs6d+4MIUSddSQnJyuVTZgwARMmTNAmNCIioibL2tq6zjG5ANCuXTuYmppqPDYXqBqMPX36dPz73//GsGHDFPbVZ7yvvjWqQdXUOImUcEOH0CjMfd/J0CFQIyJ1Vf2IgUivJEnLQdWa9cCYm5vDz88PiYmJGDdunLw8MTERY8aMqfG8Xbt24fnnn8euXbvwxBNPKO2vHu+7YMECeZmhx/syISIiIjIa2i6/ofm5CxcuRFhYGPr164eAgABs27YN2dnZmD17NoCqN8SuX7+OuLg4AFXJUHh4ODZt2gR/f395T1DLli1hY2MDoGq87+DBg7Fu3TqMGTMGBw8exJEjR3D8+HEt2qadRvWWGRERETUuoaGhiI6OxsqVK9GnTx+kpKTg8OHDcHV1BQDk5uYqzEn04Ycfory8HC+//DKcnJzk27x58+THVI/3/eSTT9CrVy/ExsYafLyvJNQZqNPMyGQy2NjYoKioSK1nrE0dH5lVET+rnoSsuTF5KdnQIRA1Kg3xmyH/jsx3Yd2mZf3ruXUXNl0X8fdNBT4yIyIiMhZc3FVvmBAREREZjYYfQ9RcMFUkIiKiZo89RERERMZCRxMzkjImREREREbDBNo93OGDoZrwyhAREVGzxx4iIiIiY8FHZnrDhIiIiMhYMCHSGz4yIyIiomaPPURERERGg4Oq9YUJERERkdHQ8pEZJ2asEVNFIiIiavbYQ0RERGQ0uHSHvjAhIiIiMhZc3FVvmBAREREZCwlavnavs0iaHKaKRERE1Oyxh4iIiMhocAyRvjAhIiIiMhYcQ6Q3vDJERETU7LGHiIiIyGjwkZm+MCEiIiIyFlzcVW/4yIyIiIiaPfYQERERGQ0u7qovTIiIiIiMBR+Z6Q1TRSIiImr22ENERERkLDgPkd4wISIiIjIafO1eX5gQERERGQuOIdIb9p0RERFRs8ceIiIiIqPB1+71hQkRERGRseAjM70xaKoYGRkJSZIUNkdHR/n+/fv3Y8SIEWjXrh0kSUJGRkaddcbGxirVKUkSSkpK9NgSIiIiMmYG7yHq0aMHjhw5Iv9samoq//v27dsICgrC008/jZkzZ6pdp7W1NS5duqRQZmlpqX2wREREBsVHZvpi8ITIzMxMoVfofmFhYQCAq1evalTngz1NRERETQIfmemNwVPFzMxMODs7w83NDZMmTcKVK1e0rrO4uBiurq7o2LEjnnzySaSnp9d6fGlpKWQymcJGREREzYdBE6KBAwciLi4O8fHx2L59O/Ly8hAYGIjCwsJ61+nl5YXY2FgcOnQIu3btgqWlJYKCgpCZmVnjOVFRUbCxsZFvLi4u9f5+IiIivZHwv16iem2GbkDjJQkhhKGDqHb79m14eHhg8eLFWLhwobz86tWrcHNzQ3p6Ovr06aNRnZWVlejbty8GDx6MzZs3qzymtLQUpaWl8s8ymQwuLi4oKiqCtbV1vdrSlIiUcEOH0CiIn7MNHUKjYPJSsqFDIGpUZDIZbGxs9PqbIf+OP/fD2tpKi3puw8ZhPH/fVDD4GKL7WVlZwcfHp9beHE2ZmJigf//+tdZpYWEBCwsLnX0nERERGReDjyG6X2lpKS5cuAAnJyed1SmEQEZGhk7rJCIiMgitHpdpOSC7iTNoD9GiRYswatQodOrUCfn5+Vi9ejVkMhmmTp0KAPj777+RnZ2NGzduAID8VXpHR0f5W2Th4eHo0KEDoqKiAAArVqyAv78/unbtCplMhs2bNyMjIwPvv/++AVpIRESkS1zcVV8MmhBdu3YNkydPRkFBAezt7eHv74+TJ0/C1dUVAHDo0CE899xz8uMnTZoEAFi+fDkiIyMBANnZ2TAx+V9H182bNzFr1izk5eXBxsYGvr6+SElJwYABAxquYURERPogmVRt2pxPKjWqQdWNRUMMkDMmHFRdhYOqq3BQNZGiBh1U/dch7QdV24/m75sKjWpQNREREdWGj8z0hQkRERGRseAjM73hlSEiIqJmjwkRERGR0ZB0sGkuJiYGbm5usLS0hJ+fH44dO1bjsbm5uZgyZQq6desGExMTzJ8/X+mY2NhYSJKktJWUlNQrPl1gQkRERGQsDDAP0Z49ezB//ny8+eabSE9Px6BBgzBy5EhkZ6t+0aS0tBT29vZ488030bt37xrrtba2Rm5ursJmaWmpcXy6woSIiIiIarRhwwZMnz4dM2bMQPfu3REdHQ0XFxds3bpV5fGdO3fGpk2bEB4eDhsbmxrrlSRJPq/g/fMLGgoTIiIiImNRPahamw1Vr/Hfv92/nuf9ysrKkJaWhpCQEIXykJAQnDhxQqumFBcXw9XVFR07dsSTTz6J9PR0rerTFhMiIiIio6GbMUQuLi6wsbGRb9WrPTyooKAAFRUVcHBwUCh3cHBAXl5evVvh5eWF2NhYHDp0CLt27YKlpSWCgoJ0upappvjaPRERUTOTk5OjMDFjXQucSw+MPRJCKJVpwt/fH/7+/vLPQUFB6Nu3L7Zs2YLNmzfXu15tMCEiIiIyFtou0Prfc62trdWaqbpdu3YwNTVV6g3Kz89X6jXShomJCfr372/QHiI+MiMiIjIaJjrY1Gdubg4/Pz8kJiYqlCcmJiIwMFCbhigQQiAjIwNOTk46q1NT7CEiIiIyFhK07CHS/JSFCxciLCwM/fr1Q0BAALZt24bs7GzMnj0bABAREYHr168jLi5Ofk5GRgaAqoHTf/31FzIyMmBubg5vb28AwIoVK+Dv74+uXbtCJpNh8+bNyMjIwPvvv1//tmmJCRERERHVKDQ0FIWFhVi5ciVyc3PRs2dPHD58GK6urgCqJmJ8cE4iX19f+d9paWn4/PPP4erqiqtXrwIAbt68iVmzZiEvLw82Njbw9fVFSkoKBgwY0GDtehBXu1eBq90r4mr3VbjafRWudk+kqEFXu7+ZBGvr1lrUUwybh4by900F9hAREREZCx0NqiZlTIhUqO40k8lkBo6kcRC3ywwdQqMg7pYbOoRGwYT/XhApqP6taIgHLjLZbYOe35QxIVLh1q1bAKomriKiByyqeSp+oubs1q1btS5VoQ1zc3M4OjrCxeUJretydHSEubm5DqJqWjiGSIXKykrcuHEDbdq0qdfEUzKZDC4uLkoTXzVFbGvT01zaCbCtTVVDt1UIgVu3bsHZ2RkmJvqbzaakpARlZdr32Jubmxt0EdXGij1EKpiYmKBjx45a16PuxFdNAdva9DSXdgJsa1PVkG3VV8/Q/SwtLZnI6BEnZiQiIqJmjwkRERERNXtMiPTAwsICy5cvr3OxvKaAbW16mks7Aba1qWpObSXd4aBqIiIiavbYQ0RERETNHhMiIiIiavaYEBEREVGzx4SIiIiImj0mRBpKSUnBqFGj4OzsDEmS8OWXX9Z5ztGjR+Hn5wdLS0u4u7vjgw8+0H+gOqBpW5OTkyFJktJ28eLFhgm4nqKiotC/f3+0adMG7du3x9ixY3Hp0qU6zzPG+1qfthrrfd26dSt69eoln5wvICAA3377ba3nGOM9BTRvq7He0wdFRUVBkiTMnz+/1uOM9b5Sw2JCpKHbt2+jd+/eeO+999Q6PisrC48//jgGDRqE9PR0vPHGG5g7dy727dun50i1p2lbq126dAm5ubnyrWvXrnqKUDeOHj2Kl19+GSdPnkRiYiLKy8sREhKC27drXgTRWO9rfdpazdjua8eOHfH222/j9OnTOH36NB599FGMGTMGv/zyi8rjjfWeApq3tZqx3dP7/fTTT9i2bRt69epV63HGfF+pgQmqNwDiwIEDtR6zePFi4eXlpVD2wgsvCH9/fz1GpnvqtDUpKUkAEP/880+DxKQv+fn5AoA4evRojcc0lfuqTlubyn0VQoi2bduKjz76SOW+pnJPq9XWVmO/p7du3RJdu3YViYmJYsiQIWLevHk1HtvU7ivpD3uI9Cw1NRUhISEKZSNGjMDp06dx7949A0WlX76+vnByckJwcDCSkpIMHY7GioqKAAC2trY1HtNU7qs6ba1mzPe1oqICu3fvxu3btxEQEKDymKZyT9VpazVjvacvv/wynnjiCQwbNqzOY5vKfSX94+KuepaXlwcHBweFMgcHB5SXl6OgoABOTk4Gikz3nJycsG3bNvj5+aG0tBSffvopgoODkZycjMGDBxs6PLUIIbBw4UI8/PDD6NmzZ43HNYX7qm5bjfm+nj9/HgEBASgpKUHr1q1x4MABeHt7qzzW2O+pJm015nu6e/dunDlzBj/99JNaxxv7faWGw4SoAUiSpPBZ/Hdy8AfLjV23bt3QrVs3+eeAgADk5OTg3XffbfT/ka02Z84cnDt3DsePH6/zWGO/r+q21Zjva7du3ZCRkYGbN29i3759mDp1Ko4ePVpjomDM91STthrrPc3JycG8efOQkJCg0arvxnxfqeHwkZmeOTo6Ii8vT6EsPz8fZmZmsLOzM1BUDcff3x+ZmZmGDkMtr7zyCg4dOoSkpCR07Nix1mON/b5q0lZVjOW+mpubo0uXLujXrx+ioqLQu3dvbNq0SeWxxn5PNWmrKsZwT9PS0pCfnw8/Pz+YmZnBzMwMR48exebNm2FmZoaKigqlc4z9vlLDYQ+RngUEBOCrr75SKEtISEC/fv3QokULA0XVcNLT0xt9l7QQAq+88goOHDiA5ORkuLm51XmOsd7X+rRVFWO4r6oIIVBaWqpyn7He05rU1lZVjOGeBgcH4/z58wplzz33HLy8vPD666/D1NRU6Zymdl9Jjww1mttY3bp1S6Snp4v09HQBQGzYsEGkp6eLP/74QwghxJIlS0RYWJj8+CtXrohWrVqJBQsWiF9//VXs2LFDtGjRQnzxxReGaoLaNG3rxo0bxYEDB8Rvv/0mfv75Z7FkyRIBQOzbt89QTVDLiy++KGxsbERycrLIzc2Vb3fu3JEf01Tua33aaqz3NSIiQqSkpIisrCxx7tw58cYbbwgTExORkJAghGg691QIzdtqrPdUlQffMmtK95UaFhMiDVW/rvrgNnXqVCGEEFOnThVDhgxROCc5OVn4+voKc3Nz0blzZ7F169aGD7weNG3runXrhIeHh7C0tBRt27YVDz/8sPjmm28ME7wGVLURgPjkk0/kxzSV+1qfthrrfX3++eeFq6urMDc3F/b29iI4OFieIAjRdO6pEJq31VjvqSoPJkRN6b5Sw5KE+O/oMiIiIqJmioOqiYiIqNljQkRERETNHhMiIiIiavaYEBEREVGzx4SIiIiImj0mRERERNTsMSEiIiKiZo8JEZGeXL16FZIkISMjQy/1S5KEL7/8st7nJycnQ5IkSJKEsWPH1nrsI488gvnz59f7u0jZtGnT5Ndfm/tIRLrBhIiapGnTptX5I69vLi4uyM3NRc+ePQH8LwG5efOmQeN60KVLlxAbG2voMJqsmu77pk2bkJuba5igiEgJF3cl0hNTU1M4OjoaOow6tW/fHg899JChw8C9e/ea1WKbNjY2sLGxMXQYRPRf7CGiZuno0aMYMGAALCws4OTkhCVLlqC8vFy+/5FHHsHcuXOxePFi2NrawtHREZGRkQp1XLx4EQ8//DAsLS3h7e2NI0eOKDz+uP+R2dWrVzF06FAAQNu2bSFJEqZNmwYA6Ny5M6KjoxXq7tOnj8L3ZWZmYvDgwfLvSkxMVGrT9evXERoairZt28LOzg5jxozB1atXNb42t2/fRnh4OFq3bg0nJyesX79e6ZiysjIsXrwYHTp0gJWVFQYOHIjk5GSFY7Zv3w4XFxe0atUK48aNw4YNGxQSr8jISPTp0wcff/wx3N3dYWFhASEEioqKMGvWLLRv3x7W1tZ49NFHcfbsWYW6v/rqK/j5+cHS0hLu7u5YsWKFwv2LjIxEp06dYGFhAWdnZ8ydO7fWNtdV34YNG+Dj4wMrKyu4uLjgpZdeQnFxsXz/H3/8gVGjRqFt27awsrJCjx49cPjw4VrvOxE1Luwhombn+vXrePzxxzFt2jTExcXh4sWLmDlzJiwtLRWSkJ07d2LhwoU4deoUUlNTMW3aNAQFBWH48OGorKzE2LFj0alTJ5w6dQq3bt3Cq6++WuN3uri4YN++fXjqqadw6dIlWFtbo2XLlmrFW1lZifHjx6Ndu3Y4efIkZDKZ0nieO3fuYOjQoRg0aBBSUlJgZmaG1atX47HHHsO5c+dgbm6u9vV57bXXkJSUhAMHDsDR0RFvvPEG0tLS0KdPH/kxzz33HK5evYrdu3fD2dkZBw4cwGOPPYbz58+ja9eu+OGHHzB79mysW7cOo0ePxpEjR7Bs2TKl77p8+TL27t2Lffv2wdTUFADwxBNPwNbWFocPH4aNjQ0+/PBDBAcH47fffoOtrS3i4+Px7LPPYvPmzRg0aBB+//13zJo1CwCwfPlyfPHFF9i4cSN2796NHj16IC8vTymhul9d9QGAiYkJNm/ejM6dOyMrKwsvvfQSFi9ejJiYGADAyy+/jLKyMqSkpMDKygq//vorWrdurdV9J6IGZuDFZYn0YurUqWLMmDEq973xxhuiW7duorKyUl72/vvvi9atW4uKigohRNUK2g8//LDCef379xevv/66EEKIb7/9VpiZmYnc3Fz5/sTERAFAHDhwQAghRFZWlgAg0tPThRBCJCUlCQDin3/+UajX1dVVbNy4UaGsd+/eYvny5UIIIeLj44WpqanIycmR7//2228VvmvHjh1KbSotLRUtW7YU8fHxKq+Dqnhu3bolzM3Nxe7du+VlhYWFomXLlvIVxS9fviwkSRLXr19XqC84OFhEREQIIYQIDQ0VTzzxhML+Z555RtjY2Mg/L1++XLRo0ULk5+fLy7777jthbW0tSkpKFM718PAQH374oRBCiEGDBom1a9cq7P/000+Fk5OTEEKI9evXC09PT1FWVqay3Q+qqz5V9u7dK+zs7OSffXx8RGRkpMpja7rv1e6/j0RkOOwhombnwoULCAgIgCRJ8rKgoCAUFxfj2rVr6NSpEwCgV69eCuc5OTkhPz8fQNVAZBcXF4UxQgMGDNBbvJ06dULHjh3lZQEBAQrHpKWl4fLly2jTpo1CeUlJCX7//Xe1v+v3339HWVmZQv22trbo1q2b/POZM2cghICnp6fCuaWlpbCzswNQdX3GjRunsH/AgAH4+uuvFcpcXV1hb2+v0I7i4mJ5PdXu3r0rb0daWhp++uknrFmzRr6/oqICJSUluHPnDp5++mlER0fD3d0djz32GB5//HGMGjUKZmaq/3NXV32tWrVCUlIS1q5di19//RUymQzl5eUoKSnB7du3YWVlhblz5+LFF19EQkIChg0bhqeeekrpnx8iatyYEFGzI4RQSIaqywAolD84wFeSJFRWVtZYR32ZmJjIv7/avXv3lGJ7MJb7VVZWws/PD5999pnSsfcnHHVR9V0PqqyshKmpKdLS0uSPuaq1bt1aXk9N1/h+VlZWSnU7OTkpjUcCIB9/VFlZiRUrVmD8+PFKx1haWsLFxQWXLl1CYmIijhw5gpdeegnvvPMOjh49qnLQdl31/fHHH3j88ccxe/ZsrFq1Cra2tjh+/DimT58uv08zZszAiBEj8M033yAhIQFRUVFYv349XnnlFaU6iahxYkJEzY63tzf27dun8KN94sQJtGnTBh06dFCrDi8vL2RnZ+PPP/+Eg4MDAOCnn36q9ZzqcTwVFRUK5fb29gqvX8tkMmRlZSnEm52djRs3bsDZ2RkAkJqaqlBH3759sWfPHvlA5Prq0qULWrRogZMnT8p7yv755x/89ttvGDJkCADA19cXFRUVyM/Px6BBg1TW4+XlhR9//FGh7PTp03V+f9++fZGXlwczMzN07ty5xmMuXbqELl261FhPy5YtMXr0aIwePRovv/wyvLy8cP78efTt21fj+k6fPo3y8nKsX78eJiZV76Hs3btX6TgXFxfMnj0bs2fPRkREBLZv345XXnmlxvtORI0L3zKjJquoqAgZGRkKW3Z2Nl566SXk5OTglVdewcWLF3Hw4EEsX74cCxculP/g1WX48OHw8PDA1KlTce7cOfzwww948803ASj33lRzdXWFJEn4+uuv8ddff8nfUnr00Ufx6aef4tixY/j5558xdepUhZ6XYcOGoVu3bggPD8fZs2dx7Ngx+XdVe+aZZ9CuXTuMGTMGx44dQ1ZWFo4ePYp58+bh2rVral+z1q1bY/r06Xjttdfw3Xff4eeff8a0adMUrounpyeeeeYZhIeHY//+/cjKysJPP/2EdevW4fDhwwCAV155BYcPH8aGDRuQmZmJDz/8EN9++22dvWrDhg1DQEAAxo4di/j4eFy9ehUnTpzA0qVL5QnVW2+9hbi4OERGRuKXX37BhQsXsGfPHixduhQAEBsbix07duDnn3/GlStX8Omnn6Jly5ZwdXVV+Z111efh4YHy8nJs2bJFXt8HH3ygUMf8+fMRHx+PrKwsnDlzBt9//z26d+8OoOb7TkSNjKEGLxHp09SpUwUApW3q1KlCCCGSk5NF//79hbm5uXB0dBSvv/66uHfvnvz8IUOGyAcRVxszZoz8fCGEuHDhgggKChLm5ubCy8tLfPXVVwKA+M9//iOEUB5ULYQQK1euFI6OjkKSJHldRUVFYuLEicLa2lq4uLiI2NhYhUHVQghx6dIl8fDDDwtzc3Ph6ekp/vOf/ygNxs3NzRXh4eGiXbt2wsLCQri7u4uZM2eKoqIildeopsG+t27dEs8++6xo1aqVcHBwEP/617+UrkdZWZl46623ROfOnUWLFi2Eo6OjGDdunDh37pz8mG3btokOHTqIli1birFjx4rVq1cLR0dH+f7ly5eL3r17K8Ulk8nEK6+8IpydnUWLFi2Ei4uLeOaZZ0R2drb8mP/85z8iMDBQtGzZUlhbW4sBAwaIbdu2CSGEOHDggBg4cKCwtrYWVlZWwt/fXxw5ckTlNVCnPiGE2LBhg3BychItW7YUI0aMEHFxcQrXbs6cOcLDw0NYWFgIe3t7ERYWJgoKCuTnq7rv1R68j0RkGJIQagwaIKI6/fDDD3j44Ydx+fJleHh4GDqcOiUnJ2Po0KH4559/GmRixpkzZ+LixYs4duyY3r/LmEiShAMHDhh8ZnWi5o6PzIjq6cCBA0hMTMTVq1dx5MgRzJo1C0FBQUaRDN2vY8eOmDx5ss7rfffdd3H27FlcvnwZW7Zswc6dOzF16lSdf4+xmj17tnwQOhEZHnuIiOopLi4Oq1atQk5ODtq1a4dhw4Zh/fr1Sq+MN1Z3797F9evXAVSNHdL1MiMTJ05EcnIybt26BXd3d7zyyiuYPXu2Tr/DmOXn50MmkwGomtLhwTfuiKhhMSEiIiKiZo+PzIiIiKjZY0JEREREzR4TIiIiImr2mBARERFRs8eEiIiIiJo9JkRERETU7DEhIiIiomaPCRERERE1e0yIiIiIqNn7f+klOegsLaX8AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ss.rmse.plot(model='SW_1', cmap='YlOrRd');" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "base", "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.12" }, "orig_nbformat": 2 }, "nbformat": 4, "nbformat_minor": 2 }