{ "cells": [ { "cell_type": "markdown", "id": "vital-string", "metadata": {}, "source": [ "## RS Oph Gamma-ray Lightcurve - Figure 2 \n", "This notebook remakes figure 2 from the paper: " ] }, { "cell_type": "code", "execution_count": 1, "id": "first-update", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from matplotlib import rc\n", "from scipy.optimize import curve_fit\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "id": "environmental-shipping", "metadata": {}, "outputs": [], "source": [ "SMALL_SIZE = 16\n", "MEDIUM_SIZE = 16\n", "\n", "plt.rc('font', size=SMALL_SIZE) # controls default text sizes\n", "plt.rc('axes', titlesize=SMALL_SIZE) # fontsize of the axes title\n", "plt.rc('axes', labelsize=MEDIUM_SIZE) # fontsize of the x and y labels\n", "plt.rc('xtick', labelsize=SMALL_SIZE) # fontsize of the tick labels\n", "plt.rc('ytick', labelsize=SMALL_SIZE) # fontsize of the tick labels\n", "plt.rc('legend', fontsize=SMALL_SIZE) # legend fontsize\n", "\n", "plt.rcParams[\"font.family\"] = \"Helvetica\"" ] }, { "cell_type": "code", "execution_count": 3, "id": "essential-photography", "metadata": {}, "outputs": [], "source": [ "def format_func(value, tick_number):\n", " N = int(np.round(np.log10(value)))\n", " if N == 0:\n", " return r\"$1$\"\n", " elif N == 1:\n", " return \"$10$\"\n", " elif N == -1:\n", " return r\"$0.1$\"\n", " else:\n", " return f\"$10^{{{N:d}}}$\"\n" ] }, { "cell_type": "code", "execution_count": 4, "id": "vanilla-elite", "metadata": {}, "outputs": [], "source": [ "lat_data = np.genfromtxt(\"FermiLAT_LightcurvePoints.txt\",skip_header=1)\n", "hess_data = np.genfromtxt(\"HESS_LightCurvePoints.txt\",skip_header=1)" ] }, { "cell_type": "markdown", "id": "detailed-approach", "metadata": {}, "source": [ "### Define reference time\n", "\n", "Choose a reference time of t=0 to be one day before the optical peak (T_0)" ] }, { "cell_type": "code", "execution_count": 5, "id": "mathematical-ultimate", "metadata": {}, "outputs": [], "source": [ "time_vmax = 59435.25\n", "start_mjd = 59434.25" ] }, { "cell_type": "markdown", "id": "alike-contamination", "metadata": {}, "source": [ "### Power-law fit to decay beyond peak flux" ] }, { "cell_type": "code", "execution_count": 150, "id": "reflected-corrections", "metadata": {}, "outputs": [], "source": [ "def fit_func(x,a,b):\n", " y = b*(x/5.)**a \n", "# y = b*(x/1.)**a \n", " return y" ] }, { "cell_type": "code", "execution_count": 151, "id": "equivalent-conversion", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[-1.31150418e+00 1.47843782e-12]\n", "[7.20752666e-02 8.90798098e-14]\n" ] } ], "source": [ "popt, pcov = curve_fit(fit_func,lat_data.T[0][2:]-start_mjd,lat_data.T[2][2:]*1e-3,\n", " sigma=lat_data.T[3][2:]*1e-3)\n", "\n", "print(popt)\n", "print(np.sqrt(np.diag(pcov)))\n", "popt_err = np.sqrt(np.diag(pcov))" ] }, { "cell_type": "code", "execution_count": 152, "id": "changed-information", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[-1.42502230e+00 5.59694817e-12]\n", "[1.81660147e-01 5.14510390e-13]\n" ] } ], "source": [ "popth, pcovh = curve_fit(fit_func,hess_data.T[0][2:]-start_mjd,hess_data.T[2][2:],sigma=hess_data.T[3][2:])\n", "\n", "print(popth)\n", "print(np.sqrt(np.diag(pcovh)))\n", "popth_err = np.sqrt(np.diag(pcovh))" ] }, { "cell_type": "markdown", "id": "adolescent-nancy", "metadata": {}, "source": [ "Identify moon break in H.E.S.S. data taking and which Fermi-LAT data points are coincident with H.E.S.S. observations" ] }, { "cell_type": "code", "execution_count": 153, "id": "administrative-publication", "metadata": {}, "outputs": [], "source": [ "#Times coincident with HESS are short 6 hour bins\n", "mask_hess = lat_data.T[1]<0.2" ] }, { "cell_type": "code", "execution_count": 154, "id": "local-hamilton", "metadata": {}, "outputs": [], "source": [ "moon_start = hess_data.T[0][4]+hess_data.T[1][4]\n", "moon_end = hess_data.T[0][5]-hess_data.T[1][5]\n", "xmoon = np.linspace(moon_start,moon_end,10)" ] }, { "cell_type": "markdown", "id": "bound-departure", "metadata": {}, "source": [ "Define ranges for butterfly plots" ] }, { "cell_type": "code", "execution_count": 155, "id": "medical-happiness", "metadata": {}, "outputs": [], "source": [ "hx1 = np.linspace(hess_data.T[0][2]-start_mjd,5.,10)\n", "hx2 = np.linspace(5.,hess_data.T[0][-1]-start_mjd,10)\n", "\n", "hf0 = fit_func(hx1,popth[0]-popth_err[0],popth[1]+popth_err[1])\n", "hf1 = fit_func(hx1,popth[0]+popth_err[0],popth[1]-popth_err[1])\n", "hf2 = fit_func(hx2,popth[0]-popth_err[0],popth[1]-popth_err[1])\n", "hf3 = fit_func(hx2,popth[0]+popth_err[0],popth[1]+popth_err[1])\n", "\n", "fx1 = np.linspace(lat_data.T[0][2]-start_mjd,5.,10)\n", "fx2 = np.linspace(5.,lat_data.T[0][-1]-start_mjd,10)\n", "f0 = fit_func(fx1,popt[0]-popt_err[0],popt[1]+popt_err[1])\n", "f1 = fit_func(fx1,popt[0]+popt_err[0],popt[1]-popt_err[1])\n", "f2 = fit_func(fx2,popt[0]-popt_err[0],popt[1]-popt_err[1])\n", "f3 = fit_func(fx2,popt[0]+popt_err[0],popt[1]+popt_err[1])" ] }, { "cell_type": "markdown", "id": "unexpected-arlington", "metadata": {}, "source": [ "### Plot Figure 2" ] }, { "cell_type": "code", "execution_count": 158, "id": "balanced-swimming", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1IAAAGhCAYAAABrv4WZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdd3iT5frA8e+bpmnS3TJaQOYBVETlCC5Q2YocQUUEVDyKOI44caGiRxQVnDgRBT3gRFRUHCCgAo6fCh5xoR6RXWZLoS1dafP8/rgbutI2SVPScX+uK1fpm7zJEwjw3r3HYxljUEoppZRSSinlP1u4F6CUUkoppZRSDY0GUkoppZRSSikVIA2klFJKKaWUUipAGkgppZRSSimlVIA0kFJKKaWUUkqpAGkgpZRSSimllFIBsod7AQ2ZZVlTgHsqHne5XHTr1u3QL0ipIPzxxx8AHH744WFeiVKqqSksLES3YVGNlWVZOByOcC9D1dL333+fboxp4es+S/8BC71evXqZNWvWhHsZSvnll19+AaB79+5hXolSqqlZv349kZGR4V6GUnXC7XbTuXPncC9D1ZJlWd8bY3r5uk8zUko1cRpAKaWUUkoFTnuklGrivv76a77++utwL0MppZRSqkHRjJRSTdydd94JwIoVK8K7EKWUUkqpBkQDKaWauOeffz7cS1BKKaWUanA0kFKqidNpfUoppZRSgdMeKaWauJUrV7Jy5cpwL0MppZRSqkHRjJRSTdw998hWaNojpZRSCmR/r8LCQjweT7iX0qAZY/jtt9/CvQxVhcjISFq2bEl8fHzQz6GBlFJN3EsvvRTuJSillKonCgsLcbvdHHbYYURFRWFZVriX1GAZY3A6neFehvLBGENeXh5paWkAQQdTGkgp1cR16tQp3EtQSilVTxQWFnLYYYdpAKAaNcuyiI6Opk2bNmzfvj3oQEp7pJRq4pYvX87y5cvDvQyllFL1gMfjISoqKtzLUOqQcLlcuN3uoM/XjJRSTdz9998PwKBBg8K8EqWUUvWBlvOppqK2n3UNpJRq4l555ZVwL0EppVQjFzl4MADuZcvCvBKlQkcDKaWauLZt24Z7CUoppZRSDY72SCnVxC1ZsoQlS5aEexlKKaVUnRg8eDBdu3b1ed/KlStxOp28/PLLPu+fOnUqTqezytsJJ5xQ7Wvv2bOHCRMm0KVLF1wuF506dWLcuHH8/vvvNa47JyeHO++8k27duhEdHU27du04//zzWb16dc1vWh0SmpFSqombPn06AEOGDAnzSpRSSqn6x+l0smDBAp/3xcXFVXleVlYWJ554Ik6nk4kTJ9KqVSs2b97Ms88+y4IFC1i1ahU9e/b0eW5xcTEDBgxg+/bt3HTTTXTq1Ilt27Yxd+5cTj75ZN577z3OOuuskLw/FTwNpJRq4ubPnx/uJSillFL1VkREBKeffnrA57322mts27aNbdu2kZKScvD4+PHj6dKlCzNmzODVV1/1ee7ixYtZvXo133//Pccdd9zB4//617/o2bMn06dP10CqHtBASqkmLjU1NdxLUEop1Ug52rfH2rXr4PdRJftTmZQUCjdvDteyDomNGzficrlo3rx5ueNxcXHcf//9ZGdnV3nuhg0bgMp9zHa7nSlTpvD999+HfsEqYBpIKdXEffDBBwAMGzYszCtRSinV2JQNovw5XleMMeTk5FQ6np+f79f5vs4F2YcoIiLC533du3cnOzubkSNHMnHiRE4++WQiIyMBuPzyy6t9vaOPPhqA8847j8mTJ3PaaafhcrkAOPfcczn33HP9WreqWxpIKdXEPfbYY4AGUkoppWrHO+I82MfX5Wj0rVu3VsoM+evAgQNVnrtw4UKGDh3q876LLrqIH3/8kRdeeIH33nsPp9NJr1696NevHxdddBFHHHFEla/Zv39/7rvvPqZNm8aQIUOIjIykR48enHbaaYwZM4ZevXoF9V5UaGkgpVQT9/bbb4d7CUoppVSdSk1N5bXXXqt0/Mcff+Smm26q9lyXy8WHH37o875u3bpVeV5ERATPPvssU6ZMYdmyZXz11Vd89dVX3H///UybNo3Zs2czbty4Ks+/++67uf766/n000/58ssv+b//+z+eeOIJHnvsMe69917+/e9/V7tuVfc0kFLhtW0bPPQQrFkDP/4IeXmwcSN06BD4cw0ZAp98ApMnw/33lx7/5BN5jXXrIDMTWrSA3r1hyhSo5h/AkNi6FSZOhGXLwBgYNAieeALatav+vEB/X776Cu69F9auhfx86NwZrr0WLrusxiUG+xM6pZRSqixfGSVvT5S/j68rDoeDPn36VDpeVFRU47k2m83nuTUpKirCGEOLFi248MILufDCCwH4888/GT16NBMnTmTs2LEHy/3KKi4uxrIsEhISGDFiBCNGjABg+/btjBs3jnvvvZcrr7xS+5zDTPeRUuG1fj0sWABJSXDqqcE/zxtvSMDhy9690LMnPPMMLF0K06bBr7/CSSdBXTa65ubCgAHw++8wbx688gr8+Sf07w8HDlR/biC/Lz/9JAGa2w2zZ8M778Dxx8P48fDcczUuc+HChSxcuDCAN6aUUkqpmqSmpvrMGnXp0oVrr72W/fv3k56e7vPcY489lksuuaTS8datWzNp0iQ8Hg8bN24M+ZpVYDSQUuF12mmwaxd8/DGcf35wz7Fvn2R9Hn/c9/0XXACPPAIjR0LfvnDxxbBwIWRnQ6BlbR06SCbLH7Nnw4YN8N57cM45cPbZsGiRBG/PP1/9uYH8vsyfD8XF8MEH8hqDB8vzn3giVLHBYFlPPfUUTz31lH/vSSmllAqAKTP225/jjclJJ53E22+/TV5eXqX71q1bR1xcXJVVIb1792bJkiU+A61169ZhWRYdgqneUSGlpX31kdstZV0vvQQ7dkCPHjB3rlz4n3aaZCtatw73KkPDFoJY/rbb4KijJGAqSZvXqFkz+eojnR4yixZJ1qtz59JjHTtCnz7w/vtQXU12IL8vhYXyPkqm+RyUmCiljDV4//33/X8tpZRSKgDeEefewRKHspwvGEuXLqV169Z079794LHi4mKWLl3q8/GRkZH0798fgFWrVuFyuTj++OMBmDp1KoMGDeLEE0/kqquuom3btqSlpbFq1Srmz5/P9OnTD5b1fffdd+Tl5dG3b18A7rrrLt5//32OP/54rrnmGrp27cru3bv59ttvmTt3LldddRWtWrUC4OeffyYtLY0hQ4bU2e+L8k0DqfrG7YYzz5RyrWnTICUFJkyAu+6S3pdrrqk/QZQxkgmpiWVBFaNBa+3LLyXrUlVZX1nFxXLbvBluvx1SU2HMmKofX9X783igbE11Ve/v118lQ1TRUUfBW2/VvF5/XXqplPBdf730h0VHy/N/+qmUE9YgISEhdGtRSimlGrDhw4czduxY5syZc/BYfn4+w4cP9/n4Zs2akZaWBshI8/bt27OsJFj8+9//zurVq7nvvvuYOnUqe/fuJSUlhV69erF06VIGl5laeNttt7Fp0yY2bdoEQLt27VizZg333nsvTz75JLt27aJZs2Ycc8wxvP7665xfplrlscceY968eRhjQv3boWqggVR9M3MmfPaZBAi9e8ux1avhySflwn7u3LAur5yVK6XfpyZ9+8KKFaF/fbcbrroKbrkFDj+85sefeCJ4N7Dr3Fl+n1u2rPrxVb2/qVPl5lXV+9u7V3qcKkpO9itT5Lfu3eX1zz1XPj8gGapZs6oPFEu8+eabAIwePTp0a1JKKaXqiWXVZMH69u1bbi+pivtK3X333dx9991+vc7//ve/Sse6devG/Pnzazx3hY/riLZt25YL6Koyd+5c5tan68MmRAOp+mbWLDj99NIgCqREa/9+6c3xlqR16CC9Nz16lD6uXz+48Ubpx+nQAaKiypd7zZ1b+viFC2WynccDBQWS5Vq2TErKqruvrJ49JcirSVxcwL8NfnnoIZlmN3myf49/5RXIypK+pUcflV6iL7+sehKer/c3fDicdRZceWXpseren2VVPhbqnxj9+Secd55kumbNkj/z99+Hf/0LnE646KJqT3+uZCCFBlJKKaWUUv7TQKo+2blTJrxdfXX544WFksWoYZ+DSt58s3ygVfZ1rrxSgoSOHeXYf/8rF/3V3VdRbKzv56/I17m1tWULPPAAzJkjwV5BQel9BQUygCIurnzJ3ZFHytcTT5TyyQ4dYPp0CT58iYuDihveORwSWPqzEV5SkmSlKsrM9J2pCtadd0oG6sMPS3u+Bg6EjAy44QbpHaum5+rjjz8O3VqUUkopH+p7b5RSwdCpffXJ+vXy1RvAgPTovPwy/O1vocvs7NghAUaLFqXHjjtOAp7q7qto5Uq5cK/pNnBgaNZd1oYN0jM2dqwEJd4bSLYpKQl+/rnq8xMTpbzP+3teF446SvqkKlq3LrT7V/38Mxx7bOXBGSecIMHU7t3Vnh4dHU10dHTo1qOUUkop1QRoRqo+8WYNymYxnntOLryPO67y40ePLl+6VzEoqHj/ihUSQBx7LJxyCrRvL1MAe/eWaXdt2lR/X0XhLO3r0QM+/7zy8f79JbgaP778tLyKdu2S7F8NZW+VlDSB+mX4cOnf2rABOnUqPf+rryQTFiqpqbIRb2GhZMy8vv1WSvuSk6s9/dVXXwVg7NixoVuTUkoppVQjp4FUfXLssZCQAPfdJ1+3bZPR3uedJ3sELV4sgYJ3l/CKpXv9+pV/vqpK+2w22bT1998lq7R4sZTJrVkjwUd195Xlq/QtGN69nLyDIBYvloxYixYyyAFkLQMHykj4f/5TAsKK79erffvy9517rgSixxwD8fHwv//BjBlgt8PNN1e9rqwsCWJrEh/vO8N0xRWyCfDZZ0vPmWXB3XdD27YyJMOr4nvz8uf3BeDaa2WvqWHDZMKjyyWj1994Q/bXKhtc+eBtZNVASimllFLKfxpI1ScxMTK2euJEGDVKLphnz4YzzpCeoLPOkov7UDniCLlddRUMGSIX394+rOruC7WKG85OmCBfy07D844i93gCf/6TToIFC+CxxyRr07atBFp33FH1oAmQ3rDaTCWMiZHJgBMnyibAxkjA9MQT0l/mVdV78+f3BWSj4Y8/luEbl18uJY9/+xs8+2z5gK0K1U0zUkoppZRSvmkgVd8MHgy//FL5+Hffhe410tKkxKxPH/k+MxM2bpSL7+ruqyv+TLHr18+/x/l6zKRJcguUv69ZnXbtJMMXzOsE8tpnnim3IETW5abESimllFKNlAZSjVnFHqmHHpLsVlGRlA9u3CibtxYVwSWXSAna5s1V36caJe/eE5deemlY16GUUkop1ZBoINVQ+Rp6ULbcq7qhCO3bwyefBH6fapQ0kFJKKaWUCpwGUko1cb52U1dKKaWUUtXTQEoppZRSSinV4A0cOJD09HQsyyIuLo6nn36aHr4mWIeIBlJKNXGzZ88G4IorrgjzSpRSSjVWgwfLYKNly9xhXolqzBYuXEhCQgIA7777Lpdeeilr166ts9ez1dkz13OWZSVblrXe3+NKNVZvvvkmb775ZriXoZRSStWJyy+/HKfTWeVt6dKlh3xNlmXx6KOPVvuYfv360a+qPTN9GDduHJZl8cADD1R6Hsuyqr2Fuk9627ZtXHfddZx88slER0djWRabqujf37p1KyNHjiQhIYH4+HhGjBjBli1bgnpdbxAFkBXKLYOq0CQzUpZlPQicCzj9Oa7qxrZtMkhwzRr48UfIy5NhgdVt7eT1ySdy7rp1MqG9RQvo3RumTCm/N25tXqO2tm6VLaSWLZNJ5oMGyRZS7dqF5tx+/WQvX1/OOAOWLPFvncuXL/fvgUoppVQD1bJly4Mb0Fd03HHHHeLVwOLFizniiCNC9nz5+fksXLiQyMhI3njjDSZPnnzwvocffpi9e/ce/P6SSy4hJSWFhx9++OCxNm3ahGwtAOvXr2fBggX07NmTU089tcpgNTc3lwEDBhAVFcW8efOwLIu77rqL/v3789NPPxETExPwa1900UWsXLkSm83Gxx9/XNu3Uq16H0hZlrUA2G2MubbC8euAG4FWwDpgsjHG33FzXwDrgfv8PK7qwPr1sk9uz55w6qkQyA+E9u6V8yZMkCBqyxaYPl323v35Zxk+WNvXqI3cXBgwAKKiYN48sCy46y7Z3/enn2Sv3tqeO3Nm5f2Z/+//ZN/k4cPr7r0ppZRSDY3T6eT0008P6XO63e6g92IcMmRISNfy/vvvk52dzb333su///1vfv75Z44++mgATjjhhHKPdblcJCcnB7SGnTt3smHDBnr37l3pvry8PFasWMGZZfazPO2009i1axcAc+bMqTKQmj17Nhs2bOCPP/6gc+fOABxzzDF06dKF559/nptuugmQYLeqLNUPP/xA27ZtD37/2muvHXzdSZMm8dFHH/n9PgNVr0v7LMs6FRjm4/jVwBPAm8CFwC/AR5ZlnejP8xpjFgOf+Xtc1Y3TToNdu+Djj+H88wM794IL4JFHYORI6NsXLr4YFi6E7Gx4++3QvEZZHTpItstfs2fDhg3w3ntwzjmyDdeiRbJN1/PPh+bcbt0kcCx7++UXcDhgzBj/1zpz5kxmzpzp/wlKKaVUI2OM4fHHH+fII48kISGB3r17V7oA79q1K3fddRc33XQTzZs3Z926dYAEaY8++ig33ngjKSkpdOjQgVmzZvHjjz/Su3dvnE4nXbp0Kfd8/pT2BeKVV16hb9++XHfddTgcDt54442QPTfAjBkzOP300/niiy/KHc/Ly+Oss85i3Lhx5OTkHDxus/kXYixatIiTTjrpYBAF0LFjR/r06cP7779/8Nh///tf0tPTfd7KBlFljR8/nmXLlpGRkRHIWw1IvQykLMu62LKsP4FVVC6/iwTuAZ41xtxpjHnPGHMpsBa4+5Avtg643XD//dCpk+yne/LJ8McfUp4WHQ3bt4d7haHh598xvzVrJl/L/nAo1K/hr0WLJLAp8+8CHTtCnz5Q5t+FkJ6blwdvvQXDhkFysv9r/eCDD/jggw/8P0Ep5VtREezYAYWF4V6JUvXKm2/a+O03iy+/tDjuOAdvvnno/3M2xpCTk1PpVljy9/XOO+/kvvvuY+zYsbz88st0796dkSNHVvr/cc6cOfzxxx/MmjWLDmX6BB555BEyMjJ49tlnOfroo7npppsYNmwYffv25bXXXqNly5ZcdNFF5Ofnh/y9paens3TpUsaOHUtiYiKDBg0KeSD1wAMPMHToUIYOHcpXX30FSBA1bNgwfvzxR5YuXUpsbGzAz/vrr7/SvXv3SsePOuqog4GqvzIzM9mxY8fB79955x1atmxJciAXRQGqr6V9G4EXS349qcJ9fwdSgNcqHP8YuNWyLDtwE3Cbj+cdaIz5MdhFWZZl/Hlcz549g30J3G4480wp4Zo2DVJSpHztrrsgPx+uuQZatw766UPGGCgurvlxlgUREXW3juJiuW3eDLffDqmpgWVjfKnqvXk8cp3kVd17+/VXySRVdNRREuxUJ9hzvRm5Sy6p/vkrWrx4cWAnKKV8Ky6Wps2sLGjZEhIT5R8KpZqwN9+0MWWKnQULLE45Bb78EsaPtwNFjB7tOWTr2Lp1K82bN690fNKkSVx22WU8+eSTPPvss4wbNw6As88+m127dvHggw8ybFhpcVR8fDzvvvsuDoej3PN069aNV155BYDOnTuzbNkyzjnnHKZNmwZAcnIyAwYMYN26dSHvyZo/fz42m42RI0cCMGrUKC699FK++eYbTjrppJC8ht1u5/XXX2fMmDGceeaZvPvuu0yfPp21a9fy6aefcswxxwT1vHv37iUpKanS8eTkZDIzMwN6rszMTEaPHk1+fj42m42WLVvy4YcfYtXhv8P1MpAyxnwJfAlgWda/Ktzt/ZP6pcLxX5HsVVtjzMPAwzRAM2fCZ5/JPzTeMtTVq+HJJ+UCf+7csC7voJUrpWenJn37Ql3u93riifD99/Lrzp3l965ly9o9Z1XvbepUuXlV99727gUf/y6QnCzXWdUJ9tyXX5b3XqZEWSl1qNntUjqwe7cEVKmp0vCoVBPgHXFe1m+/WSxYYB38f7V/f3jxRYtRo+zMmVP+59N1ORo9NTX1YO9MWW3btmX58uV4PB4GDRpUrjxt8ODB3HrrreTn5+N0SoFUv379KgVR3uNeiYmJAJx4YmnHSVxcHADZ2dnlzisqKiqXpXI4HD6fvzqvvPIKgwYNwhjDvn376Nu3L3a7nTfeeCNkgRRIMDV//nxGjRrFoEGDSEpK4tNPP+XYY4+t1fP6CnSM8St3UU6nTp1YvXp1rdYSqHoZSNUgGSg0xhyocNx7iRlfVy9sjPErpO3Vq1fgf/olZs2C008vDaJAfqi5f7/06HjL1zp0kB4aX3uMHTgArVpJiZf334ziYhm4AFJ18scfUNKDSOfO5fuK/NGzpwR4NSn5d6POvPKKXKts2ACPPgqDB0sQWpupfL7e2/DhcNZZcOWVpcdqem++fgDi778LgZ67fTssXw433CDXcYF48sknAbjhhhsCO1Ep5ZvNJv9AFBTApk0yEScxMXy1xkqFUUYGnHJK+WOnnCLHDyWHw0GfPn183rd7924AunTp4vP+vXv30rqkHCg6OtrnY+LjK19+ulyuSscqBgivvvrqwSwYwD333MOUAJqy//zzT7777juASpmdBQsW8PjjjxMRwtKgoqIicnJysNlsuN1ucnNza/V8SUlJ5SYKemVmZvrMVNU3DTGQqmk8Sp4/T2KM2QQc5u/xQ2HnTvj9d7j66vLHCwslG1EyuKRGb74Jxx4LH30kGYykJClB8+5HtmmTBGC12Z8sNtZ3EFdRXVe1HHmkfD3xRMnEdOgg0/tmzQr+OePioFev8sccDimprHi8KklJklmqyPvnEepzX31VSg8DLesD+PTTTwENpJQKuago+ccjPb00O+XU3TVU4+Uro3TccQ6+/LJ8pceXX8r/3/Vlc974+HhsNhtLlizxGXQ08/4Uuw4MHTq03ACHdv7skVLGq6++SkxMDIsWLSo34GHNmjXceuutrFixgoEDB4Zkrfn5+Zxzzjn897//5ZtvvuHBBx/kzDPP5JNPPuHkk08O6jmPOuoofv3110rH161bR7ey+9nUUw3xx2OZgMOyrIr/G3l34Eo/xOsJmfUl2wB37Fh6rLhYSrb+9jf/szsvvgg33yx7CfnIYofEypUy1KGmW4j+7volMVGya+vrwXbKRx0lvU4VrVtXfp+rUJ378ssSPAeTXV+0aBGLFi0K/ESlVM0sS37yBNLMmZ4uP/VQqomYNKmI8eMNn38ufeCffw7jxxsmTSqq+eRDpE+fPng8HiIiIujTp8/B28qVK5kxYwZRdVie27JlS0455ZSDt0ADqddee42zzjqLAQMGHNzAt1+/fkyYMAGXyxWyoRMFBQWMGDGC1atXs3z5co4//ngWLFjAwIEDGTJkCN98801Qzzt8+HC++eYbNmzYcPDYpk2b+OqrrxjeAPZyaYgZqd9LvnYH1pQ5fgSwzRjj42f5DYP3BwllsxHPPScX0P72Jf72G/z1l5ShRUfDpElw7bU1nxeo+lLaV9auXZLRu+ii0D93FZtxV2n4cLjlFik57NSp9Dm++koyZqE8d80aCbwefzywNSqlDiGHQ3665B1G0aqVjGVVqpGTgRJFjBplJyNDMlFTphzaQRM1OfbYYznvvPO4+OKLufXWW0lNTWXVqlW88MILB8vfwyUtLY1nnnmm0vEePXpgs9n466+/eOihhyrdHx0dzeDBg3nnnXeYOXNmwH1XFU2ePJlvv/2W5cuX06OkJCkyMpIFCxYwatQohg0bxsaNG8tN7nu7pG/k+5Jm9sWLF9OiRQtatGhB3759Abjiiit45plnOPvss7n//vuxLIu7776btm3bctVVV9VqzYdCQwykvgb2AWMoCaRKJvWNAupux61D4NhjISEB7rtPvm7bBrfdBuedBx98AIsXS2q8usqQF1+UPZXsdhg0SIKy778v7Y8KFV/lb8Hw9mZ5B0YsXiztBC1ayDAHkOzXwIHw0kvwz3/KsXPPleDymGMgPh7+9z+YMUPe9803B/4aZWVlSfBak/j4qjNEV1wBzzwj0/fuv19+KH333dC2LZT9d8HXe/P3XK+XX5b3feGFNa/ZF+8+FrfccktwT6CU8o9lyY7ahYWSnWrWTOq263K0qVL1wOjRnoODJepLOV9F//nPf5gyZQrTp09n//79dO3alZkzZ3JJMDXzIbR+/Xquu+66SsdvuOEGCgsLiYmJKbcRblnnnnsuixYtYvHixZztaxxwACZPnswll1xycJNfL28wtXr16krjz8+vsIHnhAkTAOjbty8rSqZ1xcTE8NlnnzFx4kQuvvhijDEMHDiQJ554Iqhx6oeaFcxUjEPJsqxNwIfGmGvLHLsFmA5MBX4GLgP6AccYYzb4eJpDqlevXmbNmjU1P9CHZctg4kQJDFq0gIcflhK9oUMlEMjKkv+HfQ2bcLvhsMPkotq7l9K+fXKBXXa/VW+P1L59Qb/FkKmqh6rsRLwVKySA/M9/4NJL5dhDD8GCBZJ9KyyUIKNfP7jjjsqDJvx5jbK8r1eTmiYSbtkif5bLlsmgiIED4Yknyq/P13vz91yQP/PWrWXfqWC3gjrvvPMA2W9BKVULBQUSIPnzn78xkJsrQVRqqpQQNEHr168nMrKm1md1KGVlZdG1a9eQP693ol99DaTqgjHm4LQ/VX/99ttvHOltuvfBsqzvjTE+0wcNNZCygFuA64DmwA/ATcaYbw/x2qYgmwOX06pVK7bX8a65vgKpd96R0q+yJXcbN0rmZvv20iqS+hRIKaVUo7F9Ozz2mGz+52/Zntstu2knJUHz5k0uO6WBVP2jgVToaCDVMNQmkKr3wyaMMR3KBlElx4wx5hFjTDtjTLQxps+hDqJK1jHFGGNVvLU+RDvmnnGGZKC8tzlzpKyvrI4doXv3mjeBVUopVUv9+kmzYr9+8OGH/u13EBkptdLZ2fJTrgMVd/ZQSilVXzXEHilFYMMPykzVBCSbpdko5TW9ZILF7bffHuaVKNXApaZKRlzwk/cAACAASURBVMmypKHx5JOl6bWmUZ2WJaV9RUWwdas0ybZoEfimcErVY00pE6WajnqfkVJK1a21a9eytjabiimlSsXHw5Il8OCDMkZ1xgz/z7Xb5fwDB+SnZdnZdbZMpZRStac/7lKqiZs/f364l6BU42K3y+7Yw4fLNByQDe5WrZLxnDVlmrzZqbQ0zU4ppVQ9phkppZRSqi4kJUFKivx64ULZx+D00yvXW/tSNju1caOMbK3nw6GUUqqp0UBKqSZu6tSpTJ06NdzLUKpxu/VW2egvLw/GjIHLL5dR6TWJjpbNA7dvl5tb+0xU3avvE52VCpXafta1VkCpJu6PP/4I9xKUathSU2HXrtLv27SRry1agLf/0LJgyBCZ6PfCC/D009C+vWSpahIRIdmpvDzJTqWkyPdVbZJXnX795Gt1m+CpJs1ms1FQUKBju1WTkJeXV6stGOr9PlL1WTj3kVJKKVVPVBfQpKX5Pr5jh+yuHh8P33wjjxsxoubgyOOBnBw5NyUFHI7A1lrPAindR6r+KSwsxO1206ZNG6KiorCCCdgVoPtI1WfGGPLy8khLSyMlJYX4+PgqH1vdPlKakaoFY8wUYErF47169dLoVCmlGjtvUFKdkSPh7bcrH2/VqvTXr78uO6rPmyfj0svutF6RzVY5O5WQEFx2SikfHCXB+bZt2/B4PGFeTcNmjNEfFNRjkZGRNQZRNdFASqkm7t///jcA9913X5hXolQT9cQT0KcPTJsG//gHjBoFd9wBLVtWfY7LJdmpXbtg/34pL4yKqvrxFcsPvYFXSgrs3Bma96EaDYfDcTCgUsFzu9107tw53MtQdUiHTSjVxG3dupWtW7eGexlKNTwrVtRcIucrG1WRzQajR8s0v6uvhnffhY8/9u+8uDgoLpZ9pzIzq57sVzaI8ue4UkqpGmlGSqkm7j//+U+4l9Ag5eRAfj4kJuoWPypE4uLgrrvgoougbVs59tFH0gc1aFDV5XtOpzxm927JTrVqVT475U8JolJKqYBpRkoppYLgdktF1MaNcu2qrQRNmHevqIqaN5d9oAL9cHTsKNG5MfDSS3DppXDxxbKpb1W82Slj5EO5d69+KJVSqo5pIKVUE3fHHXdwxx13hHsZDVJUlLSq7NwplVUHDoR7RSosdu6UAKZvX7kZI7cdOyA5GXJz5cMR6JRcy4L58+Gee2DNGhg4EKZMkci9KlFRElClp8s+Vfn59WZCn1JKNTYaSCnVxGVkZJCRkRHuZTRY3kRARARs3Sq3/Pxwr0rVC3Y7NGsGnTpBUpIEU7m5gQVUkZFw5ZXw5ZfSRzVnDnz/ffXnWBbExsrXTZskqFJKKRVyWtmvVBP3wgsvhHsJjUJkpNzy8+XaNSlJkhE6+VZht0uZX2Ii7NsnZXc2m6Qz/R1b3rw5PPww/OtfEpgBvPoqdO0KJ5zg+xyHQz6Ae/fK+b4CqqrKEpVSStVIM1K1YFnWFMuyTMWbbsarVNPldEqGKjtbW1VUBd6AqlMnCapycgIv+fMGUYWF8NxzcO65cM01VW/8681OrVkD69ZB795w2mml5Yc6+lwppYKmgVQtGGOmGGOsirfWrVuHe2lK+e2WW27hlltuCfcyGhXLguhouWVkSECVnR14i4xqpCoGVMGU/DkcsGwZTJwIS5ZIcDRjhmzU60tkpGzkW1QkadPc3NC8F6WUasI0kFKqicvLyyOvqosvVSs2G8TEyDVsWpr0/utvtTooMhJatJApffHxkqEKJKCKjoZbbpFhEoMGSSC1eXPVj7cs2aNqwQLYskX2kCouDslbUUqppsgy+iPSkOvVq5dZs2ZNuJehlKpDO3dCVpZUTfmrsFCSAQkJMoPA4ai79akGyO2WTXUzM2V6SSA9VCDNeR06yK+feQYGDIBu3ap+vDcr1aqVRPxhsH79eiK1kVA1Um63m86dO4d7GaqWLMv63hjTy9d9mpFSSqkg3HcfXHABfPut/+c4HJJ4yM2Vcr/0dE0IqDIiI6FlS8lQxcVJhiovz/8MlTeI2rsXZs2CM86AO+6Q732JjpYP5datMqq9qCgkb0MppZoKDaSUauJuvPFGbrzxxnAvo8Hp1g22bYMRI+Cf/4TffvP/XJdLMlmZmaUb+mpxgDrI4ZBpeh07ygclOzuwmtDkZPjiC9nI97XX4NRTZWNfX4GS3S7R/YEDpc18Siml/KKBlFJKBeGii6TX/447YPVqGDwY5s71/3zLkmqqqChpVdm4Ufv/VQXegKpTJ/mwBBJQJSXB1KmwdCl07y6j0/ftq/rx0dEycjItTW5ud2jeg1JKNWLaI1UHtEdKqcYvM1NK82Ji5NczZ8p+qZ07y3Wo0yl9UP5yu+UaOTZW5g9ERdXd2lUDVVAgZXpZWZJJcrn8O88YKd9r105m8T/0EFx4IbRv7/vxeXnyuJQUKTEMpE8rQNojpRoz7ZFqHLRHSiml6lBSEkyeLEEUwJQpcPLJ8Pjj0ubiD+906sJCmRmwe7e2rKgKoqJkMESHDhJEZWX5l6GyLAmiAP78U8r8+veH6dOlpK8il0t+ErBjh2anlFKqGhpIKdXEXXPNNVxzzTXhXkajMmkS9O0Ljz0m+5++9JIkE/zhdEpWKisLNmyQbJdu6KvKiYqC1q2lhyqQgArg8MNh1Sr4xz/g6adl/6l33qn8IYuIkGxUQQH89ZeMqdQPolJKlaOBlFJNnMvlwuVviZDyS+fOMHs2fPihXLfefbeU/vmr7Ia+e/ZIhko39FWVeAOqDh0kAs/Olvn6NWnVSoKo99+H1FR49FFJhfri/bfhjz9KJ6NoQKWUUkAte6Qsy3IAqUA8sBfYbYxpMsUolmVNAe6peLxVq1Zs37790C9IKXXIlO2Rqo4xsHIlHHOMDFNbvVoSCAMG+N96UlQkCQeXS6ZjO521X79qhPLzISNDAiqHw78PiscD27fDYYfJh+yRR+Bf/5IPmldWlnzY27eXiSh2uzTyxcbKrtO1oD1SqjHTHqnGIaQ9UpZlJVmWdZNlWZ8DWcBG4EdgK5BrWdYqy7JutSwruVarbgCMMVOMMVbFW+vWrcO9NKVUPWFZ0K+fBFEAL7wg49JHjgR/Z9LY7VJlVVwMmzdL64q2rahKnE5o00YyVA6HBEA1ZahsNgmiQD6QL70k49Kfe65ylspb7udwyIdw40Z5Dc1QKaWaKL8DKcuybJZl3QRsAq4E/gdcD5wPnA6cB1wN/Az8E9hoWdYdJVkrpVQ9deWVV3LllVeGexlNxrPPwgMPSNvJ2WfDZZfB//7n37lRUXIdm5sr/VMZGbqhr/LB6ZTgqH17mWKSne1fk96pp8Knn8JJJ8H990vadNmyyjWl3oAqMlJ6pzZtkoBKa0+VUk2MPYDHfo8ET2cYY76p5nEvAliW1Qu4GVgLdAt6hUqpOtUskBnd6iCbTbJCxcVyXekvh0P2ST3/fJgzR37w/9VX0LWr/8/hcsk1a0aGlBi2bFnnU6pVQ+RyQdu2UrKXni7BTlRU9bP1//Y3mDcPPv8c7rlHIv8TTvD9WLtdyvuKiiRDlZ5eWvKnH0alVBPgd4+UZVndjTG/BPwCltXVGOPnz1sbB91HSqnGrV8/+fr++3LtWFws16yBBFRee/fKdafDAQsWwG+/wXXXlZYC1qS4WDJUUVESUEVHB74G1UTk5cn0ktxcyVo5aigYcbslWo+OlmETn3wiH86EBN+P9zbzRUb6HVBpj5RqzLRHqnEISY+UryDKsqyWvh5b4bwmFUQppZqOhASZQJ2SIu0kOTmBl9olJ5dez/71l2SpeveGJ5+U692aeKusALZskW1//B21rpoYl0v2k2rXTlKq2dlVT+sDCYhSU+XXq1fDrFlS/vfaa74/6N5mPrtdBlhs3Ch/KbTkTynVSNV2/LmmXZRq4MaNG8e4cePCvYwGy2YrH1AVFAQXUAHccQcsXy6B1MMPQ58+MkLdHw6HbOibny/Xr3v26Ia+qgrR0VLy17atfJ+VVX1ABTBkCCxeLKV/t90GQ4fCd9/5fqw3oIqIgG3bpIfqwAENqJRSjU6NPVKWZQ2v5m4dwqtUA9fWezGlasUbUMXFlU6L9njkmjWQCdGHHy6D01avhmnTSrNVBQWSIKjpuVwuqdravx/27ZMKq/j4Wk+pVo2Nd7Oy9u2lHG/3bvngVlfyd/TRsHAhLFoEU6fCK69U3T8F8oGNjJQSwW3bpP60RQt5Xe2hUko1AjX2SFmWVQysBHz9q3eSMUZ38qxAe6SUaty8PVIrVlT9mOJiuS7NyAguoILSH+BbFjz0EHz2Gdx5J5x2mn/XoR5P6bY/KSn17PrVn99EdegYIx+WPXskanc6JQiC0p8MdOpU+vjcXEl/JifDr7/C0qWy/1R1m3u73XJOmYBq/V9/aY+UarS0R6pxqG2P1J/AeGNM/4o3ID2kK1VKqXrujTfkunHVKujeXb73JSICkpKk5K9FC/mhf05OYFvuWFZp4HPkkZJluvBCGDUKfvih5vNtNun3t9th61a51bStkGqiLEt2l27fXkanGyM9VFVtWBYdXToRZflyePRRCY4//LDqEr7IyNKGvq1bYfNmrLw8LflTSjVY/gRS84DmVdw3K4RrUUqFwdixYxk7dmy4l9EgvPEGTJ4s0/UKCuDpp+X7qoIpKA2oOnUKPqACGD4cVq6E++6TAWpnnSWTqf1ht0t5X3GxtKvs2qUb+qoqVAyoiovlQ1udG26At96SIOmqq2S2/7p1VT/e29AHROzYgS0tTQIqpZRqYGrskTLGTKvmvgdCuxyl1KF2+OGHh3sJ9Yq34syXX3+VIKp/f/m+f3948UXJED3/fNXnrVhRGlDFx5dWSoFUQvlb8hcVBePHw+jR8MILpWvdsUO+tmpV8/lRUZJo8PZPJSQEN7ZdNXLegKpjR/nAbtkiAyOqqg/t3RuWLJGJfg8/LJmpbjVsIelwYGJisADb9u14nE5Ms2YYp7ZfK6UaBr/3kVKVWZY1Bbin4vFWrVqxffv2Q78gpVStVRdIrVpVOvTBy+2W4OS006o+z1cbUHGxlOplZMj3gQRUFd1wg1y3XnYZXHMNJCbWfI63fyoiQgKqQ76hr/ZINSzFxbLpWUZGzXtQZWbKX4roaPnz3bABLr64/F+cEhs3biztkSosxCoowONyaUClGgXtkWocquuRCjiQsiyrjTEmLSQra6R02IRSjVP37lLO581IAXz+uexR+kvA25ULb0CVni6BTDAB1ZYt8Mgj8O67kvG65hoJqqrr+/fy7qEaFSUDKfw5p1ZSU6W2sKKUFNi5s45fXNVaXp6kQIuKJGNVU/R9222SperaVaL2xER4++2Dd5cLpLy8AVV0NCY5WQMq1WBpINU4hGRD3jL8aHFWSjUUY8aMYcyYMeFeRoMwebKU1n3+uWSiPv9cvp88OfjnjIiQnv1OneRrbm7gW+60aycB3tKl0LMnPPggPPZY+cesXi23irxb/gBs3iz7qNa0pVCt+Aqiqjuu6heXCzp0kA9rTk7NH5aHHpJ5/gUF8Ntv8Pvv8kGrjsOBiYvDKi4mIi0N2/btOiVFKVUv1dgj5UN9GZ6rlAqBHj16hHsJDcYFF8jXUaOkwqlbN3jggdLjtWG3Q7Nm0rO0b59UUdlsct3qb8ldt26ytc8330hrC8DatTIgLSWl+udxOOSWlycb+jZrJj1dteqfqq5OsrrHa7lf/WazQfPmMhJy504JqKrKTlkWnHEG9O0rqdxt2+D//k+GWdTE4cA4HFgFBdjT0vBER+NJSpLSQqWUqgeCCaS0qUqpRuT2228P9xIalAsuKB0sURfX+3a7XKMmJgYfUJ10Uumv582TARlHHCHXssXFcg1b1WAKl0uyYZmZcmvZUjJWuqGvqsTplHTovn2y/5TDITWiVT22TRtpyBs1So699Rax6ekU/OMf1X/AoqIwUVGlAVVMDJ7ERA2olFJhF0yP1G5jTMs6Wk+joD1SSjVuhzJxUlQk16kZGZIdCiSgAkhLk96pP/6Q6qq//Q06d5YhazVN+SsullLDyEjJaMXE1O69ANUvXocfNVwFBZKdys+XD0rZwKhHDwm0KoqMBLebgu7d2XfzzbhrmvJX5rUst1sCqqSkqoM3pcJMe6Qah+p6pILJSCmlGpHzzjsPgHfeeSfMK2m8vv669s9RVCQVVPv3yzWqvz+M//13mSh47LESTP3yi2S8li+XLJU/iovlOtnplHK/6ga2ldW7t3+PU41AVFT57JTdXvoh9RVEAbjd7LnzTpJnzSJl3DgO/OMf7L/mGjzNmtX4WiYqCis/H/u2bRpQKaXCRgMppZq4k08+uU6e1xs86MV0aNjtUu4XGyv7QGVnS0BV07Wj2y2T/HJypIeqb185/sUX8OWXEliNGiWzA6oSESGTrAsLZWBbXJw8pz2Y/0FSUqqe2qcaNsuSSDsmRv6Ms7NrTGPmDB1K4cCBxP/nP8TOn0/u0KEU1BRIeTmdGCgNqGJjpeRPAyql1CESzH+DW0O+CqVU2Nxyyy3hXkKDE2hJX10Ek2639DDt21d9D1VxsWSgvJVWXbqUDsr47jvJTK1cCZdfDhMmSIBUHWNkIIXHU7qhb0D9U94R5zpYovFyOOCww2QjXz+mMZrYWPZfdx3ZF1yAp3lzAGJffZWi9u3JP+WUmmtZywZU3h6qQFKnSikVpIDbh40xPetiIUqpxiMtDX76SS7UV6yQ71VoRUbKIIiOHUszTrm5lduM2reHn3+WgMvjkSDq55/l+PXXy5/P6afL+PSTT5bBFNWxLMlORUdLxdbGjZJ40PYmVY5lSZTtHR/pB28QRVER0YsX0/yWW2h+443YN23y7wmcTkxMDLaSDJVt9+46nuWvlGrqajWHybKsxDK/HmhZ1jjLsg6r/bKUUofK8OHDGT58eMieLy1N+nKOOQYGD5ZNbH//XYOpuuINqDp18h1QtWol2/78+CN89pmU8nXoUDpoomNHmDkTPvlEZgJ4q6IKCqQvqyo2m5QZRkbKn+2WLZKpUqqcipvt+sNuZ/e8eey78UYcP/9MyoUXkjBjBlZ2tl+nG6cTEx2NLS8P+9at2Pbs0YBKKVUnguqRsiyrA/A+8BcwwrKsiYB3+8dsy7IGG2O+C8kKlWoEQjFsoK506jQQCN0af/pJgqisLLl5y84++ECO13cNtafLG1AlJpaW/Hn7/Vu1Kh0scfzxvs/v3h1ee600AHv+eVi4EG6/XbYBqqq6ym6XAK6gQPZZTUiQPai0qkodVFVfXIsWVZ9jt5NzwQXknnEG8bNmEfvuuxw47zyKvLtH18SyME4nGIMtNxdbdjaeuDg8CQn64VRKhUywGakngSjgwZLvJwIzgGTgO2B67ZemlDoURo26gVGjbgjZ8+Xny8V0WQkJclzVPYdDrls7dpSMUU6OZIp69ao6iCrLGzAdeaSUAo4fD8OHyx6q1YmKkoAqN1fK/dLTpT/Ll379oN/aJwJ6X6oB27lTIvS+feWWlQX/+59fP73xJCez78472fHuuxS1awdA/LPP4vjhB/9euySgMtHR2A4ckJK/9HRpMlRKqVoKdmpfP+BGY8way7KOAtoAzxlj9lmW9QIwK1QLVKoxaKhZjmAUFsqFvHcAwZFHykV1jx5N6/ch3LwBVVKSbOq7b59krVwu/84fPBj695eeqcceg5Ej4brrJENVnbIb+u7bJ0mH+HgfGa0ePXTQRFMVFycflD17sHJy5PsaRkB6R6Lb9u8n+pNPiH/5ZXIHD2b/dddR7M/ER8vClHw4bTk52LKy8MTHS4YqmPJDpZQi+IyUB9hb8usBQJoxZn3J91GAn/9VK6XC7cwzz+TMM88M2fN16SL9OJmZktFIT5fvu3QJ2UuoADgckJoqPVQxMTIYwt9eJrsdLrxQxqTfdZeU+IEMmdi8uerzLEteKypKkhGbNkmmqjHr1690EKHyg90OrVpRnJoKbjeWnx9KT0ICuxYsIGv8eFyrVpFy/vnEzZnjf8q7JKAy0dHYcnKI2LIFW0aGZqiUUkEJNiP1DXCNZVkHgBuBjwAsy3IBVwA/hWZ59ZtlWVOAeyoeb+Xt4laqARg2bFhIn69NG/n6wQdybdOjh/TneI+r8PAGVN4MVVaWXMv6k6FyueDqq0u/f/xxeOMNuPhiuOEGGa/uS0SEJBvcbhlGERtbfVuMauR8ZCBNTAzFcXHY9u7FlpUlfU01ZKeM00nWlVdyYNgwEp55htgFC8gZNUrO9VfZDFV2Nta+fZikJDxBb5CmlGqKLBPEzNqScr4lSEnfbuBkIA3YD+QD5xljPgvhOhuUXr16mTVr1oR7GUqFlW7IW78VFAQeUHnt2AEzZsD8+ZJ1uuoqudU0ByA/X0o/x4+X11y5snbvoT7RbbGCs379eiJLSuusvDxsu3ZhGSNBTk37R5Ww7d2LJzkZiotJfPRRDowYgTvQFLgxWPn5GI8Hk5ioAZUKCbfbTefOncO9DFVLlmV9b4zp5eu+oEr7jDG/Ap2B44DOxpiNQBFwGXB4Uw6ilFKqIYiKKh2N7nJJQOVvyV+rVvDwwzJOfcAACaoefLDm85xOCbaKiiSo8pZ/KgVgXC6K27bFExeHdeBA9fP3y/AkJwNg37IF16ef0vKf/yTxoYew7dvn/4t7M1TR0diysqTkb+9ev9eglGqagv5xizGmAFhb5nsP8EYoFqWUOnQGDRoEwPLly0P6vJqJahiioqB1axlZnpEhAZXDIUFPTTp3ljHpa9eWluz9+iusWwcjRkhpX0WWJceNkV6rzEw5NzbW7wSEaswiIvA0b46JjZXsVEEBJjrarw9HUceO7HzrLeLnzCH2nXeIXraM/VdeyYERI/zPLpUt+du/H2vfPjzJyZi4ON8faKVUk1arDXmVUg3f6NGjGT16dLiXocLMG1B16CC/zs72v3+/R4/SHrj58+HGG+H002HZstJ9qSqyLAme7HbYvh22bq3bEflff12/93NT5RmnU7JTiYmSnfJzGIRJSGD/zTez69VXKTzySGLfeSe4BViWBHAuFxGZmURs2YK1b1/VM/2VUk2SBlJKNXFXXHEFV1xxRbiXoeoJp1OCovbtJaDKygoswLn3XnjuOTnn0kslM7V6ddWPt9ul3K+4WKb77dihA9RUCZsNT3IyRW3aSER+4EDVkXkFRZ06kf7UU+x5/nmw27EOHCDp/vuJ2LYt4DWY6GhwOjWgUkpVooGUUkqpSrwBVYcOUurnb0Bls8kGvitWwLRpEhyVHSrx3nuyF+s338DAgfI9lN/Qd8MGKTNsKNeqb7whJY2rVkH37vK9CiGnk+I2bfAkJUl2qrDQv/MsC09iIgCRf/yBa/lyUseMIX7mTKxA5/GXDaj27tWASikF1KJHSinVOPQrGTe2QseNKR+cTjjsMAmi0tMloIqKklt1IiPhn/+UjXy97r0XPvwQ3noLTjlF9qe67DK575xz5Kt3Q9+MDOmfatlSMla16Z9KS4OffiqdGtilS+jG8b/xBkyeLBsXe9/T+PFy3wUXhOY1FBLIJCVR5HIRkZ4OOTkQHV2683cNCo87jp1vvUXCs88SP28eMR99xP5rriH3zDMD+3DZbJiYGPB4iNi7F7NvH57ERO2hUqqJCmr8uaqejj9XDcncuXMBuPTSS8O6DtUw5OVJQJWbK5mq6gKqskEUwG+/wcKF0L9/6bHPP4fzz4euXUuPvf22fC0ulteJipKAKjo68PWmpcHvv0uAlpAAKSmyQXRt9zbzjjv/9VcJoiq+p1Gj4KijdBx6TcqOP/ebMVj790sgY7fXHNVX4Pj5ZxIffxxPbCzpTz1Vuyjd45G/FDYbnqQkCaj8DO5U46fjzxuH6saf15iRsiyrIzAH6Ai8B9xpjMkvue87Y8wJoVysUurQ0gBKBcLlgrZtSwMqfzNUAPv3S9amrFNOkf2syqrYU1VUJJmk6GhITJRsl79++gmOOUbWmZUl17gul2wYfcwx/j9PVTIyfL+njIzaP7eqgmVhEhMpcrmw7dmDFWh26uij2f3ii1jZ2WBZROzcSdxLL5F11VV4mjULbC02G8TEQHExtpI0qic5GRMbqwGVUk2AP6V9M4GFwDfAdcCnlmUNMcZkAwH+GEkpVd+4Szr7A/6psGrSygZUe/b4Dqi8mSWvgQOl9K1s9ubLL+Hww+H11+V7h6NyIGW3y62wUCb8xcdLuZ8/E63z8yUTlZVVeiwhofYTAr2Zpu7dfb+nbt00G1XnoqLwtGmDtX+/7PkUEeHf3H6QEr2EBAAcP/1EzEcfEb1sGVmXXUbO6NHyQQxERERpQJWeDiWbBGtApVTj5k8glWKMebbk15dalnUHEkwNBrQuUKkGbvDgwYD2SKnguFzQrl35gMrp9H0det110hP10kvle6RuvRXmzYP//Ed+ffbZVV97GiOvZYzsPxUfX/11amGhlPN5H3PkkZJJ69EjNHudTZ4sPVEvvli+R+qBB2r/3MoPJdmp4uhobOnp2LKzpYcpgOAl7/TT2XXEESQ89RSJzzxDzHvvsX/iRPIrphr94SugatYs4DUppRoGfwKpcgUbxphplmW5gU+BuDpZlVLqkLn88svDvQTVCHgDqtxcCVSysyU7VTag8g6UOP98Kec7/HAJnM45B774Qq4/r70WZs6EO+6QLE/F9hXLkioujwd275bnSUmRY75aXbp0kZ4ob49Uenppj1QoeAdKjBol5XzdukkQpYMmDjGHA0+rVpjsbAlgAslOAUXt2pHx6KNEffMNiTNm4Fy1KrhAyqtsQLV7t2w0rAGVUo1OjcMmLMt6F5hpjFlW4fhNwKPGGP0XoQIdNqGUasq8WaM9e+RrxQyVdwhFxdI/jwfefx8efhi2bJFs1dSp1b+W2y2vER0tAyl8XTunpUlPVH6+ZKJCObXPyzt8QhO7gQlq2ERN3G5sGRnYDhzAuFyBT9Mr9ElTegAAIABJREFUKsLKz8fExhL5669S8jd+vAySCFZRkXwA7fbSkr/aDLlQDYIOm2gcajVsAhjj66Ax5nHLst6s1cqUUmGXW7KfSnQwI9GU8sGbNfKW/O3eXX3Jn5fNBueeC//4B7z2mpThgWSd9u4FX9cjkZFyy8+XPauSkiA5ufxAijZtSgdLhKKcT9VzkZF4UlIwOTnY0tOxLEsCKn/Z7RLoAFFr1xI7fz7RS5aQ9a9/cWDYsODGnNvtEBsLRUWSoSpb8qcBlVINVo3ZJGNMgTGmwPu9ZVltytyXVlcLU0odGkOHDmXo0KHhXoZqhLwBVfv2ElSBBFQeT/XnORwwbhycdJJ8/+yzUuZ3yy0ybMIXp1MGUGRny4a+e/fW/DqqEbMsTFwcxW3b4nG5ZEJfUVHAT5Nz0UXsnjePonbtSJo2jZbjxuFYuzb4dXkDKrsd265dsrFvTo6kcZVSDU4wG/L+ALQM9UKUUuFx9dVXh3sJqpErG1CVJEApKpKyPH+quiZMkMe//LLsQzVunPRSJSX5fh2PR3qhMjNlIEVtKrJUA2e3S3YqJkayU0VFmAB6pwDchx/Onuefx7V8OQlPP03U2rUU9uhR63UdzFDt2iVZtGbNMFU1+yml6qVgAin9G65UIzJ69OhwL0E1EZYl/fdOp2y2a4xkkJzO6gOqZs3g3nvh8svh0Ufh+eclUHrySd+Pt9kOXqOyY4dkp/7+dxk4oZomExtLsdMpvVPZ2VLq58/8fC/LIm/wYPJPPRVTMizCuWIFkevXkz12bECDLcopE1BF7NiBiYqSgMrl0oBKqQYgmEERmn9WqhHZv38/+/fvD/cyVBMTESEZqsMOkwxSVpZkqKrTtq0ET8uXw803y7E//5RMla9z7fbSbNTmzVIWWFgY2vehGpCS7FRx69bgdmN55+gHwJRp9Iv64QcSZs8mdfRoXJ9+WrvyPLtdhlnYbETs2EHEtm1Yubla8qdUPacT92rBsqwplmWZirftVRXxK1UPnX322Zx99tnhXoZqgrwZqg4dJKAqLvYvoDriiNKeq3feKR2VvmiR774oh0P2m8rLg40bJZtVXBzyt6MaCBMdLb1TMTFYBw4E1TsFsH/iRPbMnIknNpZmd95JiwkTiPzzz9otzjvowhtQpaVpQKVUPaaBVC0YY6YYY6yKt9atW4d7aUr57frrr+f6668P9zJUE2ZZUt3UsaNM2CsulpI/f65vJ02CuXNlz6qrr5aJf1984fuxLpe8TmamDKTYvz90AylWrNDR5w1KRASeFi0kO1VUFHSwUtCzJ7tffpnMSZOw//UX9g0bQrO+MpMDI7Zvl4AqiAyaUqpuBdMjpZRqREaMGBHuJSgFSEAVFyfBTk5O6T5U1bWzWBYMHgwDBsggikcflYDm1FOrfnzJPqns3Cmb6KakyDHV9BiXi+LDDsO2bx+2ffswUVH+TUApKyKCAyNGkHv66TLOHIh56y0sj4ec884LrBeroshITGQkuN1EpKVhnE7poXI6tYdKqXogmIzU1pCvQikVNunp6aSnp4d7GUod5A2oOnaEkoRBjRmqiAg4/3xYtQpuukmOrVoFV10Ff/3l+/FxcfJ12zbYuhUKCio/TjUBERF4mjWjuE0b8Hik3C+IzM/BTXaNIWrNGhIff5yUsWOJ+vbb2q8xMvLghsAR27dj275dMlRKqbAKOJAyxvSsi4UopcJj5MiRjBw5MtzLUKqSsgFVq1b+BVRRUaXZpbQ0+Owz6Z+67TbJQFUUGSmv4XbLhr67dgXdMqMaOON0UnzYYXgSE6XUr6ZmvapYFnunTyf94YexCgtpcf31NLv1ViJC0T8dGYmJjcUyhoi0NGxpaVj5+bV/XqVUULRHSqkm7uabb+Zm7wg0peohy5JhER06+B9QAVxwAXz9NVxyCSxYAH36wHPP+X6s0yklhd4NfTMzdUPfJslmw5OcTJG31znI7BSWRX7fvuycP5/9EyYQtWYNERkZoVtnSYbK8niI2LZNAyqlwsQyQTYuWpblALoBcfjYW8oYs6p2S2u4evXqZdasWRPuZSilVKPk8ZT2UBUVySa8ERHVn7N5s/RP9ewJl15auiGwr72lPB7ZODgiAlq2lABL21Hqxvr164kMtCfpUPF4sLKyiMjIwDgcB8eeB8PKyTk4PCL++ecpatuW3CFDZNOzUCgsxCosxONyYZKTA950WNUNt9tN586dw70MVUuWZX1vjOnl676gOiAtyxoKvAokUD6IMiXfG6CG/9aUUvXBzpJ6p9TU1DCvRCn/2GySoSo7lKKmgKp9e3j66dLv33wTHn9c9qMaNar8PICyG/pu3y7ZqpYtw7Oh79dfy9fevQ/9azd5NhsmMZEilwvbnj1YOTnyIQsi+PEGURQVEfXdd8S/9BIx77zDvptuwn3UUbVfq8OBcTiwCguxpaXhiY7Gk5QU/EbBSim/BPujkCeAdcBA4PgytxPKfFVKNQBjxoxhzJgx4V6GUgHzBlQdO0Jqqmy2m5Pj3x5RXbtKmeCtt8LAgfDxx5UruLwb+no8ktHasUM39G2SoqLwtGmDp1kzGfBQm6kkdjt7Zs9m77//jX3HDlIuu4yk++7DFqqyP4dDeqjcbuxpadh27gQt+VOqzgQ7k7MtMMkYsyKEa1FKhcHtt98e7iUoVSs2GyQkSNCTnV264a7LVXWG6vjj4YMP4JNPYPp0uOIKyUzNmFH5sVFRcsvNlQ2DmzWDpKSaywlVI2JZkp2Kjq51dgqbjdx//IO8vn2JnzuXmHffJeuKK0K73qgoTFQUVkEB9rQ0PDExeBITNUOlVIgF1SNlWdZXwOfGmLtCv6SGT3uklFIqfDweCaj27JFfVxdQgZTwvf22jFo/7TQJlrZsge7dKz/WGAmoLEvK/eLi6rZ/qrGX9tXrHqmqGIOVnY1tzx5JW9YyODnYP2UMiQ89RH6fPuSfckpoP1gFBVhutwRUSUnykwFV57RHqnEIeY8UMAH4xLKsROAbILfiA4wxC4N8bqXUIbR1q2wN17Zt2zCvRKnQqJih8gZUVSUQ7HYoW906e7b0T51zjpT+dehQel/ZDX137Cjd0Dc6OvTvIy0NfvpJKrMKC6FLF2jTJvSvowJkWZj4eIqdTmzp6diyszH+TDypgrd/ypaVRdQPPxD77rvkn3QS+yZOpKjsh682vBmq/Hzs27ZpQKVUiASbkZoAPEXVPVbGGNNkix40I6Uakn79+gGwYsWKsK5Dqbri8cD+/RL0GCMZquoqsvbvh5kzYc4cyVaNHQs33CAZqIrcbsjLk+EULVqE7ro0LQ1+/13WmpAgwdqPP8IRRzSuYKpBZqTKMgYrJ0fK/SIiaj8tr6iI2LfeIn7OHKy8PHLOP5+sK64oHVYRKvn5WMXFpSV/GlDVCc1INQ7VZaSCDaS2A2uBO4FMX48xxmwO+IkbCQ2kVEOyfPlyAAYNGhTmlShVt4qLpWwvI6P6DJXXrl3SM/X663DGGZKpqkpengRVycly804B9JbmBeqnn+CYY2S9IFMHMzNLj/urvpcENvhAysvtxpaRgS0np1bZKS9bZibxzz2H6+uv2Tl/fugDKS9vQBUbKwFVLUa8q8o0kGoc6qK0zwW8ZIxZG/yylFL1gQZQqqmIiJAhEfHxEqCkp0uGqqqAKiVFBlFceWXpsc2bYckS2eS3bPLB5ZLvs7Jg3z7JTiUkBL/W/Hw53xtIgXyvA9jqqchIPKmpGG92yrIwtZiX70lKYt+dd7L/wAFMTAwUFZF8zz3kjBxJ4d//Hrp1O50YwJaXhy07G09cnAZUSgUg2EDqVWAE8HYI16KUCoMNGzYA0KlTpzCvpAEqKYtEyyIblEADqrJ/NT78EB58EF58UfagGjmyNPlgWfIcHo/0ZWVmSvYoJub/2bvv+Dar6/Hjn/toS16x4x0nhIYALdBQ0i+F0jDTUgjQUkJIm1JmGWG1jDJaEmaBAg20kLIhQMMK88cqe7eMElYhNE0IiZN427Jsy7Kk+/vjWvGIp6xh2ef9eukV/Mh6nhNbaZ+jc+85w+8bEAqZRC4Wz447mjhnzBj9VabxTGdlmb1TdXVm75TH03NI2XDP5/MBYN+4Eecnn1D04ou0HnAATaefTiSBs/+02w1am4QqEDAJVW6uJFRCDCLepX3nAecDnwCvAP5e36K11n00kR0fZGmfyCSyR2oEJJEaEyKRrj1UMPgeqjfegD/+0exZ2n57uOACmD176+8Lh02HP6/X7K8azvYZ2SOV+VRLi6lOwYiqU1vOFwySvWwZ2ffdhwaajz6a5qOPTnyyozUqGIRolGhOjkmoxujvKNlkad/YkIw9UtFBvkWaTUgiJTLEa6+9BsDee++d5kgykCRSY8pwEiqtTXXq6qtNy/Qrr+z/vO3t5pGXZ2ZQDfWetLLSzLoKBk0laqRd+0bj23UsJ1IARCJY9fVYfr+p+oygOhVj27SJ3L/8BfvXX1N9zz3JG2gmCdWISSI1NiR8j5TWOo4JdEKI0UgSKCEMm800iui+5E+pvhMqpeCQQ+DAA02SBPDuu3DTTXD++WYpXkxsoG8gYBK12P6pwe5/y8u7GkvIcr4MZbMRLSxEZ2VhVVWhQiFTnRrBjKhIaSn1V16Jam0Fmw3l95N/xRX4jz+ejunTExd7bJ+X1liBAMrvR+fmEs3JkYRKiE5xJ0RKqR2VUouVUlmdX39TKXWCUmoE22uFEKm2atUqVq1ale4whBg17HaTUG27rdlL1doKLS2mCtWbw2Fan4OZK/Xee2aZ3xlnQOeIti28XvO9tbWwdq1J1uJYFCIykPZ4iFRUEM3ORrW0mHWfIz1n5/Ayx9q1OD/8kKJf/Yq8q6/Gamwc8bl7iCVUXi+W34/t66+x6usT8ncQItPFlUgppX4EfAycQldVawJwA/CxUmqHxIQnhEi2k046iZNOOindYQgx6tjtMHFiV0LV0mKSqv6Sn8MOg7feglNOgaefhh/8AK69tuf3KGWSKZfLJF7r1plzinHAZiM6cSKR8nIIh01ClYBMOvTtb7P5kUcIHHEEvieeoOSII8h68MHEZ+lKmeTN5zMJ1bp1klCJcS/eitQlwHPAZK11I4DW+i1gCrAJuCYx4Qkhku3KK6/kyoE2eAgxznVPqPLyBk6oJkyAiy4yDSnmzoXsbHM8GjWvi7HZup5bv97shwqFkv93Eemn3W5TncrNNclUR8fIz5mTQ9PZZ1N1332EdtwR13vvjWj54IBiFapYQiUVKjGOxbvrcSfgaK11e/eDWutapdSfgZtHHJkQIiX2lM0XQgxJLKHKyzOzourqTELU15aXsjL405+6vn78cVi8GM46C37xC1ORAtNwzek0A33XrjVLCidMSEhPAjGaWRbRggK0z4dVU2OybK93xMlPeNttqb3xRtMkArCtX0/uzTfTtHAhkUmTEhF5l+57qJqaUI2NRCdMQGdnyxtYjBvxVqQagcn9PFcU5zmFEGnw6aef8umnn6Y7DCEyRvcKVW6uaSIx2CqtadNMq/Q//AH23htWrDBVqhiPxyz5a2oyCVVTk3l+zz1HV6OJt982D5EY2u0mUl6Ozssz1alElCW7DQN2rF6N+513KDnqKHJuuslcI9FiS/48HmyNjdi+/hrV2GjaYAoxxsWbSN0FXKqU+rVSqgBAKeVRSs3FLPuTQb1CZIjTTjuN0047Ld1hCJFxHA7TgS+WUA205G+XXeChh+D++01XwDPOgN5bE2MDfT0e2LwZvvpq8ARNjAGWRTQ/n3B5uXkTBAI9s+wRCO67L5sffpjW2bPJWbaMkiOPxPPccwk591YsqyuhamiQhEqMC/HWXhcBhcBSYKlSqgNwAAqTRP02MeEJIZLtT93XHwkhhs3hMAN3J0yAhgbz6GvJn1JmltOsWWY+VGyPVGsrfP457Lab+dqyzHMdHWb/lM9nErbhDPQVGcjtJlJWhvL7sdXXo+32rjWgIxAtLKRh0SJaDj+cvD//Gednn9F24IEJCLgfsYQqGjV/j4aGriV/yZp5JUSaxDtHKgqcrJS6CtgPKAaqgXe11p8kMD4hRJJ997vfTXcImaekBKqqur6O3S0XF5tSghiXhppQWZbp8Bdzzz1w+eXwox+ZGVSxUUAOh3kEg6Y6NWGC2UOVqSN8RuNA4FHHstB5eYQ9HqyaGlQgYMqU/U2GHobQzjtTffvtW5YPOj/8EN9TT9F06qlEJ04c8fm3Yllon08SKjGmjWg3oNb6K+DOxISSeZRSizHVuR5KS0tTH4wQcVq5ciUAM2bMSHMkGaR7EjWU42JciSVUeXkmmWps7L8pBcDRR5uhvkuXwv77m25/Z59tBvKCqUS5XNDcbPZOxRpeJODeWoxWLhfR8nJUU5PpiJeg6hSWtaW06fjf//A+/zyeV17Bf9xxBObNM51PEq13QtXYSDQvD52TI29ikfGUlsXXCTdz5kz9/vvvpzsMIYZkn86PiV+Vj4m7xD46789rr/X/3N57931cfr7jVijUlVDZ7eY+tq+Eqr4e/vIXuPtu02Di/vu3/p5o1CwFtNlMspaVNfRGb4moCMUaTcTTAKOv669evRpHppbYUiUUwqqtRbW1Jaw6FWP/+mtyb7wRzxtv0DFpEk1nn00w2d1NolHTprJzb5jOyhqzCVVHRwfTpk1LdxhihJRSH2itZ/b13Nh85wohhmzJkiUsWbIk3WEIMWY5nWbV59SpZr9Tc7O5j+z9OWZ+PixaBG++CZdcYo5t2gQ33tg1tNeyTPLkcJjZU+vWmXMNZvly+OwzeP112Gkn87XIEE4n0dJSooWF5pfd2do8EcKTJ1N37bXU3HAD2O3Y161L2Ln7ZVnmH4LTiVVba5pS+P0Ja7AhRCpJo38hxjlZ0teHwT6yH6gEIJUn0Q+n02yvy8831aemJlOh6uxUvUVsSR/A88/D1VfDXXeZGVQ//7lJoux20/0vFDLJVE6OWfLX18qs5cvNkOCHHoK99jKJ2vHHm+fmzx96/JWV8PHH5j4+FILttusZq0gipdA5OUQ8Hqy6OqxAwLQ4T9Beo/bvfY+qbiVQ77PP4vjiC/wnnGD2NCWDzWYSqkgEq7YW6uvHfIVKjD3yThVinHvvvfd477330h2GEONGLKGaOtWs1IpVqPpyzDFmmO/UqXDhhWZ53JNP9jxXTk7XQN/a2p7dpvfZx7Rav+MO2Hdfk4Ttu6/5+owzhh5zZSV88YVp4z57tqlqffGFOS5SyOEgWlxMpKgI2tu3DN5NCLt9yyBd+7p1ZD34ICVHHIHvsceS28I8llA5nVg1NaZC1dwsFSqREeLaI6WUWgv090KNGdj7H+BmrfU78YeXmWSPlMgkskcqDr279sVI1z4Rh/Z2U6Hy+/uuUIFZBvjyy/DHP8LOO8Of/9z397S2moJpYaFJsPbd1yzna2/v2e2vo8P0LnjzzaHF+PHHJony+83XU6aYfV+x4wOJbbmRPVIJFg5j1daa6pTXm/BOeI5Vq8i7/npcK1cSmj6dxvPOI7Tzzgm9Rp8iEfPJgM1GtKDANKrI0AqV7JEaG5KxR+p6YCImYboD+CNwNxAE8oFnAS/wmlLqx3FeQwiRAn/961/561//mu4wMsvmzeaude+9zUNr85AkSsTB5YLSUthmG5NE+f1bV6iUMh39/vEPuOwyc+zjj+HII+HDD7u+x+cz59u82VSonn4avvnNrROmN980x4cqGDRDh7vLzU3odh0xXHY70ZISIiUlEAqZZhQJ1LH99tT87W/UXX45Vqx7YCrYbFs2Alo1NdjWrzcVKmmOJkahePdIzQBe0Vof1v2gUuoKTBJVrrX+mVLqJuDizmNCiFFop512SncIQghMAlRWBgUFUFdnEiqHo2eFKtZsAkyy9MUXMGcOHHQQ/O53MG2auQ+NDfTdsAFOPtnsibrjjp57pK64Yujd90IhU3CNFQZ23NEsI5wxI74OfiJxdFYWEbfb7J1qbjZ7p+wJ2gKvFG2zZ9M2a9aWDXjZy5ZBKETzggXJnRIdW2oYDmNVV5s9VLEK1VBbVQqRZPFWpA4H7u19UGsdAW4FftJ56B+A3KUJMYq9/fbbvB3raSyESLtYQjV1qrlPbW7uu/Lzwx+aduS//a3pyL/ffrB4cdfzDodJqA46CBYuNPOpXC44/XSTRA2n0cR228FHH5nlfNGoSaI++sgcF6OA3W72TpWVQUdHwqtTuFwmedEa+1dfkXvbbZTMm4fnpZeSXymy282nB3Y7VlWV2UMVCEiFSowK8SZSYWBSP8+VYJb4AXiAjjivIYRIgQsvvJALL7ww3WEIIXpxuUxXvClTugby9k6osrLM8N633zaNKfLzzXGtu/Yzud1w1FHwjW/Ad78Lb7wB8+YNL5bycthhB7Oc8IUX4NNPzdfStW900V4vkYoKoj6fSTbC4cReQCkaLr6Y6qVLiWZnU3DhhRSeeir2tWsTe52+9E6o1q9HtbRIQiXSKt7a773AYqVUHfCk1rpJKeUGDgQWAY8rpfKBM4B3ExOqECIZbrnllnSHIIQYgNttEpZgsGvJn9PZc1XVxIlw6aVdXz/3nEmwFi6E444zywNtNnPPWVNjmlsMd6BveXlXYwlZzjeK2WxECwvRWVlYVVWoUMgs90vgcrjQd75D9T334HviCXJuvRUVCiXs3IOKJVThMLZNm9Aul1nyl+C/oxBDEW8i9TtMs4l7AK2UCgFOQAEfAecCuwM7AAckIE4hRJJsv/326Q5BCDEE3ROq2tq+E6qYbbeFmTPhyivhzjvN8j+tzX1m5z0oGzea1xYV9d0pUGQ27fEQqajAamjAampCu92J2zsFYLPRcvjhtBx8sCmZArnXXUekvJzAEUck9lp9sdvNjKvuCVV+viRUIqWGvbRPKWUHdgROwyRKJwOXAedgkqbvaK0bgbeAMq31vxMXrhAi0V577TVee+21dIchhBgitxsmTTJd/hwOs+Svvb3n92y/PSxbBo8+apKv886Dzz/vet5uN/unolEz0HfTJtNQIlmWL4fPPjOt2HfayXwtUsBmIzpxIpHycgiHk7MUrjOJIhzGsW4deX/+M8ULFuD6178Se53+2O1miK9S2DZtwlZZiWptlSV/IiWG/XGB1jqslHoTOFZrvQL4sp/vC4w0OCFE8i1atAiQOVJCZBq3GyoqTKv02lqTUDmdXfe1ALvvDk88YdqmX321ORYMwr//bZbnuVzm0dpqXp+fDxMmJHYk0fLlcNFF8NBDPbsGgtmzJZJPu92mOtXYiNXQgHa5eg4WSwS7ndobbsD9+uvk3XADhWecQdusWTSecw6R4uLEXquf6+usLOjowLZxI9rtNkv+3G6pUImkibfZxN3AiZ3VKSFEBrvzzju588470x2GECJOHo9JqCZPNpUmv79nhUop+NGPuhpRPPig6eD385/DJ590ncPnM1351q4150jUB/pnnGFar++7r7l333df8/UVVyTm/GKILItofj7h8nLzy01GdUopgnvvzeYHHqDp1FNxfvYZOsGDggflcJglf4Bt40asysrEdzEUolO8iVAL8H/Al0qpl4HWXs9rrfWZI4pMCJES2267bbpDEEIkQCyhamszDSX8/q6KU3fz5pmq1I03woEHwqGHmqV/U6eaZCoSMUv96urM7Civd2Rx1dWZSlR3e+3Vc6mhSCG3m0h5Ocrvx1ZXh3Y6t8yIShink+Zf/Yrm+fPNubUm//e/J/j979N64IFdA8mSyeFAOxyojg6sjRuJut3oWIVKiASJ9518FNCEaS6xP3BIHw8hRAZ48cUXefHFF9MdhhAiQTweU52aPNks0Wtu7rn/ye2Gk06Cd96BM8807cx/85uu52MDfZWCr782Q31bWsxze+45/I593/ymWc7X3ZtvmqG+Ik0sC52XR3jSJPOLDgTMhrlE60zQlN+PfdMm8i+5hMITTsDx2WeJv1Z/HA50VhYqGsW2YYOpUPU1mE2IOMRVkdJaT010IEKI9Lj88ssBOOAAabApxFji9ZpkqrXVVKjC4Z6FgJwcU4k65hhobDTH6upMl7+TTjLPO51mud+6daYTYEHB8IsXF11k9kTdcUfPPVKytG8UcLlMdaqpCVt9Pdpu37qEmQA6N5fq22/H++yz5N50E8XHHUfLwQfTeNZZ6JychF+vT04n2ulEhUJYlZVEPR50fr5UqMSIjHiPk1LKCeQDdVprGb4rRIa599570x2CECKJYgmV220qU36/+e9YQlRUZB4AL78MS5bA3XfD6aebJMtmM8WK1lZoaoK8PLPfaqgJ1fz55s8jjzSJ2je/aZKo+fNh9epE/23FsCllqlNeL1ZNjRnk6/UmfvmdZdF68MG07bMP2XffjefllxO/pHAo+kioovn5fc8REGIQcf8rUUodopR6D7M/agPwXaXUY0qp4xMWnRAi6SoqKqioqEh3GEKIJFLK3Be73SapApNQ9W55PneuGeY7YwZcdpmpID36qHm9x2OW/LW0wJo1UFU19Jbp8+fDt74Fs2bBp592JVdiFHE6iZaVEZ040WTNSVr+pn0+/AsXUvXAA6Ya1N7OxDPOwP3GG6ltWe50miV/4TD2ykqsTZuS9ncWY1dciZRSag7wOFANnN3tPP8Gliqljk5MeEKIZHvuued47rnn0h1GZnr1VfMQIoN4vTBlimlMAWYPVUe39SQ77wz332/alZeUwBtvdD3XPaEKBEyHv+EkVGKUUwqdm0ukosI0amhuNt1HkqGz/bqtpgZbVRUTzzmHiWedhX3t2uRcrz8ul0moOjpMQrV5syRUYsiUjiP776xEfaS1PkEp5QOagb201m8rpf4CfF9r/Z0Ex5oxZs6cqd9///10hyHEkOyzzz6AzJESYqzr/KfeI/fXumsPVXu7qVh1Hy+kNWzebDoBhsNmT9Xvfgd77NH1fFubudeOLfnrbzxRX9dfvXo1jkTPMxIuAGoQAAAgAElEQVSJoTUqEDDL/Wy25O4lCofJeuQRcm67DdXWRmDuXJpOPTUp+7UG1d6O6ugg6vMRnTBhRDF0dHQwbdq0BAYn0kEp9YHWemZfz8W7tG8n4Kl+nnsF2D7O8wohUuyBBx7ggQceSHcYQog0UMq0PJ8yBSZNMnuhuleoYs+DSajWr4cjjoBf/hI++8w87/VCVpZ53Zo1UF3ds8IlMpRS6OxsIhUVRF2u5Fan7HYCRx3F5kceoeWQQ3B+9lniBwYPVaxC1d6OfcMGU6HqPphNiG7iTaQqMclUX0oxrdGFEBmgpKSEkpKSdIchhEijWMK0zTYmoYpEzB6q7gnRrFmm497vfw///rcZ8nvuuaYq1T2h8vsloRpTHA6ixcVESkpMtSaJw22jEybQeMEF1CxdCpaF1dDAxFNOwfnhh0m7Zr/c7q6EqrISq6pKEiqxlXi79t0CLFZKVQIvdB6zlFJ7AhcAyxIRnBAi+Z56yhSXDzlExr8JMZYNZfVuLKGaOtU0laiuNkv3YjweOOUU0yxi6VJTNFDKJFONjTBhgkmotDYJVexYXl7S/loiFZRCZ2URcbux6uqwmpvRHg/YR9z8uW+x/VMbN2KvrKTo5JNpPeAAmk4/3SR0qeR2owEVDGKvrOxa8peOjoNi1IlrjxSAUupPwG8wQ3kVEMFUuB4BfqG1DicqyEwje6REJpE9UkKI/mhtWp5v2GDubX2+vrtiv/oqnHiimT910kmmGUXs9a2tZsngcceZ++7XXut6neyRykwqEMCqrUWBSaiSea1gkOx77yX73nvRQPPRR9N8/PEmg08DFQxCOEw0O5toXt6ACZXskRobBtojFfdHCVrrczsbS8wGigE/8IbW+qN4zymESL1HHnkk3SEIIUYppUw1KSfHVJhqaro693W/j50yBfbbD/78Z7jnHjjzTLOPyuUyyZfWpllFRwfU1ppzJquYIZJvS3Wqvh7L709qdUq73fhPPJGWOXPI/ctfcKxdm7YkKhYPWmO1tWEFAiahys2VCtU4FXdFSvRPKlJCCCHGonAYGhqgvt7cN/cuRqxcCX/8o9lL9Z3vwJNPdt3zHnGESajuvddUqPLzobZ2NR6PVKQymWptxaquRmltEqpkJznhMNjt2NesIe+662g680w6pk9P7jX7ozWqvR0ika4KVbcKq1SkxoaEdO1TSl3Q2ep8OBfOUUpdPJzXCCFS69FHH+XRRx9NdxhCiAxgt0NhoWlK4XZvPdR3xgx48EFYvtzspVLKVKFefbVnUwqfz+yfWr/eQX29RXjcbgbIfNrrNZ39srJQLS0k/ZfZWfmyb9qEY/Vqin71K/KuugqroSG51+2LUmi3G+31YrW0YF+/Hqu2VrqsjCPD6dpXCKxRSl2plBow9VdK7aKUug74cpjXEEKk2I033siNN96Y7jCEEBnE5YLycpg82XzduzP2rFlw0EHmvx9/HH7xC9Mu3e83x2IJlccTpanJ4uuvbTQ0KEmoMpXNRrSwkEhZGYTDJqFK8oqn4Pe/z+aHHyYwdy6+J5+kZO5cfCtWJPWa/VIK7fGYhCoQwCYJ1bgxrKV9SqmdgD8CBwPVwL+BWiAA5GCSrd2APOBZ4Fyt9RcJjnnUk6V9IpM0NZlpBbm5uWmORAiRibQ2iVR1tVmy5/P1XN3V0WEqVBdfbP579mw4/3zYYQdYu3YtDoeDaBSCQYXWmvz8KNnZGpstfX8nMQKRCFZjI1ZjI9rlSsk8KPuaNeQtWUJo++3xL1yY9OsNSmtUMEi4vZ0pu+5qNgVKU5WMNdDSvrj2SCmltsckUz/AzI3KxyRTa4B3gEe11mvjjjjDSSIlhBBivIlEzHK92tq+90/99KdmqG9jI1RUwPPPw1dfre3RtS+WUIFmwgRJqDKZCgaxqqvNUj+vN/l7p7Q2b0K7Hdc//0nWI4/QeNZZRCZNSu51B9ARCjG1pMTEVVAgXVYyVMK79mmtVwGrgOtHEpgQIv0efPBBAObNm5fmSIQQmcxmM/eKOTkmmfL7TSMzl6vr+fJyePppU71SCpqbLZYvz+EXvwiQlxfFssDr1USj0NBgo6FBEqpMpd1uIpMmmepUQwPa6UxuZzultiQpttpaXB98QMlRR9E8fz7NxxyD9g1rm3/iYooNVmtsNF1a8vMloRpDZP+SEOPc0qVLWbp0abrDEEKMEQ4HlJaalug2m1n2133vU36+WdYH8MEHbu67L5vDDy/hzjuzaW01VYtYQuV2Q329ja+/ttHYqHrswxIZwLKI5ucTLi83SUUK9k4BtM6Zw+aHHqJ19mxyli2j5Mgj8bz0UtKv26/eXVbWrDGfNsimwIwniZQQ49wzzzzDM888k+4whBBjjNttlvB19h8gEtn6HnqffVq5//4qZs5s55ZbcvnZz0p4+GHflu+zLPD5JKHKeG43kbIyovn5qNbWnq0ekyRaWEjDokVU33EHkeJirFink3SKJVRer5kjsGaNqVLJGzpjSSIlxDjn9Xrxer3pDkMIMQYpBdnZpl26w2H2QLW29kyott02zDXX1HH77dVss02Yd95xb7WdpndCtX69JFQZx7LQeXmEy8vRSkEgYN4QSRbaaSeqb7+dlsMOA8D32GNMuPRS01UvXcwb2iRUdXWSUGUwSaSEGOfuu+8+7rvvvnSHIYQYwyzLDOl96y2zhyoQgGCw5/fsvHOIm2+u4fLL6wHYsMHGcccV8c47rq0qVC4X1NdbWxKqFNyPi0RxuYiWlxMtKIC2NmhvT/41Lcs8AMvvx/uPf1Aydy7Zy5alpDo2YFw+n+nMEkuoGhokocogcSVSSqnDBnn+zPjCEUKk2u23387tt9+e7jCEEOOA3Q5FRaZC5XRCINBzdpRZ+WSyptpaGw0NFmedVcipp07kk0+6GhXE7j9dLnj9dTfPPOPB75eEKmMohc7LIzJpEtpuT1l1CqD5V79i8/LltO+2G7k33UTJ/Pm40t1puXtCVVMjCVUGibci9ZhS6kGlVFH3g0qpnZRS/0S6+QmRMV544QVeeOGFdIchhBhHXC6YNAnKyiJEo6YHQe97xhkzQjz00GbOOaeBtWsdnHBCERddlN9jWaBlmXtPhwNqa81g32QnVB984OKDD1zJu8B44nQSLS0lWlhoqlO9y5RJEqmooO7aa6m54QaTyI2WhMWyICurK6Fau1YSqlEu3kTqRGBf4HOl1K+UUk6l1OXA+0AFMD9RAQohksvhcPSY4yKEEKni8WgmTYowcWKUUAhaW1WPRMnhgLlzW3j00c38+tdNlJWFt+yfamzsuoWJfaDvdKYuoRIJohQ6J4dIRQXa6UQFAilLHNq/9z2q/v532nffHYCcv/2N3OuvRzU3p+T6/YolVG636e63dq3p9idv6FEnrkRKa30HsD2wArgDqAbOA24BdtBaP5SwCBNMKZWvlFrd69iZSql1SqlVSqmHlVK56YpPiFS7++67ufvuu9MdhhBinFIKcnI0FRURcnOjtLaqzqG8XbxezfHHN7Nwoem89uGHTg49tJS//CWXlpau77XZxndCldHVMoeDaEkJkaIiaG9Hpag61X1AmWppIeuhhyg54gh8jz2W/kpQ9zWs1dVmyV9TkyRUo8hImk2Egdi7PAdoAt7VWqc5je+fUupK4C3A3e3YDsD5wB5a6+2BSuCi9EQoROpJIiWEGA1sNsjPjzJpUhi3O0ogoPrtA1BaGmH//Vu5//4sLriggGee8fZIvronVDU1JqFqbh4/CVXGUgqdnU1k0iSiLpepDKUwmWk6+2yqly0jPHUqE666iqJjjsGxalXKrt8vm81UqFwuqKoyFSpJqEaFeJtN/AJYBZwE/AmYBrwO3KuUelUptdNIA1NKPaSU+msfx09XSv1PKdWqlHpfKfWjYZz2jc54u/sG8Het9cbOr18Cdo4vaiEyz6uvvsqrr76a7jCEEAIwyU9xcZTy8ghKaZqb1VZzS0tKIixa1MB991UxfXoHjz2WxbHHFm11Xxm7/3Q6obpaEqqMEatOlZRAKIRqa0vZpTumT6dm6VLqrrgCq6UF7XQO/qJU6f6GjiVUfr8kVGlkj/N19wIfAD/WWn/UeexnSqmfAH8F/g3E/c5TSv0AOASzbLD78VOAJcDVwLvAT4CnlVLf11r/a7Dzaq2fVUptA1za7djTwNOd588DzgWeijd2IYQQQoyc260pL9e0tGhqay3a280efKvbR8DTpoU57bQm/vtfB7m5USzL3FO+/babPfcMbvne2P1nOGwSKrvdVL98Pt3jfGJ00VlZRNxurLo6rOZmtMdjWj8mm1K0HXAAbfvuu2Xp34RLLyVcVkZgwQK02z3ICZIs9oaORGDTJvMzKSw0x+QNnVLx/rTPAXbvlkQBoLV+HPgmcFs8J1VK/VIp9V9Mdcvd6zkHsAi4SWt9odb6ca31McBK4A/xXK/X+X8B/AeTIErXQTFu3Hbbbdx2W1z/ZIUQIqmUgqwszeTJEfLzo7S1QVtbz4YUANtt18EBB5iqxWuvuTn77Ikce2wR777bc7+Q3W7uNe12s+Rv/XpToep9PjGK2O1Ei4uJlJVBR0dKq1Nb9k+Fw6j2dnJvu43iefPwvPgio+JNY7OZidcOh0movvrKVKhGQ2zjRLzNJq7XWvdZR9Ra+7XWC+OMZy2mCnUB0NjruV2BYuD+XsefAfZXStmVUucppWr7eHx7oIsqpe7CNMs4TGv9G6219JkU48aDDz7Igw8+mO4whBCiX5YFeXkmofL5TEOK9nZTXVq71sl//+tm5Uo31dUWs2YFufjiehobLU4/vZDTT5/I55/37Exqt5s9VHZ715K/QEASqtFMe71EKiqI+nyms1/v9Z7JZLdTf8UVVP/tb0Szsym46CIKTzkFW2Vl6mIYiN1uEiq73SRUa9dCc7MkVCkQV31UKfXyYN+jtd5vuOfVWr8JvNl5jZN7Pb1L55+f9jr+GaZ6VaG1vga4ZjjXVErNAWYC/6e1HvBjDqXUkN6Ru+2223BCECKtXnzxxXSHIIQQQ2JWMEXJyYnyxRcO1q/3sMsuitzcKHl5is8+8wBtHHxwK7Nnt7JiRRZ3353NxRfn8+CDVVuteopVqMJhqKoyS/4KCsySP6X6DEGkk81GtLAQnZWFVVWFCoXMcr8U/bJCu+5K9T334HviCbL+/ne012ue0DplMQwollCFw7Bxo6lUxZb8jYb4xqB4F5q2At2TCguYgEl2NPDwCOPqSz4Q0lq39Dre0PlnTpzn/R5QBnyout5kb2utj4vzfEIIIYQYxNtvw4YNbuxx7nlZtcrFN7+pqKuziI39yclRvPOOl6lTTbu/b30rxKWX1lNXZ/Hpp06CQcWTT/r44Q9bycvrubBml11CWxIqh6NrD1X3+8/Nm22sWuWio8MiHLaYMiVESYksYkk17fEQqajAamjAamoye5ZSsXcKwGaj5fDDaTnsMLO0Lhpl4plnEtxzTwJz56YujoH0lVAVFZkyrCRUCRXXb1trPaev40qpbOBOYN1IgurHYBNDh7RoVmv9FTCp29e/B34/xNcO6d03c+ZMqaWKjHHzzTcDcOqpp6Y5EiGEGLqODouCgijt7ZqODnPP6HKZ493Fhv4CrF7t4OWXPbz2mof992/lwANb8Xq7/i+7e4Vq82YbTqemoCCK16upqrKxbp2Hb39bkZMTpaBA8cknpgImyVQa2GxEJ07sqk61t5sKUaoShc79U7HhvXlLluB7/HEaf/Mb2r/3vdTEMJhYQtXRARs2xFpiQip/TmNcQtNmrXWzUup84FXgkkSeG1N5ciql3Frr7lPaYsNzaxN8PSHGhaeeMk0qJZESQqTSnnvC6tVBHI7BPiftWzhsUVCgtizX23bbDtavt6ioUEyfHqKvxmq77BJi773buOWWXJ591sdbb3k4+mg/P/95gI8/3rrZcCQCoZDCbtdUVzv47ncVgYBFIGD2bWVlKV5+2cf227f3G+duu/X/nBg57Xab6lRjI1ZDA9rlMhWYVF0/N5faG2/E/cYb5N1wA4VnnknbD35Aw0UXmax8NHA4zKN7QlVUJAlVAiSjR2IeZhleon3R+WfvGVU7ABu01vVJuKYQY96zzz7Ls88+m+4whBBiWKZMCfHJJxaNjWYuVEODYu1axYwZbbhcmkBA0dGx9evKyyNcemk9995bxbe+FeKdd9xbkrHee/NtNlPRsixoarITjfbscZCTo7eqgMXjgw9cfPCBa/BvHEBs2eGnn3r41788bN5sG3FcGcOyiObnEy4vN7/ElpbUNlpQiuCsWWxevpymhQuxVVcT9fnMc6Op4YPDYSpUlmUSqnXrUv+zGmPibTbx236eKgLmA/+MO6L+vY3p5HcU8H5nHHbgSDrnQAkhhBBifDDL6dp4+WUfHR0WO+ygmTKla5ldMGjmTwUCZv6UrVdeMX16B0uW1NLWplAKiosjnHtuAcce62effYJbfVCvlEV+vkJrU6GqqAjT3Aw77KDTXnXavFmWHQLgdhMpL0f5/djq6sww3VQO1HU6aT76aJoXLADLQvn9cMghcPTRcMQRo2fGU6xCFQrB+vXgdpsKVQobd4wV8S7tu7af4y2YhKd3x70R01oHlVJXAFcppZqBT4DjgKnAYYm+nhDjxQ033ADAmWeemeZIhBBieEpKIluW1fVOZsxA3wiBgGlIEY2qLdWl7jwe82l8fb1Fe7vi/PMn8q1vhTj11CZmzuw6Z1lZiM8/95CTY1ZE/e9/Nr78UjF9emvKmrb1V7VatcrFt789tGWH6U76ks6y0Hl5hD0ebDU1EAiYX1gqk5jOa9kaOyf5/OY3sGwZXHopfOc7qYtjMLFEMxSCr7+WhCoO8c6Rsvp5ZGutf6S1XpvoQDtdh5kxdTxwH6ZT4P5a6zVJup4QY95LL73ESy+9lO4whBAi4ZSC7GxNRUWE/PwIra1moG9fdtyxg/vvr+Kii+qprrZYuLCQ3/ymgEhnQaeoKEpxcRsrV2pef93iyy+jFBW1EYlAZaWtz0HBqdLRYZGT0/PiiVp2mLFcLiLl5UQLCswQ3/bUJ5Dh0lJ48klYsgQqK0116qyzIJVDhYfC6YSczubX69ebpKq1Nb0xZYhR0KOxb1rrbfo4poE/dT7STim1GFjU+3hpaWnqgxEiTk8++WS6QxBCiKSy2cxAX58vQkODRXOzhdOpt1r1ZbfDoYe28sMftvLII1nU1tq2LAmsr7coKopuaa2+yy6hLa/r6ICNG224XFEKCjRud3LmUPVXTerdeGPatA7q6tSoWHaYVkptqU5ZtbVmkG86qlNz58KPfww33giffsqWTiijZf5UTKxC1d5ukimv18yh8njSHdmoNeRESil14zDOq7XWY36dkNZ6MbC493Fpfy6EEEKMPmacjhnoa+ZPWXg8eqvRP243LFgQ2PL15587OPHEIg47rIU99mjbqvpjtpxoOjoUGzdauN1R8vP1lmWDyWYab3jIylLk5Gjq6hSffGIxZcooq3yki8tFtKwM5fdj1daajLmvto7JlJUFF14I0ahJnjZtMnunzj0XZs8eXQmVy2Ue7e2mIYUkVP0aTkXqkGF8rwbGfCIlxFhw7bVmy+M555yT5kiEECI13G4oK4vS1qapqbFobzf3iP0VKoqKIhxySAuPPebjqae8HHBAG9tu20FWVv8JVWWlhcfTVaFKpsEabwhMdSo3l0hndcpqbjZzp3p3IUm22Jusrs6UMo89FvbeGxYvhunTUxvLYLonVLEKVXn56GmaMQoMJ5E6GXhfa12XrGCEEKn3zjvvpDsEIYRIOaXA6zX7p5qbFfX1FkqpPpflFRRE+d3vGpk/P8DVV+fx9NM+PvzQxcMPb96qmgVdCVUopNiwITUJ1UCNN0Q3TifR0lJ0IIBVU4Oy2dCprk4B7LQTvPAC3HMPXHcdHHCASaoWLx5d1SnoSqgCAdP/P5WdEEe54SRSzwAHAi8AKKUcwB+AW7XWG5IQmxAiBVasWJHuEIQQIm0sC3Jzzf6ppiaLxkYLu133ufJr8uQwJ53k58ADW3E6zZJAreH1193stVdwq+KG2XLSlVB5vWbJX7IrVGIQSqGzs4m43Vh1demrTjkccMIJ8NOfwjXXmEQllkSNtv1Tok/Dqc2pzkeME7gImJzQiIQQQgghUsxuN5WnSZPCuFya5ua+B/oCTJkS5kc/MvuP/vUvF+edN5EFC4p5/XV3n537nE7TPTASUVRW2ti40SIYlJvktHM4iBYXEykuhvZ2090vHQoK4OqroXOpPR9/bJpT/DMZY1lFIo10kaP8r4AQGe6qq67iqquuSncYQggRl912a0/oUjaXC0pKopSXm/1FgQBbWqD3Zffd27nyyjrCYTj33In8+teFrFzZ99InpxOysjThsEmoNm2yCAYTFrqIR6w6VVFB1ONBNTeb5WtpigUAv9/sofrZz+Dkk03rdDEqyW4xIca5lStXsnLlynSHIYQQo4rHYwb6FhVFCYWgtbXvOVFKwf77t7F8eRXnn9/Axo12Fi3KH/Be3OUyCVVHh2LlSifvvedi5UoP//qXh82bU7y8TBh2u6lOlZRAR0f6qlMAe+0Fr78OZ59t9lHNmgU33ZS+eES/Ru0cqUwgc6TEWPDAAw+kOwQhhBiVYgN9vd7Y/ilFKKRwOrfOqOx2+OlPW/jxj1vZsMGG3W6and10Uy5HHRWgrGzrslZTk0Ug4GbnnS18Pigutvjf/zyAdNxLF52VZfZO1ddj+f1oj4c+O4okm8cDv/0tzJsHl10Goc65ZbFsXvZPjQrDfWf0tTty3O6YlDlSQgghxNhns0F+fpTsbPjf/6K0tVmEQn03L3O7NdOmmXLUf/7j5NFHs1ixIou9927j4INbyM7uukVYu9bJjBmK1laLxkaw2+1obfHEExY779yOy9Xz3NKNL0XsdqJFReisLKzqalRHh+nsl47kpbwc/va3rgTquefg9tvhkktM5z+RVsNd2rdMKbVGKbUG+LTz2COxY90e/0twnEKIJLnsssu47LLL0h2GEEKMeg4H5OdrJk6MoJQmEFADLuHbddcQl19exx57BHn5ZQ8XXljAU095t7wmHLZ6JFYOh2l4EQrZqK21UVenthQiROppr9fsnfL5UC0t6ds7BV1JXEcHfPklHHggnHee2Usl0mY4Fal7khaFECJtVq1ale4QhBAiY8SqQlpDS4umrm7ggb777BNkn32CfPWVnaVLc1izxsGuu4ZQCiIRi7w8teV122wTpr5eMX16mBkz2gkGIRxW+HxRJkyIpvBvKbaw2YgWFnZVp0Ihs9wvXUvrDj3U7Jm6/nq4+2546ikze2revPTEM84NOZHSWh+bzECEEOlx3333pTsEIYTIOEqZhhFebwS/3wz0tay+B/qCSZKuvrqeYNC8tr7e4qqrcjj88BD/938RcnM19fWKzz6zKCszjQ7MLCtNMKjYsMFOVlaUvLzoVkv+RPJpj4fIpElYjY1YjY1ol8uUENMhLw8uvRQWLDBJVGz+lcyeSjlpNiGEEEIIESfLgrw8TVZWhMZGC7/fDPTtL9mJDfptaLBwOuH667P4xjfC7LFHiL32aqW8vI2iomgfrzEJVWWlHZ9PEqq0sNmIFhSgfT6s6mpoaQGvN33Jy/TpcP/9XV/fequZPXXxxTB1anpiGmek/bkQ49zFF1/MxRdfnO4whBAio9ntMHFilPLyMA6HGeg70Jaab3wjzD33VHPZZXX4/XDffV5uvTWLvLz+l/C53eDzaYJBi8pKO1VVluyhSgPtdhOZNAmdl4dqbaXfyc2poFRXIud0wltvwX77wZVXmiFoIqmkIiXEOLd+/fp0hyCEEGOGywWlpVHa2jQ1NWb/lNvdtfqqO8uCH/6wjYkTI7zxhoeODrWlE2BdnUVBQd9JldttGlQEgxbr11tkZ5sKlXT1SyHLIpqfT9TrxVZbm/7qFMCxx8JBB8Ef/2jmTj38MFxzDcyenb6YxjhJpIQY5+666650hyCEEGOOx6OpqIjQ3Gz2T2mt8Hj63j9lt8O++7axyy6mvLRmjZ1f/rKYAw9s5de/9lNc3PdMKbdbozW0tVk0N1vk5ETJzY322ZZdJInbTaSsDOX3Y6uvRzscfffFT5XiYliyBH71K7PELyfHHB/p/qkjjoBIBN54IzFxjhGytG8ElFKLlVK692Pjxo3pDk0IIYQQaaYU5OSYhCo3N0prqyIYHPxmtqAgypFHBnj+eS9HHFHCDTfk0tTU9y2bUiah8vk0ra0WGzbYqamRJX8pZVnovDzC5eVopcySumiauyzuuis8+STsvrv5etEiOOssqKpKb1xjTFyJlFLqr0qpPRIdTKbRWi/WWqvej7KysrTF9Pbb5iHEUF1wwQVccMEF6Q5DCCHGrNhA30mTwrjdUQIBNeC2mtzcKGee2cQjj2xm9uxWHnggi/nzi2kfYOVeLKHyejUtLV0JVTq374w7LhfR8nKiBQXQ1oYa6BeWCt0rUFlZ8MQT8IMfwF//yoBvJjFk8Vak5gNvdg7fvUIpJaOVhchQdXV11MlAPyGESDqnE4qLo5SXm6V6gzWkKCmJcPHFDdx/fxWnndaEy2VWaL38sqff1ylllhXGEqr16+3U1kpClTJKofPyTDMKux2am9NfnQIzvPfll2Gvvcweqv32g/ffT3dUGS/eRKoI+DHwInAC8JFS6mOl1LlKqfSVY4QQw3brrbdy6623pjsMIYQYN9xuTXl5hJKSCB0d0NY28L32ttuGOeigVgA++cTJBRcUMG9eCf/4h6ff13VPqAIBi/XrbZJQpZLTSaS0FEpKoLXV/JLTbepUuPNOWL7cVKgmTjTHR0Oil6HiSqS01hGt9T+01r8GSoHZwBvAb4B1SqmnlFJzlJKpYEIIIYQQvcUG+k6eHCEnJ0p7O7S1KbQe+HU77xziuutqcbuj/OEPBRxzTBH//G6LW6wAACAASURBVKer39d1JVQQCFh8/bWNujpJqFJCKcjNNQmM02mqU5G+G4ek1KxZ8NxzsM025usTTzSNKRob0xpWJhpxswmtdRTY3PmoBWzAt4HHgFVKqT1Heg0hRPKcc845nHPOOekOQwghxiXLguxss+TP54vS0qIG3Qu1115Bli2rZtGievx+i0suyR90y0ssofL5wO/vSqgGWlooEsThgEmToLTUVKZGQ3UqVuvo6DCVqTvvNPun7r13dCR7GSLuREoptatS6nKl1OfAJ8BC4FVgD631ZGAq8B/gzkQEKoRIjra2NtpGw/+oCyHEOGa3Q2GhGehrtw8+0Ndmg4MOauWhhzZz4421uN0QDsO11+axbl3/022UAq/XJFTNzRbr1tmor5eEKulMC0dTnXK7we8fHQmLwwFXX20qVNOnw/nnw4EHwurV6Y4sI8Q1R0optRaYDLQAjwNnAS9qrbe8I7TWG5RSfwfuTUSgQojkuOmmm9IdghBCiE5ud8+BvsGgmfNq9fPRt9MJ221n1un9978Onn7ay6OP+pgzp4UTTvBTVNT3/pdYhUprU6FqbFTk5WlycqLYZcpo8jgcUF5ulvlt3hzLbNMdFey0EzzyCPy//wdLl0JhoTkejfb/5hNxD+T9FDgfeFJrPdBH2c8A34jzGkIIIYQQ406salRRESEQUNTVWSilcLv7Hugbs+OOHaxYsZm77srm0UezeO45H0ceGeCEE5pwu/u/liRUaZCdDR4P1NRAU5NJptL9A1cKDjkE5swx/x0Ow09+Yjr8jYbq2SgU129Ma33IEL8vAATiuYYQIjXOOussAJYsWZLmSIQQQnRnWWagr9cboanJoqnJwmbT/SZFYOZVnX12E/PnB7jllhz+9S8Xp55qntOafhOx7glVU5NJqCZMiJKdrdN+fz9m2e1m31R2tqlOhUKjozoVe5O0tMB//gMfftj1nMtl/iwuNjGPcyNZ2tdfXxkNNGL2R92stX4nzthGPaXUYmBR7+OlpaWpD0YIIYQQY5LdDgUFUbKzo9TXWwQCFi6XxuHo/zVlZREuuaSB9naTkDU3K04+uZC5c1uYM6el3+QoVg2LRqGx0UZDgyY/3yRUNlty/n7jXlaW2TtVWwsNDaZSNdAvN1Vyc/sf3FtVldpYRimlB+uz2deLlDoduAJYDawAqjFt0I/q/PMmYEfgEOAwrfWziQo4E8ycOVO/n6YhZ2+/bf7cU3olCiGEGOVWr16NYzTcMGaYYFBRW2vR3m7uuYeS4Kxfb+OSS/L55BMXkyd3cPLJfvbbr23ApYJgtsgEgwrQWypUklANTdGRR+LxeIb3omi0K3mxrP5LiKnyzgD1kL33Tl0cQ/Hqq0k5rVLqA631zL6ei3f32AzgFa31d7TWV2itb9NaXwrsDLwHlGutfwbcBlwc5zWEEEIIIUQvsYG+RUVRQiFoaRl8/lRFRYTbbqvhmmtqsdngwgsLOO64IlpbB75RtyxToXK7oaHBxtdf22hsVLJlJlksy3QcsdvNviQZljuqxbvq9XDgxN4HtdYRpdStmATqWOAfwNHxhyeESLaFCxcC0r1PCCEyiVKQnW32TzU3m4YUNpvC4+k/o1IK9t47yF57BXn2WS+ffurE6zXfX1NjUVjY/017LKGKRqG+3iz5kwrVwCrvv59p06bFf4K2Nti0yTR98PnSU50qL+//uSRVgDJJvBWpMDCpn+dKgGDnf3sAmZ0txCjm8XiGv/RACCHEqGCzQV6eZvLkCF5vlEBAEQoN/po5c1o5//xGwCz7++lPS7noonzWrx84K7Is8PlMhaq+XipUSeXxwDbbQH6+aZc+2C9WpFy8Fal7gcVKqTpMC/QmpZQbOBDTfOFxpVQ+cAbwbmJCFYOprISPP4Zg0Pxb2267gT9IEALg2muvTXcIQgghRsjhgKKiKDk5UerqLJqbLTyeoXXcmzAhyoIFzSxfnsUrr3j4yU9aOP54PwUFA1eofL6uClVjoyYvTypUCWdZMHGiaUixaRMEAqmtThUWmhbtvRUXp+b6o1y8zSYcwB3AAkyXvhDgBBTwEbAvsDuwHDhAa/3vRAWcCdLRbKKyEr74wnx4kZtr3t8ffQQ77CDJlBBCiNFJmk0kh9bQ1qaoqbGIRMy9wVBmqtbWWtx5Zw6PP+7D69U88cQmfL6h3SdGo2YlmmVBXl6UnBw97ue4dnR0jGxpX2/RqOnqV1Nj2pDHWpGnwhFHmD1bb7xhpkCPIwM1mxh2RUopZcd05DsNuBzYGygE2jBJ1Ctaa62Uegso01oH+z3ZGBUIdHXPS5WPP4ZddgG/3zwsy/wP51NPmeMivUZzF8Vf//rXANx6661pjkQIIUQidB/o6/crGhqGNtB34sQo553XyPz5zXz8sWtLEvXSSx722qttwPt2U6Ey99r19RaNjWamVVaWJFQJY1lQUGB+0Js3m+V+Pt/QsmSRFMNOpLTWYaXUm8CxWusVwJf9fJ8M4k2hYNBUovz+rmO5uea4EAMpKChIdwhCCCGSwFSHNFlZERobzUBfu33ggb5gOvxVVLQC8OWXDi68sIDi4jAnnujnoINaB1y6Z7N1JVS1tRb19ZJQJZzbDVOmQGMjVFebdZ2D/VJFUsS7tO9GYDowR2sdTnhUGS4dS/tefRV22qlrGeuOO5q5bp9+Cvvsk9JQhBBCiCGRpX2p1d4ODQ1moK/bPfBA3+7ef9/FzTfn8tlnTqZO7eCUU5qYNSs4pG06kYj5UNeyxl9ClfClfX1pbzfVqWAwudUpWdqXmKV9nVqA/wO+VEq9DLT2el5rrc+M89wiDtttZ/ZExfZI1dZ27ZESQgghhHC5oKQkSlubprbWIhAY2kDfmTPbueOOal55xcPf/pbD5Zfn89hjm8jKGvzDeKlQJZnLBZMnd1Wn7HbzSxUpEW8idRTQhGkusX8fz2tAEqkUijWUeOop86HEjBnSaEIMzbHHHgvAXXfdleZIhBBCpILHo5k0KUIgYOZPaW3mTw1UYVIK9tuvjVmz2vjqKztZWZpIBK6/Po/DDmth+vSBp910T6hqakxCVVAQxeeThGrElIIJE8wPuKrK7PPIypK9UykQVyKltZ6a6EAykVJqMabdew+lpaWpDwaTNMUaS4zm5gZidKmoqEh3CEIIIVKs+0DfpiaLxkaFzWYaUgzEbodp08yujnXr7Dz/vJcVK3z86EetnHSSn7KygQdK2WzmHj8SgepqC7vdVKgkoUoApxMmTTKJVFWV+WFLdSqp4toj1eMESjmBfKBOay3Dd0nPHqmYWLdASaSEEEKMdrJHavTo6DDd9gIBC6dTD3kbjN+vuPfebB58MJtIBA4/PMApp/jxeod2fxkOm20+NlvXkr9UjUhKtpTsker/4mapX6yz30iHe8keqT73SMWd+yulDlFKvYfZH7UB+K5S6jGl1PHxnlMIIYQQQqSewwHFxVHKyyMopQkEFOEhtBPLydEsXOjnkUc2MWdOC++958bpNEnUUD6rt9vNfb7dbipUX39to7lZDem1YgAOB5SVmUcwCK292xmIRIgrkVJKzQEeB6qBs7ud59/AUqXU0YkJTwiRbAsWLGDBggXpDkMIIcQo4HZrysujFBdH6OiAlhYzB3YwRUVRLrigkWXLqrDbobVV8YtfFLN8eRah0OCvt9vNkr/uCVUgIAnViCgFOTkwdarJVv1+hpQdiyGLtyK1CLhLa30wcHvsoNb6MuAW4KwExCaESIHtt9+e7bffPt1hCCGEGCWUgqwszeTJEfLzo7S1QVvb0JKa2Kovv99iwoQIS5bkceSRJTzzjJfIwNungJ4JVVWVxfr1klCNmN0OpaVmM30oJNWpBIp3jlQbcJTW+gmllA9oBvbSWr+tlDocuFdr7UtwrBlD9kgJIYQQg5M9UpkhHDbzp5qbzUBfl2tor9Ma3n3XxU035bJqlZNp00LcfHMNublDv/cMh83KNIfDdPnzejNnD1Va90j1Jxw2Q0ebmsDrNUnWUAUCsM02skeqm3jbn1cCOwFP9PFcKaY1uhBCCCGEyHB2OxQWRsnJiVJX1zXQt/s9+Mcfm5vrXXbpWsenFOy+ezvf/W41L73k4d133eTkmCSqqspGcfHgJapYhSochk2bbDidmokTo4O2axf9iFWncnLMIN9QyHT2kx9mXOJd2ncLcKFS6hggL3YupdSewAXAsgTEJoRIgaOOOoqjjjoq3WEIIYQY5VwuKCuLUloaIRIxBYqhLNezLJg9u42LLmpAKZNEzZ1bwjnnFLBmzdA+07fbTbt2m80kVBs22GhtlSV/cfP5THUpO9t09uuQxtvxiHeO1J+UUkWY/VGxFPYVTGL2CPD7xIQnhEi2GTNmpDsEIYQQGcTr1VRURGhuVtTXm4G+Wg+9qJGTE+W44/wsW5bNz39ezEEHtXLiiX5KS4daodJbKlQulyY/XypUcbHZoLjYJFOx6pTXK9WpYRjRHCml1BRgNlAE+IE3tNYfJSi2jCV7pIQQQojByR6pzBeJQFOTxZtvurHZYObM9iG/tqnJ4p57snn44SyU0jz55Gby8obQIrAbs4dKjcqEalTukepPJAL19VBXB2533/ugZI/UVuLdIwWA1nod3br2CSGEEEKI8SM2SLeoKILfrwgETFIzlPw4NzfKGWc0MW9egPfec21Jol56ycMeewSHNNQ3VqHq6ICNG2243ZqCgihu9+hJqDKCzQaFhWZD2qZNJmny+aQ6NYi4EqnOTn3nA7OAXLbea6W11t8eYWxCiBT42c9+BsCKFSvSHIkQQohMZbrqacrLI9TWWjQ3KzwePaSmcMXFEebMMS2516+3ceGFBeTnRzjuOD8/+UnLkJIyhwMcDpNQVVZKQhU3j8dUnWLVKZdr3FWghiPeZhO3AxcBUWAl8EGvx78TEp0QIun22GMP9thjj3SHIYQQYgwwA30jlJTEBvqqIQ30jamoiHD77dVss02Ya6+dwLx5JTz/vGfI53A4TFMKMBWqykqLtjbJpIbFsmDiRJgyxXwdCCBdPfoW79K+OcAVWus/JDIYIUTqnXPOOekOQQghxBgSG+jr9ZrlfvX1Fpalhlwd2nlnM2/qnXfc3HxzDlddNYHdd28f1v6prgqVYuNGC7c7Sn6+xuORhGDI3G6YPBkaG6G6WpKpPsSbSAWBDxMZiBBCCCGEGDssC/LyND5fhMZGC7/fwuEY2kBfpWDPPYN873tB1q61k5cXRWu47ro8DjywlZ12Cg1+EroSqlBIUVlp4fFEKSjQuN2SFAyJZUF+vtkvVVtrvhZbxPvTuA/4ZSIDyURKqcVKKd37sXHjxnSHJsSQHXrooRx66KHpDkMIIcQY5XCYgb7l5WHsdk1zsyIcHtprLQu+8Q3zzZs22XjxRQ/HH1/E736Xz1dfDb0e4HSaJX/RqGLDBrPkLxiUJX9D5nJBeTlD2vQ2jsT70/gCWKCUeg94Amjt9bzWWv95RJFlAK31YmBx7+MzZ86UjzlExth///3THYIQQohxwO2G0tIobW2amhqLYNCMLRpqkaOsLMKKFZtZvjyL++/P5vXXPRxySAunndZETs7Qbr2cTnA6e1ao8vOlQiXiE28itbTzz//f3r3HR1ld+x//rMn9QgIhyE0IIgq2WrFcWqlW2gOt4kFbteK1IGq1vLRq/XmO5fRiPdpSPFjqpWprW/pDtLUtv7aoVUEKPy16DtKieCmo3BQMkITcCCEh2eePPQNDDCFPEvLMZL7v12teyTzzzH5WhmFgZa+9dl9gTCuPO6DHJ1IiPcFNN90UdggiIpIizA5u6Ftba5SXRzBr//qpvDzHNdfUcOGFe1iwoBerVmUfSIKam9uflCmhkq7QoUTKOacCSRERERGhtDSN9euzaGyMsH9/hJKSBgYMaGrzOZEIFBT4hhRVVRGqqiKkpTmys9t3zT59mrnllipuuKGKjAzYtw9mzjyGyZP3csklte1OiGIJ1b59PqHKzW2mT5/mdschqa3dCZGZLTazU1sc+5iZZbc49mkz+++uClBEjq5zzjmHc845J+wwREQkCZWWprFlSw6nnmp8/vPNnHKKsWVLDqWlae16fno69O3bzLHH7icry1FbazQ2tv/6sT2mamsjDBjQxEMPFXLhhQNYvDiv3euwwC8B8hv7Gtu2pfPhh770UKQtQWakvgQ8HLsT3ZR3HXAmsCruvEJgbJdEJ4FNmBB2BJJspk6dGnYIIiLSTdasaUfLvADWr8/i1FON2toItbV+pik/31i+PI+RI/cFHm/fPqiqitDYCGPHNpDWvnyMvn2bmTevnLVrM3nwwUJ+9KM+PP54Pg8/vIvi4va3Tc/Kgqys2AxVOnl5zfTurRkqaV1nW2+o3YlIkps1a1bYIYiISJJqbIxQUNBMbe3BYwUFjsbGjq0CycryHf727oWGBmhuNnJz27d+CmD06AZ+9rNdvPRSNn/9aw59+/ok6sMP0xg4sO1yw5ZxKKGSI1EPQxEREZEUMWZM8FmituzfH6FvXzvQ5GHEiEbKy41Ro1ynr9XUBDU1viFFerq1e92TGZx5Zj1nnulr8yoqIlx6aX9OOaWBWbOqOOmk9tcOxhKq+vqDCVWfPs3t2gtLej41jRBJcZMmTWLSpElhhyEiIkmopKSBdesiVFYazc1QXm6sW+cbTnRWWprf0Hfo0Cays5uprTUaOjBsfn4z111Xzfr1GcyY0Z/Zs4vYujXYXEJ2tl9DtW+f8cEH6ZSWRtjXtTmpJCHNSImkuGnTpoUdgoiIJCnfnW8vy5fn0dgYYdQoR0nJ3iN27QsiIwP692+msLCZ8vIItbURsrNdu/eGzcyESy+tZerUPSxa1Isnnshn5cocFi8upX//YHH60r5DZ6h699YMVaoKmki1NqeqhvsiSezaa68NOwQREUliAwY0HWgs0dWlg/Gys2HQoIMb+u7bBzk57d87Kj/fcd111Vx0US1/+1v2gSTqhRdyGD++nl692v9fWiVUAsETqUfMLLacMPa2XWhmdXHn5Hc+LBERERGRQ8Vv6FtdbezeHWxDX/Ad/s47z//XdefOCN/5ThF5eY6vfrWar3ylNlBDiYMJVYRt2yLk5/uEKjMz8I8mSShIIvXrVo6t7qpARCQcEydOBGDFihWhxiEiItJekYhfP5Wf30Rlpd/QNz29/Rv6xpSWpjN79m6WLcvlgQd68+ST+Vx7bTVTptS1u3QQIDvb4Rzs3RuhpiZCr17N5OUFi0WST7vfIs65q45mICISjhkzZoQdgoiISIekp0NxcTO9ejWze7dPYnJy2r9+CmDo0P3Mn1/GmjVZPPhgAf/1X304/fR6+vVr//5T4GfL4hOq3bsz6NsX+vRBM1Q9lJpNiKQ4JVIiIpLssrJgwAC/fqqszG8OnJNDuzf0Bb++6xe/2MWmTen069eMc3DvvYVMnFgfaO1XLKGKRPz+WpWV0Lu3EqqeSO3PRVJcY2MjjY3t31NDREQkUeXkOI49toljjmmmoQHq6gwXoC2aGQwfvh+A8vIIK1bkMGtWP26+uZgNGzICxWLmk7n8fKipgU2bYOdO0D+5PYcSKZEUN3nyZCZPnhx2GCIiIl3CDHr18g0pCgubqauD+vp2dqKIU1zczO9+V8qNN1by5puZXHllf77znSIqKoL999k3yPAJVXU1bNyohKqnUGmfSIq75pprwg5BRESky6WlQVFRM/n5sHu3338qM9MFKq/LzoYrrqjl/PP3sHBhL5YtyyUnx09xNTe3v/U6HEyonPMJVXzJX0awyS5JEEqkRFLcFVdcEXYIIiIiR01mZmxDX0dZmQXe0Bf8DNesWdV87WvVpKfD/v1w1VXHcMYZ9Vx+eQ35+e2vH2yZUO3eDUVFPqlSQpVcVNonkuLq6uqoq6s78okiIiJJLDvbMXhwM/37N9HYCHv2+FmlIGLJ1549xtCh+/nlLwu44IIBPPFEPvsC7kUcX/JXVeVL/srKfJImyUGJVCeY2R1m5lretm/fHnZoIu02ZcoUpkyZEnYYIiIiR50Z5Oc7hg5toqiomX37YN++YA0pAAoLHXffXcGCBTsYObKR+fN7c/HFA9i+PUCbwLiYYglVZSW8954SqmSh0r5OcM7dAdzR8vjYsWMD/nUUCc/Xv/71sEMQERHpVrENfQcMaKa62qirM9LTHVlZwcY56aRG7r+/jNWrs/jLX3IZMKAJgO3b0yguDtZNIr7kr7ISKioOlvwFKUOU7qM/FpEUN23atLBDEBERCUVaGvTp4xg8eD/l5ZEOrZ8CGDduH+PG+dq+6mrjq1/tz7BhDdx5J4wfH2ys+IRq924lVIlMpX0iKa6qqoqqqqqwwxAREelWpaVprF+fxRtv5LB2bQ6RiDFwYBNNTVBbC01NHRs3N9cxa1YV27en8+Uvw/Tp8Pbbwccxg7w8n1Tt3u3XUFVUqOQvkSivFUlx559/PgArVqwINxAREZFuUlqaxpYtOZx6qlFQ0Ezfvsa6dTmUlOzl2GObqK01yssjgJGT47AA21Clp8MFF+xh8uQqXnhhGA8+CF/4AqxcCcOHB481EvEJVXMzlJf7W9++UFjoZ9QkPEqkRFLcN77xjbBDEBGRJDdmTMCWde2wZk3ABUsBrF+fxamn+lbotbU+WcnPN5Yvz2PkSP+zNDX57nw1NRHS0oLtPwV+/dQNN8Dll8Pzzx9Mop56Ck4/3SdDQSihSjxKpERS3AUXXBB2CCIiIt2qsTFCQUEztbUHjxUUOBobD656SUvzx3JymqiuNvbujZCREXz9VJ8+EFuOXFEBN93kx77+evja13y3viDiE6qyMp9QFRdDQYESqu6mREokxZWVlQFQXFwcciQiIiIHHY1Zrpj9+yP07WtEonnTiBGNlJcbo0a5w163vt4oK4tQX2/k5rojJi2NrTTtKyqC556DH/0I5s2DBQvg5pvhiisIPOPlZ9F8QrVrl0+qlFB1LzWbEElxF110ERdddFHYYYiIiHSbkpIG1q2LUFlp0VI5Y926CCUlDYd9jt/Qt4n+/ZtoaPBlf0E39AUYMQJ+/nNf4jdyJNx5J5SWdvxniSVUOTk+odq0yTen6GizDGk/zUiJpLhbb7017BBERES6ld/vaS/Ll+fR2Bhh1ChHScneA/tAHY4Z9OrlyMvz5X4VFREiEd+QIqjTToMnn/Qb8A4d6o/dcQeceSZ8/vMEanABH52hiq2h0gzV0aNESiTFTZ06NewQREREut2AAU0HGksELSOMbeibl9dEZWWE6mq/firohr5mfoYK/CzS88/72apPfxpmz4YxY4KNF4utZUIVK/mLqBatS+nlFElxpaWllHampkBERCRFZWRAv37NDB68n/R0R02NdXifpz59YMUKuPtuePddOO88uPpq+PDDjo0XS6iysmDnTr8PVVUVHSpHlNYpkRJJcZdccgmXXHJJ2GGIiIgkrexsGDiwmUGDmti/32/o25GEJTMTZsyAVavgttvgrbd8hz7o+JqntLSDCdWOHX4NlRKqrqHSPpEUd/vtt4cdgoiISNIzg9xcx5AhTdTUGDt2GHV1vglE0PVOeXm+m98NN/gNfpuaYOpU+NSn4MYbffe/oGIJVVOTT6hiXf569VLJX0fpZRNJcWeffTZnn3122GGIiIj0CJEIFBY6hgzZT2Ghn53au7djY8X2rNq713f4e/RRmDABfvITqKvr2JixhCoz03cL3LQJqqs1Q9URSqREUtz777/P+++/H3YYIiIiPUp6OvTrB8cd50v/qquh4fDd1duUnw8//jEsWwannw5z58JnPgPvvNPx+NLS/GxUZqZfh6WEKjiV9omkuCuvvBKAFStWhBuIiIhID5SZCYMH+xmknTuhpgZyczvWknzkSPjVr2D1ali40CdpAJs3+xbqHSnRiyVUTU0+ocrI8CV/+fkq+TsSJVIiKe7b3/522CGIiIj0eLm5UFLiE6mdO/3MT15e8PVTAOPG+RvAnj1w/vkwcCB861vw2c92bMxYQrV/vy/5S08/uIaqI+OlAuWZIilu0qRJTJo0KewwREREejwzv5/Tccf5zXI7s34qJjsbvvtdqKyEyy6DadNg7dqOj5ee7mej0tMPlvzV1IALvudwj6dESiTFbdy4kY0bN4YdhoiISMpIS/OJ1PDhfqaqpgb2BdsT+JCxLrwQVq6EO++Et9+Gc8+FN9/sXIzp6X42Kj0dtm/35YNaP3UoJVIiKW7mzJnMnDkz7DBERERSTkaGL8krKfEJUU0NHd7QNyvLb+D78stw773w8Y/7408/3fFNfeFgQtXQ0PG9rHoqrZESSXHf//73ww5BREQkpWVnw5AhvtRv1y6or/czVR1p9pCf78v7wI/3zW/65Ozqq2HWLOjdu2tjT2WakRJJcWeddRZnnXVW2GGIiIikNDM/8zNsmG+bXlfnb51Zm5SfD88/D1OmwE9/6vegevDBzq/LEk+JlEiKW79+PevXrw87DBEREcHPQvXp49dPFRT4cr/6+o6PV1IC998Pzz0HY8bAnDnwwQddF28qU2mfSIq77rrrAO0jJSIiqWfMmA52eOgG6elwzDFQWOjL/aqrISfHr6vqiI9/3O89tXGjT9IA7roLTjvNz1ipxXlwSqREUtwPfvCDsEMQERGRw8jKgmOP9WV+O3Z0bkNfOJhE1dXB8uXw0EMwerTfg+qMM7ou7lSg0j6RFDdhwgQmTJgQdhgiIiLShtxcv35q4EDfKn3Pns6tn8rNhaVLfYe/nTt9g4rLLoOtW7ss5B5PiVQnmNkdZuZa3rZv3x52aCLt9sYbb/DGG2+EHYaIiIgcQfyGvn36dH5D37Q0n0C9+KLf1HfTJt+gAjrehj2VmNM2xV1u7Nix7tVXXw07DJF2mThxIqA1UiLS/d59910yOrrgQyTBNTY2MmLEiKN6jYYGKC+HqipfApiV1bnxmpp8cuUcfOlL8LGPc+E4MQAAEi5JREFUwS23+LVaNTW+LDDV/sqa2Rrn3NjWHtMaKZEUd88994QdgoiIiHRAZqYv9evd25fn1dT4hhTpHfwffmzdVX29T6Iefxx+9zu49lq48squi7unUGmfSIobN24c48aNCzsMERER6aCcHBg6FAYNgsZGX/LX3Ny58X74Q1ixAiZPhvvug7VruyzcHkMzUiIpbm30k3H06NEhRyIiIiIdFdvQNy/Pl/rt2uX3pMrJ6Xhr8+OO8139br3Vl/fJoZRIiaS4m2++GdAaKRERkZ4gtqFvr15+/VRlpV/XlJ3d8TFHjPBlg3IoJVIiKW7+/PlhhyAiIiJdLD0d+vc/uH6qsxv6ykcpkRJJcSrpExER6bmysmDIEL/v1I4dvpFEZzb0lYPUbEIkxa1evZrVq1eHHYaIiIgcRXl5fkPfAQN8MtXZDX1FM1IiKe+2224DtEZKRESkp4tEoLDQb7q7e7dfQ5We7kv+JDglUiIp7oEHHgg7BBEREelGaWlQXAwFBb67X3W1b0aRmdn6+atXQ12d35BXDlIiJZLiTj755LBDEBERkRBkZsLgwbB3r18/VVOj9VNBaI2USIpbtWoVq1atCjsMERERCUlODpSUwMCB0NDQ+Q19U4VmpERS3OzZswGtkRIREUllZr7UL35D37Q0P0MlrVMiJZLiHnnkkbBDEBERkQSRlgZFRYdu6NvYGHZUiUmJlEiKGzlyZNghiIiISILJyPCt0nv3hrfe8jNWciglUiIpbuXKlQCcddZZIUciIiIiiSY7G/r392un0pU5HEIvh0iK+973vgdojZSIiIi0zgyysjQr1ZISKZEU98tf/jLsEERERESSjhIpkRQ3XLvriYiIiASmfaREUtyyZctYtmxZ2GGIiIiIJBXNSImkuLvuuguASZMmhRyJiIiIJJpt2+D116G+3jecOOEEGDw47KgSg2akRJLQxIkTMbM2bzNmzGjXWAsXLmThwoVHN2ARERFJOtu2wT//CZ/4BEyeDCef7O9v2xZ2ZIlBM1IiSWju3LlUVFQcuD99+nT69+/P3LlzDxwb3M5fFw0ZMqTL4xMREUl1f/97Njt3hh1F57z+uk+iqqv9LRKBnBxYssQfTyQTJnT/NZVIiSSh8ePHH3I/JyeHoqIizj777MBjPfvsswAdeq6IiIj0XPX1UFjok6iYwkJ/XJRIiaS8OXPmAEqkREREutInP1nPiBFhR9E5DQ1+M95IdDHQSSdBWRmMHh3ODFCiUSIlkuJ+85vfhB2CiIiIJKATToDXXvPlfIWFPol67TUYNSrsyBKDEimRFDdgwICwQxAREZEEFFtuvWSJL+cbPdonUera56lrn0gP9t5773HGGWdw4oknctppp/Hqq69+5JwlS5awZMmSEKITERGRRDd4sG8sMX48TJyoJCqeEimRHuz6669nxowZbNiwgblz53L55ZfjnDvknHnz5jFv3ryQIhQRERFJTirtE+mhdu3axSuvvMIzzzwDwOTJkwFYs2YNY8eOPXDe73//+1DiExEREUlmmpES6aG2bt3KoEGDyMjIOHCspKSErVu3HnJecXExxcXF3R2eiIiISFJTIiWSQlqW9QEsXryYxYsXhxCNiIiISPJSaZ9IDzV06FC2b99OY2PjgVmpLVu2MHTo0EPOu++++wC44IILuj1GERERkWSlREqkB9i8efNHjvXr14/x48ezYMECrr32WpYuXYpzjjFjxhxy3p/+9KduilJERESk51AiJdKDPfzww0yfPp177rmH3NxcFi1ahJkdck5hYWFI0YmIiIgkLyVSIj3YCSecwKpVq9o857e//S0A06ZN646QRERERHoEJVIiKe6hhx4ClEiJiIiIBKFESiTFxfaZEhEREZH2S7n252ZWZGbvxt1PN7OfmtkGM3sn+n1amDGKdKfc3Fxyc3PDDkNEREQkqaRUImVmPwD+BmTHHZ4ODARGAScDY4Ep3R+dSDgee+wxHnvssbDDEBEREUkqCZtImdmTZvZAK8dvNLP3zKzOzF41sy8GGPZF4J4Wx/4BzHbONeMTrBygssOBiySZRx99lEcffTTsMERERCRBTZjgb3KohEykzOxMYGorx78OzAd+C1wGvAE8bWafas+4zrm/AMtbHPu7c+5tM5sDlAG1wP907icQSR5Lly5l6dKlYYchIiIiklQSKpEysyvN7B3g/3No+R1mlgF8D3jQOTfbOfdH59wMYC3wnc5e2zl3O1AEVAO3dHY8kWSRkZFBRkZG2GGIiIiIJJWESqSATcAvgG/x0fK604D+wKIWx58B/iXaNOLfzKyslduph7ugmf2nmZ0B4JyrAf4MDOuin0ck4S1YsIAFCxaEHYaIiIhIUkmoRMo595Jzbo5zbg5Q1eLhT0S/vtHi+Jv42ashzrm5zrniVm6vtXHZeuB687KB8/BrqT7CzFx7bsF/cpHwKJESERERCS6Z9pEqAhqcc3taHN8d/VrQwXHnAz8D3gUagD8Aj3dwLADWrFlTZmZbOjPGYQwCth+FcTsjrJi647pH6xpdOW5XjDUGwMzWdD4cCVkifkaEJZlfi0SLfUz0a5d/RgwbNmxEWlraYWuLKyoqMoqKihq7+rrdcY2uHrez47333nu5AMcff3xdV8UkbWtqamrcvHnzu0c+M7BE+4wIy1H7bGqh5HAPmHOJOYFiZpuBp5xzN0Tv/wfwXedcVovzvgA8B4x0zm3o9kC7kZk555yFHUe8sGLqjuserWt05bhdMVZsFjXR3lsSXCJ+RoQlmV+LRIs9zM8IfdZ33Xj6rO85Eu0zIiyJ8J5OqNK+I9gNZEbL7+IVRr+WdXM8IiIiIiKSopIpkfpn9OvJLY6PAj5wzlV0czwiIiIiIpKikimRWoXv5HdJ7ICZpQMXA0+HFZSIiIiIiKSepGk24ZyrN7O7gTlmVgOsA2YCxwHnhxqciIiIiIiklKRJpKLmAQbcCBQD/wD+xTm3MdSous/3ww6gFWHF1B3XPVrX6MpxE/E9IeHR++GgZH4tkjn2rqbP+qM3niQvvRcSRMJ27ROR7pEIXW9EJHHpM6Jn0J+j9DSJ8J5OpjVSIiIiIiIiCUGJlIiIiIiISEAq7RMREREREQlIM1IiIiIiIiIBKZESEREREREJKNnan4uIiEg3MLMIkNvKQ43OuX3dHY+ISKLRjJSIYGZPmtkDYcchIgnlTKCmldtPwgxK2qe1z3Uzi5jZf5rZB2a2x8xWmtm4sGIUORIzO8bMFplZuZnVmNnfzezKuMdDfU8rkRJJcWZ2JjA17DhEJOEMBzbiE6r42z1hBiVH1sbn+g+BW4H7gSvxifELZjas24ITCWYxcC4wB7gKeBf4v2b25ejjob6nVdonkqKiv9H5LjAi7FhEJCEdD7zlnHsp7ECkfdr6XDezfsBNwGzn3L3RY08Bm4H/A9zQfZGKHJmZjQY+A1zgnPt/0WN/AE4BvmJmLxHye1qJlEjq2gT8Ivr9v4cZiIgkpOHAewBmlumcawg5Hjmytj7XPwdkAYtiB5xzDWa2DPhi94QnEkgR8N/AgV/mOOecme0BckiA97RK+0RSlHPuJefcHOfcHKAq7HhEJOEMB04xsy3APjPbYWa3m5mFHZi07gif658AypxzO1ocfxMYbmb65bokFOfccufcp51zu8ws28z6mdk3gdHAEyTAe1p/aURERKQ1xwMOuB3YDpyPX4+Qiy8fk+RSBFS0cnw3/hfreeiXapK45uDL+ACeBZYAEwn5Pa1ESkRERFpzM/Cyc25j9P6zZpYF3GZmc5xzdSHGJsFlHOHxvd0ShUjH/AR4Ht/w5lbg5xz5PXvU39Mq7RMREZGPcM4tikuiYpYA2fiyP0kuu4HCVo4XAjVaAyeJzDm3yTn3jHPuW8B84FIS4D2tREpEREQOYWYjzOyK6AxUvNisRss1CZL4/gn0j3bvizcKeC2EeETaFN0fakMrD72Jz2EqCfk9rURKREREWhoILASmtDh+GfBP59yu7g9JOuk5oAmYFjtgZr3w+009HVZQIm3YBJxgZh9vcfxz+F/m/JqQ39NaIyUiIiItvYxvO/wrMxsBbAHOwzecuDDMwKRjnHPbzOxhYI6ZZQAfALcA+4GHQw1OpHWPA3cAz5jZPUApcA4wHbghEd7TSqRERETkEM65/WZ2HnA38E2gD/A2cLlzbnGowUln3ALswe8xlQ/8Dficc64y1KhEWuGcqzezLwLz8B1DM/GfQzOcc7+Onhbqe9qcc91xHRERERERkR5Da6REREREREQCUiIlIiIiIiISkBIpERERERGRgJRIiYiIiIiIBKRESkREREREJCAlUiIiIiIiIgEpkRIRkYRiZivMzB3htsDMhkW/vyjEWH9sZv/RxuMTozEO66Z4TjGztWamfSJFRI4yfdCKiEii+TegKO7+r4Ed0eMx26LHzgH+0X2hHWRmJwNXAMeHcf3WOOfWmdk7wE34TSxFROQo0Ya8IiKS0MxsM7DZOTcx5FAOYWZPAW865/69jXMmAn8FjnPObe6muE4Hnotes7w7rikikopU2iciIkmpZWlftNxvrZldaGYbzKzWzBab2TFm9rCZ7TazSjO7z8wiceMMNrMnoo9VmdnjZnbsEa49EjgXWNDi+E1mtsXM9prZX4HhLR5PM7O7zGyTmdWb2WYzu9fMsqKPL4kmjvHPSTezcjP7afT+yWb2bDTeajP7s5kdFzvfOfcy8CFwVeAXVURE2k2JlIiI9CTHAz+I3n4IfAl4EzgBuB54GLgRmAZgZgXASuAUfDncjdHvXzSzwjaucy6wyzn3duyAmV0OzAdWANOBt4AHWjzvP4BbgEeAS4BfRa/5/ejji4ASMxsT95zP40sdHzezbOBZfGn+zOhYI4E/m5nFPWclMLWN+EVEpJO0RkpERHqSfOBfnXPvAJjZxcBg4Dzn3B7gt2Z2KXA68AQ+eeoLnOCcK4s+5wVgE/BV4P7DXOezwN9jd6JJzLeBPzjnpkcPP2lmucCMuOcdB/zQOTcnev+PZvYpYGz0/p+BWuBCYE302EXAFuBvwOjozzPTOfd89NrrgNlAH6Ai+pw1wAwzy3bO1R/pRRMRkeA0IyUiIj3J+7EkKmo38Ho0iYqpAXpFv/8i8DJQb2b5ZpYPVAGv4ZOtwxkIlMXdHwSMwidn8f4Qf8c5d5Vz7i4zG2RmnzWza4AxRH+x6ZyrA/6IT6QwszT8rNpvnF/U/C6wE3jEzL5hZic75/7HOfcl51xF3KXKgAxgQBs/g4iIdIISKRER6UmqWzlW18qxWBlcf3znv5oWt7FAcRvXKcLPHMUMiX4tbXHerkMuanZOtKveVnxZ37lAy4YQi4ATo10BzwL6AY8DOOdqgAnA8/guhuvM7H0z+2aLMaqiXwva+BlERKQTVNonIiKprBq/5ujuVh6rbON5tRyapMRmvI5pcd7A2Ddm1gf4XfR2mnOuNnr8L0BO3HOW4WedLsQnem86516PnpsGfOCcuy56/0Tg68A8M1vnnFsaHSM24xY/SyUiIl1IM1IiIpLKXsR31nvZOfeSc+4lYC1wB37m53A24meKYt7Gzyy13Bz44rjvjwfygIVxSVQBvrTvAOfcfuDJ6FhfBhbGPTwNqDKzvOi5G4A7o48NjTuvL9CET8hEROQoUCIlIiKpbB6+hO/PZnapmU3Hzwh9Ani6jee9DIyJdcqLJj93AZeZ2SNmdrGZPYhvShGzAT8Ddle0Rfss4BX87NaQaNOJmEXAyfgZrkVxx18A9gF/MbMro800noyOsTTuvE8CrzjnGoK8GCIi0n5KpEREJGU5594HzsD/e/gLfPvyncBZzrltbTx1MdAbODVurPnArcC/4tc/DQauiHu8Gl+uVwA8BlyN77Z3NX4GaVbcua8AHwDPOec+iDu+IzpGNvBzfBv1DOALzrmtcfF9Fniq/a+EiIgEZb4JkIiIiARhZk8DG5xztxyFsY8FNgMXOuf+FPC5Y4G/AsfFWrqLiEjXUyIlIiLSAWZ2En7j2xOdc201pggyZi98F8HL8KV9I51zTQHH+A3whnPurq6ISUREWqdESkREpIPM7E6g2Tl3RxeNNwR4D78B70zn3IsBnz8S+BMwWhvxiogcXUqkREREREREAlKzCRERERERkYCUSImIiIiIiASkREpERERERCQgJVIiIiIiIiIBKZESEREREREJSImUiIiIiIhIQP8Lfw4vLAmjWfkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#This is the final plot\n", "fig, ax1 = plt.subplots(figsize=(12,6))\n", "\n", "ax1.errorbar(hess_data.T[0]-start_mjd,hess_data.T[2],yerr=hess_data.T[3],\n", " xerr=hess_data.T[1],fmt='rs',label=r\"H.E.S.S.\")\n", "\n", "ax1.errorbar(lat_data.T[0][mask_hess]-start_mjd,lat_data.T[2][mask_hess]*1e-3,yerr=lat_data.T[3][mask_hess]*1e-3,\n", " xerr=lat_data.T[1][mask_hess],fmt='bo',markerfacecolor='white',label=r'Fermi-LAT $\\times 10^{-3}$')\n", "\n", "ax1.errorbar(lat_data.T[0][~mask_hess]-start_mjd,lat_data.T[2][~mask_hess]*1e-3,yerr=lat_data.T[3][~mask_hess]*1e-3,\n", " xerr=lat_data.T[1][~mask_hess],fmt='bo',markerfacecolor='white',alpha=0.3)\n", "\n", "\n", "ax1.set_ylabel(r\"Energy Flux (erg cm$^{-2}$ s$^{-1}$)\")\n", "ax1.set_xlabel(\"Time (days)\")\n", "ax1.set_ylim(1e-13,2e-11)\n", "ax1.hlines(0.,0.6,40.,linestyles='--')\n", "ax1.set_xlim(0.35,35.)\n", "ax1.set_yscale('log')\n", "\n", "plt.xscale('log')\n", "\n", "ax1.fill_between(xmoon-start_mjd,1e-13,2e-11,color='lightgray',alpha=0.7)\n", "ax1.legend()\n", "\n", "#Plot best fit line\n", "ax1.plot(lat_data.T[0][2:]-start_mjd,fit_func(lat_data.T[0][2:]-start_mjd,popt[0],popt[1]),'b--')\n", "ax1.plot(hess_data.T[0][2:]-start_mjd,fit_func(hess_data.T[0][2:]-start_mjd,popth[0],popth[1]),'r--')\n", "\n", "#Shaded error band\n", "ax1.fill_between(fx1,f1,f0,color='blue',alpha=0.1)\n", "ax1.fill_between(fx2,f2,f3,color='blue',alpha=0.1)\n", "\n", "ax1.fill_between(hx1,hf1,hf0,color='red',alpha=0.1)\n", "ax1.fill_between(hx2,hf2,hf3,color='red',alpha=0.1)\n", "\n", "ax1.vlines(time_vmax-start_mjd,ymin=1e-13,ymax=2e-11,colors='k',linestyles=':')\n", "\n", "ax1.text(0.38,0.9e-11,r\"$\\alpha_{\\rm LAT}~\\,= 1.31 \\pm 0.07$\",color='blue')\n", "ax1.text(0.38,1.2e-11,r\"$\\alpha_{\\rm HESS} = 1.43 \\pm 0.18$\",color='red')\n", "ax1.text(time_vmax-start_mjd-0.12,1.2e-13,r\"T$_0$\",color='k')\n", "\n", "ax1.get_xaxis().set_major_formatter(plt.FuncFormatter(format_func))\n", "plt.tight_layout()\n", " \n", "plt.xticks(ticks=[1,5,10,30],labels=[\"1\",\"5\",\"10\",\"30\"])\n", "\n", "#increase size of ticks \n", "ax1.tick_params(which='major',length=10,width=2)\n", "ax1.tick_params(which='minor',length=6,width=1)\n", "\n", "plt.savefig(\"JointHESSLAT_lightcurve_RSOph.png\",dpi=250)\n", "# \n", "None" ] }, { "cell_type": "code", "execution_count": null, "id": "restricted-graphics", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "cheap-hungarian", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0" } }, "nbformat": 4, "nbformat_minor": 5 }