{ "cells": [ { "cell_type": "markdown", "metadata": { "collapsed": false }, "source": [ "# SNIFS R-channel\n", "\n", "**Author:** Yannick Copin \n", "\n", "We construct a very simple optical model mimicking the red channel of the Supernova Integral Field Spectrograph." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "%load_ext autoreload\n", "%autoreload 2\n", "%matplotlib inline\n", "#import mpld3\n", "#mpld3.enable_notebook()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import warnings\n", "warnings.filterwarnings(\"ignore\")\n", "\n", "import numpy as N\n", "from spectrogrism import spectrogrism as S\n", "from spectrogrism import snifs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Load the SNIFS-R channel optical configuration:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "------------------ Optical configuration 'SNIFS-R' -------------------\n", " name : SNIFS-R\n", " wave_ref : 7.6e-07\n", " wave_range : [5e-07, 1e-06]\n", " telescope_flength : 22.5\n", " collimator_flength : 0.169549\n", " collimator_distortion: 2.141\n", " collimator_lcolor_coeffs: [-4.39879e-06, 8.91241e-10, -1.82941e-13]\n", " grism_prism_material: BK7\n", " grism_prism_angle : 0.301592894745\n", " grism_grating_rho : 200.0\n", " grism_dispersion : 2.86\n", " grism_grating_material: EPR\n", " grism_grating_blaze : 0.261799387799\n", " camera_flength : 0.228014\n", " camera_distortion : -0.276\n", " camera_lcolor_coeffs: [2.66486e-06, -5.52303e-10, 1.1365e-13]\n", " detector_pxsize : 1.5e-05\n", " detector_angle : 0.0\n" ] } ], "source": [ "optcfg = snifs.SNIFS_R\n", "print(optcfg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Load the simulation configuration:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------- Simulation configuration u'standard' ----------------\n", " name : standard\n", " wave_npx : 6\n", " modes : (1, 0, 2, -1)\n", " input_coords : [-0.01 -0.005 0. 0.005 0.01 ]\n", " input_angle : -0.174532925199\n" ] } ], "source": [ "simcfg = snifs.SNIFS_SIMU\n", "print(simcfg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Create a `Spectrograph` instance from optical configuration:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------- Spectrograph ----------------------------\n", "Collimator: f=0.2 m\n", " Geometric distortion: center=(+0, +0), K-coeffs=[ 2.141], P-coeffs=[]\n", " Chromatic distortion: lref=0.76 µm, coeffs=[-4.39879e-06, +8.91241e-10, -1.82941e-13]\n", "Grism:\n", " Prism [BK7]: A=17.28°, tilts=+0',+0',+0'\n", " Grating [EPR]: rho=200.0 g/mm, blaze=15.00°\n", " 1st-order null-deviation wavelength: 0.76 µm\n", "Camera: f=0.2 m\n", " Geometric distortion: center=(+0, +0), K-coeffs=[-0.276], P-coeffs=[]\n", " Chromatic distortion: lref=0.76 µm, coeffs=[+2.66486e-06, -5.52303e-10, +1.1365e-13]\n", "Detector: pxsize=15 µm\n", " Offset=(+0.000, +0.000) mm, angle=0.0 deg\n", "Spectrograph magnification: 1.345\n", "Central dispersion: 2.77 AA/px at 0.76 µm\n" ] } ], "source": [ "spectro = S.Spectrograph(optcfg)\n", "print(spectro)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Test the optical model:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "--------------------- SUCCESSFUL ROUND-TRIP TEST ---------------------\n" ] } ], "source": [ "if not spectro.test(simcfg.get_wavelengths(optcfg), verbose=False):\n", " print \"ERROR: backward modeling does not match.\"\n", "else:\n", " print \" SUCCESSFUL ROUND-TRIP TEST \".center(70, '-')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulate spectra on the detector and plot:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAArQAAAJoCAYAAACXwnKOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8leXBxvHffUZ2QhIIYROmDAUEERwMRQW3YqvVDlve\n9rVW7bK21qq1tX3do2pttUOttmqts06WgCB77z0SRhjZ+4z7/SMBAkkgwDnnyTm5vp8Pn+bczxlX\nKJ5cuc/93I+x1iIiIiIiEq1cTgcQERERETkVKrQiIiIiEtVUaEVEREQkqqnQioiIiEhUU6EVERER\nkaimQisiIiIiUU2FVkRERESimgqtiLQ4xpjzjTFzjDFFxpj9xpgvjDHDjDE3G2OCxpifHXX/XGPM\n6Lqvf22MebXesaAxptQYU1L3vwX1jt1jjNlSd2yHMeb1Y2S62Rjjr7tvkTFmuTHm2nB8/yIicmJU\naEWkRTHGpAL/Bf4AZACdgd8A1XV3KQB+boxJPsbT2KO+HmStTbPWplprM+te52bg68CF1to04Cxg\n2nHifVn3POnAH4F/GWMyTuw7FBGRUFOhFZGWpi9grbX/trWqrbVTrbWr6o6vBeYCdzbz+Uzdn6Od\nBXxmrd1G7Qvutdb+9QRyvgrEA71O4DEiIhIGKrQi0tJsAALGmJeNMROMMelHHbfAfcCPGzl2IuYB\n3zLG/KxuOUOz3w+NMW5gElAErD+FDCIiEgIqtCLSolhrS4HzgSDwIrDPGPOeMaZ9vfusAKYAv2jm\n0y4xxhQaYwqMMU/XPcc/gTuAS4AZQL4x5ufHeZ5z6tbgVgKPAlfW5RUREQep0IpIi2OtXW+tnWSt\n7QYMpHYd7dNH3e1+4Nb6RfcYzrTWZlhrM621P673Oq9bay8B0oHvAw8aYy6uOynt4IlkK+s9z9y6\nNbjpwAc0v1CLiEgYqdCKSItmrd0AvExtsa0/vh54B/gVR54E1pjG1tDWf66AtfZtYAVwurV2dt0J\nZGnW2jMauX8F8ANgzMHdFURExDkqtCLSohhjTjPG/NQY07nudlfgRmrXvB7tt8B3qJ0xPdHXudkY\nc5kxJsXUuhQYAMxvzuOttYXULon45Ym+toiIhJYKrYi0NKXACGC+MaYU+JLamdMGuxrU7VDwKtDc\nLbzqKwHuAbYDhcDDwPettV+eQNangbHGmEEn8BgREQkxY+3xPqkLcwBj2gB/BU6n9iSQSdSe5fwm\n0B3YBlxvrS12KqOIiIiItFwtYYb2D8DH1tr+wGBgHXA3MNVaexowHX2kJyIiIiJNcHSG1hiTBiy1\n1vY6anwdMMZam2+M6QDMsNb2cySkiIiIiLRoTs/Q9gD2G2NeMsYsMca8aIxJArKttfkA1to9QHO2\n5RERERGRVsjpQusBhgJ/tNYOBcqpXW5w9LSxswt9RURERKTF8jj8+nlArrV2Ud3tt6kttPnGmOx6\nSw72NvZgY4yKroiIiISctfaY+1eHWzuTaA9QFcmX3G6tzYnkC4aSo4W2rrDmGmP61m2ePg5YXffn\n28AjwM3A+8d4jkhEbdUeeOABHnjgAadjxDT9HR+2v+pd4lzZpMWdG9Ln1d9x+IX673hb6QHe2rKE\nHwwYTbI3PmTPG+30bzn8jHG0ywJwgCr+bi6M2OtNstO7R+zFwsDpGVqAHwL/NMZ4gS3UbpLuBv5t\njJlE7R6R1zuYT0QiqG38NS3ih4k4z+tyk+yJw2WcXh0nIi2d44XWWrscGN7IoYsinUVEnKcyKwd1\nTk7nBwPHOB1DRKKAfu2V4xo7dqzTEWKe/o7DT3/H4ae/48jQ37NIQ45fKexUGGNsNOcXCbUD+8pp\nm3Wsq8CKiMjxGGMcPynMGGMjvIbW8e/5VGiGViRGLFmQy0+/9w57dpU4HUVOgLUBAsFSp2OIiEQ1\nFVqRGDFwUAd+8LNRZHdMdTqKnIBy/wr2Vr2JtUGno4iIRC3HTwoTkdCIT/Ay4vwcp2PICUr0nIbX\n1Q6jM/lFRE6aCq2IiIPcJgm3O8npGCIiUU1TAiIiIiIS1VRoJWr4gwGWH8gjqJ0tREREpB4VWoka\nG4r38vyaWewsL3I6irQARYWV/OpH/2Xzhv1ORxEREYdpDa1EjQEZHfnNsMvpkNTG6SjSAsTHu+nW\nI5OU1Dino8gJKvUtwlo/aXEjnY4iIjFChVaiisqsHJSYFMctPz7P6RhyElzEYbWrg4iEkAqtiIhE\nVLJ3kNMRRCTG6FdkEREREYlqKrQiIiIiEtVUaEXCwOcLMH/2tlN+nv1VZdw17x22lOhM/tauutrP\nn574gh1bC5yOIiLS4mgNrUiI+XwBnn1kJssX7WTvnlKu/MoZJ/1cKd54zm6fQ9uE5BAmlGgU8AfZ\nu6eUygqf01FERFocY6N4k3pjjI3m/BJ7/P4gzzw8g+WLdmLag90L13/rTC6feLrT0UREpJmMMVhr\njcMZ7N/NhRF7vUl2uuPf86nQDK1ICHk8LrLap+AZY0g4C6o+hqzsVKdjiYiIxDQVWpEQS70knkGl\nAa7sUM4rianE9XA7HUlERCSm6aQwkRAblBnkykXLWHrGu9xctoWuWv4qIiISViq0IiHmffcz5ty6\ngpK9QRZe8wVmyQKnI0WVqTvXMX/vVqdjSAuwcukuPnpnldMxRCQKqNCKhNCeGUuYccsU/DW1tytK\nLFOueIHSrbudDRYBC+Zs56Xn53GqJ2ruKCtgV3lxiFJJNNuxtYAVS3Y5HUNEooDW0IqEUMXuIgIB\nFxA8NObze6jMLyS1R0fngoXZgjnb+fOTXxAIWIJBy6TbRmLMyZ0sO+m0c0OcTqLV5RNP1w4hItIs\nmqEVCaGeN17IqJd/gXHX/qflTUng4k8fpf3IAQ4nC5/83aX8+ckvCGZZ3MMMs6ZuYurH652OJSIi\nrYhmaEVCrNfXLwJg3h3PctGH/0f2uQMdThRe2R1TueS7/ViQuJb26UFKu6Uz9pI+TscSEZFWRIVW\nJAx6ff0iulw2gviM2N+DNmgt2zrs5HtJB0hYvYe3hwxgb00pnb3pTkcTEZFWQksORMKkNZRZAAPc\nnLGXTRM+Yc7XFzJx1UqyE3UFv+baV1nG2sI9TscQEYlqKrQickpqDhSw6MK32b/dj68avvzuEva/\nN93pWGFXXFTJmy8vJhAIHv/Ox/D5rvW8smFeiFKJiLROKrQickpm3PAg+7f5Dt32VcPUb/+Nsh35\nDqYKr+KiSh6+dwofv7eGPz85+5RK7bU9hnDPmRNCmE5EpPVRoRWRUzL415PwJMUfMXbardeQ0i3b\noUTh5fcHeeT+KezaU4z7DMOCOdv5xwsnf/EMr8tNWlxCCBOKiLQ+KrQSVTZv2E95WbXTMaSeDqMH\ncfHHDx8qtaffdQPDH/lfh1OFj8fj4qKrTyNuoouuV1sSL/Rw4YS+TscSEWnVVGglagSDlofuncyU\nj7THaUtzsNQO/tU3YrrMHrS9SwFXnlvF/1Rvpf8FNexKLnI6kohIq2ZO9TKVTjLG2GjOLyduZ24R\nbdslk5DodTqKtGL55bPIveV51r61j3OePYOO3/kxyd4cp2NFjSq/jzi3G5fRnIq0TMYYrLUnd7nD\n0GWwfzcXRuz1Jtnpjn/Pp0LvJhJVOndNV5kVx+V9/2VW/msffh98eftKCl/7yOlIYVdUWMlv7vqY\nTev2nfJz3b/4Qz7YvjIEqUREaqnQioicgJUPv8aKf24/dDvgh+m3vsfeuasdTBVeRYWVPHzfZLZs\nPMBjv5l2yqX2a72GcU77HiFKJyKiQisickI6jDsLb1rSEWNtBvSgzWldHUoUfv/620J255fgvdxQ\nHefjT09+gd9/8luVDW3XjeyktBAmFJHWToVWpIXTOvGWJWt4P8ZPeexQqc0c0osJUx8nPjN2C9oN\n3x1G6k1ezhvtJ+VGF9/82dl4PPrxIdIaGGMmGGPWGWM2GGN+0cjxdGPMO8aY5caYecaYAfWOtTHG\nvGWMWWuMWW2MGRGunHpHEmnB/L4Ad3z7LWZM3uh0FKnnYKntcMGQmC+zAGurdnPp0HJGvjabmwcU\ns9aV63QkEYkAY4wLeA4YDwwEbjTG9DvqbvcAS621g4GbgWfqHfsD8LG1tj8wGFgbrqwqtCItmMfr\n5rJrB9Knf5bTUeQoWcP7cem0J2K+zAIMbVtN2598yrzfb6X42o8Y367E6UgR4fcH+fyzDfqURFqz\ns4GN1trt1lof8AZw9VH3GQBMB7DWrgdyjDFZxpg0YJS19qW6Y35rbdjePFRoRVq4y64ZSOeu6U7H\nkFZs1U1Ps/a9AgByl1ax5KI/E6ipcThVePn9QZ5/fBYv/2k+Lz0/T6VWWqvOQP2PZPLqxupbDkwE\nMMacDXQDugA9gP3GmJeMMUuMMS8aYxLDFVSFVkREmrTlzc9Z/faWI8Zyl5Sw/HevOZQo/Ky1PP/4\nLBbPy8XVDWZO2cQrf57vdCyRkFpnC3nPbjn05xQ8DGQYY5YAtwFLgQDgAYYCf7TWDgUqgLtPMXaT\nPOF6YhERiX45XxnNjvfGsvXNGYfGskb25/Sf3eBcqDAzxtBnQBbLk/NIPwuKpkPfAe2djiUSUv1M\nBv3IOHT7A7utsbvtpHbG9aAudWOHWGtLgUkHbxtjtgJbgGQg11q7qO7Qf4AGJ5WFimZoRUSkSS63\nm9Gv3UOPG8YCtWX2kk8fIS4t2dlgYVY9yM+Ii/z8uMM++lweJHialhxIq7QQ6G2M6W6MiQO+BnxQ\n/w51Oxl4677+HjDTWltmrc0Hco0xfevuOg5YE66gmqEVkRblrVeX4o1zc80Ng5yOInUOltr0/t0Z\n8OPrYr7MApydFaD4sS+Z/cRmbnhlOG36ne90JJGIs9YGjDG3A5OpnQT9m7V2rTHmltrD9kWgP/CK\nMSYIrAb+p95T/BD4Z13h3QJ8J1xZVWhFpEWJi3fjjXM7HUOO4nK7GXL/t5yOETFlT77B/Ac3Yy18\n+c2FXPRGP7iyp9OxRCLOWvspcNpRYy/U+3re0cfrHVsODA9rwDomms/cNMbYaM4vIiItz5Z/TmbW\ntx6h/o+XuGQ3l87+E5mDezkXTCLGGIO11jicwf7dXBix15tkpzv+PZ8KraEVERGpJ75dBq447xFj\nnvR0PMkJDiUSkeNRoRUREamn8/jhjHvvQdzxtaU2qXNbLv38KdJ6H739poi0FCq0IieopibgdAQR\nCbODpTa1d2eVWZEooEIrcgLWrtzDXd9/ly0b9zsdRWLEkvm5TP14vdMxpBGdxw9n4pqXVGZFooAK\nrUgzrV25hyd/N52igkoe+/VUlVoJiU0b9rFy6S6nY0gTXB7tuCESDVRoRZrprdeW4ksOEPcNF5Vx\nPt55fbnTkSQGXP/NofzkVxc4HUNEJKqp0Io00413DqPNN9xcPbya9G+5uf72IU5HEhEREVRoRZot\n17+XSd320+7Hn/DtbvvYXrXX6UgiIiKCCq1Isw3z7GbnFR+z9v0Cdl76CYNcO52OJCIiIqjQijSL\nr7ySOWOeJX+DD4C9m33MHv00vrJyh5OF34zJG9m0fp/TMSRGWGuZO2srlRU1TkcRkRiiQivSDJtf\nm0r++iPL695NlWz8+2cOJYqMaZ+s5+U/zePx30xj8waVWjl1hQcqeOHpOaxYop0dRCR0VGhFmqHf\nLVcy4IcTjxy79Sr633GtQ4nCb94XW3n1xQWYHKiK9/HYA9PYs7PE4VQS7TLbJfPUXycy4vwcp6OI\nSAzxOB1AJFqMePo2ANY88w79br2Kkc/9EGOMw6nC5/TBnWg/MgXPmDLKSgy9NnegfcdUp2NJDMjI\nTHI6gojEGBVakRMw4unbyB51Bt0njorpMgtQ4a0hcVwF36zeTnlOKu9mJmOxQGx/3yIiEn1UaEVO\nUM51o52OEBFt4w2TKrYy7+oZJLd18z/TvoJLXVZERFograEVkUaVLl7NvKtnUFFs2bfFz/Jxb1NT\nWOx0rLD7+N3VFBVWOh1DREROgAqtiDRQdaCYyVc/TEWxPTS2b4uP6Vfeg7X2GI+Mbq+/tJg3X1nC\nw/dNVqkVEYkiKrQi0kBC2zbk3HjJEWPGZej7g+tidu3wu68v59P31+DqZ9idX8LD902mqtLndCwR\nEWkGraEVkUad/cStAKx+6j8Yl2HUK3fT6+sXOZwqfIaN7Mrk3LVkXehjz2bDMNuVhESv07FERKQZ\nVGhFpElnP3ErxuUic0ivmC6zAIVpFfQcX8P1BzazbmRXNpeXOh1JRESaSYVWRI5p+GO3OB0hIvq2\nMSS/v4xZty6nxwVtuPY/reP7FhGJBVpDKyICFLzxCXNvXY7fBxsnF7Plq393OlLY7csv46F7J7N/\nb5nTUURETokKrYi0evuXbGD6//4bf71zwDZ8tp9Fdz7nXKgw25dfxsP3TmbdqnyVWhGJeiq0ItLq\npffrRruRpx8x5k1NpNtXL3QoUfj9+cnZ7C8rx3uZYX9JOX96crbTkURETpoKrYi0ep6kBC7+8P/o\nMHYwUFtmL/nsUdqPHOBwsvC56bZhJN/oZtRYP0k3uvn6bWc5HUlE5KSp0IqIcLjUdrvmvJgvswBr\ngju5aVAxZ/5pBtefWco6u9vpSCIiJ027HIiI1PEkJTDund86HSMiRmcUsXDcR2xcUkmPxQWc/+lg\npyNFxI6tBbTNSiE5Jc7pKCISQpqhFRFpZQI1PhZd/Ay5S2ov77t1djlLL3uaYDDocLLw2rrpAA/d\nO4VHfz2F8rIap+OISAip0IqItDJr/vAOOxYVHzG2ZdY+Nv7tE4cShd/WTQd49NdTqXTXsG1bgUqt\nSIxRoRURaWUG/PBaulw24oixnOvH0mfSBIcShZ+vJkAwNUDy1w3eq13UBPwEArE9Iy3SmqjQioi0\nMu74OC58+4FDpTbn+rGM+ec9uNxuh5OFT0aPJNre7OGW/sVcen4VWd9OISVN62hFYoVOChMRaYXc\n8XFc+M5vWPvsuwz48XUxXWYBUr2Gb6XvYf2oD+k+pi19nrwJl9Gcjkis0H/NIiKtlDvOy+l3Xh/z\nZRYgmL+N9RM+pCA3wMrX9uL76b+x1jodS0RCRIVWRERiWk1xGVMvfICC3MChsRWv7WL5PX90MJWI\nhJIKrYiIxDQbCGKT0huMB73JDqQRkXBQoRURkZgWn5nGhGlPkDm456Gxwfd+g6G//Y6DqUQklFRo\nRQC/P8jalXucjiEiYXKo1A7ppTIrEoO0y4G0en5/kD89/gWLF+TyvR+ey3ljex7/QSISdeIz07h8\nzrN4EuOdjiIiIdYiZmiNMS5jzBJjzAd1tzOMMZONMeuNMZ8ZY9o4nVFi08Eyu2jeDmhr+cszXzJn\nxhanY4lImKjMisSmFlFogR8Ba+rdvhuYaq09DZgO/NKRVNIqWCyecYa0bxpcmpwVERGJOo4XWmNM\nF+Ay4K/1hq8GXqn7+hXgmkjnktbB43HR/roURl7k5ydd99P1BsgYnOh0LBERETkBjhda4CngLqD+\nDtfZ1tp8AGvtHqC9E8GkdRjWvoZz/zmbeYP/yzfzN9AxSZuti4iIRBNHC60x5nIg31q7DDDHuKsa\nhoRN4Im3WPjQVqrKYMHX5uCb9YXTkSLiwL5ypyOIiIiEhNO7HJwHXGWMuQxIBFKNMa8Ce4wx2dba\nfGNMB2BvU0/wwAMPHPp67NixjB07NryJJaZsfvkj5v1+BQevgFldAVO/8ncu/XIEGQNzHM0WTksW\n5PL8Y7P4+neHc8H4vk7HERFx1IwZM5gxY4bTMRq4/vZjzfWF1qRnI/ZSYWFayrWsjTFjgDuttVcZ\nYx4FDlhrHzHG/ALIsNbe3chjbEvJL9Ep98O5TL/u1wR9hy+Jmdy1HZfOeJrUHh0dTBY+Sxbk8sdH\nZ+F3BzE1cPP3R6jUiojUY4zBWhu5Ntl4Blt2x7iIvV7Ks9Mc/55PRUtYQ9uYh4GLjTHrgXF1t0VC\nrusV53Dh27/B5XUDsV9mAZbMzyWQHST+uwbTx7Dwyx3oF0MREYlmTi85OMRaOxOYWfd1AXCRs4mk\ntThYahfc+Scu+fSRmC6zAKO+0YvNa7dyU58y3s1IYnSnXhgTtb+Ui4iItNgZWpGI6nrFOVy76u8x\nX2YBXKac75htbBn2NjdtW01CXLHTkURERE6JCq1IHZe3xXxgEVZtt6xg+RWfU5AXYP63FpL1xQyn\nI4Wdzxdg755Sp2OIiEiYqNCKtCLlefuYNv7PlBXWrpmtqYIZX/+Q/fNXO5wsfHy+AM8+PJPf/fIz\nduVqNlpEJBap0Iq0IqWbd1F51PazNdVQuDbXmUBhdrDMLl+yk+LqSh6+f4pKrYhIDFKhFWlFOowZ\nzEUf/B53grd2wBjOe+En9Pn2BGeDhUlZaTU7dxbhvdSQcJOhzFaxe6cKrYhIrFGhFWllOl00jIs+\n+D2epHjOe+En9P3u5U5HCps2GYn0urUdF15Qw/8OKSX7ewn0GNLW6VgiIhJireMsGBE5QqeLhnHd\nptdI6pDpdJSwchnDsA7FuO+eyp7ZBdz06VXE6dd4EZGYo7d2kVYq1sssgA0G8d79FitfyWfvZh/b\nx3+Ap3in07HCbu7MrZSXVTsdQ0QkYlRoRSRmrbj3eZa/crjA7tviY9q4+/FXVDmYKrymfLiOPz81\nm0fun6pSKyKthgqtiMSsqsqGb3E+k4y/MjaL3rRP1vPaXxfi6gHbdxTwyP1Tqaz0OR1LRCTsVGhF\nJGaNeOoHDPjhxEO3M87IYcL0J0lo28bBVOHTNSeDhCEeOn4N4iYaOvdsQ3y8TpUQkZNnjJlgjFln\njNlgjPlFI8fTjTHvGGOWG2PmGWMGNPexoaR3OhGJaSOevg2A3Z8vY8L0x2O2zAL4OwTofrWPGws3\ns+7crmwu8+NyGadjiUiUMsa4gOeAccAuYKEx5n1r7bp6d7sHWGqtnWiMOQ34I3BRMx8bMiq0IhLz\nRjx9G77SCrypSU5HCaucVLh+xQpm37KUnFFpXPX295yOJCLR7Wxgo7V2O4Ax5g3gaqB+KR0APARg\nrV1vjMkxxmQBvZrx2JDRkgMRaRVivcwCFP3nM+beshR/DWyaVsKO617CWut0rLDy+4M8+/BM5s7a\n6nQUkVjUGah/Kcm8urH6lgMTAYwxZwPdgC7NfGzIaIZWRCQGHFi2iWn/8zr+msNjG6cV4P3R04x4\n5ifOBQsjvz/Ic4/OZOmCPJYsqP25ec7oHg6nEokOs/IK+WJnYSie6mHgD8aYJcBKYCkQCMUTnwgV\nWhGRGNDmtK60P28Qu6YuOTTmTUkg52uXOJgqvD77YA1LF+bhGW8Irrf89ZkvOW1AezLbJTsdTaTF\nG90lg9FdMg7dfmhBo59y7KR2xvWgLnVjh1hrS4FJB28bY7YCW4Ck4z02lLTkQEQkBngS4xn3/u/o\ndNFQoLbMXvzpo2SfO9DhZOFz8ZX9aP+NZEaP89PpazDhh/1VZkVCayHQ2xjT3RgTB3wN+KD+HYwx\nbYwx3rqvvwfMtNaWNeexoaRCKyISIzyJ8Yz74Pd0nzgq5ssswK6qYoYPKeaMx6fyndJNFGTtdjqS\nSEyx1gaA24HJwGrgDWvtWmPMLcaY/627W39glTFmLTAe+NGxHhuurCaaTxgwxthozi8iIifPF8hn\n6Td+xOo395GYAue/N54uF/7c6VgSA4wxWGsd3fPOGGPL7hgXsddLeXaa49/zqdAMrYiIRKXlkx5h\n9Zv7AKgsg9nXTKZo0XKHU4mIE1RoRUQk6hSt3c6q11ceMVZZZln02387lEhEnKRCKyIiUSe9f3fG\n/utejOfwj7G2w/ow6h+/dDCViDhFhVZiWllpNcGg1lmLxKKcr4w5VGrbDuvD+CmPE5+e4nQsEXGA\n9qGVmFVUWMnD906mV992/M8d5+qa9iIxKOcrY7g4NYl2I/qrzIq0Yiq0EpMOltnde0vYvbMEQKVW\nJEZ1Hj/c6Qgi4jAtOZCYtGbFbnbvLyHueoNnnGHpojwO7Ct3OpaIiIiEgWZoJSYNPrcLnUnhkl4H\n2Fruxnd+Du3a6wpCIiIisUiFVmKSL+jjqp4HKJv4HoMHp1P46wwsoAUHIiIisUdLDiQmJQVKKZv4\nHrmLKln+t91k/vZdTCu4qpzfH6S8rMbpGCIiIhGlQisxJxgI8MVlvyJ3ceWhsWV/zWPVfc85mCr8\n/P4gf3r8Cx6+bzJlpdVOxxEREYkYFVqJOb7ickoONJyN3b+5GBujs7QHy+yieTvYkVfII/dPUakV\nEZFWQ4VWYk58ZhqXTn+S9IHdD43lXD+WMa/dgzGxuYp2z64SVq3YhXeCIe4rhtzdRWxYu9fpWCIi\nIhGhk8IkJiVkpXPp9Cf55MKfkj6wB2NeuweXx+10rLDp1LUNA3/cgbYZW+mU4Gdax2z6D+3gdCwR\nEZGIUKGVmJWQlc5lM5/Gm5Yc02UWwGI5o1MBqb+aSuGqEq5891r8QT/gdTqaiIhI2GnJgcS0+My0\nmC+zAC4bJO1X77Ly1XzyllZScu27xPuLnY4Vdju2FcbsumgREWk+FVqRGLD8F8+w4tXdh27nLa1k\n1oRfEfQHHEwVXiuW7OS3d33MS8/PU6kVEWnlVGhFopy1luJdlQ3Gy0td+MoajseCFUt28sxDM/An\nB5k5ZZNKrYhIK6c1tCJRzhjDqH/cDcCW16cDkDm4J+OnPk58eoqT0cJmz64S/NlBUr9iqFhee9tX\nEyAuXm9pIiKtkd79RWKAy+0+VGqL1mxj/NTHSWjbxuFU4dNzVDu6ZLm4OXMvMzukk+zKVJkVEWnF\n9BNAJEYcLLX+skri2sTmzOxBGfEBbi7awOKL53DWNzsR/2BPpyOJiIiDtIZWJIa43O6YL7MAdv5c\nFn11NpVllqV/2onv/n85HSnsCgsqWDxvh9MxRERaJBVakTAoKa7ik/fW8MbLi1m5dJdOWAqh4g25\nTL3mRarKDo8tfn4Ta5583blQYVZYUMFD907muUdnMXfmVqfjiIi0OFpyIBJiudsKeei+yZSX1gDw\nyXtrGDWuF9+941yHk8UGb1oycR3bU1mcd2jM5XWTlNPFwVThU1JUyUP3Tia/sBRXN3jxD3NwuQ0j\nzs9xOpo29FJvAAAgAElEQVSISIuhGVqREHvr1aWHyuxBX0zbzKb1+xxKFFuSOmQyYfpTtOlXW2Bd\nXjdjXr+PnImjHE4WHkkp8WT3TCXpa4ZON1oST/fQuWvsnvAnInIyVGhFQmz9mr2Nj69ufFxO3MFS\nmzGoZ0yXWQC32+AbV8V3+hfynb0b6HC1n0CmlrCIiNSnJQciIdaufTJ524saHZfQSeqQyVWL/4zL\nHduXNjbGcFO7AlZe+BEb8nzc9LKfTkMudzqWiEiLohlakRC7fOLABmMdO6cxbERXB9LEtlgvswC+\n4hJWjvsXezf78FXDvG8vpODDz52OFXbrV+fz+kuLnY4hIlFCM7QiIXbumJ54vW4m/3cdRYWVnHFm\nR66+YRAeb2yXr62bDpC7vZDuPTLp3jPT6TgxY8ZXH2DvpsNrsn3VMOUbf+HqleeSmtPBwWThs351\nPk88OJ3qKj8+X4Bv/e/ZTkcSkRZOhVYkDIaf253h53Z3OkZEBAJBnn/8CxbNPbxH6ojzu/P9n5yP\ny60PgU7V6fd8iz1zfom/8nCp7T3p8pgtswX7y3nywenUZPjx5BimfbyerPYpXHrNAKejiUgLpp82\nInJKZk3ddESZBZg/ezuzZ2xxKFFs6Th2CBd99BCexDgA+t9xLSP/cLvDqcIns10y593Yk+yb4MzL\nA7S9KonRF/V2OpaItHAqtCJySpYuzGt0fFkT43LiDpba0++6IabLLIC1loKe+Uyq3MLYN76g78gy\n9tuy4z9QRFo1LTkQkVOSlBzX6HhyE+NycjqOHULHsUOcjhF2xhgm7lvLl1/5gsoyy3nVQTo9cqXT\nsUSkhdMMrYickgsu6YMxR44Zl2HMJX2cCSRRrXDBMr68ZgqVZbV77S56Yjub7n/B4VQi0tKp0IrI\nKTltYDa33jmK7I6pAHTolMZtPxtF79OyHE4WPnv3lPKffy7l5T/PZ9nCPKzVhQ5CZe5P/3aozB60\n4LHPKduR71AiEYkGJprfiI0xNprzi8Sammo/cfGxvZJp3ap8nvjtNGpqAofGxlzcm0m3neNgqthR\nlruXTy/4CaVb9gBgPC7Gvn4fOdeNdjiZtCbGGKy15vj3DGsGW3bHuIi9Xsqz0xz/nk+FZmgl6gQC\nQacjSBNivcwCvPnK4iPKLMDMKZvYsbXAoUSxJaVreyZ8/hSpPTuozIpIs8X+Tx+JKX5fgOcem0Vm\nu2Rtti4RFwwE2bLxQKPHNm/YT7ceuqBEKKR0bc+lM56mYNkmul6hmW8ROT4VWokaB8vs0oV5ULfS\nRKVWIsnldpGVncK+/IbbSLXvkOpAotiV3CWL5C6xuw5bREJLSw4kaixdmMfShXl4Lzd4xho+/3QD\nWzc1PlsmEi5XXX9Gg7He/bIYMCg2r9x10L78Mgr2lzsdQ0SkUZqhlagx7Jyu9PxBJgN65lPiM5QN\ny6FH77ZOx5JWZvS43iQnxzHtkw2UlVYzeFhnLps4EHP03mUxIn93CS88NYfNG/YD0P+MbG75yflk\nZCY5nExE5DAVWokaZb4azupbQNuffoo33cuOe6+iJuAnzq1/xhJZw0Z2Y9jIbk7HiIg/PDSDnTuK\nD91euzKfF56azd0PXuJgKhGRI2nJgUSNFI+ftj/5L2vfPcCKl/bQ7Xcf4jWtY9u2ooIKpyNIK7R1\n04EjyuxBa1fmc2Cflh+ISMuhQitRY9Gk/2Pte4fXzK54aRcr7nrSwUSR8dkHa/nFD95n/WptLC+R\nFQw2vUWets8TkZZEhVaigq+8kl1LGha6vHl5BH1+BxJFxmcfrOVff19ElfHzxIPTVWolonr0bkdW\ndkoj421bxa4OunCPSPRQoZWo4E1OZML0J8k4PefQWNbI/lzy6SO4vLG5htbvDzL3i624BxkSv2uo\nyfSz8MvtTseSVsTlMtz+89G0a598aKxTlzZ8/yfnO5gq/KZ9sp47v/cO35n4Gr+/5zO2bNzvdCQR\nOQ5d+laiStX+Yj698Kd4UhIZ/9mjeFNb7pnW1lr8NojX5T7p55iybS3LihZyZbsSXs1N5zsDLqF3\nm/YhTClyfMGgZcuG/bjchh6928bsjg4Ac2Zs4cWn5xwxlpjk5eE/Xk16RqJDqSTSdOnb6BObU1sS\nsxLatWHC9Cdxx3tbbJkNBIN8tGMVs/ZspNRXTYfENC7tOpCR2T1O+LmyUkq58u35LH9wPTe/fh5x\n7lFhSCxybC6XoXe/1nGRg2mfrG8wVlnh48sZW7js2oEOJBKR5tCSA4k6Ce3atNgyC/DKxnl8lLuK\nUl81AHsqS3hpw1xm79l0ws+V8uoHzL97HRUllsXXzyZhwZehjtvilJdVs31LgdMxpJUqK6ludLy0\niXERaRlUaEVCKL+yhAV7tzV67KMdqwja5p8ZvuOdGcz5+VyCgdrbVeUw7dqXKVobu+toy8uqefTX\nU3n43sls3rDP6TjSCg0a2qnR8cHDOkc4iYicCBVakRDaWnqAplZ1F1RXUFRT2fwn88TB0etvE5Ig\nRteNHyyz23ILqEzy8fgD01RqJeKuuWEw3XpkHDF2yRX96Hd6tkOJRKQ5tIZWJITSvAlNHvMYF0nu\nuGY/V7erzuWCtx7g868+QNAXICErjUs/f4r0ATkhSNry2CD4TZCErxhSsi1lb1sCgdgs79JypaTF\n85snLmfl0l3s31vGaQOz6dIt3elYInIcKrQiIdQvPZu28ckcqG54FaVhWd1I8HhP6PkOltq5tz/D\n+E8fjtkyC5CU6iXlRg/j00pon1/AP7/Rizbdm/4FQSRcXC6jJQYiUUaFViSEXMbF9weM4rnVMymu\nt7ygR2pbbug57KSes9tV59LpkrPwJDR/djfSFu/bwZSda9ldUUJWQgrjOp/GOdk9T+g5DIavdihk\n0/gPWZZbw7fehQz1WRERaQbtQytygoqLKikuqqJz1za43Y0vQ/cFAyw/kEdhdQXdUjI5LT1219/N\n2LWB1zcvajB+dffBXNat+dsc+cormDbyq+xeXQVAYqph/OT7yRgxOmRZW6K5M7fSd0B72mYlH//O\nIhIR2oc2+miGVqSZfL4ALz0/j7kztxIMWtIzErn51hEMPbtrg/t6XW7OyuruQMrI8gUDfLhjZaPH\nPs1bzQWd+pLYzGUWMyb+6lCZBagstXx25eNctXwgSR3bhiRvSzNjykZefn4eWdmp3P3gxSq1IiIn\nSbsciDTTe2+sYM7nWwgGaz8VKCqs5I+PzqJgf8P1sq3FrvLiQ/vtHq064Gdb6YFmP1fv/7ka4zny\nLanLVaNJ7JB5Shlbqjmfb+Hl5+dh+hj2lpby8H1TKCk6gV0wRETkEBVakWb6cuaWBmN+f5AFX8bu\nvrDHc7zZ1+bOzgL0uH4sY1771aFS22fSpZz3lztj9jKrPXq3JfmsOHreECTlRkPXAemkpMY7HUtE\nJCppyYFIMzW1XPsErpUQc9onptIjtS1bG5mJ7ZiYRk7qiS0V6HH9WAB2T1/KOX/6cYsus7vKi6gK\n+OmakoH36P2Cm6G6jY/ul1dz9YIllJ3TnWldMnE1sSZbRMQpxpgJwNPUToL+zVr7yFHH2wKvAR0B\nN/CEtfblesddwCIgz1p7VbhyqtCKNNOI83P49P01R4y53Ybh53ZzKFHL8O2+5/DUymlHXDQixRPP\npH7nntTz9bh+7KFi2xLllRfy8vp55JYXApDsieeq7mcwtlPfE3qeLimGK6fOZ94v15HecQ1fm94m\nHHFFRE5aXRl9DhgH7AIWGmPet9auq3e324Fl1tpLjTHtgPXGmNestf664z8C1gBp4cyqQivSTBNv\nGsz+/DIWz9+BtZCcGsfNt4ygXfsUp6OF1f69ZZSWVNM1JwOPp+EMYoekNB4860oW7tvOnsoSshJS\nOTur+wnvuRsNKv0+nl75OaW+wyevlfureX3zItrEJXJmu4YnCDYl7/l/s+CX6wgGoCAvwNIL36T9\nwitJyM4KR/QWY/bnm8numEqffu2djiIix3c2sNFaux3AGPMGcDVQv9DuAc6o+zoVOHCwzBpjugCX\nAb8HfhrOoCq0Is0UH+/hjrvHsC+/jOLCSrr1zCQu7sQ/ao4WVZU+XnhqNksX5mEtpGckMum2cxh8\nVsMN5+PcHs7r0MuBlJE1f+/WI8psfVN3rmt2oS1en8usO98lGDg8VpDn54ubH+HiTx8PRdQWada0\nTfz9j/OIj/dw16/H0btfbJd3kRjQGcitdzuP2pJb31+AacaYXUAKcEO9Y08BdwFh/whKC7ZETlBW\ndgq9+2XFdJkFePOVJSxZkHdo7XBRYSXPPTqT0pLGC11rkF9ZclLHjtbmtK4M+uXXjxhLyEpj+JO3\nn3S2lm7RvB38/Y/zMKfVrh9+7DfT2Jlb5HQskVZrVl4hv5+/5dCfU/BLYLm1thNwJvBHY0yKMeZy\nIN9auwwwdX/CRjO0ItKoebO2NhirqQmweF4uYy/p40Ai57VPTD2pY40584GbAVj223+QkJXGpZ8/\nFdOXNu43MJus0clkjC6huMxFu1Xt6NAprEvqROQYRnfJYHSXjEO3H1rQ8D0f2AnUP1GkS91YfedR\nu6QAa+1mY8xWoF/d+FXGmMuARCDVGPMPa+23QvZN1KNCKyLSTCPb9+C/21dR7m+49+64Tv1O+PnO\nfOBmPCmJdL3s7BZbZotrKpmSt47VhbuIc3sYntWdsR374DnBnR0qPT7ajy3jipmLMP3b8c7IRKzR\nlR5FWriFQG9jTHdgN/A14Maj7rMWuAiYY4zJBvoCW6y19wD3ABhjxgB3hqvMgsOFtm6x8D+AbCAI\n/MVa+4wxJgN4E+gObAOut9YWOxZUpBU6Z0wPpn2y4YixuDg3w0Y2/8SnWJPoiePHZ1zAS+vnsqui\n9i0p0e3liu5nMCzr5Ha7OONn14cyYkgVVVfw8PLJFFZXHBrbVnqA1QW7uOP0C3CdwLZqmfGWK6bO\nZ8Ev1xKXaLjpvTjcLXhbNhEBa23AGHM7MJnD23atNcbcUnvYvgg8BLxkjFlO7bKCn1trCyKd1dim\nNteMxIsb0wHoYK1dZoxJARZTe/bcd6g9S+5RY8wvgAxr7d2NPN46mV8kllVV+njxD3NYMj+39qSw\nzLqTwoY1PCksFgSDlhVLdrJp3T7ad0hlxPndiU9oeqeGHWUFVPl95KS2Jc4dmx92vbl5MdN3rW/0\n2A8GjGZw2y7Nfq6tf/wbs378r0MnwiWkGCZMfZSMs4eGIqpISBljsNY6+huXMcaW3TEuYq+X8uw0\nx7/nU+FooT2aMeY9avc7ew4YY63Nryu9M6y1DT7PU6EVCb8D+8opLamia04G7hjd+D8YCPKHh2aw\nbNHhpWHtO6Tyq/+7hPTMJAeTOeu+Rf9lb2Vpo8dGdejNN/ocfbJz46r2FfFWjxvxV9QcMZ55Rg5X\nLftri76AhrROKrTRp8X8dDLG5ABDgHlAtrU2H8BauwfQhoUiDmmblUxOr7YxW2YBFs7dcUSZBdi7\np5T//meVQ4laBo9p+v9zr6v5/x4SstIZ9Y97Dl3WuHasDWNev09lVkRCokX8hKpbbvAf4EfW2jLg\n6GlXTcOKSNisX53f6Pi6VY2PtxbDs7o3eeysYxxrTM7EUYx9436Mx0VCVhsu/fzJFnsi3KqCXfx1\n3RyeXTWDyXlrqfDXHP9BIuIoxxd+GWM81JbZV62179cN5xtjsustOdjb1OMfeOCBQ1+PHTuWsWPH\nhjGtyImx1moGKgq0zUo+ofHW4qLO/VhVuJvNJfuOGB/X+TR6pZ34RRFyJo7igrceoE2fzi22zL61\nZQlTdx6+CNKqwl18sXsjPxt8MW3iEh1MJuE0Y8YMZsyY4XSMBpLunBC5F3t2WuReKwwcX0NrjPkH\nsN9a+9N6Y48ABdbaR3RSmEQbay3l/uVU+FcTsGW4TRrJnjNI9p7udDRpQnFRJb/64X8pLTm8HZdx\nGe769TgGDu7oYLLwqq72g7XHPPktEAyyeP+O2m27XB6Gt8+hb5vYXAWWW1bI75Z+0uixsR37cGPv\n4RFOJE5pKWtog9sei9jruXLucvx7PhVO73JwHjALWEntsgJL7Z5lC4B/A12B7dRu29XgkjIqtNIS\nFdd8QYV/TYPxZM8Q0uJGOJBImmNXXjHvvr780C4Hl183kEFDY3NHh/Kyal7583wWzd1B0MKZZ3Xm\n5ltHkp7RumcgP9i2go9yG183neZN4LGREyOcSJyiQht9HF1yYK2dAzS1O/dFkcwiEgr+YCkV/rWN\nHiv3ryTFOwiXad2loaXq1KUNt9012ukYEfHCU3NYvvjwSXBLFuRRUjKT+x6O4MebLVDwGKdrWJ3K\nIdKitYiTwkRiRU1wF02fwxigJrAnknFEGig4UMGKJUdfuRI2rdtH3o4GH4S1KkOOsa/uiey5KyKR\np0IrEkLmOB96GNP0WkWRSKiq9NHUSq2qCl9kw7QwOaltOb9DrwbjGfFJXN5Na+BFWjLHdzkQiSXx\n7m4Y4rA03ObHZZKIc8XuCUYSHTp2TqNj5zR27yw5Yjw9M5Eefdo6lKrl+EbvsxmQ3pG5e7dQ6ffR\nLz2bsR37khqX4HQ0ETkGFVqREHIZL23iRlNUMx0I1jvipk3cGIxpasl4dKvwr6PctxK/LcRtUkjy\nDCDZM1hblrVAxhhu+cn5PPW76RQXVQGQnBrHrT8dFdMXz5g1bRNfzthCIGAZeX4OF0zoi8vV8N+n\nMYZhWd0YltXNgZQicrIc37brVGiXA3FK0AZxHeMqSr5gIRX+NQSCJXhc6SR5BuBxtYlgwsgp8y2h\n1LewwXiiux/p8WMcSCTN4fMFWL1sN8GgZeCQjsTHx+78xluvLuXDt4/cveCC8X349q0jHUokLZ12\nOYg+sfsOJhJiARvk09zVzNy9ieKaSjonpXNp1wEMb5/T4L5eVwZt4s6LfMgIC9oaynzLGj1WGVhH\nSnBIzBb5aOf1uhkyPPZPdKoor2HyfxvuPDJzyiauvmEQGZlJDqQSkVCL3c+XRELs1Y0L+GD7Sopr\nKgHYWVHEX9d/yZf5WxxO5hxfMB9L0ycSVQfyIphGpKGC/eXU1AQajAeDlr27Sx1IJCLhoEIr0gz7\nKkuZ10Rx/XD7SoKtdOnL8Xd10IdA4qz2HdNITo1rMB4X56ZL9wwHEolIOKjQijTD5pL9Te4ue6C6\nnKKaiojmaSm8rmzcJqXRYwYPCe6cyAaKEGsDVAV2UOnfQsBWOh1HjiEuzs313xzK0ecnXnvjYJJT\nGhZdEYlOmj4RaYYUb3yTx9zGRaK7df5gNMZFm7jRFFRPBvz1j5AWdx4u0/TfW7SqCmynuHomQQ4W\nWRfJnjNI9Y7Qrg4t1NhL+tA1J50vZ2wlGLSMGJVDv4HZTscSkRBSoRVphgEZHciMT6KguuFM7NB2\nXUn0xO4FE6r8PioDPtrEJeJqpLDFu7uSlfAVKvxr8AcLcLtSSfIMwOtq50Da8PIHiymsngLUX5MZ\npNy/HLdJIdmrzfdbql59s+jVN8vpGCISJiq0Is3gMi5u6T+K51bPpNRXdWi8e0omX+s1zMFk4VPm\nq+bNzYtYvD+XgA3SLiGFK7qdzjnZPRvc1+NqQ1rcOQ6kjKwK/xqOLLOHlftXqtCKiDhEhVYEqKn2\ns3LpLoJByxlndiIhseGMa05qW/5v+FUsPZBLYXUF3VIy6Z/eISY/ZrbW8syqz9leVnBobH9VGS9v\nmIfH5WZ4VncH0znHHyxq8ljAlmBtEHOM/YlFRCQ8VGil1du8YT9P/X46pcXVACQlebn9F2MYOLjh\nZWrj3B5GtO8R6YgRt7ZozxFltr5Pc1e32kLrcbWhOtj4MbdJUZkVEXGI3n2lVbPW8uLTcw6VWYCK\nCh8vPj0Hv7+J5tIK7GiizALklRe12m3Kkjz9aeptM8kTm8sNagK7OVD1Ibsr/sKeipcprplD0FYf\n/4EiIhGkQiut2p5dJezZVdJgvKiwkm2bDjiQqGVIj2/66klp3oRGTw5rDTyuDNLjxmGov3uDIckz\nkGTPIMdyhUt1YBcHqj+kJrgTCGKppsK/igNV/8Va/3EfLyISKVpyIK1aQqIXY6CxCcfE5NjdueB4\nhrbtyn+8Syj1NZyJG9OxjwOJwi9oLasLd5FbVkhGfBLD2nUjzt3wLTLR05MEdzeqA7kE8RHv6oTb\n1fhevNGu1LcQaPhJhd8eoDKwiSRPv8iHEhFphAqttGoZmUkMPqsLyxYeeYnWvv3b07lrukOpnBfn\n9nDbwDG8sGY2hXUXjTDAiPY9uLTbQGfDhUFJTRXPrPqc3PLCQ2Nvb13KbQPH0CO14fZjxnhI8MT2\nWmpr/fiCe5o8Xh3IU6EVkRbD2CheC2eMsdGcX1qGivIaXn1xAfPnbMcGLUNHdOXmW84mLT3R6Whh\ns2n9PhbP3UFcvIdzx/Yku2Nqo/cL2CBrC/dQ5qumZ1o72ic2fr9o9+La2Szev6PBeGZ8Er8ffhWu\nVniyl7UB9lT+ncZmaAES3aeRHj82oplEIsUYg7XW0bVVxhgb3PZYxF7PlXOX49/zqVChFanj8wWw\ntvZSmbHsg7dW8vY/lx267fG4uOPuMQw5q4uDqZxT4a/hznlvN3mi249Ov4ABGQ13vGgNCqunUhXY\n3OixzPjLiXfH3r+ZoK2mzLecqsBmLEHiXV1I8Z6Jx5XmdDSJIBXa6NP6ph1EmuD1ukNaZm31bmzZ\nKmxVbsie81QVFVTw3psrjhjz+4P862+LaK2/HFb6fcfctaGskXXErUWadwRu03B9cKL7tBgtszUc\nqPqAcv9SAraEoC2jMrCOA1Xv4g8WOx1PRI5Ba2hFQswGymHf+1C96/CYNwvaX4PxtHEwGWzZeIBA\nI9uR5e8upbiwkvTMpnc3iFUZ8Um0jU/mQHV5g2MGQ8+02LuE70FVfh8YSHA3fgKk25VKu4TrqPCv\noyawC2PiSHT3JsGTE9mgEVLhX4ffNtyyLkgVZb4lpMdf4EAqEWkOFVqRUNv/8RFlFgDfvtqS2/Fb\nzmSq0y678bPxk5K8JKXEN3os1rmM4YruZ/DKhnkNjp3XoSftEmJvB4PcskLe2rKE9cX5GKB/ege+\n2nMonZIbngjpMgmkeIeAd0jkg0ZYdWB7k8eqjnFMRJynJQciIWR9hVDVxA++mr3Yo4tuhHXLyWDw\nsM4Nxi+5qn/Mrx0+lnOze/K9fufRNTkDQ+2s7TU5g/l67+FORwu5A1XlPLlyKuuL8wGwwJqiPTy+\nYhpF1RXOhnNc0z8SDVG7tFCkVdAMrUgo+RtepOHI48UQ3ykyWZpw+89H8+E7q+t2OXAz+qLeXDC+\nr6OZwiUYtHzy3mqmf7qBinIfQ4Z34Yabh5Ke0XAHi7OyunNWK7ik7+e71lPh9zUYL/dXM3P3Rq7O\nGexAqpYhwd2DmmBe48difJs2kWinQisSSt5MandsbeIkI6/z6zHj4j1MvHEwE2+M/eLy/psrjjgJ\n7ssZW9ixtYAHn7oCl6t1zrhtLW36CnjbjnHJ49YgydOXysDGBvvvuk0KKZ6hDqUSkebQkgOREDKe\nVEhuYrP5hBxMXFZkA0WQrdqBLVmCrdiEtY3vXRpJwaBlykfrGoznbS9izYrdDiRqGdrENb2/8rGO\ntQbGeGgbfzlp3nPxurLxmHakeIbSNmFizF4NTiRWaIZWJNQyLwEMlK+jdlN6A4m9oe14h4OFhw2U\nw953oaberJY7Fdt+oqMF3u8LUF5W0+ix4sKqCKdpOUZ16N3oRSRqj/WKcJrIKPdVM23nelYU7MRl\nDEPbdePCTn0bvbSxMR6SvWeQ7D3DgaQicrJUaEVOkC9YQMCW4jHpeFwNt+EyLi+0uwybMQZ8heBp\nUztzG6sOfHZkmQUIlMLed7Gdv4tx6CpbcfEe+vTPYuPafUeMu92G/oM6OJKpJeif0YGJOUN4f/sK\nAnUz6R7j4roeZ9IrLfY+QSj31fDo8insqTy8vn17WQHLD+Tx00Hj8Lpa78mQIrFEhVakmQLBcopq\nplETPPxxdby7O+lxF+AyDbe8Mu5kcCdHMmLEWX8pVG5p/GCgBKq2QWLPiGaq75vfO5tHfz2VstLa\niyMYA1/95lAy28bmfrvBQJBli3ayZ1cJPfu0o9/p2Y3eb3zXAZyT3YOVBbswxjAosxMp3oQIp42M\nGbs3HFFmD9pSup+F+7ZzbrZz/z5FJHRUaEWaqaD6U/x2/xFj1YHtFNV8Tmb8BIdSOSxQduzj/tLI\n5GhC956ZPP7CNSyYs71ul4POdOzs7MUtwqWivIZHfz2VrZsOn/Q1bGRXbr9rNC53w1nytLhEzovR\nJQb1rSzYecxjKrQisUGFVqQZqgO7GpTZw8e24w8WN7r8IOZ5M8B4wPobPx7XPrJ5GpGYFMeYi/s4\nHSPsPnlvzRFlFmDxvFwWfLmdkaNa75ZTx1pS4DFabiASK7TLgUgz+INFxzleGKEkLYtxJUDKoMYP\nxnfFxHeMbKAIsDaALV2K3fMv7K5XsIUzsP7jzFRHwOrlje/csGpZ693RAeCsdk3vLXxWVrcIJhGR\ncNIMrUgzHG/LHrcrhk/6Op6MsYCBshVgfYALkvpA5sUOBws9a4O1OzpUbTs86NsH5WuxHW7CeJyb\npW+T3vga2MYuItGanNehJ8sK8lhTeGSxH9m+B4MyG141T0SikwqtSDPEu7riNmkEbMOTS7yuDnhd\nbR1IFX5BW0NVYCtBW4nX1Z54d8OrnBnjgswLsOnn1l4JzZ1ce0JcLKrYeGSZPShQDkVzoZ1za6kv\nvqIfSxfmYetd0yM+wcOYi3s7likSAoEgfl+A+ARvo8c9Lje3DxzDsv15rCjYibtu266BGR0xpnVe\nXEMkFhlrm7iiURQwxthozi/RxRcspLD60yNKrce0JTP+Utyu2CtwVYEdFFVPw3J4L1evqwOZ8RMa\n3dWhNbD7P4byNY0fdCVgut4e2UBHWTI/l/f/vYI9O2t3ObjuG0PofVrsbcUFtTs6vP2v5Uz/ZD0V\nFfifv8MAACAASURBVD76DmjPN//3bLrlZDgdTWKAMQZrraO/8RhjbHDbYxF7PVfOXY5/z6dChVYE\nKKguZ8n+XHzBAKdndKJrSuM/FK0NUh3MIxAsxeNKJ94dmx9ZBmwl+yr/haXhyV6J7j6kx1/oQCrn\nHbvQJmK63hbZQK3Yu28s5703VhwxltYmgcdfuKbJ2VqR5lKhjT5aciCt3vSd63lryxKC1P5y9N62\n5Zyb3ZNv9hmB66iPJI1xkeDuBjF+cnSlf2OjZRagMrCZNHte65ylTerTdKFNis2dFKy1ULUdqv+f\nvfMOj+o68//n3HunqTeEAFFF79XGxgVsg3HBNa5xuh3vpnlTnKz3l930ONnsZje9bRLHSRzHLe69\nYJtiGxsw3XQQQkICdWk05d7z+2MEkpg7o0GaOyONzud55nnQee+d88IjpHfOfd/v9whoXsieNiBa\nSl57YU/UWnNTB+++VcmSpUqKS6EYaqiCVjGkqWxt4KH973H6Of+6Y/upyCvhvLLM7j+MhSXjTe1b\nWNI/NAtaX0Xk5d/Xc13PhfzF6cnJQaQVhNrHIsXsSRreQBavRORMT19igL/N3tbY3x5KcSYKhWIg\noGS7FEOa9bX7o4rZk6w7diCluQwkDK0oZkzgRhfxVR8yFSE0GHY1FF0K3jHgHhEpZEfchjDy0p1e\n8ml8s2cxC4AFJ55HhqMHJFPJ3EXlUWtCE8xZkJltQAqFIj7qhFYxpGkL2Z/yRGKBFGYysPDpE2kR\n72LJtqhYljEDITLvR0fArKIttIWQrEcX2WQZ08gypkRdJ4QGubMirwxGSgmt22NErUjrRRpPpW/5\n5EKOHG7kaGUTALqhcesnFzJs+ND8sKVQDHUy77eSQnEGTMwbxlu19iexE/Mzczo8EYQwKPZcSWPw\nNUJWbeeqQbYxjVzXwrTm5gT+8B4ag69B53m9JVtpCh4jZB0n370kvcmlCxkGGfsDH2Z76nKxoag4\ni+/9ZBU7tlTT3NTB9NkjhrzmrkIxlFEFrWJIc3bpOF6p2kW1v+fjU5/uYsWoaWnKylk6zBBvVu9l\na/1RdE1jYckYFg8fjy56diAZWgEl3msJWQ2dOrRFaMJevH8wI6VJc+gtsGk+aQ9vI9uYOSRtjYXm\nQrpKIGRv+cwAcIHTNMHMudHayP1FhlvAvxekhKyKtBpmKBSKxFAFrWJI49YNvjT7Ep48tIUNdYcI\nWyYzi0Zy1djZDM/qX0+klBL8beD1IeL4yacSfzjEf295mcq2LqveHQ3VbDxeyWdnXIAmotvqXVoh\nkLnaniGrDkvGPm3sMA+Ro8Ww98108s+B409Fr7uKIGty6vNJAbJxPTSt49QHnIbXkHmLEIUXpDUv\nhUIRH1XQKjKWqspGnn9iJzVVzYyfVMzKq6dTVJwVdV2e28ttk87itklnJW3v0JqnCL38ILL+GGTl\n4jpvFa7LPobQ01vYvnZ0d49i9iTbGo6y6fgRFgxJb/v4souil3gmI7KnILGgcR2EG+iyNb4IIQbG\nh7RkIv0HoGnt6avQ/A7SPRyRHd1TrVAoBgaqoFVkJFWVjXz77ufo6Ihoqe7eWcuGdYf49o+vIDfP\n2cfmoTeeIPjIz7oW2lsIvfgAsrUJz81fdHTv3nj/RGXM2OYTQ7OgdWnD0ERODKkygUcfl+qUUkLI\nqqPDPAwIfPp4DM3+FF5kT0NmTQWrHYQLoblTm2gqadkSO9a6BVRBq1AMWJRslyIjefYfO04Vsyep\nP97O6y/tdXRfaZmEXvqbbSz81vNYTTH6EVNEPO96bYgeRAqhke86F7sfhznGXAwtN/VJOYiUksbA\nqxzveIzW0Lu0hjZQ1/EQzcH1Me8RQiD07MwuZgHMOPrL8WIKhSLtqBNaRUZSXdVku370iP16spDN\nDchYRatlYlXtR8svcTSHeMwrHs2BlhMxY5lKrb+FgBlmZFY+uhZduHqN8RSLq2kLbyVsdcl2eY3x\nacjWWdrD2/Gb0S5bbeEtuLWyjPw7J4y7FILVMWLDU5uLQjGEEELMT+CykJRya6ygKmgVGcm4imL2\nfRBdWI6vKHZ0X5GVC24vBDts41pBeqXAlo6czMYTlRw8raidVzya2cXRQvWDnaNtjfxpz9un/r75\nbh9Xj53NkrKKqGvdeilu/eJUp5hy2sO74saGdEGbtwDatkcky3qgR2IKhcIpXgc2EH+oYTwwLlZQ\nFbSKjOTya6bz7vrDNDX4T62Vjy3gvIujC5lkItwejLNWEF7zZFRMmzATbaSzxUIwEKalOUBhkQ9N\njz6J9OgGX551MetrD7C1vgpdRGS7FgwbgxanHSERZLADeaIGkVeEyE6/a5Y/HOJ/tr5Kc6jrw0VT\n0M/9e94mz+1lVtHQdJSKp+hgxokNBYSrCFl6PdS/CqG6yKKrGAqXIdQJrULhJBuklBfFu0AI8Wrc\nuJSxjD8HPkIIOZjzVzhLU6Of117YQ01VE+MnlnDB8on4fC7H95XBDgJ//iHmljURHUtAGzsVz6e+\niVbgXLvBPx58nxee3Im/PURRSRa3fnIhi84d69h+J5FSEnrufkKvPxaRKdMNjPlLcd/wBYQ3WlUi\nVbx+dA8P7NtgG5uSP5wvzc7801g76jueI2Adto359KkUeC5McUapwZJBglYNAh23NiLi+BYHGWoE\nJMKVuZJ1itgIIZBSpnWyQAghrYM/Stl+2ri70/537g/qhFaRseQX+LjmptTrhwq3F++nvoFVewSr\naj+ieDj6GGeno994eS+PP9g1oV1/vJ1f/febjCjPp3xMgaN7h15+kNDzf+5aMMOEN7yM7GjHe8e3\nHd07HkfbY/dLV8eJZTrZrjkEApVEG0noZLsy0863NbSZ1tBGJCEANJFFvut8vMa4mPcIl7P/bxSK\nwYIQYiXwv0QmZ38vpfzhafGvAB8m8kPFBUwDSqSUjUKIe4DbABPYCnxCyngWhCCEmE2kteBUjSql\nfKy3PJXKgUJxBsjW7cjqPyMP/wxZ/Rdk286Y12ql5RjzLnC8mAV445Vo9QbTlKxbvd/RfaVpEl5t\n/3PG3LYeq/aIo/vHY5gvp0+xTMejj6TAfQm66FJvMEQhRZ7LcGlFaczMGdrDu2kJvX2qmIVI20VD\n8CVCVrQms0Kh6EJEHmX8HLgUmAHcIoSY2v0aKeV/SSnnSSnnA/cAqzuL2bHAHcA8KeVsIgXqzb3s\n9wfgD8D1wKrO15WJ5KpOaBWKBJFNb0Hjmq6FYA0cfwYZbkbkn52+xIBQ0LRdD8ZYTxrtzciWGEWB\nlFjHDqOVpmfYbHHpeJ45vJX2cCgqtmxkZuqJ7muu4/XqPTQE2inPLmTZyMmU+qJlx3zGBLz6eMKy\nHtA63eAyk7ZQrKFoi/bwNvLd56c0H4VikHEWsEdKeQhACPEgcDUQa7r0FuCkdmUzEASyhRAWkAUc\n7WW/xVLK6X1JVJ3QKhQJIK0OaHrLPtj0FtIKpDah04jVK7voHIeNErJyIc4AmDYsfYNXOS4Pn5+x\njGHertNYj25w3bi5LBrmfG9xqnmjei8/ev8l3q49yO6mWl49+gHf3fgc+5rrbK8XQuDSijO6mAUI\ny8bYMSt2TKFQADAK6O7Ic6RzLQohhA9YCTwKIKVsAP4bOAxUAY1Sypd72e8dIUSfClp1QqtQJEJH\nlY2UTycyBIEq8E1IbU7dWHnVNA7uq2fDukMAGIbG1TfNZsoMZyezhW7gOv/qnj20nehTF6CVOVM4\nSinZ8t5Rtm+ppqDIx3nLKsjLj3aAm5BXwncWrmJ/y3ECZpgJuSV4jf4PBpqHdhF+71UIdqBPXYg+\newlCS58VrD8c4pEDG6O6YgNWmIf2vcc981amJa+BgCFyYha13dsuFIqhxur1+1j91r5kvuUqYI2U\nkf9wQogJwBeBsUAT8IgQ4lYp5QNx3uM+4C0hRDUQICLjJTtbFuKiClqFIhFEL/9VhPPqCfEwXDqf\n++oFVFc1UVvdyriJReQX+FKyt2vlbciOdsJrn4ZQAISGPnMxng/f7diev/mfNax/4+Cpr595dDv3\nfG+F7QCcEIKKvOTp/wafu5/Qc/ef+jq87lm0yfPw3vk9hCs9Tlo7G2sImPYfuA621tMQaKfQkz7F\niXSSZcygObQ2RqxPB0EKRUaw9JwKlp7TJWX57Z+8ZHdZFdD9UV9555odN9PVbgCwEFgrpawHEEI8\nBpwLxCto/4/IENlWwIr/N+iJKmgVikTwloOeDWZbdEzPBc/A0DQdMSqfEaPyk/qe0mwHyw9GPsKm\nsBeajue6f8a98jasY5WIghK0wtKk5tCdXduO9ShmAVpbAjzyl038y78tc2xfAOvo/h7F7Kn13ZsI\nv/kkros+5Oj+sehNZ2fQ6vAkgSxjBmHZSHt4B13KDgb57nNx6859nyoUGcIGYGLngFc1kaL1ltMv\nEkLkAxcSUTs4yQfAvwshvEROWy/ufL941Ekpo4XcE0AVtApFAgihI4tXQt0TPVsPhAHFl/aqaTkY\nkWYb1L8M7XsBCZoPmbcw5gCcyMpFH+/8idfunbVntJ5MwhtXx45tWp22gnZaQRle3aDD5pR2fG4x\nBRl6OhuyTLY3VNMRDjG5oJQiT3bUNUII8t3nkW3MIWhVAgZefSya8KQ+YYVikCGlNIUQnwNepEu2\na6cQ4s5IWP6289JrgBeklP5u974vhLgfeI+IbNcm4LfEZ5MQ4gHgKSJF8Mn36lW2SxW0CkUnUlqA\nZXsKCSB845EjPg6tWyDcAEYh5MzOSL1KKSUce6TLLQkip7SNbyLREPmL0pbbsFJ7ya2SYSmQ4rJR\nTEgo5jBew8WNExbw5z3vILt10np1FzdVZKZl686GGv5v11paw5HfeRqCpSMnc+OE+Qgb1ztDy8XQ\nVIuBQnGmSCmfB6actvab077+E/Anm3t/BJyJO4SPSCG7ovvbAKqgVSh6w7RaaQ69TYe5H7BwaWXk\nuhbi0aPbCISrAAovSH2Sqca/v2cx252Wd5F58xEiPUNQC84ZQ9nf86g52txjfdWHZjq+tz5jMaFX\nH44Rc1a6zbIkSGlraQywpKyCUdkFvFG9l/pAG6NzCrlwxCRKvP0v9K2GOqyDOxDZ+WgTZyO09D6R\naAl28KsdbxCwuk6kLSSvHv2Asqw8LhwxKY3ZKRSKviKl/MTpa0KIxYncqwpaxZDGkgFOBJ7ElC2n\n1kJWDfWBZynyXIFHH5nG7NJIMM7je7Mt8jJiy3U5idutc8/3VvDkQ1vY/n4NBUU+Lr1qGvPPGu34\n3vqkOejzl2FufK3Huigtx7Xsekf2bG8L8tffv8vbbx5ASjhryVhuu2MR2TnRj8zH5RYzLrc4aXtL\nyyL42C8Ir3kKrMh8higZieeT/4FePjFp+5wpb9Ue6FHMdueN6j2qoFUoBilCiD/S08ZQB5bRczDN\nFlXQKoY07eFdPYrZLixaQ+8N3YJWj3eqp4MWLZGVSgoKfXz0zvSYWXg+eg/haYsIv/sKhDrQpy7C\ndf5ViDh6vP3hN/+zhs3vdg0Vr3v9AM1NHdz9zUsc2a874dcfI/zGEz3W5PGjBH7zdXzf+DMiCRJo\nfaEh2B47FvDHjCkUigHP093+7AKWEl8V4RSqoFUMaYJmbNOSoFWNlNK2H683PnZNtC4rwJ8e/8gZ\nv1dayJ4CDatB2hhGZE9FaOmRp3ISafqhbXvkdFrPg5yZtv3RQtNwnb0C19krbN4luZyoa+tRzJ5k\n2+ZqjlW3MHyEszqqobVP267LpuOY29ZjzE1P+82YnNgWvWNyMtsoQqHIZKSUj5629KAQojdlBEA5\nhSmGOELELswErj4Vs5mA0NxQenX0SaynHIouSk9SDiKDdXD0D5Eivm0HNL8FR/+IbNud1rza2oIx\nY+1xYslCNh6PHWs64fj+sVhQMsbW0lcguLRcDX4pFJmCEGIqsCeRa9UJrWJI4zMm0mHujRnLVCzp\nJ2BWAQKPPhrNprAX3jHIUZ+OyHZZbeAuQ3id71NNCydeiKg49MCEE88jfePSdiJdPqaAkmHZHK/r\nqX9cWJzFmPHOn0Rqoydj7dsSI5a+PlWXpvOlWRfzt33vsvVEFRaSEb48rhk3h2mFZf1+fxkKYu3f\nCgi0illpa61QKIYaQogWunpoNSAL8J9cl1LG7O1SBa1iSOPVx+LTp+I3d/VYN0QRua70SVM5SUvo\nPVpDm4jIAkZOovPci21dk4TmhpzMPvGSoQYI1sQIBsG/D7KnpTapTjRNcMddS/jJ91+jvT0iC+b1\nubjjC+eix1A7SCauFbcS+PU2kD0Ne7RJc9EnOKcqUV3VxPvvVpGV42bRuWPx+aILykJPFp+ZfgH+\ncJCAGU6a1m5442sEHv45tDVFFnIL8Nx4F8ac85Py/gqFIjZSyj73UamCVpHR+MMhmoLtFHqy8ej2\n3+4FngvxmRX4w3uRhPDo5fj0iYg029k6gT+8l9bQuz3WJCGagm9iiELc+og0ZZZGZC+P7i3nH+3H\nY+rM4fz499ezecMRkDB30Sh8Wak5MTamLYTbv0Xo2fuwqvaBx4ex6BLcV93h2J7PPr6dh/60Edl5\nRvPQ/Rv56reWM2ac/Ym0z3DjM5Lz72Ee2Uvg/ntPKToA0NJI4L7voX31N2gjxiZlH4VC0RMhRJmU\nMsbJQmLXqIJWkZGELJOH9r3H+toDhCwTj25wQdkkrh0/B93G1cujl+PRy5Oyd+WzbyflfZygLbwt\nTmz70CxoXSWgZYEVY3Le26tajOP4fC7OuWB8Ut9TWkEIVAE6eEfF1BU2Zp2DMescZEc7uDwI3Tn9\n4RN1bTx8/6ZTxSxAS1OAv/5uA/d8z/khvPDap3sWsycxw4TWPYPn+s84noNCMUR5Fpjfn2tUQavI\nSP68523erj146uuAGealqp1Y0uJGh52Tttz7AEvXRheORo6PUOuHcOX4HN0/HqbVFDsmY8d6YzCr\nOgihIwvOjdj8nk72TIQr86bmZcsmaFjTpWKhZSGLlyOyYvfFCq/z9rnbt1RHDCROY9f2Y4RDJobL\nWTMP2RBbfzleTKFQ9Js5QojmOHEBxIsrlQNF5lEfaGND7SHb2Bs1e2kPO/sIuWHLftv1cKuflgPV\nju7dG7oW26ZXF5ln4ZsoInculFwJ7uGABkY+FFwAxc6fCqYa2b4P6l/pKclmtUPdUxG1hzSSX2D/\nYS8n14NuOP/rShtVESc2wfH9FYqhipRSl1LmxXnlSimj7Tu7oU5oFRlHdXsTFtGnPBBpRajztzI2\nN7aOZX/JrRhJ/eZo5QTd4yK7fJhj+yZCtjGTRtsBKEG2MSPl+TiNlGHawjvoMPdgyTAefSTZxhwM\nLXpQVmRPheypacgyxbRsjBGwoGUzFC9PaTrdmTV3BCPL8zl6pOfTguVXTk2JhJ5x3ipCa54Cf2vP\nQHYexrlXOL6/QqHoOzELWiHEdQnc3yGlfDaJ+SgU/abYkx0zpiEo8Dj7yH/GFz/Emx/7QdT6xI9d\niqfQWSH83vAZFZiyhZbQe0DEOlTgJs99Lm69/3JHAwkpTeoDzxG0uswz2sON+MP7KPZejUvLvFaC\nhAjVx46FG1KXhw2arvHVb1/C3+97j83vVpGd4+ailVO4/NrUKG1ohaV4P/cjgo/96pRcmTZpLp7r\nPoOWnzxLYYVCkXzindD+DniCSN9CLC4g0qSrUAwYyrLymVownF2Nx6Ji80tGk+92tqCd+JHlBJva\n2PL9v+KvqcfI9jLpk5ex6Ed3Orpvc9DP+toD1He0MzqngEXDxtkqO+S45pJlTCNgViHQ8OijMlLR\nocPc36OYPYkkQEvoHYo8l6YhqwGAqxBMO7tnwEh/kV9YlMU/fckZiSxptoHQEXGsm/XRk/Dd9WNk\nWzMIgchK74dQhUKRGPEK2ueklJ+Md7MQ4i9JzkehSAqfmrKE3+5cw57mrkGOmYUjuW3SWSnZf/rn\nrmHqnVfSfvQEnpI8XNnOFtE7G2r41Y43CFjhU2vPHN7Gl2ZdwjBfTtT1mvDgM5LTE7j79wPzM22H\neTBmLGAeGrq2xrnzoeOwTUCD3LkpTycVyPZ90LgGQnWAQPrGQ+GyuAN/IjumfrtCoXAQEZFcGU63\nGlVKafdDqwcxC1op5W293ZzINQpFsgmHTNrbguTme2MWJHluL1+ZcwmVrQ3UdbQyIiuPEVn5/d7b\n/GAjodcewTpWiVZajrH0Ooxp9gYMmssgZ+zwfu/Za06WxR8+WNejmAWoD7TzwN53uGuWs1a1O376\nGEu3Hoha11wGHSeuwlvc/393RfIQWRORhcsiBZ6MmDWg+aDoEoQ7vT3eTiA7DkPd43SZD0nw74dg\nLXLExxF67NNahUKRWoQQnwe+ARwDTmroSWB2b/f2OhTWWSlfAYyjZ7X84z7kqlD0iyf+voXnn9xJ\ne1uQspF5fPj2hcyeH3vwcXROIaNzkvMYNbzhFQJ/+QEnRTLNE9WYu95F3vJlXItXJmWPvrCzsYbm\nUEfsWLCDPLdzv7RbD0a3dgBYoTDtVcfTVtB69XF0mPaKEx59bEqGjFKNlCYd5kFM2YZLK8atjbT9\ne4q8BcicWdBRCUID72iEyNAZ4aZ3wG5I1GyFtq2Ql5mOgArFIOUuYIqU8sSZ3piIDspTwMeBYiC3\n20uhSCmvvbCbx/72Pu1tEdmtmqPN/PTe1dQdi9EPmESkaRJ88nf0UHwHkJLQ039AmmH7G1NAIM7e\nspd4MiiaN9F23ZWXTW7FSEf3jodXn4Bbi95f4CHXlZrWk1QSNI9R6/8rjcGXaQmtpz7wNCcCj2FK\ne8MIobkRWRUI3/jMLWYBAtF91F2x9MroKRSKKCqBPomiJ/JTrFxK2etRr0LhNK+/HC2FFQpZrFt9\ngKtvcvZbVB47jGw6bh9rrseq2o8+ZrKjOcRiUn4phtAIy2iHo1JfLiXe2KoPyWDu12/jxcv+FWn2\n3H/ml29wvHe4PRzErekYWrTgvhA6RZ7LaAtv72ZrPCqmbNdgRsowDYHnseh5Uh+yjtMUWE2R9/I0\nZTYA0LMglva07rxZhEKh6B0hxJc6/7gfWC2EeAY4JZadSFdAIgXtC0KIFVLKF/uWpkKRHAId9ieN\nHTHWk0ovUl+pcFGKRZ7by/LyaTxXub3HukBw7bg5jj9aH3nJApY/ey/vf/cvHH93NzljhzP989cy\n9Z+vcmzP9+oO8/ThrRxtb8Kt6SwuHc914+fhM3qqNQhhkOOaQ45rjmO5nMQKhdFc6Tnp7DAPRBWz\nJwlYlYStFgxtiD5Yy5kFjW/ax7JnpjYXhUIRi5M/oA53vtydL7DtGYomkZ++64HHReS3YoiIjJeU\nUmbWEYdiwLNg8WiOVkY/iViweLTje2vFZWgTZmLtj7a01cZORSstdzyHeFwzbg7Dfbmsrt5DfUcb\n5dkFrCifzrTC/mvLyvZWZFsTomg4wkYGDGDU8oWMWr6w33slwqbjlfxu15pTP+GClskbNXs52t7E\n3XOcNQVoPXSMZU/fj7SiT8P3/XUYkz6enl5qU7bGjVuylb52ig16VYe8hRCoAf+ebosCCi9EeDJL\ne1mhGKxIKb8FIIS4QUr5cPeYEOKGRN4jkYL2v4HFwFYpT28gVChSx6oPzeLg3nq2bor0xOm64Oqb\nZjNxSmomsz23foWOX9yNbOiyBxX5xXg+fLdje4ZDJm+vPcT+3ccZPiKXJcsqyM5x2157zvAJnDM8\nefacsr2VwMM/xdz8BphhRF4RruW34Lrw2qTt0Reerdxm+3F9b3MdHzQeY0qBc8oSlc+8ZVvMAlQ+\n/VbaClpDi+d8p8W1PB7sWDIEhNGE/VMUIXQovRoZOAr+gyAMyJ6CMJT6hkIxALkHeDiBtSgSKWgr\ngW2qmFWkG4/H4CvfuJhD++upO9ZKxZQSCotS96hfKy3H9//uI7zpdeSxw4jhozHmXYhwSEEgFDL5\nz2+8zO4dXVq6zz+5k6/feylFJc72xQJ0/OFbWLs3nfpaNtcTfPQX4PbiOucyx/e3w5IWh1tju1kd\nbDnhaEGre+0/TEDE2jhdeLQxGKKAsGyMivn0yegxir3BTNhqoTm0joB5CJAYoohc1yK8xjjb64Vn\nJHjSN6SoUChiI4S4DLgcGCWE+Gm3UB4nbS17IZGC9mSD7nOcYYOuQuEEYycUMXZCvBOpM0NaQWh+\nB9p2gRUE71jIPxvhLom6Vrg9uM5ekbS947H+jQM9ilmAE3VtPPPYdj7yaWen9M3DH/QoZrsTeuXv\naStoNaGRY3hoDQds407KkwGMvWYJb33+Z5j+6P3H37zM0b3jIYSgyHM5jcHXCFonJ/c1fPpE8t1L\n0paXU1gyQH3gyR6tFmFZT0PwRQrFSrz6mDRmp1Ao+sBR4F3gKuC9bustwBcTeYNECtoDna/uDbop\nQQixEvhfIvJiv5dS/jCV+ysyHylNqH2kp7RP+07w70OW3ZJWofm9u+ps1/fEWE8mVvWhmDFZewRp\nmgg9WlkgFZw3ooLnK3dErWcbbuaXOFvIeIryOP++r/HmR+/FDIROrU/73LWMWXWuY/vubz7OphOV\nAMwtLqciL/r7UtdyKfZeRdhqwJRtGFohuujfSX6s9op00x7+IEbfsKQ19J4qaBWKQYaU8n3gfSHE\nA1LKUK832NBrQXuyUTfVCCE04OfAxUQq9w1CiCeklLvSkY8iQ2nfba9TKYPQtB6GOTep3xvDhkdb\n1gKUltmvJxOtZETMmCgclrZiFuDKMbOoaW9m84kjp9ZyXR7unHY+nhhDa8lk/A0XUnbhbA4+/Dqh\nVj+jr1hM4czxju331z3v8EZNl2Tdi0d2cl5ZBbdNPMtWwcLQCjFIjpnIB797hqVP3he1njtxFFLe\nljZziq5T6GhCVi1SmpHe2T4w6IfgFIrBzUYhxOktrk1ETm+/G89wIRGnsIXA/wPG0tMpzGlt2rOA\nPVLKQ515PAhcDaiCVpE8/NGWrQnFUsCFl0zkpad30dTYJcdkGBqXXTPd8b31illo5ZOwjuyJ/9qR\nCAAAIABJREFUirkucHYozDQtwmELj8f+x5NL0/nn6RdQ2drAvuY6cl1eZhePwmWjRXsmyNYmQuuf\nxarah1Y4HOOcy2KqV/hKC5n22Wv6tV8ibD5e2aOYPcmamn3MKBzh+In03vvt1Rpb9lZRu247w5ek\nR/ZKwxMzJnCRmGeQQqEYgDwHmMADnV/fDGQBNcB9wKpYNyZynPFX4G5gK12+uqlgFJGBtJMcIVLk\nKhTJI55DUprdk/IKfHz9Byt56uGt7N19nLIRuVxx3UwqJqdI1eHT3yFw//ex9m6JLLg8uC64BuOi\nhBRUzphw2OLhP29i9Yt76PCHmDpzOB+98yxGjbaf0E+mrbFVcwj/z74MLZGhKhMIrX4Mzyf/HWOW\nc60EvfF27cG4MacL2lCLP07M3oEsFfiMyfjND2LEJmWkrbFCMUS4REo5v9vXW4UQG6WU84UQt8W7\nMZHf2MellE/2Lz+FYoCSNQVat9jHsqekNhcbSsty+dTnk19QSSsEgSOABt5y28ezWkEJvi/8GOtY\nJbK5Hm3UBESWc+L8jz2wmeef6OqN3bXtGD/65iv86NfX4HI52+IQePSXp4rZU5ghgg/+D/q0RQgj\nPQoGfjN2K5k/3Kc2szOi/NJFNG6LflLhys1K2+ksgEcfSbYxh7bw+z3WDVGckbbGCsUQQhdCnCWl\nfAdACLEIOPkLIK7aQSIF7TeFEL8HXqanysFjfUw2UaqA7scP5Z1rPZP75jdP/Xnp0qUsXbrU4bQU\nmYTwjUXmzI4uao0iyE/fyZyTyJbN0PAmyM7/zloWsuhiRIwCXhs+GoY7b17x2gvR7Q0NJ9rZvOEI\ni84d69i+0t+KtXujfaylAWvfVvQp823jTjO1oIydjTW2sWSYZvTGzLtv5NA/3qRlf8+e1QX33o4r\n11nJvLBlogkNLcZpa557MV59PH5zL1JGbI29+oQ+984qhjarV69m9erV6U4jipphQ87h73bgD0KI\nHCJGXs3A7UKIbODeeDcmUtB+HJjSee3JlgMJOF3QbgAmCiHGAtVE+ihuOf2i7gWtQtEXRPEKZNZE\naNsJVgi8YyBnJkJLqahHSpD+A1D/cs9Fqx2OP4N0FSLcpWnJy7IkHX77E8f2tqCzm0sZecUKW6az\n+8fhghETWVuzl9qOnhP9w7w5XFA2yfH9faWFrHrnl3zwm6epeXML3tICptx+BcPPm+XYnrsbj/HE\noS3sba7DremcVTqO68bNJdsV3Tfr1ofj1pOnO9y467CtE5zmMmg9vIKcMc5pHCvSy+kHYt/6Vlrm\n4Yc8UsoNwCwhRH7n193tQR+Kd28iBe1CKWXKn71KKU0hxOeAF+mS7dqZ6jwUg5f28C7aQtsxZTOG\nVkC2MQufMdH2WuGbAL7kuWwNWFrsTyLBgpZNUHxpStM5iaYJ5iwcxaZ3jvRY13XBnAWjHN1bZOWi\nVczC2rc1OpiVi17h9PxrbLIMN3fPWc5zldvZdPwIEsnc4tFcPmYG2a7+feCSbc2YuzeBbqBPXYhw\n2w9aeYrymH3Prcy+59Z+7ZcI+5rr+Mm21wjLSEEZtEzW1OzjUEs998y7FF04O+y1788v2UqVWaEw\n+x94ldn/GnWmolAokogQwgNcD4wDjJP98FLKb/d2byIF7TohxHQpZbTwo8NIKZ8ncjqsUJwRzcG3\naQtvPvV1yKqlMfgKpmwhxzUvjZmlmVBsl624sRRw2+2LqD7STM3RZiCi6PDRO8+iIAVucO5r/4mO\nn98NHd0GnYSG+7p/jlnoJYPDBxvYuaWGgiIf888ebdsrnOf2cVPFQm6qWJi0fUOvPUrw6T9AqLPt\nxJeD55YvYcy9IGl79IXnKrefKma7U9nWwObjR1gwzNkhuGBTW+xYo53urUKhSDJPEJHpeo9uba6J\nkEhBuxjYLIQ40PnmApApkO1SKPqEabXRFrYf9GoNbSTLmI4mnCtSBjRGAYSj7VEBcCVHMaCvlJTm\ncO/PVrF9Sw2tzQFmzCkjryA1lq36mCn4vvYbQm8+iXVkL1phKcZ5q9DHTnVsz7/98b0eQ3ClZbnc\n893ljtsamx9sJPiPX/Vc9LcS+NP30UZNRBuWPnvY/c0xJSbZ33Lc8YJ25PIF7PrlE/axFcn7QKFQ\nKGJSLqVc2ZcbEylo+/TGCkW6CFhVxFKYk4QJmtUx/d57w050fVAJrufNh46DNgEBOXNTnU0Umq4x\na17yCiopzcjAX/f+6LwFCCMveu/iEXiuuTNpe8dj3+7jPYpZgNqaFh7562Y+fZezVrWhtU/bB8ww\n4befx33lJx3dPx55Lg9tMWyNc13O2hoDjFl1DqNWLqLq+Q091sdefz4jLxrCT3YUitSxTggxS0pp\n0wMWn0ScwmJ7YCoUAxDRm6i6w314Axnhm4AsvAga10Tc0AA0HxRdjPBk1sCLlBbUPd7TICNUB207\nkWW3Ilz2+rapYNtmG3c6YNsm+/VkIuOcgsqm2LFUsKSsgkcObIpaN4TG2aXjHN9faBoXP/4d9t73\nAof+sQY0wfgPXUjFR5c7vrdCoQDgPODjfekKiFnQnhSyjXdzItcoFKnGo49GYCBtJOs0vHg0Z4eM\n0oWUFgGrEtNqw6UV4dbtZZ1E3nxkzkzoqIwU997RiDSbSDiCf5+925vVDk3roOTy1OfUSV6+/Wlj\nKlostNFTsPZvt4+Nmez4/vG4aNQUDrfW805d1zmKW9P5+ORzKPT0r5daWhby2GFwe9GKY0ue6W4X\nUz59JVM+fWW/9lMoFH3isr7eGO+32DQhRAzFeSBSNef3dWOFwik04SHPfS5NwTdOiwjy3EsyUqcy\nZJ2gIfA8puwaXHFppRR6VqKL6CJJaG7IqkhliqmnfV+cWLSdbCo5+7xx/ONv7/ewNQZYfoXzM7Cu\npdcSfucF8PccgBJFZRiL0nsSqQuNT01dwory6exqrMFnuJhXPKbfig7hTa8TfOK3yPpjAGjjpuO5\n9ctoZc7pGysUijNHSnlICHEeMElK+UchxDAgJ5F74xW0iUxDpE+gUaGIQ5YxDUMU0h7eQVg2Y4gC\nsl0zcGmpsY1NJVKa1Aeew5I9C5SQVUtT4DWKvOk7iRy4pNcaNSvbzb9+ZwUP/uk9dm6pobA4i0uv\nmsaFy53XltWKR+D9wo8JPfE7zA/eA01Hn3Me7qs/jfA6pygR6AjxynO72bm1hoKiLC65fApjJxTZ\nXptMW2Nz31YC930PuqknWAd30PHzr+L79/sQntQMHioUit4RQnwDWEhE4eqPgAv4C9DrcEHMglb1\nzioGIgeaj/NmzV6agn7G5BRx4YhJFMR4FOnWy2I+du8Lu37zlO36jp/9g+mfvzZp+5wpHeahqGL2\nJAGrkrDVjKFFD0FlPFmToG1b7FiaGTk6ny99/aKkvqeU4cjps9kCrlLwjkHYOG3poyrQP/MDZCgI\nmo7QnX1qYZoW//mNV9j7Qd2ptXWr93P3ty5h6gxne7dDqx/tUcyeRDafIPzuK7iWqNYChWIAcS0w\nD9gIIKU8KoRIyC4tAxvnFJnK60f38Ld9Gzjp6bStoZrXq/fy5dkXMyrb2QEfKSVbf/A3lh46FhXb\n+l4J0z57NUJLz7CZKVt6jRuceUFrp+gAg0jVwTchUri2n2apq+dAwTnpyclBZKAaav8R6RE+iXs4\nsvQ6hG4vBSb6+Sg/UTZtONKjmAUIhy2eeHALU7/jbJuDdawyTuywo3srFIozJiillEIICdBpeZsQ\nQ3fcWzGoaAsFePjARk43KG0LB3h4fyz3q+QRbvXTalPMArRXHSfQEL+odBJDi1fMCwyRvmn+dCKE\ngJJVULwSvGPBPQLyF8OIjyCMzGr/l9KMKDp0L2YBgsfgxEvpSaoblQftTTsqDztv5qENiz0Eqg0r\nd3x/hUJxRjwkhPgNUCCEuAN4GfhdIjf2ekIrhPg88BcpZXpthBRDmq0NRwlZ9i3buxpr8IeD+Azn\nTpuMbC++EcX4q6NljbzDCnDnJ9Sz7ggebQy6KMCU0YYJXr0CXXNWqD8dBMyjtIU2E7SOoQkvPmMK\nOcacqIE/ITTImRl5ZTL+fWDGcLny70OabTFPaVPB2PH2vbJjxjlv5uFaeh3mtvUgT/s4nFuAsfBi\nx/dXKBSJI6X8LyHEcqCZSB/tf0gpE/pUnsgJ7XBggxDiISHESmHXkKVQOIw8/ZdR9xhgxQ4nBaFp\nzPjih2xj0++6Ds1In3KCEIIiz+W4tNLuq3j1CeS702tl6gQd5iHqA08TsCqRBDFlM62hDdQHno/7\nfZLRhOPZssrYxW6KmLtwFJOnlfZYc7l1rrl5juN765Pm4rnta4i84lNrWvlEfJ/5T4Qv8z7sKRSD\nHSnlS1LKu6WUX0m0mIXEjBW+LoT4d2AF8Ang50KIh4DfSynj6OIoFMljZuFIDKHZ+rxPyivtt6xP\nIsz6yo3IsMm2Hz9M4HgTnuI8pt91PbPvudXRfZuCfqrbmyjyZFPqs++NN7RcSrzXErJOYMpWDFHo\n2CCYtKy09QsDtATfhqjmEwhaRwhaR/Doo1OfVLpxxxmsEu6I5XEa0XSNu795Matf2sOOLTUUFmVx\n8eVTKB/T/7yk2QEyhDBiz40Yiy5Bn78Mq2ofwuNFG+6sha5CoTgzhBAt2P1g7zJW6PUXWkJDYZ0N\nujVADRAGCoFHhBAvSSm/egY5KxR9Itft5epxc3j0NBchr+7iQxOSY0kpwyEww3FlfGb/6y3M/PIN\ndNQ14inJR3e7krK3HWHL5IG977K+dj9W58nj9MIRfHLyOeS67YX5XVoxLoptY8li5y+eSJuqg2m1\nEo7T/dRhHu5zQTuYh+CEdxTSUw6BI9HB3HkR3eE04/YYrLhyGiuunJaU95PhJqh/pdM8QyKNIihY\ngsi21/IVuo6eZuMIhUJhj5QyISWDeCTSQ3sX8FHgOPB/wN1SypAQQgP2AKqgVaSEFeXTGJ1deEq2\na2xOEctGTmGYr3/9q7K5gcDjv8bc9AaYIbTxM3BfdTt6xSzb6zWXQdbIkn7tmQiPHtjE2mM9H4Ls\naKjmt7vW8OXZlzi6d8P2gyx98j7b2J7DE9MnU9aLbbEg80wzEqb0GjjxMrTvBiwQHsibB/m9yjcO\nOqQVhJq/g9nctRiuh+NPITUXwjchfckpFIp+cZqxQgmQK6W0sX3sSSIntEXAdafr0kopLSGEEvBT\npJRphWVMK0yetqwMh/D//CvImq5vb+vAdjp++TW8X/wpevnEpO11JnSYIdbW7LeN7W6qpbK1IWnC\n83YE6pv7FHMaXWTh1soIWjW2ca+eme5nYauJDvMQIPHq4zC0aJUGoXlh2JVI0x/pmTXyEZpzTxDS\nStv2nsVsd5reiki2KRSKQYeNsYKbBI0Vem2Ek1J+I5bJgpRy55mlqlAMLMz31/QoZk8RChJ65aHU\nJ9RJU9BPwArHjNf6nZUJK1kwGXeB/cn3yIvnO7p3b+S5lyCIfoSeZczArWeeE1xz8C3qOh6kJbSe\nllDkz83BdTGvF7oP4S7J3GIWIGD/gabXmEKhGOhcC1wFtEHEWAFIqB1B6dAqhjRW5e44sT0xY05T\n6M7Co8d+gFKW5azzl5HlZcH3PhW17i0tYM7/+7Cje/eGSyuhxPshso3ZuLQyPPo4Ct0ryHefl9a8\nnMAf3kdb+P2o9bbwVvzh9H1/ph09jkVvvJhCoRjoBGVEruaMjRWUU5hiSCPyY/fCioL+9cne+eYD\np/78m/PPTAnBrRtcUDaRl6p2RcWmFZQ57owGMPWfryJvcjm7fv0k/up6hi2ezvS7riNndGnvN/cB\nKSXv1h3irdqD+M0QU/JLWTZyCnk2A3CGlkueO3luXzVvbEnaeyWT9nDsh2Dt4Z34jL5Z+A7mATgg\noivc/C62Q9GZrjmsUGQ2pxsrfJJkGSsoFJmMsegSgs/eBwF/VMy1ZFXqE+rGtePnEpYWa2r2EbJM\nBII5xaP46KSz+/3eVvUhwtvWI3QDfe75aEX2sk8jL56fshaDP+95m7XHuvqG9zXXsf7YAb46dzlF\nHmf1Qjd/+36Wvropat1Tko8ZuAndkx6VAFPG1o+NF8t0hKsYWXRJROWAblJ+vgkRNziFQjEo6Y+x\ngipoFUMakZOP947vELj/+8jm+sii4cK14sMY89JrSqALjZsrFrJqzCxq/M0UebIp9PT/cWrgsV8S\nXv1Y18ITv8V99adxXWRvHJEK9jcf71HMnqQh2M4zh7fxkSQU8fGoe8v+JDRwvInmvUcpnDHO0f1j\n4dKKMc1oBzgAQ7N33xoqiNw5SF8FtH8AVhC8YxDe2Da3CoVi4COE+BLw9zMxVDiJKmgVGY+UZpQl\nanf0yXPxfesBzN2bIOBHnzgHkRM9RZ4usl0eKlzJGXYKb1nbs5gFkBbBx3+NNmkO+ui+PcLuL+/X\n2+innoydqOIjDqeVNaqY5j1VUevC0PENd96eNRbZxiw6zAP0OIUEQJBjzE5HSo4TNGtpDW8kaFYh\ncOEzJpLjWoAmPFHXCiMH8hakIUuFQuEQucCLQoh64O/Aw1LKY4ncqApaRUYipQXNG6BlE5itSKMA\n8hYicufaXi90A2PaoqTmcKZ9s6kg/M6LcWIvpa2g1YjtqK2lwG176meu5p0v/jJqffxNS/GWpO/D\njVsfTqH7EppD6zBlxN5WE9nkuc7BrY9IW15OETSPcSLwFGACIAnTFt5KwDxKifcahFC/shSKTEZK\n+S3gW0KI2cBNwOtCiCNSyl7F19VPB0VmUv8KtHabDg83Qv3LSNOPKEjeMNFgQ7bHkftyWAosHvNL\nxvBs5fYYMeetbKd/4Tr8NQ3s+OljmP4AQtMYe/35nPvLf3F036AZ5mh7E7kuL8Ve+z5hrzEejz6W\nkHUciKg8iF4MJuIhpZ275MCgJbSBk8Vsd8LyBH5zD1lGclzGFArFgKeWiDvtCSChSWRV0CoyDhlu\nhtYYU+vNG5B5CwaEFWg60CfNxdpr/2+jTbI/vU6U/qg6jM4p5JJRU3n5NFWH4b5cLh/d/6n13myN\nhRAsvPd2Zn3tZpp2HiJ7dCnZ5c5q2j5XuZ0Xj+ygPRwCYGrBcD42ebHtAJwQGm49OeoSR555K6YL\nXMP289PWLyylRdCKbvs4ScCs7HNBO+hVHRSKIYIQ4jPAjcAw4GHgDinljkTuVQWtIvMIVGEr5wMg\ngxCsAe+YlKY0UHCdfzXht19A1vdsSdLKJ2LMX5amrCLcMGE+MwpHsP7YAfxmkKkFZSwZXoHP6LtB\ngGxpJPj4rwlveh3CIbRx0yO2xhPt+089BTmUnjOjz/slyutH9/D4wZ76srsaj/HTra/xHwsuR+vH\nCWxvHH4itinD4cfXpq2gBQHo2J3QRhjCtsYKxdBhNPAvUsrNZ3qjKmgVmYcWPTzSMx6tazpUEDn5\neL/4E0IvPIC5bT3oOsa8C3EtvwXhSv+p9fTCEUwvTE5vqDRN/L+4G3m0ywLcOrijy9Y4Tf3CAC9W\n2asqVPub2VZfzexi56b1NXfsH/vxYk4jhMCrj6fD3Gsb92WorbFCMdARQqwE/peIGdfvpZQ/PC3+\nFeDDRE6SXMA0oATIAe4HhhOZbP2dlPKn8faSUt7T+Z6lgLfb+uHe8lROYYrMwzsW9Bi6pa4ShNsZ\nY4B0I00/snkjsuF1ZOsOpLS3ztXyS/Dc+AWyvv03sr7xF9xX3YHw2dvcDmbMLWt7FLOnCIfSamsc\ntkyOd7TGjFf7mxzdf/xN9ifxQtMYd8OFju7dG3mus9BE9PeiV6/Ao49NQ0YKxdBGRBr2fw5cCswA\nbhFCTO1+jZTyv6SU86SU84F7gNVSykYgDHxJSjkDOAf47On32uy3SgixBzgAvA4cBJ5LJFd1QqvI\nOITQkcVXQN0/QIa6ApoPSi5PX2IOIv2HoO7xnn/fprXI4TcijNRM6Q80VYf4tsaxY05jaDqFniwa\nAu228VJvQrblfabsgtnM+totbP3h306tCU3j7J9+jtxxZY7u3Ru6lssw7/W0h3cQMKsQwoVPn4hX\nr0CkQO1CoVBEcRawR0p5CEAI8SBwNRBtYxnhFuBvAFLKGiKDXUgpW4UQO4FRce4F+C6wGHhZSjlP\nCLEMuC2RRFVBq8hIhG8MctTt0LoNws3gKoacGYje2hEGIdIKwfGnehazAOEmOPECDL8xPYmlmXjW\nxSK/uF/v3X0ADs68mL9o5BQePRDtTFbizXa03eAkC++9nYrbLuHw42vR3AbjbrjQ8WJWSsmuxmPs\nba4lx+Vl0bCx5Lii/z9qwkuOaz45ruQ41FW9+K7tEJzucdFRtwrvMOdtpBWKQcwooLLb10eIFLlR\nCCF8wErgszaxccBc4O1e9gtJKU8IITQhhCalfE0I8b+JJKoKWsWgImgeoz28Cws/Lq2ULGMqurB3\nzxJ6NuQ76zA1IPDvA6vDPtZxGBluRhh5qc1pAGAsuIjgU38Am5NQ13nptTVePmoqLaEOXq36gLCM\nmCaMzi7kjqlL0PsxEGZV7Se09mlkwzG0URUY512FFqOwL5wxLmUDYB1miF9sf53dTbWn1h49sIk7\npi5hTnG5o3vv/9srtutmIMTBR99k6j+l93tBoUgX6974gHVvJvVp1SpgTWe7wSmEEDnAI8BdUsrY\n/VYRGjuvfwP4qxCiFkjI51sVtIpBQ2toMy2hrg93AfMQ7aGtFHmvxKX178RtUGP5e4l3AEOvoBXZ\neXg//R0Cf/o+svlEZFF34Vp+c9oVHYQQXD9+HpeWT+dwaz25Li+jc/rnSBbe+BqB++8FK1Igm9vf\nJrTmKbyf+xF6+cRkpN1nnjy4pUcxCxCyTP5v11p+cNa1ZDs4kGh2BPsUUygynXMvmMK5F0w59fWP\nv/+M3WVVQHdZoPLONTtuprPd4CQi4obyCPBnKeUTCaR1NdABfJHIoFk+8O0E7lMFrWJwELaaehSz\nJ7HooCm4hhLv1WnIaoDgifOIWvOCqyh1uaQQKSWYraB5YuoK65Pm4PvWXyO2xh3tEVvj3IHziDnH\n5UmKqoMMBQk8/LNTxewp2lsIPvYrfF/4737v0R/W1+63XQ9aJhuPH+b8Ec4V3KOvPIcDf18dHRCC\n0auGrsmKQpEgG4CJQoixQDWRovWW0y8SQuQDFxIpQrvzB2CHlPIn8TYRQvwLsA7YKKU8qd33pzNJ\nVBW0ikGBP4aUD0DIqsG0WtC1vg3TDHbRdeEuRfoqIq0Hp5O3KCPtQmXLVmhaD2YzoCOzJkPRMoQe\n3X4yFGyNrX1boa3ZPrb3fWR7CyLL2WGzWFhS4g+HYsbbwgFH9x9/0zL2PfAKVc+902N9zr99mLyK\nkY7urVAMdqSUphDic8CLdMl27RRC3BkJy992XnoN8IKU8tQjQyHEEiIF7lYhxCYisl7/JqV83mar\nciLSYFOFEFuBtUQK3HVSyvpEcs2833SKjCSWBNVJLMJDW3Z92CpoeDMyBCcDoOdGitm85AzWDCRk\n61aof6HbigntOyF0Ajnitn7Zwg5aehUASJ9CgCYEFXnD2NtcZxuflO+sjJ5m6FzyxHc5+MjrVD7z\nFobXw4QPX8yIpf1zxlMohgqdBeiU09Z+c9rXf+K0E1Up5VoSdESRUn4FQAjhBhYC5wKfAH4rhGiU\nUk7v7T1UQasYFHj0UbSF7Y1DNJGDIVIjTZUOwlYzYdmILnJxafZ9lkIYULQMWXgBWEHQvBkpcySl\nhKYYQ7KhWvDvh6z09oumA23CLMjOh7ZoDVtt0lxEVv90hvur6rBq7Cx+su01LNnTwW9m4Ugq8vpu\nMSxDQcydG5DtLeiT5qAV27dvaIbOhJsvYsLNF/V5L4VCkRJ8RIY+8jtfR4GtidyoClrFoMCtjcKt\njSRoHY2K5boWZuSpnCWDNAVX02F2mQO4tREUuC9G1+yNI4TQQfelKsXUY/kh3Bg7HqgekgWtcLnx\n3HgXgT99D6xu1rHZeXiu+0z6EutkakEZ/zLzIp45vI19zXVkuzycO3wCV4yZ2ef3NPe8T8d934GW\nzu8HoWEsuQL3DV/IyA9zCkUmI4T4LRHjhhYi0l7rgB9LKRsSfQ9V0CoGBUIIijwraQm9R3v4AyQd\nGKKEHNc8fMaEdKfnCKcXswBBq5r6wHOUeK8fmr+0NTcIA2K1oNj00A4VjHkXoJWN6ZTtqo3Idi25\nEq2fmrvJYkrBcKYUDE/Ke0l/Gx3/9x/g76bmIy3Ca55CGzEO1/lDeEhUoRicjAE8wB4iKgpHgDin\nF9GoglYxaBDCRZ57MXnuxUhpJeVUNtAYWxKv5vX3KbtwTr/36AthqyWqmD0VkycIWkfx6H0T4B/M\nQ3BCGMjsadBq8wRKGJAd11VxUCKlBf690L4XpISsCsiabPv9r40Yh+dDn0t6DgNtCC686fWexWz3\n2LpnVUGrUAwypJQrReSUZgaR/tkvAzOFEPXAeinlN3p7D1XQKgYlyWoxOPTYm7YuQgB7i1amr6Dt\n5SlL2Groc0E76Cm8EILHIVjdtSYMKLk8YqaRQUhpQd0TPRUs2neCdwyy9LqMVLBIBNkSe+hZtiT8\nhFKhUAwgpJQS2CaEaASaOl9XEnEmUwWtQhGPUHO0i9RJgi0JmZM4gi7iSyzpon9DPoMZoXmRZbdC\nx0EIHAXNB9lTbSW7Bj1t2+3l2DoOQ8v7kLcg9TkNAPRx04klBKaNm9bv9+/vEJxCoTgzhBBfIHIy\ney4QolOyi4iOrRoKUyh6Y9SlC2PGyi9NrnbpmeDSCnFrIwha1VExXeTi0cfY3DW4kdKkPbwTv7kH\nKYO49ZFkG3MwtGiXMyEE+MZHXplM2674sSFa0GqT56FVzIro73bHcOFaHqX5rlAoBj7jgIeBL0op\no3/xJUDmjYYrFGdAwbSxTPv8tVHrpUtmMuG25WnIqIsC98UYoqTHmi5yKfRcmnGqDlJaNASepzm0\nlpBVS1g20h7ewfGOxwhZCWlqZybx9JdlbLOCTEcIgffO72FceB34ckAItIlz8H72P9E3YpGMAAAg\nAElEQVTHZl4ftUKR6UgpvySlfLSvxSyoE1qFgsU/+Rwjls1l/wOvEG4PUH752Uz6xEoMr3P+8omg\na9kM811PwKwibDWgixw8+ph+FbOhVn/vF6WBDvMgAetI1LokQEvoHYo8K9OQ1QDAOxYCMWzTfeNS\nmkoqkVJC6ARoLoRhrzEtvFl4rv8Mnus/g5QyqaofqsVAoRh8qIJWkdGYB3ci64+hjapAGz465nVj\nrzmPsdecl8LMEsejj0raANje+1+MOQRX/dpMRiybl5R9zpSAeTBO7DBSmhGN3TNkMCs6AJA7N+L+\nZp5ma6tnQ25mthvItp3Q8AaYLZGv3WVQvALhju0oNiQl7BQKRQ9UQavISKyGWgK//ybW4d2n1vQ5\n5+H5yL8i3N40ZpZeatdvjx1btyNtBa1E9n7REEToWciyW6DpLWjfA8iIcUT+YoQRf3BwMCL9B+D4\nMz0XgzVw7GHkyI9nnIqFQqFIHqqgVWQkgT9+p0cxC2C+v4Zg3u/w3PD5NGWVfrLKYovs+0YUpTCT\nnnj1cXSYe21jHm10n05nBwOmbCdk1qIJLy5tuO1JozByoXh55JXpNL1jv275oWULFJyT2nwUCsWg\nIbMmSxQKwKzah3Vwp20s/PYLyFAwxRkNHCbffjmaK/pzrKc4j/E3Lk19Qp149fG4tei2CoGbXPdZ\nacjIWaS0aAq+Sa3/rzQEX+BE4AnqOh4iZNWlO7X0EqrtW0yhUAx5VEGryDhkY5yiINiB9Md2B8t0\n8qeMZumD/46vrOs0Nm/SKJY/ey+uHF/a8hJCo8hzGbmuczBECbrIw6dPpcR7LS4t+datEf3u9NES\nepf28A7AOrVmykZOdDyLJQPpSyzd6HH0lePFFArFkEe1HCgyDm3URNB0sMyomCgchsgp6Nf7dxdd\nH4zT0GOvPY/RVy6m7p1d6B4XxQsmD4ihGiF0clyzyXHNTsr7BZvbuOTVvxO2UXY4+MgYxt9wYVL2\nOVMierv2vcySDvzhPWS7Zvbpve2G4AbNABxAzhxoeNUmICAnOd8XCoUiM1EFrSLj0ApKMBZfSnjd\ns1Ex1/JbEVpmPpgIrX+O8NqnkI0n0EZPwrX8ZvQJ9oWR5jIYvqRvRdNgofLpt2yLWYD9D76atoLW\nlO1IYre99GZ7nNHkzoNQHbR2M0wQBhRdjHAPS19eCoViwKMKWkVG4r7hLkRuIaE1T0FbM6J4BK7l\nt+A69/J0p+YIwSd+R+iVv5/62tx+AnPnBrz/9H30qZkp79QbZkfsotEKpM+UQBdeBAYSe9OE3myP\nMxkhBBRfisxbBP6DoLnANxGhp68dRqFQDA5UQavISISu477iE7gu+xgEO/5/e3ceJ1dZ5X/8e6qq\nt3Q2kpAEspKErISQAAmRrQFZBTOgbOoo8NNBGQdmRnCD+QmOijqMqIA/ZUZBFGQYQKMiCFECggQS\ntgAJGJYEskICIWtvVef3R1WS6u6qTi9Vdeve+rxfr35Z9Ty3qk6uIa/TT53nPFJNXVl8rV4MqS3v\nqmXhPTkmkmr+/c9UV6EJ7chTZ8sScXlrx9KTUacHt1verEp1iYmZGtp2c0qoLjEpgKiKz71VO5Ov\nqTm5TjGrUV3iQFXFhuS81qoGSVXBdd0AED4ktAgN91bp/SelbS+l2/jUjEz346zNf+iAxWJSbZ+C\nxlFudbOp116UkrlX+1JvviJv2imrqbwVrj77Ddah3/q0lnzxJ23Ghx87QwdeeHJAUaX1r5qrpO9o\nc6BETLUaWHOC4taz/6+SzeV7FG7Sd+jdxt+p1TfvHtveulR9qw5Xv6pZAUYGICpIaBEK7i69/Wup\ncdWewcY3pMZV8qEfkdWNCS64gFmfTprNV9VIiarSBVNmpl9+joYddZBeu+0hNW/drpEnH66x5zQo\nXh3sPTFLaFDNyWpJbVJzcoNiVqva+GiZ9fyf5MVf+LGkjhsel/9ovqZcMq8X0fbeluYn2iSzu2xr\nWaza+ChVxbpfHxv6U+AAFBQJLcKhcVXbZHa3lLT5MamCE9rYgYfI9tlX/l7HdmWJQ4+TxXv+n3l2\nRwepfFan3V2pV55R8vUXZf0GKjHrOFl9/5zXDj1iqoYeMbXEEXZNVWxwQdqStWzbqb/97H417OzY\n8uulZSMCTWjdW9WYfD3v/M7WFapiwxeAXiKhRTg0rsw/17xOnmqWxapLFk45sVhcNRf8mxp/cqW0\nY+vu8djoiaqe9w8BRlYc3rRTjT++UqnXlu4ea57/X6q56P8qMTV6hzB0xc717yqZI5mVpG1vrCtx\nNG2lN7+l8s5XdN9dAAVDQotwsM6+Io5JFs1WXF0VP2Cq+lxzh1qfWSh/f5NioycqPvmwSLYoa77/\ntjbJbHqwUU0//6biX7+zIuuF60ftq5rB/dW0aUuHucEzDwwgoj1iVquE7ZO3HVl1fL8SRwQgikho\nEQ71k6X3n8g91+fAXtUeljNPNUrblkktm6TEQKnvNFk89yY3q6lT1dxTC/r55VJikK118UO5J3Zu\nV/LFRUocelxpAyoD8ZpqTf/ieVrypZvbTpjp4Cs/HkxQWfpWHarNzQs6jMdtgOri4wOICEDURDML\nQORY1WD5wKPS9bLZEgOkfYJpkF9s3rReevuedEeHXd5/Qj70LFntyOACC9rO7XmnKvlY4+lXnKuq\n/n207Af3autrazXokPE6+Ksf15h5RwYdmuoS4yW5trU8ndkcFlNt/AD1r54r6/Tbl+7b+vpa9Ru3\nf0HfE0D5I6FFaNiAI+S1Y7Lado1Ir1jGaoIOrTg2/aFtMitJ3ixtvE8+4jOyCi2ziE+cqeSypzpO\nmCk+cWav3rucjzX2ndukVCrv5jdJmnzxGZp88RkljKrr6hITVJeYoKTvVEyJXiWyje9s1vEP3K5U\njlZlL+y3SR/4f//cm1ABhBAJLULFavaTaqJfc+dN66SWd3NPJrdKjW9VbGeHqlM/qeSrS9MHZmRJ\nzD1NsaHRW7lObXhTzffcpOQrz0juio2bpuozP6f4mMlBh9YjPe2zm+2dp17OmcxK0obHXsg5DiDa\nKnOJByh3qb3s/E41dj4fYfExk1X7Lz9QfNZxsn32VWzURFWfc5mqz43eqpxvf1+NP/yCki8/LblL\nklKvv6TGG7+o1Ma1AUcXnLph++SfG84JY0AlYoUWkVQOTdfX7dhzQtV+fS7u3otrhkuWkDzXCWCx\ndLlFRHnrFinVLFUNyltWER8xXvELriz4Z5dbmUHLogfkW3N0B2jaoZZH56vmrM+VPqgyMOSwSRp8\n6ERtevpvHeYmfeZDAUQEIGis0AJlyGK1Uv/Dck/2myFL9C1tQCXgzRvl638lrblZWnertOZm+bbK\n/vo4tfq1TuZeLWEk5ef4u6/WkMMm7X4er6vRzGsu0AHnNAQXFIDAsEILlCkbeJQ81kfa+rTU+r4U\n7yf1myn1Pzzo0ArOk43Shruk1I49g8lt0qY/yq1GVj8xuOACFBs0TMlO5nqjXE+B66q+Y4bpjKd+\npE3PvarGtzdryOGTVLNPv6DDAhAQElqgjFn/WVL/WXJPyiwedDjFs/2ltslsti1PSRWa0CbmnqqW\nh++WWtttgLKYEkeVZzeD3vJUSsnlTyn54iIpHldiZoPi46fnvX7wIRNKGB2AckVCCxRJt+tmOxHp\nZFaSmjd0Mvd26eIoM7Eh+6vmoq+p+Vf/uaeWtq6vqs/6nOJjpwQbXBF4Mqmmn12j5At/3T3W+uh8\nJY49SzUfuSTAyACUOxJaVJT3XlqpfaaNDToMtBfv5KviePTqhbsjcdARin/9V0q9ulSeSio+/mBZ\ndTR7L7c+9WCbZHb3+CP3KnHI0Z2u1AKobCS0iJyWbTt14sN3qWVrx6+wl9a/qWNvL/zu+KCFvqtD\n34OkLYslpTrO9ZvRu8DKlKeapW0vSk1rpFht+pCQmtwnXFk8ofikWQX9/HKsmU0+uzDvXOvTD5PQ\nAsiLhBaR8+5zr+ZMZiVp/V8qe9d8ubKqfeRDTpM2PdC2VVn91GhugmvdKm24M73Zb5dtz8sHHCkb\nODe4wALmrbna1GUkcx+k0FVh3wQHoHMktIic2s6arg8dWMJI0B1WP1leN1ba8Wq6D23taFn1kKDD\nKo73HmmbzO7y/uPy+omyqsGlj6kMJKbNUfOrz+eci087osTRAAgTElpEzoADR2r4cYdo/cPPdZib\n9A+nBxBRZSjEJjiL1abLDyLMPSnt6HggwG7bl0sDjypdQGUkceTpal3yJ6XWtO2/G586W/GDSGgB\n5EdCi0g69pdf1cJzv64Nj70oSYpVV2nqP51ZcQltsrlF8eqqoMNANk8qZ63wLqnefbUeZlbbR7WX\nXa/Wx3+v1hefkOIJJWYeq8QRp8hiEe/0AaBXzDPng4eRmXmY40fxvfvC69qxZqMGzzpQdUPzlyKE\nXb5NYV84IqGDv3x+iaMpvnLYBNcbvv4OqWlt7smhH5HVHVDagAC0YWZydws4Bl+7/ccl+7z96z8b\n+J+5N1ihRaj4zu1qfeGvUtNOxSfNUmzoyE6vHzR9nAZNH1ei6IKx8Zm/qeG3t+ace+31AyKZ0JaD\nXnV1GHiUtOFudViprR0t1Y7tdWzlyJvWSpsflxrflKxKqp8iDTxSFu8TdGgAIoCEFqHRuvRxNf3i\nO1JTpoOBmRJHnq7qsy+VWWh/qey15I6mTuYaSxgJuspqR8uHnSu9v6hN2y71nx3Jv8vetE5a/z/S\nroN8vVna9rzUtEY+/OOyGGUxAHqHhBah4FveU9PPvym1NGcNulof+51ioyaqau6pwQUXsCGzJ6t2\n34FqfGdzh7mRp7GRplxZ7Qip9iNBh1Ea7/9Vu5PZbC0bpR0vS33pLwugd2JBBwB0ReuSP7VNZrPn\nnri/xNGUl3h1leb88POyeNv/nPuN318Hf5Vem8WyX5+Ld/9gLxrf6tkcAHQRK7QIBd/WcfVx99z2\nHP08uynsTdfHnXucBk4erVf+6z7tXP+u9j1iqiZ++jTVDIzesbHrFnZsxxZlYd8AJ0my6rYHZrSf\nA4BeIqFFKMTGHyQtyD0XHxfxvqWtLfLtW2R9B8ji+f+THTRjvObeeGkJIwvGS9/7XzX8flHHCTNt\nfeOD6nfAfqUPKuJ6tQFOkvpOlbYsyT8HAL1EQotQiE+ZrdjYqUqtXNZ2oqZOVcefE0xQReappFru\nv00tf/mttGOrVD9AVQ1nquqkj0dy41BXbX45z1fU7tqyYg0JbTkaMFdqXC01r2873n+OrGb/YGIC\nECkktAgFi8VUe8m31Xz/bWpdskBqTLftqjrtU4rtNybo8Iqi+Tc/UevCe/cMbH9fLffdKrW2qvpD\nFwQVVuAGTh2jra+u6TBusZgGTB4VQETYG4vVyIefL+1YITWuSpcZ1E+R1QwPOjQAEUFCi9Cw2j6q\nOfOzqjnzswV/73KrmfUdW9X6+O9zzrU8cq+qPniurKauxFGVh+mXn6PV9y2SJ9v2cD3g/OPUd/Sw\ngKKKtoIca2xxqX5y+gcACowuB0AZSm14K29XBzXukG9an3uuAgw7arpOmP8NDTpkgiSpekC9pv3L\nR3XUf18ecGSllWrN0QYLACpUYCu0ZvZdSWdIapL0mqQL3X1LZu4rki6S1CrpMnd/MKg4gSDYPkMl\ni0me6jgZT8gGDO7V+4e9q8Oo0+Zo1Glz1LJ9p+K11YrF40GHVBTvvfhG3lPgVt4zUuPOPa60AZVA\nuXR16PVGOAAlFeQK7YOSprn7IZJWSPqKJJnZVEnnSJoi6VRJP7JK3gGDihQbOETx6R/IOZc49DhZ\nff8SR1QarS88ocZb/l2NP/6qWh6+R964o9Prq+rrIpvMStJb9+Xo5pCxupM5AKg0ga3Qunt2E6ZF\nknYdmfNhSXe6e6uklWa2QtJsSU+WOEQgUDUfu1xNrc1KLnsqPWCm+PQjVX12NFtzNd19g1ofnb/7\neXLZU2pZdL/qLrte1qdfgJEFJ1Gfv0460ae2hJEAQHkrl01hF0n6VebxCElPZM2tyYwBFcX69FXt\nZ7+l1NurlXpntWJDRyu2b2FaHJVbiUHyrRVtktldfN1KtTx8t6o/dGHpgyoDB5x9jBZf/mOlmls6\nzI3/xAcDiKhyUGYAhEtRSw7M7CEzW5r180Lmf8/IuuZKSS3u/qtO3gqoWLGhI5WYdkTBktlylFz6\neN651k7moq5u2CAdfesXFa/dc5qWxWKaec0FGnbU9AAjA4DyUtQVWnc/sbN5M7tA0mmSjs8aXiMp\nu5nkyMxYTldfffXuxw0NDWpoaOh+oACC1VmZfC9L6MO+AW7cecdr/xMP1Zu/fkzJphaNOn2u+o6J\nZnuyxnfyH3G9Y/276jN8UAmjQSVZuHChFi5cGHQYHVz9dGWWW/VEkF0OTpF0haRj3L0pa+q3km43\ns+uVLjWYIOmpfO+TndACCKfEjKPV8kDu3e2JGUeXOJrS8eYmpda8JuvTV7Fho/NeVzt4gCZ++kMl\njCwYr/7iobxdHV6dm9DBXzq/tAGhYrRfELvmmmuCCwY9EmQN7Q2SqiU9lGlisMjdL3H3ZWZ2l6Rl\nklokXeLuHmCcAIosNmKcqk44Vy1/+p+246MmqqrhI3leFW4tj/xazffflj7WWFJszGTV/P2XFRs6\nMuDIgrNjzcb8c6vzzwFAkF0ODuxk7lpJ15YwHAABq573GcWnHK7WJQvkjTsUn3yYEoedIKuuCTq0\ngmt99lE133NTm7HUqpfV+KMvq+6qW2SJqoAiC9a+c6Z0MscJY0AQMt+of1/pfVc/dffv5LimQdL1\nkqokvePux2XGB0j6b0kHSUpJusjdi9K1qly6HAB7RcP16ItPPETxiYcU9D3LsWa25ZF7co77u+uV\nXPq4ErMaShtQmRhz5lEaPOtAbXpmRZvxQTPGa+zZxwYUFVC5zCwm6UZJJ0haK2mxmc1395ezrhkg\n6SZJJ7n7GjMbkvUWP5D0B3c/28wSkvoUK1YSWgAosdQ7efe5djrXFdmb4Moxme9MrCqhkxdcp6Xf\nul0r7/mLlEppzFlHa8ZVn1C8pnrvbwCg0GZLWuHuqyTJzO6UNE/Sy1nXfEzSPe6+RpLcfWPm2v6S\njnb3CzLjrZK2FCtQEloAKLHYsDFKbc29oz82fEyJoykvNQP76vDvXqzDv8u3H0AZGCHpraznq5VO\ncrNNlFRlZg9L6ivph+7+C0kHSNpoZrdImiFpiaTL3H1nMQIloQW6iTID9FbV8Wer6dXnO4zb0FGK\nHzQ3gIhKI/nKM2p96UlZIqH4zAbFR+XdSgGgyNY+u0xrn11eiLdKSJqldAvWeklPmNkTWeP/6O5L\nzOz7kr4s6WuF+NBcQQAASihx0BHyj39RLffdIt/8TvpY40mHqvr8f5XF40GHV3CeSqrp1m8q+dyj\nu8daFvyPqk76mKpPvyjAyIDKtf/Mqdp/5tTdz5+59d5cl62RlN1TMNfZAKslbXT3RkmNZvao0iuy\nj0l6y92XZK67W9KXChN9RyS0CIVkU3PeucZN76t28IASRlM6uTbCsQkuGqrmnKTE4SfI31kj1dYr\nNmBwQd63HOtmWxcvaJPM7tLy4B2KHzRX8bH5uxsACNRiSRPMbIykdZLOk9S+IfR8STeYWVxSjaQ5\nkr7n7hvM7C0zm+juf1N6Y9myYgVKQotQWHnPX/I2XF9xS42mX35OaQMCCsBicVknBypEReuSP+ef\ne/phElqgTLl70sw+L+lB7WnbtdzMLk5P+83u/rKZ/VHSUklJSTe7+67E9VKlD8uqkvS6pAuLFSsJ\nLUJh87JVeefeX55/DkAZaGnqZK6xV28d5q4OQBi4+wOSJrUb+0m759dJui7Ha5+XdHhRA8wgoUUo\nDJySfxVrwOTor3AFiTID9FZ8yuFKvf5inrn2G6YBoPtiQQcAdMXYjx6jfuP26zBeu+9AHXjhKQFE\nBKCrqo6el7O0IjZpluLTo9vVAUDpsEKLUIjXVOuUP/+nFl16o1bft0iecu13/EzNuf4S1Q6J5oaw\nlXc/knP83edf06AZ40scDdBz1qev6v75+2p55NdKvrhISlQpMfNYJY46QxbrXVcHygwASJK5e9Ax\n9JiZeZjjR8+0bN8ppVxV/Yp2gl7gPJXS3eM/oW2rNnSYG3XGXH1w/jcCiKq4yuFoYzo6AJAkM5O7\nW8Ax+D88envJPu/mYz4e+J+5N1ihRehU1dcFHULRbVmxJmcyK0lrFzxT4mgAAChv1NACZah6YL1k\nuX9Rrtmnb4mjAQCgvLFCCxSJr9rTwcTGXN6t19YNG6SRp83W6vue7DA34QI2wRULZQYAEE6s0AJl\n6sj/ulxDDm/T+k9jPnK0Dvm3TwQUEQAA5YkVWkRWdsN1KXy7ofsMH6QznvyRNvz1JW17Y50Gz5yg\ngVPHBh1WUTRv2Z53bs0fF2vEySXpy11SbIIDgMIhoQWKpLtlBvkM+8A0DfvAtIK8V7l6c/5f8x5t\n/GrfEyKZ0AIACoeSAwCBa922M+9cSydzAABIrNAiwsJWYhB2vdkEN+Lkw9JdHXL0lR55CkejFgtl\nBgCighVaAIHrN25/Tf/ieR3G9507VRM+dVIAEQEAwoQVWiAiwr4J7rBrP63hx0zXa79coJZtOzXy\n1Dma8KmTlKirCTq0gtv21ttBh1CR2AQHRBcJLYCCKMQmuJGnztHIU+cUIJrytvyGX6vht3flnHt3\n6Qc06ODxJY6o+MqhqwOA6KLkAABKbNMzK/LObVzytxJGAgDRwAotAJRY/eiheef6jhlWwkgqC2UG\nQHSR0AIREbaa2bDrTVeHKZfM02u/eEieTLUZHzhtrPY7fmZB4gOASkLJAYBQuvgvd+z+CZshh03S\nsXdcpfpRe1ZqhzfM0In3fUtmFmBkABBOrNACQAAOOPtYjTnrKG1etkrVA+rVd3R0Sw12rNuUd27r\nyvXqN3Z4CaMBEEUktADQA4Xo6hCLxzVo+rgCRFPeXvvlgvxHG890zfzaJ0sbEIDIoeQAAFBUje9s\n7tEcAHQVK7SAaLgeRmyCK53sDXBS91enhx9zsF68Lnff3eHHzuhxXACwCyu0AICiGnnaHA0/7pAO\n40OPPEhjzjwqgIgARA0rtABQ5sJ+rLHFYjrxvmu1/IZfa+Xdj8hTrrFnHa0pl56pWCIedHgAIoCE\nFhBlBkBnCrEBLlFbrelXnKvpV5xbgIgAoC1KDgAAABBqrNACFYhNcOESthKDKOjNSXAASo8VWgAA\nAIQaK7QAgLIS9k1wAEqPhBaoQJQZAJ2jzAAIF0oOAAAAEGoktAAAAAg1Sg4ARB5dHcKFmlkA3cUK\nLQAAAEKNhBYAAAChRskB0A3ZzdYldkKHBWUGABBtrNACAAAg1FihRajQcB3oHBvgAFQiElqgGygx\nAACg/FByAAAAgFBjhRZAp7I3wrFCXf4oMwBQiUhoESrUzAIAgPZIaIEyxiY4oHNsggMgkdAC2AvK\nDAAA5Y5NYQAAAAg1VmgBAKFFmQEAiYQWKGvUzJYWHR0AIJwoOQAAAECosUILFMn2Sz+4+3H9DxcE\nGEl40NUBANATJLQAkEGZAQCEEyUHAAAACDVWaBFZNFwPH0oMSotNcACigoQWKBLqZgEAKA0SWgAI\nCJvgAKAwSGgRWZQZlBZdHcKHMgMAUcGmMAAAAIQaK7RARGRvgpNYoQbaYxMcEF0ktAAKgjKD7qNm\nFkC5M7NTJH1f6W/1f+ru38lxTYOk6yVVSXrH3Y/r6msLhYQWAAAAHZhZTNKNkk6QtFbSYjOb7+4v\nZ10zQNJNkk5y9zVmNqSrry0kElogIigxKC02wXVf0F0dKDMAum22pBXuvkqSzOxOSfMkZSelH5N0\nj7uvkSR339iN1xYMm8IAAACQywhJb2U9X50ZyzZR0iAze9jMFpvZ33fjtQXDCi0AAECFWfvsMq19\ndnkh3iohaZak4yXVS3rCzJ4oxBt3NwgACB2ONg4fNsEB5WP/mVO1/8ypu58/c+u9uS5bI2l01vOR\nmbFsqyVtdPdGSY1m9qikGV18bcGQ0AJAD1A3C6ACLJY0wczGSFon6TxJ57e7Zr6kG8wsLqlG0hxJ\n35P0ShdeWzAktACAvWITHFB53D1pZp+X9KD2tN5abmYXp6f9Znd/2cz+KGmppKSkm919mSTlem2x\nYiWhBUTD9TCizAAAis/dH5A0qd3YT9o9v07SdWon12uLhS4HAAAACDVWaAEgBNgEBwD5kdACoswA\n2BvqZgGUMxJaAADaYRMcEC4ktEAFYhNc+FBmAAD5sSkMAAAAocYKLQCg7LAJDkB3kNACFYgyA6Bz\n1M0C4ULJAQAAAEIt8ITWzL5gZikzG5Q19hUzW2Fmy83spCDjAxB+vuq63T8AgOgJtOTAzEZKOlHS\nqqyxKZLOkTRF0khJC8zsQHf3YKIEAJQadbMAuiPoFdrrJV3RbmyepDvdvdXdV0paIWl2qQMDAABA\nOAS2QmtmH5b0lru/YGbZUyMkPZH1fE1mDAgczdYBACg/RU1ozewhScOyhyS5pKskfVXpcgMAKCq6\nOgBAtBU1oXX3nAmrmR0kaayk5y29PDtS0jNmNlvpFdnRWZePzIzldPXVV+9+3NDQoIaGht6GDQAA\nKsjChQu1cOHCoMPooPE/k0GHEBpWDnutzOwNSbPc/T0zmyrpdklzlC41eEhSzk1hZsZesQoTdLN1\nSg5Q7jjWGOg9M5O7296vLGoM/sl5t5Xs826b/8nA/8y9US4HK7jS5Qhy92VmdpekZZJaJF1C1opy\nQRILAED5KYuE1t3HtXt+raRrAwoHAAAAIVIWCS2A8pRdYiGxQh0GlBkAqEQktAgVmq0DAID2SGgB\nAKHV/jhjVqiBykRCC5Sx7K4OUulXqCkxAACEQdBH3wIAAAC9wgotAGTQZzh8KDEAIJHQAmWNTXAA\nAOwdJQcAAAAINVZogSK4JXZCm+cXpv4UUCThwiY4AEBPsEILAACAUGOFFgAqFFZDH9sAAAtZSURB\nVJvgAEQFCS0iK7vheql3QlNi0DNsggMA9AQlBwAAAAg1VmgBFET2RjhWqLuGTXAAUBgktIgsGq4D\nAFAZSGgBABUhexOcxAo1ECUktEBEBLkJDgCAIJHQAigI6ma7j64OAFAYdDkAAABAqLFCC0QEZQal\nRVeHnsnu7EBXBwCFwgotAAAAQo0VWgChxCY4AMAuJLQA0AOUGfQMG+EAFAMlBwAAAAg1VmgBtW24\nzsaRcKDMoLTYBAegnLFCCwAAgFAjoQUAAECoUXIAiDIDlL+guzpQZgCgnLFCCwAAgFBjhRYAgBzY\nCAeEBwktUIHo6hA+dHUAgPwoOQAAAECosUILACg7QW+CkygzAMKEhBaoQJQZAACihJIDAAAAhBor\ntAAij01w4cMmOADdwQotAAAAQo2EFgAAAKFGyQHQDTRaDyfKDAAg2lihBQAAQKixQgsAEZK9AU5i\ndRpA75jZKZK+r/Qi6E/d/Tvt5o+VNF/S65mhe939G1nzMUlLJK129w8XK04SWoQGjdYBACidTDJ6\no6QTJK2VtNjM5rv7y+0ufbSTZPUyScsk9S9epJQcAAAAILfZkla4+yp3b5F0p6R5Oa6zXC82s5GS\nTpP038ULMY0VWgB5sQkufCgxAFBAIyS9lfV8tdJJbntzzew5SWskXeHuyzLj10u6QtKAokYpElqE\nCI3WAQAojPUbl2v9xuWFeKunJY129x1mdqqk30iaaGanS9rg7s+ZWYPyrOIWCgktACDUOAkO6L7h\nQ6Zo+JApu58vfeU3uS5bI2l01vORmbHd3H1b1uP7zewmMxsk6QOSPmxmp0mqk9TPzG5z908W7k+x\nBwktUMaC3ghHmQEAVLTFkiaY2RhJ6ySdJ+n87AvMbJi7b8g8ni0p5u7vSvpq5mdXJ4QvFCuZlUho\nAQAAkIO7J83s85Ie1J62XcvN7OL0tN8s6aNm9jlJLZJ2Sjo3iFhJaAEgg01wANCWuz8gaVK7sZ9k\nPb5J0k17eY9HJD1SlAAzSGiBMsZGOGDvqJsFQB9aAAAAhBoJLQAAAEKNkgMAZSO7q4NU+pIL6mYB\nIJxYoQUAAECosUKLSMputC6xaQTIha4OAKKChBZA2aCrAwCgJyg5AAAAQKixQotIosQAYRHk8caU\nGQCIClZoAQAAEGqs0AIRkb0RjhVqIDc2wgHRREILAAFiIxwA9B4lBwAAAAg1VmgBABWDMgMgmkho\ngYigbhZhEPTxxgCiiZIDAAAAhBortABCia4OAIBdSGgBACVDiQGAYiChRcWjLyUAAOFGQgsglCgz\nAADswqYwAAAAhBortAAQAmyCA4D8SGhR8aibBQAg3Cg5AAAAQKixQgtUmOyuDhIr1GFBmQEA5McK\nLQAAAEKNFVoAQFliIxyAriKhBSoMJQYAgKih5AAAAAChxgotgEhjExwARB8JLQCgLFE3C6CrKDkA\nAABAqLFCCyDSKDEAgOhjhRYAAAChxgotAERM9kY4VqgBVIJAV2jN7J/MbLmZvWBm384a/4qZrcjM\nnRRkjAAAAChvgSW0ZtYg6QxJ0919uqTrMuNTJJ0jaYqkUyX9yMwsqDghLVy4MOgQJKVPDdr1EzXl\nco+jjHtcfNzj0uA+Ax0FuUL7OUnfdvdWSXL3jZnxeZLudPdWd18paYWk2cGECIl/PEuBe1x83OPi\n4x6XBvcZ6CjIGtqJko4xs29J2inpcnd/WtIISU9kXbcmMwYA6ALqZgFUmqImtGb2kKRh2UOSXNJV\nmc/ex92PMLPDJf2vpHHFjAfhRpN1AACQi7l7MB9s9gdJ33H3RzLPV0g6QtJnJMndv50Zf0DS19z9\nyRzvEUzwAAAg0tw90P07ZrZS0pgSfuQqdx9bws8rqCBLDn4j6XhJj5jZREnV7r7JzH4r6XYz+57S\npQYTJD2V6w2C/ssGAABQDGFOLoMQZEJ7i6SfmdkLkpokfVKS3H2Zmd0laZmkFkmXeFDLyAAAACh7\ngZUcAAAAAIUQyqNvzezrZva8mT1nZgvMbGTWHIcyFICZfTdzD58zs3vMrH/WHPe4AMzso2b2opkl\nzWxWuznucQGZ2Slm9rKZ/c3MvhR0PFFgZj81sw1mtjRrbB8ze9DMXjGzP5rZgCBjDDszG2lmfzaz\nlzIHEF2aGec+F4iZ1ZjZk2b2bOY+fyszzj0OmVAmtJK+6+4z3P0QSfMlfU2SzGyqOJShUB6UNC1z\nj1dI+orEPS6wFySdKemR7EEOFyksM4tJulHSyZKmSTrfzCYHG1Uk3KL0Pc32ZUkL3H2SpD8r8+8G\neqxV0r+6+zRJcyX9Y+bvLve5QNy9SdJx7j5T0sGSjjezI8U9Dp1QJrTuvi3rab2kTZnHHxaHMhSE\nuy9w91Tm6SJJu1bBuccF4u6vuPsKpdvZZeNwkcKaLWmFu69y9xZJdyp9j9EL7v6YpPfaDc+T9PPM\n459L+ruSBhUx7r7e3Z/LPN4mabnS/xZznwvI3XdkHtYonRe9J+5x6IQyoZUkM/uGmb0p6QJJ12aG\nR0h6K+syDmUojIsk/SHzmHtcfNzjwmp/P1eL+1ksQ919g5ROxiQNDTieyDCzsZIOUXqBYRj3uXDM\nLGZmz0paL2mhuy8T9zh0guxy0KlODmW40t1/5+5XSboqUw/3fUkXBhBmqO3tHmeuuVJSi7v/KoAQ\nQ68r9xiIMHYdF4CZ9ZV0t6TL3H1bjh7s3OdeyHwbOTOzV+SPZtagjveUe1zmyjahdfcTu3jpHdqz\nerhG0qisuZGZMeSwt3tsZhdIOk3pfsG7cI+7oRt/j7NxjwtrjaTRWc+5n8WzwcyGufsGMxsu6e2g\nAwo7M0soncz+wt3nZ4a5z0Xg7lsyhz4dJu5x6ISy5MDMJmQ9/TtJz2Ue/1bSeWZWbWYHqJNDGdA5\nMztF0hWSPpwpmt+Fe1wc2XW03OPCWixpgpmNMbNqSecpfY/Re6aOf3cvyDz+lNKbdtE7P5O0zN1/\nkDXGfS4QMxuyq4OBmdVJOlHSs+Ieh04o+9Ca2d2SJkpKSnpd0ufc/e3M3Fck/R+lD2W4zN0fDCzQ\nEMscRVytPRvuFrn7JZk57nEBmNnfSbpB0hBJmyU95+6nZua4xwWU+QXtB0r/Ev/TXUdro+fM7A5J\nDZIGS9qgdLeZ30j6X6W/YVgl6Rx33xxUjGGX2W3/qNIdUTzz81Wlf8G9S9znXjOz6Upv+jKl/334\nhbtfZ2aDxD0OlVAmtAAAAMAuoSw5AAAAAHYhoQUAAECokdACAAAg1EhoAQAAEGoktAAAAAg1EloA\nAACEGgktAAAAQo2EFkDoZU4B22lmzxTp/X9pZpvM7KxivD8AoHdIaAFExQp3n1WMN3b3T4ijLwGg\nbJHQAihrZnaYmT1vZtVmVm9mL5rZ1L28ZoyZLTezW8zsFTO73cxONLPHM88Py1z3NTO71cweNbM3\nzOwsM/sPM1tqZn8ws3j22xb1DwoA6DESWgBlzd2XKL06+k1J31H6rPVlXXjpeEn/4e6TJE2SdJ67\nHynpCklXZl03TlKDpHmSfinpIXc/WFKjpA8V6s8BACieRNABAEAX/LukxZJ2SvqnLr7mjazE9yVJ\nCzKPX5A0Juu6+909ZWYvSDJ3fzDrurG9ihoAUBKs0AIIgyGS+krqJ6m2i69pynqcynqeUttf5psk\nyd1dUku71/BLPwCEAAktgDD4saSrJN0u6btdfE1Pal6pkwWAEGL1AUBZM7O/l9Ts7neaWUzS42bW\n4O4L9/JSz/O4q68BAISEpb9lA4DwMrMxkn7v7tOL+Bm3SPqdu99brM8AAPQMJQcAoiApaUAxD1aQ\ndIzSnQ8AAGWGFVoAAACEGiu0AAAACDUSWgAAAIQaCS0AAABCjYQWAAAAoUZCCwAAgFD7/38igrlc\nW1RIAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "detector = spectro.predict_positions(simcfg)\n", "ax = detector.plot(modes=(-1, 0, 1, 2), blaze=True)\n", "ax.set_aspect('auto')\n", "ax.axis(N.array([-2000, 2000, -4000, 4000]) *\n", " spectro.detector.pxsize / 1e-3) # [mm]\n", "ax.figure.set_size_inches(12, 10)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This sequence of actions is readily available under `snifs.plot_SNIFS()` function." ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.6" } }, "nbformat": 4, "nbformat_minor": 0 }