{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Get Anomaly Data\n", "\n", "## get_anom\n", "Get climate data anomaly\n", "\n", "### Parameter:\n", "+ DaArray (xr.DataArray): shape = (time, *number) original Dataarray\n", "\n", "+ method (int, optional): method of getting anomaly.\n", " * 0 is Minus the multi-year average of the corresponding month\n", " * 1 is Remove the linear trend of the corresponding month\n", " * Defaults to 0. \n", " \n", "+ Returns:\n", " anom (xr.DataArray): climate data anomaly\n", "\n", "### Example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Load Module and Date" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import sacpy as scp\n", "\n", "sst = scp.load_sst()['sst']\n", "ssta0 = scp.get_anom(sst,method=0)\n", "ssta1 = scp.get_anom(sst,method=1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### plot and see difference between different methods" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkUAAAHFCAYAAAD8Jo2EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABTeElEQVR4nO3df1xUVf4/8NcdYIbfKCIMCKL5A39gxmop1gZq4o80lWozd00rLTPrY65rue03sV3FzFx318p+rVjq1m6o/VATTdBMLTEtszRLFFIIFeQ3DMyc7x8us47MGWBmGO7A6/l43IfOfd977rl37sy8OefecxUhhAARERFRO6dp7QoQERERqQGTIiIiIiIwKSIiIiICwKSIiIiICACTIiIiIiIATIqIiIiIADApIiIiIgLApIiIiIgIAJMiIiIiIgBMilzmwIEDSElJwZUrVxrEEhMTkZiY6PI6ucIPP/yABQsWYNCgQejQoQOCg4Nx66234v3337e6fGFhIWbMmIGQkBD4+voiPj4en376qcUypaWlWLp0KRITE6HX6+Hv748BAwbghRdeQHV1dYMy//SnP2H8+PHo0qULFEXBjBkzmrUPP//8M+bNm4eEhAR06NABiqIgLS2twXLNrZfM6tWrkZycjO7du0NRlEbPjQ8++AAJCQkIDAyEn58f+vfvj9dffx0AkJKSAkVRGp3qt3HixAnMmTMH8fHx8PPzg6IoyMrKarDN/Px8/OlPf0J8fDxCQkIQGBiIQYMG4fXXX4fRaGzyvp45cwbJycno0KED/P39MWrUKHz11VdWl3333Xdx0003wdvbGxEREZg3bx7Ky8tVuS13sGzZMmzdurXB/LS0NCiKguzsbNdX6jplZWVYuHAhkpKS0LlzZyiKgpSUlAbLGY1GrFq1CmPGjEFkZCR8fX3Rt29fPPPMM1a/c4mkBLnEiy++KACInJycBrETJ06IEydOuL5SLvCPf/xD9OnTRyxdulRkZGSI7du3i+nTpwsAYsmSJRbLVldXi9jYWBEZGSk2bNggMjIyxMSJE4Wnp6fIysoyL3f8+HEREhIinnrqKfHBBx+ITz/9VKSkpAhvb28xcuRIYTKZLMr19fUVQ4cOFbNnzxZarVZMnz69WfuQmZkpQkJCxB133CHuv/9+AUCsW7euwXLNrZdMTEyM+NWvfiUeeugh0blzZ5GQkCBdNjU1VWg0GjFnzhyxY8cOsXv3brFmzRrxj3/8QwghRF5enjh48KB52rx5swAgnnjiCYv59edfWlqaCA8PF+PGjRMTJkwQAERmZmaD7X700UciKipKPPvss2Lbtm0iIyNDPPXUU0Kj0YgHH3ywSftZWFgoIiIiRP/+/UV6errYtm2buO2220RAQIA4efKkxbIbNmwQAMTMmTPFnj17xNq1a0VQUJAYNWqU6rblLvz8/Kx+FtatWycAiMOHD7u+UtfJyckRQUFB4vbbbxczZ84UAMTixYsbLFdWViYCAgLEI488Iv7zn/+IzMxM8dJLL4mOHTuKfv36icrKStdXntwSkyIXsZUUtWUXL160mgzceeedwtfXV1RXV5vnvfzyywKAOHDggHlebW2t6Nevn7jlllvM88rLy0V5eXmDMuuP8WeffWYx32g0mv8v+yGw5dr1Dx8+LE2Kmluvpmyvf//+0qQoOztbaDQa8cILLzSpXCGu/sgAEC+++GKj2/7Pf/4jTYqKioqEwWBoMP/xxx8XAERubm6jdfnDH/4gvLy8xNmzZ83zSkpKREhIiPjNb35jnldXVyfCw8NFUlKSxfobN24UAMT27dtVtS134Q5JkclkMn9/XLx4UZoU1dXViUuXLjWYX38Ov/POOy1dVWoj2H3mAikpKfjDH/4AAOYukWu7Ja7vPjt79iwURcGLL76IF154Ad26dYOPjw8SExPxww8/oLa2Fs888wwiIiIQFBSEyZMno7CwsMF233vvPXM3iL+/P0aPHo2jR4+6YpfNQkJCoChKg/m33HILKisrUVRUZJ63ZcsWxMTEID4+3jzP09MTv/vd7/Dll1/i/PnzAAA/Pz/4+flZLRMA8vLyLOZrNI6d5k1dv7n1cnR7a9asgU6nwxNPPNGk5Z257Y4dO8LLy6vB/Pp9/fnnnxstY8uWLRgxYgSio6PN8wIDA5GcnIyPPvoIdXV1AIBDhw4hPz8fDz74oMX69957L/z9/bFlyxZVbcuabt26Yfz48fj4448RFxcHHx8f9O3bFx9//DGAq11Wffv2hZ+fH2655RarXVcffvgh4uPj4evri4CAAIwaNQoHDx60WKa+u/TEiRO4//77ERQUhLCwMDz00EMoKSkxL6coCioqKrB+/foGXaj1ysrK8NhjjyEkJASdOnVCcnIyLly4YNf+26u+bo3x8PBAp06dGsxv7mePiEmRC8ycOdP8w7V582YcPHgQBw8exK9+9Sub67388sv4/PPP8fLLL+PNN9/EyZMnMWHCBDz88MO4ePEi/vnPf2LFihXYvXs3Zs6cabHusmXLcP/996Nfv37497//jXfeeQdlZWX49a9/je+++67ROtfV1TVpEkLYdUwyMzPRuXNnhIaGmud9++23uPHGGxssWz/vxIkTNsvcs2cPAKB///521amltFS99u3bh759+yI9PR0xMTHw8PBAZGQknnnmGRgMBqduq6n27NkDT09P9O7d22L+9T+6VVVV+Omnn6Tvd1VVFc6cOQPg6nlRP/9aXl5e6NOnjzneGttqjq+//hqLFi3C008/jc2bNyMoKAjJyclYvHgx3nzzTSxbtgwbN25ESUkJxo8fj6qqKvO6mzZtwsSJExEYGIh//etfeOutt1BcXIzExETs37+/wbbuvvtu9O7dG+np6XjmmWewadMmPPXUU+b4wYMH4ePjg3Hjxpm/j1555RWLMmbOnAkvLy9s2rQJK1asQFZWFn73u981up9CiCZ/f7Q0tX4nkHp5tnYF2oPIyEh07doVABAXF4du3bo1ab0OHTpg69at5r/eL126hHnz5qFPnz744IMPzMudPHkSq1evRmlpKQIDA5GXl4fFixdj7ty5+Pvf/25ebtSoUejVqxeWLFmC9957T7rds2fPonv37k2qY2ZmZrMvEn/zzTeRlZWFv/3tb/Dw8DDPv3z5MoKDgxssXz/v8uXL0jK/+eYbrFixApMnT7b649daWrJe58+fx8WLF/Hkk0/iz3/+M/r164dPP/0Uy5cvR15eHjZu3OjU7TUmIyMD77zzDv7v//6vwV/tHh4eFu91cXExhBBNer/r/5Ute/bs2VbbVnNcvnwZhw4dQpcuXQAAERERuOmmm/DGG2/gxx9/hK+vL4CrSd2kSZOwe/duTJgwASaTCX/4wx8wYMAA7Nixw/x9MG7cOPTo0QNPP/00Pv/8c4ttPfzww+bW6TvuuAM//vgj/vnPf+Ktt96CoigYOnQoNBoNOnfujKFDh1qt75gxYyy+P4qKirBw4UIUFBRAr9dL93P9+vUNWtpk7P2jqinOnz+PZ555BoMHD8b48eNbbDvUtjApUrFx48ZZdGf07dsXAHDnnXdaLFc/Pzc3F7Gxsdi5cyfq6urwwAMPWPw15u3tjYSEBGRmZtrcbkREBA4fPtykOsbExDRpuXo7duzA448/jnvuucdqt4+tpnJZ7OzZsxg/fjyioqLw5ptvNqs+9YQQDe6a8vR07ONhq17X/5Xs4eHRpG6Ca5lMJpSVleFf//oXpkyZAgAYPnw4KioqsHr1aixZsgQ9e/Z0aB+a6quvvsJvfvMbDB06FKmpqQ3islaB5rzfsmWvn+/KbTXHTTfdZE6IgP99bhMTE80J0bXzz507BwA4deoULly4gHnz5ll8H/j7++Puu+/Ga6+9hsrKSosy7rrrLott33jjjaiurkZhYSHCwsKaVF9rZdTXy1ZSNGHChCZ/f7SUoqIijBs3DkIIvPfeew53oVP7waRIxa7/a1Wr1dqcX3/b9y+//AIAuPnmm62W29gXhFarxU033dSkOl77F3ljdu7cieTkZIwaNQobN25s8APTqVMnq61B9dcdWfvr/dy5cxg+fDg8PT3x6aefWl2mKaz9devIX7GN1ev663HWrVvX7KECOnXqhIKCAowePdpi/tixY7F69Wp89dVXLkmKjh49am6F3L59O3Q6XaPrdOzYEYqiNOn9rm91unz5coMf9KKiokbfc1duyxZ7P8/19Q4PD29QZkREBEwmE4qLiy2Soutb6urfk2u75BpjbxnBwcEICgpq8nacrbi4GKNGjcL58+exZ88e3HDDDa1WF3I/TIraoJCQEADA+++/b3FhaVO1RPfZzp07MWnSJCQkJCA9Pd38xX+tAQMG4Pjx4w3m18+LjY21mH/u3DkkJiZCCIGsrCxERkY2qc7WOPOv26bU6/ptNfV4X+vGG29EQUFBg/n1yZwr/jo+evQo7rjjDkRHRyMjI6PJP4Y+Pj7o2bOn9P328fEx/5gNGDDAPL9fv37m5erq6nDy5Encf//9qtlWS6hPTvLz8xvELly4AI1Gg44dO7q6WlKt2X1WXFyMO+64Azk5Ofj0009V1ZVO7oFJkYvY85eavUaPHg1PT0/89NNPuPvuu5u9vrO7zzIyMjBp0iTcdttt2Lp1q7QlYfLkyZgzZw6++OILDBkyBMDVH6MNGzZgyJAhiIiIMC+bm5uLxMREGI1GZGVl2ZX8XatTp05W715prqbWa/DgwQ5v6+6770ZGRgZ27NiBqVOnmudv374dGo1G2lLoLMeOHcMdd9yByMhI7Nq1q9k/zJMnT8bq1auRl5eHqKgoAFfveNq8eTPuuusuc/flkCFDEB4ejrS0NNx3333m9d9//32Ul5cjOTlZVdtytpiYGHTp0gWbNm3CggULzC2sFRUVSE9PN9+R1lw6na5Fvo9aq/usPiE6c+YMdu3ahbi4OJfXgdwfkyIXqf8L9G9/+xumT58OLy8vxMTEICAgwOnb6tatG55//nk8++yzOHPmDMaMGYOOHTvil19+wZdffgk/Pz8sWbJEur5Wq3XKjzYA7N+/H5MmTYJer8cf//hHHDt2zCLer18/BAYGAgAeeughvPzyy7j33nuxfPlyhIaG4pVXXsGpU6ewe/du8zqFhYUYPnw48vPz8dZbb6GwsNBiSILIyEiL1pm9e/fi4sWLAK6OfHvu3DnziNoJCQno3Llzo/tRv3z9XUrZ2dnw9/cHANxzzz121UsmOzvbfEFvaWkphBDm7d98883mROvBBx/Ea6+9hjlz5uDSpUvo168fdu/ejZdffhlz5syxK1GsrKzE9u3bAVy9PR24evwuXboEPz8/jB07FsDV61zuuOMOAMDSpUtx+vRpnD592lxOjx49LI6rp6cnEhISLEYnX7BgAd555x3ceeedeP7556HT6bB8+XJUV1dbjFrs4eGBFStWYNq0aXj00Udx//334/Tp01i4cCFGjRqFMWPGWOyDK7elKAoSEhKsjvrtLBqNBitWrMBvf/tbjB8/Ho8++ihqamrw4osv4sqVK1i+fLld5Q4YMABZWVn46KOPEB4ejoCAgGZfI2iNs/7AqLdjxw5UVFSgrKwMAPDdd9+ZPw/jxo2Dr68vqqqqzEOOrF69GnV1debzFwA6d+6MHj16OK1O1Ia1yuhI7dSiRYtERESE0Gg0FoPiJSQkWAzQJxtgLzMzUwAQ//nPfyzmywZb27p1qxg+fLgIDAwUOp1OREdHi3vuuUfs3r27RfbPmsWLFwsA0un6gQELCgrEAw88IIKDg4W3t7cYOnSo2LVrl8Uy9cdBNl0/uFtCQkKTty9ja3v21kumfsRva9P1g0ZevnxZPProoyIsLEx4eXmJ3r17ixdffNFiEMZrNTZ4Y33c2hQdHW1erv6ca2o9AVgdhPLHH38UkyZNEoGBgcLX11eMHDlSHDlyxGrdNm3aJG688Uah1WqFXq8XTz75pCgrK2uwnKu2VVZWJgCIKVOmWC3jWtHR0eLOO++0WtfHH3/cYp7sPdq6dasYMmSI8Pb2Fn5+fmLkyJHi888/t1im/vN28eJFi/n179e1g8ceO3ZM3HrrrcLX19fimMm+T+rP76Z+ZpwlOjpaep7V74+t8xZAswdspfZLEaIF74kkImqjtm/fjvHjx+Prr782twQTkXvjfYpERHbIzMzElClTmBARtSFsKSIiIiICW4qIiIiIADApIiIiIgLApIiIiIgIAJMiIiIiIgAcvLEBk8mECxcuICAgwKGHPxIRUdsmhEBZWRkiIiJa9LE61dXVMBgMTilLq9XC29vbKWW1RUyKrnPhwgXzYwCIiIgak5eX59CzF22prq5GJx9/VMLolPL0ej1ycnKYGEkwKbpO/WM38vLyzI+fILJH5ZbV0pjn7fdKYz8Y/KUxvZ+XNGZrdI0O1b9IY4pJ/mWrlF2UxlBn/S/Xmu+PSFe5+PUP0ljpuUvSWNnPpdLYDePkD/28cipPGou4faA0pvHykNflnPxYevlZf64fAPj1kD+tXeMnf5Cups9QaczkLV+vsFb+9R4VLD/H2oJzzz0qjRmulEtjneN6W52v6xxidX5ZZTV6zXq+RR7XVM9gMKASRvwWXaB18IoXA0zYWHAeBoOBSZEEk6Lr1HeZBQYGMikih3j6yr90PG18ifob5LEAf/uSokCvSmlMMdXJY0K+Huqs16XGxn5X67TSmMlL/nVk8pTHArzlZdZp5ccr0EeewGi0Nr4abWzPy1tepr+vj3x7fvIHumpsnCsmH/l3VJWNpCgwsG0nRQE2zrMaO84JnY1zGoBLLrXwgQZaxbGkyIOjEjaKSREREZHKeSgKPBxMvjygXH0aHEkxKSIiIlI5jQJ4ONggpQGYFDWCt+QTERERgS1FREREque07jOyiUkRERGRynk4oftMfk8l1WP3GRERETWwb98+TJgwAREREVAUBVu3brW5fFZWFhRFaTCdPHnSNRV2ArYUERERqVxrdJ9VVFRg4MCBePDBB3H33Xc3eb1Tp05ZDGnTuXPnZm23NblNS9Grr76KG2+80Tx+UHx8PHbs2GGOCyGQkpKCiIgI+Pj4IDExESdOnGjFGhMRETlHffeZo1NzjB07Fn/5y1+QnJzcrPVCQ0Oh1+vNk4eH+3TcuU1SFBkZieXLlyM7OxvZ2dkYMWIEJk6caE58VqxYgVWrVmHNmjU4fPgw9Ho9Ro0ahbKyslauORERkXqUlpZaTDU1NU4tPy4uDuHh4Rg5ciQyMzOdWnZLc5ukaMKECRg3bhx69+6N3r17Y+nSpfD398ehQ4cghMDq1avx7LPPIjk5GbGxsVi/fj0qKyuxadOm1q46ERGRQ+q7zxydACAqKgpBQUHmKTU11Sl1DA8Px+uvv4709HRs3rwZMTExGDlyJPbt2+eU8l3BLa8pMhqN+M9//oOKigrEx8cjJycHBQUFSEpKMi+j0+mQkJCAAwcO4NFH5c/BqampsciSS0vlz1giIiJqDQocb8Wo7z27/tmeOp388TTNERMTg5iYGPPr+Ph45OXlYeXKlbj99tudso2W5jYtRQBw/Phx+Pv7Q6fTYfbs2diyZQv69euHgoICAEBYWJjF8mFhYeaYTGpqqkXGHBUV1WL1JyIiam311+bWT85KiqwZOnQoTp8+3WLlO5tbJUUxMTE4duwYDh06hMceewzTp0/Hd999Z45f/1A+IUSjD+pbtGgRSkpKzFNenvyp2kRERK3Bmd1nrnT06FGEh4e7fLv2cqvuM61Wi549ewIABg8ejMOHD+Nvf/sbnn76aQBAQUGBxcEvLCxs0Hp0PZ1O16JZMhERkaNaY/DG8vJy/Pjjj+bXOTk5OHbsGIKDg9G1a1csWrQI58+fx9tvvw0AWL16Nbp164b+/fvDYDBgw4YNSE9PR3p6umMVdyG3SoquJ4RATU0NunfvDr1ej127diEuLg4AYDAYsHfvXrzwwgutXEsiIiLHXE2KHB2nqHmys7MxfPhw8+v58+cDAKZPn460tDTk5+cjNzfXHDcYDFiwYAHOnz8PHx8f9O/fH9u2bcO4ceMcqrcruU1S9Mc//hFjx45FVFQUysrK8O677yIrKwuffPIJFEXBvHnzsGzZMvTq1Qu9evXCsmXL4Ovri6lTp7Z21YmIiNxOYmIihBDSeFpamsXrhQsXYuHChS1cq5blNknRL7/8gmnTpiE/Px9BQUG48cYb8cknn2DUqFEArr4ZVVVVmDNnDoqLizFkyBBkZGQgICCglWtORETkGD77zDXcJil66623bMYVRUFKSgpSUlJcUyEiIiIXaY3HfLRHbnX3GREREVFLcZuWIiIiovZK44TuM7aCNI5JERERkcqx+8w1mDgSERERgS1FREREqse7z1yDSREREZHKMSlyDXafEREREYEtRURERKrHC61dg0kRmeVfqZDGwjv4ubAmbYOu/y3SWJE2WBoL1pjs2l4HVEljIueYNKYJ7yGN1UTfLI15Vl+xOt9b6y1dJ2rYeGmsS0GONObRKUIaE9Xl0ljn6fKnc2vKL8nL9PKRxryH2/gsKPLG94rASGmsziR/lIKtnzEvG/0pIZ7yWGGJ/LMeGuT8z/p7of3k2+sif+pA7APDpLHOT/3Vrrp06NlFGtN2kn8uW5sHnNB9Jj/N6L+YFBEREamcxgktRRoH128PeE0REREREdhSREREpHpOufuMDUWNYlJERESkck650JrdZ41i9xkRERER2FJERESkeuw+cw0mRURERCrH7jPXYPcZEREREdhSREREpHoaRXF4nCGOU9Q4JkVEREQqp3goUDSOJTUKk6JGsfuMiIiICGwpIiIiUj2NhwKNgy1F7D5rHJMiIiIitfPQQNE42Lmj8ImwjWFSREREpHKKRoHi4EBDCthS1BheU0REREQEthQRERGpnsZDgcbBliINW4oaxaSIiIhI5RSN49cUKYLXFDWG3WdEREREYEsRERGR6rH7zDWYFBEREamc4sG7z1yB3WdEREREYEsRERGR6l1tKXLwQmuYnFSbtotJERERkcrxmiLXYPcZEREREdhSREREpHqKokBx8IGwioktRY1hUkRERKRyGg8NNA5eU6QR7BxqDJMiIiIilXPKLfmCLUWNYVJEqlRUVimNBQf4urAmwNmnp0tjdRVV0ljU8jekMaNRPty+v1b+15yvUX5chEb+cVbqauXreemkMd0vJ+XraX2sz/cJlK/j6S2NGfuPlMbKjPIvc4PJvkcXeAdGSWMeNropbG3OYON9vVJRJ40F6eTveZBSI415FubKK2ODJixGGqv7OkO+vYFJdm3vvsLv7Fqv7ugn8tiRbdJYtxfW27c9SZmKn/VzWlteYdd2SL2YFBEREakcW4pcg0kRERGRyvGaItfgESIiIiICW4qIiIjUzwndZ2D3WaOYFBEREamcRlGgcXCcIo3CpKgx7D4jIiIiAluKiIiIVE/x0Dj+QFgT20EawyNERESkcvUPhHV0ao59+/ZhwoQJiIiIgKIo2Lp1a6Pr7N27F4MGDYK3tzduuOEGrF271s49bh1MioiIiKiBiooKDBw4EGvWrGnS8jk5ORg3bhx+/etf4+jRo/jjH/+IJ598Eunp6S1cU+dxm6QoNTUVN998MwICAhAaGopJkybh1KlTFssIIZCSkoKIiAj4+PggMTERJ06caKUaExEROUf94I2OTs0xduxY/OUvf0FycnKTll+7di26du2K1atXo2/fvpg5cyYeeughrFy50p5dbhVukxTt3bsXjz/+OA4dOoRdu3ahrq4OSUlJqKj43zDrK1aswKpVq7BmzRocPnwYer0eo0aNQllZWSvWnIiIyDH11xQ5OrWkgwcPIinJ8lEwo0ePRnZ2Nmpr5Y8aUhO3udD6k08sn4Gzbt06hIaG4siRI7j99tshhMDq1avx7LPPmrPa9evXIywsDJs2bcKjjz7aGtUmIiJymMYDzb4mqEEZpqv/lpaWWszX6XTQ6eTPQGyqgoIChIWFWcwLCwtDXV0dLl26hPDwcIe30dLcpqXoeiUlJQCA4OBgAFf7MgsKCiyyVJ1Oh4SEBBw4cEBaTk1NDUpLSy0mIiKitioqKgpBQUHmKTU11WllK9eNhSSEsDpfrdympehaQgjMnz8ft912G2JjYwFczVABWM1Sz507Jy0rNTUVS5YsabnKEhEROUjRKFAcHLyxfv28vDwEBgaa5zujlQgA9Hq9+be4XmFhITw9PdGpUyenbKOluWVSNHfuXHzzzTfYv39/g5i1LNVWhrpo0SLMnz/f/Lq0tBRRUVHOqywREZGDNBonPBDWeHX9wMBAi6TIWeLj4/HRRx9ZzMvIyMDgwYPh5eXl9O21BLfrPnviiSfw4YcfIjMzE5GRkeb5er0eAKxmqde3Hl1Lp9OZT5CWOlGIiIjcTXl5OY4dO4Zjx44BuHqZyrFjx5CbmwvgaqPCAw88YF5+9uzZOHfuHObPn4/vv/8e//znP/HWW29hwYIFrVF9u7hNUiSEwNy5c7F582bs2bMH3bt3t4h3794der0eu3btMs8zGAzYu3cvhg0b5urqEhEROU1r3JKfnZ2NuLg4xMXFAQDmz5+PuLg4PPfccwCA/Px8c4IEXP0d3r59O7KysnDTTTfhz3/+M/7+97/j7rvvdt6BaGFu0332+OOPY9OmTfjggw8QEBBgbhEKCgqCj48PFEXBvHnzsGzZMvTq1Qu9evXCsmXL4Ovri6lTp7Zy7YmIiOznlMd8NHP9xMRE84XS1qSlpTWYl5CQgK+++qq5VVMNt0mKXn31VQBX36RrrVu3DjNmzAAALFy4EFVVVZgzZw6Ki4sxZMgQZGRkICAgwMW1JSIiInfjNkmRrWy1nqIoSElJQUpKSstXiIiIyEUUjQaKxsGWIgfXbw/cJikiIiJqrzQeTrj7rIVHtG4LeISIiIiIwJYiIiIi9XPGs8vYUtQoJkVEREQqp2iccPcZrylqFJMiIiIileOF1q7BI0REREQEthQRERGp3tXBGz0cLMPopNq0XUyK2qDCkgppzNPGU5ZNNoaCMhQXSGPajvom1et6V8or5bEa+YfXv875dbElatoD0phH7EhprO7CKWmsY4D8eXyKsU4eqymX16XwJ2nM5O0rjWkqiqQxYaiWr+djfVDUuqAI6TqV3sHSmAfk56YfauQxG48uyKmQx2x9FsL85F+NVXUmaayDjWdeBvjLG+YrjfK6lMNbGvPv1E0aKzLKKxNskp9jpqhYaawma6M09svO3dJY2Kjh0phuhPzz5Rk3Rhoz/fC5NGbMkY+oLHR+0pgS2tX6Oor1904IrbQsZ2uNEa3bIx4hIiIiIrCliIiISPU0Gg00Dl4o7ej67QGTIiIiIpVj95lr8AgRERERgS1FREREqseWItdgUkRERKRyiuKEwRsld9HR//AIEREREYEtRURERKrH7jPXYFJERESkckyKXINJERERkcppPDTQOJjUOLp+e8AjRERERAS2FBEREameolEcv/vMxvP+6ComRURERCrHa4pcg0eIiIiICGwpIiIiUj22FLkGkyIiIiKV44jWrsEjRERERAS2FBEREame4uEBjYeHw2WQbUyKiIiIVI7XFLkGjxARERER2FJERESkemwpcg0mRURERCqnaJxw95mD67cHTIqIiIhUji1FrsEjRERERAS2FBEREameolEcbyniA2EbxaSIiIhI5XhNkWvwCBERERGBLUVERESqp2g8oGgcHNHawfXbAyZFbZCPl7wBsNxgksY6edtoOKxxfqOij6laGuvi4yWNaUpK5IV21DtSJatMoT3kwbzj0pBn1ABpTPySI40pRoO8LjlfS2N1VRXyukT2lNdF5y+vi4f8fagLirA63+DbSbqOra9kTxvXO1QYtdJYpY1z+tyVSmnsxjA/acxoEtKY1kY9L1XL61JUbZTGbO27rctAKrzkX+E/Fsn3vWuQtzTmpQmUxqr6TpDGevb6lTRWEhgtjf2QL/88DwgPksbqguVlakMipTFn8iwtdcl2AAAaj6uTo2WQTew+IyIiIgJbioiIiNRPo7k6OVoG2cSkiIiISOUUDw+Hn3Lv6PrtAdNGIiIiIrCliIiISP14obVLMCkiIiJSO43GCUkRO4caw6SIiIhI5TiitWvwCBERERGBSREREZH6KR7/u67I3klpfvfbK6+8gu7du8Pb2xuDBg3CZ599Jl02KysLiqI0mE6ePOnInruUWyVF+/btw4QJExAREQFFUbB161aLuBACKSkpiIiIgI+PDxITE3HixInWqSwREZGzOJoQ2XGh9nvvvYd58+bh2WefxdGjR/HrX/8aY8eORW5urs31Tp06hfz8fPPUq1cvR/bcpdwqKaqoqMDAgQOxZs0aq/EVK1Zg1apVWLNmDQ4fPgy9Xo9Ro0ahrKzMxTUlIiJyb6tWrcLDDz+MmTNnom/fvli9ejWioqLw6quv2lwvNDQUer3ePHm40fhIbpUUjR07Fn/5y1+QnJzcICaEwOrVq/Hss88iOTkZsbGxWL9+PSorK7Fp06ZWqC0REZFz1F9o7egEAKWlpRZTTU1Ng+0ZDAYcOXIESUlJFvOTkpJw4MABm3WNi4tDeHg4Ro4ciczMTOcdBBdwq6TIlpycHBQUFFi8gTqdDgkJCTbfwJqamgYnCBERkao4sfssKioKQUFB5ik1NbXB5i5dugSj0YiwsDCL+WFhYSgoKLBaxfDwcLz++utIT0/H5s2bERMTg5EjR2Lfvn3OPx4tpM3ckl//Jll7A8+dOyddLzU1FUuWLGnRuhEREalFXl4eAgMDza91Op10WUVRLF4LIRrMqxcTE4OYmBjz6/j4eOTl5WHlypW4/fbbHay1a7SZlqJ6zXkDAWDRokUoKSkxT3l5eS1dRSIiouapH7zRoenqT35gYKDFZC0pCgkJgYeHR4NWocLCwgaND7YMHToUp0+fdmzfXajNJEV6vR4Amv0G6nS6BicIERGRmtQ/ENbRqam0Wi0GDRqEXbt2WczftWsXhg0b1uRyjh49ivDw8CYv39raTPdZ9+7dodfrsWvXLsTFxQG4eqHY3r178cILL7Ry7YiIiNzL/PnzMW3aNAwePBjx8fF4/fXXkZubi9mzZwO42tNy/vx5vP322wCA1atXo1u3bujfvz8MBgM2bNiA9PR0pKent+ZuNItbJUXl5eX48ccfza9zcnJw7NgxBAcHo2vXrpg3bx6WLVuGXr16oVevXli2bBl8fX0xderUVqw1ERGRgzQax59d1sz177vvPly+fBnPP/888vPzERsbi+3btyM6OhoAkJ+fbzFmkcFgwIIFC3D+/Hn4+Pigf//+2LZtG8aNG+dYvV3IrZKi7OxsDB8+3Px6/vz5AIDp06cjLS0NCxcuRFVVFebMmYPi4mIMGTIEGRkZCAgIaK0qExEROc6OwRetltFMc+bMwZw5c6zG0tLSLF4vXLgQCxcutKdmquFWSVFiYiKEENK4oihISUlBSkqK6ypFRETUwhSNBxQHkyJH128P2syF1kRERESOcKuWIiIionZJccI1RQrbQRrDpIiIiEjl2H3mGkwbiYiIiMCWIiIiIvWrH9Ha0TLIJiZFREREatcK4xS1RzxCRERERGBLUZukNRmksZCaIvmKtfKmWU11iTRmvHJBGhNaH2kszytCGutm+Fka8+zSVxqzpfaXHGnMo/yiNObV4xa7tmfM+Uoa09hqBq8olq/na+PZfP4dpCHF00saMxXlS2PG8ivyMv06WZ1fUmOSrhPiVSeNFVTJj8kX50ulsWAf+b59duayNHamuFIaq6o1SmOGOvn+FZXLP3s1NtazpZO/Vho7mS8/Lu//9VVpbPish6Wx2bffII0dt7G9zgHyz3qvThXSmM5D/rf54Vz5ZyHcxvkeKY24r+Y+u0xWBtnGpIiIiEjtWmlE6/aGSREREZHaMSlyCV5TRERERAS2FBEREameotFAcfDuMUfXbw+YFBEREamd4oTuM4XdZ41h2khEREQEthQRERGpn6I4/kBXRXFOXdowJkVERERqp2ickBSxc6gxPEJEREREYEsRERGR6glFA+FgS4+j67cHTIqIiIjUjt1nLsEjRERERAS2FBEREamfojh+9xjvPmsUkyIiIiK102iuTo6WQTYxKSIiIlI5XmjtGjxCRERERGBLERERkfrx7jOXYFJERESkdkyKXIJHiIiIiAhsKSIiIlI/thS5BJMiIiIilROK4oS7zzhOUWOYNhIRERGBLUVERETqx+4zl2BSREREpHZ8zIdLMG0kIiIiAluKiIiI1I/dZy7BpIiIiEjl+Owz12BSREREpHaKxvGn3DMpahSTIjdVW3hWGvOsvCKNKaY6acyk85PGjP6d5ZWx84MWXXPJrjKN576WxjyiB0pjXmHd5duzFbOTIkzSmOmKfN89ByZJY4bP/y1fT99Nvj3vAGkMfp3kMd+O0lC+0dvq/GP5ZdJ1jp0vkcY+Py0/Jr3C5PU/c7FcGsv9sUga8/CUn2PVFQZprKwwXxoL1EdIY3W1RnldPOR1MdTIP7NGQ5U01rnPUGnsu4OnpbEHN++QxgIjY6SxTl1CpDFPrYc0ZqiS71+fPvIy35txizRGZC8mRURERGrHa4pcgkkRERGR2jEpcgkeISIiIiKwpYiIiEj92FLkEkyKiIiIVI4PhHUNpo1ERETkdt5++23U1NQ0mG8wGPD222/bVSaTIiIiIrWr7z5zdGpDHnzwQZSUNBzmo6ysDA8++KBdZdp1hB566CGUlTUci6SiogIPPfSQXRUhIiIiifoHwjo6tSFCCChW9unnn39GUFCQXWXadU3R+vXrsXz5cgQEWA6oVlVVhbfffhv//Oc/7aoMERERkS1xcXFQFAWKomDkyJHw9PxfKmM0GpGTk4MxY8bYVXazWopKS0tRUlICIQTKyspQWlpqnoqLi7F9+3aEhobaVRFneuWVV9C9e3d4e3tj0KBB+Oyzz1q7SkRERPZrpe6z5v6e7t27F4MGDYK3tzduuOEGrF271t49lpo0aRImTpwIIQRGjx6NiRMnmqcpU6bgtddew4YNG+wqu1ktRR06dDBnZ717924QVxQFS5YssasizvLee+9h3rx5eOWVV3Drrbfitddew9ixY/Hdd9+ha9eurVo3IiIie7TGA2Gb+3uak5ODcePGYdasWdiwYQM+//xzzJkzB507d8bdd9/tUN2vtXjxYgBAt27dMGXKFOh0OqeV3aykKDMzE0IIjBgxAunp6QgODjbHtFotoqOjEREhf/aPK6xatQoPP/wwZs6cCQBYvXo1du7ciVdffRWpqamtWjciIiK7tMI4Rc39PV27di26du2K1atXAwD69u2L7OxsrFy50qlJUb0RI0bg4sWLiIyMBAB8+eWX2LRpE/r164dHHnnErjKblRQlJCQAuJoNRkVFQePoE3udzGAw4MiRI3jmmWcs5iclJeHAgQNW16mpqbG4pa+0tLRF60hERNSarv+d0+l0DVpb7Pk9PXjwIJKSLB9oPXr0aLz11luora2Fl5eXE2r/P1OnTsUjjzyCadOmoaCgAHfccQdiY2OxYcMGFBQU4Lnnnmt2mXZlNdHR0dBoNKisrMTJkyfxzTffWEyt5dKlSzAajQgLC7OYHxYWhoKCAqvrpKamIigoyDxFRUW5oqpERERNdnXwRscnAIiKirL43bPW6mPP72lBQYHV5evq6nDp0iUnHYn/+fbbb3HLLbcAAP79739jwIABOHDgADZt2oS0tDS7yrTr7rOLFy/iwQcfxI4dO6zGjUajXZVxlutv0ZPdtgcAixYtwvz5882vS0tLmRgREZGqCHF1crQMAMjLy0NgYKB5vq1rcprzeypb3tp8Z6itrTXXfffu3bjrrrsAAH369EF+fr5dZdrVUjRv3jwUFxfj0KFD8PHxwSeffIL169ejV69e+PDDD+2qiDOEhITAw8OjQRZbWFjYIHutp9PpEBgYaDERERG1Vdf/5llLiuz5PdXr9VaX9/T0RKdOnZy3A//Vv39/rF27Fp999hl27dplvg3/woULdm/PrqRoz549+Otf/4qbb74ZGo0G0dHR+N3vfocVK1a06sXMWq0WgwYNwq5duyzm79q1C8OGDWulWhERETnGJIRTpqay5/c0Pj6+wfIZGRkYPHiw068nAoAXXngBr732GhITE3H//fdj4MCBAIAPP/zQ3K3WXHZ1n1VUVJjHIwoODsbFixfRu3dvDBgwAF999ZVdFXGW+fPnY9q0aRg8eDDi4+Px+uuvIzc3F7Nnz27VehEREdlL/HdytIzmaOz3dNGiRTh//rz5OWOzZ8/GmjVrMH/+fMyaNQsHDx7EW2+9hX/9618O1ty6IUOG4Oeff4bBYEDHjh1x7tw5bNmyBdHR0XjqqafsKtOupCgmJganTp1Ct27dcNNNN+G1115Dt27dsHbtWoSHh9tVEWe57777cPnyZTz//PPIz89HbGwstm/fjujo6FatFxERkTtp7Pc0Pz8fubm55uW7d++O7du346mnnsLLL7+MiIgI/P3vf2+R2/EBYOLEiUhOTsbs2bNx5coV3HLLLdBqtbh06RJMJhMee+yxZpdpV1I0b94880VMixcvxujRo7FhwwZotVqsX7/eniKdas6cOZgzZ05rV4OIiMgpTOLq5GgZzWXr99TaHV4JCQku6zH66quv8Ne//hUA8P7770Ov1+Po0aNIT0/Hc88957qk6Le//a35/3FxcTh79ixOnjyJrl27IiQkxJ4iiYiISEIIYb6Ty5Ey2pLKykrzM1gzMjKQnJwMjUaDoUOH4ty5c3aV2eSk6Nrb1huzatUquypDRERE1BQ9e/bE1q1bMXnyZOzcudN8HVFhYaHdd5I3OSk6evRok5ZribEIiIiI2rPW6j5Ts+eeew5Tp07FU089hZEjRyI+Ph7A1VajuLg4u8psclKUmZlp1waIiIjIcW0sp3HYPffcg9tuuw35+fnm2/EBYOTIkZg8ebJdZdp1TRERERG5DluKrNPr9dDr9Rbz7B2jCLBz8EYiIiKitoYtRURERCrHu89cg0mRitWd/14e9PSWhoTWRx5T5I2DJu8A+Xo2tqeY6qQxjyvn5WV6ycs0+XaUxjyLcqWxuvzT8vXCe0ljtth6H4TWT76iv3x4irpvD0hjhh+/kcZ8714g354NZW/9SRrTRXWXxjwGjpDGOgf5W53/K731+QBQWF4jjT14azdpzEsjv4EjvIP8PIrtEiSNVRrkD64+d7lCGhvao580tu/URWmsg6/8MQc3dJYfsz1H5J+hDiHywXLLS+XHWust/+oP8u8jjV28VCmNGetM0pjexjlx6XKVNLb3o4PSGGbY30Xijkz/nRwtg2xj9xkRERER2FJERESkekJcnRwtg2xjUkRERKRyvPvMNdh9RkRERAS2FBEREake7z5zDSZFREREKse7z1yD3WdEREREYEsRERGR6gk44e4zp9SkbWNSREREpHImIWByMCtydP32gEkRERGRygk43tLDlKhxvKaIiIiICGwpIiIiUj0O3ugaTIqIiIjUzgmP+WD/WePYfUZEREQEthQRERGpngkCJgebehxdvz1gUkRERKRywgndZ7wjv3HsPiMiIiICW4qIiIhUj3efuQaTIiIiIpVj95lrsPuMiIiICGwpIiIiUj3efeYaTIqIiIhUjt1nrsGkiIiISOVMQjj8lHtH128PeE0REREREdhSREREpHpG09XJ0TLINiZFREREKsfuM9dg9xkRERER2FJERESkeiYhYGRLUYtjUqRiJv/O8pjWTxqrsdFv7KHIY5W18hU9IF/RR+stjSm+HeUxQ4U0Jjy00pjRVpl11dKY6Uy2fHuKjUZTrY98ezVlNupSKy+y543SWF3heWms9out8ljuD9JY0MN/kcYu/OUxaaxjtfw90vUZZHV+aGgv6ToPxMjPW1vvQYnQSWN9QvylsTJDnTRWWGGQxmrtvPgiLEBez9IaeV2KyuV1+VW/ULvqUuTrJY0F+co/X73C5MfTP1b+k3GxrEYaiwiSf0doNPLvlo0+8u3pBs2SxmqOvCGNuaurj/lwNClyUmXaMHafEREREYEtRURERKrHu89cg0kRERGRyvHuM9dg9xkRERER2FJERESkekYn3H3m6PrtAZMiIiIilTPB8bvHeElR45gUERERqZzRJGB0MCtydP32gNcUEREREcGNkqKlS5di2LBh8PX1RYcOHawuk5ubiwkTJsDPzw8hISF48sknYTDIB0UjIiJyB+K/d585MgleU9Qot0mKDAYD7r33Xjz2mPXRd41GI+68805UVFRg//79ePfdd5Geno7f//73Lq4pERGRcxmFc6aWUlxcjGnTpiEoKAhBQUGYNm0arly5YnOdGTNmQFEUi2no0KEtV8kmcJtripYsWQIASEtLsxrPyMjAd999h7y8PERERAAAXnrpJcyYMQNLly5FYGCgq6pKRETUrkydOhU///wzPvnkEwDAI488gmnTpuGjjz6yud6YMWOwbt0682utVv4IGldwm6SoMQcPHkRsbKw5IQKA0aNHo6amBkeOHMHw4cNbsXZERET2U/Pgjd9//z0++eQTHDp0CEOGDAEAvPHGG4iPj8epU6cQExMjXVen00Gv17dIvezhNt1njSkoKEBYWJjFvI4dO0Kr1aKgoEC6Xk1NDUpLSy0mIiIiNam/+8zRqSUcPHgQQUFB5oQIAIYOHYqgoCAcOHDA5rpZWVkIDQ1F7969MWvWLBQWFrZIHZuqVZOilJSUBv2J10/Z2fInm19PURo+bVkIYXV+vdTUVHMfaFBQEKKiouzaFyIiIndwfUNATU2NQ+UVFBQgNDS0wfzQ0FCbjRJjx47Fxo0bsWfPHrz00ks4fPgwRowY4XB9HNGq3Wdz587FlClTbC7TrVu3JpWl1+vxxRdfWMwrLi5GbW1tgxakay1atAjz5883vy4tLWViREREquLM7rPrf+MWL16MlJSUBsunpKSYr+eVOXz4MAD7GiXuu+8+8/9jY2MxePBgREdHY9u2bUhOTra53ZbSqklRSEgIQkJCnFJWfHw8li5divz8fISHhwO4evG1TqfDoEGDpOvpdDrodDqn1IGIiKglOOPusfr18/LyLG4+kv0GNrXh4ptvvsEvv/zSIHbx4kWbjRLXCw8PR3R0NE6fPt3kdZzNbS60zs3NRVFREXJzc2E0GnHs2DEAQM+ePeHv74+kpCT069cP06ZNw4svvoiioiIsWLAAs2bN4p1nRERE/xUYGNik38WmNlzEx8ejpKQEX375JW655RYAwBdffIGSkhIMGzasyfW6fPky8vLyzA0brcFtLrR+7rnnEBcXh8WLF6O8vBxxcXGIi4szX3Pk4eGBbdu2wdvbG7feeit+85vfYNKkSVi5cmUr15yIiMgxjg7c6IzuN5m+fftizJgxmDVrFg4dOoRDhw5h1qxZGD9+vMWdZ3369MGWLVsAAOXl5ViwYAEOHjyIs2fPIisrCxMmTEBISAgmT57cIvVsCrdpKUpLS5OOUVSva9eu+Pjjj11TISIiIhcxmQRMDt495uj6tmzcuBFPPvkkkpKSAAB33XUX1qxZY7HMqVOnUFJSAuBqQ8bx48fx9ttv48qVKwgPD8fw4cPx3nvvISAgoMXq2Ri3SYqIiIjaK5MTrilqyefBBgcHY8OGDTaXufYxIz4+Pti5c2fLVchObtN9RkRERNSS2FJERESkcmoe0botYVJERESkckYhYHQwqXF0/faA3WdEREREYEsRERGR6qn97rO2gkkRERGRyhnhhBGtnVKTto3dZ0RERERgS5GqaYPsey6cplj+VOIyT/nQ7oFaeY58oaJOGqu10STr7SPfBz9hksY8i85JY8LLWx7T+kljdf6dpbHaHa9LY6Vn86Uxn05B0liHR5ZJY7Z49LNrNXgNkceO3z9OGqsqrpbGqi+XSGMhly9Zna+L+kG6jle3vtKYqK6Qxjr5d5LGgoK7SmMXdVpprKOP/OvP38tDGquslZ+331+S74NGI38wZveOvtJYTnGlNHYyv0waqzTI2wXqTPKnkBvq5PuXZ2P/FBv7Fxks37/be8u/I8YMtPG4h+kPyGNtEO8+cw0mRURERCrHu89cg91nRERERGBLERERkeqZTAJG3n3W4pgUERERqZzRCUmRo+u3B0yKiIiIVI5JkWvwmiIiIiIisKWIiIhI9Ywmx1t6jPLRFui/mBQRERGpHLvPXIPdZ0RERERgSxEREZHqsaXINZgUERERqRzHKXINdp8RERERgS1FREREqmcUTug+47PPGsWkiIiISOV4TZFrsPuMiIiICGwpIiIiUj22FLkGkyIiIiKVqzMJeDiY1NQxKWoUkyIiIiKVY0uRa/CaIiIiIiKwpYiIiEj1OHijazApIiIiUjmjEA6PM8RxihrH7jMiIiIisKWIiIhI9XihtWswKSIiIlI5JkWuwe4zIiIiIrCliIiISPXYUuQaTIqIiIhUzihMMJpMDpdBtjEpamUXSyulsc6BvnaVqdRWS2P+Wj9pzONyrjQWHtJDGiurMUpjtTY+gx5XzktjJt8O8hVtMcnrYotnWFdpLPz+P9lXFxUZ8K/t0tgWfX+7yvT01lmdbzyRI11HG3BIGjNWG6SxoJ5dpDFhlJ9kYXFDpTFF4yGP+cg/J3WRN0pj0V3kX6nna+Vlaj3kVzL07CT/HtD7W38PACCnSP7dYqvF4HK5/H3w7xIkjeWXVEljhaXy76QtX8m/B3y18vdoRP8waWzC6welsY8eiZfGiJgUERERqRwHb3QNJkVEREQqZzQJaHhNUYtjUkRERKRydSZAcTCpqeMlRY3iLflEREREYEsRERGR6rH7zDWYFBEREakckyLXYPcZEREREdhSREREpHpsKXINJkVEREQqx3GKXMMtus/Onj2Lhx9+GN27d4ePjw969OiBxYsXw2CwHHk1NzcXEyZMgJ+fH0JCQvDkk082WIaIiIjIGrdoKTp58iRMJhNee+019OzZE99++y1mzZqFiooKrFy5EgBgNBpx5513onPnzti/fz8uX76M6dOnQwiBf/zjH628B0RERPYzmoTD4xSx+6xxbtFSNGbMGKxbtw5JSUm44YYbcNddd2HBggXYvHmzeZmMjAx899132LBhA+Li4nDHHXfgpZdewhtvvIHS0tJWrD0REZFjhBAQJgcn0XJJ0dKlSzFs2DD4+vqiQ4cOTd6nlJQUREREwMfHB4mJiThx4kSL1bEp3CIpsqakpATBwcHm1wcPHkRsbCwiIiLM80aPHo2amhocOXKkNapIRETULhgMBtx777147LHHmrzOihUrsGrVKqxZswaHDx+GXq/HqFGjUFZW1oI1tc0tus+u99NPP+Ef//gHXnrpJfO8goIChIVZPjW5Y8eO0Gq1KCgokJZVU1ODmpoa82u2KhERkdqYTMLhC6Vb8kLrJUuWAADS0tKatLwQAqtXr8azzz6L5ORkAMD69esRFhaGTZs24dFHH22pqtrUqi1FKSkpUBTF5pSdnW2xzoULFzBmzBjce++9mDlzpkVMUZQG2xBCWJ1fLzU1FUFBQeYpKirKOTtHRETkJEIIp0zA1T/+r52ubRhwlZycHBQUFCApKck8T6fTISEhAQcOHHB5feq1akvR3LlzMWXKFJvLdOvWzfz/CxcuYPjw4YiPj8frr79usZxer8cXX3xhMa+4uBi1tbUNWpCutWjRIsyfP9/8urS0lIkRERGpSv11QY6WAaDBb9zixYuRkpLiUNnNVd+Dc/3vc1hYGM6dO+fSulyrVZOikJAQhISENGnZ8+fPY/jw4Rg0aBDWrVsHjcaykSs+Ph5Lly5Ffn4+wsPDAVy9+Fqn02HQoEHScnU6HXQ6nf07QURE5Eby8vIQGBhofi37DUxJSTF3i8kcPnwYgwcPtrsu1/fkNNa709Lc4pqiCxcuIDExEV27dsXKlStx8eJFc0yv1wMAkpKS0K9fP0ybNg0vvvgiioqKsGDBAsyaNcvizSciInI3zrymKDAwsEm/i83tzWmO+t/ugoICc0MGABQWFtrs3WlpbpEUZWRk4Mcff8SPP/6IyMhIi1h9H6mHhwe2bduGOXPm4NZbb4WPjw+mTp1qHseIiIjIXQnT1cnRMpqjOb05zdW9e3fo9Xrs2rULcXFxAK7ewbZ371688MILLbLNpnCLpGjGjBmYMWNGo8t17doVH3/8cctXiIiIiMxyc3NRVFSE3NxcGI1GHDt2DADQs2dP+Pv7AwD69OmD1NRUTJ48GYqiYN68eVi2bBl69eqFXr16YdmyZfD19cXUqVNbbT/cIikiIiJqz669e8yRMlrKc889h/Xr15tf17f+ZGZmIjExEQBw6tQplJSUmJdZuHAhqqqqMGfOHBQXF2PIkCHIyMhAQEBAi9WzMUyKiIiIVE7t4xSlpaU1OkbR9UmZoihISUlx+Z1vtrjtiNZEREREzsSWIiIiIpVz5jhFJMekiIiISO2ckBSBSVGj2H1GREREBLYUERERqZ5JCCgO3j1masG7z9oKJkVEREQqJ4QTriliUtQoJkWtzMfTvme8GPOOS2Om00fl68X/RhrTnD8tjVVuTZPGOkx7WhoTivwUM0T9ShrTnv9aGjMG6qUxj/JL8vU8vKSxK4PulsZCpZG24dYFd9i1XlneL1bn++o7Sde59E2ONOan7yCNeXhrpTGNp/x9FYZqaazkyBfSWPXlUmkssNseaaz8/EVpTGuUDyesDfSVxvrcNlwei+ovjZl6yEcivmSUP++xqk5eT1u/qWeKq2zEKqWxkqpaaayo3CCN7T5eII35+8rPCXfFC61dg9cUEREREYEtRURERKpnMgGKw4M3OqkybRiTIiIiIpVT+2M+2gp2nxERERGBLUVERESqJ0xXJ0fLINuYFBEREamcySSccE0Ru88aw+4zIiIiIrCliIiISPU4TpFrMCkiIiJSOSZFrsHuMyIiIiKwpYiIiEj1+EBY12BSREREpHLsPnMNJkVEREQqJ4QTkiK2FDWK1xQRERERgS1FREREqidMwuHBF9l91jgmRURERCrHB8K6BrvPiIiIiMCWIiIiItXj3WeuwaSIiIhI5UwmAfCBsC2O3WdEREREYEsRERGR6gmTEcJkdLgMso1JERERkcoxKXINdp8RERERgS1FREREqidMJie0FJmcVJu2i0kRERGRygmjEcLoYFLk4PrtAZMiIiIilRPCCdcUCSZFjeE1RURERERgS5HL1JRdsTpfW1ctXSf9p2JpzCg6SWMTSi9LYx4H/y2NCZ23NKYN9JPGPC/+KI2ZKkqlMY2PvEx46qQhj/JL0ljt2e+kMWVQN2nMX+shr0sb5+XnY9d6BUfOW53fa1KYdB3/LvLztiK/SBr75Uv5+6rrECCNaU7nSmM+nYKkMe9OgdJYyU/W9xsANF7yr9TqohJprPBojjR2+YQ8Fh4fK43pwrtIY/r+w6Qxpa5GGhM2PpddI2+Qxn7dVX48DUb5dS6lBnlsR6i/NPbx1xekMXfFu89cg0kRERGRyjEpcg12nxERERGBLUVERESqx5Yi12BSREREpHIcp8g12H1GREREBLYUERERqZ7JZAQcbCkysfusUUyKiIiIVI7XFLkGu8+IiIiIwJYiIiIi1WNLkWu4TUvRXXfdha5du8Lb2xvh4eGYNm0aLlywHLU0NzcXEyZMgJ+fH0JCQvDkk0/CYDC0Uo2JiIic5L8PhHVkAh8I2yi3SYqGDx+Of//73zh16hTS09Px008/4Z577jHHjUYj7rzzTlRUVGD//v149913kZ6ejt///vetWGsiIiLH1T8Q1qGJD4RtlNt0nz311FPm/0dHR+OZZ57BpEmTUFtbCy8vL2RkZOC7775DXl4eIiIiAAAvvfQSZsyYgaVLlyIwUP7sHSIiIiK3aSm6VlFRETZu3Ihhw4bBy8sLAHDw4EHExsaaEyIAGD16NGpqanDkyJHWqioREZHD6gdvdGxqucEbly5dimHDhsHX1xcdOnRo0jozZsyAoigW09ChQ1usjk3hVknR008/DT8/P3Tq1Am5ubn44IMPzLGCggKEhVk+nbtjx47QarUoKCiQlllTU4PS0lKLiYiISE0cT4gcv1DbFoPBgHvvvRePPfZYs9YbM2YM8vPzzdP27dtbqIZN06pJUUpKSoMs8fopOzvbvPwf/vAHHD16FBkZGfDw8MADDzwAIYQ5rihKg20IIazOr5eamoqgoCDzFBUV5dydJCIiauOWLFmCp556CgMGDGjWejqdDnq93jwFBwe3UA2bplWvKZo7dy6mTJlic5lu3bqZ/x8SEoKQkBD07t0bffv2RVRUFA4dOoT4+Hjo9Xp88cUXFusWFxejtra2QQvStRYtWoT58+ebX5eWljIxIiIiVREmE+Bg91d999n1PSI6nQ46nc6hsu2VlZWF0NBQdOjQAQkJCVi6dClCQ0NbpS5AKydF9UmOPepbiGpqagAA8fHxWLp0KfLz8xEeHg4AyMjIgE6nw6BBg6TltObJQERE1BTCCY/5qO8+u/4P/8WLFyMlJcWhsu0xduxY3HvvvYiOjkZOTg7+3//7fxgxYgSOHDnSar/LbnH32Zdffokvv/wSt912Gzp27IgzZ87gueeeQ48ePRAfHw8ASEpKQr9+/TBt2jS8+OKLKCoqwoIFCzBr1izeeUZERPRfeXl5Fr+LsgQkJSUFS5YssVnW4cOHMXjwYLvqcd9995n/Hxsbi8GDByM6Ohrbtm1DcnKyXWU6yi2SIh8fH2zevBmLFy9GRUUFwsPDMWbMGLz77rvmN9PDwwPbtm3DnDlzcOutt8LHxwdTp07FypUrW7n2REREjnFmS1FgYGCTGguae4mLo8LDwxEdHY3Tp087rczmcoukaMCAAdizZ0+jy3Xt2hUff/yxC2pERETkOiaTEYqLH/PhyCUu9rh8+TLy8vLMl8C0Bre6JZ+IiIjUJzc3F8eOHUNubi6MRiOOHTuGY8eOoby83LxMnz59sGXLFgBAeXk5FixYgIMHD+Ls2bPIysrChAkTEBISgsmTJ7fWbrhHSxEREVF7JowmQHGwpcjYcoM3Pvfcc1i/fr35dVxcHAAgMzMTiYmJAIBTp06hpKQEwNVLXo4fP463334bV65cQXh4OIYPH4733nsPAQEBLVbPxjApIiIiUjkhnHBNUQs++ywtLQ1paWmNbP9/4wr6+Phg586dLVYfezEpIiIiUjlhMjreUtSCI1q3FbymiIiIiAhsKSIiIlI9thS5BpMiIiIilWNS5BpMiq5TfyHY9c+GcVRNmfXylLoa6TqV5fI7BYxCGkJpZbU0pvGoksaUOvn2aqrk9TSWVUhjpkp5TGPr8+lRayMoV1thY/9Ky6SxOq183+tqDXbVxV2U2nhvbamoq7M6v6xaXl5FjfxYVtbK3/O6GnnMYKNMjY27bWpt7LepTr69Kju3V21jH8prrR9LADAa5OvZeu90lfLPgldZuTRm6ztJeNp4j7Ty70yjxksaM9g4ZmUGeayqQv55rq2Sf+8487u9vqxrLyBuMcZaOLwVo33fq+2JIlzybrqPn3/+mQ+EJSKiJsvLy0NkZGSLlF1dXY3u3bujoKDAKeXp9Xrk5OTA29vbKeW1NUyKrmMymXDhwgUEBARAUZTWrk6TlJaWIioqqsHzbNo7HpeGeEys43FpiMfEumuPS0BAAMrKyhAREQGNpuXuW6qurobB4JwWa61Wy4TIBnafXUej0bRYxt/Smvo8m/aGx6UhHhPreFwa4jGxrv64BAUFtfi2vL29mci4CG/JJyIiIgKTIiIiIiIATIraBJ1Oh8WLF0On07V2VVSFx6UhHhPreFwa4jGxjselbeOF1kRERERgSxERERERACZFRERERACYFBEREREBYFJEREREBIBJkdtISUmBoigWk16vN8eFEEhJSUFERAR8fHyQmJiIEydOtGKNW8a+ffswYcIEREREQFEUbN261SLelONQU1ODJ554AiEhIfDz88Ndd92Fn3/+2YV74VyNHZMZM2Y0OHeGDh1qsUxbOyapqam4+eabERAQgNDQUEyaNAmnTp2yWKY9nitNOS7t8Xx59dVXceONN5oHZIyPj8eOHTvM8fZ4rrRXTIrcSP/+/ZGfn2+ejh8/bo6tWLECq1atwpo1a3D48GHo9XqMGjUKZWXyhya6o4qKCgwcOBBr1qyxGm/KcZg3bx62bNmCd999F/v370d5eTnGjx8Po9E9nyDd2DEBgDFjxlicO9u3b7eIt7VjsnfvXjz++OM4dOgQdu3ahbq6OiQlJaGi4n8PCm2P50pTjgvQ/s6XyMhILF++HNnZ2cjOzsaIESMwceJEc+LTHs+VdkuQW1i8eLEYOHCg1ZjJZBJ6vV4sX77cPK+6uloEBQWJtWvXuqiGrgdAbNmyxfy6KcfhypUrwsvLS7z77rvmZc6fPy80Go345JNPXFb3lnL9MRFCiOnTp4uJEydK12nrx0QIIQoLCwUAsXfvXiEEz5V61x8XIXi+1OvYsaN48803ea60M2wpciOnT59GREQEunfvjilTpuDMmTMAgJycHBQUFCApKcm8rE6nQ0JCAg4cONBa1XW5phyHI0eOoLa21mKZiIgIxMbGtuljlZWVhdDQUPTu3RuzZs1CYWGhOdYejklJSQkAIDg4GADPlXrXH5d67fl8MRqNePfdd1FRUYH4+HieK+0MkyI3MWTIELz99tvYuXMn3njjDRQUFGDYsGG4fPkyCgoKAABhYWEW64SFhZlj7UFTjkNBQQG0Wi06duwoXaatGTt2LDZu3Ig9e/bgpZdewuHDhzFixAjU1NQAaPvHRAiB+fPn47bbbkNsbCwAniuA9eMCtN/z5fjx4/D394dOp8Ps2bOxZcsW9OvXj+dKO+PZ2hWgphk7dqz5/wMGDEB8fDx69OiB9evXmy+CVBTFYh0hRIN57YE9x6EtH6v77rvP/P/Y2FgMHjwY0dHR2LZtG5KTk6XrtZVjMnfuXHzzzTfYv39/g1h7Pldkx6W9ni8xMTE4duwYrly5gvT0dEyfPh179+41x9vzudKesKXITfn5+WHAgAE4ffq0+S606/8iKSwsbPDXTVvWlOOg1+thMBhQXFwsXaatCw8PR3R0NE6fPg2gbR+TJ554Ah9++CEyMzMRGRlpnt/ezxXZcbGmvZwvWq0WPXv2xODBg5GamoqBAwfib3/7W7s/V9obJkVuqqamBt9//z3Cw8PRvXt36PV67Nq1yxw3GAzYu3cvhg0b1oq1dK2mHIdBgwbBy8vLYpn8/Hx8++237eZYXb58GXl5eQgPDwfQNo+JEAJz587F5s2bsWfPHnTv3t0i3l7PlcaOizXt4XyxRgiBmpqadnuutFutcXU3Nd/vf/97kZWVJc6cOSMOHTokxo8fLwICAsTZs2eFEEIsX75cBAUFic2bN4vjx4+L+++/X4SHh4vS0tJWrrlzlZWViaNHj4qjR48KAGLVqlXi6NGj4ty5c0KIph2H2bNni8jISLF7927x1VdfiREjRoiBAweKurq61toth9g6JmVlZeL3v/+9OHDggMjJyRGZmZkiPj5edOnSpU0fk8cee0wEBQWJrKwskZ+fb54qKyvNy7THc6Wx49Jez5dFixaJffv2iZycHPHNN9+IP/7xj0Kj0YiMjAwhRPs8V9orJkVu4r777hPh4eHCy8tLREREiOTkZHHixAlz3GQyicWLFwu9Xi90Op24/fbbxfHjx1uxxi0jMzNTAGgwTZ8+XQjRtONQVVUl5s6dK4KDg4WPj48YP368yM3NbYW9cQ5bx6SyslIkJSWJzp07Cy8vL9G1a1cxffr0Bvvb1o6JteMBQKxbt868THs8Vxo7Lu31fHnooYdEdHS00Gq1onPnzmLkyJHmhEiI9nmutFeKEEK4rl2KiIiISJ14TRERERERmBQRERERAWBSRERERASASRERERERACZFRERERACYFBEREREBYFJEREREBIBJEREBSExMxLx581q7GkRErYpJERERERGYFBEREREBYFJERNcpLi7GAw88gI4dO8LX1xdjx47F6dOnzfG0tDR06NABO3fuRN++feHv748xY8YgPz+/FWtNROQ4JkVEZGHGjBnIzs7Ghx9+iIMHD0IIgXHjxqG2tta8TGVlJVauXIl33nkH+/btQ25uLhYsWNCKtSYicpxna1eAiNTj9OnT+PDDD/H5559j2LBhAICNGzciKioKW7duxb333gsAqK2txdq1a9GjRw8AwNy5c/H888+3Wr2JiJyBLUVEZPb999/D09MTQ4YMMc/r1KkTYmJi8P3335vn+fr6mhMiAAgPD0dhYaFL60pE5GxMiojITAghna8oivm1l5eXRVxRFOm6RETugkkREZn169cPdXV1+OKLL8zzLl++jB9++AF9+/ZtxZoREbU8JkVEZNarVy9MnDgRs2bNwv79+/H111/jd7/7Hbp06YKJEye2dvWIiFoUkyIisrBu3ToMGjQI48ePR3x8PIQQ2L59e4MuMyKitkYRvBCAiIiIiC1FRERERACTIiIiIiIATIqIiIiIADApIiIiIgLApIiIiIgIAJMiIiIiIgBMioiIiIgAMCkiIiIiAsCkiIiIiAgAkyIiIiIiAEyKiIiIiAAwKSIiIiICAPx/qB0jhcTX9bUAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ssta0[-1].plot()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkUAAAHFCAYAAAD8Jo2EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABPPUlEQVR4nO3de1xUdf4/8NeZAYb7TYSBuEiJeMFcV0uxNlATpTIv3f2uaRc316yvmttm7SOx708xM9d2raytFU3NdkOtXS2xBM1VS0y/uZZmiYHKhBfkDgMzn98fLvN1ZD4DzAwzB3g9H4/zUM77XD7z4cB58/l8zucoQggBIiIiom5O4+kCEBEREakBkyIiIiIiMCkiIiIiAsCkiIiIiAgAkyIiIiIiAEyKiIiIiAAwKSIiIiICwKSIiIiICACTIiIiIiIATIrIRfbt24esrCxcvny5RSw9PR3p6eluL5M7fP/995g/fz6GDBmC0NBQhIeH45ZbbsGHH35oc/uysjJMnz4dERER8Pf3R2pqKj7//HOrbSorK7F48WKkp6dDr9cjMDAQAwcOxMsvv4z6+voWx/zDH/6Au+66C9dddx0URcH06dPb9RnOnDmDOXPmIC0tDaGhoVAUBTk5OS22a2+5ZFauXInJkycjMTERiqK0em189NFHSEtLQ3BwMAICAjBgwAC8/fbbAICsrCwoitLq0nyOY8eOYdasWUhNTUVAQAAURUFBQUGLc5aWluIPf/gDUlNTERERgeDgYAwZMgRvv/02TCZTmz/rqVOnMHnyZISGhiIwMBBjxozB119/bXPbTZs24Re/+AV8fX0RExODOXPmoLq6WpXnIuqyBJELvPLKKwKAKCoqahE7duyYOHbsmPsL5QZ//vOfRd++fcXixYtFXl6e2L59u5g2bZoAIBYtWmS1bX19vUhJSRGxsbFi/fr1Ii8vT0yYMEF4eXmJgoICy3ZHjx4VERERYu7cueKjjz4Sn3/+ucjKyhK+vr5i9OjRwmw2Wx3X399fDB8+XMycOVP4+PiIadOmtesz5Ofni4iICHH77beLhx56SAAQa9asabFde8slk5ycLH75y1+KRx99VPTs2VOkpaVJt83OzhYajUbMmjVLfPLJJ+Kzzz4Tq1atEn/+85+FEEKUlJSI/fv3W5bNmzcLAOKpp56yWt98/eXk5Ijo6Ghxxx13iPHjxwsAIj8/v8V5//GPf4i4uDjxwgsviG3btom8vDwxd+5codFoxCOPPNKmz1lWViZiYmLEgAEDRG5urti2bZu49dZbRVBQkDh+/LjVtuvXrxcAxOOPPy527dolVq9eLUJCQsSYMWNUdy6iroxJEbmEvaSoKzt//rzNZODOO+8U/v7+or6+3rLu9ddfFwDEvn37LOsaGxtF//79xc0332xZV11dLaqrq1scs7mOv/jiC6v1JpPJ8v+AgIB2J0VX73/w4EFpUtTecrXlfAMGDJAmRYWFhUKj0YiXX365TccVQoiioiIBQLzyyiutnvvvf/+7NCm6dOmSMBqNLdY/+eSTAoAoLi5utSy/+93vhLe3tzh9+rRlXUVFhYiIiBD333+/ZV1TU5OIjo4WGRkZVvtv2LBBABDbt29X1bmIujJ2n5HTsrKy8Lvf/Q4ALF0iV3dLXNt9dvr0aSiKgldeeQUvv/wyevXqBT8/P6Snp+P7779HY2MjnnvuOcTExCAkJASTJk1CWVlZi/N+8MEHlm6QwMBAjB07FocPH3bHR7aIiIiAoigt1t98882ora3FpUuXLOu2bNmC5ORkpKamWtZ5eXnh17/+Nb766iucPXsWABAQEICAgACbxwSAkpISq/UajXM/xm3dv73lcvZ8q1atgk6nw1NPPdWm7V157rCwMHh7e7dY3/xZz5w50+oxtmzZglGjRiEhIcGyLjg4GJMnT8Y//vEPNDU1AQAOHDiA0tJSPPLII1b733fffQgMDMSWLVtUdS6iroxJETnt8ccft9y4Nm/ejP3792P//v345S9/aXe/119/Hf/617/w+uuv45133sHx48cxfvx4PPbYYzh//jz++te/YtmyZfjss8/w+OOPW+27ZMkSPPTQQ+jfvz/+9re/4b333kNVVRV+9atf4dtvv221zE1NTW1ahBAO1Ul+fj569uyJyMhIy7p///vfuPHGG1ts27zu2LFjdo+5a9cuAMCAAQMcKlNH6ahy7dmzB/369UNubi6Sk5Oh1WoRGxuL5557Dkaj0aXnaqtdu3bBy8sLffr0sVp/7diouro6/Pjjj9Lvd11dHU6dOgXgynXRvP5q3t7e6Nu3ryXuiXMRdTdeni4AdX6xsbGIj48HAAwePBi9evVq036hoaHYunWr5a/3CxcuYM6cOejbty8++ugjy3bHjx/HypUrUVlZieDgYJSUlGDhwoWYPXs2/vSnP1m2GzNmDJKSkrBo0SJ88MEH0vOePn0aiYmJbSpjfn5+uweJv/POOygoKMBrr70GrVZrWX/x4kWEh4e32L553cWLF6XH/Oabb7Bs2TJMmjTJ5s3PUzqyXGfPnsX58+fx9NNP43/+53/Qv39/fP7551i6dClKSkqwYcMGl56vNXl5eXjvvffw3//93+jRo4dVTKvVWn2vy8vLIYRo0/e7+V/ZtqdPn/bYuYi6GyZF5DF33HGHVXdGv379AAB33nmn1XbN64uLi5GSkoIdO3agqakJDz/8sKVbAAB8fX2RlpaG/Px8u+eNiYnBwYMH21TG5OTkNm3X7JNPPsGTTz6Je++912a3j62uttZip0+fxl133YW4uDi888477SpPMyFEi6emvLyc+/G3V66rvy/AlRu5vc9ui9lsRlVVFd5//308+OCDAICRI0eipqYGK1euxKJFi9C7d2+nPkNbff3117j//vsxfPhwZGdnt4hf+3mbtef7Ldv22vXuPBdRd8OkiDzm2r9WfXx87K5vfuz7559/BgDcdNNNNo/b2rgRHx8f/OIXv2hTGa/+i7w1O3bswOTJkzFmzBhs2LChxQ2mR48eNluDmscd2frr/aeffsLIkSPh5eWFzz//3OY2bbF27doW40gc7RpsS7muHY+zZs2adk8V0KNHDxgMBowdO9ZqfWZmJlauXImvv/7aLUnR4cOHLa2Q27dvh06na3WfsLAwKIrSpu93c6vTxYsXERUV1WLb1r7n7jwXUVfHpIg6nYiICADAhx9+aDWwtK06ovtsx44dmDhxItLS0pCbm2tJ5K42cOBAHD16tMX65nUpKSlW63/66Sekp6dDCIGCggLExsa2qcy2jB8/vs2tY61pS7muPVdb6/tqN954IwwGQ4v1zcmcswPM2+Lw4cO4/fbbkZCQgLy8PISEhLRpPz8/P/Tu3Vv6/fbz88P1118P4Mp10by+f//+lu2amppw/PhxPPTQQ6o5F1FXx6SIXKL5r+e6uroOP9fYsWPh5eWFH3/8Effcc0+793d191leXh4mTpyIW2+9FVu3bpW2JEyaNAmzZs3Cl19+iWHDhgG4cjNav349hg0bhpiYGMu2xcXFSE9Ph8lkQkFBgUPJ39V69OjRYhyMI9parqFDhzp9rnvuuQd5eXn45JNPMGXKFMv67du3Q6PRSFsKXeXIkSO4/fbbERsbi507dyIsLKxd+0+aNAkrV65ESUkJ4uLiAABVVVXYvHkz7r77bkv35bBhwxAdHY2cnBw88MADlv0//PBDVFdXY/Lkyao6F1FXxqSIXKL5L9DXXnsN06ZNg7e3N5KTkxEUFOTyc/Xq1QsvvfQSXnjhBZw6dQrjxo1DWFgYfv75Z3z11VcICAjAokWLpPv7+Pi45KYNAHv37sXEiROh1+vx/PPP48iRI1bx/v37Izg4GADw6KOP4vXXX8d9992HpUuXIjIyEm+88QZOnDiBzz77zLJPWVkZRo4cidLSUrz77rsoKyuzmpIgNjbWqnVm9+7dOH/+PADAZDLhp59+ssyonZaWhp49e7b6OZq3b35KqbCwEIGBgQCAe++916FyyRQWFloG9FZWVkIIYTn/TTfdZEm0HnnkEbz11luYNWsWLly4gP79++Ozzz7D66+/jlmzZjmUKNbW1mL79u0ArjyeDlypvwsXLiAgIACZmZkAgBMnTuD2228HACxevBgnT57EyZMnLce54YYbrOrVy8sLaWlpVrOTz58/H++99x7uvPNOvPTSS9DpdFi6dCnq6+uRlZVl2U6r1WLZsmWYOnUqnnjiCTz00EM4efIknn32WYwZMwbjxo2z+gzuPBdRt+OxGZKoy1mwYIGIiYkRGo3GalK8tLQ0qwn6ZBPs5efnCwDi73//u9X6NWvWCADi4MGDVuu3bt0qRo4cKYKDg4VOpxMJCQni3nvvFZ999lmHfD5bFi5cKABIl2snBjQYDOLhhx8W4eHhwtfXVwwfPlzs3LnTapvmepAtCxcutNo+LS2tzeeXsXc+R8sl0zzjt63l2kkjL168KJ544gkRFRUlvL29RZ8+fcQrr7xiNQnj1VqbvLE5bmtJSEiwbNd8zbW1nABsTkL5ww8/iIkTJ4rg4GDh7+8vRo8eLQ4dOmSzbBs3bhQ33nij8PHxEXq9Xjz99NOiqqqqxXbuPBdRd6MI4cRoSyIiIqIugpM3EhEREYFJEREREREAJkVEREREAJgUEREREQFgUkREREQEgEkREREREQBO3tiC2WzGuXPnEBQUxJcjEhGRlBACVVVViImJ6dDXztTX18NoNLrkWD4+PvD19XXJsboiJkXXOHfunGWafCIiotaUlJQ49W5Ce+rr69HDLxC1MLnkeHq9HkVFRUyMJJgUXaP5tRQlJSWW1zMQOcL00zfSWGPPG6SxnaerpLEefi1fNNvMz1v+l2p0oLc0FqzTSmNetS3fvG6htV0WbXWZzfUAIM6fkcbqvz8sjVWfOS+NhQ7sJ43VFZdIY36J8u+BJiBQGjNXVUhj2vAoaUy5Tv4ePeEl/75W6eTvrKtuNEtjpVXyloVhCeHSWFdg+rZAHgyQf3ahC7C53uxtO4GoqqrG9YNu7pDXGTUzGo2ohQn/hevg4+SIFyPM2GA4C6PRyKRIgknRNZq7zIKDg5kUkVNMQfIba6Oda8tfvhsC/OU3T39veXITFORgUqS102QvuZFrlVrpLqLOXxrz8Zf/klZ85Z872N9PGvP2s/1yXgDwC5DvpwmQl9NslteJNlC+n2LnehBe8nLCV36tKEZ5UlSJBmmsq/9uMwXaTm4AAIF2vg/SpEh+rQBwy1ALP2jgoziXFGn5/opWMSkiIiJSOa2iQOtk8qWFcuXtfSTFpIiIiEjlNAqgdbJBSgMwKWoFH8knIiIiAluKiIiIVM9l3WdkF5MiIiIildO6oPtM/kgFNWP3GRERERHYUkRERKR67D5zj07TUvTmm2/ixhtvtMwflJqaik8++cQSF0IgKysLMTEx8PPzQ3p6Oo4dO+bBEhMREblGc/eZs0t77NmzB+PHj0dMTAwURcHWrVvtbl9QUABFUVosx48fd/yDu1mnSYpiY2OxdOlSFBYWorCwEKNGjcKECRMsic+yZcuwYsUKrFq1CgcPHoRer8eYMWNQVSWfHZiIiIhsq6mpwaBBg7Bq1ap27XfixAmUlpZalqSkpA4qoet1mu6z8ePHW329ePFivPnmmzhw4AD69++PlStX4oUXXsDkyZMBAGvXrkVUVBQ2btyIJ554whNFJiIicglPdJ9lZmYiMzOz3eeJjIxEaGhou/dTg07TUnQ1k8mETZs2oaamBqmpqSgqKoLBYEBGRoZlG51Oh7S0NOzbt8/usRoaGlBZWWm1EBERqYmCKzdsZ5bmlOjae15Dg/yVMI4YPHgwoqOjMXr0aOTn57v02B2tUyVFR48eRWBgIHQ6HWbOnIktW7agf//+MBgMAICoKOuXMUZFRVliMtnZ2QgJCbEscXFxHVZ+IiIiT4uLi7O672VnZ7vkuNHR0Xj77beRm5uLzZs3Izk5GaNHj8aePXtccnx36DTdZwCQnJyMI0eO4PLly8jNzcW0adOwe/duS/zal/IJIVp9Ud+CBQswb948y9eVlZVMjIiISFVc2X1WUlJi9VJgnc7OC4nbITk5GcnJyZavU1NTUVJSguXLl+O2225zyTk6WqdKinx8fNC7d28AwNChQ3Hw4EG89tpr+P3vfw8AMBgMiI6OtmxfVlbWovXoWjqdzmUXBBERUUdw5eSNzU9xu8Pw4cOxfv16t5zLFTpV99m1hBBoaGhAYmIi9Ho9du7caYkZjUbs3r0bI0aM8GAJiYiInHclKVKcXNxf7sOHD1s1Vqhdp2kpev7555GZmYm4uDhUVVVh06ZNKCgowKeffgpFUTBnzhwsWbIESUlJSEpKwpIlS+Dv748pU6Z4uuhERESdTnV1NX744QfL10VFRThy5AjCw8MRHx+PBQsW4OzZs1i3bh0AYOXKlejVqxcGDBgAo9GI9evXIzc3F7m5uZ76CO3WaZKin3/+GVOnTkVpaSlCQkJw44034tNPP8WYMWMAAM8++yzq6uowa9YslJeXY9iwYcjLy0NQUJCHS05EROQcT7z7rLCwECNHjrR83Tz+dtq0acjJyUFpaSmKi4stcaPRiPnz5+Ps2bPw8/PDgAEDsG3bNtxxxx3OFdyNOk1S9O6779qNK4qCrKwsZGVluadAREREbuKJeYrS09MhhJDGc3JyrL5+9tln8eyzzzpSNNXo1GOKiIiIiFyl07QUERERdVcaF3SfsRWkdUyKiIiIVM4T3WfdERNHIiIiIrCliIiISPU88fRZd8SkiIiISOWYFLkHu8+IiIiIwJYiIiIi1eNAa/dgUkQWx3+ulMb6Rrnn5YFdiSk0Rho7VWmSxsx2JkuzFwvwkTf89hBV0lgTQqWxbxsCpTG9t7fN9eFB8pcwG3v0lsZ0cQOkMW//MGlMaWqQxrzS5TPaKyb5frWKjzRWb5J/D+oazdLYpTr599xPkX/vfOrl+/na6U+JDpR/hlMX5NfD9RGufwvATKWXNPa7p4dLY9dljpLGfMf9RhprMhRLY17JkdKY2dvP5npNY51kfb30WK6mhQu6z+SXLv0HkyIiIiKV07igpUjj5P7dAccUEREREYEtRURERKrnkqfP2FDUKiZFREREKueSgdbsPmsVu8+IiIiIwJYiIiIi1WP3mXswKSIiIlI5dp+5B7vPiIiIiMCWIiIiItXTKIrT8wxxnqLWMSkiIiJSOUWrQNE4l9QoTIpaxe4zIiIiIrCliIiISPU0WgUaJ1uK2H3WOiZFREREaqfVQNE42bmj8I2wrWFSREREpHKKRoHi5ERDCthS1BqOKSIiIiICW4qIiIhUT6NVoHGypUjDlqJWMSkiIiJSOUXj/JgiRXBMUWvYfUZEREQEthQRERGpHrvP3INJERERkcopWj595g7sPiMiIiICW4qIiIhU70pLkZMDrWF2UWm6LiZFREREKscxRe7B7jMiIiIisKWIiIhI9RRFgeLkC2EVM1uKWsOkiIiISOU0Wg00To4p0gh2DrWGSREREZHKueSRfMGWotYwKSILs4pmgC++VC2NxYcHurEkQP32N6Ux48+l0tipjHnSWG1jozTmbeevwQh/H2nMLiF/6sRkZ+r/3mE6aaxJcsE0aYOl+xib5OU4i3BprPxikzTWaNLKY2b5dRSi85bGqo218v18Hfu1WdtoksbC/eSfITpAfj5t9XmHytIQ0FMaM//4lTSmueFmh863Wpx2aD/j5TJ57NI5aUx3+yMOna+ppsrmepNfiO31qHToPKReTIqIiIhUji1F7sGkiIiISOU4psg9WENEREREYEsRERGR+rmg+wzsPmsVkyIiIiKV0ygKNE7OU6RRmBS1ht1nRERERGBLERERkeopWo3zL4Q1sx2kNUyKiIiIVM4lL4Tlaz5axbSRiIiICJ0oKcrOzsZNN92EoKAgREZGYuLEiThx4oTVNkIIZGVlISYmBn5+fkhPT8exY8c8VGIiIiLXaJ680dmlPfbs2YPx48cjJiYGiqJg69atre6ze/duDBkyBL6+vrj++uuxevVqBz+xZ3SapGj37t148sknceDAAezcuRNNTU3IyMhATU2NZZtly5ZhxYoVWLVqFQ4ePAi9Xo8xY8agqsr21O1ERESdQfOYImeX9qipqcGgQYOwatWqNm1fVFSEO+64A7/61a9w+PBhPP/883j66aeRm5vryEf2iE4zpujTTz+1+nrNmjWIjIzEoUOHcNttt0EIgZUrV+KFF17A5MmTAQBr165FVFQUNm7ciCeeeMITxSYiInKaRgsXjClq3/aZmZnIzMxs8/arV69GfHw8Vq5cCQDo168fCgsLsXz5ctxzzz3tO7mHdJqWomtVVFQAAMLDr7xEsqioCAaDARkZGZZtdDod0tLSsG/fPulxGhoaUFlZabUQERF1Vdfe8xoaGlxy3P3791vdgwFg7NixKCwsRKOdl2CrSadMioQQmDdvHm699VakpKQAAAwGAwAgKirKatuoqChLzJbs7GyEhIRYlri4uI4rOBERkQMUjeKSBQDi4uKs7nvZ2dkuKaPBYLB5D25qasKFCxdcco6O1mm6z642e/ZsfPPNN9i7d2+LmHLNjJ1CiBbrrrZgwQLMmzfP8nVlZSUTIyIiUhWNxgUvhDVd2b+kpATBwcGW9TqdzqnjXs3WPdjWerXqdEnRU089hY8//hh79uxBbGysZb1erwdwJVONjo62rC8rK2uRuV5Np9O59IIgIiJSs+DgYKukyFX0en2LnpmysjJ4eXmhR48eLj9fR+g03WdCCMyePRubN2/Grl27kJiYaBVPTEyEXq/Hzp07LeuMRiN2796NESNGuLu4RERELuOJR/LbKzU11eoeDAB5eXkYOnQovL29O/TcrtJpWoqefPJJbNy4ER999BGCgoIs2WhISAj8/PygKArmzJmDJUuWICkpCUlJSViyZAn8/f0xZcoUD5eeiIjIcS55zUc796+ursYPP/xg+bqoqAhHjhxBeHg44uPjsWDBApw9exbr1q0DAMycOROrVq3CvHnzMGPGDOzfvx/vvvsu3n//fafK7U6dJil68803AQDp6elW69esWYPp06cDAJ599lnU1dVh1qxZKC8vx7Bhw5CXl4egoCA3l5aIiKhzKywsxMiRIy1fN4+/nTZtGnJyclBaWori4mJLPDExEdu3b8fcuXPx+uuvIyYmBn/60586zeP4QCdKipoHa9mjKAqysrKQlZXV8QUiIiJyE0WjgaJxsqWonfunp6fbvffm5OS0WJeWloavv/66vUVTjU6TFBEREXVXGq0Lnj5zcv/ugDVEREREBLYUERERqZ8LBlqDLUWtYlJERESkcorGBU+fOTkmqTtgUkRERKRynhho3R2xhoiIiIjAliIiIiLVuzJ5o9bJY5hcVJqui0lRF3SirFIa89bIp3n3sTMFfHVtnTQW6O/XtoJd41uDvJznqhqkMV875YwMCXCoLHalPywNBfvLzxdTUSONXaqXN9IG6eQ/lkE6+X5RPnZ+4Zl8pCGdsUq+X5NRGvL2C7G5/kK9fF6Tc1Xy44X6yj93TJC8/Fo7L5o8bKiWxny95PUVFyJ/H2J9k/zzhejkN60bwuTHNJrkx6xqlMcCgyKlsTOVjdJYT7P8mI2RfaQxr9NHpDHjkXxpzGfwKGlMmzBIvl+o/PPV1tVLYw12fl9VG83y80mSDi/J7816N95CPTGjdXfEGiIiIiICW4qIiIhUT6PRQOPkQGln9+8OmBQRERGpHLvP3IM1RERERAS2FBEREakeW4rcg0kRERGRyimKCyZvVJgUtYY1RERERAS2FBEREakeu8/cg0kRERGRyjEpcg8mRURERCqn0WqgcTKpcXb/7oA1RERERAS2FBEREameolGcf/rMzrsv6QomRURERCrHMUXuwRoiIiIiAluKiIiIVI8tRe7BpIiIiEjlOKO1e7CGiIiIiMCWIiIiItVTtFpotFqnj0H2MSkiIiJSOY4pcg/WEBERERHYUkRERKR6bClyDyZFREREKqdoXPD0mZP7dwdMioiIiFSOLUXuwRoiIiIiAluKiIiIVE/RKM63FPGFsK1iUkRERKRyHFPkHqwhIiIiIrCliIiISPUUjRaKxskZrZ3cvztgUtQF+drpdy6tNkpjPfy8pbEqo1kaC/RvW7muFeYr/wEN85UfVOfl3gbOykY7wdoaaSgyJEAaM0O+X4hOXi+axnppTHv5rDRmCouXxoyK/NeAyUtIY9UNtq+JSqNJuo+vt/x75+clH+9Qb5KXo7xO/g0qq2mQxpJ6+Nkpi7ycgfIfE9Q1yctpqG6Sxox2Pp8Z8liAUV7O87Xyn3UvjY80Vucl/4DlvtdLY71vjZDGyrzCpLHvfjgvjY3q3VMau1Qnr8/Y8EBpLMjB31e2KCZ7vxxcTKO9sjh7DLKL3WdEREREYEsRERGR+mk0VxZnj0F2MSkiIiJSOUWrdfot987u3x0wbSQiIiICW4qIiIjUjwOt3YJJERERkdppNC5Iitg51BomRURERCrHGa3dgzVEREREBLYUERERqZ/igjFFCscUtaZTtRTt2bMH48ePR0xMDBRFwdatW63iQghkZWUhJiYGfn5+SE9Px7FjxzxTWCIiIldpHmjt7NJOb7zxBhITE+Hr64shQ4bgiy++kG5bUFAARVFaLMePH3fmk7tVp0qKampqMGjQIKxatcpmfNmyZVixYgVWrVqFgwcPQq/XY8yYMaiqqnJzSYmIiDq3Dz74AHPmzMELL7yAw4cP41e/+hUyMzNRXFxsd78TJ06gtLTUsiQlJbmpxM7rVN1nmZmZyMzMtBkTQmDlypV44YUXMHnyZADA2rVrERUVhY0bN+KJJ55wZ1GJiIhcxhMDrVesWIHHHnsMjz/+OABg5cqV2LFjB958801kZ2dL94uMjERoaKgzRfWYTtVSZE9RUREMBgMyMjIs63Q6HdLS0rBv3z7pfg0NDaisrLRaiIiIVMWF3WfX3vMaGlq+ONloNOLQoUNW91QAyMjIsHtPBYDBgwcjOjoao0ePRn5+vuvqwA26TFJkMBgAAFFRUVbro6KiLDFbsrOzERISYlni4uI6tJxERESeFBcXZ3Xfs9Xqc+HCBZhMpnbdU6Ojo/H2228jNzcXmzdvRnJyMkaPHo09e/Z0yOfoCJ2q+6wtFEWx+loI0WLd1RYsWIB58+ZZvq6srGRiRERE6uLCyRtLSkoQHBxsWa3T6aS7tOeempycjOTkZMvXqampKCkpwfLly3Hbbbc5U3K36TJJkV6vB3ClxSg6OtqyvqysrEWmezWdTmf3giAiIvI0V74QNjg42CopsiUiIgJarbZFq1Br99RrDR8+HOvXr29/YT2ky3SfJSYmQq/XY+fOnZZ1RqMRu3fvxogRIzxYMiIios7Fx8cHQ4YMsbqnAsDOnTvbdU89fPiwVUOF2nWqlqLq6mr88MMPlq+Liopw5MgRhIeHIz4+HnPmzMGSJUuQlJSEpKQkLFmyBP7+/pgyZYoHS01EROQkjcb5d5e1c/958+Zh6tSpGDp0KFJTU/H222+juLgYM2fOBHBl+MnZs2exbt06AFeeTuvVqxcGDBgAo9GI9evXIzc3F7m5uc6V2406VVJUWFiIkSNHWr5uHgs0bdo05OTk4Nlnn0VdXR1mzZqF8vJyDBs2DHl5eQgKCvJUkYmIiJzn4OSLLY7RDg888AAuXryIl156CaWlpUhJScH27duRkJAAACgtLbWas8hoNGL+/Pk4e/Ys/Pz8MGDAAGzbtg133HGHc+V2I0UIITxdCDWprKxESEgIKioqWu1zVaufLlZLY6XVRmmsh5+3NBboI/8LIzo0oG0Fu7Ysl2sc2s/fW16WkAA/h45pT1mFvJzB8iqDr7+8Xgx2jhmik//i0jTWS2PairPSmCksXhozKvK/jUxm+a+HaqPZ5vpKyXoAMNn5dRNm53PXm+T7ldc1SWNHy+QTt94SHyqNBfvIy6KVP7eBuiZ5OY12PoO9mBnyWICdnwWDnZ/16EAfaczXS37M8nqTNNbbSz6dyUWvMGnsuwu10tio3j2lsTOX5L/nYsMDpTFXcsf9ovkcFz7fgOAAf+eOVVOLiNH/1anvbx2ty4wpIiIiInJGp+o+IyIi6pYUF4wpUtgO0homRURERCqnaLRQnBxT5Oz+3QHTRiIiIiKwpYiIiEj9XDijNckxKSIiIlI7D8xT1B2xhoiIiIjAlqIuyc77b9ErVP6etwY7c6zYm5ultk4+d47Zzrw031+sk8aG6+Xl1AU4Ni/S+Ur5fCghJvkcK5FheofOV18n/3xhWvncOtpLZ6Qxc0APeSxI/j4iTW25NKbzks9Zo9iZF8kr0Pb57M1TFB8sP5ehplEa++qM/PsTZmd+rYs18rl6vjgtrxONRn7B1zbK5+qpM8pjWjvHtDcflD0/2Znj5/13/yGNpd97uzQ2J723NHb6svya/red+Y1uCJN/b2OC5D/rh0ouS2MBduZO64pc+e4zkmNSREREpHYemNG6O2JSREREpHZMityie7U/EhEREUmwpYiIiEjlFI0GipNPjzm7f3fApIiIiEjtFBd0nynsPmsN00YiIiIisKWIiIhI/RTF+Re62puvhQAwKSIiIlI/ReOCpIidQ61hDRERERGBLUVERESqJxQNhJMtPc7u3x0wKSIiIlI7dp+5BWuIiIiICGwpIiIiUj9Fcf7pMT591iomRURERGqn0VxZnD0G2cWkiIiISOU40No9WENEREREYEsRERGR+vHpM7dgUkRERKR2TIrcgjVEREREBLYUERERqR9bityCSREREZHKCUVxwdNnnKeoNUwbiYiIiMCWIiIiIvVj95lbMCkiIiJSO77mwy2YNhIRERGBLUVERETqx+4zt2BSREREpHJ895l7MCkiIiJSO0Xj/FvumRS1iklRJ1VRUyeNRfjJv60aYZLGmry10tjlevl+1UZ5zM9b/kM4PFJeTsVYI401Vl+QxryjEqWxnsH+0hhgL+YYpalBGtNWn5fGvGKSpbFGw4/SmDmghzRmDOgpjdnT4BUsjZ2+aPvznb5cK91nV5G8TgyX66WxuHD596fwp3Jp7NT5amnMZBbSmD0VdsoZEuorjfl4yX8WKqqN0lhjg/znSwj5Z4i9cbA09r+FZ6Wx+z85KI2FxyVIY/peodJYaKCPNBbk6y2NDU6QH/O/b7leGiNyFJMiIiIiteOYIrdgUkRERKR2TIrcgjVEREREBLYUERERqR9bityCSREREZHK8YWw7sG0kYiIiDqddevWoaGh5ROtRqMR69atc+iYTIqIiIjUrrn7zNmlC3nkkUdQUVHRYn1VVRUeeeQRh47pUA09+uijqKqqarG+pqYGjz76qEMFISIiIonmF8I6u3QhQggoNj7TmTNnEBIS4tAxHRpTtHbtWixduhRBQUFW6+vq6rBu3Tr89a9/dagwRERERPYMHjwYiqJAURSMHj0aXl7/l8qYTCYUFRVh3LhxDh27XS1FlZWVqKiogBACVVVVqKystCzl5eXYvn07IiMjHSqIK73xxhtITEyEr68vhgwZgi+++MLTRSIiInKch7rP2ns/3b17N4YMGQJfX19cf/31WL16taOfWGrixImYMGEChBAYO3YsJkyYYFkefPBBvPXWW1i/fr1Dx25XS1FoaKglO+vTp0+LuKIoWLRokUMFcZUPPvgAc+bMwRtvvIFbbrkFb731FjIzM/Htt98iPj7eo2UjIiJyhCdeCNve+2lRURHuuOMOzJgxA+vXr8e//vUvzJo1Cz179sQ999zjVNmvtnDhQgBAr1698OCDD0Kn07ns2O1KivLz8yGEwKhRo5Cbm4vw8HBLzMfHBwkJCYiJiXFZ4RyxYsUKPPbYY3j88ccBACtXrsSOHTvw5ptvIjs726NlIyIicogH5ilq7/109erViI+Px8qVKwEA/fr1Q2FhIZYvX+7SpKjZqFGjcP78ecTGxgIAvvrqK2zcuBH9+/fHb37zG4eO2a6kKC0tDcCVbDAuLg4aZ9/Y62JGoxGHDh3Cc889Z7U+IyMD+/bts7lPQ0OD1SN9lZWVHVpGIiIiT7r2PqfT6Vq0tjhyP92/fz8yMjKs1o0dOxbvvvsuGhsb4e0tf/mvI6ZMmYLf/OY3mDp1KgwGA26//XakpKRg/fr1MBgMePHFF9t9TIeymoSEBGg0GtTW1uL48eP45ptvrBZPuXDhAkwmE6KioqzWR0VFwWAw2NwnOzsbISEhliUuLs4dRSUiImqzK5M3Or8AQFxcnNV9z1arjyP3U4PBYHP7pqYmXLhwwUU18X/+/e9/4+abbwYA/O1vf8PAgQOxb98+bNy4ETk5OQ4d06Gnz86fP49HHnkEn3zyic24yWRyqDCucu0jerLH9gBgwYIFmDdvnuXryspKJkZERKQqQlxZnD0GAJSUlCA4ONiy3t6YnPbcT2Xb21rvCo2NjZayf/bZZ7j77rsBAH379kVpaalDx3SopWjOnDkoLy/HgQMH4Ofnh08//RRr165FUlISPv74Y4cK4goRERHQarUtstiysrIW2WsznU6H4OBgq4WIiKiruvaeZyspcuR+qtfrbW7v5eWFHj16uO4D/MeAAQOwevVqfPHFF9i5c6flMfxz5845fD6HkqJdu3bhj3/8I2666SZoNBokJCTg17/+NZYtW+bRwcw+Pj4YMmQIdu7cabV+586dGDFihIdKRURE5ByzEC5Z2sqR+2lqamqL7fPy8jB06FCXjycCgJdffhlvvfUW0tPT8dBDD2HQoEEAgI8//tjSrdZeDnWf1dTUWOYjCg8Px/nz59GnTx8MHDgQX3/9tUMFcZV58+Zh6tSpGDp0KFJTU/H222+juLgYM2fO9Gi5iIiIHCX+szh7jPZo7X66YMECnD171vKesZkzZ2LVqlWYN28eZsyYgf379+Pdd9/F+++/72TJbRs2bBjOnDkDo9GIsLAw/PTTT9iyZQsSEhIwd+5ch47pUFKUnJyMEydOoFevXvjFL36Bt956C7169cLq1asRHR3tUEFc5YEHHsDFixfx0ksvobS0FCkpKdi+fTsSEhI8Wi4iIqLOpLX7aWlpKYqLiy3bJyYmYvv27Zg7dy5ef/11xMTE4E9/+lOHPI4PABMmTMDkyZMxc+ZMXL58GTfffDN8fHxw4cIFmM1m/Pa3v233MR1KiubMmWMZxLRw4UKMHTsW69evh4+PD9auXevIIV1q1qxZmDVrlqeLQURE5BJmcWVx9hjtZe9+ausJr7S0NLf1GH399df44x//CAD48MMPodfrcfjwYeTm5uLFF190X1L0X//1X5b/Dx48GKdPn8bx48cRHx+PiIgIRw5JREREEkIIy5NczhyjK6mtrbW8gzUvLw+TJ0+GRqPB8OHD8dNPPzl0zDYnRVc/tt6aFStWOFQYIiIiorbo3bs3tm7dikmTJmHHjh2WcURlZWUOP0ne5qTo8OHDbdquI+YiICIi6s481X2mZi+++CKmTJmCuXPnYvTo0UhNTQVwpdVo8ODBDh2zzUlRfn6+QycgIiIi53WxnMZp9957L2699VaUlpZaHscHgNGjR2PSpEkOHdOhMUVERETkPmwpsk2v10Ov11utc3SOIsDByRuJiIiIuhq2FBEREakcnz5zDyZFKlZ6uUYa8/OSD2g3muQXvr+pThrT6gKlMR+t/HyNZmkI4ZeLpTGzX4g01uQbKo15KfIGzobqCmlMFyg/nz2Nhh+lMXOgfAoKs7evNOZdJy9n06Ft8v2G3CmN2aM9/Kk0poRGSmMiqp80Fuantbnez1t+Hdm5NJEUHiCNBelsn6s1CRH+8rLY6UuoqG2UxobcGiqN/XBR/jNrjz5Ifq188YP87eJJUfK6PnVeXpZQP/krF0L8Bzp0zOp6eZ0lRQVJY6UV9dLY+/mnpLH/vuV6aawrMv9ncfYYZB+7z4iIiIjAliIiIiLVE+LK4uwxyD4mRURERCrHp8/cg91nRERERGBLERERkerx6TP3YFJERESkcnz6zD3YfUZEREQEthQRERGpnoALnj5zSUm6NiZFREREKmcWAmYnsyJn9+8OmBQRERGpnIDzLT1MiVrHMUVEREREYEsRERGR6nHyRvdgUkRERKR2LnjNB/vPWsfuMyIiIiKwpYiIiEj1zBAwO9nU4+z+3QGTIiIiIpUTLug+4xP5rWP3GRERERHYUkRERKR6fPrMPZgUERERqRy7z9yD3WdEREREYEsRERGR6vHpM/dgUkRERKRy7D5zDyZFREREKmcWwum33Du7f3fAMUVEREREYEsRERGR6pnMVxZnj0H2MSkiIiJSOXafuQe7z4iIiIjAliIiIiLVMwsBE1uKOhyTIhWrNJqksWqjIo3Zm4tCq+ikMWOdURrz0crP52snZgrWS2OKsUYaa7QzH71GFyQ/prlJGjNeOiffr7FBGoPiWIOq0WTn+xAaI41pfPyksabSk/ITlhVJQ16Dx0ljdf98XRrzrymXxmL1yTbXm0IjpPtcHxImjdlT3iAfDKFR7Pws2LkJNNr5/ti7/uzpFym/NhvtDOgor2+UxvpGy48ZrJP/Ch9wXbA0Fu7rLY1dF+wrjaVf30MaO1tZL41FBsh/71Q0yH9mo0PkZek392Np7Ls/3i2NdVZXXvPhbFLkosJ0Yew+IyIiIgJbioiIiFSPT5+5B5MiIiIilePTZ+7B7jMiIiIisKWIiIhI9UwuePrM2f27AyZFREREKmeG80+PcUhR65gUERERqZzJLGByMitydv/ugGOKiIiIiNCJkqLFixdjxIgR8Pf3R2hoqM1tiouLMX78eAQEBCAiIgJPP/00jEb5hIRERESdgfjP02fOLIJjilrVaZIio9GI++67D7/97W9txk0mE+68807U1NRg79692LRpE3Jzc/HMM8+4uaRERESuZRKuWTpKeXk5pk6dipCQEISEhGDq1Km4fPmy3X2mT58ORVGsluHDh3dcIdug04wpWrRoEQAgJyfHZjwvLw/ffvstSkpKEBNz5RUKr776KqZPn47FixcjOFg+5T0RERE5bsqUKThz5gw+/fRTAMBvfvMbTJ06Ff/4xz/s7jdu3DisWbPG8rWPj0+HlrM1nSYpas3+/fuRkpJiSYgAYOzYsWhoaMChQ4cwcuRID5aOiIjIcWqevPG7777Dp59+igMHDmDYsGEAgL/85S9ITU3FiRMnkJxs+z2JAKDT6aDXy9+R6W6dpvusNQaDAVFRUVbrwsLC4OPjA4PBIN2voaEBlZWVVgsREZGaND995uwCoMU9r6HBzgux22D//v0ICQmxJEQAMHz4cISEhGDfvn129y0oKEBkZCT69OmDGTNmoKyszKmyOMujSVFWVlaL/sRrl8LCwjYfT7HxtmwhhM31zbKzsy19oCEhIYiLi3PosxAREXUGcXFxVve97Oxsp45nMBgQGRnZYn1kZKTdRonMzExs2LABu3btwquvvoqDBw9i1KhRTidpzvBo99ns2bPx4IMP2t2mV69ebTqWXq/Hl19+abWuvLwcjY2NLVqQrrZgwQLMmzfP8nVlZSUTIyIiUhVXdp+VlJRYjbPV6XQ2t8/KyrKM55U5ePAgAMcaJR544AHL/1NSUjB06FAkJCRg27ZtmDx5st3zdhSPJkURERGIiIhwybFSU1OxePFilJaWIjo6GsCVwdc6nQ5DhgyR7qfT6aQXBBERkRq44umx5v2Dg4Pb9PBRWxsuvvnmG/z8888tYufPn7fbKHGt6OhoJCQk4OTJk23ex9U6zUDr4uJiXLp0CcXFxTCZTDhy5AgAoHfv3ggMDERGRgb69++PqVOn4pVXXsGlS5cwf/58zJgxg0+eERERtVNbGy5SU1NRUVGBr776CjfffDMA4Msvv0RFRQVGjBjR5vNdvHgRJSUlloYNT+g0A61ffPFFDB48GAsXLkR1dTUGDx6MwYMHW8YcabVabNu2Db6+vrjllltw//33Y+LEiVi+fLmHS05EROQcZydudEX3m0y/fv0wbtw4zJgxAwcOHMCBAwcwY8YM3HXXXVZPnvXt2xdbtmwBAFRXV2P+/PnYv38/Tp8+jYKCAowfPx4RERGYNGlSh5SzLTpNS1FOTo50jqJm8fHx+Oc//+meAhEREbmJ2SxgdvLdZc7ub8+GDRvw9NNPIyMjAwBw9913Y9WqVVbbnDhxAhUVFQCuNGQcPXoU69atw+XLlxEdHY2RI0figw8+QFBQUIeVszWdJikiIiLqrswuGFPUke+DDQ8Px/r16+1uc/VrRvz8/LBjx46OK5CDOk33GREREVFHYksRERGRyql5RuuuhEkRERGRypmEgMnJpMbZ/bsDdp8RERERgS1FREREqqf2p8+6CiZFREREKmeCC2a0dklJujZ2nxERERGBLUWqlhzp2OtJfjhfJY2V1zVJYwE+Wmns5MVaaSwyQP7uuIoGed59Q2iYNOZfcU4aE7oAeUzrI495+0tjOLFfGjJXXZbGvPoPl8ZCYgfIzxfgJ4+Fx8hj9kQnSUOX335eGqs+e14aC6sql8a8ys7aXO8d0kO6jxIZL48Z66SxiGC9NBYeIb+OzjfIX0ZZ1yT/s1tr58/Fkgr5G7yLa+Qxf2/5z1dSuPyatuf0ZXmdBerkv97tlaXKKG9PsHe+MD9vacxefSb1kP8s+HvLdwy/N0V+0C6IT5+5B5MiIiIilePTZ+7B7jMiIiIisKWIiIhI9cxmAROfPutwTIqIiIhUzuSCpMjZ/bsDJkVEREQqx6TIPTimiIiIiAhsKSIiIlI9k9n5lh6T2UWF6cKYFBEREakcu8/cg91nRERERGBLERERkeqxpcg9mBQRERGpHOcpcg92nxERERGBLUVERESqZxIu6D7ju89axaSIiIhI5TimyD3YfUZEREQEthQRERGpHluK3INJERERkco1mQW0TiY1TUyKWsWkiIiISOXYUuQeHFNEREREBLYUERERqR4nb3QPJkVEREQqZxLC6XmGOE9R69h9RkRERAS2FBEREakeB1q7B5MiIiIilWNS5B7sPiMiIiICW4qIiIhUjy1F7sGkiIiISOVMwgyT2ez0Mcg+JkUedqKsUhpLjgx26Jh1TfILv95OLDrQWxoL85PHqo1NdmLSEJJ01dKY2T9MGvOqOCvfTxckjwX0kMa0sX2kMZ+4gdJYZxH6myXSWNHEDGnMbOd76x9dbnO9l69Ouo/W96A0Vmu4KI0F97leGoOXjzQU03+YNCY08l9/wkt+vYdEJkpjKZH+0tjPNfK6DNFppTGtnUEOPQPk5TxX1SCNNZrkLQYNJvnviMgAeV032mmF+NnOL4KqBpM0FmSnXm66LkQa2/S/8t8RDw66ThojYlJERESkcpy80T2YFBEREamcySyg4ZiiDsekiIiISOWazIDiZFJjZ/QE/QcfySciIiICW4qIiIhUj91n7sGkiIiISOWYFLkHu8+IiIiIwJYiIiIi1WNLkXswKSIiIlI5zlPkHp2i++z06dN47LHHkJiYCD8/P9xwww1YuHAhjEbrWVKLi4sxfvx4BAQEICIiAk8//XSLbYiIiIhs6RRJ0fHjx2E2m/HWW2/h2LFj+OMf/4jVq1fj+eeft2xjMplw5513oqamBnv37sWmTZuQm5uLZ555xoMlJyIicl7zC2GdXTrK4sWLMWLECPj7+yM0NLRN+wghkJWVhZiYGPj5+SE9PR3Hjh3rsDK2RafoPhs3bhzGjRtn+fr666/HiRMn8Oabb2L58uUAgLy8PHz77bcoKSlBTEwMAODVV1/F9OnTsXjxYgQHO/YeMSIiIk8TQkA4mdQI0XFJkdFoxH333YfU1FS8++67bdpn2bJlWLFiBXJyctCnTx/8v//3/zBmzBicOHECQUHy91h2pE7RUmRLRUUFwsPDLV/v378fKSkploQIAMaOHYuGhgYcOnTIE0UkIiLqFhYtWoS5c+di4MC2vUBbCIGVK1fihRdewOTJk5GSkoK1a9eitrYWGzdu7ODSynXKpOjHH3/En//8Z8ycOdOyzmAwICoqymq7sLAw+Pj4wGAwSI/V0NCAyspKq4WIiEhNzGbhkgVAi3teQ0OD2z9PUVERDAYDMjIyLOt0Oh3S0tKwb98+t5enmUeToqysLCiKYncpLCy02ufcuXMYN24c7rvvPjz++ONWMUVRWpxDCGFzfbPs7GyEhIRYlri4ONd8OCIiIhcRQrhkAYC4uDir+152drbbP09zY8W1jRlRUVF2GzI6mkfHFM2ePRsPPvig3W169epl+f+5c+cwcuRIpKam4u2337baTq/X48svv7RaV15ejsbGxhaVfrUFCxZg3rx5lq8rKyuZGBERkaoIswvGFP1n/5KSEqtxtjqdzub2WVlZWLRokd1jHjx4EEOHDnW4TNc2WrTWkNHRPJoURUREICIiok3bnj17FiNHjsSQIUOwZs0aaDTWjVypqalYvHgxSktLER0dDeDK4GudTochQ4ZIj6vT6aQXBBERUVcTHBzcpoeP2ttw0R56vR7AlRaj5ns2AJSVldltyOhoneLps3PnziE9PR3x8fFYvnw5zp8/b4k1V2xGRgb69++PqVOn4pVXXsGlS5cwf/58zJgxg0+eERFRp3b1mCBnjtEe7Wm4aK/ExETo9Xrs3LkTgwcPBnDlCbbdu3fj5Zdf7pBztkWnSIry8vLwww8/4IcffkBsbKxVrLmPVKvVYtu2bZg1axZuueUW+Pn5YcqUKZZH9omIiDorYb6yOHuMjlJcXIxLly6huLgYJpMJR44cAQD07t0bgYGBAIC+ffsiOzsbkyZNgqIomDNnDpYsWYKkpCQkJSVhyZIl8Pf3x5QpUzquoK3oFEnR9OnTMX369Fa3i4+Pxz//+c+OLxARERFZvPjii1i7dq3l6+bWn/z8fKSnpwMATpw4gYqKCss2zz77LOrq6jBr1iyUl5dj2LBhyMvL89gcRUAnSYqIiIi6s6ufHnPmGB0lJycHOTk57Tq/oijIyspCVlZWh5WrvZgUERERqZwnxhR1R51y8kYiIiIiV2NLERERkcq5cp4ikmNSREREpHYuSIrApKhV7D4jIiIiAluKiIiIVM8sBBQnnx4zd+DTZ10FkyIiIiKVE8IFY4qYFLWKSZGHXRfo7dB+F6tqpbGkEPm3dW+1URoL0mmlsaH156Wx0/4x0lhto3wK1R+MPtJYb0X++Yw9k6Qxr4pSaUwo8t7iM/6J0liCNNI1DHjmUWnMVHFRGms4d9bmep8e4dJ9Kr4vksaC4uXvO9IE95DGFL8AaczeFL4Nhz6TH9Nbfm0GR8VLY03nbdcJACSaTdKYaJT/XPr+4jZpLMI/TBq7IUj+e6DcN1IaM9m5b5rs3JRrmuR1XVkv/+w6L8dGcRRX1EtjIbqud2vjQGv34JgiIiIiIrCliIiISPXMZkBxevJGFxWmC2NSREREpHJqf81HV8HuMyIiIiKwpYiIiEj1hNnucwNtPgbZx6SIiIhI5cxm4YIxRew+aw27z4iIiIjAliIiIiLV4zxF7sGkiIiISOWYFLkHu8+IiIiIwJYiIiIi1eMLYd2DSREREZHKsfvMPZgUERERqZwQLkiK2FLUKo4pIiIiIgJbioiIiFRPmIXTky+y+6x1TIqIiIhUji+EdQ92nxERERGBLUVERESqx6fP3INJERERkcqZzQLgC2E7HLvPiIiIiMCWIiIiItUTZhOE2eT0Mcg+JkVEREQqx6TIPdh9RkRERAS2FBEREameMJtd0FJkdlFpui4mRURERConTCYIk5NJkZP7dwdMioiIiFROCBeMKRJMilrDMUVEREREYEuR2zSWnba5Xmcn899dGiKNHTFUSmOz42uksXQ/eR4squTHtCc2QH5MpalRGmv08pPGKpqCpTHRIK+zHsZqacyESGks0Kcb/32gkX92n159pbFLh76xud43vpd0n5DrY+XFCOkhjZmrL8v3k0YA449HpTHv2N7SmOLrLy9LbZV8Pzt12VRxUR6rqZfGxFd50pjWTp1pwuTXe3hCijSmNMrLYvaV/1yaAyOkMWOgr3w/O/MJGk3yoFZRpLGvSx37XaZmfPrMPZgUERERqRyTIvfoxn8eExEREf0fthQRERGpHFuK3INJERERkcpxniL3YPcZEREREdhSREREpHpmswlwsqXIzO6zVjEpIiIiUjmOKXIPdp8RERERgS1FREREqseWIvfoNC1Fd999N+Lj4+Hr64vo6GhMnToV586ds9qmuLgY48ePR0BAACIiIvD000/DaDR6qMREREQu8p8XwjqzgC+EbVWnSYpGjhyJv/3tbzhx4gRyc3Px448/4t5777XETSYT7rzzTtTU1GDv3r3YtGkTcnNz8cwzz3iw1ERERM5rfiGsU0sHvhB28eLFGDFiBPz9/REaGtqmfaZPnw5FUayW4cOHd1gZ26LTdJ/NnTvX8v+EhAQ899xzmDhxIhobG+Ht7Y28vDx8++23KCkpQUxMDADg1VdfxfTp07F48WIEB8vf2UNERESOMxqNuO+++5Camop33323zfuNGzcOa9assXzt4+PTEcVrs06TFF3t0qVL2LBhA0aMGAFvb28AwP79+5GSkmJJiABg7NixaGhowKFDhzBy5EhPFZeIiMgpwmx2+pH8jpy8cdGiRQCAnJycdu2n0+mg1+s7oESO6TTdZwDw+9//HgEBAejRoweKi4vx0UcfWWIGgwFRUVFW24eFhcHHxwcGg0F6zIaGBlRWVlotREREauJ019lVA7Wvvec1NDR47HMVFBQgMjISffr0wYwZM1BWVuaxsgAeToqysrJa9CdeuxQWFlq2/93vfofDhw8jLy8PWq0WDz/8MIQQlriiKC3OIYSwub5ZdnY2QkJCLEtcXJxrPyQREZGKxMXFWd33srOzPVKOzMxMbNiwAbt27cKrr76KgwcPYtSoUR5N0jzafTZ79mw8+OCDdrfp1auX5f8RERGIiIhAnz590K9fP8TFxeHAgQNITU2FXq/Hl19+abVveXk5GhsbW7QgXW3BggWYN2+e5evKykomRkREpCpXus+c6/5q7j4rKSmxGmer0+lsbp+VlWXpFpM5ePAghg4d6lB5HnjgAcv/U1JSMHToUCQkJGDbtm2YPHmyQ8d0lkeTouYkxxHNLUTNGWVqaioWL16M0tJSREdHAwDy8vKg0+kwZMgQ6XF0Op30giAiIlID4YLXfDR3nwUHB7fp4aP2Nlw4Kzo6GgkJCTh58qTLjtlenWKg9VdffYWvvvoKt956K8LCwnDq1Cm8+OKLuOGGG5CamgoAyMjIQP/+/TF16lS88soruHTpEubPn48ZM2bwyTMiIqJ2cqbhwhEXL15ESUmJpWHDEzrFQGs/Pz9s3rwZo0ePRnJyMh599FGkpKRg9+7dllYerVaLbdu2wdfXF7fccgvuv/9+TJw4EcuXL/dw6YmIiJzjyoHWHaG4uBhHjhxBcXExTCYTjhw5giNHjqC6utqyTd++fbFlyxYAQHV1NebPn4/9+/fj9OnTKCgowPjx4xEREYFJkyZ1WDlb0ylaigYOHIhdu3a1ul18fDz++c9/uqFERERE7mM2m6Co+DUfL774ItauXWv5evDgwQCA/Px8pKenAwBOnDiBiooKAFcaMo4ePYp169bh8uXLiI6OxsiRI/HBBx8gKCiow8rZmk6RFBEREZF65eTktDpH0dVPi/v5+WHHjh0dXKr2Y1JERESkcsJkBhQnW4pMHTd5Y1fBpIiIiEjlhHDB02cd+O6zroJJERERkcoJs8n5lqIOHFPUVXSKp8+IiIiIOhpbioiIiFSOLUXuwaSIiIhI5ZgUuQeToms0PzJYWVnp0uM2VlXZDti5SGtq5L2b9TWS4wGorKqVxoQiP6bwapTG7M2P0aSR15XSJH+xX6Od89U3CWns6sc6r+VdVS2NNXrLy1lrkr802Fs0SWNdgbFGfr1otb7SWFW90eZ675o66T7m2nppTOMl3080yq8jjWJnvyb5Nab1kn9uxc71Z66T72eulZelyc5nb6qTfz5vrfxnVustP5/GR15OjZ2fE6VRXk5zo/znxGz2kcaMivxWY5ZXNRpN8mB1rfznsq7azu9HF/5ubz6Wvd9JLmNqhNNnMcl/HugKRbjlu9l5nDlzhi+EJSKiNispKUFsbGyHHLu+vh6JiYkwGAwuOZ5er0dRURF8feV/8HRnTIquYTabce7cOQQFBUFR5H8VqUllZSXi4uJavPm4u2O9tMQ6sY310hLrxLar6yUoKAhVVVWIiYmBRtNxzy3V19fDaLTdMttePj4+TIjsYPfZNTQaTYdl/B2trW8+7m5YLy2xTmxjvbTEOrGtuV5CQkI6/Fy+vr5MZNyEj+QTERERgUkREREREQAmRV2CTqfDwoULodPpPF0UVWG9tMQ6sY310hLrxDbWS9fGgdZEREREYEsREREREQAmRUREREQAmBQRERERAWBSRERERASASVGnkZWVBUVRrBa9Xm+JCyGQlZWFmJgY+Pn5IT09HceOHfNgiTvGnj17MH78eMTExEBRFGzdutUq3pZ6aGhowFNPPYWIiAgEBATg7rvvxpkzZ9z4KVyrtTqZPn16i2tn+PDhVtt0tTrJzs7GTTfdhKCgIERGRmLixIk4ceKE1Tbd8VppS710x+vlzTffxI033miZkDE1NRWffPKJJd4dr5XuiklRJzJgwACUlpZalqNHj1piy5Ytw4oVK7Bq1SocPHgQer0eY8aMQZXsRbSdVE1NDQYNGoRVq1bZjLelHubMmYMtW7Zg06ZN2Lt3L6qrq3HXXXfBZOqcb5BurU4AYNy4cVbXzvbt263iXa1Odu/ejSeffBIHDhzAzp070dTUhIyMDNTU1Fi26Y7XSlvqBeh+10tsbCyWLl2KwsJCFBYWYtSoUZgwYYIl8emO10q3JahTWLhwoRg0aJDNmNlsFnq9XixdutSyrr6+XoSEhIjVq1e7qYTuB0Bs2bLF8nVb6uHy5cvC29tbbNq0ybLN2bNnhUajEZ9++qnbyt5Rrq0TIYSYNm2amDBhgnSfrl4nQghRVlYmAIjdu3cLIXitNLu2XoTg9dIsLCxMvPPOO7xWuhm2FHUiJ0+eRExMDBITE/Hggw/i1KlTAICioiIYDAZkZGRYttXpdEhLS8O+ffs8VVy3a0s9HDp0CI2NjVbbxMTEICUlpUvXVUFBASIjI9GnTx/MmDEDZWVlllh3qJOKigoAQHh4OABeK82urZdm3fl6MZlM2LRpE2pqapCamsprpZthUtRJDBs2DOvWrcOOHTvwl7/8BQaDASNGjMDFixdhMBgAAFFRUVb7REVFWWLdQVvqwWAwwMfHB2FhYdJtuprMzExs2LABu3btwquvvoqDBw9i1KhRaGhoAND160QIgXnz5uHWW29FSkoKAF4rgO16Abrv9XL06FEEBgZCp9Nh5syZ2LJlC/r3789rpZvx8nQBqG0yMzMt/x84cCBSU1Nxww03YO3atZZBkIqiWO0jhGixrjtwpB66cl098MADlv+npKRg6NChSEhIwLZt2zB58mTpfl2lTmbPno1vvvkGe/fubRHrzteKrF666/WSnJyMI0eO4PLly8jNzcW0adOwe/duS7w7XyvdCVuKOqmAgAAMHDgQJ0+etDyFdu1fJGVlZS3+uunK2lIPer0eRqMR5eXl0m26uujoaCQkJODkyZMAunadPPXUU/j444+Rn5+P2NhYy/rufq3I6sWW7nK9+Pj4oHfv3hg6dCiys7MxaNAgvPbaa93+WulumBR1Ug0NDfjuu+8QHR2NxMRE6PV67Ny50xI3Go3YvXs3RowY4cFSuldb6mHIkCHw9va22qa0tBT//ve/u01dXbx4ESUlJYiOjgbQNetECIHZs2dj8+bN2LVrFxITE63i3fVaaa1ebOkO14stQgg0NDR022ul2/LE6G5qv2eeeUYUFBSIU6dOiQMHDoi77rpLBAUFidOnTwshhFi6dKkICQkRmzdvFkePHhUPPfSQiI6OFpWVlR4uuWtVVVWJw4cPi8OHDwsAYsWKFeLw4cPip59+EkK0rR5mzpwpYmNjxWeffSa+/vprMWrUKDFo0CDR1NTkqY/lFHt1UlVVJZ555hmxb98+UVRUJPLz80Vqaqq47rrrunSd/Pa3vxUhISGioKBAlJaWWpba2lrLNt3xWmmtXrrr9bJgwQKxZ88eUVRUJL755hvx/PPPC41GI/Ly8oQQ3fNa6a6YFHUSDzzwgIiOjhbe3t4iJiZGTJ48WRw7dswSN5vNYuHChUKv1wudTiduu+02cfToUQ+WuGPk5+cLAC2WadOmCSHaVg91dXVi9uzZIjw8XPj5+Ym77rpLFBcXe+DTuIa9OqmtrRUZGRmiZ8+ewtvbW8THx4tp06a1+LxdrU5s1QcAsWbNGss23fFaaa1euuv18uijj4qEhATh4+MjevbsKUaPHm1JiITontdKd6UIIYT72qWIiIiI1IljioiIiIjApIiIiIgIAJMiIiIiIgBMioiIiIgAMCkiIiIiAsCkiIiIiAgAkyIiIiIiAEyKiAhAeno65syZ4+liEBF5FJMiIiIiIjApIiIiIgLApIiIrlFeXo6HH34YYWFh8Pf3R2ZmJk6ePGmJ5+TkIDQ0FDt27EC/fv0QGBiIcePGobS01IOlJiJyHpMiIrIyffp0FBYW4uOPP8b+/fshhMAdd9yBxsZGyza1tbVYvnw53nvvPezZswfFxcWYP3++B0tNROQ8L08XgIjU4+TJk/j444/xr3/9CyNGjAAAbNiwAXFxcdi6dSvuu+8+AEBjYyNWr16NG264AQAwe/ZsvPTSSx4rNxGRK7CliIgsvvvuO3h5eWHYsGGWdT169EBycjK+++47yzp/f39LQgQA0dHRKCsrc2tZiYhcjUkREVkIIaTrFUWxfO3t7W0VVxRFui8RUWfBpIiILPr374+mpiZ8+eWXlnUXL17E999/j379+nmwZEREHY9JERFZJCUlYcKECZgxYwb27t2L//3f/8Wvf/1rXHfddZgwYYKni0dE1KGYFBGRlTVr1mDIkCG46667kJqaCiEEtm/f3qLLjIioq1EEBwIQERERsaWIiIiICGBSRERERASASRERERERACZFRERERACYFBEREREBYFJEREREBIBJEREREREAJkVEREREAJgUEREREQFgUkREREQEgEkREREREQAmRUREREQAgP8P6MqpOUA14xQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ssta1[-1].plot()" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "### Reason\n", "\n", "Removing the linear trend has a great impact on the Indian Ocean, because the warming trend of the Indian Ocean is more obvious" ] } ], "metadata": { "kernelspec": { "display_name": "sacpy", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.16" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }