diff --git a/suncet_instrument_simulator/Diffraction/diffraction_pattern.ipynb b/suncet_instrument_simulator/Diffraction/diffraction_pattern.ipynb index f9b6a17..2a110d6 100644 --- a/suncet_instrument_simulator/Diffraction/diffraction_pattern.ipynb +++ b/suncet_instrument_simulator/Diffraction/diffraction_pattern.ipynb @@ -7,8 +7,8 @@ "metadata": { "collapsed": true, "ExecuteTime": { - "end_time": "2024-04-23T04:18:42.793871Z", - "start_time": "2024-04-23T04:18:42.186058Z" + "end_time": "2024-05-09T21:04:51.944773Z", + "start_time": "2024-05-09T21:04:51.451208Z" } }, "outputs": [], @@ -35,26 +35,26 @@ "execution_count": 2, "outputs": [], "source": [ - "mesh = suncet_psf.filter_mesh_parameters(171, angle_arm = [0, 90] * u.deg)" + "mesh = suncet_psf.filter_mesh_parameters(171, angle_arm = [0, 90] * u.deg, angles_focal_plane=[90, 0] * u.deg, lpi = 20)" ], "metadata": { "collapsed": false, "ExecuteTime": { - "end_time": "2024-04-23T04:18:42.797813Z", - "start_time": "2024-04-23T04:18:42.795314Z" + "end_time": "2024-05-09T21:04:54.304942Z", + "start_time": "2024-05-09T21:04:54.300111Z" } }, "id": "4118077f978a9939" }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 4, "outputs": [ { "data": { - "text/plain": "{'wavelength': ,\n 'angle_arm': ,\n 'error_angle_arm': ,\n 'spacing_e': ,\n 'mesh_pitch': ,\n 'mesh_width': ,\n 'spacing_fp': ,\n 'width': ,\n 'CDELT': [, ],\n 'Area': }" + "text/plain": "{'wavelength': ,\n 'angle_arm': ,\n 'angles_focal_plane': ,\n 'error_angle_arm': ,\n 'spacing_e': ,\n 'mesh_pitch': ,\n 'mesh_width': ,\n 'spacing_fp': ,\n 'width': ,\n 'CDELT': [, ],\n 'Area': }" }, - "execution_count": 3, + "execution_count": 4, "metadata": {}, "output_type": "execute_result" } @@ -65,12 +65,173 @@ "metadata": { "collapsed": false, "ExecuteTime": { - "end_time": "2024-04-23T04:18:42.985057Z", - "start_time": "2024-04-23T04:18:42.983103Z" + "end_time": "2024-05-09T21:02:58.041513Z", + "start_time": "2024-05-09T21:02:58.037435Z" } }, "id": "63d883175cec901a" }, + { + "cell_type": "markdown", + "source": [ + "20 vs 65 LPI comparison to demonstrate improvement resulting from this optimization" + ], + "metadata": { + "collapsed": false + }, + "id": "97529b3bb0ede6d8" + }, + { + "cell_type": "code", + "execution_count": 3, + "outputs": [], + "source": [ + "psf_185_65_lpi = suncet_psf.psf(185, diffraction_orders=np.arange(-100, 100, 1), angle_arm = [90, 0] * u.deg, angles_focal_plane=[0, 90] * u.deg, lpi = 65, output_size = [100, 100])\n", + "psf_185_20_lpi = suncet_psf.psf(185, diffraction_orders=np.arange(-100, 100, 1), angle_arm = [90, 0] * u.deg, angles_focal_plane=[0, 90] * u.deg, lpi = 20, output_size = [100, 100])" + ], + "metadata": { + "collapsed": false, + "ExecuteTime": { + "end_time": "2024-05-09T21:04:57.802332Z", + "start_time": "2024-05-09T21:04:57.720752Z" + } + }, + "id": "10dc90fc3f0c5dfe" + }, + { + "cell_type": "code", + "execution_count": 4, + "outputs": [ + { + "data": { + "text/plain": "" + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAGgCAYAAAAD9NhnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAAArEklEQVR4nO3df2yV1f3A8U9py+0P2ssv29tCwTapoqLRgZoBERK1yyTLDIlT8Wf8RwcqlUSE4WZHpFWWEaJMHGZhLI5plrGMLftB549GwzYZDmVgcLMVGmgtSG1LW1ranu8ffHl2zqftvdxL4dxb3q/kJs+557nPPfe5z+3pOZ/znJNmjDECAIAHY3wXAABw6aISAgB4QyUEAPCGSggA4A2VEADAGyohAIA3VEIAAG+ohAAA3lAJAQC8oRICAHhzwSqhV155RUpLSyUrK0tmzZol77333oV6KwBAisq4EAd98803pbKyUl555RWZO3eu/PSnP5VvfvObcuDAAZk2bVrU1w4MDMjRo0clLy9P0tLSLkTxAAAXkDFGOjo6pLi4WMaMidHWMRfATTfdZB577DHnuRkzZpiVK1fGfG1jY6MRER48ePDgkeKPxsbGmH/zR7wl1NvbK3v27JGVK1c6z1dUVMiuXbsG7d/T0yM9PT1B2vz/pN6NjY2Sn58/0sUDRkx/f7+TPnHihJM+fPhwsF1fX+/kNTU1OemioqJgu6yszMkrKSlx0hMnTgy2MzIuSGcGcF7a29ulpKRE8vLyYu474lfw8ePHpb+/XwoLC53nCwsLpbm5edD+NTU18sMf/nDQ8/n5+VRCSGq6Ejp9+rSTHjduXLCdk5Pj5GVlZTlpO99+nYgM+iHbvwsqISSzcwmpXLArWL+5MWbIAq1atUqWL18epM/WoECyS09Pd9L6n6ZJkyYF293d3VGPNWXKlCFfN9RxqXgwmoz41Tx58mRJT08f1OppaWkZ1DoSEQmFQhIKhUa6GACAFDDiQ7THjh0rs2bNktraWuf52tpamTNnzki/HQAghV2Qdv3y5cvlgQcekNmzZ8vXv/512bx5sxw+fFgee+yxC/F2gBdnB9Gc1dnZ6aS/+uqrYPuLL75w8nRPgd3FpmNCujvO7jmIOfwVSHIXpBK6++675csvv5Q1a9ZIU1OTzJw5U/74xz/K9OnTL8TbAQBS1AWLcC5ZskSWLFlyoQ4PABgFaMsDALxhrCeQoI6ODid97NgxJ33kyJFg+/PPP3fy9M2qfX19wXa0e4hERDIzM4PtcDh87gUGkhAtIQCAN1RCAABvqIQAAN4QEwISZMdmRAbHcuz7fXJzc528aGl9n5A+rn5fIJXREgIAeEMlBADwhu44IEEnT5500vY0PSLu1DwtLS1Onk7bM3JPmDDByRs/fryTtods6+HbQKqhJQQA8IZKCADgDZUQAMAbYkJAgvRijDo+Yy/LrYddR0vr5bxjLQ0OpDJaQgAAb6iEAADeUAkBALwhJgQkqLu720nrpR1aW1uD7ba2NidPp+24j/06kcH3I3V1dQXbeulvINXQEgIAeEMlBADwhkoIAOANMSEgQbGW4bbjNfHcJ6TjPNnZ2VHfF0hltIQAAN5QCQEAvKE7DkhQb2+vk+7p6XHS9lBqeztWWufp4+r3BVIZLSEAgDdUQgAAb6iEAADeEBMCEqSXctBDp+1h17m5uU5etLTO00O09fsCqYyWEADAGyohAIA3VEIAAG+ICQEXQVpa2jnnx9oXGE1oCQEAvKESAgB4Q3cckKCMDPfno4dO27Nq62HWesZtO1/njR07Nur7AqmMlhAAwBsqIQCAN1RCAABv6FwGEqRjN5MnT3bSp06dCrb18gs6fjR16tRgu7i4OOpx9bQ+QCqjJQQA8IZKCADgDZUQAMAbYkJAgowxTlrHfbq7u4PtkydPOnnR0vbrhjruwMBAsD1mDP9HIrVxBQMAvKESAgB4Q3cckKD29nYnfezYMSfd2NgYbNfX1zt5TU1NTtructPT9Ogpf+xpe8aPH3/uBQaSEC0hAIA3VEIAAG+ohAAA3hATAhKkV0DVw6XT09OH3I6V1nn6uAzLxmjC1QwA8IZKCADgDZUQAMAbYkJAguylGkQGT8XT1tYWbOt7inTaXp7Bfp2ISGdnp5O2p/XJz8+Po8RA8qElBADwhkoIAOANlRAAwBtiQkCCYt0nZKdj3euT6L5AquNqBgB4QyUEAPCG7jggQf39/U769OnTTtpenkHnRUvrlVT1vn19ffEXFkhStIQAAN5QCQEAvImrEqqpqZEbb7xR8vLypKCgQO688045ePCgs48xRqqqqqS4uFiys7NlwYIFsn///hEtNABgdIirEqqrq5OlS5fK3//+d6mtrZW+vj6pqKhwphVZt26drF+/XjZu3Ci7d++WSCQit99+u3R0dIx44YFkkpaWdkEe0d4HSHVxDUz485//7KS3bNkiBQUFsmfPHrnlllvEGCMbNmyQ1atXy6JFi0REZOvWrVJYWCjbtm2TRx99dNAxe3p6pKenJ0jrObUAAKPXecWEzk60OHHiRBERaWhokObmZqmoqAj2CYVCMn/+fNm1a9eQx6ipqZFwOBw8SkpKzqdIAIAUknAlZIyR5cuXy7x582TmzJkiItLc3CwiIoWFhc6+hYWFQZ62atUqaWtrCx6NjY2JFgkAkGISvk/o8ccfl48//ljef//9QXm6r9oYM2z/dSgUklAolGgxAG9iTadjX/P6+o+WjjUdEDCaJHR1P/HEE7Jjxw555513ZOrUqcHzkUhERGRQq6elpWVQ6wgAgLgqIWOMPP7447J9+3Z5++23pbS01MkvLS2VSCQitbW1wXO9vb1SV1cnc+bMGZkSAwBGjbi645YuXSrbtm2T3/3ud5KXlxe0eMLhsGRnZ0taWppUVlZKdXW1lJeXS3l5uVRXV0tOTo4sXrz4gnwAwBdjTNR0PK8dqX2BVBNXJbRp0yYREVmwYIHz/JYtW+Thhx8WEZEVK1ZId3e3LFmyRFpbW+Xmm2+WnTt3Sl5e3ogUGAAwesRVCZ3Lf2RpaWlSVVUlVVVViZYJAHCJYNgNAMAblnIAEqRvLcjOznbSdhd0bm6ukxctrbuu9XG5pQGjCS0hAIA3VEIAAG+ohAAA3hATAhKkR4sODAw4aXv5b70UuH5ttH31cblvCKMJLSEAgDdUQgAAb+iOAxIUDoeddF9fn5Pu7u4Otnt7e508Pcx6ypQpwfa0adOcPD35r35fIJXREgIAeEMlBADwhkoIAOANMSFghERbETWelVWjrdA6VBpIZbSEAADeUAkBALyhEgIAeENMCEhQe3u7k25paXHShw4dCrYbGhqcvKamJidt30c0duxYJy8rK8tJZ2T872c7YcKEOEoMJB9aQgAAb6iEAADeUAkBALwhJgQk6PTp0066p6fHSZ86dWrI7VhpnaePq98XSGW0hAAA3lAJAQC8oTsOSJBeAVWn7W4zvcxDtLTubov1PkAqoyUEAPCGSggA4A2VEADAG2JCQIL0kt16aHVXV1ewbS/1PdS+dr79uqH21e8LpDJaQgAAb6iEAADeUAkBALwhJgQkyF5SQUQkMzPTSdtLMuh9o6X1Ug76uOnp6fEXFkhStIQAAN5QCQEAvKE7DkhQtGHWIiKdnZ3D5kVL268bal89qzaQymgJAQC8oRICAHhDJQQA8IaYEJCgrKwsJ52dne2kx40bF2zn5OQ4edHS9uuG2le/L5DKaAkBALyhEgIAeEMlBADwhpgQkKBY9wmdPHky2NbLM0RL268bal/uE8JoQksIAOANlRAAwBu644AE5eXlOenx48c76csuuyzY/vLLL5083aU2ceLEYHvy5MlRj6uHcAOpjJYQAMAbKiEAgDdUQgAAb4gJAQmKNUS7o6Mj2NbLM+hh13a+HqLNUg4YzWgJAQC8oRICAHhDJQQA8IaYEJCgSZMmOenTp087aTvOo+M4mZmZTrq4uDjYLisrc/IikYiTtu8pAlIdLSEAgDdUQgAAb6iEAADeEBMCEjR27FgnrZfhtud4y8/Pd/L0vUB2vp4bTi8brt8XSGW0hAAA3lAJAQC8oTsOSFBra6uTbmpqctL19fVDbouIHD161EnbUwDp4du6+23MmP/971hYWBhHiYHkQ0sIAOANlRAAwJvzqoRqamokLS1NKisrg+eMMVJVVSXFxcWSnZ0tCxYskP37959vOQEAo1DCMaHdu3fL5s2b5brrrnOeX7dunaxfv15+/vOfyxVXXCHPP/+83H777XLw4MFByyEDqcxeqkFEpK2tzUkfP3482D5x4oSTp+NJWVlZQ75uqOOGw+Fgm5gQUl1CLaGTJ0/KfffdJ6+99ppMmDAheN4YIxs2bJDVq1fLokWLZObMmbJ161bp6uqSbdu2DXmsnp4eaW9vdx4AgEtDQpXQ0qVLZeHChXLbbbc5zzc0NEhzc7NUVFQEz4VCIZk/f77s2rVryGPV1NRIOBwOHiUlJYkUCQCQguKuhN544w358MMPpaamZlBec3OziAzuIigsLAzytFWrVklbW1vwaGxsjLdIAIAUFVdMqLGxUZYtWyY7d+50+rC1tLQ0J22MGfTcWaFQSEKhUDzFAJJCenq6k7bv39H5et+MjIxz3jdWGkhlcbWE9uzZIy0tLTJr1izJyMiQjIwMqaurk5deekkyMjKCFpBu9bS0tBBABQAMElcldOutt8q+fftk7969wWP27Nly3333yd69e6WsrEwikYjU1tYGr+nt7ZW6ujqZM2fOiBceAJDa4uqOy8vLk5kzZzrP5ebmyqRJk4LnKysrpbq6WsrLy6W8vFyqq6slJydHFi9ePHKlBpKAHjqth2HbQ62//PJLJ++LL75w0vZUPfaIU30cEXeItjHGyRuu2xtIViM+d9yKFSuku7tblixZIq2trXLzzTfLzp07uUcIADDIeVdC7777rpNOS0uTqqoqqaqqOt9DAwBGOeaOAwB4w1IOQIL0iqd6RdTc3NxgW6+6qldatfPt14nIoK5s+32JASHV0RICAHhDJQQA8IZKCADgDTEh4Bz19/c76a+++spJ63t/oi3loJf3tpfw1vcFtbS0OGn7PqHTp087eXppcCDZ0RICAHhDJQQA8IbuOCAKuwtOz16th07r6XbsbjO976RJk4Y9lv26odL2UHA9G7fuMmTGbSQ7WkIAAG+ohAAA3lAJAQC8SdqYUH9/f9C/bfdr66nrOzo6nLQ9RPXkyZNOnl7Btbu720nbq8X29vZGfa1N98vrKVrsMre3tzt5etqVU6dODZuv+/s1e2VPfZ50+e18HXPQdJntYcG6TPq82edGfza9Oq+dr2Mo+rU6pmIPcW5tbXXy9DViX096OQY9FY89DFuXSQ+d1os5RlvKQQ/Rtr8fHVvS72OXQ6/mqj/r+PHjnbR9zevvXX+X9vvoMunv2f58+nvVZYr2veu8vr4+J21/zzrepYen66mRbPp717/Dnp6eYfMGBgaGPa4uk/4d2p9P5+ny6vft6uoKtvV50ezy2+dMZPBvyb7m7deJDJ6KSt8SYF8jiU4hRUsIAOANlRAAwBsqIQCAN0kbEzpx4kTQ/2j3lXZ2djr7HTt2zEnbfa56WhUdq9F91Xa+7hvVfdU2HW+ZPHmyk7b7z3V5dZ++jmPZ+bo/VvfBRosJ6ViH3a+t+5f1cXVMwj43Oo6g+5vtfnodg9NlsvOjxTJEBp8L+7trampy8nT/v32e9HQ6ug/cnopHx0V0DEjHjOzPoD/PZZddds776qUdPvvss2Bb/x50PKywsNBJ29fXxIkTnTwd67BjRvp7tuMTIu55ihZXGyrf/m51nv6e7d+hjr/o36F9Xevfg55iSf8O7TLFignZx9ZxKb2vfX3pvGjxLxE3xqj/Pmn271CfF/3d2detztPXov5929exfZ3q/aKhJQQA8IZKCADgDZUQAMCbpI0JHT58OOg/te8J0f3LR44ccdJ2n2usPnvdf27HnnTfqI4V2HSsSfeH2v3LjY2NTp7u144Wv9D3Zui+6mjj9PVnt/v4dbxF948fOnTISdufT/fZ6/Nm92vr+IU+p3a8QsdMdPwu2rHq6+udPL00gn3OdZ6Ov9j5+r4a/Vrdf27HpnQsUMet7H57XQbN/g3o2Ia+fvT72OdNxy513MfO1zEIHbu0z7m+1vRnj/a968+ur3n7t6bjLzpua1/XOv5y+PBhJ61/h/p6s+n4kp3WcZzzibPp+Jh9b5neV7PPqf77pO/7s2Od+nNHIpFhj6vLYV//0c6fRksIAOANlRAAwJuk7Y6rr68PmpF2s1p3P3z++edO2m7O66HFuhtAd13Y+bq5q7sJ7K4v3WzWXQh2E1Z3FeluAN1UjmeIdrTuOF1+uysgVjdfQ0ODk7a74/SwUt21Z0/bo/N0N4F9zvU0N7r7TXcP2V2p+hzrYdj2Odfvo8tkv1Z3M+nX6u44uxtKf/ZoXbb6mtaf3e6OizVFkc63z7HO091D9ufT156+Tu0uW32t6fMf7XvXvyV9fdn5eros3R1nX9e6Cy3W71Cfc1u07rhYQ7Tt6zTWlFf675U9NF9fT7pMdvn1+dZdanYXoc7T14g+L/Z1bN8OEO38abSEAADeUAkBALyhEgIAeJO0MaGmpqYhp8rRw671ENSLFROy6T5XPU2G3c+qy3uhYkK6j1iX387X5dXH1WWOFhPSsY5EY0I65qO/D933bp9jvUyC7te2y6RjjHo6fftYevkI/T56WLl93vR50bcW2PR1r8+FHWvSeTpOpeMM9vWlr6do36Uuk75O7d+l/l71+T+fmJBdjlgxIXu4tP496L8j+pq3Yxo673xiQvZxdUxID+/Wf6/s6ynWEO14YkL2sXSeLqOO9djfgf1dxSqfjZYQAMAbKiEAgDdUQgAAb5I2JlRUVBT0ZU6ZMiV4XvcD6z5jO/ah4y26j1X3lY7UfUJTp0510nY/q+6j12XU73Mx7hOyz+9Qx9Flvhj3CellBnRfdHFxsZO2YzmxlhG3z7nuw9dlsvvp9X1COpam7xOy8/Vn16ZPnx5s6+tJx5rs61THqXQZ9Xm0z5N+re7/t1+rz3e0ZZ/1tabP08W6T8i+rnUcR8fSkuE+If13Q59ju8wX6z6hgoKCYY+ry8x9QgCAlEMlBADwJmm748rKyoacRVs3UXUT3M7XK2Gezyza0brjdHNXd13YTWc9DDMZZtGeNm3asO85VJkvxizaeoZn3U1QVlY27LF0l0i0WbT1NRLPLNr6tdFWRNVT8ejr1u5muvzyy528oqIiJ21P26PPk75+dH6is2jr862/D7vLTX+v+vxfrFm07S5O3S0W63foYxbtkpISJ627Ju1j699ZPN1xIzmLtt3lySzaAICUQyUEAPCGSggA4E3SxoRKSkqCOIYdq9HTquj+TrtfWPfRxxqqaPfB6iGcun/WpvuBdV+73a+tj6PjLzpuEm2Itmbvq/uI9fvafdX20EqRwbGlaNPI6H5tPTza7rePZ4i2/u70a3Vftf359Peh4yTRVlbVsTN76icdE9LTQkWb2inWUhT2MGwdA7riiiuGPZYeQqs/q863++pjxYTsz6vLFO271NeaXg05niHa+pq3r0Udx9FDwe3Prn8P+prWv0P78+nfg47z2GIN0bavEZ2nh+Lr69j+fDpWptnl1+dFf3fRbk3RcU/9+7avIfuzxVod2EZLCADgDZUQAMAbKiEAgDdJGxOaOHFiEP+xp+fQ/Zu6D9ZOR4sXiQzu/7TzdZ+rfl+bnj5E94fafb9632h90TpfT2ES7b4gTZff7iPXsQ59XF1mu58+1hLFdr+9joNEizVFm7JEZPB0NHb/uT6n+vPZZdJ5OiZh5+sy6fiRvfyyiBsL0XERezkGfWy9r44n2ffs6JiPjnNGi63pz66/S7tM+nzr79n+3elrTb+PzrfjDDov2nIHOiak/xbY76tjQvqa1te8vt6isY+tj6vjPvb1FW1KH5HB17H9efXfAs0uf7T7/ETcv5M6L9r0TCLubyCe+LWNlhAAwBsqIQCAN0nbHZeRkTGoaSsSu6vFprvjNN38vVDsMuvuEe1ilSkeephmstNDznXaprtpdLeM3a0QqytVd5vZq7bqKaJ0d489NZUeOq273Oz0VVdd5eTpbhrdRWV/3ni6czV93GhDcnX5k4HuXkwF8Qx7TpSPvz+0hAAA3lAJAQC8oRICAHiTtDEh4GKIFRexYx96uLCerl7HfewpdPS+X375pZO2h7rqqXd02h6GHSsGpJ1PHAi4EGgJAQC8oRICAHhDJQQA8IaYEHCO9DQx+p4vfT9SU1NTsK3jRXraG/u+lVj3CdnvGysGBCQ7WkIAAG+ohAAA3lAJAQC8ISYEJEgvu6GXUbDnktPLhrS3tztpe84uPQedvsfIft9Yc98ByY6WEADAGyohAIA3dMcBCdLLiOjlAeyh1nqItl55MtpSDjptvy/db0h1tIQAAN5QCQEAvIm7Ejpy5Ijcf//9MmnSJMnJyZHrr79e9uzZE+QbY6SqqkqKi4slOztbFixYIPv37x/RQgMARoe4KqHW1laZO3euZGZmyp/+9Cc5cOCA/PjHP3amEVm3bp2sX79eNm7cKLt375ZIJCK33377oGGmQKrr7+93HgMDA85D59uPvr4+5xFt32gPINXFNTDhxRdflJKSEtmyZUvw3OWXXx5sG2Nkw4YNsnr1alm0aJGIiGzdulUKCwtl27Zt8uijjw46Zk9Pj/T09ARpff8EAGD0iqsltGPHDpk9e7bcddddUlBQIDfccIO89tprQX5DQ4M0NzdLRUVF8FwoFJL58+fLrl27hjxmTU2NhMPh4FFSUpLgRwEApJq4KqH6+nrZtGmTlJeXy1/+8hd57LHH5Mknn5Rf/OIXIiLS3NwsIoNnEy4sLAzytFWrVklbW1vwaGxsTORzAABSUFzdcQMDAzJ79myprq4WEZEbbrhB9u/fL5s2bZIHH3ww2E/fu2CMGfZ+hlAoJKFQKN5yA97ZS3KLDL5vyL6/Ry/nferUKScdbSkHfdxx48bFX1ggScXVEioqKpKrr77aee6qq66Sw4cPi4hIJBIRERnU6mlpaRnUOgIAIK5KaO7cuXLw4EHnuU8//VSmT58uIiKlpaUSiUSktrY2yO/t7ZW6ujqZM2fOCBQXADCaxNUd99RTT8mcOXOkurpavvOd78gHH3wgmzdvls2bN4vImW64yspKqa6ulvLycikvL5fq6mrJycmRxYsXX5APAPgyYcIEJ62HTNujPvU0PVlZWU66uLg42C4rK3PyioqKnLSeHghIZXFVQjfeeKP89re/lVWrVsmaNWuktLRUNmzYIPfdd1+wz4oVK6S7u1uWLFkira2tcvPNN8vOnTsH9Z8DAJBm9IIknrW3t0s4HJa2tjZnjRUg2fT19Tnp48ePO+nPP/882Nbd2EePHnXSdkvoyiuvdPLOdnefZQ9cyMzMPPcCAxdJPH/HmTsOAOANSzkACert7XXSevVUe6VVPROITtvDrvUKrXoFV/t9aQkh1dESAgB4QyUEAPCGSggA4A0xISBBeiqeL774wknX19cH24cOHXLy9Kwi9n1EehqrnJwcJ23HgXJzc+MoMZB8aAkBALyhEgIAeEN3HJAgPfVOdna2k7ZnCdHdZrqLzc7Xs2Tr4zLrPEYTWkIAAG+ohAAA3lAJAQC8ISYEJKijo8NJf/XVV0762LFjwfaJEyecvNbWVidtx3n0RKj6uPZKqwUFBedcXiAZ0RICAHhDJQQA8IZKCADgDTEhIEGx7hOy7/fR9wVFS+v7hPS+3CeE0YSWEADAGyohAIA3VEIAAG+ICQEJOnXqlJPWy3Dby3Trpb+jpfXy3npf/b5AKqMlBADwhkoIAOAN3XFAgmIN0baXZ9B50dJ62QeWcsBoRksIAOANlRAAwBsqIQCAN8SEgAT19fU56dOnTzvp3t7eYfeNlrZfN9Rx+/v74y8skKRoCQEAvKESAgB4QyUEAPCGmBCQoLFjxzppfd+QvQSDvtcn2j1GeukGva9+XyCV0RICAHhDJQQA8IbuOCBB6enpUdOZmZnBdkaG+1OLlrZfdy7vA6QyWkIAAG+ohAAA3lAJAQC8ISYEJEjHbvQSC/bQaj3MOlpa5+nj6vcFUhktIQCAN1RCAABvqIQAAN4QEwISlJ+f76T18gzd3d3Btl6eQU+9M2XKlGB7+vTpTl5BQYGTDofD8RcWSFK0hAAA3lAJAQC8oTsOGCHGGCc9MDAwbF60tP26c3ktkMpoCQEAvKESAgB4QyUEAPCGmBCQoLa2Nif9xRdfOOnDhw8H2/X19U5ec3Ozk+7p6Qm2Y63Yai/7MHHixDhKDCQfWkIAAG+ohAAA3lAJAQC8ISYEJCgtLc1Jjxnj/k9nL8Otl+TWr422rz6ufi2QymgJAQC8oRICAHhDJQQA8IaYEJAg+94eEXfpBhGRjo6OYLuzs9PJi5a2XzfUcfX7AqmMlhAAwBsqIQCAN3THAQnSQ6XjGTp9ofYFUg0tIQCAN1RCAABv4qqE+vr65Nlnn5XS0lLJzs6WsrIyWbNmzaAVJKuqqqS4uFiys7NlwYIFsn///hEvOAAg9cVVCb344ovy6quvysaNG+WTTz6RdevWyY9+9CN5+eWXg33WrVsn69evl40bN8ru3bslEonI7bffPmjYKZDqBgYGoj6MMSPy0McFRpO4Bib87W9/k29/+9uycOFCERG5/PLL5Ve/+pX885//FJEzraANGzbI6tWrZdGiRSIisnXrViksLJRt27bJo48+OuiYPT09zn0P7e3tCX8YAEBqiaslNG/ePHnrrbfk008/FRGRjz76SN5//3254447RESkoaFBmpubpaKiInhNKBSS+fPny65du4Y8Zk1NjYTD4eBRUlKS6GcBAKSYuFpCzzzzjLS1tcmMGTMkPT1d+vv7Ze3atXLvvfeKyP9WiywsLHReV1hYKIcOHRrymKtWrZLly5cH6fb2dioiALhExFUJvfnmm/L666/Ltm3b5JprrpG9e/dKZWWlFBcXy0MPPRTsp+9rMMYMe69DKBSSUCiUQNGB5GKMGTYdLS/WvrHeB0hlcVVCTz/9tKxcuVLuueceERG59tpr5dChQ1JTUyMPPfSQRCIRETnTIioqKgpe19LSMqh1BABAXDGhrq6uIRfuOjtip7S0VCKRiNTW1gb5vb29UldXJ3PmzBmB4gIARpO4WkLf+ta3ZO3atTJt2jS55ppr5F//+pesX79eHnnkERE50w1XWVkp1dXVUl5eLuXl5VJdXS05OTmyePHiC/IBAF/0CqiZmZlOeuzYscPmRUvbrxtq34wMZtvC6BHX1fzyyy/L97//fVmyZIm0tLRIcXGxPProo/KDH/wg2GfFihXS3d0tS5YskdbWVrn55ptl586dkpeXN+KFBwCktjSTZFHO9vZ2CYfD0tbWJvn5+b6LAwzr2LFjTvrIkSNO+pNPPgm2Dx486OQ1NTU5aTuGeuWVVzp5V111lZOeMmVKsH3ZZZfFUWLg4ojn7zhzxwEAvKFzGUiQ7kTQU+rY6Wh557MvkOpoCQEAvKESAgB4QyUEAPCGmBCQoKysLCc9btw4Jx0Oh4NtPUKos7PTSdv59utERHJzc510dnZ2/IUFkhQtIQCAN1RCAABv6I4DEhRriHZ/f/+Q27HSOo8h2hjNaAkBALyhEgIAeEMlBADwhpgQkCA97Lqvr89Jd3V1Bdu9vb1Onl6uYerUqcG2Xt5eT1LKxL4YTWgJAQC8oRICAHhDJQQA8IaYEJCgtLQ0J63jPPb0OnpKn2hpPS2PPu6YMfzviNGDqxkA4A2VEADAGyohAIA3xISABNn3AYmIHD9+3EkfPXo02P7ss8+cvKamJifd09MTbOsYUCgUctLp6enBtl7mAUg1tIQAAN5QCQEAvKE7DkiQnqbH7lITcbvruru7h83T+TpPT/mj3xdIZbSEAADeUAkBALyhEgIAeENMCLgI9FLg0fJj7QuMJrSEAADeUAkBALyhEgIAeENMCEiQvi/o1KlTTvrkyZPBdmdnp5MXLa3z9D1G+n2BVEZLCADgDZUQAMAbuuOABMWa7TonJ2fI7VhpnaePq98XSGW0hAAA3lAJAQC8oRICAHhDTAhIkB6SrZdgaG9vD7bt4dqx0vbrRAYP0dbvC6QyWkIAAG+ohAAA3lAJAQC8ISYEJCg7O9tJ5+XlOekJEyYE2+Fw2MnT8SM7336diMi4ceOctL6PCEhltIQAAN5QCQEAvKE7DkiQns1ad7F1dHQE2/EM0bZfN9Rx7SHa+fn5cZQYSD60hAAA3lAJAQC8oRICAHhDTAhIkB46PX78eCcdiUSC7dbWVievv7/fSRcUFAz5OpHYQ7aBVEZLCADgDZUQAMAbKiEAgDfEhIAEnT592knrJRbse386OzudvGhpfQ+RPq5+XyCV0RICAHhDJQQA8IZKCADgDTEhIEF66YbLLrvMSdtzvum4TkaG+9ObOnVqsD1lyhQnb/LkyVHfF0hltIQAAN5QCQEAvKE7DkhQWlqak87NzXXS9jQ+hYWFTl5fX5+TtvP19D/6uGPG8L8jRg+uZgCAN1RCAABvkq47zhgjIiLt7e2eSwLER4+As1dI1TMk6NVS7Xy9sqrujmPGBCS7s3+/z/49jybpKqGzP8CSkhLPJQEAnI+Ojg4Jh8NR90kz51JVXUQDAwNy9OhRMcbItGnTpLGxUfLz830XK2m1t7dLSUkJ5ykGztO54TydG85TdMYY6ejokOLi4pgDaZKuJTRmzBiZOnVq0JzLz8/nSz4HnKdzw3k6N5ync8N5Gl6sFtBZDEwAAHhDJQQA8CZpK6FQKCTPPfechEIh30VJapync8N5Ojecp3PDeRo5STcwAQBw6UjalhAAYPSjEgIAeEMlBADwhkoIAOANlRAAwJukrYReeeUVKS0tlaysLJk1a5a89957vovkTU1Njdx4442Sl5cnBQUFcuedd8rBgwedfYwxUlVVJcXFxZKdnS0LFiyQ/fv3eypxcqipqZG0tDSprKwMnuM8nXHkyBG5//77ZdKkSZKTkyPXX3+97NmzJ8jnPJ1Z8+nZZ5+V0tJSyc7OlrKyMlmzZo0MDAwE+3CeRoBJQm+88YbJzMw0r732mjlw4IBZtmyZyc3NNYcOHfJdNC++8Y1vmC1btph///vfZu/evWbhwoVm2rRp5uTJk8E+L7zwgsnLyzO/+c1vzL59+8zdd99tioqKTHt7u8eS+/PBBx+Yyy+/3Fx33XVm2bJlwfOcJ2NOnDhhpk+fbh5++GHzj3/8wzQ0NJi//vWv5r///W+wD+fJmOeff95MmjTJ/OEPfzANDQ3m17/+tRk3bpzZsGFDsA/n6fwlZSV00003mccee8x5bsaMGWblypWeSpRcWlpajIiYuro6Y4wxAwMDJhKJmBdeeCHY59SpUyYcDptXX33VVzG96ejoMOXl5aa2ttbMnz8/qIQ4T2c888wzZt68ecPmc57OWLhwoXnkkUec5xYtWmTuv/9+YwznaaQkXXdcb2+v7NmzRyoqKpznKyoqZNeuXZ5KlVza2tpERGTixIkiItLQ0CDNzc3OOQuFQjJ//vxL8pwtXbpUFi5cKLfddpvzPOfpjB07dsjs2bPlrrvukoKCArnhhhvktddeC/I5T2fMmzdP3nrrLfn0009FROSjjz6S999/X+644w4R4TyNlKSbRfv48ePS398vhYWFzvOFhYXS3NzsqVTJwxgjy5cvl3nz5snMmTNFRILzMtQ5O3To0EUvo09vvPGGfPjhh7J79+5BeZynM+rr62XTpk2yfPly+d73vicffPCBPPnkkxIKheTBBx/kPP2/Z555Rtra2mTGjBmSnp4u/f39snbtWrn33ntFhOtppCRdJXRWWlqakzbGDHruUvT444/Lxx9/LO+///6gvEv9nDU2NsqyZctk586dkpWVNex+l/p5GhgYkNmzZ0t1dbWIiNxwww2yf/9+2bRpkzz44IPBfpf6eXrzzTfl9ddfl23btsk111wje/fulcrKSikuLpaHHnoo2O9SP0/nK+m64yZPnizp6emDWj0tLS2D/uO41DzxxBOyY8cOeeedd2Tq1KnB85FIRETkkj9ne/bskZaWFpk1a5ZkZGRIRkaG1NXVyUsvvSQZGRnBubjUz1NRUZFcffXVznNXXXWVHD58WES4ns56+umnZeXKlXLPPffItddeKw888IA89dRTUlNTIyKcp5GSdJXQ2LFjZdasWVJbW+s8X1tbK3PmzPFUKr+MMfL444/L9u3b5e2335bS0lInv7S0VCKRiHPOent7pa6u7pI6Z7feeqvs27dP9u7dGzxmz54t9913n+zdu1fKyso4TyIyd+7cQUP8P/30U5k+fbqIcD2d1dXVNWhV0PT09GCINudphHgcFDGss0O0f/azn5kDBw6YyspKk5ubaz7//HPfRfPiu9/9rgmHw+bdd981TU1NwaOrqyvY54UXXjDhcNhs377d7Nu3z9x7770MFTXGGR1nDOfJmDPD1zMyMszatWvNf/7zH/PLX/7S5OTkmNdffz3Yh/NkzEMPPWSmTJkSDNHevn27mTx5slmxYkWwD+fp/CVlJWSMMT/5yU/M9OnTzdixY83Xvva1YDjypUhEhnxs2bIl2GdgYMA899xzJhKJmFAoZG655Razb98+f4VOEroS4jyd8fvf/97MnDnThEIhM2PGDLN582Ynn/NkTHt7u1m2bJmZNm2aycrKMmVlZWb16tWmp6cn2IfzdP5YTwgA4E3SxYQAAJcOKiEAgDdUQgAAb6iEAADeUAkBALyhEgIAeEMlBADwhkoIAOANlRAAwBsqIQCAN1RCAABv/g/U9pXugFB2sgAAAABJRU5ErkJggg==" + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plt.imshow(np.log10(psf_185_65_lpi), cmap='Greys', vmin=-12, vmax=0)" + ], + "metadata": { + "collapsed": false, + "ExecuteTime": { + "end_time": "2024-05-09T21:04:58.707355Z", + "start_time": "2024-05-09T21:04:58.649644Z" + } + }, + "id": "b319b5605ee0debf" + }, + { + "cell_type": "code", + "execution_count": 5, + "outputs": [ + { + "data": { + "text/plain": "" + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAGgCAYAAAAD9NhnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAj1klEQVR4nO3df2xV9f3H8delpbcttBWKtFQBi6lBBSMDZQMibArLJMsMiVNBxfgPDFAqiQjDzY5oqywhRJ04yIIsjmmWsY0t+0Hnj0bCnAyHMjCwIYMGbBqwtgVKS9vP9w++HM/59HJuL5zuc3t5PpIbzud8zvn0c0/P7ZvP533OuTFjjBEAAA4McN0BAMCViyAEAHCGIAQAcIYgBABwhiAEAHCGIAQAcIYgBABwhiAEAHCGIAQAcIYgBABwps+C0CuvvKLy8nLl5uZq4sSJeu+99/rqRwEA+qnsvmj0zTffVGVlpV555RVNnTpVP/3pT/Wtb31L+/fv16hRo0L37e7u1vHjx1VQUKBYLNYX3QMA9CFjjFpbW1VWVqYBA5KMdUwfuP32283ChQsD68aOHWtWrFiRdN/6+nojiRcvXrx49fNXfX190r/5kY+EOjo6tHv3bq1YsSKwftasWdq5c2eP7dvb29Xe3u6Vzf8/1Lu+vl6FhYVRdw/4n9m7d6+3fOjQoUDd8ePHA+WysjJv+frrrw/UjR8/vg96B/SdlpYWjRw5UgUFBUm3jTwInThxQl1dXSopKQmsLykpUUNDQ4/ta2pq9KMf/ajH+sLCQoIQ+rXBgwd7y/n5+YG63NzcQNlf799PEp8D9Fu9San0SU4o0Q83xiTs0MqVK7Vs2TKvfCGCAv2dP5gMHTo0UHfmzJlA2V9vByEgk0UehIYNG6asrKweo57GxsYeoyNJisfjisfjUXcDANAPRH6Jdk5OjiZOnKja2trA+traWk2ZMiXqHwcA6Mf6ZDpu2bJleuihhzRp0iR97Wtf04YNG3T06FEtXLiwL34ckJZOnTrlLTc1NQXqPv/880B50KBB3nJxcXHfdgxII30ShO677z6dPHlSq1ev1meffaZx48bpj3/8o0aPHt0XPw4A0E/12YUJixYt0qJFi/qqeQBABuDZcQAAZ/psJARcabq7uwNl/2XYzc3NgTp/vsiuty/ftttN+hgUoB/hbAYAOEMQAgA4QxACADhDTgiISFdXV6Dc0dHhLZ89ezZQ539or13v3y9Ru+SEkEk4mwEAzhCEAADOMB0HRMS+lLqzszPhsiSdO3eu19va7QKZhJEQAMAZghAAwBmCEADAGXJCQESMMYGyP5dj53nsy6799XYOyG4XyCSMhAAAzhCEAADOEIQAAM6QEwIiYudu/Hkfuy6sbOeLyAkhkzESAgA4QxACADhDEAIAOEMQAgA4QxACADhDEAIAOMMl2kBEYrHYRcthdaluC2QSRkIAAGcIQgAAZwhCAABnyAkBEbFzN9nZ2QmXJSkrK6vX25ITQiZjJAQAcIYgBABwhiAEAHCGnBAQkQEDgv+nGzhwYMLlZGW7zm4XyCSc3QAAZwhCAABnmI4DImJPo8XjcW85Pz8/UBdW9u+XqF0gkzASAgA4QxACADhDEAIAOENOCOgjhYWF3nJxcXGgrq2tLVD21/v3AzIdIyEAgDMEIQCAMwQhAIAz5ISAPtLe3u4t2zmgsLJ/PyDTMRICADhDEAIAOMN0HBCRrq6uQPnUqVPe8hdffBGoa25uDpQHDx6ccL9E7drfygr0Z4yEAADOEIQAAM4QhAAAzpATAiLS2dkZKEd1ibbdLjkhZBJGQgAAZwhCAABnCEIAAGfICQER6e7uDpQ7OjoSLks98zxh29rtApmEkRAAwBmCEADAGYIQAMAZckJAROzcjT/vYz//zc4J+evtOnJCyGSMhAAAzhCEAADOMB0HRMQYEyj7p9HsKbXL2RbIJIyEAADOEIQAAM6kFIRqamp02223qaCgQMOHD9c999yjAwcOBLYxxqiqqkplZWXKy8vTjBkztG/fvkg7DQDIDCkFobq6Oi1evFjvv/++amtr1dnZqVmzZun06dPeNmvWrNHatWv18ssva9euXSotLdXMmTPV2toaeeeB/soYE3gBV6qULkz485//HChv2rRJw4cP1+7du3XHHXfIGKN169Zp1apVmjNnjiRp8+bNKikp0ZYtW7RgwYIebba3twe+P6WlpeVS3gcAoB+6rJxQc3OzJGno0KGSpMOHD6uhoUGzZs3ytonH45o+fbp27tyZsI2amhoVFRV5r5EjR15OlwAA/cglByFjjJYtW6Zp06Zp3LhxkqSGhgZJUklJSWDbkpISr862cuVKNTc3e6/6+vpL7RIAoJ+55PuElixZoo8//lg7duzoUReLxQJlY0yPdRfE43HF4/FL7QaQNuxz/GLnfKp1YdsC/d0ljYQee+wxbdu2Te+8846uvfZab31paakk9Rj1NDY29hgdAQCQUhAyxmjJkiXaunWr3n77bZWXlwfqy8vLVVpaqtraWm9dR0eH6urqNGXKlGh6DADIGClNxy1evFhbtmzR7373OxUUFHgjnqKiIuXl5SkWi6myslLV1dWqqKhQRUWFqqurlZ+fr7lz5/bJGwDShT1tNmDAl//Hy8rKumidXW/XMR2HTJZSEFq/fr0kacaMGYH1mzZt0iOPPCJJWr58udra2rRo0SI1NTVp8uTJ2r59uwoKCiLpMAAgc6QUhHpzU10sFlNVVZWqqqoutU8AgCsEz44DADjDVzkAEbFzOQMHDvSWs7ODH7Wwsn+/RO0CmYSzGwDgDEEIAOAMQQgA4Aw5ISAi9r1Aubm5CZclKS8vr9fb2u0CmYSREADAGYIQAMAZpuOAiNiXVg8ePNhbHjJkSKDO/23Edr1/v0TtApmEkRAAwBmCEADAGYIQAMAZghAAwBmCEADAGYIQAMAZghAAwBnuEwIi0tHRESi3tLR4yydPngzUffHFF4Gy/zE+/v0StZuTk3M53QTSCiMhAIAzBCEAgDMEIQCAM+SEgIh0dnYGyv7nw505cyZQF1a2nytnt0tOCJmEkRAAwBmCEADAGabjgIh0dXUFyu3t7QmXpZ6XXYdta7cLZBJGQgAAZwhCAABnCEIAAGfICQER6e7uDpT9l1bbeZ1z584Fyv56+5Jsu10gkzASAgA4QxACADhDEAIAOENOCIiIMSZQ9ud57JyQnecJ29ZuF8gkjIQAAM4QhAAAzjAdB0TEnjbzl8PqUt0WyCSMhAAAzhCEAADOEIQAAM4QhAAAzhCEAADOEIQAAM4QhAAAznCfEBCRWCwWKA8YMOCidWFl/36JtgUyCSMhAIAzBCEAgDNMxwERsafRsrO//HhlZWUF6sLK/v0StQtkEs5uAIAzBCEAgDMEIQCAM+SEgIjYuZucnBxvOR6PB+rCyv79ErULZBLObgCAMwQhAIAzBCEAgDPkhICI2LmcwsJCb7mgoCBQ19LSEij76/37JWoXyCSMhAAAzhCEAADOEIQAAM6QEwIiYj8Pzp/Lyc3NDdSFle0ckN0ukEkYCQEAnCEIAQCcYToOiMjZs2cD5ZMnT3rLJ06cuGidFHxsj11ntzto0KDL6ieQThgJAQCcIQgBAJy5rCBUU1OjWCymyspKb50xRlVVVSorK1NeXp5mzJihffv2XW4/AQAZ6JJzQrt27dKGDRt0yy23BNavWbNGa9eu1WuvvaYbbrhBzz77rGbOnKkDBw70eHQJkEk6OjoC5dbWVm/51KlTgbqwsn+/RO2SE0ImuaSR0KlTpzRv3jxt3LhRQ4YM8dYbY7Ru3TqtWrVKc+bM0bhx47R582adOXNGW7ZsSdhWe3u7WlpaAi8AwJXhkoLQ4sWLNXv2bN11112B9YcPH1ZDQ4NmzZrlrYvH45o+fbp27tyZsK2amhoVFRV5r5EjR15KlwAA/VDKQeiNN97Qhx9+qJqamh51DQ0NkqSSkpLA+pKSEq/OtnLlSjU3N3uv+vr6VLsEAOinUsoJ1dfXa+nSpdq+fXuPx474xWKxQNkY02PdBfF4vMdXHQP90blz5wLltra2hMvS+Wno3m5r54SATJLSSGj37t1qbGzUxIkTlZ2drezsbNXV1enFF19Udna2NwKyRz2NjY09RkcAAKQUhO68807t3btXe/bs8V6TJk3SvHnztGfPHo0ZM0alpaWqra319uno6FBdXZ2mTJkSeecBAP1bStNxBQUFGjduXGDdoEGDVFxc7K2vrKxUdXW1KioqVFFRoerqauXn52vu3LnR9RpIQ52dnYGyf8rNnqoLK9tTdV1dXVF1EUg7kT87bvny5Wpra9OiRYvU1NSkyZMna/v27dwjBADo4bKD0Lvvvhsox2IxVVVVqaqq6nKbBgBkOJ4dBwBwhq9yACJi5278eR47X2Rfdu2vt/NF5ISQyRgJAQCcIQgBAJwhCAEAnCEnBESku7s7UPbncuw8j50j8tfbOSC7XSCTMBICADhDEAIAOMN0HBARexrNXw6bqrPrw9oBMg0jIQCAMwQhAIAzBCEAgDPkhICIGGMCZX+eJ5WckL2t3S6QSRgJAQCcIQgBAJwhCAEAnOn3OaGwufZkc+thc+2xWCy0PGDAl/F74MCBveuseuYC7Me32H32l5PlBvx9DOuvJGVlZXnLyfof9rUDyR4x4+9zKv23+2uXc3JyAmX/+zl79mygzu5/2FcspHL/jv1+7J/jr08lJ5Ss3UOHDnnL9u/Zfxyk8N97dnbw42+fB/5jnJubG9p/fx/DzmEpus+d/V7t92OfI37Jvi4j7PeRSv/DzmO7v/b7CWP3P6rPXSp/NxLVXwpGQgAAZwhCAABnCEIAAGfSNie0d+9eDR48WJK8fyXp1KlTge3OnDkTKPvnppPlW8KeyWXPjYbNn8fj8UBdYWFhoNze3u4t2/3310nh+Zdkj/QPy6nY8/3+OX7/8U2kpaUlUD59+rS3bPc/LMeSbG46bL7cnt+3j7G//uTJk4G61tbWQLmtrc1btvsflitIJf8lBc/NsK/+tuvtc/q3v/1toBz2e7bn7O1y2Hmbl5cXKBcUFHjLxcXFgTr7PPWfI2HnsBR+HifLcfnPC7v/gwYNCpTtc8TP/hzaeUT/7yeVnJb9+wg7j+3+259Du95/jJOdt/5jnuxz5z/GqfRfkvLz871lf//tz1wYRkIAAGcIQgAAZ9J2Ou7QoUPeUG/o0KHe+qampsB2zc3NgbJ/WG1PA9jlVC61DJuO8w9JpZ5TF/7pny+++OKidVLPqYywS87DJLt81T8dN2TIkNC27Okt/3SRPS0Q9jUEqVwqavffnprwTxVJwfdz4sSJQJ099eI/5ql842myS4/DpnuPHj0aqDty5Ih6K+yy2GSXstvTsP7zwK6zp+P80yvDhg0L1NnTV/7pl2TnRCqfu7DpRfucSPY59LP/jtjvJ+w2hDDJpkf901n28b7qqqsCZbve/zm0p2wvZzrOf8ztvxNhfzckqaioyFv2/x3xT9knw0gIAOAMQQgA4AxBCADgTNrmhI4fP+7NP/rnPz///PPAdmGXPCd7NMflzE2H5YTsPI+/bOew7G3D8ljJHh8SdumuPbfrn29ONn9r57H8v4+wR+JIl/7YoWQ5Ifuycf9ctZ3DSuUcSeURRcnK/rbtHNCxY8fUW3buxn9skl2SHXaJrd1u2CXDdp7HzqGEXbZ/OZ+7sBxX2OXCUs/Plp99/oR9DpNdoh322JuwHIud87HPU7ven8dKlhOKKheb7BzxH0f/ORB27G2MhAAAzhCEAADOEIQAAM6kbU6orKws4X1C9qM5+tt9QvajOfrDfUL23HS63ydkb5sO9wklM3r0aG951KhRgTruE0pcvtLuE/KXuU8IAIAIEIQAAM6k7XTc9ddfn/Ap2vYQm6doJ+4zT9H+Ujo8RTsZ/xRceXl5oM6eZuIp2uddaU/Rvvrqq71lnqINAEAECEIAAGcIQgAAZ9I2JzR+/PjQOd0LwvI8yeZyU7lUNOxxIna+JYydV0g2X36pj71J9vgT/zxwsv6n8iibsGOeSv+TXXpsz0373489v2/3319O9ugd//tJdqmx/XP834ia7HEo/no7B3TPPfcEyv73niyHEvZ7D7t82C7bl+bax8L/3pNduh7V5y6VRxTZkj3SK+xzdzmPHQrLe9rvJ0zY47HsPvbV341E9RfYeeQwjIQAAM4QhAAAzhCEAADOpG1OqLeS5Q7STbL7ONJRWK4gHdn3i9jlvnLo0KFAOZX7efz19ry8fbyvv/76y+pnFJLlk9KdnZNLJa+bDvpbf8Ok919sAEBGIwgBAJzp99NxQLoIu7w1lek4e1u7XSCTMBICADhDEAIAOEMQAgA4Q04IiEjY5fep5IT642X8wKViJAQAcIYgBABwhiAEAHCGnBAQkVS+LiPs6+J7+7h8IBNwdgMAnCEIAQCcIQgBAJwhJwREJOwrvJN99bS/Puyrv4FMw0gIAOAMQQgA4AzTcUBE7Cm3eDzuLSf7Jk9/2b+fxHQcMhsjIQCAMwQhAIAzKQehY8eO6cEHH1RxcbHy8/N16623avfu3V69MUZVVVUqKytTXl6eZsyYoX379kXaaQBAZkgpJ9TU1KSpU6fq61//uv70pz9p+PDhOnTokK666ipvmzVr1mjt2rV67bXXdMMNN+jZZ5/VzJkzdeDAARUUFETdfyBt2HmevLy8hMtSz7xP2Lb25dxAJkkpCL3wwgsaOXKkNm3a5K277rrrvGVjjNatW6dVq1Zpzpw5kqTNmzerpKREW7Zs0YIFC3q02d7ervb2dq/c0tKS6nsAAPRTKU3Hbdu2TZMmTdK9996r4cOHa8KECdq4caNXf/jwYTU0NGjWrFneung8runTp2vnzp0J26ypqVFRUZH3Gjly5CW+FQBAf5NSEPr000+1fv16VVRU6C9/+YsWLlyoxx9/XD//+c8lSQ0NDZKkkpKSwH4lJSVenW3lypVqbm72XvX19ZfyPgAA/VBK03Hd3d2aNGmSqqurJUkTJkzQvn37tH79ej388MPedrFYLLCfMabHugvi8XiP+XGgP7JzN/4c6ODBgwN1YWU7d0pOCJkspZHQiBEjdNNNNwXW3XjjjTp69KgkqbS0VJJ6jHoaGxt7jI4AAEgpCE2dOlUHDhwIrDt48KBGjx4tSSovL1dpaalqa2u9+o6ODtXV1WnKlCkRdBcAkElSmo574oknNGXKFFVXV+u73/2uPvjgA23YsEEbNmyQdH4arrKyUtXV1aqoqFBFRYWqq6uVn5+vuXPn9skbANJFbm5uoFxcXOwtDxs2LFDnvyLUrvfvl6hdIJOkFIRuu+02/eY3v9HKlSu1evVqlZeXa926dZo3b563zfLly9XW1qZFixapqalJkydP1vbt27lHCADQQ8wYY1x3wq+lpUVFRUVqbm5WYWGh6+4AvdbV1RUov//++97yJ598Eqg7duxYoHzNNdd4yzfeeGOg7qtf/WqgzANNke5S+TvOs+MAAM7wVQ5AROyRUEdHh7d89uzZQF1Y2b9fonYZCSGTMBICADhDEAIAOEMQAgA4Q04IiIidy/E/Eb61tTVQd/r06UDZX28/Sd5ul8f4IJMwEgIAOEMQAgA4w3QcEJHu7u5A2T+NZj+mJ6xsT7/Z7QKZhJEQAMAZghAAwBmCEADAGXJCQETs3E1nZ6e3bD96J6zs3y9Ru0AmYSQEAHCGIAQAcIYgBABwhpwQEBH7+yH9uRy7Lqxs54DS7HsngUgxEgIAOEMQAgA4QxACADhDEAIAOEMQAgA4QxACADjDJdpARGKx2EXLYXWpbgtkEkZCAABnCEIAAGcIQgAAZ8gJARGxczdZWVkJlyVpwIABvd6WnBAyGSMhAIAzBCEAgDMEIQCAM+SEgIjYeZ7s7C8/XnaeZ+DAgYGyv96/X6J2gUzC2Q0AcIYgBABwhuk4ICL2lFs8Hk+4LEk5OTm93tZuF8gkjIQAAM4QhAAAzhCEAADOkBMCImJfWj1o0CBvOT8/P1AXVvbvl6hdIJMwEgIAOEMQAgA4QxACADjDZDMQEfven8LCQm+5uLg4UNfW1hYo++v9+yVqF8gkjIQAAM4QhAAAzhCEAADOEIQAAM4QhAAAzhCEAADOcIk2EJFz584FyqdOnfKWm5qaAnUtLS2Bsr/ev1+idu1vZQX6M0ZCAABnCEIAAGcIQgAAZ8gJARHp6uoKlM+ePZtwWer52J6wbe12yQkhkzASAgA4QxACADhDEAIAOENOCIhId3d3oOy/v6ezszNQF1a27wuy2wUyCSMhAIAzBCEAgDNMxwERMcYEyv5pNPsya3uKzV9v19ntApmEkRAAwBmCEADAmZSCUGdnp55++mmVl5crLy9PY8aM0erVqwPTB8YYVVVVqaysTHl5eZoxY4b27dsXeccBAP1fSkHohRde0KuvvqqXX35Zn3zyidasWaMf//jHeumll7xt1qxZo7Vr1+rll1/Wrl27VFpaqpkzZ6q1tTXyzgPpxBhz0VeybXvbDvkhZJqULkz429/+pu985zuaPXu2JOm6667TL3/5S/3jH/+QdP7Ds27dOq1atUpz5syRJG3evFklJSXasmWLFixY0KPN9vZ2tbe3e2X7e1YAAJkrpZHQtGnT9NZbb+ngwYOSpI8++kg7duzQ3XffLUk6fPiwGhoaNGvWLG+feDyu6dOna+fOnQnbrKmpUVFRkfcaOXLkpb4XAEA/k9JI6KmnnlJzc7PGjh2rrKwsdXV16bnnntMDDzwgSWpoaJAklZSUBPYrKSnRkSNHEra5cuVKLVu2zCu3tLQQiADgCpFSEHrzzTf1+uuva8uWLbr55pu1Z88eVVZWqqysTPPnz/e2i8Vigf2MMT3WXRCPxxWPxy+h60D/dbHPA3ClSSkIPfnkk1qxYoXuv/9+SdL48eN15MgR1dTUaP78+SotLZV0fkQ0YsQIb7/GxsYeoyMAAFLKCZ05c0YDBgR3ycrK8i7RLi8vV2lpqWpra736jo4O1dXVacqUKRF0FwCQSVIaCX3729/Wc889p1GjRunmm2/WP//5T61du1aPPvqopPNTDJWVlaqurlZFRYUqKipUXV2t/Px8zZ07t0/eAJAu7Ck2/3/Y7P+8Xc62QCZJKQi99NJL+sEPfqBFixapsbFRZWVlWrBggX74wx962yxfvlxtbW1atGiRmpqaNHnyZG3fvl0FBQWRdx4A0L/FTJrd/dbS0qKioiI1NzersLDQdXeAXjt9+nSg/Pbbb3vL+/fvD9QdP348UC4rK/OWb7rppkDdN77xjUB50KBBl9VPoK+l8necZ8cBAJzhqxyAiNi5nOzsLz9eWVlZF62z6+06u10gk3B2AwCcIQgBAJwhCAEAnCEnBETEzt3k5OQkXJZ65n3CtiUnhEzG2Q0AcIYgBABwhuk4ICL2FJv/6fB5eXmBurCy/VR5u10gkzASAgA4QxACADhDEAIAOMNkMxAR+9E8gwcP9pavuuqqQN2pU6cCZX+9f79E7QKZhJEQAMAZghAAwBmCEADAGXJCQB+J6j4hIJMxEgIAOEMQAgA4QxACADhDTgjoIy0tLd7yyZMnA3VNTU2Bsj8ndPXVV/dtx4A0wkgIAOAMQQgA4AzTcUBEzp07Fyi3t7d7y2fOnAnUhZX9+yVqd+DAgZfVTyCdMBICADhDEAIAOEMQAgA4Q04IiEh3d3eg7M/l2HmdsLJdZ7cLZBJGQgAAZwhCAABnCEIAAGfICQERMcYEyp2dnQmXJamrq6vX29rtApmEkRAAwBmCEADAGabjgIjY02b+clhdqtsCmYSREADAGYIQAMAZghAAwBmCEADAGYIQAMAZghAAwBmCEADAGe4TAiISi8UC5aysrIvWhZX9+yXaFsgkjIQAAM4QhAAAzjAdB0TEnjYbMODL/+NlZwc/avaUm7/ev1+idoFMwkgIAOAMQQgA4AxBCADgDDkhICJ2Lsef57FzQgMHDuz1tna7QCbh7AYAOEMQAgA4QxACADhDTgiIiH3vT05Ojrecm5sbqIvH44Gyv96/X6J2gUzCSAgA4AxBCADgDEEIAOAMOSEgIvb9PPn5+d5yUVFRoK6lpSVQ9tf790vULpBJOLsBAM4QhAAAzjAdB/SRwYMHe8tDhgwJ1J0+fTpQ9tf79wMyHSMhAIAzBCEAgDNpNx1njJHU8+ohoL9pbW31lu3pt7a2tkDZX+/fT+KzgP7nwjl74e95mLQLQhc+gCNHjnTcEwDA5Whtbe1xe4ItZnoTqv6Huru7dfz4cRljNGrUKNXX16uwsNB1t9JWS0uLRo4cyXFKguPUOxyn3uE4hTPGqLW1VWVlZUnvc0u7kdCAAQN07bXXesO5wsJCfsm9wHHqHY5T73CceofjdHHJRkAXcGECAMAZghAAwJm0DULxeFzPPPNMj+9dQRDHqXc4Tr3DceodjlN00u7CBADAlSNtR0IAgMxHEAIAOEMQAgA4QxACADhDEAIAOJO2QeiVV15ReXm5cnNzNXHiRL333nuuu+RMTU2NbrvtNhUUFGj48OG65557dODAgcA2xhhVVVWprKxMeXl5mjFjhvbt2+eox+mhpqZGsVhMlZWV3jqO03nHjh3Tgw8+qOLiYuXn5+vWW2/V7t27vXqOk9TZ2amnn35a5eXlysvL05gxY7R69Wp1d3d723CcImDS0BtvvGEGDhxoNm7caPbv32+WLl1qBg0aZI4cOeK6a05885vfNJs2bTL/+te/zJ49e8zs2bPNqFGjzKlTp7xtnn/+eVNQUGB+/etfm71795r77rvPjBgxwrS0tDjsuTsffPCBue6668wtt9xili5d6q3nOBnz+eefm9GjR5tHHnnE/P3vfzeHDx82f/3rX81//vMfbxuOkzHPPvusKS4uNn/4wx/M4cOHza9+9SszePBgs27dOm8bjtPlS8sgdPvtt5uFCxcG1o0dO9asWLHCUY/SS2Njo5Fk6urqjDHGdHd3m9LSUvP8889725w9e9YUFRWZV1991VU3nWltbTUVFRWmtrbWTJ8+3QtCHKfznnrqKTNt2rSL1nOczps9e7Z59NFHA+vmzJljHnzwQWMMxykqaTcd19HRod27d2vWrFmB9bNmzdLOnTsd9Sq9NDc3S5KGDh0qSTp8+LAaGhoCxywej2v69OlX5DFbvHixZs+erbvuuiuwnuN03rZt2zRp0iTde++9Gj58uCZMmKCNGzd69Ryn86ZNm6a33npLBw8elCR99NFH2rFjh+6++25JHKeopN1TtE+cOKGuri6VlJQE1peUlKihocFRr9KHMUbLli3TtGnTNG7cOEnyjkuiY3bkyJH/eR9deuONN/Thhx9q165dPeo4Tud9+umnWr9+vZYtW6bvf//7+uCDD/T4448rHo/r4Ycf5jj9v6eeekrNzc0aO3assrKy1NXVpeeee04PPPCAJM6nqKRdELogFosFysaYHuuuREuWLNHHH3+sHTt29Ki70o9ZfX29li5dqu3btys3N/ei213px6m7u1uTJk1SdXW1JGnChAnat2+f1q9fr4cfftjb7ko/Tm+++aZef/11bdmyRTfffLP27NmjyspKlZWVaf78+d52V/pxulxpNx03bNgwZWVl9Rj1NDY29vgfx5Xmscce07Zt2/TOO+/o2muv9daXlpZK0hV/zHbv3q3GxkZNnDhR2dnZys7OVl1dnV588UVlZ2d7x+JKP04jRozQTTfdFFh344036ujRo5I4ny548skntWLFCt1///0aP368HnroIT3xxBOqqamRxHGKStoFoZycHE2cOFG1tbWB9bW1tZoyZYqjXrlljNGSJUu0detWvf322yovLw/Ul5eXq7S0NHDMOjo6VFdXd0UdszvvvFN79+7Vnj17vNekSZM0b9487dmzR2PGjOE4SZo6dWqPS/wPHjyo0aNHS+J8uuDMmTM9vhU0KyvLu0Sb4xQRhxdFXNSFS7R/9rOfmf3795vKykozaNAg89///td115z43ve+Z4qKisy7775rPvvsM+915swZb5vnn3/eFBUVma1bt5q9e/eaBx54gEtFjQlcHWcMx8mY85evZ2dnm+eee878+9//Nr/4xS9Mfn6+ef31171tOE7GzJ8/31xzzTXeJdpbt241w4YNM8uXL/e24ThdvrQMQsYY85Of/MSMHj3a5OTkmK985Sve5chXIkkJX5s2bfK26e7uNs8884wpLS018Xjc3HHHHWbv3r3uOp0m7CDEcTrv97//vRk3bpyJx+Nm7NixZsOGDYF6jpMxLS0tZunSpWbUqFEmNzfXjBkzxqxatcq0t7d723CcLh/fJwQAcCbtckIAgCsHQQgA4AxBCADgDEEIAOAMQQgA4AxBCADgDEEIAOAMQQgA4AxBCADgDEEIAOAMQQgA4Mz/ATA4lyvTFmDgAAAAAElFTkSuQmCC" + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plt.imshow(np.log10(psf_185_20_lpi), cmap='Greys', vmin=-12, vmax=0)" + ], + "metadata": { + "collapsed": false, + "ExecuteTime": { + "end_time": "2024-05-09T21:05:03.286272Z", + "start_time": "2024-05-09T21:05:03.243932Z" + } + }, + "id": "1850b49a64d5abb6" + }, + { + "cell_type": "code", + "execution_count": 33, + "outputs": [ + { + "data": { + "text/plain": "Text(0, 0.5, 'Intensity')" + }, + "execution_count": 33, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkIAAAHFCAYAAAAe+pb9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAADn70lEQVR4nOydeXwU9f3/X7N3NsnmPgkJpyJyGtCCcnlgg7ettVWpF7UqrQdWW39+1WoPWlupVkFrraK2KvWi2lIRT6igQgQFuSEhIeS+s/cxvz8+85mZ3Z3Z3YQkG+D9fDzy2GR2js9ujnnl9b4EURRFEARBEARBnIAYkr0AgiAIgiCIZEFCiCAIgiCIExYSQgRBEARBnLCQECIIgiAI4oSFhBBBEARBECcsJIQIgiAIgjhhISFEEARBEMQJCwkhgiAIgiBOWEgIEQRBEARxwkJC6Djm888/x2WXXYbS0lJYrVYUFBRgxowZuOuuuwbl+qFQCC+99BLOPfdc5Obmwmw2Iz8/HxdeeCHeeecdhEIhAEB1dTUEQdD9+OUvf4mPP/445j7qDz0ir2M2m5GTk4Pp06fjzjvvxDfffBN1DL/uxx9/HLb9iSeewJgxY2CxWCAIAjo6OgAA//d//4fS0lKYTCZkZmb211uZEGvWrMEvf/lLzedGjBiB6667blDXAyDq+2Y0GlFQUIArrrgCu3bt6vX5Xn75ZTz22GNR210ul/xzMhTx+Xy4+eabUVRUBKPRiClTpgz4NUVRxPPPP4/TTz8dqampcDgcOO200/Cvf/0rbL8RI0Zo/h7dfPPNA77GeFRWVmLx4sWYOHEi0tPTUVBQgHPPPRcffvih5v4HDx7E5ZdfjszMTKSlpeG8887Dl19+mdC15s6diwkTJsTc55e//GXYe2SxWDBy5Ejcfvvt8t8AAFi5ciUEQcCWLVsSfq3Jhv99/OMf/5jspQw6pmQvgBgY/vOf/+Diiy/G3Llz8cgjj6CoqAj19fXYsmULXn31VTz66KMDen2Px4NLL70U7733Hr7//e/jqaeeQmFhIZqbm/Huu+/iiiuuwKpVq3DJJZfIx/z0pz/FVVddFXWukpISOBwObNq0KWz7ZZddhtGjR/f6F5dfJxQKoaOjA1u3bsVzzz2HJ554AkuXLsXdd98t73vaaadh06ZNGD9+vLxt27ZtuO2227Bo0SJce+21MJlMSE9Px7/+9S/85je/wX333YeKigpYrdZeretoWbNmDZYvX64pht566y04HI5BXY+a3/72t5g3bx58Ph+2bNmChx9+GB988AG2b9+OYcOGJXyel19+GTt27MAdd9wRtt3lcuGhhx4CwG5oQ42nnnoKf/nLX/DEE0+gvLwcaWlpA37NW265BStXrsSdd96JpUuXIhAIYPv27XC5XFH7nnnmmVG/RwUFBQO+xni88sor+OKLL3DDDTdg8uTJcDqdePrpp3HOOefghRdewA9/+EN53+bmZsyaNQtZWVl47rnnYLPZsHTpUsydOxebN2/GySef3G/revfdd5GRkYHu7m6sWbMGjz/+OL744gts3Lgx5j9jxBBFJI5LZs+eLY4ePVr0+/1RzwWDwQG//i233CICEF944QXN5/fu3St+9dVXoiiKYlVVlQhA/MMf/tCra5SVlYkXXHBBwvvHuo7L5RK//e1viwDENWvWxDzP3//+dxGA+Pnnn4dt//Wvfy0CEBsbG+OuxeVyJbzuRFm8eLE41H6lP/roIxGA+Nprr4Vt/9vf/iYCEH/961/36nwXXHCBWFZWFrW9ublZBCA++OCDR7HaaHw+n+bvUG9ZtGiRmJKS0g8rUoj1M/TWW2+JAMRVq1bFPU9vf48GE63fpUAgIE6aNEkcPXp02Pa7775bNJvNYnV1tbyts7NTzM3NFb/3ve/FvdacOXPEU089NeY+Dz74oAhAbG5uDtu+cOFCEYD4v//9TxRFUXz++edFAOLmzZvjXneo0Ne/w8cDFBo7TmltbUVubi5MpmjTz2AI/7bz8FMkkeEUbvd+9NFHuOWWW5Cbm4ucnBxcfvnlOHLkiLxfQ0MDnn32WZx//vlh/7GpGTt2LCZNmtS3FzcApKSk4G9/+xvMZjP+8Ic/yNsjQ2Nz587FNddcAwA444wzIAgCrrvuOowYMQL/93//B4D9J61+T0eMGIELL7wQb775JqZOnQqbzSa7F8uXL8fs2bORn5+P1NRUTJw4EY888gj8fn/UGt99912cc845yMjIgN1uxymnnIKlS5cCAK677josX74cAMKs++rqankNkaGxmpoaXHPNNcjPz4fVasUpp5yCRx99VA5ZAuF2+bJlyzBy5EikpaVhxowZ+Oyzz/r8fn/rW98CABw6dCjh92Hu3Ln4z3/+g0OHDkW9xry8PADAQw89JG9Xv959+/bhqquuCnut/P3i8O/1Sy+9hLvuugvDhg2D1WrF/v37cd111yEtLQ379+/HggULkJaWhuHDh+Ouu+6C1+uN+VoFQcCzzz4Lt9str23lypUAmHN67733YuTIkbBYLBg2bBgWL14cFmYBYv8MafH4449jxIgR+N73vhdzbUdDc3MzLBYL7r///qjndu/eDUEQ8Oc//xkAc+x+9rOfYeTIkbDZbMjOzsa0adPwyiuvxLxGfn5+1Daj0Yjy8nLU1taGbX/rrbdw9tlno6ysTN7mcDhw+eWX45133kEgEOjLy0yIyJ/nROF/Uz/88EP86Ec/Qk5ODhwOB374wx/C6XSioaEB3/ve95CZmYmioiL87Gc/i/rb4PP58Otf/xrjxo2D1WpFXl4err/+ejQ3N4ft9+GHH2Lu3LnIyclBSkoKSktL8Z3vfEfTIezP3/VjAQqNHafMmDEDzz77LG677TZcffXVOO2002A2m/vl3IsWLcIFF1yAl19+GbW1tbj77rtxzTXXyHH7jz76CH6/H5deemmvzhsKhTT/WGmJuYGguLgY5eXl2LhxIwKBgOZ1V6xYgVdeeQW//vWv8fzzz2PcuHHIy8vD7bffjuXLl+Nvf/ubbJuXlJTIx3355ZfYtWsX/u///g8jR45EamoqAODAgQO46qqr5BvhV199hd/85jfYvXs3nnvuOfn4v/3tb/jRj36EOXPm4Omnn0Z+fj727t2LHTt2AADuv/9+OJ1OvP7662EhxKKiIs3X2tzcjJkzZ8Ln8+FXv/oVRowYgX//+9/42c9+hgMHDmDFihVh+y9fvhzjxo2T83Puv/9+LFiwAFVVVcjIyOj1e71//34AkAVMIu/DihUrcNNNN+HAgQN46623wl7ju+++i29/+9u48cYbsWjRorBz79y5EzNnzkRpaSkeffRRFBYWYu3atbjtttvQ0tKCBx98MGxt9957L2bMmIGnn34aBoNBvhn7/X5cfPHFuPHGG3HXXXdh/fr1+NWvfoWMjAw88MADuq9106ZN+NWvfoWPPvpI/h0ZPXo0RFHEpZdeig8++AD33nsvZs2aha+//hoPPvggNm3ahE2bNoWFV/V+hiIJBALYtGkTFixYgGXLluHxxx/H4cOHUVZWhltvvRV33XVXVPhm/fr1SE9Ph8fjwdixY3HjjTfijjvugNFo1H1deXl5uPDCC/HCCy/goYceCvsH6/nnn4fFYsHVV18NAFiyZAleeukl/PrXv8bUqVPhdDqxY8cOtLa26p5fj0AggA0bNuDUU0+Vt7ndbhw4cACXXXZZ1P6TJk2C2+3GwYMHcdJJJ/X6eokQ+fPcWxYtWoTLL78cr776KrZu3Yr/9//+HwKBAPbs2YPLL78cN910E95//338/ve/R3FxMZYsWQKA/c285JJLsGHDBtxzzz2YOXMmDh06hAcffBBz587Fli1bkJKSgurqalxwwQWYNWsWnnvuOWRmZqKurg7vvvsufD4f7Ha7vJb+/l0/Jki2JUUMDC0tLeJZZ50lAhABiGazWZw5c6a4dOlSsbu7O2xf6IQUysrKxGuvvVb+mtu9t956a9h+jzzyiAhArK+vF0VRFH/3u9+JAMR33303obVyS1bvY8OGDZrH9WdojHPllVeGhbd4aOejjz6S99GzvfVs87KyMtFoNIp79uyJub5gMCj6/X7xxRdfFI1Go9jW1iaKoih2d3eLDodDPOuss8RQKKR7fKzQWOT38he/+IVmeO+WW24RBUGQ18rfs4kTJ4qBQEDe74svvhABiK+88krM18Tfv1WrVol+v190uVzi+vXrxTFjxohGo1EOjybyPohi30Jj559/vlhSUiJ2dnaGbf/JT34i2mw2+fx8rbNnz446x7XXXisCEP/5z3+GbV+wYIF48sknx3wP+PGpqalh2959910RgPjII4+EbV+1apUIQHzmmWfkbYn+DImiKNbX14sARIfDIZaUlIgvvPCC+MEHH4g333yzCED8f//v/4Xtf+utt4rPPfec+Mknn4irV68Wr776ahGAeM0118S91ttvvy0CEN977z15WyAQEIuLi8XvfOc78rYJEyaIl156adzzJcJ9990nAhBXr14tb6urqxMBiEuXLo3a/+WXXxYBiBs3box53t6ExhoaGkS/3y+2t7eLf//738WUlBRx+PDhotvtFkUx8dAY3++nP/1p2PZLL71UBCAuW7YsbPuUKVPE0047Tf76lVdeEQGIb7zxRth+mzdvFgGIK1asEEVRFF9//XURgLht2zbdtRzt7/qxDIXGjlNycnKwYcMGbN68Gb/73e9wySWXYO/evbj33nsxceJEtLS09PncF198cdjXPMTVW1s4kttvvx2bN2+O+hiMChuOKIoDct5JkyZp/je6detWXHzxxcjJyYHRaITZbMYPf/hDBINB7N27FwCwceNGdHV14dZbb+23RMwPP/wQ48ePx+mnnx62/brrroMoilFVORdccEGYO9Db7/mVV14Js9kMu92O2bNnIxgM4vXXX5fPk8j70Bc8Hg8++OADXHbZZbDb7QgEAvLHggUL4PF4omz/73znO5rnEgQBF110Udi2SZMm9fnnnr/HkSHLK664Aqmpqfjggw+irpWIo8FDm11dXXjttdfwwx/+EGeffTaeeuopXHrppVi2bBl6enrk/ZcvX47rr78es2fPxiWXXIK///3v+MlPfoK///3v2Lp1a8xrVVRUoLCwEM8//7y8be3atThy5AhuuOEGedvpp5+O//73v/jFL36Bjz/+GG63O+7r0OLZZ5/Fb37zG9x1111hhRacWL8f/ZnEXFhYCLPZjKysLFxzzTU47bTT8O6778Jms/XpfBdeeGHY16eccgoA9nsXuV398/bvf/8bmZmZuOiii8J+tqdMmYLCwkI5pD9lyhRYLBbcdNNNeOGFF3Dw4EHdtRzt7/qxCAmh45xp06bh5z//OV577TUcOXIEd955J6qrq/HII4/0+Zw5OTlhX3P7nv9xKy0tBQBUVVX16rwlJSWYNm1a1MdgVNhwDh06BKvViuzs7H49r1aIqqamBrNmzUJdXR0ef/xxWbjy3BX+fvJYvzrUdrS0trZqrqm4uFh+Xk2873k8fv/732Pz5s348ssvUVNTg4MHD8qh00Tfh77Q2tqKQCCAJ554AmazOexjwYIFABD1T4FeONFut0fd6KxWKzweT5/XZjKZosIpgiCgsLAw6nugt65IsrKyIAgCHA6HnLvCqaiogMfjwc6dO2Oeg+fBxcsNMZlMWLhwId566y05r2nlypUoKirC+eefL+/35z//GT//+c+xevVqzJs3D9nZ2bj00kuxb9++hF4TwMJtP/7xj3HTTTeF5fEBymvWCrW1tbUBQL/+Tr///vvYvHkztm3bhpaWFvzvf/8LqyztLZFrs1gsutvVP2+NjY3o6OiAxWKJ+vluaGiQf7ZHjx6N999/H/n5+Vi8eDFGjx6N0aNH4/HHH49ay9H+rh+LUI7QCYTZbMaDDz6IP/3pT3JuCcB+0LUSPvsSvweAefPmwWw2Y/Xq1UOiF0mi1NXVobKyEnPmzOn3vCSt/0ZXr14Np9OJN998MyzBc9u2bWH78Rvl4cOH+209OTk5qK+vj9rOk95zc3P77VoAMGrUKEybNk3zuUTfh76QlZUFo9GIhQsXYvHixZr7jBw5MuzrwSp/zsnJQSAQQHNzc5gYEkURDQ0NmD59ep/WlZKSgrFjx6KhoSHqOe54RhZM9HU/ALj++uvxhz/8Aa+++iquvPJKvP3221H5RampqXjooYfw0EMPobGxUXaHLrroIuzevTvuNZ5//nm5XcXTTz8d9V6kpKRgzJgx2L59e9Sx27dvR0pKCkaNGhX3OokyefLkfv8d6Qu8YOXdd9/VfD49PV3+fNasWZg1axaCwSC2bNmCJ554AnfccQcKCgrw/e9/f7CWPCQhR+g4ResmB0BuYsf/8wdYRcrXX38dtt+HH34YZp/3hsLCQixatAhr167Fiy++qLnPgQMHoq6ZTNxuNxYtWoRAIIB77rlnUK7J/5irE2JFUcRf//rXsP1mzpyJjIwMPP300zFDd735z+2cc87Bzp07o5rNvfjiixAEAfPmzUv4dRwtib4PfB+t16f32u12O+bNm4etW7di0qRJmo5j5H/Ag8U555wDAPj73/8etv2NN96A0+mUn+8L3/nOd9DV1YWNGzeGbV+zZg3S0tLCEo214L+3kY6SFqeccgrOOOMMPP/883j55Zfh9Xpx/fXX6+5fUFCA6667Dj/4wQ+wZ88ezaolNStXrsSiRYtwzTXX4Nlnn9UVhJdddhk+/PDDsGqy7u5uvPnmm7j44osHrehiMLnwwgvR2tqKYDCo+bOt1TvJaDTijDPOkB3XRBtOHs8cfz8ZBADg/PPPR0lJCS666CKMGzcOoVAI27Ztw6OPPoq0tDTcfvvt8r4LFy7E/fffjwceeABz5szBzp078eSTTx5VhcCyZctw8OBBXHfddVi7di0uu+wyFBQUoKWlBevWrcPzzz+PV199NayEvqamRtOKz8vLw+jRo/u8lkj4dUKhEDo7O+WGiocOHcKjjz6K+fPn99u1YnHeeefBYrHgBz/4Ae655x54PB489dRTaG9vD9svLS0Njz76KBYtWoRzzz0XP/rRj1BQUID9+/fjq6++wpNPPgkAmDhxIgAWhqqoqIDRaMSkSZNkm13NnXfeiRdffBEXXHABHn74YZSVleE///kPVqxYgVtuuWXAqmu0SPR9ANhrfPPNN/HUU0+hvLwcBoMB06ZNQ3p6OsrKyvCvf/0L55xzDrKzs5Gbm4sRI0bg8ccfx1lnnYVZs2bhlltuwYgRI9Dd3Y39+/fjnXfe0e1SPNCcd955OP/88/Hzn/8cXV1dOPPMM+WqsalTp2LhwoV9PvfPfvYz/OMf/8AVV1yBX/3qVygpKcHrr7+Ot99+G3/84x+RkpICgDWofPPNN3HBBRegrKwMHR0deO211/Dqq6/iuuuuw+TJkxO63g033IAf//jHOHLkCGbOnBl1Az7jjDNw4YUXYtKkScjKysKuXbvw0ksvYcaMGWEVS5G89tpruPHGGzFlyhT8+Mc/xhdffBH2/NSpU2UR/LOf/QwvvfSS/DNttVrxu9/9Dh6PR7fjeiRdXV14/fXXo7bn5eVhzpw5CZ1jMPn+97+Pf/zjH1iwYAFuv/12nH766TCbzTh8+DA++ugjXHLJJbjsssvw9NNP48MPP8QFF1yA0tJSeDweuRrz3HPPTfKrGAIkL0+bGEhWrVolXnXVVeLYsWPFtLQ00Ww2i6WlpeLChQvFnTt3hu3r9XrFe+65Rxw+fLiYkpIizpkzR9y2bZtu1VhkJYRWZZUosuqRF154QTz77LPF7Oxs0WQyiXl5eWJFRYX48ssvy40d41WNXX311Zqvsa9VY/zDaDSKWVlZYnl5uXjHHXeI33zzTdQx/VU1prfOd955R5w8ebJos9nEYcOGiXfffbf43//+V/P9XLNmjThnzhwxNTVVtNvt4vjx48Xf//738vNer1dctGiRmJeXJwqCIAIQq6qq5DWov5eiKIqHDh0Sr7rqKjEnJ0c0m83iySefLP7hD38Ia7gZq9IOCTQw1Guo2Nf3oa2tTfzud78rZmZmyq+R8/7774tTp04VrVarCCDs9VZVVYk33HCDOGzYMNFsNot5eXnizJkzwxo6xlqrVtWXKCrf83joHe92u8Wf//znYllZmWg2m8WioiLxlltuEdvb28P260vTw5qaGvH73/++mJWVJVosFnHSpEnic889F7bPpk2bxHPOOUcsLCwUzWazaLfbxenTp4srVqzoVePVzs5OMSUlRQQg/vWvf416/he/+IU4bdo0MSsrS7RareKoUaPEO++8U2xpaYl5Xl6tp/fBf745+/fvFy+99FLR4XCIdrtdPOecc8TKysqEXsOcOXN0rzNnzhxRFPV/xyPpbdVYon9LtH6O/H6/+Mc//lH+/UlLSxPHjRsn/vjHPxb37dsniiL7Pl922WViWVmZaLVaxZycHHHOnDni22+/LZ/naH/Xj2UEURygMhmCIAiCIIghDuUIEQRBEARxwnJCCKF///vfOPnkkzF27Fg8++yzyV4OQRAEQRBDhOM+NBYIBDB+/Hh89NFHcDgcOO200/D555/3e58YgiAIgiCOPY57R+iLL77AqaeeimHDhiE9PR0LFizA2rVrk70sgiAIgiCGAENeCK1fvx4XXXQRiouLIQgCVq9eHbXPihUr5KnG5eXl2LBhg/zckSNHMGzYMPnrkpIS1NXVDcbSCYIgCIIY4gx5IeR0OjF58mS5V0okq1atwh133IH77rsPW7duxaxZs1BRUYGamhoA2rOjBqtzLEEQBEEQQ5sh31CxoqICFRUVus8vW7YMN954IxYtWgQAeOyxx7B27Vo89dRTWLp0KYYNGxbmAB0+fBhnnHGG7vm8Xm/YuIlQKIS2tjbk5OSQgCIIgiCIYwRRFNHd3Y3i4uLY42KS2cSotwAQ33rrLflrr9crGo1G8c033wzb77bbbhNnz54tiiJrNjVmzBjx8OHDYldXlzhmzJiYTbx4Iyv6oA/6oA/6oA/6OPY/amtrY2qLIe8IxaKlpQXBYBAFBQVh2wsKCuSBgyaTCY8++ijmzZuHUCiEe+65J+ZsoXvvvRdLliyRv+7s7ERpaSlqa2vhcDgG5oUQBEEQBNGvdHV1Yfjw4WHDZ7U4poUQJzJkJYpi2LaLL74YF198cULnslqtYcMfOQ6Hg4QQQRAEQRxjxEtrGfLJ0rHIzc2F0WiU3R9OU1NTlEvUW5YvX47x48dj+vTpR3UegiAIgiCGLse0ELJYLCgvL8e6devCtq9btw4zZ848qnMvXrwYO3fuxObNm4/qPARBEARBDF2GfGisp6cH+/fvl7+uqqrCtm3bkJ2djdLSUixZsgQLFy7EtGnTMGPGDDzzzDOoqanBzTfffFTXXb58OZYvX45gMHi0L4EgCIIgiCHKkB+x8fHHH2PevHlR26+99lqsXLkSAGuo+Mgjj6C+vh4TJkzAn/70J8yePbtfrt/V1YWMjAx0dnZSjhBBEARBHCMkev8e8kIo2ZAQIgiCIIhjj0Tv38d0jtBAQsnSBEEQBHH8Q45QHMgRIgiCIIhjD3KECIIgCIIg4kBCiCAIgiCIExYSQjpQjhBBEARBHP9QjlAcKEeIIAiCII49KEeIIAiCIAgiDiSECIIgCII4YSEhpAPlCBEEQRDE8Q/lCMWBcoQIgiAI4tiDcoQIgiAIgiDiQEKIIAiCIIgTFhJCBEEQBEGcsJAQ0oGSpQmCIAji+IeSpeNAydIEQRAEcexBydIEQRAEQRBxICFEEARBEMQJCwkhgiAIgiBOWEgIEQRBEARxwkJCiCAIgiCIExYSQjoMePn8V68C790PHNk2MOcnCIIgCCIuVD4fhwErn//HFcC+94CL/gyUX9t/5yUIgiAIgsrnhzy2TPbo6UjmKgiCIAjihIaEULJIyWSP7o5kroIgCIIgTmhICCULcoQIgiAIIumQEEoW5AgRBEEQRNIhIZQk9veYAQCurtYkr4QgCIIgTlxICCWJf+9xAQA8JIQIgiAIImmQEEoSgj0LAGD0dSZ5JQRBEARx4kJCSIeBbqholISQ2d81IOcnCIIgCCI+JIR0WLx4MXbu3InNmzcPyPktqUwIWQM9APW0JAiCIIikQEIoSVgcOQAAI4KAtzvJqyEIgiCIExMSQkkiLTUdXtHEvqBeQgRBEASRFEgIJYlMuwVdSGVfUC8hgiAIgkgKJISSRIbdjE5REkLkCBEEQRBEUiAhlCQyU8zoJEeIIAiCIJIKCaEkoXaEQu72JK+GIAiCIE5MSAgliQyVI+TtbkvyagiCIAjixISEUJKwmoxwCmkAAF8PCSGCIAiCSAYnhBC67LLLkJWVhe9+97vJXkoYXpMDAOB3khAiCIIgiGRwQgih2267DS+++GKylxFFwMKEUMjVkdyFEARBEMQJygkhhObNm4f09PRkLyOKoDWDfUJVYwRBEASRFJIuhNavX4+LLroIxcXFEAQBq1evjtpnxYoVGDlyJGw2G8rLy7Fhw4bBX+gAINoyAQAGb0dS10EQBEEQJypJF0JOpxOTJ0/Gk08+qfn8qlWrcMcdd+C+++7D1q1bMWvWLFRUVKCmpkbep7y8HBMmTIj6OHLkyGC9jD5hsGcCAEw+mkBPEARBEMnAlOwFVFRUoKKiQvf5ZcuW4cYbb8SiRYsAAI899hjWrl2Lp556CkuXLgUAVFZW9tt6vF4vvF6v/HVX18CJFBMXQv6eAbsGQRAEQRD6JN0RioXP50NlZSXmz58ftn3+/PnYuHHjgFxz6dKlyMjIkD+GDx8+INcBgBQ7K583hnwDdg2CIAiCIPQZ0kKopaUFwWAQBQUFYdsLCgrQ0NCQ8HnOP/98XHHFFVizZg1KSkqwefNm3X3vvfdedHZ2yh+1tbV9Xn887HbWUNEkeuPsSRAEQRDEQJD00FgiCIIQ9rUoilHbYrF27dqE97VarbBarVi+fDmWL1+OYDCY8LG9JTWVCSGz6AdEEejFayIIgiAI4ugZ0o5Qbm4ujEZjlPvT1NQU5RL1N4sXL8bOnTtjukdHS3paqvJFgFwhgiAIghhshrQQslgsKC8vx7p168K2r1u3DjNnzkzSqvqP9NQ05YuAJ3kLIQiCIIgTlKSHxnp6erB//37566qqKmzbtg3Z2dkoLS3FkiVLsHDhQkybNg0zZszAM888g5qaGtx8880Duq7BCI2l2FIQEgUYBJEcIYIgCIJIAoIoimIyF/Dxxx9j3rx5UduvvfZarFy5EgBrqPjII4+gvr4eEyZMwJ/+9CfMnj17UNbX1dWFjIwMdHZ2wuFw9Ou5j3S4kfWnUqQIPuD2r4Gssn49P0EQBEGcqCR6/066EBrqDKQQ6nD5gN+PQKbghP+Wz2EuGNev5ycIgiCIE5VE799DOkcomSxfvhzjx4/H9OnTB+waNrMRHlgAAF63c8CuQxAEQRCENiSEdBiMqjGryQAvzAAAn8c1YNchCIIgCEIbEkJJRBAE+Lkj5HEneTUEQRAEceJBQijJ+AUmhMgRIgiCIIjBh4SQDoORIwQAAQMTQn4vOUIEQRAEMdiQENJhMHKEAEUIBXwkhAiCIAhisCEhlGSCBisAcoQIgiAIIhmQEEoyISNzhII+GrFBEARBEIMNCSEdBitHKGSwAQCCFBojCIIgiEGHhJAOg5UjJJpYaCzkJyFEEARBEIMNCaEkIxqZIyT6KTRGEARBEIMNCaFkIzlCJIQIgiAIYvAhIZRszMwRQoCEEEEQBEEMNiSEdBisZGmD5Agh4B3Q6xAEQRAEEQ0JIR0GK1laMKewxyAJIYIgCIIYbEgIJRmjhYXGSAgRBEEQxOBDQijJGC12AICBhBBBEARBDDokhJIMd4SMQUqWJgiCIIjBhoRQkjFZWY6QKeRL8koIgiAI4sSDhFCSMVskISSSECIIgiCIwYaEkA6DVT5vtpEjRBAEQRDJgoSQDoNVPm+xsWRpMzlCBEEQBDHokBBKMlYuhOCDKIpJXg1BEARBnFiQEEoyXAhZ4YcvGEryagiCIAjixIKEUJKxpqSyR/jh9gWTvBqCIAiCOLEgIZRkeNWYFX64SAgRBEEQxKBCQijZSENXzUIQLg91lyYIgiCIwYSEULIx2eRPvW5XEhdCEARBECceJISSjeQIAYDXQ0KIIAiCIAYTEkI6DFZDRRiM8MMEAPB6nAN7LYIgCIIgwiAhpMNgNVQEAL9gBgD4PO4BvxZBEARBEAokhIYAAcECgIQQQRAEQQw2JISGAH4DyxPyeSlHiCAIgiAGExJCQ4CgJIQoWZogCIIgBhcSQkMA0SgJITclSxMEQRDEYEJCaAggmnhojHKECIIgCGIwISE0FJCaKvooNEYQBEEQgwoJoSGAIDlCAa8nySshCIIgiBMLEkJDAIOZDV4N+sgRIgiCIIjB5LgXQrW1tZg7dy7Gjx+PSZMm4bXXXkv2kqIwShPoQ35yhAiCIAhiMDElewEDjclkwmOPPYYpU6agqakJp512GhYsWIDU1NRkL03GaJPW4idHiCAIgiAGk+NeCBUVFaGoqAgAkJ+fj+zsbLS1tQ0pIWSy2tknAQ9EUYQgCMldEEEQBEGcICQ9NLZ+/XpcdNFFKC4uhiAIWL16ddQ+K1aswMiRI2Gz2VBeXo4NGzb06VpbtmxBKBTC8OHDj3LV/YvFxoSQVfTC6QsmeTUEQRAEceKQdCHkdDoxefJkPPnkk5rPr1q1CnfccQfuu+8+bN26FbNmzUJFRQVqamrkfcrLyzFhwoSojyNHjsj7tLa24oc//CGeeeaZAX9NvYU7Qjb40OX2J3k1BEEQBHHikPTQWEVFBSoqKnSfX7ZsGW688UYsWrQIAPDYY49h7dq1eOqpp7B06VIAQGVlZcxreL1eXHbZZbj33nsxc+bMuPt6vV75666urkRfSp8RpKoxG/zodPtRnJky4NckCIIgCGIIOEKx8Pl8qKysxPz588O2z58/Hxs3bkzoHKIo4rrrrsPZZ5+NhQsXxt1/6dKlyMjIkD8GJYxmZo5QiuAlR4ggCIIgBpEhLYRaWloQDAZRUFAQtr2goAANDQ0JnePTTz/FqlWrsHr1akyZMgVTpkzB9u3bdfe/99570dnZKX/U1tYe1WtICKmztA0+dJIQIgiCIIhBI+mhsUSIrKLqTWXVWWedhVAolPC1rFYrrFYrli9fjuXLlyMYHITkZTk05kOTJzDw1yMIgiAIAsAQd4Ryc3NhNBqj3J+mpqYol6i/Wbx4MXbu3InNmzcP6HUAKEJIoGRpgiAIghhMhrQQslgsKC8vx7p168K2r1u3Lm7S8zGFSXGEKDRGEARBEINH0kNjPT092L9/v/x1VVUVtm3bhuzsbJSWlmLJkiVYuHAhpk2bhhkzZuCZZ55BTU0Nbr755gFdVzJCYynwoctDQoggCIIgBoukC6EtW7Zg3rx58tdLliwBAFx77bVYuXIlrrzySrS2tuLhhx9GfX09JkyYgDVr1qCsrGxA17V48WIsXrwYXV1dyMjIGNBrqXOEyBEiCIIgiMEj6UJo7ty5EEUx5j633norbr311kFaURLgVWOCF11uSpYmCIIgiMFiSOcIJZPly5dj/PjxmD59+sBfTNVQsabNiXU7GxEIJl7pRhAEQRBE3xDEeHbMCQ4PjXV2dsLhcAzQRY4Ay06BXzRirPclAMCfrpyMy6aWDMz1CIIgCOI4J9H7NzlCQwHJETILQZjAQmM7jwz8aA+CIAiCONEhIaTDoIbGTMpsMRt8AIDaNvfAX5cgCIIgTnBICOkwqA0VTVYArFP201eeAgA41OYa+OsSBEEQxAkOCaGhgCDI4bGSdCaIattccavpCIIgCII4OkgIDRWkEvoiNogePd4A2l3UU4ggCIIgBhISQkMFM1NAVtGLQgcTRYdanclcEUEQBEEc95AQ0mFQk6UBwMzEDwIelGYzUVRDeUIEQRAEMaCQENJhUJOlAaVyzO9CaY4khFpJCBEEQRDEQEJCaKhg5kKIHCGCIAiCGCxICA0VeGjM70ZZDgkhgiAIghgMSAgNFaRkaQTcGC45Qgeae2jmGEEQBEEMICSEdBj0ZGkTd4Q8GF/kQJbdjJYeH975+sjgXJ8gCIIgTkBICOkw6MnSZiVZ2mY2YtGsUQCAJz7cj2CIGisSBEEQxEBAQmiowIVQwAMA+OGMMmSkmHGw2Yl1OxuTuDCCIAiCOH4hITRUUJXPA0C6zYwFEwsBADvraRI9QRAEQQwEJISGCqryec6wTLatoZMm0RMEQRDEQEBCaKggd5ZWRE9hBhNC9Z0erSMIgiAIgjhKSAjpMPhVY9wRUoRQUQYTRySECIIgCGJgICGkQ/KqxtSOEBNCDSSECIIgCGJAICE0VNAQQtwR6vEG0OXxJ2NVBEEQBHFcQ0JoqBBRPg8AdosJGSlmAOQKEQRBEMRAQEJoqKCRIwRQnhBBEARBDCQkhIYKqqGrapQ8ISqhJwiCIIj+hoTQUEE1dFVNkVRCf6SDHCGCIAiC6G9ICA0VeI6QtztscxFVjhEEQRDEgEFCaKiQNRIwmABXK9B+SN7MQ2P1XSSECIIgCKK/ISE0VLCmAcWnsc+rN8ib5WTpDsoRIgiCIIj+hoSQDoPeWRoARpzFHqv/J29ScoTcEEVx8NZCEARBECcAJIR0GPTO0gAwchZ7rNoASKJneHYKjAYBTl8QTd3ewVsLQRAEQZwAkBAaSgw/AzCYga7DQHsVAMBqMqI0m1WU7W/qSebqCIIgCOK4g4TQUMKSCgwrZ5+rwmOj89IAAAeaSQgRBEEQRH9CQmiooQ6PSYzOTwUAHCBHiCAIgiD6FRJCQ40RkhCqVvKEuCO0nxwhgiAIguhXSAgNNYafDhgtQHc90HYQADAmXwqNNTmTuTKCIAiCOO4gITTUMKcAJVLJftV6AIoj1NDlQbfHn6yVEQRBEMRxBwmhoYg6PAYgI8WMvHQrAOBgM7lCBEEQBNFfkBAaisiNFT+VN43OkxKmKU+IIAiCIPqN414IdXd3Y/r06ZgyZQomTpyIv/71r8leUnyGlQOCAehpALobAABl2UwI1bXTqA2CIAiC6C9MyV7AQGO32/HJJ5/AbrfD5XJhwoQJuPzyy5GTk5PspeljsQO5JwPNu4Aj24CTv43MVDMAoN1FOUIEQRAE0V8c946Q0WiE3c46M3s8HgSDwWNjZlfRZPZY/xUAINtuAQB0uHzJWhFBEARBHHckXQitX78eF110EYqLiyEIAlavXh21z4oVKzBy5EjYbDaUl5djw4YN0SeKQUdHByZPnoySkhLcc889yM3N7afVDyDFU9hj/TYAQJYkhNpJCBEEQRBEv5F0IeR0OjF58mQ8+eSTms+vWrUKd9xxB+677z5s3boVs2bNQkVFBWpqauR9ysvLMWHChKiPI0eOAAAyMzPx1VdfoaqqCi+//DIaGxsH5bUdFRGOUKadhcbaKDRGEARBEP1Gn3KErrvuOtxwww2YPXv2US+goqICFRUVus8vW7YMN954IxYtWgQAeOyxx7B27Vo89dRTWLp0KQCgsrIyoWsVFBRg0qRJWL9+Pa644grNfbxeL7xeZcp7V1dXoi+lfymcCEAAuuqAnmZkp1JojCAIgiD6mz45Qt3d3Zg/fz7Gjh2L3/72t6irq+vvdQEAfD4fKisrMX/+/LDt8+fPx8aNGxM6R2Njoyxmurq6sH79epx88sm6+y9duhQZGRnyx/Dhw/v+Ao4GazqQM4Z9Xv8VMnlozElCiCAIgiD6iz4JoTfeeAN1dXX4yU9+gtdeew0jRoxARUUFXn/9dfj9/Re6aWlpQTAYREFBQdj2goICNDQ0JHSOw4cPY/bs2Zg8eTLOOuss/OQnP8GkSZN097/33nvR2dkpf9TW1h7VazgqiqR1Nm5HlhQa6/IEEAiGkrcmgiAIgjiO6HP5fE5ODm6//Xbcfvvt2Lp1K5577jksXLgQaWlpuOaaa3Drrbdi7Nix/bJIQRDCvhZFMWqbHuXl5di2bVvC17JarbBarVi+fDmWL1+OYDDYm6X2L7mSc9WyHxkpZggCm8Pa4fYjN82avHURBEEQxHHCUSdL19fX47333sN7770Ho9GIBQsW4JtvvsH48ePxpz/96ajOnZubC6PRGOX+NDU1RblE/c3ixYuxc+dObN68eUCvE5NcKTTWug8mowEOG3OFKE+IIAiCIPqHPgkhv9+PN954AxdeeCHKysrw2muv4c4770R9fT1eeOEFvPfee3jppZfw8MMPH9XiLBYLysvLsW7durDt69atw8yZM4/q3McEOZKj1rIPAOTwWJuTKscIgiAIoj/oU2isqKgIoVAIP/jBD/DFF19gypQpUfucf/75yMzMjHuunp4e7N+/X/66qqoK27ZtQ3Z2NkpLS7FkyRIsXLgQ06ZNw4wZM/DMM8+gpqYGN998c1+WnjBDIjSWM5o9utsAZyuyUi2obnVRLyGCIAiC6Cf6JIT+9Kc/4YorroDNZtPdJysrC1VVVXHPtWXLFsybN0/+esmSJQCAa6+9FitXrsSVV16J1tZWPPzww6ivr8eECROwZs0alJWV9WXpCbN48WIsXrwYXV1dyMjIGNBr6WJJBTKGA521QOs+uakihcYIgiAIon/oU2jso48+0qwOczqduOGGG3p1rrlz50IUxaiPlStXyvvceuutqK6uhtfrRWVlZb/0Lzpm4CX0LYoQotAYQRAEQfQPfRJCL7zwAtzu6CnobrcbL7744lEvaiiwfPlyjB8/HtOnT0/uQnKlPKHWfXKOEDlCBEEQBNE/9Co01tXVJTs23d3dYaGxYDCINWvWID8/v98XmQyGRGgMCEuYziqkeWMEQRAE0Z/0SghlZmZCEAQIgoCTTjop6nlBEPDQQw/12+IIKCX0LfuQNYpCYwRBEATRn/RKCH300UcQRRFnn3023njjDWRnZ8vPWSwWlJWVobi4uN8XmQyGRNUYAGRLlWPt1chKMQKg0BhBEARB9Be9EkJz5swBwErcS0tLE+7ufCwyZEJjjmJAMAAhP/IENjONQmMEQRAE0T8kLIS+/vprTJgwAQaDAZ2dndi+fbvuvrFmeRG9xGgG0ouBrsPICzUBANpdFBojCIIgiP4gYSE0ZcoUNDQ0ID8/H1OmTIEgCBBFMWo/QRCSH0463sgoAboOI9PfCCAVHS4fQiERBsPx68gRBEEQxGCQsBCqqqpCXl6e/PnxzpDJEQKYEKoFUt31AMYgJALdngAypHJ6giAIgiD6RsJCSN3JeaC7Og8FhkyOEMCEEABTdx3SrOPQ4w2gzeUjIUQQBEEQR0mfGyr+5z//kb++5557kJmZiZkzZ+LQoUP9tjhCQhJC6DyMTEn8xEqYFkURgWAIAPCXTw5g5afHv4NHEARBEH2hT0Lot7/9LVJSUgAAmzZtwpNPPolHHnkEubm5uPPOO/t1gQSAzFL22Fkbd96YLxDCOcs+wWUrNuJAcw+W/nc3Hvr3TnR5KMGaIAiCICLp09DV2tpajBnDGv2tXr0a3/3ud3HTTTfhzDPPxNy5c/tzfQQQ5ghl5cVuqljX4cbBZicA4G//Y06QKAK1bS6cWpzkEB9BEARBDDH65AilpaWhtbUVAPDee+/h3HPPBQDYbDbNGWTEUcKFkLsNBbYAAH1HqLHLI3/+z8218ue1bfR9IQiCIIhI+uQInXfeeVi0aBGmTp2KvXv34oILLgAAfPPNNxgxYkR/ri9pDKmqMVsGYHUA3i6UGdsB6OcIqYVQIKS0Nzjc7hrYNRIEQRDEMUifHKHly5djxowZaG5uxhtvvIGcnBwAQGVlJX7wgx/06wKTxeLFi7Fz505s3rw52UthSK7QMANz4vRCY2ohpKamjYQQQRAEQUTSJ0coMzMTTz75ZNR2Grg6gGSUAE07USg2AcjRDY01dHo1t9eSECIIgiCIKPokhACgo6MDX3zxBZqamhAKheTtgiBg4cKF/bI4QoWDDbPNCsUJjXUzR8hhM6HLE8DpI7LxRXUbatspR4ggCIIgIumTEHrnnXdw9dVXw+l0Ij09PWz4KgmhASKVdfXOCEpCSC801smE0AMXnQqzUcD4IgfO+9N6HG53QRTF43pQLkEQBEH0lj7lCN1111244YYb0N3djY6ODrS3t8sfbW1t/b1GAgBS8wEA9gB7f+M5QiNz7bhkyjCMyE2FQQA8/hCae7TDZgRBEARxotInIVRXV4fbbrsNdru9v9czZFi+fDnGjx+P6dOnJ3spjDTmCNm8TAh1uPxRQ29FUURjFxM7+ek2AIDZaEBRBmt+SXlCBEEQBBFOn4TQ+eefjy1btvT3WoYUQ65qTAqNmT0tAABfMASnL7y0v8Plhy/A8rXyHVZ5+/BsLoQoT4ggCIIg1PQpR+iCCy7A3XffjZ07d2LixIkwm8OHf1588cX9sjhChRQaMzhbYDUZ4A2E0O70Ic2qfAt5WCw71QKryShvH55lx2dow6f7W3DJlGLKEyIIgiAIiT4JoR/96EcAgIcffjjqOUEQhkYTwuMNKTQGbycKUoCabuYADc9WdmmQEqULHLawQ789oRCvVR7Ga5WHUZptx0/PGTtYqyYIgiCIIU2fQmOhUEj3g0TQAGHLBAzMeRthZ7k+bREJ001SflCBKiwGAOecUoAHLhwPAPjzh/vg8dP3iCAIgiCAPgohNR6Pdidjop8RBDlPaLiFDVWNbKrYIHWVLoxwhADg+jNHwGIywB8U0erUrjgjCIIgiBONPgmhYDCIX/3qVxg2bBjS0tJw8OBBAMD999+Pv/3tb/26QEJFai4AYJi5GwDQFiFoDjT3AAAKM6KFkCAIyLazyfXtJIQIgiAIAkAfhdBvfvMbrFy5Eo888ggsFou8feLEiXj22Wf7bXFEBGksYTrf0AUA6PEE5Kd8gRA+2t0EADhrTK7m4Zl2FlqLFFAEQRAEcaLSJyH04osv4plnnsHVV18No1GpTpo0aRJ2797db4sjIpAqxzLFDgCAW5Xrs+lgK7o8AeSlW3FaaZbm4dmpkiOk04yRIAiCIE40+txQccyYMVHbQ6EQ/H7t0Q/HGkOuoSIgh8Yygx0AwoXQuzsaAADzxxfAYNAuj8+ShBA5QgRBEATB6JMQOvXUU7Fhw4ao7a+99hqmTp161IsaCgy5hoqAHBpLl4QQr/4KhkSs28mE0LcnFOoeTjlCBEEQBBFOn/oIPfjgg1i4cCHq6uoQCoXw5ptvYs+ePXjxxRfx73//u7/XSHCk0FiaNG/MLXWWbur2oKXHB5NBwLdG5egeLjtCFBojCIIgCAB9dIQuuugirFq1CmvWrIEgCHjggQewa9cuvPPOOzjvvPP6e40ERwqN2f2SEJIcIZ407Ugxw2zU/5ZmS8nSepPr9QgEQ71eKkEQBEEcC/TJEQLYvLHzzz+/P9dCxEMKjaX42gEAbj8TKN1eJoTU4za06EuO0Pq9zfjRi1vwy4tPxQ9OL+31kgmCIAhiKNMnR2jUqFFobW2N2t7R0YFRo0Yd9aIIHezMEbL4OwCI8PjCHaF4QqgvVWPv7WyANxDCB7sae79egiAIghji9EkIVVdXa47S8Hq9qKurO+pFETrYMgAABjGIVHiU0FiijpC9947Q7nrWvPFQq6vXyyUIgiCIoU6vQmNvv/22/PnatWuRkZEhfx0MBvHBBx9gxIgR/bY4IgJzCps3FvLDAVdUjlCaLTFHqMPlhyiKcafQi6KIPQ2SEGpzIRQSdUvzCYIgCOJYpFdC6NJLLwXAxjVce+21Yc+ZzWaMGDECjz76aL8tjohAEJgr5GqBQ3CiRwqNJZwjJDlCvmAITl9Qd39RFNHc44UvEJLP7QuE0NDlQXFmSn+9GoIgCIJIOr0SQqEQS84dOXIkNm/ejNxc7VEOxADChRBcaOmlI5RiMcJmNsDjD6Hd6dMVQu/uaMAt//gS5WXhHaoPtbpICBEEQRDHFX3KEaqqqjrmRJDL5UJZWRl+9rOfJXspR0dKJgDAIThVOUKsHD49jiMEKE0VY+UJPf9pNQCg8lB72PZDrc7erpYgCIIghjR9Lp//4IMP8MEHH6CpqUl2ijjPPffcUS+sv/nNb36DM844I9nLOHqkhGmeIySKon6ydNcRIL2IhdQkslItONLpidlUMd9hDfvaZBAQCImopoRpgiAI4jijT47QQw89hPnz5+ODDz5AS0sL2tvbwz6GGvv27cPu3buxYMGCZC/l6OFCSHBBFFm+T7dWaGzdA8CyU4DXbwD8HnmzXEIfwxHqcIU3XJwpTbOvaSNHiCAIgji+6JMQevrpp7Fy5Up8/vnnWL16Nd56662wj96wfv16XHTRRSguLoYgCFi9enXUPitWrMDIkSNhs9lQXl6uOecsFj/72c+wdOnSXh0zZJEdISZKPL4QnJGOUNV64NPH2effvAm8ciUgigASK6Fv6fHKn5sMAhZI88v+u6MBt/6jMipkRhAEQRDHKn0Kjfl8PsycObNfFuB0OjF58mRcf/31+M53vhP1/KpVq3DHHXdgxYoVOPPMM/GXv/wFFRUV2LlzJ0pLWafj8vJyeL3eqGPfe+89bN68GSeddBJOOukkbNy4sV/WnFQkIZRpcAFBNmaDh8bSbSYg6Af+9RO275hzger/AQc/Buq/AoqnyI5QIkJo6eUTUZZtR77DBoBpqTXbGyAIQlQiNUEQBEEci/RJCC1atAgvv/wy7r///qNeQEVFBSoqKnSfX7ZsGW688UYsWrQIAPDYY49h7dq1eOqpp2SXp7KyUvf4zz77DK+++ipee+019PT0wO/3w+Fw4IEHHtDc3+v1homqrq6uvrysgUMSQtlGF+BnQkgOjVnNQP3XQMchtt8VK4HVtwK73gZ2/gsonoK8dJb/09Dl0Tx9MCTKIumccfnId9jgDYQ3zzzc7h6gF0cQBEEQg0ufhJDH48EzzzyD999/H5MmTYLZbA57ftmyZf2yOJ/Ph8rKSvziF78I2z5//vyE3Z2lS5fKgmnlypXYsWOHrgji+z/00EN9X/RAwx0hgYkRt09xhNJsJqD2M7Zf6QzAmg6Mv4QJoV1vA+c8gJIsVv5epyNm2l0+hESWX83dI6vJiNvOHoPnN1aj2xOAP0BDWAmCIIjjgz4Joa+//hpTpkwBAOzYsaM/1xNGS0sLgsEgCgoKwrYXFBSgoaFhQK557733YsmSJfLXXV1dGD58+IBcq0/YMgEAGQKr4FKHxtKsRqBGEkLDT2ePY+cDRivQuh9o2oXiTJbvc6RTWwjxsFiW3QKTapL9kvkn49zxBbj4yU97NaKDIAiCIIYyfRJCH330UX+vIyaRoyASGQ+hxXXXXRd3H6vVCqvVGne/pCEJIYfAkqXdvqDSUNFiAmo/Z/sN/5a0vwMYfTaw97/A7v9g2GSWP1Tf4UEwJMIYMTKjtYeJnBzJDVKjzi/q6/eAIAiCIIYSvRJCl19+edx9BEHAG2+80ecFqcnNzYXRaIxyf5qamqJcov5m+fLlWL58ueZw2aQihcbSROYIdbh9CIRYRVi69wjQ08jmkQ07TTlm5GwmhOq3oWCWDUapL1BTtwdFGeGdorkjlJsWLQa5EIo3ooMgCIIgjhV6VT6fkZER98PhcPTb4iwWC8rLy7Fu3bqw7evWreu3qjU9Fi9ejJ07d2Lz5s0Dep1eIwuhHgBAczcTLoIApNRLay2azAa0cgpOZY8N22E0CCiUqsCOdESHx/j5ctOjhZDdYoLNzH5k2noSD48danXii6q2hPcnCIIgiMGiV//SP//88/2+gJ6eHuzfv1/+uqqqCtu2bUN2djZKS0uxZMkSLFy4ENOmTcOMGTPwzDPPoKamBjfffHO/r0XNUHeE7KILAkKycEmzmGA48iXbZ3hEB+3Cieyx4xDg6cKwrBTUdbhxuN2N8rLwXVud+qExtt2Kug43Wp1elObYE1ryDSs3o6rFibV3zMbYgvSEjiEIgiCIwSDpsY0tW7Zg3rx58tc8Ufnaa6/FypUrceWVV6K1tRUPP/ww6uvrMWHCBKxZswZlZWV6p+wXFi9ejMWLF6OrqwsZGRkDeq1eIQkhA0JIhUcRQjYTEzoAkDs2/Bh7NpBeDHQfAZp2oiQzBV8AqNNwhFqk8+VpOEIAkJVqRl2HG+0xRnSoqe9040Azy2f6rKqNhBBBEAQxpEi6EJo7dy5EqeuxHrfeeituvfXWQVrREMdsY1VgQS8ccKFZyulJs5qAzsNsnwyNKrfCCUwINWxHceZsANqhMSVHSNsRyk5lAqk1wdDYtpoO+fOvajuw8FsDK2AJgiAIojf0acQGkWSkCfQZgjPcEeqsZc9nlEQfUzCBPTbuwLAYvYSU0Ji2I8RDZok6Ql/WKOM4ttV2JHQMQRAEQQwWJIR0WL58OcaPH4/p06cneynRqCbQy8nNZi/g6WTPZwyLPqaQC6FvMCxTEkIxQmNaydKAMqusNcFeQltVjtCB5h50efz6OxMEQRDEIENCSIchWzUGqCbQO+XQWKlRqsqyZbKO0pEUSAnTjTtRnMFETl27OywsGQqJaImXLC2FzBKpGvMFQthex8RZitkIUQS2H+6MexxBEARBDBYkhI5FVI4Q1zHFQiv7RCs/CACyRwGCEfA7MczExIhTNZ4DAKpbnfAFQrCaDCjMsGmehjtCiYTGttd1whsIISPFjLPH5QMArn72c5z9x4/R6SZniCAIgkg+JIR0OBZCYxlSd2kAKBBb2Cda+UEAYDQBjmJ2uKsevKG0y6e0B+DuzfhiB8xG7R8N3lQxXmjsk73NuPEF5qZNH5GNqaWZ8nMHW5zkDBEEQRBDAhJCOgzp0JgU+kqFkuOTF2xin+gJIdVzQudh2C2sYDBMCEniZOIw/XYBcmgsjhC6763t6HD5Mb7IgfsvPAXfmz4c3y1X1tbU7Yl5PEEQBEEMBiSEjkUsaQCAVEEREwWI4wipn+s8jBSLEQCbVcb5ui6+EOKhsVhCKBAMyYnYK6+fjrKcVDhsZvzxism4bCpL5G6SkrIJgiAIIpmQEDoW4UIIKiEkNrNPEhRCdi6E/CxHKBQS8Y0khCaVZOqegidRd3sC8AVCmvuwoayAQQByImaW5TvY141d5AgRBEEQyYeEkA5DOkfIGu4IZaSYYXPVs+f0kqWBcEfIzIQQD40dbHHC6QsixWzE6LxU/VOkmOX8Ir2Eae725KRZo6bb56fbwvYhCIIgiGRCQkiHIZ0jZGFChTtCpw1Ph9B1hD0X0xGSRFJnrRwa40Joe10HAJYobdJJlAYAg0GQJ9bXtrk09+G9jfI1ehHxbc1dJIQIgiCI5ENC6FjEEp4sfVYRgFAAEAxAeqH+cVqhMe4ISfPAxhXGnwU2Jp85UvuaejSfb44xr6zAwR2h3oXG2pw+BILaoTiCIAiC6CskhI5FuCMkMMFxWp4kEGyZgMGofxwXQu42ZBhZHx/uCPF+Qhkp5riXl4VQo7YQ4iInL03fEWrs8sadMcfZUt2Gby39AOcu+wSVh9oSOoYgCIIgEoGE0LEIzxGSHKGTMyUhlJIV+zhbBmB1AFAaMLr9TAi5vOwx1aqaw9uyH/j3ncCTpwO73pE3j5WE0P7m2I4QT4xWw7e5/eHNHPUQRRG/WbMLvkAI1a0ufO8vn2Gran4ZQRAEQRwNJIR0GNLJ0pIjVJgSxG8umwB7oIttjyeEANkVKpLK7d0+Jkac0iMPmcHvBp7/NrDlOaBlD/CfuwAfC5+NLZCEUGO35iX42A8tR8huMSFdEluJJEy/v6sJW2s6YDMbcNaYXARDIpat2xv/dRIEQRBEApAQ0mFoJ0uzPB6H4MHVZ5QBbskh6YUQyguxcnseGuOPshDa+y7gbAbSi4DMUqCnEfj8LwCAMXns+kc6PejWGKLa1MVzhLTHdOT1ooT+r+sPAgBuOHMkll4+ESaDgA37WihERhAEQfQLJISORSRHiDs0ihDKjH+sJIRyg40AFAHk9HJHSAqNbX+dPU7+PjDvPvb5p48BPicy7GY51+dAszLmg8MdIa3QGKCqHIvjCPmDIWw73AEAuGLacAzPtsvdqZ/6+GDMYwmCIAgiEUgIHYtIOUII+YGAt3eOUBqrKnMEOwAoVWM8VyjVamTn2/ce23/iFezDMQzwdAK1XwBQEqb3R1SOiaKoOEIaoTFAVTkWp4R+b2M3fIEQ0m0mjMixAwB+cHopAFCeEEEQBNEvkBA6FjGrGh76nIC7g32eiBCyZwMAUoMsr8jl13CEdv8HCPqA/FOBglNZJVrZmez4ms8AKAnT+5rC84ScvqAsqrTK5wHFEYpXQr9DNfJDEFhjRp6f1Or0oT3OvDOCIAiCiAcJoWMRowkwsaaG8Hb3zhGShJA9yESGOyJHKNViAg5vYfueNF85rvRb7LFmEwCgLIeJscPtyuBXQAl3pVqM4RVoKnh36cY4jtB2jdlndosJwzLZa9erWtOj0+3HhU9swC/f/qZXxxEEQRDHLySEjlXUeUK9EkI5bFd/BwBl1pjsCFmNQPNutm/+eOW40hns8fAWIBhAmk2aXh9RAt8kJUDruUEAUJjBhBAfzKrH9sNMCE2IGAI7WicsF4/3vmnAjrou/OPzQ/AGgvEPIAiCII57SAjpMKTL5wElT8jX0zshlMIcIauvAwBzgkRRVBwhsxFo2sX2zRunHJc3jvUh8juBxu3MOQILhamRS+djCKHReYqQ0Wuq6AuEsKuBhd0mlYQLoTF5fRNCH+9hlXL+oIjd9dql/wRBEMSJBQkhHYZ0+TwgT6CHrwfwdLDPbZnxj5McIbOvA4AIty8IXzCEQIgJktRAKzufYAByT1KOMxiA4Tw89hlzjgC4fOGOUEMnc4TydUrnAWBUXioEgYWqWnq083x4orTDZkJptj3suXgjPrTwB0NYv7dZ/vprqRqNIAiCOLEhIXSswoWQt5eOkJQjZBADSIcbLl9Q7ioNAPaOfeyTrJGAOULMDD+dPdZVyo6Q+lhACXeVZKXoLsFmNmJ4FhM3eq7OASn/Z1yRQ06U5nAhdCABIdTl8ePsP36My1dsRLcqjPeVFHYjCIIgTmxICB2r8BwhbxcrawcSE0LmFLnqLEvohtsflLtKW00GGFv2sP3yT4k+ljtEbQflxovOCEeIJ08PiyGEAFX5vU7CMz+PlqDix9Z1uOXcJj2+qu3AwRannHjNQ3ZH4wglOiONIAiCGPqQEDpW4TlCnXXKtkQaKgKyK5SNbrh9wfCu0s0a+UGc7FHssa1KrgiLcoRiCBg18rwynTEdR7izlBl9nuxUC7JTLQCAgxoNHdXUd4aX6C85j4m5/U09cUWUFq9XHsakh97Dxv0tvT6WIAiCGHqQEDpW4aGxzlrp63TAGH9yPABZCGUK3XD5AuE9hJqkijEtIZQ1gj2625AWYgLG6QuEOSSH210AgJIse+TRYYyO4wjxEFuxhhACgNF5zNU62BI7PMZzli6ZUoz/3HYWvj99OAodNoREpU9Rb1izvR7dngDe39XU62MJgiCIoQcJoWMVLoS6JEcokbAYR0qYzkY3QiLQ4WLzwlItBlXpvIYQsqYBaQUAgDQXE2AhEfAGQmwpHj+6PExUDdMRMBy9ztQc7gjphdh4jlG8EnzuCI3MTcWpxawxY3kZe69WbzsS81gteAPJA73sYUQQBEEMTUgIHavwHKEOyRFKydDfNxKphD5LYDd1XvKebfYrFWhZI7WPlcJj1u4aeVOP5CjxsFiW3RzeTLFhO/DNakDlHHEh1NjlRVfE4FZRFOVz6TlCPPQW2dAxkoZO9nxRhpL4fd2ZIwAAr22pRW2bK+bxaly+gHy9eE4UQRAEcWxAQuhYRc4ROswe++AI5RnYzbxFEkLFJjZ2A5Y05fyRSALJ0K4kTPM8ISXBWRUWO7INePY84LVrgfcflDc7bGYUSENZI/N8utwBuT+RnrPErxFPCHFHqDBDOc/0EdmYNTYXgZCIJz7cF/N4NQebnbKWO9zuhsdPTRkJgiCOdUgI6TDkGyry0FhAEgJ9EEK5RiZAWqVePoUGKWdGCn9pIidMV8uT6nnlmJIfJIkOZyvw6lXKGj99HNj8rHyqnFQmhDpc4b2EDne4pOctsJmNmstQHKHYjk6D1Ola7QgBwOJ5YwAA63Y2xjxejXqumigC1a2xE7UJgiCIoQ8JIR2OmYaKnF4JIRYay4lwhPKERISQFDJrO8gm1UNpqsjDWbKLs+0fLIcpezRw1p1s22dPyafiYzp6Iqq3jnQw8RKrBJ87QnXtbt1ydrcvKOc/FUYIoZML0gEA7S5/3HEbXx/uwO2vbpU7U3MONPVNCB1ud6GmNfGQHEEQBDFwkBA6VrGkhn/dByGUDeZwcCGUiw72fFq+/rEqISQ7QlGhMUnAHPyYPU5fxISQwQS07gdaDwAAHJIQ6vaEC6E6yeUpztAXQoUZNhgElqjNc5wi4W5QqsWI9IgBsJl2MyxG9uOv192a8+Db3+Bf247gX1JytdHAGjwe7EPCtD8YwiVPfooLntggD7wlCIIgkgcJoWOVyBweXtqeCFJoLENKlm7pZkIgW2xjzycSGutpYMnVUBwhHtIqybIDAZ88qR6j5rI5ZXxw6773AADpNlbu3x2RLH2kM74jZDEZUOhgLo9enlC9VFFWmGGL6k4tCILcXJEPitWjMaIX0YxR7P17dN1enPX7D1F5qC3m8WoOtbrQ6vSh2xPAXp0eSgRBEMTgQULoWEUdGjOlAOMvTfxYqWosUwx3hDKC0qiOWI5QShZgZRVqw42tABRHiPfsKc5MAQ5vBvwuIDVf6VJ90vnsce9adhnJpemJcoRiV4xx4iVM80TpIh1nKZcLoW5tR4njSAnvz3T+qYpQPNzuxjtf1cc8Xo267H4PCSGCIIikQ0LoWEUthE69LPGu0oDsCKWHOgGIaHUyR8gRSMARAoC0PABAgUFpqggoIa50m0kJi42cDXA3ZqwkhA59Cnh72H6A3HuIw3sDDcvUH9wKxE+Y5qGxyPwgTn6CQoi/PwAwKjcVM0bnhD0fmeMUC7UQ2ttAQoggCCLZkBA6VlHnCJVf27tjpRwhI0JIh+KmpAaYwxNXCKUyxyhXSq52eoMIhkS5saLdYgSqPmH7jpqjHJc7FsgoBYI+oK5SN1m6WRIm+Y5EhZCeIxTdQ0gNF0LNMUJjoZCINkkIrV58Jt64ZSZG56XhB6cPl9sHtDlj5xipUSdYkyNEEASRfEgIHatklgIj57CQ2PAzenesyQYYWLgnTSWE7F4uhGKExgDZEcqRkqtdvgDcqp46drMRqP+afcHzggDmDBWcyj5v2aubI8SFRY40T0yPeKExHqor0BFU+elseyxHqNPtRzDEqtLGFzmQlWqBIAhYevkkPHblFADhjlE8whyhxm58VduBLdWJ5xgRBEEQ/QsJoWMVgxG49m3gey8ooadEEQQ52TpVYCLCgBCsvgRDY6lMCGWKiiPEE6YFAbB5m1jvIMEYncSddzJ7bN4tV42pHSG3LyiLquw4QojnENXrjNmoamHuS1mO9twznizdHEMIyWFDmwkWU/ivS04aO75Vp2otElEUw4RQY5cXlyz/FN99elOUGCQIgiAGBxJCJypW1keHh8ay0A2DGAQgAKm5sY+VQmMZwQ4AzBHi3aXtZiOEtiq2X0ZJ9CBYWQjtkZOl1eXzrU4mKsxGQX5eDy6U2l3RIsIfDOGQ1KtndJ52l+xEcoS4yMmVRI+a3DSLtE9ijlBzjxfdngAMgiLCOA2dsSvX1Lh9QTzxwb5ejQchCIIgtDkhhJDJZMKUKVMwZcoULFq0KNnLGRpYmBBKkxwhuZmiPSf+FHspNJYWYFVmTl8QLqknTorFBLRLQoiX2qtRCSElNKYIoXYnEzXZUggKAODtAdbcDXz8e6BDmXGWlcqO73D5opoqHmp1IRASYbcY9XOEHFwI6YsQ7ghpuVPcEXL7FUdMj0ff24Nrnv0cAFCabceEYkfY83q9kLR4YVM1Hl23F//vre0JH0MQBEFoE/tf7uOEzMxMbNu2LdnLGFpIjlCm0QuEgHxBKp1PL4x/rBQaS5WqzFzeANx+JgTsFiPAHaFsjcGtuSexR2cTHFL5vpYjlJ2qckw+fRz44hnl80XvAwXjkWVn4iQQEtHtDcBhUwQcD0GNzkuL6iHE4TlCLT0+BEOi3ChRDRdCOWnRQijVYoTVZIA3EEJrjw/2bO1fp26PH098uF/+ekRuKk4tzsBHqk7V8Zo6qtlcxd73T/e3oKnbI78OgiAIovecEI4QoYEkhM4czgRHHvh4jTiJ0oAcGkuRcoqcvoDsCDEhdJDtp+UIWdMBRwkAIMvFBJM6PyYqUdrVpozlSC8G/E55XpnNbITNzH6EO5zh4bH9TVwIRXTgVpGTZoEgAEFVZVgkPDSWoxEaEwRBDpm1xHB0IkNnk0syceNZI3HXeSfhjJGsgq8lTgk/RxRFfFnDRGtIBNZ8nXgPI4IgCCKapAuh9evX46KLLkJxcTEEQcDq1auj9lmxYgVGjhwJm82G8vJybNiwoVfX6OrqQnl5Oc466yx88skn/bTyYxwpWfosLoSEDrY9XqI0IIfGbF6loSJvqmi3GJXQWJaGIwQAecwVSu9mgskbCMEnld63RYaiNj4B+LqBwonApcvZth1vAAEmHLgr1B4xuJU7QmPytfODAMBsNCBbOl4vYZqLGL0KNr7OWCX06vDaMwvL8eM5o5CVasFPzxmLU4pYiCyWkFJT1eIMy4n664YqXPH0Rny0pymh4wmCIIhwki6EnE4nJk+ejCeffFLz+VWrVuGOO+7Afffdh61bt2LWrFmoqKhATY2SK1JeXo4JEyZEfRw5wmZDVVdXo7KyEk8//TR++MMfoqura1Be25BGcoSG2QO4de5ozC2VprxLXadjIoXGjEE3UuCB06cOjZlUjpCeEBoHALB17JM38cqxqJyc7a+zx9l3s3YBjmGApwPY818AQKauEGIVY3qJ0vJS0mPnCcUr5c9JIGGan2N4th3zTy2UZ7QBSsJ1okKo8hBzg0blpkIQWPPJzdXtWPlpdULHEwRBEOEkPUeooqICFRUVus8vW7YMN954o5zk/Nhjj2Ht2rV46qmnsHTpUgBAZWVlzGsUFxcDACZMmIDx48dj7969mDZtmua+Xq8XXq9yUzpuRZOVORGCtxv3fHsc8LYJqEdiHaotaWysR8CNXKETLm+mHBrLMToBjxRm05t/JuUJGVr3wW6ZA5cviB5PANmpFrSpHZieJqCzBoDA5pUZjMCkK4H/LQO+ehU49VJk2XnCtOKSiKKIAzw0FsMRAljTxt0N3bqVYy0xQmNsnVJozKkvZNp43pM9OgldCa0lliPEw2LnnVqANIsJT3y4H75gSO7G3Vt21XdheLY9boUeQRDE8UrSHaFY+Hw+VFZWYv78+WHb58+fj40bNyZ0jvb2dlnYHD58GDt37sSoURq5KxJLly5FRkaG/DF8+PC+v4ChDB/R4ZP62nDxYsuIf6wgyOGxPHQyR0gSQmVoYPukFYZ3v1aTIb2n3fXyDbhLyhNqk5yd7DQLUCcJ3NyTlHWdchF7rP0MEEVkaYSmmrq96PEGYDQIuj2EOIVS5VjkYFVOPEcokRJ6xeXSKsGPn2OkhjtC5aVZ+Ok5Y7H2ztkAgCMd7qjKuXi8tKkaFY9vwB2vbu3VcQRBEMcTQ1oItbS0IBgMoqAgPG+loKAADQ0NCZ1j165dmDZtGiZPnowLL7wQjz/+OLKz9cM/9957Lzo7O+WP2trao3oNQxYpNAavNOahN0IIkMNjOUIXXD4lR6g4JH1ftBKlObwyrbtBnjfGQ2NyjpBdJYRKVO5d/njAYALc7UBXncoRUoRIjdRfZ1hmCqwmY8yXUSgNZK3XGbOhVI3pOEKyEIrhCPXoV57xwa+JJEsfaO7B3sYeGA0Cpo1gP8O8NYDLFwxzxeLxZU07Hv73TgDA+7uacKjVGecIbfzBENp70VmbIAhiqHFM+OGR5c+iKOqWREcyc+ZMbN+eeL8Vq9UKq1X7pndcESmEvFIIMGEhpMwbCwZFOUcnR5TGdGQM0z82vYg9ulqQJfVu5CX0YcnSW7ewJ4edphxrtgG5JwNN3wAN25FlZ4JLnUDc5WafZ2mEoiLhQkKroWEwpHpdGiIGUEJjscZsyC6Xhquk5Aj54v5cv1F5GAAw56Q8+Vw2sxG5aVa09HhR1+GWHbJYdHv8+OnLW+EPspYBwZCIV76oxS8qxsU9NpKH3vkGr35Ri9WLz8SEYQn+7BAEQQwhhrQjlJubC6PRGOX+NDU1RblE/c3y5csxfvx4TJ8+fUCvkzSkqjF4I0JjVof2/pFIobFcqeyeh3YyQrwxY4zu1PZsedZZiZn3EmLiRS5XTzUBdV+y/YdF5HMVTWKPDds1k6V5mM2REiGEepqBXf9WkrmhTKav1xBC7S4fRJFFAnl1WiTZvUiW1hZCTEj5giF0efSbMgZDIt7aWgcA+M5pJWHPDctkryHRPKGH39mJug43hmen4JHvsPfy9cpauXKvN6z9phGBkIiPqWqNIIhjlCEthCwWC8rLy7Fu3bqw7evWrcPMmTMH9NqLFy/Gzp07sXnz5gG9TtLggoc7QX0MjRUYmZDh5efpXAil5ugfKwiyK1Rs7ADAQmN+lRjI9dUB3k42IJYPauUUTmSPDV9rJktzd4mH3QAAG5YBfxwDrLoaeH4B4A+fTN/QGS0i+IiO/HSrZrNFAMiVHaFYydL6eUY2sxHpUp5UrDyhTQdaUd/pQUaKGeecEt7raVgWC+8dSUAIVR5qx2uVhyEIwKNXTMHFU4qRn25FS48PGw+0xD1eTVOXR/6+76g7TosKCII47km6EOrp6cG2bdvkzs9VVVXYtm2bXB6/ZMkSPPvss3juueewa9cu3HnnnaipqcHNN988oOs6/h0hKTTm6wFEse9CyMQcJe6opAUTcIQAIJ05eoVSR+tuT0DONTEIgKNjN9uvYEL0yA9ZCG3X7CPEQ2Nyp+mAl3WkBpgT1V0PVL4AAChypEjH++GRhr1y9jQwkXdyob5Lpi6f51PqI+FukV7YKpE8oa1Stdg54/JhM4fnPRVLeU517fGF0JdSsvX88QU4fWQ2zEYDZo1l38vN1W1xj1ezva5T/vyb+s4Ye+rj8QdxyZP/w71v0rgQgiCSQ9KF0JYtWzB16lRMnToVABM+U6dOxQMPPAAAuPLKK/HYY4/h4YcfxpQpU7B+/XqsWbMGZWVlA7qu494R4lVj3m4g4AGCkpBIVAhJjlKWgQkg7kak+DvY8/EGt0oJ03kqIcTzbLLsFhi6pCR1rRL8ggnssb0aOWZ2fXXCbpQjtOe/rPdQejFQ8Tu27dPHAL8HjhQTUiRhEZkntKeBuRzjCtN1X0aBwwaryYBASNQdgppo5VmsEno+i4y7P2r4tkRCY9VSUvTYfOU1nT4yCwCwuao97vFq1EKots2Nzl4ka6vP8dXhTvxzS99CcwRBEEdL0oXQ3LlzIYpi1MfKlSvlfW699VZUV1fD6/WisrISs2fPTt6CjxfUydLcDRIMikCKh40JoQwDu/kGJDckxS/dTOM6Qiw0lhPiQsgvi5nsVAvQyRKDkanRvsCeLZfg5ztZU8awZGlPhCO07WX2OPn7wNSFbMRHdz2w818QBEEOj0XmCe3mjlCBvhAyGgS5aePexu6o592+INyS06SVIwQoeULNMYa/8hBU5NR6gFXHAYmFxrgQGpGrtDaYLlWgbTvcAW8gqHmcFjvqwl2gvrhCVS1sPcGQKK+NIAhiMEm6EBqqnDChsYAHcEmVXlYHYEjwR0I6Pg3hLohVmj8Ge4wcIUB2hDKD7No93oAsZrLsFqBDcoQySjQP500ZHW4mmNz+oBza6lI7Qj1NwP732TFTrgJMVmDC5ezrGtaLiidMN3QpQkIURUUIxXCEAOCkAiaE9klNHNXw3CGL0aDbtDCRpoq84WOeRhl/cWYvHKEW9v0aoeqvNDI3FblpFvgCIWw/HF/MLP9oP8783Yd4fxdLkC6QejF904c8oeoWRfzsa4x+/wiCIAYaEkI6HPehMavq5t7JqpG4y5PY8Wxfu6gIIQv8MPmlm1nc0BhzhBx+lqDb7QnAKfUSSrOZFEcoo1T7eKk8P8VdLycy84RpniOUbjMDhzcDYhDIPxXIHcuOLZHE7WFWnq9VOdbY5UWn2w+jQYg5rwwAxkqOkZYj1C4Ng81OteiWxvPrH27XDq0BiiOU74gWQiVSaKylxxeV56TGGwjiiJQUXpajOEKCIGBaGXOFvkggT+j5T6vCRNflUhXb51WtYQN0E0HtAu3XEJKJIIpin8JyBEEQAAmhExejmVVkAUCn5L4kmh+k2jclpNzIMiHdyAQjYMuMfbw03DXV1wwA6PEE4PTxeWVGZU16jpAUGhM6DyNTKpPnCdM8R8iRYgYav2H7F01WjuVCqGkn4O3W7CW0W8oPGpmbGpWcHMlJshDSd4T0wmIAcGoxE5Vf67gxoigqobE0W9TzGSlm9p4htitU2+aCKAJpVpOcl8SZNoLlCW2pjp8n5A8qSeGnFDlw+kgmot7f1YTTf/OB/N4lQlWLIv72NUULyUT455ZaTH74Pbzz1ZE+HU8QxIkNCaETGe4KdXFHKLMXx7KbtyXohAEsyTVHkG6A9uz4ITbJEbJ5JCHkVRyhLKOPJTcD+kLIMUxee6Y9XAjxHKF0mwlo3MH2U5fgO4qYkBJDQN2XSndplRDaoxUWE0Xg0z8Dy88ADn4ib+ahsQPNPVGVY7F6CHEmlWQCAA62ODWdjR5vQM4zyk2PPo8gCBiVxxyeWOElHhYry7FHuVPjJTFW1RI7TycYEmXX596Kcfjz96fgrDG5uPqMUmSkmOH2BxMSUwATeOrQ2I66Tlz//BdYumZXQsdzPpBCdP/5ur5XxxEEQQAkhHQ57nOEACUxmoeheuUIKWG0NDAXIlsWQnHCYoCcI2T2dcAKH5y+AJzSvLJioUVZj164jneu7qyTS+h5aEx2hGwqRyiyFxEf23F4M4oc0Y4Qz/cJS5R+9xfAuvuB5t3AGzey/CMAw7PssJkN8AVCUaMqEhFC2akWlGaznJ2v6zqinuduUJrVFDa5Xs3JBex94gJOC61Eac7wLHb9ug43QjptAAD2ekJSk8kbzxqJsQXpMBsN+M1lE3HhJCZu9QbYRtLY5ZUFHlufCx/tacZf1h/s1dy0/c3se/XV4Y6EjyEIguCQENLhuM8RAhRHSM4R6oUQMlkBI8tXSZcSpnMg3YTj5QcBQEqWfHye0AmnNwCX5AgViFKXYr38IIBVfgFAV50cFnJJQkruI2T0Aa0H2H685J5Tcjp7PLxFlSytCKFO6Rw8kRkN24HPnwYgNYN0NgNv3wYAMKjyiCLzhFoTEEIAMHl4JgDgq9qOqOfk/CCNijEOL/Hf06gflpKFkMYg2sIMGwwC4AuEYjZ2bJE7f1tgMob/+eAVbbGq3zj/+PwQXtxUDQAYnp0Cmzn8XF1u/S7bapj4ZD9/9Z0eNHZ5UNvmiinmCIIg1JAQOpHh3aX7kiMEyEIqXWCOUJYgiQC7/lBbGUGQmyrmo52FxiQhkxNg4TLdsBgAOIrZo68HOUZ243X7AvAGgvBK/Wgyew4AENlcNGkkiAx3hI5slcvseWhO/XmqVcoP2vMuezzp28A1bwIQgL3/BbrZ+JeTpL48B5rDHaFa1QDYWEwuYe/9ttroPCHusOTGEEIncSEUwxHigkGdKM0xGw0okkKEtTEaM3JRlqtRvZafbgvbR48DzT24760dWPExE6mjctPkFgTydWKIMTWHWp1h4cibXtyCWY98hNeluWwEQRDxICF0IsPnjfEcoUTnjHGksBV3hHINXAgl4AgBQAoTTJlCDzz+kOzkZAca2fNaPYQ4Frtcos9DaS5fUA6LAYC9nXenPjXqcOSMYY89DbAJkpDyB+WQjCyEeChqrySETv42UDBe6W5dtZ6dTko+5q+Bo4iPaBdGzRTuCB3uiAoLxeohxOGOUHWrS7dyTHGEooUQkFhjxlhr4Y5VvNBYZOPKogwbppZmhm2L5UqpOdAcnhP1lZRw/uJn1Qkdr6Y34TiCII4fSAjpcELkCPHQWEgSD712hKTu0iZ2Y8s39CI0BgApmQCADLAbNL+BZvikIbuxHCFATpjOR7QQSreaYGjSyQ8CWGhOWn+qm1UbiSJkN6lHdoSkXkR1ley4sfPZ46g57FFKmua5O7zyjZ1PSQYeqZGXo+bU4gwIAhMabRGT7Lk7Eis0lp9uRabdjGBI1CxDD4ZEWYCUaHSnBoASybWKVcbPBYpWPyMujpq6YouYyNc3qSQTd58/Ds9fPx3Tylj1Wqwhtmr4a7Wawv+UaTlWeoiiiIV/+xznLvskoaaUBEEcX5AQ0uGEyhHi9FYISY5QsY25IHm9dYSkKrUsI7vxyknBHi6EYjhCgCyU8oIslObxB1U9hExAyx62X/746GMFAchkY1ps3UoYhecZOb3sMc1qAvatAyCyEnwekhs5lz1WfQKIYlSeEsBu+N3eAAQBGJ6tcoRaDwDNe8OWk2Ixyu5T5BR6LixiOUKCIMiJ3VrhsZYeL/xBEQZBX1BxgXQ4kdCYliPk4I0hvTFzdHh136SSDDx40XhcftowZKSYMe/kfPk1JuoIcSH07QmFmutMhKoWJzbsa8GBZid++NwX+Kq2dx22CYI4tiEhdCJztEJIclQKLFICLc8RijV5Xo3kCOWa2I2X37xSvFLVWHqh1lEKkiOUIwkhtSPkSDEDXVJfGT1nKYsJIWNnDSySo+CSHJ2wHCEp/IWx5yvHls1gA1w7a4H2KtilrtEub3gVFAAUOWxKL6KW/cDTZwErzgD+9ydmQ0nwfCR1rhKgOEJaLoyak+WE6WghxJ2OQoctKsmZU8Irx2IIoViOEHdhAiExbAhuJNwRmjAsA9efOTKsT5PSZTtBISSFxiomFOLbpxbKIchEK9cA4NMDrcr5mnpwyfJPcfYfP4nZnFKP+97ajorHN1CDR4I4hiAhdCJTOjP86147Qmz/XDO76WTyqrF44zXk4zMBADnSvDJfkIWlzH4pYTglTtK1VEKfGWBVZi5fMLyHkJTILLs4kUiOEDoOyY6OR8oT4iGuVKvKWVI3ZbSkKo0ZD34Cu3Qzd6nLwVsiytVDQeBftwJ+F+th9P4vpUo0KNcCwvKcgMRyhABgrFS5drA5uhfQkQ4WFiuOkbStOEIxOlz36K/FbDTI1XGxhIg8U84eXUmXk8AAWo7LF5Bf69iCdDy9sByv3zxTOt6LQDCxIa6f7mPC+/KpwzD35DyYjQLqOtzyiJVEaery4OUvarCrvgv/3dG3nkadbj8W/u1zrNpc06fjCYLoPSSETmTGnhdeot6bERuA7ChlS1VbGehFHyFAdoSyDMqN24AQjF5JCMWrPpNCZxk+JoTc/oDc7C/PEgC80nr0nCU+2b69WhEyviA8/hB4ZCfVYgRa2GBXeUQHZ7gkhBq/kd0cl8rN4T2F5CqtrX8Haj8HLOnA6TdJ2/4h758uCaEoR0gun4/uKh32clK1E7YBxRGKJYSGqUJjeonDsarG2BqtYftp0cZnymm0FEjUEaptc+GCP/8PLl8QmXaz3IcpJ9UCo0GAKCYmpoIhEZsOMkfomhllWHn96SiX8pQONkfnWsVi7c5G2eBbt7OxV8dyPtzdiA37WuSKOoIgBh4SQjqcEMnSBiMw7Xrl6z6GxjKNbgAi0kLSjUMSOHGRHKFMQXEg0uGCADHseV0kgWP3KcnSvP/MMFMH28eSHh0C5EihMbQfgk2V49OjEiJ2TxPg62FjQ7JGhh/PK8/aDiBFyu9R5wjx0Jjct2fPGvZ41u3AnF8AggFo3A60VwNQHCF1wnUgGJLHdMRzhDKkUSNdGvO+6hIQQkUZKRAEljCuJyL4dr215CVQOSY7QqnmqOcSFUJ/+18VqlqcKMqw4ZmF02CWwn0GgyCH7ZoS6Gf0zZFOdLr9SLeaMGkY+/kflafvrMXiXZULtGF/S5SgTYQDTeyah1pdYT+HvUEUxT5dmyBOVEgI6XBCJEsDwNSFyuepefr7aSE5SA7BDRt8MKGX1WfSfg5B+c87k39uSQNMsZsQ8tCZVQqluX1B2REqNkhjHmLlGalDY1JDP7c/qCqdN8LQtp/tkz0yej1cCLXuZ84RlBwjIKKTcygIHNrEnhh9DsujKjuTfb37P+x6khBS3wDbnD6IImAQ4jdl5P2QYjlCwzL1XSWLyYBCh/4AWH8wJOf3xBdC+iKEnyNLIzSWK4fGEivBv3XuaHnWGadAStpujFO9BgAbpfygM0blyLlTo6RQZrxxI2ranT58dpANrM2ym+ELhLBhX3PCx3PUFX+xekLF4pG1ezDpofewtSaxUScEcaJDQuhEJy0PuPlT4KaPWd5Lb5AcoWEpfswcJv13LxiU0R3xkJyjdFG54WTxwa3x8oMAVgIPwOLvhIAQ3P6gXHGVD2mKuqNI//hMKSzo7UKeiTdlDIaXzrdI1V05Y6OP50KooxZ2IzuGO0KiKMo30hE5qWzmmbeTOVSFk9hx4y5kj7v+DUCqUEN4aKxJFYoyGqT5YKEQsPc94Mg2IKjs65AcoU4tIdQZ3xECIHfZ1hIRvKTdZBDkQbeR8PBdrBJ6nkitJey4IxSvfF4WUxrnyHfw1xDfEfryEBMLp4/MkrfxuW2RPYr0cHoDuPv1rxEMiRhf5MB3TmPJ+Wu2NyR0vJr9qmvuqk98eK2aDfuaEQyJ+GRv74UYQZyIkBAigMIJQPHU3h8nOUL2kAvPff8kts3qYKXpCR2fCQBIDSl//LN5vlAi4TVJCAliCGnwhCVL54QkIZQeQwhZ7KzrNIDhBiXh2qkWQq2SIxSZHwQwB82SDkBEhrtOPh5grg5Peh6enQJUf8qOKT0DMEpNGsddwB5rNgHuDjnPqEdVeaaZKP3ZcuDlK4Bn5gB/ngK42c2ch8acvmBUonC9lCzNu0frESs0JY/XSLPAYND+Hss5QjqOjiiKsR0h6XiXLxjmrkXS5tJPuOaOUFMcISSKIrZKI02mlqqEUC4T8tWtzoRGddz+6ja8v6sRFpMBd59/Mi6azJLz393RkHD1G8AcN/Wsur4IIVEU5Saeu+v75igRxIkGCSGi7/BO1N4uwCMlOPcmz0gSO/aQ8gc73yyFZFKyNA6IwGwDzCz/JlPolkJj0gT7oFQSHUsIAXKeUInIE66DsphJtRoVR0hLCAkCkDOa7dtTBYCFxliOBjuHySAgxWwEDklCiIfDANY521ECQASadyuhMVXVWJQQ8ruBT//MPjeYWPn+1/9kL9WmDGRVV555/EF55lm8UR+xhFC9FI6KlbQtzxvTcYTcfmUEipYjlGoxys0RW7r1XSE5zyhNwxHirlScEvojnR40d3thMgiYUKz83JZkpcBsFODxh2QnTQ9RFLFeCoE9d+10zBuXj8nDMzF5eCZ8wRBe/SLx6q+aNhf8QUV47azvwvs7G+UxLYnQ4fLL33utNgoEQURDQojoO1wIebrYB9DLCfaZAABLyAOzlF+ULzVXTGheGSALpkw44fIF5LBQup/3IoojhKTKs3ww4eT2BZTQmMWkqhg7Sft4SQildB8CAISk7tTczUixGCEAwKGNbP8RZ4Ufn3cye2zapVk1FtVVets/AGcTW/d5v2LbKl8ARBFmo0FuA6AOj/H8oFSLEY4U7en18nIkYaEVmopqB6ABX2ejTo4Qd4MsJmWtagRBUMSYU1vIhFR9imI5QvFCYzyHZlxROlJUazEZDXKlX7w8oVanD75ACIKAsFyl62Yygf3SZ4fgT7CMn+cHcUG7taYDi17cgrtf/yqh4wHgkEo0Vbc68d43DfjxS1tk4ZgoKz7ej39uru3VMQRxrEJCiOg7vNze2wl4OqRtvRFCyr58zEaOkYfGEnCEADmXKEvoZjlCkgBIk0rqY+YIAfI4kAyROVrq0Fi2OaAMpNUVQixPyNxxUN7EwjqSq2QxAd31gLuNVZ4VTQk/Pv8U9qh2hFQhIR7ekR2hTSvY48yfAlN+AJhsQNM3QN2X7HVoVI6pewgJccKWOTEcoYN8XEiMuWm8g3Zdu1uzj0+7U5onZ7foroWHx1p0HJ1Ot19ub5CpIYSUHKHYjtC2mg4AwNTh0T9rPGE6XuUYDznmplnlppwAsGBiEbLsZjR2efH14ehBumo63X4semEzHvwXGwkz9+T8sOd5EnYiqENrogjc9FIl1n7TiL9uOBjjqHD2N/XgkXf34Bdvfh01DoUgjkdICOlwQpTPHy1yaKy7b0LIYASsbP8MqVosx8CTpRMVQpkAmCPk8Ydkp8DmkYRQuk4zRY7U/DEjxBwtt19Jli4xSOE1W4a+QyUJIUPbQdikyjOnNyCLKbvVKJfHI6MkuvIsbxx7bNqllM9rOEJ5aVbWILLtAEtIn3IVe4/GX8p23Mb6ESmVY8o56qXwDk+ElvF7gOr/Ae4OeVOs0FhVC/vexHKECh02WEwGBEKiLMDU8NwerSRneQ2psZsq8nOk20xh4oOTn0DlGgA5P4gPvFUzMo8LodgJ03ISesR7azUZcZI08iRWg0qAJTe/v6sJDZLoPbkgTZ65BgBGg4BgArlKAFDTqn0tjz8xVwpQqtVCIvDh7iZ8WdPe5wo2gjgWICGkwwlTPn80cEdIDAFdUg+V3vYiSpGEkOQIZQrcEepdaIwLqeZuLwSEYHZxIRRnTIckhNJDSgk+z+/JN0j/yacV6B+fzUJjaN0vD151+4Nyh2m7xQi0VUn7jow+XuUIaVWNKTlCNtn1Qd44pTfSKVLlWe0XACCHvtShMf55jlp8fPki8MexwMoLgFevkkd95Mbo7Fzdwm6ysQbIGgyC3NzwUFu0mxKrhxCHN3bUEyFt8jm0xVSB5Ai1On26M8NCIRE76tj3d0rE5HsAGJvP3t99GgNs1dTH6M/E87G0BKGayDDkmPw0/OGKyfjjFayTeTDOyBI1vHeV2RjutsVKPI9EnVu04uP9uHzFRlz1188SShwniGMREkJE3zHbWbgHADqlwaW9HtORCQBwSAIoA710hCSnhpfdh0QgG90QQpIQSFAIpQZVoTHpppGDDrZPTCEkiZueBmSY2X/dTm9Anjlmt5gUR4h3slbDc4R6GuXXrq4a4wm/+Q4rcEQSQsWnKcfzsR/Nu4CAVzM01iF1cubPwdsNrLlH6bx96FNg/wfsNes4Qi5fQHYsYgkhQGkgWa3hTsSqGONwh4Y7Nnrn0BNCOakWZKdaIIqQxU4k7S6fnLRdmh0d6jupgFWO7Y2TcMwTyLWq8bg4quuI7QjxRPYUsxFXlJdg7sn5GJmbiu+Wl8jiNdEhsjWS+Jw5Ory7e2+G0O5TvWYeGmx1+uR19oZQSEx41IkeNa2uoz4HQcSChBDRdwRBKXNvl1yPXjtC7HgH2M3CIUp/hBMOjbH95NwiAIW8mWJqHmDUdx7YPuyGYQ90AAhvqJgtqs4T6/oGdo1iM1u7WyWmUi1G5b2J7EwNMGfHwfrO5LjZfpqOUJpVcYSGqVodZAxnYjIUAJp2yaExtSPU4WY3sAwuPna9AwTcLKz3rcVs20e/BkRR7src7QmEDR3lblCW3ayZl6OmNJsJpZpWDUcoRg8hDhdC2+s64Qto5RnpJ0oDLOGah5a+qNJuKsgdryy7We5KrWZMfhoEge0XqwRe6dgdXUnHna14jlCblBT+o1kj8YcrJocNoU2kU7caXjp/6dRiqDsc6LUz0EKv2qw+TgVdJJ/ub8GsRz7CvEc/DmsU2Rve39mI2X/4CMvW7e3T8QDrzt6RoKNGnJiQECKODp6D07SbPfK8oUSRHKEMyRFK46X0CVeNsf3UQmi4VfoPPJFO2dJcNJu/AwAPjTEhkhli22I6QoIgX6fIyNbu9AXh9iXoCAFAPssTyuxh86V4jpLTG5CTrvPSLNqOkCAorlDD13JTRXV3ae4IyU0Qv3qFPU7+PjBrCWBOBY5sBWo2wZFigkUSBmoHoCqBijHOiNyjc4RG5qYiI4V1Z97dEN1LpzVGM0UOr+DaXK2daBxvkK3dYpKdoliuUEKOUHtsARHL4cpLYHYbAPzj80OYsfQDWTDNOzkf//zxDDx51dSEjud4/EG5OvDSKeH5dbz6MBHe3VGPq5/9HHUdbtS2uXHlXzYl3KBSzfu72My2f3/dtyG2AHDz3ytx+m8/6FUbAuLEgoQQcXTwye58UGofHSGeI2QPSje+XjpC2YIihAos0h/9eLPKADk0ZvV1wABW9s5DU+kByU1IiyOopOcLDDy8FpAdobg5QoCcMJ3axZo3ciHEb2qpFiNSXYdZ40SjBSiYEH58kdSpul4lhDzROUKZdjMLYVZtYE9M/B5zxMaex76u2QRBEJQJ8KqbJ0+UDguLhULABw8Dz8wFmnbJm3np+SENR4jnw8RyhARBwFQpb2erVNmlhjtCOTHOMX0EE0Jbqts0c1uae5RqLz14svPeGInCPEeoSMsRknOEYgsI+T3RWEsiI0sA4F9bj8iiLNViRKbdgmkjsnGa1CiypcerO0hXzcFmJ0IiC6P+/ruTsHrxmbhgIqu8rIvjbKlZubEaAHDhpCKcWuxAq9OHv65PvHKNs00Kj9a0ueImnWshiiI+O9gGXyCEz6sSr74jTixICBFHhyP8v8a+5ghlCE4YEIIt2NuqMamPkGpeWb7ZnfhaJOdJgIhM9MDtD8mOUFpAqhpLzdc7Ouz5PAMTcS5fUM4RyjJ5AZfU00jPEZJmnlldbCSDLxCCPxgKdy24G1Q4MbryrFByhOq/gsPGk6WV8JrsCNnNUi6QCAz/ljJ0dvjp7LGWFQZwIdTqVAshKVFaEjkIhYDVNwMbHmVu0mvXAT62T5nkpNS0uaJESI30X3lJVuzGjnKekMa8rEQqz04tdsBuMaLLE8Depmghw5s1xhpky/OE9jRqOxnBkIhG6XtUrOkIMXHU7Q1ojj3htMUQdrw5ZDxHR91z6ZoZZfLn/HvpD4ox18Dh7tfJBemwmoyYMjxTFeJLzBFq7vbiC0l0/KJiHG4/hzUj3aaT86WHyxcIc+M2SXPhekOr0yf/Y7G7jyNLADactzddwoljCxJCxNHhGBb+9VE4Qg44lcnzvUyWzhCVP5g50tywhMZ0GM3ymrOEbrhVbk6KT/oPMlZoTPV8DqITrgtD0ryplCz990ZK6Da5GuVNTm8gXAg1SzkSkW4QoITGGncgw8Z+pbu0coRSLEAj61WDkmnK8SWSEDr8BSCKSgm9qrMzT8It447Q/veBr1exZPmULKB5N7DuAQAsN8ZoYJ2Z1bkt6vlr8RKu+cgLrYTptjg5QgBrisjdkM0aTkBYWwIduCO0Tyc01tTtQTAkwmQQNAWV3WJClp05dLFERCKhsXg5Qrxdwn9uOwv3Vpwib7eajEwAI7HwGM8POqlQmRfIWwMkmiP03s4GhERgckkGSrLsmCyJ2r2N3b2qXttR1wW1jt50sPdC6JAqPLu7jy0A9jR046In/oefvPxln44nhj4khIijI8oR6nuOUJY8eT49fpIzRxJM6aJ6XlkvHCFAzhPKQVdYQ0WrV/rDm2BoLFvsAAC4vAE5R6gwKAkhrURpjtT92tDdKPfF6fEG0NytGmnBq/L4oFg1OaNZBZ/fhcIgy6XQqhrLtJsVIaQWVEWTWMjN1Qq0HZSFkDrBlp+D9/jBnjXssfw64PJn2edfvQIE/TAbDXJYSB0ea+zywu0PwmgQ5MaLeozNZzfiunZ3VEinPU7VGOfUYexn8YBGU0Qe9suN4QidXMiE0J7Gbs2wEk+CLnDYlIG4EcRzU9RdsrUdofg5QqIoysnAmvPb0hLLMwqFRLy7g/28ThqWKW8vkqvfEguN/VcaNlshhdQKHDYUOmwIiUzcJMpXkgjmLR02HWhNKLynpkbVwmFXfRcauzyyW5XwOg53ICQCm6vb5d9r4viChJAO1FAxQY42NCa5NlkGJzJ7WzoPyMnSaWIPBLAKI7kXUSI5QoBcOZYlzStzeoMQEILFLYW0EgyNZUjJ1S5/EE7pD2au/wjbRy8sBigl/t31cFh4U8agHN5ypJiBTmlmlTQSJAyDUT5/jp+5SjwMEgiG5NlTmTYT0LiDHVMwXjneZFU6Xh/erNlUkYcX0mwm1nNo33vsiZO+DYw+m30ffD1AXSUAadAsgMOqROGDUp5RabZdqdQK+oH37mfz00LKTYY3lwyERPgiSqcTCY0BqiG03mgXIhFHaFRuGowGAd2egGaX6iMxKsY4PGRWpyOEOlRdsrVeTyLJ0j3eAALSSbSEUJ6GsNXik33NqGpxIt1mwgWTlI7sieY6AazYgDs3FROU1hWTStjfha96ER7jobSrziiDxWhAfacHi17Y0qtWALzaEWBhsjN++wG+95dNmkn4evAmlcGQiO11nXhm/QF81gd3ihi6kBDSgRoqJsjRhsYksZIpuOSmiLD3Rgix4w0IIZ2X4Auu3q1FSpjOkcZ0OL0BZMAJQZRuoPGqz9KYEHIEWT6LyxuAi1ee+aRwl5aTIx8vhd5CfhRZ2c2mxxuQwwhpVqPiCGWUaJ9D+j5k+lkjSR4q6VINX80INLMO4IIRyD05/Hg5T+gL+T9wdaM/Pgg23WZmYqqrDjClACNnAQYDMHI22/HgJ2w5Uhm/OhSiGRbbsAzY+Gdg3f3AP74LeNnPQKpq9pfLG/5feFsCCdcA5NltPVpCKE7VGMDmoeXIXa61htBKidIa+UEc7gjpCSFeOu+wmTTL+BNxhLhbZzUZwmamcRKtPHv+02oAwJXThstCFFCq35q7vboNKjkHW3oQDInITrXISfMA5PDYtsMdMY8HWKjwrn9+hXVSxdi3RmXj3gXjYDII+GB3Ex59b0/cc3BqdCrF4o09UaOe33bXa9vw2zW78f1nPkv4eGLoQ0KIODoiZ3n1tnye5wgJTrkpYq8cIZOVlX9DcYJSQ/w8mYmdQxJCWeiGyxdEjy+AXKFTWUtkcnIkklBKl5Kr1bPG7FKjxphiymSRw3PDTWz/Hq+q8sxsADrr2L56QkjanuYJD43xkEm61QRTi9TiIHcsYI5wMYaVs8f6r6IcoVBIlOefpVlNwN61bN9RcwGzJAJGzWGPVUwI8RuySxVKqGqOEEIN24H1j7DPDWbgwIfA50+zt8RokKfQO1ViyqNy2+IJoTSbvhDiry1W1RigOCxcbKjhoTGtijHOsDgl9Fxs5uisg4sYtTCOpD1GWEx9jlhCqKnbg/V7myEIwA9njAh7Lstulr8XjZ2xxRRvwDgqIgeMJ79/nYAQevPLw3jjy8PwBUI4uSAdp5Vm4fozR+IPV7DqSL0+R1pUS6HZFHO4QOz2JJ6rpO6HVdumfB97E6b7YFcjfvDMZ3I/JQqxDS1ICBFHhzVdnhcGSzoL0/QGyRHKElwYmy7dbHojhAA5YZqH1uxcCPXaEWJ2uSgCeVwIxQuLAbIjZJeSq5kQkhKu/bwpY67moTJSntAwUweA8O7UuUInEPSyGWORoUhOBnOEbG6l8szjD6JDCpFl2M1KWCx/fPTxOdKokPZqti8UMeX0BfgEDjYZ/eDH7IuT5ivHj5SEUO0XgM8pT5YPE0KRvYjW/4E1ghx3IbDgD2zbrnfk/ZXZa8o5eGKxySDIFXJ6pFnZ64i86QWCIbkXUSxHCICcaKw14qJenjOm7wgVxwkrxeuSnWY1yTPs9IRMuzoHTINEhFCDVHpf6LChNGKoriAIiqCLEx7jvYJG5YULoYlSaKy2zS3neOnB2wBcdUYp/nv7LLnB5MhcljdW34syfh7WmnNS+D8ivQmvHdJxlbo1BLYez6w/iE0HW/HzN77G/63ejlMffBcbD7QkfDxHFEV0e+JX/8WisctDPZUiICFEHD385tzbsBggix5LyIVbp0tuUqJzxuRzZAKAnGxtC0r/MfYhR4iTBz5nLAEhJIkla6ALZgSkPkLs5m3zSUJIElu6SHlCRUIHgHBHKDfYLO1TpJ9ELnWnNvcckTsKd7n96NRMlD41+niphB+uFqSB3Wi4iOGOitkoMGeA9wxSN3bMHsXyl0J+oOazsLlrHC6ERuWmsvJ73s/ozNuZGIIA1G+Tw4BcTKkdIfUAWXl6fct+4O2fAp8+zj6X4LPbIm8cbU4fRBEwCPFdJS4utDoTK80U9R0hPuhWK8cIUJpD6q1DEIS4JfSxEqWBxHKE4jW6jCfoOLIjlJcWtt1hM8tian+cxor8dY7KTYVBlYTOq9eauj1xR27Utrnw7o56+f295ltlEFT57PH6MnE63X5NNxAAmroSO4coitgple5XHmrH3z+rYQNtdzUldLyalz47hIm/fA+vflHT62MB5u5+56mNWPDnDdRtWwUJIeLoORohpD6Gd2DurSPEK8+kpowWf1fv1iOJlFyVECo2SedIRAilZAEGdtPNQWdYZ2mzt3dCKF9g+6u7SmcHpTwjvbCY6jmh6zDL4wH7I86djMwUiyJgtIRQSqYy983Lwmv8NXBHJc1qguBuV/oi5YxRjhcEpSS/aacciuCJyv5gSM7XGJmbysrt3W2s2q14Kqu8K/0WO373fwAAqZKYUucI8cqlsCGn/76DDZFd9wDw13mAh33v0nVCY1wQ5KRZdau9OImExrQGrnKKZCHk0WzsGKuHEKdQGiKr58a0y522++4IxRt9whPCtZpkquEJ8ZGhMQAYLVUCxhu3wUVKpFuXm2aF2SggJELu36TH4pe/xM1/Z+XuFpMBZ43NxYZ75uG3l00E0IvZbZKjlJtmjXpNeuI2ksPtbs1QHBfSiRIMiXj6Y9Z9ful/d8v/5PSGug63vJ4t1drjZ+Kxv6kbc/7wEV7poxgbipAQIo6eoxFCBqOSV9QmdZ7ttRBi100XXABEGH3dYdvjIuXn5BoUITTGzsd0JCCEDAY5ByhX6FTNGhNh9kiluvFCY9J7mCuy/Z3egCwiMryJCCEpab2zDulWJkJ6vAFl4KrdDHRIf7iyR2mfQ2qwmO5h+UiuSCFkMwGtkuPiKAGs4f/1y8Kodb/s5nAxdbjdjUBIhM1sYDf2Q5+yfYefrrhc4y5kj1J4LNUa7QjJlVrchan5HKjewHKM0ovZINkdb7DXwYVQxE2I3wTj5QcBkOeqtUfcdLyBoJxnFMsRykuzwiCw6jd100NOvNAYAIyVGjvuqtfOjVFCY9rn4OuraXNpzm5j62Dn0KvEm1iSCYCVkOshiqKuIwQAY/ISE0J6iewGg4ACSRTWx3Gm1MnQp0htEEqy7HLyeqJC6BDvn5Vjx81zRmP6iCy5tUNjgo7QN0fYWk4qSMNVZ5SiXJqDVx1HVEayfm8zjkjiqdPtx/KP98c5Ihq1G1ep0aw0ET7Y1YRDra4+u1JDERJCxNHDK8d620OIw0NYfDhponPGOJKQSoMbNvhgCEmWb6LJ0qnMrcmGUlI7iguheD2E5HMoQsjpY/k9qfBA4GtJ0BHKCjEh1ONVEq4dXqkXUSwhxL8HATcKzGztLp+SI5Rv8SvT5tOLtM4gh8fsTiaE3BGhsTSrGWiRGjvmjok+PlvKM2o9EJUj1KbKyTEYBKD6f2zfsrOU40/6Nnus/RwIBlQ5QqrQmFyyLrkw/1vGHid/H5ghDZD98kVpvdLxviCCKjdGt2IsFAI6agFVEqxeaIzn1FhNhpgixmQ0yNdp0HAA4oXGAGBcEfv51iv5VkJj2o7Q6Lw05KRa4PIF8aXOzU8ZZKt9jhmj2O/klzXtupVjDV0euHysT1SpRp8oLugSFUL5GvlbPB/rSAw3JRgS5TDYqcUO3DxntPxcIlV4anhDxrJsO743fTheu3mmnO+UqCO08wj7vk0uycRvL5uI339nknzu3iRcvywJj8nS9VdurE6oW7ia/aou6V8e6psQ4i0xdjd0xw1RHiuQECKOntIz2GPx1Nj76cGdG4+qUqs3WNl/fOmCS55iD8EIWKL/K9W+fiY7Hqp5ZWbenTpBUSaF0PKETnS5A/AFQ8iWkq9hSgEscYaVSuIkI8DCTj1ev1I15pYGTmr1EOKYrLJ7VWpUwmud0k1ymFSNBkuavmCVHCFbTy0AwBcMSX2I2B/bdJsJaNnH9s09Kfp47gi1HUQKD2tJOUJczKRapD5E3BEacaZyfPYoFioL+oD2KlWOUHRorCgzBehuBPa+C0AAzrqTiSGDmY0jadguV42xcyhiik+eD+shtOV54E+nAo9NAD78lfKW6CRLq8NighA7vFYo3by1QiG8fJ6PwtBifBH7+d6lMyKCO0J6+T0Gg4BZY5kjuX5vs+Y+8Xozjc5LQ26aBd5ACF/VapeeczeoNNsuNwZVMyaB0JjHH5RbPuSlRTttvEIvliPU6fbLWvZfi8+UGzsCivhtdfrgT+AmzkNj6gRynrOVuCPEvm+nFrPfu+HZKRAE9g9GS09ieTpHOtz4cDfLKXr0e5NxUkEafIEQ3vumIaHjOer3/qvDHfAFQmH/JCRCrTTzzRsIYWttB+5+7St83su+Sv/dXo+7X/tqyFTPkRAijp7RZwP3VAFzft634yOdm94mS9sURyjPpOohFOcGpRzPhJgdHhikpoxZRt6dOkGXi88bgzKTKAdSKCNeWAyQHaF0PxNC3R6laszmkpoyxhJCgBweKzEqZfzcESqUkrDl5o1aSI6QubtW3uT2B5UeQtZ4Qkj6z7urDukG9gfeLQkQl7r8vu0g4GwGjFalbB9gIcY8qb9R0y5VjlB0svSwTJvcvBH5p7Brp+YC4xawbTvehNVkhEXqzaMOjyml89JNv6MW+PedQLf0Pm98Euhi4pOHmzoi/vNWegjph8U4RQ79G6cyhDZWh2uHdLxXdtbUyHlgMUaOzJaqptbv0xZC8bp1C4KA00ey30u9m97BZv38IEAJjdV1uDWbXALK98ZiNMCREl0VWCiP+9AXIfz9SLeZYIrozZRtt8Ak5YW1xhAhoijiH58fwprt7OdghKonUoEjsUG4HC6ExhezvzNWk1F2tuLlXHFe3HQIwZCIb43Kxpj8dFw0iYXS3/m6PqHjOftUc/c8/hBO+r//4pxHP4bHn7ggUVecXfH0JrxWeRhXP/t5r9bxh/f24LXKw3j7qzo8/ckB/N/q7djTxxEo/cEJIYSqqqowb948jB8/HhMnToTT2bvYLJEA9uzEhUckkbk8R+EIFdt82udM4HgASJUqplJ4Cb410TwjaeaZavirnHOUSKhPcoRSvK0wIIQut+IImXt4D6FhekczpPBYIVRCiI/GEFvDrqOJ1J3a0HlIrjxz+4LhXaV5aCxHIzRmz5a/d9new/IaABbqAwC7VZVnlHcSc7LU5I1jj827ldCY6r9GniNUlJGiDKIdpqpeGzWXPdZvU9aM8IRpOUGZC6GvXwUgAqUzgOFnsFYF//sTe0t0kqWVirGIRGlRBALhIRO9m7coivJNJdYQ2jSrSQ41aQ0O7ZAdIf2xNGdJjtCOui60alSPxasaA4BvjWLhXb0p7tWtqmR4DbJSLXJS+EGNsSdAeNhSy2krlt01fUeoPcZrMRgEOTcslpD5sqYD9721A93eACaXZOD8U5V/IApkYRs/NNba40WDJIBPKVL+zozIZd/PatUsND3cvqCcmHz9mWxUz4WTmRD6dH+LpjjWQhRF7JMcIXXYsbrVpfv90DrHYY2eWIFeuEqBYEh22j7Z24yXNh3C3z+riRsyHUhOCCF03XXX4eGHH8bOnTvxySefwGqNnyRJDCJRjlBvhRD7j9kBN8Y4AtrnjIXJymZtgblKACB4eOVZgo6QtAbe3RoACs3SHxd7Ao5Qah4gGGBACDnoQpvLD48/BCOCMHqkWH5aDDcHkHOICkXmKrl8AdnJyA5KlV56fYgA2RES2mvCcnx4snSGFUoel5YjBMh5Qplu9oebW99hXbLbpHNojR1RCSE7T5aWRIzbF5TDQMWZKYojpC7jL5QG0NZ/DYiiqoReEUJKXo6VCZdtr7Anpi4E5t3HPq98HnC16fYR0hyv0VkHPDOXhdh4hR4UIRSZI9Tq9MHpC7KCuxhCCFBuors0/mtOxBHKT7fhFCnXaKPGFPdY8844Z4xkQqjyULtmbgsXMYUxXDK5cqxZ+7//5jgz4IoScoRiJ34nUkXHv7+TSjLw5q1nhnXs5o5QIqGxrTUdAFhfJV7NCSgOUyKO0L+21aHT7cfw7BScewrrQj8yNxUThjkQDIlYm2B4rLnbi25PAAaB9WhSk2iYr7nHC69Own2iHOnwyMJp7TeNqOtwI81qwjmnJFCYMkAc90Lom2++gdlsxqxZswAA2dnZMJliN2IjBpnIfj+9rhpjf+BnlJhx67ck0dHbCjaVqzQ6L1VJLE70PNIa8izKH9dCo/QfTiKhMVX1nENwolH6Q+9Q5S3FFXeSEMoNsfCH0xuUc4TSpdyj2KEx6Y+jrxsFJnYjUAuhErGJNUA0p+oLKik8lu46xNYgCSDuyNgtJkVMaQ2izZcmpzftVkJj0jmOSC5AmtUEh9WoCCF1eK1gPMsPc7UA3Q2avYS4I5KTagEObwbaDrDcpPEXs1EhOWNYntLhLbIQ6nT7w8rfoxyhriPA385jTpSzGXhjkewMFelMb+eJuEUOG6ym2I1Ix0nhMa08oUQcIYBVLQHaN714VWMAZFfK7Q+GNcrkJDK2ZGycPKF4M+CUfkbxQ2N67wd3Q5pijS2R/oEoyogeqMtzhJq6vHGTnT+vYqLzjJHhrjAXQok4QtyB+85pJWFrOftkJhwSnd/G3aCynFTcNHsUll4+UZ4Bl2gpP++srWX+JxpeU1fL8fyk808tlBtnJoOkC6H169fjoosuQnFxMQRBwOrVq6P2WbFiBUaOHAmbzYby8nJs2LAh4fPv27cPaWlpuPjii3Haaafht7/9bT+unugX1Dd4qwMw9lKo8qox0YV0kecIZervH+McU/MNeGZhudyLJuGRIdJ+BWbFOcjhQigRRwiQRVcGnLJtn21wKefXa6bIkfKU0kMsmdXlC8iJp6le3pQxhiNktsmu02gT29/tD6DHK7kwISkfIWeUfhhUCpnZu6ulNUiOkBQaYzlCkhDK1hBCPEeodR9STeyPJO8srYTFbBDaq1hyvdEa3hfJnKK4VQ1fa/YSCguN7X2XbRx3ARPDgsDCYwBw+AvWfwnMOOpSiSl5LdwR2vI8m7+WPZp9vxt3sM7ZUPoARYZSeFgsspMzNG6u3BHaGzFewhcIya8tVlgL0O7UDbAme/H6CAGAzWyA2ci+710a3Y2bExhbwsWU3o03npjiorKlR3/uWdwGk1wIxQhtdar7b0WQLzlCvmAobtUWn3R/eoQQKpO+59Ut8R0h3j8qsiUBd9cOJnCOm1+qxDV/Y3k8o/PSYLeY8IPTSzFxGPub0xAj1KjmsJQoPbkkU8614iReiRe93kunxvi7NAgkXQg5nU5MnjwZTz75pObzq1atwh133IH77rsPW7duxaxZs1BRUYGaGqWHQXl5OSZMmBD1ceTIEfj9fmzYsAHLly/Hpk2bsG7dOqxbt26wXh6RCGrR0puQFoeLFW8XGyoK9NkR+v2FozA6QwBE6Y9soqExab8sk/IHRW7QmGg7AOm1OwSnPJG8SCqFT+h9kZy0VKmzdo83IDshVrfUiyhyNlwkUh5SsTTqw6XKEcoWpRBdrDwjqUeRrasaQHQJfqrVqDTO1HKEMkrlyrHCYL20BnZsvbqBYZ2UH1Q0KVogFrKmeaj/OqqXkCiKqgRlC9AgjR3h4gdQGkMe3gyLySC7SupeQtxVkZOl9/6XPc7+GfDt37HPt7/OlqNyhNQOAneEwkrN678Clp0CvHq1UkUJ5eYdeePtcLPXIgiAIyW2UE6Th9CGn6PbE5D/M9cb08GuIcjDdPlQXzU80TmWI8RnquklKscTQtmplrhzz+KNHJFL6Hvih9cyNZpUWk1G2W2KlSfU4w1gh5QozcOKHJ5HdbC5J27VFhfdwyJm2vFzVMURQm5fEO9+0wBRBIwGAeefWiA/l0ioUQ3PDxqdl4aKiUUodNjk5qmJJo9zF4yHYXPTrJgxKk57kQEm6TGiiooKVFRU6D6/bNky3HjjjVi0aBEA4LHHHsPatWvx1FNPYenSpQCAyspK3eNLSkowffp0DB/OKm4WLFiAbdu24bzzztPc3+v1wutVfri7urRLVol+RB0K623FGKCIFW+3qgQ/s3fnkMVUpxIWE4zsptyL4x1QhJDclyiR0BggizeHKs8o3+wCAkgsXCgJLj7otcPlhz/I/sianZIQiuUIAarmkspwSDlHKMjnpsWI5Usuj6mLVZ4FQiJ8gZCc52M3GxQhpOUI8cqxI1uR76kCkCc7GHXqvJwjW9n+6vwgTtEkYPs/gYavkGY9F4AixHq8rLUBAOSkWlVjRyYox5dMZ4+HK4FQEBkpZqk5pQ9AKkRRlAVJlt3CRoI0bAcgAGPns3wzwchCgB01KJCS2D1+5iDwXB6lWZ+UXOz3AG/eBHTXA7v/DfztAHDjWsCWgRQzDxOGuyByw8wUc9wu2TzU2BPhCPHS+TSrKW6IzpFiRqvTF+UIeQNKYr5eWAtQEtRbNZpLAvGFkCAIGJaVgoPNTtS2u6LdNPSPI8Rfi5YjBLCE6XaXH41dHpxcmK65T+WhdgRDIkqyUqK6j4/KS0Oa1YQebwC7G7pwarH2P27BkCjnlkWeg8/sY7k//rAcJDX8vbYYDdj24Hny+Bv+OgDICd3xUCf333neSRBFNrLjy5qOmO+nGu4I3TJ3NKpbnZh3cn5Udd9gk3RHKBY+nw+VlZWYP39+2Pb58+dj48aNCZ1j+vTpaGxsRHt7O0KhENavX49TTjlFd/+lS5ciIyND/uACihhAwhyhXuYHAUrVl6cLcHdI5+ybI8TElCpROuESfCaEUkTlv7MMUTpPwqGxTABAlqCcI8/IHaEEBKL03lmlESP8j5tBCEHokRIqY+UIAbJoy5HcLLdfEUJpAalaKNbYEUkkCa4WAEyEuX1KTkk+2llVlsEkz0eLQgptZXmY68vzjPh/nAUOmyKmcsdGH1/IGtah/mu5aoy/Bh4Ws1uMSAl0Al2sug0FqkG0+eNZHpSvG2jeI4+u4DfHHm9AduwcNjOwdy37Yvjp7P2zOZRKtqoNsJmNcshJ/Z83r5wZzh2hj5ey0SOpeUBaAdC8C/j6n+xtlRLHXRFl59yFyY4TFlOfI7J0vS3OiA41fNBtV4QzxR0es1FARgxnKldqE9DSreMIxckRAlhzQ0C/M7NSAae9Dp7XFauzc6ebJ6Brn6Mki61he512TyVAaTMQ6QYBzJnhHaY361ThAdJctZAIk0GZOcdx2MxyGLK6RT/XSN29XC2CAOW90Gr2qSYQDOGdr45gm5SPxH9m1bPwYuVcqeEO1rhCB3596UScc0pBnCMGniEthFpaWhAMBlFQEP5GFRQUoKEhsUx5k8mE3/72t5g9ezYmTZqEsWPH4sILL9Td/95770VnZ6f8UVtbq7sv0U+o3Zs+CSHJzQn5gR7J+ehtjpDaVeptorRqDWa/ksPBc3XidpWW18Cul2tW/ijlGKU/1om8L3yAbdAJEwLyH7fhFjeEkB+AEF8ISa5SttQDSR0aS/UnIoRYvxoh5EeuJOJc/oB8jjy/qieSXi6YlIid6mM3En7j5mLGYTMDnVJonCd4q+GhsY5DyDKFJ2u3qMNiTTultZSGf68NRkXIHN4sOws8j4bnXVmMBjYZnucZ8c7YAEu6BoCq9cr1EF59xmevlWXbWV7Q1r+zJy54VOmSvWcNAMg3MJc/GJa0XSeFKoapq84+ewpYNh5YtRA4sk3ezMOEkUJI6SodX0zx8FukI8SdnJxUa9ig1EjUjlCsyrNY4TXuoNXoJBrHGzkytTQTALC3sUe39FxxhLSFEK9w4n2GtNgtVfjx60XC84a+qNYXQjwsVqiRtA0oPZv4jDctYnUv16tojOS/Oxrw01e2yq9JXeWY34u+SsGQKCdcl2m4ecliSAshTmQ/CVEU43ZzVVNRUYHt27djx44dWLZsWcx9rVYrHA4HXnrpJXzrW9/COeec06c1E71ALVp6O14DkDpISz8PHaxaqfehMbUjJAmYRBOlAflGKvhdMEJKDA50sOcSDY1Ja85T5RllG3ohhGwZ4O+DOuF6hIWH6PLiJ1xLoi1TCuuxqjF2U0jxSpVnsYSQ2Sb3XhpmVsQUz/PJ8Uk9kbTCYhy5p1KTfLz6MdVqZE0QAe0mk/ZswMK+n4WC1Kk7whHKSbWowmIaQ2h5eKyuUnY4+A2WuyGOFBN7tw9tYvuOOVc5Xi2ERBGpEbPX3L6g/B90WY4daN7DKt1MKUxQnXyBdPwGwNMpuzmiCHhUScI8Z0O+MTlbgA8eZonbu94GXrpMrl5LlXOEIhyhOF2l1ejlCCWSHwQoN2N/UJQFJUcURUUIJZBwfUhHCHXESfzOSbPK1Wtf6Lgx8VoSnH9qIYwGAd8c6dJNeObViTz8FMn0EZIQqtJuRwDoDBlWkUieUHtk3ywVXAh1q/IJtahRNVGcVpaFKcMz5a/zEwg1cuo73fAFQzAbhZjDigebIS2EcnNzYTQao9yfpqamKJeov1m8eDF27tyJzZs3D+h1CBy9I2QwKEKmlU1njpsLE4k6vMaFUB8cIYD1IjIjgFSe69NLR0juag0gizdoTEQgGozyOTKFHjk/qNQsvZ54bhAgh/Eckpvl9gVkEWHxcCEU53dPms9WbGJ/nF3eoJyXkiENdNVMlOZI67S6mRDiDgZ/dBi8SlK83vw1Kek7LyR16pYShJWRFlZW2QUAhROij+fVa+1VsiPEK4l4fpDDZga6G1gITTAqPZAAlnxttLJu1W0H5T40XMzxG0u6zcSEVrVUCTv8dNbXKncMCxGG/MC+dbCZjHKUVl31xccd8FANNj0J+F0s5ym9GHC3ya6SnhDSdYQ8nUwsqm7SvNtzZNJ2Ik4OANjMRjlpO7KxY5vTJ/eoKciI5QhJQqgtniOkL/rPkGan8fL2SPjr0ztHdqpFTvD9j44rFG+O3KSSDFiMBrT0eHXL6GXHT08I5cUXQrEG+6ZZTbJTGKuXEBeXN80ehddvmRlW6p7PqyITCI3tk2adDc+2x81pG0yGtBCyWCwoLy+PqvJat24dZs6cmaRVEf2OWnD0RQgBihDh1V6xBpRqHq9yhHhorDeOkMkCmNgfhJklZmTI/X+ExMN0PEfIoPxRzIAkhBJ9X6T9MqFY5UUmPmw1ESHE/rinh9gxPd6g3NnZ6JZK8GMlS6ueLzRyVykgi5g0t5STo9VMkSOV8JtdLMzp9AUhiqKcK5QdYAIJtgz9qj4pQTknyNbMw2otWhVjWo6Q1FwS7YdU88YiHSEz0LKH7Zc9kv0McMwpSt5R0y4ltCW9Brl0PtvO3G0+hHbELOUcJ0vjQvasgcEgwG42hp0DiHCE3B3AF39lT8z7f8CUH7DPpZBbmsYQW0DDEQr6gVd+APx+JPDUTGDzs/K+iiOkI4RiODkcJTwWHpbi+VO5aZaYSdtcCNW0OqOcFFEU4yZLA0rejpYjxM4RX0xdMIk5l+/tbNR8Psx91MBmNsrOymad8Jhm404ViThCrXG6hhfJ4bH4yeNa+V+KIxQ/NLZyYzUA7bypZJJ0IdTT04Nt27Zh27ZtANg4jG3btsnl8UuWLMGzzz6L5557Drt27cKdd96Jmpoa3HzzzQO6ruXLl2P8+PGYPn36gF6HAAvX8AGpfakaA8LGZEAwxC7x1jxeXYLfy67SEef43YIy3HdusbIuQ4K/ZpIQUjdRTBelnKNE3xcuhFSjPnK4w5RQ5Rn7A5UqVZ7xG5wVPhi4QIwVGgNkR6iACyG/EhqzeiQxFavDtSTYjM4mACKCIRHeQEjuRZTpk248GRr5QRzJEcrwMdHEXRD55mQ3Kd2f8zWEkDSAFl11yLSx7x93CXhIhwmhGLPXpC7baDsQ5Qjxhn05aVJ3a1kInaUcz0Nth5krLQ+yVQ+hlYTQ8Gw7c5V8Pey6Jy8AplzNdjrwIdBZpyqfD68ai5oztnctc5H4PxWfPAL43cprhkaOEO8hlB4/vMaFQaQjpDu2JIKSLDsEgQnkSDHl9AVlJzS2EGK/Tzvru6LcrR5vQO58rFc1BkDuwXNEYwCsR9V0MlbIcYzU5LJep0GkIoTihMaao0Uhp60ntiArTGBsSUcMh4wnS8frI7S1ph2f7G2G0SDgljmjY+472CRdCG3ZsgVTp07F1KlscvmSJUswdepUPPDAAwCAK6+8Eo899hgefvhhTJkyBevXr8eaNWtQVlY2oOui0Nggw12TvjpCatGSXtTnpoxMCPUhR0i1hgyDG5edkt77c0jOWJqq8iwtxIVQbx0h5RyZPM8okVCflM+UKuU38TyjIqO0DqMl/nn4AFqBh9eUhGuLR/rPN1a4UBJCQtArO2suX1B2hNI8POE6husnVaSlS6KJh/f4zXeY1QUEpD/8Wu5UWiF7raGAPKeNizklNKaavaZVvcaH0LYekHOE5NlrkpBIt5rC84PCumRLAq2jBvA5lcoxaR3+YEi+eZVkqUaOjDiTVTvmjGbz08QQsGeNriMkj7XgOSTb/sEez7iZiU1nE/Dli8prRowcIbUj1N0IvHIV64nkVxUASPtETl5PdJCtzWyUB9lG5glxUWc1GcLGYkSS77ChLMcOUQS+iaj84u5HvHNwd0wrt4YLbrNRkN8zLXgyNu8HFUldR+zQGM+X6vYGoubhceQQnUaOEKAMBY6VMC2HCmM0mGx1+uAP6o/fePoTlrZw+dRhmm0PkknShdDcuXMhimLUx8qVK+V9br31VlRXV8Pr9aKyshKzZ89O3oKJgaF4CmAwh5cx9wa1I+SIM5w01vF9rRpT7+/5/+29eZxV1Zk1vM6d51vDrZGCAkNQEEQBjTgRNYAkKL7mtW3T2tgq74etGOLwdkx+Ch3TcehoZzCaLyafsdP5XmJaMUbbgVZxQo0iBpwQGQJCQUFRt4ZbVXfc7x97P/ucc8dzbiFVsfb6/epX1K1z9z1316HOqvWsZz1VqkrCKxVkupoTzIp17BIhCnMEENVspG0LL5Inm4AHaXmTnOAV64WaKkcKCMWoHvwG0z+UwVCa/4J0EhEqZyB3eaUCRv6mRDIjvTHBQeHJqCkTbSEUocAg9xcSEaObQqsjLs6jwVzSIjgc0ohdl6JgR/76VBaK+t0GInRs4RpSEdohS2P5AZMhr0sOiMW4WeZzCdTJLjwc+liuQfvQER9CjvEbdkPIC3z6tlhnjr5Gu7AQ7N8iPUKD6SwyhhsWmbYbwz6gv1OPA5hzJXDmt/i/X/sJwFjFrrEGavHu2g48eDaw9SmeibReT/QnwpUfqrivgjHYCLqR7j5sLgl1WyiLEYi05b+XSv4gAnlrhtK5AgJgHGJbrrFHjnApQWIqESGfWw93LJUFJH1xJRShJkE895UjQmVKY3UBj0yZPlRkoC+BOs7+52ybtoWjgBEnQqMVqjR2lHHxw8CNW4u3Q1uBUXmx6w8CSucIVXMOyd7qfEaCSPmzOhHyZyiLyGJpTBxXoxnLazYUIV8NN/4CqEWfvEmOc4tzqlQWA+TNu47FAehlEw05aIOkCFXopBOq0Hg3eZUyUhHyDRha8EtBkGHvACcx+V1jjULlKVtCFeWx2hQnU+SVohtnxO8GDhIRKlIaMyhCpCzQe+gjIuQzjBypL1IukENotxpUJf7cT6VR2g+N5fRWeaOqRFEC+7dIRcn4XgCdCDWEvTy3iGU5mWo4Fpj5DW767v0UOLxD9wiVJEJCEXr5h7xzjf4oee0nwJ4/8bcZJAUhvzRWRBFKDQBv/396bpRAex0NLc1ThCx4ewiBEiNHusuM1zAiZFB6+vM64CoZpQn0GvlDfQGuNJG3raUMOWySI1xKESHDkOEimBTjpPKTzuKDcAFdsSq2rw6HJn1C5bxKpFjVW/CRHW0oIlQCqjR2lOF0AcFhGOiMpCVahSJULEeoytJY1YqQICqeTD808MnzXiJFthUhnUyFSGGyEimgabJsVaf1yV+irU5qwbdOhKK5OABDxoxjABr5Tip10gkiNE6M+ujqT8nmJXe/6DwrR3jF91z9+wAw9KcyyOUM4zVyggiV8yoJUh5J8tcbzCuN1btSvCsMKF4aE+NG0LcPEUdarEGlMYMiRENo6XgjqHvt4EcFN25Tx9ihj3n3mjuoD64F9HDJzg/g1ZicFUblsWyOyXJhY8QL7FjPj5/xP/lnt48ndQPA3o26IjRYpryWy+rZSv/j/wVOuAQAA/70C75vJRQh8snImz5jwB/+EXjyW8AvzuZJ3wJSESpRGrOiCIW8ZnJKkMbgCmTKTTlS0M34BL07sfx50GvEi8wso/EdYa9LljWLoRIRqkTKprbw31EfdfSV9BmVM0sDwKlf4P+fn32veL5fNsckeS4XuDlSUERI4fMBU2lsGIqQqX2+WkWop0pFqAYA4EAOIQwZOs9gvfNMEKGYofPMnxNEyGqpTxIhfbxMm4dKYxaIkDgmIkZykCLU6hHn5I0WL0cZIZSaVocwbYu5UJoGOCgNupx6KJQIRzqBCAbAGL/hEbGLpGkIbRlFSHSOBQc4ESICQiSgLScIWaipOMkM1MmfR1OWE6aBvNJY2OcCDu/gxxeLFCimCIkp36aOMSqLtZ7EYxQItZM4OcoMAV2fFPiEuvqTyDHAodHIEdFJZxxdIkeOvIWoaJ83KkKDKb2zMBb2Anve5G37vhruUSLT9q5XAcYMHiGzIrRPKEKtpAht+g/g/bXiRQ4D/76E5yRBV54O5ykpZLhuruAzAvSQyvw4gbgcn1L5hh0uoZDp8+zKqx96NEMhEbKazVRqqC/Ah/ISSStVGjsmFoLbqaEvmZHXlBFD6ayMNSiltC0WHXT/9d7+orPT+obS8g8ZRYQUFD4reA03+apKY4KwpBPAYLf5MasYrkfI7eNlCPDOManqeKPWzd/ixlvnMJTXxBBWy2RK+HcoXRowtODbKI2FMt0AGA5Rec0jiJ0V5U8oQo1aHICuOETcDFqf8AiVK415AnIv2pyH5Ro0Z8w/JFrwyylCojTmT3DCQyUpuuk1pUR4Z7GyGEH4hBpSn4o1OGFIGD1Ch60pQrLzLEmlMUPHGBmlx80yP9/h0HOSDD4huvnTzTMW8sI5eJjPOQPMXj0qtX36tql9ntQDIiNup8bN3x89xY+fspBft+NP4cbzvg6gazti1DVm6PjK5Zg+yLbGz1WldbxhBvP+iXuwUn16JhJ10OWVtfbG9XKhxGP/D/CDNuDRZUD3X+TDRArz1yg3eT4f4bwRLoRKrfOEmjKKkFTZKhAh8vgU8whRyc3pKD36xONyYHIj/0OQfDxGkBrkdGgllakzJjcg6nfjYF+yaCQBqahBjxMe1+ijHaPvjEYJlEforwxGRaia0pjx+b2i3GF3TMdwPUKAVBaiWkLPArKTkl2kNOahsR+WFSHuM6o1GK4bhPG5YpgiIMmSmyURwqBUhJqcFA5pIWlbZAk1gP9SpdTaNk+Cd0E5XLqRuBSEMvgFTxyArhY4NMCVEBJ+2dIYJ0LuPp5inW+WrkkLMlWOkAnfT32S1jCPC6lxDvGOMaB42jYpQt27EHUJEiXOw+SpoSiAlpmFa0if0GZDC715fltjxAsc2MKPq51k/v9AitD+LYi4+WvnmH4eVI6qIWPw1qf58ZSD5Pbra+x62TCBXlcwDiWSSGcZHBrQFPYCHX/mKpA3Cpz1v4HpX+cHbuUlt0CJstan+QGEu98ANq/hJGrLI8Dvl8pjyTNVEDBpw2cULtE5Vi7E0AjZNTaQMo1OAUp04hVBk+jaOlDE7EzKVG3AXXb0yVQxOPbDjsIh49If5HeXNH57XA451f6pLfsK16gw9mSkoYhQCSiP0F8ZTB6hKgblurxSjdG7xo6yRwjQJ9BrAzqZsTN2RHRbGctq7jR1ntVYWyNgHrwKADXC71ORfACAJ8jLMQDqtV6pCDVKImRdEarPCSIk1mhxG4bQVspnEoR4gosrfJTJEvK6dFXJQmnM2d8BNzJIZnLc6yCIUIhmr4XLkEOh8kQGeC5afmksRiNHAjEz+SAEGzi5ZTlZiiNFiNaI+N28xR4oXl4zGabNpTFTxxgFTOYnbddM4OeRS8N78D14xKRw2ge6ydUG3Lzr7PB2nuU12TCeiIIid74ifTPdA2nZvUb+oMawj08iFzPaMPF0riodu4h/veNFID2kK0KpPEUof/baC9/nn4/9Ku9K3bcJOMBnzOldeBU8QowB7/yGK0u/vwJI6GnU1BqfT6ZkiGEFIkSeqxzTDfQEq2ndsjRWZNaXVUImfUL7ixChCv4gwumT+e8NSo82rTE4ev1BgCJCCp8XkPLi8lkfaVGwRt6NyK6aY1KEqs0iqgEgSmN2U6UNx0ZEEKMXKTiy4i9vmx6hWkNpLGA3z0iEKsbQIwMIYw4yXFshQpyg1AoiRGWTBjG2wxI5FD6hNicnQpSTEva5ddWvnCIUjAEuPzQwNGt6lhC9H3/KwsgRQYRCg2S4NhOh2uRe03EF0DSgnhuxWzLCqyTLa/xz2JXTS1rFIgWaBBHq/MBQGuPPJaWtMWzwB9HxxnMQLfnavnfkmA0qEZrmch0USds17eb/T5MEEdr1Kmr9bpnAQOqLVLcoQZmIEM1sa57BFb70ALDzZQSoC89AHhhjBt9UANj9Jg+ZdHqARXfzUh0A/Pn/ADCUxvLIFE2el4brrU8DT1zHlaX31wKv6vMqh1sa87md8IvE8HyfUEG+Uwk0yRygQo/QYYtxAse1kCJUWBqjslYl8ziVTfNVOkAf0aGIkILCZwnKpalpr5xzUwr5ROiIKELVZRFFtYTeAl8FEQqwAbiQQYTmnWkOOYi0IoqYpd0ZMlxb3BPRXRbT9LA66TmyVBrjzw9nOBGiG1yD08aeCIISE8pWhyBT9d6sPqusnCKkaVIBaxCeq96hjCQxXjl7rYxvSnzPm+TvYyAt1Bxx4wwPiuGx5YbQCrJWKzrdqHuNziOaPgiA8T8Ciil2lJLd34kaDy+/JPJLY2FvaUUI0D1DXdsLBq/qIy0MI0fI20RoO5mXMxOdcPbtleNC6DyobNka9QOZFLBbDLIlJUnTdCKz7Vk9E8lAYroH0hgURvKWqI8nagPA1PM5QZwpRo5sfgTIZUvPXstXQDb+mn+msuPbDwED/OcZ9hYvjRWkdQNAx2bgoa8CT95gmt9WK31CZuO3VbM0EaGuRLIwz6jfWvcaKUK7uhKmES6ATtBqKpAYihPIjxIArGczjRQUESoB5RH6K0PbKcA5twKL7618bCkYiZDDLWeHWX/+EfAICSI0PpBGzEllIBtEyEC8okggQmTKG7E+6kOQynqDIuS0GzJJqpLBq1QLUoQsECFRxnOzFLxIyY6iekrJtjJyJM8v1SFKYxNcgpy5A5Xfj1Cemt38dY3pu64B0XlWThESxMQtSBMZc6kMEugvU9IiiDJhVJDCRF4LfiQp/E7RtuJ/BATquSoChhYHV8f680pjTSEncPAjfnzzjMI1qEMvvhth2UJPihCVxgyKUL6B3OU1zG/bqUcBiJuuaUL6vne48hOoBxoNpm0aP9KxWXbQJVIZadqmXKXGsJcPBd0rOukmzOWfv7iAXzf9+4E9b+pr5BEhUh8bw16gZy/wiZh1edEvuVqWTsiZbqESilBXviK09Rngl18B/vIa8PavgM2/k8dGhVqTnwx90CIRqg964HZqYKxwzIXVPKNYyIuw1wXGdFJK0DOEyq8R9JjVRiN6bPiuRgKKCJWA8gj9lcHhAM66yTyryS5o3hkANB5nX1k6Eh4hQQCunFWLv5spnmvHtO10mSbQR6WqZGMNcSyRqKCbQUsTobJIhMT7DpMiBSDKBAGxogh5wlzFAmT7O2AgVlbIIb0PUWKkX/BtwjOEcEvln7EgdGT0pptkwOOElhBm6bJEiCtCzqFuOMHnT6UyOaREO7KvTxChUqUxw/qhtCBTqQyyOSbVD2rvL+mN0zSpKjWDq0r5RKhdO8An3XvCxdcxECE5ZqNcaSxfEQL0wMjDOyQJGZQhlYY4gT1v8uPaTzeTd9lBt1WWxnIMsrXb5A/K5fRIgTaRtO3y6KSo48+6X8qgKhm715qjPuDd/5+b89tPBxqmAGes5AeKMSRhuRc6EUpnc1IBkQRk/Q+AbFIng898G+jnRJqUlvxQxUN9/OuY0Sz9598BPz6Rn5cADzQs3jlGxJ18ROVApK6kZ6qCIhSWfqnCDjjdI6TM0goKoxvUNg8AFz5g//lHUBEKsX7d8GyHxBjWiCKBZveQ6TFLEOccgvBb+A2/GK0SO7FGWNOJUEgMcrWkCDkcco2IISVbzlALWCFC/BgKlJQZMzReo5w/iCDOtVGU5GiNBm9O/xmXK40F6gDNAQ0MdejDYDprKqM4+4VXqdy4EFG+CwpPUiKZNfkwfESEyuYq8Q46StSmm91BuumL7rySqhLdwI1EqJhZuuzIEUH2urbrRuX8tG6fm4/nAMzBkABQP5mT42QP/MlD8mF6LzSOoq02wA3bQ3Gu6jYZSn0ySuC9AuM4wD016SyDpgkD+Ud/5N+gLKQvzgegAfG/AIlDRbvGiNBomiCH8d28C05zAFc+y1WlwW5g02/41tKYDUMLfS7HCktje9/hXqXuncDj13ADtwDN+8qfAL9fkrrKY0t035VZ0Sk3cNUI2s+hdM40wgUwdo0pRUhBYXRjzpVcFbr418XLA5VAyk2qHxiMi8fsKkLiBj8Y130sVfqMwtogGiQRqrH+fFEiJBLT4hO/oF1+wGnxFxkNoNX0gLaQGORq2cxOHXQGVSlM5TUbpTGa10Y3mkYmbvpWiJA415iDlwn3ixJdu08oUy5febLrcMo16oXnioaN+t1OaNSBVK4bT3Sl+ZJcQRhIZeTN2+3U4OojIlRu5Ah/r/VZvkZ/kpeUqPxSnxPEIlLCM0XZXOkEmkXnXo/wCNGNv9Gd1E3bDUWylYgIHd6hD5AV74P8RmGfWw+YpFltBJdXDsh1dn0sTcZkdja1zn8qlPzWk8zXLJGiA1sKhtgCuoISC3nh0bJ6LAHNbPNF9bLf3o1FzdLGOWNOh6bnKk2Yy/f3pMv417teBWDIEjKUxnoG08iIdvr6oBfIpnnHWjYFhMV1+9QN8o+35hKDUwsUIcaAx68FfrUQePf/8LwmgVJDecvNGTPCNMIlP5tpUJmlFRT+OvCl/wV8ezdw/P+o7vmBOu47AXiZAagiR0jc4Ae6DGSqxt4aXiIQCTS4xC9GO6qSSRFiaPHa7DozrDGriX45MvjTQnGzoggZXi9iUJUidrrXyDie7TM9XEeZSFaiAIRHqE4zq0rj3Yak7UrlNfE6ZByn0kutl+ndheXIochU8g7yUtxAKiv9QUGvS2+dj5ZRhESUQE2GVKUMugfSSGf5zTZMSdulyKHLK5WpceDHSn+PuFG2Zvbo51vsWpFEaGeBItQn57dVCJgkpenQx7qCIWevGZO2BREyzl0DdEWo8yMERTag0SxNP9+WqI+rW9kUv5ZJETOuuXdj0fb5w/15vpwPn+Sfj1vMP088nX/e8yaQTUvvjZEIEUGtCbh5AOGnb3EVyl8H/OPrnIxlU9IQLsds5HmE9ucnbX/4R+Dd/wD2vAE8vhxYf6c8Npjn2yKUmzNmhNfllNEK/SXGllQyXI8UFBEqAWWWHqMwjiewC00rHJxZLBumHOiGOHhYV4Rsl8aoLDWIOjJc2yIx/Jw9WhZepNHoISJkZ1wIP3ZqHb/RhjAIJxO/6K14hABTiY8gSY2V9nlBhDzZBFzQfzGHaT0rBFP8PGpAipDobpJJ2xYCJgURanXxNciX00YjRzRn+XMRZmlXMg4P0pwIiRtv0OMCegQBKasIcSIUSR0AwP9ip0DDqN8NVz/lKlWevdaU42vQOVDXWENSpDYXU4MAkyIUEDO69LRu0QHnyvIBr8bjjaC1D27VQxWTpAjx/RxnHDnSlvf7u2YiV32zSUQSuwCYSzmk+DVHfLzLC+AqktGrROndezcWLY1Jg3LAw7vLdm/g3zjua/xz4/H8553qBzr+rIcqGrrGKHtLhinueIl//sLZ/PcBddB9/BwAPY/IeB79yYw05TdHfUA2Azz/Pf7NlhP5540P8cdhzFXKK41JEiOIXX8n8MjfAw+czseeJPU/NEp1jlltwR8pKCJUAsosrVAV6g3DNz1h+8SKbvADh4ehCOlG5Zhr0P4anhAATawxiONJrLCjboljXal+PPets3Dj6WIRd4CPv7ACQQCbvLrc7yefkZXSWF4HHSFIs9esEExBhKKMEx/yoTQ57BOhFkmEhOrgEecRjJXv6PPXyrDPBsSRSGbkzSrq1YAeQRzKBYkKIhRMchLTl9TzkCJ+l7VcJUGEYhkiU+aOr9oBoeSUGjlSM4GTvsyg9GnR+6AbeCwtCJk3WpzskiJ08CNDqKKYmybOozHk0ctrxq4zgO9z0/EAgED3h/Jhmt9mmlW2XyRt09BZgkERCgsyVqw0Vhf0AB3vcrN1/WQ9xsDh4OZrANj1atHSGClC0ihNw3AnzeOfvyiI0CfrgFxWDpA1jguRuVk0tPW9R4Gubfz/zuVr+bWdOCjXlkNo80tj+STmtR8DH/yB507tWK/HC6B0Wrf0GRnN0rks8Mnz0jQ+klBESEHhSKJ+sv5vu/4gQFeEBoajCHECcM4kH05pcZoeswSHQ6pCy09twLkThdHSliJEc9d6MKUpjH+YJb62QmDy1mhy63K/T6ZkWyiNOZymDjqCDIe0oQiFc5yAyRulk0aOWB9C2yjIE4UYNtE8uEoKmaZJwtWoxZHM6F1J49y9QC7DM3rKZSIJguMf5K32iWRGko+w12LApCBCtYKsJJIZZLI5PWBSeJhKzvpzuuUa4xlfQypCdKNMkho0qXjJsaGwNDaQV16rQT9XW4ASAZO8POY++D7cTk2+F8BoLvYB+4UilO8XbJrO4wgGu2UgZlEiFDJ00eUbv6k89pfXZCdV3NA1ZkqVTvbpUQDHfJl/nnAqJ4sDXcDed4oOkC0oi33y3/zznCs5yTz+Iv71lkcAFC+NMcb0AbIBD5AelB1zstT3xgPcwwQgJHKVjOfBGCtsnz+8A3hoEfAfF/HPmcIwyKMJRYQUFI4kYgZFyK4/CND/Ck71cQ8AYF8REoTlSy0uQyijzTUEEbr65Ho4aVaZLUVIlASps6qaOAHxvqXhGwyeNHlqLBIqQZiMipA3Qz6jmsrPF0RIdrwJ1LM4/4clRYgTnQZhlj4oR46I87CUtE1EiPusSFUa76QogNbyg3mFIuQZOgQ3uNmablYhnz1FKDrESUzfUMbU6eQZpFyl5tJriHJXc5bIVBa5HJMlnFCiQpwA/f/qPyCvi4FUBulsDkNpXt6KJIWqFGric87yQT6hA+8XjNkg8tAS8eqKUD4RcnnkY9Hu9/npJDNy6ropVZqymWhmHIEUoT1vokYkdRsHr5oyhP6ygZPd2om6quR08zIZAGx/oWhKtonUAcCnfxKvLYzfJ1zCP3/4JJAaKNpF1zuYkTENzVEf8P7j3KAdHQ9c9CCPh+jdyxO3wdUnwFwaG0xn5bBjaZb+zyv1mISubcCGn2AkoYiQgsKRxHAVIW9U5ucA4KUEuz4jIixDPfwDGMYA2T7DGvbN0tI/UE2cgHi9ehf3foQxCAcTv2CthkwSEdKKECEr70cQoUC2Fw7oLcE1YvSHJUVIZAnVC7M0kRjqRLNk2hY+oRYHX0MqBnJsSQVVKRgToYpAk9aNvqGMVDHqvDl98GukzMBiQYSCQ5w0JVIZaZQO+1xwyFylMnsiCE5Tmis/A6mMCEXk3/b37TIdVwBfVJLPCZrudzLevP0DgtSVKhVSN1p8t6FTit/siQhNcHVzRdbhKiQxgCyB+wc+lQ8RsTSVxmSuUt4apGwN9SAmQkKNIzZIhakPefiYEEAvixHI/9T5gVTHzIqQwe/U3wl07wKg6blKbXP4tZlOAAfeLxqI2NHL16gNuHlIJZXBZi/lJe5T/hf/WowtCRYpr1HJz+3U+Hl2becz3zQn8JXV/KCXfyjOb2SgiJCCwpGEkQhVowg5HObykS9afbBjsncYPiNSdIxEyL5ZGkO9vGW3KkWIE5Va0YIfpfKWy1/8L/1ioHRp6KUxV8oGORTP18AQNawhB67a8AjVMiJCnMTI7jVLI0c4EWp15a1BBK+SQmYIVWxBF/qTGdn2Ps4pfjZOb3mCKTqn/IlPATAkklnDeA0P0CcSrsNlFCFRqopkRJ5RKitLax6nA654mY4xgli/CXEAnEwRqfO5HXD1knm8RBcdkb2evQh6+C0wIeIEZOhmWhi/Y1N4x1yJ9+Hu+1TvlEqSV0n8bAJuvf0+P2DS7ZfnQeXAXoPR2WRQJjLVepJ5DfJiHdpmUIQKO+Caoz5gj1CDGo7T/wDQNF3tOvBe0TgBGobbHPXzshh145Ga9MX5/PPedwDGECLzuIEI9RjCFDVN4/4iADhmHnD6Sj5GJZMEtr+IkYIiQgoKRxK+iH5ztJv/QzDe1OyWtIyvO9RbfRaRsbQl1ZwqFKFcGsgMVakI1QDQu7xqiYhYLYsBBWM2NOTgTNkojTndcu/qNL07xkeBflaIkBhAW5OLA9A9QmTAthQnIEpjpAjJFnyZtG1lCC337rRonMRREnOLw5CrVI50i/fqzAwigKQgU/wm1+AH73Q0HFcUghSGMrykN5jK5LXOU4ZQuaRtToRiokyYSGYlEQr73HqcQEkiJLxU6QRibk7k+pMZUxmoLtNZfg3jyBGZJcTfBylCzc5+/v9Pc5j/QCKIsSqBBCdd6SxDMiMCJo0hhqX2hMqEXZ8g4Cavk67mmBKyqSw2/hTzGsI4jgPvFy2NmeIEOj8AWJarpKS2NR3PCfRQHDi8o6jhuiBMkYjQ1Av49bb4R8Cy54E5/1C4R0cJigiVgGqfV6ga9EuvmtIYYM6UqYZMGROu+6lcYaH8YoTPWBqrQs0xdJ5VPXJEvPcg40Soxs54DYIkQnyNMAaggZnWrwiam2aYveZKCm+OFVImbv6RXByAPhojnI2b1i8LYYQmjxCVxqJykK2F8xAlq3EiA4k64ChtumLApCcou9fqtD6RRcRv+hO8QplyuMqTMrEXfqGoJZJZGaYY8Tr5bC+gfEq2IIV1Ob4XAymj8dtVmQh5gpJkj3fRGlnpqakJuOFOCHWr1J4QEYjvKQhVJCLUlNrFj6mdWFzBFIN2fb275UPkrZHZPV5NLxnlR3PUTOCBntkkalJ6uZJgIjF7hJJTQIR0v1SxIbRUXmuJ+ng6NgA0n6ATZqdbV5X2bZLKVL+JCBnCFLv/wjvpNIduto5NLsx7OspQRKgEVPu8QtUgydrODdsIU2msxv7z6Qbfs4erMYA15cIIoyJUjUfIMCKDq0pVGK7F6/lFdpBUhKogQqScyHBGl794yaMYBFEhRcjryEJLi3WsqGTi5u9iaVNKtj5yxAJJJRVElIOoNBbO2UnargEANIhIBQogjFlN2tY0U67SQCorb/ptFDAZbCwfBSDUL1+Kv6aRxLR6h7jiAJTfE3Et12T1IbR04w37jESovejTAciyFKlj/cmMbi6O+Azm8RKeKaMi5NUVoVyOSZWsLiFGhRTzGAGSCDm6d8qUbFJ0SEWJZQ9wo7TLV5jx5HDKP7rC/TvF83UCQopQU9jDPTkAH05tRJOIFzjwPoIy28ngETKSKcpVaplpXkPmKr1TlEwdklEAHj50ls7D7h9nnyEUEVJQONI49R95jP5Jl1f3/OGWxoxmaYDfrK16agrWGM7cNOMadC42jN+CCLmFubmGSlNVEKF6ESxZ7xBExM6+kiIkyNSEgCFwzorC5fbzTCnoYzYAIEA+IyulMfE+IqJMSAQkSETIRsBkTJjPqTRWmxHdXuXa7+VJm0nhPqEqtVKcQLgC4RYEx5M8DIAJjxC/6be4KegyyjuzSkEQoYggQoMGIhTyOg0Bk+Vmrwm/lKYTsrgcHuvWQx1LESGKCMgMotXDz5u66Kh7LNj7CT+m2ABaQC91de8sIBDUQVY7tEc/thjBFOWxQA8nXeksQyqTQzqbk6NcWhxxIDPI1bp8VSk2hT+e7EGtKAeaFCHjrDKKE8jPVWoVRGjfO7oiZOgak12SYZ/umcpfY4ShiJCCwpFGwxRgyc/kX3y2YSqN1dh/fr5yU+nmVAzeYZbGTGv0GNawoSqJY12pXgAMMUrJrsIjVCc6c5o9Vcxeo5u/KEOND4hf8u6A9dlrguzUk0EagDdF5TULREgO402YHg5kbOQqifdcJ8igLNFlBCGzRIQEsRQZSHsO87UaHZSrVOFaE+/VkUsjggQGhDcHAFqkqlRBKRBm6ZAcQquHQza6DX60UnlGgCRCDaIs2G/wGUV8FnKVXF6p0k1w8PPoHcrIUMew1wVnfBc/tpg/CJAeIRzeIctriVQGQ+ksUhkRBTBQKU6Aq8/e+Hb5EM+I0olIZIi66NoKA15dXrlGff82AOZQRqkIhV3AAR4VgOY8RYhM3B1/Rlj8dzCSqU5jJlKpOIERhiJCCgqjDcM2S+cRFrtlMcDcNXZEFKEq1hDvXctl8N35E3HeMUIlsBXKyNeoFR6hRvegaW1LED+PWqGCtAUo38k+qQuLDjg3MnBTPpMVRUg835dLQDO08fsyNnKVxHs2RgkAgD9rg0wJUtgsVKXdggjVi5JdxWvNrQ+pjWm9SKSysjTWSHEClYihICABQYSMs9fGC1KCYEP5BHOh9NSJIbRG8hD2unSvkoU4gTYHJ1N9Q2npmaoLeYA+3t5fsouO/lBKHERMjLHpT2ZkWczl0ODp2SWOLU+EHF3b4HWJDrhURpqtgx5n5S46YZgO9/DutITJIyTiBHKf8jK7J1TctO0JAekBNKd3y/dBMIVDdioipKCgYAXDVYRcPsBhUCqqIULGFvxq1BzA7BGqRlVyB7hsD2DZyXWYEhF/qdohMTJQkSsYVY0cyZse3+wVRMiW30kfewIAM2pEq3SlOWN5z9fAEMagfNibivN/WBo5wl8nkqcqeaWqZOE8xF40OoUiJMprdVk7cQKkjvViIKWrOfWaRWIoFE7euceQMPqMxEDYsmUxQCo9NISWn4fogPMmebYOoHeYFYNooW9moqQ0lNFTmIMeoF8YrksFTPqicj8nOQ7KNeKGSe0adYzll7QIsoV+K4Ie3WfUKzvx3JU9UyL1OtDLVSWe68TQO5SWhKYh8TE/tnlGYYnO4FWqE6ZtY2mMFKEWXwbo2W16zdECRYQUFEYb/MNUhDTNTDiGowgNdHF/AWC/NDZcRUjTdPI1GK8uHFIQoTDjN9k6J72XKtQcQWKqGkJL89+EIrTyNKG+BOrLm4sJLi83eMNg+AaDKxkX61hXhIKsz/SwJ0UktabyGuLGHRNEiEo4IVles0KEKGCyF+msYYQDs2geF9ezI5dGDfoxkNQ9QjFQ0naFMp8gOJG0ICBJXZki3xB8NbzDrBQE2WrMivltQxnp3WoMaPz/DlA+V0moKxNlOKSuCEUDbuCwKHnVlSBCRJAGu9Hk5dclV4QMZb5ukYlUighF+fugTrkc44NoSQ2K+Fzw9QlVqVSZT5DCaGq/PAcCKULjMoIEBRvtlbePAhQRUlAYbRiuIgSYb/RVeYQEEaISAWC/NFasa8w2mdJnllXVvSaIUCDXDw056Y+xRTDF+wgJJabBRT4j+7lKYQygKeLFma2i/dhKxxhBvB4ZpgNIwpET6pQNRSiQ7Tc97KSASSt7Il7HmKkEAIEUjdewrgjFRNL2fpFeHBE5SxUVIZdXH7+i9SCRykgFIkokr9KNNm8IbcLgVdLjBMqUxQDZQl+XISKUxmERpmg5TqB2IgCgFbpRmQII632aTmJKlcY8QXltjXP16u/FmM0UJyJUPlfJ0d8hH+pPZqQ/qLXGr5vHS/muopQ6ztegn0cux2TXWOOQCMtsHF1lMUARIQWF0YfheoQAM2kpN/up5PPFTZ7kfU+o0GhZCcWyiOyOCyEiWDUR4s93gCGMAX32mh2CKdUcTqLqiAhVURpbOqsOL9z4ZWgUPmjnL2MaQisUKZlr5PSUVy4IYi94KUxkKYHBkbSvCBmTugHAPUgBkxauNUH+yBMkp6RTrpKNkSONWjcGDOnUYRqoWymbSZTGvOle+DFkGkIby9Lw2ApESCgs0aQ+e43M0uNc4jxCTeUVP0Ec60AkJitHbUz09vE4AaenvMIl1mgVr8mzmQQR8rl1IlRbQhESa2u9+0xJ23JER9RKnAAnhf6EKI0Jha57IIVMjvG5wbKLbnSVxQBFhEpCBSoqjBiGG6gI5JXGLMzDykc+Yakq2LFIec224XqYipCppDQoU6ptEUzp76FRH1WUCsX7bvGleKv0MGavHRPmN5kaY6q0lTEsgug4WRo+8Bt2CIPQmA3vlSBuMhUbgENjcCSEMdjS7DVOdJqcZiIUyMTF962MHBFDaBEXA2TF9PMcqVsVCKY3IkI/gWatGz2DaUmmZJxApVwlobT6RSaSsTTW4oibzrMkiFiKIb6JlO4RaiEyVSmbSbxGs3hNo/G7xqdVDqmk95kZlCQ7kdJJXX3QazCPlw+Y9PTz45KihZ+Gx9YFPHAeEkZppQj99UAFKiqMGHxRbqIFqg9lNN5gy3kUSp5D3g26mrlp9JzhlNckEYpXR4QMx89u0jDBX0XHlyB0EYco4UhVyY7PyOCXMn6uIml7QkAE9hERsqoqecPyuooKQtjiEaTO6bWWNSVu3DLIEcB4fwZaVuyrDSJEA2MTIsCPCIU9RSiOZCYnfTV+q110miaVkGbtMLoSKT2dmsp8lUpjgmy5Uz0AuLmYiFDMahedeK9hEa7ZN6R7hIgoViSGYs8bNP6aRuN3m6NbV5VKqXVuv/w90+7m52FM/I763XzCPFCmNMYfd/bpQ2gTyYwcJ9MQ9gKHRWlMDKwdTVBESEFhtMHhBL78bWD2FeXTccvBWOYYjlmaILwM9s7BkHANAO4g4HTZW4NUisRBgGaE2fVNifP48ZJjUEMeoSpKY1HHEFacM1lXhOyQOmMcATAsRWicT8yycpPfySIRMpjPqYW+hYzfVhUyQTACWb281iaImWUyJT1CvaaH3UPCm2NpCK1QhMTNnwy5vnRcrGFhbIk4j1r04XAiJb05AeEbqqgIib1wsAzCGOSlMTJ+5yyax/NmryUMHqEG2h+LuUoxoSr1G0pjrRqV+caXV5VEajWNHDGeR8yYzVQhaVtLHETExQlU31BGb50PGYbyVtrXEYAiQgoKoxHz/jdw/o/tT54nkNLgcFenKrkDuioFAMd91f4aUhH61HxOdkA3gS49MM6+ImRI2h6M83/bMktzEuPMpXDjOROhJatQpozGcWBYilCjCIVsdNHkeTtJ2zUAdEWoiRQhq8RQEAw3SyEA0Q1kV2UTP9Na6ETIiaw+v82GItQkiBCVcdzURWfDPB7VEsjmmGzz9g5ZNH67/fz/CXjqeb9h9lokI0hdJc+UIGOBtE5AKFW6jsI3K3bRcUWoNqevQWW+ZtHRVtIfRKCkbYdYw6AqNcsuuijgDRV/vr+W/6ED4AueOADuE6I9bQtm9UiCav4w+4yhiJCCwucRdOMNNVVHpjRNn/sEAFMW2V+DbvK0TjXlNboJUBuxJ2RfVTL5jOLisRrrzzeqY6ak7WGUxuTsNfsjR2LCrN0kFSH7I0dIEZK5SlaJoTvAc6qgd461UK6SVVJHREjcuAFggk8Y0KFZjAIQIZeUNi5gaxiu2IsmN39tJvzjzqRFn5HhmDr0oW8oLQkZhT1WVoTE7LV0NwDGu8ZoWrtQeCqXxvhrRHM0e00PVKwXHW1ySGwpRKhMaCBTYo3GnHgv5UqFmiYN08d49PltpAhN9Irr3RMuTaZGEIoIKSh8HkE3pWpa5/Ph9FQ3IDH/l6/djjFAL3GQIlSNaXu4hmuHUxprMdRTXVmLuvDoJjuM0liLN4mbFkzB/InVJ22TImQ7YNIweJW61hrdNuMEqBzE+uEGVy7OGidYSKDOWneieK1ah06EXMjAkbLYNQZI8idVMQGHJEI1ldcI6EN9c0zPVfIM8Xb4yooQ3wtndggBJE1m6RCZxyv5rgQRolEpxvb5MJnHK60hSA5FBxg7zyh922qcwCQXP4/4QBqdffzaaBOt/VX5FY8CFBFSUPg8gjw9sRIDH+1g6gXVPS/aZvZ7DKc0lhLG4OEQoZ5P+SRvwH650DR7jciU/fZ5qQQNozSmDfXiunO+iBZShGzNXqsBoI8aqSpXidQYoQjVuyhg0uLPxl8LaPzWUyPI1MkxC1PnjRCvVe8ckg/VypZ+zdq5UBaRUydCYXcOGl1rNkaO1Dn0OAGf2wFnvyBClf4Q8QRlV2O91oP+Id2bQ+UyqwGTwZQgMSnd6Cw78SoRQ2Ecr8/p89voPKJp8V4qxQkIwzSNHIkPpKQiRB1tiggpKCgcPUw5D7jsMWDhv1S/xvk/Ab5wDvDVf63u+ZoGjJulf11VaSyvLDAcImQsr5WbQ1UM0uxc5cgRev5QL6/BkOrgtVNeMyhbAEBZRFUoQg1CCaoqV4lmr4HSugWZsvrzdTj0TjyRzXR8rTBc2yRCEcPctBaPoVRoRVUSRKfWoa8xTqQzWyZTYu8lKQUQC7iBhEVFCJDvuR596E9mZdeYN0lz06yVxvzpbriQMSlC+iw6a7lKNcLb1J/SiVBQmscr5SpxRagV/LzjA2nZPl/PKPFbESEFBYWjBYcTmHzu8KLsZy8FLl87vDVoMjUwPEVIrjEMItQlAt2s5NQUrDHMcSFeg18qPTAsj5AkQkM2SjgEcWy9IEKRanKVSAURihDFCtgjhvxYSuseTyTG6s9G/Dy8WQMRsquQifdcYyJThvdihUyJ12pw6mu0BwaF8qhZjBMQs9e0HvQOpWUGkOyiq0QOA/WysaEevRhIZuUa3pRFz5QgQqT+GA3X/gHq9qpAhCJcEYoJVSk+mJJddNEMBW6OPqM0oIiQgoLCZ4nWYSpC/joABrN3VcGO4nW7d/HPdsZayDUEYUkcBLI2zcEAL4FQF95QlUNojSnbxs921BxxLJXEwqzf9LgliJsqBTrKkEpbAZN8P2n2mpyZZrVkKa4DRzYJrwiHtK2QifccYnpZq5mIkNXzkLPXdCI0kcZrBOoAp7vYs8yQRKhXeow05KANWCRCDocpS8g4LNU1RETIWmnMn+mBFykc6E0im+O+LbdIi65YGhOvERJhm4f6DJEE0jxeYQbcCOFzT4S2bt2KE088UX74/X48/vjjI31aCgpjA0ZFyO6IDoB3iBlvStXMXiPyRP4gKzk1+ciPAoDGO2CsQtPM5bVhmKV1IlQFmRIqCJWDJAmoYvZaRAyhDTHKZrLfRdfmT+Mnl55kfz88YRBBpvOolQGTFozSgLyugjSWA3rJ0DIREqSr3jB7rUXEG1gmZIbSGOGYYFpP/LaUq8SJUKMWx/5e8k0xfZRLpXPx18puwAYtju2dfC/dTg1avyiNVSIxwjgezPCf5c5D/BpzaIBngNZQpbERwbHHHot3330X7777Ll599VUEg0HMnz9/pE9LQWFswGgWpWRZuzD+RVxVaawmb71qiBANoRVEyBuxNjXedB7iJt9/QCdlVZilkewFcrnqOuDEXhwTyuAfTp9o6PiqsbGGeQitj4a42vE7CcJz59faccHMVp0IWd0Ph0Mee2IjJ0R1RCRslsa8aT3PSJa4LCtCQh0zzF5rsN1Fx69HStoGgNkxuj5qAJen8hrCi9Sg9UiPUcyTgZYVnqdK5FDTJFmqQT92dfF9qPU5oNHPpiKZ4t/3Cl/S9oN8T+qCHgOZUkRoxPHEE0/g3HPPRTBoYUChgoLCkUHbKfzz9Iuqe/6wiVDec6opjRk7z4Dq/E4FI0c0GUJn7RzoNZl5flsVQ2gD2T6sOv94uOxMnicIUljjHERtwA1PpoqOvlJddFVEEqya34ZJsSAumMIVDevltRoAgDvTDwd4SUq249skQhGmqzn1Lps/F6H4yJEaAGbUiHKf1WuVSmM02gNAm9cwPsXGUN6oloCoiqFVpJgbv18SYi9c2SF4kZKZSnVBQ6p0NQOgjwJGnAi9/PLLOP/889Ha2gpN04qWre6//35MmjQJPp8Ps2fPxiuvvFLVaz3yyCO45JJLhnnGCgoKtnDZo8A/PANMXVLd84PDHEJ7JIiQVIT2VH8edJPvrVJVMgyQRc/uwnWtQPqM4vwzJW3bCpjk7/20Ni9+v/w0aNXECRi76IBhlQrb/Cm8eNOXcVyUSp9WS2M18p9kGpflNZulMTnsFUCts7pcpZhDJ0JfDA2Zvmd1jXrD2JJxXjKP19sayksZU4ChzOcJVfY7eSOAg4ed1hrKfOP8WT3+4kjkmn0GGHEilEgkMHPmTNx3331Fv/+73/0OK1euxHe/+11s2rQJZ555JhYtWoTdu/VfBLNnz8b06dMLPvbt2yeP6e3txWuvvYavfrWKUQEKCgrVwxcB2ufaLyURRoMilO8Rqsb47ctboxpVid5LXBAyT9ie94pu/oNx3sZfTeeZIDFBlsDkxpChrDUyI0fkeyBzsdXSmNMtgzLJ+B2BTSJEs9cyOgGJ2h3KKz1COplq99rtoqM4Ab2Nv9llIEJWIK6BVhkhADRJv5OF/dA0PZJA00uFE32UKh2qLlT1KMBmVv2Rx6JFi7BoUen4/nvvvRdXXXUVrr76agDAj370Izz77LN44IEHcMcddwAANm7cWPF1/vCHP2DhwoXw+Xxlj0smk0gm9Quhp4dfnL29vaWeoqCg8JkiDCSFVp92AXb/L6Y1/fkAkPNVsYZLrCFUh5zf/hpZH19j/07+ORewvwYL8efu/Zh/9gTtrZF2ivcxBHTsAESLNFIO6+vQGvFu/px4N//azhpZD3/O4UNijbhYw2l9jVyAP+fQfv6crk7+ddbGz1eLAMk+hNOHkctGwRKH0ZtkQM5rbY2MW+6nO9mNJLzAQDdfI2txjRy/LhzoRi7JyYtvsJOvoUUsngffzwZfErl+0RGYOiTWCFk8D76fdb5eeR7hTBdfg1lcA1Eg2Ylg6hByOaF0pcV+BGL2r/dhgu7bjLHyB7JRBABs7dq18utkMsmcTid77LHHTMddf/317KyzzrK19uLFi9kTTzxR8bhVq1Yx8LHK6kN9qA/1oT7Uh/r4K//Ys2dP2fv+iCtC5XDo0CFks1k0NZnrik1NTdi/f7/ldXp6evCnP/0Jjz76aMVjb7nlFtxwww3y61wuh8OHD6O+vh5atZPAi6C3txfjx4/Hnj17EIlUIZGPMaj9sg61V/ag9ss61F5Zh9ore/gs9osxhr6+PrS2tpY9blQTIUI+AWGM2SIl0WgUBw4csHSs1+uF1+s1PVZTU2P5tewiEomo/yQ2oPbLOtRe2YPaL+tQe2Udaq/s4UjvVzQarXjMiJulyyEWi8HpdBaoP52dnQUqkYKCgoKCgoKCXYxqIuTxeDB79mysW7fO9Pi6detw2mmnjdBZKSgoKCgoKHxeMOKlsf7+fnzyySfy6507d+Ldd99FXV0dJkyYgBtuuAGXX3455syZg7lz5+IXv/gFdu/ejeXLl4/gWQ8fXq8Xq1atKijDKRSH2i/rUHtlD2q/rEPtlXWovbKHkdwvTXRrjRjWr1+Ps88+u+DxpUuX4te//jUAHqh49913o6OjA9OnT8e//du/4ayzzjrKZ6qgoKCgoKDwecOIEyEFBQUFBQUFhZHCqPYIKSgoKCgoKCh8llBESEFBQUFBQWHMQhEhBQUFBQUFhTELRYRGCPfffz8mTZoEn8+H2bNn45VXXhnpUxpxrF69GpqmmT6am5vl9xljWL16NVpbW+H3+/HlL38Z77///gie8dHDyy+/jPPPPx+tra3QNA2PP/646ftW9iaZTGLFihWIxWIIBoO44IIL8Omnnx7Fd3H0UGm/rrjiioJr7dRTTzUdM1b264477sDJJ5+McDiMxsZGXHjhhdi6davpGHV9cVjZK3Vt6XjggQdwwgknyJDEuXPn4umnn5bfHy3XlSJCI4Df/e53WLlyJb773e9i06ZNOPPMM7Fo0SLs3r17pE9txHH88cejo6NDfmzZskV+7+6778a9996L++67D2+99Raam5sxf/589PX1jeAZHx0kEgnMnDkT9913X9HvW9mblStXYu3atVizZg1effVV9Pf3Y/Hixchms0frbRw1VNovADjvvPNM19p//dd/mb4/VvbrpZdewrXXXos33ngD69atQyaTwYIFC5BIJOQx6vrisLJXgLq2CG1tbbjzzjvx9ttv4+2338Y555yDJUuWSLIzaq4rW5NLFY4ITjnlFLZ8+XLTY8cddxz79re/PUJnNDqwatUqNnPmzKLfy+VyrLm5md15553ysaGhIRaNRtnPf/7zo3SGowOAeTixlb2Jx+PM7XazNWvWyGP27t3LHA4He+aZZ47auY8E8veLMcaWLl3KlixZUvI5Y3m/Ojs7GQD20ksvMcbU9VUO+XvFmLq2KqG2tpb98pe/HFXXlVKEjjJSqRQ2btyIBQsWmB5fsGABNmzYMEJnNXqwbds2tLa2YtKkSfjbv/1b7NixAwAP2ty/f79p37xeL+bNmzfm983K3mzcuBHpdNp0TGtrK6ZPnz5m92/9+vVobGzElClTsGzZMnR2dsrvjeX96unpAQDU1dUBUNdXOeTvFUFdW4XIZrNYs2YNEokE5s6dO6quK0WEjjIOHTqEbDZbMCutqampYKbaWMOXvvQl/Pu//zueffZZPPjgg9i/fz9OO+00dHV1yb1R+1YIK3uzf/9+eDwe1NbWljxmLGHRokX47W9/ixdeeAH33HMP3nrrLZxzzjlIJpMAxu5+McZwww034IwzzsD06dMBqOurFIrtFaCurXxs2bIFoVAIXq8Xy5cvx9q1azFt2rRRdV2N+IiNsQpN00xfM8YKHhtrWLRokfz3jBkzMHfuXHzhC1/Aww8/LM2Gat9Ko5q9Gav7d8kll8h/T58+HXPmzEF7ezueeuopXHTRRSWf93nfr+uuuw6bN2/Gq6++WvA9dX2ZUWqv1LVlxrHHHot3330X8Xgcjz76KJYuXYqXXnpJfn80XFdKETrKiMVicDqdBWy2s7OzgBmPdQSDQcyYMQPbtm2T3WNq3wphZW+am5uRSqXQ3d1d8pixjJaWFrS3t2Pbtm0AxuZ+rVixAk888QRefPFFtLW1ycfV9VWIUntVDGP92vJ4PJg8eTLmzJmDO+64AzNnzsSPf/zjUXVdKSJ0lOHxeDB79mysW7fO9Pi6detw2mmnjdBZjU4kk0l8+OGHaGlpwaRJk9Dc3Gzat1QqhZdeemnM75uVvZk9ezbcbrfpmI6ODrz33ntjfv8AoKurC3v27EFLSwuAsbVfjDFcd911eOyxx/DCCy9g0qRJpu+r60tHpb0qhrF8bRUDYwzJZHJ0XVdHzHatYBlr1qxhbreb/epXv2IffPABW7lyJQsGg2zXrl0jfWojihtvvJGtX7+e7dixg73xxhts8eLFLBwOy3258847WTQaZY899hjbsmULu/TSS1lLSwvr7e0d4TP/7NHX18c2bdrENm3axACwe++9l23atIn95S9/YYxZ25vly5eztrY29t///d/snXfeYeeccw6bOXMmy2QyI/W2PjOU26++vj524403sg0bNrCdO3eyF198kc2dO5eNGzduTO7XNddcw6LRKFu/fj3r6OiQHwMDA/IYdX1xVNordW2Zccstt7CXX36Z7dy5k23evJl95zvfYQ6Hgz333HOMsdFzXSkiNEL42c9+xtrb25nH42GzZs0ytV+OVVxyySWspaWFud1u1trayi666CL2/vvvy+/ncjm2atUq1tzczLxeLzvrrLPYli1bRvCMjx5efPFFBqDgY+nSpYwxa3szODjIrrvuOlZXV8f8fj9bvHgx27179wi8m88e5fZrYGCALViwgDU0NDC3280mTJjAli5dWrAXY2W/iu0TAPbQQw/JY9T1xVFpr9S1ZcaVV14p73MNDQ3s3HPPlSSIsdFzXanp8woKCgoKCgpjFsojpKCgoKCgoDBmoYiQgoKCgoKCwpiFIkIKCgoKCgoKYxaKCCkoKCgoKCiMWSgipKCgoKCgoDBmoYiQgoKCgoKCwpiFIkIKCgoKCgoKYxaKCCkoKFTE6tWrceKJJx6x9davXw9N0xCPx6t6/sDAAL7+9a8jEokMax0FBQUFRYQUFBRwxRVXQNM0aJoGt9uNY445BjfddBMSiQQA4KabbsLzzz8/wmep4+GHH8Yrr7yCDRs2oKOjA9FodKRPCQCwadMmXHzxxWhqaoLP58OUKVOwbNkyfPzxx0f0dTRNw+OPP35E11RQGKtQREhBQQEAcN5556GjowM7duzA97//fdx///246aabAAChUAj19fUjfIY6tm/fjqlTp2L69Olobm6GpmkFx6RSqaN6Tk8++SROPfVUJJNJ/Pa3v8WHH36I3/zmN4hGo7j11luP6rlYRTqdHulTUFAYeRzRgR0KCgp/lVi6dClbsmSJ6bGrr76aNTc3M8YYW7VqFZs5cyZjjM/+mTZtGlu2bJk8dseOHSwSibBf/OIXjDE+Q+iuu+5ikyZNYj6fj51wwgns97//vTyeZoF1d3czxhjbtWsXW7x4MaupqWGBQIBNmzaNPfXUU0XPdd68eaY5T/PmzWOMMdbe3s5uv/12tnTpUhaJRNjf//3fM8YY+8///E82bdo05vF4WHt7O/vhD39oWo+ed/nll7NgMMgmTJjAHn/8cdbZ2ckuuOACFgwG2fTp09lbb71Vcv8SiQSLxWLswgsvLPp9ep+MMfb++++zRYsWsWAwyBobG9lll13GDh48aHp/K1asYDfffDOrra1lTU1NbNWqVabzNb7/9vZ2+b0nnniCzZo1i3m9XjZp0iS2evVqlk6n5fcBsAceeIBdcMEFLBAIsNtuu63ke1JQGCtQREhBQaEoEVqxYgWrr69njJmJEGOMbdq0iXk8HrZ27VqWyWTY6aefbnr+d77zHXbcccexZ555hm3fvp099NBDzOv1svXr1zPGConQ1772NTZ//ny2efNmtn37dvbHP/6x5CDirq4utmzZMjZ37lzW0dHBurq6GGOcIEQiEfav//qvbNu2bWzbtm3s7bffZg6Hg33ve99jW7duZQ899BDz+/2mgaLt7e2srq6O/fznP2cff/wxu+aaa1g4HGbnnXcee+SRR9jWrVvZhRdeyKZOncpyuVzRc3rssccYALZhw4ay+7xv3z4Wi8XYLbfcwj788EP2zjvvsPnz57Ozzz5bHjNv3jwWiUTY6tWr2ccff8wefvhhpmmaHFbZ2dkpB312dHSwzs5OxhhjzzzzDItEIuzXv/412759O3vuuefYxIkT2erVq+XaAFhjYyP71a9+xbZv38527dpV9nwVFMYCFBFSUFAoIEJvvvkmq6+vZ3/zN3/DGCskQowxdvfdd7NYLMZWrFjBmpubparR39/PfD5fASm46qqr2KWXXsoYKyRCM2bMMN2wK+Gb3/ymVIII7e3tBYrMN77xDTZ//nzTYzfffDObNm2a6XmXXXaZ/Lqjo4MBYLfeeqt87PXXX2cAWEdHR9HzueuuuxgAdvjw4bLnfeutt7IFCxaYHtuzZw8DwLZu3coY40TojDPOMB1z8skns3/6p3+SXwNga9euNR1z5plnsh/84Aemx37zm9+wlpYW0/NWrlxZ9hwVFMYaXEe5EqegoDBK8eSTTyIUCiGTySCdTmPJkiX46U9/WvL4G2+8EX/4wx/w05/+FE8//TRisRgA4IMPPsDQ0BDmz59vOj6VSuGkk04qutb111+Pa665Bs899xy+8pWv4Otf/zpOOOEE2+9hzpw5pq8//PBDLFmyxPTY6aefjh/96EfIZrNwOp0AYHqtpqYmAMCMGTMKHuvs7ERzc3PB6zLGLJ3fxo0b8eKLLyIUChV8b/v27ZgyZUrB+QBAS0sLOjs7K6791ltv4V/+5V/kY9lsFkNDQxgYGEAgEABQuEcKCmMdiggpKCgAAM4++2w88MADcLvdaG1thdvtLnt8Z2cntm7dCqfTiW3btuG8884DAORyOQDAU089hXHjxpme4/V6i6519dVXY+HChXjqqafw3HPP4Y477sA999yDFStW2HoPwWDQ9DVjrMBIXYy0GN8rHV/sMXpv+SAC89FHH2Hu3Lklzy+Xy+H888/HXXfdVfC9lpaWoudDr1/qtY1r//M//zMuuuiigu/5fD757/w9UlAY61BESEFBAQC/QU6ePNny8VdeeSWmT5+OZcuW4aqrrsK5556LadOmYdq0afB6vdi9ezfmzZtneb3x48dj+fLlWL58OW655RY8+OCDtolQPqZNm4ZXX33V9NiGDRswZcoUqQYdCSxYsACxWAx333031q5dW/D9eDyOmpoazJo1C48++igmTpwIl6v6X79utxvZbNb02KxZs7B161ZbP0MFBQVFhBQUFKrAz372M7z++uvYvHkzxo8fj6effhp/93d/hzfffBPhcBg33XQTvvWtbyGXy+GMM85Ab28vNmzYgFAohKVLlxast3LlSixatAhTpkxBd3c3XnjhBUydOnXY53njjTfi5JNPxu23345LLrkEr7/+Ou677z7cf//9w17biGAwiF/+8pe4+OKLccEFF+D666/H5MmTcejQITzyyCPYvXs31qxZg2uvvRYPPvggLr30Utx8882IxWL45JNPsGbNGjz44IOWydnEiRPx/PPP4/TTT4fX60VtbS1uu+02LF68GOPHj8fFF18Mh8OBzZs3Y8uWLfj+979/RN+vgsLnCSpHSEFBwRY++ugj3Hzzzbj//vsxfvx4AJwYxeNxmZdz++2347bbbsMdd9yBqVOnYuHChfjjH/+ISZMmFV0zm83i2muvxdSpU3Heeefh2GOPPSJkZdasWXjkkUewZs0aTJ8+Hbfddhu+973v4Yorrhj22vlYsmQJNmzYALfbjW984xs47rjjcOmll6Knp0cSkdbWVrz22mvIZrNYuHAhpk+fjm9+85uIRqNwOKz/Or7nnnuwbt06jB8/XvquFi5ciCeffBLr1q3DySefjFNPPRX33nsv2tvbj/h7VVD4PEFjVl1+CgoKCgoKCgqfMyhFSEFBQUFBQWHMQhEhBQUFBQUFhTELRYQUFBQUFBQUxiwUEVJQUFBQUFAYs1BESEFBQUFBQWHMQhEhBQUFBQUFhTELRYQUFBQUFBQUxiwUEVJQUFBQUFAYs1BESEFBQUFBQWHMQhEhBQUFBQUFhTELRYQUFBQUFBQUxiwUEVJQUFBQUFAYs/i/BDwrVkM5rs0AAAAASUVORK5CYII=" + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(psf_170_65_lpi[750, 750:1050])\n", + "plt.plot(psf_170_20_lpi[750, 750:1050])\n", + "plt.yscale('log')\n", + "# plt.xscale('log')\n", + "plt.ylim([1e-7, 1])\n", + "plt.title('SunCET Diffraction Pattern for 65 vs 20 LPI mesh')\n", + "plt.xlabel('Pixels from Center')\n", + "plt.ylabel('Intensity') " + ], + "metadata": { + "collapsed": false, + "ExecuteTime": { + "end_time": "2024-04-24T05:10:56.451769Z", + "start_time": "2024-04-24T05:10:56.301800Z" + } + }, + "id": "203d573a11a98042" + }, + { + "cell_type": "code", + "execution_count": 39, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "0.18059616005355936\n", + "0.05062473026403846\n" + ] + } + ], + "source": [ + "mask = np.zeros([1500, 1500]) + 1.\n", + "mask[747:754, 747:754] = 0\n", + "print(np.sum(psf_185_65_lpi * mask))\n", + "print(np.sum(psf_185_20_lpi * mask))" + ], + "metadata": { + "collapsed": false, + "ExecuteTime": { + "end_time": "2024-04-24T05:29:32.098085Z", + "start_time": "2024-04-24T05:29:32.079725Z" + } + }, + "id": "841a93bcb180a540" + }, { "cell_type": "code", "execution_count": 5, @@ -224,8 +385,8 @@ " \"SIMPLE\": \"T\", \n", " \"BITPIX\": -32,\n", " \"NAXIS\": 2,\n", - " \"NAXIS1\": 1500,\n", - " \"NAXIS2\": 1500,\n", + " \"NAXIS1\": 2000,\n", + " \"NAXIS2\": 2000,\n", " 'WAVELNTH': 170.0, # NOTE: FITS doesn't have a keyword to store the units for wavelength, it's usually listed in the header comment instead\n", " 'WAVEUNIT': 'Angstrom',\n", " 'TELESCOP': 'SunCET',\n", @@ -249,7 +410,7 @@ "source": [ "for wavelength in range(170, 211, 1):\n", " print(wavelength)\n", - " psf = suncet_psf.psf(wavelength)\n", + " psf = suncet_psf.psf(wavelength, diffraction_orders=np.arange(-4000, 4000, 1), angle_arm = [90, 0] * u.deg, angles_focal_plane=[90, 0] * u.deg, output_size=[1500, 1500] lpi=20)\n", " header = make_header_template()\n", " header['WAVELNTH'] = wavelength\n", " if wavelength == 170:\n", diff --git a/suncet_instrument_simulator/Diffraction/suncet_psf.py b/suncet_instrument_simulator/Diffraction/suncet_psf.py index c291122..b125425 100644 --- a/suncet_instrument_simulator/Diffraction/suncet_psf.py +++ b/suncet_instrument_simulator/Diffraction/suncet_psf.py @@ -14,8 +14,7 @@ __all__ = ["psf", "filter_mesh_parameters", "_psf"] - -def filter_mesh_parameters(lambda_value, angle_arm=None, lpi=None): +def filter_mesh_parameters(lambda_value, angle_arm=None, angles_focal_plane=None, lpi=None): """ Geometric parameters for meshes in SunCET filters used to calculate the point spread function. @@ -77,11 +76,14 @@ def filter_mesh_parameters(lambda_value, angle_arm=None, lpi=None): psf_sigma = psf_80pct_px * 0.6178878 # Compute sigma from 80% encircled value if angle_arm is None: - angle_arm = [45, -45] * u.deg + angle_arm = [0, 90] * u.deg + if angles_focal_plane is None: + angles_focal_plane = [90, 0] * u.deg return { "wavelength": lambda_ang, "angle_arm": angle_arm, + "angles_focal_plane": angles_focal_plane, "error_angle_arm": [0.02, 0.02] * u.deg, "spacing_e": theta_arcsec/px_scale, "mesh_pitch": line_pitch_micron, @@ -93,7 +95,8 @@ def filter_mesh_parameters(lambda_value, angle_arm=None, lpi=None): } -def psf(channel: u.angstrom, use_preflightcore=False, diffraction_orders=None, angle_arm = None, lpi=None, use_gpu=True): +def psf(channel: u.angstrom, use_preflightcore=False, diffraction_orders=None, angle_arm=None, angles_focal_plane=None, + lpi=None, output_size=None, use_gpu=True): r""" Calculate the composite PSF for a given channel, including diffraction and core effects. @@ -188,16 +191,19 @@ def psf(channel: u.angstrom, use_preflightcore=False, diffraction_orders=None, a AIA PSF Characterization and Deconvolution `__ """ - meshinfo = filter_mesh_parameters(channel, angle_arm = angle_arm, lpi = lpi) + if output_size is None: + output_size = [1500, 1500] + meshinfo = filter_mesh_parameters(channel, angle_arm = angle_arm, angles_focal_plane= angles_focal_plane, lpi = lpi) angles_entrance = meshinfo["angle_arm"] - angles_focal_plane = u.Quantity([45.0, -45.0], "deg") + angles_focal_plane = meshinfo["angles_focal_plane"] if diffraction_orders is None: diffraction_orders = np.arange(-1500, 1500, 1) - psf_entrance = _psf(meshinfo, angles_entrance, diffraction_orders, use_gpu=use_gpu) + psf_entrance = _psf(meshinfo, angles_entrance, diffraction_orders, output_size, use_gpu=use_gpu) psf_focal_plane = _psf( meshinfo, angles_focal_plane, diffraction_orders, + output_size, focal_plane=True, use_gpu=use_gpu, ) @@ -213,8 +219,8 @@ def psf(channel: u.angstrom, use_preflightcore=False, diffraction_orders=None, a return psf -def _psf(meshinfo, angles, diffraction_orders, focal_plane=False, use_gpu=True): - psf = np.zeros((1500, 1500), dtype=float) +def _psf(meshinfo, angles, diffraction_orders, output_size, focal_plane=False, use_gpu=True): + psf = np.zeros((output_size[0], output_size[1]), dtype=float) # If cupy is available, cast to a cupy array if HAS_CUPY and use_gpu: psf = cupy.array(psf)