Skip to content

Instantly share code, notes, and snippets.

@krischer
Last active March 15, 2019 22:27
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save krischer/c3180c6fb2af4b430bf12c9ad5b95eda to your computer and use it in GitHub Desktop.
Save krischer/c3180c6fb2af4b430bf12c9ad5b95eda to your computer and use it in GitHub Desktop.
Synthetic relative calibration example
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEACAYAAAC6d6FnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xl8VPW9//HXJ5MNEgiQsIYl7IsCsghu1WrVogUXtO5WEaR4q7d7q/faen/tz+q1O3VFQETckNKKilJrq5SqQEBlR5E1AhIIa0L27/0jAUPMMiEzc87MvJ8P5jHnfOfMOZ8Pk8wn37N9zTmHiIjEnwSvAxAREW+oAIiIxCkVABGROKUCICISp1QARETilAqAiEicUgEQEYlTKgAiInFKBUBEJE6pAIiIxKnESG3IzNKAR4FS4G3n3LONvScrK8vl5OSEOzQRkZiyYsWKvc659o0t16wCYGYzgbHAHufcqTXaxwB/BALAdOfcg8B4YJ5z7hUzexFotADk5OSQm5vbnBBFROKOmW0LZrnm7gKaBYypteEA8AhwCTAIuN7MBgFdgR3Vi1U0c7siItJMzeoBOOcWm1lOreZRwCbn3GYAM3sBuBzIo6oIfIiOPcSN0vJK9hWWsPdwKQVFpRSVlFNYWkFRaTmFJVXPZRWOispKyisdFZXuS/PlldV3rD3+VDVx7Ea2NW9o+6XXqL2M7n4r0eH/XzGYThmpYd1GOI4BZPPFX/pQ9cU/GpgKPGxm3wBeqe/NZjYZmAzQvXv3MIQn4VBQWMrKbfvZ+PlhPs0/wqf5hWzfV8j+orIG32cGSYEEEhOMQIJVP1fNJwaq5hMSDDu+fNWU1Xh/1bydsM4Tt1H3e0T8rLyyMuzbCEcBqOvXyznnCoEJjb3ZOTfNzHYB45KTk0eEPDoJiYpKx/KtBSxau5t3NuazeW/h8dc6Z6TSu306lw7uTIdWqbRvlUJWejKZ6cmkpSSSlpxIy+QALZMTSU1KOP4FLSKRFY4CkAd0qzHfFdjZlBU4514BXhk5cuTtoQxMmm/bvkLmvL+N+Ss/Y19hKcmJCZzdO5NvjuzGiB5tGdSlNekpETu5TESaIRy/qcuBvmbWE/gMuA64oSkrMLNxwLg+ffqEITw5GRt2H+K3f/uYv6//nIAZFw3qyLihXTivX3vS9IUvEpWaexro88BXgSwzywPuc87NMLM7gUVUnQY60zm3tinrVQ/APz47cJSH3tjAgo92kp6SyF3n9+HGM3rQsXV4D06JSPg19yyg6+tpXwgsPNn1qgfgvYpKxzPvbeWhRRupdI4p5/Xm2+f2ok3LZK9DE5EQ8WXfXT0Ab23OP8IPX/qID7Yf4Nx+7bn/ilPp1q6l12GJSIj5sgCoB+Cd11bt4qd/XkViwPj9tUO54rRsnaUjEqN8eUGWc+4V59zkjIwMr0OJG6XllfzPgrV857mV9OuYzsL//ApXDuuqL3+RGKYegLDvSAm3z85l5fYDTDynJz8dM4DkRF/+bSAiIeTL33L1ACJnc/4Rxj/2Lmt3HuKRG4bzs7GD9OUvEid82QOQyFi+tYDbZ+cSMOP5yWcwvHtbr0MSkQjy5Z96ZjbOzKYdPHjQ61Bi1murdnHj9KW0a5nM/P84S1/+InHIlwVAu4DCa27uDu56fiVDsjP48x1n0SMzzeuQRMQD2gUUZ55+dyv3LVjLV/pmMe3mkbRIDngdkoh4RAUgjjz29qf87xsbuGhQRx6+YRgpifryF4lnvtwFpGMAoeWc47d/28j/vrGBcUO78OiNw/XlLyL+LAA6BhA6zjkeWrSRP/1jE9eM7Mofrj2NpIAvP3YRiTDtAopxv3/zYx57+1OuH9Wd+684lYQEXdkrIlX0p2AMm/rWJ0yt/stfX/4iUpsKQIx69O1N/O7Njxk/LJsHxg/Rl7+IfIkKQAx6cvFmHnpjI5cN7cKvvzmUgL78RaQOviwAOgvo5M1csoX7F67nG4M787tr9OUvIvXzZQHQWUAn55n3tvKLV9cx5pRO/OG600jU2T4i0gB9Q8SI55Zu52cvr+XCgR2Yev0wneopIo3St0QM+MsHefz3X1dzfv/2PHLjcN3OWUSCom+KKLdo7W5+9NIqzuyVyWM3jdAVviIStIgVADPrZWYzzGxepLYZ65Z8spe7nvuAwdkZPPmtkaQm6ctfRIIXVAEws5lmtsfM1tRqH2NmG81sk5nd3dA6nHObnXMTmxOsfGHFtv3cPjuXXu3TmDXhdNJSdFG3iDRNsN8as4CHgdnHGswsADwCXATkAcvNbAEQAB6o9f7bnHN7mh2tALBu5yEmPLWMjq1TmD1xFG1aJnsdkohEoaAKgHNusZnl1GoeBWxyzm0GMLMXgMudcw8AY0MZpHxhc/4RvjVzKekpicyZNJoOrVK9DklEolRzjgFkAztqzOdVt9XJzDLN7HFgmJnd08Byk80s18xy8/PzmxFe7PnswFFumr4U5+CZSaPp2ral1yGJSBRrzo7jui4xdfUt7JzbB0xpbKXOuWlmtgsYl5ycPKIZ8cWU/MMl3DR9KYdLynlh8hn0bp/udUgiEuWa0wPIA7rVmO8K7GxeOFV0JfCJDhSVcvOMpew+WMysCadzShf9v4hI8zWnACwH+ppZTzNLBq4DFoQiKN0L6AtHSsq59anlbM4v5MlvjWREj3ZehyQiMSLY00CfB94D+ptZnplNdM6VA3cCi4D1wFzn3NpQBKUeQJXisgomz85l9WcH+dMNwzinb5bXIYlIDAn2LKDr62lfCCwMaURU9QCAcX369An1qqNGWUUldz63knc/3cfvrhnK10/p5HVIIhJjfHkriHjvAVRWOn700kf8ff0efnn5KYwf3tXrkEQkBvmyAMTzMQDnHD97eQ0vf7iTH3+9PzefmeN1SCISo3xZAOK1B+Cc48E3NvDs0u1MOa833zk/fneBiUj4+bIAxGsP4NG3P+WJdzZz0xnd+emY/l6HIyIxzpcFIB57AE+/u5VfL9rIlcOy+cVlp2KmoRxFJLx8WQDizZ9X5HHfgrVcNKgjv756CAkax1dEIsCXBSCedgG9sWY3P573EWf3yeRP1w/TOL4iEjG+/LaJl11A//okn/98/gNO69aGaTdrQBcRiSxfFoB4kLu1gMmzV9C7QzpP3TpKA7qISMSpAHhgzWcHmTBrOZ0zUpl92ygyWiZ5HZKIxCFfFoBYPgawac8Rbpm5jFYpiTwzaTTtW6V4HZKIxClfFoBYPQaQt7+Im2csxQzmTBpNdpsWXockInHMlwUgFu05VMyN05dSWFLOMxNH00sDuoiIx3TkMQKqBnRZRv7hEuZMGs3Azq29DklERD2AcDtSUs4tTy1ny96qAV2Gd2/rdUgiIoB6AGFVXFbBpKeXs+azgzx243DO7qMBXUTEP3zZA4iFs4DKKir5zrMrWbqlgN9+cygXa0AXEfEZXxaAaD8LqKLS8YO5H/HWhj388vJTuWJYttchiYh8iS8LQDRzznHvX1fzykc7ufuSAdx0Rg+vQxIRqZMKQAg553jg9Q08v2wH3zm/N1PO6+11SCIi9VIBCKGH/7GJaYs3c8uZPfjRxRrQRUT8LWIFwMyuMLMnzexlM7s4UtuNlKf+vYXfvvkx44dnc9+4UzSgi4j4XlAFwMxmmtkeM1tTq32MmW00s01mdndD63DO/dU5dztwK3DtSUfsQy/l7uD/vbKOr5/SkYeu0oAuIhIdgr0OYBbwMDD7WIOZBYBHgIuAPGC5mS0AAsADtd5/m3NuT/X0vdXviwmvr97FT/+8iq/0zWKqBnQRkSgSVAFwzi02s5xazaOATc65zQBm9gJwuXPuAWBs7XVY1T6RB4HXnXMrmxO0X7zzcT7/+cIHDOveliduHkFKogZ0EZHo0Zw/V7OBHTXm86rb6nMXcCFwtZlNqW8hM5tsZrlmlpufn9+M8MJr+dYCvv1MLn07tGLmrafTMlkXVYtIdGnOt1ZdO7pdfQs756YCUxtbqXNumpntAsYlJyePaEZ8YbPms4Pc9tRyurRpweyJo8hooQFdRCT6NKcHkAd0qzHfFdjZvHCq+PlK4E17DvOtmcto3SKJORNHk5WuAV1EJDo1pwAsB/qaWU8zSwauAxaEIii/3gtoR0ERN01fRoIZcyaNposGdBGRKBbsaaDPA+8B/c0sz8wmOufKgTuBRcB6YK5zbm0ogvJjD+Dz6gFdjpZVMGfSKHpmpXkdkohIswR7FtD19bQvBBaGNCKqegDAuD59+oR61Sdlf2EpN89Yyt4jJTw7aTQDOmlAFxGJfr48ad1PPYDDxWXc8tQytu4rYvotIxmmAV1EJEb4sgD45RhAcVkFE5/OZd3OQzx6w3DO6q0BXUQkdviyAPihB1BaXskdc1awfGsBv71mKBcO6uhZLCIi4eDLAuB1D6Ci0vH9uR/yz4353H/FYC4/TQO6iEjs8WUB8LIH4Jzjv/+ymtdW7eK/Lh3ADaO7RzwGEZFI8GUB8IpzjvtfW88Ly3dw1wV9mHyuBnQRkdjlywLg1S6gqW9tYvqSLdx6Vg4/uKhfRLctIhJpviwAXuwCmrlkC7//+8dcPaIrPx87SAO6iEjM82UBiLS5y3fwi1fXccmpnXhw/GAN6CIicSHuC8Brq3Zx9/yqAV3+cN1pGtBFROKGL7/tInUM4O2Ne/jeix8wXAO6iEgc8mUBiMQxgGVbCpgyZwX9OrZi5gQN6CIi8ceXBSDcVucd5LZZy8lu04LZt42idaoGdBGR+BN3BeCTzw/zrZlLyWiRxJxJo8nUgC4iEqfiqgBs31fEjdOXkhhI4NlJo+mcoQFdRCR+xU0B2H2wmBtnvE9pRSVzJo4mRwO6iEic82UBCPVZQAWFpdw0YykFR0p5esIo+ndqFZL1iohEM18WgFCeBXSouIxbZi5jR0ER0285naHd2oQgQhGR6OfLAhAqR0srmDQrl/W7DvHYTcM5s3em1yGJiPhGzJ78XlpeyZQ5K1i+rYCp1w3jggEa0EVEpKaY7AFUVDq+/+KHvPNxPg9cOZhxQ7t4HZKIiO9ErAdgZgOB7wJZwFvOucfCsZ3N+Ue4/7X1vLVhD/d+YyDXjdKALiIidQmqAJjZTGAssMc5d2qN9jHAH4EAMN0592B963DOrQemmFkC8GSzoq7D7oPF/PGtT5ibu4OUxAR+NnYQE8/pGerNiIjEjGB7ALOAh4HZxxrMLAA8AlwE5AHLzWwBVcXggVrvv805t8fMLgPurl5XSBwoKuWxdz5l1r+3UukcN5/Rgzsv6EOWrvAVEWlQUAXAObfYzHJqNY8CNjnnNgOY2QvA5c65B6jqLdS1ngXAAjN7DXjuZIOGqjN8nnp3C4+//SmHS8q58rRsvn9RP7q1a9mc1YqIxI3mHAPIBnbUmM8DRte3sJl9FRgPpAALG1huMjAZoHv3L++/L6uo5MXlO5j61ifsOVzC1wZ04Mdj+jOgU+uTSkJEJF41pwDUNWyWq29h59zbwNuNrdQ5N83MdgHjtm/f3tfMPqm1SAZw/BLhmdWPepywbBCvBdNWc76+6Sxgb/1hBaWh2JuybH2vNZRXQ/OhzrOhGJuyXKjzrDkdjXnW1ea3n93GlouV31Ev8uwRRFzgnAvqAeQAa2rMnwksqjF/D3BPsOs72QcwLRTL1vVaMG015xuYzvVzno3l1dB8qPNsSq6RzLNWzlGXZ7C5heMzDUWeweakPBuebuzRnOsAlgN9zaynmSUD1wELmrG+YL0SomXrei2YtleCmA6FcOZZV3uw86HOsynrjGSejW3vZEQyz7ra/Paz29hysfI76ts8rbpiNLyQ2fPAV6nqNn0O3Oecm2FmlwJ/oOrMn5nOufuD3XAsM7Nc59xIr+MIN+UZe+Il13jJszHBngV0fT3tC2nggG4cm+Z1ABGiPGNPvOQaL3k2KKgegIiIxJ6YvBeQiIg0TgVARCROqQCIiMQpFYAIM7OBZva4mc0zszu8jidczOwKM3vSzF42s4u9jidczKyXmc0ws3lexxJqZpZmZk9Xf443eh1PuMTyZ9gYFYAmMLOZZrbHzNbUah9jZhvNbJOZ3d3QOpxz651zU4BrAF+ehhaiPP/qnLsduBW4NozhnrQQ5bnZOTcxvJGGThNzHg/Mq/4cL4t4sM3QlDyj7TMMJRWAppkFjKnZUOOuqJcAg4DrzWyQmQ02s1drPTpUv+cyYAnwVmTDD9osQpBntXur3+dHswhdntFiFkHmDHTli/t9VUQwxlCYRfB5xq2YHRIyHJwP74oaDqHI08wMeBB43Tm3MrwRn5xQfZ7RpCk5U3WDx67Ah0TZH4tNzHNdZKPzj6j6UH2qrruiZte3sJl91cymmtkTRNdFdE3KE7gLuBC42symhDOwEGvq55lpZo8Dw8zsnnAHFyb15TwfuMrMHiM8twCJtDrzjJHP8KSoB9B8Ybkrqg81Nc+pwNTwhRM2Tc1zHxBNBa4udebsnCsEJkQ6mDCqL89Y+AxPinoAzZcHdKsx3xXY6VEs4aQ8Y1e85BwveQZNBaD5vLoraqQpz9gVLznHS55BUwFoAqu6K+p7QH8zyzOzic65cuBOYBGwHpjrnFvrZZzNpTxjK8+a4iXneMmzuXQzOBGROKUegIhInPL1WUBZWVkuJyfH6zBERKLKihUr9jrn2je2nK8LQE5ODrm5uV6HISISVcxsWzDLaReQiEicilgPwMxSgcVASvV25znn7mvoPRWVjv2FpSdchXPsoPWJbdXPNVpTEgO0Tk2k6o4EIiJSWyR3AZUAFzjnjphZErDEzF53zr1f3xvW7TrEsF++edIbTAoY7dKSaZeWQmZaMpnpybRLSyazuq1dWjJZx9tSaN1CBUNE4kfECoCr+tP9SPVsUvWjwXNQO2ek8j/jqm7WV/OL+djkCV/V1Y3H2o6WVrCvsJSCwhIKCkvZV1jK9u1FFBSWcqSkvM7tJSYcKxhVxSIzLYXu7VrSI7MlPbPSyMlKIzMtWUVCRGJCRA8CV9+OdQXQB3jEObe0oeWz0lO49eyeIY+juKyC/UWl7DtSerxIHJ+u0fZBwX5eW72Lisov6lR6SiI5WS3pkZlGz8y048WhX6dWtE5NCnmsIiLhEtEC4JyrAE4zszbAX8zsVOdc7QEbJgOTAbp37x6WOFKTAnTOaEHnjBaNLltWUUne/qNs3VvI1n2F1c9FrPnsIG+s2X1CceiZlcbg7AyGdM1gcHYGp2RnkJ7i6xOtRCSOeXYlsJndBxQ6535T3zIjR450fj4N9Fhx2LL3COt2HmL1ZwdZnXeQnQeLgaq9Ur3bpzOsWxvO7pPFWX0y6dAq1eOoRSTWmdkK51yjIw5G8iyg9kCZc+6AmbWg6l7x/xup7YdDUiCBnllp9MxK44IBHY+35x8uYc1nB1mVd5BVeQf427rPeWlFHgD9OqZzVu8szqkuCC2T1UMQEW9ErAdgZkOAp4EAVdcfzHXO/aKh9/i9BxCsikrHup2H+Pene/n3pr0s21JASXklqUkJnN+/A2NO7cTXBnbU7iIRCYlgewC+vhlcrBSA2orLKsjdup9Fa3fzxtrd5B8uITkxgfP6teeakd04v397EgO6Rk9ETo4KQJSorHSs3L6fhat388qqneQfLqFDqxSuHtGVa0/vRo/MNK9DFJEoowIQhcorKvnnxnxeXL6df2zYgwO+PqgTk8/rxfDubb0OT0SihO8OAkvjEgMJXDSoIxcN6sjug8XMeX8bs9/byhtrd3N6Tlu+d2E/zu6T5XWYIhIj1APwucKScl5cvoMn/7WZXQeLObtPJj/5+gCGdmvjdWgi4lPaBRRjissqeHbpdh755yYKCku5clg291wygA6tdV2BiJwo2AKgU02iRGpSgInn9GTxT87nzvP78NqqXVzw23eYuWTLCVcji4gESwUgyqSnJPKjr/dn0ffPZUSPtvzi1XVc+8R7bNlb6HVoIhJlVACiVM+sNGZNOJ3fXTOUjz8/zCV/XMzzy7bj5116IuIvKgBRzMwYP7wrf/v+eYzs0Y575q/m+y9+WO/trkVEalIBiAGdMlJ5+rZR/PCifiz4aCdXPfouOwqKvA5LRHxOBSBGBBKMu77Wl9m3jWbXwaNc9vASlm8t8DosEfExFYAYc07fLF6+8xzatkzmpulLeXvjHq9DEhGfUgGIQT2z0nhpypn06ZDO7bNzWbh6l9chiYgPqQDEqMz0FJ67/QyGdm3Dnc+t5KXcHV6HJCI+owIQwzJaJDF74ijO7pPFj+etYv7KPK9DEhEfUQGIcS2TE5l+y8iqewjNW8U/dUxARKqpAMSBlMQAj980gv6dWvEfc1aycvt+r0MSER9QAYgTrVKTmDVhFB1ap3DbrOV8mn/E65BExGMqAHGkfasUZt82ioAZk2fncri4zOuQRMRDKgBxpkdmGg/fMJyt+4r4wdyPqNSdREXilgpAHDqzdyb/felA3lz3OX/6xyavwxERj0SsAJhZNzP7p5mtN7O1ZvbdSG1bvmzC2TmMH57N7//+Me98nO91OCLigUj2AMqBHzrnBgJnAN8xs0ER3L7UYGb86srB9OuYzg/nfsTeIyVehyQiERaxAuCc2+WcW1k9fRhYD2RHavvyZalJAaZeP4xDxWX86CUdDxCJN54cAzCzHGAYsLSO1yabWa6Z5ebna9dEuA3o1JqffWMgb2/M56l3t3odjohEUMQLgJmlA38GvuecO1T7defcNOfcSOfcyPbt20c6vLh00xk9uHBgBx56YwObdX2ASNyIaAEwsySqvvyfdc7Nj+S2pX7HjgekJgX4ybxVGmReJE5E8iwgA2YA651zv4vUdiU4HVqnct+4QeRu288s7QoSiQuR7AGcDdwMXGBmH1Y/Lo3g9qURVw7L5oIBHfj1og1s3VvodTgiEmaRPAtoiXPOnHNDnHOnVT8WRmr70rhju4KSAgn8ZN4qnRUkEuN0JbCcoFNGKj8bO4hlWwt4UYPIiMQ0FQD5km+O6Mronu148PUNukBMJIapAMiXmBn3X3kqRaXl/Oq19V6HIyJhogIgderToRXfPrc38z/4jHc/3et1OCISBioAUq87L+hD93Ytufcvaygpr/A6HBEJMRUAqVdqUoBfXnEqm/cW8sQ7m70OR0RCTAVAGnRev/aMHdKZh/+5iS26NkAkpqgASKN+PnYQKYEEfv7yGpzTtQEisUIFQBrVoXUqPx7Tn399spcFH+30OhwRCREVAAnKjaN7MLRrBr98dR0Hikq9DkdEQkAFQIISSDAeGD+E/UVl/Gqhrg0QiQUqABK0QV1ac/tXejE3N0/XBojEABUAaZLvfq0v3du15L/mr6a4TNcGiEQzFQBpkhbJAX515WC27iti6lufeB2OiDSDCoA02Tl9s7hqeFemLd7M+l1fGtVTRKKECoCclHu/MZDWLZK4e/5qDSEpEqVUAOSktE1L5udjB/HRjgPMWKLbRIhEIxUAOWmXn9aFiwd15DeLPtauIJEopAIgJ83MeGD8YFq3SOJ7L3yos4JEoowKgDRLZnoKv/7mEDZ+fpjfLNrodTgi0gQRLQBmNtPM9pjZmkhuV8Lr/P4duPmMHkxfsoV/b9IFYiLRItI9gFnAmAhvUyLgvy4dSO/2aXz3hQ/5/FCx1+GISBAiWgCcc4uBgkhuUyKjRXKAx28aQVFpOf/x7EpKyyu9DklEGuG7YwBmNtnMcs0sNz8/3+twpAn6dmzFQ1cPYcW2/bphnEgU8F0BcM5Nc86NdM6NbN++vdfhSBONHdKFSef0ZNa7W/nzijyvwxGRBviuAEj0++klAzizVyZ3z1/FuzooLOJbKgASckmBBB6/eQQ5mWl8e84KNu4+7HVIIlKHSJ8G+jzwHtDfzPLMbGIkty+Rk9EiiVm3jaJFUoAJTy0jb3+R1yGJSC2RPgvoeudcZ+dcknOuq3NuRiS3L5GV3aYFT004ncMl5Vz/5PvsPHDU65BEpAbtApKwOqVLBnMmjuZAYRnXTXufXQdVBET8QgVAwm5otzbMnjiK/YWlXDftfbbtK/Q6JBFBBUAiZFj3tsyeOIpDR8u46rF3WZV3wOuQROKeCoBEzLDubZl3x1mkJgW4btr7/H3d516HJBLXVAAkonq3T2f+HWfRu306k2bn8vs3P6ZSI4qJeEIFQCKuQ+tUXppyJlcN78of3/qESbNzKSgs9ToskbijAiCeSE0K8JtvDuGXl5/Ckk/2cvHvF/OPDdolJBJJKgDiGTPj5jNzePnOs8lKT+a2WbncM38Vh4rLvA5NJC6oAIjnBnZuzct3ns23z+3FC8t3cMFv3mH+yjyc07EBkXBSARBfSEkMcM+lA3n5O2eT3bYFP5j7Edc+8T4fbN/vdWgiMUsFQHxlSNc2/OWOs3hw/GA+zT/ClY++y+2zc9mw+5DXoYnEHPNzN3vkyJEuNzfX6zDEI4Ul5Tz17y08sXgzR0rKuWhgRyaf24sRPdpiZl6HJ+JbZrbCOTey0eVUAMTvDhSVMv1fW3jm/W0cPFrGad3aMPGcnlx8SkdSEgNehyfiOyoAEnOKSsuZtyKPGUu2sG1fEW1aJnHFadlce3o3BnZu7XV4Ir6hAiAxq6LSsWTTXubm7uDNtZ9TWlHJwM6tueTUTlxyaif6dmzldYginlIBkLiwv7CUlz/8jFdX7WLF9v04B73bp3HxKZ34St8sRvRoq91EEndUACTu7DlUzKK1u3l9zW6WbSmgvNLRIinA6F7tOKdPFqN6tmNg59YkBXTym8Q2FQCJa4eLy3h/cwFLPsnnX5v2sjm/agyC1KQEhmS3YViPNgzv3pZTurQmu00LnVUkMUUFQKSGnQeOsnL7flZuO8DK7ftZu/MgZRVVP/utUhMZ0KkVAzq1ZkDnVvTr2IqczDSy0pNVGCQqBVsAEiMRjIjXurRpQZc2LRg7pAsAxWUVrN15iA27D7F+1yE27DrMXz74jCPvlx9/T3pKIj0yW5KTmUZOVku6t2tJ54wWdM5IpVNGKq1Sk7xKRyQkIloAzGwM8EcgAEx3zj0Yye2LHJOaFGBEj7aM6NH2eJtzjrz9R9mUf4StewvZtq+ILXsLWbvzIG+s3U1FrXEL0lMS6ZSRSueMVDrJztJ4AAAFxklEQVS0SiUzPZl2aVWPzOPPKbRLTyYtOaDehPhOxAqAmQWAR4CLgDxguZktcM6ti1QMIg0xM7q1a0m3di2h/4mvlVVUsvtgMbsPFbPrYDG7Dx6tfq6a/3TPXvYVllJSXlnnupMDCaSnJtKq+pGekkir1CRapRxrSyI9NZG05ACpSQFaJAdITax+TkqoaksKnPCckphAQoKKipy8SPYARgGbnHObAczsBeByQAVAfC8pkPBFcaiHc46i0goKCkvZV1hKQWEJ+46UUlBYyv6iMg4Xl3GkpJzDxeUcLi5jR0HR8ekjJeWczMBoyYEEEgNGUiCh+mG1nqumEwMJJNeaDiQYgQTDDAJWNZ2QYATMSDCOT9fXnlD9/oBVryPBSKieNgAzrOoJ44v2Y/Mcn6+xXI1lqf1arXVwwvyJ66h+tTqQE9VVMuvrndW9bF3L1fP+ILdfV+Owbm1pkRzeU5gjWQCygR015vOA0RHcvkhYmRlpKYmkpSQ2WCjqcqx4FJVWUFxWwdGy6ufSCorLKzlaWkFJedV81WuVHC2roKyikrLySsorHaX1TJdVVFJaXrX8oeKq6bKKSioqHZWO6mdX6xkqKx0VNdqPLSuR8Y8fnkev9ulh3UYkC0Bdhe9LP01mNhmYDNC9e/dwxyTiCzWLh9/VVxgqq+cd4Bw4HNX/js+74/NV09T3WnU7J7TXWK7GdO11VL2LOseTqKt81X8iZB3vr2PZ+t5e57JBxgTQOaNFfYGFTCR/2vKAbjXmuwI7ay/knJsGTIOq00AjE5qIBCshwUjASNIF1lEvkpdELgf6mllPM0sGrgMWRHD7IiJSQ8R6AM65cjO7E1hE1WmgM51zayO1fREROZGvrwQ2s3xgW63mDOBgkKtoaNm6XgumreZ8fdNZwN4gY6xPOPOsqz3Y+VDn2VCMTVku1HnWnI7GPOtq89vPbmPLxcrvqBd5tnHOtW80sqqDMdHzAKaFYtm6XgumreZ8A9O5fs6zsbwamg91nk3JNZJ51so56vIMNrdwfKahyDPYnJRnw9ONPaLxtoivhGjZul4Lpu2VIKZDIZx51tUe7Hyo82zKOiOZZ2PbOxmRzLOuNr/97Da2XKz8jvo2T1/vAopWZpbrgrgRU7RTnrEnXnKNlzwbE409gGgwzesAIkR5xp54yTVe8myQegAiInFKPQARkTilAiAiEqdUAERE4pQKQISZ2UAze9zM5pnZHV7HEy5mdoWZPWlmL5vZxV7HEy5m1svMZpjZPK9jCTUzSzOzp6s/xxu9jidcYvkzbIwKQBOY2Uwz22Nma2q1jzGzjWa2yczubmgdzrn1zrkpwDWAL09DC1Gef3XO3Q7cClwbxnBPWojy3OycmxjeSEOniTmPB+ZVf46XRTzYZmhKntH2GYaSCkDTzALG1GyoMdLZJcAg4HozG2Rmg83s1VqPDtXvuQxYArwV2fCDNosQ5Fnt3ur3+dEsQpdntJhFkDlTdcfeY2N4VEQwxlCYRfB5xi3/33zcR5xzi80sp1ZznSOdOeceAMbWs54FwAIzew14LnwRn5xQ5GlVQyw9CLzunFsZ3ohPTqg+z2jSlJypuoV7V+BDouyPxSbmGbejEkbVh+pTdY10ll3fwmb2VTObamZPAAvDHVwINSlP4C7gQuBqM5sSzsBCrKmfZ6aZPQ4MM7N7wh1cmNSX83zgKjN7jPDcAiTS6swzRj7Dk6IeQPMFNdLZ8Recext4O1zBhFFT85wKTA1fOGHT1Dz3AdFU4OpSZ87OuUJgQqSDCaP68oyFz/CkqAfQfEGNdBYDlGfsipec4yXPoKkANF+8jHSmPGNXvOQcL3kGTQWgCczseeA9oL+Z5ZnZROdcOXBspLP1wFwX5SOdKc/YyrOmeMk5XvJsLt0MTkQkTqkHICISp1QARETilAqAiEicUgEQEYlTKgAiInFKBUBEJE6pAIiIxCkVABGROKUCICISp/4PBZrPr8xUDN8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEACAYAAAC6d6FnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xd4XPWZ9vHvMzNqlqxiSe6W5W6MCy4YcAgxoZmlZSkhhHTeEMhLNsmmkCx5k2yyG5JNstmQuoawJpuYBEwznVAMJGDcsI27cZEsy7YsWbLVpZn5vX/MYIQjWSNrqub+XJeuOec3Z848j8fSPefMOWfMOYeIiKQfT6ILEBGRxFAAiIikKQWAiEiaUgCIiKQpBYCISJpSAIiIpCkFgIhImlIAiIikKQWAiEiaUgCIiKQpX6ILOJmSkhJXXl6e6DJERFLK2rVra51zpb0tl9QBUF5ezpo1axJdhohISjGzikiW0y4gEZE0ldRbACID2dHWTirqmqluaGV/Qxu1Te00t/tpavPT0hHA4fCYYQZmxqAML7lZPnKzwreZPnKzfAzODv3kZ2ccv83L9pHh1fs7OTkFgEgc+ANBNlQdZeXuOjZWNbC5+hhV9a3vWSbDa+Rlhf6oD8r0YhgOR9BBMOho7QzQ3O6nuSNAINj7ZdxzMrzHw2Fwl3A4cWxwl7H8E8YUIgObAkAkRuqa2nluyyFe2HqIlbuP0NTuB2BcSS6zxhTy0bPKmFiax6iiHEYV5lCQk4GZ9bpe5xzt/mBoa6HdT2Obn2NtnTS2+cM/nSfcvnt/dUPr8bHWzkCvz5Wd4XlPUORlecnyecnyecj0ecjyeY7PZ2V0mfZ5yMp4dzmfx4PPY/i8Fpr2Gj6P4fUYGV5P+NbwdlnO6zEyPB684WU99s4P79kyklOnABCJorqmdp7edJCn3jrAyt11BB2MGZLDFbNGcu7EEs6ZUMyQ3Mx+PYeZkZ3hJTvDS3Fe1imvpzMQpKlLQEQSIi0dAY61+mn3B2j3B2nvDB6fbusMEMGGSVRZOAw84TA4Hg68GxIej3W7zPH7u4wRzhM7vn7rMg3vzHXNnXdCyLqMH7/Fuky/e8d719n9c/38I7MZWZgTjX+mHikARPrJHwjy6s5a/rx6Hy9sO0RnwDG+JJfPL5zIpTOGM21EflK+U83weijKzaSon4HUlT8QpN0fpMMfun0nHPwBhz8YxB90x6cDx6cd/kDovkDQ0RkI3dcZdATC4/6gI+gcLrw7LOgIz787HXSE59+Z7n2Z4+t07vhuteMZ5sCF51x4faHhLmPvLOrefeS7y4We693pd8fpct+763rv4+NBASByivbWNvPAmn08tK6KQ8faKc7N5JPnlHPN3NFMHT44Kf/ox5rP68Hn9ZB76hsmEkcKAJE+aPcHeGbTQZa+Uckbe47gMVg4ZSj/euUYPjh1KJk+fWgqqUMBIBKB3YebuH9VJcvWVlHf0knZkEF87ZIpXDt3NMPysxNdnsgpUQCI9KDDH+S5LQf548pKXt9dh89jXDRtGDeeNZYFE4rxeNJvF48MLAoAkRNU1DVz/6p9LFu7j9qmDkYX5fC1S6Zw3bzRDB2sd/sycCgARAgdEvn8lkMsXVXJqztr8XqMC6YO5aNnlXHepFK925cBSQEgaW3fkRb+tLqSB9ZUcbixnZEF2Xz5wslcf+YYhhfo3b4MbAoASTv+QJAXt9WwdFUlL+84jAHnTwm92184ZShevduXNKEAkLRR3dDKn1bv44HV+zh4rI1h+Vl84YOTuP7MMYyK8RmXIslIASADWiDoeHlHDX9cWclL22twwHmTSvneVafzwalD8eliZ5LGFAAyIB061safV+/jz6v3sb+hldLBWdy6cAIfObOMMUMGJbo8kaSgAJABIxh0vLLzMEvfqOSFbTUEgo73TyrhW5edxoXThunSxiInUABIyqtpbOPBNVXcv6qSqvpWinMz+ez7x3PD/DGMLc5NdHkiSUsBICkpGHS8tquOpasqeG7zIfxBxznji7l90VQuOX24rskjEgEFgKSUuqZ2HlwberdfUddC0aAMPv2+cm6YX8b40rxElyeSUuIWAGY2HrgDKHDOXRuv55XU55xj5e4jLF1VyTObDtAZcMwfN4R/vmgyl5w+nOwMb6JLFElJEQWAmd0LXA7UOOemdxlfBPwc8AL3OOd+2NM6nHO7gZvMbFn/SpZ00dDSwbK1VSxdVcnuw83kZ/v42Nlj+ej8MiYNG5zo8kRSXqRbAEuAXwK/f2fAzLzAr4CLgCpgtZktJxQGd57w+M8452r6Xa0MeM451u9r4A8rK3liYzXt/iBzygr56XWzuGzmCL3bF4miiALAOfeKmZWfMDwfeDv8zh4z+xNwlXPuTkJbCyIRa273s3xDNX9YWcHm6mPkZnq5du5objxrLNNG5ie6PJEBqT+fAYwC9nWZrwLO6mlhMysG/h2YbWbfDAdFd8vdDNwMUFZW1o/yJBVsP9jIH9+o4JF1+2ls9zN1+GD+7UPT+dDsUeRl6RgFkVjqz29Yd1fM6vHrjJ1zdcAtva3UObcYWAwwb968OH49ssSLPxDkuS2HWPK3vazae4RMr4fLZo7gY2eXMaesKC2/S1ckEfoTAFXAmC7zo4Hq/pUjA9nR1k7+vLqS+16rYH9DK2OG5PDNS6dy3bwxDMnNTHR5ImmnPwGwGphkZuOA/cBHgI9GpSoZUHYfbmLJa3tZtraKlo4AZ48fwneumMYFpw3TpZdFEijSw0DvBxYCJWZWBXzHOfc7M7sNeJbQkT/3Ouc2x6xSSSnOhc7U/d1f9/DithoyvR6uPGMkn35fOaePLEh0eSJC5EcB3dDD+FPAU1GtSFJaMOh4bstBfr1iFxurjlKSl8mXLpzEjWeNpXRwVqLLE5EudJiFREWHP8ij6/fz25d3sftwM+XFg7jz6hlcPWcUWT4duy+SjBQA0i+tHQHuX1XJPa/upvpoG6eNyOcXN8zmH2aM0P59kSSnAJBT0tYZYOkblfx6xS5qm9qZXz6Ef796Bgsnl+owTpEUoQCQPmn3B3hg9T5++dLbHDrWzjnji/n1jXOYP25IoksTkT5SAEhEOgNBHlpbxS9efJv9Da3MG1vEz64/gwUTShJdmoicIgWAnJRzjmc3H+JHz2xjT20zs8YUcufVM3j/pBLt6hFJcQoA6dG6ynp+8ORW1lTUM3FoHvd8Yh4XnDZUf/hFBggFgPydirpm/uOZ7Tz51gFK8rK48+oZXDd3ND59qbrIgKIAkOPqmzu468Wd/GFlBT6Phy9eMImbzxtPrq7KKTIg6TdbCAQd96+q5CfPbedYayfXnzmGL184maH52YkuTURiSAGQ5tZWHOHbj21mc/Uxzh4/hO9eeTpTh+sLWETSgQIgTdU0tvGjp7fz0Loqhudn84sbZnP5zBH6gFckjSgA0kxnIMh9r+3l58/vpM0f4NaFE7jt/Inazy+ShvRbn0ZW7q7j/z26iZ01TXxgcinfuWIa40vzEl2WiCSIAiAN1Dd38IOntvLg2ipGF+Ww+ONzuWjaMO3uEUlzCoABzDnHo+v38/0ntnK0tZNbPjCBL14wiZxMXZ5ZRBQAA1ZFXTPfenQTr+6s5Yzw5RtOG6Gje0TkXQqAAaYzEOTuV3fz8+d3kuH18L2rTufGs8bq2vwi8ncUAAPIusp6/uXht9h2sJFLTh/Gv145neEFOplLRLqnABgAGts6+Y9ntvOHNyoYnp/N4o/P5eLThye6LBFJcgqAFPfyjsN886GNHDjWxifPKeerl0whT8f0i0gE9JciRR1t6eTfntzCg2urmFCay0O3LmBOWVGiyxKRFKIASEHPbznEvzzyFnXNHXx+4QT+6YJJZGfo0E4R6RsFQAqpb+7gu49v5rH11UwdPpjfffJMZowuSHRZIpKiFAAp4qm3DvDtxzbR0NLJly6cxOcXTiTTpy9oEZFTpwBIcocb2/nO8k089dZBpo/K539vOksndIlIVCgAkpRzjuUbqvnu8s00twf42iVT+Nx54/W1jCISNQqAJHToWBt3PPIWz2+tYXZZIT++diYThw5OdFkiMsAoAJKIc44H11bx/Se20OEP8q3LTuPT7xunyziISEzELQDM7DTgi0AJ8IJz7jfxeu5UUN3QyjcefotXdhxmfvkQfnTtTMaV5Ca6LBEZwCIKADO7F7gcqHHOTe8yvgj4OeAF7nHO/bCndTjntgK3mJkHuLtfVQ8gzjkeXBN61x9wjn+98nQ+fvZYPHrXLyIxFukWwBLgl8Dv3xkwMy/wK+AioApYbWbLCYXBnSc8/jPOuRozuxL4Rnhdaa/mWBvfePgtXtxWw/xxQ/jJtbMoKx6U6LJEJE1EFADOuVfMrPyE4fnA28653QBm9ifgKufcnYS2Frpbz3JguZk9CSw91aJT3TtH+Hz7sc20dQb49uXT+NSCcr3rF5G46s9nAKOAfV3mq4CzelrYzBYCVwNZwFMnWe5m4GaAsrKyfpSXnGqb2vnWI5t4ZvNBZpcV8tPrZul7eUUkIfoTAN29XXU9LeycWwGs6G2lzrnFwGKAefPm9bi+VPT0Wwe449FNNLX5+calU/ns+8frCB8RSZj+BEAVMKbL/Gigun/lDEz1zR18Z/lmlm+oZsaoAn764VlMHqbj+kUksfoTAKuBSWY2DtgPfAT4aFSqGkBe2HqIbzz8FvXNHXz5wsl8/vwJZOhsXhFJApEeBno/sBAoMbMq4DvOud+Z2W3As4SO/LnXObc5ZpWmmGNtnXzv8S0sW1vF1OGD+Z9Pncn0Ubpyp4gkj0iPArqhh/GnOMkHuunqlR2Huf2hjdQ0tnPb+RP5pwsm6cqdIpJ0dCmIKGpq9/ODp7ay9I1KJg7N4+GPzWXWmMJElyUi0i0FQJS8vquOry3bwP6GVm4+bzz/fNFkfUuXiCQ1BUA/tXYE+NEz21jy2l7Kiwfx4OfOYV75kESXJSLSKwVAP6ytOMJXH9zIntpmPrWgnK8vmsKgTP2Tikhq0F+rU9DWGeBnf9nB3a/uZkRBDks/exYLJpQkuiwRkT5RAPTRxqoGvvLABnbWNHHD/DLuuOw08rL0zygiqUd/uSLU4Q/yixd38usVuyjNy+K+z8znA5NLE12WiMgpUwBEYEv1Mb7y4Aa2HjjGNXNG8+0rplGQk5HoskRE+kUBcBL+QJDfrNjFXS/upCAnk7s/MY+Lpg1LdFkiIlGhAOjBzkONfOXBDWysOsoVs0byvStPpyg3M9FliYhEjQLgBIGg455Xd/PTv+wgL8vHr2+cwz/MGJHoskREok4B0MWe2ma++uAG1lbUc8npw/j3f5xBSV5WossSEYkJBQAQDDrue30vP3pmG5leD/91/RlcdcZIzPRlLSIycKV9AOw70sLXlm1g5e4jnD+llB9eM5Nh+dmJLktEJObSNgCccyxdVckPntyKmfEf18zkunmj9a5fRNJGWgZAdUMrtz+0kVd31nLuxBJ+dO1MRhXmJLosEZG4SqsAcM6xbG0V33t8CwHn+P6HpvOxs8r0rl9E0lLaBEDNsTb+5ZG3eH5rDfPHDeEn186irHhQossSEUmYAR8AzjmWb6jmO8s309oR4NuXT+NTC8rxePSuX0TS24AOgLqmdr716Cae3nSQ2WWF/OS6WUwozUt0WSIiSWHABsAzmw5wxyObaGzzc/uiqdx83ni8etcvInLcgAqAY22dPLf5EI+t38+rO2uZPiqfpdedwZThgxNdmohI0kn5AGhq9/P8lkM8sbGaV3bU0hEIMqowh68vmsJn3z+eDK8n0SWKiCSllAyAlg4/L2yt4cmNB3hpew3t/iDD87P5+DljuXzmCM4YU6hDO0VEepEyAdDWGWDF9hoe33iAF7fW0NoZoHRwFjfML+PymSOYU1akI3tERPogqQPAOfhLePfO81sO0dwRoDg3k2vmjuKyGSOZP26IPtgVETlF5pxLdA09MrPDQMUJwwXA0QhXcbJlu7svkrGu8z1NlwC1EdbYk1j22d14pPPR7vNkNfZluWj32XU6FfvsbizZ/u/2ttxA+R1NRJ+Fzrnev7TcOZdSP8DiaCzb3X2RjHWdP8n0mmTus7e+TjYf7T770ms8+zyh55TrM9LeYvGaRqPPSHtSnyef7u0nFQ+ReTxKy3Z3XyRjj0cwHQ2x7LO78Ujno91nX9YZzz57e75TEc8+uxtLtv+7vS03UH5Hk7bPpN4FlKrMbI1zbl6i64g19TnwpEuv6dJnb1JxCyAVLE50AXGiPgeedOk1Xfo8KW0BiIikKW0BiIikKQWAiEiaUgCIiKQpBUCcmdlpZvZbM1tmZrcmup5YMbMPmdndZvaYmV2c6HpixczGm9nvzGxZomuJNjPLNbP7wq/jjYmuJ1YG8mvYGwVAH5jZvWZWY2abThhfZGbbzextM/vGydbhnNvqnLsF+DCQlIehRanPR51znwU+BVwfw3JPWZT63O2cuym2lUZPH3u+GlgWfh2vjHux/dCXPlPtNYwmBUDfLAEWdR0wMy/wK+BSYBpwg5lNM7MZZvbECT9Dw4+5Evgr8EJ8y4/YEqLQZ9i3wo9LRkuIXp+pYgkR9gyMBvaFFwvEscZoWELkfaatpL4YXLJxzr1iZuUnDM8H3nbO7QYwsz8BVznn7gQu72E9y4HlZvYksDR2FZ+aaPRpoetx/xB42jm3LrYVn5povZ6ppC89A1WEQmA9KfZmsY99bolvdckjpV7UJDWKd98lQeiXZlRPC5vZQjO7y8z+G3gq1sVFUZ/6BL4AXAhca2a3xLKwKOvr61lsZr8FZpvZN2NdXIz01PPDwDVm9hticwmQeOu2zwHyGp4SbQH0X3fXo+7x7Drn3ApgRayKiaG+9nkXcFfsyomZvvZZB6RSwHWn256dc83Ap+NdTAz11OdAeA1PibYA+q8KGNNlfjRQnaBaYkl9Dlzp0nO69BkxBUD/rQYmmdk4M8sEPgIsT3BNsaA+B6506Tld+oyYAqAPzOx+4HVgiplVmdlNzjk/cBvwLLAVeMA5tzmRdfaX+hxYfXaVLj2nS5/9pYvBiYikKW0BiIikqaQ+CqikpMSVl5cnugwRkZSydu3aWhfBdwLHNQDMbBHwc8AL3OOc++HJli8vL2fNmjVxqU1EZKAws4pIlovbLiCdhi0iklzi+RnA8dOwnXMdwDunYYuISALEMwD6eikBERGJoXgGQESn2JvZzWa2xszWHD58OA5liYikp3gGQESnYTvnFjvn5jnn5pWW9vohtoiInKJ4BoBOwxYRSSJxOwzUOec3s3dOw/YC96b7adgiIokU1/MAnHNPkVrXwBcRGbB0KQgRkTSlABARSVMKABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTSlABARSVMKABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTSlABARSVMKABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTQVlwAwsx+b2TYz22hmj5hZYTyeV0REehavLYC/ANOdczOBHcA34/S8IiLSg7gEgHPuOeecPzy7Ehgdj+cVEZGeJeIzgM8ATyfgeUVEpAtftFZkZs8Dw7u56w7n3GPhZe4A/MAfT7Kem4GbAcrKyqJVnoiInCBqAeCcu/Bk95vZJ4HLgQucc+4k61kMLAaYN29ej8uJiEj/RC0ATsbMFgG3Ax9wzrXE4zlFROTk4vUZwC+BwcBfzGy9mf02Ts8rIiI9iMsWgHNuYjyeR0REIheXADhV2w42suDOF/5u3Mx6fazXY2T6PGR6PaHbrtNdx3oaD08PyvQyKNNHbpb33elMHzmZXnKzvORkeCOqR0Qk2SR1AORl+XjfxJL3jHX3qXB3HykHgkE6AkE6/EHa/aHblg4/Da2h6eM/gdD9neFlg3382NkMcjK85GdnUDgo/JOTSVFuBgU5mRQOyqBoUNfp0G1BTgbZGd6+PZmISBQldQCMLsrhx9fNiutz+gPvDY7WjgDNHf7wbYCWdj8tHQFaOkK3zR0Bmtv9NLZ1Ut/SydGWTnYdbqKhspOGlg46Az0nSn62j+EF2QwvyGF4fhbD88PTBVkMy89mVGEOhYMy49i9iKSTpA6ARPB5Pfi8HqLxd9c5R0tHgIbWTuqbOzja2klDSyf1LR00tHRQ09jOgaNtHDrWxrYDxzjc1P53WzP52T7KS3IZW5zL2CGDGFs8iPKSXCYPHUzBoIz+FykiaUsBEENmRm6Wj9wsH6MKc3pdvjMQ5HCXUNhf30rFkWYq6lrYsK+Bp946QKDLPqrh+dlMHj6YqcMHM2XYYKYMH8zEoXnatSQiEVEAJJEMr4eRhTmM7CEsOgNB9te3sqe2mR2HGtl+sJFtBxtZsruODn8QAI/BpKGDmV1WyBljCpldVsTEoXl4PfqgWkTey05yUm7CzZs3z61ZsybRZSQ9fyDI3roWth9sZPvBY2yoOsr6fQ0cbe0EQh+mzxxdwBljCpk7togzxw0hP1u7j0QGKjNb65yb19ty2gIYAHxeDxOH5jFxaB6XzRwBhD5/2FPbzPp9Dazf18CblQ0sfmU3/qDDYzBjdCHnjC/mnAnFnFlexKBM/VcQSTfaAkgjbZ0B3qxs4PVdtby+u443KxvwBx0ZXmPW6ELeN7GE86cOZeaoAjzaZSSSsiLdAlAApLGWDj9r9tbz2q46Xt9Vy8b9R3EOinMzOW9yKQunlHLepFKKcnUoqkgq0S4g6dWgTB/nTS7lvMmlABxp7uCVHYdZsb2GFdtreOTN/XgMzhhTyMIpQzl/ylBOH5mvrQORAUJbANKtQNCxsaqBl7Yf5uXtNWyoOgpASV4WC6eU8sGpQzl3Uok+TBZJQtoFJFFV29TOy9sPs2JHKBCOtfnxeYx55UV8cGpo62Di0DxdF0kkCSgAJGb8gSDrKht4aXsNL22rYdvBRiB06Y53wuCcCcU6IU0kQRQAEjfVDa3hMDjM396upbUzQJbPw4IJxaFAmDqU0UWDEl2mSNpQAEhCtHUGeGPPEV7aVsNL22uoqAt9AdykoXmcP3UoCyYUc2b5EHKzdPyBSKwoACTh3jkZ7cVtNazYfphVe47QEQji8xhnjClkwYRizplQwuyyQu0uEokiBYAkndaOAGsr6vnbrlpe21XHW1UNBB1k+TzMHVvE2eOLmTe2iDPKCnVmskg/6DwASTo5mV7OnVTCuZNCX/JzrK2TVbuP8NquOl7bVct//mUHEPo2t2kj8pk7toi5Y4uYV17EiILer6YqIn2jLQBJGkdbOlm3r561e+tZW1HP+n0NtHYGABhZkM3c8iHMKStkTlkRp43IJ9PnSXDFIskpKbcAzOyrwI+BUudcbTyfW5JfwaAMzg+fcQyhy19vO9DImoojrK2oZ83eIzy+oRoI7TaaMaqAOWOLmFMWuuz1sPzsRJYvknLitgVgZmOAe4CpwNxIAkBbAHKiA0dbebOygXUV9ayrrGfT/mN0BELfhTCqMIfZ4TCYU1bI6SMLtJUgaSkZtwB+BnwdeCyOzykDzIiCHEbMyOEfZoQue93uD7Cl+hjrKht4s7KeNysbeGLjAQAyfR6mj8xnTlkRc8YWMbusUJ8liHQRlwAwsyuB/c65DbpUgERTls/L7LIiZpcVAeMAOHSsjTcr61kX3lL4/coK7vnrHgBGFGQzp6zo+JbC9FH5ZPl0CKqkp6jtAjKz54Hh3dx1B/AvwMXOuaNmtheY19MuIDO7GbgZoKysbG5FRUVU6pP01eEPsvXAMdaFQ+HNynqq6lsByPR6mDYydMTRuZNKOGe8LmEhqS9pzgMwsxnAC0BLeGg0UA3Md84dPNlj9RmAxErNsbZQGOyr582KBjZUNdDuDx6/hMX54WsajRmiS1hI6kmaAPi7J+xlC6ArBYDES0+XsJgxqoArZo3gspkjGVWozw8kNSgARPphT20zz285xOMbq9kY/i6EeWOL+PC8MVwxayQ5mdpNJMkraQOgLxQAkgz21jbzxMZqHnlzP7sONzM428c1c0bzsbPHMnFoXqLLE/k7CgCRKHPOsWrPEf74RiVPbzpAZ8Bx8bRhfP78iZwxpjDR5Ykcl4znAYikNDPjrPHFnDW+mNqmafzv6xUseW0vz235G+dOLOHLF01m7tiiRJcpEjFtAYj0Q1O7n6VvVLD4lT3UNrVz6fTh3L5oKuUluYkuTdKYdgGJxFFzu5+7X93N4ld20+EP8olzyvnniyeTpy++kQSINAB0oRSRKMjN8vGlCyez4qsLuW7eaP7ntT1c+NOXeXbzSU91EUkoBYBIFA3Nz+bOq2fy0K0LKByUwef+dy3/5741HDzalujSRP6OAkAkBuaUFfH4F87lm5dO5a9vH+bin73MY+v3k8y7XCX9KABEYiTD6+FzH5jA0188j4lD8/jin9Zz29I3OdLckejSRAAFgEjMjSvJ5cFbFvD1RVN4bstBLv7ZK7yw9VCiyxJRAIjEg9djfH7hRB77v+dSkpfJTfet4fZlG2ls60x0aZLGFAAicTRtZD6P3fY+bl04gQfX7uPSn7/Kyt11iS5L0pQCQCTOsnxebl80lQdvOQevx7jh7pX82xNbaOsMJLo0STMKAJEEmTt2CE9/8f187Kyx3PPXPVz+i7+ysaoh0WVJGlEAiCTQoEwf3//QdH7/mfk0tfn5x1+/xn/+ZQftfm0NSOwpAESSwHmTS3n2S+dx5ayR3PXCThb916u8vONwosuSAU4BIJIkCgZl8LPrz+D3n5kPwCfvXcWtf1jLviMtvTxS5NQoAESSzHmTS3nmS+/na5dM4aXtNXzwpyv47vLN1DTqchISXboaqEgSO3C0lbteeJsH1uwj0+vhEwvG8ukF4xhekJ3o0iSJ6XLQIgPIntpm/uv5HTy+oRqPGVfMGslN545j+qiCRJcmSUgBIDIA7TvSwv/8bS9/Xl1Jc0eA00fmc82c0Vx1xkiK87ISXZ4kiaQLADP7AnAb4AeedM59vbfHKABEunesrZNH1u3noXVVbKw6is9jLJhYwkXThnHRacO0iyjNJVUAmNn5wB3AZc65djMb6pyr6e1xCgCR3m0/2MjD66p4dvNB9taFjhiaMaqABROLOXt8MfPGFjE3sk0kAAAGoUlEQVQ4OyPBVUo8JVsAPAAsds4935fHKQBEIuec4+2aJp7bcogV22tYv6+BzoDDYzB9VAEzwj/TRxUwedhgMn06CHCgSrYAWA88BiwC2oCvOudW9/Y4BYDIqWvtCPBmZT0rd9exau8RNu8/RmO7H4BMr4fxpblMKM1jfGlu6Kckj/KSXPKzfZhZgquX/og0AKL2jdVm9jwwvJu77gg/TxFwNnAm8ICZjXfdpI+Z3QzcDFBWVhat8kTSTk6mlwUTS1gwsQSAYNBReaSFTdVHeWv/UXYeamJz9VGe3nSAYJffxNxMLyMLcxhRmMOowmxGFOQwPD+b4rxMhuRmUpybRXFeJoMyvQqKFBevLYBngB8651aE53cBZzvnTnquu7YARGKv3R+gsq6FXYebqTzSTHVDGweOtnLgaBvVDa3UNnX/DWZZPg/FuZkMyctkcFYGg7N95GX7yM8OT2f5GJydQV62j8FZPrIzvGRneMjO8JKT4X3PfJbPozCJorhvAfTiUeCDwAozmwxkArVxem4ROYksn5dJwwYzadjgbu9v6wxwuLGduuYOjjS3U9fUwZHmDuqaO6hr6qC+pYPGtk4qj7TQ2Oansa2Tpnb/e7YqemMWCpTsDC/ZPi+ZPg8+r5Hh8ZDhM3weD5ne0JjP6yHTGxrzee34eIbXQ4bXg8cMrwc8Zng8hsfAa4aZ4Q3Ph8YtPE54/L3LH18mvK5QnYaF6zUsfBua5z3zXZYLL0s394UfdfxxdFnvrNGF5GR6T+1FjVC8AuBe4F4z2wR0AJ/sbvePiCSf7AwvY4YMYsyQQRE/xjlHS0eAxjY/Te2dNLb5aesM0uYP0NYRCN12Bml9Z7ojQJs/SFtngNaOAJ2BIJ1Bhz8QpDPg6AwE8QccHf4gzR2B8Hh4LHzrDwbp8AcJOgg6RyDocA4CzhF0oelU8uJXPsD40ryYPkdcAsA51wF8LB7PJSKJZ2bkZvnIzfIByXFOgnOOoINAMBQIwS7zLhwYgXBQHF8mGA6T4wESunUQvg3Pd5nmxPvCz/3OY7qu4/iyXe53hO4cUZAT83+TeG0BiIgklJnhDe/ukRAdCCwikqYUACIiaSqpLwZnZoeBihOGC4CjEa7iZMt2d18kY13ne5ouof9HOcWyz+7GI52Pdp8nq7Evy0W7z67Tqdhnd2PJ9n+3t+UGyu9oIvosdM6V9lqZcy6lfghdUqLfy3Z3XyRjXedPMr0mmfvsra+TzUe7z770Gs8+T+g55fqMtLdYvKbR6DPSntTnyad7+0nFXUCPR2nZ7u6LZOzxCKajIZZ9djce6Xy0++zLOuPZZ2/Pdyri2Wd3Y8n2f7e35QbK72jS9pnUu4BSlZmtcRGchZfq1OfAky69pkufvUnFLYBUsDjRBcSJ+hx40qXXdOnzpLQFICKSprQFICKSphQAIiJpSgEgIpKmFABxZmanmdlvzWyZmd2a6Hpixcw+ZGZ3m9ljZnZxouuJFTMbb2a/M7Nlia4l2sws18zuC7+ONya6nlgZyK9hbxQAfWBm95pZTfiy1l3HF5nZdjN728y+cbJ1OOe2OuduAT4MJOVhaFHq81Hn3GeBTwHXx7DcUxalPnc7526KbaXR08eerwaWhV/HK+NebD/0pc9Uew2jSQHQN0sIfa/xcWbmBX4FXApMA24ws2lmNsPMnjjhZ2j4MVcCfwVeiG/5EVtCFPoM+1b4ccloCdHrM1UsIcKegdHAvvBigTjWGA1LiLzPtKXLQfeBc+4VMys/YXg+8LZzbjeAmf0JuMo5dydweQ/rWQ4sN7MngaWxq/jURKNPC32/3w+Bp51z62Jb8amJ1uuZSvrSM1BFKATWk2JvFvvY55b4Vpc8UupFTVKjePddEoR+aUb1tLCZLTSzu8zsv4GnYl1cFPWpT+ALwIXAtWZ2SywLi7K+vp7FZvZbYLaZfTPWxcVITz0/DFxjZr8hNpcAibdu+xwgr+Ep0RZA/3X37RI9nl3nnFsBrIhVMTHU1z7vAu6KXTkx09c+64BUCrjudNuzc64Z+HS8i4mhnvocCK/hKdEWQP9VAWO6zI8GqhNUSyypz4ErXXpOlz4jpgDov9XAJDMbZ2aZwEeA5QmuKRbU58CVLj2nS58RUwD0gZndD7wOTDGzKjO7yTnnB24DngW2Ag845zYnss7+Up8Dq8+u0qXndOmzv3QxOBGRNKUtABGRNKUAEBFJUwoAEZE0pQAQEUlTCgARkTSlABARSVMKABGRNKUAEBFJUwoAEZE09f8BVH0MeidqcwAAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"\n",
"from obspy.signal.invsim import paz_to_freq_resp\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"\n",
"sts2 = {\n",
" \"poles\": [\n",
" -0.03677 + 0.03703j,\n",
" -0.03677 - 0.03703j],\n",
" \"zeros\": [\n",
" 0.0 + 0.0j,\n",
" 0.0 + 0.0j],\n",
" \"sensitivity\": 1.0,\n",
" \"gain\": 1.0\n",
"}\n",
"\n",
"\n",
"modified_sts2 = {\n",
" \"poles\": [\n",
" -1.05677 + 1.05703j,\n",
" -0.0577 - 0.00015703j],\n",
" \"zeros\": [\n",
" 0.0 + 0.0j,\n",
" 0.0 + 0.0j],\n",
" \"sensitivity\": 1.0,\n",
" \"gain\": 1.0\n",
"}\n",
"\n",
"\n",
"def plot_paz(paz):\n",
" spec, f = paz_to_freq_resp(poles=paz[\"poles\"], zeros=paz[\"zeros\"],\n",
" scale_fac=paz[\"gain\"], t_samp=0.01, nfft=680000 // 2, freq=True)\n",
"\n",
" plt.subplot(211)\n",
" plt.loglog(f, np.abs(spec))\n",
" plt.subplot(212)\n",
" plt.semilogx(f, np.unwrap(np.angle(spec)))\n",
" plt.show()\n",
" \n",
" return f, spec\n",
" \n",
"\n",
"plot_paz(sts2)\n",
"desired_f, desired_spec = plot_paz(modified_sts2)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"... | 1970-01-01T00:03:20.000000Z - 1970-01-01T01:56:39.990000Z | 100.0 Hz, 680000 samples\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAuUAAADtCAYAAAAY9gu6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3XecFOX9B/DPlyIgxwGKig2wgGJBLCSISlBUUhSJSYwRUzAaFSNJRAkQG4giGFRsITaa2FBQRFoUAeWA0EEwHEXq0TnuuH63+/n9sbP7m623e+w5h3zer9fzYueZZ575zuwzu9+Zmx2MJERERERExDu1vA5ARERERORop6RcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykU8YmaDzGytmfnN7NaIecea2b/NbI+Z7TazB13zrjKzAlcpcvo4wZnfwMzeMrNDZrbVzH5TSRwdzGyl089cM2vpmtfbzFaYWYWZ9U9im1Luy8wGural1MzKXdOjnDY/MbMNZlZoZh+bWVPX8ieY2afOOteZWddKYuxvZnvN7ICZDTczc+ozzOwrM9tvZrlm9rmZnZugn2fNbJOzn5eYWeeI+X8ws+1mlm9mo83smMr2X6rM7CQzm+psT0mM+e3NLMuJcbmZXeyaNypiHJWb2Seu+XHfyzixxNuvdczsQzPbYWY0s+ZJbFfKfZnZGte2+M2s2DV9lZnVMrPnzeygc0z9Ldn9GCO+OWZW4up/esT8jma20Jm33cxuidNPwjFnZq3MbKaZ5Tlj+7oEMTUzs/edfbbVzHpGzP+LmW1x+hpjZvUq204R+Y6RVFFR8aAAuB3AdQAWArg1Yt6TAGYByARwOoD1ALrF6acPgC9d08MBTHOW7QTgIIA2cZatB2AbgDsA1AcwDMA81/weAG4AMAlA/0q257D7AtAfwJiIuhMB5AL4CYCGAN4CMM41/30ArwE41lnHAQBN4/T/UwBbAJwJ4GQAawHc4cyrC6AtAhcragG4D8CiBNv7GICznba/ctbbyJl3oTN9GYDGAL4AMLgaxtAJAO4BcCOAkoh5dQF8C+BPAGoD+DWAzQCOidPXMgB/TOa9THG/1gHwFwAdARBA80q26bD7AvA/AF0i6no723gigHMA7ARwTWX7MU7/cxBxzLrmnQxgB4CfOfEeD+CsOG0TjjkA8wEMcfrpDGA/gGZx+poAYKzz3p0PYDeA85x51wHY7ozXYxE4BoemezyqqKgcXvE8ABWVo73E+oIHsBTAT1zTAwBMiLP8IgB3u6Z3Aviha3ocgEfjLNsNwDeu6YYAigC0jGg3BpUn5YfdF2In5XcDmO6aPgtAsZN8ZAAoBXCya/48AL+L0/87AP7umr4DwOwY7WoDuBfAnhTexx0ALnVeDwXwL9e8awBsqsYx1CoymUTgxGB3RN06ANfFWL4tgBIAjVN5L6uwX5NJyg+7L8ROyhcA+LVrejCANyvbj3H6jzpmXfOGA3iyCu9h2JgD0AiAH8CxrjYzAdwVZ/n9ANq6pv8djAPACADDXPMuB7CjusajiopK1YpuXxGpuSzi9flRDczOBtAewERnuimA5gBWu5qtjLWs4zx3W5KFADY69alKZ1+J+t0IoAKBK6mtAeSR3OlqH9peM7vSzA7G6wsx9o2ZrUIgQX0JgSvEwfr+ZjY1VoBm1grAcQA2JFjPGWbWIPGmVruY4whATwBTSeY506m+l5Xu1xSks6+09Gtmr5jZKxHVLzq3u/zHzNq56jsAoHM7zU4zG2/ht1utMrPbIvqPOeaCsyNeJ4o5UdvIeaeYWZMEfYnId0xJuUjNNAvAX82siZPw/QGBPztH6glgBskDznQGAB/JIlebfKc+lgxnvlui9omks69k+024TpJfkXQnHpHto+Ij2Q6BW3/uAfCNq/5pkjdEBmdmdRG4beAZV1Ibaz3B+u/KOgAlZnaPmdW1wG8LgrcvRLoNgdsfglJ9LyvdrylIZ19p6Zdkb5K9XVX9AJwBoAWA/wCYZmbBvk5F4Lj8OQL7uw6A5119tSP5dkT/UWOO5CEE/go20MyOMbMuAH6E2O8fEPjM+LsFflNyIYCbXW1nAbjNzNo4cfZ16uP1JSIeUFIuUjMNQeD+328ATEfgSviOGO0ik6kCALXNzP1lm+nUR/4YroVTnxnRZ6h9IunsqxKJ+k11nZHtY7YlWQzgdQBvuq9yRnJ+gDgGwB4Aj1eynmB9XBb+I97pidpWhmQZAonhbwHsQuB++88QMY7MrBOApgj8DiFe/HCmC+LEmNR+jZTOvpKQtn5J/pdkAcliksOdfn7gzC4GMJpktvMXhicQuE++sj5jjbmeCPwuIQfAwwA+QOzPASDw25JjEbgf/00EbgPa4fQ9E8CzCHyWZAPIQuCvTbtT2W4RqV5KykVqIJKFJP9E8mSSbRG4f3aJu42ZdUDgR2WfuJbLRSABu9DV9CIAa5z555PMcMpWBH5Ed6Grz4YI3LO9NokY09ZXJSL7PROBq4+bEPgBbOOIJ3qEtreyvippawhcST05QWwvAjgFwO0k/ZWs51sn8YqL5JeuffqTRG2TQXIZyStIHo/ACVwbBH6v4NYTwAckS+PF734v48SYyn51x5e2vpJQXf0CgXu/g76OmGdIXtiYI7mJZDeSzUheC6AlIj4HgkjuJfkrkieS7IDAidYS1/xnSZ5F8hQEbt1ZQdKXQmwiUs2UlIt4xLmloD4Cx2FdM6tvZrWceaeZWXMzq+08Bq0XgOciuugJYFKMRO8tAI+YWSMz6wigO4D34oQxB0CGBR7fVw+Bq3FLSG5x4qjjxFgbQB0nxtrfQV9ukwF0NLNuzl8ABgGYSLKUZAGAKQAGOX+27w7gArhOVGLsm3vN7AwnkX/AqYOZXWRmnZ1bBRoCeAqBJ9esj9WRmQ0CcAWAmyISWgB4G8AtZnaJmTUG8I/getLN2af1gq/dj7ozswvNrJ6ZNQLwNAKJ2Neu+XUA3ILwv7YAlbyXMcTdr8566jlxAoD7dXX3FdnvQxZ4hGYbAHfCtd2J9qObc0vZdU4cx1jg0YrH4f8T4DEAepnZmRb4DcEAAJ/G6SvhmDOztmbW0Bnbf0HgB7fT4vR1lpk1dT5XbgVwFYDRzrwGTl/mbPsIBJ7wJCI1ide/NFVROVoLAl/ejChdnHlXI/BIuiIEHuN2VcSytRF4ykqsJ2k0QCDZKHD6uK2SODoAWIXAn93nwfWEDQRuyYiM8Q/V1RdiPH3Fqf8pAj80LEIgCW/qmncCAolKEQJ/mr/WNe8qAAURfQ0AsA+BxywOB2BO/WUAlgM4hMCTLGYAaOdabiDCnwJDBH6cV+AqPV3z/4DA7QOHnPe6XjWNo8h9utk17zkEkrw8AOMBNIlY9mfOGKmVynsZJ46Y+9WZtzkyzursC7GfvlILgXu7DyJwu9EDKezHUQBGucbbEuf9PgBgNoBLIvrqg8DxuR+BEzT3eF0THCdJjLkHnfpDCJxonh5vbCPwl5DdAAoBzAVwgWvecc56C539d091jEUVFZXDK8EvIxERERER8YhuXxERERER8ZiSchERERERjykpFxERERHxmJJyERERERGPKSkXEREREfGYknIREREREY8pKRcRERER8ZiSchERERERjykpFxERERHxmJJyERERERGPKSkXEREREfGYknIREREREY8pKRcRERER8ZiSchERERERjykpFxERERHxmJJyERERERGPKSkXEREREfGYknIREREREY8pKRcRERER8ZiSchERERERjykpFxERERHxmJJyERERERGPKSkXEREREfGYknIREREREY8pKRcRERER8ZiSchERERERjykpFxERERHxmJJyERERERGPKSkXEREREfGYknIREREREY8pKRcRERER8ZiSchERERERj9XxOoDDYWYG4BQAh7yORURERETE0QhADkkmu8ARnZQjkJBv9zoIEREREZEIpwHYkWzjIz0pPwQA27ZtQ2ZmptexiIiIiMhRLj8/H6effjqQ4p0cR3pSDgDIzMxUUi4iIiIiR6zvRVKerPLycixbtgy7d+9GCrf4oEGDBrj00ktx/PHHV2N0IiIiInK0OmqS8kWLFmHAgAHIz8+v0vJmhp///Ofo378/atXSQ2tEREREJH2OiqR83759+Nvf/oZLL70Uffr0wVlnnZV0Yk0SBw8exLRp0zBy5Eicfvrp+O1vf1vNEYuIiIjI0eSouOQ7e/Zs+P1+PPnkk2jdunVKV7rNDE2bNkXPnj1x7bXXYvr06THbpXI7jBydlixZAjPDHXfckfa+S0pKcPDgwbT3G2n9+vWYPHkyCgsLE7YLHg9lZWVYuHAhfD5ftcd2tPL7/dq/1aCwsBB79uzxOozvtQ0bNqC0tLTa1+Pz+fDYY49h8ODBae/b7/envc9YsrOzsWjRIlRUVHwn6xNvHBVJ+aZNm9CiRYvQj0HLysqwZs0arF69GmVlZQACB63P50N+fn7cg+yiiy7C5s2bo+rLyspw0kknwczw0EMPpS3uBQsW4PnnnwfJtByIOTk5KCkpSUNkNVtVPiRJIjc3N6p+z5492Ldv32HHNGfOHHTo0AEAMHr0aDz++OMp95Gbm4v58+dHnQDOnj0bDRo0QNOmTWNuQxBJHDhwIOX1AoF9etlll6FNmza4+eab8atf/Spu25KSElxwwQX40Y9+hHr16uHyyy9P24lIeXn5YX2Jr1+/Hh988AFIYv/+/Ukls/FOuFetWoVrrrkGEyZMqHI8lSkqKsLo0aPxwQcfoKKiImpsHzhwABkZGTjjjDNQXl4OANi7dy+effZZTJkyBZs3b0ZeXl5K6ywuLg6b/vrrrzFixAhs2LAhqeXz8/Ph8/lQUlISWvfKlStxww03YNKkSaF2wXiTRbLSix+pXBxZvHgxevTogQULFoTV79u3D126dEFGRgZOOumkIyYxLyoqwltvvYWJEydGjROSVT5uUv3uSfY9ePrpp9G6dWvUr18fq1evTmqZPXv2hH2H5ebm4vzzz0f79u2xe/fuqM+/Q4cO4eKLL0adOnUwePBgPPbYY+jevTvMDN26dcOhQ+EPxsjJycGYMWOS2le7du1C8+bNUbt2bbRq1SqUS0TauXNnyvve5/NhzZo1oX05c+ZMnHPOOejYsSPq1q2b0rHoVlJSgnfffRcffvghdu7cmVJMkRYsWIDHHnsMBQUFh9VPOvzvf//DqFGjDitPys7OxgsvvFDl48Tn82HevHnIyso6vBO14AfdkVgAZALg3Llz+eWXX7K8vJylpaU8dOgQSXL9+vU85ZRTCIDXXHMNi4uLWVFRwcWLF4eVVatWhU1v27aNPp+Pfr+fbhMmTOCpp57Kzp07c8KECfT5fNy6dSvbt29PAKFSq1YtLl++PBSHW2SfbsXFxVy4cCGfeeYZXnHFFWF9AmCvXr2ilsnLywv1mZ2dzalTp3LZsmXs2rUrf//733PdunX87LPPeP755xMAMzIyYsZQXFwc1e/kyZND6/71r3/NoqKiuLGTZE5ODh999FFu27YtYbvK5Ofnc+3atRwyZAhfffXVsHm7d+9mVlYW+/Tpw+HDh4faP/bYY/zpT3/KX/ziF6GYZ8yYwalTp7J58+Y89dRT2aNHD86ZM4ejRo0KvTezZ8/m7373u9AyzZs3Z0lJCUnyjjvuCNVPnz6dv//97wmALVq04GOPPcYpU6bwn//8J3fv3s2TTz6ZAPjggw9yxYoVnDx5cmgdpaWlUe8lAHbp0qXSfRq0cePGqOWHDx/O/Pz8sLoZM2aELbdw4ULOmjWLf/3rX0NtevfuzQULFnDdunVR6/nqq6/YuXPnUNuf/vSnXLx4Md99992o9X/xxRdRyw8dOjTmtgLge++9F9X+7bff5llnnUUAPO6447h06dLQvPLycn777behsfnFF1+E+urYsSM///xz5uTksF+/frz33nv5n//8h7/4xS/48ccfkyQXLFjA999/P9TfN998EzOuOnXqhPbFu+++y65du3LChAk888wzQ22WLFkSFveGDRvC+giOGZIcP348e/fuzY0bN0Ztr8/n43PPPRdab6dOndi3b1++8sor/Oabb0LtgsdocMy5S3Z2NqdPn87HH388rL558+ZRcQVLPJs3b+Y//vEPvvTSS2zevDnbtGlDAOzXrx9Jslu3bmH9zJo1K2Y/Bw4c4Pjx4zlp0qS473+wPP3006HX27dvp9/v56xZs/iXv/yFI0aMYGlpadh+2L17N/1+P6+88kqedNJJ3Lt3L0ly3759fOedd7h8+XJu3ryZrVq1ClvPgAEDOGLECC5atIjl5eVRMbvbvvbaa/T7/fzggw9ixnz88cdzyZIlodj8fj8//vhjZmdnkwx8jjz99NM8dOgQe/bsyUmTJoVtQ/D93Lx5c+j1hx9+GHacNWrUiC+88ELC74igTZs2sUOHDrzllltYXl7ORYsWhcXbr18/Hjp0iCUlJXz22WdD9SNGjGCTJk1Yv359zpw5k+Xl5dy/fz9JsqKiIhQvGTg2g8v17ds3Kobc3Fz+61//4ubNm0kGxnbws3To0KFhbX0+H2fOnMnXXnuN8+fPZ9OmTWN+nsTadr/fz8WLF4eOGwB88803mZOTE/O9atOmDefOncsNGzbwmmuuqXQ87t27l6NHjw77jPvxj38cd9+/9957bNKkSVQ/b7/9dtQ2N2/enADYqlWrsHlvvPEGx44dy6ysrJjfAXfeeScB8Mknn+T+/ftjxj1u3LiY8R08eJB9+vQJtRs1ahRfffXVmH3cdNNNHDx4cNjnV1BeXh6Li4v52WefsaioiCtXruTrr7/O999/P6yPu+66K2rZpUuXcsKECSwsLAwdex9++CFvvvlmTp8+ne+//z4/+ugj7t69O+5+dtuwYQOvuuoqvvzyy9yxY0dYjrFjx45QLCNHjkyqv7Fjx7Jly5YcNmxY1D659NJLuXPnTu7atSt0TJDRuduuXbv4+eef0+fzhX2/Bj97Vq1aFZzOZCp5bSqNa1qBk5THKt9++y1btGgRmr7mmmuikvFEZdmyZaHXS5Ys4YoVKzhs2DBmZGRUepC7S05ODslA4tixY0cC4LHHHss33niDWVlZfOihh7hq1aqwD85EpXnz5rzpppt45ZVXcsGCBaH6Dh06JB3T7bffTjLwoXH99dcnvVyXLl3CBuV7771HADz33HO5adOmsLbXX389d+zYwYEDB/KGG25gr1692KdPH3722WdcuXIl8/Lyog6UJUuWxIxn4cKF9Pv97N+/f9S8n//85ym9H8Fyxx13RH2RVUf5+OOPefnllydsM27cOPr9fm7cuJGdOnVi+/bt2b17d65YsYIHDx5kYWFhtcU3depUkoEviaos7/4i//Of/5zUMmeccUbC+Z9++imLi4vD6twJcjJl+vTpodexTnBjlfHjxyecv3TpUubl5fG6666LOf+cc87h66+/HlYXecIe68vcXdatW8cmTZqwTZs23Lt3b9reZ7/fzzFjxhAAL7nkEi5btizs8yOV8sILL/Cee+7hfffdl5bYevXqFVX3+uuvxz3BqGoZPXo0SfLJJ588rH7mzp0b9jkVeTIQLF9//TX9fj8vvfTSsPr777+/0ve2S5cunD17Nh999FE2adKEDz30EP/3v/9x1apV/OMf/5jW/RKrtGvXLmb9E088wZUrV0advAQvfsUqN954Y1TSkqhEnkAFvzerqwRPRCPLV199xQMHDtDn85EMnAgdf/zxcfs5/vjjSQZO7Js1axYV9xNPPBF2IhYsl112Gf1+P0tLS9mlS5eUYl+9ejU3bNjAyZMnc+fOnVy6dGmV90OPHj24ZcsWlpeXc9asWSktO2XKFP7oRz/iww8/zOeffz6lZTMzM/nSSy+FTniD+9vt3HPPjVquX79+Mfu74YYbCIA/+MEP+Nxzz3Hbtm3Mzc1lQUEBSXL58uUpxTd//vzQMdehQ4eUtw8pJuXGFP7cV9OYWSaApP42e80112DYsGGHtb7Zs2fjiSeeSPnPNSNGjEDfvn0Pa93p1rJlS5hZzNtxEuncuTPWr1+Ps88+G19++eVhx9G1a1f06NEDWVlZeOedd+K2e+qppzBw4MDDXl9Qs2bNUFpaGvXnS6906dIFK1as+E7uC49044034pNPPqnSsjfddBMuuuiiarlXU+T7qmXLltiyZYvXYdRIF154Ia6++mps2rQJBQUFmDNnjtchJe3qq6/GF1984XUYR6yTTz457LaaF198ES1btkT37t3T0n+vXr0wevTotPSVgsYkk37sn5LyFFQ1KRcRERGRo05KSflR8UNPEREREZGaTEm5iIiIiIjHlJSLiIiIiHhMSbmIiIiIiMeUlIuIiIiIeOyoTsr79OkTejxO586dsXz5cgDA8uXL0blzZy9DExEREZGjSB2vA/DSCy+8EHo9b9680OuLL744bFpEREREpDod1VfKRURERERqgmpJys3sBDP71MyKzGydmXWtpH0rMys2s1HVEQ8AHMn/SZKIiIiIfL9V1+0rLwPIAdAMwPUAJprZWSRz47R/DsCyaooFAJCbmwuSMLMq95GXl6fkXkRERETSLu1Xys0sA8BNAB4lWUTyIwBfA7gxTvtuAAzAf5Lou56ZZQYLgEbJxpWTk4O1a9cm2zxKSUkJvvzySxQWFla5DxERERGRWKrjSnlrAHkkd7rqVgI4P7KhmR0D4BkANwO4PYm+BwB4rCpB7d+/H48//jiuvfZanHbaaahVK7nzEZIoKCjA/PnzsXLlyqqsWkREREQkoepIyjMA5EfU5QNoEqPtAwCmkdyQ5G0lQwE865puBGB7Mgv6fD5s3rwZ48ePR926dZNZJMTv96OoqCilZUREREREklUdSXkBgMyIukynPsTMTgVwB4BLk+2YZCmAUlcfKQdXWlqK0tLSyhuKiIiIiHxHqiMpXw+gsZk1J7nLqbsIwOsR7ToAOB3Aeie5zgBQy8xakfxxNcQlIiIiIlIjWXU8TcTMJgI4AOCvAK4DMAZA2NNXzKwegKauxR4EcAKAB0juT3I9mQDy0hS2iIiIiEi6NCYZeUt3XNX1nwf1RuAq+H4A/wRwC8lcM+tpZmuAwK0oJHcFCwK3txQnm5CLiIiIiHxfVMuV8u+KrpSLiIiISA1VI66Ui4iIiIhIkpSUi4iIiIh4TEm5iIiIiIjHlJSLiIiIiHhMSbmIiIiIiMeUlIuIiIiIeExJuYiIiIiIx5SUi4iIiIh4TEm5iIiIiIjHlJSLiIiIiHhMSbmIiIiIiMeUlIuIiIiIeExJuYiIiIiIx5SUi4iIiIh4TEm5iIiIiIjHlJSLiIiIiHisWpJyMzvBzD41syIzW2dmXeO0e9bMNpnZITNbYmadqyMeEREREZGarLqulL8MIAdAMwB/BzDRzJrGaJcH4HoAjQEMA/CRmTWqpphERERERGokI5neDs0yAOwH0IrkTqduHoDXSY6rZNkdALqTXBpnfj0A9VxVjQBsT0vgIiIiIiLp05hkfrKNq+NKeWsAecGE3LESwPmJFjKzVgCOA7AhQbMBCFxdDxYl5CIiIiJyxKuOpDwDQORZQb5TH5OZ1QUwFsAzJPMS9D0UgVtdguW0wwtVRERERMR7daqhzwIAmRF1mU59FDMzAGMA7AHweKKOSZYCKHUtexhhioiIiIjUDNVxpXw9gMZm1txVdxGANXHavwjgFAC3k/RXQzwiIiIiIjVa2pNykgUApgAYZGYNzKw7gAsAfBLZ1swGAbgCwE3OVXARERERkaNOdT0SsTeA0xF4Css/AdxCMtfMepqZ+4r5owDaAsgxswKn9KymmEREREREaqS0PxLxu2RmmQg8hUVEREREpCbx/JGIIiIiIiKSAiXlIiIiIiIeU1IuIiIiIuIxJeUiIiIiIh5TUi4iIiIi4jEl5SIiIiIiHlNSLiIiIiLiMSXlIiIiIiIeU1IuIiIiIuIxJeUiIiIiIh5TUi4iIiIi4jEl5SIiIiIiHlNSLiIiIiLiMSXlIiIiIiIeU1IuIiIiIuIxJeUiIiIiIh6rlqTczE4ws0/NrMjM1plZ1zjtGpjZW2Z2yMy2mtlvqiMeEREREZGarE419fsygBwAzQBcD2CimZ1FMjei3SAAxwE4FcAFAKaZ2VKS2dUUl4iIiIhIjWMk09uhWQaA/QBakdzp1M0D8DrJcRFtdwLoQXKRMz0OwAaSg5NcVyaAvHTGLyIiIiKSBo1J5ifbuDqulLcGkBdMyB0rAZzvbmRmTQE0B7A6ot0P4nVsZvUA1HNVNTrsaEVEREREPFYd95RnAIg8K8h36iPb+UgWVdLObQACV8aDZfvhhSoiIiIi4r3qSMoLAGRG1GU69ZHtapvZsZW0cxsKoLGrnHZ4oYqIiIiIeK86kvL1ABqbWXNX3UUA1rgbOT/63AXgwkTtIpYpJZkfLAAOpS9sERERERFvpD0pJ1kAYAqAQc4jD7sj8GSVT2I0fwvAI2bWyMw6AugO4L10xyQiIiIiUpNV138e1BvA6Qg8heWfAG4hmWtmPc3MfSX8UQTuDd8JYCKA3iTXVVNMIiIiIiI1Utofifhd0iMRRURERKSGSumRiNV1pVxERERERJKkpFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY2lPys2sg5mtNLMiM5trZi3jtDvRzN43s91mdsDMPjGzFumOR0RERESkpktrUm5m9QBMAjASwHEAFgIYH6d5QwBfATgfQHMAGwCMTmc8IiIiIiJHAiOZvs7MugF4nmRbZ7ohgL0A2pLcUsmyrQEsI9kohfVlAsg7jJBFRERERKpDY5L5yTauk+aVnwdgdXCCZKGZbXTqEyblADoBWJOogXMlvp6rKukEXkRERESkpkr3PeUZACLPCPKd+rjM7HQaT+R6AAAWFElEQVQATwN4uJL+ByBwZTxYtlctTBERERGRmiOlpNzMZplZSZzyMIACAJkRi2U69fH6PA7ADABPkfyskhCGAmjsKqelEr+IiIiISE2U0u0rJK9PNN+5p/xPrumGAM4CsDZO+wwA0wB8TPLFJNZfCqDUtXxygYuIiIiI1GDpvn1lDoAMM/uDc//3wwCWxPqRp5kdg8CTWtaQHJjmOEREREREjhhpffoKEHhOOYA3ALQGsBjAb4NJuZmNAgCS95jZjxBI4osAuIM4j+TWJNelp6+IiIiISE2U0tNX0p6Uf5eUlIuIiIhIDZVSUp72/9FTRERERERSo6RcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykVEREREPJb2pNzMOpjZSjMrMrO5ZtYyiWVuNTOa2a3pjkdEREREpKZLa1JuZvUATAIwEsBxABYCGF/JMg0BPAxgTTpjERERERE5UqT7SnkXAAUk3yRZAmAwgMsquVr+CIA3AOyrrHMzq2dmmcECoFE6ghYRERER8VK6k/LzAKwOTpAsBLDRqY9iZm0A/ATAS0n2PwBAnqtsP5xgRURERERqgnQn5RkA8iPq8p36WEYC+DvJ8iT7HwqgsaucVpUgRURERERqkjqpNDazWQA6x5k9BEABgMyI+kynPrKvmwBUkJyR7PpJlgIodfWR7KIiIiIiIjVWSkk5yesTzTezbgD+5JpuCOAsAGtjNL8aQGcz2+VMHwegvZm1ITk4lbhERERERI5kRjJ9nQWevrIBgR9vvgPgcQBXkIy6um5mjQA0dFVNAjAWwNskDyW5vkwE7i0XEREREalJGpOMvK07rrTeU+7cXnIzgAcAHARwBYDfBueb2SgzG+W0PURyV7AAKAOQl2xCLiIiIiLyfZHWK+XfNV0pFxEREZEayrsr5SIiIiIikjol5SIiIiIiHlNSLiIiIiLiMSXlIiIiIiIeU1IuIiIi32t9+vTxOgSRSikpT6PBg/V/Hkli/fv3R1FRkddhVFnXrl2TbnviiSfigQceqMZoRMTtsssu8zqElJx33nmoW7duta+ndevWGDlyJA4ePJj2ffTnP/85rf0lQhJDhw79ztYnHiB5xBYAmQCYqNSvX59+v58tWrQIq9++fTv9fj9LSko4btw4ZmVl8cCBAwn7AsDf/va3fPPNN6Pqp02bRpJcsWIFb7nlFk6bNq3SvmKVU045hfv27eP+/ftJkq+88gr/9re/cezYsSn3df7554dNn3DCCSn3MXLkSJ555plV2pbDKe+88w579eqVdPsrr7ySS5Ys4YABA/jCCy+wcePGBMBzzjmHDz/8MCdOnMh27drFXX7MmDH8+9//HlX/5JNP8sknnwxN//Wvf+XDDz8c1e7YY4/lW2+9xYyMDL766qts1apVzPXs27ePJPnggw9Wab8sWLCAN998M6+77rqw+rlz51a67K233sp58+bx7rvvTtjuueee45w5c2KOhaBEyzdv3pwHDhwItd2/f3+lsQ0fPpy33XZblfZJp06d+Omnn1ba7rbbbuPnn3/O2bNn84033uCaNWu4fv36hMu0adMm7vi///776fP5ePPNN6cUb4MGDfjf//6X55xzTtQxGqt07NiRK1euJMmY7/O8efNIkmPGjGGjRo04duxY+v1+fv3111Xan7Vq1Qqbvuuuu9inT5+kl2/WrBm//vprkmROTg5fe+01Ll26lDfeeGOozeDBg6sUW2Xlv//9L0ly4MCBobpGjRrFbb9371526NAhNH3PPfdw27Zt9Pv9JMlly5bx7LPPrpZYg+Wkk06q8rK333576PUpp5xCklHHwqBBg5iVlcUhQ4YQAC+55BLeeuutSfU/aNAgLly4MO3bfPzxxzM/P5/l5eVcu3Ytjz322KSXzc7O5osvvsinnnqK33zzDQ8ePMgtW7bwkUce4c6dO/mPf/wjapm1a9cyknuMuEtWVlalMcydO5d9+/bltGnTuHPnzmodHwDYpUsX7tmzJxT7/PnzK10mMucBwMLCQu7Zs4d+vz/lGE488UTWqVOHDRs2ZN++favUxxlnnBF6PWDAgJSXnzFjRtx5vXv3Trm/4HdO165duWPHDv74xz+u0vvTqVMnvvXWW8zLy+NDDz0UOT8zpbzW68Q6HUn58uXLOXToUHbq1IlDhgzha6+9FhoAmzZtCg3kKVOm8K677uLq1aujDtCgRYsWhR14paWlnDVrFufOncsXX3yRPp8v1LZ9+/YEwFdeeSWqn7Kysphv3o4dO/jLX/4y7GALvr733nvjxkWSBw4c4C9/+UvWq1cvrM9vvvmGY8aMCSV8OTk5LCoqCi03atQoPvDAA/T7/TETsrp165Iki4qK+Mwzz9DMOGzYsNDyhYWFMRNRAKxTp07o9THHHMOpU6fy1FNPTTiAhw8fnnB+8AM01kE/ffp0FhUVMTc3l4WFhXzttde4bdu2mPsr+MUaVFFREfMD190u+Hrr1q1h+9Dn87GioiI0XVJSwm3btjE7O5uLFi1iSUlJ1DpfeumlsPWsWbMmrL85c+bwggsuSPrAjzVuy8vLQ68nTZrEP/7xj5w0aRJ37tzJBQsW8J133uG5557LHj16RO2bK664Iqz/Zs2a8cEHH4zab0uXLuWnn34aVjd58uSo+F555RVWVFRELR/c9wUFBVy/fj13795NkiwtLY1qN3LkyLA+c3NzSZIzZszgo48+ykOHDsXsP9b+BsBt27Zx5cqVLCgoiLkMSY4fPz5smc2bN4fNz8nJier3+eefD4tj//793LJlC3/961/z4osv5rp16/jFF19wwYIFLC8v5+DBg1m7dm0+++yzobESVFZWxkcffZSXX345mzZtyksuuYQ9evRgw4YNQyf7kfuyuLiYZPQYj+SOvWXLlmzfvn1YkvrFF1+wuLiYu3btYnZ2dmi58vLysDgLCws5YsSI0HI33HAD58yZw4kTJ7J9+/acPn16wjjc8QQ/Q99///1QfytWrAi1iTw5zsrKos/nY1FRUdR7MGzYMGZnZ/Pbb7+N2hfLli3jhg0bSDIq6RsyZEgoDr/fz48//jjuNmzZsiXm8di9e3e2bNkyNP3tt9+yqKiI+fn5UW2feuop7tq1K6p+/Pjx9Pl8XLp0KYuLizlz5kzWr1+f/fr1Y8+ePWOu9/e//z23b98eGtPFxcWcN29e2IlwWVkZR40axaysrLjvRWlpKVevXs3s7OzQ+seNG8epU6dy1apVnDx5clj7yJM1AFyyZAmzs7PZunVrnnfeeaHje+jQoRw5ciQ//PBDXn311WzZsiUfeOABLl68mDNnzgz7Hg1K9Nl31VVXsX///jGPh1jKysp4991388wzz+SCBQvitnOv48svvwzVv/rqqzHjaNq0Kbds2RLWR3l5ecy23bp1Y0FBAffs2RN1EaVXr170+/1cvnx51HLBpDknJ4f//ve/uXjx4qi4fT5fzPcDAO+77z4ePHgwrG1hYWFYHUlu2rSJc+fOZX5+Pu+8886ofn72s59x3bp1MT+ng2655ZZQ+2nTprFfv35cvXo1Fy5cyOHDh3PWrFl8+eWX+cEHH4Q+T3w+X+hY3bZtGy+//PKkvv/c358VFRXs0aMH27Vrx1mzZtHn8zE3NzfUduvWrayoqOCECRNinvT+6le/CsVTVFQU9tmxbt06fv3111EXk1588cXQ65dffpk+n495eXkx90tFRYX7ItPRl5TH2zFVVVFRwcLCwkrblZaWcs2aNXG/GNetWxd2Jezaa68Nzfv888+5ceNGkoFBGrwynkqMN910EwcPHpzScmVlZezTpw/r1q1LALzuuuvCks1EgoNsyJAh/POf/8y+ffuGrnaUlZWFtc3Ly+N9993H119/PbR/fD5f6Itk/fr1nD17Ntu2bcuHH36Y5eXlLC0tjdqXeXl5PPPMM9m2bdvQScfheumll9iwYUOOHz8+4QfO4fL7/ZwxYwYff/xx5uTkxG3n/otA8AojSb7xxhts2rQp33777cOKw/0hGGnXrl184YUXQmMxFZFfXOmycuVKjh49OuYXdyIVFRV85513uHHjxrgnapWJt58OHDjAgoICzpo1K2GCXxMFr5q7j/M5c+ZEnXwkY/v27Xz33XfDTgYPR2lpacz3+e233+bAgQO5ffv2qHm7du1KeWxs3bqV48aNY69evXjrrbdWejITS3CZDRs2cNmyZSQDY27FihVR8ZSVlXH69OkcP348ly9fHqr/5JNPCIA//OEPw05E4tm7dy+nTJnCsrIyVlRUJP1ZXR1ycnI4depU+nw+fvTRR2FJUjo88sgjoc+S4447jm3btuWVV17JuXPnpnU9bkuWLOFVV13FUaNGhdX7/X7eeeed7Nq1KwHwN7/5TcIxU1BQwAceeIAAeNppp8Vs6/f7436HTZw4kffff3/UCXsifr+fX331Fbdt28aysjKuXr26SuM62NeXX37JVq1asVmzZkm/t36/n4cOHUrLZ+K2bdu4ceNGrlu3jllZWbz33ntZt25dtmvXjh988EFSfXz00UecOnVqzHmrV6/m2LFjw/7iUJmCggLOnDkzlN989tlnXL9+fVLL5uXlVSkp/178j555eXnIzMz0Opy4ioqK8NVXX6FTp07IyMjwOpwjDkmYmddhVKvy8nLUqVPniNvODRs2YPLkybjhhhvQtm1br8MREUmr4uJiNGjQIKm2mzdvxqmnnvqd3CcvNVt+fj4aN24MpPg/eiopFxERERFJk6om5Xr6ioiIiIiIx5SUi4iIiIh4TEm5iIiIiIjH6ngdQDrk5yd9u46IiIiISLWpal56pP/Q81QA272OQ0REREQkwmkkdyTb+EhPyg3AKQAOeR2LHJUaIXBSeBo0BuXIp/Es3ycaz+K1RgBymEKifUTfvuJsaNJnICLp5Hqm+KFUHnkkUhNpPMv3icaz1AApjzv90FNERERExGNKykVEREREPKakXKTqSgEMcv4VOdJpPMv3icazHHGO6B96ioiIiIh8H+hKuYiIiIiIx5SUi4iIiIh4TEm5iIiIiIjHlJSLiIiIiHhMSbmIiIiIiMeUlEtamNkgM1trZn4zu9VVf6yZ/dvM9pjZbjN70DXvKjMrcJUiZ/kTnPkNzOwtMztkZlvN7DeVxNDBzFY6/cw1s5aueb3NbIWZVZhZ/yS2J+W+zGyga1tKzazcNT3KafMTM9tgZoVm9rGZNXUtf4KZfeqsc52ZdU0Q301m9j8zyzOznWb2rJnVTib+GH0l3M9m9gcz225m+WY22syOqWz/SfWIN0bMrLOZzXPG1Ywk+kk0vlM97uKODzM7y8zmO+tZZmYXHc72y/dLOsazmWWY2Vdmtt/Mcs3sczM71zVf41mOGErKJV3WA/gLgP9G1P8DwBkAzgZwGYC7zawbAJD8kmRGsADoD2A+yb3OsoMAHAfgVAC3AviXmbWJtXIzqwdgEoCRzjILAYx3NckB8DCAKZVtSFX7IvmUa1seAzDBtX33mNmJAN4GcD+AEwEcctYR9LLTdzMAfwcw0Z20R1gC4CqSjQGcB6AdgD8lGX+kuPvZzC4E8CyAHgBOB9DK2XbxRrwxUgRgFIChlXWQxPhI5birbHy8A2CW09+bACabWZ3kNlWOAoc9nhF4DvldAE4AcDwCY3usa77Gsxw5SKqopK0AmAPgVtf0UgA/cU0PQCBZjbXsIgB3u6Z3Aviha3ocgEfjLNsNwDeu6YYIfLC3jGg3BkD/SrbhsPtC4ARjTETd3QCmu6bPAlAMoB6ADAS+XE52zZ8H4HdJ7POmAD4DMDyV+JPZzwh8Kf7LNe8aAJu8HmdHY0lmjAD4A4AZlfSTcHykeNzFHR8AzgGQB+AY1/wtADp7vS9VvC/pGs8RfdYGcC+APa46jWeVI6boSrl8Fyzi9flRDczOBtAewERnuimA5gBWu5qtjLWs4zx3W5KFADY69alKZ1+J+t0IoALAmQBaA8gjudPVPrS9ZnalmR10d+bU5QE4AOBCBK7cVBq/mfU3s6nO68r283kx5p1hZg1S3Xg5bAnHSCJmtsrMbnMm446PysaDmbUws4Nm1iJWXwgfH+cBWEeyzDV/VTLxylEhXeM5VAegBMBLAIY5dRrPckTRn12kus0C8FczywLQBIErH7H0ROCKyAFnOgOAj2SRq02+Ux9LhjPfLVH7RNLZV2S/eyPqgv3Wj7POJgBA8qvg6yCnrrGZnYHAfs11rSdu/CSfjogp0X6O7CvfVV8cYxul+sR7X5vEaBuGZLsk+slAJeOB5NaI9SUaH9V1HMn3Q7rGc6jOSZ5vB7DDtQ6NZzli6Eq5VLchADYD+AbAdASuhO+I0e42ABNc0wUAapvZsa66TKceZrbG/v9HlC2c+syIPkPtE0lnX5VI1G+V10nyWwSu7jyfxHpixRR3P8foK9NVL9+tdI3LysZhovFQWV/u8VFdx5F8P6R9fJAsBvA6gDedq+Qaz3JEUVIu1YpkIck/kTyZZFsAROBHiiFm1gHAyQA+cS2XC2AXArdlBF0EYI0z/3z+/48otwJY625rZg0RuGd7bRIxpq2vSkT2eyYCf63ahMAPZRubWXNX+9D2JqGWE2Os9cSNv7L9HNmXM+9b58tPvluHO0aC4o6PJMZDwr4QPj7WAjjHzOq65rerQrzy/ZSu8RzJELh6fbLGsxxplJRLWphZXTOrj8CYqmtm9c2slpmdZmbNzay2mV0HoBeA5yIW7wlgUoxE7y0Aj5hZIzPrCKA7gPfihDAHQIbzOKt6CPxifgnJLU58dZz4agOo48RX+zvoy20ygI5m1s25cjMIwESSpSQLEHiayyDnEV7dAVwA14mKm5n9KngfpHM/fn8AXyQTfwyJ9vPbAG4xs0vMrDECT9N5K4ltlTRLNEacY60+gLoAajljsm6cruYg8fhI5biLOz5IrgOwDkB/M6tnZr0B+ABkHdaOkO+FdI1nM7vIAo9QPMY5wXwKwEEEkn5A41mOJF7/0lTl+1EQeBIJI0oXAFcD2IbA0x2WIfAYP/dytRH4dfx1MfpsgMAtLQVOH7dVEkMHBH54U4zAr/hbuuY9HiO+P1RXX4jx9BWn/qcI/KiuCIEvpKaueScAmObMywZwrWveVQAKXNMDAWwHUAhgK4ARAOonGf9AhD8FJuF+RuB+9R0IPMJxDIB6Xo+3o7XEGyPOsRY5Jse4llsDoGeS4yPueAAQvL2rRTLjA4FHoc531rMcQHuv96FKzSnpGM8IPGp3uTP+9gOYAaCdq63Gs8oRU4wkRERERETEO7p9RURERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGNKykVEREREPKakXERERETEY0rKRUREREQ8pqRcRERERMRjSspFRERERDympFxERERExGP/BxiaGIk9I0BSAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 800x250 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAuUAAADtCAYAAAAY9gu6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XucTVX/B/DPNzSEcUsPueVa8iChkqhIQjc9XRT1oJJSQgr9VFTqoaIRRRRPLokSQpLKtdzvt5G7GXdjZoyZOXM5398f55z97HM/Z+Yce/B5v17rZfZaa6+99j7r7PPda+9ziKqCiIiIiIisc4XVHSAiIiIiutwxKCciIiIishiDciIiIiIiizEoJyIiIiKyGINyIiIiIiKLMSgnIiIiIrIYg3IiIiIiIosxKCciIiIishiDciIiIiIiizEoJyIiIiKyGINyIouIyFAR2SkidhHp5FF2lYiMF5GTInJCRPqbylqISJoppTvbKO8sLyYiU0XknIgcFpEng/SjqYhscbazTESqmcpeEpHNIpIjIgND2Kew2xKRN037YhORbNPyOGeddiKyV0TOi8hcESljWr+8iCxwbjNeRFoH6eNAETklIkkiMkJExJlfQkRWisgZETkrIr+JyA0B2hkpIvudx3m9iLT0KO8qIgkikioik0TkymDHL1wi8g8Rme/cn0wf5TeJyJ/OPm4SkUamsnEe4yhbRH4ylft9Lf30xd9xLSwiP4hIooioiFQIYb/CbktEdpj2xS4iGablFiJyhYh8KiLJzvdU31CPo4/+LRWRTFP7P3uU3yYiq51lCSLyuJ92Ao45EblORH4RkRTn2G4ToE9Xi8hM5zE7LCKdPcpfFZFDzrYmi0hMsP0kogtMVZmYmCxIALoAaANgNYBOHmXDACwGEAugCoC/AbT1005vACtMyyMALHSuezuAZAB1/KwbA+AIgO4AigIYDmC5qfxhAPcDmA1gYJD9yXdbAAYCmOyRdw2AswDaASgOYCqAb0zlMwFMAHCVcxtJAMr4ab89gEMAagCoCGAngO7OsiIA6sIxWXEFgF4A1gTY33cA1HLWfcy53ZLOsvrO5SYASgH4A8C7URhD5QH0BPAAgEyPsiIADgDoAaAQgCcAHARwpZ+2NgJ4NpTXMszjWhjAqwBuA6AAKgTZp3y3BWA3gLs88l5y7uM1AK4HcAxAq2DH0U/7S+HxnjWVVQSQCKCDs7/lANT0UzfgmAOwCsD7znZaAjgD4Go/bU0D8F/na1cPwAkANzrL2gBIcI7Xq+B4D34Y6fHIxMSUv2R5B5iYLvfk6wMewAYA7UzLgwBM87P+GgAvmJaPAbjVtPwNgLf9rNsWwC7TcnEA6QCqedSbjOBBeb7bgu+g/AUAP5uWawLIcAYfJQDYAFQ0lS8H8Iyf9r8FMMC03B3A7z7qFQLwIoCTYbyOiQAaO//+EMAXprJWAPZHcQxd5xlMwnFhcMIjLx5AGx/r1wWQCaBUOK9lHo5rKEF5vtuC76D8LwBPmJbfBfB1sOPop32v96ypbASAYXl4Dd3GHICSAOwArjLV+QXA837WPwOgrml5vKsfAD4BMNxU1gxAYrTGIxMTU94SH18hKrjE4+96XhVEagG4CcAs53IZABUAbDNV2+JrXacbzXVV9TyAfc78cEWyrUDt7gOQA8dMam0AKap6zFTf2F8RuUNEkv21BR/HRkS2whGgjoFjhtiVP1BE5vvqoIhcB6AsgL0BtlNdRIoF3tWo8zmOAHQGMF9VU5zL4b6WQY9rGCLZVkTaFZHPReRzj+zPnI+7/CoiDUz5TQGo83GaYyIyRdwft9oqIk95tO9zzLmKPf4O1OdAdT3LrhWR0gHaIqILjEE5UcG0GEAfESntDPi6wnHb2VNnAItUNcm5XAJArqqmm+qkOvN9KeEsNwtUP5BIthVquwG3qaorVdUceHjW9+qfqjaA49GfngB2mfL/o6r3e3ZORIrA8djAR6ag1td2XPkXSjyATBHpKSJFxPHdAtfjC56eguPxB5dwX8ugxzUMkWwrIu2q6kuq+pIp6w0A1QFUBfArgIUi4mqrEhzvy45wHO/CAD41tdVAVad7tO815lT1HBx3wd4UkStF5C4Ad8L36wc4zhkDxPGdkvoAHjHVXQzgKRGp4+zna858f20RkQUYlBMVTO/D8fzvLgA/wzETnuijnmcwlQagkIiYP2xjnfmeX4ar6syP9WjTqB9IJNsKIlC74W7Ts77PuqqaAWAigK/Ns5yenF9AnAzgJIAhQbbjyvdL3L/E+3OgusGoahYcgeHTAI7D8bz9EniMIxG5HUAZOL6H4K//cC6n+eljSMfVUyTbCkHE2lXVtaqapqoZqjrC2c4tzuIMAJNUdY/zDsN7cDwnH6xNX2OuMxzfSzgKYDCA7+H7PAA4vltyFRzP438Nx2NAic62fwEwEo5zyR4Af8Jxt+lEOPtNRNHFoJyoAFLV86raQ1UrqmpdOJ6fXW+uIyJN4fhS2U+m9c7CEYDVN1VtCGCHs7yeqpZwpsNwfImuvqnN4nA8s70zhD5GrK0gPNutAcfs4344vgBbyuMXPYz9DdZWkLoCx0xqxQB9+wzAtQC6qKo9yHYOOAMvv1R1hemYtgtUNxSqulFVm6tqOTgu4OrA8X0Fs84AvldVm7/+m19LP30M57ia+xextkIQrXYBx7PfLts9ygShcxtzqrpfVduq6tWqeg+AavA4D7io6ilVfUxVr1HVpnBcaK03lY9U1Zqqei0cj+5sVtXcMPpGRFHGoJzIIs5HCorC8T4sIiJFReQKZ1llEakgIoWcP4PWDcAojyY6A5jtI9CbCuAtESkpIrcBeBDAd366sRRACXH8fF8MHLNx61X1kLMfhZ19LASgsLOPhS5AW2Y/ArhNRNo67wAMBTBLVW2qmgZgHoChztv2DwL4J0wXKj6OzYsiUt0ZyPdz5kFEGopIS+ejAsUBfADHL9f87ashERkKoDmAhzwCWgCYDuBxEblZREoB+D/XdiLNeUxjXH+bf+pOROqLSIyIlATwHzgCse2m8sIAHof73RYgyGvpg9/j6txOjLOfAGD+O9ptebb7ujh+QrMOgOdg2u9Ax9HM+UhZG2c/rhTHTyuWxf8C4MkAuolIDXF8h2AQgAV+2go45kSkrogUd47tV+H4wu1CP23VFJEyzvNKJwAtAExylhVztiXOff8Ejl94IqKCxOpvmjIxXa4Jjg9v9Uh3OcvuhuMn6dLh+Bm3Fh7rFoLjV1Z8/ZJGMTiCjTRnG08F6UdTAFvhuO2+HKZf2IDjkQzPPnaNVlvw8esrzvz2cHzRMB2OILyMqaw8HIFKOhy35u8xlbUAkObR1iAAp+H4mcURAMSZ3wTAJgDn4Pgli0UAGpjWexPuvwKjcHw5L82UOpvKu8Lx+MA552sdE6Vx5HlMD5rKRsER5KUAmAKgtMe6HZxj5IpwXks//fB5XJ1lBz37Gc224PvXV66A49nuZDgeN+oXxnEcB2Ccabytd77eSQB+B3CzR1u94Xh/noHjAs08Xne4xkkIY66/M/8cHBeaVfyNbTjuhJwAcB7AMgD/NJWVdW73vPP49YzGWGRiYspfcn0YERERERGRRfj4ChERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZrLDVHcgPEREA1wI4Z3VfiIiIiIicSgI4qqoa6goXdVAOR0CeYHUniIiIiIg8VAaQGGrliz0oPwcAR44cQWxsrNV9ISIiIqLLXGpqKqpUqQKE+STHxR6UAwBiY2MZlBMRERHRReuSCMpDlZ2djY0bN+LEiRMI4xEfFCtWDI0bN0a5cuWi2DsiIiIiulxdNkH5mjVrMGjQIKSmpuZpfRFBx44dMXDgQFxxBX+0hoiIiIgi57IIyk+fPo2+ffuicePG6N27N2rWrBlyYK2qSE5OxsKFCxEXF4cqVarg6aefjnKPiYiIiOhycllM+f7++++w2+0YNmwYateuHdZMt4igTJky6Ny5M+655x78/PPPUewpEfkzbdo09OjRAzk5OVZ3hfxQVcTFxeGPP/6wuiv5MmrUKDRo0AAnT560uisEICsri+97uixEJSgXkfIiskBE0kUkXkRaB6jbTUT+FpHzIrJLROpEuj8HDx7EddddF/TLoElJSdi6dSvOnz/vs7xhw4Y4ePBgpLtHRCHo0qULJkyYgKlTp1rdFfLj119/RZ8+fdCqVauotH/8+HHccccdUR0D6enp6NevH7Zt24b33nsvatvxtd1wqCoyMzOj1JuCIysrC+XLl8f1118fkfZsNhvmz5+PtLS0iLRHF05OTg4OHTpkdTeiKloz5WMBHAVwNYABAGaJSBnPSiLyAIDXADwMoASA+wGcjnRncnJyUKRIkaD19u/fj6ysLOzdu9dneZEiRXxerasqcnNz891Pf8aNG4d58+ZFrf3L1fDhwzFmzBiru5FvkyZNQrt27XDuXHT/D6158+Zhz549Ud1GKDh7GTk7duzAgAEDsHv3buzbty9g3WXLlqFp06ZYv3693zoHDhyIdBfdDBw4EKtWrYrqI4QPP/yw8XdGRkbUtmM2ePBgFC9eHIsXLw55nXvvvRfFihXD6dMR/8h0k5WVFbTO+PHjMWDAgLB+QCFUu3fvRmpqKvbv34/4+HhkZ2cHXcdms8Fut/ss69evHx544AE8+uijke4qRVmHDh1w3XXXYf78+VZ3JaCNGzdi+fLleVtZVSOa4AiubQAqmvKWA3jGR921ANqE0XYMgFhTqgRAU1JS1CUpKUmPHDmiZsOGDdMuXbqoPzabTdetW2ekjRs3+qz3ww8/aNOmTb3ymzdvrtWrV9esrCy/28ir7du3KwB1vFQXj9zcXON4ZGRk6JIlSzQzM9PiXv1PYmKicVzz+rplZ2fr2bNnjeWEhAT99NNP3cajy9dff60dOnTQc+fO5bnPqo7j6sm1H4MHD85Tmz/++KPWqFFD165dq3a7Xe12u1edJUuWWD4OXdu/4YYb9Pfff4/69szH2mazuR2XCRMmaN++fb2OVW5urrZt21aff/75qPcvElzH1JWOHj0atG5sbKzfOl9++aVRz2azqaqq3W7X7Oxst3rDhw/XW265RVNTU8Pqb8eOHaM+Ds3Ho1u3blHbjq9t3nDDDWGvM27cOL91fL2XQ7Vt2zb99ttvFYD26tVLVVXXrFmjM2fO9NuX1atX53l7gfphfk3atGkTsH5KSooWKVJEb7vtNp/lhQsXtvxc5ikxMVF/+uknPX78uNVdUZvN5va5Fkx2dnZYn+3Jyclhv+9dXK/bvffem6f1g0lOTtY9e/bkux2P82qshhFDR2OmvDaAFFU9ZsrbAqCeuZKIFALQCEB9EUkQkQMi8raISIC2BwFIMSWv/82zbNmyqFKlis/ZA1VFWlqa15X20aNHveqFY9WqVThw4ADWrFkDIH+3WFwz7qdOncLs2bORmPi//wjK81ZqYmIiatasiY8++shvOxeCqmLcuHHG/gPAbbfdhmuvvRaZmZno3r077rnnHvTu3dtr3bi4ONSqVQuHDx/2Ktu7dy/++uuvoNvPy2yW+dZlXo/VTTfdhDJlyiAhwTEMmzdvjj59+qBXr15u9eLj49G9e3csWLAAI0eOzNO2AOCtt95CmTJl/M5onj171vg7Ozsbu3fvDqndjh07Yv/+/ejQoQNatGiBVq1aub0HVBVr167Nc79DdejQIbz66qs+9+/MmTPG37t37/Z6PCLc92ww06ZNQ4kSJRAXF4dvv/0WZcqUQbt27Yzy559/HqNGjfJ6dnrTpk345ZdfMGHChIj2x+zw4cNo1aqVcfcsOTnZ592DXr164bnnngur7a1bt7otnz9/HsnJyW55nr9gZZ6RNH9fJyYmBunp6WjVqhUqVqxoPJ6hqhgwYADWrl2Lzz77LOS++Zv59BTKWEhLS8Pnn39unF8zMzNx66234vXXX/fb1urVq/Hss89G9U5NXsax5zqu5fj4eFSsWBGjRo0C4D7jbbPZsHTpUthsNp9tZmdno379+njyyScBAGPHjgUA3HrrrXj88cexYcMGn+uZz0GR4vkdMNcjUv62tWTJEmRnZ2P16tUAHO+PAwcOYNWqVUhKSvJ5vl+5ciWOHDmCzZs3u5Vv2bLFGP+nTp3Cpk2b8rwfge5oVKpUCQ888AAqVKiQ78dqXM/fz5kzBx999BGSkpKC1neNmUWLFiEmJgZlypTBqVOn3Oo9/fTTuOWWW7Bz5063u7K1atVC6dKlfT5KFR8f77bfNpsNpUuXRmxsrN/P3TNnzuC5555D165dvc49LqGeC3xZvHgxHnroIa+4DwAqVqyIOnXqYObMmV5lf/zxB5YtW5bn7YYsnAg+lASgBYC9HnnDAIz1yLsWjquIBQBKA6gKYAeAbgHaDjhTPmfOHOPqxDyT5popT05O1nXr1umGDRvcrmr27dvnNlPuWe7ib6YcpquiDRs2aJs2bRSALliwwGc7/pw8eVJLly6t99xzj9sVobl9s+eee85n/o4dO7Ro0aL69ttvB92m3W7XKVOm6K5du7zyf/vtNz1x4kTQNubNm+fVD9dyzZo13frvORPryn/88ce92nWV/f333363PWHCBAWgQ4cO1Y4dO2qLFi18ziZ7KlGihNF+enq6kW+z2UKezXat//nnn7stly1b1q3eFVdcYZTlZwbV1cZTTz3lM//FF1/UPXv26I033mjk1a9fX1977bWQ2jWnpKQk3b17tw4ePNhtPALQ/fv3B+3rli1bNCEhwW/5kSNHjHHQs2dPHTRokNarV08BaJUqVbzqV6tWzauPLj169NDatWvrW2+9ZbznxowZo02aNNFTp055teVr9jAjI0MTEhJ00aJFunnzZp/HxNf4njVrlls769ati+gsXFZWli5fvtxtJuq+++4ztjF79mzjb/Ps0/nz5438N954w6tdu92ua9as8dq/hQsXakJCgg4ePFgTEhKMsXvu3Dmfx2HFihVaunRpnTx5subk5Ojo0aPd6j344IPG302bNtVFixbpP/7xD7c6TzzxhJ46dUq/+OILrxk613s5NzdX//nPf3r14ZVXXtFWrVppdna2JiQk6NixY/Xqq6/WTZs2BTyu5nPn3Llz9d///rfP11tEjHVceR07dvQ7Az19+nRdtmyZV/6RI0f0+uuv19GjR6uq6uLFi3X8+PFebdeuXdttvdWrV+vLL7/sdlzsdrv++eefxjpjx441yt577z2tWLGiHj58WNu2bWvU6d69u9v5qlu3bgpAn332WZ/7kZaW5nPsu/6ePn26W31X/oIFC/TgwYPat29fPXDggNv+P/bYY7pixQoj79y5c9qtWzdduHChzz647Nq1y+dr8/TTT/us/+OPPxp1hg4d6ve97Non87EEoC+//LKqOsY2AC1durTbPrrupKempur8+fODzhKvXbtWBw8erAD0mWee0W+++cZr/Ji3v3379oDt+bNv3z799ttvNTY21q29Nm3aeN2pcklOTtaiRYsqAD1w4IDbep53RMxl5cqV88r3fMJg//79RpndbtejR4/qwYMHfZ6vzFq3bm3U6dq1q88+tG7dOi+HyK2N+++/328ZAD148KBx/klNTTXyzfFCsG04U1gz5dEIyhsBOOmR9xmA4R55ZZwdbmnKewnArDC2FQtAly5dqrfddpvbgViyZImeOXNGx48fr4MHD9YuXbrooUOHjMBbVTUzM1NTU1N1//79AYPynJwcTUlJ0e+//94tKLfb7XrmzJmAb/qcnJygL6DLe++9F7At84ehqhonVgDavHlz483coEEDt/o2m0179Oihb775ptcbYebMmV5tr1u3Tl977TUFoMWLF/fZV7vdbjwmNGLECK82/PW/bdu2essttxiD3ZX/0EMPGe16tjFjxgw9fvy4vvfeezplyhT9v//7Pz1+/LhXoOBK27Zt0/Pnz2t6erpu375df/rpJ6/+m+unpaUZ+ZUrVw54wvDVxmeffebV5r59+4zb957989W2zWbT7du3a4MGDfSHH37QjIwMt/LMzExj/SeffNJnP3r27KmtWrXyeUz++usvr23+9ttvfj/wzBdanqlw4cI+j4fdbtcNGza4BTg2m03XrFmjGRkZmpSUpKpqBG79+/fXPXv2+NyGv4s3X+8FX/muv3v37m3UO3XqlA4aNEirVKniFqyvXbs26PvO3zYfe+wxt2OwYcMGn+/VcOXm5mrfvn2N8dipUyfjGPvrm+tD8cyZMzp37ly3shYtWuhbb72lr776qp4/f16nTZvms40uXbroTTfdpADcLoSuv/56r+NgHpMAtFGjRiEfR89Uu3ZtBRwflMnJyXrmzBlduXKlli5dWidNmqTHjh3zWsd1jvKVGjRoEPD4XnvttSH3bfPmzW6vee3atY19Xblypaqqjhw50q3NqVOn6smTJ43tdenSxSjLzs52a3/SpEnG3zVr1jTWWb16tVu97t2765o1a7RcuXJu+a7zj/mRvK5du7oF5YHSr7/+6nV8PvroI696p06dctu//fv36xNPPKGffPKJkT9//ny94YYbFIDWqVPHaM/cF5dbbrnFyBs4cKDxuexp9+7dPvvdsGFD4z3x559/6s0336zVq1cPa9zt37/f7fPL3Mc33njDbdn1d+PGjTUrK8uYMHvllVeMvu7atUtbtmypw4cP15ycHG3SpInfbZsDZXP+8uXL/Y5bu92uiYmJPssC7Wfx4sX1xIkTumfPHr3zzjt10aJFqqpatmxZv+t89913evToUR00aJD+8ssvPo+R53bNF8MzZsww8vv27asA9J133jHyzpw5owkJCQEvUADonXfeqbVq1dJ3333XyLv77rtVVXXPnj167733+rwQ9sc8fj755BOdMWOGqjomQDy33aFDB/3qq6/czoWnT5822vJ89NXzYtiZLA/KXc+UVzDl+XumPBHuQXkv5CEoL1WqlNfB/Prrr42/a9WqpU8++aTu2bPHCLxTU1ONv7dv3+4WlK9fv15tNpvxptm5c6euW7dOhw8frk2aNNHjx49rVlaW9ujRI+ib/ttvv/U5MPbu3asNGjTQCRMmqKpjJuGpp54K2p4roLPb7frss896lbdr185tuX379l51mjZtqgC8ThjVqlXT3377zav+okWL1GazaU5Ojnbr1k3Hjx+vL7zwggLQyZMn68cff2zUdV2EBNuPnTt3ej0r2L9/fwWgH3zwgd9gMa9pzZo1Pt+YgGO2s2jRolq3bl0j7/fff9eEhAQdOnSoHj16VJ966ilt0aKF26yIuY2nn37aa5tlypTR3Nxcr/z77rtPV65cqXa7XTMyMnT37t1esxsAjA+pn376yS3/iSee8LkvPXr00GbNmvk9Bj169NDNmzdrnz598n08VVVfeOEFbdeunXGBZZ59dKWHH344T+0/+OCD2qBBA83OztaTJ0/67YOvsWaeAY6JidHZs2frZ5995lanRo0aumTJEp/BXqA0a9Ysr7s/5g/XjRs3Gvn169c3Zk4Bx4dLZmamnjp1SidOnKjJycmanp6up0+f1jZt2rjNPn733Xc+x4M5uPNM69evV1XfdxXM6e233/Z5Xggn/etf/8r3GAqWSpYsafwd6CLRV6pXr57be8TzItd1sRNKmjNnTsBzWqCy5s2b65AhQ9zyfJ0rXKl69epqt9vD2t+4uDgdP358vo71H3/8oatXr9aJEyfqiRMnfNZxXawBjufYfdXxfG8sXLjQLTgDHHcfzLOh5vTbb7+pqurSpUv16quv1nfffVfj4+P99tt8zs5LatasmdvnlyvZ7XYdMGCAsewrYPNM//rXv7RWrVohb9scfHuWrVixwufs9iuvvKKA40LOUyjbdH32A8j3mBkxYoTXBSbwv/Oyr3OYr/T666+7zT6Hss6dd96pqu4TkK5Jn8TERK/Z7NTUVG3UqJHbec98B9t8YRwsVa5cWSdOnKgvvfSSFi5cWHft2qW33357oHWsDcodxxSzAIwHUAzAgwCSAJTxUW8YgPkASsLxOMs2AF3D2E6svwNh/jIHAG3VqpVb4B1qOn78uPH38OHDjcceXLfagyXPAEpVddSoUW515s+fH/KAqFq1qpYvX95nEBfN1KtXL/3+++99lplnSooWLeo2+15Q0pgxY1RVjavcYPVff/11v2UPPPCA1yMd/pKv28DmYxVo3TvvvFPT0tK8Ljrvuusu3bJli9aoUUMbN25s5D/zzDMX7Hh+8MEHxt9ffvllyCfTcNPNN9/sdRfMlVS9ZxOtSK5Zp7/++iukC2vPZP7y4r59+7RChQp56se6det006ZNQetVr17d7x2VSyXVrVtXs7KydMCAAfrYY48Zx/nuu+/WP//8U6tWrRpyW64ve/orN98dyW+67rrr9Oeffw5rnYEDB3p93oWb/L3HLnR64oknQrpjHKnk747JpEmTdNCgQVHd9rJlyzQpKcl4TMZfGjt2rD777LN6//33u+V7Ph5oxev10ksveeW5dOjQIay2/v77b83JyQm5/quvvuqVZ74QKFWqlP7666/64YcfWnJsTKlABOXlASwEkA5gD4B7nPmdAeww1bsSwAQ4vrSZCGAoAIlEUO6Z8hqUm5M5KA8nJScnq6rqokWL3D6AL7Z01113Wd6H/KTPP//8gm/T8xnbSzVF+wPMV7Lqg8hXqlOnjuV9+Ouvv7RIkSKW96OgJPOz9/lNrjuD0U5Vq1b1enae6dJMgSZ+Qkn//ve/dfz48bpo0SIdO3as5fvjSsWKFdO777477PUjRDGxAAAgAElEQVR69+6tI0eOtLz/UUhhBeXiDG4vSiISC0dAH1SrVq0wfPjwfG3v999/x3vvvRf2t6M7d+6M+Pj4gL/vS0Thee211/DJJ59Y3Q0iIsqnl19+GevXrzd+NecSUkpVU4NXc4jWfx50UejduzcmTZoEAGjZsqXxc0ebNm1Cy5YtI7adadOmMSAnijAG5EREl4YxY8ZcigF52Apb3QErjR492vjb/L8vNWrUKO//GxMRERERUZgu65lyIiIiIqKC4LIJyi/mZ+eJiIiI6NJ22QTlZ8+ezXdgnpKSwuCeiIiIiCLusgnKjx49ip07d+Z5/czMTKxYsQLnz5+PYK+IiIiIiHD5/CRioUKFUKVKFdxzzz2oXLkyrrgitOsRVUVaWhpWrVqFLVu2MCgnIiIiolCE9ZOIl01Q7hITE4MiRYqEtR273Y709PSw1iEiIiKiy1pYQfll95OINpsNNpvN6m4QERERERkum2fKiYiIiIgKKgblREREREQWi0pQLiLlRWSBiKSLSLyItA5S/zoRyRCRcdHoDxERERFRQRatZ8rHAjgK4GoA9wKYJSI1VfWsn/qjAGyMUl+IiIiIiAq0iM+Ui0gJAA8BeFtV01V1DoDtAB7wU78tAAHwa6T7QkRERER0MYjG4yu1AaSo6jFT3hYA9TwrisiVAD4C0D+UhkUkRkRiXQlAyUh0mIiIiIjIStEIyksA8PxNxlRnvqd+ABaq6t4Q2x4Ex++Su1JCXjtJRERERFRQRCMoTwMQ65EX68w3iEglAN0BDAuj7Q8BlDKlynnvJhERERFRwRCNL3r+DaCUiFRQ1ePOvIYAJnrUawqgCoC/RQRwzKRfISLXqep9vhpWVRsA43/+ca5HRERERHRRE1WNfKMiswAkAegDoA2AyQDcfn1FRGIAlDGt1h9AeQD9VPVMiNuJheMxFiIiIiKigqSUqno+0u1XtP7zoJfgmAU/A+BjAI+r6lkR6SwiOwDHrLeqHnclOB5vyQg1ICciIiIiulREZab8QuFMOREREREVUAVippyIiIiIiELEoJyIiIiIyGIMyomIiIiILMagnIiIiIjIYgzKiYiIiIgsxqCciIiIiMhiDMqJiIiIiCzGoJyIiIiIyGIMyomIiIiILMagnIiIiIjIYgzKiYiIiIgsxqCciIiIiMhiUQnKRaS8iCwQkXQRiReR1n7qjRSR/SJyTkTWi0jLaPSHiIiIiKggi9ZM+VgARwFcDWAAgFkiUsZHvRQA9wIoBWA4gDkiUjJKfSIiIiIiKpBEVSPboEgJAGcAXKeqx5x5ywFMVNVvgqybCOBBVd3gpzwGQIwpqySAhIh0nIiIiIgockqpamqolaMxU14bQIorIHfaAqBeoJVE5DoAZQHsDVBtEByz667EgJyIiIiILnrRCMpLAPC8Kkh15vskIkUA/BfAR6qaEqDtD+F41MWVKuevq0RERERE1ischTbTAMR65MU6872IiACYDOAkgCGBGlZVGwCbad18dJOIiIiIqGCIxkz53wBKiUgFU15DADv81P8MwLUAuqiqPQr9ISIiIiIq0CIelKtqGoB5AIaKSDEReRDAPwH85FlXRIYCaA7gIecsOBERERHRZSdaP4n4EoAqcPwKy8cAHlfVsyLSWUTMM+ZvA6gL4KiIpDlT5yj1iYiIiIioQIr4TyJeSCISC8evsBARERERFSSW/yQiERERERGFgUE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWYxBORERERGRxRiUExERERFZjEE5EREREZHFGJQTEREREVmMQTkRERERkcUYlBMRERERWSwqQbmIlBeRBSKSLiLxItLaT71iIjJVRM6JyGEReTIa/SEiIiIiKsgKR6ndsQCOArgawL0AZolITVU961FvKICyACoB+CeAhSKyQVX3RKlfREREREQFjqhqZBsUKQHgDIDrVPWYM285gImq+o1H3WMAHlbVNc7lbwDsVdV3Q9xWLICUSPafiIiIiCgCSqlqaqiVozFTXhtAiisgd9oCoJ65koiUAVABwDaPerf4a1hEYgDEmLJK5ru3REREREQWi8Yz5SUAeF4VpDrzPevlqmp6kHpmg+CYGXelhPx1lYiIiIjIetEIytMAxHrkxTrzPesVEpGrgtQz+xBAKVOqnL+uEhERERFZLxpB+d8ASolIBVNeQwA7zJWcX/o8DqB+oHoe69hUNdWVAJyLXLeJiIiIiKwR8aBcVdMAzAMw1PmThw/C8csqP/moPhXAWyJSUkRuA/AggO8i3SciIiIiooIsWv950EsAqsDxKywfA3hcVc+KSGcRMc+Evw3Hs+HHAMwC8JKqxkepT0REREREBVLEfxLxQuJPIhIRERFRARXWTyJGa6aciIiIiIhCxKCciIiIiMhiDMqJiIiIiCzGoJyIiIiIyGIMyomIiIiILMagnIiIiIjIYgzKiYiIiIgsxqCciIiIiMhiDMqJiIiIiCzGoJyIiIiIyGIMyomIiIiILMagnIiIiIjIYhEPykWkqYhsEZF0EVkmItX81LtGRGaKyAkRSRKRn0SkaqT7Q0RERERU0EU0KBeRGACzAcQBKAtgNYApfqoXB7ASQD0AFQDsBTApkv0hIiIiIroYiKpGrjGRtgA+VdW6zuXiAE4BqKuqh4KsWxvARlUtGaBODIAYU1ZJAAn57jgRERERUWSVUtXUUCtH+vGVGwFscy2o6nkA+5z5wdwOYEeQOoMApJgSA3IiIiIiuuhFOigvAcDziiDVme+XiFQB8B8Ag4O0/yGAUqZUOW/dJCIiIiIqOAqHU1lEFgNo6af4fQBpAGI98mOd+f7aLAtgEYAPVHVJoO2rqg2AzbRuCL0mIiIiIirYwgrKVfXeQOXOZ8p7mJaLA6gJYKef+iUALAQwV1U/C6cvRERERESXikg/vrIUQAkR6er8UuZgAOt9fclTRK6E45dadqjqmxHuBxERERHRRSOiv74COH6nHMBXAGoDWAfgaVdQLiLjAEBVe4rInXAE8ekAzJ24UVUPh7itWDi+8ElEREREVJCE9esrEQ/KLyQG5URERERUQFn6k4hERERERBQmBuVERERERBZjUE5EREREZDEG5UREREREFmNQTkRERERkMQblREREREQWY1BORERERGQxBuVERERERBZjUE5EREREZDEG5UREREREFmNQTkRERERksYgH5SLSVES2iEi6iCwTkWohrNNJRFREOkW6P0REREREBV1Eg3IRiQEwG0AcgLIAVgOYEmSd4gAGA9gRyb4QEREREV0sIj1TfheANFX9WlUzAbwLoEmQ2fK3AHwF4HSE+0JEREREdFGIdFB+I4BtrgVVPQ9gnzPfi4jUAdAOwJhQGheRGBGJdSUAJfPfZSIiIiIia0U6KC8BINUjL9WZ70scgAGqmh1i+4MApJhSQl46SURERERUkIQVlIvIYhHJ9JMGA0gDEOuxWqwz37OthwDkqOqiMLrwIYBSplQ5nP4TERERERVEhcOprKr3BioXkbYAepiWiwOoCWCnj+p3A2gpIsedy2UB3CQidVT1XT/btwGwmdoPp/tERERERAWSqGrkGnP8+speOL68+S2AIQCaq2pLH3VLAihuypoN4L8ApqvquRC3FwvHYyxERERERAVJKVX1fKzbr4g+U+6cyX4EQD8AyQCaA3jaVS4i40RknLPuOVU97koAsgCkhBqQExERERFdKiI6U36hcaaciIiIiAoo62bKiYiIiIgofAzKiYiIiIgsxqCciIiIiMhiDMqJiIiIiCzGoJyIiIiIyGIMyomIiIiILMagnIiIiIjIYgzKiYiIiIgsxqCciIiILkplypSxugtRV7hwYdx6661Wd4MugMs6KN+wYQM6derkt/zzzz+/gL0hAPjvf/9rdRcM7dq1s7oLFMTcuXNRuXJlq7sRlv79+0et7Q8++CBqbV9MKlSogNjYWDz77LMXbJsNGza8YNsKR8eOHa3uQtQ88cQTQc/Ta9aswTXXXHOBehQdEyZMwKxZs/yWh/saT5kyJaz6hQsXDqu+L6tWrcp3G5cFVb1oE4BYAJrX5PLoo496lXXq1EnVsZGQUvfu3fPcD1fq37+/pqSkaG5ubr7b8pVuuOEGt+VGjRpprVq1/NZv0qSJPvHEE1Hpi7906tSpC7o9z/Sf//zHbXyE+tofPHjwgvWxevXqlhybypUr52m9sWPHRq1Pqqo2m00/+ugjv3VeffVVBaBvvvmmTp48OSr9OHv2bEjjZfHixW717rnnnpDaf/jhh/2WVa1a1e14jBgxwme9Rx55RLt27WrJ2IlEat68edA6JUuW1NmzZ2tOTo7abDZNS0vTkSNHetWbN2+eV56IaI0aNfLcvzvvvFN79uxp+XHyTOF8jrnqnzlzJmi9O+64Q8eMGaM//vijkVeqVKkLum/x8fGalJSknTp18lneq1cv43P+iy++iOi2r7nmGn3mmWf03nvv9VunTJkyEdnWV199paqqM2bMiNhrHE797Oxs4+/bbrstT/sQbJt33XXXBRs3S5Ysicjn0rRp0/T9998PVi82rLjW6sA62kH5Z599FnCQqKo+/vjjXmX//ve/VVW1Zs2aQV+YunXranp6uj711FNeZYULFzb+Hj9+fMB2+vXrZ/TJs8zfB22g1KhRI92wYYPu3r1bt23b5tVu06ZN9c033/S5rvlkFs42zcHplVdeqa1btw5rfbvd7nbMPNPo0aO1UqVKfstvueUWv2WJiYleF2Cff/65jhgxQnfv3q02m0337dvnNj7MdXv37q3XXnutV7vjx48POlb++c9/hv36mff55MmTxvIzzzyjM2bM0E2bNgVcz/MkZ7PZQtreI4884jP/9OnTbssvvPBCSO2Zj+PUqVP17NmzOnr0aD127JhbvZ07d7ot9+zZU6+//vqgbQcLIn766Se12Wxe43nRokVedUP9wDlw4IDPfkybNk3vv/9+tzLzmNm5c6dbH4YMGRL0wnfBggWalpamDz74oM/yZs2a6datW/XQoUMB37O5ubmqqrpv376QLiJXr16to0aN0vbt2+d57K5atUq/+eYbYzkzMzPkdZ955hldunSp3nLLLTp9+nQdN25c0HVc++jJs96RI0d0x44d2qVLFyNv8uTJunPnTr3nnnt05cqVRv7HH3+s5cqV027duunAgQP9bts1qeKZ/69//Ut79eoV1nEL9vrcfffd+sUXX2h8fLxXWcWKFfWdd95xG5vmcs91WrVq5TWWz549G7SPmZmZxvHdvXu3rlq1Shs2bBjS/rVp00Zr1aqlLVq0CPmY/P777zp37lxjuXbt2j5f4yFDhmizZs20Z8+ebuXm8fPTTz9p/fr18zSmXa+z3W5XVdXExESvOr/99pv27dtXhw4dmqdteKYVK1YY+zFq1Civz0hV1RUrVhjLzZs318WLF/ttL9TP9WrVqunnn3+uqqqvvfaaPvDAA5qWlqajR492qzdv3jzdvHmzpqSk+G070DYrVaqkNptNRSTkY/Kf//xH77jjjjwdT1XVJUuWhFy/bNmyXhf3x44dM16TIBMG1gblAJoC2AIgHcAyANUC1C0C4GMAJwCkAlgR5ra8gvL333/fbdZaVdVut/t9YVRVX3nlFa+y4cOHq6rqnj17vMo8TzzNmjUz2po5c6ZbWWZmpo4YMUI3bdqkqqrLly/3++L16dPHaMcVPE6bNk3PnTunqqpNmjRRANqgQQNjnbffflsB6IsvvqgbNmzQvn37eu2fmXl7jRs31jNnzmi9evV02LBhbuWjR4/2uY6vYzhmzBifb7wffvhBjx07po0aNdKqVatqiRIlFPD+EHClSZMmqarqrl279I033vAq37x5s6qqZmVl6Zo1a4z8v/76y/h7ypQpbsFojx49jL/PnDmjWVlZbsHsnDlzvI7RxIkTjXxzkLZ161a12+367rvvuvXrxIkTqqp69OhR/eKLLzQ1NVUTExO1Xr16Rp3s7GydM2eOsXzllVcqAL3qqquMvLVr12qVKlW89jshIcHtuPbv39/oq2v9rl27qqrq4cOH9e6779Y5c+Zobm6uTp061e21GT9+vFasWDHga+o5hgHoyy+/7DUWzAFInTp1dPr06Tpx4kTt3Lmznj9/Xj/66CPduHGj23pTp071Ox7Nyw0aNFBV1bp16xp55teiefPm+sMPP7i19eGHH/rcn3nz5rnVK1q0qNH/li1butX96KOPdOTIkVqjRg1dtGiRduvWze+43759uwLQEiVKBHyfxcfH6yuvvKKdOnUyPsi3bNmiw4YN04yMDK+TvTl4OHPmjNHm4cOHjfx169a5vc8Cbd/f+SDQXTLP+sHGDOC4I5Genm7M3Ddp0kRVVdevX+/zNXYlzwuvOXPmGMfJbOnSpV7rxsbGBtxHX8fDda5zuemmm9zex57rzJ8/360/I0eO1L59+7q9v9977z1NS0tTVdXMzEz9/vvvjbJHHnlEVVWvuOIKv8eucePGxt/vvPOO39fQlcz92bt3ryYlJRll//jHP9zuHJnbqlSpkqr+L0CdMGGCqqrbeVJVvS4u4uLi1G636969e/XFF1/U6dOn+zzOjRo1MtZZu3atpqWlGed9c7r//vuNfVi4cGFIwfzSpUvd9qVp06Zu2963b59++eWXbhfgZubzr6pqdna2dunSRe+77z5NTU3V999/X+Pi4vxe/LrSzJkz3dr1nFwwj8Np06YZeStWrNBhw4Zp/fr19fDhw8bkR5s2bQJOXn3zzTde+3Lo0CGf23PdRXXFDDNmzNBbb71VP/nkE/3kk0+0RYsWOm7cuKDjy5Wys7N9Hkvz+h06dPBb5tlHX/nTpk0zLvDsdrveeuutPutt2LDB+Lt9+/ZGfV91165dG3C/VFW3bt0a0jEAoKVLl9ZPP/3U5zFXVY2Liwu0vnVBOYAYAEcAdAdQFMBwAMsD1P8EwEwA5QEUAtA4zO15BeW5ublutyVc1q5d6zVz6JKUlKQPPPCAzpw5U1evXq1DhgxxmwXIysoy1unXr5/efPPNbu2MHTvW7QWaPHmyNm3aVI8cOeI1WFUdJ1EAWqtWLbeZ2d69ext1cnNz3Wa+VFVPnjypY8aMcftwPnz4sB4+fNjtJP3dd98ZFwGeXnvtNWPdN954w6vcVfbtt9965bmS62TSqFEjo86ff/6pR48eVVXVHTt26JQpU7w+WLOzszUpKclnwA1A//jjD599qVatms6YMcOrr8eOHdOsrCxVVZ09e7b26dNHc3JyNDc3V1esWKGLFy/W7OxsHTJkiL711ls+216yZInP4+SSnp6uy5cv1/3797vlnz59WnNzczU9Pd3vuuZHE1TV7WJg/fr12rp1a7cTzcaNG/XAgQM6aNAgt+PiGkezZ8/Wxx57TJOTk41t7N27V4cNG+aWZzZ//nyv8X7u3Dkj7+uvv9bChQu7BYIzZszQfv36ufVh1qxZqqrG7f3KlSurquPD2/xe8adXr15av359zcjIcMv3fD+6bgWuXLlSVR1j2VVunj30fG+oOt6n48aN01WrVumUKVOMuosWLXKrl56erklJSarqmGE338EaMWKEW92cnBy3PrZv317//PNPo3zPnj3G7JCZOcgIJisrS+Pi4nT79u2amJioqo7X6Pz58151Z8+erb/99pvRt4MHD/ps098HkZmvCQd/9T3vEN1www3arFkzY3nu3Llu42D79u3Ga71x40a3Nr/55hu9/vrr9dVXX9W3335bjx8/bpTffvvtAY/VV199pX/++adRv0aNGnr48GFt1qyZxsXF+V1vwIABxoe1p5ycHCOI8XUMPS/qXFyPujVs2NBnuWv9jh07qqr6vQPYpk0bTUlJ0QceeECbNWtmjKdixYopAC1evHjQ19K8vWuuuUY/+eQTt7pt27ZVALpr1y6jvusiQlX1xRdfdKtv/sxzvfdDYb64cLHZbJqbm+vWpiuo8tX/u+66S8uXL699+/bVdu3aGfm///67W73GjRuH3C9Vx2dqv3793D7b/PH3nrjvvvu87saYz6ee+56Tk6P9+/f3OfmjqsZ7xHyeMd8xLFasmN8+uj5Hv/zySyMvOzvb53nDF88L7e3bt2tOTo5+9dVXRp7r89WXRYsW6UMPPWR87pv5O5+4/i5UqJDfsWx+VMbf+vfdd59Rf+7cuXrDDTfo+vXrtV+/fjpw4EC/ffDc5scff+yWv2vXLp/1S5Uq5XV3wCwnJ0d//vlnt/Ly5cu7/rY0KG8LYJdpuTgcM+bVfNQtB+A0gNJhtB8DRyDuSpV8HeysrCzt0aOHfv/9914v+O7du/W2227T+fPne5UF4mq/T58+bieeZcuW+b1lGsipU6eMK3rX86LmE2YgaWlpxvb9fSj7k5ubq3PnztUxY8b4DChnzpypvXv31pycHCPPta077rhDExMTNTU1Vb/88ks9efJkWNt2SU5O1vbt2+uUKVN0xowZGhMTo0OHDvWqF+xDMT/GjBmjL7zwgs8ZuUhp06aN27i02+366KOPul18qTpmLGNjY92CGtct/zZt2uSrjwsWLPA6iZjHz4EDB4zZEFfe9OnTNSsrS5cuXarx8fFuM5d79+7V5557Tvfs2ZPnPpk1btxYRcTtosLzgyAhIUHtdrtXv4N5++239cEHH3Qby/7cfvvtCvgO9s0XLKH69ddfFXC/q3EhuS4I4+PjtW3btn5nNl13HzwfvStdurRbPfP3CeLi4jQzM1NzcnL0v//9r/79998B+2K+GPXn0KFDOnjwYLdbwoG42rvuuutCqq+qIV08mrVv314rVqzoN8jJzs7WX375xe8F8TPPPKMAdM2aNaqqxt0xV3r00Uc1PT3d7/t706ZNev/99+vWrVv1/PnzbndtfDEHq6NGjXKrm5ub67efqur2eI3Lli1bdMOGDX7X8aVp06Yh9dEcVLlMmTJF+/fvr3a73e2YuNZxTaC4ljt37hxW38Lx3XffGXfUQglQPb9TFi7z98jMQXlMTEx+diOgjIwMr+9QqbpfrPu76xCMa33z3ThV1Q4dOmj58uX177//1rvuusvnZJuqYzKpTZs2+vXXX7ut/9xzzyngeCwu1D4A3t+HMjtx4oQ+//zzxp0dVx3zo88lS5b0ehTaF9dserly5fTo0aMFIijvC2CmR942AO181L3TWRbnDM63AegYpP0hnkG464MPcDzCES2uZw/37dunmzdv1jJlyrg94pEfdrvd50yNP+fPnzf23XMGNxpcH/Cu24cXimsf586de0G3Gymu2/bBAjPXl9KiwXz17mJ+ttx8YeXK83zEJJpyc3O9Zs/9SU9Pj9q4z8nJ8TnjreqYnc7Lh2047+lIs9vtIR3X1NRU/fXXXzU7O1t79Oihb7zxhk6aNEkPHz7sVs/zy6ThsNlsevXVV2vdunXDXtcfV1+qVq0asTY92e32kC7oAq2fmppqLMfExLh9bj377LNhtXfy5Em98cYbve7muKxYsULvv/9+3bt3r9vjfaGYOHFinl9fM9eXqv2147pTu379+pDbdLXn+pL06tWr9YUXXtDTp0/nq6/BeH6vIBDPH2cIl91uN+6MZGRk6N13360AvJ6LjzTz4zUuR44ciVhQHhcXpx9++KHxGKPdbg/4SIynzMxMrV69utvjKub3VCCrV6/Wp59+WhMTEzU3N9ft+2b+ZGRkGHUOHDhgfIl30KBBXl8M9cc1SWt6BMzSoPwtABM98lYBeMxH3aecHR4Cxwx4cwDJAOoEaN/nTPnZs2d1y5Yt+TqBBmO3291mWvIyOx4pubm5WqJECS1UqFDUgjmznJwc47b6hWR+c1yswp2di7RffvnF50nk22+/NZ7hd3EFXqHOVl5o5i8J7t2794JtNycnR7t06aKfffbZBdtmQVOtWrV8BW1ZWVkRPT+7+hLOTLnVpk+f7vah7ppBj5Y//vgj5HNnTk6Ojh492u9jj6FKS0vTDz74IOBd33DPia7j9csvv+Srb3mxbds2t+91BOL5XbZwZWRkGBfyKSkpOnv27JAnLPIqIyNDb775Zq+7t6+//roOGTIkz+26jsOOHTvy20XNycmJyB1t87PqgXTv3l0ff/xxtdvtmp6ersuWLdPs7Gy3X+8xP8LozwUJygEsBpDpJw1GeDPlHQFkAShsypsJoFcY/YkF4HeG61KWmZkZ8HnmS0FCQoJu2bLF6m5c1MzfwA/GZrMZz1oXROYZ/kg9PkOhMf8MZ0HwzTffaIUKFXw+I16Qpaamak5Ojh4/ftzqrlw0XM/UR3tmPL9OnjyprVu39vuo2OXk1KlTun37dqu74WbDhg1arFgx40c8wuVvgsufvAbloo7gNiJEpC2AT1W1rnO5OIBTAOqq6iGPurUB7ABwlarmOPNmAViqqmND3F4sgJSUlBTExsZGbD+ILhXr169H06ZNAQCRfK9bIScnB0WKFAEAxMfHo06dOhb36PJRq1Yt7Nu3D0DBGUeqChGxuhsUZenp6cjIyEC5cuWs7gpd5HJzc1GoUKE8rauqiIuLQ/369dG6deug9VNTU1GqVCkAKKWqqaFuJ9L/o+dSACVEpKuIxMAxe77eMyAHAFX9G8A6AINEpLCI3AqgDYAlEe4T0WWrcePGeOONN/DFF19Y3ZV8u+KK/52uypYta2FPLj8dOnQAAFx77bUW9+R/GJBfHq666ioG5BQReQ3IAcf5pk+fPiEF5PkR0ZlyABCRpgC+AlAbjqD7aVdQLiLjAEBVezqXqwKYBOBWAAkA3lTV2WFsizPlRJeRlStXIiMjA23atLG6K5eVjIwMTJ06Fe3bt0elSpWs7g4RUYGW15nyiAflFxKDciIiIiIqSArK4ytERERERBQmBuVERERERBZjUE5EREREZLHCVncgElJTQ35ch4iIiIgoavIal17sX/SsBMevthARERERFRZcRZoAAAerSURBVCSVVTUx1MoXe1AuAK4FcM7qvtBlqSQcF4WVwTFIFz+OZ7qUcDyT1UoCOKphBNoX9eMrzh0N+QqEKJJM/3nJuXB+8oioIOJ4pksJxzMVAGGPO37Rk4iIiIjIYgzKiYiIiIgsxqCcKO9sAIY6/yW62HE806WE45kuOhf1Fz2JiIiIiC4FnCknIiIiIrIYg3IiIiIiIosxKCciIiIishiDciIiIiIiizEoJyIiIiKyGINyiggRGSoiO0XELiKdTPlXich4ETkpIidEpL+prIWIpJlSunP98s7yYiIyVUTOichhEXkySB+aisgWZzvLRKSaqewlEdksIjkiMjCE/Qm7LRF507QvNhHJNi2Pc9ZpJyJ7ReS8iMwVkTKm9cuLyALnNuNFpHWA/j0kIrtFJEVEjonISBEpFEr/fbQV8DiLSFcRSRCRVBGZJCJXBjt+FB3+xoiItBSR5c5xtSiEdgKN73Dfd37Hh4jUFJFVzu1sFJGG+dl/urREYjyLSAkRWSkiZ0TkrIj8JiI3mMo5numiwaCcIuVvAK8CWOuR/38AqgOoBaAJgBdEpC0AqOoKVS3hSgAGAlilqqec6w4FUBZAJQCdAHwhInV8bVxEYgDMBhDnXGc1gCmmKkcBDAYwL9iO5LUtVf3AtC/vAJhm2r+eInINgOkAXgFwDYBzzm24jHW2fTWAAQBmmYN2D+sBtFDVUgBuBNAAQI8Q++/J73EWkfoARgJ4GEAVANc5952s4W+MpAMYB+DDYA2EMD7Ced8FGx/fAljsbO9rAD+KSOHQdpUuA/kez3D8DvnzAMoDKAfH2P6vqZzjmS4eqsrEFLEEYCmATqblDQDamZYHwRGs+lp3DYAXTMvHANxqWv4GwNt+1m0LYJdpuTgcJ/ZqHvUmAxgYZB/y3RYcFxiTPfJeAPCzabkmgAwAMQBKwPHhUtFUvhzAMyEc8zIAlgAYEU7/QznOcHwofmEqawVgv9Xj7HJMoYwRAF0BLArSTsDxEeb7zu/4AHA9gBQAV5rKDwFoafWxZLI+RWo8e7RZCMCLAE6a8jiemS6axJlyuhDE4+96XhVEagG4CcAs53IZABUAbDNV2+JrXacbzXVV9TyAfc78cEWyrUDt7gOQA6AGgNoAUlT1mKm+sb8icoeIJJsbc+alAEgCUB+OmZug/ReRgSIy3/l3sON8o4+y6iJSLNydp3wLOEYCEZGtIvKUc9Hv+Ag2HkSkqogki0hVX23BfXzcCCBeVbNM5VtD6S9dFiI1no08AJkAxgAY7szjeKaLCm+7ULQtBtBHRP4EUBqOmQ9fOsMxI5LkXC4BIFdV0011Up35vpRwlpsFqh9IJNvybPeUR56r3aJ+tlkaAFR1petvF2deKRGpDsdxPWvajt/+q+p/PPoU6Dh7tpVqys/wsY8UPf5e19I+6rpR1QYhtFMCQcaDqh722F6g8RGt9xFdGiI1no08Z/DcBUCiaRscz3TR4Ew5Rdv7AA4C2AXgZzhmwhN91HsKwDTTchqAQiJylSkv1pkPEdkh//sSZVVnfqxHm0b9QCLZVhCB2s3zNlX1AByzO5+GsB1fffJ7nH20FWvKpwsrUuMy2DgMNB6CtWUeH9F6H9GlIeLjQ1UzAEwE8LVzlpzjmS4qDMopqlT1vKr2UNWKqloXgMLxJUWDiDQFUBHAT6b1zgI4DsdjGS4NAexwltfT/32J8jCAnea6IlIcjme2d4bQx4i1FYRnuzXguFu1H44vypYSkQqm+sb+huAKZx99bcdv/4MdZ8+2nGUHnB9+dGHld4y4+B0fIYyHgG3BfXzsBHC9iBQxlTfIQ3/p0hSp8exJ4Ji9rsjxTBcbBuUUESJSRESKwjGmiohIURG5QkQqi0gFESkk8v/t3L9qVEEUgPHvuIkbQRELW30BURGENIKCQbBIp4VptPUJQkAwjZ34CimCIIJ2YmcarcRUFrHzXylYrIqFjMVM4Lrcu15xdUz8fjDF7uydnR3OZc8uMycWgGvAnbHLl4AHLYneOnAjIg5ExDywCNzrmMIGsL+UsxqST8w/Tym9LvObKfMbADNlfoO/MFbTQ2A+Ii6Uf25Wgfsppa8ppRG5mstqKeG1CByj8UOlKSIube+DLPvxl4EnfebfYtI63wUuR8SpiDhIrqaz3uOzasomxUi51+aAWWBPicnZjqE2mBwfv3LfdcZHSmkL2AKWI2IYEdeBb8Cz31oI7QrTiueIOBG5hOLe8gPzFvCRnPSD8aydpPZJU9vuaORKJGmsnQXOAW/J1R1ekMv4Na8bkE/HL7SMuY+8pWVUxrjykzmcJh+8+UI+xX+00XezZX5X/9RYtFRfKc9fJB+q+0z+QjrU6DsMPCp9r4Dzjb4zwKjxeAV4B3wC3gC3gbme81/hxyowE9eZvF/9PbmE4xowrB1v/2vripFyr43H5FrjupfAUs/46IwHYHt715E+8UEuhfq0vM8mcLL2Gtr+nTaNeCaX2t0s8fcBeAwcb7zWeLbtmBYpJSRJkiTV4/YVSZIkqTKTckmSJKkyk3JJkiSpMpNySZIkqTKTckmSJKkyk3JJkiSpMpNySZIkqTKTckmSJKkyk3JJkiSpMpNySZIkqTKTckmSJKmy7+MSTLf0Lx3wAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 800x250 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Generate some random data.\n",
"import obspy\n",
"random = np.random.RandomState(13)\n",
"\n",
"sampling_rate = 100\n",
"seconds = 60 * 120\n",
"data = random.rand(sampling_rate * seconds)\n",
"\n",
"tr_1 = obspy.Trace(data=data, header={\"sampling_rate\": 100})\n",
"tr_2 = tr_1.copy()\n",
"\n",
"# Add two different instrument responses.\n",
"tr_1.simulate(paz_simulate=sts2)\n",
"tr_2.simulate(paz_simulate=modified_sts2)\n",
"\n",
"tr_1.trim(tr_1.stats.starttime + 200, tr_1.stats.endtime - 200)\n",
"tr_2.trim(tr_2.stats.starttime + 200, tr_2.stats.endtime - 200)\n",
"\n",
"print(tr_1)\n",
"tr_1.plot()\n",
"tr_2.plot();"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"from obspy.signal.calibration import rel_calib_stack\n",
"import scipy.signal\n",
"\n",
"# Recover the unknown instrument response.\n",
"freq, amp, phase = rel_calib_stack(tr_1, tr_2, calib_file=sts2, window_len=100.0)\n",
"\n",
"a, b = scipy.signal.ltisys.zpk2tf(modified_sts2[\"poles\"], modified_sts2[\"zeros\"], 1.0)\n",
"desired_spec = scipy.signal.freqs(b, a, freq * 2 * np.pi)[1]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA3YAAAHZCAYAAADKeKmfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XeYVOXB/vHvc6bsbGeXXbo0QaoaYWn2qCgq2FNIXt9YovE13V+KsWsw9hijJrFEscUSVBQMYsfYXWw0qVLWpZftU8/z+2OWZWnCwu6emd37c11zzcyZ55xz78DDzs2ZOWOstYiIiIiIiEj6crwOICIiIiIiIvtHxU5ERERERCTNqdiJiIiIiIikORU7ERERERGRNKdiJyIiIiIikuZU7ERERERERNKcip2IiIiIiEiaU7ETERERERFJcyp2IiIiIiIiac7vdYBvUlRUZHv37u11DBEREREREU/Mnj17g7W2eE/jUrrY9e7dm9LSUq9jiIiIiIiIeMIYs2JvxumtmCIiIiIiImlOxU5ERERERCTNqdiJiIiIiIikuZT+jN2uxGIxysrKCIfDXkdJC6FQiB49ehAIBLyOIiIiIiIiLSTtil1ZWRm5ubn07t0bY4zXcVKatZaNGzdSVlZGnz59vI4jIiIiIiItJO3eihkOh+nYsaNK3V4wxtCxY0cd3RQRERERaePS7ogdoFLXBHquRKQlhaMxtmxeT6RyI9GqjUSrNxKp3EC0toJIXS3RcB1uLExuwCXTxPDZGJG4S9y1GONgjIPjS15bXwYJX4i4k4ETzMIfyiKQkU1GZg7B7Hz8WQWQ2QF/VgeKiwrJ0FvMRUREGqRlsRMRkdZRG4lR/vVKfBu+hI1LqF2/AndLGYGacvKiaylyN9LFJPa4naj1EyZADD9g8BkLWIy1GCwOLkFiBPdiWwAJa6ggi2qTTa2TQ8SfSzhYQDxUiJtZhMkuxuQU4c/tRCi/Mx27dKdTcWd8Pt/+PSEiIiIpSsVuHxx++OG89957LbqPyZMnc+KJJ9KtW7e9Xmf58uWMHz+euXPntmAyEWmronXVlC94n81LPiZcPp+cyiUckFhJP1PTMCZmfaylkC2BTpTnHEJ5fnec3M6QVYCTVYg/pyPZHTqRl19Afm4u/mAm+IIEjMFxLQB+364/BZBwLbFEDGJ11NXUUF1TRW1NNeG6aqLVm7HhCpxIJSZcQbhqEyZSgT9WRSBWRTC2hc61i8mrriCf6l1uP24dNpg8qnwdqA0UEM3oSCKzCHI64c/rTKhDF7I7diO/qDu5Hbtg/BnN/ySLiIi0EBW7fdDSpS6RSDB58mSGDh3apGInIrK3KsMxSj//gnVfvEqHDZ/QN7qQvu5KehuX3sBmclkf6sPKDuNY1Wkg4YL+uIX96HFAH7p1yKaH07S3eRvA7/vmdXyOwecEIRAkkJVPXvG+/WzxaJiqzeuo3byGyJa1hCvWUrd5LdHKdTi1GwhENpIV20RxuJyCLVvINpFdbqeSHCp8BdQGCohkFG1XAjMKupBT2I0Oxd3JKugKgdC+hRUREWkmKnb7ICcnh+rqat566y2uu+46ioqKmDt3LsOHD+fxxx/n5Zdf5uGHH+aZZ54B4K233uKOO+5g2rRpvPLKK1x77bVEIhEOPPBAHn74YXJycujduzcXXHABr7zyCpdccgmlpaX88Ic/JDMzk/fff5/58+dz2WWXUV1dTVFREZMnT6Zr167Mnj2bCy64gKysLI488kiPnxkRSVXWWpYuX8mKj17ELH+bfrWfcZxZB0CNyaEsaxAfdTgO54ASOg0YTe+efSjYzZG1VOcPhijo3JOCzj33ONZ1LZsqtlCx4WuqNpRTt3k1sYq1uNXrcGrXEwxvICu2ifzwAgq3bCHX1O1yO9VkUenrQE2gMFkCs4ogu/5IYH0JzC/qTqigCwSzm/tHFhERSe9id/20ecwvr2zWbQ7ulse1E4bs9fhPP/2UefPm0a1bN4444gjeffddxo4dy09+8hNqamrIzs7m6aef5nvf+x4bNmxg0qRJvPbaa2RnZ3PLLbfw5z//mWuuuQZIfufcO++8A8CDDz7I7bffTklJCbFYjJ///Oe88MILFBcX8/TTT3PllVfy0EMPcf7553P33XdzzDHH8Nvf/rZZnwsRSW+1kRiflb5L5Zz/0G3tWwxxF9HPWKpMNmsKh7Oy/yV0+9ZYsrsMZYCTniVufzmOobCggMKCAug/9BvHJlzL+opKtqz7muqN9SWwci1u1Tp8dckSmBndRF54MQVbPqZDo7ewNlZHaLsjgfHMIkyjt4PmdOxGXnE3gnldICMXdBIsERHZC2ld7FLByJEj6dGjBwDf+ta3WL58OUceeSTjxo1j2rRpnHPOObz00kvceuutzJo1i/nz53PEEUcAEI1GGTNmTMO2vve97+1yHwsXLmTu3LmMHTsWSL5Vs2vXrlRUVLBlyxaOOeYYAM4991xmzJjRkj+uiKQ6a1n95Yd8Nesxeq6eyeFmPQArgv1Z0OsndC05nY79R5Hr6CQiTeVzDMUF+RQX5AODv3FsLOGytrKazetXU7W1BFaswVYn3w6aLIEbyQt/RfGWzyikCsfYnbYTIUiF04HqQCF1wUJioSLcrGJMbmcC+Z3JyO9CVkEXcjp2Jze/I06aHmUVEZH9l9bFrilH1lpKRsa2D9f7fD7i8TiQLGn33nsvhYWFjBgxgtzcXKy1jB07lieffHKX28rO3vXbc6y1DBkyhPfff3+75Vu2bNHXGYgIADWr5rBq1qN0+Go6XRPlFFkfi3OGs2Tor+g5+kx6FXT3OmK7EvA5dC7Io3NBHjDgG8dG4gnWVNayef3q+iOBa4hVrsGtSr4dNCOykezYJnKryymunE8hlfiNu/N2rJ8Kk0u1k0edL5dIIJ9YMJ9ERgfITJ7cJpBTSDC3I5l5ReR0KCK3oJisnA6YdnrEVkSkLUnrYpfKjj32WC688EIeeOCBhiNxo0eP5qc//SlLliyhX79+1NbWUlZWxkEHHbTT+rm5uVRVVQEwYMAA1q9fz/vvv8+YMWOIxWIsWrSIIUOGkJ+fzzvvvMORRx7JE0880ao/o4h4x3UtHy/8ivXvPkH/1S8wILGY/tbwie8Q5vQ7n0NO+B8Gd9HJl9JBht9Ht8JcuhXmAjv/PmjMWktNJMbmDWuo2lhOpNHbQZ3a9Th1m/FHtxCMVdIh8jXZtV+Sa6vJ2s0JYiB5ptMqk021k0utk0edP4+oP5d4IAc3mIebkQsZ+TihXHxZ+fiz8glmFZCR24FQTgFZuR3ICQXJCvr0n40iIh5SsWshPp+P8ePHM3nyZB555BEAiouLmTx5MhMnTiQSSf6SnTRp0i6L3Xnnnccll1zScPKUKVOm8Itf/IKKigri8Ti/+tWvGDJkCA8//HDDyVNOOumkVv0ZRaT1rausY9bLU8j/8mmOTnxAyMRYGejDWz1/TcGoH1AyoL9eXLdhxhhyQkFyevSEHns+OcxWkXANVZs2ULVlPXUV6wlXbSRWtZFE7WZs3Sac8BYCkQpC8Qry4hvIjKwgy9aQTS0B9vzdglU2k9VkUkM2dU4WdU4OEX8OMX82cX8ONpgNgWxMRjZORg6+UC7+UC6BzFwysnPJyMojMyePrOx8srOzyA76cZp45lURkfbOWLvze/pTRUlJiS0tLd1u2YIFCxg0aJBHidKTnjOR9Lfs69XMmf53Di7/N31NOTUmh3W9J9D12IsI9RymE2xIy7AWG6slXL2F2spNRKq3EKneQrR2C4naLSTqKrHhCkykEhOpwherwh+rIhivJiNRTWaihgxbR4joXu8yan3UEqKOEHUmk6iTSdSXSczJIu7Pwg1kgj8TAlk4gRAEMjFbLxlZ+IJZ+DKy8Acz8Ydy8GdkEgzlEAxlEczMJpSVQ0YwqP8AEZG0YYyZba0t2dO4VjtiZ4zJBv4GRIG3rLV636CIyB6sWvgJy2fcxWGbX+Z0E2ZV9mDWH345xaO+Rx99d5q0NGMwwWwyC7PJLNyPz2m6CYjWYKPVRGurqK2pJFJTSbimklhdFbG6KhLhKhKRGmykCqK1EK3GF6vBF68lI1FLbmIdwVgdGTVhgkQJ2QhBE9+nODHrI0yQiAkSIYOoySDqZJAwQRJOgISTgesEcX1BrC8IviDWlwH+jOS1LwP8QRx/EPwZGH8IE8jACYRw/Bn4Ahk4wUx8gRC+YAaBYAhfIEQgI4Q/GCIQDBIMhggGMnTCGxFpNvtV7IwxDwHjgXXW2qGNlo8D7gJ8wIPW2puBs4Ap1tppxpinARU7EZFdsZZ1n81gy2t3cFBNKZ1sgIXFY+l50i85oP9or9OJNJ3jg1AeJpRHRh5k7HmNvRKLxQjXVROpqyVaV00sUkssXEs8UkM8UkciUkMiWocbqcXG6uovtRAP48TqMPEwTqIOJxHGnwjjuFGCiSi+eA0+N4bfxgjYKH7iBG2UIDGCxHd5BtN9FbcOcXzEjZ84fuL4SBg/CfwkjI+4CZDAh2v8JIw/ee0EcOtvWyd57ToBrOPHbr02AawvgHUC4PjBF0ieJMcXwBgfxufHONuuHZ8f4/hxfNsuxvHh+AM4vgCOz4ev4bHkfRw/Pp8P4wQwvuS2fD4/xufHqV/XOD58/kDyvi+A4zj4HINj0FFTkWa2v0fsJgP3AI9uXWCM8QH3AmOBMuBjY8yLQA9gTv2wPb9hX0SkvUnE2PDhU0Tf/gvdwkvAFvBGj0s45LRfcEhnndVSZEeBQIBAoIDcvIJW22ci4RKORYlE6ohHw8QidSSiYWLRMIlomHg0jBuLkIiFcWNhErEINl5/qb/tJmKQiGLjMXBj2EQM3DgmEU1euzGMG8excRw3lry2cXz1ywLxOnw2gWPrq6BN4COO3zZUQ/zE8dvELs+g6iXXGuI4xHBwcXAxWAwJHOzW+8ZgceqXJW9bDAnjNNy2Jnnt4msYv+3a2e4+xuAaH9QvAwfXOFB/sZjktfE1Wmfrfnz1jxnAQKNtmvqsQHJb9Y/TMHb7ZaY+N422ZQCMD2vA4Gx7W73jwwC2fr3kNpz6Mly/bZNcN7nI1K+fXAYmuXuS62zNb7fL2OjnYtvPsf2Yxrb+TDQa32hM/c9jtz5Go21st63kNhqeOwzG1K9nGpYkxzQ8T9u2tXX7ptF2bf02Gv8s1pgdfoRd57eNBvXs2ZeuXbqSrvar2Flr3zbG9N5h8UhgibV2GYAx5ingdJIlrwfwGVv/REVEhERdJcte+RsFX/yTosQ6ltjuTO19BYeffgnHFeZ7HU9EGvH5HHy+EKFQmrwV2nXBjeHGoyQSceLxBIl4lHg8TiIeJ5GI1V/HcRMxEokEiXjy2o3HcBMx3ESi4TGbiOMmEuDGwCbAjWNdFxJxsAmsm7zeeh83gXETyds2gXHjGDeBtS7GuljrJpdbN7kdLMYmwLr1F4ux7rZlJO/TsCxZ78zWsSQwW9chOcap31fy/rbxZmuttNsqZsP1dsu2LadR1dxaMZLbbHyhYazTaFlzHumVlvHR0Gvpes5lXsfYZy3xGbvuwKpG98uAUcBfgXuMMacC03a3sjHmYuBigJ499/6MXyIi6Wbzpg0snnYHA756lP5U84kZzH8HXc7h4yZyRn6W1/FEpC1wHHAycPwZOEDA6zztnbVY6+K6yWtrLe7WUttw22KxuG7ycRqNtW4C11osbv1ysDa5jPrxtn4f2PptYBsKcHL/yRwNJRlbX5ppWNZwcsWGkyza7W+T/PoVGrqq3cWYxvtKll273fhtORqvaxvGs9M+t1uv0WO2fplpNK7xvhv20Xibdudl/fqN2MUfWvpoiWK3qzdMW2ttDXD+nla21t4P3A/Js2I2czYREc+tXlPOvOduYeTaZxhpapmdMYq6Mb9m1FEnMUwnUhARabuMwRgf+qdeWkJLFLsy4IBG93sA5S2wn7R13nnnMX78eM455xyvo4hIK6rYsJp5z97EIeXPcIKpY16Ho8ka+weGDz3c62giIiKS5lqi2H0M9DfG9AG+Br4P/KAF9pMSth7ydhz914uI7NqXX61k3cu3UrL2GUbbKJ/nH0vXCVczpP9wr6OJiIhIG7FfbcQY8yTwPjDAGFNmjLnQWhsHfgbMBBYAz1hr5zVxuxOMMfdXVFTsT7wWs3z5cgYNGsSll17KsGHDeOyxxxgzZgzDhg3jO9/5DtXV1QDccMMNjBgxgqFDh3LxxReTyl8GLyLN79VPF/P4LT+l2+RRHL32MRbkHs5X33udwy6bSheVOhEREWlGJpXLRklJiS0tLd1u2YIFCxg0aFDyzozLYc2cXay5H7ocDCff/I1Dli9fTt++fXnvvffo168fZ511FjNmzCA7O5tbbrmFSCTCNddcw6ZNmygsLATg3HPP5bvf/S4TJkxo9bdibveciUiLq6is4K3Hb+KotY9RaKpZUXwsBadcR16fw7yOJiIiImnGGDPbWluyp3Et8VbMdqFXr16MHj2a6dOnM3/+fI444ggAotEoY8aMAeDNN9/k1ltvpba2lk2bNjFkyBAmTJjgZWwRaUnxKMteuZfcj+7idDazvGA0eWffSK8D9vhvsYiIiMh+Se9it4cjay0pOzsbSH7GbuzYsTz55JPbPR4Oh7n00kspLS3lgAMO4LrrriMcDnsRVURamrWUf/Bv/G9cR9/Y13zhDGbLuPvoP/Ikr5OJiIhIO5GSZ/xI9c/YNTZ69GjeffddlixZAkBtbS2LFi1qKHFFRUVUV1czZcoUL2OKSAv5et67LL31aLrNvIjKqOG5QXfS//f/VakTERGRVpWSR+ystdOAaSUlJRd5nWVPiouLmTx5MhMnTiQSiQAwadIkDjroIC666CIOPvhgevfuzYgR6f2FhyKyvdjG5Sx76vcMWP8yG2w+M/r8nlFn/4qzcvXF4iIiItL60vvkKbJX9JyJNKNwBRtevpm8zx7AtfBWx+8x/IfXU9yxyOtkIiIi0gbp5CkiIs0oEYtS/sY/6PjxHRTFt/CSOZrsU65n3MhhXkcTERERUbETEdkday0vz1nNqg+f54Sv76UvX/O+O5h3+t7Kj79zJgXZQa8jioiIiABpWuystRhjvI6RFlL5rbYiqWxLbZR7npjCcSvv5mTffNYGe/LhYfcw8KjvMiYnw+t4IiIiIttJyWJnjJkATOjXr99Oj4VCITZu3EjHjh1V7vbAWsvGjRsJhUJeRxFJK5/Oncfq567gisQsIqF83BNuo3PJ+XT2BbyOJiIiIrJLaXfylFgsRllZmb4Tbi+FQiF69OhBIKAXpCJ7Equt4JMnr+eQlY/hM5aKQy6k+JQrIJTvdTQRERFpp9rsyVMCgQB9+vTxOoaItCEbK2v4Yto9HLL4b4xiC7Pzj2fgD2+nuHNfr6OJiIiI7JW0K3YiIs1lQXkF7778FEev+CvfNmUsDA7lq2MfpOTwsV5HExEREWkSFTsRaXfCsQSPvfASA7+4lR87c9iQ0YPV376fAaO/C/rsroiIiKQhFTsRaVc+n7+Asueu4sLY64QDOdQdcyNFh18Mfn11gYiIiKSvlCx233RWTBGRfVFbXcFHT1zPyPLHGWwSrB50Ad1PvxoyC7yOJiIiIrLf0u6smCIiTeHG43w6/W/0/OzPFLOZuR2Oo+/3byOri/7jSERERFJfmz0rpojI3prz9lSyZ13L8MRyvvQPZP3YBxg6SidGERERkbZHxU5E2pylcz+mevofODT8MeWmMx+PuIPh4y7A8TleRxMRERFpESp2ItJmLFyylLUvXssRFdOpMZl82O/XfOuc39EtlOV1NBEREZEWpWInImlv/so1LJp6Cyds/Bd9ifF5t+/S75w/MqpjZ6+jiYiIiLQKFTsRSVvL1lXy1pR7Gbf2fs4wm1jc8Vi6nHUzw3oM8jqaiIiISKtKyWKnrzsQkW8SiSeY/sIzDPjiFi4wX7EmdxDVEybTf8AxXkcTERER8YS+7kBE0srnn35M9UtXcET8Izb5O+Mbey35IyaCoxOjiIiISNujrzsQkTalYkM58/51JSM3TiViMlh6yG84cMJvIJDpdTQRERERz6nYiUhKi0dq+eLZWzho0f2MtGE+73wmgybeyIEFXb2OJiIiIpIyVOxEJCVZ12XuKw9T/OHNDLPrKM0YRcHpNzF88HCvo4mIiIikHBU7EUk5S0pfxc68koNjC1ns9KXs6NsYfuzpGGO8jiYiIiKSklTsRCRllC+dx/rnL+fQ6rdZSyHvHvxHRp5+KQG//qkSERER+SZ6tSQinitfU85XU65hxPrnyMfPf3v+hG9990qOyM33OpqIiIhIWlCxExHPrNtUwSfP3saYsocYTS2fFJ1G73MmcVTXnl5HExEREUkrKVns9AXlIm3bhqows55/gBFL/8o4s46FeaPocPrNjOg3zOtoIiIiImlJX1AuIq1mS22UaS+9wJC5tzLMLKI8oy/+cTfS6bBTvI4mIiIikpL0BeUikjIqwzGmvPJfus2+hXPNB1T4C1l39O10O+oCcHxexxMRERFJeyp2ItJiaiJxnpz1GaH37uB/7ExcJ8CG4f+PohP/HwSzvY4nIiIi0mao2IlIswvHEjzxzkKq3v4bF7jPkmPCbBn0fQpPvZZQbhev44mIiIi0OSp2ItKsPl+5iWlP3M154UfpYTZQ0fM4nAl/orDTIK+jiYiIiLRZKnYi0iziCZepzz/NgC9u4SrnK6oLB8OEh8jve4zX0URERETaPBU7EdlvZYs+o3zK7zgn+iGbg52oPelv5AyfCI7jdTQRERGRdkHFTkT2ma1ay+JnrqTvymfJNxksGHwZg878HQQyvY4mIiIi0q6o2IlI00VrqZ71F3zv/ZU+bpQ3cidw6A//xKCuPbxOJiIiItIuqdiJyN5zE8Q+eYLYq38kJ7KOme5Iao68kjNOOAbHMV6nExEREWm3UrLYGWMmABP69evndRQRqVc1dyaR/1xBUe0S5rr9eKbwt1z4gx/Qr1OO19FERERE2j1jrfU6w26VlJTY0tJSr2OItGvlCz+mdvoV9Kv6iBVuJ6Z3upjDTjqPMf2KMEZH6URERERakjFmtrW2ZE/jUvKInYh4b86CBVT95zpGV86kkixe7PIzBp72a37avcjraCIiIiKyAxU7EdnOe/O/ovylWzi1+ll8xqW020T6nHktp3Xq4nU0EREREdkNFTsRAWBR+UY+/PcdnLzpMQ43lSzrMo4uZ/2JkZ0P9DqaiIiIiOyBip1IO7e+MsyrU/7B4Sv+xrlmLeUFJUTPupm+PUd4HU1ERERE9pKKnUg7VRdNMGPaM/T/4jZ+YJayNrMvVaf+i25DTwGdFEVEREQkrajYibQzrmt5fdabZP/3j5zlfsKmQCfWHXMnnY/4ETg+r+OJiIiIyD5QsRNpR0q/+IIt06/j+Mgb1DpZrBx2OT1P/hUEMr2OJiIiIiL7QcVOpB2Yu2QFK16YxAmVz4OBpf3P48AzryEnu9DraCIiIiLSDFTsRNqwBSvXMWfq7Zy48XEGm1oWdz2VXufcSP+i3l5HExEREZFmpGIn0gYtXr2F96b+nePXPMh3zQZWFB5O4Iw/MaDXYV5HExEREZEWoGIn0oYsX1/NzBef4OgV9/AjZyVrcgZSPf5+eg063utoIiIiItKCVOxE2oBVm2p57qXpjFj8F37izGNzZneqx95Pl2HfAcfxOp6IiIiItDAVO5E0tqYizL9efouD5v2FXzrvUxPsQNVRN1JwxMXgD3odT0RERERaiYqdSBpaUxHmkVc/puvn9/Az51Xw+aka+Wtyv30ZhPK8jiciIiIirSwli50xZgIwoV+/fl5HEUkpayrC/PONuWR+cj8/dV4k0xelbuhEck68imBeV6/jiYiIiIhHUrLYWWunAdNKSkou8jqLSCpYWxnmvjcXEit9lJ85U+js20Jt33H4Tv4jOcUHeR1PRERERDyWksVORJLWVob5+5tLWFf6HJc5T9LPV06k6wg4eRJZPUd7HU9EREREUoSKnUgKWlsZ5u9vLWXBR6/yW+dflPgWEutwIJz0BBkDTwVjvI4oIiIiIilExU4khayrDPP3WUv54MP3+bXzJNf5S0lkdYLj/kLgsHPBpykrIiIiIjvTq0SRFBBLuPz19cU8//ZsfmqmcJX/LUwgE468Ct+YSyGY7XVEEREREUlhKnYiHlu1qZbf/+tdRq95nNcDLxM0CUzJj+GY30F2kdfxRERERCQNqNiJeOjlz1cy+/k7udtOoaO/EgafBcdfDYV9vY4mIiIiImlExU7EA+FonOf/9TdGLbuXcc4awj0Oh1MmQffhXkcTERERkTSkYifSyso+e42qaVcwMbGQ9dl9iZ3+NKEBJ+lMlyIiIiKyz1TsRFqJXfcl5c/+gR5r32AdhSwY+ScGjbsEHJ/X0UREREQkzanYibS0qjXEXv8Tvs8eI89m8FTe+Xz7vGsY1LHQ62QiIiIi0kao2Im0lEgVvHc37rt/xcSjPBIfS/yI33DBSSPwOXrbpYiIiIg0HxU7keaWiMEnj2DfuhlTs54Z7mgmh87ltz86hZF9dJRORERERJqfip1Ic7EWFkyD16+HjUtYlHEIv4v8nE6DjuCBcw6hQ1bQ64QiIiIi0kap2Ik0h5UfwCtXQ9lH1Ob352r/FUyrPpgrxg/iR4f3xuiMlyIiIiLSglTsRPbH+kXJI3RfTsfN6cxrB17JT+cPpEfHPJ47/zCGds/3OqGIiIiItAMqdiL7omotzLoZZj9Cwp/J610u4g/lR7JxQ4CzDuvODWcMJSdD00tEREREWodeeYo0RaQa3r8H++5fsfEIM0KncvXmU6ip68Bph3bjf0b34tADOnidUkRERETaGRU7kb2RiMEnj5J48yZ8tet5zYxhUvg7mKy+XHpqL84Z3kMnRxERERERz6jYiXwTa0ksmE54xjVkVy1jtjuQm+I/p3jgEUwa04sjDizC0XfSiYiIiIjHVOxEdqNy0TvUTr+CLpWfU+525+/+y+kx5kzuHdWLbh0yvY4nIiIiItJAxU6kEWst8+Z8QuLV6zi06m3CtgP3dfglPb79Y24e2oOg3/E6ooiIiIjITlqt2Blj+gJXAvn7FRJdAAAgAElEQVTW2nNaa78ie6M2Gmfmh18QeOc2xoVfJkKQN7r9mANO/R0/6dHZ63giIiIiIt9or4qdMeYhYDywzlo7tNHyccBdgA940Fp78+62Ya1dBlxojJmyf5FFms/mmij3v/YFOZ/cx494kZCJsbTXd+lxxnUcV9jV63giIiIiIntlb4/YTQbuAR7dusAY4wPuBcYCZcDHxpgXSZa8m3ZY/wJr7br9TivSTGIJl8ff/4qvXnuQn9l/0clsYVOvcWRPmMSAov5exxMRERERaZK9KnbW2reNMb13WDwSWFJ/JA5jzFPA6dbam0ge3RNJSW8tXMeLL/ybC6rv53xnOXVdhsH4Wyg8YKTX0URERERE9sn+fMauO7Cq0f0yYNTuBhtjOgI3AocZY/5QXwB3Ne5i4GKAnj177kc8ke0tWVfNfVNf47hV9/Bn38fU5XTFnvxPMoeeDUZfWSAiIiIi6Wt/it2uXgnb3Q221m4ELtnTRq219wP3A5SUlOx2eyJ7q6I2xt9f+YSOs//Kjc7LmGCA+FFXknnEzyGgry0QERERkfS3P8WuDDig0f0eQPn+xRFpPvGEy5MfLGPFa//gEvdpCn1VRId+n9BJ10FuF6/jiYiIiIg0m/0pdh8D/Y0xfYCvge8DP2iWVCL76e1F6/nP1Cc4v/oBznXKqOk2GmfCLYS6fcvraCIiIiIizW5vv+7gSeBYoMgYUwZca639pzHmZ8BMkmfCfMhaO685QhljJgAT+vXr1xybk3Zk2fpq/vn8TI5bdTc3+z6lNvcA7KmPkj3oNH2OTkRERETaLGNt6n6MraSkxJaWlnodQ9JARV2MB2eWUjz7TiY6r+H6M/Ed+1v8Y/4P/BlexxMRERER2SfGmNnW2pI9jduft2KKeC6ecHn6w2V8/eo9XOw+Q56vjvCh55I19mrIKfY6noiIiIhIq1Cxk7T17uL1zJz6CD+qepADndVU9zgK57Rbyeo82OtoIiIiIiKtSsVO0s7yDTU8PPU/HL/iLm7wzaE6rw92wtPkHHSSPkcnIiIiIu1SShY7nTxFdqUyHOPhmR/RqfQOrnHeIJaRS+y4m8kZ9WPwBbyOJyIiIiLimZQsdtbaacC0kpKSi7zOIt5LuJYpHy5lzat/4fzEFLJ9USKH/ZissVdAVqHX8UREREREPJeSxU5kq/cWr+f1qQ/xv1UP0stZR2XP4/GdfgtZRf29jiYiIiIikjJU7CQlLV1fzT9feI1xK27nat8cKjv0x55+P3kHftvraCIiIiIiKUfFTlLK5poo9746l7zZ93Ct8yIEg8SOv4W8UT8Gn/66ioiIiIjsSkq+UtbJU9qfaNzl0feX8/HrU/iD+yC9fWsJDzyL0Kk3Q25nr+OJiIiIiKS0lCx2OnlK+2GtZea8tTzw0rucX30/9/k+IFrQF05/gVDfY72OJyIiIiKSFlKy2En7MKesghunz2HQqqd4NDCFzEACjrmS4BG/BH+G1/FERERERNKGip20utUVddw2cyHLPp3FzRkPMTCwHPfA43FOvR0K+3odT0REREQk7ajYSaupicS57+1lPPX25/yap7gj43Vsbhc4+VGcQaeBMV5HFBERERFJSyp20uISruXZT8q4/eUvObL2NV7PfIoctwoz6lLMt/8AGbleRxQRERERSWspWex0Vsy2472lG5g0fQHRNfN5KOcxhgbnQNeRMP7P0OVgr+OJiIiIiLQJKVnsdFbM9LdsfTV/+s+XvLNgJVdmT+OHoRcxvhwYdxcc9r/gOF5HFBERERFpM1Ky2En62lwT5a7XF/P4Bys4KfApH+U/Rl5kNXzrhzD2Bsgu8jqiiIiIiEibo2InzSIad3nsgxX89fXF5IZX82Lxvxlc+V/IGwTjH4Jeh3sdUURERESkzVKxk/1ireWV+Wu56T8LKNtYyaTOb/Nd3xM4dQZOuB7G/BR8Aa9jioiIiIi0aSp2ss/mfl3BH6fP58OvNnFm4Qpe6vwQ2RWLYcCpcPLN0KGn1xFFRERERNoFFTtpsjUVYW6buZDnPi2jb2Ydrx/4Igd+/QLk94SJT8GAk72OKCIiIiLSrqRksdPXHaSm2mic+2Yt4763l2Jdl78NnMO41X/HrK6GI38NR/8WgtlexxQRERERaXdSstjp6w5Si1v/BeO3zVzIuqoIPxlQw2WR+8j4qhR6HQGn/hk6DfQ6poiIiIhIu5WSxU5Sx3tLN3DjSwuYV17J6O5Bph30Op3nT4bMAjjjH3Do98EYr2OKiIiIiLRrKnaySys21jDppQW8On8t3fND/PvodZQsuBUzbzUMPw+OvwayCr2OKSIiIiIiqNjJDsKxBPfNWsa9by0h4Bj+eHQWP9hwN76PXocuB8N3H4UDRngdU0REREREGlGxkwb/Xbyea16Yx1cbajhjaEcmdXqDnI/uAicA426GEReBT39lRERERERSjV6lC2srw0x6aQHTPi+nd8csppzTkZLS38CSOTDkTDjpT5DXzeuYIiIiIiKyGyp27VjCtTz6/nLueGUR0YTLr47vx0/z3yMw8wrwh+D7T8LAU7yOKSIiIiIie6Bi1059tmoLVz4/h3nllRzVv4gbT+pBz/f+AO++AH2OgTPvg7yuXscUEREREZG9kJLFTl9Q3nIqamPcOvNL/vXRSopzMrjnB4dxat5XmGdOhOo1cML1cPgvwHG8jioiIiIiInspJYudvqC8+Vlref7Tr/nTfxawqSbKeYf35rLj+5L74Z3w/G3QoRdc+Ap0H+51VBERERERaaKULHbSvJasq+KqqXP5YNkmvnVAByafP5Kh2RXw5Bmw6gM4dCKcchtk5HodVURERERE9oGKXRtWF01w9xuLeeC/y8gM+PjTmQfz/REH4Mx/Hh79FVgXznoQDvmO11FFRERERGQ/qNi1Ua8vWMu1L86jbHMdZw/rwR9OGUhRIAbTfgafPg7dS+DsB6Gwj9dRRURERERkP6nYtTFfb6nj+hfn8cr8tfTvlMPTF49mVN+OUP4ZPHshbFwKR/0Gjr0cfAGv44qIiIiISDNQsWsjYgmXh975ir+8thiL5ffjBnLhkX0IOsB7d8Nr10N2MfxoGvQ5yuu4IiIiIiLSjFTs2oCPl2/iqufnsnBtFScM6sx1pw2mR0EWVK2FqZfA0jdg4Hg47W7IKvQ6roiIiIiINDMVuzS2sTrCzTO+5N+zy+jeIZMH/reEsYM7Jx9c/CpM/T+IVMGpf4aSC8AYbwOLiIiIiEiLULFLQ65reaZ0FTe//CXV4TiXHHMgvzi+H1lBP8Qj8Np18MHfoNOQ5FsvOw3yOrKIiIiIiLQgFbs0M7+8kqumzuGTlVsY2aeQSWcM5aDO9d8/t34hTLkQ1s6BkT+BsTdAIORtYBERERERaXEpWeyMMROACf369fM6SsqojsS589VFTH5vOfmZAW7/zqGcPaw7xhiwFj55BGZcDsEsmPg0DBjndWQREREREWklKVnsrLXTgGklJSUXeZ3Fa9ZaZsxdww3T5rOmMszEkT35/bgBdMgKJgfUbYZpv4T5L0DfY+GMf0BeVy8ji4iIiIhIK0vJYidJKzbWcM0L85i1aD2Duubxt/8ZxrCeBY0GvAfPXgTVa+CE6+HwX4DjeBdYREREREQ8oWKXgiLxBPfNWsa9by7B7xiuHj+YH43phd9XX9oScXj7Vnj7NijoDRe+At2He5pZRERERES8o2KXYt5dsoGrp85l2YYaTj24K1ePH0yX/EYnQNm8Ap67CFZ9CIf+AE65FTJyvQssIiIiIiKeU7FLEeuqwkyavoAXPy+nV8csHrlgJMccVLz9oLnPwrRfAxbOehAO+Y4nWUVEREREJLWo2Hks4Voe/2AFt89cSCTu8ovj+3PpsQcSCvi2DYpUw4zfw2ePQ48RcPaDybdgioiIiIiIoGLnqS/KtnDl83OZ83UFR/Yr4obTh9C3OGf7QeWfwbMXwsalcNRv4NjLwRfwJrCIiIiIiKQkFTsPVNTFuH3mQh7/cAXFORncPfEwxh/SNfmddFu5LnxwL7x2PWQXw4+mQZ+jvAstIiIiIiIpS8WuFVlreeGzcia9tIBNNRF+NKY3l514EHmhHY7AVa2FqZfA0jdg4Hg47W7IKvQmtIiIiIiIpDwVu1ayZF0117wwl/eWbuTQAzow+fwRDO2ev/PARa/A1P+DaDWMvxOGnw+Nj+SJiIiIiIjsQMWuhYVjCe55Ywn3vb2UzICPSWcMZeLInvicHcpaPAKvXgsf/h06DYFzXoJOA70JLSIiIiIiaUXFrgW9+eU6rnlxLqs21XHWYd35wymDKM7N2Hng+oUw5UJYOwdG/gTG3gCB0M7jREREREREdkHFrgWUb6njhmnzeXneGg4szubJi0Yz5sCOOw+0Fj55BGZcDsEsmPg0DBjX+oFFRERERCStqdg1o1jCZfK7y7nztUW41vLbkwZw0VF9CfqdnQfXbYYXfwELXoS+x8KZ90Ful9aOLCIiIiIibYCKXTMpXb6Jq6bO5cs1VRw/sBPXnTaEAwqzdj14xXvw7EVQvSb5tssxPwdnF+VPRERERERkL6RksTPGTAAm9OvXz+soe7SpJsotM77k6dJVdMsPcd+5wzlxcOftv5Nuq0QcZt0C/70dCnrDha9C92GtnllERERERNqWlCx21tppwLSSkpKLvM6yO65r+ffsVdw840uqwnF+cnRffnF8f7IzdvOUbl4Bz10Eqz6EQ38Ap9wKGbmtG1pERERERNqklCx2qe7LNZVc9fxcSldsZkTvAiadcTADunxDSZv7LEz7NWDh7H/Cwee0WlYREREREWn7VOyaoDYa585XF/HQu8vJC/m59ZxDOGdYD5wdv5Nuq0g1zPg9fPY49BgBZz+YfAumiIiIiIhIM1Kxa4K4a3nx83K+M7wHvx83kILs4O4Hl38GUy6ATcvgqN/AsZeDL9B6YUVEREREpN1QsWuCvFCAVy87hrzQNxQ014UP7oXXroecTnDedOh9ZOuFFBERERGRdkfFrom+sdRVrYWpl8DSN2DgeDjtbsgqbL1wIiIiIiLSLqnYNZdFr8DU/4NoDYy/E4afD7v6ygMREREREZFmpmK3v+IRePVa+PDv0GkInPMQdBrodSoREREREWlHVOz2x/qFMOVCWDsHRl0CJ1wPgZDXqUREREREpJ1RsdsX1sInj8CMyyGYBROfhgHjvE4lIiIiIiLtlIpdU9Vugmm/hAUvQt9j4cz7ILeL16lERERERKQdU7FrisrV8ODxUL0Wxt4AY34OjuN1KhERERERaedU7JoitwsMmgCHfA+6D/M6jYiIiIiICKBi1zTGwMm3eJ1CRERERERkO3ofoYiIiIiISJpTsRMREREREUlzKnYiIiIiIiJpTsVOREREREQkzanYiYiIiIiIpDkVOxERERERkTSnYiciIiIiIpLmjLXW6wy7ZYxZD6wA8oGKZtrs/mxrX9Zt6jp7O35P44qADU3Ybzprzr8f+6Olc6TKPNjX9ZuyTnPNA2g/cyFV5gG0n7mgeZCaUmUutJd5sK/ray60rFSZB9B+5kJLzoNe1triPY6y1qb8Bbg/Fba1L+s2dZ29Hb+ncUCp139urXVpzr8fqZwjVebBvq7flHWaax7Uj2kXcyFV5kFrZEmVuaB5kJqXVJkL7WUe7Ov6mgste0mVedAaWVJlLrRGT9jTJV3eijktRba1L+s2dZ29Hd+cz0m6S5XnoqVzpMo82Nf1m7KO5kHTpdJz0V7mguZBakqV56O9zIN9XV9zoWWl0nPRXuZCa/SEb5TSb8WUfWeMKbXWlnidQ8RrmgsimgciW2kuSFuWLkfspOnu9zqASIrQXBDRPBDZSnNB2iwdsRMREREREUlzOmInIiIiIiKS5lTsRERERERE0pyKnYiIiIiISJpTsWuHjDFnGGMeMMa8YIw50es8Il4xxvQ1xvzTGDPF6ywirckYk22MeaT+d8EPvc4j4gX9DpC2RsUuzRhjHjLGrDPGzN1h+ThjzEJjzBJjzOXftA1r7VRr7UXAecD3WjCuSItpprmwzFp7YcsmFWkdTZwTZwFT6n8XnNbqYUVaSFPmgX4HSFujYpd+JgPjGi8wxviAe4GTgcHARGPMYGPMwcaY6TtcOjVa9ar69UTS0WSaby6ItAWT2cs5AfQAVtUPS7RiRpGWNpm9nwcibYrf6wDSNNbat40xvXdYPBJYYq1dBmCMeQo43Vp7EzB+x20YYwxwMzDDWvtJyyYWaRnNMRdE2pKmzAmgjGS5+wz9J6+0IU2cB/NbN51Iy9I/5m1Dd7b9zyskf2F3/4bxPwdOAM4xxlzSksFEWlmT5oIxpqMx5h/AYcaYP7R0OBEP7G5OPAecbYz5OzDNi2AirWiX80C/A6St0RG7tsHsYtluv3neWvtX4K8tF0fEM02dCxsB/eeGtGW7nBPW2hrg/NYOI+KR3c0D/Q6QNkVH7NqGMuCARvd7AOUeZRHxkuaCyPY0J0Q0D6SdULFrGz4G+htj+hhjgsD3gRc9ziTiBc0Fke1pTohoHkg7oWKXZowxTwLvAwOMMWXGmAuttXHgZ8BMYAHwjLV2npc5RVqa5oLI9jQnRDQPpH0z1u724yciIiIiIiKSBnTETkREREREJM2p2ImIiIiIiKQ5FTsREREREZE0p2InIiIiIiKS5lTsRERERERE0pyKnYiIiIiISJpTsRMREREREUlzKnYiIiIiIiJpTsVOREREREQkzanYiYiIiIiIpDkVOxERERERkTSnYiciIiIiIpLmVOxERERERETSnIqdiIiIiIhImlOxExERERERSXMqdiIiIiIiImlOxU5ERERERCTNqdiJiIiIiIikORU7ERERERGRNKdiJyIiIiIikub8Xgf4JkVFRbZ3795exxAREREREfHE7NmzN1hri/c0LqWLXe/evSktLfU6hoiIiIiIiCeMMSv2ZpzeiikiIiIiIpLmVOxERERERETSnIqdiIiIiIhImlOxExERERERSXMqdiIiIiIiImlOxU5ERERERCTNqdg10QNvL+OrDTVexxAREREREWmgYtcE6yrD3PX6Yk68cxaTps+noi7mdSQREREREREVu6bolBfijd8cw9nDevDPd7/i2Nve5NH3lxNPuF5HExERERGRdkzFrok65Ya4+exDmP7zIxnYJY9rXpjHuLv+y5sL13kdTURERERE2ikVu300pFs+/7poFPefO5x4wuX8hz/mfx/6iEVrq7yOJiIiIiIi7UyrFjtjzDhjzEJjzBJjzOWtue+WYIzhxCFdeOXXx3DVqYP4bOVmTr7rv1w9dS6baqJexxMRERERkXbCWGtbZ0fG+IBFwFigDPgYmGitnb+7dUpKSmxpaWmr5GsOm2qi3PXaIh7/cCVZQR+/OK4/Pzq8N0G/DoyKiIiIiEjTGWNmW2tL9jSuNRvHSGCJtXaZtTYKPAWc3or7b3GF2UGuP30oL//yKIb3KuDG/yxg7J2zeHnuGlqrQIuIiIiIyN6x1hJPuETjLgk3vV+v+1txX92BVY3ulwGjWnH/raZ/51wmnz+SWYvWM2n6fC55fDaj+xZy1amDGdo93+t4IiIi0s5Ya0m4loS1uC64duvt5HLX1i+rv28tJOzWx2zDY66bXO42Wjdh68c33LYkXHZet355w7r1167dNjaZNTnGUn9d/1rbdZPLGj+OTa5vsfXLk7ex29bdaTsNy+vXTW4I1/2G7eywX7t1Gzs8bhs/3vDc1183+rNgp8fs9vcbLd9xfXaxzo7bbrzvbevt/NhO+9iXrI3G7pxx58f29DPv7ufYKSs7rvfNz0fy72/93ztLw+3G+73prIOZOLIn6ao1i53ZxbKdarEx5mLgYoCePdP3iQU45qBijvjlUTz58SrufHURE+55h+8M78FvThxAp7yQ1/FERETaLWstsYQlmnCJxV1iCZdowiWesA23Y/W3Y/Ed7tePi7sucdc2rJNwLXF32+1YInkkIO7Wj03Y+vEuMdeSqN9GLGHrx2/bxrZCRKMXo40LETsUpu1frO5YztoqY5IvMI0xOAYMyQVbbxsDjjHJF6Fbb9evs/U21K9bv07ydvJlq+PsvB1T//j222i07Ub7No2Dsu3+1nW2Zt/1Y9sGmYaxW+8nbzSMNdvGm4adbn3M7HK/Wx/blqNhjd2M3Xkfu8u647a3rb/rx2i0nV3tt/E+tj7HTc3qGIPPSV6c+r8vW29vXX5wmh+Aac1iVwYc0Oh+D6B8x0HW2vuB+yH5GbvWidZy/D6Hc0f34rRDu3Hvm0t4+N2vmP7Fai499kB+fFRfQgGf1xFFREQ8E0+4ROLJSziWqL+dIBxzidTf37Z82+3GYyOxbdfhhvu72F58+3Vb8lMSjgG/4+D3GfyOwe9z8DuGgM/B5xj8PkPASd4O+JKP+xxD0O+QWf8i02cMxhh8Do1ejJpGL0aTL1ad+rGOoeG2z2m0bv2YXa67dZlT/0K38facRi+GG7LsYt3GL5IbsiTX32ndHba77bq+gDQqULssTY1v79hgRNq51jx5ip/kyVOOB74mefKUH1hr5+1unXQ7ecreWL6hhptmLGDmvLV075DJ78YN4LRDu+kfJxERSWmua6mKxKmpv1RH4tREEvXXcWqiyWW1kQS10QS10Xj9dePbyeu6aKKhXMX382hS0O8Q8jtkBHxk+B1CO1xvf9tHKLBtbIbfIeh3CPiSl6DPIeA329/3OQR8hoA/ed/vSz6+tawF6kuYv36cz0kWNsfR73URaR57e/KUVit2AMaYU4C/AD7gIWvtjd80vi0Wu63eX7qRSS/NZ155JYf17MDV4wczrGeB17FERKQNi8ZdKsMxKuqSl8pG15XheHJ5bWz7MeHksqpIfK+OcDkGsoJ+MoM+soM+MoN+soK+hkt2/WM7FrCG8hVwCPl9ZAQaFTH/jmUt+XjQpwIlIm1fSha7pmrLxQ6SH+J89pMybpu5kPVVEU7/Vjd+N24g3Ttkeh1NRERSWDTusrk2yqaanS/bilrjApcsbXWxxDduNxRwyAsFyM9MXvIyG90O+cnLDJCT4Sc7w092RrKkZWf4/397dx5eV3XYe/+7zjnS0TxL1mR5nm0wthhsIAaCA2YIJZgQSmlIekPSTE/et2lvU1puLuW2Ic2bvG9zc5smucHklhCCQzCTGZwEwlTAI54nPEm2NVrzeM5Z7x/raLRsy0Y6k36f51nP3mfvtc9eR9a29NPae63+bRl+HylJHt2FIiIyhhTs4kh7d4B/e+0gP33jQwDu/9h0vrRiBun+SD4CKSIi0WStpa07QG1rN7Ut3dS2dg0sB29r7aa1K3DG98lM8Q0LZ75B4SyJ7LSB9axhdfw+PfctIhJrFOziUHVTJ4+s38Oz245TlOnnmzfMYfWSct1mIiIS50IhS31bN8ebuzje1Mnxpk5ONHdxsqWLupZuasIhbqQeNb/Pw6SsFIoy/RRl+SnKTCE/PZnc9OTTljmpSfi8kZyiVkRExpuCXRzbfPQUDz23i63HmlhYlsU/3Dyfy6fnR7tZIiJyFu3dAY6d6uBoQwdHG4eWqsZOeoKhIfVTk7wUZ/cFtvAy0z8kxBVmppCV4tOtjSIiE5iCXZyz1vLstuM8sn4Px5u7uHFBMd+6aS5T8tOj3TQRkQnJWktdWzeH6to52tjBsXBoOxJer2/rGVI/0++jIj+Nirw0KvLTKM9JpSQ7ldKcVEpzUshOTVJgExGRc1KwSxBdvUF++scP+bfXDxIIWu67cipfvW4mWSlJ0W6aiEhCCoUsx5s72V/bxsHaNvbXtHGgro39Na20DHq2zWOgNCfVBbe8NCaHl1PCYU7BTURExoKCXYKpaeniey/vZe3mKnLTkvm/Vs7m7ksn61kKEZGPoLmjlx3Hm/mgqpl9Na0cqG3jQG3bkGfdCjKSmVGYwcwiV6YXZjAlL42y3FSS9H+wiIiMMwW7BLWjupl/fH4X7x5qZPakDP7+5vl8bHZhtJslIhLzmjt72VndzPbqZj6obmZHdTNHGjr695dkp/SHt1lFmf3reenJUWy1iIhMdAp2Ccxay8s7T/JPL+7haGMH184p5IGb5zGzKDPaTRMRiQmtXb3sqG5hx6AQd6i+vX9/WU4qF5Vns7As2y1Ls8lVgBMRkRikYDcBdAeCPPb2YX74uwN09Ab5s8sr+Mb1s/XLiYhMKNZaDta18eb+erYea+KD6mY+rBsIcaXZKSwqz2ZRWTaLynNYVJatXjgREYkbCnYTSENbNz/YsI9fvnuUDL+Pr398Fn++bCrJPj37ISKJJxSy7Dzewuv7annv8Cm2HWuiubMXgOKswSHOLQsy/FFusYiIyIVTsJuA9p5s5eEXdvHG/nqmFaTzdzfN4/p5RRqVTUTiWihk2VvTynuHGnn3UAPvfthIQ7ubWmBucSaXVORwcXkOV84sYHJeWpRbKyIiMrYU7CYoay2v7a3j4Rd2cbCuneUz8vn7m+czvzQr2k0TERm12pYuNh89xev76nlpxwlOdbgeubKcVC6flsfVswu4elaheuNERCThKdhNcL3BEL989yg/2LCP5s5e7qqczF99Yg6FmfolSERiT0dPgO1Vzbx1sIHf76lhR3ULAGnJXlbOn8SK2YVcNi2P8lz1yImIyMSiYCeAm6PpX3+/n8fePkxKkpcvXzuDz185jZQkb7SbJiIT3Kn2HnadaOH5D47zm83V9ARCeAwsqcjl2rlFLJ+Rz/zSLPw+/X8lIiITl4KdDPFhXRv/9OIeNuyuoTw3lb9dNZebF5Xo+TsRiaj9Na28squGDbtr2HqsCWvB7/Nwx9Jyrp9XxCWTczWyr4iIyCAKdjKitw/U89Dzu9hzspXKKbn8wy3zuXhyTrSbJSIJrLmjl2e3VfPkxmP9t1heVJ7Nx+dOYumUXBaVZZOdlhTlVoqIiMQmBTs5o/joBFEAACAASURBVGDI8tTGY3zvlX3Ut3Vz+yVl/M2NcyjJTo1200QkAQRDlh3Vzbx1sJ63DzTw3uFGegIh5pdk8enKclYtKmFSVkq0mykiIhIXFOzknNq6A/yvPxzgZ28ewmPgix+bwRdXTCct2RftpolIHNp9ooXfbqnmmS3V1LZ2AzBnUiZXzSrg9kvKWFiWHeUWioiIxB8FOxm1Y40dfOelPbzwwQkmZfn5mxvmcvslZXg8ev5ORM6utqWLdVuP8/SWanafaMHnMVwzp4hbLy5h+YwCjcQrIiLyEcVUsDPG3Al8G5gHXGatHVVaU7CLrI2HG/nH53exraqZi8qz+daqeVwxPU8DrIjIEB09AV7ZWcPTW6p5c38dIQsXT87hjiVl3HJRKXka/ERERGTMxFqwmweEgH8HvqlgF7tCIcu6bdU8sn4vJ1u6mFmUwZ1Ly7l9SRlFmXomRmSiCgRDvHeokd9srualHSdo7wlSlpPK7ZeUcfuSMmYUZkS7iSIiIgkppoJd/8mMeQ0Fu7jQ2RPk2W3VPLWxio1HTuH1GK6ZXcidleVcN3cSyT5PtJsoIuMsFLK8suskv3r/GBsPn6KtO0Cm38dNi0q4fUkZl03N0y3bIiIi4yxug50x5n7gfoCKioqlR44ciVDr5Ew+rGtj7aYqfrO5ipqWbvLSk7ltcSl3Lp3M/NKsaDdPRMZYbUsX63ec5JfvHmVvTSvlualcM6eQ5TMKuG5uESlJmjBcREQkUiIe7IwxG4DiEXY9YK1dF67zGuqxi1vBkOWN/XU8tamKV3fW0BMMsaA0izuXlnPb4jJNKiwSxwLBEL/fU8sv3zvK6/vqsBbmFmfyl9fM4OZFJfi86qUXERGJhrjtsRtMwS52NXX0sG7rcZ7a5CYcTvZ6uH5+EXcunczVswr0S6BIHLDWsr26md9uqea5bcepb+thUpafT1dO5pMXlzJrUma0mygiIjLhjTbYacIyuSA5acl8dvlUPrt8KrtPtPDUxiqe2VrNi9tPUpTp51NLyrmzslwDKojEoJPNXfxmcxVPb67iYF07yV4PH59XxO2XlHHd3CL9YUZERCQORWpUzNuBHwKFQBOw1Vp7w7mOU49dfOkJuFu51m46xh/21hEMWZZU5HBn5WRuuaiEzJSkaDdRZMLqCYT43e4afr3xGK/vc1MUXDY1j9uXlHHTwhKy03R9ioiIxKKYvBXzfCnYxa/a1i6e2eJG1dxf20ZKkodVC0u4c2k5V0zP10h6IhHQ2RPkzQP1/G53Da/sqqGxvYfirBRWLy1n9dJyphakR7uJIiIicg4KdhITrLVsq2rmqY3HeHbbcVq7ApTnpnLHEveL5eS8tGg3USTh7DzezBPvHWXdluO0hqcoWDGnkNVLy7l6ViFe/WFFREQkbijYSczp6g3y8s6TrN1UxZsH6rEWlk3P587KclYtLCE1WUOoi3wUO6qb+X9e2csf9tbh93m4eVEJn1pSzmXT8jT3pIiISJxSsJOYVt3UydObqnhqUxVHGzvI8Pu45aIS7qwsZ0lFLsaoR0HkXE42d/H7PbXsPdnC7hOtvHe4kezUJL64Yjr3XDZFz82JiIgkAAU7iQuhkOX9w408tamKF7efoKMnyPTCdFYvLeeOJeVMykqJdhNFYsrxpk7e2F/Hs9uO8/bBBqyFDL+PWZMyuHpWIX9x1TSyUxXoREREEoWCncSdtu4AL24/wVMbj/H+4VN4DHxsdiF3Lp3Mx+cVkZKkWzVl4rHWsqO6hd9sruK1vbUcbugAYEp+Gn+yuIxPLi5lekG6erlFREQSlIKdxLVD9e2s3XSMpzdXc6K5i7RkL9fOKWLVomKunVNEul9TMEriau8OsOVoE+8dbuSlHSfYV9NGss/D1TMLWD6zgOUz8plbnKkwJyIiMgEo2ElCCIYs//lhAy9uP8HLO2uob+vG7/NwzZxCVi0s4bp5RWRpfjxJAIFgiA27a3n83SO8fbCBYMjiMbB4cg53LC3nlkWlemZORERkAlKwk4QTDFk2Hm5k/Y6TrN9xgpqWbpK9Hq6eVcCqRSWsnDdJv/hK3LDWUt3UyeajTWw+cqr/e7okO4XbFpexbEY+SypyyNQfLkRERCY0BTtJaKGQZcuxU6zffpL1O05S3dSJz2NYPrOAmxYWs3L+JPIz/NFupsgQzZ29vHWgnj/uq+OP++o43twFQEqSh2XT8/nTy6dw7ZxCfF5NTSAiIiKOgp1MGNZaPqhq7u/JO9LQgcfAFdPzWbWohBsWTKIoU6NrSuSFQpbt1c28vq+O1/fVsfVYE8GQJdPv48qZBeFeuVzmlmSSpDAnIiIiI1CwkwnJWsuuEy28tOMkL2w/wYd17RgDl07JY9WiYm5cWExJdmq0mykJrLa1izf21fP6vjre2F/HqY5ejIFFZdmsmF3Ix2YXsnhyjoKciIiIjIqCnUx41lr217bx4vYTvLTjJHtOtgJwSUUONy0s4caFxUzOS4tyKyWeBYIhqps62XK0iXcPNfLeoQYO1rUDUJCRzMdmFbJiTiFXzSzQrcEiIiJyQRTsRIY5WNfGS+HbNXdUtwBwUXk2Ny4s5qKyHIqy/BRl+slOTdIw8jJEV2+QvSdb2XWihd0nWjhU387Rxg6qT3USCLn/QzNTfFw6NY9Lp+Zx9awC5pdk4fHo+0hEREQ+GgU7kbM42tDB+h0nWL/jJFuPNQ3Zl+zzUJTpD5eU/sA3sO6WeWnJ+sU9gVhrOdrYwbaqZqpOdVB1qpPqU50cO9XB4fp2wvmNDL+PGYXpTM5LoyJcFpVnM7c4C6++H0RERGSMKdiJjFJNSxeH69upae2mtqWLutZualu7qW3toqbFbWvpCpx2nM9jKMjwDwl7/QEw022flJVCfnqyRjmMMdZaGtt7qDrVybYqdxvl+4caqW3t7q+Tl55MWU4q5bmpzJqUyfySLOaXZFGem6pALyIiIhEz2mDni0RjRGLZpKwUJmWdfdTMrt5gOPB1UdvSTU1LVzj8uVJ1qoPNR0/R2N5z2rHGQH66vz/sFWW6wFeU6adwUCAszPTj93nH62NOSIFgiIN17XxQ1cT26mYO1bdT3dTJ8aZOunpD/fVKslO4Yno+l07LY2lFLlPy00j3679HERERiR/6zUVkFFKSvEzOSzvnYCs9gRD1beHA1xf+hoTALnYdb6G+rbv/1r7BctOS+nv/Cgf1/k3KGtojmJo8+gBoraU7EKKrN0gwZMlKTYr7ERmttbR0BTjV3sOpDlca23tp6uihsd2V/bVt7Dze3B/g0pO9zJyUybziLD4+t4jSnFRKc1L7e+H0XKWIiIjEMwU7kTGU7PP0B4azCYYsDW0DYa+2xa0P7gk8WNtGXVs3vcHTE2Cm30dhlp9J4ZDX1RukszdIV68LcH2lb9twGX4fOWlJrqQmn7ael54c7sl0oTIzJemCvybBkKWutbu/p+x4UycnW7roCYQIhmx/CYQsQWsJBsPL8LZQyBIIhQgELc2dveEg10twpGQMeD2G3LQkpuan86eXTeGi8mwWlWczLT9dt1CKiIhIwopIsDPG/AtwK9ADHAQ+Z61tOvtRIonL6zEUZaVQlJUCZJ+xXihkOdXRM9DjFw5+dYNCYEtLL6lJXtKTfeSne0hJ8pKS5CU1yUtKkofUJC/+8DavwfV0dfTQ3NFLU6fr5Tre1Nm/PlJeSk/29t+yOinLz6TsFLJSkugOhOgOBOnuHbx06y2dAY43d3Kyuat/5MjB75eS5MXrMfg8Bk946e0vHrwe8Ho8brsxJHk9zCzKICctmbz0JHLTkslNSyYvPZnc9GRy05LITU8m0+9T75uIiIhMOJHqsXsV+Ja1NmCMeQT4FvBfI3Rukbjl8RjyM/zkZ/iZVzL+5wuFLK3dARrbe6hp6eovJ5u7qWntoqa5i01HT1HT3E1P0PUE+n0e/D4XKP1JHvw+L36fh3S/G/6/NCfF9WJmu57MkhwXCkVERERk7EQk2FlrXxn08j+B1ZE4r4icH4/HkJ2aRHZqEtMK0s9Yz1pLTzBEstej3jERERGRGBCNERQ+D6yPwnlFZIwYY/D7vAp1IiIiIjFizHrsjDEbgOIRdj1grV0XrvMAEAAeP8v73A/cD1BRUTFWzRMREREREUlYYxbsrLXXn22/MeazwC3Ax+1ZZkW31v4E+Am4CcrHqn0iIiIiIiKJKlKjYt6IGyxlhbW2IxLnFBERERERmSgi9Yzd/wQygVeNMVuNMT+O0HlFREREREQSXqRGxZwZifOIiIiIiIhMRNEYFVNERERERETGkIKdiIiIiIhInFOwExERERERiXMKdiIiIiIiInFOwU5ERERERCTOKdiJiIiIiIjEOQU7ERERERGROKdgJyIiIiIiEucU7EREREREROKcgp2IiIiIiEicU7ATERERERGJcwp2IiIiIiIicU7BTkREREREJM4p2ImIiIiIiMQ5BTsREREREZE4p2AnIiIiIiIS5xTsRERERERE4pyCnYiIiIiISJxTsBMREREREYlzEQl2xph/NMZ8YIzZaox5xRhTGonzioiIiIiITAS+CJ3nX6y1/wBgjPk68CDwpQide+xYC7+4DXImQ+E8KJrrllmlYEy0WyciIiIiIhNURIKdtbZl0Mt0wEbivGOupw1sCPa+BFv+Y2C7PxsK54SDXrgUzYPMEgU+EREREREZd8bayGQsY8z/AP4caAautdbWnaHe/cD9ABUVFUuPHDkSkfadt/Z6qNsDtbvDyz1Qtxs6GgbqpGQPDXp9y4xJCnwiIiIiInJOxphN1trKc9Ybq2BnjNkAFI+w6wFr7bpB9b4FpFhr/9u53rOystJu3LhxTNoXMW11LuD1Bb26vS78dTYO1EnJCYe8uUNv6cwoUuATEREREZF+EQ92o2WMmQK8YK1deK66cRnsRmIttNcN6t0btOxqGqiXmjso6A3q6UsvVOATEREREZmARhvsIvKMnTFmlrV2f/jlJ4E9kThvzDDG9cZlFMH0FQPbrYW2mqG3ctbugR2/ga7mgXqpeUNv5ewPfAWR/ywiIiIiIhJzIjUq5neMMXOAEHCEeBwRczwYA5nFrky/ZmC7tdB68vRbOrevhe5BgS+tYNAtnX2hbx6k50f6k4iIiIiISBRFalTMOyJxnoRhDGSVuDLjuoHt1kLridNv6fzg19A9aODR9MKhvXt962l5kf8sIiIiIiIy7iLVYydjwRg3Z15WKcz8+MB2a6Hl+NAevto9sPUJ6GkdqJdedPqALUVz3bN9IiIiIiIStxTsEoExkF3myszrB7ZbC81Vg3r39rrQt/VxNydfn4zi8Dx8w57jS82J/GcREREREZHzpmCXyIyBnMmuzFo5sD0Ugpaqob17dbth8y+gt2OgXkYxFM6Ggjku+BXMdkvNwyciIiIiElMU7CYijwdyKlyZ/YmB7aEQNB8b2sNXvxe2/WroLZ3+7EGBb9AyZwp4vJH/PCIiIiIiE5yCnQzweCB3iiuzbxjY3jdoS91eqN83sNz/Cmz9j4F6Xj8UzBro2SuY5UJf/kxISon85xERERERmSAU7OTcBg/aMuPaofs6T0HdPtez1xf4qjfBzt8CNny8x/XmDb6ds6+XLyU74h9HRERERCTRKNjJR5OaCxWXuzJYbyc0HDi9l+/g7yHYM1DvtOf4wr18mcV6jk9EREREZJQU7GR8JKVC8SJXBgsGoOnIwPN7fb19Iz3HVzBroJcvfwbkzYDcqZCcFtGPIiIiIiIS6xTsJLK8PhfS8mcANw1sP+NzfK+66RkGyyyFvOmQP90t86a70Jc3DZLTI/pxRERERERigYKdxIazPsfXBI0fhsshaDzo1veuh/a6oXUzis8c+vyZkfs8IiIiIiIRpGAnsS81B8qWuDJcV3M47H04tOx/FdpqhtZNLwrf0jndBb3+0DcdUrIi81lERERERMaBgp3Et5RsKF3synDdrcNC30H3+uDvYeuJoXVT8wbm9sud4kbx7HudU6FbPEVEREQkpinYSeLyZ0LJRa4M19M+KPQdhKajcOqIm5h938sQ7B5aP61gUOjrC3x9AXCyGyxGRERERCRKFOxkYkpOh+KFrgwXCkF7rQt7TUfdKJ6njrj1Ex/AnheGTtkA7jbP3GG9fH3BL7tcE7SLiIiIyLhSsBMZzuNx8+hlFsPky07fHwpB28mB4HfqiAt/TUehejPsWgehwNBjMktGCH3h9ezJ4EuOzGcTERERkYSkYCdyvjyegRE8K644fX8o6KZuGBL8wj1/x96DHU+DDQ46IDwi6OCglzM5vKwI9/jpVk8REREROTMFO5Gx5vG6MJZdDlOWn74/GIDW48NCXzj4HXkHWtYOC35AeuEIgW/Q69ScyHw2EREREYlJCnYikeb1DfTOTb3q9P3BgOvxaz4GTceg+Wh4eQxqdrnBXQJdQ4/xZw0LfsMCYEaRmytQRERERBJSRIOdMeabwL8Ahdba+kieWyRueH0umOVMhikj7LcW2uuHBr7ByyPvQHfzsPf0ux7EM/X4ZZW584qIiIhIXIrYb3LGmMnASuBopM4pkpCMgYxCV8qWjlynq3lY6BsUAve97Eb9HPKeHsgsHRr8cqdA7lRXssrcLaYiIiIiEpMi+Sf6HwB/A6yL4DlFJqaUbCjOHnk6B4DeLmiuGrnX7+h/wo7fDH3Oz+NzgS936tDAlxNeT83VrZ4iIiIiURSRYGeM+SRQba3dZs7xy58x5n7gfoCKiooItE5kAkpKgYKZrowkGICWKjh12A3wcuqwK01HYPdz0NEwtL4/G3Irhga+vGmQP9MFQvX2iYiIiIyrMQt2xpgNQPEIux4A/g74xGjex1r7E+AnAJWVlXas2ici58HrGwhpI+luPT3wnToMdXth/6tDB3fxJkNuOOTlTw8vZ0LeDDdXoHr6RERERD6yMQt21trrR9pujFkETAP6euvKgc3GmMustSfH6vwiEkH+THeb50i3eoZC0FYDjR9C40FoOAgNB9zywAYIdg/UTUqH/BnhMnOgFMyGlKzIfR4RERGRODfut2Jaa7cDRX2vjTGHgUqNiimSoDweyCpxZeqVQ/eFgtBSPRD0Gg64cnwr7FoHNjRQN6sciuZC4VwomueWhXNcqBQRERGRITS+uYhEjsc7MIffjOuG7gv0uNs56/dB3R5XavfAoTeG9vJlTw6HvUGhr2g+JKVG9KOIiIiIxJKIBztr7dRIn1NE4oAvGQpnuzLvloHtoaALfLW7oW63e46vdg8c+uNA4DNed/tm8SIouQiKL3LraXlR+SgiIiIikaYeOxGJbR7vwHN4gwNfMBAOfLvg5HY4+QEcfhO2/3qgTnZFOOgtcmGv9BJ3i6iIiIhIglGwE5H45PUNTNkw/5MD29vr4cQ2F/RObocTH8CeF4DwILtZZW5i9/JKKKuE0sWQnB6VjyAiIiIyVhTsRCSxpBfAzI+70qe7DWp2wvHNULURqt6H3c+6fcYLk+a7kFdeCeWXuts6NQ2DiIiIxBEFOxFJfP4MqLjclT5tdVC9yYW86o2w4zew6VG3L60ApiyDKVfBlOUwaYEmWRcREZGYpmAnIhNTRiHMudEVcPPv1e+DY+/C0Xfg8Fuw+zm3LyUbKpa5kDflSihZ7G4FFREREYkR+s1ERATc/HtF4WkUln7WbWs6CkfegSNvwZG3Yd9Lbrs/C6Z9DKZf46ZtyJuuWzdFREQkqhTsRETOpG/OvYvvcq9ba1zIO/Q6HPw97Hl+oN6M62D6tS7waZoFERERiTBjrY12G86osrLSbty4MdrNEBE5nbXQ+KELeB++5ubV624B44Hyy9wtnrNvdJOoqzdPRERELpAxZpO1tvKc9RTsRETGQDDgBmM5sAH2v+ymXADImeIC3uwbYOpV4PNHt50iIiISVxTsRESiqeU47HvZlQ9fg0AnJGe4gDf/Npi5EpLTot1KERERiXEKdiIisaK3092quecF91xeRwMkpcGsT7iQN/sGTZIuIiIiI1KwExGJRcGAG4Bl1zo3nUJ7LfhSYdb1MP9PXMjzZ0a7lSIiIhIjFOxERGJdKOjmzNu1DnY9C20nweuHWSvhortcyNMzeSIiIhOagp2ISDwJhdzk6LuegZ2/hbYaNzH6gtvhos9AxRUaXVNERGQCUrATEYlXwQAceg22PemeyevtcKNrXnSXKwUzo91CERERiRAFOxGRRNDd5sLdtl+5idFtCMqWwuI/hYWrITUn2i0UERGRcaRgJyKSaFpOwI61LuTV7HCDrsy/DZbcC1Ou1K2aIiIiCUjBTkQkUVkLJ7bC5l/A9rXQ3QJ5M+CSP3M9eZnF0W6hiIiIjJGYCnbGmG8DXwDqwpv+zlr74rmOU7ATETmHng43quaW/+OmUTBeN5rmJfe6efK8vmi3UERERD6C0Qa7SP7E/4G19nsRPJ+ISOJLToPFd7tSf8AFvK2/hL0vQsYk14O35M8hb3q0WyoiIiLjyBPtBoiIyBgpmAkr/zv837vgM7+E0kvgrf8P/vUS+MVtbhqFQE+0WykiIiLjIJI9dl81xvw5sBH4K2vtqQieW0Rk4vAmwdybXWmuhq2Pu+fxnroP0gpcL97S+yB/RrRbKiIiImNkzJ6xM8ZsAEZ6Yv8B4D+BesAC/wiUWGs/f4b3uR+4H6CiomLpkSNHxqR9IiITWigIB/8Amx6FvevBBmHq1S7gzbsVfP5ot1BERERGEFODpww5oTFTgeettQvPVVeDp4iIjIPWk7DlP1wvXtMRSM0LP4v3WSicHe3WiYiIyCCjDXYRecbOGFMy6OXtwI5InFdEREaQWQwf+yZ8fSvc+1uYdjW8+2P40aXw81Ww7Uno7Yx2K0VEROQ8RGq6g/8DLMbdinkY+KK19sS5jlOPnYhIhLTVutE0Nz8GjR9CSg5c/BnXizdpfrRbJyIiMmHF7K2Y50PBTkQkwkIhOPyGC3i7n4NgD0y+3AW8Bbe76RVEREQkYhTsRETko2mvh21PwKbHoGE/+LPhok/D0s9C8aJot05ERGRCSNhg19vbS1VVFV1dXVFqVXxJSUmhvLycpKSkaDdFROKVtXDkbdi0Bnatg2A3lC11I2ou+BT4M6LdQhERkYSVsMHu0KFDZGZmkp+fjzEmSi2LD9ZaGhoaaG1tZdq0adFujogkgo5G+OBJF/Lq9kByJixa7UJe6eJot05ERCThjDbYRXKC8jHR1dXF1KlTFepGwRhDfn4+dXV10W6KiCSKtDy44i/h8i/BsfdcwNv2hJsfr2gBXHwXLFwN2WXRbqmIiMiEEpHpDsaaQt3o6WslIuPCGKi4HG7/N/irvXDT99zAKq8+CD9YAI99ErY8Dl0t0W6piIjIhBCXwU5ERGJIag5c9gX4Lxvga5thxX+FpqOw7svwvdmw9vOw72UI9ka7pSIiIglLwe4CLF++fNzPsWbNGo4fP35exxw+fJiFCxeOU4tEREYhfwZc+y34+hb4i1fhknvg4O/hl592Ie/Zr8GB3ynkiYiIjLG4e8YuFrz99tvj+v7BYJA1a9awcOFCSktLx/VcIiLjwhiYfJkrN/wzHHgVdv4WdjwNm38Bqbkw9xaY/ycwfQV4NXKviIjIR6EeuwuQkeGG9n7ttde45pprWL16NXPnzuWee+7BWsv69ev59Kc/3V//tdde49ZbbwXglVdeYdmyZSxZsoQ777yTtrY2AKZOncpDDz3EVVddxRNPPMHGjRu55557WLx4MZ2dnWzatIkVK1awdOlSbrjhBk6cOAHApk2buPjii1m2bBk/+tGPIvyVEBEZBV8yzL0Z7vgZ/PVB+MwvYeZK2PkMPH4H/MtMeOYrsOdF6OmIdmtFRETiUlz32P3353ay6/jYPpg/vzSL/3brglHX37JlCzt37qS0tJQrr7ySt956i5UrV/LFL36R9vZ20tPTefLJJ7nrrruor6/n4YcfZsOGDaSnp/PII4/w/e9/nwcffBBwc869+eabAPzsZz/je9/7HpWVlfT29vK1r32NdevWUVhYyJNPPskDDzzAz3/+cz73uc/xwx/+kBUrVvDXf/3XY/q1EBEZc0kpLuTNvRl6u9xtmruegd3Pwtb/AF8KTL8G5qyC2TdCZnG0WywiIhIX4jrYxYLLLruM8vJyABYvXszhw4e56qqruPHGG3nuuedYvXo1L7zwAt/97nd5/fXX2bVrF1deeSUAPT09LFu2rP+97rrrrhHPsXfvXnbs2MHKlSsBd6tmSUkJzc3NNDU1sWLFCgDuvfde1q9fP54fV0Rk7CSlwNybXAn0wJG3YN9LsPdFtwQovQRmr4JZ10PJYvB4o9tmERGRGBXXwe58etbGi9/v71/3er0EAgHAhbQf/ehH5OXlcemll5KZmYm1lpUrV/LEE0+M+F7p6ekjbrfWsmDBAt55550h25uamjSdgYgkBl8yzLjWlRu/A7W7BwLea/8Mr/2Tey5v2gqYcZ2rl1MR7VaLiIjEDD1jN06uueYaNm/ezE9/+tP+nrgrrriCt956iwMHDgDQ0dHBvn37Rjw+MzOT1tZWAObMmUNdXV1/sOvt7WXnzp3k5OSQnZ3df/vm448/Pt4fS0Rk/BkDk+bDx77pplD45n6443/DnJvcpOjPfR3+30Xww6XwwjdhzwvQ0RjtVouIiERVXPfYxTKv18stt9zCmjVreOyxxwAoLCxkzZo13H333XR3dwPw8MMPM3v27NOOv++++/jSl75Eamoq77zzDmvXruXrX/86zc3NBAIBvvGNb7BgwQIeffRRPv/5z5OWlsYNN9wQ0c8oIhIRGYWwaLUr1kLdXvjwD+75vK2Pw/s/dfWK5sOU5VCxDKZcCVkl0W23iIhIBBlrbbTbcEaVlZV248aNQ7bt3r2befPmRalF8UlfMxFJWIFuqN4ER9525di70ONGGyZ32kDQK78UCmbpGT0REYk7xphN1trKj9/YWgAAC39JREFUc9VTj52IiMQvn9+FtynL3etgAGq2DwS9vetdrx5AciaULoayJVC21JWsMnfrp4iISJxTsBMRkcTh9bmRNEsvgWVfgVAIGg64Xr2+8s7/glCvq58xyQW8koth0kIoXgg5UxT2REQk7ijYiYhI4vJ4oHC2K4vvdtsC3XByx9Cwt3c9EH40wZ8FkxYMBL1Ji6BoHiSnRe1jiIiInIuCnYiITCw+P5QvdaVPT7ubYuHkdqjZ4YLftl/B+63hCsZNr1A4Bwpmu+f1CsLr6flR+RgiIiKDRSzYGWO+BnwVCAAvWGv/JlLnFhEROavkdCivdKVPKARNR1zQq9kF9Xuhfh8cegMCnQP1UvPCgW8W5M2AvGmQO9UN3pKSFfGPIiIiE1NEgp0x5lrgNuAia223MaYoEucVERG5YB6PC2l502DerQPbQyFoPgb1+wfCXt2+8Hx6DUPfIzVvaNDLnepKTgVklYI3KYIfSEREElmkeuz+EviOtbYbwFpbG6HzxqT77ruPW265hdWrV0e7KSIicr48Hsid4sqs64fu62qGU4ddaTw0sF69CXY+AzY4qLJxg7dkl7nRObPLw8syyCp3y4xJmqJBRERGJVLBbjZwtTHmfwBdwDette+PVNEYcz9wP0BFRUWEmnfhrLVYa/F4PNFuioiIRFtKthths+Ti0/cFA66n79QhaK6C5mpoCS/r9sCB30Fv+9BjPD5IL4KMIsgsdsuMScNKeF9SamQ+o4iIxKQxC3bGmA1A8Qi7HgifJxe4ArgU+LUxZrodYXZ0a+1PgJ+Am6B8rNo3lg4fPsyqVau49tpreeedd/jGN77Bj3/8Y7q7u5kxYwaPPvooGRkZPPTQQzz33HN0dnayfPly/v3f/x2jIbRFRCYmr2/g1s6RWAtdTeHAV+1CYHM1tNVCW43bdnwLtNeBDZ1+vD/Lhby0AkjLh7S8YcvBJQ/82a73UUREEsKYBTtr7fVn2meM+Uvg6XCQe88YEwIKgLqPdNL1f+tGMBtLxYtg1XfOWW3v3r08+uijPPTQQ3zqU59iw4YNpKen88gjj/D973+fBx98kK9+9as8+OCDANx77708//zz3Hrrred4ZxERmZCMgdRcV4oXnrleKOie5Ws9ORD62moG1jsaoOmoC4Ed9RDsOcP5vAOhLyXH9TamhpdnLOH9/iwXVEVEJGZE6n/lZ4DrgNeMMbOBZKA+QuceF1OmTOGKK67g+eefZ9euXVx55ZUA9PT0sGzZMgD+8Ic/8N3vfpeOjg4aGxtZsGCBgp2IiHw0Hm/4lsxRjENmrZvKoaMhXBoHrfeVevdsYNtJd0toVzN0t4zcKzhYcmY45GW6UUX9GZCcMfA6OWNgW//6sLp9230pmhReROQjilSw+znwc2PMDqAH+OxIt2Get1H0rI2X9PR0wD1jt3LlSp544okh+7u6uvjyl7/Mxo0bmTx5Mt/+9rfp6uqKRlNFRGSiMsYFJ3+GG+xltEIh6GlzIa+r2d0i2r8+qHQ2uRDY0+7qt9aE11uhuw1CvaNspweS0txzgkmpbt2XMmzboH1D6gzblpQCvlTwJbv93mQ3d+HgdW+ygqSIJJyIBDtrbQ/wZ5E4V6RdccUVfOUrX+HAgQPMnDmTjo4OqqqqKCpyf0ktKCigra2NtWvXahRMERGJDx6Pm4MvJQuYfOHvE+hxga+7dSD8DVlvcyGwpx16u6C3A3o73TyBveHS1eRuO+3b19vp1kcbGs/E6x8Ieb6+db8LhH37+rcNq+tNdgPbeJOGrnuS3C2q3uTwetKwfUlD63nCx5+1XpJCqIiMim6Q/4gKCwtZs2YNd999N93d3QA8/PDDzJ49my984QssWrSIqVOncumll0a5pSIiIhHmSwZfnnuWb6wFA4MC4ODQFy7Bbgh0u2cMA33rg7d1ueAZ7HbLQNew9R7oaB+hbpc7d6j3zM8vjjXjccHPeN3SM/y11xUTXvbvG+XrIe/jGzjfiHU97rXxDCpm6GvP2fYP3zdo/2nHDT/+TPs97mtypn2YgXBszKDXZuC9h28bsvSc/fghdc72PoPaITIOzFjcETleKisr7caNG4ds2717N/PmzYtSi+KTvmYiIiLjwFo3mE2oF4K9EAq4ZbAnvC0wbF9PeL13aDg8Z72AmwMxFHC3yQ55HXRlyOuAe0ayfz04UO+M7zX4uMBZjg2Fn7+M3d8fY9+ZQt+5AibnF0L7zgWDtpkhizPvH+/XY3X+sXqP8OtlX4XZNxBrjDGbrLWV56qnHjsRERGRC2FM+NZL38SbR9DacAkNKsFhr4ftDw3fP7jOSPsG7R/x2L79Z9oXcuGVcFvPugydYR8DbTjX+5xWh1HUOVM7uIC29r1/37/RsJX+zpxIv+Y864/y9UdtQ9/38OD6oSDxTMFORERERM5P/22FmgtRJFbE5dUYy7ePxhp9rUREREREEl/cBbuUlBQaGhoUWEbBWktDQwMpKSnRboqIiIiIiIyjuLsVs7y8nKqqKurq6qLdlLiQkpJCeXl5tJshIiIiIiLjKO6CXVJSEtOmTYt2M0RERERERGJG3N2KKSIiIiIiIkMp2ImIiIiIiMQ5BTsREREREZE4Z2J5dEljTB1wBMgGmsfobT/Ke13Ised7zGjrn6teAVB/HueNZ2P5/fFRjHc7YuU6uNDjz+eYsboOYOJcC7FyHcDEuRZ0HcSmWLkWJsp1cKHH61oYX7FyHcDEuRbG8zqYYq0tPGcta23MF+AnsfBeF3Ls+R4z2vrnqgdsjPa/W6TKWH5/xHI7YuU6uNDjz+eYsboOwnUmxLUQK9dBJNoSK9eCroPYLLFyLUyU6+BCj9e1ML4lVq6DSLQlVq6FSOSEc5V4uRXzuRh5rws59nyPGW39sfyaxLtY+VqMdzti5Tq40OPP5xhdB+cvlr4WE+Va0HUQm2Ll6zFRroMLPV7XwviKpa/FRLkWIpETziqmb8WUC2eM2WitrYx2O0SiTdeCiK4DkT66FiSRxUuPnZy/n0S7ASIxQteCiK4DkT66FiRhqcdOREREREQkzqnHTkREREREJM4p2ImIiIiIiMQ5BTsREREREZE4p2A3ARlj/sQY81NjzDpjzCei3R6RaDHGTDfG/G9jzNpot0Ukkowx6caYx8I/C+6JdntEokE/AyTRKNjFGWPMz40xtcaYHcO232iM2WuMOWCM+duzvYe19hlr7ReA+4C7xrG5IuNmjK6FD621fzG+LRWJjPO8Jj4FrA3/LPhkxBsrMk7O5zrQzwBJNAp28WcNcOPgDcYYL/AjYBUwH7jbGDPfGLPIGPP8sFI06NC/Dx8nEo/WMHbXgkgiWMMorwmgHDgWrhaMYBtFxtsaRn8diCQUX7QbIOfHWvtHY8zUYZsvAw5Yaz8EMMb8CrjNWvvPwC3D38MYY4DvAOuttZvHt8Ui42MsrgWRRHI+1wRQhQt3W9EfeSWBnOd1sCuyrRMZX/rPPDGUMfCXV3A/sMvOUv9rwPXAamPMl8azYSIRdl7XgjEm3xjzY+ASY8y3xrtxIlFwpmviaeAOY8y/Ac9Fo2EiETTidaCfAZJo1GOXGMwI284487y19l+Bfx2/5ohEzfleCw2A/rghiWzEa8Ja2w58LtKNEYmSM10H+hkgCUU9domhCpg86HU5cDxKbRGJJl0LIkPpmhDRdSAThIJdYngfmGWMmWaMSQY+Azwb5TaJRIOuBZGhdE2I6DqQCULBLs4YY54A3gHmGGOqjDF/Ya0NAF8FXgZ2A7+21u6MZjtFxpuuBZGhdE2I6DqQic1Ye8bHT0RERERERCQOqMdOREREREQkzinYiYiIiIiIxDkFOxERERERkTinYCciIiIiIhLnFOxERERERETinIKdiIiIiIhInFOwExERERERiXMKdiIiIiIiInFOwU5ERERERCTO/f+AKL9uxt/RMwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 1080x576 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(15, 8))\n",
"plt.subplot(211)\n",
"plt.loglog(freq, amp, label=\"inverted\")\n",
"plt.loglog(freq, np.abs(desired_spec), label=\"real\")\n",
"plt.legend()\n",
"\n",
"plt.subplot(212)\n",
"plt.semilogx(freq, np.unwrap(phase), label=\"inverted\")\n",
"plt.semilogx(freq, np.unwrap(np.angle(desired_spec)), label=\"real\")\n",
"plt.legend()\n",
"\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment