IOSACal mixed calibration curves

IOSACal mixed calibration curves.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from iosacal.core import CalibrationCurve, RadiocarbonDetermination\n",
    "\n",
    "from iosacal.plot import single_plot, stacked_plot\n",
    "\n",
    "import pkg_resources"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# needed to plot results in the notebook\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# this pkg_resources call is going to be unnecessary with the next version\n",
    "curve_filename = pkg_resources.resource_filename(\"iosacal\", \"data/shcal13.14c\")\n",
    "curve = CalibrationCurve(curve_filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve.mixing('marine13', 0.5, 0.05) # these are simple numbers, real values will be different"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "r1 = RadiocarbonDetermination(2233, 25, 'test')\n",
    "cal1 = r1.calibrate('shcal13')  # calibrate using only SHCal13\n",
    "cal2 = r1.calibrate(curve)      # calibrate using the mixed curve\n",
    "cal3 = r1.calibrate('marine13') # calibrate using only Marine13"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# reservoir effect with same proportion values as above\n",
    "curve_filename = pkg_resources.resource_filename(\"iosacal\", \"data/shcal13.14c\")\n",
    "curve = CalibrationCurve(curve_filename)\n",
    "curve.mixing('marine13', 0.5, 0.05, 50, 1)\n",
    "cal4 = r1.calibrate(curve)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAIZCAYAAACI8tb1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeVQUV94//velafZN2REQwQ1kBzEYxR0X1Gg2Jyb5JhMTM/lNJjPPM/HEzGKczJJnYmY0mcRMzEQTl4lEEyWJRnEjRqMiuKC4AoIiKouy093QfX9/AH00AgJCNwXv1zl9pLtuVX3qNuCb27eqhJQSRERERERKYGHuAoiIiIiI2ovhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIqIkQ4lMhxF+6aFtSCDG4E+sFNK1r2RV1EBH1NgyvRNTjCSHGCCF+FEJUCCFuCiEOCiFGNi17VghxwNw1moMQYqkQYr256yAiMiX+ZU9EPZoQwgnAtwBeAvAFACsAYwFozVkXERGZB0deiainGwoAUsrPpZR6KWWdlDJVSpklhAgG8G8A8UKIaiFEOQAIIZKEEMeFEJVCiCtCiKW3b/C2kdzypuXP/nSnQghHIcQ+IcR7opG1EOIdIcRlIcQNIcS/hRC2t7VfJIS4JoQoEkI819YBCSHShBBvCSHSm0aTU4QQ/Vtp6yOE+LppxDlHCPFC0+vTAPwOwLymYz/ZkU4lIlIqhlci6ukuANALIT4TQkwXQvRrXiClPAvgFwAOSSkdpJQuTYtqAPw/AC4AkgC8JISYAwBCCH8A3wH4FwB3AJEATty+QyGEK4A9AA5KKV+RjffR/jsag3QkgMEABgBY0tR+GoBXAUwBMATA5HYc1/8D8BwAHwANAN5rpd3nAAqb2j0K4G9CiElSyh0A/gYguenYI9qxTyIixWN4JaIeTUpZCWAMAAngYwAlTSORnm2skyalPCWlNEgps9AYAMc1LX4SwO6mkdx6KWWZlPL28OoD4HsAm6SUfwAAIYQA8AKA/5FS3pRSVqExOP6saZ3HAayRUp6WUtYAWNqOQ1t3W/s/AnhcCKG6vYEQwq/p2F+TUmqa6vwPgKfbsX0iol6J4ZWIejwp5Vkp5bNSSl8AoWgMmCtaay+EGNX0kX+JEKICjaOzbk2L/QDktrG7JAC2aJyO0MwdgB2AzKapBuUAdjS9jqZ6rtzWvqAdh/XT9urbamzmA6A5LN/edkA7tk9E1CsxvBKRokgpzwH4FI0hFmgckf2p/wL4GoCflNIZjUFUNC27AiCojV18jMZgul0IYd/0WimAOgAjpJQuTQ9nKaVD0/JraAzFzfzbcSg/bV/ftJ/bFQHoL4Rw/Enbq01ft3TsRES9GsMrEfVoQojhQojfCiF8m577AXgCwOGmJjcA+AohrG5bzRGNI5YaIUQcgPm3LdsAYLIQ4nEhhKUQwlUIEfmT3b4M4DyAb4UQtlJKAxpD7XIhhEdTHQOEEFOb2n8B4FkhRIgQwg7AG+04tKdua/8mgM1SSv3tDaSUVwD8COAtIYSNECIcwIKmY2g+9gAhBH+XE1GfwV94RNTTVQEYBeCIEKIGjaH1NIDfNi3fCyAbwHUhRPPI5f8H4E0hRBUaT6r6onljUsrLAGY0rX8TjSdr3XGyU9MJWgvROEqbIoSwAfAagBwAh4UQlQB2AxjW1P47NE5j2NvUZm87jmsdGkeQrwOwAfBKK+2eABCAxlHYLQDekFLualq2qenfMiHEsXbsk4hI8UTj72giIjIVIUQagPVSyv+YuxYiIqXhyCsRERERKQbDKxEREREpBqcNEBEREZFicOSViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBTDsiON3dzcZEBAQDeVQkREREQEZGZmlkop3Vta1qHwGhAQgIyMjK6pioiIiIioBUKIgtaWcdoAERERESkGwysRERERKQbDKxEREREpxj3DqxBioRAiQwiRUVJSYoqaiIiIiIhadM/wKqVcJaWMlVLGuru3eNIXEREREZFJcNoAERERESkGwysRERERKQbDKxEREREpBsMrERGZXHl5OVauXNnp9VesWIHa2to7XnvrrbewYcMG4/OIiAg88cQTd7T59NNPUVRU1On93o/x48ff941+Ro8e3el1HRwc2lze2fdk06ZNCA4OxoQJEzpbWrt8+umnePnll+94rSv6tLP+9re/dXidtup99NFHkZeXh9raWiQlJWH48OEYMWIEFi9ebGyj1Woxb948DB48GKNGjUJ+fr5x2bRp0+Di4oKZM2fesd0FCxYgIiIC4eHhePTRR1FdXQ0A+Pbbb/HGG290+Bh6AoZXIiIyue4Ir6mpqUhMTAQAnD17FgaDAfv370dNTY2xjTnD6/3Q6/UAgB9//LHb9tHZ9+STTz7BypUrsW/fvjteb2ho6KrSeqTOhNfWZGdnQ6/XIzAwEADw6quv4ty5czh+/DgOHjyI7777DkBjX/fr1w85OTn4n//5H7z22mvGbSxatAjr1q27a9vLly/HyZMnkZWVBX9/f7z//vsAgKSkJHz99dd3/RwpAcMrERGZ3OLFi5Gbm4vIyEgsWrQIALBs2TKMHDkS4eHhxhGhmpoaJCUlISIiAqGhoUhOTsZ7772HoqIiTJgwwTjaV1lZCZ1Oh+ar4vz3v//F008/jcTERHz99dcAgM2bNyMjIwNPPvkkIiMjUVdXh4CAAPzud79DfHw8YmNjcezYMUydOhVBQUH497//DQCQUmLRokUIDQ1FWFgYkpOTAQDXrl1DQkICIiMjERoaih9++AFA4wjnb3/7W0RHR2PSpEm4/TKTmzZtQlxcHIYOHWpsr9frsWjRIuOxf/TRRwCAtLQ0TJgwAfPnz0dYWJhx283efvtthIWFISIi4o7RuWaXLl1CfHw8Ro4ciT/+8Y/G16urqzFp0iRER0cjLCwMKSkpLb4nrbW73ZtvvokDBw7gF7/4BRYtWoRPP/0Ujz32GGbNmoXExMRW+y4tLQ3jxo3D448/jqFDh2Lx4sXYsGED4uLiEBYWhtzc3PZ+Kxl9/vnnCAsLQ2ho6B2h7pNPPsHQoUMxfvx4vPDCC8bR29zcXDzwwAMYOXIklixZ0urI9Pr16xEXF4fIyEi8+OKL0Ov1WLx4Merq6hAZGYknn3zyrnVeeuklxMbGYsSIEe0a3dywYQMeeughAICdnZ3x+9rKygrR0dEoLCwEAKSkpOCZZ54B0DhSu2fPHkgpAQCTJk2Co6PjXdt2cnIC0Ph9XFdXByEEAEAIgfHjx+Pbb78FAHz99ddYsmTJPWvtEaSU7X7ExMRIIiKi+3Xp0iU5YsQI4/OdO3fKF154QRoMBqnX62VSUpL8/vvv5ebNm+Xzzz9vbFdeXi6llHLgwIGypKTE+PqXX34p//jHPxqfDxkyRObn58udO3fKWbNmGV8fN26cPHr0qPH5wIED5cqVK6WUUv7mN7+RYWFhsrKyUhYXF0t3d3cppZSbN2+WkydPlg0NDfL69evSz89PFhUVyXfeeUf+5S9/kVJK2dDQICsrK6WUUgKQ69evl1JK+ac//Un+8pe/NO77f//3f6WUUm7btk1OmjRJSinlRx99JP/85z9LKaXUaDQyJiZG5uXlyX379kk7OzuZl5dnrNfe3l5KKeX27dtlfHy8rKmpkVJKWVZWdlcfz5o1S3722WdSSinff/9947r19fWyoqJCSillSUmJDAoKkgaD4a73pLV2P3V7n65Zs0YOGDDAWE9rfbdv3z7p7Owsi4qKpEajkT4+PnLJkiVSSilXrFghf/3rX9+1nzVr1kg3NzcZERFhfNjb28ujR4/Kq1evSj8/P1lcXCzr6+vlhAkT5JYtW+TVq1flwIEDZVlZmdTpdHLMmDHG9yMpKUn+97//lVJK+eGHHxr753ZnzpyRM2fOlDqdTkop5UsvvWTs05baN2s+/oaGBjlu3Dh58uTJu/rqdgkJCTIrK+uu12/duiUHDRokc3NzpZRSjhgxQl65csW4PDAw8I6fg3379smkpKS7tvPss89KDw8POX78eOP3jJRSrl+/Xr788sutHoc5AciQreRRjrwSEZHZpaamIjU1FVFRUYiOjsa5c+dw8eJFhIWFYffu3Xjttdfwww8/wNnZucX1d+zYgenTpwMAjh49Cnd3dwwcOBCTJk3CsWPHcOvWrVb3PXv2bABAWFgYRo0aBUdHR7i7u8PGxgbl5eU4cOAAnnjiCahUKnh6emLcuHE4evQoRo4ciTVr1mDp0qU4deqUcdTLwsIC8+bNAwA89dRTOHDggHFfDz/8MAAgJibGOF8xNTUVa9euRWRkJEaNGoWysjJcvHgRABAXF4dBgwbdVfPu3bvx85//HHZ2dgCA/v3739Xm4MGDxjm/Tz/9tPF1KSV+97vfITw8HJMnT8bVq1dx48aNu9Zvb7ufmjJlirGe1voOAEaOHAlvb29YW1sjKCjIOOUjLCzsjrmct5s3bx5OnDhhfMTGxgJofM/Hjx8Pd3d3WFpa4sknn8T+/fuRnp6OcePGoX///lCr1XjssceM2zp06JDx+fz581vc3549e5CZmYmRI0ciMjISe/bsQV5e3j374IsvvkB0dDSioqKQnZ2NM2fOtNn+2rVr+Om19BsaGvDEE0/glVdeMU4nkE2jrLdrHklty5o1a1BUVITg4GDj6DcAeHh4KHIaDcMrERGZnZQSr7/+ujGU5OTkYMGCBRg6dCgyMzMRFhaG119/HW+++WaL66enpyMuLg5A48fH586dQ0BAAIKCglBZWYkvv/yy1X1bW1sDaAydzV83P29oaGgxMABAQkIC9u/fjwEDBuDpp5/G2rVrW2x3e7ho3r5KpTLOCZVS4l//+pfx2C9dumQMcvb29q32V3tCS0ttNmzYgJKSEmRmZuLEiRPw9PSERqPpdLufur3m1voOwF19ffv70NH5sq3tp639t3e7zzzzjPG9OX/+PJYuXdrmOpcuXcI777yDPXv2ICsrC0lJSffsN1tb27vaLFy4EEOGDMFvfvMb42u+vr64cuUKgMZwW1FR0eIfLi1RqVSYN2/eHT8LGo0Gtra27Vq/J2F4JSIik3N0dERVVZXx+dSpU7F69WrjmdBXr15FcXExioqKYGdnh6eeegqvvvoqjh07dtf62dnZGD58OFQqFQwGAzZt2oSsrCzk5+cjPz8fKSkp+Pzzz1vcb3skJCQgOTkZer0eJSUl2L9/P+Li4lBQUAAPDw+88MILWLBggbE2g8GAzZs3A2iceztmzJg2tz916lR8+OGHqK+vBwBcuHDhjpPMWpKYmIjVq1cbT7a5efPmXW0efPBBbNy4EQDuuApDRUUFPDw8oFarsW/fPhQUFAC4u29aa9cRrfVdVxs1ahS+//57lJaWQq/X4/PPP8e4ceMQFxeH77//Hrdu3UJDQ8Mdwe2BBx4wPm/up5+aNGkSNm/ejOLiYgCN/dzcD2q12vie3a6yshL29vZwdnbGjRs3jCdbtSU4OBg5OTnG53/4wx9QUVGBFStW3NFu9uzZ+OyzzwA0zuGeOHFim3/ESCmN25VS4ptvvsHw4cONyy9cuIDQ0FAAwJYtW/D666/fs9aewNLcBRARUd/j6uqKBx98EKGhoZg+fTqWLVuGs2fPIj4+HkDjiUnr169HTk4OFi1aBAsLC6jVanz44YcAGkelpk+fDm9vbyQlJWHatGkAYBwJHTBggHFfCQkJOHPmDK5du4Znn30Wv/jFL2Bra4tDhw61q9a5c+fi0KFDiIiIgBACb7/9Nry8vPDZZ59h2bJlUKvVcHBwMI682tvbIzs7GzExMXB2dr7jY9qWPP/888jPz0d0dDSklHB3d8fWrVvbXGfatGnGj82trKwwY8aMu85+f/fddzF//ny8++67eOSRR4yvP/nkk5g1axZiY2MRGRlpDDM/fU9ee+21Ftt1RGt9d+7cuQ5vqy3e3t546623MGHCBEgpMWPGDOMJUL/73e8watQo+Pj4ICQkxDj1ZMWKFXjqqafwj3/8A0lJSS1OSQkJCcFf/vIXJCYmwmAwQK1W44MPPsDAgQOxcOFChIeHIzo6+q5LtEVFRWHEiBEIDAzEgw8+eM/6k5KSkJaWhsmTJ6OwsBB//etfMXz4cERHRwMAXn75ZTz//PNYsGABnn76aQwePBj9+/e/I3SPHTsW586dQ3V1NXx9ffHJJ59gypQpeOaZZ1BZWQkpJSIiIow/QwCwb98+vPXWWwAaT2BrPrmrpxMdGVKPjY2V5rqeGhERUUumTJmCtWvXwtvb29ylAGgM3s0jyGR+1dXVcHBwQENDA+bOnYvnnnsOc+fORW1tLWxtbSGEwMaNG/H555+3eEUFU6irq8OECRNw8OBBqFQqk+zzxo0bmD9/Pvbs2QOgcX728uXL75p7ay5CiEwpZWyLyxheiYiIug7Da8/y6quvYvfu3dBoNEhMTMS7774LIQR++OEHvPzyy5BSwsXFBatXr8bgwYPNVufOnTsRHBwMf39/k+zv6NGjUKvViIyMNMn+OorhlYiIiIgUo63wyhO2iIiIiEgxGF6JiIiISDEYXomIiIhIMRheiYiIiEgxGF6JiIiISDEYXomIiIhIMRheiYiIiEgxeHvYVqxateqebV588UUTVNIzffTRR+Yu4a7+72hN5nj/flrjwoULu3wf7fnebUt31ERERNRVOPJKRERERIrB8EpEREREisHwSkRERESKwfBKRERERIrB8EpEREREisHwSkRERESKcc/wKoRYKITIEEJklJSUmKImIiIiIqIW3TO8SilXSSljpZSx7u7upqiJiIiIiKhFnDZARERERIohpJTtbhwbGyszMjK6sRwiIiIi6uuEEJlSytiWlnHklYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GVqAtJKZGVlYXU1FSUlZWZuxwiIqJex9LcBRD1FufOncPWrVtx8uRJaLVa+Pr64o033oCrq6u5SyMiIuo1OPJKdA9SSmi1WpSXl6O6uho/vSudlBJfffUV/vKXv+DEiRMICwvDzJkzcfXqVaxZs8ZMVRMREfVOHHkl+gmNRoMff/wRFy5cwLVr13D9+nXU1dXBYDBACAF3d3eMHj0aiYmJcHR0xMaNG7Flyxb4+flhzJgxsLRs/LGKiorC8ePHkZubi6CgIDMfFRERUe/A8Ep0m8uXL+Nf//oXLl26BCsrKzg5OcHFxQUDBgyAWq2GXq9HUVER/vvf/yItLQ0xMTHYsWMHgoKCEB8fDyGEcVshISE4ceIEdu/ezfBKRETURRheiZoUFRXh7bffRmlpKRITE+Ht7d1iu4iICNy6dQt79+7Ftm3b4O3tfVdwBQC1Wo1BgwYhIyMDzzzzDGxsbExxGERERL0a57wSoXGqwAcffIDi4mLMmjWr1eDarF+/fnjooYcwduxYTJo06a7g2mzw4MGoqKjAyZMnu6NsIiKiPofhlQhASkoKzp07h4kTJ8LJyald61haWiIwMBAqlarVNt7e3hBC4NSpU11VKhERUZ/GaQNdqLS0FIcPH0ZFRQX8/f0RFxcHa2vrbtmXTqdDXl4eLl++jLKyMtTX18PCwgK2trawtbWFnZ0d7O3tYWdnBzs7O6hUKgghIISAwWBAbW2t8eHo6Ijg4GBYWVl1S6093eXLl7Ft2zYMGjQIAwYM6NJtCyHg5uaGnJycLt0uERFRX8Xw2kUOHTqE1atXo6yszBgU3dzckJCQgClTpnTJtT6llDh79iwOHDiAjIwMVFRUQKfTQaVSGUf/dDodpJTG1ywsLIyhtfmh1+thMBiMD5VKBWdnZ0RHR2PkyJGIjIw0njHf2xkMBqxbtw51dXV44IEHumUfHh4eyMnJQU1NDezt7btlH0RERH1F30go3SwjIwMffPABbGxs8Nhjj8He3h43btzA6dOnkZycjB07dmDKlCmYM2dOp8KLlBKZmZlISUnBxYsXodfr4e/vj/DwcHh6et51IpBer4dWq4VWq4VGo4FWq4WU0vgQQsDa2hrW1tawsrJCRUUF8vLykJaWhn379sHT0xMPPfQQJk6cCAuL3j2zJCMjA1lZWd06Su7h4YHs7Gzk5+djxIgR3bIPIiKivoLh9T7dunULq1evhlqtRlJSknHE0svLC15eXqiursaxY8fw5Zdf4tChQwgODoazszNcXFyM/7q6usLT0/OOuZMGgwH5+fk4c+YMMjMzcebMGajVasTExGDIkCFtzrNUqVTG6QLt4eTkBD8/PxgMBhQWFuLkyZP46KOPcPjwYSxYsOCeJy8plcFgwNatW2FlZYVhw4Z12348PDyg0+lQUFDA8EpERHSfGF7v05YtW3Djxg08/PDDLX7U7uDggISEBISEhCAjIwNHjhxBXV0dtFqt8aN9Kysr2NraYtCgQfD390ddXR2ysrJQVlYGrVYLZ2dnxMbGYujQod06EmphYQF/f3/4+fkhJycHhw4dwtKlS/Hcc89h1KhR3bZfc8nIyEBOTg5Gjx7d6tUCuoKVlRXs7OyQl5fXbfsgIiLqKxhe78PVq1eRlpaGoUOH3vMMdTc3N0ybNs34XK/Xo66uDtXV1aisrERpaSkKCgqQlZUFCwsLeHt7Iy4uDn5+fia/PqgQAkOGDIG3tzd27dqFFStWICkpCU888USbI75KYjAYkJKSAmtrawQGBnb7/tzd3XHx4kXjtA0iIiLqHIbX+5CamoqamhpER0d3eF2VSgUHBwc4ODjAy8sLQ4cO7YYK74+DgwMeeughHDlyBFu3bkV1dTVeeOGFXhFgs7KykJubi1GjRpkkTHp4eODs2bMoLy9Hv379un1/REREvVXvPhunG1VUVOCHH35AYGBgr75zkoWFBeLj4xEZGYldu3bh448/hl6vN3dZ92337t0AYJJRVwDw9PSETqdDfn6+SfZHRETUWzG8dtKRI0dQXl6OsLAwc5diEhEREcYAu3r1ahgMBnOXhPLyctTV1XV4vStXruD48eMICQkx2dUU3NzcUF9fz/BKRER0nzhtoJMOHjwIJyenPvURcGRkJKSUSE1NhZ2dHebPn2+W+ZuXL1/Gd999hx9//BEAMGTIEAQHB2PYsGEYPHhwmyPhDQ0N+O6776DRaBASEmKqkqFSqeDo6IiCggKT7ZOIiKg3YnjthKKiIly4cAGhoaHmLsXkoqKioNFo8M0338DGxgZz58412ehlTk4Otm3bhqNHj0Kj0RivvpCXl4fMzEyo1WrY2dlh2LBhCAsLQ1hYGPz9/Y3rHjlyBIcPH0ZxcTGGDh3abdd1bU3znbZ40hYREVHn3TO8CiEWAlgIwBgE+rr09HTU1tZi8ODB5i7FLB544AHodDokJycjLy8Pzz//fLePQG/YsAHbt29HfX09QkJCEBoaekf4rK+vx40bN1BYWIjz588jPT0dtra28PDwgJWVFQoLC6HVajFgwACMGzcOvr6+3VpvS9zd3ZGdnY2ysjK4ubmZfP9ERES9wT3Dq5RyFYBVABAbGyu7vaIeTkqJAwcOwNPTs903AehthBBISEjA+fPncfjwYeTn5+P5559HVFRUt+zvxx9/xLfffgtfX1/Ex8dDrVbf1UatVsPX19cYSmtra3HlyhUUFBSgqqoK0dHRCAwMbHFdU/Hw8MDx48dx+fJlhlciIqJO4rSBDsrNzcWVK1cQGxtr7lLMSgiB4cOHw8vLC3v27MGyZcswY8YMzJkzBw4ODl22n8LCQqxevRr29vYYM2ZMu6coNE8f6M47Z3WUq6sr9Ho98vPzO3V5NSIiImJ47bD09HRotVqTXWKpp3NxccGcOXOQnp6Or776Crt27UJwcDAiIiIQFhaGAQMGdHp+p0ajwcqVK1FZWYlHHnnEZHNru0vzSVuXL182dylERESKxfDaAXq9HgcPHoSfnx+srKzMXU6PoVKpEB8fjyFDhiAvLw9nz55Feno6bGxs4Ofnh5deegkBAQEd2qaUEmvXrsW5c+cwderUXjNFo3///sjNzeVJW0RERJ2k7KEsE8vOzkZJSQmGDBli7lJ6JDc3N8TFxeHRRx/F448/joiICFy+fBl//etfcfLkyQ5tKy0tDXv37kVERAR8fHy6qWLTc3d3R3l5OcrKysxdChERkSIxvHZAeno6DAYDr7rQDvb29ggJCcGcOXNQX1+Pf/7zn0hLS2vXuhcuXMC6devQr1+/bjsJzFw8PDyg0+l4vVciIqJOYnhtJ41Gg8OHD2PgwIFQqVTmLkcx7OzsMHv2bNjZ2WHVqlX48ssvIWXrF624ePEi3nnnHeh0OkycOLHXfbTu6uqKhoYGhlciIqJOYnhtpxMnTqCiogJDhw41dymKo1arMW3aNPj4+CA5ORkrV65s8WPz5uBaW1uLWbNmtXmnLKVSqVRwcHBgeCUiIuoknrDVTkeOHIFKpYK3t7e5S1EkCwsLJCQkICsrC3v27MHx48cxY8YMzJgxAzY2NsbgWlNTg1mzZvWaE7Ra4urqypO2iIiIOonhtR0qKipw/PhxBAUFMWzcByEEIiIiEBQUhMOHD2P9+vX4/vvvMXHiRGzfvr1PBFeAd9oiIiK6H5w20A5HjhxBVVVVj7rgvZI5ODhg8uTJSEpKQr9bigAAACAASURBVHV1NdavX99ngisAeHp6QqvVIicnx9ylEBERKQ5HXu+h+XawTk5O6Nevn7nL6VU8PDzw0EMPoaioCK6urr1yjmtL3NzcIIRAdnY2HnjgAXOXQ0REpCgceb2HgoIC5ObmctS1mwghMGDAgD4TXIHGY/b29sbJkydhMBjMXQ4REZGiMLzew6FDh6DRaDB48GBzl0K9iL+/P8rKynjVASIiog5ieG1DWVkZ9u7dC19fX1hbW5u7HOpF/Pz8UFdXh+zsbHOXQkREpCgMr62QUiI5ORllZWUYNWqUucuhXsbW1hbOzs44deqUuUshIiJSFJ6w1YrMzEwcOHAAERERcHR0NHc51Av5+vri3LlzqKqq6rbvMYPBgMuXL+PSpUvo378/AgIC4Ozs3C37IiIiMgWG1xZUVVVh3bp1sLKyQmRkpLnLoV5q0KBBOHv2LE6cOIGxY8d22XYrKiqQlZWFM2fO4NSpU7h58yY0Gg3UajWsra3h4eGBwMBABAQEwNfXF+7u7nBzc4Nare6yGoiIiLoLw+tPSCnx2WefobCwELNnz4aFBWdWUPdwc3ODlZUVjh492mXhNT09HatXr0ZpaSksLS3h5eWF2NhYeHt7o6qqCsXFxSguLsahQ4ewe/duWFpaGh+urq6IiorC3Llz+WkDERH1WAyvP/Hdd9/hwIEDiIqKQv/+/c1dDvViQggEBgbi5MmTKC8vh4uLS6e3VVNTg02bNiE1NRU2NjaYOXMm+vfvf8cd4RwcHO64vbFWq8WtW7dQWVmJyspK3Lp1CykpKUhPT8ejjz6KhIQE/vFGREQ9DsNrk6qqKhw8eBDJycnw8PBAeHi4uUuiPmD48OE4c+YMDhw4gJkzZ3Z4fSkljh07hvXr1+Py5csYPnw44uLi2hU6ra2t4eXlBS8vL+Nrt27dwg8//ICVK1di3759eOSRRxAWFsbbIhMRUY8hpJTtbhwbGyszMjK6sRzTqq+vx4kTJ3DkyBFkZmaisrISrq6umDZtGuf/kcl88803cHR0xNtvvw2VStXhdT///HOo1WqMGzcO7u7u912PlBJ5eXlIT0+HlBJRUVF45JFHEBQUdN/bJiIiag8hRKaUMralZX1q5LWhoQEFBQW4cOECLl68iOzsbJSXl0MIgUGDBmHcuHFwdXU1d5nUx4SHh2P//v04dOgQxowZ0+71zp07h02bNsHd3R0TJ07sso/4hRAICgpCQEAAzp49i8zMTJw6dQrx8fGYM2cOfHx8umQ/REREndEnwmtBQQG+/fZbnDhxAtXV1dBqtbC3t4enpyciIyMxYMAAfixKZuPv7w87Ozt89dVXCAwMbFc4vHHjBt5//30YDAaMGzeuW+amqlQqhIaGYtiwYTh58iT27t2Lw4cPY+zYsZg5c+Yd82eJiIhMpddOGzAYDDh16hTS0tJw9OhRaLVaYzDw9vaGra2tuUskMiopKcHOnTvRv39/LFiwADExMa3+QXXt2jW88847uHLlCh566CE4OTmZpEaNRoNjx47h/PnzcHBwwOjRo5GYmIhBgwbxjz8iIupSbU0b6FXh1WAw4NKlSzh+/DgOHz6Mq1evwmAwYPjw4QgPD4eVlZW5SyRqVXV1NXbu3AmdTodRo0bhscceg4+PD+rr61FbW4uamhrk5+dj3bp1KCsrw4wZM8wyzaW2thZZWVk4f/48rK2tMWjQIIwePRpRUVHw8vJikCUiovvWK8OrlBIVFRUoLi5Gfn4+Ll26hDNnzqCkpARarRbu7u4ICQnBwIEDebkfUgy9Xo/Tp0/jxIkTsLGxga2tLbRaLQwGAwwGAxoaGmBvb48pU6aY/VqsOp0OOTk5uHjxIm7evAk7Ozu4u7tjxIgR8Pf3x4ABA+Dj4wMXFxcGWiIi6hBFh9cLFy7g8uXLqKmpQXV1NcrLy3Hjxg1cv34ddXV1qK+vh06ng42NDdzc3BAQEICBAwdylJUUTaPR4MyZM9BoNLC2tjY+bG1t4ePj0+P+IKusrERhYSEKCwtRUlICnU4HKysrqNVqODg4oH///ujXrx+cnJxgb28PKysr40OtVt91lQUpJfR6PRoaGu74t/nrhoYGWFhYQK1WG7dx+/aEEG0+btdasG4rcN++rL3tmnl4eNxxeTIiIrpbl4VXIUQJgIKuKqyXcQNQau4i+iD2u+mxz82D/W567HPzYL+bXk/s84FSyhav/9ih8EqtE0JktPYXAnUf9rvpsc/Ng/1ueuxz82C/m57S+rxnffZIRERERNQGhlciIiIiUgyG166zytwF9FHsd9Njn5sH+9302OfmwX43PUX1Oee8EhEREZFicOSViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBTDsiON3dzcZEBAQDeVQkREREQEZGZmlkop3Vta1qHwGhAQgIyMjK6pioiIiIioBUKIgtaWcdoAERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBsMrERERKVZ5eTlWrlzZqXVXrFiB2traO1576623sGHDBixduhRCCOTk5BiXLV++HEII491GZ8yYgfLy8s4X32Tp0qV455137ns7fQXDKxERESlWV4fX1NRUJCYmAgDCwsKwceNG47LNmzcjJCTE+Hz79u1wcXHp1L5NRa/Xm7uELsfwSkRERIq1ePFi5ObmIjIyEosWLcKyZcswcuRIhIeH44033gAA1NTUICkpCREREQgNDUVycjLee+89FBUVYcKECZgwYQIAoLKyEjqdDu7u7gCAOXPmICUlBQCQl5cHZ2dn4zIACAgIQGlpKY4ePYrw8HBoNBrU1NRgxIgROH36NAC0WA8A/PWvf8WwYcMwefJknD9/vs1jzMnJweTJkxEREYHo6Gjk5uYiLS0NM2fONLZ5+eWX8emnnxrrevPNNzFmzBi8/fbbiIuLM7bLz89HeHg4ACAzMxPjxo1DTEwMpk6dimvXrnXqPTA1S3MXQERERNRZ//d//4fTp0/jxIkTSE1NxebNm5Geng4pJWbPno39+/ejpKQEPj4+2LZtGwCgoqICzs7O+Oc//4l9+/bBzc0NALB7925MmjTJuG0nJyf4+fnh9OnTSElJwbx587BmzZq7ahg5ciRmz56NP/zhD6irq8NTTz2F0NBQpKam4uLFi3fVY29vj40bN+L48eNoaGhAdHQ0YmJiWj3GJ598EosXL8bcuXOh0WhgMBhw5cqVNvvFxsYGBw4cAAAkJycjLy8PgYGBSE5OxuOPP476+nr86le/QkpKCtzd3ZGcnIzf//73WL16dYffA1NjeCUiIqJeITU1FampqYiKigIAVFdX4+LFixg7dixeffVVvPbaa5g5cybGjh3b4vo7duzAz3/+8zte+9nPfoaNGzdi586d2LNnT4vhFQCWLFmCkSNHwsbGBu+9916b9VRVVWHu3Lmws7MDAMyePbvVY6qqqsLVq1cxd+5cAI2htD3mzZtn/Prxxx/HF198gcWLFyM5ORnJyck4f/48Tp8+jSlTpgBonF7g7e3drm2bG8MrERER9QpSSrz++ut48cUX71qWmZmJ7du34/XXX0diYiKWLFlyV5v09HR8+OGHd7w2a9YsLFq0CLGxsXBycmp13zdv3kR1dTXq6+uh0Whgb2/faj0rVqyAEKLdx9QSS0tLGAwG43ONRnPHcnt7e+PX8+bNw2OPPYaHH34YQggMGTIEp06dwogRI3Do0KF21dGT3HPOqxBioRAiQwiRUVJSYoqaiIiIiNrF0dERVVVVAICpU6di9erVqK6uBgBcvXoVxcXFKCoqgp2dHZ566im8+uqrOHbs2F3rZmdnY/jw4VCpVHds39bWFn//+9/x+9//vs06Fi5ciD//+c948skn8dprr7VZT0JCArZs2YK6ujpUVVXhm2++aXW7Tk5O8PX1xdatWwEAWq0WtbW1GDhwIM6cOQOtVouKigrs2bOn1W0EBQVBpVLhz3/+s3FEdtiwYSgpKTGG1/r6emRnZ7d5jD3FPUdepZSrAKwCgNjY2JbjPxEREZEZuLq64sEHH0RoaCimT5+O+fPnIz4+HgDg4OCA9evXIycnB4sWLYKFhQXUarVxdHXhwoWYPn06vL29kZSUhGnTprW4j5/97Gdt1rB27VpYWlpi/vz50Ov1GD16NPbu3YvExEScPXv2rnqio6Mxb948REZGYuDAga1OY2i2bt06vPjii1iyZAnUajU2bdqEwMBAPP744wgPD8eQIUOMUxNaM2/ePCxatAiXLl0CAFhZWWHz5s145ZVXUFFRgYaGBvzmN7/BiBEj2txOTyBaG45uSWxsrGy+thkRERFRbzFlyhSsXbtWMfM+ezshRKaUMralZZzzSkRERH3erl27zF0CtRPDKxEREVEP8Mtf/hIHDx6847Vf//rXd10Boa9jeCUiIiLqAT744ANzl6AIvMMWERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBi+VRUQmtWrVKrPt+8UXXzTbvnuKjz76qEPtb++zjqx7r77uaB0LFy7sUHsi6r048kpEREREisHwSkRERESKwfBKRERERIrB8EpEREREisHwSkRERESKwfBKRERERIrB8EpEREREisHwSkRERESKIaSU7W4cGxsrMzIyurEcIiIiIurrhBCZUsrYlpZx5JWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYiIiIgUg+GViIiIiBSD4ZWIiIiIFIPhlYioh2loaEBVVRUMBoO5SyEi6nEszV0AERE1OnfuHNLS0nD8+HHU1NTA1dUVU6ZMwbRp02BpyV/XREQAwysRkdnV1dVhw4YN2Lt3LxoaGuDr6wt/f38UFhbi008/xdmzZ/HKK6/A2tra3KUSEZkdwysRkRlVVFRg+fLlyMrKQlhYGKKioqBSqQAA4eHhuHDhAg4ePAiDwYBf//rXsLGxMXPFRETmxTmvRERmUlVVhWXLluHUqVOYPHkyYmNjjcG12dChQzF27FgcOXIE//nPfyClNFO1REQ9A0deiYjMoKamBv/4xz9w5swZTJ06FT4+Pq22HTx4MGpra7F//34MHToUiYmJJqyUiKhn4cgrEZGJaTQa41SByZMntxlcm4WFhcHV1RVffPEFiouLTVAlEVHPxPBKRGRCUkp88sknOHbsGMaPHw8/P792rSeEQEJCAm7duoV169Zx+gAR9VkMr0REJpSSkoK0tDTExMRg0KBBHVrXwcEBI0eOxNGjR3HkyJFuqpCIqGdjeCUiMpEDBw5g8+bN8PPzQ1hYWKe2ERISAhsbG2zatAk6na6LKyQi6vnuGV6FEAuFEBlCiIySkhJT1ERE1OtkZWVh9erVcHBwQEJCAoQQndqOEALx8fHIz8/Hrl27urhKIqKe757hVUq5SkoZK6WMdXd3N0VNRES9Sl5eHj744ANIKZGYmHjX5bA6ysfHB15eXti6dSvKy8u7qEoiImXgtAEiom50/fp1rFixAjU1NZg2bRrUanWXbHf06NEoLS1FSkpKl2yPiEgpGF6JiLpJeXk5li9fjuvXr2P69Omws7Prsm07OTlh+PDh2LVrF65cudJl2yUi6ukYXomIukFdXR3effdd5OTkYOrUqXBycuryfcTExECj0SA5OZmXziKiPoPhlYioDXV1dSgtLe3QOhcuXMDf/vY3nDx5EhMnToSHh0e31GZtbY3o6GhkZGQgKyurW/ZBRNTT8PawREStyM7OxqpVq1BWVoagoCDExsYiOjoaPj4+LV4toLKyEps3b8aePXug1+sxadIk+Pv7d2uNI0aMQHZ2NjZu3IgRI0bA0pK/1omod+NvOSKin9Dr9di6dSu2bt0KABg6dCgKCgpw6tQpbNy4EQMHDsSgQYPg4uJifFRUVOCrr77C9evXERwcjJiYmC47OastFhYWiI+PR1paGr7//ntMmjSp2/dJRGRODK9ERLcpLS3Fxx9/jGPHjsHf3x9jxoyBpaUloqKiUFtbi7y8POTn56OgoAB1dXUAGgMkAPTr1w9z5sxBv379TFqzv78/+vXrhy+//BJRUVHo37+/SfdPRGRKoiOT/GNjY2VGRkY3lkNEZD4XL17E8uXLUVxcjAcffBCDBw9us72UEjqdDjU1NdDr9XBzc+v0zQfuV0VFBbZs2YL4+Hj89re/NVsdRERdQQiRKaWMbWkZT9giIgJw9epVrFixAhUVFZg7d+49gyvQeLcra2tr9O/fH+7u7mYNjM7OzoiLi0N6ejrvvEVEvRrDKxH1eTdv3sQ///lPlJaWIikpqVsua2UKwcHBcHNzQ3JyMgoLC81dDhFRt2B4JaI+raamBitWrEBBQQGmT58OR0dHc5fUaUIIjB8/HlVVVfj4449RX19v7pKIiLocwysR9Vn19fV4//33kZ2djcmTJ8PV1dXcJd03W1tbJCQk4PTp08arJRAR9SYMr0TUJ2m1WqxatQrp6ekYM2YMBgwYYO6SukzzpbxSUlJ48wIi6nUYXomoz8nNzcWf/vQn7NmzBzExMe06OUtpRo8eDQsLC/zjH//Atm3bYDAYzF0SEVGX4HVeiajPaGhowLfffostW7ZAq9UiMTGxV4243k6tVmPOnDlIS0vD2rVrceHCBSxYsMBkJ6MZDAacOnUKJ0+ehLW1NVxcXODs7AwXFxd4eXnBxcXFJHUQUe/D8EpEfUJhYSE++eQTnDp1Cn5+fhgzZgysrKzMXVa3UqvVmDx5Ms6ePYsDBw7g8uXL+MUvfoFhw4a1uZ5Wq0VFRQVqa2tRU1OD2tpa6HQ6eHl5wc/Pr81+Ky0txYEDB7B//34UFRWhoaEBQONdy1QqFVQqFezt7fHII49g2rRpUKlUXXrMRNT78SYFRNSrGQwG7Nq1C8nJyaiursaYMWMwaNAgc5dlcmVlZdi1axeEEPDx8UFAQAD8/f3h5+cHW1tbXL58GQUFBcjLy0NhYSEaGhpgMBiMDykl1Go1bGxsEBAQgEGDBsHV1RU1NTWorq5GdXU1ysvLcfHiRdTU1MDDwwMhISHw9/eHSqVCfX09ampqUFNTg+zsbFy/fh1hYWF47rnn4Ovra+7uIaIepq2bFDC8ElGvVVxcjDVr1iAzMxPu7u4YN24cbG1tzV2W2dTX1+P06dMoKSnBzZs3UVtbC7VaDQsLC+h0OlhZWaFfv35wd3eHo6MjrK2tYWNjA2tra1haWqKsrAzFxcUoLi5GeXk5dDodLC0tYWVlBSsrK1hbW8PT0xPBwcFwcHBos5bc3FwcOnQIdnZ2eOSRRzB9+nSOwhKREcMrEfUpUkr88MMP2LBhA27evIkHHngAQ4YM4S1Tf6K+vh5lZWXQ6XRwc3ODnZ1du9eVUsJgMNxX4NRoNNi/fz+uX7+O+Ph4/OpXv4Jare709oio92grvHLOKxH1KgaDAevWrcN3330HJycnPPzww/ccBeyr1Go1vLy8OrWuEOK+R0ptbGyQmJiIs2fP4uDBg7CyssJLL73EEVgiahPDKxH1GlqtFh999BH279+P4OBgjBw5kqOtChAcHIyGhgakpaXB1tYWzz33HN83ImoVwysR9Qrl5eV47733cPLkScTGxiI0NNTcJVEHhIWFQafTYefOnbC0tMTTTz8NCwteipyI7sbwSkSKV1BQgH/961/Iz8/HhAkTMHDgQHOXRJ0QExMDvV6Pbdu2Qa/X49lnn2WAJaK7MLwSkaIdPHgQa9asQXV1NWbMmAF3d3dzl0T3IS4uDhYWFtixYwfc3d0xa9Ysc5dERD0MwysRmZ1Op8OVK1dw5coV9OvXDwEBAXB2dm5zHa1Wi+TkZOzYsQMODg54+OGH+/RlsHqT2NhY3Lp1C1999RXCw8M5kk5Ed2B4JSKTk1IiOzsbJ06cQG5uLvLz86HRaIzXDW2+XmhQUBAGDhyIgIAAeHt7w9nZGUIIFBQUYNWqVTh37hyGDx9uHK2j3iMhIQGbN2/Gxx9/jCVLlvT6u6ERUfvxOq9EZFI3btzAxo0bceTIEdTX1xsviu/p6Qk3NzdUV1ejuLgYN27cwK1bt4wX0ler1XB2doafnx/Onz+Puro6jBs3jndn6sUKCwuxa9cuTJ8+Hc8//zyvQEDUh/A6r0RkdhqNBtu3b8e2bdtQUVGByMhIhIaG3nVNT0dHR3h7exufa7Va4x2hbt68iQsXLsDFxQVJSUmwtrY29WGQCfn6+iIiIgJ79uyBr68vpk+fbu6SiKgHYHglom4lpUR6ejo2btyIy5cvw8/PD4mJie2+m5O1tTV8fX05wtpHRUVFoaysDBs2bICTkxMefPBBc5dERGbG8EpE3SY/Px+bNm1CZmYmrK2tMW3atE7f0Yn6JiEEJkyYgB07dmDVqlUQQmD06NHmLouIzIjhlYi63LVr15CSkoKDBw9Cq9UiJiYGwcHBnLNInWJpaYlp06bhu+++w8qVK1FWVoaZM2fy+4moj2J4JaIuU1ZWhm+//RZ79+5FTU0NQkJCEBERwTPF6b5ZWlpixowZSEtLw/r163HlyhU8/fTTcHR0NHdpRGRiDK9EdF+klMjLy8PevXvx448/orKyEkOGDEF0dDSvu0pdSqVSYeLEicjKysKePXtw9uxZzJ07F7GxsXBycjJ3eURkIrxUFhF1iJQSZWVlyM3NRV5eHs6cOYO8vDzodDoEBgYiMjISDg4O5i6TermbN28iLS0NVVVVcHBwQEREBCIiIjB8+HB4e3tzSgGRwrV1qSyGVyJqkZQS1dXVKC0txbVr13D16lVcu3YNeXl5KC0thU6ng8FgQL9+/RAYGIhhw4ZBrVabu2zqY8rKynDx4kVcunQJWq0Wtra2cHV1RVBQEAYMGABfX194enrCzs4OdnZ2sLW1vevybETU8zC8EimElBJ6vR46nQ46nQ719fWor6+/47mFhYXxoVKp7vi3+evm7bT2MBgMqK+vh0ajQV1dHerq6oxfl5eXo7S0FNevX0ddXR0aGhpQX18Pg8EAR0dHuLi4wMvLC56ennB1deWdrajHqKysxNWrV1FUVIRbt26huroaFhYWsLS0vOPnw8rKCg4ODrCzs4O9vT3s7OxgbW0Na2trWFlZGR+tvdbSc2tra/4sEHUh3qSAqAeoqanB3//+dzQ0NECv1xv/bf66+bmUEgaDwfgv0Bhqmx/NH4cKIe74aLS1j0mb/0BtXv/2r2/fj6WlJSwtLWFjYwN7e3u4u7vDyckJTk5O6N+/P5ycnPifM/Vobm5ucHNzQ0REBABAr9ejvLwcFRUV0Gq10Ol0xtsQa7VaVFVVoaysDDqdzvhHXfPPYkNDA4QQxtDb/HXzz93tX9/+2k+D8u1/bAohoFKp7linmYWFBX75y1/C09PTXN1HpBgdGnkVQpQAKOi+chTNDUCpuYvog9jvpsc+Nw/2u+mxz82D/W56PbHPB0op3Vta0KHwSq0TQmS0NrxN3Yf9bnrsc/Ngv5se+9w82O+mp7Q+52eARERERKQYDK9EREREpBgMr11nlbkL6KPY76bHPjcP9rvpsc/Ng/1ueorqc855JSIiIiLF4MgrERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBsMrERERESkGwysRERERKQbDKxEREREpBsMrERERESmGZUcau7m5yYCAgG4qhYiIiIgIyMzMLJVSure0rEPhNSAgABkZGV1TFRERERFRC4QQBa0t47QBIiIiIlIMhlciIiIiUgyGVyIiIiJSDIZXIiIiIlIMhlciIiIiUgyGVyIiIiJSDIZXIiIiIlIMhlciIiIiUox7hlchxEIhRIYQIqOkpMQUNRERERG1S3l5OVauXNmpdVesWIHa2to7XnvrrbewYcMGLF26FEII5OTkGJctX74cQgjjDZtmzJiB8vLyzhffZOnSpXjnnXfuezt9xT3Dq5RylZQyVkoZ6+7e4l26iIiIiMyiq8NramoqEhMTAQBhYWHYuHGjcdnmzZsREhJifL59+3a4uLh0at+motfrzV1Cl+O0ASIiIlKsxYsXIzc3F5GRkVi0aBGWLVuGkSNHIjw8HG+88QYAoKamBklJSYiIiEBoaCiSk5Px3nvvoaioCBMmTMCECRMAAJWVldD9/+3de3iU5Z3/8fc9OZ9JSIBwCCc5mQRICCAiAgIBBFFExVOra7fYrm7r7lW2um6tP213tdrWUs9rQaVuQbAVDwgREAQEMZyDQDmDJBggCTknk8z9+yOTaYAkgJDMDHxe1zVXJs9hnu98ZyCf3LnneaqrqR+su+WWW1i0aBEA+/fvJyYmhoYDed26dePEiRN89dVX9O/fn8rKSsrKykhOTiYnJweg0XoAfv3rX9OnTx/Gjh3L7t27m32Oe/fuZezYsQwYMID09HT27dvHypUrmTx5smebhx9+mDfffNNT11NPPcV1113Hb37zG4YMGeLZ7uDBg/Tv3x+AjRs3MnLkSAYNGsT48ePJy8v7Tq9Bawv0dgEiIiIi39UzzzxDTk4OW7ZsISsri4ULF7JhwwastUyZMoXPP/+c48eP07FjRz7++GMATp06RUxMDL/73e/47LPPiI+PB2DZsmWMGTPG89jR0dF06dKFnJwcFi1axPTp05kzZ85ZNQwePJgpU6bwX//1X1RUVHDvvfeSkpJCVlYWe/bsOaueiIgI5s2bx+bNm6mpqSE9PZ1BgwY1+RzvueceHn30UaZOnUplZSUul4sjR44025fQ0FDW7mcc9QAAIABJREFUrFkDwPz589m/fz89evRg/vz53HHHHTidTv71X/+VRYsWkZCQwPz583n88ceZPXv2Bb8GrU3hVURERC4LWVlZZGVlkZaWBkBpaSl79uxhxIgR/OxnP+PnP/85kydPZsSIEY3uv2TJEv7pn/7ptGV33nkn8+bNY+nSpSxfvrzR8ArwxBNPMHjwYEJDQ5k1a1az9ZSUlDB16lTCw8MBmDJlSpPPqaSkhKNHjzJ16lSgLpSej+nTp3vu33HHHbz77rs8+uijzJ8/n/nz57N7925ycnIYN24cUDe9IDEx8bwe29sUXkVEROSyYK3lscce48EHHzxr3caNG1m8eDGPPfYYmZmZPPHEE2dts2HDBl555ZXTlt10003MnDmTjIwMoqOjmzx2QUEBpaWlOJ1OKisriYiIaLKeF154AWPMeT+nxgQGBuJyuTzfV1ZWnrY+IiLCc3/69Oncfvvt3HrrrRhj6NWrF9u3byc5OZl169adVx2+RHNeRURExG9FRUVRUlICwPjx45k9ezalpaUAHD16lPz8fHJzcwkPD+fee+/lZz/7GZs2bTpr3x07dtC3b18CAgJOe/ywsDCeffZZHn/88WbrmDFjBk8//TT33HMPP//5z5ut5/rrr+dvf/sbFRUVlJSU8OGHHzb5uNHR0XTu3Jn3338fgKqqKsrLy+natStff/01VVVVnDp1iuXLlzf5GD179iQgIICnn37aMyLbp08fjh8/7gmvTqeTHTt2NPscfYVGXkVERMRvtW3bluHDh5OSksLEiRO5++67GTZsGACRkZH8+c9/Zu/evcycOROHw0FQUJBndHXGjBlMnDiRxMREJk2axIQJExo9xp133tlsDW+//TaBgYHcfffd1NbWcu2117JixQoyMzPZuXPnWfWkp6czffp0Bg4cSNeuXZucxlBv7ty5PPjggzzxxBMEBQWxYMECevTowR133EH//v3p1auXZ2pCU6ZPn87MmTM5cOAAAMHBwSxcuJCf/OQnnDp1ipqaGh555BGSk5ObfRxfYJoajm5MRkaGrT+3mYiIiMjlYty4cbz99tt+M+/zcmeM2WitzWhsnUZeRURE5Ir36aefersEOU8KryIiIiI+4KGHHmLt2rWnLfvpT3961hkQrnQKryIiIiI+4KWXXvJ2CX5BZxsQEREREb+h8CoiIiIifkPhVURERET8hsKriIiIiPgNhVcRERER8RsKryIiIiLiN3SqLBG5Yr3++uveLuG8PPjgg94uodW99tprXjnujBkzvHJcETl/GnkVEREREb+h8CoiIiIifkPhVURERET8hsKriIiIiPgNhVcRERER8RsKryIiIiLiNxReRURERMRvKLyKiIiIiN9QeBURERERv2Gstee9cUZGhs3Ozm7BckRERETkSmeM2WitzWhsnUZeRURERMRvKLyKiIiIiN9QeBURERERv6HwKiIiIiJ+Q+FVRERERPyGwquIiIiI+A2FVxERERHxG4HeLkBERC5cWVkZGzdu5NChQzidTtq1a0dKSgpdu3bFGOPt8kREWozCq4iIn1m/fj1z5szh5MmTOBwOAgICqKysJDQ0lHbt2jF8+HAmTJhAdHS0t0sVEbnkFF5FRPzIypUreeONNwgJCWHSpEm0bdsWYwyVlZUcOnSIAwcO8Je//IWsrCxuvvlmJk6cSGCg/qsXkcuHLg8rIuIntm7dyvPPP09kZCTjx48nICCg0e2KiorYsGEDeXl59OzZk2nTpjF48OBWn05QVlbG4cOHcTgchIaGEhoaSkhICJGRkQrUItKs5i4Pq/AqIuIHCgoKePzxxykvL+eWW245r/B3+PBh1q9fT3V1NVdffTXTp0+nb9++F3xsay0VFRWEhobicDT+OV9rLaWlpRw4cIDdu3fz9ddfs3//fiorKwFwOBwYY3A4HERFRZGZmcnYsWOJioq64HpE5PKn8Coi4sestcyaNYtVq1Yxbdq0C5rLaq1lz549fPXVVwQEBDBixAj69OlD27ZtiY+Pp23btgQFBVFYWEhBQQEnT570fD158iTHjx/n+PHjVFRUEBgYSGxsLHFxccTFxREUFERRUREnTpzg5MmTVFVVUVVVhdPpJC4ujk6dOpGYmIgxhurqaqqrq3E6nRw9epSjR48SHx9PZmYm48eP1/xcETnNRYVXY8wMYAZAUlLSoEOHDl36CkVEpElffPEFf/jDHxgwYACpqanf6TFqamrYvHkzO3bswFpLYGCg5+ZwOKipqaGmpoba2lpqa2sJCQkhNDSU8PBwoqOjCQsLo6qqivLycsrKyqisrKS2ttazTWRkJBEREbRt25YOHTqcc2S4oKCAjRs3cuTIERISErj77rsZNWrUd3puInL50ciriIifqqioYObMmVRUVHDzzTdf9LzV+ikAxcXFlJSUUFxcjNPpJCoq6rRba81JLSoqYs2aNZSUlPDAAw8wduzYVjmuiPi25sKrZsyLiPiw5cuXk5uby+TJky/JB66MMYSHhxMeHk6HDh0uQYUXp02bNkycOJFPP/2U2bNnAyjAikizdIUtEREfVV5ezkcffUSHDh1ISEjwdjktJiAggHHjxhEbG8ucOXNYuXKlt0sSER+m8Coi4qOWLVtGfn4+gwcP9nYpLa4+wEZHRzN79mw2bdrk7ZJExEcpvIqI+KCqqiqWLFlCYmIi8fHx3i6nVdQH2ICAAF566SWOHDni7ZJExAcpvIqI+KAvvviC/Px8Bg4c6O1SWlVQUBA33ngjp06d4oMPPvB2OSLigxReRUR8jMvlYsmSJURFRfnEh6paW2hoKFdffTXr1q3j6NGj3i5HRHyMwquIiI/ZunUrBw8eZMCAAd4uxWtSU1OpqKjgk08+8XYpIuJjFF5FRHxMVlYWxhi6d+/u7VK8JjQ0lN69e7N69Wry8/O9XY6I+BCFVxERH3Lw4EG2bt1KSkoKDseV/V/0wIEDKS0tZcmSJd4uRUR8yJX9P6OIiI9ZvHgxVVVV9OvXz9uleF14eDjdu3dn1apVFBYWerscEfERCq8iIj5i7969fPHFFyQnJxMcHOztcnxCWloap06dIisry9uliIiPUHgVEfEB1loWLlyI0+m84k6P1Zzo6Gi6dOnCihUrKCkp8XY5IuIDFF5FRHzA5s2b2bp1K+np6QQFBXm7HJ+SlpZGQUEBy5Yt83YpIuIDFF5FRLyspqaGBQsWEBgYSN++fb1djs+JjY2lQ4cOLFu2jPLycm+XIyJepvAqIuJln3/+Ofv27WPo0KFX/BkGmpKenk5+fj4rV670diki4mX6X1JExIvKy8v529/+RkxMDElJSd4ux2fFx8fTtm1blixZQlVVlbfLEREvUngVEfGipUuXkpuby7BhwzDGeLscnzZo0CDy8vJYvXq1t0sRES9SeBUR8ZKTJ0+yePFiOnXqREJCgrfL8XkdOnQgJiaGjz/+WKOvIlcwhVcRES+w1vLWW29RUFDA0KFDvV2O3xgyZAhHjhxh+fLl3i5FRLxE4VVExAtWrFjBl19+yZAhQ4iOjvZ2OX4jMTGRdu3a8dFHH+m8ryJXKIVXEZFWdvToUf7v//6PuLg4XQb2O7jmmmvIz8/n448/9nYpIuIFCq8iIq2oqqqKl19+mZKSEkaNGqUPaX0HsbGxJCUlkZWVxfHjx71djoi0MoVXEZFWYq3lnXfeYdeuXYwePZrw8HBvl+S3hgwZwqlTp5g7dy41NTXeLkdEWpHCq4hIK/niiy9YtmwZycnJdO7c2dvl+LWIiAiGDh3KF198wZw5c3C5XN4uSURaSaC3CxARuRJ89dVXzJ49m4iICAYNGuTtci4L/fr1o7Kykk8//ZTw8HDuvvtuTcMQuQIovIqItCBrLUuWLOGdd94hJCSEzMxMXQL2EkpLS6OqqooPP/yQ4OBgbrvtNgVYkcucwquISAupqqrirbfeYsWKFSQkJDBmzBgCAgK8XdZlZ+jQoVRXV/Pee+/hcrm44447FGBFLmMKryIiLSA3N5fXXnuNnJwcBgwYwMCBAxWoWogxhhEjRrB27Vree+89ioqKuO+++wgNDfV2aSLSAhReRUSaYK1l3759FBQU0K5dO9q1a9fsGQJqamrYtGkTq1atYuvWrVRVVTFmzBi6dOnSilVfmYwxDB8+nIiICLKysjh8+DAPPfQQHTt29HZpInKJGWvteW+ckZFhs7OzW7AcERHvKywsZP369axZs4YDBw5QXV1NUFAQQUFBxMbG0rFjR2JiYggLCyMsLIzQ0FAqKipYs2YN+fn5BAQE0KdPH/r27avTYXnB0aNHWbFiBW3atOHee+9lxIgRGvUW8TPGmI3W2oxG1ym8iojU2b9/Px988AGbN2+mtLSUNm3a0KdPH9q1a0dxcTFFRUWem9PpxOl0UlNT4znPaKdOnTynwVJY8q6ysjJWrFhBcXExQ4cOZcSIEaSkpGgqgYifUHgVEWlGcXExixYtIisri+rqanr37k2fPn2IiYk5r/2ttbhcLn0Yy8dYa9m2bRvbtm0DIDo6mrS0NNLS0khNTSUqKsrLFYpIUxReRUQaUX+O0I8//pj8/Hx69+7NkCFDCAoK8nZpcgm5XC6OHTvG/v37OXz4MDU1NYSFhdG9e3f69etHnz596N27NxEREd4uVUTcmguv+sCWiFxRysrKyMnJYfv27WzZsoVjx47Rvn17brnlFmJjY71dnrQAh8NBx44d6dixI9ZaTpw4weHDh8nNzWXHjh0EBAQQFhZGjx496Nu3L7169VKYFfFhGnkVkcuatZbc3Fy2b9/Otm3b2LVrF2VlZQAkJiaSnJxMhw4dvFyleEttbS35+fl888035ObmcvLkSQICAggJCaFbt27079+f9PR0rrrqKs1jFmlFmjYgIlcUp9PJ7t272bZtG1u3biU3N5eKigqioqLo1KkT3bp1o127dgojcpb6MJubm8vRo0cpKCggJCSEDh06MGTIEDIyMujZs6feOyItTOFVRC5rNTU1HDhwgD179rBr1y527txJcXExLpeLdu3akZSURFJSEpGRkd4uVfyM0+nk4MGD7N27l2+//Zbg4GA6duxIeno6KSkp9O3bl5CQEG+XKXLZUXgVEZ9iraWkpISCggIKCwupqKigtraW2tpaXC4XQUFBhIaGEhISctqtfllJSQmHDh3i8OHDHDx4kN27d1NWVkZVVRUREREkJiaSlJREp06ddAYAuWSqq6s5dOgQe/fu9ZzPNyIigpSUFJKTk+nVqxdJSUkEBurjJCIXS+FVRFpNTU0NJSUlnDp1ynOrD6knTpwgPz+fEydOUFlZ6Qmr1lrPzRjjuTkcjtPu13/vcrlwOp3U1tYSFRVF+/bt6dSpE4mJiTqPp7SK2tpa8vLyOHz4MN988w1lZWWEhoYSHh5O9+7dad++PfHx8cTHxxMXF0d0dDRRUVFERkbqFyqR86CzDYj4iZ07d3L8+HFPkANOC3aNLWvse5fLdc7tLuaxXC4XVVVVVFdXU1VVRVVVFRUVFZw6dYqSkhJqampwuVzU1NR4RlRDQ0MJCwsjMjKSDh06eH6QR0ZGEhISclo4ra2t9VwEwOl0Ul1dTU1Njed+cHAw7dq1IzY2VkFAvCIgIIDOnTvTuXNnoO4sFseOHePbb7/lwIED7Nixg8rKShwOBwEBAafdoqKiaNOmDWFhYYSEhBAcHOz5GhwcTFBQkOffQ1O3M+fcNvX9+WzXcFl4eDgZGY3mBRGfcUEjr8aY48ChlivHr8UDJ7xdxBVIfW996rl3qO+tTz33DvW99fliz7taaxMaW3FB4VWaZozJbmp4W1qO+t761HPvUN9bn3ruHep76/O3nju8XYCIiIiIyPlSeBURERERv6Hweum87u0CrlDqe+tTz71DfW996rl3qO+tz696rjmvIiIiIuI3NPIqIiIiIn5D4VVERERE/IbCq4iIiIj4DYVXEREREfEbCq8iIiIi4jcUXkVERETEbyi8ioiIiIjfUHgVEREREb+h8CoiIiIifkPhVURERET8hsKriIiIiPgNhVcRERER8RsKryIiIiLiNwIvZOP4+HjbrVu3FipFRERERAQ2btx4wlqb0Ni6Cwqv3bp1Izs7+9JUJSIiIiLSCGPMoabWadqAiIiIiPgNhVcRERER8RsKryIiIiLiN84ZXo0xM4wx2caY7OPHj7dGTSIiIiIijTpneLXWvm6tzbDWZiQkNPqhLxERERGRVqFpAyIiIiLiNxReRURERMRvKLyKiIiIiN9QeBURERERv6HwKiIifq2oqIiXX375O+//wgsvUF5eftqy//mf/+Gdd97hySefxBjD3r17Pet+//vfY4y54CtOPvHEEyxbtuw719mU+++/n4ULFza7zZtvvklubu4FPe7x48cZOnQoaWlprF69+mJKPKeAgAAGDhxISkoKN910E0VFRefc59prr23RmprzXd9z3bp148SJE2ctt9Zyww03UFxczJEjRxg9ejT9+vUjOTmZP/zhD57tCgoKGDduHL169WLcuHEUFhYCcPLkSUaPHk1kZCQPP/zwaY89YcIEBgwYQHJyMj/60Y+ora0F4MUXX2TOnDkX/Bx8gcKriIj4tZYIr1lZWWRmZgKQmprKvHnzPOsWLlzI1VdffUHHqK2t5amnnmLs2LHfuc6L8V3C6/Lly+nbty+bN29mxIgRp62rD0CXSlhYGFu2bCEnJ4e4uDheeumlc+7zxRdfXNIaAGpqas5ru4t9z51p8eLFDBgwgOjoaAIDA/ntb3/Lzp07Wb9+PS+99BJff/01AM888wxjxoxhz549jBkzhmeeeQaA0NBQnn76aZ5//vmzHvvdd99l69at5OTkcPz4cRYsWADAAw88wKxZsy7Zc2hNCq8iIuLXHn30Ufbt28fAgQOZOXMmAM899xyDBw+mf//+/PKXvwSgrKyMSZMmMWDAAFJSUpg/fz6zZs0iNzeX0aNHM3r0aACKi4uprq6m/vSQt9xyC4sWLQJg//79xMTE0PDUkT/+8Y/JyMggOTnZcyyoG2V76qmnuO6661iwYMFpI6TdunXjl7/8Jenp6aSmprJr1y5PjQ888ACDBw8mLS3Nc9yGrLU8/PDDXH311UyaNIn8/HzPuqeeeorBgweTkpLCjBkzsNaycOFCsrOzueeeexg4cCAVFRWNbtfQli1b+I//+A8WL17s2ScyMpInnniCoUOHsm7dOpYvX05aWhqpqak88MADVFVVeZ7bf/7nfzJs2DAyMjLYtGkT48ePp2fPnrz66qvnfD2HDRvG0aNHPd839loCREZGArBy5UpGjhzJHXfcQe/evXn00Ud55513GDJkCKmpqezbtw+oG0meNm0agwcPZvDgwaxduxaAJ598khkzZpCZmcn3v//9s+pp7PiNvecauuWWWxg0aBDJycm8/vrr53zO77zzDjfffDMAiYmJpKenAxAVFUW/fv08/Vi0aBH33XcfAPfddx/vv/8+ABEREVx33XWEhoae9djR0dFAXTCvrq7GGANAeHg43bp1Y8OGDQC8+uqr5/X6+ARr7XnfBg0aZEVERHzJgQMHbHJysuf7pUuX2h/+8IfW5XLZ2tpaO2nSJLtq1Sq7cOFC+8///M+e7YqKiqy11nbt2tUeP37cs/y9996zv/jFL6y11v7yl7+0zz33nJ06dardvn27/dWvfmXffPNNO3LkSPvVV19Za609efKktdbampoaO3LkSLt161bP4z777LOex73vvvvsggULPOtmzZplrbX2pZdesj/4wQ+stdY+9thjdu7cudZaawsLC22vXr1saWnpac/3vffes2PHjrU1NTX26NGjNiYmxvO49bVYa+29995rP/jgA2utPa3e5rZraM6cOfahhx7yfA/Y+fPnW2utraiosJ07d7a7d++21lr7ve99z/7+97/3PLeXX37ZWmvtI488YlNTU21xcbHNz8+3CQkJZx3HWmsjIiI8PbztttvsJ598Yq1t+rVsuM9nn31mY2JibG5urq2srLQdO3a0TzzxhLXW2hdeeMH+9Kc/tdZae9ddd9nVq1dba609dOiQ7du3r7W27jVOT0+35eXlZ9XV1PHPfM+dqb6/5eXlNjk52Z44ccLTm4bvtXpJSUm2uLj4rOUHDhywXbp0sadOnbLWWhsTE3Pa+jZt2pz2/ZmvWb3MzEzbpk0be9ddd9mamhrP8l/96lf2+eefb/J5eBOQbZvIoxp5FRGRy0pWVhZZWVmkpaWRnp7Orl272LNnD6mpqSxbtoyf//znrF69mpiYmEb3X7JkCRMnTjxt2Z133sm8efN4//33mTp16mnr3n33XdLT00lLS2PHjh2eP/ECTJ8+vck6b731VgAGDRrEwYMHPbU/88wzDBw4kFGjRlFZWcnhw4dP2+/zzz/nrrvuIiAggI4dO3LDDTd41n322WcMHTqU1NRUVqxYwY4dOxo99vlu11BAQADTpk0DYPfu3XTv3p3evXsDdaOAn3/+uWfbKVOmAHVTLoYOHUpUVBQJCQmEhoY2Op+1oqKCgQMH0rZtW8+8zvp+NPZanmnw4MEkJiYSEhJCz549T5vyUd/bZcuW8fDDDzNw4ECmTJlCcXExJSUlnnrDwsLOetzzPf6ZZs2axYABA7jmmms4cuTIOfcpKCggKirqtGWlpaVMmzaNF154wTN6+l0tXbqUvLw8qqqqWLFihWd5u3btLng6iS8I9HYBIiIil5K1lscee4wHH3zwrHUbN25k8eLFPPbYY2RmZvLEE0+ctc2GDRt45ZVXTlt20003MXPmTDIyMk4LEgcOHOD555/nq6++IjY2lvvvv5/KykrP+oiIiCbrDAkJAepCYf1cS2st7733Hn369Gn2Odb/6behyspK/uVf/oXs7Gy6dOnCk08+eVotF7rdmUJDQwkICPDU2Zz65+ZwODz3679vbF5p/ZzXU6dOMXnyZF566SV+8pOfNPtaNna8M4/Z8Hgul4t169Y1GlKbep2aOn59IG7MypUrWbZsGevWrSM8PNzzS0hzAgMDcblcOBx1Y4pOp5Np06Zxzz33eH7JAWjfvj15eXkkJiaSl5dHu3btmn3chkJDQ5kyZQqLFi3y/HJQWVnZaD98nUZeRUTEr0VFRXlG0ADGjx/P7NmzKS0tBeDo0aPk5+eTm5tLeHg49957Lz/72c/YtGnTWfvv2LGDvn37ekJavbCwMJ599lkef/zx05YXFxcTERFBTEwM3377LZ988slFPZfx48fzxz/+0RMON2/efNY2119/PfPmzaO2tpa8vDw+++wzAE9Aio+Pp7S09LQzEDR8js1td7769u3LwYMHPWdhmDt3LiNHjrzgxzlTTEwMs2bN4vnnn8fpdDb5Wn4XmZmZvPjii57vt2zZcs59mjr+me+5hk6dOkVsbCzh4eHs2rWL9evXn/M4ffr0Yf/+/UBdYP7BD35Av379+Pd///fTtpsyZQpvvfUWAG+99ZZnnmxTSktLycvLA+rmvC5evJi+fft61v/9738nJSUFqDv7QMP++DKNvIqIiF9r27Ytw4cPJyUlhYkTJ/Lcc8+xc+dOhg0bBtR9sOfPf/4ze/fuZebMmTgcDoKCgjyjqzNmzGDixIkkJiYyadIkJkyY0Ohx7rzzzrOWDRgwgLS0NJKTk+nRowfDhw+/qOfyi1/8gkceeYT+/ftjraVbt2589NFHp20zdepUVqxYQWpqKr179/aExjZt2vDDH/6Q1NRUunXrxuDBgz373H///fzoRz8iLCyMdevWNbnd+QoNDWXOnDncfvvt1NTUMHjwYH70ox9d1HOvl5aWxoABA5g3bx7f+973Gn0tL2TEsd6sWbN46KGH6N+/PzU1NVx//fXn/IBSZmZmo8fv2bPnWe+5ehMmTODVV1+lf//+9OnTh2uuueactU2aNImVK1dy1VVXsXbtWubOnUtqaioDBw4E4L//+7+58cYbefTRR7njjjv405/+RFJSkufMAVD3Qbn6Dxu+//77ZGVl0bZtW6ZMmUJVVRW1tbXccMMNp71Oa9eu9XwIbdeuXRf9/m0t5lxD/w1lZGTYCz2vnYiIiL8YN24cb7/9NomJid4uRa4geXl5fP/73+fTTz9ttWNu3ryZ3/3ud8ydOxeAyZMn89e//pXg4OBWq6E5xpiN1tqMxtZp5FVERMStNcODSL3ExER++MMfUlxcfNEfzjpfJ06c4Omnn/Z8f+YIvy/TyKuIiIiI+JTmRl71gS0RERER8RsKryIiIiLiNxReRURERMRvKLyKiIiIiN9QeBURERERv6HwKiIiIiJ+Q+FVRERERPyGLlIgIiKXpddff93bJTTpwQcf9HYJPue1117zdgmXxIwZM7xdwmVPI68iIiIi4jcUXkVERETEbyi8ioiIiIjfUHgVEREREb+h8CoiIiIifkPhVURERET8hsKriIiIiPgNhVcRERER8RvGWnveG2dkZNjs7OwWLEdERERErnTGmI3W2ozG1mnkVURERET8hsKriIiIiPgNhVcRERER8RsKryIiIiLiNxReRURERMRvKLyKiIiIiN9QeBURERERvxHo7QJERETEd9XW1lJVVUVlZSVVVVUEBwfTtm1bb5clVzCFVxERkctAbW0tOTk5ZGdnU1hYSExMDFFRUURHR9OmTRuMMVRXV1NVVUVNTQ3GGBwOBw6Hg8DAQEJDQ3E6nRw/fpzjx49z7Ngx8vLyqKiowOVyYa3FWovD4WDYsGHcdttttG/f3ttPW65ACq8iIiJ+rKysjBUrVrBs2TKOHTuGy+UiIiLCE1Rra2sJDAz0hM/6Wz1jjCfIWmtxuVyEh4cTGRlJTEwMHTp0ICgoyHMrKipi+fLlfPXVV0ycOJEpU6YQFhbmxQ7IlUbhVURExA+5XC5WrFjBwoULyc/Pp3379lx33XUkJSXhcPzjIy21tbWUlZUBEBjRXRcdAAAWMUlEQVQYSGBgIAEBAZ7HsNZSW1uL0+nE4XAQERGBMabZYycnJ7Nhwwbmz5/Pl19+yfe+9z3S0tJa7smKNGAa/vbV6AbGzABmACQlJQ06dOhQa9QlIiIiTSgvL+eNN95gzZo1xMXFcc0113hlHmp+fj6rVq2itraWO++8k8mTJ58z+IqcD2PMRmttRqPrzhVeG8rIyLDZ2dmXrDARERG5MHl5efzhD39g586dDB06lJSUFK8GxtraWlauXElubi4TJkzg/vvvP23kV+S7aC68atqAiIiIn9iyZQsvv/wyBQUF3HjjjXTs2NHbJREQEMANN9zAxo0bWbx4MTU1NfzgBz/wTE0QudQUXkVERHyctZaPPvqI+fPnExAQwLRp04iMjPR2WR7GGDIyMggICCArK4uAgAAeeOABTSGQFqHwKiIi4sOqqqp44403WLVqFR07dmTkyJEEBvrmj++0tDRcLhdLly4lNjaWW2+91dslyWXIN9/9IiIiQklJCbNmzWLjxo2kp6fTv39/nx/NTE9Pp7S0lAULFhAXF8eoUaO8XZJcZjSjWkRExAfl5+fz61//mo0bNzJq1CgGDBjg88EV6qYQjBgxgpiYGGbPns22bdu8XZJcZhReRUREfExhYSG/+c1v2LdvHxMnTqR79+7eLumCOBwOMjMzMcbw8ssvk5ub6+2S5DKi8CoiIuJDKioqeOGFFzhw4AA33ngjHTp08HZJ30lQUBATJkygoKCAWbNmeS6UIHKxFF5FRER8RE1NDS+//DLbtm1jzJgxxMfHe7ukixIZGcm4cePYu3cvr776KrW1td4uSS4DCq8iIiI+wFrLm2++ydq1axk+fDhdunTxdkmXRIcOHRg6dCgbNmzg3Xff9XY5chlQeBUREfEyay0LFiwgKyuL/v3706dPH2+XdEn17duXHj16sHjxYj7//HNvlyN+TuFVRETEyz788EMWLVpEt27dSE9P93Y5LeKaa64hIiKCt956i7///e/eLkf8mMKriIiIF2VlZfHuu+/Svn17rrvuOr84HdZ34XA4GDt2LNXV1bz44oucOHHC2yWJn1J4FRER8ZJVq1Yxd+5cYmNjGTly5GUbXOuFhISQmZnJsWPH+OMf/0hlZaW3SxI/pPAqIiLiBevXr+eNN94gMjKSMWPG4HBcGT+S27Rpw6hRo8jJyeF///d/qamp8XZJ4meujH8pIiIiPiQ7O5tXX32V0NBQMjMzCQgI8HZJrapLly5kZGSwcuVKXn31VQVYuSCB3i5ARETkSrJp0yZefPFFHA4HEydOJDDwyvxRnJKSgrWWzz77DJfLxY9//GOCgoK8XZb4gSvzX4yIiIgXbN68mT/+8Y84HA4mTZp0xYe11NRUHA4HK1euxOVy8dBDD13xPZFzU3gVERFpBVu2bGHWrFkATJo0ieDgYC9X5BuSk5NxOBysXr0al8vFww8/rN5IsxReRUREWti2bds8wXXy5MkKZ2fo168fxhjWrFlDeXk5P/7xj2nbtq23yxIfpQ9siYiItKCtW7fywgsv4HK5NOLajL59+zJq1Ci2bNnCk08+yc6dO71dkvgohVcREZEWYK1l6dKl/Pa3v8XlcjF58mRCQkK8XZZP6969O1OmTKGoqIhnnnmGpUuXYq31dlniYxReRURELrGqqireeOMN5syZQ3R0NFOmTFFwPU+xsbFMnTqV6Oho5syZw+uvv05ZWZm3yxIfojmvIiIil1BeXh6vvPIKOTk5DBgwgLS0tMv+ylmXWlBQEOPGjWPLli0sXbqUnJwcbr/9dq677ror5mIO0jRzIcPxGRkZNjs7uwXLERER8V/r169nzpw5FBYWMmrUKLp06eLtkvzeyZMnWb16NaWlpSQnJ3PvvffSo0cPb5clLcwYs9Fam9HoOoVXERGRi+N0OvnLX/7C0qVLCQsLY8yYMURGRnq7rMuGtZa9e/fy5ZdfEhwczKhRoxg/fjydO3f2dmnSQhReRUREWkheXh6vv/4627Zto2/fvgwZMkR/2m4hTqeT7Oxsdu3aRVhYGOnp6YwdO9ZzsQO5fCi8ioiIXGLWWtasWcPbb79NUVERo0aNIikpydtlXRGqqqr4+uuv2bFjB9ZaunfvzoQJE7j22mv1wbjLhMKriIjIJVReXs7bb7/NypUriYyMZOzYsURERHi7rCuOy+Vi3759bNu2jZKSEjp37sxNN93EyJEjdZlZP6fwKiIicons37+f1157jT179tC/f3/S09N1NgEfkJeXR3Z2NoWFhXTp0oUpU6YwbNgwjcT6KYVXERGRi+RyuVi6dCnvvvsulZWV3HDDDXTo0MHbZckZjh49yldffUVxcTFxcXEMHz6c4cOH06NHD/2S4UcUXkVERL4jay1bt27lk08+YfPmzbRv356RI0dqRM+HWWs5duwYO3fu5NChQwQHB9OtWzdGjRrFNddcQ5s2bbxdopyDwquIiMgFcjqdrF+/nqVLl7Jv3z4A0tPT6dOnj0bw/IjT6WTv3r3s2rWLkydPEhMTw5AhQxg9ejRXX321XksfpfAqIiJyDk6nkyNHjnDo0CEOHTrEli1byM3NJSwsjAEDBtCzZ0+djsnPFRUVsWvXLvbs2UNAQADdu3dn2LBhpKSk0LVrV72+PkThVUREpIGamhq++eYbT1A9cOAABw8epLKykurqagICAoiPjyc1NZVOnTppdO4yU1tby/79+9mxYweFhYWEhoYSHx9P//79SU5OJjk5mejoaG+XeUVTeBURkSteaWkp27dvZ/PmzWzevJnS0lKcTicOh4M2bdqQkJBA+/btSUhIICIiQoH1ClFZWcmRI0c4fPgweXl5uFwuwsLC6NmzJ7179+aqq66iZ8+exMbG6j3RipoLr4GtXYyIiEhrKS4uZsOGDWzevJmcnBxKS0txOBx07dqV5ORkEhISiIqKUii5goWGhtKrVy969eqFtZaTJ09y+PBhjhw5wvbt2wkICCAkJIS4uDi6du1Kp06d6NixIx07dqRz584EBwd7+ylccTTyKiIilxVrLXv27GHNmjWsXbuWoqIiwsPDSUpKokePHiQkJCisynlxuVycPHmSvLw8jh8/TmFhIcXFxRhjCA4OJiwsjB49enDVVVfRvXt34uLiiIyMJCoqioiICM2hvQgaeRURkcuatZa8vDy2b9/O559/zoEDB3A6nXTr1o2RI0cSGxvr7RLFDzkcDhISEkhISPAss9ZSUlLCyZMn+fbbbz0f7jPGEBAQgMPhwOFwEBAQQFRUFLGxscTFxREbG0tsbCxt2rQ57atG/i+cwquIiPidsrIyvvnmG89cxR07dnDs2DEqKiqIjY0lLS2Nnj176hKhcskZY4iOjiY6Opru3bsDdSO0RUVFlJeXU1FRQUVFBZWVlZSXl1NYWMjRo0epqqqisrLSE3IDAwNxOByEhITQvn17OnfuTIcOHUhMTCQxMZEOHToQGhrq5WfrmxReRUTEZ9XPQTxy5IgnrB46dIhvv/0Wp9NJVVWV58wA/fv3p3PnzvqUuLQ6h8NBXFwccXFxzW5nraWiooLS0lLKysooKSmhtLSUU6dOceTIEc+c7MDAQIKDg4mPjycpKYnExEQ6duxIYmIinTp1Ijw8vJWemW9SeBUREa+pra2lqKiIgoICCgoKKCws9NwvKCggLy+PkpISqqurcTqdREVFERMTQ/fu3YmPjyc+Pv6K/0Eu/sMYQ3h4eJPvWZfLxalTpygsLPTctmzZwurVqwEIDg4mODiYhIQEkpKSaNu2LTExMbRp04aYmBjP/aioqMt6vq3Cq4jIFc7lcuFyuaitrT3rVr+8qfXnulVXV3tu9X82LS0t9cwZLCwsxOl0UlNT4zlGUFAQ4eHhhISEEBkZSdeuXYmPjycuLk7TAOSy5nA4PHNjG7LWUlpa6vml7uTJk2zatMlzXuL6Obb10xGCgoI8I8H1wbb+a1hYGCEhIYSEhBAcHExoaCjBwcGefRtOafDVubgKryIiV5j169czd+5cnE6nJ7hC3Q/I+jPQnPn1zPUN9znz65n3rbWeH4iBgYGeH5rh4eF0796dqKgozye0IyMjdeohkUbUX0jhTLW1tZSXl1NWVua51U9LyM3N5cCBA5SXl1NVVeX5MFl9MK3/Wn8DTrsfHBzMiBEjuOeee1r1uZ7LBZ0qyxhzHDjUcuX4tXjghLeLuAKp761PPfcO9b31qefeob63Pl/seVdrbUJjKy4ovErTjDHZTZ2PTFqO+t761HPvUN9bn3ruHep76/O3nl++s3lFRERE5LKj8CoiIiIifkPh9dJ53dsFXKHU99annnuH+t761HPvUN9bn1/1XHNeRURERMRvaORVRERERPyGwquIiIiI+A2F1yYYY7oYYz4zxuw0xuwwxvzUvfxpY8w2Y8wWY0yWMaZjg30eM8bsNcbsNsaMb7B8kDFmu3vdLOOrl6zwAc30/TljzC537/9mjGnTYB/1/SI00/Pb3d+7jDEZZ+yjnl+kZvoeZ4z51Bizx/01tsE+6vtFMsbMNsbkG2NyGiwbYIxZ5+7hh8aY6Abr1POL1ETPBxpj1rt/lmYbY4Y0WKeeXwJN9H2+u+dbjDEHjTFbGqzzn743vAKKbv+4AYlAuvt+FPB34GogusE2PwFedd+/GtgKhADdgX1AgHvdBmAYYIBPgInefn6+emum75lAoHv5s8Cz6nuL97wf0AdYCWQ02F49b9m+/wZ41L38Ub3XL3nfrwfSgZwGy74CRrrvPwA8rZ63eM+z6nsG3AisVM9bvu9nrP8t8IQ/9l0jr02w1uZZaze575cAO4FO1triBptFAPWfeLsZmGetrbLWHgD2AkOMMYnUBd51tu5d8DZwS6s9ET/TTN+zrLU17s3WA53d99X3i9RMz3daa3c3sot6fgk01Xfq+vuWe7O3+EcP1fdLwFr7OVBwxuI+wOfu+58C09z31fNLoImeW6B+hDsGyHXfV88vkSb6DoB79PQO4C/uRX7V90BvF+APjDHdgDTgS/f3vwa+D5wCRrs360RdqKr3jXuZ033/zOVyDmf2vYEHgPnu++r7JdRMzxtSzy+xM/re3lqbB3UB1xjTzr2Z+t5ycoApwCLgdqCLe7l63nIeAZYaY56nbgrjte7l6nnrGAF8a63d4/7er/qukddzMMZEAu8Bj9SPulprH7fWdgHeAR6u37SR3W0zy6UZjfXdvfxxoIa63oP6fsk01fPGNm1kmXr+HanvPuEB4CFjzEbqpnBUu5er5y3nx8C/uX+W/hvwJ/dy9bx13MU/Rl3Bz/qu8NoMY0wQdT9U3rHW/rWRTf6Pf/x56Rv+8ds61P1ZO9e9vHMjy6UJTfXdGHMfMBm4x/3nC1DfL4nzeK83pJ5fIk30/Vv3n+pwf813L1ffW4i1dpe1NtNaO4i6H+j73KvU85ZzH1D/nl8A1H9gSz1vYcaYQOBW/vEXTPCzviu8NsE9H+RPwE5r7e8aLO/VYLMpwC73/Q+AO40xIcaY7kAvYIP7z38lxphr3I/5fer+NCWNaKbvE4CfA1OsteUNdlHfL1JTPW+Gen4JNNP3D6j7wY7766IGy9X3FlA/NcMY4wD+C3jVvUo9bzm5wEj3/RuA+j9fq+ctbyywy1rbcDqAf/Xd258Y89UbcB11Q+PbgC3u243UjZLkuJd/SN0HW+r3eZy639h30+DTeECGe599wIu4r2ym2wX1fS9wpMGyV9X3Fu/5VOp+664CvgWWquet0ve2wHLqfpgvB+LU90va978AefxjLt8PgJ9Sd7aHvwPPNOyfet5iPb8O2EjdJ9y/BAap5y3fd/fyN4EfNbK93/Rdl4cVEREREb+haQMiIiIi4jcUXkVERETEbyi8ioiIiIjfUHgVEREREb+h8CoiIiIifkPhVUR8njGmgzFmnjFmnzHma2PMYmNM73PsU+r+2tEYs9B9/35jzIsXWcsjxpjwC9xnlDHmo4s57oVy1/l99/03jTEHjDFbjDG7jDG/bLDdSmPMbmPMVmPMWmNMH/fyeWec11pExCcovIqIT3OfGPtvwEprbU9r7dXAfwLtz2d/a22utfa2Czme+2T1TXkEuKDw2trcV9B5gLqrANabaa0dCAwE7nOfiLzePdbaAcBbwHPuZa8A/9Ea9YqIXAiFVxHxdaMBp7W2/qpHWGu3WGtXG2MijTHLjTGbjDHbjTE3n7mzMaabMSanwaIuxpgl7tHGXzbYZqcx5mVgk3ubV4wx2caYHcaY/+fe7idAR+AzY8xn7mWZxph17hoWGGMi3csnuEc511B3KcazuI+72r3vJmPMte7lDmPMy+5jf+Qeab7NvW6QMWaVMWajMWapcV9K9gw3AJustTWNrAt1fy1rZN3nwFXu+6uBse4gLCLiMxReRcTXpVB3JZ7GVAJTrbXp1IXc37pHapszBLiHuhHI240xGe7lfYC3rbVp1tpDwOPW2gygPzDSGNPfWjuLustajrbWjjbGxFN3OdGx7hqygX83xoQC/wvcBIwAOjRRSz4wzr3vdGCWe/mtQDcgFfhnYBiAMSYI+CNwm7V2EDAb+HUjjzucs3v2nDFmC3VX2plnrc1vZL+bgO0A1loXdVe2G9BE7SIiXqHfqEXEnxngv40x1wMuoBN10wmONbPPp9bakwDGmL9Sd5nK94FD1tr1Dba7wxgzg7r/JxOBq6m7lGtD17iXr3Vn5mBgHdAXOGCt3eM+zp+BGY3UEgS8aIwZCNQC9fN4rwMWuAPksfpRXuoCdgrwqft4AdRd/vFMicDOM5bNtNYudI8MLzfGXGut/cK97h1jTAVwEPjXBvvkUzfS3NQvDyIirU7hVUR83Q6gqTmr9wAJ1F0X3WmMOcg//izelDOviV3/vefP6O75oD8DBltrC40xbzbxuIa6MHzXaQvrwuj5XHv734BvqRvddFA3klz/uI0xwA5r7bBzPG5FE/VirS01xqykLiDXh9d7rLXZjWwe6n4sERGfoWkDIuLrVgAhxpgf1i8wxgw2xowEYoB8d3AdDXQ9j8cbZ4yJM8aEAbcAaxvZJpq6MHvKGNMemNhgXQkQ5b6/HhhujLnKXVe4+ywIu4Duxpie7u1OC7cNxAB57hHW71E3kgqwBpjmnvvaHhjlXr4bSDDGeKYRGGOSG3ncnfxj7upp3HNYhwL7mqipod7U/fIgIuIzFF5FxKdZay0wlbrQuc8YswN4krq5p+8AGcaYbOpGYXedx0OuAeYCW4D3GhtxtNZuBTZTF9xmc3rAfR34xBjzmbX2OHA/8BdjzDbqwmxfa20lddMEPnZ/YOtQE7W8TN0n/9dTFxTrR3/fo25uag7wGvAlcMpaW03dKPSzxpit7udwbSOP+wlw/RnL6ue8bqNuXutfm6gJAHdorrDWNjYtQUTEa0zdzwUREfElxphI95/42wIbgOHW2ubm8p65/9+A/6ifd/sdjv9vQLG19k/fZX8RkZaiOa8iIr7pI2NMG+o+BPb0hQRXt0ep++DWdwqvQBF1I9QiIj5FI68iIiIi4jc051VERERE/IbCq4iIiIj4DYVXEREREfEbCq8iIiIi4jcUXkVERETEb/x/7rh0C+gL3bMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x576 with 4 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "stacked_plot([cal1, cal2, cal4, cal3])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}