{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Pupils\n", "\n", "ost any physical optics model begins with a description of a wave at a pupil plane. This page will cover the core functionality of pupils; each analytical variety has its own documentation.\n", "\n", "All Pupil parameters have default values, so one may be created with no arguments," ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "from prysm import Pupil\n", "p = Pupil()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Pupils will be modeled using square arrays, the shape of which is controlled by a `samples` argument. They also accept a `dia` argument which controls their diameter in mm, as well as `wavelength` which sets the wavelength of light used in microns. There is also an `opd_unit` argument that tells prysm what units are used to describe the phase associated with the pupil. Finally, a `mask` may be specified, either as a string using prysm’s built-in masking capabilities, or as an array of the same shape as the pupil. Putting it all together," ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEOCAYAAAB2GIfKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO29fbRlVXUn+vtZfGjSRvkQm4A02FSSxo/4cUORTtrnCB+WDmMlHUGI0dLgY5CWZ9J5nSG0LyEDzRvQyYut0Zd0CaTBVhBQm0oamy5AXj6LUBAClmgokGghAyxAJFHAqjvfH3vvqnVnzTnX2ueee84+567fHWfcvdfHXHOtffZv/dbc++xNEUFFRUVFxXzgOdN2oKKioqJifKikXlFRUTFHqKReUVFRMUeopF5RUVExR6ikXlFRUTFHqKReUVFRMUeopF5RUTHXILme5FdJ7iB5vpH/OpJ3ktxN8q0qbyPJ+9rPxiT9tSTvaW1+lCQn0ZcSVFKvqKiYW5BcA+DjAN4I4AQAZ5E8QRX7OoB3Afi0qnsogAsBrANwIoALSR7SZv8hgHMArG0/61eoC71RSb2iomKecSKAHSLygIg8C+BqABvSAiLyoIjcDWBR1X0DgC0i8riIPAFgC4D1JI8E8EMi8tfS/HrzSgA/t+I9KcQB03Zgkjj88MPl2GOPnbYbFRVzjTvuuGOXiLxoOTbesP6V8tiufyxs72vbATydJG0SkU3t9lEAvpHk7USjvEtg1T2q/ew00geBVUXqxx57LLZt2zZtNyoq5hok/2G5Nh7b9Y+4bdsHi8oewF96WkQWPHeMtNJno3h1l2NzxVHDLxUVFYODQLC4uKfok8FOAC9J9o8G8M1CN7y6O9vtUWyuOCqpV1RUDA8iWFx8puiTwe0A1pI8juRBAM4EsLnQixsBnEbykPYC6WkAbhSRhwE8RfKk9q6XdwK4frSOjh+V1CsqKgYHgWBRdhd9QjsiuwGch4ag7wVwjYhsJ3kRybcAAMmfILkTwOkA/gvJ7W3dxwF8EM3EcDuAi9o0APgVAJcC2AHgfgBfGPcYjIpVFVOvqKiYFQgkQ9jFlkRuAHCDSvutZPt2LA2npOUuB3C5kb4NwMvH4uCYUUm9oqJigBgfqa82VFKvqKgYHkQgi5XUR0El9YqKimGiKvWRUEm9oqJigFiE7PnetJ2YSVRSrxgMtqw7A+TS33CINL/z6NJFiFNvu2bivlVMFiI1pj4qBnNLI8mXkPwiyXtJbif5q236oSS3tE9J29I9UIcNPto+Je1ukq+Zbg8qxoWOyEUIUkDK3rSK1QIBFneXfSqWYEhKfTeA/1NE7iT5fAB3kNyC5ulpN4vIxe1jM88H8H40T13rnpC2Ds1T00qf6VAxJdx00ul7yRrAkm398FKt2ru0m046fUlaauOUrdeugNcVE4dIJewRMRhSb3+l9XC7/RTJe9E8JGcDgNe3xa4AcCsaUt8A4Mr2KWlbSb6Q5JGtnYopoyNeTeAAlyjvbjsl8HQ7Vev7bMAs27VrlalkP4Oo4ZeRMBhST0HyWACvBnAbgBd3RC0iD5M8oi3mPUFtCamTPAfNc49xzDHHrKjfFftghU40cXvEnqbrPE+9W2VqyGZ2QVkEdz+dL1ixHwZH6iT/GYDPAvg1EflO8EKRoieltY/g3AQACwsLg3mS2jxBq+OGVOmq7Bwxa6WeQqv6pW3uX06k8WPLujOWlKsXXIeOGn4ZFYMidZIHoiH0T4nI59rkR7qwSvtw+kfb9OU8fa1ijLCI1lLLFhnn7Fmx93Tfujsm3fcuspb6UjEtCFjDLyNhMKTePu3sMgD3isjvJ1mbAWwEcHH7//ok/TySV6O5QPpkjaevHNKLkymhdqq8S0+hwycW4eryOt0L0Wi7UZ6n/NOLttq3GoOfMgRA/rG6FQYGQ+oAfgrAOwDcQ/KuNu0/oiHza0iejeZdgh273ADgTWiekvZdAO+erLurG5E612W8uLp30TMXQ0/reRdZPUK3FH1E/BXTgoA1/DISBkPqIvIXsOPkAHCyUV4AvHdFnVql6H4E5JFdLkbeh8zTMhbppnUt1R/5Fql73Ya2Z41BVe+ThFSlPiIGQ+oVk0d68VCTpUW0KSKlm5bR4RGPjDtbFrzJI7Wd5ltKW1+EtVYFkW8iXHJBuF5wXWGIgLuzL8CoMFBJfRWjJM6tywM+0etyJfVyd7NY4ZXchGOV89qMQki5C8D6Tp2KMUKqUh8VldTnHJa67OCRo3eR0crL3XKYu4DppWn/rPppXtq/3J04uUnGshH1JX1mje5LDdmMDlZSHwmV1FcBvPBCFK4oQWnIJhd7H8X3ND8KoXhhlpwi130pUeTeyqBiFFSlPioqqc8p9PNROuTIyVKj6ballnN3mug0j2Cj+Hzqv3fHiy6T2i8JF0UXab0x0mW0je44VMXeDxQZm1InuR7ARwCsAXCpiFys8g8GcCWA1wJ4DMDbRORBkm8H8BtJ0VcCeI2I3EXyVgBHAuieD3yaiDyKAaCS+pzCiweX3IqY1vfIq48azZGwLhO1nbsoW3JtILWZi5t79XQbXgipYkSIgLufXbYZkmsAfBzAqWh+sHg7yc0i8uWk2NkAnhCR40meCeASNMT+KQCfau28AsD1InJXUu/t7btKB4XBPHq3Ynzo7mrJQd/dkQufpNtenFrXj8hfK+LSelF+1Ke0rTTN6p81BlaYSfuflk33S49JRYLFPWWfGCcC2CEiD4jIswCuRvMwwBQb0DwsEACuA3Ay938+yVkArlpmjyaCqtTnCPqxtoB/kTFS4RqRSi6Nk2vi1O1Hd6R4ZaJQjXfXjJcWrSa8vGglZE0i3QXVGoopgYCLi+MwZD34Tz+ie28ZEdlN8kkAhwHYlZR5G/afDP6Y5B40jzb5UPvbmamjKvU5gqV00ztMdNnuv6V6o4ufXR1LnXZpll0rVm+p5HTfKxepeqtvpRdvLXtempWXi8/Xi6iFEPRR6oeT3JZ8zkkslTz4LyxDch2A74rIl5L8t4vIKwD8m/bzjlG6uRKoSn0O0ClA64GWFnlFF/YikovUdHSxMrrrxEKJzymx5y7aWunWvjce1oVXT+Fb/usy9eJpCXrd/bJLRBacvJIH/3VldpI8AMALADye5J8JFXoRkYfa/0+R/DSaMM+VpQ6vJKpSnwOkJOdB51sq1iJii+wsEtX+6G2t4i3S9exaSjyKvXt1LXJOxyL9ePajcfFWDBWjgbJY9MngdgBrSR5H8iA0BL1ZlekeGggAbwVwSxdKIfkcNM+bunqvX+QBJA9vtw8E8GYAX8JAUJX6DEPfthiRcgdLRUY2oniyV84iZmsF4MXDo3i2t1ooVejeWETxea9e2obli7ff2eguntbHDBgQAXZ/fwxmZDfJ8wDciOaWxstFZDvJiwBsE5HNaJ4O+0mSO9Ao9DMTE68DsFNEHkjSDgZwY0voawDcBOATy3Z2TKikPsOIwgdpvhV20PVT5EI1ej+KW+fCNCVtR6Qckbkum6aXkH/ORprm2Y8mobRehYIIMJ4LpRCRG9A81TVN+61k+2nse/qrrnsrgJNU2j+huad9kKjhlxlF+tN0IB+GSMvomLAXltF29QVQi8iiGLN1ATXXtg73WH31Jiwv1FN60VPvp4q8hIxLriGQUm93dMDFPUWfiqWoSn1GUapMuzxLheqQiQ4jeEo/8idSpRYhpvajlUAutt/lRUrcatP6n1P+Jfme/xZKwz6rC+NT6qsNldRnDNbP/zWxRfnpvhfzjgg+2rYQxd1zcf+0T5p8rTh1NBYlMXM9UXhq3JpYrEkjF7NP0+o97AqCSuojopL6jECTuRf66JC7MGhBhzii2Lp3YdRrK1Lb6bY14UQXHUti9t6k5l1r8FYWenKxxsnzzRsb3Wa9eNqhKvVRUUl9xmAp1ZILel65HPGUlLWItCSUYsG7mJhbVVi2S8M23uSSW11YNixVbvndZyxXIygCjuHul9WISuoDR6fQcyQL5FWqRYwlBKTLW2mWgvbSo+1IlXt99PI8FR5d5IwmjDQ9RU61l45BOl71B0qoSn1EVFKfAViEECl1XS61ExF7Scw7zfdUt5cW1bHqWpNC5JsXt07b19CEXLrK8fzNTUpRnL2q9AQ1pj4yKqkPFNEPi6wLmRoWUVtEqNO17egCX0rSXuhDx8WjyUPvW/+9ScEKm+i+6u1cmdwkFx2PKM0b19Rmh+4hbasvxl5j6qOikvqAMWoMeZQwi6dirdCBRYIlE4VnP9fvdL80Vl2i4KPVgxf28UJMGl44StvwQjxe/1YNBMDiKuz3GFB/fDRAdHdAWARqkbQmJosoLNKNYtZeyMcrm350Wi6kEk0aHsFHoQvti1XPW5V40BOLR+hpuyW29QTsTcrem6zmFwLs3l32qViCqtQHhH0nbvn91+l+REppfhQ7Ll0d6HBEzrYmUG81oftkqekoLu+tSHJlIrUdTTTeaiHy0fPPmww6e6vq4mlV6iOjkvrA4KnsEqLqUBLz7RvSSetZ7Xu2vb5Z5TwFb4UyPLtWP/pMDJG6T7e1Itfkbh1DL+wS2dV5qyoUk38CY4WBSuoDwL43FtmqNxdqsNLT+iWwiD+1EanWiGSj+LhV11t1WCEgK5TjIRqbkpWL5U9ukvLi5RFJW+OdfgdEVsvTHaUq9RFRSX0A8OLEuTqeKkzhkZ4mRI/MLB8t+x7pe2RutW+RmdXfyD89JiVqV7eXi397bVu+eGGqqP+6Pa+Pc40afhkZldSniC5Gar2xKCICq5ynCi2UhDKiEIrVrkVsXqgj2ta+lRCk5WNOCes0az+34vAmoFyYy+uH1zdvwp372x0rqY+ESuoDQKTQPYKK1JpHUt6+pfotPzxbuTBCabhIp+X6nPPZI+Cc0o1WL16budVS9790NeaNQzRZzhNEANk9n31badRbGqeETmUBZUQWhVH0Et0KEaT5fbc1GUZK1wov6LBGZ8vy3bITrTw8X9M0vW/5pOGtkDx/df/Sflp2vb7nVlk6fAXM6e2OAmCx8FOxBJXUpwQvjmulWSSkScFb1nsTgd72SEunW357RGn1KUVK7p4/Oj1nOzd+2mY0WeoykXq26nptWxOkNaaWj6NMBDOLSuojoYZfJoz0h0UlsdUUUfzZspnWsdqKbFjte8QWTS7e5OP1qavj+aQRtW31Qdf1fNNjZB2vtFxUL/Lbqq/7ZU1kaR4p8xlfr9GXkVCV+oSwZd0ZuOmk080leSlhaXjxYYsgPDLNhQh0mVzIQbeb1i9R1bq8R6JWjDxH6JYS9vyyYIVxrHRrUo7UulVOHx9r8rds3nTS6fPxejwBZJFFnxxIrif5VZI7SJ5v5B9M8jNt/m0kj23TjyX5PZJ3tZ8/Suq8luQ9bZ2PktbtDtNBJfUJIVK9Jcrci8F25SwC9trRfmh7Fll6pBRNLN6EkfrsTRgecXlhGGsVoevnJjIP3orGS0996D65FUO0n+uHdezmAmMIv5BcA+DjAN4I4AQAZ5E8QRU7G8ATInI8gA8DuCTJu19EXtV+zk3S/xDAOQDWtp/1I/RwRVDDLyuMfeGW/fNKTkDrhPeUaV94hKjhKU69AtDkHIVdtN1IxXv+eH5HE16uXt90z18vvKL7oOtbq6uSFVlaToSz/3o8AWT3WDTniQB2iMgDAEDyagAbAHw5KbMBwG+329cB+FikvEkeCeCHROSv2/0rAfwcgC+Mw+Hloir1CUEv/UuUnt7Xy3CrvhUnttItBZ9TiBqaiDy7uRWJ5Wvkh6V+rZCP9ietZ7XhhWFSYrbCRJ5vVn+9lYkHj8ytNkonxdkAgcXCD3A4yW3J55zE0FEAvpHs72zTYJURkd0AngRwWJt3HMm/Jfn/kfw3SfmdGZtTQ1XqK4T0h0VRCMQjH09VlqhqTzFa7aY2PRLS5fT/tJyVHhGp14eoXxGB6zajicRTx3rbqm+Fu6wVTcnEE4WkLKRj4YXdOsz0Q8DKJ6VdIrLg5FlG9EB5ZR4GcIyIPEbytQD+O8mXFdqcGiqprxA8JewRkFU/zbdIJVWKEelZPnk2U3jKVefnFHZEbNqONwlZ45Uq6HQsrEnHU9jeaicajxzpeuVLwzhWn0vLR+p9ptBeKB0DdgJ4SbJ/NIBvOmV2kjwAwAsAPC4iAuAZABCRO0jeD+BH2vJHZ2xODYMKv5C8nOSjJL+UpB1KcgvJ+9r/h7TpbK867yB5N8nXTM/zpeji6JZi7rb1ia4JCYhj2RqenWhC0fW7Nj0iTNOsUEakhq2ylm8loQlvcikJa3k+pna0f9YkYIWzotBNWs47rqltbUenp3WtdH1MZvIHSovPKfvEuB3AWpLHkTwIwJkANqsymwFsbLffCuAWERGSL2ovtILkS9FcEH1ARB4G8BTJk9rY+zsBXD+eTi8fgyJ1AP8V+19FPh/AzSKyFsDN7T7QXM3urjyfg+Zq9FRx00mnL7lt0SO8KN1SdFoZWoToKUcr35pQPF885egpzoiwdLrnu0emluL34JG7LuOtAHQ73jHxxsRbiXjK2mozmuAs9R6tTkQ4W7c7CiF7nlP0Cc00MfLzANwI4F4A14jIdpIXkXxLW+wyAIeR3AHg17GPY14H4G6Sf4fmAuq5IvJ4m/crAC4FsAPA/RjIRVJgYOEXEfmz7h7RBBsAvL7dvgLArQDe36Zf2S6RtpJ8Ickj21l0arCIwYv16hPRI0XrZPXq67ZL/LWIJt32VG7apmcnl2/laf+sPlvw8i0F75GrVc8be8uONU6WbyUK3vtORN8Za0VRsvoZJPIqvAgicgOAG1TabyXbTwPYbykjIp8F8FnH5jYALx+Lg2PGoEjdwYs7ohaRh0ke0aZ7V7WXkHp7JfwcADjmmGNWxEH9PPSl7fvEZp3ouq6Gp/itJXhEVJZdz88+aZ6PmsxKCSsaC51u1dN9W87xsfxP63nkGU1ounzJBOuFZ6w2rRXDLNzuKOOLqa86DC380gdFV6BFZJOILIjIwote9KIVd8qKfQJ+yMKqo09UiwRSFWbZtghY18/5pOvp/NIVgUdEHumk+V5owvLRIm9LEVv9TOuk4xv127KdI3CrnoWUwNN+WhNIyaSQ9rWk/emDfW5prEgwC6T+SHuzP9r/j7bpJVe1J4aI4KKTKCLV0no5ks4RpUWEadk0zSKu1EbfvkQTi2fLsqHtWWpdl82pX73t7eeOoeW3NbZpueiYaF+tvGgl0GdCnibSyTX6VCzFLJB6emV6I/ZdZd4M4J3tXTAnAXhymvH09GTSCisqn9bR6Wm+bsPKT8vklL+u75GI50tnJ3dS9T3ponHQk4NFviUhG6uONcZR/ZKVTwlx5iZtPSHn+mJN2LmJaZAQjOvul1WHQcXUSV6F5qLo4SR3ArgQwMUAriF5NoCvY98FjRsAvAnN1efvAnj3xB1ukVOP3X4u1pqWs5b2VnxVt28pMc9O5J9FcpHC66MArb7oPGslkEvzxkErWc9Xj9xLVgdeP738aGLN9cvzVZfxxng2QCxm7mypsDEoUheRs5ysk42yAuC9K+tRGbwTzSKPEvLzlJlF5stRZhG5WnW9UEDJkt/rl+d7brs0dq3zLJK2SD8i7dwk6fmWm3y9NI2SlVrartW/waNT6hW9MShSn2VoleURRURKHaITPi1vEUWksL0Tui9peiGBXKggUv25/qbp0Rh442j1K9pO24rGMkf43lhEY+LVKykb1dfjNnTUu19GQyX1ZaC7Naxv3BjIx4Bz6qvUlqVSrckg3fZI0lOllr8euZRMMlG/Ip89Mta+RvZSv/uGubQvfVdUUd89/3SaHpuorSHf2ijofz2mokEl9WXC++JZBJNTtaVq2CPvbl+f9H0nDk2SGqVKNdd/7bdXP7ey0T5FBBfZ0P6k5GzZ9up46bn6ui/Wvlc/KlMyfoODsIZfRkQl9WUiUmxe+Q5eSCCqb9nwfIkmi1x7EQFYStybtLzJrbRdb/LwJhPPD6v9KByTtuNNsFGatheNUW6llbMdTYzeKi1aBQ0FNfwyGiqpLwOekrbyPQUXlbPKW215hJWm5foR1bHqe8Tg2fH6HNn0xsPqaxTasMbF88Mr401C3jhYE0l0HLzJJtcP7ZN1vKPJN/fdmBqEkD1rpu3FTKKS+jKQU3ae8tJL5egk80jZmiS88iVpHqF77fdZoaR9juxrMisldCvNm/Cibct/zz+PKK1jWzqB6fZzYRQr7OMdq6h/Q0VV6qOhkvoykCOG6ASLSCu16Z2wJQorIl3PxxI1mbMR+epNBhYBer5EirOkX33GtEThR3l9JlSdH02AuTBWDiWrh2lCMOwJZ8ioVyJGxJZ1Z+xVnh1ShZb+T/OjGKlGLsaa2izxJVL3Xrt92tTtWOW8E1WPRV9/o4/VF089W/56vkXHuWvDm7z7EJZVNlq9Rcfc2h/ks9alUeoln4qlqEp9RHhq0VNZXX6KXFmLLCIVmJ60XixW27PsR+o0Nyl5ClTXy/WlJFyj2/EI2wuhWP3WyB1nbaPUP+tY6f1ovK0073hF36+hKnWAEKmacxRUUh8RnvqKiDuFDmP0OcFypKQJIwpJ5Gzp5X5azyK7nL8WcXnt5HxN7WsfIkVu9SeXrxVuSQjHO/ZWWd1Pb1/7aNnzfC+Z2IeE3AswKmxUUh8RHlEDNunqMhE8xazb9khH+2L5ltax0ix/Ujs5IvRIM1L3uTCD7rdVP1ph5CagiOj7HDcLqQ+lKyTPpndM0zyrXOqH9mlwkHqhdFRUUh8R1nLbUszRiRedbGk5qz1vArHIWrcZKbeS/qYoUX9RHc93r11vnDwyzo2XtaIpIbpSf/Rxy6l7y29vdZObDFL/or4NUa1LDb+MjCypkzy0wM6iiHx7DP4MHvvecmSrxtyy3kNqs1Rt6nwNjzgiFVpix1OG2tcS8o0UdLSdjpPnc8kqx/OlRDFHytojz2gitVZfKWlHE7vXF2/S12ndxdIhPTJgXEqd5HoAHwGwBsClInKxyj8YwJUAXgvgMQBvE5EHSZ6K5imxBwF4FsBviMgtbZ1bARwJ4HutmdNE5FEMACVK/ZvtJxrhNQBW5l1xA4N1Ynox0NwJluZF9b3yUbjC8tNTbSXL+dzqQfvotWmViyZJr18lPkSrn5IJQZfzVHfpeHgTmjfe0WQbTUzRuFniY5CQ8awgSK4B8HEAp6J5sc7tJDeLyJeTYmcDeEJEjid5JoBLALwNwC4APysi3yT5cjQvrz4qqff29l2lg0IJqd8rIq+OCpD82zH5MxOICNNCRAjRiaXtWydl6QRg7XuIyKLP0t0bk1y4JK1bsqqwCCuyFSlcz7+IFHMqPpqErT6lq5w0rdv2jmHuWHj7uRXEtDAmf04EsENEHgAAklejeWl9SuobAPx2u30dgI+RpIikvLYdwHNJHiwiz4zDsZVCCan/5JjKzAVKl+cl5Jkj2YhAS5W/TrMmlxJi9ZRkascKF+h+eCGGaNLwSFj7lluJaJve2Edq3PIxInSrXDQhWPuRXU956/G3tkf5zk4KAmKx/DEBh5NMFfMmEdnUblsvqF+n6u8tIyK7ST4J4DA0Sr3DLwD4W0Xof0xyD4DPAvhQ+46HqSNL6iLy9DjKzBsspQ74RO2d+FY4Q8MjPO/k93y1SKpPOEETduRratvrm9VuRDSeuizxQ29bJGfV1fVyx1lPbGn9kkncG9uI6HXbUd9ydQYD6RVT3yUiC06eZUQfhLAMyZehCcmcluS/XUQeIvl8NKT+DjRx+amj+PIyyQWSnyd5J8m7Sd5D8u6VdG6I0ApIq1hvqW0Rs64fEWvadmRbE6lF4mlbmiD0RxOubkeTVUR02k6ap8fCKmv1Q4+LHl/L/5KViNeO9k1PWtYkGPXHW4VYdbxxto6zHgPLF2uCGxLS4xp9Mih5Qf3eMiQPAPACAI+3+0cD+DyAd4rI/ft8k4fa/08B+DSaMM8g0OeWxk8B+A0A9wBYXBl3hg9PTVknt7Wvy0dpudCE5Ve6rU/aknxPzVv+lawUon7nlLxuO/LNs23Vj8YqTcutoqLJXBOwPoa5FZNVR4+X11dvDHPtDg1jWj3cDmAtyeMAPATgTAC/qMp0L7f/awBvBXCLiAjJFwL4HwAuEJG/7Aq3xP9CEdlF8kAAbwZw0zicHQf6kPq3RGTzinkycOy7lbHZzy2P+yjXCKN8sXNL7hxy5KDJx1PtqS99iWcUxa/LaKJO/bDsRcc08ssbF10mN7nmwiDWJKXHIee351+X1n3PT73tGtePiUDG81yXNkZ+Hpo7V9YAuFxEtpO8CMC2ltMuA/BJkjvQKPQz2+rnATgewG+S/M027TQA/wTgxpbQ16Ah9E8s29kxoQ+pX0jyUgA3A9h7sUBEPjd2rwYKj7SiEzinZj173nI8R56RDzlS8vyKSNwaG4tYI3vROFh+Wn5ECtcbH6u/0Zh5E4ClqK1+WD5GbVskXzJOGt73qXRVOQ0IgMXF8TxPXURuAHCDSvutZPtpAPs91UxEPgTgQ47Z147FOQDj/i1QH1J/N4AfA3Ag9oVfBMCqIfUcAUcE1iEiptRuV9ZSox48QtX5uqz2xarr+VNCAF7oRNu2bOZWHd7EUtpOycTihU70GFq2dP2SY6T3rWPiEX9ugtb+W74OgdQBYHE84ZdZwFh/C9SH1H9cRF7Ro/zcwSJCi+SiEzFH2J7i9YjRIpUciebUdGfLasfzPSKT3Mqj1LdoHHLj5vmqJzXtV04pe+o8/Z9TylEbVp7ud+Rf9H3RyImGiUJW1WN1x/pboD4PV9hK8oQe5ecKKdFZJ7a3bO3qRQSrl8BeOb0fEYVVJ02LlFvqV7rtTWSej12+LqNVbTR2Uf8tP3V9j7StCTb1NSqv+5U7vlYbVt88MvaEgtW33PdI50UT6DQhwH7j5X3mAGP9LVAfpf7TADaS/BqamDoBiIi8soeNmUP3TIzuAqmn1FKUqMMSWIquVIHqk92bkKyTwiKLqH5EuJHq1H3Uad6YWJOQtZKwyNLzJbUbqVlrcosm2KgP0YQTjYW3iogm8ZwA0Ptb1p0BAFO9YDonhJ1F9zsfkgsAPgDgX6Dh5r0c2+e3QH1IfX0fR+cNUYgBiE8+j0lzrnMAACAASURBVIh0uahdjwQsUkjb8Wx7IYeI3LxJw5tsvFCH157lX9oHj0i9/peQv4bXpjd56XLWaiH10SrvkXOO1PSx0vVLy2t/vZXTpLFaSD3BWG4bLyZ1EfmHURuZZVjkVXICWeSWW3ZrErLI1VLZOdLVfuVUvUeClhrWiEhN1/OUvDfO0aogR/zajm4nmiSjMdfpng3r2HrjrfM9X9O63uokV96zO22IEHvKHxMwLxjLbePFpB4tDZbrxKwgWhZ75GEp1ZwajtRkSZ6nfD3fvCW6tTqIJjavvJXnEYo3YfaZUNM8z36Jmh1lvCLy9RS/7qOXllPXHryJIvo+R+1NCtNufwoYy23j9RelhfBUaqRySgjZI2jLtkU4Flmm9jzFbhGkReTeisPL033xyntKUvct9c3zX08gfcYhUv9WmmXbq+eNjwc9Iep2Ozue795xsnzR7aU2u7xpK/dVSOrvxhhuG6+/KO0BSy1aKiincKOTDyiLHVtquMRvj/DS/yW29H5OfZaoyUitW+OR7kfk6Y1RbiLMkWOuj9Hx9vod+aARCQJrPD2B4LUxVVKVVXWfeoex3Dbe55bGC0leSvIskv+2+yzXgSFjy7ozQiJKT5DcUtbK805aazJIT9BSUrDspD5YYQbvpPcUnWXH2/d8t8g1TbP8jBRm5Js1bukYaV+stnP2+hByOi5WWW+CssbFQ9c/y5b1vU3/d3d/TRrN6+xWzS2NHcZy23j9RWkhcoqxy/PUmUdeVnkLlnpNty2yySk4y340AVjpaV1Pxefa0v3TffB802W8sbBIPvLTU+M539J0z3dLzUdtWP54vnmTttU36zsYtTkNzBlhl2Ast43XX5QGKPlS5xS6deJHJ7I+wSxlb9WNlvhp+dyJby3RrRPfUqO5pX1uLNLylj/apxJ49S2/dX6aFk2Q1qSm+xH5ZfW5xEbqn+WPN7bWKsg7htMk9j2Lq+7F02O5bbwPqW8leYJ6t99co6/yBGz1bdXPkVea5pGxpaJL1bkma48IIkL1VibaT2+SshCpSGvSSMulfbPasvpcOhHnFK43SaR5egxzZXV66l/UL10mtxL0juE0ITIMPyaJcd02Xn9RGsA6AXNqxgotWPa8JXEJWUQhiG5bk7H2zetrrs0Sko9sW322bFl5JZOD549HZta4WKsdazy135Fit8rmJklrAikRCl5/rXGzfBkGuXPVXSgleTCa1+Ydi4SbReSiPnbqL0oD6JO65IS0YJFAmu6dsFpR5kjTU6Ceii9RlVYb2veo/1HfdBuWHU99WmVTUrX8t/zS29qWbss7xrqv1kRhtWVN5rmxjY6RVV6n5VZllriYBlabUgdwPYAnAdyB5D71vqi/KA3gfemt/UhhAv1jwzlC0jYj8vPKl/oS2YnCBX0mvoiILeLuMwbRJBfZ023m+uSNbx/fdLu58bCQWyVYE4vVZsl3YiWxCkn9aBFZtnjOkjrJO0XkNcstM0vQbzmyvlwWcedUb2rH27bqR2St60YEq216y+xIGWo1HPU39TFVuZYPegJMbVhhCKuv1mol6ov2bxR4E4DOtyYGq77VF893q82SibdkP93ubu09Zeu1pq8rAZFVeaH0r0i+QkTuWY6REqX+rzIvmCaaF7XOFSICStPTslF6rlxErKVKriMz7b/lQ9+VgEZuKW/5mAthlIxbZN9rLzeW1liUHkedXrr68spHbUXjEKnvPuOZ+jRdpb76Yuporlu+axK3NP5YQZk9fRodJ0iuB/ARNG8GuVRELh6HXSvOmVPUuf2IyD37USxWw1NbkeLLLbc9koxgEaI1aen+ef0t6buV7k3GVp4+RrkVgdVmTo1b4xe1n9r1SDZS75ENzwc9NqXHfNwQABK+CKgcOY5oL1BeieYVdY8BeJuIPNjmXQDgbDQc9z4RubHE5oh44xhs5El9yLF0kmsAfBzAqQB2Arid5Obl3nYZnTxWKCIiaEuZlpzU2lZuYoiUlWdT50dqO7JZ2m6K3OrECqV4PuVWC1Y5bwVTUlaXjybjkn55/3Vdq62SPqfHNkfoVhuTJvRxtlvIEWcDeEJEjid5JoBLALyt/XXnmQBeBuCHAdxE8kfaOivAO+Ph2lkPWp0IYIeIPCAizwK4GsCGlWjIIkJrOy1vQZ9kkU2LrEqW0h7JWOo49cMqY8Ga1Kz9lLx0OW/MPD9z422lW5NZup225411eryilVDuv0fyus103+uzddzSctFkqtuzyusJJTdxrBQWhUWfDEo4YgOAK9rt6wCcTJJt+tUi8oyIfA3AjtbeWHmH5J3jKNOhzy2NQ8RRAL6R7O8EsC4tQPIcAOcAwDHHFL23FYC/xM4pq6VtxwQcKbY03/ItN2lYaRZx5YjZW8Zr+3psrHY9e3pMrDxvlZJTtpYfXpslffWQI/KS45tT0RHZRmMR9ctqv0+/Vw69nutyOMltyf4mEdnUbmc5Ii0jIrtJPgngsDZ9q6p7VLuds9kHY71u2ed56s8D8FY0b7Te7+n1fW+QHxOso77k29ge3E0AsLCwUPRNjUIG3f+cKkzT0zql5KlJIi3jTRbeSWnZ8nzLhRUiWx65RqrUgxWmyCniXL+tctFxtI5VyWRu9TUaDy9Pl/Hsen3U5UuUd26imBR63v2yS0QWnLwsRwRlvHTLseUM0FivW/ZR6psB7AJwJ4bzPPWdAF6S7B8N4JvLNVqillJEatRTjTqvj/LW6ZFvlirUxBQt8y3S033PbUd1PFjKNupPml+ikksIOSVBrfJLVGQf5R+p6j7lSr53up9euWkReofF8VwoLeGIrsxOkgegUcWPZ+qOjXfGfd2yD6kfJiKnjrPxMeB2AGtJHgfgITQXNX5xuUZL1XFJaMQjyUi5R8pTE42nRq2+WOVyy/KcPavvFsHkSMVS+ZGyjXzTxBdNstEqJbfqiI6/t2rwylv+WHVzY9jHrlfOy5skBGO7QFvCEZsBbATw12iiEbeIiJDcDODTJH8fzYXStQD+Bo2CHzvvjAt9SP3PxnFj/DjRxr/OA3AjmpDQ5SKyffl2R7ubI/2vy0eKq0/b3qrAUrCpb7n2Syctr05f8vfsRqQT9dtqO1LmnnL3VgWWPxEBl6wk9ASeli1R9dp27jtrjWMkLvQ4TBbjuU/d4wiSFwHY1r745zIAnyS5A41CP7Otu53kNQC+DGA3gPeKyB4AWAneGRcoUjYTk9yOZqb6Gmb0gV4LCwuybdu2fMEWW9adsWTfOlE90iklqb7KN5dmnYyRGtbtRGlWH0raKkHkg7UdtWUp9XSscuNW4lOJn7k+RJNqzjfLfppWMj4lNkWIU2+7xi1rgeQdQYy7CGt/8DD58I+V/WL+Z+/89LLbGwpIvggARORbo9roo9TfNGojs4roRLNUY6R2I3UXwVPgnmLX/lgKPfXHUqzaX69umqZt5spbvnr9StNL1G2koCNEqluPVemKJfIzrZMTB1H9yDftU24SGGXcVgKC4Vy4W2m0t09eCOA8NGL5OSR3A/iDUW5AqQ/0ChCRrz5RvLIWAXmq2jupPOKMiN3yM0ceOs2z6534OdK3JhjLrteW9j9S17kVh554tP3cWFr7o6wIvLHS/Um/Y1HfvLHyhIc3QVm+TBSyqp798msAfgrAT7T3w4PkSwH8Icl/LyIf7mNs1YzacqC/4B55awJLicNT8Z6S99qzylsna59+aGLzTuxIuWoCsSY6j4yjEISG51u0OrFgjZtF/HpVpvOsul77Xv9S39OPN3Fb/bZ89Y5N+j+yPy2V3kHAos8c4J0AzuoIHQBE5AEAv9Tm9cKs//hoRRGdWDocYKkiy462r2167af7kb9WvUjdespN2+vy0v/Wtqf6olWFp/Q9UokUaq7fJX54Y5Oml05Yus3cysnK9/zxFLt3HHL9LxmTSUFW1wO9DhSRXTpRRL5F8sC+xqpSD2B94YFYVev6qZ0cGeo6VlnLtrdCSP31fIiUm+VLV0ar3NxKRJex2rHsaAXqEU6kqL2+WGMR5afbljqO7PUJfej+RKsnr5+WvfS/953Qq52SVc9KYVHKPnOAZ0fMM1HyPPWnAPPXUgQgIvJDfRudFXhK2SPoVLlrlJwkngK1TlRvQrFOfqu+p4qjFYPlp2UjUsVeG5YPerLT/fWOjeV7NAFZ5bzxsPrgrUwsv6NVhKe40zY8RN8HS1x4Y5kTE5PEnIRWSvDjJL8D7O1w94UjgOf2NVbylMbn9zU6Lzhl67V7X5ihEZ3QpWXTNE/Vp8opWmLrtq2T3CPhNL+kXlrWKme1YfmWI2QP1gTj+Rr1SRNctPLyfCshQe+45VYx6f9oZeOFZSyf+kxGAHrfzjguNI8JWB2kLiL7PXZlOagx9Qw88vJOwkill5BPatOCdwJ75BXZs0jeUolWnXTfU9XeKkO3pfth9c3qQ7R6yI1ntBKyyE2TqzUxemQajUVuUtI+R2res2v1Kfe9ib5Lk8KYHhMweJB8LoBzARwP4G40P2baPaq9kvDLX4jITydhmHSk5zr8AuRVT/c/p6J1vVSt6cnAOsHSepZflgqM1F33X5eNCCOaiDxy81YTkf+lZb1xjWCp8uhY5Pqv+xkRqLan9/WkaPWrZFKxxsEqF00snkCZFATTD/9MEFcA+D6AP0fze6CXAfjVUY2VhF9+uv2/KsMwEeFYylCne6Sm83UdC56qjPz0JhlPqUb99vphqT6rTmTfI0JvtWCRVrRCyPW95FhGPlh1cz6n9TxforHVtnNtRGOhy1n7k8WquvvlBBF5BQCQvAzN82VGRp9H7z4XwL9D8x49QTOr/JGIPL0cB2YFJUpWl9Np3XZk20JJe16dXFnvxI38jJbsubpeOzl17OVZpJeWs/a9MShRwx68Put2tI8lilr75H2XosnKs2XZnqZCTzEMLyaC73cb7bNqlmWsT0z9SgBPAfiDdv8sAJ8EcPqyPBg4vBOzy4sUW1o/d6LklK3Xnq5r+Va6LE/LeKrRI1HLhudbt1/alrbt+WL1x5sEtD3LfqSKo4kn1xfLJ2vlU4LcCqdELGh/LT8nDQFWk1JP734RAD/QhruBEULcfUj9R0Xkx5P9L5L8uz6NzSJO2XotAOCmk5q5y1uWpyd7mp8LqXR1uzKly+iSyUSf2Nq2p1C9MpaaKyHNSE2XTjS6T7pfXlu6rBd6sMYm11bkc25stL3ceFv+emlRe14b6b5I/4d4rQT2rBJS7+5+aV9EdDyax97cP2oUpM+Pj/6W5EndDsl1AP5ylEZnER6ReWSQU5xeG95+jgTStiKS1NsW0VlldF/1JJLaLiGfEkLPjWFJP62JKxqPkkk4shettCLfu7olqw+vHSu0ZI1zNBEMBSJje0fp4EHyAJKXAPg6gP8K4L8B+AbJ/7TSvyhdB+CvSD5I8kE0D5T/30jek3m/3tzAIrKSE9Ej6+gE9CYKq23Pntd+SV7aVrftEXhfIrPCLX19z/VJrzh0PWulFdnKlbfydEjFG8+SicRbjXjfyXTbmyi9OqUiZKUhhZ85wO+ieSfqS0XktSLyagD/EsALAfxeX2N9wi9lDzeeY0TL9xLVWRqySPPS8qXtRSQRqTbdri4fTTgeohCCbq/Et1w7Ub62GY17bmWQTrJ65aEnlNxxKu2T/u6VhoC8yX/oin0eVHgh3gzgRyR5uYWIfIfkrwD4Cnre3tiH1DdaiVN64fRUEKkzj/C9+KWX3uVZCtMiA8tO5EvphBIRTx8yjEIIUVvexGEdAy/k5Kn1qH1NfhZ55myOMgmWTmh6gtR99L5jaXm96ovEyrQgWD3PUwcgKaEniXs4wsHoQ+r/lGw/F83scm/fBmcV3YUj67EBmri8sEoH68TThBwt+SOCsIg9UowWIfQhFN3fiChyE4SemKwy1n5Emt5katnR6daEq8tZfcz5q/3sO9FHK58cuecmo+7GgOnDfkDbnOLLJN8pIlemiSR/CY1S74U+L8n4f1SDv4fmha2rCjkVB+RVlnVSpsrJqqO3vUlCE4JHKLrtEsXuKb7S/pbU88YoIjGtOvX46H1P4VrlrGPsTT5pW2ld73uStpdb/UQTbzS2WplbE0nJambSEEzm7heShwL4DIBjATwI4AwRecIotxHA/9XufkhEriD5AwCuRRP/3gPgT0Tk/Lb8u9DEyh9q63xMRC513HgvgM+R/GUAd6Dp/k8AeB6An+/bp+U8++UHALx0GfVnEt7JlOZZKstKz6msiNQ0mXkKMVqORysKr472L7KtSdEjlT6TS66Ot+LJ2bT6UxoSsWxGq6QS1e1Nyn3GIbLv2RwSJvRY3fMB3CwiF5M8v91/f1qgJf4LASygIdw7SG5G867m3xORL5I8CMDNJN8oIl9oq35GRM7LOSAiDwFYR/Jn0DwigAC+ICI3j9KhPr8ovQf7LjavAfAiAKsmnh4ht6RNkQtxeHU9BWZte+14hBD5GynEiHA8GxFRWkSWI3yLyKy61n4OOULW8FZIpQQfhV9K/IsmC2+CKyH9aWFC3mwA8Pp2+woAt0KROoA3ANgiIo8DAMktANaLyFUAvggAIvIsyTsBHD2qIyJyC4BbRq3foY9Sf3OyvRvAI8t5ktgso+TE9pRSRFJeyMBTxbo9bwIp9blvPzvkVL0ul7Od8yM3blG7lr9emgVvnLy8Lj0i1c6Gd3xTG149K69ktaV9HAq6+9QLcTjJbcn+JhHZVFj3xSLycNOmPEzyCKPMUQC+kezvbNP2guQLAfwsgI8kyb9A8nUA/h7AvxeR1MaKoeQpjeljIe8BcNlqJfMU0QmwnGV7attS3zo9p26jNG07p8o9ezrdCwlFRGuVt/pu1U0ntWhiyY2b55euEx0Py4Y1CXt98sp6ZTwCT/PTSUWneW0NAT3uftklIgteJsmbAPxzI+sDhfatk33vYJE8AMBVAD7avlsUAP4EwFUi8gzJc9GsAn6msL1loUSpX4F9j4V8I4ATsIzHQs46IlXoEUVu+Z87oSLlX6r4o2V2ycmcU3vW0j9a6pegZJWR+ubZjiaBEuRUbToJWpNc1O9oIrV88MI1Vt8sgrcm4KGpdABjvVAqIqd4eSQfIXlkq9KPBPCoUWwn9oVogCbEcmuyvwnAfSLyn5M2H0vyPwHgkhFcHwklpD7Wx0LOOrpbG7esO2Nvml5eWySvlRKw/4lsndgeWVrE3pe8PdWoy1h1PTvpvhca6qusrXKe6syNUdS33PjkVivWsfTq5trRiMbM62+alxvLITzrRWP/O7dXBJvR/Abn4vb/9UaZGwH83yQPafdPA3ABAJD8EIAXAHhPWqGbKNrdt2CCt3+XPCZgyWMhV9CXmYM+QSxS8MIG1gnukYo+adP/kQKMQi3e0tvqRxQa0GkaJQRrkarVX8tfa9/zwepTiYqO0rVKt9rWPngTkFXGC7VY46f7E01iadulK5ZJYxEs+iwTFwM4leR9AE5t90FygeSlANBeIP0ggNvbz0Ui8jjJo9GEcE4AcCfJu0h25P4+ktvbhx6+D8C7lutoKUqUevdYSKCJLT0vfUzkvL/5KEKO5EoUZGrHIn2PFKyluE7LEZmV76nkyJaur/M8xe6p59wKQ8NbyXj1opVDWqbL06TrTdTWMUnTLQLP+ZpbXaT9idS4Vye1OSQIJnNLYxsmOdlI34ZEfYvI5QAuV2V2wo63Q0QuQKvmJ42SNx+N9aWo8wRvyZw7QUtIK1pmWyTm2c7ZL1Ftke3UL8/v3H/dtyiUkfPHm1AsRGWi/VEmJs9W6aojGpPc8Y18yPk6TUwo/DJ36POUxgoFK/RQcoLklukdPCVvhSO8UIXnZ6mvli0vdJPajWx7fS3xMwr3lB6PVC17ZUrH0VPuUWgszff65vns2Upt9vm+5NqdHspCL6vl5dR9UEl9RJx62zXuiaZPLE26aT1vWd/ldTZSlKpGi+w16VjkotuPyGIck5il7NNx61u/hKQ88vVs6WPUdyKwwjjWpO2Fgqz+esfX+r54Ew8pg71IuqfwU7EUy3lMQAXsWHmXHm1bS2dN9npf28itDqKwQKlC9CYDTR6pbS9cUlIvLR+Fe7r8XOikRA1H4afoeOn63jhY/0fxz2o/mlitPP0dGKZKbzChxwTMHSqpLwMly/Hlnoze5OC1ZRFKWtcikdITP/LZI/F0PyK7yEcP0USSkrRuz/M/miysycfqe1TWO16Wfx65W/3TtqzyuYlyiJgdT4eFSurLQHQCalKJlLaHUhXlkYE10UTEaZGwtVIo9SUX0knLe/Y14Vnk6/VBE59FnDot1+eSVYI1kVtteIo5EgJdujVxepODtRLITXLThmBVvSRjrKikPgZ4S/YoDJHbj05M3ZZHclH9SNVGfUvLRyo07Zflp2XP2vb8sNS/l6f7mJbP9SV3zKw+RSswK82b/C1Y49SHmNN2osl0CKh3v4yGSurLQKSCuv8WYVknYenEYKFkxVDiv9duNMFEZOYRouW/NTZRm7mwUdSG9iUXusiFYqKy2ge9HfVL++Mp93Q7Oj6WL0NW6vUi6GiopL4MpI8M8JQuYMd3dV667yEKO6TpfcjOC8uUhBdK2/DKRCGpqB2dr5WuDkVEbaT/vTZKCT1qwysfpVuwJqZoRWP50v0fzluObKyi19mNFZXUx4CSE9IjyygsEu17pNu11ZXpq8QiNemFMCxYfcz1zUvL+Wf5oom9T9t9Vh8e6Y9C0p4/fX20bEffl0FC6t0vo6KS+hhQsqyPTkJL/eWUprUyKCVcy57lj+V/35h9rm9acUc+lqhQnV4a/y/x1+p/ujrI+eatPrQfOUXv+a7LWYQ+5Bh6Cmk/Ff1Rf3w0BqQni3fS6PyODHRc1FN9Hrl5RK59strSaal9y1bqi6VYtZ3Sycobo8iuNRa6jia26H8UztL5uYlG20j/W3n6WFjtesLBWq1Yx17nz4JiX5SyT8VSDILUSZ7ePtFskeSCyruA5A6SXyX5hiR9fZu2o3234FShVZF1wnsE6Kl5fcKmhGiRss5LffJOfk9h67zUtlU+GhNrfCz7ue3cJGGNu0fG2qZ1/PQ45I6ftyLT9rxxsMY4UvlRv3U57zswZIiUfSqWYijhly8B+LcA/kuaSPIEAGeieRnrDwO4ieSPtNkfR/OozJ0Abie5WUS+PDmXl8I7oVNYBGUpaKu+RWZpuqc8NTl5pOfZtep7Plo2o75E0H55/bfsWmETXc6yY/kbEbVlM2o7slVyvHOTZISS7+eQIAB2D9vFwWIQpC4i9wIAud+JugHA1SLyDICvkdwB4MQ2b0f36iiSV7dlp0bqpaSbklR0MkaEGK0CujKWyrZUp6fgdb9yfffCNbn+pGkeyUXhFA1PTad2LfuWPau8Pn7aP2tb+6ttWvu6L96Eq33x2rLyho7Z8XRYGASpBzgKwNZkP33hq34R7DrLAMlzAJwDAMccc8wKuLjv1rCbTjo9XEpby/79/d2fGHSaR0glqrSkbAkiex5p5sIEVl3PV51XMhHl2k3LRETokbiuH008OWJP2/bGLZqgo+Mz9FsZAUzseerziImRevTyVxGxXiEFwH3hq3UtwPwKtG8V3wQACwsLK/o1yS29gf2J00vv4JF8zo8S36LYbUTwubCCl25NdDrPUuhWm14/Lb+jScJSxTpdr3Q8BW35F5G8pbgtH/S4RX2MwlIzhRovHxkTI/Xo5a8BdgJ4SbJ/NIBvttte+tSQ/hipj/LzwhdpPS8MoLdz9tLy3sRj+WaRmdV21K5nz8qzxilaweTGyLJh9bUrY41RbrKIiNabUDWZpxNHNBHkxlH7RcpMKPQU9cdHo2EQd78E2AzgTJIHkzwOwFo0L76+HcBakseRPAjNxdTNU/RzCdLnU0fk3qVHhO4pSQ2tgC0y8mKwmsS8kIr20SOPKIRg9c0KY+iyufHxVK3lv+63N7bWuHVtRMfTmhy8SUjb0GSv+5I7xmm7qe+zRuhd+GWlb2kkeSjJLSTva/8f4pTb2Ja5j+TGJP3W9i68u9rPEW36wSQ/096ddxvJY5fnaTkGQeokf57kTgA/CeB/kLwRAERkO4Br0FwA/Z8A3isie9oXYJ+H5i3f9wK4pi07GGiy0QQdLZM1MVoka5FZ5Etu6R/5ldqISFTb0aSrP5bKz42LRXoWuXljmNrxiNrray60EZG2NyFYedFEoSc1PVFbx9o7pkPHhF6ScT6Am0VkLYCb2/0lIHkogAvRXLc7EcCFivzfLiKvaj+PtmlnA3hCRI4H8GEAlyzb00IM4kKpiHwewOedvN8B8DtG+g0Ablhh10ZGevE0p870idqlWaTlqWPLtia7yF6ubBS2sPyxlLxG1K/cJKUJTIdmLJ+j1YVlwxqTdN+Dpc698fP6afXPm8itFViXN8S3GpViQjH1DQBe325fAeBWAO9XZd4AYIuIPA4AJLcAWA/gqozd3263rwPwMZIUWfleDUKpzzssdaa3I0TlLRJL073wgmUzUsoWgVqhAY/MUl+iUEGUZ7VjqVVvfCyfLFg+W/7pcfFseKupyB9r8o9CMVG9WYSgiamXfAAcTnJb8jmnR1MvFpGHAaD9f4RR5ijsf7fdUcn+H7ehl9/kvvuy99ZpIwtPAjish18jYxBKfZ7RKfYt687Ym2apYC+tK6/TLEQq1Sqr8z3SiWxY+V5oxlLTli8lKlWXLw09eGq7RJnr/mh4xy9qNx0Pffx0f3IrirTOrMXQLSyWi9pdIrLgZUZ33hXat2bIzrm3i8hDJJ8P4LMA3gHgykydFUUl9QnCO8k9JRfV0WEaq4wum+aXhEc8tZer55XNEbJHzCWTYDR2OvzhteP5HRG0zo8mCssvb0LUk4fV32jCmgeMqzfRnXckHyF5pIg8TPJIAI8axXZiX4gGaO62u7W1/VD7/ymSn0YTc78S++7c20nyAAAvAPD48nuTRw2/TAin3nbNXvWkT0Zrie+FVSxVmluS33DdRgAAE99JREFUAz756zas8IK2bZGkBctvL4ShbUWE5cWWrTZzfS8Jzei6XjjG88/qjzdJWW16q4M05NR9Ttl67VyodBHBnsLPMrEZwMZ2eyMA6zczNwI4jeQh7QXS0wDcSPIAkocDAMkDAbwZzSNPtN23ArhlEvF0oCr1iSMNx1ghCY8kI5VmkU5O3Zbkeb5Z5OuRvKc2+7Th9S2qZ014uq+6ft82vGNRMsnqsbHKl443MNsXRC0IJvaL0osBXEPybABfB3A6ALQPFjxXRN4jIo+T/CCaW6kB4KI27QfRkPuBANYAuAnAJ9oylwH4ZPtok8fR3HY9EVRSnxL6hBQ8gkrtpPuaBHIhFy8EoLe9cIfnnyY/3e9c+CBHolY/PaVt2bHGxDoO2n+NkuPljUVJPd2HFN4xnQdM4sdHIvIYgJON9G0A3pPsXw7gclXmnwC81rH7NNoJYtKo4ZcpIV0iR4otUuSaKCLiiEhK17H8SPdzZJWmW2GEHCJys7a9sbD66YVEvAk2F2LRtqwJNpokPVWfs5/2dR7CLRZEpOhTsRRVqQ8AkTIEytVYLjTgqdQoZOGRokfOEXGmdbWtSEGXtuv5Gals7Vtqw+pXtFqx6nhtWPWiFYM1VvMMQX1MwKiopD5F6BdXA/FJrtOs8p761HZSW96EYpFHLgwSkbIXErHqRSuIEnL0wiq6D96YlExMUR3LRjRxe3Ujm/Oq0DtUFT4aKqkPADkC82LQgB1CsMpHpBaRr+dLul+innXfPEVvqVPLVqTktbLV+9qHSJlb/bR8iyaQtI1oJZPuWysDr715hADYXUl9JFRSHwDSRwp08FR0l9aVsQg5Fx7IkXvUvp40tK1IfUd+p8iRufbDgkfSXvlo0vRUs2fTm5gsG9YKJvIBmL87XTxIfU3GSKikPkDkVJhFiCVqLkcoVrmSEFCuntWm53epXUtJeyGYSElH7eT8jCaWSI17fdK+rQZFHqHG1EdDJfUBQSv26GJZuqTPqVYLJTaikI1ly6tvka+u1yFaCeQmmzQ9GherjVzoKq3n2YrGORcu0vas+qtFoQONSl+sSn0k1FsaB4hTtl6bDbt0+1pJpqRQEoMvJbEodKF9sdrpYBGztu2Rq4dIKVt98yYT3bbln972YE2C3sSYHq+07dUYctkLaZ79UvKpWIpK6gNFFArwiEnXLSEUTdxW22l9jwh1O15fSkg6mpTSMlEsusuP4vmaOK2JKLXj+WCNkU7z+qDH2jo+qxVS+FexFDX8MlBYz2P3Yryeso7gLflzFya1HyX5uQnK89Xz0VLVkeq2ylr+W/b0RKrb0MSu24pWDBa0j/N+26IHAbC7RtVHQlXqM4SIqC316SlZTeB6srDUY2kYJm3LUtwlKrzP6kATqOWbRey639Y46LpeehSO8fqpJ5EKjVKdXsdOoyr1gcO63RHYn6iiMIMmxYgYvVhwjuRK6kQkqmHZ1WSs1bBHtFGYxutLmpdT3R6iMIzle7q/6mLoCgLUC6UjopL6DCE98SMFaNWxbHl1ohi4l2fZsMISnl1dJxdWsdrwJp4uzSLYaGIpGbvoWoXVN8+uN6GuWhBYZA2/jIJK6jMCHWNPFXoKL+6cwiPatLylXksJLErXE1OO+KNQjLbb1bfGxmpDk6xFtl671vh6E5VlxyP21RpDt1CV+miopD5j0M9jB+L7tLs0K7ySCyF4YQ6vrNV+ro0oHBNte6TolbHi6iVKPzdu0UXSiMBTiKyue9BLIBDswZ5puzGTqKQ+B4hCLhYZWSQdxb8jQiqNXUdxeYtkc/U10VqhG6+v2v+SycQidismXhKy0X2vsFHDL6Oh3v0yo+iUXRdmyIViUhLyFLRFgCnJ60+OmLwLmZ1PaRnta0R2nrJO83P1tS9WeKQ0hBNdB9DtWiuTqtL3h0CwWPhXsRRVqc8wLMWr89LtSEl3aZ6yT9M8xe0pb23fCmHk/LNWGlF/vfKe8vbspum6rdwEGdmsyKMS9mioSn2Gcept1+yn8jRxWZ+uXFqnQ1omqqPLalXrqfLUx+jj2Yv8i0IpJTF4L+wS1bPSrdCRNUbpy8grNKR9+kv+sxyQPJTkFpL3tf8PccptbMvcR3Jjm/Z8kncln10k/3Ob9y6S30ry3mPZXQlUpT5HyF2M0+W8C4a5kENq0yOyiGh1XS/eHBGyByve7a0sNKLwk6X0rTZ1H7w+V8QQTCymfj6Am0XkYpLnt/vvTwuQPBTAhQAWWtfuILlZRJ4A8Kqk3B0APpdU/YyInLfSHdCopD4HsB4pUEpkXZncBb9cCKHkwmo0WVg+6DpR37zQUVo+quf1vySMEtmx+l7VeQkEe/D9STS0AcDr2+0rANwKReoA3gBgi4g8DgAktwBYD+CqrgDJtQCOAPDnK+tuHjX8MmfIxco9ePHxNK9Lt9Ssp7q1H144xyJ5a2Lx1HBJ/6LViueDFeLRZaM+ePkVMSZ4ofTFIvIwALT/jzDKHAXgG8n+zjYtxVlolHl6kH+B5N0kryP5kuU6Woqq1OcI3vPYgfgipaVgS+7o8MIUet9T0VY7qU+5OpYv0Z0nVn+1bW3DUty58JDVBll/WNQXPQj7cJLbkv1NIrKp2yF5E4B/btT7QKF960ulZ+gzAbwj2f8TAFeJyDMkz0WzCviZwvaWhUrqc4hTtl6737NiOlh3o+TK6PIW2Wt1baltr57Vri7v9cXzqYSgI3gkbbWv07RvQA259Efz86NC7BKRBdeSyCleHslHSB4pIg+TPBLAo0axndgXogGAo9GEaTobPw7gABG5I2nzsaT8JwBckuvEuFDDL3MKLyZshRO0qvYuFloTgqfyvXKRvdSHNFxiXfiMiDxS+Z09Xc9awWgf0zpWvu5H+r+iHwSYVPhlM4CN7fZGANcbZW4EcBrJQ9q7Y05r0zqchSS+DgDtBNHhLQDuXa6jpahKfU7R3erYKXZLUZaqz65+7iKqRWiamHUdr33dnkXC2j8LFilbcfk+YZvcCietU39YNDqWe7tiIS4GcA3JswF8HcDpAEByAcC5IvIeEXmc5AcB3N7Wuai7aNriDABvUnbfR/ItAHYDeBzAu1awD0tQSX2VIRfOSP8D8QXFKB4dEacXsvAmAY9EvT55Fz89H3JhG8uOFzay+l3RHzKhu1/aMMnJRvo2AO9J9i8HcLlj46VG2gUALhifp+WopD7niGK5+omPuYuQgE92Viw9reOpf23fstvt53xL28qFZUp988ql+1WNrwQEi1If6DUKKqmvYlgEp/ej8IPOy8XctRJPy3XbOSKOLtTq8qUEbsXSPb+8/lWMHxMKv8wdKqmvYngK03qsb3Sh01Pc0cVHr2xfRW+1Y01EOt2rb8XwuzpVkU8Sve5+qUhQSb1iP1jktWXdGXu3LRLM/df1LHWv7ecuYnZlozreHTYdNIF35SuBTxcCYFGqUh8FldQrimAp9ejipBWq0KTrTQzaRi6+rdvJKXsdVtF9jCaRiglBBIsykccEzB0GQeokfxfAzwJ4FsD9AN4tIt9u8y4AcDaAPQDeJyI3tunrAXwEwBoAl4rIxdPwfbUguuCahmsAmxSj0In+ry+0djZzd7RYdq1QS9dG/UHQcNE9JqCiP4by46MtAF4uIq8E8PdobwUieQKan9++DM0DdP5fkmtIrgHwcQBvBHACgLPashVThEfG1n8gvqhpqWjrrpr0rhvPD2sFUS9wDh8ii0WfiqUYhFIXkf+V7G4F8NZ2ewOAq0XkGQBfI7kDwIlt3g4ReQAASF7dlv3yhFyuSGDFn9PHFHjxc0t9R3eo5G5njMI4NUY+a6gXSkfFUJR6il8G8IV223s6WslT0wAAJM8huY3ktm9961sr4G5FDt6th7n7w3MhF2tl4F18rZg9VKU+Giam1KMnpYnI9W2ZD6D5We2numpGeYE9GZnr6fZpbZsAYGFhoa65JwQvXp3+4Amw74X3QirWRc0UVZHPE6Tepz4iJkbq0ZPSAKB9RdSbAZycPJN4J4D0OcRHA/hmu+2lVwwYuV+4Av4tkyU2KuYDAsHiYr37ZRQMIvzS3snyfgBvEZHvJlmbAZxJ8mCSxwFYC+Bv0DxYZy3J40gehOZi6uZJ+12xctC/LvXCLRXzi0m8o3QeMYgLpQA+BuBgAFtIAsBWETlXRLaTvAbNBdDdAN4r0jwQguR5aB5/uQbA5SKyfTquV4wLVYFX7IWgxstHxCBIXUSOD/J+B8DvGOk3ALhhJf2qqKiYFmpMfVQMgtQrKioqUggAqU9pHAmV1CsqKgaI5t1HFf1RSb2iomKAECzK7mk7MZOopF5RUTFQVKU+CiqpV1RUDBP17peRMIj71CsqKiqWQiZynzrJQ0luIXlf+/8Qp9z/JPltkn+q0o8jeVtb/zPt72bQ/rbmMyR3tPnHLsvRHqikXlFRMVAsFn6WhfMB3CwiawHc3O5b+F0A7zDSLwHw4bb+E2geE472/xPt7dofbstNBJXUKyoqBgiZ1AO9NgC4ot2+AsDPmd6I3AzgqTSNzS8lfwbAdUb91O51AE5uy684VlVM/Y477thF8h9GqHo4gF3j9meCmGX/Z9l3YLb9H9X3fzGGtm8Edh9eWPa5JLcl+5vaB/mV4MUi8jAAiMjDJI/o4eNhAL4tsvc2nfRpsXufJCsiu0k+2ZZf8e/CqiJ1EXnRKPVIbhORhXH7MynMsv+z7Dsw2/5P03cRWT8uW9ETYpdr2kiTgrwVxaoi9YqKitWH6AmxJB8heWSr0o8E8GgP07sAvJDkAa1aT58W2z1hdifJAwC8AMDjo/WgH2pMvaKiYjVjM4CN7fZGANeXVmwfEf5F7HtTW1o/tftWALckjxRfUVRSL0NpfG6omGX/Z9l3YLb9n2XfS3ExgFNJ3gfg1HYfJBdIXtoVIvnnAK5Fc8FzJ8k3tFnvB/Dr7as2DwNwWZt+GYDD2vRfh39XzdjBCU0eFRUVFRUTQFXqFRUVFXOESuoVFRUVc4RK6glI/i7Jr5C8m+TnSb4wybug/cnvV5N4Gkiub9N2kJxY3MwCydNJbie5SHJB5Q3ef40h+wYAJC8n+SjJLyVp5s/O2eCjbV/uJvma6Xm+19eXkPwiyXvb782vtukz04cKAyJSP+0HwGkADmi3LwFwSbt9AoC/Q/PKveMA3I/mNXpr2u2XAjioLXPCFP3/VwB+FMCtABaS9JnwX/VlsL4lPr4OwGsAfClJ+08Azm+3z0++Q28C8AU09y+fBOC2Afh/JIDXtNvPB/D37XdlZvpQP/t/qlJPICL/S/b9OmwrmvtOgeYnv1eLyDMi8jUAOwCc2H52iMgDIvIsgKvbslOBiNwrIl81smbCf4Uh+wYAEJE/w/73Hns/O98A4EppsBXN/c1HTsZTGyLysIjc2W4/BeBeNL+EnJk+VOyPSuo+fhmNKgGSn/y26H4O7KUPDbPo/5B9i7DkZ+cAup+dD7o/7VMEXw3gNsxoHyoarLpflEY/GRaR69syHwCwG8CnumpGeYE9Ka7oPaIl/lvVjLSp+N8DU/uZ9QphsP0h+c8AfBbAr4nId4LnTg22DxX7sOpIXYKfDAMAyY0A3gzgZBHpvrDdT347pD8H9tJXBDn/HQzG/x6IfB4yvJ+dD7I/JA9EQ+ifEpHPtckz1YeKpajhlwQk16P5hdhbROS7SdZmAGe2D74/DsBaAH8D4HYAa9sH5R8E4My27NAwi/4P2bcI3s/ONwN4Z3sHyUkAnuxCHNNC+yjYywDcKyK/n2TNTB8qDEz7Su2QPmguIH4DwF3t54+SvA+guRvjqwDemKS/Cc1dA/ejCYFM0/+fR6OmngHwCIAbZ8l/oz+D9a317yoADwP4fjvuZ6P5qfjNAO5r/x/aliWAj7d9uQfJ3UlT9P+n0YRP7k6+82+apT7Uz/6f+piAioqKijlCDb9UVFRUzBEqqVdUVFTMESqpV1RUVMwRKqlXVFRUzBEqqVdUVFTMESqpV1RUVMwRKqlXVFRUzBEqqVcUgeQekneR/BLJa0n+wIh2/irZ/kcj/0dJPkhSSH6T5B8YZY4l+T2Sd43iQ8a/57X9fJbk4eO2X1Gx0qikXlGK74nIq0Tk5QCeBXDuKEZE5F9nijyB5nG2J4jID4vI/+GUu19EXjWKDxFE5Hut3fpMk4qZRCX1ilHw5wCObxVz+taf/0Dyt9v0r5C8on1DznWdsrfUucLhaH66/u1SZ5L2Lm1XEp8ieQrJv2zf3nNiSZlRBqKiYmiopF7RCyQPAPBGNM/+iPCjADaJyCsBfAfAvyuxLyJfBvA5ADtJPkLyhELXjgfwEQCvBPBjAH4RzbNN/gOA/9ijTEXFTKOSekUpntfGsLcB+Dqap/tF+IaI/GW7/d/QkGcWJF8N4N0AfkpEXtySfAm+JiL3iMgigO0AbpbmwUb3ADi2R5mKipnGqnueesXI+J6OYZPcjaXC4LnJtn5SXOmT404BcJU0r0vrg2eS7cVkfxH7vuclZSoqZhpVqVcsB48AOILkYSQPRvNykQ7HkPzJdvssAH9RaPOrAP41yeeN0c+KilWDSuoVI0NEvg/gIjTvtfxTAF9Jsu8FsJHk3QAOBfCHhTY3A/gigNtI3k7Su/uloqLCQH2eesXY0b7E+E/b2x+XY+cHAXxFRF6i0sdiP9P2g2heArFrpdqoqFgJVKVeMWT877BfYbcHwAtW8sdHAA5EE2uvqJgpVKVeUVFRMUeoSr2ioqJijlBJvaKiomKOUEm9oqKiYo5QSb2ioqJijlBJvaKiomKOUEm9oqKiYo5QSb2ioqJijvD/A35UTXirfb+qAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "p = Pupil(samples=123, dia=456.7, wavelength=1.0, z_unit='nm', phase_mask='dodecagon')\n", "p.plot2d()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "`p` is a pupil with a 12-sided aperture backed by a 123x123 array which spans 456.7 mm and is impinged on by light of wavelength 1 micron.\n", "\n", "Pupils have some more advanced parameters. `mask_target` determines if the phase (`p.phase`), wavefunction (`p.fcn`), or both (`both`) will be masked. When embedding prysm in a task that repeatedly creates pupils, e.g. an optimizer for wavefront sensing, applying the mask to the phase is wasted computation and can be avoided.\n", "\n", "If you wish to provide your own data for a pupil model, simply provide the ux, uy, and phase arguments, which are the x and y unit axes of shape (n,) and (m,), and phase is in units of opd_unit and of shape (m,n)." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "x, y, phase = np.linspace(-1,1,128), np.linspace(-1,1,128), np.random.rand(128,128)\n", "p = Pupil(x=x, y=y, phase=phase, z_unit='um')\n", "\n", "# below examples will want something nicer...\n", "p = Pupil()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The wavefront quality can be evaluated by classical metrics:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0.0, 0.0)" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "p.pv, p.rms, # in units of opd_unit" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "or Strehl ratio under the approximation given in [Welford]," ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.0" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "p.strehl # ∈ [0,1]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The pupil may also be plotted. Plotting functions have defaults for all arguments, but may be overriden" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWgAAAEOCAYAAACkSI2SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAf60lEQVR4nO3deZBddZ338fcnHVAcQYSwRCAGhqhPcAkPPTCOC84Aj8FyQOdxAbegjJEHeeaZx9ESl1ILaqoQYRgfx2KMSgnqyKYMKcRRjMu4RemIw46ERQikAgEEFQU6+T5/nNNy+vbt2+f0Pffc373386q61fcsfc63O51v//p7fosiAjMzS8+CfgdgZmbtOUGbmSXKCdrMLFFO0GZmiXKCNjNLlBO0mVminKDNbOhIOk/SfZKuL+zbTdJVkm7NPz4z3y9J/0/SRknXSvrv/Yt8OidoMxtGXwBWtuw7FVgXEcuAdfk2wNHAsvy1Gji3oRjn5ARtZkMnIv4TeLBl97HA+fn784HXFPZfEJn1wK6SFjcTaWcL+x1Ary1atCiWLl3a7zDMhtqGDRu2RsQe3Vxjl12Xx+QTvyt17u8fvesG4A+FXWsiYs0cn7ZXRGwGiIjNkvbM9+8D3F04b1O+b3OpYHpo6BP00qVLmZiY6HcYZkNN0q+6vcbkE7/jec9/f6lzr/nZu/8QEePd3jOnNvuSmANj6BO0mQ0GCbSgXa6szRZJi/PW82Lgvnz/JmC/wnn7Avf2MpCyXIM2szQItMNYqdc8rQVW5e9XAZcX9r8t783x58DDU6WQfnML2swSITRWTwta0leAVwCLJG0CPgqcAVws6UTgLuD1+elXAq8CNgKPAm+vJYgaOEGbWRpqLHFExPGzHDqizbkBvLuWG9fMCdrMklFXC3pYOEGbWRoELPBjsSInaOupiO83ch/p8EbuY72jGmvQw8IJ2rrSVAKey1xxOIEPAMGCHdyCLnKCNrM09L4f9MBxgrbSUmktz0e72N2qTowEY25BFzlBm1ky3IKezgnaZhjklnIVs32dbln3h+Rudq2coEfcqCTjKlwO6RPJDwlbOEGbWTrcD3oaJ+gR5FZzdcXvmVvTveMa9HRO0GaWBtegZ3CCHhFuNdfHreneEHILukVSBR9JKyXdkq+ue2qH814nKSTVtaKCmfWbyPpBl3mNiGRa0JLGgE8DR5GtcHC1pLURcWPLeTsDfwf8tPkoB4dbzM1o/T67Rd0FD/WeIZkEDRwKbIyI2wEkXUi22u6NLeedDpwJvLfZ8NLnpNx/Ln90wyWOVin9upptZd0/knQwsF9EXNHpQpJWS5qQNHH//ffXH6mZ1S9/SFjmNSpSakF3XFlX0gLgHOCEuS6UL7++BmB8fDyJ1Xl7yS3nNE39u7glXY4AuR/0NCkl6LlW1t0ZeD7wPUkAewNrJR0TERONRWlmvSFYUEOJQ9JzgYsKuw4APgLsCrwTmPqz+oMRcWXXN+yhlBL01cAySfsD9wDHAW+aOhgRDwOLprYlfQ9476gmZ7eaB4fr0uVIYqyGh4QRcQuwIr/mGFk+uYxsMdhzIuKsrm/SkGQSdERMSjoF+CYwBpwXETdIOg2YiIi1/Y0wDU7Mg81lj87qaEG3OAK4LSJ+lf/lPVCSSdAA+Z8bV7bs+8gs576iiZjMrCE1lThaHAd8pbB9iqS3ARPAP0TEQ3XfsE6uyA+IiO+79TxE/O85U7ZmrEq9gEVTPbXy1+oZ15N2BI4BLsl3nQv8KVn5YzNwdjNf2fwl1YI2s1EmFpTvxbE1IuYaSXw08POI2AIw9RFA0meBjt11U+AEnTC3sIafHyA+SYIF9fZxPp5CeUPS4ojYnG++Fri+zpv1ghO0maVBMLawnqqrpKeRTRvxrsLuMyWtIBtfcWfLsSQ5QSfILefRNOo9PIRqe0gYEY8Cu7fse2stF2+QE3RCnJgNRjhR96YXx0BzgjazJIjaa9ADz93sEuHWs7UauZ8JVepmNxLcgjazRFTqZjcSnKD7bORaSVbJKNWjVWMvjmHhBN0nTsxWxagkategp3OCNrMkyL04ZvDfE33g1rPN13D/7JR7QDhKSdwtaDNLglvQMzlBN2i4Wz/WlKGtRwvGdhjrdxRJcYJugBOz9cKwJeo6h3oPCydoM0uDSxwz+CFhj7n1bL02TD9jfkg4nVvQZpYEAQsGcN3AXnKC7pFhatVY+oaiHu0SxwxO0GaWBCF28FDvaZygzSwJkkscrZyga+bShvXToJc6XOKYzgnazJLhBD2dE3RN3HK2lAxiSzob6l3XtXQn8BtgGzAZEeOSdgMuApaSLRr7hoh4qJ479oYr8maWCDGmcq+S/jIiVkTEeL59KrAuIpYB6/LtpDlB18CtZ0vVIP1sSrBw4YJSr3k6Fjg/f38+8JpaAu8hJ2gzS8LUQJUyL2CRpInCa3XL5QL4lqQNhWN7RcRmgPzjno19cfPkGrSZJaPCQ8KthdJFOy+JiHsl7QlcJenm7qNrnhN0Fwbpz0cbXQPzwFD1zbMREffmH++TdBlwKLBF0uKI2CxpMXBfLTfrIZc4zCwJAsYWqNSr43WkP5G089R74H8A1wNrgVX5aauAy3v31dTDLeh5cMvZBlHqLWkJFo7V0mbcC7hMWa16IfBvEfEfkq4GLpZ0InAX8Po6btZLTtBmlow6ShwRcTvwojb7HwCO6PoGDXKCNrMkqMYa9LBIqgYtaaWkWyRtlDSjE7mk90i6UdK1ktZJenbTMbq8YYMu5Z/hCt3sRkIyCVrSGPBp4GhgOXC8pOUtp10DjEfEC4FLgTObjdLMekV4RZVWySRosm4wGyPi9oh4HLiQbOTPH0XEdyPi0XxzPbBvwzGaWa+onl4cwySlGvQ+wN2F7U3AYR3OPxH4Rk8jMrPGCNXVi2NopJSg2/1ajLYnSm8BxoG2/YXyoZ2rAZYsWVJLcCnX7cyqSrLLXY2z2Q2LlL4dm4D9Ctv7Ave2niTpSOBDwDER8Vi7C0XEmogYj4jxPfbYoyfBmlm9Ks7FMRJSakFfDSyTtD9wD3Ac8KbiCZIOBj4DrIyI5Idpmlk1I1ReLiWZBB0Rk5JOAb4JjAHnRcQNkk4DJiJiLfAJ4OnAJfkoobsi4pjexuXShg2vlEodUy1oe1IyCRogIq4ErmzZ95HC+yMbD8rMGuMEPV1SCdrMRpcEO7jGMY0T9Cxc2rBRkkqpw/l5OidoM0uCGK0eGmU4QZtZMpygp3OCNrMkSC5xtHKCNrNkLHSGnialkYTJ8ANCG1X9/Nn3SMKZ3II2s0SMVvItY84ELWm3EtfZHhG/riEeMxtRddWgJe0HXADsDWwH1kTEJyV9DHgncH9+6gfzwXHJKtOCvjd/dfrWjQH1TBtnZiOrphb0JPAPEfHzfHXvDZKuyo+dExFn1XGTJpRJ0DdFxMGdTpB0TU3xmNkIq6MFHRGbgc35+99IuolsvvmBU+Yh4YtrOsfMbFYiG+pd5gUskjRReK1ue01pKXAw8NN81yn5mqbnSXpmE19XN+ZsQUfEH+o4ZxC494ZZ/4Z9ZzXo0k3orREx3vl6ejrwVeDvI+IRSecCp5MtBHI6cDbwji5C7rnSvTgkjZNNlP/s/PMERL6Aq5lZ1+rqxSFpB7Lk/OWI+BpARGwpHP8scEUtN+uhKt3svgy8D7iO7MmomVltsn7QNVwnmyz+82TPz/6psH9xXp8GeC1wffd3660qCfr+fNJ8M7MeqK0f9EuAtwLXSfpFvu+DwPGSVpCVOO4E3lXHzXqpSoL+qKTPAeuAP64FOPXng5lZN6R6hnpHxA9p3y24kT7PdY4dqZKg3w48D9iBJ0scAQx8gvbDQbOZmn9YGGg4qqe1jR2pkqBfFBEvqHC+mVkl0lAk6NrGjlSZLGm9pOUVzjczK00EC9hW6pW42saOVGlBvxRYJekOshq0u9mZWa2GoQU9NS6kU9fksmNHqiTolVUDNTMrb2hq0FO67ppcOkFHxK/mcwMzszJEMKYn+h1GnbrumuyRhGaWBg1HiaOg667JHkloZskYshJH112TR34koftAm3UW8f2G+kLHsLWgu+6a7JGEZpYMpd+Fror1kpZHxI3zvYBHEppZEkSwYLha0F13TfZIQjNLRLBAk/0Ook5dd02ukqC7bq6bmc1GDFeJo46uyR5JaGaJGK6HhJKeAvxPYCmFXBsRp5W9hkcSmlkyFgxXN7vLgYeBDRQ6VlThkYRmlowh6we9b0R01bCdczY7ST+v4xwzs06kYIGeKPUaED+W1PN+0P9N0rUdjgt4RjdBmJkxnN3sTuh1N7vnlTinlkevklYCnyRbbeBzEXFGy/GnABcAhwAPAG+MiDvruLeZ9V9dJY65cklDju72AnMm6KZqz5LGgE8DRwGbgKslrW3p1nci8FBEHCjpOODjwBubiM/Meks19eIomUt6ro7cWWVFlV47FNgYEbdHxOPAhcCxLeccC5yfv78UOCJfYt3MhoDYVuo1hzK5pGfqfG5XpZtdr+0D3F3Y3gQcNts5ETEp6WFgd2Br8SRJq4HVAEuWzLkuo5kloVILepGkicL2mohYk78vk0t6qbbndlXmg94JeB3ZSrRjrcerdL6e7RZt9sU8ziH/h1oDMD4+PuO4maVJ20sP9d4aEeOzXabNvibzQG3P7aq0oNeStVR/Tm/mg94E7FfY3pds6fJ252yStJDst9CDPYjFzBoXELWkljK5pGfqfG5XJUHvHhFH1XXjNq4GlknaH7gHOA54U8s5a4FVwE/IWvPfiQi3kM2GQQD1/Hcuk0sGQpUE/Z+SXhAR1/UikLymfArwTbISynkRcYOk04CJfLGAzwNflLSRrOV8XC9iMbM+qaEFPVsu6frCfVAlQR8FnNzLyZIi4krgypZ9Hym8/wPw+rruZ2YpCdheT/W0XS7pB0l7AETE/fP5/CoJ+lXzuYGZWWn11KD7Ku/6+1HgFLKG7AJJk8Cnqnam8GRJZpaGCCjfiyNlfw+8BPiziLgDQNIBwLmS/m9EnFP2QikNVDGzUbd9e7lX2t4GHD+VnAEi4nbgLfmx0lIaqGJmIy3q6sXRbztExNbWnRFxv6QdqlzICdrM0hAMRQ0aeHyex2aYM0FL+g3tR+FM9eLYpcoNzcxmNRwJ+kWSHuHJEY1T+VPAU6tcqMxsdjtXi83MbB4iiG2D/5AwImZMhTFfI1/ikA4n4vv9DsMsWdLhDd2ptqHefSXpqcBJwIHAtWQDZeb1m6dMieOHEfHSQqmjOBGJSxxmVp/heEh4PvAE8AOy8SMHAf9nPhcqU+J4af7RpQ4z660haEEDyyPiBQCSPg/8bL4XqjLd6FOBk8nW2Qqy3w7/mg+/NjPr0nCUOMhaz8Af5wWZ94Wq1KAvAH4DfCrfPh74Ip4bw8zqEAzCIJQyir04AnhaXiKGimXhKgn6uRHxosL2dyX9V4XPNzPrYDiGek/14sgXOTmQbP782+ZTbagy1PsaSX8+tSHpMOBHVW9oZjariHKvhElaKOnjwF3AF4AvAXdLOrPqSMIqCfow4MeS7pR0J9mk+YdLum6O9bfMzOYWMSxzcXyCbK3UAyLikIg4GPhTYFfgrCoXqlLiWFnlwoNkqp+n+0ObPam5/s8Fw/GQ8NXAc4qrPUXEI5L+F3AzFbrcVUnQq9rtrGGxWDOzzHAk6Gi3FF9EbJNUqT5TpcTxu8JrG3A0sLTKzczMZjU8JY4bJc2YVlTSW8ha0KVVmbD/7JabnUW2iKuZWQ0CJnvfi0PSJ4C/JptZ7jbg7RHxa0lLgZuAW/JT10fESfO4xbuBr0l6B7CBrKvdnwE7Aa+tcqFu5uJ4GnBAF59vZvakALY30kPjKuAD+SCSjwMfAN6fH7stIlZ0c/GIuAc4TNJfkQ3zFvCNiFhX9VpVRhJex5PT5o0BewBDVX/2w0KzPj0cnNJA+SIivlXYXA+8rkf3+Q7wnW6uUaUF/erC+0lgy3xnaDIzm2GqBl3OIkkThe01EbFmHnd9B3BRYXt/SdcAjwAfjogfzOOatSkzm11x6rzrgM87MZtZT5RP0FsjYny2g5K+Dezd5tCHIuLy/JwPkTU2v5wf2wwsiYgHJB0C/LukgyLikdLx16xMC7o4dd7RwHLmOXWemdmsor6HhBFxZKfjklaRVQWOmOoSFxGPAY/l7zdIug14DjAx64V6rEyCrm3qPDOzjhp4SChpJdlDwcMj4tHC/j2AB/P+ygcAy4Dbex5QB2USdG1T55mZzaq52ez+BXgKcFWez6a6070cOE3SJNlYj5Mi4sEmAppNmQQ9NXUeZN1FdipOpTeMK6p4GSwbVX3twUGlh4Tzv0vEgbPs/yrw1Z4HUEGZFVVqWwDRzKyTGI6h3rUZ+UVjzSwR1brZjQQnaDNLx+S2fkeQFCdoM0uDW9AzOEHPwsO+bZT09+FggRP0NE7QZpaG4Vk0tjZO0GaWCJc4WlWZsL9nJO0m6SpJt+Yfn9nmnBWSfiLpBknXSnpjM7Edns6ff2Y1S+rne2qod5nXiEgiQQOnAusiYhmwLt9u9Sjwtog4iGx9xH+WtGuDMZpZr22Pcq8RkUqCPpZsUibyj69pPSEifhkRt+bv7wXuI5uT2syGwVQNevCXvKpNKjXovSJiM0BEbJa0Z6eTJR0K7Ei2XE0j3KvDhkkyZY1pXINu1ViC7jQ/a8XrLAa+CKyKWcaFSloNrAZYsmRJxUjNrG+coKdpLEF3mp9V0hZJi/PW82Ky8kW783YBvk620sH6DvdaA6wBGB8fH52CldkgC4ht/u9alEoNei2wKn+/Cri89QRJOwKXARdExCUNxmZmTYiAJ7aXe42IVBL0GcBRkm4Fjsq3kTQu6XP5OW8gm6/1BEm/yF9drb47H2nW7szKS/VnOIDYHqVeoyKJh4QR8QBwRJv9E8Df5u+/BHyp4dDMrCkBuMQxTRIJ2swsS9CjU74owwl6HtzlzgZRqqWNJ41W+aKMVGrQZjbqpkocZV5dkPQxSfcUnmW9qnDsA5I2SrpF0iu7/ZK65RZ0F9yStkGQfss5FxDN9dA4JyLOKu6QtBw4DjgIeBbwbUnPiYi+rSLgFrSZJSL6PdT7WODCiHgsIu4ANgKH9upmZThBm1kaqpU4FkmaKLxWV7zbKfmsmOcVZs/cB7i7cM6mfF/fuMRRA+lwlzksSQNT3shVeEi4NSLGZzs4x9QS5wKnk/1KOB04G3gHoHYhlQ2oF5ygzSwNNfaD7jS1RJGkzwJX5JubgP0Kh/cF7q0loHlygq6JHxhaSgat5QxARCMPCafm/ck3Xwtcn79fC/ybpH8ie0i4DPhZzwPqwAnazNLRzECVM/NpIgK4E3gXQETcIOli4EZgEnh3P3twgBN07dyStn4ayJZzLqJSDbqL+8RbOxz7R+Afex5ESU7QZpaI7gehDBsnaDNLQzBS6w2W4QTdIy51WJMGubRR5An7p3OCNrM0RMATfX0mlxyPJOyxYWnZWLqG5mcsPGF/K7egzSwdLnFM4wTdANejrReGpuU8paFudoPECbpBTtRWh6FLzAV+SDidE7SZJSFitOrLZThB94Fnv7P5GubWMwHbmpuwfyA4QZtZEgKI3k3GP5CcoPvE9WirYqhbzlMiXINu4QTdZ07U1slIJOYC16Cnc4I2szSEe3G0coJOhB8cWquRaz0HbJt0DbrICdrMEuEadCsn6IS4Hm0wei3nP/JIwhmcoBPkRD2aRjYxFzhBT+cEbWZJiIYeEkq6CHhuvrkr8OuIWCFpKXATcEt+bH1EnNTzgDpwgk5YsUXl1vRwcqu5KBoZqBIRb5x6L+ls4OHC4dsiYkXPgyjJCdrM0hCwvcGh3pIEvAH4q8ZuWpET9IBwXXq4uOU8UwRsb7YG/TJgS0TcWti3v6RrgEeAD0fED5oMqJUTtJklo0INepGkicL2mohYM7Uh6dvA3m0+70MRcXn+/njgK4Vjm4ElEfGApEOAf5d0UEQ8Uv4rqJcT9IBxS3qwueXcQbXpRrdGxPjsl4ojO32ypIXA3wCHFD7nMeCx/P0GSbcBzwEm2l6kAUkkaEm7ARcBS4E7gTdExEOznLsL2ZPWyyLilKZiTI0fIA4OJ+XyGhyociRwc0RsmtohaQ/gwYjYJukAYBlwe1MBtZPKorGnAusiYhmwLt+ezemAM5LZsGl20djjmF7eAHg5cK2k/wIuBU6KiAfruNl8JdGCBo4FXpG/Px/4HvD+1pPyutBewH8As/55M2pc9kiTW87VRIO9OCLihDb7vgp8tZEASkqlBb1XRGwGyD/u2XqCpAXA2cD75rqYpNWSJiRN3H///bUHa2a9EGzfXu41KhprQXd6qlryEicDV0bE3Vn3xdnlT3PXAIyPj4/Mv6br0v3nVvP8BeAFVaZrLEF3eqoqaYukxRGxWdJi4L42p70YeJmkk4GnAztK+m1EdKpXj6zWROGE3RtOyDUKJ+hWqdSg1wKrgDPyj5e3nhARb556L+kEYNzJ2Wy4OEFPl0qCPgO4WNKJwF3A6wEkjZM9Sf3bfgY3DFz+qI9bzb0RAZPb+h1FWpJI0BHxAHBEm/0TwIzkHBFfAL7Q88DMrDGuQc+URIK2Zrk1XZ1bzQ1wDXoGJ2gzS4YT9HRO0COuXctw1FvVbi33h0scMzlB2wyzJahhS9xOxIlxiWMGJ2gzS0IETE72O4q0OEFbaYNcDnFreTBEjMzA31KcoM0sCa5Bz+QEbV2Zq2XaVAvbLeQh4Br0DE7Q1lNOnFaFE/R0TtBmlgSXOGZygjazJLgXx0xO0GaWBtegZ3CCNrNkjNBiKaWksuSVmY24qRp0mVc3JL1e0g2StudTGhePfUDSRkm3SHplYf/KfN9GSY3NQ+8WtJmlobkSx/XA3wCfKe6UtJxste+DgGcB35b0nPzwp4GjgE3A1ZLWRsSNvQ7UCdrMkhA085AwIm4CaLO26bHAhRHxGHCHpI3AofmxjRFxe/55F+bnOkF3a8OGDVsl/aoHl14EbO3BdXtlkOIdpFhhsOLtVazP7vYCd/DYN9/MLxeVPP2pkiYK22vyxaK7sQ+wvrC9Kd8HcHfL/sO6vFcpQ5+gI2KPXlxX0kREjM99ZhoGKd5BihUGK96UY42IlXVdS9K3gb3bHPpQRMxY83Tq09qFRftndY08zhz6BG1moycijpzHp20C9its7wvcm7+fbX9PuReHmVlmLXCcpKdI2h9YBvwMuBpYJml/STuSPUhc20RAbkHPX7f1rqYNUryDFCsMVryDFGtPSHot8ClgD+Drkn4REa+MiBskXUz28G8SeHdEbMs/5xTgm8AYcF5E3NBIrJ5/1cwsTS5xmJklygnazCxRTtAlSdpN0lWSbs0/PrPDubtIukfSvzQZY0sMc8YraYWkn+TDXq+V9MaGY+w4fDZ/WHNRfvynkpY2GV9LLHPF+h5JN+bfx3WSuu4X3I2yQ5MlvU5StA55tjQ4QZd3KrAuIpYB6/Lt2ZwO9HuxvjLxPgq8LSIOAlYC/yxp1yaCkzRGNnz2aGA5cHw+1LboROChiDgQOAf4eBOxtSoZ6zXAeES8ELgUOLPZKJ9UMl4k7Qz8HfDTZiO0spygyzsWOD9/fz7wmnYnSToE2Av4VkNxzWbOeCPilxFxa/7+XuA+sifbTTiUfPhsRDwOTA2fLSp+DZcCR6jN+NwGzBlrRHw3Ih7NN9eT9ZXtlzLfW8gaEmcCf2gyOCvPCbq8vSJiM0D+cc/WEyQtAM4G3tdwbO3MGW+RpEOBHYHbGogNsiG0rcNn95ntnIiYBB4Gdm8kulniyLWLtehE4Bs9jaizOeOVdDCwX0Rc0WRgVo37QRd0Gh5a8hInA1dGxN1NNPRqiHfqOouBLwKrIqKpKdNnG1Zb9ZwmlI5D0luAcaCfizF2jDdvSJwDnNBUQDY/TtAFnYaHStoiaXFEbM4T2n1tTnsx8DJJJwNPB3aU9NuI6Mn8sTXEi6RdgK8DH46I9e3O6ZFOw2pbz9kkaSHwDODBZsJrG8eUtkN9JR1J9svx8HxGtH6ZK96dgecD38sbEnsDayUdExHFCYisz1ziKG8tsCp/vwqYMeFKRLw5IpZExFLgvcAFvUrOJcwZbz5s9TKyOC9pMDYoN3y2+DW8DvhO9Gdk1Zyx5iWDzwDHRETbX4YN6hhvRDwcEYsiYmn+s7qeLG4n58Q4QZd3BnCUpFvJJu4+A0DSuKTP9TWy9srE+wbg5cAJkn6Rv1Y0EVxeU54aPnsTcHE+1PY0Scfkp30e2D2fl/c9dO450+9YP0H2V9Ml+fexkbka2ikZrw0AD/U2M0uUW9BmZolygjYzS5QTtJlZopygzcwS5QRtZpYoJ2gzs0Q5QZuZJcoJ2mYlaVs+6OJ6SZdIeto8r/Pjwvvftjn+XEl35vMS3yvpU23OWSrp95J+MZ8Y5ohvp/zrfFzSorqvbzZfTtDWye8jYkVEPB94HDhpPheJiL+Y45SHyObYWB4Rz4qI/z3LebdFRO0jHSPi9/l1Z8yvYdZPTtBW1g+AA/OW7PVTOyW9V9LH8v03Szo/X1Xk0qkWd7tWc4tFZLOt/bpsMIX7fS5v4X9Z0pGSfqRsFZlDy5wzn2+EWVOcoG1O+UxyRwPXzXHqc4E1+aoij5BNvzqniLgR+BrZrHVb2q3+MYsDgU8CLwSeB7wJeCnZRFUfrHCOWZKcoK2TnfKa7wRwF9nkRZ3cHRE/yt9/iSwRzimfCe7twEsiYq88YZdxR0Rcl89hfQPZEl9B9otkaYVzzJLk+aCtk9+31nwlTTL9F/tTC+9bZ94qOxPXkcBX5jEfdXHO5e2F7e08+bNd5hyzJLkFbVVtAfaUtLukpwCvLhxbIunF+fvjgR+WvOYtwF9I2qnGOM0GnhO0VRIRTwCnka0EfQVwc+HwTcAqSdcCuwHnlrzmWuC7wE8lXS1ptl4cZiPF80FbLSQtBa7Iu+R1c50/AW6OiP1a9tdy/TnufScwHhFbe3UPsyrcgrbUvJOZS18BbAOe0cuBKsAOZLVpsyS4BW1mlii3oM3MEuUEbWaWKCdoM7NEOUGbmSXKCdrMLFFO0GZmiXKCNjNL1P8HUGX5gkCJyhkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "p.plot2d(cmap='RdYlBu', clim=100, interpolation='sinc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "`cmap` and `interpolation` are passed directly to matplotlib. `clim` is made symmetric if only a single value is given. A figure and axis may also be provided if you would like to control these, for e.g. making a figure with multiple axes for different stages of the model" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcwAAAFmCAYAAADgXKoOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9e7D2WVUm9iybi5ggAt2OpJuexgIzAa/FJ00qxlhCYzOltlNyaUVpJ1hdOkOZya2CZZQEnRSoiTMaywkKCnhBwTF2IhbpBkkmEyHdCAOCGlqckQ8otWnCMBrA5tv547y/7t3re55n7d973nPOe87ZT9Wp87vsvfZaa+/fetba7y1aa5iYmJiYmJjw+JyTVmBiYmJiYuI0YBLmxMTExMTEACZhTkxMTExMDGAS5sTExMTExAAmYU5MTExMTAxgEubExMTExMQAJmFOTExMTGyFiLgxIv4oIu6OiJeQ+18bEb8XEfdFxHPSvVsi4gObv1u660+NiPduZP5kRMRx2DKCSZgTExMTE6sREVcA+GkAzwbwZADfFhFPTs3+FMB3Afjl1PcxAF4K4HoATwPw0oh49Ob2zwC4FcCTNn83HpEJqzEJc2JiYmJiGzwNwN2ttQ+21j4D4PUAbuobtNb+ZWvtPQAupb7fAOD21tq9rbWPA7gdwI0R8TgAn99a+9128K06rwXwLUduySAectIKTExMTEwcDb7hxi9vH7vn32zV953v/JP3AfhUd+mVrbVXdudXA/hQd34RBxXjCFjfqzd/F8n1vcAkzImJiYkzio/d82/wjrt+eKu+D4nv+FRr7YJpwl5bHP2uVdX3MDKPHJMwJyYmJs4oGhouXfrsUYm/CODx3fk1AD6you/Xpb5v21y/ZkuZR44zT5hXXnllu+66605ajYmJiYlVeOc733lPa+2qw0lpaO2+3Sh0Oe4E8KSIeAKADwO4GcC3D/Z9M4D/rnujz7MAfH9r7d6I+GREPB3AOwC8EMBP7VjvrXHmCfO6667DXXfdddJqTExMTKxCRPyrk9bBobV2X0S8GAfkdwWAV7fW3hcRLwNwV2vttoj4agC/AeDRAL4pIv7b1tpTNsT4wzggXQB4WWvt3s3x9wL4BQCPAPDbm7+9wJknzImJiYlziwa0dmRbsmitvQnAm9K1H+qO78SDt1j7dq8G8Gpy/S4AX7pbTXeDSZgTExMTZxQNDZeObkv23GES5sTExMSZxZG+hnnuMAlzYmJi4sxiEuYuMQlzYmJi4qyiNbRLkzB3hfnVeBMTExMTEwOYFebExMTEWcbckt0ZJmFOTExMnFnM1zB3iUmYExMTE2cWDbj01yetxJnBJMyJc4k7nv5ctBaIeOB7nZfz/L+/D+BB13rk+0pOxjPf/oZD2zMxwdDarDB3iUmYE+cSC7n1ZOaIMKLdT36KOHuiVOO5NhMTu0cD5rtkd4ZJmBNnAnc8/bn3H7uK7oF7QQnQEakiyr4ta8Oq1P5ernazPnm8WZFOTJwMJmFOnBkoomFbrT0yuWXS7CtMhixTVZFVRZorWWfjxMQQ2qwwd4lJmBOnBrdf/zwAvPLqyS3fU69DutcsmSz3GmUei1Wa1WuZPVS/XJHmcWb1OXEZ5muYO8MkzIlTg0xUC1hFyCpMR7S9/H68fsuW6cPIMt9nW7uu0lU6KcKemFAINMSsMHeGSZgTpx5qqzOfu61PVmWOvMvVkS+Txe4xmey+q5oncU5QzC3ZnWIS5sReYdluBMaIYwF7444jq75PrvTya5bbbrOq1zJVn15+VU2r11aXbWvmi7ldex4xCXOXmN8lO7FX6Mklb13226P5OPfJW7Pu9UZVufUy1UdLsi7MDvZGHtVfbeuO+oYlBrnvxMTEdpgV5sSxo3/zDpCrKE1sPVhf9q7WTGbVG4Pca5yOdNQ7WyuCVjYynUf6uXvLR29Y8nDDO35Nypk4zWiI+aafnWES5sSJwBFebtPfcx/VYK8bqrHcOdsaVWP0fXsZihyZTso/lb5sPGYLs2k0EZg45WgALn32pLU4M5iEOXGk6L9QYEGE/3C+et3SvV7IiCNXk4wo2Bt7GDGrd66y1xdVv9y+P1fkrggxy2bjMLJkfSLag147XtrO1zzPAua7ZHeJSZgTRwoVrDNRjJKhq+7YmPl1vur1SjU2s6WXU22ZsopPVXjVm4SyP1S7Si9G8Or6xGlFmxXmDjEJc+JI4QhKbbPm+xXJVdcYSbMqi5271znV66yO/Ef0H3ltdUReRbz5WlXFT5xCzI+V7BSTMCd2hvwNNFUAdsE7V4Vse5H1W+658SsS6sdXxOn0UNugatuXtVHVrEsulJ5uW1cRf2uB269/3mV95jbtxHnGJMyJnUFt8VVtlnYs6Gc5mQT6ayPVam43upXK2rrXL9VWLNPHbSkz+5ld6lr12ifTdznvfTRxehFzS3Zn2CvCjIgbAfxjAFcA+LnW2stFu+cAeAOAr26t3XWMKk5ssFQfD3698IH7KogzUlPEN0JUCq6yy/eZLu5NREyHisiUHup+hqqw3dbx0o/pPpI4MF+xXYRZde4x2nwNc5fYG8KMiCsA/DSAGwBcBHBnRNzWWnt/avdIAN8H4B3Hr+XEArZV6LZC1ety/f+KIFm73MZVtCNyVXulo3o901Wt+do29ldb3uq1WQW2VVz1G0lcJk4es8LcHfbpm36eBuDu1toHW2ufAfB6ADeRdj8M4EcBfOo4lZu4HEtA7gOzql4y1JtrGNlkMnpwVet/dov9d3qw9lX15cZwtlSyc3vlLya/kun8XSUFSr+JfcSmwtzmb+Iy7E2FCeBqAB/qzi8CuL5vEBFfBeDxrbX/NSL+CyUoIm4FcCsAXHvttUeg6vlD/paYfNxjhKgy3NZhPxYb31Vo1diOSBipsyqv11/ZprZsWXVabc0qW6sqkrUZ6eN8ybZp57cG7Q+itVlh7hD7VGGyiHP/ExoRnwPgJwD855Wg1torW2sXWmsXrrrqqh2qeH6RiWLtG0L69uw1zExMjHTZFi6rEPsxna5uSzjbmqvavn1/vRqzb6/s7/Vj1V+eA2Z3Hntp7wi+kufkM7smJs4a9qnCvAjg8d35NQA+0p0/EsCXAnhbHLy75IsA3BYR3zzf+HM06L/zdXlDDwu61Wtw/T0VsEcrnVxxsorMVX/qNccMVm26rcksq9q+dCTNkoQsq+/nqtVKL9emv+6q9ty+342Y1eYeYFaYO8M+VZh3AnhSRDwhIh4G4GYAty03W2ufaK1d2Vq7rrV2HYC3A5hkeYRg1QwLsP2x2n7M56PEkmWrsVlV54iN6eeqsFwBVkTYn1fVKZOZr2UyZDY4e9XrwA7V65ZsR2C+prln2GzJbvM3goi4MSL+KCLujoiXkPsPj4hf3dx/R0Rct7n+goh4d/d3KSK+cnPvbRuZy70v3KFHDoW9qTBba/dFxIsBvBkHHyt5dWvtfRHxMgB3tdZu8xImDov+g+oHxOFJh70elo/Z62WsQh15HY5VVfkea8/GUH1UXyVvpOrMMkZ9qrabXfKSK2Gns7quKvNlHOVHJj9/l/D8CMoJ4IgqzMFPNrwIwMdba0+MiJsBvALA81trvwTglzZyvgzAb7bW3t31e8E+FkN7Q5gA0Fp7E4A3pWs/JNp+3XHodJ4wWnkwKAKrCJCR5ppxRqvITFT5NcBMfG4L0slX1eaIb9h2M2vXj8fmjFWlFZGrXYCqMmWVPrNl4qTQEJcuHZXw+z/ZAAARsXyyoSfMmwD8N5vjNwL4HyMiWmv9ovg2AL9yVEruEvu0JTuxR6gCbL7Ggqm63stRVVYmHlZVLv3dlq3bqsxbikpPRqy5vdp2dRWbS1AU2bLzrCfbSmb3ss6uUsy2q50FV/VOnAAaDvOxkisj4q7u79YknX2y4WrVprV2H4BPAHhsavN8XE6YP7/Zjv3BiNibxbNXFebE8aP/xp5lWbrXoxQ5jWylVluoKlArfVS7vL04WjWyik4Rg7NN9WPVGKviHEZfD1Vjj5CvmwOGbAur5vMWbWvzDUGnAPe01i6Y+/aTDSNtIuJ6AH/VWvv97v4LWmsf3nxJza8D+E4Arx3U+UgxK8xzDlZBLMfL/XyNvR7Xy+hl931U5dhXLVkHV50qubmd0zNvparx3bZo/mN+yG37/orMlG1VxVjZ0stRPlPXmB8d8Tv9Jo4D7Si/uKD6ZMOD2kTEQwA8CsC93f2bkarL1tqHN/8/CeCXcbD1uxeYFeY5Rf4iAlUR5YCugjer7nK7RTaruHJFyHRSlWomCLcFq7YZe7DAz+zL4zDiyInHGrJQVR3zhdJZkZgjS5Y4VTqqPsxPy/9ld2O+EehoEe3IXsO8/5MNAD6MA/L79tTmNgC3APhdAM8B8Nbl9cvNZ+ufC+Br79f1gFS/oLV2T0Q8FMA3ArjjqAxYi0mY5xSZlFRF6Qii79fDbbey60sfto3qZCoCVPeYHLedrAjHbW9m4qhISxExI3FVWVZtnZ1VQqTg1kROjtwW7aw2jxhH+OXrg59seBWA10XE3TioLG/uRHwtgIvLm4Y2eDiAN2/I8gockOXPHokBW2AS5jkA+11D4PKPH7jtNdZmuQfU74ZkQbmq7nJ718cRAmurdGVy8rZttkmRl9tqVeTKdB8hYEd2rgpX+jB7KpIeucf0yh89AebHT3aKo3uXLKpPNrTWPoWDKpL1fRuAp6drfwngqTtXdEeYhHlOUG1XKtJU22m5YlDXVPXHxu77rqlSe9nOXlWN5j6jZDAqO19XSQOrOrOM0XFUnzy/I+RWVc6qalTjOx2rxGtiJVo7UsI8b5hv+jlnYITAArbaomVy3BalGmcBIwB1XhFtHktVUY5kGFw1moleJRhZr6VvPu+J1JFXT3x5Thk59j6okhGVHCkb2H2WNE0ynDjtmBXmGcbyKxJA/W7Y3KbaKlRkWG01VgSSoSoQRayqgnE2KrmjVaojVKZ731f5wM1XnjtHaOw+87cieKVL7suSk9E1wHSfbwjaHeavlewOkzDPMHIQAx6c/eegxkgp93dbnxWRVKQ2Qj79eSacTP5525MRYB53hDSdHzJGiF3NU9VG+YL1Z+3XJAMuCapsqnzS6612Aia2xdyS3SUmYZ4xPPgNFOsqPYYcCNXWYz+GqqxUgHYBuKoQGcEw0lAEzs5ZtZWJMicaCqqCroib2Zt1rao2Ze8oEbmEgm3hqrGVX6vKGXjwG9ZmtbkFGiZh7hCTMM8YcravtjsrwlsTENm1UdJjUP3U9nG+14/fy+vbrdkuZOOyJMBVyCP3nZ9Y8qH0YxW1S1icbqM7B0zvKjlSSQ6TObEtZoW5S0zCPGNgW3Es4Fek1fdT23tqK66qzNgxq34cSSj78rjsfhWU3Zas2pplldcomB97XZ2PnF/yGEzvkW1UlTyxuVbrhLXv+7H2apyJFWiYv4e5Q0zCPCPov7kH8FuI/Xm+xrZVVbXIjnP1oLY3HYEq0lc6uKpxhDiZDdkPbhy2LVtVbJmUq4TCzRfzRbajH4dVm0w3tk3KyM3NF+s7UmkzvZc1PrdmJ04KkzDPENR2bH/MAvzaLVJFWu5/7sfIr9dlDakqvVyl6Pw2YjPTlVXb+VglEEubqqpn9rq5VMmFWhtKR6aTWgO5rdNpGT/LUPJmpbkOcbQ/73XuMAnzlKN/U4QiKoURMnGBeyRg9+3yudMvkw1rO6LbiF6K6Jke1biKCFWl63R0eimbnT7bbImOEGDWl8mudhEySTtbZqW5EpMwd4ZJmGcALJCxwOO295TcvH239HdVZn/ej8cqF1eBMl2Y3i64Vluaub2qyJSdauvREWdViTObRirvrEO/Bcv85SpyVY0qqASlSo6qJICN08ufKNAwCXOHmIR5CvHgqvKB6+x1pIM2nmhYds+CvgruqqJlgZoRa98+t8v2sO06Jltt1ypdmS9GyNxViKqSZFXsKMGMVvXZdnad2eWIvCfhLF8lOlVypBKOrA/r3/+W6/xtTYU2CXOHmIR5CuECHmvHKp/RykuRVr7G2uT/rHrJ+irC68dxwZzJYOMw3RnBqKowj+nGUYlJPs5jqmpqzZz311VVripLZVOVcDgyZDpuW9lW/p+Y2DXmd8meQvRZPruX/ytS7AMX2+YbqWJGtvXWbLupyjK3YbpW4y12qkDNxsg6KAJlc7JcG0kw+nFVBc18VPk1kzOb/5GKla2rrPNy3ZFlf858rypaNZ9r1ta5RANwqW33N3EZZoV5irC82WHZhmXBeDQIseOlLYOq+PJ4rtpSBOWqTdZXtWN6MIJjbVmgd+OosZiNKhlRc8LgyK1KKBwBuoRBzVe1BpQ8dS9fV9Wo8lNEu/97k+fWLMHckt0ZJmGeUoxsReUgk7dGc1vAb/GxQKkqCtZ2JFCy6k0F6CzDJRCqzcj2Yta/ItdeRtYh27ucMx84MAJyVWn2rdIjt2fnrL/zYZ5PZUfWkfnZ+bvy2flEm4S5Q0zCPAXo3+SjAqqqslxVl+WwYKbIwpGaIy1VPbBzpkuvQ6WLq8Rc5bkGjjwdMSqbqqpwJPGoruXrjKCy3u6c6enG7sfMdo0kJy5Zm5VmwrIlO7ETTMI8RXCv1fTBRGX3TI6qQllF0o+zHK8JroygM5mrStVVG6rCqQhEVV0sAWBbqJlo2NjZ11kHhaqiY7ZkO1l13N/rz51s55e+j7PP+cwRJUO1lmelmdBmhbkrTMLcU+SPjuSA6KpBtf3V91NV6NLGkSaTsVx35MH6KlLMOrsqoz931aYisb6fC9RK9z5J6e/355mQ2T0nn8lw1boiXEZUrAJWYyvyVnPG7qmkhcGtZ5fMLZVmxPxNzYndYRLmnmIkgAD6zSNMFgucblvNVaeMLFVfVbExO1UAVAFbkTUjGSc3y1LVpPIdszPrOVIp5vGzLEa8LLlRemSblT7Zv85uJTfrrYhyJJmp7M/6KyI/f2hzS3aHmIS5p1DEsBy7zJz1USTbt2MyGTEpkmO6qLFYe5cksPGUPo6Qc58RsnIk0Z87WRWJsfFcssGqQXacx1FropddVdkjpOxsYm3Z2mPVshrTrelzjYZJmDvEJMw9Q/7oCOCDuApe+V5Fso4I83ijcBWsC+pu+9LpXOmm7FcVpdPFbUEyfZzP2Ty6+VS2ZYwQoIIiW5eIjMqqSL0frx/Xje/aLC9vnNut2UmYO8MkzD2D2mZb7rFtKVWBZJm9DNWebcv1153srKNqx4iw2oZ0VVIeP7fJ+ld2K9kjhMzsG7VjxH5WObH1MUK8VYK0ZhzXjo3PbFS2s75V9Z71PK9oDfM9PzvEJMw9wFJVAgeVZfWAM9LMx33bBSPEmPvka5k4XYXkAneWqSoZ1Te3zeMqP1TVurNZ9R8hqnzNJRD5fsYosWV9mZ/UsarGXcLmEqER0lK2MBuqpIQlbP1zdq6qzVlh7gzzq/H2DG7LjWXcC1SG3f9nQUfJ6UlnJPtnOih9HZkpPdYQSh/Ac9AfIen+WFWFVVXf6850ZokFS0acjSMkksfudWK+URU6s4G1qxKmkWTI2c7GUwlB1peto4mJNZiEuYdwVYAK2Lm/q3bcOQuYbruL9VcBURELk52DuguUua+Tm3XvrytiVKSvSFbpqMjCzUcmQqaLSgj6djmJyDLU+Gw+q/4MIySpbKwqVEeK/fqpkpAziQbg0pZ/E5dhbsnuGaqA11/Lx31AyEFNVVZZVk9qrOJQQT7boIjRVVaOAEcqE9aul6P0yj5nCUdVKbIKeOQek8n0HSFz5geXdDFZbs6yTix5YOusst/Nf+V7Nk9urTN9zjwm+e0MkzD3CGqrbIGq+pZAwYKUIk4lt5fX33NBjgVMRYQjbUZ0GbFH2TDqB0foa7ZPezuYLkzHioC2SVyqY+a3ith7HdRayeuzIvPR9cjsd/4dnaszh3Nm7lFiEuYJov82nwUuaFfk5wJLlq3kVeSw3GMBMfdXwSvLctdVgGW+yjKYb1RSoUi0H8dV8mwb0FW6Izax+WB2KNvdPDHfqDEcWS3t1JjMLqYru6cSpCrZYONlecsbgM78m38a0C6ds4r6CDFfwzxB9A8w2+bq2/TX+raZbKoMeqSi69v2MpUurqLIumUdsi2qoukDoSIj54NMVFnXPB67l/2S54r5LdupkogsL/tegc0DqxbzmI6o2XwqO/s2LsnK/fJ8Ovt6m/IzwOZoW73OLI7wNcyIuDEi/igi7o6Il5D7D4+IX93cf0dEXLe5fl1E/H8R8e7N3z/p+jw1It676fOTEbE3kzUrzBPA7dc/DwCwfITEZcMLGAkt1/v/rCId2YJSRMXaMVLPwb0i0Mo+V/GxMVnQdNWR8qMae4TwWWBWc6AIyVWOrA9DNffsHuvDKjeley9nuc58rdYIk8HWoNK5IkDl7+VZnL9ush4RcQWAnwZwA4CLAO6MiNtaa+/vmr0IwMdba0+MiJsBvALA8zf3/ri19pVE9M8AuBXA2wG8CcCNAH77iMxYhVlhngAqAmNByFUDfZsqoLAKsiJite1WVXtZL7ddxuSO6KiIiVVvrC2rcJncXo7zU77O+lQVEWszMr5bF6qiVeP0/5UPFBkznVTCkm1XayXr6gherS+VRI0klKcal2K7vxpPA3B3a+2DrbXPAHg9gJtSm5sAvGZz/EYAz3AVY0Q8DsDnt9Z+t7XWALwWwLesNfmoMAnzBOCqp5EgpNqyoO3IjBGUkpOvjVaCWR4jdEcuSscMRgjOZ7kqZD5gdqjqbKRPlr2mcutlrU24ej1HyEWRXC/fEW1/rGxU42ffsPXMSHeE9CpCPZNoB69hbvMH4MqIuKv7uzVJvxrAh7rzi5trtE1r7T4AnwDw2M29J0TEuyLif4+I/7Brf7GQeWKYW7LHBPZtPpksVHX4QL/LA5gL4KpvrqwUQWZdXAWgdB8JzFkms1HpWyUJqopkOjH9qnnKYzD7XD83ZvaD80seg/khgyVJrp3yXZ4LNZfZJmd/1ovpnOdT+Vz5spe/PJ+tnbUfnx6uFhnuaa1d8MIvQ15Aqs1HAVzbWvtYRDwVwP8cEU8ZlHlimIR5AnAZ7UhwccGN9V2OVTCsxmT9RoKeI2Smr5Jb6T6aLKjKSNmmxh8Zk/UfJTDVdrR/TnCyDo78ne4sUcnkWBHzyLp1vnXrrgJb607XM4Ojq6AvAnh8d34NgI+INhcj4iEAHgXg3s1266cBoLX2zoj4YwBfsml/TSHzxDC3ZI8R/cM6ErxZQMr33dZYrgTZVtZoVq62SnObKuhXAYnZ21eHTBemuyN5199B+cZVMKrvci1Xv72dowmM6sf65EqbraF8LfvN+Zbp0yPbWo3H9FdkqdY4G4PZM0K6pw6H25KtcCeAJ0XEEyLiYQBuBnBbanMbgFs2x88B8NbWWouIqzZvGkJEfDGAJwH4YGvtowA+GRFP37zW+UIAv3l4R+wGs8I8YbCM223xVRVYJqosR8nM/ddk4KyvskvpoWQzG1S1VLV3FZdLDFzwV37JFVhPUo7Qsm5uzP68qooVlC6qwq2OHRix9WOoBFGtFbbO8jGriJXcM1ldHiFaa/dFxIsBvBnAFQBe3Vp7X0S8DMBdrbXbALwKwOsi4m4A9+KAVAHgawG8LCLuA/BZAN/TWrt3c+97AfwCgEfg4N2xe/EOWWAS5pGj/whJhnq4M1imra4puWqbi5EOy+BdVVEFQBZ8mQxlN5Pv9GZJgBpfJSF5fIYRu/IxO88YTV6ynlln1tfNpRvHyVJryxEfO8/jKEKuiHqkUswJU78OztzvZ146uo3E1tqbcPDRj/7aD3XHnwLwXNLv1wH8upB5F4Av3a2mu8EkzGNARRI58DCiUdVUlscC9UjFx+T17UYqRlX5uWSA6V3JUJWismWEIJj+ikwZKVWV6Rp/OAJnslh1phIM5Wunh1p7zK8jPnHHzn51LdvYt1NVJtNT+e5Uo8Vh3vQzkTAJ8xigqo2qggH0Gz7yFqILVL2cvq0ijRF7qqqSjZMJxAXUkfGzXUy/qkqrKivma3esiFGRsNOxqtBdezXfzCY2T9kHzM6RJC6DyXcE6cDWW5ajbGH3nT2nGWfJlpPGJMwjwgPfE6tfe+zB7qsKiJ331/q+jBjzfSUr65Vlr8322bhMB0ZGrlLMUIHU+c8RVtZd6czarvHRaPXjKs8RO/q2+dzNhVuDjKQU+TjSZ2B+V3ap/ipRdc/Gmfm+2SPckj1vmIR5jHBVoMt82f0sb+kzQg79GFUbpoPKxEeTAUX0juD6c0bcqpIbCa4jFTMjxlHyUoTmdFb9WCLD/vf31by7tcJ8pGS4xC3rxOzaxt9VwqX+s3nINjm7Txtaw/zy9R1iph5HBEUCLJtV5MOy9urcBaO+LQvcarylvSO+ingdWfUBzQVrVdGw8Zk/sl3MBlWJMVJmujFbMwFmO1mQHiEhZbMaM/8fIcZsb3+u7FVjuoRArZ/eNkW2imBVIlXpPJIETpxPzArziDBS3bBrmaiqgJ3bsfHZOWuvdHMViiJHRmwV2fXH2VZHdOy6ImE1tqs4qoDMSIDpuTb4jqyh3o5qbFbNuXl3iYFbgyzxquZe+UnNo6r8nHzXlyWUTN7pw3zTzy4xCXOHeOAjJA3A5UFaBZXl3tJWZf8sC3cZOAMLmnmMipwUqVXHDIxcWYLAxnXy1HUWmBUp9/YzfV2/KvFh/VkSoogmj8OqNuaD3EeRqiNqtuaUT/PYefz+nOnOKslMoioxy75h/hh5du54+nPvv3cavzbPPS8T6zAJc4dwgUQF3eWeytxZ++UaOx7VkfVVRKmCP+tfVV9VJVIRr7LXEQc7ZmMqWyo5inj766pCZD5WSVWlMztXJMRsdPcY2Vb6snU9OqaazwzVnvlixBamf+53qtAw3/SzQ+yVJwd+jPQ/i4j3R8R7IuItEfE3T0JPBVat9ccjhFc95P041daS0ovdc9WCC6KKfJnOzuYc6BQJZd17nUf9y4Krq7SzXiywuyRJ+SXbzHzJ+ipSXI6zfSNrKsty5MDmwyU/IyTW33PzqWRVxNrDVZ7uWTytldoRfjXeucPeEGb3Y6TPBvBkAN8WEU9Ozd4F4EJr7ctx8NtqP3q8WtbIQViRG7vXZ+O5bW6X0bdhMplXTKcAACAASURBVFh1xchxdFwWIBkU+bJxqv7ZNhcgmd+zza6vmi9VKfXXVLWs5pXpkteHS5Iy2FyPjNn3Vfoogl8DRfDZv6xP9UyxREaRe76nnpHTjQd/R/Gav4nLsTeEiYEfI22t/U5r7a82p2/Hg7/Vfq+QH97+YcwPpdr6coGLnWdZKsAr8mKBa2nnAm22jenHAl0ep2rvqjEXdPtrmRBye0Z2zFY3pksm1Jzm+1VVmO1friniZvPC/O3WEPOR8uUo+ai2+Zi1UXZlP7Bnj81zllklVhPnD/v0Gib7MdLrTfsXQXwp7+aHTm8FgGuvvXZX+kksbwpg2e2D9dKvffXnOQDm9hVYfxXM2ZZUJmAXgFyA7tsqnfI4yle9n5w9I35iVaDyHQueVSLh5o2RDAvMbJ2w9mwclggxGS7pUPqofspfTmfnJ0WwGWouRwiPJZLVM3fqvmd2voa5U+yTJ9mqppEvIr4DwAUAP8but9Ze2Vq70Fq7cNVVV+1QRQ6W+argtdxzQVsRRRVQFSGrcyZHkWjuk8lQ2c2CM/MXG0/Zpfqp+1UiU+mkAnxvcz+W8jNrw9qyNdLLYLb3FWLVlunB/MVsc1Vilq+wtmLLa1QRoVr7ysf5nvJrvn/aKs75GubusE8V5siPkSIingngBwD8R621Tx+TbhYjD6cLYEqG6seuKaKq2lfZd2WHC/zsetapuj9yj1UTyi9KL9XWEZSbz6yXI6aqylH+cPZUpM3sUP2Zzkp2tf7ZebU+1by4cdg6deuGyRyt4PcZDXNLeZfYpwqz/DHSiPgqAP8TgG9urf35CehYwi1OFthVexfAXBbMxlHZ9HJNbcUt91Q168Zn4yh72Dae0qvym9NH6ZLbumpT3asqYKdvVbWodaOqTHZcjePWAZOfz9l8Zv0ZOY8Gc2ZX7l8lHsz+kbW+Rs+9Q4uDLdlt/iYuw95UmIM/RvpjAP5tAG84+DFu/Glr7ZtPTGkBlm1XhNFfW1O1VA9yVflVAUKN4bavqupVVY+qOlT6jfh5NNA5n247nxmVH12Fm9sqG1jy0Y+tKjk23sg6G630mF35uErynN3Mtt52pWd/f+R5WtqfJszt1d1hbwgTAAZ+jPSZx66UwPJLBsDYA6SIaTkfCdK5r6pi1JZXfzwS1Ee2sXqdRuzK15leiqyd7uxeDoIukOb+LolxSUYes6pqVbBnVZnyG9Ol10nZwdoxXXJ/tTbVHORxVbWc/ef0cP2ZXUpnpmPlU+DBz/+peQPQxKEx6+4twbJ+FRAWjGTLSzsWVPL4jGBYVaACggoeOWjlMbP8RQe2XeYqL1ZJqUou+8ORtgrm+Vj5Ivugb8t8MFLls3nM/lNEwuxT/mH3mZ+ZL5ivmC2sXdZ1JClTtmR7GJH1a4UlbC6pGdU/t8nyTkulmZ/P0b+Jy7FXFeZpRPVA9+369lmGyqpZ5cF06MdUQdAF+qoSYdVPHjMTg9J5pCJwMpRvnJ5qHNaX6ZcTFCZXJQBuPnNf5scMJVONm21l7Zn/Rv3O5tzJHpkvN34+dv7M99nzN0IOzmd7jRbz9cgdYhLmDqBIkAWVTI4LXIDM/ftrTNYoCam2jMCZbipQ93KzD3JwrWQym5lcF8yzTKZPrlRGbXUBX5G1uu/8nG1xJJPlMqLv/alInOmu5pvZ4UissnM5r8jVrSNmZx6HJYq9jOq5Pg2Yr2HuDpMwt0R+cNmDzQJGlWUv5646UrooIlZZ/6hcJZu1UTJcQHYE6exQ8lUwHwl+zjdrAv5IwGYYrbCUDWr+15LZiH7b2JP1UtfcveX+Wp+OJpZqnNHncZ/QcHq2jk8DJmGuBHuzjwpYbiuIVZz5fJtMVsmqyEOhCtgj/VnyUBFVVUVVskYDdVXB5L5LH0fqla1KLzZOPmfzxpIOtw7VOCNQtoyu9V7nvl9FfswXVbLA2qt7Su+R+Vxiwl6++afNCnOXmJvbK7EQmasil3bVQ87IrZfTE6fTJ4NVFKq6yn+OeJQsB9Y/X886s/ZuS3FEN0UOynYVJKsKUM1HHs+16c+XMfNaUXJyG7YOmC+U//r2yv+KHBXx920zWSl5+ZlzY+Y5cTor3ZR92bbq+Zw4W5gV5kqo7FcFE0Wq+V6faVfjLe3zOG7sHPRU1eIqhRGdmC+Urap/P4artrI+yoeqTZ4bFxgVVH/m515mVd1lPSs5rD+zwyUwKkFytrkqWM2R0kOt4ewDltwouWptuMSwR+9nVdHuN2EGWpt10a4wCXNLqAcnB/zl/nJNPVyjZKnG6Mdh4+dA6wIY0031ZUGxClqq+mDHFclWCYUjxNxeJQ8qIXF2V1C2qqrQ2T6iC7NDzYtbX45cHHkzmep4VMaIf3JipEjZHbMEk/l+bzG3ZHeGSZgr4bLnKtvM96pqR2XFVSXDAltVcaj7rH9VHVfVnKpoVVDqZWYows8VSfaNskXJUnbmMRiqxCfrwfylKrgRAmPVmUtKXDXmCKca3yUdLhHLOvbXmM5sPFcZumdB6enW1L5hn3U7bZiEeQioAANc/jD3YG2ZXNaHBdFeRn64swx27AKECtAu8CofMBsUOY+M4WzNurPrijjV+Pkaq0JGfJDHzDa46sWRmap83LphqPzG5sKt84osVYWnbFd2MJnO30on9Xyq9cpk7Q0a5pt+dohJmCtRZfgsmKmsfW3l4AhvJBBURMBsZXZl2/I4LACqakbZw2xzeitZo3oz//dYQyAuOWD9lD7KL2o+nW+qtVX5jflZ3XP6u4RiZD6rJJA9U0pnpQs7zv9VUrlvaPM1zJ1iEuYgHviR6LqtykhVRcSCgyPAXoYiIKZPJc8FmG10YcTtSKMicKZDpY8iXyavItrRSlgF1hEdWX8mT5FLtmu0vbKV+SbrXxGVI94q2WP6KVtZ/zw/zsbRBIHpH9FO349LT6zGJMxDQlUVayqP/l4P94Cz8dSWEcvQXYXFArgi+ywrQ1VquSJSNvf6ZTnKLmYHs7WSmWUpPzOCUwnRGsJj4zhSX2OjSrbYmnMVHvOba+t0Xu7nc3ac9XPresQnymb3TDmZ+4S5Jbs7zFp9BdiDph7UDBZ4FfqgnCuOTMKKSHudFQG6rLzvmysJpkevOwvy+S/r7PSuKiwX8LNdTO9MglWgV8SY5aoxHNj4LLlyZMvWopPJ5kCRSl6TjkDzX26v1mVuo9aeIrEqEcly8/hqDlnCsvdoD6yXtX8jiIgbI+KPIuLuiHgJuf/wiPjVzf13RMR1m+s3RMQ7I+K9m/9f3/V520bmuzd/X7gjbxwas8I8BEZJMj9oKmt3Aau/z0ihqijUedazl6v6jhCoy+wdIY3o5mxXpKYqlxF/5/GYndmuTOhqDNXPBazKFlWRZbtYEqdsdgkJI1tVfVUkw+afkWyV5GVb+3Yjvs1jVf2YTvuCUfJbi4i4AsBPA7gBwEUAd0bEba2193fNXgTg4621J0bEzQBeAeD5AO4B8E2ttY9ExJfi4HeQr+76vaC1dteRKH4IzApzC1TBm2XWgN5mWpu55vasX9WGtWWBRpFzhnoolQ6OnNf4hlUZSs9M7kym8huDGpvJq8g2I88JmyNla05mqqRGJStOL3ZdtWfHai7ceKPz4sbMMvs2breA6VP13xe0S7HV3wCeBuDu1toHW2ufAfB6ADelNjcBeM3m+I0AnhER0Vp7V2vtI5vr7wPwuRHx8B2Ye6SYhDmIPtOttrlyP9XWVaX9mFW7XC2MyHVjuApYtctBoyJfdpz1rmzqx1W25DYMbBxG4Ep/ZwMby/l60XekSsx6V5Upk6FIM8tTckeIVu0osDbKLtaX2a0qUmYb082tS/ZMqPv7guVdstv8AbgyIu7q/m5N4q8G8KHu/CIeXCU+qE1r7T4AnwDw2NTmWwG8q7X26e7az2+2Y38wYuStlseDuSU7iLUPQ59B5+CngpTKaFVluozDZCm9FXG7jFuRIguyquLJpJn7qATA2cGqJwY15mgFnm1jAXJkffQ65/WxRjaT099zst16UGvAkbAao1q3DmxNZlvV/WybW4/VmnRrYE1VfopxT2vtgrnPFkZ2iG0TEU/BwTbts7r7L2itfTgiHgng1wF8J4DXjql8tJiEafDAR0nqDJqRAgserF9/vmAksDFCY8GXyWP9M1TAGNnOYn2cLx3Ru4Cr7rHgpnQZJePcXwVa5hvlt2zDWnszkTBdFImxsRmYX93aqPzSX2M+GyWiKtlk7bJNyk+s3yih9r9odOIfMWk4ynfJXgTw+O78GgAfEW0uRsRDADwKwL0AEBHXAPgNAC9srf3x/Sq39uHN/09GxC/jYOt3LwhzbskOwG3TAPw1ob49y3xZkFoewopUWLadM2ZXsalsu7eN/Tnds21VBp99WumgfKKSEmaXShqUPSPts7/VvCrdWTu1XtQcMx1Y8sR8xebB2TyaPGYbma1qLNaXjavGyTIYsak1Xa1vJofpPSrjOKCe5+pvAHcCeFJEPCEiHgbgZgC3pTa3Abhlc/wcAG9trbWI+AIAvwXg+1tr/3xpHBEPiYgrN8cPBfCNAH7/UA7YIWaFaeCyz6q9q7zWBh1HBnksFYRcBeL0yZWzqxpGqoKsQ9VGXWMVkqocFXFXOrN56PWu7B8N9KzqrZISV4m56tGN56rNao1knV2S5qD0Uc9E1abXy613J8utezbnzP6TxFHp0Vq7LyJejIN3uF4B4NWttfdFxMsA3NVauw3AqwC8LiLuxkFlefOm+4sBPBHAD0bED26uPQvAXwJ484YsrwBwB4CfPRIDtsAkTAO3jdTfZxntNg+OyoadHEcSqo0bn8lgY6ugk++5cZleqjrIJDHi5zV+r2xS9/M9pecanfq+jMhGSdhVdSrwVwkV8zebDzYmk+uqwcpX1RpUBF4lG3nuqmRRxYXRJOFI0Ybf8bqd+NbeBOBN6doPdcefAvBc0u9HAPyIEPvUXeq4S0zCHITKKkeqjZGqIaOq4Fy1tTzkSr8+aOb+ebyqsh6pvJU9/ZhMTjVmbjNC1Kpf9gnTjV1XOihSHyFdRiiV7yobR+71to6QpbNbjeWqNpcUsGenWnsu0RtZ62yNuIRCrZeTRMPRVZjnEfM1TIMlw2eVi8peVZWTHza1VcYITy340eu9DqoaVbayKoDpnMfI/UfkOJ2zTvmeq0wrqODP7mU/9fZXvlP6qP6qD/MjS5KyzJFzp3OutpT/q/np2+Rj1pa1yW2zzSqZYTq5dZTtZzYyXSdJnU2UFWZEPGZAzqXW2v+7A332ClVWyYJpvu7as35Kh6Vdf51l//2YSndW+bg2rtJ0gbC/70hJjcOSjEoX5RN1nMdilY7yBQukTH7Wz9lakXfWSa0bNrfKV6q/s5np0kOds37KD8p+hxF/ZfuVHSrZ68dxOu8Lac5fK9kdRrZkP7L5c7N/BYBrd6LRHqB/WzgjOPawjJDUAhX4cwa8XBupfhixsIc/H/dylE7ZVjWuQxVYRuW5cXJAdwTtZCsf9+fMFkeAzEalu9LJrasqcVE2uUqsv87850iTrTVmg6ri3NrIdjAfjyaIykalP9PXrbPbr38eAOCGd/waTgqX9oS4zwJGCPMPWmtf5RpExLt2pM/egGWWLsioqkfJVXCEsm3gzX1cFq8yekairF0e0yUcfd/RbNyRF2vDdMrHSi47V/4Z0TPrkmUokmBkqJKf/hojkaxDlQww4lUJWbaN2VCtOQW2ZllS5/orO5xspuPoWhm17UjRjvZNP+cNI4T57++ozalBtcXEjl2Gy+6PjpuhsnnWLxPtcq3/3x87m0d1yX0zIYwElD6AsSC/9FOkovyj5obJH6nGnT8cobG2I7rl+2w81WY5zjpU8rNfmMxqTplMpS+zp+9XJWsj99XcVGvR2ZPlj87nUaNhf7aGzwJKwty8LfjQbU4TVBa9QG3tsHYVqYwEs5Hslz34Tr6zhclzJK22r5gtavy+rQpqI+Mv565yY8GMVfGVDsp2V9Xl+8pXbL5dBThCnlUg75MU5RdmsyLWXn9XdTKdFOmOzInyRx678o3zd7ZX9T1JslwwCXN3GP5YSURcAPADAP7mpl8AaK21Lz8i3U4MrArrUWWVTBa7loODk9Pfz8dZryqr7vVQgT63VYFWBatej5Ggmu9VpM0SChak2PVtqz/Wtx9b+USdqzl0CZmqwKoKMLet/MsSCLeOlL2KcBUxMj1Vu9EEqpLJ5prFAEW+WQYbc+JsYM3nMH8JwH8J4L0ALh2NOvsBRzwZKpt3gYBluGpsRqrsvrIh69Wfu+pEEWtlHyOyimBZJs/6qzEdqanrWZ8czFXyUgVupR8bSxFnJpc8hqtqWHKg1oJbI24e1tifCTfbnW3qdWE2ZF+tsTfr49Z6f6ySIVd5Mt1OCpOwd4c1hPkXm686OhdQ2SYLoKMLsspGXZXC2iid+nbsgVWBiN1n/dhYLii7YMsy+ZExnayKlLOOqupgtqlr1T0WtPMYqkpWOir7RuY+91XnmWAcgTAZypasE5sLtdYZCTNyV8/ZWsKuEg41VyopPla0wKX5sZKdYQ1hvjQifg7AWwDc/7tlrbV/unOt9gCM3Fjwr4I0k1VVLyyQLvercXLAUZVLlqUIRI2hdO7ty35heq6pZNgcKF/ltirIsUBZVaKqysh9FJT/FKGPEE5/zGyqxs9+ZTYx36lnwp3n8Vyl55K6EWKsEh62DlVbJ2MZzyVGJ4EGnOt3ye76ewTWEObfBfC3ADwUD2zJNgBnkjCBdVuRKmCOZpfsoRsZUxEUI0pHrCxQ5rZZxojdjlxz+yobV2SZs/k1VaEjyBF93HmWr9o5vdUcLdcyFMk53V0blyBWYyj5LuliRMigyDInHIyE1TxXc8WSGdbGJTAngZMe/4Sx0+8RWEOYX9Fa+7IV7c8cHJGwdq4q6Y+roJD7qDbsnAW1PL7SWyULOThVwc1VK9smIyNVy9rAtWYu1PyrNaH6qfl1Ps2BX/Xtz1m/Ed2qNaySL1Y9ukSGEZ+yI8tzz6Ujc/bcurnr5TufKR9OnAh2+j0Cawjz7RHx5Nba+1f0ObVwQZEFeBYwlnMmr7+XZfTtM1jAyTqp6jHLYZWj0zWPwcbNvsjyquDK5DC7HDmrc5U0ZJnMl2ruRubJ+cr5mhEOk+lIKMvLtvcyHCmqecnjZvvV+so2qaTMrffq+Ry1QSVeTi/1PDBfnDT2SZcTwE6/R2ANYX4NgFsi4k9w8BpmADiTHysBxrYYl3NHKqMBfIQQ1DhO9yq7dUGwCiqVflWfiqiyni4JGNFPVV29TEfolV8Z4SgyyjLUWhiZyyx7TWKWZbiqiNk2ksTk/jl5cslpHk+RG+tXJWH5fMTXoz6tfHFsaOf7q/GW7whwH4tc8z0CawjzxjWKnkbcfv3z6APMwKqAfE/1yeesAmDtFqjKYxudVfDK4zq5GU73/nrWmQU+1s7pkUlRVQqsD9OVVd4jtqm1kWWyfkr3vu8oObnKZzSZYvKYP0fXiyIbtn6V/o4gnX1Z/0pnV3Wzc9a3tcAdT38uWotj/07ZhpN909EeYScfixwmzNbav9p2kNMClTGzIKxIq9qacbKqh7Xvz6ogN04eLxNl1qnv5/TZBs5/WY9eT6eLCwpsPquqaNTPbPyqSspyl3NHxooIqnWT57kfVxG0Gl8Rh6pm1VpnOigZTDcGlVy4xMOtraX9yPMwMte7ena2wSRMADv6WOT8pp8OqrJRmW9V9bnqqSLJKkBnfVnfvo166BVps2ydESzTgxFzVbkx25R/c7/czsll8l2wVnYo/TMpsHEzFCkrwmF9nY2MGPJ8On8o//Y6sgRD6c58VCVtzM5Kl36s/l6eT9VeJSlKp7VJ1nHhPG/JdtjJxyLnN/0IqOxzudfDBQBGRH2f0fZOz6yvCzjsIXYZcBVEq6BSBV9HVhUcoaoAzvzqCCGP5+SqpIqRUt83H6txqurJ+bcK7lkH5g+XoDh9WEJVyWNQ83vYta4IvfL50rdKILOeEyeCv4sdfCxyftOPAAt2LgCqe6p6cuP0bVTAUf0qnZUuI0FjJICM9GN2VQTAqocRHXLyUfmL3et1VmOvnX+VnChSZcSgfFn5lMGtGyc3+5bJUteUjQsJ9XPuyNCtdUey1XobSVBY29F1cFyYRA1gRx+LnN/006EiGRV4XZv+Wm7jAucC9RCrh74apyIEJqsaR/mksq/3t9OByWDjLroq+1ngZfqMJCnsnkoMer2UDkpP5Rt1jdlf6ZSvM5+wNsqebKsjffX8VGtp5DlgY1QJiLvW93PPYrXOjhOtTcLcYCcfi5zf9EPgFnufAY/KycgyXAatHvqsA8twl+OKtHq91tijSKgK1NkPalxFuE4vRSr9f1dJjPRhOqokR+mpCEAd58prJBBX5MSuVQlatmc0wVS6Op8627KNvT5ujGo+1Zi97GqNsfV8EmS5GXm+hnmAnXwscn7TT4cqM64C40iAr6odRR7sgT1MAFeViAuKuf2I3VW16eD8oqoCVsmpCkTNn/NN71sXHBl5uMSoIrNRglTtWMXGZDk7HSm6yt6B6cH0dH2z/pX8/p5KtnLCmav4as6qZ+44MStMADv6WOT8pp8OLqj1/1UAUcG2WrBVNah0YHq7MVV17IK7stnZwIIz800VxFTwyQSeoZIFFsiVTaytIvAsj8EF0JHgOkIO1dpZ2vQ2q4qyIvy+DUtYKrvZM6MIVpHyyPPJbK7m0+nv1qSTdZKkNQkT2NXHIuc3/SSooKOyRpehs4owB/vRyqvKtnuwANnr6rJgFgQY2Tm9XQLBgjLzS9YlXx/1W4WqIstgCRG7PjrHal5d5cbgfJdlublnCQRbx71cRbhZD5Wwjczj2vl25KjGHH22R8h3l2t04vCIiIcD+FYA16Hjvdbay9bImd/0k8CCiqpAVIBYm+m7Nss1FQAr3VzGP6pnL5MFvSowVMTTyx/Vk9nI9Or/jxKiqr6ynqpSYudZtiO0/pqqorIPRgO0qrqWa9k+5i+1Rvu2GZW/e1mVvkqeOnfrQxHdqG4sGWX3lYyjRsP8HOYGvwngEwDeie5Nq2sxv+knQVUcjKRUW3Vt9AFiD+1IgKgCQO6jgpgKThVZVPo4O0cqAOdzlcAo2cyOfhwmv9J/tEKp5lPtJOS+I2uQtWUkO5KojewqZFvV/Dr5zj6WsI36lek8msTmtZ71cEnZiVaabW7JbnBNa+3QRV/5U9wR8Xu7aHMasTwYLvPv2+W+/f+qX1VdqYDT6zUaEPs2eXyXMavsuSLeUbsyKSs/5XNHCMrfSo+RwNbLdTozfd08qsqMkQST5eS45CTfd0mMW7NZL+Z/1i/rW5EM68/6sLXC/jPfOvsYRhLKk8HBu2S3+Ttj+L8i4lg+h/nvRcR7zP0A8KjDKrJPUAGqvw+sr4KyXFUNsQe1qijYw8/IksmrAm22S5Egg9NFkdQI+Y70GUl01lR7o3ZW91mFxOSPJiJunSq9RqpF5p9tgj/rm9dhpdeaKs0ljCNJHpOlEueq/1rdd40GoNnfTj43+BoA33UcHyv5WwNtPrtmUIWIuBHAP8bBL2D/XGvt5en+wwG8FsBTAXwMwPNba/9yF2M/eBy9zaLOWZa8yHDBnGWmbmwV5FWQYxWVqhJGAm/GCOn311gWzqpDdqz0YvopElQk4MZ2c+nGYb4Y9TGbR5UsKf9UiVsFlWSotiwpUvrmeR1JNFXfvLZzW0fwbN26RG6E0N1aPQkcZXV7mJgdEd8P4EU44I/va629eUTmlnj2DmTUhHlcr11GxBUAfhrADQAuArgzIm5LH2N5EYCPt9aeGBE3A3gFgOfvSgf3sKlgubRT2fNoAMljVIF4WzuyTopcXeDINjByd0F+JABXgc8lBspPVbBT5F3NPfNVf935dqS6VcmP04URNrOdBXmnnyNtRZrMR8of2d6sh/IV80ceP9ub26lx+35uvGyT0vms4DAxOyKeDOBmAE8B8O8AuCMivmTTp5K5GrvisfI1zGPE0wDc3Vr7YGvtMwBeD+Cm1OYmAK/ZHL8RwDMiYmfpEwse7CFn7fuHwz0gLkN3GWoec/R+RVJOn+U+s5uRlJOhxmKymB2uYmX2MXsVeSg9R4OzSwB6/zFdXQB2PmXXlZ8r/XI7Rbj9f3asxsiyXVuV1KhrKknt+1T6uLXD7meZ/XNfze1J4AhfwzxMzL4JwOtba59urf0JgLs38kZkDmPX78FZ87GSo8bVAD7UnV8EcL1q01q7LyI+AeCxAO7pG0XErQBuBYBrr712WIHRrYuq8mDVipOTZR7Y4IOzC/jMHlatuX4uo2YyqyrB+czZwIKfSmLy9SrI5raVz13Vm2W4Kqm3U80nI6zsn8oPKtlRxKrmkNnJ2qtxXD/mk+xLNYZbC0wPtQ7dc9zrVT3XrIp1ycrx4FA/IH1lRNzVnb+ytfbK7vwwMftqAG9Pfa/eHFcy12Cn78FZ83uYjwDwHADX4mBv+UFY+wFQNgS5llflSBtsJvWVAHDhwoXhNG9bEnSBIxOhG7v/r+SMnFfVQT+OIo9errpX2cqIV5Hqcs8Fm3zf6VzJcTqw8bKd2XfKXjUmG0vdX+CINJNkJrqK5FQiUSVPKumryJQdMz3YGEs7tSazDmq8Ub8zndk8qDlhYx0XWjvU5zDvaa1dMPcPE7PVdbbreRjH7fQ9OGsqzNtwUMn9Ho7m9zAvAnh8d34NgI+INhcj4iE4yAzu3aUSOQtlcBkwe3iz3EVGL6uXneW4oKCyZUbyyq6KDB2RLvfV+GzMkSDZy1I6j/hG+SmPk2Wy+3k8RiZ5DHbsyC4HeKf7CCEpHUYDeG7Dqjpmc/ZZ/7/SO9ub55PpwsZWuqv5zH1ccjUyx0qv48YRjnuYmO36VjKHsev34KwhzMe21m7Y5eAJdwJ4UkQ8AcCHcfCC8LenNrcBuAXA7+Kg2n1ra21naRsjQlaBLEFiOc/9l+P+Nd3HEwAAIABJREFUv8p0q2BcEVXWQclh/avKLI+hfMFkjWbTrHLsiW6E+KsAyYKVuqd8oIi/SqrytewfR4jZt0qfbZIVNj+VXGajsnvEh0y+IkY1XtZf6cN0ZmuM6ezGHk1wRxOUo8Clo/tYydYxOyJuA/DLEfE/4OBNP08C8H/joPKsZJ4Y1hDm/xERX9Zae+9RKLLZ334xgDfjYMv31a2190XEywDctfnx6lcBeF1E3I2DLOXmXergKocH9PRbmv3xSIaqMnNHCCrIs74VwVWZOwsUa4K3k6euMZ2YT50ePfrrlS6qQmGJA7Mjz6uybU2CwXQaIWNmk1qvbKxqPvtztdadD9ncrfFRP26e80oXl2Tlfmvvq7msCP+04TAxe9Pu1wC8H8B9AP5+a+2zAMBkHrdtCmsI8wYAf+8ov3y9tfYmAG9K136oO/4UgOfuarzLx798G4kFypEqo8qaR9tUhKuqkKyzkpftzkFLkZK6nuVlH7lsX9nB/MGIk43Dko5tSEMFQhfsmd6OeHu/ZB8rKALt5VS2sT4uaVza5Xv9NTafeR1Wz032p6vynP2KWNfMrUpAqz7Oh8eBhqMl6sPE7NbaPwTwD0dk7gIRcdVG/l9sK2MNYf7tbQc5LVABVz3A7GFTJKLG6/v047H/7D57MF2fPDbTRQUC16/v63znoAJM7qcSFkcQKlgr3zC7lM65HVsDo/+z7SPz6UhLzTmbj2qts7EcEfXIa92Rqkv6HFkxUhvRhcElYZWOI8/+8eFMfs3dMDYfYXkpgBfjoMj7nIi4D8BPbfNG1fnl6wQsACzXc5v80FfBjl1nD7B7OHu9XCavgp0LtMoHuR8joawT01UFzIp4lG67CLau0lpkqHFccuB0cbZksHkabevmahQj/lJjO5+N6OzWBSM892yxJKJaH4pU3XPJcJIEWiWqZxz/AMB/AOCrN5/3RER8MYCfiYj/tLX2E2uE7dPnMPcCigT6axWBuqxWyauyU5cluwC+TbDrdXSBlhGoCwprAiZrl/tUfmNJgOqvkg9ld2Wn87EL0qova8faM/+5NZv7MZ+qftX8qPauSnTPEbMry1PrR9lSkZ7ybR6f6c2Sw+NGw9F8pOEU4YUAbmit3f9Z/dbaByPiOwD8bwAmYR4GLDNWVUTus6AKfrlqYQFqdIunevir6lT9r4hkNACsIVJljyJDRzprwYhIjc3GrRKL0fXh/J7bsSo7j6kq5SxXJRDbJFyOEN31te1Gk4B8r3omchvW3813tV6PFe3cV5gP7clyQWvtLyLioWuF7dNX4504+iCTj1nbBTlY5ev5b7mXgx4j5Dx+Hme0X29f1pHZn8di19i5GjvL6P8UoeZkopehCCAHUWYP65f9wZIEpofzQ98/68d8kXVj16t1pOxm+rv1nW2q1mAeR9nmfJUTN+evbAtrw0icJTnqOWFJiPKPa1PZPXGk+MyW9yjKCjMiPgnQb1oIAK219vlrB91nsCyd3a8qj3ydkTGT4477B09VFFkPpZurWJhPWFBW/ZTt2Q/Zp2w8R/BZJzUfCopEKll5THc/Q817P76TkX1TJUWKNCqSzX5h/dh6YTa7OWJtqrZsPLX+1fph/su2qSStSpCZHf2148Z5ftMPgK+IiH8N3P9h1GUSAsDnrhU28mslj1wr9DRDVRMs41zuMRnsWD04VZBXhOTInMlw1WSu5hQ5MllsHJUUOFtZUlDJZ6TZQyUPI4StUBGVShiYnL7vmuCq7q+ZL5bAKTkV+bBEjyUgI35l/9WcsTVWrR1mk0uW8jjZP0x+1vMkcZ5/D7O1dtnXuB4G8zXMBJZFrwlk7KFVQWC5t0Blrv091s+1YzqzvipIqIpX9cnyR6rgEd/kPllmH7Srazkw5vZZLht/VGfnB0b0le1qLnI/RXiOdJyO7L7Tq/K3gkrqHOG7Z9IluGwORhJBp0OWX62Lo0abHyv5XADfA+CJAN6Dgy9CuG9beSNbsv9na+1ruq3Z3vtnakuWBdn+nmrvqgc3jnqgquw2y2APLLuWs+t8XAXbPBYL2K5SUeSpCFxB+dZVQI78HWGxZKNKEkaqGzZWlsUSN9Zejd/LcZVR308RYdYhX6uqK1XdqcQy6+DaMTuyHxhhV4kj01GtH6an89tx4tLJDb0PeA2Avwbwz3DwXQJPAfCfbCtsZEv2azb/z8XWLCMLBhcgVPbvAq17kEcDsMuClX4sw1770FdVAevLyEpVJ44MHBTRZ517Pdh8sPbK/mxHvu+CcH/f2ch0ZeM7vUfnuBpDJQG9rr3OuY0joDyWm8/cT81dVbUq3Zlc10et5ZPAed6SBfDk1tqXAUBEvAoH31e7Ndb8vNfnAvh7AL4GB5XmPwPwTzZffXRm4AJ9FeQYUSxtKxKtyE7po9qwqlERp9NT6av8NiKD2aN8k6uUEV36tmq8fM6qvCrRGCERVWUpfdV8sPms7HCJ1sha66+5hION6YjQJaGV3ur5yqhI1T2/arwsQ60nVX1OnAj+ejnYfPftoYSteQ3ztQA+CeCnNuffBuB1OMLvdj1usKCkHhz1sI6QpQswqk9/PZ+rYMtsy2MyQurPFQky8siE4QIas7OymREY05e1rQhF+YuNNZpU9XIZQSmd87yqe0pmtoGRMLObVUX5OvMV61MRhHu2su5Oz3xt1JasC9NLnTsbnE0ngYb5LtnuXbINwOdtXl4EgNUvKa4hzH+3tfYV3fnvRMS/WDPYvuOGd/waAOCOpz/3sqCk/gNjBLi0ywQ0Uh2wAKGCE7vugpIjjVFSyDqOEmx/j42v9KjmQvk++6LyiUsKRvznqr0qyPZjuOSL+Ubdc+uLXVNtWBLC1r2SWa3LHi4ZYwnUCNG5sdT47JyNx3QGgGe+/Q3l+EeCBuzuBxBPH5Z3yUbEI3Dwxp9LAP54253RNV9c8K6IePpyEhHXA/jn2wx6GuACsws+gK9KFBG7/jlI9Q9pLzOj0pEFEVZZ9WO4ADdCYlle1o35rrfBVR6ZMNQYqo2S24/N5jDLqRIAJrPSi+net1UJnBvX6cn8qnzI1iA7dwkF82MeU8nK+ql1pZ6zXmZOVFwCwHzOdK7m46hxCbHV31lARDwkIl4B4E8B/AKAXwTwoYj40W2+6WdNhXk9gBdGxJ9uzq8F8AcR8V4clLY7+5mvfYB62F1GucBVAlVWroKsCxwjlcFIHxcgFOG6DJz5iunI5Dub1LWM0YpJtXGJSm7D5pXZxXTMhK3kVXNbBXAlx/mkv1dVa85uljCxZDGTXratTxJZe0VYOdFx+lYJITuv1tFJoeHkCfuE8WMAHgngi1trnwSAiPh8AD+++Vv1jtk1hHnjGsGnHYpAMtRirPpv81BW+jACZ1VMrlxZAHPZen9eVUnKD0yHEf2Y/1jwZfpuWwWMEjOToRIWZXN/nufdnff/1Xy6ZEn5mt1n9uRj5gMmx5GzWweZAJ0O2c9qrVdjqzWl5tOti+PD+f4cJoBvBPAlrT2wMd1a+9cR8b0A/hBHSJi3sIvb/KbYvqMiLYBv/7nz/pp6YFUW7AJhHqu/px70vh+DCwwKKnAweUwus4GRTf8/+5tdY/5T+uZ+uY/TmbUdsZXpqWSpOagCtEqMMoGtSQxGfL9cU+SWdWX3KyIf0aG/7hIrNU4P9cwqAp44cbSeLLuLn40tJmrNa5h/2f19FsCzAVy3dsDTAPUQVQ98f56Dfv7PxmFycmDr+1ZBt7/HgmivJ9OJtXcysw7sfs7wR4K0qkiyDiphUZVN1Sbb4PztgiULqnl++75OVpUwKbJ0trC+bO0qXdT67NsxuUqPbFP2Xz7OfTPWPGf9eJXc3H9k/JNA2/LvjOD9EfHCfHHz815/uFbYmh+Q/u/TgD8O4La1A54msEprAcvQ+3t9++rhcxVMP9aonqyPqqJU8HFEzWzL+oxWR/keG9dVZ6p/hSrRYVWOOmfzz9bDCBGOnDtSZNfUeH3b0T59osZ8zhKAfK8nPaW/8sOI75gPss2qKmRJsUoolN9G1+BxoOHcf6zk7wP4pxHxHwN4Jw5c8tUAHgHg76wVdpjvkv08AF98iP6nEuoB68GCgurbt+/7K1kViboKj+mQH3BXWbC2jOCz3rkvC2KMENQ4zH5HYEoGC3T9uMx3jMRY1ZQD76juyu7cTiUYa6qwimidzVkflYCsgUpKGLmxNVsRF+u3Rne2hkcSu218sUuc5x+Qbq19GMD1EfH1OPhavADw2621t2wjb803/bwXD1TqVwC4CsCZe/0S0EQ2EsjzfyZnZHxFmq6dIuTcVhFTD0fijvBUAsH6quCcA5kK7Fle33Yk43fzPFIhuMpHEasjHzYfao7ZehrVpz+vgrwib6enSljy+Cx5Y74YPXbynL5KT2UDS/KYjeracWMfdDhptNbeCuCth5WzpsL8xu74PgB/dphvfd93MOLLgaWquPprKgC7yssFoF63rI+6P1KpKNtH+qmAq/qxMVifarwqw2dyqgqEBV8VTN38975zwbUiC6Zj396tlYrIWNLHSMElEo6s8/XK9+qZYLId4btrI75SYHqN+uC40dq535LdKUZ+raT/eZT3AnjVWSZKYF2AqQLjSBauMmcmQ+lQEVcvq6qgRqvF/lzZrWQ5PVSmzpKCXhbTT13LOrBgnsmI+bmXrxKN0aCs/D1KVMpWtU5Ye5Y4OdJUCYbTecRup1++XpGTS+TcM5nbV/7Ia716ziZOH0YqzNfggZ9HeTaAJ+MQP49y2qAqDVWNqMycycxgstj/kcrDwRF01llVRqwq6du7zF7pzIJhX3EpnfN9J3vEfpeAVLJYG6ZPj4pAmA1VgpTXqxtbran+/hpfqnZZHpM/Mg6Tw3ypEp/+vF9fTEa2S81Bb2O+p57348Kk7N1hhDB3+vMopwHL9z4u3ymrSEVVM1VVtYARjiIWBpU1VzoqOb1Oa2xj16pqhgXZUT2yrUw2I1k2LtO1SjxGyF/5L9vvfJf7qmRGkXa2xRF/lUCpeVXEldu6tcD69G22JeSR58YlQv11Nq+sTdZv+X7qk8Tckt0dRghzpz+PctqgHj5XBVTZ+nLcX+9lKJnV2GwsFTgcebskgY2vrqmApXzT93Pjs/YVWTE5IwGcya+gyEOdu7WQybG6N5IwZaJ1/XOfPH7lkxE92P018+SeUWXLiE6KqFl7lhyMrpejRMP5fpfsrjFCmMvPowBAAHhE/3Mpa38e5TTCPQwjGXRVbfTkUFUca4JURerMnlG9+2BSEftIIKmuVRWgqnwcWOB0FXTfduSekueuj/rGoaqs3DVXWbl5dD5g8t1cjbQZtZUlhRVRMtnqeXfrYD/Av8R/YjuUhLn8PMp5RCaFKstc2uZr/bkjEVe5VVmyIkFGCiogLfeqijbLU77o5bI+IySibHFVkauc+3bKp0oWa1vZ5MaobKiOWULl1hjTj9mVE7TeFrc2qvWkSGVkbbN77plkstnzMPIMuHXKMJKsHSdmhbk7rPlqvHMHlkXnhyQHc9aW/XeBoyIZpduIDfkea1eRY2v6p6eUTBaAl3NXneQgp9rkOVAJjNNntIJRAV0lPFme87NLXHKbSk8mw9mTfcH6qfGdP1lf1p6RYT5XiYzyqbOLXVfJghqr/5s4+5iEaeCCO2unHrT8YLmg47JWV71lHVkg6uXksXrd2XHfpvJBlpH79PcrgmE2j+jlEpj+XNnC5oIFcFdNsMoqX+/lsqqGtWV2jiRZTD9mG5On1rgjRQZGbGr+nM/YOO7cEbPSQ5GwS4R73ZwfjgsNlz/zo38Tl2MSpkG/8Bkh5ICpHviKhNiD6ioGJ58FWZatq2prtKJh+vd92DGzLfuGBRtGBmwM1o75Jo/V980VICMuNYdZ9yynmqPcjumqbGBrRhGI8m32V/azW+uONNgzxOBImj0bjmDzfbZGlG7uWc26MT32DZfadn+HQUQ8JiJuj4gPbP4/WrS7ZdPmAxFxy+ba50XEb0XEH0bE+yLi5V3774qIv4iId2/+vvtwmq7DYb5L9sxj+XgJANx+/fPstkt+ONdkaSzQOeJygZIF10yOLrD0cnJ7RrBZN0VoilSzPo58XV/lGybDtc12sLEr4mN2rQnmaj7UmlA6ORvyOEy+updtXEsYzi9Zxohebj7ZGJVf3JpUa6q/38eNfcAJbRa/BMBbWmsvj4iXbM7/q75BRDwGwEsBXMCBmu+MiNsAfBrAj7fWficiHgbgLRHx7Nbab2+6/mpr7cXHZkmHWWGugMqQWSY68mDlQO0qhZGMtgouTB9VBVT6L/ccGeVxR4Iru6+SD5aYKBJg/mVjsUqk0peti0rOSELl7GbHSl9WVVVjsIQp654TIbZ2mL+d3Y6cVbtqTakqstfPVee5jVrjlW0ngeWr8bb5OyRuwsGX3mDz/1tIm28AcHtr7d7W2scB3A7gxtbaX7XWfudA//YZAL8H4JrDKrQLTMIchAq0/T2W5bO+/TEL1I5scp/R4MnA5FSEksdSBDsSJJk9zi4VmEZ9oyqI5X6WqYK3ks2CLtN9ZD5d8lPBVVT9mJk0ej2rte5kZ6gqmLXJY6nngyVGKiFzOrp1yezo5bH/287ZUeLSln8AroyIu7q/W1cM+zdaax8FgM3/LyRtrgbwoe784uba/YiILwDwTQD6Xxf51oh4T0S8MSIev0KnQ2NuyW4BtYWnHnAWfPuH3BGteuizHlk/pTeTm+1xgZ9dVzaqdiNBhVUwWQZLLEb0VYFZ+Z/NmyPDPE4meqWnWyvs2K233D7fZ+2reWQ+Y35ixFKtsTyfa9dgPneVqrPX3R/x9xnDPa21C+pmRNwB4IvIrR8YlM8C1f3OjIiHAPgVAD/ZWvvg5vL/AuBXWmufjojvwUH1+vWD4x0akzAPCRWc8kOmsv7+oVRBhwUSd5x1YLqygKTOFXFn/SviyX1cf0aWym/ZLuZLNheZlJ0PsnxFDsp+do8F+uzj7Bs2H4u9Slbvj7yWKjKrSHTEH6qdSgqU/RnV/FT98hrL/dn6db7ZV/JkftqN3PZMdS8i/iwiHtda+2hEPA7An5NmFwF8XXd+DYC3deevBPCB1to/6sb8WHf/ZwG8YgvVt8bckl0BtvD6YKUeFhWUXRac+zGZSqfcJh+zIMHas3MXvJ1+KkFQOqign+9lKH+wOerl9/Y6nbPNzE95HEZYzFd923ycxxpJhFjb3h7mHzaumzOW8FWJSh6nb+cSst6WkXXf6+V0Z/PC1tua53Nf0HCoLdnD4DYAt2yObwHwm6TNmwE8KyIevXkX7bM21xARPwLgUQD+Qd9hQ74LvhnAHxxe1XFMwtwBqqpnaaMCr6qGljYs01Xj5/P84Ds5a2wdCd6qfw6QOZiN9M3jMQJTpNbfW45VwFRjKb3ZNZV4sDZMV5bsuHmtKicGNn5FeEoG8527ztZrXstqjVe+UQkNO3brhaFKlPcFrW33d0i8HMANEfEBADdszhERFyLi5w70avcC+GEAd27+XtZauzcirsHBtu6TAfxe+vjI920+avIvAHwfgO86tKYrMLdkB9H/gsmCqvrr24w8tKr6zNVWf809rCobZ4HWkagKlirrV5WZs9klGE5vZh8bZ/TeqG491Pz0dvTtRirl3DePXemUydMR35rrFXGtmQu2vlnb6r67zqpI5pu+/dJGJapMFoC9+GUShkv0pcKjxWbr9Bnk+l0Avrs7fzWAV6c2FwGudGvt+wF8/06VXYFJmFuiypT7cxUAWEDrZTGiyw+3atfLUA9+1k1VTYwAXTDJ8rJfHMGwfupa1kkRkiK+HExH21cJRPbLCAlnv6jdAOZTVykxAlPzyPRw9jBd85jKZ2vmlumfj0dkq+tq3tQ6Zf51ic9JouHwX0Iw8QAmYW4JFbBcsHBBsqool+MsS5FmHiPr3V9XAYOdK1msrdKZ6Z0DVeUbRi6V7fma86Ga02wbu5+TiWreXEKjxmbyFhksAXCk4GxWflG2u3mpkqIRP7sES/Wt1osiz76tSixH1uDE2cEkzC2ggp7KiBmxsDa54nEZt3uQ+/ajY6p7Ixm0Ik7VZ6S6cnb1MlyWPzqOIixGZq76WuNjV6VlfVWfEZscyTDfsHvbrO1MuEyGOmZE2j8XvVw2bqWnkp3vOzJ0Cda+YQevR05sMAlzJVT15wJ+f01t81Sy3EObg0nWbdSG5d6IDKVblc0zGyodVSXiKgMlw5FNLzPr5uZIjeuqlTyGIzB1n91z9o6uCUXSKmlxOqndguwjNgdrE0ZXsW77PLBrbP3tL+JEXsM8q5jvkt0SKkBXD5B6oEf6uaDliLyvFF0wV9Wp21rLuigS7fuxsfp7TG7fjhEf68uqQqeXkufaucowz4eqvtzYTDZrV80xO3ZjqHXa31dz5pIktkbZemU6q8Qgj937VyVJSs/sz/z8qOp7b9EOKswTeJfsmcSsMFciZ9wqcKsHy1WF7FgFMhYwqopqpLodraKYTipwjlQRzJ+OWFy1s6a6UxWLm4uqInE6KDtd5cT65vl1laRKkipdXEXoxlHzPOJjppdKNFw16Oxi/mR2jMhh+u8TGuYPSO8SkzBXgn28BLg8EG9LNi7Qsmtq3CpougpCEXrWOetdVcAuyOe2aox87qqELDtXm9neXElUQZAFXVXducQq+4b1zT6riJ4RqEpsnF7sejWOkjGSnChZy7XeV1mOmgP2PDh7mU1qXUS0vf04ycTuMQnzEFBVkspKRzN1dc3pwcACBiMpFghz/94OdlwFYUWGTAc3VpZX+ZIlFMq+Kjiy9su1qrpiuvUyFRxxrLWX6aHWpOujxnGJorLTVdRqXKaHI1Alc6S9Sjycz/YN82Mlu8MkzC2hSLAiv5EApPozcnZVmSPSKotWOjqdXPBQla+SnQN9VWUrWaPE63TubRv18Zp1kdspwlEVFJvPKiFTurLrTEZFJkzP7LsRomb2qLaV3/o2rHpm42Yd8liVb/cB+6/h6cEkzB3AEVN13wVzBxZA1BZh1rMap+rfQ1U02S4XzBTZq+qMBefcJleGOblxNrPrLGiPkGLWK+utZCuyUSTJxnRtXQLg5oXZk/srm52uTpbyzYhNyjaVCLn5PA3kmNGAXfy25cQG812yWyIHpIwcsPs/1qeqdPoxVQDJslygVLr0fVWwZnqNBmRmW29XVaVmnTN59pVLb5MKjNnuDEZUWRemX76e/48kWVVV3stya0vpp+aVrW1HUnmts/8ZauysA3uOsny2VnudK5uzjZmsWVun+75hvkt2d5gV5pZY3vwDHLwBaKRCrCoJVR2OEmAlr5elqkU2ngoSrh07z3r17RipMbsrvZ2P+nPmI6aTGivbOFKBVv5bjrO/HXmO+KUf1+nCCL1KrLJ8ppdKFBi5ZRlZFyZPPTNMB6e/WwsZ/fO/75jvkt0dZoW5A6iHPf9f0245HqkuHBRJqyDGqqRFj/6Ptc33WPXX33OEzmxQwctVEL3OrArpz3O7Xn62mRFVtj/rnturhKjyR1XZsPlUejId+36qPdOXJTpVsqPmxCWL6rnI+lRVLhuX2ZnX78T5xV4QZkQ8JiJuj4gPbP4/mrT5yoj43c1Pu7wnIp5/EroyjJKWC9pKnsr6c/tMZgyjRDwiI4+rZLBAPeIHFZxU4pH7uWqMja2In42vqjRVsamkSc1nTigURgO4s4XJzOupIk4nU5HXSMWb242MpapPlQg52aNrfZ/RMLdkd4m9IEwALwHwltbakwC8ZXOe8VcAXthaewqAGwH8o4j4gmPUcQgsaI6Q3kgAZsdKFhuHBaeRQDKyhZfR32PVYa5C2NabsqeqVHI1wGyuwAJvtjUfu2pF+XjkXNnc25rXGUOeA7YL4OaVVVlsrbj5ZAmPqnCzXKWDqgrZPad/9iN7hk5jpXlCPyB9JrEvhHkTgNdsjl8D4Ftyg9ba/9Na+8Dm+CMA/hzAVcem4SBUYGPB1D3oqt8Ct9WlyJUFXha8qqy/P6/IuZfnAqLT0ZFg1o3JzcGaQfmw/9/Ph/J9hZFqvKrkquRpdP5HK89RXZjuI/Kd3Iqg1P1MfmycXtdqnS9tTl2V2Q4+h7nN38Tl2Jc3/fyN1tpHAaC19tGI+ELXOCKeBuBhAP5Y3L8VwK0AcO211+5Y1cuhvv1nQc5iR0hsASOCKoBUwd8RpyMcN1ZVYWXbmS96ua4P227L4yjZapzKB8r20WqjCsJZbq+f01mNo6p6t95GyJDpnOUwrPWTS7ZGfMcq0+U6m+cRn0S0U/VmH2CzJXvSSpwhHBthRsQdAL6I3PqBlXIeB+B1AG5prdGdg9baKwG8EgAuXLhwrOuFkYgiwJzh5od5hCyz7BwIHKmNjpnHUTYzOWxcprcas9KPtVF6q2DsKtfKFjW20ysHalV95/bsvLcrX8v+ciS/Rues5+gaVz7M97KMfK/yZa+HShZHEiql02nDrBZ3h2MjzNbaM9W9iPiziHjcprp8HA62W1m7zwfwWwD+69ba249I1UNBkdNotp4f8P46C0oVMbFxVOWW5eRKjcnKxKV0Ubo5slF6Or2znrldJklFBs5/ri3zI+vjgrYjNdbXkVmWWSU1LBFRSQKz0xGnItDKfvcsKJ2Vf91ad+sm2zZxPrEvr2HeBuCWzfEtAH4zN4iIhwH4DQCvba2din0RltVX1Yw7X5CrBlcxuPv5nAWZqrJj4ykblD054DsSU2RZXV/GUUFUkVS2L/u194eqblhlx8Zi/RXRM58wm5hsRXBZliLiET/lsVSl6KrGPLbTma1L9Qzktc70yeOe1uoSAOa7ZHeHfSHMlwO4ISI+AOCGzTki4kJE/NymzfMAfC2A74qId2/+vvJk1PXoM2n3oOUHNxNhruyYvP4aq8BYvyro9bIru9ixu59ls3tq7P5/lu0qUTUeC87ON07X/k8FdOYn5Rc1LrPNEWBPAJVsZY+zVcllSRCbH5foKPJWc1fpq5KD3C/mHBaKAAAgAElEQVQfux2IfUfDfJfsLrEXb/pprX0MwDPI9bsAfPfm+BcB/OIxq7YK7Nt/FqgAl+9VbdU1JWtpv9xzlQ0jH0YorF2Fygamo9I3y2Ak4chBXXNEpPTMfmJBm/msIjnlV9dWVc9rbB7VQ7VdOycqwRsBSx6Y/KUtmwM33tLutP9813wNc3fYC8I8qxgJVss1FVTzeQ6yjuyybCanb6MIIEPdc9t7rF2uMrL+anssQwXobL+6PkJCTFenFyMCtm2Y71c+ZOe9bsruqv/ae8zGKpFj608lDSOJjEv+Fhm5HZsDtQ7PCiZf7g6TMI8QI5VCld2OZOCZ7DJyUFM6ugqTje0q5SowVdWZ0nNNIGPBUiUTWe/c31Vx7jrTX43LbHeVeH/O/KIIb7TCYvcc2Y0mENmWER2U35y+1VpxCZtah6cNDbPC3CX25TXMMwkX1Fi1UgU9J6vHQqB9xcKIoCIHJjsjj8HGdJWTC9qu+lTVg/Nh336ELJzN+VxVcszP+T7Tt7/X+zX7JMt05Jp1Ve1HiCaPn+0dkTW6tliywZ4JJt+tK5ecjeo4cb4wK8wjQv9lBixgu6y2DwqsIstgxMeC6prqpqrmVEWkKhlX/fYB0OnHgmRlDzuv7HBVs0o2lL1Vxex0YlUtq5Kznj1GSKnSw82jqtqdfsxfzvf5nltjbK2MVL6s/Wn7kgKKBsx3vO4OkzBPEIoQMzGoYNkjk7GrAkbINMtSY7NgqhKAbcZh5MNIqZczUjmM+nnU5y5pYUF5DVSyUI3lEjR1jVXiat2M+EatbbY+sm1qTGY38wtDVXXnZOgsYL7jdXeYhHkCUMFLBQNVMTK5KiAw2SrgqGxfBTs1tiK7SmdFMFUAZv7o7R+pHJlOjFhdYHdJRN/G+SX3WUu4lT1r5oXp2x+rOciy1bxV86kSSWWzWuu5Tb+Omd5nAQ3zNcxdYhLmEaPfmh0lQRfce7gMP4MRk6rUVIXmAv4aXbIsJo+1ydfdmMyHrD0jD+V/V9kx/Zw/WDBfro/2Hzl2ycJa/6k2rCJ1Ow5Z5ohuau7cmHkMNn95TOCMbMV2mHy5O8w3/ZwAWKDM1R3LoJfg1wfB/j+ToarJBSyQqDZON2ZHDqT9X9bPBWdWaa7xYdajb++CMBtbVb7KD6wfu5eJp9LBIVdZTIYaJ18fWZfZPnYvj92vKVct5vZM56qSZOtgOVckfpZwEr9WMvIbx5t2t2zafCAibumuvy0i/qj7kpov3Fx/eET8akTcHRHviIjrDqfpOkzCPEaMkNPae/mYkY8jhiyH6aeqBVbBsKqhD1h9PyVj6aNIJ9uo2uYAynTK/dfMiSNPV3m6/s6nh0VOthiJu+TKETer7kYJl629qh0jO0aaLBFh/RgpTxwK5W8cR8RjALwUwPUAngbgpYlYX9Ba+8rN3/L94i8C8PHW2hMB/ASAVxylERmTMI8ZVbB39zPZbJP9syDGgrQKrCOoKqKRoOmCqPKZ8hsjxL49SzIyWamAqwI3q9aY7S6xUHYqP2SwZEFVnIpYXJLCjlnilKtANaaq9tjaGElMXOWaK9RdJyf7hBP6LtnyN44BfAOA21tr97bWPg7gdgA3rpD7RgDPiIhjy3Lma5jHBPa1eTmoMKwhKxbwVEDIAY9Vk052lsu2xnJwzO2UDNanP89Q1UEOrEo35wNW+aoEgJGEst2BzSHziSLbaszKr5W9WS8lb5SE3Brqx2LXmE5s3GqNnbXXLRc0HOpdsldGxF3d+Ss3P504gpHfOL4awIe684ubawt+PiI+C+DXAfxIa631fVpr90XEJwA8FsA9g3odCpMwTwgj20Q5AIxWfCxAueDjxlTjuuCkSLsi0cp25zelF/MDG9cF9VyNMBudHUoW079KPFjSoqo9pncPlQRlOUrfvr+aKzc+G7fSUaFKYtSYlW5nAZe2Lxfvaa1dUDd38BvHbGIXZV/QWvtwRDwSB4T5nQBeW/Q5ckzCPEGMVJa5SnRZPXv4M0E5klZVqJKhCL2yJY/B9GLVQ+Ub11YlHlkHJquqpJheam6dD5bjESJndvW6ZX3yOH1bNc+9PkoXN0aV5DAZjMjcHDK7XSWvErOzjKNikx38xvFFAF/XnV8D4G0b2R/e/P9kRPwyDl7jfO2mz+MBXIyIhwB4FIB7D2/NGOZrmCeAZ779DZdtAfWBqT8G/Gtm+Tz3rQJwDpo9VNv+ngqyuU/fXpGfsitXMllH5rsc6BWhZh1UtdZfy2Sr2qtqLeudSdT5Pevfy2EElf3Ty1vGqsZk9lZJVN9vTRWn5qCX1fuMrW/1PORqOKLhmW9/w6n/NRKH1k7mXbIY+I1jAG8G8KyIePTmzT7PAvDmiHhIRFwJABHxUADfCOD3idznAHjrZqv2WDArzBOEI0XXTgWn/j/bUmMEpCoJFbRYZl9VFMpule0zedm+KoD395j+ioSYHEbKzB5nax4rB33nm6oqHan82LphlZ+qkKvqTvm8gktAchulT2672MLsrvpP7AwvB/BrEfEiAH8K4LkAEBEXAHxPa+27W2v3RsQPA7hz0+dlm2v/Fg6I86EArgBwB4Cf3bR5FYDXRcTdOKgsbz4+kyZhniiWzLZ/E5B6kNWDn+9VgYAF2hyQsi7s/tJfBVEV4FRFxsZjsp197F4mJuc7BxbQFeEw+7JOzG42jrKHzZ3zEUsQVHKi7FFwZMfajWCkMmX6j4y3VJbnAw3tBL66YOQ3jjfnrwbw6tTmLwE8Vcj9FDbkexKYhLknyBnxSGbsAm0ly1W0TDdVRTiCcgFcJQVMVkWWjMyZPFe1jVTRrrLPlaK6V8nozx1puLUxQkxMHwe2bpRflD2Vj5keKpEcmW9m73mrKhvmV+PtEpMw9xCqqlMEUG3PVdtvbPy+b9+ejeEq0hzAXfBf7lfVxagdSj/Wd4TcXfWl9Kh8w8apSD9Xzcw3rAJlctgYozqzdVrpwcZXa92tGyWX6axknBfML1/fHSZh7hlysFBE2AckFfAWVCRaBZgqO1e6uqCW5VVEurRR23CuenGBMrevfJnbOsLoyV+N17dVFaKr7tX4/Vj9MdOf6ex8o+Yzj83aMt8w348kaqqPIkmm53nAMb4n5sxjEuYewAU9dV8RaH8vB1pXdTAyqioiRmhM9zx+VTWxYNoHexXY1diMGPqxHRG7SqsK9lX1ymwdqTwZETJCZX3V2MpvuT2TX82t0oPZPJK4ZVTrU90/D2iYFeYuMQlzD9C/rf32658HQBOMIjRVgShUZOrIyFUprE3fTlVYzD4li42jqiUm1wVgRQysslO6KxJw1Y4iHEbSbI6r6smRRuXbvj9rP2LXSPKidBmpvPMYeeyI8/RGn4mjwvwc5h6CBWNVheYgytqqvlWwGgmQyznTQ8FVWa5dtpHZ5SrNfG0keOfgqwinqoqUbEbIjqR7PVxC5ORnP6i+6v4awnPVupPpKmO2DpjezJfnEa21rf4mLsesMPcM/UdNFrBM2mX8C6oq0FWwfT8VnFnV6CoMVz1VY7vqtNKj2u5juqnxmZ+YHupa7pPHYzqx61U1pexj/UZ0rkhszVrZRjclW+na3z/vleXckt0dJmGeEuSAoarIvtJzgajvk8fIGburnNR9F2xdxVpV0qzSYuMre1g7Fbir8ZXdiiB7VH7IbRhhMB8ov2e7VFXGqtIqmcg6M5vUfCo9mW2j1Ww/5nmuLIHNa5izWtwZJmHuMVQFks+r7baeRPvziiDyeQ5AjDxGSEQFbiaTkfeIL5i+rmrr749Ug9lHjIBYgpJ1q4iY6cHkKWLswZKEkap7NIHIYzCbexuczSM+UXaqBOq8op3AFxecVUzC3FPknwMDdIDOcBk7Ox8N1mw8RQ55LHW/t6XaslPjsIqEVXDqugr+VaXoiCLb1rdhiUuViLhq0lXSVfWW/ZB1yMcqgVE+UDawMbPt7DzboNqe923YHnNLdneYb/o5JRipAvN1V3H211SFpSrS/h6rLNW4CipoMhlV9bX0V8Fb9cnXcpXL/KL65ePctyebXjdFepUt7JhV18t9dS/bovyY+zgf9mP217MOqvplUPbmexMTu8asME8Blmy5/85ZtfUF6K1QlqW7rby+vwr4rgpk51XFxfTodel1qHRhVRuzv6okHZQPlH2qihvRpyKx0SQhX8+JAasm3TnT043dj5ntUuuNVZXs3qwsH4yGhktzS3ZnmIR5iqC2GQFPCKx9bjcakNYQiwr8rIKoyJ21yWMwslQkxoJzbqP8nX2jxmEEM1K9923XklNFqqrvKLGyMZXO+Tq75pIDNga7NqtKgzbf9LNLTMI8RVg+cnL79c+jBAnoIKbIiBGjquYyqnYjlV+W0+vndHfylB+YrxwROf0Z8So/MJJn/lM6KL2YrWoO+r5u7eRjl2BVyVbWz5GkmgfXppcxK0uN+aaf3WES5imEqsIY6anK0G0hZqituBwoWXDNuuWxqgpZ6cL0Zv1dRcr6K38pwmZ+d3rle2wsRkqV/5htjjzVuK4fm0+WPKh1oapFVbU7n+VxJjgaMLdkd4hJmKcYLsD310arJdZ3VAdV+VVbioo4VBWsAixLEno5rl8eJ4/PbMp+ZwnKSPJS6cN8NTpO5VcGZ08PR5pq7KVftsHpqRKefq5GqtPzjEmYu8MkzFOI/JETl8Gz65kY8nFFfi7oOqLrwcbMFYbSjZEVk53HyTqyaiYH5yy3qlb766payjo5Yqx0ynD65P7MFmaj06lKLHKfTHR5bYxs1y5t+u9gnpg4DkzCPOVQVZKq7BSBqSw+92Vt2dgq8KqxXKWrEgJHdjmg52qEEXTWiemixmbblcyOUb9m3ZU9rN3IzoPSr5fDrmW9mX9VdciuMX+5+VbrekKhzdcwd4hJmKcc+Y1AC1T2zo57qOCriDgHOxWwRwOdIpw8ZtY3jz/qByYjtxkhAEb6zBcjVSrzLSMiNZ9OPquuXZKlkiBFrsr/a5KtHuz6/9/e+cdaVlV3/PuliGABZRgZKDJODQRDDcX4gsXWSgooNgbaFKZQbWcaKCXE+oehEWuTGvgHaYw2bWNL1Uhtg/wIxim1TocpGGqFAJGCCHTQUqBQBAdsiBQYZvWPey7s2bPW3vved+697937/Uxe3jln77322ue+Od+z1t7nXC3waceglGyfSDDnBO8CmBKJXbqfU7tA5v3nbSI/vf1IeGuURLwmZJGtkrim9SKRiUSydsMSid6wbST0kUDlUWLeJvo7SX+n58MT3rx+RHQz50Xmpb/T0t+WcCCwm3rXT19IMOeE9OUGwN4XmVIkVks/RnXysigNG4lWWq8kwrntSNgi+9H4hiJUi7ry/koC4p2PSDTHOccppVRt67ksjckTs/zc1M5FNJYoYi7Z0ZzleCjC7A8J5hxSSgN6Kb50P0qd1QQ2Sqd526U2UZ2SrVqaM/LfszM8HkU8pe3WFGc0zqg8Go8XIbbcnIzaf24/7Sc9Vov2vfF45yiyJUbHunf9iH7Qu2TnjFbhqtWJUnV5Gs2rE+17tqKyEqU0oFeWjjmNLD0b0Y1G7q8nEl6UnNf1bNXaRtFja2Set42EdpSI1zv3pcg6bVe6UcnrRGMUYhYowpwz0rTVMD1bukgB9YgtKqtFnVGkVrswtkQgeZ8tPrfOiZWEftwLtydSuW+1z6EmaLWbpbxOafzeDVT6uyR6Jfs1scz7Uxp2+Sgl2x8SzDmnlOryjnsRZBSxlOx4EWnuT34h9qKRYb1RxKSWFiz1kwtDVDcfQ0TJTu5nfu5KNj1/S/1HY6v14/nstY8+n9KNQEv0KZaPFv30hwRzjskXAgFtc2Q18Rn+jqID71iaBi1Fo+l2LnBeWrEkpqXoLLrYlyKiSDwi+16knI83umnwBDXyqSR83g1LdKxkM4oAS+ej1m+UTQC0wKcvBjOYEsy+kGAuAJE4eJFfJJali2gkdF5ZenyU1GBLmjQS2ZK4twqnN45c6NJ2LdG3Zz/d9qLD3JdSX8Oy2g1K2pdnP/IzOpbfVNWyErktRZn9IsHsDwnmAhBdgFovXvkFsDUKq5W1+tbqvxd1liK8ku1S1OVFRi2Cnv/OhamW/o3G3VonEqMo8ovs1s5Nq6/R31ItVSxGQatk+0SCuQB4b0bZ9s6Nr2yPErm0XIhrUZ/XPrddEo8onZlGY56/pcipZZyeIEfk9WrnKGrv2Un986K70hxjTTSjm6JS36XzUfJhuK8394jVggRzgSmlWYdEYuYd89KVUX+lSKSUTs3nBnPbeQqzlOKL6njpV+8mIpp/8/pI7ea2S+29dlHfnvDm9VujN+/ce+e/Nu7o7ynaF/1i0KKfPtFzmAvKabdf+8qdfSRCQNuFvBb9eCKazm3lNluilbRuKeqKUrutglMq9wSoFPGlP8MyLyouzQlGx0pimdf3bgjydrmPLXifd/Q5A4PIUot7Js1g0c84/5YDyTUkt5Hc0f0+JKi3qauzg+Sm7thBJO9Ofp4m+dmubDPJp5Ky85fl6IgowlxwogjLi0xaUrWlqKnU3otISmLnjSHqo5Te9SLjyLZHaYxemrg1Qq1F6yXyz8trXzvf3ude24589G6oFFlOD8PLs+j2EgDbzexykpd0+x9LK5BcA+BPACxhEAzfRXKLmT0D4ISk3l0AbkiaXmNmH570ADwkmAtO/t2aQDmSSSldNFvSva1lQ/ulNpEgpGWlCLcUqUZjL/meC3BryrZ2Dmt+ROcpF81IHNOySGjzdpHYe8c1XzldZvhYyZkATu62rwJwCzLBBPA+ANvMbCcAkNwG4HQAVw8rkDwGwGEAbp2su20oJSv2IE+jRZFXKUrL8dKOUVsvjZenMkv+RuVDO94cXK1OrZ9SdJr+bpljzO2l4/Y+m1KaN23j2Y76qqWmc/v57/QzzX/E9FlGSnYtyTuTnwtG6HadmT0BAN3vw5w6RwJ4NNl/rDuWci4GEWX6n/83SN5D8nqSR43g07JRhCma8NKykWhGZaUoMO+rVC+KYDwbniAObea2PaEsjT+14/kdHasJRy1azsW9Fs3n5V6kGY3dE9uon9boVqwanjazpaiQ5E0ADneKPtFo3/uPkP+hnAPgt5P9fwBwtZm9QPJCDKLXX2nsb9msCMHsctnXANgA4GEAG7s8tlf3YAD3A/jqrPLY80qeLrvpF84uXvBqqb7Wi6R3MY/m93LhaOnDEwXvQh/1XRI4L7WZH/d8z8ed+xGlND2b0dygd1OQ7rekpKP0bH6DFNlRCnbW2MTmMM3s1KiM5JMkjzCzJ0geAeCHTrXH8GraFgDehEHqdmjj5wHsa2Z3JX3+KKn/NwA+NZ7347FSUrLDCeJjAGzv9iMuA/DNqXi14EQX4mFZdHGObJVSst68YSQoXh0viozGkh/z2qRpxJbUai217J2rXIxKNx6tqdRa+ju1l5/3SERLn0PtpknMFsOyUrLLYQuATd32JgBfc+psBfBekod0q2jf2x0bci6S+UwA6MR3yBkYBE9TY0VEmGibIAbJdwBYB+AbGKysEhMkXfKfvuhgiHex9YTHq5sfS21GKdxSH56o5G3y+bW8bjS+fNvru5SyraWuo35K587zI/I9La+lVUuRp+dTPi5ALyJYaczoTT+XA7iW5HkAHgFwNgCQXAJwoZmdb2Y7SV4G4I6uzaXDBUAdGwH8amb3IyTPALALwE4Amyc4hr1YKYK5xwQxyb0miEnuA+DTGOSzT5myfwtPlHpsTUPmdSNbwzali3bUrydEXlSa1/V8i8pLY/Da5OfA8zm3FZ2/Wpq3do5K5zSKOtMy78Yjt1W6ARGzwLB7Bo+VdKnTva7TZnYngPOT/S8C+GJg4y3OsY8D+Hh/no7G1ASzhwniiwB83cweJWuLJngBgAsAYP369aO4KQK8+c1oXgsoL0hJ63pzZ/l2ixCUhKGlbe5H6aagJJ61CLAmwLUIPT9Wijoj372+anaizxFQRLmSMcwswpxLpiaYPUwQnwTg3SQvAnAggP1IPmdme813mtmVAK4EgKWlJU2oTIBRIoo8AvGiMc/2sG1uy6tTSyOWhKgk2F6EGtn3hCyK0nJfSpFuVDfyOafVpneDU7opUDQpFo2VkpIdThBfjmCC2Mw+ONwmuRnAkieWYjq0zKuVLrQl0fSixiHR8ZY0aLrt+e8JdRQRRunRlii1NjYv+vTSpbXPoOVY1KcnnqNEpGKlYNhtM3nTz1yyUgSzOkE8S+fE3nhpuOGbgkZJVbYKQC4YUWRaSnmWxDm3GUVbUZ+RkNai11LEPM4NQUl0I99TG7XPTenX1YdSsv2xIgSzdYI4Of4lAF+auGNiLErC483FpXgCFUWjpfZ5Gjjtu5YurYlV3j4/7tVpGaNXVks353Zb6ue+l9KxeTux2pjcc5iLyIoQTDEfRFHnKEJSSgPWUrclofPqerZyu5G9UgQXCU0prZrv11Kh0XlsSePWzvMQfZPI6scA7DZFmH2xUl5cIOaU0kKSNA3YMr83/Bnup3W8aC/qP6UU2XnHU3+jiDdKzbZEaF4fuY2W+dSSKObnPT2vabkQYk8UYYqJUovCSpHnKPNwpVSkl6bN90sRZlQnH0vJZ69tSyq0NaVaikJrNkuRuFjtmOYwe0SCKSZKa5p2SEu6Mq3nXew9UcvreFFsqd/cbillmvsXiawXHUZC15KqzcdRap/vazHPnGKAaZVsb0gwxUwozTXm5SleFJjvRxGSV8f7nUd9nr1SBBb5mNoq3Rh4PpciyZo/0TkQ84/N7vsw5xIJppg6pWgmfzQlmnvMiQQwLQfKz3hG9qJjkSh7tryylki3FCGOOieqKHIxMS366Q0JpljReEIxSso1PZ7aKYlfbT40StHWIsl8PJG4lqJvr70XlUbnTywaeqykTySYYkVRioK2vXPjHgtxauI4JFrQU0tN1hbulCLESHRLKWTPnjc3mUeXaVtFkUJMDgmmWDVEIldaqDNsl9tJ66bHIgEtpWejqDHqNxe90hym1y7qU1Gk8FBKtj8kmGJV0ToHmZd586Gtc4BeX6VjeduakJZWvkZ9loRaiFfRYyV9IsEUq4aWxUIppcU5Xt2auJUeD8m3o7bD47X0a2nRj97AI1ox6LGSPpFgirkhiiDTFGhp3rC2anV4vGXeNJr/bF0pm0eT3tylEHVMKdkekWCKuWDUxS7DiLQljVqb5yylh2vPdGqhjpg0Ssn2h94lKxaWaIGNlxotzTt6xz1RzaNcIcTqQhGmWEgU1YmFwLRKtk8kmEIIMbdolWyfSDCFEGJO0SrZfpFgCiHE3GKAIszekGAKIcQcoznM/tAqWSGEEKIBRZhCCDG3aNFPn0gwhRBirpFg9oUEUwgh5hnNYfaGBFMIIeYWpWT7RIt+hBBirtk95s/4kFxDchvJHd3vQ4J63yD5LMkbs+M/S/L2rv01JPfrjr+223+oK9+wLEdHhGbz/U5Lkk8B+K8Zdb8WwNMz6nsWLNp4gcUb86KNF5jdmN9sZm9cjgFyHyNfM1ZbsxfvMrOl8frlFQB2mtnlJC8BcIiZfcypdwqA1wH4fTP7QHL8WgA3mNlXSP4VgH83s8+RvAjA8WZ2IclzAPy6mf3mOD6Ow9wL5iwheee4f3CrkUUbL7B4Y1608QKre8zkPsYxZ94MLy1HMB8EcLKZPUHyCAC3mNmxQd2TAVw8FEySBPAUgMPNbBfJkwB80szeR3Jrt/1tkvsC+B8Ab7QpCZnmMIUQYm6xrYaX1o7ZeH+Sdyb7V5rZlY1t15nZEwDQieZhI/R7KIBnzWxXt/8YgCO77SMBPNrZ3UXyx139qWQAJJhCCDGnmNnpk7JN8iYAhztFn1iuaeeYNZRNHAnmZGm9G5sXFm28wOKNedHGCyzmmKuY2alRGcknSR6RpGR/OILppwG8geS+XZT5JgCPd2WPATgKwGNdSvb1AHaON4LR0SrZCTJC+mIuWLTxAos35kUbL7CYY+6BLQA2ddubAHyttWE3H3kzgLOc9qndswD8y7TmLwEt+hFCCNEzJA8FcC2A9QAeAXC2me0kuQTgQjM7v6t3K4C3AjgQwI8AnGdmW0m+BcBXAKwB8B0AHzKzF0juD+DLAN6OQWR5jpn9YGrjkmAKIYQQdZSS7ZHWh3W7ugeT/G+SfzFNH/ukZbwkTyD5bZL3kbyH5NSemeoTkqeTfLB7YPoSp3ymD1T3TcN4P0rye91nup3km2fhZ5/UxpzUO4ukddGSWCAkmP1yCYDtZnYMgO3dfsRlAL45Fa8mR8t4fwLgd8zs5wCcDuCzJN8wRR+XDcmfAvCXAN4P4DgA55I8Lqt2HoBnzOxoAJ8B8KnpetkfjeP9DoAlMzsewPUArpiul/3SOGaQPAjARwDcPl0PxUpAgtkvZwK4qtu+CsCveZVIvgPAOgD/PCW/JkV1vGb2H2a2o9t+HIPVcst6e8kMOBHAQ2b2AzN7EYO5lTOzOum5uB7AKd0D2KuR6njN7GYz+0m3exsGKxlXMy2fMTC40b0CwP9N0zmxMpBg9sseD+sC2OthXZL7APg0gD+csm+ToDreFJInAtgPwPen4FufvPKwdEf6IPVedbql8MMHqlcjLeNNOQ/AP03Uo8lTHTPJtwM4ysz2eO+pWBz0HOaI9PCw7kUAvm5mj66GAKSvh5O7Z7G+DGCT2ar7vqGWh6Vn+kB1zzSPheSHACwBeM9EPZo8xTF3N7qfAbB5Wg6JlYcEc0R6eFj3JADv7l4ifCCA/Ug+Z2al+c6Z0cfDySQPBvCPAP7YzG6bkKuTZPiw9JD0Qeq8zkweqO6ZlvGC5KkY3Di9x8xemJJvk6I25oMAvA3ALd2N7uEAtpA8w8zS18eJOUYp2X6pPqxrZh80s/VmtgHAxQD+dqWKZQPV8XZfy/NVDMZ53RR965M7ABzTfeXQfgDOwWDsKTN9oLpnquPt0pN/DeAMMxvlLS4rleKYzezHZrbWzDZ0/1e4S2sAAALkSURBVHdvw2DsEssFQoLZL5cDOI3kDgCndfsguUTy8zP1bDK0jHcjgF8GsJnk3d3PCbNxdzy6OckPA9gK4H4A15rZfSQvJXlGV+0LAA4l+RCAj6K8QnpF0zjeP8UgQ3Jd95nmNxCrisYxiwVHLy4QQgghGlCEKYQQQjQgwRRCCCEakGAKIYQQDUgwhRBCiAYkmEIIIUQDEkwhhBCiAQmmEEII0YAEU6x6SL7cPTz/XZLXkXzdmHb+Ldl+zik/luTD3XchPk7yz506G0g+T/LucXyo+HdAN84XSa7t274QoowEU8wDz5vZCWb2NgAvArhwHCNm9q5KlWcweD/scWb2M2b2B0G975tZ728zMrPnO7t7vddVCDF5JJhi3rgVwNFdpPfd4UGSF5P8ZHf8AZJXkbyH5PXDiNSLKjPWYvANFs+2OpP09/kuAv57kqeS/BbJHSRPbKkzzokQQvSLBFPMDd23hLwfwL2VqscCuNLMjgfwvxh85VoVM/segBsw+EaSJ0ke1+ja0QD+DMDxAN4K4LcA/BIGL9//oxHqCCFmiARTzAMHdHOGdwJ4BIMXoZd41My+1W3/HQbCVKX7ho7fBfCLZrauE9AW/tPM7u2+B/Q+ANu7bzK5F8CGEeoIIWaIvg9TzAPP53OGJHdhzxvC/ZPt/BsHWr+B4FQAV4/xnZ7pd0XuTvZ349X/gy11hBAzRBGmmFeeBHAYyUNJvhbAB5Ky9SRP6rbPBfCvjTYfBPAukgf06KcQYpUgwRRziZm9BOBSALcDuBHAA0nx/QA2kbwHwBoAn2u0uQXAzQBuJ3kHyWiVrBBiDtH3YYqFguQGADd2j6Asx85PA3jAzI6ahP1K3w8DWDKzpyfVhxBibxRhCjEevwdgi3P8ZQCvn+SLCwC8BoO5TSHEFFGEKYQQQjSgCFMIIYRoQIIphBBCNCDBFEIIIRqQYAohhBANSDCFEEKIBiSYQgghRAMSTCGEEKKB/wfxUR3M5wLS6AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from matplotlib import pyplot as plt\n", "fig, ax = plt.subplots(figsize=(6,6))\n", "p.plot2d(fig=fig, ax=ax)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A synthetic interferogram may be generated," ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWwAAAEOCAYAAACtoy3oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de7QdZZnn8e9PkEuPXAJBRCAGBwKCKOhp8IKtImh0HMBuRFCXwYbJMEjjNINLaG3F2KwVbZXWHps2AzSojCDYalpRRC6iYpCD0oRrEhEhhAZiAuJAgCTP/FG1odjsS9XZtevUrv37rLXX2VW79rufk5w858lb70URgZmZ1d/zpjsAMzPLxwnbzGxEOGGbmY0IJ2wzsxHhhG1mNiKcsM3MRoQTtpk1jqTzJD0o6ZYur0vSlyStkHSzpFdlXpsnaXn6mFdd1P05YZtZE50PzO3x+tuBPdLHfOBsAEnbAZ8EDgQOAD4pacZQIy3ACdvMGicirgXW9LjkcOCrkVgCbCtpJ+BtwBURsSYi1gJX0DvxV2rT6Q5g2GbOnBmzZ8+e7jDMGu3GG29cHRE7DNLG3LlzY/Xq1Xk/71ZgXebUoohYVODjdgbuzRyvTM91O18LjU/Ys2fPZnJycrrDMGs0Sb8btI3Vq1fn/rcqaV1ETAzwcepwLnqcrwV3iZjZOFoJ7Jo53gVY1eN8LThhm1ltRESuRwkWAx9IR4u8BngkIu4HLgfeKmlGerPxrem5Wmh8l4iZjY6yVg+V9A3gTcBMSStJRn48P/2MfwYuA94BrAAeAz6YvrZG0qeBG9KmFkREr5uXlXLCNrPGiYhj+rwewIe6vHYecN4w4hqUE7aZ1YbX5+/NfdhmZiPCFbYN1bJlyyr5nDlz5lTyOTY8Jd5QbCwnbBtIVQm5n35xOKGPBifs3twlYmY2IlxhW251qaanolPsrrrrxxV2b66wzcxGhCtse45RrqSL6PZ9uvKePq6we3PCHnPjkpyLcPfJ9HHC7s1dImZmI8IV9hhyVV1c9s/M1bZNFydsM6sFT5zpzwl7TLiqLo+rbZsuterDljRX0p3pTsan9bjuSEkhaZAdJ8ysZipcD3sk1abClrQJ8GXgUJJdH26QtDgibmu7bivgZOD66qMcHa6oq9H+5+yKezDjnIzzqE3CJtlSfkVE3AUg6SKSnY1va7vu08BngVOrDa/+nKSnn7tLbJjq1CXSd7diSfsDu0bE93o1JGm+pElJkw899FD5kZrZULhLpLc6Vdg9dyuW9DzgLODYfg2l290vApiYmGj8364r63pq/b240ray1KnC7rdb8VbAy4FrJN0NvAZY7BuPZs1RVoXdbwCDpLMk3ZQ+lkl6OPPahsxri0v+FgdSpwr7BmAPSbsB9wFHA+9tvRgRjwAzW8eSrgFOjYjJiuOsBVfVo8P92vmU1d2RZwBDRPx15vq/AvbPNPF4ROw3cCBDUJuEHRHrJZ1EsqX8JsB5EXGrpAXAZETU6jfddHGiHm3uJqlE3gEMLceQ7Kpee7VJ2AARcRnJ9vPZc5/ocu2bqojJzKpT0g3FTgMYDux0oaSXALsBV2VObyFpElgPLIyI75QRVBlqlbCtO1fWzeJKe2Az06TasigdbAB9BjC0ORq4NCI2ZM7NiohVkl4KXCVpaUT8poSYB+aEbWa1UaDCXh0R3QYc9BvAkHU08KG2GFalX+9K75XtDzhhW2+uqpvPNySfraQukZ4DGFok7QnMAH6ROTcDeCwinpA0E3g9yUS9WnDCNrNaKGuUSIEBDMcAF8WzP/RlwFckbSQZ9rywfXmM6eSEXUOurMeT+7XLW0skzwCGiDijw/uuA/YtJYghcMKuESdqAydu684J28xqY5zXCcmjTlPTx5qra2s3jj8TXvypN1fYZlYb45yM83DCnmbjWEVZfuPUnz3u1XMeTtjTxInaihiXxO2E3Zv7sM3MRoQT9jRwdW1T1fSfHd907M1dImZWG+OcjPNwwq5Q06sjq0aT+7OdsHtzwq6AE7UNQ9MS97h3d+ThhG1mteGE3ZtvOg6Zq2sbtib9jPmmY29O2GZmI8JdIkPSpKrH6q8p/dnjXD3n4YRtZrXhhN2bE7aZ1cK490/n4YRdMneF2HQa9a4RJ+zenLDNrDacsHvzKJGSLFu2zNW11ca4/zxKmivpTkkrJJ3W4fVjJT0k6ab0cXzmtXmSlqePedVG3psrbDOrjTIqbEmbAF8GDgVWAjdIWtxh9/OLI+KktvduB3wSmAACuDF979qBAyuBK+wSjHMlY/U2aj+bJU2cOQBYERF3RcSTwEXA4TlDeBtwRUSsSZP0FcDcKX9DJXPCNrNayJus04Q9U9Jk5jE/09TOwL2Z45XpuXZ/IelmSZdK2rXge6eFu0TMrDYKdImsjoiJLq+pU9Ntx/8GfCMinpB0AnABcHDO904bJ+wBjNp/N208jdJQv5JGiawEds0c7wKsavuc32cO/w/wmcx739T23mvKCKoM7hIxs9ooqQ/7BmAPSbtJ2gw4GlicvUDSTpnDw4Db0+eXA2+VNEPSDOCt6blacIU9Ba6sbRTVvdKOCDZu3FhGO+slnUSSaDcBzouIWyUtACYjYjFwsqTDgPXAGuDY9L1rJH2aJOkDLIiINQMHVRInbDNrnIi4DLis7dwnMs9PB07v8t7zgPOGGuAUOWGbWW14pmNvterDzjE76RRJt6VDca6U9JKqY3R3iI26Ov8MewOD3mpTYeecnfRrYCIiHpP0P4DPAu+pPlozG4amJeN05mQ/GyPi4Tzt1SZhk5mdBCCpNTvp6YQdEVdnrl8CvL/SCM1sqJqWsEmGE66i8/julk2AWXkaq1PC7jTD6MAe1x8H/GCoEZlZZcoaJVIzt0fE/r0ukPTrvI3VKWHnnmEk6f0ki7O8scvr84H5ALNm5frF1Ved+/3MiqrrEL8GVtivLekaoF43HfvOTgKQdAjwMeCwiHiiU0MRsSgiJiJiYocddhhKsGZWvqbddIyIdWVc01KnCvvp2UnAfSSzk96bvUDS/sBXgLkR8WD1IZqZTU26wNQ+wMuBfYF9eqyH0lFtEnbO2Ul/D7wAuEQSwD0Rcdgw43JXiDVZ3bpGRql6zkPSfwfmAXsDmwPfB24hmSp/ZtH2apOwIdfspEMqD8rMKjFq3R05nU4y9Hg1sBDYkqQYvWcqjdUqYZvZeGvgKJF3RsQt6fN3S5oL/Juk84EvRkShb7hONx1rZdz3xLPxUpef9wbedLyl7fiHJHNOtgN+XrQ9V9hmVhujlIynKh3d9reSvlb0va6wzawWCm4RNhIk/arbaxGxrN817Vxhm5kNz8sk3dzjdQHb5G3MCdvMaqOBNx33ynHNhryNOWF3UIebL2bTYdmyZdM6JnuUujvyiIjfldmeE7aZ1UbTEnbZ+ibsstdzNTPrpMwbiul45y+SzJo+JyIWtr1+CnA8yZ6ODwF/2aqGJW0AlqaXDn02dRF5KuxS13M1M+umjIRdwmYoj0fEfgMH8uyYBLwPeGlELJA0C3hRRPyySDt5Enap67mamXVTUoVdx81Q/gnYCBwMLAAeBb4F/GmRRvKMwy51PVczs24KjMOeKWky85ifaabTZig79/jY9s1QtkjbXCLpiJK+tQMj4kPAOoCIWAtsVrSRvhV22eu51plHh5hN3wp+BXecWd1jadJBN0OZFRGrJL0UuErS0oj4Td7Aungq7aqJ9HN3IKm4C8k9SkTSBMnGAS9J3ycgIuIVRT/UzKyTkrpEim6G8sbsZigRsSr9epeka4D9gUET9peAbwMvlHQmcCTw8aKNFBnWdyHwEZK7p40b3W5m06+khD3lzVAkzQAei4gnJM0EXk9yQ3IgEXGhpBuBt5AUu0dExO1F2ymSsB9KNxEwMxuKMhL2gJuhvAz4iqSNJPf4FraNLhkkrjuAOwZpo0jC/qSkc4Argex/H/51kADMzCBJ1hs25J6l3a+tKW2GEhHXkWzfVSpJ7wZ+GBGPSvo48Crg7yIi98JPUCxhf5BkXvzzeaZLJICRT9i+2Wj2XNNx87HBMx3/NiIukXQQ8Dbgc8DZwIFFGimSsF8ZEaX/5jEza2lwwm791+G/AGdHxHclnVG0kSLrYS+RtHfRDzAzy6tp62Fn3CfpK8BRwGWSNmcK+xEUqbAPAuZJ+i1JH7aH9ZlZaQqOwx41RwFzgc9FxMOSdiIZdVdIkYQ9t2jjZmZFjGj13Fe6ZsnVJMMN/yw9XXjCYe6EXfa6rmZm7ZpaYUs6HvgwySSem4DXAL8gWVskN890NLNaaHiXyIdJFnpaEhFvlrQX8KmijXimo5nVRlO7RIB1EbFOEpI2j4g7JO1ZtJGxn+noMdhmvVW5bViDE/ZKSdsC3wGukLSWDuub9OOZjmZWG01N2BHxrvTpGenNx2149pKuuXimo5nVwgiPse4rHXf9F8Bsnsm7+5FsZpCbZzqaWW2UtZZIDX0XeAS4kUwPRVFFEvYSSXuXtXKVmVm7plbYwC4RMfBcFs90NLNaaHKXCHCdpH0jYmn/S7vzTEczq42mjcOWtJTkXt+mwAcl3cUABa9nOppZbTSwwn5nmY31TdiSfhURrxr0GjOzXho60/EB4ARgd5JJh+dGxPqpNpanwn6ZpJt7vC6SMYVmZgNpYIV9AfAU8FPg7cDeJNPUpyRPwt4rxzWljMWRNBf4Isk+bOdExMK21zcHvgq8Gvg98J6IuLuMzzaz6VdWhT1ILpF0OnAcSV47OSIuHyCUvVvDoSWdC/xygLb6J+yq+q4lbQJ8GTiUZJv6GyQtbhtGeBywNiJ2l3Q08BngPVXEZ2bDVdYokUFySbpJy9HAPsCLgR9LmhMRUy1Kn2o9STcHnmIzicI7HgzRAcCKiLgrIp4ELgIOb7vmcJL/YgBcCrxFg/4JmFltlLTjzCC55HDgooh4IiJ+C6xI25uqV0r6Q/p4FHhF67mkPxRtrMiwvmHbGbg3c7yS525Q+fQ16W+rR4DtgdXZiyTNB+YDzJo1a1jxmlnJClTYMyVNZo4XRcSi9PkguWRnYEnbe3fOG1S7iNhkqu/tpMh62FsCRwKzSPqFniUiCs2J7/QRHc61/+3luYb0L24RwMTEROPuYpg1UUQUmZq+OiImurw2SC7JlWOmS5EKezFJJfsrhrMe9kpg18zxLjx3+cHWNSslbUoyOmXNEGIxs2lQ0iiRQXJJnvdOmyIJe/uIOHRokcANJPud7QbcR9Lx/962axYD80i21jkSuCoaOA7IbFyV9M95yrlE0mLg/0r6AslNxz0YcGRHmYok7GvLmAvfTdqPdBJwOUmXy3kRcaukBcBkunnCucDXJK0g+W149DBiMbPpUUbCHiSXpNd9E7gNWA98aIARIqUrkrAPBU4c5uJPEXEZcFnbuU9knq8D3l3W55lZvZT1H+ZBcklEnAmcWUYckk7p9XpEfKFIe0US9juKNGxmVkRDV+vbKv26J8kmvK1tFv8rcG3Rxrz4k5nVRtPWEomITwFI+hHwqoh4ND0+A7ikaHt1GodtZmOugRV2yyzgyczxkyTbhRXihG1mtdDQ1fpavgb8UtK3ScZ1v4tkLZNCnLDNrDaaWmFHxJmSfgC8IT31wYj4ddF28qyH/SidZ/q0RolsXfRDzcw6aWrCTtcp2RvYJiIWSJol6YCIKDTGO89qfVv1u8bMrAwN7hL5J5IZ4gcDC4BHgW+RjBzJrU6r9U2LOXPmTHcIZrVW1b+RvCv1jWgVfmBEfAhYBxARa4HNijaSp0vkZxFxUKZrJLs4irtEzKw0I5qM83gqXac7ACTtwBTWZMrTJXJQ+tVdI2Y2VA1O2F8Cvg3sKOlMkvVLPl60kSLLq24BnAgcRPJb4qfAP6dTPM3MBtbUhB0RF0q6EXhLeuqIiLi9aDtFhvV9laSj/B/T42NIxhZ6bQ8zK0VTE3Za8L6DZFjfRmAzSb8tWvAWSdh7RsQrM8dXS/r3Ih9mZtZNwQ0MRk2r4P1SejylgrdIwv61pNdExBIASQcCPy/yYWZmvTS1wqakgrdIwj4Q+ICke9LjWcDtkpZS8jKrZjaeGpywSyl4iyTsuUUbHxWtcabLli2b5kjM6qPqOQojPMY6j1IK3iIJe16nkyVsvmtmBjS6wi6l4C2SsP9f5vkWwDuBwsNSzMy6aXDCfoBnD4v+GXD20EaJRMTns8eSPsczuyeYmQ2sirVEJG0HXEyyHvXdwFHpVPHsNfsBZwNbAxuAMyPi4vS184E3Ao+klx8bETf1+dhShkUPsrzqnwAvHeD9ZmZPq7AP+zTgyohYKOm09Pijbdc8BnwgIpZLejFwo6TLI+Lh9PWPRMSlBT6z2lEirc7x9HATYAeSVacawzcfzaZ3QbSKEvbhwJvS5xcA19CWsCNiWeb5KkkPkuS8h5maykeJvDPzfD3wQESsL/qBZmbdFOgSmSlpMnO8KCIW5XzvjhFxP0BE3C/phb0ulnQAycp6v8mcPlPSJ4ArgdMi4ok+n1nNKJF0SuUJwO7AUuBcJ2ozG4YCFfbqiJjo9qKkHwMv6vDSx4rEI2knkr7meRHR+m1yOvAfJEl8EUl13q+3obJRIhcAT5Es9vR2kl0TPlzGh5uZtZS5p2NEHNLtNUkPSNopra53Ah7sct3WwPeBj7e6MtK270+fPiHpX4BTc8Tzu0LfQBd5EvbeEbEvgKRzgUJb2piZ5VVRH/ZiknklC9Ov322/QNJmJMuhfjUiLml7rZXsBRwB3JLnQyXNAPYgGRYNQERcWyTwPAn7qUzj65MYzczKV1HCXgh8U9JxwD2kQ+skTQAnRMTxwFHAnwHbSzo2fV9r+N6F6QYEAm4i6TLuSdLxJD0Tu6TveQ3wC5Itw3LLk7BfKekPrc8FtkyPG7sJ75w5czxSxMbSdG+ZV0XCjojf88y61Nnzk8Dx6fOvA1/v8v5CSTb1YZL9G5dExJsl7QV8qmgjeXac2WQKwZmZFdLwtUTWRcQ6SUjaPCLukLRn0UYGmThjZlaqBifslZK2Bb4DXCFpLbCqaCNO2GZWG1VMTa9aenPy5HSW5BmSrga2AX5YtC0nbDOrjSZW2BERkr4DvDo9/slU23LC7sLT1G2cTPfNRmh8H/YSSX8aETcM0ogTtpnVRoMT9puBEyTdTbJUdWuUXaGdupywzaw2mpawJf1PkkWeDiMzp2WqapGwB12fdpjcNWJNVoeukKwG3nTcBfgisBdwM3AdSQL/RUSsKdrY88qNbcpa69PuQbr6VYdrWuvT7kOykMo/pMNkzKwBWn3YeR6jIiJOjYjXkSxE9TfAGuAvgVsk3Va0vbok7MNJFpki/XpE+wURsSwilqfPV5Es2LJDZRGa2dA1LWFnbEnSO7BN+lgFXF+0kVp0iVDO+rRD5a4Ra5K6dYW0jGgy7krSImAfku3BrifpEvlCe5dvXpUl7CGvT9t+zXxgPsCsWbMKRmpm06VpCZtko4LNgeXAfcBKpr5rTXUJe5jr03b4rEUkC4szMTHRuJ8As6ZqWsKOiLnpTMd9gNcB/wt4uaQ1JDceP1mkvbr0YbfWp4UprE9rZqOvtYFBnscoicQtwGXAD0hGifxnprARTF0S9kLgUEnLgUPTYyRNSDonvaa1Pu2xkm5KH/tVHWhd+/7M8qrzz3DTbjpKOlnSRZLuBa4l2Rv3TuDPge2KtleLm46Drk9rZs0wSsk4p9nApcBfZ7YWm7JaJGwzszL3dKyLiDilzPacsKfAQ/xsFNW5K6SlgRV2qerSh21mVkkftqTtJF0haXn6dUaX6zZk7pctzpzfTdL16fsvTgdEVMIJewBz5swZiarFxtso/ZxWNEokz1IYAI9HxH7p47DM+c8AZ6XvXwscN2hAeTlhm1ktVLiWSN+lMLpJx1QfTHIjsfD7B+WEbWa1USBhz5Q0mXnML/Axz1oKA+i2FMYWadtLJLWS8vbAwxGxPj1eCew8hW91SnzTsQRz5szxDUirpVHpCmkpUD2vjoiJbi+WtBTGrIhYJemlwFWSlgJ/6HBdZXdKnbDNrDbKGiVSxlIY6aqgRMRdkq4B9ge+BWwradO0yt6FKex+PlXuEinJKN3YseYbxZ/HiGDDhg25HgPKsxTGDEmbp89nAq8HbovkN8rVwJG93j8sTthmVhsVjRLJsxTGy4BJSf9OkqAXRkRrw4GPAqdIWkHSp33uoAHl5S6RknlSjU2nUauqs6paJyTnUhjXAft2ef9dwAHDjLEbJ2wzq42mTU0vmxO2mdWGp6b35oQ9JO4asSqNcldIlhN2b07YZlYLrVEi1p1HiQxZUyofq68m/Yw1cceZMrnCNrNaGLXdZKaDE3YF3J9tw9CkyrrFCbs3J+wKOXFbGZqYqFvGubsjDydsM6sNV9i9+abjNGhyhWTD1eSfnQrXEhlZrrDNrDZcYffmhD1N3J9tRTS5sm5p4q7pZXPCnmZO3NbLOCTqLFfYvTlhm1ltuMLuzTcda2LcKinrb9x+JnzTsT9X2GZWG+4S6c0Ju0bcn20wfpV1lhN2b07YNeTEPZ7GOVGDR4nk4T5sM6uN1gJQ/R6DkLSdpCskLU+/zuhwzZsl3ZR5rJN0RPra+ZJ+m3ltv4ECKsAVdo1lKy5X28007lV1u4oq7NOAKyNioaTT0uOPZi+IiKuB/SBJ8MAK4EeZSz4SEZdWEWyWK2wzq4VWl0gF62EfDlyQPr8AOKLP9UcCP4iIxwb94EE5YY+IOXPmuBprEP99dlZRwt4xIu4HSL++sM/1RwPfaDt3pqSbJZ0lafNBA8rLXSJmVhsF+qdnSprMHC+KiEWtA0k/Bl7U4X0fKxKPpJ2AfYHLM6dPB/4D2AxYRNKdsqBIu1PlhD1iPIJktLmq7q1Awl4dERM92jmk22uSHpC0U0TcnybkB3t8zlHAtyPiqUzb96dPn5D0L8CpeYMeVC0SdtqpfzEwG7gbOCoi1na5dmvgdpI/xJOqirFufENydDhJ51PhFmGLgXnAwvTrd3tcewxJRf20TLIXSf/3LcMKtF1d+rBbd233AK5Mj7v5NPCTSqIys0pVMayPJFEfKmk5cGh6jKQJSee0LpI0G9iV5+abCyUtBZYCM4G/GzSgvGpRYZPctX1T+vwC4BrahtkASHo1sCPwQ6Drf4fGjbtJ6smVdXFVrBMSEb8H3tLh/CRwfOb4bmDnDtcdPMz4eqlLhd33rq2k5wGfBz7SrzFJ8yVNSpp86KGHSg/WzMqXt7oe5+nrlVXYJdy1PRG4LCLuTbqOukvvFi8CmJiYGJu/XfdrTz9X1YMZ52ScR2UJu4S7tq8F3iDpROAFwGaS/hgRvfq7x1Z74nACHw4n6HI5YfdWlz7svndtI+J9reeSjgUmnKzNmsUJu7e6JOyFwDclHQfcA7wbkru2wAkRcXyvN1t/7i4pj6vq4WhtYGDd1SJh571rmzl/PnD+0AMzs0q5wu6tFgnbquVquzhX1dVwwu7NCdvMasMJuzcn7DHXqXIc96rb1fT0GPcx1nk4YdtzdEtYTUvkTsz14y3CenPCNrPacMLuzQnbchvl7hNX0/XnLpH+nLDNrDacsHtzwraB9Ktcq6rAXUE3gxN2b07YNlROpFaE+7B7c8I2s1pwH3Z/TthmVhuusHtzwjaz2nCF3ZsTtpnVhhN2b3XZIszMxlxEsHHjxlyPQUh6t6RbJW1Ml3Dudt1cSXdKWiHptMz53SRdL2m5pIslbTZQQAU4YZtZbVS0p+MtwJ8D13a7QNImwJeBtwN7A8dI2jt9+TPAWRGxB7AWOG7QgPJywjaz2qiiwo6I2yPizj6XHQCsiIi7IuJJ4CLgcCUbyh4MXJpedwFwxEABFdD4Puwbb7xxtaTfDaHpmcDqIbQ7LKMU7yjFCqMV77BifUkJbVxOEl8eW0iazBwvSjffLsvOwL2Z45XAgcD2wMMRsT5zfucSP7enxifsiNhhGO1KmoyIrv1fdTNK8Y5SrDBa8dY51oiYW1Zbkn4MvKjDSx+LiOfsGdupiQ7nosf5SjQ+YZvZ+ImIQwZsYiWwa+Z4F2AVyf9OtpW0aVplt85Xwn3YZmbPdQOwRzoiZDPgaGBxJHc8rwaOTK+bB+Sp2EvhhD11ZfaXVWGU4h2lWGG04h2lWIdC0rskrQReC3xf0uXp+RdLugwgrZ5PIulXvx34ZkTcmjbxUeAUSStI+rTPrSx2D1Q3MxsNrrDNzEaEE7aZ2Yhwws5J0naSrkino14haUaPa7eWdJ+k/11ljG0x9I1X0n6SfpFO071Z0nsqjrHj1N/M65unU39XpFOBZ1cZX1ss/WI9RdJt6Z/jlZLKGJc8Zf3izVx3pKToNUXb6sMJO7/TgCvT6ahXpsfdfBr4SSVRdZcn3seAD0TEPsBc4B8kbVtFcH2m/rYcB6yNiN2Bs0imBFcuZ6y/BiYi4hUks+A+W22Uz8gZL5K2Ak4Grq82QpsqJ+z8DieZhgo9pqNKejWwI/CjiuLqpm+8EbEsIpanz1cBDwJDmWjUQcepv23XZL+HS4G3pFODq9Y31oi4OiIeSw+XkIzPnS55/mwhKSw+C6yrMjibOifs/HaMiPsB0q8vbL9A0vOAzwMfqTi2TvrGmyXpAGAz4DcVxAadp/62T/F9+pp0mNUjJMOoqpYn1qzjgB8MNaLe+sYraX9g14j4XpWB2WA80zGj13TWnE2cCFwWEfdWUQiWEG+rnZ2ArwHzIqKqLT/yTPGd1mnAGbnjkPR+YAJ441Aj6q1nvGlhcRZwbFUBWTmcsDN6TWeV9ICknSLi/jTBPdjhstcCb5B0IvACYDNJf4yIXv3d0xkvkrYGvg98PCKWDCPOLrpN/e10zUpJmwLbAGuqCa9jHC0dpyNLOoTkl+UbI+KJimLrpF+8WwEvB65JC4sXAYslHRYR2QWVrGbcJZLfYpJpqNBlOmpEvC8iZkXEbOBU4KvDStY59I03nXL7bZI4L6kwNugy9bftmuz3cCRwVUzPTK++saZdDF8BDouIjr8cK9Qz3oh4JCJmRsTs9Gd1CUncTtY154Sd30LgUEnLgUPTYyRNSDpnWiPrLE+8RwF/Bhwr6ab0sV8VwXWb+itpgaTD0svOBbZPpwCfQu+ROdMd69+T/Ec5HrcAAAJfSURBVK/qkvTPsf2XT2VyxmsjyFPTzcxGhCtsM7MR4YRtZjYinLDNzEaEE7aZ2YhwwjYzGxFO2GZmI8IJ28xsRDhhW1eSNqSTQG6RdImkP5liO9dlnv+xw+t7Sro7XZd5laR/7HDNbEmPS7ppKjH0iW/L9Pt8UtLMsts3K4sTtvXyeETsFxEvB54ETphKIxHxuj6XrCVZI2TviHhxRPxVl+t+ExGlz8SMiMfTdp+zPohZnThhW14/BXZPK91bWiclnSrpjPT8HZIuSHddubRVkXeqqtvMJFlN7uG8wWQ+75z0fwAXSjpE0s+V7LJzQJ5rpvIHYTZdnLCtr3SlvLcDS/tcuiewKN115Q8ky832FRG3Af9KsirfA512R+lid+CLwCuAvYD3AgeRLLz1NwWuMRsJTtjWy5Zpn/EkcA/JYky93BsRP0+ff50kMfaVrnT3QeD1EbFjmsDz+G1ELE3X8L6VZEu0IPnFMrvANWYjwethWy+Pt/cZS1rPs3/Rb5F53r6SWN6VxQ4BvjGF9biza05vzBxv5Jmf7TzXmI0EV9hW1APACyVtL2lz4J2Z12ZJem36/BjgZznbvBN4naQtS4zTrHGcsK2QiHgKWECy0/b3gDsyL98OzJN0M7AdcHbONhcDVwPXS7pBUrdRImZjzethWykkzQa+lw4BHKSd/wTcERG7tp0vpf0+n303MBERq4f1GWaDcIVtdfPfeO5WYQAbgG2GOXEGeD5J37ZZLbnCNjMbEa6wzcxGhBO2mdmIcMI2MxsRTthmZiPCCdvMbEQ4YZuZjQgnbDOzEfH/AfCmglcZ6gOFAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# this one is empty because there is no phase error.\n", "p.interferogram(passes=2, visibility=0.5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Pupils also support addition and subtraction with other pupil objects," ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "p2 = p + p - p" ] } ], "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 }