{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Scalar field Z" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2021-09-28T10:57:56.182717Z", "start_time": "2021-09-28T10:57:56.168939Z" } }, "outputs": [], "source": [ "from diffractio import degrees, mm, um, nm\n", "from diffractio import np, plt, sp\n", "\n", "from diffractio.scalar_fields_X import Scalar_field_X\n", "from diffractio.scalar_masks_X import Scalar_mask_X\n", "from diffractio.scalar_sources_X import Scalar_source_X\n", "\n", "from diffractio.scalar_fields_XZ import Scalar_field_XZ\n", "from diffractio.scalar_fields_Z import Scalar_field_Z\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Definition of variables" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "radius = 3*mm\n", "\n", "xin = np.linspace(-radius, radius, 4096)\n", "xout = np.linspace(-150*um, 150*um, 128)\n", "\n", "wavelength = 550 * nm\n", "\n", "focal = 250*mm\n", "z = focal\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "t0 = Scalar_mask_X(xin, wavelength)\n", "t0.lens(x0=0, focal=focal, radius=radius)\n", "\n", "u0 = Scalar_source_X(xin, wavelength)\n", "u0.plane_wave(A=1)\n", "\n", "u1 = t0 * u0\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABdnklEQVR4nO29eZhdVZ0u/K46VUlVZhISpgQZEmgCAoGIIDwNKCgOUZSrrY1AX/Wh7b5+SA92i4B9uZLP9hPwE1uv4vBdVHBAAbngpWUUNYIiBgiEJCSBJJCBJCSpTFU5dX7fH+esk1Wr1rzWPlOt93nqqbPXXns6Z+/17vc3LUZEyMjIyMjIcEFXs08gIyMjI6N9kEkjIyMjI8MZmTQyMjIyMpyRSSMjIyMjwxmZNDIyMjIynNHd7BPwAWOM5s6di76+Pu9t9+zZg+effx7Tp0/HoYceig0bNmD69OkYO3as134GBwfx7LPPore3F8cffzyef/557NmzByeddBK2b9+Ol156CUcddRQOOOAA7Nu3D11dXSiVSt7nK6JSqeCZZ57BuHHjcMwxxzhvVy6XMTQ0NOwa//SnPwEA5s2bh66u6jvD7t27sWfPHkybNq3eb/PmzXj55ZcBAKeeemq9vb+/H8uXL8eECRNw7LHHDjvexo0bsXnzZhx77LHo7va/tXbu3IlKpYJJkyYZ+61evRpbt27FgQceiDe84Q1ex9ixYwdWrFiBgw8+GIcddpj3OTYK/Hc6+eSTve8fvu1f/MVfYPz48dp+RIStW7di4sSJGDNmjPc57tixA2vWrMFRRx2FcePGKc9BvHeGhoawePFiZfuWLVswderU+n2zYcMGvPLKK2CM4ZRTTqn33b17N/r6+sAYczrHffv21Z+d4447zvsaZQwMDKCnpwddXV1YtmwZdu7ciRNOOAEbN27Ea6+9hnHjxuHwww/HCy+8MOI6XbF9+3bs3bsXBx10kPJ79AF/VgBsJqLpQTuRQURt8weApk6dSiG45557CAC9853vpG984xsEgP7pn/7Jez/PP/88AaCjjz6aiIhmzZpFAGjlypX02c9+lgDQF77wBdq9ezdNmjSJ5s2bp9zPvffeSzfffLNyXaVSoW9961v0pz/9iYiIVq1aRQBo5syZXud6xhln0NixY2nbtm31/QIgANTf31/vx9t+85vf1Ns+85nP1NtFXH755cp2IqKzzjqLANDDDz88Yt3WrVvphhtuoA0bNmjPl+93z549xuvi/SZMmGDsp8KXv/xlAkCXXXaZ97aNBL/GRYsWBW974403GvvdfPPNBIAmT56s7fP444/Td77zHeNxZs2aNax9YGCgvm5oaKjevnnzZuW9c+GFFxIAOvfcc+ttf/VXfzWi74033kgA6KMf/ajxukR8//vfH7Gf22+/nR599FFl/61bt9LVV19NK1asGLFuyZIlBICOPPLIYdf/la98hebMmVNfvuuuu7TPiAv4tn/4wx/qnyuVStS+ADxJicbhppmnGGO9jLE/MMaeZow9xxi7zmW7Gmt6o1wu1z/v3r0bQFU1iHj66adx0UUXYeXKldr9DA0Nafcrft62bRt27NiBtWvX1tctXryYkx+uuOIKXHHFFdi0aROICNdccw1+/OMfAwCef/55/O3f/i2uvPJKANW3JRnPPvssHn/88fryypUr8alPfQrr1q2rt61ZswYDAwP170w8d9XbKz9XoPpGpYKuHQD27t2rXXfrrbfin//5n/Gtb31LuV48t0qlot2PCJsiUWHjxo3e2zQTMee7fft24/o//vGP1n6nn346PvGJT2DRokXaPv39/cOWRRUg3i89PT31z+Jv/OCDDwIAHnnkkXqbSvnwe+eHP/zhsG0/85nP1O+fcrmMX/ziF9i8efOIYwLVZ+Kv//qvcc455wCovjR/5jOfwV133QUAuPLKK7Fw4UK8+c1vBgDs2rULzz33HADgoYceAlB9exdRLpeHKWtXlX3nnXfi4osv1j5T/BoA92eiEWimT2MAwFuJ6CQAJwO4gDF2eujO7r//fpx99tnDBj4iwvXXX49777132KAkD/wct912G+68807cfffd9e1liMQgL+s+A8CNN96IefPm4Wc/+xmAqikGqD5Ua9aswcKFC3HttdcC2P8Qc3JTkca73/1uvPWtb63fcN/73vfw9a9/HT/5yU/qfWRSlJdliA9qatLYsmULgP3XJGPXrl31z66mh0waZsiDuQyT6UrGihUrtOv27NmjXSfeEzoyURGEqk0mAAA4//zzccMNN9SJ5Bvf+AYuvPBCnHnmmcptXn/99WHLDzzwAG644QZ84AMfAAA888wzAPa/nL7pTW/CCSecgMcee0xLBr6kwceViy66CLfffju+/e1vK/vJYxYfy0TcfffdOOeccxp6XzeNNGrqaWdtsaf255Se/sQTT9TfTjhuv/12PPbYY3jkkUewZMkSLFq0CKtWrcK1116Lq666yok0+KBaqVRw//334+CDD8YjjzwCIsKLL76ISqUSTBpr1qwBALz66qsj1vMHi7eZjsGxadMm7Nmzp/7wqQhGJhsV+YgQH1QdwYSShmlgAfRkYkImDTNSkobp9zPdE7p14v3lShomvwtX2A888AAAYPny5QBGkoa8vGPHDuP6pUuXAgDuvffeYNLo7+/H+vXrAQDXXHMNjjjiCGzbtq2+nr9QqfbL8dvf/hbXXnstFixYgHK5jLvvvhtbtmzB+9//fvz617/G1VdfXe87ODiIH/zgB/VjpkZTo6cYYyXG2GIAmwA8QERPKPpczhh7kjH2JG+76KKL8O53vxvbtm3DF77wBSxZsqT+BRMR3vve9+Jtb3tb/W1hcHBQO6A//PDDOO200/D8888Pa//1r3+NTZs24Xe/+x1++ctfYs6cOfjSl74UTBqmZdvgrhrsXfr4Kg3RYa572E3EELoOGK40XBFCGhs2bPDeppmIIQ15QJQhO69NUJGG6s1fhotiVQWjuCoNDn5v20hCHtBtyxwyMZjWyf2mTZuGQw89FFu3bsXChQuxZs0a3HrrrcO2/973voezzjprmKlQfLkVTVVf/epX8f73vx9nn312vW379u2444478MMf/hBf+cpXcOmll2L+/PnK841FU6OniGgIwMmMsSkA7mKMnUBES6Q+twC4BahGT5VKJbz++usYHBzEvffei89//vNYtmzZiC947969dRMQMFLqcdx111344x//iAceeEDbh9sw16xZUwhpyOtshFCpVEbYOFWEIG/XTPOUTWmEkMbEiRO9t2k3pRFDckUrjd7eXqt6dbmPXAnCpDT4edhIw5dUxP27rpP9hfzcli1bVm+Tx5qPf/zjAICbbrpJ2UfE/fffDwB1Xwvv+6EPfQgAcO655wLYb9VIjZbI0yCibQAeBXCBrW+pVKp/mXyg2bt377CBl6/XkYDLQK/77OoIt5GG+LD59FUtq9qIqK1IoxHmqaGhIbz22mvex2kmWllp9Pb2WrdLaZ4yKQ0dacgDvbieiIzrRZiUxtDQkJNPw2WsEF+exHbRv6oKYjGNPanRzOip6TWFAcZYH4DzALxg204kDd1griINW1/TZ9MP4tpPHsBjzFMupKG6cXxIo9E+jUaYp7Zs2dJSUSguaGWfRgxphDjCUygNnpsEVBV7CqVRLpeHDeQxpKEbg8QgAhVp6FRJEWimeeoQALcyxkqoktdPieheyzYolUr1wck24Lv8GL4/pNgum4lclUalUhn25uBjygLS+TiA4dcm3uy+Po1yuWx8w0lFGuL35mueajfTFFCs0hCTZInIGLWm+v1cEmMbpTRcfRoiVETgqiZEuEZPuYwpLoph1JIGET0DYJ7vdl1dXfWBw0YEqvVElIRMiEhrqpKPIe8jVlm4kIaKIMQ2/h26PNQu7SYFAqQjDdGM5Zt13k6k0d3djXK5jP7+fuuAroNNaYiDz969e42VFlIoDZHwG+XTEO8ROYS+XC47k4zNES5+l7oM/hjLh8081UjSaAmfRihUX7D45m/7MeSB3+VNQNcuE4VIIHLfUB8G35ev0uDbuRKJri/gpkBU+S0+Pg3V9hy2gdAE0alsOkYrYMKECfXPodds2078DmykrVKYOqUh7lenTMX7SzVQu7Zx6EhDJFt5YPVRGibzlKxCRBOY+F2Iz6TLy6qOCDJpRMDmx0hFCK5+i5B1Koe1TDgh/gobQZjaAH9HuC2klpOGbrB2VRo2k4sJ7aQ0RPieNx/AbD4sEbZABNXv22ifhklZ6khD1UdcdvVp+ITc6ogqJhjH5tMo2vktoqNJI5QQTAoklDREteBLEi4qIkSNAHGKwqWdI5V5KpOGHSGhyLbvv10d4XIf2ZcoO8ZTkIa8jiOVeUp1LF3fItAxpMF/BNWPIZNAjKJIrTQAu5IIcYS7qgpfRRGrNHTIpKGH73mHJD3alEaRIbcpkvtClYaIGGLwjZjyjeC0+TREs3yI/8sHHUMatogp38goleowkY/pGD7rgMZliOvaTO1FkYZrnkYmDTtEpeH61tlqSkMchPmA6KI0TCYslyhFk99CHKxlxaJzfqcK9ReRfRqeEH8s1WDuKwF1JqlY85RPXkYRpFGUT4M/wKKsB+ykkaqMSCYNO8SBz5WMm6k0VGQgvi37qIgYpWHLxdD5KmxkY/vsMk6J455KSWTScITqy46xG/qapBplngpREUX5NDjkAcNECkSUzDyVKnqqnVBkrgZHiNKIydOwKQ0R/L5N4dMwLatCcE19xc9FJPS55GyotisabU0aodFTMXbGIkgjlkRUba4E4ZunwSEPGCbSsBEQULx5qlKptF0JEY4is8I5Wk1pqPqmVhqq58pl8JeXY0nDdzxSIZOGI2xKI1W5EBcfiGkbwFxrKoXSSG2esr3Z+JCGTWUAxZuntmzZ0tAHKyVaVWmIA6XuftHdFz5Kg9+jpn4+xKJbdh385eWifBo+pJFDbh0RUjok1nch+kDkH8qVUBrh07ApDVPCnwjdg9pupNGu/gyg2Eq3HCGkIdrWXcjBJyNclWgaqzRcFL4pI1y37OrTiFEd2TwVAZUjHNj/pfqWAkllnhL3JX82bQc0royIajvdQ80RQhryflKShjgI+sSjy6TR6hnh4vk1QmmEJPfpMr915CBClxGu+l1UhCD30/XRZWSrllMpDV1dOd+xRlfdQoVMGo7w8WkU4QhPEXJrqlPlk7hncoS7lhFRQefsDFUausG66DIi7aw0du/ePWxuGBvE78/0fbmUEeGDsI1kY/J9xEg81eDH71FTP36/u/QxLZvKiOjIQC4jIsLFtO071qh+i0wajgglDdsPE1KTyrRfcVmVEW4jBZMy0RFLrE+DIzVp6NAo89TkyZODtm8WeKJeKOnFKg0XZzfgZp4S4ZMf5KKaQ0y2qme0GY7wHD3VQLj6NFJVtjV9tq0zOcJ9o6kakdzH4WqeMkVIFWWe8gGPnJo+fXrQ9s3CQQcdBCCcNGJ9GqbKtyJ8ScOHTFza5JcmsY2jXC6PWC+rh9SkEZMCoPusytNo5DwxbUcatmJgMWxuc5zrKtmqlkMJJWRZ1Rbi01BBJAfxAUutNFxDbmNJY8aMGUHbNwucNDZt2hS0fazSaAXSCH1J8lUatkq2un25kk0qE7nNPJVrTxkQmhHuokZSK40io6cqlYrVfgvEm6fE85QfkhjSKJfLzhVZR5vS4OcbmmMSqzTE399kAvEtZOmTVFqUeSqVIzzEpxFauQJQk0IOuXWE6sv2/WF8HVWxPg3bOpflkAcGiCeNVpjqdbSRxoEHHgggnDRilYZLWK1pXaPNU6Y2l2W56i1HM3waObkvEXxCbmOipxqZER5rnoohCFWb+B2LPo2QQcO2DvCb6tUnikhEu5JGs5WGCBP5d4ojXNffVkbEd7Y+35fSTBqJYIue8v0xYn5gUz95OebmdVn2aVM91OLxxYehKKXhU1Qv1OE3WkkjVZ4GkJY0Yn0aIUrD97nyeQl0USgxL65iey5Y6AlTrLSuTWx3nWDJ1sc2D7j4uVKpGKVuqykN3wcdaIx5KjTcdnBwENu2bUOpVMIBBxwQtI9moV2URgpy0LW7EInKph/i09At+6xLlY+h+5zzNDyhM0+F5GmETM7k4tPwibJSLRcRTQW4lxHRzRtuyh6PyQgPDbd1jRLZvHkzAGDatGn1N8J2yQjn0V6toDTkFwOXucBdyoiI4Peebb563mabP8MWgivC9FyFmptTR3PqkEnDgBSkkSo3Q/6BVTeoTZG4TP8qL/sWJ3TNJOf9QiZaSmme0g3oo80JDoQpjZQZ4SJczVMu5KAjE9O9rHq+YufPECGTjPys26wctnUxFg4bKeToKUeE+jRsn0PmCI9xfNtIwWfZp9SIrq1I85SKEHJinx78nDdv3hykjmKVhnjMZjnCTfdyTCl0n2XZlxby7MeojkaSgg1tRxoqpVGpVJSFAVPlYPgSiGmdSnnEksRo8WmkII2i509OjfHjx6O3txd79+71Ck3m6O/vdyKbctmeKxNCGikc5CYiMZmnVM+NKSPcpxZVLGmkMk/poraKREeQRhGp+jGEIy/HvvEU4dPwbQM6gzTaDYyxKGc4ETl/vyETMXE0K+RWLnPjU9U2VS0qV9JIlQ4gqh552mUgZ4QbkSKhj382RTj5Os5123Okzgh3JQjXfs32aegwGkkDaJ0IqmaF3JraYiddsq0PIY1GJhirUHQdqqaRBmNsFmPsEcbYUsbYc4yxT7tsp1IaMV+6inhC9hkargfE523EEEQrmKeKnup1tJNGqgiqENIYGhpSDnKxPo2iHOHyclFKI5UCaUb4bTOVRhnAPxHRcQBOB/DfGGNzbRupIhhCv3SdL0Tun9oR3goht75tQOuYp1SS3ITRShr8e2qE0jDdG6p1zVAaIc+da00p34mXYvwezY6kahppENF6Inqq9rkfwFIAh/nsIzbM1sTeNlOV6bNtXepkvnZ3hPuapyZOnOjUn2O0kgafi6OZSgMohjRaQWm4zAseQwiuY42MTlYadTDGjgAwD8ATinWXM8aeZIw9CdiVhu2L9q2IKy+7OMhN5xGyXGQZkUb4NGy1pzJpmBFKGvx7SqU0Qn9j14goXXtRPg0XEtHtf9++fU6kUXR2uAodTxqMsQkAfg7gSiIa8UpERLcQ0Xwiml9brq+z+TRsn11+XHnZ5bNMKClr4Lj0B9wT/myk4ZIRPjQ0pMwV4UitNPgbtGuUiIo02iUjHIhXGq6kUZTS4Ovke8bk6zBFQYltohKQHcC2kFqXlzWT0jBNDavaZ4poTvmzCh1NGoyxHlQJ4zYiutNxm/rnWPNUrKnJhRhseRm2jHHVvn0d4b4lQ3x9GraJnFJnhPPB0AVDQ0PYunUrgGoZkXZDTMgtVxo685T8Pfv6NHRlROT9+pQMMakSVea4OB6kePkyzeQXMsuf60um6nNodnjHkgar/trfBbCUiG5y3W78+PH1zz6kofrSXZWGb3icvOxKCrr5vnXL8vZinxj/hS852MxPtozwIs1TW7ZsARFh6tSpxkSwVobvnBr8O06hNGIzwoH4RD7TS1HqvAx5OURpuJqkxO1Vpi7XscZ0/CLQTKVxJoBLALyVMba49vcu20Y+pUNUfgzRye2asu8bNmfaXwpl4ZIRnjqiCvDP+uVIHT3lQxp8mlT+tt5uGeFA4xzhReRpAOknXALc7m+XZV3UE18fQhquJimRKFT7ctmPCuL6CRMmGPuGoGmvXkT0WwDeT7APacT2rZ1nw0NuUzjCGxk9FUsavnkaPqTRrnODi2h3n0aqkiG2foODg1bl4aNMio6eErPZVftyeam15WkUoTqa7gj3hY0IbF+0b3ibrFKKDrlNoURUbT5JVr4+DdOAMTQ0ZK1pVKR5qt0jpwBgypQp6O7uxs6dO60ELcLm05DRyDwNXXsjlUZR5inXaCiRKHxIw0dpZNKAXxSUra/LZ9fqlqYb0OdmVf3IKUJuVX2ANCG3MeG2QCYNGxhjdb8GnxvEBa2uNGLUcMgUsCEKPjVpcPMoETXEp5FJA/YoKBs7+zjK5XbTOtMN6aMc5BtXroml6uPyUPmQRkqlYTNNAW6ksW/fPgwMDKCrqwt9fX3W/hydQBpAmImqkUqjWY5wU1sKpSEO5j6RVRzyeKHzY6gIKFRp+PQNQduRhoiQ6KnQSCvVulTmKZMqCVURLsSia0/p03AhDRdHuKgyfJzZo5k0WkVp6NbFEESI0gghFdPc3yEhtzqTlI/SGLUhtykQ4tx2dVKp2sVl00RNpn2rln1UidxftQyMfKhUD5nPXBxAWJ6GjTSIqD5Y9fb2avuF5GgAo5s0XJUGd8g2MiNc157SOR6rPEzrTY5w03ikUxc2n4YpUMe0XRFoa9LwUQ820ggJv3UNq/WJlmqU0tARBqDP/NYlcMUojYGBARARxowZY8yjGK0lRDiKVBrjxo0DEKc0+O+oWweMHOBcp3aN8WmYQmp9SEJebzJd6V4CZXViC7l1HZNM2xWBtiaNkOgpX0e4/AP4kotqcDYt225k2/Yccqa3y5tZaEa4TBquhALsH6j4wCVvz6EKt3UpBcIdx9yR7LNtMyGfXxFKgx+DJ8zGZIQD/kpURSa+c4SLSGGOciUZm+mKw9WnEao0VMikYUBIGRFfn4ar0rBJWd060dEtD+6uiXsub10pHN4pfBryQCOShslXEerT2LJlC4D2LCEC7E9GDImeEgsWmkiSk0ZMRjigrjEFFJOn0aiMcN32IRYI2TzlozQyaSSCq0/DpbJtrIPbtZ/NdOVSIsRELESkJBsf0mhk9BR/uxXLw6gQYp4iohGk0Y4Z4cB+0uDXYwK/T7q7u9HX12ed8jVUacgoepY+W1tqR7hPwIrrWOKjNFzHGhmZNAzgb+i2L1csHeISkRDi4A6NnvJ9+1FtY9uHqg+QjhxSzKUhmqdUCCGNnTt3YnBwEH19fdb9tzo46fkoDcaYk1/D1adRLpeNA1aKKV9bqYxIaHCL6bOOKFLmaeQ5wh0QGmare3MwqYRYM5ZtXQpHeOwMfY1UGkWSRrubpkTwa3BRGiJcIqhclQYQFnYbG3Lr02YrhW6rNSX31603PdOml1Ef0gjN0ygaHUEath8sxo/henOEkoZv2J+t1LpqGUiTj5FJo3mIJY2dO3dq+7gqDSCMNFpZaYhQPVcu07jK+zL5NHREYavAnEkjMXyc377mKTkmOoRcYqKlVDWkZLhkf3eK0vDJ0+hU0vAxP/AqpybzFM+w37t3r3VASkka7ZDcp1tuhHlKRNGlQXzQEaShm7jEVpvK1xHuWofKJ/Ev9sZWtcU8eKb2GJ+GzgHtShohFW47iTTGjRuH3t5eDA4OOpeSB9yUBmOs/v3bSL5opaFTHzJRNqqMiG45NHpKJAefoAzdeCKPSY1AR5CGDiqfRowj3FV1mEqPyMsxYX+6tlilUYR5SlcvqsjoqU4iDSDMROWiNAC7X0NUIzoUZZ6SX8J0/WIzvn22N/lHTC+pNkUBqMkkK40GIST5T/4c4rfweaNplNJQJfK5KA3+YAwNqUurA3GkkX0a7gghDRelAez//m2kEaM0QjPCVX1TlBFptNIol/UJgeJnFWlkn0aDkKLeVOyNYnNcx779uPTREYQuI9zFbGDK+lYlhBVBGjbbPh9c5Wxwl22bDdX5NUJp6JzhKtJoVEY4MNIp7foc+E665Bpd5TMW6MxQprFFRlYaDUKKhL5Y0gDCk4RcSpAAbsUJY30apqKEPkrDlBHOoRowQzLCO0FpiNeaWmmI37NOafA+PkrDlUxcwmZ1fV1Ki7gElYRu7zMWuMzw55PlbSOYotHRpMGhcxz5+jHkH0tc53pDycrDVjbERiKupUdifRqhpMHXuSgN1zIiruCJcCJptGtGOFCc0mCMBSkNGb6BErp7SjUoupinbIO+am6aUP+i6zo5N0OEDxE0myhEjArS0EUYuJiqxBvPdW6NVElCqn2l8mno2nTtLqQxduzYEev4IKMre54d4e5wLSUivj03wqfBiThFyK2uPURNq54V0/NVVMitTmmYxgEZzYiS0mFUkIYOvmarkHoztuUiHOHNUhoqYkjl0+Aht6M1TwPwVxqMsUJ9Ghz8d08RPQWEE0KsGinKEe5inpItCq2MTBo1uITihvo0fG6+okjDlUh0TvNYpRFDGkRUf1Pmg6ALRjtpAI1RGqlJwyUUV0csMeXTfWpTuVoPymX9DH+t5Nz2wagmDZ1/I4QMXMlFjjmPjSUH3B4o1zY5goSjmUpj165dICL09fVZyy1wDA4Oor+/H11dXZg8ebLTNq2OkKKFraw0YsxT8oDO20zLuv1w+JiKfXwaLvNuZKXRhvANy3XdBojzYdjWq9pizFO+CX+AG2noSMGFNEL8GVu3bgUATJ06ddhD285oBaWhIoZGmKdU4eHy8xji0/DxRaY2VWel0eZohTLpLuuLdoSHzAMeozRcHOExTnBVjka7otXyNDj4756CHHTtKcqG2PrIKtsnIzz1mNHKaCppMMa+xxjbxBhb0szzAPzzNEz1pULzMlTri/JpEJGT7VhO+FM59Yo2T+XIqSpaQWmEmKdC8jdcssdtL0suz4H8/OnMR/Jy6hfIrDTc8b8AXNDkcwCQdupXU15GbIhtiNLwyQi3vS3yAUKVEa5aVyRpmGy/qhwN121bAarzmzJlCrq6urBjxw7lfaBCM30aY8aMAWBXGq7lRXxLi4T4PUL8FvKya5h+KyXs+aCppEFEjwHY2sxz4HD50UN9GqHmKZe5M4C0pdFtPg05QqpSqdS3SRU9JT/ocoVblwS9TlEa4rV2dXXhgAMOALDfZ2ODqDTk79UlI5zDVEZERxrd3d3o6upCpVJRDoopSqabzEmqZD4f0kiVER7Sr5XRbKVhBWPscsbYk4yxJ4s8Tmz0lOsNZboRXcqGuGSEF+nTkIlBbFcN5j5lRHRkkHIujXbOCAf8TVRjxozBmDFjUC6Xtb+tKSPcpYyITBp8G8ZY/X7xyf0JLYPj6/hWLZue41DScPF9ZPNUQhDRLUQ0n4jmF3kcF9no6tMw3Ww2EjBtqysjEhJhBaRRGqJpSh6Qy+VyPeRQFSpLRPU32+zTcEOMM9xl9r7UPg0TafhE6oUm/Pm+hIVaDFK/aLYyWp40mgHfOvmmfj7rgLibnCP13MuAXWnYnOCqN/zBwUFUKhX09PSgp6dHeVwgk4YIl1Iisprj35vJr1FUnkYqpVFUaREf85Tri6GrSbtdkUnDghR+Cx+fRhFlRQB/85QcJeWiNGTkuTTSw0dpcKJ2cYY3Umnw84pJ7gPSO76B4pVGO5mhdGh2yO2PAPwewLGMsXWMsY8383xUSK00fBKIVMsuORkpiMTkuxDRCNLgZpXRnqcBFBd2W1Sehoo0TOoDSOvTiHkJk/0R8vPvYpGwPc/tCLeaDAWBiD7SzOO7IDb81rauCKWRwhE+ZsyYYQOIjkyKrjsFINedElBUgl9RGeH8vhDXjRkzBnv37jUqDVnpFuXTCA1aAdzGBhO5tCuyecqC1GG1scqiaPNUUUpD5dPIpOGPkPpTjVIaruap7u5uMMa0UwjHmKd8+4Sap+RtTduZEgbbEZk0POBaaLCRGeGhBGGrZssTs+T2lOYp17k0fEmjUqnUSWPq1KlO27QLmq00VKShUhMcKtOVLRTXNSM8hXkq1Qtep5ukRGTS8ECI38KWEd4oR3ioT0MuIyJnfauip/i61I5wmTR0Wd3bt29HpVLBxIkTR5Bfu6Mon4ZICqrvNYXSEPfLf5dQVaFqSx09lUlDjUwaHjD9+Kl8GnIuiC/JAGnKiPiap1L6NORBS1YatgQ9F9NUO5YRAdIqDfEYPOxZzPAXoSIVOfFPJg2dopDbXWpShZQRUfXxeZ58/B+uY0MnIJOGB1KQRuqMcNXA76sqxByJInwaKhUCjCQNHRn4Rk91kj9D/k6KUBr8GPx3EP0a/Pfq7u5GqVRCpVIZcX/xciFDQ0MjIo50ZiiV0uAJoKnKiKjaYojAR5WYjtHuyKSRCKFlRGId4akipUxtYrtPGZFmRU+Jc2nI6JQyIlu3bnVWS65FC1WkwcEY05qoGGNaE5WONFTtviarRoTcmvq6Ko1OQyaNRDCF1fnciLFKxKeNP3QiQTQyT6MoR/jrr78OoPOc4ED1++/r60O5XNY6rXUZ4a7l0VV+CyAsV0PnJFcRhIk0QgjBpY+JCGzhsuKyXBixk5FJowCkLCMS4tPwNU+lUhoy+EDR6PnBTUqjE8Cvy1bp1icjHNj/O6UMu22k0nAhG9/naTQ5uF2RSaMAxCQMpcjTMD1g4tSnzVYaLqQxMDCAcrmMnp4e50gorjR4GfFOA78ufp02+CqNlKRhUyA20uD3azPKiMjLmTSqyKRRAFIqjRCfhq8pKoXSKCq5LySxLyuN4Ujh0wDMWeE6M5SPIzzWZBUSlhuawKc73mhAJo0CIL/tFEkariG3vJ8tUipUaRRVRiSk7lRWGsPRTKXhY57ibSrScA2njVUaPiG4oxWZNAqGPO9GI8xTqho5RIRSqTSspk8KpVG0IzwrjZFottJIQRomVeHq51Dd+7HPT2guxmhCJo2CERtimyK5T5WPISb8mdSHa0Z40UrDhzSy0hiORkZPmUjDlr8hEgTvmzoMN0bpZ9KoIpNGwYghiZBkP1WbqCi4b4H34YlZ/HgpHOGxZUTEwcVEGro8BRfSaNeMcMBfaegmYZKPYYue4r8x/13F7VWk4ZIRbiMIFZFwhPgwYpRHJo0qMmkUDJ+J7VNlhMuRUjYnuOikThlymyIjXFV3ypag10nmKdW1+ioNPv/6nj17lFVlXTLCAb3SEJP7fPM0Qnwa4n2dOlrKlreRkUmj4UjpmFMti23c7OQaOSW268hB7m+aI7wVzVPtnhEO+CsNxpjXPOG+GeHAcKXhUkbENeRW1cbv69QZ4ECOkHJBJo0Gw3c+DZu5yiVpz+S70CmK1CG3KhQRPbVnzx7s3bsXY8aMsZYnaVf4Kg0A0aQBpM3T8CWIVGG4Mcl9GVU4kQar4qOMsc/Xlg9njJ1W7Kl1JtpRaZiq38aE3BYRPSWqjE5QFSrYlIbKH6Lza4iIIQ3fPI2ilEZKH4ZqOcNdaXwDwBkA+PSs/QC+XsgZdTgaOXOfbA/2Ca0V23n8e1dXV70SKUfRGeG5hMhIuCoNkTR9lIYtespnytdUyX3ifW7K3UhtnsqkMRKuc4S/mYhOYYz9GQCI6HXGWGfNbtMgNCIjXDZH+ZinVErD5Ox2rT3VqIzwTg+3Bfx9GoCb0mh27SlX9aEyWXH4Ko2QkPbRDlelsY8xVgJAAMAYmw5g9JR1TIiiMlbFQVlWFj7mKVO5dBMxFKU0dLP26TAalMbkyZMBVGcodI3oaZZPw0YaMT4N8VmQQ8k5fENwZRJp9dDsZsCVNG4GcBeAGYyxhQB+C+D/LuysRhFsjm/XNyH+IKlKhqjMU41QGqqQW/G8y+UySqXSMAUkIyuNkSiVSpgyZQoAYNu2bU7bFO3TsJGGS8htjHNc5+fIyiI9nMxTRHQbY+xPAN4GgAG4kIiWFnpmowTim2KlUhmW16F6i9Td9D09PfWKsEA1aY+/fdnKhYjHlElA5wTnb2AmpWEqjS46wU0Oa9/oqdFAGkD1+rZt24bXX3/daYbCRikNkxlKF4rLP5sywm0+jZ6eHgwODkabp3Iuhh2u0VNHA1hNRF8HsATA+YyxKUWe2GiES4y47iEw+S/4g9zT0zOCSGQisCkNeYBXhdwSEYjIaJ5KMWufynTgap5qdbOD7fxSZYWLsJFGSEa4ikwYY84lQ1wzx3VKw5dEMuxwNU/9HMAQY2w2gO8AOBLA7YWd1SiFi5Nb18c3vNYnI1yX8AfoQ27L5TIqlQpKpdKwiCv+4JtIw1ZGxKRMOk1p6K7VN1dDpTR0ZUREJeGbEZ66yi2/d/j9BAwnCJMicVnO0VH+cCWNChGVAXwAwFeJ6B8AHFLcaY1O+Mx5zB8mk9KQnYNF+DR0GeEmsxWgJg3VAJk65LZTcjdSKg1TGRGxj2uehktGuE0t8H0wxrS+OZNzXH4+OHLGdzx8oqc+AuBSAPfW2vTeS0cwxi5gjC1jjL3IGPts7P7aHTqnuFjO3Ja4J69X9dG1AX7RU6akvxSRU4B/9FSnKQ0dUigNGTGTMKVQGiqCAPT3tuzTMPXh9aqyeSoerqTxX1FN7ltIRKsZY0cC+GHMgWshvF8H8E4AcwF8hDE2N2af7Q5dcUP+1mSKjHKJlLIRSaVSqROVyifiaoIC7ErDJRscyMl9OjTTp5GSNGR17RKKq3OEi222ZQ75mcuww4k0iOh5IrqCiH5UW15NRP8eeezTALxIRKuIaBDAjwG8L3KfHQWVarA9FK4EEeLn8Klwa4qcAtyUBhG1fPTUfffdhzPPPBOzZs3Chz70IaxcubIhxzUpDZUT3UVpiMl9qn2knE9DFXKra9fd26YwXFfSyPCHa/TUHMbYzxhjzzPGVvG/yGMfBmCtsLyu1iYf+3LG2JOMsScjj9d2kG94MSSXm6x0D4lIBq4+jRRzaXDIORqyL8GFNAYHB1Eul9HT0zPCjKZDI5XG1772NbznPe/BokWLsG7dOtxxxx049dRTsXjx4sKP7aI0xO/cRWn09PSgp6cHlUpFGYQRojT4fC2iigXiyou4hOHKz4POx5HhD1fz1P8H4H8CKAM4F8D3Afwg8tgqj+SI1xsiuoWI5hPR/MjjtR14zLjqrUnn5E6lNESEVLgNcYTL8DVNEVE92a1opfHoo4/i05/+NABg4cKFeO6557BgwQJs3769/r9IFDV7n0lNhORpAO7lQXTtuns7RGnkPIx4uJJGHxE9BIAR0ctE9N8BvDXy2OsAzBKWZwJ4NXKfHQldpJTYZipO6EokvmaomAq3RZBGf38/hoaGMGHCBGOWeSwGBgbwsY99DESEa6+9Fp/73Ocwd+5c3HHHHTjttNOwbt06XHXVVYUdH2jOPOEhSkNcZ6szBbhFRsX4NDLi4UoaexljXQBWMMY+xRh7P4AZkcf+I4A5jLEja8UPPwzgnsh9dgTkG9z0AIT6NEz5GGKbS1n0opRGaN2polXGt7/9baxevRpz587FtddeW28fO3Ysvvvd76JUKuGWW24p1L9R9DzhKUkjRGmY+rqUEbGRhhy4keEOV9K4EsA4AFcAOBXAJQAuizlwLe/jUwD+E8BSAD8loudi9tkpkG9ofsNzR6KLY9zm4DaF1ro4wk1kQkRJoqdsSkN21jbCCb5v3z588YtfBABcf/31IwajE044AZdccgmGhobw7/8eGyuiRzOURk9PD7q6ulAul0eEh+tqTKnWyXOHqxzk/J4Vw3Dle1alNGRiUT0zquUMd7hGT/2RiHYS0Toi+q9E9AEiejz24ET0SyI6hoiOJqKFsfvrFNhucJWKiImekonEZJ5yTeCLMU/xB10XOaVL0PNxgoeWEbnnnnvw6quv4rjjjsOFF16o7HPVVVeBMYYf/OAHXrPr+ZxfTJ4G37fqGLYEP9W0rrqMcP47uWZ/y+VFOHSOcFONKp9nKMMPrtFTxzDGvs0Y+xVj7GH+V/TJjVa43PCyI9yUER5jshLbQ0JuZUIxlRGRyaCICrexGeHf/OY3AQCf/OQntfs65phjcN5552FgYAC33XZb1PF0xxg/fjx6enrq09va0NPTg7Fjx2JoaGhEf/EYMmnIxGKb1tXXPBXi03BJAJTv4Uwa6eBqnroDwFMArgHwGeEvIwF4tiqHzjylWtbNBx4SKRUbcisOPp04a9+aNWvw4IMPore3F5dccomx7yc+8QkAwHe+851CzoUxFuzXcJmISVVfyrReJAY5YU5luoqZBrZUKqGrqwtEZIwwFGFblp/BDD1cv6kyEf1PIvoDEf2J/xV6ZqMIPiQhL7s4/kx+DlWNqXYNuS3ap/Hzn/8cALBgwQLrMd73vvdh8uTJePrpp7F8+fJCzifUr8H9SSqEZoV3dXVp1UNsyG1MiC2H/IzZljP0MJIGY2wqY2wqgP/NGPtvjLFDeFutPSMBbDd4CGk0Umn4zugHuDnCQ+tOFaU0OGlcdNFF1r5jx47FggULAAB33XVXIefjqzT4d93oOTVShdy6mKxUpligqk5EJZzNVeGwKY0/AXgS1Uipfwbwu9oy/8sIhFiEMKXSMPkvUvs0fBzhIRnhoeapIpTGq6++ikWLFmHs2LF417ve5bTNBz7wAQDFkYZOaeic6EUqDdM6FUHwQbxSqQxLuHMhCMBfachtNpOw+HxmDIeRNIjoSCI6CtWCgv8B4GkAiwF8DcDxhZ9dB0O8gW0kYSIV1RwDQDFKw5VMAHfzlM7nAYTP2leE0rjvvvtARHj729/ufD7veMc70NfXhyeeeALr169Pfk42pSETdQqlkYo05LBbjlClwZfl50GE2CaG8qr6Z3OVHq4+jVsBHIfqXOFfq32+taiTGg0Qb8oYpcFhKiPiSiRFKg0ZfFCKydOQUaTSePDBBwEAF1xwgfM248aNwznnnAMAeOihh5KfE79O13nCi57yVZer4VPpNtanIa83tYnLtmcwYz9cSeNYIvoEET1S+7scwLFFnlinQ7xJfUnCRBquZURcZ+gzJfHFKA0bqQCt4wivVCr1Qf/888/32pb3f+CBB5KeEwBMmTIFgD9pmMxTJlIA0pZHN2V/29SHb7QUYHZ+Z6XhDlfS+DNj7HS+wBh7M6r+jYxAmN5yfJWH2GbzX8i5GzEZ4aZcjEZETzUqI/zPf/4ztmzZgje84Q2YPXu217bnnXcegCpppJ6XnJNGIx3hMaQhZoSL7SaCcE34sz0jsjlK7pOjqdzhShpvBrCIMfYSY+wlAL8HcDZj7FnG2DOFnV0HwectJ0Rp8DbVJEqtmBFuWs8HV130lC7hzafCrc8AzlXC+eef750ceMIJJ+Cggw7C+vXr8fzzzztv53J+oUqDk4ZvRjgwXImIGeHASBKQ22UykVWFSBA2IjHd76plVZvPM5hJZD9cSeMCAEcCOLv2dySAdwF4D4AFxZxae8MWnRFzw4rZ3kD1IZKjPWyqotMywiuVSr0c+aRJk6z9ffDrX/8aAPC2t73Ne1vGWH27xx57LGh7HVKZp1JkhJvW+STyqfwcKp+GfI/Ky11dXSOeQR9HuOkZHO2JgK61p142/RV9ku0In2SimDwNYOQDI/eJ8WmIfX0ywk0z91UqFav5CvCLnurv7wcRYeLEica3Ql+lUKlU8Pvf/x4AcNZZZ3lty3HmmWcCABYtWhS0vQ6+jvAY85ScEe5SmJBDl7/hGymlIhK5j3i+Ps+Vz/M62lXH6KbMAmEb+H2ip0Lq6MTM3JdCacgz94kQCcP01uajNPjAyd++U+GFF17A9u3bMXPmTMycOTNoH295y1sApCeNWPOUCq4+jdApX23tJtIQK+vanglVm495KpOGHpk0EsI17luWzrGOcJcHBnA3T8kkoFMgofNpuORoAK1BGnyg5wN/CN74xjdiwoQJWLVqFTZs2JDq1FouesqHHMR2V6UhIjbEVl6Wzbs+L3mjjUQyaSSE643V3d3tVdIgJgRX19YopaHKCOfrTJFTQOeQRqlUwumnnz5sfyngSxqNKiOSIuTWZIoS+5n6+CgNxpiz+s+kkZEMrqTh68OICcHVtaWc7rUIpTEwMIB9+/ahu7tbSYAyWpk0xO1TkgZ3+O/YsWNYZVlbGZGiQm5DlYareUqEr1/PZTnV89vpyKQRCd0biI9NNJZEbKpC1cbtwyIR+CqNUNJwURqiynBxXhdBGlu3bsWyZcvQ29uLk046KWpfZ5xxBgDg8cej5y6ro1QqYdKkSSAi7NixY8R6+Xtrdu0pXcitiiBU9ahExIbY2vYRSiijgUAyaURCd8P43Fg2EnB5Y5IHiBCTla/S8DFPiSii7lQRpPHnP/8ZAHDyySc7qR0TTjnlFADA4sWLR8w3EQMfE1XR5qmUSkOES1useUpeDjVPjYbyI5k0IuFKGj4x4UX6NGRysfk0dGVEbFnfKZWGC4ogjaeeegrA/gE/BjNmzMDMmTOxa9curFixInp/HD6kUbR5yjarn0vIrfwCogqdBfyVhi0jXF5OQSixs0S2KjJpeMJVtqaKCbftS7Uc2qYqIyKThq28CFDNbTBN0GQrIUJETqQh2u59ScMl45orjXnz5jnt0wZOPpyMTHDNWA8hDW6eUh1DlfEtQjVHOIdrRrhLGRERqhyMopVGCkLp1PLqmTQ8Id8IunkxYvIyUpUVsZmsSqXSsMlpVCTAczy6u7tH5FTo1IQ4OZN8DkSkNV2JfU2koXqDK1JpNIM0OGxvqynMU3IkX3d3N4aGhrBv374RxKAyT/HtY/M0VGQApHOEyyoglU9DHANM8+R0CjJpeML1LaNIR7hMCK42X5NpClBniZv8FrKa4PtPUazQd9Y+V9JwNRns3LkTy5cvR3d3N0444QSnbWwIIQ0bfIoWctLYtWuXUcmoTFQ2YgD8fRqmkFtVP1ObS0a4adkWchuyLpNGBoA0pBEbUhtinnJ58LiqUM0bbqovJa+LKVbI0WyfxjPPPAMiwvHHH6+89hCIpJGq4q2P0iiVSujt7R2m9lQw+TVS5mnwdjEwwEUxqNpUfXxfzjJpuCGThidCbxjZJtro5D5XZ7ncHlKUMGVZdN/oqcmTJzv1tyG1aQoADj30UMyYMQPbtm3DSy+9lGSfjZ6IyUVp+Cb3pWiLzQiXl0NLjmSfRsYIyDdCESG3RZQRcX0YAXeloVunUhKqjPCUSoNXuE2lNFI7wYHqd3DiiScCAJYsWZJkn42uP5UyT0N1T4X6NIoOuc1KYz8yaXiiCPNUI5SGD2m4KA0x5NYn6Q9IMwGTjNTmKT6o80E+Fbh/pEjSMJm+RL+GDiYTVEiehi3kNqRNFdRhC/zweRnzCZnPpNEAMMY+yBh7jjFWYYzNb8Y5hKIIR3is0giNLolRGmJUlay+Gj3VqziXRgrzFBHVJ0w6/vjjo/cngu/vueeeS7I/k9JQOf2bqTTELG+xXYSrTyPEZOWjPEzrQq0NnYJmKY0lAD4AwH9WmiZAF0bHGNNWq42Jnoq9+V32oeqjaucPuo+asJFC6uipnTt3olKpYMKECUke1HXr1mHnzp048MADMX369Oj9ieBKoxGkoUIjfBo6h7euv60tlXPcR3mY1smTnrmSRqdM3tSUqyCipUS0rBnHDkFqNRFbayqkjEio0tC1hWaDA35lRJpR4ZYP6KlVBgDMnTsXALB06dIRb94hCK10G1p/SpxPw3VWPx1ppHSEA/aQWh8S8bEG6MjGto92RctTH2PscsbYk4yxJxt5XPGtwOWtQl4OtZe6rI9xhNvyO3i7LqdDzsXwrTslrjcpDT6o8UFOBT5opSYNbpriA3xKTJo0CYcffjgGBgawcuXK6P01WmmUSqX6vSWW/wDsGeEcunbALfsbSO/4lvuHOs1dx4d2Vh2FnTlj7EHG2BLF3/t89kNEtxDRfCIq3PfhYmpqBGmkcoTbCII/tDJBxCoN06x9RGRVGkRkJA35fENIw+QsLlJpiPs1OcOLKCMCDCcN3TFE0rCVGgHsiX+88oCMUIJw7eerPEITcF3HB13meLvVqCqMNIjoPCI6QfH3i6KOGQuXcgBFmKeKyAgP8XuY2l0q3LpmhOuUhnhNnDRSm6dcHtAilQbg5wy3ne+kSZPAGMOOHTuczF0q85R8DFOpEEA/5asuSgqI81+EEElsRrjP7JqueRqd4jBvX41UAFx+fFNtGZOa8LF1NjPkVvdWWIQj3MWnYTJPcaQ0TxUZOcWRMuy2q6tr2GRMNhRZ6VaXjyGuE6G6/1QDaKrcjRhfYgrS0Jm/2i0JsFkht+9njK0DcAaA+xhj/9mM85Dhqy5iEoDEN7fUPg5VW4yDEXCfNxwIVxoiXHwaHClJ45VXXsGOHTtw4IEHYsaMGdH7U4ErmGXL0sSCNLo8ui5Xg//eqvlCXNQrEF7BVtWWMuDE9VkXxxBZrXRKYcNmRU/dRUQziWgsER1ERO9o5PF1Tijdjx8rP4E4R7hvCK5MSqo+qv0CfuGRgHv0lHg+LkqjKPOUDdxkVJRpCgDmzJkDAFixYkWSCZlCixbqEKs0VHA1TwEj78tSqTTimW0HpeHaT6c0WtXXMSrNUzpm1xFFSPRUypDbIsxVPm2AH2mkiJ5qlnnqhRdeAAAcd9xx0fvSYcqUKZgxYwZ2796NV155JXp/PvWnijRPmWY3dDVPufZNXWsKiLMc6Pr5jimm82kVjBrS8JWJru06M5PPDRh7s4f6L1IojRDzlMskTLxPo0mDz6rH1UBROOaYYwAAy5cvj96XbJ4yRV75kIauEq7OPMUYi07kA8ILFMYmxZqWTU5z1+feZUxRjVOtFp7bWmdTIFx+sCJ9GqY3i0bkbbSi0jCtJyL09vY6OQnbkTSOPfZYAMWQBofKvOFTeypmylcZPkojZP4MIK2ykJdDlIYrudg+t5riGDWkkYooQpxjskMstfnJxeldlE/DV2kMDQ0Zp4IF/ExTQHuSRpFKw4SUjnDXKCldu4/SCAn+iH35cjUpuyYFuoxBqvGo1aKrRg1ppFIUqZ1jLssupGGanwNwzwhX9QPClIbJSd7X16d19LlGTqXOCB8cHMTLL7+Mrq4uHHXUUVH7soGTRooIqph5wlVw9WmYJmKSIedQyG2+fVOUEfF57lzXmSwMvuMO759Jo0mIcXK7vCHI/eS3D58QW3nZFO9tWg4hCJt5yoVMuJJQmZ9SlEUvKiN81apVqFQqOPzww5PN1qeDTWn4zOwXOk+4S0a4CnISn/h7qJI95XZdm4lMbIl7QNqM8NDpX10jKn3Go0waTYKv89v3DQFwf/vwdcCpwmdNA30jHeG6sFp5HYcp3JZvX6R5yhTGyE1TfEAvEkcffTS6urqwevXqETWcRLiEXcaap+RjuJYRUcHHpxHjCLeVIAnJCA81T8VaH2zEkh3hTUIRzu9GmKdUTrAQR3iMcxzwM09xhCoNn8S+lHNpNMqfAVS/zyOOOAKVSgWrVq2K2leI0nAxT9nKiKjg89Lh4wgvwqeRyhHuGuQSqjrEUuytkLsxakjDN1mvaNLQ+SBsTnNVW8rs7xSOcI6QyCrAL7Ev5VwajSQNIJ1fI5Q0dImFrtFTPutU90l3d7fyLbpR0VM+y6HPeqqX0lbyb4x60khlqpKXQ5WGzZSlagsNuY1NuALcBxBOkqmVRjtGTnHMnj0bAKJLpPuQRqlUQl9fH4hIm4fhGj2lgu/LReh92S7JfTF+0lb1b4wa0oiRibFvFaFJRHzZVhKkVZWGanDxmerVhTS4aaodSYNHaDXSPAXYTVQ2n0YqpaFrD03ua/UyIin8qrFqOgVGJWmkNFWJD1VIuRF52SWSypdEdG1F+zRMjnAXpdHIulN79+7F2rVrUSqVcMQRR0TtyxWpScOl9hRgz9Xo6elBd3c3hoaG6nPBiwhRGjGJfPycfPukJJFGKY1snmogbA4iFwkYYsIS5zDQ+SlspBHal6PVyohwmBzhLj4NH/NUrBN81apVICIcccQRyu+4CHDSWL16ddR+Jk6cCMYYdu7ciXK5bA3XTTVPuArNVBq2Pr7Plu5lzhSOK/poTHOJ+7SLxzOF/DcKHUMati8wxgll+qyb+MbHTxFqyvLp08gyIqZ1rVqskA/cRSf1iTjyyCPrx46pdqubU0P3IhUaQcWRKnpK1+5yX6oq3zbLp+FKKDG5Yiql0SwCGTWkEZODYcq/KJfLyuP5mKfk/Zuyu3WkkSoj3CeJD/AnDVPiH0czzFMvvfQSgP0DeSMwYcIETJ8+HQMDA1i/fn3UvrjS4j4e23EBs9Lgv08ReRquSadi3yKmd/UxT6X2d4Q6xV3NYkVi1JNGrM1RRxomIvC5OUOUhoogYogE8CeTUqlkNPOkUBpEFEwa8mDIlUYjSQPQ+zV8MsKB/aSRevY+FWRiCM0IB9wJwqWfL4mYnsPQjHDTuhR+VZd9FI1RQxq+zO76hiCSho4YYqKndDe2rxppdEa4SYEAcRnh4nFS+TQ4aTTKCc5h82u4JnP5KA2xlIjuGD6k4bIu1jxle1mKVRIqEklNGjEBOCrzlG67otExpGFjWl9ycCWQUPOU6zwcIYQQ47/w9WmEmK2AdNFTviG3ukG4GeYpIF0EFfdppDJPyb+P+L0VHT0V6ufwLSPi86yFkkYKohA/28xXjSg50jGkoWLa0B/Gp08oaehCdZvtCE9FDqaBBUiX3JeqhEizzFP8eLGkUZRPQ4WUSiPmZSZ1Brirs1teJ2/nkq8VasLyIZWi0NGkEZOs5/qGEEoaOrNWSp+GbTsfpREbVSW/4adO7uNv2iHYvn07Xn/9dYwbNw7Tp08P3k8IUimNEPOUS/SUCq0QPQUUWzYEaE/zlLi+KNUxakgjJh8jtXmqp6fHaTve17Ts2qcdlYaLeYo7fmOUBjdNHXHEEQ0vCJeaNFI5wk2/X6oqt0A4Qaj6pVQW8rIt4tFluxQWDtuYxhir38NFqY6OIQ0VqxZdpNCVNFQ3qyrTVtfXFgWVsoyIr9LQDRIhjnCORpunmuUEB4CZM2eiu7sb69ev19aCckGrmKdS1J4qSmn4vpCJy/L4ktKn4dPfh2yKiqjqGNJQoegihd3dYcl9rmQDhKkI1zIirnMc9PT0aKVuEY5wTqiNMk81ywkOVO+tN7zhDcPOIwQ+jnD+vZoSCk2/T0qlEWOesgWAyPe8KUPbtixGV7mG48r9QnPFxPP22V9WGhaozAopfjDb59TmqWY7wn1MVkC4ecq2Hmi80mgGaQD7Fc6aNWuC9yEqDdcyIiZ0otKQ22LKjOgsCSZCiTFbqaKjMmkUgJTkoPvMH1CZtGxvIibzlItjXHxLlG+OmKldVUl8JtLwNVtxmN5k+fFtxFIul7F79250dXU5DYQ6NNM8BQCHH344gHSkwWErI2KC6fcZM2aMdt8+c4TL7aa+sT4N1X58VL0rGZjWpSATxlidOGx9O4o0GGNfZoy9wBh7hjF2F2NsShHHiQmzDflsKlgo52W4mqd0N7ZIOqLzS+yTao5wUwkIndkqVmmMHz/e6pTmTt9JkyZFObCbaZ4C0pKGjyPcBNPvwxjTJnX6ZoSr+ocSiU05iG188BXvX505SrXeRBqpiqOG+jc6kjQAPADgBCI6EcByAFcVcZAiEvpc9mnaHkhjntIpFSAu5NZnLme53ScjXPUmK25vehvm/Xg2eIg/gytEImq6eUpFGqFlRHwc4RwqwrWRum8AROpJmFRtLqShe55sc4TLyymVhg+B+Pg3OsoRTkS/IiI+aj4OYGbIfmxxyDE+jRgFIi+rbkaX8iOqbV1Iw0QQouNe9Sbi49MISfLisA1KRYXbygPk1q1bsXPnTkycODHJRE4hMCkNVwUV4gg3wZQRDqRxeAPpfBryoM/Vt+15sr3cufRtpHlK/GwjlU5TGiI+BuD/6FYyxi5njD3JGHtSXid+qao3s5joqRgCcVnnGnKbSmnwtsHBQe12um1DSCMmTwNonBN87dq1AKoDd6NzNDiK8mnoEGueMq1vVu0p1+1i6r41Qmn4vND6mK9Sohj9AoAx9iCAgxWrriaiX9T6XA2gDOA23X6I6BYAt9T6D2MG25cS47uIVRquDrfQOQFClYaNNGLMUyJiMsKBxoXbctKYNWtW8D5iwY+9du1aVCqVoExecT4Nm2JKQRqtoDRcQ83F/Bef59S1r2m7VD7TViopUhhpENF5pvWMscsAvAfA28jXgFtDCtJI5feQzyXGdGVaVkVPyTCZolpBaTDGjIQDNC4bvBVIo6+vD9OnT8drr72GjRs34pBDDvHeR3d3N8aPH49du3YZy4MAdpXn0sdXacRO9xqiIlz6mJ5FV/OU67Pvqy58HN0d7whnjF0A4F8BvJeIRs4n6QjxS7HlaYQ6nkz9Y8xTunWx9XLENhVB6BzwHEUpDRF9fX1WU1CjzVPNJA2gsSaqUqkUHd2WSmmkTO4LIQ0fpaF7bl2zxVNZO0YlaQD4DwATATzAGFvMGPtmyE58lIbrD65q9/1RTdvIyy6lDkyOctPD5Ov30LWndoS7vOm6kAYnwEwaVXATlel357ApudTRU+J9LirgmNpTRSsN12falPmdwo8hZ7I3mzQKM0+ZQESzU+wn1Dyl286lJIjuzcQkUXt6eoY56kMTjFyWxTaTKUq1HZBOaYQWvOPwSdaL8WmsW7cOQGeQhg95jh8/Hq+99pp2fah5SndPiC8+4n2Z0qcR6wh3VRoySfDlcrk8otxIkT6N0aY0kiCF0hDhknDnus9U5ilfUhHbTKTRCUqDI4XSmDkzKOo7GXj9qUaShgmpzVMixPvSZp6KmSNcbOP7CQ1YsfWN8WOE9rEdJ5OGAinyNET4Fh8sgjS6u81VbU3JSRy+5ilV5q5Llm+7k0alUulopWHyGxVFGowxZaUBEQMDA/XPRc4RrmpzNU/p1ITLfor4bFMSmTQcoVMMHC7mKdFsJJKG2O7i3+juHl7xVvZDmEqMiEhpnjKRhkgEOtOZqq+p3XVK0FjzlDwYhWaEv/baaxgcHMTUqVOdiKxINMM8ZYKPeco18Y9DJI2UeRo+feRSJLysiMl/GPIimDqSatRmhKeCLXpKx8RiXzF0VTRPmQhA99lk3hKPI5KQa4FC12XALSdDfDDFfqrvsVHmKdcyIjJCM8JbxQkO6EnDJ+EwlDRcyoj4EIPufnGtQJAqI1zuo1o2vbT5ZI+nipjy+Zx9GgHQKQbVeheTlO6zbp+upCHXmnJ9m+Hb+iyLbSGkoUIKR7irCuFohHmqlUhj+vTpGDt2LDZv3ozdu8Oi0H0Ul+377enpMb6ppppTo9GOcB8fYqjSSB1aq/usUyIqUkmJtiYNm/zS/QA6k5Q4sIc4xU3bhBYoDHGE84eHX6fK9yM+YKK5QIVmOMIbET3VSqTR1dWFQw89FADwyiuvBO0jpXkKaMycGqo28Z7nCt2FEFSDZIxyL9I8FZMr5mqeyqShQGjIrQgdUfhOrhRKGkUqDQ7bW1oRSiOmLhUw+pQGABx22GEA2p80YpWGqEq5X870ciRu5/u8mJ6/IqOnUiiTbJ4KQGjIrXhT6kxSvqRRKpUKURoxPg3dstxWhNIwDTipo6c6QWkArUcaoVO+xmZ/i+AvNKmc441WGjGmKlv4bSaNAKQgDZ26CAm/NZGGLpIpVlmExq+nUBqh5qmUyX3jxo3TJiraMJpJo8gpX32Uhi2vw0QaKaaATRkq7/PZp0yRS0Z4Drl1hE/0lK8jvFHmKZeb2Ce6StUWqzR0b4OhGeEplUa7V7gVEUsaKR3hQHqfhuvAL4KThsrP4aKyi3CEm3I4dJ8B//IiPj6LHHLrCBuT6n4Y0RHeCJ+GHD2l2x/vG7MM+EeVtKLScCWNUH/G0NAQXn31VQD7B+tmo1PMU7E+DRH83lT5OXyUhmsIbuqQ2xizVc4ILwA+pBETcqvLsTCRRsyc4CF5G6YEvRDzVMi8z+I5mFRIyuipUNLYsGEDhoaGMGPGDKfSF41AatKIyQgH4pWGfHyfEjUcKrOuyTmeyjylUhOu/g4f05NpW1u7OBZln4YjfHwa4hesUxqxORspoqdClYaOsFTbAMX4NEQzl6nES6zSEK8v1DzFB+ZWURlA6ykN16z+mKlgdeVsOFT3Jm+LIQ1VRrjcP7VPI0WYrejf4Fnsur7ZPKWAODCpJiXSDVwuJqkQp3goaaTICA+d/hXYT6K6m8yVNPbu3as9BxE6pSEOEK7ObV+lwb/rjRs3AkA9N6IVwM+Fm818MXbsWGfVZMsIB8yKMFWehqt5SoTKPOVCAoA9xN0UcpvSPJWCZHTzgWelYYAuSY9DNwi6KIoQp3hI9JTvTV0EaZjaAPdBwEQaLhnhthnnVIipcAu0FmmMHTsWBx54oPYlxQWuystXafiUEfFRGuL9rbruopSGy3JqpaHbt69/w0QQOiJJibYmDRGqAdsloc8lIzyWNBrhCPedOyMFaciDgzgXswm6N1jX7UXERE8BrUUaQLy5zJVEU5qnZPgoDd08Gxwmn0bIbH4+Po7UeRq+vg6f/qq2TBoK6JL0OGLCbFOSRkzIbVHmqRSTLfkoDRE2n0aj5tIAMmmYEGqe8p0nnMOkKlRtRTjCGxU9FUMsruG3mTQsCCWNGEe4/EM2I7nPtzihbxvQOKXB4VN3KpY0DjnkkKjtU6OVSKMRtadE+JJGiuQ+H2IIIZQi/BuZNBLB5tPwLR3SLKURmhGewqfRKKVhI42sNMLRCuap0Fn9VEmmJpNViNLwcYT7rAtVCC7b6jLCxc+ZNAKQWmnoSEMknyJCbkMd4a3g02iGeSr7NIbD9fsQiVsVeSj3kdGpSiNETZgywlMVKbQpkEwaAfAhDReiCImeMpmxYsxTLiG5vkpDRSS+D3tR5qlGKY2uri7MmDEjePsi0CilIT4but+tEbWnRKTwaegc4SGTMpkS/+SxxXewt+VYuH5WOcLFSCqfSbxc0dak4RNy6zKHhku7bv+ympAH+naLnjJlfuvabTWsOGxKo1E+jYMPPriwt7FQNIo0ROhCnU3k3t3d7Z3XE5OToWor0hEuDuqqviYV4eKXENfZ8i1cPrv0TYWOJg1fpeGbp+FjnnIhHr4ckuzXqDwN8dxCo2RSKo0Y81SrmaaAtKTh+papmynQRu7899flcLiUERHhqjRCak9xiAO1uJ4vm8xOKZ3f8rKtnpTNv+HiO0mFtiYNG3RVcF3IwSV6SryJbNFTJke4iSRSJPcV5dPwKWQoXr8t29tGKuL3FaM0WpE0pk6d6pwNr0JKpWEjDd36UJ+GSqmqnhuV0nDNCI8JcU8ZZqtTKj55GjaFojKjpUBbk4ZPOXTfelM6R7ipeGGKkNvu7m6v4oMplYb4UIvXKR5TbNeRjIo0RLu57Q24SPOUeOxWC7cFqud38MEHB28forxczFOq30ynNFP4NEzE2aiQW3k5NqFPNnnZlEaIecoUaZUKbU0aKphsiBy+YbY6knGtZGtSIaqbWDyGbFdtlHlKfOsTr1PXLkL1BuqT7W0zT4mDi0vxQx1aUWkAcWTWSKXhSw4+pGEyZYU4wkOXUyb0ycu2yrUupMTbO96nwRj7AmPsGcbYYsbYrxhjQU+vTWmI68U35BjzVKgJyscR7mNqCk3us0VP6ZzaLmG1qsFEZzdXwUYa27dvr3+OiQ7JpFFFavNUijwNE2nw56NUKjm/VHHEhLinJg2bzyLUEe7y4hyDZimNLxPRiUR0MoB7AXw+1Y51pCHCRUWEVLw1kYZPwULdflT9Q30atoqjOtJwiZCKVRo289SOHTuc92VCJ5JGUeYpFZqtNIC0eRlAGvNUSBKgj3PbNU8DsH/nIWgKaRCR+NSPB0C6viKmTZuGp556StzPiD46P8bFF1+ME088EQsXLsTRRx9db/c1VU2ZMkV5bkceeWSy+TRkAhBVkvzmUFQZkXZRGjHIpFFFavNUitpTLkpD1S9l7Sl5vWkdf9Hp7e2t+6SmT58elM/R19cHxhgOO+wwnH322bjooouCQ24XLlyIa665BimhDrRuABhjCwFcCmA7gHMN/S4HcDkAHH744Zg3b55xvzo5NmHCBDz99NMAgP7+fhARLrnkEpx//vn1Pi4mKb4fPrAtWrQIv/vd77BgwQJMnToVF1xwAW655ZZhx7aRhkhupVJpBGmYlIdLGRHX6CkX85SL0sikEYdGkUZXVxcqlYr2haOZ0VOtojRM5qnzzjsPjz32GGbOnImPf/zjGBoawmWXXYYZM2bgN7/5DY499tgg5/nNN9+Mbdu2Yfz48Xj00UcBAA8//LCyr02V/M3f/A0A4Prrr0cqFEYajLEHAajCQK4mol8Q0dUArmaMXQXgUwD+TbUfIroFwC0AMH/+fKsicQlXnDhxIm666SYAw8lh5syZ2LRpEyZMmGA0SZ144on1z2eccQbOOOMMAMBZZ52F/v5+MMZARJg1axYqlQq6urqcI6sYYyMIYPLkyTj//PNxwAEHaK/XV2nMnj3b2E9nmnBRGq1snhIHpgMPPDB4P0UihjTE7842L8f48ePR39+vXR9qntK1T5061bg/8R7+6Ec/iuuuuw7vfOc7jf10tdg45s2bh6OPPhrvete7lP19SOT444/HuHHjMH/+fADAd77zHRARGGM45JBD8G//tn8IO+usswCMfJkTB3R+rl1dXcOOc+mll464Zlv4bdHObxGFkQYRnefY9XYA90FDGq6477778Hd/93f46U9/itNPPx1AVR7aIA7mP/3pT/Gv//qv+PznP4+LLrpI2ccG7kdhjGHVqlX1Zfnm+dSnPoWf/OQn+OAHPzji7X3WrFkj9vmrX/1KeTwXn8Zpp51W/7x69Wps2rQJM2fOHNFv7ty59c9vf/vb8clPfhJ/+Zd/OazPm970pmHHVUHljG0VpbF+/fr6Z9OUtM1EDGmI12QiBKBKMKY+ojJQEZBOaUybNg29vb0jSrScccYZuOKKK3DSSScNaz/iiCPw0ksv1QdjALjmmmvwlre8BW95y1vqbR/5yEfwox/9CBdffHG9zWaemjFjBl588cX6Mn9R4GZmWQmceeaZmDJlSv1FUFw/e/ZsvP7669YpakWYzFP/8i//gt/97neYM2eOdaDnY9nUqVOH7WPBggXYuHHjMAtM4VUOiKjhfwDmCJ//LwA/c9nu1FNPJSKik08+mQDQP/zDPxCq/hASwdueeuop+v3vf0/PPfcc6dDX16fcx3vf+14CQNOmTaOPfOQjyj4+uOGGGwgAXXfddfW2SqVCRET79u2jKVOm0OzZs+vrHnroIVq6dKl2f/x8HnjgASIiOuWUU0ac49q1a+nBBx+0nttTTz1FX/rSl6hcLjtdy5NPPkmvvfbaiPYf//jHdOqpp9LatWtHrDv33HOt3yFf/8gjjxiP//d///cEgI466iin8xXxuc99Lvq3LBqvvPJK1DnybV988UVjv9mzZzv/JitXrhyx7uc//zmdcsoptHr16hHrnn32WeV9oMLLL79M119/PW3dutXYb9euXXTPPffQ7t27R1zD8uXLiYjoxhtvJAD06U9/Wruf+++/v/7drF27lsaPH0//+I//WF/Pn0siogsvvJDGjBlDGzdudLoWFRhjBID27t1Lc+bMIQAjvpu5c+cSAFq8eLHxvNeuXUvf+MY3CABde+219XVbt24lADRx4kS6//77iTFG3/72t+vrATxJqcbvVDvyOijwcwBLADwD4H8DOMxlO04a5XKZdu7cSStXriQA9MEPfnDYl/vEE0/Q97//fe2XL2LBggUEgM4555xh7Rs3bqQrr7ySli1bRl/84heTDDRr1qwZdkOKGBwcdB60iYiuu+46esc73lHfZtmyZfT2t7+dnnjiiahzLAp8oB87dqy2z6WXXkqzZ8+mvXv3Gve1Y8cOuummm2jdunXe5/HQQw8RAHrjG9/ovW2jUC6Xo+633/zmN073/6JFi2jy5Ml02223afvw8wj5rhuBr33ta/ThD3+YhoaGiKhKVqeeeio99NBDzvvYs2eP9rkcGBigDRs2RJ3jpEmTaMyYMTQ4OEh33XUXXX/99SP6nH766QSAVqxYYd3ftm3b6Ktf/eqw8yqXy3T00UfTWWedRURE27dvH7ZNStJgpIhAalUwxvoBLGv2eRSIAwFsbvZJFIhOvr5OvjYgX1+741gimphiR02LngrEMiKab+/WnmCMPZmvrz3RydcG5OtrdzDGnky1r9b0BGZkZGRktCQyaWRkZGRkOKPdSOMWe5e2Rr6+9kUnXxuQr6/dkez62soRnpGRkZHRXLSb0sjIyMjIaCIyaWRkZGRkOKNlSMM0xwZj7CrG2IuMsWWMsXcI7acyxp6trbuZ1XL6GWNjGWM/qbU/wRg7ogmXNAyMsS8zxl6oXeNdjLEpwrpOuL4PMsaeY4xVGGPzpXVtf30mMMYuqF3bi4yxzzb7fFzBGPseY2wTY2yJ0DaVMfYAY2xF7f8Bwjqv37GZYIzNYow9whhbWrsvP11r75Tr62WM/YEx9nTt+q6rtRd/famyBGP/AEwSPl8B4Ju1z3MBPA1gLIAjAawEUKqt+wOAMwAwAP8HwDtr7X8vbP9hAD9pget7O4Du2ucvAfhSh13fcQCOBfAogPlCe0dcn+G6S7VrOgrAmNq1zm32eTme+18COAXAEqHt/wHw2drnz8bcp02+tkMAnFL7PBHA8to1dMr1MQATap97ADwB4PRGXF/LKA3Sz7HxPgA/JqIBIloN4EUApzHGDkGVaH5P1Sv/PoALhW1urX3+GYC3NfvtgIh+RUS88uHjAHjFwE65vqVEpMrW74jrM+A0AC8S0SoiGgTwY1TPv+VBRI8B2Co1i9/9rRj+m/j+jk0DEa0noqdqn/sBLAVwGDrn+oiIdtYWe2p/hAZcX8uQBlCdY4MxthbAxdg/m99hANYK3dbV2g6rfZbbh21TG6i3A5hW3Jl742OoMjrQmdcnYrReX7viICJaD1QHXgC8VG3I79gSqJk356H6Nt4x18cYKzHGFgPYBOABImrI9TWUNBhjDzLGlij+3gcARHQ1Ec0CcBuqc2wAVckkgwztpm0Khe36an2uBlBG9RpN59qW16faTNHWktcXiHY61xiE/I5NB2NsAqoFUq+UrBkjuiraWvr6iGiIqlNmz0RVNZxg6J7s+hpae4rC5thYB0CcYGImgFdr7TMV7RC2WccY6wYwGSNleHLYro8xdhmA9wB4W00KAh10fRq0zfUFQnd97YqNjLFDiGh9zXSxqdYe8js2FYyxHlQJ4zYiurPW3DHXx0FE2xhjjwK4AA24vpYxTzHG5giL7wXwQu3zPQA+XIuoORLAHAB/qEmvfsbY6TV796UAfiFsc1nt838B8LAwSDcFjLELAPwrgPcSkTgjUUdcnwGdfn1/BDCHMXYkY2wMqo77e5p8TjEQv/vLMPw38f0dm4bauXwXwFIiuklY1SnXN53VIjAZY30AzkN1zCz++podBcD/YJhjA8DVqHr7l0Hw7AOYX9tmJYD/wP4M914Ad6Dq7PkDgKNa4PpeRNWmuLj2980Ou773o/rWMgBgI4D/7KTrs1z7u1CNzlmJ6nTGTT8nx/P+EYD1APbVfruPo+o7egjAitr/qaG/Y5Ov7SxUzSzPCM/cuzro+k4E8Ofa9S0B8Plae+HXl8uIZGRkZGQ4o2XMUxkZGRkZrY9MGhkZGRkZzsikkZGRkZHhjEwaGRkZGRnOyKSRkZGRkeGMTBoZGRkZGc7IpJGRkZGR4YxMGhkZjmCM9THGfs0YKyXY1xjG2GO1MikZGW2DTBoZGe74GIA7iWgodkdULaP+EIC/ij6rjIwGIpNGRoYAVp3t7fza5+sZYzcLqy+GUJeHMfYoY+zY2udpbPgMeHcwxv6DMfZbxtjLjLGzGGPfZ4wtZ4x9t9bt7to+MzLaBlkaZ2QMx78B+B+MsRmozsHwXqBqTkK1BtZLQt/ZqNb4Aaq1gJ4V1r0RwO+J6FOMsf+BavG8cwBsQbUS6d+jWu/nTcVdSkZGemTSyMgQQESP1ap9/iOAcwRT1IEAtvF+jLE3AHiFiCq1phNRLR4HxlgvgCkA/t/auj0Avku1yXEYY7sBDBIRMcYGGWMTqTq7XEZGyyObpzIyBDDG3ojq/NID0kC+B9Xquxwno0YSNZwqLB8P4CmBUE5CddY4MMZmAniV9lcKHQtgb8pryMgoEpk0MjJqqE1acxuq8ynvYoy9g68jotcBlGoqAqgSQW9tuzm1bbh56o0AnhZ2XVchte24IpkG4DUi2lfIBWVkFIBMGhkZABhj4wDcCeCfiGgpgC8A+O9St1+hOk8DUFUaXYyxp1Gdz34p9k9+80ZU52/gpqq+GukAwwnkXAC/THwpGRmFIs+nkZHhCMbYPAD/SESXMMZeBDAvxhfBGLsTwFVEtCzZSWZkFIysNDIyHEFEfwbwCGNsMoBKJGGMAXB3JoyMdkNWGhkZGRkZzshKIyMjIyPDGZk0MjIyMjKckUkjIyMjI8MZmTQyMjIyMpyRSSMjIyMjwxmZNDIyMjIynJFJIyMjIyPDGf8/g20kz6oq58QAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u1.draw('phase')\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Scalar_field_X\n", " - x: (4096,), u: (4096,)\n", " - xmin: -3000.00 um, xmax: 3000.00 um, Dx: 1.47 um\n", " - Imin: 0.00, Imax: 1.00\n", " - phase_min: -179.70 deg, phase_max: 180.00 deg\n", " - wavelength: 0.55 um\n", " - date: 2023-10-07_17_30_14\n", "\n" ] } ], "source": [ "print(u1)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## CZT algorithm" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "CPU times: user 3.3 s, sys: 219 µs, total: 3.3 s\n", "Wall time: 3.31 s\n" ] } ], "source": [ "%%time\n", "zs = np.linspace(focal - 50*mm, focal + 50*mm, 1024)\n", "u_z = u1.CZT(zs,0,0)\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Scalar_field_Z\n", " - z: (1024,), u: (1024, 1)\n", " - zmin: 200000.00 um, zmax: 300000.00 um, Dz: 97.75 um\n", " - Imin: 4.36, Imax: 261.70\n", " - phase_min: -179.73 deg, phase_max: 179.85 deg\n", " - wavelength: 0.55 um\n", " - date: 2023-10-07_17_30_18\n", "\n" ] } ], "source": [ "print(u_z)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Scalar z" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "DOF = 6.64 mm\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAD4CAYAAADy46FuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAqzUlEQVR4nO3deXxV9Z3/8dcngZAExCSsYd8SxqWFQkpBOmjRh1VstbVOpY+f1rZWLDp1qMu0dDpDS7W2Lp12xnHBtrZuVQQdW5e6IKhFKRMQUPaAELayhS1kIcn9/P6494aEhJDk3puEe97PxyOPHL5nud97uLnv+/me5Zq7IyIicjIp7d0BERHp2BQUIiLSJAWFiIg0SUEhIiJNUlCIiEiTOrV3BwB69uzpQ4YMae9uiIicVpYtW7bP3Xsl+nE6RFAMGTKEwsLC9u6GiMhpxcy2tsXjaOhJRESapKAQEZEmKShERKRJCgoREWmSgkKkMaWlMGsW9OoFKSnh37NmhdtFAqZDnPUk0qGUlsL48bBpE1RUhNv27YN77oH582HJEujWrX37KNKGVFGInOjee+uHRFRFRbj93nvbp18i7URBIXKiBx9sGBJRFRXw0ENt2x+RdqagEDnR/v2xzRdJMgoKkRP16BHbfJEko6AQOdFNN1Fh1vi89HSYPr1t+yPSzjrGWU/r18MFF9Rv++pX4aaboKwMpkxpuM43vhH+2bcPrrqq4fzp0+Hqq2HbNrj22obzb7sNvvjF8GPfeGPD+T/6EVx0EaxYATNmNJz/s5/BeefBe+/BD3/YcP6vfgWjR8Obb8Kddzac/8gjMHIk/PnPcP/9Dec/8QQMHAjPPtv4mPi8edCzJ/z+9+GfE73yCmRmhsfb585tOH/RovDv++6Dl16qPy8jA159NTz905/CggX15/foET77B2DmTHj//frzBwyAJ58MT8+YEd6HdeXnw5w54elp02DDhvrzR48O7z+Aa66B7dvrz58wAe6+Ozz9la80HAq68EL4938PT196KZSX15//hS/A7beHp0983QEHJk1ihzvDgYw67W6GpaRAbm64Qa89vfbi/NqL+X0vQVRRiJxg3ZYtjAf+PGIEdO6MA8eAA2ecAWPGhKsKkQAxd2/vPlBQUOC6e6x0FDfddBMPPfQQ99xzD3fccQevvfYal1xyCZ/5zGdYsmRJe3dPpJaZLXP3gkQ/jioKkRMsW7YMgLFjxwLwmc98BoAVK1ZQXV3dbv0SaS8KCpE6qqqqWLlyJQBjxowBICsri8GDB1NZWcnGjRvbs3si7UJBIVLH5s2bqaysZPDgwWRlZdW2jxo1CqA2RESCREEhUsf69esBGDlyZL32aFB8+OGHbd4nkfZ2yqAws4FmttDM1prZajP7l0j7j81sh5mtiPxMqbPOTDMrMrP1Zvb5RD4BkXjaEDldMj8/v1579N8aepIgas51FNXAbe6+3MzOAJaZ2RuRef/p7vfVXdjMzgamAucA/YA3zSzf3Wvi2XGRRDhZRZGXlwcoKCSYTllRuPsud18emT4CrAX6N7HKFcAz7l7p7h8DRcC4eHRWJNFOFRRFRUV0hFPKRdpSi45RmNkQ4FPA3yJN/2xmq8zsd2aWHWnrD2yrs9p2GgkWM5tmZoVmVrh3796W91wkAU429JSTk0N2djalpaXs3r27Pbom0m6aHRRm1g2YD8xw98PAQ8BwYDSwC4jeC6Cxm+Q0+Ajm7nPcvcDdC3r16tXSfovE3ZEjR9i9ezddunRh4MCBDeZr+EmCqllBYWadCYfEU+7+PIC773b3GncPAY9yfHhpO1D3r2wAsDN+XRZJjK1btwIwePBgUlIa/mkoKCSomnPWkwG/Bda6+y/rtOfWWezLwEeR6T8BU82si5kNBfKApfHrskhiFBcXA+GgaMzw4cMB2LRpU5v1SaQjaM5ZTxOBa4EPzWxFpO2HwNfMbDThYaUtwI0A7r7azOYCawifMXWzzniS00E0KAYNGtTo/GiAbNu2rdH5IsnqlEHh7n+l8eMOrzSxzl3AXTH0S6TNRYeeThYU0eMWCgoJGl2ZLRJxqqGnaIBElxMJCgWFSMSphp7qVhShUKjN+iXS3hQUIhGnCorMzEx69OhBVVUVe/bsacuuibQrBYUIUF1dzY4dOwAYMGDASZfT8JMEkYJCBNizZw81NTX07t2bLl26nHS5aFDogLYEiYJCBPj73/8OQN++fZtcLnqcQhWFBImCQgTYtWsXALm5uU0u179/+LZlO3fqZgMSHAoKEZofFNH50eVFgkBBIYKCQqQpCgoRFBQiTVFQiHD8jf9UB7MVFBJECgoRjp/1dKqKIicnh86dO3Pw4EHKy8vbomsi7U5BIULzh55SUlLo06cPcDxcRJKdgkICz92bHRR1l9HwkwSFgkIC78CBAxw7dozu3buTmZl5yuUVFBI0CgoJvJZUE3WXU1BIUCgoJPCae8ZTlIJCgkZBIYEXvWV4S4NCB7MlKBQUEnj79u0DoGfPns1aXhWFBI2CQgJv//79APTo0aNZyysoJGgUFBJ4La0oevXqBcDevXsT1ieRjkRBIYHX0oqiblC4e8L6JdJRKCgk8FpaUWRmZtK1a1eOHTvG4cOHE9k1kQ5BQSGB19KKAqB3796Ahp8kGBQUEngtrSjg+PBT9NRakWSmoJDAU0Uh0jQFhQRaWVkZ5eXlpKWl0bVr12avp4pCgkRBIYEWrSZ69uyJmTV7vWhFoaCQIFBQSKC1ZtgJdC2FBIuCQgKtNQeyQRWFBMspg8LMBprZQjNba2arzexfIu05ZvaGmW2M/M6us85MMysys/Vm9vlEPgGRWLS2otDBbAmS5lQU1cBt7n4WMB642czOBn4ALHD3PGBB5N9E5k0FzgEuAR40s9REdF4kVq2tKHQwW4LklEHh7rvcfXlk+giwFugPXAH8IbLYH4AvRaavAJ5x90p3/xgoAsbFud8icaGKQuTUWnSMwsyGAJ8C/gb0cfddEA4ToHdksf7AtjqrbY+0nbitaWZWaGaF+mOT9hJrRbF3715CoVDc+yXSkTQ7KMysGzAfmOHuTd3gprFzDBvcOc3d57h7gbsXRP/oRNpaayuKLl260L17d6qrqzl48GACeibScTQrKMysM+GQeMrdn4807zaz3Mj8XCA6WLsdGFhn9QHAzvh0VyS+ohVFS4MCdIqsBEdzznoy4LfAWnf/ZZ1ZfwKui0xfB7xYp32qmXUxs6FAHrA0fl0WiZ+6F9y1lE6RlaDo1IxlJgLXAh+a2YpI2w+BnwNzzex6oBj4JwB3X21mc4E1hM+Yutnda+LdcZF4aO3QE6iikOA4ZVC4+19p/LgDwIUnWecu4K4Y+iXSJlp7MBtUUUhw6MpsCazKykpKS0vp1KkT3bt3b/H60YoiGjYiyUpBIYFVd9ipJTcEjIpWIRp6kmSnoJDAiuX4BOgYhQSHgkICK5bjE3XX09CTJDsFhQSWKgqR5lFQSGDFcrEdqKKQ4FBQSGDFcrEd1K8o3BvcpUYkaSgoJLBirSgyMzNJT0+nsrKSo0ePxrNrIh2KgkICK9aKwsx0nEICQUEhgRXrwWzQcQoJBgWFBFasp8eCznySYFBQSGCpohBpHgWFBJYqCpHmUVBIIFVVVXHo0CFSUlLIyspq9XZUUUgQKCgkkEpKSgDIzs4mJaX1fwaqKCQIFBQSSLGeGhulikKCQEEhgRTrxXZRqigkCBQUEkiqKESaT0EhgRSPU2NBFYUEg4JCAikep8YC5OTkYGYcOHCA6urqeHRNpMNRUEggxauiSE1NJScnp942RZKNgkICKV4VRd1t6DiFJCsFhQRSvCoK0HEKSX4KCgmkeJ0eC6ooJPkpKCSQ4nV6LKiikOSnoJBAimdFEQ0KVRSSrBQUEjg1NTUcOHAAoPaMpVhEqxJVFJKsFBQSOAcOHMDdycrKolOnTjFvTxWFJDsFhQROPI9P1N2OKgpJVgoKCZx4nhoLqigk+Z0yKMzsd2a2x8w+qtP2YzPbYWYrIj9T6sybaWZFZrbezD6fqI6LtFY8L7arux1VFJKsmlNR/B64pJH2/3T30ZGfVwDM7GxgKnBOZJ0HzSw1Xp0ViYdEVhTuHpdtinQkpwwKd38HKGnm9q4AnnH3Snf/GCgCxsXQP5G4i+epsQCZmZlkZGRQWVlJaWlpXLYp0pHEcozin81sVWRoKjvS1h/YVmeZ7ZG2BsxsmpkVmlmhSnZpS/E+mA06TiHJrbVB8RAwHBgN7ALuj7RbI8s2Wou7+xx3L3D3gugfmUhbiHdFATpOIcmtVUHh7rvdvcbdQ8CjHB9e2g4MrLPoAGBnbF0UiS9VFCIt06qgMLPcOv/8MhA9I+pPwFQz62JmQ4E8YGlsXRSJL1UUIi1zystSzeyPwAVATzPbDswCLjCz0YSHlbYANwK4+2ozmwusAaqBm929JiE9F2klVRQiLXPKoHD3rzXS/Nsmlr8LuCuWTokkkioKkZbRldkSKO5OSUn4bO+4BUVpKZ9//332AHf/4hfQqxfMmgU6VVaShIJCAuXQoUPU1NRwxhlnkJaWFvsGS0th/HhGv/46vYic9rdvH9xzD4wfr7CQpKCgkECJ+7DTvffCpk10qqqq315RAZs2heeLnOYUFBIcpaWk//zn7AE2bdkSnyGiBx8Mh0JjKirgoYdav22RDkJBIcEQGSLKfeIJehF54cdjiChyBlWr54ucBhQUEgyRIaLUY8fqt8c6RHSqIaw4nlkl0l4UFBIMiRoiuukmSE9vfF56Okyf3rrtinQgCgoJhkQNEd1xBwwf3jAs0tPD7Xfc0brtinQgsX9hcDysXw8XXFC/7atfDX9aKyuDKVMarvONb4R/9u2Dq65qOH/6dLj6ati2Da69tuH8226DL34x/Ng33thw/o9+BBddBCtWwIwZDef/7Gdw3nnw3nvwwx82nP+rX8Ho0fDmm3DnnQ3nP/IIjBwJf/4z3H9/w/lPPAEDB8Kzzzb+aXfePOjZE37/+/DPiV55BTIzw5+k585tOH/RovDv++6Dl16qPy8jA159NTz905/CggX15/foAfPnh6dnzoT3368/f8AAePLJ8PSMGeF9WFd+PsyZE56eNg02bKg/f/To8P4DuOYa2L69/vwJE+Duu8PTX/lKwzf5Cy+Ef//38PSll0J5OXTqBCeemVRXaurx12BLX3vZ2dCnD1XFxXRypyYzk0633w5f+xp84QsN19drL/w7KK+9ur7wBbj99vD0ie95EPv7XoKoopBg6NcPUk7yck9JCc9vrdRUGDKE9T168Dbwf9OmwU9+Al27tn6bIh2IdYRv5CooKPDCwsL27oYks8hZT5Vr19IlFDreHh0iWrIEunWL6SFuvPFG5syZw4MPPsh0HZuQNmBmy9y9INGPo4pCgqFbN1iyhN/26MEewM3C11H867/GJSRA93uS5NUxjlGItIVu3bizUyduBrYVFzNgwIC4bl53kJVkpYpCAsPdE3Ln2ChVFJKsFBQSGKWlpVRVVZGZmUlGRkbct6+KQpKVgkICI5HVBKiikOSloJDAiH6zXaKCQhWFJCsFhQRG9A08+oYeb3Urio5w2rlIvCgoJDCiQRHP78quKzMzk8zMTI4dO0apvrBIkoiCQgIjeuwgUUEBx6sVHaeQZKKgkMBI9NATHA8hHaeQZKKgkMBI9NATqKKQ5KSgkMBoi6EnVRSSjBQUEhhtMfSkikKSkYJCAqMthp5UUUgyUlBIYOgYhUjrKCgkEEKhUMKvzAZVFJKcFBQSCAcOHCAUCpGVlUXnzp0T9jiqKCQZKSgkENpi2Knu9lVRSDJRUEggRD/hJ/KMp7rbV0UhyeSUQWFmvzOzPWb2UZ22HDN7w8w2Rn5n15k308yKzGy9mX0+UR0XaYm2qiiys7NJSUnh4MGDVFVVJfSxRNpKcyqK3wOXnND2A2CBu+cBCyL/xszOBqYC50TWedDMUuPWW5FWaqugSElJqT1YHj14LnK6O2VQuPs7QMkJzVcAf4hM/wH4Up32Z9y90t0/BoqAcfHpqkjrtdXQE+g4hSSf1h6j6OPuuwAiv3tH2vsD2+ostz3S1oCZTTOzQjMr1HiuJFpbVRSg4xSSfOJ9MNsaaWv0G1zcfY67F7h7QVt8ypNga8ug0FeiSrJpbVDsNrNcgMjvPZH27cDAOssNAHa2vnsi8dEW93mK0leiSrJpbVD8CbguMn0d8GKd9qlm1sXMhgJ5wNLYuigSu7a4c2yUKgpJNp1OtYCZ/RG4AOhpZtuBWcDPgblmdj1QDPwTgLuvNrO5wBqgGrjZ3WsS1HeRZmvLoFBFIcnmlEHh7l87yawLT7L8XcBdsXRKJJ7cnd27dwPQp0+fhD+eKgpJNroyW5JeaWkp5eXlZGRk0K1bt4Q/nioKSTYKCkl6dasJs8ZOzIsvVRSSbBQUkvTactgJVFFI8lFQSNJrr6DYu3cvoVCoTR5TJJEUFJL02joo0tPTyc7OpqqqipKSE+9+I3L6UVBI0mvroADIzc0FYOdOXW8qpz8FhSS99giKfv36AQoKSQ4KCkl6CgqR2CgoJOkpKERio6CQpNeexyh27drVZo8pkigKCkl6qihEYqOgkKRWVlZGaWkpaWlpnHnmmW32uAoKSSYKCklq0Wqib9++bXL7jigFhSQTBYUktegxgr59+7bp40Yf7+9//7uuzpbTnoJCktr27dsB6N+/0a9uT5j09HRycnKorq7WPZ/ktKegkKS2Y8cOoO2DAjT8JMlDQSFJrSMERbSqETldKSgkqbVnUAwaNAiAbdu2tflji8STgkKSWnsGxeDBgwEoLi5u88cWiScFhSS1jlBRbN26tc0fWySeFBSStNy9QwSFKgo53SkoJGmVlJRQWVlJ9+7d6datW5s/voaeJFkoKCRptWc1EX1cM2PHjh1UVVW1Sx9E4kFBIUmrvYMiLS2N3NxcQqGQrqWQ05qCQpJWewcFaPhJkoOCQpJWe92+oy6d+STJQEEhSevjjz8GYOjQoe3WB1UUkgwUFJK0tmzZAsCQIUParQ/RoIiGlsjpSEEhSasjVBQjRowAoKioqN36IBIrBYUkpWPHjrFjxw7MjIEDB7ZbP/Lz8wHYsGFDu/VBJFYxBYWZbTGzD81shZkVRtpyzOwNM9sY+Z0dn66KNN+2bdsIhUIMGDCAtLS0duvHwIEDSUtLY+fOnZSWlrZbP0RiEY+K4nPuPtrdCyL//gGwwN3zgAWRf4u0qY5wfAIgNTVVw09y2kvE0NMVwB8i038AvpSAxxBpUkc4PhGVl5cHaPhJTl+xBoUDr5vZMjObFmnr4+67ACK/eze2oplNM7NCMyvcu3dvjN0QqS/66X3YsGHt3JPjxyk2btzYzj0RaZ1OMa4/0d13mllv4A0zW9fcFd19DjAHoKCgwGPsh0g969aFX4pnnXVWO/dEFYWc/mKqKNx9Z+T3HuAFYByw28xyASK/98TaSZGWigbFP/zDP7RzT1RRyOmv1UFhZl3N7IzoNHAx8BHwJ+C6yGLXAS/G2kmRlqiqqmLTpk2YWe2n+fY0cuRIANasWYO7imc5/cQy9NQHeMHMott52t3/Ymb/B8w1s+uBYuCfYu+mSPNt2rSJ6upqhg4dSkZGRnt3hz59+tCrVy/27t1LcXFx7dXaIqeLVgeFu28GRjXSvh+4MJZOicSiIw07AZgZo0aN4s0332TlypUKCjnt6Mps6ZDKysr47//+b775zW9y9913s3///mav+9FHHwEdJygARo8eDcCKFSvatR8irRHrWU8icbdjxw4uvvhi1qxZU9v2wAMP8Je//IVPfOITp1x/2bJlAIwZMyZhfWypUaPCxffKlSvbuSciLaeKQjqUsrIyLr/8ctasWcPIkSP59a9/zYQJE9i5cyeXXnpps74pLhoUY8eOTXR3my1aUUT71hJ79uzhoYceYvbs2fzv//4vx44di3PvRE7B3dv9Z+zYsS7i7v69733PAR82bJjv27fP3d3Ly8v9s5/9rAN+2WWXeSgUOun6u3fvdsC7devm1dXVbdXtU6qurvbu3bs74Nu2bWvWOqFQyH/5y196ly5dnPDFrQ54Xl6eL1myJME9ltMBUOht8B6tikI6jFWrVvFf//VfpKSkMHfuXHr06AFAeno6zz77LGeeeSYvv/wyTz/99Em3UXfYKTU1tU363Rypqamcd955ALz77runXN7dueOOO7j11luprKzksssu4/bbbyc/P5+NGzdywQUX8Prrrye62yKAhp6kA5kxYwY1NTXcfPPNDYaN+vXrx/333w/ArbfeyqFDhxrdxuLFiwH49Kc/ndjOtsKkSZOA5gXFfffdx/3330/nzp157rnneOmll7j33nv58MMP+fa3v01FRQVXXXUVH374YaK7LaKhJ+kYFi5c6IBnZWX5gQMHGl0mFArVDkHNmDGj0WXGjx/vgL/00ksJ7G3rvPvuuw74WWed1eRyixcv9tTUVAf82WefbTC/pqbGr7766tohuiNHjiSqy9LB0UZDT+0eEq6gEHe/4IILHPDZs2c3udwHH3zgKSkpnpqa6h9++GG9eSUlJZ6amuqdOnXyw4cPJ7K7rVJZWelZWVkO+Lp16xpdpqSkxAcNGuSA33HHHSfdVnl5uY8aNcoBv+GGGxLVZeng2iooNPQk7W7RokUsWrSIrKwsbrnlliaXHT16NNOnT6empobvfve74U87Ec8//zw1NTWcf/75nHHGGYnudoulpaVx+eWXAzB//vwG892dG264geLiYj796U9z5513nnRb6enpPPnkk6SlpfHoo4/y0ksvJazfIu1eTbgqisBrbjURtX//fu/Ro4cD/swzz9S2T5w40QH/zW9+k6iuxuzFF190wM8999wGZ2/9z//8jwN+xhln+KZNm5q1vfvuu88B79ev30mH7CR5oaEnCYK333679tjEwYMHm73enDlzHPC+ffv6li1b/LXXXnPAs7Oz/dChQwnscWwqKiq8b9++Dvirr75a27506VJPS0tzwP/4xz82e3vV1dU+YcIEB/xb3/pWIrosHZiCQgLhwgsvdMB//OMft2i96upqP//88x3wnJwc79q1qwP+s5/9LEE9jZ9f/OIXDnh+fr6XlJT48uXLayukm266qcXbW7t2be21Fn/5y18S0GPpqBQUkvT++te/OuDdu3f3kpKSFq9/8OBBP++882ovRLviiiu8qqoqAT2Nr4qKCj/77LNrh5miZzhNmTLFKysrW7XNu+++2wEfOHBgh66oJL7aKih0MFtaZcOGDcyYMYNx48YxduxYbrjhBpYuXdrs9d2dH/3oRwDccsstZGdnt7gPZ555Ju+++y7vvPMO7733Hi+88AKdOnX825d16dKFV199lTFjxnDkyBEApk+fzvz580lLS2vVNm+//XYKCgrYtm0b3//+9+PZXRFVFNIy1dXV/pOf/KT2U/CJP1//+tebdaxh/vz5tcNG+/fvb4OedzyhUMiLiopqb1USq1WrVnnnzp0d8Lfeeisu25SODVUU0tFUVFRw5ZVXMmvWLEKhENdffz0LFy5k8eLFfP/73ycjI4PHH3+cT33qU3zwwQcn3U5paSm33347ALNnzyYnJ6etnkKHYmYMHz689lYlsfrEJz5RW6Vdf/31lJaWtngboVCIRYsWcdtttzF58mTGjBnDxIkT+da3vsUf//jHVm1TkkBbpNGpflRRdHxHjx71iy++uLYKeOONNxoss27dOh8zZowD3qVLF//tb3/bYJlQKOTf/OY3HfBPfvKTp8UxhdPJsWPHai/Eu/LKK72mpqbZ6z366KM+ZMiQRivF6E9WVpbPnDmzVceUJP7QwWyJl1Ao5Hv27PGtW7e26mDpkSNHaq916NWrl69cufKky5aXl/u0adNq31iuv/56Ly8vd/fwrSdmzpxZGyQnXlkt8bFmzRo/88wzHfDvfOc7Td5Ft7q62p966ikfMWJE7f/Z4MGDfebMmf7yyy97YWGhL1y40O+///7a26MA3rNnT3/kkUeaHUSSGAoKidnGjRv9pptu8t69e9f+gWdkZPjFF1/szz//fLNuw33o0KHaC9lyc3N9zZo1zXrsxx57zNPT0x3w/v37+9e//vXaaiMlJcXnzZsX69OTJixYsKD2uozLL7/cd+/eXW9+VVWVz507t/bsKyK3L3/66aebfPN///33fdKkSbXrjBs3zpcvX57opyMnoaCQVisrK/NZs2bV+x6D7t27e79+/Rp8r8Hjjz9+0uGfbdu2eUFBgQM+YMAA37hxY4v68cEHH/i5555b7zH79OnjL7/8cjyeppzCwoULayuLzMxMnzp1qs+cOdOvueaa2ov+AB80aJD/5je/afYwYCgU8meffbb29ZSSkuIzZsyI2/21SktLff369f7RRx/5mjVrWnQhZtAoKALo4MGD/sADD/i1117rU6ZM8alTp/pdd93lS5cubfLLeup65ZVXfNiwYfXOQlqxYkXt+nv27PFf//rXPnTo0Npl8vPz/fHHH/eysjJ3D5/n/9hjj3mvXr0c8KFDhzb7lhInCoVCvnjxYn/00Uf9hRde8KNHj7ZqO9I6mzdv9ssuu6zR4w35+fn+8MMPt/rajUOHDvmMGTM8JSWltnJ87rnnmv1adQ+/PjZu3OgPP/ywX3XVVd6nT5+THhu5+OKLffbs2S36e0h2CorTRCgU8uLiYl+8eLG/9957vnnz5ha/iJcvX+433HCDZ2ZmnvQgYn5+vt95552+ZcuWRrexbt06/+IXv1i7/LnnnuvvvPPOSR+zqqrKH3vssXqh0qVLFx8+fHi9flx00UW+d+/eFj0f6XjWr1/vDz/8sM+ePdsfeeQRX758edzebJcvX+7jxo2rfc2MGTPGH3vssZOe9lxcXOxPPfWUX3/99T548OAGr/W0tDQfNmyYn3XWWZ6Xl9fo38WgQYP8e9/7ni9evDjQx0naKigs/Fjtq6CgwAsLC9v0MQ8ePMjmzZs5fPgwGRkZDB06lN69ezdrXXdn5cqVPPfcc8ydO5eioqJ687OysjjvvPOYPHkykydPZtSoUaSk1D8T+ejRo8ybN4+HHnqIv/3tb7XtkydP5uqrr6Zv376UlJSwdOlSnn/+eXbv3l27zKRJk5g8eTL9+/dn//79LFy4kNdffx13p1u3bvz4xz/mlltuoXPnzqd8LlVVVTz55JM88MADfPDBB0RfD6NGjeLWW2/lmmuuadB3kRPV1NQwZ84cZs2axd69e4Hwt/qNHDmSIUOG0LlzZw4cOEBRUVGD7z3Pycnhwgsv5KKLLmLy5MkMGzas3mvO3dmxYwfvvfceixYt4sUXX6y3jX79+vGVr3yFq666iokTJ7b4mw0PHDjAhg0ban+KioooKSnh0KFDVFRUkJ6eTkZGBr1792bw4MEMGTKET37yk4waNard71JsZsvcvSDhj3M6BcWRI0d47733WL16Nbt27SIlJYXs7GzOPvtsRo8ezcCBAzGzRtetrq5myZIlvPLKK7zyyiusXLmywTK5ubmMGzeOCRMmMH78eAoKCujatWvt+kuXLuWll17iueeeqxcO2dnZ5OXlYWZs2bKl3ps6hP8Qzj//fAoKCjh27BirVq3itddeo6ysDAgHyze+8Q2+853vMHLkyEb7/sYbb/DEE0/wwgsvUFFR0WCZtLQ0rrvuOmbPnk3fvn1PuS8bc/DgQfbs2UPPnj0De22DxKa8vJynn36aZ555hkWLFlFdXd1gmaysLCZOnMikSZO46KKLGD16dIs+jIRCIZYsWcK8efOYP38+xcXFtfN69OjBxIkTGT9+PCNGjGDQoEF069YNM+Po0aPs2rWLnTt3sm7dOlavXl37XtIaZkZeXh5jxoxh7NixjBkzhjFjxpCVldWq7dVVWlrKzp072blzJ4cOHaKyspKKigpCoRCZmZlkZmbSrVs3Pve5zwUnKLKysvy73/0uo0aNYtSoUfTv35/U1FR27drFqlWrePvtt3nnnXdYvnw5oVDopNvp2bMnY8eOZezYsQwdOpTU1FS2b9/OsmXLWLhwIYcPH65dNj09nREjRpCTk8PRo0fZsGFD7e0UolJTUxk+fDgAxcXF9d6ge/XqxZVXXslXv/pVJk2aVO/WEdu3b+ftt9/mrbfeYsGCBWzdurXR/k6YMIFvf/vbTJ06lczMzGbtq8OHD/Pqq69SWFjIgQMH6Nq1K6NHj+byyy+P24VbIvEQ/bsqLi6urXaHDh3K0KFD41alujuFhYXMmzeP5557jo8//rjF28jIyGDkyJHk5+eTn59PXl4evXr1Iisri/T0dCoqKigrK2PXrl1s3bqVoqIiVqxYwerVq6mqqmqwveHDhzNmzJjaoBowYADdu3cnIyOD1NRUysrKKC8vp6SkpDa4du7cya5du9ixYwc7d+5s8F7UhOAEhZk1qxOdOnWioKCAgoIC+vfvD8CePXtYtWoVH3zwASUlJU2un5eXx5QpU7jsssuYNGkSXbp0qZ0XCoXYtGkTS5Ys4f3332fJkiWsWrWKmpqa2mVGjBjBJZdcwpe//OUG4dCUjz/+mLfeeov169eTlpZGfn4+kyZNYsiQIc1aX0ROzd3ZvHkz77//PoWFhWzdupXi4mLKy8uB8D22cnNzyc3NJS8vj3POOYdzzjmHIUOGtCq4KisrWbNmDcuXL2fZsmUsW7aMlStXUllZGfNzSU9Pp1+/fvTr14/s7GzS09NJT0/HzCgrK6OsrIwjR47w7rvvBicohg0b5lOnTmXlypWsWrWKffv2UVVVRd++fRkxYgT/+I//yPnnn8+ECRNqh4JO5O5s3bqVwsJCli1bxt69e6mqqiI3N5eRI0cyefJkBg8e3KJ+HT16lC1btmBmtZ8KREROpqqqirVr17JixQq2bt3K1q1b2b59O0ePHqW8vJzq6mq6du1KZmYm3bt3rw2DaID179+f/v37k5WVddJh9Lp0jEJERJrUVkGh01lERKRJCgoREWlSwoLCzC4xs/VmVmRmP0jU44iISGIlJCjMLBX4H+BS4Gzga2Z2diIeS0REEitRFcU4oMjdN7v7MeAZ4IoEPZaIiCRQor5guD+wrc6/twOfqbuAmU0DpkX+WWlmHyWoL6ebnsC+9u5EB6F9cZz2xXHaF8c1vJVDAiQqKBo7AbjeebjuPgeYA2BmhW1xitfpQPviOO2L47QvjtO+OM7M2uS6gkQNPW0HBtb59wBg50mWFRGRDixRQfF/QJ6ZDTWzNGAq8KcEPZaIiCRQQoae3L3azP4ZeA1IBX7n7qubWGVOIvpxmtK+OE774jjti+O0L45rk33RIW7hISIiHZeuzBYRkSYpKEREpEkxBYWZDTSzhWa21sxWm9m/RNpzzOwNM9sY+Z1dZ52Zkdt6rDezz9dpH2tmH0bm/ZdF7rFrZl3M7NlI+9/MbEidda6LPMZGM7sulucSqyb2xb1mts7MVpnZC2aWVWedQO2LOvNvNzM3s5512gK3L8zsu5Hnu9rM7qnTHqh9YWajzWyJma0ws0IzG1dnnWTdF+lmttTMVkb2xU8i7R3zvTOWL9wGcoExkekzgA2Eb9lxD/CDSPsPgF9Eps8GVgJdgKHAJiA1Mm8pMIHwNRivApdG2m8CHo5MTwWejUznAJsjv7Mj09lt8UXjLdwXFwOdIu2/CPK+iPx7IOGTHLYCPYO6L4DPAW8CXSLzegd4X7xe57lMARYFYF8Y0C0y3Rn4GzCeDvreGVNF4e673H15ZPoIsJbwVdlXAH+ILPYH4EuR6SuAZ9y90t0/BoqAcWaWC3R39/c9/EweP2Gd6LbmARdGEvPzwBvuXuLuB4A3gEtieT6xONm+cPfX3T36xcFLCF9TAgHcF5HZ/wn8K/UvwAzivpgO/NzdKyPz9kRWCeK+cCD6rWBncvyaq2TeF+7upZF/do78OB30vTNuxygiZc2nCCdjH3ffBeEXB9A7slhjt/boH/nZ3kh7vXUib7iHgB5NbKvdnbAv6voW4cSHAO4LM7sc2OHuK09YLHD7AsgH/jEyJPC2mX06slgQ98UM4F4z2wbcB8yMLJbU+8LMUs1sBbCH8Bt3h33vjEtQmFk3YD4ww90PN7VoI23eRHtr12k3J9sXZvZvQDXwVLSpkdWTdl8Qfu7/BvxHY4s20pa0+yLyuuhEuOwfD9wBzI182gvivpgOfM/dBwLfA34bXbSR1ZNmX7h7jbuPJjzKMM7Mzm1i8XbdFzEHhZl1Jvyf/pS7Px9p3h0piYj8jpbVJ7u1x3aOD8nUba+3jpl1IlyaljSxrXZzkn1B5GDRF4D/FykPIXj7YjjhsdWVZraFcB+Xm1lfgrcvINzP5yNDEEuBEOGb3QVxX1wHRKefI3z3aUjyfRHl7geBRYSHfzrme2ccDsg8DvzqhPZ7qX9A5p7I9DnUPyCzmeMHZP6P8Ker6AGZKZH2m6l/QGauHz8g8zHhT2XZkemcWJ5PgvbFJcAaoNcJ7YHbFycss4XjB7MDty+A7wCzI9P5hIcCLKD7Yi1wQWT6QmBZAF4XvYCsyHQG8C7hD5Md8r0z1if7WcIlyypgReRnCuFxsAXAxsjvnDrr/BvhI/briRydj7QXAB9F5j3A8avG0wl/yigifHR/WJ11vhVpLwK+2V7/6afYF0WE3wSibQ8HdV+csMwWIkERxH0BpAFPRp7bcmBygPfFZ4FlhN8I/waMDcC++CTwQWRffAT8R6S9Q7536hYeIiLSJF2ZLSIiTVJQiIhIkxQUIiLSJAWFiIg0SUEhIiJNUlCIiEiTFBQiItKk/w9MlJt0VSSNXAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "dof = u_z.DOF(percentage=0.5, remove_background=True, has_draw=True)\n", "print(\"DOF = {:2.2f} mm\".format(dof / mm))\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEJCAYAAACDscAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAArMklEQVR4nO3de3iU9Z338fc3R0BCOIUkEBQ5KKciVDwhVWy1ou3W2uqqZa3Xo091rbaeul2ta9vH1W6lrbb2oXVxdeui1FotFYuilgd1YVUIRwWVQxGMhHAIEI4JJN/njzkQYJLMJHNnZpLP67rmyuSe+77nm9+VzCe/3+8+mLsjIiISr6xUFyAiIplFwSEiIglRcIiISEIUHCIikhAFh4iIJCQn1QW0Rd++fX3QoEGpLkNEJKMsWbJku7sXtXb7jA6OQYMGUV5enuoyREQyipltbMv2GqoSEZGEKDhERCQhCg4REUmIgkNERBKi4BARkYQoOEREJCEKDhERSYiCQyTJtm3bxpNPPsm6detSXYpIIDL6BECRdHPw4EEmTpzImjVrKCwsZMWKFZx00kmpLkskqdTjEEmiZ555hjVr1gCwe/dupk6dmuKKRJJPwSGSRE888QQAd911FwDPPfcchw4dSmVJIkmn4BBJkp07d/Luu++Sm5vLj3/8Y0aMGMH27dtZuHBhqksTSSoFh0iSzJ8/n4aGBiZMmED37t258MILAfjv//7vFFcmklwKDpEkWbx4MQCf+9znjvqq4JCORsEhkiQrVqwAYOzYsQBMmDABgPLyctw9VWWJJF1gwWFmA81svpl9YGarzOy28PIfm9mnZrY8/Li00Tb3mNk6M/vIzC4OqjaRICxfvhyA0047DYD+/fvTp08fdu7cSUVFRQorE0muIHsch4G73H0EcDZwi5mNDL/2iLuPDT9eBgi/djUwCpgM/MbMsgOsTyRptm3bRmVlJd27d2fw4MEAmFk0RFauXJnK8kSSKrDgcPdKd18afr4H+AAY0MwmlwHPunutu28A1gFnBlWfSDJFhqnGjBlDVtaRP6sxY8YACg7pWNpljsPMBgHjgHfDi241s5Vm9qSZ9QovGwB80mizCmIEjZndaGblZla+bdu2IMsWidvq1asBGD169FHLhw8fDhA9KVCkIwg8OMysO/ACcLu71wC/BYYAY4FK4BeRVWNsftyMortPd/fx7j6+qKjV91oXSar169cDMHTo0KOWR77XdaukIwk0OMwsl1BoPOPufwJw9yp3r3f3BuBxjgxHVQADG21eBmwOsj6RZIkEx5AhQ45aruCQjijIo6oMeAL4wN0fbrS8tNFqlwPvh5/PBq42s3wzOxkYBiwKqj6RZGoqOMrKysjLy2PLli3s3bs3FaWJJF2QPY5zgWuBzx9z6O1UM3vPzFYCFwB3ALj7KuA5YDUwF7jF3esDrE8kKRoaGtiwYQNA9IiqiOzs7GiYRMJFJNMFdll1d19A7HmLl5vZ5kHgwaBqEgnC5s2bqa2tpaioiIKCguNeHzp0KB988AHr1q2LHp4rksl05rhIGzU1TBWheQ7paBQcIm3U1DBVRORGTps2bWq3mkSCpOAQaaPI5UQGDhwY8/XI8k8++STm6yKZRsEh0kaffvopAAMGxL4wQiQ41OOQjkLBIdJGLQXHiSeeCKjHIR2HgkOkjSJDVU0FR1FREXl5eVRXV7N///72LE0kEAoOkTaK9DjKyspivp6VlRV9Tb0O6QgUHCJtUFdXx9atW8nKyqK4uLjJ9TRBLh2JgkOkDSorKwEoKSkhJ6fp82kVHNKRKDhE2qClifEIBYd0JAoOkTZoaWI8orQ0dG3PLVu2BF6TSNAUHCJt0NLEeISCQzoSBYdIG8Q7VFVSUgIcmRMRyWQKDpE2iDc41OOQjkTBIdIGVVVVwJFgaErjHof7cXdEFskoCg6RNti6dSsA/fr1a3a9E044gYKCAmpra9m9e3d7lCYSGAWHSBvEGxygeQ7pOBQcIq1UX1/P9u3bAejbt2+L60eGsxQckukUHCKttGPHDtydPn36NHvWeESkx6EJcsl0Cg6RVopMjMczTAXqcUjHoeAQaaXI/EZzFzdsTD0O6SgUHCKtlMjEOKjHIR2HgkOklRINDh1VJR2FgkOklRINjsiQVmQ7kUyl4BBppUSDI7Letm3bAqtJpD0oOERaKdHgiJzrsX37dhoaGgKrSyRoCg6RVko0OPLy8igsLKS+vp6dO3cGWZpIoBQcIq2U6OG4AEVFRYCGqySzKThEWinRHkfjdRUckskUHCKtsH//fvbu3Ut+fj4FBQVxbxfpcejIKslkgQWHmQ00s/lm9oGZrTKz28LLe5vZ62a2Nvy1V6Nt7jGzdWb2kZldHFRtIm3VuLdhZnFvp6Eq6QiC7HEcBu5y9xHA2cAtZjYSuBuY5+7DgHnh7wm/djUwCpgM/MbMsgOsT6TVWjNMBQoO6RgCCw53r3T3peHne4APgAHAZcBT4dWeAr4afn4Z8Ky717r7BmAdcGZQ9Ym0RSQ4IkEQLwWHdATtMsdhZoOAccC7QLG7V0IoXIDIv2wDgE8abVYRXnbsvm40s3IzK9cfn6RK5D4ciQaHJselIwg8OMysO/ACcLu71zS3aoxlx92c2d2nu/t4dx+f6B+tSLLs2LEDgD59+iS0nSbHpSMINDjMLJdQaDzj7n8KL64ys9Lw66VA5C+oAhjYaPMyYHOQ9Ym0VluDQz0OyWRBHlVlwBPAB+7+cKOXZgPXhZ9fB7zYaPnVZpZvZicDw4BFQdUn0hYKDunMWr7fZeudC1wLvGdmy8PLfgD8FHjOzG4ANgFXArj7KjN7DlhN6IisW9y9PsD6RFqtrcERuV5VVpZOpZLME1hwuPsCYs9bAHyhiW0eBB4MqiaRZGltcOTn59OjRw9qamrYtWsXvXv3DqI8kUDp3x2RVogER+SKt4nQcJVkOgWHSCtEDsdNtMcBCg7JfAoOkQS5e6uHqkDBIZlPwSGSoH379lFXV0fXrl3p2rVrwtvrXA7JdAoOkQS1pbcBRx9ZJZKJFBwiCWprcDS+haxIJlJwiCRIPQ7p7BQcIglSj0M6OwWHSIIUHNLZKThEEpSs4NDhuJKpFBwiCVKPQzo7BYdIgtpyuRGAHj16kJuby759+zhw4EAySxNpFwoOkQS15XIjAGYWDZ1ICIlkEgWHSILaOlQFGq6SzKbgEElQMoNDE+SSiRQcIglSj0M6OwWHSAIOHTpETU0NWVlZFBYWtno/OntcMpmCQyQB1dXVAPTu3btNt31Vj0MymYJDJAHJGKYCBYdkNgWHSAIUHCIKDpGEJDs4dFSVZCIFh0gC2nrWeIQmxyWTKThEEqChKhEFh0hCkhUcke23b9+Ou7e5LpH2pOAQSUCygqNr166ccMIJHDp0iD179iSjNJF2E1dwmFnvoAsRyQTJCg7QBLlkrnh7HO+a2R/N7FIzs0ArEkljjU8AbCvNc0imijc4TgGmA9cC68zsJ2Z2SnBliaSnZPY4dGSVZKq4gsNDXnf3a4D/DVwHLDKzN83snEArFEkj6nGIQE48K5lZH+AfCPU4qoDvALOBscAfgZMDqk8kbbh7tMeh4JDOLK7gAN4GZgBfdfeKRsvLzeyx5Jclkn72799PXV0dXbt2pWvXrm3enybHJVPFO8fxL+7+r41Dw8yuBHD3h2JtYGZPmtlWM3u/0bIfm9mnZrY8/Li00Wv3mNk6M/vIzC5u5c8jEphkzm+AehySueINjrtjLLunhW1+B0yOsfwRdx8bfrwMYGYjgauBUeFtfmNm2XHWJtIukjlMBZocl8zV7FCVmV0CXAoMMLNHG73UAzjc3Lbu/paZDYqzjsuAZ929FthgZuuAMwkNkYmkhcjEuHoc0tm11OPYDJQDB4EljR6zgdYOJ91qZivDQ1m9wssGAJ80WqcivOw4ZnajmZWbWbnGhqU9JbvHoeCQTNVsj8PdVwArzOwZd2+2hxGn3wL/Cnj46y+A64FYJxXGvICPu08ndE4J48eP10V+pN2oxyES0tJQ1XPu/vfAMjNr/CFthE7vGJPIm7l7VaN9Pw78JfxtBTCw0aplhHo7Imkj2T2OyH6qq6upr68nO1vTepIZWjoc97bw1y8n483MrNTdK8PfXg5EjriaDcw0s4eB/sAwYFEy3lMkWZLd48jJyaF3795UV1dTXV0dnSwXSXctDVVFPuS3AwfcvSF8qZHhwCvNbWtmvwcmAX3NrAL4ETDJzMYSGob6GLgp/D6rzOw5YDWhSfdb3L2+lT+TSCCS3eOA0HBVdXU127dvV3BIxoj3BMC3gM+FJ7PnEZowvwqY0tQG4cuTHOuJZtZ/EHgwznpE2l2yexwQCo41a9ZonkMySrzncZi77we+Bvza3S8HRgZXlkj6CarHAZogl8wSd3CEL2Y4BZgTXhZvb0WkQwiqxwG67IhklniD4zZCZ4rPCs9HDAbmB1eWSPpJ9iVHQGePS2aKq9fg7m8RmueIfP834LtBFSWSbtw92uPo1atXC2vHT0NVkonivaz6KcD3gEGNt3H3zwdTlkh6qampob6+noKCAvLy8pK2XwWHZKJ45yn+CDwG/Aegw2Sl00nmDZwaU3BIJoo3OA67+28DrUQkjQUxvwEKDslM8U6Ov2Rm3zazUjPrHXkEWplIGgmqxxGZHNdRVZJJ4u1xXBf++k+NljkwOLnliKQn9ThEjoj3qCrdU1w6tSBO/gPo0aMHOTk57N27l4MHD9KlS5ek7l8kCHENVZlZNzP7FzObHv5+mJkl5cKHIpkgiJP/AMws2uuIhJNIuot3juM/gTpgQvj7CuCBQCoSSUNB9ThAw1WSeeINjiHuPhU4BODuB4h98yWRDimoHgdoglwyT7zBUWdmXQnflc/MhgC1gVUlkmbU4xA5It6jqn4MzAUGmtkzwLnA/wqqKJF0E2SPQ8EhmSbeo6peM7MlwNmEhqhuc3f9lkunoR6HyBHxHlU1z913uPscd/+Lu283s3lBFyeSLtTjEDmi2R6HmXUBuhG6/WsvjkyI9yB0b3CRDq++vp6dO3cC0LNnz6TvX5PjkmlaGqq6CbidUEgs4Uhw1ADTgitLJH3s3r0bd6dnz57k5CT//mXqcUimafavwN1/BfzKzL7j7r9up5pE0kqQ8xug4JDME+/k+K/NbALH34/jvwKqSyRtBDm/AQoOyTzx3shpBjAEWM6R+3E4oOCQDi+oCxxGRPa7fft23B0znVsr6S3eAdvxwEh39yCLEUlHQV1SPaJbt25069aN/fv3s3fvXgoKCgJ5H5FkiffM8feBkiALEUlXQfc4QEdWSWaJNzj6AqvN7FUzmx15BFmYSLoIuscBmueQzJLIJUdEOqXIh3mQPQ4Fh2SSeI+qejPoQkTSVeTDPDKcFAQFh2SSls4cX+DuE81sD+Er40ZeAtzdewRanUgaiMw7RD7cg6DgkEzS0gmAE8NfdZiHdFrt0ePQ5Lhkkngnx0U6LfU4RI4WWHCY2ZNmttXM3m+0rLeZvW5ma8NfezV67R4zW2dmH5nZxUHVJZIId49+mCs4REKC7HH8Dph8zLK7gXnuPgyYF/4eMxsJXA2MCm/zGzPLDrA2kbjs3r2bw4cPU1BQQH5+fmDvo+CQTBJYcLj7W0D1MYsvA54KP38K+Gqj5c+6e627bwDWAWcGVZtIvNpjfgMUHJJZ2nuOo9jdKwHCX/uFlw8APmm0XkV4mUhKtcf8BmhyXDJLukyOx7qqW8zrYpnZjWZWbmbl+iOToLXH/AYcOSu9urqa+vr6FtYWSa32Do4qMysFCH/dGl5eAQxstF4ZsDnWDtx9uruPd/fxQQ8fiET+OQn6dy0nJ4devXrh7tG7DYqkq/YOjtnAdeHn1wEvNlp+tZnlm9nJwDBgUTvXJnKc9upxNH4PzXNIugvycNzfA28Dp5pZhZndAPwUuMjM1gIXhb/H3VcBzwGrgbnALe6u/rqkXHv1OEDBIZkj+TdQDnP3a5p46QtNrP8g8GBQ9Yi0Rnv2OCLhpOCQdJcuk+MiaSkVPQ4d9CHpTsEh0gzNcYgcT8Eh0gzNcYgcT8Eh0gz1OESOp+AQaUJtbS01NTVkZ2fTs2fPwN9Pk+OSKRQcIk3YsWMHEOoJmMW6uEFyaXJcMoWCQ6QJ7XWdqggNVUmmUHCINKGqqgqAkpKSdnk/9TgkUyg4RJqwZcsWAIqLi9vl/QoLC8nPz2fv3r3s27evXd5TpDUUHCJNaO8eh5lRWloKQGVlZbu8p0hrKDhEmtDePQ4gGhybN8e8OLRIWlBwiDShvXscAP379wfU45D0puAQaUIqexwKDklnCg6RJqSyx6GhKklnCg6RJkSCQz0OkaMpOERiOHz4MNu3bycrK6tdLnAYoeCQTKDgEIlh27ZtuDt9+/YlOzu73d5XQ1WSCRQcIjFEJsbbc34D1OOQzKDgEIkhFfMbAH369CE3N5ddu3Zx4MCBdn1vkXgpOERiSFWPQ2ePSyZQcIjEkKoeB+jscUl/Cg6RGD799FPgyGR1eyorKzuqBpF0o+AQiaGiogI48iHengYOHAjApk2b2v29ReKh4BCJIZXBceKJJwIKDklfCg6RGBQcIk1TcIgc49ChQ2zZsoWsrKx2P6oKFByS/hQcIseorKzE3SkuLiY3N7fd31/BIelOwSFyjFQOUwEUFRWRn59PdXW1biEraUnBIXKMVAdHVlZW9L0/+eSTlNQg0hwFh8gxUh0coOEqSW8KDpFjRE68U3CIxJaTijc1s4+BPUA9cNjdx5tZb+APwCDgY+Dv3X1nKuqTzk09DpHmpbLHcYG7j3X38eHv7wbmufswYF74e5F2F5lXGDBgQMpqOOmkkwDYsGFDymoQaUo6DVVdBjwVfv4U8NXUlSKd2ccffwzAoEGDUlbDsGHDAFi7dm3KahBpSqqCw4HXzGyJmd0YXlbs7pUA4a/9Ym1oZjeaWbmZlW/btq2dypXOYv/+/VRWVpKTk5PSoSoFh6SzVAXHue7+WeAS4BYzOy/eDd19uruPd/fx7XkvaOkcGvc22vOWsccqKSmhe/fuVFdXU11dnbI6RGJJSXC4++bw163ALOBMoMrMSgHCX7emojbp3P72t78BMHjw4JTWYWYMHToUUK9D0k+7B4eZnWBmBZHnwBeB94HZwHXh1a4DXmzv2kQiwTFkyJAUV6LhKklfqTgctxiYZWaR95/p7nPNbDHwnJndAGwCrkxBbdLJpUuPAxQckr7aPTjc/W/AaTGW7wC+0N71iDS2fv16QMEh0px0OhxXJOUiH9IaqhJpmoJDJKy2tpZ169aRlZXFqaeemupyojV8+OGHNDQ0pLgakSMUHCJhH330EfX19QwZMoQuXbqkuhz69u1LaWkpe/fu1RnkklYUHCJhq1evBmDUqFEpruSI004LTQeuWLEixZWIHKHgkA6lpqaGn/70p1x55ZXce++9bNmyJe5tV61aBcDIkSODKi9hCg5JRym5Oq5IEDZu3MhFF1101GTy448/zquvvsq4ceNa3D4dexxjxowBFBySXtTjkA5h3759XHbZZaxdu5YxY8bwxBNP8PnPf55t27bxpS99KXqPjeYsXboUgM985jNBlxu3sWPHAkdqS8SWLVt49NFHuffee3n++eepra1NcnXSabl7xj5OP/10F3F3v/322x3wYcOG+c6dO93dvba21idNmuSAX3rppd7Q0NDk9pWVlQ549+7d/fDhw+1UdcsOHz7shYWFDvimTZvi2qahocEffvhhz8vLc0IXFHXAhwwZ4u+8807AFUsmAMq9DZ+96nFIxnvvvff49a9/TVZWFn/4wx/o2bMnAHl5ecycOZPCwkJefvllnn322Sb38e677wJwxhlnpPTihsfKzs5mwoQJACxcuDCube6++27uvPNO6urquOyyy7jvvvsYPnw469ev54ILLmDevHlBliydgIJDMt6dd95JfX09N99883FzGaWlpfz85z+Prrdnz56Y+3jnnXcAOPvss4MtthUmTpwIwIIFC1pc95e//CVTp04lJyeH5557jj//+c/cf//9LF++nOuvv54DBw5w+eWX8+GHHwZdtnRkbemupPqhoSpZsGCBA96jRw/fsWNHzHXq6+v9rLPOcsC///3vx1xn4sSJDviLL74YZLmt8uabbzrgo0ePbna9xYsXe05OjgP+zDPPHPd6fX29X3nllQ748OHDfe/evUGVLGmONg5VpfzDvy0PBYd88YtfdMDvu+++ZtdbtGiRA56bm+tr1qw56rXq6mrPzs72nJwc37VrV5DltsrBgwe9oKDAAV+3bl3Mdfbs2ePDhg1zwL/73e82ua89e/b4qFGjWlxPOra2BoeGqiRjvfPOO7z22mt0796d22+/vdl1zzjjDK6//noOHTrEHXfccdRrr7zyCvX19UycOJHCwsIAK26d/Px8vvzlLwMwa9asmOvcdtttrF27ls985jM89NBDTe6re/fuzJgxg+zsbB599FHeeuutQGqWjk3BIRnrgQceAOA73/kOvXv3bnH9n/zkJ/To0YM5c+bw8ssvR5f/53/+JwBf//rXgyk0CSK1/eEPfzjutWeffZYnn3ySLl268Pvf/77Fy6WMGzeOH/zgBwDReQ+RhLSlu5Lqh4aqOq8lS5Y44N26dfOtW7fGvd0vfvGL6KGpO3fu9HfeeccB79Kli1dXVwdYcdvs27fPe/fu7YAvWLAguvyjjz7y7t27O+C//e1v495fbW1tdMjqn//5n4MoWdIYmuOQzuhrX/uaA37XXXcltF1tba2PGTPGAT/11FO9pKTEAb/77rsDqjR57r33Xgf8jDPO8IMHD/qmTZv85JNPdsCvuuqqZs9TieXtt992M/Ps7GxfunRpQFVLOlJwSKfz/vvvO+D5+fm+efPmhLffsGGDDx06NHpi3Nlnn50RRxjV1NR4WVlZ9ETHXr16RYOkpqamVfu87bbbHPDPfvazfujQoSRXLOmqrcGhOQ5pV3v27OGxxx7jm9/8JlOmTOHhhx+mqqoqoX3cf//9AHzrW9+itLQ04RoGDRrEsmXLmDFjBjNnzmT+/PmccMIJCe+nvRUUFPDSSy9RWlrK2rVr2blzJxdffDFz586loKCgVft84IEHOOmkk1i6dCmPPPJIkiuWDqstqZPqh3ocmeWFF17wPn36HHUZDMC7du3qDz30kNfX17e4j4ULF0a3+eSTT9qh6vSzb98+f/XVV33p0qUJD0/F8sorr0TnedauXZuECiXd0cYeh4X2kZnGjx/v5eXlqS5DWuDu3HfffTz44IMATJgwgWuvvZbc3FxefPFFXnrpJQAuvPBCZs6cSVFRUcz91NfXc84557B48WLuu+++aM9D2u7aa6/l6aef5oILLuD1119v1WVXqqqqeO2111i5ciXbtm3jhBNOYOTIkXzhC19g+PDhAVQtrWVmS9x9fKt30JbUSfVDPY7019DQ4HfeeacDnpOT47/61a+O+y95zpw53q9fPwd84MCBvmjRopj7+rd/+zcHvH///r5nz572KL/T2LZtmxcVFTng9957b0Lbfvzxx/6tb30retZ6rMc555zjf/rTn5LSQ5K2Q5Pjkq4aGhqik6+5ubk+a9asJtetqKjws88+2wHPy8vzxx9//KjXZ82a5VlZWQ743LlzA668c/rrX/8abeNj2z+WTZs2+T/+4z96bm6uA56VleUXX3yxP/DAA/7EE0/4I4884lOmTIme9Q74xIkTffHixe3w00hzFBySdLW1tT5z5kz/+te/7kOGDPEBAwb4pEmT/P7770/o0t633nprNAheeumluN7329/+dvRD5itf+Yo//vjjfsstt7iZxXVpEWmbRx99NNr+U6dOjTnvtHHjRr/11lujl203M58yZYp/+OGHMfe5d+9ef/TRR71v377R9W+++ea0Pm+mo1NwSFK98sor0XMDYj2ys7P9pptuavYw2Lq6Or/pppuih8zOmTMnoRp+97vfeZcuXY5636ysLP/Rj36koY52EDlJkvChvv/+7//ur7/+uj/55JN+xRVXeHZ2djQArr76al+9enVc+929e7d/73vfiw5p9e3b15988sm4DoqQ5GprcGhyvANwd9577z2WLl1KTU0N3bt3Z/To0Zx++ulxT3Ju2bKF22+/PXpJi+HDh3Prrbdy/vnnc8IJJ7By5UpmzpzJ888/T0NDA127duXmm2/mrrvuon///tH9bNiwgRtuuIH58+eTn5/Pn//8ZyZPnpzwz1RZWcmMGTNYtWoVRUVFTJkyJa7bv0pyzJ49mxtvvDHmodLZ2dlcddVV3HPPPYwePTrhfa9atYpvf/vb0etknXvuuUybNi16f/V4NTQ08N577/HGG29QXl7OmjVrqKys5ODBg+Tm5lJSUsLJJ5/MmWeeycSJEznrrLPS6l4rqaTJ8U5s//79Pn36dB85cmTM3kG/fv38pptu8oULFzb5n/rhw4d92rRp3rNnz+hhrlOnTvW6urqY669evTp61jbhYaiLLrrIb7nlFr/00kuj/00WFxcfdWkMyTx79+716dOn+1VXXeUXXHCBX3HFFf7LX/7SN27c2OZ9NzQ0+IwZM7y4uDjak73ttttaPMR669atPmPGDP/GN74RHfqK99GvXz+/4YYbfM6cOV5bW9vmnyGToR5HZvU46urqmDt3LrNnz6a8vJxNmzZRX19P//79Of3005k0aRKTJ0+mrKysyX1UVlYybdo0HnvsMXbs2AFAcXEx559/Pv369aO6upq3336bDRs2RLcZOnQo3/zmN/nSl75EWVkZO3bsYN68eUybNi16U59LLrmE3/zmNwwaNKjFn2PZsmX85Cc/4YUXXqDx71BWVhbf+MY3+NnPfkZJSUkrW0k6i127dvHDH/6QadOm0dDQgJlx7rnnMmHCBE466SS6dOlCdXU1H374IUuWLGHFihVH/b6VlZVxwQUXMHHiREaMGMHAgQPp2rUrtbW1bNmyhdWrV/Puu+/y6quvHvX3UFhYyFe+8hWuuOIKvvjFL7Z4YcjmuDs7duxg7dq1rF27lo0bN3L48GHMjMLCQoqLiyktLWXEiBEUFxdjZm1qs2Roa4+j0wTH4cOHqaiooLKykqysLHr16sXgwYPJyclpcVt354MPPmDu3Lm8+uqrrF69ml27dlFQUMDAgQM566yzmDhxIuedd17MD0t3Z+HChcyYMYM//vGP7Ny5s8X3PO2007j00ku58MILKSsr48CBA6xcuZIXXniBv/zlL9TX10fagDvuuIMrrriCvLy8o95z5cqVPPPMMzz99NNUVlY2+V4nnngijzzyCJdffnnCv9Rbt27ljTfeYMuWLfTt25dJkyYdNXQlEo/IPyKzZ8+mrq6uyfXy8vKYNGkSl1xyCZdccgmnnHJKXL+z7s7777/PrFmzeOGFF1i5cmX0te7du/N3f/d3TJ48mdNPP53hw4cfN6TV0NDAtm3b2LBhA2vXrmXdunXRoFi7di27d++O6+fs06cPo0ePZuzYsYwbN45x48YxYsQIcnNz49o+UQcPHqS+vh4zIysri/z8fMyscwfH4MGD/cUXX2T48OHHNfyhQ4dYsmQJb7zxBm+++SYLFy487rahXbp0YfTo0Zx55pnRR//+/cnOzmbz5s2Ul5czf/585s6dS0VFRVw1nXLKKZxzzjkMGDAAgPXr17NgwQI+/fTT6Dpjxozhmmuu4bzzzuOUU04hKyuLjRs3snDhQl5//XX++te/sn///ibfIzs7m69+9avccccdTJgwocU/nPr6eubNm8fTTz9NeXk5VVVVFBQUMH78eL72ta9x5ZVXBvaLK5KInTt3smDBAsrLy9m6dSsHDhygd+/enHTSSZx++umMGzcuKZeHWbNmDc8//zzPP/88y5YtO+q1nJwc+vXrR2FhIQ0NDRw8eJDKyspmA62goIBhw4YxbNgwBg8eTF5eHg0NDezatYuqqioqKipYtWpVzIDJz89n9OjR0SAZN24co0aNoqCgIObf9qFDh6iurqaqqootW7ZQWVlJZWVlzOd79+49atvs7Gx69uzJjh07Om9whA/RJC8vj1GjRlFWVhb90F+9evVxjda/f3/Kyspwd6qqqti0aVPc79WvXz8uvvhiJk+ezDnnnEOvXr3Ys2cP69atY+HChbz11lv8z//8D/v27Yu5fVlZGVOmTOEf/uEfWpxQPHjwIG+99RZz5sxh8eLFVFVV0bVrV4YMGcKkSZO45pprNAwkkiTr169n1qxZvP322yxZsoSNGzfGXC8SYMOGDWPo0KHRoBg2bBhFRUUt/gPn7mzevJmVK1eybNmy6GP9+vUx18/Ly6NPnz7k5+dH5xZqamrYtWtX3D9bbm4uOTk5uDsNDQ2Nw6/zBkevXr28T58+TTb8qaeeyvnnn8+kSZM477zzor2AiF27drFs2TIWLVrEokWLWLp0KTt27ODQoUOUlJQwcuRIzjvvPC666CLGjh1LVlbz14Q8dOgQS5cuZfny5VRVVZGVlUVpaSlnn302I0aMaHF7EUm92tpatm7dSk1NDdnZ2eTm5lJaWkq3bt0Ceb/du3cfFyZr1qxp8gZbWVlZ9O7dm6KiIkpLS6OPkpKS454XFhYeFWh1dXXs2rWL4uLijhUcZjYZ+BWQDfyHu/+0qXUjcxw1NTW89957bN++ncOHD1NSUsLQoUMpLi5ut7pFRJJp//790X9kw/MS9OjRg549e7b5n9C2znG0PDPcjswsG5gGXARUAIvNbLa7r25uux49enDuuee2R4kiIu2iW7dugfVy2irdxk7OBNa5+9/cvQ54FrgsxTWJiEgjadXjAAYAnzT6vgI4q/EKZnYjcGP421oze7+dakt3fYHtqS4iTagtjlBbHKG2OOLUtmycbsER67CEoyZh3H06MB3AzMrbMk7XkagtjlBbHKG2OEJtcYSZtenM6XQbqqoABjb6vgzYnKJaREQkhnQLjsXAMDM72czygKuB2SmuSUREGkmroSp3P2xmtwKvEjoc90l3X9XMJtPbp7KMoLY4Qm1xhNriCLXFEW1qi7Q7j0NERNJbug1ViYhImlNwiIhIQtI6OMxsoJnNN7MPzGyVmd0WXt7bzF43s7Xhr70abXOPma0zs4/M7OLUVZ88zbTDz8zsQzNbaWazzKxno206XDtA023R6PXvmZmbWd9GyzpdW5jZd8I/7yozm9poeadqCzMba2bvmNlyMys3szMbbdNR26KLmS0ysxXhtvg/4eXJ+9xsy12ggn4ApcBnw88LgDXASGAqcHd4+d3AQ+HnI4EVQD5wMrAeyE71zxFgO3wRyAkvf6ijt0NzbRH+fiChAys2An07a1sAFwB/BfLDr/XrxG3xGnBJePmlwBudoC0M6B5+ngu8C5ydzM/NtO5xuHuluy8NP98DfEDo7PLLgKfCqz0FfDX8/DLgWXevdfcNwDpClzHJaE21g7u/5u6Hw6u9Q+i8F+ig7QDN/k4APAJ8n6NPGu2MbXEz8FN3rw2/tjW8SWdsCwd6hFcr5Mh5YR25LdzdI/eUyA0/nCR+bqZ1cDRmZoOAcYTSs9jdKyH0CwP0C68W65IlR19LPcMd0w6NXQ+8En7e4dsBjm4LM/sK8Km7rzhmtU7XFsApwOfM7F0ze9PMzgiv1hnb4nbgZ2b2CfBz4J7wah26Lcws28yWA1uB1909qZ+bGREcZtYdeAG43d1rmls1xrIOc7xxU+1gZvcCh4FnIotibN5h2gGObgtCP/u9wA9jrRpjWYdti/DvRQ7Qi9DwxD8Bz1nopgydsS1uBu5w94HAHcATkVVjbN5h2sLd6919LKFRiDPNrLm7xyXcFmkfHGaWS+gX4Rl3/1N4cZWZlYZfLyWUqtCBL1nSRDtgZtcBXwameHjAkg7cDhCzLYYQGptdYWYfE/p5l5pZCZ2vLSD0M/8pPGSxCGggdIG/ztgW1wGR53/kyBBMh26LCHffBbwBTCaZn5upnsiJY5Lnv4BfHrP8Zxw9yTM1/HwUR0/y/I0OMOHVTDtMBlYDRccs75Dt0FxbHLPOxxyZHO90bQH8I3B/+PkphIYhrJO2xQfApPDzLxC6ZWpH/70oAnqGn3cF/pvQP5dJ+9xM+Q/ZQgNMJNRlWgksDz8uBfoA84C14a+9G21zL6GjAj4ifDRFpj+aaYd14Q+FyLLHOnI7NNcWx6wTDY7O2BZAHvA08D6wFPh8J26LicCS8Afju8DpnaAtxgDLwm3xPvDD8PKkfW7qkiMiIpKQtJ/jEBGR9KLgEBGRhCg4REQkIQoOERFJiIJDREQSouAQEZGEKDhERCQhCg6RBJlZ1/DFA7MD2Heemb1lZjnJ3rdIsig4RBJ3PaFrQdUne8fuXkforN6rkr1vkWRRcIg0wcz+X/jOccvN7KCZXRl+aQrwYqP1/mhm/9fMFpjZRjObaGb/ZWZrzOyJRNcD/hx+D5G0pEuOiLTAzG4mdFe9a4BsYJO7lzR6/UNgurs/bGb3E+otTAJ2AFVAibvXJrBeNrDF3Yva7YcUSYB6HCLNMLNvApcQumx9PaHLk+9q9HoXoCfwy/CiA8ATHrojXR2wH6iLdz0I3UshvE1BkD+bSGspOESaEB6amgL8vbsfCi8+AHRptNooYKm7N4S/P43w3RnNrAzY7KFufbzrReQDB5P/U4m0nYJDJAYz+zLwbeBr7h79AHf3nUB2uAcB8BlCl+yOGEPoctYQCoeVCa6HmfUBtjUKK5G0ouAQie0pQndCWxieHL+h0WuvEbrPA4QCYTlEh626hsMFjg6HeNeD0HzKy8n8YUSSSZPjIgkys3HAne5+bUD7/xNwj7t/FMT+RdpKPQ6RBLn7MmB+UCcAAn9WaEg6U49DREQSoh6HiIgkRMEhIiIJUXCIiEhCFBwiIpIQBYeIiCREwSEiIgn5/5CEdxGULWEAAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u_z.draw(z_scale='mm')\n" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEJCAYAAACDscAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABUv0lEQVR4nO2dd3gVRffHvycJIUjvRUoApYmCBESQpkiVImJDEMRCERUQVPAVRVFEUOyKvC/8KCogVSkCASGAEHrvHekIhIQWSHJ+f9zC3bl77+7euzeFnM/z3CeZ3SlnZ3fnzJw5M0vMDEEQBEEwS1hGCyAIgiBkLURxCIIgCJYQxSEIgiBYQhSHIAiCYAlRHIIgCIIlIjJagGAoUqQIR0dHZ7QYgiAIWYqNGzf+y8xFA02fpRVHdHQ0NmzYkNFiCIIgZCmI6Ggw6cVUJQiCIFgiS484BCGzkZCQgGnTpuHChQto1qwZateundEiCYLtyIhDEGxi9erVqFSpEnr16oV3330XderUwRtvvIG0tLSMFk0QbEUUhyDYwMGDB/HYY4/h3LlzqF+/Pnr27InIyEh8++23GDp0aEaLJwi2IopDEIIkLS0NXbp0QUJCAtq2bYu4uDiMGTMGv//+O8LCwjBs2DCsWrUqo8UUBNsQxSEIQfLrr78iPj4eJUqUwOTJkxER4Zg6bNmyJd59910AQJ8+fZCampqRYgqCbaS74iCi8UR0loh2eBwbSkQniGiL89c6veUShEBITU11m6I+/fRT5M+fX3N+8ODBKFu2LLZt24ZZs2ZlgISCYD8ZMeKYAKClzvEvmbmm87cgnWUShICYM2cODh48iAoVKuD555/3On/HHXdg8ODBAIDhw4dDPmMg3A6ku+Jg5hUALqR3uYIQCr799lsAQP/+/REeHq4b54UXXkDx4sWxZcsW/P333+kpniCEhMw0x/EaEW1zmrIK+opERD2IaAMRbTh37lx6yicIGo4cOYK4uDjkypULXbt29RkvKioK3bt3BwCMHz8+vcQThJCRWRTHjwAqAqgJ4BSAL3xFZOaxzFybmWsXLRrwViuCEDQ///wzAKBDhw7Ily+f37guxfHbb78hKSkp5LIJQijJFIqDmc8wcyozpwH4L4AHMlomQfAHM2PSpEkA4He04aJSpUpo2LAhrly5ghkzZoRaPEEIKZlCcRBRSY9gBwA7fMUVhMzA9u3bsX//fhQrVgxNmzY1laZLly4AIIpDyPJkhDvuFABrAFQmouNE9BKAkUS0nYi2AXgYQP/0lksQrDB37lwAQJs2bdzrNox4/PHHERYWhtjYWFy6dCmU4glCSMkIr6pOzFySmXMwc2lmHsfMzzPzvcx8HzO3Y+ZT6S2XIFjBpTjatm1rOk2xYsXQsGFD3Lx5E/PmzQuVaIIQcjKFqUoQshJnzpzBunXrkDNnTjRr1sxS2o4dOwKALAYUsjSiOATBIgsWLAAz45FHHkHu3LktpW3Xrh0AYMmSJbh582YoxBOEkCOKQxAsMn/+fACO+Q2rlCtXDlWqVEFiYiLWrl1rt2iCkC6I4hAEC6SlpWH58uUAgObNmweUR4sWLQAACxcutEssQUhXRHEIggV27NiB8+fPo3Tp0qhYsWJAebRs6diqTRSHkFURxSEIFli2bBkA4OGHHwYRBZRH48aNERUVhY0bN+Ls2bN2iicI6YIoDkGwgKfiCJRcuXKhUaNGAByT5IKQ1RDFIQgmSU1NRVxcHIDgFAcA92pz13yJIGQlRHEIgkm2bt2KhIQEREdHIzo6Oqi8mjRpAkAUh5A1CUpxEFFuItL/CIEg3Ga4GvlgRxsAUKtWLeTJkwf79+/HyZMng85PENITS4qDiMKI6Dkimk9EZwHsAXCKiHYS0Sgiujs0YgpCxrNixQoAjsntYImIiECDBg0AwG3+EoSsgtURxzI4vpsxGEAJZi7DzMUANAQQD2AEEXWxWUZByHCYGfHx8QCAhx56yJY8xVwlZFXMbet5i0eZ2WufBGa+AGAmgJlElMMWyQQhE3H06FGcOXMGhQsXDnj9hopr5CIjDiGrYWnEoac0AokjCFkN12jjwQcfDHj9hkpMTAxy586NvXv34tQp2RBayDqIV5UgmMBTcdhFjhw53PMcrvkTQcgKiOIQBBOEQnEAcCuOv//+29Z8BSGU2KY4iKiEXXkJQmYiOTkZmzdvBhGhTp06tubtmmgXxSFkJewccYyzMS9ByDRs3rwZN27cQLVq1ZA/f35b837ggQcQHh6OrVu34vLly7bmLQihwjbFwcyP2ZWXIGQmQmWmAoDcuXOjZs2aSE1Nle9zCFmGgBQHEa0nonFE1I+IHiGionYLJgiZhVAqDkDMVULWI9ARR3sA0wFEAugF4AgRHTWTkIjGE9FZItrhcawQEcUS0X7n34IByiUItiOKQxC0BKQ4mPkkMy9k5pHM/DSA2gD+ZzL5BAAtlWODACxl5rsBLHWGBSHDOXXqFI4ePYq8efOiatWqISnDpTji4+ORmpoakjIEwU4CNVWV9Qwz824A95hJy8wrAFxQDrcHMNH5/0QAjwcilyDYjWvewTWJHQruvPNOlCtXDomJidi5c2dIyhAEOwnUVDWNiI4T0Uoi+oGIRgOoEoQcxZn5FAA4/xYLIi9BsI1Qm6lciLlKyEoEaqqqx8ylAXQHEAtgJ4A2dgrmCyLqQUQbiGjDuXPn0qNIIRsjikMQvLG6yaEGZj4A4IANcpwhopLMfIqISgLw+SFmZh4LYCwA1K5dm20oWxB0SUlJwfr16wEAdevWDWlZojiErERm2XLkDwDdnP93A/B7BsoiCACAHTt24OrVq6hYsSKKFg2tx3n16tWRL18+HDlyRD7sJGR60n3LESKaAmANgMrOeZKXAIwA0IyI9gNo5gwLQoaSXmYqAAgPD3eXI6MOIbOT7luOMHMnZi7JzDmYuTQzj2Pm88zclJnvdv5Vva4EId1JT8UBiLlKyDrIliOC4ANRHIKgT6DrOBoR0Qrnt8Z/JaIH7BZMEDKSCxcuYO/evYiKisJ9992XLmXWrVsX4eHh2Lx5M65cuZIuZQpCIAQ64hgPYBiAJgAmAfiKiJ62SyhByGjWrVsHAKhVqxYiIyPTpcw8efKgRo0asuGhkOkJVHH8y8yxzHyOmRcCaA7gfRvlEoQMZfXq1QCAevXqpWu5DRs2BACsXLkyXcsVBCtYUhxENImI+gFYRUTvE5FrHUgygOt2CycIGcWaNWsApL/iaNSoEQD5lKyQubE64hgHgAEUgmM/qQNEtATAHgAL7BVNEDIGT1NRRo041qxZgxs3bqRr2YJgFkuKg5njAHwL4AAz1wJQEUB/AB8AyGW/eIKQ/uzatQtJSUkoV64cSpUqla5lFy1aFFWrVsW1a9ewcePGdC1bEMxieY6DmdMANHX+n8rM25n5Z2Z+y3bpBCEDyCgzlQsxVwmZnUAnx7cQ0QdElFm2LBEE2xDFIQj+CXSTwzIA7gXQm4jWAtgGYBszT7dNMkHIIDJacbjmOVatWoXU1NSQfQdEEAIl0G3Vn2bmqgDKAfgQjh1yZRGgkOXxXPhXo0aNDJGhTJkyiI6ORmJiIrZt25YhMgiCP4IyNTFzMjNvYuaJMsch3A641k/UqVMn3Rb+6eEyV8XFxWWYDILgi0C3HClMRL2JqDsRPUBE4lEl3BYsW7YMAPDwww9nqByu8pcsWZKhcgiCHoGOOGYDKApgOIBRAC4R0R7bpBKEDGL58uUAgCZNmmSoHC1atADgUGTXr8vaWiFzEajiyMvMHwE4w8yNAXQC8H/2iSUI6c/58+exdetW5MyZM8Mmxl2ULFkSNWrUwNWrV7Fq1aoMlUUQVAJVHK4uUDIR5WLmmQBa2ySTIGQILvfXBx98EFFRURksDdCyZUsAwKJFizJYEkHQEqji+JyICgGYBmA8Eb0O4E77xBKE9Gfx4sUAMn5+w4VLcSxcuDCDJREELYG6485k5gvMPBqOParKAGhvq2SCYJFz587hnXfeQbVq1VCmTBl06NDBtJmHmTF//nwAQOvWmWPwXL9+feTJkwc7duzAsWPHMlocQXBjaQEgEREzs+cxZp5sFEcQQs2qVavw5JNP4syZM+5jx48fx5w5c/DWW2/h008/9buQbseOHfjnn39QvHhxxMTEpIfIhkRGRqJly5aYMWMGpk+fjgEDBphOm5iYiLlz52L79u0IDw9HrVq10KpVK9xxxx0hlFjINjCz6R+A5QBeB1BWOR4J4BEAEwG8YCXPYH4xMTEsCKtWreI77riDAXDjxo05Li6O9+zZw0OGDOHw8HAGwL169eK0tDSfeQwbNowBcPfu3dNRcmNmzpzJANjss56amsqjR4/mfPnyMRw7Wbt/xYoV43HjxvmtByF7AGADB9H2WlUcUQBeBfA3gJMAdgE4DOAogP8CqBmMMFZ/ojiEf/75h4sVK8YAuFu3bpySkqI5Hxsbyzlz5mQAPHToUN080tLSuHLlygyAFyxYkB5im+bq1aucN29eBsD79u3zG/fSpUv82GOPuRVF/fr1+cMPP+QhQ4bw/fff7z7eqVMnvnbtWjpdgZAZSVfFoUkI5ABQEkCBYARQ8jwCYDuALWYuTBRH9iYlJYUfeughBsBNmzblmzdv6sabM2cOh4WFMQCeMmWK1/l169YxAC5evLjPPDKSrl27MgB+9913fcZJSEjgBx98kAFwoUKF+Pfff9ecT0tL40mTJnGePHkYADdr1oyvXLkSatGFTEqGKY5Q/JyKo4jZ+KI4sjdff/01A+CSJUvyuXPn/Mb96quvGADnzJmT4+PjNedcDXP//v1DKW7A/P33326FcPnyZa/znkqjXLlyfPDgQZ95bdu2zT1Ca9u2rdcITcgeZBrFYcfIQxSHYJbDhw+75zXmzJljGD8tLY179uzpHlkcPXqUmZl37tzJERERHB4ezocOHQq12AHjUgwjRozQHL948SI/8MADbqVx+PBhw7x27drFBQsWZADcp08fmfPIhqS74gCQG46dcLsD+ALAQgAnAFwMRhBn3ocBbAKwEUAPH3F6ANgAYEPZsmVtr1Ah85OWlsbNmjVjAPz000+bTnfjxg1u2rQpA+AqVarw9OnT+d5772UA3KNHjxBKHDwLFy5kAJwrVy7eu3cvMzOfOnWKa9euzQA4OjralNJwsWLFCo6MjGQA/N1334VIaiGzkt6T40ecSiIWjk/IvuqcJC8WjBAe+Zdy/i0GYCuARv7iy4gjezJ+/Hi36ebMmTOW0l64cIHvuecejbdR5cqV+eLFi6ER1kaee+45BsAlSpTg3r17c5EiRRgAV6hQwT2CssIvv/zCADg8PJyXLFkSAomFzEp6K47vnJPXT3scOxyMAH7KGgpgoL84ojiyHqtWreIBAwZwx44duXv37jx+/HhOTEw0nf7UqVNcoEABBsCTJk0KSIbExEQePHgwN2jQgHv16sWnT58OKJ/0JjExkRs3bqxRek2bNuVTp04FnOfgwYMZABcsWNDQa0u4fcgIU1U0HOs1NgBoBeBQMAJ45Jsbjs0TXf+vBtDSXxpRHFmHY8eOuc1L6q9o0aL8008/mbK1d+zYkQFwq1atsqVtPiUlhWfPns0jRozg2NjYoOsgNTWV27dvn6VGXkLwBKs4yJGHdYjoHgAfA2gBoDUzLw8oo1v5VYBju3bAsaL9V2b+xF+a2rVr84YNG4IpVkgH9uzZg2bNmuH48ePInz8/evTogZiYGJw5cwa//vor1q5dCwBo164d/u///g+FChXSzefXX39F586dkSdPHuzcuRNly5ZNz8u4bUlKSsJDDz2E7du3o0WLFpg3bx4iIqx9VXrHjh2YMWMG4uLicPjwYSQnJ6NQoUKoWbMmHnvsMTz++OOyaj0TQUQbmbl2wBkEo3WcSqcugCUAVgSbl9WfjDgyP+vXr3fb4uvXr+/lNpuWlsZTpkxxm5/Kli3La9as8cpn69atbi+qMWPGpJf42YbDhw+775MVt+StW7fyE088oTuS9Pzly5eP33//fRnRZBKQidxxH7UrL7M/URyhIyUlhTdt2sSzZ8/mhQsXBjQP8Ndff7kXnLVq1crvgrPDhw9z3bp1GQBHRETwiBEj+Pr168zMvHr1ai5atCgD4K5du2ZLE1V6EBcXxzly5GAA/OWXX/qNu3XrVrfZEABHRUXxK6+8wnPmzOEDBw7w8ePHecOGDfz111+77yuc255MmzZN7mEGk2kUR0b8RHHYz40bN3j48OF85513evUaGzRowLNmzTL10s+ePdu91cdzzz3HN27cMEyTnJzMb775pru8QoUKabbKaN68OV+9etWOyxR88L///c9d3z179uRLly5pzsfHx2sURs6cOblv37588uRJv/muXLmSGzRo4E7Xtm1bPn78eCgvRfCDKA7BNnbv3s0xMTHul7tcuXLcpk0bbtSokdtMBIAbNWrEGzdu9JnP2LFj3Vt89OnTh1NTUy3JsWDBAvf6CgAcGRnJ77zzjnsEIoSWSZMmudd45MuXjzt27MjPP/88V61aVaMw3njjDT5x4oTpfFNTU3ns2LGcP39+BsAFChTgCRMm2DL6SElJ4UOHDnF8fDyvXLmSt23bxhcuXAg639sVURwCb968mT/66CN+5pln+Mknn+S+ffvyjBkzTDe0aWlp/O2333JUVJRbYSxcuFDzQicmJvI333zjtoMD4Oeff56PHTvmjpOUlMR9+vRxn3///fcDbhTS0tL44MGDvHr1ak5ISAgoDyFwNm/ezI0aNfIadRYoUIDfeecdSwpD5cSJE9ymTRt3nm3atAkovyNHjvCXX37JDRs2dCs69VexYkXu2rUrz549WzZ29EAURzZm/fr17k3+9H5FihThQYMG+fXzP3HiBLdo0cKdplu3bl7mCU8uXrzIAwYMcL+oERER3LhxY27fvr17gjtHjhw8bty4UFyykM7s3r2bJ0+ezOPGjeOVK1dycnKyLfmmpaXxhAkTNKOPH374wbCzs3fvXh4+fLh7xbznr1SpUlynTh2uV68eV6lShXPlyuU1Qf/GG2/IehUWxWFHHaYraWlp/Ndff3GvXr04JiaGS5UqxXfffTc/8cQTPGbMGFNeJ2fPnuWXX36Zici9eKt37948ceJEnjZtGn/44Ydcs2ZNjVmhR48evH//fnce165d42+//db94hYuXJhnzJhh+joOHjzIzzzzjPt7F65fvXr1eMOGDYFUjZANOX78OLdu3dr9/JQuXZoHDRrECxcu5F27dvHWrVt51qxZ3L9/f42pDADnzp2bn3rqKZ4yZYrue3Pz5k3etGkTjxgxgmvVqqVJ26pVK69RdTBcu3aNd+/ezYsXL+Y5c+bwb7/9xvPnz+d169ZlygWmojjSkQMHDvBPP/3EAwcO5D59+vAnn3zCf/31l6mJ37S0NF60aJHGw0TvFxUVxd26deO1a9d6PdTJycn8xRdfuBv7HDly8FtvvaW78jotLY1Xr17NHTp0cCsYIuJ77rmHGzZs6M7DZSowmtz0xdmzZ3n+/Pk8Y8YM3rVrV0B5CNmbtLQ0/u2337y2gtH7FShQgLt27cpz5syx7CixefNmfvHFF91OG4Bjz7Iff/xRd9dhf/Lu37+fJ0yYwC+//DJXrVrV/Y75+hUrVoxbtGjBn332GW/evDnDvcqCVRwBLwDMDJhdAJiWloZly5YhPj4eZ86cQe7cuVG1alU0adLEcBEZM2PJkiX44osvsGjRIt04+fPnR4cOHdCpUyc88sgjmsVTzIylS5fiww8/dH//umjRonjllVfQqlUrlC1bFomJiYiPj8fUqVOxdOlSd9patWqhZcuWKFWqFPbv349p06bh9OnTAIAWLVrg66+/RuXKlQ2vf+/evRg1ahR+/vlnJCcnu4/ff//9eO+999ChQwcQkWE+ghBK0tLSsGLFCsyZMwebNm3C2bNnERYWhvLly6NGjRpo2bIl6tWrhxw5cgRVzr///ouxY8fihx9+wIkTJwAAefPmRZs2bdC2bVvUqVMH0dHR7vf4/Pnz2Lt3L9avX4+VK1di1apVmk8UA0B4eDjKlSuHcuXKIU+ePIiMjMTly5dx7tw5HDx4EJcuXdLEL168OB577DG0adMGjz76KPLmzRvUNXnCzDh9+jQOHjyIq1evIjk5GUSEAgUKoGDBgihZsiQKFy4c1ALALK047rvvPt62bZvP8//88w8mTJiAcePG4ejRo7px6tWrh06dOuHJJ59EyZIl3cevXr2KKVOm4JtvvoGrjDvuuAOtW7dGTEwMcuXKhcOHDyM2Nha7du1ypytWrBg6duyIu+++G2fPnsW8efOwY8cOAEChQoXw9ttvo0+fPsiTJ4+uPAcPHsRPP/2EcePG4cKFC17n77nnHowcORKtWrWy3Nhfv34d27Ztw5UrV1CxYkVZeS1ka27evIlZs2bh66+/xpo1a7zO582bF9evX8fNmze9zhUpUgQNGjTAQw89hAYNGqBWrVqIjIzULYeZcfToUcTHxyM2NhaLFi1yKyzA8W35Jk2aoE2bNmjTpg3Kly9vSv6UlBQcPnwYe/bswe7du92/PXv2eCkqHbKv4iAirlmzJtq3b4+2bduiYsWKSEhIwMqVKzFlyhQsWrQIaWlpAIBy5crhiSeecPfwN23ahNjYWFy9etWVF2rWrIno6GgkJCRg3bp1uHLlCgCgRIkSeP3119GrVy/d7TD27t2LKVOmYMqUKdi3b5/X+cKFC2PAgAHo06cP8uXLZ+rarl27hkWLFmHdunW4ePEiSpUqhaZNm6JevXoyOhAEmzl48CBmzZqFuLg4bNq0CadPn4arbcyTJw8qVaqEe++9Fw0aNEDDhg1RqVKlgN9DZsaOHTswb948zJs3D2vWrIFnO3z33XejRo0aqFq1KooVK4a8efMiOTkZV65cwfHjx3Ho0CEcOHAA+/btw40bN3TLKFCgACpVqoR8+fIhZ86cSEtLQ0JCAi5evIjjx4/j8uXL2VdxhIeHs0sx6BEZGYkOHTrgpZdeQtOmTREWFqY5f/nyZfzxxx/47bff8Oeff3rdhLp166JPnz54+umnkTNnTkN5mBmbN2/GggULcO7cOeTLlw/16tXDo48+6rM3IghC5iMlJQWXL19Grly5EBkZGdLO2rlz57Bw4ULMnTsXixYtQmJioum0pUuXRtWqVVG1alVUqVLF/X+xYsV8yszMCAsLy76KIyYmhj/++GPMmTMHS5cuxZkzZ5AvXz5Ur14d7dq1wzPPPIMiRYqYyispKQnbtm3DyZMnkTt3btSqVQslSpQI8RUIgiDc4saNG9i5cyd27tyJvXv34vz580hKSkJUVBRy5cqFUqVKoUKFCqhQoQIqV64c8NxIsJscZmnFIbvjCoIgWCdYxRFmHEUQBEEQbiGKQxAEQbCEKA5BEATBEqI4BEEQBEuI4hAEQRAsIYpDEARBsIQoDkEQBMESmUpxEFFLItpLRAeIaFBGyyMIgiB4k2kUBxGFA/geQCsA1QB0IqJqGSuVIAiCoJJpFAeABwAcYOZDzHwDwFQA7TNYJkEQBEEhwjhKunEngH88wscB1FUjEVEPAD2cwWQi2pEOsmUFigD4N6OFyCRIXdxC6uIWUhe3MP6Qjx8yk+LQ28rRayMtZh4LYCwAENGGYPZbuZ2QuriF1MUtpC5uIXVxCyIKapO/zGSqOg6gjEe4NICTGSSLIAiC4IPMpDjWA7ibiMoTUSSAZwH8kcEyCYIgCAqZxlTFzClE9BqARQDCAYxn5p0GycaGXrIsg9TFLaQubiF1cQupi1sEVRdZ+nscgiAIQvqTmUxVgiAIQhZAFIcgCIJgiUyrOIioDBEtI6LdRLSTiPo6jxciolgi2u/8W9AjzWDndiV7iahFxklvL37qYhQR7SGibUQ0m4gKeKTJVnXhcX4gETERFfE4lu3qgohed17vTiIa6XE8W9UFEdUkongi2kJEG4joAY80t2tdRBHROiLa6qyLD53H7Ws7mTlT/gCUBFDL+X9eAPvg2IpkJIBBzuODAHzm/L8agK0AcgIoD+AggPCMvo4Q10VzABHO459l57pwhsvA4VxxFECR7FoXAB4GsARATue5Ytm4LhYDaOU83hrA8mxQFwQgj/P/HADWAnjQzrYz0444mPkUM29y/p8EYDccq8vbA5jojDYRwOPO/9sDmMrMycx8GMABOLYxyfL4qgtmXszMKc5o8XCsfQGyYV04T38J4G1oF45mx7roDWAEMyc7z511JsmOdcEA8jmj5cettWG3c10wM192BnM4fwwb285Mqzg8IaJoAPfDoTmLM/MpwPGwACjmjKa3ZcmduM1Q6sKTFwH86fw/29UFEbUDcIKZtyrRsl1dAKgEoCERrSWiOCKq44yWHeuiH4BRRPQPgM8BDHZGu63rgojCiWgLgLMAYpnZ1rYz0ysOIsoDYCaAfsyc6C+qzrHbytfYV10Q0X8ApAD4xXVIJ/ltWxdwXPt/ALyvF1Xn2G1bF87nIgJAQTjME28B+I2ICNmzLnoD6M/MZQD0BzDOFVUn+W1TF8ycysw14bBCPEBE1f1Et1wXmVpxEFEOOB6CX5h5lvPwGSIq6TxfEg6NCtzmW5b4qAsQUTcAbQB0ZqfBEtmvLirCYZvdSkRH4LjeTURUAtmvLgDHNc9ymizWAUiDY4O/7FgX3QC4/p+OWyaY27ouXDBzAoDlAFrCzrYzoydyDCZ4JgH4Sjk+CtoJnpHO/++BdoLnEG6vyS69umgJYBeAosrxbFcXSpwjuDU5nu3qAkAvAB85/68EhxmCsmld7AbQxPl/UwAbs8FzURRAAef/uQCshKNzaVvbmeEX6efiG8AxXNoGYIvz1xpAYQBLAex3/i3kkeY/cHgE7IXTk+J2+PmpiwPORsF1bEx2rQsljltxZMe6ABAJ4GcAOwBsAvBINq6LBgA2OhvGtQBiskFd3Adgs7MudgB433nctrZTthwRBEEQLJGp5zgEQRCEzIcoDkEQBMESojgEQRAES2Sa73EEQpEiRTg6OjqjxRAEQchSbNy48V9mLhpo+iytOKKjo7FhQ1CfzhUEQch2ENHRYNKLqUoQBEGwhCgOi1y+fBkXLlwIaRlJSUlITk4OaRmCfVy8eBFXr17NaDFCypEjR5CQkJCuZc6fPx9z584NaRkJCQm4ceOGO3zp0iUcPRpUZ9wyaWlp+PfffzXH1GUS27Ztw6VLl9zhS5cu4eLFi+kinx6iOCySN29eFC5cGNevXw9J/tevX0e+fPlQvHhxS+lOnDiB1NTUkMjkwmjNz/z58zF+/Hh3eNWqVejevbvmgffHokWLUKVKFWzevDkoOdOTa9euoVChQihQoEBGi+Lm8uXL2LFjh235nT17FuXLl0fBggV9xklNTcW+ffsMnxEXf/31F8qUKYPly5cDAG7evIn+/fvjr7/+AuB41tq0aYN27dr5zHPlypWIjo52p7HKxYsXUbBgQVSoUMF9rGDBgoiOjsbp06dN5zN//nxs2rQpIBkAoFWrVihatCi2b98OAFi8eDEKFiyI+fPnAwDWrVuHGjVqoEqVKu40BQoUQKFChdzv/MWLF7F69WrT9R80Gb3KMZhfTEwMp6WlsYsDBw7wn3/+6Q6fPn2av/nmG05MTGS7gGN1Kh87dsx0msWLF/PBgwdNxT1w4IC7DLPExcUxAG7durXpNFaJj4/nYsWK8YwZM9zHrl69qql/l9wnTpzQhPv162eqDFf8u+66y17hTXDz5k33/wcOHODevXv7vceu67Z6v06ePMk//fQTX7t2LSA5r1+/zi+88ALPnTvXZ5y77rqLAfDKlStN5blmzRru2rUrnzt3jpmZz549y926deN169YxM/Py5csNr/HZZ59lADxx4kRTZbryCw8PZ2bmMWPGaMpITU11h1NTU3XziIqK8pLLV1w99K7LFV6yZImpPI4ePWr5fVVxpX/vvfeYmTk8PFyT54gRI3zKefnyZWZmLl26NAPgBQsWmC1zAwezOj2YxBn9i4yM5CZNmnhV5vr165mZuWbNmgyAu3XrZqoyzaAqjtmzZ3Pz5s3533//ZWbmpKQk/v777/n06dPMzLx582avm56SkuKz4di/f7/lB/HFF18M6OFNSkoyHbdcuXKaMs6ePcsAuFmzZu44rvN79uzRhJ9++mlmdjTOsbGx7oddxRW/bNmylq7Dk7Vr13LTpk15165dptO88sorHB4ezqdOnWJm5goVKjAArlu3rm780aNHc968eXnPnj2W71d0dDQD4EGDBumev3r1qt9OxjfffGNYntoQGeGK/8ILLzAzc+fOnTVlLFu2zHSZ9evXt1RmWFgYMzMPGTJEU8bNmzfdYU+l7klERIQmzZkzZzh//vw8YMAAUzLoXZdVxREfH++Vx9y5c3nYsGGaTpU/jBTHp59+aqg4XOE333yTmZl37tzJzZo1c7eFOmUGpTiytKnqxo0b7qGuJzt37gQAbNmyBQDccZgZGzZswOXLl73SWMWxSzXQoUMHLF68GB999BEAoF+/fujTpw9atWqlkcWTKlWqIFeuXLaZuxzPgZYff/wRU6ZM8Znmq6++Qt68efHrr78CABITE9G+fXvMmTPHVJmLFi0CAMTGxnqdc9WNyogRI9CsWTN07NjRVBmBUL9+fSxdutRSGf/973+Rmprqrq9Dhw4B0L93APDmm28iKSkJ77+vt4u7N1u3bnXbo48cOQLAYWbRIyYmBhUrVtSYPm7evOn+X7WFmyE5ORmLFi0yfN6OHz+ukTEYTp48iUcffRR//vmn33i+npW0tDTd/z1Rn/vx48fj0qVL+OKLLyxKa14uIxkAoG3bthgyZAji4+PTrUyVxx9/HLGxsWjUqJGldGbJ0orDKn/88Qfq1KmDhx56yFI6MxWemOj4PMaSJUsAwK+d/sCBAwAQskm4Cxcu4NVXX8Vzzz3nPvbDDz+gefPm7kn3/v37AwD69u0LABg1ahT++OMPdOjQwXZ5XC+Eq2F2KZ1AOHXqFDp16oR169bpnnfZfM+dO2c5b/XF9QyvXbvWa3LYzIu+adMm1KxZU2NH98fu3bsBwN3gTp06FZGRkfjtt98AWHv5XfL17dsXLVu2RO/evU2nDRRXmf3798fSpUvRunXrgPIJRHHoMWDAgHS5bj1cTjTbt29H7969TT+TRtfled5X3DNnzgBwzMG5ZClZsiTeeustUzIYke0UB+DwUHCxbds27N+/3x2eMGECnnrqKXcv7/fff0fRokWxYsWKgMr09xC4XrIDBw7gm2++wY0bN3zGv3DhgqbnefPmTVy5ckU3rp6HT58+fRAbG+s1CnHJoOcxs2/fPnfv087eiidxcXE4fPiwrkyAoyH98MMP3df06quvYurUqahbt66p/BMTE9GxY0f3vQeA2bNnG/aEPYmNjcWDDz6I++67z+ucXr2cO3fOPWJwjSwC9Ujq1KmT5m8gjBvn+HbRhAkT/MYz2+M1g971xsbGat49f2V6Kotgnr3Ro0djzJgxmnfF1bmw45k2837XrFkTY8aMQZ8+fdznUlJSfCUzxJ9SdZWpyjVp0iScOXMGn3/+ecDlenJbKg5fD6NamVeuXEGNGjVQqVIl97Hu3btjxowZmDZtGgDHkO/8+fOGpo9gXrrKlSujb9++GD16tK68p0+fRuHChVG9+q2PeFWqVAl58uTB1atXLb0Avtx81Txu3LiBypUro3z58qbzdmG2/nft2oUmTZp49cY901erVg1Dhw51mwJd5hTPPF9++WWMHTtWN4/PP/8cs2bNQvv27QE4lOoTTzxhqSe8ePFiAMA///xjENNBsWLFEBMTg507dxrem+TkZIwZMwbHjh3Tld9X2JMLFy6gc+fOpjs3zIw5c+aYNkl5XoOrJ5uamop33nlH11SppgEcdde8eXPUqFHDVHzPBnHZsmXuY02aNPE5gjDzHrzzzjvImTOnV2dFTetZNz179sSTTz7pjrN582Z89NFHGjdef7iuxdVB3bJlC3LkyIFPPvlEE89Xfp6y3bhxQ1M3np3J9OS2VBwqvl46f71AtTev5vHzzz/rpgukF+N6EFxzMi5cPZS///4bgGME4ML1YHuOluzE12jGhb/rVHu2v//+u248l1lG5ciRI16uxa75BrXcFStWYNy4cejZs6duXqqvu51u1OozkZSUpAn729XAlfbTTz9F7969ERMTY7l8Vx299957+PXXX9G4cWPdeGqdxcbGokOHDj47Bf7u7eOPPw4AmDFjBkaOHInmzZubkvXkSf0PyqWkpHg1focOHdI0jo899hgAR0cjLi4OY8aM8crHrCv6yJEjkZqaiu+++05z/IMPPtCEX3rpJff/Y8eOxcyZM92mp1q1auGDDz7wysOoEXfd86FDhwJw3DdPPvvsM680qnnrgw8+0NSN6364UJ/BUHFbKg5PswRgz2QfoH2hBg8e7CemfhozN9Uz/o8//uh1zIgTJ074PW+HOULtGanD7uHDh2vCgTTWLp92F77kNlp4Z6Xuli5dqgm77pe/PDzPBTI/tHr1agDmJ709y/vyyy8BOOZ99Bg2bJjucV/zb3qjB3U9g2vC11eZf//9t+XOk1rvL730ku68hpqvZ3jhwoWG8T1RnydfdeWP/fv3a8rw1XmximeeqvPPnDlzNHWzatUqzXl15K2Xpx1YVhxElJuIwm2VIkgmTZqkCc+cOdMrjt7ElJUHC/A9SQcA//d//+d17ObNm5oyevXqZViGVdQ8Bg4caOkhcdWLv5dSDU+ePFlz7j//+Y+psozK8CRYD5NAJsfnzZvndczlyKCH2qiqDaC/a1BfeBUz1//LL78AsLdRUJX8qFGjLOd//vx5r2P+8lBHCwkJCaYUhydqB+Lo0aOGE+zB1pt6j9Q2QDWHBVqGPzOeL8xMoAeDoeIgojAieo6I5hPRWQB7AJwiop1ENIqI7rZdKot069bNMM7JkyeDqkAiMnXDPFHNWepIyDWPEiye16W6Gqvbo+jVgdoAGrkrJyYmavL54YcfLMkIeJuQzGLlHuqZR6w+A/483/Rcwa2U589U+u6772rCwewKYEVhq+UE0rlRG7uUlBRLnkJAYA295/nu3bub8swKJa+99lpI8rV6LcF4MfrCzIhjGYCKAAYDKMHMZZi5GICGAOIBjCCiLrZLFmI85wv00HtIrWp6I0+aIUOGaMJ6vQtfsphFtUGrox7Au2f0/vvvW+6xWG0YXGa4YDDaz+vo0aN+G02z9WrWVGWEaqpMTU21lP7QoUNB9x6vX79uONKzu4f63//+11J8s500TzlVBXfq1ClNHur8YSAEayFwmQjTuwPrcmhwYce8qBnF8SgzD2PmbczslpiZLzDzTGbuCMCernM6MmDAAE1YnQxWG9izZ88G5EJn5SE5ffq04UNhNGmvvvgbN240Xb6L7du3+21grZrx9M4b1Ytaxty5c73SjB071lL9XrlyRRNfz7xohDoXYeUldq2ZCRTVmcCX66U/jDYNtEOZqvdOna+y43kyg2ceTZo08SrTaodIT26rSsCMW7ZRnoEoVc/ws88+a5jeCEPFwcw3nYXHEVE+5/+9iKgfEUV6xsnMGDWwLndNf/hadBZomSrLly83fChUhabniRHIS6c+aJ55TJkyxbChMCqTmf0+zOq8iZmR19mzZzVh12InXwwdOlQjp6fnjC/Ue6g2/v46E5s3b9bIrG7GZ+e6CRdmdkXwdy/11gWo8VVzlt5u0f5Gem+//bZf+dLS0rzkUJ0+FixY4LdMIvJr3jOznuHChQsaOUaNGqU573LTtoI66tRzL/dE7ezs3r07aKUayO4DKlYmxwswcyIRxQB4BUBBANbGoH4gojJEtIyIdjvnT4Lrnung72FWJzb1UN3t9Dx6gjUDqA+FOonq2iLExS+//BLQ5Jkqlz85unTxtkQalak2kpcuXfJ77S+88IIm/OSTT+qW6c88YeR8sGXLlqDrRvXQ86c4vvrqK8P8Pa9HVYRGqIoJcLj3+itDvR7VkUStHz3vKXVbGrXejZ5xtdH+5ptvNGm2bt3q9a598MEHmjguF10XcXFxXuVYvdeqI0SPHj00eaidNH+OE75Q2wCXa64LdWW3+l4A3nOSrt0q0hMriuMmEUUA6ArgM2b+AMA9NsqSAmAAM1cF8CCAPkRULdDMVBNNWlqapUlCPdQH0bUozV9+nscSExMtm3Ti4uIsyblw4UJbhvn+emvqZFtSUpJXmer8TdOmTTVh9QVSGzS9l9LoulTvOsD/SEoP9dqMlJWqOKyY5NRzH3/8sV/Z9Doeak9cdQg4ePCgXxnUtQQbN27UxHd5bnmijmrWrl2rCbsW7JlFb/JW7YkbOQeoayoA4+dFrRfXeikX69evt/wu+Vqv4gv1uszss6WO8NQ5JJebtgsri0nNYkVxfAtgK4A2AFyG0jxBS+CEmU8x8ybn/0kAdgO402x69UFTJ4XVlcKBTDKp8dU9+Pfu3asJu/aD8hdHJdhvaty8edNyHnv27DFsYD3P6/nMq5OP6sOpbjWh5mkGO1wqje55y5YtNeGnn37ab3y1rqdPn25aFlX2QFYBGyk2z03u9M6rYVXB68mpHlPzeOqpp/yWYYZAJpDtNumYNcF6EooPT1m1KLz55puaNOr7bQemFQczTwRQF0B1Zr5GRHcBWGOrNE6IKBrA/QDW6pzrQUQbiEizLHfixImaeOok1KlTp4I26RhVvq/tFDxRTR2qXd5oAl7vJTSyQauoC/i6d+/u9aBZfWFUzw2jNL62qvCH1fulZyqymsf58+f93nfVI+3gwYOasL8GVt1yQkX1fjG7c7Endrgkq5hRBOpI26rXnR2dBF+LFK1g9T1Q60Z1VbdDifran8oX6oJcOzCtOIjoEQBfA/iYiLoDyA/A27czSIgoD4CZAPoxc6J6npnHMnNtZq6tpLNUzh9//GHZbGQU398mhb7yUG2cquIwoyjUsJqHa3WyCyPvnmPHjln+dK3Vl9qqOQNwzAF4lmO0UeG0adMMR1JmPOX8XZvnDsRW+eqrr/wqFtUTyKpbqx6qOSZUeF6Xr+15/GG1cdRDXStkVQkEmsYTdSsTdbI7kOsKdiRlx67cVkxVPwOYB8fajQoA3geg/8GCACGiHHAojV+YeZadeeuhNhrqQ6B6f1jZUdVXnupNV19kvYbMMw8zvv9qHqrcenv9qFjdJsQOs4ARtWtr+grYsGGDpZ7smjVrvOR07RxrlmB661a3vjEzWrDa8Kj2bz3MzAUG2/gZEYiCN8pDdb82M8qx2rk0ol69epbi6xGsGc8OrCiOA8w8m5mnM/MQZm7PzHfZJQg5nr5xAHYz82ij+CbyM4yjmnRUFz91okpdtBYK263qcbJmjbc1UP1AjDpBasd6E/WFMLK/q/HNNHqeYT2Z1WN6boRWXoikpCQvhWi0EFRFHb3ZSSDPU0JCgqYOzMyxWG301Xv7zjvvaMJmFKLVhku9T5MnTw7a3KW6z6or/9Vn9ujRo15thOoIYGSiVetXndg2s5ZIrQt180w964DdikLFiuKII6L+FIruhYOHADwP4BEi2uL8md77WhXLqJcEeCsOq1/sCmQ/JKMbqk6u6e2fpD6sqtnHjm0j1JdOdRU16pXasU+PHd8OUOVSXzK7H2cj910VT/kCce+8//77NWGrX7fctWuXYZw9e/Zowqo7qB5WG66tW7dqwqq7LWDcu1cbWKP4qu1f9QQEvOtH9XxTF9OFonlU5VS3o9Gb+wq14oiwEPceANUBvENEGwFsAbCFmc27kfiBmVcBCLjW1QbBzBAy2L3s9XYZNZpvMNq91gxGD4Ude/QbTbCrPSczJhCrDaPeXl7BvhDBLo60yvfff286biALyoL5IJAvVOXTr18/y3kYmYFV9Bb0qVhVBHY0nqrzipFiUN2b9Qh1ow7Y+/kAPUwrDmZ+AgCIKBduKZEHAdiiOILlww8/NIyj9pxC8dKpeb7xxhua8Ouvv64JW/VQAYxfoM6dOxvmqaJ6oRnVzcCBAzVhq6M1u7DbHVcPK8rejglsf6RHo+PrW+tW5FAXrtoht9V7py7QTY+JaHVkHjoDjTXs3uTRtOIgosIAngZwHY5J8WnMPMFWaUKMahs048FkFfXBMdL8gXx4xeghUNe0mPGQcn3T2kUolKpKekzyqXt7qSMp17fi/aG3Tb8vrJrorF7PTz/9FHQe6aF8VOzee8oM7dq1C7rMYOvK6m4AocJMx9oKVuY4ZgMoCmA4gFEALhGR/ifcsghWV3mGAtW2awarL9DIkSMtl6G6CYcC1VVRRU+RX7p0SRM2UrzqymdVIYZ6hGCEma3ZPQl0O/qMJpA1KCrBKh+jj36Fokw9MkJxB3Lt/rAyx5GXmT8ioieYuTERdQRgm1dVRvDMM89owlZf4oxCHR2EAnWvqVAQyPbOqg/6888/byl9ZukBZiTffvttupf59ddfB52H1bVFKoGs6lYVR7AyADD9bfhMjWvLAqMfgDXOv2sB5HL+H2c2fSh+ADg7/N59990Ml0F+mfdXt27ddC+zevXq6V5mwYIF073Mpk2bZvj9DdFvQzBtL5kdNjlHGMsAvACgDoDVAPrauZbDKkRkTnhBEATBk42s7L5hBdOKQ5OI6HkA9wKYyMy2rh63KIcoDkEQBOsEpTjMfHP8A+ff+kSUFwCYeTIzv52RSkMQBEHIGMxMjrtWJfUDUN35TY5dALYB2MbMM0IkmyAIgpAJsWyqIqKccCwAvBeOLdbfMkgSMsRUJQiCEBBBmaoMRxxE1ArAVmY+CQDMnAxgk/MnCIIgZDPMmKo6AhhGRMUB7IFjj6qtzr+7mTm4T9YJgiAIWQrDyXFmftk5pPkCwD4AhwE8DGAdgOC/CCIIgiBkKaysHO/OzDVcASL6AUCGzW8IgiAIGYOVvaoSiSjGFWDmjQAq2S+SIAiCkJmxMuJ4EcDPRLQLwEY4vKqC//CDIAiCkKUwPeJg5v0A6gNYAKA4gN0ATH+hTxAEQbg9CPR7HHMBbGfma/5TCYIgCLcbwX6PY4//JIIgCMLthhXFkZeZPwJwhpkbA+gE4P/sFIaIWhLRXiI6QESD7MxbEARBsAcrisP1DdRkIsrFzDNh4xwHEYUD+B5AKwDVAHQiomp25S8IgiDYgxWvqs+JqBCAaQDGE9FqAHfaKMsDAA4w8yEAIKKpANrDsaGiIAiCkEkwNeIgojAAlZn5AjOPhsOzqgwcDbtd3AngH4/wcegoJiLqQUQbiGiDjWULgiAIJjE14mDmNCJqCsfEOJh5cghkIb2idWQZC2AsILvjCoIgZARW5ji2ENEHztFHKDgOxyjGRWkAJ0NUliAIghAgVpRAGQDPAjhJRL8T0TAiespGWdYDuJuIyhNRpLOsP2zMXxAEQbAB05PjzPw04PUhp7oAptshCDOnENFrABYBCAcwXj5NKwiCkPkw/AIgEREbRDITJxTIHIcgCEJABPUFQDOmqmVE9DoRlfU8SESRRPQIEU0E0C1QAQRBEISshRlTVUs4dsadQkTlASQAiILDnLQYwJfMvCVUAgqCIAiZC0NTlSYyUQ4ARQBcY+aEUAllFjFVCYIgBERQpiorK8fBzDcBnAq0MEEQBCHrY9odl4juJqLxRPR9KAUSBEEQMjdW1nFMhsP1tiEAEFF1IpoUEqkEQRCETIsVxRHGzH8CSAUAZt4BoHpIpBIEQRAyLVYUx0mnVxUDjrUbAHKFRCpBEAQh02JlcrwfgP8BKEFE3eFw090RCqEEQRCEzItVd9wIAI8DqAGHd9V4Zr7uN1EIEXdcQRCEgAjKHdeS4shsiOIQBEEIiJBvOQIAIKI4Isrn/L8XEfVz7mIrCIIgZCOsTI4XYOZEIooB8AqAggD+GxqxhGrV5HPrQubizTffDHkZ4eHhmvAjjzwS8jJLly4d8jLsoG/fvhktghsriuOmc46jK4DPmPkDOLZXz3aEhYXqW1a3SI8XJhRUqFAh6Dx69+5tgyS3N3/8EfpP1dSvX18TLlasmO1ltGvXThOeO3euJpwjRw7DPPLnzx+UDK+99pomHBmZOQwpP/zwgybctm1bTbhRo0bpKY4GKy3gNwC2AmgDwHV389guURagSpUqltO0adMmBJJoqVmzZsjLuP/++zXhd999VxN2eGlrqVq1qqUy8uTJlo+VJfLmzRvyMgYPHqwJ693bYHn00Uc14dy5c1vO4/PPP7dLnIAZOHBg0HmULFlSE1aVptphDcX9MItpxcHMk+D4cFN1Zr5GRHcBWBMyyUzSq1cvwzitW7e2tUwzN+z777U7s3Tq1Mn2MlS+/fZby2msojYmZkxqRYsWtVRGRjlsqC+qleemT58+dovjFzOj3hIlSljKs0OHDpbLUHnuuecsxVfvtR2jebVnboQdDXAgCs8qZurGjhG/GaxMjhcG8DyAZ4noAQAnmLl7yCQzSUSE8VKUe+4JvUWtTp06mnChQoWCyi+Qh9lMXah88MEHmnBUVJTf+IG82GrjkB695eLFiwedR5MmTUzHDbX5Uu14mCkv2AbRjnudHmWq9O/f31J8tZ4ysifvDzNyqvU/fvz4kMhi5S7NBlAUwHAAowBcIqLdIZEqxKSHDd3oBahYsaLtZaqK48UXX7Scx7x58/yeN3p4zbx0kydPtiyXimp/VzG6x1Z7pXoUKVIk4LRGPfMCBQpowvXq1dOE02OezY4yfv3116DKDKQRT4+6UbFD2Zj4GqvlPFVng9dff90WxxsrNZyXmT8CcIaZGwPoBGBC0BIESSCV2blzZ004FHMDoXh4X3/9db/n1YfETN2oD2uuXP53kbHjuozkUq8DAFq2bKkJDxo0yG8e6nWptnQzox6jF9nKC/jggw+ajqtXdiAmHTXNf//r3wlSjW9HT9yquUzvuu68806/aYzkDmQS2Y4Rq0r79u0txbfDZKvWZ+vWrbFz587g87UQ17VCPJmIcjHzTAD2Th6kE2rDVKtWLcM0Q4cOdf9vZohoR8+pWzftF3nvvfdeTVht/PQaXKuoedx1112asJkG69VXX/V7Xq2Lt956SxMeMGCAVxqjxsGI999/3zCO0Ys6fPhwn/GN5DFTvhX07oPRKMZqnel5Kw0bNkwTfuyxxzRho/dA7fyUL1/eUMaVK1dqwkYj6WB7/0SE1atXa47FxMQElSfgPdc3YsQITfijjz7ShJ955hnLeQb7npjFiuL4nIgKAfgNwHgieh2A/66ASYhoFBHtIaJtRDSbiAoEk59qD1YrUzXpBNIzN5qECuSGefaMiMjrBVHzVL2P7FAcat088MADfmXQ4+uvv9aEjcwslSpV0oTNTKareagulSpmRmM//vijJqze88aNGxvK5Qu1PDvmAl566SXNsR49evhNo5b5xRdfaMKqV4/eM64eM6p39T6pc3+qtyER4YUXXtAcy5kzpyasmgibN2/uVwa9e+2psIgIP//8syasXqeeovb0FgykDalbt64m/Morr2jC+fLl04SJSKN4iQhr167VxFHnWkOFFa+qmcx8gZm/ALAAQBkA1sZevomFw1vrPgD7AAw2iO9G74aNGzdOEy5VqpQmHMgksoraC1JvstHEqplRi1VUryAzD7NqmjKqm6ZNm3qV4dlgEZFXuZ6jNcAec5eR8jHqLOjx8ssva8Lqddjx3Phi+vTpfs+b6Umqx4xGIOpI++OPP9aEw8LCvEYERqhejkb3Ws8c5qnkzTzD5cqVsyChOTnMoF6bkYdnKDzI1LmwMWPG+I1v1wjEilfVI0Q0joi+gGNX3WkADtghBDMvZuYUZzAegOmlnESEJUuWaI6pjaFqOgm2Z05EKFiwoObY2LFjNeG8efN6PRiffPJJUOUaERkZabgGQrXLq26kRo1jrly5vBSUOtGvXrd6P+zw1rH6AgTS6Ks9QH95EJFX58Ef6vWovWj1fNeuXb3KU1Hr1Z9pTS8PdTSgZ2416tyoo7JA5klCMUkcbJlmaNasmSYcHR3tN34o5gvVe2iHiU0PK5L/DGAeHA17BQDvAwh+lsWbFwH86eskEfUgog1EtMF1zGh1qTrUNXPDjOKo5/UaT7XnpDY8RpN+RqgPTVhYGB566CG/aVSvKaN5Er0X3+jFNTpftmxZS/EBbwcGo/ujmiv1zJNGbptq50CtG7WxUVc9q+V5YmS+mzFjhl9Z9FDzUJ97I5n08gu2EbdjfYPVRp2ILE9Eq+nNHDNi2bJlfs8H0oE1cppQCWSxshmsKI4DzDybmacz8xBmbs/Mdxknc0BES4hoh86vvUec/wBIAfCLr3yYeSwz1/bc2dHqg6W+UHoPxV9//aWWqwlbHYLrER8f77MMvd6dmRddJTY2VhM2Mj2EhYUZlmNlUliPu+++23Iadb2JujZHlUN9YfReUqvuymr9vvfee35l8Ic6ClbzVnuvelh9JtUV/IE8T1axug+UXaMJT/MVEXmNeoN9hs3koY44GjRooAlbLTczrTcxfDKIaBIR9QOwhoi83V1MwsyPMnN1nd/vznK6wbGdSWcO8dLhiIgILxc91WSjDrm7dOmiCZtRHP56B0Tk5apYuXJlTVhvnUDDhg19lqlnWrBiPnHl8fjjj2vysIKZ+GFhYYb7C6mjyDvuuEMTLlCggKHrsCd2OA6o19aiRQtN2MpjG4iDhor6vBjloTZcRujVmdXJ17CwMFSvbu0L056jFCIKaI8s9V5s2bJFEzYz+a1i5H5tdP9VL0krz6+LwoULW07jSXrOcUx0xisB4HkiOkZEfxDRMCJ6yg4hiKglgHcAtGPmqxbTWi4vPDzc66apNmQVz16yGXONihmb9KhRozRhvRXvTz2lrXJPFz4zpgUj0mO1sF7PSW0ctm3b5jePsLAwU+Ybz/ieCjEQ/N1zI/u/3b17IvLaLdXqmgkjwsLCNA4DROTlhGAmD0+IyGsdlacyIiIv1+UcOXJolA8RWd7iRfXUmzhxoias92yooz7VCjFkyBBLMqh1oacQ1Q7FihUrNGHVVT0Q5WMHhk8rMy9l5tHM3I2Za8Ixv/EeHN5Pdf0mNs93APICiCWiLUTk3zUgSIIdgptVHEaLiNQ89BpCtWeujlo8zSV61xVIo66WsXz5cr95du/uvfOMPx90vbpTzT5GtllVSZppmD3rx8xEr4qem6aVtG+//bbP84HsyKqaXHPmzOm11sdIJn+EhYUZLrQkIi8nAjUPtZ6feOIJTXj+/PmasN5IW20gVTdez50AzMzLqSY0z/xdcT3n1YjIa77m2Wef1YSN5pRUwsLCvEZw6jyZamFQ5ySjoqIyRHlY8apqRUTxcEyIvwfHnEfwW0ICYOa7mLkMM9d0/ox3LvRAb/sJdddWT6w2OnqYSfPnn9o5fnWdhlEP1uoxPfOXVfTKUIfHnhPsvkwJRh5ldo+MiMhrG2o1vpkyjbbIsLIC3Gj7Fk/MuL0+/PDDfs+HhYUZ7jWm8sYbb2jCnnNiqoOHL/yN/Mx00jxHA4Ha8X/5RTstauX9NluGv2shIq/dDYLJz1cZemFV+fjbYSI9TVUufgAwAMCDAMYCGEVE1rZ8DQFEhBw5cnjZ8tXevuf+SHZtB6LuJqqu/FR7f6onjRksfhNes8W02YdkwYIFltKotn1VBj3U3pmKVUWipwj8edfpjTD0Rg9qHHUbeV9y6nUEjFZVexIeHm64GeesWbM05anoPddGJlhVJk/3zUDfk+3bt2vysHvKUu/a1Z64VZOaGYzqw6rLt13tkFq/Tz75pCYciu92WJH8DDP/zcwXmXkJgBYA/mO7RAFi1OB52lWDNen4KktVJGbSWMWoJ6XacvUaR9WG3KpVK01Yb9huN57+5mbL8FQMZnqmnl4tevFz585tONm4aNEi03L623gxkDkOtSFSF3vp5aE+x1a3OffETB3rKWRPjyY9t3Q7FInRfJK6GC+QMu3YvsOf+SoUSlWPmTNn2p6nFcVxhIg+plvfGb8JIMl2iWzC300nIi87q95LabSDbShvuq+HVN0WwiiN3vYdZcqU8Vv2yJEjDaQLHvULdmbqcvPmze7/icjw41ieq/uJSNfsYuR1ZvZbIhEREYiIiDBs3D3xnJDVUxzTpk0znZevPNR6tboq2wz+3ILDwsIwadIk03n5ksnqu6Z2Mvw14GYVpD985aF+v8YTu0xVKkaLS+3AiuQM4AkA/xDRKjhWjS8nIutO+SHA6oPlOaFHRF6ujYD3DqyuvWXs+Jqf66b76hG6branvz8ReQ1DzZThT9nooU4cBuItZHQ/VDOemUZMtaV/9dVXfuN7XgcR2aIQ1ev64YcfUL16dfdGjf6uWz03YcIE9/+Bmi08N6E0k4e6eV8g+JvkV+UICwvzWrxpZRLZ5SSh3murczmqO7cnLmXviuPqLKjvuZnP2Kqo93z9+vXu/8PCwrw6cf6eH3USXN1KyQzpPsfBzJ2YuRqAcgD6AfgQQG4A/yOif2yRJgACrQi9l8yogf3zzz8xY8YMy254/lAflAULFqBLly5uV0t1N1J/jYO/l8MsevXpbxsTV3xX2S4F/OabbwK4tQmekfJ54YUXULlyZc2+VkZKUvVyMVJWZuY0rHZAevfuje3bt5samagvutrABmIu9XTXDAsLc9eZrw0Z1U0rzZT50Ucf4a677nIvmFR3da1du7Ym7PkM6TkxtGzZEm3btvXaZBG49Q5u27YNgwYNcj8PnjsiEBEaN26MZ555xrDz4InqoTR9+nR07drVvUZrxYoVaNq0qduzSa3DsLAwr3ko11yqkeOCC8/0YWFh+P7779GpUyf3ZoVNmjRBmTJlNGs+xowZg9dffx01atQA4FD+nTt3ds9nuq5L3ZEhlFjewIeZrwPY4PxlOK7eQrly5bB9+3Z3r9T1oAXjqqa+VAULFkTHjh3dYX+rUVXMmjxatWqlmXMw2h8JcLxk169fD4lbHhGhfPnyGDFihO6D6ZJh/fr1+P77793zJ3Xq1MHly5dNbzmRL18+7NmzR3PMaJGgKocVjzJfHQ7XTqt6O64SEXr37o2pU6danjsgIrz66qs4fPiwey5MVRxWcPV+PUdhYWFhGDBgAGrVquX2/jJSDH369MF3333nVvB6z9CQIUM0nSXV7Pv000+Dmd1lqudVIiIivEyVe/bswYULF9weevfeey8+/fRTn3KHhYVh6tSpuufMmguffPJJTeckJibGa987ldKlS2u+Z7Fnzx4cOnTI1KcZAO97XrJkSY0nX65cuXD06FFNvfXs2VOTR7169TSONkOGDEF0dLTbGSM9PmQVuu0+04ECBQq49xuaNWsWBg0a5N6aonbt2li3bp2um6MvW6HV3uZnn32GvXv34p133gGgby9fsGABNmzYgCZNmgBwzJscPHjQ8At2evhq7Ix894PtVQNwX6OLwYMHY+HChe7hfLVq1by+s+6pNFy97WAeaiO5H3vsMQwZMsTv6nqV5s2bY+zYsW5voqJFi+LatWs+zSmNGjXCuXPnfE6q+5IxMjISOXPmxDfffOM+pjYiL7/8Mt58802vRZ6efPrpp9i7d697+xA1j4iICFPblQCO56lKlSpITk52ryOJiorCihUrLJmTiEjjNReIFUDPVGyVmTNnYt++fV4rvO0wz/jKo0CBAhqlMXXqVDz77LM+56fMdBasyhsVFaUxvderVw+NGjXyGl3aSZZWHBUrVnSbUe666y6vjeF8bY9gtIOsWcqUKYONGze6w+XKlcOXX36psVuqI4jly5dj6tSpht9OCIb33nsPw4cP1/0gEuDw/urZs6fmZX/uuedw/Phx03sLDR8+3Gv3VX/kypUL58+ft9QgBeKiq7pEu/D1Mo4ePRoxMTEajzgj+7mVycbRo0dj/fr1uptPqiadvn37omHDhrjvvvsA6F+/Ou9m1BC53gH1g1yeqIsPrSheX/kNGjQooEWNZvB1L1WHFxdWt93xh9Ez+cwzz6Bjx45uS4E/J51QjQwiIiIQFxene87Mly9NlWFLLlmE48ePIzk52XTjZXUlKAD069fP7/nSpUtj4MBb6yZD4Zk1bNgwDB061O1B8+677+KNN95wzzsUKVIEycnJmsk+dQGVCzvdcdUtn62W4c91NhB3V8AxKjJS4u3bt8f+/fv9Nr4uSpcujV27drnNbP524FVNOkSkmS9o2bIlSpcu7XdhmVFDVKRIEVy8eNGWXWqt4MvMFCr3bj3+/PNPjBw50v1hsfRwfQW05uWePXvi888/dzu5BGOeDIbJkydj69atto1CspXi8LWNuctM0alTJ3z55Zdo164dAEfvYfLkye5wKBg4cCCmTp1quMW3JwULFsTFixf9Du893S5fe+01tGnTRrOuwchDpEmTJli+fLmpz1cGists52siV1UUUVFROHv2rG5P1tco8q233sKpU6e8JqatfFN6zpw5YGZTjd7s2bPx1ltv+Rz5eGKUX+7cuXHs2DG/8cw0RL5s/qE04ahUr14dO3bsQOvW6fe16ZYtW2qUrpEbeigoVaoUkpKS3PfGaP4nVHTp0sVro9agYOYs+4uJieFgOHjwIP/vf//jlJQUZma+fv06z5s3j5OSkoLK1yrXrl3ze75Xr15ctGhRTkhIYGbmXbt2cY8ePfjEiRMhkyk5OZn37dsXsvyZHfV96tQpn+cTEhK4Q4cOPH/+fJ9xJk2axM8++yzfuHHDVJlLly7lLl268KVLlyzLGwoGDBjAgwcPDioPOFzlOS0tzVT8UqVKMQD38xQIuXPnZgB85swZU/GvX7/Ohw8fDrg8ZuZPPvmE77zzTtNlqpw5c4a7dOnC8fHxptMMHDiQAfC8efOYmXnChAkMgJs2bRqQDMzMw4cP548//jjg9HYAYAMH0fYSp9PwLRTUrl2bN2zIFM5dIYdN9niF7Mf58+fBzKbnXm7cuIGrV69aWqyocvnyZSQmJga0liAYMuI9SEhIcNcVM2Pz5s2oVq2a5bUkmQki2sge3zSySrYyVWVlRGkIvrD6jYbIyMigJ67z5Mljm5OJFTLiPfBUsERk2vX2dib9ZmcEQRCE2wJRHIIgCIIlsvQcBxElAdib0XJkEooA+DejhcgkSF3cQuriFlIXt6jMzAEv6sjqcxx7g5nguZ0gog1SFw6kLm4hdXELqYtbEFFQXkViqhIEQRAsIYpDEARBsERWVxxjM1qATITUxS2kLm4hdXELqYtbBFUXWXpyXBAEQUh/svqIQxAEQUhnRHEIgiAIlsi0ioOIyhDRMiLaTUQ7iaiv83ghIoolov3OvwU90gwmogNEtJeIWmSc9Pbipy5GEdEeItpGRLOJqIBHmmxVFx7nBxIRE1ERj2PZri6I6HXn9e4kopEex7NVXRBRTSKKJ6ItRLSBiB7wSHO71kUUEa0joq3OuvjQedy+tjOYHRJD+QNQEkAt5/95AewDUA3ASACDnMcHAfjM+X81AFsB5ARQHsBBAOEZfR0hrovmACKcxz/LznXhDJcBsAjAUQBFsmtdAHgYwBIAOZ3nimXjulgMoJXzeGsAy7NBXRCAPM7/cwBYC+BBO9vOTDviYOZTzLzJ+X8SgN0A7gTQHsBEZ7SJAB53/t8ewFRmTmbmwwAOAAjdtxPTEV91wcyLmTnFGS0egOvzfdmuLpynvwTwNhzbjLvIjnXRG8AIZk52njvrTJId64IBuD4BmB/ASef/t3NdMDNfdgZzOH8MG9vOTKs4PCGiaAD3w6E5izPzKcDxsAAo5ox2J4B/PJIdx60G5bZBqQtPXgTwp/P/bFcXRNQOwAlm3qpEy3Z1AaASgIZEtJaI4ojI9Q3l7FgX/QCMIqJ/AHwOYLAz2m1dF0QUTkRbAJwFEMvMtradmV5xEFEeADMB9GPmRH9RdY7dVr7GvuqCiP4DIAWA6/uv2aou4Lj2/wB4Xy+qzrHbti6cz0UEgIJwmCfeAvAbOfYjz4510RtAf2YuA6A/gHGuqDrJb5u6YOZUZq4JhxXiASKq7ie65brI1IqDiHLA8RD8wsyznIfPEFFJ5/mScGhUwKElPb8NWRq3hqVZHh91ASLqBqANgM7sNFgi+9VFRThss1uJ6Agc17uJiEog+9UF4LjmWU6TxToAaXBs8Jcd66IbANf/03HLBHNb14ULZk4AsBxAS9jZdmb0RI7BBM8kAF8px0dBO8Ez0vn/PdBO8BzC7TXZpVcXLQHsAlBUOZ7t6kKJcwS3JsezXV0A6AXgI+f/leAwQ1A2rYvdAJo4/28KYGM2eC6KAijg/D8XgJVwdC5tazsz/CL9XHwDOIZL2wBscf5aAygMYCmA/c6/hTzS/AcOj4C9cHpS3A4/P3VxwNkouI6Nya51ocRxK47sWBcAIgH8DGAHgE0AHsnGddEAwEZnw7gWQEw2qIv7AGx21sUOAO87j9vWdsqWI4IgCIIlMvUchyAIgpD5EMUhCIIgWEIUhyAIgmAJURyCIAiCJURxCIIgCJYQxSEIgiBYQhSHIAiCYAlRHIJgESLK5dw8MDwEeUcS0QoiirA7b0GwC1EcgmCdF+HYCyrV7oyZ+QYcq3qfsTtvQbALURyC4AMi+sv55bgtRHSdiJ5ynuoM4HePeNOJ6DsiWkVER4moARFNIqJ9RDTOajwAc5xlCEKmRLYcEQQDiKg3HF/V6wQgHMAxZi7hcX4PgLHMPJqIPoJjtNAEwHkAZwCUYOZkC/HCAZxm5qLpdpGCYAEZcQiCH4ioK4BWcGxbnwrH9uQJHuejABQA8JXz0DUA49jxRbobAK4CuGE2HuD4loIzTd5QXpsgBIooDkHwgdM01RnA08x803n4GoAoj2j3ANjEzGnOcA04v85IRKUBnGTHsN5sPBc5AVy3/6oEIXhEcQiCDkTUBsCrAJ5gZncDzswXAYQ7RxAAcC8cW3a7uA+O7awBh3LYZjEeiKgwgHMeykoQMhWiOARBn4lwfAntb+fk+Ese5xbD8Z0HwKEQtgBus1Uup3IBtMrBbDzAMZ+ywM6LEQQ7kclxQbAIEd0P4E1mfj5E+c8CMJiZ94Yif0EIFhlxCIJFmHkzgGWhWgAIYI4oDSEzIyMOQRAEwRIy4hAEQRAsIYpDEARBsIQoDkEQBMESojgEQRAES4jiEARBECwhikMQBEGwxP8DpofH7n7OEeQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u_z.draw('field', z_scale='mm')\n" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.10.6 64-bit", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": false, "sideBar": true, "skip_h1_title": true, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false }, "vscode": { "interpreter": { "hash": "31f2aee4e71d21fbe5cf8b01ff0e069b9275f58929596ceb00d14d90e3e16cd6" } }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }