{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from sklearn.decomposition import PCA\n", "import os, sys\n", "from astropy.io import fits\n", "from lightkurve.lightcurve import LightCurve as LC" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "directory = '/Users/AdinaFeinstein/Downloads/'\n", "files = os.listdir(directory)\n", "targets = [i for i in files if 'TIC' in i]" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "fluxes = []\n", "cdpps, mags = [] ,[]\n", "for t in targets:\n", " hdu = fits.open(directory+t)\n", " mag = hdu[0].header['TMAG']\n", " q = hdu[1].data['QUALITY'] == 0\n", " flux = hdu[1].data['CORR_FLUX']\n", " fluxes.append(np.array(flux)/np.nanmedian(flux))\n", " cdpp = LC(time=hdu[1].data['TIME'][q][0:600], flux=flux[q][0:600]).estimate_cdpp(transit_duration=2)\n", " mags.append(mag)\n", " cdpps.append(cdpp)\n", "fluxes = np.array(fluxes)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6gAAAJCCAYAAADN6ep4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHBhJREFUeJzt3WGI5ft91/HP7+zZqRhhHOZGtE0yN9OEYCglMEN68IHmSTGWbooVpWEfKHZcK+Sh2JaACuKjIn0gi7Ku4SKMKUEQcuuVtsiWIPRA9kiqCWnqcugxty3GTA8jojB7PD8f3E28udl7d+/O7P6/c+b1gjyYs+fu/Sz8We47/9//nNZ7DwAAAAxtNPQAAAAASAQqAAAARQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBLGQw9Ikpdeeqm//PLLQ88AAADgOZjNZt/uvb/3Se8rEagvv/xy7t+/P/QMAAAAnoPW2uJp3ueILwAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShg0UFtrN1prd05PT4ecAQAAQAGDBmrv/dXe+63t7e0hZwAAAFCAI74AAACUIFABAAAukdlimdv3HmS2WA495cKNhx4AAADA05ktlrl5d5qz1Tpb41GOjyY52NsZetaFcQcVAADgkpjOT3K2Wmfdk4erdabzk6EnXSiBCgAAcElM9nezNR7lWkuuj0eZ7O8OPelCOeILAABwSRzs7eT4aJLp/CST/d2NOt6bCFQAAIBL5WBvZ+PC9Dsc8QUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAACe0WyxzO17DzJbLIeeshF8DyoAAMAjs8Uy0/lJJvu7T/yu0dlimZt3pzlbrbM1HuX4aLKx30/6oghUAACAvPvgnM5PcrZaZ92Th6t1pvMTgXpOjvgCAADk8cH5Tib7u9kaj3KtJdfHo0z2d1/Q0s3lDioAAED+f3A+XK2fKjgP9nZyfDR56iPBPFnrvQ+9IYeHh/3+/ftDzwAAAK64d/MMKk+vtTbrvR8+6X3uoAIAADxysLcjTAfkGVQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAACgvNlimdv3HmS2WA49hedoPPQAAACAdzJbLHPz7jRnq3W2xqMcH01ysLcz9CyeA3dQAQCA0qbzk5yt1ln35OFqnen8ZOhJPCcCFQAAKG2yv5ut8SjXWnJ9PMpkf3foSTwnjvgCAAClHezt5Phokun8JJP9Xcd7N5hABQAAyjvY2xGmV4AjvgAAAJQgUAEAAChBoAIAAFCCQAUAALikZotlbt97kNliOfSUC+FDkgAAAC6h2WKZm3enOVutszUe5fhocuk/SModVAAAgEtoOj/J2WqddU8ertaZzk+GnnRuAhUAAOASmuzvZms8yrWWXB+PMtnfHXrSuTniCwAAcAkd7O3k+GiS6fwkk/3dS3+8NxGoAAAAl9bB3s5GhOl3OOILAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAADwXM0Wy9y+9yCzxXLoKRQ3HnoAAACwuWaLZW7eneZstc7WeJTjo0kO9na++2vT+Ukm+7vffY2rTaACAADPzXR+krPVOuuePFytM52f5GBv5x3DlavLEV8AAOC5mezvZms8yrWWXB+PMtnfTfL4cAV3UAEAgOfmYG8nx0eT7zvK+51wfbhaf0+4crW13vvQG3J4eNjv378/9AwAAOAF8gzq1dFam/XeD5/0PndQAQCAQRzs7QhTvodnUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABcKrPFMrfvPchssRx6ChdsPPQAAACApzVbLHPz7jRnq3W2xqMcH01ysLcz9CwuiDuoAADApTGdn+Rstc66Jw9X60znJ0NP4gIJVAAA4NKY7O9mazzKtZZcH48y2d8dehIXyBFfAADg0jjY28nx0STT+Ukm+7vv6njvbLF823/unX6NF0egAgAAl8rB3s67jsh3enbVc611PJcjvq2197TWZq21n3wevz8AAMC78U7PrnqutY6nCtTW2udaa99qrX31La9/srX2jdbag9baL7zpl34+yRcucigAAMCzeqdnVz3XWkfrvT/5Ta39+ST/K8m/6r3/yKPXriX53SQ/nuT1JF9O8ukkP5jkpSR/LMm3e++/+qTf//DwsN+/f/9Z/wwAAABP5BnU4bTWZr33wye976meQe29f6m19vJbXv54kge99/mjf+GvJPmpJH8iyXuSfDTJ/2mtvdZ7X7+L7QAAABfunZ5dfZbnWrl45/mQpB9K8s03/fx6kh/rvX8mSVprfyNv3EF9bJy21m4luZUkH/jAB84xAwAAgE1wng9Jao957bvnhXvvr7zT8d7e+53e+2Hv/fC9733vOWYAAACwCc4TqK8nef+bfn5fkj843xwAAACuqvME6peTfLi19sHW2laSn0nyxYuZBQAAwFXztF8z8/kkv5XkI62111trP9t7XyX5TJJfS/L1JF/ovX/t+U0FAABgkz3tp/h++m1efy3Jaxe6CAAAgCvpPEd8AQAA4MIIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAgA0xWyxz+96DzBbLoafAM3mqr5kBAABqmy2WuXl3mrPVOlvjUY6PJjnY2xl6Frwr7qACAMAGmM5PcrZaZ92Th6t1pvOToSfBuyZQAQBgA0z2d7M1HuVaS66PR5ns7w49Cd61QY/4ttZuJLnxoQ99aMgZAABw6R3s7eT4aJLp/CST/V3He7mUWu996A05PDzs9+/fH3oGAAAAz0FrbdZ7P3zS+xzxBQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAuEJmi2Vu33uQ2WI59BT4PoN+zQwAAPDizBbL3Lw7zdlqna3xKMdHE19HQynuoAIAwBUxnZ/kbLXOuicPV+tM5ydDT4LvIVABAOCKmOzvZms8yrWWXB+PMtnfHXoSfA9HfAEA4Io42NvJ8dEk0/lJJvu7jvdSjkAFAIAr5GBvR5hSliO+AAAAlCBQAQAAKEGgAgAAUIJABQAAoIRBA7W1dqO1duf09HTIGQAAABQwaKD23l/tvd/a3t4ecgYAAAAFOOILAABX3GyxzO17DzJbLIeewhXne1ABAOAxZotlpvOTTPZ3N/p7Q2eLZW7eneZstc7WeJTjo8lG/3mpTaACAMBbXKVom85PcrZaZ92Th6t1pvOTjf2zUp8jvgAA8BaPi7ZNNdnfzdZ4lGstuT4eZbK/O/QkrjB3UAEA4C2+E20PV+uNj7aDvZ0cH02uxHFm6mu996E35PDwsN+/f3/oGQAA8F1X5RlUeBFaa7Pe++GT3ucOKgAAPMbB3o4whRfMM6gAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKGHQQG2t3Wit3Tk9PR1yBgAAAAUMGqi991d777e2t7eHnAEAAEABjvgCAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUA4FKZLZa5fe9BZovl0FOACzYeegAAADyt2WKZm3enOVutszUe5fhokoO9naFnARfEHVQAAC6N6fwkZ6t11j15uFpnOj8ZehJwgQYN1NbajdbandPT0yFnAABwSUz2d7M1HuVaS66PR5ns7w49CbhArfc+9IYcHh72+/fvDz0DAIBLYLZYZjo/yWR/1/FeuCRaa7Pe++GT3ucZVAAALpWDvR1hChvKM6gAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlDBqorbUbrbU7p6enQ84AAACggEEDtff+au/91vb29pAzAAAAKMARXwAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKCEQQO1tXajtXbn9PR0yBkAAAAUMGig9t5f7b3f2t7eHnIGAAAABTjiCwAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQwqCB2lq70Vq7c3p6OuQMAAAAChg0UHvvr/beb21vbw85AwAuxGyxzO17DzJbLIeeAgCX0njoAQCwCWaLZW7eneZstc7WeJTjo0kO9naGngUAl4pnUAHgAkznJzlbrbPuycPVOtP5ydCTgAE4SQHn4w4qAFyAyf5utsajPFytc308ymR/d+hJwAvmJAWcn0AFgAtwsLeT46NJpvOTTPZ3/UcpXEGPO0nh7wJ4dwQqAFyQg70d/zEKV5iTFHB+AhUAAC6AkxRwfgIVAAAuiJMUcD4+xRcAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABgI0wWyxz+96DzBbLoacA8IzGQw8AADiv2WKZm3enOVutszUe5fhokoO9naFnAfAuuYMKAFx60/lJzlbrrHvycLXOdH4y9CQAnoFABQAuvcn+brbGo1xryfXxKJP93aEnAfAMHPEFAC69g72dHB9NMp2fZLK/63gvwCUlUAGAjXCwtyNMAS45R3wBAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoYNFBbazdaa3dOT0+HnAEAAEABgwZq7/3V3vut7e3tIWcAAABQgCO+AAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAANsBsscztew8yWyyHngLPbDz0AAAA4Hxmi2Vu3p3mbLXO1niU46NJDvZ2hp4F75o7qAAAcMlN5yc5W62z7snD1TrT+cnQk+CZCFQAALjkJvu72RqPcq0l18ejTPZ3h54Ez8QRXwAAuOQO9nZyfDTJdH6Syf6u471cWgIVAAA2wMHejjDl0nPEFwAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhVgILPFMrfvPchssRx6CgBACeOhBwBcRbPFMjfvTnO2WmdrPMrx0SQHeztDzwIAGJQ7qAADmM5PcrZaZ92Th6t1pvOToScBAAxOoAIMYLK/m63xKNdacn08ymR/d+hJAACDc8QXYAAHezs5PppkOj/JZH/X8V4AgAhUgMEc7O0IUwCAN3HEFwAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAOcwWyxz+96DzBbLoacAXHrjoQcAAFxWs8UyN+9Oc7ZaZ2s8yvHRJAd7O0PPAri03EEFAHhG0/lJzlbrrHvycLXOdH4y9CSAS02gAgA8o8n+brbGo1xryfXxKJP93aEnAVxqjvgCADyjg72dHB9NMp2fZLK/63gvwDkJVACAczjY2xGmABfEEV8AAABKEKgAAACUcOGB2lr7s621f95a+zettb9z0b8/AAAAm+mpArW19rnW2rdaa199y+ufbK19o7X2oLX2C0nSe/967/3nkvy1JIcXPxkAAIBN9LR3UF9J8sk3v9Bau5bkdpK/lOSjST7dWvvoo1/7VJL/mOQ/XNhSAAAANtpTBWrv/UtJ/ugtL388yYPe+7z3fpbkV5L81KP3f7H3/ueS3LzIsQAAAGyu83zNzA8l+eabfn49yY+11j6R5KeT/ECS197uH26t3UpyK0k+8IEPnGMGAAAAm+A8gdoe81rvvf9mkt980j/ce7+T5E6SHB4e9nPsAAAAYAOc51N8X0/y/jf9/L4kf3C+OQAAAFxV5wnULyf5cGvtg621rSQ/k+SLFzMLAACAq+Zpv2bm80l+K8lHWmuvt9Z+tve+SvKZJL+W5OtJvtB7/9rzmwoAAMAme6pnUHvvn36b11/LO3wQEgAAADyt8xzxBQAAgAsjUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUMGigttZutNbunJ6eDjkDAACAAgYN1N77q733W9vb20POAAAAoABHfAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShg0UFtrN1prd05PT4ecAQAAQAGDBmrv/dXe+63t7e0hZwAAAFCAI74AAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQwaKC21m601u6cnp4OOQMAAIACBg3U3vurvfdb29vbQ84AAACgAEd8AQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKAEgQoAAEAJAhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACghEEDtbV2o7V25/T0dMgZTzRbLHP73oPMFsuhpwAAAGys8ZD/8t77q0lePTw8/FtD7ngns8UyN+9Oc7ZaZ2s8yvHRJAd7O0PPAgAA2DiO+D7BdH6Ss9U66548XK0znZ8MPQkAAGAjCdQnmOzvZms8yrWWXB+PMtnfHXoSAADARhr0iO9lcLC3k+OjSabzk0z2dx3vBQAAeE4E6lM42NsRpgAAAM+ZI74AAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACQIVAACAEgQqAAAAJQhUAAAAShCoAAAAlCBQAQAAKEGgAgAAUIJABQAAoASBCgAAQAkCFQAAgBIEKgAAACUIVAAAAEoQqAAAAJQgUAEAAChBoAIAAFCCQAUAAKCEQQO1tXajtXbn9PR0yBkAAAAU0HrvQ29Ia+1/JFkMvYML81KSbw89givFNceL5prjRXPNMQTXHRdpr/f+3ie9qUSgsllaa/d774dD7+DqcM3xornmeNFccwzBdccQPIMKAABACQIVAACAEgQqz8OdoQdw5bjmeNFcc7xorjmG4LrjhfMMKgAAACW4gwoAAEAJApVn1lr7XGvtW621r77ptb/aWvtaa23dWvOpb1y4t7nufqm19juttf/cWvu3rbU/OeRGNsvbXHP/6NH19pXW2q+31n5wyI1slsddc2/6tb/bWuuttZeG2MZmepu/5/5ha+33H/0995XW2k8MuZGrQ6ByHq8k+eRbXvtqkp9O8qUXvoar4pV8/3X3G0l+pPf+o0l+N8kvvuhRbLRX8v3X3C/13n+09/6xJL+a5O+/8FVsslfy/ddcWmvvT/LjSf7bix7Exnslj7nmkvxy7/1jj/732gvexBUlUHlmvfcvJfmjt7z29d77NwaaxBXwNtfdr/feV49+nCZ53wsfxsZ6m2vuf77px/ck8YEOXJjHXXOP/HKSvxfXGxfsHa45eOEEKrBp/maSfz/0CDZfa+0ft9a+meRm3EHlOWutfSrJ7/fef3voLVwpn3n0OMPnWms7Q4/hahCowMZorX02ySrJ8dBb2Hy998/23t+fN663zwy9h83VWvvjST4b/0cIL9Y/S/LDST6W5A+T/JNh53BVCFRgI7TW/nqSn0xys/v+LF6sf53krww9go32w0k+mOS3W2u/lzceY/hPrbU/PegqNlrv/b/33v9v732d5F8k+fjQm7gaxkMPADiv1tonk/x8kr/Qe//fQ+9h87XWPtx7/6+PfvxUkt8Zcg+brff+X5L8qe/8/ChSD3vv3x5sFBuvtfZneu9/+OjHv5w3PggTnjuByjNrrX0+ySeSvNRaez3JP8gbD9j/0yTvTfLvWmtf6b3/xeFWsmne5rr7xSQ/kOQ3WmtJMu29/9xgI9kob3PN/URr7SNJ1kkWSVxvXJjHXXO993857Co22dv8PfeJ1trH8saHcv1ekr892ECulOYkHAAAABV4BhUAAIASBCoAAAAlCFQAAABKEKgAAACUIFABAAAoQaACAABQgkAFAACgBIEKAABACf8PxYCmcxInUnIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(16,10))\n", "plt.semilogy(mags, cdpps, '.')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXd8VfX9/5/n7uwQspgJhLD3EkQFRXHWrWjVVmtr/bZW26/VDq3i6Fc7XNWidSuKRVFUQEX2XmGEBEL23ju5ufuez++Pc3NH7s0AQ234nefjwYN7zueMT+54nc/nvT6SEAIVFRUVlbMLzffdARUVFRWV/kcVdxUVFZWzEFXcVVRUVM5CVHFXUVFROQtRxV1FRUXlLEQVdxUVFZWzEFXcVVRUVM5CVHFXUVFROQtRxV1FRUXlLET3fd04Pj5epKamfl+3V1FRURmQHDp0qEEIkdDbcd+buKemppKRkfF93V5FRUVlQCJJUmlfjlPNMioqKipnIaq4q6ioqJyFqOKuoqKichaiiruKiorKWYgq7ioqKipnIaq4q6ioqJyFqOKuoqKichaiinsfKGu0sDO//vvuhoqKikqf+d6SmAYSF/xtKwAlz175PfdERUVFpW+oI3cVFRWVsxBV3FVUVFTOQlRxV1FRUTkLUcVdRUVF5SxEFXcVFRWVsxBV3FVUVFTOQlRxV1FRUTkLUcVdRUVF5SykV3GXJOltSZLqJEnK7qb9NkmSjnn+7ZEkaVr/d1NFRUVF5VToy8j9XeCyHtqLgYVCiKnAU8Dr/dAvFRUVFZXvQK/lB4QQOyRJSu2hfY/f5j5g+HfvloqKiorKd6G/be53A1931yhJ0j2SJGVIkpRRX68W4lJRUVE5U/SbuEuSdCGKuP+uu2OEEK8LIWYLIWYnJCT0161VVFRUVLrQL1UhJUmaCrwJXC6EaOyPa6qoqKionD7feeQuSdJI4DPgDiFE3nfvkoqKiorKd6XXkbskSR8Bi4B4SZIqgMcBPYAQ4jXgMWAwsFySJACXEGL2meqwioqKikrv9CVa5tZe2n8K/LTfeqSioqKi8p1RM1RPASHE990FFRUVlT6hivspoGq7iorKQOH/C3F3udpxOtu+83VUbVdRURko/H8h7tt3TGfHzhnf+TqyOnRXUVEZIAxYca9ptfH2ruL/6D1VbVdRURko9EsS0/fBPSsyOFbRyiUTkxgRF/4fuadQDTMqKioDhAE7cm+xOAFwy/85wVVH7ioqKgOFASvuSr7UqTk55e/4IFDFXUVFZaAwcMXd8/+pxJ5f8sL273RP1SyjoqIyUBiwNnepc+gOCCFTVPwig2LnceToHWTb/g+3fgYPLE5n+dYCJnoOLazv+E73VEfuKioqA4UBK+4Ayxc/SHPtz4g3XkJJyT8p4Z8ANNSvZ0VOJLFhep7bmMdbS/rnfqq2q6ioDBQGtFnGqHXSUrccIVwhj3ly3XEk5H67p1p+QKW45J9s3pKGLNu/766oqPTIgBV3JJ/QChFawC8csZM3l/y63275HwzMUfkvpajoeeX/4n98zz1RUemZASvuGr8ReUPDlqD2SH07t09YHbQ/q6KVzTm1p3dTVdxVPJSWvkZd/YbvuxsqKt0yYMVdK7m9r0vL/tWlVZAUHnqN1h+8sou738sI2VZbW8tTTz1FS0tLyHY1WmZgU9Zo4Zvs6n67Xnb2/f12LRWV/mbgirumN1u61Et7MBkZGbjdbnJzc0O2qyb3gc3lL+3g3g8On/J5suwMafqTJG1/dEtF5YwwcMXdb+QeitA67Nv7xzVZ7ClsCGz1qLd/mGXv11QZKHQ4ev7OdMfWbePJzVsWtF+S9N+xRyoqZ44BK+46TW/iHizQWsk3+lq5v5Tb3tgbeE4v4n6mqkK2trZSXl5+Rq6t0nfa6us4smFdwL7OEXtl5YdBx0vSgP35qPx/wID9dkaLRhwOU8g2SQJEsEDrNU7v6x+OXx0USeMv7rW1taxbtw5Zlv3a+6HjIXj55Zd56623zszFBxBOt8znRyq/t5DTtS88w5a3X6O90TejE6L7QYQkDeg0EZWznF7FXZKktyVJqpMkKbub9vGSJO2VJMkuSdJv+7+LwZSVlTGqqZT9+24KKbgaZCYMDrab6/zEffHInUBg7Lq/uH/44YdkZGTQ1uZb5ONMOVRdrtBx+v+/8dq2Qn696ijrs/rP6RmK7moMuRwOANa++Cy1RQVA92G2oIq7yn83fRm5vwtc1kN7E3A/8Pf+6FBfePvtt72vHY6woPbzh+/jhvR1Qfv1GhfXpH3Fr2b4R9f4j8yVH31HRwdutzJi02r9nGaq0f2MUteuJAY1dzjO6H26M68ZIyIAqM47yer/e8yz1y8qq/T1gONVh6rKfzO9irsQYgeKgHfXXieEOAg4uzumP+lwBtaH6Wp8WZV7Le8dXxry3L8vfJyr075hesJx7z7/aXenuG/ZssUr7v4j+/5MYrJaK3C52vvvggOQFosjoGSzxvNhnukyzu5uxN3fhq7RKK/9vx8FhX8JON5ur6ak5LUz0EMVle/OgLO5V7RXBGzLIvBP+Lb0InZULujz9QrrfWaXzMxM7+tQ4i4QtFh8o0qbzXbK9mGHWzl/z96FHDh4dUCbv32/L5SWlrJz585TOue/hQ67i+lPbuTp9Se8+97bWwr070O0rt0WZIbp7m32d6RrPDO2nswyAIVFf/tuHVRROUP8R8VdkqR7JEnKkCQpo74+dJJRb+i1geFnQoT+Ez7OvQan3LtN9Kp/7KC4oSNIWJ1OZSLS3u4bXR+raOVXz7zG8/98ndbWVp599ln27vVF3Miy3KP9fHv5dmZ9MIuDJQfZueMOysoC2/37cPToUT7//HNkWfbub2trY9myZSxbtoyCggLeeecdNm/e3Ovf+N+IxROW+MXRqqC2PYWNfHwwMHrI5ZZP+UFa22Zj7p838+Lm/ID93UY9+Ym7uakRWXb36FD1Xk920dp6lJO5j6v1h1T+a/iPirsQ4nUhxGwhxOyEhITTuobUJQpGlkP/CRtKF7O3ao5yX1mLCBE9Y7FEMVrTyEfvvkFVVbDIALz55pve12/tLCJN20RbfRVHjx4F4OTJk972Dz/8kKeffpq9e/eyY8cOHA4Hzz//PIWFhQBsr1DqyR8tUM6tqxsdcK/jlS00mBW78+eff87Ro0dZvnw5Tz/9NEBAH48dOxayvwOFToFt8tjX39tT4m3blFPLw58eo9Wz2lZ1q5Uxj3zNJxkVQdfppLrV6n3vOqlptQGw9WRdwP7uzTKB35HiIxmIPhSea2vP5NDhpVRWfoAQZ9ZfoKLSVwacWcbpCjTtdzdyB7A4FWdr7urXqMm4I6j9UMbVnKOtxGFu4Z1P1vZ67yMlvtnG1q1bASVyp5NOEd+wYQNbtmwh42QJbW1tbNq0CQDZM8XvfEBJXTy0N726m8teDDSzNDQ0eEfuDkfPwuGWBZtO1PZp9Phx7sdsLz+9xUuEELy5s4h2W9/cLAXNBUF9crh8ovmXb07y+JfHu57GtCe/xWx3kV2pmM6+OV4TdEzq79fzyJos5j+zhdlPb+pb/7vRa4020EGq1emhDyN3rcbn1JdlNfJJ5b+DvoRCfgTsBcZJklQhSdLdkiTdK0nSvZ72ZEmSKoD/BR71HBN9pjpsdwSOznoSd7tsQJY1mKMKaC6ZH+II37llTdZe720g9A/948Mfc/XnVwftf3qVItSyLHP8+HGEyxNq6XEDS5LAKfsEUoMIGn128ubOIrKzfdGo/qPMTvF/e1cxP30/g6+zg0WwK0/te4r7ttzX63H+ONwOZn8wm+f3ruLp9Tk8te5Er+fsqNjBdV9ex/ri9QH7rU7fe/nqtsJuz7c63F4/R2y4zyRnd7nZXaDEo3+4vyxg/7Ivj9Pc4fA+OjVdJm3djdzpMnLXGQx9MssUFb9Ep2u/u/LTKmcXDlvvevF905domVuFEEOEEHohxHAhxFtCiNeEEK952ms8+6OFELGe1229Xfd0sTu7iHs3ZhlQfm7V1elYI6qwRJb2eF25x1o0ghRNE/P0oa9x4ssTWKuDP+yxWmWkH1O7l62fvIa0xSPqnSN3SWb+St9DR9ODCeDp9Tnk5eV5t3OqWr2v3W43LS0tVLZYAKhqUfricDt4M+tNnG4nLrfM/32VQ2M3D4+urCtax76yfOwun7g12Zqwu+wcOrweE04+zqjgi6OVQefanG7vefnNir07rzkv4Jj7PzoCQIQMCX56GOGwEmM3s2zvWzyUsZJjf13OY58cYUxLBeOLMtlX1IgQgoc+OcZtb+4PuvfazGre3VPCXzfkek0/gsDY9u5s7qKL32XVst/36lAFaGjYhBBOz7X/I0FjKqfB+w//isNffxmwr66kiJ0fvYcQgk1v/pNjm77p9TpNVZW8/OObOL69Z3/XkQ3r2PHhO97t2uJCVj3xe5w22+n9AafIgMvCsFoDRdQtdx9rrJXcfuKv/KCdbh06javrIA19CJt8J5fqcxmi7TlscUFtcISOUVJU64d8AcAyfgNAa4UizAWto7BH+gpZafUt4EziSFlz0LXuqN0OKeHe7Y66MvSeLufm5rJ69Wp0I2eiQcLpEdbXv3mdD6s/RCfpGKG9nNd3FFHd2Mp5hlL0bj1OrZOM44W0VBdz0UUXecP/bnh1N3nhf0B2RbIk8lX+euNUHvj3EaaPaSHWEcvYhjTC9cVsco7lgX8fZVtuPS8sne7t27mPfUp6hJ1Vj95Bm9mKJCS0XWLCT9a0M7M2l8WGqSBpONi+m23DZ7D6qz8B4NKacOhNrNWmc6O2kB9uU2LMlzYm8OAl6XyZGdpHsr9hPVETXuGjjCeYPCwatB0UhT/GHSsfAqKA7pOYhCxjiohk5NQZ5O3tTHI7tXo0QlZt7v+NNFVVUF9azNZ3X2faJVdQlp3JqOmzWP/SX2mqqmDaxZeTufFrAKZe3H1ajxCC3R9/AMCxTd8waeFiAJx2G7VFBQyfMBmA5poqtrythMnGJg/BYbFQkLGfypPHKc/JYvSMOWfyzwUGoLhXl1QQgQUnOhwYyMs9l4TxhxgWWY1R50AruXAL5c9SwiR9I/t6Sxxr9t3FhPByzp/3ccB1R7p0OLp5N3oTdoDBNNNKFC6/t9QYwoyT0p5Ca4Mi7o3WuIDEKJ2hHpxJPP/IctLSAs/T+gk7gN6vTs7q1Urd+pjdX/OjEUnU7stjE7NpONjAxdLFVLa1UNKgmDA09QXktBVwNVfz6ahPWffJCgBknZEli84H4FBpI1ETQKMzsyO/nvRHvkYXcZydzhUkyIojXOd3/zVHKgPE/WrjCfCMxls3NHBu2Lm0O3zvYXGDkqsQ7rKBUfl8fnlyK5eX7PMes3v+07j0BmA3Wp3v3J+3mzB/Vg6xvvfCKEOMLBEtS2TUKr4TSdfGI2uy0UaWIWkdHG5dDdwFdG+WEbLM4BEpzFx4nlfc68pCVwjtDllWR+7fN06bjbUvPsvin9xLTGIyjRXlvPvg/3jbX7lrKa4u5t3Nb7/qfZ1/cC/7P/uY2/78HEgSx7dvZujYCcQNHUZp5mHvd8Pa7ps9b3pzOSd2bOGuF15Do9Hy9gP3eNs2vv5KwL0s3ZQU728GnEO1sLiah/gX97ASAIcjgj8feJA3shWHqX9BMbfQ4nQaPVsSha2jmKCrA4eRutrUgOs6TL56ImZhJBSzyWQZL6Dvkq+1gIP8ine5nq8D9mskQde01tkNs72vR2ubmNUwy7s9NDyXadpK0tJOTyAaRyQpL+zt7NquOHy1QsuOIxvIyfwEgJjq0Nd+e3MWAJnlLeAn3E6P4/N2t4UFtQvQeHwcbiERI1nRC8FCp5Xyqlp2rc7nzb8c8J6b+nslSzjZmsyq3FUcqD7A+mPVXPGS8uPYNWwabo0NWXJyYM6joA3jsStm8ty1Wty6sD6Xe7ipw8iPzSausxiZn/uDwMYQM7KuA/fPj1SyLbcOIQQajYZP//6sty03/6E+9cF7O9Us871TcGg/xUcy2LnyPQDamwKrv3YVdoCSzCPe11/+/c/UFuVj6zCz/7NVbHj1Rd75zc/J27+bT5953Htcc3WV1/ZeV6z4jd75zb289cDPeuxf5qave2zvLwbcyN0Unwi1EI9iuoiMqQMbHG8Yz8mmdLSSC1DE2SVrqaia5D23xRZDhOf1jblH2EwE+5gZdA+XkEKWg1+AsshHLK3UE+/ZK7iEXQCMpizonHn4vjSplFPCiID2VHOq9/Wc9pGgD21u6AkJmVlkUUAqLcQEtc+vPweAI0IQ7dZh7jzPT/gmuBysXLuWF4+FATJTGqcwtm0sH8sCSVJ+DMnWZEqiSgBI0Ji5zphNQ8do4qOKeOv1bBJqLkDWOCBRuaZWG+h62V6xnX9+lsYiXTlGvYsNznE0JR70tq+I/DHyuGe4pDWOOHMpDe6h3rbm2HQGtQTGq3cyxO0bo0SZBwW9Owo+RT9e2cqw2DDe2FFEbZuNN3cVA/CMQUan1+F0+Y7VRZg5FWTVofq90xls0BmhtWNF70X5ZHfw59bR3MSJXdu822uffybomD0ff8DxHVuxtffdzVhTkEdbQx3R8Yl9Pud0GHAj9yk6X5ddQkO4x5brkI38LeNXWFwR3vavSy4JOLejdCyghCDqcXMZXUMBBRezk0QaCEWnySUaXwkErX9tmhBPhEn4HIl3spoYQn8JkqljLEXe7Ug6WMYLTMUXkRKGlTgC7fE6XKRQyVVs4XK2hbx2Jz8MywC/6CKt8NnBhaGN4kP7uMP+EXcajzG2TXmvLI4OwmN8K1fp3UrEisEzug/TB75XjYk+00pYZKBJY2X2Sm6K2Euqtpkh2nZMQz8KaJ/hkvlD1i382ZzJrfG/Bkmgx8kyzXOY582iIyzee+zExmKeyXyIhzPeDXw/3MpnsLiss8/B4n7PikMA/PmrHD7eXkKa01NqQJahhzK+RmNyt22dqDb375/2Bk/IshDs+eRD6stKTus6HS0tOK2WHo85tP6LUxL2TtobQmtMfzLgxD3SroyO2yQDkWgZ3DC1T+fJGgeJGkWUNd2ENEZi4Twy+IX0Uch2rec8PU7OJYPJ5KIj8Im/hO38mE+8210jYH5D6FHEvXzodbwCXhGfzTH0OJEQ/JL3uZ93A867m1Xcicfm3s2Dwx+H0VcmSOefwSvJXMFWfqVZzVB8oZSXmo6SGOvLwo1yBka5hum6H9lGJXzrfR3tiCalPYUIv4evPiYwEWtYWDHXxf3Juy0QRHvmGUtin2f/Ob4p8SJDKrnn3Ixt0k0B19B4/C235HVGMijivrR+KnN1yndHr20hv0GZBdzVbuT6DmWmZ7Naqc4/SXcY9IO7betEjZb5/umMUGmqKmfv6tC/5b7w7esvYzOf2sytr+z99PT71VcGnLhXD55NsW0OFhHBeHKB4OnUrMSjQfvsYb4EJG03IYdJKMfESmauYlO3oYk6XCxhJzfylVfwQRm5n8thRqFkUs4gm9iQgiuIpZVQZSancJIJ5CM8H80wanmEV7iSzUSijCJG4svUHIIv+9LVByubw9Toff3LlkYm+s0sOh8Ot/IFt7GGxeziFfEq51Zd6D1mbFt6wPU0PSx3OLdurvf1JZWXML1pekB7jLVzWiqQkHEZAx1N9yfewkyyvNvDDb6HgU6rhJNZw6sxSW2kGDMA4RX3ZlMECMEFLUO44diDAEzUKQuj36bP551XVjDeocXgEf/48BwaLdYew9T6ItxVlf+moWFrr8ep9B+NFeVUF+Tidjk5uds3G28oDzaTngpt9bW4nGdmJhY3dPgZua4/A07c0QiaXcNIpplrpS+YIG0hqUvRynB9YLikBhkJmSvYwhBquxXtO1jjfT2bLIZSG/I4nZ+g38qXIY+5gH1cw0YiCI5/T6KBX/M253IoqO0GvmYpvnLFnQ+i2X4it5Dg+G5QxD2ZOnS4PKalnh2Sd7Y3cjPricBCAg2MRqnnEoWFdEo4n4PoJScmTeAXfDrHWcYLjKIME92L4SBHV/t3ILcefQSAcznE47xEHMFRBAv83qMhhpMYpHaE5AZJZiSVLNBs5+6kH3PVoD+TZtyDhI5dJVVo5o4kQkjM7khgkM1X6uJck/KAMKDhKosec1QhTl07V8lmOpJG9tjfvoRFVtd8Ruaxn3Is65e9HqvSM0IITuzcitvV80P1o8d+y8pHHiRj7RrW/+O/r5CbzmDwvh46djwAYVFnLM/Ty4AT9wkxE7DLqb5tzXb2m+5jjuSbTltdygpNF+mVqfdjvMRS1jKXTH7OyoCRu4RgLkdCilSiJYy5fg7RTjodqwAj8C0sEe53jYsIXMLPn0EoIVTplAChk5e6mnv8SaMspGlJh4t7+ZC7WcUvWMFF7PG2jaSSQX7imYBvBH8vHwT8TV2JJnBqeq7n2B/zKXezCoCHeI1zkx8FYCo5jKP7rNNO9JKZZbzAEpTomc6ZU3fMjVzFz5J+RHzSewjJzU/4mEV6XyJJpLaRiYZdxAgXN5g2c1P4QWymWlrifDO5sfgiJYTkJjwim+FxHzObTO93ICXKF3p58pNRyC5ldD8m7bch7e56fVzQvvr63pNhvjecVtj4OLbmep5bepX3X+esZcFHC/jTbp957NjmDZQcDR6I9IWOlmYayko4tP4L2hrquj3u7d/cy46V7wbsKzx0gK9feY49HwcvceiPvUMxtx7d+NVp9bEvXPaL35z2uS6/siFVeYpOGcLCuzu83xhw4l5TUEm6yffG6D2RHFM0SsTDrRh4oFGJgBmpbaFz9Drez1npL6Z382+uYFtIZ+TV4e9wBdvQ42QC+cR5RDmh+/L2fWICyio/WpxMJI/Fnmgbf3oSd1Di6ruOzN2ej7PTVDONHEBxzv6Ej3kAX7bcEnZ4X0fR0aNJZzyFLGQvBpQvqcEvFLTzvYjAygyU+jDX8w238qXH9NQ958V8ELDtP2PpiVtYS3z4waD9Aoklsb6Y4qv5FmvscYS2u5GfzK94lxs167iKLVzJFkxaJ1em+XwOtiYTDrPiRI6Nncd5C3aTmHhFwFXGpD3cp373hiw7yTn5R2y2U4+YOiX2vgK7X6R58/KA3e1NygO/zdHG5wWfe/dvfP1lPn3mcY5sWMdzS6/C3Vn51GmFqmATqD9v3nc37z10H9vef4M3fvkTAGS3m13/XoHV7HuINldVcPCL1QHndjoz2xqCH/pVeTk8t/QqWmqqGTFxCgDmxv5zUv5gfGDgQk1h6Eit02X6pVf26/VCMeDEXTaDIDgO/W7Pn/JLTIxyRXmlqmtxLoAwvxH2cI/z0Ej3afkXsL/PwtMXOkVXaDq4mfUBpgdfH3suE/BLVrCMFwP26bs8EDofEBeEMON0zho6mUxe0DGdXMoOLmQfizyzESPd2yEvx2dvTiV0FccpnORPvMRE04Zur9MbtxveCbE3MFppJsd5lJe9sxQ9TpbxAst4gQUcREiBM6ZwrGiQCdMJxkX7CYrnK6TRKN+qKZNfZsTwOwFIH/MIg+Mv5HT56quvyM1VoopaWg5QVbWKnJN/7PU8IQQbSjbg6q5QWf5GsHeTfOfZ3zVL22mz4rT7fhsv33kzVXk53u0dK5QV0Lx1Vdb9Bl5fCO012C0Wnlt6Ffn79wRcs6vNete/V1B89BD716xi6zvKimgFB/cRCq1eeai6ncEP5xM7lO9ZceYhpK7Fg/qBSGtJwPbRDesIj+3ZzNgV4fe/jMSR6KmsSb6a5Sk/Y2fhdxsg9oUBF+feoZUQIZKMIruYKf7oCqOM0CaPX7AiaJ+b7ssYnE/wKPFUkJHQhHjImHqo3ng9pz6t7yq6nb6BJk86p5lwxlFIDcHllk29PEwAzuUwsbQFPXii8YnIOfQ8kvsBG5lFyOV4T4lQTvHzo0NHIv2S91nLYsrxxc1fwi5OKLXvvLjRkmCyQEctrU6f+UX2hEr6r5nqTbCSJKTvMEY6cOAABw4c4LcP3YZG48nPcJmRZTsulxmDoUuEzokv4eM7+Oamf/Jwxl/4tTaJu6PGwzV+WZDNpfDhjTDhB7BUmR1ZSo5iemchmmm3QJhHpLqsjfDBH36NKTwcFkFaRQQOq4V9n63ytncKtctup6mtld2b8rl8sITO1kp1pVJjaO/qlQyfOJn60mIaK4IdmvvXrOKGO68FoDo/l3/+9IcBoYRfvfx3LrzzHipPnmDtC0oyWSibe6bHBNOZ4t/fVFhjg/ZZWoLLgvjToQ1n3vBqdsfeynttKRjdNuxaU8hjixs6WDSuX7raLQNO3MeNjkTkBb9hkhQ4grne+BV7yQtwfvaEfAYnMXmMYqwoplWkMEhT4t2fJPo3zKqruCuRPMIbz99MDLfyJR3y4NOes030mJT8+V/eDHGk8sAYTwEnGYMeJ2HY+kXYT4eL2MOaLksBGzWB7/9YislkGrva7qLO6QBPUpo7hLj7091+p7MVvV5JKjt48CDr16/n4YcfJjw82N760cpHWLz4XABkt4WsrPtoaNzC4osCfReuzPfQAU0Niu22prmAnEMlfLXyKu555H6iUibBCY9JJWctVB3BMWgcr/7uUaYPSmNO0Sfkpk1iV84ChrYH5iFUDbbi1CmmkPOPKTkF1vbg0X9j3lH2rV1LZUMYeQ3ncUtBIZ++qJh46stK+OAPv6atvnv7utj8NDCZltrghdBzdm3DYbNhafUJqctv5G63dLDi9w90e+2+cPWwE7jCEviqoPs1JXbWjeq2rdI4hEi3mb2DziHVUopVG0ZG7CxsWhNvA50Bct0J++M/mMhdC7q/fn8x4MRdNOSHNMuAi3i/aXms/g0uP4XrTuHUaoj4c8h8PbMiP+u2vYNw9rZ9QVrYZjC8dNr36Y3OUMlOdLgDTDedzt9wqecRSH/R+VCxYKKVKIaEcJi2ERnksPXHhbbPD+iekBABYasAUYaSoOPOTa7k86ZlgM/+67Z1irv/E9E366qzhnYEHz5yG+fMVcx5hw8rBeKam5tDintl5USKS15W7idbaWjcgtOi5evlLzBu/vmMmj6LutqvyE46ypxKLRrh60VWizLLaFrxM6IiWjjUOJTMlplcuuA4iW9dyF9b52FCS15bPGETmijKdiBFG6YzAAAgAElEQVSjoaKiMaAPQxvDKE0K/A41VQaPvj998UXCtQ5AiQL594uBtvuehB1gY/WYHtsLM/YhaXzvdVVeDnUlRSSmjiZ3705aa3svaS0ITjK/b/xejifdxHsVY6iWkigcnsiIUSk0NLZQ32Im2V7H2I584hzNZMTOYJCzhd1x84l0dWDTGLHoIki21VBj8s3q8iPT6Y0ZI2N54upJTBmmPOi7Lgpzphhw4i7X5COLsUH7JZw8RliIM8484dJ51LGNRI9zMZ/UAJt2lGwkyahBkrrPdvuW871RI2camXC0PQhqfxOOLSCSyJ9Qwn6M8UxFGZm+xE94kDf6pQ9dw1aXRgWHzQ0zeJxykm/U1bZnECkLnuTnx0t4ftwIInQ+E56ExM6K3SEMXdDSUkBpaR4pKWPR6ZSfWk/LMArRaQdXhO34irHAZk5s34xGo2HBZVkwQoM5QoezxszQehPC2I4mysU8RzEdbVqIgG11acRPbuLk2Cj+HKFlp1zJuSKO6dVGhBvqjnafjJVSG45feSYc1tB1yy1ufcj9faHdFXpE649/+WWX3c6K393Pg6vWodEGS5YASsNGohMuwtxWNiYspt6ofCJxjiaaDEo008v2/8FbIUQG9FBSYQfCwBhGgzGB7OhJQddvMvgGk/7CHoovDI+SMnoCxec8yYQxozDpA829zZ98gmn8BMKmTO71PfiuDDiHqoZ2ZIJjRAUOJvax0FR/E6fTYXD7PqwdnEO7PMW7PcKT8u90X8oR91Uhr1FJ76ntXXHIab0fFALxHR6CjSLltM/tCwWkBJjIZDTUEk8TMWxjXr/eq7sZQbVW5tOkxYxLaiFSZydJ18Erx8v4oq6FtJ1ZnOyw0jlyl4VgRc7KkNfJzrqYd95ZSV3dRvDkO2RmZvK3v/2NTd9+S2tdYB5F58IzJy0u3qg3BLTJskyROR6bLZwcjZ76NQUsOZiEDESc00rcIRlpuxGzx+k7fEEtQsABaxiX7kuiLMmKxeAm79BInGHB9mR/frShL5/xqY0+7ZIBu8YQtL9VF4VZGxHijNBotcG+sW2DL2Bt8pWsGXINK4ff4hV2wCvsXYnEwmSp2Lsdpwn9EJvs7j5K5sP5DsZ0FDK6Q4nEG2ktZ5qmiNiS9czIfT5I2AFq/vQYJTfdFLT/TDDgRu4YdGRbYXiXh3+jNo8xuusAqLX/o19u1e76AVG67pffa3ddj4QFp0gB17WgVcILjXICTY4niTIp/Wl1/dhzhokE572gDY68aeTUPPEAFvcFNDvvRic1ojG8waA+lB8AECLsVH+bXjLdN7JAswGjpv9s5xUke6OW1rHYm6RVK08kRR7Lqzql4mcEFhahRFZ056QG2M453SZ69YUPohSH8b/0d3JrxwqKzXG8mfUnhl+4DYAnCqp42Kj8dDY0tlNmriHUx9fertitj2X9Dx0di4DhXvPMrj17yHzrHzDBVyW0dNswUi6o4IsjNqbl+hKqBOCOjKGsZRSWA4P5dNSn3Olpq7RGY3MBnvDUWhGFrNUhOyWyV6Rzq1MRmCFNJt69vJSbsxbjHhlLRH4mmh6Sg9xoMOsiiHG145R06EMURJORaNHHMMjZQqM+jsHOJu/XyiVpsWlMaIWLirDhfJO4BIDrqr9gqK2aLfGLGGktY4Nnf4K9nutqvsTYQ20e2e1m16oPgvaHGm0DXCCOYHGGkWKo48eGDaRJVXxuv5IEeyKLot/AILkp+DIRp0XH+MdmIOWtR1z1D3IizyPr3fuJ1duYTg7/LF+EU9IR72ziwQk7WeueR8n5L3DuxeP5zbp57NTNoIjRgEAIaMiOIjq+FqPLoSzT2NEAGi1EDw3ZzzPFgBP3JslKoV3mo4aXuDX+AdxCi1ZyEyf5RkFJxvv75V6trp8HiLtDTkYjudBJytTdKYZhcV/qee0zFV3oOA+Bb9rqEqkB17W5p2HSZlJrf5kk469oJwIzEdTanyPJ+GDAsWXWFeglN0NMdwb1T8JFmWMyQ/QaYng1qL079Jru07LXcRGDZZn5mm0h260k0uy8l2Rj8BJ9H3INRpzciC+Z5KiYyg5pFvcTKnRR4U1u4Y/iAwxSA1ZMdI6Ka+SpTHGP5IROCam0eGYcR5hELG2M8mTUruB63GhIpIFWoskljUNM6dbR21eyGEJWwsPoZQdypZGLGvfRrI8mx57A5carmEAqOS2TYeQFdIg7iCC02U0ICY0meJbgjAk0j7QUxtBYOYTznHbcYZHIOjMalwN3ZAzWEemd+k1yg29kE1+QRNEhI0XTYGJlPfvFUhgLx94OFuM7v07BnGbCIbS0GZOJFs0YZTuHYmcwpe04hRGjsGlMxLlaOBmRTkl4qvdcjXAjS1rGmfPIjQw2iwIk2WqpNSX1+J6uGXKN93VO1Hjv63pjAsVhKaR3FLJv0FxGWCvQIDPUVu19iB/+6gva6n2/884H/DBrJZVhw5hvz2SvcZq3/T7rJ0R96yAs3k7qxYp/oaXlDlqAJTFKlI3TokigcCvjHWnd/Uz8dRYTTZ7cEy38fvly6p6a5i0M+APtPlgyEYArh+VyuFF5uGiMgtwhkYhVUbQUHyddswjq/NYGXtZz3kd/M+DEvdVjt25yjcQltOgk5Udj4MwXbCp3PIGJOKTyF0gck4vNPTfkcYIIrOZ68PwGqxwyQw0aXA156OLH0uB8CpwCrSf+utPRJwh2tO3piOK8SB1OeRhuMRiT1ldfRZJcdEZThoXyIPnxDQtDVMEMJtZyFUM0o7Eb6hDCAGgwaX2rRU10jUFIvkzXMoYykipOiKnkS6MBAsS9jBE46N4+qyxvKFHpeIw842rsGInyVN00S+EBA2KBxDP8Aid6rkEpSvYWS70hjv7llNs8qy4BNBHjTUDrjU3umZxr07HH5BNHp8bAt8YFrM3+HQLIGT6Y1+Lv4OMYZWZ2u3iHehKJoISYsgtpHbmVUrufaUmW0GqDQ01tQ0cxjRNUkUg98TiiB9GRPBo9biXE0ukgsuAYsk55/9pkI5GSg0sOJONGokUfi7bDRWd4xolhijlCCMgaOo9shpNiLedQ7EzGaOsZbi7FKUez2z4WEgNLXWfEzqI7xktlFDEEB9puhR3oVdh7Y0/cPDYmXgzA4dgZ3v1T2rJJtNfx8m4jt+lj2Rm3gAnmXDYk+qq+ptor+SjmL/xvYjyflSkz945dJqJwYLXpuG1IEs9stoDHSuPzb3i2rX5RQS/6TKoA4eFGUiO7lMaoPgYGxZzkP4PUWZVnsMsiIWqPI7sktHqlPas2Ex3gjhIIIc64Y3XAiTthPtt0p7BD70k/rfowYpzfbVFbYTcjjMPoKNYhj1AWAjDvfonIBUpoll2e7DFXSDi3Pgke01qW1U10/peI/K+xXPkqSXpPiVmP+nfGbMsi0BbeISKIKlgP06+h1vEvyh0y86MDF+LudDs55DnotIp4F+y8iMHzxjNIvxynPBS9pop8UkOK+zHGM9Yd7hXwFqeW+tp8dCP/CoCU/QDDZivhnBG2vyAfX4dF2MFjTVjLxdwu1qBx3AfGnIBrO4WB40xEDlFfpxNXZ36BSGEvisAcELMYTRnxrvNo2/sKLPIJkd0TKfUVF1LGUMoZAsAwdxyVWsWhfY19Dl8YfbkJXRPZbBiD4vrdaCiTE9GiYYFNzyy7jpdjfE7grOjJCAEHZ8bSHiVxDR9wW/u7NERGkNxqxmrUYDdpaR25FVlAa9E80tP3otM50GrdjB23n9radJKT8zF3xNLclowOmeskJZHrCvmvnBg0B+wQL5kZranlfu0nPJdwBSIqnB+5v+VB571KQblRs+lKpKsdl6RHa5PowNi5pAENHvtzgTuBgrAEDDhJkyopFMNCfh4GbR0Ot6/OeALNfGP8Patci/idS1ldaEn5bvYlTybO1Yg9PIJqESjqEa4OHq5+mxq3kayoSewavACd7MSl6dkJ26GLDLk/K9rnzzoWPZmy8JGUhQfWATLYrVib9Gwc5RsgaTpHPmYtz68305gXBYs8u+pMRAz2fQdqyk6QrJPQer4rnY92pySRl/sFnfMBhw4qh4SRmvkR0kFlZih5foU6t8ywEhvHrg7HdERD/tFYzGEmmu9xMajZwStbl3LvJTLt17kxFPyZsemP9vh+fFd6FXdJkt4GrgLqhBBBLl5Jefy8BFwBWIA7hRCHux7XX+iMoR0kvbEzcSpXVZ6aHdZVl4Nz+BD0GiWEsOPYv3HWNYHTSsfmx5EiEjFbjhMhu5A0OvIsTxCnUx44ktv3xZm740Haz+8gMl9Led4WSlIWcU6kDtkj7m55EDc75wdZkO0YGF2iiDuA+ehKCloTiV7yEGHaPbS7rvee0+R8gHCPuLv9YoRt8lTKna/TaNzMX/k5P7LPJNn4cwCqSaTceR+0RTB1sGcJOiQMTbk4pVacZXsJ0/uuJQjHkb8Bfcp5AMjCRL00mJflu7lTDAUCxR2RRFNEFBZtBLllcxgXF5wM5kbLzfb56NAQL0fRoGmnnHSq7Gs5VPMJDYNDz8jsGDnEVBY4x5Ekx3JE53OORXkekm9yC0OpZX6XDOBPuZzb8KXXv8ZttBPJtdI6kAQliTrC7IIftRt5P8r3OX40ewIFLaNI0xSTHFGPOVKHCTuzMpVZwd2JP+fWcZ+ikWDMmMP848jPSY6oY4IrnkZbHJ/a5kJJ4GzvCc/s7oTDVyWwQUTygPQx52mOsSFmJk/p3wXgSecdtPjNSPwx60LvB0h01FFnUAT7Kd07LNVt4//qrmb7uArmNy9k1sx1tEp6ZKGhsSKduPgyjjVO5PKxX6M1a+AI/CByB/FzFT+L1C4x7qNMWpKiCYsaQULuTkZn5fPxXWacFdFMKo9menY5myaPwuBZdjDMbaO9F3HvCxZt6JosetlFfWY04yc7fCmHOoErTuAeLGgoi6AmbQbJs99DyDqavogkxxiDdq6bsAwNYp+bhsgoEqe18c7Qa3l62ALOqVvN85UZTFv9c+99DC4YVW6Fcl/45yCDMgjQJzk4nhyJebiM+TKZ6M+MtF2vPCaaBxm4DSftymQPydb70p3flb6M3N8FXgHe76b9ciDd8+8c4FXP/2eE2sHRPLU0jh/kFdHRpCfC3jdzTCUjyYpsZIo5OAmnO14fUs69zr8x0ng7AKsXnI+m2YzW7aYyMRkhSeyedT33WV6m6NhSqoYmo4uPYWqLm+ips8gaP4Lpb5ehdds4PGoySUU2hpetp0TU887MGUQmTWF6x6NUxe/CnWNktlHCfzxl0+vJHTeS+iFaTG6JIY6TPHH3L1jaNpnUjokYtIJmg8ASpWFsC2TrDEx2OZDNtdiPfwbT4cggLTfPjOLn28FCOFX19Vh5gFHDXyLaGkVyaS1NpWtBcR3w/AWx3FNYyu9njWH7Tx/jxoNf8YrlL/xzzM00VR5j2JU3YApLIn3QLEqNSewLc1ITPZg1hghM5tmMNLvwmMJ5P/4pVkxRkjXSDBWMawoWd5vWyJ7kaB6dGkZa/XSMzjbaIuKJKNdzOCaNovjZJKIn1txAOmNwO5v5ZjDUJw5jfHUp4e5R5EdpGVFdT73GRrjdxsZBRjIj0sh3DCfW0s5sSyb+Ca3RI4vxXzTraEQstSKGKyw6XDoNa87RMtjZQmnYMNKLixhRVMUe1yj+6FdMKyW6jHpLPBZXOCWmHwKwuWwhxa0pFLWmeo/Lb0ljZ+X8kN+vOD8HuAaZaVIhKaPz+LzwSm+doNTIcjonGRMjighDw6y87fwibQ832R/joFDs1tHONtr00czSlTNJW4uMhKnqELp2icK0SBqk8wmTbczRKCGmf0z8knOmxIMnW7tTMkekVyMkGB5fCoBRp/y+JL+Rh4gSTL+ns1DfARpcY4mvqkdjcnNioov7vmrC4JbRyDKd/hP/uJHprZkkOBrYmKAsLj2xPYdGfVyfzDqFEaEjxPSyC1uLjk+qakj17Gt+yo7L4K8P+71L7xY8CHjMoS13urG9YcRQKLHzTiO1dTsoyvgcl11DVf5QFlwmuLHDzH1tbTgMGopSwpmQ046t2kh7hx7rDEHn+j6tw32S2nZ99/kZCXmlMKXb5n6hV3EXQuyQJCm1h0OuAd4XyppW+yRJipUkaYgQIjj9rB9oMwwDO6wdO5oR5XfzWFHf0o+Lw4bxzKif8VHWb5nfeqz3E4C3L76Bt4EajzVjzYILaDAEh0U8EfEYLPZtbxiihwm/VTb+4HdgSBP9Jco/T8La8pzFVBiTub/8Q6x6A/f++i++Q2coCVAbQ/Y2gtiOZ0mz6Tn06lQWNh1kVdZv+SZJSRj5Yvp5NEbE8JrHfru0JozNcfM8f8/F1GxfCIDZZOLB3/imi6vnXMGn4jKEpIHQM3kf4cMpTMQr7o9M8WXhlSQkQBN8krSEm2p9i3g4tDoenaqMtAsTBgOKk/HpGAClT4HLZwzxvtobPob9HXZciVGQOFkxpMqAQYPUNha0GoROIq2ykP+t9JWc2CjO4d0pV/Nszpuku6r5tnE+MhokA4yUa1iTeT9TzAUkL9xO/qjR/GbYE+zZ+n8BvShtC10e2F/YQ3Hz2DWE66y8e+KHHDb5yh/8Uvs5D+pXczAhhqtGbyD+uA6a4NxhB+iseffLma9jC9OSNMYGufCieIlna68ixVKGVu8mUhdNfJqTOe7jvCXdjNGuR6e1cI3JQJruFdIdFRS5hwT0J7oepuc1oncLds6L4/x9imkre3wUtSIKTaTi45IEdBCB0WFHZwh01safl0fDeXAjcCNWLHdL2M0SsxLzaM5RzKjhse20dCizi/OblPozneK+qGEHn/k5Wk8HnezGbdeSb/PFpEueJ5IbDRIiwDYugDZiiPH4Ypp/pvxNg4AUi46ilHAavk2m4oqJPDQohwydzI4pym9H6xZkTYnGPNSII11gr/GJO4CuBozZGjoulonocGHKlmg8R3m8JS7To6uTKLprTIgFPvuX/rC5D8P7cwagwrPvjIi7s6mdB8Nf5jnpj97RTW9s1c3BUCDxw/p17NJN4bPhS/hbxd8B2J2Sht2g4aL8fBrCI4m3mMmOGs038QuCriP/BzLLfjHhMea0ZnF/ec9lTkPREjGLQ56Q4e1xc7hixnIORyle/ZqY+IBjVyWHzt/1/xvDhZlJZHFQms9PxXKymM55bKeYUXwm3QLATHGQw9KcgGuUhSUx0hoYw233TMlNbjtbYs6hJGIoP6laE3AMsmCQaGByw2HaBw8iQduIrslKdusEJKubuopYTEYHGkmJ3LBZlR+yrrjnhKyXuZS3WMRx090AZBeMpUgM5QaeIF2q9MbVT9SUKg+Hzst5vG5/1T/mvZaQAkexd2t9zmPX0HC0dVZco6KQ3AJNgw05Uo97eDgYtIgwLe9rPFFGI8DfBTLJU5bi08a72StP5qmwvwE11DaPY4KnxHJb1RgK08LQy+WAmWGmVq6adZytmvFkj5/JZazjuh1KOO7ssVtona9hsKOeSzN2dw5SA6p/PuB4neyc673bncIOMPlkO7fPe4cERzMb+Rkmh8w90vtghFVP3Muq667BEOaiYmQMi5w7qWEIwyPy0eNg89RLcaEji2kQZofjggpNChpP0EP63Aqubi8myzGPbfIM8tImcV71HlabfH05VYrHjOL+xY+QRC0/yNrDy4ZXuL/+d1SGJ7N70EzGWEq5v/J9lqddz0mNL3QyXtQxnHKOMZ1ERxOXOr/hg+RbiaORuT/byz4W8CPeZqHfUpMTcttJanCwae4QLBhwegIGBBJDf6HE8rtjBUgwT684YnN+ERgGWWwyDwhxD6V4IQOQJUm6B7gHYOTInhdG6I6LHbupDj/Eb8WfiXYrBssmXTRxruAYbwc6DLgoM89CYwJDq5YOEonwS14x7mrHtdRFW7WOvPF6MsPjgTbOKfiGVdnf0jzJ5+R8UjxMhGgjm6mkUUAtyVgIZzpHkJGQEMhoKGEUoyjmGNPYw/lczRqGU04u4xhEC1G0eitTlpBKCaOZw36KGU00rQgUx5JOuPhQ3OBNpS4lhWFUUEcSD0lKqvoD4m+UkkoZKdgxMZ4TTOA4X3ADZVFDkRDcJFbSQizx1JPPeA5KSjLQNWI1kznGBHzhWo86HiPW4Eunj2ySMccp4nchSu30mWRwg1CWEqxiGIdtM4mz1fNkzKMYZSvHJhqoFrF8KG7wXselU2Y8wxy17NFMwhTTBFUgyTLTDmwjt1lJ47YCB1FcOydDTBVs9lClJ3pGoKHDL3GrM6StjUgOie6rN92a9Q0fTb2cEikNE0phLPuSYWB3g0EDLsGf9vjirsen1pM9xS+aJL3vCzJcplXEoyY6jrzYUVxwTPFfLGr2icrllQcZPeprbhYbuKCzrIQIZ/OoaymXhlBtS+U+TynnK3Oz+TZuPr8b+wiXstR7jYmaUu/rQ/t9+0MxyNlGqq3Suz2j7QRHoiey9PHXlJleM0xI/ZLVhxRzVfLCYIe9VupAT0tAuuRXw2/n6pynuEG7k23yDNZdovRj/tffkKmdhCxOPbeyNnoIlcmxZDGZTTlKXsk/Cv7i7dfynKeYas5n+dDbwC9nqkFKpE2OpWrnRd597y28g3qSWI9S4OwFfhdwrx2WO0iijCfkZ8mXUtFIVgw0cYLJXPiqb/m8ez5bybzBSuXL9vAIDkxU3LJzT2TSmBO8ME1/0x/iXgF+MWgwHAhZkFoI8TrwOsDs2bNPK53UPSkZKiFpfyQpVgtowOqMBSlQ3A8ylaHUMoxahpe1MaLqADsWLfS2b2Mei9hHxG4dYQe15F4wGG21gCtlOl7TM/SYIgDf/Krd8xFD8iET6dHFTNdsI7HBjiQXBswdXFqJ7QsGk0Yh0W1Opkcc4ebincgaidLhYYwz5DL+ZDvFKeHYw5RpWqpH3gEme1ZbqmG095qzj7TQEa6lYqiJqXIBrTF6hlLFk+J3gOBH+w+idQsOTY2hI1KHvU1PWm07s+pPYjk31IpTa5UBqWcr/YSZ9hqT95sgNo7kaJ0yyt859QrM0Qau3bYaW3sYGJ24LFo0GoF+jAtTjB2zIwpTYw0W4OC4aRS1pHKwdiaXp27E7jaQHFGPUWunoHwsl/ICbsnG32tuYHxtGT81bgUn5Jp7r8/RHXcnfMCIMRXsLZ+PaBxMqv04wye2M7Z9HB9oO5gYl8dfMwLzHrp+8c43ZrLTPo2uvND8LB91qVB09YaVZE2Ywy0HnGyZGXilTYd/xv+M+xN1zCOqIZswh53r2jdyktE0XGChiDRsBUPZlR58r06GN9WRTvd5CAmOZobYfQ/fSxv3cGnjHs6b8z5DbYE1XcZYy4h2dR+pZOxl2cAodwdvnFjm3b64PoMj0RMDjnk150nfhhDohQunRs+fTh7gl7UP8U7MVTzBDz3VM2W0yN4olq7F+vZefhm6nBZ0ZR2cMn5DTHuX0Nu3di1niM0JerjiRA3OtASub97Gj4u3s/XIYL6+ItAk9OvS99k+aDZCODkaE/xZefvfZTH1cGELKN33+vU/5Mntirhf/VxgzsWv95z6zPxU6Q9x/xK4T5Kkf6M4UlvPlL0dwDj4fA7ueYHE9ijinYMhAu+0yB/ZY2cDGFJzBLM50Lmxjfkk/Vv5EWkcEpGbFLGNXg9lfoVC/KfgyIKkBt/Htz4ynG1hYTymGUpUbRY6t2DxjtALBoyotGI3aAizywyrs1OwNpHIUVZM51px6SQyJ8f47ul301WDr+Suii8ZWmtny+Dp/FF3Pz8u/Zzx43YjybB21B2k5X2G3a7l6IcTkDSCxMsLqZ4SS65lFEtzDrDm0Hic6NALN2FocZrsmAbbMVdGcAQYnd7MJM8qTYeiRlDsiKBNF82J9hRiG1sob01mQ+IltOhjvXHCWDz//FiV65tWf11yCV25UXqM7FrFDt/VpCYh89zCP1F3PI6crHTsiSYmz81ld9U53D35g6Da4xVfDaGpLoa08aVERFmJS9lMTscFTJt3EL3eQQfH+FFBGKZaF3WmgxjjK7120c6R+5KwbGZN2M+o+jxG78lSQgK6cKt4nx22iyHMjMnZwbj8E4wrPoEjNokZ5SOCjr8l91ssziHs19dgxM4SdjKPI2w4OpuUqGbKq6YwuaoYoQ3tbHu6+h9sqg42CXby5pHHmeIMTonfdfBHQftGWytJsp/+AhZfHv1VwPaFFfvQVbi5S+Mrkuc/s9j89UomRbzOE/pfsdR2CLRwV+s69DorB90TuFa/kwu1mdyHUq9ePt006W545ZsHGao149YFfreudK/C7ant87DlDzzcuWKlHtImXsvbDf8TcPzvS97i9yVvUVOYxJqIl8lJ1bJw7ws4zreRW7WEIdGK+WrJYQvUfMC/zruBCVIpX1v/QNmuZLY0P8Ty88eT7g7M4p5QVczJ5BSERkNT9OlF/Z0KfQmF/AglOjRekqQK4HFQ1FQI8RrwFUoYZAHKz/2uM9VZALThfKiJ4HK9nnlOpftatxz0lwSIR4i1L0cVFgXt68R/zWf/GGmHFFjzIMXhomnyBegXv8LRLyYwuNHBiGobXP1P2P4XiBsFxcpU9fDUGGZn+hJpxtxkB1s7XSsGzDjWyoZI37xx9Ijd/CtvNo4iIwuuzuamD95hY9x8vi2YzIymw9RqBQ8l/YtJJ04wPKKSKtMQPjpxE5MacmjWx7JHmkZWquKWH91RhFkXiVvSkmyrpXhEChZdBLjgep0S8bFJnkn9IF/tkRZDLKuG9V4Lw+S2YtMGxukbZDsaIVji2MuQsFb+Ja71tkl+C2Vk/2ksBw9egyzbiZnZTvpMn+lg2pBm/FYp48TKNK574k5c5ufRnnARPlwZmUZEtDB7TmBhsJiwJpzRGkaYskjObgDP8/PCyFoM5OFwhjE8sgDtICtLXz0SUtxvP7CTGFsky3gJszDwL08g2GALhBpfLmQ/f9bPAfTe5LRozNzUto29bTMYydKlo4MAACAASURBVEEKSWG2u3un/sXs7rYtlLD3xMdZvz2l43tCj4vZZBEthy5zMTZM8T8kOFto1NbQ6eW5XbeZ252bvSEzd7mVUWynuL8m7uRe6V0ARHTv4ZJd/R4Aw1yVTKm4kogRc6jmf4PO0UqhC7bNjPw85H6A5LRa0lt2IpdcRIJrMGkbtzN80lQiPJ/85Fo7EzevR5/YyqV6ZTnOke4a7ox+kJsOX0SEdov3Wt/87Uk+nzuVi7L2MuvAEUaNdhAYbdH/9CVa5tZe2gXwH1sNeEzsGG4fdw+a0gqE59uikYPFWxm5e/ooaTDafdPT8I4O5h48SOFF6aRtUX4sW6ZKVI7Wc8fnDoyyL93Tf2xhFn4FlyZdz+TyA7x12bsATL725P9r77zj46jOhf2cmS3qzZJl2ZItufcqgyu2sTHGBhwwl95CgCQkIRCTXEK4hJBGviQ3CQm5kIT0ECAQSggJxXRCL24Y44ptbNlyt1W3nO+PmdXO7s5qV9JqZ1ec5/eztTt7Zuad9s4573kL+nfMSZOpF9M+8UIO3/85KsxIuCPFblbN7cfClw6wv3A07s+9wLOP3MEtu7dxzq6j1Dx6lE+Xvsj/bZvF0bGGz7Nsk1z5tJknpxbuWAO5g1toMZXoewPDYQfrW8ayvn94yPx8eWyewq35YXPPgegiEKFzFdWbcgfb8Wkeytv2c8Lht9icP4yPcwbyqYbHOs5Pqe8wQkhyylvwN7mQQuBvMh7SqhP3UTn5AE1bSvgUbmYUryW4s401TfUd7okF+SMoKBjN0aOrY+SZOOEu3nr7HAAObS5i0LCZVA+6kL17H8dVYB+3MGDAWTQ0PMyuQbnmqXuPwtVhZ7xcDtLuzzfll6DD4UsCtPl0vP4Ae/p7WXNkKqe2vcq01g/5B0ZJtAJ3OydWLGP78fXMH3Aev2FV7M6BKaznDSbHVMYawH7qeJdTbMoqZgOJ0mJrwvDdHsA+RhG/81TXaMROBISOFFDIMX4qP8dXuBNXlcAu5X/7pDI8q40ec/sJ5XhfjxyRVMgm8muMif14OYe6w6KSn7N93ySm1T9DvmhCtIerRp1Z9iPun/IVbpI3xKRgzNUiayiPnb6GBvysyRlNzf4DbK51MyZlUtqTdRGq/1q7l/97ZChD9RJOM282v4idgHHjZyN1Rj3RYz48Ph/n3Xc/zXm5uHzGQ/cftnCsTjB5m+SuZTrLAqOBNehBuGlFLfvyCzjZUjj60fbp1LGJ7/kvZO26iZTlX8KOG//JtScP5+Dho5zbPIjtrmoKPtjLdfe9x9HW5QwX09gsqzGj5anLa2Nboxe+vQoYB4zjd2aH97bWz0EVjDn8ETd6bQccHYq9M/L8zTS7wsEew49vZmDbHnbmVLMrdxA+Mzvf9ENvMeb4RtxBH5jzimcFX2LDoVLGHduAJ8om6y5sZ1jjtohl3uI2Bs/bzbZnqmlpzKOwvIK29l2UjjxM5fggORVGigWPp5kzap7Ed8SFr0aPiXeaPOkeNm68lX79TuL9DeGapJrmpV+/eRw48AKzl32f/pULSYS0O3GWd5YQbeRpHo4FPB3LA8WS+45fyOK8B1jrncOxOBO3k0reYFTRNNqD8T108mnGhb+j5myIuginsr6HT7SgA2fTefnEssPGfeWvyqdthNEhqqCRP3IuB0Q5X83/Fl9u/Rt/C8yjwXSNlfmWKlglXgIDctEbjA7b2dqL/O/Bu2gQt7NFVJKX4uywpQVvky+M3rouwtfUre3gSLF9Bk1NRM51eEQjS3iBYXzEW2d+AX+/+AXpU0XWKffptYbXxdZAGQfMJzMgYlNrFtBEw9tFrN9di9Zk6bU3hz/7Nfj+eRpCGulWtzePpFXfwItV03g3cA4chafbZvJt8S0AXm4fwSq+a6zcFuC4ORq441kjMOrP2g+N3ujvwxduswxHHgJsa07s7RFtixzQ2sDglp00eiqYdHQNPs3NPyuNib5Ld/6ZvEAL+z39GNC2t1Mr5qSj6wgiaNO8HHUVUtkeW2RC04JMORppMsgf0Exx3TF2v1rJsNM/omlvLg1vGhGPbUe8bHqsjv51w9i3bQtnXH8jOaV+3ll9Jrm5tYTSgWsWW1deeRt6nj/CruF2lzJ+vOHHH63cQwl03F4vuumnLzo50sLCcezdG2misfasNCnBUwAt4Qe1faTk6PO1PHHsRnaKo4wiXAGpwJJzvtD1dwr5O77gYCDsDWRlEht6lJWyJzzBApZa6timk65a0Gtbd4MuQBoFUfScIOXsZyzbWel+kIX6u3yq/dsxG9elH9/EUkSzH+2oj/m6MeLbrb/PS66PuNGmWlhPGG4JhKp0R277ov6fj27eKSPYzl9du5gguuct2BWyTrm37tjIkJadfJRbw3tCMg/wW8wyLzGdubzJJuooaDqI1hSbQnTTJbMZ8adXWFV8IU07Cgi0DEHozbzmL+asMyJzhje0jepIABZNoix4Z43I472DsO1Ac0f7ivb9BIVGVesecoJtaDJIXqCZF/vNobR/JQO2vEBFrQYBEBpcs+1u23qh9YffRkhJsd8YCle17Y1p0yHnlP3sfdewgGpIcoOt5LYbrpiDZjXw8X/C+Xrsag2PWG7YwHe/WknhoGYKBzUzdsItNG7fytpnjSFJxeA69m3bQk5BIUUlpQihM2zYSgL+JjZ8cCOatXQQlhqkNsycsYpXXzN66KG6okBkpqeoGdYTT/w3x49/QFHhRHJza9iz50GamjZRWDieY8fW0To2CCGvvojateHPAkGbaZGz+oPfYFMsxK3t4Br+YCt/SZKplxPxay7gLP5NOclXznqf4Ukp959wOdfz+07bbGYIw/mIhzmVsxL0xqHzwul2zDpi1NsNBFxwPBdyzJB889Lmm+7CgX7eqDeHaevUzc6d+eb2aHsYbrlu92uncl6w+0XYQ0xyxc8uWqJ1fcL6JF5jv3tZT0RKiqxT7ttXv8Op+57iV0M+wxaXofQCQufPnMVRCthHOc8xkyA6c6Xh8yuBR4bN5cSG9yluO866+VO4rWgzBz6e3OHxIf3FMfsao22hUIZv2PNy32XKzqfZeizSXh1AQydIu3DjkT4aPf0obz+A2Gb4hfqFjiaDndoCz9nzCANyhtFwfAtHRXhmz06xA8w89EYypwuAqhMaqTqhkffujrXyVUw4FKnc9VgZt/yzBleeYco6vK2AkrrjTDnVKDoy678uQtN13Lm5jJt3MqUDjGH2yQs+7Fjf5SrizUPXkkM4SCp6QsxKXl4t7vYgPo+GprnRzYRSmnBb2gzj0KGwXbMgfwQF+eHzNuPEfyOl5ODBF3lv9RX4hsmwcgdGDCrjnc0N6LpxXB5RYEYqGCQTIBeqvPUsMzmZVxO07jqNlPEPFvFp/tZpO2vWS3+Sj/RREhfIOGTOQB8k9tlIFcOCGwkGNXAXoJuF1gOmJg/VAD6hZht18j0eZCkAubRyHDdu2gkQdqkcrz2NNfnVh6Gw7074NtdSz2rm8xrrGMV0kote7wkVHGRbXu8HRGadcvfk5uENtjNv/4u4zeR17W4PmzsySkDQnGiVSL47/RJ2FlbyUdEAfjXB9Gf9B8BXY7Y9fsRWbnrszwwbsofC0lZyacMXDJt89P27YxQ7hBVwyEZd0R5Zm9JlZwO2oWGbYQqwKr6cslZaDyYuS9YdimuNh+nMlTfBE4aZRwQln7v7T9z1WaNAxsW3/4yn7r6DfR8asm1/qpqv3Bc2eRSUhc9HzbiJtvvp3/9USktn0MLLtBzw4i32ESfOrQN30IWPICAYPeo2CvJHUVoaztEyYvhNlJcvYPXqz8TdhhACKWNfjgLJ0vo6Zp+2nPdWG2kJZs5cxVMv/y/9fV13UcvrJOtlMvzdfRJ5PhGTtdOHG828t7ZRwzaqmdB+MnvajzKx4Ocd7e7hfL6K4U/dhoedVHGEQsbzIfGxVy6v+Wez3jWI4+TTTC4fMpSdDGQDwxhjMVWliu8GbqSVHJoJdLxCQv7jQ7UGrtEf42sf3M+q2uEdyv1/uJl3ZT3/ZD4BNAJxgp6khFeC/8VsLf7LMYDO60zldaaiEUiLcgc43BBrSk41WafcSyqN3BgjmrZwVNQCkb0sKeHDgOEp8vsTO3cD81T8G2/58x3rTRx2BrOmfMAxn4cc2hECGloLGJxvujAG01XGL7yfVCj2tiNuU6FGUneqUQRjxAmzCKVgLykrJ6+4hDFzFzBu3kIq64Zxye2GLfzH550OCITNHEciygbVsuG1d6gcNpQ2mfgBmrzgOfbtewKv13iD19VFFgfRdS/l/eYn3E6epeBEB1LiysmnpCSsyF055RxyHaXMZ3hERacJ7oxke8udEWw/n3vEcBrdHm40C69IRIdy9yNYfeBqHqYAHcFPLcrdmi9fonEPRmoIq3J/m/FMi3BDCSv3l6lHJ8hM3iEg3ey0RAZvMnu/D7GU5TzVo0LydlTvbmVzVQ5FbS2424OUHPUxP/hex+9fc98PQGlbeIKmml1Us4snzNxDcf3lJfxLq+YY8zhEMfubpjHL+wTTXPbzIUF0msglv4cv60SMYxOPtqieewxFFYbizg22ovmPgBv8Umd7sJQC0c7j7WM7XT9n4H3IQA563kfoOUasVf2GUnZWNnOJGMHbBwby/L5hDC04QE3eEV7YN5SVY4zC1SKOiST19Hw/cy+8nJfu/T0AH/69Dt2b3OhBExpCCJZ+cWXMb6d9cWXHhGZXGT7sa+TlDqam5tM8+9yIDqUVj9zcaoYMuTrhdmec+BTBYPxc/nl5dbjdZfh84bwpwuWGIbPRgAGVn6LYUqjiiGbY6Q7b1OmNx/PMZA7Jez+8yURypZ83CsZxRZPRqxztr2ZpXtTkXFBDNyei/TIHj6+EA14fZUG4s+FhLi3/OYWuZ/Hj4h7O7Sg27gscY8z7ayOyDkafbylFh35fxRxOMssXHpb23lh+XDzE0hjlHl3YvQVvRG2FoyKPIhm/MPzw7c1sriph0QdhL6yJxJbS04X9CAxiI11DaAhcwGtmFheP7uLNwExGudZSEKdq1iZqmcwGnmF2p/EGPSUYUMo9BpelEnmbbvRqPwhU8rx/uG37P9bsp+Lnt/PoiYI/L3BFRH+GGLetkPHbinjhzb/SFDBSim493o+txyNNMJ3ZiVOJrknsajdPWryM1U/9s+N7blExQyZM5oNXIofzU5cuJ6/IGOROXLiENav+TaA9trddmfMtKkdEFuburE8+du6C5A8iCpcrn8GDDRPKwKpz0QLJTxJ2Rn6cFLCR+y6IVO6a2wyHh3Hjfhxe7gsnVztECT/ksx3mjs7oas/9/cAFjPPX0C7D5i27nHRl+6fzWEEVF+Y9yoviErMhrPP4GeNz8a32K+nnGo1ERPS2B75+iGk7dvLzsRfwJf2vsRs2txNCInidKZRxmLvkyUyzzx5iy3+oZxLvU8kBXmEa7zCBL1kmat/yjOLktnc73Ub13sTpBjrrWCU7ymrP2U8Dgh/xWQbSYJZ0jOQfLOJlpnOQUupZQwnHjOhuS8zCH1jBZNazl4qIF1s8YkdN0J6fXGerJ3Q9Q4/D6O5wz7Gy3cilET35VShamajv5pcb7+Hf642KQprEVrFD2K2uKRBbnd2KkZ+6e0xYuCRm2aDR9qMMIaPkFIJlX/4aiz4T2bO75td/Ydm1X2Xl/Y8z6ZSlHctn/deFjJk7n+lnruCki69g5f2Pc+2fHorZj0vWUFGxKHJXabglxoz5PkL2vs0xhIw6n1q87J7+yEIQTeTxCvHLz3WH1U3LWOibwABZYpErOnTM4NcFkte0gdS23ktr8UyzJWx3S3490M+DGtzbZlRlOhr0ovnyyT86lJ2Fxuh2a6CGp5jLrzkfV0xvIfI6t5LDw5zGIRE70ZqoU7MNw63vGAUcs0zUBqXggJZ4Mnb81sQlEGsPxm/jMaNP7yEyEZr9WTXYzQAOURKzPICL/fQjiMb/YbxQ/84S/mHJ6X2QEuNcdTLR3EwOP+YqbuV6/sEifs5l/CtUBgoYGx3o0QtknXJ3ucMKeEyTYVMMlWqr1g4zy7WNFd61THV/TNPSEzve6anodXtsPEmSZeQJkQUbasZNxJNjPwS21la85jf38uU//Z3Rs04C4PTrbrRdZ8oSw3vl0h/+Am9ePrrLzUkXfRpvnqGw3B4vX/r9A0w4eTHTlxsRn3nFsbnp9d4fLQJd94nuGTYmiWQISJ5hTqdN/oiRTyd0ZxxusK8UFOLlY1dyd+k/+SjQjLTclNEP4iNlV3BMC/8uzWjbLWalrzvON2qMtktj1LAjWMq/2saxxT+AB4bNQd7xa47rufyHej6miqeZy1rCGTBfxr5jYefRJQW86D2Rt3zVfMP3aT4K9o9u0fG/zzKKOSIKOajFrxDVFYpa7d0sp4hN1Ik97NQq2EkVT1gUaFdol7Hmxja83Mr1bGAEbzORjWZCv1B93t1EukG/ro3hdSYD8CIncoxQ2UDBAcp4g3ASskoZmeStN8g+5e4JK3eXaYer1o4wUd/NIs8mRrrCfqftw6oJPcfxlHvtns4fRisiSa8X23W1yJ7qOTd/m7ySWOV6xle+jkaoxqogt7AIl2W0MmrmHE5Yfg7nfOM7Eev1qx7Myvsfp2JwbVwZPLl5LP7stcw9/1LOu/V2Bo+P9W7RtPT0qENK6TfSPq98Khkw4FMR3+On5Y+8Sdqk36hZ2glbiYxQ/OjtCv5GeBS1gUhz4fbSdTxS+QSfG39Dx/4EsS+7d0tP4/tnh43m7oocfljSwj6Xsc5JIyu4edkY7v38XO5rnczb/mr2uCT/zvMR1DTGLp7DydPCrq9HKeQhi1yjRoTNWRdffDErVhgBWadOqeWCi2OTkF04ZxTrAlX8JXAK89qNdMP7zCxyYdkFL7SHj/d3nEvA9Hx5R+t+5s8wknxaGP/+UQqP+dnYcikPe7/JVG0zfqEDgjeYYmmffBdibzCxHvgby/gRV3WMCI5QxBpGW6QzRnuArcnHei8FUlByMBFZZ3P35oWHfUd9hv29XctjqvZxTFsZlGw3OxmbB8Ze6HFbi5j+QayCjYcWCHYpsHnA8JEUlpWz6Y3/4PKG5woKSsvQNJ35l15JefVgRs6cy6+/8GnOXHmT4bnyzP2dbnfuhZd3QYpYhKZRPSamHC4Aupae971fc1Pbei94MHME9h5D665jyODPwouGD368F330HdLc3Ehu3iCebb8Jd97TzCW2TGAIaXrJi4CMMAe84J/PGFc4qvHfo2MDogD2+CQ1loHc6ZNrmD5hMGdOGkiLL0C+x8XGhmM8tjpsD79yrtGT/OD2T1F7Y3gu5htLDaU+tKqMjWHHkwgKPGHlMny4oZCHDRtGnjnSOxb0UqgZE6PXXHMNL+yInLQe13oPF+S8R3vA2zFRI4GtwTJOb/sOu2U/fv1fo3jtOSM3UpPoudeXiwDztdVU7m+PyM4KlkLrFobV9ufD7cnlTU/Gbu/HxfGO3njsehLBK9TTRC5rEmSOWZ0/m0Wdtug5WdtzDwiND45W8MjuK3iXcfaNBayv1fjC53VeGRd7qCN22ldbj0sSrpBFFeEha05BIYuu+gLzLr6CgSNHc8kP7ohom5NfQP0ZZ1NUXsHK+x83FLspt/EnXa6XVtLzvhfJmkZSsS+h4XKFOwVanBdYdMCY39/CyI1l7POPYlUC80zHvoKSDy35+Jtawykenj8SLrT8o72NWF8n29sjw9qnDzXs5vleF+UFXnI9OndcMIWt31vKpu92Ptq5ZKYxmpg8eTITJkxg9OjRMW1yvbHXOaTYAdwW+1x5eTlnTIqsJNRELn9urWfz8Ym8zHS2MJg1jOHRL8xhnRzKQYqZNm0arqBxruN5tHQFN36OxAm+sssvNXpULfXT42d9tDJpUPcCtaKVeACdt5nYqb0f0mOWzDrlHsIIdBB85BuM3any7N9DhfmANJYY+SvGbiskt1Wnu4Z4ERUQU9CvPKbNoiu/wLD6E9F0F3MvuIy8omLqzzgbIQQ5BYatLuFendDpJiJNPfc0VCyMixbnuutRy30BHS3ooyG5GuwA3HaeiPCeeWaU0dtvDpSwvsWYVC8LBDi1uSXi5X1m2TejhLQ3j2mawK3Hv0bbb19GjttY1+v1smLFCmbPNvLDL168uKPdysUjYVT8EPj+hUYnKj8/H03TyHHrzBwa6T3WiptvXj+bvf5h/IkVtJLDpJqoSUrzEAMpUGdfcT3I11332v5m57Hk8XrJzT2W1LYLc7pnjtxEHT8rMOawEnXGGoPhF1P0vdYbZJ1ZBmD49JlsfOt1QEcGDmOn3L2NH6NpGlWNORzL96MH4YQNZdTszcMVFFQc7nq5tlBOkuHTZ7D5zdcYOGI0Hx54OapJkE999X9sVw+5J84575JOd+MJ2ueeTgfpsrk7+QaLN2oIWny/AfZTjBb0QYJRxoT3NmDOo7G1UtDgX8XMg+X04xhBc9WgxWxw/8cNIUnMv5IcLUoJeVMX8l9TU8Ott95qfDGzkxbmuOG8P4NNBC+EPYyuvjoca3D3pdP4z+b9vLRpP3953Sh0M25wCf/vO58Pbx9YMbWaU8aGRrDGMaai536V64m4v/ltAutcrp6rt+HDh7N5s30ishMnzuH1NS9H3Mr/bBvDH06fxSP/eJj2nLD77R9aDa+mX7hW4HZ5yclNnP6hp2Rlz335Dd8g3x0E4QXZHpUIKpJT36zkrBcH4vYbh+r2d1OxQ4dZZsHln+W8W29nyeev6/hp6RdXUjelnurRcUxEGCallfc/zvgFsVWKrLQGcng6MJVvcnH35OwBms3wtjfIDRo201At2XQi4gwbXDGBMi60oM/WPfT3nMNvMYqYuI6u4S2fMfkZEJLD3sPcUJXHp6sqO9S3Vbm5Yl5sNvJ0YwQ1oaumBU0D3V4BhpS79VwV5bhZMr6K7541wXadED8+dxJLxhuR5DJKudvZxlOB3Xa9uZFm14EDBzJ37tyO7wsWLKCuro7CwkLOOOMM2+1Onjy54/OYMZEmmCVnLWT58jxcltumURYwdnYNyy8KW9QvuOACpDEjw+NyDnuoJLew9/L1hMjKnjsYboltfg+B9rUgp3Y6zteDghwziCdoeWaacvysmtbIma9UJbVPl8fLl373IABFZjGMyacu48i+vYyZu4AxPQjysSIRXOW7gQq9+yXSuouwSwvZC4RuPHcXMwmmgnh+7jPmzGftf1YhZZDi4mLW6V42B9rw2uQu2W6WDS4+fJibL3Px2zuaeaL/LAJ1uyLa6aYib9f8tLiMCFJdAvkV9p2SU74N+7ruA7321sV4XKl7MRcVFXH8+HF0PTllPHDgQHbvjg1+Ch2hz+xZP8McTvBvYXvbNKbmP5oqcW177jk5/Rg//hd4vR4OH2rh9NMNd2Gfz8c777zDvHnzmDdvXsx6kdswJoILCgo477zzIkYoQgimTPkarz51Vcx6o8eOAjO0ZNSoUWCmIQ5FABeWRbuTpp6sVu7hgUfkQyIOGg/Y9sfDZa7cfuOB7m/ptft0ydQPYwMZ4qFpsQ/jwiu6ls85KRy0R4s4PblUE29SMx3EywV/9imzOfuU2fj9foQQvHrvb9gc3MX4vUb77/Albiac08XdVsziJx/gt1N0vnZZG1pwJ9EXb2vLRP5WeIx7i7xsHvoNwLxrp10Ob9lUK5p9beyyJCjMSa1r3YUXXsi2bdvIz0/OfHD55ZfT1habBsLac7+V6wEonva/DBhaDA+nTrkHbHruuttF//6nURmlR5csWcKSJbFBhSFWrlzJsWPH2LFjB56QA0cg7Abdv39/zj333Jj1rJOo8V6KDcEiHmqbwPzBIzs9nlSQtcrdrUlCMfq5OzbRUmt4BOzyv8novQCCoy+u7Wg/Z03s5GdJk5uSps4fiqmnnQnbjRBjLU1at7NCFL1NulwhHXx/xY9QNQnZav9y4RVcd/NPCUlrnbQTATfFhyYahT+A/cX223z06LOsqLYZGWqZ/egVFBQwYULn5hcrHo+nQxFGYHNaZp5l+sI/bPzZIyuoErGFY7qCX+hcc801fPzxx/DoTwDQupgH6brrrqO5uZnCwkIKCwsZOHAge/YY+af8fmMe7JZbbkEIEWnaS2L6aPN3T8MflDQcaWX+j55n2/4m5o9KvF5PyEqbO0CrX+DTjOxtrpbjaM3GkHfErkJb5agHu6dOrPbxovLeH0pBt8ytKSNdPXfdQXeZZC1Pbt3NuKJYN0KA3OaBSb2Ecy2ml9JAgJUHDlEcDMZ6w1z1HJz+k+QEyybMww85h9hZou6Wl3FXzuU92o0fnf79+zNlSjiIyeXqmn2/pKSEgQMjXT5zcw0zis9nuExpmhZ3zqYzXLrhcVRbns/Mof1wdeLxlCqS2oMQYokQYqMQYrMQIib+XQgxRAixSgixRgjxvBCi2m47qeRQsxabg4Wu+4a7NT8X18YmNsotNNwWrcFHExYujmnXG+hOepIkaWPtKZ258/U2ehdyLMwcWYkWDL/wGjBGgHogh6o9iQt0WGOaZ7W0cvnRULWhqPM8aCrUX5G0XE4zf1QFX1xgn6zPSrS/d7wiKL6cnnVj7bxxXN3MYGrF6vsfD+vTevms2oTt/3r1DC6ZYV97NZUkfMKEkbz7TuA0YCxwgRAiOjHFj4A/SiknArcB30+1oNHMG9YS6cYV9irrErqQ5LkiJ/Uun/QxupnDRtddnN72Hb7quxqXp3eKZsTipFkms80FqaAr+ehdXhcCnfKG2Wj+HH7FRXyba9GCHsZsjE1N+zt9MA/u2tPxvcK01c5oaeGW/WHXOCrDEcLOBKv1jN9/+gRuODWxQrbmqlnrH8CqYKw3mRCC0vLkUyzbYZcvSE+BW2/I1BTqwdvu26JGbz0zfHwLFixg0qRJdqukhWSe5BOAzVLKrQBCiPuA5cD7ljZjwZwtgeeA5MLCekB9dQvPff1IgwAAHxpJREFUbQtP9nj37aK1qg6trWuJ9l0iSIFFuZ9atZF+uV6mLDmDl+79PTkFBayTQ1kXGMqVaVJ8Tgb4pMvmrqXJK8cOvQsvT91jXHOBjkB09DyF1NlVX8O7WqR3TD254AtHPeVLydptOyI3+oU3oWIkks7TTPQphORtfw25bht/9GA3XZOt25CxsSGpmrS/7LLLKC2Nn6YkXj3gRJ44vU0y2moQsNPyfRdwYlSb1Ril4H8GnAUUCiH6SSkP0FvIIAHdhWb6S7uaj1OwZW2ClWJZMXgdQsD5Q96jxNNKvssHchDTz1zB9DNXRNjXUhEUkQxOToToaTpGJyeNXV1wGdTtJgkBITX+deV4nv4o3EvPd+cnV2elwvSUyL4Oe5cpbTPMULnB+PEMp85bzgcNydcEtsOulGWqAvLq6hLVYnXSPSA+ydzldpJH35Y3APOEEO8C8zBKEce8SoUQVwsh3hJCvNXY2LPZcWSQ7bXx7VY7+jdz38KdMcs/6t/EwQLjhXDNyFcp9xoVWQblHTMUu7ntmBlx0jfZ6GTPPX0vMAd77l3IJR+MaBuWeeKiwTz90dMd36+ZdA0vnZe4cEMEmakTUkooX49uV33GpH7RSILBIH/mLNpE12vYgn3u9nSl0pBRfzOFZI5+F5gRGwbVEFmqRUq5W0p5tpRyCvANc1lMdn0p5a+klPVSyvoKs1xet5EBWnJjxW8sbuOvi3bybH0jrd7YbtSQffmUHfcwY/F8cvU4Yf5jzrRd7O7i7Ht3cbLnrnWjPmr39uOgcu/ChGrAGvVmeXrbio9HtPM0H8K9dx2Z94g7S/TZuHahferfQCDAZmrZMvTOhNtMNsq1O14t3SJDX9LJ6JE3gRFCiDohhAc4H3jM2kAIUS5ER9z614HfplZMG2SQNVXPxyz265I2T+dj42EF+5kxr5MKO6d+z3ax5u79HMwAMksUX49w8BhdXRiu+wOxj4gW8PLbDfdELBOv/xJ+Nb/TVBhx6c46WYK1R7399mV8fr59WcSOICFXGW3uztWSXS/dPgtjeu6xTL16CZW7lNIPfBF4EtgAPCClXC+EuE0IEerizgc2CiE+BCqB7/aSvGH6j8Wvx4au+9yRiv35yY28NjbspfDlUS+zvHoD+t8ui7/tOOYXLV3Rm2lMhxuNrnVeajBViHS9RGzQulAooW6ytdqOIXPRoXEc9kfWgO04GnPS/bkdu/j3TrPGwKB6+407eJ3TTaIjDSl3TWi2ivol19SOz5mnTIXl/8whKW0lpXwCeCJq2S2Wzw8CD6ZWtAScdRfcNytm8avjDkZ83z7QsKnPeN+w5blCKQTaEtdtjCZdk41O3iYiTf7n6UpQZrtvV/Lnt7h/2M/ZulZAC5v0XBJObDEnDPOMtLjlAUsn4+IH4Qe1NlvPPDXlFEGzPnE8D5cdmrWQe7LXL10990xT6wZZG6GK1wgyWjPUUNI7+jfz7ojDtOTYT9zMv/QqFp3/Xz3apZYuTxIH3QS74knSE9KWWdgGXXT3OoaDKfxa2N3x3fZ+jG03v/ttvEISudA6OYOeIYQiQwvi+LtbfcntXonKLBNL1kesvDP6MO+MTlxKa9qy5bDvA8Nps5uIdPXc05DIPx56Gmo7ArZpdNNFKkYN7S5rPIVFiWx6ymaHlnM65ys93nc2kihQa+nSpdTX11NVNZA2m8Bfq7E1Gfc9RR9Q7l2ih0Uwkk1/2lOc9AHX0nWMDnbde3qMrw/5B8e9hs19Qmsb7Hkn/GP78dgVrD33ReGKS3kBI4uiW3ah1FMfxe12M2jQoPgNLGmX7V8UNq6QaXuMjB1l2gsme80y3cEm0KErpCtNrZOjdLcrPROqXie9Zbqr3M0J0I+LPuxY9LN9ScRrxHmR5QeM3n9uMDZVbt8hNdfZmjsm+XQNaTLLZJpWN/nkKPfWoxDsmXJPV1ItJ5wo7vKfTkCKtBnDncyn0tPgFuvIypPMYcR5kWXqRFxqSc11tuaOsauX6qR+lRnqLfOJUO617T544Qc9V+7pCvBx4Da53X8hw9r+kkbTk3NmGdHtHEGhhzh8fdx23bbcUlj2Y/jcK/DZxFGrGdrxSwmd3ck/4HP8Pz6X1HasL8J7iC2UYb+nTFO36eUTYXP/5d5GGHA8ObPMf38U9yeRJvc9J3t0Wpq8ZZxNsdC9F0tD/i7KfUW0a+H4Clvl7i2E6Vcm3N4nQfV09uJqwT7Tor03TJhGYgvv5GjpL9cYIlNHYJ8I5e6WEgL+5HruufHL7qUqEVEi0hY2bYOmp8lbxoFj/FTbbczR1jK/m9fxxYHPUNpeyjHPsY5lPXmAQgorM1VDakjVqMTqCokk5qT1020mstNlc3fQdbkzstos89LY5OpN+gWw+90eT6imy//cQRdwSJdZxgHl/p4czi8CZ+F2d00lt87azlHvAaQmOZgTGSRnexRd1mh91zCTqqtsPUPe1uQqoqXrHsvUQOOsVu4l7nwWNjXH/X1EqZGkqDAoYd96eD9OQd4rnkxqf+m6WdLjr2KPnoLqNcnhZJm9ru3bP7qRe6feFrM8Nxgnh5GTEVoZRniysWcvMKvpQxYnjmtxgkwzz2S1cieBDfxHJ/2Id6beYtSsBGjcaN9w8IzkdudgyHy60NJUkMTlYG4Zl7vnxSE65fy/9O72P4FYzTKf+cxnYn+3fWFnlrJNN9mtrYRGu81FnVo6irWXrWVoyVDc1gvcwyCmpCsr9xAnC2Sna17B6zVHCA5YJPQuThpfOOZCw+MqWSpjS8kpus8bvpEELbaPysrKTlpbSc/zmqmvkOxW7ggabDwfqnIss+kBy0PZQ1fItOUgl85dFr0LSbV6gsvBSaiuxivUFNbwj4/3JG7YbTJVPfScbk0aRzVuHjifmpyC5HbkAGXtRwHwyswKRstubxkhOGDzoEbUyAxalPuunpXySmQGShXCwdwyottJtbpGXmHPCiL3hFSNTmKuUs0MmHJxSratCJOXl0tTa6LayHbpB9I00jbvBC3DJsazu+cuNNpsLqDLqoQDqfN/TVf6gaCDPrtpM8v0tt27E/QevqTvbNgH2Cj32dfC1Eu6scXMUgqZhq67E4621BmMJct77hqtNsp9Runo8JdAJ3b2ggFQatZhnfffUJegWnnaqnY5d6tqaTKXOOnLr3en6Mo5v+XB9iPsfuYmxrUZL98Y74g0ZdTMRrpyT0e31HQvGp333J30VHHyee2MLFfugoCNkjitwqx8s28DPPvt+OvfYPGeWXBT3Ga3zfoee5oqEZzWXUmzhnRlhdTSVBTEdt+ebhzj+BWMAkY9dj1HzBFcjH9zZwFg162D1jgufH3X5E4q+tTC7UlinsR5Beu8BJFkuXLXmNDaxtqcqCF+0A9/Ogu2PJuS3QwqaGBQQUPaYua9aSp1Z0falK6ek5792O26J+a1wTPRd74G2DzMnSn3khoi68x/Uui8nnEyuD25sQVP8vpB84HwdwcjifIwRnI5NmU/nSTrlfvdDftocLk4u7oqvPzwjpQp9ojdpUm553cUBUl/XyBt5e/czin3HqVYmHwRmqncYybYu2uWybQuXyrp5JG59tp4EeaRK7k8+bHxF599CX4yNqn99DZVehv4oMKbWcFV2T2hiqBQSkb4fIxps7w1H73GOZFSQMcA1IGHvqfpcJPej4MPY0/rxLrMZGFd6rl3Jk+PpMleysrKKCsrS9gur6Aw9poVD+JAXri+rZM2d72lGAB3c6ljMtiR1F0uhFgihNgohNgshLjR5vfBQojnhBDvCiHWCCGWpl5UG3qYKyZTcXmde+emu3qNE4juTKiG0N0dD01PlXtf7rD3hJgJVZcH3Z3g3Dp4MqXfMAtr7XkJWqaXhFpEGEnM7wROA8YCFwghxkY1uxl4QEo5BTgf+GWqBbXF4uaYqTPW3SFtphEb0pWz3tHMlz3x5R+/An3Wl5lXNYtfLLwTSoaEf9OdzAqUoaTgMmtuNyLhi/OTOv6JTzJa5ARgs5Ryq5SyHbgPWB7VRgKhqJRiYHfqROwEGzfHX5g+yHFZdGuviJJKnAxiStczImXPJ9q6S/eLdQC6G7H4Nn6x+G7mVs+F69ZAsTlRmqa8PNlEt26nqMlRzZPTs2uWNjLrBZPMGRsE7LR83wWcGNXmVuApIcSXgHxgkd2GhBBXA1cDDB48uKuyxmIToNQvkEBpzLkeJl0A/swKFY7EQZNF2tICCJwaS/daucQuB4Cl//jfmFKC2x9kSpr325ORtebx2npxyTifX60vwdMumdbtPfYNkum52z3t0VfqAuD3UspqYCnwJ2GTQlFK+SspZb2Usr6ioqLr0kYT6KaCLrQEL2UgwuVc9Ga6zCW6ZlThCfrT/yLrvVKCXTsWJ17hH/kG05BfmLb9BVPgoqh5c9ESFm4P76c5z8XhkvQFlKXLCaGrJCPVLiIddKuJNbt8BngAQEr5KpADNrWwUo0lKVjo0vYFy7t02ZcfSwfpSmusC8MVMuhL/4Oh9WRCtTO6+GJ0wsNj/bqFvP6aXQ3S3mGPMNTAi9qEpNeJfoZFTh66nXJ30uXKSobm70/myXoTGCGEqBNCeDAmTB+LarMDWAgghBiDodwbUymoLRETqn2ILHYTzAZSHoW76FYQOuQnVyEoljR2SeIVGOklDosialvv5T96V9IgRz4AujsHLWHdW+ceGtnhhJBZXcuEXRgppV8I8UXgSQwX7N9KKdcLIW4D3pJSPgasBH4thLge4wgvl9KucnCqCV/QEtPW7sqwE6zIPFI+OplwjvGvyxj3bDrVUsGm1XEKW/QO3dlTTM9dd6MnNFU699y3ae0QhLZgZs3jJTU+lVI+ATwRtewWy+f3gdmpFS0Jpl8JT/8PAN9rPMATBXmM7kpRBUUMaSslaOaNL28/kKBl6slUG2k6EMFAWl8mIiW+kDp6Am8Zq4lr9erFeD3NLDy557tOBp+7HdqgJdCUnh0mSTb4F8XHEw4aKA0GueioXQX07MNJH/B0UeTVOHv3I1S0NwKfT+u+M0259+Wx5hCP8fKe3W9N0usE/B7AYnLVXEkEnoXP4tEjyVZqShWZ+bxm1l2eoWx/ZhD7N5SkbX+i2LlCFml7r+T1Y1DbHrw4EGWcIbVwQ+6BmakaUsMAl497Fl/L6MKPkl4nEIiyrwuB2yZCNZ4rpMIgM+7yDOfwliJ2vViVuGGKcPZhT1fSeifTD2TGbe9kPpR0IbrhChm9hhAiCZu7IprsNst0laHznZYgKQryqoB2ylsPOi1Kn6RHuWVSSvr7m2dcfyPHD6Z/nqNrxJ4XlzuRzV0RTabc5b3D0Pmw9fnw93P/5JAgXaO8pB8rdv+SysAh4Cvp3Xm6OpPpcKZSxDByxpz07rA7Ec9BImwKQgg028Rh9tvObWqhLbrGQy8SNPPeOJhS3pbMGJ/2hJroTAgW9KgL7E1fZF5PGdjWgFt2UiIwxXhJr5dRx6SxA+aZXgti6iIZpgt6hZycAgBcXThaEYx88QshEtb2tfYVXDvWk//hO8kL2UOaPYZeCWaYOs2Mu7wnXPx3aGqEHa/CI1FeF9ZMcif/T+ZEtGUgy73rOSa9kKZSgnlFRg7snPyCtOzPSqZ5y/RlSsUADgAlR5Kfs7Ib0+nuBOkHnAz8c27XnZL9d7m3AMrqYPKFsb9ZU7CedEP6ZMpC8oSPSi19rqQ5hYZH0NyLLk/bPsNkxuO4QdSyPjiEh12znBal1/Aer6Hmzf8m74P53d6GECJxUXNHrXxxirc4TPb33DsjN33ui32GNOk9l9vNyvsfT8/OTM7Z/RA782oQ2pK07jceo3Ib+HHpEhbmfeC0KL1GwN1GYeMYdrd2v+ylpmm4EqQfkA6q1iLNiDTWMky9Z3/PvTNGnOq0BFlIZvRqe4Oqtn2ccOhtMuUYdQEXjXmIfFer06L0GgFvOw9/9DP2Nq9Peh0ZNTNp9NwT+bmHvz1w8i4emZuekhJg7SEr5Z51zFhxPsX90xf1pqYGPhlkliroHVyV/WkPtqL365f0OtH+/0KIhJPgAUvJzeacAIcL0+cgUJxryJZpyrRvm2VS9PjMPvdiZp97cUq2lanoUiMggpnSqe1dMuXtGTR6owF/3w3QcVcZE6nu6uourGWj3F2Zq6q00AR9ptxXJpl7xlJBMACXPgoth52WJOM5p30Gh0Wz02Kkhwx5CAP7h7Blcxti+3inRel1upJATGLk+j9IMUE0yjQtYfUsazeuWA9yJJBp/ej007eVuwxmTVSq0xTKXAqlc0VCPom0FQXYvXsMpWRWqlin0YQxknmIpXzMAL4pRLh3nAQ3DWjFn06bV4ba1/r46y1Dz7pCAfQrMpSYph1zWJLMwtbmbpN+ICIvvWUS1qtBvgPFkTIt6LpvKfcvr478LtNbdSZVOHqTZIbF4hNBKDI/JTnP+xSGWhKEn1+R0BVSncNo+pZyL62Fix8Kf89S5e4kn4Rc8hlDnllm2FKXoK8xYNgIAMbNX5j0OqEwfqvfeOL0A87dt+4hkwDY5ZnhmAx29C3lDjB8EYw72/icaeOkJCk03cbyilUQVp8mt9T46853Vo5epKi8Pyvvf5xh0zrJARWDMP8PP78JOx1B55R7sHAQdzXcxzbPIsdksKPvKXcIF2PI0p57KNTa5UmQT0OR1ZQPqQWgX81gZwXJOGyUewb33AECeDNu1KuUewYSCtjIKcieLJaKruPNNcwxTiRPy2SkDNnck++5Z8IYPRNksJKUK6QQYgnwM0AHfiOlvD3q958AC8yveUB/KaVzNoUO5Z5ppzs5ivtXsvCKzzP8hJlOi6JQpB0pYm3udpk8I57uTEumngEkVO5CCB24EzgF2AW8KYR4TEr5fqiNlPJ6S/svAVN6QdbkCb3ls7TnDjD51GVOi6BQOILIzYNAF23uSrnHkIxZ5gRgs5Ryq5SyHbgPWN5J+wuAv6ZCuG6T5WYZRe9wUe27zKnY5rQYikRoRp8zwhXSNojJotAdnFAN4bwEkSRjlhkE7LR83wXYTn0LIYYAdYBtfk8hxNXA1QCDB/fiJNJJX4XGjTDm9N7bhyLrGJB7nAG56ctZr+ge0sYVMvE6zqnWXKM0AQUV+x2TwY5keu52Zy3eWT8feFBKS4o260pS/kpKWS+lrK+oqEhWxq5TVgdXrQq7mikUYMRBZBBVZlKtUaNGOSxJhiFiJ1TtkJaevZOzmSUDNGoXf4ua+vecE8KGZHruu4Aay/dqIF6y5POBL/RUKIWiV7hyFRzMHLNMZWUlN998M64MznjoDLGukLZYfvY1O5gXSUpySnYhtCHOyWBDMj33N4ERQog6IYQHQ4E/Ft1ICDEKKAVeTa2ICkWKyC+HmulOSxGBUuyxSGH4tCdS7u8OsNRlzYgJ1UyQIUxC5S6l9ANfBJ4ENgAPSCnXCyFuE0KcaWl6AXCflFnqf6hQKDKCHfkrANjJwE7btbjd7KrKSYdISZJZqi+pboOU8gngiahlt0R9vzV1YikUik8qh7xTuLPhYY4NeDHpdTKrz5wZ9M0IVYVCkb0kGcYf2U/OBPWeCTKEUQY/hUKRUcw9bwTuHJ3RA+ey/0BmuRdmE0q5KxSKjCK/2MvCS8cAY5JfSSgjRDTqjCgUiqzEagTRNZVBNRql3BUKRXZiMbo7q9wzy0smhFLuCoUi6xHKLBODOiMKhSLrceFARewOMstLJoRS7gqFIuvRHO25K7OMQqFQpIz21nB0qrRNCZxmVJk9hUKh6Dm73h/H4UMDgCSKeaSDDMu8opS7QqHISoJBnZZWI5m6fTGPTzbqjCgUiqwllDnSSW8Zr7cSgKLiSY7JYIeKUFUoFNmJBM0s2CHzCx0To6BgFCec8E8K8kc4JoMdSrkrFIqsRTeVe47X2dS/hQWjHd2/HUq5KxSKLEV29NzRwxGqM2c8Q7vvoEMyZQ5KuSsAqLyhHn9Dk9NiKBRdQsco1+wp6texLC+vjjzqnBIpY1DKXQGAuzwXd7mDdSgVim4Q6rnnFPdL0PKTh/KWUSgUWUvI5o6m+qnRKOWuUCiyFjubu8IgKeUuhFgihNgohNgshLgxTptzhRDvCyHWCyHuTa2YCoVCEUvI5q567rEkPCNCCB24EzgF2AW8KYR4TEr5vqXNCODrwGwp5SEhRP/eElihUChCaMosE5dkeu4nAJullFullO3AfcDyqDZXAXdKKQ8BSCn3pVZMhUKhiEaEbe7KLBNDMsp9ELDT8n2XuczKSGCkEOIVIcRrQoglqRJQoVAobNE0/snJ7Gofgq8484KInCYZ5W6Xbi06/ZkLGAHMBy4AfiOEKInZkBBXCyHeEkK81djY2FVZFQqFogMhNHYzgPsOXwVuZyNUM5FklPsuoMbyvRrYbdPmUSmlT0q5DdiIoewjkFL+SkpZL6Wsr6io6K7MCoVCEUEmZPzNNJJR7m8CI4QQdUIID3A+8FhUm0eABQBCiHIMM83WVAqqUCgU8RCa0u7RJFTuUko/8EXgSWAD8ICUcr0Q4jYhxJlmsyeBA0KI94HngK9KKQ/0ltAKhUJhJSOKdWQYSfkPSSmfAJ6IWnaL5bMEvmL+UygUil6nSHPRDGiucmWWsUFFqCoUiqwkp6rK+CBcyixjg1LuCoUiK3FbnDJUzz0WpdwVCkXWo3rusSjlrlAosh41oRqLUu4KhSJraAvm2S5XPfdYVLYdhUKRFfxx3920y1yutPlNddxjUcpdoVBkBceC/XF5wsYGqylGmWViUcpdoVBkBRd/eyaeXL3ju8tlqi8p8OYpVRaNOiMKhSIrKK6IrPG7bNkySkpKGV09hZJKe1v8Jxml3BUKRVaSn5/P4sWnOC1GxqK8ZRQKhaIPopS7QqFQ9EGUclcoFIo+iFLuCoVC0QdRyl2hUCj6IEq5KxQKRR9EKXeFQqHogyjlrlAoFH0QYVTIc2DHQjQCH3Vz9XJgfwrFSTfZLH82yw7ZLX82yw5K/lQxREpZkaiRY8q9Jwgh3pJS1jstR3fJZvmzWXbIbvmzWXZQ8qcbZZZRKBSKPohS7gqFQtEHyVbl/iunBegh2Sx/NssO2S1/NssOSv60kpU2d4VCoVB0Trb23BUKhULRCVmn3IUQS4QQG4UQm4UQNzotTzRCiBohxHNCiA1CiPVCiC+by8uEEE8LITaZf0vN5UIIcYd5PGuEEFOdPQIDIYQuhHhXCPG4+b1OCPG6Kf/9QgiPudxrft9s/l7rsNwlQogHhRAfmNdgZjadeyHE9eZ9s04I8VchRE4mn3shxG+FEPuEEOssy7p8voUQl5ntNwkhLnNQ9h+a984aIcTDQogSy29fN2XfKIQ41bI8M3WSlDJr/gE6sAUYCniA1cBYp+WKkrEKmGp+LgQ+BMYC/w+40Vx+I/AD8/NS4F+AAGYArzt9DKZcXwHuBR43vz8AnG9+vgv4vPn5GuAu8/P5wP0Oy/0H4ErzswcoyZZzDwwCtgG5lnN+eSafe+AkYCqwzrKsS+cbKAO2mn9Lzc+lDsm+GHCZn39gkX2sqW+8QJ2ph/RM1kmOC9DFizETeNLy/evA152WK4HMjwKnABuBKnNZFbDR/Hw3cIGlfUc7B2WuBlYBJwOPmw/jfstN33EdgCeBmeZnl9lOOCR3kakcRdTyrDj3pnLfaSo5l3nuT830cw/URinILp1v4ALgbsvyiHbplD3qt7OAv5ifI3RN6Nxnsk7KNrNM6OYPsctclpGYw+QpwOtApZRyD4D5t7/ZLBOP6afA14Cg+b0fcFhK6Te/W2XskN/8/YjZ3gmGAo3A70yT0m+EEPlkybmXUn4M/AjYAezBOJdvkx3n3kpXz3dGXQcLV2CMNCD7ZM865S5slmWku48QogB4CLhOSnm0s6Y2yxw7JiHE6cA+KeXb1sU2TWUSv6UbF8Yw+/+klFOAJgyzQDwySXZM2/RyjGH/QCAfOM2maSae+2SIJ2/GHYcQ4huAH/hLaJFNs4yUPUS2KfddQI3lezWw2yFZ4iKEcGMo9r9IKf9uLt4rhKgyf68C9pnLM+2YZgNnCiG2A/dhmGZ+CpQIIUIF1a0ydshv/l4MHEynwBZ2AbuklK+b3x/EUPbZcu4XAduklI1SSh/wd2AW2XHurXT1fGfUdTAndE8HLpKmrYUskd1Ktin3N4ERpveAB2MS6TGHZYpACCGAe4ANUsr/tfz0GBDyArgMwxYfWn6p6UkwAzgSGtI6gZTy61LKaillLcb5fVZKeRHwHHCO2Sxa/tBxnWO2d6TnIqVsAHYKIUaZixYC75Ml5x7DHDNDCJFn3kch+TP+3EfR1fP9JLBYCFFqjl4Wm8vSjhBiCfDfwJlSymbLT48B55seSnXACOANMlknOW3078YEyFIMD5QtwDeclsdGvjkYw7I1wHvmv6UYttBVwCbzb5nZXgB3msezFqh3+hgsxzKfsLfMUIybeTPwN8BrLs8xv282fx/qsMyTgbfM8/8IhvdF1px74FvAB8A64E8Y3hkZe+6Bv2LMD/gwerGf6c75xrBvbzb/fdpB2Tdj2NBDz+5dlvbfMGXfCJxmWZ6ROklFqCoUCkUfJNvMMgqFQqFIAqXcFQqFog+ilLtCoVD0QZRyVygUij6IUu4KhULRB1HKXaFQKPogSrkrFApFH0Qpd4VCoeiD/H8RP7sDtPOanAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(fluxes.T)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 90, "metadata": {}, "outputs": [], "source": [ "pca = PCA(n_components=4)" ] }, { "cell_type": "code", "execution_count": 91, "metadata": {}, "outputs": [], "source": [ "commons = pca.fit(fluxes.T)\n", "transform = pca.transform(fluxes.T)" ] }, { "cell_type": "code", "execution_count": 92, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([-0.18309613, -0.17031454, -0.18458051, ..., 0.09266309,\n", " 0.10965268, 0.0886739 ])" ] }, "execution_count": 92, "metadata": {}, "output_type": "execute_result" } ], "source": [ "transform[:,0]" ] }, { "cell_type": "code", "execution_count": 93, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 93, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD8CAYAAABzTgP2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnXecFFXW93+nexJBYMiZISuCiI5ERUVWQV1RV3fFhOlV13WDGxSzqz6Kq4/u45rXgGt2dXdhFUElqAgCg4Ek4AAjGQaGDBP7vn90VXeFW6m7uqt7+nw/n4Hqireqbt1zz7nnnkNCCDAMwzCMSijoAjAMwzCZBQsGhmEYRgcLBoZhGEYHCwaGYRhGBwsGhmEYRgcLBoZhGEYHCwaGYRhGBwsGhmEYRgcLBoZhGEZHXtAFSIS2bduKkpKSoIvBMAyTVSxdunSXEKKd035ZKRhKSkpQVlYWdDEYhmGyCiL60c1+bEpiGIZhdLBgYBiGYXSwYGAYhmF0sGBgGIZhdLBgYBiGYXSwYGAYhmF0sGBgGIZhdLBgYJJmxvJtqDpUG3QxGIbxCRYMTFLsPFCNm974Gje+tjToojAM4xMsGJikqK2PAAC27D0ScEkYhvELFgwMwzCMDl8EAxGNI6I1RFRORJMl20cT0ddEVE9EFxm2NRDRt8rfdD/KwzAMwyRO0kH0iCgM4GkAPwGwGcASIpouhFil2W0jgKsA/FFyiiNCiOOTLQcTDEIEXQKGYfzGj+iqQwGUCyHWAwARvQ1gAoCYYBBCVCjbIj5cj2EYhkkhfpiSugDYpPm9WVnnliIiKiOir4jofB/Kw6QRoqBLwDCM3/ihMciaBi8Ghu5CiK1E1AvAHCJaLoRYZ7oI0fUArgeA7t27J1ZShmEYxhE/NIbNALppfncFsNXtwUKIrcr/6wHMAzDEYr8XhBClQojSdu0cExAxDMMwCeKHYFgCoC8R9SSiAgCXAHDlXURExURUqCy3BTAKmrEJJvPhwWeGaXwkLRiEEPUAbgYwC8D3AN4VQqwkovuJ6DwAIKKTiGgzgIsBPE9EK5XDjwFQRkTfAZgLYIrBm4lhGIZJM77kfBZCzAAww7DuHs3yEkRNTMbjFgAY5EcZGIZhGH/gmc8MwzCMDhYMTFKwuyrDND5YMDAMwzA6WDAwDMMwOlgwMEnB7qoM0/hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcAwDMPoYMHAMAzD6GDBwDAMw+hgwcBkHe8s2YiLn1sQdDEYptGSF3QBGMYrt72/POgiMEyjhjUGhmEYRgcLBiYphAi6BAzD+A0LBoZhGEaHL4KBiMYR0RoiKieiyZLto4noayKqJ6KLDNsmEdEPyt8kP8rDpA8BVhkYprGRtGAgojCApwGMBzAAwEQiGmDYbSOAqwC8aTi2NYB7AQwDMBTAvURUnGyZcpnKAzVYvnlf0MVgGCaL8UNjGAqgXAixXghRC+BtABO0OwghKoQQywBEDMeeBeATIUSVEGIPgE8AjPOhTDnLuL9+jp8+NT/oYjAMk8X4IRi6ANik+b1ZWZfqYxkJuw/VpvV6PPjMMI0PPwQDSda5bS5cH0tE1xNRGRGVVVZWui6cW6rrGvDAB6twqKbe93MzDMNkE34Ihs0Auml+dwWw1e9jhRAvCCFKhRCl7dq1S6igdry28Ee8NH8Dnvtsne/nbsywwsAwjQ8/BMMSAH2JqCcRFQC4BMB0l8fOAnAmERUrg85nKuvSTm1DdPijIcJNHcMwuU3SgkEIUQ/gZkQb9O8BvCuEWElE9xPReQBARCcR0WYAFwN4nohWKsdWAXgAUeGyBMD9yrq0E1EEAsmMW4wlIsBBhiCvzTCNGV9iJQkhZgCYYVh3j2Z5CaJmItmxLwN42Y9yJIPaxIR8lAxCCDz28RpcclJ3dGvd1LfzMlGEYEHOMKmAZz4rRISqMfjX0qyrPIin567DDa8t9e2cmUaQfXbWFxgmNbBgUFCHFvzsgKrnrGswTt9g/IBNSQyTGlgwAPjX15vx5OwfAHBjk03wm2KY1MCCAcBTc8tjy0/OKcd3m/b6ev7G3IAFKUcjLMQZJiWwYJAwv3yXL+fhcdHUwnKBYVIDCwYJj85aE3QRsghunRmmscGCAUh525YL4xZBuI3mwGNlmEBgwZBC1MayMbdfQTbOnAuCYVIDC4aUwqMMqYQ1BoZJDSwYGF8IopFmucAwqYEFA7iBSYZAZz6zysAwKYEFA7iB8YNABp/Tf0mGyQlYMKSDRtyCBTr4zJFGGCYlsGCQcGKPYt3vW9/7DiWTP/R8Ho78mVrYK4lhUgMLBpg79MZQC++WbU5fYbKMIBtntgAyTGpgwSCBk7hlB/yaGCY1sGCQEPFZMnADlhrcOA3MXLEN3/ocFJFhGju+ZHBrbPgVtTPIIQYhhK9Jh6yvk/JLWF/bxT43vv41AKBiyjmpLQzDNCJYY4C5cWsMpqRcsL/nwj0yTBCwYIB5ANVvU1IQpOsOOFYSwzQ+WDBIsDIlJToRjifQpQZ+rAyTGlgwSGiwaHGySZFIlzBqzO6qby3eiJLJH2LfkbrUXohhMgwWDBKsGhyvg9KBxhEK8NrpItVC6ZUvNwAAtu07ktLrMEymwYIBwKYq/YdvJQAaPKoMgdrf03TtTL7HgzX1vlyHOHw6k2OwYJBgJQASdWMNou3MhYFZpzu86NkFyZ1fc4EF5bvYpMTkDDktGMoqqqQxkHbsr8b+anMj4FVjyA2DTnA4jaOs3n7Al+scrKnHpS8uwg2vlcXW7a+uw33TV6K6rsGXazBMJpHTgmHqggrp+roGgdF/mWtaH/EYzTOTzSyNAfUeIxGBWSu3+z7grp6ttj764tfuOBjb9uSnP2Dqggq8vXijr9dkmEwgpwVDXYN1S7/3sFljyKbB51xi6oIK3PDaUlyYpOnIiCpoVLNcSDPUUK9ojw38knOWFVv2YfnmfUEXIyXkuGDw9lVbubE6EUjaywwYfBZCYOqXG3DIp0Fgq2urXkPfbPQ3JpIwLcQlA4dUz24Wb6hC+c7kTI3n/m0+fvrUfJ9KlFnkuGDwZhvyOiM612cFz1m9E/f9dxUe/uj7lJy/8mBNSs4bw/AIQzbCYM32A1ixpXH2HhsjP39+IcY+/nnQxchYfBEMRDSOiNYQUTkRTZZsLySid5Tti4ioRFlfQkRHiOhb5e85P8rjFtV27BavGkMmNM6pxu4eD9VGB2b3SMxyfnDTG0tN6+obIqiptx8QnrlimydXVrU/ELJRE8766+c492+Ns/fI5B5JCwYiCgN4GsB4AAMATCSiAYbdrgWwRwjRB8ATAB7RbFsnhDhe+bsx2fJ4wWuP06tTUq4PPqfa2rJjv/n9nfPkfPS/a6blMeU7D+DG17/Gbe8vi607WFMv1QbVNWqHQCYXONwJo2XAPTPxxCdrUe/RGpFp+KExDAVQLoRYL4SoBfA2gAmGfSYAeFVZfg/AGZSOmNAObN3rbUZrosH1gtAcMqG5SnUZ1BqkbZvX7LC3Gx+siWoTm6oO4+fPLcQ/FlZg4L2z8Pgnay2PaVDc0bQaA096Y4w0RAQO1zbg/2b/gD//d1XQxUkKPwRDFwCbNL83K+uk+wgh6gHsA9BG2daTiL4hos+I6BQfyuMar/MSZq3c7mn/YDWGNMVKcnGZVDWhiZz3rUVx99LFFVW4Z9pKAMB7S83pW9VnqDopBN+VYbxSXdeAd8s2peV70I5Z/nfZ1pRfL5X4IRhkn4vxLVjtsw1AdyHEEAC/B/AmEbWQXoToeiIqI6KyysrKpAqs4tUr6cEPv8cnq3a4VhMDDTAX2JXTj9t73bm/Gu+URfswywxuhvWGSSq7DtagYvfh6DYHweDViYFJH098uha3vrcMH6/akdLrrNiyDzV18XqQrBzavq86UDOlH4JhM4Bumt9dARjFZWwfIsoD0BJAlRCiRgixGwCEEEsBrAPQT3YRIcQLQohSIURpu3btfCh2Yvy/f5Thuc/Wudo3F8zPdreYafZ3O+cBoyNCxa5DseV6iSlJZfX2AxwqI4OpPBAdh9qfwne0vvIgzv3bfNz/gT/moxVb9mH4w7PxZoCTJ/0QDEsA9CWinkRUAOASANMN+0wHMElZvgjAHCGEIKJ2yuA1iKgXgL4A1vtQppSyxePYRGOex9BYqDeYFQvzwrHlW975FoBcMLy3dLN0MiSTGXgZC9p7uNayM7P3cC0Wb6iSbqtWNIXZq+NaiRAC32/bj9veW+Z5bHJdZXSG/Vfr5ddLB0nnfBZC1BPRzQBmAQgDeFkIsZKI7gdQJoSYDuAlAK8RUTmAKkSFBwCMBnA/EdUDaABwoxAibU+juGl+ylwpAydtE9ysL5Rq/wIighBCOlbkde6ESTDkx/tM6ib1bozxtVI1gY9JH5UHanDS/3yKW8b2w2/H9jVtv/ylRVixZT/WP3Q2QoYJLQV50d/aDoIQwPWvlWFT1REM790a63Yewh/P6p/am/CRpAUDAAghZgCYYVh3j2a5GsDFkuPeB/C+H2VIBNkYw4hebbBw/W7b49z2xnNpglsQrpwE4PFP1kpjXj3/mTfFM6zcwKqt+3H2k1/g1H5mcyWPMWQvTjVx35FaAMC0b7dIBcOKLfsBRMPihAxaiLFToV5PHZu65Z3vAMCzYMj2MYaspVbyQd95zjGx5UdnrcbgP39s2se1YGjE2c1i13FxvVRpDvURgZfmb5Buyw+br2lnVjhS14AD1XUx76TP1podHKzuQ1uP3lzEQfVSTUNE4N2yTa6cQNxWvYJw1HR4uNZ+cqRsnKpe0sEUQtg6t0QiAkt/3KNbd9d/lkujPQdBzgqG6Iuzr1hPz10nHVj02uCzvT91WH3IzQv1yrAQwnGm+5KKKlTbzJq2aoi0DcAd/15uew0med5eshG3vrcspimWTP4QkzUTFqU4fINqg3+4Vm8WXL55n24mvSzCsiy4poDZ003Lqwsr8LNnF2Demp2xda9/pe9UHKgOzkSZs4KhPiJS3mAHa0rKbYy9uIXrdmP0o+ZQ6lqumVqG2d9buzXWWAiWOo+hVZjkUBvMnQfiM9/fXrLJandXqJMY1YHkhet2Y/z/fYGfPjVf12DLhIDUlCTkmoRK+c7oAPOmqsOW+3y2thJHHDSYVJGzgsFqcpsb1dO9KSk4GvMEt44tihz3OWAYEL70xUWuzi0Ls6FilZTH2DNcX3lQuh/jD/nhaLNlp/GrnkBu657auKvvcuLfv8L326LjCg9o3FB37K/GnkO1uOCZL2PrZG1JRGKR0H6TYWUAW3as1mRp1GDShS+Dz9lIomk6vRDk4FG6r2wnUP0sS3VdA7bvr/bxjF6uLW+Iag09wy9+2IVe7Zqno0g5SYEyfmQlGHYfrMGJD36KR342KLbOyfyr9u6dPEvHPv6ZaR8rzcCoSfS8fQYuHdYdD10wKOb6nKn5PHJWY7Bqs5sWOMvKTHmXlQdqMPn9ZY7RRIGoH/Y7S1IxMJrepyFT29NFTX0DHpu1xrR+tdKzVElHpyOXyVM1hnr5c1ZNTP/78VpHbzkhBE544BO8tvBHV9eWVT+rMQaZNqA6J6gaQ6Lx11JNzgoGq4+3Z9tmaNOswJdrpPqVP/DBKry9ZBNmrjDHcDLe3h/e/Q63vb88ph6nEz9NSUHOMu7fsQWemltuWv/MPP1M+Az91hsNMVOSxeCu+m1babHa9yMEUHWoNhYqJRFknRU7R4en55bHBYNDJyKoWKM5KRgqdh3CoPvMbqgqPx3cGQV51o8mU+Yx2IajMGzdpYQYtxpA9UJ9QwSrt0cFjH0Gt6QvZWLUlDn+n9QlJ5UUu9ov00KBNDbyY6Yk+XNuiI0vxBtV7SvR9uQTzcqov563b2rqgoqY0JJd36jRBEFOCoZPbTxPVOwir7p3V03tS7WtNBab/KhoUz5ajXF//UIXT8iORDs9izdUYeeBaizeUJURs4v/4drckNuCoa4h4jkJlhdiGoPFNWKCgfTC4cfdh1Ay+UOd1vyWD/GIvM5vFCKuUfxl5ho8aBNjKSjtM2cHn+0IEdmH5Pb4slI+A1jS8pqu6KNK+vXG6MScXZpER6m4xZ8/vxDNC/NwsKYeY4/pgBcnlfp/kRSQ6xOhR06Zg8oDNaiYck5Kzm/nlbTvcB2eVUx7xho/Y3nU5KoNsa6GXU8GrxoDEfDKlxWx3y/O34C7zo3nNouwxhAMsob0mlE9sfqBccp2f64TpCnJz2PSeT4VtUelpuAMYmwkUWQaw+tf/YiF6+xDrTQWKg+kNhe3GqpI9pzvnrYiFmJb+50LOI89JIpXhwindkErZ/wwdSVCbgoGybrC/BCK8sOW2xMhXa9UmuxCyPfxo56prnbaU/n9sRkn9gghbHtP3Vs3TfhaVwzvkfCxMmTlvOs/KzDx71/5ep1cR9YeH6jWOydo62WqPICswrJY4WRq1G4PypSUk4JBhrZdc2rkXI8wpPql2lq79Bvj95R8obQpNVN1j4fr9GMKAvbX6tCiMOFrtWqan/CxKgVhczRWJjWoz9f4mL/btBdz18RjXGm/YyFS916+2bjX0/5OMZ5e1piZgnJnZcEgwclFzKvdL9WvVu6rbdjHYn1C11PO5vQcvMaUOlLbgM17oiECDtUYNQb759i2eSF+Orizp+up5IWS/wy0QftCPmtPjB613hnr34Snv9T9Nn4XkdhxqSubG/Y7xEDSmk2DKmtOCgZZQ6pdt8gh7LZbUj1wpDa8DRFhSh5kvHKi/tAzV2xDyeQPsVsz0KxKmWgv3vke3V754Y++x8mPzMWugzWm8BMCwlYFDxHh0YuOwytXn4Tbxx/t8opRwiHg09+P9nSMkXyNe7Pds15Qvgs7DwQzc7uxoNaCqMZqXScIZNAaovtu25c9zz8oD7ecFAwytG5tfs2uTdcrfWjG9xg1ZY7OS8gKr2VS1dofdprj//hdZ78s3wUgGljM6BUWNQXYNAIEFOWHcXr/9gldu0/7o6Shut2iT9JiXc5LX1yEi59baLn9w2XbAgucFgT7q+tQMvlDT3b6mMYAYWse0gkFiNhArht39bbN/ZnkmiwsGAJGW4nCDraATBljUM+vBn7bc6hWsy268cwnPsMFz3zp6+zjrYp24mQq8nr/2kFt4wfhNMagfWdeH7t6Xm06z0Q4uU9bAM7uqj/ulkfU/HrjHvzqza/x5/8m70KZLexU4l699MV63PBamau0uer7cuwsAPhIExXAbX+vtj6Cfh2OcrezDxTkhbBiyz7pNq1rbTrJScHg1EjKcvtq8drgJSMgyncexGOz1rgy2Rin+gPA2h0HdYNjVqe5b/pKPP7JWldl2rxH+Xgd7P52/OebLTjfYBPWmQgM+wthH5/e6Z254ZhOyTUG6jyLRF0MDyr35zWneDajPqqt+6oxa+UOPDTDOSWr+nQjwj50PhHpNDm3A7n97voIC9LoWlxbH8G5f5sv3fbMvHWo0nT40kVOCgYZ2mbFSWPYc9jdi/Ijg9uVLy3CU3PLXfmGO5laALmZ4+uNezB1QQWenP2Dp7I53Z3RvqsdVPvdO9/i2016bw61bM/OW2cq566D0Zy87q7lUDAL/n5lKV68MvFJdEX5YYTIedzFaXA66MFRIw0RgQc/WIVt+/wXWMZbddN4aweR7ea3aB+zk3aRyZzwwCdpv2ZOCgbp4KBmnZNg+OKHXfYzo1USrIefr63EPdNWANCkjXTheaQtk9kryfqeLnxmgXVhlPMsUOz/xuvLvrX6hgjmrN6h2/bylxUY/39fYElFlW7fiKTMn36/w7NrYdgHjaFV0wKMHdAhoWN7tW0WLUfIYdY8rL2gYsI7Y+L3Rln64x68OH8D/vDud76f264OA9He9HbDYHHMlASzJ5IO3RgDuxF7IUcFg2SdZtmNu6Fd2j6VROvhlS8vdh2XR3e91MkqPDnHHFXUqgF7dt46XDO1LJY3mYiwUrGhGu3rDUJg694jeO6zdbqz2Q3QytCakgZ3benqGCdZMtmDd1OLJvmxcjiZkqw6HnbCO0jUnnYqQp4b65C2V79t3xH0u+sjDH94ts5LLW5ytC+P9mlGIvZebZmOGoYmXeSkYJChbSTc2Ku9hEdJpgeo1mVZmYzn1TZIpmtqJqUlijGfcXQsIHrCTVVxM8N6JbiezjZqYcpqiAj8+q1vMOWj1djgMiifDG0nfKQyCGzF368sxbRfjcJNp/W23W/iSd1dX19t7AWAr9ZX2e6b5+TckGHtl5fyjPnfeUldqyEiMG/NTiypqMLc1fHJatrc3kJjSrJDaxmoa4gkNVlsxm9OSfhYP7DV6lMAB9GT4GRKAtwNMPrxgauncNOX1Aff0m+LTXAz9tA8fCxvLtqIhy6IZ8VauH43+ku8N9TgZmoDWNcQ0TWcxjIv/TH53pCXeRp1DREM7tbK0W3Ry7w39f5q6yP4btNeLNu8F8d1bSXdt1mh/LMjH4S33xyqqUe1i0RQKusr48J95optaFaYh1P6ttPtU1sfwYZdh9C/41FmU5IArnplCYC4eQ6IvrN9h+sw5n/n4dzjOgFw1n7LNS7WdQ0C//jKmxZ+02m9ce3JPVGUH7Z8Z8nSr0NzrN2Realgc1JjkDUhWjXelWBQYsG/s2Qjlm+Wu5q50RT2HanDfdNXWuYTjvtsy7bJ91X3//2738Z+q41OTX1E5/Xy8EfOXiBWPDtvHX73zrem9WqqQzXT1gfLtmHat1vjBdPwVQKTCW86rTdO7KHPjWA1xjDrd+aJa4ddzhPw4ulk1AKMM7e1tGwiD8FhJbyD5Nh7Z+FqpaH2yo2vf40rXlqMMsO40oMfrsJZf/0cm/ccNtVhbUdlvUaDrK2PYOnGKuw+VItXFTOrF9NQbX3Ek8At7VGMW8cdjTbNC2NCwWvYlcuGOWucmWo+zEnBIMPr2KWqMdz2/nL89Cm5q5mbivjYrDWYuqAC09WG0wI3H4HWf14IgX99vcW0z2/e+gajpsyJfYDTHK6bCOr4i8xkIiB0qUivmVrm+fzNCvNMn5OVLO/f0azRyLQcGV4Eg3HXZoXWcyKOKrLofWZmG+EaK5v/nsP6wHZLKqIaYsWuwyZnBKuB+5r6BlMj6qWh95onXNY5XHTHWE/nsMoRriVTxz3YlKTg9Zt045Xk5pWrCWhkFVEIoRlocz6/m0qmzgWICIEQyBf/fy0791fH5jlIc+EKYM8h/9NzGk1J/75pJGZ/v1O67yDD4LTVU/NiSvLS88tzmGGdoW2FI1YZ1cIWz/HaV5eYMgpamWj//N9V5tnwHsq23aOrrdM7cqJjiyKc0rct3v/afoKaLKJAJpCbgsHgmtoQEZ794F25q7o4n2q/Lcw3fz0NkfgEHjeNvlYNNxbP2HBt31+NrsVNXZnNtNh5ghyqqcfQh2bHfssaCgEkPWFHVobLh+vV9iHdizGkuzkVZ692cbt1shMd7VA9eA7V1Jvu1+oRxoITInqPj3+yFucP6YLe7ZonXI5U8peZq/HfZVsRiQA3j+kTs/0bMT5H9f3J0sxajXl98YPMXdr9N3jEwlRrRbIdpl+e1jujTIJeyUlTkpNrqpvX6W7w2XmfGkXdLJKEY4gIew+M9ZUHTfurjH38M902Yz3/xfPR3ABeA4vaycNXF1bofs9Zbe6xR4RwPUHQCiH099OiKA992jubh16+qhRvXz88vsLh4/cyN8IUyVN5UJe/tAin/GWubptVp0J7jl0Ha/G3OeW4/MVFrsuQbp6Ztw6bqo5gy94juP1fyy3ziXvpfHgxrSyzGNuT4eQpZsRrh8lI9zZNfc/kN1MT3iPV5KRg0KKaILxGH22wUJu1uKni6sdUkGd+FfqEHeazravUu3faCSvj7akD0E6Nn7HXYzd/Y98RZxOREEjKLRUwC6dPfn+q4zEdWxRhzNEd0P6oItfXSabTqGoMslj9e52ek4h3WKwa2yD4cfch286OpWDw8CBdfFZpwU2ZT+/fTrp+xm9Owen92+P4bmavtCcnDokt/7y0q6cy3fj6Uk/7J0POCwapxuCi17K/2rkRfHbuOsd91IFYmWDQxgxyZ96ybkSsbOAhjz0jOxPam4ucE6sLRLOZJYOA0N1Phxb2jf28P56Gmb8z+6EXKea7QsmzB6KdhUQ7jnbPqdzCrqz1SlJ7rF5Mln6yY381ftx9yLCuBlMXVOCNRT/G8mZoqbUQDFZ1TDYI7zV/sp901NQjN9/F3y3Cpwzo3AIA0Kd9c1Pe6/MGd8YNo3sBAHplqIkQyFHBoIukGtMY7I8Z1EU/YHnu3+ZbRkRUWax4XNh92rU2XaSIRjK4UbHrE+hueQ0lYTf71S7InYoa6iMZhIAnb4GSts3Qqqk5jPI1o3ri12P6YNLIEstjE03i4zRLePX2/RBCxJwPgLjWqn3VQWXwGvbQbJz66DzT+s/XVuLOf6/AeU+ZQ1FYdahe+Hw99h02d6RaFJnddoMaeL/z7GNw46m9Yr+dJiECcXdsAFjx57NcXyuUhNBPVtt2iy+CgYjGEdEaIionosmS7YVE9I6yfRERlWi23a6sX0NE7p+uT+RZuUwYkNkcV251b+O0Qv2YZA1/xMErSba/FVbtv9XAYHVdA7bvi3sYqRx338fOBbHBl0l/PrUeRflh/OHM/rFc3zLU9/7rMX1sz2U0RTr1fK+dWoaet8/AsffOig1OyxwgknFnnLliO178Yr1p1u/yzfsSfobqaWQOBFbt3JzVO/HYx2tM61tI5nMEIRhO6dsWV48qQVjTFrjVpJ+57AT8orQbmhW4D9munjoRoX/6Y/M8H5MISQsGIgoDeBrAeAADAEwkogGG3a4FsEcI0QfAEwAeUY4dAOASAMcCGAfgGeV8aUPtGWhNE7LKKZPUby/Z5OoadpVd/fD/+ukPmPjCV4Zt8Ubw+c/X4eRH5theJ5FYNsYPQD3FzW9+jeEPz87IbFfpbDtUt8XrTukl3V5skS/aSXvTTjK4eYVfAAAgAElEQVSsOqSPnCsQrxeJhvAGojbpBz/8Hn3v/Aj3To/mePhq/W789Kn5nhPYq9gpmEaHB+NxxlzHXhrTVHLJSd2RFw4hX/MtWGnSRk3i7EGd8MhFx3kao2zXPDpRzphrvFNLvUn07nONzWj68ENjGAqgXAixXghRC+BtABMM+0wA8Kqy/B6AMyj6JCcAeFsIUSOE2ACgXDlfSpHNcjZmezIiG1h1mwTcrtendiwXb6jCQsMs4MF//hiHlFm6by3epOu9y9RQO9XUquIaO0bqOT61mAOQCQgRtyR5TePplQfPH4i2zQvRrCCMId2jg4kvTYrblh//+fEAzJYtL2YCVWuLjTFoNEW/LEmvKeEg1NDZy23MoBUpMFf8Y+GPGGIIHy3NVe77lZ1Rhb+2k2TllfT5racnfb0rRpTg0YuOw6XDeujW3zymD85T8pafM6gTrhlVguk3j7Ic5E4lfgiGLgC0XefNyjrpPkKIegD7ALRxeSwAgIiuJ6IyIiqrrKyU7ZIQMcHg2xnNqLGDauob0PP2D/EvzaSXRE0FMr/shMYYTBpDhriF2KAto3Gymh88d/kJuGBItBpOOL4Lyu4aqzM5Gnt6QLyR+8c10X6NF+1N1gip9+j3GIM6ZnKktkGX8U+LXY6DZFDHoNTXJxOetR7iMvmFmtI1z4Vg6NyqiefzG3N8hEOEi0u7ma5REA7h7EHRuSB1DREQEY7r2ipp19lE8EMwyEptfONW+7g5NrpSiBeEEKVCiNJ27ZKToA7pGHy3c6qCYdveaggRNRslei31Y5LlBd7pIpmPEaMnSVBeMF4Q0LyvFBR33MBOeOIXx5vWx3v01sd2LY42HIloDOohIvaPf4JabVvUxu/jVTtMPfhZK7fjt29/Y/tI/fw2ZMLT6IKdDlRhqRX+iQgAK9zm+AgRxYRUncbs5tWV3g/8EAybAXTT/O4KwBiAJ7YPEeUBaAmgyuWxvqN9zNJQEz43NvUN0ZwDaryWFk3ibnqm3MYOF1crjCzo3iMzV1sed7hG7jE00OBtlYxNO5X845qhuHJEVPWOCOEpBIVfGBtwGWoj40kwhAjlOw9ixvJtANRsY3C8lh07DbGB1F6n0dlCiGjsqnWVB3HDa0sdY2clWzuiZrLoWRLRcFNBnkRjcHI2SAXhEMXej1Zo+pGEyit+hMRYAqAvEfUEsAXRweRLDftMBzAJwEIAFwGYI4QQRDQdwJtE9DiAzgD6AljsQ5k8k8qGpj4iMHJKfOC4RVE+auobUJgXNgkGpwlNdQ0RFOWHXUcIVSmThLZesWWfKdn4tG+2xMY1MolWTfPRqaXSi9M8snQ2Ldr0qG/9v+EIEWIeTX0Un/Sw0sj84Z/fuXYtFELoBm4Fko+w+qf3lhnKHi3X64bQ0zX1Edz1nxW6evCnf1pnakvWI0wbJsVNsqt0kK80xqrw/MmADrF16YQIsQFwrcYQhCkpacEghKgnopsBzAIQBvCyEGIlEd0PoEwIMR3ASwBeI6JyRDWFS5RjVxLRuwBWAagH8CshRFpbJfUDtE1y4zO7Dtag/10z8cjPBpm0kxqHiIx1DQL//W6r59gvMmQJyO+etjLp86aCvFAoZg7RmpLSqeBoYxmN6N0mtv7N64ahtKQ1AKCJxvX1qbnmrHfFTfNN0UZl9yDTFJZv3oeF63fh+tH2CYYAfcMCxE1JalY9ldqGiCltq13HQBazyAuHa+tj+QcyRmNQHs7QktZoXpgXm4CWbkJEOLGkGD8Z0AG3jesfX5+NggEAhBAzAMwwrLtHs1wN4GKLY/8HwP/4UQ63aDWzuMquEQwJ1Nft+6rRsaW7cAuqHfXD5dtNIuhwnf0ksbqGCH791jfeC5jl5IUp9t6CmvRllQVPmzGuZZN8hMjaBBSWTJgz1TchpD1zNbx7cdMCXFzazbTdDqugcEdqG7A1jS7J2qRMqUgVmgiqdlDcrMDTRDW/CRGhMC9smlGdZKDXxMqS/ktmJrqk9JLtTma+4Q/Pxudr3XlLhTQmCaMp6aqX7ZOiWIUdyGZG92snjSujJS9Euh57ELj5PsMhQrFklnV8u3mdsQ58t3mfdNaxitFMBADTvt2C7zZZu09bCYZhmmi46eDDZdtiy8Z5DUGRbIhtLarzgYyemox0MqysV36HxndDbobdluBH72Xl1v0Y3c/ZYyrawMlnPK/ZccD2WD9MSJnCgE4tsGrbfpzevx0iAvjWpmHLC4V05qPYchrFhDpTN9+hIWlaGMZui+EF2UCi7A5MuQcc1Njfvh3NpFcx5RxU1zVgwTr9nJgArBFS/vVNPHlUOjUVOxINe2Lkw9+cHB8HM7D2wfGO78DS+4g1htSy70gdDlTXGWY5qw20ZscE2xrti9d6DQ226Q1vqvKWQMRNPKIgUZ+Bm4ZIjT9VmBd2bPiipiTVKygYr6S//Ow4TB5/tCmtqBE7LxKZvdhVHCyXHZd9R+qkCYqCGMBMhvEDO6btWk6C3i3Hdm6J1s3k2mJBXsgx/I6VZpDvk+DyQk4Jhp89uwCT31+uWxeLSaMzJZk/QrdmBJVn58Ujq5oaPWW3RAbyDriI6hok6p2qH0FJm6aW+6qRZYvyQ46NY9SUJLleGu1Kxc0KcOOpvR39yu0GC2Uff7KJobSdkAc+WKULzKdCRFkxR0UlFeEgBnRqIV3vNl5aqrGqNrIJlakmM55ImgiR8oFJXoDWK+mRnx1nykblpscVIsLzn63DHf9erjP5GDOZJdM/OWgxHyFTUB+jeo92H/j1o3ujQ4tCnNqvnWPjmBcOxcIZD+neKqlcCanGTmOwSuFqx1WvLLYVnEffPTO2fKimXvpsqg7VovcdM8wbMpQmhsCGA7uYG/V+HbyFrT53cCd0aFFoWp+fIdqUlYCSRQZONTkmGMj0gXVqFfUk0vamerVrjqcuPcF0rPP5gYc/Wo03F23UNf5+hhiY7jABKWg6GzyzVD//343ta9p3QOcWWHTHWLRpXuhovcsLE4b3aoOFt4/BhOPjUVMysQ9s14m4dGh30zpjjCwj89ZUuu7tR4TQBejLBvp3MGffM0a8/eDX8Xwa6vNVv8kHzx+IX56md+E9trNZkISI8MZ1w03rE9EYUmGaswoq2Lww/cEGc0owVNc14ONVO3CrxqtDDVrl5ALppofq1t84maxcH6/akfCxqWZw15b45y9H4v8uiYeTKMoPoWLKOfjd2H66fa83+Iq7MSUBiA3unaoM8nez8QIJCrtOxKmSgGj3uJg74nZCY8QQciUbmHXLaNM6WfKkL249Hd/c/ZOYRqY+5z7tm+O2cfpgin84s5/p+BBF3Yll670w83enYOHkMd4OckHTArkvkMzFOdXklGCo2G3OOqVWLqdQEO40Bs0+CXYo1Kximc6IXm3MK4nQpVUTTDi+S6yhL5TksgbMbreOpiTDx3HtyT2x+M4zMjILlt13nGg/05jW0cr85FeuCrecf3znlJxX1snq1ropipsV4MmJQzC4a0sUKt+K7NuUNaYhInmcNI9v5eiOLdDeIWtgIjQvlAsGP91p3ZIdrVAKUVVCJ43BjWDQqpert9m7nVpR7TDzOQj+pslT6xbV9GEl6GoNPuxDXMxj0EJEnvI3B4XRrz3RgGjGEO9W7X+6x5f9DDan8vmfTtf9njRCH5563MCOmHbzybZjObIMbCEijeYZrzuUIa1gEwtTkptscn6TIY8kONSHbqcxXHhCF10Mfiu0788YeiDTsYv5bpyYM6JXG4w5ur1pP231VRuods3ljbdRYxjZp600L7NKEGEBEkXbA+3eWu+V5ddtNAiBuoYIpn65Qbc+3WHTUzH5qtDQmfjzhIG215bds6xcIYoO5D7xi8GYenU87Usm1Ky2zQstvY+yNex21iDtRcTyr1of9/jPj8eQ7va+60Aw4XFlfPr7U13tV6Cx475ytTw/0qI7zsDALi3xscYO/NwVJ6LIIfvWVUoe5ZYWlf0qSZ5luxnD2YqxSvhVRxoiAq98uQH3/XeVbn26NYZUtFlFFuZH07WV6qsKhsuGxQf2ZeYX9Vu/YEhXdNSYgjLhuy27a6xl4L4gAvrllmCQVRayNiV1LW4SGwRzo84FMXVdRp/27uzubu6pg/IB9dN4jlgdpr39+847FhVTzpHut+r+s0zhvo3Hq7x+7TA88YvBjuXMJLT3YawTfjWktQ0R7D5oTrTjdozhzrOP8acgKajzViYVI/FvN/r7yhElsW2yuk0WY4CZ8dVawxpDipHNIAzbDD5/9qfTser+cQDcmTL+aBOuOBVcabC9AsBd57j/4J0q3JQLB0nXh0g+XDdUiTDqhJUAlZ21Z7tmuGBIV1fnzRTsnqrVQOeZLpO5qIz938/w/OfrTevdTpqcOMzsNpsIbpusLzykxCyQeCTJMJqStNWqbXPzfAVtdSeL5UyExxhSjJ16KdMYwiHK6FAC4yRhA7zMbnWqcJdIfO4BfcM+QeOV8qez+st2N2H1TDP9A/UDq3v0ai5IJFufyoPnD7T0mfeKWy3ZyuNG5ZKTvEWLBbTxspTfmm3tjjILBt2xmuVM0fStMKXfTYO9MKcEg5Hbxh0di2iYqZnL7CiQNCZe7iMcCmHGb07BezeO8HRd7Xek9b12O1HIyptEHvIi+96LFuPAqJXm6Ve8HjeEQ+SbXd2p3/T0pSfgi1tPR76DFjDlZ8clcG19p057S0X5Ydx0mnXeikwXBlqG9myNwV1bxuJHpcPBIKcFw3nHd47PY/AohTOhXkmTyBvuY+wxZu8hlRBFZx+XujQBxY/T+oN7r6RePIysgpJlNJrKYUxSZnXn6YzX42eqSKd3WdwsH91aN3Ul+G4Z2w8Th7rXHGKu5hYN5Q2GhEba3bLJlNS0IA/Tbj45Ni6Xjk5sTgkGYy+J4Fy5rMiEHodWMKi9I6N861rc1HLcIdFbSJV1zVj2l68qtZwNmsno3XYNGoNVBM00agxuBfO9P3UOZCfTWrWoQqgwL4w3rhsWWy8LT//bsX3x8IXuNQf1NuKdOvv70j577VhPEJF6E6FTyyLHyL5+kVOCwQiRZvDZhcaQAbJAh1YwqDZq2X1ce3JP6fGJfhDRwefosX52XoxRbTNB+CaCttiVB2sst2lJp0ui20tdPSpebxbfcYZpe9fiJrhC4gChRSuERmky3V13ck/888YRePt6c+witzx4/iD8vLQrTusf1YpNz9bwW3vf2aQxqFx4Qle8/8uRltEE/CS3BQMoPvjsooHb8LDc/TJdaCND3nn2Mbqoi1aajxDC0p6c6AdBlJqPyShkslYwaJbXV+oz9ljdkl/JYtxQEHbfsBzdMeqmLBvMzQsRivLDeNJmZrzdOzyppDWGy0KrKFwxvIcu7paRji2L8JeLBse8mJxqi05jsHEpZnI8g5tWY/A60p/uqrT6gXEgAvrfFQ2xfMWIHqg6FPdjN6vVUWR3pWZOs7qHT24ZbZspjojiCXlChIlDu6GLD6ER2jYvxCl92+LL8l2IiOxLLuOGTDAlebnWuzeOwK4DNdLOhfp+zhvcGaU9ivHrt77R5XQGrM2Obr62B86Xz3i2wlhGY+hurVei3pTEGMlxjSFeub2m9ky3r0xRflinQuaF4gPA+WHCVaN64sITupgG3FTe/+VIDO0ZHWRu0zyqaVhpEn07HIXjutrHLjp/SBdcOaIHbj2rPx6+8DjcPMYcVtsr4RDhtWuHxQbDs7UnZ1eVrAef0ygYlB72m/9vGH4zpo/tvi2K8i0DFWoFt1XMJKOJrK8y+TIV3mbGJ1iQF8L6h86O5Vax0hiytJqllJwSDLL3H0pw8DlowqG4nb8oL4zmhXl4/OfHx0JQqB+Delsn9iiOhRhX8fpBaOP+FOaFcf+EgSlJIqJqb9mqMdjlrrbWGNJpSopea2TvttIZ6G5xM0ZlFHidUhB0T0X2aEMhin0D2vqkFxLZWc9SSU4JBhMUT4Lh1S0y6KpEGpdRY9AxADHtQDugaxR9Xr+Hf900Eu/e4G3OQyKUKEH7ZLHzs51MGHx2o4kd19VZYBidBWRnTed9WQkq1byqddMN+vvNdHJaMBAIJ3QvxsMXDsKDLu2ZH/z65NjyOYM62eyZOMdY5KY1okYolXkp2FV8NVid11ATbZsXxgROKnlgwkBMvfok9O9ozuyVDdxjk87UqneazrAHejOK/LqyjHtGrJTsazTeTOlMZG8l7xpiITPkpiTGTG4LBopWlolDu+OoIne9U20v9unLTsD6h852LVSM/OMafURTNUzB/ROOdXV8TUwwWL9G7cc7um/UXXDSyBJ8f/843OLi4/eTB84f6EqYNikIx1wQsxGj9jnu2HjoEqsGyW3gOL8xFuem03rjxlN749R+1s9fnSn/VwuPofGD4vdrNXaSTsOtqjFohS+bj+zJbcGQwDHGyUGhEOHy4fa+3FaM7B111TtNyYVw+tHtUTHlHJzkciZy51ZFCJFFjCKl4ms/wB5tmqFiyjk4sUcxmhSE0/5xXDG8B56+7ATnHbMc9b2qhMPOJgynWEJ+oi2DOsZQ0iY6ftS6WQEmjz/adnyntKQ1Kqacg2M7O5ubjIIhiOa4IcvHrIIgp9xV/YiN72fVCocI5f8zHuEQYc7qnVKf7scuHowBFqalpgV5WG8xt0Itp5cx9fm3nY4d+6vdH8BIad+iCKP6tMGX5bsBuAtB0cyDYOjYogjbbd7TeYM7Y/p3WwFEhdRZx3bEvdPleaU7tixCxZRz8OAHq/Di/A2+O2E4zYz2E6vHrN5TNiV7ChrWGLwe4+Ig2XR/+bkIeeEQiAhnHNNB1zg8fOEgvDSpFBed2BUDOrsbc5CX0/2H3rW4KU7skfoxhFxAO8CrNWFYvQ03GsOVI3rgneuHY+HtY6SJjlSMk8KM5z66o7k+qY2mH3JBe450xoCy6ujdfe4ADO/V2nVYeCbHNAYjiVhS0hVXZaJFyOtHLzoO6wyzaWWkImQF4x5tI6XtqVppD25NScMUrdLOi0577V+e1hvb9sa1C6vkSeohfkd0thxUT0G9tPoy+3U4Cm9fn3pvusZETmkMRje9RBp59RR2QkUIgdP6t8Nzl5/o+fxOXFzaDZPHH+24n/pBsl01GLRPfdbK7bHlZoV5eHLiEJx/vH5OiZvBZ62Zx+17PaVvO5xnuJYMu/zJWga7cGPVYnRXjedQSMEEN67qvpGUYCCi1kT0CRH9oPwvDf1HRJOUfX4gokma9fOIaA0Rfav8pdQVxQ/Pj1jFtqnXQgBTrx4qTaSTLiYM6YyrRpbg1rOchQjjP9p2+0B1PU7pGw8gd97gzqY5Glb+/lrXZauw0TJGaMar7LzWjOW1m5G86I4z8JbHoHdGAZbKtjtboqRmA8makiYDmC2EmEJEk5Xft2l3IKLWAO4FUIqoArmUiKYLIdSgKpcJIcqSLIcrTEnGU2RKSkVvyCuFeWHcd547t1fGf4za6WvXDtP9rq7TJ2qwMrlcOqw77v7PCgB6M4/TJLXXrxuGeiUZhBsnC7Ve25mS1Pzfbrj25J54af4G0/o7zzkGh2obbIPnJQprDP6RrClpAoBXleVXAZwv2ecsAJ8IIaoUYfAJgHFJXjchiowaQwLttxtTkjE5CwD85oz0zhlggsWpkaqp1wcptPKYOVhdH1tuqqm/Tp5O4RB5Cs8ccqEJO6Et0t3nDpCOZ/RpfxTevWFESvJssFzwj2QFQwchxDYAUP6XmYK6ANik+b1ZWafyimJGuptS7FjfxBA6IpGevZsCys57/ehenq/FZC9OVVmdnKhiNWQwd83O2PLvf9Ivvr/PY0eXD++BId1bYeIw77mXVSaPPwa92jXDwC7eveh8gSWDbzgKBiL6lIhWSP4muLyG7HWpLedlQohBAE5R/q6wKcf1RFRGRGWVlZUuL63nV6frI0km0jtSVXi7Y2XquJ/pFJnMx6ndrjaENbdKFHXH2fHse1p3Zqvze83frdK+RRH+fdMotD/KvbnIyIk9ijHnD6cFlnWPxxj8w1EwCCHGCiEGSv6mAdhBRJ0AQPl/p+QUmwFouyFdAWxVzr1F+f8AgDcBDDUdHS/HC0KIUiFEabt27uYJGDmlbzuU3TU2fs4EzuGmfT9Fk6lKRRsyplPLxD8+Jjtw8hq6bJh+tnxdg7w2Hm0RL0o2FwFA1saXYjKLZE1J0wGoXkaTAEyT7DMLwJlEVKx4LZ0JYBYR5RFRWwAgonwA5wJYkWR5HGnbvFAXPtoraq/ESkAsvH2MSTMB9Bm6Zt0yGl/dbk6VyDQenAaHxw7ogNUPxIfarOSIlUfRiN7ywdsCFx5IjRVWyv0jWZ1vCoB3iehaABsBXAwARFQK4EYhxHVCiCoiegDAEuWY+5V1zRAVEPkAwgA+BfD3JMvjijeuG4ZZK7d7DrUNwNGO2amlPN68+uGHQ4QWRflo4TJoH5OduJlnoHVR7dGmGR67eDAGdmmBI7UNuOCZBQC8h21JZzTTTCMb5UJ+mCy1xSBJSjAIIXYDMHV9FffT6zS/XwbwsmGfQwD8nwHmgm6tm+K6UxIbDE48TzLhj2f2w5ijOyR2Aiar0FYTK3OQUXhcdGI0DPqB6jrd+qtGluD4bvYZ9VRyOR5QNkZMXXj7GTrPs0whp0NiJEIy6Sb9SH/JZB9eZ8Ab6xjPR3FH9omFqGm7bfPCoIthggWDR7Kx8jHBclSRt88sFbmu7zrnmFiCpsZKFioMGQsLBo9w5WPcoLUa53m0+6eijiVqOs0mctmM5jcsGDwS80pyqTtMvfokrNy6P5VFYjKcsEUWMyu485EY6UyP2thhweARrx/taf3bZ3WaSiZ5vDZYiZiSurRqgi17j3g+rjGRChNcrsKCwSNWYYO7tW6CNs0ybxCJCR6vDVYiDdzHt4w2hdnINVhj8A8WDD7xxa1jgi4Ck0FoQ6Z4zYmh7t61WD4nxsiwnq3RrDAPud4v4dwj/sGCIUE4LgtjhyoXfnV6b9sG67dn9EXVoVrdOiLCc5efgMEu5y488YvjnXfKAbJxHkOmwoLBI6qHyQ2nNn4vDyZ5erVtbrv9Fk3EVC3jBnZyfY1cDoPBpAYWDB4Jh8gyby7DGEmHeYP7yYzfcFeDYVJIOqwbbEJh/IYFA8OkALvcyX7DYoHxGxYMDJNC0tGbZ4WB8RsWDAyTAg7XRjO0pcODkk1JjN+wYGCYFHCxEkJ7ZG9zNj+/YbnA+A17JTFMChg/qFPavNdYLjB+w4KBYbIcNiUlxotXlmLb/uqgi5GRsGBgmCyHxUJijB3A2RSt4DEGhslyWGFg/IY1BobJcjhuV5x/3zQSLZrkB12MrIcFA8NkOawxxBnSvTjoIjQK2JTEMAzD6GDBwDBZDmsMjN+wYGCYLIfHGBi/YcHAMFkOawyM37BgYJgsh+UC4zcsGBgmywmxysD4DAsGhslSmhaEgy4C00jheQwMk6VM+9UozFtTiVA6YnszOQULBobJUvp2OAp9OxwVdDGYRgibkhiGYRgdSQkGImpNRJ8Q0Q/K/9L56EQ0k4j2EtEHhvU9iWiRcvw7RFSQTHkYhmGY5ElWY5gMYLYQoi+A2cpvGY8CuEKy/hEATyjH7wFwbZLlYRiGYZIkWcEwAcCryvKrAM6X7SSEmA3ggHYdRbOLjAHwntPxDMMwTPpIVjB0EEJsAwDl//Yejm0DYK8Qol75vRlAF6udieh6IiojorLKysqEC8wwDMPY4+iVRESfAugo2XRnkteW+dgJq52FEC8AeAEASktLLfdjGIZhksNRMAghxlptI6IdRNRJCLGNiDoB2Onh2rsAtCKiPEVr6Apgq4fjGYZhmBSQrClpOoBJyvIkANPcHiiEEADmArgokeMZhmGY1EDR9jnBg4naAHgXQHcAGwFcLISoIqJSADcKIa5T9vsCwNEAmgPYDeBaIcQsIuoF4G0ArQF8A+ByIUSNi+tWAvgxwWK3RVRbyUayuexAdpc/m8sOcPmDJJPK3kMI0c5pp6QEQzZCRGVCiNKgy5EI2Vx2ILvLn81lB7j8QZKNZeeZzwzDMIwOFgwMwzCMjlwUDC8EXYAkyOayA9ld/mwuO8DlD5KsK3vOjTEwDMMw9uSixsAwDMPYkDOCgYjGEdEaIionIqtgf4FCRN2IaC4RfU9EK4not8p6aRRbivKkck/LiOiEYO8AIKIwEX2jRtK1iqBLRIXK73Jle0mQ5VbK1IqI3iOi1co7GJEtz56IblHqzAoieouIijL52RPRy0S0k4hWaNZ5ftZENEnZ/wcimiS7VhrL/6hSd5YR0b+JqJVm2+1K+dcQ0Vma9ZnZLgkhGv0fgDCAdQB6ASgA8B2AAUGXS1LOTgBOUJaPArAWwAAAfwEwWVk/GcAjyvLZAD5CNLzIcACLMuAefg/gTQAfKL/fBXCJsvwcgF8qyzcBeE5ZvgTAOxlQ9lcBXKcsFwBolQ3PHtEYYxsANNE886sy+dkDGA3gBAArNOs8PWtE5z+tV/4vVpaLAyz/mQDylOVHNOUfoLQ5hQB6Km1ROJPbpcALkKaXOALALM3v2wHcHnS5XJR7GoCfAFgDoJOyrhOANcry8wAmavaP7RdQebsiGn59DIAPlA95l+Zjib0HALMAjFCW85T9KMCyt1AaVzKsz/hnrwiGTUoDmac8+7My/dkDKDE0rJ6eNYCJAJ7XrNftl+7yG7ZdAOANZVnX3qjPP5PbpVwxJakfjoptJNdMQFHvhwBYBOsotpl2X38FcCuAiPLbLoJurOzK9n3K/kHRC0AlgFcUU9iLRNQMWfDshRBbADyGaPSBbYg+y6XInmev4vVZZ8w7kHANoloOkIXlzxXB4CmSa9AQUXMA7wP4nRBiv92uknWB3BcRnQtgpxBiqXa1ZFfhYlsQ5CFqGnhWCDEEwCFYJ54CMqj8ii1+AqJmis4AmgEYL9k1U5+9E1blzcj7IKI7AdQDeENdJdktY8sP5Eg7PT4AAAH2SURBVI5g2Aygm+Z3xkZyJaJ8RIXCG0KIfymrd1A0ei1IH8U2k+5rFIDziKgC0fhXYxDVIFoRkRrFV1u+WNmV7S0BVKWzwAY2A9gshFik/H4PUUGRDc9+LIANQohKIUQdgH8BGInsefYqXp91Jr0DANHBcADnArhMKPYhZFH5VXJFMCwB0Ffx0ihAdMBtesBlMkFEBOAlAN8LIR7XbLKKYjsdwJWK18ZwAPtUVTzdCCFuF0J0FUKUIPp85wghLoN1BF3tPV2k7B9Yb0kIsR3AJiLqr6w6A8AqZMGzR9SENJyImip1SC17Vjx7DV6f9SwAZxJRsaI1namsCwQiGgfgNgDnCSEOazZNB3CJ4g3WE0BfAIuRye1S0IMc6fpD1LNhLaJeAHcGXR6LMp6MqCq5DMC3yt/ZiNp/ZwP4Qfm/tbI/AXhauaflAEqDvgelXKch7pXUC9GPoBzAPwEUKuuLlN/lyvZeGVDu4wGUKc//P4h6umTFswfwZwCrAawA8BqiHjAZ++wBvIXoeEgdoj3naxN51oja8suVv6sDLn85omMG6rf7nGb/O5XyrwEwXrM+I9slnvnMMAzD6MgVUxLDMAzjEhYMDMMwjA4WDAzDMIwOFgwMwzCMDhYMDMMwjA4WDAzDMIwOFgwMwzCMDhYMDMMwjI7/D2CzOYrCrGafAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(transform[:,3])" ] }, { "cell_type": "code", "execution_count": 94, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([0.58566725, 0.22671932, 0.09035797, 0.01797021])" ] }, "execution_count": 94, "metadata": {}, "output_type": "execute_result" } ], "source": [ "pca.explained_variance_ratio_" ] }, { "cell_type": "code", "execution_count": 95, "metadata": {}, "outputs": [], "source": [ "A = transform#[:,0:3]\n", "f = fluxes[1]\n", "ATA = np.dot(A.T, A)\n", "invATA = np.linalg.inv(ATA)\n", "A_coeff = np.dot(invATA, A.T)\n", "weights = np.dot(A_coeff, f) ## same as xhat" ] }, { "cell_type": "code", "execution_count": 96, "metadata": {}, "outputs": [], "source": [ "science = f - np.dot(A, weights)" ] }, { "cell_type": "code", "execution_count": 97, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 97, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnXd4VMX6xz+TXkgBEiAQMPRepfdeReyAHfWq14rl+kOxIHq9XHvviuXasKMgCAgonSC9BwgQeigJJKSf3x97zsnZ3bMt2WQL83mePDk7Z3b3PbPnfGfmnXdmhKIoSCQSiSS4CPG1ARKJRCLxPlLcJRKJJAiR4i6RSCRBiBR3iUQiCUKkuEskEkkQIsVdIpFIghAp7hKJRBKESHGXSCSSIESKu0QikQQhYb764qSkJCUtLc1XXy+RSCQBybp167IVRUl2lc9n4p6WlkZ6erqvvl4ikUgCEiHEfnfySbeMRCKRBCFS3CUSiSQIkeLuBpuyzrBm3ylfmyGRSCRu4zOfeyBx6ZvLAcicMcbHlkgkEol7yJa7RCKRBCFS3CUSJyiKgtzQRhKIBLy4K4pCflGJr82QBCl3ffE3jR+d62szJCbMSj9I2yfnUVomK18zAl7c316yhzZPzudUXpGvTZEEIb9tOeprEyQOmDZ7K3lFpZwvLvW1KX5JwIv7LxsPA3D4zHkfWyKRSCT+Q0CL+7HcAsJDLZewcs9JH1sjCWYmzVxDUUmZr82QGNCGQqRbxpyAFffsc4X0eG4Rmw/lAPDvudt9bJEkmFm88wQtHv/N12ZIDChYRL1MirspASvuZ/Ltfez7svOk711SLfy84RB5hYE3kD9z+T6+W5flazO8gt5yl9FMpgTsJKZSkx7yoBeXAHKykaTq+HRFJsfPFvDW4j1c0aUBL1/Tyaf2rNt/ivqJ0aQkRLuV/+lftgFw1cWpVWlWtSLdMuYErLjL8EdJVbLlUA5t68fbpT81e6t+fOi07wfxr3xnJVHhIex4ZpSvTal2NEmX4m5OwLplzhc5Dn/KyS+uRkss3P3F3wx9eSlgGQ/wBEVRWLPvVKUny+w9cY5ps7dKH2QFWZ6RzeasHJbuOsElbyzjxo/XOM3vL6VcUOybgd6jOQW+DWSQA6pOcSnuQoiPhRDHhRBbHJwXQojXhRAZQohNQogu3jfTnnwn4n7wdH51mGDFnM1HyDh+jp/WH6LrswvZePCM0/yn8oooUX1LczYf4Zr3VvJteuV8oXd/uZ5PVmSy+/i5Sn2Ot2nz5Dzu+fJvq7SC4lImf72eIzkVb/0WlZQx5KUlLN11orImAnDdh6sZ++Yydh7NBeCv3dlO8/ti5urpvCJ+3XSYKd9vsmrgrNt/irQpc8gy3PtnC6wbOcZzaVPmVLoRNO6tZUz8YFWlPsMZ/Z7/g8veWs602Vu56p0VDvOVSZ+7Ke603D8BRjo5Pwporv7dDrxTebOck1dYwm2fOd7oo7DEsfAv3nmck260rBdtP8aqvdatklcW7LJ6fejMedKmzOGad1fqacsyLIKw69hZh599vqiULs8s4Nk5lgif/SctD92+k3lW+crKFKuW0dbDOeTkF7Pz6FmrB9WYHzy72XcczaWgEpNA8gpLOJNfRGmZwuIdx1EUhcNnztPk0TlsUSOZ8otK+XXTEav3zdtylJ82HOa5uTsq/N1Zp/PZcyKPaQZXiTc4muNez+vwmQKX7sFdx87q91vW6Xy7e8odjuSc13uDt3y6lnu+XM/Xaw/y3d/ljYGPl2cCsG7/aXLOF9P0sbm0n/a7vpqpoij0/e9iq8/NOOH4HnWHY7kWm5w9b0UlZR73ZDUOnjrPhoNn+GRFJun7TzvMV+Jhy33xzuMUmw3aBRkuxV1RlD8BZ+vdjgM+UyysAhKFECneMtCMjVnOW8Uvzt9lmr7tcC6TZq7ltUW7OXAyn7Qpc/ht8xHTvLd+ms6E961bJa8t2q0L1l+7T9Bnxh8ArMksL56fNxwCIDREABYhX7rrBM/8uk3Pk3O+2CqvIz5evo+JH6zi3aV7WLrrBGNeX8aED1Yx4tU/rR7UlXtO8tzc7XpLTeumfrRsHw9/u9HqM7cfydXjtU+eK2Tkq3/x+E9byC8qIdempTfm9b9ImzKHDLUnUFamWLX2TucV0fap+XSavoCmj81l0idrmb/1GAu3H6NMgU9WZDq8Nq0CCrUUE3mFJU5FwoziUstnhKllbcuJs4X6JLeDp/JJz3Rv2eaZK/a5le/QmfNc9+Fqp3mGv/InI179C4CBLyyxu6fAMhHPWeRNr//8QddnF5I2ZQ7rD5Tf+0b3m/a7xEaEseVQjn4PbFKfFU2IjUSGherH+UUljH9vJZuzcvjPb9s5X1TKlkM5HHJjcuDZAse2PzhrA12fXVglrkLbUMiHZm1kzOt/OX3Pij3ZTJq5llcXmmtEZSkrUypcmXkbbwyoNgAOGl5nqWnmqukFnN1MACv3nuRIznmu/WA1n9/andSaMQCcUAt9x5GzbDpkuel/2XSYhrViiAgLoUXdOJfffckby4iLCnMYbaAJTmiIYMPBM1z21nL93ONjWnMst5A5aoVyOr+Y/SfzdAH6dEUm/zeylfo5ZbqozvitvHW7/UiufpyZnUdaUqxd11hriWsVyotXdwTg0R828dWag9zSpzFPjm2jh42u2XeKAS8s4cTZQqtIo62HLd/11ZoDPHFJG95cnMHLC3ax7vGh1K4Ryb++22R3/esPnOa9P/cC8N26LB4e3tIuz+CXlrD3hKWXEqIKc9un5tOqXhzzJve3y7/9SC7FpWV0SE20StdaX2Gh1m2UktIyQkMEN89cw9bDuQxomUy/5y2VoTuRVJ708o1i6wjtYTdrYW45lMO9X60H4K6BTXlE/f3dwfh5+09ZytO2R/uf33YwoEUyw1750+79xaVljHjlT27r15jT+UWs3neKsW8uA+C9pXv1fD2b1GLV3lOsmDKY+on2UTlnC0pIqhHJX7tP0Dolno0Hz/DS77t4ZXwnffmGE+cKqRsfBcCy3dmczi9i9/Fz3D+kOaEhgvUHTtO0Tg3io8IdXm9uQTHnCkpYuP0Y1/e4yC4U8vu/Hbs1i0vL+H3rMT1vZrZz1+2lby5jcKs6TB7awip94bZjJMaE07BWDHXjo9h48AwpiVHUiYvi/q/Xs2DbMfKLSlnz2BDqqNfrK7wh7mbNJtPHQwhxOxbXDY0aNarwF97x+TqXeSbNXMu+7Dzu/3oD6/afZvVjQyhURa9UUdCei7mbjzJ3s+UG3PnsSCLDQl36Us8WlDBT7QY74v6vN9il/brpiP4gawx4YYl+nF9UykOzNjq9SY0MfHGJqVjZrrVxvqiUotIyvlpjqYO13kf2OYu4h4UITpx13Nr4aNk+Hhvdmjmqa+VITgEHT59n4fZjdnkXbLNO6/mfRfrxe0v30LZ+gi7sACFCsPOoxT2w4+hZ9mXn0Tgp1uozRr1maY01q1ODh4e3YGS7FF6cv5P4aMvtGxFafgvmnC+m49O/M3V0a33socDJ+Iw3SM88ReuUeMJChVVr2IjZ7Na8whIueWOZ/vrtJXvokJrIyHb1AIsgHTlT4PB7S8vKP/OMA/95aZnClB82m577ddMRdh47a1pJG1m119Lj6T3jDwa3qsP+k3m8eHVHIsNCKCwpI/d8MdsO53LDR2uoFRuhNxpGvFpeofR4bhHpjw8lPiqc6z8q7+3M2XSYtNqxLNpxnKGt6/LmtZ15ZcEuvYFgpMO03/Xjdg0S3IqWWZGRzVtLMlieYXGHjW5fTz/31M9bmLP5COmPD2Pb4Vy+WXuA63teRI2oMDZl5bApK0cX96KSMsJChFXlmTljDOPeWk5iTDjrnxjGzxsO6+f+3J3t83BTb4h7FtDQ8DoVOGyWUVGU94H3Abp27VqloyA7VMFYp/rqHv9pC+M61QcsLgGzbuKcTUe4oksqg19aqqfNXO5eF90dbIXdDHeFXePtJRl2ae//uZcejWvrry954y8+vKmb/joy3NLSPaC29vZml4tt7/8s4nBOAXufG231mVe+s4KIMMv7rnp3hcMIDeNn2fKf3+z96yWlZbwwvzx9kIMKCyDj+Dnu/N/fbHl6BG8uLr/u8NAQ8gpLiI0M4zl1HGNW+kFdUB0Nvk+bvZX5W49yS5/GvPHHbod2u+IqdcylRd0a/P7AANM8xtmtZWUKISHCdCD4zv+tY8czI4kKD+X+r9frDQ8ztF4iOO/NOhK/j5Z5fm//seM4AJe/XT7AOe6t5bx8jaV36GwSYddnF7Lu8aFWaXtO5LFHrewXbj/Gkz9vYZYbgQX/990m/bpsr+/mmWvo3zyZnUfP8tuWI+QaykYrzzkGd+zWwzmMed1SyX660n7v6TP5RXSavoB+zZOs0u9WgwTO5BfbBTE8/O3GoBD32cA9QoivgR5AjqIoVeaSqejg34Jtx2iSbGkRlinlfm8jWoN9n0Gg3lmyp0LfV108P2+nXdpfu7Ot7N5zIo8p35e3zjJP5pE2ZY7p5x3OsbQUmzxmvcztBkP0jzdD737aYN8OUBQFIcz96GBZ58XIzqNnafvUfF68uiPfpFt6J8aH7WReea/Esj47FJWW6WMCxqUrWtStwa5j5e/914iWNE2O5c7/WUf7mKG97+CpfGZvPExclPnjVVRaRlRIqMNw3lZPzHP5XeBcSI1scBG55Q3+O8+9gfGLn13o9Lw7wg7Wv29pmcL8reWV4JKdJ1iy0/0Iqq/WHHB6/ovVlvO20VNzDEECw03cXj9vOERcVBiDW9V12xZv4lLchRBfAQOBJCFEFvAUEA6gKMq7wFxgNJAB5AOTqspYgA8M3bVmdWrofml30Fo3iqJYTUbROHAqn89XZlqlHXfirvBnXrEZMFpt2AP24CnfT75xxqm8ImrXiORcYYneUjSyNtM6cuKsOhhpO3is8cL88gpw17FzfLYyU39gbYmNLH8kVj46WJ/5mVQjQndj1YgM45yDAdCS0jLdv++IwpIyosJDKSmrXCVZkZZ3VWE2YFtdvL1kj5070BP+t8qxuL++aDfxDippV2iu2ZiIUK7sksozl7Xj102H+X3rMV6f2LlCn+kJLq1WFGWii/MKcLfXLHKB5hoAqBsf6ZG4F6vddEfd1NcWVbxrLvEerlp3nqL5jMHaD2xG/cRofZDUOKX/5t5pvPj7Loa0qsM/BzYl82S+aWXiTmt66o+b7UJDqwLNJ+6PXNkllXoJkby1uPI948oIuyteXlD5qJr8olI+X7Wfiy+qyeRvLIL/2OjW1Euo2gHXgJuhGmLors+4ogPX93R/YDZfdelUNCqrV5PadmkrHx3MN7f3ZMYV7a3SHxzWwi6vkZt6XcTkoc354a7eetrU0a0rZpiX+PzW7rwyvqNV2uShze3S3MWsvDS6p9Vi4YP2kTG+5sWrzK/1nsHNyZwxho9u7kbXtFpc2aWBab7uzy0yTTdSHcK+69lR7HzWekmCW/o01o/v6N/Erc9pkxJP0+RYBrVMrrRNPZvU0o+v79nIKpKqdmxEpT/f39GEHWD/ScdjU94i4MR9UKvym6xhrRievaw9r47vxIIH+jP3vn5O36v5yIzhhJ7QrXH5zfnrvX354MaupCRE06NJbSZ0b8S711/M+zdczNanR3DfkOb0bmovbmEhgswZY3h6XDsmD21Bl0Y1AbikQwodUhMqZNd7N1xsmh4VHkLmjDEseMA9Ea0ZE8HlncsHgb6+vSeTh7bg8s6pfHhjV2ZO6maV//4hzWmqjmP8o19jq3OZM8bw1e09rSovgKcvbcu0sW348OauNKsTx8IHywcgX7q6YpWIN4mOCGVcp/oufwtnYwLe5AFDKJ5tGTvir0cG6T3cvx4ZRJzqaurdtDYxEaHc1rcxU0aVh1w+OKwF9w1uZvc5z1zWjl/u7cuihwbygslvM21sG/b9Z7Sd8L8yviMDWthXBhO7lzfEakSGIYQg/fGhPDC0BZ/e0p0ujRLt3uNthrSqox/XiYv0+P1LHh7IdT0qFulnHIM5mus4CspbBJy4N6sTR+aMMVYRFZd1bkDzunG0MVnoyZsMaplMUo0IZk7qRrsGCQxrYz1QMrJdPYa3raf7bd+5/mLevd4ivO0bJNClUSI/3d3H7nN3PDOS1yZ0JibC3Et27+BmPHd5e2zn6jRJjmXTtOG0qmcdn19D/f7wEMvP27xuHJP6pAHw2oROfHtnL/b9pzwapr/6INrGMPc0tLyHtqnLoJZ1GNm2PJSsR+NaJKsPyLhO5S3ZtVPLIyK6NKrJnudG868RLdn41HBu6p3GzX0a6/HMyYYHbHT7FMa0t5//tvRfA01KpZy4qDB9fsDHN3fV4/ptcTTZyZbXJnRm9j193cpb1QxpXS5GU8e0YcnDA63OP3FJGzo3SuSi2jF6WsNa1scDVUFrlRLHtukjefySNggh9PtkfLeGPGgzH6Fjw0Ru6HmRPhkvqUYk47s2tMpzKq8IIQRTRrXmyi6pek9tXMcGXNLB8jte0bn8vri0Y339uIYqdEk1Irl/aHPaNUjQv8sVYzqkkFrT9SqYxntJE+SPbu5GhDov4vNbe5i+TxP9EW3tB0LTkmI5ec7c9dasTg2n9hgjmqpjxYSAXRWyOln92BBqxUYQHhpC+uPD3H5fQnQ4I9vV43+39qBdg3gSY8y7nlHhltjo6AjzuvbBYS0QQjCmQwqKotBp+gIA/nhoIFA+Q++ugU25tW9jcs4XM/ilpRgbl09e0oYn1Yfals9u6e4yQkXj9YmdeW/pHi7r3ICGtWJoUS+OP3Ycp12DBN66tgs9m9Sidg3rFlFoiODuQfYtQ7CUkUZ0RChvXdeFMZuPcNcXf5NaM5qFDw7Qywdg45PDSd9/isU7j1MnLorLOjWgdo0IYiPDuKVvmh5nbuYPN5tE1CAxmkNnzvPpLd1pYhNf74pXx3ey6mp7k1/u6ctri3bTom4cr4zvyPwtFr9yWlIsV12cqq/JHhMRyo93WRoMz/66jVnpB+0+679Xtmdit4b6ZD5bIgyTwJokx/LCVR1olmw/oe++oc31aCSAwa0t4teyXhwvXdOR80WlHMstICRE6EJ3cVpNflhvmYkthCA0RFBapugVi5EQF/ffq+M7MbBlMokxETz242a+XH2AsR3r88vGw3RPq2U1Uxzgreu6MGfKHBokRvPvy9vz78strtPI8BCKSsuoZXAFbZ8+ktZPzqNmTDgrpgxm57GztK2fQH5RCW2enM+1PRpxZRdLr7Zb41rM22ofoto0OdbSK1LnFUSEhejhuOmPD6WrOpb0rxEtuayzuVvPmwSduI/pkKK7X4a3qcvvlRxs8cba8H1t4mMdoQlTj8a1eO6K9gx5aSnX92yki64mhM9f2YGlu8tDvRJjIvj7iWEkRocTEiL0yR1GsTYT7gUP9NcHD43nv7ithz7RyZaIsBDuHdJcf51UI5Jr1BbdmA4VW3WiT7Pa+iQTQA9ZVZTyiu/5KzsQHiZIiAlnSOu6DGlt36oym0Dk6B6YfU8fUhKiEcKylISZG8EVYzvWNxX3YW3qVnqQr31qAh/e1BWAyzunWrnLXry6I0UlZczeeJio8HJhfvySNjx+SRu7z4qJCKN3M/t7MC0phi2HcglXXTirHxtCjcgwq4ghIw0So1kzdQin84rJLyqhU0NrN0p0RChpagXZuVFN/vzXIBrWimbqj+VrDk4b24anf9lGrEkv1bbl/sxl7difnceHy/bx9nVdGG1oiU8b25Zb+jQmNjKUM/lFTB7anCvfWWn7kcyb3I86cdYDl29d24W3FmdQK9bSC996KIfoiFCmjW1D3+ZJhIWG0LZ+gl52Wi9Xe0Zu6ZNGq3pxvLJgF+n7T1MvPoqjuQXUiAxnQvdGjO1Yn/UHztChYYI+8SqpRqReAbnbg6wsQSfub13bhaxTy9iYlYOtnnVMTWBjlrlo+QMNa8Xw+a3d6XpRLaIjQtk2fYSpYF3TrSHXdLPuIhtbIZpv784BTZ1+X/O6cTQ3WXKhT7Mk+piIQVXx8c3drGK+Y8It9hujmmyv1xWbpg0nMiyE7HNFHMkp4OObu9Ht35aW099PDLMqL6NwekJoiGB0+3p2E41SEqK4uXea07V1NB4a1oKXKhCR0btpbWZvPEzrlIq7Ij+Z1J30zNN6K7quG9Pl68RF2YmlIxoZXEUpamTIDb3SuKFXmml+Tdwjw0KYPq4t47s1QlEU7hzYlCSb3mBEWIjeO9DcK1d2SWVwqzr65CKAVvXsy6d/i2TdFTmoZR0GtbS4rW7uYz6mYdswEkLQp1kS24/kkr7/NMfOWvznWkUbGxlm2qC7tV9j1mSeorM6zlbVBJ24A3x4UzcWbj+mT2vXCAsN4dd7+1pN+XZEg8RoehgGUKuLfs3LW5COfPCuiAwLDajdqCLDQq0qMW0GbWW2T9N8+g0So/nl3r4257x327993cV2E8JChGDKqFas2JNtNSHKlibJsdw9qBmXdqrPpyv2cyTnPCPb1TN1WdgyvltDhrapayd6npBUI1Jf6qAq2T59JCFujO49Nro1j/6wmS9u66H3HoQQbl/jS+os2bu/rLCpHqFVHInR4ZzOLyY81P4iP765K9vUNZpGtK3H2qlDrcaZqpKgFPfkuEgmdm/Emfwiq9ZTeKigXYPyKAitO2XGZZ3r868R7i/iJPEemjBrPk5vY7vQmLeYdUcvFm4/xt0DmxEVHsrvDwzgke82snrfKX1ZZ41f7+2r34sX1Y7lybH27hRneCJ6viY6wny9HVtap8SbBhz4K32bJ/HpLd3ZejiH5+ftNB0zGNyqrtUM1eoSdghScdewHcDUatZZd/QiPjqMVvXieX7eDt62WWLg7kFNHQ4ASqoezSUV5WARLn+lXYN4utv09p5X4+a11v2NvS7is5X7aehgcFMSWAxokayHVoeFVo8v3V2CWtxt0QYybB9AIzVjwmWL3Q+oqEvKl7hTGU0f147p49pVgzWS6kIbG3I3lLO6CLwnyEMy/j2K3IIShr/yJ/cZojw0UgxTgDs2TGTmzd3s8kgk7hDiZw+3pHoocbFpjK8IenEPCw2hVmwE6TZLjWpc1+Mi6sRHUScuktYp8VYx1RKJO3RLq2m3mJkt8yf397uWncQ7aAP//vb7Br24uyIkRDCibdVHDEh8T0pCFEdyvD/t+/Nbe9htUWhLy3qud/mSBCbapimh1bQkhbtc8OIuuXBY9NAAiku8P+87KjxU9vguYK7p2pBv1h7kcgeLyfkKKe6SC4aYiDAI/sUHJdXMRbVjPVqWpLoIuIXDJBKJROIaKe4SiUQShEi3jEQiCSpGtq1HTKQcA5HiLpFIgop3HWxec6Eh3TISiUQShEhxl0gkkiBEirtEIpEEIVLcJRKJJAiR4i6RSCRBiBR3iUQiCUKkuEskEkkQIsVdIpFIghAp7hKJRBKESHGXSCSSIESKu0QikQQhUtwlEokkCJHiLpFIJEGIFHeJRCIJQtwSdyHESCHETiFEhhBiisn5RkKIxUKI9UKITUKI0d43VSKRSCTu4lLchRChwFvAKKANMFEI0cYm2+PALEVROgMTgLe9bahEIpFI3Medlnt3IENRlL2KohQBXwPjbPIoQLx6nAAc9p6JEolEIvEUd8S9AXDQ8DpLTTMyDbheCJEFzAXuNfsgIcTtQoh0IUT6iRMnKmCuRCKRSNzBHXEXJmmKzeuJwCeKoqQCo4HPhRB2n60oyvuKonRVFKVrcnKy59ZKJBKJxC3cEfcsoKHhdSr2bpdbgVkAiqKsBKKAJG8YKJFIJBLPcUfc1wLNhRCNhRARWAZMZ9vkOQAMARBCtMYi7tLvIpFIJD7CpbgrilIC3APMB7ZjiYrZKoSYLoS4VM32EPAPIcRG4CvgZkVRbF03EolEIqkmwtzJpCjKXCwDpca0Jw3H24A+3jVNIpFIJBVFzlCVSCSSIESKu0QikQQhUtwlEokkCJHiLpFIJEGIFHeJRCIJQqS4SyQSSRAixV0ikUiCECnuEolEEoRIcZdIJJIgRIq7RCKRBCFS3CUSiSQIkeIukUgkQYgUd4lEIglCpLhLJBJJECLFXSKRSIIQKe4SiUQShEhxl0gkkiBEirtEIpEEIVLcJRKJJAiR4i6RSCRBiBR3iUQiCUKkuEskEkkQIsVdIpFIghAp7hKJRBKESHGXSCSSIESKu0QikQQhUtwlEokkCJHiLpFIJEGIFHeJRCIJQqS4SyQSSRAixV0ikUiCEKEoim++WIgTwP4Kvj0JyPaiOdVNINsfyLZDYNsfyLZDYNvvT7ZfpChKsqtMPhP3yiCESFcUpauv7agogWx/INsOgW1/INsOgW1/INou3TISiUQShEhxl0gkkiAkUMX9fV8bUEkC2f5Ath0C2/5Ath0C2/6Asz0gfe4SiUQicU6gttwlEolE4gQp7hKJRBKEBJy4CyFGCiF2CiEyhBBTfG2PLUKIhkKIxUKI7UKIrUKI+9X0WkKIBUKI3er/mmq6EEK8rl7PJiFEF99eAQghQoUQ64UQv6qvGwshVqu2fyOEiFDTI9XXGer5NF/ardqUKIT4TgixQ/0NegVY2T+g3jdbhBBfCSGi/LX8hRAfCyGOCyG2GNI8LmshxE1q/t1CiJt8bP8L6r2zSQjxoxAi0XDuUdX+nUKIEYZ0/9QkRVEC5g8IBfYATYAIYCPQxtd22diYAnRRj+OAXUAb4Hlgipo+Bfivejwa+A0QQE9gtR9cw4PAl8Cv6utZwAT1+F3gn+rxXcC76vEE4Bs/sP1T4Db1OAJIDJSyBxoA+4BoQ7nf7K/lD/QHugBbDGkelTVQC9ir/q+pHtf0of3DgTD1+L8G+9uoehMJNFZ1KNSfNcnnBnj4Y/QC5htePwo86mu7XNj8MzAM2AmkqGkpwE71+D1goiG/ns9H9qYCi4DBwK/qw5htuOH13wCYD/RSj8PUfMKHtser4ihs0gOl7BsAB1WhC1PLf4Q/lz+QZiOOHpU1MBF4z5Bula+67bc5dznwhXpspTVa2fuzJgWaW0a7+TWy1DS/RO0mdwZWA3UVRTkCoP6vo2bzt2t6FXgEKFNf1wbOKIpSor422qfbrp7PUfM/lntOAAAgAElEQVT7iibACWCm6lb6UAgRS4CUvaIoh4AXgQPAESzluY7AKX/wvKz96jew4RYsvQ0IQPsDTdyFSZpfxnIKIWoA3wOTFUXJdZbVJM0n1ySEuAQ4rijKOmOySVbFjXO+IAxLN/sdRVE6A3lYXAOO8Cv7Vf/0OCzd/vpALDDKJKu/lr8zHNnql9cghJgKlABfaEkm2fzWfgg8cc8CGhpepwKHfWSLQ4QQ4ViE/QtFUX5Qk48JIVLU8ynAcTXdn66pD3CpECIT+BqLa+ZVIFEIEabmMdqn266eTwBOVafBNmQBWYqirFZff4dF7AOh7AGGAvsURTmhKEox8APQm8Apf/C8rP3tN0Ad1L0EuE5RfS0EkP0agSbua4HmavRABJZBpNk+tskKIYQAPgK2K4rysuHUbECLBLgJiy9eS79RjSboCeRo3drqRlGURxVFSVUUJQ1L2f6hKMp1wGLgKjWbre3aNV2l5vdZq0VRlKPAQSFESzVpCLCNACh7lQNATyFEjHofafYHRPmreFrW84HhQoiaas9luJrmE4QQI4H/Ay5VFCXfcGo2MEGNUGoMNAfW4M+a5GunfwUGQEZjiUDZA0z1tT0m9vXF0i3bBGxQ/0Zj8YUuAnar/2up+QXwlno9m4Guvr4G1a6BlEfLNMFyI2cA3wKRanqU+jpDPd/ED+zuBKSr5f8TlgiMgCl74GlgB7AF+BxLdIZflj/wFZaxgWIsLdhbK1LWWHzbGerfJB/bn4HFh649u+8a8k9V7d8JjDKk+6UmyeUHJBKJJAgJNLeMRCKRSNxAirtEIpEEIVLcJRKJJAgJc52lakhKSlLS0tJ89fUSiUQSkKxbty5bcWMPVZ+Je1paGunp6b76eolEIglIhBD73ckn3TISiUQShEhxl0hsKCwpZe+Jc742QyKpFFLcJRIbps3eyuCXlnI6r8jXpkgkFUaKu0Riw5p9luVZTuYV+tgSiaTiSHEPInILipm35aivzUBRFGalH6SopMx1Zj8kNtISZ5BXWOpjSySSiiPF3Y85knOe6z9cTc75YrfyP/jNRu783zr2n8yrYsucM2fzER75bhNv/rEbgCU7j7Pr2Fmn71m2O5vtR5ytjGxPQXEprZ74jbmbvbPWV1mZwpZDOcREhAJwtqDExTsk/o6iKBQUX5iVtBR3P+aNPzJYlpHNLxvdW0E067RlEbtzhdUnSusPnGbR9mN0fXYhA19YTHFpGWfyLZXRiXMWn/XNM9cy/JU/HX6Goihc/9FqRr32FwdO5jvMZ0vW6fMUFJfx4vydVukFxaUoisLx3AJ+3nDI7c97ds52LnljmW6/qwopUFEUpVrvEV/y0u+7aPXEPPIukOs1IsXdTziWW+Dw3OM/bSE98xQLth3j0JnzdH12AYNeXEJ6pvXS3aEhln0DyrzoDflp/SG+X5fl8Pzlb6/g1k/TyT5XSObJfL5cfYDHf7LsN7z+wGleW7jb5XfMNlRe/V9Y7LZtmtsnIqz8Ns7JL6bVE/N4848MJn2ylvu/3kBOvns9n4+X7wMgPNTyeUdyzrttSyDx6YpM2j01n8NnKnd9hSWlfj/o/E26ZZMkTyuzD//ay4o92VVhUrURkOL+8u87+eDPvb42w4rCklLOFznv/mUcP8vlby9n+5FcFu88rqf/uesEPZ5bxMJtx6zyGxfsvOrdlfzjs3TGvrGM7HNF7MvO46p3VwKwPCObGb/t0MW9VH3j4TPnOZNf8Ydv1d6TTP5mAw99u9Eq/Zp3V3Ldh6tM3zMrvXzHsR1Hz/LKwl0uv2dtpnv7S5w8V0hZWXmhFJVaxF0TY4BsdRD0h/WHOJpjqTALSyvWLS8sKePwmfOs23+6Qu+vCgpLSjl4qrx3c7ag2G5so7RM4eNl++zcEcdyCygoLuWnDZbK9EiOfYNCURSWZ2Tjzmqxt3+2js7PLHCZr6S0zGeuEe06zLZLcsazc7Zz7QerXWc0cOBkPifO+s8gfMCJu6IovP5HBv+eu93Xplgx/JU/af3kPB79YTPLdmfz5eoDdnn++b+/WX/gDKNe+4tJM9dSqgrV/32/CYBNh3Jcfs8pk5bSdR+u5t2lewgRllu4RBW93jP+YMALSzh5rpDjJj2D0jLFaetkwvvmAr4m8xTLM06yL9vet3/Ig9bgibOFHDpznojQULtzB0/l672ZHUdzGfP6X1z87EJGv/4Xxer1aZWpEOgPlVYGiqIQolZ2hcWedWU2q79DYXEZ/Z5fzJXvrPDo/VXJw99uot/ziykssVx7+2m/c8NH5SJUUlrGt+kHmf7rNq58Z4Uubi/O30mP5xZx31frdaGNDrcv9582HOK6D1fzrZPemsbSXScAXFYE//zib1o9MQ+Apo/N5eXfdzrN7020tkCZYqncjuScJ23KHJZneL9V3v+FxXR/bqHXP7ei+Gz5gYpSXOpf688/P28HDWpGs1/1FX+15gBfrbEI+7U9Gun5Dp7KZ/dx64kx5wpKSIgJ11tQkap7QVEULBvxeHatWsu90NCSyzlfzMXPWm64xJhw2qTE8+U/egLw5ZoDPPHTFt68tjMDWiSz9XAuPZuY76+sKAp/7s62aplrlZKRM05cIIt3HKe0TGHuliM0SYrlo2X7OO0gf7/nLe6ZzBljGPnqX3r6jqNnaT71NzY+NZz8IktXe1NWDj2eW8i6x4fpLbQyBcLU8jhbUMLXaw7Qs0lt0pJiHdpnyzfpB11nqiYURaGwpIw/tlt6dwXFZUSGWcR59b5TrMjI5toPV5McF6lXdFsP53LZW8uZ2L0Rby7OAOD3bcdoWCsagBCTpt3BU5bK2dHYx9JdJ7j7i7+t3BxFpeW2mLHA0CMtLbM0zh4c3tJhfm+hKIreGJqVfpCXF+zitr6NAfhi9X76NEvS8xaVlJF9rpD6idEefcdP6w/Rsl4crVPi1e/0kvFeIODE/fu/XbcoqoLlGdl8sXo/b0zsoovo8dwC3l6yx+n7CopLiQoPNRXCjtN/Z1DL8vV/Dp85zx2fpzN/6zH+fmIYv2x0HQUybfZWu7TCklK9V2DkTH4xK/acpKxM4WReke5Lv+fL9bStH8/Ww7n89cggGtaKsXtv40fn2qWFeNjXnfTJWrfyaS1CwPQ6QHUN9SyvPMsUWH/wNEdzyrvFWiv+P79t56/dlpZa5owx+vmc/GK2HcmlV1PzCs2f+HLNAab+uEV/XVhcCtHh+usf1lsGjm3dAhuzctiYtdkqTRPwEkND6cDJfI6fdTzuo/H6ot12/uuCIou4HziZT3JcJNER5kKv9SiriyWG++jlBRb34I6jlkHy80WlLN11gh6NaxEVHsqjP2zm+7+z2PHMSP35Bpiz6QhjOqTor9OmzOHWvo154pI2AEz+ZgNgfV8NfXkpX9/ek6QakVV3cW4QcOKuRYRUN9d9aOn6JtXYyh0DmvLLxsOMalfP6Xs2ZZ3h0jeXc+eApqzYc9I0z+Kd5TfgFwZXThc3fJkAn6zI1I81ISwoLuNsgeMW9DtL9zBvy1Hd/QCWVh5Y/N81YyOoEen61li1t2r2Yr7p4zX6cdPH7CsVgJ3HzvL0L9us0m75pHwhuuLSMt0Xv/OofdRL9rlCuqo9mis6N+Dl8Z0qbXdlOJVXRHrmKYa3tb+nVmRkWwk7QNaZ89SJj9Jfl1WgyVhiqDi1gWxtcHrbkVzSpszRz39ze096NKltel+cLy4lriyM/i8sZlDLZGZO6q6fyzS47p4yaYh4k3OFJRSXlJFXVELf/5oPzC9T3TGLd55g8c4TRIaF8N2dvVmwzTI/pKC4lDDDGM7dX/7NmA5jOJVXRM0YS2X60bJ9uribkXH8HL9tOcoNPS8CLPdiiBBWlUZ1EHDiHmrWl6xGvll7kA0Hz7ApK4e02s67+O8u3WP1v6rZcPAMANN/2cZTYx3ffMszsq2E3ciDszYCG/nXiKrvNlcWR616sAwWRoVb7hXjQHdZmcKP6w/RvG4NPe2H9Yf4WY3YqRsfybHc6h0UUxSFWz5Zy4aDZ9jw5DASYyKszn9hMn5zxdsrmH1PH/31D3+7H/KpUVJaxlM/b+HTleWLDGqDs3/sOG6V98mft/Li1R2telUa54tLOae6yLTzU3/cTO+mSTxsGIw3Xse5whK9oigtU2j62Fz+NaIldw9qxv6TeQx4YQk/3d2HTg0T3b6e/s8v5lReEeM61Xf7PYUlZYx9cxkJai+o0/QFrHp0iFWe03lFLhtbR20Gp9MzTzFn02ESoyOYt/UoA1ok8+kt3R28u2oIuAHVsGqu/WwpLCmP43Y1PX3uZt/MFj2aW8A/v/jb4fkok4E0W16YX3WDXmm17d0+VUGBOpCab4jU+HH9IR76diMzl2da5dUqivFdG1aLbUau/WC1XjGbzepVHIy9ZJ2uXChjUUmZlbA7Y+exs4x9c5npuexzhXq4aZliEcMvVh/g7i//5ryDKJkHVXfGqr0n9bkIb/5hGRfQXGja+M6U7ze5jI6776v1un/d0TiOM4wTBXv+Z5HVuexz9s95QXEpM9XQWbCEThr5ecNhVu09xbytFg3QKr2yMqXa3FMuxV0I8bEQ4rgQYouD80II8boQIkMIsUkI0cX7ZpZj7No4a7lVJQfUUDTbrnKgYNsqq076t0hmyb8G0axODdeZVerEVc53abxPNPfTgVPm7r0yBZ4a24ahretU6jvd4VxhCWlT5rByb7nLzuyW3n3MfIVKrWfiCREGl8O1H3oW6ueIq99daSWO7oRH/n3gDLkFxUx4f5XaW7REPQH6DGGtx/X12oMuo+OMcyX+NOldVIZcExfnibOFVm7BD5fts8tjxsPfbqTZ1N+8Zpsz3Lk7PgFGOjk/Cmiu/t0OvFN5sxxjbLlvPew6dFBS/YztWJ/Hx7Q2PXeROlh7c+80AHo0ruXy81ZMGew127SJSqcdxP8Xl5YxqU9jXry6o9e+0xGnztnbsC87j0d/2MyZ/CJKSi1x9rZRVhpbDnm2XAOUzw3wNrluLpGhkX2ukA7TfrdKKy4tY+byfXqQwI/rD7m99EZVYmaDFs3lKdrAd3E1tN5diruiKH8CzkbOxgGfKRZWAYlCiBQn+SvFqHblH+2o9eUNthzKIcPBQ+UNxnWqz8xJ3ars873Fjb0uYkI3c1dFZFgIDUxCxyZ2b+hSRLRWWqQbLiLjAJct1/ZoxB8PDXD5GbbsPWG+/o42oOgstM9bmA0fTfxgFV+tOUCn6QuY9MlaRr7qeNkGLQLEHzBr3XpKcanC079sI9ewps+stdbhqGv2nSJtypxqmz0cIuC4l8ZgMo6XD+yfNKnYvY03fO4NAOMvkKWm2SGEuF0IkS6ESD9xomJdp0a1Y9j45HDAfhDDm1zyxjKGvryU4tIynp+3w2neN6/t7PHn39G/KYNa1uHHu3pX1MRqYfq4dsy4sgPf/7OX3bn6idFM7G4v/MWlCsUl5f6FzBljmKYO8GqiroUpRoU5vwW/vdPyvVNGtbKLUJh7Xz+eu7w9TZLdd/G4oqY6mFkRl4enhLkIDvhrd7aV0FWWNmosti0PDmtR6c+uqha27e/wxWrLGMEqgyurojrgTqhivfgorzXyhr5cXlGbTUb0Nt64g81GOE2d4YqivK8oSldFUbomJ7vc39UhcVGWUXZv3viO2HjwjMtY9ks62I/O1413fONc2rE+bepbHjTtvzOc+ad7NanN9Wq8d0qCJTSusTpRJzREeNyqXfZ/g/Tjdg3Kbbv4olq8c10X7ujfRE8b0qoOdw9qxj/6Nbazd2IPi+jfqk4aGdjS4sO+sksqUH7TGFvuz13e3upz2jWIp1uaxW1z54Cm+mdp1K5hHVUC5feGLWY9DDMmqJWVEILR7cvDEt2Zju8p1R341djBBK6asfbl6Cn/9/1m15kqwDmbZZe1MTdjjL4jF5srbO9bMxJiItz2p3vC+eKq1y5v3F5ZgLH5lgq4t4xhBQkJEUSFh3C+qGoK6I8d5TPqHI32a2ghg1ufHqGn9WuexOrHhpIYUz7JJDy0vA4Uhuow3PCEPzOurdVnT+zekM3ThjsNa/zq9p48NbYto9vX4/0burL16RF8poZchYcKmiTX4LHRrZxeg5GUhHIR/P6f1r2KUe1T9PC1id0bMWVUK4QQdGpYU8+TOWMMDRKjqRMXReaMMXprOy0plswZY+iohrZpZSCA1yZ04pIOKSTZiLWx9a/x78vbMeOK9jwysqXVQOsr4y0+8pt6pVlNKAGYN7kffz4yyGXs/pj2KcRElOcxtqz9aeZhRYmPNr/+eAcVojsYZ2Fr1PJCZaHxX5teszbmZhwkv/tLx5FhzkgwTAIzQwg8Woa6X/Mkh+da1o2zem1baVUF3hD32cCNatRMTyBHURTvLLDthJiIMPJdLNRVUYyTYW74aI3DfDVjwrl7UDPAemXC9264GIBQg4rXjo2kQ2qCJd0wKBxiOE41mRkaFxVuNdnEjPDQEN6+7mLapyYQGxmmRxtoFcft/ZuSOWOM3c08om1ddjwzkpWPlg9YGrthZn7nSLWbnFwjQveFj+mQwjPj2vLAUPe798LwTeM6NeDNa7swrE1dXhnfUY/fLjFZ3vK6HhcxoXsj7hrYTF2iQf2Mjg146eqO3D+0uVX+yzs3oFW9eEJDhMuVARvZhGiGGSrkikwSckkVVRit6sXpPSQjcVHmYlYnLsrqtTbYbeTzW+1jtB8e3sJ0QbCaMc5FE+AyD2LRjWjzXKb8sJn7vloPOB4/cYVZeRgjrT39yZPjInl9YmceGtbCrqEysJW1pyK/GpYgdicU8itgJdBSCJElhLhVCHGnEOJONctcYC+QAXwA3FVl1hqIiQh1uQpjVdGsTg3WTB3C3Pv76WnGKB6t9Xe1IWY6KtwyE+7GXhcxdbR9JMmlHeubLuQEUGyIfXYnhFDbSWiiTatq/uT+Vq//NaIVUeGhVq11ISytsQgHg5g39Ezjlj6NuX1AU+v0Xml2wuoUE2eeEILLO6fqD52rSs1ISIjgyotT9VmpP6hjGVd3tRc5QF9fRWN0+3rcO7iZVdqRM+W+3NIqEPeq6gzMm9yfl67pyE9397GaRR1n6LkYe3MNEqNZ81j5xJ1pl7bl13v78vHNXfW0psk16J5mHdlUXGq+EYarFjHAvUM8uFcMGJ+z2RsP2w24ekK9BHvXqQJW7jhPiAoP5dKO9bl3SHP+eHig1blwGx9cXjVolzvRMhMVRUlRFCVcUZRURVE+UhTlXUVR3lXPK4qi3K0oSlNFUdoripLu6jO9QUxEaJW03M1m4Gnco7bSBZbWjrUoWm66iy8qd1E8MqIl26aP4PExrflkUnciwkKYPq4dtW0GcnY+O5JXx3fSW9y2dEurRXxUGN//szcLH3TtQ48KD2X79JFMGWntjqmXYN1Ci420/z4hBM9d3p5d/x5l+tnREaE8ObaNW8sTOENfvdHkXJiJX9VTujSqSeaMMfRuat5V/ub2XrwyviOpNS2/4ZSRra1cMmCJmNLwhrZnZuexxLDUs6ef2TTZ4jO/onMDq7EPR3RqmMg7119cPgZj6In8o1/5+8PDBHXio/j9gf76AHa7BgkMblWXd6/vQve0WtSJi2SqTXhrSVmZPlEMLG7AjqkJvDahs13j4PLO1jEW7kykM8N2Cv8jNms2XeRggtxv9/fjGpuKvqlhIL6NYeGvGVd24L9XWo//aCx8sL9pOkCUoacbHxXO/27tofcEbO2ujolMATdDVSMqPJSCEu+Lu3FdE1u0yiREmDQ7gSUPD9T93WBpTcZEhHFbvyZOVyOMDAslJMTx2hM1YyPYNG2EXnEMaWUZnDRWJLZER4RauXw0vruzPOqlsgJdGZzNM9Za31UVkw2W3s3lnVP1CtXsJz1r6Dp7wy0z8MUl3Dxzrb6Wv6vxHKPvevv0kdyp9ZaEea/ml3v6MusO+6im3urCaA1rlgufEEIffNaEuEXdOH0AW2NkuxRm3dmLsNAQ3a2oUVKmWIn2GxM78/M9fWlYK4bNTw/nmq6ptFUDBlrVs/Y5u4qScoSrGerGn2n9E8P049Yp8fz3yg48NbYNt6sVY0xEmD5Zbdql5eNd8VHhjO9m3evVxtZqxzoOlIi0iezp2zyJtVOHsvDBAfqYW2rNaPY+N5oJ3e3HKrxNwIp7aIio9hmq2lZdDrSdtKRY3SVSlXxwY1c2TRvO17f39Pi9XQ0PryMfbHVgHFC1JTkukv4tknl9guchpu6iicRHN3XjwWEt9Ba8Iyp6rx3NKbDb4m3STMvqmINeXOL0vV0aWQaf7+jfhOiIUL2XEyIEYzvWJzxUWEU3tU9NoLvJpLCnxrZl/uT+DGtT1/qE+oGuQjI1hBC8dW35BPSSUoWxHeuTOWMMO54ZyUjDHJTIsFCev6qjXlmE27Tk3ZnfYIajFr8W2WVcqsE2CkgIwaQ+jXlsdGsyZ4whIiyEF67qyJRRreiWZmkoXWdwZa42uKruGtiUXc+OsvpM2zVszOqd2jUiaVanBkPVsn/7ui6mja6qIOAWDtMIFd4V9+O5BRx3sYvKNd0a8k36QauBPG/Sul48dwxoQliI4K3FjsMvQ0IE8ZUQ5l/v7evzmX/OijA0RFj1gKoCbe38hrViuM+B//enu/vwj8/SOXG20HRZAHfo+Z9FtKhbg98fsHan7T3hOnb62cvaUy9hN3ep7sDaqrA0rBlDp4aJ7P73aMDyezobn4gIC6FlvTi7WZG9mtbm923HrIIBXDGmQwqvLqzB7uPnGN/NOKbkXKxtf29j9Ji7PPnzFtYfOGN6LibcImWxNq61SX3S6OVgjwKwVABaj2jnsyOtfON146OYPq4t7/+5FyEEEWHWNnsy5teqXrxdFFdVE7DiHuLFlvs1761kzT7Xy9e2VLuW99kMvHmLkBDBo6Nam+7i5E3aNUhwnamKceZzr2o2TRvudNarRqeGidw1sClP/7KtUnHuu0zWhhn80lKX76uXEMWzl5X7fge3qsP7N1zM4FbW6964+3vauhNfm9CZrNP5Dtdfd8QCN8Z9nDF9XFurSKxhbeqyau9JzrqYt/KZk0XONBdevM1g7lNj25plN8UsOuzGXmnc2CvNKm3y0Oa8vmi33kDq3bS2wyW9fUnAintYiKj0+gyr954kt6DELWEHi4+6umtfiffxpM2oCWJ1uwBfm2C/vrwQwnS9d3ex9QZER4TS3Cb+2tsYK8VXx3eied0atK1vXRkJKr/WSvM6Ndhw8AzXdm/k9vNcUSYPbcHkoS2Yt+Uoq/edok1KvF+Ke0D73D0JlTNj/Pur+Mdn1RLcUyGCYeKMP+KJW03zj1altk8ZZT/JbFwn0xU8KkVVuRPd5bLODeyEHSy9t8pun5mWFMve50ZzWWfvl5sjRrarZzV/RHi8DXfVEtDiXuajJX8l3sMXj4NnLXfL/6pYfkDD19ux+QNaz2i47aCvm4SFiGobqAwUAlfchWDPibwqjxd1Z0naqsLHDa2gxZNy1dwyths4uCKvsITDZ8pXLlxo2CTazh6PPjn4eNoQhvjcFebx5a6o7i3sjPhrEzNgxX3NvlOcKyzhu3VVu2H2u9dfXKWff6HjiwfDk+6ztoREmeJZ6338+yvpPeMP/fVtTtx/F3olXt+wqJtttEsgoN0W/vY7Bqy4axNMNmZVbMOOjQfNQ6pssR19lwQ+njyExry5591fD8SdjTS0xaT8TRR8SUWXWvb1eII/ErDirlHR3pi7G334srsn8T3G8MH527y3J+6LV3fU1wky9iR+uKs3c+/r5+htQYu2WFlFRVo+pfYEXh/IBkdLAbjC1XTyDqkJDGpZ9ftoOuNCiJbxyYCqB19qXBnyke82cY2LDbT3nDjH71sd+9c1jDNijfZ0aeR4SYlgZtqlba2WAPAUX7bBtFmx/lbBBIG4V+x9rsT9wWEt9A0mbunTmFN53tlqS2KNv/vcIz1cA2Xi+6tcznQG++n4c+/rVy0bOFQ31eUukW4ZewJe3H/bcpSnx7Xz+H3OgmzGdEjRhR3gSSebZUgCD090wNM1UFzNstSoGROuN/WEEG7tyBWMpNisVFpRpLbbE/Di7k4ryQxnMfKRbkxNrw7kDVs1eFKsUR5ulK246Is0Torl35e1s9r39UL9mddOHerx0geO8Isy9LMH1j9UrJooLVNo9thcvlpzwOnmC1WxMYPEf/CkC+/JolrukH22kN7NrNeYD+a7zVn4aHJcpPeWnfYzYfUHgkLcB76w2G5ZVTNO5hVSUqbwwvydTpcu6Kzu8ympevx9hqqt9KZNmWMqWD+tP8TqvSddDoIb13CXcuQefz0yyGUeX5alv7YFg0LcM0/mszwj22W+47kWF07t2AjmbHK8h/dNJntI+gJ/vWm8iU8GVD1QArM2gNYw2HIoh7Qpc1i99ySTv9nA+PdXubyeyq6HdCGSHGe/PEOsjTvHHxrufmCCFQEr7i9d3dHq9cm8InYePcv+k443y9WW6KwRFcaqvY5XjvP1yLs/3KjBjCe/b/sGCfpuQhravqHL1AbFoh3lW+fVi7cfIJx5c7eKmClRMdt9aev0kVavKxoS7Q38tboOWHG/8mLr/RDLFIURr/7JgBeWOHyP1mo6eOq8wzz+wIXQYg8UosJD+eK2HlZpS3Za77NrlBXb3X8ABhnWXzfToKpclCwYsBVuMzeNP7SH/K1RFrDibos7z4e2yFj2ucCIWfe3m+VCxXa1wXu/Wg+U33Pv/blXP+dqWQvj7FNf9xADBWMxJUSH07CW/SbYsijtCSJxd63u0t8pqQhmboHcgmJ+r8ByBK1TLsx49sogRPm2i46ec5+upe6nPa+AFnd9N3gcb6awau9J0qbM4XhuASWV3BBAcmFi5s99aNZGh/t5Gqll4qaReI62xpNDHfWDlru/bdYR0JOYjHONHC0n8NnKTABW7zvlcskBsB+o9SV+2iC44DBruWPj2VYAAA0dSURBVK9ysq3amPYpPDCsOQnREV6bpHOho1Wwjp5h/5JV/yCgxd3YojK23BVF0f2ZEWoNUFTielOPJkmxdgO1vkD6D/0Ls5VBzzqZVxEdEUqzOq73JpU/s/voO2I5OO/L8Qt/bYMFtFvGKO7P/LpNPzb61rUZhg99u5F8wwQSicRdPBUO2ePyPqEhzlvu/rAyt781yoJG3I0Yd6o3rr63yUEkw6KHBnjXsEoixSGwGdOhnkf5L4Tfu7K6p1WwjsrK34TVHwhocXe0vleJA3EPtGgZecMGJoNbubfJs/x93cdVWflyMNNfK+eAFndH3eVSQ1SMcT3ujOPnqtwmyYXNwJbJvjYhKHEl3f5QUfqBCVYEzYCqkeKy8sFTYwWw+ZD5fqtCz+s107yCv7YIJOa8d8PF9Gpa29dm+BVt6ycAWC1xXBG0Z10+Eu4T0OLuSIyNPvctDgQdLNExV3RpQOOkWG7r25hrezTytokVwt8qGQk8NbYNUeGhPPrDZod5mibHEh/leEN1R8sHu1oDPpC5umsqnRol0qKu6+ghZ7h0y/g0WsY/f7+AFvdiB+GNxaVl/OOzdBZsc76X5R8PD9SPH79E7rYkccykPo0pKilzKu6RTjb22Pr0CDuB8oUcdUxNoK7J4mZVhRCi0sIOrhcGk9Ey9rjlcxdCjBRC7BRCZAghppicbySEWCyEWC+E2CSEGO19U+0pdCDus9KzXAq7ROIpEWEhXNapvv7atqfnbGOP2MgwYiJ835b6+Z6+vH9jV1+b4RZvTOzMT3f3cSuvv80O9QdcirsQIhR4CxgFtAEmCiFsm7mPA7MURekMTADe9rahZhSWmMetv75od3V8fZUhfe3+i7H7f5XNhDdPN9PWkL+3OT2b1KaTunGO3nI3lJWxpezLVrO//n7u3I3dgQxFUfYqilIEfA2Ms8mjANqKSAmA450wvEhBsetZpxKJN7mma0MA7hnUzG7HLk+35JOrQjrHE/H2h5L0t9/TnbuxAXDQ8DpLTTMyDbheCJEFzAXu9Yp1LiiowIzTF67qAFh2Y/JX/OwekRjo1bQ2mTPG8PCIlgghWDFlsH4uwk82Vg8WjI+B2TNh9MP7tOXuu692ijtOQLNis72eicAniqK8JIToBXwuhGinKIpV01oIcTtwO0CjRpWPTHHkc3dEfFQYV3dtyLA2da0mN0kkFaV+YrR+HCbvKa9iFO/yUMhy6bEWJtkissWduzELaGh4nYq92+VWYBaAoigrgSggySYPiqK8ryhKV0VRuiYnV36yhyOfuyMSYixhaokxEcR6a9d1iaSS+KvP1tdYuWU8yFvdNK4dC2C6iYgvcUfh1gLNhRCNgUNYBkyvtclzABgCfCKEaI1F3E9QxXjaco8Ol8uvSrzPy9d0ZN4WzzfukG1N5xgjYMzWlrG05hXDsW+4umsqjZNj6XpRTZ/ZYIbLlruiKCXAPcB8YDuWqJitQojpQohL1WwPAf8QQmwEvgJuVqphY8ixHeq7zmTAWRyyRFJRruiSGjDhhYGAvsSyqwFVD1r2VYkQgm5ptfxuQNUt34SiKHOxDJQa0540HG8D3AtI9SJXXpzKQ99udDt/RUPVLgRm3dGLnPPF1fZ9UWovKj5Kusck1oSHCkrLFCtBd7X8gJ/pql9wQT1ZkeFS3B3RvXGtav2+Ya3r8viY1kzs7h9LPvgS6XK3JjwkhALKrFwwZtod4qplf4FzQYl7o1qxvjbBLeolWKaHpyUFhr0VISREcFu/Jr42w7dIQTIlIiwECqGk1LgAoH0+K5+8LEw7glrcE6LDrVwNT40NjPVjBrWsw5e39aBHE7nCoOTC48t/9OT7v7OsNhfX3TIOd+uoDssCi4D3U9wzqBl14iJNz6UkWC+QFBVA0TK9myWZ7t0pCT6qIfYgoGhZL47HRrd2OUA5uFUd/Vg+KfYEvLg/PKIla6YOtZopqBFpEPPN04ZXp1kSiUukK8F9QkwaOi9d05EodRzNl6GQ/krAi7uG2W/bzRB3GudknW2JROLfmEl3VHgojdSJQ1Lb7QkecTf5+cd3a2iSUyKRBBouQyFlL8iO4BF3s4WFpM9aEgBIj7trXO/EVD12BBJBEy1j9tuGCsGsO3pRQ64jI/FDpCC5j8u1ZarFisAiaFTPbGmBECGqfXKORCLxPmZry1hnqD5bAoWgccskxITz2S3drdJCgubqJJILG0e9HE3sZbSMPUElf/1bJDN/cn/9tfzBJf7M/41sxRWdG3BpR88WwLsQcfUsyyfdnqASd7BMgNAmNUlxl/gzyXGRvDy+U0BNrvMVrtdzl8+6LUEn7lC+ZGiZnPknkQQFjhpqWrLUdnuCUty1G6G0TIq7RBIUBMAG2f5GUIp7fLScjSqRBBMyzt1zgiYU0siHN3Xll42HSa0Z7TqzRCLxe1yPn5Wfv75nI2IiglLaPCIoS6BBYjR3DmjqazMkfsbL13Tk4KnzvjZDUgEcSbs2rGbU/mcva1/l9gQCQSnuEokZV3RJ9bUJkgriqOVeUFIKQLSMOLIjKH3uEokkuHDklTlXUAJAnNyL1w4p7hKJJGA5V2gRdxlEYY8Ud4lE4vc4csto6TXkAKodskQkEonf48gt89PdfVix56Rc3tsEKe4SicTvcdRyb50ST+uU+Gq2JjCQbhmJROL3yHa550hxl0gkfo+cgeo5UtwlEonfI1d99Bwp7hKJJGD450A589xd5ICqRCIJCDJnjPG1CQGFbLlLJBJJECLFXSKRSIIQKe4SiUQShAjFR1vRCSFOAPsr+PYkINuL5lQ3gWx/INsOgW1/INsOgW2/P9l+kaIoya4y+UzcK4MQIl1RlK6+tqOiBLL9gWw7BLb9gWw7BLb9gWi7dMtIJBJJECLFXSKRSIKQQBX3931tQCUJZPsD2XYIbPsD2XYIbPsDzvaA9LlLJBKJxDmB2nKXSCQSiRMCTtyFECOFEDuFEBlCiCm+tscWIURDIcRiIcR2IcRWIcT9anotIcQCIcRu9X9NNV0IIV5Xr2eTEKKLb68AhBChQoj1Qohf1deNhRCrVdu/EUJEqOmR6usM9XyaL+1WbUoUQnwnhNih/ga9AqzsH1Dvmy1CiK+EEFH+Wv5CiI+FEMeFEFsMaR6XtRDiJjX/biHETT62/wX13tkkhPhRCJFoOPeoav9OIcQIQ7p/apKiKAHzB4QCe4AmQASwEWjja7tsbEwBuqjHccAuoA3wPDBFTZ8C/Fc9Hg38hmXJ6p7Aaj+4hgeBL4Ff1dezgAnq8bvAP9Xju4B31eMJwDd+YPunwG3qcQSQGChlDzQA9gHRhnK/2V/LH+gPdAG2GNI8KmugFrBX/V9TPa7pQ/uHA2Hq8X8N9rdR9SYSaKzqUKg/a5LPDfDwx+gFzDe8fhR41Nd2ubD5Z2AYsBNIUdNSgJ3q8XvAREN+PZ+P7E0FFgGDgV/VhzHbcMPrvwEwH+ilHoep+YQPbY9XxVHYpAdK2TcADqpCF6aW/wh/Ln8gzUYcPSprYCLwniHdKl91229z7nLgC/XYSmu0svdnTQo0t4x282tkqWl+idpN7gysBuoqinIEQP1fR83mb9f0KvAIUKa+rg2cURSlRH1ttE+3XT2fo+b3FU2AE8BM1a30oRAilgApe0VRDgEvAgeAI1jKcx2BU/7geVn71W9gwy1YehsQgPYHmribrdjvl+E+QogawPfAZEVRcp1lNUnzyTUJIS4BjiuKss6YbJJVceOcLwjD0s1+R1GUzkAeFteAI/zKftU/PQ5Lt78+EAuMMsnqr+XvDEe2+uU1CCGmAiXAF1qSSTa/tR8CT9yzgIaG16nAYR/Z4hAhRDgWYf9CUZQf1ORjQogU9XwKcFxN96dr6gNcKoTIBL7G4pp5FUgUQmhr/xvt021XzycAp6rTYBuygCxFUVarr7/DIvaBUPYAQ4F9iqKcUBSlGPgB6E3glD94Xtb+9hugDupeAlynqL4WAsh+jUAT97VAczV6IALLINJsH9tkhRBCAB8B2xVFedlwajagRQLchMUXr6XfqEYT9ARytG5tdaMoyqOKoqQqipKGpWz/UBTlOmAxcJWazdZ27ZquUvP7rNWiKMpR4KAQoqWaNATYRgCUvcoBoKcQIka9jzT7A6L8VTwt6/nAcCFETbXnMlxN8wlCiJHA/wGXKoqSbzg1G5igRig1BpoDa/BnTfK1078CAyCjsUSg7AGm+toeE/v6YumWbQI2qH+jsfhCFwG71f+11PwCeEu9ns1AV19fg2rXQMqjZZpguZEzgG+BSDU9Sn2doZ5v4gd2dwLS1fL/CUsERsCUPfA0sAPYAnyOJTrDL8sf+ArL2EAxlhbsrRUpayy+7Qz1b5KP7c/A4kPXnt13DfmnqvbvBEYZ0v1Sk+QMVYlEIglCAs0tI5FIJBI3kOIukUgkQYgUd4lEIglCpLhLJBJJECLFXSKRSIIQKe4SiUQShEhxl0gkkiBEirtEIpEEIf8PcC4Y2MAci08AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, (ax1, ax2) = plt.subplots(nrows=2, ncols=1, sharey=True)\n", "ax1.plot(science)\n", "ax2.plot(f)" ] }, { "cell_type": "code", "execution_count": 100, "metadata": {}, "outputs": [], "source": [ "np.save('pca_components', A)" ] }, { "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.6.5" } }, "nbformat": 4, "nbformat_minor": 2 }