jason-neal/eniric

View on GitHub
docs/Notebooks/Incremental_precision.ipynb

Summary

Maintainability
Test Coverage
{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Incremental RV and Quality\n",
    "\n",
    "The incremental functions split up the spectrum into small sections an calculate the precision on each.\n",
    "This enables the change in each precision/quality accross the whole wavelength range to be observed.\n",
    "These were inspired by the plots in Bouchy 2001 and Artigau 2018.\n",
    "\n",
    "Here we show how to use easily use the **incremental_quality** and **incremental_rv** functions.\n",
    "\n",
    "The split the specturm in sections with a width of X% the current wavelength. This is logarithmic and the section size increases with wavelength. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from eniric.precision import incremental_rv, incremental_quality\n",
    "from scripts.phoenix_precision import convolve_and_resample\n",
    "from eniric.snr_normalization import snr_constant_band\n",
    "from eniric.utilities import load_aces_spectrum"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Settings\n",
    "band = \"J\"  # analyse this band only (for speed)\n",
    "snr = 100   # SNR per pixel normalization at center of band.\n",
    "R = 100000  # Resolution\n",
    "vsini = 1   # Rotation speed per pixel\n",
    "smpl=3      # Sampling"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "wav1, flux1 = load_aces_spectrum([3900, 4.5, 0.0, 0])\n",
    "wav2, flux2 = load_aces_spectrum([3500, 4.5, 0.0, 0])\n",
    "wav3, flux3 = load_aces_spectrum([2800, 4.5, 0.0, 0])\n",
    "wav4, flux4 = load_aces_spectrum([2600, 4.5, 0.0, 0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Convolution and resampling\n",
    "wav1, flux1 = convolve_and_resample(\n",
    "    wav1, flux1, vsini=vsini, R=R, band=band, sampling=smpl\n",
    ")\n",
    "wav2, flux2 = convolve_and_resample(\n",
    "    wav2, flux2, vsini=vsini, R=R, band=band, sampling=smpl\n",
    ")\n",
    "wav3, flux3 = convolve_and_resample(\n",
    "    wav3, flux3, vsini=vsini, R=R, band=band, sampling=smpl\n",
    ")\n",
    "wav4, flux4 = convolve_and_resample(\n",
    "    wav4, flux4, vsini=vsini, R=R, band=band, sampling=smpl\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEVCAYAAAAPRfkLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnXd8VMX2wL+zm7LpnSSkkBB6F0IRkI40GyoIdhQFFHtXROzY9dme6LP/sBfUh2J5WFBRigoqIKiIKL1Jh2Tn98fN9ru7d2s22fl+PpDduXNnZu+dO+fMmTPnCiklCoVCoVAAmOq7AQqFQqGIHZRQUCgUCoUdJRQUCoVCYUcJBYVCoVDYUUJBoVAoFHaUUFAoFAqFnQYpFIQQTwshNgshfjSQt58QYqkQokYIcbLbsbOEEKvr/p0VuRYrFApFw6BBCgXgWWC4wbzrgLOB2c6JQohc4CagJ9ADuEkIkRO+JioUCkXDo0EKBSnl58B25zQhRJUQ4gMhxBIhxBdCiDZ1eddKKZcBVrdihgEfSSm3Syl3AB9hXNAoFApFoyShvhsQRmYBk6WUq4UQPYHHgEE+8pcAfzp9X1+XplAoFHFLoxAKQoh0oDfwmhDClpzs7zSdNBXzQ6FQxDWNQiigmcF2Sim7BHDOemCA0/dS4NMwtkmhUCgaHA1yTcEdKeU/wO9CiDEAQqOzn9PmAUcLIXLqFpiPrktTKBSKuKVBCgUhxEvA10BrIcR6IcS5wGnAuUKIH4CfgOPr8nYXQqwHxgBPCCF+ApBSbgduBRbV/bulLk2hUCjiFqFCZysUCoXCRoOcKSgUCoUiMjS4heb8/HxZUVFR381QKBSKBsWSJUu2SikL/OVrcEKhoqKCxYsX13czFAqFokEhhPjDSD5lPlIoFAqFHSUUFAqFQmFHCQWFQqFQ2FFCQaFQKBR2lFBQKBQKhR0lFBQKhUJhRwkFhUKhUNhpcPsUgqWm1sqYJ77mu3U7ARjRoYh/DhzmyzXbAMhPT+K0ns14ZP4aZk/siSXRzOFaKy0LM8hKSazPpisUCkXUaHCxj6qrq2Uwm9e63PIhO/cdDqrOy4e24v6PfuHYzk1594e/veZbO3MUANv2HKTbbR97zbfohiEUZGive/jmt22cMmuhz/q7lGVzVMt8pgyowpJgxmTSexWEQqFQeEcIsURKWe03X7wIhdbT3udgjfsbOcPP99OH0uWWjyJeD8BPNw/jUI2VWinpdccn1Fi938sMSwK7D9T4LO+6EW2Y1L+Kmlort89dwTNfrmVUx2L+u3yD4TZdNKgFO/Yd4sWF60hKMDGsfRFT+leRkmTGkmhiz4EamuWlsX3vIbJSEklKMGFWQk6hiDhKKLjx/vINTPm/pRFoUePCJMCHbGkQpCcnsOegQwB+euUAFqzZyrS3f3TJd9aRzehclk1hpoXctCTSkhIoz0t1yVNrlVilZMvug+SmJZGcYMIqUYJM0eAwKhTiZk1BYYyGLhAAF4EAMODeT3XzPff1H/C1oXAwAXPL8e0Z1bGYO+au5I2l673mq8hLZe22fS5pfVrksXbrPvYfrqUo08KUAVXM+2kjX/+6jc5l2Uw/ph3F2RYSTSa27j3IT3//Q9OsFBLNgvyMZCwJZpISTFitklopSTSbkFLi9Kpaj+8KhY24mSnMXb6BC9RMQaEIK19fN4jCDAsfrdjEpBeWANC5NIv89GQ+WbnZI/+0UW05XCvpXZVHy8J0UpMSsFolh2qt/PjXLjqVZlNjtXLwsJUaq2TN5j0UZVlISzZjtUJ2aiJJZhMmk0BKiZTYPysh5xs1U3DD2sCEn0LREDjyzv95pP2wfpfX/Lf9d0UkmxNWHhrXhUFtmpBh8e19eKjGyuK121n42zbW79jP0HaFXPbq9xw4bKUsN4XHT+vGp6s20zQ7hU6lWWzfe5jnvl5LfloSWSmJHNEsh27Ncti+5xArNvxDt2Y5JCeaAUgym0hJMkfh1zqIG6GgZIJCoQiES17+HoALBlRx9fA2AGz65wA97/gEgJ9vGcau/Yc9BOOb3/1l//zn9v0c8/CCsLSnc2kWM0/qRNvizLCU5424MR8drKml9bQPItAihULR2Pn9zpFc+8ZyXln8Z722o0lGMt/eMCSoc42aj+JmR3NygpmmWZb6boZCoWiAPP3l2noXCAB7D/p2Kw8HcSMUFAqFIlhufe/n+m4CAPl1m14jSVwJhWgZyrqWZ0epJoVCEU8kmiM/ZMeVUPi/iT2jUk9Sgu/LesngllFph0KhaFxEY89kXAmF5gXpYS3v6HaFuum+pPnvd47kiAYyk7AkxlX3UChiHrNJzRRiGm/2vWQfMwUhBAkGbmyTKNgO/fHv07vVW93NC9I4rWd5RMo+NULlKhSRJgrWo/jZpxAo5bmprNu+z2ceb968/sxHURD2YaF5fvAzq6yURHbtDzwq7cS+lRRnp3BCl6Y88flvQdfvi44lWREpV6GINOYo7NpuIMNT9Km1SkqyU3zm8XZ//C0G6d3YYzs3dfkerUXxm49rH5Fyg93/0qEki3P7VpKXHpmZ0uVDW1GRlxaRshsKr00+sr6boAiSaARiVELBC6GExUjyIhS+uX4wgKEYLbG+qTDUxfKyXH2BG05F6KiW+Sy8brBL2sWDW3JkVV74KmmAVBWk07nU+2zpxXN78t5FfaPYosbHUS3zI1KuEgr1SG0I4UITvZiPCjO1zXN6A597UrRkQrDCp7fTwDqpf/OAz+9SlqOb7k2g2lg+42jDdbxwbs+we2tMP6ad3zw/3GS8jTmprnF1elTk+sw/dWALw2XrUZKdQqYlgTlT+/Kfs/Q3t/Ztme8RSuGtC3pT4RZWPBBGdSwO+tyGyPPn9PCbZ3yPMvvnVIPxjUzKfFR/GJkpeF1T8DOwGbmtI/08RJX5vk0gtrfAhZPuFY6BvFszx+f+rQo88vq7et6E0dHti7yeM75Hud/gZL54Y4qn2aS5n+vozKwzunFO30rdY87rSIG8vrUoy3XG9ND4LlwwoMpr/iuHteaqYa1d0u4d05lf7xhJea7/QfvLaweRUNc/fb3Bz2wSPH9OD/59eje+unYQR5Tn8OlVA7k4iBni8+f04IFTugR8nh7vTu1LJx+znFjBiDVgQOsm9s8PGrw+lw1tFXSbjKKEAvoDQygzBX8LzXr9ZUBrx8BakZfK9GPbuQy87vRt4X16unhacLFRnGnmphW+cG4PXpvcm4fGdeGdqX1cOn2XMmMutv68iXpW5nqdHq+dOYo7T+zoklaZn+bRTg/qistPT6ZbM99auD9sAuuLqwd67Hn59vrBeqcYZmi7Qr6fPpTirBR78DXQF+7uAlWgDeKzz3Nt0+uTj+TeMZ291pngZxrVr1UBwzsU0dTP2po/jmqZT1KCiSObh262y0lL9OgHAC+d1yvksoPF39qjN4a2LeT6kW348eZhDG6r797ujq8xIVzErVD4v4k9efCULrx4bk8m6Gh/gcoE5+lfSqLnVHByf2ftz/VhrG6Ww4ldS3l2QncAynJT7S9G8Ya357k4y0K+wUXaL64e6PWYc9WlOSkc1VITWsd3KaFTqUMICAGpSTpObH6u39jqMo+06grXDu/t99uE7vwrB/DZVfq/4Y7RHQ21w9C0zY2y3FT6uAnl7NQk3rqgN8/5MBuMPqLEQ2GwVX/J4JZkpyb5rNcmeLz1TXfttCI/zafyoDfremdqH59tqMwP3ISkpzX3DnJdp9Yq7f3NWSHwt07UrVkOa2eO8lgLsz1zzrxyfmAC5pMr+vPshO48cuoRAZ1nMgnO71dFenKC4bWCaHgfxa1LqvND/feu/S7HRnUq5uzeFVxaFzrXO46ns1VhBt//uRPw1Jy/nz7U5YF3v69pyQl16a4HfI1n3qanMwx6EwmhDW7esEoZ0qJvQWYyu7e4Bu9yLq+fm8nppfN60aPSmCb/2VUD+GvHfq/HP7miP1V1GxVT667t4DZNdPPeObojp8xaaKhefxxR7luLG9+j3ON914Fc4yaZvgM66s0gipyCQA5x00bd+6kl0eQi8PU4oUsJS//YyQsLQ3tj3cndSvnq1226xzqWZPHLpt0e71SvzE+jKMvChp0HAq7Pdm0uG9qKXzbt5v0fNwJQU+u4Zr/eMRKTMGb6cSbRbLKbgu7+YBXtIhjaOhrvEYrbmYI3HhrXhUdP7Up3Pwt+zvRvVcDJ3Urt392lvrsG6O2+2tKNrP166xyZOtpfhsUh+3+/c6ShOqQ01g5vmsv/Tezp4WZrQ2/BsmVhumFtqTgrhWof96fKaed6enICX183iNtGd9DN29OASeO+MZ25xsmkExJhcCAw6hvgvij5lM7Csq/1Cz2EEEGbMJyb42twu2pYa92+MP/KASQnOGbhgfhIOGe9+XiH4uQ82zCbRFBvb3Nu6udXD+TfZ2ibPudc6HvWFQzReLtcRIWCEGK4EGKVEGKNEOJanePlQoj5QojvhBDLhBAjI9kef4zsWMTxXUoCOMNxg5xDXvi7b95ubIe6TVXnHqWZs2yd/oFTOrssrg1u0wQRgN3jw8v6ea07wWlRPD3ZITyK/YQZN5sEFw9qwdteOn5xVgoT+lR4pM8+ryevTe7tkvbMhO4+TV6DvGj5RinOSgk6kNik/s05qVspUwIcPJ29ipzNiXrXBAIb4Iy6S9tudV5aEhO9LJDXl+Ozr/7rbw0j1HGxSYaF3+8cyZJpQ2hZmBFaYXh/njv7WGub1C9wj71oETGhIIQwA48CI4B2wHghhLs/3zTgVSnlEcA44LFItccIurZxn2iPlHufEAgu9+El4OF+Wvc3Ny2JtTNHMbBuKmpLb5aX5qE5BeJqWZzl+ZBdOqSlvU4bzuakJ87o5vfhu/zo1nZB9s7UPtw/1nVRs6uOOaV3VT4FbiE8Brb2Pej3am5s1paenGDYtU+Ps3tXeKRlJAdnYZ17yVG66deOaGOfrUFwA5z7QG4rw11W2AarJTcOZZoBV1rj9RsTJX1b5HPTsYHX26JJeGKU3X1SJ68L0EIIQxskO5VmMbZaswIE4lXmjnPYmrHVpVw3sm3QZUWaSM4UegBrpJS/SSkPAS8Dx7vlkYDNAJcF/B3B9gTMxKP0tSt3LVrY/6v7LrRNUs+d00N30cp9IPC2oJqerA1wel4iocRVXztzFJcO8RRaztUEuqO4U2k2J3Yt9ZlnXPfIxhxaeuNQvp9ufI+AO22KMjw2HQU7XS/OSqFl3eDmajbRN1EYHWgDIdCmBzL7NMKLE3syoY/jGQqkPbZHYt6l/bzncbtmiWZHBef3a87Y7mX2hf1g9/28M7Wvfc9GMF3BFvwyVA+uaBJJoVACOL+qaH1dmjMzgNOFEOuBucBFEWyPV4a0LaQiL5XJbpuwnDu0DSGgfVOtk9hc0bqW57g8ULZP/VsV6NqsjT58D5zShauGtaZjSZbHGeO6e3rvhEokd0uunTnKPquIFEkJJr/uwHo4r32c43bPfQ0Ezm7EetjGoSisDeoSjY1OgRCM0CnSMWN6K2fmiZ3sn61W/Vl8MNiei2Cu50l1ipLzsxVu4RtuIikU9H65u7weDzwrpSwFRgIvCCE82iSEOF8IsVgIsXjLli1hb2huWhKfXjWQFk0Csy+2Lsrk48v7c+HAFh7aoC+M9q0mGZa6soXHhfO18SgQnEuJtUEkWru6U+pChAvhqX36eoCfndCD+VcO4PHTuvosPxqLg7r1GsgTq9FUgrlko48ose85CufPst0/5+cjUhF8ARZc491VPBpEUiisB5zV2VI8zUPnAq8CSCm/BiyAh2O1lHKWlLJaSlldUOBbO/PHmxf0ZtYZ4QkJLaWkRZN0TCbXoSMqbmNeKgnFDBGJdl89vLXPBTcj1Kdm5U/2VuanMSKEEA623+ZrcJ578VHcfbJDCzZ6NWJMxgeEEWHlnsdkEvaw6O7HQhEStj7g3BcC3VkcyK0ozQk+nEg4iOQ+hUVASyFEJfAX2kLyqW551gGDgWeFEG3RhEL4pwJO6C1+Bo7nLXb24vHXATzXFAKrXRKdNzDZ6wvhibpgQAsuGBBavJ5oIaWnS28oA6ttrSjUW9WuaSbtmnr3fffWRiMzv2gKjmDq0o0TZqAcm3JkzxpgJ57Yt5J3ftB0WJPOTMHorNpWq9H8RjeeRpKICQUpZY0QYiowDzADT0spfxJC3AIsllK+A1wBPCmEuAzt+p0tYzw8aJ1e55Hu7Jng13wUxDDhfka4TD3Ogc8asGLpshcjUK4Z3gYp4YQjSrAkmnn01K5cOHspgOHwA3rYe0mYL2ykBnKj5cbKE6rXDttzYTsWrOlu2jHt7B5bNgXMxUQcaIEG9mj8ePMwv6FHokFEdzRLKeeiLSA7p013+vwzEP4dHlHCvcMVZiaz6Z+DBvYphFhvaKe7UJGfxkvn9WL8k+HZ1RsIncuy+aFuF3gozL9yQEjugnnpydzjFCNoVKdiLpytfQ42ro0zvu6X3Z005Fo8CVRxiOZgX9+CLbAyPU18kVh/Sw/S/TncxEYrGgm2GUCg5qNwTauDxeYZEe0F0dkTe7J976GQy/EXMTaWieQVN2RmiRGt3xvBXh93g0MoP9O2a99lvS5gd9+GgxIKQeDtQXJMM8NvPooG0W5VWnKCPe6THjE+XvknRBOGN4z2n0BrjbROEA2lw1aDuztwKMLP9vpc15lCgO2KzUdeFyUUAsS5b7nfZ1unj4T5yP2BMjIwnNi1hMO1DX5obfAY09hD8Brz4sFUX66wsUQ4LoGeqSjQaxtr7t6+UEIhQJxvrbfH2N/tj9Za0v1jA3uxSQPqt42GDiVZ/LB+l9+w2XqU5qSwfsd+WheFHr8HGpaJwx/h/C32NQWntFCe4Vh/zpRQCBAtTIH/PH5yBFyvR1jkCKwppAQc+0nhCyM7mqcf246Tu5UGtC5iu/cndCnh/P7NdSPjBoPRuUpOEAIMghuo/T1Ld4zu6PV9386Esn/HJgCky5JCYK5aDWmmoEJnB4iU0q99MtCF5vqma3k2lw9txX0+3tJVH8TYZQoaXwNbcoLZ73sYfOEsEKLVr/yF9zCCe1DEYBACTu1Zbn8BFERmHcq+0Oz04AccV8olf2z3bCUUQsD91kqDWkE4NKZwT48vHtwyLA9qJIhEsLhoEq0hwGNNIcDzje+UDu6dCrYufN5RlfSu8v5GOAj9ngezBuevLFfzUWDlOUcjNiK6Xj6/F9NG1U8kVSUUwojdXBCyeSm2qK/megsJXR9tiCUcHjYNU1j6EwjO+Lr8vvqFp0tqOMxHQcwU6jKmJycwI4Aw4r2a5zHxqPp554ISCgGiF5wu4DLC1A5dIjBOxLove6xiD3MRoR3N9SosgzgnkPaGy207HOWYQpkpOP3oRHsE3xjUNJxQK4sh4K1fRHpHsyJ08tKS6N8qdNu4L5z0yojW46gvcClRL7MNIy66Btrl6zlymHxC/316+xSC2QPSUJQrJRTCiD3eip8uE47YR15LiMg2//CXaazeyFW85MahESs7HohWn4gFBcoR5iL4hWbX8kJtUWRR5qMQ8Cb5Iz1T8KlwNELzUX1WH4u7z/UWPvXzBVduQ8ZzsT3032TWud5Gr1UDmRy4oIRCGLFNVaPhk9wInl+/hCNEQX3iiNRZv+3wRSwKvXASFvNRJKMWxiBKKASIrylktAaBaD/Ggf6esMWEb9zjVUwQ7KDpTZjMntgzlOZ4radVYbo9LRglIRThZwpBJjTELqzWFCJA4ItQQawxRFH9DOQhnHNhn7C/pLyhul5GmlicQfVu4d3dNNjm/nLbiJBDw4T2RkLPNYXg2tAwUEIhQHwNxkb3KTRmQn31ZqwRjnsZ7u5gtE0RVRxCWWgNMH9SgqtBI5CfFR6XVO1vMIO63jmxPjwooRAE/jWGWL/tgRHPQi4UfGmnS28cal/ADFt9IaqijfE2h8clVbsy1hAj2Xav0HaBD20X/Nv8ooESCgESjjUFNcgaw8hL7RsCevfbNexBcDRUs1qorfa5kznEsvUIl+NIm6JM1s4cFZayIolaaA4Bz84ZOe+jeBQkXcqyAHy+tD4e6Vyqmei6hhBIT48BbZoYzhvp7mhEEfC9ea3ubzjNRw1TBgdM/MwUpIQN38PerbB/J6TmaH9/nQ/fvwjHPQxpTeDv7+CoywEBq/5LmdjBn9LLdM9aC1arfcujY/MaUFsD+3dAap5jS2QUaEz9dniHYhZcM5DSnNR6a0MoQ0qkBpE+LfJZdMMQjwCGoVR3eq9ybjzGeGyeUGho+o1emIvGTPwIhftaw55N3o+/c5Hj82cz7R+/SIYJh67imaR7OCTNJIla+AOwAK+5FrEEOJKHqXikqc+mlAJN+RejzAu5IXE2rANm+G7+a2h1vlrTnx5//w1zB5FJNYnUMMy8mDsS/6NlfFH//IXJORSJHT7qEYCkjZjJYcwwI4sy4P2kctruX6edl2CBmgPQeTz88JJnEWW9oOPJ8Pm9sGej9x8zaBosfBw6joEdf0C/K2Hbr7BxGbQ8Gg7+A/+7Hcp6UHr0bbDiPUi0wJ4t0GYkbF6p5ak9pOWXVkhIhtrDmiBOyYG9WyAlV0u3xRjYvUE7Jq2Q1HDf6wzhCT3tTEl2KskJZsP5I//qztDOtwlktU8hcOJHKPgSCH54JukeAE0g+OFry0V+8wB8Zbk4qLaMTfgMDgHfrmaZ5QnD5xWJHX5yaD3+g+RrXVLbmtY5vtQc0P7qCQSAPxdq//zxv9u0v9/8W/v7y/uOY18/4vi8ZQUsfc5/eRHg5aS2PFkzkoRbTw26jDtN3VlsLueyQ29oSoT7T8kohqrBmiBrd5wmtP78Fpr14WzzKn6XxbBoA+S10Ga4v87X+vHY5+DAP7D5JxAmOLgbmrQjaesOSsUWdsk0dpMCezZr/3Ir4dA+TSBu/hmySmHVXCg/0rPRh/aByay1yZINO36HJm21UdZ5pA5iGhSIS2csmWocMiGGGhVBRKi+t9GmurpaLl68OPATZ2SFvzEKhYIJh67i8cQHsYjDAZ33q7WYF2uHMNq8ACuCLqbfQmtIZilbyMSy63cyxH4ADmY1Z/KWk7g47SOOSN8FJd3gx9eh/WjoOQX++JJ/zfuBn6wVPDEkEdodDwkpkJQKmdqMf8WGfxjx0BeYTYJaqzZerp05iopr/2v/bMdqhT+/gexyqD3Eq99v4doPN3Fmj1JmnHhEwD/JVoeNUBaqhRBLpJTVfvPFjVC4vz38sz78DVIoFI2a9SOeo+9biZgECFnLr5YzAGhz4Bky2Mciy4URrX+ttZDDJPBc7dHcdtUVkF0WVDlKKLjz3YswJ7I3T6FQNE46HZjFm0kzaGH6u34bUtkPzno3qFONCoX4WVMQxhfRFAqFwplllvPruwka676JeBXxs08hHh39FQpF48Ic+qZHf8SPUMgsqe8WKBQKRWhIa8SriB+hUHlUfbdAoVAoQiMKFo/4EQoAGb43lSkUCkVMkxjesPR6xJdQEPH1cxUKRSMjwRLxKuJrlDz7vahUc7hpd5/Htw+c6ZnY5pgItUahUDQaTJH3oowvoZBbGdbi/m4+Rjf9YOsTvJ90/d/UNOnomZ4Q3lg2DRZzMlz9e2Tr6D4xsuUrFA2Y+BIK4caUqJssfcWBTEpD6l52t3PSY+BFHE27Rr/OEx6D1FyYsQt6Bxcfyi9FOkI5nuh0Sn23QBE0aqG5/khK95ulOMuLdu/HQ8DQHvJo7TTvONb7sRMeD75cS4CxpkwJ0HOyFnvGhtpbEhmG3VHfLWi8HH1bZMtv3j+y5RNhoSCEGC6EWCWEWCOEuNZLnrFCiJ+FED8JIWZHsj2BIbSgVr5z6Kf7ffVajMji4XdpAcK8kewmGG/cFrm2XPgtjLgLzE6zr6JOnvmu+tX1+yleYoWDFlbbnapBWsCzuEZAUob3wyc8DgOuj15zAuH4R+u7Bb7pfRGc9gZc4Cda8BinkLntTzRe/nCd9cgwE7HRSQhhBh4FRgDtgPFCiHZueVoC1wF9pJTtgUsj1Z6ACUVL9XbupC+0w0ZeuhNK/UNuNpYvIcDdkWanqChZTkG5znzHM2+gEx09V7uOJ7t+73Y2pOU7vh91BbT2EjVyyAw47TW4YpX23ZINN26FM97yLDcQWgz1b1YraGO8PPe8fhQRwFWQ95xs/3joKF29y5Xu52nmuStXwTkf6ufpcioMuMY17dyPA5/92ehxPly+Irhz3Wk+kNqUfP/56pOWQ7Rw496wZEP7ExzPTU4z42U3cJfUHsAaKeVvUspDwMvA8W55zgMelVLuAJBSbo5gewIkAqaLYpvma6DsihA22+U2936st9P7Hrqc7j1fToXvOs6cE1CTACjv7fg84QPXY/7MZTN2wbEPaZ+z6x6iwdO9v9XO/bclJDtmIc4eHPmtfNfrzOQFcPrrcP58aDXcR74vjZfp7mJYNcj1+0VLXb9P2wzJTq8nHXGXdm1m7CKpq4F3P4y6V1M4ktKgvKexNp71HpR11wRsrwuMnePMyHsgsykrknVmfoFirWFHb51ZzFFXQJP2xstpPUrbt3Ttn57Hxnl5X4g3bIqHUVoM1v427w/jX4GBN8DUIIJ8RohICoUSwPmKr69Lc6YV0EoI8aUQYqEQQvdJE0KcL4RYLIRYvGXLlvC3NFnnHcC9pxKsYLCm5Pk8rrvQXNzZ8TkxTXs9aKkP11ZfQiOvyvuxQTc6PvuaKVitrmYum8vsjF1w007/QkNvquBcXzP3F7wEMLWY9Dlc/J3345evhPQC7XN6IXSboM0aQsV5gfrEJ73nMwcQZ9J9RjjiblcNPq9K09JtJCRj+Fpds1Yzy/lgZ1Kx94M3boXp2x3RABJT9J8Vozg3O1ilR9ZyoLSvZ/rg6dDvCu/nlXSDaU5jx/jZcMUKsDj9nqNv12bzbUbCpct932NnMoo0ZWvE3a7pl/3k5QSne956uKas5LeEUfcFHQE1nERSKOiNqO69OQFoCQwAxgNPCSGyPU6ScpaUslpKWV1QUBCe1g2erk2lxzwHfdy8XGbsgv5XB120NcfN9bX3Ra6aq9tAsNOcC0dO1WyRoGlwSd7fS7wspYf+W7NAewVlYXutU+vh7jHlzUwlrdrJ+qF2AAAgAElEQVSbwQZPh4u/h3H/5/0cvTIq+3mm5dYJK70ZSkoAL6FPyfY9G8p0GuiEgGMfdBW64cDiNjgmZWiDiM7ivNXnGkbdtZvwAdywURv03TX4MjfloKyXsTaak6CgteP7yHv9n+OssZoTdfzig3eAkM79pIuPWc2o+/U3afW+CHIq7bPKjcJtLHC3zbdzcg0XJk0pGX6X9+vQe6pjNp9drr0u1ujM6OjboOck17SsUmPn2ug+Uf+5iTKRDJ29HnB+G0Qp4B6MfD2wUEp5GPhdCLEKTUgsimC7NI5y0ioW/cfxecYux+fuE+EjJ83ag7oHpMVQOPICeGG09t19IdnNI8HdZfXPxEqyTSan9/5J1796eBvMWw7V/nqzTbuf5/yu4sp+8Pvn2ufB07W8R3nRvpx/Y1OdN0qd9BSseh9en+Ba9w2bPCM9Ol/zaDPoRnj1jNDKuGipZidO8zJDtF3zITfDxzfp5zEn+rYXT/wEuwDxOhNx6y+2e3T5SkjO8HQcAFfVLTFN01h90XkcLPy3NsDuDWXW7mMW3v1cmHeDZ7rHc+RepKNMKUwIvevZa7JnmtcmChh+J7uy25P1gdu7WIQp+OB0seJo4oVItm4R0FIIUSmESALGAe4rkm8DAwGEEPlo5qQQ38kXBLab5O6/7j6D8DzR8bHAeWHJpGnsXqvzctltJpmqgXUJvrQyLw9V53GOz3oake3BsS2W2lxvCzto9k17OX582YWw27JJzvAUHokp0EHHqyLR4n0dQI/LfoIrfjGeP1DaHec/zyU/+PYmyavyFAjOmq7tNpb18BSIJz2lzZyKu/huQ2k1lNoWmA2aNW39OrNYXyAAae1HOuU3UG5uc7hunaax1xv+ZytCWl0VoxBcvK16MdMCHdibD4Tqc+vOjW1X64gJBSllDTAVmAesAF6VUv4khLhFCGF7EucB24QQPwPzgauklBH0e/SC7QaHYmJw6iRSmOAaH7ty3TqF/VtelbZoZdu0VWZwIdAZ587f4zztr7vb6XnzNS8cZ3IrfZqs/DJ4umbDdiclF/pdrc2mjrrS9djwu/x722SVQkY9buQb95ImrH15k+hxwdeOz873291NNq8KTng0sHUI26DYaZzvbAaER6IlzW8e/aINDmxOXmotCpzq8nt+3W88+vbA2uVMv6ug72XBn2/DpNNWo0LBpnzlNHNaI4xToQAgpZwrpWwlpaySUt5elzZdSvlO3WcppbxcStlOStlRSvlyJNvjlbQ622SGj0U3Z2yeJzb3yILWrp3E9vnK1bohG9wXml3MSRlFjgdm6C0w5SuY+Am/J7gtHnt7qFLdNNbLV3guXpV01ezy4UZvXeCa32HQDZrXTqbb9e012fvaR0AIzZ5syYbkIN0mveHuDWSU3ObswNYWp3t10n+0mUc4yKrPd4R46X/OXmkXfweTv7B/TUpwXp8wODBWT9BJdDrXku22luB0zJzo5LIcymZQ7XmtlU5lnzjL2KlHnFbXLCdzkz+BmFGsOUcE4sUWRuLndZy+aD0CxjwbeFC6pl01X+PyI6FmvyPdJhTSm+ieph/mQgdzorZoDBwSbmaHI06HT+90TTt7LjR1M0NkRjlc+LA7gpvhhMKMndrf2prgzm8+AH77FPpcCltXw6r/Oo6FY6rvXEaixYDnlh9sSkxaGJwuwm3KaD7A8dmXM0AgXPw9JDrPYp0G+Gv/cM0rhOtsOdz3z0brkZ5periYrWyf/bTpipXGyo4Q8ScUzvkQ9mxyTRMC2o8OrjzbtnNnd0u/O5oD76jCXdNJL/LMVNEn4HLDzpEX+s8TKQIywThhM3HkVsLQm2HTT/C4bU+Fj3t12uuum+ncEYD0EwsrGI44U9uV7b4Jz91uHknbdahlGz5f+AhkqVeGN2+6EGYKeqaiYN75bmtDjK8pxJ9QMLphxwhuN3cbWeSxS1tTMH5aWOr2GbZA4Qc3Da7Q4CYom6eXP/RueF5LT1OfUUwm/44AQRFIx4yFgU1noO9wIix/DXrY3EOF97yGCWFNwQWDM4V6Jv6EQjjwonXYtfmIuJy5dySn7zds1DmuMExxF/juxfCZOzzQuTcXRWIHa5SCKIaFCPXX1LqZWyChI/wg6haapW3qB4F50GmlOM0Uwta0iKCEQki4exFpN134EQpSBt4rfPlkhxQPxX1vRLww8l7Hzu/uE6FZHyhs5/ucBkcMm4+yy/zn8YKjp/rQ4MPan8Ok5NliR6XprzXGCkoohBFTRGcKboTLLtlqhLZoPXBaeMprKNjcdUG7lo1OIBjEefBs5x6azBch9L+ht0K5wV3ZgWJXcmyePmEsM1TaHa9Fee2o/3KuWEEJhTAigrQZGluIjJDWl5AU++GI65NIea/EGkdO1SLLRoPi0ALjCUM7/d3yhHuhOahyhKaAxTixvd+6gWETClFZaFY0IKJ0w4PxPrLlSct3fZeF0fMiTcD1uM0UwtKEOFEM6lBCISj0tQ7DrofuO5oD7C9hd3FURBQZywNCNNeSwnQdfLY4EmsKocwUGuBanRIK9UCo/cRjz4IixolhoWAnwDaGMsCHaaDUVY7c1xTCs6gQhjIaDkooNECyUgJ8Y5oiBIIfEIJdYwofkaw3RgdKu1YfIzOFBohaaA4Fb9pSBJ4XZxNEeV4QQcxO+k9gNmNF+Ihl81Gw1OtvMhA92GM2EspCc3zNNuJLBIYLL9Nfo7b+4PqYT+9s/3Q8OUCXQ0Wo2PrDzo7naAn5rX3kbmgE0wsDOKf5wLpTPMNJ2IKW6j5HEVhTsC00x8tanpophES4Fs4CdUlVawoNiT1Vx5A/4rrIVxTN2EeR4Jo/HGsBY56BXX/pvjK2abYWHK8gPdmzDI81hTCgzEeK0PH9MDZAh4T4pdG7I8ZQZ3QO556YAvktdLPZ1moSfL7nIIwzhZCUvxi6vgYxJAKFEB7bPYUQA8LemkZC5LtBLA8yCncahNkhUMEVC4JO336k/Ym1mUIsXC+DGP21rwohrhEaKUKIh4E7/Z6lMESTDJ1psE8anvYRj9S791EgA1Gg09e2xwaWPxLotTki+xRCuH+2iLv2t67FPkaFQk+gDPgK7d3LfwMxELy/vvDd4fx1IfcdkkWZFi85FY2CmNYSg2xbtwnhbUa4sEe5CN9MQQYcEdWJir5w6Y/QaWzY2hNpjP7aw8B+IAWwAL9LGc75WQPF68Pu70ELVYuJ5UFG4U70zEfB9Ksg+6IQYArWxTlMWrze8+e+pmCuW6hOzgylohDOJaSIsPWBUaGwCE0odAf6AuOFEK9HrFWxji3CY6brO3KNdvVgZrYNwi7dKGnsC811BNXGADuytzosAb4v3OcD5LamUNgejr4dTn46sDqcS2wI9y+MGPU+OldKaXsryEbgeCHEGRFqU+zT9wrtZeG2ePyB4t7J4qzTxQsmYQIJloQouTRG+33cgXLEGdq7sAvauKZftAT2bg2iQAPvUxACek8NomydMuMEo792sxCi3Pkf8FkkGxbTmEy+BYKfQd4jdpGBqUNasvIebmhkJGsbr2x+9QEx9BbILA3snKQ0mLErsHM61L3nudXwwM4DAp5FdTxZa5+78ErLhyZt9M/xid5Cc4xGSQ0Fc3TD2hgdaf6LdgcE2ppCJbAKMPgyW0WolOakwO4oVdYkTl84E2bsY0kwg0qfS7R/kaZpl8AFiZ0Y9ILT2QEdepma7lwvv/asdyGnEh7sELUqDQkFKWVH5+9CiK7AJC/Z4xbjdv/AzUemaK0pXPaT47WBihBNe8osGFl0rm/1ObB1FfS7Kny1OPeBBIun+SuSVPaLXl11BGWTkFIuFUI0HMfbGCMojSNaU9isAE0WCu/Ut9kh4sTg70tKheMeDm+ZzmsK0zaFt+wYxJBQEEJc7vTVBHQFtkSkRY2CGHxYFPVAY+8HMWg+CoYB10NyhvfjjV64u2J0puB8xWrQ1hjeCH9z4oP46mJxTJwNJg2WAdf4Ph5n3kdG1xRujnRDFAGgBpvoEWex9BV6xNf98ykUhBDv4mOOKKU8LuwtigtCnHarMKsKRdSod5fUKONvpnBvVFrR6PA9aAc3psdXx2wc2APx1GsrGi3RGqyV+ciF36WU66LSkkaBwU7q0ZnVgN8osW+kUkIhrNg2c2UUR6U6UScUTHEi3P2JwLdtH4QQamE5TIQsAuJsOttwUTOFiJBdBqNnwdgXolShqPs/Pu6jv5mC8+jTPJINUQSA0jwV8U7nU6JXV5zpYP5mCtLLZ4Uv1KCtUDQa7FbAOJEO/mYKnYUQ/6DJypS6z9R9l1LKUIKUxy0h72hW5qOGQWNfU7Bkwb5t9d2KqBEv5iOfMwUppVlKmSmlzJBSJtR9tn33KxCEEMOFEKuEEGuEENf6yHeyEEIKIaqD+REKRWzSyNcUzv2ovlsQJeJLCYuYr5UQwgw8CowA2qG9mMcj/KYQIgO4GPgmUm2JFkYf/ZC7WGH0IiYqQqCxzxTyqqIbHK6eiRfzUSQdcHsAa6SUv0kpDwEvA8fr5LsVuBs4EMG2RJWAh4BAzEGDb4IB1wVag6JeaOQzhTjB9ngq81HolAB/On1fX5dmRwhxBFAmpXzPV0FCiPOFEIuFEIu3bIndOHxR0SNKuoFZvXBHoYgacbaGF0mhoHcl7aJWaDtCHgCu8FeQlHKWlLJaSlldUFAQxiZGBn9dKD70DUV8EF8DZjwQSaGwHihz+l4K/O30PQPoAHwqhFgL9ALeaciLzUEN9u1Hw4i7w90URSxgtx41ZjWgMf82DRHRYTL2iKQdYhHQUghRCfwFjANOtR2UUu4C8m3fhRCfAldKKRdHsE0RJaj3ro15NsBaGv9D2HhQawqKhkfERKCUsgaYCswDVgCvSil/EkLcIoRolNFVLQna5UxNTvSdMRgbZZzZNRsF9eV91HNydOsLM1v2bWHZlmVhLXPDng3UWGuCOtffo7d6x2oWbVxk/75t/zbW/bOOpZuWcqDmAFJKaq219uNb92+1fz9sPcyOAzvsx3Yf2s2Bmvr1uYnoiqWUci4w1y1tupe8AyLZlmiQkmSGGodw8EqjNifEPlv2baEgNfbXpoJmxF3avxhh6OtD2bh3I88Of5YkUxIdCzrq5luxbQUmYeLkd08G4LyO53Fx14t1827et5lDtYcozXC8PvaXHb8AUJlVyW87f6NVTive++09nvnpGVbvWA3AtT2u5fvN37P9wHbu7nc3b65+k4kdJ3LYepiXVr7EqW1P5anlT/HY949xebfL6d20N9mmFAqd6t57eC87DuzAVvOJ75wIwJDyIew+vJtvNvj3rh/XehyDygdx/kfnA/DE0Cd4aOlD/LztZwBObXMqiaZEJnacSLYl22954US5sYQVY9p8UCKh+hz4/XMoaBvM2bo8/N3DtMhuwYjKEUGXsXHvRtIT00lPSren7Tu8D0uCBVOAIYcP1x7GbDLbz/t+8/d0yO9Arayl+sVq7u1/L8MqhgEw/cvpHLYe5s6j7nQpY/66+cz/cz639LnFJX3HgR3sq9nHhA8msGHvBq7ufjVntDvDcNsO1h7kse8fY3LnyaQkpHgcl1KGHHd/6/6tAOSn5PPDlh9YtmUZX//9NT9u/ZH/jf0fCSbvj+tFn1zE7sO7eXb4s/a0WmstTy1/itPanma/P8e9fRy/7/odgHM6nENVdhXHVTkm7h2f0wbs5Wct50DNARZvWkzfkr729Dv63sExzY+x/1ZrXW+23enDtYeZVpDHpdt3Ugws2bSEjXs3AnD2B2cDMKxiGOv+WcfdCQlU1Di097HvjXX5TU8uf5L+Zf1pl9cOszBzoOYA076cxkd/ODbN9W7am/sH3M+Vn13Jgr8W+L3GM7+daf98wpwT2HlwJ8/8+AyTOk/i3sX3cu9ix9sC7l9yP/cvuR+AW9LT6Lv3AFu3rbC384kUC5OKmtjzf7zuY7/123h51cu8vOpl+/dJH01yOT575WwAnvv5ORafvpgfk5NZkGrhzAM7yLHkGK4nGIRsYFprdXW1XLw4+GWHQ7WHMAkTCaYEDtYe5I1f3mDXoV0cX3U8TdObep4wI8sj6eCNmxEIkupC+G7au4lfdvzCt+9M5PLNG9l3yTJeWPc+j37/KLf2uZVuhd3YX7OfXQd3cc68c3i8+n76vnYye4UgYdomtu/fjkRiEib+2vMX7fLasWbHGhb8tYDBzQbTKqeV/o+xtW3GLq+/V0rJsq3L6JTfyWPQch4ApJT867t/cWzVsTTPas4vO37h1VWvcn3P6+2DtJSSNTvX0DKnJYs3avdgwrwJALx7wrtUZFXw5LIn+dd3/6K6sJp/DfoXc3+by+M/PM7Lx7zMgZoDnDb3NBaMW0CtrGXRxkU8sOQBBpcPpmthV86Zd473G1fHuye8y7FvH2v//vqxrzNv7TzW7V5HoimR937TvJufGfYML696mXlr53kta/lZy71X5HRtdx3cRd+X+wIwpfMUOhV0YsrHUyhJL2HW0FmMemuUvS2b9m2iX2k/7dyHusCO3+GipdpGLycWbVzk8ntv73s7Nyy4AYDPTvmM/q/099vebzd8y7kfnsuyM5fR6flOAFxZfSVntT+L2Stmc+e3msDMTMpkROUIxrYey0nvnORR7lntzmLqEVOZt3Ye076cBsDwiuF8sPYDAJqkNGHz/s32/H1L+rLr4C5a57bm9V9eB+DxjZuZ4jRAArx27GtM+XiKXdjp8dXaP9l52XJ2HNzBaXNP85ov0rTIbsGanWvqrX6jlGWUMffEuf4z6iCEWCKl9OvIE3dCwTYQ+mJc63EcU3UMZRllyLcvIG/l+wAcAt5NT2NGQR4AR5UcxRd/feFy7r82beHiQv+miW/W/knPijJMwoRVWn3mtQ0Gl396OQPLBnJsVd2gWDdwPT/mYTbu28jV3a/2OHfOmjlM+3Ia9/bXNKBuhd3IT8nnq7++YtLHmnby5fgvOVhzkEGvDQK0h3nMu2MAmD1ytn26f/zbx/Pbrt+4oMsFPPb9Yx51fXTyRwx9fajf3x4rfDLmEwpSCvQ1fCehcN0X19mFTVFakV3zBajIrGDtP2tdTr2g8wX0K+1H+9mnw/bfWDfhPQ5mFnHiOyfy/onvI4Rg+BvDvbbrjHZn8MLPnmGhB5YNpGdxT2Z+O5M5J8zh+Lf19oLq98tYp21uW1ZsX1HfzYh50hLTWHjqwqDOVULBjR0HdtDvlX5B1bn893VIoLpZGYdM0V/wPaX1KUzrNc1Fs1+0cRG5/xlO+eEaulaW29N/2/Ub2/ZvY/WO1YxqPsqu4fYr7cfn6z8HYMG4BfZ0gGk9p3HbN7d5rf/Z4c/SrbCbIYHaENGdMTgJhWB/9/JdiWzZ9TuDykv9Z1YoDJCakMo3pwUXEcioUIibNYXRc0aHdP68tNR6EQgAr6x6BWfhveCvBUz5eAqUupq79hza46I9zv3dMc20CQTARSAAPgUCOGzBjZWv/vqK3iW9XRPzW8HWX7jr2xAWbIXgp6Tk0BqnUDixr2ZfxOuIG6Gw7UBoIX6vapLvP1MEefWXV+2fp3w8RTfPgVpXV7YftvwQ0TY1FiZ9PIkrq68kLTGNtrltufaLa7nz+AfpYDXx4qfnB13uW4m1JNWTIqFQBEvcCIVQ2GZqGDsaH/nukfpuQoPF2esEYPz/JvPgwAdDKnN6Si2k1K8yoWhc5FpyI15Hwxjt6pkBzRqGTfiN1eo12uHk0vmX1ncTFAoXzMIc8TqUUFAoFIoGQqB7f4KqI+I1KBQKhSIshLpB0ghKKCgUioiSfEhSsNO46/up82t59c4azLXhcZcXUvLyzBq6rfa9H8hG0mHJ4O+thsPRXPhuLf9+uAYhJSar5Mifrdz6fA0P/buG896vpWKjpOdKKykHJKVbJD1WWUnbr5VdsFPS+2dHu/r8ZKViY/1uE1ALzQqFIqK8cJ8W/O2MK8wcTMRrhLnCHRIJnLBQGxQf/nctF1yoP0QV7JQk1sLfeY6yirdp521wShu61Mp587RB95rXrfxWaGV/MjTdBjecZeb4hVaeHmrCJGHgMsknnQUv3qu1d9L7cMt4EzVmwcpSR7tzdkuabpf81EzTqfv/qNX7ysxatqdD7h5HO4t3SIZ+7wiGZ+PnMnj4ODOPPq4d+6OJ4IEnXfPtSINJF7v+fhGF91cooaCIK0q2Sh54spZbxpv4sSKAibKU9PhFsqilQEbRzTRntyTlkOvgp0fzDRIh4demvvOlHJCUbYXbXqjl3EvM7E51zV+yVbI7Bf5Jc02/8N1a1hQL5lX7v2bmWkmtCRCCss0OrfeF+2r5OxeuO9tM2gHYmuVax8P/dh0U8//Rfv+ODEe+9n9YuWm2Q7O+c4yJH5oLrCbBQ7O089/oLTjpK8nMk012gWCj+SbH56nv1tLuT1ifJzFJmPCxleFu+2Knv6SdP7+T4PU+Jtr+KZn6npZ2/ZmCO553bbOzQPBFuz/h8Ucd57oLBICcvXDvkzVcOdEc1SjJcbOjOdq7cTP3Skq3Sn5u5niIzLWSkYsk71cLahLq33/9/lk1/NhM8PSw4D0a3AcRISWpB2CvhYA7srlWYjWBdDtv5tM1NN8Ezw42MbeHieYbJJfMqWVZpeA/bm2/49kaWmyAsde56juZeyX7k+H/7nE8fO55fHH0EisTP7TyzBAT71cL0g7A3pTw3ENhlUgBJ30pSaqRvDTATNp+yaT3rfRapT2fZ1xh5mCS9/pevVMLLHfeRWZ2pWv5svZKHnyilltONfN7kXDJ58w5l5jZk+p6fOx1CQ7ziRD29EUtBQ+cYLL336FLrawpFvxerH0v3iZ5aFYt73UXPD/ETM+VVq54y7vZxnYPEmuky72xsaRK8NERgm2Zgr4/Wzl+YWTGK5sgiXXWjerCsPteCupctaPZD0JKyjdD1QZJ8mFY0F54aE2BkLVH0nGtZEEHTQjc8VwtTXbBreNgTbHAJOHUT60M/V6Su0fw3BDPgThrjyShFqSAlENwIBG2ZRIZLUFKSrdB6TbJ01rgUSo2Sv7Og0OJAqQ2PT/sRXh1XW1FCrjuNSsHEuHMK7Wu9MpM7cGe30nw+Cj/wqZgp6ZV7kqDl+7Wzl3eTFC8Q/JTueCnZsKu3Z39iZV53QQzn9XyFe+Q/GcYpO2XFO2AtYXQYoOWt+k2yaAfrPReIcn/J9iLpJG9RzLxQ21gy9stGf+ZZPTXku+aC2aOMXHDK1Y6rZVMmmp20Wr98fLMGg6bIbkGNmZD0U4tfW615MmHXQfIF+6r9RBiwipJO6D1Fxs9fpF81FXQcr3kmEVW0g7CFW/Wcs9JZpd8zlz5Zi0zTk9wsaGP/spKj1VWqjbCzJMdik331ZIHZ9Xy8HFmMvdJuyZ++pVm+vwsmTJX+37MIkmb9TX2++ENk1Xr84U79Y93+1XS7dfID9YNQSAAlH8R+aB9cScUznu/luy9ULZF2h9CgAkfw+t9BCd+JTntKjO1Zu8Pd9kWbSCzTeknfFjLiCVap1pRJrjqDU0gANz4sqeWNGqR5Lkhmla2pghabPTI4kIgGq1Rcnc7PifUSBJr4O5ntIHolGvMnPOhlWHfSRcNtcNaK9NfsnLxJDPXvu74XZbDMPg7K58c4Rg8Bi6TfNXWyg2vaPkummSmwzpN+7VpppUbJHfVDfBft3Fc745/aNey/4/Sbq+1cesLroPlyV9YGbvA8xpPfafW73U1ypErHG0w18Jx32jfj/hNMmiZpNNa7fvwJVbmVpvsmrovBn1vxSQ1gQC49EV/JoiCnZKtmfDvR2rJ2QtvHemo77x5VjbmuPa7Jrvgnqdr+bSjfrva/ek5gxj/meN853ttK8/9Ptjs8M74EwgAL9/leZ7CB1HYSBt35iO96bMep15l9mriMVqGL2YNN3H+B8a8IWwa6IwXa/imtYn3u7t2jIE/WMnYD+/08uwwrf+U3PpiLVdPMNPyb8nSKsG2LMElb9fSp26wu3CKmV6rJGf8T2vP120ER67Uji2pEtw1VtP4Zz5TQ/ONsDkLu9Bz5tW+Jt0BOlZxNrX44p6nami2Rf/Yhhwo3uGadsMZZn4rxq5YZOyTVGzStOirzjGzz+KYUemxpIWg2xrP5/L9boLdKYKxC6z80hRa/a1zsqJxk5NN26+/DupUo+ajuHJJLdlqXADOvqeWlIOu+ZttkmERCIBhgQAwrk5ra/enthjmzpS5Vk6fr19er5Vaes9VmhfG44/VgpR2gQDarMEmEAC7QABcpu7N6zRvPYEANCiBAHDmJ1aeetD//fQmEMBTIADc/kKt3RQG8J+HarnxZSv5/8AzD9Zy6qe+r5OeQAAYsUTar7ESCHFKFJwc4kYonF5ygu4Kvy+eu7+Wyg2OB/Sep+tnqjtwuesgkbZfcuxCK91XeQ4u7dda6fSbleGLtWNpB7X0Zo53pHh4ZNz6YnxO4Y/6WZK5P7J1uCsWQMQWSxWNH7Mp8mEu4mZNYei9wb105K5ntQW+cG2kCZZeKxwDee4eOKNuZrCi1JGeckBy00uO70taCAbUCZRqJ+1z6HeB/ZbbnqvhllMj3xnrCyGlh8eTjbxdod33tuuUAFCED+uevRGvI26EQsraTf4z+WBIgANpuLn8bcdgf99TDs2+7XpHHpNbE28O0wyg1d/6C4mNhVdm1nL+RWYOJcA+i6twePyx0H63WcmEBoPMzKR28iRkWVlUFnSDQghWrPD9hjqLxUJpaSmJiYlBVRE3QiFUjv4u9u3laa6vUyB/t34+hSez6lxAbx1n4saXrZx1mZn9ltDst91XWT0EtSJ2qZ08ibwuXchOTIxKjKFgECYzlrZtvR6XUrJt2zbWr19PZWVlUHXEqDiMPcq8v3s8Zmi9Xo1AoWJz5XzugVoy9oV2Pa9604op9nUJRR2yrCymBQKAvygXQgjy8vI4cOCA74w+UELBAPW9nmAUPc8kRfDc9UzoJrPL5nfZwu4AAB/KSURBVKh70mAwmWJbIBgk1N+ghIIBKkJbjoga6cErBwodQt0JrVCEHRU6OzYoDWB/g0KhUARLaseOnHvddfbvNTU1lPfrx4kXXghoawYXX3wxLVq0oFOnTixdujTsbVBCwQAX/leZABQKReRJS0nhpzVr2F+3JvDJ119T3KSJ/fi8zz9n9erVrF69mlmzZjFlypSwt0F5HykUCoUbdyzezsrth8JaZpvcJK6vzvWb7+i+ffng888ZffTRvPb++4wdMYIv62YEcxcu5Mwzz0QIQa9evdi5cycbNmyguLg4bO1UMwWFQqGIIcYMH85rH3zAgYMHWf7LL3Tv1Ml+7K9NmygrK7N/Ly0t5a+//gpr/WqmoFDEISX/eoi/Lr7EIz2pWTMO/fFHPbQotri+OhdTSgrW/RGOg6JDx9at+eOvv3h17lyGHXWUyzG9AKbh9phSMwVFzGFKS6vvJjR6TCkpuulN774rpHIzjj46pPNDIaGgwOV71gkn1FNLQmfUgAFcf999jB0xwiW9tKSEP//80/59/fr1NG3aNKx1K6GgiDmSgtyJaZTcs86MaPkNgdQePXTTLZ06kTlyBObsbJf08uefM1Ru6b8eCrltwdLyi89dE8za8FZ0880Bl5VcVRWOJgXNmaNHc93kyXRo1col/dhRo3j++eeRUrJw4UKysrLCup4ASigoYpSU6m4RKzupef0+8KGQ0rlzSOdb6uzTwqwf4FAIQcn999Nq4de0+OxTe3piYaFLvryJ54bUjkApvHEaIshYPoGGmzZnZnqdSUWL0qIiLjz9dI/0kSNG0Lx5c1q0aMF5553HY489Fva6lVCIMsJiQZSEV7IHS/HMOyNSbpMrrwjpfJGcTPmTT5J9yil+81o6d/J5POv441zbds01ZB17TEjtq08krjbltitXkOjFfJDavbtHWsXLL9H6++8QCf6XE90FgTMiNZWiGTf5LcNGk2uvIalFYMI4wan+3NNOo83yZV7ztlq8KKCynWn24gvgdD0Sw6x565GQl6ebvuXbbz3S+nXvzpuPPgpoQvvRRx/l119/Zfny5VRX+31nTsAooRBlksrKMFeW13czNGoD23+ROWqUbrr74JM3cWLATUnr3ZuEpsVkDB9Oyb33YEpJoWj6jR75EpoWU3TrLZgL8l3SzTk5uuU2ueoq++e2K1eQN+FsTKmpAbevPkg9spfL98LrrgWdfZQ2QZF/gavPeunjj2HKyHBJEyYTJotFtz5v91erxLPinHHj2FxgTHtPKi2l6r33/OYz5zvua7PnnvU4XnzHHfbPueec4zgvPd0jb97EiSS1qCJjyBBH2qRJVH38EeXPPgNA1oknklpdTdsfl9vzuM9ITPmOAdzsdj0BEouK7J+TysqwtGnjV7AEK3hEUlJQ5wVCRIWCEGK4EGKVEGKNEOJaneOXCyF+FkIsE0J8IoRoFsn2xApJp4+p7yYAIK2BxfZpeucdnolC0OyF5z2S3Qcjf5Q+8jAt//c/Sh98wP7AOJs4ymY9Qf4FF9B8zhxyxoyh7OGHtQMSSh64n4rXXrPnTWrm2o3arlxB25X64YZbL11iqH0JPrRmI1S+M8dvHnOuqw97/mTXQT73rLN0B2cbWSee5PJdCEFLJxOQO0W33Ezzd9+xfy+5716PPM3+70Uq337LYzAK1uMlrXdvr8eqPvrQ5XtSRYVHnuwTR9Ni/v/ImzKZJldd6bUsU3o6yZWVVL33HglOCkOTyy4lqbSUtF69aL3sB4pvu9V7GSma8uCsRLj3LYAEJ0FmzspCJCR4nQkA9lAVwQzw0YjNFDGhIIQwA48CI4B2wHghRDu3bN8B1VLKTsDrwN2Rak8sIfI9O0yC067FqOFnplDx6isu3/U6saV9+5Cb0XblCq/ae8Gll1I4/UbS+/Wj4OKL7JqabZBO79uXzBEjSCotcZzkHAvfz0Nk5MFsvewHsk8+Wf98g3Zui9uCYVqfPgAupp+Ciy92PclJAORfNFWrLzmAgcRs9jkryhk7luSWLX0Wkdqtm6b5Nm1K8cw7yR5jXKFJ7dnT0ZSsLABMWZke+Zq/P5cWn84nqayM9AH9tfw+BtXE4mKaXHKJzwGy8u23Xb6XPfUU2ePHuaSZkpIQPt6bkFBUSHJVlb2P2O6GOSvLw9NJJCRgcputeLu2iSUluumxQiRnCj2ANVLK36SUh4CXgeOdM0gp50sp99V9XQiURqoxWSedGNbymt59l4dmZxidzqxn/w2W4ttvN9aMxERdbcxGSiff9vpms2dT/p+n9A/qaLRpvXtT8sD9htpmI3/yJHJPPdUjPbG4mBaffUb+1As9jtke2LQ+fXxrbAYxJSW51JPQ1DH1L3/macPlOK+1lD78L9o4mSwAsk8Za78f5c8+g8mSbD+Wf8EFAJTc66nN2xACKuc4BkPbgOecpkfxnXcackXNPuEE+0CfZMA7p/TRR2j2wvOUPvqIz/6dXFlpN8EU33orOaedRuWbb/otX4+K116j+PbbXZUEIL1vH4pv8r0GYsrMdBHSwmTSXXBOKivzWG+xtGlDstuzZEpOxh1hNpPg5tkVa0RSKJQAfzp9X1+X5o1zgff1DgghzhdCLBZCLN6yxcdb1H0QDo3Wmazjjgt6wDEyBcydMMHvIqo3TJnGTDdZxx9H9tixAZVd9uQs++fUrkfYNUB3cs86Syc1vIEFEwub6Gp6pY89RuGN0yh76kn/hRicjgshSOnaFYCSu++m1cKvKbjiclK6GfeSyps4EWEbKEwmREKCfT2g+ftzEULYNeTEkhIsnTtTOP1GWn2z0N5nEouKaPXtNy7lNr39diwdO5LQpAmW1q0dB+qujaV1a6o+/tirCSt79AlkHXec7jF3MkeNpHLO22Ta9yN4v37m9HRSu3cnY/Bge5q/vi+EoOjGaSQWGps5t/5uKa2XLLZ/T+nYgewgFUBzejoJeopemEw2wmx2cwrQyq1v91d3IikU9K6k7qgghDgdqAbu0TsupZwlpayWUlYXuE3bYpHSp//jO4NOJzPn5dq1qeSWLTWtMtgx1Ifd2dk2LhISArZrWtq0AXxP7wEKLprqYce3dOyEpaMm6NKdBopwY05PI/e003ybF+a8TYtP5wf2wNuuqxCYs7PJP++8wG28TmWAw8RgM4uVPvQgRTffTFJZGUIIck891UPwmjNdTTBpRx5J5WuvepqynNZkkkpLPExYwSCEcBU8QdLUx4wnEEwpKRHf7BiKSEhq1oyE3FwsbdtiadvW9V7aCo6xdzhEUiisB8qcvpcCf7tnEkIMAW4AjpNSHoxUY7yt9rsPbkUzZhguU0p9m3xK27aUPeUwq2QMG2b/bEpPR+h0syaXX07e+ecB2vqCMJspvOZqw23xhvvAXO6mPbsPON4eVtt1Mefnk3PqeI9y/FE4/UYKLryApNIS2q5cQdmjjwR0vhFaLV5Mq8WL/WdE056dvUaMYNvw5W5P9sbGJprAFaVO2qGbUCi57z7Kn3navliZkJ9PzinGZm8WP+a9WHxhTOYxxwKaRm+EvMmTfC4oR4UQrqM5I4PEpk297gtxx1/o7Dlz5tCpUye6dOlCdXU1CxYsCLpt3ohk7KNFQEshRCXwFzAOcDEOCyGOAJ4AhkspN0ewLWQMHKibnlxVxb5t2+zfM0eNZKNRwWD1rpGn9+1j/5wzfhwlDz7A9meeJfOYUWzavs4178CBdW6CdZ2vbuBI9WGaaPHZZ6zp39+z3sGD7TOMlCOO8Djuvhjm3OGTKirIGKRdp8IbbgC0TUq1/+wmZ9wpddkFRdOne22XN/TWBUzp6Vj37Am4LG+Y04PQGJ1+f/rAgeyZP99r1oKLLyL7pBNJKtd3Ka76cB4iMZE1AwfVlV3XrksdLroiKQl5+LD9HpvT00k78siAm9162Q8+F0mjRW5KLrDRcP6MQQO9eoLp0eTSS4NoVcNA4GkMcA6dnWKxeITOHjx4MMcddxxCCJYtW8bYsWNZuXJlWNsVMaEgpawRQkwF5gFm4Gkp5U9CiFuAxVLKd9DMRenAa3VazToppTHjZoQwZ2RQdNN0Nt58i9+8ztP1xGblHP5jnfe8QpB3zgTty3Yv+QzsvBRJSVS+9aaLzbXFZ5+ypv8AQHPt3F+nLad01YRC5nHHYkq2YGnfzqtrZUJhIc1m/x+m1FSXh7bJlca0tFaLvuXQ2rWIJLfFtcRECi7Qj/le+fbbHFwV3g4dKM7adOkjD7Ph+uvZt2gxh//2mNQizGavAgGwH0tu2YKDq9c4DjgN3s1ems3ujz/WXYQMBFMU/NWNkJqYRngDTMcI718LGzVHAJO0krSvLjDeIofikbR3r0eaT4o6woiZTgmuSqANX6Gz052Uur1790ZkNhjRKKlSyrnAXLe06U6fh3icFEGSWlRxaM2vAFS+9SbCYmHjTTM88uWMH88/H8xj3zffeBwDaLngCwDSevXi4KpVAJhSnTqGTdPv0YN9OjsUPajLn9azJ9mnnEL+5Enef0Ozco+FqcTCQkoffYQDK1cihCC1e3dK//046XU+4SV3Ozx9bVEWhc1Vsa5PpXQ9Qn+RzSDmjAxSOnb0SG/rYxdqUmmJh5dIfZA3eRIZg4cgzGaa3qV54axo0xaAXB3vJn80e/FFDm/cxMbzNfdN4STsLa1ahcW2r2gE6I3nQjBm+HDufOIJRvTvz/JffuHM0aPtQgHgrbfe4rrrrmPz5s3897//DXuz4ip0dtV777Fu4nnsX7oUS9u2HseTW7awfy6b9QTW3btZ3dc1dC0JCY7NKk5SOm/C2Zizs9m3aLF9d21yixbs+/Zbj9227msKOWdoMU5EQgLFN88I6rdlDB7s4uWRMWCAbj4hBIXTbyStVy/7dyDcjkENCl8mioxjfezyraPZS7NdBKo5KwtzVpbjkor6N/MoAsRJo5eHD3No1SokkNrBsRZy6McfAUjpYGx9xJ2EoiIO//W3wyMNQAh76Ow3Fy1iuI6JePTo0YwePZrPP/+cG2+8kY8//jio+r22K6ylNQC8LZCWP/O0y2YbU3Ky3yl+avdqtj/7LOXPP0da3SJker9+9uNNrr2GjGHDfHprBGJfBSDYoGBOuNj3Y3AxMpYwmfwvEKbqrN24FhLda1w5Zw77De7UVtQf5vR0zK31Z42jBgzg6htu4OM5c9i8erVunn79+vHrr7+ydetW8vPzdfMEQ9wJBXcSy8vg228xZWYGvHCXMXgwLb/+ymUbvTOmpCTSeuqHKDZKxeuvIxLM7P7wQ7Y+9jipXcMcPVTo2zUVdYRDy4/ygrCldSssXgabSFH+3HN219Dktm05uCJAZSeWqQe96czRo8lv1YpOHTvysZNQWLNmDVVVVQghWLp0KYcOHSIvDBs0nYl7oVA0bRoZgweTEuTmNm8CwReBDL8pHbR2Hd6oeXj4C0sQOEoo+CQcWn4cmI8ScnPsfbPipdn18sayyBF9qVBaVMQlQ4YgpXQxTb7xxhs8//zzJCYmkpKSwiuvvBL2xea4Fwomi8Wru2oskTFgABWvvYalQ3h3ZttRQkEXk0H/ct+FNF4TXUKTAg79+uv/t3fvQVKVZx7Hvz+GGQaGiyMIKCMBFZfLMg7CskAsQlYrinFlNRghCWqI5R/W7mIMBLSyZCuptdZKFpPdyl5c0BJXwZUQ40ZXodBZWS8oqNzEgBfUSYyQcTCy4gV99o/zTtM0PdMz06f7dM88n6qpOX1u/czMmX7OeU/38xz3Trxe1dVtVmJ17Tt8+DBHwr0KiO75nT9nDjPCm0uWLl3K0qVLCxpDj08KiehiZu/oB346o/+sLzDgwgsZ9t0luVfuieI4C+vGVwojVqzg8GOPt1tDq+z1sPtunhSSUEJn5b2qq6n76U+SDqN05fOC0LppRfdNCr1ra7tca8iVJk8KndT/vPOSDsEVkTrw7qPcO+lZZ5rdTeuYfSH/jH3GjCmZk0VPCp1w5oZH4+l74C8SZUNl+O4jVyjH/99WjRwZ298230+4x8mTQie0V+Kgc0rjjMDlpjiGfjwpdA8Z/7aZ1Wq7Cz9anWtPHi/oxz7RHEskLik97Mrek4Jz7cln+Kj1tcQvDLuFigH9c6+UJ0knlM4+5ZRTUqWzARobG2loaGDChAl8IUsZjHz58JFz7cjrJDEkA/NTr7ImiT5nn416F/7lsqam5ljpbGDjxo2MSOvpfOjQIa6//noeeeQRRo4cyYED8Xcc8KSQgAFVA3g/6SBcx+QxfNTLx4/aNXr9zzmye3fSYWR167O38vK78ZZ1H3vyWJZOzf3Bs9bS2V+bMoU1a9Ywf/58Gh+Oik3fe++9XH755YwM9zeHxvHGlwx+DpOAvpWFbR/oYhTHTeIeNibdUdXjx1N7xRVJh1Fy5l15Jes2bODDDz9kx44d/Glaoc69e/fS0tLCrFmzmDx5MqtXr479+f1Kwbn25PGCrtYrhW5c5qK76sgZfaFMvuAC3ly2jDVr1nDxxRcft+zo0aNs27aNTZs2ceTIEaZPn860adM4O8YeHZ4UnGtHKfY5dt3fpZdeyuLFi2lsbKS5uZle/fpROWIEdXV1DBkyhJqaGmpqapg5cybbt2+PNSn48FGCKk87LfdKrux5YnGdtXDhQpYvX87E0M1QFRX0rq1lzpw5bN68maNHj/LBBx+wZcsWxmVpGJYPv1JIkPl7FUtfDC/oViLlC1z5qKurY9GiRSfMHzduHBdddBH19fX06tWLa6+9lj/uYue3tnhSSICfOJYR/2O5Ijp8+PAJ82bNmsWstPa6S5YsYcmSwlU19uGjJITa85Wn+vBRyfOk4HoYv1JIQOXQoYy4bQX9pk1LOhTnnDuOJ4WEDJw9O+kQXEf4lYLrYTwplLHh319On7POSjqMbi2fdw5V9aoEPqZPRemURXYuF08KZax2/vykQ3DtGFozjE94jWH94i9F4Fyh+I1m5wrEP5/gypEnBeecKxGSWLBgQepxa+nsSy65BICWlhYuu+wy6uvrmTp1Krt27Yo9Bk8KzjlXImpqati1axdHjhwBTiydfcstt9DQ0MCOHTtYvXp11g+45cvvKTjnXIbf3XILH+2Jt3R2n3FjGX7zzTnXmz17Ng899BBz585Nlc7evHkzAC+99BI3hSY8Y8eOZf/+/bzzzjsMGzYstjj9SsG5LHoPH550CK6HmjdvHmvXrs1aOvucc85h/fr1ADz77LO88cYbNDU1xfr8fqWQQ+2CBbTcfXfSYbgiG3XffXz08p689qGqqjDhN5zLTUfO6Aulvr6e/fv3Zy2dvWzZMhYtWkRDQwMTJ05k0qRJ9I65I5wnhRyG3XyTJ4UeqHLYUCqH5fdW0rp//CcO3X8/VWeeGVNUrqfILJ3dauDAgdx5551AVGhx9OjRjB49Otbn9qSQg7+t0HVVVd0Ihn77hqTDcGVo4cKFDBo0iIkTJ9LY2Jiaf+jQIfr160dVVRUrV65k5syZDBw4MNbn9qTQUaGInXPOFVpbpbP37NnDVVddRUVFBePHj2fVqlWxP7cnhQ4YumwpNTNmJB2Gc66by1U6e/r06ezbt6+gMRT03UeSLpL0a0mvSFqWZXkfSfeF5VskjSpkPF01+JprqI6x3Z1zzpWqgiUFSRXAz4DZwHhgvqTxGat9C2gxs7OA24BbCxWPc8653Ap5pTAVeMXMXjOzj4G1wJyMdeYAd4XpdcD58ju7zrmEdIfWqfn+DIVMCiOAt9IeN4V5Wdcxs6PAe8DgzB1Juk7SVklbDx48WKBwnXM9WXV1Nc3NzWWdGMyM5uZmqquru7yPQt5oznbGn/nb7sg6mNntwO0AU6ZMKd+/mHOuZNXV1dHU1ES5n3hWV1dTV1fX5e0LmRSagNPTHtcBv21jnSZJvYFBwLsFjMk557KqrKyM/YNg5aiQw0fPAWMkjZZUBcwDHsxY50Hg6jA9F3jMyvnazTnnylzBrhTM7KikvwQeBSqAO8xst6QfAFvN7EFgFXC3pFeIrhDmFSoe55xzuRX0w2tm9jDwcMa85WnTHwJXFDIG55xzHadyG62RdBB4I+k4gCHA75MOIguPq3M8rs4r1dg8rvZ9zsxOybVS2SWFUiFpq5lNSTqOTB5X53hcnVeqsXlc8fAmO84551I8KTjnnEvxpNB1tycdQBs8rs7xuDqvVGPzuGLg9xScc86l+JWCc865FE8KzjnnUjwpZJB0h6QDkna1sXyspKclfSRpccayb0vaLWmXpDWSul6qsPNxfV3SjvD1lKRz0pa12+woibgknS7pcUl7wu/sxN6DCcSVtrxC0guSflUqcUk6SdI6SS+H39v0EootyWN/TojpxVBN+by0ZVdL2he+rs62fbHjktQQXkN2h+VXxhlX3szMv9K+gJnAucCuNpYPBf4E+Dtgcdr8EcDrQN/w+D+Ba4oY1wygNkzPBraE6QrgVeAMoArYDowvgbhOBc4N0wOAvaUQV9ryG4F7gV8V+fhqMy6i3iPXhukq4KRSiK0Ejv3+HLs/Wg+8HKZPBl4L32vDdG0JxHU2MCZMnwa8HfffMp8vv1LIYGZP0E6lVjM7YGbPAZ9kWdwb6BsqvvbjxKqwhYzrKTNrCQ+fIapKCx1rdlT0uMzsbTN7Pky/D+zhxH4bRY8LQFId8GVgZVzx5BuXpIFEL0Krwnofm9mhUogtSPLYP2zhFRao4Vj5/QuBjWb2boh7I3BR0nGZ2V4z2xemfwscAHJ+0rhYPCnExMx+A/wYeJMo879nZhsSCudbwH+H6Y40OyqW9LhSQm/uScCWIsfTKjOunwDfBT5LJpyU9LjOAA4Cd4ZhrZWSapIL7VhspXDsS7pM0svAQ8DCMDvxY7+NuNKXTyW66nu1mHG1x5NCTCTVEp2Bjya6JKyR9I0E4vgi0T/s0tZZWVYr+vuQs8TVOr8/8HPgBjP7Q9JxSboEOGBm24odS3txEZ2Jnwv8i5lNAv4PiPX+UFdjK4Vj38x+YWZjgb8AftgaarZVixdVm3EBIOlU4G7gm2aW9AlIiieF+FwAvG5mB83sE2A90Rhs0UiqJxrymGNmzWF2R5odJREXkiqJEsI9Zra+mDG1E9fngUsl7ScaavszSf9RAnE1AU1m1no1tY4oSRRVG7Elfuy3CkM6Z0oaQgkc+23E1Toc+BDwPTN7JomY2uJJIT5vAtMk9ZMk4HyicfKikDSS6J9xgZntTVvUkWZHRY8r/I5WAXvMbEWx4skVl5ndZGZ1ZjaK6Hf1mJkV7ay3nbh+B7wl6Y/CrPOBl4oVV3uxkfyxf1Z4XiSdSzQc00zUy+VLkmrD1cyXwrxE4wr/h78AVpvZ/cWKp6MK2k+hHElaA8wChkhqAr4PVAKY2b9KGg5sBQYCn0m6gehdM1skrQOeB44CLxDjx9tzxQUsBwYD/xyOw6NmNsXaaHaUdFxEZ+QLgJ2SXgy7u9miHhxJxlVQecb1V8A94UXlNeCbpRBbCRz7XwGukvQJcAS4MtzgfVfSD4lOjAB+YGaxtfvtalySvkr0poHBkq4Ju7vGzF6kBHiZC+eccyk+fOSccy7Fk4JzzrkUTwrOOedSPCk455xL8aTgnHMuxZOCc865FE8KzjnnUjwpuJIm6bbwAcHWx49KWpn2+B8k3Rjzcx6OeX8nSbo+7fGotmrwZ9m2r6T/kVQRQxxVkp4IlUydy8qTgit1TxHq6EjqBQwBJqQtnwE8mUBcnXEScH3OtbJbCKw3s0/zDSKUTt8ElFZTF1dSPCm4Uvckx4qrTQB2Ae+HejZ9gHHAC5IekLQtdLO6rnVjSbdmnKX/raTvSPqGpGcVdcX6t2xn4m2tE87090j69/B8GyT1Dcv+RlFntI2KOpAtBv6eqBjai5J+FHZfkW37LL4O/DItpsbW+keSBivqdDYqPOfK8PgeSRdIelJRx7Gpaft7IOzTuaw8KbiSFpqQHA3F2GYATxP1XZgOTAF2hDPghWY2Ocz7a0mDwy7WcvyZ8VeJalddCXzezBqAT8l4oZQ0Lsc6Y4CfmdkE4BDwFUlTiOrdTAIuD7FAVOL6VTNrMLMlbW2f+bOHGkdnmNn+tNlnAfvCdD2wM23+T8O8scDXgPOAxcDNadvvIuoc6FxWPrboykHr1cIMYAVRo5QZwHtEw0sQJYLLwvTpRC+6zWb2gqShkk4j6m7VAkwEJgPPhcJufYm6X6U7P8c6r6cVMNsGjCIa2vqlmR0BkPRf7fxM2bbPNIQoYRD29zngN2m19+uBHWn72xnW2w1sCsXXdqbv28w+lfSxpAGh451zx/Gk4MpB632FiURnum8B3wH+ANwhaRZRTf/pZvaBpEYgvXH8OmAuMJzoykHAXWZ2UzvPmWudj9KmPyVKGtmaurQl2/aZjnD8z9HAsSQAUdK6L8v+Pkt7/Bkn/p/3AT7sRKyuB/HhI1cOngQuAd41s09D+eOTiIaQngYGAS0hIYwFpmVsv5aoN8JcogSxCZgraSiApJPDWXi6jqyT6X+BP5dUraij3JfD/PeBAZ39oUNf4QpJrYnhHEKSkDSGqNvZzjY2zyoMq7U2w3HuBJ4UXDnYSTSU8kzGvPfM7PfAI0BvSTuIWh4e18kq9I8YQDT08raZvQR8D9gQttkInJqxTc51MpnZc0QNjLYTNaPZGmJsBp4MN4F/1N4+sthAdG8AoiuFXpK2E/U22ANc3cn9fRGIpWeF6568n4JzMZLU38wOS+oHPAFcZ2bP57G/ScCNZrZA0ivApHzuBUhaD9xkZr/u6j5c9+b3FJyL1+2SxhMN89yVT0IACDfKH5c0CPgsz4RQBTzgCcG1x68UnHPOpfg9BeeccymeFJxzzqV4UnDOOZfiScE551yKJwXnnHMpnhScc86leFJwzjmX8v9a61IdOFBhfgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# The 4 spectra examples\n",
    "plt.plot(wav1, flux1, label=\"M0\")\n",
    "plt.plot(wav2, flux2, label=\"M3\")\n",
    "plt.plot(wav3, flux3, label=\"M6\")\n",
    "plt.plot(wav4, flux4, label=\"M9\")\n",
    "plt.legend()\n",
    "plt.ylabel(\"Flux\")\n",
    "plt.xlabel(\"Wavelength ($\\mu$m)\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Only dealing with one band so band=band\n",
    "norm1 = snr_constant_band(wav1, flux1, band=band, snr=snr)\n",
    "norm2 = snr_constant_band(wav2, flux2, band=band, snr=snr)\n",
    "norm3 = snr_constant_band(wav3, flux3, band=band, snr=snr)\n",
    "norm4 = snr_constant_band(wav4, flux4, band=band, snr=snr)\n",
    "\n",
    "flux1 /= norm1\n",
    "flux2 /= norm2\n",
    "flux3 /= norm3\n",
    "flux4 /= norm4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This should now be 100! if band == J\n"
     ]
    }
   ],
   "source": [
    "print(\"This should now be 100! if band == J\")\n",
    "norm4 = snr_constant_band(wav4, flux4, band=band, snr=snr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "x1, q1 = incremental_quality(wav1, flux1, percent=1)\n",
    "x2, q2 = incremental_quality(wav2, flux2, percent=1)\n",
    "x3, q3 = incremental_quality(wav3, flux3, percent=1)\n",
    "x4, q4 = incremental_quality(wav4, flux4, percent=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEOCAYAAABB+oq7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XuU1XW9//HnWy4OoKBcVGREINQRBdEmRXMZZr9SJLHy1kUMbNGpTtopTXClnKw8urrX8VikknhSUDMxxJJDUWaKAiooF0FCGUS5g8pFRt+/P76fgc2w9549M9/vd1/m9Vhr1uz9+d7ee2Dv9/58Pzdzd0REROJwQLEDEBGRyqGkIiIisVFSERGR2CipiIhIbJRUREQkNkoqIiISGyUVERGJjZKKiIjERklFRERi077YAaStZ8+e3q9fv2KHISJSNubPn7/B3XsVsm+bSyr9+vVj3rx5xQ5DRKRsmNmrhe6r218iIhIbJRUREYmNkoqIiMSmzbWpiDTX7t27qaurY+fOncUOpVWqqqqorq6mQ4cOxQ5FKpiSikgT6urqOPjgg+nXrx9mVuxwWsTd2bhxI3V1dfTv37/Y4UgF0+0vkSbs3LmTHj16lG1CATAzevToUfa1LSl9SioiBSjnhNKgEl6DlD7d/hIRKTObp93PthkzmnXMgcfXcMT11ycU0V6qqYiUATPj8ssv3/O8vr6eXr16MXLkSCBqM7nqqqsYOHAgQ4YMYcGCBcUKVVKwbcYMdi5dWuwwslJNRaQMdOnShRdffJEdO3bQqVMnZs2aRZ8+ffZsf+yxx1i+fDnLly9n7ty5fOUrX2Hu3LlFjFiSVlVTw9H3TCl2GPtRTUWkTJx33nk8+uijANx333189rOf3bNt+vTpjB49GjNj2LBhbNmyhbVr1xYrVGnDVFMRaYbv/vElFr++LdZzDjqyKxM/eUKT+1122WXcdNNNjBw5koULFzJ27FieeOIJANasWcNRRx21Z9/q6mrWrFlD7969Y41VpCmqqYiUiSFDhrBq1Sruu+8+RowYsc82d99vf/X2kmJQTUWkGQqpUSTpggsu4JprrmHOnDls3LhxT3l1dTWrV6/e87yuro4jjzyyGCFKG6eaikgZGTt2LDfeeCODBw/ep/yCCy5gypQpuDtPP/003bp1060vKQrVVETKSHV1NVdfffV+5SNGjGDmzJkMHDiQzp07M3ny5CJEJ6KkIlIW3n777f3Khg8fzvDhw4Go/eS2225LOSqR/en2l4iIxEZJRUREYqOkIiIisVFSERGR2CipiIhIbJRUREQkNkoqImWgqanvp0+fzpAhQxg6dCi1tbX84x//KFao0sZpnIpIGWhq6vtzzjmHCy64ADNj4cKFXHLJJSwt0fU2pLKppiJSJvJNfX/QQQftmUDynXfe0WSSUjSqqYg0x2Pj4Y1F8Z7ziMFw3i1N7pZv6nuAP/zhD0yYMIF169btST4iaVNNRaRM5Jv6HuBTn/oUS5cu5eGHH+aGG24oQoQiqqmINE8BNYok5Zr6PtNZZ53FK6+8woYNG+jZs2fKEUpbp6QiUkbGjh1Lt27dGDx4MHPmzNlTvmLFCj7wgQ9gZixYsIB3332XHj16FC9QabOUVETKSK6p73//+98zZcoUOnToQKdOnZg2bZoa66UolFREykBTU99fd911XHfddSlHJbI/NdSLiEhsVFMp0Bs338yuJc0bTNZ15EgOvfSShCISESk9qqkkZOfSpWybMaPYYYiIpCqxmoqZ3QWMBNa5+4mhrDswDegHrAIucffNFrUo/hwYAWwHvujuC8IxVwDfCaf9vrvfHco/CPwW6ATMBK52d0/q9Rxx/fXN2v/Vy0cnFImISOlKsqbyW+DcRmXjgdnufgwwOzwHOA84JvyMA26HPUloInAacCow0cwODcfcHvZtOK7xtUREJGWJJRV3/zuwqVHxKODu8Phu4MKM8ikeeRo4xMx6A58AZrn7JnffDMwCzg3burr7U6F2MiXjXCIiUiRpt6kc7u5rAcLvw0J5H2B1xn51oSxfeV2W8qzMbJyZzTOzeevXr2/1ixBJW1NT3wPMmTOHoUOHcsIJJ/CRj3ykGGGKlEzvr2yjtLwF5Vm5+yRgEkBtbW1i7S4iSWlq6vstW7bw1a9+lT/96U/07duXdevWFTFaacvSrqm8GW5dEX43/M+vA47K2K8aeL2J8uos5SIVK9/U9/feey+f/vSn6du3LwCHHXZY1nOIJC3tmsojwBXALeH39IzyfzezqUSN8lvdfa2Z/Rm4OaNx/uPABHffZGZvmdkwYC4wGvhlmi9E2qZbn7mVpZviXfyqpnsN153a9Gj4fFPfv/zyy+zevZvhw4fz1ltvcfXVVzN6tHogSvqS7FJ8HzAc6GlmdUS9uG4B7jezK4HXgIvD7jOJuhOvIOpSPAYgJI/vAc+G/W5y94bG/6+wt0vxY+FHpGLlm/q+vr6e+fPnM3v2bHbs2MHpp5/OsGHDOPbYY4sUrbRViSUVd/9sjk3nZNnXga/lOM9dwF1ZyucBJ7YmRpHmKqRGkaRcU99XV1fTs2dPunTpQpcuXTjrrLN44YUXlFQkdRpRL1JGxo4dy4033sjgwYP3KR81ahRPPPEE9fX1bN++nblz53L88ccXKUppy0ql95eIFCDX1PfHH3885557LkOGDOGAAw7gS1/6EieeqIq8pE9JRaQMNDX1PcC1117Ltddem2JUIvvT7S8REYmNaioJ2rl0aaITS2pqfREpNUoqCemaMX1GEnYujcZKKKmISClRUknIoZdekugHvqbWF5FSpDYVERGJjZKKiIjERre/RMqAmfGFL3yBe+65B4imZenduzennXYaM2bMYPPmzYwdO5ZXXnmFqqoq7rrrLo1TCTZPuz/Vpb3begca1VREykDm1PfAflPf33zzzQwdOpSFCxcyZcqUrAMk26ptM2bs6diStJ1Ll6aawEqRaioiZaJh6vuLLrpoz9T3DbMUL168mAkTJgBQU1PDqlWrePPNNzn88MOLGXLJqKqp4eh7piR+HXWgUVIRaZY3br6ZXUvi/dZ74PE1HHH99U3ul2/q+5NOOomHHnqIM888k2eeeYZXX32Vuro6JRVJnW5/iZSJfFPfjx8/ns2bNzN06FB++ctfcvLJJ9O+vb4zSvr0v06kGQqpUSQp19T3Xbt2ZfLkyQC4O/3796d///7FClPaMCUVkTIyduxYunXrxuDBg5kzZ86e8i1bttC5c2c6duzIHXfcwVlnnUXXrl2LF6i0WUoqImUk19T3S5YsYfTo0bRr145BgwZx5513FiE6ESUVkbLQ1NT3p59+OsuXL085KpH9qaFeRERio6QiIiKxUVIRKYC7FzuEVquE1yClT20qIk2oqqpi48aN9OjRAzMrdjgt4u5s3LiRqqqqYodS8ZJenK/hGlU1NYleo6WUVESaUF1dTV1dHevXry92KK1SVVVFdXV1scOoaEkvztegqqYmtWs1l5KKSBM6dOiggYRSkKQX5ysHTbapmNmPzOyENIIREZHyVkhD/VJgkpnNNbN/M7NuSQclIiLlqcmk4u53uPuHgdFAP2Chmd1rZmcnHZyIiJSXgtpUzKwdUBN+NgAvAN80sy+7+2UJxiciFSqtFRlLuadUJSqkTeUnRLfARgA3u/sH3f1Wd/8kcHLSAYpIZUprRcZS7ilViQqpqbwIfMfdt2fZdmrM8YhIG5LWioySnkIa6j/fOKGY2WwAd9+aSFQiIlKWctZUzKwK6Az0NLNDgYahxF2BI1OITUREyky+219fBr5BlEAWZJRvA25LMigRESlPOW9/ufvP3b0/cI2798/4Ocnd/7s1FzWz/zCzl8zsRTO7z8yqzKx/GAuz3MymmVnHsO+B4fmKsL1fxnkmhPJlZvaJ1sQkIiKtl+/210fd/S/AGjP7dOPt7v5QSy5oZn2Aq4BB7r7DzO4HLiPqXfZTd59qZr8CrgRuD783u/tAM7sMuBW41MwGheNOIKpN/Z+ZHevu77UkLpHWSKt7bNeRI9v8NCACD7z8ADNXzmzWMTXda7ju1OsSimivfA31Hwm/P5nlp7X989oDncysPVG7zVrgo8CDYfvdwIXh8ajwnLD9HIumih0FTHX3Xe7+L2AF6o0mRZJG99idS5emkrik9M1cOZNlm5YVO4ysctZU3H1i+D0mzgu6+xoz+xHwGrADeByYD2xx9/qwWx3QJzzuA6wOx9ab2VagRyh/OuPUmce0CWlMsa1vxoVLunvsq5ePTuXfHPTvXg6O634ck8+dXOww9pPv9tc38x3o7j9pyQVDT7JRQH9gC/AAcF62SzQckmNbrvJs1xwHjAPo27dvMyMuTWkM5mr45q0Pl9KQ1gA+/btLa+Tr/XVwQtf8GPAvd18PYGYPAWcAh5hZ+1BbqQZeD/vXAUcBdeF2WTdgU0Z5g8xj9uHuk4BJALW1tRWx/F0aU2yn8Y1YCpfWtOr6d5fWyHf767sJXfM1YJiZdSa6/XUOMA/4K3ARMBW4Apge9n8kPH8qbP+Lu7uZPQLcG6aRORI4BngmoZhFRKQATU7TEgZBXknUy2rPWqTuPrYlF3T3uWb2INHYl3rgOaJaxKPAVDP7fii7MxxyJ3CPma0gqqFcFs7zUug5tjic52vq+SUiUlyFzP11D9GEkp8AbgI+DyxpzUVDJ4CJjYpXkqX3lrvvBC7OcZ4fAD9oTSwiIhKfQub+GujuNwDvuPvdwPnA4GTDEhGRclRIUtkdfm8xsxOJGsr7JRaRiIiUrUJuf00K3YBvIGo0Pwi4MdGoRESkLDWZVNz9jvDwb8CAZMMpXbc+cytLNzVvxPSIASO4+NiszUEiIhWpkN5fWWsl7n5T/OFUjoYpFJRURKQtKeT21zsZj6uI5v1qVe+vctTcidjG/CnW2W1ERMpCIbe/fpz5PMzb9UhiEYmISNkqpPdXY51pw20rIiKSWyFtKovYO1FjO6AX0SBIERGRfRTSppI5NWo98GbGFPUiIiJ7FJJUjgEGhcfz3H1NgvGIiEgZy7eeylFEMwW/RbSIlgGfMbMdROuhXJ4xhkVERCRvTeU24Bfu/tvMQjMbTTQNPYCSioiI7JEvqdQ0TigA7j7FzG4GTkksKikZaS1fmyQtjSuSnnxJJWt3YzM7ANjh7uuSCUlKRVrL1yZJS+OKpCtfUvmjmf0G+Ia7vwNgZl2AnwIz0wiu3C3btCzRkfVJzy2W1vK1SSr3WpZIuck3+PHbwFbgVTObb2bzgFXANuDaFGIrayMGjOC47scldv5lm5Yxc6Vyu4iUlnxr1O8GrjGzG4CBRL2/Vrj79rSCK2cXH3txorUIzS0mIqWokLm/dgCLUohFRETKXEvm/hIREclKSUVERGKTb0R93nEo7r4g/nBERKSc5WtT+XGebQ58NOZYRESkzOXr/XV2moGIJCWNWQF2Ll1KVU1NotcQKQeFzFKMmZ1INFNxVUOZu09JKiiRuKQ1K0BVTU1FzEAg0lqFLNI1ERhOlFRmAucB/wCUVKTkVcKsAMWg2p20VCE1lYuAk4Dn3H2MmR2OZicWqViVVrt74OUHUp19Iunpk0pdIUllh7u/b2b1ZtYVWIfWqBepWJVWu5u5cibLNi1LdNqkBss2LQNQUmnCPDM7BPgN0WJdbwPPJBqVSJlJ69twJX0LTutv1pBQJp87OfFrafqkJgY/mpkB/+XuW9z9V8D/A65wd/3lRDI0fBtOUqVNIprG3wzguO7HMWLAiMSvI5G8NRV3dzN7GPhgeL4qjaBEylHS34Yr8VtwWjUISU8ht7+eNrMPufuziUcjzZL0ei1QWbdbRNKQxvsyrTailigkqZwNfNnMXgXeIZoC3919SKKRSV5pVOfV6CjSPGndZivlW3qFJJXz4r5oaPi/AziRaMqXscAyYBrQj2gxsEvcfXNo1/k5MALYDnyxYd4xM7sC+E447ffd/e64Yy1VSa/XApV5u0UkSWm8L0tdIbMUf9/dX838Ab7fyuv+HPiTu9cQjYFZAowHZrv7McDs8ByipHZM+BkH3A5gZt2BicBpwKnARDM7tJVxiYhIKxSSVE7IfGJm7QgN9y0RxrqcBdwJ4O7vuvsWYBTQUNO4G7gwPB4FTPHI08AhZtYb+AQwy903uftmYBZwbkvjEhGR1suZVMxsgpm9BQwxs23h5y2iwY/TW3HNAcB6YLKZPWdmd5hZF+Bwd18LEH4fFvbvA6zOOL4ulOUqFxGRIsmZVNz9v9z9YOCH7t41/Bzs7j3cfUIrrtkeOAW43d1PJmr8H59nf8sWXp7y/U9gNs7M5pnZvPXr1zc3XhERKVAht7+eMbNuDU/M7BAzuzDfAU2oA+rcfW54/iBRknkz3NYi/F6Xsf9RGcdXA6/nKd+Pu09y91p3r+3Vq1crQhcRkXwKSSoT3X1rw5PQ/jGxpRd09zeA1WbW0Mn6HGAx8AhwRSi7gr232B4BRltkGLA13B77M/BxMzs0NNB/PJSJiEiRFNKlOFviKWgdljy+DvzOzDoCK4Ex4Tr3m9mVwGtAQ7+8mUTdiVcQdSkeA+Dum8zse0DDoMyb3H1TK+MSKWlpDKwDDXqVlit0QsmfALcRtVl8nWhiyRZz9+eB2iybzsmyrwNfy3Geu4C7WhOLSLlIa7CbBr1KaxSSVL4O3EA0MBHgcfYOOBSRlKQ1sE6DXqU1mkwq7v4OMN7MDnL3t1OISUREylSTDfVmdoaZLSZqTMfMTjKz/0k8MhERKTuF9P76KdHo9Y0A7v4C0Yh4ERGRfRSSVHD31Y2K3ksgFhERKXOFNNSvNrMzAA9dgK8imgBSRERkH4XUVP6NqEtvH2ANMJQcXXxFRKRtK6T31wbg8ynEIiIiZa6Q3l8DzOyPZrbezNaZ2XQzG5BGcCIiUl4KaVO5l2g0/afC88uA+4gWx5IKl9a0IEnSlCMi6SmkTcXc/R53rw8//0uOKealsowYMILjuh/X9I4lbNmmZcxcObPYYYi0GYXUVP5qZuOBqUTJ5FLg0bCcL5rEsXJVwnrb5V7LEik3hSSVS8PvLzcqH0uUZNS+IiIiQGG9v/qnEYiItC3rtu1iwzu7uPTXTyV+rVFD+/C50/omfh3Jv0b9h8zsiIzno0PPr1803PoSEWmpDe/sYvuu+sSvs3jtNqY/vybx60gkX03l18DHAMzsLOAWomnwhwKTgIsSj05EKlrnA9szbczpiV4jjZqQ7JUvqbTLaIS/FJjk7r8Hfm9mzycfmkg80ugWvWzTsrLvKZe2Q9/bSNf3t8Dk8xO9zo0bt/Jkp7OBZJOXRPImFTNr7+71RCsyjivwOJGSkdZqicd1Py61a1WKru9vocp3Jn6dfrtXJn4N2StfcrgP+JuZbQB2AE8AmNlAYGsKsYm0WiV0iy6GNGp3r3Z4j6N3V8GYRxO9zqqbz0z0/LKvnEnF3X9gZrOB3sDjYa14iBr3v55GcCKSvrRqXEfvbseHt3dM5VqSnry3sdz96SxlLycXjogUW1q1u5dUg6hIBS3SJSIiUgglFRERiY2SioiIxEZJRUREYqPxJuVq3mRY9GDy1xl8EdRqpl8RKYxqKuVq0YPwxqJkr/HGonQSl4hUDNVUytkRg5MdOJbw9BkiUnlUUxERkdgoqYiISGyUVEREJDZKKiIiEhslFRERiY16f5WpN9/ayYa3d3FTgqva3bhxKz0POpDDE7uCiFSaotVUzKydmT1nZjPC8/5mNtfMlpvZNDPrGMoPDM9XhO39Ms4xIZQvM7NPFOeVFMeGt3ex/d33Er3G9nffY8PbuxK9hohUlmLWVK4GlgBdw/NbgZ+6+1Qz+xVwJXB7+L3Z3Qea2WVhv0vNbBBwGXACcCTwf2Z2rLsn+0lbQjp3bMe0Lye3ROpLN7dL7NwVRzMciABFSipmVg2cD/wA+KaZGfBR4HNhl7uB/yRKKqPCY4AHgf8O+48Cprr7LuBfZrYCOBVI5n7QY+ObP4JdHwBtR8MMB0cMTu4aDf//9H9KSlixaio/A74NHBye9wC2uHt9eF4H9AmP+wCrAdy93sy2hv37AJmLiGUeU3z6AGh7NMOBSPpJxcxGAuvcfb6ZDW8ozrKrN7Et3zGNrzkOGAfQt2/fZsW7x3m3NG//yedHiSWhD4J+u1eyqsOARM4tItJSxWio/zBwgZmtAqYS3fb6GXCImTUkuWrg9fC4DjgKIGzvBmzKLM9yzD7cfZK717p7ba9eveJ9NbkMvijRWyGrOgzgyU5nJ3Z+EZGWSL2m4u4TgAkAoaZyjbt/3sweAC4iSjRXANPDIY+E50+F7X9xdzezR4B7zewnRA31xwDPpPla8qodk+itr4auxOMSu4KISPOV0jiV64CpZvZ94DngzlB+J3BPaIjfRNTjC3d/yczuBxYD9cDXkuz59d0/vsTi17c165hRQ/vwudNaeLtNRKQMFTWpuPscYE54vJKo91bjfXYCF+c4/gdEPchKzuK1UQJSUhGRtqSUaiolbeInT2jW/pcmONJdRKRUae4vERGJjWoqCVq8dltiNZbFa7cxqHfXpncUEUmRkkpCRg1NdhzmoN5dE7+GiEhzKakk5HOn9VUjvYi0OWpTERGR2KimIhIDrW8jQHqzVbfEEYObP91UCyipiMQgzfVtlFSar9/ulelMyPnqP6LfR5+Z6GUavsQ0x1s7NjIsoXgyKamIxETr25SmhjnymjfSrIWOPjOVJS+u+vVTLN7UvB6gg47sqqQiItJaszuPYHbnEUwbk1zCL4ZBvbsm+iWmpdRQLyIisVFNRUQqXpIDkRtr6xPJKqmIyL5S6sGU1kJzaQ4S1kSySipS6SrsAzIVix6MVi1NcJE52LvQXNIN6GkORNZEskoqUukq7AMyla6xDX+vMY8mehktNFeZlFSk8lXIB2RqXWOPGBx1i5UWSaP9ppQnlFVSESkTldo1tpKk1X5TyhPKKqmIyD7unfsa059fk/h1SvnbdktpIlklFWlC4vfwUxh9LM0z/fk1qXzgl/K3bWk5JRXJKfF7+G8sin4rqZScUh2tLaVPSUVySvwefhoT/IlIqjRNi4iIxEZJRUREYqOkIiIisVGbihTXG4uSbVtJYTR9mtr6wDopfUoqkleSH2LnbD+FUV12JruSYQWNDtfAOikHSiqSU9IfLL/ceiaze2uEeKE0sE7KgZKK5JT0h5hmdBWpPGqoFxGR2CipiIhIbHT7SyqaJkcUSZdqKlLRGiZHTJp6TIlEVFORiqfJEUXSk3pNxcyOMrO/mtkSM3vJzK4O5d3NbJaZLQ+/Dw3lZma/MLMVZrbQzE7JONcVYf/lZnZF2q9FRET2VYyaSj3wLXdfYGYHA/PNbBbwRWC2u99iZuOB8cB1wHnAMeHnNOB24DQz6w5MBGoBD+d5xN03p/6KpMWSHiGutg6RdKVeU3H3te6+IDx+C1gC9AFGAXeH3e4GLgyPRwFTPPI0cIiZ9QY+Acxy900hkcwCzk3xpUgrjRraRwtBiVSYorapmFk/4GRgLnC4u6+FKPGY2WFhtz7A6ozD6kJZrnIpExohLlJ5itb7y8wOAn4PfMPd83XPsSxlnqc827XGmdk8M5u3fv365gcrIiIFKUpSMbMORAnld+7+UCh+M9zWIvxeF8rrgKMyDq8GXs9Tvh93n+Tute5e26tXr/heiIiI7KMYvb8MuBNY4u4/ydj0CNDQg+sKYHpG+ejQC2wYsDXcJvsz8HEzOzT0FPt4KBMRkSIpRpvKh4HLgUVm9nwoux64BbjfzK4EXgMuDttmAiOAFcB2YAyAu28ys+8Bz4b9bnL3Tem8BBERycbcszZDVKza2lqfN29escMQESkbZjbf3WsL2VfTtIiISGyUVEREJDZKKiIiEps216ZiZuuBVxM6fU9gQ0LnjotijIdijEc5xAjlEWeSMR7t7gWNx2hzSSVJZjav0MasYlGM8VCM8SiHGKE84iyVGHX7S0REYqOkIiIisVFSidekYgdQAMUYD8UYj3KIEcojzpKIUW0qIiISG9VUREQkNkoqBTCzu8xsnZm9mGN7jZk9ZWa7zOyaRtv+Iyyb/KKZ3WdmVUWK8fNhOeaFZvZPMzspY9u5ZrYsLNk8Pon4WhNjriWoSynGjO3tzOw5M5tRijGa2SFm9qCZLQ1/z9NLMMZU3jMFxjkqxPh8WD7jzIxtqSxn3tIYzWxo+Fx6KWy/NKkY9+Hu+mniBzgLOAV4Mcf2w4APAT8Arsko7wP8C+gUnt8PfLFIMZ4BHBoenwfMDY/bAa8AA4COwAvAoBKLsTdwSnh8MPByqcWYsf2bwL3AjCL+f8wZI9Gqql8KjzsCh5RSjGm+ZwqM8yD2NhMMAZaGx92BleH3oeHxoSUW47HAMeHxkcDapP69M39UUymAu/8dyDkDsruvc/dngd1ZNrcHOplZe6AzOdZ8SSHGf3q07DLA00TrzwCcCqxw95Xu/i4wlWgJ55KJ0XMvQV0yMQKYWTVwPnBHErG1NkYz60r0AXVn2O9dd99SSjEGqbxnCozzbQ+fykAX9i4EmNpy5i2N0d1fdvfl4fHrRGtUJb6glJJKgtx9DfAjoqn81xKtBfN4caMC4ErgsfC4VJdlzoxxD9t3Cepiaxzjz4BvA+8XJ5ysMmMcAKwHJodbdHeYWZfihbbHnhhL8T1jZp8ys6XAo8DYUFxS75scMWZuP5WoZvpK0rEoqSQoLB42CuhPVP3sYmZfKHJMZxO9ia9rKMqyW1G7BGaJsaG80CWoE9c4RjMbCaxz9/nFjCtTlr9je6LbKLe7+8nAO0BibWiFyPJ3LLn3jLv/wd1rgAuB74Xiknrf5IgR2LOS7j3AGHdP/AuPkkqyPgb8y93Xu/tu4CGie8lFYWZDiG7NjHL3jaG44GWZ05AjxlxLUBdFjhg/DFxgZquIbiF+1Mz+t0gh5vu3rnP3hlreg0RJpihyxFhS75lM4TbUB8ysJyX2vmnQKMaGW56PAt9x96fTiEFJJVmvAcPMrLOZGXAOUXtA6sysL9Eb9HJ3fzlj07PAMWbW38w6ApcRLeFcMjGGv122JahTlytGd5/g7tXu3o/ob/gXdy/KN+w8Mb4BrDaz40LROcDiIoT+CiZrAAAEwElEQVSY7/9jybxnAMxsYIgDMzuF6BbSRkpoOfNcMYb38x+AKe7+QFrxFGM54bJjZvcBw4GeZlYHTAQ6ALj7r8zsCGAe0BV438y+QdQ7aa6ZPQgsAOqB50ho1GtTMQI3Aj2A/wn//+rdvdbd683s34neEO2Au9z9pVKKkRxLULv7zBKKMTWtjPHrwO/CB85KwvLcpRJjmu+ZAuP8DDDazHYDO4BLQ6N4asuZtzRGM7uEqGNGDzP7YjjdF939eRKkEfUiIhIb3f4SEZHYKKmIiEhslFRERCQ2SioiIhIbJRUREYmNkoqIiMRGSUVERGKjpCIVzcx+GgajNjz/s5ndkfH8x2b2zZiv+XbM5zvEzL6a8bxfrrU1shzbycz+ZmbtYoijo5n9PcweLJKVkopUun8S5o4yswOAnsAJGdvPAJ4sQlzNcQjw1Sb3ym4s8JC7v9faIMLSCLOBdBZ7krKkpCKV7kn2Tkh4AvAi8FaYs+lA4HjgOTN72Mzmh1XyxjUcbGa3Nqol/KeZfcvMvmBmz1i02t6vs9UEcu0TahpLzOw34XqPm1mnsO0Gi1ZlnGXRqofXALcQTRL4vJn9MJy+Xbbjs/g8MD0jpjkNc3+ZWQ+LVlfsF655R3j+OzP7mJk9adGqhqdmnO/hcE6RrJRUpKKFxYnqwwSGZwBPEa3FcjpQCywM38DHuvsHQ9lVZtYjnGIq+34zv4RonrdLgQ+7+1DgPRp90JrZ8U3scwxwm7ufAGwBPmNmtUTzOJ0MfDrEAtH09K+4+1B3vzbX8Y1fe5jfa4C7r8ooHggsD4+HAIsyyn8eymqAzwFnAtcA12cc/yLRKqciWeneqLQFDbWVM4CfEC2mdAawlej2GESJ5FPh8VFEH9ob3f05MzvMzI4kWjVvMzAY+CDwbJgMsRPRqnqZzmlin39lTOw3H+hHdGtuurvvADCzP+Z5TdmOb6wnUcIhnO9oYE3GmhpDgIUZ51sU9nsJmB0mJVyUeW53f8/M3jWzg8MqnCL7UFKRtqChXWUw0Tft1cC3gG3AXWY2nGgdj9PdfbuZzQGqMo5/ELgIOIKo5mLA3e4+Ic81m9pnV8bj94iSTraFn3LJdnxjO9j3dQxlbxKBKOlNy3K+9zOev8/+nxMHAjubEau0Ibr9JW3Bk8BIYJO7vxemKD+E6BbYU0A3YHNIKDXAsEbHTyVaI+UiogQzG7jIzA4DMLPuoRaQqZB9GvsH8Ekzq7JolcvzQ/lbwMHNfdFh7fR2ZtaQWE4iJBkzO4ZohcVFOQ7PKtwWbFhAS2Q/SirSFiwiuhX0dKOyre6+AfgT0N7MFhItxbrPCnlhfZmDiW4drXX3xcB3gMfDMbOA3o2OaXKfxtz9WaIF0l4gWsBqXohxI/BkaET/Yb5zZPE4UdsIRDWVA8zsBaL1TJYAVzTzfGcDsa9jI5VD66mIlBAzO8jd3zazzsDfgXHuvqAV5zsZ+Ka7X25mK4CTW9MWYmYPARPcfVlLzyGVTW0qIqVlkpkNIrpNdXdrEgpA6GjwVzPrBrzfyoTSEXhYCUXyUU1FRERiozYVERGJjZKKiIjERklFRERio6QiIiKxUVIREZHYKKmIiEhslFRERCQ2SioiIhKb/w+I3iDFSeHo/QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(x1, q1, label=\"M0\", drawstyle=\"steps-mid\")\n",
    "plt.plot(x2, q2, label=\"M3\", drawstyle=\"steps-mid\")\n",
    "plt.plot(x3, q3, label=\"M6\", drawstyle=\"steps-mid\")\n",
    "plt.plot(x4, q4, label=\"M9\", drawstyle=\"steps-mid\")\n",
    "plt.legend()\n",
    "plt.ylabel(\"Spectral Quality\")\n",
    "plt.xlabel(\"Wavelength ($\\mu$m)\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "xrv1, rv1 = incremental_rv(wav1, flux1, mask=None, percent=1)\n",
    "xrv2, rv2 = incremental_rv(wav2, flux2, mask=None, percent=1)\n",
    "xrv3, rv3 = incremental_rv(wav3, flux3, mask=None, percent=1)\n",
    "xrv4, rv4 = incremental_rv(wav4, flux4, mask=None, percent=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEOCAYAAACEiBAqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XmUVeWZ7/HvT0QKERwApaQk4IhgYdldzl6DmpWrSKBNjFM7RJJrxquJwUbtG00nHVe8ncRk3bbTIUaiRlHjEAwOiSEShxi0MAQQcCamCAqiiAPK4HP/2LuwxBpOVZ29z6lzfp+1anHOPnvv96mi6jznffd+n1cRgZmZVa9tSh2AmZmVlhOBmVmVcyIwM6tyTgRmZlXOicDMrMo5EZiZVTknAjOzKudEYGZW5ZwIzMyq3LalDqAQQ4YMiZEjR5Y6DDOzXmX+/PmvRMTQzvbrFYlg5MiRNDU1lToMM7NeRdJfC9nPQ0NmZlXOicDMrMplnggk9ZH0Z0mz0+c/l/SCpAXpV0PWMZiZWfvyuEZwAbAUGNRq20URcVsObZuZtWvjxo00NzfzzjvvlDqUHqmpqaGuro6+fft26/hME4GkOuBE4DvAhVm2ZWbWVc3NzQwcOJCRI0ciqdThdEtEsGbNGpqbmxk1alS3zpH10NAPgX8B3ttq+3ckLZR0laR+bR0o6TxJTZKaVq9enXGYZlaN3nnnHQYPHtxrkwCAJAYPHtyjXk1miUDSRGBVRMzf6qVLgNHAwcAuwLS2jo+I6RHRGBGNQ4d2ehusmVm39OYk0KKn30OWQ0NHApMkTQBqgEGSfhERZ6avvytpBjA1wxgsAzfNe5FZC1Zk3s7khuGcceiIzNsxq3aZ9Qgi4pKIqIuIkcBpwO8j4kxJtQBKUtg/AYuzisGyMWvBCpasXJdpG0tWrssl2ZiVmiTOOuusLc83bdrE0KFDmThxIpBcAzj//PPZe++9GTduHE888UTRYyjFzOIbJQ0FBCwAvlCCGKyHxtQO4pbPH57Z+U/9yaOZndusnAwYMIDFixezfv16+vfvz/3338/w4cO3vH7vvffyzDPP8MwzzzBv3jy++MUvMm/evKLGkMuEsoiYGxET08fHRkR9RBwQEWdGxJt5xGBmVq5OOOEE7r77bgBmzpzJ6aefvuW1WbNmcfbZZyOJww47jLVr17Jy5cqitt8rag2ZmWXt3379JEv+XtwhzzG7D+LyT4ztdL/TTjuNb33rW0ycOJGFCxcyZcoUHnroIQBWrFjBHnvssWXfuro6VqxYQW1tbdHidIkJM7MSGzduHMuXL2fmzJlMmDDhA69FxIf2L/adTu4RmJlBQZ/cszRp0iSmTp3K3LlzWbNmzZbtdXV1/O1vf9vyvLm5md13372obbtHYGZWBqZMmcJll11GfX39B7ZPmjSJ66+/nojgT3/6EzvuuGNRh4XAPQIzs7JQV1fHBRdc8KHtEyZM4J577mHvvfdm++23Z8aMGUVv24nAzKyE3nzzwzdOjh8/nvHjxwPJ9YCrr7460xg8NGRmVuWcCMzMqpwTgZlZlXMiMDOrck4EZmZVzonAzKzKORGYmZVQZ2WoZ82axbhx42hoaKCxsZGHH3646DF4HoGZWQl1Vob6uOOOY9KkSUhi4cKFnHLKKSxbtqyoMbhHYGZWYh2Vod5hhx22FJl76623Mlla0z0CMzOAey+GlxYV95zD6uGE73a6W0dlqAHuvPNOLrnkElatWrUlYRSTE4GVrSUr12W+UpnXRbZy0FEZaoCTTjqJk046iQcffJBvfOMb/O53vytq+04EVpYmNwzvfKceall32YnAgII+uWepvTLUrR199NE899xzvPLKKwwZMqRobWeeCCT1AZqAFRExUdIo4GZgF+AJ4KyI2JB1HNa7nHHoiMzfoL0uspWTKVOmsOOOO1JfX8/cuXO3bH/22WfZa6+9kMQTTzzBhg0bGDx4cFHbzqNHcAGwFBiUPr8SuCoibpb038BngR/nEIeZWdlqrwz17bffzvXXX0/fvn3p378/t9xyS+9aoUxSHXAi8B3gQiXRHwucke5yHfBNnAjMrEp1VoZ62rRpTJs2LdMYsr599IfAvwDvpc8HA2sjYlP6vBloczBY0nmSmiQ1rV69OuMwzcyqV2aJQNJEYFVEzG+9uY1dP7wyMxAR0yOiMSIahw4dmkmMZmaW7dDQkcAkSROAGpJrBD8EdpK0bdorqAP+nmEMZmbWicx6BBFxSUTURcRI4DTg9xHxz8ADwMnpbucAs7KKwczMOleKEhPTSC4cP0tyzeBnJYjBzMxSuUwoi4i5wNz08fPAIXm0a2ZmnXPROTOzEuqsDDXA3LlzaWhoYOzYsXz0ox8tegwuMWFmVkKdlaFeu3YtX/rSl7jvvvsYMWIEq1atKnoM7hGYmZVYR2Wob7rpJj75yU8yYkRScmXXXXctevvuEZiZAVc+diXLXi3ugi+jdxnNtEM6nxXcURnqp59+mo0bNzJ+/HjeeOMNLrjgAs4+++yixulEYGZWYh2Vod60aRPz589nzpw5rF+/nsMPP5zDDjuMfffdt2jtOxGYmUFBn9yz1F4Z6rq6OoYMGcKAAQMYMGAARx99NH/5y1+Kmgh8jcDMrAxMmTKFyy67jPr6+g9snzx5Mg899BCbNm3i7bffZt68eey///5Fbds9AjOzMtBeGer999+f448/nnHjxrHNNtvwuc99jgMOOKCobTsRmJmVUGdlqAEuuugiLrroosxi8NCQmVmVc4+gwtw070VmLViRaRtLVq5jTO2gznc0s17BPYIKM2vBii2LsmdlTO2gXBaXN7N8uEdQgcbUDuKWzx9e6jDMrJdwIjDLUB5DdS0mNwznjENH5NKWVRYPDZllKI+hOkiu2+SVcKzyuEdglrE8hupO/cmjmZ7fsiOJM888kxtuuAFISkrU1tZy6KGHMnv2bF577TWmTJnCc889R01NDddee63nEVgZaJoBi27Lvp36k6Hx3OzbMSuhzspQX3HFFTQ0NHDnnXeybNkyvvzlLzNnzpyixuChIeu6RbfBS4uybeOlRfkkG7My0FEZ6iVLlnDccccBMHr0aJYvX87LL79c1PYz6xFIqgEeBPql7dwWEZdL+jnwUeD1dNfPRMSCrOKwjAyrh3Pvzu78M07M7txmbXjpiit4d2lxy1D32380wy69tNP9OipDfeCBB3LHHXdw1FFH8dhjj/HXv/6V5uZmdtttt6LFmeXQ0LvAsRHxpqS+wMOS7k1fuygi/HHPzIyOy1BffPHFXHDBBTQ0NFBfX89BBx3EttsW9607s0QQEQG0FNHom35FVu2ZmfVEIZ/cs9ReGepBgwYxY8YMACKCUaNGMWrUqKK2nek1Akl9JC0AVgH3R8S89KXvSFoo6SpJ/do59jxJTZKaVq9enWWYZmYl114Z6rVr17JhwwYArrnmGo4++mgGDSpuiZdME0FEbI6IBqAOOETSAcAlwGjgYGAXoM3VICJiekQ0RkTj0KFDswzTzKzk2itDvXTpUsaOHcvo0aO59957+dGPflT0tnO5fTQi1kqaCxwfEd9LN78raQYwNY8YzMzKUWdlqA8//HCeeeaZTGPIrEcgaaikndLH/YGPAcsk1abbBPwTsDirGMzMrHOd9ggkNQL/A9gdWE/yxv27iHi1k0Nrgesk9SFJOLdGxGxJv5c0FBCwAPhCT74BMzPrmXYTgaTPAOcDLwDzgaeAGuAoYJqkxcA3IuLFto6PiIXAQW1sP7bnYZuZFUdEkAxQ9F7JTZrd11GPYABwZESsb+tFSQ3APkCbicDMrNzV1NSwZs0aBg8e3GuTQUSwZs0aampqun2OdhNBRFzdSeOeDWxmvVpdXR3Nzc309lvUa2pqqKur6/bxhVwj+L/Av5NcH7gPOBD4akT8otutmpmVgb59+xZ9clZvVMhdQx+PiHXARKAZ2Be4KNOozMwsN4Ukgr7pvxOAmQXcLWRmZr1IIRPKfi1pGcnQ0JfSWz/fyTYsMzPLS7s9gpaJXxFxMXA40BgRG4G3gcn5hGdmZlnrqEdwraSdgbkkF4kfBoiIt4C3sg/NuuO4t+/hyPUPwIwds2vkpUXJegRmVhE6un30hHRxmfHAScD3JL1IkhTua28imZXWkesfYOTG52ljLl/xDKtPlpG0srJk5bpc1i6e3DCcMw4dkXk7lp8OrxFExDukb/wAkkYBJwD/KWlYRBySfYjWVcv77snYLFcPs7IzuWF45zsVwZKV6wCcCCpMwdVHJQ0iWV7y5vTrwyXzzKwkzjh0RC5vznn0OCx/hUwo+zzwLZK7hloKWkRE7JllYGZmebhp3ovMWrAil7bKdVitkB7BVGBsRLySdTBmZnmbtWAFS1auY0xtcVf92lo5D6sVkgieI7ll1MysIo2pHcQtnz880zbKeVitkERwCfBHSfOAd1s2RsT5mUVlZma5KSQR/AT4PbAIeC/bcMzMLG+FJIJNEXFh5pGYmVlJFJIIHpB0HvBrPjg05OJzZp3IZaZ3i/qTofHc7NuxilNIIjgj/feSVtsC6PD20XRW8oNAv7Sd2yLi8nRS2s3ALsATwFkRsaGrgZv1BrnM9Iak7Ac4EVi3dJoIIqK7qza8CxwbEW9K6gs8LOle4ELgqoi4WdJ/A58FftzNNszKXi4zvWecmO35raJ1tHj9URHxcAevDwJGRMTitl6PZDXlltnHfdOvAI7l/V7GdcA3cSKwtry0KNM3uMvWvM4j/Y8hKa5rVr066hF8Kl2m8j5gPrAaqAH2Bo4BPgJ8vaOTS+qTHrs3cDXJnIS1EbEp3aUZaLNISnpd4jyAESPKbwKGZSyHonbJkI2ZdVR99GtpGeqTgU8DtSRlJpYCP+mot9DqHJuBBkk7AXcC+7e1WzvHTgemAzQ2Nra5j1WwxnMzH+9efsVRmZ7fbGvdqRA7ZvdBXP6JsRlFlOis+uhrwE/Tr26LiLWS5gKHATtJ2jbtFdQBf+/Juc0slfFQGuQ3nJZn/Z88yktAfhViu6Pg6qNdlS5puTFNAv2BjwFXAg+Q9DJuBs4BZmUVg1nVyGl9iLyG0/Kq/wNJeYk83qTzqhDbHZklApKhpOvS6wTbALdGxGxJS4CbJf078GfgZxnGYFYdchhKg3yH0/Ko/5Orphmw6LauHzesHk74bvHjaSWzRBARC2nj5umIeB7wgjZmVl0W3Va2y7wWlAgkHQGMbL1/RFyfUUxmZpVpWD2U4eqBhSxMcwOwF7AA2JxuDsCJwMysAhTSI2gExqQTxMzMrMJsU8A+i4FhWQdiZmalUUiPYAiwRNJjfLD66KTMojIzs9wUkgi+mXUQZmZWOoVUH/2DpN2Ag9NNj0XEqmzDMsvH2xs2Z7qW7NQNm9l+uz6Znb9S5bqOQ17K9NZRKOAagaRTgMdI6g2dAsyTlM80RrMMDdmhX+Zv0ttv14chO/TLtI1K9P46DhVkWH1uM8C7qpChoX8FDm7pBaSlI34HdGOKnFn52G1gDbsNrOGWczOcvVpJn2hzlss6DgYUdtfQNlsNBa0p8DgzM+sFCukR3CfpN8DM9PmpwD3ZhWRmZnkq5GLxRZI+BRwJCJgeEXdmHpmZmeWioFpDEXE7cHvGsZiZWQl0tGbxwxFxlKQ3+OAqYiJZkjj7QuFmZpa5jpaqPCr9d2B+4ZiZWd4KmUewl6R+6ePxks5P1yA2M7MKUMg1gtuBRkl7k6wmdhdwEzAhy8BKIY91Uic3DC/b5erMrDoVMh/gvXSh+ZOAH0bE10iWoaw4LeukZmXJynW5LchtZlaoQnoEGyWdTrLQ/CfSbX07O0jSHiSL1wwD3iO57fRHkr4J/C9gdbrrpRFRNvMSslwnNcuaNmZm3VVIIjgX+ALwnYh4QdIo4BcFHLcJ+HpEPCFpIDBf0v3pa1dFxPe6F7KZmRVTIRPKlgDnt3r+AvDdAo5bCaxMH78haSkwvPuhmplZFtq9RiDp1vTfRZIWtvpaJGlhVxqRNBI4CJiXbvpKeq5rJe3czjHnSWqS1LR69eq2djEzsyLoqEdwQfrvxJ40IGkHkjuPvhoR6yT9GPg2ySS1bwPfB6ZsfVxETAemAzQ2NnZrveR/+/WTLPl74Rd/l6xcx5haz5Mzs+rS0YSylenDbYCVEfEOgKT+wG6FnFxSX5IkcGNE3JGe9+VWr/8UmN290ItvTO0gJjd49MrMqkshF4t/CRzR6vnmdNvBbe+ekCSSeQdLI+IHrbbXtkoyJwGLuxRxF1z+ibFZndrMrGIUkgi2jYgNLU8iYoOk7Qo47kjgLGCRpAXptkuB0yU1kAwNLQc+37WQzcysmApJBKslTYqIuwAkTQZe6eygiHiYpEDd1spmzoCZmRWWCL4A3CjpapJP8c3A2ZlGZWZmuSlkHsFzwGHp3T+KiDeyD6tEmmbAouyWYr5szes80v8YIMM1cs3MuqiQ6qO7SfoZ8Mt0YtgYSZ/NIbb8LboNXlqU2elHbnyeI9c/kNn5zcy6o5ChoZ8DM4B/TZ8/DdxCckdQ5RlWD+fencmpl19xVCbnNTPriUKqjw6JiFtJCseRViLdnGlUZmaWm0J6BG9JGky6XKWkw4DXM42qWO69uGtDPS8tSnoEVj1eWgQzTsz2/P6dsjJXSCK4kGQxmr0kPQIMBU7ONKpSGVYP9ZX5rVkb8vi/9u+U9QIdJgJJ2wA1wEeB/UjmBTwVERtziK3nTui0SGru3t6wOdN1CaZu2Mz22/XJ7PwVpfHc5MusynV4jSAi3gO+HxGbIuLJiFjca5JAGRqyQ7/M36S3364PQ3bol2kbZlZZChka+q2kTwF3RES3qoBaYreBNew2sIZbzs1wHsGMHbM7t5lVpEKvEQwANktaTzI8FBHhes1mZhWgkJnFA/MIxMzMSqOQHgGSPgkcRXIL6UMR8atMozIzs9wUUmLiv0gKzy0iWTvgC2kBOjMzqwCF9Ag+ChzQcqFY0nUkScHMzCpAISUmngJGtHq+B9ClxevNzKx8FdIjGAwslfRY+vxg4FFJdwFExKSsgjMzs+wVkgguyzwKMzMrmUJuH/1Dd04saQ/gemAYSeXS6RHxI0m7kJSxHkmyZvEpEfFad9rolVzkzMzKTCHXCLprE/D1iNgfOAz4sqQxwMXAnIjYB5iTPq8O9Sdn/ybtImdm1kUFzSPojohYCaxMH78haSkwHJgMjE93uw6YC0zLKo6y4iJnZlaG2u0RSJqaDu/0mKSRwEHAPGC3NEm0JItd2znmPElNkppWr15djDDMzKwNHQ0NDQf+KOlBSV+UNKQ7DaSL3t8OfDUi1hV6XERMj4jGiGgcOnRod5o2M7MCtJsIIuJrJPMHvgGMAxZKulfS2ZIKqj8kqS9JErgxIu5IN78sqTZ9vRZY1ZNvwMzMeqaz9QgiIv4QEV8kmUj2Q+BrwMudnViSSBa4XxoRP2j10l3AOenjc4BZ3QnczMyKo9Cic/XAacCpwBrg0gIOOxI4C1gkaUG67VLgu8Ctkj4LvAh8uqtBm5lZ8bSbCCTtQ/LmfzqwGbgZ+HhEPF/IiSPiYZK1C9pyXBfjNDOzjHTUI/gNMBM4NSJcZM7MABi58flsJ0WmbSzvu2embdj72k0EEdHm/4KkPsBpEXFjZlGZWVl6pP8xAIzNuJ3lfffkkf7HZN6OJToaGhoEfJnkNtK7gPuBrwBTgQWAE4FZlZmz/QT+3+tHMWZDtivVLtmwjjGDB3Fepq1Yi46Ghm4AXgMeBT4HXARsB0yOiAUdHGdmOfvl07/knufvybydN3d+l9o+B8LmozNtZ0ztICY3DM+0DXtfR4lgz4ioB5B0DfAKMCIi3sglMjMr2D3P38NTrz7Ffrvsl2k7r21azn51/ZhxfHVUhakWHSWCjS0PImKzpBecBCwveX3CnbDnBD69b2XcwbzfLvsx4/gZmbZx7n2ulVWJOppQdqCkdenXG8C4lseSCi4VYdYdLZ9ws/TUq0/lkmzMyl1Hdw31yTMQs61l/QnXn27NElmuR2BmZr2AE4GZWZXLbGEaKw1fZDWzrnIiqDB53EbY9HITTS83ZZpw8rgV0swSTgQVKOuLrHn0OvbbZT8m7Dkh0zbMLOFEYF326X0/7WEhswrii8VmZlXOicDMrMo5EZiZVTknAjOzKpdZIpB0raRVkha32vZNSSskLUi/fFuImVmJZXnX0M+B/wSu32r7VRHxvQzbNSsbeU3w87wL64nMegQR8SDwalbnN+sN8qiiCp53YT1TinkEX5F0NtAEfD0iXitBDGa5yWOdALOeyPti8Y+BvYAGYCXw/fZ2lHSepCZJTatXr84rPjOzqpNrjyAiXm55LOmnwOwO9p0OTAdobGyM7KPLXh7jxR4rNrOuyrVHIKm21dOTgMXt7VuJ8hgv9lixmXVVZj0CSTOB8cAQSc3A5cB4SQ1AAMuBz2fVfrnyeHF5eerVpzJdqcw9tPKX151dUL7l2zNLBBFxehubf5ZVe2ZdlUfPqRJ7aFknzxZ5vWnmUbod2DIaUFWJwKzcuYpq1+WV1PJ+08yjp17Oa2Q7EZhZwfJKnuX8plmJnAjMrCzlNQTl6zhOBGZWhvK8rlKJ13G6yonAytJrt9zKutntTjMpmkETJ7Lzqadk3k7W8vp5QT4/M1+/yZfLUFtZWjd7Nu8sW5ZpG+8sW5bbm2fW8vh5QWX9zOx97hG0kvX9xB6L7Jqa0aP5yA1bF68tnr+edXZm5y6FrH9eUHk/M0u4R9BK1jN/PRZpZuXIPYKt9PaZvx5bN7OuciKoMC1jxTWjR2fWxtuPP87bjz+eacLJ+nsws/dVdCK48rErWfZq4RfQKmUMP+ux4jx6HTWjRzNo4sRM2zCzREUngq7yGH5hdj71FA8LmVWQik4E0w6ZVuoQzMy26M5s6dG7jM78vayiE0G5yWNIxWPrZuWpnEcbnAhylMeFXI+tm5Wncp4t7USQszwm/Vjh3lm2LPNJUr7V1sqdE4FVrTx6TnncagseErSecSKwqpXH3U95TfDzkKD1RJZrFl8LTARWRcQB6bZdgFuAkSRrFp8SEa9lFYNZqflWW+sNsqw19HPg+K22XQzMiYh9gDnpczMzK6EsF69/UNLIrTZPBsanj68D5gK+2d+sF8njAjv4Inue8r5GsFtErASIiJWSds25fTPrgbyuQ7SsreBEkI+yvVgs6TzgPIARI0aUOBozg/yueeS57kGlre7WHXmvR/CypFqA9N9V7e0YEdMjojEiGocOHZpbgGZWXby6W/49gruAc4Dvpv/Oyrl9M7MPqfbV3bK8fXQmyYXhIZKagctJEsCtkj4LvAiU1XzrrLuInvRjZuUoy7uGTm/npeOyarOnsq4F5Ek/ZoXL6+4kf0Ar44vFxfDSFVfw7tLCx/5afiFcC8istPL8wOQPaBWeCLrKvxBm5cEzsvNV0Ylg2KWXljoEM7Oyl/fto2ZmVmacCMzMqlxFDw2ZmZWT7twJ1W//0ZkPczsRmJnloJxvRHEiMDPLQTnfCeVrBGZmVc6JwMysyjkRmJlVOScCM7Mq50RgZlblnAjMzKqcE4GZWZVzIjAzq3KKiFLH0ClJq4G/ZtjEEOCVDM9fDI6xOBxjcTjG4sg6xo9ERKeLvveKRJA1SU0R0VjqODriGIvDMRaHYyyOconRQ0NmZlXOicDMrMo5ESSmlzqAAjjG4nCMxeEYi6MsYvQ1AjOzKucegZlZlavoRCDpWkmrJC1u5/XRkh6V9K6kqVu99jVJT0paLGmmpJoSxfjPkhamX3+UdGCr146X9JSkZyVdnEV8PYlR0h6SHpC0NP1ZXlBuMbZ6vY+kP0uaXY4xStpJ0m2SlqU/z8PLMMZy+ZuZnMa3QFKTpKNavXaOpGfSr3OyiK8nMUpqSN+TnkxfPzWrGD8gIir2Czga+AdgcTuv7wocDHwHmNpq+3DgBaB/+vxW4DMlivEIYOf08QnAvPRxH+A5YE9gO+AvwJgyi7EW+If08UDg6XKLsdXrFwI3AbNL+PvYbozAdcDn0sfbATuVU4xl9jezA+8Pe48DlqWPdwGeT//dOX28c5nFuC+wT/p4d2BlVv/Xrb8qukcQEQ8Cr3bw+qqIeBzY2MbL2wL9JW0LbA/8vUQx/jEiXkuf/gmoSx8fAjwbEc9HxAbgZmByOcUYESsj4on08RvAUpI3jLKJEUBSHXAicE0WsfU0RkmDSN5YfpbutyEi1pZTjKly+Zt5M9J3UmAA0PL4fwL3R8Sr6fdwP3B8OcUYEU9HxDPp478Dq4BOJ4T1VEUngu6KiBXA94AXSTLy6xHx29JGBcBngXvTx8OBv7V6rZmM3mS7qHWMW0gaCRwEzMs5nrZsHeMPgX8B3itNOG1qHeOewGpgRjp8dY2kAaULbYstMZbb34ykkyQtA+4GpqSby+pvpp0YW79+CEnv77msY3EiaIOknUk+XY8i6Z4NkHRmiWM6huQPb1rLpjZ2K+ktYG3E2LJ9B+B24KsRsa4UsbWK5QMxSpoIrIqI+aWMq7U2fo7bkgwz/DgiDgLeAjK7JlSINn6OZfU3ExF3RsRo4J+Ab6eby+pvpp0YAZBUC9wAnBsRmX9AcSJo28eAFyJidURsBO4gGRstCUnjSIYtJkfEmnRzM7BHq93qyKgrXoh2YkRSX5IkcGNE3FGq+NJY2orxSGCSpOUkw2vHSvpFiULs6P+6OSJaelO3kSSGkmgnxrL6m2mRDtHsJWkIZfY302KrGFuGAu8G/k9E/CmPGJwI2vYicJik7SUJOI5kfDt3kkaQ/FGdFRFPt3rpcWAfSaMkbQecBtxVTjGmP7ufAUsj4geliK1VLG3GGBGXRERdRIwk+Rn+PiJK8km2gxhfAv4mab9003HAkhKE2NHvYzn9zeydxoCkfyAZXlkD/Ab4uKSd0x7Mx9NtZRNj+rd8J3B9RPwyr3i2zauhUpA0ExgPDJHUDFwO9AWIiP+WNAxoAgYB70n6KsldLfMk3QY8AWwC/kxGMwA7ixG4DBgM/Ff6e7MpIhojYpOkr5D8IvcBro2IJ8spRpJP22cBiyQtSE93aUTcU0Yx5qaHMf5v4Mb0jeJ54NxyirHM/mY+BZwtaSOwHjg1vTCuFtzJAAAEIElEQVT7qqRvk3yIAvhWRLR7QbcUMUo6heTGgMGSPpOe7jMRsYAMeWaxmVmV89CQmVmVcyIwM6tyTgRmZlXOicDMrMo5EZiZVTknAjOzKudEYGZW5ZwIrOxIuiqd3Nfy/DeSrmn1/PuSLixym28W+Xw7SfpSq+cj26tN38ax/SX9QVKfIsSxnaQH04qgZm1yIrBy9EfSOjWStgGGAGNbvX4E8EgJ4uqKnYAvdbpX26YAd0TE5p4GkZYonwPks8CJ9UpOBFaOHuH9gmVjgcXAG2mNmH7A/sCfJf1K0vx0NafzWg6WdOVWn8a/Kenrks6U9JiSVaF+0tYn7vb2ST/RL5X007S930rqn772DSUrh92vZGWuqcB3SQqJLZD0H+np+7R1fBv+GZjVKqa5LXWGJA1WsgLYyLTNa9LnN0r6mKRHlKy+dUir8/0qPadZm5wIrOykC3JsSgucHQE8SrKOweFAI7Aw/aQ7JSL+Md12vqTB6Slu5oOfgE8hqSl1KnBkRDQAm9nqzVHS/p3ssw9wdUSMBdYCn5LUSFI35iDgk2kskJSJfi4iGiLiovaO3/p7T2sJ7RkRy1tt3ht4Jn08DljUavuP0m2jgTOAo4CpwKWtjl9MshKfWZs8bmjlqqVXcATwA5IFRI4AXicZOoLkzf+k9PEeJG+0ayLiz5J2lbQ7yepOrwH1wD8Cj6fF0vqTrP7U2nGd7PNCq+Jf84GRJMNWsyJiPYCkX3fwPbV1/NaGkCQJ0vN9BFjRqib9OGBhq/MtSvd7EpiTFi5b1PrcEbFZ0gZJA9OV4sw+wInAylXLdYJ6kk+0fwO+DqwDrpU0nqQG/uER8bakuUDrxdJvA04GhpH0EARcFxGXdNBmZ/u82+rxZpJE0dZiJ+1p6/itreeD30cD77/xQ5KobmnjfO+1ev4eH/7b7ge804VYrYp4aMjK1SPARODViNiclgveiWR46FFgR+C1NAmMBg7b6vibSdYXOJkkKcwBTpa0K4CkXdJP260Vss/WHgY+IalGyUpsJ6bb3wAGdvWbTtfS7SOpJRkcSJoYJO1DsgrYonYOb1M6ZNayYIzZhzgRWLlaRDJM8qettr0eEa8A9wHbSlpIsszfB1ZyStdmGEgyrLIyIpYA/wf4bXrM/UDtVsd0us/WIuJxkgWB/kKyYEtTGuMa4JH0Qu5/dHSONvyWZKwfkh7BNpL+QrIWwFLgnC6e7xig6GtAWOXwegRmPSRph4h4U9L2wIPAeRHxRA/OdxBwYUScJelZ4KCejO1LugO4JCKe6u45rLL5GoFZz02XNIZkCOe6niQBgPRi9wOSdgTe62ES2A74lZOAdcQ9AjOzKudrBGZmVc6JwMysyjkRmJlVOScCM7Mq50RgZlblnAjMzKqcE4GZWZVzIjAzq3L/H4Ihl8gem+Q4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(xrv1, rv1, label=\"M0\", drawstyle=\"steps-mid\")\n",
    "plt.plot(xrv2, rv2, label=\"M3\", drawstyle=\"steps-mid\")\n",
    "plt.plot(xrv3, rv3, label=\"M6\", drawstyle=\"steps-mid\")\n",
    "plt.plot(xrv4, rv4, label=\"M9\", drawstyle=\"steps-mid\")\n",
    "plt.legend()\n",
    "plt.ylabel(\"RV precision (m/s)\")\n",
    "plt.xlabel(\"Wavelength ($\\mu$m)\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}