GS_plot_traj.ipynb 75 KB
Newer Older
Patrick Wagner's avatar
Patrick Wagner committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "### "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Once deleted, variables cannot be recovered. Proceed (y/[n])? y\n"
     ]
    }
   ],
   "source": [
    "### %matplotlib inline\n",
    "%reset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [],
   "source": [
    "## Load required modules\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import xarray as xr\n",
    "from mpl_toolkits.basemap import Basemap\n",
    "import matplotlib.lines as mlines"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [],
   "source": [
    "## Load data\n",
    "\n",
    "path=\"/gfs1/work/shkpwagn/ARIANE/VIKING20-K301_Turtle/\"\n",
    "figdir=\"/gfs1/work/shkpwagn/ARIANE/VIKING20-K301_Turtle/FIGURES/\"\n",
    "out=figdir+\"gs_traj_050607\"\n",
    "d=1\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/h/shkifmwr/TM/software/miniconda3_20180110/envs/py3_std/lib/python3.5/site-packages/mpl_toolkits/basemap/__init__.py:4995: MatplotlibDeprecationWarning: The is_scalar function was deprecated in version 2.1.\n",
      "  elif masked and is_scalar(masked):\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<xarray.Dataset>\n",
      "Dimensions:     (nb_output: 147, ntraj: 284)\n",
      "Dimensions without coordinates: nb_output, ntraj\n",
      "Data variables:\n",
      "    traj_lat    (nb_output, ntraj) float64 ...\n",
      "    traj_lon    (nb_output, ntraj) float64 ...\n",
      "    traj_depth  (nb_output, ntraj) float64 ...\n",
      "    traj_temp   (nb_output, ntraj) float64 ...\n",
      "    init_t      (ntraj) float64 ...\n",
      "    ts_gb       (ntraj) float64 ...\n",
      "Attributes:\n",
      "    history:                    Fri Oct  6 13:25:50 2017: ncks -A -v ts_gb /g...\n",
      "    history_of_appended_files:  Fri Oct  6 13:25:50 2017: Appended file /gfs1...\n",
      "    NCO:                        4.6.7\n",
      "<xarray.Dataset>\n",
      "Dimensions:     (nb_output: 147, ntraj: 393)\n",
      "Dimensions without coordinates: nb_output, ntraj\n",
      "Data variables:\n",
      "    traj_lat    (nb_output, ntraj) float64 ...\n",
      "    traj_lon    (nb_output, ntraj) float64 ...\n",
      "    traj_depth  (nb_output, ntraj) float64 ...\n",
      "    traj_temp   (nb_output, ntraj) float64 ...\n",
      "    init_t      (ntraj) float64 ...\n",
      "    ts_gb       (ntraj) float64 ...\n",
      "Attributes:\n",
      "    history:                    Fri Oct  6 14:05:54 2017: ncks -A -v ts_gb /g...\n",
      "    history_of_appended_files:  Fri Oct  6 14:05:54 2017: Appended file /gfs1...\n",
      "    NCO:                        4.6.7\n",
      "<xarray.Dataset>\n",
      "Dimensions:     (nb_output: 147, ntraj: 427)\n",
      "Dimensions without coordinates: nb_output, ntraj\n",
      "Data variables:\n",
      "    traj_lat    (nb_output, ntraj) float64 ...\n",
      "    traj_lon    (nb_output, ntraj) float64 ...\n",
      "    traj_depth  (nb_output, ntraj) float64 ...\n",
      "    traj_temp   (nb_output, ntraj) float64 ...\n",
      "    init_t      (ntraj) float64 ...\n",
      "    ts_gb       (ntraj) float64 ...\n",
      "Attributes:\n",
      "    history:                    Fri Oct  6 14:05:54 2017: ncks -A -v ts_gb /g...\n",
      "    history_of_appended_files:  Fri Oct  6 14:05:54 2017: Appended file /gfs1...\n",
      "    NCO:                        4.6.7\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVIAAAD7CAYAAAA4sRx9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzsnXd0VOXWxn9TMum99xBaaCFAQKQTBAQU6YoKggroBS7KRewQ7Cj3ExUUUaSogAIKSBHBhBYQCIQWQhIICem9l6nn+2OcY4YUEpKQwM2zFiuc/p6ZM/s87y7PlgiCQCta0YpWtOLOIW3uAbSiFa1oxb2OVkPaila0ohUNRKshbUUrWtGKBqLVkLaiFa1oRQPRakhb0YpWtKKBaDWkrWhFK1rRQLQa0la0ohWtaCBaDWkrWtGKVjQQrYa0Fa1oRSsaCkEQWvQ/iUQidO/eXZDL5fX6K5FIBCsrK0EulwsdOnQQAEEikQgymUzo3r27IJVKBXt7e0EqlQrOzs4CIJiYmAhSqVTo3r27IJPJBHt7ewEQXFxcBKlUKjg4OIjHy+Vywc3NTZBIJAIgSKVSwdXVtVH/yuVyARAcHR0FiUQimJiYCBKJxGg/iUQitG/fvt6fz936K5VKBalUWu24PTw8jPa3trYWAMHBwUGQSCSCjY1Nk3yu98pfw3NX03bDc3sn5wcEU1NT8Xlv7uekur/dunVr9nEA2XWxU/Jms+B1hIuLC+fPn6/XMdnZ2UycOJHXX3+d7du3c/PmTSZPnsyiRYuws7MD4MiRI0ycOBEHBweys7OxsrLCxMSE/Px8pk2bxoIFC4iKimL69OlkZ2djb29PeXk5a9as4fnnnwdAIpGI12zbti2TJ09GoVA0yn2Hhoai0+mwsbHB09OT3NxczM3NGTFiBF27dhX32717N3PmzGH27NmNct3GRGpqKv7+/lhbW1NeXs4LL7wgfmYFBQVs3bqVrl27snTpUsaOHcuVK1fIyspCrVbzyCOP0LVrV2QyGcOGDWvmO2kenD59mpKSEkJCQqrdrtVqSUpKwt/fv87nLCgoYOXKlQAolUpiY2Pp0KFDo4y3sbF69WoyMzN55513mm0MEokkqU77tfRa+549ewrnzp2r1zE3b97E19cXAF9fX/r378/Zs2fJyspi/fr1BAYGUlJSQm5uLl988QU5OTls27aNoqIi2rdvj52dHbm5uUilUsrKyoiMjMTFxQVXV1fs7e3F6xiMgrW1NcXFxcyePRs3Nzek0oZ7TC5cuICnpyerVq0CwN3dnTlz5gBQUlJCWFgYGo0GjUbDtGnTeO211xp8zcZGcnIyPj4+9OzZkzFjxiCTyYy2C4JAbGws+/fvZ9u2bQwfPhxBEAgKCuLixYt4enoyffp0TE1Nm+kOmhdqtRqdTtco9y8IAitXrqS8vByVSgVAUlISPj4+DT53U6GsrAxBELC0tGy2MUgkkrOCIATfdr+Wbkjd3NyEjIyMeh/XtWtXunfvjre3N+bm5gBcu3aNP//8k/LychQKBVlZWQCEh4czZMgQACIjI8nJyeHhhx++7TViYmLIysoSjwU9M502bVq9x1sd1q9fT1JSEkuWLEGn03HkyBGys7MRBIHJkyczZswYDh06REhICP369WuUazYmiouLsbGxAaB///60adOGdu3aoVKpjJh7eHg4wcHBrFixgszMTNzc3FAoFAwePJiKiopWRloDI60rysvLycvL48qVK0RERHDgwAECAgJatBGFFsNI62RIW/zU3sPDo97HCIKAi4sLJSUlohEFaNeuHe3atWPz5s3ExcWJ69u3by/+Pzj4tp+ZiE6dOuHm5sbq1auZO3cuoGfDjQVXV1eSkpKqPEjm5uZMnDgRNzc3nn32WVxdXRvtmo0Ja2trli1bxubNm4mIiCAiIkLcFhoaCkB8fDw2Njbiy8jgxlGpVGi1WgYMGHC3h91i0KNHD3Q6XYPPk5uby7fffgtAnz59GDFiRIPPeTcwc+ZMWjrRM+CeZKRqtZqUlBQqKiqqPUYQBG7evImZmVm10yK1Wk1ZWZm47O7uXmffplarRaPRVDlvWVkZOTk5yOVyLCws6nSu20GpVKJSqUhJSWHt2rX06tULW1tbTp8+zdmzZ8WpX+/evdmxYwdOTk6Nct3GhlKppG3btqSmpgJgZmbGrFmz+P333xEEgU8//ZRHHnkEgPfee4+3334bgDZt2uDl5dXKSBvASM+cOcOFCxdIS0ujf//+/Pjjj3h7ezfiKJsOrYy0EWFgpAUFBZiZmWFmZkZKSgrW1tb4+fkZBXwqo6ysDGtra6ytrWs8t06nIyMjA51Oh5eXF5aWlrX6N0tLS4mJiQGgY8eORmzXMMa8vDxUKhVWVlYN9pUKgoBOp8PPz4/y8nLmzZsHgKmpKaNGjSIoKAidTseff/5Jt27dSE9Pb9D1mgoKhYIRI0YQGRlJz549iY6OZvPmzSxevJgFCxYYvZRcXFzo2rUrly9fJjg4mHbt2jXjyJsXDWGkWq2Wa9euERkZSWZmJl999RXjxo3Dzc2tkUfZdLiXGGmLzyNNS0tjw4YNODs74+joSGxsLBUVFTg6OtZoRAHkcjnFxcUolcoa95FKpVhZWaHVaklISODChQtGTPVWGLZZWlpiYmJSZbudnR1t2rTB3NycwsLCetxlzdcrLS3FwsLCKDKrVCrZuXMn27Zt49SpU1hZWZGRkUFJSUmDr9kUkEgkfP311/To0YPt27eLhnLx4sVGRlSpVDJnzhwuX74MwLZt2/jll1+aa9jNjqioKCN3SH1w48YNtmzZQmZmJu+++y4vvPDCPWVEQR8jWL58eXMPo05o8VP7oKAgIScnh1GjRpGZmUlcXBw7duygS5cutR6nVCq5dOkSAA4ODpiZmd32Wvn5+djb21frl83Pz+fGjRvodDq6du1a5XyGKb+JiQkZGRnk5eXh4OBQjzutCkOOmlQq5dq1awwePBgAJycnZDIZISEhODs7o9VqUSgUfPjhhy0+wv13bnC129RqdRUXy+OPP06nTp3uxtBaHO40ah8VFcWuXbsAWLlyJQsWLGiK4TU57qWofYuf2l+/fh0fHx+8vLzw8vIiIyODnJycWn+QoJ/+Ojs7k52dTV5eHu7u7jXur9PpKCwsRKfTGaU3GVBaWsqNGzewsLBAIpFUa5TT09PJzMzUJ+fK5Wg0GkpKSqioqMDe3r5K6k9dUFZWhlarxcbGRhx727ZtuX79OqCPxn722Wf1Pm9zorbv7Nq1a1XWxcXF/c8a0qioqDvyke7evRuAhQsXikHQexHr169vdh9pXdHiDalMJjNKQH/44YdRqVRkZ2fj4uJS67HZ2dmAPvpd2w+4cjCrOoNXWlqKRCKhoqICPz8/QJ8jOX36dDIyMpBKpUyaNImnn34apVLJ7NmzSU1Nxc3NjQ8//BCVSoWrqyuhoaGEhYVhbm7Op59+Srdu3QDw9vYmICAAAE9PTzZs2ACAhYUFgiCIU3wPDw/RiIJ+6ns/4caNGwA8+OCDnDx5EktLy3smwtwUuBMfqSE9DmDFihW1PvctHa0+0kZEWVmZEbWXyWSYm5uTkpKCRqOp8ZjIyEhxWa1W13oNV1dXMeJdVFRUZbuTkxNubm506tRJrIySy+X897//JSYmhr/++otNmzYRFxfHypUrefDBB4mLi2P06NFs3LgRgN9//53r169z7Ngxli9fzuuvvy6e38zMjIMHD3Lw4EHRiBruo6SkRHyY0tLSjMb1+OOP13pf9xoMhvTkyZO88MILDBkyhBMnTjTzqJoPdfWRCoJAamoqBw4c4Msvv+Ttt9+mvLz8njaicG/5SFs8I5VIJFV8JDKZDLlcTmlpKba2tlWOMUT4KyoqsLa2vu20WiaTIZPJsLKyory8vMp2qVSKu7u70Tp3d3dxnbW1NR06dCA1NZXw8HCOHTuGqakps2fPZuDAgSxcuJADBw4wfPhw8vPz6dWrF4WFhWRmZtaaA2pgpFKpFFtbW2bPns3atWsBCAkJYcuWLbXeV0uHIAg88MADnDlzhnXr1hkVNlRUVIhZCfcrwsLCCAwMrDFtrTZGmp2djVQq5dKlS+Tm5nLp0iVkMhnR0dH3jSvkXmKkLd6QqtXqap3NcrkclUrFSy+9VKUWv6KiAkEQbmtAdTodarVaTFMSBAELCwt69uwp1iPXBYmJicTExPDEE0+wePFiMVjl7u5Obm4ubm5uotE0RPvd3d3JyMjA1dUVpVLJqFGjkMlkzJs3T6yqquwjBYyMfFhYGKampoSEhDBnzhxCQkJabB5pTaicHvbcc8/RuXNncfnixYtkZ2dTVFR0X+aRFhcXc/ToUeLj45kyZUq1vvnqfKTl5eWEh4dz+vRpcZ2/vz/R0dF4enpWSyzuVbT6SBsRgiBUGyGXSCQ1pjbpdLo65XAapvw6nQ65XI6pqWm9cz9LSkqYOHEiK1euFA2eVqs1MuIuLi5ikn7lezFMvU6fPo2bmxtJSUlMmTKFgIAA/Pz8REZqQNu2bcWKoNDQUNzd3QkLCyMsLAwTExPS09NxdHSs1/ibA4Ig8NVXXwHQvXt3+vXrx65du7hy5Yq4z6BBg7C0tLxvGamlpSUWFhakp6dz+PBhxo8fX2UfKysrUlJS0Gq1ODk5ceTIESIiIsjMzGTQoEGUlJRw9OjRZo1qNyVaGWkjwtTU1MjXU1RURElJCZaWlpiZmYnMMS0tDSsrK2xsbLh48WKdIuXl5eXk5+cDEBgYWG/lJrVazcSJE3nqqacYMWIEiYmJ2NnZcejQIdq1a4dGo8HFxQWJRIKvry8ajYacnBzs7OxIT08Xp/WG/D5fX18efPBBLl++jJ+fXxVGWhndunUT07sAPvroo3vCiH733Xc899xz4rJKpcLW1hZHR0fRB/z6669jamrKmTNn7ltGKpVKmTlzJqtXr+bChQtYWFjQo0cPXFxc+P333ykqKuLGjRtUVFRQXFws+o9XrlzJs88+W2uhyf2CVkbaiLh1ypObmwvoWZ9hGlNcXExGRgaCIBAYGIhKpaoTszQzM0Mmk6HVauvNfARB4LnnnqNTp04sWLCAS5cuYWFhwciRIzl8+DD+/v6sXr2aMWPGUFxczIABA9iyZQtDhgwhMjISMzMzXF1dKSgowNzcHFNTU/Ly8jhz5gz/+te/AKow0sro3r27aEjlcjn/+c9/6NChg+giaKm41djHxMRw8+ZNBg0aJN7PtWvX6NKly33vI83JyRH/f/LkSU6ePAnoyYNhtrVu3Tp+/vlnxowZwzvvvFOtC+B+RSsjbUQUFxcbLet0OszMzOjRo4doLA16omVlZeTn5yOXy+sUsZRIJCgUCsrLy+ud9BwREcH3339Pt27d6N69OxqNhjfffJP58+fzwgsvsHXrVlxcXFi3bh2xsbF06tQJT09PxowZg4WFBa+++iqCIBAfH89rr72GRCJBEATmzZsn6kPWxkjbtWtHaGgoRUVFIit/9NFH6dWrF5s2baJ9+/bs27cPBwcHBg4cWK97a0potdoq6+RyOQUFBeLytm3b6NSpE+fPn79vGSnov8OAgADi4uIIDAwUff0GIxoQEEBUVBS///57cw6z2XAvMdIWX9nk7u4uvPDCC+JyYmIi06ZNM3LAZ2Zmkpubi1qtRqPRYGdnV6UOviao1Wry8/MJCgq6o/EJgkB0dDSWlpZVjHFWVhaWlpaoVCpMTU0pKytDrVbj7OxMVlbWbbUAKlc2JSUl8ccff1S7n8FvWhP27NnDmDFj6n1vTYHy8vI6i7rMnj0bR0fHFl+tdadQq9UcOXKE48ePU1xcjLm5OWlpaSgUCuzt7dFqtWi1WqysrJp7qM2Ce6myqcXnkd7KSL28vNBoNKI4LeiZpVarRa1WY2FhUWcjCno2JAiCESOqD4qKipBIJNX+2GUyGfb29qjVamQyGba2tlhbW6PRaHBwcKC0tLTWqYshj/R2WLJkSa3bT5482SKmyBkZGfVSxtq5cyfHjx9vwhE1H5RKJWvWrOH48eP06dMHKysrZDIZ3t7euLq6olAo7qk8yqbAvXT/9xwjBRg4cCBdunQRgzVqtVoUM7kThXqVSkVeXh5dunSpN/tJSEhAq9VW+9bMz8/H3d0drVZLTk6O6It1dnbGzc2N2NhYTExMajT8dWWkhmsdOHCAq1ev1rhPc3/XUVFR9OzZs9pthnJeA3r06MGgQYOwsLC47xhpXFwcYWFhPPLII6xatQq5XF6tX7u8vLyVkbYy0sbBrYwU9LJs6enpoh6piYkJPj4+mJiY3JF0nUKhwMrKioSEhHozN4VCUeMxEokEjUaDs7MznTp1omvXrgQGBop1/zY2NrVWXdWVkYI+KPfEE08A0KtXL3F97969xf9HRUXV6VxNhejoaKPl/v37s2DBAvz8/MjOzsbZ2VncNnbsWK5du3bfMdKKigo2b95MRkYGn332GaampjUGB+8lRtYUuJfuv8UHm6rzIcpkMiwtLbl+/ToBAQGiwWpISZyVlRV5eXmkpaXh5eVVr+Oqk97T6XRiak9NsLCwELMQatpeXxb58ssvY21tzeDBg8XPacyYMYSGhvLpp5+yadOmep2vMWEQdgZ48skn6dChA1qtlilTpmBubo5EIiE3N5fw8HB0Ot19GbXPysrCz8+PuLi4aqUYK2PmzJnVBuf+V9AatW9EVMdIATGIk5aWRklJCaWlpQ2+lpmZWa36pTWNo7S0FEtLS9GQa7VaCgsLcXBwqHVaam5ujlqtrlHJqraofU0wGG7DMVu2bCE2Nhb4h5FevXqV7Ozsux7Nv379Or6+vnh4eBATE8PmzZuNtvv5+TFjxgwmTZoEwLlz5+67qP2+ffswMzNDLr/9T2/9+vWkp6fz7rvv3oWRtTzcS1H7Fm9Ia4tqm5iYkJmZCfwzPW8IVCpVvZPaTUxMMDExMQpWqdVqHB0db9tvyjDmwsJCUQylMu6Ekd6Kyue9fPkygYGBXLp0CVdXV27evNlo7aPrgpSUFJKSkkhK+qfD7ejRo2nbti3R0dFVKtjuN0aalZVFRkYGTz31VJ1mT62MtJWRNhpqYqRgnJMol8vrJN5cG1QqVZ0rRirL6EkkEmbOnMncuXPJy8tjxowZJCUl4efnx88//4y9vT2CILBgwQL27duHhYUFGzZsoGPHjqIf9MUXXyQtLQ2JRML333+Pt7f3HTHSW2HolAr6+n5D0vuQIUOQyWRkZWWRmZkpSvo1Jfbv3y/+/4033jAy4oMGDaqy//2WR2pQFrtdupoBrYz03mGkLT7YVJths7a2xsHBATc3t2oZnVarpaCgQEwzquntplKpxAhhXVOnKsvonTp1im+++YbExEQ+//xzHnroIeLj4xk2bBgfffQRoDci8fHxxMfHs3btWmbNmkVcXByWlpa88sorvPjiixw5coS9e/eK4iMWFhYNZtmGiOeSJUuMfLk+Pj689957BAYGEhgYKMr93S0UFxcTExNjZFxBX+kUGhpKWloaQUFB900X0YsXL/LDDz/w+uuv17kP1cyZM3n11VebeGQtF/fS/d/TjFQqldbKQg3TftCzAUEQqlXKLy4uxsHBAS8vrzoHrG6V0evUqROpqans2rWLw4cPA/DMM88wZMgQli9fzq5du5g+fToSiYS+ffuSl5eHUqmksLAQrVYrMrLKqR6NwUjHjRvH+PHjkUqlRvqtn3zyidF+M2bMwMfHh6FDh97xtW6Hs2fPihkFX3zxhbi+Q4cOtGnThk2bNpGYmAjA2rVrGTJkCFqt9r5gpAZ3xssvv1znY1oZ6b3DSFu8Ib3dVHvJEhuuXKk++qnT6afUGs0/LgCFQoHBVgoCf+d2WmNpaSUa0aAgqIeKHomJiURFRfHAAw+QmZkpGlh3d3dxap2amiq2wVWpVDg5OREXF0dBQQHW1tY8//zz3Lx5k4EDB/LGG28gk8nq5CPV6XRotdoaI8CVgxqDBw9m3759gD4oZWjQ16VLF6Kjoxk9erRYcdUU6NmzJ4IgsHjxYtGQ29nZYWpqyocffoharaZnz56Ul5eTkZFBQEDAfVFbHhERwdmzZwGMUrxuh1YfaauPtNFQGyOtDP3nLQB6Y6jTaZFKZUilIJFIq+RrGlodC4JQRWGqPqhORq/68ekfCKVSyZUrV5DJZEgkEsrLyzl16hQHDx7E09OTF198kZ9//pmpU6felpGGh4dz5MgR4J+Id23o06cPLi4u5ObmcubMGZRKJQ8++CDh4eGAPtK/Z88epk6degefRN2wcuVKIzY8a9Yscbljx448+uij4ndx5swZoqOj72lGmp6ezsGDB4H6t4a51xhp5QKSxkArI21E3I4dvfOOfsqenp6OiYmJKOisVquxsbERSxLz8vKoqKhAoVBgbW0tquh7e3tjZXVnkevKMnoTJkwA9G1L0tPTcXd3Jz09Xewr5eXlxY0bN7C0tMTS0pLs7Gy6du2KnZ0dgYGB+Pr6AjBy5EjOnTvH1KlTa2Wky5YtM9qWmJiIUqm8bRWQn58ffn5+BAUFiZKEZmZm7N+/n6ysrCY3Wrc26zMY0eq6hbbkqH1NmrdarZbi4mLRZ294ftu2bSumddUVTcVIKyoqOHXqFOnp6djb2zNy5MgGn/PEiROMGzeO8vJyfvnlF4YPH97gc95LjLTFB5vqwkglEglyuZx27drRsWNHfH198fT0pLi4GK1WS3p6uqgIpVKpyM3NpU2bNnTq1OmOgzmVZfQWLlworh87dqwYuNm4cSOPPfYYAKNGjeLrr7/GwsKC2NhYbGxscHV1JSgoiIKCAjExPyIiwkj9qabKpuoMpiGwpdPp+OCDD0TGfeTIkSodOjMzM/nss89Yt26dGPAJDAy8bUPBhiIlJYWXXnqJhQsXiv5gV1dX2rRpU2Xf8+fPt8jKppSUFN555x2++eYb0tLSjDoX5OTksHLlSn777TdCQ0PZv38/Q4YMIT4+vt7XaarKnvnz5zNkyBCmTp3Kww8/TF5eXoPPmZ2djUKhYMKECYwfP56TJ0+KGSpyufyOXoitlU2NiLr66ywtLYmNjcXCwoJ27dohl8vRarVikKmyQXJ1da02yl8fVJbRMyhHffDBB7z22mtMmTKFdevW4ePjw7Zt21CpVPj4+ODr68vIkSMxNzfn//7v/wB9ldaSJUt4/PHHEQSBbt268eSTTwK155GOHj2aX375BYAHHniAadOmMW/ePM6fPy9WVUmlUr777jtu3rxZRYH91ii9jY0N//73vxv0mdwOGo0GjUZj1MZl2rRptG3bttr9WyojNbiJUlNTxaCYod+UQf/B4BONjo4mNjb2jlxHTcFIU1JSjFq6gF4j9ttvvzUS3L4VgiBw+fJl0f9vY2ODQqEgPz+f2NhYHBwcyM3Nxd/fHw8PD5YuXSq6NJYvX35H0/17iZG2eENaVx+pQTG/sLCQ69ev07FjR4KCgsjKykKtVmNiYkJ5eTkKhUIM+jQEAwYMqPFL/vPPP42Wk5KSMDU1rRIpN2DQoEEcOnSoyvrqfKRarZaMjAz8/f1xd3enV69ebNmyBa1Wy759+9i5cycAHh4eFBQUiAGs7t27G527cgXX1q1bmTRpUpMLQhtU3g2YP38+jo6OKJVKMjMz8fHxMdreUvNIKxdtyGQynJycyM3NJT093cjw/f777wwfPvyOfYaN7SNNTEyswvw7d+6MiYkJzz//PO3btxezR3755Re2bdtGaGgobm5ufPXVV7z99ts8+OCDpKamkpyczM6dO9m3bx/r1q2joqKC4GC9tkevXr04cOAAAFOnTmXx4sV3NN57yUfa4tWfTExMhClTpojTXYARI0aIPsXK0Gg0lJaWYmNjg4eHh8gCCgsLiY+PRyaT0b59+7uuplNWVkZcXBxOTk71+lFVp/5UWlrKJ598goODA6NHj+bXX39l4sSJIsN84YUX+Prrr+nRo4dYEjp9+nT8/f0BKC0txdzcnP/+979iWa1cLicqKoquXbs28p0b49bW2qAXNza4HW5NVFer1eh0uhal/pSVlcXevXuNqrNAbzwMLBT0zQkbmkrW2OpPhmfjVoSGhoqf/alTp+jTpw8hISHExsYycOBALC0t2bZtG3379iU+Pp7ExEQGDRpEZGQkOp2OcePGkZaWRseOHfH09AT0z+7JkyeJiIhg7ty5dzRFb1V/akRoNBo2b95MXFxcrfsJgkBOTg6Ojo4IgiDK6oF+GmJubo5MJjMSzrhbsLCwwM7Ors7s2oDqfKSGgoG8vDx++OEHSktLjZjcG2+8wZw5c0hISKBt27aMHTvWSITlk08+4aeffiIgIEBct3HjRjp27Hgnt1YvWFhYiFkGBly7dk2scCopKTGqxGqJPtK0tLQqRhT0U+b27dsDemPbGPm4jekjzM/PF1uZVAdDxsepU6cAfVCyW7duhIWFERsby+DBg+nfv7/oIsrPz+fpp5/mmWeeISAggJCQENGIgj5u0a9fPx544AGKioqqFfa5He4lH2mLN6QAb731ltgvvqbgiyGXUi6Xk52djVqtJjo6GpVKhUQiwc3NDZVKRXFxMampqfUWJ2koPD09RTGTuqK6yqbKjHbWrFkkJycbTf18fHz46quvKCwsRBAEzp8/L7ogdDodTk5OxMbGGrGnwYMH31aJqLFwayloz549WbBgAd7e3qxYsYIvv/yS0NBQwsLCWkxlU3l5OZs3b6asrExkx/7+/kYtkWNjY4mLi0MQhHrlitaGxqzsee+997h48aK4PHnyZECfclZUVMSGDRsA/QwhJSUFW1tbSktL6d69OxEREaKQuq2tLVZWVly6dAkPDw/RZ1oTYmNjWbt2Lc8//3y9x3wvVTa1eENqZmbGF198wf79+/nzzz9ZsWJFFWOq0WjExHbDX0dHRxQKhWi4HBwcxNLL9PR0oyqfuwGZTEZAQAAVFRV1DiDUFLVfsGABAN988w3e3t5GPxDQswFHR0cSEhK4efOmuD0hIcGo4RroRYYrM4m7gcrVPefOnSMyMpLk5GT69+/Pc889h5eXF0ePHiUqKqpZGalWqyUmJobly5cTFxfHxx9/zE8//QTop8m9e/cmPz+f/Pz8JunnP3GBAAAgAElEQVQn31iMbOfOnWJw0wBDTmtBQYGYR2xY9vb2ZuXKlchkMvr168cTTzxhlJr28MMP19mfPnjwYEJCQti+fTsDBgwwEu++He4lRtrifaTu7u6Cr68vhYWFHD16lH79+vHVV1/RqVMn0QeanZ2Ni4sLbm5uxMfHi4rj2dnZBAQEUFhYiL29vTgtNvRQagxcuXIFd3d3bG1tq/g/BUGgoqICQRBISEgQWXB1PaUMAbFbj69OIb+8vNzoAXvuuef49ttvjY49evQogwcPFpcfeOABPD09xUj/woULiYmJYe/evQ3Scb1T7Nu3r0ofqbfeesuoEqu5faS7du2qIoYtl8tZuHAh77zzTpOPq7F8pIGBgVy9erVWEfH27dtXSdF68803G2WmYsiaCQ8P59y5c+h0ujo9c/eSj7TFG1IrKyvh5ZdfZt26dUycOJGRI0dibW2NnZ2dmICv0Wjo1q0bUqmUmJgYzM3NUSgUYh6pSqVCp9PRuXPnBitEVYYhiCWRSEQWaGlpSXFxMTY2NuTk5Ih+UQsLC6ytramoqMDU1LSKHmVGRgZmZmbY2tqKD1lpaakYtb+11YihoMCQO1rd9xgSEmLENhYuXCgyk5bwvUulUnEcr732WpXvprn72icnJ7Nu3Tq+++47Zs6ciUajqZOOaGPhyy+/bJSofVpaGn/88Qd//vkn9vb2zJs3D39/f8rKyli8eLEYgPL39ychIQEXFxdmzJhRr/5at4MgCKxZs4bMzMw6P3urV69utqi9VquvljQxqZshbfHpT9bW1sjlch5//HFOnjzJb7/9xtdff41GoxHfsFKpFJVKhZmZGRqNRmSGpqam5OTkiA3n0tLScHJywtLSssGpPsnJyTz11FMkJyejUCiYOnUqU6dO5dq1a7zxxhtiddO3336Lubk5ZmZmLF26lLCwMMzNzfn000/p1q0bERERhIaGotVqEQSBxMREvvzyS0aNGlVrHqmdnR1arVbsgb5s2TKWLl1qtM+OHTtwcHAQl8+dO9ege25sjBo1Sqz9Ly0txczMjKioKHbt2sWECROaLY9Up9OhVCpFJmiYjt5NIwqNl0fq4eHBjBkzqpQQ7969m/Xr1zNo0CDMzMzo2bMnH330EZMnT25UIwr6YJxBRKiwsLBOrpDmyiP98Ud4+un6HdPifaQGRmdnZ8eoUaOYOXMmNjY2uLi44O7ujoeHBwqFQvQlarVa0ZCamJjg4OCAlZUVVlZWKJVKbty4cdsMgLpAo9Ewd+5cdu7cya5du9i0aRPp6ens2LGDkJAQ/vrrL4YNG8aaNWuwsLAgPDycGzducPz4cZYvX87rr78O6PsWHTx4UHyozczM6NKlCzqd7rY9m2QyGY888gigT2G51f9kb2/PmjVrxGWDKlVLgYFNGyAIgqiXamFhcdej9oIgEBoayjvvvMPy5cvZu3cvAC+99NJdG0NlNJWPsKKighMnTjBt2jRUKhVHjx7l5MmTmJmZ4eXl1eBilepQOX3R4F66HZrLR3onE4B7gpFWBwOjNLAHW1tbcZpvmBpLJBJxumhiYiLuk5WVRVFR0R3L0xUVFZGfn0+PHj3EFiPt27cnIyODAwcOsH37dkAfGZ00aRKvvvoqBw4cYNKkSUgkEnr16kVhYSGZmZm4urpSXFxMcXExe/bsYdiwYVhbW1NSUoK1tfVt38iVH1AXFxeKiorQ6XTiG3/27NmMHj2an376iVdeeeWO7rep0K1bN1avXs3cuXNFWb3g4GASEhKws7PD19f3rjLSypkc48aNIzw8nN27d9/VLgKV0ZiVTWq1mqioKE6fPs38+fPF9U5OTlhbW4sBxzuJrtcFr732mvj/hx9+uE7HNAcj3bgRYmNhzBjIyZnK39lgt0WLN6S15V4KgsDRo0cpKiri0qVLaLVasXqpNuh0OiIiIlAoFMjl8ipBIjc3t1q/bIP4iWHql5yczOXLl+nRowc5OTlimaChxbChxUTlVBF3d3cyMjJEQwpw8OBB/vWvfyGVSikrK0Mmk91Wj9TU1JS3335b9KMZ9l22bBlLlixBIpHg7e3NokWLePTRRwkICKjSzbM5YYiCG+Dl5UVkZCTh4eHk5eVhbW0tlsw2NQwv57lz57Jq1aq77hO9FY1Z2XTs2LFqfc3du3dv8t5d0dHRohunY8eOt02ZMuBuVzZFR4PB+zFiBMyZswEzs611OrbFG9Laau2LiopQq9Wi4azsH60NUqkUExMT1Go1SqUShUJRrwjszZs3RR9SaWkps2bNYtmyZeJYDcpAldlURUUFRUVFYhQfEJmztbU1N27cICEhgcDAQLHiSCKR1CliW52/d+nSpSxZssRoXceOHVtEkKkyZsyYwdGjR5kwYQIlJSViiavB2Fdm3E0NExMTFi1axLp16wgODr6tLGFTozEZ6dChQxk7diylpaW8+OKLXLp0iWXLlvHnn3/Sp0+fO25lXhcYUq3mzJnDV199VefjZs6cSUmJQGIi+Pk1ydCM8N//6v8mJoK+cLLuNqHFG9KaGKnhAXvooYfEacnly5exsrKq11RMp9ORm5uLo6Njndswe3p6kpeXh1qtZtasWYwfP57Ro0cD+qnS5cuX8fX1paSkRMxddXFxISUlRVTaSU1NFZmrpaUlhw8fZvDgwZiZmaFWq9FoNGi1WkpKSurkgggNDeWLL77A29ub8+fP06FDB+zs7MjKymq2qWld4Obmhr+/P126dKkisQf6Cpq7CSsrK/r378+RI0ea3ZA2JiOVSCTs3LlTfHlPnDiRKVOm0KVLFz788ENA346msY1p5RfBf//731rTnsrLy/njjwiuXk0gJ2c4P/+8j+TkTOAdwsOhUjZfk8AQYsjIMBjSuqPFB5tqYqTFxcWYm5sbTRMUCkW9fWpSqRR7e3uysrLqzNZsbW0pKytj/vz5tG3bljlz5gB6Rjx8+HD27NlDaWkpGzduFCtzBg0axN69e8WAiq2trWhIpVIpf/zxB1OnTsXe3l4MpBmCZHXF/PnzGTduHH369CEuLo7CwkKR4bVUtG/fnry8PGJjY+natatRtRDoex1lZGTclbFoNBp0Oh2xsbFs3br1jsoaGxONXdlzqxHr3LkzL774orgcFxcnVjA1FgwVdLt37642H1Sng/XrlSxYcIT33vuMqKgIsrOLsbTcxMCB/ixc+Ap+fjBnDvxd8d1k+PRT/d9b6lvqhBZvSGtipCqVCnd3d6M3qIWFRa1JxzXBoFVa17JRc3NzMjIy+O233zh69CjDhg1j6NCh/Prrr0yaNImzZ8/yxBNPEB0dzRtvvIGDgwPDhw/H09OT8ePH8/777/P222+L50tOTiY9PZ0HH3xQXCeRSG4bta8Jbm5u4v8NUf2WCj8/P8rLyzl8+DDz5s1j7Nix4jbDS/RW49oUyM7OZu3atXzxxRfEx8dTUVHR7HX+dyNqXVmR7Oeff2604N7OnTsJDQ1l3759/PDDD9U+h2q1ms2bTxAT8zkODofRaBwwMQnE29seHx8f8vNXERf3FJ99VkhsLPxNnJsMfn7g5gZ/x4rrhRafkK9QKIQ5c+aQm5tL27Ztyc3N5cknn8TJyclIeAP0Qh5paWlGuZO3wpCcX3naLQgCubm5ODg44OzsjFKprNN0WqfTkZaWhkKh4ObNm3h4eODg4FBj0r+hUiklJYWsrKxa9628/62VTbeDQSEK9A9rcwZM6oIBAwYQERFR7bY5c+bUOTjREKSlpbF582aCg4NJTU3lu+++a/Y6/8ZWf7oVWVlZuLi4iEz17bffbnB+9cWLF7l58yaRkZEArFq1irlz5xrto9VqiYiI4PDhwwiCgFYrRSYzNuCG8tRjx45hbm7OtWuPsWFDRy5cgFsaKTQqFi7UM9PycjAzu48qmyQSidEAFQoF+/fvp1+/fkZGqLCwkOTkZGQyWY3JvjqdjqysLGxsbMT0IoVCIfo7K8Pf379Wg9xQGB40S0tLrKysqn2Aa6tsuh2ysrL48ssvgX8i+C0Z/v7+ol7poEGD6N27N0ePHkUulzdKK4y6wJCClpuby/Dhw9m8efNduW5NaKzKpupQUFCAvb09kydPZtu2bTz11FO0adOG/Pz8OxZduXTpEn/99RcbN25k4MCBNfrmN2zYQFJS0t8Exo7ERHvS0j4nLS2NrKwskRWPGzeOtm3b0qFDB9LT04mK6kN+/nAOH5Yjkehbd7/88susX7++0V62w4ZBZiacO6f6u1HmfVLZBHqBiDVr1jBjxgxGjhyJm5ubkRFNSkoiPz8fGxub2zI8qVSKm5sbWVlZotiJlZW+g6ipqSkVFRXk5+eTkJBQ78BVfRAYGMjFixcpLS1FrVZjZWVVZex16SJaExwdHQkODiYyMrJJBDUaGydOnBB/DIZqtBEjRty1PNKCggLi4+N57LHH+Pjjj2uVnLtbaMouogY31pUrVxg1ahRFRUV89dVX5Obm4uHhwdChQ0VZwJogCAIqlYoLFy6IFWqPPnoJX9+u1PSzUat1FBcXI5fLycgoZdWqlzh2DAYMeAZBEPjxxx959913iYuLo0OHDixZsgRTU1N+//0QcIrMzCS6dl3DoUP7efPNNzlw4ABvvvkm3333XYM/k337UggP92TKlHg6dBjB+vXr63xsi2ekNjY2QlFRkdG6mJgYAgICkEgkFBYWkpCQgLOzszgVNjSBu3VKa2j21qtXLyQSCaWlpcTFxWFnZ2d4+wD6oINBF7Nnz55NlhaSmZlJcnKyuGxnZ2dUmmdgpNbW1iQmJvL444/TvXt3CgoKiI2NpWPHjgwbNqxaNnv16lW2btXnwPXo0YOzZ882izhJXVFZxGT+/PlYWlpy6dKlu1Zrv2fPHiIjI/H29sbPz4/r1683i3ZtZTQlIy0pKSEkJAQzMzMcHR1xcHDA1taWq1evij28apvq79y5k/PnzyOXy0UlNVvbeAoL2wGweTN06ACdO0NkZBGff36Jgwe1eHtbMmnSHgCSk2HDhiWUlEgwM9MbZkdHR4qLS+jbN4bjx3+nd+9M/vjjHf79bzh9+irjxm1GKtWXh69YsQIvLy9SU1M5e/YsPXr0uOPPIyUlhbZtY9Bo+qDT+dG5swffffcdffv2vT8Yabt27aqsMzMzE1OWlEolMpmMnJwcMX9TJpMhCILoW9JqteTm5iKTyQgMDBQNilarRaFQUFZWRlFREdbW1piZmVFaWoq1tTV+fn5NZkRB39+nrKxMbHxXUFCARqMR/bMWFhbodDpKS0spLi4W/bqGt/+JEyfo0qVLtTJ4lfMvo6KiRIWslorKzfkMVU5vvvnmXWGkOp0OX19fUc4vOTnZKJrdXGhKRmplZVVjEG/Pnj08+uij/PHHHwwZMsRIqUypVBIXF8f58+cBqKhQ8t57UkJDobLUrqGGQiLRYW+vpqysFxpNOjk51wGD/9+a9u31RlS/r4QPP/yQF164zvHjbYGZnDkjYG+v3x4cHM9XX33GxIkT8fX1ZdKkSXzwwQf069ePhQsXEhYWdkdkobi4mJ4956JS7UIme4P333+FRYsW1Ws22uIZqZubm3Br+otarSYlJUXU9iwsLDSqNCooKDBqTVxRUYFEIsHJyanGN2xZWRkFBQVotVpkMhlubm7VGtGysjIxF9Ta2hobG5sGMT2D31an04mBIUOaiEqlQqvVolarOXfuHDdu3GDr1q2MHTuWH3/8UTzHW2+9JQpYV37oK7fuSExMrLY9S0tC5c9x4sSJ5Ofno1arm5yRnjhxoor/OTMzs9lfPE3JSG+H7Oxspk+fzvHjx/H19RVblxcXF/Pnn3/i4+PDyZMn0Wo9qNxqy9NTha3tDYqKLOna9TJdupzF2lrPWE1NTcnIKMbeXm+gtFoJsbGv8NNP/zyz8fECHTpU4Op6nbZtD3PiRBZgqGyS8uqri8WmkoMHD0atVuPi4sJLL73Ezp07xa69tcEQ8D19+jSnTp1i7dq1FBaOBH5iz56bjBnzzw3dN8Gmnj17CvVVLVq+fDm7d+9mxIgRlJSUsGbNGk6fPl2lb/qt0Ol07Nmzh759+9b4I3r//fdZuXIlI0eOZP/+/ezYsUPsIHknqKioEI3fwIEDCQkJEQ1KdXqcWVlZ7Nixgx07dohRZUOV1sCBAxkyZIj4sqhsSCMjI+nVq9cdj/Nu4PLly3Tr1s1o3euvv96kup+CIBATE8PPP/8srrt27VqNnU3vJpo6an87ZGZmsnTpUqKiokhNTSU/Px+lUmnEkvv168fFiy9RUXEOJ6cf6No1gC5dumD/N41s164dAQEBKJVKCgsLOX36NBqNFLlcP9MQhLGEhvb4+//w0ENw7FgpanV7/PxMCQ4+yfbtboAtUERmZianTp1i7Nix+Pn5MXHiRCwtLTl79ixxcXFcvny5Ria5du1a9u7dy+nTp8XcZEMZ9kMPvc2hQ++wbx+MGvXPMfdNsCktLa3ex9ja2pKXl8f58+c5ePAgCxcuvK0RBX1ifOU8xurw2muvER4ezunTpzExMeHYsWMNMqRmZmZkZ2fj7OxMQUEB165do7CwkODgYC5cuEBBQQEPPfSQuL+LiwtmZma8/vrrREREIJfLeeCBBwB9PfWxY8dEA2ptbS3m4d6aKtYSUV3zvb179zJhwoQmu2ZERITYwdXe3p7c3NwW40tu7C6i9cGuXbuYNWsWRUVF9OjRg6CgIBwdHXF0dMTJyYmEhATS0tJwdnamX7/TWFhYIJXqBU/c3d3p06cPAQEBVQKoFy5coLxcQC5XUVhog0JxHkHogU4HL74IYWHwxRcKUlOf4aOPPmLIkDepqPiW5cv/w3vvfYKr66OAvnw4MTGRNWvWMG7cOHr37o2rqyv/93//x/jx4/H09DR6AQmCwH/+8x/kcjljx46lT58+lJaW8uqrrzJ9+nTefXcZvr56v+2doMUbUg8Pj3ofM2nSJLZs2UJhYSGHDx8W28TeCq1WS1hYGN7e3nU2NDKZjM2bN7No0SL69u3LU089Ve/x3QonJye++eYbZs2aJcrIBQcH071792p9hI8//jgff/wx/fv359133xXzAQ0IDQ0VNSWLi4v55ptvmlVlvD6o3NHS3Nycixcv0qFDhybrcHrz5k3x//n5+Rw7dqxKX6nmQlP6SGuCRqPh008/ZfHixfTo0YP33ntP1LTQarVoNBqUSiU+Pj5iW3MHBwe8vLzw9fUlICCgVi1TqVSKRCIDVKSnm9ChQzIDB8bx7LMKvvnGD4CNGxPo1i0QkNC9ezc2bvyGtLQnUasNKXwvYWHxDWVlZZSWlrJ582b69OlDcHAw5eXlYpDVxsYGLy8vHntsEh9/LEEiiUYuX8qGDd9x8uRJpkyZwqBBg/jss8+wspIgk+lISTnPsWOlDBgwoF4v1BY/ta/OR9oQHDhwgKeeeor27dsTExODhYUFRUVFvPvuu0a9hJoDlduDtG/fHj8/Pw4ePFhtDXRaWhpr164FYMWKFXTv3p01a9awY8eOKud99NFH2b17d9PfQCPh6tWrRjOIl156qUk0Mg8fPmyk0Tp16lS+//77BielNxbulo+0sLCQgwcPcvLkSfbs2UNcXBxTp05l6NChpKWlYW9vj6WlpdjCx8TEBBcXF7y8vPD09KyXCPS6detITc1AEDTodBJAQKMxIS7Ok/j4eKKjPdBopgL2wDPApr+PHA8YdEzdgX9sgkQiwd7enrfeeovXXnuNefPmMXnyZFJTU8nJUbJ27ZOcOwc2NmUUFVlgYbGOsrK5DB3aj8GDB1NUVIS3tzcJCWY4OurFp/v2DaFbt/Z4err/7/pIa4MhudvQIM/Dw4M9e/bw5JNPNrshBf0UZM+ePWzfvp1NmzaJ6ydNmlSFlSmVSlFwIioqiqCgIARBYOXKlSxcuFDcLzw8vEHuh+ZAZTZwa7/7hiIhIYGLFy+KkWfQK2NdvHixRQm83A0f6SuvvMJnn32GRCLBw8OD3r17M2DAAPLz87G0tGT48OFGmS4NwV9//YVSqeTo0aPodDoEASQSKCy0wsREjYWFEpVKRmysBadPP0RysjnQC5gBLAJ8gMWAvmrP39+fdu3aERcXR3p6OiqVCicnJ/Ly8oiMjCQoKAhXV32N/vr1MHKklocfjuL48WAghs6dV9Knj4SSEjmdOztRVubM/v2D6dgxjsDAS2i1Et59d+n9YUgbk5FGREQwcuRI5s2bJwZ4rl27xh9//MGJEyfqHWBISkoiLCyM4cOHi2/txoIgCCxYsIC//vqLM2fOADWr8xgMzaFDh8QI9+rVq5k3bx49evRocS1G6oINGzYwc+ZMHnjgAUZV9v43Aq5cuSIGl4KDg9m/f79YLtyS0JSMtLi4mKSkJA4cOEBRUZHRcyWXywkODhbVyBoLXbp0EfO2u3XrRq9evbGy0rNZQQCtVkpxsRVmZkrkci1ffCGhqKgCeAcTk5Oo1ct44gkHtm7dAuilASv3JKsMiUTyt3vrR374Qc6FC9Chg4atW7cybdoPwJeAP5ALOGJlFc2ePV2IjASlch9qtf43Fxoaen8Y0sZipFqtFj8/PwYOHEjHjh0BfRLu5s2bOXDgQI111efPnyc/P5+hQ4card+6dStPP/20mDAvCAI///wzw4YNq5bVCILAt99+y6lTp/jmm2/q9IY3MJKgoCCuX79O586dcXZ2rjIWQRC4evUqP/30k6gZoFKp+OCDD5g6dap4v/cSNmzYwJIlSxg5ciQWFhaNWq4rCALLli0Tlw0C2y0Njc1Iy8vLCQ8PJyEhQcxdNjU1xcfHBy8vL5ydnXF2dsbBwaFJ8qfPnLnChg2ncHcvISSkFzExMZw7dw5BEHBwcMbERIrBHGk0MoqLJaxbF0Rp6UeADlPTb1AqnwEib3stc3NzysvLMTX1RCKJw8MjC0EYxo0bCQDI5W3QaM4BdsABQD9TlUp1vPnmCqTSciSSuhvSFq/+dCdR++oglUrFHM3y8nIyMzP57bff2Lp1a41GdNu2bQwcOJDHHnvMqATt7bffZurUqTz33HOEhobyn//8hxEjRjB79mxsbW1p3749EomEJUuWkJWVRXh4OE8++SRvvfUWBw4cwNnZuU6KRt9//z3vv/8+165dY/Xq1Vy5coUjR46wZ88eSkpKxECUoYMp6AVABEFAoVAQGhp6TxpRQIzUf/vtt3z++eckJiY22rklEgmPPvqouFzZqLYkNKb6U1FREevWrePs2bOiGtmsWbNYvHgxTz75JIMGDaJTp044OTk1WRFKUFBnLl9+EKWymL/++ouHH36Yt956i0GDBmFhYYogCEgk+um+iYmW+PgTdOq0CpnsU+BnlMrOODqurVa39tdff2XLli1ivnhSUhKLFi1Co8mgouI/JCT4odXqn6nOnXuj022iT59ohg/fg97nCiDQu/dpJJJyQM+Q64r/GUYKsGjRIlauXImlpSUKhYJVq1bx+OOPV7uvWq3G2dmZ8ePHI5PJ+PXXX8W3eJ8+fYiOjmbx4sVVjjO08AV9Oo1BmLhv374MGzYMuVzO+fPniYuLY+XKlQwdOrTG3uG3MpKCggKmTp3K77//Lu7TtWtXRowYgZWVFdu3b+fKlSu8/PLLYtvlexm3VmOZmpqKTQMbA8nJyWzfvp2ffvrprgmj1AeNxUhLSkrYtm0b6enpPPXUU81amLFxI7z6aiZz525BIilBpxOQy82oqGiPWm1Djx75nDlzGisrC1EbNT29LevXzwQkwDxgtXi+kJAQtmzZwi+//MIrr7yCWq1m06ZNTJkyBdDnXc+a9QK7d7+FtXVbios74uT0C6NHx+PvnwjAt98+S3q6PWPGHKJnzwt/uxkkyOVCnRlpi09/aixGCvoa7qCgIAYOHIhOp6NNmzY17muogffz8yM1NRWfSuUboaGhLFiwoNrjXF1defrpp8XS1oSEBPbv38+IESPEN31QUBB5eXmMHj2a/fv3M3z48GrP9f3333Pjxg0xoGRnZ8f+/fvFDqPbt29n5syZXL58mVGjRtGvXz+uXLlyV2Tn7gZ27NhBp06diImJAfTBtbS0tDtKiasO+fn5olZDS0RD80hVKpWY8ywIAuPHj2/26rbp0+H771356KNZTJmyi4ICCw4cGIlabU6bNnHcvOnO008P4uTJVZw/f/5vTQAFzz77HevXP0vv3oPo21fGoUOHuHbtGmFhYQwcOJD4+HhCQkJYu3Yt/v7+4vVcXFxYuvQtdu+ei5/fZjp2fB8vr9NYWeUTHn6WoUN74+sbw+DBSbRvnyYGwOTy+hHM/ylGWh9EREQwYcIEXnzxRQ4cOED//v1ZsWIFoO8PP3r06AbVYxcUFLBy5UpA3552/PjxVfapCyMRBIEzZ85w/PhxPvroIyoqKvjkk09E1f57GeXl5ezcudOo+d2CBQvEqpmGQKPR8Msvv3DlyhVR9auloSGMNCEhgd9++42CggJ69OhBv379WkxAraREr/kZFwe9eoG3N4SECKxevYbCwjK0WgnW1iVoNHqFKkPMobjYnM8+06DRvEl0dDRubm5s27aN7du3ExcXx/Hjx8XcVgO0Wi1jxowhLCyMN998U1y/adMmEhKu8dxz3+HunopMZuijZjzWVkbaQAQFBeHl5cWRI0c4f/682MAL9Ar0RUVFRvX89YVBxkwulzNhwgT27t3LqFGjjIJQtzLS6iCRSOjTpw99+vRh2LBhREVFMXz4cIqLi8nNzaWsrIy2bds2aZllU2Hjxo1G6ljAHbfQvhXHjx/nypUrgF58uCWiIYw0MjISmUzGjBkzmp2F3gorK6jUIALQi59ALlZWWvRTeKESIw0BwNq6nDfegBUrBtClSxemT59Oly5d6NixI++++24VIwr6z/DAgQOEhn4NpIvrXVxeJCFBwsGDDzFz5oYqBrS+aGWkteDSpUsEBgYyY8YMI21CQRB47rnnWL9+PTNnzqvkGQEAACAASURBVGzQg2pIvQJ958/KOZMNYSROTk6iT7dTp06i0biXYLj/oqKiKgpXS5cuvePcxoKCAtasWSOy0AsXLhAYGNjg8TY2GvL9l5eXY2Ji0uK7I1RGbm4uly9fxtLSkosXL3LtWiZhYT3p2/codnYK0dhpNBJWr64gP/97IEU8/siRI1Wq0vRpgP/hjTeWoFD8I96+f/8ITp16kLlzv8DZOa/GMbVG7RsB3bp14+DBg3z++edG6w05aqAX2mgIAgMDxTfprVNWQ9T+TrBs2TIxpedud+JsLGzcuJEPP/yw2ilpQ0onDfKEBrREIwoNi9qbm5u3WCOam5tbbS8yR0dHBg8eTHBwMM8++yxHj+Zz6tQxbG1NjRijXC7w4osWmJnF0r//JkaP1mdgVJ66G9C5c2dgFHFx/8hxqtUybGxKcHZOx8kpt1HuqcUb0sYKLNwpHnrooWo7mRqm5oY2zHcKNzc3UZXpVim3adOmVftw1AVz584VXwC//vprg8bYXHjmmWd4/fXXUSgUbNy40Wjbe++9R2lpaZVjjhw5QmhoKFevXq3R2H7wwQdiRLiyIExLQ2N3EW0J2LNnD25ubtjY2NCxY0fmzp1bbZPBwsJCLl8+i6fn00ilAoKgzz7IzNTnmioUWmbNWoerqyX79+9k8eL3iYiIYMsWJf/+N+zerZe83LZtO337zqFtW31AUaORodHICA6OZO7ctTXOauo7UW/xU/vGrrVvLBw9epQJEyYwf/78Bp9Lq9WKfrDAwEAmT55MUFAQ27Ztw83Njeeff5527drd0VS2pKSk2WTYGoo1a9aQnJwssnKtVmvEsnx8fHj22WfF5evXr/P9998bnaOyq0QQBD7//HMjhn7mzJkaRW2aG82pR9rY0Gg0HDp0iCNHjqBQKJBKpaSkpPDbb7+RmZnLmjXX0WiusX79e0ilGqCCM2fOYGX1LIsW/ZMxU1CgwtTUCnNz/YtQpZLz8cevEBCQz+XL45FIIpDLJWg0AgMGLCI4WI6tbd21AG5Fa2VTE0Kj0dC3b1+8vLwa1N6gMlQqFVevXkWpVLJ3715x/aJFi1ixYgWHDx8WBU3+V1Cdj1ClUnHw4EGxva+trS0mJibk5ORUOd7Z2dmog6UgCPz6669crNS4vKKiosUG4ppbj7Qx8csvv3Dp0iXUarVR3rROJyE2tgNarRQLizIuXfIlICCdCxcCuXLFhaCgq4wbFy6mJVWHvDwbVq16CZ1Oi4mJihkzbnDxYjEPP3wAlUrG9evt6Nw59o7G3Rq1b0J8/vnnxMXFcfbsWRQKBV26dGnwORUKheirM/xdtWqVmHJ16NCh/zlDaojaV/YTKxQKsbcT6KeANSE7O9toWSKREB8fLy537ty5xRpRaF49UgMEQSA9Pb1GF1tiInzyCfTrV063bvF07dqVnJwcfv31V7Htj6urq5gLXNmIfv/9k7Rpk0SvXmcxN9cH/vz8kpBKoWPHeH77LRHIBB6oNapua1vMuHG/sH//A0gkmykt9WDUKH1LE1NTLZ063ZkRrU+XmzobUolE8gSwFL0ESwYwQxCEYxKJZBj6UgMf4NTf65P+PqYd8BPgAPxLEIT9f6//f/bOO66q+v/jz8O97L0UkBC3koaoKCq50NzZMlNLXC01v65sJ/3KLM1yZFbmTs3cmqsstygOQFAUARVQRGSPC3ed3x+Xe+QyLwhG2evx8HG9Z3zOOZdz3uf1eY/X+wgQALQQRTGpeFlf4CdRFL1LHvfv9pGWBwcHB0kwuS4k1/QP9/jx4/n999/JyMiQGNijhODg4Ar9nKdOneLMmTMEBgZKU/MdO3aUEYHesmULhYWFDBw4kNTUVBQKhbSuMiNcH/Cw9UjDwnQqSc2bg4dHNhcv/srx47+QlZXGTz+tICCgC4WFkJcncupUOB9/nE1UVCBarRyFYgdxcde4cuWKJFZua2tLVlaWQQpbYSHFje4gISGV+PguZGQ48PTTullYyerU/v09gKqff5lM5IknonniiWjUagG5PN5gvbEescTERFxcXCRZwOpUyhplSAVB6Ad8CYwAwiguThUEwQWdSOBEYA/wKTrDGVC86yfAVOAasBvYX2LYfOAj4LXKjl0fGen48eNp1KgRAwYMqNMk5/j4eDIzMwkKCpJU8B8llMdI9ejatStdu3Y1WPbss88ye/Zs5s+fLy27dEmnpl5ermjJarX6iIfNSM+fh++/BxMTNT16nMfHR0lQUA9MTODAgQOsXHmRs2f9yc+3oaCgIbdvtwfU+Pltp3Hja1hauhETE4MgCIiiSFZWVpljpKQ0okmTW0REtEcUdRoXvXodLfd8SueRGgN9Yn11cfr0aQ4cOICX1wDGjw+oeodSMJaRfgL8nyiKp4u/3wIQBOE14JIoiluKv4cA9wRBaC2K4hVAVnwMGVC6oHwJMEsQhPmiKMZRAeojIwWk1KLaShAvD76+vly5cuWBhEfGjRvH0KFDefbZZ+tNCw1jURkjrQjvv/8+BQUFZQxnQEAAp0+flr7b2dlx/PjxWjnPusLDZqRPP53FF19k8OKLm7GxURqsEwTw9LyNp+cuadnRo0+SnOzNkCG6HGWFQhcU1rdFt7CwYMCAAZw8eVLyYTdpomtxXVhoCoi8+OJm7OzKpkKBriimuqjJLa5QmKFv+Z6YeID8/HbVlsSskrwKur4AnQBXQRDiBEFIFgThW0EQLIHHgUj9tqIo5gPxxctB1/5vCbop//8ZjswtYAUQUtnx6yMjBV27Y29v7zKVN7WJyMhIlEol33zzTY32379/P2vWrOH555+v07bSdQV9Hml1YG9vz9KlSxFF0UDwpKQRBXjqqafqjRJ+RahN9aeqsGPHDgYO7MG4ceuwtlYWC3dAUZGcmJgWREQ8QV7e/W6fogg9ex5n0KDvkcm0FBSUvb8KCwvZvXt3uYHArl3PMmfO/1UaBIqIiCg3NepBUF5s3cJCKRXFAFhZWaPRVM8iG/N0NUTHJl8AngTaA37Ah4ANUNrRlA3Y6k5avCyKoq8oil6iKO6iLOYBQwVBqDBaY2NjQ1FREfv3769Xn0ePHkUQBKmM89q1a6jV6gf+jI6OZtmyZaxZs4bMzEwSExNRq9U1Os/SOa5paWl/++9Wnc+XXnqJjh071nj/u3fvGlx/yRJCW1vbv/36qvocNWoUfn5+tTruvn37uH79Ovv27ZOW5+Tk8Oqrr9K8eXPpPoyLu0ZenimJiTG0aBFDREQKX331BuvWOaFS6dar1Wru3buNUqnm1q2rBvdxXFwcGo2G2NhY1Go1sbFl73eVSlPucv1nu3btcHd3r5XnquR1lV6u0agNnhWN5v52xsIYQ6r3zi8VRTFFFMV7wNfAICAPKD23tQNyjTm4KIppwLeUZasSbt26RVFRERs2bKh3n/ocUhMTEy5evIharX7gT30K1I0bNzh58iSgywWt7vnpmUxJ1rVixYp68bsZ+7l69Wo+/vjjGu9/5MgRqVps9OjRREdHM2rUKC5cuEBBQcHffn1Vfa5cuZI5c+bU2niRkZF8+umn/Pjjj3zwwQdERUWxYcMG1q1bh7W1Na1atSIqKgqNRnc/mpvncvHiRZRKLba23zNo0GZSUg6j0ai5cOEySqWWixcvotVWfV9HRZVdrtWqyl2u/7xw4QKHDx+uleeqqk/9zNLU1KqY0ESiVBrvVjEqj1QQhCTgA1EU1xV/fx4dI10OBIui2L14uTWQBnQo9pFWNN4R4GdRFH8SBMEBSEAXzHqzdNS+PuaRgs7l8PLLL3P48GGmTJlSJ0EnlUpFbGwsW7ZsQavVGu3jzM/Pl1T7X331VUmdH3T+q38K/k15lDVBbV5/VlZWGUFkmUzGiBEjGDhwIKNGjUImk5GW5oSLS4bka9TfLoKgbwciA0QKCy2wsSmQtinv1iy5b0XbVZYfqq8+q80+WpUdLyIiEicnRykIGR3tw9atL9ZqHulq4C1BEA4AKmAa8BuwA1hQbFj3Ah8DFyszoqUhimKWIAgL0XW1KsNk66uP9MCBAxw+fJguXbrUWeQ+MjKS2NhYQPcgGCsfZ2pqiqmpKZMnT8be3t5gXVxcnKSVWt9RWdT+UUBNo/Y5OTlkZmaSk5NDdnY26enpREZGltnOzMyMjRs3YmPTAJlMhkZjgqtrhoEfURBArTZBLtclVcpkGoqKTLGxKZCMUkWGqfTy8rarjBvUJGpfFSoz5O3b+xqsc3Aom3VQ4bhGMlJTYDEwCigEfgVmi6JYWJz/+S3QmPt5pDeqGO8IxYy0+LsNOlZa8E9hpHK5HI1GU2FDutqASqXi+PHjeHh48PPPPxsdHFEoFFIu3Pjx4/Hy8uLPP//k+PHjeHp6kpiY+I+I4P/HSKt3/aIocujQIYPASWmYmJgYCLbo9rvPOOv6tqjOMeqCkVYXtar+JIqiShTFSaIoOoii6CaK4lRRFAuL1x0SRbG1KIqWoij2qsqIFu/TS29Ei7/niaLYoLQRhfrJSNPS0qS0lLpMT4mMjCQ+Pp5ffvkFuVxe5gGoCOUp7uurr5KTkxkwYECtnmddoSZR+38TqhO112q17N69m1OnTtG+fXtefPFlGjcehP4Rd3R8ggYNhldoRKFujagoVt9Q10XU3liEhnbm0CHjmfB/tfY1gL+/P+fOnWPGjBl1mkeqUqlITExk/fr1CIKARqMxiknqt3nhhRdo27YtgNRVVI/k5OQyGp/1Df8xUuOvf9u2bURHR9OtWzfy8xWsWGGNv7/OCKWkNOTy5Ta4uNyjXbto6nnWl4S6YqSiCPn5FtjYlN8V4fffPTl16k1ACQj/1drXBWJiYjh37hxvvvlmnRpR0DHS3377DbivXGRnZ8eYMWPKnearVCp27twpfS+Zv1f6d/T09Kz3gaf/fKQV+0jVajUymQxBEFCpVERHR+Pv74+/vz+LFy+mc+f723p4pOLhkfrA51MdRimKOsm6Q4eacubMSOztcxg16mcaNiybU1oRauojrcpVIYpgbV1U4f4ZGU2B6skr1ntDWt8qm3x8fHB3d6/zPugKhUIyoqAT6t24cSO3b9/m6NGjrFmzxmB7rVZLy5YtpbbFHTt2NKgM8fb2ZsqUKdja2krT5fj4eJo1a1an1/EgqEll078JFVU23b59m59//hk7Ozt69eoldWgoLCwst1WxHg/qA63Ovj/8MI7U1EY0bnwHEMjOtsfSUlHlfiVRk8omMDzP+HhvPD2TMTe/nxN68mRXunYNQy4v+9uKInTrlsDzz3fD1PQTSqgwVop6b0jrGyMFSElJQavV1mm1kJmZGQ4ODmRlZWFubs7s2bMRBIHCwkIWLlzI8OHDDVSQTpw4IQmpgE5wujRrdXFxobCwkN69e3P48GHOnTtXrw3pf4y0LCO9ffs269evx9zcHIVCwebNmyVXTlRUVKXj1WVA6fbthri63sPUVGec3nhjNVotZGdbk5npglyuws6urBB3ZXjQqP3x47mcPfsWHTp0xNu7MTdu3CA+Pp7k5BDu3fMtt+GkIICXV3LZwapAvTek9Y2Rvv/++3z++edcvHixxm9MYyCTyZg8eTIRERG0adNGelgsLCxQqVQMGTKEZ555hlWrVuHk5ARAUFAQf/75p7R/eVi6dKmkLP/SSy/h7e1dbwRR4uPjOXz4MD4+PnTr1u0/RlqKkaakpLB+/XrkcjkODg7cvHkTuJ8bHBfXg61bA5g5cz6mpZUt6hju7qllDLSJCTg65uPoaLwBNUxHMv75ysnJYefOnXh6ekqG96+/FgJw9OgRjh4Fb29fWrduj5mZVblVSzrfqSWWlgpkMjhxopvRx6/3hrS+MdK5c+cybNgwBg0axI0bN+jbt2+dBUMiIyPJzs6ucPydO3dKPtEBAwbQsWNH/vzzz0p1S0u35yivjcrDxLVr1/juu+/Yt2+flDML0L17d5544gkcHBwMgmSPEkoz0r/++ovCwkI8PDwkIwogk9mi0eSSnZ3DmDGfYGpqWdGQdcJGa9PVXvL8jGWkkZGRUjudhIQE+vTpY/D7jBw5kkaNGmFjY4NCoSAmJhxf385l2LkggLX1ffdDRITxhrzeK1nUN0YK0LlzZ+Li4vDw8GDv3r11xpp8fX158sknDZbpfaB6NGjQgA8//JCAgABMTU0ZO3Ysvr6GicUlMW7cOF555RXpe22IUj8I5s+fz6JFi2jatClLlizh0qVLLFmyhOTkZJYvX86OHTtYt24dKpWq6sH+JcjKyuLAgQOkp6djZWXFrl06mQp/f39Apxmgz8YA0Gh0Lp2OHSPw8KjYiNYWShvOypLyH2Ts9u3bExgYWOl2aWlp5fYk0yt7DRs2jFatWmFjY0NKSgobN26kUaNGvPnmZDQamcHxSl/XxIkrjD7vem9I6xsj1cPBwYENGzYQExPDp59+aiAYXFuIjIw0kHrLzc2Vgkz6G2zMmDEGfYy8vb2lqX55cHNzK9PXKDAw8G8TOW7atCkAmzdv5q233sLHx4e33nqLa9euMW7cONLT0wkODqZZs2YsXrz4gQ2qKOr6pddHw6x/0BcvXkxYWBiJiYn8/vvvRERE8Mcf+TRv3hIfHx9Onz7NsmXLOHq0rI6nsbnGD4LyqoNqi5WWZqQV5ZHeLzs1PHDz5kdJS8skLk6nzOnn54dGoyE0NJQ1a9YwfPhwXnzxRVq27ENBgaXkN9aPKQhQVKQTVrewMP4eqfdT+/rISPVwcHAgJSWFFi1akJeXh6Vl7bIBX19fgwfDxsaGDz/8UDKcQUFB1a5QMjU1ZcSIETRr1ow1a9Zw+/ZtTp48yb59+xg5cmStnr8xaNeuHYCUAwm6KW14eDiDBg1i7ty5RERE8MUXXzBt2jSuXr3Kd999V+Pj/fjjj7zxxhsMHjyYrVu3YmFhUSvXAbrUuG3btrFz50569uzJwoULjdrv7t27rF69msLCQuRyOU2aNKF169YIgjVbtuwgNdWJ/v0tadUK3Nwy8PfX0q5dP1JS3IBEgGJlePGhyiVWVSJa2T7GwBgfacnrHTt2Ak2bRvPxx/czF+7cuYNSqeTgwYPS9kqlEo0GsrIcsLHJQa1WG+SqWlgUVfvFUO8NaX1lpHq4ubnx6aefsnz5ckaNGlWrY0dGRpKVlSW1DBYEwYB91qTM08TEhDZt2gDw2muvkZWVxaJFi8pVM38Y0BvSqKgounXrxvz583nnnXeQy+UsXboU0PkGjx49yjvvvMP8+fPx9fXl9ddfr3Ls7OxswsLCyMnJITc3l+joaBYtWoSzszN79+6lWbNmdO7cGZVKxezZs+nRo0e1zl3Pbrdt28a2bdu4cuW+xERFU9LysHz5cun/arWa69evc/36dcLCwsjNzSMoqA/Tpi0mKupxUlO7cufOBZo2vYmTU6K0n1wuUlQkN0jzqUvUNPpf0T7ljVeVjzQ31wJbW2f8/f2xtbVlzZqV+PnpindatGjB008/jbW1NeHh4dI+ixYtwtzcnM8++4yiItBqFZiYmNC/f386l0i+re611XtDWp8ZqR4TJkzgo48+QqlU1moVRmlGWtu4ePEi27dvB3Q33t8BvUshPj6eb775hnfeeYeXXnqJlStX8uabb7Ju3ToKC3UVKJ9//jnR0dFMmTIFHx+fMv5jPWJiYvj2229Zu3ZtmeAaQGZmJlZWVty+fZuDBw+iUCh44oknyjWkeknD+Ph4EhISpH/67wUFBZiYmNCzZ0+mTJlCYmIi8+fP56mnnuLEiRO4u7vj7u4uaR9UhqCgIA4fPiz9zf38/KT/29vn0L17qPSAa7UCIBo88MYa0dpIgaps/9KqT1VBoxHKbRFSFSM1M8tBJjOhX79+UkBSbzSvXbsmBVL1LhAnJyeCg4Oxs7MjJycHc3NzRNGGPXssOHfuKwNDWl3Ue0Na3xkp6CLffn5+3Lx5s1YNUmlGWtvw8vLCycmJjIwM+vXrR35+vlEPfG1h165dvPzyy7i7uyOXy5kxYwYvvPCClOLTokULTE1NWbhwIQMGDEAmk7Fx40a6dOnC888/z9mzZ6Vk9Li4OHbt2sXOnTs5efIkHTt2ZNKkSQQGBuLt7c20adMIDQ2VFLsEQWDLli3SLMLBwaHM+YWGhtKjRw+DVBkzMzMaN26MtbU13bt3p3///gQHB+Pi4oIoivj6+tKpUyeioqKkF4B+/IkTJ5bbwsLb25v09HQCAwNxc3Njw4ZNgJbw8HDy8vLo1asPJiYY+PNMTGrulKxrYZLqji+TidVmpKII5uZm0ousPGg0GmQyGa+//jpZWVl4eHhIs7j7qmhqMjKW0qGDzmjn5FhjaanA1FRbrS6i9T7Y9E9gpAA9e/as9bYj5UXtaxMODg5MnTpV+o0dHByk3jV1CVEUmTdvHs8++yxt2rRh6tSpzJs3j2HDhrFx40bJfTFz5kzmzZtHbGwsixYtIjExEXt7e3bt2oVSqaRfv37MmjULHx8fWrTwYdasb9Fq2/HOO/MYMmQIVla2nDt3ifh4U0nUOSAgQHqYnnnmGUlSsHS5b2pqKv369SuTbyiXyxk9ejRdu3bl3r17LFiwgF69evHWW29x9OhRoqKimDBhAubm5lIuryAIZGVlkZCQUO7vMWbMGKZPnw5A8+bNadFC9+L08/Oje/fuBt0sazNCXl9QEUOuLGqvk/dTY2VlhYWFBe+//36ZbVJSUgCwtramUaNG5brCdG6fDDp16oJKJcPOLh9TU23xeRlvHus9I42MjGTChAkEBQURFBRU56WZNYWTkxNFRRXX79YEdc1I9Zg4cSIRERHs3r2bqKgounfvXmfHKiwsZOLEiWzYsIGRI0cSGBjI5MmTGTx4MJs3bzboe7527Vpu3LjBK6+8Qk5ODqtXr8bJyYmIiAjAhGvXxrNwYS/MzN7FxUWkX79DtGoVS2amBXv2DObOHS+effZXZLJfpdShkpg7d67k1/T09JSWq9Vq2rdvT35+vtQRE3RGVM8oGzZsyLBhwygsLOTHH3/k22+/Zdu2bVhYWDBy5Eg2btwouS30+x84cIQuXdrxySfw9tsYHG/Pnj1SZVJamheurkiMtHfvPv8641kSFQstV+4jLclEzczM+PDDD/nqq6+kmUBpLd7yEB8fj59fF+RyobhpnxmiqEIQRGQy4ylpvTekTk5OUlLy+PHjWbx4sVGBhocNS0vLavV4MQZ17SNNSUkhNDQUa2trQkND+eSTT+jQoUOdHS8jI4OBAwcSFhbG3LlzGTt2rKRAtW7dOszNdWknoihSVFTEoEGDSE9Pp6ioiJ9++gk/Pz86duyIlZUVvr6LSUvrSJMmJ2ne/ChOToWAGY0b9+XFF7vwyitypkyBnTuDefbZ/9G5c2dMTU1RKBRYWloSGhrKZ599hrW1NS4uLgb+0WHDhnHnzv2OmKDLdli9ejV5eXnSOtA9zGPHjmX58uWkpKRgb29frDSfVub6c3MzUShg9mzdd70xzczMNCjvdHXVBZH0PtJ/sxGtDFX5SEtP6dVqtfT3sre3Z+HChYwdOxZvb+9y909PT+fIkSP88MMPfPVVZ9q1c+fbb7V89dVKlMpcTE3/RelPCoWCLl260KVLF9LS0pgxYwY3b97kww8/fGj+vAsXLvDss89y48aNCiPllpaWtZ6bWNeMND8/n4sXL0rfx40bV+spXCVx8uRJwsLC6NKlC/369TOQ8Rsy5CSXL4ejUHyNUqljctbW1nTo0AG1Wk1ycjLJycnExMTw9NMv0adPApCAUqkkKSmJ8+cTiIiIoEWLq2i1rkRFtUAuV6JWv8mKFavp0KEDQ4YMYcWKFYwaNYpXXnlFqlf/5ZdfsLOzQ6PRMH78ePbt2wdAkyZNSEtLIz8/n8mTJxMXF2dQqltUVERhYSHW1tb07fsif/3VCYXiMAsWLCjXb1fSr/npp7nIZAfIyYknJ8cKKysTaUqph56R1qZC/D8J1a21z8rKkmaF+tlAWlpahYb06NGjmJubM2LECLp1g0aNdGWtL788hPXr11FUZHxL5npvSEuWMLq6uvLGG2+we/dufv31V86ePWt0+40HQceOHQEdO05PTy/3IbGwsPjHMVKdP64F165dA6jVogJRhLAw6NgR9BlbQ4YMYf78+cyePZtXX321xNazCQ0dCgzF0vJNWrb8idjYD7GzsyMwMJDff/9d2vL69YZ8991M2rTZx717H5CaehtHR0dSU71wdNxPRMST6BrZzkWlWkSHDo25cAEaN27Mxo0befnll5kzZw7x8fEAfPXVVwQEBBAXF0dgYCCpqTq5OUtLS65fvwn0pkkTRxwcHMjIyMDR0VESjwFd616l0gx4DlPTJxk82BUTkytotdoy+Ykg4u6ehCDIGD9+Jbm5OrZpb1++S6hk1P5RRHW1LBo2bEhMTAwvvviixPBLT+/VajWbNm2S/v76bgIlGszStKkuj1e/jTGo98GmkopGoAsKDB8+HGtra3744Yc6P75eXBZ0b7yKsghiY2PLZXObNm0iJCSEmohTl65sqgsMHTpU+v+kSZNqLZ80Ph4CAsDUFPLydMsEQeDtt99m69atXLlypfj36gV8iUyWSHDwWRSKG8TGvgdE4uTUguPHjxMdHV08qicazVZMTVWEh79FWpou2drV9f9o2vQX+vW7Sbt2P+Hn9yzwEWZmuahUKry9vbG1tSUuLg6ZTEZeXh6CILBixQpmzpxJYaGaoKB3uXvXAWgLzEShWAekAofo0GEEGo1JcYaDgFp937jJ5TIsLTV07BjOtGlL8PG5gijqpp2lU+EEQWDUqOX4+kYgijpmq1JVrLIcHh4udZJ9FFFdhXw7Oztat27N28U+k0aNGhEZGWnwDBcUFBgYyIqybHr27FmtGWa9V8h3d3cX33jjjTLL79y5w/fff8/58+fr1K/Xp08fDh8+LH0vr5unKIpS/XPXrl0N1n333XfY29tz7do1fH3LW+gJPQAAIABJREFUl+6qCCqVCq1WK/kO6wKiKLJz506pOVpISAhz5syphXGRos1t2sCQIfDMM1BcvFT8G5oC3YE16Fp+gc54OQNyevQ4QLdux1m61Ib8/PeALECOmVkPlMr7SdZy+XLef/8uJiYioaGhHDx4kJYtW2JmZiZVTJ07dw6lUkmvXr3o2bMnaWlpLFu2jD//hFdeySYlxR5r63z8/MLx8wtHo5Fx9Wor0tJceeaZnZiYiCQmehpIrOXl5VUqWFNRNDo31xpBELG2LqjU//kw/v71GdVVyLexscHc3JwDBw6wd+9eafmQIUPo1KkTe/bsIS0tDWdnZynf9NatWxVmBu3evZthw4bVXs+mvxOlGakebm5u+Pj40LFjR4YMGcLChQtrPXXnzJkzBkZULpeX6yP95ZdfKCoqolOnsr/3008/LeUoVree/WEwUkEQeOaZZ6TvISEhnDhx4oHV8wUBfvxR9/+YGFiwALp3h6FD4eJFuHAhHA+PcARhCzAfaImFxRTgDxwcdEwgPPwYx46d4OWXnYod/w6ANUrlH8At4E9sbW3p2FEp+R87dOiAlZUVGRkZUkT+1KlTzJ8/n1GjRnH58mUEQeDMmTP8739R9O0LWVlKnntuOzNmfE3fvn+i0cgwNy8iMPAEzz23QxrbyytZyi28fNmzRkYUwNY2HxubAoO80PLwIIy0Nuvf/y5Ul5Hm5eWRnp5uYEQBfvvtN+7du0ejRo1ITEwkPDwcExMTWrRoQdu2bVmyZEm545WcrVWFfywj1ePevXvcu3eP6Oho4uPjmTp1KtOmTcPNze2Bj33mzBkCAgKk7+W9vSIjI2nfvj2vvfZapTmvKSkpODs7V6vy6WEyEpVKRWZmpkEd+/r163n55ZdrPKYowqZNMHkyZGWBpyfk5EBuLkyfDl9/XXLrK8DbwG/Y2s4mL28uoqiiefMrvPzybsLC/Nm3b1DJ0QGBkJBPDI4ZHh7Onj17MDU1NUhHGzJkCH/88QceHh4EBwezYcMGvLye5uTJz3F2vsOrr67mzp2GuLndb8khilBYaEZmpg22tkrMzIowN9cZ+SNHfPHw2IKtbRfc3TNr9NtUFY3/j5FWzUhL/o6iCJ06dWT69OmSaAnoKsb8/PykF1/JXlgZGRmsW7eOQ4cOSbGQkhAE43o21XtDamNjI75dMumuEmRmZhIaGsr169f5/vvvee655x74+Ppqn/KYqFarpUmTJiQmJhJibE+CauDcuXMPJY+0JDIyMqQ39IABA9i/f/8Dj3nrFkyYAAcPQqdOYGsLeqIfHAzDhqkZNSqZwkJvBOEwopgIBAPfA0mMHdsCMzMlly49zpNPHsfCogiNxqQ4yKPE0vK+z3Ljxo1S/qepqSk3b940CNgEBQXRrVs3vvzyS2xsbJDJlpOR8TzvvvsFMpkaEFCp5CQlmePuno1a7YAo3sPGxor161/G1jaXF17Y/sC/iTEICwt7pKP21b3+5ORbeHk9xt69ezlz5gwnTpxg6tSp2NjYVDrG5cuXJT0HHx8fg3XGGtJ/VNS+Kjg6OjJo0CCSk5OZMGEChw4dYsmSJcjlcvLz88nNza02Uy2vpA90vsVp06YhiiLvvfdetcY0FnUdtS8PTk5OvPnmmyxfvpwDBw4QFRVFy5YtH4gVNWoE+/fDDz/AzJmGTOzUKVizRs7t2y40aPA/1OrPsbGxLg5QBQMa7O1XkpbmSlDQn+TlWWNhUYRMpsXSUhc5V6nkmJrqMiZKC8dcv36dX3/dhULhA3SjWTMHbt1SoVROIiMjDzOzfsjlSonZaDQm3LzZmJYt49BoTJDJcklNdWHnzme5edOb7t3vTzVVKhG5XKiSWda0tv2/qH3lUXu1Ws3WrVvp378/jo6ONGrUCI1Gy5kzZwCdOPjQoUMlPdeK4OPjQ0xMDOfPny9jSI1FvTekFflIK4OnpyevvfYa8+fPZ+3atcUVC6KUsjJy5EgKCwvp3LkzTz75JN26dauWklJubi5jx45l+/btTJ8+vdpGRqPRcPr0aVJTUxFFkcGDB5cr51adPFKtVotKpaqVaWDJqdQTTzwBgLu7O2vXrqVfv341GlMQ4I034KmnYMoUnWEFGD5c9+noaMOsWVZ88UUUeXl5gBqZbAUWFnE4ODhjZlaEIIjY2eVJknF6yOX3085u3brF+fNXyc1tgoNDK/r3b8prr81h27bnycuzxc1tCceO9UKXLQBKJTRtGo+ZmW6MuLgmtG4dhyiCTCYgipCR4YGjYxYdOmzn8ccvSccyNTXunqlpQv1/eaTl55FevXqVTZs2Sd+vXLmCs7MzLi4uWFpaYm9vz5NPPklkZCTz5883qoDHxsbGoNCiuqj3U/uqfKSVIT09HZlMhpWVFWZmZqjVai5fvoxGoyE3N5eYmBhSUlLYsmULL7zwQpXjFRYWsnz5cj766COaNWtGr169KhVRLg8nT57kjz/+MFg2efJkXF1dy2xbHR+Z3rUwc+bMB24fsmrVKhITE8tdp79fVCqVQTlndSCKsGMH/O9/kJysM6oDBkBycjJeXtHY2TUlL68h3bufoHv3k9L1Z2fbIperDdpBlBxTEODw4TiOHv0WsAcK8fT8kxEjzmNre/8+X7VqLYmJmYAtsApLy+7Mnj0fQdAdw84uF5VKRmqqHY89dt//qVBYYGFRyM2bligU4bRp07pG128s/vORlvWR6rN1SsPP70VatBA4ePAAI0aMwMPDg1WrVuHs7MywYcOk7c6ePUt2djatWrXisRLJo1FRUaSmpnLs2DEDUvWvmdrXhJHq4ezsbPBdLpdLDAugR48e/Pzzz5K4QUVQq9WsXbuWGTNm4O7uzpgxY2pc869PEPb29qZJkyYEBARU+KDUpLLp6tWr5WYPiKJoNOsOCAjAxsYGZ2fnMlkDgwYNIikpiejo6BqnngkCPPecjp36+MA33+gMqaenJ15eN0lMlGNiksG9e/s5ceIiQUFBiCLY2+cWX0v5LO/YsZscPfo9DRvm4eiYwpUre0lOnsUPP9jQvHlzqWgiKekGumBVFjAXpfJzaVxb2zzUap2r4PbtZjz22DlpfL0rwdtbgVbbmsTERnh53ar29RuL/xipISMVRVEyov/73zccPboWmUxGmzZtaNasOVeutKJ5c0u2bt3KhAkTSExMlAhBRkYGgYGBUkRfn7Koh4+PDxERERw7dqzSnmcVod4b0rpuzubk5GQgd1Yax48fp0ePHnh7ezN8+HCDH78maNu2rUG/ncpQHR+pvkLp3r17gE7V5tixY/Tv35+DBw9y7tw5owNiPj4++Pj4kJWVxfHjx5k1axaWlrob9I8//pAquPr16yfVM9cENjYwcSLMmQMJCdC0Kfj7N+fmTTmieBEvr574+elUmUQRlEpD4WKtVpBSk5RKOY0a9cbKCl55ZT3W1vmsXj2WxMRPyc/PoXHjxmRlZREbG1sqtWsJGk07FIpjWFkVIggiKpVASoo7TZrc5O5dV2xtc7C0NKw+MjGhTo0o/OcjLe0jLflSz8uLNkjbA2jT5ioeHi6cP5/BkSNHaNy4MTdv3iQ2NhY/Pz9UKhWtW7emWbNmZURsZDIZLVq0oFevXqSlpeHi4lKtc633hvRBGKkxuHHjRqXK6PqSweeff/6hd9ysDiPVl3mePn2aAQMGEB0dzblz5/D3969xtZKDg4OB8R0xYoT0f72yvr6ZXkXtn6vCCy/oDOnJkzpD2rp1Q0DD0KGxqNV7OHEih6CgIExM7gsX69mo3oiKImi1JjRrlsDs2V9JY3fpsoDExCD8/T3x89O5YFq2bMmBAwdISkrC3T2IlJR2mJp+TVJSI1q2jEcQwMxMxa1bHvj7n+fgwac4fTpAYsHm5kU4O9+lU6fz+PpexMSk7nrF/8dIDRlp69at+euvvwDdc1seqVGrdS+8s2fPMmvWLC5cuMDq1WuwsrLE2toaR0dHg6l+Sbi7uwM6nd4FCxYwefJko8+13hvSujReBQUFZGVlVTo91ftO169fz6RJk+rsXMpDdRhpq1atuHr1KgC///47HTp04Pfff2fr1q2SEtFff/31wA9lYWEh+/btk8ROunXrVmMjCqDfVV+Pb2enAWSoVHb4+bVHq9VSUGCOmZnSIMBUEoIAFhbKMgZNoXiM8eN9eeyx+2W9Hh4eTJgwgdTUVHbvdgW0qFQLiI5+k1at7o+nD2AlJzdCLlcTEHCahg3v4Op6D2fne8jlWoPj1wUedUZqZmbGsWPHuHbtWrGrSRdHkMnalakg1KOkQLeZmRk9evTgxInHGDVKN53v3Lkzhw8fpnPnzpibmxvYl8aNGxMSEkJGRgYff/wxgwcPNvpc670hrUtGqm85UZUhGDBgAAcOHGDXrl0Vvs3qAtVhpAMHDpQM6alTp6SqnpJybseOHePYsWOAzj9srFG9ceMGYWFhXL582WC5tbU1V65cITs72yjtx/KgL4O+/yfIBFy4fTsPQYggJ0fHSKtCeZ0tO3UqBG6X2k63YcOGDXnlFViy5DoFBXeIipIzZIiAmZkSQQAfH921Zmc70KTJFYKC/iIry57MTAdcXR+OcXvUGenOnTsBsLAYybVrX3P9uq7t+YgR7Qx6l+mRmJgo6TIEBgZKQSp9ZoggCCxevJgOHToQFhZGixYtGD16dJlxnJyc8Pf3p1mzZkafa70vEa1LRurh4YGlpaX0B6sIelm1kk20Hgaqo5Dv4ODA+PHjpYBYRkZGmW1KBt/0BrUqZGVlsWbNmjJGdNKkSbz99tvY29szY8YMo8YqD489BtbWuig+QJs2et/UdUkhvSaJJVWVXwKYm8PEib9ha3sSGENRkTkajX6droLJxESLpaVu4dq1wZiYaB6aPqheIf9RRZMmTQDw87vB9OnTee+90YSEhNCyZctyt1+1ahVhYWFs27aN3r17AzpXTsnkEj8/P1as0PWrL88Y69GhQ4dqzQbqvSGtS0YqCAIBAQGsXLmyyu30/dcf5lSrurX2Xl5e2NnZGUjUOTo6MnnyZPr374+JiQnu7u74+fnx7rvvVjneoUOHWLRokcGyWbNmERISQoMGDQBdwvuqVauwtramadOm1a7Rd3CAqVNh82ZdTX6vXrp+4oLQiIiICI4fN77WujSqMniCAHZ2uUyYEIaZmQuCIEouBv1ljBr1M3Z22cXnmknjxskPzZA+6upPFhaPA7B9+3ZJMLsy+Pj4EBwcLPX3At2sszT0XXR79epV7jhRUVH8+OOP5aYkVoR6P7Wv6wBPXl4erVtXnQ/YvHlzEhIS0Gq1D613eE0rm0rePGPGjMHR0RFXV9cK/UoVQe8H7dmzJ1euXKF3794GQh2XLt1PTi8oKOD69euYmJiQlJRk0LqjKkybBl9+Cb/8Ap98AsOHqzh3rjvt26fUuUK8XK7BwSGbV1/9iU2bXuS553bi7JwhMVpX13R69tSx96Cgv9BoSroh6haPso9Uq9USF6frgWZsOlLDhg3Zs2cPy5YtQyaTodFoGDhwYJnt9Ey0Iv3gbdu2Vft8H2lGCmBlZVXuW6s09K1aH5YRhZqrP5VMYC6thGMs9u/fL6lpFRYW8vrrr5d54eg7ePbs2ZOQkBApMPfiiy9W61gNGujk9fbtKyIzM5OpU60oKHCQGOnDqBlxdU3n1VdX4eKSUZwFIKBUyosDT1pEEdzdb1H6z1+X5/YoM1JBEFCpdNKOSqVSEsquDA0aNCAjIwMTExPee+89pk+fLrkHQGcg//zzT0mIqKSwSUmEhIQQEhJCdQqB/tWVTcYgNjaW+Ph4nnrqKTIyMjAzM2Pu3LllRJrT0tJo3Lgx77zzTp2dS2k8SGWLKIrk5+djZmZWLcUpwCDlqaKqKz3y8vKwtraWgjgHDx4kNDRUUpM3Ft9+CwcO7MXf/xwvvTSS556z4Nlnf0SjEbGwqN75Vxepqa5kZTnQsuU1Ax1VPUorDOnZal1P8R/1yqbIyEvs2LHFYNkHH3xQYUVdeno6S5cu5fbt21IqE+iIgJ+fnxSA1cPT05OJEydWeg4hISH/VTYZg6ioKNLT0/nmm28IDAzk+vXr3Lt3j7Vr1xpUAmVkZKDRaMjIyKh2WWhN8SA9mwRBqFQvszK88847yOVyo0pASx5Do9EQGhoKGC/Gq8eECbB9u06G8JdfNjFu3Ivs2xeLSnWHvn2rjto/CBo0SMPKSoEgQEGBOYcO9eOppw5JlUwlDab+/w/DT/qoR+2LivLx8AggLe08KpWKpk2bVuondXZ2xsvLi8uXL0uGVBRFJk6cSH5+PrNmzcLa2hpRFFEqleXqW9QU9X5qX9c+UplMRnZ2Np6enpw5cwZXV1d+/fXXMlHqli1bolQqWbp0aZ2eT0nUdV/7imBpaVmjOvo8fU8RKlbNqviY8Nxz7dBodLekQrGVhg370b17xcUStQVB0JWGAlhbF/H007+RlWX3twsjP4pR+5K/efv27enTJxiVSsVjjz1GYGBgla61oqIig5f4woUL2bBhA0lJSRw/fhxBEDAxMTHKiF6/ft3o8673hrSuGekzzzzD448/jkaj4cCBAzz++OOYm5tz65Zh+Z+enT5MV8jDUMivLeTk5LB9+32dzvLSr6rCq6/KOXtWFxzQarXk5Ozi1Cnj0rRqE4IA7u53//Y2yI+ij7Tkb37hwgV++WU6AElJSaxbt66M4E9JpKamYmpqavDyadu2LZMnT8bb21uS1zMWJf2rVaHeG9KHUZb51FNPUVhYSGhoKN9//z0ZGRk89dRTZbb79ttvAcr4WuoKfxcjrQlEUeTmzZvSd30pX3Vgbg7vvtuBu3ddUaks8PRsQWBgIBkZdrV5qv8YPIqMtCTc3T1Qq+/rYHTs2BE/P79ytxVFXb+ukSNHGrDWAQMG8O233+Lu7l5jtTJjUO8NaV0zUtClQzg7O0tCzRX5YfS1tw+LJf6TGKm9vT0zZ86Uvm/cuLFG4/TsaYK1dX9MTQs5fDiUo0dP4OCQU2m3zX8rHkVGqkd09BVWr14FwPDhwwkJCWHo0KHliokcOXKEzz77jIiICNq1a2ewThRF/vjjDy5dukRwcHCdnW+9N6QPSyhEX2/fu3dvkpOTy91G3/e9m74VZh2jPEaak5NTrTaxDxO2traMHz8eoFwpP2MREtKMGzea8/TT/dFqdVH0goKyra7/7XiUGenWrb9I/y/P356dnc3cuXMJCQnhyJEjaDQaWrZsafC87N27l0GDBjF69Gj69OlDo0aN6ux8670hfRiMFHRVQR999BEAgwcPLtORVKPR0LZtW5o0acLjjz/+UM5Jz0iLioqYN28e27dv5+uvv2bu3Ln1MlFbpVJJyuWzZ8+u8ThWVtC69QAuXw4nMzOXzMx87Ozy/vbgz8PGo8xI9c+Yq6srhw4dkpZHR0cTEhLCN998g0ql4sCBA2g0GkRR5OrVq3h7ewNw9OhRhgwZQm5uLq+88gpNmjQpo8erVColdbcHRb03pA9Tuk4mk9GrVy9MTU0NknEvX77M008/TUJCQrm+07qCnpFev36doqIiqdII7svm1RecOnWKZcuW0a9fP7Zt2yYJqNQUkyY5U1j4Bs89N5DTp3UBp787+POwoWekj9oLBJBKkNPS0khOTiYtLY0LFy6UqTrq0KFDuZF8PdE4efIkX3/9Nd988w2HDh2SZpUKhYJFixbxo75n+APikc8jLQ1BEOjbty/Lli2jsLAQjUZjwEAPHjxYYyFjPbRaLQqFosoUocrySGuq0F9byMzMZPHixbRp04bMzEzu3LnD8OHD2bJlC7t37wZ0UVT9A1Fd2NqClVUCv/9+nief7E1oaAB+fhFSbuejgMrySLVa3Yvl3/py6dmzJ+Hhl8jKugvAsmXLymxTmQBzr169ePXVV4mNjWXq1Kk8//zzhIWFkZWVhY+PD7/++itgqLH7IKj3hvRhiymDrkIiLy+PoUOH4unpibu7OwMHDmTVqlX4+vo+0Nh5eXl89dV98eGPPvqowuBWyVp7BwcHHBwcuHHjBqCTDCupvfgwoHcxlETnzp25cuUKJiYmXL58GTMzM6nXTk2NqB6LFk2gb99eODhsQ6sVmD//bYYN20779peq3vlfgMpq7R9ipfLfhkmTJrJ3bxaFhXfw8GiIvb0fN2+msWRJMAEBjahEvAlBEAzYZnx8PHv27GHevHnSi/7555+XBEweFPX+z/GwGSnojPezzz6LRqNh06ZNtG7dGi8vLwRBqHY755JISEgo07irsvrhklH79u3bk5ycLCUU16XjvDxoNBoWLFggfd+wYQOiKLJ69WqaNm1KdnY2w4cPl3pilfRr1RSHDq2lqGgTMTFN6dr1DNbWeezZM5TCQl3C9b99yluffaQP47ePiIjA3j6al156gtat29G8eTaDB2sZNMgOc3Nd3y9jEB4eTrNmzfjggw+wtrbG1NSUtm3blonwPwiMYqSCIBwBAgD9U39LFMVWxetGAfMAF+APYLwoihnF67oAawEtMEYUxXPFy28AlkBTURTzi5dNBF4WRbFXyWP/HYzU3t5eYp4BAQFS/tmcOXNqPKYoiqxbt85g2Zw5cyptSFeSkfbs2ZOkpCTi4+OZPHlymcZ+dY2kpCTUajXDhg0z0G9NT09n06ZNdOrUCYVCQXZ2Nk2aNKmVzIbg4GCSkzWsX59OmzYJtG17idOnu7Jp0zCCg7eSkWGHi0v2Ax+nrvCg9fj1Wf3pYbgULCxeJDS0N2FhDVAqLfD01CAIyUydaouJCZQop68QoihKGTljx46tUDeioKCAZcuW4e/vj5eXF/b29mUCzpWhOlP7KaIo/lRygSAIjwM/AIOBC8CPwHfAS8WbzAWGFP//O6BkpEYO/A/4vLKD/h2MtCRqSzBCP93VoyojCoY+UkEQJBX66jbmelCkpqayZs0aAgMDDYxoQUEBrq6uuLu789RTT5Gamsq1a9fYv39/GdGXmmDt2rXExSWRmqrLprCyuo4gdOHmzZ2cOtWNwMD6ydb0eFBj86jX2uflbSMzMx74PwD27JHRq1djrl+HDRuMG2Pbtm3Y29szfbquQiouLo4bN26gUCiwsrLC09OTLVu24OrqSn5+PkeOHAF0MQh9BoAxeFAf6WhgjyiKxwAEQfgIiBEEwVYUxVxAVvxPKOdYC4DZgiB8J4pihd3Z/g5GWhc4ceK+QLExRhTK6pH26dNHkvN7WNAb0YEDB0qpTSXXgY61y+VydhTL3A8YMKBWjh0cHIyVlYbNmy3IzbXG0bEIJ6cY0tP7ceRIe3x9I7Cxyf/XBlyMZaQPQ4nKGNT2eXTq1A5f3wZ8+mmRRGh++w2eeAK+/lrXONHLq+L9x40bx5o1axg3bhxKpRKZTMbPP/8MwJdffsm8efMk19kHH3xA//79ad68OfHx8Xh7eyOTyYxuYV4dQzpPEIQvgKvAB6IoHgEeB07pNxBFMV4QBCXQEjgPzAH2AhpgbKnxzgFHgFnAhxUdtCQjzc/P58yZM5iamhIYGGj0RdYHlKxQMva8S0ftbWxsKlR0UigUXLlyBSsrKxISEiT1JlEUsbOz4+7duyiVSuzt7fHw8OD69eu4ublx9+5dnJ2dadWqFaampty6dYuEhAS8vLyIiori4sWLfPTRR3z44YdlFJ22bt1Kx44d8fX1JSkpiYyMjDJiLw+CtWvXcvhwEjAXrdYVL68M0tOjgeGo1Sbcvu1Bq1b1Kw2sNmEsI63JY1AXxre2x4uIiCAz8wyJic/RooVOGaxbN/j+e5g5U/evIl9tVFQUa9aswcnJib1793L37l1p3e7duxk6dChTpkwhOzvbQHIPqFavJj2MDTa9AzQFGqGbvu8RBKEZYAOUdlJlA7YAoigeE0WxuSiKrURRDC1n3I+BtwRBqFDw0tLSEo1Gw71799ixYwfnzp3jzz//5LfffpOW/xM+9dPyPn36GL1f27Zt8fHxqXK7/Px8vvzyS3bt2sWxY8d45pln6NKlC4GBgfj4+DBy5Ei6du3K2LFjcXNzIzk5GX9/f1544QX8/f1JS0tjwYIFfP7556xdu5bjx49z6dIlLl68yIwZM/joo4+4ceMGarWa2NhY6XPz5s14enqi0WhYuXIlLi4u+Pj4cO7cOU6ePFlm++p+jh49milTXkAQ1OTn22Njk4mrqwM6l3sszs4pf/vftS4/n3jiCVq3bl0n46en//3XZ8z1u7l15ZtvtjFlyky+/PI3Vq7cygcfpAKxgJorV8reN6+++qoU9MzIyGDHjh2sXr0amUxGq1atGDBgALGxsZiZmZGbm1vpfWgsjDKkoiieEUUxVxTFIlEU1wIngUFAHlBaUcIOMMqxKYpiNPAbUGEDoYyMDJRKJVu2bCEuLk6a6p8/fx6lUsnOnTv/EZ9arZbevXsTGxtr9H7h4eFs2rSp0u127NjBmjVrcHZ2JiIigscee4zXX3+dU6dOMWPGDMLDwzlx4gSLFi1i586dnDlzhu+++46wsDAGDhzI6dOn2b17N40aNeLq1asEBASQkZGBhYUFGRkZnDp1iry8PIKDg8nNzZU++/Tpw/nz5zl37pzUYeDevXsA+Pv7ExgYSHZ2dpn9qvO5YsUKXnppCO3b5/LLL/MoKlJSVDS9+PYKxtQ062//u9bl5/nz56v8+5f3WVRU+XqFQl2j/f6O6z90aDF2drfYtetXMjNPEBIyHTe3azRuHExSUi7jxpW9b7p37y7ZiXbt2tGmTRt++OEH0tLScHR0LPd+rujTWNRIIV8QhP3AfsANaCyK4uji5U2BK4BzsY+0ov1vABNFUTwkCEJzdIGqhUDv0lF7vUL+1q1biY6Opn///mg0GhISEhgzZky1z/2fBGMU0s+f14ne7tu3r0KXwRdffMF7772HhYUFXl5eJCUl4ejoyIwZM/jf//5XaTfFitCsWTMSEhKYM2cOR44c4ejRozg5OUnyeTExMVKUHZhpAAAgAElEQVRrElEUa+SGUSgUaDQa3n/fhpUr4f/+bxlnzjRgy5bhwDYCAxvQt2/1Vab+KahLhfzSqv8qlRwzM3WVU/6S63/+eRSjR28s0z3AmH2NgT5AW9KlJAhvYmrqyKxZplQl5mRsf7WlS5cSExODSqVizJgxBvX6giAYpZBf5VEEQXAQBKG/IAgWgiDIBUEYDfQADgIbgKGCIDwpCII1uvDa9sqMaGmIohgHbAamlrde/1bQ53wdPHiQQ4cO1biP+j8Jxqg/Xb16lVdeeaVSQ/Xuu++SmZnJc889h4WFBbNmzaJv375899139O/fv0KRlopw6tQpEhISAJ2/V19QUFJdp02bNgiCwNmzZzExMWHYsGH88MMPZbIXKsPatWuZN28eHTtCQQFkZ5thZqbb/4knMuq1Ea2NPEt9Hqko3h+v9LhVHaei9aVV/83MjJvGltxv1ChDI1pYWP4LWRTh1i23ahlRtVogIiKCjRtvsmbNu/zxxyTatfuIjz9uwHvvVW1Ewfj+apmZmWzdupWffvqJ3377rUaZQlUy0mL/5T6gNbqg0RXgI1EU/yhePwr4AnAGDgHj9HmklYx5g2JGWvz9MeAacLoiRqrVajly5IjUj3348OEPTTzk70JFjCQuLo5Lly6Rl5fHnTt3uHnzJo6OjhQVFXHy5ElSU1NJS0vD1dWVAQMGSL2T8vLy6NKlC3l5eXTr1o0mTZpw8uRJwsLCmDFjBnPnzjXqvHr27MmxY8fw9/dn8ODBiKIoRZfXr18vVV+Vhre3NyYmJuzdu9eozq16Rrpjhw1jxsD776/jzh01+/YN4I03Vhh1rv9kKJUqrl1rwuHDQ2nXLpqePXUv1er0jTJmG7VaQC4XUankyGRqg3Yqd+86U1Bgjb19No6OZXN29epclR1Pq4UVK17j9deNq2vXamHLlhdIS7Ph3j0nikMuAEyapOs4+/XXMG4cPPaYUUNWitOnT9OvXz/69u3Lzp07sbCwoFOnTjRq1IjNmzcbxUjrffM7Gxsb8e233zZYdu/ePXJycqQH89+Kc+fOlam1X7VqFYmJiQQFBTF16lRat26Nra0ty5cv5+uvv8bR0RF7e3usra3Jysrixo0b9O/fn/Xr12Nvb49arWbr1q0sWbKEiIgIbGxscHFxITY2lpycHKysrCo9p127dvHMM88AMGrUKFq2bFnudvoGej4+Pvj4+NC8eXNOnjzJ8ePHuXDhQoUCvSXx/fffk5SUxGOPzeXNN2HatM0UFWWQkNCErl2rp3b+T4Iowo0bjfnll3SKirKAT+nd+yA9e56WmJ+FhbrcKXV1p88ljXLJxn/Z2TbY2eWh0cj47rs3mTr1W2mfnBzdupI4daor7duHY2VVVguhqEhApbLCxia/3GOXxiefzGLbNmvGjl2GWp2Kk9P/Ud6k6ZNP4OOPjb/WyvDBBx9w6NAhevfujUwmIzExkby8PHbu3PnvaH5XOo80KSmJlStXArqk2d69e+Pm5oZWq8XMzKzGDd/Kw8qVK0lKSiI4OBhvb+8a+fny8vLIysqqsM+7RqNBrVaX6wcrmUeqY2Y7SExMZObMmXz22WdcvHiRF154gWvXrtGuXbtyKzdUKhVz587FwcEBhUKBhYUFL730Ei+99BJKpZLLly9z4sQJ3Nzcquxjk5eXx5dffil9r6xL6IgRI1AoFFJVib5nDmCUEQWdq0Cj0RATo3vg0tIsadw4EycnnVunogdRrRaQycR6kVtZGsb4EVeuHEtycmNAgW4SCL6+urSy0NAAzp71Y/Lk5cjlkJtrxY8/voadXS5t20bTpcvpEm1xKj6Wfp0ggEYjYGJi+HuZmOjuO5lMw+jRG0lJccPF5R4pKW1wc7tiMIZCYUpSUiOaNLmJldXtMsfNz3cgK0tD6UezonMbNOhPvvnmaXbsGEeXLiLW1nD1KqxdC/PmgY0NNGoEs2ZV/DtWhBMnTjB9+nRatmzJhhJZ/Y6OjshkMqmQpHnz5gAGBSiVod4b0pL+isLCQsmIdu/eHU9PT44dO0ZeXh4mJiYUFBTw1ltv1VoSf1JSEqDz1UHlAiOloVQqWbBggSTC/Prrr5fJVwNYtGgRubm5DBo0qEyyfck80mPHjiGTydi+fTtyuZzevXtz5coVOnfuzKBBgyqsJDI1NcXOzo6cnByuXr1qILpiZmZG+/btad++fZXXk5iYKPWxf/fdd6s0uqXFIPR/N3Nzc5YuXcrdu3f59NNPKx1j7dq1JCUlMXfuXKZMgQ0bujNu3DWaN48HKn4Q66sRBdixYyjDhu1FJis/0X7//n7FRhRgNZBCq1ajsbfPQRQhJkbFpElLkcvlaLVga1uAr+9FTpx4ku7dPRkwwJ7ff/8dgJUrx3PvngsuLqlMmLC2DGvNy7NCpTItM2UXRd24oJtm29vnUlhoSkyMD23bRhmMAZCd7cCIEVvJyrJn+/Zn6NbtFG5uurzNXbsGER7eib59/6Bp0/IyIMvC3z+cBQv6EhS0mqFDU9m9+/9o1Qo+/1z3r/TxK4MoiixcuJClS5eSmJgoLT937hzDhg1j8uTJUrYJQFBQUI2Cr/V+aq/3kZ46dUq6QQAmTZpURl1o586dNGjQoNYU7PPz81mzZg1mZmbcunWr0p7apVGyN7we06dPLxMk27BhA+np6Wi1WqZNm2awTu8jVSgULF68mC1btrBhwwaOHTtG3759adWqVZWuDVEUSUxMZPXq1dy5c6dG8ntKpbIMY67OSwV0N+5vv/1W5twqg95HamNjQ34+dOkCKSkZ9O79F23axCCTact9mHRtSQRMTY3PA3xYiI1tTsuWcWX8nPqp9YIFs1Eo9C9FBdOnh2JndwJB0P1W+kh0UZGca9c64u0dTUJCU7Zvf47sbJDJ8vnqq6/QauH2bQ+2bp2Ik5NARkYu8+cf5PbtS2i1OgYKkJnpwKlTATRvHl9ucYNCYYWlZUGF16NWy1GrZRw61JcLF/zQamUIgpY2bS4RFHQYjUbO8uVvMmLEdlq3jjb6d4qMbMeOHQOYMEHkp5+q15G2JKZMmcKyZcvw8/PD0dERR0dH2rZtS2xsrEGlnpubG/b29jz//PMGWQL/qr72Wq2WU6d0BVS2trZMmjTp/9k787ioyv2Pv5+BYd/dUUDELfd918g1Tc0ltdwry8gWu9du2630drv9WswytVJLM69Z7gupuV9NTERE3BUBERCQfR+YeX5/HGacYR0QEXU+rxevYZ5znnPO98zM93ye71oqA+vUqROBgYH06tWrWmynjo6OzJkzh+zsbC5dulSpJ1WTJk1KeMPDw8Pp16+fydi4ceO4dOlSqcvko0ePEhYWRqNGjRg2bBivv/46fn5+vPTSS2Uq9KysLJKTk/Hx8UGj0XDgwAFu3rzJ4sWLq1zDVK1Ws2nTJlasWMHZs2cZPHhwpZQoKB0hjbFo0aIKw1OMGamjI4SGwrZtHkyZ8hS7duXQs+dxWrW6ZGA/Fy605MSJfsTGelFYqKN+/UT69TtMu3Y106zQHPj5XSM62gtHxyzq1k1FCNBoFCdPREQLIyUKvr6LCAnZw6hRo8jONrUvHjvWl0cf9ScqKhEPD8W3q9XCmjWO5OXZYWeXR5MmcXh7h3PmTAfAmRdffIrWrTvg5RVNeroLqanuXLvWDK3WmuDgngwfvoeePY+bnMfePsekNKIxhBBYWxfyxx+DOXu2GzqdkrK5ZYuKuXPbc/AgPPXUZmbMuEDjxrmGebm5dmzZMg5Pzwj8/Uu3dbdvH86WLXHUqZOIPte+Kjh58iSgREC8+eabhhrArVq14m9/+xsuLi5VDs8zRq1npA0aNJCenp6kp6czbdq0cn94UkpWr16Nt7c3PXv2NCyry7Pl3U3ExsbSqFEjkpOTcXZ2NquXtjGSk5P55ptvACU99PHHH6ddu3Yl9tPpdHz11VeGajX6QiK3bt3ikUceYcWKFZVqLVschYWFzJw5kz///JPJkyeXSBUtjoyMDFxcTPM0cnJy+OyzzwCYNWsWK1eu5L///S+TJ08u8zjGjNQYUsK6dfD11xAZCSkpyhJ0zhyYMAH8/UGlkjRu/H+kpW3hjTeeqHX1O8+fb03r1hfJz7fD2rqQ6GgfNm0ab1CkQki+/PJbcnKUMhT673LLli05fTqK2NimODg8w/Xru+nRI5gNGyZw7VprtNpCZs36gTp1bqFSSdLS3FiyZA5KuYuyYW0N165J9u1bbbIE7tixI+fOnaOwsBA7Ozvy8m47k7p06cqpUyEsWTKHW7fqolZDWprSKkarBVtbHQEBy7C3z+XWLQ+8vW+g1VqxdGkALVrU4fXXs7l06YsSn42epaemgkbzBp07u2BnB0OHQmW5wLp165gyZQpApVaUepjLSGvZ16skkpKSUKvVTJgwoUKWKYRg3LhxxMbGsmTJElauXMnSpUvvWSmyxo0bo1KpqFevXqWVKCj1S/Ve9DFjxpSqREFZetvY2BhssDqdjs6dOxMXF0dycvIdV2JatmyZQempVCouX75Mfn4+oGSe6VPprly5wvz58/nyyy+JjY01zNdqtQYlCrBypVJErKLarvo40uIQAqZMgRMnICkJ8vOhWTNITAR9zWxX1/8jN/dLpkzpSnx80h3JXxrulH888shF0tLc+fLLN0hKqkNISGcTNvrmm0ns3RvIiRMn6NmzJ0IIvL29GTFiBA4OGpyd01Crf+bmzYbExzdiwoQNtGhxlqefXk+jRjexsSlEo1FTt24yTZtGlXkdjRpB//7w7bfg5SUYPXq0ye8sLCzM8PnqlaibmxuvvPIKarU1hYVqkpPr4OgIK1cqShTAygoCA1X88sszFBZa07hxPADW1loGDjzIiBGncXLaj0bjbLiXWi1kZt5exl+5coLDh5/lxRdh+nS4cgXi4+OJj483+z63atXK8P/HH3/Mnj17OHPmTKXimc1BrWekDRs2lAEBAVWaGxoayoULF8plPbUVUkp+/PFHYmJieOSRRyrVEmHFihXExsYalmS9e/c2mEaqgr/++ot+/foxaNAg9uzZYxhXqVRmPaTs7e0NvXKKQ6vVlvmALIuRFkdBgdKaZNgw2L4d1OoMrK0bcejQQTp16sGSJYvJyEitdgdUdRb+kBKWLJlDcnJdPDxS0Gic+emnA5w8eZhmzZrh5ubGpUtWZGXFYGeXajI3L88GjcYGZ+csw/VcudKMyMhm9O9/lPj4Rqxf/zQajQ2ffgqvvgpbt8LkyXDmDBSvb3zmzBkOHDhAenrptV49PT3x8/MjNDSUunXrMm3ajDIZ/65d8Mwz6UyZ8l/q108q9Z5JCTqdyuCA0++jV3Z16jQgNzcbISQFBQV4enoya9Yss5fjISEhLFu2jPT0dJOeT56enkyfPr1ckmMuI631irR4HGl6ejo7d+7kypUrNG/enKlTp5Y5Ny4ujrVr1zJp0iSDx/l+gL4b57Vr17C1teX111+vML7TGLGxseTn55Obm8uGDRvIy8u7ozRDjUbD22+/TX5+PvHx8UybNg0nJyfOnj3LTz/9RFhYGKCYFJKSSmd/vr6+REZGlhiPj48vk5nq40grShQ4elRhVbcxgfff741GM4WEhP1otRI/v6tmyVoR7mbJuoUL/4adXR4BAd9y40YTUlMzqFPnHO3atUOj0ZCaqgV0RY4iSVqaB05Omdjb55lcV3a2PVu2jOPq1eZ06xbME0/8Tm5uHbZvH8zGja2IixM8/bQSMxodfZtFFkdubi4nTpxAq9Xi7e2Nn58fISEhnDp1ivj4eOrUqcOYMWPKDO3TIzISnn++ECur/1GvXiKPPBKNTqewWzs7exISfLhwwYVevU6YzDtx4oSh+pWdnR2dOnXCysqKHj16lDAdmQspJU2bNjUxX5TmGDba9mAoUr3XHiA6OppNmzaRm5tLQUEBzz//PF4VpDacPn2aoKAgqspq7wVOnjzJiRMnSElJYebMmRV+UUtDTk4OX3/9Nfn5+RV6x6sT8fHxfP311ybxpq6urixcuJBZs2bh7u5Oly5dSE9Px87OrtwUWHMZ6cKFt2MKPTxSSEmpQ9++V+jUaQ8eHumoVIUGr3dVoVdUWq0So1p83Jy5iYl1qV//Vpn7HTjgz8CB0YDywMnOtqNJE3vUajU5OTlkZt5mnHl5tvz22wQiI5vRtGksXbuG4+NzjawsX9zdr2Fnl0x+fhtmzBiGjU0yW7duJyMjjYwMDzZsGIuraxO2bIGicMlKQ29Oqgz0kQlCSJKSkrCxscHOzhV7e8GLL36Pp+fNEucAGDlyJK6urrRp06ZqF1sM+lXUwIEDOXz4MI8++iiPPfZYqfs+MDZSfZmrY8eO8dtvv5GXl0dBQQG9evXCw8OjwvlxcXE4OzsbbHq1HdeuXWPnzp0kJibSoUMHLl6smsf5hx9+MMgcHR3Nt99+y4ABA8ya+8svvyCEKDPVszyEh4ebNPdbunQpaWlpnDt3DiEEqampZGVlcfLkyVLjao1Rlo20OIzMscyZY0WzZr2xsoqjXr1bzJgx6Y4fJHoFkJnpaAgbKg1llV2zs3NHSkoo0eKXNXDgIfRKVEo4cSKMwMBAMjIyyMrKQqdTExLSlbNnX+W9994iKMiPlBTBtWtN2LBhOF98MYfvvhvBRx+9xGOPPYaDwyXWr19EYGAg3t5NAG+cnFKZMOFPgoKqrkSBSitRUB4mKpXiy6hfvz5ubm7Y2Qk2b4bDhx8tsf/p06c5evQoLVq0qDYlCopJSqVScfDgQV555RUOHz5cpunJXNT68CeNRsP//d//0bhxY1544QX27NnDhQsXOH78OOHh4UyfPt0Q1iOlJDExEa1WS0FBAVFRUdjY2HD16lUOHDjA8OHD77E0pSMvL4/Lly8TGhpqWP5OnjwZX1/fKjvK5syZQ3BwMLt27eL69eu8/PLLZs2LjY017Pv999+bpciMUVhYiI2NDR4eHgwdOtRwrG3bthkU2qVLl7CxsaF46m9x6DObKkJRoX4A4uOjsLK6raBXrTpiVoGL8qDPAHJ2VkKQ9ApQXzEpM9ORlBQnfHwSSp1fUJCOEL2A4ybzAVJSXJFSUKeOaZMIIaBfv2aoVM1p0KABfn5+vPFGT0JCbPjgAyW7pyyibm1tzYABA2jXrh1nz54lLi6u6KGYhUoFLVqkljm3ppCWloZWq6VOnToMHZrNuXNbMf6onZycmD59On379r1rrXWEEEyZMoU1a9ZUyRlsjFqvSK2trfn73/9uiOEcPHgw/v7+NGjQgJMnT7J161aaN29Os2bNCAkJITo6GkdHR0Oco/41Ojqa5OTkGm8aVxEiIyMNmVMAr732moFpl5Zrby5UKhUdOnRg165dTJw40TBev359vvvuO8aNG8eVK1fQaDS4uLjg5eWFlJKAgAA6duxIdna2WYy/OEaMGEFOTskA7j/++AONRkObNm1IS0vD39+frl27lnss4zjS8pCZCS4ukJEBK1fGAE+Sn6+YfKytoystgzHKKhBSWKgyVEyKjW2Cu3vpdXp0OitAS1YWhIb2oWfPYDQaGyIimlG/fgp2dlb88MNMOnQ4zbBh203m/vijYPJkLTNnzkRK0C9OzK0e6eHhYbIKycjIIDY29p6378nNzWXJkiXY2Njw2muv8csvv6DV5uPo6IidnR0DBgzAycmJXbt2sXz5cv71r6rHkVaETp064eXlRWhoqCGduSqo9YrU1dXVJBDeWBF26dKFrKwsCgsL2bFjB76+vrz66quGZcfVq1cJClLS0hITEwkMDGTAgAF4enpWaWlSndBqtezdu5fjx4/j7OzM3LlzSwS5F+/ZVFnY29vTvHlzrl69SsuWLbl8+TJJSUmMHz++xL6xsbEcOXKE48ePM336dBYuXEjTpk3Jzc2tlkZ2+vYNubm5qNVqVCpVhV5XcxhpZibs2KF/J3F1bUXr1il06xZYLY4h40pIxmNq9e3PJSSkM8OG7TCZpz+3SqUteh/OgQPz2L9/CHXr5tCyZTj9+kWRnGwNCNLSOvP6674EBQXRuHFjPvusPfHxefTooczfvBmys6FdO6hCJwwAXFxcquykqU7Y29vTokULLl68aLClOzs74+/vT2JiIu3bt0cIwbPPPnvX7ft2dnasW7eOUaNGcfHiRXx8fOjbt2+lj1PrnU2lVX+qCBqNhrCwMOrWrYunpydZWVl88803eHp6IoQgMTERJycnpJRIKWnRogVnzpxh9OjRtGnTxqysKCklCxYswNPTkxdeeKFSmRGXLl0ypKf17duXIUOGlLrfnTBSPeLj4/n+++8RQiClLDVL5cyZMzRq1IimTZsyatQo1q9fb9jm4uJCTEzMPfkBVuS112qVQHIFEn//Q/j7/6/Grk+Pc+fa0LZt+b2qpISbNxuSk+NA06ZRWFnp0GobsXlzH86da8fgwbB37+399+yBxx9fRtOm8Vy79hFNmkBcHISEwB0Qp1oDKSUHDhzgwoULuLu7M2rUqBLfsaVLl5KQkHBXGakeeXl5DB06lCNHjph48R+YFFFzliGHDh1Cq9WSmZlJYmIiGRkZuLu7o9VqSUpKMmSFxMXF8cEHH7B//37+/FNp5atWqw3tPzZu3Ago2SPjx48vN2RICEH37t0JDg5mwYIF9OrVy6zumYWFhQYl+sILL9C4ceMy971TRgoYYgHr1KnDrVu3GDduHL6+via2z6SkJPbv34+np6eJEgVlOVhaplJNoCJGqleirq5pTJ26lnr1ku/q9Wi1Ap1OYG2tfCZCwJo1Uxk/fmOFc62sbMnKcqJt20yOHu1BSEgnEhJup+m0aGG6/9Ch4Ob2LFFRWoMSffTRB0OJgvL7GTRoEIMGDSpzn5pgpHrY2dmxceNGGjRowI0bNyodKVPrFWlF1apjYmIMvagBpk+fjpubG25ubpw/f56NGzdSv359evTowSOPPIJKpTLJOS8oKDAo2i5dutCyZUvCwsLYs2cPTzzxRLk55U888QReXl5s3ryZ48eP079/f0Mur5SS+Ph4Tp8+jY2NDf7+/sTHxxuqIH3wwQcVMt/iXUQrC31LFlBquM6YMcPEHvvUU08hpWTQoEG4ubmZ3Gtra2sKCwv5+OOPy1X2dxPl2Uj1Ga8qlZZx4zbddSWal2dDfr4aZ2el/XNhoeDnn5+hc+e/cHTMIzDwSUaO3G744et0EBvblIED3YiIOENamiNHj05h9Wo4fNjUQTZggBLCZQwhYMSIVaxbF09c3EfUqweBgXdVxFqHVatW1RgjhdsZdxVFk5SGWr+0N44jLQ6dTscPP/xAkyZNkFLi5eVlaEkCsHbtWpydnWnSpAmHDx+mVatW9OnTx5B7n5eXx9GjRw095/Xl4dLS0vjqq6/o3bu32QU6jJcDnp6exMXFlbnvwIEDzQpFupOePdevX2fHjh2GAPnx48dTp04dbG1tSU1NxcPDgytXrhjiTM+dO2fSsrZHjx588803JUr71STKiiPt3z+funV3cflyC/r0OUaTJmXf6+qA/idibL3R6ZRBlUrQsOEA3nijJ6+++g35+ba4u4/Aze0gaWk3i+YLtm0bzeLFnejfH9asAaOuLCxdqlR+L47c3Fzi47UkJzvRsSPcY7N+jSMnJwcppYGc3G3cuHEDLy8v1Go17777LkKIB2dpXx4jDQ8PR0rJsGHDTNhdfn4+v//+O1evXqV79+507dqVli1bEhQUxPLly+nWrRuDBg3Czs6OwYMH8+ijj3Ls2DECAwPp06cPjRo1okuXLgQFBeHt7V2itmZpeO+999i2bRvp6enExMTQqFEjJk6ciLu7O5mZmSwsohxt2rQxu0WKOYw0NTWVyMhI2rZtS1paGmFhYWRkZJCamsrs2bP597//DWCSGufj48Ozzz5Lz549AcWmbBwzumnTJsaNG2fWNd5NlMZI586F1NSrDB4cRqdOYTVyHXl5aqytJWp1Ifn5Sixnr15BqFSKZh01qgN//LEbW1sNa9ZMZ8iQRvz0Uwu2bcvg2WchP1/NO+84GLKvjFezn34K06aVft5Vq1YRHx9fYd3WBxU1zUitrKxQqVQUFBQQExODt7e32XNrvSIty0aq7080dOhQEyWal5fHl19+iUajYezYsYZCxs7OzgwdOpTGjRuzYcMGBgwYYKgEo1arUavVhIeHk5OTw/jx4xk9ejQ3b97kyJEjZilStVrNU089ZTIWHx/PkiVLDIVjn3vuOerVq8enn37K448/Tq9evco9pjk20gsXLvDHH3+wfbsSOjN16lRatmxJQECAwSOqdzQ1b96cYcOGGeLy9Dbk7777zuSY5th6awKl2Uizs6OYMKFim2RlUZaHX0o4eHAgI0YoNQY2bRrL1asNmTfPhZs3r5KQkMC5czfo0CEcIR4lO7sRtrbKPR8zxpW0NKW4yrvv3j5m48bmFT159tlnzYqjfVBRkzZSgKioKBo1akRsbCw//vhjuamjxVHrFWlZjHTHjh0kJiayf/9+fH19sbKyQqvV8uOPPxq80qXZOvSBt8Vri/bu3ZsjR44QERHBhg0bePrpp/H29ub48eNs2bKFsWPHVuq6V69ebWB5HTt25IknnjAkB4DSibMiRWoOI9Xbd7/++muGDx9OixYtCA4Opl69erRu3Zo333wTBwcHrly5gru7O7a2tobulMnJt+2KH374YaW+ODWB0hhpvXpr7sq5ygq6SE9PIyVlBjCZjAwnWrfWsn+/O1lZ7hw7ptifjx3bSlpaQxYv7o9OB08/fXv+zJnKX1VgYaQ1y0gvXLhASkoKL730ksG3YC5qvSItjZHm5uZy9uxZHB0diY+PR6fTYWVlZejxHhAQwKFDh0pU0AdlWVtaDKMQgnbt2hEcHExkZCSrV682xE9W1tmyfv16oqKiSu10mpamZLCYU8W/OCP9+eefDV5FPz8/Dh48aNg2bdo03N3dee6551i1ahWg9N2l1UwAACAASURBVE0SQpCXl0dgYGCp1Xw2btxYalxpbUBpjLRhQzdSU1PLmFH9cHFx4caN69jZ2fHDD1+TkPAubm7J/PDDr4AS59yrVy8aNepI3bpWeHnBwIHVc24LI61ZRvrhhx8a4pwr22Wj1ivS0hiptbU1/fr1Q6VSkZKSYlii5+Tk0LZtWxo0aFBm2bmMjIxSg/HT09M5d+4cM2fORK1WExUVxenTp5k+fTrNmjWr8Dp1Oh0xMTHs3r2b+Ph4Ro4cWaotVJ8tZE7WkDEjTU5OJiIigogIpV+RvnrN559/ziuvvIKdnR0XLlwwKNF33nnH8LDYsmWLiRJdvHgxr7zyyh1XBS8LsbHK8vVOURoj9fNrzsmTwXd+8AogpaSwsBC1Ws1LL73EmTNn6N27N/v27SM1NZUnn3wSDw8PEzvae+9V7zVYGGnNMtIGDRrg6elZpezHWq9IS2OkarWawYMHs3v3boPQeXl5nDp1qlx2JaVky5YtdO3alZs3b1KnTh2DEt61axfdu3enadOmgMJCy8tw0Gq1nDp1isBSYlKGDx9Ot26lO/r06ZPmVKzv2LEjcXFxpKenm9QTXbhwITExMXz22WcmFb/1FfKN+ymlp6dz6dIlg/x5eXl3nFdcFt57TynJ9s9/Qq9eSkrj3LlKy9yq6OzSGOnevbeoiYYHQgjDvXVzcyMnJ4f69evTpEkTLl68WGYSRXXCwkhrlpFmZGRUmLZcFmq9Ii3LRpqfn094eDjTilyeSUlJ2Nrallt39NdffyUxMZGYmBhCQ0Np1KgR06ZNQ6fTcfHixRLOImPk5uYSGRnJjRs3iI2NJT4+vkSGkLu7O6mpqWVWsgcIDlbYlDktD4KCgjhw4IDh/Zw5c1iyZEmZ+6empmJvb09mZiZubm6kpKSwePFi4LYCvxtKNC9PsQtu23Z77HhR65/585W/t95S4iWHDzdfqRZnpLm5ubi7l6xpai5Ky5uXEgoKXLCxyShjjiQkJISoqChatmyJn58fx48fp1OnTiVaX1c3LIy0fEZaHb2WjHEntThqvSIty2ufm5uLSqUyFAW2s7Mr0SDMGDqdjoiICHx9fUlPT+eFF15g8eLFJCYmGrzYZTW3i4yMZMuWLdSvXx8vLy9Dvr6UkqCgIC5dusT48eO5cOECycnJ5ca9tW/fnujoigtppKammijRRYsW8dprr5U7p1WrVuTm5vLVV1+ZFFJOSkqqlnx5Y6SnQ3g41KsHrVtXvP+nnyp/+/aZhv+Uh+KMNDf3zgIpy8qbt7HJwMnJiaysrFLmCLp162aywtDpdCxbtow+ffowePDgu2YisTDS0hlpUFCQwYb56quv8vHHH5OWlkZOTg5WVlaEh4dX2jms0WhIT0+vcn+3+6IeaWlITEzEzc3N8D4rK6vcNEaVSsXo0aOxtbVl6tSpWFtb4+PjY9JbqDTcvHmTDRs2MHr0aKZMmUKbNm1Qq9VkZ2fj4ODAoEGDePnll8nMzCQkJKTCTqPdunXj7bffLncfwFDwePLkyWi1WubOnVthJpSvry+bN28GMCjR06dPV2sZMimVNhU+PkpV+rKU6MKFMGIEFK+7PWYMLFhg3rmK1yM9e1ZFdnbVHghKOwsoKChd6ZWmRIv/iNPT002avx07doyVK1dWqW6rOVi1apVJgeyHDaXJf/78efr06UO9evXw9/dn8+bNuLi48Mgjj9CnTx9at27NjBkzEEJQt25dXn/9dTZv3szixYvp0qULQggmT57M8uXLmTVrFgkJCRQUFLBo0SKaNGlSpZ72cB9nNm3dupWGDRsaQoh0Oh3/+te/GDVqlNl2jqCgIFJSUvD09GTfvn2l1scMDAzk+vXr+Pn5GYoTOzs7k5qaipOTEzY2NuTn55OTk8PYsWPNckxVhKioKHbv3s3WrVtp27ZthRXiS5sfGhrKmDFjqpUtFRbqGDcujx07THtTTJwIU6fCqFGKspJSaX62b98+kpKScHWtw5NPDqaw8PaDwMkpk7ffFrz3XtmyFc9sev/9y2i1G7C1rXy/+vPnW9CmzRXS051xdS0/7VgP/dIxPT0drVaLh4cHMTExJboyuLi48MYbb1T6miqCuR0CHlQYZzb9+OOPPP/884ZtEydONBR7TkxMNETo6KNcMjIyyMzMJCoqiujoaFxcXMjLyyMmJoacnBy0Wi329vZIKQ0NKkeMGFGiZfkDn9kkhDBhaCqVimbNmlUqNKZRo0YcOHCAK1eu8Oyzz5a6j7+/P+Hh4Vy/fp309HTeeust7O3t0el0XL9+3eSDqExfpfJw7NgxPvroI8LCwti6dWuliys3bdrU4DS7E8TFFTBmTCo63X4iImJIS5sBNADygP/g6HiAmTOn8PnnMzl9+jTr1kXSrFkzWrRowblz5wyFYb79dpZBiXbv/hcdO4Zz8mRX/vnPzsyfn8LHH8fx5pvtSthOf/rpJ65fv87YsWM5ePAgVla5VJYw6HRKVfbU1DrAFezsFEZZXok9e3sHcnNzDIkMWVlZrFixgtdff72E2Wb37t2cPXv2rihSi41UsZEGBATw/PPP0717dx577DHs7e1NCIJxmKNeJ+jrbXh5edHftKEXoChclUpFfn5+lWvvGuO+ZaShoaFERESYOIj++usv4uPjGTNmTIXHlVKSkpJCcnIyPj4+2NraIqU0LNmzsrLYu3cvCQkJJCcnG2xVAQEBJZ5a1YmMjAwWLVpEbGwsbm5u94SRpKen85///IelS2+Snf2T0RYdDRocJzX1BTSa22XjHB0dmTt3bgkHWny8B40apZCd7UB8fAN8fK6jVuvrc0JISCeOHOlPeroHDg5JLFlShylTVNjYQEJCArt37+bSpUtFn03lPf/GczQaa9Rqhcmaexzlp6FnpVa4uOjnC44fP052djb+/v589NFHFBQUVHlZWBYsjFRhpOnp6YZY7nnz5tXo/XjgGamdnZ0hq0cPrVZrtlf65MmTBAYGGuym58+fZ8+ePeTn56PT6Qwe+enTp1O3bl127NhBbm7uXa+wX1hYiJSSd955hz59+hAZGVlpRlr1c2tZvvx7VqxYgZQ9sbb+GCEkUuYBPzN27EE6dmxNePgj2Np2Ys+ePSQnJ5Odnc2aNRuYPn0S+fk2HDnyGCqVIDfXjjFjdqBSafH1jcS49osQ0K3babp2PW0Yu3YNPvzQCrVai5WVEuGQkZHBoEGDqhQ+lZlph7NzXpFDyXxzgNKjSRT1Z1JO7OqqpbBQi0ajISsri2eeeYYZM2YYbNHZ2dm4urpW/iLLgYWR3vbaP/PMM/zyyy8cO3aMoUOH3utLK4Far0jL8trn5OSU8ETHxcVV2FUUlJJy+vjPQ4cOUVhYyO+//84zzzxDy5Ytyc/PJzo6mqtXr5KcnMzmzZtp3bo148ePr3bWURz6cKXvvvsOKeVd9drevAk//6zjl180XLki8PS8RL9++UaM/gcKC63IynLCxSUDlUrxLLVr156IiKYkJ/8EKNoxJkawdes5JkzYiJUVHDvWBx+fKADs7PLLVITG49bWYG19W95OnTrdUT1WZ+c8k/fp6S64upYe5lTymkqu1KytrcnJyeHWrVtMnz6dffv28fPPPwN3R5FavPa3vfb/+c9/DIp0yJAhdy1Soqq4b732NjY2Jp7W1NRUIiIiDEVKysOWLVsAxUaakZHBgQMH8Pf3p1WrVgghsLW1xcHBgcTEREJCQnj66acZOXLkXQtkN4a+YMh7771H7969eeedd+74mKmpcPo0RERAYSGcOQPDhsXj6anlrbcgL+8PJk36gsmTt+DmpmHbtlF8//1UAgOHERzcDZAm3TOFgObNo3j//Y+ZOPFXGjdWIh/OnWvLtWu+9Ox5AiEkfn4Rhv2rAn0Xyaqi+HmPHetFdLQ3+fnK55iba0dERDNOnuzGxYutivZSM2zY8+Tn3w61ysuzZdKkSbRv35727dsbKnktXbqU5kWtOEvz+t8pLF772/Lrbf7W1ta1TonCfcxIW7VqZcgft7OzIzAwkG7dulWo7PLy8gwhTzdv3jRpFaLRaAgPDyc4OJj8/Hx69+5Nt27dzGo9Uh24fPky3bp1Izg4mC+//BJQYkC/+eabKh3vf/+D2bPh0iXjakPKclWttmfo0K/o2jULGxvlQbRz55+cOnUKnU6Hp+dUevUKwsMjo0xFaGUFbdpcok2bS+h0EBXlg1ZrjatrBnZ2OXTocIaCAhVqta5KNs47ZaRw+5zK9XXlr796M3BgNrdu5RIdXRd95uyIETrS05ewf/9jDBzYhLCwLtSte5OGDRO4fLkFvr6taV0s1svZ2ZkNGzYwf/78u1L82sJITeNI3d3dmTJlyj28orJR651N5fVs2rdvH+fPn0ej0dC8eXOeeOKJCjOGQkND2b17N/369aNBgwb4+vqi0WhITk5m48aNNGzYkB49etCsWbMaU6Cg2HdLs4XZ2dlVqef2r7/C5MmKAmncGFq2TOby5RukptrTr981evYMxspKUVKFhdYcPFiP5s0lXl43sbK686Zx1dF4zthGeqfQ6SA724mFC/9u6DgKSpyrUs2+9HlC6GjYUBAeLqjpBrTLli17qG2kxXs2+fr60r9/f0MjxZrAA+NsKq9n08CBA/H19cXV1ZU6deqYRfk7deqElZWVIXDdzc2N3NxcbGxsGD58uCE27W4jKysLBwcHMjIyuHDhQomyXQMHDiQgIMCsWqjF8ccft0u5zZ6tITPzOdat+y9+fn68+upT2NnZFYX2KMrO2rqQIUPiq0MsA6pj9VUdjFQPlQqcnLKoXz+BxMTbURcxMbf3GT0atpt2RObwYRW9e1PpsKvqgIWRmjJSV1dXQ/W02oZKfT2EEC2AcGCjlHJq0dhk4BOgLrAXeE5KmVK0rSfwE6ADpkspTxaNRwH2QDMpZXbR2CxgqpTS3/ic5VXIV6lUlX46CSHo0KEDnp6eZGZmUlBQgJ+fn1ntRKoDycnJBAUFcfLkSUDJuS8oKMDLy4tFixYxe/ZsgxNt+fLlrF27tlJe+4wM0Be+at36MD/99DgtWrTgH/94BwcH05YltdDUZILTp09XGyMFRd5evYLYvt00PM7BAXJylKIrW7bAN98oxVbatYNSQhBrDBavvWmu/cyZM1myZEmVC4vcTVT2ObsUMNQwE0K0Bb4HngBOAcuBZYC+tO3HwMii/5cBxnEL1sDrwH/KO6E5XUSrgrp161Zr6qQ5OHPmDFu2bMHd3R1XV1fS09MpKCjgypUrBqeFMaZNm1ZpRjJ+PKSlSWxsjtG3789MmvRWrTTOm4PqZKR6dO4cxunTnbl+/XZxG50O3Nyga1eFub7+utK1s4Ls4bsOCyM1ZaS9evXiww8/vIdXVDbMNgIKIZ4G0oD9RsNTgB1Syv9JKbOA94FxQgi99rMq+rOmpNL+HJgnhHCjHFTURfR+wfz589m8ebMhEUBfH3TdunWlKlFQCjmX1dO9NHz1lVIUxMXlPG+9tQcvL6/7VonCnXvtS4MQ8PTTpi2n8/OVYszGi5JPP4W1a6v11JWGxWtvKn9oaGiNltWrDMxipEIIF+BfwCDgeaNNbQFDoUwpZYQQQgO0BEKAD4FAQAvMLHbYk8AhYB7wz7LOfbcYaU2iuLPok08+YeDAgXTp0qXcuNTKMNJTp+Bvf9Mya9b3NG6ciBA1Y6q4m7gbjBTA1jYPffQCKI4xf/9qP80dw8JITRnpsWPHDNXeahvMZaQfAT9IKWOKjTsBxftXpAPOAEVMtbmUspWUMqiU434AvCqEKLOw44PASO3t7fH09ASU1NS3336bHj16VBjcXxYj1Wg0/Otf/yIkJASA9euPMGhQFO+//2+aNEm6r1moMe4GIwXIyHBDr0T1qEJgxF2HhZGayr927doyC6bfa1SoSIUQnYDBwKJSNmcBxWvXuQBmaT8p5VlgJ1BmXbnCwkLy8vJYv379ff2ak5PD8OHDGTt2LOnp6Wa9jh49mrNnz5YYf/fdd/nwww/p168f9eo1YcWKGBo2HI5Gc+/lrM7X1q1bk5CQUO3HjYuzA8aiPPOV1+3bzf9caup13LhxnD59+p5fx72UPzQ0lPT0dJ588kkAjh49WmPfP+PmkBXBqqLOkQsWLHgaGA/MWbBgwTzAH2i7YMGC8UAQ0Gr+/PmbAYQQzYC5wFvz58/XlHFIFixYMBfYPn/+/GsLFiw4g+KkOgv4zp8/f7Xxvh9//PH8fv364eLigpub233xam9vz/bt2+nYsSMNGjTAzc2N4OBgXnzxRcaMGYOvry9eXl4Vvh4+fJikpCSeeuopw3jjxo2ZPXu24SHTsmUAgwal1gq5q/s1KiqK3NxcHnnkkWo7bkJCFw4depbCwmaAL+AF+PLtt160amXe51JTr/v37yctLY1x48bViuu5l/LfuHGDy5cvM2jQoBr5/jk6OrJ//35SU1Pj58+fv7wiRVphQL4QwgFT1jkPaAoEAPWLlKnea/89YC2lfJpyUBT+NEtKua/o/QpgHBBePPyprOpPtRGRkZH89NPtakkffPABKpWKtWvXcvXqVZKSksqMFJBSUlBQQG5uLo6OjlhbW5da/UffMRVgyJDJ9OnTstaHMVUVBQUF6HQ6bG1tK965DOhjZbVaFSdO9ODAgccoKDCttP/ZZ1BGzsc9haX6k1L9ydramqlTp5KdnU3Pnj1r5Nxnz54lOTmZgwcPVk9AvpQyB8jRvxdCZAF5UsokIEkI8RLwX6AOsA8ovbBn+fgXMK20DfeLjfTkyZPs3LkTAD8/PyIiIli8eDH169fn6tWrPPfcc6UqUY1GQ1hYGCdOnODWrVuAEng8fvx4du3axbVr1+jZsydnzpyhZcuWXLhwgffem2/IPnpQlShUTxypvkfTt9/O5tat23Ur27eHhARITFRCxmojLHGkShypg4MDGzduZN68eTV27rS0NNq3b2/S8rw81PoU0fuFkX755ZeGLp4ffvghW7Zs4cyZM4CS2lY8cwmUp97u3bvJzs7G09OT1q1bY2trS1BQEGlpabi6upKcnHzXK07VVlQXIz12rBf79w/l0UcFf/4JGg0EBMCiRXDuHHTpUo0XXY2wMFKFkX7yySf88ccfuLm5ldvZtzpx8+ZNtm3bRnx8vFmM9L6t/lRbkJ+fj1arNak6lZOTQ6dOnQDo0qUL58+fLzHv6NGjbNq0CUdHR2bOnMmsWbNo0qQJ+fn5+Pr6otNJ9u/fz6FDh2pKlFqH6vDax8Q0Ye/eYQweLDhwAPbuheeegxdeAFvb2qtEweK118u/b98+goOD2bt3b42du2HDhnTu3Nns/Ws91antcaT69E3jroU7duxgUlGeZlRUVImKVCdOnGD//v20b9/e0FcpMHAXISGGpDFUKkHHjh3vShxlTUK/4Km8CULNyJEjcXJyon79+jg7O6PT6QgNDS3hTS1eIEVKCA9vTYcOF/H2voEQt/jlFwfAgQEDlLbQ9wMscaRKHKmrqyt//fUXoCjXstoCVScyMjLM6varh4WRVhO2bt1q+P/ixYvk5+cDSt1UY4SHh7Nr1y5atmzJwIEDiYpK4N///q+JEtUjLCzM0E30XkKvDDUaxclVnm7Xbzt6tDsff3yDBQtWsmDBP0lJcTU5Vnmwt3fkpZdmAXDhwgWGDh1K79696du3L6+88gpz5swhOzvbUPhap9MaHjiFhVrOn2/C5s0TDcfz9HyZ69cvV1Lqew8LI1Xk//vf/25ov1wZ5VZV5OTksHjxYiZMmGD2HIuN9A5RWFjIqlWrSrR1bty4MbGxseTn5xuUqU6nY9GiRbi6ulbYBhqqx0ZYHTD+ipTHLAsKlHYee/akcPLkeGAAMIJmzdYTE+NEnTqZvP++Izk5S8nJyS4z3W/y5Mm0aNGiQhvh999/z8svv4ynpyfDhg3j6tWrXLt2DWdnZ86fP0+9eqMZObIVGzd+xx9//GHoOHu/wGIjvd1FFMDHxwd/f/9q6dRbHg4cOICzszM7d+5ECPFglNGr7YzU2tqaFi1aGBTj6NGj2b59O7GxsYSEhJgw0suXL5OVlUW/fkOIjd1S4bHDwsJIS0tj8ODBd+36zYHe861/zc0VODhILl/WkZl5ni5d2pGXp2bPngOcPv2/ollrgD9xdNzJxIk25OXB++/XwcMDcnICUKvV/Pnnn1y/fp2cnBzS09NRq9VIKfH29gaULqIxMTFl1hsYMWIEOp2OGzduMGnSJAYPHmzI6oqJiWH16tV88MEHAPelw87itTet/mTcS+1uobCwkPDwcDZs2FCpeRZGWg3Izc0tsQQbO3asoeapHr/88gvR0XF8/fUb3LypYsaM1Xh5xVBQoMbePr/Ece8WI61K0eW0NBeCg/dy4sQ+3N3dGTVqFOvXryc7Oxt39+6kph5AKQw2BKW0Ajz66AzOn1+Nmxv89Re4u5d9fH0jQysrK0NBbXMY2ffff8+QIUPKZCk6nY6zZ8/Svn37+y511sJITRnprFmziIuLuyuxpBqNhgMHDnD9+nXi4uKIiorCx8fHbEZqsZFWA65cuVJibMmSJSbvMzIyuHz5CgcPduLmTRU2Nvn4+kaTmFi/VCUK1WsjNXd5XhoyMhz5+uu5/PnnX0gpSUxM5IcffiA7OxuA1NRgYA7wGL6+awzzNm78gs2bISrqdo3UsqBWq1Gr1SZdCX766acKa7HOnj273KWeSqWiQ4cO950SBYuN1Fh+KSUbN27E29vbYBIqKCggLCyMr7/+mtOnT5d3qHJRUFDAunXr8PLyYuHChVy4cAEfH5+KJxqh1q93arvXHpSYM2MMGzbMUKREj6Cgs4AkNFQJqdBobEhJccHT82aZDLE6vfZ3okeuX/dGSsG4caPYsOFH5s2bx6JFpqUX3nyzAVFRGWzYMBnIBv7B9evX6devLhMnwoYNSuO9yqywZ8yY8dB7rR92+fVKMzIykvT0dNasWUP9+vWJjo6mYcOGht/e1q1bcXd3r7QCBLh06RI+Pj6sX7++yu2FLIy0GjB06FCTXtstWrQosc/WrQrrTE52RYgkmjW7jJtb+bJVByNVerSb5y0vC6GhSkXyzZtXERio4osvFrJ3714TJvj555+zYUM9hFiIktx2kV9/1ZCbC4GBMG5c5dt1mMNIH2RYGOlt+Rs0aMD69eu5ceMGa9eupWHDhqxevRqtVouUkm3btrF69epSbajJycmEhIRw9OhRwsLCuHDhAsePH2fPnj38+OOPhIWFMXv27Dvq0WaxkVYj8vPzDT98/X09duwYdnYerFx5EJWqkKVLXwEEr732NQ4OWdjZFZZ5vMrYSKWE7Gw7wsPbEx7eFi+vOHx9o6lfPxFb23wuXGjN8eO9aNDgJq1aXaJJkxu4uaVT0XdHSvj00zfJy3MwjE2aBOvXK/GyTz75JC+//DJCCIKDgwkPD8fevjdSLiclpRlt2yrZQwcPVr7m58NuI3zY5S9uIy0PUkpUKhV16tTh1VdfBZTOuNHR0ezbt49x48ZRt25dYmNjycrKonHjxqSlpTFp0iTi4uJ45ZVXSj2uxWt/D2Bsp7lx4wZ2dnaGbIwGDWD79uGAwN4+G0fHNIqFmJaAuV77zEwHVqx4kYwMV8NYXJwPf/3VGwAnJzh6FCZPzuLcufacO9ceyMfaWtCjRzB9+hzFySmv1GMLAV5eMURFJVFQ0A9QOpS6usJ3340kLi7OpNiulBIpJRqNipkzlX0BHn20fFlLQ0Ve+wcdFq+9qde+PAghmD59OmvWrCEjI4OsrCx+/fVXhg0bxpo1a0wSZu4GLIy0mlC8HGFAQACPPvooFy9eREprhCjk66/nMGXKSurWLd25VBzFGWlpGTyHDg3g8GF/ihcqBqV1Rt26SnEOOzvw8dnDpUsrUas/oaDgKmp1C6ysGpCXd5SnnsqiXbtzJsfWnysioieffdaFI0fsmTjxts36zTeVykllITdXUaBz5sCMGWaJXGz+w83IHnb5K8NIQXmIz549m82bN5OcnMzTTz/NL7/8UqVzZ2dnY29vj5WVlcVrX5OoU6zp+bfffmuILc3NVQPw2mtLzVaiUNJGqo/j1OO33yZy+PBjDBwoaNy45HwpwcND+T8vDy5dGsaIEWqaNRuJSvUEBQXHycvLAp5k06Zx7NgxgtxcOzQaazZtGkNqqhsJCQ0YPNibhg3rM2GCM8bk6PPPoTyyZG8PJ05UTYmCxUZqsZFWTn4hBNOmTTPY7hs0aFDBjJI4d+4cb7zxBn5+fqX6Oso8t4WR3jnOnj3Lpk2bTDJ13nlnKfb2yXfkddczUhsbW4SA8PC2XLzYmgkTNrF//2McOVJ+0rivL0RGKv9PmKB4zvW2So1Gw5df3uKddzz56KPf2bZtOCdPCoTQYWWlpbBQzT/+oTSBK44//oBRo5QqSgCvvgqLF1dZzDLxsDOyh13+yjJSgKlTp3LmzBm8vb1JSUlh9+7duLgUb+IB165dKxE2d+LECZ544gmaNWuGt7c3zs7OrFq1ysJIawpXrlwxUaJTp07HxibpjkOXQkPDOXjwGAUFatavn0hg4FP88ENbRo4cxxdfdGPVqrIdOLNnw++/g4+PUnvz5ZeVcf0l2djY8Nxz+hCtEfz1lyA4GFq2VFFYqMbHp3QlCjB0KGRmgj4u+ptvICzsjkQtFRZGamGklZFfSsm+ffsYOHAgnTt3pqCggK5du7Jjxw62bdvGrFmzePnll/n444/x8/Nj5cqVJCQkIKUkMDCQkSNH0rVrV0aMGEG7du0qFUpV651N90McqfE1jhgxgubNq5YLbFwpqbBQ0Llze2JjG7By5UQSExswaBD07CmA9gD06AHTp99uIywEjBmjhBmFh0O3buDoqDh80tKUfYwLUdWvD927Kwr3pZfg7bfh0iWYOROWV9BcwcYGjh+HL76AJUsU51N1wxJHaokjrcyKOTAwEAD3ohS6kSNHsF4nTgAADb9JREFUEh4ezuuvv46Ukvr162NlZUVQkNKH84UXXsDW1tZQYWrgwIG0bt26Stda6xVpbWekCQkJJjUz7ySDRmmJIUhN9cDBwQdb23AiImJJTFRsPaV5+fVLd9A7nyArC/r2hWeegfffB29v2L5d2ad4JNWIEYqds1s3uHkTVq1SFKm5mDdP+bsbsHjtLV57c732UkoCAgLo37+/yXj79u1p3759if2ffPJJCgsLsba2Jjs7GwcHhzuKI631irQ2M9KwsBS2bPnW8N7f35/u3btX6Vg6Haxf/zSXL7cEBL17w759g/nHP7SsW6fYJUuL4PjtN9P3DRsqKZn//CcMHKgo55QUZfkNUMwnxogRsGCBknV05IjCUGsLLIzUwkjNZaSHDx8mJyeHNm3amLW/EAK1WnECV4cNutYr0trKSNesiSIubj2PP/44PXv2rBQTLS0lVKWCmzfXAQsACAqCZct+Jjk5kk8++YQXXyz9WM8+Cx07wn//C+vWKWxy8WIYPBi8vKBzZzh7FmJilOV606am87t3hy1boE8fZalfm2BhpBZGai4jff755+nXr98dsco7Qa1XpLWJkep0kJOTxY0bN7h2bT39+/evdI3L3FwlLCg8vB1+fhE4OOSSkJDAsGHDyMh4B2hHr14OLF78BO3aTauQkTRsqLBKf3+FgU6fDq+9BmvWKB760FDFZnrokKIsi0NvV62NsDBSCyOtiJFKKVm+fDm3bt2iQ4cONXRlJVHrFWltYaS5udacOnWYI0eO0LhxY9q3b1+l7pZqtYrExEI2bTqKtfUveHrasnLlSoYMGcK1ax1o1qwZ3bu/SvfuT7B8+c9ERkaa5bl2cIDnn9efA158kTJZ7P0CCyO1MNLyGOm1a9cICAggPDycKVOm3NOas5Y4UjNx8uRJjh49yrPPPotrFVzUBQVWqNVaYmIas25dNrm5cwFYtmwZAQEBJvtKKRFCPPRxhBb5H275y4ojPX78OF9++SU7d+6kT58+9OnTByt96Eo1Y/78+ZY40upCVFQ6O3fupFevXpVWosbPKY3GGg+PFITogre3NxcvXiyhROG25//nn39+aNkYWOJILXGkpvJLKZk0aRLDhw8nMzOTuXPn0r9//7umRCuDWr+0v5c2Up1OR1jYVbZtWwdg0nLZHGRkOOLklI0Q8NtvE5g06TdiYxuSk/M2Tz75JK1atSp3/rRpFdtIH2RYbKQWG6nxivnzzz/nyJEjvPLKK7WudYyFkZaBW7du8fvvu9i2bR2jR4/mn//8J/b29mbPlxJcXLJRqSA/34bISD+OH+9Hy5aRzJgxnU2bfmfLlvL7NlkYqYWRWhipIv9//vMfPv30UyZNmlTrlChYbKQlUFhYyL///W/D+xEjRtCjR48y9y/et11fRNk4CuPcueZs2DAFIXT07XuMjh1D2LjxU7y9m3DixIkyj/2w28gs8j/c8uttpImJibRt25aAgIBS8+bvJsy1kdY+1V4MNc1IjStsv/vuuyX60hdH8YpMQpSMET1/XumpLuWfHD3an6NHQ7GxSebo0f9RHn7+2Xyv/YMIi9fe4rVPSEjgmWeewcbGpsaVaGVQ6xVpTdpIY2NjWbFiBQDvvfeeIfOhOPSKUykAIlCpymb1UkoiIgYCocBgIB8Hhxx27txN8+bNy70ei43UYiN92OWXUjJr1izatm17ry+nXFhspEbQK9GXXnqpTCV65Yod/fsPws/vbUaNGsvEiXPLbSwnhCAvTwPk4uGh5AF//vlbPPbYYxVej8VGarGRPuw20lmzZrF79278K9unpoZhFiMVQqwFBgGOwE3gMynlyqJtg4ClgDfwFzBTShldtK058CvgAbwspdxVNH4I6AW0kFLGFI0NBlZKKZsan7umGKmxrdi4dYYxTp3SsXnzPKNwiw4lOoiWBnv7W+Tn98HWdhcNGsDo0eZdk4WRWhjpwy7/+++/z8iRIys0sd1rmMtIPwGaSildgNHAv4UQXYUQdYHNwPsoyvIkiuLUYwHwGtAT+LDYMbOL5pWLmmKk165dAxTlVRxSQkKCmm3bFpSIWUtMTDT87+7ujqOjI+fOnWPJkiXs2bMHnU7H8OGPodM5kpWlZtcuaNLEvGuyMFILI33YGWlmZiaNGjW615dSIcxipFLKc8Zvi/78gK7AOSnlBgAhxHzglhCitZTyImBVdA4roPhaeTEwTwjxmZTyalnnLo2RSim5dOlSlWsHlgY9s/T29i6xrV49b+bPf7bUeR06dKBVq1bldvoMD1fy3V9+WSkiYi4sjNTCSB92+d99913y8vIqVSX/XsBsG6kQYpkQIge4CMQDvwNtAUNtdCllNhBRNA7wLxSF+VfR/8aIBVYA88s7rzEjzc7OpqCggIsXL7J+/XoyMzM5deoU8+fPJzY2lgMHDhAREWGuSCbw9fUtMSalZMyYMcyZU7oS1aOidsnt2yvFlSvbSdPCSC2M9GFmpCtXriQjI4O8vNI73NYmVCqOVAhhBfQG/IFPge+AJCnl20b7/AmskFKuLuc4h4C1wDbgKtAHaEQpNlInJyfp7u7OjRs3zL5OgP79+5sUFdFqteWmkhUUFPDpp5+SlZVlVh/5msDDHkdokf/hlj8nJ4elS5eybNkyZlam2ng14q7k2ksptVLKo0ATIADIAooHd7kAZhk2pZRJwBJKslUDsrOzSU9Pp0OHDkRFRfHiiy8SHByMg4MDoDDJs2fPMm/ePK5fv86wYcMAOHLkCLGxsezdu5crV67w0UcfGeyW+fn57N271+T10KFD1KtXjzlz5pCVlcVbb711z19XrFjB4MGD7/l13KvX5cuXM2TIkHt+Hffq9bvvvmPo0KH3/Drupfw//vgjDg4OJX6vNfVqLqqU2SSEWIniLDoHzJBS9i0adwSSgC5FNtKy5h8C1kopVwoh3IBrKAw3oDgjFUIkAdGVvkgLLLDAgjuHj5SyXkU7VehsEkLUBwYCO4FclKjyZ4DJwDHgcyHEeCAQ+AA4U54SLQ4pZZoQYiHwD0phsuYIYYEFFlhwL2HO0l6iLONvAKnAF8BcKeW2oqX5eODjom09gaercB1fAw+ve9ICCyy4r1Hri5ZYYIEFFtR21FiKqBDiESHEASFEuhDiqhBirNG2QUKIi0KIHCHEQSGEj9G25kKIECFEpBBieNHYM0KI88WOv7eMsbe5xxBCZBX70wohvjHa/kDLD4pdXAiRZ3QPLhltmyyEiBZCZAshtgohPIy29Sy6N+eFEN2Kxt4RQvxe7PhXyhirygqp2iGEWCuEiBdCZAghLgshZhlte6A/fyHEK0KIk0KIfCHE6mLbHgzZpZR3/Q/FFnsZ+BtKcP5AFGdVS6AukA5MAOyAz4HjRnP/C/QF6uvHAU8Uk0M9o+MnAlHFxrKAPjUhYyXuhWPRdQ0oev9QyA8cAmaVMt4WxTY+AHAC1gHrjbbvA5oX/f1RNNa36J5ZFb1vWCT7zWJjEvC817IbyWlb9H/romvt+jB8/sA4YAzwLbDaaPyBkb2mGGnrohuwSCohVAeAP4FpKDf5nJRyg5QyDyVAv6MQQp+2VCI7SkoZh+LpH1C0TxeUCILDxcZUKGmrtQlPoXzwR4reP2zyF8cUYIeU8n9SyiyUtOFxQgh9SpsVt++B3jkajHIvOhW9HwAcBC4VG4soulf3HFLKc1JKfTyNcXbgA//5Syk3Sym3AsnFNj0wsteUIi2tPpIA2lH17Kj/cfvGDUBRTEeLjR2XUmqoXZgBrJFFj04eLvk/EULcEkL8KYTwLxorLn8EoEFZrYBSoyEQJXnjvaJ9NCj3oyL5yy/4WsOo5uzA+/HzL44HRvaaUqQXUVjYm0IItRBiKPAo4ICynEsvtn864AwgpTwvpewopfSWUm4z2sf4CdQf5WYeKTZ2+G4IU1UIIbxR5P7JaPhhkf8toBnQGFgO7BBC+FGx/P+TUjaXUraSUgYZ7XO/yY+U8mUUufqjFPvJ5+H5/EvDAyN7jShSKWUBio3kCRTb0N+B31BCqqqaHfU/oIMQwh2lJF+QVOJXGxWN9aOWMRJgOnBUShlpNPZQyC+l/EtKmSmlzJdS/oRi2hnBncnfr0jWelLKKyhxzX2KxtpRi+TXQ1ZfduB99fmXgQdG9hrz2kspz0gpH5VS1pFSDkNhJydQ7BuG9pxCyY7yKxov73jXgDjgReB6kX0NIKhozAk4Xu2C3BmmY8pG4eGS3xgSxbxTXP5mgC2Kc7I8BAGuKLL+CSClzOD2PYkr9sCqbbDm9uf8MH7+8CDJXlNeLaADimfOAZgHRKL8YOqh0PnxRds/xchzV8Ex1wEJwJdGY38vGvuzpmQz81r7oEQqOBcbf+DlB9yAYUXyWaM4mLKBVij2sAyU5ZgjSjGb9WYe91iRrK8ZjX1TNPbfey230TXVR0lUcUJxnAwrkv/Jh+Tzty6S7RPgZ6PvwQMje03ezM9Rsp+ygF1Ac6Ntg1HsqLkoYTJNzTzmbBRmM85orEfR2Cf3+gtU7Fq/B34uY9sDLX/RDyYYZcmWhsIWhhhtnwxcL1Iu2wAPM4/7SZGsXYzGJhaNzb7XcheT/3CR7BlAOPDCQ/T5z+d2pIL+b/6DJLsls8kCCyyw4A5R65vfWWCBBRbUdlgUqQUW/H87dSwAAAAAMMjfehB7CyKYRAowiRRgEinAJFKASaQAk0gBJpECTAHJmSyFD90hywAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x2aaafc6ffb70>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Setup basemap\n",
    "def setm(ymin,ymax,dy,xmin,xmax,dx,res,ax):\n",
    "        m = Basemap(projection='merc',llcrnrlat=ymin,urcrnrlat=ymax,llcrnrlon=xmin,urcrnrlon=xmax,\\\n",
    "                    resolution=res,ax=ax,suppress_ticks=True)\n",
    "        m.drawcoastlines()\n",
    "        parallels = np.arange(ymin,ymax+dy,dy)\n",
    "        meridians = np.arange(xmin,xmax+dx,dx)\n",
    "        m.drawparallels(parallels,labels=[1,0,0,0],fontsize=12)\n",
    "        m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=12)\n",
    "        m.drawlsmask(land_color='grey',ocean_color='white',lakes=False)\n",
    "        return m\n",
    "\n",
    "def plottraj(year,d,basemap,color):\n",
    "    file=path+\"DATA/GS-\"+str(year)+\"/ariane_trajectories_qualitative_GB_initt.nc\"\n",
    "    data=xr.open_dataset(file)\n",
    "    traj_lon=data['traj_lon'][:,::d]\n",
    "    traj_lat=data['traj_lat'][:,::d]\n",
    "    xx, yy = m1(traj_lon.values, traj_lat.values)\n",
    "    for i in np.arange(0,xx.shape[1]):\n",
    "        m1.plot(xx[:,i],yy[:,i],color=color,zorder=1)    \n",
    "\n",
    "fig,ax=plt.subplots()\n",
    "m1=setm(30,70,10,-90,0,20,'l',ax)\n",
    "artist=[]\n",
    "years=[2005,2006,2007]\n",
    "colors=[\"k\",\"blue\",\"gray\"]\n",
    "\n",
    "for i in range(len(years)):\n",
    "    plottraj(years[i],d,m1,colors[i])\n",
    "    artist+=[mlines.Line2D([], [], color=colors[i], label=str(years[i]))]\n",
    "plt.legend(handles=artist)\n",
    "plt.savefig(out,facecolor='w',dpi=300,format='png',bbox_inches='tight')\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "SSH shkpwagn@hdata2.hlrn.de python_py3_std_/gfs2/work/shkpwagn/NB_WDIR",
   "language": "",
   "name": "rik_ssh_shkpwagn_hdata2_hlrn_de_python_py3_std_gfs2workshkpwagnnb_wdir"
  },
  "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.5.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}