Histogram.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "8484e32d",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.0077531107247804515\n",
      "0.09332304232182348\n",
      "3.140645713444366\n",
      "4.567539198743687\n",
      "[[0. 0. 0. ... 0. 0. 0.]\n",
      " [0. 0. 1. ... 0. 0. 0.]\n",
      " [0. 0. 0. ... 0. 0. 0.]\n",
      " ...\n",
      " [0. 0. 0. ... 0. 0. 0.]\n",
      " [0. 0. 0. ... 3. 1. 0.]\n",
      " [0. 0. 0. ... 0. 0. 0.]]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyEAAAKrCAYAAADmsoL5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACQcElEQVR4nO39eZRc12Hf+/52DV09FXpCA2iMBEESpEiKpEhRo21aHq5sK1IG3US5uY6c6zzlZcUrzkryEjt/OOt63T+Se99zkhvnxVfL9rKc+HmIPEj2lWzLsiibtiaS4gSSIDEQUzeARk+onqu69vsDrYSh928DXWpUF9Dfj1YvgXv3PmefffY5VafPqV+FGKMAAAAAoF0KW90BAAAAANsLFyEAAAAA2oqLEAAAAABtxUUIAAAAgLbiIgQAAABAW5W2asVdoRK71bdVqwcAANi2QghtWU8uhTUU/N/CY7O5qf2oaeZKjHF0Uxd6E/wP390Xp6bXtmz9z7yw8gcxxg+2Y11bdhHSrT69K3zPVq0eAABg2yp0dbVlPc3VVVtX7K/aurVabVP78Ufx02c2dYE3ydT0mr7+Bwe3bP3Fsdd3tmtdPI4FAAAAoK227E4IAAAAgP8mSmpqcx9F61TcCQEAAADQVlyEAAAAAGgrHscCAAAAOkLUWtwej2Pd8EVICKEo6WlJF2KMH3pL3Y9I+j8kXVgv+tkY489vVicBAADwF+VSrnLJVLm6zVSs+gSsVrS6veg8G7kT8uOSXpG0w9T/eozxx779LgEAAADbz7UPpvvvVrmd3NBnQkII+yX9kCTubgAAAAD4ttzoB9P/raR/JmUzw/5aCOGFEMKnQwgHUr8QQvhECOHpEMLTda1ssKsAAAAAbgfXvQgJIXxI0uUY4zOZX/tdSXfEGN8u6QuSPpX6pRjjJ2OMj8UYHyur0lKHAQAAgNtVcwv/1043cifkfZI+HEJ4Q9KvSfpACOE/v/kXYoxTMcZv3dr4eUmPbmovAQAAANw2rvvB9BjjT0r6SUkKITwh6Z/GGP/nN/9OCGEsxjix/p8f1rUPsAMAAAC4QVFRa3F7fDC95e8JCSH8tKSnY4yflfQPQwgfltSQNC3pRzanewAAALePXMSs02r0bGlkxNY1pqY2vLxW+r5Wq224TbYPmcjfZgvbhK2zoYuQGOOTkp5c//dPvan8v94tAQAAAIAcvjEdAAAA6BB8TwgAAAAA3ATcCQEAAAA6QJS0xp0QAAAAANh83AkBAABogUuLajXNyilmEqFyYqYfbpmhhQQsyadt5RK1Whmn5ianbWHrcBECAAAAdAg+mA4AAAAANwF3QgAAAIAOEKVt843p3AkBAAAA0FZchAAAAABoKx7HAgAAADpEc6s70CZchAAAgG2hNDKSLHfxslI+YrZgYm6bmeWFew7buvja6XSb4UHbpnHmnK3LRfuumahbN0ZSfpxaihHe5Ajh3Lij8/A4FgAAAIC24k4IAAAA0AGiotb4nhAAAAAA2HzcCQEAAAA6QZTWtseNEO6EAAAAAGgv7oQAAIBvWytJTK0uL5eQlEtwiiaNKZeAFSqVDa8rmxR1ZtzXGXF61tbl+p5jU7AqfnnZfVLtT1es+ASsHLe8WJtvaXnoPFyEAAAAAB0gavt8TwiPYwEAAABoK+6EAAAAAB0haE1hqzvRFtwJAQAAANBWXIQAAAAAaCsexwIAAAA6QJTU3CbfE8JFCAAAtzEX4do0cbWtiisrts7GwUpqthDfm2tT2jvmG5q42ObbDtsmhRMX/LruPpIsjxcv2zZheNDWxYlL6TYu/lb5cc9FGbu44rVMxHFubBvjE8nyXKxvbl64KN5W4p7RmbgIAQAAADoEH0wHAAAAgJuAixAAAAAAbcXjWAAAAEAHiOJxLAAAAAC4KbgTAgDANuRSsySpkEk0clpJuZKkUKlsuE1hbLfvh0mYkiTdf1eyuJhJwHIpUpIUFhbTfcilbS34NKtC3adF2TY7N95Gkv1be25euHQxyadg5RK6cnPGpbfl+rfZiW9bpRm5EwIAAAAAm46LEAAAAABtxeNYAAAAQAfgg+kAAAAAcJNwJwQAAADoAFFBa9vkHsH22EoAAAAAHYM7IQAAtFkrMaMuAlWS1jJRpy5uNxc9m+NiVbOxvhW/ve7p91x8a6g3bF2uH7G+lq7o7/XrmrdVUjn9Nqp4btL3oZZZ4PBgus30rG0S9uzyy8u0c+vK/XW6lRjm3NxsRTbS+TaJ6N0uuAgBAAAAOgTfEwIAAAAANwEXIQAAAADaisexAAAAgA7A94QAAAAAwE3CnRAAAG6CXJpVXFmxdaVDB5LlzYlLLa3LJT+Feiah68qUrzMJRIVMwlRz2PevMG3StjKJWi6VSpLU5/sRXzudLA9ju/3yMv2IVTO2Kz69K+T6d/FyuuLQXttGE1d8XabvzTPnkuXZ9KmMghvDzLxtZV25Y+f2ELQWt8c9gu2xlQAAAAA6BhchAAAAANqKx7EAAACADhAlNbfJPYLtsZUAAAAAOgZ3QgAAAIAOQUQvAAAAANwE3AkBAEDXibk1cnGhYc8uX5dbpolpddG4klTKxQG7GNlMtGu457CtWxvtT/fh+Lhts+nqPgI3Jz5yb7rixAW/qvvSkcmSVD5nooxz/cvEC4dqemxVW7Rtmgcz8cIZBdePTN9t/+QjpAu5uZmZ07j9cRECAAAAdIAYO/97QkIIb0iqSVqT1IgxPhZCGJb065LukPSGpL8eY5zJLaeztxIAAABAp/nuGOPDMcbH1v/7JyR9McZ4t6Qvrv93FndCAAAAgA7RvDU/mP4RSU+s//tTkp6U9M9zDbgTAgAAAECSdoYQnn7TzycSvxMl/WEI4Zk31e+OMU6s//uipOt+WIk7IQAAAAAk6cqbHrFy3h9jvBBC2CXpCyGEV99cGWOMIYR4vRVxEQIA6FiFri5bl0uL2vSkq0olWZ5L/tH0rF9Xpu/BbHPpkE9pajUtyjrjk64qKyb1K5P6VJj3Y+vE4R22Lqz47V0b7vP9WNh4P8qXan5de4aS5cVzk7ZNc9jPmVbGqTDt+6d5n6plVfwxl+OOx5CZ67ljLnd8386ipLUOf1Apxnhh/f8vhxB+W9Ljki6FEMZijBMhhDFJ6Zi/N+nsrQQAAADQEUIIfSGE6rf+Len7Jb0k6bOSPr7+ax+X9JnrLYs7IQAAAEBH6PiI3t2SfjuEIF27jvj/xRh/P4TwDUm/EUL4UUlnJP316y2IixAAAAAA1xVjPCXpoUT5lKTv2ciyOvpSCwAAAMDthzshAAAAQAeIkprb5B7B9thKAAAAAB2DOyEAgA3JxebmImubNR8l6trlomyLJjZXkkK1P12RiZFVtdcvb/pquqLPt8lFxYbVTKTuiXPpNvVMdGqmH43BnmR5oW+fbdMs+79RlmaX0hUVP7bLezPzoiu9rp4L87ZNLBdtXU6opSNrlx49bNt0j2fmrYv8zcyzUF+zdc3+9JyOXX55xWlbpebBTBzw2UvpddX8uLv4aMnHYttj8TrLy0X7OmuZc8ytZC3ekt+YvmE3fCckhFAMIXwzhPB7ibpKCOHXQwgnQghfCyHcsam9BAAAAHDb2MjjWD8u6RVT96OSZmKMd0n6N5L+9bfbMQAAAAC3pxu6CAkh7Jf0Q5J+3vzKRyR9av3fn5b0PWE9QBgAAADA9UUFramwZT/tdKNr+7eS/pmkpqnfJ+mcJMUYG5LmJI289ZdCCJ8IITwdQni6LvMcJQAAAIDb2nUvQkIIH5J0Ocb4zLe7shjjJ2OMj8UYHyvLf6AQAAAAwO3rRtKx3ifpwyGEH5TULWlHCOE/xxj/5zf9zgVJBySdDyGUJA1Imtr03gIAtlwuASunmUm7cUk4YXjQtonTs35lLi0q1weTnCTJph3Fi5d9m2GfuLQy6tOseqYH0+vKpHe5BCxJWh1MJxD1vjFn2xRXfHqXS3DKbVMuYcqN+9qeIdsmp9FbtnUuU6vnpQu2TRzeYeuafemxKLo0NeXnWePAX3iIRJJUqLsHUfIpbMWLM7ZOlfS8yD1Lv3aXT1QrnptMljev+LeDhUMH/MrMsXW7JGDlNOP2+AaN625ljPEnY4z7Y4x3SPqYpD9+ywWIJH1W0sfX//3R9d+Jm9pTAAAAALeFlr8nJITw05KejjF+VtIvSPpPIYQTkqZ17WIFAAAAwA2KUts/IL5VNnQREmN8UtKT6//+qTeVL0v6HzezYwAAAABuT9vjUgsAAABAx2j5cSwAAAAAmycqaC1uj6/a404IAAAAgLbiTggA3OYKXekozutxkboxF3ObWVcpF8dZ95GwrazLRfHmYl+XR7ttnYuzXXxwzLbpO+6jSVcPZuJsTRTv/N2Dtk3v+JJf3uRystzFy0pSY7Tf1lXOz6bLbYt8LG3dxNIWFzLzrL5m65oDfl64mOPlvT52umvO98NZObrX1jXL/q/cbl25udk16/tXyMQ6h0wMs1OcXvCVJsa6MLbbt8mcS5zc+SxUMrPQpyZ3nOY2uUewPbYSAAAAQMfgIgQAAABAW/E4FgAAANABYpTW+MZ0AAAAANh83AkBAAAAOkJQU9sjopeLEADYIi7lJZfwEldWbF02GaaF5bk0qzg9u+H1XFefT/FxVt6+y9YVTHpSod60bUoLPnEpTKejdXK9Xjrsk7gGvn7B1l15Yn+yvG+inlmb1yynH3oozfpErdybg/rudJJUs1y0bWpvH7B1BRPSVD3j+7A66BOScmlRLmGse9IfB278JGl5NH3MVV+6YtvkEqFceptLOLueXAJWrKT3cpjOJFatZOoq6X3SvOJT4go708lorVqr1TZ1ebi5eBwLAAAAQFtxJwQAAADoAFF8MB0AAAAAbgruhAAAAAAdYm2b3CPYHlsJAAAAoGNwEQIAAACgrXgcCwBukIvUbWYiN1sRzHqu2254MF1R9zGdYY+PuW32p+NHs3+9yvTdRbtKPjq3eHHGtqmcn7V1ru/1wR7fpuyz+ePwDlvnNPp8ZG1unHb++eVkeW78Gr1lW1c7mF5XebHbtuma83OmtJiuc+uRpOEX5myd2ydzR/y+Gjjp44VzcjHMTi7Wuf/12WS5m3/X+LriCRPd3J8Jg87MpeaZc7Yu3HM4XeHOI5Iar5+0daW7jyTLC5lY3+xxVZv3dbexqKBm3B7fE8KdEAAAAABtxUUIAAAAgLbicSwAAACgQ5COBQAAAAA3AXdCAAAAgA4QJTW3yTemcxEC4JbmEqtCxSfQxJUVW5dLuipUTTpRrWbb5PrhUrDi2E7bZnW039Z1Pf1auuLQXttmea9PXOoeN9uVSeOZe3jU1u14xSddNUxC0urR3bZNoR5tXWmxniwvz/pUpZVRn0AUy+mkq9m3Ddg2gy/7RKjcOHVf2XjfV/f5eVE9m57TufQul4AlSauDJiUuky525dFBW+f61zPlk6zK56Zs3dzj+2xdvTfdx6FX/XGfSx5b3Zuet92T/hxTO5RJ/epKvy0Lq35/5OZtYf+grXP7q+clk9AlqejOgZI0PZtez0F/DOvYCVvlzrdrmfMtbi3b41ILAAAAQMfgTggAAADQEYLWxPeEAAAAAMCm404IAAAA0AG20wfTt8dWAgAAAOgYXIQAAAAAaCsexwLQVi5SNxuNm4mELYyZ+Me6j7R00biSVKj4OpXTp8zC8KBvY2IrJR/FG6av2jaV2qKtW3vb4WT56oDfpmwUq4kDzkXFDnzdx3uOf/iArdv99Xlb59hIYklXP3h/srzQ8LG+SyM+snZpNB3TWj3jx6JuYoclafIR/zfAA19Ilxfmfezrwp5BW9d3MV3e6PF9KC34sVgZSNeVlvzY9kymY4clHwXd7PPx1rV3+Njprjk/pwuNdN/nD/qY2+oJfzwW6pnzhTH0/LStc1HQhbOXbJvuuo/0bimCu8+PRVjx52mZ82DuOPDhx1I0rwmlkRHbpjHlo5tvJXwwHQAAAABuAu6EAAAAAB0gxsAH0wEAAADgZuAiBAAAAEBb8TgWAAAA0CHWtsnjWFyEAMja7DQrp9XEk4JJwWruHvJtpk0qjCTl+u62OTMWa3fts3XFc5PJ8vqdJvFLUqPXn7Z7jqcTdErlzNhmlrcwZhKhbAupkklwGjidSeIyCV5dc35sl99z1NatdaXTZXLpWNWzfl0uRWxxrNu2yaU07XvSr6tZTvf9/If8vBh63a/LLa/7ik+sKtTX/PLMlCn4Lmh1wM+z1YH0sdp/1ifBFep+P+Z0T6bn5+qgP+6vPDpo6wZOppfXNenT3sJKJrnP1C2854ht0/vGnK3Lpd+55Kyel3zCXfZ8dnEmWZ5LsdM96UQ/SVp76Xh6PdXcGQi3Ei5CAAAAgA4QJTWJ6AUAAACAzcdFCAAAAIC24nEsAAAAoCOEbfPB9O2xlQAAAAA6BndCAAAAgA4QJTXj9vhgOhchAFqKx83F8IZKZcN9aB708aOlysYjf1vVPHPO1q09/rZkeflUOhpXkorTC7YuDu+48Y6ta3b5G9hre9JRp4V607apTPoY1NrBwRvu17dMP77Lr2vOx766CNdcxOjlx/y82PlCul1pwfdh7oift4V6Oq64vOijYht9RVvXO75k6ybe25/uQyYCt2vWR/66GOHVQf8WoFD3fXdRxpMP+/HrvWyrWpLb3hwXxds9uWzbNMs9tm5hb3p55Vk/flfv8/HhA8+lY7tz/XNRu5LU88xpW1du4byaO5/ZePNDe22TMHHF1rnXmLjiY8Cz8fC+GbYIj2MBAAAAaCvuhAAAAAAdYm2b3CPYHlsJAAAAoGNwEQIAAACgrXgcCwAAAOgAUYF0LAA3X6sJU7l0EKe56tNkmrWarSvtHUv3oTZv24RqOt0nJ8xntqmcOVXV05FBhczyVu7MJDiZ5UnSSm+6H6vvOGDb5PRcSI/h3BGfxjNw0qcqzR/sTZZXT1y1beqDfl27vpxO/Wr2+7nZPOjndLPsX1gnH0nfmB/s9f2rTNsqu65apn8NvyoNXEynas0d9nOze9onZ538aJ+t2/fljSd7NczclKSVgfTYDpz0x0guKcy9dRg55s8xC2PpdDHJ76sL3+nPIyOv+OPUHVeSVDDjNPmIX1f1vF9Xoyfd94vvH7RtRr/pE+nmHh5NlrvULEmqrGZi0/rT54SszDmwvtsncZVfSScLZt9OZ/oXMq9ZTiuvjdg6XIQAAAAAHaK5TT4tsT22EgAAAEDH4CIEAAAAQFvxOBYAAADQAWKU1rbJB9O5EwIAAACgrbgTAgAAAHQIInoBbEgubtfF4+Zic0tju21dmPcRj06h4vuXi9ttjE+k+/Cut9s2xYsztq45nI54zEXqXv4+H4E7+tWp9Hr6MjGymajY+p1+3JdG0zGjxVUfxdo15+Mu1/rS+6RytWnbTD3gIy1HXkrPi9pdO3wfuvxYLOxNRxlX5nxU7PKwX16zXLR1o99Mb3MuXnjhDh8X2juejjJeGvFRrM3MK6KbM4W6b5Mbp91f9w8iuCje89/t53Rp2fejbBK4l3f62NxVP7Q2enjuTn+OWfKp2KqeSS9v/x/N2Ta588X0435l7lgtLfljOHe+cHMmF8ObOw4KjXQ/VvYP2jarg5mYaHOOkaTSbPoYaWRiu8vn0ufblmXigF3Ue+71CrcWHscCAAAA0FbcCQEAAAA6wLVvTN8e9wi2x1YCAAAA6BjcCQEAAAA6xJq2xwfTuRMCAAAAoK24CAEAAADQVjyOBSS4uN1Q8RGZazWTgympNDKSLI+ZiN7GmXN+eXvH0hV9Pr5VCz4yMuzxkZal4cF0xYkLfl2ZOOC6i3/MxELufGbW1kUTd9kst/Y3lkLdx+P2TKbzWHPRpKUl34/CXHpd3ZM+frS04E/bLvozF8PrtknykcRT9/k+7P+jWb+8fT4e98rb08usZaJJu/whp8uPpNe180W/f5eH/L5aGkmPbc+UX97c4Ux0qom5zdnzdR9nmltXoztdvjTit3fodb8uNxZ9l3wkcaHhY2mHX0hH8c6+bcC2ycUf5+Z011z6nDtzb59tU6j7fTV0PH2sLo/614rc8b2yI30uceWSNPTqgq0rLPh1hZX0Pi4u+NelWPWvMXHiUroPVZ/33Nw9ZOsKZ9PLy73W3g6its/3hHAnBAAAAEBbXfdOSAihW9KfSKqs//6nY4z/8i2/8yOS/g9J3/rT6M/GGH9+c7sKAAAA3M62T0TvjTyOtSLpAzHG+RBCWdJTIYTPxxi/+pbf+/UY449tfhcBAAAA3E6uexESY4yS5tf/s7z+s/GHWQEAAABAN/iZkBBCMYTwnKTLkr4QY/xa4tf+WgjhhRDCp0MIB8xyPhFCeDqE8HRd/sNSAAAAwHbUVNiyn3a6oXSsGOOapIdDCIOSfjuE8ECM8aU3/crvSvrVGONKCOHvSfqUpA8klvNJSZ+UpB1hmLspuOlcypUkNTPJVK6utDOdciVJhZWNX1iHTP9sApaktQOj6eWt+kSbQmZ7YyWTuDScTo0Jqz7xZGXUJ6j0nJ5Jltd3++UV5v3Yuna5dJqcRp9P6ar3pv9uUz2fG3efntToTadPrQ76/ZFL6ikt+n3suAQsSWr0pF+Qei/7Prz+8R22buBV/3evptnknkwQTt/FTILTrvQCrzzo+1BJT01Jftxr+/3yculdOXNH0nM3c8hpOX1KkCTtOGXa+GAim1YmSbueTs+z1QHfZuhVn843/sRgstzNCUlaNPtXkqrn/THX7Ervr4KfSlnNcvoYcecKSTJhZVmVq36bcuf91VGfSNfoSyeW9b6RTiuTpFDLpCya18c47M8JuXN7s4UUrNxrPn/77jwbiuiNMc6GEL4k6YOSXnpT+dSbfu3nJf3vm9M9AAAAYHuIUVojoveaEMLo+h0QhRB6JH2fpFff8jtv/pPthyW9sol9BAAAAHAbuZE7IWOSPhVCKOraRctvxBh/L4Tw05KejjF+VtI/DCF8WFJD0rSkH7lZHQYAAABwa7uRdKwXJD2SKP+pN/37JyX95OZ2DQAAANhetsv3hGyPrQQAAADQMbgIAQAAANBWG0rHAjbCReXlonFbWV6ritVM3qWz4vteGNtt62I1HVmbiztUZntdJGMu7jAnrGTyKfvScaFrfb5/3V85butWHrsnWV6or/k+ZOZMozd9Guua9dtUO9jaXMrFZDouBlOSVnZs7t+BLr8jPc+qZ3ykrosYlaSFvenyyoxvc/enfLyni2KVpF3PpvfXyoAfv/Pf68dv9BvpbS6u+rFY6/LbNb83XVdatk2yEbi9dV/nonhXMsvrnvR1zsAbfj7nxsIdP7lI3ZUBH9vtInX7z/rz49SDPnq23pvp+/50J0eOZSLbTayv5KObc2PhzlmSVD2b7kfX5HyyXJIW7xjwK8voe/ZcsjwXqVs/4GPqy+emkuXxtdO+ExUfpV4wkb+FzOtwzL2/uEUieqOCmqRjAQAAAMDm404IAAAA0CHa/c3lW4U7IQAAAADaiosQAAAAAG3F41gAAABAB4jStvlgOhch+K82O83KtSuN+HSNxlQ6XUOSQiZFw7ap+gSVXNKV+k2SS4sJXWH6arqinDkEM+Me6ul29d0tJH4pn9bSPV5Llq/u9euqffiBDfehembJ1q3t8bFALn1qYbdPVcoZOu4jVFpN1XJ6pjKJYEYuLWrg5MaTely6jyQVGukXwtzyZt/mk3pcCpIkTd2XXmghE9zWcyGXZpUuLy37NiU/BdV7OT1OjR6/vN6rfmxz7erm0Nr/JT83547486ObZzN3tXaMdM+ky5d2+TY9l31d11x6J+cSsFyCmCR1T/u60efSYzj1gB+/PU/N2jqXEpibt6sD/gByYxHLfl/1feWkX5l7LZO0dmA0WV6cXrBtmpl+uBTIQtcB3+ainxixlk4Eiyv+OGj1/Qq2BhchAAAAQIfgG9MBAAAA4CbgIgQAAABAW/E4FgAAANAJIt+YDgAAAAA3BXdCAAAAgA4QtX2+MZ2LkNuQi9q9HhuB22LknY38raUjX3Ntrie00C5mtmv1zjuS5ZXzs7ZNs99HPLZyyzFW/OG5bOJxXZxurs31LN6RjlztfWPOtuke3/jyclYH/P5dGklHRlbP+4zMwqqPih3/Dr8fKyaaNBftWlryMa0zd6f3cZffjdnlLe5Kv3Dlltfo8XW+D75ufq9/8VzORLjGcnqfNAfrtk3fMb+vlnelx6n7su/f4vt9NGn5+b5keaPbNlH1rF/XgokQzsnNzbz0MeKidqV8/0rj6bHd81W/r+bu9Mewi77uu5g5hus+KnZ1hx/3yYfTY5ib0/VBf5C4aOTKnD/HdF/x4+TGorTox2/l0cO2Lif3euFUTmWylucXN7y8MDxo6xpnziXLi1X/WlbMRfmbpHxsHR7HAgAAANBW3AkBAAAAOgQfTAcAAACAm4A7IQAAAEAHiOJOCAAAAADcFNwJ6XC5tKimSXdy5ZJUGhmxdY2pqWR5LoliLZN0pfvvShaHE+nECymflJFN3qiYcSr7Kb74XffZup4L8+mKBd+HlcNDfnmX0jE0tXcdtG16x31cS8/xS+mKzHwp1H2qUrPs/+riUrBWR/ttm/Ks73vXbHp+Lo5lYoYymmYX57ap0ePnRfWMH6eCCeuZOerXVWj4urI5fNw2SVKjxy/PJfwsZVKphl9Zs3ULu30CkZNL4uqeyYyTCQzqmfJ9uPJ2v6617vR+rFczc/2pdAKWJM0fTC9v7M99CtLMXb7vjR2+XfVk+u+DuXnh5qbkk66KmUSovkzCXc9kemflErByc9Bt19U7/QYf+pw/F5//QK+tWzOnmf6ztomWd5ZtXc9U+vjpmvU7pNHn58XAyfROaZZ9m8rkxlOpJJ9U2PfihG9Uz0w09zqck1leae9YsjzWzOuzrvOeBB2HixAAAACgQ9wKj2OFEIqSnpZ0Icb4oRDCYUm/JmlE0jOSfjjGmP2OBx7HAgAAALARPy7plTf997+W9G9ijHdJmpH0o9dbABchAAAAQAeICmrGrfu5ESGE/ZJ+SNLPr/93kPQBSZ9e/5VPSfrL11sOFyEAAAAAJGlnCOHpN/18IvE7/1bSP5P0rQ+3jUiajTF+60M+5yXtu96K+EwIAAAAAEm6EmN8zFWGED4k6XKM8ZkQwhPfzoq4CAEAAAA6RFMd/cH090n6cAjhByV1S9oh6d9JGgwhlNbvhuyXdOF6C+IipE1yUbuhUrF1ubg5t8xCJlK32cLy4sqKbZOL/NVZEyObi+FtNf6vLx3J2Oz3Y9v3rI8KjsM7kuVrB0Ztm+7xTDSg6V82xrHXx0IuP5iOLnTxt1I+ZnJlwD+Z2XM63ceFvX5/NI74cXf6LmaiGhcyMbJ70uvKtWl2+e1d6/JRmMXVdExrJRM92z3tI38X9qXbjT7nj7k3fjQT7fq1nmT54Ou+zdKI396Zd6f7sf8z/qXj0jv98kLm8F7dZfZXt9+Pu7+Q2VdL6X287A9hzd3jx6l3Ir286fs2HmMsScUxn4+7mD68pZM+QjinezJdnhuLpX1+LFar6WN/LT39JEmVdEr5tXWZfgwe98fO1AM+hjcXL1yZS29X11zm/LPo6ybek+5H90wmUtf0QZJKC+l2uVjfynl/vqjv9u8HuieXk+Vre3zcfPGi35Gxmh6LsOLHL1687Osy7z2wdWKMPynpJyVp/U7IP40x/q0Qwn+R9FFdS8j6uKTPXG9ZfCYEAAAAwLfjn0v6xyGEE7r2GZFfuF4D7oQAAAAAnSDeGt8TIkkxxiclPbn+71OSHt9Ie+6EAAAAAGgr7oQAAAAAHSDq1rkT8u3iTggAAACAttr2d0JyqVVOc9UnELnl5dqoheXllpm7sswlcQW3rlwq1Upmu/pNekkuASunnJmuC4vJ4lDxbep37t5wFwp1n2rSGPTRMI3edD9yqSvNsk9DWRlI160M+D5U5nzKUL3X/9Vl+vFdyfK+cb/vVwf9uE+8Nz1DK3Ot/U3EJet0XfVpYMMvzNm6qfsGbV2hkR6n3ss+xWe16sc21NPlZ37QH3OlU355KybUZmXIj+3yaCYR6rX0+eL8B3ybrkwK0spes8GSivPpOT30tJ9Lk+/0/XBjm1Nc9mNbmU6Xrwz75S3v8v3reXrjSVfRT2n1+JAhzTyQ7sfoN/z2Nnr8nGm404wP/FIhsz/2fD19HnTnOUla9aFP6s7Mwbk70tvVVcu81pZ93ejz6fPg6oCft/1n069XkhRW02PRNenP30uHfZpVZdKvyynMZ9K2DvhEzPIrJnEy8x7Cvu/I1OVSPnFr2fYXIQAAAECn4HEsAAAAALgJuBMCAAAAdICowJ0QAAAAALgZuAgBAAAA0FY8jgUAAAB0iLhNHsfiIsTIReoWqz4bMK74aDundPcRW9d4/aRvt3csXdFnonElBRNlK0mxNp9u46J2r8dF8Wb61+zPRAifMPF/ksLwYLp8+qptU675sai9Y6+tc7on/b5vdqVvOjbrPoKydtBHF/ZNpPMuc7GQ2bjLHf6EVz2fjoacesDvq6Hj/vhZ606vq7Y/d2PW11XPpssXd/ltWn3/oK1b8ynHGnk2Padn/x8+MvLqjJ/v1efTY9gY8tHNxSWf07p0MD0vui77Nt2TfmxdzG2uTSUTj7qa2Sehnq7LRbE2/ZRW7E+P4Y5jfiyamVdEF8WbiwLuO+PH6er9vmH/8UwWr1Fa8jHRTm57awd93VpPel0Fsw8lqd8cp5J0/nvT4zTwql9eVyal9eqdvm7fkybavu7Hr9HnJ9rCWHpf1Q5lorQH+m1d38WNR9jn4pSbXX5drt3A1y/YNsWFTCy/eR3OykTvNycubXx5uKVwEQIAAAB0iKa2x50QPhMCAAAAoK24CAEAAADQVjyOBQAAAHSAGPnGdAAAAAC4KbbFnZBCl08Z2vR17RxJl2faNM/41KfSoQO+4Xw63cmlXElSyI3FIZMINXHFt8kkZ8XqxlO1lvb5JI/S4D0bXt7qoJ/ipYV06pMkVV/KbLNRe2CnrVsaSaerLA9l+rfs11W5ml5eoeETXmbv9rOwnEmaufTO9Lp6J3ybXEpX9VS6H0Ov+1SYy+/wy1seTZeHTMhMIVO31u3rJt6b7seO3x2ybfp8lRb3pvdXNs3qIR8/VfpaemUrQ35e9E74v7gtmgC+td6mbVOv+uUVaz5laK2aPh6v3u/7Xpz3yxt4Oj2GC5ngu9y8cMfISmb/Nu5csnVdZ30Mm+tHPZPcNn1/7i+n6TGcO+Lb9I37pc3dm15e/0m/vEam705uX8XMu5cdp3zd5cfSr4H7v+hfNxf2+tfNnqn0vB1+wS8v9zrXLKfHsGvWT87Cqj8eu8f9yb3Zl07nW7lzl21TWszEwSmdEJlNvaxn0rHOpJO4cgmlyiSb3kq2S0Qvd0IAAAAAtBUXIQAAAADaals8jgUAAAB0vsAH0wEAAADgZuAiBAAAAEBb8TgWAAAA0CG2SzrWbXMRkovhbbYQ2ZaNgMtw8bhheNC2KeTWZWJ4c3IxvM2aj+sLZpyii+6VFKav+nWZ+L/CQjrGT8rH5raie9Kvq9Hrp//acF+yfOpBH63Yd9FHKI68mJ4Xy6OZPNiMRk/6JmbXnO/Dzhd81KmLhZSkQj0dgzp03I/tqb/qI2YHXkuXzx3O7I/MMLko3sYOH1uZuwm8ustHUO78anq7Zp7wUaw9z/ts0nItPe4Vn8IrTfhM2MX704NR6PfbNDeauSFeS29vYcjv+7V5v++7Lvt9XFxK163u8nO6OObHfe1w+ny2enbAtuma8WMRnkjvlNUJf/7uOpXb97ZKdbPItfsXbJvep9LnLEnqnk7Ps5mjvg+rmZel4mJ6nK4e8eeYyow/x1RPppfXaO30qGbmnU3P5XT5/EEfKd8zmTl+7ky/3jZ6dtg2A1+/YOsWHkznYl96zMfc7n7aH4/zdw/auuqz6RzmMNza+x8Xj1vwacX5ryjYmx6L3NcQ4NZy21yEAAAAALeyKL4xHQAAAABuCi5CAAAAALQVj2MBAAAAnSBK0X+86rbCnRAAAAAAbXXL3QlxKVitJGBJPgUrrvi0idy6SiMj6TYTlzbWsXXZ5Kz+dJpHrPqUj0LFJ2ctPbAvWd497mNc1g6M2rpGr0nJceWSKudnbd3K/kG/rr50glPPBZ+eUz43Zetq70gnguVSVwqrPo1pdSA97isD6X5L0mrVfzCt71I6Rax2MLN/d9kqBR/+otJyujyX1tJ31i9vziQ47f5TPxa1TOqOS+rJWevxywsNP+52DK/4sVi4359LCqX0nCm87FOVcgaOpSfo3KM+dc71QZKKtfTYDt6ZSafx4V2aWk6fHyWp//Bcsrxx2qdZ1U2iliStTaTHMDdbcklcH9h7Oln+xWPvaGl5pWXf9+qZ9PycGvbzYiUz7lJ6TrtjW5IWDvp5MfBaehSXRv2xU/SnYhXM+ac706Z20NeNPu9fo5dG068/c3f4mdF72b9mFcwuziUVru3J7qykg7+Snn+StHLUJ1j2vz5r65omBWutz7+OlC9lYt1Mmmdc9X0oHDrglzedbpd7f3a7aJpj9nbDnRAAAAAAbcVFCAAAAIC2uuUexwIAAABuR1Hb5xvTr3snJITQHUL4egjh+RDCsRDC/5r4nUoI4ddDCCdCCF8LIdxxU3oLAAAA4JZ3I3dCViR9IMY4H0IoS3oqhPD5GONX3/Q7PyppJsZ4VwjhY5L+taS/cRP6CwAAANymAt+Y/i3xmm9FoZTXf94a3fERSZ9a//enJX1PCGF7jCAAAACADbmhz4SEEIqSnpF0l6T/EGP82lt+ZZ+kc5IUY2yEEOYkjUi68pblfELSJySpWz5GNtuXSjoKs2jKpXycWytRb6VMpJyLvGt+08fButhhSVImUtdF8YYVHw2Y03N6Jlne7Pdj2yz769hCPR0LWvz6y7bNyvse3PDyJKl7Mp3xWLtrh1/XwKCtG3o1HTXYO+5jK2fu7bN1fRPp/uUif0tLPka23pse99VMorOL/ZSk8qLfrisPpqNzh1/xbWbv9vOiWEsvb+ao/7tF9eTG/6Yxd4+vKy/65e14zu+UmYfTx1Zh2W9vbPi60tn0sbV4jz8vHTlw2dadfH1sw31Yy9SNPXoxWT4978/fS+P9tq485jNX52vdyfLDD16wbU6+nI4Vl6TKgXSM8Mo53799h6/Yui88+UiyfG3Yn297z/q51EhvriRp9gfT55+/de8zts1/+c3vtHVdH5xM9+FJH7HePennRfd0+lxSqPvjavoxnwMeltLnhN1fzR33vm7qfv+6WTfnyNz50Z2/JanZlR4nFxt/je9fs5zerlwcftekj8yefLePxR55Md1uedS/5tcO+ZjowZfT7cKJc7aNcl+vMDyYLG5Otfh+6vZP9r3l3FA6VoxxLcb4sKT9kh4PITzQyspijJ+MMT4WY3ysLD/JAQAAgO0oxq37aacNRfTGGGclfUnSB99SdUHSAUkKIZQkDUjyl6oAAAAAtq0bSccaDSEMrv+7R9L3SXr1Lb/2WUkfX//3RyX9cYztvp4CAAAAbm0xhi37aacb+UzImKRPrX8upCDpN2KMvxdC+GlJT8cYPyvpFyT9pxDCCUnTkj5203oMAAAA4JZ23YuQGOMLkv7CJ/JijD/1pn8vS/ofN7drAAAAAG5Ht9w3pq/Vasny0t50Ksxmr0eSmmd80kNxOh29Ubr7iF/ZQjoJ5XpCzbSb98tbevSwrev+yvH0eg7ttW1ySRl94+nUi2YmAWt10E/JZqls67qvpNNLFvb4Jw6HjvtUDpcOsjTiE0/6LvqUnEJ9408nDr26YOvm7vJJXM7iLn+btVny27W8N71dtZrfV2V/+Gh5V3osBl7ybaYzURjlWnq7ho75fd/0U0nz3+3Hfd9gesOu/v4e22Zxrx/30nK6vDHrOzheHbB1hf70cVA65Y/Txg6fcjb5tfR25favjvjj4J7dPtnr1FQ6xeeNZ/f7Vb3jvK2bnE+nYC32+5S92vLGA1O6x/1xsHBnJhGqkXns4XJ6f/36cZ+A1bjHJ4/pN3Ymi+/4O6dtk/O/6V8r3LlkZch3oeesn9MFM2UWd/nl7XzRz9uVgUxqo1mXK5ekmaM+ccklheUSB2sH/fJcKmIufbHQ8HW7vnzJ1tV3p9+v9L3hD/DY5ed7Yd7ET5mUK0mK07O2LuSSroxmLm3rFnHtA+Lb41suNvTBdAAAAAD4dt1yd0IAAACA2xXfmA4AAAAANwEXIQAAAADaisexAAAAgA6xXb5pjzshAAAAANrqtrkT0hifsHXFajqGTpLiSjpSrjSSjouUpMbUlO+Ii7NdyeT/5axk4ubK6d1Xv++AbdI97qP3QjUdabkymi6XpOKqv1xf3pmOZCw0fJuF3T4qtmfKRx46Lj5Rksbf7+P/Bk6m2zUzR8zM3b6yXk3XFTOpmpKPXWz0mFja1/08q+3PRCvmpufJdLvey35sc3HA1ZPpv31Urvro1ELdz4tVE026dNBvby4utPurftwvHE1HuPb6BFwV6n4sGt2+nfPA7ou27o25dEbq5K5MJnFG9+G5ZPnVCX9OdTHBknTsVX9uUt38TWzUnwMnvuiX58Z2cNJ3YWnUZ8ze8b50NPvJ1308fHHez9viov8bYLfp4+Jef8w1G355lavpc+fJP/YxvMv3+PNtcTk9p4tLfq6v5ea6iaoOmfPtpcf9uvrP+nYDp9Mnu4U9fmV7f9/H3F7+rt3J8sqc31c7n5m1dVceHUyWj7w4b9sUp32seLPfx04X6ul9nIvhza3LaWS+1qCQieFt5f1ZzEX0XvVV2Bq3zUUIAAAAcKvje0IAAAAA4CbgTggAAADQAaICd0IAAAAA4GbgIgQAAABAW902j2PlErBc6lNOLmEhl8zQfO10ug9j6QQNSTbl6rp19XTKR/HrL9smq+970NY196bHsOe4TwYpz/pxXxntTZYX6j41JFfX+0Y6qUeS5u8eTJb3TPqknuKqH9ulkY1fn9f9UKjncrp8NdOmdmjjt2NnHvZ1A8d83dKuDa/KJnRJ0lomLar7fenon+k/G7Vtumb88haG08lPPeN+/y6P+uSftW6/78vT6WUujm08uU2SYo9PBHOeecMnQh0eu5Is7zqcLpekiy/4c9PV7vQx7MZBkurdfpv6TNKaJN3xQ+lz57E3TOKgpMV70uk5klSYTc+L9/zAc7bNl8/eZevGZweS5bltKr7fT9z37k1vryT94fH70ut6xh9YazM+BenSO9PlRZNKJUk9e30aU+P4jnT5nT7ur+d53/eF+9P7seus36ZDn/Ov0af+hk8lGzqePlaHXl20bVwCluQTE6+83c+LXavp40qSKnPp/k0+4t/HDB336Xddk34/LtyRfgGqvuTPF9m0rXmTZrXXJ8jF4fRckiSdGU8W5xJKc2lbt5Jt8jUh3AkBAAAA0F63zZ0QAAAA4JYWiegFAAAAgJuCixAAAAAAbcXjWAAAAECn2CafTOdOCAAAAIC22rI7ISGElqLUQsXHwzmN8Qlbl4uOa0VhOB2fGSt+qGPZxwkWLvmIx+buoWT5wrsO2jb9r8/aOqd+wEcSL45127odr6T77uJ0JWllh78uXhoZtnU9UxuPOr3yoF/X0PGN/xmi0PAfJGuYdMqBN3y06+zdmajYmqvJRKdm4oDrVb+9wy+lyyff48c8F+G6/EfpKN7lQ34sui/7seg2UbxrfmpKVR/dvDbk+9H1WnpHhobvXy62tDmfjtbMxb4uHPHz7I1n96f7MJSO85aku95x3tZNfDF9PlsZ8vOl3OPXtfSQH9sDfenzxbFlH0ncM+pjVZdq6bH94qmjtk2p5Of06ilzAO3yY9GY9LGqvz9zv60rTqZfG2sP+Uji7lP+tXGtNz3u5Xk/b3s+76NTZ4+mt7n7ZR/DmzvXlZbTfc+dsyYf9ts79JxvtzSa7nuz5GNua/4lVQf/MD0HT/wNfwJqdvlxLzTS/Rt5yc/13Ouw5Oegi8RfOpx+byHlI/tjNR09HDJx+GH6qq1r1NIvdKVD/pzQOHPO1t1K+GA6AAAAANwEXIQAAAAAaCs+mA4AAAB0iMgH0wEAAABg83EnBAAAAOgAUdvng+lbdhESY1RzdXXD7YotpGNlE7BWTB8qG0/uysklQDQPpNOCJPn+SSpMp5MjCvt8GkYuiSvU08kwy6N+zLvmfBKOS8HqmvVtumZtlS495vuxuCs9lYdO+LSbnS/6tJaVgfRNwmbmiMmll5z9/nRqSLPkb0bm1jV/MH2vtrjkT1y55XVf9u3mjqTLcwlOqz5cRQsmBSvUWzvpHvnA6WT5sVd9gsq+MZ86NzE5sOF1vfLcIdtm4CmfGLS0K12+cNDPzZzCwfQcfGD3ZdvmpRPpRC1J6n00fd6qPOOTk5bH/Jyu9PhUsj/8s4eT5QMn/fJWn7BV+qff838ny//d7/wl22Zpr0+fcmfOQmbe9g75ZLTFGT8v4thyumLGnwOX7/R9730t3W7se3yS0InTe2xdZTydJOVT+6SL7/N1PRfMeqZ9m5y+S/687xIYq2dz70f8+4G5I+n9OPCanxdds2b/ZpbXPWmbaK3Lr6vywhu2rnxwd7K8MO/nkkvAkqRmX3qehdfS500pn3jqUrDi9Kxtk01d9ZuFLcLjWAAAAADaisexAAAAgE4QJW2Tx7G4EwIAAACgrbgTAgAAAHQIInoBAAAA4CbgIgQAAABAW23Z41ghBBullovujSvpjLVczFtOs2YyBQ/e5RsdO2GrwrCJrqz4SMvCgs+NW3iPyUeV1HNhPlne+8acbZOL6G0MpqMBe8d9zGSjNx3VKEnNUvqDVQt7fYRez6SP8Oy97O9PNno2/iGueq9vU6in1zX0qh+L1UG/XY0d6cjVyoz/O8DASb+9c0fSfc9FZFbP+9jX8e/JRMJ2m7jLs/6Y2/XoRVt3/tyIX5dRmfHzzEXMdk3709uFusnGlVSu+X1y+uXDyfLmQR87XfpLU7buQ3tPJsu/cvkO2+bChM8/rl9OH8OnSn7M9x/w/btweme6YszPl3LZ160d8+fBw+9Lx8Weu+qjlg8M+nPdb088nCxf3eXPMf3H/JwumF3sIqclScf99vb7lGg1350+zzTO+ljf5vLGY+UnvujHNmT2cTSHVsmfHjXwqj+uFvamy7sy57Oumj8/uhheSVoeTp87Gz2ZqNiljZ+Lh477Notj3baubyI9P3OvL7nXTfX7SN0lE+ffPenPnWHVn+us+/37qXDJHwiNM+lzQi6Gt5WvfuhIPI4FAAAAAJuPixAAAAAAbUU6FgAAANARgiLfEwIAAAAAm487IQAAAECn2CYfTN+yi5AYo00xyCUfuBQsl5olSSGzvMLY7nTFvF+eDvlEkXhmfMN9yOlb8UkU9d3VDS+v0et3uUuEyrUpLfr+rZltrsyZtKXrWNjnb08Ovp5Ocikt+HXV9vuksMUxU77LJ4249BxJ2v9H6cq5w/5mZHHVn4WKSyZ57FEfT7O41yeyFPp9usrIF9LtaodsE118wRxXkorl9HYVF/1YDH2XT9uqP7MnXb7XH8Ojo1dt3fyyT8k5PJJOknrlfLoPkjQ56ROSvtm9L1meS8AqlDJpTEPpbV5Z8uliFyZNApakWErvqyceftm2Ga2kU/sk6dNLj9i6s1eGk+X1Xf7AmpxPp/tI0pJJi6qM+7G444dO27pjb6QjnHL7I3emu+O9fk6//ifpFDYdXbBt+rp9KtDqM+n5tHTQH/ddl/04OVPv98dcXyZ5rHp2w6vS9P3+9aArkzzmUhZzCYvVs7nEpfQ8yyVWNbv8ua7R51+XnEKme2t7/Lmk78WJZHmcnrVtQtUfc7GWPvbjXf49k1Z850sj6VS/mEnAKlUz74uu+CpsDR7HAgAAANBWPI4FAAAAdIIoPpgOAAAAADcDd0IAAACATrFNPpjOnRAAAAAAbcVFCAAAAIC26sjHsQqZiLVmrZZuszMd5Sb52DhJCq6q4mMhs8sz8XUrR9PxjpJUWvRRfjnFhXRMXaj7YMjyJb9d04/vSpb3jfs4vOVRH7vYNLOr0eOvfSff66dkfybGsdBI37ucfMhHI7uoRkkq1NMfCque9+M3d9j33dXN3e+XN38wF9WY7nv3yz22xcKdfp4NVJdtXe1QepnLd7YWgTs1kz5GBof8cTXzZR+B22U2q2vGz83J+31s7r4xn+/5ynPpXOKx+y7bNrVM5K+LpVXNx6M2q34/9lQ3HtFbWPbH430Pn0mWvzzlI5inZo7YunKPn++rtfSxGrr9+WzhtUFbt1ZNtxt9dNK2Ofaiz50u19Lj1O0Xp8b7/XFw7FUfW1q4Mx21nTsjLB3zUax1c6zu+pI/P17+Dj/PdhxLz6e1Xj/P6pnk1DVz2qpX/Tl68Lj/4K6Lm5ekgZPpsZ16wMevL+z149Qsp/tRmVy0bdb6/PIWx9KR6P1n/fKWR338elj1x1xzOL1T1g7491OFuo+kLk6bMTxxzrZRJvK3eSUdie6+qkHy7xFvPXwwHQAAAAA2XUfeCQEAAAC2JT6YDgAAAACbj4sQAAAAAG3F41gAAABAp+BxLAAAAADYfB15JyQXseai2XKxuVn96Ui55sQl3+b+u2xVmE9HIXZN+v7Fsg9enH3bgK2rnklHDebk4vUqc+lIy8rxcdtm6q8dtnXV8+nlrezw1767v+HjOOu9vp1bZv+4/3OCixCWpAWTqFxotHbIdF1N96NY8/t+rdv3fehYentnHvZxjMVuX3f1rJ9nXWaTu0/5mMSF10Zt3dh3XUyWXzi907bRQd/33rPpDh7+/tO2zWuX0nHUknRhwked7jg8lyw/f85HWv7AQy/aui+eOpos33ePj/zNrWvlXDrust/0W5KuNvxx1VtKx7TOZ2KH/58P/6mt+4NL99m6ye503+sNf4w0MpG/a0vpeTH7vJ+b5SMLtm7wznRE6uQpvz/Kx30U9OhDPtv38sRgsrzrtUzUsh8K1U3k846/fd734bSPxXaRus1MhnDFJ1/b5e37sn89yEWid2deGl3keFfNn2/da+M16Y1uJYY3p1n2x2lpwfcvuhO4fHxv+ZR//7N2wB8/sZJeVxgetG1y3Nc15N4jNlf9VwrcMqKkSEQvAAAAAGw6LkIAAAAAtFVHPo4FAAAAbEeRD6YDAAAAwObjIgQAAABAW3Xk41guEUGSokk+CF0+iUK5ZIaFdOKJS+GSpHApE/PRl07bCrX0eiSpccCnqwy+7FNtXPpGM5O21ej1u3zm7nRd3/H0Nkk+9SlnedinPuT63kqaVU6PD6fRzhfTKWK1/f66Pde/Rrfb5sz4lX2S2cy706lFA8/4eVt7V7qNJKmaqZtJLzM8dNU2qWdSgabn0/Np9IA/rnJpTM3x9LpeOe/TfZqZRKjyuF/X0bvOJMufW/bnn2cuH7B1axPppJ4LZswlqbjs+374wQvJ8pPnfBpYzIzFM0+nkwDX+n0az8/90ffaupG7p2zd3GQ6Hat3yEcdPbwvvb2SNFJJJxJ+cfwdLS3vG99Mj0V5l+9fmOmzdZPnfArb4LH0yeQ9f/tZ2+bJ3/Xb1T2Z3scTX/RzM2QS6ZzeCT+Xakf98kaeTm9vo8cvr+7fJqhe9a8xg6+nz6tLI35dFf8yrPKiP0873Vf8+XZ1MD0Wy6P+nJBLx5o/6F+/eyeW0xW7/dwsXsy8/2lBnJ71dSvptNHbIgHrengcCwAAAAA2X0feCQEAAAC2Jb4nBAAAAAA2HxchAAAAANrquhchIYQDIYQvhRBeDiEcCyH8eOJ3ngghzIUQnlv/+amb010AAADg9hXi1v200418JqQh6Z/EGJ8NIVQlPRNC+EKM8eW3/N6fxhg/tPldBAAAALDVQgjdkv5EUkXXriM+HWP8lyGEw5J+TdKIpGck/XCMMRtldt2LkBjjhKSJ9X/XQgivSNon6a0XIZumWavZusJOE2e7ktlOE8MrSSqnhyDkYn1zccBmXWsHRm2T2qF0TKckFVe7bV3/2fS6Fsd8m1w04OhzJg6v30cDLu7yH56qmATXgTd8nGC919+cK676S/S+S+m63PLmjvi+V8+nl7fikwtVMmmHkrQ4lo5x7NmbjhGVpPCMj7ktLqXnbS62cuApP89Whn27ejo5VYszfnl9frNUKqX3fy6Gd2ncdEJSsZreV8VMxHEuojfn5c/ck67IzAs95M8/xbF0vKsPqvaxvpI0PjuQaWn60O2jU5ul9Bj2VtPnCknae9jnmZ6e2Gnruqobj9187sI+W9c8m44mXdvr+/7iH5r9Kz/PdNLH8D70gdds3fEp/5ow153u+xeefMS2ketfxo6T/hy4/Db/WtEspc8/a/6lR12X/duNylx6nk0+kolEL/vtPfwZH5vsYurn7vCv67moYGd10C+vtt+PxdDx9PzsmvQn1dVRf34sNPw4lS+Z91otRuDGqvmKgokrvo2J4ZX81zU0p3zU920hqtMjelckfSDGOB9CKEt6KoTweUn/WNK/iTH+Wgjh5yT9qKT/mFvQho6sEMIdkh6R9LVE9XtCCM+HED4fQrjftP9ECOHpEMLTdfmJBwAAAKCzxGu+dVVcXv+Jkj4g6dPr5Z+S9Jevt6wbvggJIfRL+k1J/yjG+Na/cT8r6VCM8SFJ/17S75iOfzLG+FiM8bGy/F89AQAAALTdzm/dMFj/+cRbfyGEUAwhPCfpsqQvSDopaTbG+K1b6+d17amprBv6npD12y2/KelXYoy/9db6N1+UxBg/F0L4/4YQdsYY/T04AAAAAG8Stvp7Qq7EGB/L/UKMcU3SwyGEQUm/LeneVlZ0I+lYQdIvSHolxvgz5nf2rP+eQgiPry/3Nn9oDwAAANieYoyzkr4k6T2SBkMI37q5sV/Sheu1v5E7Ie+T9MOSXly/9SJJ/0LSwfUO/Jykj0r6+yGEhqQlSR+LMXb2x2oAAACATtPB76BDCKOS6jHG2RBCj6Tvk/Svde1i5KO6lpD1cUmfud6ybiQd6ylJ2ftCMcaflfSz1+/6jbEJWJKaV9I3WELFf8YkVDJpVvV0Mkys+SSKxkNHbF1p1qRtrfoEmqHnp21ds89vV7O88cSOQt0nU7k0j0uPDdo2Qyf88rpm09u8vLNs29QO+ak2+Lo/Kl0KVm55Ayf98pZG0vlEjR0+cWnouF/X4j3ppJk7R/wNw8vv8+EN0y+nU4YKdd+HlSG/vY0hPz+DSUja8byfm7mUrqXl9Dw7uNMfBycySVzvf+fxZPnTv/mA70Qmzerh7/SJRiOV9HnhD4/fZ9tMnvLns9E70/t/VyZe7NjlQ7auvzs9ZxrHfdLaj//l37V1l+rptK1fefGdts1kySf1NOf9sf8dD7+aLH/ypaO2zcCoH6fK29KpZAt/5lOplnf547u4uPHz7alf8mlbVx/y6+qaSa8rl8DXzLyaF01Y1NUj/pyQS9NbNcf36nBmm6b9+NX2p+uGj/n+XXnYVmnqgXRKkyR11dLLLPhToOq9/rza6EnXjbzkd9bQgn/dnHogfV7tnvbvY6pnfBpY77hPOVsbTie7Fc9lEkX7/Ni6FKxGJs2qtHfM1rn3e4VMQmmzxWQvbMiYpE+FEIq69uTTb8QYfy+E8LKkXwsh/G+SvqlrT1Fl3dBnQgAAAABsbzHGF3QtKfet5ackPb6RZXERAgAAAHSKDn4cazO19o1dAAAAANAi7oQAAAAAnYI7IQAAAACw+bgIAQAAANBWHfk4VmN8wta5aLYwPGjbxOlZWxdc1NuhvbZN6fmTtm7xu9JRnX3HfUTd3MM+MrJ3wsf8NcvpGNnSko9JrB3ysYsu5q/3cno9klTKRA3WDqbHdnGXjzus+JTWLLfMXGzl3BHfDxc1WT3pr9sXd/l17RubSZYfe9HHrQ4d8+uK7zeRjGf9/i0u+e39K9/7jK37rS+nwy66v3fSttnV7eft6Yl0vHBOsdvnZ55fSMfIrmZieOvDfnm9JR/x+IUn/0IgiCSpOerbuBjenN6Sj9U88jb//U+T8+l43PpeH/f8qdPvtnUfP/zVZPnasj+wRsd8bO5cJr73xFx6Xjx+72nb5hvfvMv3w4x76V3pY1GSdDo9lyQ/Z3aM1WybqWGfVV2s+fNqNMObi9necdIf3zMPp/u+/w98m8uP+P5Vz7oaf87KxQu7SO+Zo75/O/zLsBb8y7dWltPL7Lvgx7Zy1b+mutdhF3l/PUPH0+eSyqSPzS1M+zmY5b6iYNhHeue+NqBQSU/cUub9mTLvz+zXNaz4820zEwd8y4ja6m9MbxvuhAAAAABoq468EwIAAABsR4EPpgMAAADA5uMiBAAAAEBb8TgWAAAA0Cm2yeNYHXkRUqz6RJFQNekqJuVBksIeH1tU351eV/lSJm0is7y+Z88ly5u7fVRPLs1qcazb1u14JZ3yUlr06RVds/7mV6M3PR0aPT6lodHnE1Rcu/5xf3TN7/Xr6qr5ut7L6WU2Mm36Lvo5M3N3eixqR/y+6jvjx/bCa+k503fWt1n7AZ/iUzyWnk8Pf+drts3zf3yPrfv0U+kELMnfLp163SSXSJpd9Ns18lA6VeuNZ/fbNs2qT2E7/eK+ZHlXJo3nnnvTx6kkffNSenmSVDiYTqjp6/ZpLbWv+fS7ogk5e3aXTxDb8/ZLtm5pOZ3I033KnxOmZv1+/HfzH0iWd1X99p48l4mJy6gtp/u4WCnbNgU/LTQ1k36tCBP+nNrs9uemQn86sezqTK9t0zPqE41Wa/51zqWZFWb9WDTL/lxXPZ4+n11Oh71Jknp8+J2mHkqfB4smeUqSKjO+zh0HwYfEadUPn3Z900+Mem/63JR77WnloZGFPf7tVfWsP35Ki/51abO59yVL+3yKXd+LPr3UJZHm3k/bhFJJsZZO2lur+fdnLkFVkuRDArFFeBwLAAAAQFtxEQIAAACgrbgIAQAAANBWHfmZEAAAAGA74ntCAAAAAOAm4CIEAAAAQFtt2eNYIQQbpZaLbNNKOtqumYls05UpW1WqHE6Wrw332TaFl0/busZDR5Lly6M+IrPvDd/3+mCPrZu/ezBZXqj7+3guelaS9j45myzvGvWxkKUFH4XYLKfjeytzvk3fuK974yO+H0Mvpa+nc2Mx8V4/Ft0mnjLUfYzj0Os+WnGtJ72upX0+8lc1HyU69ujFZPk3vnmXbZM5qhR7/LiPHUgfP5dnfUbmrkE/py+cTsfPFsp+X5WHfLZi/XL6GFm7f8G2+etjT9u639Bjtu6V2p5k+dUZPxb73pfeV5L04PB4svzPx9PnJUm6+MJuW3fHO84ny08c9Hs/F7eb249O906fq5qLYb7aMGM4MGvbNH1CuIomircx5I/Tuw77fdVK9PDSuI86VSYOuDiZ3l+5mNvZh/y47/rT9PlnZcgvb+5Rf8yNfindv9mjfnm5uF0bsZ6Jh8+ZeK//++qhz6Xne3nRvx7kXqOn3z6QLB95Ztq2ca/dktTsSve9Z8Efp81hf/6JXX673HuZvvnMXJ/3sdNheDBdkfkKBZV9/+JEetxLIz5WvDHl3+/dUmJrc/9Ww50QAAAAAG3FB9MBAACAThC1bb4xnTshAAAAANqKixAAAAAAbcXjWAAAAECn2CaPY23ZRUiMUc1Vk3SVSTdwqQiFQwf8ysx6JCmeSafT6G0+nSbs8ckR5VOXTI1PtFkZ7bV1C2OZRKhX0+k/jV7fZuANn7iwtC+d5FJc9UeDS/LIqe3PTTtfN/SSb3X1znR5/9nWEiYW9278DJDbrqar2ukTaNaW/fKm/iyd0pQJQlH3QzO2rvLkkK1b/dP03C192CfGnD/n00vcHiku+rm0WvPpTr1755Pl/+C+L9s2/9tzP9jSumySVL+P/rkw4cd2b9/VZPkn7n7KtvkPje+ydWevDCfLw5KPkaqX/KSZLqXPTUs1n/Y3OpreJkl6/3ccs3V/eiKdLHhqys+lffdctnXT8+m+v3/fWdvmyZeO2jpnYDQ9/yRJftfr6kQm0WgoPZ8Whvy+Ks34OpdaVc50vXDK7+NZM0z7nvSvtVP3++NqcVe6f2s+HFLlTHDb4HF/3l/emX59nLvDn3+6ZjMdMaYeTR+LklQ968eptJje96HuEwyX9/q51HPcvSeRZN7L1Hf75ZUqG3/bGE1SpiTp2ImNLy/zng63Fu6EAAAAAB2Cb0wHAAAAgJuAixAAAAAAbcXjWAAAAECn4HEsAAAAANh83AkBAAAAOsU2uRPSkRchhS4f5aeKqZuetU2aB308rltXYcFHp+asHRhNlpfP+djhxtFM/xp+Xcuj3cnyhT1+t1bmmn6BRvcVHz+6sNfvqxUTT1k944+u8qKvm3zE37jbcSpd3ncxM4CZ6b8ynI54rExnlvaXJm1d+Y/S86J5ykc/di/7dbnI38adS7bNaL/P4zxxfzqeWZJqjfRYFM75NplARoV6enk99/sI4YdH/Ng+/8f3JMs/WXq/bVO/7Md9fyb21cXtxoafm4/fe9rWjVTS++Rn/vCHbJtHH/ORlouNdPzoqZKPuW0c32Hrlvamy7vO+vjWqZKfFz27ztk6F0ndP+TnbS7++MiB9H788qt32zZdl328eW5+OpWSP/+MPJ2J4H5/Oo61Mu7719jhz+3F2sajynOvPT0m5bjR54/8kWM+VtVF0dcO+X53XfWvFY0e3653In1ire33UflzR/x8H3kxk3NsNMv+fBFW0wPfcgxv5v1UvJg+RkqmXJJC7v3Z8GC6zcQV32anPzfF2sbHFrcWHscCAAAA0FYdeScEAAAA2G5C5HtCAAAAAOCm4E4IAAAA0Cnixj/HdSviTggAAACAturIOyGFqk+B0Eo6YaNZq/nlmSQPSTZtK6z6oWkOZ/rXQptC3T/81zfuE0Wa5fSVcmnJL6+23193DryRTmQ59Vd9IsvI8/5qPZeC5czd4ftXH/ZxLY2J9P668F25nCbfv8LBxXTF8322zcqfpROwJGn50XRq1cBTPqVp5mG/vR9919PJ8t/77XfbNhMvH7B1T3z4JVv39QsHk+VLNZ8YUyj5pJ47xtJJKSdO77FtvnF2wNb1mhSxj9zxom3zx93pRC1Jujzrj9WPPvTNZPmDvT716f/z6vfauqszh9IV3X783pjziVBzz6TnYC456fC7fN+vLqcT+C43Bm2bu8z+laTPP/+grRsYTSfhTM34tC23PyTpzMJwsvxcJtlr54t+nC7sSvfjiQeO2zYn5nbaust32iqbPrZ60Kc2lsf9dpXNy2PJh+lpaZevc+emhcuZZMYZfy52/ahkAslySYpDry7YukZv+vVs7M99EtPiXn+edmlWK6M+bavnpQu2rrk7fXx3zfn3Amt7/Dkhl/QZTJpV/YBPrCpf8u+1cimlmymutJZeis7TkRchAAAAwLbEB9MBAAAAYPNxEQIAAACgrXgcCwAAAOgQfE8IAAAAANwE3AkBAAAAOsU2uRPSkRch2bjdQ+mY0UK/j8NT3Uedqi/TzljrS8f6SlJxIR2jl4vrWxr1EbirVR+BO3Aqva7yoo+ZrJ73y1vYnY5QvPO3fBzepcd8LGTv5XT59AO2iY2SlKSh5/x0nX6sniwffN6PbddVf5Q3j6ejeK9mYjXXev24d72WjnjMxfCWp/32fvqpx5Pl8WB6HCQpNPy+/9MTR2xdcya9j2NPOtJZksKkP0Y0Ztos+QjP3LqGvis90X75G+/xfcj4f733923dJ19/f7L8t197l23z/u84Zuueeu3+ZPk73vO6bfPMGz5qWXems07/8Tv+yDb5f3/9+21dbKRvlvcf98fVGzP7bV3uBWepauLSJ9IxwZL0h9Wjtm7+dDrWuXjUx7fO1XwEd8/Z9Fh8uXGvbZPjR1Da8/X0eWFu0p9vF81xJUmLezf3HU3fyfSe7B/367nysF/ezufS5bnXv4U9/kGOpREf6zzyUjp+ff6gf43ORcd3zabP7bno/dz7jsJ8+vU2Vvw5MKfZl4lSN+XlU5dsm7UDPoq+cNG86Oes+uhhp9lCG3QmHscCAAAA0FYdeScEAAAA2HYiH0wHAAAAgJuCOyEAAABAp+BOCAAAAABsvo68ExIqPs0htpC+ELoyST3l9BDE6VnfJLOuWEkvr/srx22bwmP32Lqesk8H6R5PR0kt3pFOhZHyiR0Fk4RTOT/r2zy829aVF9PrGvtz34e5O3xCUlfNtzv0O+lxuvBdvk1lxo+tS+lqlv3yRu6esnXLf5ROFFld9n8HyKVtjd6ZXtfkuSHbxiXaSNLqkK8beWhyw+tqDPnUr3N/lk536l22TbRwp0+Gmf2jPenlvfuqbbM07tNzcmlR3adMUtg96VSq6/mfvv9Pk+Vnl4Ztm2LZz4tdg+mJ+6vn3mnb/F/f8cu27u997keT5UsP+e199I5ztu65C/tsXfdX08lU4YkZ2+bqWX+ui/3pOePS3iSpfjCTcDeTPlYLmWM4J3d8T7w3fTzWh/1x1XvWH8Mr5lAtLm38HChJXaZu5qhfXu+4X16zlD6v5s75i7syyZGnM4mYxsqA34+jz/s0pvJs+lioD6ZTsyRJC+mELklaeEf6/Ng96U+QhQWfYJkTJq601G7D66n6823zin/dzL0XxO2hIy9CAAAAgG2Jx7EAAAAAYPNxJwQAAADoEET0AgAAAMBNwEUIAAAAgLbiIgQAAABAW22Lz4Q0D/oY2cLZS8nysGeXbbNwdMTWuRi9+Q/eb9t0zfk4wdUBv4sK9XTs3cIe32bZp6pq99PpmL+lw77RwBs+OvXSO9Nxu93pxFdJUjMzIytzfl0zd6cb9p/1y5u734/70HPp5Q0e97GQy+PpGF5Junp/PVnef9wHPsdMFvSU0nOwesb/XaFe9ctbO+gjVycnd6Qruv3+qJz10Yr/8KO/myz/mWe/17YJJj5akpominfv4Jxts9zvIzLPn/PH9/Kd6WOkcMVv79d7Dtq6E/07k+UXJjIH6rKPsb6wlJ63+8Z8zO3PT3yXrSsup+d78ZSPH33xNR85Xh/1sbQHP5SO9j39oo/17Tkwb+uclXM+LrR6MjPP7LnJt4mZ81numNPZ9PiWp/0Ci5nF9dfS+7GVKFvJRwgPH2stUtfFufdO+FjaypyP3m/0+H3S++WTyfK+wftsm9KiH6e1vnQ/uiYzc9N8NYAk9b04kV7PHn9OCDV/Pgtlf75QxYxhpn+5OOAwPJiumM/0LxPDu1ZLZ0EXMl+70Fz1ccroPNviIgQAAAC4JfDBdAAAAADYfFyEAAAAAGgrHscCAAAAOkHke0IAAAAA4Ka47p2QEMIBSb8sabeufVTmkzHGf/eW3wmS/p2kH5S0KOlHYozPbn53JR3am+7ndDohR5LCJZ8Mo/7eZHEuAau04FOBXHJE9xUfdVRaTCcnSVKzy18nrg6md9/Qqz6JYube9PZKUrOcTi9ZGs3ENLWglElxKWTCWnKJJ2smrKfrsl/ewDE//RtmeQuHfLpPMxNCMvh8egxzaWCLu/yfQiomFWh50SRZSarv9akm3a/5tKOlg+n5edfhi7bNyZpPNPrVc+9Mlld6/HEQn/fb9UN/5ZvJ8t/68uO2jXb6sfgH7/6Srfu5574jWe72hyQt1Xz6y8RSel7ETBpY/0k/aRbuT5+bJiYHbJvVhp+473jP68ny5/84k4BVzfwJr+yPn5Ovj6WbLPqxyI1t37F0Xf/7/etB/bJPIFo8mD45dY9nEqvS4T7X1uVOMhk9mWTBrqt+3Ce/O50Y1Oj249d72S/PpWD1TPpjuNHjE41KS+l5MfmIf72qnvevwzkzH34gWZ47F/e+4V+0wkp6XqwN9/k2mcSqwnz63FScXrBtVPcvnCujfgx73HujXJpVps6JLSZW5VKwcHu4kcexGpL+SYzx2RBCVdIzIYQvxBhfftPv/ICku9d/3iXpP67/PwAAAIAbxeNY18QYJ751VyPGWJP0iqS3/pnzI5J+OV7zVUmDIYT0n7UAAAAAbGsb+mB6COEOSY9I+tpbqvZJevM3TZ1fL/vvvnUnhPAJSZ+QpG75W4QAAADAtsSdkP9eCKFf0m9K+kcxRv8BjIwY4ydjjI/FGB8ryz+LCgAAAOD2dUMXISGEsq5dgPxKjPG3Er9yQdKBN/33/vUyAAAAAPjvXPciZD356hckvRJj/Bnza5+V9LfDNe+WNBdjnDC/CwAAAOAtgq59T8hW/bTTjXwm5H2SfljSiyGE59bL/oWkg5IUY/w5SZ/TtXjeE7oW0ft3rrfQEEJL8Ws2infFR8Ct3eXjQl3sXd/xKdumvrtq61xcn4u/laTl0W5bl+vHlffuSpav7PDRj6UlP8MWxtJxoZWrPlaztt9HDfaay9CldLclSQMnff9cjKMkxUfTEam1Yz7atXrW98NZHvNRjdWnfCTj/MH0djWGfLRiWPJj2zie3q76sF/ekQM+r/jCqQO27vF7TyfLv/HNu2yb6in/942Ly7uT5a//rZ+zbR786t+3dZ9+Kh3FW6z7Y+6v3P+8rfuDS/fZuv7qcrK82u0jfyczc7BxZ3o+hZKf66s+RVY91XQ/lsb7bZvZs/588UxvOqq8JxM9e+0lNK24lInF3pHe5ly0tNteSVrelX59KT/jB3DZxPBKkrrTkbC7MmH0V97utzd3rlseTo9h3b/0KDfuI0+lH3/OnYtzFnel15WL4c0prKb3ffe0H6NmyW9v15zfj11z6fcKzbI/Z4Waj6VduXPjg1he9f1bOpyenz2nfbT0ytH0VxdIfnuvNWwtOnej4oo/TkPFP5rv6tZq2RMQbiHXvQiJMT6l3Nnt2u9ESf9gszoFAAAAbEt8MB0AAAAANh8XIQAAAADaakPfEwIAAADgJtmCD4hvFe6EAAAAAGirLbsTEmNUczWdzFDMpCW0onhu0leW00OwtscnqJRmfUJSsz/d98qkT9eoD/p0mlwSV/VsevwmH/bjN/qcT6lw5o745ZX8UKhnKp0mU6j71KdsPw776RqeSScQrQ35PydcebdPKNn51XRSWE+3TxOpHfH70Rl+Or0eSZq53yckjX4zXTf+nf7vCidf9ilx2uvH4tXfvCdZHu+v2zYLj2bm2ZX0fLrzD/8X22Tf9160dXtK6X7s75uzbf7LM4/aur5Tfp84D/4VH5H05w/542dpOZ0mNDbq+z4x69N4ur6aPg7iUb+v6lVfVzApXfP9mWMxk+w1Ouq/5/byxGB6ebaF1PVkJnlsLF2+nJnrOd2n0vtxyoepqTLt61wCVs6up/35Z+5On0zVNLtreZffVztf8HV9E+m6Rp8/t+fSrEqL6X3Sk3npbnb5c51LesxZGvVtimP+3Fk9kZ7Ta31+f8Quf/x0j5vkpwX/HqJy3lYpXvSpiKqmU/Ma4/4bForVTDpoC4mnueSsws50Ol8h08a9r7zlcCcEAAAAADYfFyEAAAAA2ooPpgMAAACdgsexAAAAAGDzcScEAAAA6BBE9AIAAADATbBld0JCoaBifzrqba1mIuoklYYHk+XNK1O2TWH4gK1z8XXNA+loOElSn4/cbPSmY/4W9vrouuoZn3ObiwNumGjfgTfS0biStLzTxxCuDKSvSUtL/pLcRT9K0vR9G4/i7bvkYyFz18x9F9N9nLvDt+k66cdiyaSghqd8dPPeE37cL70zPRbTj/l41P7jvn+Tj6TL+87aJup9wuddzj0zautqD6XjELt6Wos6Hbw7fay6iFZJujDhx/277n09Wf7lb7zNthk46CNwD9w1a+uOvZo+l3zx1FHbpr7kD5LCbHofX5jx55gjD16wdeMHBpLlxXPpKE5J6r6ciSY1p5+FQ/44DXW/vQuv+Xm2w6wrFwOe02tSRueq/pwQe/wx7MYiZtJgF/b6uoGTmfNqOR1nO3PUv44U/KnE9mPoJT8WSyM+Urd6Nj1OXbP+nJB77Zl8pDdZPvKSj2LtOe1f85sP7LR1Li6/tOjHNhvLb94PlM/5/sXpWVsX9pgXn770GEmSMrG0wbxnkiSZSN3SIf+eSfWNn/cL8rG+jSk/TjLv626bGF7wOBYAAADQMXgcCwAAAAA2H3dCAAAAgE4QxZ0QAAAAALgZuAgBAAAA0FZb9jhWbDZtClZpJJNMZRSqPn0hxyVRLI512zbVE1dtXXMgnTZRXvRpMo1evxuaZZ9qM/VAOpWjet4nvKx1+cSTvovp1Iup+3z/Vn1okUa/md7m2n5/7Ttzl0/U6r3s70/OHUkvc3nUj/vqLp94Un6+b8N9mHhvJnWnnO5HaPj9UciEkBSX0+W5/bFwyh9X5aMLvh/1jf+tYi3TZvJcupMDYz4Vb+G1QVv35+MPJMt33D9j24z2z9u6l07st3X9J9PHwsL9/pjbMZRO45Gk4lfTY7H3Y2dsm1z/BkbT2xWn/f4Y+55ztu7K76RTckpX/fLWevwxsrzLH4/dl9PLXMnM6dUhv7yC2SWjT/tjbuWv+nmxMDOYLB887pd39U5bpYV9vt3iWHq7qqcyaX9X/bjvfSodndXsyqVj+XNxLunKL8+vyyWg5VIlG33+fFZa8MdjYTp9nilM2yaq59IyjTDs35Nkz6gu+amFVCpJijU/p4NJx1LFj3s2idS8D8smYGVs5xQsvicEAAAAAG4CLkIAAAAAtBXpWAAAAECn4HEsAAAAANh83AkBAAAAOgQfTAcAAACAm6Aj74TETCxbnPAxnk7uSmvtrn3pNg1/GVof7LF1lcl0HGezy0ftFuqtXfLueWo2WX7mLw3aNtUzfl0uknHkFR8N2Ojxo3vlwXTdwEnfh1xs5eoOX1c0EY89e3084e7/y8cwX3x3unzqkdb2VdOkXd5/r49HfaV/j60LE+m+F5f8GKmajumUpPceOm3rvvzcvcnyscNXbJsLE5lc1RbEMZNJLKn4Wvp4XHvK9+HskK8LVR/7unD/SrL8nXf5SN1vnDhk6/Tu9PLW5gZtExfDK0lXZ3qT5YVhv02nX0yfAyWpeU+6XbmWi6Nu7Rhx8b2h7ud0Pm53Nlk+XR+0bXb/5wFbVx5Ir6uZeRUNmVTVvgt+nPoupNe1mE6UlyQ1evxYnPu+dKTuzuf88prlzLnEPLDeNec3uLTk50xhNb3vcxHCPRf8cVCYTx9XktTMROc6pVkf5+40Mu8TCvM+AjdevJwsD1X/HkJ96eNeksKKfz/lonMLLrpX+a9DaJqvXShm2rivasj1YztH995uOvIiBAAAANiWeBwLAAAAADYfd0IAAACAThDFnRAAAAAAuBm4CAEAAADQVh35OFYuBcLWZRIgVMkkPSykUzS6J/31WS4pwyVi5JI81vp8/3JcSteuZ31CSS7xxKV0lRbWbJuVARP7JGnoeHp5F7/b968ynk5xuR6XULN6zs+lqfv9Ph5+JZ3WMn2fb1Ov+vunsT+9zecyKUguAUuSCiYxqLHDpyD1HavYuq+cfcCv6870fM8lYBVKvh9NpefM1QmfoNIzmk6dk6TVanqclkf9/sgd36Hh68pn02P4DfkErJhZnpvvV22L/H7UwfS4d037Piwd9Klprn8x88rRf3jO1uUSy7pn0uWrmTCjy9/h+37gl3Yky4u7/TlwaSSTzmf6MXTCnx+bZX9+XB7OpU+l7X7apz6Nf4efF0c+vZAsX9zrE5zK/pBT/+uzyfKr9/n923/WL3B5NH0M977h59LqqD+3VzLpWC45a224z7ZZGfXpU+51s/LCG7ZN7j3J6mP3pJscH7dtmmd8ymKo+HlRGhlJLy+TWJVLL7VtVvz+aIXrt+QTv24lYf1nO+BOCAAAAIC26sg7IQAAAMC2xAfTAQAAAGDzcRECAAAAoK14HAsAAADoEIHHsQAAAADgmhDCgRDCl0IIL4cQjoUQfny9fDiE8IUQwuvr/+/j8tZ15J2QxviErctFs1mZ+N5oIhSLC75Ns89H3rlIxlLdxzjO3Ovj/6pnNx6Hl4vhzcXtXnx3OjawuOSnSZdP8tPKcLr80O/kIjL95f/cEd9urSfd7s5PL9s2Z7/fj3uzlL4+L2e2d63b96/nbDrqdOmyP0Zjr4+5XSunt7fZ7dssPZSJSbzi5/Tacnr/F+d9/GhOeTE9tvWq7/vqqUxOq1vPLh+lve9t07auu+gjpN/4vw+nK5b9WAyM+UkTnk/v/4GTPsJz9qitstvcnPHxozu/uvFY7KVdvq7xNT+nF4/4sU2HyEp7v+j/VtYs+b5P35cur/hdr9JS7s+P6eO7a9Zv08qDfl5UTCTxtX6ky5d3+u3d81X/WjH5SDrOduCUb7M64M/7b/zV9Otwrg8uhrdV5UxUfi5u18XyF89N2jbdC+m4Z0kKNRM9PDxo2ygTc2ujeMt+f+RieHNfeWCXmYnozXH92OyI3tshhve6OvtOSEPSP4kxPhtCqEp6JoTwBUk/IumLMcZ/FUL4CUk/Iemf5xbEnRAAAAAA1xVjnIgxPrv+75qkVyTtk/QRSZ9a/7VPSfrL11tWR94JAQAAANB2O0MIT7/pvz8ZY/xk6hdDCHdIekTS1yTtjjF+61Gmi5J2X29FXIQAAAAAnWJrH8e6EmN87Hq/FELol/Sbkv5RjPFqCP/tkdUYYwzh+h+v53EsAAAAADckhFDWtQuQX4kx/tZ68aUQwth6/Ziky9dbDhchAAAAQCeI1yJ6t+rnesK1Wx6/IOmVGOPPvKnqs5I+vv7vj0v6zPWW1ZGPY5UOHfCV8yaJIqfik2acwrRPh6gf8Aldq4NuSH1Cxc5nZm+wV/+92bcNJMvLiz5lqNnlrzv7TChHPjHGcyldSyM+MaZ2yCdM5ZKp+sbT7S69M5c85pc384RJGZr36TTd4/5warZwpOWSrmIpvU8q475/pbO+E4v3+PSS0mT6+Om7Z9a2iU/6hCSX/LO0y8/N5b0+gciNRc8xn5Bz9n5bpWKmXZ9JNCod82M7W8qk0zxUTxZ3Xc7MMx/io8Wyme8mPU6SGj2ZtLrRdHku2Smn76Qfp4LZxZOP+OWV531dLgVro33ImTvik4kGTvpxz56nS36fOIW6X1fPVHpdS6N+nvVMpufmNelzQs9pPzEWjvrXTZfa2BjssW1yCZaNXr9dFZeClUmf0hnz4qjWnpoJmeSs5pV08lM2AavLv8eJNX+QuNSqZia9q1T1SYXRtMv1PZcUlusHttT7JP2wpBdDCM+tl/0LSf9K0m+EEH5U0hlJf/16C+rIixAAAAAAnSXG+JRcZrn0PRtZFhchAAAAQKfo7O8J2TR8JgQAAABAW3ERAgAAAKCteBwLAAAA6BA3klJ1O+BOCAAAAIC22rI7IaFQULHfRL3lYnhd3G6fj2LNRcC1onwuHaEnSYV6Opp07i4f+zmw6KMQVwd89F71TDrrdOI9fiz6x3NRnenyQiap0cXwStLQ8fS4X3jCb1MucrO07OtctKaLg5WkVZ80qNHPd5s2fnuvHvFj22sihBf3+jYjT2diX4+m261lolh9mIU09FUfobgynC6vdvtY3wtHfdZpoT89oYpnfRxnz6g/JyxNZo59Ixz3x2N46Kqtu3Iw3cfB530k6K6xWVt3eWIwXXFkwbZZafi+7ziZ3se5ubk4Zqt8P77q+7Dmd2PWwkPpg3XgqdYW2DS7pOtqa39idFHlPVPpeFlJWhnwceQrO/zfAMuL6XXVe/0xXHvYH8PuPDh03B/DlfOztm5XVzpuN1b8Oavngj+5F+bT/VjZP2jblC/5OOBKps7KvO/IReDKxO3Gi/572uL0rK0rjO1OV7Ty9QTX4aJzc3+dbkz59z9OITN+xUzk71otncufW95tE+vLnRAAAAAA2HxchAAAAABoKz6YDgAAAHQIPpgOAAAAADcBd0IAAACAThDFB9MBAAAA4GbYujshMSqupGP5XLkkKZ3Ylr+aysT3FqbTC2xOXLJtVp54u63rmkvHwxVXMxGZe30EZS7GsXo2vcxcDG95sWnr3Cjm+l7b7yMjx9+fjtErZqJ2q2f8ui69OxNn+3y670u7/LqKmfjemaPp7cpF4A6/5JdXO5Qu777sxy8XIezarT7mo12ViaWdu9fPi+JYeqCW/8se26bnB3zMbempHcnylXS6tSRpqebjR7sup09jD37/a7bN1189bOsqx9L9k6Q+M3cX3+/HvfHkqK0r7UrPp+JZH/m7lk6PvrY8EyNbqPt5lovFLpoo3ty+yukb93WDn3H72M/NuTv8+bH3shkLnx6txV1+nHw7H8Pbanxv/9l0HOvyqN/5Ayd9NGmhnh7DmXv9OUEatDVdsxuPQV3a15+pTdf1vjHnm6z4PtTvO2Dryqf8a7sTM7GvwcTtZmN93VcNSDaKN9uHzLpydU0TgZvTSjxuLja3lb+E3zYxvOBxLAAAAKBj8DgWAAAAAGw+7oQAAAAAHSCIiF4AAAAAuCm4CAEAAADQVtd9HCuE8IuSPiTpcozxgUT9E5I+I+n0etFvxRh/+nrLjTG2lHBQGhlJljevTNk2hXJmM03SQ3zkXttkYcwn1zgu7USSCgs+Daxwx4CtKy3WTY1Pr1jY7RNZRr+ZjsmZu8snqHRlwjUWDqUTWfZ+ubX7jNVTvu9Dr6bHt9GTSUZzwyefrDP5WCblbJe/pndJXPNHfSdKMxt/WnLwC35fzTzh48CajczfIybS6W3LmYSkRsPvq5E30olB5+/0yUTdp3w6lkuLeubpu2ybvssb31eS1DS7pHjMj/vAGz7daa0n3Y/lUd9m6Jjv+/yH0wdkyCR+LR70cVG7/zS9H5slv69y21ta8nVLI+l15dKsckmALtUvl6g18opf2cR70zt/6Lg/hmsH/bl44KQ/77sUrJ4LPsosdvnzhVve8As+fao+6FMbS7PpgySW/XHf++VXbF3zbem0ujDtU/Y0PGirSs+f9O3ca34ufarqk71iLb1Pwp5MNOOCfz/gEqsKVR+XmOt7Lm3UvQfLJWBtNpKujG3yONaNvMv5JUk/K+mXM7/zpzHGD21KjwAAAADc1q57ERJj/JMQwh1t6AsAAACwrYW4PW6FbNZnQt4TQng+hPD5EML9m7RMAAAAALehzYjofVbSoRjjfAjhByX9jqS7U78YQviEpE9IUrf8s/oAAAAAbl/f9p2QGOPVGOP8+r8/J6kcQthpfveTMcbHYoyPleU/aAoAAABsO3GLf9ro274ICSHsCSGE9X8/vr5MH1UFAAAAYFu7kYjeX5X0hKSdIYTzkv6lpLIkxRh/TtJHJf39EEJD0pKkj8X47X2iprR3zFeauN1CJq4vVvxmro6mo/dWB32b6lkfKdc1mY7rc+uRpNq9Pt6zb8LHPzZ601HBfeO+f81MZKSL4nVRl5JUO+SjOu/4TLrvC3t9H6pnfD7qwm7/CN/Ee9J1uQjPxV2+7/MH03UDr/k2DRMVK0lLD6W3a+eTPgazdtAvr2yikZcyqZClU35dq0OZSNiX0n+ryK1rRyYq+Pz/kI5BHTjm4z1Lmdjcwb95Lll+9sqwbbNS93Npx7TfxzMPp/s++hXf90uPZ6KHL5vyydb+PtT/2XSMZ+VqOhZZki4/4s91i2Yf56K5Z+/2fW+WfF31TPpYzR2nQ6/7SN3a/vR2jX3Fx6Muj/o79Lu/kR7DZpffpp3PzNq6nNJiuu+52NxGn5+DfcfTfxdcG/bHaeX8rK1zr6mFs5d8m0N7bZ1tY+JvJSnu9hnhrUTMxrsO2LqQ2S6neSZ9XpKkUPHzzNW56F5p82Nuic3detvlG9NvJB3rb16n/md1LcIXAAAAAK6Lb0wHAAAA0FabkY4FAAAAYDNsk8exuBMCAAAAoK24CAEAAADQVlv2OFYIwSZY5BIx5BI2zozbJs23HbZ1lVMmnuZOH/3jErBaNXDSR/80y5nEoEWfnOUUGj41xKVqzRz1SR4upUmSLr4n3a7HDLkkLY75iKlmOgxMkhRNXS5Zp54OEpIkDR5P3wudfsC3aQz5pJ6j/z6drPPGD/nlrR5csXV9X02PbS6xqmfS18VMalFpKT0WlUyKVM7I0+nTztRjfvxiKZPQ9mfpVJtmj2/TO+77PnevTwqrHk/3fWGfbaJef2qyqV8NPxTqqvntmjuS3q7GZX8eyc2LrqvpdeVS+yar/hwzYBKwJKm8mB73ninfZmnEb1f1fHoQZ+71yWh9F/3Arwyk1zX48pxtc/H9g7Zuz1Ozti6spvuRe8Euz/rXkfru9MmufMmfwHOpktG8LoV+P7ZhxY9tMMvLJVbp2AlfN7bbVjUn0klX4UQmESqTthVX/Hm6lTatJFPl0sBIuro1bZd0LO6EAAAAAGgrPpgOAAAAdAruhAAAAADA5uMiBAAAAEBb8TgWAAAA0AkiH0wHAAAAgJtiy+6ExBhtdFyx4iNh42unk+WFnSO2TXF6wXekbqIQM/G3zT7fv0uP9yfL9/5+OhZQklb2D9q6Qj0d7ZqzOuDj+nJxkqf/Srof+7/k4wTnjvix6LuYjtwsNPwlfqHu67ozkbBLo+m6PV/xfV/e6TN/Z+9OX5/nIon7z/rlLe7d+PX+6Jf8flw0UbzNzBHd8OnHKi77ukvfkZ6DLq62VeVpv7zeCd/ORS3Ho4u2zdpMn63rvuz3VcGkjHZf8PN2edjP26Hj6fk59YA/rnJWh9PHXLm28W2SpJ6pjZ9/Rl7xC2yW/Vis7Ej3MRdhXlrwEb3u+B442do5wZ07C/N+ebkY3vpgj62rnE+3ixW/P7IRuC28jrgYXkkqXJpJV5Qz54RMVGzhkj9WrarPWI/Ts7bORtZm+lca8e8vOiECtxP6gE3GnRAAAAAA2HxchAAAAABoKz6YDgAAAHSAID6YDgAAAAA3BXdCAAAAgE4Rt8etkI68CAnVdMKUJK3V0vFEhZXNTYcIqz5pZK3PpxaNfjOd8hErfqhzCVi1Q5kElbmNJ56c+UuDtu7gH6b73sylpGTSrEpL6aSeXEJO96SPaWr09Nq6gZPp8lzKUPe077tLwZo/6Nv0XPbbdemd6THsMiEz19M/nu7HlYd9m+pZ37+aD5rRwLH03M2lbS0c8cfP6FfcfPL9y6V+OcVjPgFr5wu+f3OH/cpcPxo9vu+Le/2cmVlKz0+XLCdJSyP+Bvb+P0qfEwqZ81nOwlg6LWroVX+czt3lx71nMpM6WPJjaNtkziX9Z9Pns1z/cudUlzBV3+0Pnkavn0tdc/41q9mfnhe5JK548bKt06G96fKFTCrVmXO+ziVTZV6Hm+a1W/LplrE2b9vEFT8WOUXTd/feQpIaU1MtrQtAHo9jAQAAAGirjrwTAgAAAGxHfDAdAAAAAG4C7oQAAAAAnSCKb0wHAAAAgJuBixAAAAAAbbVlj2OFEFToSkfdNsYnbLvSSDrKLxv/1797Y52Tj2OUpNLskq2LJs421HwUYjETgVs9Y6u0PJqOcSwt+L4Pv+KjP2fuTUfg7nxm1raZ+kg6wlOSui+np1fwKZ268qCPZ66e8fcnXURqLkZ2edjHey6Npst7x32bvkt+3Hum0u1WBvzfAVZ3+HWVzBQs+qmphk97ViUTFezGonrWtykc96eW1R3pcrdNkjTykj9+XIx19YxfoJvrklTPxBW76OZCJgF3x0m/H6vn0w2XRjKx2C2k7U4+5GPFB97w54ThF+aS5Uv7csepH/fFMX9Ads2lN8yd56T8ua6wkI5wHf79C7bN2l37bF1jMD3Pcq8H5Veu2LqF9xyxdT0X0tG0uaj3YF5PpcyTHZlI3YKL4ZWk/vTxE6dnfZuM3Gu+494/XE+r0b5AOwV/Wr6tcCcEAAAAQFtxEQIAAACgrUjHAgAAADoF6VgAAAAAsPm4EwIAAAB0iO3yjelbdhESY1RzNZ3MUcylcgwPJouDWZYkKZOisbI/vbxCJh1rddAvz6Wa1A+kU70kqXxuytYt3LHX1jVL6dSd7nGfFNY9bqu0vDc97rkknJydL6TTbnKJNjkX3+NTcoZfSS+z0eNThnIpSGs96TNAKZOOVdvv1zVwOj0WVx70NyNzKUjV8+ntLS355dUO+b67fSVJc4c39zRRyKSjOZfe6dOsRo6lj/2pB3yblWG/rl1P+3PJzNH0sd930Y/fwh4/foXVdARKs+zb9Ez546d7cjlZPqJMKtWs3966SYTqe9GnGdXe4c9Z3Vf8zm/0pY+f/tdnfRvTPymTVGheQyQprPr96M5bq6P+/Nh1xp9w3WuFJBUupePq4rCJlruOMJFO6cqlSoZ7Dtu6+NrpDffBvd63arOXB6D9eBwLAAAAQFvxOBYAAADQCaKkuD2ex+JOCAAAAIC24k4IAAAA0CG2ywfTuRMCAAAAoK24CAEAAADQVlv2OFYoFFTsz2SkGs0z55LlhVysb0bleDpCceWoj5l0sZqSFEy0b/mVdL8l2dhhSeqeXLF1y6PpyNrFOwZsm1wcZ7MrfU3aLPto19Fv+LpGz8aXt7LDXxdXz/j7k0sj6TjO7unW7mlWz6TLa4d8m4ZPQVWzlD7U9nzdR4LW9vvD021voe63t7hkq7Qy4OOFXRzwwm7fJhcv3Cyny5dGfZtd3/SxtGd/MD1OA6/65eXm0uqAH/eh4+njZ2nUbJSk0pJfV2kxPVCNHh8DnuMiZnvfmLNtcucLF/mbixzPReo2+3zMdnk2PUHDip9M5Us+YjZW0vvRxbJfT9dkOlLXvYZIkqo+vjdM+77LRDSvvXTcNilkouhdnG2uTSfE8ALbFo9jAQAAAMDm44PpAAAAQAcI4oPpAAAAAHBTcBECAAAAoK14HAsAAADoBDFum29M37KLkNhsaq2WTgcpZpKuCjvTqSyxlk4ukSRdvGyrgkmmqpyf9cvLWVhMl2cSsBaO+qSZngt+uzJhTNb4d/bauoOfnUqWzzw0bNssD/ukq6Hj9WS5S+GSpJ4pn4LkEqFyqmd9WsvCXp8M0+hJb9eer/i0srkjPvlnaVe6PJeAlUtVWtiX7l/fBdskm1iVU+9N76/5g77Nvif9uLskqep5nzpXWvDzYuDV9H4cOt7aviqu+nFv9KXn4MCJBdtm/qA/5lYH030vZZLM3P6QfDrf6qhPaep71if3xeEdyfJGZnmNwR5bV1zw86Ju2lUuzdg26vNjG2rpc3G5nEl1m/dzRtOz6fJ+34fo2kgKuTQr89rYKpeCRZoVgK3EnRAAAACgQ/DBdAAAAAC4CbgIAQAAANBWPI4FAAAAdAoexwIAAACAzcedEAAAAKBDbJcPpm/ZRUgIwcYG5jTGJzbcpnTogK+sm9xSV34dSw/sS5Z3zfkoxK5ZX5eLjJx920CyfOj5adtmz1d8PGXMRFc63dP+SFkdSE+v0pKPYm30+JtzLjZXklZNqnPPoJ/iq9VcvHB63Kce8NGuzczRFNJpxdnY3PKiH9vB19N1C3v8+NV98rW6rvq6Zjk9Tr2ZQ9HF8Eo+1nlxl59/A2/4dbkxdHG6Un5f5eaniwqeubfPtsnFThdW0+va9eVLts3cw6O2LqykB6O0aCagpPqdu21daTadFVw25ZIUu/zguhje7DLLfnmx4uuCiUtvKYY3Z97Esks2hl7KR9GHijnPZCJ1idsFcKvhcSwAAAAAbcXjWAAAAEAniJKa2+N5LO6EAAAAAGgrLkIAAAAAtBWPYwEAAACdYns8jbV1FyExRpvmkbs94xK1CmM+4SVmEk/C8GC6IpOOtXbAp9NUJtNJKblUmFa5FKz5uwdtm0I9N7PTfRw4sWBbLO7125Vf18b1TPnUoq5aOnFpYbdPSOq75FOL5o6k02lGv+mTcGqHMsk/i+m+T9+XS3DyR4Ibi5GXfPJPoZ7bXt/3gZPpZboxkqTiqt/31fPpvrvkKek6aVtD6fKdz8zbNn0v+tSi+oERWzfxnt5kef+4397cdrl90uz3Y9t9JZN0tTsTgWaUT/kkrji8I12eScBqlv28rZyf9R1pId0pTPtYN5ekWNo7Zts0M2lWraRP5RKwcslZALAdcCcEAAAA6BDb5XtC+EwIAAAAgLbiIgQAAABAW/E4FgAAANAp4vZ4Hos7IQAAAADaijshAAAAQIfYLh9Mv+UuQgou8jATqRuq/bYuVtORm41MpG7p+ZMbXldzNL2e62klcjMXCeoihCVprS8df7w82m3b5GJ4l0bS8bN9Ez5itHfWx2Aujvl+uEjYXV+bs21WMvukasZw8hHfpuA3S6tVMxYXMlG2ZzceCdro85G/C3t87GtlLhN/PJmOuq2W07HIkjT5UHouSVL3TLp8YNYfw8NPnbd1je87kCx381mSrjyabiPlx/3A56eS5Yt3DNg25dklW1eYT8cf5yJ6S4uZiWYUL5pBl7TwDj8WvW+kj5/CQiYKOtOP5plzvp05t8dMNG5cyfTDxLk3r6T34c1ADC8AeNd9HCuE8IshhMshhJdMfQgh/J8hhBMhhBdCCO/Y/G4CAAAAuF3cyGdCfknSBzP1PyDp7vWfT0j6j99+twAAAIBtJm7xTxtd9yIkxvgnktJfzX3NRyT9crzmq5IGQwj+K2kBAAAAbGub8ZmQfZLe/KDv+fWyiU1YNgAAALAtBElhm0T0tvWD6SGET+jaI1vqVmsf1AYAAABwa9uMi5ALkt4cr7J/vewviDF+UtInJWmgMBJtekkmDUUmbaQw7BNeckItnRbV8sCU0y1zqVSzb/PJOn3j/mq4cn42Wd44PGTb5BKDnEaPf2qveuKqretJzgKpdtcO2yY3TtUTfl647cqNbfWMTy1q9Kb3Y/f05v51omfSJx1NPeATkrqupvtR8AFT6pnyqWldmWSqsJKua3b5eeESsCSp72J6eS6FS5KWju62dTufmU2W59LPcmOR2y6XmteVSXVzCViST8EqTGdSlcx5M7e85rBP2et71idWObHm91WOTTeUT8HKJWDlZF9HAABbbjO+rPCzkv72ekrWuyXNxRh5FAsAAADYqOYW/rTRdf/gH0L4VUlPSNoZQjgv6V9KKktSjPHnJH1O0g9KOiFpUdLfuVmdBQAAAHDru+5FSIzxb16nPkr6B5vWIwAAAGCb2i4fTN+Mx7EAAAAA4IZxEQIAAACgrdoa0QsAAADA2IJvLt8qW3YREmO0EYqlkRHfzsUuTs/6lVV8pKWLmgyZGEzt2eXrjFD3kaB94z5KstFX9HVH0+O0NOLbDJz0UazNcrrdwB+/btusvP0OW1eeTUfg5mJ9rzw6aOuGXvXxvbVD6ejU8qKPepi5d+MRrsVVf2YYeG7S1sVK+lBb3puJTr3o+959JR3tuzroD+lC3fd97oiPA945m67rnly2bUoLZVvnrI7227qe0z7zd224L1neNeePq+5xH/vqYrslSXVz/PT5ueRicyVJx04kixe+/2HbpO9FHz4YTlz263Jykb8mEj0Xf+ui1yVpzSwv146oXQC4PXEnBAAAAOgIUeKD6QAAAACw+bgIAQAAANBWXIQAAAAAHSLErfu5bt9C+MUQwuUQwktvKhsOIXwhhPD6+v8P3ch2chECAAAA4Eb8kqQPvqXsJyR9McZ4t6Qvrv/3dXXkB9MbU1O2ziWohOHBltYVqiaRJ5PIEqs+CSespNNzokmekqTK+Vlb5xKwJJ9O1O1DmrS4N50ida1dOjGoft8B26a0mE5pkqTCfHp5ubSgnc/M2rrY5afr4MtzyfJc2pZLwMpxqVSSVN/tk65cElcrfZDyKVhO3/HMcbXq/2hRmE4nGk2/f79t0zPpx8m2OX7JV7pUKknL5hjJpUjlEqFymrvT4xROnLNtCu4cI0nV9Jzp+cPnfJud/pzgzmcuBVDKJA6qtWSqVtOsSMECgM4XY/yTEMIdbyn+iKQn1v/9KUlPSvrn11tWR16EAAAAANvS1qZj7QwhPP2m//5kjPGT12mzO8b4rb/8XZS0+0ZWxEUIAAAAAEm6EmN8rNXGMcYYwo18uoSLEAAAAKAzRCn47yruVJdCCGMxxokQwpikG/rmXD6YDgAAAKBVn5X08fV/f1zSZ26kERchAAAAAK4rhPCrkr4i6WgI4XwI4Ucl/StJ3xdCeF3S967/93XxOBYAAADQKbb2g+lZMca/aaq+Z6PL6siLkKKJrcyJ07O2zsbwykduFi7N+OVNX/Ud6UtHsbq4WkmKFb8bcrGqNkY4Ez/a1Vu2daXZJb8uo9nn43a1sJgsXstE2TYzUca5OGAnF/kb6j4eN9TSfc/FMzcGffzx0Kvp5TV6/b7vmvORpcWL6fm5dNQHUiwd9jG83ePpGF5JWtuTbldc9SfJ3PLc2DaH/bzIReA2yyFZXj+QibJ96jlbVxrx7VxEuPbssm2ycd8mOreQieFtjGeihwEAuEV05EUIAAAAsC117o2QTcVnQgAAAAC0FRchAAAAANqKx7EAAACADhE6+IPpm4k7IQAAAADaasvuhIQQfNJMK8trcVkuBSsO77BtcolQhYV0ClbMpD4Vpn2SUC7hp7iQTt1ZGfUJTpXJdDKR5FOLNO/bNO87YOvm33coWd59xadcVU75L9nMjUVpIZ10ld1XZy/ZurW79qXbmP0r+f0hSbErfajlErBy63KJVdlUqpWGrcspmnnRW/Z/w8iuq56uc2MkSYXhQVtXfXZ8Q+uRJB3y8zbbLpPC58RcOtZKeh+v1fx+BADc5rgTAgAAAACbj4sQAAAAAG3FB9MBAACAThAlNbe6E+3BnRAAAAAAbcWdEAAAAKADBEUiegEAAADgZrjl7oS4SMtsRG9545uZixgt5uJHF9JxpqHPx+bmIkHLr5zz7Srpbe5pIXpWkmI13cdW4493vJKOP86Zfv9+Wzf0/PSGlxdW/djmxqLRW06WFzKxtMujmTjgevqvGt2Ty7ZNLtbZ9a9ybtK2qd+529aVT/k5IzN3s21yx9xKOrK28PJp36ba7+tyx5bRnMj0fZM1MxG9AABsV7fcRQgAAABw2+JxLAAAAADYfFyEAAAAAGgrHscCAAAAOgWPYwEAAADA5uvIOyGtpjFZ8+nEKkmKJrnGpXBJUmHniF+XSwWanvVt+n26TzbR6FJtY32QVJxe8P1wKT6Z/VGaXfLLM5b3Vm3d8O+/Zuvi2E5b59LMciln9d2Z7VqsJ8uLJy7YNoW6T9tyKV1rfb4P5RM+walyybTL7Hs7XySbWJVdZittJDvfQz1z3GfW1bwylV5exaeVkVgFAOhIfGM6AAAAANwcXIQAAAAAaKuOfBwLAAAA2I4CH0wHAAAAgM3HnRAAAACgU3AnBAAAAAA235bdCYkx+pjMWiZK1Gg11jdU+9Ple3b5Rgs+8tfG2aZXI0laOupjeCuTmXW1EjOaiwq2nfAxt0E+XjhW0tMrt025GN5YLtq6xmBPsrxQ9zl35VfO2bq1u0zcbiZOuXhu0ta5iNlCZnktycZRz/p2h/baqjB99dvo0F/UOOPHfVMRwwsAQMficSwAAACgI0QexwIAAACAm4E7IQAAAEAniOJOCAAAAADcDFyEAAAAAGgrHscCAAAAOoUP9rytbNlFSAhBBRNnGyoV2y6urKTLW43jtDHB87ZJNg64Lx25Gqs+irXn9IxfXguRujZeVlKzPGLrypfS0cjN/sz+6PJTqHgxvV258cutq3D2ku+H2eaw6uOFs3G70wvpihajoJtXppLlmx7Rm5Eb98ZLx21dsVpNlq/lorSn0tsryR73NrIbAADcdrgTAgAAAHSIwAfTAQAAAGDzcRECAAAAoK14HAsAAADoFDyOBQAAAACbb8vuhMQYfRpOJiXHJmpV+1vqh0stCvcc9m1eO23rgum7S/WSJGXSwHLb1dw9lCwvnpu0bYomvSunlQQsySeChYkrtk1Bg355Yztb6odVzyRnOfOLvi6TdFU4dCBdkZnrzUz6VEHpxKpcSlxuDroELOk6KVgtIAULAAAjSmpyJwQAAAAANh0XIQAAAADaig+mAwAAAB0h8sF0AAAAALgZuAgBAAAA0FY8jgUAAAB0im3yONaWXYSEEGzcbktWfOxnYyodwyv5yN+Yi+FtMVK3pTaZSN1QX9vwujQ9a6tcJGzhot/eeGivrQvTV5Pla3fts21y8cJrr0/YupKLwG2RjRducXmN108my3PRuNnlZeZ0S4jNBQAAbcSdEAAAAKBTbJM7ITf0mZAQwgdDCMdDCCdCCD+RqP+REMJkCOG59Z+/u/ldBQAAAHA7uO6dkBBCUdJ/kPR9ks5L+kYI4bMxxpff8qu/HmP8sZvQRwAAAAC3kRt5HOtxSSdijKckKYTwa5I+IumtFyEAAAAAWhUlNXkc61v2STr3pv8+v172Vn8thPBCCOHTIYTN/ZQwAAAAgNvGZn0w/Xcl/WqMcSWE8PckfUrSB976SyGET0j6hCR1y6c+5RRcmlDFJ23lUrgKY7uT5c2JS7ZNyCyveWXjqUVheNAv78w5W1fYObLhdbkErOzyMsljMglYkmyyV/HEBd8mI5skNb+YLM5tb06hnD40YiZdrBVxZcXWNUmsAgBgm4lSbG51J9riRu6EXJD05jsb+9fL/qsY41SM8Vvvpn5e0qOpBcUYPxljfCzG+FhX6G6lvwAAAABucTdyEfINSXeHEA6HELokfUzSZ9/8CyGEsTf954clvbJ5XQQAAABwO7nu41gxxkYI4cck/YGkoqRfjDEeCyH8tKSnY4yflfQPQwgfltSQNC3pR25inwEAAIDb0zb5npAb+kxIjPFzkj73lrKfetO/f1LST25u1wAAAADcjvjGdAAAAKATENELAAAAADfHlt0JiTHaCNJcFGtjKh2Bm2tjY33lo3gLh/xXnWRjc03kbzba1cTLSvm+x9p8snwtE0tbym1XJpbYCasVWxdNXHGo+Da5uOI44bfLbXMunjnXDzcWmx2bSwwvAADYjngcCwAAAOgU2+SD6TyOBQAAAKCtuBMCAAAAdAruhAAAAADA5uMiBAAAAEBbdeTjWLl0J5d2FDIpSM3c8lxaVC7NKqNhkrNyKU1xZcXW5dKx3DiVRkZsm1wSVy4tqhWFnel+uFQv6TopYi3Ipk+RTAUAADpK5HEsAAAAALgZOvJOCAAAALDtREnN5lb3oi24EwIAAACgrbgIAQAAANBWPI4FAAAAdAo+mA4AAAAAm++2uRPSmJqydcVMzG3TROrm5GJfbYRwJv42VPv9uq747crF/trlZeKKN5uLEM71u5WxBQAAwK3ltrkIAQAAAG55PI4FAAAAAJuPOyEAAABAR4hSkzshAAAAALDpuAgBAAAA0FZb9jhWKBRU7E+nVrlUJcmnTOWupsLwoK+cnk0Wx5UV2ySX0lQwSVy59C5ltjeX7OX6GDMJU63IJVa1c3mb3Q8AAICOEqUYm1vdi7bgTggAAACAtuKD6QAAAECn4IPpAAAAALD5uAgBAAAA0FY8jgUAAAB0Cr4xHQAAAAA235bdCYnNZjaK17bLROc6zYlLts5F/rYcI2uieHOxvq2yfSTKFgAA4NYTo9QkohcAAAAANh0XIQAAAADaig+mAwAAAJ2CD6YDAAAAwObjTggAAADQIeI2+WB6R16E5JKkWkmtyi3PJXRtdh9ybXLraiUNDAAAAOhkPI4FAAAAoK068k4IAAAAsP1EPpgOAAAAADcDFyEAAAAA2orHsQAAAIBOECU1eRwLAAAAADZdR94JaSUCt1it2joXw7vZfZDycbubvS4AAADcZuL2+J4Q7oQAAAAAaCsuQgAAAAC0VUc+jgUAAABsN1FS5IPpAAAAALD5uBMCAAAAdIIYt80H02+bi5C4stJSO5dm1WpiFUlXAAAAQB6PYwEAAABoq9vmTggAAABwq+OD6QAAAABwE3ARAgAAAHSK2Ny6nxsQQvhgCOF4COFECOEnWt1MLkIAAAAAXFcIoSjpP0j6AUlvk/Q3Qwhva2VZXIQAAAAAuBGPSzoRYzwVY1yV9GuSPtLKgkKMW/PhlxDCpKQzW7Jy3Kidkq5sdSdwQ9hXtw721a2DfXXrYF/dOrZqXx2KMY5uwXo3JITw+7o2RlulW9Lym/77kzHGT37rP0IIH5X0wRjj313/7x+W9K4Y449tdEVblo51K0yE7S6E8HSM8bGt7geuj31162Bf3TrYV7cO9tWtg32VF2P84Fb3oV14HAsAAADAjbgg6cCb/nv/etmGcRECAAAA4EZ8Q9LdIYTDIYQuSR+T9NlWFsSXFSLnk9f/FXQI9tWtg31162Bf3TrYV7cO9tUtLMbYCCH8mKQ/kFSU9IsxxmOtLGvLPpgOAAAAYHvicSwAAAAAbcVFCAAAAIC24iIECiF8MIRwPIRwIoTwE4n6HwkhTIYQnlv/+btb0c/tLoTwiyGEyyGEl0x9CCH8n+v78YUQwjva3UdccwP76okQwtybjqmfancfcU0I4UAI4UshhJdDCMdCCD+e+B2OrQ5wg/uKY6sDhBC6QwhfDyE8v76v/tfE71RCCL++flx9LYRwxxZ0FVuID6ZvcyGEoqT/IOn7JJ2X9I0QwmdjjC+/5Vd/vZUvosGm+iVJPyvpl039D0i6e/3nXZL+4/r/o/1+Sfl9JUl/GmP8UHu6g4yGpH8SY3w2hFCV9EwI4QtvOQdybHWGG9lXEsdWJ1iR9IEY43wIoSzpqRDC52OMX33T7/yopJkY410hhI9J+teS/sZWdBZbgzsheFzSiRjjqRjjqqRfk/SRLe4TEmKMfyJpOvMrH5H0y/Gar0oaDCGMtad3eLMb2FfoEDHGiRjjs+v/rkl6RdK+t/wax1YHuMF9hQ6wfqzMr/9nef3nrUlIH5H0qfV/f1rS94QQQpu6iA7ARQj2STr3pv8+r/RJ/a+tP4bw6RDCgUQ9tt6N7kt0hvesP6rw+RDC/VvdGUjrj4M8Iulrb6ni2OowmX0lcWx1hBBCMYTwnKTLkr4QY7THVYyxIWlO0khbO4ktxUUIbsTvSrojxvh2SV/Qf/vLBYDWPCvpUIzxIUn/XtLvbG13EELol/Sbkv5RjPHqVvcH3nX2FcdWh4gxrsUYH9a1b9R+PITwwBZ3CR2GixBckPTmOxv718v+qxjjVIxxZf0/f17So23qGzbmuvsSnSHGePVbjyrEGD8nqRxC2LnF3dq21p9Z/01JvxJj/K3Er3BsdYjr7SuOrc4TY5yV9CVJH3xL1X89rkIIJUkDkqba2jlsKS5C8A1Jd4cQDocQuiR9TNJn3/wLb3n2+cO69hwuOs9nJf3t9SSfd0uaizFObHWn8BeFEPZ869nnEMLjunYu5sV3C6zvh1+Q9EqM8WfMr3FsdYAb2VccW50hhDAaQhhc/3eProXfvPqWX/uspI+v//ujkv448g3a2wrpWNtcjLERQvgxSX8gqSjpF2OMx0IIPy3p6RjjZyX9wxDCh3UtmWRa0o9sWYe3sRDCr0p6QtLOEMJ5Sf9S1z7spxjjz0n6nKQflHRC0qKkv7M1PcUN7KuPSvr7IYSGpCVJH+PFd8u8T9IPS3px/fl1SfoXkg5KHFsd5kb2FcdWZxiT9Kn1BM6CpN+IMf7eW95b/IKk/xRCOKFr7y0+tnXdxVYIHJsAAAAA2onHsQAAAAC0FRchAAAAANqKixAAAAAAbcVFCAAAAIC24iIEAAAAQFtxEQIAAACgrbgIAQAAANBW/3/9WtEy5aIt9AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x864 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# try a 2d histogram for radiance\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "size = 100000\n",
    "x = np.random.random(size)\n",
    "x = np.arccos(x*2-1)\n",
    "y = np.random.random(size)\n",
    "y = np.arccos(y*2-1)\n",
    "y = y + x/2\n",
    "print(np.min(x))\n",
    "print(np.min(y))\n",
    "\n",
    "print(np.max(x))\n",
    "print(np.max(y))\n",
    "\n",
    "\n",
    "fig=plt.figure(figsize=[15,12])\n",
    "(h, xedges, yedges, image) = plt.hist2d(x,y,bins=100)\n",
    "print(h)\n",
    "plt.colorbar()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "fe825f3a",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/joel/.local/lib/python3.8/site-packages/cupyx/jit/_interface.py:161: FutureWarning: cupyx.jit.rawkernel is experimental. The interface can change in the future.\n",
      "  cupy._util.experimental('cupyx.jit.rawkernel')\n",
      "/home/joel/.local/lib/python3.8/site-packages/cupyx/jit/_interface.py:161: FutureWarning: cupyx.jit.rawkernel is experimental. The interface can change in the future.\n",
      "  cupy._util.experimental('cupyx.jit.rawkernel')\n"
     ]
    }
   ],
   "source": [
    "# try a 4d histogram\n",
    "\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import optics_cuda\n",
    "\n",
    "source = optics_cuda.MonochromaticLambertianSource(0.001, 0.001, 555, 5e7, 0.001)\n",
    "photons = source.make_photons(10000000)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "a685b04b",
   "metadata": {},
   "outputs": [],
   "source": [
    "import cupy as cp\n",
    "import scipy.constants\n",
    "bins = 10\n",
    "\n",
    "points = (photons.r_x, photons.r_y, cp.arccos(photons.ez_z), cp.arctan2(photons.ez_y, photons.ez_x)) \n",
    "wavelength_m = photons.wavelength_nm * 1e-9\n",
    "frequency_hz = scipy.constants.c / wavelength_m\n",
    "energy_per_photon_j = scipy.constants.h * frequency_hz\n",
    "energy_per_bundle_j = energy_per_photon_j * photons.photons_per_bundle\n",
    "\n",
    "(h_j, edges) = cp.histogramdd(points, bins=(bins,bins,bins,bins), \n",
    "                              range = ((-0.0005,0.0005),(-0.0005,0.0005),(0,np.pi/2),(-np.pi,np.pi)),\n",
    "                                       weights=energy_per_bundle_j)\n",
    "\n",
    "# hack to see the total\n",
    "waves = 40\n",
    "h_j *= waves"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "dad5fc73",
   "metadata": {},
   "outputs": [],
   "source": [
    "# bins are 0.00005 m squares,  which should yield 2.5e-9 m^2 areas\n",
    "\n",
    "bin_area_m2 = cp.outer(edges[0][1:] - edges[0][:-1] , edges[1][1:] - edges[1][:-1])\n",
    "bin_area_m2_stretched = bin_area_m2[:,:,None,None]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "92a9cf31",
   "metadata": {},
   "outputs": [],
   "source": [
    "bin_area_sr = cp.outer(cp.cos(edges[2][:-1]) - cp.cos(edges[2][1:]) , edges[3][1:] - edges[3][:-1])\n",
    "bin_area_sr_stretched = bin_area_sr[None,None,:,:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "f1afd6be",
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "#bin_area_sr_s_m2 = (bin_area_m2 * bin_area_sr)\n",
    "#print(bin_area_m2) # 10x10 1e-8\n",
    "#print(bin_area_m2.shape)\n",
    "#print(projected_bin_area_m2) # 10x10 1e-8\n",
    "#print(projected_bin_area_m2.shape)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "3bf4635a",
   "metadata": {},
   "outputs": [],
   "source": [
    "#bin_area_sr_s_m2_stretched = bin_area_sr_s_m2[None,None,:,:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "a232ad2e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0.07853982 0.23561945 0.39269908 0.54977871 0.70685835 0.86393798\n",
      " 1.02101761 1.17809725 1.33517688 1.49225651]\n",
      "[[[[0.07853982]\n",
      "   [0.23561945]\n",
      "   [0.39269908]\n",
      "   [0.54977871]\n",
      "   [0.70685835]\n",
      "   [0.86393798]\n",
      "   [1.02101761]\n",
      "   [1.17809725]\n",
      "   [1.33517688]\n",
      "   [1.49225651]]]]\n",
      "(1, 1, 10, 1)\n",
      "[[[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]\n",
      "\n",
      "\n",
      " [[[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]\n",
      "\n",
      "  [[9.96917334e-09]\n",
      "   [9.72369920e-09]\n",
      "   [9.23879533e-09]\n",
      "   [8.52640164e-09]\n",
      "   [7.60405966e-09]\n",
      "   [6.49448048e-09]\n",
      "   [5.22498565e-09]\n",
      "   [3.82683432e-09]\n",
      "   [2.33445364e-09]\n",
      "   [7.84590957e-10]]]]\n"
     ]
    }
   ],
   "source": [
    "theta_bin_centers = (edges[2][:-1] + edges[2][1:]) / 2\n",
    "print(theta_bin_centers)\n",
    "theta_bin_centers_stretched = theta_bin_centers[None,None,:,None]\n",
    "print(theta_bin_centers_stretched)\n",
    "print(theta_bin_centers_stretched.shape)\n",
    "\n",
    "print (bin_area_m2_stretched * np.abs(np.cos(theta_bin_centers_stretched)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "302df5f2",
   "metadata": {},
   "outputs": [],
   "source": [
    "h_w = h_j/source._duration_s\n",
    "\n",
    "intensity_w_sr = h_w / bin_area_sr_stretched\n",
    "#radiance_w_sr_m2 = intensity_w_sr / projected_bin_area_m2_stretched\n",
    "\n",
    "radiance_w_sr_m2 = intensity_w_sr / (bin_area_m2_stretched * np.abs(np.cos(theta_bin_centers_stretched)))\n",
    "\n",
    "#radiance_per_bundle_per_angle_w_sr_m2  = h_j/bin_area_sr_s_m2_stretched"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "ae970de3",
   "metadata": {},
   "outputs": [],
   "source": [
    "max_radiance_w_sr_m2 = cp.amax(radiance_w_sr_m2, axis=(2,3))\n",
    "#print(max_radiance_per_bundle_per_angle_w_sr_m2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "4cbefa9b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0IAAALNCAYAAAACz9HNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAA8QElEQVR4nO3de7xuV1kf+t/DzoVLAgTCNRcSJKDgBTREekAEEYhUD9qDNejBaFFaFQsttdw+FYv1c6S0aC0iRqAgokElQloDISIIqISENBCSNLJJgCRGYhIgAbntvZ7zx5qBdy3XZSfZe82VNb7fz2d+9vuOOeaYY77rzc569vPMMau7AwAAMJI7zD0BAACArSYQAgAAhiMQAgAAhiMQAgAAhiMQAgAAhiMQAgAAhiMQAgAAhiMQAoADoKpOqqq/rqr3VtUfVNXBc88JgK8TCAHAgXFlku/p7scm+USSp847HQAWCYQAbqOquriqHjf3PA6Eqnp9Vf2n6fWOvc4Dobuv6e4vTm+/kmTp5n1V9cGqetg8MwMgEQgBA6iqT1TVV6rqyFXt/7uquqqOuy3jd/fDuvs9t2WM24PtdJ1V9cKqevuqto+t03bKwvv7V9VVWzXP6ZwPSPKkJP9zofm/JHnpVs4DgJUEQsAorkjy9JvfVNW3JLnzfNPZWlV10Nxz2M/em+T/qqpdSVJV90tycJJHrGp70NT3Zk9J8o5bcqLb8tlV1V2TvDHJT3T3Vxd2nZnk8VV131s7NgC3jUAIGMUbk/z4wvtTk/zuYoeqekFVfbyqbqqqS6rqh6b2b6iqG6rq26f396+qv7+5TGzKOH3vwjifqKpfqKqPVNUXquq1VXWfqnr7NPafVdURU9+uqgctHPu1UrRbOtZq07HPr6qPJPlCVR203jUuHPOIqrpg2v/mJHdcNd7ida471tT3303z/lxVvbmqFsc6pqrOmD7H66vqlQuf7Vum9iuq6l+v8/M8L8uBz8On99+V5N1JLlvV9vHu/tuF456S5KzpXM+vqqun+V9WVU/Y6LO7pT+P6ZjTk/zH7r5scfLd/aUkH0ry5HWuD4ADTCAEjOIDSe5aVd80ZQxOSfJ7q/p8PMu/PN8tyX9M8ntVdb/u/niS50/v75zkfyR5wyZlYv9PkicmeXCSH0jy9iQvSnKvLP/du94v+Pt7rKcn+adJ7t7de9a7xiSpqkOSvDXLQeM9kvzRdO71rDvW5J8nOTnJ8Um+NclPTOfZleR/JflkkuOSHJXk9Kq6Q5bLxz48tT0hyXOr6h8FC939lSTnJnns1PTYJO9L8v5VbV/LBtXyqm2PTXJOVT0kybOTPLK7D89yQPKJTT67m+3rz+PpSb4zyX+oqvdU1Y+sGv/SJN+2+toA2BoCIWAkN2eFnpjlX0KvXtzZ3X/U3X/b3Uvd/eYkH0ty0rTvd5LszvIv3/dL8uJNzvXfu/vT3X11ln9BP7e7//eUCfiTJI+4BfO+LWP9RndfefNN+xtdY5JHZTnL8uvd/dXu/uMsZ17WtMlYN5/7b7v7hiwHOA+f2k9Kcv8kv9DdX+juL3X3+5M8Msm9uvul3f2V7r48ye9kOWhdy1/k60HPd02fzftWtf3FQv/HJvlwd9+UZG+SQ5M8tKoO7u5PTAHvup/dgn36eXT3G7v7nt39uGl786pxbkpy93WuDYADTCAEjOSNSX40y5mJ3129s6p+vKourKrPVtVnk3xzksUFFn5navvv3f3lTc716YXXX1zj/WG3YN63ZawrF99sco33T3J1d/fCIZ9cb+B9+Lz+buH1PyzM85gkn1yVZUmSByS5/83jTWO+KMl91pnCe5M8pqrukeUA6mNJ/irL9w7dY5rP6vuDzkqS7t6d5LlJfinJtVV1elXdf9X4V2Zt++tne3iSz+5jX4DZVNXrquraqvroPvb/51PJ9MVV9fsHen63lkAIGEZ3fzLLiyY8JckZi/tqeWWv38lyudQ9u/vuST6apKb9hyX59SSvTfJL0y/a+8M/ZOWiDfv75vmvBTWbXWOSa5IcVVW1cPyxaw26D2Nt5Mokx9Y/XoTgyiRXdPfdF7bDu/sp64zz11kuy/vpJH+ZJN19Y5K/ndr+truvWOj/tUBo6vv73f2YLAdgneRlq8bvHFjflOUyQIDt7vVZLnXeVFWdkOSFSR7d3Q/L8j86bUsCIWA0z8zyQy6/sKr9Lln+xffvk6SqfjLLGYWb/bck53f3TyX50ySv3k/zuTDJj1bVrqo6Ocl376dx17LZNf51kj1J/nVVHVxV/ywrS91uyVgb+WCWg65fraq7VNUdq+rRU/tN0yIFd5o+k2+uqkeuNchUsnZ+kn+b5RK1m71/alu8P+j4JId296XT+4dU1fdU1aFJvpTlTM5Stsi0cMR3JDlnq84JcGt193uT3LDYVssLCb2jqj5UVe+rqm+cdv10kt/s7s9Mx167xdPdZwIhYCjd/fHuPn+N9kuS/NcsBwOfTvItmbIMVfXULP9L2M9M3f9tkm+vqh/bD1N6TpZvuP9skh/L8mIFB8RG1zjt/0qSf5bl0sEbkvxIVmXO9nWsTeaxN8vX/KAkn0pyVZIfmdq/P8v3El2R5Lokr8ly1mc9f5Hk3lkOfm72vqltsSzun2YhG5Tl+4N+dTrH3039X7gv899PfiDJe1ataAdwe3Jakp/v7u9I8u+SvGpqf3CSB1fVX1bVB6Z/5NuWamUpOADsPFV1VpJXdvdZm3beAlV1bpJndvc+1dsDzK2WHz7+v7r7m6dy8b/P8iMLbnZod39TVf2vJF/N8sqhR2f5H6W+pbs/u8VT3tROe8AeAKzlPVl+ztC20N3fOfccAG6DOyT5bHc/fI19V2V5Nc2vJrmiqv4myQnZYBXSuSiNA2DH6+7/vMYy2ADcCtPCNFdU1Q8nSS27+blob03yuKn9yCyXyl0+wzQ3JRACAADWVVV/kOV7Qh9SVVdV1TOzfF/rM6vqw0kuTvLUqfvZSa6vqkuynIn/he6+fo55b8Y9QgAAwHBkhAAAgOFYLGEDuw6/Sx90r7vPPY1ZHXL9vjwbcWerJVnTJNlzZ/9uUnvmnsH87vBV/z0kSR/k78aoKCFJ7Z17BtvD52+6+rruvtfc89jIkx9/l77+hnl/YB/6yJfP7u5ts5y2QGgDB93r7rn/f/q5uacxq2N/z1fkoC/67TdJ/v5b7zT3FGZ3xxu27Hmb29Zdrvnq3FPYFr505MFzT2F2u77ivweSg28UCSXJe971wk/OPYfNXH/D3nzw7GNnncOu+33syFknsIp/4gUAAIbjn/sBAGCH6yRLkcldJCMEAAAMRyAEAAAMR2kcAADseJ29rTRukYwQAAAwHBkhAADY4ZYXS/D8r0UyQgAAwHAEQgAAwHCUxgEAwAA8R2glGSEAAGA4MkIAALDDdTp722IJi2SEAACA4QiEAACA4SiNAwCAAXiO0EoyQgAAwHBkhAAAYIfrJHtlhFaQEQIAAIYjEAIAAIajNA4AAAZgsYSVZIQAAIDhyAgBAMAO10n2tozQIhkhAABgOAIhAABgOErjAABgAEtzT2CbkRECAACGIxACAACGozQOAAB2uE5nr+cIrSAjBAAADEdGCAAAdrpO9koIrSAjBAAADEcgBAAADEdpHAAA7HAdzxFaTUYIAAAYjowQAADseJW9qbknsa3ICAEAAMMRCAEAAMNRGgcAADtcJ1nyHKEVZIQAAIDhyAgBAMAALJawkowQAAAwHIEQAAAwHKVxAACww3WUxq0mIwQAAAxHRggAAAaw1DJCi2SEAACA2VXVMVX17qq6pKourqrnrNHnblX1P6vqw1Ofn1zYd2pVfWzaTt3sfDJCAADAdrAnyfO6+4KqOjzJh6rqnO6+ZKHPzyW5pLt/oKruleSyqnpTksOSvCTJiVm+JepDVXVmd39mvZMJhAAAYIe7PSyW0N3XJLlmen1TVV2a5Kgki4FQJzm8qirLwc8NWQ6gnpzknO6+IUmq6pwkJyf5g/XOJxACAAC2wpFVdf7C+9O6+7S1OlbVcUkekeTcVbtemeTMJH+b5PAkP9LdS1V1VJIrF/pdleUgal0CIQAA2OE6lb3zLw9wXXefuFmnqjosyVuSPLe7b1y1+8lJLkzyPUm+Ick5VfW+WzOZ2T8NAACAJKmqg7McBL2pu89Yo8tPJjmjl+1OckWSb0xydZJjFvodPbWtSyAEAADMbrrv57VJLu3uV6zT7VNJnjD1v0+ShyS5PMnZSZ5UVUdU1RFJnjS1rUtpHAAADOB28ByhRyd5RpKLqurCqe1FSY5Nku5+dZJfTvL6qrooSSV5fndflyRV9ctJzpuOe+nNCyesRyAEAADMrrvfn2y8tF13/22Wsz1r7Xtdktft6/mUxgEAAMOREQIAgB3u9vAcoa0mIwQAAAxHRggAAHa8yt6WA1nk0wAAAIYjEAIAAIYzayBUVSdX1WVVtbuqXrDG/kOr6s3T/nOr6riFfS+c2i+rqiffgjF/o6o+f8AuCgAAtplOspQ7zLptN7PNqKp2JfnNJN+X5KFJnl5VD13V7ZlJPtPdD0rya0leNh370CSnJHlYkpOTvKqqdm02ZlWdmOSIA3phAADAtjdnaHZSkt3dfXl3fyXJ6UmeuqrPU5O8YXr9x0meUFU1tZ/e3V/u7iuS7J7GW3fMKUh6eZJ/f4CvCwAAtp29qVm37WbOVeOOSnLlwvurknznen26e09VfS7JPaf2D6w69qjp9XpjPjvJmd19zXIstblDP/HFnPATF+1T351q133uPfcUZrd0xF3nnsK28NXD7jT3FGb3xXtvv7/Et9pNxx469xS2hcM/tTT3FGZ38E099xRmt+vLvgeHfviKuacAt9oQy2dX1f2T/HCSx+1D32cleVaS3DF3PrATAwAAZjFnIHR1kmMW3h89ta3V56qqOijJ3ZJcv8mxa7U/IsmDkuyeskF3rqrd071HK3T3aUlOS5K73uEe/rkLAIDbvW7PEVptzk/jvCQnVNXxVXVIlhc/OHNVnzOTnDq9flqSP+/untpPmVaVOz7JCUk+uN6Y3f2n3X3f7j6uu49L8g9rBUEAAMAYZssITff8PDvJ2Ul2JXldd19cVS9Ncn53n5nktUneWFW7k9yQ5cAmU78/THJJkj1Jfq679ybJWmNu9bUBAMB2s7QNFyyY06z3CHX3WUnOWtX2iwuvv5Tle3vWOvZXkvzKvoy5Rp/Dbs18AQCAnUGhIAAAMJwhVo0DAICRdZK9ciAr+DQAAIDhyAgBAMCOZ/ns1XwaAADAcARCAADAcJTGAQDADtdJluRAVvBpAAAAw5ERAgCAAeztmnsK24qMEAAAMByBEAAAMBylcQAAsMN1KnvlQFbwaQAAAMMRCAEAAMNRGgcAAANYajmQRT4NAABgODJCAACww3VisYRVfBoAAMBwBEIAAMBwlMYBAMAO16ns7Zp7GtuKjBAAADAcGSEAABjAkhzICj4NAABgOAIhAABgOErjAABgh+tO9rYcyCKfBgAAMBwZIQAA2PEqS7F89iIZIQAAYDgCIQAAYDhK4wAAYIfrWCxhNZ8GAAAwHBkhAAAYwF45kBV8GgAAwHAEQgAAwHCUxgEAwA7XqSy15wgtkhECAACGIxACAACGozQOAAAGYNW4lXwaAADAcGSEAABgh+skSy0HssinAQAADEcgBAAADEdpHAAA7HiVvfEcoUUyQgAAwHBkhAAAYIezWMI/5tMAAACGIxACAACGozQOAAAGYLGElWSEAACA4cgIAQDADtddFktYxacBAAAMRyAEAAAMR2kcAAAMYK/SuBV8GgAAwHBkhAAAYIfrJEuWz15BRggAABiOjBAAADC7qjomye8muU+Wk1indfd/W9XnF5L82PT2oCTflORe3X1DVX0iyU1J9ibZ090nbnQ+gRAAAOx4dXtYLGFPkud19wVVdXiSD1XVOd19yc0duvvlSV6eJFX1A0n+TXffsDDG47v7un052bb/NAAAgJ2vu6/p7gum1zcluTTJURsc8vQkf3BrzycjtIEvH3vnXPbiR8w9jVkd9BlfkfuctzT3FLaFY865ce4pzO5L97rT3FOY3SE3fnXuKWwLH/8ZNxzf4ao7zj2F2R3yOd+Dw+/34LmnsD383twT2FwnWerbz3e2qo5L8ogk566z/85JTk7y7IXmTvLOquokv93dp210Dr/lAgAAW+HIqjp/4f1pawUrVXVYkrckeW53r/cvsT+Q5C9XlcU9pruvrqp7Jzmnqv5Pd793vckIhAAAgK1w3WYLGFTVwVkOgt7U3Wds0PWUrCqL6+6rpz+vrao/SXJSEoEQAACMbO82Xx6gqirJa5Nc2t2v2KDf3ZJ8d5L/d6HtLknu0N03Ta+flOSlG51PIAQAAGwHj07yjCQXVdWFU9uLkhybJN396qnth5K8s7u/sHDsfZL8yXIslYOS/H53v2OjkwmEAACA2XX3+5NsuqJDd78+yetXtV2e5NtuyfkEQgAAsMN16na1atxW2N6FggAAAAeAjBAAAAxgSQ5kBZ8GAAAwHIEQAAAwHKVxAACww3Uney2WsIKMEAAAMBwZIQAAGIDls1eSEQIAAIYjEAIAAIajNA4AAHa4TmWp5UAW+TQAAIDhyAgBAMAA9sZiCYtkhAAAgOEIhAAAgOEojQMAgB2u4zlCq8kIAQAAw5ERAgCAHc/y2av5NAAAgOEIhAAAgOEojQMAgAEseY7QCjJCAADAcGSEAABgh+tO9lo+ewUZIQAAYDgCIQAAYDhK4wAAYACeI7SSTwMAABiOQAgAABiO0jgAANjhOpUlq8atICMEAAAMR0YIAAAGsBQZoUUyQgAAwHAEQgAAwHCUxgEAwA7XicUSVpERAgAAhiMjBAAAA1hqOZBFPg0AAGA4AiEAAGA4SuMAAGCn67JYwioyQgAAwHBkhAAAYIfrJEuREVokIwQAAAxHIAQAAAxHaRwAAAzAYgkryQgBAADDkRECAIAdriMjtJqMEAAAMByBEAAAMBylcQAAMAClcSvJCAEAAMOREQIAgB2uUzJCq8gIAQAAwxEIAQAAw1EaBwAAA1iK0rhFMkIAAMBwBEIAAMBwZg2EqurkqrqsqnZX1QvW2H9oVb152n9uVR23sO+FU/tlVfXkzcasqjdN7R+tqtdV1cEH/AIBAGA76OXnCM25bTezBUJVtSvJbyb5viQPTfL0qnroqm7PTPKZ7n5Qkl9L8rLp2IcmOSXJw5KcnORVVbVrkzHflOQbk3xLkjsl+akDeHkAAMA2NudiCScl2d3dlydJVZ2e5KlJLlno89QkvzS9/uMkr6yqmtpP7+4vJ7miqnZP42W9Mbv7rJsHraoPJjn6QF0YAABsJ51sy6zMnOYsjTsqyZUL76+a2tbs0917knwuyT03OHbTMaeSuGckecdtvgIAAOB2acTls1+V5L3d/b61dlbVs5I8K0kOutsRudNVY99KdL+//PLcU5jdIdd+Ye4pbAuf/q4j5p7C7G56wNwzmN/x//Orc09he7j20LlnMLs7HO/vxmN/+StzT2F2n/yhe8w9BbjV5gyErk5yzML7o6e2tfpcVVUHJblbkus3OXbdMavqJUnuleRfrjep7j4tyWlJcsf7H9P7fjkAALB9KY1bac7SuPOSnFBVx1fVIVle/ODMVX3OTHLq9PppSf68u3tqP2VaVe74JCck+eBGY1bVTyV5cpKnd/fSAb42AABgG5stI9Tde6rq2UnOTrIryeu6++KqemmS87v7zCSvTfLGaTGEG7Ic2GTq94dZXlhhT5Kf6+69SbLWmNMpX53kk0n+enm9hZzR3S/dossFAIDZdLbnEtZzmvUeoWklt7NWtf3iwusvJfnhdY79lSS/si9jTu0j3g8FAACsYdYHqgIAAMxBlgQAAAbQSuNWkBECAACGIyMEAAADWIqM0CIZIQAAYDgCIQAAYDgCIQAA2OG6k6WuWbfNVNUxVfXuqrqkqi6uques0ecXqurCaftoVe2tqntM+06uqsuqandVvWCz8wmEAACA7WBPkud190OTPCrJz1XVQxc7dPfLu/vh3f3wJC9M8hfdfUNV7Urym0m+L8lDkzx99bGrWSwBAAAGsN2Xz+7ua5JcM72+qaouTXJUkkvWOeTpSf5gen1Skt3dfXmSVNXpSZ66wbEyQgAAwPZSVccleUSSc9fZf+ckJyd5y9R0VJIrF7pcNbWtS0YIAADYCkdW1fkL70/r7tNWd6qqw7Ic4Dy3u29cZ6wfSPKX3X3DrZ2MQAgAAHa8fVuw4AC7rrtP3KhDVR2c5SDoTd19xgZdT8nXy+KS5Ookxyy8P3pqW5fSOAAAYHZVVUlem+TS7n7FBv3uluS7k7xtofm8JCdU1fFVdUiWA6UzNzqfjBAAALAdPDrJM5JcVFUXTm0vSnJsknT3q6e2H0ryzu7+ws0Hdveeqnp2krOT7Eryuu6+eKOTCYQAAGAAt4NV496fZNNJdvfrk7x+jfazkpy1r+dTGgcAAAxHRggAAHa4TrbDYgnbiowQAAAwHIEQAAAwHKVxAACw03XSPfckthcZIQAAYDgyQgAAMIClzVemHoqMEAAAMByBEAAAMBylcQAAsMN1kvYcoRVkhAAAgOHICAEAwI5XWZIRWkFGCAAAGI5ACAAAGI7SOAAAGED33DPYXmSEAACA4cgIAQDAACyfvZKMEAAAMByBEAAAMBylcQAAsMN1K41bTUYIAAAYjowQAAAMYElGaAUZIQAAYDgCIQAAYDhK4wAAYADdc89ge5ERAgAAhiMQAgAAhqM0DgAABuA5QivJCAEAAMOREQIAgB2uUzJCq8gIAQAAwxEIAQAAw1EaBwAAA/AYoZVkhAAAgOHICAEAwE7Xls9eTUYIAAAYjkAIAAAYjtI4AAAYgdUSVpARAgAAhiMjBAAAA7BYwkoyQgAAwHAEQgAAwHCUxgEAwADaYgkryAgBAADDkRECAIAdrmOxhNUEQhupZO8hY+cQd31579xTmN3lP3LE3FPYFo5/201zT2F29zvrc3NPYXZLh99l7ilsC0t397/Po3/74LmnMLs73PSZuacwu6P//I5zT2FbuGzuCXCrKI0DAACG45+0AABgp+skSuNWkBECAACGIyMEAAADsHz2SjJCAADAcARCAADAcJTGAQDACJTGrSAjBAAADEcgBAAADEdpHAAA7HiV9hyhFWSEAACA4cgIAQDACCyWsIKMEAAAMByBEAAAMBylcQAAsNN1LJawiowQAAAwHBkhAAAYgcUSVpARAgAAhiMQAgAAhqM0DgAAhmCxhEUyQgAAwHBkhAAAYAQWS1hBRggAABiOQAgAABiO0jgAABiB0rgVZIQAAIDZVdUxVfXuqrqkqi6uques0+9xVXXh1OcvFto/UVUXTfvO3+x8MkIAALDTdZLe9stn70nyvO6+oKoOT/Khqjqnuy+5uUNV3T3Jq5Kc3N2fqqp7rxrj8d193b6cTEYIAACYXXdf090XTK9vSnJpkqNWdfvRJGd096emftfe2vMJhAAAgG2lqo5L8ogk567a9eAkR1TVe6rqQ1X14wv7Osk7p/ZnbXYOpXEAADCAnn+xhCNX3btzWneftrpTVR2W5C1JntvdN67afVCS70jyhCR3SvLXVfWB7v6bJI/p7quncrlzqur/dPd715uMQAgAANgK13X3iRt1qKqDsxwEvam7z1ijy1VJru/uLyT5QlW9N8m3Jfmb7r46WS6Xq6o/SXJSknUDIaVxAADA7Kqqkrw2yaXd/Yp1ur0tyWOq6qCqunOS70xyaVXdZVpgIVV1lyRPSvLRjc4nIwQAACOYvzRuM49O8owkF1XVhVPbi5IcmyTd/eruvrSq3pHkI0mWkrymuz9aVQ9M8ifLsVQOSvL73f2OjU4mEAIAAGbX3e9Psuka39398iQvX9V2eZZL5PaZQAgAAEaw/Z8jtKXcIwQAAAxHIAQAAAxHaRwAAAygtv9iCVtKRggAABiOjBAAAOx0ndvD8tlbSkYIAAAYjkAIAAAYjtI4AADY8cpzhFaREQIAAIYjIwQAACOwWMIKMkIAAMBwBEIAAMBwlMYBAMAIlMatICMEAAAMR0YIAABGICO0gowQAAAwHIEQAAAwHKVxAACw03WSrrlnsa1sGghV1R2TfH+S70py/yRfTPLRJH/a3Rcf2OkBAADsfxsGQlX1H7McBL0nyblJrk1yxyQPTvKrU5D0vO7+yAGeJwAAcBuUxRJW2Cwj9MHufsk6+15RVfdOcux+nhMAAMABtWEg1N1/usn+a7OcJQIAALjd2KfFEqrqxCQvTvKA6ZhK0t39rQdwbgAAwP6iNG6FfV017k1JfiHJRUmWDtx0AAAADrx9DYT+vrvPPKAzAQAA2CL7Ggi9pKpek+RdSb58c2N3n3FAZgUAAHAA3WEf+/1kkocnOTnJD0zb99/Wk1fVyVV1WVXtrqoXrLH/0Kp687T/3Ko6bmHfC6f2y6rqyZuNWVXHT2PsnsY85LbOHwAAuH3a14zQI7v7IfvzxFW1K8lvJnlikquSnFdVZ3b3JQvdnpnkM939oKo6JcnLkvxIVT00ySlJHpblh7z+WVU9eDpmvTFfluTXuvv0qnr1NPZv7c9rAgCA7cpzhFba14zQX03Bx/50UpLd3X15d38lyelJnrqqz1OTvGF6/cdJnlBVNbWf3t1f7u4rkuyexltzzOmY75nGyDTmD+7n6wEAAG4n9jUj9KgkF1bVFVm+R2h/LJ99VJIrF95fleQ71+vT3Xuq6nNJ7jm1f2DVsUdNr9ca855JPtvde9bov64+dCl7jv/SPl3MTnWHL+/ZvNMOt+cu/vkkST77kMPmnsLs7nb+R+eewux23efec09hW7jPu+4+9xRm94X7zT2D+V31+KPnnsLs7nRtzT2F7eH9c0+AW2NfA6GTD+gstpGqelaSZyXJriPvNvNsAABgP2mB66INA6GqOqy7P9/dn9ysz60499VJjll4f/TUtlafq6rqoCR3S3L9Jseu1X59krtX1UFTVmitcyVJuvu0JKclyaEPPEoqAAAAdqDN7hF6W1X916p6bFXd5ebGqnpgVT2zqs7Orc8WnZfkhGk1t0OyvPjB6mcVnZnk1On105L8eXf31H7KtKrc8UlOSPLB9cacjnn3NEamMd92K+cNAAC3L70Ntm1mw4xQdz+hqp6S5F8meXRVHZFkT5LLkvxpklO7++9uzYmne36eneTsJLuSvK67L66qlyY5f3qA62uTvLGqdie5IcuBTaZ+f5jkkmk+P9fde5NkrTGnUz4/yelV9Z+S/O9pbAAAYECb3iPU3WclOetAnHytsbv7FxdefynJD69z7K8k+ZV9GXNqvzzLq8oBAACD29fFEgAAgNuzbVieNqd9fY4QAADAjrFhIFRVZ1XVcVs0FwAA4ACpnnfbbjbLCP2PJO+sqhdX1cFbMSEAAIADbbNV4/6oqt6e5D8kOb+q3phkaWH/Kw7w/AAAAPa7fVks4StJvpDk0CSHZyEQAgAAbie2YXnanDYMhKrq5CSvyPIDTL+9u/9hS2YFAABwAG2WEXpxkh9eeCgpAABweyQjtMJm9wh911ZNBAAAYKt4jhAAADCcfVksAQAAuB3brs/ymZOMEAAAMBwZIQAAGEHX3DPYVmSEAACA4QiEAACA4SiNAwCAEVgsYQUZIQAAYDgCIQAAYDhK4wAAYACeI7SSjBAAADAcGSEAABiBjNAKMkIAAMBwBEIAAMBwlMYBAMBO1xZLWE1GCAAAGI6MEAAAjEBGaAUZIQAAYDgCIQAAYDhK4wAAYARK41aQEQIAAIYjIwQAAAOwfPZKMkIAAMBwBEIAAMBwBEIAAMBwBEIAAMBwBEIAADCCnnnbRFUdU1XvrqpLquriqnrOOv0eV1UXTn3+YqH95Kq6rKp2V9ULNjufVeMAAIDtYE+S53X3BVV1eJIPVdU53X3JzR2q6u5JXpXk5O7+VFXde2rfleQ3kzwxyVVJzquqMxePXU1GCAAAmF13X9PdF0yvb0pyaZKjVnX70SRndPenpn7XTu0nJdnd3Zd391eSnJ7kqRudTyAEAAA7XS8/R2jO7ZaoquOSPCLJuat2PTjJEVX1nqr6UFX9+NR+VJIrF/pdlX8cRK2gNA4AANgKR1bV+QvvT+vu01Z3qqrDkrwlyXO7+8ZVuw9K8h1JnpDkTkn+uqo+cGsmIxACAAC2wnXdfeJGHarq4CwHQW/q7jPW6HJVkuu7+wtJvlBV703ybVP7MQv9jk5y9UbnUhoHAAAj2P6rxlWS1ya5tLtfsU63tyV5TFUdVFV3TvKdWb6X6LwkJ1TV8VV1SJJTkpy50flkhAAAgO3g0UmekeSiqrpwantRkmOTpLtf3d2XVtU7knwkyVKS13T3R5Okqp6d5Owku5K8rrsv3uhkAiEAABjBLVywYKt19/uT1D70e3mSl6/RflaSs/b1fErjAACA4QiEAACA4SiNAwCAHa5yy5/ls9PJCAEAAMOREQIAgBHICK0gIwQAAAxHIAQAAAxHaRwAAOx0bbGE1WSEAACA4cgIAQDACGSEVpARAgAAhiMQAgAAhqM0DgAARqA0bgUZIQAAYDgyQgAAMADLZ68kIwQAAAxHRmgDh1xfOeZNY39Ed7jxxrmnML97Hzz3DLaFTz9h19xTmN3eQ/7J3FOY3V0/+ZW5p7AtfPWwmnsKs7v7x3wXPvDy1849hdk9+rn/au4pwK029m/5AAAwCqVxKyiNAwAAhiMjBAAAO11HRmgVGSEAAGA4AiEAAGA4SuMAAGAAniO0kowQAAAwHIEQAAAwHKVxAAAwAqVxK8gIAQAAw5ERAgCAAVgsYSUZIQAAYDgCIQAAYDhK4wAAYARK41aQEQIAAIYjIwQAADtdR0ZoFRkhAABgOAIhAABgOErjAABgh6tp4+tkhAAAgOHICAEAwAgslrCCjBAAADAcgRAAADAcpXEAADCAUhq3gowQAAAwHBkhAAAYgYzQCjJCAADAcARCAADAcJTGAQDACJTGrSAjBAAADEdGCAAAdrq2fPZqMkIAAMBwBEIAAMBwlMYBAMAIlMatICMEAAAMRyAEAAAMR2kcAAAMwKpxK8kIAQAAw5ERAgCAEcgIrSAjBAAADEcgBAAADEdpHAAADMBiCSvJCAEAAMOREQIAgJ2uY7GEVWSEAACA4QiEAACA4SiNAwCAESiNW0FGCAAAGI6MEAAA7HAVy2evJiMEAAAMRyAEAAAMR2kcAACMQGncCjJCAADAcGSEAABgANVSQotkhAAAgOEIhAAAgNlV1TFV9e6quqSqLq6q56zR53FV9bmqunDafnFh3yeq6qKp/fzNzqc0DgAAdrrO7WGxhD1JntfdF1TV4Uk+VFXndPclq/q9r7u/f50xHt/d1+3LyWSEAACA2XX3Nd19wfT6piSXJjnqQJ1PIAQAAAOonndLcmRVnb+wPWvduVYdl+QRSc5dY/c/qaoPV9Xbq+phC+2d5J1V9aGNxr6Z0jgAAGArXNfdJ27WqaoOS/KWJM/t7htX7b4gyQO6+/NV9ZQkb01ywrTvMd19dVXdO8k5VfV/uvu9651nloxQVd2jqs6pqo9Nfx6xTr9Tpz4fq6pTF9q/Y7oRandV/UZV1UbjVtWPVdVHpmP+qqq+bWuuFAAA2FdVdXCWg6A3dfcZq/d3943d/fnp9VlJDq6qI6f3V09/XpvkT5KctNG55iqNe0GSd3X3CUneNb1foarukeQlSb4zyxfxkoWA6beS/HSWo78Tkpy8ybhXJPnu7v6WJL+c5LQDcVEAALBt9czbJqbkxmuTXNrdr1inz30XkiAnZTmeub6q7jItsJCqukuSJyX56Ebnm6s07qlJHje9fkOS9yR5/qo+T05yTnffkCRVdU6Sk6vqPUnu2t0fmNp/N8kPJnn7euN2918tjPuBJEfvx2sBAABuu0cneUaSi6rqwqntRUmOTZLufnWSpyX5marak+SLSU7p7q6q+yT5kylGOijJ73f3OzY62VyB0H26+5rp9d8luc8afY5KcuXC+6umtqOm16vb93XcZ2Y5aAIAALaJ7n5/ktqkzyuTvHKN9suT3KLbXw5YIFRVf5bkvmvsevHimymC2++rmq81blU9PsuB0GPWO25aYeJZSXLQ3Y/I1d899noSB5+4Viw5loe89NNzT2FbWLrrneaewuw++5Cx/z5Ikst/xGKjSfLAN39l7inM7ktHHjz3FGb3lG987NxTmN2Xnr7h76xsM/v/N+7btwP2f/Xu/t719lXVp6vqft19TVXdL8m1a3S7Ol8vc0uWy9neM7Ufvar96un1uuNW1bcmeU2S7+vu6zeY92mZ7iE69OhjfF0AAGAHmuuf9s5McvMqcKcmedsafc5O8qSqOmJaJOFJSc6eSt9urKpHTTdK/fjC8WuOW1XHJjkjyTO6+28OxAUBAMC2ts0XS9hqcwVCv5rkiVX1sSTfO71PVZ1YVa9JkmmRhF9Oct60vfTmhROS/GyWszu7k3w8X7/nZ81xk/xiknsmeVVVXVhV5x/g6wMAALaxWQrep9K0J6zRfn6Sn1p4/7okr1un3zffgnF/anFcAABgbO78BQCAna4tlrCa5X8AAIDhyAgBAMAIZIRWkBECAACGIxACAACGozQOAAB2uIrFElaTEQIAAIYjIwQAACNoKaFFMkIAAMBwBEIAAMBwlMYBAMAALJawkowQAAAwHBkhAADY6Xra+BoZIQAAYDgCIQAAYDhK4wAAYAC1NPcMthcZIQAAYDgCIQAAYDhK4wAAYARWjVtBRggAABiOjBAAAAygZIRWkBECAACGIxACAACGozQOAAB2uk7SauMWyQgBAADDkRECAIABWCxhJRkhAABgOAIhAABgOErjAABgBErjVpARAgAAhiMjBAAAO1zFYgmryQgBAADDEQgBAADDURoHAAA7XffyxtfICAEAAMOREQIAgAFYLGElGSEAAGA4AiEAAGA4SuMAAGAESuNWkBECAACGIyMEAAADsFjCSjJCAADAcARCAADAcJTGAQDATtdJltTGLZIRAgAAhiMQAgAAhqM0DgAARqAybgUZIQAAYDgyQgAAMADPEVpJRggAABiOQAgAABiO0jgAABhBq41bJCMEAAAMR0YIAAAGYLGElWSEAACA4QiEAACA4SiNAwCAna6nja+REQIAAIYjIwQAADtcJSnLZ68gIwQAAAxHRmgDB3+hc78P7J17GrP64hG75p7C7D7/0HvOPYVt4eAbx/5vIUm+eljNPYXZnfC7X517CtvCtd9+57mnMLujzrxy7inM7v/80kPnnsLsHvKyj889BbjVBEIAADCCpbknsL0ojQMAAIYjIwQAAAOwWMJKMkIAAMBwBEIAAMBwBEIAALDT9TbYNlFVx1TVu6vqkqq6uKqes0afx1XV56rqwmn7xYV9J1fVZVW1u6pesNn53CMEAABsB3uSPK+7L6iqw5N8qKrO6e5LVvV7X3d//2JDVe1K8ptJnpjkqiTnVdWZaxz7NQIhAADY8TrZ5osldPc1Sa6ZXt9UVZcmOSrJusHMgpOS7O7uy5Okqk5P8tSNjlUaBwAAbCtVdVySRyQ5d43d/6SqPlxVb6+qh01tRyVZfNLzVVPbumSEAACArXBkVZ2/8P607j5tdaeqOizJW5I8t7tvXLX7giQP6O7PV9VTkrw1yQm3ZjICIQAAGEDNXxl3XXefuFGHqjo4y0HQm7r7jNX7FwOj7j6rql5VVUcmuTrJMQtdj57a1qU0DgAAmF1VVZLXJrm0u1+xTp/7Tv1SVSdlOZ65Psl5SU6oquOr6pAkpyQ5c6PzyQgBAADbwaOTPCPJRVV14dT2oiTHJkl3vzrJ05L8TFXtSfLFJKd0dyfZU1XPTnJ2kl1JXtfdF290MoEQAACMYPuvGvf+JLVJn1cmeeU6+85Kcta+nk9pHAAAMBwZIQAA2Ok6qaW5J7G9yAgBAADDEQgBAADDURoHAAAj2OaLJWw1GSEAAGA4MkIAADACCaEVZIQAAIDhCIQAAIDhKI0DAIABlMUSVpARAgAAhiMjBAAAI5ARWkFGCAAAGI5ACAAAGI7SOAAA2Ok6ydLck9heZIQAAIDhyAgBAMAOV2nLZ68iIwQAAAxHIAQAAAxHaRwAAIxAadwKMkIAAMBwBEIAAMBwlMYBAMAIlMatICMEAAAMR0YIAAB2uk6yNPckthcZIQAAYDgCIQAAYDhK4wAAYABlsYQVZIQAAIDhyAgBAMAIZIRWkBECAACGIxACAACGozQOAAB2vFYat4qMEAAAMBwZIQAA2Ok6MkKryAgBAADDEQgBAADDURoHAAAjWJp7AtuLjBAAADAcGSEAABhAWSxhBRkhAABgOAIhAABgOErjAABgBErjVpARAgAAhiMjBAAAO10nWZIRWiQjBAAADEcgBAAADGeWQKiq7lFV51TVx6Y/j1in36lTn49V1akL7d9RVRdV1e6q+o2qqn0Zt6oeWVV7quppB/YKAQBgO+nlxRLm3LaZuTJCL0jyru4+Icm7pvcrVNU9krwkyXcmOSnJSxYCm99K8tNJTpi2kzcbt6p2JXlZknceiAsCAABuP+YKhJ6a5A3T6zck+cE1+jw5yTndfUN3fybJOUlOrqr7Jblrd3+guzvJ7y4cv9G4P5/kLUmu3X+XAQAA3B7NtWrcfbr7mun13yW5zxp9jkpy5cL7q6a2o6bXq9vXHbeqjkryQ0ken+SR++MCAADgdmUblqfN6YAFQlX1Z0nuu8auFy++6e6uqv3+U1k17q8neX53L023E62rqp6V5FlJsuuII3L1Y8deT+IuV2/8eY3gHpd+ae4pbAvXPPqwuacwu4M/738gVz/2znNPYVtoD5/Itd9z9NxTmN1R79079xRm9/dP+Ya5p7A9vG7uCXBrHLC/yrv7e9fbV1Wfrqr7dfc1U6nbWuVqVyd53ML7o5O8Z2o/elX71dPr9cY9McnpUxB0ZJKnVNWe7n7rGvM+LclpSXLoMcf4rQcAgJ1BRmiFudIdZya5eRW4U5O8bY0+Zyd5UlUdMS2S8KQkZ0+lbzdW1aOm1eJ+fOH4Ncft7uO7+7juPi7JHyf52bWCIAAAYAxzBUK/muSJVfWxJN87vU9VnVhVr0mS7r4hyS8nOW/aXjq1JcnPJnlNkt1JPp7k7RuNCwAAsGiWKufuvj7JE9ZoPz/JTy28f13WqLqc+n3zvo67qs9P3PIZAwDA7VgnWVIat2jslQAAAIAhWfcGAAB2vE56ae5JbCsyQgAAwHAEQgAAwHCUxgEAwAg8R2gFGSEAAGA4MkIAALDTWT77H5ERAgAAhiMQAgAAhqM0DgAARmCxhBVkhAAAgOHICAEAwAhkhFaQEQIAAIYjEAIAAIajNA4AAHa8Vhq3iowQAAAwu6o6pqreXVWXVNXFVfWcDfo+sqr2VNXTFtr2VtWF03bmZueTEQIAgJ2ukywtzT2LzexJ8rzuvqCqDk/yoao6p7svWexUVbuSvCzJO1cd/8Xufvi+nkxGCAAAmF13X9PdF0yvb0pyaZKj1uj680nekuTa23I+gRAAALCtVNVxSR6R5NxV7Ucl+aEkv7XGYXesqvOr6gNV9YObnUNpHAAAjGD+xRKOrKrzF96f1t2nre5UVYdlOePz3O6+cdXuX0/y/O5eqqrVhz6gu6+uqgcm+fOquqi7P77eZARCAADAVriuu0/cqENVHZzlIOhN3X3GGl1OTHL6FAQdmeQpVbWnu9/a3VcnSXdfXlXvyXJGad1ASGkcAAAwu1qObl6b5NLufsVafbr7+O4+rruPS/LHSX62u99aVUdU1aHTOEcmeXSSS9Ya42YyQgAAMIL5S+M28+gkz0hyUVVdOLW9KMmxSdLdr97g2G9K8ttVtZTlZM+vrl5tbjWBEAAAMLvufn+Sf3Tjzwb9f2Lh9V8l+ZZbcj6BEAAA7HidLG37jNCWco8QAAAwHIEQAAAwHKVxAACw03XSvTT3LLYVGSEAAGA4MkIAADACiyWsICMEAAAMRyAEAAAMR2kcAACMoJXGLZIRAgAAhiMjBAAAO113smT57EUyQgAAwHAEQgAAwHCUxgEAwAgslrCCjBAAADAcGSEAABhAWyxhBRkhAABgOAIhAABgOErjAABgx2uLJawiIwQAAAxHIAQAAAxHaRwAAOx0nWRJadwiGSEAAGA4MkIAADCC9hyhRTJCAADAcARCAADAcJTGAQDADtdJ2mIJK8gIAQAAw5ERAgCAna7bYgmryAgBAADDEQgBAADDURoHAAADsFjCSjJCAADAcGSEAABgBBZLWEFGCAAAGI5ACAAAGE51u2lqPVX190k+Ofc8ZnZkkuvmngTbgu8Cie8BX+e7QOJ7cLMHdPe95p7ERqrqHVn+ec3puu4+eeY5fI1AiA1V1fndfeLc82B+vgskvgd8ne8Cie8Bt29K4wAAgOEIhAAAgOEIhNjMaXNPgG3Dd4HE94Cv810g8T3gdsw9QgAAwHBkhAAAgOEIhAZSVfeoqnOq6mPTn0es0+/Uqc/HqurUhfbvqKqLqmp3Vf1GVdW+jFtVj6yqPVX1tAN7heyrrf4uVNWPVdVHpmP+qqq+bWuulLVU1clVddn083vBGvsPrao3T/vPrarjFva9cGq/rKqevNmYVXX8NMbuacxDDvgFsk+2+Hvwpqn9o1X1uqo6+IBfIPtsK78LC/t/o6o+f8AuCvZFd9sG2ZL85yQvmF6/IMnL1uhzjySXT38eMb0+Ytr3wSSPSlJJ3p7k+zYbN8muJH+e5KwkT5v7M7DN811I8n8tHPt9Sc6d+zMYdZv+m/x4kgcmOSTJh5M8dFWfn03y6un1KUnePL1+6NT/0CTHT+Ps2mjMJH+Y5JTp9auT/Mzcn4Ftlu/BU6a/LyrJH/gebJ9tq78L03EnJnljks/Pff22sTcZobE8NckbptdvSPKDa/R5cpJzuvuG7v5MknOSnFxV90ty1+7+QHd3kt9dOH6jcX8+yVuSXLv/LoP9YEu/C939V9MYSfKBJEfv16vhljgpye7uvry7v5Lk9Cz/3BYt/hz/OMkTpqzfU5Oc3t1f7u4rkuyexltzzOmY75nGSNb/rrH1tux7kCTdfVZPsvwPKf4O2D629LtQVbuSvDzJvz/A1wWbEgiN5T7dfc30+u+S3GeNPkcluXLh/VVT21HT69Xt645bVUcl+aEkv7VfZs/+tKXfhVWemeUsEvNY7+e6Zp/u3pPkc0nuucGx67XfM8lnpzHWOxfz2MrvwddMJXHPSPKO23wF7C9b/V14dpIzF/5fAbM5aO4JsH9V1Z8lue8au168+Ka7u6r2+5KBq8b99STP7+6l6RYSttA2+y7cPKfHZzkQesz+Ph9wu/CqJO/t7vfNPRG2XlXdP8kPJ3nczFOBJAKhHae7v3e9fVX16aq6X3dfM5U3rVWudnVW/gV1dJL3TO1Hr2q/enq93rgnJjl9CoKOTPKUqtrT3W+9xRfGLbbNvgupqm9N8pos3090/a24JPaPq5Mcs/B+8ee3us9VVXVQkrsluX6TY9dqvz7J3avqoOlfkdc6F/PYyu9BkqSqXpLkXkn+5X6YP/vPVn4XHpHkQUl2T78b3Lmqdnf3g/bPpcAtozRuLGcmuXnlr1OTvG2NPmcneVJVHTGt+PWkJGdPKewbq+pRU13wjy8cv+a43X18dx/X3cdluab4ZwVB28aWfheq6tgkZyR5Rnf/zYG4IPbZeUlOmFZzOyTLNz6fuarP4s/xaUn+fLq348wkp0wrSB2f5IQs3++x5pjTMe+exkjW/66x9bbse5AkVfVTWb7v8OndvXSAr41bZiv/TvjT7r7vwu8G/yAIYlZzr9Zg27oty/W870rysSR/luQeU/uJSV6z0O9fZPmGx91JfnKh/cQkH83ySjCvzNcfyLvmuKvO/fpYNW7bbFv9XchyJugzSS6ctvPn/gxG3rK8gtffTD+/F09tL03yf0+v75jkj6af+weTPHDh2BdPx12WabXA9cac2h84jbF7GvPQua/fNsv3YM/UdvPfAb849/Xb5vkurDqvVeNss243//ICAAAwDKVxAADAcARCAADAcARCAADAcARCAADAcARCAADAcARCAIOqqmOq6oqqusf0/ojp/XFr9L1TVf1FVe26BeM/u6r+xX6cMgDsN5bPBhhYVf37JA/q7mdV1W8n+UR3/39r9Pu5JAd193+7BWPfOclfdvcj9t+MAWD/kBECGNuvJXlUVT03yWOS/Jd1+v1YkrclSVU9bsoOva2qLq+qX62qH6uqD1bVRVX1DUnS3f+Q5BNVddIWXAcA3CICIYCBdfdXk/xClgOi507vV6iqQ7L8JPlPLDR/W5J/leSbkjwjyYO7+6Qkr0ny8wv9zk/yXQdm9gBw6wmEAPi+JNck+eZ19h+Z5LOr2s7r7mu6+8tJPp7knVP7RUmOW+h3bZL777eZAsB+IhACGFhVPTzJE5M8Ksm/qar7rdHti0nuuKrtywuvlxbeLyU5aGHfHafjAWBbEQgBDKqqKslvZbkk7lNJXp417hHq7s8k2VVVq4OhffHgJB+9TRMFgANAIAQwrp9O8qnuPmd6/6ok31RV371G33dmeTGFW+rRSc7ZtBcAbDHLZwOwqar69iT/prufcQuOeUSSf3tLjgGArSIjBMCmuvuCJO++JQ9UzfIiC//hAE0JAG4TGSEAAGA4MkIAAMBwBEIAAMBwBEIAAMBwBEIAAMBwBEIAAMBwBEIAAMBw/n86AtkFvGbESAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x864 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "fig=plt.figure(figsize=[15,12])\n",
    "plt.imshow(max_radiance_w_sr_m2.get(),# vmin=0,\n",
    "           extent=(edges[0][0].item(), edges[0][-1].item(), edges[1][0].item(), edges[1][-1].item())) \n",
    "plt.title(r\"Maximum radiance $\\mathregular{W/sr\\, m^2)}$\")\n",
    "plt.xlabel(\"X (m)\")\n",
    "plt.ylabel(\"Y (m)\")\n",
    "\n",
    "plt.colorbar()\n",
    "plt.show()\n",
    "# the expected number is like 2.5MW/sr/m2"
   ]
  }
 ],
 "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}