From 4e02bc8e650eb6384772775ab508de4f360b2402 Mon Sep 17 00:00:00 2001 From: djperrefort Date: Wed, 31 Jul 2019 16:00:21 -0400 Subject: [PATCH] Updates ew notebook with new plot --- notebooks/equivalent_width.ipynb | 78 ++++++++++++++++++++++++++++++-- 1 file changed, 75 insertions(+), 3 deletions(-) diff --git a/notebooks/equivalent_width.ipynb b/notebooks/equivalent_width.ipynb index a5a0605..1ce0102 100644 --- a/notebooks/equivalent_width.ipynb +++ b/notebooks/equivalent_width.ipynb @@ -17,11 +17,11 @@ "source": [ "import sys\n", "from pathlib import Path\n", - "\n", + "import numpy as np\n", "import sncosmo\n", "from astropy.table import Table, join\n", "from matplotlib import pyplot as plt\n", - "from sndata.csp import dr1\n", + "from sndata.csp import dr1, dr3\n", "from specutils.analysis import equivalent_width\n", "from specutils.fitting import fit_generic_continuum\n", "\n", @@ -39,6 +39,7 @@ "source": [ "models.register_sources(force=True)\n", "dr1.download_module_data()\n", + "dr3.register_filters()\n", "\n", "# Load models for different masses\n", "m102 = sncosmo.Model(sncosmo.get_source('CMFGEN', version=1.02))\n", @@ -71,7 +72,7 @@ "data": { "text/html": [ "Table length=8\n", - "\n", + "
\n", "\n", "\n", "\n", @@ -481,6 +482,77 @@ "fixed_width.head()\n" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "We plot the pw values at peak for feature 7 (Si II 6355) verses feature 6 (si II 5972). Observed values are taken from the CSP DR1 paper." + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "def plot_peak_pew(model_pew_table):\n", + " \"\"\"Plot peak pW7 vs peak pW7 for modeled and observed values\n", + " \n", + " Args:\n", + " peak_pew_table (Table): A table of modeled pw values at Bmax\n", + " \"\"\"\n", + "\n", + " fig, axes = plt.subplots(1, 3, figsize=(18, 6))\n", + " for axis, typing_system in zip(axes.flatten(), ('Type', 'Wst', 'Bst')):\n", + " csp_type_table = dr1.load_table(1)['SN', typing_system]\n", + " is_typed = ~csp_type_table[typing_system].mask\n", + " csp_type_table = csp_type_table[is_typed]\n", + " \n", + " csp_pew = dr1.load_table(6)\n", + " for subtype in set(csp_type_table[typing_system]):\n", + " is_subtype = csp_type_table[typing_system] == subtype\n", + " obj_ids = csp_type_table['SN'][is_subtype]\n", + " plot_data = csp_pew[np.isin(csp_pew['SN'], obj_ids)]\n", + " \n", + " axis.scatter(plot_data['pW7'], plot_data['pW6'], label=subtype)\n", + " axis.errorbar(plot_data['pW7'], plot_data['pW6'],\n", + " xerr=plot_data['e_pW7'], yerr=plot_data['e_pW6'], \n", + " linestyle='', label=None)\n", + "\n", + " axis.scatter(model_pew_table['pW7'], model_pew_table['pW6'],\n", + " label='CMFGEN', color='brown')\n", + " \n", + " axis.legend()\n", + " axis.set_xlabel(r'pW7 (Si II 6355) [$\\AA$]')\n", + " axis.set_ylim(0, 70)\n", + " axis.set_xlim(40, 190)\n", + " \n", + " axes[0].set_ylabel(r'pW6 (Si II 5972) [$\\AA$]')\n" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAABCUAAAF9CAYAAAAz7dKvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xl01OXd///nRRJIWAoSQbOoIFWEBAiLC2KVgkI0YkFagVpB9Aie3ljgVlDL0ri1WPUGbP3+bqUt2tbKVsCGWLlRoLaKshNQRCqyJUEwFBSSISS5fn8kM2aFmcxMPrO8HufkJHN9Zj6f93Dal5+55lqMtRYRERERERERkabWzOkCRERERERERCQ6qVNCRERERERERByhTgkRERERERERcYQ6JURERERERETEEeqUEBERERERERFHqFNCRERERERERBzR5J0Sxpiuxpjt1X6+NsZMMca0N8asMcbsrfp9QVPXJiISLZTFIiLOUg6LiFQy1lrnLm5MDJAPXAv8F3DcWjvHGPMYcIG19lHHihMRiRLKYhERZymHRSSaOT19YzDwubX2APAD4LWq9teA4Y5VJSISXZTFIiLOUg6LSNRyulNiNPBG1d8XWWsLAap+d3SsKhGR6KIsFhFxlnJYRKKWY9M3jDHNgQIgzVr7pTHmhLW2XbXj/7HW1plDZ4yZAEwAaNWqVd+rrrqqyWoWEfHGli1bvrLWdnC6Dm80JouVwyISDsIli3VPLCKRytscjm2KYhpwK7DVWvtl1eMvjTFJ1tpCY0wScLS+F1lrXwFeAejXr5/dvHlz01QrIuIlY8wBp2vwgc9ZrBwWkXAQRlmse2IRiUje5rCT0zfG8O0wNYC/AeOq/h4HvNnkFYmIRB9lsYiIs5TDIhLVHOmUMMa0BG4BlldrngPcYozZW3VsjhO1iYhEC2WxiIizlMMiIg5N37DWFgOJtdqKqFx5WEREmoCyWETEWcphERFn15QQkRBx9uxZDh8+jMvlcrqUsBIfH09qaipxcXFOlyIiEUBZ3DjKYhEJFOVw4/ibw+qUEBEOHz5MmzZt6NSpE8YYp8sJC9ZaioqKOHz4MJ07d3a6HBGJAMpi3ymLRSSQlMO+C0QOO7nQpYiECJfLRWJiosLXB8YYEhMT1ZMuIgGjLPadslhEAkk57LtA5LA6JUQEQOHbCPo3E5FAU674Tv9mIhJIyhTf+ftvpk4JEQkZc+fOJS0tjfT0dMaMGYPL5cJay4wZM7jyyivp1q0bL774ouf5b7/9Ntdccw1XXXUVGRkZjBo1ioMHDzr4DkREwt/8+fNJT08nLS2NefPmAbB06VLS0tJo1qwZmzdvBmD16tVkZGSQkZFB69at6dq1KxkZGYwdO9bJ8kVEwt6RI0cYPXo0Xbp0oXv37tx222189tlnpKenO11aUGhNCREJCfn5+bz44ot88sknJCQkcNddd7Fo0SKstRw6dIhPP/2UZs2acfToUQB27drFQw89xN/+9je6desGwN/+9jf279/PpZde6uRbEREJW7t27WLBggVs3LiR5s2bk5mZSVZWFunp6SxfvpyJEyd6njt06FCGDh0KwMCBA3n++efp16+fU6WLiEQEay0jRoxg3LhxLFq0CIDt27fz5ZdfOlxZ8KhTQkR8tnJbPs+t3kPBiRKS2yUwbWhXhvdO8fu8ZWVllJSUEBcXR3FxMcnJycycOZO//OUvNGtWObCrY8eOADz77LP8/Oc/93RIANxxxx1+1yAiEi6CkcW7d+/muuuuo2XLlgDcdNNNrFixgunTpweiZBGRiBKMHF63bh1xcXE8+OCDnraMjAz279/vZ7WhS9M3RMQnK7fl8/jyneSfKMEC+SdKeHz5TlZuy/frvCkpKTzyyCNceumlJCUl0bZtW4YMGcLnn3/O4sWL6devH7feeit79+4F4OOPP6ZPnz4BeEciIuEnWFmcnp7Oe++9R1FREcXFxbz11lscOnQoMEWLiESQYOXwrl276Nu3b2CKDBPqlBARnzy3eg8lZ8trtJWcLee51Xv8Ou9//vMf3nzzTb744gsKCgo4ffo0f/7znzlz5gzx8fFs3ryZBx54gPvuu6/Oa4uKisjIyODKK6/k+eef96sOEZFwEKws7tatG48++ii33HILmZmZ9OrVi9hYDawVEaktWDkcjdQpISI+KThR4lO7t9555x06d+5Mhw4diIuL48477+SDDz4gNTWVkSNHAjBixAjy8vIASEtLY+vWrQAkJiayfft2JkyYwKlTp/yqQ0QkHAQriwHuv/9+tm7dynvvvUf79u254oor/D6niEikCVYOp6WlsWXLFr/OEW7UKSEiPklul+BTu7cuvfRSPvzwQ4qLi7HW8u6779KtWzeGDx/O2rVrAfjHP/7BlVdeCcD06dN55pln2L17t+ccxcXFftUgIhIugpXFgGdB4YMHD7J8+XLGjBnj9zlFRCJNsHJ40KBBnDlzhgULFnjaNm3axIEDB/w6byhTp4SI+GTa0K4kxMXUaEuIi2Ha0K5+nffaa6/lhz/8IX369KFHjx5UVFQwYcIEHnvsMf7617/So0cPHn/8cX73u98B0KNHD+bPn8/YsWO56qqrGDBgALt37+bHP/6xX3WIiISDYGUxwMiRI+nevTvDhg3jpZde4oILLmDFihWkpqayYcMGsrKyPLtuiIhEq2DlsDGGFStWsGbNGrp06UJaWhrZ2dkkJyezZ88eUlNTPT9Lly7161qhQpMERcQn7hWFg7H7xhNPPMETTzxRo61Fixbk5ubW+/ysrCyysrL8vq6ISLgJZhb/85//rNM2YsQIRowY0eBr1q9f7/d1RUTCSTBzODk5mSVLltRpP3v2rN/nDkXqlBARnw3vnRKQwBURkcZTFouIOEs5HBiaviEiIiIiIiIijlCnhIiIiIiIiIg4Qp0SIiIiIiIiIuIIdUqIiIiIiIiIiCPUKSEiIiIiIiIijlCnhIhEhYEDB7J582YAbrvtNk6cOOFwRSIioad169Y1Hr/66qtMmjSJ9evX079//xrHysrKuOiiiygsLGzKEkVEIlpMTAwZGRn06tWLPn368MEHHwCwf/9+0tPTHa4uOLQlqIhEnbfeesun55eXlxMTExOkakREQt+NN97I4cOH2b9/P506dQLgnXfeIT09naSkJGeLExGJIAkJCWzfvh2A1atX8/jjj/OPf/zD4aqCSyMlRMR3eUtgbjpkt6v8nbfE71Pu37+fq666inHjxtGzZ09++MMfUlxczJYtW7jpppvo27cvQ4cO9Xwj9+9//5ubb77Z04v8+eefs379em6//XbPOSdNmsSrr75a51qdOnXiq6++AmD48OH07duXtLQ0XnnlFc9zWrduzezZs7n22mvZsGGD3+9PRCTggpDFDWnWrBk/+tGPWLx4sadt0aJFjBkzJmjXFBEJeUHO4a+//poLLrggoOcMRRopISK+yVsCOT+DsyWVj08eqnwM0PMuv069Z88efv/73zNgwADuu+8+XnrpJVasWMGbb75Jhw4dWLx4MTNmzOAPf/gDd999N4899hgjRozA5XJRUVHBoUOHfL7mH/7wB9q3b09JSQlXX301I0eOJDExkdOnT5Oens6TTz7p13sSEQmKIGVxSUkJGRkZnsfHjx/njjvuAGDMmDFMmDCBRx99lDNnzvDWW28xd+7cRl9LRCSsBTmHXS4XhYWFrF27NgDFhjZ1SoiIb9598tvwdTtbUtnuZ6fEJZdcwoABAwD4yU9+wi9/+Ut27drFLbfcAlROo0hKSuKbb74hPz+fESNGABAfH9/oa7744ousWLECgEOHDrF3714SExOJiYlh5MiRfr0fEZGgCVIWVx82DJVrSrjX47n66qs5deoUe/bsYffu3Vx33XVR8Q2eiEi9miCHN2zYwNixY9m1a5c/lYY8dUqIiG9OHvat3QfGmBqP27RpQ1paWp3pE19//XW9r4+NjaWiosLz2OVynfN669ev55133mHDhg20bNmSgQMHel4THx+vdSREJHQFMYvPZfTo0SxatIjdu3dr6oaIRLcmyOH+/fvz1VdfcezYsYCdMxRpTQkR8U3bVN/afXDw4EFPB8Qbb7zBddddx7FjxzxtZ8+e5eOPP+Y73/kOqamprFy5EoAzZ85QXFzMZZddxieffMKZM2c4efIk77777jmvd/LkSS644AJatmzJp59+yocffuj3exARaRJBzOJzGTNmDH/+859Zu3atZ1qHiEhUaoIc/vTTTykvLycxMTFg5wxF6pQQEd8Mng1xCTXb4hIq2/3UrVs3XnvtNXr27Mnx48d56KGHWLZsGY8++ii9evUiIyPDsy3Sn/70J1588UV69uzJ9ddfz5EjR7jkkku466676NmzJ3fffTe9e/c+5/UyMzMpKyujZ8+ezJo1i+uuu87v9yAi0iSCmMXn0r17d1q2bMmgQYNo1apVUK8lIhLSgpTD7jUlMjIyGDVqFK+99ppn9O6ePXtITU31/CxdutSva4UKTd8QEd+458i9+2Tl8LS2qZXh6+d6ElC5uvv//u//1mjLyMjgvffeq/PcK664ot6Ff37961/z61//uk77+vXrPX/v37/f8/ff//73ems5deqUl1WLiDggSFlcO/vuvfde7r333hptO3bs8OsaIiIRIUg5XF5eXm97p06dOHv2rF/nDlXqlBAR3/W8KyCdECIi4gdlsYiIs5TDAaHpGyISEjp16hTxKwuLiIiIiEhN6pQQEREREREREUeoU0JEREREREREHKFOCRERERERERFxhDolRERERERERMQR6pQQkZAxf/580tPTSUtLY968eQAsXbqUtLQ0mjVrxubNmz3PffXVV5k0aZJTpYqIRCxlsYiIs44cOcLo0aPp0qUL3bt357bbbuOzzz7DGMNvfvMbz/MmTZrEq6++6lyhAaJOCREJCbt27WLBggVs3LiRHTt2sGrVKvbu3Ut6ejrLly/nxhtvdLpEEZGIpywWEXGWtZYRI0YwcOBAPv/8cz755BN++ctf8uWXX9KxY0fmz59PaWmp02UGlDolRKLMgXvGcuCesX6dI3dfLkOWDaHnaz0ZsmwIufty/a5r9+7dXHfddbRs2ZLY2FhuuukmVqxYQbdu3ejatWu9rzl06BCZmZl07dqVJ554wtP+1FNPcdVVV3HLLbcwZswYnn/+eb/rExEJJGWxiIizRr28gVEvb/DrHMHI4XXr1hEXF8eDDz7oacvIyOCSSy6hQ4cODB48mNdee83v64SSWKcLEJHwkrsvl+wPsnGVuwAoPF1I9gfZAGRdntXo86anpzNjxgyKiopISEjgrbfeol+/fud8zcaNG9m1axctW7bk6quvJiur8vp//etf2bZtG2VlZfTp04e+ffs2ui4RkVCkLBYRcVawcnjXrl3nzMvHHnuMW2+9lfvuu6/R1wg1jnRKGGPaAb8D0gEL3AfsARYDnYD9wF3W2v84UZ+INGz+1vme8HVzlbuYv3W+XwHcrVs3Hn30UW655RZat25Nr169iI09d0TdcsstJCYmAnDnnXfyr3/9C4Af/OAHJCQkADBs2LBG1xTJlMMi4U1ZHBmUxSLhK1g5fD6dO3fmmmuu4S9/+UvQrtHUnJq+MR9421p7FdAL2A08Brxrrb0CeLfqsYiEmCOnj/jU7ov777+frVu38t5779G+fXuuuOKKcz7fGFPnsbXW7zqihHJYJIwpiyOGslgkTAUrh9PS0tiyZcs5n/Pzn/+cZ599loqKCr+uFSqavFPCGPMd4Ebg9wDW2lJr7QngB4B7csxrwPCmrk1Ezu/iVhf71O6Lo0ePAnDw4EGWL1/OmDFjzvn8NWvWcPz4cUpKSli5ciUDBgzghhtuICcnB5fLxalTp8jN9X9uX6RRDouEP2Vx+FMWi4S3YOXwoEGDOHPmDAsWLPC0bdq0iQMHDngeX3XVVXTv3p1Vq1b5da1Q4cRIicuBY8BCY8w2Y8zvjDGtgIustYUAVb871vdiY8wEY8xmY8zmY8eONV3VIgLA5D6TiY+Jr9EWHxPP5D6T/T73yJEj6d69O8OGDeOll17iggsuYMWKFaSmprJhwwaysrIYOnSo5/k33HAD99xzDxkZGYwcOZJ+/fpx9dVXc8cdd9CrVy/uvPNO+vXrR9u2bf2uLcIoh0XCnLI4IiiLRcJYsHLYGMOKFStYs2YNXbp0IS0tjezsbJKTk2s8b8aMGRw+fNiva4UK09TD64wx/YAPgQHW2o+MMfOBr4GHrLXtqj3vP9baC851rn79+tnqe2WLNJXxb48HYGHmQocr8Z17tffL/vRHT9vu3bvp1q2b1+fI3ZfL/K3zOXL6CBe3upjJfSYHde6cr06dOkXr1q0pLi7mxhtv5JVXXqFPnz5BuVZ9/3bGmC3W2nOvDOcg5bBEgnDOYVAWB5qyWFksznDvXrF4Yn+HK/FdfbVHWg43JX9y2ImFLg8Dh621H1U9XkblXLkvjTFJ1tpCY0wScNSB2kTEC1mXZ4V04E6YMIFPPvkEl8vFuHHjgnYTHMaUwyIRQFkc9pTFImEu1HM4XDR5p4S19ogx5pAxpqu1dg8wGPik6mccMKfq95tNXZuIRIZIWo04GJTDItIUlMXnpiwWEankyJagwEPA68aY5sA+YDyV61ssMcbcDxwEfuRQbSIi0UA5LCLiPGWxiEQ9RzolrLXbgfrmlgxu6lpERKKRclhExHnKYhER50ZKiPgsFBY1cy9Mxt0xjtUgIuIU5bCIiPNCYXFJdw0igeDElqAiUW382+M9N/YiItL0lMMiIs4b9fIGdW4IoE4JkYhwMieHvYMGs7tbd/YOGszJnBynS/KZMYaHH37Y8/j5558nOzu7SWu49957WbZsWZNeU0QiSN4SmJsO2e0qf+ctcboin4VCFouINNbKbfkMmLOWzo/lMmDOWlZuy3e6pEZ55plnSEtLo2fPnmRkZPDRRx+xatUqevfuTa9evejevTsvv/yy02UGjDolJCzk7ssl71gem7/czJBlQ8jdl+t0SSHjZE4OhbNmU1ZQANZSVlBA4azZYdcx0aJFC5YvX85XX33VqNeXlZUFuCIRqU45fB55SyDnZ3DyEGArf+f8LOw6JvzNYhEJrpXb8tl28AQffXE8rD90B8PKbfk8vnwn+SdKsED+iRIeX74z7P6NNmzYwKpVq9i6dSt5eXm88847XHzxxUyYMIGcnBx27NjBtm3bGDhwoNOlBow6JSTk5e7LJfuDbEorSgEoPF1I9gfZPt0QR/JQ3YIZM7EuV40263JRMGNm0K4ZjJEZsbGxTJgwgblz59Y5duDAAQYPHkzPnj0ZPHgwBw8eBCpHNvz3f/833//+93n00UfJzs5m3LhxDBkyhE6dOrF8+XKmT59Ojx49yMzM5OzZswA8+eSTXH311aSnpzNhwgSstX7XLxLJApHDEMFZvDAL3pwEZ0tqtp8tqWwPkqbOYhFxlvtDd2l5BdC4D92RPGVi+rI8Ss6W12grOVvO9GV5QbtmMHK4sLCQCy+8kBYtWgBw4YUX0qZNG8rKykhMTAQqO5C7du3q97VChTolJOTN3zofV3nND92uchfzt853qKIQU1rqW7ufgjky47/+6794/fXXOXnyZI32SZMmMXbsWPLy8rj77rv52c9+5jn22Wef8c477/DCCy8A8Pnnn5Obm8ubb77JT37yE77//e+zc+dOEhISyM3N9Zxv06ZN7Nq1i5KSElatWuV37SKRTDnshfIzvrX7yYksFhFnPbd6T70fup9bvcehikKLu7PG23Z/BSuHhwwZwqFDh7jyyiv56U9/yj/+8Q/at2/PHXfcwWWXXcaYMWN4/fXXqagIzvtygjolJOQdOX3Ep/ZgK3IVhdQQ5tjkZJ/a/XV07rx6R2YcnTvP73N/5zvfYezYsbz44os12jds2MCPf/xjAO655x7+9a9/eY796Ec/Iibm21X4b731VuLi4ujRowfl5eVkZmYC0KNHD/bv3w/AunXruPbaa+nRowdr167l448/9rt2kUimHD6P8bnQ9pL6jzXU7icnslhEnFVwosSn9mD76tSZkJpKktIuwad2fwUrh1u3bs2WLVt45ZVX6NChA6NGjeLVV1/ld7/7He+++y7XXHMNzz//PPfdd59f1wkl6pSQkHdxq4t9ag+mIlcR+0/u93sIcyB1nDoFEx9fo83Ex9Nx6pQ6zz2Zk0Pxjh0Ub9rU6CFmZYWFPrX7asqUKfz+97/n9OnTDT7HGOP5u1WrVjWOuYe6NWvWjLi4OM9zmzVrRllZGS6Xi5/+9KcsW7aMnTt38sADD+Cq9R8UEalJOeyFwbMhrtaNb1xCZXstkZLFItK0khv4cN1QezB9deoMX3x12q+pJIE2bWhXEuJqbhedEBfDtKF1pzkEYm2OYOZwTEwMAwcO5IknnuC3v/0tf/3rX4HKL9mmTp3KmjVrPG2RQJ0SEvIm95lMfEzND93xMfFM7jO5SesY//Z49p/cTwU1h0q5yl3Mfr/uTWdTaTtsGElPPVk5MsIYYpOTSXrqSdoOG1bjee4hZu5pHY0dYhablORTu6/at2/PXXfdxe9//3tP2/XXX8+iRYsAeP3117nhhhsafX53B8SFF17IqVOntNuGiBeUw17oeRcMe7FqZISp/D3sxcr2asI5i0XEWb586A6mUS9vYN+x01TUWpIr2Os3nM/w3in86s4epLRLwFA5QuJXd/ZgeO+UGs8LxNocELwc3rNnD3v37vU83r59OxdddBHr16+v0XbZZZf5dZ1QEut0ASLnk3V5FgCz359NaUUpSa2SmNxnsqe9KdW+EXZzf2PnlLbDhtF22DAO3DMWgBNLlnJiydIazynesaPOOhPuBTHLXvotZ/Z9Uee8LS7vXKet49QpFM6aXWO4WkMjMxrr4Ycf5re//a3n8Ysvvsh9993Hc889R4cOHVi4cGGjz92uXTseeOABevToQadOnbj66qsDUbJIRFMOe6nnXdDzLk8Ws3UVUHPNmnDOYhFxlvvD9fRleZSWV5DSLoFpQ7vW+dDdFBpaIjxY6zd4a3jvFIb3TvEs5vnGxoO8sfFgjedsO3iiTp3uDpXf/SCJ5sdO1Tlvlw6t67QFK4dPnTrFQw89xIkTJ4iNjeW73/0u8+fPZ+LEiUycOJGEhARatWrFq6++6td1Qok6JSQsZF2exbLPKr/RXpjZ+A+k/liYuZC/zcugtKLuomVJrQLzzVRQBWhBTPcIjKNz51FWWEhsUhIdp06pMzLDV6dOffsfgIsuuoji4mLP406dOrF27do6r6kdxtnZ2Q2es/qxp59+mqeffvq85xORbymHAySMs1hEnDe8d4rnQ/biif0dqWHxxP5cOePv9XZABGv9hkAK1IKYwcrhvn378sEHH9Rpf+utt/w6byhTp4SID1LbpLD/5P4abU4MYW7IZX/6Y4PH9g4aXLk6cC2xyclw4YX1fhPXEPfIDBGRphbqOQzKYhGJfJe0T+CLr2pO4XBiKklDztVhM2DOWvLrWRw0pV0CHdq0qHdUREOUw4GhNSVEfJAYn0intp1o3qw5UPnNXPb12Y4MYfaVLwtiioiEqnDOYVAWi0hkuLB1Czpf2IrmMZUfJxtavyEUhcraHPItjZSQsOHUcOHaEuMT6dmhIxA6NXnD3YtbMGMmlJYSm5zsGWJWsHu3w9WJSDgIlcwL1xwGZbGI+M+paRu1Xdi6BRe2rtz1LFRq8sa51ubYvftrh6uLTuqUEIkibYcN8yyAea7hxSIiEjzKYhERZ4XC2hzyLU3fEBERERERERFHqFNCRERERERERByhTgkRCRnz588nPT2dtLQ05s2bB8DSpUtJS0ujWbNmbN68GYDVq1eTkZFBRkYGrVu3pmvXrmRkZDB27FgnyxcRCXtHjhxh9OjRdOnShe7du3Pbbbfx2WefYYxh1qxZnud99dVXxMXFMWnSJKBy2+WUlBRPNj/22GMAlJWV8fOf/5wrrrjCc+yZZ57xnCcmJsbTnpGRwZw5cwAYOHAg/fr18zxv8+bNDBw4sAn+BUREnBdtWaw1JUQkJOzatYsFCxawceNGmjdvTmZmJllZWaSnp7N8+XImTpzoee7QoUMZOnQoUBmWzz//fI3AFBER31lrGTFiBOPGjWPRokUAbN++nS+//JLLL7+cVatW8dRTTwHfdhhXN3XqVB555JEabTNnzuTIkSPs3LmT+Ph4vvnmG1544QXP8YSEBLZv315vPUePHuXvf/87t956ayDfpohISIvGLNZICRHx2WcfHeG1n7/PSw+u5bWfv89nHx3x+5y7d+/muuuuo2XLlsTGxnLTTTexYsUKunXrRteu2qJJRKS2L1atYuXNN/OX9HRW3nwzX6xa5df51q1bR1xcHA8++KCnLSMjg0suuYSEhAS6devmGbG2ePFi7rrrrnOer7i4mAULFvCb3/yG+KptUNu0aUN2drZX9UybNo2nn366cW9GRKQJBDqHITqzWJ0SIuKTzz46wrrXP+XU8TMAnDp+hnWvf+p3x0R6ejrvvfceRUVFFBcX89Zbb3Ho0KFAlCwiEnG+WLWKjb/4BcWFhWAtxYWFbPzFL/y6Id61axd9+/Zt8Pjo0aNZtGgRhw8fJiYmhuTk5BrH586d6xn6u3r1av79739z6aWX0qZNmwbPWVJSUmPI8OLFiz3H+vfvT4sWLVi3bl2j35OISLAEI4chOrNYnRIi4pMNb35OWWlFjbay0go2vPm5X+ft1q0bjz76KLfccguZmZn06tWL2FjNMBMRqc+OefMod7lqtJW7XOyoWo8nGDIzM1mzZg1vvPEGo0aNqnN86tSpbN++ne3bt3um2FW3cOFCz7d97k5n95Bh90/t886cOVOjJUQkJDmRwxCZWaxOCRHxiXuEhLftvrj//vvZunUr7733Hu3bt+eKK67w+5wiIpGo+Ej9o9MaavdGWloaW7ZsafB48+bN6du3Ly+88AIjR4487/m++93vcvDgQb755hsAxo8fz/bt22nbti3l5eVe1TRo0CBcLhcffvihd29CRKSJBCOHITqzWJ0SIuKT1u1b+NTui6NHjwJw8OBBli9fzpgxY/w+p4hyfm3vAAAgAElEQVRIJGp58cU+tXtj0KBBnDlzhgULFnjaNm3axIEDBzyPH374YZ599lkSExPPX2PLltx///1MmjQJV9W3ieXl5ZSWlvpU14wZM/j1r3/t02tERIItGDkM0ZnF6pQQaYSFmQtZmLnQ6TIc0f8HXYhtXjM6Yps3o/8Puvh97pEjR9K9e3eGDRvGSy+9xAUXXMCKFStITU1lw4YNZGVl1TsMTUSiTzTnMECvKVOIqVqwzC0mPp5eU6Y0+pzGGFasWMGaNWvo0qULaWlpZGdn15ivnJaWxrhx47w+5zPPPENSUhLp6en07t2b733ve4wbN85zztrzmN3b11V322230aFDh0a/LxEJnsUT+7N4Yn+ny3BEMHIYojOLNWFbxEsnc3Io3rEDSkvZO2gwHadOoe2wYU6X1eSuvLay93fDm59z6vgZWrdvQf8fdPG0++Of//xnnbYRI0YwYsSIBl+zfv16v68rIuFBOfytzrffDlTOaS4+coSWF19MrylTPO2NlZyczJIlS+q079q1q07bvffey7333gvQ4CrucXFxzJkzx7PnfW0NDR2une3nGsosIk1r5bZ8th08QWl5BQPmrGXa0K4M753idFlNLlg5DNGXxeqUEPHCyZwcCmfNhqphTmUFBZWPISpviK+89uKAdEKIiHhLOVxX59tvD8jNr4iIt1Zuy+fx5TspLa9c9Dz/RAmPL98JELUdE8ph/6lTQiLagXvGAjD6+KeVj18f26jzuL+Zq866XBTMmMmJJUt9Old9tVz2pz82qi4RkVDnzmHwL4uVwyIijTfq5Q0AfFL4dY3HvnKPkKiu5Gw505fl8cbGgz6dq75aonUqSLTTmhIi3mhoIRgfF4gREZFGUg6LiDiudofE+dpFvKGREhLR3N98Zb89HqDRi6LtHTSYsoKCOu2xyck+fbuWuy+X2e/PprSilKRWR5jcZzJZl2c1qiYRkXBQPSP9yeJoyWFrLcYYp8sIK9Zap0sQCXnuEQjuUQmNHZEwYM5a8k+U1GlPaZfg0zlXbstn+rI8SssrOPyfkpBal0I57Dt/c1gjJUS80HHqFEyt1XVNfDwdp3q/um7uvlyyP8imtKLyW73C04Vkf5BN7r7cgNYqIhKJoiGH4+PjKSoq0odsH1hrKSoqIr7W/zZEJDimDe1KQlxMjbaEuBimDe3q9TkaWpdi5bb8gNbaGMph3wUihzVSQsQL7kXUCmbMhNJSYpOTvVr1fXzVt4IAecfyPDfCbq5yF7Pfn82yz5Y12dZ2mjctIuGoKXIYGj+izlf1ZXFqaiqHDx/m2LFjTVJDpIiPjyc1NdXpMkSigns0g3uUQ0q7BK9GOVRfN+J861I01boS9V1HOdw4/uawOiVEvNR22DDPYmqN+WBf+0b4fO3R6MiRI0yZMoVNmzbRokULOnXqxLx58+jatSszZ87kqaeeAuCrr74iKSmJiRMn8tvf/pbs7GwWLFjg2Ts5MzOTOXPmUFZWxuzZs1m6dCmtWrUC4Ec/+hEzZswAICYmhh49eniuP3r0aB577DEGDhzIqVOn2Lx5MwCbN2/mkUce0fajIg6L9ByOi4ujc+fOTpchInJOw3uneBa1bEwHQiivS6EcdoY6JSTquFeCb4oRA9W/cRuybAiFpwvrPCepVVKTfTMXyqy1jBgxgnHjxrFo0SIAtm/fzpdffsnll1/OqlWrPJ0SS5cuJS0trcbrp06dyiOPPFKjbebMmRw5coSdO3cSHx/PN998wwsvvOA5npCQwPbt2+ut5+jRo/z973/n1ltvDeTbFBGUwyIiocDf9SV8Uf0agVqXQiKH1pQQaSKT+0wmPqbmXKv4mHgm95nsUEWN98WqVay8+Wb+kp7Oyptv5otVq/w+57p164iLi+PBBx/0tGVkZHDJJZeQkJBAt27dPCMXFi9ezF133XXO8xUXF7NgwQJ+85vfeOa4tWnThuzsbK/qmTZtGk8//XTj3oyIhKRIymERkXAViHUpJLKoU0KkiWRdnkX29dk0b9YcqPxmLvv67JBa9d0bX6xaxcZf/ILiwkKwluLCQjb+4hd+d0zs2rWLvn37Nnh89OjRLFq0iMOHDxMTE0NycnKN43PnziUjI4OMjAxWr17Nv//9by699FLatGnT4DlLSko8r8nIyGDx4sWeY/3796dFixasW7fOr/clIqEjUnJYRCScDe+dwq/u7EHzmMqPointEvjVnT1CZvcNaXqaviHShLIuz2ryxdQCbce8eZS7XDXayl0udsybR+fbbw/adTMzM5k1axYXXXQRo0aNqnO89vSNvLy8GscXLlzI/PnzKSoq4oMPPvCMwGho+gZUTv94+umnefbZZwP3RkTEUZGQwyIi4c7fdSkksjgyUsIYs98Ys9MYs90Ys7mqrb0xZo0xZm/V7wucqE0kVBy4Z6xn3nUoKT5yxKd2b6WlpbFly5YGjzdv3py+ffvywgsvMHLkyPOe77vf/S4HDx7km2++AWD8+PFs376dtm3bUl5e7lVNgwYNwuVy8eGHH3r3JsKMsljk/EI1iyUyKIdFzm/Uyxtq7N4hkcfJ6Rvft9ZmWGv7VT1+DHjXWnsF8G7VY5EmM/7t8TW2jpP6tbz4Yp/avTVo0CDOnDnDggULPG2bNm3iwIEDnscPP/wwzz77LImJieevs2VL7r//fiZNmoSramRHeXk5paW+rbI/Y8YMfv3rX/v0mjCjLJaQoRyWKKUclpCiTgBpaqG0psQPgNeq/n4NGO5gLSLSgF5TphATX3OhuJj4eHpNmeLXeY0xrFixgjVr1tClSxfS0tLIzs6usXZEWloa48aN8/qczzzzDElJSaSnp9O7d2++973vMW7cOM85a68p8dhjde/7brvtNs9Wo1FCWSwi4izlsIhEFafWlLDA/xljLPCytfYV4CJrbSGAtbbQGNPRodpE5Bzc60bsmDeP4iNHaHnxxfSaMiUg60kkJyezZMmSOu27du2q03bvvfdy7733AjS4o0ZcXBxz5sxhzpw59R5vaBrH+vXrazw+17SSMKcsFhFxlnJYRKKeU50SA6y1BVUhu8YY86m3LzTGTAAmAFx66aXBqk9EzqHz7bcHdVFLaTKNymLlsIhIwOieWESiniPTN6y1BVW/jwIrgGuAL40xSQBVv4828NpXrLX9rLX9omxItTRAc5BFGqexWawclvooi0V8p3tiCSStBSHhqsk7JYwxrYwxbdx/A0OAXcDfAPdk8XHAm01dm0Sv3H255B3LY/OXmxmybAi5+3KdLkkkqJTFEmqUwxJtlMMSilZuy2fbwRN89MVxBsxZy8pt+U6XJFHAiekbFwErjDHu6//FWvu2MWYTsMQYcz9wEPiRA7VJFMrdl0v2B9mUVlTuylB4upDsD7KByv3so4W1lqr/X4qXrLVOl+APZbGEDOWwRCnlsISUldvyeXz5TkrLKwDIP1HC48t3AjC8d4qTpUmEa/JOCWvtPqBXPe1FwOCmrkcin/vbt9KKUoYsG8LPvvqKuGaxZFcNM3Yfq85V7mL2+7NZ9tkyT9vCzIVNWndTio+Pp6ioiMTERHVMeMlaS1FREfG1diIJF8piaWrVs3jLlxU0j2nhcw4DZDdVwSJBphyWpuYeBVFaXsGAOWspOn2GuJhmnikf7mPVlZwtZ/qyPN7YeNDTtnhi/yatWyKfUwtdijSJ+r59c5WXA99+kKx9I3y+9kiUmprK4cOHOXbsmNOlhJX4+HhSU1OdLkMk5NXO4goqcJW7KHIVkRifqBwWEQmy+kZB1Fa7Q+J87SKBok4JiWjzt87HVe6q1WqpsBWekQ9Dlg2h8HRhndcmtUqK6NER1cXFxdG5c2enyxCRCNVQFrvKXCzMXOhTDh94fWwQKxURiUzPrd5Dydm6W6FXVHw78mHAnLX1dlaktEvQ6AgJKkd23xAJlPMtjHbk9JF6X1f927fJfSYTH1NzCH58TDyT+0wOfMFeOpmTQ/GOHRRv2sTeQYM5mZPjWC0iIufT2Cx2t4diDoOyWETCx/kWqCyop7MBao6CmDa0KwlxMTWOJ8TFMG1o18AX7CUtvBkd1CkhYauhhdGq3wxf3Oriel/bvFlzz99Zl2eRfX22py2pVRLZ12c7trjayZwcCmfNhtLK91VWUEDhrNm6GRaRkORPFrvbQy2HQVksIuGjoQUqq3+AT26XUO9rm8d8+3FweO8UfnVnD09bSrsEfnVnD8cWufTmfUlkUKdEhIjG/eFnvz+7znBg98JobvV9+9aMZqS2qRmuWZdn0bNDT/pd1I//++H/OXYjfOCesRTMmIl11Xxf1uWiYMZMR2oSEe9EYw6Pf3t8g1k8f+t8z+OGsrj6SIhQyWE4dxYfnTvPoapExBujXt7gWbgxWkxflldnaoZ7gUq3+kZBNDNwSfuanRXDe6fQ+9J2XNu5Pe8/NsixDolRL2/w6n1JZNCaEhK2vFkYzX1TO/v92ZRWlJLUKolObeNIjE9skhobpbSBhd0aahcRcVBDWVx9ykbtLG7erAWpbVLoE8rbfTaQuWWFdde+EBFxkjcLVLo7F6Yvy6O0vIKUdgnExzXjwtYtmqTGxtDCm9FDIyUkbFWfglFdUqukGo9rf/sWyh0Sl/3pj8QmJ9d7rKF2ERGnLMxc2GAW156yUT2Le3XoGb5ZnJRUb7uIiFOqT8GoLqXduUdBhHKHxOKJ/evU79ZQu4QvdUpI2Eppk4LB1GjzZ2G0hZkLQ2K3jY5Tp2Diaw5zNvHxdJw6xaGKREQaFsgsDpUcBmWxiISPS9on0KxmDPu1QOXiif1DYreNUFx4U4JD0zckbLm/Zcv/Jt8zNWNyn8lBnYf86fFPAbgsaFeAtsOGAVSuIVFaSmxyMh2nTvG0i4iEEqeyOPvt8UHtwFAWi0i4cI94OHS8xDM1Y9rQrkFdD+KTwq+Ddm63+qacBPt9iTPUKSFhLTE+0XND7M3NqXt7N0pL2TtosE83mLn7cjlz9jQVVDBk2ZCg3nS3HTaME0uWApXDiEVEQpkvWVxWVETxjgONymGAIlcRp8+e9mw/qiwWEansmHB3TngzysG91WZpeQUD5qz16cP+ym35nHKVYcHn1/pqeO8U3th4EPDufUl4UqeERI2yoiIKn6u7vRtw3hti95Z3rrubAc2gass7wNEV4kVEwklZURGlX3wBpWWVj33IYajK4h8cxVVeOfu0UFksIuKzr06dqXerTeC8nQvubTpt1WNfXivSEHVKSMQ5cM/YOm2jj39KxenTWFfN1XrdW226vwlzq/6N2Pi3x5N3LK/OCvPu7Ud1IywiUld9WfzjQ59Dhe85DMpiERFf1bc16ieFX3tGOVTn3mrTPSrBrfrohFEvb/CMrqjvteqUkMbSQpcSPSoa2D7Ii602vdl+VEREzsOPHAZlsYhIINTukHDzZqtNbdMpwaCREhL2as9frm/eb/bb4ynZsYNfLCypcyw2Ofmcc4UXZi5kyLIhFJ6uuzd97e1HRUSilTdZ/OizfbBnSsn+S3mN9vPlsPv8ymIRkYbVXnOhvjUYGhrtAJVbbZ5r3YbFE/szYM5a8k/UvZ/WNp3iD42UkKgRl5ra6O3dJveZTHxMzdf6s/2oiEg0iktNpfa+db5ss6ksFhHx3yXtExq91aa26ZRgUKeERI3YxESSnnoSmjevfJycTNJTT3q1uFrW5VlkX59N82aVr01qlUT29dmawywi4oPYxESad+7cqBwGZbGISCBc2LoFv7qzB81jKj8KprRL4Fd39vBqTYjhvVMa/VqRhmj6hpzT+LfHA95ttxkO/NneLevyLN34iogjIimLYxMTadmrHdC4bTaVxSLiBPeikZGyLaU/W20O752iTggJKK87JYwx7b14WoW19oQf9UiEy92Xy+z3Z1NaUUpSq6Sg7i8vIiJ1KYdFRJy3cls+05flUVpeQUq7BKYN7aoP+hK1fBkpUVD1Y87xnBjgUr8qkoiVuy+X7A+yPaukR+v+8pHwTaeIhCflcCXlsIg4aeW2fB5fvtOz2GT+iRIeX74TIKo6JiJl1In4z5c1JXZbay+31nZu6AcoClahEv7mb52Pq9xVo81V7mL+1vkOVSQiEl2UwyIizntu9R5KztbchajkbDnPrd7jUEUizvJlpIQ3XVnq7pIGHTl9xKd2p5zMyaFgxkwoLSU2OZmOU6d4vQhbIDVmrrW3DtwzNujXEJHQoxz2XbByUjksEr0K6tlS81ztTgmVKSbBHFERaWuFhCuvR0pYa10NHTPGtDvfc0QubnWxT+1OOJmTQ+Gs2VBaObS5rKCAwlmzOZmT43BlIiL+Uw6LiDgvuV2CT+1OaGiKycpt+Q5XJpHovJ0Sxpi+xphfGGMuMMZ8xxhznTHmfmPM/xhjVhtj8oEvmqBWaWK5+3LJO5bH5i83M2TZEHL35fp1vnDYX/7o3HlYV82+NetycXTuPIcqEpFoF8gsVg6LiPhu5bZ8th08wUdfHGfAnLV+fzCfNrQrCXExNdoS4mKYNrSrX+cNJE0xkabkzfSNl4GJwEHgG+Bj4FNgNzAayLDWHg1aheKIYCyG5n5dKK/6XlZY6FO7iEgwBTqLlcMiIr4JxqKU7teFwtSIhoTLFBOJDN50SnwATAO2Ai2BBdbaJQDGmGnqkAgP7j3uvZV3LM9zE+zmKncx+/3ZLPtsmdfnqb3CeajvLx+blERZQUG97SIi/gqFLFYOi0g0c68h4K1tB094OiTcSs6WM31ZHm9sPOj1eWqvWTC8d0pIdULUltwugfx6OiBCaYqJRI7zTt+w1v4MuM9aexMwFLjOGLPBGHMrYINdoDij9k3w+dojRcepUzDxNYc2m/h4Ok6d4lBFIhLNojGLlcMiEkpqd0icrz1ShMMUE4kcXu2+Ya0trvp9HPhvY8xlwNPARcaYgdba9cErUQLB1z3ZhywbQuHpukNlk1olRfT+7u7V3Y/OnUdZYSGxSUmOrvouIpFFWXx+ymERCSZfd1kYMGdtvSMGUtolRPSODe5RHM+t3kPBiRKSQ3CKiUQOX7YE9bDWHgDuMca8AMwxxjxprb0xsKWJkyb3mUz2B9k19rMPtcXQgqXtsGG6+RWRkBCtWawcFpFQMW1oVx5fvrPGoo/RMmIg1KeYSOTwekvQ+lhrt1trM4FfBKgeCRFZl2eRfX02zZs1Byq/lcu+Pjtk5iEHemcQEZFQpCwWEXHW8N4p/OrOHjSPqfzYlNIugV/d2SNkPqwHemcQESd4PVLCGLPVWtunvmPW2nXne46En6zLszwLqYXSMOFg7AwSTU7m5FC8YweUlrJ30GANixYJccriyKMcFgkvw3uneBa1DKUpG8HYGSSauDt0SssrGDBnraanOMiX6RvdjDF55zhugLZ+1iMCnHuF+kCtRh+NTubkUDhrNpRW/vuVFRRUPgbdEItIDefbKURZ3DjKYRHxxbl2CwnUziDRSB06ocWXTomrvHhO+fmfIuKfaFyN3l8H7hkL4PlmrjrrclEwYyYnlizlsj/90YnyRCQMKYt9d+CescphEQmYaN0ZxB/uTp7zdeiE0oiYaOB1p0TV4pYiTeJcQ5Qbsxr9+b7xixqlDXxYaKhdRKLW+aaKKIsbSTksIj4414fjxuwMcq6RF9FEHTqhxa+FLkWcMLnPZOJjau5hHw2r0fvjsj/9kcv+9Edik5PrPR6bnKxv50TEJ8pi3ymHRSSQpg3tSkJcTI22aNkZpLEWT+zP4on9SWmXUO/xSN/qNVQ1aktQiR6htKiam3sBtdnvz6a0opSkVklM7jNZC6t5oePUKRTOmo11fbu9oImPp+PUKQ5WJSLnoyyOHMphkfAUih9U3WsfTF+WR2l5BSntErRYo5eieavXUNToTgljTCvAZa3VOhIOc2/JVlpRypBlQ6LipjBUV6MPde5F1ApmzITSUmKTk7Xqu0gARGMOg7K4MZTDIsETjbsphOrOIKFOHTqhxZctQZsBo4G7gauBM0ALY8wx4C3gFWvt3qBUKQ3Slmziq7bDhnFiyVIADRUWCQDlsPhKOSwSeNpNQXylDp3Q4ctIiXXAO8DjwC5rbQWAMaY98H1gjjFmhbX2z4EvM/IEarGvYG7JFqnfeukGUETcApHFwd4aMxKzWDksIm6BWngxmNtjRuoH1kh9XxJ+fFno8mZr7VPW2jx3hwSAtfa4tfav1tqRwGJvT2aMiTHGbDPGrKp63NkY85ExZq8xZrExprkPtUUtbckmIo2lHA4M5bCI+ENZHBjaTUEkfPmyJehZY8xVQArwkbX2lPuYMSbTWvu2tfasD9eeDOwGvlP1+FlgrrV2kTHmf4H7gf/Ph/OFlUB989WYLdlERKpEdQ5DYLJYOSwiforqLA7Ut/WN2R5TREKD1yMljDE/A94EHgJ2GWN+UO3wL325qDEmFcgCflf12ACDAPc419eA4b6cM1ppSzYRaQzlcOAoh0WksZTFgaPtMUXCly9rSjwA9LXWnjLGdAKWGWM6WWvnA8bH684DpgNtqh4nAiestWVVjw9TOSKjDmPMBGACwKWXXurjZSOPtmQTkUZSDgeIclhE/KAsDhDtpiASvnzplIhxT9mw1u43xgyksmPiMnzolDDG3A4ctdZuqToHDbze1vd6a+0rwCsA/fr1q/c50UZbsomIL5TDgaccFhFfKYsDT7spiIQnXzoljhhjMqy12wGqRkzcDvwB6OHDeQYAdxhjbgPiqZw/Nw9oZ4yJreoZTgUKfDinSEAduGcsoBXiJWIphyUsKIslwimLJeS5d0dRJ48Eky+dEmOBsuoNVWE51hjzsrcnsdY+TuW2olT1Cj9irb3bGLMU+CGwCBhH5foVIhIEusGPbsphEecph0VZLOI8dbaEBq8XurTWHrbWHmng2PsBqOVR4L+NMf+mcj7d7wNwThER8Z5yWETEecpiEYkqXo+UqNp9Y7m19nCgLm6tXQ+sr/p7H3BNoM4tIiLnpxwWEXGeslhEopkv0zeeAh4zxnwOvAEstdYeC05ZIoGlhedERJynLBYRcZamK0go8nr6BrCPysV2ngL6Ap8YY942xowzxrQ590tFRERERERERGrypVPCWmsrrLX/Z629H0gG/h+QSWWHhYiIiIiIiIiI13yZvlFj32Rr7Vngb8DfjDEJAa1KRERERERERCKeLyMlRjV0wFpbEoBaRERERERERCSK+LIl6GfGmC7GmEeMMfONMS8YYx40xrQNZoEiTe1kTg7FO3ZQvGkTewcN5mROjtMliYhEHWWxiIizVm7LZ9vBE3z0xXEGzFnLym35TpckEcrrTomqLUH/F4gHrgYSgEuADcaYgUGpTqSJnczJoXDWbCgtBaCsoIDCWbN1Mywi0oSUxSIizlq5LZ/Hl++ktLwCgPwTJTy+fKc6JiQofFlT4gEgw1pbboz5H+Ata+1AY8zLwJtA76BUKBJAB+4Ze87jxTt2eG6C3azLRcGMmZxYsrTB1132pz8GpD4RkWigLBYRcdaolzec8/i2gyc8HRJuJWfLmb4sjzc2HmzwddpyVBrDlzUl4NtOjBZAGwBr7UEgLpBFiTim1k3wedtFRCTwlMUiIo6q3SFxvnYRf/gyUuJ3wCZjzIfAjcCzAMaYDsDxINQmEnDn+xZt76DBlBUU1GmPTU7WN3AiIgGiLBYRcdb5RjQMmLOW/BN19zJIaZeg0RAScF53Slhr5xtj3gG6Af9jrf20qv0YlZ0UIk1qYebCgJ+z49QpFM6ajXW5PG0mPp6OU6cE/FoiIpFAWSwi4qxgdBJMG9qVx5fvpORsuactIS6GaUO7BvxaIr6MlMBa+zHwMYAx5gbgGmCXtfb/glCbSJNrO2wYAAUzZkJpKbHJyXScOsXTLiIiwacsFhFx1vDeKQBMX5ZHaXkFKe0SmDa0q6ddJJC87pQwxmy01l5T9fcDwH8BK4BfGGP6WGvnBKlGkSbVdtgwz0JqGiYsIuIMZbGIiLOG907xLGqpKRsSTL4sdFl9McsJwC3W2ieAIcDdAa1KRERERERERCKeL9M3mhljLqCyI8NUrSWBtfa0MaYsKNWJiIiIiIiISMTypVOiLbAFMIA1xlxsrT1ijGld1SYiIiIiIiIi4jVfdt/o1MChCmBEQKoRERERERERkajh0+4b9bHWFgNfBKAWEREREREREYkivix0KSIiIiIiIiISMOqUkJB14J6xHLhnrNNliIhELeWwiIjzRr28gVEvb3C6DJGg8Xv6hkgkuuxPf3S6BBGRqKcsFhFx1uKJ/Z0uQaKARkqIiIiIiIiIiCO87pQwxnxjjPm6np9vjDFfB7NIiT4nc3Io3rGD4k2b2DtoMCdzcoJ+vb2DBrO7W/cmuZ6ISKhr6hx2X1NZLCLyrZXb8tl28AQffXGcAXPWsnJbftCvN2DOWjo/ltsk1xMB37YEbRPMQkTcTubkUDhrNpSWAlBWUFD5GGg7bFjQrmddria5nohIqGvqHK5+TWWxiEilldvyeXz5TkrLKwDIP1HC48t3AjC8d0rQrldytrxJrifipjUlJOQcnTvPc1PqZl0ujs6d1+gb03Mt1Fa8Y4fnxrv69QpmzOTEkqX1vkbznEUkkgUjh0FZLCLii+dW7/F0ELiVnC3nudV7Gt1JcK4FM7cdPOHpAKl+venL8nhj48F6X6M1JyQQvO6UMMZ8A1jAVGt2P7bW2u8EuDbxwcLMhU6XEDBlhYU+tfut1k3wedtFROqhHPaTslhEAiCSPiQXnCjxqd1ftTskztcuEiiaviEhJzYpibKCgnrbG+tc36btHTS4/uslJ+tbOBGJSsHIYVAWi4j4IrldAvn1dEAkt0to9DnP1WkzYM7aeq+X0i4hojp7JPRo9w0JOR2nTsHEx9doM/HxdJw6JSKuJyIS6pzIRWWxiEhN04Z2JSEupkZbQlwM04Z2jYjribhpTQkJOe75ygUzZkJpKbHJyXScOsvLX5YAACAASURBVCVoC525z3t07jzKCguJTUoK6vVEREJdU+dw9Wsqi0VEKrnXjZi+LI/S8gpS2iUwbWjXoC066T7vc6v3UHCihOQgX0/EzVhrna6h0fr162c3b97sdBkSJO4F0TRsV8KNMWaLtbaf03U0BeVwZFMOSzhTFkukcC9OqSkUEm68zWFfFrr8l7X2hmoLXnoOoYUuRURERERERMRHvix0eUPVby14KU1C38yJiDhLOSwi4jyNkJBIp4UuRURERERERMQRXndKGGOuNsZcXO3xWGPMm8aYF40x7YNTnoiIiIiIiIhEKl9GSrwMlAIYY24E5gB/BE4CrwS+NBERERERERGJZL5sCRpjrT1e9fco4BVr7V+Bvxpjtge+NBERERERERGJZL6MlIgxxrg7MQYDa6sd86VzQ0RERERERETEp86EN4B/GGO+AkqAfwIYY75L5RQOkbB04J6xgFaZFxFxzMKsyt/jc52tQ0QkmimLxSFej5Sw1j4DPAy8CtxgrbXVzvGQt+cxxsQbYzYaY3YYYz42xjxR1d7ZGPORMWavMWaxMaa5929DJPBO5uSwd9Bgdnfrzt5BgzmZk+N0SSIBoyyWsJC3BOamQ3a7yt95S5yuSCRglMMSNpTFEmS+7L5hrLUfWmtXWGtPu9uttZ9Za7e6n+PFqc4Ag6y1vYAMINMYcx3wLDDXWnsF8B/gfl/eiEhjnMzJoXjHDoo3barR8XAyJ4fCWbMpKygAaykrKKBw1mx1TEgkURZLaMhbAoc3wYF/1bzZzVsCOT+Dk4cAW/k752e6GZZIohyW0KEsFgf5sqbEOmPMQ8aYS6s3GmOaG2MGGWNeA8ad7yS20qmqh3FVPxYYBCyran8NGO5DbSI+c3c8UFoKUKPj4ejceViXq8bzrcvF0bnznChVJOCUxRIS3De75WcqH1e/2X33SThbUvP5Z0sq20UigHJYQoayWBzmS6dEJlAOvGGMKTDGfGKM2QfsBcZQ2aP7qjcnMsbEVO3YcRRYA3wOnLDWllU95TCQ0sBrJxhjNhtjNh87dsyH8kVqKpgxs8GOh7LCwnpf01C7SDhqbBYrhyVg3pzU8M3uycP1v6ahdpEwpHtiCQnKYnGYL2tKuKy1/89aOwC4jModOPpYay+z1j5grfV6W1Brbbm1NgNIBa4ButX3tAZe+4q1tp+1tl+HDh28vaRIXVUjJGorKywkNimp3mMNtYuEo8ZmsXJYAsb9rVxtJw9D29T6jzXULhKGdE8sIUFZLA7zZaSEh7X2rLW20Fp7wp+LV71+PXAd0K7alqOpQIE/55bwceCesey5+hrPLhhNJTY5uf72pCQ6Tp2CiY+v0W7i4+k4dUqjr3fgnrFN/h5FvKEsFnAoi9te0kB7KgyeDXEJNdvjEirb/bEw69sV5kVChHJYAMa/PZ7+f+nP+LfHN+2FmzqLlcNSS6M6JfxhjOlgjGlX9XcCcDOwG1gH/LDqaeOAN5u6NnFGWVERFadO1VlsMtjO1fHQdtgwkp56srLjwhhik5NJeupJ2g4b1iS1iQSbslhqcySLz3Wz2/MuGPZi1c2yqfw97MXKdpEIoByW2opcRZw+e5rNX25myLIh5O5roq05lcXisNjzPyXgkoDXjDExVHaKLLHWrjLGfAIsMsY8DWwDfu9AbdLETubkUPrFF1C1w6x7sUkg6B0A7vMXzPj/27v38Liu+t7/76WbZctGxg6OR7Hj2DSE5AQ7cR0ggV84takJP1kN8AuGnj4JNZwnwK8pkuFAE5woinMeDj2cHtvpaYH8aN0k5RJHzQXZLaR1aHlyITcHK4FcDL4ktuQ4cSIRW5ZHl/X7Y8+WZkZ7j+ayZ/aemc/refRYs+a2PLY+WrNmre+6EeJx6lpaJiYk3Os1CSEVTFksE0LLYndQ+8B1zvLh5sWTg2D3eg18pXIph2XCrv27ODh4EJvYrdN/sp+uR7sAaF1W5FUFymIJWdaTEsaY9wHPW2t/l5jNvR5YCfwa+Ia1djCbx7HW9gIXe7Tvx9lLJ1Xk2JatE4Ngl1tsshQTAs1tbQzsuAeAJXfdWfTnE4kKZbEkCzWLl6+Hp+9wvt9Qok8FRSJAOSzJtu3ZNjEh4RoeG2bbnm3Fn5QAZbGEKpeVEn8PrEh8vw0YwjlHeQ2wHfhEsF2TahCFUy40GSEi1S70LNYAWESq3NGTR3NqLwplsYQkl5oSNUnHE62y1nZYax+21t4CLCtC36QK6JQLx2BPD/tWr+H58y8oaV0NERFQFk/o3QFbLoSuuc6fvTvC7pGIVImFTQtzaq9YyuGqlMukxHPGGLcU7F5jzCoAY8y7gJHAeyZVYcHGDjAmpa3QUy7KzWBPD/03dTLa1wfWTuzl1sSEiJSKshhn4NvzJRh8BbDOnz1f0oBYREqifWU7htQcbqxtpH1le0g9CoFyuGrlMinxX4EPGWN+C1wAPGaM2Q/8f4nrRHz5HYfZ3NZGw9KlE4Phajzl4tiWrdjh4ZQ2dy+3iEiQlMUZ7N4MI6dS20ZOOe0iIgHZ8JMNnkd+ti5r5ZzmcyYmJmJNMbou6ypNPYmoUA5XrawnJay1g9baPwUuAq7FOUf5Umvth6y1e4vUP6kAgz09DO3d63nM3GBPD/HDh8HaKadfVIvQ93KLSFVQFk9j8HBu7SIiOdq1fxe9r/V6Hvm5a/8ujrx1BIsl1hSjfWV7dU1IgHK4iuVy+kYH8DDwS01CSLbcrQnE40DqMXOA73XVNBiui8WcrRse7SIiQVAWZ6F5UWLJsEe7iEiBdu3fRdejXcTHnaxNPvIT8L2uqiYmlMNVK5ftG4uA24Bjxph/N8Z8wxjTaoyZV6S+SQXo23Sj79YEbVtwaC+3iBSbsjgL9TPB1ExtW9PpfXsRkRx0PtLJ8Fhq1rpHfm7bs833uqqyplM5XKVy2b7x36y1lwELga8DbwCfxSmA+esi9U/KXeKTt3Sj/f1lvW3Bb192Pib2cjc0gDHVuZdbRIqrQrOY7a3OVxCaFsD834PmxYBx/my7DZavD+bxRaSquasg0h09eTQax4HmK8gcXr7eyeHaGSiHq0vW2zeSzATeBjQnvvqAZ4PslFSOupaWjFsTqmrbwlH/H5O6+fOpmz+fJXfdWcIOiUi1UBYnbG+Fo72wcLn39U0LYMMu7+tERAoQa4rRf3LqZK975Gem6yrO0V7/65oWKIurUNYrJYwxtxtjHgHuBi4FHgU+aa1dZa2dWkJWii7IT+uD4NWfBRs7MI2NKW3u1oRM11WTTMXnRCT6opTFfn1RFk+jdwccfhIOPQxbLtTxcyJlxu9Ei7B49ad9ZTuNtalZ6x75mem6qqIsrlq5rJQ4G5gB7AOOAIeBgWJ0SiqHuwXh2JatjPb3UxeLTanq3rfpRojHK7ri+2BPD2MnToC17Fu9ZuLvman4XCW+DiISDmVxwsljcPqtyQGvu0+550swdtr5fvAV5zJoybCIBMYtWLltzzaOnjzKwqaFU07Y6Hykk/h4vLJP3+jd4eQwdjKHl6932pXFVSvrSQlr7RXGGAP8J+Ay4CvAhcaYN4DHrLU3F6mPUgbcT/uJx1PedIMzGPYb3Da3tTGw4x6Ait264E48NDRawJl46PvaX/D6t7/jHMGXttfbDg/Tt+nGynxDICJFkymHQVlM7w44/hvAyWIGX4H7Pu8UVRsfTb3tyCl44DoNhEUkZ+6xn/HxOGu716ZMLrQua/WdaGhd1kr3S90AbL9ie8n6W1LuxEN6Dv/8W/DmockJCZeyuGrkcvoG1vEc8M/AvwCPAO8EqmxtkSTz+7S/Urch5Lrdwq16H3+rlvhbtU6jtZ4TEhP82svIgZ07uf/DH+YHF17I/R/+MAd27gy7SyIVq9pyGMh9me8D14EdT22z41MnJFzpg+MypBwWKS2/Yz937a/Q+gj55PDIqdQ2O+49IeFSFleFXGpKfMkY8yNjzCvAz4F1wIvAJwAdC1rFMh01V2nyGvj7TTCMjFDX0uJ5lV97uTiwcydP3HwzQ/39YC1D/f08cfPNCmGRIqm6Yz39lvlmGhDnOrBtXpx//yJAOSxSepmO/aw4QebwWNw/c5XFVSGXlRLnAN3Ae621y6y1V1tr/9Zau9fa9I8epKpkOGqu0uQz8HcnGBrmjNEwZ2yyPbGnuxILzO3dupWxtNdpbHiYvVsr9A2SSMjK+ljPfOzePPXTtpFTTrsfv4HtzHlQPzO1rX7mZL2JMqUcFim9TMd+Vpwgc7h5kZO5yuKqlfWkhLX2y9babuBNY8yXjTH3GmP+yRiz0RjTOO0DSMXy/bQ/7Ti5KFWozyRTP/MZ+GeaeGhuayN262ZoaACc1zJ26+ayrycxdNT7l69fu4gUxu/4Tq/2SshiBg/n1g7OwNakDXvqZ8JH/xLaboPaGU5b82LncpnvYVYOi5RerMk7i9OP9ozaaSG+trc6X17yzWG/iYfl65XFVSynmhIJd+IUu/xr4P8A5wN3BdkpKS+V+mm/l1wG/q6JiQdjnNumTTw0t7Uxa8UKZl1yCec+tLvsJyQAZi30Plfbr11EClOOOVzQ5EjzotzawRnYzv89wCRumzTgXb4eFl0CSz4IG58r+0EwKIdFwlB2R3tmmnSYTr453HYbnjnsXq8srkr5TEqcZ639nLX2Z4mva4F3Bd0xKR/ZfNo/XXHI4RdeKItP7vId+De3tVE7eza1c+Z4TjwsuevOiqp4v6Kjg9q016m2sZEVHdF9gyRSzrJddVUpWZz3Mt+mBTBjjveAd8Mu56tCKIdFSq91WStdl3XRUONkcawpRtdlXSknbrinczz16lOs7V47pQjmC2+8UB6rKPLN4eXrnRye8TbviQdlcVXK+kjQJM8YY95vrf0FgDHmfTincEgVy3TUnF9xSPd+S+66szwGwTDxd+zbdCPE49S1tEw5dk9g6bp1gLOPbujoUWYtXMiKjo6JdhEJXqYchmhl8XTHl07LHcQ+cJ1TOK158eTyXwGUwyJhyXTsp9/pHO79tl+xvXQTEu7JGWOnnZMzcs1Q5XBWlMXZyWdS4n3ANcaYlxOXzwaeN8Y8i3Nq6PLAeicVIVNxyHJ8M9/c1sbAjnsAKmp1Q9CWrlunwBWJEHcyNVkYWTzd5EjWlq+Hp+9wvq+gT9WCpBwWiZZte7b5ns7hN5FRFH4nZ0DuExPK4Wkpi6eXz6TEFYH3Qipa1VWFFxGJooiclBSVyRERkVLzO4Wj5KdzZDo5QysdJAQ515Sw1h7K9FWMTspU0+0LjpJ8ikP6KZeq8SJSHcoqi7M8KSkbBWVxRCZHRKQyTFejIUrST+GYrj2Tgk7wyOfkDJEiyqfQpYTMb+lrVAfD5VQVPts3GJVWmFJEcqcszk82kyNZT/ZUWEE0EcmNX42GqE5MROZ0jmxOznBrThx62Kk50bvD+z7KYQlAPts3pIiy+eTJLQ6WzA4P07fpxolaB37CeCM9XXHIggueBSSwfc4iUvaUxcXL4gUbO+i/qTOl1lDy5IiyWESArFYB9L7WOzEh4RoeG6bzkU66X+rOeN/tV2wvqH/5cOtGdD7SSXw8TqwpRvvK9ol2d9VHfDzO2u61KdcFak2nU0MieQtH8skZQdWcEMmSJiXKkc/SV9/2CPArDlmKwWe2S4wLeYMBKnopUnWUxTlJzuL6lhbiBw6AtdDQQH1LCwM77mFgxz1lN9kjIuFJn5CYrj0KWpe1TkyYJE+MTHcyRyC2Jz1O82J4fR9goXaGc/npO5wv91SOZCOnnJM23MKWmWjlhORIkxIRk82Aat/qNYz29U1pr2tpifSAzKtvfgXP+jbdWPpPxMrwDYaIFEc1ZfGhq6/xnQgoVvHJuvnzqZs/3/tKZbGIkN1KhrXda+k/ObUeTawpFspKiGx59c1dPZHMXfVRlNUSTQucLy/pExLTtYsUSJMSZWi6pa9lJY/B5+jx486bgf5+6mKxaZcYZ/vmoFzfYIhIOKohizMVn1QWi0jY2le20/VoV8oxm6HUaAhAXqs+Tr7m1HsYPOzUg1jTmXl7RbYrGLZc6GzZSNe8WKsgpChU6LIMNbe1Ebt1MzQ0AM5ALXbr5rLca+tX8MzMnTulbbCnh6E9e4jv3+8MWK0NtLBcVIrAiUh5qJQsXnLXnTmdzKEsFpGoaF3WStdlXTTUODkca4rRdVlXcVYWFFmsyfskpOaG5iltu/bvovfYMzz11kHWzhljV9PMyboPfgUpc7Gm06kxkSy55oRIwDQpUaaa29qYtWIFsy65hHMf2l12g2DXgo0dUOexYOfkyZTB7cR+57GxKTd1lxgXqlLeYIhI6VRKFs/+0OVZtSuLRSRqWpe1svwdy1l15ioevOrBspyQAGfVR52ZOiYeGh1KOU1kovaEHeO8+AjdR/q5fOgUx2tqnLoPuzcX3pnl66HtNqfWBDgrJNpuU5FLKRpNSkioBnbc4xQ6S2NHRlIGt8e2bE1ZIp0uqPPtK+UNhohItg5dfQ0D/3Sv53Xp7cpiEZHi6H6pG8vUMfHI+Ajb9mybuLxtz7aJ7Sov1zuTGBY47H7IN3g4mA4tXw+LLoElH4SNz2lCQopKNSWy5FYN137WIvD4xA1SB7fTDXS9lhhHzZL/4r08WkSyoxwuoizr+1RCFrNwufZEixRAWVw8Y9Z7THz05FHP70/V1PDiDGdV2ecWLqD34CtObYmoW7g87B5IxGilhBTVoauvyXgkZ7Z7mTMNdP32Gg/29LBv9RqeP/8C9q1eE8heZxGRclNQDqe1K4tFRPKz4Scb2PCTDb7Xb79iu29diYVNCz2/B/hs7Ew+GzuThaNj/nUfenc4xSu75jp/BlF3QiRAmpSQ0C3Y2AHGpLSlD269Cp8B1M6d67nX2N33XIwibCIilSbb4pLKYhGR4mlf2Y4hdUycfppI+8p2GmtTc7hxfJz207XedR96dzgFMAdfAWywBTFFAqLtGxINNTUT2zjM3LnENn09ZXDrft+36UaIx6lraZk4fu7Q1dc4tSmSDO3dO2XZsR0epm/TjVNuW5Llh7074PCTzvnOWy6c/sgmH1oyKSLF0NzWxtCePQz88EdOQ20tzR//2JRJhmJlcTnlMADbE4X0tA1ERAJWY2omtnE0NzRzw/tuSCne6X7f+Ugn8fE4saYY7SvbnfbtrfD0HakP6OZespFT8MB1U29bikwLKouVwxVFkxISKs9K7j5F1Jrb2iYGsdMOYLPcH10S7gy1+wvBnaGGnEJ4sKdnYoC/b/WaiTcCIiKFGuzpYfC++ycbxsYYvO9+Zq1c6TkxUXZZHFAOTzxWUJMbIiIJ7qkayXUlTqdPJiS0Lmul+6VuwNn2kZHPY/i2F1NQWawcrjjGepx8UC5WrVpln3rqqaI/z2BPj+enQmGL3Kfm21OPYBo9/jov/9MbYC0L33+KhkWLqJt/Rspt9t/+ElhL/K1aABrmJILYGJZd+64pT3HoB33A9EUjh/b+Eusx6DUNDcxacVHWf6VAeM1Qg3PM0qJLsnqI0eOvEz94kKNPzQHgzJWDmJoaGs45Z8prWtYKmO2O0s+DMeZpa+2qsPtRCtWewxCt/3tAShZnk8OHftDH2IkTlZ3FAeQwACePwfHfgB2fbDM1MP/3oGlB4f2MgkI/dYzQp5fK4uBFNYvdWg3TvkEvlbQx8fHh43yhbgBrLTcNnGTRnEXMb5yfcptPjv4Wa+3EiRpnj4wCYIzhnrp3TnmKDeZV56nsmZn7ElT+BSGIvlRDDkNhGVqGOayaEtOY+CQ/MajSftjsuG+i3eM+bTxO/OBBRo+/nnpDv0mxAifLGhYtwtSk/vc2NTU0LAqhInEAM9Txw4ex4+PE36ol/lYtr+5pdi4fDujYJ5EIUw7nJ+schsrP4qA+KXzzUOpAGJzLbx7Kr18iZURZnJ/jw8c5OHgQ94Pg+Hicg4MHOT58POV2fh8UF/wB8tuXOG/ak5kap73Ugshi5XBFKvn2DWPMYuBOYCEwDtxurd1mjJkH3A2cAxwE1ltr3yx1/5IduvqajPthozAzHCmJ2TjndTsM8XkTV738kPPJXF3LbM59aHLWbuyuNU4BtITTAzWJ27V4z+49lKggvyHzJ5J1wMmeHo5t2cpofz91sRgLNnYwK4x/sy0XJooLpWlenPUM5qHzL0h5c3B6oMZ5TY3h/K3hz4KGTVtbclNJOTyw457orFCIikSu7Fu+wjOHaRjm/N7J3FiyAfatrvAsDiCHAadyvZexeCQ+kQqdllTnpJKyWL9z0yTyYMNPNtD72lHis94xcdVnY87KhljT23jwqsnc+F33WvpPTh677B71GWuKwVUe+eKe5JHN6pDeHbB7Mwwedo4MDetnM4gsVg5nVqY5HMZKiVHgK9ba84H3A39mjLkAuB7Yba09F9iduBy+qOyHLTc+r0/6GffZVnzPR3NbG+c+tJvzn/815z60O7xfmGs6nSOakvkd2eTD7xi+TMfzVQt9cpMX5XA1yOF1q/gsDiCHAWcwn0t7NfHbK64K/5koi6tAfNz79Tl68mjKZc9TNdJO3sjb8vWw8TnoGnD+DOtNahBZrBz2V8Y5XPJJCWttv7V2T+L7t4DngbOAKwG3BOwdwMdK3bd0uZzdLpMyvm5pb6Kb29qI3brZub0x1LW0eB4rl/zYZfep6PL1zhFNtTOcy82LvY9syqCYbxjK2aGrr6Fv043YtOKo7ic34q2Scrjs8qCEcvn9VfFZHEAOA8FNblSa7a1OJf+RU6ntI6ecT2fFUyVlsXjbfsV2Z6WDh4VNC1Muty5rpeuyLmJNMQyGWFOMrsu6Uk7eSH/syNTQyFYQWawc9pYphx+4Lpw+5SDU0zeMMecAFwOPA2daa/vBCWljTCQqlSzY2EH/TZ0pb3qi8mYwygPCXF635ra2yl/2t3z95LFLeSwty3QMX9XTJzcFUQ4XLqpZnOvrVvFZXGAOTzwGOAO8sdPOgLpMlsYWnd+e8EHVPsqGsrgwUX5z3r6yna5Huxgem3zd/FZAtC5r9Z2EqBiFZrFy2F+UTlrJUWiTEsaY2cA/AR3W2t8ZY7K937XAtQBnn3128TqYoDeD+dHrFryKf8OQhyV33TllL7xLn9xMTzlc2fS6Fcny9Rr8ptuwK8NecS2pno6yuLK5kwydj3QSH48Ta4rRvrK98icfikk5PFXGHF5c+v7kKJTTN4wx9Tjh+31r7b2J5leNMbHE9THgmNd9rbW3W2tXWWtXveMd7/C6SeCa29qYtWIFsy65JNzaBCVy6OprJo64K0S1vW4SDm1tyY9yOPqCyOJqfN0kJFpSnRdlcbRt+MmGieNGC9G6rJXl71jOqjNX8eBVD2pCQoqjjHO45JMSxpn+/TvgeWvt/0666sfAZxLffwZ4oNR9E5Hy4+6Fp8GpUj3dXnhRDotIEQRVt6OKKItFJFBlnMNhbN/4AHA18Kwx5peJtq8D3wR2GGM+B7wMfDKEvolIHg7s3MnerVsZOnqUWQsXsqKjg6Xr1pXs+bW1JWfKYZEKFHYWa0l1zpTFIhVGOZyfkk9KWGsfBvw2y60pZV9yEdVCZlGn163yHdi5kyduvpmxROGrof5+nrj5ZoDShrBkTTlcXfS6VQdlcflRFleXKBfjlGAoh/MXSk0Jia7Bnh6G9u5l6Mkn2bd6DYM9PWF3SSJu79atE+HrGhseZu/WrSH1SKT8KYslV8pikWDt2r+L3td6eerVp1jbvZZd+/M8tUeqhnI4f5qUkAmDPT3039Q5cZTiaF8f/Td1lvVgeLCnh32r1/D8+RdMGdhnuk6yN3T0aE7tIpJZpWXxtFnbu8OpGN411/mzd0c4HS1zymKR4Ozav4uuR7uIjzs53H+yn65Hu8p7YiJT1iqHA6Eczl9oR4JK9LhHPCWzw8Mc27K1LPfruwN79zxtd2Dv8ruuHP+uYZq1cCFD/f2e7SKSm0NXX8PQ3r0Vk8WZcri5rc0Z+PZ8CUZOJe7winMZynJPbJiUxSLBcY/vTDY8Nsy2PdvK8+SMTFkLyuGAKIfzp0mJKuEeK5dxH2DaINg16vHDVUzZHoE3/MILGa8fP3ECrE1ps8PD9H3tLxIXPK7bdCMDO+7J+rkb3/3urPoKwNE+58+Hrqmo/ZgrOjpS9s8B1DY2sqJDR3KKpCubLN6e/aD70A/6fK8bO3GChkZLfKQWgIY5Y8AIx7/xZQZu64TTb4GdBcxKud+S8evg6Tumf/Kjvc6fC5dn11n39ttbnTPdK4iyWCQ77hGfmWo8pE9IuI6eLPEn3tlmsZttfk6/BaSOexk5Bfdem7jgcd0DRcrh5PtUWBYrh/On7Rsyoa6lxbs9Fpv4/tDV12Q9aRC6tEmHlHa/63zeDIi/pevW8d5bbmFWLAbGMCsW47233KKCPiJ5WHLXndNmccXkMIAd975+7HRx+lPBlMUiwYk1xTzbFzZNfuK94ScbJiY4os8ni7H+1ymHc6Yczp9WSsiEBRs7UpbZApjGRhZsLO3sXlCrCPatXsNo39RP8NwBv9912Tx/Vp92pnNnuzdUzioJ19J16xS4IgGJRBbn8MnVkgxj8vQcPj3gfBZS19LCsn/d7exdHnxl6h2bF2fXh4lczbK/ud6+zCiLRYLRvrKdrke7GB6bzOHG2kbaV7aXtiNBZVWmrIXS5nC+9ykTyuH8aKWETGhuayN262ZoaACcQWPs1s1lt4fZtWBjB6axMaXNHdgv2NgBxnheJyISpkrK4kw5DED9TDBpQ5H6mbCmExGRsLQua6Xrsi4aapwcjjXF6LqsqzzrSYCTqfUzU9vcrF3TqRyW0GmlhKRobmubqKlQ7nUP3AH8sS1bGe3vpy4WY8HGjon217/9HeKHD8PIyJTrRETCVClZPF0O07TA+XPkFAwehuZFzkBYxdVEJGSty1rpfqkbyFx/oiy4mbp7s3fW/vxb8OYhGIsrhyUUmpSQsjbdlZgKpwAAIABJREFUNormtjbPiYbBnh5nQiIep66lRRMSIiJ5yjeH6d0Bh5909i03L4ZP3K5BsIhIvqbbErF8vXfG9u5ITEgkslgTEhICbd+oAoM9PQzt3cvQk096nxFfZdwj6tyilu4RddX+uohIcSmLk7jH07mF1Nwj6Hp3hNsvEalou/bvove1Xp569SnWdq9l1/7Kq2mQE2WxRIRWSlQ4vzfggO/KgKCXCudVFLKI+jbdOOWUDTs8zLEtW7VaQkSKIuwsjloOs3uzs2Uj2cgpp12f0IlIEezav4uuR7smjvvsP9lP16NdAL61IoLetpHNcaQl9cB1U0/ZUBZLCDQpUWHSj4kb2rvX8w1436YbJ/Yru6YbrLqf8hGPs2/1mvLd8uBz7Odof3+JOyIilcjruM5ss7hqcnjwcG7tIiI5Sj+us/e13okJCdfw2DCdj3RO1I5wTTdp4K64iI/HWdu9lvaV7eVZBNPv2E9lsZSYtm9ERNHOnfd5A+7b7iPfLQ9RXK7sHgk6pT3mfSa1iFSPKGdxJeUwzYtyaxeRqlGsHE6fkJiu3Y/fiovptoJEcuuIeyTolHZlsZSWVkpEQJCffKV/ypZ+RryrrqUlp2W8+Wx5yGe5ciks2NhB/02d2OHJs6d1HKiIBJXFXtkaRBZXUg6zptPZt5y8hUNH0IlUvSDHxOmrHdZ2r6X/5NRVsbGmWE7bKTof6fRccbFtzzbf1RL5bB0pCWWxRIRWSoSs2EUXpz0jPlt5bHk4tmVryht/mBxAh6m5rY3YrZuhwTl7uq6lhditm8tzCfQ0ivapr0iFKYssrqAcZvl6aLsNamc4l5sXO5crcQ/z9tbJqvgi4qvYOdy+sp3G2tQcbqxtpH1le06P47ey4ujJo7732bZnG8NjqVnsTmSESlksEaFJiRK676/2cN9f7UlpK/aAMag34PlsefAbKEehdkNzWxuzVqxg1iWXcO5DuytyQkJEpvrUdx/jU999bEp7OWRxpeUwy9fDoktgyQdh43OVOQgWEU9eH5oUO4dbl7XSdVkXDTVODseaYnRd1pXzSoVYk3fmLmxa6HsfvwmLTBMZJaMslgjQpETIMg0Yg/qUO4g34Pl8yuc3UI5K7YYld90ZnUr0IhIq3yzu6wtstVGhWVyJOcyGXc6XiFS9UoyJW5e1svwdy1l15ioevOrBvLZO5LPiwm/CItNERkkpiyVkmpTIgddKh0JFfsCYkM+nfJ4D6Pp6xoaGeP78C6YUXAviF462K4hUNr+VDoXyzdxE5kVBYDnc2MjsD13OvtVripbFWiYrUtmKMd4qlzFxPisuvCYy6mvqGRoZYvkdy6cUvtzwkw1TTg/JmXJYyogmJUIWWM2HEsj1U770AXTt3LlYa7EDA2Bt4HsFI0+z0CKR5ZfFDYuiVYG80Byua2mh+eMfY/C++53Cm9WWxcphkcgqpzFxrisu0icy5s5wxsSD8UEsNusTPCqGsljSaFIiINOtonjp8aMcPTBI374B7vj6I7z0uLOHrBhFF6O0WqC5rY3ze/dy/gvPY2bNgtHRlOvdvYL5HFmXfp/+W24pybF3kTxeT0SAzCsp7n/mCM+8PMDjB97gA998iPufOTJxnV8W182fn3dfopLFyTl87kO7OfEfP/fdt51rvnnevncHHH4SDj0MWy50LhdDqZ5HRHIyXfb55UwxxsSBrDgISOuyVp6++mme/cyzzKybyahNHRO7hS/zOjo0PQ93flk5LGVFR4KWwEuPH+Vn33+B8VELwIk3TvOz778AwLvet5DmtrYpgRvkkUhRkWnPdq5H1nlVaB744Y+mPmaGx8hHZI/X81CJ/4dE8nX/M0e44d5niY+NA3Bk4BQ33PssAB+7+CyAKVlcqT9DQWWxZx5u2gSXvEnz4tOJG73iHDcHwRZP693hPO5YkZ+nUO6Afey0M2Bf0xmt/omU2HTjqGoZE/sVuHRXTOR0dKhXHj71d5PXV3sOg7K4DBhrbdh9yNuqVavsU089VbLnc1dCfPwrKz2ve/3wW5yxaM6U644eGJyYkEhWU2dYuLR5Svvqdx2h/6bOlE+yTGNj1rPF7ux0MYo4FvKpn/sLJWvGUDN79sTFxne/O7/Hqq2lZtas7J83g/ETJ8DrZybR1+Q+FqqQfz/3l36+/4ekMMaYp621q8LuRymUOofdVRB3f/5Sz+vf0/VTAC6IvS2l/ZmXByYmJJI11NZw8dlzp7SPHT/Of99xU0E/Q8XK4kJXXwSVxcmPY2rHaZgzlriFpbbeIydr6qChKd9upzr9Ftip/56YGpgxhyX/xfukkrzku8TYHbCPnJpsq59ZucftRZCyuHgy5duhq69h+IUXPMdEvvnT0MCsFSumNM9d/8mCxjPuKontV2yf9ra5mnjs/mM533fva3uJj8f5szPfAcDZI6mrJl6ur0tpN8Ywu34yh7fbMydv7L7ZzoapDTaH8XofaWDGHFi4PJjngcK2eiiLQ5VtDmulRADcrRnjo5ajBwZpPmMmM+dMFkfzmpDI1J7pSKR83lAO9vTQt+lGiMepa2kJbYa5YdEi4gcOeL+p95LpdrkMqIPk16cITO4lv1Hx+qVvh4fp23QjAzvuAYozaSUSpvufOcKJ4VEsziTE4nkzOWO2c/a614RExvbDhwPNYajALPbNYYP3QDVAXhMSmdpLJbmonNcbhZFT8MB18PQdzmXtqZYKk7yqYWjvXhoWLUrdBueXGz7tQY+Jd+3fRecjncTH48SaYrSvbM/rBI5CLZqziIODB6dMRvjJ+CFythMSgfPrU/hjYmVx+dGkRIHSt2aMj1reemOYVR89h3e9zznm546vP8KJN6YGxux5MzxXXTz/vfzPlU9f4jb7Q5czeN/9BW03SJ4NL/SN7GBPD8e2bGW0v5+6WIwFGzucy319U25b19LCuQ/t9nycfavXeN4nl8fIh9/zBv08Bcvxl75IuXO3Z7hDofjYOH0Dw/z56nP52MVn8YFvPsSRgVNT7nfW3Jmeqy6eP/9zns+TTQ5D8FkcZA67/Ss0i5Pz0I7VcHrAKVNVN2uUZVccT71x82LY+FzB/Z6w5UJnqXC65sWw8ZngnqcQfm8UQnsDIVJc6VsziMcZ6evjjC9+YSLnMo2jvLLt+fMv8HyubLLYrc0QH4+ztnstly+6nAd+80BuWyPSJK+8KGT1xXzgF/t38bs92zh68igLmxbSvrKdbXu20X9y8u/24ozJEz66r3rQ+8H88jBdSXM4wOcplLK4LGhSIg/uNo7XD79F/NTYlOtH4+Psvut5fvWwE7r1DbXUNdQwGp/8BKeuoYZLr3yn5+PXxWLegT3NkUjT1VlwuZ+Yh/EJnddeQcBzaV6massLNnZMuU+6YlRs9nreqFSGTv5lnusvfZFy427j+HX/7wA4PTI+ZdXDqZExvtbdyw+feJnG+hpm1tdyamQys2fW1/LVj5zn+fj55jBUTxZ75mFDPQsuPpF6w/qZzv7dIK3p9F6OG/Tz5Cr507ZMA3Z9KicVwp0wHX7hBc8trumrNGsaGzGNjVnnTL5ZvGv/rim1Ge5+8e4ptxseG6bzkc5QVku0Lmv1fN6uR7sYHpt8fRprG2lf2e7/QF55mK6achiUxWVIkxJZmFgWf1H2bzyTt2bMnNPAqo+ew+67nmd81DJ73gwuvfKdEysp0mX7xjd9X3FO+4Q9blfMWhSZuAPjXJY1e91n9ocuZ+Cf7i3q0uh8+hqGKE+eiOTr0NXX8Kf9v+MfPvkXU66bbnvGGbNn8Oerz+Vr3b3Ex8Y5a+5MvvqR8yaKXKbLN4ehwCx2l5yGMFDKNd98b78ksSx27LQz6CtGQTH38Yr9PIWI8oBdJE8Za9pksdWrbv58zvjiF7LOmWyzOP2EDXeFRDbSb1fMOhTTcScpctpi4pWH566FZ/5ROQzK4jKhSYk8uFsu7vurPb5FLL22ZrgrJ7y2bCTL+41vDkvz61oCLAIWgOa2tpxrHXjdJ3bzzcXp4DTPGzXlMnkiki93y4W7YuLwm6ey2p7xwydeTrm/n4J+hqooi31vX4pB6fL10d4PXC4DdpECuD/3h66+xndC1muVZrY5k28WZzshAc7WiChpXdZK90vdQA4TI155uO5/F6F3WTxv1CiLy4ImJfKUXNwyXaatGdnKZmCY3u5bZ8GYlNnrqH5iHtU3+FFzYOdO9m7dytDRo8xauJAVHR0sXbduyu3KYfJEpBD3P3Nk4lSNuTPrqa81jIwlrVLLsD0jG/nkMCiLq0VWWVwOA3aRAiXX0ClGzmWTxelv3td2r02pzeBn2q0RIQljlUY5ynZMrCyOvpqwOxB1btAOPfkkp/buZfT48SnFLZPNnjeDP/iTd3tuzfj4V1ZOu0qiEAs2dmAaG1PaTGMjcz/9KWhwCuXUtbR4HqOU/Pfct3oNgz09BfVlsKeH55ev4Pl3nx/I44njwM6dPHHzzQz194O1DPX388TNN3Ng586wuyZSVG5Gjf3uLU7t3cvuv76TG+59dmJ7xsCpEbBQV2MAZ4XE//jEe6Zsz7j785dOu0qiUHlnsXuO+qGHnT2wvTsK7ouyuDiUxVKNkseKQ4kx8ZTilkkTEn5jTgiuaK+f9pXtNNam5nBjbSOfOu9TNNRMFo/suqwrZWuEWxzzqVefYm33WnbtL/wN7K79u/j9u36f99zxnsAeU5TDlUYrJTJIDtp4/WwGGluwx2v41+2/wjnyLFVNneEz3/hA6TuakGmJW6ZtDV5F2XI9oaOYjyeT9m7dylhacc+x4WH2bt3qPTMsUgFSMqURbDzOvO/8LQ1/+HVONc6ZuN3IuKWhtob3LX170SceMskri91z1N1q4IOvOJch7yWmyuLiURZLtfE6WSO+fz99X/2a9x0aGkI9lSxTbYYb33+j5328imPmejpHKR5THMrhyqJJiQzcs5HHTQ0146M0nXyVE3MW+d7ea+VEKSQXqMx2yX5yoSKvPYDp1ZJz6Uumx9NAuDBDR4/m1C5SCZLPqW85+ToAjWMjLBx6g8GkSQmAkbHxUCYk0gsFZ5XFuZyjfrQPFr4np/4oi4tHWSzVJjmHszIyUrzOZJBcpDLb2gzufbyKY7qnc7iPk2tfMj2mJiUKoxyuLJqUSOMe9wkwtGA9LLBYagDD0Mx3ZLxvTd3U1ROFyGZZW/I+vn2r1+S3b8+vKFsOxdqK+ngyYdbChc4yNY92L9obLuXILV7pGlp2JZ8/8QAtJ19n3Bj2N7fwtQ9+kYbx0Sn3ra8Ndldivjmc85v+YpyjriwumpyyWPuXpUwlf4DlWScnk/r6QPuSTRa7Wy/i43HWdq/NuVaEX3HMXIpmluIxxZHrmFhZHG2qKZGBSez9HaudwVjtDMbrZmS4MTQ21fPt637G33zhIe74+iO89HhxZ+r8luaOHj8+7X1T9vIl/p7p6lpamLv+kznVmlhy152+1eSjVmW+HK3o6KA2ba96bWMjKzqiVyxPJCimoYHvLr+SvqYzOPi2ySrp8ZrUefUaA4vnzeT+Z47wrk3/wjnX7+ID33yI+585UrS++eVwVrUbNuya+Dr0Hy0c2j1/6m2aFzN4xhcYenU8p5o/yuLiUhZL1fEZK3oyhtq5zSWtZ+O3TeL48PRj4u1XbGf7Fdsnak2kizXF2H7Fdq5611U51ZvYfsV235M9onbiRzlSDlcWrZRIk1yIcrDnCP03dTJQv5Dxmnqeudj7P/nseTM458L5vPCLoxNbOE68cZqfff8FAM+il4XKtDQ3fuAAdfM9Brc+GhYtYqSvb8oZ0LM/dHle+5GzPVNacufukcuq0rBImUrffjHY8zoHrt/EjLERfvu2FpYN9rH5se+x7aKreHHeEsApbOmetJFcAPPIwCluuPdZgClFL4Pg1o1I5m6RmLViRfYP9PZz4Pi+1Lb6mQw2XJl3XQhlcfEoi6UaJK9OcCdgp9vCUdfSwuwPXc7gffeXrJ5Npm0SBwcPMr8xuzHxWXPO4uDgQSyT27Hd0znyrQ3RvrKdrke7GB6bfN2ieuJHuVEOVxZNSmTgBufpzm8w1JAaaHUNNcyZ18jMOQ18/CsruePrjzAaH0+5zWh8nMce+G1RJiUA/yW41ua0bL9u/nzO+OIXphRl89o/aIeHObZla8ZfKvmeKS3ZWbpunQJXqkpzWxvHD77JGd/+FuCskNh20VU8vuy9vHNuI2fMnjExkfGBbz7EqZGxlPufGhnjWz99sSiTEpm2SOS0faopsT2w9ncp56gf6/huXjkMyuJiUxZLNZnIk+tvgLHUjDWNjdS3tFA3fz5L7rqTfavX5J1b+fLbDmGxWR+v6U5eHHnryJTimGu716ZMLIAz6bFtz7aMkxKZCm5K4ZTDlUOTEtNobmtj7o57aHgrTk2dYXzUMnveDC698p386uHJ/XUn3vDe9+vXXig39L32+GW7NDd5H3T8yBEaFi2a+IUC0Pe1v/C836jH/q10zW1tJRn4qmaCSHVY8+fX8PRDuxh67SQvzlvCiWXv5n985Dx++MTLKbfrGzjleX+/9kLVtbQUlMOQmsVDZ9fTsGgpdRsfB2C0v8vzPtnkMJQmi5XDItXBLeA7evw48cOHUyY7kwuj++VTtrmVq+1XbGdt91r6T059/Gy3SSTXo2ioaWBp81J+/LEfT1x/9KT3lmy/9mSty1pLMwmhmglSxlRTIkuz5jSwcGkzLefO5TPf+MCU1Q+z53nXm/BrD8KCjR2YtL1U2S7N9doHHT9wIKUeRV3MO8j92kVEiumM2TOY01jH+5bO45HrV3uufGiZO9Pzvn7thSokh2FqFtt4nPjBgxP7r5XDIhI1dfPnM2vFCmZdcgnnPrR7ysRnGLnVvrKdxtrULM52m0T61oz4eJyDgwdTakYsbPJe9ezXLiK50aREFlKKQib5+FdWTtSguPTKd1LXkPpy1jXUcOmV7yxav5rb2ojdutn5RM4Y6lpaiN26edpPxQ5dfQ19m26cui/QWsaT2godbIuIBGnJXXfyD59MXcF19+cvTalB8dWPnMfM+tqU28ysr52oNxG0fHMYgO2t1O/+QkoWv/zQGRz6t3kc27IVUA6LSLT4jYmT28PIrdZlrXRd1kWsKYbBEGuK0XVZ17QrFDb8ZAOdj3RO2ZphsWzbs23iciGTHiIyvZJv3zDG/D2wDjhmrb0w0TYPuBs4BzgIrLfWvlnqvhXCXTnx2AO/5cQbpye2eBStnkRC3ktzffZBJy+t035kkcpVqVnsrp741k9fpG/gFC2JAphFqSeRUMgWCTtNFiuHRSpXpeawm0/HtmxltL+fulisJLmV7zYJv3oUyVszVBtCpLjCqCnxD8D/AZKnWa8Hdltrv2mMuT5x2bugQYS9630Liz4JEZiGBs+JifSlde7+QdC+YZEK8w9UaBZ/7OKzijoJESTjdyRzUhYrh0Uq1j9QoTlcqtpiQWioafCcmEjfmtG6rJXul7oBsi6eKSLZKfn2DWvtz4E30pqvBO5IfH8H8LGSdqqKJBdUw5iU67QkWKR6KItD1rsDDj9JQ+MQJu3jAWWxSHVQDocrubhlOm3NECmtqJy+caa1th/AWttvjFngd0NjzLXAtQBnn312ibpXGdILqmEnz2HWkuDcRPnTykNXXwNEu48SWVllsXK4QL07oOdLMHaaukaIrXqTvifmwrhRFucqqtXmtyeWdEe1fxJlGhOXQHpxy2TampGjKOecsrhsRGVSImvW2tuB2wFWrVplp7m5JDm2ZevU4pYADQ2c+9Du0ndIRMqScrhAuzfDyOQRpc3nnGLgt7OgbgZLlMUikiVlcf627dk2pbglOFs5HrzqwRB6JFLdonL6xqvGmBhA4s9jIfenIvmeD+1TaE1Eqo6yuBQGD3u3j54ubT9EJIqUwyWQXMQymV/RSxEprqhMSvwY+Ezi+88AD4TYl4rle250S0uJeyLF4tYMGXrySfatXsNgT0/gz3Fg507u//CH+cGFF3L/hz/MgZ07A38OCY2yuBSaF01pWrLmOEs+MSuEzkjgEvVCOPQwbLnQuRww5XBFUw6XQHoRS1esyXusLGVIWVxWSj4pYYz5IfAYcJ4x5rAx5nPAN4E/NMbsA/4wcVkCpvPuK1t6zZDRvj76b+oMdGLiwM6dPHHzzQz194O1DPX388TNNyuEy5CyOERrOqF+Zmpb/UynXcpbUr0QAAZfcS4HOBhWDlcO5XB42le201ibOiZWccsKoiwuOyWvKWGt/WOfq9aUtCNVKN/z7lUwsTx41Qyxw8Mc27I1sKJ5e7duZSztOcaGh9m7dStL160L5DmipJKLhiqLQ7R8vfPnA9c5A6bmxc6EhNvuoRL/D1aktHohgHN59+aM/765qLYcruRCdcrh8LhFLDsf6SQ+Hs+6uKWOAi0TyuLgFTmLy67QZZg+/pWVYXehYDrvvnL51QzxrSWSh6Gj3nsw/dpFgnb35y8NuwvBWL4enk6c+leBb7aqll+9EL/2PCiHJQoqYQzZuqyV7pe6AU02VBxlcdmJSk0JESmQb80Qn/Z8zFrovQfTr72claI+h4hUGI96IRnb81BNOVyKPeEiUoGUxcEqQRZrUqIKLbnrzoqY4ZZUpagZsqKjg9q056htbGRFR2XVJSlFfQ4RNuzSKolKU4J6IdWSw6XYEy6y/YrtWiVRiZTFwSlRFmtSQqRCNLe1Ebt1MzQ0AM6pKrFbNwdWTwJg6bp1vPeWW5gVi4ExzIrFeO8tt1Tc3rlM9TlERHwtXw9tt0HtDOdy82LnckB7mKF6cjjjnnARkUyUxcEpURarpoRIBSlFzZCl69ZVXuCmKUV9DhGpUCWoF1INOVyKPeEiUsGUxcEoURZrUkKkwmhrTuHqYjFG+/o820VEpqVtOYVrXuQsE/ZqFxHJhrK4cCXKYm3fEBFJU4r6HCIikkEJ9oSLiMg0SpTFWikhIpLGrcPRt+lGiMepa2lhwcaOQOtziIhIBu7e7weucwqsNS92BsEB7gkXEZFplCiLNSkhIuKhFPU5REQkgxLsCRcRkWmUIIu1fUNEREREREREQqGVEiIiPrRCQkQkZFohISISviJnsVZKiIiIiIiIiEgoNCkhIiIiIiIiIqHQpISIiIiIiIiIhEKTEhFy6OprOHT1NWF3Q0Skqn3qu4/xqe8+FnY3RESqlsbEItVFkxIiIiIiIiIiEgpNSoiIiIiIiIhIKDQpERGDPT0M7d3L0JNPsm/1GgZ7esLukohI1bn/mSM88/IAjx94gw988yHuf+ZI2F0SEakqGhOLVB9NSkTAYE8P/Td1QjwOwGhfH/03dSqERURK6P5njnDDvc8SHxsH4MjAKW6491lNTIiIlIjGxCLVSZMSEdC36Ubs8HBKmx0epm/TjSH1SESkunzqu4/xte5eTo2MpbSfGhnjWz99MaReiYhUF42JRaqTJiWiIDEbnHW7iIgEzl0hka5v4FSJeyIiUqU0JhapSpqUiIC6lpac2kVEJFh3f/5Szpo70/O6Fp92EREJlsbEItVJkxIRsGBjB6axMaXNNDayYGNHSD0SEak+X/3Iecysr01pm1lfy1c/cl5IPRIRqS4aE4tUp7qwOyDQ3NYGOPvoiMepa2lhwcaOiXYRESm+j118FgBf6+4lPjbOWXNn8tWPnDfRLiIixaUxsUh10qRERDS3tTGw4x4Altx1Z8i9ERGpTh+7+Cx++MTLgLOlQ0RESktjYpHqo0mJCFHwioiET5MRIiLh0phYpLqopoSIiIiIiIiIhEKTEiIiIiIiIiISCk1KiIiIiIiIiEgoNCkhIiIiIiIiIqHQpISIiIiIiIiIhEKTEiIiIiIiIiISCk1KiIiIiIiIiEgoNCkhIiIiIiIiIqHQpISIiIiIiIiIhEKTEiIiIiIiIiISCk1KiIiIiIiIiEgoIjUpYYy5whjzojHmN8aY68Puj4hINVIWi4iESzksItUkMpMSxpha4G+AjwIXAH9sjLkg3F6JiFQXZbGISLiUwyJSbSIzKQG8F/iNtXa/tTYO/Ai4MuQ+iYhUG2WxiEi4lMMiUlWiNClxFvBK0uXDiTYRESkdZbGISLiUwyJSVerC7kAS49Fmp9zImGuBaxMXTxtjnitqr/J3BvB62J3IIMr9U9/yF+X+RblvEGz/lgT0OGGYNovLKIch2v/v1Lf8Rbl/Ue4bRLt/QfetXLNYY+LSinL/1Lf8Rbl/Ue4bhDAmjtKkxGFgcdLlRUBf+o2stbcDtwMYY56y1q4qTfdyE+W+QbT7p77lL8r9i3LfIPr9K6Fps7hcchii3T/1LX9R7l+U+wbR7l+U+1ZiGhOXUJT7p77lL8r9i3LfIJz+RWn7xpPAucaYpcaYBuDTwI9D7pOISLVRFouIhEs5LCJVJTIrJay1o8aY64CfArXA31trfxVyt0REqoqyWEQkXMphEak2kZmUALDW/jPwzznc5fZi9SUAUe4bRLt/6lv+oty/KPcNot+/kskxi6P+ukW5f+pb/qLcvyj3DaLdvyj3raQ0Ji6pKPdPfctflPsX5b5BCP0z1k6pmyMiIiIiIiIiUnRRqikhIiIiIiIiIlWkrCYljDG1xphnjDE7E5eXGmMeN8bsM8bcnSgGFFbf5hpjuo0xLxhjnjfGXGqMmWeM+ddE//7VGPP2kPq20RjzK2PMc8aYHxpjGsN87Ywxf2+MOZZ8dJXfa2UctxljfmOM6TXGrAyhb99K/Lv2GmPuM8bMTbruhkTfXjTGfKTUfUu67r8ZY6wx5ozE5ZK+bpn6Z4z588Tr8ytjzP9Mag/1tTPGXGSM+YUx5pfGmKeMMe9NtJf8tSs3Uc1i5XBO/VEOB9i/pOtCzeIo57Bf/5TF+YlqDif6oizOvj/K4oD6lnSdxsQ59i0SOWytLZsv4MvAD4Cdics7gE8nvv8O8MUQ+3YH8F8T3zcAc4H/CVyfaLse+MsQ+nUWcACYmfTCxPJIAAAK+UlEQVSa/WmYrx1wObASeC6pzfO1Av5v4F9wzux+P/B4CH1bC9Qlvv/LpL5dAOwFZgBLgd8CtaXsW6J9MU4xrEPAGWG8bhleuz8A/g2Ykbi8ICqvHfAg8NGk1+vfw3rtyu0rqlmsHM6pT8rhAPuXaA89i6Ocwxn6pyzO77WMZA4nnl9ZnH2flMUB9S3RHnoOZ3jtIpHFUc3hslkpYYxZBLQC30tcNsBqoDtxkzuAj4XUt7fh/AP/HYC1Nm6tHQCuTPQr1P7hFDSdaYypA2YB/YT42llrfw68kdbs91pdCdxpHb8A5hpjYqXsm7X2QWvtaOLiL3DOC3f79iNr7Wlr7QHgN8B7S9m3hC3A14DkAjElfd0y9O+LwDettacTtzmW1L+wXzsLvC3xfTOTZ8CX/LUrJ1HNYuVwbpTDwfYvIfQsjnIOZ+ifsjhHUc3hRF+UxTlQFgfXt4TQczhD/yKRxVHN4bKZlAC24vwnG09cng8MJP1gHMaZAQ3DMuA1YHtiKd33jDFNwJnW2n6AxJ8LSt0xa+0R4H8BL+ME7yDwNNF57Vx+r9VZwCtJtwu7r5/FmTGECPTNGPNHwBFr7d60q0LvW8K7gP8rsSzyP4wxlyTao9C/DuBbxphXcH5GbohQ36IsqlmsHC6ccjhPEc/iKOcwKIvzEdUcBmVxEJTFeYh4DkO0szj0HC6LSQljzDrgmLX26eRmj5uGdZRIHc4ymG9bay8GTuIstwpdYh/alTjLgVqAJuCjHjeN6jEskfl3NsZsAkaB77tNHjcrWd+MMbOATUCn19UebWG8bnXA23GWfH0V2JH4RCcK/fsisNFauxjYSOJTHaLRt0iKeBYrh4snKv/GkcthKIssjnIOg7I4JxHPYVAWF1Nk/p2jlsVlkMMQ7SwOPYfLYlIC+ADwR8aYg8CPcJZZbcVZQlKXuM0iJpealNph4LC19vHE5W6cQH7VXeKS+POYz/2L6cPAAWvta9baEeBe4DKi89q5/F6rwzj7w1yh9NUY8xlgHfAn1lr3hzHsvr0T5xfr3sTPxiJgjzFmYQT65joM3JtY9vUEzqc6Z0Skf5/B+XkAuIfJpXJR6FtURTmLlcOFUw7nJ+pZHOUcBmVxrqKcw6AsDoKyOHdRz2GIdhaHnsNlMSlhrb3BWrvIWnsO8GngIWvtnwA/A65K3OwzwAMh9e8o8Iox5rxE0xrg18CPE/2C8Pr3MvB+Y8ysxGyc27dIvHZJ/F6rHwPXJKq/vh8YdJe0lYox5grgL4A/stYOJV31Y+DTxpgZxpilwLnAE6Xql7X2WWvtAmvtOYmfjcPAysT/x9Bft4T7cQZMGGPehVPw6nVCfu0S+oAPJb5fDexLfB+V1y5yopzFyuFAKIfzUAZZHOUcBmVxTqKcw4n+KYsLpyzOURnkMEQ7i8PPYVvk6qNBfwH/mclKw8tw/tF+gzOrMyPEfl0EPAX04vynezvOHr/diX/Y3cC8kPp2C/AC8BxwF05119BeO+CHOHv5RnBC43N+rxXOsqG/walE+yywKoS+/QZnP9UvE1/fSbr9pkTfXiRRtbaUfUu7/iCTlYZL+rpleO0agH9M/N/bA6yOymsHfBBnL+le4HHg98N67crxK4pZrBzOqT/K4QD7l3Z9aFkc5RzO0D9lcf6vZ+RyONEXZXH2/VEWB9S3tOtDy+EMr10ksjiqOWwSTygiIiIiIiIiUlJlsX1DRERERERERCqPJiVEREREREREJBSalBARERERERGRUGhSQkRERERERERCoUkJEREREREREQmFJiVEREREREREJBSalBARERERERGRUGhSQqSMGGP+2hizxxhzSdh9ERGpRsphEZHwKYsriyYlRMqEMaYJWAB8HlgXcndERKqOclhEJHzK4sqjSQkpKmPMFmNMR9Llnxpjvpd0+a+MMV82xpxnjPll0tfvku+XdPuZxpj/MMbUJi5vMsb8yhjTm7jf+5Ju+6jH/U94fe9xu7nGmG5jzAvGmOeNMZcaYxqNMU8YY/YmnvOWpNsfNMY8m+jDU2mPNeU6Y0yDMebnxpg6j+c+xxhzyhjzy+R2a+1JIAb8O3Bb0uvxS2NM3Bhzht/fR0Sql3JYOSwi4VMWK4vF35R/fJGAPQp8EthqjKkBzgDelnT9ZUCHtfZF4CKARLgeAe7zeLzPAvdaa8eMMZfizI6utNaeTgRQg3tDa+1lBfR7G/ATa+1VxpgGYBZwGlhtrT1hjKkHHjbG/Iu19heJ+/yBtfZ1n8dLuc5aGzfG7AY+BXzf4/a/tdZelNxgjJmf6MdbwFjicU4BFxljDub9NxWRSqcc9rhOOSwiJaYs9rhOWSyglRISkMRM5gvGmDsSM7TdxphZwCM4IQvwn4DngLeMMW83xswAzgeeSXu4NTgBdMjjqf4EeCDxfQx43Vp7GsBa+7q1ti+pT76zvtP8Xd4GXA78XeJx49baAetwH7M+8WXzeY6E+3H+Ptm6EfhfwK+ACwp4XhGpQMrhvCiHRSRQyuK8KIurnCYlJEjnAbdba5cDvwP+30QgjhpjzsYJ4seAx4FLgVVAr7U2nvY4nwZ+mP7gidnZZdbag4mmB4HFxpiXjDF/a4z5UEB/j2XAa8B2Y8wzxpjvGWfvGsaY2sQSsmPAv1prH0/cxwIPGmOeNsZcm/Z4ftc9B2RVnMcYcw7O63c38DzOLzMRkXTKYeWwiIRPWawslhxoUkKC9Iq19pHE9/8IfDDxvTsz7AbwY0mXU/a4JUL2j4B7PB7/DGDAvZCYof194FqcwLzbGPOnAfw96oCVwLettRcDJ4HrE885llhCtgh4rzHmwsR9PmCtXQl8FPgzY8zlSY/neZ21dgyIG2PmZNGn/w5sttZaFMAi4k85rBwWkfApi5XFkgNNSkiQ0pdtuZcfxQnb9+DMhP4CZ1b4MpxwTvZRYI+19lWPxz8FNKY8gROI/26tvRm4Dvh/CvobOA4Dh5NmfLtxAjn5eQdwiutckbjcl/jzGM6+v/cm3db3OmAGMJypM8aYi4BPAH+T2Cd3EwpgEfGmHFYOi0j4lMXKYsmBJiUkSGcnCu0A/DHwcOL7R3CK77yRCMw3gLk4IfxY2mP8MR7L1ACstW8CtcaYRgDjVCc+N+kmFwFee+5yYq09CrxijDkv0bQG+LUx5h3GmLmJ554JfBh4wRjT5M7sJpa0rcX5RcM0180HXrPWjkzTpb8E2qy151hrzwFWoAAWEW/KYeWwiIRPWawslhzo9A0J0vPAZ4wx3wX2Ad9OtD+Ls8zsB0m3fRaYbZOq7yaKAP0hzpnDfh7EWQL3b8Bs4K8ToTgK/AZn2VoQ/hz4fmLp3H5gA7AYuMM4lZBrgB3W2p3GmGXAfcYYcH6mfmCt/Unicc7McN0fAP+cqRPGmNVAk7V2t9tmrX01EezzEr/MRERcymHlsIiET1msLJYcGGc7jkhhEkVndlprL5zmpoU+z8XAl621VxfzeUrBGHMvcIN1jn5Kbj+HHF/LxBK2Vdb/+CURqXDK4dwph0UkaMri3CmLRds3pKxYa58BfpaYmS1bidnm+9PDN2EMaE5UNJ7ucWYmblcPjAfcTRGRKZTDUx5HOSwiJacsnvI4yuIyppUSIiIiIiIiIhIKrZQQERERERERkVBoUkJEREREREREQqFJCREREREREREJhSYlRERERERERCQUmpQQERERERERkVBoUkJEREREREREQqFJCREREREREREJhSYlRERERERERCQU/z/AM4skq5KW4QAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "peak_pew_table = ew.tabulate_peak_model_pew(model_list)\n", + "plot_peak_pew(peak_pew_table)\n" + ] + }, { "cell_type": "code", "execution_count": null,
feature_idlower_bluelower_redupper_blueupper_redfeature_name
str17int64int64int64int64str3
Ca ii H&K3500390038004100pW1