{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Units and Labels\n", "\n", "prysm tracks units using [astropy.units](https://docs.astropy.org/en/stable/units/), but purposfully avoids `Quantity` objects for arrays, since these are hard-bound to `numpy`, and prysm's backend system allows api-compatible replacements such as [cupy](https://docs-cupy.chainer.org/en/stable/) to be used. The label system used by prysm allows flexibility override of its default choices for axis labels when [plotting](./plotting.html).\n", "\n", "We begin by performing some imports:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%load_ext autoreload\n", "%autoreload 2\n", "\n", "from prysm import Labels, NollZernike, config\n", "from prysm import wavelengths as wvl\n", "\n", "from astropy import units as u\n", "\n", "from matplotlib import pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To include a wavelength, we need a unit that describes one wavelength of light. To facilitate this, prysm includes a `mkwvl` function to make wavelengths of arbitrary size as well as a `wavelengths` sub-module populated with common laser wavelengths:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "built-in wavelengths are:\n", "['CO2', 'NdYAP', 'NdYAG', 'InGaAs', 'Ruby', 'HeNe', 'Cu', 'XeF', 'XeCl', 'KrF', 'KrCl', 'ArF']\n" ] } ], "source": [ "print('built-in wavelengths are:')\n", "print(wvl.__all__[:-1]) # last element is mkwvl\n", "\n", "telecon_wvl = wvl.mkwvl(1550, u.nm) # == mkwvl(1.55)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVwAAAEVCAYAAACsQV2IAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9fbAty1Uf9lvdc8657+kjiRFVUPqwIBA+YoMRT0joCUkgZJ4lRUICIfERZCClsh0IDiEOhAqk5KLKNk5syiYVvwAukShgY0OhCkKyJPQiISGhJ/FlkKnIEMMTJEg4GIR4957da+WP1at7dc/M/jhnn332vXdW1Tl7ZvbsmZ6Z7t/8+rdWryYRwWKLLbbYYldv4boLsNhiiy12t9gCuIsttthiB7IFcBdbbLHFDmQL4C622GKLHcgWwF1sscUWO5AtgLvYYostdiBbAHexxRbbixHRDxPR7xPRv9zDsb6YiH7J/T1KRF++j3Jep9ESh7vYYovtw4joOQA+BuBHROTP7fG4fwbAhwA8SUQ+vq/jXoctDHexxRbbi4nIOwD8W7+NiP5DInoTEb2fiN5JRJ95gUN/JYCfud3BFlgAd7HFFrtaexDAt4jI5wP4dgD/0wWO8SoAP7rXUl2TDdddgMUWW+zONCJ6LIBnAfhxIrLNZ/m7lwN47cTPPiwiX+aO8ckA/jyAN19taQ9jC+AutthiV2UBwB+KyF/ovxCRnwDwE1sc46sA/KSInO+7cNdhi6Sw2GKLXYmJyB8B+C0iegUAkNrn7niYr8YdIicAC+AutthiezIi+lEAPw/gM4joESL6JgBfC+CbiOiXAfwagJfucLynAngygP9z/6W9HlvCwhZbbLHFDmQLw11sscUWO5Dd8U6zJzzhCfLUpz71uoux2GJ3tL3//e//qIh84mWO8WUPfI78wUc/tuX5fuvNIvLAZc53HXbHA+5Tn/pUPPzww9ddjMUWu6ONiP7NZY/xBx/9GN778N/cat+Bvu4Jlz3fddgdD7iLLbbY7WECAXO67mJcqS2Au9hiix2HiYD55nWX4kptAdzFFlvsKEwgYFlddzGu1BbAXWyxxY7EBLIA7mKLLbbYIWwB3MUWW2yxw5gIhBfAXWyxC9vvvfpzACaIBIgQhDVrlIh+cv4UHo/BocBlOZCOiCT7DAIiAREDQfDJr/uVK72OxQ5kC8NdbLF5+/DXPg3CBBFCSlFBNQOrLj8eIgFc1ivgsrTgO2cGsoGkAm5eDkFB90Nffr9uK0AsiDGVbU98/Qeu8C4sth9jSPrT6y7EldoCuItdyj7+sXtboBVC4sxmO6Bl2wd1Gagsd8qM2QIKsoEEBKnLHfASCWLgst2Ad7HjN5FFw11ssWIffNHzJkD1MUgpFEBlCWD3vQFrEr+ux+MMvrY8ZUVKgDhZobLdmEG2rAdGIEEgrtsi45f/4gtGYPxZP/3Qld+zxXYxARYNd7HF1G7eOlsLrolDAdbEuk8BWijQGgAzOvDFDOCiBdnQsFuU9ZgBNoa6buA6B8KLHZnJAriL3YX2gedrTpCxFHCKxBVkmQMYGWi5gmuy75GBF7qtAHU+tkHenLTQO8pCXg7wDFdBOAb7ZETyoJsBF4IQuIBvTIz3PvfFI2kCAJ72tjdd9S1ebM4WSWGxO9l+4XkvanRSBb040loNaDmzWANYY7G2XrYDmfFSWZZ8fA+4FWynyxcyBhvQNoALZDAFojHZhAq0/tMx4EiMYMAbnCSR2fLUPfmCh356z3d+sd5IGLR69LqLcaW2AO5dbsZIATR6qrHRIhdIaJjsFMgmIawcwxWgrOvxUABYBFlWyOWYKZ9x3lDkAxSADQQQB0QSDCSIHDKzDRhY2e8IfJkb5hu5lR0MfFlUN17skLZICovdgfauL3pp6T77HPRePjA220sGPdCeGxAb0DagS5XhCiEVDRcjWWETw/XslqgCbszbYmGvCr4rMsANiMI4IUIUQSRCDLotcl4ODOZQWK+XGcze/ZyXlHt0/zt/an8PYzFnAlokhcXudPMaqoxYrWO3+btzkxUyuK64ZbfnTA6AkWUHPc9KkKWFMcsF9LvefVbZLSmrBTAQZTZLKg0UQBWcOOCNJBgASP5+CAQRhgTdJmCIADHoW0Ay22XUV5EH3sWu0ATAkp5xsTvB3nH/yxDzyC0PINuC7bnTbo3VnhsAO0Z77oDWGO5KkAFXJYYCuKjMdiPDJQLldQVcBeBICr7KcIFIhCSCExKkDLwshBQYJyTKrCkDPWV+HQAwAyEggpE4VNDtwJZI8K4v0nkQEwc8510/edlHs1gxAS2SwuGMiB4A8P0AIoAfFJG/NbPfVwL4cQBPF5FlOoctzDuHzHxUQB8zO2K2E2C74lAY7kqAFVeWq8st0KYMsitGA7rAFk4zOLANKlFEUkaaSIF2EMIQBEwKqsZsEViHDgcGODQAG0RG2yIYLKQM192vcu8yiMewhJbt12RhuIcyIooAfgDACwA8AuB9RPQGEfn1br/HAfgvALz38KW8feztX/gV6hQqsagycgIFQhMp4ONjy1+v4Tpm24PtudNtzwVIXIFWpQUU0K2SggHvNOISKEsKFosLpAREAhIBiSrwSoDmbSDkZQDBAJVBDnRDBtiAAIoMEgKBQAIQEUhoBLb+/kUSCBHe/ZyX1BhkDvjin//ne3yKd5mJgFb7SUBORP8+gB8E8OegVe4bReTn93LwS9jRAC6ALwDwIRH5TQAgoh+DzmH/691+fxPA3wHw7Yct3vHbm57+SgyRMYSE0yGNRmKZFWYrwFgxrVYGMKDG0qbCgqtcUJjtBNjq9xVoE7cg64F3ykpkAgEkFXwlaDeIJYMrSFfKMkAOdCkz1iCElD+J6oAMAa27FSVkDMDoXopQiXR46Fkvx4ojVinggff9kzVPa7GRyV4Z7vcDeJOIfCURnQK4d18HvowdE+A+EcDvuPVHADzD70BEnwfgySLyfxDRLOAS0WsAvAYAnvKUp1xBUY/TzoYVhphy8H9lt0DVbUVU77QIgQApYWHeJIMs0IaJ1YgDBVhjxR6Ap8C2SgsKtB549XySz6XrJiUwCJDqOIshb2OASZlu9XAp6FJepvw7yWUKwYBW2msSApGFxY2Z9lTkQvnM9y6Kar9EFnoWL/wc72ajPQAuET0ewHMA/GUAEJFbAG5d+sB7sGMC3Cl+UWo/EQUAfw/5Jq4zEXkQwIMAcN99993xLua3PvMVGELC2cmqAG0IXBK6eLOEMkBmilJHWRHJCG/EhXGVUK4MvAq+9c/WTUZg6b4XIIk0wMvdKTmP9bWYYOvGCyrQKtiRC2hT5SDlX5BIBVqpcoQCfAZac6CJSSpj51j5w3SinP6+BmZEDkgUGrb7pe/58W0f5V1uOzHcJxCR9988mNs9AHwqgI8A+MdE9LkA3g/gW0XkT/ZX1ovZMQHuIwCe7NafBOB33frjoHrMQ0QEAJ8E4A1E9JK73XF2NpzjNK4QA2MYEiiPpAKmur/q3RcJSCnk2FZR5xMy493inGXYL3zsbgYwZPnAljeALYuMJAUDXh1VpscOJFgJYQCQUEFXKAMrtefjLEMwGbCi0WXnzPI3hG7Yb4xcMpKNZRpGjNB43hQROSCQ4DwtjrVtjUR2YbgfFZH7Zr4bADwNwLeIyHuJ6PsBfAeA/24PxbyUHRPgvg/ApxPRpwD4MIBXAfga+1JE/h2AMhc9ET0E4NvvVrB9yzO+SoF2WOGek4RhSAghjVITmknWMMswXQZiNOB0I6zy/nOjrIzdTpmxUDbAsz8owK4DW5GZ85GipQddBiGIOsrYHF1AHh2mv9Hz0QhgjeVOmb/m0LHbGBkhpCbjGFClhZjvrQTGEBMSBwwcMawi3v2cl2CVIm6uTvCC9/7Tmbu3mDrN9tLzfwTAIyJijvV/BgXca7ejAVwRWRHRNwN4M7T+/rCI/BoRvRbAwyLyhust4XHYG+97FW4MK9x7eo6TqJrtkLVbLyP0wfo+P62ZRkFRDoGShrkl1G3rrI828Nv1vC3wlt+gBdtZp5nICHQtRKwsA8VnFqmew5gtKxmeNWP5AFrpwIb8hspse7lm7j5HJjAnhBARVoxViIiB8bPP/Eo8uhrwwod/bO19vWttDxquiPw/RPQ7RPQZIvIbAJ6PsfP9WuxoABcAROSNAN7YbfvumX2fd4gyHZvde3oLN04UbE9OzgvQxvzpZzwAupCvpHKCNxEd5iqJRmBjRk7iDCTah5+xuXjacj7HIqeYbZESXDklg65ptgJlthd1S9m1mbbrt/vrL4nOg8s0lu81kSDEaaekzXhhw4WHmLBKEfG8HmOxKRMQ7+3efAuA1+cIhd8E8A37OvBl7KgAd7F5e+szX4Gz4Rz3np7j9GSFk+EcJycrlREsBCwDQW82/Q2RgFPovktgrmFNiUMJxfLAG1CHujbOqg57Q9ZTNxlLC7yz+7nzcc6FsItZSFkpX/cJeKBV/dZnEVPATQ3YhsAFbDfd80FWSKuIgZXhxhUjEOOdz/5y3FydLA41b4K9hYWJyC8BmNN4r80WwL1N7N7Tm6rZnt7CEBNOTs4Rh4SQw8BCdugAUJc90EzeSCRI3IOtjpaSmJSVJWoYXknyPRUOlb9DDqkKsh3YUg69Unmggvl8HK5bJotc2Ay6Fq9bzjtzDf4a/awRZLl1Y9XFPdjafe8nswSgWgaAOCRwCggxgVNECNYbWSHe4mWk2sh2ilK4LW0B3CO2Nz39lbhxco57Tm7h3rNznJ7cUhnhZIVhWGnDD0lHSgUHKDaUFSYlBIAJxBFMglV+7DbrgcWmepZbHEb+jwRBqHwCxggBE0lDBlMDX0unWKIH3PXpwDB12pUxCo7B1uQxbQiYfRdsQETJt1DjdfV3UmaJqL8R91n/TD7onY6hMNv8OSQF3nzvETppobv3gQkxBTBH8BAQh4TV+aBgfq5J0P/0/BSPnp8sAyWgOXHvZFsA94jtMWc3cc/JLZye3sLpybmC7em5RiQY4ObGD/RJaXQb59AvTqGA1oAVVhiKnBAye9PQrpquMAZBktrdrnGpakUDlfyd03YN/Axok9uuDrMMrlmfNdDtzcDWJ69pQNZt89IBNcdw4A0fYwv3ommn7zF2WwY9BL3PHmwp5p5FZr31GaT8DKg8A04xs938zIrzLZUh2He9iQCr8+suxZXaArhHaG995iuU1Z7ewo2zmzg9vYWT03MMwwrxtEoJITJoSA27AtAwLKIIEQaFAA4RyPU5iAJszJ/GcnXer5wbVkKZJSFlQDR2GEmQSMOzCBX4KEdiEeofoDIAkWhqAyjIpw50p8zAFnZ8N+DBpAXq/wrTbcsWC+N11xGqfGL5cIMNh/a6rXu5hZAQTtwzIFn7HIIQZKWAy0NCWGU9N7NlY9Xv+qKX4k/PT+9eXVdEQ2fuYFsA94jsDU/7Wjz29CYee3YTN05v4exUwfb07BaGk3OEk4Th9BxkDT0mlRJMOyTOiQUShNVTLh27BQCcAzFozKgEQoyZ0caUma0bpkqMBHWkRWIkosICjSUOOYmM6bjGYu0zBmtHhIQsH+Qhxp7pTlmblpEKqMcsW8TQMtzgwZYsjWM7NDfk6zI5weebKHOeZXZrUkKNUujAduOzgD6HYYVgLPd8ULZ8vqq6cJYtLHTsY7fO8JIPvH7/lezIbR9De4/ZFsA9Invc2aO49/QmbpzdxNnpLZzduKkywtktxBNltyEywsmqNHB0w0wlN3JigqQICZLZrcAetzAhCBd2y5xBRUhDyzLoRmYkCkXbXImBsI5Gizle12ZgSARE5GGzOdGMxcpafKxJDBo3q0xX42fHQy2KdOHA1ksIkSqrtdSNOssDSjn9lDz2vb0oql6LmnsC4rr7PJISyDPdkxVCdBp6TMVxNnomTBDOoBsZ4ZyRYtTfkjTON58o5+6yheEudiB76FkvV832xqO4cXYTJ6fnymzPbmE4PdfGnf/UUaONmwLX0VQWc8sBsgra+Fe5W4sMtACEGSJJAVYSQgojlpukAm1kQqJQpQSgOM5iEKREBQDFsdyIPPgAHmg1pMGKzJbXdiZ5jckGdvxg7BYto43Ustvol0OehsdNOBl6oHVOQ89uNd62ZbN1OYNtTKAhD44YuEgK/rkIK+hyyr89WSGcD1keypEPLsaXSPDQs16O5737J/Zf2Y7VBAvgLnYYe+zZozjLgHtypo6y4cYt1Wwd4NLAjt1y7cJCARUcAFIgllVsowI4z2KQs3yBk+q1kQvLJScllLm+bJ0ICRpKFnMCmGTSAJATgdtsClLAtBYAmuUra7huE5CHBPv92/nMWg1XZ3eof57dGqBGsmnU3dxnecp0P5NvnV7dx94qgIaY9wlVxyXHdg08KaiOi8xW554N8Up7HqugzzGmvG8GbDcrRz/bxJ1vC8Nd7IrtHfe/DPeeZmZ741Gcnt3CyY2bGo1w4xbi2TnIwDYmdc5EruzWG6tYSnmkE0hHSQkFN4BgADLDLV7zhuUyOOYBEBl4ogR1mpkDLVVA8M6zSMpobVlItQOLSiAoUBroEnQfksqMe2uAtjBWD7LI0+wo2NZJJsfOsgq61VEWQ3WWWYKa6MLAShRIqEzXyzrNS3DIL0KS0fMhaM9DpZ4EiVF/NzDCUAdSUAF4/Xz4i1+Ij986uyum8iER0BKlsNhV2BvvexUee3YTjzl7FPfcUHZ7duMmTu65iZOzW4hntxBOzwvgGtDSkIChDeAvVpIaELBqB76SqBOHYtJk3EIKJB3LrSOrAlha51mkgJTZ4UCUp7fJ84ZRPbVklmvOsggUKit5mcQ2VeCdsqkZHwxsB+cwM+eZsVr/F5pPY7hjZ5k5xUxLHbPbHBWSmakBpn8R0pB13DDzfADQiiCRIUmPIzEBVMPLmpjq/CJ4x/0vw8dunt35ORgWhrvYVdhjz24WZnt2dhM37n0Uw9ktBdsbNxHPzhHOVEowwMUgmd2iHUYFaJ5Y1pEGEjQewFtg0iG87EaTxcpyKVVnkckKIcsK5jwzKWElAZHqVORMAg5ASoSBdKobi0pgcpljgiYkN8DUmRcUKGuy87yrkxMAx1od2FZJoUoJMehsvcZkDXS9E61OPVQ13NDJCcHmgHPstsTv5miEUIC3fRlSZrj9cyLkF2KE7mfAa85PY8e2v5ck1uSvuGNs0XAXuwp7x/0vK8z2xj2PKrO9cRPDjZsYbtxCODtXhpvZbQHaATqbQQ+4XD9lJSAWCEIFXQ7AkBRs7U8IlMEWnBBjyMlW8jThLkRMgZaRmPL04xoeNpAUlhuKlkvFWQZYfG0LuiYvALpZUHEmdrhiAGvfeSfZCGy77T5m2MLDBuIcsVBZrQ8FM6cZBRvK28XaRs9wcyiYB9v8QqTBPaP+WTEq9R8UbOOQVIIITgM2lhsq233H/S+7g+WFRcNd7ArMM9sebOONm8psTzOzPWFtyAP0aflhU2YRKHN3ARAGCKwRr6xAIUzZsaNs1tgbJwstkxIGRTkUbOxEUy03Zl1XJlgus/acNWKBMBQAraBrYVxFYoAD4M78rL1DqL/1IDu40C8bLTbFbiPVobv6PdchzI7VVqeZuzc+fCtLCTAQ9mBrL8UNz4okM+CVHlfOQ5kluJ9NAshON9zhTNfE/jvYFsA9oD30rJfjsU6zvXHvowVsT+69iXCqrDacnYNOFWzpFPqUCr3LLbhM+mXxVXmSMIjG3ENBF0NeHpImmmFSdsYElMQqWVZgS9KSnWdBHLut0QqJAgZinYo8a7llyhqiHJ0mdebcAqdaqsJ5AxqH2pSVobxoPz3Yeimh2ebYrerOVUbwiWpizm/rQ8GKnkoWgseV2Wa9lYYEmAPTwNaY7TbPKwpoEM3KbnkZZrRsMxHCw1/8Qnzs5o07MGRMgNXqugtxpbYA7gHsDU/7Wjzu7FE85uymAu0NlRGGs1tVRji9hXjPLWW1NxzYFgoXagP27vzSgAlYmW4rCmZZetCuqjJS0x8pJ12RZJ85DKxznpHF3wbdrtOS24gzOHarSVskO9Ay3COP5UVJbkM2gaQDXkxHKABjoLWJJAkt2LYSQwbZwI7looR+lTA3F53QO8uiCwuze+XvHcq9ZFDMTjKTegxshzgGXHtmTDmhBDehZAgrdRCGFnhFCGdsOrcmjX/7F34F/vjmjTtnRNrCcBfbhz32tDrI7rHQrxyN4GWEEdieWKPtG6/rpzJn3Svo01ylIi8gABRz9ipr1EMCUqxgkdmcaYdNl9o5zyy/QhSNXoghQERZbiLVbQvTLaFfAEAN6CoIa5faRykAVV4ws+0GtD5KwSSECrbi1lsJoWe3PhQsENcBDj5TWHFk1XtUtpmMUJxjzkFGHdiWl2T3zBD1c5WjGlbJAugAKMsL7EG65sdgDjXR+Z0mMSzZwha7jL31ma/Q3Aj57yTH2Vo0gkkI4ex8GmztzzdaCrVihiycrlYA8uiyQYXUhuVmSUGYQMQQYoCCgkfOKOa70oGti93mVzB2G4khASVigXPyGxsEgRKpAHjQNRVXQ8Eq8JpZtEJw8FtjcV1cLbAWbMmcZvn3Pbv1eRPIBkAUJtumaCzgSjnncNFyUwHawm4HB7CBgGFowZaCMlnA5jhyPZMWdO1XNmUPoFquzckmQhAQ3vrMV9whCW9kYbiLXc7uObmVE9Foxq/TUx2qG890UIPG265azbYH28E8MXCN1yFagO6zWuWua3DsS73bllTRd4sLow1VVgiZpdrMv5SniSEPuiVtI4qWqw4zBZOSrczPzZNBl/KU5CGz3aw+lL1iP125XS0ZeOYBDiSdhICxlNCxWxvo4JOMl/AvG+lFUkHYxcT2983CKtSxBsduXXCwge3cs4sOfIs50OWV5h0WKgHOQwqQFCAccJaizr7cJZa/bW3PkgIRRQAPA/iwiLx4bwe+hC2Ae0X2pqe/Eo85u6m67alGI5xmCUE1W8dsz843g23fLfXG+V/IDJcz85JMhq3baw4ZH/cJNLKCj8k1YDLnWcq5cpMo0+Uc88sgMEsZvNB+VtAN0M6yJbVJLpph3YwPFh5lQKv6reVC6MC2cZRl1puB16ISors2C4PzkoLF3hrTBVDkhDKrQ3MfUUEXGPdGpqQg65kw9Bk3lkGXM+gi58dgwlDkhFCkhZQC3vnsL8ef3Dy7/ZOY75fhfiuADwJ4/D4PehlbAPeKzGZqODvVvAgl61cePRbzoAZEdbpUB9kE2FqDnAJc5va7wJVh5aiFmgDWxvhzkRUUaFVT9LICBY3JZc4ZrLLzzELEJJCORCNL3aiAW9Mu6qeagm4OWsu6rWTg1SnNLbLNmx/0YEBbHGeEHFdbWS7ByQ0FfFvtNjYsViqrLey1ygl2PzVBkMt1kF9ODciWJLwuOqEH2KnnF4L2THrQ5QQ6VdAV1qHEynIDODPc0xTKZJUidNszXRGNI9+HEdGTALwIwPcC+La9HHQPtgDuFVmZqcHy2easX/E0O8fs7zQP1R2Q9T8ag+2cs6xZL/9a69ogBSnRCrqhlRZMz2SqQFvzyarzzOeXNS2XWXVddZgFcODi6MknApewMO8wkzIJMHUOIAuR8kDr0y2OhvJ2UkJwZexn5PUjyxrgtWV3T4p80EUnzN1jd7OdnjvxDM1MDjLQZda6IKLxvZIAG7CSIuJpADggrSJOsxMtccA9tzngFvl6O3sCET3s1h8UkQfd+t8H8DcAPG4/hduPLYC7Z3vrM1+Be09v5jnIznWmhpNzRJdeMVhuhDIkFLm/HKqDjCYYkjfrkk5ZCEA3tFe1Roa4WXt9wuyp8fs+ZSEHQeC8bB7yzHKFCEwMpogIBXTNn9CCLhGK40w64AVyiFk3tBdogdacYSXl4iTYmtzAheUau/Ujy7yc4P/KPej0294mc9auk37WPcMQcnRJrgMsNU6XWUcKpoBwsoKsAiRFDKfnABNOUkRKEcwR737OS/DxW2e3rxNte8D9qIhMzspLRC8G8Psi8n4iet6eSrYXWwB3z3Y2nOvsujbh47BCOEltisWcR7UyW0x3RTt2JEET0qzNih9C1nAd6G4gPl5KQNe17mUFgozYoo3c0nnUQpuIxoEu5bIEUbbrgdesH9prxS9RB2SOMwe4GWQHt23IINvKB3kZE3JCkVLacLCN5u/tHNB21jzHqV5L43jLIXVJQKuaDlJSRDwfwCcJw2qFk5NzZbrpNnei7UdRuB/AS4johQBuAHg8Ef1vIvJ1ezn6JWwB3D3ZG+97Fe49vYV7cySCTfgYT/WzpPQbuGW3jXd7DaPd1owxzbHfzNiahCn+a8/yguTRZ6GNUbVuORRoWbhMuW7SQoNEGXRLmJNdcgbeKf22/NQ+qYIuZXAloAFbY7Ml5AuV3YbuRdHG3dbctY1+CzQgPJVgZnzv41ju2dZCAJJjvE7CQKwsl1IswDvkWZkHDkVaSBzw0LNejo/fOr29sotJHcJ8qcOIfCeA7wSAzHC//RjAFlgAd292Y1jhxsk5zk5vYYgpT2WeNCmJ5U0tiU8yu/Xebc9wabqxNsy2B1QfMN5/t6GbZo4gcuetOm4Fl5BlheAAq8z6GxjMUcE2BGXFcBUscJ3FFjV2ocy+OxPAb+dvWW6bI2Eoy94xhm4kGZd9SsxtD8AeaLMEs9H6XZiBzGAh3YvHjr2uh1J0X3N+ov4NleXSoC9wTjr1+jAk8MkKw2rA6ckKnKMYbju7s8c9bOpsLraNveUZX4Wzk3OcxBWGmKqUMKya6VmKFmh6qfdyr7MRgDoWW0aaTex3ATPHWa/nFmcS1bSNnuVSxyb9AAQfrjUK29rhb+44VM6FUp6Y5wbz7LbVa+u9GjnPJvTanW3qGU09n3UvTrMm5MzipuvMEyEmhGGFYVBp4WTQunh2co63POOrLn8thzKB6tNb/G19SJGHjiUGF1gY7l7sLFfwk5NznJys2nmqHNj6sKy1Jty+6accZL4BiwddqX+WkBwuy9QOQ0ELuyw6rvfytyw3BoGwQHI2sTLcDKokD9AY3JiPK5Jn47XLmSlDG6kgY1mhAV8uUoIf5FDia2lKVrgAuLqRXiXHrQ1ftnsfkNluvg/rXobrANlZcXx6rT22dS0OCSe8KtLC2XA7cWJ62PIAACAASURBVKo8wOMOtgVw92CngzLboeRPrZ9l2hVzxMw18AKUfbjXTCPsWe1MoxVGBV0bATZTqUvZvOMshKLjGkCJSHGeSWa5krvrTARCnkrdeUAMdFs5QcqLYK6ZFcB1QGuyggdbP217Ix3AAS0qyPaMttFq/b3oTJjqdDlI028K5goc9rJcF1Uy+TKdH+ZawtLsOQ0JYdXWPauPw+027fidlhuiswVwL2FvfeYrcDac456TpBJCUIYRA+fptF0O1TwGf2TZz1LXuTKjdYSHuXY/V6u2kdoIJv/7kda4XjOd2q66LoOIQB1TVDAzhxkcwNSXh4FsEGqAlwE3SKK1Oq6gBVoDXz91Tgn9yhpu6xhrQdaOs8s9AFP3rADEnH/YLsiumfOVWXxt4lltftRT2UaD9y8LkjKpZQwMHhLCijEMK5yyjkK7uTo5/lCxPTnNjtkWwL2EncYVTqI6LIaYyhTXNgZ/o7H7LKzIjzwyZtT9zoC20QelXS7zm6GCK+fJJaWub7JRbKrJDFQZbgwMSZoXV0jHlInYgIYKujqljk67EwDLXtOw28J4HejVSIV5sO3Z7OT6xHX0UQmzZmBbEgC5e8idrAD3LA10Q9DVKeD1ktBq1T1LO//mIpY0kitNOjRwAA8JJ5xun1Cx26WcF7QFcC9gb3r6K3E2rHB2kjCE3JXLMwZES1BtbCoY/UFprOCg2bz6hrpK+YmE2kgBTLa2RkrIDXSVUJIUcJUTxCU/QScrrGMUgaScWXVcY41VVvAhVpyBiIqO68ufz+sibwt4W/RCBmB/fvj9YEw3O8gc2JbZHFBlhYbddnKCXVN/rikTJlCscgI4aJgW5/njWCrLNRnFP0u4ZxVmui79i9OepXtxilDJqVCH5/XZzbQOphBhM3gMIeF00Py5N1fD8eZbEGoG5tyJtgDuBWyIrIyWWEO/chKUMlqrH/opASJel815Doh0DjLdSb+z/KiBipMNQAVf390sTMh/QhvnCpXdctVvgbGWO6XtNrGoHfBRjrM18BOn5Wp5gxI5EtcDZwQQVggIkuNLkfO5ZmCfkhVaSUG/t8gDkzmiA1VLVtOz27bs42vqlw1k+3XhoGF9GfyIk8a2sUBW0FGDlJ/n6FlSfX7985x6lkmAlT1Lansr0LSNkDFAmYxlI+lC0IlAh5iQJvY/KlsY7mK9DSG5mV9d8PwEYCiYpbq80uTfs+yIoflsYSPGSqBqPahvpLbumG35W1Flt1lOEFuWDmRnKnoPssZ6DexEOq+/oEgLACMhFNBlAAPUscZCBXg5Tzw5Zwa0xmqnwNbnTggYs1sPumEGfKesJG+PCnD6vPQeElF9jivUUWEQN2Jq6lnm9blnuarMdqqXIqs4+QI1R5pdk9XJktNYGEM4Xiea3AUa7lG9TojoASL6DSL6EBF9x8T330ZEv05Ev0JEbyOiP3voMr79C78CJ5nd9hMPtgwpFD3SGmhpMKyNBqvchTIWs0JmNQLcSsqQtvpzMkJhRLmRrvK5TE6QCrrGjvysAetmEPCs0HfNzdo4VzTs0lIlGlOtQ29rPoSA9q+dHkd/008S2YZ7tXHD/UAGX+a1ICsTLyEJrQZe7mVwL7b2/pdnue1zvOWepasTWkcMaA10XXlQ5QZ/vV7uCTmhz0lMePsXfsXstV+vuXay6e82taNhuDlZ8A8AeAGARwC8j4jeICK/7nb7RQD3icjHieivAvg7AF55yHL2uVS17F3j5RbAhAmSoiaAARqdkpBKOu4yO4PFcZoLH0CfjCYfPJ8PDRtqAMCAYRVL4mqRUNnuDNjOAW9/rV5WSHnWhyA5VSNUz7VMkUH0+kQkM1zTiU2PnDylY7gVWD3LJafp2vby2+5FuM5JZrHB/vpFdEp5YSpDsBtgQ46+4PwcWYAhywshX7SdbtOz7J8jd2Cbgn6uNPG4JAXg5ll1YOT1ap/T+FjtjpsyqLOjAVwAXwDgQyLymwBARD8G4KUACuCKyNvd/u8BcPDx0Zb6z4cX1fJZI1V2yymAKGoXOwiwCjq1TWo70ISUtUHVf+vkjwAsn+2UuWgiH5FQmC0HIFVwFQ6FGaGAbtjYaHur3VUayQoBKilYJjDk9AAedJlIQZRy3lzIRoCvMsAYbGM3TLcZ3NDkUNgCaFzol3XjhbKsEELWakOeBl2loXw1oMCQIcsLo+eIVhZqzpnP1z/HlM/vwdZAdhXKi1xWUXPkipOK3P2s4Xyh5DI+ShMsGu4B7YkAfsetPwLgGWv2/yYAP3OlJZqwTWFENrEf53nCRFiXs5NEFb1uKmih7PXOGlzIgwJCB7Z+uQ8X4spqwR3YFkmBRuxWcvfYg+9F7skUy2XKRD1HTjVMl6DaLYkS9XxPbb6uVg6wy++Yrh/U0LHb3lm2i0kGWSEBUYAEGbFcyhckq6hJZQx0V7mnIpT3ES11AV13Iu6WNzxDBdcAPh90OYUKtOyWN8hCF7knhzFtN3eyHRPgTtWSyZpBRF8H4D4Az535/jUAXgMAT3nKU/ZSuHfc/7LsKBs3ZD+DKosmigZTAdoAVDkBOUBBWGdlFY1f1QkJM5MakiYvdIA76ZArWqMdODfSxrlSlznFSXY7Z7ym8frwMJMlepabfYPKhCEFdI3ZGuFjqueYjlRogRaAY65owNaz276s/vf9fWSh2Y6ETLBcRtQwtSEp6AbJI9CoRKlIAJDyWydt8QxNj80vSP+ytOdnuXAVbBWEOQOvabopxclZfb22/a4veikSBzznXT85c9XXYAvDPag9AuDJbv1JAH6334mIvhTAdwF4rojcnDpQzvz+IADcd999e3mdF4dQ9w6whioSwElAZLlOYwZakw8GCBOiY7dsWfwHboHXx6Na450ROIueaMDpYn1Ns1Xwj7mbGsfs1jVO65Zuis9NXZdVlKaW72qMLYqeq8CZwZbqe2IbScHOY0Cr37W67Uh+wDTAzlkBVpHmfhB3LBcRhFRAVy/Swr0YJFLvXxkqvOEZds+vgm6oIGvPLwWkzHK5AG8Ac0TiKmeZvMDdvbWY5GPUcu/0KIVjAtz3Afh0IvoUAB8G8CoAX+N3IKLPA/CPADwgIr9/6AL2eq1qtdo4LejfpAROed9z93vk0VbMCrRB9FNsviypwGsxuAluAAQah1vj8LIBFbbsG6t52EujDY1227DdXs/tTGUSGq0byywhY1BArd+hOMeM7UaqswnPNX0PnrZeJAYHsj277b/rr2HOJPdQhKqsoIXPYFxKW0GXmGtLygxXC6jgC2Tn17pnBxRpx/dORMaAyyk0YMvnOtuDzXVm0gJz9xLtXpJpzXO+DhMsTrODmYisiOibAbwZ6rb4YRH5NSJ6LYCHReQNAL4PwGMB/HgeWvnbIvKSqyzXu77opQDaRlomSswgOxo22Um0OM/sibOMwAzKQKTzVEnNKJZiBV9jt7573DlexDEjXa9dS2lCmUIDttoFdaDrGiYwDjOaMnVCVQcgobJczmWy2FzTcwE4tlvDTz34tsdT61mtLre6bS85jMs6b8KhOPvK9ZdeAEAU3Hg5vYqabFzvmz03sYEvqc5+Mfvc8n2ae24Nsy3SQgbaPGW6gW1KEWkV88wPsSQkN/DV+92+LK1+3//On1p7fw5iUknDnWpHA7gAICJvBPDGbtt3u+UvPXSZ1oUQ6aey3CnQFSEESYiBEIQhknIjTggGuB5482g18YDb5WeVKU+3aX9cyzQCWsd6C9gaqDh2a8fhcn0t451iIERVx/XMqRkanEFXxAC5st1yGd2lTTrOJsC212ZHicQ7M9btr4+yrkyZyRJzYbmWPJ24BV1hAQXKIV0GvKK/RWiyr0kfCeZDB+dejh3QmlTAHctNXLdNga2x29lnd0S2SAqLFTPtsaRABZCyV9UzKBFCjA5gAyNmkAUnBbSUk0ingBCVDWnSGwYoNOALoGm804VrQXcKaE0+gNf5Gu922zBboF3PPDzQeknBZw4w0AVatluOMWKl7fHt2LZvLx/0zsyd9FvJaRfL9auWyynUa4lcQNdelBAeAa+BraR80VODu7Z4XvZiLBJBBlUwNayW8z7MESnvm/KUOyKhbCuTf4JGuu5RmOgL5k62BXBn7Bee96La8JxpRdUuscWfihBCMMBKxRnBTIiSIIEK2w0UEFiTinAGWQNfSRlEBs06pvkYIsq8Vhh3Tb2NEtLMAG3Vbh3AOj1XjzHOrzBlIQhSqkCrmDnPpDRPD9X7RjUcjCZONcVWeyfZ6Fz94Ie53LYSQKQvQIqi1xxTjjSAst1g8bih6NMkpPG4gUfACwe2NAe26J6Vk288m5Xs7LSXowdadgBsrFalg8x4OYJZ17cB2vc+98UgEnzBQz89XeAD2cJw71LzjG2qK1biTWEeekBEdV2JymZDBhVjuwawBXgjQ8fYS/mOAutEgSURTs15CmDjbBHNGHsHtH69j9ts4nM9E+quf1uHhkkJDcv17LeEg9VyrmOjHmj9up2rZ7nbmF5LdfjZtZXwMEaRFrgk6oE6MvNFEQcIowVeoEQsjGQEf/7umUy9DP3zmQJaBVhltQasHnyN6RrQbju45LpsX04zInoygB8B8EnQp/WgiHz/pQ+8B1sAd0ubqrBEkoPjNQzImJU14piBlzkghIAYGcyMYOCaBDHmxhwZknTgQ01cXid29KC73eSGXYMuy23DNrBl3+107LYAkWMeu3ZHp6IWRJzjbMPhvNQwFZWw6dzezOEJ5B5INAaLwnI5sJY3BYQI1WSDrhfNXdy8b8Z4gRxZssHxs+HZlJha91ySY7peKpgC2qLdOqCdq7/XDbKNyd4Y7grAfyUiHyCixwF4PxG9pUsTcC22AO4aE7SVta+0RKJRW45dWTLueeDlCryByzbinH0sSJYWGDbZJJVPy27lxvKv8X4XZjrBcAurddvYNfqpyAV/TG99qBjQsVxpnWi906rXbqdsNJBhhuVuY+35q8pcrqEDXQqaEwLERUbyjJbEZR1L9QXZnHPmuax7Jomzxp6dmrsAbXJSgneANvcQtQciU3rOwa2taxc1Efk9AL+Xl/+YiD4IHcm6AO6xmlVSq7QlvMbrk4Im9pRIoxVskARLQCB2QJvU6cJ5UsO8jfN8YeyYbXRSQnHINdLCZk98WXYsqjhrNoGtZ7QudGzOfAPW7IQTvQFxIWMYA+/ccfvl3oHW7N+NOJuzytyBWGY197FrFXRraEKOXgAV4AUyw+3L6+SEbZ+FSQae6RZpIDvHfF0sDjMhnRbd1dVeFuLueViCoRjYAtyOwpG2QwLyJxDRw279wTzgqTEieiqAzwPw3ksXbg+2AO6EfeD5DwCII7BNuVJ7M4YLqaFKHnQD5eXMVAxYQ0jNNmO9Ji8wVVZbAvqb9IIVhNdaw6xq8HsjF0yA7RS7vYgVIJbWUWbrWi5aC45TkQcefC8b2uRZbhP35UDXei2UE78Y4wVQoxHWmQ/ZQ9vb6O+/jxrxUQdes/X6rAdYTwz88x1dMwkiuNRVi8D5wPMfwNPe9qZL3M1L2G6SwkdF5L51OxDRYwH8cwB/XUT+6LLF24ctgDtjpfJbRZaA5ONTHWOgbhSUOtxCSYRtlTqJZmvSdQPfVKQFY72BldWGPGuuB996nnHXtQeenjV6dlWugx0IT4BtO1qpdaj576fMgNUv96Dbl1smtvnj+e/mYm/7e6Az21RQbQcwYPS930FSzCP/dPh1sMgUUAFf+826+7/u3tt9N8mAyzLNygZJXPhX3ocbCQxr66oQgBAQwSM2fF0me5IUAICITqBg+3oR+Ym9HHQPthFwiejPbHEcFpE/3EN5jsIaHYxDAdvkum2jLlqXqYosbaBkBssujy4JWLiRG0KQwnoD1RkLgg1+cCzXABhooxbWaZjcNP42hV+jI9p6B7aXsV7jnQoJ89/1tm4gw+XZbZUWPOhWRpsBKsWi2Zd7zhe75+tebD6Mqw/tYgkFYHvg1W00ikiYq6cCG6wRAAY4cPFZXKftw2lGOgz1hwB8UET+x0sfcI+2DcP93fy37k5EAPtJy3UEVphFBlYD21RYxDhqoe3iasVOUmeVTVlaMIA11mtyQ8wRCTGy5i61rh4rmy3gyxV8AZS0j+S00U2arl9vvOQFcMfM1bPbuXvmjTAeqtuzWvuNB62p9Iz9PZ79HmOGOXcvjMR6WWPEdKGLymZdCJkLzpcN93ruPk/1HvzIsCnZoGezuj9GRKBht/2LjpDnkpPs0WRQpFIfrxVwZT9hYQDuB/CfAvhVIvqlvO2/zSNZr9W2AdwPisjnrduBiH5xT+U5CuOuERjYJuvCoXbbvPlx/sYkYhAkGMtVVhcbaSEzX1Hmm5gLw2XHbnvwBVAAWM8970iasklQyCBg2zzj5W7/bbt+m6SEvoHNMcZNOu42Zlqt/W4qJaPezpgHS7AD2XmpY/pc1H2O720Psj2bbZisY7QmbRV2m0mB73mtq58BblbloFJ0JB495+uwfQCuiPwc1hPEa7NtAPcL97TPbWMt0OryuVsvWln3O5NTzXlGBGW5uYKnnHE/iSCSjqyKwg34eiebsV4iKeAbQp1lYgqA9fzbp90rKfy4B4gJzbFjuHzB7t+UfjvX0HqAu2jcKDOV++N1YgPdHvw98IrUe7oNIEzdU3/vjKX2LLd3grUg27JZBVuMCEDjLOvKFZxjl0FgGzYOQQp6juEaZ4MQ5DSid7BtBFwReXQf+9wO9sEXPQ83b51BcFoqfJLMbJ2Om4rcMDZjTHUWAmUPKwlZXsgsNwNrYkIM03JDyr8JZI40Jy04XdfOxcVxE7ZmYsAc22093FOSQ//bTeyolxE2abe9bUpK442lasReWrCsXlOgO2Up2W+3A4K5F9aclDAX0uVlA9NmrZdlgGsgmySU+OYSzjh1//I9jFT12yACkoDEAg56nl96wZfh7PQmPuunH9rqmvdm+xv4cLS2dZQCEd0HTfz9Z/PvCICIyOdcUdkObiKEZNmX2AFsXj530oJVcGCsOxJJyQamSf/zAAgSJFCWGQRRdNRTkvVyQ2G9FmPagS9Tq+vOxalORQP012/Wa4AebNf9bhfrQXfTvhcxD7bGcnvQBVyCnYnfTWn1/Xa/3r+kpkC2ccx28sKUbODlLOtlJQ/aqNr0pnrJJDrYjRkBQXtizCUPg+nH12HX7bS7atslLOz1AP5rAL+K6RfobW++QSiTzUDLoYDtuTWOKTbhmJJGKlRGYR7uGBgpVdar04VXuYGFkSiDcgbfwsLycWybzbSQcwMX5isTgGt5HeasB745jRcYyw+bbJ1uOxWp0P+2t03OsykzAJ0CXV+mpmLLOMdDf829tm2fU5qtDzEs9WwU5tXWP5MNVuIkBZiUgAZ0xZXHX0dwwBvNuZr1WwtnSxQa59sCuFdjuwDuR3IS8DvWvHZW5YPqNDuX2q3rWa63huEKIWXwVbCNymqJ6mcGXgnodF6dYdUAwUDWIh4s/MyDr8kJvZ475U1fZy3gtkDrv+/Z1EXtKsb020vKg2oPunnP2WNsmhVh6mU0KSegBdgefH3YYQ+0SdqXvNU91XxRttmVrKuTbJJCdpqFoH6GJFL8FVbOg5vQIik4+x4i+kEAbwNQ5hI7pqDiy5o2hC72lrUynwthxQErqV251UwYS51lVkfvFDmhgK2ur0gwkA6xFEJhvBL0GEyU43WlAK+xWw/A5XMEvmjYby3f5g5KLx2sA9v+Hm6yqeiEbW0brdfLAf02D7oAOuAFtum8+f2nnI09yPptU0C7idEmIazmQLcct5Z8rk4OhDILu7AmVg8G9kTgYFES1yMpCADmu9xp5uwbAHwmgBO4ZwvgtgfcR7768/GnH78HIo8psbfWCM4lRyhksFXQzXqv1MBym61AE5RTiVJQzzCqdivUgK4U4M1ThnfA66UGD7w2iq1kKytOupqnV8tjOqSXGWqlnssVOxWBMCUz6LXv3jgvAroXGeQw5TxrHGbddYYwzez6/dY6F2eA1o8GuwjQjhlur93O1EcgT+apdVGINO0nB5zn+xBFslxhZQ74v176bNz72I/jia//wM73/aJ23WFpV227AO7nisifv7KSXKOpsywWh1kfEtYyXGBVnBa+oemx6lQwAEQ9i0Ou6ClRIydwB7zmzLDGwSCVFShAyKaB4ZyVLHeNZ4AXAqfxStEjPdgAaJOHd/dk7l55m9Iwt7VdQHdXsG2cZR3orjs3dwR3zjFmx7VtPrn3tkDbMlY3mhFopCsPtCmz2SppQUF5VCa7b/ppsyQHyf8ahhuQiJ2DWB1nyWZ4PpQtkkJj7yGizz6GnJL7ttFAhyIltJXewLbKCmgYBgAFOnhZASoRSFvpLS948RqTgHNo2EACRgCzgHMYGRMj5tFOMegsETFkoJ4AXkugbqwXQMN8e4fQVDfcts/ds3X79aPM5mwb0N16sIGOoWrK5EEXwAh4+3PMsXh/DPtOCptsY2yngNYPvfXRBlW39ZEIKMveX2Dyln3P+l4tTNfK1TrM6j1mAiKhVELy53USWo2uOKzzTLA4zbw9G8Criei3oBou4Q4JCyvp8KTqax5sV7mxrJiy46xWfKv0fgJEkxUoLxuzJQJSrvSMup1D1dXEsVRbjsRgihDhHK8blfmyabo0CbwAWtabpwYyUDJArjdi/cwLc929yzaSucEPF5EQeutfJP4a+heOt7VRGzNs1n43lUhmziHWywfeQWbrPdBazypJBVirh8C4LpKoL8F6TpYcLQBYESF60HX5Gqznd0hbALfaA1dWimu2GrZjeloGWMcAzh3wGtvVCj9OEJKkzlKgwyl1nyAZbPPfkEGQGXmGWClst8yLJgQOQIRGIDCrzIAQIKIMWGWH7CDL3ueq6brpgMjyLThg7TBtl3i/WdnhgqMq9wKwHcsFWhnB2zq9cI6xrwPaUo9cnoOpIbjbslqrZ+Yv8D2rCsAZbPOlzdVFMnYLfdtT7nVFaUmFEQ5znB0aAI8NcPedvGtrwBWRf7PtvrebmZRQGW6VE1ZCOGdqGsG5AImrdtZP7w0AIWsKJAqs3AGvMQ2TGMSAuGO7OiqoTtVdwDZJHgyRAQCiOVqRwRZSkkxPAW9ltx2rvSTmXRRs92lToAvs1pj9dfSDQNYBbS8fTA3DXafV9qzWHLNJekkBVcudYbdmgQgq44uyW+h6EOCcCZEIg6vz9RoOG61wHYx6C9tr8q5lpBlQUuVV/XaqMSBLChVsk5MT+ooeMsCq74wa4I0ZeJlRdTUo+I2StSKH8mRtFwgAM5ioBK0b29Whwihsdw54+1FVZZ3GGb6mgGvyHh4B0Hqz8lym/FNAa+tzQFsyyXVarYYRTksIFk87BbYrV/+qw6x1lvGGephrlv7n7GTIIxutfps04duAyOGdWMfGcLHn5F3LSDOYpOASPqM2glL5OcsIGWxXhWFMs4okFXANfA14dWoTwOBYhArbVXBtJYYSqF4sZCDRJOdglRTAlkOBS4KS3LbqTLlA0WrngBdwzqU9AOlVhvr4CIQp27X8U1EXPdDa8hzYTkUgeAlBOrDtJQT/kjeHWF2vQGt1b8qPYGb1MJKxYAVdCkBkQqJcx0N+QUT/Mjl8PO4RAu5ek3fd9SPNfu/VnwORx2tqPOtWFZ22ZRfGODRaAQ3DnUqFZ8yisA+noym7VQ12cJ7j2nPxyyiyQjl+c7bKWyOApNG4KjEENKBrTqRmJhl3FC17m+PgMqPAboe4yjkH2ZyUUNZnJIR9gq33F3iGmzqg3VQPOfsPUCQFY7cyIWVYAh0t7++9+nPwya/7lSu7//XGH199mUvMRUSfBOBVIvL3d0netYw0Y8duuTo3eo/wyoMu10pvoWGjrhxq5Y+hAq9ZYRtlK+0EugOaeQrRgy4jR0hkiSHqzNu6p2e7qIm4/VBYoDqZ+kawLQAfqvH0IV/b7t/bXIztFKu1dYurPQTYVt22BdxUJ6CYrYdlUt7sMBsEWAkwCIEhRRu2l0j9cxXniu0YZpxYZ3mOtK8E8LVQzfZHdz3GXT/SrJneBDULk5+Mz1fw1P1ppZdRRRcoIFJms7VLp78Zik579aDLqjKASkxvm7uXhab13A54SxGOtFFctFxzQNuvr9NrDwG29rL3coKvg1re7tqAIi0B+SWce2rKcuvw4JoJz0VWZEJyKDs2wCWiCOBFAL4OwAsA/O8AvkdE3n2R4931I81qSJi90S2ovVZo2+5ZhZQKL4Xltsc1VlEBMwcRIJI2HGwA3YAc85nDeqxHaOVbQR+gIIM+CCGzcxtW7EF3TtcFbD+Xa8Cxxjngvd1tCmj9dh+NsA3Ysl/2dUnqEN1St+yvq1/bgm3i+vuUu1J9PQx53xqDndfz5yDG4MXV61ruXkI5hF1L0pz19lPQHv2bADwBwDMAvIuI3isiae0vJ2yXq3sPEX32ric4dvOT+VkoDBeWQkVO8I4KdVy0YGsVP+XK6+N0y3bob0379dJEGjVEY0QWn9nrym5EEwALPWqGkHqA6FRfDyBlvWN3wBiIjo2BXMT669gFbJvjdGDrp7nxqRTbcMPp0K+Ve5YMA9b1YGv1ysC6OnJ1+yrXTauD5vQt9a9j19z9lRfEoSQF94La9LfJiOgBIvoNIvoQEX3HJYp1j4h8hYj8LyLyJQD+MoAvAfBBIvobux5sGWkGY4xUAKuAHnIlQAVIA0WBVvTVpKQg5b8xSmW3yjDNYxUz8yRk1lFYrzJei9s1uAwOlAnKXGw7k0afZjKDmMvHRFmhMOkgs1zTZ3NccBn+iZbpAuNIgNuV8U411G3B1u/fT9JoAOuB1+qR3VMD3jaczOqZe5kbyHKNQCharXvZC6RhuVOylgkNDNKeEEk514qBGF1vDb7O5+s5eKgfNc/jwkdRGeAHoBLAIwDeR0RvuEhaAhF5frf+LwH8Z0T0iQBes+vx7vqRZr6R9XqdZ699rKOv+G3guZRBWpOiNQAAIABJREFUDywu+TZZPgWt9OSA1mIlOXuPVTqgDMabpQXD6SlpIUByJBBNSgt+WG0PsOtAt//NsdsuYDv1mzkpQb+rvxtpvp6ZYVpKYDumUAXZDIy+zlkKRqtzPdiuq4MroeIwY1FHanNMxy7biSgP26vZ07m+AMCHROQ3AYCIfgzASwHsLQ+MiHwEwPfu+rtlpBlaRuO7ddVjWyu+WAVHW9FZalwri+T5o6gF3WwGmEUAsqGWeXtuG7AYym1YLjL4EQghD/ktzBZ1uTcRTWpSE9zQKGLhdgfdXcG2kRtmPOfcgZNnuqWXtCO7NSlBHMCOImKKbDAG2411EJR7ZnV0mmfmU2U9tIN0B8B9AhE97NYfFJEH8/ITAfyO++4RqPa6sxHRB0TkaZfdx2wj4O77hMdmTcPLDQZAqZTWKEqlxrirZxXdcs4C2pVHZhoedBVQPQjauccsN2Zmug3LhWOqcIDJogly7FpHLHfLjGG3K+huC7abftez26njTQH2tuwW8Iyzslspx2mlrDmwXVcHyc5hParycpBSJjvX1NDmqzYb/rylfVRE7pv5bqrAF62gn0VE64KQCcC/t+3BtmG4ez3hMVr/JjewBVpwtfUCyhhX9NK9zyzTV3jLEmZygkkJVvmLxABd0IYqNTfumjqjgQ1UGG1E/TQtd+73m1jutvfw2EB3l+7pLuy2/b4u+09jsfO/dcDagV1f56pUUCUAW9+lDvrjwx3DW+9MPKyTdG+M+hEAT3brT4LmQ7iIfeYW+2wdrbAN4O71hOuMiB4A8P3QMNIfFJG/1X1/BuBHAHw+gD8A8EoR+b8vc05xb1R2XURg/EosrALmjkCp6L5xTVX4SJR/kyWKCS03urpmjCh2x52TFaLbx//GGK8tM1AYr13zPlju7WIXZbf9es9im++75Sk5YVQutPWtMFAHssZuAVf/JsDWfxZ5CVXfV2bckgo7n/62IyAHCtVS8rEXwH0fgE8nok8B8GEArwLwNRcq056l1I2AeyjtdkvP4jcB+P9E5NOI6FUA/jaAV+67LPt4yzZhtmhZSpg5vO9C+uOQkwXWntPtJ06uGJ+HiqxAeyQwx8Ry983M5uSEfh+zbesQTwDwpjvYM9+yfaszbi7Pddo+npuIrIjomwG8Gco/flhEfu3SB96D7RKlcNW2jWfxpQD++7z8zwD8QyIi8cLVFVjR07quGHen9ethAsnKePaZ9bJ9Qj6Y27eU0THZ8Xm3A+z+eLvKCrerTckJu/52yjYdbyrZTJ/Pdm7fqd/q9uk6yKIO1b4OHduT3Rfgi8gbAbxxLwfbox3TsI4pz+IT5/YRkRWAfwfgE/oDEdFriOhhInr4Ix/5yBUVd7HFFtuvbTfo4bC6sisd0Sfm+NsL2y75cO+BJm54ClqZEAAgIq+9TEGwnWdxK+9jDg95EADuu+++S7/ETWMlqmezdIubWK23XkqYkxbCBO+Y27eUcQ0LvQhD7Wf9vZPNojkuMpOw5ZSdMpu8c/63Gt7VbIMgTVTzft+p3+r26fMFosk6dL0CQms7RikcxIiIAHwPgG+G+bqJVgD+wUUwb5erewOAFwK4BeBPJv4ua9t4Fss+RDRAoyP+7R7O3dg+QKa/sVaxNwLnxHG2LY/fbx32E3TWh33qt3rO4wHnfZdFp0qq88HN7WO2yzObqyvzv5nebx9Qdd0vWEsgtenvgPbXAdwP4Oki8gki8h9AY3rvJ6L/cteD7aLhfoKIvGDXE+xg23gW3wDg1QB+Hsq2f/ay+i25HLNa+R1odftqQvE6PBY5PlYHGbQhOUDOzoXKOmw6HCrHqcnJLVm5P3ffgAJqKkWbKqdvsKPG24CwjBpUvz4HVH77dTfKy5gfhDLHaPvt/cCVkrzdHccz2WZQS97fR5UEoVFYj5+JA7Yftf4DRnZKitWlHO6VIxFE6owefR3UY+bRhaBS54Ba98o62npCYR/uuM2m13lMnBsA8PUAXiAiH7UNIvKbRPR1AP4FgL+3y8F2eSm+g4iuLFtY1mTNs/hBAP9URH6NiF5LRC/Ju/0QgE8gog8B+DYAl0lKUawHIl3XZQ+Itm5MJvcvMgjqximwbStzezxCrewFeO078oMc1oOcH6Zr+zYTWa5lZt2L5gJywjGxW7NdymTX2rxYIGuP4aM8xvdcygty+rf598gvYHLPisZ1RMuY/9z6LnXQHx/uGN76F+thn+s4gc7c3wHtxIOtWR7ae7LrwXZhuC8A8NeuMnnNlGdRRL7bLT8K4BX7Oh/QdcNdAzLQtYZRgBO14lvsomcZ5bhowdZYhWe3vgtagBfSAbzU3zVlckN7aR4sRwzWyQn7YrfHCLZmUyz2wizXejcC+P6I6bjN5Jzu+AQN7Sv5Lgg6YhA1PjblUYDK8nKMtjtLTbdJBcy5nHO7OugBu2yDlDrtQd9LJ4fs0Rwhw711we8mbRfAfeGuB79dzL/JrYteGAgJyAMnWeOpQ2MNdLUR1FE9PdjGjsXG/DeEdrsBe3Tsdp2cEJpyS2G0nvFOyQnAZnZ7O4Ot2bagO+U8C9CE7ZO/h4JgMhCUKi8JCJKBV0F2LCuoXJBHBJLJBtYjUfAdnNaQULebWcjfVnUQ+mks17NrX8esbs3Vmasyuz9HZp9LRH+E+u6zG0IAbux6sLs+eU0DtA5krDvvGa1nCjE3hJVU0AWA6BmGA9uAtuL32m3Pbiuobma3oZR3zG5Llzc7fDy7nWLG/vd3Atia7Qq6/W/qjMeV5WqkAJUXXM9y7cVZRgRuYLl6HBqx3CktVxML6YF5izo45Do7VZett1XqE3z9OaCscIRRCiKy13nbj2ngw7UZZTCyhlMZbmWcBGUbKaE0PCLBAAVdnUInA1V+GfrKrAxXG1cMjr064LXtkWQndlv1QkEM4ph5x3QnpAT/vWc068D2dgJabz4Vpdks6KIO/rCufNk/g26SltlOsVwmqWx2C5aLPPwbATp7M6p0VSaB5FyeIJryEyjgDIzr4VB6TlRe1kOo9boQC/h61soKh7LD5+Bdb0R0A8BfAfBpAH4FOmptddHj3fWAS0EaIIuBR12qgQhDZi6RACndPEKCgq52hVxDtuMb0OZKHUML4CYrRMeIrFs3ZOCNpEAayzb9K+ALjADWA+9UOJMxmV3A9nYF2t564N0adDtpgWBJ4lVaQNA8yQMYq8xVg1DOmzERvO7qkcGkzyymowOp7LtiXbPeFbLUlaDls7L09dDqoLFce5kPJBjyCzj6F3nXa6JwmOfucwwfkb0OwDmAd0Jl1f8YwLde9GAL4I665dXRMWQwNB03EpAImnCGDFS1QfVYVD3MrX5bKnsHtpHsfAaum6WEwTWKSO1LogfbXkqoLHhas70TgbY3D7xlBgyn4Zbwr2wGulPSgoGu2HJ+welEotrQVlgvLcRyEqCVDCvoIrNeBLg6l/ViYLIetg5bVw9z/bK6WesWujZxuOc/N2T5Gu2zbS5HIvohAL9wmYNtkw/3jzE95JqgUQqPv0wBrtuIGCFUEArIDDIwiENer2AbCc2s0X13rmwvgFxZbQ+868DWWOzglq1cxnTJbVeWwoWtTIFtzPGUU2C7ScNtr227VnFotnLRcvXAa6ALoEQDlPUNoFtZbJu+KME1tsBAo1Uqwy3LG0DXQMnqVMr7svtFvSee5TpHbek9odQXq1ta51nrUBAQHc6VdWySApTdAihJcS51sG2yhT3uUmc4dssVSisXO+CtYGiyQqIcYO48x5YzZjzgoHVOeMCtoHv1YBuDY74uRnRKQpgD2ot6qvsBA1dpu5Sx39fLCUB9lga0Zd1+LwByr0dEU2wmPgzoEuqgCvsFhfy5oR5aRIyxW5MTrK63kTn2x8ChJAXRWYuPzHyUggC4N5NQ4AKE866XFD75db+CD335/QihAtdYL6XyqY4PpQ1B6rTTUy/mUWSD0868g0y13f2Bra1vkhB2DQG7iB0CdC9bVv/7nt0CDmhh0mmWHBzbRdDIAgIhIewVdDUGVx1xxmYtl4Ilri8ceKYe2kt+aNht/Sx1PdT6Yz2/T37duvkH9msHHra70SxKIeeS+TToA/3XeUzAzraNpPBzIvJsJy34O3LbSwoAypvc65/KEDMQigLkIMZurWFlViNj7akPv6kMt2pm4wo/dpBpGbhpFAa4Qy6zOchi4I167TqgnQz92tJTva4reJWguwlsdy1/A752Drceqc5uDExLDABDiAAOIALAjFBcop0+l0E3kICEVN9lp9kH5OQ4meEGgCRHSABlMkiLZFhXD4uUEHKvKlSQtXpW6r5rE4cycxoek+WcLd8L4BsB/Db0sTyJiP4xgO8SkfN1v+9tG0nh2fnzjpUWLFIhhsxuM8DFwIhS2e3QtC5tAPYG6vOM+mGZBq6FaaCVEFQzm49GqEy2MvACsITKaCckhCn5YK2D7IKhQPa7Y9Hgdr2Ofn+Bn0pezQNvAGBzx1FOTN5LDDEwSAgIQWdkYCAQ5W5zKMcMgduMY130QshkN7AOstB11Y4L8MLmwBtfmwGuge0JAUPXo1KwRdMGDhmhoHaUUQrfB+BxAD5VRP4YAIjo8QD+bv7bKWJhl/SMNwD8NQDPhj7fdwL4ny9KrY/JDGwNqAqwUSihM6rf2ojNyjRsvjAA8LM5NHkMOqAlx2rr8rSEUBpBJyH4sC8vIRhLmQPaEbudCBfb1iaHxEKuHXTnwHYntuZ2jZSvieqknDYP3hzwAqFIDHo4AkVWsA0BJArGSYwJZ4bLIYOwgixl2crYLolqnQF5HjqxCSC1rHN10ZxmfT2rEoPV+TZE8tARKscXpIAXA/iPfJIsEfkjIvqrAP4VrgpwoXOJ/TGAf5DXvxrA/4o95za4DqvgxMU7G4NgYMGKBENgrCRAqI3BNTnBZk81ltuzW4uTNQdco91OsNoqGdSGoE4OHrFa+65ntRaRsAlod3GQrfPuN9svCLqj41ygsU+B7UWiLUZ+ebHf1O+ngDdxyDG5XNhukIQkobDdICmPqNIomFBkA/2d6rXUSAx++HBqgFiB12SFOFEXfT20+jaQ4CTX7VLXgn+B88FDwlQWufqXNRF9H4D/BJoL4V8D+AYR+cO5Yk1lJBSRRBe4ObsA7meIyOe69bcT0S/vesJjNArcAFPMwFtZgODEulZcmYZVeGAcKgageLiNxYYJoK3RENuzWtNqe/kgOJY+KSGgbUDbDNftR2WZ+YZhTqSL2qbJHC/T6Ld9oYz2G+2QWawbxWWsl4SashoIMwcnMyi7DRmUKTKIJQMxO3AN4Am2GyH6eyiYJgFiroMi7aCMubro/QVDqMBbGW+u9763F/cyP+zWNpfMfc/2FgDfmcO8/jaA7wTw38zs++tE9PUi8iN+Y07P+K92PfEugPuLRPRMEXlPPuEzALxr1xMeo/kA76phVVnBKqTUlE1FVysT+rn26ofPmnZGjmFUh5nXZjstrQPYPgLB2MgU0DaRCBizWy3jerCdAzs/OgtoQ6r6UVjbsNxtgHrbySnXySPbDk+27VMvkei3kRudBoIQNXIDSQVXkxlMv6WgQBuCRR4EBCIEpsJ2Tds1tstCRWbgLDEwAYMIVlIVdJ0d2F9PvX7zFwxhQsKiWu/bsLADMlw5DMMVkX/hVt8Dza09Z/85gJ8gom8E8H4oEX86gHsAvGzXc+8CuM8A8PVE9Nt5/SkAPkhEv4o9p2k8tNmbPEbOgMu1Aoo6zk5IahhPkKyjVTZhwNvrZj7XgWe0HmiJNssH2wKtH8Y75SAbMbkZpwhPxLp74Fs3LHZb23XfXRr/HNjSzHLZt7kfjtFz+wISB66V/Wa9lqiw3hHjjbo9sOm+Gsmg9SmACnMNRVaIWesNbtkYL5Pur8OQ19dHGzU5EDAExonvWRnBsPofGTGmAzvNdtJwn0BED7v1B/PUWrvaNwL4J7PlEfkwgGcQ0ZdAh/USgJ8Rkbdd4Fw7Ae4DFznB7WBP+tH3AwB++S++oAkLi8Q4IYIEDQMSDkD2PKccH1nYhMMOH8dpmb+8lOBjZvcNtD3LLWUJHmg2jxwKoe1UM7dg68H1ImFfF5EgLjIN+xTYNmDcAcrcvYkuEYKIDcx1QJyvx8DVWG8BXjeKrQfeKKzSAkGXuQJvlKCRCB5s7TNqeJhk1qtFsZdAW37/0q+aLeMkME5cPathYYxP/6mf2+le78N2qEcfFZH75r4korcC+KSJr75LRH4q7/Nd0BHXr990MhH5WQA/u23h5mwXwH31TEFee9lCHIuZ4yw0DFf/jOF6Xc1YxeRx8nIsYIhJoJ2KqfUOMe/IaHTmwGtB1oDEAGT3rmEFHpGAEKRheT3o+mvfBKaX1Xs3Xcs6uWQKbD3IbiVdUHLLBrT6xJlpBL6SoxdYKMsHUrYThQK8GnGg6R6JTZpAA7YG4CYzDPlYpn1uqpNeuqoM19W7UB1mhzbB/vLhisiXrvueiF4NjUB4/mWn6drFdgFcP1HkDWhhP7jf4lyvGYiVLhYzdJQZ6YCgzHC9rgbqEpzYpwPZhuE6jXYbRlujJlpWG6iuexbbA+wmJ5nZHLNQUMg5GIIDlQ50Lzu4Ye63lxlFNienrAPaXaM2Kuu331cA5uzksqm92d2nwoQNWCkgisoEURiJApKEwnhXEgrjZWjeWAYyw1VALl7+NXXSXvInVs/J/QVp2sDh7TBToBPRA1An2XNF5ONXfkJnuyQg/x/8OhH9XeikjneMkQO0WgkDToRrtDsHzQDlmJyBRRMZ4D5b0LXohHlnmAfZUp4cmeDZrDFZGxHUO8fWMb21DiGMmYYBr2e722q2Hkz7/TeBdH9va1mq024bUJ4D2xEYT0gQU8cxDddbWx8YIWR2m8GXhAqLNdZrAGzAGolUv3VSQySp32ft1uSEwfThfH4elWNcLyvDVUmhOojHL/FDmuBgUQr/EMAZgLfkF+R7ROSvHOLEl8mlcC+AT91XQY7B1HnGiMl174UhIY/xZs6e4lzJJyplw3BLBa+xsj7vQR/eVbt0Ldj2bLYH2SnNdpdudx/S5K/DA68H3f5YF2EmuzDibYG1t2nQrGA795La5Fybu78ildXaOssYfI31BmKwBAXdwDmOV6NjlOFmgIUmT9KczJS/U/02ZdZrrLlnuGb28m+jEiq7Dc5hdh2ACxwmPaOIfNrVn2Xadhlp9quoTsQI4BMB3DH6LQB81k8/BAB473NfnIf1MiLnkB/iLNwyLLZyrlKbc8tAdmoYroGwVXgvG1Qg5hHIBkux6AB4XVd5nVXgrPrs1PW0XWfK5QgjlrtJVmgYdLffHGD7hu9Bd5OWO8VUe8nFg+3cC2tXZ6OWrUovgL6gq6zAiBFgDg58dZBE4lBAlwMjsOq7KYTMfpXxqmarkTLGei0cjLFd3aystkYnmE/gL7zlzVtd51XY9cD84WwXhvtit7wC8P9eZqqJYzZrjCYpGF8IIggIYFCp4O3v9DMUwK3yQS8bkKvoHmyNzYbADZP14Gtl3KTZrjMFrFSWYzTdsd9HJkEXQONI28XWSQxT5QTqdV2U6bbAOQ+26/TwfnneFEBD8Ew3a7q2HirIGvMNzGAOhfUa+LaartN5s9zA0PheBdqaj2CqfnoCMBUSdl3MFjhcHO512jbZwvycPr8K4IfuVKA1M5ZpjEMERb+18fCRWrBoGzAKyPZs1od1Wcxjz2ZjTCOQndNrKYyBdqrR9MBWgKBhqTwJvEBluvUcY2lhymrjb3XFqTL573tQndJtbdv2Om4LKL1eOfUSW3efp9b766r32e5tKy2EDLwefFOKlfWSDpxIrCxXB1eEDLJSJIVIFXA9cE3V0daHwG1P6hoBF9hflMKx2jYM93Woc/r8JQCfjUvM6XM7WNVOa/cLrDnzSAghUM7b6RpvN9igH3o7JRv0bDaEVM9tGm4HspbZzM5Fzpu8bRSC5GQpIgSKGWybmm6SQcswSxoJTOc9aIBmYpTZHNhOsRq/bUpGmAJZQR16a2Wasl73LtvcvbZt/f0GcOF77nVd4bGuG0JmuBmEmWNhvfZSNrnBSwwxO+IS1/Cz8oLrnkNPBAx8w8Q9ObQd0Gl2bbYN4O51Tp/bwagAJmujYVT9FhpU7gHFGnlht8YWOjbbA6+x2WYgg5MSQuBmkktr6E03eBS4P3aGAcpQhW0klI2MyoM5gAK8UiSCCrr9sXa1OXY9BbRTQ4o96+2124sMhvA213uYe7Htcu89/5d8I4tjCxWE7b57XZc5QDKwihBCitnJVnVdY8U2mCKGGm7WsNyZuupjb30UzHXa4SJir8e2Ady9zulzu5jvYgkJItRpVuvjOF6zOsvawQm9bBAjFzbrt5Uk6JEbkA2TIJCRcM3QSwKaTCaSGytBQFnzQ0zaQG3mAU9jO9AdfT3BhL15QN1VWpjqCq+TEuZkhSkw7R2Oc2Dr778dy0sLjROtB9/RzahabkRluQUgc3SCCIGT03XzsrHeyBEp1fjcnvUW5xzIgftENI17+V9nKFhvxzbjw75tG8C1OX0ArUf3+Dl+7oQZH3p72tvehF943ouaUV3qDdYIBW9T84TNOcF6bXYkJRjQ2m8L8FZwnWJXa3XEWJktxIA2NIxW5Wk33UsHuhhP8L21NYyuA9sprbM3L1V4kPXM9vIsdz3YTr7sdn0W+RZ6putB2IBWk9TU9RBaaUECl+9DipNOtspyN0sLPvwQ0Lp/XSY4TFjYddo2Mz5cvLXdxqYhTvkTKPICyZh5+RFgvms2JRuEEudYJQP9zN260uBro24a+4Qn2Td44TGICYWyLpYAwmYdyOA6BbqSfDRCe8zemkY+wVI82G4DtP33HngZVV6w5bK/6j1lnMo6HdfYbdm2AWzXPRN055p7JnZN1qPwTFc13VBehCkGCAcwcyM5hBSK1sscipPN67r+eVi+3r7X4Ostuvt4XbZICnexlZFdqA0Y1LIZX2l3AVpjr8Z6Q/7OGnSI3DRk38CxZtglRcBPTKggK245NMBLrNs4hRHoVgCpLJdIgBlN1ptv7P22/jebupEBMgJeD7QGKD34zlnrLGsjQXqw7Z9LAV/34mueDTB6PtRTFg4AUgHWGskQipwAJgSp6ylFBGFwCtD8CwEpae5cA17KOq8NF/YOtP5FZ/ehvFSOIgKWFklhMbVA0oAt0EcktGFd64A2xFQkgxDSCGQbgM2s184HtF7y3oQDELk2LtZx+kAFXwNeTgESAPr/2/vWYNuyqrxvzLn2OY1RjNiCbaBtU6IJMVQn3LTy6AaETpBS2geSGI1NSqWU0jwsYijjj0TLqo7GUitJmXShJWiigoJQoBBtaR4N3dIYfKMNCRGlA7Yhptp7z9l7zTnyY4wx51hrr/04956z7znnzq9q1177ueZea69vfeub45EzQsSAdDlFWa/6syEwkile4iPNJo8P+PFl7iaVawLRiHdAulitYut3D0PByrL3Y0PdtqvIdnL/AJe3j9zv5kxADmDOCNEm0WT/UFTyjUK2ORJyilsRbwk1c1XLTjMkxO1qj+Jk0Qh3BVZ1NBjYCC7ExvuxMeQS1lXIdwXRmtINWnuUiEFdKuuikIEBIawhlyx9szhrZVZTtIHluVwtBAQpouNJlwIVe6GoYhjBbZ/CO1ZV4wmyKaJdV7ymELQRL48m0lCvQKz1zSZUP3bagyeLx93iZOj3E4At9lWqESPaKK9ELIzINyt5InHxczPxEvGKApZ0cyKW9+Vq4/gi6X7b+5PVaUg6uOY93GsVt9z7FgDAe257EYDVRCukWhXsOBJByFUmzMZEG5ydALMZHMFSYMCC8O3gXXUQZwJFI10AHOTgJRZ1o8RritcrxEK6WlWdE4GcyjUvd3K1I894yr8tr02Q7ZCMh58hmg4Jy6CyGSYnqdTHNQLMmQb1bOv3D7dtsXAKoQ7J1k6OU0R71P1FnnCBsr88AecUirI18h0Tb6IIypVkbZskYhAFvTKpnq61Bhrr7y9+x5snx7lrnHO+bYR7FKwiWx/iVcjVqdpoB6uRrt6bmiVikClcPWgp2AGfh0pp6gDOVKIRSDPFkBkgAoVl4s2A1PWNeUC6Qhah+LRDZb8+WmFJxY5Urg9V8u/3JDtWWMEdfkbUpagO1U3h1a332a3lzTYYK91BiN6IbFedGP0+AzaQrj9RWeRIIdwAzlmKlKcI6gjcx2orKPESRSHXHGQ5cVG7ZTVZiqcns4sAgIFAp0PRekiUwuka03GjEe4GTF162WN/YE6Rbdf1RdVSzEOyHROttTPxJBuGl6rA8mUqK9mCayIDZeu/NkG8qOREmnFmXSyYGUiyjsBcylCOLZVViQzAcuztKrItbb09UY/Usb+8NPLNqGFhY9KNE0Q7FTI2nDSrXmwYWwmhTo6NyZZiXn1ytP0GLNkMg/1mY7CohUxL5Gt2kHnvRrJEUe2DAFK12684nHMGSjjjFhOeVxMtSmEHIKLHQfoK3QTgIwBewsyfHL3nZgA/BuCxABKA72fmlb2IjgvPfNcbAQD33XpHOXhKGE3gSbKNXRpYCGGWKtG6+0K0XS4HK3Vp4Nn6CRnZEG6ZqVCUeX/FwzVfsA+VeBFhAbaFdFl820IqlMEcy+Uyp9WqdlxLwQhzPEE2DMLXqDRbdp9Zhk9s0M0BHqhdI10/kRaJl9J8mQN8p4alNY3C7YZRC9NkSzFVou2GqlcGuHrf+fDC4q/Debp9FKLtMrgPA+JFCHoyCMghltSkDv0k6fpJTz/hyUzl/30awGiTZrvCKwHcw8x3EdEr9fG4bfFFAN/IzA8R0ecAeD8RvW1NP/ljh1dOXt1aPOf4EnQV2YZZX8nWCNcO1s417lMfEVg9611Sc8nLPaqhYR1AmZGloio86QrZstaIqCqXeHXH3VXqdsk+GM2Kj1WthSyt+t7yubJkxbqHatcrWbMUbJm0jfk4UmFdokTxbk3d+nCwCbINRrr62oBkt9wYg+YOAAAgAElEQVR/EsqX61VKH6sF1EcwZSHcFMWHR1ftE4OSbuBaAS4Elv1Z/qvLls5pQytesxvcAeA5uvxqAPdiRLjM/Idu+WNE9AlITd6dEG6pHDbONCsE7GJrg7cRKsmGKM8ZyYaYqqI1onVRCghcjqrJySEmIVH7l2bSwjSOfEkMO1G0lXRZJ1o45GIteJVrEQvbwAjUlm1sXt2aqp1SwWtjL5kcwYr3mFjUq4WJeWvBOuiayj2SjxuGJ1FvJch24SWy9VcktlwshaJ09funxlHOQqi+rpGt9957Rg4M9OLFiroV8gVQMgZj0DThEhoYSvJOWSVqyNipArcohV3hCcz8MAAw88NE9Ph1byaiWwDsAfjwLgYHALfd9wYAErUwlXHmY3DNxwvRT5ZVsg2zXlSt2QpdAmwiRknWlNZAxozaMJClVEUI2QaAcpLlPsr3umUjXVgImKnaFcRqdSTGKBNieahix+o2lwObBqp2OQ53+N227vKcnuQCK9kS66SPbIsYqrUwyKry6jcTiKYJ2CZBV8GrXai/u/ZkOd6Hft/55YxqeyjhMhMo9OUqhQIXxRsCl0lP+xojWs5Ztk2Q/6CpeE/ytm/yKbMSDIwWpXBsWNe2+IjfcwOAnwJwJ/M0VRDRywC8DABuvPHGI450PXx++vS6l+M2hXQnyLZTEvYHadRA2ABQ6dezYmUZlSwzg3SCjANASKJ2EYAuyUEbGJSzyETmocplqrbCaqtz5TaxA7k8dr7tmGwzlifNxtvUfEaJmNBkC0hH28BV7WZWUQnhKFKitZOi2QpHwYBgx+rWTaoNyNZOmB0PT5ambt3yEixt2vZhBjgFUE4S3UYMpPrhgL6SrpVvDJo0kXhaSSs2/X+vNnapcInoFQB+EMBnMfMju1jnzgh3XdtiIvo4Ed2g6vYGAJ9Y8b7HAngLgO9h5vvXrOtuAHcDwIULF451Fya9/I4sBbgz5xIsVUOLHNF2ySnZ6v8NyLbjZaK127rqbEEv0DMq+WYG9QB3qAcsghAAZ6CzsLGhyh2GBPBaO0EOWpevX2yDZR93XKPVyPaolkJRtmBkogHpJjU0A1DWZSrXgv1Jx7zRJ3GTXAPv1qlb89xLSFg5aep/oOPl/ShfOL1OKR8mUXm6HylkoCcgJLn3yATqMsBB/kssnSTETnKhffZ2t78Sh9NnJTjsKkqBiJ4E4HYAf7SbNQpOi6XwJgB3ArhL75eud4hoD8AbALyGmV+32+FV+JJ4wVXxH2DqoI1JZ7NXkG03QbQBwijle4UoC2yRoOoV4B7ynZnBEOXDmaovDI2pdSrXj3Vq2WNq0mwqk8wmycZkm/KQaDclPtj7AnE9D7BsI0+6zBK9QJrG6mNx/SZcp/Cmfv+Uui2TZ8W/zctXJx2G+xHYvC/dfkQGODCo1/0Id9nBJIkTMYHyaAJwHH5m2zq7jhKnlHAZQL87hfvDAL4LE1xzkjgtW/4uALcT0UOQs85dAEBEF4joVfqelwC4DcBLiegDert51wN97nt/AYsUpZlfDlrBiZYO5HFxE1NGgwkWT7Yd5PQXCZgFYC8CnbvtzYb3gxvJ5/RAt4Pd1FadzHFhS6bYsB3RAqOIhJF/W6yDchs+l3jYd6vnoI/FPkhG1O5mz0uzRP2MKbVsrznbgoc+8tIE3kQltSlMbgMXmzvYlrYf1dfw2x+RZN+M99eqfbkXZd9HGvwn5D8y9ordeGAn9krg/j9p/9PEAYsU8dz3/sLK3361wVveAFxPRA+628u2XQcRvQjAnzDzbx7v6DfjVChcZv4zAM+beP5BAN+syz8N4Kd3PLRJ9DkiZmns11k5PQ6TB7E/CChkDGay7SA1VRtJVFAIw3uDqaKgJkZmVUnuUpkAMT0hB2kPcICsj2W9rF6uKGAWzzbkgU9oqCQbdJWjMDCMDmxUsivkyMsqF8AwGWLFtpa6WmonqIqF83EBuSK3KUFiDZHLoSRIDG0FbznIZ1aSr0v1BezKoKpdBLcfO+fZRnd1cjn7sstFylOAToay7Ec7YfLwfwWgJE/Y/mH337SeaH0+vdVWGUfycB9h5gurXtwwZ/TdAP7uUcd3HDgVhHvW0KeASBExZPR9lII1JUuISngW7DKQMpa6AxgJemUbVA0NDtIRCdpj7bEms9ks930CIJNM1Km9EFguQUMGK6FSGJItKdlS4DLkcfTB8D6Ux9nfYM0OKxEnR76l9YtNqI2+fym1l2q2G5TcCVJIJ5OFi8m3SLIdiZ+LALj6xWXyDMskO/6d9miqvoJ/7K2EMuFpytbIdrwv/f4bL+c8sS8BSw3hLNsAOpGHJH4tT6RcM5PWz5X90fdyRdaniH7ipHpqwMfn4a6aMyKivwng8wD8pmbfPRHAbxDRLcz8v49n7avRCPcy8IL3SYLbO5/5VYiUtTJTKsWgY5eGFoNlZPlJGHKKiCbIttNdE3y4goLzyAM0ZYThgTqhjuDIlTWNd1VUwiBN1yknsxPG6tasBDbiHZFtsQBQydpGX34a1+iEEiJm5EhKPiBEZJ1sc+q+dFnmMiavcm3sIVRfs5y/mFb6a0vVv4x0112lTJ04bWVT+xOQ1/t+uC8zKqFn3X8h1HA9DsPeaCz/QbMRco5CuqfcSjCcdOIDM/82gBJ2SkQfAXDh3EUpnEfMUycRC31GCBEhZHTcS3WnUg3KKd0xfNhXOSiVbD3RjlUuQlW45SC1qXpVfe5AlQBWKrZCOVidol0Ff+LwywNl6ybKxo/Xka0tAxMHmpGtbj1TuSLqWWboS9Sm2gkaFkagYi14lTuIXpj4TWvhEkmKnWDPj0+cfl+GUE+eNHHFYj8+Bo0i6SrpBgZCkg2whTDlTIOC5X2K6PuIeeowT6f/UD+ipXAmcfr3winG8++XYIn7n/0VCH2HLiakPkqFsKTl9bgWJBkXDgGwfCDZweoPzqmD1D6c/UGq0QcIcqDy8K3Dj06T7SryqSFfw0nCsbpNOrElk1qVbAcTXPZ9bn2rUk3t5CDK1s4jpI09/U+wk1DQn5eWVa5TgZlzUc5bYxQBsBRjW06cNCTbrfYlyjlTrl7czpn6iO0Lm4hMAZzUq+1F4fZ9hz5H3PruXzza77yK2HUtBWa+aZfra4R7DJj3HQJl9Cmiy1EKQ+s99xHc9cXXLcXBx/AH6eD5FaoIqASTsjugeXCwmq0A9WrZ2QljlApW5kfnUE4Wfqbf7IQpNVtshWInVLLtJ1SuHV8rCdd+h5KsJ14Ag20ZQGInQMZCnAcq1ywQsxXKxFLO+ltzyfYqdYVHWPJyp3bN2ELw9+v25eDzdvKc2CY2R2D7qR/+53JWdZsi5v3ZOsRbtbCGjTjsZ4ghIy5yqaOQu6C1SwNCiuBc0zUlA2yNW0Vh9UELDJSc1N4LelWdHXFTtRUgl+R8xKwrH8M5ULRF7VZ1a+m7NnEmIVzLZJsc0Zra9esaY5DGikq8EUAP+QMPVa6EUZGOKVBVuTZu8r9lRXHytfAxujY4ouG2B5ZtoU37EtAsjmnIiSAAnmy164PdUh+xWHRYLGZYpA6H/eyIP+7qwZ9IzytO8ZTl2cHtD7wWh/oHF9+sQ993yH1X2qBwirWif4lmwHb/sFWRCutI+TLhFS0wvOT3Sjfl4AiWiqpMuWaTsSPbNLpNPbfptup7BrG/oGJnsC7XRIyhrVB+8+j5cZfdy8K2+2j8eDyhBgz/J/ofyik6so3I+p/re/kPLlKHw8UMtz/w2iv/LTtEZt7qdlbRFO4x4aDvtHGk1sZdaFnGPiKmiLzoEGY9OGVwSqCeaqVETcktsZjBzLwhWGM2KScsZSqRm60Po4m6DTxsqmmcyDBO6Mh5aB2Yus1sWUwY+LcDksxBM4nGtoKub8XYaoQClu0EJgw75FYfN3CSNOxQIxaSFRliLr8ljJS7iV05IW5xApvy4A1j8tTXBvtxDM4aImb/B9RbT2A9eSMT8qIDp4jUy+RYv+iw6DsczvdwsJjh4IzZCcBw2uE84uztkVOKFz74swCAdz3rKxHnuRYlDxm9Fq7Jiw4UEzhGsGYPIWdJSLA4TDvkx5ea28ImXMaka+Cqsk0tDV52BOvV7jijzqvbxKFYCUbIY7JNSrJ18gzroxRsPHqfwTpZVn3bAAgpjkiXWFKvKWYkDgiZQXFZ5eYcwCEXH3f8+2UbaLsbt718pbilbe/vj4pBqq8SLzM4Q4i21zmBPiAvOuQU0M9nSPMZ+vkMi8UM8/keDvsZnvOe11/eGK4imBnpDKvXbdAI95hhfm7X97X7w6JHirPa6aEbqdwAIDqFm7PWHBxOpEwqojHsM1M1F1agkJD3Bk29ZnepzXWyzBOvebjm22ZgiWwL4eZxlAIGxWtKZwd9SiKiLOWhEi+CU0OOdImso4G0GgphYrwY/qYyieYnzrRz7ubtPd72mw3hwX60/e33m1e4PfRW1W1edMiLDmk+Q16IhbVYzETd9rMz5dt6tLCwhiPDQsUeePaXl/q4UWvhRqscpuUZByo3aVB7n+osdcZwcmw82bIKE68VQvXPOYVbsuQw9jRVueZqHaQclqwDsRkwJFlHtklVbRoRbWZMRiokRm2vA42f9Y1z8ijKTkk32D1Y1C0zyFQuKQkTI1AumVgh5CX7xCIVyv2Er8s8bONTtn2eiDax11btw5yrndAnuU8M7ofq1sg2LzqkRYd+McNiPsN8McPB4T6e8c43La/3DOG8T5o1wj0hXFrsSeSCWguDduhdBsckNQ26hFLk1eqoWraYxdd2K3aTV0er1Njo6dqG27LgRsRSUnCXY27z4LXq5w6shBVk2+skWhoR7aD1zmjoRn5EQraeeAFC1vv6W4U85R4gzmVcUk7TR1mIguXRby2TgzmXbWLMXnvFib9bKnitErachz9qqUKYLve9vLfvlWxZymapskVP4EVXbmYh9Id7mB/uYT6X26XF3vR/4AyBm6XQcDk4WMwK2ZaW6Ye131mtrZqAsADmmsHUj2wDI91ViskvDyZbXGJ6mXipYWmSiBFqmBGG1oLl4i+TEQ3UbYkM8LYBkzYErGRrVb96riFhFuRu9sIYlmBFWg83OuKFql040rUGiWYrSJ1YnUBzKjcqEY9tBW+lFB/Xrgo4KBEHafiY8zB6gG2Cy1lDgFhD43013meebO2WGDwHOAnZ5kWHdDiTm5Jt8WyVbA8WZ9NKMFwLYWGNcE8IVm/hvlvvKA0mB40nSwHrrBPdCyBkvUSdIN2M1Zekpm4HXqDcuBBtVWrjCTPv31pqaI2vDUu1VAvplipUy2Fb/Yhsk6lbvTeSNc9ula7JgKT0knAaEZR4vagU0rXJLH8fMiEg6EUEIzIv1TQ2WyHngMC5hu6pjzueOKuDs/dpDVt7ydsKY5U7te/GZLtQsl0E8OEM2W7zGfqDPfQHomwPD/ZxON/Hpfn+mcomW4emcBuuCGYtEOWBnysql4uVEIiBIOm5k6RrsZwp13AjT7R9P1RXRXENkwq4KLmhlSBEO7QTSthXtpjWYZKDKNawZCWsIlubMDM/1yZJBj3NTK06CiZyab2qapnVhiksx5VolUBTsRZQVG3KhEhUSHfwG7meVILZKSMf1262LmYqbXFg3XcR1BZQO8jOkUUNj06Qq8h23iHPO6TDPVW3e+jnomoX85mEf50TKwGAhg02wm24Ajz//tfh177kxdJldlyUXO+tdKMclxOkm3OtPgUsq9xiJdgBzE7dyo37WP1btRPAJM/7SbKsWUs2SZZimSwbFqMJJfbWVsNcFa6pXk+21XKoRGt+LjBUuQk1BTpAlI8RL9zEmSddC18OkPhd0m61UsBcVG5kDRXT0LFAtVVS1hNPCErCKUifsD4CMRVboXjgfdSGj7KtyTySThUuQrWDgouT9hEJxbMdke1BV5Xt4QzpYB+Lwz0s9P7gcB8Hh/u4eLhfJmrPA/icR+I2wt0B/mKxV5tLug6/AAYxpKX8nyfdnKRrQGYh3dFnauaAO3i9uk2h2Ak20z0O/RqEg9nEkvMybbIs56F3u6RuN5Btr6o25fHEmfyE8aFWtauLVGC2DF55hyNdictVfjMPd6RyzcuNFJBDRuaAlIKEj/lY4zW2gnXRpcDVrjGVa5XR+6RHl8VFT+wzi0aQS4EB2fKiK0TbH+wJ2V7ax/xwD5cOrsOlg+twcb6Pvzgn6tbQPNyGK8ZXvP+/AgDufcZXD9pxU+AlxSueI4DcAzmD9pRAI1cP1xc18YQrcVfL6lYnwUpEQkkPDQM7Iam6rRNJsUYjaCjYWN1Wkg2j1FsskW3iCTsB1QmZQqA6LWbEmzNcq/NKukEfi0tDxVoITIgjlWuJGlFDwnKOJQkip4AUxVYwlZtTAMWhpWDRCkhCvpyta5Ej1FLGcbTP/P7qMbQRime7p57tPnr1bA8OrsPh4T4uzvfPRH3bo0AmRJvCbTgmWPbPg8994WTfrM5NyEhK60IOyj1VUx2j9AQfw66dNcIIvVO3FsdpPqXm4SNLHr5dPpviTckUbY279WFgXt2O6xsYsfZO6XqyLXYCgN7mjbCacP1rXajE65eMdMXaZqmFy1StBSZEDgOVmzIjktgKiWTyLOtvtyQII9mcIkLMst0CgztRvGIpMGBV2CxQTcchJ4KJZJXxvpqTRCMczpDnQ2XbH+xjcbCPg4tCtJcOrsOFt//S9MY662Cc6ToJ26AR7lXAxfn+krL1MJ8wAMCsV4WblVSUcD3p5nrPWhaXU5D4TbMRmMB9kJtaBxaRYMV0rENAXQ6uVY41fRzG2E75thYmloElsjXCNXEHbCbcchrKVfEKHP2SrMOiGaxiZXRjJB17JJaxBo3Phfu9RfGOtk0KoE4Vbh/AVDsio48gJO2um6XLBoBSc3i8r/LopKhxtlknxvLhbEC2hwf7hWwvzve3+YudWTQPt+HYcdt9b8A7n/lVNavJhWchB0TnH4Y+gmY90ElsJlwr7gFcrC0zVbJV20CW3a2vyjbniKTVzIqqdZNlReXCFKuQkalbK7NYOuwukWwl2z5Xgk08DBUGKvHaFbi9bi3TCzmX3y+kK1cEXBowAISo8btV5crYI2dEawNEoUyepSRdO4LeJ1W2SFwmz5ikCScHBvc8KHdZSDdr+/k1+2lwQrTsMYuz1Ymx/nCvKNuLlx6Di/N93HbfG472ZztDYAD9OXdxG+FeJTx6KErFWwv7PsbTTdxQJlAKoD4L8XY8aUkIYUNI18i2tzJ+sZaItMkwqzrlYm/Hk2Xm3dpMv/dux+pWFK2fNJsm22GUglkMw9+TChGLajWxyJA43F4jEgSEwDpxptYCqYq1CbSo4zQvN5FYCqbYY6iZZ7YdAmeNWnC2gvPBc4oIoysUQpKC76zhYlP7qB/tH01qyHNRuDZBNla29p85v+CmcBtOBlZdzJTuOIe/SwFdlvAtShFh1ksdBiVe1qQJj9Im20jBVK23ElTdlqLVWZTseLLMyHWQXVbsgjBQt1O1a7Ne1pcaLI5sUx4S7TgW10MKp/soBQChki7p5buFJyftZ+bbuVWy1TRfzVaT38QuRCwgjibPcpTtgcRiK1BVuQDASw0hqcZXr9g/gxOhpeuajTCfYXEg0QgHLhrhPCtbg+zn3RAuEX0HgG+HhAS9hZm/axfrbYR7lWEH0oPPfWEhN0AuOe2SP+4JYYZZD0pRiDcwWLPVDMWaKJ6t2gRaN9VK+lnRalO3Wbu8JuuJleokmYWCWbeEQbUvR64WX8tuuXfq1pNtPyLasZ1gCGrRSoZZJV7NgEAk3wOLQMzVw2V53QrmFNKFRk6QJUBQCRGzCcMQAkIKINJtQ1xUrnTpDSCNhZah9PVER5qunTCIl7aIhnrSC4OqX0mzyOZzUbiXNBrh3E6QTYGATCdvKRDRcwHcAeCpzHxIRI/f9JnjQiPcU4JHD68r1bKYCbMUS+4+co0sEKUbSudYCrn6iMUTVkXLep/k8tU6T1gXirG69XZCKUrj+pX5urfj2glDtevrKbibKlsjW7MWgGopeA83MUDajMioKxAPSBeo7zeiDapyk6pc83KFbANi5KJufYhYmUQrk2bLKpdIIhMSgIhexgxIBbhMoC7X/eD2yyDWOQ2rfvlCNIu5VP26FibIprAjhfttAO5i5kMAYOZP7GKlQCPcU4PnvOf1ePvTv6Yc7PsuXCungNhHxBQQkxTWNpULl6kGaIQDoJeu6tmqivX9r3x7Fh8KllIsYWA1IkEzyhzJlomyXGsiJOfJekL1zQsS86SlAIwyzdiiE0S11pMRNB1aW6FDLQUlZ2IWsWnrJBljjFpDIaaieGuHiqGfWyfNqn+bScg2E4NCkBZyJfYXkrkXU9n+S/vEkS2ngGTKdtGVQjRSG0GyyM5jnO0mMBhpqmj+NK4nogfd47uZ+e4tP/sFAG4lou8HcADgFcz8viMM9bLRCPcU4dH5PqzId/FT9fK+kuYCFFMtZh642gqU5cAHaoKD66lWCNYaDnqlqxNFvm6CEW0tWlPTd0W91tq2dRKtqluLt/UTZZmXyVbeO61sAollYU0grQ9mSfu1OgsswlKiEiRAK4KFePX5btAuHSVErBJvra+Qc0AmRs4RlBiJomYIBuRgZXM6cCaErPsiU9kXvGFf9FY8fDErJRYP5/s4mO/h4uE+Hr0G1S1wJEvhEWa+sOpFIvpVAJ898dK/gvDeZwD4EgB/B8Brieiv8g4q5zTCPUWwjLRf/ZKvrUW/9eCfpYiu70XtdtovLWZQl0r1MaufZYXG2Xxcdl7tIpZJMftur2593QSzFjLqZFnp1luIyxo5DtWtv9nkWWIuynNMtisPM+Yl0pXmQRKVYJNlpOu2e1O2FCRJzyv0qcmzyFx+uw8RS+rlUhaFTDkCizq8ACDnIJltJbZ6uC9sP9g+SH0Uz7bvpHi4lVic7+PSYu9c1UY4CiTT7Hg8XGZ+/qrXiOjbALxeCfbXSS5Hrgfwp8ey8jVohHsKYQfcu571lTppJWp0NlugywFdlxC6Xmrs9kq8ozAxm3yzDDJRVqpmjXxN5bpwMF83YTBRZlYCnLXA5sk61epUrd0ADfgp718m21XiIhMtky7Ej7VC5OP1eTvBj81bIeV3aPWw8pu1vgJzLio3JVO3YikAABamXLXGsW5vvx/8Pii+ed9Jw0frQaY1bS8t9s5NicUrwXER7gb8IoAvBXAvEX0BgD0Aj+xixY1wTzH+4nAfKQc8xpUQzDkgz3p0KQjhBin3aHUZLFTMq9ysKb5Gtt5KKNEG6t0WhWvtz90leFnW8ZmdYPCUOfZux89Pke34ULNiNWPS1TpcpYyCWQxZIxPGYzFbgTWKYOo3pUyIoar3lDUcT8nUWwt+fIlJkiJ0e/vtX2J23Ymu7zvprqs9yA4O989F8fDjgWjcHeAnAPwEEf0OgDmAO3dhJwCNcE81rIg5ALznthdhv1sgqb2QulishahdJEi7SBTk6gWzZpBxDkh9DQMbJzv4CTNWK6HcRhNmAMrzRc1pFbAyEbaCeOW9y2RrXm6gGp1gpDvOkKUx+RZla/4tFaULPZSHv6HeSgqzK2hTTnJmFVCABiYgsxA0BQYndtveui5TiXJIuVoJi8UMfYqlu+5Z70F2nGDsJiyMmecAvuHEVzSBRrhnBBfn+zUZYW8uFkPuEUJC7hKol4aVqzLQko96UIuiRCFYZAKqojU7wbxPn5hh93W5XrZvQg3/qm9ed4j56K/MjEhUrIR1qERbt4Eko2nxcuLyO+03Zh6SrqlcK6eZkpJpDwSV0oEt62/UoNNtcznBRSwWHfoUcahke1a7654cGMkb5OcQjXDPCMzXvfcZX43MAbPYI+WALiaEPqPreiStXuXr7XqCtIgH822TK1ZjRGPKzy63fSfdPLoHlkl22zbXYg+omsVq0vU0Fmg9yU7BFLbx8/Rv8CcOZykEJeUckSY8cmZy6pcHz1sKdd93yDmgT6JuF6nDwWJWKsc1VBznpNlpRSPcM4aL8z2J052FQrhd10vuf8jaP216Ek2SG8iRbaiNIp01MCbaYTbYetKzhIVVqI0flWiJBhNmU6RKIzuBNqjbTbDf4D3oceRF5moz5ExiJzgEJeNN27lPUXzbFLFIHQ4XMxz07bBbhUa4DacKVoPhV774JdjvAvZYyDZ3CSEk7Z826iqBodKtPm2o5DtBthnVNliFQLXLuyfbkilG6sOiEndQvzXoZFgmKz6zDCPbQEKzVkUs6LKtH+41Ww7Q9W/gZ1ZVD2BpOySdDBOC5fKYtajDeBsPrhZylIiEHDHvOxz2M9z+wGvXD+aahqQ+nGecCsIloscB+DkANwH4CICXMPMnV7z3sQB+H8AbmPnbdzXG0wZ/4N536x3qO0Z0XQKRtGYHMOgw4cmEucbfFiXnVe0KFRmISwjUGDKBJdEEjJqUYJfqRIwOUmNBYmqHpDu5PgzJlkDolGQJNCDectPRT/nZYbRNPNip77KtyKI4agaZKdvxtgVQUoL7XnzxeepauNeWkBC+861wp3oHXA28EsA9zPxkAPfo41X4PgDv2MmozgguLWQSZp46HGp852IxQ0pyKWs3KzCeUiwFavwEGYCi9Dah9GhD7dMWChHKLdqyEmEkApG0aCMtuyhESuj0ZiVky+MVZBuJBuuLbr1e1QYd41jxr4NXu/4klUrccqx2gd5S6sp2P5zrvmgTY0eGNNnZfDurOBUKF1K55zm6/GoA9wL4l+M3EdHTADwBwFsBrEzru9ZgE2pvf/rXoItJCmznrPaCKbOJSTSnbn3UgYcQVS13aMousJBj1FhYq18QyYq2UKlda7DC4EziLWTnxlr0gj1jhOkfBwAxjIlXbkHvO2J5jPqaPTZ1KyeKajmsiuzgkco1dZISgWg5asPSoZOS8bVWC+FKwS1KYWd4AjM/DADM/PBUuTQS5vghAP8IwPPWfRkRvQzAywDgxhtvPP7RnlkhfUkAAAyeSURBVFIc9h0SB3QhYRaTpKzqBNrU5e/YRvAgiDJMIKdghWgjsdQZCJC+loVo9ZMBQGa5RHchCEyoRBtcxhlXVRlHw/E+rKniSrxKssGRbRAFHYkRnbqNgd1zGP4me9+KSlWZqVQtg9ohJe7Y1K/GOy+SeLZ9Oi0Xj2cJjMzNwz0WbCgmsQ1eDuCXmPmjtMLvM2jVoLsB4MKFCzvJIDkN8IkSb3/610gAf6h+45hQVvm0BiMqTkJIhWghhJYBsMlfq9WSARSiZangFfQtDCVaSwOm0mF8VaV/T7KVeIeq1pPtzMg2MDpidIEHBEzEAzIWQp5WuUvbSt+S7ISFYexuyqGp2ivEWbYLtsHOCHdDMYmPE9ENqm5vADBVn/LpkJJqLwfwqQD2iOhRZl7n916zGB/49916x6TS9TYCESMwS/orWf0DlgpaICBnZKlLCATtXBuyPhbSpaDlebVwDDHXNuUsf7heS8Zaiq+NBahhWuPIg+rj1u69lXRZSRfoVMl6su1CVkLOmOl9pCyqN+TiQw/U7oiAM9Oo07lsu2e+641XspsaBmhRCrvCmwDcCeAuvV/6FzPz19syEb0UwIVGttujZE9NxI4CKL4mEyMiS7eCrA0rc5CUyxCQMzuSzSC7V8Is3WWIEUubdBayZYlQiNHXUxgmJtjIvPYOTtkW0lVv2YjWVG10anam5Dsbke3MyJayfs4RsLNgpkgXqL6thYg1HA/k/9AU7i5wF6Qm5TcB+CMAXwsARHQBwLcy8zdfzcGdB/ieWKZ2AW1dY2FaQElpXSJdQOvMMhY5lOaMC8h9Yqm61eegBEtIJCUZOyZtvyPkWy2FWnvBsJJw9dK/Eq5aHqTkPkG0nbMPZkHJlbIsh9Vka7HMU+FjTdWeIJiRuU2anTiY+c8wMRHGzA8CWCJbZv5JAD954gO7RmDkaxZsDFnUXAYoEijLJBpxKH6ndG4g7Y4bSp3Zjrg0mKzde63urJBvKdPIJEE+G1z26uGy83ItGsGR7QTRxilVW9TtyFYIeTCZtm0YWcPxoKX2NpxLeIX2wLO/vCwzSZ3ZMhFEUpQ7UkYK0oom5YxE0g8sESEyI1EuXX17JlG0OYADCvmKujWSXVa3q6jNlO5A5ULIt6pcIVkCinKtpLueaL2FMFa2fpLxi9/x5uPa/A1rwM1SaDjPmJpIk4aMVJRuytJcMVEAB2swSUrC0mRR1Gyo9/ZZtt5nQa0E1/XX7rdQuFXdurAuQNWpKNywgmQD1lsH0ZStmyybCqNrOGm0SbOGc45b7n3L0nO/8bwXAHAVsUItcDMoUl4KlVei5WIn1CLmYie4rDZUwrX1AMvEVpSmszwGhOvCu4bEaq8tWwZGrjHkQrD2Xbaev33PW499Ozdsh6ZwG645rCKcD9z+90pKsKS9VvK1gjiFdHNtr27FYaxx47BQjHz3qpTiYBloVAnYiJWKzSDkGtTHtcemXo1kA1WPNsaMm3/lbce85RquDDvr+HDV0Ai3YWvs7x0OqmH5IjjlOVdPthCzkuy44DdQky9WXbaP/dRAXMi31HJwk13eHiiWgXbD8DZCmxA7fWAwcj75KAUiuhnAfwJwHaSHx8uZ+ddPfMVohNtwBPz1t9y79NxDdzyrdDaYIuPs+qKVQjDjtjwbPFJvLXiv1VsBvg7wuERljAlPfuO7j317NBw/dqRwfwDAv2HmXyaiF+rj5+xixY1wG64Ij/mUS5VMfT+0QQPKIfGOiXZTzd3xRBatIFpposluObfwrrME3pmHywAeq8ufDuBju1gp0Ai34QrxxJ95/9rXH77zqUvNLIHlwjm8ImvLuuACGJBnJVapnnPDq3/rSn9Kw1XHkTzc64noQff4bq2hsg3+GYC3EdG/g8zFPuMIg7wiNMJtOFE0ImzYFlKsfuuwsEeYeWWJ1g3Fsp4H4J8z8y8Q0UsA/DiAlbVejhONcBsaGk4JpOfDsXzT+mJZrwHwT/Xh6wC86lhWugVa9Y2GhoZTAkbmfqvbFeJjAJ6ty18K4KEr/cJt0RRuQ0PDKcJOJs2+BcCPElEH4ADarGAXaITb0NBwerCDKAVmfjeAp534iibQCLehoeGUoGWaNTQ0NOwQjXAbGhoadgA+98VriPl8Z+EQ0Z8C+F8n8NXXA3jkBL73pHCWxnuWxgqcrfGe1Fg/l5k/60q+gIjeChnfNniEmV9wJeu7Gjj3hHtSIKIH1wVenzacpfGepbECZ2u8Z2ms5xEtDrehoaFhR2iE29DQ0LAjNMK9fGxbKOO04CyN9yyNFThb4z1LYz13aB5uQ0NDw47QFG5DQ0PDjtAIt6GhoWFHaIS7JYjocUT0K0T0kN5/xpr3PpaI/oSI/sMuxzgaw8bxEtHNRPReIvpdIvotIvr7Ox7jC4joD4joQ0T0yonX94no5/T1B4jopl2ObzSWTWP9TiL6Pd2O9xDR516NcbrxrB2ve9+LiYiJqIWK7QCNcLfHKwHcw8xPBnCPPl6F7wPwjp2MajW2Ge9FAN/IzH8DwAsA/AgR/eVdDI6IIoD/CODLADwFwNcR0VNGb/smAJ9k5s8H8MMA/u0uxjbGlmP97wAuMPNTAfw8pE/WVcGW4wURfRqAfwLggd2O8NpFI9ztcQeAV+vyqwF85dSbiOhpAJ4A4L/taFyrsHG8zPyHzPyQLn8MwCcAXFG20BFwC4APMfP/YOY5gJ+FjNnD/4afB/A8IlrfAO1ksHGszPx2Zr6oD+8H8MQdj9Fjm20LiDD4AUiJwoYdoBHu9ngCMz8MAHr/+PEbiCgA+CEA/2LHY5vCxvF6ENEtAPYAfHgHYwOAvwLgo+7xH+tzk+9h5h7AnwP4zJ2MbsU4FFNj9fgmAL98oiNaj43jJaK/BeBJzPzmXQ7sWkcrXuOwoQ/SNng5gF9i5o/uQogdw3jte24A8FMA7uTdVQ+Z2kDjGMVt3rMLbD0OIvoGABdQOwpcDawdrwqDHwbw0l0NqEHQCNdhQx+kjxPRDcz8sBLUJybe9nQAtxLRywF8KoA9InqUmdf5vVdzvCCixwJ4C4DvYeb7T2KcK/DHAJ7kHj8Ry+2q7T1/rNX5Px3A/9nN8CbHYZgaK4jo+ZCT3bOZ+XBHY5vCpvF+GoAvAnCvCoPPBvAmInoRM/tOuA3HjGYpbI83AbhTl+8E8MbxG5j565n5Rma+CcArALzmpMh2C2wcLxHtAXgDZJyv2+HYAOB9AJ5MRJ+n4/gHkDF7+N/wYgC/xlcnU2fjWPUS/T8DeBEzT57cdoi142XmP2fm65n5Jv2v3g8ZdyPbE0Yj3O1xF4DbieghALfrYxDRBSLaWdfPI2Cb8b4EwG0AXkpEH9DbzbsYnHqy3w7gbQB+H8Brmfl3ieh7iehF+rYfB/CZRPQhAN+J9ZEhV3usPwi5qnmdbsfxyWNn2HK8DVcBLbW3oaGhYUdoCrehoaFhR2iE29DQ0LAjNMJtaGho2BEa4TY0NDTsCI1wGxoaGnaERrgNDQ0NO0Ij3IaGhoYdoRFuw0oQUdIg/t8hotcR0adc5ve8xy0/OvH6FxLRR7Qu68eI6N9PvOcmIrpERB+4nDFsGN9j9HfOiej64/7+hgZDI9yGdbjEzDcz8xcBmAP41sv5EmZ+xoa3fBJSI+EpzPw5zPwdK973YWY+9kw4Zr6k37tUH6Gh4TjRCLdhW7wLwOer0vwde5KIXkFE/1qf/yARvVq7Hvy8KeIpVTvC9ZBqVv9328G49b1KFfh/IaLnE9F9JF0ubtnmPZezIRoaLheNcBs2Qit1fRmA397w1i8EcLd2Pfh/kHKVG8HMvwfg9ZCqYB+f6k6wAp8P4EcBPBXAXwPwDwE8C1I46LuP8J6Ghp2gEW7DOjxGPdMHAfwRpJjMOnyUme/T5Z+GENtGaKWtfwzgmcz8BCXgbfA/mfm3tYbv70JaCjHkxHDTEd7T0LATtHq4DetwaeyZElGP4Yn6Orc8roS0bWWk5wP4mcuox+trzmb3OKP+t7d5T0PDTtAUbsNR8XEAjyeizySifQBf7l67kYierstfB+DdW37nHwB4BhE95hjH2dBw6tAIt+FIYOYFgO+FdHp9M4APupd/H8CdRPRbAB4H4Me2/M43AXg7gAeI6H1EtCpKoaHhTKPVw204FhDRTQDerCFkV/I9fwnAB5n5SaPnj+X7N6z7I5BW54+c1Doarm00hdtw2vAtWG61AwAJwKefZOIDgBnE221oOBE0hdvQ0NCwIzSF29DQ0LAjNMJtaGho2BEa4TY0NDTsCI1wGxoaGnaERrgNDQ0NO0Ij3IaGhoYdoRFuQ0NDw47w/wHwNXjy60t79wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "z = NollZernike(Z25=1e-6, wavelength=telecon_wvl, z_unit='waves')\n", "z.plot2d() # units visible on plot" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The wavelength can be checked via its `represents` property if you lose track of what it is:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\mathrm{1550\\,nm}$" ], "text/plain": [ "Unit(\"1550 nm\")" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "z.wavelength.represents" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Turning our attention to the labels on this plot, those too can be adjusted by creating a `Labels` object and passing it with the `labels` keyword argument:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEOCAYAAACUxJyzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO29e9QlRXX3/9mCaDSoIKhEJGgkGi9odAT8hZCooJj4egWUoIDBxTLGpb+44i/6GqNBTVBjjPH2iogBvETACyRRCaCYeBl0EBQRDGiCDmIicskYXpSZ2b8/unumnj17V1ef55zznEt9Z531nK6u2rWra+p7vrW7ulpUlYqKioqK+cWd1tqBioqKiorVoRJ5RUVFxZyjEnlFRUXFnKMSeUVFRcWcoxJ5RUVFxZyjEnlFRUXFnKMSeUVFxcJBRA4Xke+IyLUi8irn/CEi8nUR2SwiR5hzx4nINe3nuCT9sSJyRWvzb0VEptGWElQir6ioWCiIyE7Au4GnAg8DjhaRh5ls3weOBz5iyu4OvA44EDgAeJ2I7Naefi9wIrBf+zl8Qk0YjErkFRUVi4YDgGtV9Xuq+nPg74FnpBlU9T9U9ZvAVlP2KcAFqnqTqt4MXAAcLiJ7AfdQ1a9o8xTlGcAzJ96SQuy81g5MGnvssYfuu+++a+1GRcVC49JLL71RVfdcjY2nHL6//uTGnxbW9+9XArcnSaeo6int9/sDP0jObaRR2CXwyt6//Wx00mcCC0/k++67Lxs2bFhrNyoqFhoict1qbfzkxp9yyYY3FOXdWZ5/u6qui9xx0kr3IonKrsbmxFEUWhGRu4nIo0zaPiIyM79IFRUV8w1F2bp1S9GnBxuBByTHewM/LHQjKrux/T6KzYmjNEZ+B/AJEbl7knYqsNf4XaqoqFhKqLJ168+KPj34GrCfiDxQRHYBngecV+jF+cCTRWS39ibnk4HzVfUGYJOIHNSuVjkWOHe0ho4fRUSuqncAnwSeC40aB/ZU1RqzqKioGAsUZatuLvpk7ahuBl5KQ8pXAWep6pUicpKIPB1ARB4nIhuBI4H3iciVbdmbgDfQ/Bh8DTipTQP4AxoBey3wXeAz474Go2JIjPxU4P3AaTS/Rh+ciEcVFRVLCkV7SLrYkuqngU+btD9Lvn+NlaGSNN9pNDxn0zcAjxiLg2NGMZGr6tUigoj8KnA0cPDk3KqoqFg+jI/Ilw1DV618gEaZf7NdY1lRUVExHqiiWyuRj4KhRH4W8A7gpAn4UrGAuODAo3ZIE1FUZdv3FDY9PfbKqAoiyqHrzx6/8xXTR1XkI2EQkavqbcA9J+RLxRziggOPcom1Q7obRUfETXpDzCWEHpF7mnbhQUeGtoBK9HOBreiW/7vWTswlFv6BoIrJIiVXEV1BvN35lNxXkvyOPwA2LbWRlrH15mxWzAdUa4x8VBQTuYisA14D/HJbTgBV1f0n5FvFjMFT37CSgD1V7JFpSrLd94h4bV6bltYT/Wh0vqd1HXbJWYOvQcUkoVBj5CNhiCL/MPBK4Ap23GimYglh1XeXBj45e+ETL1xibdp0q9BLfghsfRUzCK1EPiqGEPmPVbX06aiKOUaqXsEn7BSeuk7JOyXaLs1T0115e2xJOvoekXqaN/XP3ojt8tV4+hqihlZGwhAif52InApcBGx7RlZVPzF2ryqmhgsPOhKwBLrjzUT7PU3zSD4XF8/FtaMYd0TsfeEVj+htXdb/9JpUcp8eRLcim2/vz1ixA4YQ+QuBhwJ3ZntoRYFK5HOMXEzbI+A03d7gtN/TOrwboWn9HnJhk1woJWqXNzOwbbP+9vlYMU7U0MqoGELkj1LVR07Mk4qpIb1p2S0P7IsfezHpFBHJegRp7XnH1lZaPhdD936EvB+Zkh+Q1E4XhhGpa9YnB0VqaGUkDHlD0HrndUkVCwIbR07J056PVHf6PRdCKVG4tv6uvI2z55Y6ltwctf7aOHpqu94onTAU2Lql7FOxAkMU+cHAcSLy7zQx8rr8cI7QxX1hpQpvjv012TkF3ZfPi0WneS2JRgq/b+VL7qZpFEuPQjzeD5On5NNrWdX5OKFIDa2MhCFEPvEXjYrI4TRbAOwEnKqqJwf5jgDOBh5Xt9ItQ184pO/YU6MRSabn079p3lz8O4ple22wNqIfhpwij27i5spVdT4JaFXbI2LI7oerfpVTDsmbrw+jeRvH10TkPFX9tsm3K/Ay4JJJ+jPv8JYQWuTCFjlF7aV5xJe7WWoVefq9JHSTQ0TuXttz53I/Ll1aGjuHqtBXBVVkc+9LIyoczNKTndvefN3W1735+tsm3xuAtwB/PKZ6FwZe+KT5HpN4pHYjou7ODVGno6joiNxzK0+i9nix+Vy8PnfDNfpx7NK6PV/qTdERoFWRj4pZerKz983XIvLrwANU9R9FJCRyETkROBFgn332mYCrs42IsDzy8hRx34oO76+XJ3fj0RJtVHf0I5HaierM+R9dj6ge2wbvuO/aVfRDKpGPhFl6sjP7lmoRuRPwduD4PkOqegpwCsC6desWPphp9xGBflUJ/QSWIhefjsrlfiCsP7lwSi48lNpK7Vj/oh+IXOglapv3Q+mFXboZUlXmpaiKfFTM0pOdfW++3pXmNUsXN+8+5X7AeSLy9GW/4TlEEeZCKbkYeanizdXp+Vuary8EEoV+vPO5czn/vRh+FMIZeo0qQFTHpsj7Fk6IyF2AM4DHAj8Bnquq/yEix9BEHjrsDzxGVS8XkYtpXjjf7bX7ZFX9r7E4vErM0pOd2958DVxP8+br3+tOquqtwB7dcXtR/3hZSTx9qMeLOVt44RRvRUaHyF5Up60jCj9EPnlK3NoovTHp+Zk759nI3QfoC0FZv9O90uuOixmoIpt/vmozhQsnTgBuVtUHi8jzgDfTkPmHacLIiMgjgXNV9fKk3DGzyDkz82Snqm4Wke7N1zsBp3VvvgY21A27GkRv3IkQKdmIAL3Yda6OXJjGC2fkfgCi8EhJ20pj3l44x7Yv+tGLwkO5G7BpHB+o4ZY+jEeRlyyceAbw+vb7OcC7RERUNf2PfjTw0XE4NGkMIfL1IvIwuxxwnOh787VJ/+1J+THriMIEuSm+t0LEu1lo7eaUdbRyxPppidUjbls++uEoCWXkftg8cvWIuu9HwfvRia5hFPuvsFBk61jWUfQunEjztCLyVuDewI1JnufSEH6KD4rIFuDjwBsN8a8Z6pOdc4JOyYnhgChkkRJkmm7P94UHPHix9OimoVcm8jciPltXWqYE0czA8zuKh3t2cm33bKb1V2XuQBmiyPcQkTTEcUq7yAF6Fk6U5BGRA4HbVPVbyfljVPX69lmWjwMvoImzrzlm6snOijxKlF8OuXKRavTsRmo1ipVbMuz7kfBmDta+zROV99oetS0K++RCMH2zg5IQVkWHQatWblTVdcG5voUTaZ6NIrIzzbuIb0rOPw8TVlHV69u/m0TkIzQhnPki8u7JzvbBoG+q6urvSlRkkT7gAzEZeWSSIqd0+25y9sWeo3j3EOK07fOI1CKKWUd+pmWitkUkXWI/CjP19Vm693m9EQqiYwmtZBdOtDgPOA74CnAE8LkuTNIudT4SOGSbXw3Z30tVbxSROwNPAy4ch7PjwKCXL4vIXsCXgd8HPjQRjyq2IacQbR7vu0WkmL1yHuFGYQdPtXt12HNR6Cfyv8Rvmz/6kbLXIzqXlrU2baw9slXaZ0sPVdh8xxjMFC2c+ABwpohcS6PEn5eYOATY2N0sbXEX4PyWxHeiIfH3r9rZMWEQkdP8gp0OvIhK5BND+ph3ChueyJFA7nwpaaV15lB6s68k7BGFNHI/DLZNfTOKKHRSEgLqCzWtps+6fl9aZa4K47nZSd/CCVW9nUZ1e2UvBg4yaf9Ds+Z8JjGUyF8A/BbNgzi/oqrfnYBPSwvvCU3wlwHmpv4RmfQpRS9PjphsvVE9UQjElilpW6qCc0o2InPPX68eL2/Oxjj7bJkJvT6iPxqKXywhIk8ArlbVG4EP0iyorxgzcgoyJaYu3R6nfzvkQgh9+b16vHx9arTLH9mwdmx9EYH3hVssuZbGx73r4vlky4+7z5YLrSIv+VSswJA3BJ1AE1cC+BhwZHtToGIM2L680J/Op99TFWePvXIRaXnlIuVs64uUvM2bElpf+0rabtuQwrbNmxHkZieRLa/tuVmL5/cofWZvdi88lErkI6KIiEXkXjQxo88AqOp/A+uB35mca8uHSOF6Ci79bpVi3/Q9pwY9hdyl5+rzwhTWTqRMPX8iorY2ox8rW0dpyKr7nrtWHolPqs+WC1WRj4qiGLmq3gI82KS9YCIeLRnsnikdPCVqv6fwiCqHiIijeHdKhp4Nz75HiH2+55RulM/67c0ObHl7zvPRy5e7BpPos2WKl4sqMoZVK8uIoTc7K8YE+2aZ7ntf7DaKQ5fCI6KI6PoI1SuT2oxIt4Rk+wjXI24vZNH3Y+O1M7omkXqfRp91/18WntCr2h4JlcjXCDlC9QgjPfYIMhfv9WzY4z7i9Ug9IjmPtD1Cs75HZNmnvm1dab5I9Q5JS8/NQp8tLLoYecVg1JuVa4B0B0Mv5hoN3BIFaZFT1Wn9Xp5cnDdH/DmfLGl7JB2FRrzv9gfEI/YhxJ+2yap8e13Wos+83S8XBzVGPiqGvLNTgGOAB6nqSSKyD3A/Vf3qxLxbUHgD1SMuL3zhpfURQmcrF8boK2/9yf3NtclT3X2hGc9eV75k5mB98NR1rt6ojWvZZwsJBbYucPsmiCGhlffQvFDiicBJwCaaHcAeNwG/FhLbn9jsn4qnx32q2JbrU6ApkUTq2pKUJSOLiICsz14bU5SGdbxrkftB8a5PdJ0i8p2VPutukC/ezokKmzevtRNziSFEfqCqPkZELgNQ1ZtFZJcJ+bVQ2P6GmB3XMnvH3XdPgXaIYsle3j7lHSm9iKBKlHyOrDx/c+mltkpVa3SdPFuWkGepzxZuK9yqyEfGECK/o32FkgKIyJ5sf+VbRQ9yU29PwUUKMwpL9IUhvDIWuVh3Lm5b0ua+ME1f/en1iXywPzK5cE5OoXuKPbUza322UBjP7odLhyE3O/8W+CRwHxF5E/BF4C8m4tUCIX1fI/TfGLNhD++vRwjpea+eTl1GIYk+JRuVtaq1L2zghWwitZsL89hPes7WH7UzCn/0XUt7Xda6zxbnCVBtFHnJp2IFhuxH/mERuRR4Ek2M4JmqetXEPFsQeIppCIFa5JRkarsv1OCpV2+aH4Uz+hS6pzZLfLe2cgo81ybvOPIvdx2srVnus7lHDa2MjEHryFX1auDqCfmyUIi2orXwVKQXzsiFHbrvqY2+sqkNa8ue6wtd9JFcX7n0WkTIlbPlS0NHlgxLSXoW+2xh4uWVyEfCkN0PT2/3XOmOdxOR0ybj1vwjNxA9JZWLf5aSS8m50tCBR85eKMXa7yOhKPxhfbK+55S+57cta0NBnlJfhD6bZ6iCbtaizyJARO4mIo8yafuIyP2H2hoSI9+/3XMFAFW9Gfj1oRUuC3KDLKdQc/HhaMofxZhzZTyC9HwoIW0PfaRScj6XJyJka9va8QgwPZezMw99NtdQmuUTJZ/FwB3AJ0Tk7knaqcBeQw0NCa3cSUR2awkcEdl9YPmlQLoJFsRTaU9J5uARaukUPPXDI5pcDLokpGHPRSEXa9v62Rfz9UIxUazYq680lDPPfTb3m2wtDkn3QlXvEJFPAs8FTmsfstxTVTcMtTWEiN8GfFlEzmmPjwTeNLTCZYAlNasII9hBXlouR1CWXD2fSuvxykftjNIi8vP8TFV3Lt48CtnasoveZ3ODBWpKIU6leffnacCxNC/tGYwhq1bOEJENNE92CvBsVf32KJUuInJ7YHhEZKfm0eCNkKo1L92W94g3ze8pvpwKtGo4UptDFKiHqH3RtYxmD5FfuXrnuc/mUpkr6NbxhIhE5HDgHTQvSj5VVU825+8CnEHzHs6fAM9V1f8QkX2Bq4DvtFnXq+qL2zKPBf4O+AWa94G+XFVX9dOjqleLCCLyq8DRwMGj2Bm6adYNwFeBbwB7iMgho1S6iOgGaF/80lOt3uAtKZc77+Xv899Ttjml3Fc2FweO6u5sRz8i0bWK2m1h/V6WPpsbjCFG3j64+G7gqcDDgKNF5GEm2wnAzar6YODtwJuTc99V1Ue3nxcn6e8FTgT2az+Hj9BCDx+gUebf7ELXQzFk1cqLgH8Bzgf+vP37+lEqXTSs5oEMj/Qi1ekRaFfWDt6+MIQNW/SRr1d3aejAlo3aa23bUEfJD4GnXj0FuxrMU5/N1W6JCrr5TkWfHhwAXKuq31PVnwN/DzzD5HkGcHr7/RzgSe3GgC5EZC/gHqr6lVaFnwE8c5RmOjgLeBTbX6U5GEMU+ctpNsi6TlWfQLNi5cejVrxoSJWkR3BRiCE3hbf2PfWY2sihJBTikXMunJKSSeSzteuRztC227x9KjjKsyx9Nj8Q2Fr4aSICG5LPiYmh+wM/SI43tml4eVR1M3ArcO/23ANF5DIR+YKI/GaSf2OPzZGgqrep6j1V9cJRbQy52Xm7qt7exnPu0sZ2HjJqxYsEb4Dav925Lv8QhWjVra03IsUoj2fX8zPXvshH60POp1z7+mxFPxJ9JOtd+0Xvs7lCud83quq64JxnxP6iRXluAPZR1Z+0MfFPicjDC22uGYYQ+cb2gaBPAReIyM3ADyfj1nygW2qYbkvroUTlQtlKhFy6nXpHajA99gjKI0KrXEvI2irIiCBz7fPqtCSZU8wRAfeR8iL12dzc+Bzfzc6NwAOS473Zkau6PBtFZGfgnsBNbdjkZwCqeqmIfBf41Tb/3j02B0NE1gGvAX6Zho+lqVr3H2JnyKqVZ7VfXy8in6dp+GeHVLZo8JSd9zdSljmy8OqxNmwdfTY85MrlSDgi9SjNs1PS7ghRqCHyIyL1ZeqzucDWsby07GvAfiLyQOB64HnA75k85wHHAV8BjgA+p6ra7up6k6puEZEH0dzU/J6q3iQim0TkIOASmqWC7xyDrx8GXglcwSpW0Q96oEdEdqNpmAK3AAfR3ABdKqQviCidbvep1jQNYuLpkBvkns2I8Gy9ke1cWc/3yF9PTeeuTWQz9+PRF95Z1j6b+RdSqKBbVk/kqrpZRF5KsyBjJ+A0Vb1SRE4CNqjqeTQ3Fs8UkWuBm2jIHuAQ4CQR2QxsAV6sqje15/6A7csPP9N+Vosft/6sCkNe9fYimhueewOX05D4V2jWlS8VIjKK4KnAaICWElRfeMGzWVpfVNbGpaOQgEVfWz1SzPmUQ1+oJc2zrH020xiPIkdVP02z1jtN+7Pk++00DzXach+nefOZZ3MD8IixOLgdrxORU4GLaEM6bV2fGGJkiCLvVq2sV9UniMhDaZYhLh3s1LhL8+Ky3sC2x1EZm7dUyeb8TNO9c169to2WuHPhA+9733Xy4s+eHzkizF27Ib54thahz2YROr4Y+TzhhcBDgTuzPbSiwMSIfOKrVgqexnoF8CJgM83Sx99X1evG6UMfmulp/6qExt9+orF5I4Vr8+di1FFYIiVIj2ijwe6RX5ru1e/5HLU1ytOnYL3zns0c4Xrti7AofXbBgUfN6I3PbUsLlwmPUtVHrtbIkHmMXbVyLmNctVL4NNZlwLr2ju45wFvGVX8pvEFnz0dkkp6PEJFNVC4ipK6MZ7ckv+dz3zTfO+d9OoUYkVCO7LqPd96r11PPkZ+1z9YeaR/nPguE9Q7PDcYsrVrZ9jQWgIh0T2Nt289FVT+f5F8PPH+M9RcjRwjpQIrIww7QKLTgqdFI8UV12bRoEEdTc8/vSBmX2Ij8KkGp7971qn02u+S9DcrYYuRzhIOB40Tk32li5JNdfggw4VUr3tNYB2byn8B47hoPhhcHtWRgz6dl+wZVpFCtD6Vqy6rFPr89ePXljj1yH/W69F0zzw+PeG3dtc9mDcLWMaxamTOMZb+WWVq14v2PdP/XicjzgXXAbwXnT6TZ3IZ99tlnLM71LTlMyaPvfIchU+dcfNfLk+az5BQpuYiIh6R7tnL2S/32Yr59IZPUdu2zlendHiwzFStfTkX+I+A5wL6s5OOThhiZpb1WSp7GQkQOpXkS6umq+jN7HkBVT1HVdaq6bs899xyLc9GU2JuCe7FVLy6aIyqvrLVtp+9eXZ6faV7vvEeSnlL17OUINf1Y0s+VixARaI48a5/taG+WoFul6LNAOJcmhLwZ+J/kMwiztGql92ksEfl14H3A4ar6X2OsexC8AZWqunSQR9P3aJCmsOSRm1J7x546i+qI7EdEGtnsU8nWdo58o3qsos213dqpfTa7UPI/3AuKvVV11eGVmdlrpfBprLcCvwic3e44+X1Vffq4fPBg91Ppm8pbWDJI08Enrpya8tKjgR+FASKS8VRyLgSQs+O12VPjlrxyPxBRmyw5W39qn8V2ZmofFpVlDK18WUQeqapXrMbILK1aoeBprEPHWd9Q5AaIzeMRQApPEUb1eXXafBHB2bI5QvSIwbPl1Rn530dwHrmn/kR+eHXl0muflfXZWmPBwiYlOBg4fpFWrcws+pRPibLy1GCknnIE4cVSPZ/S9Gjw903hre+pLeuvne57be+r0yPXEnjlap8N67OZgAq6Zae19mLaeOo4jMzSqpWZQhdSgWHKKxciyJGMl57a7/LZAdg33U/t5hRdnwLN+Wevh83r+ZgLsXhtS/3PXY/IH5tW+2zl9ZiVEMuyKfLuyfR2O9tvtm80GoxZWrUyc4gUWhcO6JCGCCKUhhGistH0OVdvpBY9YrB2rE3vnHc9bBttu0tCDRHRWf/7FHPko/UrwrL22VpBWdmu3GeR0L5G7svAyO/lG0Lkt7c7htGtWgEW9g1BOdLpzvele+Rh06PB2Pef1RvwVtGl9iO7XTtzxJPzIQoNeHWk/lhlndorIeiSNnj5+tJrn60hlGVcfgjNvuin0+wjNRJmZtXKrMMLFaTnbHpEPF6IIC2b2sqFGWw9liCjMIDXLs++JRGvfq+Ntn19IY5ceCUil9LrXftseJ+tLQTVpVu1AvACmocbzxORX1HV7w41MFOrVmYZffHIaBD12UzLeYrSDsDo2Nbp+eaFHLz8Uf2pz54v9nsfkfXVbev1kCtb+2x4n601xvFiiXmCiDwBuFpVbxSRD9JsPfK/h9oZ6aqp6hdU9bxRA/OzjgsPavab74vHeaqpj5C6PF6Z3BTdI5DUdh8JeoM7UqteGZs/p0y776nijNpi6/eUaR9R2jbXPhu9z7pH99cEyxlaOYHmbUUAHwOOFJHBvDzq8sO7dmmqunDLD3NTYU8t2UHSlbF5PERKtk8x9RGhFwrw/OjzPSKiCFF7rdrMKeLc9fds95WpfZZHX3unBV2y0Eobqj6IJrSCqv63iKwHfgf4xyG2iq9au/zwX2ievPzz9u/rh1Q2b+iIp0RJeVPgblBZlWgHfAqPrDyysPWVTMVLCMAj3IiU0/IlpGGvQ9quPmUdkaTXxtpn4+uzaWNcilxEDheR74jItSLyKuf8XUTkY+35S0Rk3zb9MBG5VESuaP8+MSlzcWvz8vZzn1W1VfUWVX2wqmqS9gJVHUTiUJcfhkgHsw0RRIPCywc7TsOHKsQUnjqzU+p0MHv5U/veVNzz09pI68mFF2y6R2jWn+iaeCTlkWPts9X32ZpA/T70PjkUvqTmBOBmVX0w8HbgzW36jcD/at/acxxwpil3jKo+uv2s2X5PFnX5oYN0QNlpZ45IcorMqh+PgLxBmcIOvkhNRdNzbxBYv0qUp6fkPAKy5yNEitaSW47gap+Nv8/WAuMgcpKX1LT38bqX1KR4Bs2SP2jeNvYkERFVvUxVu9V4VwJ3FZG7jKl5E0NdfujA+49iVWWqbOxg7vLbfPac993Wmdqw+fv89hRuSRsjFeopUlvWU7bpdbNtyKlDr+1evZEKrX02Wp+tFRRha/kj+nuIyIbk+BRVPaX9XvKSmm152g37bgXuTaPIOzwHuMxsl/1BEdkCfBx4YxoWWUvU5YcJ0jv20X9uSzwlJGPTowHuTbVz8IipD5bEPJXmTb9T+x5ZpfZtW2x6rk4LT6X3/WjUPltdn63ZSyeUIStSblTVdcE5z4i90Nk8IvJwmnDLk5Pzx6jq9SKyKw2RvwA4o9ThSWLIzc67isgrROQTwMuAXxlSfl7g/efPDfQcidiBFSm8tFxOTdl0W1869bQEY8t7ZBep1hRe+9PjaPqe2o/q99oVtdWrI/Wv9hkr6hilz9YC6fXIfXpQ8pKabXlEZGcaYXpTe7w38Eng2PThHFW9vv27CfgITQhnVZAGzxeRP2uP9xGRwXaHEPEZwMOBdwLvAn6NHW8EzD2iQd0dR+dsejQl9ga9l8/WF+WNlFrJtDk6lxKDp0Cja2JteGScy5/zJ6dga59Nts+miTER+baX1IjILjQvqTnP5DmP5mYmwBHA51RV2/DxPwGvVtUvdZlFZGcR2aP9fmfgacC3Vt1geA/weODo9ngTzY3aQRgSI3+Iqj4qOf68iHxjaIWzjGgg5BRd33/8vrCBN3htmTTN5vXKeGXT/LmpuT3v2SmZrnuEY8m4pE57nWye2meT7bOpQsfzsE/hS2o+AJwpItfSKPHntcVfCjwYeK2IvLZNezLN69fOb0l8J+BC4P2rdhYOVNXHiMhlre83tz8+gzCEyC8TkYNUdT2AiBwIfKmnzNyhRJV5asebPuempzlS6yOaSL32TeHT9nmqLcoblctN1z0/I+Lw2pbzL0ditc/G12fThgJbt45nP/KCl9TcDhzplHsj8MbA7GPH4txK3NEul1QAEdkT2DrUyBAiPxA4VkS+3x7vA1wlIlcAg99oMUtIH8kHX/146FNg6cCPprvewIr88OqOlGrp9NoSZ6TK7Hl7vaLr4JWzZSMSyqnUEsWau26er7XPVvZZNy4OXX922KZxY+tazQbWDn9LE4+/j4i8iSbM89p8kR0xhMhX/YLQeYA3QDp4AzCa/nuDLVJllgRy0+A0LRqAuemz9S/KH5GTZzNHxLm6PWXa13Zbjy1b+2wyfTYV6MLto9ILVf2wiFwKPAkQ4JmqetVQO0NuduwfK+4AACAASURBVB4A3NS+0eIFNE9D3VtVr+vecrEI6PtP7w2yPvUWDY7cgI3S0kEbqVD7Ny1vycQjAtteTxVaHyJ1mpJZ1K4cYeYUrrVh06LrZP310mqfxX02KWhSZ99nUSAipwM/UtV3q+q7gB+JyGlD7Qwh8teq6iYRORh4Cs1TUe8dWuEsIzfoo+OSc96g80grLWcHq1eXzW//k/eRkUcaXnqOTHOqOj1vCcTWE12PEpKqfTb5PpsWlo3Igf1V9ZbuQFVvptn+ZBCGEPmW9u/vAu9V1XOBwXdXZxHRYMqprFw+O3VO/6ZkFA0gbyBa5WXVWmo7gj2fq8+bekc2LdIBF9mNrmPnY0owHknVPpten00TS0jkd2p3lQVARHZn4K60DCxwvYi8DzgUeHO7/8BCPBAUDUBPrUTl+8iixKZVgda3KAxg67TqzVO5JSRQQja5QZXWH5FxRMqWzCPb6XHts8n12TSgKmwpf0R/UfA24Msick57fCTwpqFGhhD5UTQ3PP9KVW9pXxj6yqEVziLS//w5pdjlSf92yB3bQdE3DY4GaG5qb217dfal2fZGg9lTqtaHSHlb/yKFa+vLhRRqn022z6aJBVPbvVDVM9o9Y54ICPBsVf32UDtD9lq5DfhEcnwDcMPQCmcVkfoqJaJoWusRSmrXlo1IIjeoc+2J2paq3rSM1wbbjhysQvVIKkfiEbna6xu1J83bV2/tsx1tRH02LSwbkbe4AfgqzQt79hCRQ4a+sGchQiOrQbdWFvqVW6RkcuUtEdgwQponVUSRsvQGs4dolhEpSe97REDe9YkGoGc/yl+iCqNrV/tscn2WjpGJQpt15CWfRcG4Xtiz9EQeTVPtwLToBkL3SY9TG+l3bzB7AzbK4/nqEZZ3zvpk7Xn+RCTm1Wevh6cs7TXquwbpdevSIvVd+4wV+cfdZ9OAUnajc8FU+1he2DP47uiioU+deerHG+B2YHhEZAdHTjHl/ItIIjc1jhRrlCdHFt5x5HOfXZvmqcqcjdpna9Nnk8KCkXQJblfV20WE7oU9IjL4hT3FRC4iAhwDPEhVTxKRfYD7qepXh1Y6S7DE4amt9DhSRJ7K8sggyh8d5+rw1Js36DzCSNtqB483tc9N96M6I5ved69tJdfC1lX7bPvxOPtsWtiydemCBGN5Yc8QRf4ems1cngicRLPd4sdppgVzDe8/cJruDbocEeQGUG6qnfOtpExEDqk/1k/bfm/wegTs+RCdt9N2L29JyMC7Dl5d0bWofba6Pps0VJdPkY/rhT2DNs0ax3aLs4icIvLyeQQUDdI+wkm/R0RmVWqEHBlE36393F/rY5Q3d208UipR0X2E5LU/zVf7bGV9o/TZ5LFYNzJL0T4QtB+gwC3AQTQ3QIsxhMjHst3irMEqoug/sjctjQaDZz+11aWXDJCI7HIDMjdAPR+ivKkizJGrV29UV9QOmzeXp/bZ2vbZJLFsirxdtfJymrcYXU5D4l+hiXwUY0hAqttu8b7tdotfAv5ySGWziL7pPfgD0htI0blcujfobZqdYts8nlK0dVgCyKlAW78d0F4d6VTcHluCiZBT6bYdaT21z9auz8aNzoe+zwJhuqtWzHaLAE9X1auHVjhr8AZgTh3lps42T5cW1ZWWSfNae5GKiwZ46l+UL1Ks3vXJqVvrt9dez04JGfe1v/bZdPts0lBdypud01m1IiKbYMUbqCU5p6p6j6GVzhJyU88ceXjIlfGIIhrc0VQ3R0y2Tda+bWsfeUUklOaxbY98suTm+daHnP3aZ9Pps8ljKWPkY1m10vvzp6q7quo9ks+uyWesJC4ih4vId0TkWhF5lXP+LiLysfb8JSKy7+rr3D59tYMJfMWWpnffcwrTU0GRIkyPrT8eeUVtSuuL2mGVY+qX9cnaTst719D6k7OR+ph++tpR+2x6fTYNKDQPBRV8+rAaLhGRV7fp3xGRp5TaHKnNqs9S1VtU9fU0bwb6APDMoXZmZh7T3kh9N/BU4GHA0SLyMJPtBOBmVX0wzYst3jyeunfcTD83IDt4g7w7tuesuknJKhpIaT2RqixRdta/tH0eKdnrEpGiR7qR0u67ft418Oq0vllf0rpqn42/zyYN62P0yWE1XNLmex7wcJpNAt8jIjsV2hwJIrKbiBwAK1atDEIxkYvIXUXkFSLyCRH5uIj8kYjcdWiFGRwAXKuq31PVnwN/DzzD5HkGzQstAM4BntQ+qDQyogGZU242j7WTIyCvDquyPILzBrFXf27wRSrNO/YGjh340d/0u0ce0fX0roslnByJ1j6bXJ9NC2Paa2U1XPIM4O9V9Weq+u/Ata29EpuDsRZ7rZxB8yv1TuBdwK8BZw6tMIP7Az9Ijje2aW4eVd0M3Arc2xoSkRNFZIOIbPjxj8tuAHsqyFM3HlnYAWvzeYMltWFJIJ32euWs37aMJa5IwXiDNDfgcyo0UoupH7Z9Xt2WXCNyTs/XPpt8n00Hg/Za2aMb4+3nxMTQargkKlticxRMfa+Vh6jqo5Ljz4vIN4ZWmIH3v8WOhJI8qOopwCkA69atK5IT0bTYIwuvTHfsTU89JRYRW0RAkb2IfCI/PRUXEastY+30kXJ0jdK6re99yjUiytpnfvlx99kkMXDVyo2qui44txouidI9x8ZxYcayamWIIr9MRLbFbkTkQJq15OPCRuAByfHe7Hj3dlseEdmZ5nHWm1ZTaclgj8p5U3aPQKJjb3Dn8lvk1F9OiaV128EdtdGz7Sm/kmuX1l16vXLXp7Te2mej99m0sBUp+vRgNVwSlS2xOQrsqpVzR7FbsvzwCppfnjsDx4rI99tT+wCD32SRwdeA/UTkgcD1NDccfs/kOQ84jubJpyOAz6nq2P4XRkow/W6VmHc+tRVNW7060vN9iiuq19YR1dun1KIwhhf66Juie2o5Ck2UEKDnl5ev9tn4+mwa0PHVNTKXiMh5wEdE5K+BX6J5dP6rNEq9z+ZgTHOvlacNNToKVHWziLyUJti/E3Caql4pIicBG1T1PJqlOWeKyLU0v57PG0/dO07d00GfUz1W4XgDOafCrL0uv6cSo//kQ+uzeSOS9eroU8bpuYhYSn3zyqTptp7aZ34d4+yzyWI868hXwyVtvrNoROpm4A9VdQuAZ3PVzq70+wujlu0lclW9rvuebO6Srla5bodCI0JVPw182qT9WfL9dpqXk04EkVrJDZQcOUWxXW967k3dc8qsT9GmJOWVs3Zz03RP+UX+RErP/rVttdess5tTnFE9tc8m12eTxrjqWg2XqOqbcF6A7NkcBzxeHfqqtyH7kY9lc5dZg0cmdqBYYvAGtRcysIhspMd96ik37Y2IwWtfVH+unkiZeW3PEZIloMg/r47Il9pnk+2zaUBZgF34BmItNs0ayzKZWUPfVNgbpH3T01y6HfheOU8ReYPd89uShTfwc9Pw1L7np0dIXtu6716arSc9jq5xjixrn02vzyYKbVatlHwWCFNffjiWZTKzjGiqn1OQJVPxEpWU2vHIp0TFpvYjdeapOs9OVJ9Hkn0hEouScp56zRGdtVH7bLx9Ni2UPH6/YJjuq94Y0+Yus4aSwZ47l05Jo8E2lGgsosFsB7AlA48Q7PEQpeoRSEm4wSufXrdcXanfHrnXPtuxXut3qV99fTZpaN00a3KvehORnVV187iWycwiIpVSqvyiqW806C3RRGX7lJRHLiVKOD3OKc8ShRqpyqiOvuvbp5yjenM2U9Q+W12fTRpbp1fVTGCayw+/CjzGVD7yMplZhDfg+/4zW7VkEamb0ulzpERzyjWq19qKVF16zmubpxq7cl5aybXIXYMob3Qtap9Nts+mgSUMrTCtV70t9JWNBkyOVLq/Q//je4MtV6dHRLbuPmWYCx+UkKG1G7Un8jFHrFE7IoUaKfTaZ9Prs0mieUR/oelmB4xr1UoJke8pIq+ITqrqXw+pcNZw6PqzufCgZjmp/Q+c+4+tuvIVXNHg6PJ2aTklFw0oT8F5Sis676nE1LccaaT2Ini+R9cxvQ7e9YoUbWnIoPbZ+PvssEvOCu2MGwWP3y8aulUr61X1CSLyUJpdEAehhMh3An6RBVbm9j9yhxyB2IGSU1te3hwxlCrZIcopR0Kp31G7ojoiwoh861Ox1iePjD1SjOzXPlt9n00LynRvrs4IprZq5QZVPWkEB+cG3iBJYdVRl5ZTXjZfn7Lqm3rn0vt8zdnN5feIxNqzJBYRTqSCc22KFLY9V3odap+N1mfTQ121MrFVKyywEu8wRJ14//lLyMjaSY+7833TYo/YPFKy/nv1leS318P64Nn2/Ext2DZHx9ZWWj5HmrXPVl4P64Nn2/PT2pgWpj8PWFtMc9XKk4YanTfkpuGlSiYarKl977hvulxiS9V/0UFqv4/APAXotTVSyCV19qnXFFHbI4WZa3vts9X32TSgsHSKvH3L2kuAg2kuwRcZ4RWcJZtmrWq/73lCnxryppveQMsNdI9AbJ196snWl7Np0yK/Ir8jxdyXLyIySxTRFD4i59x18dpS+2z1fTZNbFkyIqd589ommjevARxN8+a1QZsDDnmyc+ERKTX7nzqnwLwB4A0sbxD3TYUjf9I6SojBlik59uqMSDI9b9vl+RmRUR9p2by1zybTZ9OC6vIpcsb05rWF2n1mVKTTXO9c9zciCk/R5Gx7U2ZPsfX5U+q3RwieDe+757uX5p33iMleqz5StOeGtr322er7bJrQws8CYSxvXit5RH8T/rUTQFX1HkMrnTV062S79eR2KupNc71pqPefPpqGe9875FRe3zQ/qtOiT+1G0/BSH3L1l6ho61f6Q+Bdt9pn4++zQ9efnc0zCSyhIj+QHd+8dlX3ZjZV3b/ESEmMfNfRfZwvlEzlU3j5cgomJQGvrvTYU3+eQspN90t8tOl9KiwiTi9PTm1boo3qjq6/Vfi1z8bfZ9OGsnz7kQOHj8PIoBj5ON5kMeuwCjEXAkjL2LIpvLTcAI4IIc3j+VuSL6rbuw59bci1y9Zj22T98a5dRKCer7XPxt9n00ccdlpgHAB8VlU3icif0uxr9UZV/foQI8Ux8nZPgH+heWfdn7d/Xz+kslmHN6CitD6CyNXhDVCP6HIKyiOAPmWW5rXhDG8aHineoYiII62/L19kt/bZ9jyT7LNpQGlWrZR8VgMR2V1ELhCRa9q/uwX5jmvzXCMix7VpdxORfxKRq0XkShE5Ocl/vIj8WEQubz8vKnDntS2JHww8BTgdeO/QNi39G4JSpIMgRW5A2HOeDW96ntaZG8yp/ZyN1FbUNq/OqKwlGhsO8fy2RBJdj1x+r0zU7pz92meT6bNJY6uWfVaJVwEXqep+wEXt8QqIyO7A62hi2AcAr0sI/69U9aE0HPgbIvLUpOjHVPXR7efUAl+2tH9/F3ivqp4L7DK0QfUNQQm6m54XHHhUkZrxQgPdsc3TfbfT6g7p4Mupvxwx5Xy0+aNwhVWNJWGNqJ6IHG1IwRJTjkQj+7kwilem9lncZ9PcJGuH+qdTzTOA326/nw5cDPyJyfMU4ILuORoRuQA4XFU/CnweQFV/LiJfp9m5cFRcLyLvAw4F3iwid2GE1YRDCtg9Ac5lAd4Q5MGbktswQJrPTmctCXllU3gDz5axg9UbjHbK7U3FvWm3pwo9wutTzvY6Re32ro/Xfg9efV7Z2mej99laoVtHXvIB9hCRDcnnxAFV3VdVb2jq1BuA+zh57g/8IDne2KZtQ8uH/4tG1Xd4joh8U0TOEZEHFPhyFE2Y+nBVvQXYHXhlcUtaFCvyYE+AzwytcN6Rm4JbpekNzr6ptrWRlk9tpvV7PnpElFO+FqUklvM7KhP566n11J9Igfeh9tlofbYWGLBq5UZVXRedFJELgfs5p15TaN+7yNsulIjsDHwU+FtV/V6b/A/AR1X1ZyLyYhq1n91XXFVvAz6RHN8A3FDo4zYUE7mInA68XFVvUdUvtPGi9wG/P7TSeUFELN40PCIaO+XPTecjMvHqivzMIVcuV5dHhGmeXLpHhn2Ek163SJHaPFF5r+21z/rT1wLK+B7RV9VDo3Mi8p8ispeq3iAiewH/5WTbyPbwCzThk4uT41OAa1T1b5I6f5Kcfz/w5hFcHwlDQiv7t9IfAFW9mSbYv3CwgzuXJ1Jw3gBL4U3x0/weSdjpuRea6FO80Tk7lY+m+l47vOvlhRIitZnaiurPtclew9pnq++ztYJq2WeVOA84rv1+HHCuk+d84MkislsrWp/cpiEib6SJSPy/aYH2R6HD04GrVu1pIYYQ+Z3SZTrtXd2l2aslNz3vjrt80bkob0QQ3pQ+mt573629iFBK6/SI2RKqV79H4Ll071pEfufIp/bZaH22ltiKFH1WiZOBw0TkGuCw9hgRWScipwK0NznfAHyt/ZykqjeJyN404ZmHAV83ywxf1i5J/AbwMuD4PkdEZAfV7qX1YQgRvw34soicQzMLOgp409AK5wGHXXLWtpUr0VQ2HTzRlN/Lk6JvChwpWFuPh9x5ey4irbTOHHKK216PvrBHRK6W3CKSrn22uj5bi8fyOyhjWVrYX08TAtlhe25V3QC8KDk+DTjN5NmIHz9HVV8NvHqgO4ex44qZpzppWQy52XmGiFwKPIGmIc9W1W8PqWxeERED5GOMfaousuGd9xRWZD8ik5KwRR9hRMSaqyendr1pfa5tVplGqH22o62SPltrjCFsMhcQkT+g2Yf8QSLyzeTUrkxi06wUqnolcOXQSuYRkbLrm7LagRYpsWi6nVOLtpynUiPFa9uQU5+5kIN3jXJTfnvNcoSXI5VIkXs/GLXPRu+ztcVYwibzgo/QrPr7S1Y+kLRplHdAlOx++EVVPdjZBVGAhdj90EOkcLxBlZvOW1vWXno+N62233PqLsprffPaVzp993yMFF+u7dYPSyo5crOofTaePlsrqMKWJVHkqnorcCvNiyRWjZLdDw9u/y7NLojQfyMsIoa0rP0eqSs7CPsUZORb32DOqbyIIPrabNvuwSPnXFqkKj1l65Fk7bN8m23bZwnTiJHPEtonOZ8D7EvCx0NfeD9kHfmbVfVP+tIWBd1NnwsPOtJVLdF0uEOfskxteGUiu5FqyqlDD1GoIiIijxz7yFsk3qgq8r1PXeZCEWm9tc+G9dla3uRMsWQ8Ds3Sx1uBS4GfjWpkSIx8LHdX5xElZJRTah6RRErQU29pHd4Uf2hb7MCOiK6PUL0yqc2IdEtIto9wPeJOfzhqnw3rs1mAspQvlthbVVe9J3lJjHysd1fnDel//mgADElLz0WD0yOM9NgjSK+8h74pfI6YbfmI5DzS9q6H9T0iyz71betK89U+K+uzWcGyrFpJ8GUReaSqXrEaIyWKfKx3Vz20Dxd9jCZO9B/AUe2To2meR9Ps03sPmq0f36SqHxtH/TnYHRFLlVduwHhT3RIlmVOdJdN3b2oeEbDnZ0RKOZKz7egjv0gV2/O2bbkfj9pncZ/NSkgFaB/RX2svpoPuVW40HPxCEfkeTWilW0RS9Iq3DiU3O8d6dzVAtz/wySLyqvbYhmxuA45V1WtE5JeAS0Xk/HTbgEmiRGmliJRSWjZSkx7ReQTjkVjkc246HqnBiCSs/157ojZFJOn5Wqq4ozZG7al9NlsqPMUSvertaeM0NuRm511pQiwH0/ySfJFmI/Tbx+BH7/7Aqvpvyfcfish/AXsCUyHy7XX70/yIlKKB3DcV98rnFG+fqk3Leyo1p669NlhEBJRTsh5KQwTetcj9oNQ+y/fZTECXZ9WKql4HICKvcE7fKiKXqurlpfaG3Ow8A9gEvLM9Pho4EzhygI0IK/YHFhFvf+BtEJEDaN6i8d0x1F2Ebgp6wYFHrUj3VJk9b6fE0bQ9IhPPvlcm8qVPxUXT/IigSlRhjqw8f3PppbYiNZ2zv6x9tpYvj4ig7WfJsK79/EN7/Ls0e7u8WETOVtW3lBgZQuQPUdVHJcefbzeHKcIY9gfu7OxF8wNynKq6M7F2k/kTAfbZZ58h5gvqz68B7gs7RNNsjzBySjRSZzaf9T0tY+Gl96ndHGnm6otILUfQUR6vHi8skeatfTabWBZFnuDewGNU9acAIvI64BzgEJoliWMn8stE5CBVXd9WeCADVq2MYX9gROQewD8Bf9r5EdR1Cs1+waxbt26s/zUisvAGp5c2SlghmpLnSK5E/UZtiXzL1W2/59RlRLC50EFqKyI8S7R97Vz2PptFLOGqlX2AnyfHdwC/rKr/V0SK15UPIfIDgWNF5PuJA1d1d1+H3mU16PYHPplgf2AR2QX4JHCGqs7ErfZI/aVpucFYGgrw6ky/9xFJn5KLCDOqLzqfHkcE1+e7tZVT4Lk2eceRf8veZ7MCBTbPpmuTxEeA9e2rM6F5ddxHReTuQPGmhEOI/NltfqVZ/nfbgLJ9OBk4S0ROAL5PG3cXkXXAi1X1RTTb5h4C3FtEjm/LHT/khsA4cNglZ2172jMajI3vZQM+p9qsco2m0N331EZf2dSGtWXP5dSs9bfEz4jsS5Vv1J6+87XPmu+ztOTQYtl4XFXfICKfAX4DEBq+29CePqbUTskDQTsDf0HzSrfraF5GsTfwQeA1qnrHQN93QMn+wKr6IeBDq61rHMiFFdLjaCB653PEM5RcvHM5m6Xqr/tu89h2lLY9R3LW91Kl77XLqt00bRS/rW9eG2a5z2YVylLGyDue29CbMYOSNwS9lebNzg9U1ceq6q8DvwLcC/ir1VQ+7/DIKD2XwiNA77tVW32DMC1TOnCjsILnQ/rdtitHRH2EUXI+l6c7HxFaZKf22QwzpTYx8im86m3NISJfbP9uEpH/Tj6bROS/h9orCa08DfhV1e2XT1X/u310/2rg5UMrnXd0S7cuPGj7ysvcNNibStvBWkp8qY3SKXjqh0c0kQr2/C9pZxRysbatn7kfBq9+T2l7NnPl++pZpD6bxSWHFsvyQNC4d5UtUeSakniSuIXlC2m5sAPWm7KnxxFSIvTKpR+vfq8Obyre2Y787UNkv7OZqmWriPtCCbaNVu3adpf8SHjqfpn7bFbRhVZKPquBiOwuIheIyDXt392CfMe1ea4RkeOS9ItF5Dvt+zov7557EZG7iMjHRORaEblERPYt8EVE5Pki8tr2+AHtczKDUELk3xaRYx0Hnk+jyJcWh64/m0PXn92rJKOB6akmL0+firVqzNbbpUchhohkciowbbMlUs+2JTUvpOAh8s3WkSrd6EehNCxk/U7T5rHPuv+n84AtWvZZJbotQfYDLmLlHlJAQ/bA62hW6x0AvM4Q/jGq+uj20y2XPgG4WVUfDLwdKHmJ8nuAxwO/1x7/FHj30AaVhFb+EPiEiPw+zQJ1BR4H/ALwrKEVLiq8ARhNx3PT+0hl5sIK3vkc8Ub+50IhEfFEZT1ijUIKtoxFRJq2nlwdJbYWuc/mBVOKf/duCQI8Bbig2xhQRC4ADgc+2mP39e33c4B3iYh4EY0EB6rqY0TkMgBVvbldaj0IJZtmXQ8cKCJPBB4OCPAZVb1oaGUVO8IjvSgGGg1Uj5xyZGNt5eKyEUn1kYQXZ/bOe+0tafuQa+O1czWYtz6bFyiDYuR7iEi60uOU9kHAEpRsCXJ/4AfJ8cY2rcMHRWQL8HHgjS1ZbyujqptF5FaaJzdvzPhyh4jsRBumFpE9GeFWQfE6clX9HPC5oRUsAw5df/aKPVi8QdQ32HNEBnn1FtmO6vFUYU4tRzbS9JzPaV19Kr2k7TZvSQjIy5OLi9v29qXPap/Nww3OFFvLJfmNqrouOjmGLUG8X8HOuWNU9XoR2ZWGyF9AsxdVrkyEv6V50PG+IvIm4AjgtYU+bsOQB4IqCuANvmjKXaqYInXbN5WP8nh2c6q3j/Bsm3LhAmsval+frTRfH8F5Nrw2LWqfzRPGFVkZw5YgG9kefoHm2ZmLW9vXt383ichHaGLoZ7RlHgBsbJ+/uSeQfWeDqn5YRC5l+3M0z1TVq/pbuBKVyMcEb0liB0/Rpd+7QWfJxUMu3U6rLcmlSOuMpuqRbetHREY27BARZK59Xp2eOvWuRTSDSPNEmPc+mzclDqCqbJlOkLx3SxDgfOAvkhucTwZe3RL0vVT1RhG5M83y7AuN3a/QKOvPRfFxEXkn/u/WQ0QEVX3ZkAZVIp8AcsSTHufIwqZ7NmwdfTY85MrlSDgi9SjNs1PS7gjeD4VX1gtleOS7aH02j1Cm9mRn75YgqnqTiLyBZktZgJPatLsD57ckvhMNib+/zfMB4EwRuZZGiT8v40Ma3/9zmhUyI0PyN1STjCJC8+z/g1T1JBHZB7ifqn51NQ5MGuvWrdMNG1b19OtISJV5bpqem1ZH03wPOWXYV19EbrmyJSGByM/S6xC1vS+2HPkahWYiG+n5eemztVpm2L4IIYxZl+BeO91HD7n7c4vy/sOmd626vlmBiFzWPjE/MkrWkXfo1jt2r3zbxAjrHZcRuam/l9dTid5UOv3e5Ylsen7YtJziTMMA6XHqS1R/pBA9v9I6cjYj9IVaojyRb5ENe26W+myeoapFnwXDqhs0aBvbcax3XCZE8V17bJVclHeIOvTqtOneOa/etFwa042I2KZFYRavHd4PhEWOVCM/vGtR4otna5b7bJ6hLM8j+uPGECIfy3rHZUE3xfVufkI8UNPBnhu8uRh1NL23N8ts3VEc2CO/NN2r3/M5amuUp09Fe+c9mznC9doXYZb7bB5vbnpYQLXtQkQ2sV2J3y3ZKEsAVdV7DLE3hMi79Y73Wc16x2WFHegRgUXkFZXLkZ0l1dxUP60j8tnDKOf6iDNHdpH96HtOPdv8UWx9nvpsnqHA5iUh8nFtltVhyANB6XpHYcT1jssG+7AQxIQBKwdyjhy6dGvDG+glijMKp3h+R4RfYiPyqwSlvnvXK6rP+h1dx1nus0VR4wBa4GWiHgAAFrpJREFU9+EbCcU3O0XkdOBHqvpuVX0X8CMROW1yri0OrKrrIzpLll7+iFRzaV1827Pb51NOyXp1p6rSfu8jOK/uHLzr5RFxWmefml+EPptHbC38VKzEkNDK/qp6S3fQ3uxc1ZKZZUEXL7/gwKNCVZempfDIxFN1nr3cjbf02NaV2oqUaE6herYj+6V+e3H6nIK29eSUb8n5qJ2pv7adk+6zednRsBSKsrUq8pEwhMjvJCK7qerNAO02j/WBogHopsAdoVui8tQirFS3HXJl+qb9kXq15OKR9hCC9trglfEIs89eZDP1ORfGiGYD6TXwys9Cny1SKGUFdNBeKxUJhhDx24Avi8g57fGRwJvG79Liw1OYOfLzVFkUOkiPcwTplfPsR0o3stmnkq3tSHnm6rGKNtd2aycl+VSJe222hDorfbbIqDHy0VAcI1fVM4DnAP9Js8nMs1X1zEk5tsjoXkYRqdvcja6UkKyCtp++MICn+iyxRmEHa8eSc2rb89XaifKn7Y7Ub2TbkqlH+h48G9YHm39afbZo4ZQUCmxma9GnYiWGPNkJcAPwVeAbNPsBHzJ+l5YDOaXWoSRkkbPhqVyPSKxPXv1WAfeRl+eXrS8i4NSfEoXu2bbpljg9mxGBp/547Yr8m0SfLTa0+F/FShSHVkTkRTQvWt4buBw4iGaXrydOxrXFht0tsTQO7cVc0+/elN5TeEOm8J6tKKZrz9nypTHvIaGF3A9Vmt53jSO/ozzT7LNFVuIdFOrNzhExJEb+cppXvK1X1SeIyENpdu2qWAXSJ0D7bhR6YY5IrdmYb27anov35sIitj6vbOSjR5q5tqX+565H5I9N8wi5j4yHziYiH4f02TIQ+DYIbJUaNhkFQ4j8dlW9XUQQkbuo6tUi8pCJebZkyE2bvVBA383ItFwUSvDqjkjGi617xOrF3W1a5GtE7lHe3A+TPc751edT5Ps0+2xZUBX5aBgSI98oIvcCPgVcICLnAj+cjFvLhy7UUjLYIY7r2vId4UV207hxX4zY8yEKDXh1pP7kyLUvpOH5EcW/bb6+dFu/9dvLM6k+Wyo1ThMh31L4r2Ilhjyi/6z26+tF5PM0rzH67ES8WlJE6rY7hrzijEIokcpO07wwRu6850vqdy4eHv1gRATf53PkSy6GXWrDa5P1fZJ9tmyooZXRMNIDPar6hXE7UuHvmBjFcb24rkcWETF6oRNP0UZk5fliv/cRWV/dtl4PJfHtKH/0Y5WDF8fOXbNR+mxhH/jpQfNkZyXyUTBk1cpdgZcAB9PcYP4i8F5VvX1Cvi0t0kf6O+SUbu6GmyWOiFy8kEtk3/rgKX5PmdrypXVZWDVr2xyVi3zt+5FI64x+pIbOTrz8y0rgKSqRj4YhMfIzgIcD7wTeBfwaUB8ImiD6QiIemXnkE+VN83i207r74r/WZxuz9hS9V59XVxSb9vJ4Zbyy0Q9fCcFHsw87M4nsl/4QLx+63Vb6P6uBiOwuIheIyDXt392CfMe1ea4RkePatF1F5PLkc6OI/E177ngR+XFy7kWrcnQAhoRWHqKqj0qOPy8i3xi3QxXbkd7suvCgI0OySEkyFz6JSDQ9133vzkVqu0SRWkSqvySkEsWVSxR86azDXofctegLDfURfPe3qvDtUKYWI38VcJGqniwir2qP/yTN0O4l9TpgXevapSJyXrvX1KOTfJcCn0iKfkxVXzrpBlgMUeSXichB3YGIHAh8afwuVeTgKWqrfFNiieK3XV4vfhuRWmrDEmpnw4MlxujHxhKirb8kxhzNQCIy9fKlfuZ+6KLZjHe97LHXNxXKFu4o+qwSzwBOb7+fDjzTyfMU4AJVvakl7wuAw9MMIrIfcB/gX1fr0Gox6J2dwLEi8v32eB/gKhG5AlBV3X/s3lVsg/fgUG7qH8W1u+8dPDKOwgNWfZeEIbz6LSKC7lP5OSXuxbOjcl7Yx7tWkYKP2pCm2/PLtrSwBFO82XlfVb0BQFVvEJH7OHnuD/wgOd7YpqU4mkaBp53+nHbrkn8D/khVf8AUMITID+/PUjFp5MgtIl6P0KIYbvq3RJHasp6y7ZCLLXs+RWGJ3A9RdD0sSVvFbW1F57zvtk7P76q+yzCAyPcQkQ3J8Smqekp3ICIXAvdzyr2m0L6nOmwnPg94QXL8D8BHVfVnIvJiGrU/lS1MhhD5AcBnVXWTiPwp8Bjgjar69cm4VuEh3aMlUns5WBLz1Gqk0kvi4hGhR3F8W6eFp577fjSiH4TUd6+cPd/3I2TTSwi7KvEcmkeCCnGjqq4LLakeGp0Tkf8Ukb1aNb4XzW6uFhuB306O9wYuTmw8CthZVS9N6vxJkv/9wJv7GjEuDCHy16rq2SJyME386K+A99KEXFaF9sbCx4B9gf8AjupeYOHkvQdwFfDJtbipMCtICaFbpugp0JK4bno+UpQdPDWdyx8RvT0Xqdic+vXsRLHwnCLPtcm7NrYuLzST5qs3NMugTG354XnAccDJ7d9znTznA3+RrGh5MvDq5PzRwEfTAt2PQ3v4dBqemgqG3Ozsfip/l2b9+LnALmPyo7uLvB9wUXsc4Q1AfSApgado+1Rq7lxKULlYe2rDIv0RKc2f8yc36/B+wLz6cj9i0fVLv3vlPfUfnavoxzSWH9IQ+GEicg1wWHuMiKwTkVMBVPUmGq75Wvs5qU3rcBSGyIGXiciV7Wq+lwHHr9bRUogWvlpJRP4RuB44FHgs8H+Br5oliaM5IfId4LeTqc7FqrrDhlwi8ljglTRbA6wrUeTr1q3TDRs29GVbCHSvkOuL50Zk2Fe270aeTbd15vJ5iMqW/mhFeXJKusS3SKlbH5dJiYvIpblQRwnuvNOuusfdHt2fEfjRT7+46voWCUMU+VE0043D25cw705DquPAirvINEt6VkBE7kTzurneOkXkRBHZICIbfvzjH4/JxdmHd4POfs/FkdO8uViwTYtUcEk83iKKy/fNCqIwkBeb74vvR7C2vOub1lkxFMpW3VL0qViJIZtm3Uay8L0l3BviEisxhrvILwE+rao/EOkbcHoKcAo0irzUx3mHjZuXhESiG6ZRuCVSuZ69qJwt6yno0tlBNMsoDeHkYtxRWCkXT18mFT4JjCFsspQYadOsUTCGu8iPB35TRF4C/CKwi4j8VFVz8fSlhSWUdCMu8NWyR1Ldd1vO2rL5onh0asOLbfeFK2w9tmyfyvaI2/PLI3NPgVfiHicGrVqpSDA1Iu9B711kVT2m+y4ix9PEyCuJF8LeePSI2qblwinpd0uAferd+pUjzFwM3drw0nKhpJwij2LkJf5UjAYFtmpV5KNgyO6HLwU+HC0LXCVOBs4SkROA7wNHtnWuA16sqlPbfGZRkSrHTp17IZAuvYOnnlPkVHVa3vsB8H5McqtCvNBKRKyRai89l/OxrgWfEFTZqqt+/H4pMUSR3w/4moh8HTgNOF9Ll7z0oF1I/yQnfQOwA4mr6t8BfzeOupcZHomn57xQhYVH8jlFHSnw3I+A52vJypHIjuezp8b7VH3FeDHFR/QXDkNudv6piLyWZmH8C4F3ichZwAdU9buTcrBi/Mg9TNR3M9MiR26egk6P079erHpIiKMkLOLlt+2IZhhVhU8HWkMrI2HI8kNaBf6j9rMZ2A04R0TeMgHfKqYEL56dC13YZX2evTSvtRup3iiE4/0QpD7Y77betC1eiMdeh1woqWKSaG52lnwqVmJIjPxlNDcibwROBV6pqne067uvAf6/ybhYMUl4qy7sChfof3ing40te6q8JJRi7fWFUnKrS3JlUpsdqvpeO1RFPhqGxMj3AJ6tqteliaq6VUSeNl63KtYSEZF1a9OtYs2FMzwijRTvUGLuyqR/vfSoTF06OGvQuo58RAwh8r8Ani0ix6TlVPUkVZ3a5jAVs4Ho5mAuH/gxdS+UEy0LtPmjuiN1X0MmswtF2bq1rloZBUOI/FPArcClwM8m407FLKM0DBOpY++cTStR3H3nvZuYNVwyH6iKfDQMIfK9VbW+XKJiB0QhFZvHKuz0uxcKKVlOaOuqKnyOoTVGPiqGEPmXReSRqnrFxLypmDv0Kd3uBRjeAze52HeXZr9HIZoa714E1Bj5qOgl8u6dnG3eF4rI92hCK0KzIrG+q7MiRA1pVJRCAa07G46EEkVeV6RUVFRMAQpVkY+EXiLvlhuKyJnAvwD/qqpXT9qxioqKZYOyVTevtRNziSFPdn4Q2At4p4h8V0Q+LiIvn5BfFRUVS4mthZ+KFEP2WvmciHwBeBzwBODFwMOBd0zIt4qKimVDXbUyEooVuYhcBHwJeC7wHeBxqvrQSTlWUVGxbNCpvHxZRHYXkQtE5Jr2725Bvs+KyC3t+4rT9AeKyCVt+Y+JyC5t+l3a42vb8/uuytEBGBJa+Sbwc+ARwP7AI0TkFybiVUVFxZJiKqGVVwEXqep+wEXtsYe3Ai9w0t8MvL0tfzNwQpt+AnCzqj4YeHubbyooJnJV/SNVPQR4FvATmpj5LZNyrKKiYtmgqG4t+qwSzwBOb7+fDjzT9Ub1ImBTmibNC4OfCJzjlE/tngM8SfpeMDwmDH1D0G8CjwWuo3m5xL9OyK+x4dJLL71RRK7rzzkYe9DsBDkvmCd/58lXmC9/J+XrL4/BxvmweY/CvHcVkQ3J8SntS9dLcN/25fG07wm+zwAf7w3corptec1G4P7t9/sDP2jtbhaRW9v8E/+/MeTJzl8A/hq4NGnEzENV95yEXRHZoKrrJmF7Epgnf+fJV5gvf2fZ13FuASIiF9K81cziNas17aRpwbmJYsiqlbdO0pGKioqKcUFVD43Oich/isherRrfC/ivAaZvBO4lIju3gnZv4IftuY3AA4CNIrIzcE/gptFaMAyD3hBUUVFRsQA4j+YlObR/zy0t2L4l7fPAEU751O4RwOfG9V7jPlQiHx2l8bhZwTz5O0++wnz5O0++TgonA4eJyDXAYe0xIrJORE7tMonIvwJn09y03CgiT2lP/QnwChG5liYG/oE2/QPAvdv0VxCvhhk7ZEo/GBUVFRUVE0JV5BUVFRVzjkrkFRUVFXOOSuSFKH2st817DxG5XkTeNU0fjQ+9/orIo0XkKyJypYh8U0SeO2UfDxeR77SPNO8QT1zLR54dX/p8fYWIfLu9jheJyDjWVY+MPn+TfEeIiIrITC5JrChDJfJylD7WC/AG4AtT8SpGib+3Aceq6sOBw4G/EZF7TcM5EdkJeDfwVOBhwNEi8jCTbc0eeU5R6OtlwLr2RSvnAG+ZrpfbUegvIrIr8DLgkul6WDFuVCIvR9FjvSLyWOC+wD9Pya8Ivf6q6r+p6jXt9x/SrKedyANUDg4ArlXV76nqz4G/p/E5xZo98mzQ66uqfl5Vb2sP19OsL14rlFxbaATHW4Dbp+lcxfhRibwcKx7rBXZ4rFdE7gS8DXjllH3z0OtvChE5ANgF+O4UfIPkceYW6aPOO+RpH77oHnmeNkp8TXEC8JmJepRHr78i8uvAA1R1xc5+FfOJIY/oLzzG8FjvS4BPq+oPpiEcx/UYcvt025nAcTq915iXPM68Zo88GxT7ISLPB9YBvzVRj/LI+tsKjrcDx0/LoYrJohJ5gjE81vt44DdF5CXALwK7iMhPVXUiDwaM4zFkEbkH8E/An6rq+kn4GaB7nLlD+qizzTP1R54DPzp4viIih9L8iP6Wqv5sSr556PN3V5rtqC9uBcf9gPNE5Omqmm5EVTEnqKGVcvQ+1quqx6jqPqq6L/DHwBmTIvEC9Prbboj/SRo/p/26+68B+7Wb9O8CPI/G5xRr9sizQa+vbajifcDTVXXI3h2TQNZfVb1VVfdQ1X3b/6vrafyuJD6nqERejqLHemcIJf4eBRwCHC8il7efR0/DuTbm/VLgfOAq4CxVvVJEThKRp7fZ1uyR5xF8fSvNLOzs9jraH6WpodDfigVCfUS/oqKiYs5RFXlFRUXFnKMSeUVFRcWcoxJ5RUVFxZyjEnlFRUXFnKMSeUVFRcWcoxJ5RUVFxZyjEnlFRUXFnKMS+ZxDRPYWkXPbfce/KyLvaJ/my5W5V7uNQJr25VX48FNzfO/kAaMftXuzd8eubyKyr4h8q8R+RUXFSlQin2O0W7p+AvhUu+/4r9I8XfimnqL3otngaxtU9f8Zl1+q+hNVfbSqPhr4P8Dbu+N2W9WKVSL3YpMp+rCLiNx9rf2oWDIiF5Hni8hXW2X4PhHZSUQe177V5a4icvf2bTmPiNKNvbuLyD+JyDdE5Fsi8txWWV4tIqe35c8Rkbu1+T8lIpe2tk5M7Bzb5v2GiJyZ89c06YnA7ar6QQBV3QL8EfD7InK3jC8nA7/S2n1rW9dPk3ojf1z/R+iHyM7O3nUzZbPXxOuTNn2Ufploentuh2td0O8AG0TkIyLyxPYH3V6nrr2nttfhwyJyqIh8SZrZ2wFD8hnbvyYibwO+QyMeKtYaqroUH+DXgH8A7twev4fm7TgAbwT+iuatKq9OyrjpyfnnAO9Pju8J7EuzZehvtGmnAX/cft+9/fsLwLdo9tZ+OM2A2MPkCf1N6nsZjdq1fl0G7B/50qZ/y5T5afvX9SfyPy0bXPfXd+3P2em5bj8dcE126JP276B+mVL6Dte6pI1t+k7A02hmZFcB/xv4peT8vsBm4JE0gu3Sts1C85KJTw3Md3fghcAXgS8BLwJ2XetxXT/NZ5kU+ZOAxwJfE5HL2+MHtedOotlYah0rX9EVpXe4AjhURN4sIr+pqre26T9Q1S+13z8EHNx+f5mIfINmt7kHAPvRqOpzVPVGAFXttmnN+dtB8PfFTtMjXyJE/kT+j4LITp+vJdck6pOc/cifSad717qkjajqFlX9R1V9Ns3GZw8Cvm8U9L+r6hXa7DF/Jc2r/7S9RvsOzHcDzQszXqSqv6Gqp6rqJutXxdpgmfYjF+B0VX21c253mtjynYG7Av/Tkw40r0qT5tVuvwP8pYj8M3AGO5KrishvA4cCj1fV20Tk4tZmjowjfztcSaNAtxdq9hd/AM2bfu7j+ZKx19W7Q56M/4PQY6fP195r4vWJqp4U2Y/8mXR60p7BbdyWUeSewHNplPIdNET7zSRLuif61uR4KyvHfkm+I1r7nxSRj7Y+XtfnY8V0sEyK/CLgCBG5D2x7y3z3pvNTgNcCH2blC36jdFobvwTcpqofognBPKY9tY+IPL79fjTNdPSeNC8Svk1EHgoclPh1lIjcu/OrwN+0TXcTkWPbPDvRvGru73T7+yM9XzbRvFwguk6eP5H/Q5Gz4/lqfctek0yfRPYjfyad3rXHXuuSfkdEPgR8nUaJH6uqh6jq6ao6kfdvquo/q+pzaWYxtwLnisiFIrLvJOqrGIalUeSq+m0R+VPgn6V51dUdwB+KyG8Bm1X1Iy0RfllEnkjzVpUd0lX1c4nZRwJvFZGtrb0/aNOvAo4TkfcB1wDvBbYALxaRb9LERde3fl0pIm8CviAiW2ji28dH/gLbVJCqqog8C3iPiLyW5of50zTx0g47+NKSypekWe73GVV9ZWLT9Qf4rOf/CMjZ8a7bNpRcE+I+iey7/ZLxc1zp7rVW1eML2ghwFs3/k81MEar6E+AdwDvaMM6WadZf4aPuRz5mtArlH1X1ET1ZJ45Z8mWtUa9FxSJjmUIrFRUVFQuJqsgrKioq5hxVkVdUVFTMOSqRV1RUVMw5KpFXVFRUzDkqkVdUVFTMOSqRV1RUVMw5KpFXVFRUzDkqkVdUVFTMOf5/Wqo1Wb39CFgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "label_pack = Labels(xy_base='Optical Table', z='Height',\n", " xy_additions=['exx', 'why'], xy_addition_side='left',\n", " addition_joiner=' space ', unit_prefix='>', unit_suffix='<',\n", " unit_joiner=' spaaaaace ')\n", "\n", "z = NollZernike(labels=label_pack)\n", "z.plot2d()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Most of these have sensible defaults, a more reasonable usage would be:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEOCAYAAACUxJyzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO29bdglRXUufN+CYHKCyocaXj4ccyQaNAblCeB7jDEKiu+VgHkDCEEZcvBwqfEy0Ss50ddEc0DPQU00x6g5GZWIxA/ALyZRQwYUT+JxCDNCVDQKfjJgEkaQYBR15lnvj+4e6lmzVlXtvXvv3Xs/655rX093ddWqVV1Td9+1uruaIoJAIBAILC7uN28HAoFAIDAZgsgDgUBgwRFEHggEAguOIPJAIBBYcASRBwKBwIIjiDwQCAQWHEHkgUBg6UDyZJJfInkLyZcZx59M8jMkd5E8TR3bSPLm9rcxST+W5Odam28iyVm0pQZB5IFAYKlAch8AbwHwTABHAziL5NEq2zcBnAvgParsQQBeBeB4AMcBeBXJA9vDfwbgfABHtb+Tp9SEkRFEHggElg3HAbhFRL4qIj8E8D4Ap6YZROTrIvJZAKuq7DMAbBGRO0XkLgBbAJxM8lAADxSRT0vzFuW7ADxr6i2pxL7zdmDaOOSQQ2TDhg3zdiMQWGps3759p4g8ZBIbzzj5cfLtnd+trO9rNwG4N0naJCKb2u3DANyaHNuBRmHXwCp7WPvbYaQPAktP5Bs2bMC2bdvm7UYgsNQg+Y1JbXx753dx3bYLq/Luy+fcKyIrnjtGWu1aJF7ZSWxOHRFaCQQCg4BAsLq6u+pXwA4ARyT7hwO4vdINr+yOdnscm1NHEHkgEBgGRLC6+oOqXwHXAziK5CNI7gfgTACbK724CsDTSR7Y3uR8OoCrRORbAO4heUL7tMo5AK4cr6H9I4g8EAgMAgLBquyq+mXtiOwC8CI0pPxFAJeLyE0kLyB5CgCQ/HmSOwCcDuDPSd7Ulr0TwIVoLgbXA7igTQOAFwB4O4BbAHwFwMf6PgfjYulj5IFAYFEgkAJJV1sS+SiAj6q0Vybb12NtqCTNdzGAi430bQAe24uDPSOIPBAIDAT9Efl6QxB5IBAYBkQgq0Hk4yCIPDBVbDn+jL3SSIEI92yn0OnpvlVGhCAFJ269on/nA7NHKPKxEEQemAhbjj/DJNYO6WoUHRE36Q0x1xC6R+5p2tUnnO7aAhBEvxBYhez+/rydWEgEkQcmQkqupKwh3u54Su5rSX7vC4BOS22kZXS9OZuBxYBIxMjHRRB5oBqW+gbWErClii0yTUm22/aIV+fVaWk93kWj8z2t66TrLh/5HASmCQEiRj4WgsgDY0Or7y4NsMnZCp9Y4RJtU6drhV5zIdD1BQYICSIfF0Hkgb2QqlfAJuwUlrpOyTsl2i7NUtNdeb2vSdrb9kg9zZv6p2/Edvkinj5HRGhlLASRr3NcfcLpADSB7n0zUW+naRbJ5+Liubi2F+P2iL0UXrGIXtel/U/PSZD77EBZBXfdW84Y2AtB5OscuZi2RcBpur7BqbfTOqwboWn9FnJhk1woxWuXNTPQbdP+lnwM9IkIrYyLIPJ1iPSmZfd4YCl+bMWkU3gkaxGktmfta1tp+VwM3boIWReZmgtIaqcLw5DxzPr0IGCEVsZCEHkAgP+MtiZBTyHnnkLJPQfuwSLc3KxA+5se82YMnp+pfW8/MAUIgPIStQEDQeTrBF3cF1irwpt9+5nsnIIu5fNIUZNmLtbd/c09+ZK7aerF0r0Qj/azy6vJPj2Xoc77hIARWhkLg1rGtvTl6yTfaSSFpPeFkICCpYItgtPHvbypTSufJkXvQpCLm6czAK8N+sLhhWhSH6x4u9V2bxYSMfNpQRpFXvMLrMFgFHny5euT0HyN43qSm0XkCyrfAQBeDOC62Xu5OLAeIdTIhS1yitpKs5R77mapRaqlG6bWtgUvdGK1PXcsN5Po0tLYORAKfSKIgLuKH40IGBgMkSP58jUAkOy+fP0Fle9CAK8D8DuzdW/4sMInzbZP4l5YxSPq7liOJDUsorRsWPXnYvRpPv3XCsdYFydvpqL3rYuQdaHp1nyJm6JjQCTU9pgYUmjF+3r1HpB8PIAjROSvc4ZInk9yG8ltd9xxR/+eDhwd4XlhidJ2Wl4TnUecVh7Ldu4mpVW3d5FI7Yyi1tOLQlrG+1m+eXVH2GVycHV31S+wFkNS5NmvVJO8H4A3Aji3ZEhENgHYBAArKytL/7iBXkcEKKtKoPwMdwovBFGKq6f5cqGcXDglFx5KbWmSzsXXvZupNW2zlL11cehmSKHMaxGKfFwMSZGXvnx9AJrPLF1L8usATgCwOW545m8cenl1OMMKnVg3KcdRnB4h1t5szRG457+VViL1nP9eOEeTesnngA+K9KbISw9OkNyf5GXt8etIbmjTzyZ5Y/JbJXlMe+za1mZ37KE9n4KxMSRFvufL1wBuQ/Pl61/vDorI3QAO6fZJXgvgd9rv6K07pC/11DxVodWiFffNxcHTdKtOXYd18zPnk6XEtY3aG5OWn7ljlo3cfYDcebb8TtdKjxUXMxABd/1wYjOVD06cB+AuEXkkyTMBvBbAs0Xk3QDe3dr5WQBXisiNSbmzh8g5gyFyEdlFsvvy9T4ALu6+fA1gm4hsnq+Hw4D3xR0P3s08jwC92HVNrLpUt0eOmqh1eKSmbbkwilXGU9me3Vx4KHcDVsfRI9xSQD+hlZoHJ04F8Ift9vsBvJkkRST9j34WgPf24dC0MRgiB4DSl69V+lNm4dMQYanDNN1So5qQNMl4hJ5T1hZxWX5qYrWIW5f3LhxWe63ZhAeLXC2iLl0UrIuOdw692H9AQ8DV1T4MWQ9OHO/laUXk3QAOBrAzyfNsNISf4i9I7gbwAQCvVsQ/NwyKyAM+OiVHxQFeyELHn7t0fbwUHrCgQzLpdo6AS/56xKfrSsvUwJsZWH5b9xC82UGu7aUbz6HMDQhGUeSHkExDHJvahxyAwoMTNXlIHg/geyLy+eT42SJyW/suywcAPBfAu2odniaCyBcINcovh1w5TzVadks3Da0bgLk4t+eTp3itPF55q+1e27ywTy4EU5od1ISwAh1Gemplp4h4DzqUHpxI8+wguS+ABwG4Mzl+JlRYRURua//eQ/I9aEI4QeSBPNIXfACfjCwySZFTuqWbnKXYsxfvHoU4dfssItXwYtaen2kZr20eSdfY98JMpT5L1z6PG6EApZfQSvbBiRabAWwE8GkApwH4eBcmaR91Ph3Ak/f41ZD9g0VkJ8n7A/hlAFf34WwfCCIfMHIKUeextjU8xWyVswjXCztYqt2qQx/zQj+e/zV+6/zeRUqfD+9YWlbb1LF2z1Ztn617iAC7ftSDmaoHJ94B4FKSt6BR4mcmJp4MYEd3s7TF/gCuakl8HzQk/raJne0JQeQDRPqadwodnsiRQO54LWmldeZQe7OvJuzhhTRyFwbdptKMwgud1ISASqGmSfqs6/d1q8xFgH5udqL04ISI3ItGdVtlr0Xznkqa9u8Aju3FuSkgiHxAsN7QBOzHAHNTf49MSkrRypMjJl2vV48XAtFlatqWquCckvXI3PLXqsfKm7PRZ5+tZ0KP1+/HQxD5wJBTi11amrdDSfGmhFZSn54q1kjtpvu6LR7B6zKeHa+t2ndPjVuEap0Xy+/c7KRG4etyo/TZ+kN/iny9YUiv6K9r3Pd4oT2dT7dT4tL7VjlPhVvlPOWs6/OUvM7b1ef56PnhtV23IYVumzUjyM1OPFtW23OzFsvvcfpM3+xeeggaIq/5BdYgFPmAYClGve+pvzR/bvruKfMaUs/Vl4ule/FsL82bCVjKNnexqvHfCtPk1Ll3rqfVZ+sLocjHRRD5nKHXTOngEYs3wK3peQ5eqKIUfvDCBJb9VHnm8tZcVHL5tN8Wievy+pjlo5Uvdw6m0WfrKV5OEbCHp1bWI4LI5wT9ZZlu2wuDWPtd2iiwiMgjuhKhWmVSmzWxZq/uEuFaxG2FLGpmClYZyy9Pvc+iz7r/L0tP6KHIx0IQ+ZyQI1SLMNJ974Zarfq19kvEa5G6R3IWaVuEpn33yLKkvnVdaT5P9Y6Slh4bQp8tLboYeWBkxM3OOSBdwTANWwB7E2GKGgWpUYrBdvVbedJjFoF5xJ/zSZO2RdJeaMTa1hcQi9hHIf60TVrl6/Myjz6zVr9cHkjc7BwTocjnAGugWsRlhS+stBIhdLZyYYxSee1P7m+uTZbqLoVmLHtd+ZqZg/bBUte5er02zrPPlhICYHWJ2zdFBJHPEPe9sVmeiqf7JVWsy5UUaEoknrrWJKXJSMMjIO2z1cYUtWEd61zkLijW+fHOk0e+Q+mz7gb58q2cKMCuXfN2YiERRD4D3PeFmL2fZbb2u21LgXbwYslW3pLy9pSeR1A1Sj5HVpa/ufRaW7Wq1TtPli1NyEPqs6VbCjcU+dgIIp8RclNvS8F5CtMLS5TCEFYZjVysOxe3rWlzKUxTqj89P54P+iKTC+fkFLql2FM7Q+uzpUI/qx+uO8TNzikj/V4jUL4xpsMe1l+LENLjVj2duvRCEiUl65XVqrUUNrBCNp7azYV59C89puv32umFP0rnUp+XeffZ8rwBKo0ir/kF1iAU+ZRhKaZRCFQjpyRT26VQg6VerWm+F84oKXRLbdb4rm3lFHiuTda+51/uPGhbQ+6zhUeEVsZGEPmU4C1Fq2GpSCuckQs7dNupjVLZ1Ia2pY+VQhclkiuVS8+Fh1w5Xb42dKTJsJakh9hnSxMvDyIfC0HkU0JuIFoDNUc8o5KLdSxns1b9dds6j25HbdtzJKd9r1X6Vru02k3TxvFb+2a1YV59tsgQAWRXEPk4iBj5lOCp227bI8lcfNib8nsx5lwZiyAtH7xYb7ftkUeJVGqO5/J0x0uEpu1YF5D0WM7OIvTZQkMArFb+AmsQirxnpItgAf5UWg+8WuJLbdROwVM/LKLxVLDlf64uTz1boQjLz1LM11LEXqzYqq82lLPIfbbwi2wFSY+FIPIpQJOaVoQe9CCvLZcjKGsqnlONpXCALu+100srhRI08eZCFaVwSa4+XXbZ+2xhsERNmSUitNITthx/hrsOhjcwLaVr5SmpWE1+ut4u3QsxeCGFnArUJKjbacWp9U/XYcHzTdeREqpVp7ZVUv+L3mdXn3D64q3LIoCssupXAsmTSX6J5C0kX2Yc35/kZe3x60huaNM3kPw+yRvb3/9KyhxL8nNtmTeRHEw8K4i8J3gEpQejpVqtQV1TLnfcyl/y31K2OaVcKpuLA3t1d7a9i4h3rrx2a1jEuB76bGHQQ4yc5D4A3gLgmQCOBnAWyaNVtvMA3CUijwTwRgCvTY59RUSOaX/PT9L/DMD5AI5qfyeP0cKpIIi8B0zyQoZFelpZpumWerOIrhSG0GGLEvnm1KHXrtRGLoZtqX2tPi0b+txY7dft7AOL1GcLpcoFkF33q/oVcByAW0TkqyLyQwDvA3CqynMqgEva7fcDeFpOYZM8FMADReTTIiIA3gXgWeM0cxoIIu8JqZK0CM4iLC/dm5Z7IYkakkrrSW1pQrH882ykZFIbSrBIZ9S267wlFezlWS99tjggsFr5Aw4huS35nZ8YOgzArcn+jjYNVh4R2QXgbgAHt8ceQfIGkp8k+QtJ/h0Fm3ND3OzsAdYA1X+7Y13+URSiVre6Xo8UvTyWXcvPXPs8H7UPOZ9y7SvZ8i4SJZK1zv2y99lCod7vnSKy4hyzjOgrmpfnWwCOFJFvkzwWwIdJPqbS5twQRD4BukcN02VpLdSoXKDuSYRcup56e2ow3bcIyiJCrVxryForSI8gc+2z6tQkmVPMHgGXSHmZ+mxhHklsb3b2gB0Ajkj2Dwdwu5NnB8l9ATwIwJ1t2OQHACAi20l+BcBPt/kPL9icG4pETvKgCjurIvKdHvxZKFjKzvrrKcscWVj1aBu6jpINC7lyORL2SN1Ls+zUtNuDF2rw/PBIfT312UJgtZdo7/UAjiL5CAC3ATgTwK+rPJsBbATwaQCnAfi4iAjJh6Ah9N0kfwrNTc2visidJO8heQKA6wCcA+BP+3C2D9Qo8tvbX+5SuQ+AI3vxaAGQfiCidrpdUq1pGuATT4fcILdseoSn6/Vs58pavnv+Wmo6d248m7mLRym8s177bPAfpBBCdk9O5CKyi+SLAFyFhpsuFpGbSF4AYJuIbAbwDgCXkrwFwJ1oyB4AngzgApK7AOwG8HwRubM99gIA7wTwYwA+1v4GgRoi/6KIPD6XgeQNPfmzEPDIyIOlAr0BWktQpfCCZbO2Pq+sjkt7IQGNUlstUsz5lEMp1JLmWa99Nmj0o8ghIh8F8FGV9spk+14Aez1uJiIfAPABx+Y2AI/txcGeUUPkT+wpz9JAT427NCsuaw1sve+V0XlrlWzOzzTdOmbVq9uoiTsXPrC2S+fJij9bfuSIMHfuRvHFsrUMfTZESH8x8nWH4uWvvXKB5ArJD5H8DMnPtm84fTbNMykq3sZ6KckvtPVfQ/LhfdQ7CrrncnUc0iKfEtF4sVqtotL8aT0WkWo7uVBAWsarV9dtpWtbVjnLRk4hpnmsbd0OHd5Ioc/Zeu+z4T5bPtLjh4EEozy18m4Avwvgc5jC0jbJ21gnoblDfD3JzSLyhSTbDQBWROR7JF8A4HUAnt23L3k/bSWWHgfKN8s8AvOm31653HRdq1tNcqUwRo2KG+dYiTg9ssv5623n1LPOH302fww+9DNQjELkd7Q3CaaFPW9jAQDJ7m2sPUQuIp9I8m8F8Jwp+uPCU5JaIVkDEcjfjNM2rONpOW3DGug1itObmlt+e+RRY8Pzqwa1vlvnK/psuOS9B4LeYuTrDaMQ+atIvh3ANWifswQAEflgT75Yb2Mdn8l/HuZ019iKg1pTe28glwaVpxi1Dzk1a6XV+m0hp2StfYvcxz0vpXNm+WERr647+mxoIFZ7eGplPWIUIv8NAI8GcH/cF1oRAH0RufU/0vxfR/I5AFYA/KJz/Hw0i9vgyCP7eSqy9MhhSh6l4x2sgeZNnfV03LJRQ15Wuq4rteX5XGpLibBG9VsrVCtkYqG2T9Zbn3Vx8kG9JBSKfGyMQuQ/JyI/OzVP6t7GAskTAbwCwC+KyA/0cQAQkU0ANgHAyspKLxIkNyX2ps7pILTKa6Ky1GBq0yMzy4fctF8TrhcWsAhgFIK22mCVscinZM+zmfocfVbXZ0NCPLUyHka5/G01loLsE3vexiK5H5oH9NfE5Ek+HsCfAzhFRP51ir5kUVJLOYXVbXdEkYtzdsdTAknLW3lTAtGK04ImK014ntL17JXK1hBoqR7dJutceMo1+my4EKw9J7lfYC1GUeRPArCR5NfQxMgJQETkcX04Uvk21usB/ASAK9oVJ78pIqf0Ub8HvZ5KaSqv4anAHHHl1FRJzaXwwgAWmem8Oi0XYrDseMq3NO3PXSC8NqUkq+uOPsv32aDWYRFGaGVMjELkU19EveJtrBOn7UMOuQGi81gEkCIlLm9Q1wz8Ll9OhVokUyIMrw0ekXj+lwjOIvfUH88Pq65cevRZXZ/NGxFaGQ/VRC4i35imI0NGSfnUKCtLDXrqKUcQmoC8QZqme4M/RzQ6v7al/dXTfavtpTotcq2BF6KJPqvvs0FACNm9z7y9WEhUEznJFTQ3GR/elus1tDI0dCEVYDTllQsR5EjGSk/td/n0ACxN91O7OUVXUqA5//T50HktH3MhFqttqf+58+H5o9Oiz9aej6GEWEKRj4fBvNk5RJTIZhR1Y5XTdeXKerZqynnTd6tuj2Q0CXd/LYWp6861W7fJI1srb02oIfps7XZNn80Lgvy5CfgY0pudg0LpP3WNyvHIwxpEaZmaKbRHlJ561fZr2lFLdFa7LF+7+i2/vTTPb89+yd9SevTZHCGhyMfFkN7sHDRy8VBrMJTUZGm6XQoz6HpShZoLA1jtsuxr9WjVb7XRIibrHFjbnnou+ez5En02ep/NF4RIPLUyDob0Zueg4U2tSckOopJNPeX11KwX+yyRokeMnurL1Z/6bPmit0tEVqpb12shVzb6bPQ+mzf6+LDEesSQ3uwcDK4+oVlvvjTALdVUIpwuT47wRlW63kC0Bro3oDVJ5Ozn/E/LeOfEU65eXRpazeo2e+U8X6PP7iuz5fgz5nfDM0IrY2NIb3YOBpby0X9LRFCTJ03PkYY3vbbIzLLj5U3zeL7rvyWkdWvysdqm67Pq8s67lccqY5WNPtu7Pd7MZlaQNrRS8wusxShn5EkAbmw//LDmwxLLCk1GQD6Omubp0tIBadnwCMhSpdYUuStTMxX3CED7mk7vrcFt+V8zPfeI0qsjd8688xd91m+fzRqyyqpfCRUfqdmf5GXt8etIbmjTTyK5veW37SSfmpS5trV5Y/t7aI9NnwiDerNzSLDUjDWlt6a5ab5uP81rDThvwFpqLLdfIgeL4Lw2lUIGetsjhpzKzRG0Piel8xF91l+fzQVSP4vIofIjNecBuEtEHknyTACvRfORmp0AfkVEbif5WDRLhhyWlDu7/XbnoBBvdhrIkU1uUFqDRadZ5XLkYfnhHdd+5HyzbOq2eiTlEWEOFulpu5a/Obve+dTtiD4br8/mgT6IHBUfqWn3/7Ddfj+AN5OkiKQfkr8JwANI7u+ttDouSB5UkW1VRL5TY69I5CQ/IyJPmDTPIsH6z2SpNz1ALfVYUpajDEKdv+S3pSBr2mgpN4tUSsSm81ht0ISUwmq7Va93kYg+G6/P5gUBsVr/iv4hJFNlvKldvhqo+0jNnjztgn13AzgYjSLv8GsAblAk/hckdwP4AIBXi8i4V8Pb21/upO8DoOqDCjWK/GcKsXACeFBNZUNH+lFa7z93SQlZJKPTvQGeptdMeWvVnlWmNAW32qn9tsp75JCm5+rUsBRq6aIRfTZZn83toxOCUZ5a2SkiK84xy4g+0dk8JB+DJtzy9OT42SJyG8kD0BD5cwG8q9ZhhS+KyONzGUjekDueoobIH12RZ3dthUOHN6j1APHUnaU2LWWkB5SVrybkoFWjzmcpMkt55Qa8p0yt82D5kcI75qnYnPq17ESf9d9ns0RPs4Kaj9R0eXaQ3BeNGL0TAEgeDuBDAM4Rka/c55vc1v69h+R70IRwxiXyJ/aUB0AFka+n2Djg/6dP9/V2Ckv1WIOqlK+z5R3zSMabantKTrfZspUjJW/QW+SW7lv5PcIuKe7os+n02TzQU717PlID4DY0H6n5dZVnM4CNAD4N4DQAHxcRIflgAB8B8HIR+VSXuSX7B4vITpL3B/DLAK4e10ERubePPB1GeWpl6eENoBpiyKlFb8qdU32eOtV5rTJW2TS/RQ7eQLfs5Hy1yunzmLNplbWIRfdH9Nl0+mymkLpHC4tm6j5S8w4Al5K8BY0SP7Mt/iIAjwTwByT/oE17OoB/B3BVS+L7oCHxt03qa1+rygaRK+RUWXfcUzvedNZCjtRKRJObEnvTdo98cu2uJQePVDxysurUbcv5lyOx6LP++mzWEACrq/2sR17xkZp7AZxulHs1gFc7Zo/txbm16GVV2ZqnVj4K4IUi8vVxKxk6Sq/k56aiOQWWm5JbAyk3Pfbq9pSqV5+XzyIpqz6PPLzzYJXTZWvCAeMo1tx5s3yNPlvbZ924OHHrFW6b+sbqvGYD80Mvq8rWKPJ3AvhbkpcAeJ2I/GjSSocMa4B0sAagN/23BpunyjQJ5KbBaZo3AHPTZ+2fl98jJ8tmjohzdVvKtNR2XY8uG302nT6bCaSf0MqCoZdVZWtudl5O8iMAXglgG8lLkUwBROQNo1Q4dHj/6XPKx1NmpcFZoxrTtLRMzp6VR0/9rUGv06z26oGf+uspRct/3VaPMHMKV9uw0qLPJu+zWUHmUOcA8BvoYVXZ2hj5j9AE+/cHcACW9AtB3iD1FGDtsZxisv7jekor519uCu61x1KGuk1eOCA97vmdlrfIpOZ8lNoVfTabPpsV1iGR97KqbE2M/GQAb0DzuM4TROR7k1Y6NFgD1RuI3pQ5haV0PIVoDUyr3pxyTPPk4CnGGkWbs1lDvtpuKYyi/bLOUWo/+mx6fTZLrEMi30ryaLUOzMioUeSvAHC6iNw0SUVDhjcAS6pK5++gp8S1NnMq0NrOEV2a3xq0HkFYU+scOeQGnqcyLcLU56+rszbEEX22d1qffTYLiBC761/RXxY8CcBGkl9DEyMngP4fPxSRXxjPv8WBpbTSdGugWUTi7etBkTum/bHyjTMQLTIqEYs3mDVZeOpW+1qrwq36vDzRZ9Pvs1liHSryXlaVjefIW3jqq5aIrIGsy3jKyZp+Wza8QZ1rj9c2axruEZ5uRw5aoVoklSNxj1z1+fXak+Yt1Rt9trcNr89mhXVI5P+MZnGuDVjLxxeMYmTdf2qje1YWKCs3T8nkymsi0GGENE+qiDxlaQ1mC94sw1OS1rZHQNb58QagZd/LX6MKvXMXfTa9PkvHyFQhzXPkNb8lwpVoltTdheaBku43Eta9IvcUm5engzXYcoqt2+7sWbYttaj90r56yi2tr0vzVG3ufKTts9rv+Wzl0za9c+Cdw5K692xqX6z96DM7n7Y5TTSfelsqkq7B4SIycXilWpGzwXNIvrLdP5LkcZM6MG946sxSUFotpvlKROMNVssP7Z/O76lWrf70oNa2vDw6v0eUHqy6rGPW8RrFHX02+z6bFdIZUO63RPg/JKf/+GGCt6J5fvypaOI396BZk/fnJ3VintAKp1bReDY8NZcOEk/RjVqHRRbWoLPUYdpWizwtsvDOjVenZ9PattpWcy50XdFn9+332Wezwu7VdRftfRKAc6f+1EqC40XkCd1i5yJyF8n9RqlsqPBUYG5amyOC3ACypqw5xeNN8a0yHjmk/mg/dftzSi7XrtxxK0Sg83pEWDoPVl3euYg+m6zPpg2RdXmz85l9GBmFyH/UftRUAIDkQ7Akb3jmFJGVzyIgb5CWCCfd9ohMq1QPOTLwtrX93F/to5c3d24sUqpR0SVCstqf5os+W1vfOH02fSzdjcwi+vrewyjzmDeh+QsM0pwAACAASURBVGrGw0i+BsCnAPyPPpyYJ/Qg9f7z6zLWgMnZT2150+lceY8w0r+e2k2PWwPcG8zatjf1Lqlcbcc7t2leq63a1+iz2fXZrNDVXfotOkh+po88HaoVuYi8m+R2AE9rk04RkX+qLT9UeFPvFDnVpe3UKrXUtqXK0r+WevIIoEaV5hSwpQ513VYdFglpO15Z7aflo3U+cvaiz/rts1lhGUi6Er1+C7lmrZV7gDUfLmVyTETkgbWVDRE1A1Cnd/n1QLJUkB4c6d+0TJpX2ysNTo9QSkq4RJy6nTqvpYit9lp2SkRW0/7os9n22bQhsq5udvb6LeSaV/QPqDW2iLD+w+sBaJGHhVwZiyi8wW3tp+kWMek2afu6rSXy8kgozaPb7vmkyc3yrYSc/eiz2fTZ9LF+YuR9fwt5UJc/kieT/BLJW0i+zDi+P8nL2uPXkdwweZ33xQRrpt86vdvOKUxLBXmKMN3X/ljk5bUprc9rh1aOqV/aJ207LW+dQ+1PzkbqY/ortSP6bHZ9NgsI0LwUVPErYRIuIfnyNv1LJJ9Ra3OeGOWFoAeQfCnJD5L8AMmXkHxAX460T8S8Bc3jOEcDOIvk0SrbeQDuEpFHAngjgNf2U7fsNcBzA7KDNci7fX1Mq5uUrLyBlNbjqcoaZaf9S9tnkZI+Lx4pWqTrKe3S+bPOgVWn9k37ktYVfdZ/n00b2kfvl8MkXNLmOxPAY9AsaPVWkvtU2pwbRlHk70LTuD8F8GYAPwPg0h59OQ7ALSLyVRH5IYD3oVmDIMWpAC5pt98P4GkkJ5IN3oDMKTedR9vJEZBVh1ZZFsFZg9iqPzf4PJVm7VsDRw9872+6bZGHdz6t86IJJ0ei0WfT67NZoae1VibhklMBvE9EfiAiXwNwS2uvxubIILmXGLXSShiFyB8lIueJyCfa3/kAfnrUCjM4DMCtyf6ONs3MIyK7ANwN4GBtiOT5JLeR3HbHHXdUVW6pIEvdWGShB6zOZw2W1IYmgXTaa5XTfusymrg8BWMN0tyAz6lQTy2mfuj2WXVrcvXIOT0efTb9PpsN6tR4688h3Rhvf+cnhibhEq9sjc1xcJKRNvJLQqO8EHQDyRNEZCsAkDwezbPkfcH636JHQk0eiMgmAJsAYGVlpUpOeNNiiyysMt2+NT21lJhHbB4BefY88vH8tFScR6y6jLZTImXvHKV1a99LytUjyugzu3zffTZNjPjUyk4RWXGOTcIlXrrl2NgnhuQLALwQwE+pxxAPwBi8WvP44efQOHx/AOeQ/GZ76EgAE32eSGEHgCOS/cMB3O7k2UFyXzTPWd45SaU1g90r5ynCkj1rIFsDruSPV14fs+xYxGcpUcuObrd3LjyUpuslIok+m32fzQqrFTcyKzAJl+TKlmyOgvcA+BialyrTG6f3iMjInFajyH95VKNj4noAR5F8BIDb0Nxw+HWVZzOAjQA+DeA0AB8Xkd7+F3pKMN32BoClyLQC6uDVkR4vKS6vXl2HV29JqXlhDCv0UZqiewTk1Z1rp5UefTb9PpsFpL+6xuYSkpsBvIfkGwD8XwCOAvAPaJR6yWY1RORuNOGcs8a1kaLmOfJvdNskD0TTsPRplW/sVWgMiMguki8CcBWAfQBcLCI3kbwAwDYR2QzgHQAuJXkLmqvnmf3UvffUPR303tS7207TPDWUltXwBlvOji4/Sn06r0eyVh0lJZse84il1jerTJqu64k+s+vos8+mi36eI5+ES9p8l6OJNuwC8JsishsALJuT+kpyfxhfCBKRkb4QVB0jJ/k8AL+FZkpxI4AT0FzNnjpKhTmIyEcBfFSlvTLZvhfA1D5X4qmV3EDJkZNFNqndNL81dc8ps5KiTUnKKqftesdSoitN51OfrXOp21RDoDnF6dUTfTa9Pps2+qprEi4RkdcAeE2NzR5wJRplvh3NMrZjYZSbnb+FZu3xrSLySyQfDeC/jVvxUGCRiR4omhisQW2FDDQ8G+l+ST3lpr0eMVjt8+rP1eMpM6vtOULSBOT5Z9Xh+RJ9Nt0+mwUES7Kc6miY7ReCANzbXsVAcv92waxHTerAvFGaCluDtDQ9zaXrgW+VsxSRNdgtvzVZWAM/Nw1P7Vt+WoRkta3bttJ0Pem+d45zZBl9Nrs+myqkeWql5rdEmPkXgnaQfDCADwPYQvIuTHbXdnDwpvo5BVkzFa9RSakdi3xqVGxq31Nnlqqz7Hj1WSRZCpFo1JSz1GuO6LSN6LN++2xWqHn9fhmQPA24L4DfIPlVTPMLQST3FZFdIvKrbdIfkvwEmsd1/mYk7weImsGeO5ZOSb3BNirRaHiDWQ9gTQYWIej9UZSqRSA14QarfHrecnWlflvkHn22d73a71q/Sn02bcg6WjQLPT8NWKPI/wHAE9IEEflkn07MG55KqVV+3tTXG/SaaLyyJSVlkUuNEk73c8qzRqF6qtKro3R+S8rZqzdnM0X02WR9Nm2szq6quaJ7GpDkS43Dd5PcLiI31tqrCTYt/SXSGki10+Jc+fTX5bMGsUUqWrGmeS3f00HnDf5U3VnQCrD7q9PS9nhqOfWjpPS8c6CPW/55/lv1Rp9N1mezQF+rHy4QVgA8H/ctA3A+gKcAeBvJ/1prpEaRP8S5agAAROQNtZUNEVrtlaaweqqb2hmlrpo6PSLSSsvyMafYrHzax5pptkdmqe2csvPakSN0XW/02Wz7bJpoXtFfKpKuwcEAniAi3wUAkq9Cs4jXk9E8kvi6GiM1inwfAD+BZg0A67fQOHHrFe7001JDWq2ledNfCksFafVjDXRLYWnfclNjL6/2zduvUYRWGd1Oqx3aF92OdD/9a52r6LPp9tmJW69w7fSNVbDqt0Q4EsAPk/0fAXi4iHwfIzxXXqPIvzXqW0aLhpQkUtSSRrff5dHwBq0mMo+sPCU7inLy6kpt5trl1WEpvNTHmrq9bcuGReyWP9Fne297bdM+1vg2DQjyF58lxXsAbCV5Zbv/KwDeS/I/YIS1rGqIfOnPbG5qCuytWLo0rZxSWzqfRRoe4ZRCCzlYvubs5vJbRKLtaRLzCMdSoaU2WYrcOlZ7HqLPxuuz2WFdPbUCABCRC0l+DMB/QsO1zxeRbe3hs2vt1BD508bwb6Ewijqx/vPXkJG2k+7rqbBn0yI2i5S0/1Z9Nfn1+dA+WLYtP1Mbus3evraVls+RZvTZ2vOhfbBsW35qG7PC7OcB80dL3NuKGTOoWTRromViFwG5aXitkvEGa2rf2i9Nl2tsidgfOkjtlwjMUoBWWz2FXFNnSb2m8NruKcxc26PPJu+zWUCAdaPISf69iDyJ5D1Ye/0iABGRB45ib5Q3O5ceJTVkTTetgZYb6BaB6DpL6knXl7Op0zy/PL89xVzK5xGZJgpvCu+Rc+68WG2JPpu8z2aJ3euEyEXkSe3fXh4YWapFCyaF/g/tqZJUTVmD0BrQnjrS09gcYWhy0XV0fz3S0LAIzrLp1Wmdny6fJg2PVDxSTcvllHz02fT7bFYQ6e2bnQsDNngOyT9o948gedyodmpe0dfSf88hjDEFGCJy/2G9abhVvjSQLZueD54iqpm6W/m8PKkNPdW37Ofao497KthS5Om2dZ5z5yrX9uizfvpsVliHMfK3oln08akALgTwXQBvQbPSbDVqYuQL/6x4CSdddzkA4OoTmuWJtWq0prnWoPEUlN63Bp838HLTbouEvDo1SlN0b9DX+pCr3yNaq35N7h7ZRZ/132ezfH68wzKp7UocLyJPIHkDAIjIXST3G9XISDFy6wtBIvK/R610qCgNFA0rX07BaCWqbaT7lvqzFJI1OD3Cso5birjUBuuvlcdTfTpkoknGsmXVkU7/o8/677NZQ7Au1yP/Ecl90E5GSD4EY5yGQX0haAjQCtELAegyumwKKy03gD1CSPNY/tbk8+q2zkOpDbl26Xp0m7Q/1rnzCNTyNfqs/z6bPWYflx8A3gTgQwAeRvI1aL4f+gejGhnlZmf3haBviMgvAXg8gDtGrXDIsAaUl1YiiFwd1gC1iC6noCwCKCmzNK8OZ1jTcE/xjgqPONL6S/k8u9Fn9+WZZp/NAoLmqZWa3yQgeRDJLSRvbv8e6OTb2Oa5meTGNu3HSX6E5D+RvInkRUn+c0neQfLG9ve8YptF3g3gvwL472i+7/AsEbl81DaNElq5V0TuJYnuC0EkF/4LQSk85Zebdnqq0xpknurMDWY9GK282rZlKzflt/zy7FntzSnVnJq28ltlLIXt5fX8zB2LPrP9ngdmtIztywBcIyIXkXxZu/97aQaSBwF4FZrVCQXAdpKb0ax/8kci8ok2ln0NyWeKyMfaopeJyItKDpD8U9j3dh9FEiLy4lEaFF8IStDd9Nxy/BlV/6mt0EC3r/N027nB1qmpnPqz/MmpJy+/R4BaNdaENbx6LOWbIxMvPFLTtlwYxSoTfeb3WTcO5oEZzQNORbNULABcAuBaKCIH8AwAW7oXIkluAXCyiLwXwCcAQER+SPIzaMLNoyJ9k/O/oblojI1qIl/WLwRZyE33dT49UK1BXyIYi8C8weqpw7SMpUZ1ee27V19OWVvEpc+P1W6LuHV577xb7bWUqYfos7o+mwe658grcQjJlAw3icimyrIPE5FvNXXKt0g+1MhzGIBbk/0dbdoetML2VwD8zyT510g+GcCXAbxERFIbeyAilyR2fjvdHwdjvdm5bF8IGgWeyky3vcGZU5q5KbNlM63f8jEd0LoOyx+NGsVYmup7ZTx/PeLTdaVlahB9Nl6fzQMjPK6xU0RWvIMkrwbwk8ahV1Tat07ynhNFcl8A7wXwJhH5apv8VwDeKyI/IPl8NGq/5mGQiTug+mYnyUvaK1C3fyDJiyd1YMjwBnA3oKxBYZXNKUdLkeaUbC6E4R336suFI1JbXvut9LRsKWxhwVPGul4vT/RZP302Dwj6u9kpIieKyGON35UA/oXkoQDQ/v1Xw8QOAEck+4djbSh5E4CbReRPkjq/LSLdGuJvA3DsKO2fBKM8tfI4EflOtyMid6F5cmXpYA1uL4+n4ErkpQeQtu1Nv61y1nFtz1NyVnjDIjJPtXlkWJriW+fFOn+eYkzbpM9h9NnkfTYviNT9JsRmABvb7Y0ArjTyXAXg6a1gPRDA09s0kHw1mtDyb6cFuotDi1MAfNFzgOQ9JP+N5L8BeFy33aWP2qBRQiv3I3lgS+Bo7+qum0W3LGWWm5LXTK314LPIoURs3rESEXntsBSdJgaLAL2wQ1qHPpYj9pyKzfmt80efjddn88KMvv5zEYDLSZ4H4JsATgcAkito1gN/nojcSfJCANe3ZS5o0w5HE575JwCfIQkAbxaRtwN4MclTAOwCcCeAcz0H+n5jfhQi/mMAnybZvbd7OppnH5cOJ113+Z4nV7TyswaPRwhWnhSlqa5HdBZZeLZz9ZZIK60zhxwx6/Nh+WVtW+dQ+++1Kfps/D6bx2v5HQSzefxQRL4N4zsL7brgz0v2LwZwscqzA3b8HCLycgAv79XZSozy1Mq72rvEv4SmIb8qIu7UYZngEQOQf7qhRtVZNjwlmbORHvPIpKQUrTw5YvL8t2xY+XWZmrZZ4QkL0Wd726rps3mjh7DJukTN6od6AXQmx5Zi9UMLnrKzBmea35tGW+kpatSiLudNtz1SqVWfnir2zpFHBtY5yxFejlQ8RW5dMKLPxu+z+WLpPqw8M9SsftjrAuiLAk/hWIMqN53XtrS99HhuWq23c+rOy6t9s9pXO323fPQUX67t2g9NKjly04g+66fP5gURYHco8rEwyuOHr61JWxZYU2dLIaX5rSk/uffTG9qmHoQ6PVVQlrpLj+v60/2cH1ZZT7lZ9eVUXXrurPNp+eMRtbZpkWT02eR9Ni+sSt0vsBajPH54kpH2zL4cGRpO3HoFTtx6hataPGLQA1kPIG0jhR6Ell1PNY06MK1BbRFObspeIm/tuwd9Lkvq0lPdnp/RZ+U+6/6/zxtS+QusRU2M/AUAXgjgp0h+Njl0AIBPTcuxIUGrIgt6EFlK0Zvip/V4KtPaL/nktUUPbI/ocoRqhSi0TU8J15CsRbiWz7oOT8laiD4bFiUK1uWHJXpBzVMr7wHwMQD/A80qYR3u6RaUWWZoZekNzNq09Jg3OC3CSPctgrTKW8ip3BzxeirYIjmLtK3zoX33yNIj9dwFIvUl+qzcZ0NBPLUyHmpudt4N4G4AZ03LifblossAbADwdQBndC8eJXmOAfBnAB4IYDeA14jIZdPyqYNeEbFWeeUGTG46nBvYOdVZmr5rEtJTc03Alp8eKeVITrejRH6eKtbHddtyF4/oM7/PhhBO6SCIm53jYpSbnQ8g+VKSHyT5AZIvIfmAcskqdOsDHwXgGqxV/h2+B+AcEXkMgJMB/Em69su0UVJaGt6UPy3bpacDVA9ESxGmaZZvls/edDynxlL/0ro1KVmEm/qn86R2cmRvpevz5Clyq23RZ35fDQWrlb/AWoxys/NdAB4D4E8BvBnAzwC4tCc/TkWzUhjav8/SGUTkyyJyc7t9O5qFbh7SU/3VsKbv6V/ruDWQNQFoe7q8tm2RRFomR4TaX12HRVJpOU2oHhFbxF4iI0s163Z5F4yaC4qVL/psIJB4amVcjPKK/qNE5OeS/U+Q/Mee/KhZH3gPSB4HYD8AX+mp/iK6KeiW489Yk26pTn28pMoslWYNaovgcqSnt73B6ylUj6BK9qzQhGUv54+1X7Kl/fMQfYa5fjzCg7S/wOgYhchvIHmCiGwFAJLHY4SnVnpYH7izcyiamcBGETFnWSTPB3A+ABx55JGjmK+of2+S8qbg1l+tkrQSS49b9nPqzCMI7VsptGCleWq3NM3P2bXaliNoL49VjxWWSPNGnw0TobbHwyhEfjyAc0h+s90/EsAXSX4OgIjI43KFReRE7xjJfyF5aKvGvfWBQfKBAD4C4Pe7C4pT1yY06wVjZWWl1/8aHllYg9NKyw0ia2Cn5Wun2FY9Fkl4bfF8y9Wtt3Pq0iNY77xp1eoRnibaUjvXe58NEfHUyngYhcj/3za/oHlq5Hs9+tGtD3wRnPWB2w+dfgjAu0RkELfaPfWXpuUGY20owKoz3S4RSUnJeYTp1ecdT/c9giv5rm3lFHiuTda+599677OhQADsGqZrg0fxZifJfUm+DsDfAXgnmrDG/wbwmwBuF5FviMg3JvTjIgAnkbwZzRukF7V1r5B8e5vnDABPBnAuyRvb3zET1jsyutiiF8ft4qqeOtTwVFu67w3u1Ie0Xq+s/mm/rWO6jtx2iQRzNrry3vnKlUvL545HnzXbQ3rkUEMqf4G1qFHkr0fzFucjROQeYE+I44/a329N6kTN+sAi8pcA/nLSuvpALqyQ7ntkYB3PKbAaBZibMufUX6lua1vn0e2obbtFcp6CrlX6JWKMPhtwWAURIx8XNY8f/jKA/9KROACIyL8BeAGA/2daji0CcsqtpMS8AafVVmkQpmVqB64XVrB8SLd1u3JEVCKMmuO5PFrNera1neizATOlNDHyGXzqbelQo8hFZO9TJyK7Oej/FdNDF165+oTT96R5p8JSRhYB1hJfaqOk/Dz17BGU1RaPwHLt1PmtUITlZ+7CYNVvKW3LZq58qZ5l6rMhPnKoES/7jIcaRf4FkufoRJLPQfPdunUPPWCtKXu67yElQqucjofq+q06rKl4Z9vztwTPfmczVctaEZdCCbqNWu3qdtdcJCx1v577bKjoQivTfiGI5EEkt5C8uf17oJNvY5vnZpIbk/RrSX4puVf30DZ9f5KXkbyF5HUkN0zmaT1qFPlvAvggyf8MYDua8/3zAH4MwK9O0bfBo7tpdPUJp2fVWS6+qkMEtdPgVK1Z6bq8Jjut2CzFl1OBWg17anMUBWrBa593Lr3Zg+dXrt5l6bMh39zUmNFaK92SIBeRfFm7/3tphnb9p1cBWEHDedtJbk7WgDq7vYeX4jwAd4nII0meCeC1AJ49zYZ0KCpyEblNRI4HcAGaBa2+ieaL0seJyG1T9m9hkFOSOUWqB/Yo5XLHc4rP898KI3gE4ynYdDsXB/bqthR0iehqFL91fD312aJgRjHy4pIgAJ4BYIuI3NmS9xY0azzV2n0/gKeRnMnd5VE+vvxxAB+foi/rEjWqFsgPVIvocmEIjwDSfOlUfhyS0OrYO55Tvrm2j3JurHZOgkXrs0WBYKQY+SHtx+A7bGpfBKxBzZIghwG4Ndnf0aZ1+AuSuwF8AMCr2/uIe8qIyC6SdwM4GMDOSr/GxigvBAUcnLj1ijVrsFiDqDTYS9PzXKigFL7QeawQhFc+F05J03M+WyEJLwRS03adtyYEVApDLGufLcINzhSr9XJ7p4iseAd7WBLEugp2zp0tIreRPAANkT8XzaKCuTJTRRB5z7AGnxViGEUxeerWs6vL6jyW3ZzqLRGebpNHXpY9r30lW2m+EsFZNqw2LWufLRL6Yr0elgTZAeApyf7hAK5tbd/W/r2H5HsAHIeGyHcAOALADpL7AngQgJl8fGeUZWwDGZx03eWu+rEUnbXdKcfcIOyIK/2ldtJ4bZrfi8fqv3rbsq39KJGRRTaeitfpVghChxIsG1b7tP85LHqf5f4/DhUigt2VvwnRLQkCOEuCALgKwNNJHtg+1fJ0AFe1b7ofAgAk74/mPZvPG3ZPA/Bx69HtaSAU+RTgEU8HHV7Q5Sx7lg1dR8mGhZrYsyYW7Y9Fcl5IxgthjOJz53fuQqL9sNS2ly+1b9VX6/s8+2wRIZjZm50XAbic5HloHt44HQBIrgB4vog8T0TuJHkhgOvbMhe0af8BDaHfH8A+AK4G8LY2zzsAXEryFjRK/MyZtAZB5L0jfSQxhadaLQUK+MTTITfIPVWbi8mW1Gopnqt99/ytjfV6IQ7Lx1oSTvPXkN0i9tkiPWpoYRYvBNUsCdLuXwzgYpXn3wEc69i9F+1FYdYIIp8BLBVYc6MqR1Detle/tp2rzytrTf212rZQaqtFijmfcvDUuXcB8bDIfbbImFEkYukQRD5FePFdvW+RopW3Vsl20HZ0unXMqjct18WDU+LOhQ+s7VybUwL1iLakVC0/rHNR44tla8h9tsgQxCv64yKIfErwQiwdvIHqKVzrxpdFpNb0XpODFWpI00tEY92cyx3XsEjbylNSsNZxy2aOcK32eRhyny3ajU0PocjHQxD5jKAHukdgHnl55XJkp0lVk1wpjFEKmVhlao6ViDNHdp59bzunnnV+L7a+SH22yBAAu4LIx0I8fjhlnLj1CncgWwpShy+6/Olxq5wX97XsW4O/qzMHL+SSEk8p5JH6MSrStmvCLSlsT/3r8+3NJobcZ4t+gzOFVP4LrEUo8hnAUnm1+a39nJ1cmhW/TdNyRJ5Tsta+FSrw4se59pSOeX5Y6lrX3YdvOR9q2zGtPltERIx8PASRzwCdYtpy/BlFgvCmzrn4bi40kNrwyEHXldrywg1eumUrZ7/Wb0sB678WcrHztHzpuNfO1F/dzmn32TIpcaBR46uhtsdCEPkM0d2Q6ghdE5WlBoG1oYsOuTI6Pac2reO5uPIoBG21wSpjEWbJnmcz9TkXxrDIM033yg+hz5blxuZekJHWWgkkCCKfA0aJgXbH9b4Vw9X7OYK0yln2PaXr2SypZG3bU565erSizbVd20lJPlXiVptzYYx59tkyI+Lf4yFuds4BJ269wiUCj4D0flo+JaX0VwoDWKpPE6sXdtB2NDmnti1ftR0vf9puT/16tjWZWqRvwbKhfdD5Z9VnyxZOSSEAdmG16hdYiyDyOaHmhlVNyCJnw1K5FpFon6z6tQIukZfll67PI+DUnxqFbtnW6Zo4LZsegaf+WO3y/JtGny03ap9ZCdWuEaGVOUF/wLk2Dm3FXNNta0pvKbxRpvCWLS+mq4/p8rUx71FCC7kLVZpeOsee316eWfbZMivxDgLEzc4xEUQ+Z+jvfpaUoKeEvRt/3bHctD0X782FRXR9VlnPR4s0c21L/c+dD88fnWYRcomMR51NeD6O0mfrgcD3gMAqI2wyDoLIB4LctNkKBZRuRqblvFCCVbdHMlZs3SJWK+6u0zxfPXL38uYuTHo/51fJJ8/3WfbZekEo8vEQMfKBoAu11Ax2wI/r6vId4Xl207hxKUZs+eCFBqw6Un9y5FoKaVh+ePFvna+UruvXflt5ptVn60qNo4mQ7678F1iLUOQDgqduu30grzi9EIqnstM0K4yRO275kvqdi4d7FwyP4Es+e77kYti1Nqw2ad+n2WfrDRFaGQ9B5AOCtWKiF8e14roWWXjEaIVOLEXrkZXli94uEVmpbl2vhZr4tpffu1jlYMWxc+dsnD5b2hd+Cmje7AwiHwdB5ANE+kp/h5zSzd1w08ThkYsVcvHsax8sxW8pU12+ti4NrWZ1m71ynq+li0Rap3eRGnV2YuVfrwSeIoh8PESMfMAohUQsMrPIx8ub5rFsp3WX4r/aZx2zthS9VZ9VlxebtvJYZayy3oWvhuC92YeemXj2ay/E6w/daivl3yQgeRDJLSRvbv8e6OTb2Oa5meTGNu0Akjcmv50k/6Q9di7JO5Jjz7PsTgOhyAeM9GbX1Sec7pJFSpK58IlHoumxbrs75qntGkWq4an+mpCKF1euUfC1sw59HnLnohQaKhF89zdU+H0QzCxG/jIA14jIRSRf1u7/XpqB5EEAXgVgpXVtO8nNInIXgGOSfNsBfDApepmIvGjaDdAIRb5gsBS1Vr4psXjx2y6vFb/1SC21oQm1s2FBE6N3sdGEqOuviTF7MxCPTK18qZ+5C503m7HOl963+iYg2I0fVf0mxKkALmm3LwHwLCPPMwBsEZE7W/LeAuDkNAPJowA8FMDfTerQpAhFviCwXhzKTf29uHa33cEiYy88oNV3TRjCql/DI+iSys8pjDCc2AAADgRJREFUcSue7ZWzwj7WufIUvNeGNF0fX2+PFtZghjc7HyYi3wIAEfkWyYcaeQ4DcGuyv6NNS3EWGgWedvqvkXwygC8DeImI3IoZIIh8wZAjN494LULzYrjp3xpFqstayrZDLrZs+eSFJXIXIu98aJLWilvb8o5Z27pOy+9Q33UYgcgPIbkt2d8kIpu6HZJXA/hJo9wrKu1bqkN34pkAnpvs/xWA94rID0g+H43af2plfRMhiHzBkK7R4qm9HDSJWWrVU+k1cXGP0L04vq5Tw1LPpYuGd0FIfbfK6eOli5BOryHsUOI5NK8EVWKniKy4lkRO9I6R/BeSh7Zq/FAA/2pk2wHgKcn+4QCuTWz8HIB9RWR7Uue3k/xvA/DaUiP6wiCIvL2xcBmADQC+DuCMNi5l5X0ggC8C+NA8bioMBSkhdI8pWgq0Jq6bHvcUZQdLTefye0Svj3kqNqd+LTteLDynyHNtss6NrssKzaT54oZmHQQze/xwM4CNAC5q/15p5LkKwH9Pnmh5OoCXJ8fPAvDetEB3cWh3T0HDUzPBUG52dneRjwJwTbvv4UIAn5yJVwsCS9GWVGruWEpQuVh7akMjvYjU5s/5k5t1WBcwq77cRcw7f+m2Vd5S/96xQBmzePwQDYGfRPJmACe1+yC5QvLtACAid6Lhmuvb3wVtWoczoIgcwItJ3kTyHwG8GMC5kzpaC8oAPq1E8ksAnpJMda4VkUcZ+Y4F8LsA/gbASo0iX1lZkW3btpWyLQW6T8iV4rkeGZbKlm7k6XRdZy6fBa9s7UXLy5NT0jW+eUpd+7ielDjJ7blQRw3uv88BcsiPH1POCOCfv/v3E9e3TBiKIl9zFxnNIz1rQPJ+AP4YDZFnQfJ8kttIbrvjjjt6d3aosG7Q6e1cHDnNm4sF6zRPBdfE4zW8uHxpVuCFgazYfCm+70Hbss5vWmdgVAhWZXfVL7AWM4uR93AX+YUAPioit5KlASebAGwCGkU+ip+LDB03rwmJeDdMvXCLp3Ite145XdZS0LWzA2+WURvCycW4vbBSLp6+nlT4NNBD2GRdYmZE3sNd5CcC+AWSLwTwEwD2I/ldEcnF09ctNKGkC3EBtlq2SKrb1uW0LZ3Pi0enNqzYdilcoevRZUsq2yJuyy+LzC0FHsTdJ0Z6aiWQYBBPraDiLrKInN1tkzwXTYw8SLwS+sajRdQ6LRdOSbc1AZbUu/YrR5i5GLq2YaXlQkk5Re7FyGv8CYwHAbAqocjHwVCI/CIAl5M8D8A3AZwONHeRATxfRGa2+MyyIlWOnTq3QiBdegdLPafIqeq0vHUBsC4muadCrNCKR6yeaq89lvMxngWfEkSwKhO/fr8uMQgibx+kf5qRvg3AXiQuIu8E8M6pO7bksEg8PWaFKjQsks8pak+B5y4Clq81T454diyfLTVeUvWBfjHDV/SXDoMg8sBskXuZqHQzUyNHbpaCTvfTv1asepQQR01YxMqv2+HNMEKFzwYSoZWxEEQeyKpeixBzj+F16RYx1qpwqz4vT46QrTCJ1W5rv+ZiEOgbcbNzXASRr3NYT13oJ1yAemLTsWVLldeEUrS9Uigl93RJrkxqs0Oo7/khFPl4CCIP7AWPyLpn07UKLql3TaTezdNRibkrk/610r0y8ejg0CDxHPmYCCIPjAXv5mAuH1AOnXhPi9Soe8+m9iUwTAgEq6vx1Mo4CCIPVKM2DJOLPZduXNYo7tJxK2Ye4ZLFQCjy8RBEHpgYuRuiaR6tsNNtKxRS8zihritU+AJDIkY+LoLIAxOhpHS7D2BYL9zkYt9dmt72QjQR714GRIx8XASRB6aKCGkEaiEAJFY2HAtB5IFAYCAQIBT5WAgiDwQCA4FgVXbN24mFRBB5IBAYEEKRj4Mg8kAgMBzEUytjYSifegsEAuseMpOPL5M8iOQWkje3fw908v0Nye+Q/GuV/giS17XlLyO5X5u+f7t/S3t8w0SOjoAg8kAgMCCsVv4mwssAXCMiRwG4pt238HoAzzXSXwvgjW35uwCc16afB+AuEXkkgDe2+WaCIPJAIDAQCERWq34T4lQAl7TblwB4lumNyDUA7knT2Hww+KkA3m+UT+2+H8DTWPrAcE9Y+hj59u3bd5L8xhRMHwJg5xTsTguL5O8i+Qoslr/T8vXhPdi4Cth1SGXeB5Dcluxvaj+6XoOHici3AKD9TvBDR/DxYADfEdnzeM0OAIe124cBuLW1u4vk3W3+qf/fWHoiF5GHTMMuyW0isjIN29PAIvm7SL4Ci+XvkH0VkZP7skXyagA/aRx6xaSmjTSpODZVLD2RBwKB9QcROdE7RvJfSB7aqvFDAfzrCKZ3AngwyX1bVX44gNvbYzsAHAFgB8l9ATwIwJ3jtWA0RIw8EAisN2wGsLHd3gjgytqCIiIAPgHgNKN8avc0AB9v808dQeTjozYeNxQskr+L5CuwWP4ukq/TwkUATiJ5M4CT2n2QXCH59i4Tyb8DcAWam5Y7SD6jPfR7AF5K8hY0MfB3tOnvAHBwm/5S+E/D9A7O6IIRCAQCgSkhFHkgEAgsOILIA4FAYMERRF6J2td627wPJHkbyTfP0kflQ9FfkseQ/DTJm0h+luSzZ+zjySS/1L7SvFc8cZ6vPBu+lHx9KckvtOfxGpJ9PFc9Nkr+JvlOIykkB/lIYqAOQeT1qH2tFwAuBPDJmXjlo8bf7wE4R0QeA+BkAH9C8sGzcI7kPgDeAuCZAI4GcBbJo1W2ub3ynKLS1xsArIjI49C81fe62Xp5Hyr9BckDALwYwHWz9TDQN4LI61H1Wi/JYwE8DMDfzsgvD0V/ReTLInJzu307mudpp/IClYHjANwiIl8VkR8CeB8an1PM7ZVnhaKvIvIJEfleu7sVzfPF80LNuQUawfE6APfO0rlA/wgir8ea13oB7PVaL8n7AfhjAL87Y98sFP1NQfI4APsB+MoMfAOS15lbpK8675Wnffmie+V51qjxNcV5AD42VY/yKPpL8vEAjhCRNSv7BRYT8WZngh5e630hgI+KyK2zEI59vYbcvt12KYCNMrvPmNe8zjy3V54Vqv0g+RwAKwB+caoe5ZH1txUcbwRw7qwcCkwXQeQJenit94kAfoHkCwH8BID9SH5XRKbyYkAfryGTfCCAjwD4fRHZOg0/HXSvM3dIX3XWeWb+yrPjRwfLV5A8Ec1F9BdF5Acz8s1Cyd8DADwWwLWt4PhJAJtJniIi6UJUgQVBhFbqUXytV0TOFpEjRWQDgN8B8K5pkXgFiv62C+J/CI2fs/7c/fUAjmoX6d8PwJlofE4xt1eeFYq+tqGKPwdwioiMsnbHNJD1V0TuFpFDRGRD+391Kxq/g8QXFEHk9ah6rXdAqPH3DABPBnAuyRvb3zGzcK6Neb8IwFUAvgjgchG5ieQFJE9ps83tlecxfH09mlnYFe151BelmaHS38ASIV7RDwQCgQVHKPJAIBBYcASRBwKBwIIjiDwQCAQWHEHkgUAgsOAIIg8EAoEFRxB5IBAILDiCyAOBQGDBEUS+BCB5OMkr27XHv0Lyf7Zv9OXKPLhdSiBN+z9j1v9dtX9w8oLRP7drs3f7pl8kN5D8fG0dBX+OIPk1kge1+we2+w9X+TaQ/D7JG2ttj+DDj7Xt/SHJQ/q2HwikCCJfcLTLun4QwIfbtcd/Gs0bhq8pFH0wmkW+9kBE/u8+fBKRb4vIMSJyDID/BeCN3X67rOpUISK3AvgztG+ztn83icg3jOxfaf3s24fvt3b3WpMlEOgbQeSLj6cCuFdE/gIARGQ3gJcA+M8kjyb5TyQvab9c836SP96WuwjAf2xV4+uB+1QvyXPa/P9I8tKuIpIfJrmdzReFzh/X4YydfR1f07LPIfkPrd9/3n5EwcIbAZxA8rcBPAnN8sIlvza05+vtJD9P8t0kTyT5qXa2c1xNnrFOSiAwCUQkfgv8Q/OFlzca6TcAeBya5Uv/U5t2MYDfabc3APi8KvNdAI8B8CUAh7RpByXHD2r//hiAzwM4uCuX8e8Puzpzdlp/TF8T334GwF8BuH+b9lY0Xzjy6n5Ga/Mk5/iac9Du7wLws2hEzvbWD6L5MMOHa/KoOr7encv4xW9av1Dkiw/CXhu7S79VRD7Vpv0lGnWaw1MBvF9EdgKAiKTLxr6Y5D+iWS3vCABHjemzZ6fk69MAHAvg+jau/TQAP5Wp55kAvoVmydZafE1EPifNuuw3oflcngD4HBoSr80TCMwMsR754uMmAL+WJrRrjB8BYDf2JvnSKmnmhYHkUwCcCOCJIvI9ktcCeMCozhbslHwlgEtE5OUV9RyDZtXHEwD8Pcn3SfvFpALSdcRXk/1V3DdeavIEAjNDKPLFxzUAfpzkOcCeD+/+MYB3ovm48pEkn9jmPQvA37fb96D5wIBl7wySB7f2DmrTH4TmQ8jfI/loNAQ5DnJ2PF9T304j+dDON/0kSptONDc7f1tEvolmidk/GtPfQGDwCCJfcLRT+l8FcHq79viX0XxM9/9rs3wRwEaSnwVwEBqCg4h8G8Cn2ht2r0/s3YTmiZdPtuGPN7SH/gbNzcjPovlo77hfE8rZMX1NfPsCgN8H8Ldtni0ADjXq+C8AvikiW9r9twJ4NMl5fn4tEJgaYj3yJQbJDQD+WkRGiRGvG8zi/JD8OoCV7p5DIDANhCIPrGfsBvCgab4QBOD+aGLngcDUEIo8EAgEFhyhyAOBQGDBEUQeCAQCC44g8kAgEFhwBJEHAoHAgiOIPBAIBBYcQeSBQCCw4AgiDwQCgQXH/w/t0WiWM8zz3wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "label_pack = Labels(xy_base='Optical Table', z='Height')\n", "z = NollZernike(labels=label_pack)\n", "z.plot2d()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finally, across the library the default units and labels for the various classes can be controlled with the config class:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "show_units \t\t True\n", "pupil_labels \t\t X/Y Pupil Z OPD\n", "interferogram_labels \t\t X/Y Z Height\n", "convolvable_labels \t\t X/Y Image Plane Z Irradiance\n", "mtf_labels \t\t X/Y Spatial Frequency Z MTF\n", "ptf_labels \t\t X/Y Spatial Frequency Z PTF\n", "psd_labels \t\t X/Y Spatial Frequency Z PSD\n" ] } ], "source": [ "config_items = vars(config)\n", "for key in config_items:\n", " if 'units' in key:\n", " print(key, '\\t\\t', config_items[key])\n", " elif 'labels' in key:\n", " print(key, '\\t\\t', 'X/Y', config_items[key].xy_base, 'Z', config_items[key]._z)" ] } ], "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.7.4" } }, "nbformat": 4, "nbformat_minor": 2 }