{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 74,
   "id": "361279f7-fb2c-4b7a-9529-1cec58739dc3",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB3s0lEQVR4nO3df1wUdf4H8NfyQwQFylKQ/BEWqWG6pWXaD+iHlJX9sLurxNQs+2F2kXX2tbpCr6CsMyouKzWV1Ky7tOuuUqwEM7Ms5TIyMyV/BZGJoEigMt8/Znd2dnb29+zOzO7r+XjwAHaX3c8OszOveX/m8xmLIAgCiIiIiMj0YvRuABERERFpg8GOiIiIKEIw2BERERFFCAY7IiIiogjBYEdEREQUIRjsiIiIiCIEgx0RERFRhGCwIyIiIooQDHZEREREEYLBjoiIiChC6BrsDh06hIKCAvTu3RuJiYkYPnw4Nm7cKN0vCAIKCwuRkZGBxMRE5Obmorq6WscWExERERmXrsHujjvuwOrVq/HGG29gy5YtyMvLw+WXX459+/YBAGbNmoXZs2ejtLQUGzduRHp6OkaMGIFDhw7p2WwiIiIiQ7IIgiDo8cItLS1ITk7Gv//9b1x99dXS7VarFddccw3+9re/ISMjAwUFBXj44YcBAK2trUhLS8MzzzyDu+66S49mExERERlWnF4vfOzYMRw/fhwdO3Z0uj0xMRHr1q1DTU0N6urqkJeXJ92XkJCAnJwcrF+/3m2wa21tRWtrq/R7e3s7Dhw4gJNOOgkWiyU0b4aIiIhIQRAEHDp0CBkZGYiJCU8nqW7BLjk5GcOGDcPf/vY39O/fH2lpaXjzzTfxxRdfICsrC3V1dQCAtLQ0p79LS0vDrl273D5vcXExZsyYEdK2ExEREflqz5496NGjR1heS7dgBwBvvPEGJk6ciFNOOQWxsbE455xzMGbMGGzatEl6jLLKJgiCx8rb9OnTMXXqVOn3xsZG9OrVC3v27EFKSor2b4KIiIhIRVNTE3r27Ink5OSwvaauwe60005DZWUlmpub0dTUhO7du+Omm25CZmYm0tPTAQB1dXXo3r279Df19fUuVTy5hIQEJCQkuNyekpLCYEdERERhF85TwQwxj12nTp3QvXt3NDQ0YNWqVbjuuuukcLd69WrpcW1tbaisrMTw4cN1bC0RERGRMelasVu1ahUEQUDfvn3x448/4i9/+Qv69u2L2267DRaLBQUFBSgqKkJWVhaysrJQVFSEpKQkjBkzRs9mExERERmSrsGusbER06dPx969e9GlSxfceOONeOqppxAfHw8AmDZtGlpaWjB58mQ0NDRg6NChKC8vD2tfNREREZFZ6DaPXbg0NTUhNTUVjY2NPMeOiCgKHD9+HEePHtW7GRQlOnTo4HYqEz0yiK4VOyIiIq0IgoC6ujocPHhQ76ZQFImJiUFmZiY6dOigd1MAMNgREVGEsIe6bt26ISkpiZPSU8i1t7fj559/Rm1tLXr16mWIdY7BjoiITO/48eNSqDvppJP0bg5Fka5du+Lnn3/GsWPHpDECejLEdCdERETBsJ9Tl5SUpHNLKNrYu2CPHz+uc0tEDHZERBQxjNAVRtHFaOscgx0RERFRhGCwIyIiIgDA999/j/PPPx8dO3aE1WoN2esUFhaG9PmVr5WWlgaLxYJ3330XEyZMwPXXXy/dn5ubi4KCgrC0JRwY7IiIiHRUV1eH++67D3369EFCQgJ69uyJUaNG4eOPP3Z5bFFREWJjY/H000+73Ldw4UJYLBb079/f5b63334bFosFp556qse2PPHEE+jUqRO2bdum+vpaeeihh0L6/HZbt27FjBkz8Oqrr6K2thYjR47ECy+8gIULF4b8tfXCYEdERKSTn376CYMHD8Ynn3yCWbNmYcuWLVi5ciUuueQS3HvvvS6PX7BgAaZNm4bXX39d9fk6deqE+vp6fP755063v/766+jVq5fX9uzYsQMXXnghevfuHfDo4ra2Nq+P6dy5c1hGL+/YsQMAcN111yE9PR0JCQlITU3FCSecEPLX1guDHRERkU4mT54Mi8WCL7/8En/4wx9wxhlnIDs7G1OnTsWGDRucHltZWYmWlhbMnDkTzc3NWLt2rcvzxcXFYcyYMU7Bb+/evaioqPB6nXWLxYKvv/4aM2fOhMViQWFhIQBgy5YtuPTSS5GYmIiTTjoJd955Jw4fPiz9nb1rs7i4GBkZGTjjjDOk17355pvRpUsXdOrUCUOGDMEXX3wBwLUr1v4czz33HLp3746TTjoJ9957r9MVRGpra3H11VcjMTERmZmZWLp0KU499VSUlJSovp/CwkKMGjUKgDiJsH2Qg7IrVqmtrQ3Tpk3DKaecgk6dOmHo0KGoqKjwuOyMhPPYERER6eDAgQNYuXIlnnrqKXTq1MnlfmVVaf78+bjlllsQHx+PW265BfPnz8fFF1/s8ne33347Lr74YrzwwgtISkrCwoULceWVVyItLc1je2pra3H55ZfjyiuvxEMPPYTOnTvjyJEjuPLKK3H++edj48aNqK+vxx133IEpU6Y4dWd+/PHHSElJwerVqyEIAg4fPoycnByccsopeO+995Ceno5Nmzahvb3d7euvWbMG3bt3x5o1a/Djjz/ipptugtVqxaRJkwAA48aNw/79+1FRUYH4+HhMnToV9fX1bp/voYcewqmnnorbbrsNtbW1Ht+73G233YaffvoJy5YtQ0ZGBlasWIErr7wSW7ZsQVZWls/PoxcGOyLyqqqqCgUFBSgpKQnbCc9EQRsyBKirC//rpqcDX33l9WE//vgjBEFAv379vD62qakJ77zzDtavXw8AGDt2LC644AK89NJLLtcgtVqtOO200/Cvf/0Lt956KxYuXIjZs2dj586dXpqdjri4OHTu3Bnp6ekAgLlz56KlpQVlZWVS+CwtLcWoUaPwzDPPSGGxU6dOmDdvnjSn22uvvYZff/0VGzduRJcuXQAAp59+usfXP/HEE1FaWorY2Fj069cPV199NT7++GNMmjQJ33//PT766CNs3LgRQ4YMAQDMmzfPY9Dq3LmzFI7t78ebHTt24M0338TevXuRkZEBQAyIK1euxIIFC1BUVOTT8+iJwY6IvKqurkZlZSWqq6sZ7Mg86uqAffv0boVbgiAA8G0etKVLl6JPnz4YNGgQADG89enTB8uWLcOdd97p8viJEydiwYIF6NWrFw4fPoyrrroKpaWlfrdx69atGDRokFNF8YILLkB7ezu2bdsmBbuzzjrL6VqpVVVVOPvss6VQ54vs7GzExsZKv3fv3h1btmwBAGzbtg1xcXE455xzpPtPP/10nHjiiX6/J082bdoEQRCk7mS71tZW01zRhMGOiDzbsgUoLta7FUT+87FKo9frZmVlwWKxYOvWrR7P+QLEwQ/V1dWIi3Psttvb2zF//nzVYJefn49p06ahsLAQ48aNc/o7fwiC4DZ4ym9XdiUnJib6/VrKy3FZLBap69YegtXap6X29nbExsbi66+/dgqZgFgBNAMGOyJSVVVVherqauCuu1De3AwAKF+xQro/Ozub1TsyNh+6Q/XUpUsXXHHFFfjHP/6BP//5zy7h6ODBgzjhhBOwZcsWfPXVV6ioqHCqgB08eBAXX3wxvv32WwwYMMDlua+99lq8/fbbeOWVVwJu45lnnolFixahublZat9nn32GmJgYl6qW3MCBAzFv3jwcOHDAr6qdO/369cOxY8ewefNmDB48GIDYlX3w4MGgn1vu7LPPxvHjx1FfX4+LLrpI0+cOF46KJSJVBQUFGDt2LMY2N6PMdlvZO++It40dG1ETehLp5eWXX8bx48dx3nnn4Z133sH27duxdetWvPjiixg2bBgAcdDEeeedh4svvhgDBgyQvi688EIMGzYM8+fPV33uhQsXYv/+/T6dw+dOfn4+OnbsiPHjx+Pbb7/FmjVrcN999+HWW2/1OBjjlltuQXp6Oq6//np89tln2LlzJ9555x2XaVh81a9fP1x++eW488478eWXX2Lz5s248847kZiYqOklvc444wzk5+dj3LhxWL58OWpqarBx40Y888wz+OCDDzR7nVBisCMiVSUlJVi8eDEWAxhnu23cVVeJty1e7HaKASLyXWZmJjZt2oRLLrkEDz74IAYMGIARI0bg448/xpw5c9DW1obFixfjxhtvVP37G2+8EYsXL1adO84+PUkwkpKSsGrVKhw4cADnnnsu/vCHP+Cyyy7zer5ehw4dUF5ejm7duuGqq67CWWedhaefftqle9MfZWVlSEtLw8UXX4wbbrgBkyZNQnJyMjp27Bjwc6pZsGABxo0bhwcffBB9+/bFtddeiy+++AI9e/bU9HVCxSJo3UFtME1NTUhNTUVjY6PLyCEi8oHFgiUAxgJYPGMG8h9/XO8WEbn4/fffUVNTg8zMTM139GRMe/fuRc+ePfHRRx/hsssu060dntY9PTIIz7EjIt9F9nEgERnYJ598gsOHD+Oss85CbW0tpk2bhlNPPVV1Lr9oxmBHRF5lA8gBkJ2ZqXdTiChKHT16FI888gh27tyJ5ORkDB8+HEuWLHEZTRvtGOyIyCsrgAoAMMGs60QUma644gpcccUVejfD8Dh4goh8x65YIiJDY7AjIt8x2BERGRqDHRH5jsGOiMjQGOyIyHcMdkREhsZgR0S+Y7AjIjI0Bjsi8h2DHRGRoTHYEZHvGOyIIt7ChQtxwgknhPx1jhw5ghtvvBEpKSmwWCw4ePAgTj31VKfLFVosFrz77rshb0skYbAjIt8x2BFFFGWQAoCbbroJP/zwQ8hfe9GiRfj000+xfv161NbWIjU1FRs3bsSdd94Z8teOZAx2ROQ7BjuKArW1tSgsLERtba3eTXEhCAKOHTsW0tdITExEt27dQvoaALBjxw70798fAwYMQHp6OiwWC7p27YqkpKSQv3YkY7AjIt8x2FEUqK2txYwZM8IS7ARBwKxZs9CnTx8kJiZi0KBB+Ne//iXdX1FRAYvFglWrVmHIkCFISEjAp59+igkTJuD66693eq6CggLk5uZKv+fm5mLKlCmYMmUKTjjhBJx00kl47LHHINg+x7m5udi1axceeOABWCwWWCwWAOpdsXPmzMFpp52GDh06oG/fvnjjjTec7rdYLJg3bx5uuOEGJCUlISsrC++9957b952bm4u///3vWLt2LSwWi9RutQqi3L59+3DTTTfhxBNPxEknnYTrrrsOP/30k9vHRyMGOyLyHYMdkaYee+wxLFiwAHPmzEF1dTUeeOABjB07FpWVlU6PmzZtGoqLi7F161YMHDjQ5+dftGgR4uLi8MUXX+DFF1/E888/j3nz5gEAli9fjh49emDmzJmora11G2RXrFiB+++/Hw8++CC+/fZb3HXXXbjtttuwZs0ap8fNmDEDf/rTn/DNN9/gqquuQn5+Pg4cOKD6nMuXL8ekSZMwbNgw1NbWYvny5V7fy5EjR3DJJZegc+fOWLt2LdatW4fOnTvjyiuvRFtbm8/LJNLxWrFE5DsGO4pQ8mCzadMmp+8A0L17d3Tv3l3T12xubsbs2bPxySefYNiwYQCAPn36YN26dXj11VeRk5MjPXbmzJkYMWKE36/Rs2dPPP/887BYLOjbty+2bNmC559/HpMmTUKXLl0QGxuL5ORkpKenu32O5557DhMmTMDkyZMBAFOnTsWGDRvw3HPP4ZJLLpEeN2HCBNxyyy0AgKKiIrz00kv48ssvceWVV7o8Z5cuXZCUlIQOHTp4fG25ZcuWISYmBvPmzZOqiwsWLMAJJ5yAiooK5OXl+bxcIhkrdkTkOwY7ilCvvvoqBg8ejMGDB2PSpEkAgEmTJkm3vfrqq5q/5nfffYfff/8dI0aMQOfOnaWvsrIy7Nixw+mxQ4YMCeg1zj//fCkEAcCwYcOwfft2HD9+3Ofn2Lp1Ky644AKn2y644AJs3brV6TZ5JbFTp05ITk5GfX19QO1W8/XXX+PHH39EcnKytKy6dOmC33//3WV5RTNW7IjIdwx2FKHuuusuXHvttQDESt2kSZMwd+5cnHPOOQCgebUOANrb2wEA77//Pk455RSn+xISEpx+79Spk9PvMTEx0rlydkePHtW8jXbycAiI5wYqb4uPj3f5G/t71EJ7ezsGDx6MJUuWuNzXtWtXzV7H7BjsiIgo6ql1tZ5zzjlSsAuFM888EwkJCdi9e7dTt6svunbtim+//dbptqqqKpdwtWHDBpffs7KyEBsbCwDo0KGD1+pd//79sW7dOowbN066bf369ejfv79fbQ7WOeecg7feegvdunVDSkpKWF/bTHTtij127Bgee+wxZGZmIjExEX369MHMmTOdEr4gCCgsLERGRgYSExORm5uL6upqHVtNFMVYsSPSTHJyMh566CE88MADWLRoEXbs2IHNmzfjH//4BxYtWuTxby+99FJ89dVXKCsrw/bt2/HEE0+4BD0A2LNnD6ZOnYpt27bhzTffxEsvvYT7779fuv/UU0/F2rVrsW/fPuzfv1/1tf7yl79g4cKFeOWVV7B9+3bMnj0by5cvx0MPPRTcAvBTfn4+Tj75ZFx33XX49NNPUVNTg8rKStx///3Yu3dvWNtiZLoGu2eeeQavvPIKSktLsXXrVsyaNQvPPvssXnrpJekxs2bNwuzZs1FaWoqNGzciPT0dI0aMwKFDh3RsOVGUYrCjKNC9e3c88cQTIel+Vfrb3/6Gxx9/HMXFxejfvz+uuOIK/Oc//0FmZqbHv7viiivw17/+FdOmTcO5556LQ4cOOVXU7MaNG4eWlhacd955uPfee3Hfffc5TQA8c+ZM/PTTTzjttNPcdmdef/31eOGFF/Dss88iOzsbr776KhYsWOA0tUo4JCUlYe3atejVqxdGjx6N/v37Y+LEiWhpaWEFT8YiKDvpw+iaa65BWloa5s+fL9124403IikpCW+88QYEQUBGRgYKCgrw8MMPAwBaW1uRlpaGZ555BnfddZfX12hqakJqaioaGxv5jycKhPw8mvffB666Sr+2ELnx+++/o6amBpmZmejYsaPezTGE3NxcWK1Wj/PCUfA8rXt6ZBBdK3YXXnghPv74Y+nSJf/73/+wbt06XGXbcdTU1KCurs5pCHNCQgJycnKwfv161edsbW1FU1OT0xcRaUTDE6GJiEh7ug6eePjhh9HY2Ih+/fohNjYWx48fx1NPPSXNg1NXVwcASEtLc/q7tLQ07Nq1S/U5i4uLMWPGjNA2nIiIiMiAdA12b731FhYvXoylS5ciOzsbVVVVKCgoQEZGBsaPHy89zpdh1nbTp0/H1KlTpd+bmprQs2fP0LwBomjDc+yITKOiokLvJpAOdA12f/nLX/B///d/uPnmmwEAZ511Fnbt2oXi4mKMHz9emo26rq7O6STW+vp6lyqeXUJCgsv8P0SkEQY7IiJD0/UcuyNHjiAmxrkJsbGx0nQnmZmZSE9Px+rVq6X729raUFlZieHDh4e1rUQEBjsiIoPTtWI3atQoPPXUU+jVqxeys7OxefNmzJ49GxMnTgQgdsEWFBSgqKgIWVlZyMrKQlFREZKSkjBmzBg9m04UnRjsyOC0vNIBkS90nFxEla7B7qWXXsJf//pXTJ48GfX19cjIyMBdd92Fxx9/XHrMtGnT0NLSgsmTJ6OhoQFDhw5FeXk5kpOTdWw5UZQy2AaMyK5Dhw6IiYnBzz//jK5du6JDhw5uz8Um0oogCPj1119hsVhcrvqhF13nsQsHzmNHFCT5zvGdd4DRo/VrC5EHbW1tqK2txZEjR/RuCkURi8WCHj16oHPnzi736ZFBeK1YIvJdZB8Hksl16NABvXr1wrFjx7xe/5RIK/Hx8dK1d42AwY4oCPYpekpKSmC1WvVuTugx2JHB2bvEjNItRhRuuo6KJTK76upqVFZWorq6Wu+mhAeDHRGRoTHYEZHvGOyIiAyNXbFEfqqqqhIrdGvXonzRIgBAeXm5dH92dnbkdssy2BERGRpHxRL5KTc3F5WVlW7vz8nJiaxL+chHxb75JmC7UgwREXnGUbFEJlBSUiJW7MaORTmAMgDjxo1DXl4eALFiF7Ei+ziQiMj0GOyI/GS1WsWu1rFjAYjBLi8vD/n5+bq2KywY7IiIDI2DJ4jIdwx2RESGxmBHFIRsADmI8O5XOQY7IiJDY1csURCsACoAIFJHwSox2BERGRordkTkOwY7IiJDY7AjIt8x2BERGRqDHRH5jsGOiMjQGOyIyHcMdkREhsZgR0S+Y7AjIjI0Bjsi8h2DHRGRoTHYEZHvGOyIiAyNwY6IiIgoQjDYEZHvWLEjIjI0Bjsi8l17u94tICIiDxjsiMh3rNgRERkagx0R+Y7BjojI0BjsiIiIiCIEgx2RFqKlkhUt75OIyKQY7Ii0EC2BJ1reJxGRSTHYEWkhWgJPtLxPIiKTYrAj0kK0BB5Od0JEZGgMdkTku2gJsEREJsVgR6SFaAk80fI+iYhMisGOSAvREnjYFUtEZGgMdkRaiJZgFy3vk4jIpBjsiLQQLYGHFTsiIkNjsCPSQrQEu2h5n0REJsVgR0S+Y7AjIjI0BjsiLURL4GFXLBGRoeka7E499VRYLBaXr3vvvRcAIAgCCgsLkZGRgcTEROTm5qK6ulrPJhOpi5ZgFy3vk4jIpHQNdhs3bkRtba30tXr1agDAH//4RwDArFmzMHv2bJSWlmLjxo1IT0/HiBEjcOjQIT2bTeQqWgIPK3ZERIama7Dr2rUr0tPTpa///ve/OO2005CTkwNBEFBSUoJHH30Uo0ePxoABA7Bo0SIcOXIES5cu1bPZRK6iJdhFy/skIjIpw5xj19bWhsWLF2PixImwWCyoqalBXV0d8vLypMckJCQgJycH69evd/s8ra2taGpqcvoiCrloCTzR8j6JiEzKMMHu3XffxcGDBzFhwgQAQF1dHQAgLS3N6XFpaWnSfWqKi4uRmpoqffXs2TNkbSaKOuyKJSIyNMMEu/nz52PkyJHIyMhwut1isTj9LgiCy21y06dPR2Njo/S1Z8+ekLSXyEm0VLKi5X0SEZlUnN4NAIBdu3bho48+wvLly6Xb0tPTAYiVu+7du0u319fXu1Tx5BISEpCQkBC6xhKpiZbAw4odEZGhGaJit2DBAnTr1g1XX321dFtmZibS09OlkbKAeB5eZWUlhg8frkczidyLlmAXLe+TiMikdK/Ytbe3Y8GCBRg/fjzi4hzNsVgsKCgoQFFREbKyspCVlYWioiIkJSVhzJgxOraYSEW0BJ5oeZ9ERCale7D76KOPsHv3bkycONHlvmnTpqGlpQWTJ09GQ0MDhg4divLyciQnJ+vQUiIPoiXwsCuWiMjQLIIQ2XukpqYmpKamorGxESkpKXo3hyKJfBDPr78CJ5+sX1tCSf4+H34YePpp/dpCRGQiemQQQ5xjR0QmwYodEZGhMdgRaSGyC98O0fI+iYhMisGOSAvREnii5X0SEZkUgx2RFqIl8LArlojI0BjsiLTAYEdERAbAYEekBQY7IiIyAAY7IvLd8eN6t4CIiDxgsCPSAit2RERkAAx2RFqIlmDHih0RkaEx2BFpIVqCHSt2RESGxmBHpIVIDnYxss0Egx0RkaEx2BFpIVqCHbtiiYgMjcGOiDxjxY6IyDQY7Ii0wIodEREZAIMdkRaiJdixYkdEZGgMdkRaiORgZ7E4fmawIyIyNAY7Ii1EcrBjVywRkWkw2BFpIZKDXWys42dW7IiIDI3BjkgLkRzsWLEjIjINBjsi8oyDJ4iITIPBjkgLkVyx4+AJIiLTYLAj0kIkBzt2xRIRmQaDHZEWoiXYsWJHRGRoDHZEWoiWYMeKHRGRoTHYEWkhkoMdpzshIjINBjsi8kw+eIIVOyIiQ2OwI9JCJFfseI4dEZFpMNgRaYHBjoiIDIDBjkgLkRzs5NgVS0RkaAx2RFqI5GAnf2+s2BERGRqDHZEWIjnYybFiR0RkaAx2ROQZK3ZERKbBYEekhWip2DHYEREZGoMdkRaiJdixK5aIyNAY7Ii0EMnBjl2xRESmoXuw27dvH8aOHYuTTjoJSUlJsFqt+Prrr6X7BUFAYWEhMjIykJiYiNzcXFRXV+vYYiIVkRzs5FixIyIyNF2DXUNDAy644ALEx8fjww8/xHfffYe///3vOOGEE6THzJo1C7Nnz0ZpaSk2btyI9PR0jBgxAocOHdKv4URKJgx2VVVVyM3NRVVVlecHsmJHRGQacXq++DPPPIOePXtiwYIF0m2nnnqq9LMgCCgpKcGjjz6K0aNHAwAWLVqEtLQ0LF26FHfddVe4m0wUMaqrq1FZWYnq6mpYrVbf/ojBjojI0HSt2L333nsYMmQI/vjHP6Jbt244++yzMXfuXOn+mpoa1NXVIS8vT7otISEBOTk5WL9+vR5NJlJntopdYyOwYoVvj5W/N3bFEhEZmq4Vu507d2LOnDmYOnUqHnnkEXz55Zf485//jISEBIwbNw51dXUAgLS0NKe/S0tLw65du1Sfs7W1Fa2trdLvTU1NoXsDRHYmCXZVVVXiOarz56N8zRoAQPm//iXdn52d7bl6x4odEZGh6Rrs2tvbMWTIEBQVFQEAzj77bFRXV2POnDkYN26c9DiLxeL0d4IguNxmV1xcjBkzZoSu0URqTBLsCgoKUFlZ6XRb2bvvouzddwEAOTk5qKiocP8ErNgRERmarl2x3bt3x5lnnul0W//+/bF7924AQHp6OgBIlTu7+vp6lyqe3fTp09HY2Ch97dmzJwQtJ1IwSbArKSnB4sWLsRiA/dBp3Pnni7ctXoySkhLXP+LgCSIi09C1YnfBBRdg27ZtTrf98MMP6N27NwAgMzMT6enpWL16Nc4++2wAQFtbGyorK/HMM8+oPmdCQgISEhJC23AiJZMEO6vVKna1jh0LACgDkHfmmcjPz/ftCRjsiIgMTddg98ADD2D48OEoKirCn/70J3z55Zd47bXX8NprrwEQu2ALCgpQVFSErKwsZGVloaioCElJSRgzZoyeTSdyZpJgp8pbWOPgCSIi09A12J177rlYsWIFpk+fjpkzZyIzMxMlJSVO1YNp06ahpaUFkydPRkNDA4YOHYry8nIkJyfr2HIi88sGkAMg23bKg09YsSMiMjSLIJi51OBdU1MTUlNT0djYiJSUFL2bQ5FEPoDns8+A4cP1a4u/5G0vLQXuvdf9Y3v0APbtE3/u3Bng5OBERD7RI4PofkkxIlNSHg+Z+fjInyocK3ZERIbGYEekBTMHO3/azmBHRGRoDHZEgYikip23tnPwBBGRaTDYEWnBzMGOXbFERBGDwY4oEGYOckqs2BERRQwGOyItmDno+VuFM/N7JSKKcAx2RIGI1nPsAFbtiIgMjMGOSAuRHOyUeJ4dEZFhMdgRacHMwc7foMZgR0RkWAx2RIFgVywRERkQgx1RtGPFjogoYjDYEQWCFTsiIjIgBjsiLURysFNixY6IyLAY7IgCEUkVO29BTfneGOyIiAyLwY5IC2YOdv62nV2xRESGxWBHFIhoqtgF+3giIgobBjsiLZg52HHwBBFRxGCwIwqEmYOcEgdPEBFFDAY7Ii2YOehx8AQRUcRgsCMKRCSdY8fBE0REEYPBjkgLZg52rNgREUUMBjuiQLBiR0REBsRgR6SFaAp2rNgRERkWgx1RIMwc5JTYFUtEFDEY7Ii0YOagx65YIqKIwWBHFIhoOseOFTsiItNgsCPSgtmCncXi+NnfoMaKHRGRYTHYEQXC7BW7GNlHnxU7IqKIwWBHpAWzBbtgKnYMdkREhsVgRxSN5MGOgyeIiCIGgx1RINgVS0REBsRgR6QFswU7Dp4gIopIQQW71tZWrdpBZC6s2BERkQH5FexWrVqFCRMm4LTTTkN8fDySkpKQnJyMnJwcPPXUU/j5559D1U4iYzNbsGPFjogoIvkU7N5991307dsX48ePR0xMDP7yl79g+fLlWLVqFebPn4+cnBx89NFH6NOnD+6++278+uuvoW43kb7MXrHzZ/AEK3ZERKYR58uDioqK8Nxzz+Hqq69GTIxrFvzTn/4EANi3bx9eeOEFlJWV4cEHH9S2pURGZrZg509XrBKDHRGRYflUsfvyyy8xatQo1VAnd8opp2DWrFk+h7rCwkJYLBanr/T0dOl+QRBQWFiIjIwMJCYmIjc3F9XV1T49N1FImS3IKbErlogoIuk+KjY7Oxu1tbXS15YtW6T7Zs2ahdmzZ6O0tBQbN25Eeno6RowYgUOHDunYYiIVZgt6HDxBRBSRfOqKlZs6darq7RaLBR07dsTpp5+O6667Dl26dPGtAXFxTlU6O0EQUFJSgkcffRSjR48GACxatAhpaWlYunQp7rrrLn+bTqSdSDrHjhU7IqKI4Xew27x5MzZt2oTjx4+jb9++EAQB27dvR2xsLPr164eXX34ZDz74INatW4czzzzT6/Nt374dGRkZSEhIwNChQ1FUVIQ+ffqgpqYGdXV1yMvLkx6bkJCAnJwcrF+/nsGOjMXMwY4VOyKiiOF3V+x1112Hyy+/HD///DO+/vprbNq0Cfv27cOIESNwyy23YN++fbj44ovxwAMPeH2uoUOHoqysDKtWrcLcuXNRV1eH4cOH47fffkNdXR0AIC0tzelv0tLSpPvUtLa2oqmpyemLSHNmr9hx8AQRUUTyu2L37LPPYvXq1UhJSZFuS0lJQWFhIfLy8nD//ffj8ccfd6q0uTNy5Ejp57POOgvDhg3DaaedhkWLFuH8888HIHbxygmC4HKbXHFxMWbMmOHv2yIKjtmCnT9dscr3xq5YIiLD8rti19jYiPr6epfbf/31V6k6dsIJJ6Ctrc3vxnTq1AlnnXUWtm/fLp13p6zO1dfXu1Tx5KZPn47Gxkbpa8+ePX63g8grswU5JX+6YpVYsSMiMqyAumInTpyIFStWYO/evdi3bx9WrFiB22+/Hddffz0AcXqUM844w+/GtLa2YuvWrejevTsyMzORnp6O1atXS/e3tbWhsrISw4cPd/scCQkJSElJcfoiCjmzBb1ggh0rdkREhuV3V+yrr76KBx54ADfffDOOHTsmPklcHMaPH4/nn38eANCvXz/MmzfP63M99NBDGDVqFHr16oX6+no8+eSTaGpqwvjx42GxWFBQUICioiJkZWUhKysLRUVFSEpKwpgxY/xtNpG2zH6OXTBdsazYEREZlt/BrnPnzpg7dy6ef/557Ny5E4Ig4LTTTkPnzp2lx1itVp+ea+/evbjllluwf/9+dO3aFeeffz42bNiA3r17AwCmTZuGlpYWTJ48GQ0NDRg6dCjKy8uRnJzsb7OJQstswS6YwROs2BERGZbPwS42Nha1tbXo1q0bADHgDRw4MKgXX7Zsmcf7LRYLCgsLUVhYGNTrEGmOFTsiIjIgn8+xE8y24yIKJ7N9PniOHRFRRNL9kmJEpmS2IKfES4oREUUkv86xW7VqFVJTUz0+5tprrw2qQUSmZLagx0uKERFFJL+C3fjx4z3eb7FYcJwbfYpGZgt2rNgREUUkv7pi6+rq0N7e7vaLoY6iRiQNnuAExUREEcPnYOfpMl5EUc9swU5esWNXLBFRxOCoWKJAmL17MjbW8bO3oGb290pEFEV8Dnbjx49HYmJiKNtCZF5mO/DxJ9gpsWJHRGRYPg2eaG5uxoIFC3x+0ubmZnTq1CngRhEZntmrWKzYERFFJJ8qdqeffjqKiorw888/u32MIAhYvXo1Ro4ciRdffFGzBhKZAit2RERkAD5V7CoqKvDYY49hxowZsFqtGDJkCDIyMtCxY0c0NDTgu+++w+eff474+HhMnz4dd955Z6jbTaQvs4+KlQc7XlKMiChi+BTs+vbti3/+85/Yu3cv/vnPf2Lt2rVYv349WlpacPLJJ+Pss8/G3LlzcdVVVyEmhhezoChktrAj/5yyYkdEFDH8mqC4R48eeOCBB/DAAw+Eqj1E5hBJFTt/g5rZQiwRURRheY1IC2YLO8F0xbJiR0RkWAx2RIEwe8UumK5Ys4VYIqIowmBHpAWzBTv5lWQ43QkRUcRgsCMKRCSFHQ6eICKKGAx2RFowW8VO3l52xRIRRQyfRsV+8803Pj/hwIEDA24MkWmYvWInb7+/bWfFjojIsHwKdlarFRaLBYIgwCI/N0fFcW70KRqZrWInx4odEVHE8KkrtqamBjt37kRNTQ3eeecdZGZm4uWXX8bmzZuxefNmvPzyyzjttNPwzjvvhLq9RMZg9lGxvnbFqr0vHrwRERmWTxW73r17Sz//8Y9/xIsvvoirrrpKum3gwIHo2bMn/vrXv+L666/XvJFEhmfmKhYrdkREEcPvwRNbtmxBZmamy+2ZmZn47rvvNGkUkeGxYkdERAbkd7Dr378/nnzySfz+++/Sba2trXjyySfRv39/TRtHZBpmq2JxVCwRUUTy61qxAPDKK69g1KhR6NmzJwYNGgQA+N///geLxYL//ve/mjeQyBTMVrGT87dix2BHRGRYfge78847DzU1NVi8eDG+//57CIKAm266CWPGjEGnTp1C0UYi44mWrlg17IolIjIsv4MdACQlJeHOO+/Uui1E5mXmKha7YomIIkZAwQ4AvvvuO+zevRttbW1Ot1977bVBN4rI8CKpYucpqHHwBBGRqfgd7Hbu3IkbbrgBW7ZskSYtBiBNXMwJiikqma2KxcETREQRye9Rsffffz8yMzPxyy+/ICkpCdXV1Vi7di2GDBmCioqKEDSRyIDMXrGTO3bM/X2s2BERmYrfFbvPP/8cn3zyCbp27YqYmBjExMTgwgsvRHFxMf785z9j8+bNoWgnkbGZrYrFih0RUUTyu2J3/PhxdO7cGQBw8skn4+effwYgXp1i27Zt2raOyKgiqWLX3u6+/azYERGZit8VuwEDBuCbb75Bnz59MHToUMyaNQsdOnTAa6+9hj59+oSijUTGZ7Zgp2zv8eNAnI+bA1bsiIgMy+9g99hjj6G5uRkA8OSTT+Kaa67BRRddhJNOOglvvfWW5g0kMiRlMDJb2FG2/9gx34MdK3ZERIbld7C74oorpJ/79OmD7777DgcOHMCJJ54ojYwlijpmq9gpuQtrvPIEEZGpBDyPHQDs3bsXFosFp5xyilbtITKHSKvY+VOFM9t7JSKKIn4Pnmhvb8fMmTORmpqK3r17o1evXjjhhBPwt7/9De1BbPCLi4thsVhQUFAg3SYIAgoLC5GRkYHExETk5uaiuro64NcgChmzVezUumJ9eRzArlgiIgPzO9g9+uijKC0txdNPP43Nmzdj06ZNKCoqwksvvYS//vWvATVi48aNeO211zBw4ECn22fNmoXZs2ejtLQUGzduRHp6OkaMGIFDhw4F9DpEmomkUbEAK3ZERBHC72C3aNEizJs3D/fccw8GDhyIQYMGYfLkyZg7dy4WLlzodwMOHz6M/Px8zJ07FyeeeKJ0uyAIKCkpwaOPPorRo0djwIABWLRoEY4cOYKlS5f6/TpEIWW2sMOKHRFRRPI72B04cAD9+vVzub1fv344cOCA3w249957cfXVV+Pyyy93ur2mpgZ1dXXIy8uTbktISEBOTg7Wr1/v9vlaW1vR1NTk9EWkOVbsiIjIgPwOdoMGDUJpaanL7aWlpRg0aJBfz7Vs2TJs2rQJxcXFLvfV1dUBANLS0pxuT0tLk+5TU1xcjNTUVOmrZ8+efrWJKCBmCzu+VuzUsGJHRGRYfo+KnTVrFq6++mp89NFHGDZsGCwWC9avX489e/bggw8+8Pl59uzZg/vvvx/l5eXo2LGj28cpp1ARBMHjtCrTp0/H1KlTpd+bmpoY7kh7Zq/Y+ToqltOdEBGZit8Vu5ycHPzwww+44YYbcPDgQRw4cACjR4/Gtm3bcNFFF/n8PF9//TXq6+sxePBgxMXFIS4uDpWVlXjxxRcRFxcnVeqU1bn6+nqXKp5cQkICUlJSnL6IQs5swU7JnyocK3ZERIYV0Dx2GRkZeOqpp5xu27NnDyZOnIjXX3/dp+e47LLLsGXLFqfbbrvtNvTr1w8PP/ww+vTpg/T0dKxevRpnn302AKCtrQ2VlZV45plnAmk2kXYibR47fwZPmO29EhFFkaAmKJY7cOAAFi1a5HOwS05OxoABA5xu69SpE0466STp9oKCAhQVFSErKwtZWVkoKipCUlISxowZo1WzibTBih0RERmAZsEuFKZNm4aWlhZMnjwZDQ0NGDp0KMrLy5GcnKx304icma2KxYodEVFEMlSwq6iocPrdYrGgsLAQhYWFurSHyK1oGTyhhsGOiMiw/B48QUQqzBbslDjdCRFRRPC5Yjd69GiP9x88eDDYthCZR6QNnuB0J0REEcHnYJeamur1/nHjxgXdICJTMnvFjoMniIgigs/BbsGCBaFsB5G5RFrFjoMniIgiAs+xI9KC2Sp2wQyeYMWOiMiwGOyIAmH2ip0SK3ZERBGBwY5IC6zYERGRATDYEQUi0uax82e6E1bsiIgMi8GOSAtmDzuc7oSIKCIw2BEFItIqduyKJSKKCAx2RFowexWLgyeIiCICgx1RIFixIyIiA2KwI9KC2YMdK3ZERBGBwY4oEJE2jx0rdkREEYHBjkgLkVqxU2P2EEtEFMEY7IgCES0VO7XAyoodEZFhMdgRaSGaKnZqf09ERIbAYEcUiGgZFevufbFqR0RkSAx2RFowe1esvxU7s79fIqIIxWBHpIVoq9gx2BERGRKDHVEgomXwhFaPJyKisGCwI9KC2St27IolIooIDHZEgeDgCW3bE6Cqqirk5uaiqqpK76YQERkCgx2RFsxewTJpxa66uhqVlZWorq7WuylERIbAYEcUCFbstG1PoFpa9G4BEUU4s/UMxOndAKKIYJAKVsAMXLGrqqpCQUEBSkpKYLVaUVVVJVboampQXlgIACj/4APp8dnZ2bBarWFrHxFFNnnPgBm2LQx2RIFgxU7b9nig3KgWFBSgsrLS6TFlS5eibOlSAEBOTg4qKirC1j4iinAmO3BnsCPSgsk++KYaFdva6vRrSUmJWLEbOxblAMoAjDvzTOQ98ggAsWJHRBQMqWdgzRqUL1oEACgvL5fuN3LPAIMdUSDMXrFTMljFTtqofvwxyhcuBKCyUbX9XAYgr3dv5Ofnh7RNRBQ9VHsGyspQVlYGwNg9Awx2RFowW7AzeMXOp42q/E6zLX8iMjTVnoFx45CXlwfA2D0DDHZEgTD7lSd8PcfOnRC/X9WN6q23Iu+KKwDYNqqDByO7vR05ALJPOCGk7SGi6GK1WsWu1rFjAdh6BvLyTNEzwGBHpAWzV4zcVex06opV3ahedpnzRjU2Ftb2drFy16VLSNtDRGSW7TznsSMKBCt22rXFV8o2xsg2X2Zb/kRkGtmA2DPQt6/eTfEJK3ZEWjDJkZxbBqvYyUkb1TPOcL4jNlaX9hBRdLECYs+Agc+rk2OwIwqE2UfFmqhiZ4Vto3rmmc53yIMdK3ZEFGomOYBkVyyRFswWLHwdFWuAip1E2UZ5V6xJNrhEZGL+zh6gE12D3Zw5czBw4ECkpKQgJSUFw4YNw4cffijdLwgCCgsLkZGRgcTEROTm5vJi32QMZq/YKRm4YifxdI4dgx0RhRqDnXc9evTA008/ja+++gpfffUVLr30Ulx33XVSeJs1axZmz56N0tJSbNy4Eenp6RgxYgQOHTqkZ7OJXEVqxc4do1XszLb8ich8THIAqWuwGzVqFK666iqcccYZOOOMM/DUU0+hc+fO2LBhAwRBQElJCR599FGMHj0aAwYMwKJFi3DkyBEstV0Tkkg30VKxc/e+9AhSymDHc+yIKJxYsfPP8ePHsWzZMjQ3N2PYsGGoqalBXV2dNMszACQkJCAnJwfr1693+zytra1oampy+iIKObMFC4NfeUIVgx0R6YnBzjdbtmxB586dkZCQgLvvvhsrVqzAmWeeibq6OgBAWlqa0+PT0tKk+9QUFxcjNTVV+urZs2dI208EwHzBwtdRsUYaPMFz7IhITybZzuge7Pr27Yuqqips2LAB99xzD8aPH4/vvvtOut9isTg9XhAEl9vkpk+fjsbGRulrz549IWs7RTGzT1CsZPaKnUk2uERkYiap2Ok+j12HDh1w+umnAwCGDBmCjRs34oUXXsDDDz8MAKirq0P37t2lx9fX17tU8eQSEhKQkJAQ2kYTKZktWJixYseuWCLSk0mCne4VOyVBENDa2orMzEykp6dj9erV0n1tbW2orKzE8OHDdWwhEYKf4NdozFCxY1csEenJJMFO14rdI488gpEjR6Jnz544dOgQli1bhoqKCqxcuRIWiwUFBQUoKipCVlYWsrKyUFRUhKSkJIwZM0bPZhO5MlvFKNhgarSKHYMdEYWaSbYzuga7X375Bbfeeitqa2uRmpqKgQMHYuXKlRgxYgQAYNq0aWhpacHkyZPR0NCAoUOHory8HMnJyXo2m8j8FbtgrzxhhHPsOI8dEYUTK3bezZ8/3+P9FosFhYWFKCwsDE+DiAJl9mBhhoqdp65Ysy9/IjI+kwQ7w51jR2QKrNhp2x5fcPAEEenJJNt5BjsiLZjkAy+JhHPsOHgipKqqqpCbm4uqqiq9m0JkDKzYEUWwaJnHziwVOwY7zVVXV6OyslK6djdR1GOwI4oiZgsWwVbsjDDdCbtiQ6u1Ve8WEBmLSbbzuk9QTGRK0VKxc4ddsRGpqqpKrNC1tKB8yhQAQPmSJdL92dnZsFqtOrWOSGcmqdgx2BFpwWzBItgrT7ArNiIVFBSgsrLS6bayDz9E2YcfAgBycnJQUVGhQ8uIDIDBjiiCRcuoWHc43UlEKikpESt2S5ag/MMPUQZgHIC8xYsBiBU7oqhlku08gx2RFgRB/LJY9G5JYMxesWOw04TVahW7WrdsAWzBLg9Afn6+zi0jMgCTVOw4eIIoEGqBx0zhwowVO155InzkoZmIRAx2RFHGzOHCDBU7T12xJukiMY3YWGQDyAHAzlciG5NsZ9gVSxQItcBz/DgQHx/+tgQiEip2HDwROrGxsAKo0LkZRIbCih1RlDFTuFAGO0HwrwpntHPszLTszYBdsUSuGOyIIpjZz7FToxaO3HXFGqFiJ6+OHj0a3rZEOgY7IlcMdkRRxkxVI3ddyZ7oPQrVU7BrawtvWyIdgx2RK5McQDLYEWnF7BU7taNReQDUO9gpN6ryqWVMssE1DQY7IlcmOYBksCMKRCAVLyMJpP16j0JVXrtU/h4Y7LTFYEfkisGOKMqYPdgZvWKn3Kgy2IUOgx2RK+XBpUEx2BEFIloGT8jpPQpVuVGVL28GO20x2BG5YsWOKMpEYsVOjhW76MFgR+SKwY4ogpn9HDs13qY70bti5ynYmb1aajRxnLueyAW7YomijJnChRkrdp4GT5C2WLEjcsWKHVEEM3vFztf2m6ViR9qK4a6ByAWDHVGUMVPFTo3RK3begh2Dnna4LIlcsSuWKIKxYqd9m7zx1hVrpuVvdAx2RK6fA1bsiKKMmYJFsOfYGbEr1iQbXVMwe/WZSAsm3cYw2BEFIhrnsZOPlNTjYtjeKnYm6SYxBVbsiEy7jWGwI9KKmSp2arxdeSI+3vFzuN6r/HqwyqNlZZA2ydG0KTDYEbFiRxRVWLELbVvs5KMz2RUbPmZfl4m0YNJtDIMdkVbMUrFzV43xVrGTB7twvVf563vrijXJRtcUWLEjYlcsUVQx86hYdzttf64VG66Knbyt3ip2JtnomgKDHZFpDx4Z7Ii0YvbuK38qdkYMdibZ6JqC2ddlIi2YdBvDYEcUiGis2IV78IS3ipxJN7qmwIodkWl7BRjsiLRi9iqHtypcuCt23oIbg13oMNgRmXbkPYMdUSCipWKn5+AJf4OdSY6mTUG5QzP7QQtRIEx68MhgR6QVswc7o1fs2BUbPgzNRKb9HOga7IqLi3HuueciOTkZ3bp1w/XXX49t27Y5PUYQBBQWFiIjIwOJiYnIzc1FdXW1Ti0m8sDsVQ1/KnZ6dcXKbzNpN4ldVVUVcnNzUVVVpXdTXCmX/e+/69MOIj2Z9OBR12BXWVmJe++9Fxs2bMDq1atx7Ngx5OXlobm5WXrMrFmzMHv2bJSWlmLjxo1IT0/HiBEjcOjQIR1bTlEvErti/anY6dEVKwjObTTpRteuuroalZWVxjxQVYZmk1QqiDSl3MYcPWqK80/jvD8kdFauXOn0+4IFC9CtWzd8/fXXuPjiiyEIAkpKSvDoo49i9OjRAIBFixYhLS0NS5cuxV133aVHs4nUsWKnLbXl2dbmGJ1r0m4SyTffiN/1uO6uN6zYEamHuLY2ICEh/G3xg67BTqmxsREA0KVLFwBATU0N6urqkJeXJz0mISEBOTk5WL9+PYMd6YcVO+3a4467jWqnTur3m6BiV1VVJVbofv4Z5bNmAQDKX3xRWrbZ2dmwWq06ttDG7KGZSAsMdsERBAFTp07FhRdeiAEDBgAA6urqAABpaWlOj01LS8OuXbtUn6e1tRWtso1QU1NTiFpMpGD2YGe0a8WqtVMeMEwY7AoKClBZWel0W9mmTSgbOxYAkJOTg4qKCh1apqCslrJiR9HIXbAzOMMEuylTpuCbb77BunXrXO6zWCxOvwuC4HKbXXFxMWbMmBGSNhJJ1D7wZu+KNdqVJ7xtVE0Y7EpKSsSK3fLlKF++HGUAxgHIW7wYgFixMwRW7IjUt+km+CwYItjdd999eO+997B27Vr06NFDuj09PR2AWLnr3r27dHt9fb1LFc9u+vTpmDp1qvR7U1MTevbsGaKWE8lEesVOfq1YPbti3d1vgg2u1WoVu1p37QJswS4PQH5+vs4tU+A5dkSmrdjpOipWEARMmTIFy5cvxyeffILMzEyn+zMzM5Geno7Vq1dLt7W1taGyshLDhw9Xfc6EhASkpKQ4fRFpzswVO3/OsZM/NiYGsFfKjdgVa4JgJ4kzxDG1exwVSxox9LQ+3jDY+e/ee+/F4sWLsXTpUiQnJ6Ourg51dXVoaWkBIHbBFhQUoKioCCtWrMC3336LCRMmICkpCWPGjNGz6USuzFKxc8db+y0WRyAxQsVOGT6OHg1te7QUF4dsADkADNL56owVO9KIoaf18cbbwaVB6XrYOGfOHABAbm6u0+0LFizAhAkTAADTpk1DS0sLJk+ejIaGBgwdOhTl5eVITk4Oc2uJZKJlVKzysXFxYoBixS44cXGwAqjQuRlumXnZknHs3g18+KHerQicSSt2ugY7wYeJ/iwWCwoLC1FYWBj6BhEFwyxdse74UrGzn2en5zx2diYcPCExW1csK3bkB2lan4ceQrltdovy8nLpfsNM6+MNK3ZEUSSaK3aAfl2xJp/uRGL0YMeKHQVBdVqfsjKUlZUBMNC0Pt6wYkcUOlVVVSgoKEBJSYlxj/TMUrELdFRsuCt2ETgqVsJgRxFMmtZn7FiUA+K0PuPGSRcbMMy0Pt6wYkcUOvITcA0R7MxcsXPHl7DGip02jB7s2BVLQZCm9bFNvF0GIC8vz3jT+njj7XQQg9J1VCyRTw4cAMwwVN4swc7dpMPerhUrf7wRK3Ym2OBKjB7sWLEjYlcskdakE3AfewzlP/0EwEAn4Jp5Hju5uDhHSPMW1ozQFRspo2LdXDnHMDjdCWlEmtanXz+9m+I/dsUSact0J+BGU8XOCPPYmbli58OMALriBMWkESts0/qceaau7QgIK3ZE2iopKUH1N98A48cb7wRctQ98OKpYWvAn2MkZYboTecAwc/gwenWXFTvSWlsbkJiodyv8Y9KKHc+xI8OyWq3Iv/lm5EO8nibgOAE3Pz/fGIMo5MxSsZOTBzszTHfCil14mLmbW4WpL2sVKcy4Dpm0YsdgR8Zm1MqGmUfFBlOx03vwRKSMijXqem0XYaNiTX1Zq0hhps+nHSt2FE5RcwRq28FIJ+Aaef4js3TFynmr2CnpPXgiUuaxM3qwM/OyVfPLL3q3gMy4Dpm0Ysdz7EzKcPO6hYptB2iF7QRcI79XswQ7+cYqPt7xsxm6YlmxC48IOMdOGlW/dy/K/+//AADlb78t3W+ay1pFCjN9Pu28nedrUAx2ZiQIQEOD3q0ID6PuACNl8ESHDo6fjx71/HfywRPHj4vPE8ppOyL5HDujrtd2EdAVqzqq/r33UPbeewAMOKo+0pnp82nHih2FmnQEWlqK8g0bABhoXrdQMfoOUM4swU5OXrFTC3buKnaAGO5COdFuJM9jx8ETISdd1qqwEOU//iiOqj/5ZOSVlAAw+GkdkciE65BZz7FjsDMR083rpgVlsGtvB2IMcGpotFbs9A52kVqxO3rUOWTrLQKCnXRZq5dfBmzBLi8lxXyXtYoUZvp82rFiR6EWMRdW9odyB9jWBnTsqE9bvDFLsJOThwm1DZZyw2bvigXE95uQEJp2Ad7Pb4mkYNfSYqxgFwFdsRL5gaCZegAijQkPDhjsKOSsViusgwaZ/8LK/lCbhNYIwY4Vu9APoPC2UY2kCYqPHAFSUvRpi5oIGDwhsVgco+qNsO2IFmY+8LJjVyyFhdHPzdGaWsXOqMwY7OQVN2/BDnCt2IVSNI2KPXJEn3a4EwFdsZKYGMeoerNd+cDMzLTtdsekFTsDnKxEfrFVSUwxr5sWlFUho+xgzFyxkwtm8IQewc7TOXZHj5qnq83owY5dsRQsM1fU7Vixo7Aw07xuWjDTxsEswc7dPHZm6Ir1VLEDxPcQyvP+tKJse3OzPu1wJ5K6Yhns9KF2HqnZqK0vrNiR5qJtw2TUYGfmip287fKw5kvFTu+uWE8VO+X9Rmb0il0kdcXK51qMtu2nnpTL2mgHL74wacWOwc5szHI9Uq2Y6TwNswQ7OYvFUbWLhIqdCTa6AIwf7CK1KzbazlHWU6QGOyPvg2wY7Mwm2o44WbHTnrLtnoIdK3ahYfRgZ+bzF5XYFauPSA12RtkHecBgZzbRtmEy+g5QzizBTk5esfMWipQVu1C/X2/ntzDYhY5Jd2iq2BWrj0gNdibYxjDYmU20d8UaZeMQSRU7+1x2vkx3wq5YbRh1vbZTC0Bm7Y5lxU4fRl/HfcFgR2ERbRsm5fs9fFifdvjCjMHO2zl2Ru6KNemINQCu780MFbtICHY8xy58IjXYmeDgkcHObJQflkjfUBl142Dmip1cpA2eMEuwY1ds+LBipw+jbrv9YdJtDIOd2Sh3ppHeNWumip1Z/heRPHjCLOHD6MEukrpieY6dPiIh2Hm7XrVBMdiZjXJFM2OVyB9GDXas2BnvkmKAKY6mARg/2LErloIVCcHOpNsYBjuz0TjYVVVVITc3F1VVVUE9T8iYaeNglmAXTMVO767YSB0Va7T12szVUCV2xerD6Ou4L0z6OWCwMxvlzjTIMFFdXY3KykpUV1cH9Twho3y/rNgFz93giWPHPFc0LJbwdsWq7YTl86mZdKMLwPiDJyKpK5bBTh+RGuxMcPDIYGc2WlbsBAHYsye49oSaUbti1Zgl2MlZLI7pTgDv3bF6V+wARxtNutEFwK7YcOI5dvqI1GB37Jjh16M47w8hQ9Eg2FVVVYkVujffRPn77wMAysvLpfuzs7NhtVqDaaV2zLRxMEuwc9cVC4ihSR70PHXF+jLvXTDcBbvWViAhwdwVO6Ov12Zetko8x04fymVttHXcF+7Wl7Y2oGPH8LbFDwx2ZqNBV2xBQQEqKyudbisrK0NZWRkAICcnBxUVFYG2UFtGrdiZuStWTt4VC3gOa/48VgueNqru7mfFThvsiqVgGf3gxReeDi4Z7EgzGlTsSkpKxIrd2LEoB1AGYNyttyLviisAiBU7wzBqsFNjlmDnrWIX6GO15mmjCph7gmKjB7tIrdgx2IWPclm3tYnbyDgTxQ5vB5cGZaIlTAA0CXZWq1Xsah07FoAY7PIuuwz5+fnBt09rRg12Zq7YuRs8AXiv2PlzPl6wAqnYmSV8mDHYmbVix3Ps9KG2rJubgdTU8LclUO7WF4MHOw6eMBtlV6xWO1ejhhIzlfONugw98RbsPFXsQr1x81axM3NXrNHPP4rUrlieYxc+7oKdmXjbBhmUrsFu7dq1GDVqFDIyMmCxWPDuu+863S8IAgoLC5GRkYHExETk5uYad1qOcNF4HrtsADkAsk8/PajnCRlW7LQXaFjjOXbaMWPFzuA7M5+wYhc+kRzsDL6d0TXYNTc3Y9CgQSgtLVW9f9asWZg9ezZKS0uxceNGpKenY8SIETh06FCYW2ogGgc7K4AKANZ+/YJ6npAxarBTY8Zg52/3aji7Yt3thD1V7MwSPoxeiY6krlj5e2GwC59IDnYG387oeo7dyJEjMXLkSNX7BEFASUkJHn30UYwePRoAsGjRIqSlpWHp0qW46667wtlU49B4gmKJUY9A1IKdIDifN6MHM1fs5MzYFRuJFbujR8Uv+fLVUyR1xcq3mQx24aO2rI18YK6GFTtt1dTUoK6uDnl5edJtCQkJyMnJwfr163Vsmc5Cda3YUFdfAqV8v8ePG/dDZZZgF+hIV6N0xUbCOXZqO72WlvC3wx0zV0OV5Mua59iFDyt2ujFssKurqwMApKWlOd2elpYm3aemtbUVTU1NTl8RJVTBzqg7RKMe9UVrxU7eFWvEip3BN7gSo+/0IqkrVr6sQ/kZffNNICcHWL06dK9hJkZfx33Bil1oWBRdboIguNwmV1xcjNTUVOmrZ8+eoW5ieIVqVKxZKnaAMYKdmuPHzVERiJSKnX3dCGf3sFbU3puRBlBEUldsuILdmDGoWrsWuXl5qKqqCt3rmEWkBTt57jD4AaRhg116ejoAuFTn6uvrXap4ctOnT0djY6P0tcfo10L1l9qkj1ow6g5R7XqkRgh27kJHqK+fqgV/5rHzVLELZ7BTqxTa75fPAG/wDa7E6Ds9M1dDleSfyRAfeFUDqAR0mb2hqqoKubm5xgmVZjood0f+HhISHD8bdX9pY9hgl5mZifT0dKyWlbXb2tpQWVmJ4cOHu/27hIQEpKSkOH1FFOWHRauNLSt22jBbd6wy2Pkz3Uk4u2LVwpv9fhNtcCVq67WRKnaR2hUbwaqrq1FZWWmcKcGMfvDiC3fbIINvZ3QdFXv48GH8+OOP0u81NTWoqqpCly5d0KtXLxQUFKCoqAhZWVnIyspCUVERkpKSMGbMGB1brTNlRUirYGfUFdWowS5SKnaAf92r4eyKlf/vO3YE7OfLKtdVE21wJUYPdpHaFRsCVVVVUpgqt91WXl4u3Z+dnS1e6SeUjh4FPJx7rotIC3byA0iDV691DXZfffUVLrnkEun3qVOnAgDGjx+PhQsXYtq0aWhpacHkyZPR0NCAoUOHory8HMnJyXo1WX+s2Bkj2Llj1OXojrd57IzSFavcqLo7kjbSyFJPjL7Ti9Su2BAoKChAZWWl021lZWUoKysDAOTk5KCioiIkry2FysJClNuKJGEPle4YfR33hbdeA4PSNdjl5uZC8HDOg8ViQWFhIQoLC8PXKKMLVbAzaqXDqMHOpMPgAWg3eCKcXbFJSY6ff//d+b5OnRw/myXYGX3wRCR3xba3O19mLEglJSWo3rwZmDgR5RCvvT1u3Dhpqq7s7GzNXktJz1DpVaQFO/k2yODbGV2DHQUgVIMnjFppMtvGwWzBziyDJ+Qb1SNHnO9LTHT8bPANroRdseGjts2UV1+CZLVaYe3TB5g4EYAY7PLy8pCfn6/Za7hTUlKC6g0bgHvuCXuo9MqoB+X+YLCjsNDqHDvlDpsVO/8oR2zal5/Zdn7KrlhP65Oe050oq3Ly+2Jjxa7a1lbDb3AlRj9giaSuWLVeDg2DHQDXQVNhGrBhtVphTU8H7rkHQHhDpVdGX8d94W4bZKSDMBWGHRVLbmjVFRuq+fC0ZtRgJ2eicy8AuO60PbXfKJcU81Sxi4lxVO0MvsGVGL1ip1YRNdtBi51yWxeK9da2/cwGkAMg+/TTtX8Nd4w6Ej+Sg53BDyAZ7MwmVMHOxBU73edvMtFoKRcWi+/BzttAC615OlqWrxcWiyP4GXyDKzH6Tk/ePrMHu1CdlyxnC1dWABUArP37a/8aXl5bCpV6dr/KGX0d94X8PbBiRyGjVVes8ijPxBU7XeZvcjdaygw7P2UVTh5Mle33dD5eOCt2nrpiLRZH+DBzsDPSzkJt/TbbQYtdOIKdcvsZzmWlDJUDBoTvtT2JhGBn0nPsGOzMhhU752AnCIAeVxeRt0v+gTfDzk8ZijwFU38GWmjN09Gyu2BnpHDkiZlGxdqXvZHa549QDTiTUx4oh3N7qnztQ4fC99qemOE0Gm9Meo4dB0+YDYMdcPiwY/6m995D+dtvAwDKV62SHhLy+ZtMOr+RC3+DnRFGxapV7ORdsYLgfF1HIzJ6NUPePvu8oc3N5li2SqGa1F1OGa50qNhJmpqAE08M3+u7Y/R13Bcmrdgx2JmNVhsp5fMYdUV1E+xU52964w2UvfEGgDDM3xSpXbHezrEzQlesp4odIP4P5L8bkZm6Yjt3Fr8fPy7+z+Xrixno0RWrZ8WusTF8r+2JWlXazMGOFTsKmVBV7Iy6oqrNGn/4MEpeflms2I0d65i/6cYbkXfDDQDCcAJxtFbsYmPF74IQ3oqdcq46T8GupYXBLljudmjNzQx2avSs2Ck/h7/+Gr7X9oQVO93wHDuz0ep8EZNX7KxWK/Lz85EPIM92c94FF4i35eeH/jI67kKHGYKdp+lOvAU7edUunMHOPlcdoF6xM9FGF4Dxd3r29lksjoodYOhzpNyOjg/HOXYGGDwhMco1Y9XW8aNHjXvajxrlAbz9NAQjHYSpYLAzG3bFut+5hHODESldsRaLf12xgCPYhbMrVu08Ovl98nAdwEY37FPmmGXwhMXiWrEzKLej45XbulB8Ro00eOKXX8L32p64m6TZwOuQC+V8mSaZVonBzmyUH5ZAN1IRFOyk+Zt69gxLkwBEb1cs4BhAEc6KnXISYnf3AQGty2GfMsfoFTv78o2JMUew27oVePZZ9fuUyzoU78FIFTsGO+0o58s0yeh7nmNnNsoPS6BdI8qNgVFXVB92gFaI8zehR4/Qt8fOXVes2Sp2gP/BLlxdse4mIfY0QTEQ2EFKuA9s5O3v2FFc7kb6DMq7Yg0c7KTR8Q8+iHJboClfvly6Pzs7G1attpme2vH99ygAUAJxexTW7YCZumIBw61DHnnqNTAwBjuzUVbaAp2zyMwVuyNHxPbHxjrfHs73YOaKXbBdsfaKncm7YqVQcPAgygsKAADlb70l3R/SKXOU04kYLdjJu2INfI6d6uj45ctRZgt3OTk5qAhDsKveuROVAKphC3bhDC9hrNhVVVWhoKAAJSUl3j8bkRjsWLGjkFB+WKIx2AHiB8s+v5Ydg51vPAU7I1Xs3G1Uf//dfRcJ4PN6oBoK/vMflP3nPwACmzLH5x2fMtj9+quxQpNJumJLSkpcR8dffjnyJkwAYBsdf8stzn8Uigl8lZ/7cP4vlZ/DEAY7+SkLfgU7e1UaMNZ67g0rdhQWWh19GiDY+bQTVO7A7R+0w4eNE+zM3BVrsYg77w4dxAqcP8FOr4od4HzE7Ok+D6RQ8NJLKP/iCzEUdOyIvHnzAAQ2ZY7POz75en3CCeL35max+hIX2GbZr2qKNybpirVareJ7HTsWgBjs8mwj5iUhqthJFd+1a1H+2msAgHL7nZWVyB4yJPSj8wFzdMXaq9KA4dYhj9wdXB49GtRnNdSM2apocPQosHSpeF7YZZf5/neh6orVobTs005QvnFISXFMvqm2cTZIxU7THWwoyJdpjG38VEKCGNT86Yr1oWIX1LLwtbs1wIqdFApWrQJswS4vLs45FPjr/ffF797mEpO/N/lVApqagC5dAnppv6op3phwVKxEuW3TqpdDQbXia/vC4sXI2bMntJOk2ymD3a+/qp+qEiApwP74I8qfeQaAj1f5kS/31FTHZ8Iolzzzhbf5MpXFBYNgsPMgpDvo119H1d13iyfcLl8Oq21iXa+UG6mWlsCOHAxQsfOJu2CntnEwQrD7/XdUL1yo3Q42FNTCWseO4jLVuCs2qLDhqWInDxjBDp6QX01DuZP0gbTj27cP5W++CQAof+opoGtXAG52fPL1Wh7sDh4MONhpyl1XrIG70aTR8d27O98RooqdajcwbPNqXnstsmfM0OR1vFKus8ePA7/9BnTrpsnTB3yVH/lyl6/TBw9q0q6w8NZrwGBnPpoeASvdfTeqAfGE2/nzAw92gLihsnfn+Eqnip20E1y/HuXLlgEAysulDgzXnaBy47Bnj/iz2sbBCF2xK1cC+/aJP3/zDRBM9SdU1Cp29nDqKdjZHyvvivV27dBvvhG/B7IzDbRi5++6HBfnCAVq88t5obrj278fZbbuQdUdn7tg5+floKTPkyCg/IknAADl774r3R/wABB3ExQ3Nxu2Im2FbXS8MtAot3UaBTvVbmAA+QCQmQmEa9moHYzU1WkW7FQD7AUXIO+eewB4OGUh0oKdclolAw+gYLAzAncDBNSoXWLr0KHgg11Tk39/HyDVnWBZGcrKygCo7ATly+akkxw/Hzjg+uTh/KApKnZVEEfEYd8+6Tyb8rlzgYEDAQQ3wrJq1SoUFBaiZM4cbXakalU4+wAKf7piAXE9UlSLpbBx6BDKZ80CALELxxYQfF4Wyo2qvHIkr9haLM5Hzv7uuOPjHaFA7fPlhbTj27QJ5bNnOyo3ixcDcLPj0yjYqX6e/vUvlP3rXwCCuGayu67YxkZUL19u7Iq0srs4RF2xfrUhlNwFO9u2J1iqAbZfP++nLLjbdjc0aNKusFCe4604yDEqBjsFaacEoPzRR8XvtpIzEPwUCE7Pb7utvKYGWLLEt+f350oMnig3BgcPis9tr8qEiOrR37hxyMsTLwzmshOU72hPPtnxc0ODPhtsO3noSE5GAcTqq1xZQ4Pnqo0vvv8e1SNHolIQUP3RR9oHO38qdmrBrrXVJdipho2aGv+XhXKjmpLi+F0egGJinO/z9yAlyK5YacfXtStgC3Z5gOcdn/29xcSI5x/Z+RnspM/T22+j/L33fAuVvnDXFfvSS4E9Xzh5C3Yh6E6WKr7u2hBKaqdEhHoAhS+fk0ir2Cm3QWEqhgSCwU5Bdae0ahXKbCeLBryD9vT833/v+05PLdgFEmiUlQlBEFdUfyt/flI9+svLc78TlL9febA7cMB1g+bnTjEoihODS2Cr2AGOwJqcjLw5cwAEsYOdOtWxcSkrAx56KLDnkVMGJsC/YKc8mV7+O2Rh4/vvUf7kk4GHDeVrywOQfOcQ7AZXHuyC4U+1T97VGUSwkz5PO3YAtmDnNVT6275OnRwVacgOSH05gV4PylAV4IAzf7qcrbBVfN21IZTcVexCQAqw8m2xO5FQsWOwiwzSTqm6GuXFxdodAduff+ZMvDZyJOYcOYILAawDMO6MM5D3+OO+Pb+7rlh/qT3PgQMhD3YA/Lssmrtg19DgOt2GypFgyM4Hkn/gO3SAtWNHWGXvowxAXkJCwDvYqqoqfPDBB3hjzRqcbrutfN8+3yu7vrZdPioWEJe3fDCOWrBTTlirOJdHChuffgrYgl1AYcPTRlXLYBdAlQ5QWbdaW10rN+7IK2LuAqs/5FVULci7Yt1VpH05gV4jfn2ONarYBXWOdTgHmaitv4q57LTaDlrh5jxGNazY6YbBTsFpp2QLdpocAduf//33ccGRI5gD4CyIwS4vI8P35w9VxQ4I35GU7Twun67x6umoz4eKXcgGwKh94JUBNcDAADhXdr+33VZ24EDwXbuA54odIL4Pe3jzJdgpSDuRm28OrH12njaq8v91sBvcACeVdlm3fv/dtXLjjrwrVn4wFWjVOSHBJVQGtTNXdBWXQKUifeONyLMN+gr2gNcbvz7Hyv+/cpvpyzKurgZsU3sERO0c4FDxoWKn+XbQl/k61eZqBMwd7OQHYQx2JqTxrPrSuXWffip1ZWyxff+sthbZVVVOHzi3G2W1YBfIRkTPYGcbvWqFbScYSLDzpSv20CFg7lzx5wBGO3qkFjrq6wHIAqun0aJelJSU4LV77sGcDRscld3UVOT94x/iawSzI/V0jh0QdLCTdiKDB/tewfLWTm/BTj54wt8NbqATLe/f7/y7P5NTa3SOncRicQmVmkw1Y7EAHTuqV6TPPVezA16P9u0D3nnH98crt4fKbebBg27neZO20w89hHJbOHK59qx8Wbrbrnibx1BLYeyKlfhSTJAv9/h4cT1vbIycrthwnvrjJwY7d1pagtspKaidW7fO9n3Otm34rqDAqQLjdqOsFsh++83/BukZ7JQ7QE87Yk9dscpgZzsSlDbOb72FctsyL3/lFSmYaHI+kPIDLwsWVvhYtVEhtb2mBs0bNjjd1wkatV2tYufusmL+Brvt28VzAQGgrS2oZeHXOXYJCeJXa2vwFTsPg4ik/8+GDSgvLQUAlK9cKd756afIhu1aofbnlS9X5WsA2gU75Xv48Ufg738P7LkA52AHiCN3a2udHxPiUejSsrZdGQQAyj/8ULrf7WdBGezUzic+eND5QNHGp2vPyivl7goA+/d7nwpIK26CndPUUrarqXicWsof/gY7e2W6sdHcFTt2xZpcS0twOyUF6dy9+fNRvmaN82SWAwYgu6TE+Q/cHQmqVeyUlQNfGCnYedqZKScojo0V237ggGulxfY8qhvnzz5D2WefAdDofCBPH3i7AKq+Hg8AGhtdDgAColaxczc/k4/BTtqJTJ+Octtcg+WbNkkPywZg9RR01Ch3DPJlLF9X7ZWXlBSxUhJssPMwSbDqurV4Mcps5+DmQLbNaGgA0tPVX1M+OEGLbirbZ6oKti7TkSNR/uOPAIDyf/5TepjPO3N58AScgp10wGubgBkIzbmsqst6yRKU2c4zlT7Hym2i8kBXbZv522+qwU511P6llyJv4kQAKpVyd1Xa1lbxXD/5ZyVU3IyK9XtqKX/40kuk/PyeeCKwa5f4ubCF3nDMiajJKQlA8KPvw4jBzh3lZLdtbUGdoGxfoQq+/14Mc5BNZhkXB1itjp3jsWMoLygAAGkme8C2UXa3kfKXu8ET4aDcGHoKlPL3Gxsr7nB//VW9YnfgACAIjo1zQQHK9+8XN86nnII82zkzmpwP5EuwC2BknNT2zz5D+Zw5KANwNYAdAMampuJq5QFAINQqdu7ChY/BTnUnsn49bLU7MfAcPAikpfneTvk6GhvrfvBEsMFOeYCwf7/bYKe64x8+HHmTJwP/+Q+y33rL8eD6evfBTj54Qv5agXyWAekzVQDbIAdbqAOAsn//G2X//jcAP3bmahU7Gyts4VXW7lCcy6q6rC+5BHm33w5A9jlW2w7I+bHNVB21f+aZ7rucPZ2f+euv4Ql2ahW7AwdQMmsWqrdvd15+f/gD8q6/HkAA28HOnR0Vel+6mtUqdoD4/zpyBOjUKbQXAbAJ6jU8bYMY7ExI2c3Q0ODfTklFdXU1KmtrcabyDtsIJtWd4/vvo8x2/cmcnBxUDBni+sSB7AzUNga2c8RCThmaFSO4nKgd9f36q/ielRv0tjbg118dG+enngJswS4vNVXb84E8dMU6PcbPuQGltsfEALZgdwtsBwCJidrMZq8W1uST5Mp3jD4GO9WdcFYW8rZvB2A7naGuzr/PkDLUu6vY2Zev/f6mJv+6wZQ75/37gTPOUH2o6o7/jDPEdWvHDucHe/o8yStinTqJ/9uWFs+fBU9sn6kSqAxyuPBC5N19NwA/duZqFTulQHoK/KC6rLt2df0cK4O5rxU7X3ka4eop2O3fL16BItTk2/ITT5Q+G9b0dFjPO895+Z1zTuDbQfk5ib6sp8oDSHmF1L7urF4dWFvc0LwCqAx2HDxhcsrw8fPPQQc7rBM71U6B4ty9+nqgvd2xc/zgA5QvXao+1Yr9/CU5rbpif/7Z/+cJhLJipzx3R04Z7NLSgB9+EM/xUPtg7dnjGIov714M9AR5d3yp2AFiSPJlzicltWqfVifrql1STF41kocmH4Od6k64Y0c47ULq6oBBg3xvp3wdVZ6L5q4rFhADf2ur84AQT9Qqdv6w79yV67Wvwc5iEdfrn34K/ODK9tpWyM7xgyJ4+sO+7O3LVmUapKrNm1G9ZAnQ0oLy6dMBAOX/+Y90f0jmttu61fU25f/vyBHn+RUD7J2Qupzl1wdV8hbswkEe7DIzHZ+N3buBXr2cHxvMNkT+Xn1ZT5XbGds1fKsAVC9cCHz7LcptV0jR5DJ4cFTn3n//fceFAOyX2rNdwtKv1/B0cMlgZ0LKYLdvH3D22X4/jdS9euCAeAI/gB8ATILjyNp6/DjQ0ODYOTY1AbZg5zLVysKFri+iVVesXsHO0wgu5cYhI8Px+08/uT5+715g8GDx58REx8ZZZQRcUHwNdnV1gQW7w4ddB++0tIiBNtgLT3ur2AUQ7FQpd2yeArwa5dFyQoI4uu7oUecdszLYAeJnyNdgp9w526/164X0/3F3OTZfKtH2YN2tm7g+//ab8zyCvvI0IjeQHZAy2KlU7Ao++ACViu1R2VtvoczWHa3Z3HYdOiC7rU1c1tu3u54Wo3bQ9tNPgL06GWDFzgpbl7OnbYenYOdme6p5VUke7M44A7Cf27pjB3DBBQBk62qg2w5BcF7HfFlPlefy2k5LKABQWVjo9FBNLoPX3Ax8/DEAYPHixfj++++d7i77739R9t//+vcayoNL+baPo2JNSC3YBUC1e9X2BchOtq6tdZSqPZ10rxbIAhlar2ewUy5bfyp2p5zi+F0t2NlO3AcgThwM2/LVePoal8DjLrzV1QEDBvj//IcPqw/e+flnoG9f/59PTq1ip0Wwsy1jaSei3HkGG+zsAfq331zvA1yDna8XQVcbUeoDK2z/H/syUoarvXvd/7F88ATg6A0QBPHzbKtu+Ez5mYLs/xDIZQLty9e+41YJdiW9e6N69mzg2WdR/r//iT0MSUnIe+018fW1OJf12DHn0dW//w588QVw0UWOx6iFq507PQc7fyqju3a5v89ToHbzd5qfVyYPdvJtw44dTpXcCsD3gx0lte3nvn1A797u/8ZNxa4EQPWECcDChc6nbdgqa/6uN1LxZMUKlNumxTmtZ09ccsklwN692POf/+C/CPBCA8rtTGys41xDVuxMSLmh9LSR9kDqXv3mG5TPmuU8Ghayasy2bY4A0NTkXK2Rny8k/7DYz8vZv19c0fw5UdddsAvHEH3lxlB5bpKc8ohJXrFT23Du3u34WX6eZF2dtu9NGXhOPVX9cdu3A5df7v/zuxt4sXdv8MEuVBU7e7csbDsRZSWlpsa/diq7QQCxO1YZGJXn2AH+jfBWttN2XqBHHTo4/m7nTvG7cr1WVAyc2HfG9suZyU/z+OUX/4OdSsCwwvZ/COZa0h4qdtaDB2HNzwf+8x/AFuzyLBZtz2VV+xx8/LFzsFOr2Nm3KfbzXAGn88+k/5ka5bZR7QDSzlPFzt3faT37gKdgp5yWxNP79kQtwH7/vedgJ29XbKxUsbMCsPbqJXYT794trjcNDcgfMyag7bNa8eT91aul8/f62W7LA5BvtToCvy/cHUAePmzoaVtCe8V3M1MGu+pq9cd5YbVakZ+fj/xBg1xGw+ZDdj6M/NyRxkZpo2wFnCtp8g9Lv36Onz2FIzVqwa6lJTwrq1qwcxdk5EeK8fHeg922bY6f5Tu0I0e0HRzia7D79tvAnt/dzvi77wJ7PrlgK3byc93kVThvc1tt2eL5fiVlqAfUR5naN7jydcOfA7FAKnbyHZA9CCqfR+18MDv7em2viCmDnb88VY5++MH/5/OhKxY7dkijGyUqlcOgqH0ObN1tEk/BTr4OnX66Yz3yFHCU26KaGvfTT8n/5ykpwPr1jt9lwa6qqgpLlizBkquuQvl99wEQ55RbsmQJlixZgqqqKvft8Ua+jczKcvy8Y4f6PJOBUFu/PK3fynbFxzsfrNTWOl+jef9+52Xnh5KSEixevBiLIRZNAGDcueeKt91xB26VP3jBAv+eXG1baX8fv/yifU+QRhjs3FFuoGRzcgXE28m68g+JssQrD5XyjVj//o6ffew+ksg3ePLzVbQIDt4ol60guA/O/gY7+fMoN2pavje9gl0wOwC7YCt2PXs6TqbfuNHxGG+VoS1b1LvF3FE7WlZbzvb75NUDT91nSspg8OOPnisx9sEZdnV14gGRcudXU+P+gEVZsZN3GwdyAOIpUNXU+H/ZNGWwU6sgtrcD33wDNDc7ehja2/37H3ujtk5t2OB8EKEW7OwHEfLtR6dOjqvceDoQVv7PDh1yX22WL9eCAmDYMMdpGbK/KSgowNixYzH2ww+l03DKysrE28aORYFtequAyA/2k5OBHj3En7dtc92XeKoie6K2/njbnnoLdg0Nzj1Ttm5Sf0nFEzh6wvK6dRNv69sXV8lf4403/BtI52kb1N4ecE9eqDHYuaPcUO7aFdwopwMHXE+Gl1fcvv7a8bPywyi/T/5hkQc7f4/E5Cv3eec5ftYiOHijdvQnf49ynoKdWlfHzp2O/51yA712rV/N9EjtHDu1c7qqqgK7Zqy7kBTsAQagfhQqn4pAfpCgFuxiYhzrzC+/OEKUt4pdc7NzRdUbX4Od/T3Ig52n7jMl5U7r6FHgf/9z/3i1/82nn/q3XtvXCXvFTl6JlJ8n6itPFbvjx8Uw5A/lOXZupn/B118DBw449zD4e5DpidqyPnYMkHe9qYWOTZvE9Vy5/ejTR/y5ocH9AAq1ML5unettyte2n79mr5rt3SvtM0pKSrB44ULnqtLVV4tVpcWLURLM/JTy7UtcnGPk+cGDrlXyPXv8+2zYqa1fX37p+W/kyz4uDuja1fFZ3bcPOHjQeb1Ztizg6za7sFdkletmfT1gG0DhE7VtUKAHkGFkimD38ssvIzMzEx07dsTgwYPx6aefhv5F1S6X89FHgT+fcgUDxCrdsGHiz9u2ObpclaNtZFMIOH1Y5KN0P//cv/bIP0BDhzp+1iI4eKO2kVB2r9gpN8w9ejg2DmrPIwiODY5ypxDM/0/tdewsFvHrnHNcH9fU5H7n7om87WPGOELj//4X/HUo1cJaYiJw7rniz1u2ODZYao8FgOHDHT+/95743ZfLDLn7P6tRO8dO7ZwetdDnzwZXbWfiKQiphY3Vq9Wfx93IO2XFTn6QFsjBlbfr1MouJeUTZcXO3fWcV6927Y2wX2JNC/KQJe9mfPttx89qFZhDh8RtqnL7cdZZjt/ddf2pBbtPPlF/rPx/bh8dLf9s2F7DarUiPyfHuap0wgliVSk/P7hBFMr3KJ/rdM0a18cHMnec2rr9zTeezxeUB874eHFdsgfrzZvVr+Er39f54/hxwGJxFE927hT/j2rts11ezSfetkGBhOQwMHywe+utt1BQUIBHH30UmzdvxkUXXYSRI0dit/wk+VBQ69rw50LUSsqN3yOPiN8vvdRxm212eJeK3eefO0ro8g+x1eoo+1dU+Fdiln9QzzvP0R37wQfq599pSW0n9NFH6jtM5UYrMVF9LrQLL3T8vGKF+B6Ur7Nunf8n8LujFnjcjX6VzdHkM/nOZcEC4LbbHK8byPPJqVXsAOC66xw/v/ii4/Xs5MHuD39w/Dxnjvh/8uUkfX/arnaOnaeu2IwMR4Vp3Trfr/yhttPy9FlXe59vvaUebN99V/38LOU5dv37O+Zd9FYJUeMt2C1e7PvnWj7gwL5s3U358dFHrju3efO0646VL+s//MFxCsC//uWouKmFK/tjlNuP3FzH7+4Cjtp6s2KF+sG+fLnbX9s2xQgAQHZtW5duu8pKbba18nYlJTkO0AD1A6k33vD/NdwdRNsP6tQolz3g3DtkJz8NpLg4sHXHdtUhK2zFk9ZWYNEi9VOgPvwQ+Oor357X2zZI4wmWtWL4UbGzZ8/G7bffjjvuuAOAWNJetWoV5syZg+LiYt+f6IcfxBOUfR11I0/6ycniRvudd8Qd3oUXiudrxMeLG2ZfnlPetTNvHnCr7ZTOG24Qr5AAAIWFYpeMWvXt9tvFx8mfp0MHccTlsmViGLzzTuCOO8TJZjt29DzHkLzqc+KJwFVXiTuhujrgvvvEiWZPOkl8DvskqvbvwY4slZ9wf+qp4o6hsRGYNEl87RNPFN9bfLzzztK+cRg+XDzik7vpJnEahKNHgVdfFU+UVnPbbcCMGeK6EB8vvif5+/L1vcnbZf+bO+8EnntO/PmWWxw7lpdeEkPfWWeJ/xf70aun17KfvxIXJy6LG24AbJdEw+OPi90ap5/uWE7enk9OfkqB/G9uuw148klxI/788+J98hPv5Y/Nzhb/D+vXiwcd117rPCG0UmamGKo//ljceOfliaNrO3Rw33a1y4YNHy6ul/JuNPt9cXHAqFHiTnj/fuCKK4A//1ms9Jx8svprCIIjPJx3nvh/3bpV7LZ/+mngssvEwSL2z7nFon4+aH2949y4Dh3E//emTeJ6+vjjYmju3Fnc+cfGulbs7F1oGzaI1ca77gKuv14Mq/bPclyc+4Blr/LLR37K7d4tfrYmTBAfY1/mal9q3U+A+D+zV/4uuUSsBjU3uwahLVuA++8Xt3EpKeJ7jI8PbLshP/f45JPF53zpJfHA+/bbxQNk+WPuvBMoLRX/r88/7zyCOz4eyMlxXG/69dfF/9OQIeL2vGNHsY1qgwKamsTXu+8+5+2ifF2wv9all4rP19wMzJ8vBq2zz5Z6DKSq0t694mfunnvE/0lCgmMZefsuJ/9/JyWJ67HFIi4DtQLIp58CM2eK2/zUVNd1wb49lJMfEF9wAWC77jb++ldxf3PqqeLy69DB8bfyz6882C1d6vzcN94oBq2qKvEzM3Wq+H/u3FlsW0yM47u7dUjtf/bII869XzNmALYpVTBunLiedO/uvB+wf9nJ/97+Wbj4YnG5NTaK7yU1Vfx89+qlPseoLz0ZGrMIgrvhPvpra2tDUlIS/vnPf+KGG26Qbr///vtRVVXlMsQZAFpbW9EqO4JrbGxEr169sAeAmylkPevcGZg+HXj00UD+2pXFIu505KHrxhvVuwlTU8UPiruut127xPNZLrssuDZ98IG4Mx89OrjnCcSiRcDEib4dudbVieHh88+BK690vm/xYvEIeO5c17+75BIxEIfqWrhr1zqqiB9+KHa9Tp4sbjz9HYWlZL/SBuB+PQnGU08BU6Y4fn/2WTHcqZk3D/jjHx2/V1WJ657aOYSXXeaoFmRkiJ+fe+8NvJ3vv++oyr7xhnObp01zfD5//FHc8AZwnV5ccon4GbCNWvTJFVeIR+3yKkPXrmLAHz/e+9+fe67jf7p0qbiTD8aAAWKb/v538fd77wX+8Y/An2/4cEfVqb5eDKlnnikG+SFDwjsqcM4cMZgNHux+sMjs2WJwWbHC9b4//EEMWlOnit99MWECsGSJb+9z6VLg6qvFnx9+GLBNSB9WDQ1iMMnPdz2X7Oqrxc9RsB5/XAz4/py3+fXX4oHo3r3iOi+vMj78sPjZHjUq+La5Excnnup00UX+necr9/HHjm7uefOABx/06c+aAPQEcPDgQaTKZxQIJcHA9u3bJwAQPvvsM6fbn3rqKeGMM85Q/ZsnnnhCAMAvfvGLX/ziF7/4ZYivHTt2hCM2CYIgCIbvigUAi6L8KgiCy21206dPx9SpU6XfDx48iN69e2P37t3hS8tRqqmpCT179sSePXuQ4u4SW6QJLuvw4bIOHy7r8OByDh97r2EX+fW4Q8zQwe7kk09GbGws6hTXEq2vr0eafEJPmYSEBCTIT6C1SU1N5QocJikpKVzWYcJlHT5c1uHDZR0eXM7hExPIpf0Cfa2wvVIAOnTogMGDB2O1YuTJ6tWrMVw+pJyIiIiIjF2xA4CpU6fi1ltvxZAhQzBs2DC89tpr2L17N+6++269m0ZERERkKIYPdjfddBN+++03zJw5E7W1tRgwYAA++OAD9PZ08WGZhIQEPPHEE6rds6QtLuvw4bIOHy7r8OGyDg8u5/DRY1kberoTIiIiIvKdoc+xIyIiIiLfMdgRERERRQgGOyIiIqIIwWBHREREFCFMEezWrl2LUaNGISMjAxaLBe+++67T/RaLRfXr2WeflR7T2tqK++67DyeffDI6deqEa6+9Fnv37nV6noaGBtx6661ITU1Famoqbr31VhyUX8g4CmixrHNzc13uv/nmm52eh8va+7I+fPgwpkyZgh49eiAxMRH9+/fHnDlznB7D9do7LZYz12nfeFvWv/zyCyZMmICMjAwkJSXhyiuvxPbt250ew3XaN1osa67XvikuLsa5556L5ORkdOvWDddffz22Ka45KwgCCgsLkZGRgcTEROTm5qK6utrpMeFat00R7JqbmzFo0CCUlpaq3l9bW+v09frrr8NiseDGG2+UHlNQUIAVK1Zg2bJlWLduHQ4fPoxrrrkGx2UXnx8zZgyqqqqwcuVKrFy5ElVVVbj11ltD/v6MRItlDQCTJk1yetyrr77qdD+Xtfdl/cADD2DlypVYvHgxtm7digceeAD33Xcf/v3vf0uP4XrtnRbLGeA67QtPy1oQBFx//fXYuXMn/v3vf2Pz5s3o3bs3Lr/8cjQ3N0uP4zrtGy2WNcD12heVlZW49957sWHDBqxevRrHjh1DXl6e07KcNWsWZs+ejdLSUmzcuBHp6ekYMWIEDh06JD0mbOt22K5KqxEAwooVKzw+5rrrrhMuvfRS6feDBw8K8fHxwrJly6Tb9u3bJ8TExAgrV64UBEEQvvvuOwGAsGHDBukxn3/+uQBA+P7777V9EyYRyLIWBEHIyckR7r//frd/w2XtSm1ZZ2dnCzNnznS67ZxzzhEee+wxQRC4XgcikOUsCFynA6Fc1tu2bRMACN9++61027Fjx4QuXboIc+fOFQSB63SgAlnWgsD1OlD19fUCAKGyslIQBEFob28X0tPThaefflp6zO+//y6kpqYKr7zyiiAI4V23TVGx88cvv/yC999/H7fffrt029dff42jR48iLy9Pui0jIwMDBgzA+vXrAQCff/45UlNTMXToUOkx559/PlJTU6XHkDO1ZW23ZMkSnHzyycjOzsZDDz3kdNTCZe2bCy+8EO+99x727dsHQRCwZs0a/PDDD7jiiisAcL3WirflbMd1Ojitra0AgI4dO0q3xcbGokOHDli3bh0ArtNa8WVZ23G99l9jYyMAoEuXLgCAmpoa1NXVOa23CQkJyMnJkZZTONdtw195wl+LFi1CcnIyRo8eLd1WV1eHDh064MQTT3R6bFpaGurq6qTHdOvWzeX5unXrJj2GnKktawDIz89HZmYm0tPT8e2332L69On43//+J13zl8vaNy+++CImTZqEHj16IC4uDjExMZg3bx4uvPBCAFyvteJtOQNcp7XQr18/9O7dG9OnT8err76KTp06Yfbs2airq0NtbS0ArtNa8WVZA1yvAyEIAqZOnYoLL7wQAwYMAABpWaSlpTk9Ni0tDbt27ZIeE651O+KC3euvv478/HynIxV3BEGAxWKRfpf/7O4x5OBuWU+aNEn6ecCAAcjKysKQIUOwadMmnHPOOQC4rH3x4osvYsOGDXjvvffQu3dvrF27FpMnT0b37t1x+eWXu/07rtf+8WU5c50OXnx8PN555x3cfvvt6NKlC2JjY3H55Zdj5MiRXv+W67R/fF3WXK/9N2XKFHzzzTculU/AdVn5spxCsW5HVFfsp59+im3btuGOO+5wuj09PR1tbW1oaGhwur2+vl5K2Onp6fjll19cnvPXX391SeHkflmrOeeccxAfHy+NyOKy9q6lpQWPPPIIZs+ejVGjRmHgwIGYMmUKbrrpJjz33HMAuF5rwZflrIbrdGAGDx6MqqoqHDx4ELW1tVi5ciV+++03ZGZmAuA6rSVvy1oN12vP7rvvPrz33ntYs2YNevToId2enp4OAC5VNeV6G651O6KC3fz58zF48GAMGjTI6fbBgwcjPj5eKi8D4ujOb7/9FsOHDwcADBs2DI2Njfjyyy+lx3zxxRdobGyUHkMO7pa1murqahw9ehTdu3cHwGXti6NHj+Lo0aOIiXH+iMbGxqK9vR0A12st+LKc1XCdDk5qaiq6du2K7du346uvvsJ1110HgOt0KLhb1mq4XqsTBAFTpkzB8uXL8cknn7iEY3t3tny9bWtrQ2VlpbScwrpu+zUURCeHDh0SNm/eLGzevFkAIMyePVvYvHmzsGvXLukxjY2NQlJSkjBnzhzV57j77ruFHj16CB999JGwadMm4dJLLxUGDRokHDt2THrMlVdeKQwcOFD4/PPPhc8//1w466yzhGuuuSbk789Igl3WP/74ozBjxgxh48aNQk1NjfD+++8L/fr1E84++2wuawVvyzonJ0fIzs4W1qxZI+zcuVNYsGCB0LFjR+Hll1+WnoPrtXfBLmeu077ztqzffvttYc2aNcKOHTuEd999V+jdu7cwevRop+fgOu2bYJc112vf3XPPPUJqaqpQUVEh1NbWSl9HjhyRHvP0008LqampwvLly4UtW7YIt9xyi9C9e3ehqalJeky41m1TBLs1a9YIAFy+xo8fLz3m1VdfFRITE4WDBw+qPkdLS4swZcoUoUuXLkJiYqJwzTXXCLt373Z6zG+//Sbk5+cLycnJQnJyspCfny80NDSE8J0ZT7DLevfu3cLFF18sdOnSRejQoYNw2mmnCX/+85+F3377zelxXNbel3Vtba0wYcIEISMjQ+jYsaPQt29f4e9//7vQ3t4uPQfXa++CXc5cp33nbVm/8MILQo8ePYT4+HihV69ewmOPPSa0trY6PQfXad8Eu6y5XvtObTkDEBYsWCA9pr29XXjiiSeE9PR0ISEhQbj44ouFLVu2OD1PuNZti63RRERERGRyEXWOHREREVE0Y7AjIiIiihAMdkREREQRgsGOiIiIKEIw2BERERFFCAY7IiIiogjBYEdEREQUIRjsiIiIiCIEgx0RERFRhGCwIyIiIooQDHZEREREEYLBjoiIiChC/D9ayah2qHUdVAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "from netCDF4 import Dataset\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "year1  = 1700 \n",
    "yearn  = 2000 \n",
    "\n",
    "# read CAM forcing\n",
    "offset = 500 # forcing data shifted by 500 years from -499-2000 to 1-2500 \n",
    "nc = Dataset(f'CCSM4_volcanic_forcing_0001-2500.nc','r')\n",
    "colmass = nc['colmass'][:]\n",
    "MMRVOLC = nc['MMRVOLC'][:]\n",
    "lev = nc['lev'][:]\n",
    "lat = nc['lat'][:]\n",
    "date = nc['date'][:]\n",
    "nc.close()\n",
    "year_ccsm = np.int32(date/10000) + np.int32((date - np.int32(date/10000)*10000)/100)/12 - 1/24 - offset\n",
    "\n",
    "# read eruption data\n",
    "file = open(f'CCSM4_volcanic_forcing_0001-2500.txt','r')\n",
    "year = []\n",
    "hemi = [] \n",
    "load = []\n",
    "for line in file:\n",
    "    year.append(int(line[0:6])+0.5-offset) \n",
    "    hemi.append(int(line[7:8]))\n",
    "    load.append(float(line[8:15]))\n",
    "file.close()\n",
    "\n",
    "# compute global load from colmass\n",
    "area_earth = 4*np.pi*6371000**2\n",
    "weights = np.cos(lat/180.*np.pi)\n",
    "weights = weights/sum(weights)*area_earth\n",
    "colmass_global = (colmass*weights).sum(axis=1)*1e-9\n",
    "\n",
    "# plot\n",
    "fig, ax = plt.subplots(1, 1)\n",
    "h=[]\n",
    "h1, = ax.plot(year_ccsm, colmass_global,'r',linewidth=2,label='CAM forcing file')\n",
    "h.append(h1)\n",
    "for ind, y in enumerate(year):\n",
    "    h2, = ax.plot([y,y], [load[ind],load[ind]],'k+',label='eruption file')\n",
    "ax.legend(handles=(h1,h2))\n",
    "ax.set_xlim((year1,yearn))\n",
    "ax.set_ylim((0,90))\n",
    "ax.set_ylabel('Load (Tg)')\n",
    "plt.tight_layout()\n",
    "plt.savefig('verify_colmass.png')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "2a860dc6-c1c8-4296-9c04-af4db3ffff17",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0.0 0.0 0.0 ... 0.0 0.0 0.0]\n",
      "[1.2150000e+03 1.0115000e+04 1.0215000e+04 ... 2.5001115e+07 2.5001215e+07\n",
      " 2.5010115e+07]\n"
     ]
    }
   ],
   "source": [
    "print((colmass*weights).sum(axis=1))\n",
    "print(date)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "id": "40fa4a79-a56f-435f-bd7d-3be69240523a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1.0986122886681098\n"
     ]
    }
   ],
   "source": [
    "print(np.log(3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e6e177dd-7178-44ed-b0a3-fde3229f9462",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.15"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
