{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Inelastic Radiative Transfer Calculation with MC\n", "\n", "In this example we perform a full radiative transfer calculation that includes rotational Raman scattering. The \"inelastic radiance\" (with full Raman scatter) is plotted alongside the \"elastic radiance\" (which neglects the wavelength shifts of Raman scatter) in the presence of deep Fraunhofer lines. " ] }, { "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.InelasticRayleigh(), 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 = 1000 # 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 and fix the rng seed for reproduceable results\n", "engine.simultaneous_wavelength = True # enable simultaneous wavelength mode\n", "\n", "# calculate the elastic radiance in parallel\n", "engine.secondary_output = 1\n", "# wavelengths where the optical properties are cached\n", "engine.optical_property_wavelengths = np.arange(385., 407., 0.5)\n", "# incident solar spectrum\n", "engine.solar_irradiance = sk.SolarSpectrum().irradiance(engine.optical_property_wavelengths, fwhm=1) \n", "\n", "# Choose some wavelengths to do the calculation at\n", "engine.wavelengths = np.arange(391.0, 400.0, 0.5)\n", "\n", "# And do the calculation\n", "engine_output = engine.calculate_radiance()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcAAAAERCAYAAADsR3f5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAC1RElEQVR4nOydeXxU5dXHv2f2ySSTPQRIwiaL7CKgKIKKiFYRd1utxdpq1eqr1b6tVWut1lat1tZXq7VqrZbWrSqoKCICAoKy7/seyL4ns8887x93EkIyk0wmCQlwv5/PfDJzl+c+k0zm3HOec85PlFLo6Ojo6OicbBi6egI6Ojo6OjpdgW4AdXR0dHROSnQDqKOjo6NzUqIbQB0dHR2dkxLdAOro6OjonJToBlBHR0dH56REN4A6Ojo6Jxgi8pqIFIvIphiOnSQia0QkICJXR9jvFJFDIvJ858y269ANoI6Ojs6Jx+vARTEeewC4Cfh3lP2PAYvbP6Xuh24AdXR0dE4wlFJfAeWNt4nIABH5TERWi8gSERkSPnafUmoDEGo6joicDvQAPj8W8z7W6AZQR0dH5+TgZeAupdTpwM+Bv7Z0sIgYgGeA/z0Gc+sSTF09AR0dHR2dzkVEEoGzgHdFpH6ztZXT7gDmKqUONjrnhEI3gDo6OjonPgagUik1ug3nTADOEZE7gETAIiK1Sqn7O2OCXYEeAtXR0dE5wVFKVQN7ReQaANEY1co5Nyil8pRSfdFCpm+cSMYPdAOoo6Ojc8IhIv8BlgODRSRfRH4E3AD8SETWA5uBGeFjx4lIPnAN8DcR2dxV8z7WiC6HpKOjo6NzMqJ7gDo6Ojo6HYqI2ETkWxFZLyKbReS3EY45V0SqRGRd+PHwsZ7nSZ0EYzAYlN1u7+pp6Ojo6BxXuFwupZRqyYHyAucrpWpFxAwsFZFPlVIrmhy3RCl1aefNtGVOagNot9upq6vr6mno6OjoHFeIiLul/UpbW6sNvzSHH91uvU0Pgero6OjotBWTiKxq9Li16QEiYhSRdUAxMF8p9U2EcSaEw6Sfisiwzp50U05qD1BHR0dHJy4CSqmxLR2glAoCo0UkBfhARIYrpRo3514D9AmHSb8DfAgM7KwJR0L3AHV0dHR0Og2lVCWwiCbNuZVS1Uqp2vDzuYBZRDKO5dx0D7AJfr+f/Px8PB5PV09F5xhgs9nIycnBbDZ39VR0dE4YRCQT8CulKkXEDlwAPNnkmGygSCmlRGQ8mkNWdizn2akGUEQuAv4CGIFXlFJPNNk/BPgHMAZ4UCn1dHi7DfgKrVedCXhPKfWb8L5HgFuAkvAwD4TvHhCRXwE/AoLA/yil5rV1zvn5+SQlJdG3b19O1P53OhpKKcrKysjPz6dfv35dPR0dnROJnsA/RcSIZtjeUUp9LCK3ASilXgKuBm4XkQDgBr6rjnFheqcVwoff+A5gKpAPrAS+p5Ta0uiYLKAPcDlQ0cgACuBonEIL3K2UWhE2gLX1xzYaayjwH2A80Av4AhgUjkNHxOFwqKZZoFu3bmXIkCG68TtJUEqxbds2Tj311K6eio7OcYOIuJRSjq6eR3vpzDXA8cAupdQepZQPeItw6516lFLFSqmVgL/JdlUfGyb2FNoZwFtKKa9Sai+wKzyHNqMbv5MH/W+to3Py0pkGsDdwsNHr/PC2mGglhfZOEdkgIq+JSGpbricit9an7gYCgVino6Ojo3NC4Ha7efnll3n//fe7eipdTmcawEi31jHHW5VSwbB0Rw4wXkSGh3e9CAwARgMFaIKNMV9PKfWyUmqsUmqsydQ9c4DOOuusuM+96aabeO+999p83qJFi/j6668bXr/00ku88cYbcc8j1mteeqnWBGLOnDk88cRRS8Qopah2+wnp/Wp1dNpNYWEhDz/8MHl5efzkJz/h3Xff7eopdTmdaQHygdxGr3OAw20dJJxFtAgthXaTUqqofp+I/B34uCOv1x1obIiOFYsWLSIxMbHB+N52221xjaOUQimFwdC2e6vLLruMyy677KhtLl+QSpcPq9mA1WSMaz46Oic7Gzdu5Nlnn2XWrFn4/X6mT5/Ovffey6RJk7p6al1OZ3qAK4GBItJPRCzAd4E5sZwoIpnh4kkapdBuC7/u2ejQK4D6wso5wHdFxCoi/dAKKr/tiDdyrElMTAQ0o3Tuuedy9dVXM2TIEG644Qbqk5ZWr17N5MmTOf3005k2bRoFBQXNxnn00UcZN24cw4cP59Zbb20497nnnmPo0KGMHDmS7373u+zbt4+XXnqJZ599ltGjR7NkyRIeeeQRnn5ayzPatWsXF1xwAaNGjWLMmDHs3r37qOvs27ePU089lTvuuIMxY8Zw8OBBbr/9dsaOHcuwYcP4zW9+03DsZ599xpAhQ5g4ceJRIZjXX3+dO++8E4CPPvqI8WecwRnjTuf6Ky8l/5D23h555BFuvvlmzj33XPr3789zzz3XcP4bb7zByJEjGTVqFDfeeCMAJSUlXHXVVYwbN45x48axbNmy9v1hdHSOE5RSzJs3jwsvvJCRI0fy1ltv8eMf/5ht27Yxe/ZsJk+erK9/04keoFIqICJ3AvPQyiBeU0ptbpwGG64DWQU4gZCI3AMMJUoKbXjop0RkNFp4cx/wk/B4m0XkHWALEAB+2lIGaCz89qPNbDlc3Z4hmjG0l5PfTI+948/atWvZvHkzvXr14uyzz2bZsmWcccYZ3HXXXcyePZvMzEzefvttHnzwQV577bWjzr3zzjt5+GGtwfqNN97Ixx9/zPTp03niiSfYu3cvVquVyspKUlJSuO2220hMTOTnP/85AAsWLGgY54YbbuD+++/niiuuwOPxEAqFms1z+/bt/OMf/+Cvf/0rAI8//jhpaWkEg0GmTJnChg0bGDRoELfccgtffvklp5xyCtddd13E9zxx4kTmffkVtd4gb/3rdZ5++o+8+H9/AWDbtm0sXLiQmpoaBg8ezO23386OHTt4/PHHWbZsGRkZGZSXlwNw991387Of/YyJEydy4MABpk2bxtatW2P+3evoHG94PB5mzZrFn/70J7Zs2ULPnj35/e9/z6233kp6enpXT6/b0amLYOH6vLlNtr3U6HkhWqiyKRuA06KMeWML13sceDyuyXZTxo8fT06O9isaPXo0+/btIyUlhU2bNjF16lQAgsEgPXv2bHbuwoULeeqpp3C5XJSXlzNs2DCmT5/OyJEjueGGG7j88su5/PLLW7x+TU0Nhw4d4oorrgC0wvFI9OnThzPPPLPh9TvvvMPLL79MIBCgoKCALVu2EAqF6NevHwMHat2Ovv/97/Pyyy83G2vv/gPc87N7KSkuwu/zkZPXp2Ed8JJLLsFqtWK1WsnKyqKoqIgvv/ySq6++mowMrYlEWloaAF988QVbtjRU3VBdXU1NTQ1JSUktvmcdneON4uJiXnzxRV544QVKSkoYNWoUb7zxBtdddx0Wi6Wrp9chiEhaDIeFwp1nYqJ7ZoF0E9riqXUWVqu14bnRaCQQCKCUYtiwYSxfvjzqeR6PhzvuuINVq1aRm5vLI4880tDd5pNPPuGrr75izpw5PPbYY2zeHF0AOtY6UYfjSEnQ3r17efrpp1m5ciWpqancdNNNDdduLeyilOKuu+7iR7ffxSWXTmfZkq/44x9+RyCoeZ3Rfh+Rxg2FQixfvhxd8krnRGXLli08++yzvPnmm3i9Xi655BLuvfdezjvvvBMxxHk4/GjpjRmBvFgH1HuBHocMHjyYkpKSBgPo9/ubGbF6g5ORkUFtbW1DZmgoFOLgwYOcd955PPXUU1RWVlJbW0tSUhI1NTXNruV0OsnJyeHDDz8EwOv14nK5WpxfdXU1DoeD5ORkioqK+PTTTwEYMmQIe/fubVhD/M9//tPsXG8gRFVVFb17axUs7/znXw3bozFlyhTeeecdysq0Lkr1IdALL7yQ559/vuG4devWtThvHZ3jAaUUX3zxBRdffDHDhg3jX//6FzfddBNbt27l448/5vzzzz8RjR/AVqVUf6VUv2gP2thKTTeAxyEWi4X33nuPX/7yl4waNYrRo0c3yxxNSUnhlltuYcSIEVx++eWMGzcO0MKl3//+9xkxYgSnnXYaP/vZz0hJSWH69Ol88MEHDUkwjXnzzTd57rnnGDlyJGeddRaFhYUtzm/UqFGcdtppDBs2jJtvvpmzzz4b0MKnL7/8MpdccgkTJ06kT58+R52nlKLC5eO++x/k1pk3MOOiKaSlaesWHn90Azhs2DAefPBBJk+ezKhRo7j33nsBLdln1apVjBw5kqFDh/LSSy9FHUNH53jgo48+YtSoUUydOpW1a9fy2GOPcfDgQV566SWGDBnS1dPrbCZ00DENdFortOOBaK3Q9LZYXUOtx0+Fy4/FdPR9mVKKQEjRO8XeKXe2+t9c53jgr3/9K3feeSdDhgzhF7/4Bd/73veOWhI4lnRlKzQRGQs8iNZG04QWElVKqZFtHUtfA9TpFgRDISrdfkzG5gZORFBK4Q+GsOj1gDonGUopHnvsMX7zm99w2WWX8dZbb53s69qzgP8FNgLRQ0MxoBtAnW5BtTsACgwteHi+gG4AdU4uQqEQ99xzD//3f//HzJkzeeWVV+iuHayOISVKqZhqylvjpP9N6nQ9vkCIWm8AszH6krTBILj9IRIjV2Ho6Jxw+P1+brrpJv79739z77338sc//rHNHZZOUH4jIq8ACwBv/UalVJubm+oGUKdLUUpR6fZhEGhpec8o4A0Eo5Y86OicSLhcLq655hrmzp3L73//e+6//379c3+EHwJD0FSC6kOgCtANoM7xhccfxOsPtej9gb4OqHPyUFFRwfTp0/n666/529/+xq233trVU+pujFJKjeiIgXQDqNNlhJSiwuXHaJAWvb/G6OuAOicyBQUFTJs2je3bt/POO+9w9dVXd/WUuiMrRGRoY3H1eNEDyscRffv2pbS0tM3nvf766xw+fEQY48c//vFRLcI6g8bNtB9++GG++OKLZsfUeQMEQwqjITbrV78OqKNzIrJ7924mTpzInj17+OSTT3TjF52JwDoR2R7Whd0oIhviGUj3AE8CXn/9dYYPH06vXr0AeOWVV+IaJxAIxJWB9uijjzYfKxiiKkrZQzT0dUCdE5X169czbdo0/H4/X375JePHj+/qKXVnLuqogXQPsBvyr3/9i/HjxzN69Gh+8pOfEAw2F7W4/PLLOf300xk2bFhDQ+lgMMhNN93E8OHDGTFiBM8++yzvvfceq1at4oYbbmD06NG43W7OPfdcVq1aBWjyRGPGjGHUqFFMmTKl2XVef/11rrnmGqZPn86FF15IbW0tU6ZMYcyYMYwYMYLZs2c3HPv4448zePBgLrjgArZv396wvbFIb71E04gRI/jlPXc2NPW74pILeezhB7novImcNWYEK75e2vCeHnnwfs6dMJbzzx7Pa397EX9QxSQHpaNzPLB06VImT56M2Wxm6dKluvFrBaXU/kiPeMbSPcCW+PR+KNzYsWNmj4CLn4i6e+vWrbz99tssW7YMs9nMHXfcwaxZs/jBD35w1HGvvfYaaWlpuN1uxo0bx1VXXcW+ffs4dOgQmzZpEon1UkfPP/88Tz/9NGPHjj1qjJKSEm655Ra++uor+vXr19BDsynLly9nw4YNpKWlEQgE+OCDD3A6nZSWlnLmmWdy2WWXsWbNGt566y3Wrl1LIBBgzJgxnH766c3GuvPOO7n/gQcpqvbys9t/zPzP5nLhxZcAEAgG+GzhUr74/DOeeeL3vDtnLm/+41UO7t/HF0tWYDKZKC4ppc7tjUkOSkenu/PJJ59wzTXXkJuby/z588nLi7mP80mLiFiBq4C+NLJhSqnmoaZW0A1gN2PBggWsXr26oXen2+0mKyur2XHPPfccH3zwAQAHDx5k586dDB48mD179nDXXXdxySWXcOGFF7Z4rRUrVjBp0iT69esHHJERasrUqVMb9imleOCBB/jqq68wGAwcOnSIoqIilixZwhVXXEFCQgJAM3X3er788kt+/8STuFwuqiorGHLqqQ0G8JLpMwAYNfo0Dh7QbuiWLPqSH9x8S0PoNS0tnU1btsYkB6Wj052pb2I9evRoPv30UzIzM7t6SscLs4EqYDWN6gDjQTeALdGCp9ZZKKWYOXMmf/jDH6Ies2jRIr744guWL19OQkIC5557Lh6Ph9TUVNavX8+8efN44YUXeOedd1r0imJdS2ssdTRr1ixKSkpYvXo1ZrOZvn37xix15PF4uOOnP+WjL5bQJy+Pp5/4HR7vkc+vxaL1NTQYjQSCgYhzNBrAEwi0Kgelo9Odee6557j77rs577zz+PDDD3E6nV09peOJHKVUh6wD6muA3YwpU6bw3nvvUVxcDGjSPvv3Hx3erqqqIjU1lYSEBLZt28aKFSsAKC0tJRQKcdVVV/HYY4+xZs0agKhSRxMmTGDx4sXs3bu34VqtUVVVRVZWFmazmYULFzbMbdKkSXzwwQe43W5qamr46KOPmp3rcrlRSpNoctXV8vHsD1q93uTzL+Cfr/2dQEAziJUVFfQfMIjiVuSgdHS6I0opHn74Ye6++24uv/xy5s6dqxu/tvO1iOh1gCciQ4cO5Xe/+x0XXnghoVAIs9nMCy+8cJR00EUXXcRLL73EyJEjGTx4cIMS+6FDh/jhD39IKKSVCtR7kTfddBO33XYbdrv9KK8pMzOTl19+mSuvvJJQKERWVhbz589vcX433HAD06dPZ+zYsYwePbpBgmXMmDFcd911jB49mj59+nDOOec0O9dkT+R7N97EBWePIzevD6PHNF8jbHa9mT9kz+6dnHfWOMxmMzfM/CE33vwT3pz1Fr/833upqqoiEAhwzz33MGxY1wsY6+hEIxgMctddd/Hiiy9y880387e//U3v6xkfE4EfisgetBBo3GoQuhySLod0TAgEQxRWeTAapcWG17GNpbCaDWQkdowUjP431+lsfD4fP/jBD3j77bf5xS9+wRNPPHFcl/J0hRySiEwAVhBF8T2eTFD99kPnmFDl9oO0rPYQK0YDeP16PaDO8cMdd9zB22+/zZNPPskvfvGLrp7O8cpM4AVgB/AZ8JlSqmV17lbQDaBOp+P1B3H5gq32+4wVESGkFP6gwmLSDaBO92bTpk289tpr/OxnP9ONXztQSt0GICJDgIuB10UkGViIZhCXKaWaF023QKcmwYjIReF2NbtE5P4I+4eIyHIR8YrIzxttt4nItyKyXkQ2i8hvG+37o4hsC7fA+UBEUsLb+4qIW0TWhR8vxTvvkzks3NGocL9Pg8Te7zNW/MH2t0XT/9Y6nc0DDzyA0+nkoYce6uqpnCjUKaWeDWeCng8sBa4BvmnrQJ1mAEXEiOauXgwMBb4nIkObHFYO/A/wdJPtXuB8pdQoYDRwkYicGd43HxgeXvDcAfyq0Xm7lVKjw4/b4pm3zWajrKxM/2LsIFy+IP5gKOZ+n7FiEMHtb9PNXjOUUpSVlWGz6SKDOp3DsmXL+Oijj/jlL38Ztc5Wp818WP9EKeVWSs1VSt2llBrbwjkR6cwQ6Hhgl1JqD4CIvAXMABq6MCulioFiEbmk8YlKsz614Zfm8EOF933e6NAVQId2jM3JySE/P5+SkpKOHPakRClFlduPiNDRgUoVHr/YbmmXZ2mz2cjJyemweeno1KOU4v777yc7O5v/+Z//6erpnEisEJFxSqmV7R2oMw1gb+Bgo9f5wBmxnhz2IFcDpwAvKKUiubc3A283et1PRNYC1cBDSqklEca9FbgVwGKxNBvQbDY3dEbRaR+fbixgwZ5ieqd0jod1qNLNzy/sR3ay7sHpdD/mzp3L0qVL+etf/3pUMwmddnMe8BMR2Q/U0Y4yiM40gJHuy2OOK4YXM0eH1/g+EJHhSqlNDYOLPAgEgFnhTQVAnlKqTEROBz4UkWFKqeom474MvAxaGURb3pBO7JTUePlyezE9kjqmVCEaB8vrdAOo0+0IhUL86le/YsCAAfz4xz/u6umcaFzcUQN1ZhJMPpDb6HUOcDjKsVFRSlUCi2gkgSEiM4FLgRvC4VKUUl6lVFn4+WpgNzAozrnrtJO5GwswiWDqoMzPSNhMBrYWNu9wo6PT1fznP/9h48aN/O53v8NsNnf1dI45LSUyNjpGROS5cJLkBhEZE+Pw44HycN3fjcCzQFwLrJ1pAFcCA0Wkn4hYgO8Cc2I5UUQyG2V32oELgG3h1xcBvwQuU0q5mpxjDD/vDwwE9nTc29GJleJqDxvzK8lop/cXCIVazPR02s3sKKohFNIdeZ3ug8/n49e//jWjR4/m2muv7erpdBUtJTLWczHa9/RAtGWpF2Mc+9dKqRoRmQhMA/4JxJX132khUKVUQETuBOYBRuA1pdRmEbktvP8lEckGVgFOICQi96BljPYE/hk2aAbgHaXUx+GhnweswPxwEfSKcMbnJOBREQkAQeA2pVTrzS11Opwlu0oxGgztKnpXSjF73WGCIcW1Y3MjHmM2GvD4Q5TUeunh1MOgOt2Dl19+mb179/Lpp59iMJyc7ZZbSmRsxAzgjfCxK0QkRUR6KqVaE/esT/++BHhRKTVbRB6JZ56dWgivlJoLzG2y7aVGzwvRQqNN2QCcFmXMU6Js/y/w37gnq9MhVHv8fLOnjMx2tinbVlhDfoUbALcviN1ijHrswXKXbgB1ugW1tbU89thjnHvuuUybNq2rp9OZmERkVaPXL4fzKxqIIZExUqJkb7R8jpY4JCJ/Q4sMPhnWB4zrTkPvBKPToazcW45StGvtz+MPsmRnKQ6rkTpvkPwKFwN7JEU81mYysK2whrF99Rorna7nz3/+M8XFxcyePftEb9MXaK3urrVERuJPlLwWLSfkaaVUpYj0BP43xnkfxcnpn+t0Ct5AkIXbi0l3NC8vaQvLd5fh8Qe5dGQvzEbhYNgTjESSzcy2wmp9HVCnyyktLeWpp57i8ssvb1Bo0YmcyBimTYmSIjJBREQp5VJKva+U2hkev6BJfXjMxGQARWSiiPww/DxTRPRCOZ1mbMqvwu0LYTVHD1e2RlG1hw2HqhiVk0K200bvFDv5Fa6ox1tM2jpgaW27hKF1uhGVLh+ednb56Qr+8Ic/UFdXx+OPP97VU+lyWkpkbMQc4AfhbNAzgapW1v9mAqtF5C0RuSmcQ9IuWjWAIvIbtKzL+pZjZuBf7b2wzolFKKT4fEsRKfb4U75DSvHltmISLEbOHKCFNHPTEqhw+an1BFo8tyUjqXP8EAiGeGXJXlbtq+jqqbSJAwcO8MILLzBz5kyGDm3a8fGkpCewUEQ2oFUEzFdKfSwit9UnQqLlh+wBdgF/B+5oaUCl1G1KqTHAI0AqWjPs5SLyexGZVF8F0BZiWQO8Ai0hZU14EodFJPKCjM5Jy66SWkprveSkJsQ9xqZDVRTXeJk2rAdWk/ZZzg2Pd7DCxak9IytnW8P1gGP66OuAxztrDlSys7iGBIuRiQMzuno6MfPII48c9fNkRykVMZGxSRKkAn4ax9jb0LzJZ8Pe5XlozbD/BLSpH2gsIVBfeKIKQET0nj46R6GUYsHWIhIs8edUuXwBvt5dRk6qncGNEl4yEi3YTAYOtuDhOW1mthfW6A3Mj3PcviAfrT9MbmoCe0vrjpsw6JYtW/jnP//JT3/6U/LyImq16nQS7W2GHYsBfCeccpoiIrcAX6C5qzo6AByu8rCruJbUhPjDn0t3leIPhjhvcNZR2XMiQk5qAvkV7qgGzmIy4PYFKdHXAY9rvtpZgtsfxGHVbqQOlh8fYe2HHnoIh8PBr371q9YP1okLEZkqIn8XkdHh17d2xLitGkCl1NPAe2g1doOBh5VS/9cRF9c5MVi6swSz0RB32vehCjdbC2oYk5dKWoQM0pw0OzWegKYqHwWlFIdayBbV6d6U1/n4YmsRWeHuQSKws6i2lbO6nhUrVvDBBx/wv//7v2RkHD8h2+OQO9BKHb4vIuejdZdpN7EkwfQDliil/lcp9XNgqYj07YiL6xz/VLp8rNpXTkache/BkGLh9mKSbCbG94u8hpfXsA4Y3cBZTEa2FVRH3a/TvflsYwFGEczh+tFku5n1+ZXdOqxdL3eUlZXFz372s66ezolOiVKqMmyDLgTGdcSgsYRA3wUaN2QMhrfp6PDN3nJA4ha8XX+wkrI6H5MHZTZ8+TUlJcGMw2okv4WQmNNuYpu+DnhccqDMxaoDFUd1D7KbjZTV+Vr0+ruazz//nMWLF/PrX/+axMTErp7Oic4n9U+UUvcDb3TEoLEYQJNSytfo4j6gfZXOOicEHn+QxdtL4vb+ajx+Vuwto296Av0zoudWiQi5qQkcbGEd0Goy4vIFKa31Rdyv0z0JhRQfrjtEgtmIwXD02q8A+8u65zpgvdxR3759ufXWDlmO0mkBpdTs+uciMhY4X0TWhFUkNobLLdpMLGl7JSJymVJqTvjiM4DSeC6mc2Kx7kAlvkAIiym+hkJLdpYSUnBuk8SXSOSk2tlWWENZnS+qwVUoDlW6yexkDUKdjmPz4Sr2ltaRm2pvts9iMrDpcBWjclOO/cRa4Z133mHt2rW8+eabEYW1dTqVWWjrgRs5OjrZZmIxgLcBs0TkebTebQeBH7TnojrHP4FgiPlbC0mJM/Nzf1kdO4trObN/GskxFM/X1wPmV7ijGkCLUVsHHN0NvzB1muMNBPlw3SHSEiwRb4CS7Wa2HK4mGFJxh9g7A7/fz0MPPcTIkSO5/vrru3o6JyMl9Q5Ze2nVACqldgNnikgiIEopXYFUh+1FNVS6/HEVvgeCIRZuLyHFbub0vNSYznHazSTbzRwsd0U1cNo6YDVKqRO9EfEJwTd7ylv8DJmNBnyBEIXVHnqnNPcQu4pXX32V3bt38/HHH5+0ckddzG9E5BVgAZruIABKqffbOlCrBjAsNXEV0BdNAqP+Yo+29WI6JwZKKb7YUkSiNb7C99UHKqhy+7l8dK82qUbkptrZUVRLKKSOWi+qx2oyUlbrazFMqtM9qPb4+XRjQQzhasWektpuYwDr6ur47W9/y8SJE/nOd77T1dM5WfkhMAStLWd9CFQBHW8AgdlAFZquk15prMOBchcHyl1xfSlVuf2s3FfBwKxE+qS3ralQTmoCmw5XU1zrJTuK/p9CqyvUDWD3ZsHWYoIh1dDyLhoOq5mN+VWcMzDzGM2sZZ577jkKCwt577339ChD1zFKKTWiIwaKxQDmKKWayljonMQs3lGCxdT2wnelFIu2F2MQmBTHF1pOOFEiv9wV1QBajAa2F9Z0y8QJHY3CKg/LdpVG/Rs2JslqYl+Z1hbN1g6VkY6gvLycJ598kunTp3P22Wd36VxOclaIyFCl1Jb2DhRL/OlrEekQa6tz/FNa62XDwUoyHG33sPaU1rGvzMWZ/dNJtLU9fOqwmkh3WFosiHfaTGwtqNbrAbspSik+3nAYizG22lGDQVB0j7ZoTz75JNXV1brcUdczEVgnItuPRRnEROAmEdmLFgIVtEbeI+O5oM7xzfLdZRgMEnENriX8wRCLd5SQ7rAwKicl7uvnpiaw6XBV1MxAa7iAurzOR7oeBu127C6pZfPh6ohlD9EwCOwoqmFgj64TocnPz+e5557jxhtvZMQI3R/oYjosIhmLAby4oy6mc3xT5w2wbFdpXOtr3+4tp8YT4OrTc9qV0p6TZmddfiWFVR56t/AleqjSrRvAbkYgGOKDNYdItpvbFD532sxsyK/iOyN6dtm626OPPkowGOS3v/1tl1xf5whKqf0dNVYszbD3hy/oRssxaJBGag0RuSjspu4Skfsj7B8SFjT0isjPG223ici3IrJeRDaLyG8b7UsTkfkisjP8M7XRvl+Fr7VdRKbFMked2Fm9v4JASEVtWRaN8jofaw5UcGrPpHZn8+Wk2LVi1BbkkUxGYXuhXq3T3Vh7sJKCKk9MdZ+NqW+LVunqmrZo27dv57XXXuP222+nb9++XTIHnSOIyD/r1ebDr1NF5LV4xoqlGfZlIrIT2AssBvYBn8ZwnhF4Ac2DHAp8T0SaSiWXA/8DPN1kuxc4Xyk1Cq3r90UicmZ43/3AAqXUQLQ6kPvD1xsKfBcYhuYi/zUehWCdyPiDIRZsLYqo1tASSmnNrs1GAxNPaX+3fKvZSGaStVV9QH0dsHvh9gWZs/4w6Ylt75oiIqBgfxetAz700EPY7XYefPDBLrm+TjNGKqUq618opSqIIL4bC7Hcyj8GnAnsUEr1A6YAy2I4bzywSym1J9w/9C1gRuMDlFLFSqmVgL/JdqWUqtdCMYcf9d9mM4B/hp//E7i80fa3lFJepdReYFd4DjodwJbDVdR6A9jbmIm3vaiG/Ao3Zw1Ib5dgbmNy0xIorPLgD0bugmQ1Gaj2BKjoIo9BpzlLdpbg9gXj/gxYzQY2Harq4Fm1zsqVK3nvvfe47777yMrKOubX14mIoUnkL43YlvOaDxTDMX6lVFn4ogal1EJi02LqjdY2rZ788LaYEBGjiKwDioH5Sqlvwrt6KKUKAMI/6z+VMV1PRG4VkVUisioQCMQ6nZMapRSfbynGaWtb6MobCLJkZylZSVaG907usPnkptoJKThcGTkbVFsn0vUBuwtNtf7iIdmuefXB0LH16l988UWSk5O59957j+l1dVrkGWC5iDwmIo8CXwNPxTNQLAawMtwG7Su0nqB/AWKxHJFWq2P+9Cqlgkqp0UAOMF5EhnfE9ZRSLyulxiqlxppMHeORnOjsKa2joMpNUhtLF1bsLsflC3L+kCwMHZi80CvFjkFa1gc0GwzsKNLXAbsDn20qwNBI6y8e6tuiFVQdu5sapRRffPEFF1xwAU6n85hdVycyIjJBREQp9QZwJVAElABXKqXejGfMWD6RM9ASYH4GfAbsBqbHcF4+kNvodQ5wuK0TDMd6F3Ek9bVIRHoChH8Wd+T1dJqzcFsxdrOxTRl4xTUe1udXMrJ3Mj1iKHhuC2ajgexkW4u1YU67mS0FVfo6YBdzoMzFqv1Ha/3Fj2JvSV0HjBMbu3bt4uDBg0yZMuWYXVOnRWYCq0XkLbTlrfeUUv/XnoL4WLJA68LeWEAp9U+l1HPhkGhrrAQGikg/EbGgJajE1MFbRDLrs3xExA5cAGwL756D9osg/HN2o+3fFRFrWMV+IPBtLNfrCoqLi7nvvvsYNGgQy5cv7+rpRKWo2sPWguo2Jb8opVi4rQSb2ciEAemdMq/c1ARKarx4/MGI+60mA1XuQJdlDupoWn+zI2j9xYvDqpVDHCsWLFgAoBvAboJS6jal1BjgESAVeD1cRfB7EZkUT9JjVAMoIkvDP2tEpLrRo0ZEqmOYbAC4E5gHbAXeUUptFpHbROS28NjZIpIP3As8JCL5IuIEegILw9X9K9HWAD8OD/0EMDWcmTo1/Bql1GbgHWALmqf6U6VU5G/HLqSsrIxf/epX9O/fnz//+c9UVFRwySWXsGnTpq6eWkSW7irFZDC0KYS5paCawmoP5wzM6LT2VbmpCVrfzxbXAaPv1+l8thRoWn9tzRyORpLVxL7yuqg3PR3NggULyMnJYeDAgcfkejqxoZTappR6Ntyi83xgKXAN8E3LZzYnqgFUSk0M/0xSSjkbPZKUUjEFxJVSc5VSg5RSA5RSj4e3vaSUein8vFAplRMeNyX8vFoptUEpdZpSaqRSanhj5QmlVJlSaopSamD4Z3mjfY+HrzVYKdVqqcaxpLKykocffph+/frx5JNPctlll7FlyxZWrlyJ3W7nwgsvZO/evV09zaOo9vj5Zk8ZGW1IXQ+FFN/uLaeH08qQ7M7r3NEj2YrJIOSXRzdwJoOw/QRdBwwGg3z44YdUVlZ29VQi4guE+HDtYVKjaP3Fg8EgKKU1Y+9sQqEQCxcuZMqUKXrT626GiFwjIvVfLvcBPwb+oZQa29axWvIA01p6xDn3k46amhp+97vf0a9fPx577DEuvPBCNmzYwL///W8GDx5M3759+fzzz/F6vUydOpXCwsKunnIDK/eWE1K0SbJoR3EN1Z4A4/umdeoXh8lgoFeKveV6QLuZrYdPvHrAnTt3MmnSJK644gq+973vEQq1SxS7U/hmbxmVLl9cPV9bwmiAncfgpmb9+vWUlZXp4c/uya+VUjUiMhGYhlYO92I8A7X0zbYaWBX+WQLsAHaGn6+O52InE3V1dTz11FP069ePX//610yaNIm1a9fy3nvvMXz40Qmtw4YNY+7cuRQWFnLRRRd1i7t6byDIwu3FZLRx7W/V/grSHBb6ZbRN6igeclPtlNX5qPNGTkq2mQxUuv1UuU+MdcBQKMTzzz/PqFGj2Lx5MzNnzuSzzz7jueee6+qpHUWNx8/cjQVktKPsIRpOm5n1+Z2f3KSv/3Vr6mPglwAvKqVmA3HF2VsKgfZTSvVHW8ObrpTKUEqlA5cSh/DgyYLb7ebPf/4z/fv355e//CXjxo3j22+/Zfbs2YwePTrqeWeccQYffPABW7ZsYfr06bhcXdv9fmN+FW5fEGsb1vD2lbkoq/Uxtk9qh3l/LX3R5aRpSuL5UcohRARpYf/xxP79+5k6dSp33XUXkyZNYtOmTfzjH/9gxowZ/PKXv2Tt2rVdPcUGFmyLTesvHuxmI+W1nd8WbcGCBQwZMoRevXp16nV04uKQiPwNuBaYGxZtj6vGJpaTximl5ta/CK+tTY7nYicU5Xtg9RuwdwmU7sJbVcILzz/PKaecws9+9jNGjBjBsmXL+PTTTxk3blxMQ06dOpVZs2axbNkyrr32Wvz+rvFcQiHF/C1FpNjbdlO1an85iVYTgzqoa783EGRbYQ2+QOQQX1aiFYvJQH4LYVCjQdhZfPyuAyqleO211xgxYgTffvstL7/8Mp9++ik5OTmICK+88goZGRlcf/311NUduxKBaBRVe1i6oyQGpff4EBGQzm2L5vP5+Oqrr3Tvr/tyLZpjdlG4TC4N+N94BorFAJaKyEMi0ldE+ojIg0AsZRAnNrXFsGs+/tVv8veHb2HQKf258667GJCVwKJ/Pc0Xs57jrFGDoY2hmmuuuYaXXnqJTz75hJtuuqlL1nd2FtdSWutt0/rN4Uo3hys9jMlLaZfaQz1KKYqrvQzISowawjQYhJwUeyv6gGa2Hq45LtcBCwoKmD59Oj/60Y8YM2YMGzZs4JZbbjnKu87IyOCNN95g+/btXd6tpF7rz2wyYDLEX/TeGlZT57ZF++abb3C5XLoB7GY0KoR3KaXeV0rtBK0jmFLq83jGjOUb7nvAb4APwq+/Cm87qQkEAsxals+j721gT0E5Z5yayyv3Xc4Fw3sggd2w/AXtQKMZUvIgYyAk50JiD3BkgjH6r/7WW2+lrKyMBx54gLS0NJ577rljlommlGLB1qI292xcvb8Cm8nQYS3Pyup89M1wcM7ADGZ9E139JDctgT2ldVS7/TgjqAzYzAYOV3modgdITmhbK7euQinF22+/zR133NEQUr/rrrswRDEqU6ZM4Re/+AVPPvkk06ZN48orrzzGM9bYXVLH5sPV5LRD8eNguYulu0q54rTeUUtoGrdF64ibraYsWLAAg8HAueee2+Fj67SLmcALIrIDrdTtM6VUu7IGW/2WC5cZ3N2ei5xorF27lu9efQM79hxgzMBefPz7H/CdMwY3MlKNir9DAagtgrI9NHRmEwFnb0jrrz0SMzWjaDmSOHL//fdTVlbGM888Q0ZGBr/5zW+OyXs7VOlmd2kdvZNj795SVutlT2kdZ/RLa1e7q3r8wRC+QIhrx+ZiMxsAQSkV8SYgJ6wJeLDCxTB7c+PbsA5Y6SI5oeP6kXYWpaWl3HHHHbz77rucccYZ/POf/2Tw4MGtnvfoo4+yYMECfvzjHzN+/HhycnKOwWyP5sutRTgsprhv1gKhEF9uK6bS7WdfWR1DsiNXWzVui5aTmtCeKUdkwYIFjBkzhtTU1NYP1jlmKKXq68eHoKkMvS4iycBCNIO4rK21360aQBHJBH6BJjPU8K2olDq/LRc6kcjNzSUjNZknfnEhl0+b3PI/vMEEthTtUY8KgbcG9n8NexcDom2zJmneYlp/xNmLPz50D2UlxTzyyCOkp6dz5513dur7qvUG+HRjIWaDtOlLbPX+CkwGYVRuSofMo6jaw0XDs8kOG+GMRCtuf2QlgXSHBbvZSH6Fm2G9Ihs4g0HYWVQbdX93Yfbs2dx6661UVFTwhz/8gZ///OfE2q/WYrHw73//m9NOO40bb7yRL774AqPx2KmBldR42VFUS6+U+NverT1QSaXbj9Eg7C9zRTWA9ewtqetwA1hbW8uKFSu47777OnRcnY5DKbUNrTPYs+FOYeehFcL/CWhTLWAst+uzwhfrB/wWTQ9wZVsucqKRkZHBsvdf5ooz+8V3tysGzdglZWth0eQc7ac5ASoPwLZP4Nu/I4uf4O+X2bn8rMHcdddd/PtPD0DhJqgpgmDHKVkopdiYX8mTn25jZ3FNmxIYqt1+thfVMLx3cpulkiJR4fLRw2lj8qAj0jMjejujrgOKCLmpWj1gtHU+p83UrfUBKysrmTlzJpdffjm9evVi1apV3H///TEbv3oGDhzI888/z6JFi3jyySc7abaRWbWvHIMQt/dX6wmwcl85/TMcnJKVyP6y6H9PgESrifX5lXHONjpLliwhEAjo63/dmI4shI/lPyxdKfWqiNytlFoMLBaRxW29kE4riIDJqj3sR0IvJqX4z6+u5OIH/sXMXzxJSsUGvjO2jxZNTeoBKX0gtR8kZmlhVHsKGGI3RJUuHx+uPcSG/CrSHJY2t61ae6ASgDF5KW06LxKBUIhab4AfT+yPxXTk3mxQjyQW7SiJel5OWgI7imupdPlJjTB/u9nIoUo3B8vd5KV3fMisPXz++ef86Ec/oqCggF//+tc89NBDWCzxtw6bOXMm8+bN4+GHH2bKlCmcccYZHTjbyPgCIZbtLiUtDrHbepbsKiGkYNKgTAqq3GwvrKG4xhu1kXqi1cT+MhduXxC7peM83QULFmCxWDj77LM7bEydDufXSql3GxXCP41WCN/mD3tMeoDhnwUicomInIamtKBzLBDB5khi9u9/yKgBPbnqqXks3R8AZy8I+KBgPax/C77+P5j/MHx0N3z5O1jzBuxeDEVbtIzVJh5jfcuyJz/bzrbCGnJS7TisbfM43L4gmw5XMTg7iaQ2agVGoqjay7mDMpsZqdxwvV8oihZcbngd8ECUcggRwWIysGhHccT9XUFtbS23334706ZNIzExkeXLl/Poo4+2y/iB9l5ffPFFcnJyuP7666mubrVtb7vZVliNxxeMu+4vv8LFjqJaTu+TSrLdTF74772/LHqpQ31z7ZY6AcXDggULOOuss0hI6F43SscbIpIrIgtFZKuIbBaRZnkkInKuiFSJyLrw4+EYh++wQvhYvvF+F15ovA/4P8AJ3BPPxXTix+mw8ekTM5n4Py9z6QNvsPjPtzBqQM+jEmcACAXB74ZDa+HACo6SSUzMgpRcKq29+PyAYnOFhZSULCzW+NZt1udXEggpTs9rf7JAjceP02biwmHZzfbZzEb6pjkorfVGzPRMtptJspnIL3czKicl4vgZDisb8qsoqfF2Wo1arCxdupSZM2eyd+9e7rvvPh577DHs9vgzJ5uSkpLCrFmzmDRpEnfeeSdvvPFGh43dFKUUi7eXtPnmqZ5QSLFoRwlJNhNj+2ifowSLiawkK/vK6hjfL3rXRYMBthfWdFjdaWlpKevWreOxxx7rkPFOcgLAfUqpNeFw5WoRmR9BumiJUurSNo5dXwg/FXiyUwvhlVIfK6WqlFKblFLnKaVOB8pbO0+n48lMSWT+H28mKcHKtF/8g12HIpRjGoxgSdSMXcP6Yg44exEK+Dm8fRU7F7xO/11vcHX1P7jkwBNM3Pccow6/Q9/yZWTWbifJU4AlUKsl5kTBFwix/mAl/TMcpLdT6y0UUlS4/Hx3fF7U1PeROcnURml5JiLkpNrJb2Ed0GAQDAJf7y5t11zbS0FBAVOnTgVg8eLFPP300x1q/Oo5++yzefjhh3nzzTeZNWtWh49fT2G1h31ldSRHuDGJhQ2Hqiir9XHOwIyjMoj7pjsorPK0qPzgtJnZ2IFt0RYuXAjo7c86gnBt3prw8xo0RaDeHTR8fSH8tHAhfCpxFsJHvW0Laytdizbpz5RSm0TkUuABwA6cFs8FddpHXo8UPn/qh5xz98tM/d/XWPbcT+iV0bo4R403yObDHqrcZhyO3hgNosW2lcIU8pLq3kdW3TZAoTAgKgRiwG1OxmVOp86cTq0lE6/ZideUxLICwRMIMbZv+72/ohoPZ/RLY2BWYtRj+rbSWzQ3NYGtBTWU1vqieniZiVaW7y5jyqk9SIzTY2kvTz/9NH6/n88//5wBAwZ06rUefPBB5s+fz+23386ECRPo379/h1/j273lmAyGuJJfXL4AK/aUkZtm55TMo//2fdIT+HZfOQfLXQyM4uHZzUYOV3qocPk7RHJpwYIFJCUlxdy56STHJCKrGr1+WSn1cqQDRaQvmr2IJFc0QUTWo4mX/zwsa9caQbSKhGtEpPE/cpuL4Vv6FngVTWH9W+A5EdkPTADuV0p92NYLnWhUuHyE6vykJCuMx1gu5dQ+WXz6xE2cf98rTPvlP1j87C2kOSOvWQRDin1ldewursVkFJw2k5ZwU48IAaONgNGGp+nJKoQp5CPJW0Cqex9G5UdhwKeMvHT4EoZZ67i6Yh6u2jRqLRm4zOl4zMl4TYl4jUn4jQlHXysCLl8Aq8nAJSN7tvgl2ivFjsVkwB8MRaw1bFwPGM0AmowGAiHFyr3lnDckK+IxnUlJSQkvvfQS119/facbPwCTycSsWbMYNWoU119/PUuWLMFs7rhmAB5/kBV7ykiPM/nl691l+IMhzh2U1exvn+20YTUZ2FcW3QBqbdEU+8s6RnNwwYIFTJ48uc3ZtycpgViyLkUkEfgvcI9SqumC9Bqgj1KqVkS+A3yIJmTeGrOByvD53rZMuikt/aXHAiOVUiERsQGlwCntrbw/USiq8lBWVEOgroQBWYlkO22d2v6pKeOG5PDhY9/nO7/6J5c+8Abz/3gzjia9OyvdPjYdqqLOGyTRamq7KrcYNOOI7ahP2eLqbEqDDn6UtRNRIZI9h0hz78WgAmhrjoIQIigm6swZ1Fp7UG3tgduchsecjNuUTMBgQwGltT5mTujTahKN0SAMzk5iR2FNxJBrks1MSoKZg+UuxrSwJpnusLBwezFnnZLeKc2aW+LZZ5/F7XbzwAMPHLNr9unTh7/97W9897vf5dFHH+3Q9a1Nh6rwByLfkLRGYbWHzYerGZOXEtF4GQxCXloC+8vqojZBALCajGw6VM1p7VyHPnDgALt27eKnP/1pu8bROYKImNGM3yylVDMBhcYGUSk1V0T+KiIZSqnW1ilywmK47aYlA+hTSlsEUkp5RGSHbvyOxmIQggibD1Wzs6iWAZkOeqXYj5khnDLmFP7z4HVc8+h/uOqRWXzy+5kYjQYCoRC7i2vZV+bCYjJETByJl5CCjyryyLPUMspRQUA07zESooKYQ24ya7fRq3o9SjTjiAoRMFopCqVwekYfRgVGQ2E6JKSBPQ3Mkccb3svJhoOVUeeWm5rAtsKWW2TZzEbKan1syq/i9L7HTtayoqKC559/nmuuuYYhQ4Ycs+sCXHfddcybN4/HH3+cCy64gMmT29/LXinFwu0lcX22lFIs2l5MgsXYYpJL33RHuC9t9LB2st3c6t88FnT5o45FtDuWV4GtSqk/RTkmGyhSSikRGY+WkxJLn+mvRWSEUmpje+fZkgEcIiIb6ucKDAi/FkAppUa29+InAmaTAbPJQCAYYmtBDbuKa+mfmUjvFHuHtAVrjSsnDef5u6Zzx1/m8PaiDUw781Q2HarCGwiRaDNh6ODw7Jq6DA75HNyZvbm16CZKjPiNCfiNCTRtV638PhKCNYwz70Y2bKfec4SglsST1DOcvNMb0vpBUk/y0h0g0dui5aba2XioiuIaDz2ToyeWOO0m5m8t4rS81LZ7xXHy3HPPUVNTw4MPPnhMrhfp+kuXLuX73/8+GzZsaHebr4Plboqq3fRq4fccjS0F1RRVe7lwaI8WvfA+6fXlEHVRDWB9W7TDle6Gcpl4+PLLL8nKymqm1akTN2cDNwIbRWRdeNsDQB6AUuol4GrgdhEJAG7guyq2jKaJwE0ishctBBq3TWrJAJ7a1sFOZkxGA067gWAwxI6iGnYX19I3w0FOqr3TQ20/mT6eFz/6hgdf/YLkzCwcNlOH1OU1RSmYXZFHlsnNmYnRC9NjGagmIAzrnY0tJaHZPkJ+rXaxYj8EvVrnHJuT9Jxx9MdGtS8Xa4TSjfq2WAfL3S0awCSbmfwKF7tKajsshb4lqqur+ctf/sKMGTMYObJr7hsTExP597//zYQJE7jlllt4991329VgfcWesriSX7z+IMt2ldEz2caQ7JZ/9w6ricxEK/vKXIxtwVtXCvaW1sZtAJVSLFiwgPPPP/+YNZ0/0VFKLeWoGqyIxzwPPB/H8BfHNakItCSIu7+lR0dN4ETDaDSQZDNjNRvZXVLLkp2l7CyuwRNoU4/WFlEo3P4gxTUethVW8/Wecq6cMpZ9heWsWL+r0wzuVncKuzzJXJp6AKPEn3ru8gdJdVjoFUk1QASMFrAla63iUvpo5RwGC7J7IZfUvMt5u59iZME79KjZrJVrhLFbjGQkWmIqjk6wmPhiS9ExaY/217/+lYqKCh566KFOv1ZLjB07lscff5z//ve/vPrqq3GPU+sNsGp/eVzJLyv2lOP2Bzl3cGZMxqZPegIFVW68Lfz/JNlMrM+PXx5p27ZtFBQU6OHP44Sw/UkBpocfKfHapE6N0YnIRSKyXUR2icj9EfYPEZHlIuIVkZ832h61i4CIvN2oc8C+evc6rFfobrTvpc58b61hNAhJNjN2s5F9pXUs2VHKtsJq3C3UNUVDKUWtN8DhSjcbD1WyeHsJS3aWsO5gJfkVbkJKceG4QQzKy+K1j5cT6EBj25g5FXk4jT4mO+NfCg6FFMGQYmhPZ9vCs2Y7OHuTkNmXGmMy6XV7GFnwHufueYYzDvyd3MpvcHiLyU21U1DlIRBsWUcxNcHM7pJaDlV2rlp8XV0dzzzzDBdddBFjx7a5VWGH8/Of/5zzzz+fu+++m23btsU1xvqDFYQUbV7rLq31sv5QJSN6J5OVFFvzhb7pDkJK8+qjkWgzcSDcFi0e9PW/44uwPZgFZIUf/xKRu+IZq9PyfcN1hC+gVevnAytFZE6TTgDlwP8Alzc5PWoXAaXUdY2u8QzQ+NZvt1JqdMe/m/gxGIREm5lQSHGw3MWBchc5qXb6pDtwRNHcC4Y0g1fl9lNa66XC5SMYbgNmNAhWkwFrhHT2W2acxf/+34d8unwL088Z0aHvY7/XwXpXOtel78FiiF+kt9YbYFCPxLjr8JLtZkJiwm1O1bxFpbAGqhlSPA9Q4BvA2tAZ1BTvI71HLiFD5FBwfXu0r3aUcP0ZfeJ+P63x97//ndLS0i73/uoxGAy88cYbDaURy5cvx2qNvZFBKKRYvL2UlDYmv2iJLyVYjQYmDEhv/YQw2ck2LEYD+8vqOCVKnahBBIVWAhNPSHvBggX069ePfv36tflcnS7hR8AZSqk6ABF5EliO1qmsTXSmBzge2KWU2qOU8gFvATMaH6CUKlZKreRIv9H67a12EQhnGV0L/Kfz3kLHUW8IE60mDld6WLazlI2Hqqj1BvAHQ5TXedlTWsu3e8v4clsR3+wtY1thNVVuPzaTkSSbmSSbmQSLCWOUO++zRvbn1L7Z/OPjFfg72AucU94HuyHA1ORDcY/h8QdJspm0ZJY4sYRDzL56D08EnymJGls2NbaeDLZXYyCEN3895+9+MmKotJ50h5U1+yspr/PFPZ+W8Hg8PPXUU5x77rndqrly7969efXVV1m7dm2bk3L2ltVRVudtc+uzncWatz1hQHqbVEOMBiE3zc6+VtQhjCJsL6xp05wAgsEgixYt4vzzT1p1t+ORcLZcA0FaWW+MRkudYDbSoODanBgybnoDBxu9zieObt0tdBE4By2Fdmejbf1EZC1QDTyklFoSYbxbgVuBdjcejgcRwWE1oZSiqNpDQaX7qD+dxWjAYY1PVFREuOXys7j3z+/zybJNXD55VIfMuchvY3ltFpekHsRhjE+GKaQUvkCIMXmp7W4ckO20srO4FkuEtU6L2cQAWw1rfLlcYSkg3bWHHjVbQcBtSqXMMYByez+qrdm4zakYBJbvLuWSkb3aNadI/OMf/6CgoIB//etfHT52e5kxYwa33347zzzzDFOnTmXatGkxnbdsV2mb15h9gRBLdpaSmWRleO+2azL2TXewu6SO8jpf1LZ7yXYzG/KruLSVhgpNWbNmDZWVlXr48/jiH8A3IvJB+PXlwGvxDNTSbVx9g9L6ytA3wz9vAGJpwR7pU9imjINWugh8j6O9vwIgTylVJiKnAx+KyLCm54Xb9bwM4HA4ukwgrt4QUn9X20HZZ2cM68uIAb14/eMVfOesYVjM7Y9yf1yRh1EUF6ccbP3gKNR5A/TLcMTdM7IxkSSPGjPMXsGcijxcyooyp4GZhpZvvarXk1O5CkTwG2yU2PqyZ20e7p7nYE/LBWPHrAr4fD6eeOIJJkyYwHnnndchY3Y0zzzzDIsXL2bmzJls3LiRzMzMFo+vcvnZmF8VVaIoGiv3lVPrDXDx8Oy4ynKOlEO4ohpAm9nA4Sp3m9ui1a//6R7g8YNS6k8isgitHEKAHyql1sYzVqtZoMDZSqlfKKU2hh/3o2kwtUY+Wiu1enLQ+r3FREtdBML9364E3m40X69Sqiz8fDWwGxgU6/W6DJEOM37acJoXWFxRy0dL2l0nSmXAzOLqbM5JKiTNFF+o0BsIYjMZ6Z8Zf+izMUk2MwYRQlFCYsMSKglhYKs75cjGcMu3OksGNbae1Fiz8RkdZHn2MqZ8LlWfPgYf3wNL/gQ750PpTvDFL7Xzr3/9iwMHDvDQQw9129R6u93OrFmzKCoqikkxYs2BChS0qeC8wuVj7YFKhmQnRc76jYEkm5l0h4V95XVRj6n/He8rjX5MJBYsWMDw4cPp0aNHXHPTOfaIyJNKqTVKqeeUUn9RSq0NrwO2mVjWAB1h4cH6i58FxPJNthIYKCL9RMQCfBeYE8ukYugicAGwTSmV3+iczHDiDSLSH62n3J5YrneicfqQPEYPyuGfn3yD1xdZST1WPqvMJaAMTE89ENf5Sim8/hDDezs7rEOOUYR0hwWfP3IyziBbFWYJstmV0uI4IYMZtzkVr6MX21zJBB09oKYINn8IS/8Mc3+uaStu/gAKN4Kr/IjH3gKBQIA//OEPjBkzhosv7rCSpU5h9OjRjBw5kjlzWv7XDARDfLWjpM09N7/aUYLRIEw8JaM906RPegKHKzz4AtETsGwmI5sOx14O4fF4WLp0qR7+PP6YGmFbXP9oscR7fgS8FtYEVGhZlze3dpJSKiAid6LJVhiB15RSm0XktvD+l8KtcFahaQyGROQeYCgwkghdBJRSc8PPv0vz5JdJwKPhrgJB4Dal1Ekp2yQi3DLjLH76x3f4cPEGrpt6elzjuIJG5lf1YnxiCT0t8ZUL1HkD9E61k+boWA2+LKeVklovNiKsAxpCDLRVs9kdW7cTk9GAy+enpNZPdnIK2FO0HUqBrw52fQk7vwAUWJ2Q3h/SB2qixEnZ2rZGXt7bb7/Nrl27eP/997ut99eYyy67jN///veUlpaSkRHZUO0qqaXG66e3PfZi8z2lWju+c07JiFsvsJ4+6Q7WHKgkv9JF/4zI2aBOu5ltBTUEgiFMMXRhWr58OR6PRzeAxwkicjtwB9C/UZcygCRgWTxjtvqpDIcTR4mIExClVMy3WGGDNbfJtpcaPS8ksrp8i10ElFI3Rdj2X7SQqQ5w2uBcTh+Sxxtzv+WySSOxW9u+9raguheukDlu788f0L6IOqPbSoq9lXXAhAreLetPTdBEUgyJO1aTgd2ldfRw2o4YLRGwJmoP0Axi0AvF2+HwOsIdmDRR4tT+kDGQUFJPHv/dowwfPpwZM2ZEuVr3YsaMGfzud79j7ty5/OAHP4h4zJKdJdjakPyieYylpCaYGZWb0u459kqxYTYK+0ujG0CzUVMLKajyxNQVZsGCBRiNxg7pjapzTPg38CnwB6BxXXlNvM5Oq7dJItJDRF4F3lZKVYnIUBH5UTwX0zm23DLjLCpqXHywaF2bz/WFDMytyGW4vZwBtranl6O0bjVDezqxdEJP1ASrEbPRQDAUOSQ23F4BwBZXbF6g1WSg1hOgwtVCyFgETDZwZDQSG84FswPK98Dm9/ng2fvYum0HD16ch2H5C7B9HhRthroyiDLXrmbMmDH06tUrahi0tNbL9sLaVpOPGrPmQCVVbj+TB2W2q0l1PSaDgZzUBPaF1SFaYk9p85KXSCxYsIBx48bhdLaup6nT9YSF2fcppb7XpDNZ3JG+WL6ZXkcLY9bnie8A7on3gjrHjpEDe3PGsL7869OVuDxtS2BZUtODyqCVGWnxeX91viA9nFaynB0b+qxHEHokWfFEWQfsb6vBKgE2xRgGRQSzUdgb45fnUZiskJCGcubwu/c3MCgng2umjIWqg7B1Diz/K8z/NXxyr7a2uPUTKNgANYUQjK+spCMxGAxMnz6dzz77DI+nmSokq/dXYBBizuCs9vhZua+cAZkO+rSj5rMpfdITqPYEqHRHv0lJtJrYcLD1IFV1dTUrV67Uw5/HKSIySkTuDD/irveKxQBmKKXeAeqlkerX2HSOA3484ywqa92892XsWcIhpZU+9LNWMyzsSbWFUEihFAzOdiLx1afGREaSlVAosjdgEsWp9iq2tJII0xi72UhprY8aT3yJQ5+s2M66XQU8cMO5GC02sKeGvcSwp2hPg+oC2PEpfPM3+PIx+Ohu+PJxWPcf2L8cynaDNw4j3E5mzJhBXV0dCxcuPGq7VsPXtuSXpTtLUcCkgS2XVbSVvmFjur8senZuos3E/vLW26ItXryYYDCoG8DjkGPdCq1ORNIJ1/CJyJkc3X5MpxszrH9PzhrZn3/PW8VV543GYW/dI/u2NpNCfwL3ZG+Kq0LD5QuQl+5oU8ePeEi2m+uFUCKWkgxLqGBd6SmUByyxlXCIYDQIB8pdDOvVtoJtpRSPvfklfbNTuX5KlBtSo1lLsKlPsgFQIa3c4uA3sO8rtHvSkJZYk9IH0gZAci9I7AEJ6WDonN/peeedh8PhYM6cOUdlru4oqsHjC5HuiO26B8td7Cyu5cx+aR2qQwna3zslwcy+sjpGR1lX1LxUxYFyF4NbUJtYsGABNpuNCRMmdOgcdY4Jx7QV2n1o5QsDRGQZ8AZa/06d44QfX3YW1XUe3vliTavHKqU1vc42uxgXh+RRvUfWpx3abLFiNRlJtJrwB6PUA4a9180xrgMCJFiMHKp042lj0/IvVu/i22353P+9SZjb0ilFDFqSTWIPSM4LayDmgskO5Xu1EOqKl+CL32re4sI/wPq34MA32rqjtzam0ozWsNlsXHTRRcyZM4dQeK2yXrjWYY3t/QRDikU7SnDaTJzep316g9Hom+4gv8LdYrNzoxjYXtTyuvWCBQuYOHEiNlvbivp1ugWd3wqtHqXUahGZDAwOX2S7Uqp9xWU6x5QhfXtwzugB/Ofz1Vw95TSSEqL/029yp7LX6+SWrG3Ek7vg8gXITUvA1sneXz1ZTht7S+owm5rfy/Wx1uIw+NnsTuUcZ1FM44loQdv8SjenZEbONozE7/61kN4ZTm6aFl/JSZNJaOuKJivQSAcvFARvjRYq3bMYDAbN+JkdkNxb8xhTcjVP0ZGhCQu3wYW/7LLL+O9//8uaNWsYO3YsRdVe9pbW0TvGAvYN+Vpf1UtH9oypDCEe+qQnaCoole6GkGhTku1mNhysZHqUtmhFRUVs2rSJG264oVPmqNPpNG6FJmg9pju8FRoAIrIb+GPj8gUR+VgpdWkLp53Y+D3YKndSp7o+gSFWfjzjLGb+9k3enr+aH8+I3ph5TnkeKUYv5yS1XfIopBQKOjTxoTXSHBb2lEReMzOI5gVudqVGi5JGxG4xsb+sjj5pCZhj+CL/av1evtqwj7/ceSnWKAofHYLBCNYk7VFPvYBwdYG2fqjqb4ZVg4QUqWFNxXrD2KRusZ7vfOc7GAwG5syZw9ixY1m5rzxm0VuXL8CKPeX0SUugf0bn/f1zUuwYDcL+MldUA1jfFi1a79Avv/wS0OWPjleatEKDdrRCi+W/1Q+cJyJnAD8JKzv0buWcE5uizfRd/gB9AI/LicucjsucRp05DZclDa8xqUPbm3UEA3OzOO/0Qbw1fw3XThmDM7H5Xf12dzKb3GnckLELs6HtYTWXL0hOqr3T1/4a47SZwsuAKuIX9bCESr6ty6LYb6OHpXmGYySMBiEYUhRWe8hNbT2U+7t/LSQr1cEtl4xr6/TbT72AsN1y9NoiQNAPtcVQsU97Xu8xGi1HDGNKHiRkQEIaGanJTJw4kdmzZ/PAr3/D17tLYxa93ZBfhS8YYtKg2IRu48VkNJCTamdfWR2TiZxkU3/9/WV1EQ3gggULSElJYcyYMZ02T53OQ0SswBAgEc2GTReR6UqpR9s6ViwG0KWUuk5EfgEsEZFraWNT6xOO9AHkn/ZzZMuH2AxBEvxlpLn3NgShg2LCZU47yii6zGkEDZ1TEhArN182gUVrdvDvz1dz25UN3e3whQzMqchjdkUfko1eznfG3LK1AaUUSqmod+WdhclgINVhocYTiBh2HZagrQNucqfSw1IQ87g2s5E9JbX0TrG3mP7/zdaDzF+9i6duvSiuZgOditEMxmSwNUnoCfq1tm5V+RBaqK1DKoAQlw008PNXN7Dqrd8zwJeOgx54jUn4TIl4jYkEDNZmN3fBkGLjoSr6pCe0uVVaPPRNd7B4RwlVbn/U5upaW7RqxvRJa7ZvwYIFnHvuuRiNx+5GTadDmY2WiLka8LZnoFgMoAAopZ4SkdVoNYHNP1UnE/YUarPPoGbPJgJOrd+3IeQnwV9Bgr+84ZHu3kuPuiOq215jIi5zKnXm9Aaj6DYla19Ax4ABvTOYMm4w7y5Yw3VTx5CalMBGVyqvFQ+i0J/A2UmFfD9jNwnGtle5uHxBeqXYSejMEGAUeiTZKK+rjmgAe5ldpBi9bHClMSU5dgNoNhqoDgsSt6Re/rs3F5LmtHPbZW1W+uo6jGbtYWtSAK4Ul0008fNXF/HRR59w3cT+WLyCCufKGQgRFDMuc6r2sKTjMqezsjIJlw/G9rRqWa3t/DwfqnRhMxmjKj8cUYeoY2ROSsRjkqO0RduzZw/79u3jvvvua9ccdbqUHKXURR0xUCzfVg/XP1FKLRCRacDMjrj4iUTIYKbWmkWtNevIRqWwBF1hg1hGgr8ch7+cZE8+hrATHcKIy5yieYxho+gyp+E3dk4W5c3TJ/Dlyh28NncN1sm3sKwmm2yzi1/1WsdIR9tr/kDz/kJK0e8Ye3/1pCRE97xE4KykYj6r7E2J30qmOfYbRqvJwO6SOjKTrBHrGdfuPMzHK7bx2A8vICmha737DkGEgX16MTg3k8/XH+byC8/G28TbExXAFPKS6j5ARt1uDAT4R+FUeppszKz6D4ZqcJtTqDNnUGvNxGXOwGNKwmNy4jU5I3qQjamo85FkM1Prib6+nmI3k2w3s6/MFdUAmowGfMEQ+RVu+jZak6yXP9LX/45rvhaREUqpdsvdtCSIO0QptQ04JCJNg+Uft/fCJwUi+EwOfCYHlfYjylCigtj9lQ0GMcFfTornEFmuI9q+PoO9wRi6LFoo1W1OQUn7PKy87HRGjBnD+wvXkjfEzFU5e7ks9QAWQ/xtuly+ID2Tbe1ueBwviVYTpvC6XaS2W99JOci8yt58UpHHTVk7I4wQGatJ8wIrXX5SE5qH9h6ftRCnw8qdV5xYtWQTR/Xn9bkrqXF7m2UMKzHhN5oabtD2ehLZ6svkxoyd1NmyQYUwhXwkew6R5tqLQQVABIVgIERALNrSgCWDGmsWbnMaXlMiXpMTlyRS5wvwk8kD+PuSPYRCCkOEv6eI0Cctga2F1QRCoagqIxajgQVbi7h5Yr+GdcEFCxbQs2dPhgwZ0sG/NZ3OppFIuwn4oYjsQQuBhtMAWhVpb0ZL31j3AbcAz0TYpwBdQTJOlBi18JElndJG201Bd0P41OEvJ8FXTnbdFgy1WkhSIbhNybgs6draYvjhMzpiSrrZ60nk1eLBlIwZBGtuZ8Tmv3L1qImtntfie1GKUEjRL0qD4mOBiJCZZKW42ktCBCOcbvYyMamIhdU9uTJtH05TjFU8IpgMwr7SOlLzjjaAm/cW8d+vNvPgDeeSEiGh6HjFFwwxfFAewY+/YcWmfUwd37Kh+LyqN1YJMtkZzhoWAwGjjYAxctjYEApgDrnJcO0iu3YzKIUSA4LC5QuQlp7JKTv6cb5ysrOyJ6T0IWhofvPRJyOBDYeqOFzpIS9KzWlGooUtBdUcLHeTl55AKBTiyy+/ZNq0aceFSodOMzq88iCqAVRK3RL+2T3lrLuYYAcUHzclYLRTbexNta1Rkq0KYQtUH2UUE71FZLh2HzlPLA3h07pwRqrLnErIoIUG3SEj75X149PKHJKMfu45NZ/FZ53K50tWc8t3RpOZEr/xcvmC9Ei2kdhF3l89mUk2CqqiZ3lOTzvA4pqezKvqzTXp+2Ie124xUlzjpc4bOMrD/cN/FuOwWbjn6uglJccjRdUeBvfpQWpSAkvW7WrRANYETSyr6cGkpEIcMShuAIQMJkKYmoX4/YEQQWuI0zKdUFvEGM8W0spqcZSZKbf3oThxCJX2PGotmSgxkpuagFGE/WV1UQ2giGAzG/l0UwG3TurPpk2bKCkp0cOfxylhgXZExIYmizQRzRlbCrwYz5gthUCvbGUy77e0/0TmUKWbmfMCXGgczhRHbVxJIzEjBjzmFDzmFMrp37DZGPIdSbjxaT8z63aSrbYA2qfCY3JySHoy33MqgaCda5MquCCjkkRTkH6XnslnK7by5txvuPf6+L4QlFIEQ4p+nVj3FSsp9dmAUQr+eltcjHWUMK8yh+mpB7EZYvubiQgGEfaXuxjaU0sa2Zlfyn++XM+9V08kI7nr33tHoZRiX2kdDquZs0f1Z9HqnfgDwaidbRZW9cKvjFyYcqi9F8btD3JaXgpmqw2wk5iVSFV1KcpsINFXTFrxPhAIiplSx0CKEweTm5zE/jIX5wyMPnS6w8L2whr2lbn09b8ThzeAGo60Pvse8CZwTVsHaum2fXoL+xRw0hpAty9InyTh38Vj+MAVYGryIS5OOUhKrKG1DiBosFBjzabGmn1ko1JYgzUk+MvBU01FrYceoSJ+YpiNwajAD8FCLYuvnzmdq8f35v2vNvCDC0eTkZHe5jm4/UGynFactq5P/7eaDdjNRvwhhdkYObx1WeoBVtVl8mVVT76Tmh/z2AkWI4cqXAzIdGA1GXniP4uxmIzcd237wsfdjUq3H5cviNNuZuKoAXy8dBPrd+Yz9tQ+zY4NKZhf1Yuh9gryrHXtum6dL0hmkpXMpCOJRAkWIzaTEZ8CFU6gAS3bOr1uFz1qN3O+GsJrdWPIzJ+PSu1DlbU3ftPRNyQiQoJF8wIXLFjAwIEDyc3NRee4ZrBSqnHD3YUisj6egVoKgf4wngFPBk7JSuTxs02sX/wB7/vG81FFHp9W5nCus4DpqQfJNMdWcN3hiFBnTObd6uH8t7wvguLq9H1c4txLUqC8IRvV4S8nw7Wbp87w8P43Iea+O4tnZ2RqBf3hhBuXOR2PKXLHEACUIhBUUcVJjzWC0MNpY395XdTuLQPt1QyxV/JJZS4XphzCJLGFsQ0GQaF5/saAnzc+X8vtl51BdlrHC/12JQfKXQ1JROOH9sFiNvHV2t0RDeCaugxKA3a+n7G72b62EAwrhwzJTjoq01YQeqbY2Fd69N8zZDDjtmhVWKcm+6AKCkpKudi9HJSi1tqDosShlCf0ocbak6DBQprDwvZDFSxatJgbb/x+u+ar0y1YKyJnKqVWAISbtHSOInz4ApcAw4CGle14qu5PNPqby7knfTMFPjsfVeTxZVUvFlT14uykYqan7ifXGl22pTPY5k7m1eJB5PsSGesoYWbmTjLMXsBErTFCiUZ2HVef/QWvfr2HW6YkMySpilTPASRcohG5oD+doMGCyx8kI8mK0961a3+NSU+0sK+sZW9kRup+njw8imU1PY4kbsRAgsXIvtI63vl0OSLC/153Tnun263wBIIUVXsa1nJtVjPjh+axdP1ufva985oljcyr7E2aycPYxNJIw8VMnTfAkOykiPWjGYlW9pZG/3v2trhJN3lY6c1jYnpNQ9lR//JFDCgDJUK1rRelCQPI31VOXV0t55+npzScAJwB/EBE6sVK84Ct9VmibckGjaUX6EtAAnAe8ApwNfBtm6d8AtPT4ubWHtu5Km0fcytzWVDViyU12Yx1lDAj7QCn2Ko77doBJezzJvFlVU8WVvci3eThvp4bGJtY1vKJIvhMiVw3fQr/Xb6PR5ab+NXMazCEAtgDFST46ks0ypoV9HuMiVQbU3EaeyElPcGRpbXhOkYF/dGoD8VGa4sGMCqhnDxLLR9V5HFOUmHMDb+NBgOF5VX849PV/PCiMeRmpXTQrLsHheEEosa/t4mjBrB0/R5255dySu6RtmOHfAlscqdxXfoejDF60ZHw+IMk2UzkRElicdrNGERaKIfQ/p7La7MIKMEkNJQdAYgKYfNX0b98CfNXrkcExtR9Ad9WQsapkNIbknqBWVeEOM7okCJ4iM0DPEspNVJENiilfisiz3ASr/+1RLrZy42Zu7g8bR/zKnP4rDKHVQczGWqvYEbqfkYkVLS7RWhFwMJOj5Od7mR2epzs8SbhV0YMhLg09QBXpe2LOcEDICstiRmTR/L+wnXcePF4crJSqLNkUmfJpEEMSSkswbqGpBubpxRnqALr4W/hcPgL0GDWGi0nZmkG0ZGlPTcduwJxs9FAit2MyxfEGqUfqQhclrqf54uGsaYuvfUbhUZ8+OVaAqEQt86YgEJ1qtjvsSSkFPvL6pr1cJ04agAi81mybtdRBvDzyt6YJMR5cbTMq0cphS8QYkxeKsYo/xTGcHlLWa0Xe5QOQ6MdZXxZ3YudHien2o+WKVViwGdKxGdKZPmuSgbkZFIesNG/eBtyeC2aGpyCpJ6QNUTTXnT2Bkem1jdVp1tSnw3aEcRiAN3hny4R6QWUAf06agInIknGAFen7+PS1IMsqOrFJxW5/OHwaPpZq5mReoBxiSUxeR4BJRzwOtjpSWZH2OCVBLSaM5OE6GetYWryIQbZqhlkryI1FtHXCNx48Xhmf7WR1z9ewUM3R7i5CnuLPlMilbZcqs0BxvdLw24VcJVBXbHWdLm2GEq2Q0Gj9WhbimYIE3uEjWIPTc2gk+qwejhtbC+siWoAAc5MKuHtMjdzKvpwuqMspql4vH4+WbaJKeMGU+KBzYerGdQjCUsnyf4cKxSKQxVuPP5QMwHbtGQHQ/v1ZMn63fxwulbs7woa+ao6mwmJxSS3I+mr1hsgLz0hai/Peno4bRRVR19TH2avwEiIdXXpzQxgPR6vn427D3PN+aMp9wplyklGctigKwX+Oti7BHYv1LYZLZDWD7KGas3Cnb2OVuDQ6RJEZI1SqsUO5rEc05hYDODHIpIC/BFYg5YB+kosg4vIRcBfACPwilLqiSb7h6BpO40BHlRKPR3enouW6poNhICXlVJ/Ce97BK1Av95BeUApNTe871doasFB4H+UUvNimWdnYTMEuST1IBcm57OkJpuPKvL4c+FwepnrmJ56gInOoqMSMaoDZnZ4nJrB8ySzx5OET2lf5KlGL4PsVUyzHWKgrYp+1pq4FBsikZGSyFXnjeLt+Wv4wXfGk5cdvdWr2x8i1WEOtx8TSMrWHvUoBb5aqC3SDGJdsfa8dMeRY0w2zRA2NowdpHaekmBpVRrTKIpLUw/yj5JBbPckMyTKF2djlqzfjcvrZ8Y5I3HaTBRUuimp8TKslzNqq7TujscfZGthNUXVXhyWyL/7c0YP4KX3l1JSUUNmahJLarLxKBPTUmLPom2KPxDCbDTEpLeYkmDWvnGilLckGIMMslex3pXG99gTcYwNuw7hDwQZO7QPVpOBHUU1pDssWrhXRNNNtDSaS9APVYe0mzkECGkSUil5YS+xl/bZTcgAY/dZAz8JOFVENrSwX4DkFvY3IxZB3MfCT/8rIh8DNqVUq98YImIEXgCmAvnAShGZo1S4UE2jHE1d/vImpweA+5RSa0QkCVgtIvMbnftsvbFsdL2hwHfRknV6AV+IyCClVCcW6cWG2aA4P7mAc50FfFObyeyKPvyt+FTeK+/HZGchJX4bOz1OCv3aWoiREH2ttZyffJhBtmoG2qpIN3k7VWHphovG88Gi9bz20QoeueU7kQ9SCn8wxCmZydG/8EWOaNaln3Jke8ALdSVhoxg2jofXQihcQC3GRiHUeuOYpRnLNpBkM2FsYd2onsnOAv5b3pfZ5X0Y0rul/ymNeSu2kpWayOhBOYgIiTYzvkCINQcq6Z1iY1CPJKxtUYLvQhSKwioPWw5XEwKSbaaoHvk5o0/hpfeXsnT9Hi6fPIrPK3szwFrNAFvLiuvRL67V/I3OTYlJa9FmMpJkN+ELhLBE+f2OSijnrbIBVAQsEaMgq7YewGQ0MGpgjtbezhOgpKUm50Yz2FO1R3jOBH1Qvg+KtnBEDCd8A5jWD1L7hW/qMqPqLeq0m1j617Xp+z6WJBgjcAnQt/54EUEp9adWTh0P7FJK7Qmf8xaacm+DAVRKFQPF4SxTGm0vAArCz2tEZCuaBmFj49mUGcBbSikvsFdEdoXnsLy193isMAhMSCrhzMQS1rvSmF3Rh/fL+5Js9DLIVs35yYcZaKumv7WmXb054yHNmcBV543mP5+v5q5rJ5MeocDbEwiRYjeTGo/kjckKyTnaox4V0mR5GrzFIk3UtbBRj1tb8tHh08SsFr9gDCJktLJuBGA1hJiWks+7Zf3Z73XQp4VatooaF99s2sv3Lhx7lFG1mAyYjUJhlYeSGh/DejnJcnZvb7Cp19eacnvfnmn0zkxh6brdDBh3Lof9Du7o0dK/YcvU1/xlOWNfG8522tlZXBPVAI52lPFW2QA2uNIiZvau2naAYf17NshVWU0GdhbVkpFobVHqqgER7fNrsnKUEI4Kgd8F+atg37JwElhYiDg5D9L7a593R5a2rmjqfKmoE5mOXPurJxb//SPAA2xEC0fGSm/gYKPX+Wjpq21CRPoCpwHfNNp8p4j8AFiF5ilWhK+3osn1mgn3isitwK0AFkvXfCBFYLSjnNGOcmqDJhyGQLe4Ybz4rGHMmreKRat3cNX5px29M5y0MDKnBe+vrYhB8/ocGdBjWMN1tBBqOHQaRwg1K8lKUQtt0eq5MPkQH5Xn8VFFHndmb4163IKV2wmGFNMmDG3+FsLeoD8QYu3BSrKdVob0dGLrZt5gY69P0bLX1xgR4ZzRA/jvwnVkFKfhNPo4M7E4rjlEq/lrjfRECztbuGSepY4Uo5d1dc0NYHWdh+37i7h5+pGG5VazkWq3n5IaLz2c7cgAFUPz8Clo3mL1YSjbeUQeSoU0I5jaV/MWk8KfW3vqCZlw09IyVqNjBG2J7DuAC7hJKbXmWM4zFgOYE0+XbSKvxLRp0UpEEoH/AvcopeprCV4EHguP9Rhas+6bY72eUupl4GUAh8PR5cK+iTH2UDwW9O+dQb9e6XyxcnszA+gJhHDaTKRFUEXoUI4KoQ44sj3mEGoP0m0ZpPitYMkkaIzuaSQaA5yffJjPKnO4Ln1v1AYG81ZsZWBuJgN6Z0Qdy2wykGwUSmp8lNWVMrSnk+xkW7fwBtvq9TXlnNEDeGv+alZsPsD1Z4biXnuu8wYYHKXmryWSrCbMRkNUtQ8RGOUoZ1VtBkElR5VmrN1+EKVg7Kl5R51jMxvYWVxDZlKMXmBbMFrAbtFKg+pRCvxuLbqRv1KbtFJgMGmZp2n9IKWPFkJ1ZGpGtTvcFcdPa8tYABcDA8OPM9C+24+psGYsn8RPReRCpdTnbRw7H2jccygHiDlvWkTMaMZvVuO+o0qpokbH/J0j0kztup6OxgXjBvPKnK8bkh7q8QVCDO+d3HVd9GMOoe7E4t/AeIAqrR9qnTmsnmFJp86cfpR6xndS8plXmcPHFbn8MIJU0sGiCjbvKeDOaya1PkcREm0m/MEQ6/OrKKjyMLSnM6JQ77EgXq+vKSNO6Y01wYF71zdcMC0+b6W+5i83Ss1fS4gIPZKsHK7yRJXcGpVQxuLqnuz2JDHIfqTudtXWA9itZob263nU8RaT5gUWV3vITj4Gah4iYEnQHo0JBcFTqYVQgwuPhFEtDi2MmtYfksOlGcdRGDXGZawZwBtKKQWsEJEUEekZPjcqYc/xBqC/UupREckDspVSba5Pj8UArgA+EBED4OeI9pKz5dNYCQwUkX7AIbQEletjmVT4Db4KbG261tjkF3QFsCn8fA7wbxH5E1oSzED0gv02c97YQfx99tcsXL2Tay/Qsom9/iCJVhMZ8az9dSYthFAPHtyLt7KQ1FAFjnAxfz0Bg7Wh3VumOY3rE828Xz2CqyJIJc1bsRURWpUFaozZaCDZJpTX+Vi6q5RTezrplWw7pjcPHn+QrQXVFNd4SYjD62tMUEzY+o/Dvftbkg2jgbaNVV/zd1peStSav9bIdNrIr3RH3T8ioQJBsd6V3swAjh6YE7Ght81sZEdxLVlOW8d7gbFiMGpr2tZGX6f1STdVB6GkvgFFOBvVkREOo/bXEnASw2HUYz9/k4isavT65XB0rRlRlrEg8jJZb8KGswX+ihZWPR94FK0x9n+BcbFOvp5YDOAzwARgY9hSx4RSKiAidwLz0MogXlNKbRaR28L7XxKRbLR1PCcQEpF7gKHASOBGYKOIrAsPWV/u8JSIjEYLb+4DfhIeb7OIvIN2hxEAftodMkCPN/r2TOeUnEwWrNx+xAAGQgzt5Tw+NNTCIVRb9mC2+rJICneHMYR8JPg1Y+jwlZHgLyerbitGFeQxvuJhi5HSwgxM9iRNUsqSTq0plXkrtnL6kLyjvOFY5+GwmggEQ2w8VEVhlZtTezrbHP5rK029PmecXl9jvq7JwjzgTKo2LWLj7sOcNiin9ZMaUesNkJeWQIo9/huoerWPaF1+Eo0BBtqqWVeXxjXp2s1OSUUN+wvLmX7OiIhjWsKCx0XVHnoeCy8wVhon3dRnooJmGANuKNykJd5oB2tZq8k52pJBcl44jJrV2R1uAkqpsa0dFGUZq2F3hFNisTFnKKXGiMhaAKVUhYjE9eGK5b9xJ7CpLcavnrDBmttk20uNnheihSqbspQo1VxKqRtbuN7jwONtnafO0UwZN5i/fbCUovJqUp0OHFYjGYnHrqNLR5DcRB4pZLBQa+1BrbXHkYPCWosOfxmbyg2kB0s4w7OLTNcuAFbkBzhU4uJXkxPIq/yWOoumteg2Jcfc9s0U9gYrXH6+3l3G4OwkclLsnXIz0ZFeXz1KwbyqHE4ZnEGVycjSdbvaZADra/4GZLWvabrZaCA1wUKtNxA1pDzKUca7Zf2pCphJNvlZtU1zLpqu/zXGbjayo6iGLKctbu/0mCEC5gTt0ZhQAOpKoWI/qIC2Hq6CWp1iSl8tGzWpp+Yt2lKOWdJNtGWsRsS7bOUPVyeo8HUyaVuCZgOxGMACYJGIfIomPw8QSxmEznHK+WMH8bcPlrJw1U6+c85IRuemdF2IKE4sRgNJNjO+QDBq+nxjrcU6g5N7D57Ojc6dXObcRYK/jL8uWIbN5ObKISZ61GzA0NAk3IjbfEQ1o65Rk/DI19G8wWAwxJbD1RRWeRjUIwmb2YDZaGj377YzvL56dniS2e9N4sc9DsGQXJas282d10yOzYArhcsf5LTclA7pmNMz2caWguroBjChnHfL+rPBlcY5ziJWbz1ASqKdU3IyIx4PWvJStdtPYZWH3indyAtsCwaTVi5ka1QDrhQEPFC8BQ6vPpJ0YzSDM+wtZg6BHs0zmzuClpaxGjEHLaP/LbTkl6rW1v/CPAd8AGSJyONo/akfimeesRjAveGHJfzQOcHJ7ZHKoLws5q/cxtVTTjtKp+14IttpZWdxbXQD2IhTbNUMtVc0SCV5TNnMWVfC2aMHsq/vdParIHZ/JQ5/GQk+TVIqzb2PHnXbG8bwGBMbEm1c5jTqLOl4jUfavhmNBpw2odrt55u9Wg9SpcBkECwmAxaTAavJiNVkwGY2YDEZMRsFk0GrNzQZtZ9GgzRkl9Z7fUU1HhwWU4d4fY2ZV9kbh8HP2UlFVI0awNOzFrC/oJy+vVrXj6zzBclqY81fS6Q6LC12helnrcFp9LHelcbEpEJWbd3PmCG5LTZEALBbjOwsqiHbaYuYZXpcIqLVI5qbGPVQQGtfWLEPDn4LFz8R8fQO4GwiLGOhKTfURwLnopVA7EIrg4hJgk8pNUtEVgNT0CKFlyulotcxtUCLBjDsZg5USukiWicZU8YN5sX/LsFK4Ljz/upJbWPJxvTUAw1SSaa9y6iq9TDtTO0OWYkRl0VbG6S+R0CjJuGaYdR+prr3N8TvA2LWJKUs6Q0eo8GcSshwJEQbUlpDao8/iMsXJBRSBButONSP1XgNwmLUDKbLpy1zJ9vMHZ4IUR6w8G1tJtNS8rEZQkwMG8Al63e3agC1mj/V5pq/lkiwGLGZjQRCClME0WODwMiEcta70thfWEFxRW2L4c96zEYD1T4/BVVuclLbnqV6XFHvLZrtWhF/J6GUirqM1egYBfw0zvG3AdtaPbAVWjSASqmgiGSKiEUpFV+nZZ3jkkmjT+HF/y5h0ZqdjB/Ys/UTuiFJdjMiQkipmIz4qIRy+lhqmFOeR8KKV0hOtHHG8L7RT2jcJNx+5IvWEAo0KGfUe4wZdTvJVlqGqQI8puRGOouacfQZY6v9UmGjGQgqbGYDxk5a0/myqhchhKnJ2rJMVloSQ/r0YMm6Xdx48fgWz4235q8lGovkJkbxdEcnlLG0JpvPNmjVUmOHtG4AQTOuO4tryU62YToBC9NPJETkn8DdSqnK8OtU4Bml1M1tHSuWT+c+YJmIzAEa+kXpa4AnNqkpSYwa0JN3F23kF9fFUAPXDTGKkO6wUOXyY4vS7LkxIjA97QDP7e9H4brdTJ84LGL6fGuEDCZqrc0FiK3BmgYvMcFfjsNXenR5hlhwWY6sLWpixI28xYZ5CkahU8N1ASUsqOrFqIRysi1Hyg8mjh7Aq3O+pryqjrQI7fJAC8smxlnz1xrprYjk1pdDfLFym9bGLUbdRpPRgMvn53Clm7y0yO9Lp9swst74QUMW6GktHB+VWAzg4fDDAOiaICcB/mAIq8nA9ReM4pd/+4w9h8vp3yu6QkR3JivJSmmtFxuxGbIzE0v46679+P0Bpp3RgQkCInhNTrwmJxX0bdjcuDwjwVceLs/YiTHcMKO93mK8fFObSWXQyrSUo6NME0cN4JXZX/P1xj1cOrF5eUFDzV9u/DV/LZFsa1kk12ny09O1i+X79nDblRPblG2bYDGyq6SOXil23Qvs3hhEJDXcAhMRSSM2W9aMWNQgfhu+SJL2UtXGcyGd4we3L8iwXk76TR7BL//2Ge8u3sgvvze5q6cVF21dBzSKQrbPx5SSjaHnqUDrUkntIXJ5huYt1tcr1odSI3mLDYYx/GjqLcbL55U5ZJtdjEwoP2r7wNxMeqQl8dW63RENYEPNXye1zDMaWhfJZdsXAEwcF7n+LxomowG338+hCjd90nUvsBvzDLBcRN4Nv74G+H08A8WiBjEceJNwG3QRKQV+oJTaHM8Fdbo3gWAIi8kQXgtJ4IxTc3l74fFrABOsxnAfyVBMa2UllbXs3bWTzLOv5ePKPgxztC6V1OE08hbLG2lPR/IWG68tgtb6zWVObchEdZnT8JiSYq5bBNjrSWSHJ5kbM3Y2E26ub4790dJNeLx+bNYjBrejav5aoyWRXKUUe9atwJo3giJbX/rTtsbddrORXSW19EqxxyTXpHPsUUq9Ee5Cc35405VNeozGTCxu48vAvUqphQAici7wd+CseC6o030JhRS13gDDeyU3hICuO3cE9744l535pQzMid4MursiCFlJVgqqPDisrX+hzf9mG0rBZRMGMd+V3qpU0rEkurdYq3mKYY/R4S8n1X0AaahbNOEypzYYxPpHwBi5U8i8qhysEogoLQRaGPS9L9exausBJo4+0rDc7Q8yMie5Q2r+WqIlkdytewspKimj53e+y/q6NCYktc0Aal5gkPwKN/0ydC+wOyIiVmA0WgcxE3B1WKLv0baOFYsBdNQbPwCl1CIR0T8ZJxJKUecLElKKU7IS6dWoIPjqycO598W5vLt4Ew/ccG7XzbEdZCZZOVQRvY9kY+at2MLQftlce4qfJXsDrUoldTkieE1JeE1JVNj7NGw2hALYAxUNnmKkukWv0dHEKKZSbMjg65osJjsLcURRKjltcC4Ou4Ul63Y1GMD6yEFUkdkOxGYykmgz4Qs2F8n9bMUWLGYTE04bwnpXGiFFMy+2NRIsRvaU1tI71d7pxlwnLmajrU2splFzlniIxQDuEZFfo4VBAb6PVhivcwLgCwTx+ENkJlkZnJ2Eo8m6Sm5WCmcNy+PthRuOWwOYbDeDRO8jWc/uQ6XsPFjCz753HonGABckH+aTylyuTd9LVhSppO5KyGCizpJJnaVRFxSlMIfcYU+xoiEbNdlzCEO4k1QQAx+b52InAaqdDeuLjZNuzCYjZw7vx9L1exqSUdz+IAOzko5ZIXnP5OYiuYFAkC++3c45o/ozKbOW1QVWltX04BxnUQsjNcdoMBAMBcmvcNE/o3PDuTpxkaOUuqgjBorFAN4M/BZ4H62w8StirNg/kQlYnPjFgtMbDhMpRchgxm+w4TfaUdK9BFGbEgwp6rwB7BYjp/dJJT3RErVg+brzRnL38x+z7UAxQ/KyIh7TnbGajPRw2iip8UaV0wH4fMVWjAZhyjhN+eHi1IN8WpnDJ1Gkko47RPAbE6iyJ1BlP9KCUVQIW6ASu6+cZaV2hhkO0CewD1vVERHihoL+cCh16tB0Fqzczpa9BQzt3xMBeqV0vvdXTySR3BWb9lFZ6+aiCUMZ6yilj7WG98r7MSGpGJOoyANFIcFiYk9JHTmpCSekF6jQxImP03f2tYiMUEptbO9AsWSBVgD/094LnWh4Ugfxfs976e8MYfeXk+CrwOk9jNN7mCRvEQYVRCGIChE0mAkY7PiNti43jCoc7gQY1COJ3LTWU76vOmcY97zwCe8u3sSvbzy/xWO7K/0yHBRWeaK20QqFFJ9/s5Xxw/qS5tTq19JMPiY5C1lY3ZMr0/aR3EQq6URBiQG3OY0l3kE84xvBz3puJJRYijHkw+6vaAihJvjLyXDtwaS2cWsPxUMG2LF0Dlck5oEjE2tZjqZZl5CudRzpRJKsJkwGOUokd96KLaQk2jljWF8MAtel7+Gpw6NYWNWTqSltkwY1GrQGCgfLXQzIPL68QIUiGFJ4AyF84YcnEMTlDeLyB3H7Avi9XlLMfuIqnut6JgI3ichetBBovURfm4Xbo35KReQfRJemUEqpH7X1YiccjdZfKu19OMxobXujome7v4Jk72GcnsMk+ooxqAAgiFIEDBYCDR5j59+LefxBfIEQPZNtDOyRhD1GodbemclMHN6HtxduOG4NoNNmJjPJSkWdj4QIXuD6nfkUlddw+1XnHLX90tQDLKruybyqHK5NP7Ej//Mqe5Nu8nC6Q+tTGoySdKO1f6tg/ID5fLTdza+nukiq3AgV68IHCSSkaWoE9UKujkxNIb2DPuciQrbT1iCSW+vysmTdbi6bNAJTOCw6OqGcwbZK3i/vy2RnIRZD2wQDEiwm9pbWkZNqxxpHQ4TOQjNu2v+yLxjC6w9R5wvg8gVx+4N4fMGGVnrCkXwho0HrI1vfe9YbiEtAoTtwcUcN1NJt2scRtuUB90CMVcUnOJ5A6Kg70AaaFD0frr/PUiGsgRoSwnfVTu9hkj2HSfQWIw1qHoLfYMNntBOSjunvGAiGcPm07hyjc1PaXBsHcN15I7jzuY/YvLeIYf16tH5CN6R/poNvarwRvcDPVmwlwWpm0uhTjtrey+JmnKOEzyt7Mz31AHbDiSkxeciXwCZ3Gtel78bYUriwUfu38aeP5s9vLeQLywVcObYvuMs1WZ66kvCjGEqPJN1gMGneYb1BTMjQtOssSXF9zhuL5C5cvQNfIMhFZx5pXiAC383Yw2/zxzCvqjfTUw9GGyoi9V7ggXIXA7OOfQ8QhcLrD+HyBaj1Bql0+ahy+3H7tc9gvXFDgSFs3IwGwWY2ttoAXKnjs78vgFJqf7j92UCgcdx9f1vHimoAlVL/rX8uIv3ROnlPAp5Ak7k4qRmcncSkgZms3FeOABlJ1tbrhsSA15yM15xMBX05hCY4K2GlgQR/GYneYlLdB0j2HMISKkNp2RsEDVZ8RjsBgy3mLwultHU+gwhDeznpnWKPu7H1VZOGc9f/fcy7izcetwYw2W4mzWGhxuM/qoja6w+wcNUOJo8ZeFRdWz3T0w7w7cEsvqzqxSVt/BI9Xvi8sjcmCXG+MxY1Go2Jowbw57cWsnnnQa4cP+CIYePUIwcF/Zr6QF3xEeNYsR+KNh05xmgNn9vEY2yqZNCExiK581ZsJbdHKqf2yz7qmCH2KkYllDGnvA9TnIdJMLbtBsYR9gKDIUWCxYQ1rNphMRmwGA2YjNIhzb6DIYXLH8DlDVLt8VPp8lPt8RMMHbkZMRkEs9FAkrXj5K6OR0Tkx8DdaPqB64AzgeUcqQuMmdbUIE4FHkSTs/8jcJtSKnJu9ElGks3MteNymTq0B8v3lPHVjhICIUW6wxJVrywajZUGSh2D2AeNQk1lOLylpHjySfEcJMlXSP29n8KIz5jQPOlGKdz+IP6gIi8tgf6ZjnaHcLLTkpg8qi9vL9rIb2ZOOT7U4ZsgCAMyE1m5rxx7Iy/w6w17qHV7mXbmqRHPO8VWo0klVeRyYXI+ZkObtaG7Na6gka+qs5mQWISzDeucmalJ9O+dwYJVO/l1tAxhoxmSsrVHY/xuzSC6SsKGsRhKtkHBuiPHWBKbG8WEDG1Mjojk7i4oZ832g9wy46yIn8vr0vfwwMFxfFKZyzXp+2J+f6B5VjaTkfwKN6GQAmnkeQEGEWxmA3azEbvFSILFhM1sbFDrsJgMmA3SMC+F1iquzhfE5Q1Q4fJT7fZT5zvytSrh92Yzte7JnaTcDYwDViilzhORIWiJmm2mpTXAd4GxwNPAz4Ag4Gz4QypVHu3ck4lUh4XvjOjJ5EGZrN5fwRdbiyit9ZGSYCaxhYzDVjlKaaAPhzgdIJyYUI7DV4bTe5gU90Gc3gIMSruzDQQVNSELSY5kBvdNx2nrmNZYANedO5Lb/zybTXuLGNE/u/UTuiGpDjMpdjMuf7DhRmXeiq2kJzs4vQXpnBmp+/nD4dHMq8rh0hPMC1xSk41HmZiWcqhN53kCQS6dMITn319GWZWL9OQ2NL822yElV3vUoxT4ahuFUMPG8fBaTceuHntqgzHsTzKvfq1lq047M3Lv1n62Ws5ILGZuRS7Tkg+1yciDJpprNkWO7iilJZzUeAJUusIeW9hm1RtKQctEtpgNuH1B/MEja29G3auLB49SyiMiiIhVKbVNRAbHM1BL39Dj0P5+PwfuC29rLE3WP54Lnqg4rCYmDcrkjP5pbMyvYt6WQg5WuEiymkgOy/J0BFpiQja11myKkoZpG1UIk7cCb2UB2aESzs2opFeoEIO3MFwmKmBOAGsiGOPv0XjlOcP46XNzeGfRxuPWAArCgKxEVu+vwGYyUF3n4esNe7j6/NNabJU2IqGCYfYKZpWewk63k5lZO0kzHf8KYSGlhT9PsVUxwFYT83n1a9/fPXc4z723lLnfbOfGC9uZUygC1iTtkdbo60WFwF3ZxDCWQOlO0lSIr76pY2KekYuCX+AqSzuqXMNndIAI16Tv5dvaTGZX9OHGzF3tm+dRUxZMRqGlAItSilBIW8+zGA1tjhDpNCNfRFKAD4H5IlKBJtjQZlpaA+wbz4AnO1aTkbF90zgtL5XthTV8vrmQAxUu7GYj6Q5LhxjCkFK4fUFqvYGGO06jOJg0agLnDcnCbjFqd9PuCqgphOpDULIdKvZqoSfQbmEsCVqYyRSbYndWaiLnje7P24s28OgPLzguw6AA6Q4LSTYT3kCIL1fvIBAMMW1C5PBnPSJwf+/1zK3I5b3yvmzan8r1Gbs5z1nQ5k4j3YVDvgT+VjSEw34Hd6W3rbWv2xegb4aDARmZ9ExPYs7XW9tvAKMhhnBWaRpkNrrRDwVYs2ELW0vf4vdXDsBntOL0FpDpOmLgAmLBZU6lnyWdRxwHmVs9jCqnkGw9dmFsEcFoFD1zsINQSl0RfvqIiCwEkoHP4hmrU4t1ROQi4C9oWaOvKKWeaLJ/CPAPYAzwoFLq6fD2XOANIBsIAS8rpf4S3vdHYDrgA3YDP1RKVYpIX2ArUJ92tkIpdVtnvr+WMBq0xJNTeyaxt7SOL7cVs7WgBrNRyEi0xtwxQymtnqfWG8DrD2lZ5Aqyk22MyEmmX4aDbKeN9KZjihz50ugxFAZO1YyipxJqiqAqH8p2Qvkebf2lPmBjTtDuwKN4itedO5Jb//QB63cXMPqUXu39NXUJIsIpWYmsPVDJvBVb6dcrnUG5rRf4m0RxWdoBxiWW8ErxYF4pHsLSmmxuydpGL0tsrda6AwElzKnI44PyvtgkyO09tjAhMfaemUppRdQ5qQkYDAYuO+tUZn2xDq8vgLUDBXBbxWDiX8sOYjYZGXzmFLalaZmaxpC3IdO6cdPwmWoLMy3vQzF4jIlHWsBZNI/RY0o5JuVIOh2HUmpxe87vtE+riBiBF4CpQD6wUkTmNOnaXY5WZH95k9MDwH1KqTVhGabVIjI/fO584FdKqYCIPAn8Cvhl+LzdSqnRnfWe4kFE6J+ZSP/MRAqq3Hy1o4RV+yu0zNHE5pmj/qBm7Fy+YEOqS0qCheG9kzkl00F2sp0spzW+pBYRbf3EngpZQ2DgBWGjWAW1RVB9GEp2QNmOI56i0QJWJ5i07NMrzhnK7X+ezTuLNh63BhC0331VTS3rdx5qs25cT4ubh3qvY1F1NrNKT+H+A+O4Im0/01MPtLnjyLFmtyeJvxUN4aAvkQmJRfwgcycpbVwTc/uCZCfbGupIL5twKn/76FsWrtvDReMHdca0IxIIBvnPl+uZNn4QiQlHsuGDBis11mxqrI3C9OGkshWlVlwuF9+1rSM9UEKK5yCGcEpLCANuc0qzpuE+Y4K+PtcNEJGlSqmJIlLD0TXq9YXwzraO2Zm3a+OBXUqpPQAi8hYwA2gwgEqpYqBYRC5pfKJSqgAoCD+vEZGtQG9gi1Lq80aHrgCu7sT30KH0TLZz3bg8pg7NZvnuUpbsLCUQVBgMmh0CsJoN9MtwMDArkd6pCfRw2tqXTNMaIlqBsj1FCy8NOA9CIW2NpSofSrZC0RaozgeEDDEwZXQf3lm0kcd/dOFxGwY1iLBm0x4App7RcvgzEiJwXnIhpznK+WfJKbxT1p/lNVnc2mM7p9iqO3q67cYbMvBuWT/mVuaSYvRxX88NjE0sa/tAShEIqaP08s4f0x+HzcKcr7ceUwP4xerdFFXUctOFp7Uokgs0JJX1zzBzz77zWROYwP/03NKoBKk8/KhoHkY1WI/SXKxfY+wo7UWd2FBKTQz/7LCizJi+WcPeXI/GxyulDrRyWm+gcbpcPnBGWycYDm2eBnwTYffNwNuNXvcTkbVANfCQUmpJhPFuBW4FsFg6R7SzNdIcFi4Z2YvJg7NYva8CTyBIbmoC2ck2UhM6LmEmbgwGSOqhPXJOPxI6rTwIZbu47rx8fvTMh6xZu5bTB2QdSVw4jsJHSilmL93MqIG9yUiOv9VVisnH3T23MLG2iNdKBvHwwTFMS8nnuvS92LpJ0fxmVwovFw+h2G9nivMQ12fsbnM9XD3eQIhkuxmn/chXh81iZtq4gcz5eisv3H3ZMfv8vjl/LalJdi49cwhbi2oor/NFF8kNk2Lyc3FKPh9W9OUy7wH6WmsbSpAaYwp6mrWBy6rbgTGsvaio11482jB6TE7dWzyOiEUQ9y7gN0ARNLQrUUBrfdcifQraFB8SkUTgv8A9SqnqJvseRAuVzgpvKgDylFJlInI68KGIDGt6nlLqZTSNQxwOR5fGqxKtJiYPzmz9wK6mcei050guf+BcfvKXj3lnXwqnX3opFG6Eyv1Hei5ZHGBzdno/yPawZudhth0o4anbL8HtD0ZNc4+V0xPLONX+LW+X9WdeZQ6rajO5OWs7pzm6rlqoLmhiVukAFlb3oofZxUO91zIsobJdY/qCIU7t6WxW/H3ZWafy/pLNLFy7h/PHDIhydsdR4/LywdItzLzwNKwWE9nJdoqrY1PGuTT1IPOrevNOaT9+0TtyP+WA0UaNsSc1tp5HNjZqceho8Bg1man634amvdjcW4ymvagTOyJyb0v7lVJ/auuYsXxD3Q0MVkq1NV6SDzQq8iGHNqSqiogZzfjNUkq932TfTOBSYIpSWvBQKeUlnPSvlFotIruBQcCqNs5bpxXS0tKYOnUq73z8BU889zIy+GIIeLVs04oDULQRSndptVsqpK0fWp1atmk3uTv+1/x1WMxGfnjhaWwqrMUfDLVbATzBGOSHWTs5O6mIl4uH8NThUZwVXmc71o20V9Zm8FrxIKqCFi5NPcDVaXuxtrEXZlMCQS2NPyOpedbw1ZOG85vXv+DO/5vD2r/d1enJMB8s3Yzb6+f7U7XM0/quMNGanTfGYQxwWeoB/lM2gO3uZAbbq2K7aJMWh/UYQgHsjeSlNKO4lx512xqOaay9WGdOw2VJx21KPq6iJt2A+tDnYLQyvTnh19PRVIraTCyf0oNo4oNtZSUwUET6AYeA7wLXx3KiaDGUV4GtTa16OLP0l8BkpZSr0fZMoFwpFQy3bhsI7Ilj3joxcN1113HTTTexcuVKxo8frxm3tP7aY8C5EApqiTVVh6B4q7aWWF2C1kojXO9lSeySL4D65IlLzhhMRnICA4KwpaAas71j5jLIXs0fclcyu6IPH5b3YYMrje9n7mJSUmGn2//KgJnXSwbxTW0WeZZaft5rY5vq+1rC7QsyODsJY4Q34bBb+OvdM7jkgX/yx7e/4qFObpr+/+2deXhU1fnHP++dmcxkDwk7iKCygyCLiiJKRYoL4oKg4varraXaKlar1daKWvcq2iK41X2NiCCuKG5VUZHVfSMIYQsEsmf29/fHvQkhhiSTTDY4n+e5z8zce865Z84k8857znve75NvreKALpkcMdBOXuDzVIjkKgnuugd5fEYurxd057n8A/hHt5WN+lyilptSbwdKvdW1F8tIDu6osr64g3R/bpWgGxdlngzKPFmUJTiG0ZNlvMU9oKo3AIjIYmCYqhY7r2cCLzSkzXoJ4gLvicirVFHfrcvddKI0/wi8ib0N4hFV/UpEpjvX7xeRztgeWhoQFZEZwADs6dVzgS9EZJXT5LWq+howG/Bib4CEXdsdxgA3ikgYO2vNdJOtpumYNGkSHo+H7Oxs2wBWx3JBWlf72G/krmjToo2Q/xPkfQ2FG5xJ8Sh4kmvdfhFPlqywgyfOGTcUgC4ZPn7IKyEcieKOk/abx1ImZ63j8JQ8Hsrrx/1b+/NhUSd+2+k7OjWBuK4q/K+4E09s601AXUzJWhvXqNRo1BYT7pKx5/ycJxzelynHDOafT73HlGMOps9+7eNy7+ps3FbIkhU/cd25Y3dbb9wlklv3Z+izopya+TOPbuvDmrJMhsR7qlqEkCuZgsRkCnbTXtwVdGNPo+aT4d9Ax7Jd2otBK2k3g1jmyaTck97iUmqtiB7Y2+AqCEIVlzwGRLX2fxARub6m8xXWuC2TnJyspaWlLd2NNsvEiRNZvXo1P//8c8MCH8IBe+tFwQbI+wq2f2+fA7A89rSpJzHu06bn3foCi5Z+w+Z51+BLsKfOfs4v5bstxaQmxj+yL6qwpLArz+YfSESFyVk5jEndQoorXLvyQj3ZFvLxcF5f1pRl0sdXyEWdvqVbQlndFWOgxB9mv8xE+nWuPdJ8y45i+p0/i2G9u7LkrgubJCDmX8//j7888DrfP/FnenffZWQLy0N8mpNPaj3T/4VV+PO6w0hxhbl5v89bdHbeHSmvNIgVa4yJoZ1YTthF1S0apZUeYxZhV+0Jw/eEhoNY4TIOv+ThBtUXkTJVTa67ZPxx4j+mAC9h/4Q+FchW1Vtibas+grht3tAZmoapU6fyyiuv8MknnzBq1KjYG3B7IbOXfRwwxt5+Ubbd2X7xnT1tWpiLnWrF2aiv2OuKYtlBNi5PlUfneS3fZKXlQeb/7yvOPnZIpfED6JqRyI/bSohEo7WmRGsIlsBxGZsYlpzPo9t688z2g3hmuy27lGyFSHVVOaq/do6UikcrVJl5JqqwuLAbz22304Zd0OF7jkvfGP/MNKqoKt3b1Z3rs3NmKrdf9Gumz1rIE4tXcv6vh8W5M3b052H999vN+AGk+n4pklsbblEmZ+Uwd+sAPivpwGGp2+Le1/oSdiVS6OpGoa9b5TnRKL5wwW7TqOmBTbtt0QhaiZQlZFV6i6WeTPyevXtDv6reLCKvAxXinf+nqisb0lZ9okDfpYboTVVtm8qohrhx8skn4/V6yc7ObpgBrI5lQUpH++jmfHEGyyBYCuFye3N+qMx+9BfbWzMCRfbUarDEfh6q8HxklyFUBaIgHhZ+tJZSf5Bzxg3Z7dYel73/8se8ElJ9TfPlkeUJcEWXL/mqPINNwWSKIh6Kqxz5YS/rIikURzyEtObpLkFJtkKkuUJEEbaEkhiSlM+FHb+jg6d+UZCxUh6K0D7VW+/9qL87cSRPLF7JFXNf48TD+9I+PX6OwpqfNrNm7RZmXzrxF9esaiK59WF06lZe3rk/2fm9GJGyPS4eebxQsSj3ZFLuydztvDvit1ViQjtICtqPXfxf7uYt2rlQsyhN2DWNupetLeZg2y8fkCoiY1Q15kCY+vyVXFnluQ84HXv7gWEfJy0tjQkTJvDCCy9w1113YcXZcwKcfKUxqAxEo46x9O8ylqEy+yjewlO3LaFHxzRG93Db3iVQkbO/u9dLrvrRiCCupoliFIFBSQUMqmM7QiBq7WYciyMeiqMeiiMJFEfcFEc8lEXdnJa5jtGpW5t0+i4UUXpm1f8zsCyLB/58CodcNJsr73+dx66OX66KJ99ahdtlMXVszbuwOqTuEsmtD5bAlKy1zNo8mA+LO3F02pZ4dbXJCLt8FLm6UVTNW6zQFK0wjNXXFisiUUs9WZS60iizWmQGs9E0mx4g2FsKqp36SEQalX/NsPcwdepUFi5cyMcff8zo0aNbuju2F5mQbB/svrl569atLF72HVdddRXWpBttz9Ff4CgN5JFQuIkuZT+QvzWXFI+iWAjqCBInELa8hC1vs2QA8VpRvFaA9k3k1dWXYDhCstdFu6TYgpMG9erMVWeO4Zan3+O84w6Jy97ASCTKM0tWc/yhffboVWYk7RLJre/648jk7RzgLWJefk+OSNnaJvUeVSw7p2lCJturnPdEyuyAm2C+E3iTXxmJmp/QvcX620iaXg+wAhGp6n9bwHDsJNUGAyeddBI+n4/s7OzWYQBr4fnnnycSiXDOOeeAyw3JWfZRhY6DQzz0yld0SwyTrKX4wkX4wgWkBPJICW4jJbgNd8gPWAhRIuImZCUScvlQab0b/xuKPxRlcLf0BgWz/P2csTz/7hqm37OANQ9futuaa0N4d9VaNuUXcc9xJ+6xTIVIbkkgXG/ZIRFbNPfWTUNZUtSVCTHqIrZmQq4kCl1JFPp2GTvRCIn+PCQSqPYTsc3QLHqAFSxnl65jGHvu9cKG3Myw95GamsoJJ5zAvHnzmDVrFi5X6w3VfuqppzjkkEMYMKBm4VSA9EQPo3t34MMfttMlowvFdNm9gKrzq9pOk5Xm30RaYBNpgS1YGgIEUSVsJRBy+QhbiW02ICESjeJ2CR3T6ieXVZ1Er4f7Lz+F4/7yCLc8/R43/t9xjerPk2+tJC3Zy8Qj+tVarku6j683F8Wkuzc4aScDEnfy0o6eHJO2GV8jkwa0ZlRclHqysCS+kcLNSNPrAVagqr0a0rBh32Hq1KnMnz+fDz/8kKOPPrqlu1Mj3333HcuWLeOuu+6qs+xRvTvYicqjUdzV1zVFCLmTKXQnU5jYnc1pzlpUlTRZtmHcSLp/EynBPHsa1S5EWLyEXImELW+rzwJSFoxwYIeUX45BDIwbfhDnjBvKbc9+wJljD2ZAz04Naqcienfq2IPr9CTbJTvTtfXIClNBhRd4fe5w3ijozimZdaU6NrQETpKUS1W1gKbUAxSRkcAGVd3ivD4POwDmZ2Cm2WRuqODEE08kMTGR7OzsVmsAn376aSzL4swzz6yzbLvkBA47IJNl63bSOa2ekXPV0mRtTLejWEUj+MJFjjZdPun+jaQFNpEa3KW/pwghVyIhK7HVKAxU7A/uVsvG9/py98Un8Npn3/H7WQt4f9bvGhQstfDjrykpD3LucUPrLJuU4MLndhGOKm5X/adu+yQWMSx5O6/s7MG49E2kuEysX0MRkUew01XmqeqgGq4fAyzEnlEEmK+qN9bVrqqqiCzAXoprtB5gbX+JD+DstheRMcBt2CK1hTjJpA0GgOTkZE466STmzZtHJNI6FBCqoqo89dRTHHvssXTtWj8Nw6P7dCQSiRKJNi4gQsVFuacd+UkHsiHjUL7sfCof738JSw68lo97/IFVXaayNnMMxd7OJERKSQtsJjWwhVT/ZpKC+XgiZfa+x2amLBihS3piTNOIe6JDRgp3/v54PvziZ/77WvWYuvrx5Fur6NExg6MG96yzrCB0TvfhD8X+tzglay1lUTev7OzRgF4aqvAYMKGOMv9T1aHOUafxq8InjoPWaGqbAnVV8fKmYquyvwi8WCU9mcEA2NOgL7zwAu+//z6/+lXr2iK6dOlScnJymDlzZr3rdEj1Mmz/dqzJLaRTfb3AGKiaP3JbirN+74i27vIWc8nwbyQ1uA0nAwBgEbJ8BF2JRMXTNMnFVYlElR6ZMWw/qYP/mzCcJxav5KoHX+fkI/rRKbP+km5bdxSz+PMfuPrMMfX2HtunelmXH3uWp/29pRyRupU3CrozISOXDHew7kqGX6CqHzhSdk3BWOD3IvIzUMouQdy6FIp+Qa0GUETcqhoGjsXR0KtHPcM+yPHHH09ycjLZ2dmtzgA+9dRTJCYmcuqpp8ZU71f9OrHi551EVbGaI0+WI9oadKdQkNiDTem20oFomKRQAUnBfFKCeWSUbyDdvxFPNB8QZ5uGl6AryVlbbFxf/eEoGUm7a/41FhHhgctP4eDf/ZvL57zKM3+veyq6gmffWUM0qpzrKD/Uh3SfratZq0juHpicuY6lxR1ZsGN/Luj4Q0x1DTExSkRWYwewXKmqX9VVwVkDnI69FNdoavsLfxZ4X0S2A+XA/5wOHETD1CEMezFJSUlMnDiRF198kdmzZ+N2t47fSHl5eTz++ONMnjyZ1NTYhKQ7p/sY3D2D77YU0SG15bJoqLgpTWhPaUJ7trHLW/REykgO5ZMc3E678vWk+zeQGtgCIghRwpJAyJVEyPLFFHATDEcZ2PWXmn+NpW+PDlx79jHMfHwJ540fVm/1+CffWsnwPt3ov3/Het/LZQmdUr3kFQfqnRWmgs4J5YxN28zbhV05sd0GOjRB8vLmIKxCQTiBHWFv5ZEf9rIz7CU/lED/hDwOb3jzbhGpKjX3oKO1Wl9WAPuraomInIAd0dm7rkrOGuAsVR0eW3drZo9/GU6+tSVAF2Bxhe4e9rrhn+Jxc8PexdSpU3nuued49913Oe64xoW8x4vbbrsNv9/P3//+9wbVP7Z/R9bkFsS0sbpZcKJRC9zJFCT2qAy6cUUDJAdto5jh30BG+foqATdKFLcdcONKqnF7RigSxeexaJ/SsK0PdfHXs47m2XdWc/G9C/nyv5eR5Kt9g/3X67ay4odN3HPJnvf+7YleHZLZUuhv0Gd3WtY6PijuzIv5PZne+du6KzQzgajFjgpjVsXAVT0KIwlotR8xCRIh0x0g0+Un2WqURmVYVUc0tHJVoXJVfU1E5ohIe1XdXls9h09EZKSqLmvo/Suo9aeRqn5Sw7nvayprMEyYMIGUlBSys7NbhQHcuHEjc+bM4fzzz6dPn/p5G9Xp3i6Jfl3SyNlWSocahGBbGxHLS5GvK0W+rpVbNKxoyE6RFcwnzb+RzPKfSQ1scbZn2Blvgq4kgq5k/MEI/bqkNtmUrzfBzQN/PoVjLn+YG594h9suqj1O4qm3V+GyLM7cQ+qz2kj1eujWLpFNBX5SfLF5gZnuIOPTN/JawX5MzFwfd3WN+lIU8bA+kMyGQArrg/bj1lAiJdFfRgsnWyEy3QHauQPs7y2pfJ7lDthGzx0g2QojsksNoqVwpPC2Oh7dodiOVX1F15tlDdBgiAmfz8ekSZOYP38+c+bMweNp2ZD+m2++mWg0yj/+8Y9GtXNc/078Z/MPrc8LrCdRy0OJtzMl3s5sTR0I2NszkkI7SXKMYlZ5DqnlG3FHInQhAGV+W47KFf/P8OghB/Cb44fzr+wPOfvYIRx8YJcay0WjUZ5+ezXjRxwUU9BMVQ7okMKmgvIGrQWenLmed4q68tDWvhyRmkc7x4i0cwfIcAXjqroRigobQ8lsCCSzPpDC+mAK6wPJFER2/ehKcwXpkVDC4al5ZLn9lUYty+lTa9q8LyLPAscA7UUkF7ge8ICtBQtMBv7g6LeWA2dWmWWsi+Pj1s/633Pvw+gBxp9FixZx8skn8/rrrzNhQl1R0E1HTk4Offv25be//S1z5sxpVFuqytz3f2JzgZ/M5KYX7G0p8nYWcWz3KOM6+6vpMyq4vLZBdPviEnm6o6iMfhfM4oAumXz079/jqkGI+L1Vaxn754d55m9TOevYITW0Uj9+yCsmZ3tpvXUCq7K4oCtPbOtNpNqOMUHJcAUrDaL9GKx8nukO0M4VINGK7DZcqpAf9lYauPWBFDYEk9kUTCLq3MMjEbollNEjoYQe3lL2Syihh7eEDHejpix3oy3rAcYT4wEa4sr48eNJS0sjOzu7RQ3gTTfdhGVZ/O1vf2t0WyLC+AGdmfPej7RL8rRJL7AuoqoE8TBkcD9I9e7SZyzNs1Uz8r6BvK+haDuVUlPeNEhIaZBBzExL4u4/nMC5t77A/Ys+5ZJTfimn9dTbq0hJTGDSkf0b9d72z0pmfX4ZkUi0RkNbG+MzNjEufRNFkQR2OAElO8Nedka8lc+3hBL5pjyD0hqmJb0SJtMxjBEVNgSTKatSrr27nB7eUkYkb6eHt5Qe3hI6e8pblSxTa0FEiqlBmo9dU6C1qzXXgDGAhrji9Xo57bTTyM7O5sYbb6R79+bPOP/999/z+OOPc9lll9GtW7e6K9SDAzsk0yMziR2lwZiVEdoC+SVBBnZN232d07IgtbN9dB9huy/+Qtsg5v8IW7+Cog1QEWiRkALeVFuUuB5MGzeUJxav5JqHF3Pq6IF0bb/r+6s8EOKF97/g9DGD6gyUqYsEl8VBHVP4bksxqYmxZ6GxBDLcQTLcQQ6gZI/lAlGLneEKw5jAjoi3ymsvgnJEat5unl2Sq/UljmitqGrD5sFrwUyBminQuJOTk8PAgQP59a9/zUsvvdTs9z/77LNZuHAhOTk5dOxY/9D5uvh2SxEPvP8TndJ8eN2tN+l3rKgquQXlXDL2IA7skBJb5VA5FG6EnevsadP8tRB1puost20QPcl79BJ/2pjPoAvv5cTD+zJv5rTK89nvrWHqjc/x9r9+w7HDDmrgO9tFOBrlwx+2Y4ngcbfuHKzNgZkCtTF/CYa406tXL2bOnMmCBQua3QB+8cUXPPfcc1x22WVxNX4AfTulcvrw7hSUhdha5Gdv+fFYUB6iS5qPA9o34PvMkwjtD4Le4+DIy+Cku+HYf8DIC6HnaHvNsCjX9hoLfobS7RD2294kcGC3LK4791e8+MFXLPr4m8pmn3prFV2z0jhmyAFxeY9uy6JPp1TKG5AezbD30qQeoIhMAO4FXMDDqnpbtev9gEeBYcDfVPVfzvn9sPOOdgai2Jss73WuZQLPAz2BdcAUVd3pXLsGW6opgp0x/M3a+mc8wKYjFAoxYsQI8vPz+frrr0lLi3l6vkGcdtppLFmyhJycHDIzM+uu0AB2lgZZuHojqzcUkpmUEHOIfWshHI2ytchPqtfDeaP2p1es3l99CfmheJPtKW77DvJ/gECxc1EISiLDZjxJUVmArx+dQXkgRNczbuPyyUdyx+/jFvBHVJWlP+UTikTxxiHHaVumrXqATiaY7qq6IR7tNZkHKCIu4D7skNUBwFkiUl2IbQdwKfCvaufDwBWq2h9b7v6SKnX/CixR1d7AEuc1zvUzgYHYSVjnOH0wtAAej4eHHnqITZs2xSUQpT4sX76cl156iSuuuKLJjB/YahHnj+rJ747qhaJs3FlOONJ6QtDrw47SIFsK/Yzp3YGrJvRrOuMH4PFB5gHQ6yg49Lcw4Tb49c0w6o/Q/yQSMrvy4PSj2JBXyD/uf4nn3/iIcCTKucf0g0iw0ltsLJYIfTunEghH49amoXlxtkosiFd7TeYBisgobNmkXzuvrwFQ1VtrKDsTKKnwAGu4vhCYrapvich3wDGqullEugDvqWrf6u2LyJvO/ZfuqY/GA2x6Lr30UmbPns3SpUs57LDDmvReJ5xwAp9++ik5OTnN5nH6QxHe+TaPd77Nw+u2yEpOaNVRooFwhK1FAbq3S2TKiP3YL44JrxtFNMIffvcbHnzsKfbrlEl6UgKrH5gOgSIIVvyPOtGnqG3AVEFc9l5Fl8dec3R5wPKAVfNvX0VZlrODkkCYxIS26bnHg7bqATr3vg94rMkzwTSSbkBVNzUXiPkb0MkofgjwqXOqk6puBnCMYMVCTzegauaaXOdc9fYuwknsnZCw90XztTb++c9/Mn/+fC666CI+//zzJtsc/9FHH/H6669z++23N5vxA/B5XJwwuAtD98tg3vJc1m23M8bEQ0YonqgqecUBFDh9WDcOPyALd4xbApoUy8Wtd93LgtcW8/PmLdx5551w7JX2tWjUXjcMlUOo1H4MOo/+QijfaT8GCsFfBGX5EI04+U8FexUFcHkRTyJ9Ovr4dF0xvjaa2MDQNjLB1PSXFZO7KSIpwIvAjKq54xpzPydh64Nge4Cx9McQO2lpacyePZtTTz2Vu+++m6uvvrpJ7nPdddfRqVMnLrnkkiZpvy66ZiRyydiD+PznHSxcuYmCsiAdU30xZx9pCkoDYfJLgwzqlsaph3RvtZv5MzIyuP/++7n00kuZNm1XRCiWBQlJ9kFW3Q2p2lOnwVLbgyzfCaX5ULQRireQUbKFAzw7KCiN4PO4sIgQkQTClpew5SXSVDJThngRt4XhpjSAucB+VV53x5a9qBci4sE2fk+r6vwql7aKSJcqU6AVmX4bdT9D03HKKadwyimncMMNN3DGGWdwwAHxieyr4J133uHdd9/l3nvvJTm55SKzXZZwWK8s+nVO49U1m1i2bicZiR7SElsmJVwkqmwp8pOc4OLC0T0Z2DW91Xs8kyZNYtKkSY1rRATcXvtIyoR2PXe/rkrajp089+Zy9vMFSYoUkRLcRkpwq51IPJKPIiiCRZQoLkIuHyErkajVsun9DKCqP4tIO2z1iKoyLTFLJDXlGqAb+B5bS3AjsAw4uybNp+prgE6kz+PADlWdUa3snUC+qt4mIn8FMlX1KhEZCDwDHAp0xQ6Q6a2qe4x7NmuAzUdubi4DBgxg1KhRvPHGG3H7IlZVjjzySDZs2MAPP/yAz9dyskVVUVV+2lZC9ue55JcE6JTmw9OMU447y4KU+MMc1bs9vx7UmaR9eL1rT8xbvoFlOTvonJ6423krGsYbLsLnHMnB7aT7c0kLbMYd9QMWaJSIlWAra8QoN9UaaONrgL8FLsN2clZhB0ouVdWYhUib7L9CVcMi8kfgTextEI+o6lciMt25fr+TEfxzIA2IisgM7IjRg4FzgS+qqM9fq6qvAbcB2SJyIbAeOMNp7ysRyQa+xo4ivaQ242doXrp3784tt9zCn/70J5599lnOPvvsuLT7+uuvs3TpUh544IFWY/zATp92UMdUrhjfhw++38bir7fisYT2Kd4m9cKC4Shbi/10Sfdx4ehe7J/V5vcqNxnj+ndiWc5OQpHobj9Oopab8oRMyhOqRRKrkhApJSm0g6RgPun+XDL8G0kNbqPqakvI8hFyJRKRBDOV2jRcBowEPlHVsc52uhsa0pDJBGM8wGYjEolw5JFHsnbtWr799ttGb1VQVYYPH05hYSHffvtti6tP1EZekZ/5K3P5bksJ7VMS4u6RqSrbSoJEolFOGNyFIw9q36weZ1tl8VdbWPz1VrplJNZdeA+IhkkMFZIU2kFKYBvp/vWk+zfhjRRT4S1GpUKHMRFtBbuz2rgHuExVRzrO0WGqGhCRVao6NNa2zLyIodlwuVw8+OCDDBs2jKuuuoqHH27YP18FL730EitXruTxxx9v1cYPoGOaj4uOOpDVGwqYv3Ij+aVlWAgul+CxBLfLwm0JbpfgtixcMQTPlAXDbC8JMKBLOqcO69ZkYrZ7I6N7t+f977fhD0UaHLmr4qYsIYuyhCy2J/cGjgDAHfHb3mJoB6n+TWQ406hWxcSUaqVRNN5iTOSKSAb2fsC3RGQnDYz3MB6g8QCbnauvvpo77riD9957j6OPPrpBbUQiEYYMGUIkEuHLL7/E5Wr5X9X1pSQQ5rstxZQGwhT5QxT7wxT7Q5QEwpT4w5QGI0Qian8fSkWMtz3JpqqVBtLjEspDEbxuF6cP68aQ/TJafZBLa+SjH7cxf8VGurdrhj2RGq1cV0wO5tGufAPp/o34wsVExUKIEhEPIVdSk64ttmUPsFo/jgbSgTdUNRhzfWMAjQFsbsrKyhg0aBAJCQmsXr0arzd2j+WZZ55h2rRpPP/880yZMqUJetlyqCqhiFIeiuB3Dvt5lPJgmCJ/mKLyEEX+EFnJXsYN6ESK10zmNJRgOMptr3+DJUJyC42jO1JuR6AGt5Me2EhG+QZSgttQQIiiWE7ATXwiUfcWA9hYjAE0BrBFePPNN5kwYQLXX389M2fOjKluOBymf//+JCUlsXLlSizLrHUZGsfK9Tt56pOfm8cLrCeiYZKCO+0p1MBmMso3kB7YhDvir8yIExYvIVcSYcsb0xRqWzSANegBOmmBGq4HaAygMYAtxrRp05g3bx6rV6+mX79+9a73yCOPcOGFF7JgwYLG7xkzNBuhUIjc3Fz8fn9Ld+UXqEKRP4SqYrXyaWRBEY0gGsXSMJZGEKrnorX3Me7RKKoCije19sQCPp+P7t27/2KN3XiAewHGALYseXl59OvXj8GDB/Puu+/Wy5MLBAL06dOHjh078tlnn5k1rzZETk4OqampZGVltcrPzR+KsK04gMdltb14FI1iaQRLw7iiIdwaxIqGkCoOk2KhYtmGEQWNktKx556bVCU/P5/i4mJ69eq127WWNoAiMgQ4ynn5gaquaUg7Zu7I0GJ07NiRO++8kw8++IBHH320XnX++9//sn79ev75z3+2yi9Rw57x+/2t1vgBeN0WXo9FJNq2lD0AEIuo5SHsSiTgSaM0oT3F3s4UeztS6snC704jbCUAiktDWBquu0kRsrKyWp3HLiKXAU8DHZ3jaRH5U4PaMh6g8QBbElXlmGOOYc2aNXz77bd06tRpj2XLy8s58MADOfDAA/nggw9a7RepoWa++eYb+vfv39LdqJWgo5bRJr3AeiIaRaJhVJWOmel1lq/pc2vhfYBrgFGqWuq8TsbOBBNzMmzjARpaFBHhgQceoKysjMsvv7zWsnPnzmXz5s3G+zM0GQluF0kJLsJt0QusJyoWYcvjeIRtEsEWPa8gQs1iCHViDKChxenXrx/XXHMNzz77LG+88UaNZUpKSrj11lsZN25cg/cOGgz1IS3RwxkTJ7Byxect3ZUa+eh/H7Ds0z3KnO4LPAp8KiIzReQGbBm8/zakIWMADa2Ca665hr59+3LxxRdTVlb2i+v//ve/2b59OzfddFML9M6wL+Fx2YkGIjFkEo7EUriRfPzhByz79JMar4XDda/ttXVU9W7g/4AdQD7wf6p6T0PaMrtnDa0Cr9fLgw8+yNFHH80NN9zA7bffXnmtoKCAO++8k5NOOonDDz+8BXtpiBczZsxg1apVcW1z6NCh3HPPPTVeKy0tZcqUKeTm5hKJRLjuuuuYOnUqS5Ys4corryQcDjNy5Ejmzp2L1+u1U9EJRFX5658vY9WK5fj95Zw06VSuuvY6AEYM7stZ55zP+++8zW9+N51TJu9KyPDySy9y1+234HK5SEtLY8Hrb/Pc00/y+isvEwwEWP/zOk49YypX/vVvAMx7/lkevv8+QqEQw4aP5La778XlcvHO24u59YbriUQjZGZmcffs+3nikYdxuSxezH6Om++4i2eefJx27drxxZrVDB4ylJSUVJKTk7n4UntJ4ejDh/Pk87ai3Nmnn8yhhx/B8s8/Y8DAwVz8+99y/fXXk5eXx9NPP82hhx4a188knojISGCDqm5R1RUiMhg4HThIRH5W1R2xtmkMoKHVMGbMGC688ELuuusuzjrrLIYOHQrA3XffTUFBATfeeGPLdtDQZnnjjTfo2rUrr776KgCFhYX4/X4uuOAClixZQp8+fTjvvPOYO3cuM2bMQIDkBDfhiHLNdTNpl5lJJBJh8snH8/WXXzBg0GDA/uH28pvv/OJ+d99xK8/Nf5kuXbtRWFBQeX7l8s9575PPSUxMYsLY0YwbP4Gk5GQWzp/HosXv4vF4uPrPl/Fi9nP86rjxXHnpxbz02tvs37MnO3fsoF1mJuf95re7Gbhnnnycn378kRcWvobL5eLOW/6JYutBqqqTWche08xZ+xNzHnmSW+7+DycfN4ZnnnmGDz/8kJdffplbbrmFBQsWNOXH0FgeAMYBiMgY4FbgT8BQbJHzybE2aAygoVVxxx13sGjRIi666CKWLl3Kzp07mTVrFpMnT+aQQw5p6e4Z4sSePLWmYvDgwVx55ZVcffXVnHTSSRx11FGsXr2aXr160adPHwDOP/987rvvPmbMmAFAUoILS4SFL83j6cceJRwJs3XLFr7/9ptKAzjptJq/cw89bBSX/eEiTj71dE44eVeyhjFjf0Vmpr35/ISJk/jsk49xud2sWbWCCWNHA+AvL6d9hw6sWPYZhx8xmv179gSgXWamnRNWlWgUQpEoqhCNKidMPIWICpFwlKiTONZtCS5nOjcj0YMlQs9evRg7agSWJQwZPIhjjz0WEWHw4MGsW7euCUY+rriqeHlTgQdV9UXgxSqyeTFhDKChVZGZmcmsWbOYNm0ac+bMYcOGDZSWlnLDDQ2S+zIYAOjTpw/Lly/ntdde45prrmH8+PGcfPLJtdaxLGHn1lzm/ude3nz3QzLatePSP/wOfyBQWSYpueadAHfc8x9WfP4Zb735BuNGH8bbH34K8IvoZRFBVZly1jn8beau9W1VeOPVRSh2rtKqOb8EWzUkxevGbQlej0XnrAw6p/uwRGiX4sPrddExzdbHDAYCJDk5Tn1eL25HJsuyrMo8vJZltYX1Q5eIuFU1jC20flGVaw2yZSYIxtDqOOussxg/fjzXXnsts2fPZtq0aQwYMKClu2Vow2zatImkpCTOOeccrrzySlasWEG/fv1Yt24dP/74IwBPPvnkLyKMQ/5SkpOSSUxOZdOmzbzz1ptEIlF7StHxvqKqVN9OvW7tWoaNOJSr//YPMjPbsyk3F4AP3l3Czh07KC8v541XFzHisFGMHnMMixa+xKbNWwiGo+Rt207OunWMOOwwPv34Qwrycmmf6sUbKaNbu0S6dGhHNFhORlICKT4PbsvC47Yqg3d69erFihUrAFixYgU5OTlNP8DNw7PA+yKyECgH/gcgIgcBhQ1p0HiAhlaHiDB37lwGDRpEMBjk+uuvb+kuGdo4X3zxBX/5y1+wLAuPx8PcuXPx+Xw8+uijnHHGGZVBMNOnT9+t3iFDhzJixDDGHzWSXr16ccQRR5LkdZPqdSNie4mq7LZvUICZf7+GnLU/oqqMHjOWfgMH88Wa1Yw8/Aguueg3rFu7lkmTpzDw4ENwWcK1183k/DMmoRrFk5DA7NmzGTSkHw8/9BDnnTWVaDRKx44deeutt5g4cSKTJ09m4cKF/Oc///nFez399NN54oknGDp0KCNHjqyc4m3rqOrNIrIE6AIs1l1ZXCzstcCYMZlgTCaYVsuiRYvIz8/nggsuaOmuGOJAW8gE01BUlajaUaMRxyuMRpWwc0SiyjNPPcGXq1Zw9z3/xuN2BJAtCysG8eOWoLVlgoknxgM0tFomTpzY0l0wGOqFiOAScCHsSVi+XVICiQlu0pPabAaWvQ5jAA0Gg6EZuOCCC8xsRiujSYNgRGSCiHwnIj+KyF9ruN5PRJaKSEBErqx27RERyRORL6udf15EVjnHuorwVxHpKSLlVa7d35TvzWAwxM6+vOTSFtnbP68mM4Ai4gLuA44HBgBniUj1UL4dwKXAv2po4jFgQvWTqjpVVYeq6lDgRWB+lcs/VVxT1enV6xoMhpbD5/ORn5+/13+p7i1U6AH6fL6Y6+7JgalyXUTk345ztEZEhjW6ww2gKadADwV+VNW1ACLyHDAJ+LqigKrmAXkicmL1yqr6gYj03FPjYm+omQL8Ks79NhgMTUD37t3Jzc1l27ZtLd0VQz2pUIRvAI8Bs4En9nD9eKC3cxwGzHUem5WmNIDdgA1VXucS3zd4FLBVVX+ocq6XiKwEioC/q+r/4ng/g8HQCDwezy+UxQ17J3U5MNjO0BPOVoZPRCRDRLqo6ubm6aFNUxrAmmJ74zn3cRb2xsgKNgM9VDVfRIYDC0RkoKoW7dYpkYtwMggkJJhoLIPBYGgAbhGpqhf1oKo+GEP9mhykbtjf481GUxrAXGC/Kq+7A5vi0bCIuIHTgOEV51Q1AASc58tF5CegD7CbqJfzIT0I9j7AePTHYDAY9jHCqjqiEfWb2kGqF00ZBboM6C0ivUQkATgTeDlObY8DvlXV3IoTItLBCbxBRA7AnlteG6f7GQwGgyF+NJmDFAtN5gGqalhE/gi8CbiAR1T1KxGZ7ly/X0Q6Y3toaUBURGYAA1S1SESeBY4B2otILnC9qlao/p7J7tOfAGOAG0UkDESA6XXpQ5WVlamIlMfj/bYC3ECrz2bbgpjxqR0zPnvGjM0vSWxk/ZeBPzrBkYcBhc29/gf7eCq0vQkR+byRUxJ7NWZ8aseMz54xYxM7VR0YYCtwPeCBSudHsKNEJwBl2Krun9fcWtNhMsEYDAaDIa6o6ll1XFfgkmbqzh4xckgGg8Fg2CcxBnDvIZYQ5H0RMz61Y8Znz5ix2Usxa4AGg8Fg2CcxHqDBYDAY9kmMATQYDAbDPokxgG0EEfGJyGcislpEvhKRG5zzQxxJqS9EZJGIpDnnjxOR5c755SKy1yYNj3VsqtTrISIl1aW49jYaMj4icrBz7SvneuySAG2ABvxfeUTkcef8NyJyTcu+A0OjUFVztIEDO3VQivPcA3wKHI6dcedo5/xvgJuc54cAXZ3ng4CNLf0eWsvYVKn3IvACcGVLv4fWND7Y26PWAEOc11mAq6XfRysZm7OB55znScA6oGdLvw9zNOwwHmAbQW1KnJce51CgL/CBc/4t4HSn/EpVrUgt9BXgExFvM3a52Yh1bABE5BTsVHlfNV9PW4YGjM94YI2qrnbq56tqpBm73Gw0YGwUSHbyEScCQWz1GUMbxBjANoSIuERkFZAHvKWqnwJfAic7Rc5g9/x6FZwOrFQ7YfheSSxjIyLJwNXADS3Q1RYhxr+dPoCKyJsiskJErmr2DjcjMY7NPKAUW7VgPfAvrSPloqH1YgxgG0JVI6o6FDtx7KEiMgh7euYSEVkOpGL/Iq1ERAYCtwO/b+buNisxjs0NwKwqv/z3emIcHzcwGpjmPJ4qIsc2f6+bhxjH5lDsXMNdgV7AFU7yfUMbxKRCa4OoaoGIvAdMUNV/YU9ZISJ9gBMryolId+Al4DxV/akl+trc1HNsDgMmi8gdQAZ2Ina/qs5ugS43K/Ucn1zgfVXd7lx7DRgGLGn+Hjcf9Rybs4E3VDUE5InIR8AIjPJMm8R4gG0ER+4pw3meiCMJJSIdnXMW8Hfgfud1BvAqcI2qftQSfW4uYh0bVT1KVXuqak/gHuCWvdn4xTo+2AouB4tIkrPWdTTwdbN3vBlowNisB34lNsnYATPfNnvHDXHBGMC2QxfgXRFZgx2h9paqvgKcJSLfY/8TbgIedcr/ETgIuE5EVjlHx5boeDMQ69jsa8Q0Pqq6E7jbKbsKWKGqr7ZEx5uBWP927gNSsNcIlwGPquqa5u+2IR6YVGgGg8Fg2CcxHqDBYDAY9kmMATQYDAbDPokxgAaDwWDYJzEG0GAwGAz7JMYAGgwGg2GfxBhAwz6FiMwSkRlVXr8pIg9XeX2XiPw5jvd7TEQmx6u9Ku1eW+V5TxH5sp59yRGR6XHqw+Uisl5E9to9lIa9G2MADfsaHwNHQOUm5/bAwCrXjwDaQuKAa+suUiN/UdX76y5WN6o6C/hHPNoyGFoCYwAN+xof4RhAbMP3JVAsIu0ctYz+wEoR+YeILBORL0XkQSfzR38R+ayiIcfzWuM8Hy4i74utvfimiHSpfuM9lRGR90TkdrF16b4XkaOc80kiki0ia0TkeRH5VERGiMhtQKKT3OBpp3mXiDwktqbdYierSa04HuG/ReRjEVlb4amKyDFOP7Od/twmItOc/n0hIgc2dPANhtaEMYCGfQpHIiosIj2wDeFSbA24Udg5HdeoahCYraojVXUQtuzNSar6DZAgu5IfTwWyRcQD/AeYrKrDgUeAm6vetx5l3Kp6KDADuN45dzGwU1UPBm4Chjvv4a9AuaoOVdVpTtnewH2qOhAooIr0Ux10wU54fRJwW5XzQ4DLgMHAuUAfp38PA3+qZ9sGQ6vGJMM27ItUeIFHYKf86uY8L8SeIgUYK7YMUBKQia0buAjIBqZgG4upztEXW3T4LREBcGHL5VSlrjLzncflQE/n+WjgXgBV/bLC29wDOaq6qoY26mKBqkaBr0WkU5Xzy1R1M4CI/AQsds5/AYytZ9sGQ6vGGEDDvkjFOuBg7CnQDcAV2MKmj4iID5gDjFDVDSIyE/A5dZ8HXhCR+dh6qj+IyGDgK1UdVcs9pY4yFVqNEXb9X0oM76mq1mME22uNtZ7s4Xy0yuso5nvDsJdgpkAN+yIfYU/57XC04HZgyyKNwp4SrTB220UkBaiM4nRkpSLAddjGEOA7oIOIjAJ7ulNsHcaq1KdMdT7E9jYRkQHYBruCkDOtajAYGogxgIZ9kS+woz8/qXauUFW3q2oB8JBzbgF21v+qPA+cgz0dirNmOBm4XURWYysoHFG1Qn3K1MAcbKO5BlvBfg32NC3Ag8CaKkEwBoMhRowahMHQShERF+BRVb8TebkEOxglWEfVPbX3GPCKqs6LYx8vwJ4q/mO82jQYmgszl28wtF6SsLXqPNjrc39oqPFzKARuEpH28dgLKCKXA9OBFxvblsHQEhgP0GAwGAz7JGYN0GAwGAz7JMYAGgwGg2GfxBhAg8FgMOyTGANoMBgMhn0SYwANBoPBsE/y/8obYIw8z618AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# interpolate the solar spectrum to the radiance grid for normalization\n", "wl = engine.wavelengths\n", "solar = np.interp(wl, engine.optical_property_wavelengths, engine.solar_irradiance)\n", "\n", "# plot inelastic radiance\n", "rad = engine_output.radiance[:, 0] / solar\n", "stdev = np.sqrt(engine_output.radiance_variance[:, 0]) / solar\n", "plt.plot(wl, rad, 'C0', label='inelastic radiance')\n", "plt.fill_between(wl, rad - stdev, rad + stdev, color='C0', alpha=0.5)\n", "\n", "# plot elastic radiance\n", "rad = engine_output.elastic_radiance[:, 0] / solar\n", "stdev = np.sqrt(engine_output.elastic_radiance_variance[:, 0]) / solar\n", "plt.plot(wl, rad, 'C1', label='elastic radiance')\n", "plt.fill_between(wl, rad - stdev, rad + stdev, color='C1', alpha=0.5)\n", "\n", "plt.xlabel('Wavelength [nm]')\n", "plt.ylabel('Sun Normalized Radiance')\n", "plt.legend(loc=2)\n", "\n", "# plot solar spectrum for reference\n", "plt.gca().twinx()\n", "plt.plot(wl, solar, 'k', label='solar spectrum')\n", "\n", "plt.ylabel('Solar Irradiance [photons/s/cm$^2$/sr/nm]')\n", "plt.legend(loc=4)\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 }