{ "cells": [ { "cell_type": "markdown", "id": "6a3135bf-2319-4261-b1e2-34b77b17dbc5", "metadata": {}, "source": [ "# Various recipes (e.g. resolution reduction, etc.)" ] }, { "cell_type": "code", "execution_count": 1, "id": "04a6560c-00b8-4411-aa6a-10becc89e839", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ ":241: RuntimeWarning: gvar._svec_smat.smat size changed, may indicate binary incompatibility. Expected 248 from C header, got 464 from PyObject\n" ] } ], "source": [ "import orcs.core\n", "import matplotlib.pyplot as pl\n", "import numpy as np\n", "import scipy.fft" ] }, { "cell_type": "code", "execution_count": 2, "id": "acbd4be6-3051-4780-ab3a-79f7b4c70cdc", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "\u001b[0mdev.0624|INFO| CFHT version\u001b[0m\n", "\u001b[0mdev.0624|INFO| Cube is level 2.5\u001b[0m\n", "\u001b[0mdev.0624|INFO| shape: (2048, 2064, 847)\u001b[0m\n", "\u001b[0mdev.0624|INFO| wavenumber calibration: True\u001b[0m\n", "\u001b[0mdev.0624|INFO| flux calibration: True\u001b[0m\n", "\u001b[0mdev.0624|INFO| wcs calibration: True\u001b[0m\n" ] } ], "source": [ "cube = orcs.core.SpectralCube('data/M1_2022_SN3.merged.cm1.hdf5')" ] }, { "cell_type": "markdown", "id": "f8896551-a4dd-4341-b0bf-51ba2d6cd6c6", "metadata": {}, "source": [ "## reduce the resolution of a spectrum" ] }, { "cell_type": "code", "execution_count": 3, "id": "6f693c23-75e6-4ae6-bfa6-e499898f53e7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGsCAYAAADQat0+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAACaZklEQVR4nOzdd3wUZf7A8c9sTQIklBBCCb33XgWxgIL9RLAclrP+9M6CnoodPUU9C/Z2KmeFU1QsqGABpCjSey+hJJRAerJt5vfH7M7OZDdNCSSb7/v14sWW2c082+Y73+f7PI+iaZqGEEIIIUQMsZ3sHRBCCCGEON4kwBFCCCFEzJEARwghhBAxRwIcIYQQQsQcCXCEEEIIEXMkwBFCCCFEzJEARwghhBAxRwIcIYQQQsQcCXCEEEIIEXMkwBFCCCFEzKn1Ac7ChQs577zzaNasGYqi8MUXX1SLv7dp0ybOP/98kpKSqFevHoMHDyY9Pb1K900IIYSIFbU+wCkoKKBXr168/PLL1ebv7dixg1NOOYXOnTszf/581qxZw4MPPkhcXNwJ2UchhBCiplNksc0wRVH4/PPPufDCC43bvF4vDzzwAB9++CHZ2dl0796dp556ipEjR1bJ3wO49NJLcTqdvP/++3/6bwghhBC1Ua3P4JTnmmuuYfHixcyYMYO1a9dyySWXcPbZZ7Nt27Yq+XuqqvLNN9/QsWNHzjrrLFJSUhg0aFCVd50JIYQQsUQCnDLs2LGDjz/+mE8++YThw4fTrl077rrrLk455RTefffdKvmbhw4dIj8/nyeffJKzzz6buXPnctFFF/GXv/yFBQsWVMnfFEIIIWKN42TvQHW2cuVKNE2jY8eOlts9Hg+NGjUCYPfu3bRp06bM57nlllsqXOOjqioAF1xwAXfccQcAvXv3ZsmSJbz++uuceuqplW2GEEIIUetIgFMGVVWx2+2sWLECu91uua9u3boANG/enE2bNpX5PA0aNKjw30xOTsbhcNC1a1fL7V26dGHRokUVfh4hhBCiNpMApwx9+vQhEAhw6NAhhg8fHnUbp9NJ586dj9vfdLlcDBgwgC1btlhu37p1K61atTpuf0cIIYSIZbU+wMnPz2f79u3G9V27drF69WoaNmxIx44dueKKK7jyyit59tln6dOnD0eOHOGnn36iR48ejB079rj+vZYtWwLwz3/+kwkTJjBixAhOO+00vvvuO7766ivmz5//p9srhBBC1Aa1fpj4/PnzOe200yJuv+qqq5g+fTo+n49//etfvPfee+zfv59GjRoxZMgQpkyZQo8ePY773wt55513mDp1Kvv27aNTp05MmTKFCy64oNJ/TwghhKiNan2AI4QQQojYI8PEhRBCCBFzJMARQgghRMyplUXGqqpy4MAB6tWrh6IoJ3t3hBBCCFEBmqaRl5dHs2bNsNnKztHUygDnwIEDpKWlnezdEEIIIcQfsHfvXlq0aFHmNrUywKlXrx6gv0CJiYkneW+qls/nY+7cuYwePRqn03myd+eEk/ZL+6X90n5pf+y0Pzc3l7S0NOM4XpZaGeCEuqUSExNrRYCTkJBAYmJizHzAK0PaL+2X9kv7pf2x1/6KlJdIkbEQQgghYo4EOEIIIYSIORLgCCGEECLmSIAjhBBCiJgjAY4QQgghYo4EOEIIIYSIORLgCCGEECLmSIAjhBBCiJgjAY4QQgghYo4EOEIIIYSIOVUa4CxcuJDzzjuPZs2aoSgKX3zxRZnbz58/H0VRIv5t3rzZst2sWbPo2rUrbrebrl278vnnn1dhK4QQQghR01RpgFNQUECvXr14+eWXK/W4LVu2kJGRYfzr0KGDcd/SpUuZMGECEydOZM2aNUycOJHx48fz22+/He/dF0IIIUQNVaWLbY4ZM4YxY8ZU+nEpKSnUr18/6n3Tpk1j1KhRTJ48GYDJkyezYMECpk2bxscff/xndlcIIWo0VdX479Ld9GvVgJ4t6p/s3RHipKqWq4n36dOH4uJiunbtygMPPMBpp51m3Ld06VLuuOMOy/ZnnXUW06ZNK/X5PB4PHo/HuJ6bmwvoK636fL7ju/PVTKh9sd7O0kj7pf3m/2Pd56sOMOWrjQBse2x0rWt/SdL+2Gt/ZdpSrQKcpk2b8uabb9KvXz88Hg/vv/8+Z5xxBvPnz2fEiBEAZGZm0qRJE8vjmjRpQmZmZqnPO3XqVKZMmRJx+9y5c0lISDi+jaim5s2bd7J34aSS9kv7a4Pv9tgIVR7MmTPHuL22tL800v7YaX9hYWGFt61WAU6nTp3o1KmTcX3IkCHs3buXZ555xghwABRFsTxO07SI28wmT57MpEmTjOu5ubmkpaUxevRoEhMTj2MLqh+fz8e8efMYNWoUTqfzZO/OCSftl/bXpvav/34rPx3YDcDYsWNrXftLkvbHXvtDPTAVUa0CnGgGDx7MBx98YFxPTU2NyNYcOnQoIqtj5na7cbvdEbc7nc6YedPLU5vaGo20X9pfG9pvt9uNy+b21pb2l0baHzvtr0w7qv08OKtWraJp06bG9SFDhkSk2+bOncvQoUNP9K4JIUS1UkYiW4hap0ozOPn5+Wzfvt24vmvXLlavXk3Dhg1p2bIlkydPZv/+/bz33nuAPkKqdevWdOvWDa/XywcffMCsWbOYNWuW8Ry33XYbI0aM4KmnnuKCCy5g9uzZ/PDDDyxatKgqmyKEENWexDdChFVpgLN8+XLLCKhQHcxVV13F9OnTycjIID093bjf6/Vy1113sX//fuLj4+nWrRvffPMNY8eONbYZOnQoM2bM4IEHHuDBBx+kXbt2zJw5k0GDBlVlU4QQotqTDI4QYVUa4IwcORJN00q9f/r06Zbrd999N3fffXe5zztu3DjGjRv3Z3dPCCFiiiI5HCEM1b4GRwghRMVIBkeIMAlwhBAiRkh8I0SYBDhCCBErJIUjhEECHCGEiBES3ggRJgGOEELECEngCBEmAY4QQsQIGUUlRJgEOEIIESMkgyNEmAQ4QggRIyS+ESJMAhwhhIgRksERIkwCHCGEiBGKRDhCGCTAEUIIIUTMkQBHCCFihCRwhAiTAEcIIWKEDBMXIkwCHCGEiBGSwREiTAIcIYQQQsQcCXCEECJG2CSDI4RBAhwhhIgRUoMjRJgEOEIIESOkBkeIMAlwhBBCCBFzJMARQogYITMZCxEmAY4QQsQICW+ECJMARwghYoQkcIQIkwBHCCFihMQ3QoRJgCOEEDFCanCECJMARwghYoTEN0KESYAjhBAxwhzfaJp20vZDiOpAAhwhhIgVphSOxDeitpMARwghYpAqEY6o5STAEUKIGGHpojppeyFE9SABjhBCxAhzkbFkcERtJwGOEELECPNq4hLfiNpOAhwhhIgRMkxciDAJcIQQIkbYpItKCIMEOEIIESOki0qIMAlwhBAiVkgGRwiDBDhCCBEjZJi4EGFVGuAsXLiQ8847j2bNmqEoCl988UWZ23/22WeMGjWKxo0bk5iYyJAhQ/j+++8t20yfPh1FUSL+FRcXV2FLhBCi+jMvtqmpJ3FHhKgGqjTAKSgooFevXrz88ssV2n7hwoWMGjWKOXPmsGLFCk477TTOO+88Vq1aZdkuMTGRjIwMy7+4uLiqaIIQQtQY1gyO5HBE7eaoyicfM2YMY8aMqfD206ZNs1x/4oknmD17Nl999RV9+vQxblcUhdTU1OO1m0IIEROsE/2dvP0Qojqo0gDnz1JVlby8PBo2bGi5PT8/n1atWhEIBOjduzePPfaYJQAqyePx4PF4jOu5ubkA+Hw+fD5f1ex8NRFqX6y3szTSfmm/+f9YpwYCxmWvz4cvGPHUlvaXVNve/5Jisf2VaYuiaSem1F5RFD7//HMuvPDCCj/m3//+N08++SSbNm0iJSUFgF9//ZXt27fTo0cPcnNzeeGFF5gzZw5r1qyhQ4cOUZ/nkUceYcqUKRG3f/TRRyQkJPyh9gghRHXz+2GFD7bbAXisn59E10neISGOs8LCQi6//HJycnJITEwsc9tqG+B8/PHHXHfddcyePZszzzyz1O1UVaVv376MGDGCF198Meo20TI4aWlpHDlypNwXqKbz+XzMmzePUaNG4XQ6T/bunHDSfml/bWr/F6sP8M9Z6wFYcvep1I+z1ar2l1Tb3v+SYrH9ubm5JCcnVyjAqZZdVDNnzuTaa6/lk08+KTO4AbDZbAwYMIBt27aVuo3b7cbtdkfc7nQ6Y+ZNL09tams00n5pf21ov91uNy47HA6cTv16bWl/aaT9sdP+yrSj2s2D8/HHH3P11Vfz0Ucfcc4555S7vaZprF69mqZNm56AvRNCiOrLnI+XImNR21VpBic/P5/t27cb13ft2sXq1atp2LAhLVu2ZPLkyezfv5/33nsP0IObK6+8khdeeIHBgweTmZkJQHx8PElJSQBMmTKFwYMH06FDB3Jzc3nxxRdZvXo1r7zySlU2RQghahQZJi5quyrN4Cxfvpw+ffoYI5wmTZpEnz59eOihhwDIyMggPT3d2P6NN97A7/dzyy230LRpU+PfbbfdZmyTnZ3NDTfcQJcuXRg9ejT79+9n4cKFDBw4sCqbIoQQ1Z45pJEMjqjtqjSDM3LkSMqqYZ4+fbrl+vz588t9zueff57nn3/+T+6ZEELEthM0fkSIaqva1eAIIYT48yS+EbWdBDhCCBEjzFkbCXBEbScBjhBCxCBVIhxRy0mAI4QQMUIr5bIQtZEEOEIIEYMkgyNqOwlwhBAiVphiGolvRG0nAY4QQsQkiXBE7SYBjhBCxAjz7MUy0Z+o7STAEUKIGCRdVKK2kwBHCCFihHWxTYlwRO0mAY4QQsQgiW9EbScBjhBCxAjrYpsS4YjaTQIcIYQQQsQcCXCEECJGSA2OEGES4AghRAyS+EbUdhLgCCFEDJIMjqjtJMARQogYYZ7oT8IbUdtJgCOEEDFIEjiitpMARwghYoRmWWxTIhxRu0mAI4QQMUjCG1HbSYAjhBAxwjLRn6y2KWo5CXCEECIGSXgjajsJcIQQIlaY6m5kmLio7STAEUKIWCTxjajlJMARQogYYV1s86TthhDVggQ4QggRgzRJ4YhaTgIcIYSIEdbFNk/efghRHUiAI4QQMUgm+hO1nQQ4QggRI8xBjYQ3oraTAEcIIWKQZHBEbScBjhBCxCCJb0RtJwGOEELECBkmLkSYBDhCCBGDpItK1HYS4AghRIyQYeJChEmAI4QQMUkiHFG7SYAjhBAxQmpwhAiTAEcIIWKQlOCI2q5KA5yFCxdy3nnn0axZMxRF4Ysvvij3MQsWLKBfv37ExcXRtm1bXn/99YhtZs2aRdeuXXG73XTt2pXPP/+8CvZeCCFqFutEfxLhiNqtSgOcgoICevXqxcsvv1yh7Xft2sXYsWMZPnw4q1at4r777uPWW29l1qxZxjZLly5lwoQJTJw4kTVr1jBx4kTGjx/Pb7/9VlXNEEKIGke6qERt56jKJx8zZgxjxoyp8Pavv/46LVu2ZNq0aQB06dKF5cuX88wzz3DxxRcDMG3aNEaNGsXkyZMBmDx5MgsWLGDatGl8/PHHx70NQghRE8kwcVHbVWmAU1lLly5l9OjRltvOOuss3n77bXw+H06nk6VLl3LHHXdEbBMKiqLxeDx4PB7jem5uLgA+nw+fz3f8GlANhdoX6+0sjbRf2m/+P9YFAgHjst8fqHXtL0naH3vtr0xbqlWAk5mZSZMmTSy3NWnSBL/fz5EjR2jatGmp22RmZpb6vFOnTmXKlCkRt8+dO5eEhITjs/PV3Lx58072LpxU0n5pf22w8YAC2AFYvXo1jv16Fqe2tL800v7YaX9hYWGFt61WAQ6AoiiW66E0q/n2aNuUvM1s8uTJTJo0ybiem5tLWloao0ePJjEx8XjsdrXl8/mYN28eo0aNwul0nuzdOeGk/dL+2tT+jMW7mb1nKwC9evdmVNfkWtX+kmrb+19SLLY/1ANTEdUqwElNTY3IxBw6dAiHw0GjRo3K3KZkVsfM7Xbjdrsjbnc6nTHzppenNrU1Gmm/tL82tN9mC48bsdttRptrS/tLI+2PnfZXph3Vah6cIUOGRKTS5s6dS//+/Y1GlbbN0KFDT9h+CiGEEKJ6q9IMTn5+Ptu3bzeu79q1i9WrV9OwYUNatmzJ5MmT2b9/P++99x4AN910Ey+//DKTJk3i+uuvZ+nSpbz99tuW0VG33XYbI0aM4KmnnuKCCy5g9uzZ/PDDDyxatKgqmyKEEDWKDKIStV2VZnCWL19Onz596NOnDwCTJk2iT58+PPTQQwBkZGSQnp5ubN+mTRvmzJnD/Pnz6d27N4899hgvvviiMUQcYOjQocyYMYN3332Xnj17Mn36dGbOnMmgQYOqsilCCFHtSVAjRFiVZnBGjhxZ5lwM06dPj7jt1FNPZeXKlWU+77hx4xg3btyf3T0hhIhZEuyI2q5a1eAIIYT447RSLgtRG0mAI4QQQoiYIwGOEELECHO3lCzVIGo7CXCEECIGSXgjajsJcIQQIkZoEtYIYZAARwghYpHEOqKWkwBHCCFihJTdCBEmAY4QQsQg6a4StZ0EOEIIEYMkmyNqOwlwhBBCCBFzJMARQogYYZ77RhI4oraTAEcIIWKQdFGJ2k4CHCGEEELEHAlwhBAiRliWapBOKlHLSYAjhBBCiJgjAY4QQsQIc85GanBEbScBjhBCxCCJb0RtJwGOEELECMnaCBEmAY4QQsQiiXZELScBjhBCxAjzyCkJb0RtJwGOEEIIIWKOBDhCCBEjLPPgSApH1HIS4AghhBAi5kiAI4QQMcI6D46kcETtJgGOEELEIAlvRG0nAY4QQsQKydoIYZAARwghYpDEOqK2kwBHCCFihFbKZSFqIwlwhBBCCBFzJMARQogYJKOoRG0nAY4QQsQIiWmECJMARwghhBAxRwIcIYSIEZbFNiWbI2o5CXCEEEIIEXMkwBFCiBhhWWxTBoqLWk4CHCGEiEHSRSVquxMS4Lz66qu0adOGuLg4+vXrxy+//FLqtldffTWKokT869atm7HN9OnTo25TXFx8IpojhBDVksQ0QoRVeYAzc+ZMbr/9du6//35WrVrF8OHDGTNmDOnp6VG3f+GFF8jIyDD+7d27l4YNG3LJJZdYtktMTLRsl5GRQVxcXFU3RwghagQJdkRtV+UBznPPPce1117LddddR5cuXZg2bRppaWm89tprUbdPSkoiNTXV+Ld8+XKOHTvGNddcY9lOURTLdqmpqVXdFCGEqNYsNTgS4YhazlGVT+71elmxYgX33nuv5fbRo0ezZMmSCj3H22+/zZlnnkmrVq0st+fn59OqVSsCgQC9e/fmscceo0+fPlGfw+Px4PF4jOu5ubkA+Hw+fD5fZZpU44TaF+vtLI20X9pv/j/WqWrAuBwIBGpd+0uS9sde+yvTlioNcI4cOUIgEKBJkyaW25s0aUJmZma5j8/IyODbb7/lo48+stzeuXNnpk+fTo8ePcjNzeWFF15g2LBhrFmzhg4dOkQ8z9SpU5kyZUrE7XPnziUhIaGSraqZ5s2bd7J34aSS9kv7a4Mde2yEEvObt2xmXv4moPa0vzTS/thpf2FhYYW3rdIAJ0RRFMt1TdMibotm+vTp1K9fnwsvvNBy++DBgxk8eLBxfdiwYfTt25eXXnqJF198MeJ5Jk+ezKRJk4zrubm5pKWlMXr0aBITEyvZmprF5/Mxb948Ro0ahdPpPNm7c8JJ+6X9tan967/fyo8HdgPQuVNnRg1pUavaX1Jte/9LisX2h3pgKqJKA5zk5GTsdntEtubQoUMRWZ2SNE3jnXfeYeLEibhcrjK3tdlsDBgwgG3btkW93+1243a7I253Op0x86aXpza1NRppv7S/NrTfZguXVSo2m9Hm2tL+0kj7Y6f9lWlHlRYZu1wu+vXrF5EemzdvHkOHDi3zsQsWLGD79u1ce+215f4dTdNYvXo1TZs2/VP7K4QQQojYUOVdVJMmTWLixIn079+fIUOG8Oabb5Kens5NN90E6N1H+/fv57333rM87u2332bQoEF079494jmnTJnC4MGD6dChA7m5ubz44ousXr2aV155paqbI4QQQogaoMoDnAkTJpCVlcWjjz5KRkYG3bt3Z86cOcaoqIyMjIg5cXJycpg1axYvvPBC1OfMzs7mhhtuIDMzk6SkJPr06cPChQsZOHBgVTdHCCGqLfPIcE3GiYta7oQUGd98883cfPPNUe+bPn16xG1JSUllVko///zzPP/888dr94QQIuZIfCNqO1mLSgghYoRkbYQIkwBHCCFikIQ6oraTAEcIIWKEJHCECJMARwghYpAEO6K2kwBHCCFihGUUlXRSiVpOAhwhhBBCxBwJcIQQIkaYu6Wki0rUdhLgCCFEDJL4RtR2EuAIIUSMkLobIcIkwBFCiFgkfVSilpMARwghYoTENEKESYAjhBAxSGIdUdtJgCOEEDFIsjmitpMARwghYoidAC2Uwyd7N4Q46STAEUKIGGHT/Hzjuo9F7ttIzVt3sndHiJNKAhwhhIgRaXlr6GzbC0D9ovSTvDdCnFwS4AghRIxwqsXGZZumnsQ9EeLkkwBHCCFihEI4qLFp/pO4J0KcfBLgCCFEjFBMWRtFC5zEPRHi5JMARwghYoTNFNRIBkfUdhLgCCFEjDBnbRQJcEQtJwGOEELECIVwgGOXAEfUchLgCCFEjLB2UUkNjqjdJMARQogYYS4ytqmSwRG1mwQ4QggRI8xz39iQDI6o3STAEUKIGGEpMpYMjqjlJMARQogYYc7aSA2OqO0kwBFCiBgh8+AIESYBjhBCxAhLkbFkcEQtJwGOEELECGuAIxkcUbtJgCOEEDHCRjiokQyOqO0kwBFCiBhhGSau+U7inghx8kmAI4QQMeKP1ODke/zkFEkwJGKPBDhCiJMiK9/D24t2cbTAe7J3JWZUdqkGTdPo/vD39Joyl2KfdGmJ2CIBjhDipLjx/RU89vVGbv5wxcnelZihVHKYuF/VjMvpRwtL3e7f32/m/z5YgWraXojqznGyd0AIUTst33MMgF93Hj3JexI7bIS7qOwVyuCEL5eVwXnl5x0ArEg/xoDWDf/4DgpxAkkGRwhRLamqxraDeZI1qARz1kapQAZHI/zaevxqGVvqfIHytxGiujghAc6rr75KmzZtiIuLo1+/fvzyyy+lbjt//nwURYn4t3nzZst2s2bNomvXrrjdbrp27crnn39e1c0QQpxAT363mVHPL+SZuVtO9q7UGIolg1OBAMcUO3p85QcvNkX5Q/slxMlQ5QHOzJkzuf3227n//vtZtWoVw4cPZ8yYMaSnp5f5uC1btpCRkWH869Chg3Hf0qVLmTBhAhMnTmTNmjVMnDiR8ePH89tvv1V1c4QQJ8ibC3cC8Or8HSd5T2qOyo6isgQ4/ujba6aNJMARNUmVBzjPPfcc1157Lddddx1dunRh2rRppKWl8dprr5X5uJSUFFJTU41/drvduG/atGmMGjWKyZMn07lzZyZPnswZZ5zBtGnTqrg1QghRfdkrWWRs7qIqLiWDY+4htEl8I2qQKi0y9nq9rFixgnvvvddy++jRo1myZEmZj+3Tpw/FxcV07dqVBx54gNNOO824b+nSpdxxxx2W7c8666xSAxyPx4PH4zGu5+bmAuDz+fD5Ynv+h1D7Yr2dpZH214z2l7d/f3T/a0r7j5sSw8TLa7/XGw6CCj3eqNuZa3PUQKBGvZa17v0vIRbbX5m2VGmAc+TIEQKBAE2aNLHc3qRJEzIzM6M+pmnTprz55pv069cPj8fD+++/zxlnnMH8+fMZMWIEAJmZmZV6zqlTpzJlypSI2+fOnUtCQsIfaVqNM2/evJO9CyeVtL86tj/88zNnzpw/cH/FVc/2H3/1CvKMy35PkdHu0tpf7IfQ67x81RpcB1ZHbOMNhLdZunQJGeuP4w6fILXl/S9NLLW/sLD06QxKOiHDxJUS/baapkXcFtKpUyc6depkXB8yZAh79+7lmWeeMQKcyj7n5MmTmTRpknE9NzeXtLQ0Ro8eTWJiYqXbU5P4fD7mzZvHqFGjcDqdJ3t3Tjhpf/Vt/21L5xqXx44dW+n7K6I6t78qrN75IeTol+OcNkaNGlVm+/OKfdzz+88AdOjclbFDWkXZxg/LfgJg2LBh9GqRVHUNOM5q2/tfUiy2P9QDUxFVGuAkJydjt9sjMiuHDh2KyMCUZfDgwXzwwQfG9dTU1Eo9p9vtxu12R9zudDpj5k0vT21qazTS/urd/vL27c/ue3Vv//FiKzGKKtTm0tpvN5Xp+DUl+jamHgGX01EjX8fa8v6XJpbaX5l2VGmRscvlol+/fhHpsXnz5jF06NAKP8+qVato2rSpcX3IkCERzzl37txKPacQQsQay1INVGyphpDShon7VdPILBlFJWqQKu+imjRpEhMnTqR///4MGTKEN998k/T0dG666SZA7z7av38/7733HqCPkGrdujXdunXD6/XywQcfMGvWLGbNmmU852233caIESN46qmnuOCCC5g9ezY//PADixYtqurmCCFEtaVUei2q8OXiUoaJB0zDqDSZc1HUIFUe4EyYMIGsrCweffRRMjIy6N69O3PmzKFVK72vNyMjwzInjtfr5a677mL//v3Ex8fTrVs3vvnmG0sf/NChQ5kxYwYPPPAADz74IO3atWPmzJkMGjSoqpsjhBDVlrmLqkIBjulyaRmcgCmqUSXCETXICSkyvvnmm7n55puj3jd9+nTL9bvvvpu777673OccN24c48aNOx67J4QQMcEc1JhnNS6NuYuqtAyOP2DK4PyJfRPiRJO1qIQQIkbYLDMZlx/gmCfxCwSihy+qZHBEDSUBjhCiWpJZcytPqXQXlRb1spnfUoMjAY6oOSTAEUJUS3aJcCrNXskuKnNMU9qi7aoUGYsaSgIcIUS1VNrEnaJ0CuZRVBWowTFfLiV4MWdwSguChKiOJMARQlRLdglwKs1Sg1OBeXDMNTWldT8FVKnBETWTBDhCiGpJeqgqzxzUKBXJ4JjildJCFwlwRE0lAY4QolqySYRTaeagxl6RmYxNl0sLXsxdVDJOXNQkEuAIIaolKTKuvIiRU+VkcSpSQGwdJv6Hd02IE04CHCFEtSQ1OJUXUXejlp/FMTYtLYMTkC4qUTNJgCOEqJZkFFXlRYycUv3RNwyqSA2OpRD5D+6XECeDBDhCiGrJLr9OlRaRwSlnsj/LRH8VqMGRDI6oSeQnRAhRLdkkg1NpkRmccmpwzBmcUmKXgOk5ZCZjUZNIgCOEqJYkwKm8iNmLy8vgVGCdqYDpKSW+ETWJBDhCiGrJJr9OlWYvGdCUU2RckZmMzRkcGUUlahL5CRFCVEsyiqrylIhRVOUVGZc/BNycwZEaHFGTSIAjhKiWZKK/youowSlnHhxrvFJakbG5BucP7pgQJ4EEOEKIaklqcCqv8qOowkpdTbwC61UJUR1JgCOEqJaki6ryIkdRlVdkbL5ckYn+/vCuCXHCSYAjhKiWJL6pPFslR1FVZBkG82Kbmkz1J2oQCXCEENWSeS0q6RqpmIi1qMqZB6dCq4nLWlSihpIARwhRLZlrcAJyZK2QyLWoyhlFVYGZjC0ZHAk0RQ0iAY4Qoloyj6LyS4BTIZGjqCpTgxN9m0AFVhwXojqSAEcIUS2ZR4lLgFMx9mAGx6fZ9RsqUWRc+kzGshaVqJkkwBFCVEvmUVSBgBxYy2Wqt/HhAECp1GKb0bexBjh/Yv+EOMEkwBFCVEvmUVSSOagAUzDjDQY4lZnor7TXWFYTFzWVBDhCiGpPDqwVYCoo9uKMuC3qQ7TyMziqWoGhVkJUQxLgCCGqJfMBNyABTvnUKBmcyiy2WUr0olYgyyNEdSQBjhCiWqrIStfCxJzB0UJdVJUpMi5lG6QGR9RMEuAIIaolc7ZA5sGpAC2yyLi8if6owDw41h4qeR9EzSEBjhCiWqpIAawwMXVHGQFOuUs1RL9spslMxqKGcpzsHRBCiGgsK12Xl4gQRhdVQFPwh85dy5vJWIOJ9rkkUcAv2tXRn1ZWExc1lAQ4QohqyZo5kANruYLZGj92VCPAKacGJ+DjMed0ADyeTsApkdtUYLZjIaoj6aISQlRLMoqqkoLBjIqNQOinvZx5cPB7jYuDvUujP610FYoaSgIcIUS1JF0jlRTqosIWzuCUU4OjqB7jskPzRd1GanBETSUBjhCiWtI0UFCxoRKQGpzyBbM1AWwEtIp1USmmDI5Di16vYx2uLxGOqDmkBkcIUS1paoAvXA/hxoca+OVk7071FwxmAtjwE1xss9wMTjjAcWreqNuospq4qKFOSAbn1VdfpU2bNsTFxdGvXz9++aX0H6vPPvuMUaNG0bhxYxITExkyZAjff/+9ZZvp06ejKErEv+Li4qpuihDiBGmiHqKXbSedbXtRCg6d7N2p/owuqooXGRMIBzUuSglwpAZH1FBVHuDMnDmT22+/nfvvv59Vq1YxfPhwxowZQ3p6etTtFy5cyKhRo5gzZw4rVqzgtNNO47zzzmPVqlWW7RITE8nIyLD8i4uLq+rmCCGOk/bKPl50vkQ7ZX/U+9MC4d8ImyfvRO1WzaWFMziBPxLgaJ6om1hqof7cHgpxQlV5F9Vzzz3Htddey3XXXQfAtGnT+P7773nttdeYOnVqxPbTpk2zXH/iiSeYPXs2X331FX369DFuVxSF1NTUKt13IUTVecP5PO1sGQywbQFuiLi/VWCPcVnx5JzAPauhTF1UFS4yNgU47lK6qCx/QjI4ogap0gDH6/WyYsUK7r33Xsvto0ePZsmSJRV6DlVVycvLo2HDhpbb8/PzadWqFYFAgN69e/PYY49ZAiAzj8eDxxM+O8nNzQXA5/Ph80UfORArQu2L9XaWRtpffdvfzpYBQFPlaNT9a6hmGZfVomN/qA3Vuf3Hm+Lz4AACWjiDo/r0oKW09qveQuOyU/NG3c4fCAdJfn+gRr2Wten9jyYW21+ZtlRpgHPkyBECgQBNmjSx3N6kSRMyMzMr9BzPPvssBQUFjB8/3ritc+fOTJ8+nR49epCbm8sLL7zAsGHDWLNmDR06dIh4jqlTpzJlypSI2+fOnUtCQkIlW1UzzZs372Tvwkkl7a9+7R+tuYhX9APwnDlzIu63+wqMyzvXL2fzUecf/lvVsf3HW8P8rQwn1EWlALB500Zo3KzU9ucf2MiA4GWnVhz1fdi120aommHLli3MKdhcBXtftWrD+1+WWGp/YWFh+RsFnZBRVIqiWK5rmhZxWzQff/wxjzzyCLNnzyYlJcW4ffDgwQwePNi4PmzYMPr27ctLL73Eiy++GPE8kydPZtKkScb13Nxc0tLSGD16NImJiX+kSTWGz+dj3rx5jBo1Cqfzjx8gaippf/Vt/5GVdYnnKABjx4yBEr8JC9a+BcHkQceWTWg5emyl/0Z1bv/xpuxJgm2hif70UVSdO7Zn1zFKbf/aBcVwUL8ch5+xYyNf41+/3Mjig/sAaN+hI2NPa1d1jTjOatP7H00stj/UA1MRVRrgJCcnY7fbI7I1hw4disjqlDRz5kyuvfZaPvnkE84888wyt7XZbAwYMIBt27ZFvd/tduN2uyNudzqdMfOml6c2tTUaaX/1a3+eFk/TYEzj9OdDgrUbOo5wt7LTn/+n9r86tv+4s+kvph+70UW171ghO3NLb7+D8Nw3LrzRXyPFZrpoq5GvY614/8sQS+2vTDuqdBSVy+WiX79+EemxefPmMXTo0FIf9/HHH3P11Vfz0Ucfcc4555T7dzRNY/Xq1TRt2vRP77MQ4sSIMw9L9hVF3G8etmz3VPysrdbSwks1hIqMZyzbwwsbHARKmYLYUmRcyjBxTdNooRymi7JH5sERNUqVd1FNmjSJiRMn0r9/f4YMGcKbb75Jeno6N910E6B3H+3fv5/33nsP0IObK6+8khdeeIHBgwcb2Z/4+HiSkpIAmDJlCoMHD6ZDhw7k5uby4osvsnr1al555ZWqbo4Q4jipp5iCGn/kHFZu07Blm1cCnHKZRlGFZjK2Bwd2lx7gmJZqIAABP9ithwVN1Vjkvg2A1z29gY7HeceFqBpVHuBMmDCBrKwsHn30UTIyMujevTtz5syhVatWAGRkZFjmxHnjjTfw+/3ccsst3HLLLcbtV111FdOnTwcgOzubG264gczMTJKSkujTpw8LFy5k4MCBVd0cIcRxoKkq9TAVC0bJ4MSZAxx/5P2iBDW0mnh4FJUNffmG0pZYUNQSI1L8xWCva7kpyXvAuFyvOPqcRUJURyekyPjmm2/m5ptvjnpfKGgJmT9/frnP9/zzz/P8888fhz0TQpwMqrcAhxJeYErzFVFy2IHbVINjzjSIUmiRq4k7glXapa3Gbu6iAvQAx20NcJoVbTUuu3wy4aKoOWSxTSHECad68q3Xy83gyDIs5Qou1WAuMrYFg8jSVgG3lQxworwPKZ5whj3Bf+w47KgQJ4YEOEKIE071WQMWzRs5t4W56FUJSIBTLjWyyNge7KJSS6vBUUtmcCIzZW41PB9RnE8CHFFzSIAjhDjxSgY40TI4mDM40kVVLk0PZswzGYcCnNK6qCIyOCWvA65A+L1J8B49HnsqxAkhAY4Q4oTTSgQsWomAB00j3hTgZGYdK7VQVgQZq4lHFhmX1kWlBEoUGUepdXKppgBHMjiiBpEARwhxwpXM2ERkcErU3Lg0D6v2ZlfxXtVw5mHiwZmMK99FFSWDYwpw3IH8iPuFqK4kwBFCnHARGZuSAU6J627FR7Gv7JWxaz3NHOBUtIuqRMYmSheV2xTgOFSphRI1hwQ4QogTLqKLquQoqRLX40qZZVeYGF1U9ogi49J692wl58Epp4vKKcP1RQ0iAY4Q4sQrGdB4y+6iisOLEjFTjrCIOpNxMINTSheVTS0RsETpojJncJyqTLgoag4JcIQQJ1zEMPGIGhzrgVYyOBUQHEWlRikyLnWiv4gMTnkBjnRRiZpDAhwhqpA/oLJw62Hyin3lb1yLlKzB0Xwl5sEJZnAKNTcADkXFpslrWBa/X399/FHmwSltBJotYhRVlABHkwBH1EwS4AhRhd5YuJMr31nG36b/frJ3pVrRfOUMEw8eaHOoY9wksxmX7T8LtwF6F5U/FOAowW4rNfpjIruook30Fw5wXBLgiBpEAhwhqtCM3/Vp7n/fLfOHPDd3C8Of/omsfE9EUXFpw8TztHhULVh7UzIIEhY5+frrE9Aii4xLGyYeWWRcIoOjqsRp4dddMjiiJpEAR4g/yR9Qmb16Pxk5kQWYdkUKY0Ne/Gk7e48W8c7iXZFFxhEBjp5J8ODEgxOIvl6VCDPqbaIME1fLmcnYp+nz5kQEOCVWcXdqXlBLSQcJUc1IgCPEn/T+r3u4bcZqRj+3MOI+m00CnJI0DSPAKdJc+o0Rw8T1AMeLk2L0bSLqdP6kQ3nFvLlwB8cKYqOA2QhmKlFkHKprKiBOv6FkgOMtIIJfAk1RM0iAI8SftHDrYQDyPP6I+2ySwYngsClGAJNNXf3Gkt1PwYDHo5kDnOPXPeIPqIx4+meemLOZO/63+rg974lWXJgHm+eAt8Cot7EWGeuBTWlJl1AGJ594/YaSw8S9+szFoWJv/bbjG2gKUVUcJ3sHhKjp3A57qfdJF1Uku81mBDA5Wh2aKkdRSmYFggdeL06KNScoHNcanHtmraPYpx/15285fNye90R6b+luHN/czuWOn6HnpZYMjr/kUg2hDI6mwbpPoVlvSO6ALbhUQ74Wr7/GJSfyC2Zw8olH0TTiFS8c50yaEFVFMjhC/EluZ+lfI+miiuSwhzM4oVFSSqCUDI65i+o4do3MWrnvuD1XVTucF3324Idmb9CDG4C1M8KT+mFDDU6KaCc4iioU4OxZDJ9dBy/3B00zApzyuqgKNDeFBLM4UgslaggJcIT4k9yO0r9G9hj7hi3bdZSHZq8nP0p3XFnM87DYbQp5+XrXR64WDHAianBCGRwHnmCAUxsPrNMX72LA4z/wxoId5W4bCnD82CNmMjZGUWXvDT9g/wrswVFU+Vr0Lqo1O/YDUEQcRUaAc/wzOL7SxrEL8SfE2M+vECeeuYuq5IRqsVaDM/6Npby3dA8v/LC1Uo/z+MMHsMN5HlbvygTCGRzzHDfLdx8lNz9PfxwuI4NTUJBf6oR1ldVd2ckzztfpruykOifZHvlqIwBTv91c7rZlFRkbo8Q1UyCRtcMYJp4XqsExdVEdyffw2ry1ABTgpjhYEK55CygoJcAtbTh6WX7efIhuD3/PJ8v3lr+xEJUgAY4Qf5LLlMGZ+ftey498tADnf8v3snRH1gnZt6qy43CU0TVl8PjCB9ZNGbm40Q+sOcEMji3YRbVuXw7jXl/Kf37WD+gezWkcWL9btZvbZ67+s7sOwP85vmScfSFfux+ImSybzdJFVcowcfOoKH9RuIsqlMExdVEdyfdQh9CM0nEUBQPNd+dvotvD37PxQK7l7y/efoT+j//Aa/PLzzaZXffecrx+lX9+urZSj6tJcopkFu6TIUa+2kKcPA57OIi597N1fLoiXN9hN6UHNE1jx+F87v50LZe99SuFXv0s2OtXmb54F7uOlB80rEo/VmpNxolU2uKNpSn2ByyPjQhwghmcZbuPAuBS9Pu9OIwMTpziZfbqA39ux4PaKRnG5SRF7/rKyCliyY4jx+X5TwYnodXETfPgKCpJ5GMrDBZSe/PCD/AVYS9Rg6Oauqg0DeIVj3F/YXCbZdv0bqvXSnSbXfGf3zha4OXNhZULcOzVOYV2HHz0Wzq9pszlv0t2n+xdqXUkwBFVLtb71wMB68H+l+3hg2Sils9g20bcePH4VfKKw6n9HzcdAuCuT9bwyFcbuaecM9hV6ce46NUlnP7s/HL3yetXWbM3u9KBSFnMz1XaxHGlKfaFA5xiv4o7GMCEuqjsqgd8xQze8hSn2tYYAZDeRaVP9He8FtxUVY0kJRxMNrbpmYgzn13A5W/9VoOyaxoBLRwcNFD04CWg2QgER1E5CPCV634Gzx4Bx3ZbMzi+QuyatYtq7e6Dxt0BVQtncIgzMmnx6EGPwxyYbJ3LePvP+j5V8jPnipUUWinu+3wdAA9/ueEk70ntE9ufLHHSrd2XTY9Hvq902romCQVw3ZTdNOYY/lBAp2k8kXUbM1z/4m7HTDw+Fa+pFiX9qF6s+eUaPSsRyl6UZtE2PXAyB0mleWj2ei54ZTFvLNxV4XbsPVpo2T+zdfty6PPoXON6pTM4pi4qjy9gBDDZWnitqeyfnqPb3o/5r+spXEaA4zAOrO7jFOAUFBeTQnjpjBRFD3AKvHoQNm/jwaiPA8gu9EatMwmoGv/7fS8Hc8seyv7VmgP83wcrKPD4ycwp5rX5O8gu/GPtisOLXQnvS91gMBLAbmRwEiimpe0wNi0Ai1+wBjjF4S6mAk3PzhzOzg/f7QuQoBQH73cbRcYJwayOkXnRNPjoEp52vsUo2woSXBWffUTTNEsGtCp9vfYAz8/betzquET1JwGOqFJ3f7qWYp/KU9+VXyR5vBwr8PLFqv0UeQPlb/wnpGcVctuMVazZl0NfZSvfuO9jofsOkov36BsUZ9M0oAcvI22rKfYH2JKRQxslAxe+iELNhnVcZf49t9NGN2U3P7ruZPUrEy33ZRd6LT/cM37XCzZfWbAT0GsASv6wH8ot5tm5W8jIKeK3nVkMf/pnrpm+LOrfPu/lReSaAquyMjgbDuQwdc4mck0rqJszOF6/Gu6iMi2mqe1ZGm5r8H6vaaK/OKVygUBA1Xjsm83M3m1j15EChkz9kbcX7aIga78lMGik5FgeV1rAsSUzj96PzuOuT9dE3Pf12gPcPWstw578qcx9+sfHq/h2fSb/+WUXf/9oJU99t7lydUWefPj4MlgxnXpYR5WFupPMXVT1FNM2xTn644O04nC780kAMAJL0IPSUNCkj6IKvg8lMzgF4Yzl+fYlOB0KmqaRlV92V6ovoHLBK4vJLjwx9Sl//2gVL/y4jV93ln0iIWKHBDiiSpkPbCfK1dN/5/aZq3lizqYq/Tu3zVzF7NUHWL03m1PtevdSvOKlbWGwq+nYHmPbFsphTntqHr459/Kz+05muh6jsEQAVt55rKrB4863aWfLoPfhLyku1g8+i7cfofej83gkSgrc61eZtt5O/yd+5l/fWF+P22as5qWftvO36cv5aFl68Lmy+H330XLPcktL4BR6/Zzz4iLeWLiTr9eE61w8fpULbYt40fkSqrfQyMYUanHGOkhOT/jAkxwMOopNo6jcVO5A+MGve3jv13R+yrDx77nbyMgp5rGvN+I5aq3jSVZyLRmp7BIFoTsP53PJ60u47K1fAfhs5X42ZeRyyJStWbNX31+/qrH9UD7lOZxfzPI9ehapMhMNrpzxKGyZA1/dRpxiDSASCAc4oSLjuuYgyO8xZiYGyDwUzlTlBzM4LsJBbLEvQJ3g4/O0eIqCsxnHB987I/OSG645SyYXl93GU99tod+/fmD+lkOltmVVejZr9+WUej/oGZ7jvZTGobyTv2BoRk4RX605cFy7kEUkCXBElSo6QQHOyz9tY9RzCzhW4GXN3mwAvlp7fApSS7M1M1yw2UfZZlxu5Ase2LPTjdvcip80dR+X2Bfo29u248jfb3m+0E9dyS6Qr9Yc4Nm5WzhW6DUO/AA7NuuZhKeD2bH/Lt1DNLvy9APR24us3VVLd+q1JpsycnE7bPRUdpBAMZe8vpSv1uptOJRbzL5jkfOelDYceM3eHBzBg+QR0xl8sS/ANNernG9fyoW+OUYNjnkiv7jc8P63DRYBF+E25sGpbA2OOdAwF7L68q0BRTI55Hv8nGJbx632z8gtsGZG7pm1lt93H+Oo6UA75oVfGPjEj6bnD29vXnT1WIGXtxbujDio2hWlwsW1vZXtDLZtZOq3mziwfbVxeyLW9yWeKBkc8zbeAkuAc/SIHnz4NLvxGjuVcIAzb+NB6gYzQPnEG11UoUyRwxZsdE74c5ys5GC3KbweLEC+65PS68qOFpRfLD/z9730eWwu3/8Wfp456zK46f0V5BX/scxPdQgqTn9mAf/4eJVxYhGSle/hxveXM3dD5knas9giAU41oGkas1fvZ3cFRtGcDEt3ZPH6gh1/qO+6qruJQp6Zu5Vth/J585edxm1VvUxCvEvPPCio9LaFa4ySjQDHGnCcY//V0mWQlrPC8mN7tMDLM99v4ZSnfuL2GauM2//x8Spe+mk7CzbspRnhAtisnSsB8AZKf1/ceLnAtoiWin62bmTUvIWMtK8zMildi1fxpftBZrkewY2XWSv2EVA1Bj7xI6c89XPE84Zmxj22dQnvf/QBOw7rB87kpY+x1n09vZXtmD8uAVNQ0Vg9FLWI2OELB4xtbfoPfJHmwqOFiowrd0Az1xMlxTuNyzmHrYFvklJAXpGXD1xTmeT8lGH531vuP5LvxY2X6+1f006xBqUhWQVeGpPN/9m/xHZwvX7j7sUcevF0fvj2M/7vg5X4AirxFOvzDCsKqYlxxuO3HcyL+rwJFPOF+yFmuP7FVwuWkaqEa4c6lNiXaF1U5q44zVvI0WPhxyeo+nvmw4E3uGpPqItq95ECZi7fa3RR5WvxRhdVKFNkBGi54f1orGTjN32mj+R72Hu0kLs+WcPDs9dzz6drjc98tKkGSgbO9362jvNtSznr2+H85/H/48Ev1nPzhyv5bkMmL/ywLeLxFVEdApzQid8vW63B9v99uJLvNxzkhvdXWG7feTifa6f/zsr0Y5TlQHYRZ09byIwSgRPo668d74y6xx9g/f6calvXJAFONfDV2gxum7Gakc/Mr9D2voDK6uM8QqYsl731K09+u5nv/8BZhbm49HjIzClm5u/ppX5RzQFVVS+TEApw2ikHSFTCZ8opfv0A6sm2vl4X2JZYrjcu3BFR1Pvyz9s5kFPMF8Hh0Ob2KNnp2EwHrMN7NnL6M/PZlGGdj8T8Y/O881VecL3K1677SSXLCERY+DTTnVOZ4foXoNEtR18JvYstnZG21WiEi6CjUTXAV0Tcx39h4tZbePHNt0DT6LD9XRIUD6+7nrdk7+KOhLvPmmkHTQGOMzxTcRRFuI0Da3yJLpmtB/OYNHM1e7Kinxh4TEPTzZ+XnKwMy3bxaj5fLvzVuN7Ds9Jyv9th4x+Oz7nf+RHvOp825psBuPKdZdzw3nKy8r2863qae5wz6LbsHv3Ohf+mk2cdM92PsWbPYQp2r2CV+0Yec7wLQMsEH+faluLGy4e/6QekQ3nFlu/1MNt64/JI+xpSlXA3Xmeb9SAWCjxUbMZMxmYFBblkHA7Xy4QCHK8pwAkNNc9Z/x1tlQOmDE6c0UUVyqQZ3678cFdXfaUA/F762bbylONN6lHIPz9dw6cr9vHfpXuYGZwDKrvQy7+/3xKxjze8v4KLXl1sqd/6m2MOANf5PuKzX8O1fNsPl98VGPJHRgBqmsY/P1nDiz9WLJBauy+7zAL1aEoWWC/bFb0+6LK3fuXHzYe4zXTiE80TczaxOTOPez9bF3Hf+DeWMuDxH0qdpPGPmPS/NZz70qKITFR1IQFONfDbzsoNS53y1QYufGUx0yo5m+yftSer/CnaP16WbvmSessZIq5pGmv3ZVf4zOIvry7mnlnrePmn7VHvN/+QOao4wKkTHC0yyKb/6BYEDwBJqn6WtWid/sPoDdaYtLLpXQJ5wUnVUrzppY5aAvhl22FL91AT1frjGcjaw84oWb/QquY9lJ2MtetFw4lKIePsCzlW4GPdvhx2rNSzMn1s22lGFmn54cLZ7rbdaJrGlszoWQUInmkfWEW8ph8ATy+eCznhmWgbk83rC7Zzy0d6sODKC2ezmnM4HOCYJvKLphiXsdJ1aMhy6O+Pfn4hn63aX+oIPfNnL6fIFwxMNLbv2g1ArqYX1iYE8ti4LJylah2w/li7HTb+Zv8OgJa2w5agY+HWw8zdeJDd+w7Q3aY/b/3cLZB3EO/e8Fl4mnIY18+PEKf4uMLxI/7ifC4v+C8vu17iTedzHM4r5j+/7GTg4z9yn+ngNNAWPqAPsG0mmXAXZSvF+nkIdVH5TaOozDRvoaWbL17VPztenHiDWTIXfgL7VtFr/t/4yX2X0cVVsovqYttCblx3KWTtgCJrViFRzWaW6xEmOOZzi2N21EzN1oPh4KSFcpj66J+1HzYdZFV6Nq/P32EEqE7Cvw3tlHD2rTKT5/lNnwV/BU8MV+3N5pMV+3hunv47uz+7iPFvLOW79dFP9M5/eTHXv7eczZnWE46yEsmlzXRe8uaDufp7ezCn7G498xIqq/dmc8g0Z9bK9Gzyiv3ljtasjG+CXdlvLdz5h2axrmoS4FSRfI/fchZSlsr2pHzwq/4D/FIpB/mqUt6yA7/vPsrkz9Yx/o3waJiByiZ+dN3JJfb5UR/z0bJ0zn95MVOCU9KHMjTms2+zAzn6Qa60MyW/Gv4hq+oJxNxOPXA5w6YfxGcHhgJQV81D0zR8+foPyQq1k+VxXwWGAJDq21dmADjx7WXc8b/VxvUmivVA0kKJXpx6NF8/iF3tsHa1DLevI9/j47yXF5FQED6Id7LtpWFROADpruxC0yizSNsXUMlYv8C43kbJgMzwgd+uaKSQzTdrMyjyBnAWhotNGyj5lhqcorIyOJqbvGAgUk8ppBlH2Pnzf7n8rfBnLBTQkZsBe8OjwHw+H52UdECjfv421ruvZb5rEh0UvSg2w9Ec0IO/JqbMSHMO4TN9/hraCo2h0QBDbZHF3C2Lra9Vzu8zcPnCB7p2ygG8meFgpVn2cs726O/Pqfa1KDn7jSLwUO2Yxx+wBDFdlHTilPBvSmqJz4NbCU/0p0b5abf7i4zXHfTADvQMji+YwXErPtYvDwd7nWz6a2XuoorHy7Ou12ns2YP23b0s32yt7erhD78+rZVM6sU5SCWLrspu9Ll7NI4V6t1+rzmfZ5H7Nha476CzEv5Mvjp/B9sO5mNDpb0pqDG/HqvSs5m9OnqXodmWHIUBU8NtKu1AvD+7iPX7wwGkOYAKqBpTvtzAsl1HuemDFdEebli/3xrglNVVbj4JMwdh5rmBCr1+2ioH6KSk06VZYpl/2+zCVxYz7Gn9O2qei+zPlA28+OM2rn53WcSJ2e6sQs55aVG16P4zkwCnCqiqRq8pc+n5yNyIA/WmjFwmf7a23PkyKuJELHNk/jEo7+9Fy/D8z/0Y7WwZ3Ob4LOpj7v9cPyh+HExxXjP9d+6ZtY5n55adndKI/kXym+pRogU4qqpx5//WGEWQFaFpWtQfRY8vQCclnVNtevbj08CpgH4mnZ2bRxL6mevPai/L4z4OnA5AqpqJ1+tBQaW1kmEU55qZfyxT0Q/C+50tgVCAozFA2WzUTtzw3nL2HC0kmRzOtelBwN2+6wH9ALklI486FNHUdEAfYVtrGR7c3babg7nFZXZRbc7MY8mvi43rrZVMDu+1vmdtFP1Md/qS3Wh54bPexqZCaXORcTRFuMgNDmFurhzhO/e9tF1wK0l7vuchx3s843ydBsGMWeD9i+DtUWSumQfAmcdm8r37Xq6xf0fv/EUkKB5a2w5yml1/v4640wBIooBkJfw6xyte8rPCB9WWWLu0htg2RuxnX8XajbFt/geW691tu6jvD3cPNUifZ3QHATTPDxfRFnoDvLVwJ/3/9QNpSjgw7GyzrtVkDsrMVM0WNYPjL863ZHBCnzevZu2iUo9FrgmVT3gUVWMlO/ycOZkU5lhnf+7sCWegUpUs3HYbM12PMcd9H7fYZ7NmbzbPzd3KA44PGGP/HYAkpZC7HTMsz3PLRytppRy0BGVpyiFc+Hja8QaX2X/kthmrLY/Zfiif1xfs4ItV+43v7Ksb7RT5ys/gDHvyJ859aREHsvWspPkgXuwLWLIhJWtOzNdLTjNQz1bEX+3zqEvk98ncjZ5b7Ke9so+b7F/SwR7+vhw9ls3nroeY45rMQO9vgP479sScTXy7zvrZLK33zVwqEBq9eazAy80fruCVnyt+ovzcvK3M33KYb9dnRLwGmzJy2XWkgNcX7GD7odKzvyeSBDjHmaZp3PnJGiOSPZRrTSmOf2MpHy/byx2muS/Mn5PKFGudiBlAC01dR+VlcMy7o6oaAU/4C91COYKmll+PE6on+ei3ivfpml+zgOlytDOnhdsOM2vlPp40LV6oqhrzNh6MugRCQNW49r/LGf70zxGjNvI9fu5xzMCuaHwTGMhKrQP+YO3D9j17qa/oafhNWisOa/qZ1x41hfVaawo1Nw4C5GZs5yXnS8x338kzztcj/r4DP32VrTjxGxmcfXV6AtBMyeI6+xw+cT/Ks87XAJi78SD//N9qbnbMxq34Wam25/PAcLyanUSlkLzMHTRXrAekUXb9jPSA1pCAptBYyaHoaOSZsYJKe2WfUZjc0nQ2nagUsdF05g/QOlgo/NR3m8nKiD7Cy4OrzC6qIsIZnOZKllHrNMnxKX9zfMc4+0KGHv4fePKxH9YzINu+01+LCTnvAPCw831aeCJ/xDPj2gN6kXEjrGfdxYfCxeopfv212K02AfSuv77KVv7rfJLnna+QSD59bXqAs0NtCkB/mzXYG2mzzp1zjv03y/U2xdas0ONzNpFX7KOlUvow61SORb3dT3gmY7M4PFFHonlx4iHcRZVUuDtim3wt3ugq7GsLv5Za/kHjcx7SzRbO6LRTMnAUZBrds6fbV/HcvK0cObiPCXb98/KI70oATrWtoTHZNOYYo22/czDrGB0Va7DVSjnEGbaVjHcsYKrzbform9E0jUO5xew9WsiZzy3gyW/1uYW+Wnsg6u9ptJnVzUPR9x3TAxzzArFFvoDlhGn40z9b5vYyB025JbrO7rLN5F/Od3nJ+RIBVbPUwJh/o7ILivnY9Tj3OmcwRXstvL9bvidJKcSuaJxZ8A2gTw765sKd/N+HK9E0jVs/XsV1/11eymmftR4t9PfvnrWWOesyo9ZClScr38uxKPMXPTFnE09+u5lzX1pU6eesChWfclJUyPr9uXy+KnxwKPllCs1Cu6TEdPB1KaQIN4XeAHXcFXtbTkSAY/4yRivOO5RbTG6xj/Yp9bDbTGlVXwDtcDr1TNt6cjKJa9As6t9Jrms9yOVXohDO/ENUXgYnMycyc/bx7+nc//l6XHYbb13Vn1M7Njbum7Mug5826z/OWzLz6N+6Yfj5i49xun01AP/2TwAUcqhDI/LYvW8fQ4PLAWRrdbnV9w/G2JYxKzAcDRu7tFS6KXuY/fEb3OvUD3YX2pfwlv9ctmvN8OJAQ+E/zmcZaV/De/5RRoBztH43vMe+xaUEeMD5IQDn2X/lWf8l7NZSubP4ZSY45gPwiv8CfDjYqTWjs7KX+LydNFasZ1ctggHPZrUl+Uo8HZX93Kl8wCD3JpapnbnbdyN1KeJt1zP0sW1ntdqO8d6HaF2iBqRj8VrLZD5tTOs9lexeA1A1BR/2cjM43ig/U6GuE4AmxTvRMlYbf7pewR4osH6/RmAtHAZId+kBTiKFNFSsAY73yE5Az8glF+sH2aVqVxwOBy3U/XzmfsTYti5FRoDzUeB0HrR9aNx3VKtLQyXfMsoOMAK1Is1FvOKlsz9yIsxG5FJH8aBqCgFsOBVrNthccG5mXk3czKUEcBHZPWHuonLhI+dgesSpbz7xRqBp5iw8RCNF/07sUJvSzpZBd2W3cX89pYhmBRsIvcVpwW7VobYNuJQAG9VWTA+czUX2RfSy7eR0+yom2ufR3bab+YFerFQ76G3SFGyKRlMly1KXdIp9PYXeABe9uiQiK37bjNV8161JxD6/s2g3czcc5K0r+9MgOLnmVtMottDPan6RPuP1IRpQ5A1YgpF9x4p4bf4OBrRugNNu4x8fh4t/S86jdIVN74o8zb6Gv7+3iK83hz9r5iJj/+5fjexmP9sWfAXHcNZpgCM9PDihi38TqCobDoSzoAXegDEbevuUuhHtBet0HaHpDsyjdjVNQynnBNbchVbg8fPd+ky6Kzs5qiVygGQA47fyeA8u+aMkg3OcbT9sPXgs3n6EtfuyAUpdyK9JwWZ+d9/M0843K1Xhbl7FWlU1th7Mq9DjNU3j/z5Ywf99sKLcjJH5+aIVAk/8z2+898L97PnfPRDwWR5XfNSahZk1dwHLSxS42VAZYtuA4q1cStO82+YAx9wHHApwCnxw1bvLmb16v6UuKrTtl8ERS96AylXvLOMd03wx5hE6JZdI6OLTuyq2q83Yreln7tma/gOTdTiT+uhntjnUYanajYf817BG0w+qO4Pb3+KYbXnOb9z3sdF9Df9zPUo3ZQ8jg90pl9gX0DrY5aMkNme/lhzxmvRVtjHatpwJjvkENIUnfJfxo9oPgL1aCgAJBftoHCxUzdQaWB6/R2vC72pnAC6yL6aZcpQL7Uv4wPUEs10P0id45t7btoPL7D8ZP8abVb2rJ9TttdLRBwh3UQGkRAlw9KyBUuYoqmLNzbVn9Cr1foBmvnRy0sP1P01tRy0Fz6Af3Eva62oD6EsPhPb9WPD900yTNKb49GCqd+++rEkeG/E8o+wrqacUcVSryxeBUyz3LSjRPblBbWW5Huqu7KbsIQ4P4+wLmOz4kHbKfiN7k0FD4/NSEaUVGZdGLzIOd1ElYi0K9ms2CnEb61WZKWg0V/Rgcrum1zQ5FOvBbbgt3GWVomRTl0IjSFmqdgXgx0BfAB5x/Nco1h5pX8OZweziQlXPWjZVjhqfQwh2ux7MY392kZFF+at9Hh84HyeVLL7dEFmrl5mrT7L4n0XhLN0205xJoYNz5w3TWBZ3CyNtqyj2BbBFeUn/Nn05E99eZpmN2TIzc6H1927HFuvoJnNWPGHbl5b7vLv0kX1H0sMBXSIFkLWdI/nhjNNR0+XCUn7/Q21y4+VYsAvN/MtfctLRaMzZ/AJvgNVrVvG562G+c99LI3KIp5hb7F9wrX1OiWc/eSTAOc42B0eepJLFNfZveXH2Ii54ZTFfrjnA5W/9FvUx5+57nnjFy8X2XyjIz2Ptvmxe+Xm7EZWbBVSNZHJIIh+nKYNz/xfrGP38Qs57eREv/bgtojtFVTXj4L47q5Bv12fy7frMcIEmMPP3dM6etpC9Rws5VuAlu9BLgSf8oS45aZ8/oNLmyE886vwvrTa+zpKv3jHuyyv28+sqa1r+tzXruOL1Bcxetddoy7PO1/jY9Th3qv/l0a82kuCKTK2DHpSVzIbtySrg81X7LDPKmlOxobOjFUcUluw8ym0zVrP/WHgemvNfXsSGAzkRgdujX2/kp836D6P5x8ocHHn8Afqid4ksCwYF/zyrE0qCHjT4cjKoEyxM7dCqRUR7NgUPdKF5cWb6Rxr32RWNAbatvOt62rgtXvHSzqZnROzJbdinhbNMIb1sO/ir/QcA3gqcy5uB8wBo6NZo0KwdAEmeA0YNxe9qJ2OEF8BuLZXB51wd8bwDbVtIsx0mXW1sFFKHAjOtbipbtDTL9vsaDAQwAjIHfhoHMyShrjrA6BYpL4PjSKhf6v2gjy7LOxDuEmpEtmVulpB01fqaHdXC+cVQMLZa1V+n5atXGt0WTQP6czlTOrIsZTxfBIbydWAwp3he4AX/RcZzzFd7k0USh7Uk47aGAy+1/M0f1H74lHB7vw4M5pDWAKcSYLrraZ5xvsGNjm+Y4XqM3sEDebraxPJ+71IjsxJmqj7TTpnbmHk1UwZHCVgWIgXYqzVGw2bUQpWmRcfeUW8v2V3XSjlIt2AQs0rVA/4fVT3AiS+xHEfPYHfXgmCAk6octWQGuyh72LNuMTNdj3KabRWNyeZfznc5xb6Bqc7/mJ4p8oBr/m3LKfIRh4e7HTNQM9ej+Yrou/e/ANzjmEmRL4BSzmt6u+NTnnW+yrHcPH7dmaVnlI7utGyTVqLL0W7Tl7W48j9LcG39CsCY9ylwZDvHCrzUL7YG61m711gm0TycH/79i9ZtBPrJ6VDbeta7r+XSzbdCwBeRkSlPoeX18tIufwVOJUCiUsh4+wKuts/ln87/8aDzA061lb1w8IkiXVTH2cGcYlz4+NQ9hRbKEf5q/4Gx3qn8d8nuUh9T3xc+y8jbtojzvw3/AAZUlfN6NsMRDGbyjuznR/edFOPiBvtzxnbZyz/lVedSXjpyEc/OK2D9gRzemNjfuP/BLzfyvxX76d48kQYJ4ecv9gZIjNO/UPfM0s8unpu3la/WHMCvarx/7UBj2yKvNcA4Wui1nJ2N9C9mBvr2ecU+Nm/dynnh+dXobtvNo87pZH2eyJ7kuaSkNOUiu16oerb9d+5dvIt4Z/iAq6oaNpvC7iMFXPzaEv7St7lxny+gcuq/5wNwTo/w2a15ttlQStmU6GLB1sOkcIy/2H9hWUZnHpptt6RTr7J/z02Or5j83rWc/sQDDNzzBje6P+Np/6Us2Z7GjkP5zFq5n8b13Dxi0wOc34IBzrD2yWhrG0AxxOXp2SsVhUKlDpBtee1WqB0t118KXMhmLY1z7b9iQ6OPbTspivUxoKfq4xq3NTIyADvVVNraMhlm22DUxcwMjDTuz/NBUR39tWvkzyRfaQTAMUcKu9TmdELf13NHDqPdoHP5ZsXfKDq4nad9E+hp28m19m/ZpTXhSf9ltFCOcIF9ibFvSkpn9mdbR3bktjgNDr9BJ9s+htrWsyd4QPZpdnZpTY1gJ9TlUayFPyRFcSnEF4cPAsW4aNYoEY/mMEYJ3ey9lb87ZrNKbc+F9kXUUTzYDoRHtjhQITNyHpB1Whvm+3tzqf0n7vNfhyegkKslkKiER0mtVttzmn0NzbXDvL5wB92aJjIycAAUUJLbY8uJ53bf3wE4tWNjXtlan77KNpKVXJ7xjdf/jtqG0+2r8dniOHXMpXiW32Hs+w61GdnxrWhcqHdprdfasNHRmZTAUgbbwiOxGiu53Of4CIB0LcUy0myj1oo2lD7fSqCULqrS6BP9hd+D5BLddQXBzE20LqoQTbHRqlNfiFK7X7I4urlyxBgZtS2Y9dmgtWKflmx0ly50jWCEd6HxmFAGx7K2Fvqw/eLf76OjbQ+DXJt5xneJcd9Q2wa6Krt51vkaBcRzmfcBI5AD2JyZy/7sIprXj8fjV7nW/i03O76EH77k3b2vc01wu4ZKLru9AQq84SDAhsr9jg/ZrKXxSWAkfZWt3B4cSLEqvQOXvplHvNPOhnG5lnci1EXXhKMcox4BVWPrwXy8OxeT4somW6vDrMAIrnV8i3Z0F0dyC2gdfE0WBbpxin0D//3iO7bVudx4zn2mk7ZoM8ermp7Buco+F6cSoGvRCtjxE/me8G9tnsdPSsQjrQoL87ne/jXfBAZzOK8JZxSHs6a9bdstczRdbF8I3FfOM1Y9yeAcZ8U+lVNs64wvajtbBpfYF7AiuO7MFfYfWOu+jjscn+pZFr+X+gHT7LS7rFmPO2auYeiTP5FT5OPV+dvZM/cVkpRCmijZXOCZg6ZprN9ziCecbzPWvoz3XVMBje+DqdljhV4+2KpQZ/V/eNU5jQP79/HLtnBXWSg1mbV/O5fZf6SFcpiDucVGutc8X0WRzxrlH87zWNLFQ2wbjEnQsvKtywoAXGhfTH2lgHa2DA7/8i7e7PCZWALFKMHHNuMIDcjljYU7ySv28ejXG8kq8PLWL+Guo92mEVu7jhRwq/0zFrhup3t2uNA11KdsHtG491gRkxyfcK9zBp+5H2H9noNsCfa/K6jc7/iApspRprv+zaxftzL68HQaK7lcY/+emcv3Mv/n75lU8BwdD3xu1BosU7sAek2UGqdncOoX690ahUodNCXya7bLHQ5wflM7s09L4d3AGC72TuEq7z3GnDoAr/vPNS5n0oAGifX4Sh1i3PaYX194s73tAC4lwE41lV2mLg2fquCrp2dZUtVDRgbnstP747GHD1oD+g8Cm43tnW/mLt9NHKIBP6j9uMz3APf5ryeXumwMHogMKV0ZOSgcSOdqCbiadzcKcj9yPcFM92MAHKI+x0xZkzwSaF4/Hld8uG4gvmln47JHcxDATlrDBCNAAPhWHchY71Tu919LhqYHa41zrAGNul+vuZkfCHcRfRMYzEP+q+nheZtPA6fSv3WDiKzEck1/XzrY9rHzUD6PzlhAolKEqik4G7XF7QgfFG48tS1vXDOUib7JjPFOpWe37twwoi35Z01jY6u/olz4Kjhc7CRcd7ZLS2Vx+ztZGujKE77L8OJkU3xf4/7PAqdwm/dmAKPmZo/WBKV+S2ObjWpryuLHZsxrUxFenFHrnEJC61nlR+miClHjG6IkplpuC2XDQkJBV0/bThKVQgKawm4t9BiFf3j/wQa1FY/5/sriJpdbHrdDa2bMWwR6QfwBTa/96Ui4O/Eu5yfGZbfiZ5brEbrY9tLfttUYVRjy686jxuKoXr/KYNPouPgN4RFdTZRsvIU51CnYSw9Fz8iMtK3mWse3/Nv5JqlkcYY9XON1sf0XQA82vIetxe0tlMNcYf+B3+L+zivOF/D6VZbsOMI5Nr076vvAACPosx3bRcHhdJxKAI/mZJHaQ2+vbR+Zpqy1OcCJxq9Csc9nORn9dvZHePOPMty2FgXVksE5lFfMuNeWMH3xLss8XImLn+B+50e87fo3hV4/qf5wlvRU2xp62cLZquG2daiBE78OYUknJMB59dVXadOmDXFxcfTr149ffvmlzO0XLFhAv379iIuLo23btrz+euToklmzZtG1a1fcbjddu3bl888/r6rdr5Rif8CYXC2g6QfYmxxf4cBPE45yv+NDEpVCbnN8xuH9OyF3PzZT+jR9+/qI5zyU5+Gfn6zh6e+2ENg617i9c2Ar3284yDNvvEmD4EiGZCWX3op+GqVpGi/+tAPHsc087HyfsfZl/MNhfZ2KfAGOFnjZ+8YlTHW+zfvOJ0ANf9i3ZOZyke0XXnS+ROMc675l5RbQXgkXeyYphXRR9B+bI/keGgXPBHeq+o+YOSMRl7+XQHb4zM6lBGjKUXqxhZ/dk/jafT8vfbeae2ets3Qrgb520GOOd/iH/TNAI7+omL87PqeV7RCPev9tDH8Nvap+85IBqkYX0wyw/Uzp8+7Kbkutxs9fvmdcDo0eesj5Phfbf+Fp51s4FJV0tTEZ6AfZpAQnSrz+oxvKpBTYrdmNtIb6QeKSIR0Z73mQLwJDedB3jWWbXOpwu+8WNqktudX7d2PeHIAMUqjjdrBU7cpU32U85ruCn9XeZJkCh5/UPiTXdXHZQP2gODRFJZCkX05TDpESzCY5k5qyUtF/NHO0BGjQGoA2jcMrfNdPKHmgVPgy2E0FQLvT6dpnmHF1v5ZM06Q4vlXDmb9QsH9Yq2/UuADkEY+madhdpgNnSlfjohbsEmheP55NargbrF1K+DXNCB7kStbY5OzQv4PL1Y486xvHh/4z+F4dQKjm576xnbnltPbkaOG2ejQnGxzd8GgOkpVcGnkP0DrYHbJfS6Zu3brEOcM/mXVcjuBkjwqg0LR+HPeN7cL5p/Sh6zWv4Oh5MQBbtNbGYzZqrdjfYACX+R4wuhBXNjqPF/0X8pH/NO73/Y3v1IHGZJAAS9RuXHTpdRx2NacoqR3z1d6WtuaWyKyo2CwZmZK8Nmug4jHNZBwS0BR29p0MwCN+fZRTALsl8A4FGAC2uinY6lm7zkoGOFscevAaOtDu1VIs9VertA6c453K+5yLr3FPIzP6ReIVgMJuLfz8u9VUo5s3mtBvjrnL6zy7NcDppWznMcc7HMzYi8cfoKOpcP38Etvaju3m+aL7+cr9AJ2VdAaZsm1j7ctoZ+o266nsIClYf7d5g37CujfYPZqmHOJxp96VP8q+EtXvIf1IPmcHjxnfqIPYr+gnJ46cXfiO6L/lhxypbEP/DndQ9mEnwL2Oj/jI+S8KDuijoBqQa5kIMiSggZK9zzKXU4e8ZbzmnMb7rie50/EJ+ab6wvs+W8/yPcd45KuNnPLUz8YM9g026cXzXWx70fzFJKvhE+XQ/Ex71cbkafE0UPIp3rc6Yl9OtCoPcGbOnMntt9/O/fffz6pVqxg+fDhjxowhPT36MOBdu3YxduxYhg8fzqpVq7jvvvu49dZbmTVrlrHN0qVLmTBhAhMnTmTNmjVMnDiR8ePH89tv0WtcTiSft5hRtuUAXO27h8NaEi2UIzzueId/O9+wfMg2LpzFwt+tk0aFCzM17KYRD9sP5xOHh+5KOIvRhZ288MNWupQYSjnarv99j19l79EiRtnCf+MC+2IjUwJ6BmfVymX0DkbfbWwH6XD4B+P+XQcO8aTzLc63L+XyA4/re6Zp3PXJGv71vj6Sp1iJY00d/SA3y/UIF9sW6gFOcOjtWq1txOuUtXcL786xDiVsaTvEdXyBW/HTXMniSvtcvlmXYWRYQC+Se8v1LBMdP3Cn81POsi0npWin5QAXmqfE51fxB1QKfOa+c80yG6p5Vto+NutcJqFhrKAfQPvZttJbsZ6RLdO6GJcbJDix19F/9EP1J0X2epa++y9uHsZ/ruzPHWd2pF6nU7nd93e2lqhhAZin9meM90m+VIeyw9GWLapex/OufRxNk+JoXj+BRU2u4Os6fwEUvg8MMB77XWAADeu4ePi8rvxnYh8uaq2SlKq/Bw2VfGP4NnVT+MB+If/2jecsz1PG49smhw/65vWSQl73n8sGtRXL3EOg/ZnQpLtx3x6tCWkNEnjOfwmXeB5itRp+7zO1hmRjCnC0eDT0kUThF7G1cTF00I1z2vk/3+1sUVvwgO8aOqXW47vbh/PL3acZGZyQUOYoNDfOQRrwUuAv3O+/1jJ0+oYR7Yhz2sk1BThHSOSCfm3YEAxImu75nDbB12qXlkpinMOSwanjtltqxhrViV5L9DLj+SHQh/GeBwlgt0zu5rApJNVN4Dn/eO7zX08RcXhw8bpfD372acms1dqS2KILje/bSPwdK9lbosD8iGYNogPYyszI+Bx10JTwfntxomEzVnUHPcjO7nUD/xn6M/PVPsbteaaM169qOBhV6qZgNwU4Xs1uFB2H7IrXtw+d6W/Xoo+odNgVmiTFcaP3Dv7pu4FFqXqAtc30fHuc7chIGW5c/8Q/wrhcqLl5zh/uqlKDJ5pDbBtJ4Rjn2pbSVdnNbPdDTHT8QPG8x1G9RTQxnYAllFgWxLXxU1KD68D91T7PMiquvbKftqbfFLuiGb/BgSx9u19U/TvS0zSEHqBJ4Tbq5myhsZJLvhbHErUbrsb6d8adtw9blv6blOVqzkVnnQnooxPH2Rdyk+Nrhto3cuqWRznTtoLf3Lew0H17cELFML8Grmz9dysjOBVEe9sBhtn1qQn+7phNfrGH2av3859fdvLLNr0brZuyi8mOD7n34yUQ8KOo4fqeFsXbSdb0Linz5+83rQu/BjPa/m0/cbJVeYDz3HPPce2113LdddfRpUsXpk2bRlpaGq+99lrU7V9//XVatmzJtGnT6NKlC9dddx1/+9vfeOaZZ4xtpk2bxqhRo5g8eTKdO3dm8uTJnHHGGUybNq2qm1OuzgUrSVIKOaTVZ7HanSnBOR4mOOYzwr4Oj+bgh4D+g5G3YxlfLdSDslAdQivlIG2UDOa7JvGN6z7ig9PTH87z0EvRD+RHtbr4NRsNlHy82fvpEDzzCB0EzwwGNPuOFRqrJIc0VPLpbAqIin0BEg5Y10i62fcul9p/woWPtEPzje6BVN8+Alm7+H33MT5dsY+Wmp6iPOxuyZqU8wE9kn/a+Qb2g+uMDE6oC8csTTlE0RHr3CgpHKOnKc15uf1HEigmlIv5q30eS9z/sHSLnW9fTDfTzKkAl9l/or+yGZ/fz4T/LCMrYweTHP+jnbKfphylrhJO7w4zzUrbq8RQ3uF2a8bqL/ZFlhEiqqbwif9U43q8046rnn7wCf1YehzWg0+jum7O7NoEh93G21cPIJr/G9nOcvnVK/pxjfduxnkeYpWrH3FOOwv+OZLPbx5mFJq/5L+ItWobnvaNZ7nWmVM7NibOaefUjo1x2aFH2xZGcBEa9ULdJuQTzyuBC8kkHCh0aZrIOT2aMqB1A8b3jwy+cqnLOd6pPNPgYX32R2ccKxPP4JBWn3/5/0pqUhw+HPyudeaDwCjjcUvUrsYoMwgfLIvM3Snx4axAaGkAgBduHsdZ3qf5IDCKpHgnnVMTSWuYwAGsAc46m7W26WCJkWJ/6ducn+4Mv2c5hAOcLC2R5g3iKe6pd/nd5vicJx1vAfqoN4fdhttU0JXgclimdGhYJ7y/ZodsKVzn+yfLtC40TYpjwoDwa+q026IGRq8ELuAcz+OM8zwc0cWZX6Jb7QhJlusB7EYBdzSqPQ7FGX6O0Agqc1CUqyUQ73KgxdW3PNacbf7V/L2uk4KjTvi9C2DnnFP6GdcLNDcZdcOBMMCOUgIcgCaJcWRTj08CI3G54vjfjUMINOxg3L/L0ZbDbS5gn5bMBrUVD/qvYU0wmP6KEcxT+7FKbY+qKfyf7zYytIbEK16Wxd3Cy66XmOMO14c02/kpCYXRZ0QOHbw7HPjCuK2XbYdxAgN6d2bo+kd+fVTc3+zf4sZrzLwc6l4qWVfXtGgraTn67/Xvaie+m3QGLVt3wKM5sWk+EjP1rqucuBacO3wgHlsCLiXAg473jecYZNvMf1zP4lICJCgeHna+Z/kbX6fbULKCy03U6c6x5L6UpGRt57YZq/nXN5vw+FXsBPjC9RA3Or7hJj6FvAxsWvgksq13Mw5Fxa/ZmBsId1HPDfRjcTCY07b9YCmGPhmqNMDxer2sWLGC0aNHW24fPXo0S5YsifqYpUuXRmx/1llnsXz5cnw+X5nblPacHo+H3Nxcyz/Qp3I/3v/OKNQXhvs2MAAVG1+rQ3g66QF+VbuwXO3Ijb47+DKgZzva2Q4Yk64tCn4oWiiHud0xi9a2g3S27eUSuz7Vdl6xn/42PRW5RO1OerDINMW7l47BbqK3Aufg12x0tO0nTTnImc8tZH/6Djra9hPQFFYGRywMtW2gnbKfdsp+8go91MnSA6AvHGeToTWkiZLNk87/8JzzNQZjrW144uXXjaUYQqsZH41vxYHk4fzTdwOgn8E0O/C9EeAUNOoR8Z40V45E1Oj0tO00fgCKNBctbYfZGPc3fnD9k7ed/+ZfzndppORxTKvLs75xgH5WFuo7XxjQ/85Y+zI+dT/KlTmvsXZfDu+4/s2tji942vkm7W36PodS7T2UnQy2beRM2woj8/NdwBp45Gt6FuMsmz7r6qp6p/J/3tv4m++f/GbK4Pj9fmwJ1oONx5mIeQRHyc9LNO0ahbsQbjutLS2S3BwgmeVaZ+KdNnw+H5oaQNECxlD3DBpxvvdxXg1cyA3DW/OPkW0tf0NTAxxxWGskfO6GlrmNQturAT/Txvfgo2sH0DGl9MLSgKoaj3mnyf0M9rzMPq0xdlOG8MvAUHaqqeRq8XwbGFgig5OAqmkUu8NBit8V7morCL7uPp+PrqnhQKR380Tj7+5Uw7VGWt0m7MfaxkxTN8r5PZvyxAVdSavvNh5v7t7J0hKxK3Ck7V/4j38MEF6Ne4naTX8tNdN0+jYNly38+iXF2aL+JpjnY1pw53DqOMPXnXaF5vXDgdETF3Zl4uCWgMIGrQ2ZNOLZcT0sz6diI9fUhXVEKxHgaGVncAJ2N5oznJkLFd6aC3BzqINDUSPm2jmk1Tcub1fDGZVAncYETBN5ZmoN6N8t/N3YrjUnr641k1sywxOiaRqNEsL74rJDnxb1aNDjLIo0F/laHBvcvfi/UT3ZOG4+uVd8y0c3jWC89yFGe57ig4a34MHFxd5H6Od5je/VgZZarJKc+Di2WS9o9pdYpDQ0zL++aWRZD9tuS7ZngG0rLiVAkebief84srR6dLHt5RXnC0bB9lpnTyOTZNbIs5dmhXp3V+Oup9KqgZtBbRsav+/NsvRjWn6dNHx+P0cSuwFQVylG1RTmBqxBJOgBzwW2RYyzLyCRAn49ZGP3Fr2r7KC7JdpgvUg+V4vnaPCEw3nIOuppkG2TUQN2kX0xmzdbTyI7efSTv4M04BX/BexWm7BdbUarQRewvs5QVE0h6eCv3PXqzCo5zlZUlY6iOnLkCIFAgCZNrH2zTZo0ITMz+oJlmZmZUbf3+/0cOXKEpk2blrpNac85depUpkyZEnH73LlzSUgoe+hjZTTK28Qp/t8IaAofBs40bt9o68yr3geN692CKcR2ygHSFb0dG2jPKWwkTvFwgT0cqJ1r/5X3AmcBMMKufwh/VzsRZ/PQlkw6KPtoHww0lqsd+V3tzBD7Rv7P/hWP+f9qPGad1pY17gH09W3nQecHTNY+QkXh3UUFdMrSt9nu6MgzBWdzo/1rJjp+4Fx7eIXldWpreth209m3EdDTwaGhnts9jTiwZwdfBkbiIMBU59ukZK+loU3/cnsCdlaoHehn28Y+JZUWWiZuxW+ZJwXCffM71KY87Z/AC85XiFN8tLcdoD16CvjfvvG8GTgXDbjR8TUNlXxjuveX/BdRQJxx/QplLl8rA436pH62bUZh8C9qTzoqe2lrywyuqK3zaE5+jj+Ls72/G7d9ExjMBMd8YwRHnr0R36qDIt7/OXPmkJCVjnnGkiNFCiNSDrHSbufMZipz5swp8ajIr6AvfTVOm53kOPj+u2/J8Ya38xTmW56jsMgOJYavdvNv5+cfrF1p8+bNo6HakPbB637szPl5CUXFTuPxkfsGe/Ki7yNA1tFjxmP2HQivf6TfFsoKOLnQ+yhufBymAZ0bJxgDyvKJp7iomHqt2hm3LV61hW51OpFcsIX3AqMt+3VHd9iVp+A8sJo5GasB60HyqJZIui8Rc/IilMFR0Dijzl6+/87anWvO4BzSGrB10waOuOE1/19prWRypn0VuVo8i9XuzJkzh3UZCgS7uhb8OC9Y36W3dcPqlfh3Rw5HHtDAxtxCGx0SVb799tvgrfpjtICPhIPruKydQmq8Rp2Da2lges2vaBfAvm8Vc/aFJ5K7vpNC7u46JAaLfyMCHGxo2CwjzwKaYgRr+R4Vp6IZLQ/V65jrdnK1BFYtXMCW7HB7AZ7zj+Md1zP8EujOLi0cTP5yrAk5c+aQEOjLKPtK7vNfx8TftxKaNeiwlsTmLOsozPVqm4jXCsAfCLB19a/Ga7A3fQ9z5uxic7bCrZ7XUbGRWuTk++/017IIyPHqs2Jv1dJoX1DA7d0DTFvv4Bh6BubDwBlcaF+MEz/X++5koG0za9R23OH4lE62fcZyK4vV7pxqDx/sV6odjKLhksyvL+hZvsPU5xbfbXzofJwz7fp7dpBGHPHHk25PMSbHDI0Yc2XvJEU5BDZYf8zJjjlz0DS9q7cD+40FZnfl2pkzZw51tNY0R/9tWqp25S7fjXykPEEdirjddwt/tf/AJY6FvOB6FYC19rlc4H3MmGJie1Ei+el2PvQ8RobWiNsds7jC8SOFO38HwgHoINNkiilKNs9/NZep5hGxgU2g6CcQ+2nMGd5naBIP95DOrrp1+bmoN2fYV3Fv4bN893U8qq30qSAqq7Cw/EWfQ07IMPGSMySWN2titO1L3l6Z55w8eTKTJk0yrufm5pKWlsbo0aNJTKz44mXlKhjIpy+uIs+r8ffx5/DV2kxuGtGGxTuymJ+hH3CmX92PbfvSYBE0UvKMmprCui3JJo3UIuuBqZ+ylYbk0tW2h0G2zcGovT/NlCOcySpOta0lTvFRrDlJ15rwTuBshtg3crnjJ8bYlxkH9231T+GSCyfCex8D4cm4zlV/ICWgZ4Ca9z6Tp9p15oq3U7CjcrlD70Mt0Ny84r+Q113TONO+kteU53EQMH4U3O1PYWDTHnyZvpHVwSzRELueDQlgI7VDT/75240Mtm2iUd8LuH7tBBKVIjoFu8o2qWl0se01ZqfdqTXle3UgvT29aEQuVzu+p6dtJ+/7R/G1afTQMrWzMZOwR3PiSenNLQc70sp/kFecL9LVtoe7nDMtr+dfgj9W27VmZGmJtLVZg6zv1f5c8JdLYUY46PlR7cME5hvXOw8+G76IfPvHjh1L/q4U+Ojfxm31m7biukvGcvVfVGOov9n7B5axfE+2cf2Mzo0Zd34fzhrtx2FTiHfZKfT6eWiF/l6kNGrA2LHhAt5/rZsPXuvcIWPHhiej8/l8zJs3j1GjRvHLzq8gRy9mzHMmM/acc5my7mcInhGZHxey9WAez61fGnE7QP0G9Rk7Vg/0diXsZPWP243nuW1puBg+15S16dCtNwSXsCrWXLjj47hi/Fn4tyehHNnM0EG3QPEVXPPkm8bZc7T9Crl76dfG5QbtBnBweX3jerHmJIc6dG+g8si4gfRq2TDi8U8tCy9IukNrSu9ePemSWo/XNv3K3323ckZgFZu0lhQQz9ixozmwaDfs1tP95587loCqcc8yfe2rU08ZQt+W9SP+xpl+lV+2H2Fg6wbUC07JEHp96ibEc8G5I7jAtP3RAi9vb/uFNsl1ePDKQREzco8Ftj92vxHXZkWpwQE9YHEHC+6PUo/GwZo4V536JDiccORIcDuH8XqFnjOXOpxz9ihcmw7xv53hs/ef1L582mc6Tywt5hiJ/KvBv7j33J4Ma6kXnvdZ6qWhP490rQnvnT+ad7f/xl8L32dWYAS9e3Qnd36CMYvz5ii1ZwCKzc5fLxzNY6v019VZP5WxY3vT4WA+r23ST/46pqUwdmxv4zF5xeHvSLMmydwyoR/T1oc/g+u1tozwTKOuUsQurSk/B+uKztOW0Il9xrwt6VoKv6sdGWDbytxAPyOTEhKaVgBgo9aaVI4ak0SGutx+VbsyMzCSyx16Dd/mOgNx5DtZr7amdXB4/9eBIdzk+IoutnSaBJfcaDfsQnp31bNem458B/vCI7Pa9RnBmacMZcOWdhz45EeSyWGa/2Jyqcu53ieM7Y4F6jHSvtqYiqGnbRcjbWuMk2Bbiz5ceO4Y7vldf21DhduBfOu0AwMU68za5lpFCHe1HQp2EQewE1enDmPHDtOD2tyeaG+fRsc2/Wg3ehS46nC8hHpgKqJKA5zk5GTsdntEZuXQoUMRGZiQ1NTUqNs7HA4aNWpU5jalPafb7cbtjuwfdzqdOJ0VH05ZrvpNeZxrOeb3MrdZfc7vo3+BV+4NvyE9WjTgYJ7XiOBDc0T46rXgmNrOCHAWBHrSSMmlu203n7keprVN/wCuqnsKGZ5GeOu3g3w4I3iWsENrRrzLyTxvfx72XcWNjq9oFvziBTSFrmdeSd02/Vjn7EEPX7jbqfkRvdA3W6uDK7ktwzror+GrgQu4xL4ApxJgdmAo89VeFGpuGpgyJhCcbr9JT+rF6xH6dq05fs1mBFBHbY3wY2en1oydgWY8nNqSI2uSSFSKjDZt0NrQhfCZdWjoaDFu9tOYx/1/jfpy/6L2MAKc1Vo7hnRuxtqDO4M/YL3oatvDgBKTjHUIdlFtV5uzTmtDGyWDjVorflF70k3ZzTuBs9nQsQWZAyeTumwqD/uuYk2J0SBJnYYB1uH8oH+e4htZ6wrc9ZsHP2dRm8BLl/fl8W82cdXQ1vgDGl2a1sPpdNLQ9IBEh6nOo67b8pm9amhrnimxMGm0z7TT6SQvrjmhQRaH49vSwOm0zAgd7XH14sNdGfeP7UJA00zreCnGY24a2Z6ABmd2aVLmd6pdz2FGgJNJg/Df7TIGGKPnClwpxgGotP0K8eDiJf+FtFMOMPbsx+lSsNCYi0X/HClc1ylAr5YNoz6PeaTPdq05Hex2erdqxJTzu9EmuQ5XvhP+3XA6nWimkyin04kTGNi6IQfziundqiFOUxFyeDs4u4e1O6ZniyTW7svhkn4tIvarSX0nv953Bk67jThn5PMB+ExZlSysAU44IxP+3GRr9YyDnt+RgOLUIrbPpQ6gF5jmaHWoG+8mzhV5iDjaoCdH0T8DGY2H4WgXrunIJ4H8YLef0+nkm8QJPHHsdHw4OK1uHDf67uBF50vMCoxAM1VIuBw2Y3FLr1/F7Q6f8Rf6AjidTpo3CgfKzerHW163JHt4P5smJeB0OplyfjeemLOJS9v4+O82O4epz+FgF1vrRgn0admAfq5BsGaZUVS8X0vmP76xPNG/kFt/b26Z26VAczNX7c84u96dtUltSYHipqld32azaaTf8/5LGGFfR2Oy2djkXJQCvbvrXPtvZGn1+ChwOjc5vjLq4bK0etRJbmW0SWnUDoKDugKaQmrrzvpr0KoDp3ieowF5lrq5kHStCWd6nqGlcoi/2H/hGsf3/NX+A42UPFRNwZPU3vLapgcDHPOaZ078Rp3jMa0uDZR8I7seysSHZNnDBe8BTQu/J41awY2/oCQ2K3cNw8qqzDG7SmtwXC4X/fr1Y968eZbb582bx9ChQ6M+ZsiQIRHbz507l/79+xsNK22b0p7zRNInjVOIM/3QmWfAbVTXTWKckx2q9UBYr1lH9iWHp3n/We3ND8HZPUOBwKzAcJpd+Q7rp5xFl+79LI/frKUZBY//DZzFcM8LvOo/nz1qCg/6/0aLdt1AUXip0UPc5buRXsVvGkMpQZ+czB0cETKoTUP2aY251ncXD/uuYor/Kopx80VgGCXNUQfhTkwhpV6wXgKHaW4LOOZItkwQmODSf2jMVqgdLNfNjw9x2iO/JLMCI4xahHf9Z9O7Rfh55wd6W7ZdpnayXN+uNWOH1pzLfQ/wL/9EFqi9eDVwAcW4sdsUUsfcw9Od/8d/A6OhXqox5L9Ac+OqHz5YxTltXDOsNd/epo/ocCc1NbYFiG8Uvc4gpGlSPC9f3pcBrRsypF0j6idEpnLNmcmSBak3ndqOD68bxCX9ImdLLimvTngulSMJ+oH9ikH6UNvhHSKXfgAsw6Ib1HFx06nhgECzbGfnztGd6JVWP+rzzL5lGIvvPZ3UJqlwwwJe85/HnMCg4zKh+7P+8dzsux0SGvL3CeeE98nt4sJeTSnr93WdaYTfTq2p8VpfNbQ1Izo2pl8rPQhrlqR/vqOtKj/jhsH8OOlUywir8vz3moG8ekVf/n56h6j314tzlhrcgHX2Z/NQdwgHNuZC42OmLJrXWQ+c4Roej+bgyiGt8DnD9U+5JOB22CyzpYfYFIV3ru7P8A7J3D82cgCBmdtpM2p7EuOdLFW7McDzGk/69XluPr5+MF2bJjLjhsERj502oTfN68dzz9n6cPHEuHAQU/K1MWe5UoPv1VVDW7PqgdPplBT5np3WOYXnJ/SmURtrfeB+LZk9WipH219MMW7LDNKZWkPLdAUbtVZsNE0BYB5Vdpj6nOJ5gW6edyhMHQAafBYYziTvTZzteYp0LYVC05D7jWor6sWH36+WHcIF2Tu0ZrRK0T+HyXVdeHAZwc2cW8MjyUJyqMs6rS3zgsu0hE6Cd2mpuBLqWbbdEwxwQsXQdztmsC3uSuIVL1laPeYFa3ySglmrtXZrofgxU4BjXgdQf1BzyvzynQBV3kU1adIkJk6cSP/+/RkyZAhvvvkm6enp3HTTTYDefbR//37ee0+v/L7pppt4+eWXmTRpEtdffz1Lly7l7bff5uOPPzae87bbbmPEiBE89dRTXHDBBcyePZsffviBRYtO7gqmmqZRHFwqwHxguKRfGm/9spPze+lBTb04J1u1ZpyKnhbN0RJo3qwlBwqS2LQxjQNaMh8FzqAuRVxgW0xL5RBP+y/ljcB5bGrQiHiXncS0LhAukWGD2obTO6cw4/fgMgjYedp/KU+jTxU/JXQmltCQTwP6KJI56iD+btOn3J/HIG5qqX+J3pzYnx82HeTOT2Ah4eK8R/xXsVztyA6tGYe1+pxhX8kXgVN4Lc7J4LYNefyi7uw+UsC2X5sbNTM5zhQ6p4a/VG6HPWJY68oSAY75TCgp3okvoNK9WRLLSqxjlUsdLvM+QCMll6VKb6a2DZ/RrNSsz/mx/3QGusKr5u4sYwQHAIrCLRedQat2B+iUmsi9r1/PFfYf+EQ5m8cVhTcm9uP+z9fx3PjejDAtzqk4XByhvpF2TkppWdpf+EMalghwHHYbw9on0715EscKfVzUp/SA6mCDfmxVm9PRtp999fVC6tvO7ED/1g0si4iaxZmGQZvXPgPKXces5H43rx88qDbrzVP+y/TnL+cpoi2YWpY4lxN6XwGrP6T1JU/y7zY9mDNnb6nbe3DxN+9d1Cef3VrTiMn4X7qsD28s2MHVw/R6kZGdUnhm7lZLoGmzKdgqsTQC6MHi2B5Ny9+wFM/4xjPV+R+22dpEjJgKBTheU5eTefSa11EPnOFuzQLiaVjHpd8eHCiTq9VBUZSoAQ7A6Z2bcHrn0peMCA2Fd5gWcEoyDuDh12pIu0bMCZ4cpDWMZ+/R8JxXF/ZpzoWmz7M50O/QxHqgNutous9pt+GI0oRQgORMsZ74hNZ4O7t7Kmd2aULfVvVJt71C/R/v4t/+CezUmuLVZqBhY7HanXmam3H2BRThNoLlW09vz4s/bQcU/DhIa6BntALY+UwND2cPLbgLsElrRXfTiLyk9uGu+MVqd64OzkelKAoN67g4WuCldaOEiN8Ds9/VzuRp8Ubt4AatNXVLLOS8J9gFV18pYKRttT6Tc9C3gYGWGdMBMpN6ox77wig+t9VvQXDKn6irtJ9sVR7gTJgwgaysLB599FEyMjLo3l0v1mvVSj9zzMjIsMyJ06ZNG+bMmcMdd9zBK6+8QrNmzXjxxRe5+OKLjW2GDh3KjBkzeOCBB3jwwQdp164dM2fOZNCgyMLPE8kbUI2Uv9t0htGyUQKrHhxtBD2J8Q7WmQrsdtOM4Z0a88MmjTHe8HwkR3EyxvskdSk2hoLGBw84zZulkavFkxj88BYnd+P+c7pw2cCW/PU/vxlrTA1q04D+8eEJmUI1AABv+s+hm7KbQc2c3D7+PhoGD0BJCU4u7teCTqn1jGXv7z67E5k5xew60pQ1wZmQ3w8WgdaLc6AoClcMasW8jQdZsaSd0Y2V7WrKlUNbke/xc3rnFDJyisgoMXS3WZsuYFp2a4PprOjnu0bi8Qd4/JtNRLNBawMaTOjfwlgZGPQfk4nee5nq/A9P+i4zJg0DvdagkDgGtmnIjkP5ZBV4oz01ddwOJgxoyeE8D58ERvJJYCTJdd08DpzVLZWzukVmmgAOavWN1bOTUqLXGVTWRX2a88Omg1w9tHXU+5Pinfznqv5R7wtxxSVwofcxuiq7GdZQ/wF12m2M7FT6JO3mpTPsJc7Goq0uH/L1P07hlZ+38+16vSu59E3LjnAqG+AAMPYZGHILNOlm1BeV5Sc1cthsSLP68Uy5IHzW2r15EvPuGEFKlPmBTqRlWhfO8D5Ls6Q4xjfZjHnlhtAsxuYuKvO6W157HXCGA4lcLYG/dErh2IpwJihUfO0qJYNTnsb19OyEeemAxvXKLjR98dI+XPL6Um48NXLerJB3rx7Ab7uORg3kHzmvKxsO5HJ2d+v3MmqAE8q2JVunFdirNeaqIa1w2m3G92lPVlN6zqlPKDDr53kDjfBw/TM8zxBAz1Q1SHAyaXQnbjuzIw/NXk9yXTcX9mnOY19vpKQtWhrdgrMwL1G7cY0pQ0VcElzxKVnzX2HYyH9ZgrtXr+jL7NX7GdevBalJcTxxUQ/quO089vVGyyKcPhwsUHtyrl2fiuRXtSu9SwQ4HiWOQ1p9UpRsbjYt/Puu/yz+7Z8QUXtTt1knso6F67lOHXEa//5I//D5o2Q3T7YTUmR88803c/PNN0e9b/r06RG3nXrqqaxcuTJyY5Nx48Yxbty447F7x415TSNzBgfCgQlAYpzTmCsAoEPv4SQkxVvm2DCeEzfFRNYPpSbFMz0wkmsd3+LRnDzxj2vB4aRXWn3LApof/G2AZXRMwzqmkRLU5behb3LamM5Rl9ELpXoBujdL4uaR7Sn0+hkyVV86IqSu6YtZx21ngdqLe9GnOt9UfySjHHbuGKX/kBwt9Fqm+s/S6oEzgaX1zmJI3vesVdsY7Y132o0zlPLS/02SIg84v6g9OcXzonE91H/8euA8zuvVjPvGduaClxeX+byhNoVUJON61JShciRGD4Iq67nxvfAG1Ep1g5RUx2WnkDiWa505rYzuDzPzGXzJYKOsBE735km8ekVfbvpgBXnFfmMG58pylhPgKEqU/XAl6MHNH1CR97es7MGJ5nLYuP2s7mCa+sRnGsEWkk14n1ukpoKp27xz6+b0bJHEt7Zwlic0fD5aBqcir1Gz4MmSeSHcxHhrpumW06y1bX1aNmDdI2dZfitLOq1zCqd1jh6Qh7JsJUXp3TYW4sWVgA8HzmAx9tQrR3FGF+vz62v16dsn13VzJJixGNC6AQPbNOSVn/Wir7/0bc7NwTms7DaFxy+KnB7D7A3/uZxj+5XD1OfiS/4a+Vp3GEWjDqMiKm0Gt23EYFO2+vJBepbYZbfxfx9aj5vP+8dxim09OVodZgWGMywY4Hx03SDmrM8gNTGOPfNTSFGyGRichuRm763MUfUuw5Kr2Hfr0o2s9UlGPVf37n0BfTRbdczgyFpUx5EneLaiKNHPfEIS450cogEz/CM5mjyAhLMfBqz9yi9c2rvMv2W3Kbzg/wtv+cdynvdf4AifHd16ht49c+OIyDMhc5bjpztP5d4xnSO2MbY1L8oZbFuCy8Gie05jw5SzjPuamM5mnXYbm7SWvOs/i7f9YziSZO2zdTtsln7t/VoyLoeNz1rcw63ev3Ob7+9cd0ob6rkdloyE21n2R9XcP1+aO3w3c6f3Jl70/4UXL+1N06R4y6yypYmrZFDx38BoNqlp+sJ/9uNTxK4oyp8KbgDLrLtlfT5L06iu9Qy8vB4qRVF4Y2J/Prp+cKkjHMt7jtCBsjTXnaIf1MZ0Pz6BZFL88Xm/qtrpwYP8tae0AYc1uA9lbswrip/SO1wrU79hMiSED5A3ju6LoigU2U0BTjCD44gSHVQkpxZ638xrJiWYCpbvHNWRf54V+dtTVnDzZ3x4bX9Gdw13qZk/jq8kP4CqKdzuvZnkeu6Iz2o9029Lk8TwyeZTF/dkYJvw6/ivC7vTPiV68Gv+mF81RO+92KK15CzvU7zR4Q3O61P6shMVNaZHUxbfezoTTJNz7tCaM8IzjTHeJ/HgMk7WhrZP5l8X9iApwWUUGof8qnY1Mre7SgQ4Azo05R7f9Xg0J78knAmmLsiIGpxqQFYTP45CGZw4h73MYfChg/G9/hsYdfWZEK9/acwZnMZ1w1+kG0e05dedWZb+aNB/hKKNMPrH6e05tWMyvVrUR1Ota/SYg5a65QQF5jP27s3Dc22EurkW/vM0iv0BS7+unr5WmOK/CoDHUq31NnFOe0SAs+tIAUPbJfOJqheJTxzSivvGdsFm+vuhfmwz88iLxLjyD0zpWhPjyxx6f6IN3S7JvB/Rsmwl/az24WevPgrornK3PnHiTQeYaEXbpXn64p7sOJJP/1bWrsXj8XNW2nN8eN0gXvhxG09c1L2ULXT/PKszp3RozIDWDcrcriIuHZDGaWV011Unr17Rl60H8+jeLAkyj1nuC2VuFNOr27nHQAj1NrgTLQEOcfp3VDOFLluCdXDRMzjlf3aGBDMMlqy26btzonszBrZuyND2KbSZHDnX067kkbTd1wFQuDHK99tht3HZwJZk5XtoXM/NhgN69sLttDOwdUPcDhuN67krfCI05YLu3HdOFzo98B27taZ0Mk2e+Gc1rx9vjjmAcLAKRNTguB021qltjLl+NqqtOEoiI9s25K7RnXjs6408k34Jdzk/oXj009R1O1ijtWeQ52UGtWmNucTZr1a/DI4EOMdRtALjaBx2G1/9/RRUTaORKZAxZ3CS64Vv79y0HpPLGa1g5rTb6NdKLxr1lQhwzB/wskZphKx44Eyyi3xRz6RbNooMOnqn1eeU9sks2n6Epklx/KVEUOay2yxnDDu0ZiSWGDFSL85pCSoArhnWmvUHchjVpQm3z1wNgHmTeqUEazal7B/TAa0bkn600HiOS/qVXTPTosEf62qpDswZHGcFArWQ8QOivyaVKTKurGHtkxnWPvrILjOXw8appiLvP+PJi3sel+c5EeKcdnqGRg3arV3YoSHk9QjX2WiNTcW0znioY3pt4/STFxfhbqt3brtQvy1qDU7p+/W/G4ewbFcWlwY/M49f1J37P1/Poxd0s5xMaMclPK6c0gIz/aRPv69kIX3I1L/o3U2PfBmeE8jtsBHvsrP8AX1S15K/WWYtGsSzOdO0pp4pGCr0Hu9Vt0vfjzpRApyf1d48Euzj/DE4PUMdt4PuzZN4++oBrNjdFl/i7cSldjNSX9nUw+GyZnR9ksGJbe0a12XFA2dW6I3u0SIp4jbz9yPZFPgkRJmLAvQMS6CSp0Lmg1xFshGN6rotQVh57DaFD64bRJE3QEDTIr5QHn+APBK4ynsP006LI+PwAB45tRvfrAuvxmuueQmJc9p55XK9GDQU4AB0b57Ipow8hpZyMFQ1iLNrFAcU+rVqwIo91rPdh87rSuN6bi7q05x2jeuUm9FpESWTVJqTPEIywp/toiqpa9M/P0lmVQZJtYYj/P30anZjfhnLukd1TV14viJwhbujcOvv48/1zqNl9m/MCJzOk8H31umI8iEu44M9sE1DBrYJj8i7YlArzu6WWqnfkBOtfUr4tajM9yIUDNWrQPb41Sv68vCXG/j7ae0j7jv+AU5Y+8Z16F8vlxk79e9+tAzOHi2VF/0XkkghL/svBPR6vdD2p3ZqAkSOmCuZsUquW3YR+ckgNTjHkd2m0Kiu21KcWxnmIi1zLUDJD2VIXCXOwkPMwdLxOMiVJt5lj7rf3Zsn0bZxHZydRtFg1CQev3w43ZsnWWoxyqs1Cc1JMrxDYz77v2GsfXi08XrVT4j8sXm4b4Afbj+FO4OFznVMB/qkeCf3julMp9R6ZQY3fYIz1E4cXPG+8srW7lS1rs0Sje7RygRqJX39j1O4ZlhrHjqva/kbl0PCm+PAHOCYCostK2IrCvT5KyQ2h45nWWvDggFOjjOFC73/YkbgdOOu6PPgVG73ogU31SmuPa9neMqIaL8fZuaRgxU5QQxp27gu7187iEFtS5YMW+t6jgdz/PntrcPomxze55JZ+9Bv7XP+8TzivxpPcH6lkiem0YTqImfeMJjeafV59+qB5TzixJMMTjVltymM6tqE9KxCBpQyR4nbaaegktF/n5b1SWsYT7Ok+Ar1pR9vboedH+44NSKdW5mU9cwbh/DZyv1cOaQVLofNklb++h+n8PPmQzw4O5xKTnBAq0YJtE918t7fBtKyYeUP7h9cO4jM3GLaNa5b/sZB5RVGn2gp9eJYOvkMMnKKSi2GrIjuzZMsNVl/RnU60NVYpiLjQFnnrBe8AqqqF4Ymm7qsgkUbKfUiD7TOkgUdgFLJOX+qo4amWsSkBCdzbh1OkS9QbjbGHOD82RPE968dyPTFu3nw3D9/omBW8t1x2+GfozvgVxVj+L5xXylBWmkn1WahE7hBbRvxxS2Rk8BWBxLgVCPD2ifTs0USPYIHjzcn6rNIlhaI9Emrz4+bD0W9rzRxTjs/3znyuE+fXRlR+6orcaBLa5jAbWdGnwG2RYMEJg5pzabMPD76LT24MvNO4/4Rf7Beo47bUeHg5tyeTfl6bYYxZLQ6qeN2/Kng5niTLqrjwB4+WJsLhT8NjGCcfSFbm5yHMYA6FLAkt4crPoU64e/DP8/qxIGcYssonGhdVIPaRj/hqoiJg1vxzboM/lqJTOjx9O9xPVm8/Qh/6Wud+btrs4p1t5orAv7sCeLwDo0Z3uH41I+V54bhbaIucVBazVFFslPV7QQuGglwqhG3w86Xfw8v11DeF+jJi3vy7NwtxnT7FVWRkUMn2vm9m/HGwp10OQ51HQAPn9eV83o2o2ezuvwwd2f5DziOnh3fi+uGt6XnccpyxKKnx/XkkS838PrEfuVvLMpm6qIyj5y63/c3LrhyEps3HA0HOGYdRlmuNqrr5r2/WbsZzF1UP0zSZ0CvTBazpMcu7M4j53f7YxM4HgeX9E/jkv5/fPLNaEt1VDeVibvM5QCdU+sZhdAVOUaUtrxLdSIBTg3WuJ67Ro38KEu3Zkn8cvdpESnUP8rtsDOkXSN8FZjJ9nhzO+z0LmVNJqEb3z+Ni/u2OGkHuphiyuCYX00PLmh9CtrGyKHRFWXuokqMcxyXGZxr8nte1uzd1UVluhDNWZi3rx7AsCf1FdnLmh9syb2ns/OwPrVHdScBjqg20v5AbYyouWryga5asZyyH98DsLkEpzpOxX+ixdpLYP7kmKfaKCvAaVY/vtwJOKuL6tdXIYQQ4g9xcnyHHJu7MBKqaJbhmqQmZHA6Nql4F6I5Nq5jngj0D4zQrY4kgyOEEDEiXom+cOwf5XLYmH7NALx+lfoJ1W+ekxOtBsQ3XDawJUcLfAxrHzkkvaR2jetyXq9mJNd1WTKq1W2Kiz9KAhwhhBClKmu1+dqmJmRwHHabMcq0vBpERVF46bI+xvWrh7bm151ZjO3RtIxH1RwS4AghRE3nqgfevPK3E39Kn7T6zF594GTvRpV55PxuJ3sXjisJcIQQoqZr2hP2LD7ZexHz/jq4FXa7zVhMVFRvsVFJJIQQf0BoMsZx/VqUs2U1d9HrpNfrw/XeSSd7T2Kaw25j4uBWlvWrRPUlGRwhRK115+hOjO6WSrcKzmRbbdVvyYddX2PeghM7qaUQ1ZkEOEKIWstuU2JnUsbqX/8qxAklXVRCCBEDJL4RwkoyOEIIEQPM6yR9cO0gWifLzOCidpMARwghYoA5gzOobUOcdttJWYtNiOpCuqiEECIGmOegs1dmSWkhYpQEOEIIEQPMs+zaZCFTISTAEUIIIUTskQBHCCFigFYD1kkS4kSSAEcIIWKAKvGNEBYS4AghRAzQZCYcISwkwBFCiBggPVRCWEmAI4QQMUC6qISwkgBHCCFigkQ4QphJgCOEEDFAVU/2HghRvUiAI4QQMUCVIhwhLCTAEUKIGCDhjRBWEuAIIUQMkASOEFYS4AghRAyQeXCEsJIARwghYoBkcISwqtIA59ixY0ycOJGkpCSSkpKYOHEi2dnZpW7v8/m455576NGjB3Xq1KFZs2ZceeWVHDhwwLLdyJEjURTF8u/SSy+tyqYIIUS1JmtRCWFVpQHO5ZdfzurVq/nuu+/47rvvWL16NRMnTix1+8LCQlauXMmDDz7IypUr+eyzz9i6dSvnn39+xLbXX389GRkZxr833nijKpsihBDVmoQ3Qlg5quqJN23axHfffcevv/7KoEGDAHjrrbcYMmQIW7ZsoVOnThGPSUpKYt68eZbbXnrpJQYOHEh6ejotW7Y0bk9ISCA1NbWqdl8IIWoUmclYCKsqC3CWLl1KUlKSEdwADB48mKSkJJYsWRI1wIkmJycHRVGoX7++5fYPP/yQDz74gCZNmjBmzBgefvhh6tWrF/U5PB4PHo/HuJ6bmwvoXWI+n6+SLatZQu2L9XaWRtov7Tf/H8vUQHimv5Ltrg3tj0baH3vtr0xbqizAyczMJCUlJeL2lJQUMjMzK/QcxcXF3HvvvVx++eUkJiYat19xxRW0adOG1NRU1q9fz+TJk1mzZk1E9idk6tSpTJkyJeL2uXPnkpCQUMEW1WylvTa1hbRf2h/rDmTYCFUdzJkzx3JfbWh/WaT9sdP+wsLCCm9b6QDnkUceiRosmP3+++8AKIoScZ+maVFvL8nn83HppZeiqiqvvvqq5b7rr7/euNy9e3c6dOhA//79WblyJX379o14rsmTJzNp0iTjem5uLmlpaYwePdoSOMUin8/HvHnzGDVqFE6n82Tvzgkn7Zf215b2f5u7BrIOAjB27FigdrU/Gml/7LU/1ANTEZUOcP7+97+XO2KpdevWrF27loMHD0bcd/jwYZo0aVLm430+H+PHj2fXrl389NNP5QYhffv2xel0sm3btqgBjtvtxu12R9zudDpj5k0vT21qazTSfml/rLfffOJYsq21of1lkfbHTvsr045KBzjJyckkJyeXu92QIUPIyclh2bJlDBw4EIDffvuNnJwchg4dWurjQsHNtm3b+Pnnn2nUqFG5f2vDhg34fD6aNm1a8YYIIUQM6d+6Id+ur1j3vxC1QZXV4HTp0oWzzz6b66+/3hjCfcMNN3DuuedaCow7d+7M1KlTueiii/D7/YwbN46VK1fy9ddfEwgEjHqdhg0b4nK52LFjBx9++CFjx44lOTmZjRs3cuedd9KnTx+GDRtWVc0RQohq7cohrYh32hnSrvyTQiFqgyoLcEAf6XTrrbcyevRoAM4//3xefvllyzZbtmwhJycHgH379vHll18C0Lt3b8t2P//8MyNHjsTlcvHjjz/ywgsvkJ+fT1paGueccw4PP/wwdru9KpsjhBDVltNu4/JBLcvfUIhaokoDnIYNG/LBBx+UuY159s3WrVuXOxtnWloaCxYsOC77J4QQQojYJGtRCSGEECLmSIAjhBBCiJgjAY4QQgghYo4EOEIIIYSIORLgCCGEECLmSIAjhBBCiJgjAY4QQgghYo4EOEIIIYSIORLgCCGEECLmSIAjhBBCiJgjAY4QQgghYo4EOEIIIYSIOVW62GZ1FVrQMzc39yTvSdXz+XwUFhaSm5uL0+k82btzwkn7pf3Sfmm/tD922h86bpe3MDfU0gAnLy8P0FcmF0IIIUTNkpeXR1JSUpnbKFpFwqAYo6oqBw4coF69eiiKcrJ3p0rl5uaSlpbG3r17SUxMPNm7c8JJ+6X90n5pv7Q/dtqvaRp5eXk0a9YMm63sKptamcGx2Wy0aNHiZO/GCZWYmBgzH/A/Qtov7Zf2S/trq1hrf3mZmxApMhZCCCFEzJEARwghhBAxRwKcGOd2u3n44Ydxu90ne1dOCmm/tF/aL+2X9tfO9tfKImMhhBBCxDbJ4AghhBAi5kiAI4QQQoiYIwGOEEIIIWKOBDhCCCGEiDkS4FQzCxcu5LzzzqNZs2YoisIXX3xR6rY33ngjiqIwbdo047ajR4/yj3/8g06dOpGQkEDLli259dZbycnJsTy2devWKIpi+XfvvfdatklPT+e8886jTp06JCcnc+utt+L1eo9ncyP82fYDjBw5MqJtl156qWWbY8eOMXHiRJKSkkhKSmLixIlkZ2dbtqmJ7d+9e3dE20P/PvnkE2O7mvr+X3311RH7PXjwYMs2Ho+Hf/zjHyQnJ1OnTh3OP/989u3bZ9mmpr7/5bU/1r//FXn/Y/n7X177a/r3/3irlTMZV2cFBQX06tWLa665hosvvrjU7b744gt+++03mjVrZrn9wIEDHDhwgGeeeYauXbuyZ88ebrrpJg4cOMCnn35q2fbRRx/l+uuvN67XrVvXuBwIBDjnnHNo3LgxixYtIisri6uuugpN+//27i+kqTeMA/hX26ZCJtOcpuaK6K/TTIVlFwUGpjgMogiSJRRWF6KFeeGFSEUXdlFEd5F0VVkEQmJZCQbBVsbZLFYESkNTtMVY20irWc/vIjx03HRHf+v365yeDwz0Pe/76pfHh17OORDhypUrMUob7t/mn1VXV4ezZ8+K3yclJUmuHzp0CGNjY+jt7QUAHDt2DFarFd3d3QCUm3/16tWYmJiQjF29ehUXLlxAZWWlZFyp9a+oqMD169fF73U6neT6yZMn0d3djc7OTqSlpaGpqQkWiwWCIGDZsmUAlF3/hfL/Df0frf6Auvt/ofxK7/+YI/bHAkBdXV1h42NjY5SdnU0ul4uMRiNdunRpwX3u3LlDOp2OQqGQOBZt3f379yk+Pp7Gx8fFsVu3blFCQgL5/f7FRlmSpebftWsXNTY2zrvvmzdvCAA9e/ZMHLPb7QSA3r59S0TKzj9XYWEhHTlyRDKm1PrX1tbS3r17513z6dMn0mq11NnZKY6Nj49TfHw89fb2EpGy6x8tfyRq6n85+dXc/0upv1L7Pxb4EZXC/PjxA1arFc3NzcjLy5O1xu/3Y8WKFdBopDfs2tvbkZaWhsLCQpw/f15y+9Fut8NkMknuEOzZswdfv36FIAixCbMEcvPfuHEDK1euRF5eHk6fPi3+D/LAz2wpKSkwm83i2Pbt25GSkgKbzSbOUXL+WYIgYHBwEEePHg27psT6A8CTJ09gMBiwYcMG1NXVwePxiNcEQUAoFEJ5ebk4lpWVBZPJJKmtUusPLJw/EjX1PyAvv1r7H1hc/dXY/4vBj6gUpr29HRqNBg0NDbLme71enDt3DsePH5eMNzY2oqioCHq9HgMDA2hpaYHb7ca1a9cAAJOTk8jIyJCs0ev10Ol0mJycjE2YJZCTv6amBmvXrkVmZiZcLhdaWlrw8uVLPH78GMDPbAaDIWydwWAQsyk5/686OjqwefNm7NixQzKu1PpXVlbiwIEDMBqNcLvdaG1tRVlZGQRBQEJCAiYnJ6HT6aDX6yXrMjIyJLVVav2j5Z9Lbf0vJ7+a+3+x9Vdb/y8WH3AURBAEXL58GQ6HA3FxcVHnBwIBVFVVYcuWLWhra5NcO3XqlPh1QUEB9Ho99u/fL57qAUT8GUQk62f/DnLz//pc2WQyYf369SgpKYHD4UBRUREAedmUmn/W9PQ0bt68idbW1rBrSqw/ABw8eFD82mQyoaSkBEajET09Pdi3b9+865ZSW6XnV1v/A/Lyq7X/gcXVX439v1j8iEpBnj59Co/Hg9zcXGg0Gmg0GoyMjKCpqQlr1qyRzA0Gg6ioqMDy5cvR1dUFrVa74N6zb+IPDw8DADIzM8NO6j6fD6FQKOxk/19ZTP5fFRUVQavVYmhoCMDPbB8+fAib9/HjRzGbGvLfvXsXU1NTOHz4cNS9lVD/SFatWgWj0Sip7bdv3+Dz+STzPB6PpLZKrH8kc/PPUmP/RzJf/l+ppf8jWSj/39D/0fABR0GsVitevXqFwcFB8ZOVlYXm5mY8fPhQnBcIBFBeXg6dTod79+4hMTEx6t5OpxPAz4YBgNLSUrhcLskb+Y8ePUJCQgKKi4tjnEweufnnev36NUKhkCSb3+/HwMCAOOf58+fw+/3irVw15O/o6EB1dTXS09Oj7q2E+kfi9Xrx/v178fcuLi6GVqsVH0cAwMTEBFwul6S2Sqx/JHPzA+rt/0gi5Z9LLf0fyUL5/4b+j+p/ermZzSMYDJLT6SSn00kA6OLFi+R0OmlkZCTi/LlvwwcCATKbzZSfn0/Dw8M0MTEhfmZmZoiIyGazifu+e/eObt++TVlZWVRdXS3uMzMzQyaTiXbv3k0Oh4P6+vooJyeH6uvr/+j8w8PDdObMGXrx4gW53W7q6emhTZs20bZt28T8REQVFRVUUFBAdrud7HY75efnk8ViEa8rNf+soaEhiouLowcPHoRdU2r9g8EgNTU1kc1mI7fbTf39/VRaWkrZ2dkUCATEPU6cOEE5OTnU19dHDoeDysrKaOvWrYqvv5z8au5/OfnV3P9y//6JlNv/scYHnD9Mf38/AQj71NbWRpw/9x+4+dYDILfbTUREgiCQ2WymlJQUSkxMpI0bN1JbWxt9/vxZsvfIyAhVVVVRUlISpaamUn19PX358uU3JV/495ebf3R0lHbu3Empqamk0+lo3bp11NDQQF6vV7LO6/VSTU0NJScnU3JyMtXU1JDP55PMUWL+WS0tLZSTk0Pfv38Pu6bU+k9NTVF5eTmlp6eTVqul3Nxcqq2tpdHRUcke09PTVF9fT6mpqZSUlEQWiyVsjhLrLye/mvtfTn4197/cv38i5fZ/rMUREcXyjhBjjDHG2P+N38FhjDHGmOrwAYcxxhhjqsMHHMYYY4ypDh9wGGOMMaY6fMBhjDHGmOrwAYcxxhhjqsMHHMYYY4ypDh9wGGOMMaY6fMBhjDHGmOrwAYcxxhhjqsMHHMYYY4ypDh9wGGOMMaY6/wDhIFjEq2lFkQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# reduce the resolution of a spectrum by half\n", "\n", "REDUCTION = 0.5\n", "\n", "spec = cube.get_spectrum(1000,1000)\n", "spec.plot()\n", "\n", "# inverse transform\n", "zp_spec = np.concatenate([spec.data, spec.data[::-1]])\n", "spec_ifft = scipy.fft.ifft(zp_spec)\n", "spec_ifft = np.roll(spec_ifft, spec.data.size)\n", "\n", "# apply window to reduce resolution\n", "w = np.zeros(spec_ifft.size, dtype=float)\n", "ncut = int(REDUCTION / 2 * spec_ifft.size) # remove half data (2*0.25, on each side)\n", "\n", "w[spec.params.zpd_index + ncut:-(ncut + spec.params.zpd_index)] = 1\n", "spec_ifft *= w\n", "spec_ifft = np.roll(spec_ifft, -spec.data.size)\n", "\n", "# forward transform\n", "spec.data = scipy.fft.fft(spec_ifft)[:spec.dimx]\n", "\n", "spec.plot()\n", "pl.grid()" ] }, { "cell_type": "code", "execution_count": null, "id": "ab17aad2-63df-4946-8403-1f735ed84775", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10.9" } }, "nbformat": 4, "nbformat_minor": 5 }