{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Simultaneous Wavelength Mode with MC\n", "\n", "In this example we perform a full radiative transfer calculation using the simultaneous wavelength mode of the MC model. Radiances at different wavelengths are calculated using common ray tracing, reducing computation time and greatly increasing correlation between wavelengths." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import sasktran as sk\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from sasktran.geometry import VerticalImage\n", "\n", "# First recreate our geometry and atmosphere classes\n", "geometry = VerticalImage()\n", "geometry.from_sza_saa(sza=60, saa=60, lat=0, lon=0, tanalts_km=[20], mjd=54372, locallook=0,\n", " satalt_km=600, refalt_km=20)\n", "\n", "atmosphere = sk.Atmosphere()\n", "atmosphere['ozone'] = sk.Species(sk.O3OSIRISRes(), sk.Labow())\n", "atmosphere['air'] = sk.Species(sk.Rayleigh(), sk.MSIS90())\n", "atmosphere.brdf = 1.0\n", "\n", "# And now make the engine\n", "engine = sk.EngineMC(geometry=geometry, atmosphere=atmosphere)\n", "\n", "engine.max_photons_per_los = 300 # cap the calculation at 300 rays per line of sight\n", "engine.solar_table_type = 0 # calculate single scatter source terms on the fly; no cache\n", "engine.debug_mode = 1234 # disable multi-threading, fix the rng seed for reproducibility\n", "\n", "# Choose some wavelengths to do the calculation at\n", "engine.wavelengths = np.arange(332, 337, 0.5)\n", "\n", "# And do the calculation in simultaneous wavelength mode\n", "engine.simultaneous_wavelength = True\n", "simultaneous_output = engine.calculate_radiance()\n", "\n", "# Repeat without simultaneous wavelength mode\n", "engine.simultaneous_wavelength = False\n", "separate_output = engine.calculate_radiance()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# calculate HR for reference\n", "engine = sk.EngineHR(geometry=geometry, atmosphere=atmosphere)\n", "engine.wavelengths = np.arange(332, 337, 0.5)\n", "radiance = engine.calculate_radiance()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAACPZUlEQVR4nOyddXxUV/qHnzMedxeSQHB3L9SoAKUKddu6bN237bbb3W5ta9tft+4tdYMaFAoUpzgEQoLEXScZP78/TpIGiExCQkK4z4f5MDPXzp3M3PeeV76vkFKioaGhoaFxKLquHoCGhoaGRvdEMxAaGhoaGk2iGQgNDQ0NjSbRDISGhoaGRpNoBkJDQ0NDo0kMXT2AjiQ8PFwmJSV19TA0NDQ0jhk2bNhQLKWMaGpZjzIQSUlJrF+/vquHoaGhoXHMIITY39wyzcWkoaGhodEkmoHQ0NDQ0GgSzUBoaGhoaDRJj4pBNIXT6SQ7OxubzdbVQ9HQOAyLxUJ8fDxGo7Grh6KhcRg93kBkZ2cTEBBAUlISQoiuHo6GRgNSSkpKSsjOziY5Obmrh6OhcRg93sVks9kICwvTjINGt0MIQVhYmDa71ei29HgDAWjGQaPbon03Nbozx4WB0NDQ0NBoOz0+BqGhKC8v56OPPuLGG2/s6qE0yRVXXMHMmTM577zzOnS///znP3nggQcA2LdvHzNnzmTbtm0degyNLiB7A+z/HQxm0JvV/4c9N4LO8OejpdeNn2uzugY0A9GDcLlcGAxN/0nLy8t55ZVXuq2B6CwaGwiNHkLBdlj3OpgC1GvpBun58+HxAFI9RJ2TROjq3hMN/6lt6xumyT+f6411D9PhjwbjU//cB+JGQmDMUTr5o4vmYupkrFYrZ555JsOGDWPw4MHMnz8fgA0bNnDCCScwatQoZsyYQV5eHgDTpk3jtttuY+LEiQwePJi1a9cCsHbtWiZOnMiIESOYOHEiu3btAuCdd97h/PPPZ9asWZx66qlUV1dz0kknMXLkSIYMGcI333wDwH333UdGRgbDhw/n7rvvBuDpp59mzJgxDB06lEceeeSwsX/66afccccdALzwwgukpKQAkJGRweTJkwF47LHHGDNmDIMHD+baa69FSsnOnTsZO3Zsw3727dvH0KFDWzzvxrT02dx7772MHTuWvn37snz5cgBqamq44IILGDp0KHPnzmXcuHGsX7+e++67j9raWoYPH87FF18MgNvt5pprrmHQoEGceuqp1NbWAvDiiy8ycOBAhg4dyrx589rxl9Y4KpTuhdX/B77h4BuqHn4R4B8FATEQGAfBCRCcCMG9IChBPQLjIDBe/R8UV/c6DoLi6x512wQlgF8kmIPUbERKcNZCbSlU5UJJBuRvg+z1sHc5pC2Apf+CvC1d/cl0CsfVDOLv321nR25lh+5zYGwgj8wa1OzyH3/8kdjYWBYsWABARUUFTqeTW265hW+++YaIiAjmz5/Pgw8+yFtvvQUoo7Jy5UqWLVvGVVddxbZt2+jfvz/Lli3DYDCwaNEiHnjgAb744gsAVq1axZYtWwgNDcXlcvHVV18RGBhIcXEx48ePZ/bs2Tz55JNs27aNTZs2AfDzzz+Tnp7O2rVrkVIye/Zsli1bxtSpUxvGPnXqVJ5++mkAli9fTlhYGDk5OaxYsYIpU6YAcPPNN/Pwww8DcOmll/L9998za9YsHA4HmZmZpKSkMH/+fC644IJWzxtodR2Xy8XatWtZuHAhf//731m0aBGvvPIKISEhbNmyhW3btjF8+HAAnnzySV5++eWGc963bx/p6el8/PHHvP7661xwwQV88cUXXHLJJTz55JPs3bsXs9lMeXn5EXwjNDqNqnxY+RKY/MHk1znHEAKEHtB7v43DCqv+CwPPgr6nga7n3HcfVwaiKxgyZAh33XUX9957LzNnzmTKlCls27aNbdu2ccoppwDqrjYm5s8p6oUXXgioC3RlZSXl5eVUVVVx+eWXk56ejhACp9PZsP4pp5xCaGgooHLrH3jgAZYtW4ZOpyMnJ4eCgoLDxvXzzz/z888/M2LECACqq6tJT08/yEBER0dTXV1NVVUVWVlZXHTRRSxbtozly5dzzjnnALBkyRKeeuopampqKC0tZdCgQcyaNYsLLriATz/9lPvuu4/58+czf/58du3a1eJ5A62uU3/cUaNGsW/fPgBWrFjBX//6VwAGDx7cMFtpiuTk5AYD0ngfQ4cO5eKLL2bOnDnMmTOn2e01uoiaUvj9RUCAJairR3MwJj/lctrxDVRkwYhLweTb1aPqEI4rA9HSnX5n0bdvXzZs2MDChQu5//77OfXUUzn77LMZNGgQq1atanKbQ1MfhRD87W9/Y/r06Xz11Vfs27ePadOmNSz38/vzburDDz+kqKiIDRs2YDQaSUpKajLPXkrJ/fffz3XXXdfi+CdMmMDbb79Nv379mDJlCm+99RarVq3i2WefxWazceONN7J+/XoSEhJ49NFHG441d+5czj//fM455xyEEKSmprJ169YWz7t+XC2tYzabAdDr9bhcroZtvKV++/p91LuYFixYwLJly/j22295/PHH2b59e7PxHI2jjMMKq18BRw0ERHX1aJpGb1QurfwtsOwpGHdD9x1rG+g5c6FuSm5uLr6+vlxyySXcdddd/PHHH/Tr14+ioqKGi6DT6WT79u0N29THKVasWEFQUBBBQUFUVFQQFxcHqLhDc1RUVBAZGYnRaGTJkiXs36+UfAMCAqiqqmpYb8aMGbz11ltUV1cDkJOTQ2Fh4WH7mzp1Ks888wxTp05lxIgRLFmyBLPZTFBQUIMxCA8Pp7q6ms8//7xhu969e6PX63n88ceZO3cuQKvn7e06hzJ58mQ+/fRTAHbs2MHWrVsblhmNxoNmW03h8XjIyspi+vTpPPXUU5SXlzd8LhpdjMsOa1+Hqrzuf8EVQsU57FWw9J8qmH6Mo90idTJbt27l7rvvRqfTYTQa+b//+z9MJhOff/45t956KxUVFbhcLm677TYGDVIznJCQECZOnEhlZWWD7/2ee+7h8ssv57nnnuPEE09s9ngXX3wxs2bNYvTo0QwfPpz+/fsDEBYWxqRJkxg8eDCnn346Tz/9NDt37mTChAkA+Pv788EHHxAZGXnQ/qZMmUJWVhZTp05Fr9eTkJDQsM/g4GCuueYahgwZQlJSEmPGjDlo27lz53L33Xezd+9egFbP29t1DuXGG2/k8ssvZ+jQoYwYMYKhQ4cSFKTcENdeey1Dhw5l5MiRPPHEE01u73a7ueSSS6ioqEBKye23305wcHCzx9M4SnjcsPF9KNqlgsfHCr7h4KhWLrHB50CfU47ZuIRoy/S8uzN69Gh5aMOgnTt3MmDAgC4aUduZNm0azzzzDKNHj+7qoRwzuN1unE4nFouFjIwMTjrpJHbv3o3JZOrqoXnFsfYdPSpICVs/hT2/KtfNsVib4HZCZQ4kjIHhF4PRp6tH1CRCiA1SyiYvONoMQuOYp6amhunTp+N0OpFSNszSNI5hdv8EexZB0DFqHODPuETOHyoDa+x14N9kZ89ui2YguhlLly7t6iEccwQEBGitZnsS+1bCti+VW0nXhnTT7ogQ6jysRbDkCRh/A0T06+pRec2x6RjT0NDomeRvhT/eVZXJ+h7UI8MvAoy+sPw5SF/cqIK7e6MZCA0Nje5BaSasflVVMhssXT2ajsccAIGxsHW+MoLO7i/zrhkIDQ2NrqcyT1VJmwN6TJFZk+hNKi6RtQZWPAfWkq4eUYtoBkJDQ6NrqSmFlS8qiYvuViXdGQgdBCVCdZGqlyhO7+oRNYtmILqIv/zlL+zYsaND9pWUlERxcXGDYqvGwTz66KM888wzHb7f559/npqamobX/v7+HX6MHo+9Gla9rNwtvmFdPZqji3+kmlEsfxYyl3bLuIRmILqIN954g4EDB3boPjUDcXQ51EBotBGXHda+BtWF6mJ5PGIOBP9o2PQRbPpAfSbdCM1AdDLNyX1PmzatITXT39+fe++9l1GjRnHyySezdu1apk2bRkpKCt9++y2g5DVuvvnmhv3OnDnzsJTYQyW9m5P+3rdvHwMGDGhS9jojI4PTTjuNUaNGMWXKFNLS0gDYv38/J510EkOHDuWkk07iwIEDgGr001hio/4uOi8vj6lTpzJ8+HAGDx7cIM1dz9q1axuE97755ht8fHxwOBzYbLYGWfHXX3+dMWPGMGzYMM4991xqamqoqKggKSkJj8cDqBqIhIQEnE5ns2NvTHPrXHHFFdx6661MnDiRlJSUhnPyeDzceOONDBo0iJkzZ3LGGWfw+eef8+KLL5Kbm8v06dOZPn16w/4ffPBBhg0bxvjx4xtEEj/77DMGDx7MsGHDDhJDPK5xu2DDe1C8GwJiu3o0XYvBrOIS+1bCiueVy62bcHzVQfxwn0qj60iih8DpTza7uCm570OxWq1MmzaNf//735x99tk89NBD/PLLL+zYsYPLL7+c2bNnezWUQyW9m5P+BpqVvb722mt59dVXSU1NZc2aNdx44438+uuv3HzzzVx22WVcfvnlvPXWW9x66618/fXXzY7lo48+YsaMGTz44IO43e7D7rRHjhzJxo0bASUlPnjwYNatW4fL5WLcuHGAUm695pprAHjooYd48803ueWWWxg2bBi//fYb06dP57vvvmPGjBkYjcZmx96YltbJy8tjxYoVpKWlMXv2bM477zy+/PJL9u3bx9atWyksLGTAgAFcddVV3HrrrTz33HMsWbKE8PDwhr/j+PHjeeKJJ7jnnnt4/fXXeeihh3jsscf46aefiIuL06TEoa5K+jPIWX/sVkl3NEKn6iUq81R/iXHXQ1jvrh5V5xoIIcRpwAsocfU3pJRPHrK8P/A2MBJ4UEr5zCHL9cB6IEdKObMzx9pZNCX3fSgmk4nTTjutYX2z2YzRaGTIkCENctTtoSXp76Zkr6urq1m5ciXnn39+wz7sdjXlXbVqFV9++SWg+j7cc889LR57zJgxXHXVVTidTubMmdNwrHoMBgN9+vRh586drF27ljvuuINly5bhdrsbPqNt27bx0EMPNYjnzZgxA1AaT/Pnz2f69Ol88skn3HjjjS2OvZ7W1pkzZw46nY6BAwc2fE4rVqzg/PPPR6fTER0dfdBs4VBMJhMzZ85s+Ex/+eUXACZNmsQVV1zBBRdc0DBrOm6REnb9ABm/QkiSZhwaI4QSJLRVwLJnYMQl0Gtil35GnWYg6i7u/wVOAbKBdUKIb6WUjSOzpcCtwJxmdvNXYCcQ2CGDauFOv7NoSu67vsFOPUajsUHiW6fTNUhS63S6Bklrg8HQ4FYBmpTwPpSWpL+bkr32eDwEBwc3zEBaon68jcclpcThcABKBXbZsmUsWLCASy+9lLvvvpvLLrvsoH1MmTKFH374AaPRyMknn8wVV1yB2+1uCChfccUVfP311wwbNox33nmnwaU2e/Zs7r//fkpLS9mwYQMnnngiVqu11bG3dn6NP5N6jbK2aJU1/js2liN/9dVXWbNmDQsWLGD48OFs2rSJsLDjLCBbz/6VsONr1b1NaB7uJrEEKbfTH+9C+QEYfC4YukY6pjP/QmOBPVLKTCmlA/gEOKvxClLKQinlOuAwPWYhRDxwJvBGJ46x02lK7rs9JCUlsWnTpgZp6vpWpI05VNK7Oenv5ggMDCQ5OZnPPvsMUBfHzZs3AzBx4kQ++eQTQBme+pajSUlJbNiwAVCxhHpp7f379xMZGck111zD1Vdf3eR5T506leeff54JEyYQERFBSUkJaWlpDcqtVVVVxMTE4HQ6+fDDDxu28/f3Z+zYsfz1r39l5syZ6PX6Fsfuzfk1x+TJk/niiy/weDwUFBQcFPc59PNujoyMDMaNG8djjz1GeHg4WVlZrW7TI8nbAhvfUzEH3fHl3W4zBotKhd37m0oBri3vmmF04r7jgMa/hGxgXBu2fx64BwhoaSUhxLXAtQCJiYltG+FRoCm57/YwadIkkpOTGTJkCIMHD2bkyJGHrXOopPe9997bpPR3S3z44YfccMMN/OMf/8DpdDJv3jyGDRvGiy++yFVXXcXTTz9NREQEb7/9NgDXXHMNZ511FmPHjuWkk05qaF60dOlSnn76aYxGI/7+/rz33nuHHWvcuHEUFBQ0BG6HDh1KZGRkw134448/zrhx4+jVqxdDhgw56GJc35Co8QW7ubF7c37Nce6557J48WIGDx5M3759GTdu3EFS4qeffjoxMTEsWbKk2X3cfffdpKenI6XkpJNOavF4PZaSDFjzv7oqaXPr67dEaSakfa/kwHV6ZWz0RvV/40dT7zX1vt4Auqa2P+T9o+3q0emVkSjPUnGJ8Tcot9xRpNPkvoUQ5wMzpJR/qXt9KTBWSnlLE+s+ClTXxyCEEDOBM6SUNwohpgF3eROD6Aly3xrdj+rqavz9/SkpKWHs2LH8/vvvREdHd9j+e/x3tDIXlj2t7orNR+gtLj+gZMB9QlSA2+MCj1MZC7ez7nWjR8N7brWe9LR+jOYQ+kZGQ6/+NwdAn5M7v4ajthwcVTDyMkgY16HGqqvkvrOBxl0+4oFcL7edBMwWQpwBWIBAIcQHUspLOniMGhqtMnPmTMrLy3E4HPztb3/rUOPQ46kphd9fAGE4cuNQmQvbPlc++qEXtk+SQ3rqDEe9YWlsTFyNDM6h7x2yrP79igMqVtD3NIjs2Lqmg/AJBqMF1r+lZhSD5hwVMcPONBDrgFQhRDKQA8wDLvJmQynl/cD9AI1mEJpx0OgSNAn2dmKvUv5zl+PIC+GqC5TIndEXhs5rv16T0KnqZX0HBX3tVbDjG9j5LVRkQ+8TOy++Uh+XyFisGhGNvrLTpUk6LUgtpXQBNwM/oTKRPpVSbhdCXC+EuB5ACBEthMgG7gAeEkJkCyE6JmNJQ0Oj63DaVMyhpuTIjUNNCWyZDzoTDJun3DrdBXMADLsQ4sdC7h+w6UOVptpZ1MclSjNh6b+Vy60T6dRUAinlQmDhIe+92uh5Psr11NI+lgJLO2F4GhoanYHbBRvegdIMCDzCXtK2ctiisucYNg8swUc4uE5Ap1czh6A4SFsIG96G/rM6r9BNCAiIgdoyZSRGXQnxozoliK4lImtoaHQcHo9yBeVuVMbhSC5a9irY/Am4HTB0bvcX8wvvB6OuAEsgbPtMpageSVC8NXxC1Gey7jXY/lWniP1pBkJDQ6NjkBJ2LYDM3+oK4Y7AODhq1MzBWQND5oJ/VMeNszPxCYHhl0LMMDiwSp2Do7rzjmf0US6n9J87RcNJMxBHgUNloBsL7z366KPExcUxfPhwBg4cyMcff9wVQ9TQOHL2LYed3ylNoSOpknbZ1CzEVgGDz1Nd2I4l9Eboezr0O1NlXm14u3NjBTp9p/Xu1gxEN+D2229n06ZNfPPNN1x33XUN1cgaGscMOZtg44dHXiXtdighP2sRDDxbzUSOVaKHwIjLQG+GzR/DgdXdsudDS2gGohuRmpqKr68vZWVlXT0UDQ3vKd6j/OD+R1gl7XHBti/UXfeA2d1CzfSI8Y+EkZdDRD/YuxS2f3FM9KKu57gSRLntttu8EqJrC8OHD+f5559vcZ3a2tqD1ExLS0ublPD+448/SE1NJTLyOG2eonHsUZGtOsJZglSNQnvxuFWgtXw/9J8JEa3LwhwzGMww4CwIilcqtn+8rWZHAd2/4PK4MhBdhY+Pz0GG6Z133qGxJMh//vMfXn/9dTIzM/nxxx+7YIQaGu3AWgwrX6qTnDiC8iXpgbTvVFps6qkQNbjjxthdEALiRisX3I6vYeP7SqIjZni3ljw/rgxEa3f6XcXtt9/OXXfdxZdffslll11GRkYGFoulq4elodE8tko1c3A7lABfe5ESdv8ARWmQMh1iDxeh7FEExqq6hZ3fQfpPagbWd0bHVXZ3MFoMohtxzjnnMHr0aN59992uHoqGRvM4bbDmVbCWHLlxyFikujz2mqRE6I4HjD4w5HxImgqFO5SWk7W4q0fVJF4ZCCHEZCHElXXPI+r0lTQ6gYcffpjnnnvuoOZAGhrdBrdLCcaV7VPVvEfC3mWQswHix0CvyR0yvGMGIVS3uKFzwVmrjETB9q4e1WG0KvcthHgEGA30k1L2FULEAp9JKScdjQG2BU3uW+NY5Jj6ju7+CbZ9eeS9pPevhH3LlA8+dUa39sN3OvWCf5XZEDsCep/U9lThymw4+THwa3u1eUty397MIM4GZgNWACllLq008dHQ0OiBVBcq33lg7JFd0LPXK+MQOUgFpY9n4wCNBP/GKYmSTR8oDapugDcGwiHVNEMCCCH8OndIGhoa3Q4plS6STn9kAdW8zSruEJYK/c/U+lLXo9ND7+kw6ByoKVNihyV7unpUXhmIT4UQ/wOChRDXAIuA1zt3WB1LZ3XN09A4Uo6Z72bOH8pH7ncEmkiFO2H3jxCSDAPP0oxDU4T3VVlOliDVHClzaecK/rVCq44uKeUzQohTgEqgH/CwlPKXTh9ZB2GxWCgpKSEsLKyh17GGRndASklJSUn3T2m2VyupCL+I9ruDitNVrUNQnLpLPsKmOm6PxOpwtbpes6NtYoFofu0Wd+Fj0qPryGuLTzCMuBT2LIKs1ao50ICzwOzf6qYdTat/pbqMpeX1RkEI4SOESJJS7uvswXUE8fHxZGdnU1RU1NVD0dA4DIvFQnx8iy1Rup6d3ylV1fbKbZftU8Vh/pEw+PwOaZW5M6+SnPJadB12Xa7fUdtmdB4JUYFmBscFYdB14IxIZ1BtTIPiVWLAhrdh4GyVHHAU8caMfwZMbPTaXffemE4ZUQdjNBpJTtaycjU02kVJhnJzBLWz8U9FttJX8g1Vst1HotVUv8taJznltQRYDF3vFZCSwio7Gw+UMzwhGKO+g91mUYOV1Pn2r1QMKHkqJIw/aoF9bwyEQUrpqH8hpXQIIbpn2Z+GhkbH4XbCxg+UP7w9ctJV+UqZ1eyv8v2NPkc8JIlkd0EVQZ4K/BxuRMMdv0TUPZUCJDqk0NX9L1p4LRreB9H2C68QBJgNlNc4Wb+vlBGJIViMHSy97RehBP92/6iaEFVkK72qDvg8W8MbA1EkhJgtpfwWQAhxFtA9y/40NDQ6jszfoCqvfW4Na7HqI20ww9B5YOoY/3lxtQNnRR5G32A2Rp6BQKKTLnTSXfdwofPYMXic6KUDvceJweNA53Ggl66G9/TSic7jxCCd6KQLg8eJkG6kEH96mYQOZZKU0fjTdKj3RF1yp80YiL/FD6vdxdp9pYzqFYKfqYNVjAxmpXAbFA8Zi+tcTmdD4BEWK7Z2WC/WuR74UAjxMuozygIu69RRaWhodC3VRXVxg3YojtaWqU5qOp3K77cEdciQ3FKSdWAvHnMQ6+Ivx27smP02ID0HG5o6wyOkG710IQ5apv43uG30LVmM2V0IpghqnW7W7lVGItBy5LGWgxAC4kapCvYdX6t6id4nqeK6TsKbLKYMYLwQwh9VeV3VaaPR0NDoeqRUd/9C3/aYga1SZTx53DD8ItWCs4MozttPmceHXclXdLxxABA6PEKHh7Zd2Iv9+zKo4FsirLvRmSKpcQvW7i1lZK8QQn07wRtfL/iX9j3s+RkqsiBmaMcfB++ymMzAuUASYKgPCkkpH+uUEWloaHQteZuUgF5bXUsOq5o5uOxq5uAX0WFDclbmsavCwOZelyE6wzgcAXZDAH/EXkhi+Vr6Ff2ESe9LpT6A9ftKGR4fTGRgJ6QxG31UO9as1UrTqioXpj/ULqmNlvAm5P4NcBbgQslt1D80NDR6GvZq2PQR+IW3LWDrrFXGwV6llEo7shlOdT4HrGZ+Cb0Q0YEzkg5F6DgQMp7Vidfi0psJcxfgYxD8kVVOdllNJx1TQOIEFeMJ7gWmjhe58CYGES+lPK3Dj6yhodH92LVAzQSCQr3fxmWHrZ9CTakyDkEdWNdRlU+NIZB39TMICArvuP12ElWWGFYlXEv/4p+Jr1iP0RDKtpwKHG4PyeF+bS7G84qQXqDv4MypOryZQawUQgzplKNraGh0H0r3QsYS1fXMW9xOJQlRlQ8D50BIUseNpyof/MJYEHABDkNgx9cYdBJuvYXtkbPYFDMXX1lDlL6S3fmV7C6oOnakVerwZgYxGbhCCLEXsKMymaSUsnOiIhoaGkcft0vVPJgCvK958LhUAVdFlkrBDE/tuPFU5YFfONkDr2H18gJigzo/579DEYKCgEFUmmMYUvAV8TUH2F/kweGSDIwNRN/VBX5e4o2BOL3TR6GhodG17FumNH+8DUxLD+z8Fsoyoe/pEDmw48ZSlQ9+4ciJt/L16iJ8jHp0HaepcVSpNYWyLv5ykkuW06d0KcUlNja7PQyN72Bpjk6i1RFKKfdLKfcDtagSkgbpbw0NjR6AtVjNBLyteZBSxSqKd6s8/JhhHTeWqjzwj4BJf2VHmSCz0Eqo37Et3CCFgczw6ayLvxJ/ix5nWTYb9pXicHf/rpHepLnOBp4FYoFCoBewExjUuUPT0NDodKSELZ8CwruaBykh/Wcl/Z00VbUL7Sgq8yAgEibeitPoz9eb0gj2NR6mt7Rv5ya+euUfuBx2jCYLBpMZo9mM0WjGaK57XfcwNP7fbDn4fbMFo9HU9DZmC0ajGb3x8OO3lzLfJFYnXs+AwgUElW7hjwwHw5Kj8OloaY4OxBsX0+PAeGCRlHKEEGI6cKE3OxdCnAa8AOiBN6SUTx6yvD/wNjASeFBK+Uzd+xZgGWCuG+PnUspHvDslDQ0Nr8nbrOoegpO8Wz9zKeRtVIJxiRM6bhyVuRAQBRNvBUsg6/YUU1LtICHE96DV0jeu4s1HbsQvMJi43gNwOu24HA5qqyqpdNpx2e04HXZcDjtOhw2nw4HH3boseHMIIeqMhwWDydRgOIxmZUiCwiI5/YrbiYjzzjXnMPizOeZ84vz6kJzzHdt2VzGgT2/8zR0szdFBeDMqp5SyRAihE0LopJRLhBD/bm0jIYQe+C9wCpANrBNCfCul3NFotVLgVmDOIZvbgROllNVCCCOwQgjxg5RytTcnpaGh4QWOGtj0Ifh62eeheDdkr4HYkZB8QscpilbmQmA0TLgFLIFY7S4WbM0jwv/gGc2ONUt557FbCI/txfVPvkVgWKRXu3e7XXUGw47TbsfltOO0K+PhdNhw1RsUux1n3TKXw6HWr1vudDQyPI32kbZ+BTtWL2X2dfcy4cx53s02hI6coFFUWOLpe2A+u3Ztp0/vfgT5db++IN4YiPI6mY1lKE2mQlTRXGuMBfZIKTMBhBCfoAruGgyElLIQKBRCnNl4w7oWp9V1L411Dy3uoaHRkexa6H3Ng9MG6T+BX6SKO3SUcaiqMw4Tb1W9mYHfdhfhcHmw+P/petm4dCEf/vtu4nr359p/voFfoPcFc3q9Ab2PAbNPxxeSlRfl88mzD/D5i4+ybdWvzLvjCa8NV7U5is29rycu9yfYs4bkxF6EhrSh/uQo4I2BOAuwAbcDFwNBgDcyG3EoYb96soFx3g6sbgayAegD/FdKucbbbTU0NFqhbL/qWBboZVFbxmI14xh8fvukv5uiMgcC42DizQ3GoaTaztJdhUQG/Dl7WP3DZ3z2/N9IHjSKvzz+PzxGCz9uz6fW4W5YRyKp+1f/RuP/GpY3fg9USOWgfRzyXsP6stHyRgt8TXrm/e3/2Lb4c757/SmeunYW5936CMNPOMOrj8CtM3EgfhYlvr2xZX1JP1sl4dG9jlq/h9bwRqyvsazGu23Yd1Nn6PUsQErpBoYLIYKBr4QQg6WU2w47iBDXAtcCJCYmtmF4GhrHKQ01D/7eXexLM6Fgq4o5dJSERmWOaj864ZaDWmn+uC0fgcBQVxS37Kt3+fr//km/0ZO58uGXMVl8WLSzgPSCKiID/nTJNL6eikZP6iuXdYiDcjbFIcvFIRvXv27sMvrzvT/3k11Wy2d/5DBn+nn0HTGBj566l/eeuJ1tKxdzzs0P4xvgnW6UNXQga32iKcn6knHOXUTGJSP0R95c6Uhp1kAIIVZIKScLIao4+MJeXygX2Mq+s4HGbajigdy2DlBKWS6EWAqcBhxmIKSUrwGvAYwePfrYd0NJ2W3uHjR6KPtXqOI2b2oeXHbVqMY3DHpN6pjjV2arDnUTbj7IOBwoqeGPA2XEBvsgpWTRR//HD+++wJDJp3Lpfc9iMJkorLKxPbeSEYnBTE3tODHA9lJSbefrTbl8viGb2cNiueX5j1n88f/4+cNXyNi6jnl3/ot+o7z83HxCSUu+nLzc3zgzay1xUZFdrj3VbB2ElHJy3f8BUsrARo8AL4wDwDogVQiRXNeBbh7wrTeDEkJE1M0cEEL4ACcDad5se0yTsQQW3AkrX4b0RVCYBraKw+e8GhrtpaZUtQD1j/Ju/b1LwV4J/c5QfZKPlMpsCEo8zDh4PJJvNuXgY9QjgO/feJof3n2B0SefxWUP/geDyYSUkuW7i7EYdYxL6h6++jB/MxeMjsfPrOerTTnsK7Vx6iU38dcXPsHs48f/7r+KL//7OA5brVf7MxqNVMafyDv680gvtuGpyFZFiV1ESzOIFv8CUsrSVpa7hBA3Az+h0lzfklJuF0JcX7f8VSFENLAeCAQ8QojbgIFADPBuXRxCB3wqpfze+9M6BqnIga2fg08olB+AwkbJXuZACEuFiL5qWh4Qc1TaDWr0MKSELZ+hah68yJgp3w+5G1WtQ2DckR+7KlcZh4k3H6Y8ujOvkr3FVmKDzHzx0qOs/P4TJs26iLNv+hu6uorjzGIr2eW1TOsXgbkb1Q4EWIycPzqBbzflsmBLHicNiGRQ3yHc8cpXLHjrWZZ/9R67N/zORfc+RWK/1hWKDHoduqhU3iwNY47jd8aW7kYfGN0lv3nRnHhUnfaSRLmUEoGyuufBwAEpZfJRGqPXjB49Wq5fv76rh9F2XA5Y9hTUlCmZ5cZICW67kmF21aDspUdVvYb3Vfo3AdHqjlDfwR2sNHoWeVtg1X+Va6k1N6bbAevfUs9HX31k3y0pVcwhpBdMuOkw4+BweXj6pzQcTicLXn6YDYu/5cS513DmVXc2xABcHg8frD6AXie4eGxit5TecLo9LNiSx/7SGib2DmN0rxCEEOzeuIpPnr6PytIiTr7oBk656Hr0htY/TykluWU1TPLZyyz3rxiMpuZ7bFRmw8mPtasfhBBig5RydFPLmp1B1BsAIcSrwLdSyoV1r09HuXw0OordP6oZRHATQXZRd7fX+I5PSnDWwIG6ZiFCp9YLToDw/hCapBQ5fcNU20cNjYaahzDvYlx7l4OtHIZd1EHGIQkm3Nhkz4J1+0ooKK9mySsPsm3lIs648nZOvvD6g9bZklVBRa2Ts4bHdkvjAGDU65g1LJZfdhSwMqOEWoebKanh9B0xgbtf+46vXvkHP3/wMjvXLuWie54iKrF3i/sTQhAb4suqyt4U+kRzuflXLOUHVEe5jnD3eYE3RxkjpWz4a0kpfxBCPN6JYzq+KMlQ+ehtmcILoX5ojX9sHjfUlKvURVmX/qc3QWgKRA5QQcHAGOWu0oLgxx+7f1TNfIISWl+3Igdy1qlex03dtHiLlFBVbxxuApPvYatU2118tS6TH569gz0bV3L2jQ8xZc6lB61T43CxZl8pvcJ8SQrr+FqGjkSvE8wYFIWPSc/GrHJqnW5OHhCFj38gF93zFIMmnMTnLzzMszeezcyr72LyWZc0uNCaQghBTJAPeysFLzCbGxJ3Ebh/kTL0danBnYk3BqJYCPEQ8AHK5XQJUNKpozpecNbC+rfBHHTkdwQ6PVgC1aMetxMqsqFoFw2JaOYAFc8I7wdBsepuRItn9GzKD0D6L+pv3RoeF+xeqG4kkqe1/5gNxiEFxt/QpHEAWLh+Dx///Xry0rcw785/MnbGuYetsyazFKfb0y2ylrxBCMHU1HB8TXo1k3C6OXNIDEa9jmFTZpA8aCTzn3uQr//vCbav/pV5d/6LkMiYFvcZFWihuNrOs9kDuGlYKpG7PqyTJ4np1Bu+ZmMQDSuoYPUjwNS6t5YBf28tSN0VHHMxiE0fw77l3t3VdQRNxjOk8muG91WPwBgtntGT8Ljht6eUYuuh8a2myPwNslbBkAvU7LM9SKl84qG9WzQOu/blcMKJp1B0IJ1L7n+G4VMP7yxQUm3nwzUHGBIfxPR+3lUodye25VTwa1ohUYEWZg+PbRDmk1Ky+ofP+ObVf6HT6znnpr8x6qTZrUp1lNU4cLol140LJ2n/l5C/RRmJmqJOiUG0aiCOJY4pA1GwHX5/QWV1dFRlanuQUs1kHFUqMIlQdySDzoU+07tuXBodQ+YyFXvwJjBdlQ9/vAvRg6HfmS2v2xz1xiGsjzIOzcxOc3NzGTt5GgU5B7ji4ZcYOPaEJtf7emMOeZU2rpiQhI+p+2QutYWMomp+2JZPkMXInBGxBFj+vPkqzj3Ax0/fy97tfzBsygzOvfVR/FuRPqmyOam2u7hiQiIDbZtg62fqt3vavzvcQLQawayrSXhaCLFQCPFr/aPNo9D4E1slbHgHfMO71jhAXTzDV80aghJUP2G/CNWAPndT145N48ioKVXtQAOiWzcOHrfq8WDyg5ST2nc8L43D3r17mTBpMoV5OfzlH681axz2FVvZX1rDuOTQY9Y4APSO8GfO8Fiq7S4+XZ9NqdXRsCw8NpGbnvmAM6++k22rfuXp62azY83SFvcXYDES5GPkzd/3s14/DKbdD7HDvZNrbyPepLh8iCpSSwb+DuxDFcFptId6/X1n7VEJMrULvQn8I2HdG1C2r6tHo9EepFQFcVJ6V/NwYBVYiyB1BhjboSpabxzC+7ZoHNLS0pgyZQrFJaVc8vfX6Tt8fJPruT2SZelFBPsYGRYf3PbxdDPiQ3w5b1Q8Hin5bH0W+RW2hmU6vZ6T5l7LbS99jn9QCG/87To+ff5h7LXWZvfnazIQEWDmwzUHWFpgRo6/6aDCw47CGwMRJqV8EyX7/ZuU8ipUfwiN9pC9AbLXKr9hd8boC0Y/lTdvLe7q0Wi0lYLtkL3OO+2k6kI4sFK1DW1PX+kG49APxl3frHHYtGkTU6dOxeZwcvZDrzF4+Khmd7ktp4KyGieTU8PRd9O01rYSEWDmgtEJmI16vvgjm30lBxuAuN79uf2lL5h+/tWs+eFTnrlhDnu3/9Hs/swGPTFBFr7ZmMvCrfl0RrjAGwPhrPs/TwhxphBiBEpXSaOt1JTCpg+UO0ccA/UJPsHgdsPqV5QstMaxgbNWfc98w1r/nkmPyloyWKBPO8qbDjIO1zU7+1i1ahXTpk3DYrFwyeNv02/A4GYDsjanm9WZJcSH+JAS3r3TWttKkI+R80fFE+Jr4rvNuaTlVx603GAyMeuae7jpmQ+Qbjcv33kx37/5LC6no8n9GfU64kJ8WJxWQHmNs8l1jgRvrlL/EEIEAXcCdwFvALd1+Eh6Oh6PUtCUHnV3fqzgHwFVBaqq1t3xX0CNTmD3zyrOZfZCMi1rrQpO9zml7d9LKaEyCyIGtGgcFi9ezCmnnEJERATPvfc1+pBY/C3Np3Wv3VuKzaXSWjuq3Wd3ws9s4NxRccQG+fDT9gI2Hig7bJ2UIaO563/fMvbUc/h1/ms8f8v55O7d1eT+9DqBQSc6pWFOqwZCSvm9lLJCSrlNSjldSjkK1QlOoy3sWwEF27xvDN+dCIyD/O0qW6IHZb31SCqyVVGcNy7MmhKVZh3eFyL6t+04DcZhIIy9plnj8N1333HmmWeSnJzMD7/8yh+lhsM6xTWmrMbB5uxyBsUGEhHQ9XLXnYXZoOes4bH0jvBjWXoxKzOKD3MRWXz9mXvHE1z99/+jsrSI/9x8Lr9++gYet7uZvXY8zRoIIYReCHGhEOIuIcTguvdmCiFWAi8ftRH2BCrzYOunSv7iWLwjqpfxyFgCexZ39Wg0msPjho0fqhhAa4WXUsKuH1S9S59T2/a9rHcrRQ6Ccdc2axw++eQTzjnnHIYOHcrSpUvZWaHH6ZYtCu2tSC/GoNMxIaXt6ZrHGga9jjOGxDA4NpB1+8pYnFaIx3P4DdigCSdyz2vfM3DsNL5/42leufsySvKymthjx9PSDOJN4C9AGPCiEOJt4BngKSnliKMxuB6B26lSWvWmTklDO2oInUqB3foZ5Gzq6tFoNMX+Vaq5j68XBXG5G9RFvvdJbct+kVL1koismzk0851+4403uOiii5g4cSKLFi3CY/Jn2e6igzrFHcqB0hoyi62MTgrBz3x0tIa6Gp0QnNg/krFJoWzPrWThtjxc7sPlvf2DQ7ni4Ze48K4nycncyTPXz2bNj593SmC6MS39FUYDQ6WUHiGEBSgG+kgp8zt1RD2N9EVKNjmoB3S7a0h/fR187oLQbifoe/xSWwbbPlMuzNZmA7XlqmI6NAWiBrftOFV5rRqH//znP9xxxx2cdtppfPHFF/j6+vLeyn0YdH92ijsUj5QsTy8i0GJgREJwq8Mor3FQ62yfq+WwS2q9ZnXj/w+l7v1DF0mp9DBjAi3tjpcIIZjQOwwfk57fdhfx9aZcZg2LwWzQH7bemFPPpvewsXz8zH3Mf+5Btq1azAW3PQ66zgnmt2QgHFKqThVSSpsQYrdmHNpI2T7Y+e2x61pqCqMvmB0q/XXafd7JN2h0LlLCti9VIkRrNQxSwu4f1Pcx9bS2fS/t1coojLqsSeMgpeTxxx/nkUce4dxzz+Wjjz7CZDKxr9jKpqxy4kKa1/zakVtJcbWDMwZHN2tE6nF5PFjtLq6YlIyvSY8QqnGoqGshWn9KQtDkssbtRutfH7pevaBAU8sa79/plny0Zj+78quIC/Y5oqD68IRgfIx6ft6Rz+cbspkzPK7JmVRoVBw3/Ptdln/1Hgveepanr53Fydc+BKcPaPexm6OlfhA1wJ76l0Dvutf1LUdb73xxlOlWUhtOGyz9l5JZ9vWi+5XHpeS7jT4qd90vssv0kOwud11mRAs/VGsh+ITA5Ds6pUBHow0U7oQVzyvl1dbSWvM2qSB26mmq+tZbpEfNhCfeAtFDDl8sJXfffTfPPvssl19+OW+88QYGgwGPR/LSr3sorrYT6mdqctd2l5t3V+4nxNfIeaPiW73I5pbXMrVvBLOGeSE+eBRwuDx8tiGLDftUu9QjrdvYX2JlwdY8fE0G5gyPJdi36c8NIH9fOh8+dQ/lJUXsTNtFQmTbW5S2qx8E0PHm6Hgi7XtVgOStXHLmUshpbNyEujv3j1YG4ygYDYkkv0L1/I0KNDM4Lqihqfth+EUqnf/1b8P46zVxv67CaYON7ytj3ZpxsFVCxq9KlylmWNuOU5kDKSc06ZJyu93ceOONvPbaa9x00028+OKLDRLW23Ir2FdiJaGF2cO6fWXUOt3M7hvbqnGodboxG3Sc2L/7CPeZDDrmjUnEz2Tgt91FxAZZWp0FtUSvMD/OGRHPN5tz+HR9NnNGxBIZ0PTMMDoplb++MJ+09D34+XW8m6mlhkH7O/xoxwtFu5S8srcqrSUZyjjEjoSE8VCdr3LTq/KhNAMKttatWGc0AqL/NBwdZDTsLjdp+VXkVdjwNerJKbcRG+xDmF8LgfWAWCjcrqRDhl/Uc9xoxxLpP6v4Q2sxLikh/Sf1f9/T2/a3qi1XRZMDzz5sO6fTyRVXXMFHH33E/fffzxNPPNFwkXe4PHyzKYdQP1OzF/6KWiebDpTTPzqA6MDWJT6Kq+zMHZvQ7YLYep3grOGx+Jn1LNyaR1Sg5bAYQluIDrJw/qgEvt6Uwxcbcpg1LIb4kKbrVAxGE2FxnRMP7F6fck/AXq2KynxCvRPis1crkTS/COh9okpPtASq3HRQP2h71cFGoyQD8pszGjFqX14aDYmkqMrOttxK3G4PQRYDCIEPkm05FUzqE968q0kIZQT3LlXH7HuqV8fU6CAqclQ8wd+LmofC7epmo/fJ6mLvLR6X6ix3wj2HyXbbbDbmzZvHN998wz//+U/uv//+g5av2VtCRY2TuGYubAC/7ylGCJjYu/W01vIaBzHBFkYltt2NcjQQQnDygCgCLUbmr8si3N98RCKDoX4mzh8Vzzebcvl6Yy6nDY6mT+TRdedqBqIjkVKpZ9qrVUqoN+vv+l5J9Q64qOncdSH+bAR0kNGorDMaBa0YjZi6mcbhRsPp9rCroIrsshp8jQZ8GskQmwx6qmxOMous9I1qQVRQ6CAwQQnD+UVAnJYBfVTweGDzx2Dwaf1mwFGtOg0GxkNc8/pHTVKZA/3PhLCD22NarVbmzJnDokWLeOmll7j55psPWl5lc7Jwax7hLaS15pTXkl5Yzbjk0IMksJvCIyVVdheXT0w6IvdNZyOEYFxKGL4mPe+u2k+Qx9hi1XhrBFhUXObbzbks3JrH9P6RDIkL6sARt4xmIDqSvM2wf6Xy8XpD9lqV6ZR6WtuygYQAS5B6hPdT7x1mNPKgZM8hRiOiYaZRYQhlc7kJm1tHoMXYpAvAz2RgX7GV6CALgS39gPVGpS+17g0t/fVokbUaitNb/65JqdxQbif0a6NryVqsZoh9Tzvo7fLycs4880xWr17NO++8w+WXX37YpkvSCnF5ZLNuFikly3YX4W82MKpX6zOCoio7IxKCST5GtJmGxAdzwzQDbyzPxOXxtBhobg2LUc/ZI+JYuDWPX9MKqXW4GZMUclRkSJo1EEKIrTSdEQxAd8xi6lJqy+CP91RMwBshvso82PubmhW0NWDYFK0aDfWQxemI/C0EAZMR1BpDsJrCqTZGYDWFYzWFIoX6Wujq8ta351YwLjkMXUtfSKOPki9f9V844V6l4aTROdSWq7iPf1TrF/ziXVC8W7UP9W1DdbLLDi4bjLoCDH9e3IqKipgxYwbbtm1j/vz5nHfeeYdtWlhlY3l6MZGBzc8e0vKrKKyyM2NgFMZWZgROtwePlJwxJOaY0mbqHeHPzSem8r/fMiix2luO57WCUa9j5tBYFu8sYFVmCTUOFyf07XytqpZmEDPr/r+p7v/36/6/GKjptBEdi3g8qn2ox6UarrSGyw47vwGTf9sDhm3hEKNRXutga3Y50lZJpCjD31GCv7OYkNoDRFp3q1NBYDVFkB42HbshEB+jjspaF9llNSSGtnJuliCV/rr6vzDlLi39tTOQErZ/pb5rrfUSd9ao2YN/NCSMbdsxqvJg6FwIimt4u7KykmnTppGZmck333zD6acf3iJUSsnCrXkYdLpmY1dOt4ffM4qJCjTTL7r1niiFVXZmDIoirAUNp+5KXLAPt5yYymvLMyissjWbjeQNep3glIFR+Jj0/HGgnFqnm1MHRneqHHqrWUxCiElSykmNFt0nhPgdeKzTRnWskbVa5ZcHJ3m3/p6fwVYBwy5q/UfeAbg9kr0lVjIKqzEZdFj8QigjhDLfup7DUmJyV+PvKMbPUUyUdScDi35gW+QsnHpf/Mx6dhVUE+5vxtfUileyPv113Zsw4UYt/bWjKd6t6mW8SZ/es0jNAobOa5u8fHWhEu9L+bPTm5SS66+/nrS0NH7++WdOOqnprnP7SmrYkl1BXHDz3+v1+8uw2t2cMbj1GYHV7iLAbGBK6rE7I40IMHPz9FTeXJFJXkUt0UdYdT0lNQJfk4EVe4qxOXM5c0jn9Zbx5lvjJ4SY3GiAE4FjwxF4NKguVMFC/xjvZgIF21Qzl14TlQBeJ1Nlc7JmbwkZhdX4mw1YmhJKEwKHIYBS32SygseQFj4Do9tK/6Kf0Hkc6HU6BMotIL0RFQ6IhaKdsGW+pv7akbjs3tc8FKdD4Q5InKjkUbzFWQt4YOSlB2Xhvffee3z88cc8+uijzRoHj0fyzcYc/E2GZt2RVTYnf+wvo2+kP7EtGBFQRqnE6mDOiNimv7fHEEE+Rq4/oTcp4f7klNcesYbSqF4hnDwgkqyyGr7cmI2tnbIjreGNgbga+K8QYp8QYi/wCnBVp4zmWMPtUnEHofeuTWNtmZryB8ZDr0mtr38EeKRkb7GVVZkl2JweAn2M6Lycilabo9gddjJ+zhL6FS9CSDd+Jj2FVXYKK+2t76Ah/fU3db4aHUP6L2AtUa68lnDZVM2DXwQkTvB+/1KqeNXwiw9Kmti9ezc33XQTU6dO5YEHHmh286055ewvrSHYt/lZ48qMEiQwqU/rSRmlVgfJ4X4MiQv2/hy6Mb4mA1dOTmJYfDBZZbVNKre2hUGxQcwcEkNxtYMlaUXUOFwdNNI/8aYfxAYp5TBgKDBcSjlcStl8H7zjicyl6k7Nz4s7NI9b6TIJAQNmdWpHOavDxbp9pewuqMLXZGhXLna5TyIZIVMItufQu/Q3AHyNenbkVeJoQm3yMBrSX7+EbO3rcsQUpsHO77zr85Dxq+oA2O9M72px6qnKhbjRB8Ur7HY78+bNw2w28+GHH6LXN70/u8vNN5tyCWuhKC6/wkZafhUjEoIJ9GnZ9ej2SGqdbs4eEef1jc2xgNmg5+LxvZjcJ4zs8pomlVvbQkqEP2ePiKNXuG/r7t920OpVSggRJYR4E5gvpawQQgwUQlzd4SM51qjIhu1fQqCXQnz7lqvAX9/TW78DbCdSSrLKali5p4Rqu4tAi+GIAlhF/v04EDSaiJoMepWvwWjQ4XJ72FNY5d0O6tNf17+pZKg12kdhGqx8Scl4txbTKd0L+VsgYZx3/ajrsVertqPD5h70fX7ggQfYuHEjb731FvHxzdf2rM4opaLW2WyFs5SSZelF+Jr0jElqXZussNLGuJQwEkKPoe6LXqLXCc4ZGc9pg2PIrbDhcB2ZkYgL9mFgjBfdA9uBN7ex7wA/AfXKWLs53luOuhxKg8joqySwW6Nsnwpkxwxve+cuL6l1uvnjQBk7civxMenU3UQHZEflBAwnz38gsdVbiancgp/ZQFZpLaU1TffIPYyG9NdXoLroiMdz3FG4UxkHS3DrWWEuu6qs9gmFpMktr9sYjxtqilVKa6Oblx9++IHnnnuOm266ibPOOqvZzSttTn7cltdiB7j0wmryKmxM6B2GydDyZcfucqPTCWYMOga7L3qJEIJTB0Zx3qh48ittnRZDOFK8MRDhUspPgXrpbxfQPc/maLHrB5Wp401euaMG0r5T6/ZuOrh3JEgkeRW1/L6nmLIaZ92soQPdV0KwL3gCxT7JJFWsIaJmD2aDju05Fbg8Xt75WIKUGujq/6o7VQ3vKNwJK1/2zjiAivnYK6HfGa13lGtMVS6kTIWoQQ1v5efnc/nllzNkyBCefvrpFjf/Na0Qt2y+KM7l9rBiTzER/mav7nSLquycOTSGoFbcUMc6Qggm9QnniolJlFgdWO0dH0M4Ury5kliFEGHUFc0JIcYDFd7sXAhxmhBilxBijxDiviaW9xdCrBJC2IUQdzV6P0EIsUQIsVMIsV0I8Vcvz6fzKcmA3QtVn+bWkFLpLDltMGB2h6d82lxuNmeVszm7HLNBp6b3nVFTIXTsCZtOhTmG3qW/EenKpcbhZn9JG8ph/CLAWqTSX93Ojh9jT6OtxqH8AOT+oWII3si81FNbrmYcg85p+O54PB4uu+wyqqur+eSTT/DxaT7bqLDSxor04hZnD39klVNlczElNbzlYkvUbCTU38S45J7fcrSeYQnBXH9CCtV2FxW13eu34Y2BuBP4FuhdV//wHnBraxsJIfTAf4HTgYHAhUKIgYesVlq3r2cOed8F3CmlHACMB25qYtujj6NG+dMtwd7doeVsqBNIm6588R2ERFJYZWPlnmKKqhwEWYydrk8jhZ5d4adSawyhX8kiokUZGUXVVLflricgTqVebvlES39tibYaB7dTuZYswZA81fvjeFxgr4DRVx5Uj/Pss8/yyy+/8PzzzzNwYPM/Oykl32/Jw6hvvneI1e5i/b5Sekf4tRpPkFJSUePk3JHxrbqhehp9IgO4aXof3B4PpVYvMgWPEl5lMQEnABOB64BBUsrNXux7LLBHSpkppXQAnwAHOTKllIVSynWA85D38+ozpaSUVcBOwItb9k5EShWUri1XP8TWqC6AzCUQ2hti2yiQ1gIOt4ftOZX8caAcg06nhMA6cNbgkVDgaDpl160zsTPidJw6CwNLfsLfU8mO3Ao83l7shVAFXnuXw+6fOmzMPYq2GgdQCRC1ZSoBwpuYWD2VOcod1UiIb926dTzwwAOce+65XHPNNS1unllsZVtOBREtVDivyizB7ZFM9iKttbjazoDYQPq1JA7Zg0kI9eWWk1IxG1RKeXfAmyymDOAvUsrtUsptUkqnEOJ7L/YdB2Q1ep1NOy7yQogkYASwppnl1woh1gsh1hcVdWIQtGAb7F3mnWvJ7YAd36i7sv5ndtgFvMRqZ+WeYvIqagm0GDB28F1WnsOHx7JHcNv+CaytbvoH7dT7sjNCSSyMLP8Fa3UFueW13h9E6FQdyPYvIXtDRwy751BvHHyCvTcOlTmQvQ5iRkCIlyKR0KQQX2VlJRdeeCExMTG8/vrrLVb7uuuL4iyGZtcrrFLNp4YlBLcqVudye3C4PMwe1nrToJ5MZICFm0/sQ5i/ibyKIy+oO1K8ucI4gelCiLeFEPV/ZW8u9E39ldt0tkIIf+AL4DYpZWVT60gpX5NSjpZSjo6I6KRyfFsFbHgXfCO8q1/YsxhqS6H/TJXpdIS4PB525leybl8pQoB/M+qr7cUj4cfyOO47MIZshx/RxhreKuxLtbtpN5rNGKyqrT21jKlcRHpuSduyMPRGpQ+0/i0V09E42DiYvDQOHhfsWqiyxFKmeX+sZoT4brrpJvbu3ctHH31ESEjLCqtbssvJLqsluJlAspSS5buLsRh1jPMirbWgys4J/SKJ8qJpUE8n2NfEDSf0pleYH7nlti41Et4YiBop5VyUm2e5EKIX3l3os4HGWhLxQK63AxNCGFHG4UMp5ZfebtfhSAmb56sflDd3dUVpkL8ZEiZASFLbjmWvhIoDSnnVXg3SQ3mNg5UZJWSV1hBoMWI6gi5VTVHgtPB4zgjeLerLQJ9ynuq1llujt1PlNvJ+UZ9mt6s2R7I77GR8XWUMq/iVXXll3slw1FOf/rr6FSVXcjxTsEOlsrbFOICSlq8pUbMAg5dCdvVCfEPOPUiI7/333+eDDz7gkUceYfLkllNkbU43327KbbFTXGaxlezyWsanhGFuRSaj1unG0s3aiHY1fmYDV09OZlBcINnlR1513V68MRACQEr5FPAAqibCmzSJdUCqECK5buYxDxXsbv2A6lv3JrBTSvmcN9t0GlnrVDvQAC8apNsqVApsQGzb8tBB3dXVlisFzV4TcBksHNi/h127dxJgyyNGlGJ2VyNkx6TCeST8XB7LvfvHsN/uz/VRO7kndguhBgfJlmpmhxxgWVUMm6zN3/2V+ySQETqVUGcekTmLKKq0tW0QliB1wVr9iuqadzxSsANWvaz0ldpiHKry4cAqiBoCoSneb1cvxJf8pxBfeno6N954I1OnTuXBBx9sdRerM0uotDVfFOfyeFieXkyon4khsa0XhRZX2Zk9PK7btRHtaixGPZeO78X4lDCyy2u9TyvvQLz5izxc/0RKuVgIMQM4vEPIIUgpXUKIm1EGRQ+8JaXcLoS4vm75q0KIaGA9EAh4hBC3oTKehgKXAluFEJvqdvmAlHKh12fWEVhLYPOHyh3SmktHepSUBlKltLZF4kBKlYs+/CJImUZOeS0f5u+nyLealLAaAt1lBNryCLHtJ9CWh67eSEhw6i049L54hNHrWEeR08L/CvqzvTaEob4lXBu5izDjwUGxs0P3s7Y6gjcL+/FUr7X46Jp2IRX79cXkrqVXxVqyd/2Cc8SZGNsyy/GL+FP9dfyNB7k8ejztNQ4et3ItmXzbVlvjrAUhYeRlDd9Ph8PBhRdeiNFo5IMPPjhMSkNKic3pocrupNrmotLm4sdt+S2mtW7JqqCi1smc4bGtymSU1TiIDfbxqmnQ8YhBr+O8kfH4mw38sj2fmGCfVvtndOjxm1sghOgvpUwDcoQQIw9Z7E2QmroL+sJD3nu10fN8mp6NrKDpGMbRw+NWypkS7yS5962oa884q209fwEqsyFhHK7EKfy2s4AftuXja9QTGxqIjUBsRFPoP0CtKz1YXJX4OYoJsBcSbNtPcG02Jk8JEgFIXMKMU++LS2c+yGhICYsrY/mwqDdCwDWRaUwPzGvSrph0Hq6LSuPR7JF8XJzCVZHpzQ4/N2AoRncN8dXbKNwVSOSgaW07/4BY5Zrb/AmMuAQ6stCvu9Je4wCqKt9aqGoXvBGJhD+F+Eb/BfzC6t6S3HPf/WzYsIH/vv0hJQSQuauQ4mo7xdUOymoclNc4cXkkQqgfpJQSH5Oh2aK4GoeLNftK6RXmS6+wlkWfPVJSbXdxxcSkTu1pcKyj0wlOHxyNv9nAVxuziQywHDV125ZmEHcC1wDPNrFMAid2yoi6C3uXqcChN+1Dyw/UTfcHH1SN6hU1JeAXjhw2j++25vPbriJigy3N3yUIHTZjMDZjMCV+fYCJdf0crPg6S/CzFxFSe4BgWzYB9nxlIKSkwBPIC6Xj2FIbwWCfUq6NSiPC2HIqXV+fSmYEZ/NjeQITAgoZ4NNMfaQQ7A8ej9FdQ2TxamoOBOObONz7z6Be/XX/CiVN3e+01rc5ljkS42Atgv2/Q8SAP3uUN4GUEofbox4uD7Iil5LAfmzOj6Y0M4NSq4ONK5fy+X+eY/ip57PXfxCZq/eDBJNB1/AI9ze36eK9JrMUp9vDVC/6N6g2oiHHTBvRrkQIwdS+EfibDXy4Zj8hvqaj4pJrqWHQNXX/T+/0UXQ3KvNg6+dKObM1t42zVklp+ARDn1PadhynTQW/J9/O1gIny3YXEh/i2/a7KSFwGPxxGPwp9+lFTvBoAPQeO772YnbmlrMgR4+UkhtD1nG6/x6EG6RH4NT74tD7NLQZPZS5YZn8UR3OawX9+XfiOky6ZvygQpARNg19YS0he3/C7eePPqz5IPfh2zdKf/UNh4TRbfsMjhUKtqu2rD6h3nUfbIz0wK6FSIMZe/KJOGxOHC5lBGodbmocbmqdLmodHuwu5RIUCMzuGnTSyY/6CcicCkx6HfbKEn585W9E90pl3i0PYTIfefZQSbWdrTkVDIkPItSvZVfhn21Eo4/rtNa2MrJXCH5mPW+u2IvbI1tVxT1SWnIxndPShl2aWdSZuJ2w4W2lbNlaZoiUqoLVYYURl3mfSQLqx16VB2OuolgfzsfrdrX5bq01yh06vtgF+0tMxAf7cPLAKEyWVFY5yvBzlhBoyyXYdoAgW26juIbEqffBbghECh0WnYdronbxRM5wPitN5uLw5tNSpdCzJ+JU+ud/h//2r1VMJdCL4H499emvG94C35CDCrh6BEdgHOwuN0U7lhNflcfWgKnk77UisAJ1RelCqYTqdQKDEPibVX2CkB4C7ZWsj72UMH9V2+LxePj4+Qex1VRzw7/f6RDjALA8vRiTQcd4L2QyCittzBgcfUy2Ee1q+kUHcvN01ca01Opo1RgfCS3NUWa1sEwCPdNA7P4Zyg54V3SUt0m1gEyZ3jZpZVBy4SlTccSM4YOlGQghOkzPXUrJzvwqfttdhMcjmdY3gqHxQQghkIDVHIHVHEGhf52ybKO4hr+9kLDaDCKq06m0xCCFnsG+ZUwPzGVBWQLj/QvpbWk+40hVW5/GkMJv8dn6GWLEJd6JGtZj9AFzoLqQTruvbd3QujNHMnMAsnNySCpbR4klEWtgKgFexmn8HQUcCBpHiX9qw3u/ffE2uzas4LxbHyUmuXk3VVvYV2xlf2kNU1LDW+0/YrW78Lcc221Eu5rEMF9uOTGV//2WQVEnVl235GK6stOO2l0p3avcRd7c9VqLIGMxhCRDfBsawtdvGxQPQ87nx235ZJXVEOkj+On9l9EbjMT27k9sSn+CwiLbPP222l0sTitkb7GV2CALpwyMarWK9dC4xn45gdTiX0gpW0GFORaEjkvC97DJGsb/Cvrzz8T1GETzedkeox+bgmcwpnwBhi3zESMuVTUP3mIJUp/Rqv/CuOshsPN67h4VjtA4VFmrCcv+GSn07A2d7HUQ3+yqpNYQzO7wPzOdDuzawoK3nmPIpFOYcOa8No+lKdwe1esh2MfIsPjgFtetbyN6xcRex3wb0a4mKtDCLSem8vryTLLK2iCc2Qa8umUVQpwJDAIa5qJSysc6ZURdhdOmXEvmwNZVV91OldKqN7VdSsNRoypgx17D9gIbS3YVERtk4ZOn7uaPJQcnh/kFBhOT3I/YlP7Epqj/o3r1wWg6fFoupWRXQRVLdxXh8kimpoYzPCG4ff5dIUgPPxmDx0FixVoqzLH46uHqyF08kzeUr0t7cV7YvpZ34RvCOvfJjK/6CbH1M+VuMrTBleEXofpHLH4M4kZB31OVjtOxxhEaB+msQbdlPn6uEtLDTsKp924fQrowu6pZk3A1br363G011bz/rzsJDAnngtv/0WG+/205FZTVOJk1NKZVF2mp1UHviJ7TRrSrCfEzceP0Pvy4LQ+fTjC4rRoIIcSrgC8wHXgDOA9Y2+Ej6Wp2flN3Z+/FRShziVp3yAVtz1+3FsD4GykhiA/X7Cbc38SKr97hjyXfc/oVtzF59sXkZu4id+8ucjPTyM1MY9XC+TjtqghNp9MTkZD8p9FI7kdwfCrriySZxTXE1M0aQlqbNbSG0JEWeTp66SS2chOV5lhG+Zcw0b+Ar0t7Mc6/iASztYXtBW6/KDaL6QyvWITY9iUMvaBtfQr8I1SspmCrKlaMHgL9TleztmMhsFmwvS5bKaxdxgFHNe6NH2NxlLE77GTKfJO83jTAXsie0BOo8PlTzODLlx+jND+bG596D7/A4LaPpwlsTjerM0uID/FpNRupvo3oWcN7VhvRrsbfbOC8UQmtr9gOvPm1TpRSDhVCbJFS/l0I8Sw9Lf5gLYaMpd4J8RXvVrr78WPaVsEqpap3SJ2BM2ooHy5Vwd7sbWv57o2nGTp5BidfeD1CCHoPHUPvoWMaNvW43RTn7q8zGLvIzdzF3m0b2NhoxqHzCSSiV19MAwexp3QAMSn9iE7sg8HUfkMhhZ4dkbPQexxEVe+k0hzD5RHpbK0N4X8F/fl7wh/oW3A1GfQ6io2xZEdOI6HgV9VTeeBZbevHLXQqcC0lFO+B/KdVemf/MyE8tfsaiiM1DrYK5JZPEPYqtoWeSo2v9xcAH2cpleZo9oVOaXhv/aKvWb/oG0695OaDvltHytq9pdhcKq21tRlJUZWN8T20jWhPxRsDUS/VWSOEiAVKgOTOG1IX4HaoytLWqp/tlaqC1T/qIKkCr6guVNLfA2fz0/Z89pdYsdhKeO+J24lK6M2Fd/+r2R+YTq8nMiGFyIQUhp9wBqAKkn7emMmuHdswV2YRas+n5MBuVn7/CS6HvW47A5EJKcQm922Ia8Qm9yMgtPUfcz0enYGt0eegz5tPWE0GmKK5IiKdl/IH8UN5PDNDslrc3s9sYLstkdD4qfhlL4M9i1Q6cFsv7EKogLWUKsC/4jmldTVgtpKO6E7FdfnblHxIe41DbRls/hjptLE+8FQ8ft6LIOs8DgweB1ujz8FTN1srytnPFy/9nZTBoznl4hvaPp5mKKtxsDm7nMGxgS1WVsOfbURP7cFtRHsi3hiI74UQwcDTwB+oDKY3OnNQ3RLpUXfAHjcMOKttrhJ7lbqAjbmaHQU1LN5ZQLhF8sr9NyOlhysf/S9mH+8vJOmFVSxJK8Lh0nHi9GmMSgxpmLJ73G6KcvaRm7mLvDo3VebW9QfFN/yDQolN6U9MSn18oz9RiSkYjE3PNjw6I5ujz2Nk7scE2XKY4C9Y6RfJpyXJjPIrJsbUvNy3EAKLUcd6dx+mxNWiy1mn3HK9Jnp9vofsEPzCQYYpKZSVL6j6iQGzlAuqLRInnUH+NtVatb3GwVoEWz5BejysC5qB0zcSr89ISgIcReyIOBOrWWV/uZwOPvjXHegMRi6+72n0+o4rrlqRXoxBp2N8SutZakVVduaMiOvxbUR7Gq1+W6SUj9c9/aKuD4RFSulVy9EexYHVUJGlGqz4ti5f3IDbqaqlJ99GqfTnwzW7CPU18fmz95C3dxd/+cdrRMR5p+Nf63SzdFchuwuqiQwwc8rAKMIPySPX6fVEJfYmKrE3I6ad0fC+tbK8wWAoN1Uav3/7IS6no247A1GJKcSm9Cex/zAmnjkXveHPH7Nbb2Fj7DxG57yPv6OQqyJ3c/f+sbxW2J+/xW2kJZey2aCnyuYkI3QMqZE1sG+ZunjGDPP+czwUIdTfwSdEdUVb86oKbA88S+23g9u7esWRGoeqPNgyH3QGMuPOorLGgn8bZkZ+jmKKfVPIDv6zQdUP7zxP1u5tXPHwS4REtqEmpRUOlNaQWWxlYu+wVit6K2uPvzaiPQVvgtR64EwgqX59IQRdrrJ6NKnIUV27Igcq9UxvkVLpMw08C2dYPz76LQOPhLXfvsum3xZy5lV3MmCMdy0iM4qq+TWtEJvTzYSUMEb1CmlTUZ1fYDB9ho2jz7BxDe+53S6Kc/aTm5HWEBTfs3kNGxZ/y57Na7jsgecOMhIuvQ9/xF7E6Ox3iXflckl4Bq8V9mdxRSynBLes5O5nMrC3pIao5FMIdNbA7h/VRbQt1dZNIYTqvmYJVi7AdW8oozFgtsp+Olrif0fqVqrIgq2fgdGH6n7nsSfHTYDF+7t9g9sGQrA96iykUHOOXetXsOSzN5k4cx5DJ5/a9jE1g0dKlqcXEWgxMCIhuMV1pZRU1Dq59oSU466NaE/Am2/gd4AN2Aocfb3ZrsZlUxlOlkBIPbVtvvOqfKXP1HcGv2wvYG+xleqMDSx461mGTT2NE+e23NIRVJbI0t1F7MqvIsLfzJzhca36e71Frzf8OduYfiYAHo9k4Sdv8Os7z/D+P+/k0geePchIOAz+bIi/lDFZb3OGZTsrfSL5qLg3I/xKCG9B20mnExj0OrbnVzNuwFnotnwCO76GoRce1JfgiDAHqoejGv54T8l29J8JCeO8F7VrD/nbVCqrbzuNQ+le2P4FmAORw+ayM8+NQefxPg1VSvycxWyKPh+bMRiAqrJiPnr6XqJ7pTL7uvvbPqYW2JFbSXG1gzMGR7faC7242s7A47iN6LGONyY9Xkp5jpTyESnl3+sfnT6y7oCUqneyvUrdkbYlj99WASYfGHUZaQVWftlRgNFayPv/upPoXn2Yd+c/W70AZBZX8/7q/aQXVDEuOZS5YxI6zDg0RZXNSXZ5LVdedzPTLruTLSt+4oMn78LtOqhlOHZDIBviL8OjN3JL6DokgjcK+9Fa4ysfo47KWhfZVW4Ycr4qntv2mcoi60hM/qpmwmCBLZ/ATw/A7l+UJEpHk7/1yIxD8W7Y9rmqkxh+MUUOMyVWB76tVCM3xt9RRF7AEAoCBgN1UhpP30ettYpLH3iuw6Q0QAWbV2aUEBtkoU9kyyneLrcHh1sy6zhvI3os442B+EEI0XHz02OJgq1QtBOSpniXAluP2wG2chh7LeUeHz5cs58AvYv3Hr8FgCsfaTkobXe6+XlHPt9tzsPHpGfumATGp4R1miSy2yPJKa9FSrj+hN5cPC6Rh++/m6mX3sHmZT/y4b/vxu0+uFFRrTGE9fGXE2ms4ZKQ7WyuCWN5VVTLBxICP7OeXQXV1GCGIXNVGuvWTzunYZDRV9W1mANgx1fKUKQtBFuT3WvbTv5WWPVK+41DwXbY/pXKiht2EW6DL2l5VfgYdV7PVI1uKy6dmbSI0xu2Wfblu6StX85Z193XYVIaUkrsLjdr9pZS63QzpW/rmXAFVXam9YvQ2ogew3jjYloNfCWE0KH6U9fJwsvATh1ZV1NTAum/qDvRxPHeb1cfdxg6D1dwMh8ty8Tp8vDTyw+Tv38P1/zjNcJjmy/G21dsZXFaIVaHi7FJoYxNDu1UrfxKm5PyGgeT+oRzxpCYBj2oKakRnHPZ9Rh1sPjd5wBxWBaM1RTO+vjLmHngbVZY43mvKJWhvmUEGxzNHk+v0yHwkJZfxYjEYMSQC2DzR8pIDL+4bbM0bzFYlJy4y66kVHb/AL1PVH2cfdrZqKbBOISrxj1tJXcTpP+oDNjgc8FgJqvEis3pJsDLTB8hPfg6y1kfeykOg7qbz9q9jQVvPcvgiSczceaFh21TLwNuc3qwOd11Dw82lxu700Ot043d6cbmOnx5/Qyxf3QA0a1c9GsdbnyMOqb36yFaWscp3hiIZ4EJwFbZld2zjyYel4o76AyqAVBbCrsqcyBuNKRMY/HOQvYUVrP75w/Y9NsPzPzLXfQfPaXJzdweyZJdhWzPrSTUz8TMoQmdeufl9kjyK20EWgzcNL0PfSIP9hEb9DrmjU0gu/wyTHodP7z1DEIILrr3qYOMRLU5mo0Jl/JX15f8Ne9U3i5K5faY7S0e28+kp7DKTmGlnajAaNX4ZuunsO0L1XK1LSnEbcFgVhdkt1MZ//TFkDwF+pykUme95UiNQ9ZayPy1ri5mDuiN2Fxu9hRW42v2zrXkkaC3lbHBZxLbPfHYS6xUVlbx+WN/xRQQQsJZd/DzjoLDDEDjC31TmPQ6zEYdFqMei1FHuNnc8Nxi1ONr1LfqWgIVe5g3NlFrI3qM481fLx3YdtwYB4DM31Rh26Bz2yYyV1OqZKpHXEx6kZUft+VTlb6OhW8/x/ATzmD6+X9pdtMVe4rZnlvJqF4hjE8ObTX4dyRU1jqpsDmZkhrBaYOim1XfjAq0MHNILG73xeiQLHjrWYQQXHjPvw8yEhWWeHS9ZnOhdT3vVQxjTVUh4wKKmh+AEPga9WzPqyTEz4QpJEkFk3d+qx4D57TNKLcVvVHNKDwulZ229zdInACpp7SuynskxkFK1exn/woI73dQa9rMomqklOibSWv9srQXG61hVLuNVLsN1HiMeBqaLqpixeIF/8FalEvUvCfYUwEWY23dxV1PoMWAudGF3mLUYzE0em7UYTboO2S2qrUR7Tl4YyDygKVCiB+AhjSVHpvmWpIBOesgdqSScvAWlw1cNTDpVipcJt5btQtZkcdH/76L6OS+zL3jiWZ9trvyq9iUVc7w+GAm92nDnWwbcXsk+RU2QvyM3DS9D70jWr8TnJwazpaccnxmqjbkC95SDQYvuucpdI36F5f5JpHc10nvTaW8XZTKIN8y/PWuJvcJYDTosNuc7CmsYmBMkEohdliVQm76L23PGGsPOoOKLXnc6q5+/0qIHw19Zyi13UM5UuOQuQSy16rMtn5nNBjBSpuTrNLaZtNat9UE81lJCinmSpLNVQToHIRSTmnoSKRvKBajjsxVP7Jw22KmX3gDZ1x+dpe18PRIiVVrI9pj8MZA7K17mOoePRdHNexaoAquereho6r0qC50o67AFRDHx8v3Yq2uZv4//4oQOq565L+YfZq+oBRX21m0s4CYIAuTUzvPOFTUOqmsdTK9fySnDIzyWmpZrxPMHZ3AMz/vYvJ5VyOlh4Vv/weh03HhXU8eZCTKA1M5re8OXt2u5/3CFG6I2d3ivv3MBrJKa4kO8iHU16T0rRzVkLUGzP7Qa9IRnbPX6PRK4l16IG8LZK9TxXZ9T1NyHkIcuXFI/0n1D4kdeZDUiESyK78Kg140eQPhkfBeUSoRhloeid+ISech0JbLntATyAxPAqA49wCL3/wnyYNGcsZlN3fphbmoys5wrY1oj6FFA1FXJJcqpbzkKI2n65AS0r5XGUgDLmqbH7wyB5ImQ68JLNlZQFp+JUv+9wgFWRlc9883CYtpWmjN7nLz/ZY8TAYdZw5pXSq5PbjcHvIrbYT7m7n1pFSS2vHDjQy0MHtYLF/+kcNJ865DSskP7zyPEIJ5d/7rICNhiB7I1NIdLMmLY2JVAcMCmi+6F0JgNujYnlvBhJQwDDodJE9TM4l9y1XzoJgRR0+QT+iUi0l6oGgX5G1WOk+xw2Hz/HYaBw+kLYDC7ZAwXml4NTqfoio7pVYHgc3MHn6tjCXL4c9t0dsw6Tz4OMqoNMc2CPG5nA7e/9cd6HR6Lr7vmQ6V0mgrWhvRnkeL3yYppVsIESGEMEkpm09N6QkU7oCyfZA6o20By5pidVEZegF7iqz8sDWfXT++x5YVPzPrmnvoO7JpzSEpJT9vL6DK5uSckfGdEswrr3FQbXdxYv9ITh4YhdnQfp2iCb3D2ZxdQU5ZDadcdAPS4+HH914EBPPu/OdBRmJg//5sL9vNm0X9eMpnDRZD8xcLs1FPZa2T/SU1yuUlBPQ9HZw1kP6zunNPGKcUXDszLtEYoVOpp1JC+QEoTFOzyrYaB49LxVSKd0PS1MP0p1weDzvzKtVsrokLao1bz6clyfT3KWesfxF6jwODdLA1ek6DEN8P775A1q6tXP7QC4RGdVDBYTvR2oj2PLy5Ku0DfhdCfAs0VBr1qBhE/lbI3aguQjHDvd/OWaNmHGOvpcKl571VeyjcsYqf3nuBEdNnMu28q5rddN3+MjKLrUxNDScu2OfIz6ER9bOGyEAzV09OITHsyOWV9TrBBaMTePqnXdhdbk695CYkkp/eewmh0zH39n80GAmDTse0QYl8tiGbzwvjuDQmq0H+oSn8zQYyiqqJCrTgbzYol8+gc9TfJXutqri2BEPCWCV1crR0loRQhqE93hK3U1Vyl+2F3ier2MYhZJfVYnN6mm08/2VpEtVuI5eFpyOQ+DmK2NlIiG/X+hUs+fQNxp8xl2FTT2vHIDsOq91FgMWotRHtYXhjIHLrHjqg59XL26tg4V3KndH3dO/dGR43VBfA2Otx+0Uxf0Um+Vl7+eLZe4lJ7sfcFjp27S+xsiqjhL5R/gxvRcumrZRZHdQ4XJw6KJpp/SKOaNZwKBEBZs4eEctn67OJD/FhxiU3g4Sf3n8JIeCC259AV5eFExPsy7CEYBZk9WNqxT56BQlkMzMAnU6g1wm251YwJikUnRDKxRc7QsUCineruET6z7BvhfLjx41Sf7PuiMuuKsQrstV3qglRQptTpbX6NZPWmufw4cfyeKYF5pFsriTQnkeh/4AGIb6qshI+evpeohJ7M+f6jpXSaCv1bUSvnJSktRHtYXij5vp3ACFEgHopqzt9VEcTvVllzJQf8P6CU9/8p8/JEDeCpWmFbMnM59un70Cv13PVo//FZGl6X5W1Tn7cnk+Yn4mTB0R1mK/W6fZQUGEjOtjCNVNTOq0py7jkMDZllZNVWktEgJkZlyrJ8p8/+C8guOD2fzQYiYm9w9lbZOWZiun8z/QFdp/IZt1EvkY95TVOcstriQ9pNHahU3GA8H5K0C5rjUoTzVoDMUMhbgz4BHfKubYLZy1sna/SpAecBZEDmlwto5W01g+Le2PSeZgbuocgex7ZQSPZGXEmUujrpDTupba6kuv+9Waz37WjRX0b0cGxQV06Do2Op1WnrhBisBBiI7AN2C6E2CCEGNT5QztKGEww5U7lc/YWa5FqezlwDpnFVhZszmHp649SmJ3JpQ8+T2h0EymSKNfPgq15eDxw5tAYjB1U61BSbaewys5pQ6L560l9O7Vjl67O1eSRErvTDcCMS2/hlItvZO1PX/DZCw/j8ShNR6Nex0kDoihymPg/+6kE2vNotkpLCHxNetLyq7DV7ffQ5QQnKg2n0Vcro5G7Edb+D3Z8o4QRuxpHtaoKry5SLrJmjENFrZPsstpm405ba0LYYI3g7OBMEt37yQidyvbIWQ1xh+VfvUfa+uXMvu4+YlP6d9rpeEN9G9E5I7Q2oj0Rb1xMrwF3SCmXAAghpgGvA+3s+HKM46hWYiNjrqbSJXh35T42f/8O21cuYva199F3xIRmN126u4jCKjuzhsYcec9o1Kwhv9JGfLAPc8cmdngsoznC/M2cPSKO+euyiA/xQQjBaZfdivR4WPTxqwgE5/317+h0OhJCfRkcG8iiXBibPJaJ9nVUmmOadOUZ9DpsLg9p+ZUMSwhG0MwFxy9CtRxNngrZ6yFvo9LMCu6lMoXqU1OPJrYKJQxor1ZGLCSpydWklOzKr2w2rdUtBe8X9SHKUMNc3/XsiJhJVvCYhvPJTt/O928+w6AJJzJp1kWdeUZeUVhlY0LvsINnfRo9Bm8MhF+9cQCQUi4VQhyfSc4el1IenXQrbp8w5v++l+2rl7Lko5cZeeIsTjj3imY33ZZTwfbcSsYkhZDiRYFaS9T7fB0uDzOHxjAlNaLDZiPeMiYplM3Z5ewtshIZaEEIwelX3AbAoo9fBSE479ZH0el0TE4NZ19JDa8WDiIpvorYmrRmjYS/SU9+pZ3YKjuRAa1IjZgDoPd06DVBpaRmr1PuHb9IlfkU0f/odJirKVXGwWWHofNalC8vrLJTWuMkqLm01ooYshz+PBi2jO2x51EQ8Odk3V5r5f1/3oF/cKhXasCdjd3lRicEpwzU2oj2VLy5qmQKIf4mhEiqezyEKpw7vpBSNQ7qfyZEDWLZ7kJ+X7+FBS/eT2zvAVxw2+PN/mDzK20s3VVEYqivV+0ZW8Lh8pBVVkuEv5k7Tu3Lif2jjrpxAOVqOm+Uqu+odwnVG4kT517L6oXz+eKlR/F4PJgNeqb3j6DE6uRN+4mU+KYQ4Cho2t0kBD5GPTtyK3G6vWw/YrAogzDuelWhLD1KlG/t/5TRcHdihra1CDZ9CB4nDGu5t4XLo2ZHPs2ktVrdBj4rSWKIuRBj6vSDjAPAl//9B8W5+7n4vmfwC+x6GYuiKjtnDo3R2oj2YLy5slwFRABfAl/VPb+yMwfVLanOh8j+0O8M9hZb+WJ1Ogv/cyd6o5ErH3m52UBhrcPNgi15+Jr1nDY4WmXotAMpJcXVdkqsds4aHsvNJ/YhJqhrg5OhfibOGRVPYZWdeqkuIQRnXnUHJ869hlUL5vPly48hpSQl3J9+UQGs3V/BooA5VJjj8Hc0rddkMuhwuD1kFLUxH0JngOihKkYx+DzV5Cljser0tve3ju8HUZWnjIMQMOziVnWcskprsDk9zXZW+7YommqPifEDelHml3LQsg2/fse6n7/k5ItuoM/QsR12Cu2lstZJuL9ZayPaw/Emi6kMuLU9OxdCnAa8AOiBN6SUTx6yvD/wNjASeFBK+UyjZW8BM4FCKeXg9hy/w7BVKiXQUVdQ5ZS8syKTRa8+QknOfq578q1mC5Q8UvLD9jxqnW7OHxWv7hzbgd3lprDSTlK4H/PGJBDZjfT1RyWGsDmrnD2F1Q3qs8pI3In0eFjy2ZsIITjn5oc5oW8EB0pr+HFXOeEj5jEm9338HEVYTYfnzvubDewrsRLia2q7qq0QqpVpWB9V5Z61Bg6sUnpL0UMgfmzb+oo3RXmWSmU1+qiueK1kUqm0Vmuzgekyay0LqnozNNqCb9jBSQ5/LPmeT569n6SBIzj1kpuObNwdQH0b0etO6K21Ee3hNGsghBBvA80puEop5dUt7bhOpuO/wClANrBOCPGtlHJHo9VKUcZnThO7eAd4GXivpeN0Om4H2Mpgyl14zEF8unIfiz/+P3at+ZWzrr+f1OHN94pYlVFCVmktJw+IbJd0d/2sweWRnDMyjgm9w7udAJpyNcXz7x93Uet0NxhBIQQz/3I3UkqWfv4WCME5N/2Naf0i+GFbPmtzHejjLmFM9jv4OkqoMR18JyqEwM9kYFNWOYNjA4lrbxA0ME5lFNWUqqK7/K1KEym8r3JLtaURVD2NWoQybJ76vxX2FFUD8vC/n5QEOAp4quIU9Ho9Y1NjGy2S/PLhK/z43ov0HjqWKx5+qUulNOqpbyPaN+rIYmka3Z+Wvm3fN/FeInAbakbQGmOBPVLKTAAhxCfAWUCDgZBSFgKFQogzD91YSrlMCJHkxXE6DymhMle5K8L7sHxXEd9+8w2/f/p/jD75LKaefXmzm2YUVbN+fxmDYwMZ1I788PpZQ59If84f3bmtRo+UYF8T54+K5/3V+0gI8W2IxQghmHXNPQAs/fwthNAx54YHSAn3Y1VmCSkRiWyIu4Qx2W/j4yij1nSwX92g1+Fnhq25lTjcHpLC/ZrPbGoN31Alvpc0BXI2QO4fqgAvMF4ZirA+3mU+Fe9WabW+Yap3hRed5CpqneSUNaHWKj0E2vNZwig2WsOZ3CesoWGTy+ng0//8jfWLvmb0yXO44PbHMRg7XytTSolHqtmvp/65RyLr3nNLqbURPY5o1kBIKb+ofy6ESAEeAKYCTwJverHvOOqF6hXZwLj2DbN5hBDXAtcCJCY236mtXVTlqGrePiexv8TKOwt/56f/PkR86iDO++tjzf5Aymoc/Ly9gKhAMyf0a5v0gJSSomoHHik5f3Q8Y5M7r9VoRzIiMZhNWcHsLqg6aLZUbySkx8NvX76DEHDyVffywZoDLNpZwHkj41kfdzljs9/C4izHZgw+aL96nY4As2BXfjVOtyQ10v/ILkwmP5Uemzi+Trl1rZoN+IYp11PUoOaFGgu2K0HHgBgYcgEYW58VSilJy6/EeEhaq5BuAu35HAgYzlv7hhDkA8MS1I2EtbKcdx6/hYzNazntsls55eIbEUJQ63Dj8nj+vIB71AX80Iu6QFBvR+uPKOreqncJSOpbQ6pXEoGUEiEEJr3AoNdh1Osw6gVmkx6T4c/XA2ICtTaixwmtqbkOAB4ERgBPA9dLKZsX+T9k8ybe6/CmQ1LK11C1GowePbrj9l9bBpYgGHEx1U4P//tlK989ewdGk1kFpZtpBO9weViwJQ+9TnDGkBilUOol9bOGvtEBnDcqnvBjSPRMCMG5I+P5949p1DhcDXfC9ctmX3cfUkqWffUuCMGU2TexOK2ILTkVDIsPZX3cZYzNfhuzqxK74WCXjU4nCPQxsLfYitPtYUBMYLuD/Q3oTUofKXYEFKVB9hrVinTfctURMHb4we1P61uEBieqRlIG7/42hVV2yg5Ja9V5nAQ4CsgMmczntjGU1hQzc6j6rhTnHuD1h66htCCHi+97hlEnzgKUhAoCogMtGPUCo16P0SCw1F24TQY9Jr3AZFDP9TqBoU7CxKBXzw06XcNrfd3r+mX1r3UCbWag0UBLMYjPgNHAM8DtgBsIrP/ySClLW9l3NtBY5zoepenU/XHZVUHcCffiMfrz6e+ZfPrMfZTlHeD6f79NSGRsk5tJKVm8s4BSq4M5I+IItHif/uf2SAoqbZwzIp5JfcKPyarUIF8j54+O572V+0kI1R98xywEZ11/PxLJ8q/eQwhBwtjL+H1PseodYIlsMBLSVd3QY7nx9gEWA9lltTjdHgbHBbXJ+DaLTq9mDZEDlZpv1hrYuxQOrFTCjfGjoXCnavbTqEWoN9Srtfo2SmvVe+z4O4rYGXE6u/3GsHrHfuJDfEgJ92Pv9g289ehNICU3/PsdUgYrgT+PlFQ7XNx6Yvvk2jU02ktLM4gxqDv+u4A7696r/8VLIKWpjRqxDkgVQiQDOcA8oOtLP1tDeqAqF0ZcCiG9WJlexGsvPE3Ght84+8aH6DOseS/ZpqxydhdWM7F3GIltkLuQUpJTXsuJ/aOYnBp+TN/BDU8IZktiOTtyq4gOOniWJYRgzvUPID2SZV++ywQ3yIHn82taIWcNi6XKEsP6uEsZk/0OUuhw6n0P2z7QYqCwys6mA+UMSwjuuBoQISA0WT2q8pXrKXudeiBV0V3/WW0qvDtQWoPD5SGgrk7A6K7Bx1nO5ujzyA8cyprdRdhdHqamRrBx6QI+eeZ+QiJj+cs/XiMirlfDfoqq7IxKDNGMg8ZRp6UYRNKR7FhK6RJC3Az8hApqvyWl3C6EuL5u+atCiGhgPRAIeIQQtwEDpZSVQoiPgWlAuBAiG3hESulN7OPIqMxR/YmTJpNVWsPTr33Ami/+x5hTzmbyWc33Tcopq2X5nmJ6R/gxuo29eAsqVVbI6YOP/UYrQgjmjIgnveBwV1P98rNvfBCQrPjmXQbWOtk3/ELS8qsYEBNIhU8Cf8Rdwqjs95AIXHqfQw9AgNlAWY2D9fvLGJkY3KGKtYCqZxgw+08pD6GHlBPa1I+i1ukmo8iKb11aq9lVidFtY0PcJZT49aHU6mBLdjkDYwPY9N1b/PDO86QMGcOVj7x0UBGcasIDpw+J6dhz1NDwgk7NmZNSLgQWHvLeq42e56NcT01te2Fnjq1JrMXgHwnD5mF1uHl6/mJ+fPkhEvoO5ry//r3Zi3e13cXCbXkE+Rg5ZWDbFFpLrXZC/IxcNDYRQxdURHcGQT5G5o5J4O3f92EJ0R8WL1BG4iGkx8Pv331ETK2LpYbLSQz1xc9soNQ3mY2x8xiZ+zEOjwObIfDgDCMh8LcYsdpdrNtXysjEkMMMUYdgCVaKve1gT2EVSNVHw8dZBgjWJlxJpUWl1S5PL0Iv3eR9+xwbF3/DqJNmM/f2JzCYDs5UKqyyM2NQFKF+Pbvbr0b3pGdckToClw3cdhh7LR69hfd/28EHj9+CyeLDFQ+/jNHUdFDS7ZEs3JqH0+1h5pCYNt3NWu0uXB64alJyp3SU60oGxwUxMjGEwkpbk8vri+cmzryQvOWfUvjrWyxJK2xYXuzfj3XxV1BjCiXQnoefvQghD1Z59TMbsLs8rN1bSpXd2ann0xYqap3kltvwM+vxcxTh1FlYk3BVg3HYV2IlI6eQqq//zsbF3zDjslu46J6nDjMOVrsLP5Nea8Kj0WX0rKvSkaA3wcjLITCW33cX8tR9N1NRmMMN/36HkMjmp/fL04vIq7BxehtbLTpcHkprHFw/NaVbVUZ3FEIIzhoRy66CKnWha8IA1hsJkKz8/hPWI+gbfS99o1QWU5lvEmt9riLAnk9i+Rpiq7aAhBpTKG6dupj6mgzYnG7W7i1lVK8Qgn269k67Ia1VJwhwFFJtimRj3EXYDarXltsjWbR6C0UfPISrsoCL7nmK0Sef1eR+Sq0OLh3fS2vCo9FleDWDEELohRCxQojE+kdnD+yo4hMKIy+D+NFkldZwz/0PsXfjcubc8AC9h45pdrO0/Eo2Z1cwIiGYvlHeN9vzeCT5lbWcMyKOvtGtV+EeqwRYjMwbk0CJVdV1NIVOp+Ocmx9h/BlzqVz9OZ+9+gy1jkaZ1EJQZYlhe/QcliXfxp7waZjc1QTaczG5qkFKLEY9eiFYt7eU4mr7UTq7psmvtFFmtRMlCyjzSWJ9/GUNxgFg8W/LSXvtVrBXccOT7zRpHADKapwkhqqufBoaXYU3DYNuAQqAX4AFdY+mqqyPXYwWSBxPjdPNvc++zuovX2fsjHNa1NsvqrKzeGchccE+TOoT7vWhpJTkVtQysXd4m7Y7VhkYG8iYpBAKmnE1gTIS5936KMNPPpeS3+fz5gtPNoj/NcZuCGRv6FSWJd/BlqhzcOnMBNrz8HWUYjYIzAY9G/aXkV9R25mn1Cwuj4fdeeXEykJyA4axMXbeQUH2NYu/56enbsTkG8DtL84nZcjhfapB3UBYHS7OHqk14dHoWrxxMf0V6CelLOnswXQlUkpe/GwJnz/3AIn9hnLuLY82G2y2O90s2JqH2ajj9MHRbap0LqyykxLhz1nD4475jCVvEEIwe3gcafnNu5pAGYlL7voHFbUOMn95n0/8zMy74a6mm+roTOQHDiU/YDAhtQfoVb6KyOpdeIQevSGETdkVDHRLEkJ92i/N0Q4OFFXgZ8slK+oE0sNPRgrlGpJSsviT/7Hw7f9gjh/I9U+8SmRM8x0Mi6rtjE4MoVeYltaq0bV4YyCygIrOHkhX88vGDJ6862osPr5c8fBLzQalpZT8tKOAKpuTc0fGtym4XFbjwNek59IJvY4rFUx/s4F5YxJ4fXkmPkZ9s3fFOp2O6x58kqcfcrLu6zfwMUhOv/QmzD7NXCiFjjLfJMp8k/B1lBBfsYGE8nUE4GRPdi1OdwQpEUeg39QGamutlOXtZW/0GWSHTWrIunI5HXz2wiOs+/lL/AZNY+pVfyOxBeNQn9Z6mpbWqtEN8ObqlgksFUIsABocvFLK5zptVEeZ0iobN/zlSqqK8rjx6XcJjmhe13/dvjL2FluZ1jeC2Da0+KxxuLC7PFx3UmqbKqx7CgNiAhmXHMYfB8pa7GNhMhi48t5/8uoTHpZ9/iZrF85nzKlnM2nWRUQmNF+bWWMKY3fEqWSGTiG6aivJJcspydkDtaEkx8ej64iq6+ZwWMnL3s+asLOwh//Zq6GmqoJ3HruFPZvXkHjyZZjHzmVyv5Yv/AWVNk4bHK2ltWp0C7wxEAfqHqa6R4/j739/hMyNKzj3lkca5A2aYl+JlVWZJfSPDmBovPcKrS63h+JqO1dNSj5qfaO7G0IIZg2LJS2/imqbC/9mWm4CxIb4cfK1f2PVqlMw7f6Zld9/wvKv36fviIlMmn0xA8dPa1b22qX3ITt4LDlBIwmpziA8dykiI41ekcHoA6I6vgWprZyqqgo+N89BHzW4IahXkpfF6w9dS0leFqfc8Bi7A0cyPiUMH1Pzx7faXQRYDFpaq0a3QTQVDDxWGT16tFy/fn2btikpKaFvv/7EDZvCFff+q9m4QEWtk4/XHiDAYuCC0QleSzxIKckqq+WMIdFa714gLa+S/y3LJD7Yp8UArEdKlqQVsi23kt7+LkzpS1i18BMqigsIiYxlwpnzGH/6+fgHt9z4x+PxYC3ezzTDdiaa9mAQgF/EwUJ87aWmBI/08I6cyV53ZMNd/74dG3nrkRvxeDxc/vBLrKqJRAKXjOvVbLyq/nty2YRejEjs+naiGscPQogNUsom74xbNRBCiCU0ocIqpTyxY4bXcbTHQACs37abj7ZWkBDe9KzA5fbw6YZsKmqdXDgmgWBf7ydSOeU1DI8P4aJxiVpGCupC+PmGbNbuLW3VRSelZPXeUtbuLSUl3I9TB4Sza+1Sfv/2Q9I3rUZvNDJ86ulMnn0Jif2HNmvc6zPH+gZ5uCw2D8u+X8FZo9R6zYdUaXtLVT6Y/diaeDlvb6klPsQHIQQbly7k46fvJTgihr/843/keIJZll7M7GGxSpSwGUqtDsL9zdxyYh/te6JxVGnJQHjjYrqr0XMLcC7greT3MUF8QiKGtF1NLpNS8uuuQoqq7MwaFtMm41BUZScu2IfzRsdrP/o6hBCcOTSGnXmVVNmcBLQQjxFCMCElDF+jnqW7i/h2i5tZ405i6ORTyd+/h9+/+4j1v3zFhsXfEp86iMmzL2b4tDMPk2IXQhAb5ENGpY2XZB+umTqV4MpdkLYQKrLUbMIv3DutJSmVXldgDLbRN/D5kryGmcOij19l4dv/IXnwKK585GX0PkGsWbWPxFBfksKaF2+sT2v9y4hk7Xui0a1ol4tJCPGblPKEThjPEdHeGUR+hY3nftnVZPB0a04Fv6YVMjY5lAkp3jdor6x14paS20/uS4gWcDyM9IIqXlm6h7hgX6/ShNMLqvhpewHBvkbmDI9riGHYaqrZsOgbVnz3EQX79+AbENxQwxIWk3DYfgor7fhZ9Fw7NYVIfzOUZMCeRZC3WcUn/CJUVX1TSKkMSkQ/GHsNP6dX8fOOAqL89Hz+4iOs/elLRk6fybw7/4XBZGLJrkK25lRw8djEFqvs8yttDE8I4sKxvZpdR0OjsziiGYQQorGTVweMAo4LZ3p+hY3fdhXRK8yXccneN7m3O91U2V3ccmIfzTg0Q2pUAFNSI1iVUeJVNlhqVAAWo57vtuTy6YYszh4eR4ifCYuvP5NmX8zEWReRsWUtK775kGVfvstvX7zNgLEnMGnWRfQbPaUhiyky0Eyp1c7Lv+7huhN6ExfeB8L7QHUh7FsBmUvB7VQd5hq3E/W4lXFIGAMjLqPEJlm0s4AAYef1B/9K+qbVnHrJTcy49BaEEJRU29maU8GQuKAWjYPD5QEJpw3S0lo1uh/exCD28meHQhewF3hMSrmi84fXNjpyBlHjcPHx2ix0Ai4cm+i1Ho7L4yG3zMbF4xMZneS9UTkeqXW4eebnNDweCPTxLvW3sNLG15tykUjOGhZ3WM8JgPKifFYtnM/qhZ9SVVZMWGwik2ZeyJhTz8EvMFitU+PA4fZwzZQUUiIaNSdyWCF7A+z+EWpLwRQA5gBlHFJPgcHngk7PB6v3s+KP7Xz971spzj3ABbf/gzGnzAGUW/LrTbkUVNq4fGISPi18d3LKapgxWEtg0Og6jihIfSzRUQbC45F8tSmHvAobF4yOJzLAu4yX+kyUkwZEcuaQmOOiUvpI2VNYzX+X7CEu2MfrivTyGgdfb8qlxuHizCExzVYcu5wOtq74hRXffsDe7X9gNFsYOX0mk2ZfTHyfgVTZnFTZXFwxKYlBsYckKHjcqpPc7p+geDcMOhv6zgAh2Fts5YFXv+C7Z27H43Fz5SMv02fon/UPmcXVfLc5jxP6RjC8BS0lq92F2yO59/T+miCfRpfRLgMhhBgDZNX1bEAIcRkqQL0feNSLlqNHnY4yECv2FLNhfxmnDIxiYIz3Ynp5FbX0jQrgyknJbZLfON75elMOK9KL21QjYrW7+HpTDqVWB6cMjKJ/K6KHORlp/P7th/zx63c47LUkDRzBpFkXkTr+ZCodapbY5IxPSrBXgUXt3+2RXPXwi3z09L0EhUdxzeP/O6iAz+2RfLB6P0LAxV6ktV4+oRfDtbRWjS6kJQPRUtrG/wBH3Q6mAk8C76FkN17r6EF2F/YUVrNhfxlD4oLaZBxKrHbC/E1cNC5RMw5t5LRB0YT4Gqmo9b6ng5/ZwHmj4okN8uGn7QVsPFDW4vpxvftzwe2P88jHyzjr+vupLi/lw3/fzTNXnszGL1/l1QVr+G1X0eEigUI0GAcpJbfd/yjvPXEb8amD+OsL8w+r7t6cXU55rZOpqREtfg/Kapwkh/sxND7Y63PW0DjatGQg9I1mCXOB16SUX0gp/wb06fyhHX1KrQ5+2VFAdKCFqX29V1qttrvweODKScmd09msh2Mx6rlwbCKVtU5cHo/X25kNes4aHkvvCD+WpRfz+57iJlVgG+PjH8gJ51zBfW/9yLX/fINe/Yey9NPXeOvWmVx96Tz+9cZnuN2Hj8HpdHLl1X/h5aceY+gJZ3DDv9/BP+jgGUeNw8WavaX0CvNtsX90fVrrnOGaWqtG96ZFAyGEqL/anQT82mhZj7sKOt0eFmzJQ68TnDEkGoOX2j0Ol4fyGgdXTU7yOlahcTgpEf5M7x9JQWXb+jkY9DrOGBLD4LhA1u8vY9HOQjye1uNqOp2O/qOncPVjr/Lgu4uYfv7V5O36gwevnUuv1H688OJLVFZWAlBeXs7pp5/Ou2+/xfhzr+XyB55rUsxxdWYpTreHqa1IZRRW2RmTFEJiC7URGhrdgZaugh8DvwkhvgFqgeUAQog+9DB1Vykl6/aVUVbj4PTB0S0WbzXG45HkVdZy9sg4+kR63zBIo2lOGRhFmJ+J8hpHm7bTCcGJ/SIZlxzKjrxKvq9rAestodHxzLz6Lh75aBlz7/wXLp2Z2/56K3Fxcdxwww1MnDiRZcuWcfpNj3H2X25vMvmguNrOtpwKhsYFtSi053CpcWlprRrHAs3OBKSUTwghFgMxwM/yz7m7DrjlaAzuaPHxuiyyy2qZ3CechFDv7urq5Rsm9wlnUu+e3/jnaGAx6rloXCIvLk7H32zA4KXeFahq6XoxvKW7ivhqYw6zh8W2KTvIaDIzbsY5jD31bP5Yv56MZV/y9ttv4+Pjw53PvYeMGdikBpeUkmW7izAbdIxvpZiyqMrG6UNitPoYjWOCFl1FUsrVTby3u/OGc/Qpr3Hw5opM4kN8GJkY7PV2hVV2ekceP41/jha9wvw4qX8Uv6YVEBfSdhfMsPhgfI16ftpewOcbsg+quvYWIQQjR48mpu9gZl17LxNSQvhut5W4gKYL3vYWW8kqq2Va34gWDZJSazUeF50ENXoGx0/XmmYI9jXx2qWjGZMU4vWFvqzGgb/ZwCXje3mt6qrhPScNjCQi0ExZG11N9aRGBXDW8FiqbC4+3ZBFqbXt+6nXb6qUFr5Oq8bfbEDXxPfD5fGwLL2YUF8Tg+Oal4CXUlJidTBnRNtmNRoaXYl2dQN6R/h7faGvb/xz1eTk47Lxz9HAbNBz4ZheWO2uNsUSGpMQ6su5I+NwuSWfbcgiv6L5ntgtERVoITrQQkgzIo1bsiqoqHUytW94i2mtpVYHyeF+DIkLbtc4NDS6As1AtAGn20NxtYPLxvdqUzc5jbaTGObL7OGx5FXYsDvd7dpHZKCFC0bHYzbo+eKPbPaVWNu1n+bu+OvTWpPCfFvsH+3xSGqdbs4eoaW1ahxbaAbCSzxSkldhY+bQGAa14ErQ6DhO6BvJxeMSKaq2Y7W3T2E+2NfE+aPiCfE18d3mXNLyKztsfKsyS3B5vE1rDfU6AUJDo7ugGQgvyauwMTophOn9Irt6KMcVo5NCuWZKClV2F5VtqLRujJ/ZwLmj4ryuuvaGoio723MqGRof3GJGksPlAQGnDdbE+DSOPTQD4QWFdY1/zh2pNf7pCvrHBHLz9D64PB5KrW0rpKunvuq6T4Q/y9KLWeFF1XVzSClZll6E2ahrVQa+qNrGqQOj2tRoSkOju6AZiFaorHVi0guumJikZZ90IQmhvtxyUioWk57CqvYFnA16HacPiWZIXBAb2lB1fSiZxVayy2oZnxLWalproJbWqnEM06kGQghxmhBilxBijxDiviaW9xdCrBJC2IUQd7Vl26OBzemm2u7i6skpWmFTNyAywMLN01OJDDCTW1HbrhmATgim94tod9W1y+NheXoxYX4mhhwqEd6IP9Na47QbC41jlk4zEEIIPfBf4HRgIHChEGLgIauVArcCz7Rj207F5fFQWGVj3thETTOnGxHkY+T6E/qQGulPTnktnva1zGV8ShjT+0Wwt9jKVxtzsHmZKbW5Lq11Smp4i+7G+rTWwS0YEQ2N7k5nziDGAnuklJlSSgfwCXBW4xWklIVSynXAodHHVrftTKSU5JXbOGVAdJuqqzWODj4mPVdMTGZ0r1Cyy2pwt8NNBDA0PpgzBkdTWGnn8w3ZVNlaDoJb7S7W7i0lOdyv1bTWGoeW1qpx7NOZBiIOyGr0Orvuvc7e9ojJr7QxKC6IGYOjNRmNborJoGPumARO7B9FTnltuwvqDqq6Xp/dYtX16rq01impLccUCqtsjEsJ09JaNY55OtNANHVl9fZWz+tthRDXCiHWCyHWFxUVeT245lCNf8xcODZBa/zTzdHpBDOHxnD2iDjyKmxeu4kOJSHUl3NHxeGRzVddF1XZ2ZZbybCE4GarqkGltQohmDEoql1j0dDoTnSmgcgGEhq9jgdyO3pbKeVrUsrRUsrREREtFyy1RrXNhZRwldb455hBCMHUvhFcPqEXJVYH1e0sqIsMsHD+qKarrqWU/La7CB+jnnFNtSVtRFG1jVMHaWmtGj2DzjQQ64BUIUSyEMIEzAO+PQrbtguHy0N5rYMrJyUT0Yxqp0b3ZXhiCNefkILV7mpzP4l6Gqqu/Q6uus4ospJTXsv4lFDMLWQkVWtprRo9jE4zEFJKF3Az8BOwE/hUSrldCHG9EOJ6ACFEtBAiG7gDeEgIkS2ECGxu284aK4BeJzhnZDx9Iv078zAanUifyABuOSkVUA182oOf2cC5I+OIDVZV1+v3l7I8vYgwP1OLGUlSSkqtDs4eEYfZoKW1avQMRHurSbsjo0ePluvXr2/zdjanm63Z5YxOCtWC0j2A4mo7byzPpLzGSWSAuV1/U5fbw087CthTWA3A2SPiSGwh6FxSbSc6yMKN0/pomUsaxxRCiA1SytFNLdMqqVFqnWOSwzTj0EMI9zdz0/Q+xAb7kFtha1dBnUGv4/TB0YzuFcLIxOAWjYO7Tq31rOFaWqtGz0IzEBo9kgCLkWunpjAwJlAV1LWjVkInBJP6hDOlFbXWoiobY5O1tFaNnodmIDR6LBajnssm9GJcShj/3979xkhVnXEc//7YRRbWCgioKOGPsBhcwH8UlRYratUXprEt1TS0TdP2hbVqtbWtmmi1xkabVvvHGqPG2hdtIlZrUjVR06hJkaql4ALWtlFssWoABQRhWXf36Yt71p3dXNidYXZmd+b3SW64c/bOnWceZufZe+6957y5Yy+dJd4rcSAdnd2AOL/Vo7Va7XGBsJrW2DCK5SdP47zWI3lrZ3v6Qi+fLbvaOX/+UYwf59kFrfa4QFjNGzVKnNd6FF9YNI133m9nb0dpN9T1t3tfJ+PHjmbJnEll2Z/ZcOMCYXVBEktmT+Zrn5jJ9j0dA467NBBf1mr1wAXC6sqCaRO4dNkc2j/sZnuJN9QBvPtBB7OnNNPq0VqthrlAWN2ZNbmZK85uoXGU2Lqr+BvqurqDdl/WanXABcLq0lHjm7j8rBYmNI/mnSInH9rq0VqtTrhAWN2a2HwI3zpzDjMmNQ/6hrqPRms93pe1Wu1zgbC61jymka8vncXCY8azefveAScf2rKrnfNafVmr1QcXCKt7YxobWHHaDM6YO/mAN9TtbvdlrVZfXCDMyEbzvfDEY7hg4VTe2tnOvs6+90pEBNv3dPC5k31Zq9UPFwizRBLnzDuSLy6eztZd+9jT0Tv50LsfdDD7iEN9WavVFU+bZtbP4lmH87GmBn6z6g06u4LmMY3pstajPeKv1RUfQZjlmDd1PJcta6Gzu5tN23Zz2rGTmDbRl7VafXGBMNuP6ZPGcflZLZw0fSLn+rJWq0PuYjI7gCMOa+IbS4+tdhhmVeEjCDMzy+UCYWZmuVwgzMwslwuEmZnlcoEwM7NcLhBmZpbLBcLMzHK5QJiZWS4VM5PWcCdpK/CfEp8+GdhWxnBGMueiL+ejL+ejVy3kYkZETMn7QU0ViIMh6W8RsajacQwHzkVfzkdfzkevWs+Fu5jMzCyXC4SZmeVygeh1T7UDGEaci76cj76cj141nQufgzAzs1w+gjAzs1wuEGZmlqsuCoSkJkkvSnpZ0kZJN6X2myW1SVon6SlJR6f2T0taI2l9+ves6r6D8iohH4tT27r0nM9W9x2UT7G5KHjedEm7JV1dnciHRgmfjZmS9hZ8Pu6u7jsor1I+H5IWSlqdtl8vqal67+AgRUTNL4CAQ9P6aOAF4DTgsIJtrgDuTusnAUen9fnA/6r9Hqqcj3FAY1qfCmzpeTzSl2JzUdD2MPAQcHW130OVPxszgQ3VjnsY5aMRaANOSI8nAQ3Vfh+lLnUx5Whk/1O708PRaYmIeL9gs2Yg0vZrC9o3Ak2SxkTEvkrEO9RKyMeegvamnvZaUGwuACRdCLwOfFChMCumlHzUshLycS7QFhEvp+e/W6lYh0JddDEBSGqQtI7sr9+nI+KF1H6LpM3ACuCGnKd+HlhbK8WhR7H5kHSqpI3AeuCSiOisQthDophcSGoGfgDcVKVwh1wJvyuzJK2V9JykpZWPeGgVmY+5QEh6UtLfJX2/KkGXS7UPYSq9ABOAZ4D5/dqvBW7q19YKvAbMrnbcwyEfqX0e8CLQVO3Yq5EL4KfARWn9Rmqsi6mEfIwBJqX1U4DNFHS/1NIyyHxcDWwiG6NpHLAaOLvasZe61M0RRI+I2AE8C5zf70e/JztaAEDSNOCPwFci4rVKxVdpg81Hwfb/IOtamT/UsVXaIHNxKvATSW8AVwLXSbqsMhFW1mDyERH7InWjRMQasj+o5lYuysoZ5OfjTeC5iNgWWdfsE8DJlYqx3OqiQEiaImlCWh8LnAO8KqmlYLPPAK+mbSYAjwPXRsSqykY79ErIxyxJjWl9BnAc8EYlYx4qxeYiIpZGxMyImAn8HPhxRNxZ0aCHUAmfjSmSGtL6sUAL2fmZmlBsPoAngYWSxqXfmU8Br1Qw5LKqi5PUZFfe/DZ9kEcBKyPiMUkPSzoO6CYbJvyStP1lwBzgeknXp7ZzI2JLpQMfIsXm45PANZI+TD+7NCJG+hDHPYrNRa0rNh9nAD+S1Al0kZ2feq8agQ+RovIREdsl3Q68RHbi+omIeLxKsR80D7VhZma56qKLyczMiucCYWZmuVwgzMwslwuEmZnlcoEwM7NcLhA2Ykm6Q9KVBY+flHRfweOfSfpOGV/vAUnLy7W/gv1eV7A+U9KGQcaySVJZLr+VdJWk/0qqmXs67OC5QNhI9jywBEDSKLLhDVoLfr4EGAk3Ol438Ca5vhcRZRleOyLuIH8sMqtjLhA2kq0iFQiywrAB2CVpoqQxZONGrZV0g6SXJG2QdI8y8yS92LOj9Jd7W1o/JQ08tyYdlUzt/8L720bSs5JuUzaHwL96Bq9Ld9auVDaHwIOSXpC0SNKtwFhl8wr8Lu2+QdK9yuYTeCrdwXtA6Yjil5Kel/R6z5GOpDNTnCtTPLdKWpHiWy9pdqnJt9rnAmEjVkS8BXRKmk5WKFaTjdd/OrCIbNjlDuDOiPh4RMwHxgIXpDGlDknDQwBcDKyUNBr4FbA8Ik4B7gduKXzdQWzTGBGLycZq+mFquxTYHhELgZvJBrYjIq4B9kbEiRGxIm3bAvw6IlqBHeSMibUfU8nuer8AuLWg/QTg28AC4MvA3BTffcDlg9y31aF6GWrDalfPUcQS4HbgmLS+k6wLCmCZsmGXxwGHk83x8SdgJXAR2ZfpxWk5jmwgwqclATQAb/d7zYG2eST9u4ZsQh3Ivrh/ARARG3qOVvZjU0Ssy9nHQB6NiG7gFUlHFrS/FBFvA0h6DXgqta8Hlg1y31aHXCBspOs5D7GArItpM/Bd4H3gfmXTPd4FLIqIzZJuJJv0COBB4CFJj5BNAvNvSQuAjRFx+gFeUwNs0zN3SBe9v2Mq4j0Vzj3SRXbUU+zztJ/27oLH3fg7wA7AXUw20q0i61J5LyK60kBxE8i6mVbTWwy2SToU+OgqpDSMexdwPVmxAPgnMEXS6ZB1J0kqPPE92G36+wvZ0QqSjicraD0+TN1WZsOKC4SNdOvJrl76a7+2nWlM/h3AvantUbJRNgs9CHyJrLuJdM5iOXCbpJeBdfSeCGew2+S4i6yotJHNSNdG1g0GcA/QVnCS2mxY8GiuZhWQhoseHRHt6cqhP5OdLO4ocX8PAI9FxB/KGONXybrianICJCue+x/NKmMc8EzqShLwzVKLQ7ITuFnS5HLcCyHpKrI5DR4+2H1Z7fARhJmZ5fI5CDMzy+UCYWZmuVwgzMwslwuEmZnlcoEwM7Nc/wdbRomJXND9ngAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "wl = engine.wavelengths\n", "\n", "rad = separate_output.radiance[:, 0]\n", "std = np.sqrt(separate_output.radiance_variance[:, 0])\n", "plt.plot(wl, rad, 'C0')\n", "plt.fill_between(wl, rad - std, rad + std, color='C0', alpha=0.5)\n", "\n", "rad = simultaneous_output.radiance[:, 0]\n", "std = np.sqrt(simultaneous_output.radiance_variance[:, 0])\n", "plt.plot(wl, rad, 'C1')\n", "plt.fill_between(wl, rad - std, rad + std, color='C1', alpha=0.5)\n", "\n", "plt.plot(wl, radiance, 'k')\n", "\n", "plt.xlabel('Wavelength [nm]')\n", "plt.ylabel('Sun Normalized Radiance')\n", "plt.legend(['separate wavelengths', 'simultaneous wavelengths', 'HR'])\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 2 }