Skip to content

Instantly share code, notes, and snippets.

@monodera
Created May 2, 2015 17:53
Show Gist options
  • Save monodera/0ee9e8600b2b802c62c3 to your computer and use it in GitHub Desktop.
Save monodera/0ee9e8600b2b802c62c3 to your computer and use it in GitHub Desktop.
Test Simulate Broadening
{
"cells": [
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline\n",
"\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def gaussian(x, a, x0, s):\n",
" g = a/np.sqrt(2.*np.pi)/s * np.exp(-0.5*((x-x0)/s)**2.)\n",
" return(g)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"sigma_intrinsic = 5. # [km/s]\n",
"sigma_total = 200. # [km/s]\n",
"\n",
"sigma_total_rand = sigma_total * np.random.randn(100000)\n",
"\n",
"x = np.linspace(-3*sigma_total, 3*sigma_total, 1000.)\n",
"\n",
"gauss_intrinsic = np.zeros((sigma_total_rand.size, x.size), dtype=np.float)\n",
"\n",
"for i in range(sigma_total_rand.size):\n",
" gauss_intrinsic[i,:] = gaussian(x, 1., sigma_total_rand[i], sigma_intrinsic)\n",
" \n",
"\n",
"gauss_total = np.nansum(gauss_intrinsic, axis=0)"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAtcAAAEMCAYAAAAPn5osAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8FEXeP/DP3DOZIzcJEJJwBAKCCHIoingg7goeqBg8\nVllgV0V0cUVFH1mCgqCrj7cuHqi7j3jgtcoqiiiwiIIuyH0TCCSQzH3PdM9M/f7Ir3sTksAMSWbS\n09/365UXZK6u6pl8prq6ukrBGGMghBBCCCGEtJky1QUghBBCCCEkXVDjmhBCCCGEkHZCjWtCCCGE\nEELaCTWuCSGEEEIIaSfqVBdAznieh8/nQzAYhNfrhd/vRyAQgNPphNvthsfjgdPphMfjQTAYRDAY\nBMdxCIVCCIfD4DgOPM8jGo0iFosBAJRKJTQaDfR6PYxGI3Q6HbRaLSwWCywWCzIyMmA0GpGZmSn+\nZGVlwWg0ivdlZGRAoVCkeO+0L8aYuD/9fj98Pp+4z/1+P4LBIDweD1wul3ifz+cDx3GIRCKIRqOI\nRqNN/i8Q9pVCoYBKpYJGo4FarYZarYZWq0VGRgYMBgNMJhNMJhMsFgsyMzNhNpuRl5eHrKwsZGVl\nIT8/HwaDQfL7njEGh8MBu90Oj8cDj8cDt9sNl8sFm80Gl8sFr9eLQCAgfp6Fz3YkEhE/y8K11kql\nEiqVCkqlUtynGo2myY/BYIDZbEZ2djZyc3PFz3F2djZycnJgMplgNBphNpuh0WhSuXuIhFFmJw9l\ndvJQZrc/hVRnC/nkk0/w7LPP4ssvv8SSJUuwf/9+GI1G5ObmoqCgADk5ObBYLNDr9cjMzESXLl2Q\nnZ0NrVbbruVgjIlBK3wgbTYbTpw4AYfDAa/XC5fLhbq6OtTV1cHn88Fms8FqtcLr9ca1Da1WC4PB\nAIPBAK1WC71eD71eD61WC7VaDZVKBZVKBcYYYrEYeJ5HKBRCIBBAOBxGOByGx+NBKBSKu15GoxE5\nOTnIy8sTP/x5eXno1q0bcnJykJWVJf6r1+thNpthMpmg1+thMpmgUqnOdJe2iOM4uFwu2O12uN1u\nBAIB+P1+2Gw2uN1u+Hw+WK1WMWh9Ph/sdjvsdrv4ZRfvRz0jI0Osz8n7uPH/GweqsO+j0Sh4nkck\nEkEkEgHHcQgEAgiFQvD5fAiHw6fctsFgQF5eHvLy8pCfny8GuMlkQnZ2tvg5zs7OhsViEYO/vT/b\njDHxs+P1emG321FfX4/a2loxgIV9XF9fD5vNJn4J2u32U9ZTpVLBZDIhIyND/Gw3/jwrlUoolf89\nqSZ8McZiMXGfCvuY53lwHIdgMAifz9fkC7Q1mZmZyM/PFz/bWVlZsFgsMBqNsFgsyM/PR25urhj8\nOTk5MBqNMBqNMJlM7RL0jfPLbDa3+fWkgDKbMpsymzJbqpkNJJbbku25ttvtWL9+PdxuN44fP47P\nPvtMPLqtrKzEsWPHAADz589v8jy1Wi1+SPR6PXQ6HTQaTbM/SIVCIf7xCR8Q4UMh9FaEw+G4Phwm\nkwldunRBYWEhsrKy0Lt3b/H/FotFPEIWPgxZWVnIzMwU/wB1Ol277DOe5xEIBJCZmYnKykrx9v79\n+4s9AT6fT9yPwh+f3W7H9u3bUVdXB6fTGde2hD88rVYLnU4HvV4vHrEK9wlhzhgDY0z8gxSOioPB\nYJN9fTpms1k84jWZTMjPz8eAAQPE/ZmZmSn+0ZnNZpjNZrHnRziKzszMhEqlarJ/Gv+/rSKRiPiF\n7vF4YLVa4Xa7xf1stVrF8LPb7aiqqoLdbofX6wXP86d8bSH0TCYTDAYDVCqVeJtGo4FSqRS/XIQv\nFeFzLbz/HMeB4zj4fD5EIpFWt8UYQ2VlJUwmE0pKSrBhwwb07NlTbCx1794deXl5MJvNYq9PVlYW\n8vLyYDKZOqSnR+jpcjgc4ufG4XDA6XTC5/MhEAiIDan6+no4HA5YrVbs379f/OIJBAKn3U5GRoa4\nj4XPt9DbKGTHunXrmnxuVq1aJX7ZzJ8/v0l+yaVx3Vpm+/1+8fMkaJzblNmU2fFkNoAOyW3KbMrs\naDSKxx57LKHclmzP9YcffoiKigrs2LEDZ511VpP7Gn8IvvnmGzidzmZHyMKRaTgcFgNYOH0knOIQ\nXks4tSEEunD6QqfTiUehFotF/FDm5eWhoKAA+fn5MBqNUKs73zHMyUfw8QqHw3C5XOLpIrfbLR7h\ne71e8f9+v188DRoOhxEKhcDzvNhDw3GcuJ+Fsgi9DDqdTuz1Efa1xWJBVlaW2AMg/LEIvws9Fqne\nPx1J6I0QwkcIF6/X2+RUtBC6Qo9BMBgEz/PNPtfC/tZoNOKXlfDlKnyuhc94Tk4OunTpgm7duiE/\nPx+ZmZmdch+1Fc/zYu+Z1+uFw+EQc0NowAi/C19sQm8Rx3GIRqNgjGHNmjVN9s+ll14qnia95557\n4Pf7W82vdNVaZgcCgWbDGj744APK7JNQZnfcPuoolNkdLxmZrVKpcM8998Dj8cSd250vQeKUnZ0N\noKE35FQuv/zyZBRHNnQ6HQoKClBQUIB+/fqlujiyIvTclJaWprooaUuj0aCwsBCFhYXt+rqrV69u\n8vuqVasAnD6/0klrmZ2RkdHssTfeeGNSyiQHlNmpQ5nd8ZKV2UBiuS3Z2UK6dOkCAKivr09xSQgh\nJDFyzC851pkQkj4SyTDJNq5zc3MBAA6HI8UlIYSQxMgxv+RYZ0JI+kgkwyTbuBYGk/t8vhSXhBBC\nEiPH/JJjnQkh6SORDJNs41q4GjueK5IJIaQzkWN+ybHOhJD0kUiGSbpxrVar4fF4Ul0UQghJiBzz\nS451JoSkj0QyTLKzhSgUCphMpha75+fNm5eCEkkL7aNTo/1zerSPTu1U++dU+ZWuTldn+jydGu2f\n06N9dGq0f07tdPsnkdyW7DzXAFBQUIBrr70WS5YsSXVRCCEkIXLMLznWmRCSPuLNMMkOCwEaVjfi\nOC7VxSCEkITJMb/kWGdCSPqIN8Mk3bg2GAwIBoOpLgYhhCRMjvklxzoTQtJHvBkm6ca10WiE3+9P\ndTEIISRhcswvOdaZEJI+4s0wSTeuNRoNeJ5PdTEIISRhcswvOdaZEJI+4s2wpDSuXS4XduzYgWg0\n2upj/H5/i0cD1dXVOHz4cIvPUavViEQi7VVMQghJGjnmlxzrTAhJH/FmWIc3rpctW4aBAwfid7/7\nHYYNG4Zjx441uZ/nedx7773o3bs3SkpK8PDDDwMAYrEY7rzzTlx44YUYO3Ysfv/73yMWizV5rkql\nOmWDnRBCOis55pcc60wISR/xZliHNq7tdjseeOABrFu3Dlu2bMEtt9yC+fPnN3nMhg0boFKpcOTI\nERw8eBDPP/88PB4Pli9fjr179+LgwYPYu3cvjhw5ghUrVjQtvFIJCc8kSAiRMTnmlxzrTAhJH/Fm\nWIc2rlesWIGrr74avXr1AgBMnjwZq1evbvKYMWPG4Nlnn4VOp8OOHTtgNpthNBrx4Ycf4oEHHoBG\no4FKpcKkSZOaPZcQQgghhJDOpENXaKypqUFJSYn4e3Z2Ntxud7PHMcbw9ttv45FHHsE777wDlUrV\n7LlZWVnNlpxcu3ZtxxWeEEI6EGMMCoUi1cVIKjnWmRCSPuLNsA7tuc7MzITT6RR/t9vtKCgoaPIY\nxhjuuOMOvPXWW9iwYQN+85vfiM91uVxNntulS5cmz62srIRCoWj1p7KysuMqRwghp3C6fDp8+DCU\nSklP2JSwWCwmuzoTQtJHvBnWoSl37rnnYt26deL4lG+//Rbnnntuk8f8+OOP+PHHH7F69Wr07Nmz\nyXMb90yvXr262XMJIUSqotEoNBpNqouRVDzPy67OhJD0EW+GdeiwkBEjRkChUGD69OkYMGAAFi5c\niK+++gpAQ0N7zJgx2LNnDxQKBf7whz/A6/UiNzcXzzzzDKZNm4YLLrgASqUS1dXV2Lp1K957772O\nLC4hhCRNJBKBwWBIdTGSKhgMomvXrqkuBiGEnJF4M6xDG9dKpRLfffcd/va3v6G2thbr1q3DwIED\nAQAzZ87EqlWrMGHCBAAN47FNJhPsdjs0Gg169+6NTZs24dVXX0VxcTH+85//QK/XN3n9999/H5Mm\nTcKHH37YkdUghJCEVVZWnnJoWnl5ObRabfIK1AlwHCe7OhNC0ke8GdahjWsA0Ov1mDVrVrPb9+zZ\nI/5/6tSpLT63uLgYixYtavW1o9Eo1OoOrwIhhLQ7OeaXHOtMCEkf8WaYpK8sCQaDsjutSghJD3LM\nLznWmRCSPuLNMEk3rn0+H4xGY6qLQQghCZNjfsmxzoSQ9BFvhkm6cU29IIQQqZJjfsmxzoSQ9JH2\nPdccx4HjOFgsllQXhRBCEiLH/JJjnQkh6SORDJNs41pYnCY7OzvFJSGEkMTIMb/kWGdCSPpIJMMk\n27gWlkI3m80pLgkhhCRGjvklxzoTQtJHIhkm2ca1z+cDAJhMphSXhBBCEiPH/JJjnQkh6SORDJNs\n41o4gqDxe4QQqZFSfkUiEezatQtut7vJ7YFAAKtXr8bRo0eb3M7zPL7//nscOHCgye1SqjORLp7n\n8euvv+Lnn38WG0OEtIdEMkyyjWsav0cIkSqp5NfmzZsxYMAA/P73v0efPn3wxRdfAAC2bduGwYMH\n48knn8SYMWPw+uuvAwCqq6sxZMgQPP7447jyyiuxYMEC8bWkUmciPYwxMMYAAD///DN27tyJPXv2\nYN26deLthLRVIhkm2aWyhErm5OSkZPsOhwPvvfce/H4/tm/fjvr6ehgMBpSWluLGG2/EqFGjzvi1\na2pq8Oijj+K1116DRqNpx1ITQjqDVOdXvFatWoU333wTo0ePxrJly7BkyRJcddVVuO+++7BgwQJU\nVFTgxIkTGDJkCG666SY8+uijmDZtGu677z74fD70798fU6ZMQVFRERQKBYDU1ZkxBqvViry8PCiV\nye1X8nq9mDJlCl5//fVO/55LzY4dO7Bz504YjUaMGDEC1dXVsNvtCIVCUKlUcDgcyM3NTXUxSRpI\nJLcl23NdX18PACn7o1m9ejVmzpyJv/zlL9i7dy/0ej18Ph8++OAD/M///E+bXru2thZvv/02Dh06\n1E6lTX+bNm1CJBJJdTEIiUuq8yteDz30EEaPHo1YLIYVK1Zg0KBB8Pl82LVrF2688UYAQGFhIcrK\nyrB9+3Z8+eWXuOuuuwA0jEu88MILsXHjxiavmao6f/TRRygoKIBOp0NeXh5yc3ORmZmJnJwclJaW\n4uyzz8bDDz8c12tt374dXq837m0fPXoUn3zySbMhNKRtPB4Ptm3bBqfTCbvdjm+//RaxWAxut1sc\nEmKz2VJcSpIuEsltSfdcazSalF153qVLFwANIVtWVibezhhDNBpt02sLpxyCwWCbXkdOLrjgAvz9\n73/HTTfdlOqiEHJaqc6vRBw/fhxTpkwBACxZsgR1dXUoLCwUe6KBhjGIdrsdSqUSer1evN1sNiMQ\nCAAAwuFwSk/RX3XVVfjoo49w5MgRBAIBPP300xg+fDjGjh0Lj8cDl8uFfv36xfVaFRUVuOmmmzB3\n7ty4Hi8c+KtUqjMuP2lu7969ABo6pFQqFfr06QPGGILBIBhjiEQiza4VIORMPfTQQ5g7d25cuS3Z\nxrXP50vpF5PwxXJyWCoUCqjVie3WWCzW5DSl0KjOyMhIuFyxWAwKhaLJF9+pHD9+HD6fr8kBQmtq\namqg1WqRn5+fcLlqamqQk5MDg8EAnucBIO4hL/E8V6fTwW63J1wuQlIh1fkVrwMHDmDs2LG4//77\ncffdd0OpVILn+Wa9gXV1dSgtLUUoFEI0GhVzsa6uDgUFBQAAq9WKyspKzJ8/v9XtzZs3D5WVlR1S\nF71ej+uvv178/b333sPgwYMxZ86chF9Lq9UmlDeJZh45NcYYNm7ciEOHDiEQCIDnefA8D5fLBZVK\nJR7EcRyX0BkGQk6XUY07D05FssNCrFZrSk+p+v1+AM2nZDm5Z2b16tW48sorUVpaiquuugq//vqr\neN+JEydw5ZVXQq1WY8CAAXj//fcBNDSQgYaG+rx583D++edj+vTpp73y+Y033kBWVhYyMzMxZcoU\nOBwOsUwffvghRo8ejeHDh+Pjjz8Wn/PEE0/gvvvuw6ZNm1BeXo5zzz0Xa9asweWXX46ioiK89NJL\n4mNvueUWPPTQQ3jggQcwcuRITJs2TTxNIuyThQsX4pxzzsGVV16JHTt2iPfdeOONWLRoEV5++WVk\nZ2eje/fuTU4X79y5ExMnTkRxcTFGjhyJd955R7xv0qRJWLx4sfjcoqIibNy4EYwx9OjRA/n5+QgE\nAnj44Yeh1WoxY8aMU+4nQlIt1fkVr8WLF2PWrFm45557xA6A7OxsGAwG7N69G0BDr2FNTQ369u2L\noUOHYs2aNQAaDiA2btyIc889FwCwb9++lNShNTk5OXC5XC3et379eowfPx6DBw/GM888I2bymDFj\nkJOTg127duH111+HTqfDFVdcIT7v66+/xhVXXIHS0lJce+212LlzJwBqXLc3h8OBQ4cOweFwoK6u\nDkqlElqtFnV1dU2+k8LhMDWuSbuKt9NTso1rl8uV0gtDhJ6bP/3pTxg3bhwGDRqE/Px8aDQa8ctl\n6dKlGDt2LCwWC+bPn489e/aIV8/HYjFcc801qK+vx8cff4yKigq88MILAP7bQL/55pvx0ksvoW/f\nvnj77bfxj3/8o9XyLF++HHfccQfmzZuH1157Dfv27cNPP/0EAPjrX/+KiooK5OXlIT8/HzfccIN4\nX319PWw2G8aNG4eBAwdi27ZtuOSSSxAOh3HOOedg1qxZ4heQ3+/HW2+9hU8++QRDhgzBypUrceut\nt4plvuaaa7Bw4UIMHz4c1dXVGDt2LEKhEABAqVTitddew+zZszF37lwUFhaKDfcTJ07gwgsvBAA8\n/vjjGD16NP7whz+gqqpKfO6SJUvE5xYUFOCll16CQqHAyy+/jNmzZ8NoNOKiiy7Ciy++iN///vft\n8A4T0nFSnV/xqqqqwueff47LLrsMo0aNwsKFC6FQKPDwww9j4sSJqKysxJgxYzBz5kxotVo8+OCD\nmD59OubOnYuLLroI119/vXgQsWfPnhTXpqm8vLwWe59XrlyJiy66CC6XC+eccw5mz54tzoby2GOP\n4ZFHHkG3bt0wcOBAPP/883jggQcAAK+88gp+85vfoEuXLnjsscfw66+/YtGiRQD+22GS7Asp01VN\nTQ0YYzhx4gR8Ph+6deuG4uJi8cxKcXExFAoFOI5DMBhs81BNQgQGgyG+BzKJGjJkCLvyyitTtv1n\nnnmGKRQKdv7557OJEyeyO++8k82bN48tWbKEeTweFgwGWU5ODrv77rublPmhhx5ijDH28ccfM71e\nz6qrq5u99ubNmxkAlpuby6qqqhhjjF1++eXsvvvua7EskUiElZWVsQcffLDZfW63m1ksFvbII48w\nxhjjOI716dNHLNcdd9zBALA777yTxWIxVlxczO666y4Wi8XYwYMHGQC2ceNGxhhj48aNYyNHjmSB\nQIAxxtinn37KFAoFq6+vZ6tWrWIA2Nq1axljjO3atYsplUr2z3/+kzHG2HnnnccAsOXLlzPGGJsz\nZw4bNWqU+P/hw4ezaDTKdu7cycaPH8+ys7OZ3W5njDE2cuTIZs+94IILmtSzZ8+ebMGCBa2/YYR0\nIqnOr3jt3r2brVy5km3cuJFt27aNbd68Wbzv559/ZosXL2Zr1qxp8pxdu3axJ598kv3rX/9isVhM\nvL2z1XnatGlszJgxzW4fOnQou/LKK8WyT5w4UcwqweWXX86mTJki/u71epnZbGazZ88WbysvL2fz\n5s1jjDG2du1aBoDV1NS0f0Vk6Ntvv2VvvPEGmzNnDnv55ZfZu+++2+znL3/5C1u8eDF79913mdPp\nTHWRSRpIJMMkO+baZrNh8ODBKdt+JBJBfn4+NmzY0OL9X331FRwOh9irATT07gqnFD744ANUVFSg\nR48ezZ7L/n/P9ZNPPonS0lIADacThQnMT/brr79i//79uO+++5rdt337dng8HkydOlV8nfHjx2P/\n/v0AGk5x9OnTBy+++CIUCgXMZjNKS0uhUCjQs2dP6PV6HDp0CCNGjIDZbIbZbBaP3AYPHgzGGKqr\nq/Hvf/8bvXr1wkUXXQQA6N+/P/r27Sv2PrtcLgwdOlQc88gYQzgcBgD88MMPOPfcc1FRUYFPPvkE\n48ePx4YNG8SevZOfG4vFxOcKDAaD2EtOSGeX6vyKV3l5OcrLy1u8b9iwYRg2bFiz2/v374/+/fs3\nu72z1dlkMjW72M3v92Pz5s1Yvny5eN3KNddc02wWkZPzZtWqVfB6vZg9e3aTxwl5L+RVvOM1Seui\n0Sjsdrs4NNNoNLb4OJ1OJ+73qqoqDBw4kIblkDZJJMMk2bhmjMFms6V0zCLP86f8Q+U4DsB/ZxUB\nGqZvEU4PVlVV4YYbbmjxuULjevTo0eJtwnydLamqqkJubi4KCwub3Sc0yDMzM8XbQqEQsrKyxPtL\nSkrEizDz8/NRV1cHoGHMd1FREQ4ePAig4YtBmOexcR01Gg3cbneTbQjbEWY+cblcqKioEL+wYrGY\neAV9XV0dfvjhB1RUVGDLli04++yzm7yOy+XC5MmTW3yuQK/XU+OaSEJnyK9k64x1bumAXBif21pe\nCk7OG47joFAokJeXJ97WOO+FGVOocX1m2P9fJEapVGLPnj2IRCLw+XzQ6XStfg/r9XrYbDbEYjHs\n3r0bu3fvhtFoRFZWFs4555xm31eEnEqiGSbJAWDBYBDBYPCMZq1oL2azGX6/v9WppXr37g0A+PHH\nHwE0vDFarVa8qGfQoEH46quvxEZiOBzGqlWrwBgTG606nU58va5du6KmpqbFbQ0cOBB2u13cFgBs\n3boVR48exfDhw6HRaLB8+XIADWPVPvnkE4wZMwZAw5RgjVcbKigoaDIXa/fu3XHkyBEADQ18j8cD\nnucRjUbx7LPPIicnB2VlZbjwwguxY8cO8SKnf/3rXzh69CguuOACAA1fLicfaBw7dgwAMHz4cFx+\n+eVYunSp2LDesGGDOK1eMBgUZxwAGsZKCs8V6PV6mrqQSEJnyK9k64x1tlgszS5oLCgoQJ8+fbB8\n+XIwxhAIBLB06VIxLwUn503v3r3BGMOmTZsANM97oXGt1Wo7skppieM4fP311/joo4/w448/YuvW\nrfB4PPD5fKdchtpkMoExhuPHj6O6uhocx8Hj8eDw4cP44YcfaOVGkpBEM0ySPddWqxVAahdgGD58\nOFwuF+655x4UFhaC53lx2p/58+fjrLPOwqWXXopJkyZh/Pjx2LdvH3755RcYDAZ4vV488MADGDVq\nFAYMGIChQ4di/fr1iMViqKqqgs1mg0KhaNIY7devHz744AMwxppNs1deXo4pU6bg4osvxm9+8xv4\n/X6sXr0ab7zxBqZNm4Z77rkHM2bMwNKlS7Fv3z4UFxfj9ttvB9BwkFBbWyu+Vl5eHrZu3Sr+XlRU\nhOrqagANPdnr16+H2WyGSqVCMBjE0qVLYTAYcNVVV2HgwIEYMWIE+vXrhy1btmDGjBniQYbFYhFP\n4wFASUkJ7HY7fD4f5s6di4svvhilpaUYMmQI6urq8Ouvv+LRRx8Vn9t4ppSSkhLYbDb4/X7xlCDP\n83SxEJGEzpBfydYZ6xyNRpsNKVAoFJg7dy5uv/12bNiwAXa7HTzPi50TgpPz5txzz8WoUaNw9dVX\n48orr8Tu3buxefNmbNmyBS+99JJ4pjPRaVpJwywzdrsd4XAYVVVV4DgO1dXV0Ol0p2zoGI1GZGdn\nw+l0inNfK5VKZGRkiGdbTz4jQUhrEs0wSf6lC70N8azv3lFGjhyJxx9/HJ9++imsVqsYnt26dUMg\nEEBubi4+//xzvP/++/jPf/6D2267DZ9//jkWLlyIQCCA8vJy7N69G3/7299w8OBB3HfffZg6dSp0\nOh0uvvhifPXVV02uSr3uuutOubrX0qVLMWHCBKxatQpqtRrz588Xe42ffvppjBgxAt9//z0mT56M\nO++8U+wVnzVrltjbDDT0kO/atUv8/ZprrmnSGzNs2DBMmzYNgUAAV111lTg/tk6nw/r16/HKK69g\n7969mDNnTpM5ZSsrKzF8+HDx99GjR+PWW2+FVqtFv379sH//fixfvhwHDhxATk4Oli1bJo71nDdv\nHkaMGNHsuY1PB06ePLnJMBpCOqvOkF/J1hnrfPfdd4urTDZ22223oaSkBB999BEsFgvuvffeJmfO\nAGDChAlNek0VCgVWrlyJZcuWYevWrZg6dSq++OILLFiwAKFQCJdddhmWLFnS4XVKR1arFaFQCHa7\nHUVFRfB4PGCMoaSk5JSL8gjDGrt27YpwOIza2loEg0FxqE59fT01rkncEs0wBZPguZG1a9fi4osv\nxqpVqzB27NhUFyetRKNR8Dzf4tjA2267DQ6HAytWrEhByQhJD3LMLznWmbSPzz//HMeOHYPVakVZ\nWRmOHDmCcDiMvn37JvxaNpsNx48fR3l5OcrKyjBq1KgOKDFJR4lmmCR7roU5pjvT+L10oVKpWu0N\nEFZ/JIScOTnmlxzrTNouFoshEAiA4ziEw2FEo1H4/f4zvhhRmL0lGAzSir4kIYlmmCQHqQpXdJ/q\nYgbS/rRaLV3tTkgbyTG/5Fhn0nahUAixWAw8z4MxhqNHjyIajcJsNp/R6+n1eigUCgQCAXi9XroI\nnsQt0QyTZM+1MDcpTaWTXM888wytdEVIG8kxv+RYZ9J2wiwrwvLxXq9XXG/hTAgXNHq9XhQWFmL/\n/v3Npn4lpCWJZpgkG9fC3M3UC5JcneliJEKkSo75Jcc6k7Zr3Li2WCwoKio65UWM8cjMzERtbS3c\nbjd27NiBnJwcFBUVtUdxSRpLNMMkOSxEmDyepjUihEiNHPNLjnUmbSc0riORCDQaTZsb1kBDJ5HB\nYMDRo0fBcZy4WjEhp5Johkm2cW0ymVJdDEIISZgc80uOdSZtI1y8KKzI215LlyuVShQXF4MxBq/X\nC6vVSgt7K5juAAAgAElEQVTKkNNKNMMk2Y0QDAbpwjpCiCTJMb/kWGdy5rZv346dO3ciFouJKxa3\nV+MaaLg4X61WIxgMikupn+k4biIPiWaYJHuuA4GAOKUOIYRIiRzzS451JmfG5XJh+/btcDqd4Hke\n4XAYAMSFz9qLTqcTX1sYT0tIaxLNMEn2XHMcB61Wm+piEEJIwuSYX3KsMzkzVVVVYIyhtrYWRqMR\narUaSqWyQxrXQqNamGaNkNYkmmGS7LmORCJ0YQwhRJLkmF9yrDM5MzU1NfD5fIhEInC73bDb7cjK\nyoJS2b7NFZ1Oh0gkgmg0ii1bttCiMuSUEs0walwTQkgSyTG/5FhnkjifzwePxwOv1wulUokePXqg\ne/fu6Nq1a7tvS+gJD4VCYIzhm2++oQY2aZUsGtc8z7frxQ2EEJIscswvOdaZJK6urg5AQyPbaDQi\nKysLOTk57d5rDUC8OC0QCKCmpgaRSASHDh1q9+2Q9JBohkmyKyEajbbLfJeEEJJscswvOdaZJM7h\ncCAajSIcDiMrK6tDt6XRaJCRkQGr1YpoNAqLxQKbzdah2yTSlWiGSbLnmjHWIUeyhBDS0eSYX3Ks\nM0mcx+MRZ/AwGAwdvr0ePXqIcxeHQiG43W7EYrEO3y6RnkQzTLJpp1AoUl0EQgg5I3LMLznWmSSm\nceM6GbPLaLVaFBcXQ6PRIBwOIxaLwefzdfh2iTQlkmGSbFzTakqEEKmSY37Jsc4kMdFoFMFgEBzH\nQaFQJHXqRr1ej1AoBKBhnm1CTpZohkmycU0IIYSQ9OH3+wH898KxZJ7pEBrXjDGaMYS0C8k2rqkn\nhBAiVXLMLznWmcQvEAgASM3MMgaDAYwxBAIBHDhwAPv376ex16SZRDJMko1rlUqFaDSa6mIQQkjC\n5JhfcqwzSYzQuOY4LumNa5PJBKVSCZvNhkAggJ9//hm//PILHRASUaIZlrSp+Hw+H5YtWwatVoub\nbrqpxaVMt27dikAggPPPPx/19fX47rvvwBhDJBJBJBLB1VdfjdzcXJSVleHXX39NVtEJIaTdyLGh\nKcc6k8SksudapVKhoKAAx48fh8fjQWFhIQCgsLAQxcXFSS0L6ZwSzbCk9FzX1dXh7LPPxqZNm/Dd\nd9/h/PPPB8/zTR7z4Ycf4rzzzsP3338PADh69ChmzJiBtWvX4qeffsKePXvEo0i1Wo1IJJKMohNC\nSLuSY37Jsc4kMX6/H5FIBIyxlCw4lJeXh969e8NisaCurg6hUAh79uxJejlI55RohiWl5/q5557D\nLbfcgscffxwAcP311+Ozzz7DpEmTxMf4fD5UVFSIRwbRaBRDhgzB7NmzUVNTgyFDhsBisQAAjhw5\nQr0ghBBJkmMvrhzrTBLj9/vBcRwApGw1z4yMDHTv3h179uyB2+2GXq9HOBxu8Uw7kZdO2XP93Xff\n4dZbbxV/P//887F58+Ymj5k6dSry8vLECd2PHz+OjRs3YsKECXjmmWcwdOhQWK1WAMCAAQMoqAkh\nkiTHhqYc60wS4/V6xcZ1MqfhO5larYbBYBCHqdDsIQTopI1rp9MJs9ks/t7aUaDdbkdBQQEAIBKJ\n4JJLLsHWrVvx+eefY8KECXj//fcBAAMHDqRTjIQQSfrss8+wbdu2VBcjqWhYCDkVnufh9/sRDoeT\nPsd1S/R6PYLBIACa95o0SDTDktK4Fi4UEFRXV4sXDDTmdruRm5sLoGHoyBdffCE2xIuLi1FfXw8A\nMJvN4DgOlZWVUCgUrf5UVlZ2fOUIIaQFlE//pdVqxV5JQk5WV1cHoOGiRr1en9Ay0x1Br9cjGo2C\n53m43e6UloV0DolmWFI+wZdccgk++OADAA1jqz/77DNcccUVzR6n0+nEwtfU1ODIkSMAgFgshhUr\nVmDw4MEAGpZIpSVKCSFEGkwmE2U2adW+ffvE3mthaGgqCZ164XAYHo8nxaUhnUGiGZaUCxrvuece\njB07Ftu2bcOhQ4dw2WWXoby8HLt27cLKlSvx5z//GXPmzMG3336LH3/8EYWFhfB4PJg2bRquu+46\nbNy4EQaDARMnTgTQcEGj1+tNRtEJIYS0kdlspswmLbLb7Thx4gRsNhsAIDs7O8Ul+m/jOhQKwe12\ngzGW1BUjSeeTaIYlpXGdn5+PX375BWvXrkVBQQEGDRoEoOEUkDCu6aGHHsJdd90FACgqKoJKpcKa\nNWvw9ddf47e//S0uu+wy8VRRJBJBIBDA3LlzZXVqlRAiHZWVlS3m07x58/DYY49h7ty5yS9UiphM\nJgQCAUSjUahUqlQXh3QiwmqIDocDmZmZnWJmDrVaDbVajXA4jEgk0ml61EnqJJphSVtERqPRYOzY\nsU1uGzZsGIYNGwag4Wj15CPW0tJS3HHHHc1eKysrC0DDGO2cnJwOKjEhhLQ/OeaXHOtMTi8SiaC6\nuhputxuxWKzTfDYUCgV0Oh1CoRCAhs8tNa7lLdEMk+Ty58J813SakRAiNXLMLznWmZxebW0tIpEI\nXC4XdDodMjIyUl0kkV6vFxvXTqczxaUhqZZohkmycW0wGABAHFJCCCFSIcf8kmOdyekdO3ZMHHaR\nmZnZqcY16/V6xGIxcBxH0/GRhDMsacNC2pPRaATQsKITIYRIiRzzS1jnQE51Jq2LRqPYuHEjDh8+\nDI/HA8aY2DPYWTRuTFHjmiSa25LsuRYqKaygRAghUiHF/KqqqsJVV10FnucBAAsXLsTkyZNx/fXX\n47rrrsPDDz8MoKHR9Le//Q2jR4/GnXfeiRMnTgAAzjrrLADSqjPpOIcPH8bhw4fhcDhgtVphMBig\n1+tTXawmdDodFAoFQqEQvF6v+Nkn8pRobkuy51roBaHxe4QQqZFafh05cgSjR4/G8ePHxavkd+/e\njV69euGSSy6BUqlEv379AABPP/00Vq5ciSeffBLr1q3D9ddfj/Xr16OkpARarVYydSYd69ixYwiH\nw6ipqYFCoUBxcXGnGhICAEqlUryokTEGl8uF/Pz8VBeLpEiiuS3Jnmvhql0KakKI1Egtv7p27Yr1\n69dDo9EgFosBABwOBywWC1auXIkffvgBZrMZjDG8+OKLWLZsGUaNGoWHHnoITqcTe/bsgUqlwoAB\nAyRTZ9JxGGOw2Wzi6fWysjJxCEZn03gZdIfDkeLSkFRKNLcl2biW45hFQkh6kFp+abVadOnSBVqt\nFmp1w8nOEydO4M0330T//v1hs9lw8803o66uDkajEV27dgXQMJ1Znz59cPjwYQBAt27dUFFRkapq\nkE7C5/MhHA4jGAxCrVZDq9WmukitysjIAM/z4DgO1dXViEajqS4SSZGzzjoLjLH0HnMtzIdNR5KE\nEKmRYn7Z7XYUFBSIv0+cOBHff/89pk+fjueeew4//fQTQqGQ2LMt4Hm+yenUyspKKBSKVn9oUbD0\nJ3zug8Eg9Hp9pxsO0pjQW2mz2WC1WrFjx44Ul4h0tNNl1Pvvvx/X60iycW2xWKBQKOB2u1NdFEII\nSYgU8ysQCIiNZACYO3cuioqKmjwmPz8fTqdTvOCHMYadO3eipKQEAFBXV5e8ApNOq76+HrFYDKFQ\nqFPNa90SnU6H3Nxc2O12eDweVFVVgTGW6mKRFBLmPj8dSTauFQoFLBaLpL6cCCEEkF5+McbA8zyc\nTiesViuAhsa1w+EAYwxPPfUUhg4dCqPRiHHjxuGpp54Cz/N44YUX0KNHD/To0QMAxOcS+WKM4fjx\n4/D7/WCMiUOkOrOuXbtCqVTC5/MhEAjQXO0yFw6H43qcJBvXQMOpVZp7khAiRVLKrwMHDmDKlCnI\nzs7G/fffDwAoKCjAoEGD0L17d6xYsQKvvfYaAODZZ5/F1q1bkZOTg08++QRvvPGG+DqZmZk4dOgQ\nGGOt/tCwkPRms9ng8/ngdruhUqk6fc810HAwbDAYxB5LqfzdkjNTWVnZaj6VlpaKMyOdjiSn4gOA\n3Nxc6gkhhEiSlPKrrKwMv/zyS5PbZs6cienTp8Pv9yM3N1e8vaCgAJ9++mmLryOlOpOOcejQIcRi\nMbjdbmRlZUGplEb/XuNpJGnGG/lKJMOk8cluQW5uLpxOZ6qLQQghCUuH/NLr9U0a1qeTDnUmZy4a\njeLYsWNwu92IxWLIyspKdZHiptPpEIlEEIvF4PP5Ul0ckiKJZJhkG9cmk4k+5IQQSZJjfsmxzuS/\n7HY7wuEwvF4vNBqNJIaECITpAjmOo8+wjCWSYXE3ridPnox58+Z1mtN6FNSEEKlKdn4tWrQI06ZN\nw+7du5O2zZNRZsvb8ePHwRiDz+eD2Wzu1FPwnYwa1wTooMb1okWLcOzYMfTt2xczZ87EoUOHzriA\n7SErK0tS88QSQogg2fl19913o0ePHrjoootw7bXXYsOGDUnbtoAyW95qamrg9/sRjUabTOsoBULj\nOhwOw+fzNZvPnchDIhkWd+O6Z8+eePPNN7FlyxZ8++23KCsrw6RJk/Ddd9+l5IOWn58Pr9cLnueT\nvm1CCGmLZOeXxWJBZWUlDh8+jFgshgsuuACjRo3Ce++9l7SpxSiz5cvn88HlcsHr9UKpVEpiCr7G\nVCoVNBoNwuEwotGoZFZXJe0rkQyLu3EdiUTwwgsvYOTIkSgvL8fGjRsxfvx4zJkzB+Xl5UmfWD3R\ndd4JIaSzSEV+ffnllxg5ciSOHj2KL774ApWVlVi2bBl69OiBvXv3dvj2KbPlq7a2FkDDe280GqFS\nqVJcosTpdDpxOr5///vfdJAoQ4lkWNxT8V177bUIhUL44osvMGLECADAsGHDMGXKFPj9/qSPn8rM\nzAQAOJ1O5OTkJHXbhBDSFsnOr8WLF+Ott97CggULcP3114tToI0bNw6BQAAGg6HDy0CZLV91dXXg\nOA7hcDihGWY6k4yMDFitVvA8D4fDgQ0bNmDMmDGpLhZJokQyLO7G9aBBg7BgwYIWjzhTcYonLy8P\nAGhqJ0KI5CQ7v7KysvDVV1+hV69eze5L1qwNlNnyxBiDzWZDIBAAkJr2QnvIysqC1WrF3r170aVL\nF6hUKnAcJ47HJukvkQyLe1jIK6+8glAoBMYYgsFgypcAFf5A6cpdQojUJDu/VqxYgX379oExJl6U\nleyhfJTZ8iS0F4LBIJRKJXQ6XaqLdEZ0Oh169eoFxhgCgQAYY3SgKDOJZNhpG9f/+7//i+zsbHi9\nXuTl5cFkMqGoqAg9evTAsmXL2l7aMyT0tghHw4QQIhXJyq9169YhPz8f33zzDa6++moYjUYUFhai\ntLQUf/rTnzp02yejzJYnt9sNAAiFQtDr9ZKagu9kGRkZyM7OFsdeezyeFJeIJFMiGXbaYSHCMre9\nevXCvn37Os1YucZjXwghREqSlV8XXngh9u/fj9tvvx1Tp07FNddc06HbOxXKbHkSLv4Kh8PiBWFS\nptVqwfM8YrEYXZwrM4lk2Gkb11qtFlqtFgsWLEB2dnbbS9dOhIsiKKgJIVKTrPxSKpXIysrC7bff\njvPPP79Dt3U6lNny5PV6EYvFwPO8ZIeENNZ4QRlqXMtLIhl2ysZ1IBDAa6+9hj/96U8YNWoU1qxZ\nA7PZDKVSiUgkgqFDh0KtjvuayHaVlZUFAHC5XCnZPiGEnKlk5dezzz6L2267DcOHD8e2bdtgsVig\n0WjAcRzKy8vFnphkoMyWJ7/fD47jACAtLv4TDhA4jqP5rmUmkQw7Zcs4Go3CbrcDAP75z3/ip59+\nEo9CdTod3nrrLRQXF7dDkROn1Wqh0+lozBMhRHKSlV8ulws8z2Pjxo34+9//Do/HA57nodVq8eCD\nD+K3v/1th26/McpseWrcuNZoNCkuTds1Xq1RaA8JU1uS9JZIhp2ycW02m/H4448DAObOnds+pWtH\nGRkZdHEMIUSSkpFf8+fPBwDccMMNuOGGGzp0W/GgzJafQCAgLriSDo3rk1dr9Pl8sFgsqS4WSZJ4\nMyyuw63vv/8epaWl+PTTT3Ho0CFMnDgRU6dO7RTT8dFpGUKIFCUrv+x2O8466yw8+OCDCAQCuPPO\nOzF+/HhUVVV1+LZPRpktL5FIBOFwGDzPQ6FQpGwYaXvT6/Vi+8dqtaa4NCSZ4s2wuBrX3bt3h0Kh\nwLXXXovp06dj+PDhOHHiBJYvX97mgraFwWBIeQOfEELORLLyKysrC06nE/feey8WL14Mh8OBwYMH\n44knnujwbZ+MMltehB6+SCQCjUYj6Wn4GjMajQiFQuA4DkeOHEl1cUgSxZthcR1G9u3bF0ajEf/4\nxz+wefNmfP311/D7/aitrW1zQduCgpoQIlXJyi+VSoVrrrkGzz33HL777ju89tprsFqteOmllzp8\n2yejzJYX4b3meT4thoQIMjMzUVdXB6fTCa1Wi0AgkLSVTklqxZthcY/Cf+GFF7Bw4ULcdddd0Gg0\n+OabbzBs2LA2FbKt9Ho9wuFwSstACCFnIpn5NWfOHGzYsAE5OTkYMmRIyvKbMltehJ5rnufTZkgI\n0HBhm9FoFC9sO378eIpLRJIl3gyL+9N+6aWXYu/evQAaTvE88sgjGDFixJmXsB2oVCpEIpGUloEQ\nQs5EMvOrpKQEGzZsAAAwxnD11VejZ8+eSdl2Y5TZ8iKMTeV5HmazOcWlaV9GoxF1dXXirGq9e/dO\ndZFIEsSbYXH3XD/44IMoLCyExWKBxWLBn//8Z/zzn/9sUyHbSqVSIRqNprQMhBByJpKZX8uWLUNp\naSkyMzORkZGBm2++Gc8//3xStt0YZbZ8BAIBHDlyBJFIBLFYLC3muG5MGAYSDAbhcDhSXBqSLPFm\nWFw914wxvPbaa9i+fTuKioo6zUUJnaUchBCSqGTm17vvvovFixdj8uTJYIylLDsps+XhxIkTWLt2\nLaLRqHgKPd0a13q9HgAQCoXgdrtpvmuZiDfD4vokKBQKzJ07F9OnT0/J9E2tUSgUYIyluhiEEJKw\nZObXAw88gCeeeALr169PyvZaQ5md/hhj2LJlCwKBAKxWq3jxl9AYTRdqtRoajQahUEic75qkv3gz\nLK7GtdfrxTPPPINvvvkGvXv3Rp8+fdC3b198/vnnbS4oIYSQjrVo0SLs3LkTo0ePRlFREQYMGIBH\nH3001cUiaai2thZOpxNWqxV1dXWw2WzQ6/VpNVuIQKfTIRQKAYhvSWwiH3ENCzGbzThy5AiUSiW8\nXi/cbjcAoKioqEMLdzrRaDTtTjURQuQhmfn1r3/9CwqFAsFgEC6XC+FwGN27d0/KthujzE5vjDFs\n27YN4XAYLpcLjDHodDp069Yt1UXrEAaDATabDYwxOJ1OFBcXp7pIpIPFm2FxDxA6ePAgnnzySWg0\nGkQiEUycOFFc0jRVOI6joCaESFIy88vv92Px4sWoq6tDt27dcN1116Vk8QvK7PTmcrngdDphs9mg\nUChQXl6Onj17QqfTpbpoHUKv14MxhlAoRBc1ykS8GRZ34/oPf/gDdu3ahTlz5qB3794YNmwYvvnm\nmzYVsq1CoVDajeMihMhDMvNr8eLF2LlzJyoqKqBQKDB79my8+eabSdl2Y5TZ6a2urg5Aw1BSi8WS\nlkNBGms8Y4jNZkMsFktxiUhHizfD4m5cHz16FK+++ipWr16N48ePo1u3brBarW0qZFul26pPhBD5\nSGZ+HT16FFOnTkXfvn2xYsWKlOU3ZXZ6czgc4HkePM/LYsVCjUYDjUYDv98PnufhdDpTXSTSweLN\nsLgXkbnxxhsxc+ZMXHvttaisrMQPP/yADz/8sE2FbCsKakKIVCUzv2688UY8+uijuOuuu7Bo0SL0\n6NED48aNS8q2G6PMTm8ul0u8wE8OZygUCgWMRqO4WE59fT1yc3NTXCrSkeLNsLh7rp944gmMHDkS\nW7ZswcGDB/GXv/wFAwYMiLtAX3/9Na6++mpMnz69xbF+0WgU8+bNw7fffivetmHDBlx33XW49dZb\nsXv37mbPoVOMhBCpSmZ+XX311Zg3bx4+//xzqNVqDBo0CFOmTEnKthujzE5fsVgMXq9XVo1roGGl\nRp7nEQ6HUV9fn+rikA4Wb4bF3XO9bNky3HnnnZgxY0bChVmxYgVmzZqF5557DocPH8a4ceOwY8eO\nJq3/iooKfPHFFygvLwcAbNq0CZMmTcILL7wAn8+HK664Atu3b0dmZqb4nHA4nLYXShBC0lsy8+ur\nr77C0KFD8fHHHydle62hzE5ffr9fXDRGo9FApVKlukhJIQx/CQQCdFGjDMSbYXH3XM+cOVM8Ik3U\nX//6VyxduhQTJkzAzJkz0bt3b6xatarJY+bOnYuKigpxWclnnnkGzz77LK6//nrcfvvtuPzyy/HR\nRx81eQ5deU4Ikapk5tdLL72EX3/9NSnbOhXK7PQlTNErtwMonU4HpVKJUCiEYDB4xu0kIg3tNlvI\nkiVLMGDAAAQCAfTr1w+FhYXIzMxEZmYm3nvvvbgKs2PHDlxwwQXi7+Xl5Th48GCTxwwePBgulwvZ\n2dkAgC1btmDMmDHi/f369cOhQ4eaPMftduP555+PqwyEENKZcBzX4Y2QH3/8EWeffTb+/e9/43e/\n+x26deuG7OxsmEwmzJo1q0O33ZJk1JmkhjCvdSgUktV7rFAomiwmIxxkkPQUb4adtnE9adIkfPbZ\nZ8jMzMS3336Lffv2wWazwWazoaKiIq7C8DzfZD32aDQKk8nU7HF2ux2FhYXiYxqLRqMwGo1Nfq+s\nrIRCoWj1p7KyMq7yEUJIeztdPoVCoSaZ1hGGDh2KDz/8EKNGjcJTTz2FXbt24cSJE3C5XHj66ac7\ndNsni0ajSakzSQ273Q6O4xCLxWAwGFJdnKTS6XTgOA5AwzSEJD0lkmGnbVzn5OSgb9++uPfee1FW\nVibOXanRaKBUxjeqpHv37jh8+LD4+65du1pcySgQCMBsNgMAiouLceDAgVaf4/P54to2IYR0Vi11\nMrQnnU6H8vJy3HDDDbjsssuQlZUFnU4HtVoNtTruS27ahZDZHV1nknyxWAxWq1WcNUMO0/A1JjSu\nGWPUuE5jiWRY3Ol63XXX4YknnoDX60UgEMCRI0dw8803Y/Lkyad97rXXXounnnoKL730EtasWYNd\nu3bhoosuavIYnudhMBhQU1ODvn37YuLEiXj66acxfPhwbNu2DStXrsSLL74oPp7GNRFCpC5ZMypc\nddVV+L//+z/U19cjEomguroa/fr1w4IFC5KyfQCym0VCTtxuNziOg9/vh0ajkd24eqG+wj4g6SmR\nDIur65kxhosvvhiMMZSUlGDgwIH44x//iMsvvzyuAj388MPweDzIysrCjBkz8M4770Cj0eDFF1/E\nn//8ZwANw09sNhvuuusuHDt2DHfccQe6du2K3NxcXHfddXj11VeRlZUlvqbP58P8+fPxzjvvgDHW\n4g8NCyGEpEplZWWr2SSclRPO1HW0P/7xj9i+fTtKS0vRu3dvVFRUYOrUqUnZtkDo9UlWnUny2O12\nAA1nnzMyMpoMA5UDYeYzjuPorHoaSyTD4uq5VigUGDRoEPr06YObbrop4QJZLBYsW7YM0Wi0yfQ8\n48ePh8vlAgB89tlnzZ73yiuv4MUXX4RSqWz2xyqceqFTjIQQqUl2fg0aNAgulwt33HFH3MP5TvbL\nL7/g3HPPFbPY4/Fg06ZNGDhwoHitDNCwFPSGDRtQVlbWZCgfZXb6crlciEaj4DhOnJRAThr3XFPj\nOn0lkmFxpSzHcRgyZAhuvvlmTJgwAc899xxefPFF7Nu3L6GCnTzvZa9evTB06NDTPqelo2A6xUgI\nkapk51dpaSnefvttDBs2DIsWLcIrr7yC77//Pq7nRqNR3H///Rg+fDgikQgAYP369Rg8eDCefvpp\nnHfeeeL82Tt37sTgwYPx5JNP4uKLL8arr74qvg5ldvryeDwIh8MA5Pn+qtVqKJVKcBwHjuPEfUHS\nSyIZFlfPtTAc5Nlnn4XX68WxY8egUqnEq2NTQejxbryoDCGESEGy80uv12PBggUIBoNwuVyw2Wwo\nLS2N67lHjhyB3+8Xp59ijGHGjBlYunQpLrnkEuzduxdXXHEFrrnmGtx///2YN28ebrnlFlitVpx9\n9tm45ZZbYLFYEAwGAVBmpyOv1yu2B+Q0DZ9AoVBAq9WK+2Dr1q0wGAwoKyuT5cFGukokt0/buH79\n9dfxu9/9rtmcqGvXrkVZWdkZFrHthJWQ8vLyUlYGQgg5E8nKr6+++gplZWW49dZbm9y+Y8cO5OTk\nxPUavXr1wssvv4ylS5dCrVajqqoKAHDJJZcAaFiDQK/XY+/evdi8eTO+/PJLAEB+fj4GDBiArVu3\nYvTo0eL0qpTZ6SUajSIQCIDjOCgUiiYrL8tJ47muDxw4AMYYHA5Hk/U6iLQlktunHBbCGMOsWbNa\nHEM0c+bMhIeFtCfhCKLxRY6EECIFycqv5557Djt27Gh2+8svv9zidS6tcTqdyMvLg0KhQF1dHbp2\n7drkfovFgqqqKnTp0qXJmG6z2YxAIACgYSw2Ywz9+vU7w9qQzigQCIAxBo7jxOERctR4Or5gMAir\n1YqamhrwPJ/qopF2cvPNNwOIL7dP2XPt9XoRDAbx7rvvwmw2Q6vVIhqNIhKJwOVypfRDIxxByPHi\nCUKItCUrv06cOIE1a9bAZrNBp9OJjaB9+/ahvLw87tdxuVxib012djZsNluT++vq6tC/f39x1ojG\ntxcUFIj/r6ysxPz581vdzrx582iWJ4kROt/kvrS9Xq8XG9Y+n088qHQ6nejSpUuKS0fidbqMAuLL\n7VM2rgOBAPR6PY4ePQq3241QKASNRiOGtNPpTKzU7ai2thbZ2dmy/mMmhEhTsvIrEAjA5/Ph559/\nRjAYhEKhQEZGBtxud0L53bg3sk+fPqipqREbzjt37oROp0OvXr2QnZ2N7du3Y9CgQairq0NVVZXY\niN+zZw/NFJKGGjeu5fz+CgvneL1eeL1ecRiU2+2mxnUa0ev1ceX2KRvXOp0OOp2uxWVyA4EA6uvr\nz/phhqAAACAASURBVLyEbWSz2cQeEUIIkZJk5ZdOp8OsWbMwcODAJre/8cYb2Lx5c1yvUVtbi4UL\nF+LQoUOYP38+5s2bh/vvvx/jxo3DxIkT8cYbb4hLvT/yyCO44YYbMHnyZLz//vu46667xAu6qqur\nMWDAgHavI0ktn8+HWCwGnudl3dml0WiQmZnZpF0Ui8Voar40E8/S5wCgYIyxUz1g586dOOuss5rd\n7nK5kJGRkbI/pssuu0ycT5UQQqQkWfm1Z88elJWVNZsGleM4BIPBuK569/l8WLNmDWKxGLp06YLz\nzjsPALBu3Tr8+OOPuPTSSzF8+HDx8Zs3b8aqVaswYsQI8aJHgDI7Xa1ZswaHDh3C/v37UVxcLOvZ\nYGKxGFwuF/x+P1wuF8rKytCnT59mK1ITaUokw047W0hLDWsg9RcSOp1OdO/ePaVlIISQM5Gs/Gpt\nXLVWq427Y8RkMmHChAnNbr/oootabDQMHTq0xfULKLPTjzA8VJglQ47T8DWmVCqRk5MDg8EAl8tF\ni8qkmUQyTLKX9dpsNrqYkRAiSXLMLznWOd35fD4Eg0EEAgEolUrZN64FwoFrOByG1+vFaQYIEIlI\nJMMk27h2OBw0XyohRJLkmF9yrHO6O3bsGICGi/iMRmOLqynLkUqlglqtBsdxiEaj8Pv9qS4SaQeJ\nZJgkG9eMMfj9fllfmUwIkSY55pcc65zOYrEYNm3ahC1btogLyMh5rHVLdDqduAy6x+NJcWlIWyWa\nYZJsXAvzR8Z71SYhhHQWcswvOdY5nR07dgwHDhyAy+VCXV0dVCoVLBZLqovVqVDjOr0kmmGSbFwL\np1ioF4QQIjVyzC851jmd1dbWIhKJ4OjRo/D5fOjatWuzGWnkTq/XIxKJgOf5lK4JQtpHohl22tlC\nOiPhCMJgMKS4JIQQkhg55pcc65zOHA4HgsEgAKCkpIR6rVsgzO8eCoWocZ0GEs0wSfdc0ylGQojU\nyDG/5FjndBWLxeDxeMTp9+g9bZler4dCoUAgEIDb7QbP86kuEmmDRDNMkj3XLpcLQHzru5P4xGIx\nBINBZGRkNLvimzGGaDQKpVIJhUIh3h8MBqFSqcRph2KxWJNlkhvjOA6xWEx8fiQSgVarhVrd8BF0\nu90wmUxQKpXgeR5KpRKxWAwajQYcx0Gj0Yj3hcNhaLVaRCIR2O125OXlQaVSieUKBAJQq9VgjCEc\nDot18vl8UKvV0Gg0Yll8Ph8YY8jMzEQkEhEDked5MMZgMBgQDocRCoWgVCqh1WrF50ajUfH11Go1\nHA4HzGazuM/C4bD4OGG1U4VCAY/HA61WK9ZRKK/RaEQkEkEsFkMkEhHr3Hj/CvtOrVYjFos1ORUb\ni8XEbQs/KpVKfF8avzfBYBBarRaMMfF9Ff4vvF9qtbrV97PxNoX3kjEGjuNanY4rGo3SqWPIM7/k\nWOd0JazIGA6HodFo6G+6FSqVCjqdDsFgEIwx2O12FBYWprpY5AwlmmGSbFy73W4AoFNRcXC73TCb\nzadtJIXDYbhcLuh0OrHBK+B5Xuyl0Gg04mmR48ePQ6VSoaSkBNFoFD6fr9Urxv1+PzweDzIyMqDR\naOD3+5GRkYHs7GxEo1EcPnwYRUVFUKvVsFqtKCwshMvlQk5ODhwOBxwOB7Kzs8EYg9frRffu3bFn\nzx7U1taif//+CAQCsFgsiEQiCAQC0Ol04DgOHMchOzsbWq0Wu3fvhslkgsVigUajgUajQU1NDWKx\nGEpKShAKhZCbmwu3241AIIBgMIji4mJ4PB7YbDYolUp07doVPp8PNpsNGo0GPXv2FBvGe/fuRX5+\nPkwmE7xer/gYpVKJgoIC8DwPo9GIw4cPIzMzE4WFhThx4gRsNhuysrKQkZHR5CDEZDIhEomgW7du\n4j7jeR6BQAC5ubngeb7J34DT6YRSqRQPHoCG+VaFhnJGRob42Lq6OigUCqhUKpjNZmg0GvH1FAoF\nQqGQeEBwKsI+1mq1iEajCIfDLT5HeF/ob1ae+SXHOqcr4eI8ISNI64xGI5xOJxhj4vcakaZEM0zS\njWua+uf0hB7MeEQikVZfQ2isRaNR8Xahh/nkbbU016nQwyk8p/HrxGIx8Dwv3sZxnNhbLvwEg0EY\nDAaoVCrwPI9IJIJIJIJQKASO4xCJRBAOhxGNRhGJRKBQKBAOh8WeYKER2LhnWKlUimVqXL7Gry28\nRuPHR6NRsWdb6I3lOE7cvlAGYVvC84Teb+H2WCwm9vYKDdRwOCz2JAuvJbx/sVhMvE34vfG+F3q7\nG98nvHcnP7bxeyc8r/H7JzTwT0d4bOP3sjWnuk9O5JhfcqxzuhJ68MLhMJ2JOI2MjAzY7XaEQiHU\n19enujikDRLNMEmOuaag/n/svXl8VNX9//+aJdtksi/su1BBqWKVAkZQ6wLuuKOiouBCrfULtm61\nAipWLFarrVvrQhWLbd0LFFRAQJAPIAjIvichIZkks2T2mffvj/zO9czJubNAwmQy5/l45DGZe889\n933O3HnP677v+5wTH3x6QLxl9RCFESsvpiLoEQqFNDHJBCJ/LC8iRUHIC1F2HC8y2bGsHNvPXvk2\nMAHNjmfCGoAmmIPBYKu2sXLsRkBsC//HBCcT2Ow9O4a3j/+fnZfvL5nI5e2SfR68GGfpHnwdfF/x\n1wj/GfL1xEIU8nokcqPXmUlH/5WObe6sNDY2wufzIRwOa4P2FHLY00K32w2bzaYCDClMWohrNTgm\nfuIVM6Iwk+0XxZFYPtqxQKTAFqOdTGTydohiNJowZft5Ac4LSF4Q83bxQp8/n3i8bLsYDRZtFO3l\n7WP7efv16hU/I74uvb6O53OXCexYn2UsotmlhHUL6ei/0rHNnZFwOIyjR49qMyfwqWaK1mRmZiIj\nIwNutxvBYFATaIrUI1EflpLiWn2x4yORyDUrn2j98Z4rloDjRaFMYPOCXIxo8wJZTxyzulnkWSay\n+TK8WGavLJ2D2cNSRHgbRJEvtkHvPZ/qIUaPxT6VRZr1+ldvuyjkj0f8ttf11VlJR/+Vjm3ujNTW\n1sLn88HpdCIjI0PlXMeBxWLRrv/6+vokW6M4VhL1YSkprh0OB7Kzs7X8UoU+iQifaGX1hJtelFV2\nvCwaDEAqjsVziFFpXuyKApzlPPMiXKyLF8GytBO+LJ/Wwf74FA6+DjGqztcn2iIey3Leo0WnxZsF\n2b54PnPxBkb8jI9FMCd6DaUr6ei/0rHNnZHt27cjEAjA4XCgoKAgrnEZ6Q4bqB4MBpW4TmES9WEp\nK67VqPP4SFRcR0s1kP0fj7jSE596NgCIEH68EGbwOduioNaLPLNX9oPAIs9EpM3SwUQzi27r3VCw\nXGa+fbKoNP8nSz8Rc8/5KLbYJ7Lodzyfk957WX/zZeMV6fFeA0pYt5CO/isd29zZaGhoQE1NDWw2\nGwCguLg4yRalBmx2Lbfbjbq6uiRbozhWEvVhKSmuo82lq/iRRKOQidRzLHXqCT5Z/WIUViay+Qi2\nWK8oHvljouVOi/Xwgw0BtErf0Iuyy1JCRDv48/D16PVXtD7Ta7PeZ3csUepoKOEcP+nov9KxzZ2N\n/fv3IxwOo6GhAQUFBerzjJOcnBwYjUa4XC64XC44nc5km6Q4BhL1YSkprvXm0lW0JpHItZ7IZPtj\nib1YKQF8OdlgQJnNskg2X4cs8svXL0aKZZFfvWP02i+mqYj2ymwUxb3sPX/zwI4VI9SyiLWsj0W7\nZTOOxEoLiTdyLesjvbJKhKen/0rHNncmiAhVVVVwuVwIhUIoLCxMtkkpg9FoRG5uriaqDx8+nGSL\nFMdCoj4sJcU1W4BE0XZEi2LqRVOjiXHxeJmYlZ1bFuEVbRLFL4te8yKdF6niXMxiPbJ0Ff78LGot\ni0SL5+f/9OqMJqBleeJiX+r1i3gDEAu9zzXRiLaejQo56ei/0rHNnQmn06lFXplYVMRPQUEB/H4/\nXC6XtviZIrVI1IelpLj2er1qfs04OBaRE0/0OVrZWMeLUVx+n0wc6glR8Vj2KkapxRxmUZiLOdLi\nFHl6udSydA+Z8JWle/BlZDnoorhNpJ9lwl9Wnm9HtPpioSfQ47U3HUlH/5WObe5MVFVVAWgR2Var\nNeaKv4pICgoKkJGRgZqaGtjtdnz99ddobGxMtlmKBEjUh6XkNyQQCKhR5wmQyOP9eIVcotFKUTzr\n7dMTqaKgBtAqjUI8ThTIvNhm0WZZXQAiZhvht+mJV9FG2XaZOBfbw+xnRMsHl71GK8Nvl0XRY9UX\njXg+/2OptzOSjv4rHdvcWTh06BA2b94Mt9sNv9+PvLy8ZJuUchiNRvTo0QM+nw979+5FMBjEnj17\nkm2WIgES9WEpKa7D4bC6c44TPUGbaDm9/Xq50np16EVoo4lvPZEtE+WiwJZFlnm7RaErE/F8feKi\nLzLBLFuYRhatZ3bwx/Dnk/WrXj1if8j6WFY+WplEbsriLa9IT/+Vjm3uDITDYWzYsAEulwuHDh2C\n2WxWq2weI3l5eejduzdCoRCcTieqq6uVz0whEvVhKevt1PyabUu8kWv+fz0BKDteFKKy+mUDCsUc\nZPEVQKt8Z72yotjmU1T4pdTZftnNgF70mZWTCWRZhF0U9LJott5gSb3+lkW8o30uog3HSjyiPN6b\nvHQhHf1XOrY51WlsbITH40F9fT1CoRD69u0Lk8mUbLNSFqvVCpPJhObmZjQ3N8PlciXbJEUCJOLD\nUlJcG43GCBGjkHMsYiaWSJZtizd6KROYYh16glMvAi3Li5aJaX7gInvPv+qdgx/0KLYvlpAX26PX\nD2J79QZ88m3Wq1fW5/HYLSurItftQzr6r3Rsc2egoaEBQMsczfn5+dqczYpjw2AwwGKxaEtpq0Vl\nUodEfVhKimuTyaQcdQIkIpLijVwnGomUiUmeWAP/gNYpKHpCmT9OJrb1xL2eKNYT7Xp9oRfJlZ1P\nXMFRbIdeH+udQ7wxiFZW7/hjoS0i4OlCOvqvdGxzZ6CpqQmhUAiBQEANSG0jLBYLfD4fwuGwdvOi\n6Pgk6sNSUlxnZGS0WjhEcfzEmr5NJtJkkd9o9ctSG0Sxq3duvgw/vzW/n7eBT/XgBacobnmRyw94\nlB2n1w490c3naPO2RBO4elFwti8W0WYe4cuIKTfxiHA9m2OVV+L7R9LRf6VjmzsDDocDPp8PANQ8\n5W1EdnY2iAher1fNGJJCJOrDUlJcZ2Zmwu/3J9uMlCDeCDMRRSzpLatDzDcyGAxx5yCx42Xn4OsW\nhatsVg0xTYMJVqPRGCFKY9nH18vKseP4FRb5MmKUXJxBhLeBr1M8l140G/jxpkAcPCGLzOu1h+9T\nWR/w2443F1Y8Xgno6KSj/0rHNncGlLhue9gTAK/XC7vdnmRrFPGSqA9LSXGdlZUFr9ebbDPSCl4o\nMsRobTRxzurghS9DFJisHrE+UeSLkW9eEPNl+MgxO594Dl7ssmN5ccvKiZFjUcCLgyvFvpPNICL2\nE2+nrB490SwTzLKIsdhu/jMRI++xkN046JVTtJCO/isd25zq+Hw+eDwe+Hw+GAwGZGZmJtukTkFG\nRgaMRiN8Ph98Pp/6XqQIifqwlBTX2dnZ6oKMg2iCV4SJyljzKovEW79Yjyw9hH8v/q8n+njhKYrk\naHbrpWeIAyRlqR+yKfjEMqwvxcg3bycv/KP1Cw9/gyLrz1iDIfntbS14lYCOj3T0X+nY5lSHpSyw\nxTPUbC9tg8FgQFZWlvZEwOFwJNkiRTwk6sNSUlyrR4ztQ6yooyzaLJbRQ4xYi9Frve18vXppEbKI\nsnisKD5lkV0+Ys1HsPlz8QJZTO+QRdnZ/9FytVgqCauTRcz5+uO9adDrO1m/iE8UZH0WDzI79c6p\nSE//lY5tTnVqa2tBRHC73WqWkDZGievUIy3SQnJycuDxeJJtRkqQSLQhHoEs1h3PADtWlhelYgRW\nJnZlQlgUgbywE3OUeQEcCoUiRKxMeMuOE9NFoolXXozzqSd8G0KhUCsxK35G0aLP0Z4w6KXuyGBt\n5us6lsiULN1Ej0SfcnRW0tF/pWObU5lQKIR9+/bB5XIhFAqpVRnbmKysLPj9foTDYZV3nSIk6sOU\nuO7kxPv4P1oKiSzCyR+TqB2iMBQHMTJYBJmdm8815oW3mKbBn0Nv4B8/aJEvJwpUsb1ivrYooGWR\ncFGoy/qFPy5W9FdvUKisjay83j69iHMiQlsNaEyMdPRf6djmVKa2thYejwcNDQ3IyMiA1WpNtkmd\nCn5QY1NTU5KtUcRDWojr/Px8+P1+7bGKQp+2zpOTCSmZuBTLiGV5eAEdLc+Yr4uvWzb1nCh09dI8\nZHWx+oxGY1ThyCLizD4xn1o2K4gsX5o/t94NAmuTrA/4ftIT3XrT84npOKLQj0W85ZXg/pF09F/p\n2OZUprq6GuFwGE6nEwUFBWrp+jaGiWuPxwObzRb3E2BF8kjUh52wb0ysH9d4I6wA0LdvXwBQc0TG\nIFGRBMQ3jzLPsQ5olCHmVMv+l9URa9BhtHxsvTJMZIuDEnnBy+dwy5Y9l13T0X6kePv5qD2PXsRZ\nvDEQ26miyx2HoqIiAOnlv9KxzanM0aNH0dzcDCJSKSHtQEZGBjIzM9Hc3IxgMAibzZZskxQxSNSH\ntbu43rlzJ0aPHg2LxYJx48bhyJEjrcq88cYb6Nq1K0pLS/HHP/4RAHDw4EH06tULPXv2RJcuXdC1\na1esXLkSANCnTx8AaiBAWxPvNGpMTMpye2OlNMjeMzHJBCov2FlkmP3PiGeBFV7Iiqs58oKa5WLr\n5VOLEWk+JYSPbouzdYjCV9zG3wTopW7I0kZ4u0TEVBTxGJkN4rZEiecG63jq72zk5+cDSC//lY5t\nTlX8fj/sdjvcbjcMBoMazNgOGAwGWK1WuFwuEBEqKyuTbZIiBon6sHYV10SE66+/HlOmTIHL5cJ5\n552H6dOnR5RZvXo1/vjHP2LDhg3Yu3cv3nvvPXz99dfIzc1FKBTC9u3bceDAARw6dAjnnHMOAGjT\noaiBANFpayFzvAMaWR1iegZvJy+m+fMyMStOoSeKSCaUeXv5aDC/yAvbx9uiF+WWCVLxHGJUWUyz\nkKW3sP3ijQUv4PWI90ZGlk4SbeaXWPXrlVHCOT4KCgoApJf/Ssc2pyo2mw1ELbOEZGdnw2QyJduk\nTkleXh5CoRBcLhe2b9+Oo0ePJtskRRQS9WHtKq7Xr1+PvLw8TJw4ESaTCffffz8+/fTTiB/zN998\nE48//jh69uyJwsJC3HHHHfjwww9RV1eHrKws3H333aioqMDjjz+uTWdWWFgIAGogQBwkMs+1HnpC\nlm3jz6V3vCiq+fLiIEFxUROZMOcFqFgX/yqL2jKhzR8nvoqiW0+M8ucTI9+yfhdtlonoaKkc/Awf\n4jGijbGeRLDzyuo71gGN0a63ts7/T1XY48V08l/p2OZUpa6uThPXFosl2eZ0WqxWK8xmMxoaGgAA\nX3zxBfbv359kqxR6JOrD2lVc7927F6eccor2Pjs7Gzk5OXC5XLplunfvjqNHj6KmpgaVlZU4//zz\nsWDBAmzduhVvv/22VgYALrzwwvY0P62IV/jIhJRe2oKI3r5oU8GJ6RH89mjRZrF+MQWDr082dZ7M\nDr1BgXwbotkpE/ri3NuiTbIofVtMmRetnYmmcCQS4VbRbaCiogJEhObm5mSbcsJgIi2d2pyqVFdX\nw+PxIBwOIzc3N9nmdFqMRiNKSkrgcDiwY8cOuN1ufP/998pHdlAS9WHtKq5zcnLgdru19+FwGIFA\nIGKAhFjG6XSitLQUQ4cOxRdffIHJkyfjpJNOwq233oq1a9cCgPaFnzFjRkT0TfybMWNGezavw5OI\noDlW8XM8cxfrCVRxnyiEoyEOOIwmIHnBLKZs6PWHzEYxks2LZla3eEw02/nIugxZWote/+hFwRNJ\n6YkX9aMQSSz/NG/evGSbeMzccccduPHGG3H11Vdj/PjxePPNNwEALpcLjz32GCoqKvDEE09oU1eV\nl5cDUOK6o9PQ0ICGhgbY7XYtL1jRfpSVlaF79+4IBAKor69Hc3OzSg/poDDd2SHE9eDBg7Fhwwbt\n/ebNm9GnT5+IH3qxzJo1azB48GCUlpZizJgx2vaqqiotLK9ywOInniinLCc4Vnl+cCHbFs+ASDGd\ng2EymSL28SkisnbwdoqDFcUIrxilFVNI2P+yGwVeoIuClk/v0JuNgx+MKfYzn0uul74i68N4Ith6\nbefrihbBTjQtRJY6JJ5TpYVEksorFq5evRoXX3wxpk6divvvvx/jx48HAEyaNAlHjx7F888/j4MH\nD+LRRx8F0CIiAEQ8tVR0PL7//nuEQiE0NjaioKBA/da2MwaDASUlJSgpKYHT6UQ4HMahQ4eSbZZC\nAgsKx+vDzO1pzKBBg1BYWIjf/e53uOiiizBt2jRMnjwZAHDkyBGUl5dj4sSJuPzyyzFgwAAcOXIE\nH374IWbPno2VK1fib3/7G55++mn88MMPeO6557Bo0SIAKnczGeiJsFjT5OnVIUsn4aO9YpqEXiQ5\n1oA/FgGWLTTD2xEKhXRn+mDIVoBkBINB6cBB/iaEH5DJYAKb1c1WkuTrEaPl0fo4mvAW+5S1Sa/d\n8XzPxDz4WKjo9o+wMSSpSFNTE4xGIz777DMMHDgQo0ePxoEDB/Ddd99h165dMBqNGDRoEPr164e5\nc+ciKysLRqNRLSTTgTly5Aiqq6tx9OhRhMNh7YZI0f7k5+fDZrPB6XSisrISZ555ptI5HQw2a068\nPqzdp+L7/PPP4XA4MGvWLNx5552YOnUqAGDEiBH46quvMGzYMMyfPx+vvPIKli5diqVLl6KkpAQj\nR45E37598Ytf/AJPP/00/va3v+H0008H0LLGO4BWy1mLfyotRH9hF5FoA+D09okR23giqaw+fhub\nIUQmOsVtDCZSxQgtnxYirvDIEIUj/6ieF9DRZhRh28RFZMRzyISnWAd/zfKznfDl+f/5qQP10kKi\niXCx38TBnOL/0ZClm+idN91+LGbMmCH1SyxiPWrUqCRbeGwEAgHU1dXhrbfews9+9jP873//w6xZ\ns7B161aMGDFC+w4VFhYiJydHG7BVWlqKBx54IJmmK3QIhUJYv349fD4fbDYbCgsLtYVOFO1Pbm4u\nzGYzHA6HtjKmomNRVFQEIuoYU/EBQHFxMf785z/jiy++wL333qs9Ztq7d682IHH06NH4/PPPsWDB\nAk1Am81mzJw5Ezt37sSKFStw6aWXanVmZGQgLy9PjTxvQ+IVxgxZ9FfcLtvP3uvlIcv26aUvyAY0\n6tkhLirD18cPZBTznsVIrmyxGbZdZid/kyAbOCnObhJP3rcesaL4sj7Si+TztsRTZywSbUtnJtX9\nFxFh+vTpWLx4MW699VbMnTsXn3zySavvItAixPm8XTVOpmNy4MABOJ1OHDlyBEajEV26dEm2SWkF\ny29n+bx1dXVJtih9ieWjWAZFLJK2pqnZfHwZKWVlZaipqWkjazonbSl+om0/lvMArQfw6aU/8FFg\n2cwboliUbY+VWhGtDbyAltkmq0t2bjHSL7ZbFLuy1AxZhDqedhxL1DjRY5Rwjp9U9l+ZmZl47rnn\ntMgm+3706dMHu3bt0srV1NTAbDZr5dTguI7Lvn374PV64XQ6UVZWhoyMjGSblHZYLBYEAgH4/X61\nYmMHJt6c66SJ6+OltLRUXYBxEE9aiEyQ6r1n28TcYXYuvfp5gawnGGURbJbeAUQOCpQtaCNGoqMR\nTRTriV3eLhFZagVfViaG+Si63oqXsoGJsrxqfnCkaHuiN0zHQjx9n25pIdFIZf/ldDrxu9/9ThMC\nM2bMwLhx4zB06FA4nU588skn8Pv9eOSRR3DDDTdox5WWlibRaoUeXq8X9fX12gwhbOIAxYmF5fR6\nvd64l9hWnHj42e2ikbLimi0dqmg7YokfUdCKaRP8qx56olaW78uEmphfLR4bK5VCL+2BP16WliET\n8HwEW5xJhC3hLlvGnb9ZEAW32PZo4l6vD8W2RruJibbgDX8jJPZdtPPFQqWFRJLK/is3NxeNjY3o\n06cPunfvDiLCY489BoPBgH/961+YPXs2SkpK4PP58MQTT2jHWa1WLF26tNWNH/+n0kJOPFVVVSAi\nOJ1OWCyW436qrDg2srOzYTAYtCcIqTzgOZXRGytDRDj//PPRtWvXuOpJ2W9Rly5d8O233ybbjA5N\nIoImWlk+6iwSbzRSL7Kpl2LBz7Ihs43PW+bhZ9tgM4DIoqnhcFg7hygmeXHMzyIiq4e3IRwOa2MK\n+HKynHDxvGJdsr6TIfadGOGOdgxDFnWPFxWNPjZS2X8ZjUb85S9/wTPPPAOgZaYDxk9/+lPddqVy\nmzszhw4dgt/vh8fjiVs4KNoeo9GIzMxMeL1e7WZHPUXoWCTiw1I2ct2lSxc12foJQjb4MN5cX7ZP\nJgJlx0ab2o+ftk6MDovlRMEsimNRgIqRa7aP/YkRbpmo1+sf8RyytkZrw7FGe+OJOvN57KId8Z5X\ndsN0rDalC53Bf+Xn50cI61h0hjZ3NhwOB2pqarTBtQUFBUm2KL3JysqC1+sFEP8y24oTRyI+LGXF\ndWFhIVwuV6vFTBQ/opdKoVeWf5Xtl9XHp0/Eyrfl6xDLxjPThZiSIZ6T1c2nq+hFrWX2ydortpu3\nQUzhEIU+HzkXzy2Lwov79dqg18+yG41Y6SP85xerT/Tq4NuW6Kwz6Uo6+q90bHNHZ8eOHQiFQrDZ\nbMjLy9OmuVUkh+zsbPj9fhCREtcdkER8WMqK60RXy1HoE02I6qUpJBKBFEWqmGYizsbBjmFpGQy+\nfLQoskxgiiKY1a0XsWXnYmX0xCjfLn6ObFnKh3hDINrJ33SIgxbF8mJUWxzQKL4X0ctfP1ZEWxXR\nSUf/lY5t7siEQiEcPHgQdrsdwWBQLRrTAcjOzgYRqUGNHZREfFjKiuvi4mIAaj7IWLRltFAmBqoY\nqQAAIABJREFUmsQ0imjCihe0oqgVI8Ji5FacDYS98lFqfpo+McodbS5s2aIweqkqol3sVVxKXSaG\nRfiBmnzaidh+8UYh2tMB8X9ZFJuvQ/b5xZviIdYZq9zxpLh0NtLRf6VjmzsytbW1CAQCsNvtyMzM\nhMViSbZJaQ+/CqDNZlNPeToYifiwlBXXbOBFbW1tki3puMgEVDSiRR71RLq4JHgi54iW+iBGt4Ef\nRTI7Tq9dovjlzy2L+MpsY+9jLe8uDryMNhMHKx9tYKRscKFYVpY7LqaxRBPIvA3RotvxoFI9jo10\n9F/p2OaOTFVVFcLhMFwuF/Lz89V3uQOQkZGBjIwMNDc3IxAIYMGCBVi+fDkCgUCyTVMgMR+WsuKa\nzZlaX1+fZEs6Pu3hNPlZOQD9RU3YNjEnWZaHzB/Pi1TZ3XusVBDRHpnwlEWbZakpoiiXRWtFG8V0\nD7Zflq4hps3otSlWH0Tbr/e5sFd+QGa80WhZXbFSi1TkuoV09F/p2OaOSiAQwKFDh+BwOEBECQ1M\nVbQfBkPLSo1OpxPhcBiNjY2orq7Gzp07k22aAon5sJQV12xUc7zrvKcjiYqjaI/towk9FglNdPAk\nI9qxspQIWX62uD+eumX7+TrYecRVIWURdlYXIJ/xhD+PLB1FzMOOFTHXa4ve5yQbsMmno4j9EM/n\nqYTysZOO/isd29xR2bx5M3w+H2w2G7KyslRKSAeiqKgIoVAIe/bsQWVlJZxOJ/bt26f8bQcgER+W\nsuKa3Wk7nc4kW9KxiWeQWaL7Es3LFXOQxbQNWU6zbGo8Me1CrI/t519l5dg2sR5mU7T0Ctk0f7Jo\nt14fiHbx+2V9K/YLj156C1+33mdzLPnV0VDR6cRIR/+Vjm3uiOzZswe7du1CfX093G43ysrKVEpI\nByI3Nxc9e/bUFvNpamqCy+XC9u3blW9NMon4sJQV12zUpoqCtA2JpBiI2+MRVuKARDENQzagkf9f\njBDHEqNsnyw/WZbCIm6TTX/H9vODKBn8vNuyKelkM4ToiWg+oizLaY/Wx7L/xeP4QZPxDoLUI9YT\nD9E+9ePQQjr6r3Rsc0fD5/Phu+++g9PpRE1NDQoKClBYWJhssxQCRUVF6N+/P0pKSmC32+FwOLBp\n0yYcOHAg2aalNYn4sJQV1zk5OTCZTMpRRyGR3FkxyivuO97IBp8aIYs2i+fmo7yiPbJUDf54mYDV\nE9RiOV7Yinbw52GvbC5r0X69HGxejIrLussi03rbZe2MFqkX+6+taOsp/dKFdPRf6djmjsbu3bsR\nCARQU1ODzMxM9OzZU0WtOzDl5eXIysrCwYMH4fF4sH379mSblNYk4sNSVlyzxH81Z2ps4nGesfKS\nZQJXlpIQKxorRqHZfhahFSPUssiqbBYP3h5Wl96sGuKUfbJ6ZMcajcaI6LpeG0WxzQtd2T6xHvYn\nm9IvmjDWi5bzbZD9z9usl7YiQ1ZnLLGthHcL6ei/0rHNHYlgMIjdu3fD6XTC6/WivLw84dmeFCcW\ns9mM/v37w2g0orGxEU1NTermNIkk4sNS+puVlZUFn8+XbDM6LPGKGr1UiHjgnXM0gc6LT9l0c7LB\nfuKAO3HhGJlIFNsqRoxFUc/Xx9I9xAh0tHPJziETzfy5+VxzcZ5u8Rx6Nyt6AljvGL3Phk9hEW1U\nEa32JR39Vzq2uaOwf/9+eDwe1NXVISMjQ80QkiKYTCZYrVZNVFdWVibZovQmXh+W0uLaYrHA7XYn\n24wOTaICSRaRFeFFKSCfKo9HFMTioi6y84kDBvXSPWS5w3wUlk2DJ6Z28O2VzT6iJ16jrQwptoPf\nHwwGI9JZxBsZ2bSAYj/I7Ba3ieeXRefFeqNFx+O92Yp1IxdvNDydSEf/lY5t7ijs3r0bHo8Hzc3N\nKC0tVVHrFCI/Px+BQAAej0eJ6yQTrw9L6W9Xbm4umpubk21GhyURIROvYOThI72JnJMXWsFgUCrw\n+AGC7FUUnHqD/cQp9GRtleVEi6Jf70aDF/yy6fT487A28WkvYn0MNqUha5ueMBZtYcdGi2bzr7I2\nHk+UOp4nHioS3pp09F/p2OaOgMPhQFNTExobG2E0GlFUVJRskxQJkJeXB4PBALvdjvr6epValUTi\n9WEpLa4zMzPh9/uTbUanIdo8yQy93N9E83Oj5ebqiTWZkOdTR2TRWj7tgp+lhD+feF69tBA+8s5g\n5xFXXjQajdJoOt+OUCgkFflim/Qi/XxfMPEqi17H6mfZtIWJkugNliI9/Vc6trkjUF1dDaBlCjGr\n1QqTyZRkixSJYDabYbVaYbfbAUDNGpJE4vVhKS2uMzIy1LKgbYyeMBKjxgx+6rlYgyKBFtHJC08m\nJsXBgryIFUUkAN1UDz5NRIw0y8qx/dHSLkTByrdbFKd8fbIIMl8vi2gzoS0T8LEi1zJiHSvaKMub\njyfKrNe+aOWU8P6RdPRf6djmjkBVVRW8Xi/8fr82nZgitSgsLITf74fL5cK+ffvadNYnRfzE68NS\nWlybzWYEg8Fkm9Gh0cvZ5dETznr74hF/PPzS2kDrqLc47R1vuyi4+Tr0RLheHrGe6NWDnzJP1h/8\n1IIivC287Xq51bGQpXxEGzwZS+Cz//XSZhJFrFcRm3T0X+nY5mRBRPjhhx+wePFi1NbWagPilLhO\nTfLz82E2m2Gz2eByubB69Wp1o5oE4vVh5hNgS7vBRzsVrUmkb+JJ69ATcaFQCBkZGTHrMRqNWo41\ng8+tjhalFdMkZPAR9Fg3DLwN4XBYe0waCoUiHpmyfWKkWxT27FyyAYvMLj7vWszBlkXnYw361OsD\nPcRUF700mXjPxdsbbznFj6Sj/0rHNieLffv2YdOmTfB6vcjKykJjYyNyc3M1X61ILYxGI0pKSlBb\nW4uamhrtd+Tss89OtmlpRbw+LOXFtYqS6ZNIJDGRyCMRJRx9EkU1s0+MgrPcZV70BYPBiGisGNFm\nNunZLUaIxYivOOOILNVDPAef0iETp/z0gWLahzitoKwNsgizeB5ZKgo7dzzfDd6+4/0e6aWEKPRJ\nR/+Vjm1OBkSEHTt2wOPxYM+ePdqj7O7duyfbNMVxUFpaCr/fj7q6Oi2QdPLJJ6OkpCTZpqUN8fqw\nlE4LUT/m8ZFIP+ldNLwwZHXy+c+x6mCRWjE9QxZh1kt34G0Qo8P8q5h6wbaL0WJe4OrNliKKYLF+\nXojzkWG+Tr1BhuIc0yJiSom4Tw9Zyoqsr2SDIPk6YiGzQeVcx086+q90bHMyaGpqgt1uR0NDAwAg\nOzsbvXv3VikhKY7RaETPnj1RUFCA+vp6BINBbNy4UaVanUDi9WEpLa7VD7U+iYgxJuJi5Vzz6QrR\n0i5ExEco/JRzehFPmdjlB0LKcqhlgxR5Ia73KEdWp8FgiEgPYXervCDm7eOFtSxXnM87529K+JQX\nvei82A+ySDvbHi3fWq+/xTYca841I9bMJooW0rFP0rHNyeDQoUMgIjgcDhQWFqJv374oKChItlmK\nNqJLly4gIhw5cgR1dXVYsmSJmoXnBBGvD0tpcR0KhdRE+FGI9yKIJ7IYLUp5PDmUvFDXE+16IhFo\nWb1KL91E770owPXSLWTpJPxxYg44H53nX0UbxFfxOP5/WT5zrBslvfKx0klk9cdzDYliPVYkXAms\nFtLRf6Vjm080RIQDBw7A5XIhGAwqUd0JycrKQnl5OZqamnDgwAE0NjZiy5YtyTYrLYjXh6W0lxMH\nnykiiffxhRjNjCV+9NI1ZPtFIat3DLNXFq0Vo7pAZNTXaDRGpKaIedvRorh828U0EbZfjDDLouV8\nXez8fMqH3vza4v96i9iIbRD7Q+8zkLVbb15rvX6Kh3gFuhLWP5KO/isd23yiqaqqQnNzMxobG7X5\nkRWdj7KyMnTr1g1OpxN2ux379u1T6SEngHh9WEqL62AwCLM5pcdkthuJRDvZ/mhRUj3hFWslQUCe\nTy2mRgDyVBBWXjZ4kLdZ770sSs0LWb7NLO1EFMXRlgPnZ97QE6p8e/nj9aYQZPtkEXDxWPF/MfKv\nd9PE31CIbRaP0+NYRLPKuf2RdPRf6djmE0koFMKmTZvg8/ngcDhQVFSknhR0UgwGA0pKSrSZYAKB\nAKqqqpJtVqcnXh+W0t86v9+PzMzMZJuR8shElx5McLIZPMTjY0VKZcuF8+/ZTBesLj5aLNrMi0Ix\nCs0LZr0VDvm2sv18lJqvJ55IcLRcZlk9eqKZ9QNrn2xQpQyxP0R7ZVF82Y1ItJsJWfvaS4h3dtLR\nf6Vjm08ku3fvhsPhwJEjR2A0GlFaWppskxTtiMFgQGFhIZqbmxEIBLBv375km9TpideHpbS4DgQC\nas7OOIglZvhBerLVBmWiURaBBVoLLVE8y84tK8/EJy9ceZHJz8ghHiNG4WUrSIqCmIgipvzjxSzf\nP7IIPlstUmyH2B7xWJY+ws/1zR8rbov2hCBaSoh4nHjjcTzo5d2LfaQi1q1JR/+Vjm0+UdTU1GDb\ntm1wuVxwOp0oLy9XTwnSgMLCQgBAQ0MDjhw5gsbGxiRb1LmJ14eltLj2+XzIyspKthkdnmOdVk3c\nLwrlRISVKDBlx4rHiUKaF6Kycux/XmgDP4pIfuaOcDisiVsxpYTfzrdBbJPegEXeDtYucUCjmIoh\n3gjw7dZrb7RZOaKlpMSKRsvs0isbDypSLScd/Vc6tvlEsGfPHnz11VdwOp2orq5GZmYmiouLk22W\n4gSQmZmJ/Px82Gw2hEIhbNy4ET6fL9lmdVri9WEpLa69Xi+ys7OTbUaHRJZGEa0sE3t6qQ1A5KwW\nYrQ12nkS2RYtf5hhMBhaDdzQm9lDNsUdfy4iajUHt3huMZVD7C9xP983YkRdFOBiyoysj/g6Zcii\n9dEi3LI0FfF84ja9uqK919unxHYL6ei/0rHN7U1NTQ3Wr18Pp9OJXbt2IRAIoGfPnirXOo0oLy9H\nOBxGTU0Namtr8emnn6Kuri7ZZnVK4vVhKf3tc7vdsFgsyTYj5YmVFsK/xkrtkIltcT9/Dl408tFg\nWRqFuCiLeC69Oaj5tA29vGOxHl5oikJeFNNsPxsQyZcXo/Z6Nw/xpHLIFryJJpLFJwWyc+j1yfGi\n0kJik47+Kx3b3J74fD6sXr0abrcbhw8fRnZ2NgYNGoTc3Nxkm6Y4gWRnZ6O0tBQNDQ3YvXs3XC4X\nvvnmGzV7SDsQrw9LaXHtcrnUNEM6JBJVFPfHikCK0exE8m3FyC8vfHmhzS8xKg7m44WmnlAWI9Q8\negMEZQu9yAZYxoo+i7OjyG4++O167RLtY/3Ct1VvlpJoTwv4/pf1l94MI7L6ZHnieudVRJKO/isd\n29yebN++HV6vF4cPH4bRaETv3r1VTnua0qVLF/Ts2RPBYFCbjnHFihXaKp2KtiFeH5ay4trv9yMY\nDCpH3QboRXzZ/6IIY4KSCb1oKQtifWLetLhPHDzI28ALTzGKzWyS5QnrCXBeiLI/fmpBPlIsy3Pm\ny0Wrn7VZFMNi+ohYPz9oM1YON48sXUUWQedt10tPaYtUj0QEe7qQjv4rHdvcXgQCARw8eFCbHcTr\n9aJbt25qJpY0xmAwoKioCF26dEFzczOqq6tRVVWFL774AjabLdnmdQoS8WEpK67ZiFg2UlahTyID\n08Sy4j62LZoAE48Xo8xitFiMxgKtUzP4/bIp52RikT+nbLo8mfiUpZDw5cW6ZSI8GAy2Eu38flFw\ni/XrRd/F/pDdBIn7xfaKUW69mxZ2/mMVwnrXkBLWLaSj/0rHNrcHjY2N+Pzzz7F69Wp4vV7U1tYi\nOzsb+fn5yTZN0QEoKipCSUkJbDYbdu/eDbfbjbVr1x73zFCKxHxYys7Tc/ToUQAtqxQpWhNN8MrK\niqkbsuNkgjGR6dxkkVZZxDmayBTFudgGvoxoNz/Nnt6gQ74u/rwmk0lX1IrnFFM+xD4Nh8Mwm82t\nBluK7RAj0Hx9etFu8Xi9qD1vp16/ieVF+PLximaVe91COvqvdGxzW+J2u/F///d/qKmp0VJBPB4P\niAj9+vVT3y0FgBYf2717dxQVFWHfvn04cuQIMjMzsXfvXgwcODDZ5qU0ifiwlBXXdrsdQMtdmiI2\n8UauxfeiGOOjw9EWaOGP1xPrbNlymag2Go3aYAwxys3EsSwvm69HTHmIlhsupkzwfzKBzLedF+pi\n+oXsHOLNgiikWVTbZDK1Ev6x0nf47dEi1LInBOJrIk88ol1DsXLy05F09F/p2Oa2oLGxEZs3b0Z9\nfT08Hg8aGxvR0NCAYDCI4uJiFBUVqRlYFK3IyclBWVkZamtr4XK5sGHDBjQ1NeH0009XefnHSCI+\nLGXFtdfrBQA1Z6oOoiiOJmh4ASkOwJOJJj2RqHeMmJLB1yHLvxaFKZt3Wtwn5oqL5dirmBMtrsLI\n28n3hyhAZVFmsV/48qw/+br58uL0fHyKjNhOUayKNxZ6AlmvL/T+EkkF0bsGZNeMiqpFko7+Kx3b\nfLx4PB589dVXcLlcaG5u1gR2VlYW+vTpo2YFUUSltLQUTU1NOHz4MMrLyxEKhXDgwAEMGDAAp512\nGkwmU7JNTCkS8WEpK64dDgcAIC8vr9W+GTNmSP9PZ0TBw/fLvffeq1tWPE4U36yM3jGiEOcjwkCk\nIGb7RSEbT2Q9WqSUF8vsXPzMILy4JSLccccdEf1TWlqKjIyMiLxudiwfuRaFqijSA4GArgjnB4my\n7czxiRF0IDIdJ5qYFVNEWN2yXHdZlDzaU4k//elP2v+TJk2SpqyI9sk+z85KND8UzX91VmK1Wfnt\nFpg/8fv9WLt2LbxeL2677TatT7KysjBw4EB1wypw0003RVw3gwYNSp4xHQSj0Yi+ffuisrISv/rV\nr7T+2bZtG+x2O0455RRtwZmCgoK0z9uP5YMS8dspK67Z9DIlJSWt9s2cOVP7P12ddKzINd9H99xz\nT0TZWDnXMvEa7Vz8sawM2xYMBiMENT/ojxetQMuPjtlsjph5g7eBCXW9engbWeoFE81ifXz/vPTS\nSwCgpapEy6c2GFrmueYj0CzHWpzTO9qNg17EnU/L0YuUi7O5iDceogjWE9CxxDXfR7fffruuaJaJ\n/XQgmh+K5r86K7HarPw2sG/fPqxfv177/gYCARw5cgRAZP+89957yTKxQ6P6qDWZmZno168fgNb9\nw64tRv/+/TFs2LC0fboUywcl4rdPiLg+dOgQdu/ejVGjRiEnJ6fVfiLCmjVrkJubi9NOOy1i34YN\nGxAKhXDWWWdFiJ/a2loAanCMHmKkM5pIkqVXMESRLtYlE2iimGNiWBR7MiHKR7rFNA32v1iGCWiW\no8zXxf/Pi3p+VhC+bdFuDPi28PWLs4MwEcxyysUoPXsfCAQi+l2MivOv4gBMdnPA+kSM+MvqYPWI\nKSDi/2Lb4yXaoNh0FNfR6Mz+a9OmTfB4PBgxYoTy2TEgIuzYsQMHDx6EwWCAzWaDy+WCz+eD2WxG\nbW2tWspacdzInnIcOHAARqMRPp8PBoMBBQUFICLU1tbizDPPRGlpadqKbD0S8WHtLq5feuklzJ07\nF4MHD8aePXuwcOHCiBGrTqcTl112Gfx+P9xuN8444wy8+eabCAaDuP7667F//35kZ2ejuLgYn3zy\niZaIb7PZkJubKxXriviiyXzZaI/0xbxaMRIsi5zyx4s28UI0GAxqnykTjEyE8lFiXiiyMhkZGa3E\npxhVFoWxLK+a/QUCAd0+kgl11la/3y+NBjOb2CsT4WzmETG6HggEkJGR0WqQp6xd/GfDBn7yKS9i\nJF7sJ367eHPEnyteIRzteguHwxG5ffHOLtOZ6Yz+KxQK4ZZbbsHWrVuRn5+P7OxsLFy4UPuB7oxt\nPlZYZHr37t347rvv4Ha7kZOTg2AwiEOHDmnfU5ZbrVC0NR6PBwaDATk5OQiFQqipqYHdbkffvn2x\nYsUKmEwm9O/fH6WlpSgoKIDFYkn7gbOJ+LB2FdfV1dWYM2cO1q9fjy5duuCNN97AzJkz8e6772pl\nXnzxRQwePBivvvoqgsEgRo0ahdWrV2P37t0Ih8PYuHEjDAYDrrzySnz44Ye44YYbAADNzc1qMEcU\nRMEcbaEXsaw4S4conPWWLJe9F2fTEAW23+9HZmZmK/HMbBa3m0ymVrnLvCjnz8tHlkWhqSfa+Txk\nvrwoQplQZ6I8KytL28f3Hy/ymXjn7WT/87OPMFvE8/F1sag4Ox+/+A07TjZ/OLOBHxTJf0ZiKorY\nDj3EzyNa5FqJ687pv/7973/DZrNh06ZNMBqNmDBhAubPn49JkyYB6JxtTpTa2lqsX78edrsdVqsV\nHo8HTqcTBw4cQHZ2tvY96tOnD0wmE3JycqQ+SaE4XgYPHhzx3uFw4PDhw9i1axdycnJQUFCAQCCA\n3bt3a2WGDBmC008//USb2mFIxIe1q7hetGgRrr76anTp0gUAcNlll+Gpp56KKPPJJ59g3rx5LcaY\nzRg3bhxWrlyJtWvX4le/+pXmWC677DKsXLlSE9cOhyPtk++jYbVaW0UP9URSbm6u1s8mkylCoPHC\nnAkjJh75uvkyvHgMBoOtUimYiGT7eNtYBJcX+awePtrNR7b59BBeNGdmZsJqtSI7Oxtmsxn5+fkw\nm80wGo3Izs6GyWTSXs1mM7Kzs6V35nfddReMRqOW6sFeRfjIvqydrG2sv/r37w8igs/nQ3NzM7xe\nr3YDwQt1k8mktdls/vErywS33+/XhDI/vSF/s8DqY48A+Qg8H9Hnnxbwn2+sUeW8eDabzRE3Aqz/\n+TLpTmf0Xx9//DHuu+8+7Vq57LLL8NVXX2niujO2OVFKS0sxbtw47ckXELn4FAscAC03woFAoJW/\nBYCBAwdqfiUYDGp/fr8fPp8v6lO4jozZbEZWVhYyMzORkZGhvZrNZu3VbDZrPjsjI0P7E7nuuuu0\nfo7ls/nfD+afWf+yz4Ctzufz+eD3+7W+5v9Sqd9NJhMyMzO1Pu/VqxeGDh2q/f6azWbk5OQgLy8P\nOTk5yMnJSfsVQBPxYe0qrmtra9GtWzftfX5+Ppqbm2OWaWxsbLU9Ly8Pbrdbe9/c3BxxR8XDR9A6\nMoFAAC6XS4teNDc3w+12o7GxEXa7HQ6HA42NjXA4HPB4PPB4PPD7/fB6vdoXPBAIRIgjo9GIjIwM\nZGdnIzc3V3NU+fn5yM/Ph8ViwfTp0yPsqKmpQW5uLiwWi/ZqMBhgMplQXFzcyulkZWUhEAjAbDZr\n4o99SVmElhduGRkZWnS3ublZ+1yY/QaDAYFAQBOuPp8Pubm5mtAsKytDTk4OLBYLLBYLMjIy0LNn\nT2RmZiInJ0dzwrm5ucjMzNRs4YUoEWn92dzcDJfLBafTidraWjQ3N8Pj8cDhcKCpqSliUAMA3Hjj\njZo45oUygxeNJpMp4kcgMzMTFosFOTk5sFqtsFqtyM/PR0FBAfLy8lBaWorCwkJ06dJFayef3sH+\n+HxvNscmu4HIzMyEyWTSPhuWpsLnZ7MfD6/XG5Haw4S80WjUbqyysrK0HyM+9SYa/fv3j+jrhoYG\nVFZWwuFwwOFwwG63o6mpCfX19WhqaoLT6YTb7dauZ3Zty1KBmG1Go1HrU/5HNSMjQ/sRYKuTseu4\nqKgIxcXFsFqtyM3NRV5e3gmb4zVav73//vu46KKLTogdJ4rj8dlAakzZ2FF8Nu+reZ/NCIfD8Hg8\n8Hq9cDqd2v/sj4lF9p0LBAIR3z/m60ShyV/TvGhl/o59P5lgy87ORk5ODrKysiK+q1lZWcjKyoLF\nYokQ07F8dn19PZqbm1v5bJfLhblz50b00dVXX92uPru0tLSVz2bim/dtbrc7ou9ZX7Nrg91Esd/N\nePqbt9lkMmm/wey3gLXBYrEgLy9PGjRigUrW1w0NDbDZbJrPbmhoSDufzbdBxvvvv4+xY8fGVU+7\niuuioiLs3btXe19bW4uuXbu2KlNfX68tJ1lbW4vevXujqKgINpst4lgWAQdaVquaMWNGKyEksmrV\nKhQUFKC8vBxFRUVtfudFRJqjZSKivr4eNTU1aGhogNPpRFNTE2pra7XJ3Ovr61FXVwen0xnXOZiI\nZHeOLMLKnBj7crGoLouCut1u7Y7a4XBoczSKjvrUU09tdc7c3FwUFxejtLRUu/hLS0vRvXt3FBcX\no7CwUHtljzPD4TCys7NhtVo1Z1NcXIzi4mLdtv30pz+Nqw/8fj+amppgs9lgt9vh8XhQX1+P+vp6\n2O12uFwu1NXVaY7W5XLBZrPBZrNpP3bx5g+L19S+ffu0Pub7Wy+vm0U62A8Yc65soFI0cnJyUFpa\nqjnuwsJClJWVwWq1oqioSLuOi4qKkJ+fD6fTifz8fBQWFsZ1bZeWlsbVByyiTkSoqamBzWbD0aNH\nUV1dDZvNhmnTpkWU//nPf679CNpstqjtNJlMsFqtsFgs2rXNX89GozEiyuT3+zUxwvqU9TH/Q+Vy\nuaKmPzEKCgpQVlamXduFhYXIz89Hbm4u8vPzUVZWhpKSEs3xFxcXIzc3F7m5ubBarXE7+lj+6Ykn\nnoirnlShrXw2AGzZskX57Hb22Xl5edqTPavV2uZzHos+m92EtIfPtlgsrcR1ZWVlh/PZpaWlmlhv\n62ub+Wyfzwen06n57D179iifnYA4byu/3a7iesSIEXjttde0yNjChQsxatSoVmWWLl2Kk046CeFw\nGIsXL8bbb78Nm82GJUuWYPTo0SAiLFq0CPfdd592nMy5yDjnnHMi3pvNZu0iyc7O1u6oRafHp0bw\nj+z8fr8WrfD5fHFdHFarFeXl5ejatSsKCwsxYMAA7f/8/HztDpldDIWFhdqck0VFRW02YjcQCMDt\ndqO6uhpTp07VvoiXXXaZFglgixWwH5TGxkbYbDZs2bIFtbW1aGxsjOtc7IuXmZmpRTDYHSvbx5w5\nn04RCoW0u2KPxxPR17HIy8vT7nitVivKysowZMgQrT8LCgq0L11eXh7y8vK0yA+7iy4NhHGjAAAg\nAElEQVQoKAAQ+QVqy2nBgsGg9oPucDhQV1cHu92u9XNdXZ3mFG02G/bv3w+bzQan0xnzkSNzelar\nFTk5OVoEg0X32ZMBABERqUAgoH3+LJrlcrmkj6MZ06dP1x4b5ubmoqCgAP369UNubi5KSkrQo0cP\nlJaWIi8vT4v6sGgPf/PVlrBIV0NDg3bdNDQ0oLGxES6XC263WxNSR48eRUNDA+rq6rB79244HA6t\nTCwsFovWx+z6ZtE35jsMBgPOO++8qPV89NFHnWrKuREjRmDJkiVaZGfhwoWYMGGCtj9enw1E3nQr\nn618diyfzexqD7+tfHb6+Ox4Bg//5z//ievaMlC8t4bHABHhsssug9FoxODBg/HGG2/g66+/xtCh\nQzF//nyMHz8ee/fuxdixY3HLLbdgw4YNAIAlS5aguroaI0eOxDXXXKMl2W/YsCHiDiTWHcbtt9+O\nm266CY2Nja3ukNmdKcuTYo83ZFO0sUcbfPoBe5zF7kLz8/O1i7K0tFR7zJ+bmxvxqCvV8fl8aGpq\n0h4X2e127Q7f6XRq/zc3N2uPQX0+n/Y4LBAIwOv1wu/3t5oBg0UZsrKytKgP62sWnWURAPZlYe+t\nVmunzwdj0QjmfJhzcTqdEY+imdNlEQOPx6PlejPY41B2XbMfK/bjyq5rdo0XFxejvLwc3bt375TX\nNSMQCGjRM6fTqT0a5QUMe89+2JjgYdEa5lLPPffcqP5pzJgxWL58+QlqWftTV1eH4cOH4/LLL0dd\nXR2+++47bNq0KeKRtPLZJx7ls5OH8tntT1v6bCC23z7vvPPw1VdfxbSrXcU10JI3+sEHH6C6uho3\n3ngjevToAaBlpOqXX36J7t27o7q6GvPnz0fPnj1x7bXXahdAY2Mj3nnnHRQUFGDChAlpPw2MQqFQ\ndGTsdjveeecdWCwW3HTTTbBYLMk2SaFQKE447S6uFQqFQqFQKBSKdCGlJ9DcvHkzrr/+eu290+nE\nI488goqKCsycOVMbDOLz+fDMM8+goqICDz74IOx2O4CWqPpf//pXnHPOObjnnnu01Xc6G3/9619x\n44034pprrsHVV1+tLXdORHj//fdx7rnnYuLEidi/f792zKJFi3DhhRfi2muvxdatW5NlelIIh8P4\n+9//jtGjR2Py5MmoqqpKtklJgYhw66234sYbb8T48eMxfvx4bY56h8OBhx56CBUVFXjqqae0/Eqv\n14unnnoKFRUVeOihh+IeAJaKvPXWW/jDH/6gvd+1axcmTJiA888/Hx999JG2vbKyEnfeeSdGjx6N\nN998U3sEabPZcP/996OiogLPP/98XAN7Up36+npcf/31qKysBNDyXXvjjTcwevRoTJkyJWI55g8+\n+ADnnXcebr755oiB8UuWLMFFF12Eq6++Gps3bz7hbTgRLFu2DDfccAOuu+46jB8/HldccYX2u/Xt\nt9/i8ssvxyWXXIKVK1dqx2zbtg3XXXcdLrjgAnz++efJMj1pLF++HOPGjcOVV16J9evXJ9ucpPH0\n009H/N4/8sgjAFr0zquvvqrpnZqaGgDRdUBnY82aNbj99tu1942NjZg2bRoqKiowZ84cLWfd7Xbj\niSeeQEVFBR5//HEtrzsQCGDu3LmoqKjAAw88EH08A6Uo69evp9LSUurevbu27aqrrqK77rqL1q5d\nSxMnTqTp06cTEdEvf/lLuvbaa2nNmjX0wAMP0I033khERE8//TSdf/75tHr1anrqqado9OjRSWlL\ne/Ob3/yG7rrrLlqyZAktXbqU9u3bR0RE8+bNo9NPP51WrFhBr776Kv3kJz+hQCBAS5cupX79+tHS\npUvp/fffp169epHD4UhyK04cL774Io0aNYpWrlxJzz//PJ1xxhkUDoeTbVZS6N27N7377ru0dOlS\nWrZsGTU1NRER0aWXXkpTp06lb7/9liZMmEAPP/wwERHddddddMMNN9DatWvpvvvuo4kTJybT/Hbj\n2WefJYvFQvfccw8RETU0NFD37t3pjTfeoK+//poGDx5My5YtI5/PRyeddBI999xz9M0339DIkSNp\n/vz5FA6HacSIEfTII4/Q2rVr6dJLL6Xnnnsuya1qXxwOBw0YMICysrI0HzR37lw655xzaNWqVTRn\nzhwaPnw4hcNhev/99+nUU0+l5cuX0xtvvEEnnXQS+f1+Wr58OfXp04f+97//0YIFC6hnz57U2NiY\n5Ja1PZ988gmdffbZtHjxYlq6dClt3LiRiIi2b99OXbt2pQ8//JAWLlxIvXr1or1791JNTQ117dqV\n3nnnHVq2bBkNHDiQ1qxZk+RWnDjWrVtHPXr0oP/+97/00UcfUY8ePai2tjbZZiWFm2++mR555BFa\nsmQJffHFF1RZWUlERH/4wx/o3HPPpdWrV9Ps2bPp7LPPJqIWHXDaaae10gGdjWXLllFhYSGdcsop\n2rZzzz2Xpk+fTmvXrqWrrrqKnnrqKSIiuummm+j222+ntWvX0uTJk+nee+8lohYtdfnll9OaNWvo\noYceoiuvvFL3fCkrrh0OB23atIm6du1KRER79uyhgQMHUigUIqKWH7uioiJqbGykbt26kdvtJiIi\nv99PpaWl1NjYSN27d6eamhoiIgqHwzRo0CDauXNnchrUjtx55530yCOP0G9/+1t67LHH6MiRI0RE\ndNppp9HmzZu1cr/4xS/oiy++oEsuuYQWLVqkbZ84cSLNmzfvhNudLPr27UsHDhzQ3p9xxhm0YcOG\nJFqUPAoLC2nevHn0q1/9il555RUKhUK0fft2Gjx4sHbDUVdXR6WlpWSz2ah79+7k9XqJiMjr9VJx\ncTG5XK5kNqFdOHjwIL388ss0ZcoUIiJ64YUX6IEHHtD2v/vuu3TzzTfTv/71L7ruuuu07V999RWd\nd955tHr1au3HjYho586dNGjQoBPXgCQQCoVo165dNGjQINqzZw+Fw2Hq1asXHT58WCszdOhQ+v77\n7+nMM8+k9evXa9vHjh1LixYtoiuvvJI++eQTbfsdd9xBf//7309oO04Eb7/9Nl1zzTU0a9YsmjZt\nmtYX9913H/3lL3/Ryj3xxBM0Y8YMmj17Nj366KPa9tdff127NtOBm2++md577z3t/f33308vvvhi\nEi1KHuPGjaNnnnmGpk2bRjNnzqSmpiYKh8PUo0cPqq6uJqIWvXPyySfT9u3bdXVAZ6OhoYFWrlxJ\ngwcPJiKiDRs20Jlnnqn9jh04cID69OlDlZWV1KdPH+0Gw+l0UkFBATkcDuratasWaAwGg9StWzdN\nQ4qkbFoImxidzaG8ZcsWjBw5UptrsaioCGazGevWrcMpp5yirQXPFiBZt24dCgoKtHlYDQYDTjrp\nJBw4cCAp7WlPampq8NZbb2lL6l5yySUIBoM4fPhwxJRXrP2bN2/G2WefrW0fOHBgp+wXGXa7HaFQ\nKGJKns56XcTC7XajqakJ8+bNw1lnnYWPP/4Yc+bMwZYtWzBq1ChtxoDS0lKEQiGsX78ep512mjYN\nWVZWFrp166alAHQmevfuDY/Ho/mf77//HhUVFdp+ds2I29l3Sdzev39/HDp06MQ1IAkYjUYMHDgQ\nDQ0NKC4uRkNDg+aPGQMHDsS+ffuwZ88enHHGGdr2aP3ZGR9j19TU4NNPPwURYeDAgbj88stRU1Oj\n65v1rr90Qe97lo7U1NTg73//OwYPHoz6+nrcdNNNqK2thdVq1RZ5MhgMGDhwIPbu3YtDhw5F6IDO\n2ndME/I+++yzz9Z+x3r37o2amhps2bIFw4cP1ybWsFqtKCwsxJo1a9C/f3/k5eUBaJkpp1+/frp+\nOyXmZXnmmWfwv//9D263GwaDAb///e9x6aWXor6+XluUhrilrxnBYBC5ubmttgcCAd3trONSESLC\nNddcg7q6Ovh8PmRmZuLf//43LrjgAsyZMwdDhgwB0JLLePjw4VZzYrL2i30ZCASiLgTTmZBdR6l+\nXRwrBoMBv/nNb/Dkk08iKysLp512Gu666y5MmzatVR/R/7/MeTr1XTT/E+27JNseCoWQm5t74oxv\nZz7//HM8//zzaG5uRjgc1nI/w+EwmpubUVhYCJvNpnu9JOKbUv36euKJJ7B8+XJ4PB4YjUY8/fTT\nOPXUU/Hvf/8bV1xxBQBgw4YN+PLLL3Xb7/F4Ol2/JEJnvC6OlfHjx2PSpEno2bMnwuEwysrK4PV6\nE/6udUbq6+u1gKp4zbDfMLaiMU8gEIDVak3o9y0lxPX999+PCRMmICsrC2azWVvb3ePxwGq1AgD6\n9u2LXbt2acdUVlYiJycHAwYMwO7du7WldT0eDyorK/HTn/4U9fX18Hg8yMnJARFh27Zt6N27d1La\n2BYYDAa89tprcLvdyM7OhtFoRGlpKR544IGIcuFwWLuDq6qq0qZH3Lp1K+666y6tL8866yxt+8SJ\nE094e5JBQUEBfD4f7HY7CgoKQETYunVrXJPLdzZycnIwZ84c7T1bDEr8rh08eBAFBQXo378/du3a\npX3XmpubUVtbi7KysmSY3+7w/qdPnz7YuXOnto9dM3369MHGjRul2xcvXtxqe2fh4osvxqmnnoqM\njAxtwQagZeU2tlx0cXExnE4nXC4XrFar9l3r27cvunXrhoMHD6Jv374AWvrnlltuQZ8+fbBr1y5t\npdGtW7fi6quvTlYz24QHH3wQkyZN0n7fCgoKWs3/LH73hg0bBuDH6yYnJwc7d+7EpZdeGrE9XWDX\nBX+98BH+dOLxxx+PeE9EKCsrQ0NDA9xuNywWi6Z3+vXrh5KSklY6YMqUKckwvd3xeDyaGO7Tpw8W\nLFig7du+fTt69uzZypfbbDZ4vV4MGTIE+/fvRygUgslkQiAQwJ49e7R+a0XbZrWcOMLhMC1cuJBG\njRpFDQ0NFA6HacCAAfTpp5+Sz+ejiRMn0rRp04iI6Oyzz6Y33niDAoEAPfjgg9qAxuuvv55mzZpF\ngUBAG1jTGXnmmWfo4MGDRNSSy3fKKadQOBymadOm0T333ENer5f++c9/Uv/+/SkUCtHcuXNp/Pjx\n1NzcTF9++SWVl5dTc3Nzkltx4rjjjjvooYceIr/fT6+99hoNGzYs2SYlhYaGBvr9739PwWCQvF4v\njR8/nmbOnEmhUIj69u1LCxcuJK/XSxMmTKCHHnqIiIiGDx9Ob731FgUCAfr1r3/daQc0+nw+uv/+\n++nxxx8nr9dLmzdv1nL1Dx06RIMGDaJly5ZRdXU1devWjbZs2UL19fX085//nN577z1yu93UrVs3\nWrlyJblcLrrkkktozpw5yW5Wu2Oz2chisWi5nxMnTqTHHnuM/H4/vfzyyzR8+HAiInr44Yfpzjvv\nJI/HQ//5z3+oT58+FAwG6c9//jNdfvnl5HK5aMWKFVRWVkZOpzOZTWoXPvjgA1q2bBkREW3atInK\ny8upsrKSPvroIxo+fDg1NDTQ5s2bqby8nA4ePEjffvstDRw4kCorK2nfvn3Ur18/+vbbb5PbiBPI\nm2++Sb/4xS/I4XDQ2rVrtTEg6cjvfvc7stlsFA6Hafbs2XTBBRcQEdGNN95IM2bMIL/fT3/605+0\nMR/Tp0/XdMCCBQs0HdDZCIVCNH/+fBo7diw1NTWRz+ejHj160FdffUVut5uuuuoqevLJJykcDtOQ\nIUPogw8+IL/fT1OmTNEGrl944YX04osvUiAQoMcff5yuuOIK3fOlrLhetWoVnXHGGTRs2DD6zW9+\nQ0RE3333HZ111lmUl5dHN910kzaKfM+ePXTeeeeR1Wqlyy67TBs9e+TIEbryyivJarXSueeeSzt2\n7Ehae9qT+fPnU69evahHjx50+umn06ZNm4iIyG6302233UZ5eXn0s5/9THPGXq+XHnjgAcrPz6fB\ngwfTkiVLkmn+Caeuro6uu+46slqtNGrUKPr++++TbVJSCAQCNHnyZOrevTuVlJTQLbfcog0M3rBh\nA51xxhmUl5dHEydOJLvdTkREu3btotGjR5PVaqUrrrhCE1GdjWeeeYaGDRtGw4YNo88++4yIiN55\n5x3q1q0blZeX05w5c7SBMp9++in169ePiouL6eGHH6ZgMEhERCtXrqRTTjmF8vPz6e677yaPx5O0\n9pwImpqaaPjw4TRs2DDtB//o0aN0zTXXkNVqpYqKCtq6dSsRtQxYv+OOOygvL4+GDRtG33zzDRG1\n3NRMnz6dCgoK6OSTT6aFCxcmrT3tycaNG2nIkCHUt29f6tWrFy1YsICISBNMxcXF1KdPH3r33Xe1\nY15//XXq0qULde3aNe0G8wWDQXr00UepsLCQTjrpJProo4+SbVLS+POf/0zdu3enbt260dlnn63N\nzFNTU0NXXXUVWa1WGj16NP3www9EpK8DOhsLFy7UNOOsWbOIiOibb76hoUOHUn5+Pt15551aEHHr\n1q00cuRIslqtdN1111FdXR0RtQx6vOiii8hqtdLYsWMjJj4QUYvIpAl+vx92ux2lpaVaAr9CEQ92\nux1Go7HT5uEpFB2RcDiM+vp6FBcXd8plqxXth9frRXNzM0pKSpJtStqixLVCoVAoFAqFQtFGpOxU\nfAqFQqFQKBQKRUdDiWuFQqFQKBQKhaKNUOJaoVAoFAqFQqFoI5S4VigUCoVCoVAo2gglrhUKhUKh\nUCgUijZCiWtFyhMKhZJ6fnFJ1ETx+XytlqBVKBSKzkqyfTYJS18fy/Fut7sNLVJ0NpS4VnQoVq1a\nhaVLl0Zse+655yKWI+XZvn07Lr744oTOEQ6H8dhjj8Hj8QBocZSJzEj5zTffYPLkyZg1axYA4NZb\nb8WKFSsSsoFn6tSp+Oabb6T7Nm/ejClTpuDhhx8+5voVCoWivThy5AheeeWViG2ff/45Pv74Y91j\n+vbtm5DPBYCXXnoJ69ev194ncvzRo0cxefJk3HPPPXC73Xj99dcxe/bshM7P88477+CFF16Q7nM6\nnbjnnnswZcoUNDQ0HPM5FKmNEteKDsWhQ4fw2muvae/37duH5557Dl26dJGWb25uhs/nS+gcBoMh\nwjH/6U9/wl/+8pe4j9+yZQuampowZswYAC1O3mazJWQDo7a2FqtWrcLZZ58t3d+zZ09ccMEFePfd\nd4+pfoVCoWhPgsEgZs6cqfnUUCiEX/7yl7o+G8Ax+Uuz2awtZLVu3TpMmDAh7mPr6+uxdu1aXHTR\nRcjKygIRHbPwJSK88MILuue3WCy4+OKLsWrVKjQ1NR3TORSpjxLXig7Feeedh2XLlmmPDd99913c\nfPPNKCwsBNASqV66dCnsdjsAIBAIICMjQzu+vr4eq1evjkiz8Pl8+Oc//4l169YBaBHXU6dORU5O\nDhYvXowtW7ZgxYoV+OCDD3Do0CGtbgCoqalp9UOQm5uLIUOGRIhrtoLaf//7X6xZswZff/01Nm7c\niC1btmDLli2or6/H4sWLEQgEIup67bXXcPvtt8NkMgEAdu7cidmzZ+OSSy5BMBhESUkJrr/+euTk\n5Bx/5yoUCkUb06tXL+Tl5WHr1q0AgK+//hplZWUYOXIkgJbI9uLFi1FZWRlxHFsp2OfzYeXKlRF+\nl4iwaNEiLFq0SPPll19+OQYNGoT169dj2bJl2Lp1K95//30cOHAABw8e1I71er3Yu3dvxLksFgvK\nyspw9dVXw2QyIRwOaz57x44dePfdd1FVVYUPPvgAR44cwbJly+D1erFkyRLU19dH1LV8+XL06NED\n/fr1A9ASFX/ttddw6aWX4ocffoDJZML48eNRUlICi8Vy3P2rSE2UuFZ0KLp164bevXtj3bp1ICJN\nXAPAM888g4kTJ+Kdd97B8OHD4ff7EQqFNGG6cOFCnHrqqZg6dSouvPBCEBECgQDOP/98LFy4ENOn\nT8dHH32E2tpanH/++QCAZcuWYcOGDfj++++xaNEi/O1vf8PUqVM1eyZOnNgqTSUcDkcI+kAgALPZ\njM8++wyTJ09GUVERPvzwQzz77LMYP348zjnnHJx11lmYPn06fvvb32rHOZ1OvPXWW7j77rsBAO+9\n9x4qKirgcDhw8OBBHD16FADUcvUKhaJDM27cOCxatAgA8I9//AO33HILAGDp0qUYM2YMFixYgDFj\nxmDLli0RTw3379+PwYMH4//9v/+HU089VfN5Dz74IObOnYu3335bS4k7//zzcfToUWzZsgVLly5F\nVVUVFi5ciEWLFuHCCy+E1+sFADz77LOtUjZEnx0MBmE2m7F7925ceOGFyM7Oxtq1a/H666+joqIC\nEyZMwKmnnopHH30UY8eOjajr6aefxoMPPggA+OGHH/CTn/wE69atg9FoxLZt27RyZrP5uMfjKFIX\nc7INUChExo0bhyVLlgAAMjIy8LOf/Qx2ux3vvPMOli1bhpdffhlbtmwB0OIks7KyAADTp0/Hv/71\nL4waNQpDhw7Ftm3bsHHjRgwcOBBvv/225tR3796tRYKfffZZDBgwAD/88ANeeOEFuN1uDBgwAKtW\nrUJmZiZ27tyJa665JsI+v98f4ahDoRDmz5+Pjz76CKtWrcLJJ5+M3NxcrFu3DmvWrMHIkSPx6aef\nIhwO495779WOe+mll3DLLbeguLgYAPDZZ5/h4osvxpVXXomnn35au2lQKBSKjsy4cePw/PPP4777\n7sPHH3+MP/zhDwCAxx9/HP/5z3+wdu1azaeGQiHNf86ePRuTJ0/Go48+iilTpuDDDz/EBRdcgMWL\nF+O7775DZmam5rftdjssFgsmTZqEESNG4JZbbsE//vEPAC1PDF944QXcd999eOmll1qNYRF9djAY\nxKZNm1BRUYFZs2bh2muvxaJFi7Bq1Sp88cUXeP311zFixAjcfffdKC8vh9/vR2ZmJlasWIFQKIRz\nzz0XQEtwZtCgQRg/fjwuvPBCZGZmaufIyMho9aRSkT6oyLWiwzF27Fj897//xbx583DbbbfBYDBg\n165dcLlcOPfcc1FcXIxvvvkGmZmZ8Hg8yMrKQiAQQFVVFc455xyYTCYMHToUjY2N+PLLLzFx4kQA\nLRFgg8EAj8cT8bguJydHG9xosVjw2GOPYdasWXjxxRcxbdq0CKcMtDzGZI8U2fsdO3ZgyJAhWLt2\nrXau++67D127dkVhYSEKCwtx0kknYd++fQBafijefPNNLQICtAzc7Nu3LyZNmoRf//rXEedUERCF\nQtFRGTNmDNatW4d58+ahoqIC5eXlICLs3r0b1113HbZt24bly5fjJz/5CbxerxYQ2bJlCy699FIA\nwGmnnYaGhgYsW7YM119/vSZU2ZM7j8ejBUV4nw0ATz75JJ5//nm8+uqrGDNmDAYNGhRhn8/ni/Dj\nfr8fGzduxMiRI7FixQoQEQwGA8aOHYuKigoUFBSgsLAQJpMJffv2xcGDB0FEeOKJJzB79mzNpkmT\nJuGGG27AH//4R1RUVMDv92vnyMjISHg8kKLzoMS1osMxatQoHDx4EP/85z+1lJDu3bsDAFasWIHp\n06dj27ZteP/99+F2u5GVlYWMjAyUlJRo+c0bN27EkCFDUFZWpkW5q6qqsHfvXphMpoipoLKzs7VH\nigAwZcoU7Nq1C0uXLsXkyZNb2ef3+7UfBwBwOBx45ZVX8MEHH2DmzJnYsWMHGhsbtYh0eXk5qqqq\nkJubC7/fD7fbjZdffhkTJ05EQUEBgBbxvHfvXjz55JNYsWIF/vGPf0Q8Pk10ZL1CoVCcKCwWC8aM\nGYOHHnoIt912G4D/r737CWn6j+M4/vytzY1ygptryZzbMPNUFDskI1hkKB2UOkgQdDQKBh0EG3Rp\ndIgIoksjBsEQBv0j6NIlElERDykePEgHo3QHMfwDaobuuw7DL7+vW7/T4mf1ehz3+X73/Zxee/P5\nvD/flYpin89HJpPh0aNH2Gw27t+/b2Y2QEtLC+Pj42xvb/P27Vui0Sg+n4+ZmRkMw2BjY4PJyUkA\nS27vzexTp04Rj8dJJpOW1rtdezN7bW2NgYEBnj9/zuzsLLlczpLZfr+ffD4PQCAQYG5ujtHRURwO\nh9lLDjA1NUUikWBoaIj19XVmZ2fNMRXXfzcV17Lv2O12ent7OXfuHIFAACgF3J07d4jFYoRCIfr7\n+zl+/Dgul4sjR44ApTaLs2fP0traysDAAF6vl1u3bvHs2TOOHj3KmTNn+Pz5Mz6fj5MnT1qeuRuq\nAE6nkwsXLtDX10dtbW3Z/L5//27Z/jt8+DAej4dIJMK9e/fIZDI0NTWZW4Jer5fFxUUAQqEQ+Xye\nbDZr6e3e2dnh9u3b1NfX09rayo0bN9RrLSK/jStXruB2u+nu7jY/e/r0KYlEgkgkwvnz5zl27BgH\nDhwgGAwCkEqlePjwIQ0NDTQ1NdHZ2UlPTw+HDh0iFArR1tbG2NgYAB0dHdhspZKlWCxaMhvg4sWL\nxGIxTp8+XTa33baOXR6PB7/fj9PpZHBwkEwmQzAYrJjZ4XCYfD5POp2mv7/f8r0vX76koaEBj8dD\nMBikra3NHKupqbGsZMvf5Z+ilsRkH9rY2KBQKFBXV1c2ZhiGGbKAuaUHpSLVMAxLkAKsrq7idrsr\n9jEXCgUKhYJ5z+bmJpFIhOnpaRobG8uuv3v3LjU1NVy9etVcUa80RwCbzcaHDx9obm42V7AbGxv5\n8uUL4XC47L61tTWcTiculwsobYUuLi4Sj8ctJ+JFRPYTwzBYWlqq+Aq+/8rsYrHI1tZW2RuR1tfX\nsdvtZhbu9e82EYCuri6uX7/OpUuXyq4dHh7m8ePHPHjwgGAw+NPzLLsH5BcWFlheXubEiROsrKzg\ncDj4+vUroVCobNHj27dvbG9vm79VhmEwPz/PzZs3SSaTtLe3V3yW/NlUXIvs8eLFC3K5HG/evKk4\n/urVK1KpFHa7nYmJCct2Y7Wl02mePHmCz+fj/fv3v+w5IiK/q3w+TzQaZWFhwXIeZtfHjx+5fPky\nOzs7ZLNZotHoL5vLxMQE165dw2638/r164qLKPLnU3Etsse7d+84ePDgT//YRbyxP4kAAABSSURB\nVERE9o9Pnz4xMjJi9nuL/N9UXIuIiIiIVIkONIqIiIiIVImKaxERERGRKlFxLSIiIiJSJSquRURE\nRESqRMW1iIiIiEiVqLgWEREREamSHyPMnOy0eXYaAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x160432750>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"with plt.xkcd():\n",
"\n",
" fig = plt.figure(figsize=(12,4))\n",
"\n",
" ax0 = fig.add_subplot(1,2,1)\n",
" for i in range(0, sigma_total_rand.size, 500):\n",
" ax0.fill_between(x, gauss_intrinsic[i,:], color='gray', alpha=0.1)\n",
" \n",
" ax0.set_xlabel(u'Velocity [km/s]')\n",
" ax0.set_ylabel(u'Intensity')\n",
" ax0.set_xlim(-1050,1050)\n",
" ax0.set_ylim(0,0.22)\n",
" ax0.text(0.1, 0.8, u'Each component', transform=ax0.transAxes)\n",
"\n",
" ax = fig.add_subplot(1,2,2)\n",
" ax.fill_between(x, gauss_total, color='0.3', alpha=0.5, lw=3)\n",
" ax.set_xlabel(u'Velocity [km/s]')\n",
" ax.set_ylabel(u'Intensity')\n",
" ax.set_xlim(-1050, 1050)\n",
" ax.set_ylim(0,240)\n",
" ax.text(0.1, 0.8, u'Total', transform=ax.transAxes)\n",
"\n",
" plt.savefig('velocity_dispersion.png', bbox_inches='tight', dpi=300)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.9"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment