{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Scripts for Calibration of Cernox Sensors\n", "---------------------------------------\n", "\n", "Make a copy of this notebook for an other run" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import math\n", "\n", "from zcalib import read_curve, convert_res, compare_calib, make_calib, logrange, Sensor, CalibRun, nplog, npexp" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('lsdat', 'READ /afs/psi.ch/project/SampleEnvironment/SE_internal/Thermometer_calibs/2012/73027 Cernox 5/X75610.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c1.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c1.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c1.dat')\n", "{'selected': 2, 'averaged': 46}\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c2.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c2.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c2.dat')\n", "{'selected': 2, 'averaged': 46}\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c4.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c4.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c4.dat')\n", "{'selected': 2, 'averaged': 46}\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c6.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c6.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c6.dat')\n", "{'selected': 2, 'averaged': 46}\n", "calibration file has 199 points\n" ] } ], "source": [ "run = CalibRun([\n", " Sensor(3, 'X75610'), # the reference sensor must be the first\n", " Sensor(1, 'X133979', 'CX-1050-SD'),\n", " Sensor(2, 'X133978', 'CX-1050-SD'),\n", " Sensor(4, 'X133928', 'CX-1050-SD'),\n", " Sensor(6, 'X133981', 'CX-1050-SD'),\n", " ],\n", " t_points = (1.0, 1.2) + logrange(1.4, 310, n=195) + (330,),\n", " caldate = '2018-10-25',\n", " logT = False,\n", " logR = True,\n", " calib_data_file = 'calib_data/calib%s_p%d_c%d.dat',\n", " outputpath='%s/%s.340')\n", "# use slightly lower smoothtst as we measured less points (1e-7 for 60, 0.8e-7 for 48 and 0.4e-7 for 24 points)\n", "run.make(diflim=0.001, smoothref=1e-7, smoothtst=0.8e-7)\n", "print('calibration file has %d points' % len(run.t_points))" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('z340', 'READ 2018-10-25/X133979.340')\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c1.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c1.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c1.dat')\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD/CAYAAAD12nFYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4FNX+BvD3ZDeEmoTQIULoTSQooKhAaIKIIIo/QUVA\nRAWxcC0gcgW8egVsXAXFQrWAqKAogtIiRTqCdAJSQ09ogZCy+/7+mE3IhpRNb+/nefKwM3POzNll\ndr57yswxJCEiIpLAK68LICIi+YsCg4iIuFFgEBERNwoMIiLiRoFBRETc2PO6AFlljNGwKhGRTCBp\nUlpfKGoMJPPt3+jRo/P1/jOTPyN5PEmbVprMbMvpzzyv/890TuicyMr+Bw8mJk1K+/d0oQgM+VlI\nSEi+3n9m8mckjydp00qT2W35mc6J9NPqnMi5/cfEAD4+aac3ZMFuiTHGsKC/B8leY8aMwZgxY/K6\nGJKP6Jy45tFHgc6dgcceM2BhbkoSSaqg/mqUnKNz4hpPagwKDFLo6CIgyemcuMaTwFDgRyWJiKQk\nKCgIhw8fzuti5LkaNWrg0KFDicsKDCJSZB0+fBjqfwSMce9GUFOSiIi4iY0FihVLO40Cg4hIEaIa\ng4iIuFFgEBERNwoMIiIFxIABA/D6669j9erVaNiwYeL6ffv2oVmzZvDz88OkSZOyfBwFBhGRAubO\nO+/E7t27E5cnTJiA9u3b48KFCxg6dChCQ0PRvn17+Pv7o1atWhnevwKDiEgBd/jwYTRu3DhxuVSp\nUhg4cCDefffdTO1PgUFEJJ/666+/cMstt8DPzw+9e/fG1atXAQB//PEHbrjhBgBAhw4dsGLFCjzz\nzDPw9fXF/v370aJFCzzyyCOoWbNmpo6rwCAikg/FxcWhZ8+e6NevHyIjI/Hggw/ihx9+SNyecFPa\nsmXL0Lp1a0yePBkXL15EnTp1snxsT+5j0J3PIlIkmRSfK5pxmbm5et26dYiPj8dzzz0HAHjggQfQ\nokWL7ClQGuLjrX/t6Vz5FRhEpEjKy6dlHD9+HNWqVXNbV6NGjRw/rifNSEA2NSUZY7oYY/YYY/YZ\nY4ankuZDY0yYMWarMSbY07zGmBeNMU5jTEB2lFVEJK9VqVIF4eHhbuuOHDmS48fNtcBgjPECMAlA\nZwCNAfQxxjRIluZuALVJ1gXwFIApnuQ1xgQC6ARAj0gUkUKjVatWsNvt+OijjxAfH4958+Zhw4YN\nidvTevgfScTExCA2NhZOpxMxMTGIi4vz6Li5WWNoCSCM5GGScQDmAOiRLE0PALMAgOR6AH7GmEoe\n5P0AwMvZUEYRkXzD29sb8+bNw/Tp01GuXDl89913eOCBBxK3J30iavKno65cuRIlSpRAt27dcPTo\nUZQsWRKdO3f26LieBobs6GOoBuBokuVjsC746aWpllZeY0x3AEdJbk/+wYiIFHQ333wztmzZkuK2\npM1Ky5cvd9vWtm1bOJ3OTB0zNwNDZqR5pTfGlAAwElYzUrp5+vfvj6CgIACAv78/goODE2dsCg0N\nBQAta1nLRXBZrgkNDcX778/A6dPAmDFBaaY1WZ3IwhhzG4AxJLu4lkcAIMnxSdJMAbCC5Leu5T0A\n2gKomVJeAAsBLAVwBVZACAQQDqAlydPJjk9NxiEiyRljNFEP3D+HTZuAp54CNm9OXJ/iD+7s6GPY\nCKCOMaaGMaYYgN4AFiRLswDAY65C3gbgPMlTqeUluYNkZZK1SNaE1cTULHlQEBERz+VaUxJJhzFm\nKIDfYQWaqSR3G2OesjbzM5K/GmO6GmP2A7gMYEBaeVM6DNJpfhIRkbTlah8DycUA6idb92my5aGe\n5k0hTcYfISgiIm5y9QY3ERHJ/xQYRETEjQKDiIi4UWAQESlAcmNqTwUGEZECKL2pPd999100adIE\nvr6+qF27doZmcvNkLgZAgUFEJF9LPrUnAHz55Zc4f/48Fi1ahEmTJmHu3Lke7Us1BhGRfCyzU3u+\n9NJLCA4OhpeXF+rVq4cePXpgzZo1Hh1TgUFEJJ/Kzqk9V61adV2NIjX5/SF6IiJ5yozNnocpcHTG\nn8eUXVN7jh49GiQxYMAAj9LHxABlyqSfToFBRIqkzFzQs0t2TO05adIkfPXVV1i9ejW8vb09yhMT\nA5Qvn346NSWJiOSyrE7tOW3aNEyYMAHLly9HlSpVPM6nPgYRkXwqK1N7fv3113jttdewZMmSDNcy\nFBhERPKprEzt+e9//xuRkZFo0aIFypQpA19fXwwZMsSj43p6H0OWJ+rJa5qoR0RSool6LEk/h4ce\nAnr2BHr3zvmJekREpABQU5KIiLhRYBARETcKDCIi4kaBQURE3CgwiIiIm5gYPXZbRESSiI1VjUFE\nRJJQU5KISAGiqT1FRCRF6U3tOXHiRNSuXRu+vr6oXLkyHn/8cURFRXm0bwUGEZFCIPnUnj169MCm\nTZtw8eJF7NmzB4cPH8Zbb73l0b4UGERE8rHMTu1Zs2ZNlC1bFgDgcDjg5eXl8aO3FRhERPKprE7t\nOXv2bPj5+aFixYqoWLFi4kxwaXE4AKcTsHswPZsCg4gUTcZkz18mJJ3a02azZXhqzz59+uDChQvY\nt28fdu3ahYkTJ6abJ+EeBk+KrMAgIkUTmT1/mZAdU3sCQO3atTFixAjMmjUr3bSe3sMAKDCIiOS6\nrE7tmVRcXBxKliyZbjpP+xcABQYRkVyXlak9p06dijNnzgAAdu3ahXHjxrnN/pYaBQYRkXwsK1N7\nrlmzBk2aNIGvry/uv/9+9OvXD8OGDUv3mBkJDJraU0QKJU3taUn4HHbssKb23LnTbb2m9hQRKarU\nlCQiIm5yPTAYY7oYY/YYY/YZY4ankuZDY0yYMWarMSY4vbzGmAnGmN2u9D8YY3yzo6wiIkWRp3Mx\nANkQGIwxXgAmAegMoDGAPsaYBsnS3A2gNsm6AJ4CMMWDvL8DaEwyGEAYgFezWlYRkaIqt+9jaAkg\njORhknEA5gDokSxNDwCzAIDkegB+xphKaeUluZSk05V/HYDAbCiriEiRlNtNSdUAHE2yfMy1zpM0\nnuQFgMcBLMpySUVEiqiMBAYPHqeUIzx+wIgx5jUAcSS/SS1N//79ERQUBADw9/dHcHAwQkJCAACh\noaEAoGUta7mILdeoUeO6ewCKooQnr27dGootW2agf38kXi9Tk+X7GIwxtwEYQ7KLa3kEAJIcnyTN\nFAArSH7rWt4DoC2AmmnlNcb0BzAIQHuSMakcX/cxiIikY+pUYM0aYNo0azmn72PYCKCOMaaGMaYY\ngN4AFiRLswDAY67C3AbgPMlTaeU1xnQB8DKA7qkFBREp4q5eBZ59Fpg7N69Lku/lalMSSYcxZiis\nUUReAKaS3G2MecrazM9I/mqM6WqM2Q/gMoABaeV17fojAMUALHFVB9eRHJLV8opIIXHlCtCmDXDx\nIhARAfzf/+V1ifK1XO9jILkYQP1k6z5NtjzU07yu9XWzo2wiUkht3mz9O28ecP/9eVuWAiBX72MQ\nEckTu3YBTZsCDRoAx48DFy7kdYnyNc3HICKF386dQOPG1lyVwcHAli15XaJ8Tc9KEpHCLyEwAEDz\n5sCmTXlbnnxOgUFECr+dO4FGjazXt9yiwJAOBQYRKdwiIoDoaCDQ9aQc1RjSpcAgIoXbrl1WbSHh\nzuZ69YBTp9QBnQYFBhEp3JI2IwGAzQYEBQGHD+dZkfI7BQYRKdySdjwnqF4dOHIkb8pTAOg+BhEp\n3PbsARo2dF9Xo4YCQxp0H4OIFG4HDwK1armvq15dTUlpUFOSiBReTidw7JgVCJJSU1KaFBhEpPA6\neRLw8wNKlHBfr8CQJgUGESm8Dh+2+hOSU2BIkwKDiBRehw5ZQ1OTq1YNOH0aiIvL7RIVCAoMIlJ4\npVZjsNuBSpWA8PDcL1MBoMAgIoVXaoEBUHNSGnQfg4gUXgoMmaL7GESk8EqtjwFQYEiDmpJEpHAi\n064x6O7nVCkwiEjhFBFhNZT7+qa8XTWGFF29av2rPgYRKXzSqi0AeixGKiIjgXLlrj2lPD0KDCJS\ncKTVvwBcqzGQuVWiAiEyEggI8Dy9AoOIFBzp1Rj8/AAvL+D8+dwrUwGgwCAihVd6gQFQB3QKFBhE\npPDyJDCoA/o6CgwiUnil18cAKDCkQIFBRAovT2sMGpnkpkgGhuPzv8zrIohITrt40XquQ7lyaadT\njeE6RTIw+PTtj2PzZ+V1MUQkJx0+bDUjpTcYX53P14mISD+eJlUoAsP6iS+hRN8BODp/Zl4XRURy\nyqFD6TcjAaoxpKBI1hi6PjEeGye+jJJ9H8eh+dPzujgikhM86V8AgCpVgDNncnfCnn/+AVq1AipU\nANq3t5q98pEiGRgAoMsT47DlfyNQpu8TOKjgIEXFsWOI/uwTnB/7KmJmTrv2UJzCyNPAYLcDlSsD\nx47lfJkA4MABoF07oHdvYOdOoGFDoHNnICoqd47vgSIbGACg08C3sO3DkfDtOxAH5k/L6+KI5Bjn\n4UM4eG9rXGgQhIWTnsXs0I+w5o1BONqwGiI3r8nr4gEHDwITJwKvvGI1cGeHhD4GT+TWyKS4OKBX\nL+DFF4HnnwcqVgQmTbJqLZ98kmIWhwN4/XWgRw/g8cdzJ34UycDwwgvWGweA9o//B9s//Df8+j6B\nffO/yNuCFTaLF1tfgjp1rLP611/zukRFD4mTH/4XF26si+X4B+tWzUG3LVEYsCgKe0aexeymwYhp\n1xbhB7blXRkjIoA77gB27bKaVJo2BVatyvp+Pe1jAIB69YB9+7J+zPS8+641neizz15bZwwwcqQV\nIOLj3ZLHxloVi7VrraAAAB07AufO5VwRY2OtimSZMhnIRLJA/wHg00+TFSqQ//sfGRtLkuTKqaN5\nppThnnmfU7IudtKnPOd7A9+76WP26/Arf+s3g/HVg8inn772oee1uDg6jp/kxW3/0HnmLOl05nWJ\nspXz5EkebNuU26vYOHv2a3Q6nXQ6ydmzyerVybvvJl96iZxUP4S/1irN3Qcj8qagffqQL7zAA5EH\nODZ0LL8a/yijy/ry3LKFWdtvhQpkeLhnaSdMIIcNy9rx0hMWRpYrRx46lPL2228nv//ebdXzz5P3\n3ENevWotO53kkCFk9+45d7qePGl9dMlZl/9UrqupbSgof643x+3bybvuIuvVIxcssD7k1dPG8kwp\nw50/TMnsZyokIz6ezePFAnnrwLb0eaMkK7xVg7bXi7PcI+25rX4wo+/vRcbF5Vp5nE4n90fs5y/r\nv+KPY/rw926NuK26L6/YDU+XMDzkB0YWByN8vLgksArfaP4Yn+xzgG+9Rc6dSx44kOMFJP/5hzx/\nPtt2eWL25zzr78Npd1Xk7qNbSVrXpU6dyJtuIleuvJY27kosd1SpwMHtbku8AOWaxYvJOnX41uLX\nWG58Ob6w6AWOWjaKo0e04ulShsNea84pG6dw8/HNPHL+CE9cOsEzl8/wXPQ5RsdFp77f8HDrIuzp\n1XPBAitS5qTu3cnx41Pf/u23ZEhI4uLq1WSVKmREsnh99Sp5443kl196fuiMBJFdu8j69a9fn+OB\nAUAXAHsA7AMwPJU0HwIIA7AVQHB6eQGUBfA7gL0AfgPgl8p+3T6shQvJBg3IDh3IbdvItdP/w9Ol\nDHfkcHCIvHCKiz95iQtH9OKysQN4/uj+HD1ebnEcO84I73JsP6gmhy3+Fy/HXiZJRsVE8e2FXzHw\nufb8raaN05v5M/Ctpqz77s1s+UlrPjbvMc7dMZfxjvjsKYfTwcVhizn0y4f54gNluLy+D6NK2Lk5\nuCb/e2NP9rn1HU78fC33HjvFq3FXeer8Ba5YtpyzXxnMxbcE8mxJw+ltuvHBbudZvjz52GPWL6ns\ntvx/w3isrJ2nynjxSglvxtzVkVy7Ns08Jy+dZOjBUC4OW8yN4RsZfjGc8Y54xsTHcMeeVVzTqQEP\nlfXizA8H8mrcVe7dSw4fbl0n33kn5Qrb5WVLeNDPzi5Dvs3+N5mWjh254j8DWe+jejwaeZJXrlzb\nFL3oZ0YH+HLWwBbs+EZdtnu5Ip941I+j7inJhx8tyRJji/HOaXfyw3Uf8vjF4+77/eknsnNnz8sR\nFkYGBaWdJiqKfPVV8uabyfbtraqXp1fcJUvIWrWYZuSNjiZLlyYvXGBMjHVx/uGHlJNu2pRChSg2\nlvz9d7djhIeTzz1Hli1LliljFf3jj8nLl10JfvqJ7NnT+tuxg6QVkG6//fpj5mhggNVPsR9ADQDe\nrgt/g2Rp7gaw0PX6VgDr0ssLYDyAV1yvhwMYl8rxr3vDsbHkRx+RFSuSgwaRyz56k6dLGYY+040H\n1/+e8v9MJl04d5K/PGb9GtpXpyw3drqR61tW4wUf8O/OzegIP+bZjqKirP/hnj2tE65cObJyZbJJ\nE/Khh6xfJmvWpH0iZjenk/sad+eEFjU4cumoVJNt3RDOE5Vq8cM7R7HZPRsZcMtylm7zGSuOvJV1\nP2jMA5GZ/4nucDo4Y8t0PvJsNS5sUZZXSxfnpfu68sqsuRwx9BIrViSnTUv/+7x73S9c1rICT/rZ\n+dfYF/mvFxxs0IA8cSLTRbvOonFP8HRpL/793WTuPbOHw74bxGG9yvBy5XJkr17kvn2JaWPiY/jV\ntq/YZnob+r3txzun3clOszrx5k9vZqV3KjFgpI0v3+XFU742/tauBYc9v5edO5NVq1qnxQsvpN6C\nkeBM61v5TCdffvX9uex7k2nZvp0xlcqz4hvlWee23bTZyJIlyYcfJv/+25Vm61brs6hUiaxTh7zv\nPutKd+utdNSpzSUrZ7DvvL70H+fPttPbcvKGyTx56ST5+uvkyJGelyUujixePMkVMxmn0/pe3Xcf\nuWqVVZW85RayRw/yXDqfV3Q02bAhOW9e+uXo0IFcsIAffJB+BebFUZFs8cgCTt7wMb/d8S0vvfaK\ndRErV46cPZvbt5M33ED+61/kwYNWhXTpUvLee8m6dZw8PORtOgMD6Zg1k/zvf63PNzKSP/1Edut2\n/fFyOjDcBmBRkuURyWsNAKYAeCjJ8m4AldLK66pFVHK9rgxgTyrHT/WDjoy0PsRy5chxT0/n0g4N\neMLXi2GVi3HZ/7XkXz9MZlxs5i60V66Qa2e9y0MBNm66rQYjtvzptv3wkR38smsgz5fx5uXJE1O9\ncjkuR/PiqPGMC6jAs+278Z8P/sdzWzeQp09bPw+2bLHqmM89R2ezmxlfohRP1GrFlfUH8tP67/Cl\nxl9z+B2f872ek/nry5/xxPSv6Vy40Gpby2Lb//EpP3FviUps8UGL9H/579lj/eRZt856/4fJka85\nWSpkMgPeDOTuM7szfPztJ//my0/V4t4bSjKqbhCdH3xA55mz/Oknq029Xz/rY/KU0+nk8tlvc1OQ\nD8Nq+fPdQbPZqJH1Pc+qsPWLGVHSMPyPX9zWrzmyho3eCeKshxoypqwvI1q34KLnuvK+wQHs/V4r\nLtjwFaNPhVvtW8uXk2PGkK1b01mqFMPb9ua91bcyOJh8802rdeTQoQw0I2zaxLP+JVnuvid55kzW\n32N6HAMf54QOlel755ecO9cqZ2SkVaupUIEcNSqd3zWjRllXMKeT0XHR/GnPT3zkh0fo97YfV93o\nx49f7cTBvwxm/x/7s/f3vdljdg/e+829fOrnpzhl4xQevXDUbXex9Rtz6XtbuXo1efFismO9/z55\nyy2MvxzFA5EHuO7oOq4OW85LTw6gs359qykwNcOHk/ffTzqdjIqJ4rhV49huRjvWnFiT9825j9O2\nTOPVONcbfestRj/1PCtUSPwB7ybeEc9f9/3K++bcR9+3fVl6cEe2ffdJPvNOO54pZfj0lG48tmIB\nY/zKs3XZ7fz665SLtObhkdxZJoANhlWjz3982PSTpvyn771kr16cPt2qISeXVmAw1vbMM8Y8AKAz\nySddy48CaEnyuSRpfgbwNsk/XctLXLWAmqnlNcacI1k2yT4iSV434MoYw/TeQ1gYMGoUsHIlcOF8\nPDo3+Bp3eU1Fm/CNqBgVg8031QK69cRtTw6Df/mqKe7j+HFrEM4vvwDrVp/Hy74d0Pv0VrxU9Q34\n3P4auncH7roLKF36Wp44Rxze+2wA7n37B9QKaobiEz7AztK3YtkyYN1vF9B44wz0jXgf2wIqYPRd\nTvxd8wBw1R9On3Owx1ZA2ajbUJW3onR0I0SeKokD5/ajauBO3Oa9CTWjDyDwzAnUjHTAB8XgoA3x\nTgfiTBx8nA7UPO+FgBgv/NGqLvY91g033doW7Wu2QwnvEh79vyI6GifKN8CTD5zHO++vR4PyDdLP\n89NP1uiMjRutkRoAtm0DQl6YDq/2Y7Fr2HpUKl0pzV2QwNq1xIwvRqHf0gkoF10VCxpOwtZq3RAX\nb7B+PVCqFPDxx9bQ8cy4GnsFS97oj+Yf/YCVgZWxvONAPD+qNxqWbwjj6dyHyQq9uVklXL6jJdpM\n/uW6zZdjL2Pmtpn4dt1UtNl9BZ0PeqHZWW+UOhVpDUex260JZgIDgdtvx9U7O+LF+Xdi4YqS+PRT\n67zKTLEAILZnd4yNWoHQsiuwcnZz2GyZ20+6zpxBVI0aaHB/S6wYvQJ167oX+PhxYMgQa6DQp58C\nrVunVNhY4JZbgBEjgEceSVx9NS4aXtUCMfuLFxBVOQAlvEughL1E4rkcfjEcfx77E4v3L0YNvxq4\no8K9+Pv7bvjX9+Owu9GD+LF4b+zZA3S7l+ja+wi8/Zaia9eh6PtqfSxy7kWFkhUQUKwSLkfZEX75\nHwxa68Crq+IwpsEP8G7dEU88YQ2sAgCEhgIPPQT8/Tfmnv0Dzy9+Hm1rtMWA4AGo4V8DW05swaxt\ns7Dj9A68dPtLeCo+GOd7Posx92/Hp59ee6vhF8Mx7a9p+OKvL1CxVEU8efOT6H1jbxzaVwYdOgA7\na9yN0o92xHvBVzF+5Qfo810bvH9+F0rv2uR2kYm4EoFFb/RFuym/4dN3B2PFosGIOhaEIeNWY9yG\nAdj15jlMe/kI9p2rgA8+cP+4jTEgmeKZlVeBYSmAV5CxwBBB8rqnfRhj2K9fPwS5xjf7+/sjODgY\nISEhAIDQ0FAASFz+9ddQHD8O+PmFICwM2Lh0MvyOfo++V7ajxdkIPNOwK2p2fxkhISGIiAAWLAjF\n+vXA2bMh6NwZuMFnDJr8/F/UqVsRdeeswG9rwrFhA7B7dwjWrQMaNgxFy5ZAz54hKFUK+OOPUPy6\n73s0PjETry71xrp4G+JLl8UDMcexr0V9DKp8CFfr1MAbA95Ap1qdsHb1WkRHO4GKlfFH2HosXT0f\nEY4jKFHPB02q1kHx43YE+gbi3s73ok5AHWxfvx3GmMT3t2xZKHbtjcWpqIY4vnk1yq15H3ee3Iad\nTSrjzbvP4caLXTCiy6O4v1uPFD+fhOWqn63ArrVT8cPQjhh4S/9UP8/rlvv2BbZvR8jGjYC3N0JD\nQ3HiBDD0+1CUaPIbpt45Fj72YinmP3sWuK/Hb2hw4jm8feoAtj04Giur3Q6nsaFRoxAYA8TFheKG\nG4B27TwsTxrL0RGnML13L1QNXY/iDbzxcdviKFa9HXrU74G+Pfp6vL+DC2bi9tlfI+ifSKxdvylx\n+/HjwBdfhCIiArjpphAEBwMXLoTCyyv1/X35ZShGjwZatQrBlCnA5s2Zf38AEDpjBmKefQb9BtyA\nfr5bcHfHDammdziAlaHLYUCEdOiQoeP5TP4Fu/ZNQuSQqWhRv1qK6UlgzJhQfPop0KhRCB57DPD2\nDkWFCkDbtiE4cwZYPuljVP7kdZyccBThkSWwbVsoyjnP4L3FQ3Bh32ns3vNHquW5dDkeDw6ZjOX7\n/0TZztsxbPk+nIguhWVtb4AtyAf/nD2KK3vi0WeTL3qbAMzqPhHndl7Grm0lEBMTgjZtgKpVQ+Fd\nbQeKn/wGL0xbhw9qN8Gcky+iVUhfTL1/ITYMegQRLz2DOQ3CsP3UdjxX6Tk0qtDouvKUqVcG49aM\nw++LF+PrmVdxbO54VGhQA6Ghodhycgt2l9qN3jf2RrOrzVC3XF23/Ds3R6P/Sw/iiW4nsXX/JsR4\nn0StZ+fiqS+W4LKzMqJfGYb6zRvgxz0/YsvkqRgZGocWSxejwm3tsWJFKH7+GZg1KwQ9+h5D0E+1\nEFX+YZS+dwZatw7FjBkzAABBQUEYO3ZsqoEhu5qSFidZ9qQpaQ+uNSWlmBeu5iZea0rancrxU65b\nZcLfv/3KE6W9+HLIMwwJsZobX3uNDA0lL5w6zVUDO/FMKcM/3niCTofjuvznz1sDEZ58kmzdmmzR\nwho18swz5JhPtrLp/5qz+fCyfHpEE5YfW4pdv+7KNUfWZFv5U3XsGK90vo8RgQ14V89etL1ShZ8t\n+jPV5BeWbeRpb382e7U6r8ReSTVdiuLjyS5drhsqePKkk2UHPcTar/RhbOz1bSFLlpBNqv3DtdUD\nuPnG8rx4IONNT5k1algUZ906ibHVA7m3eS22GebP//zxH8Y5PBhpFRvLI5VK8PdPXiZpNUtNm0YG\nB5MBAdbQxKefJvv2tZqla9SwWkzCwtx343CQX31Fli9PfvJJ9g5ddD72GH+4txF9+wzmyy9fG0CW\n0Lc5dNBVvlXpfzyCQMbDi+dsAfy79WBe2JpGc0oS8+ZcZXgJH7797rMepY+Ntd7r//2f1Wbu7W39\nVa1qfW4rK97P2be8w+HDybffJr/s9SM3VexCf3+ybVvys8+uH9mzcKHVNderF3nU1aJ0dcZUnuvR\nhasOr+KGYxsYfjGczuhoxlWswh9Gb+N775FTplj9Hyl8nbl/Wyj/bFeH54sbXvCx80hJXw4d0YoB\n4wM4ZsWYxIEYaen28BH+Wb8xp/yrLXvM7sEhvwzh9L+mMyomKvVMv/3GK83v5PffW90YCWX7c89S\nHq1elp/dU5ld/teS817pzvjy5awmyGSOHSNfeYUc0HQAF1UN4PTp159QyOE+BhuudSAXg9WB3DBZ\nmq641vl8G651PqeaF1bnc0KQyFDnc1aELZ3LY2XtXHdTAFeMfJir3nySf9zblKdLe3HVbVV5cGto\nlvZ/6NxxIWrXAAARt0lEQVQhLjmwhOeic6lDMIHTafVVVKjAZT36sNiL5fl/b3513RfCeSmKR33r\nsl/PclwUtihzx4qIsL6l06a5rT4VcYW+L7ZkzYGvcds2awz+/v3k44+TXQN/4nG/Ypz/4E2MjcmG\nRv8MuHjRGka4YXUM+eGHjC9fjtN6VOfT855IN+8/41/l6ro+PHk6lqNHW32qXbpYozZTutj89ZfV\ncVyxonURHDbMWm7YkGzWzBqdku1OnKCjYgV2HxbI+v3eZ40a1ijKgADy/pvCeKrKTTx/x92MXbeZ\nl87Hc9XXh/l9w1E8a8rxmxbv8881zlQD1cKF5KA6Q7i6binGxMdkqngOR7JAuHOn1SmRMNx30CDy\n3//m1avk/PnWxd/X1xqR07072bix1c+6eHGyHW/aZI0DTeqLL6z/oAw4dWwff1z5Obv89z+s2vF7\n7j161qN8f/5JBgaSMW9OsDrYPfXKK1Zne0oOHLAiqo8Peccd6Y54i7twjpeKe3Hppu+u25ajgcHa\nP7rAGlYaBmCEa91TAJ5MkmaSKwhsA3BzWnld6wMALHVt+x2AfyrHTu9jzrDYy5e4ZlQ/ru5Un2tu\nr87lT3fh3pU/Zvtx8kR4ONm1Ky/Ub8A2j1Zhjcdf485dVsdy/PlLDAtsy+kNbuDzC1/M2nF277au\ntrNnu60+HHGC1d64mT6PPMRilcNYttIZTu3Yi2dKefHbtx+lw5nC1TQXTJ1KtmrlukAdO8a49iFc\ndmNpTlv7SeqZoqJ4JqAEhw0axoAAcuBA6217Ii7Ouvdg3DjyvffIZcty+H68r79mTIO6bPZ2EAfO\nHs5vfznFM5/Pp7NCBZ55Zyzn75rHN0Lf4IAfB7D77O7sNbcXh096lGFVb+DX5R7iTY3iOHGi1YF6\n6ZL1PocNIyvW+4sHAry4Z87k7C3vc8+RbdqQ33xjjTQ4634xvnTJuiZ+/z25eTMZk1JMio21zsGE\nXl+HwxozmsIvbE8NH25dj9Mb1+F0WufTjBm0jpfSeNHUNG9O/vFH2mmuXPH4hPn7tpr8dWzf69bn\neGDIy7+cCAyFntNJzpjBuMBA/lnTnyNaBvGjxkO4r2RdzmoUyDuntMn0rz8327dbX8wRI9y+udFx\n0Ry8YCgfHFKJG6oZ/t2oPDesnJP142WBw2H9Ak0c9RETw4v3dORvDby58cDqFPPseqgfv2nkw45d\no9IcxJIvOJ3kiy8y7oZq/HpgS84J9uYRfy+GDC7Bqu9V5d1f3c3hS4bz882fc/7u+fx2x7f8cN2H\nHDqnH1fU9ebSGwPYadA3rF0vhsWLW7cIPPKvrRzdszyP3dEk+8vrcFiR1mazhmln1r//TQ4dar3+\n8UerfTcLEdjhsAZOvfBC2um++86qEToctGo+pUp5dhNoZKR170M2Dktf3aMZl7zQ/br1CgySspgY\nOr6dw3U923Dm7eXYs18Zvrp4JGPjs/ERFydPWp01VapY9/5PmGA1srdoQdauTeesWfnm0RWrVllt\n3lEJzb+xsTx6Vysub1icp88eSUwXF0dOemIzTxb35jP/HZdfiu+Zn38mBwzg1XfH8+TBHR41acZF\nX+bhrndwaz1fNhjpxy5fdWG7Ge3Y4tXyvBJQxrqTNCc4HNYw6Kw4fNi6G2z/fivyf3d9k0pGRUZa\nLaUzZ6a8/cwZ6zxyq5jUq5fkZo40zJ9vdUxmozX3t+Bvz1zffJZWYMjyqKS85slwVckH9u61hrOe\nPg0UL249s751a8DbO69L5qZ3b+v5a2+84VoRF4e/72oKc/QYav+6DhfKNMLI7jvw5u4OmNjTjrdn\nHobdy56nZc4VTicwdiwcH09G2CNd4CxTBg0+nwevUf8Ghg7N69KlrXt3YNky4OmngQkTkB1jdnft\nsoZLz55tncoJnE7gnnuAJk2sQyV65BHraXkDBqS941dftb4fo0dnuYwJ1j50B86XLYG7pyx1W5/W\ncNU8/8Wf1T+oxiDZ6MgRqwN56dJr6xzxcfzmsZt5vqQ3//K5kRdK+vGpPmW45MCSvCtoXtm82XoS\nXL9+5IoVeV0az5w6RR4/nn66DFqxwhpFlvCMo/Pnyd69rRGJ1/VBvP++VWNOT6dOVq0uG617uC0X\nPt76uvVIo8ZQBH7qiHjuhhusX4F9+lg3MzZvDsDYsa/OOowKGYygGvNxwc+GiYN/wZ3V78zr4ua+\nm2+2/gqSihVzZLchIcDy5daT6EeOBC5ftu59mzo1hYpw8+bAnDlp75AENm2ybvLLRl52bzjjMzbp\ngwKDSDLt2llzzNx7L1C3rjVBV5063lg1+wtUrvIZ4hxx8LH75HUxJR9o0gT4+2/gxAnAy8uaHyhF\nzZoBO3ZYE/uk1nx68CBQsqQ1yU82Mt7eYHTGpjlVYBBJQe/eVmBYuhS48UagVq2Ex1J4KSiIGx8f\nDyaWK13aSrRzJxAcnHKaTZtcVdTsZdUY4tNPmDRPtpdCpJAoVcqaqK527cw/q0gkUfPm1sU/NTkV\nGLx9wPiM1RgUGEREckN6gWHz5hwJDDa7N6gag4hIPpRWYHA6rcCQzR3PAGDzLqYag4hIvtS0qXUD\nREzM9dt27wbKlQMqVMj2w3p5FwPiHRnLk+2lEBGR65UsaQ1z2779+m1r1wK3354jh7XZiwEONSWJ\niORPt9yScnPSn3/mXGAo5gOoj0FEJJ9q3hxYv/769WvXAq1a5cgh7d4+oENNSSIi+VP79tbNMUmf\n7xYZCYSHWzfM5ACbt4+akkRE8q369a05vnfuvLZu3TqgZUtrfQ6wefvAqPNZRCSfMga4+25g0aJr\n61atyrFmJMBqSoLDmaE8CgwiIrkpaWCIjwe+/NJ6El8O8S5WHEZ9DCIi+Vj79sDGjcClS8CCBdYz\nlJo2zbHD2b19MhwY9BA9EZHcVKoU0K0b8PDDwPnzwDPP5Ojh7MWKqylJRCTfmzkTqFbNeqb7Aw/k\n6KHs3sXh5XAmTGzmEQUGEZHcVqwY8MknwP791uscZPP2hp1AvNPzIasKDCIiecEY6zEZOc1uRzF6\nIcaRwjOaUqHAICJSmNlssMMLMfEKDCIiAlg1BqdRjUFERFxsNtjphVhHrMdZFBhERAozux12GjUl\niYiIi80Gb6opSUREEthsqjGIiEgSdjvsTqjGICIiLjYb7IQ6n0VExMVuh92ppiQREUlgs8FGNSWJ\niEgCux02J1RjEBERF5sNNidVYxARERebDTZnLnY+G2PKGmN+N8bsNcb8ZozxSyVdF2PMHmPMPmPM\n8PTyG2M6GmM2GWO2GWM2GmPaZaWcIiJFlt1u1RhysSlpBIClJOsDWA7g1eQJjDFeACYB6AygMYA+\nxpgG6eQ/A6AbyaYA+gP4MovlFBEpmmw2eOVyU1IPADNdr2cCuC+FNC0BhJE8TDIOwBxXvlTzk9xG\n8qTr9U4AxY0x3lksq4hI0ZMHNYaKJE8BgOtCXjGFNNUAHE2yfMy1DgAqpZffGNMLwBZXUBERkYyw\n2eDlyFiNwZ5eAmPMEgCVkq4CQACjUkju+aSiKXPLb4xpDOBtAJ2yuF8RkaLJboeXkxnqfE43MJBM\n9aJsjDlljKlE8pQxpjKA0ykkCwdQPclyoGsdAJxMLb8xJhDAPAB9SR5Kq4z9+/dHUFAQAMDf3x/B\nwcEICQkBAISGhgKAlrWsZS0XzeXYWLR2OLF/8370/7I/ACReL1NjyMz/yDfGjAcQSXK8a7RRWZIj\nkqWxAdgLoAOAEwA2AOhDcndq+Y0x/gBCAYwh+WM6ZWBW3oOISKEWHw+nTzE89/MQTOo6KXG1MQYk\nTUpZstrHMB5AJ2NMwoV/nOuAVYwxvwAASQeAoQB+B7ATwBySu9PKD+AZALUBvG6M+csYs8UYUz6L\nZRURKXoSRiXFXfU4S5ZqDPmBagwiImlz2rzw+LcPY0avrxLX5WSNQURE8jnabIiP0yMxRETEhTYv\nxMd63pSkwCAiUth52RCnwCAiIglot8GhpiQREUlks8ERp6k9RUTEhXab+hhERCQJmw3xmsFNREQS\nGLtdTUkiIpKEza7OZxERucbYbHDGez5zgQKDiEhhZ/eGI1Y1BhERcTF2Oxzx6mMQEREXY/eGM05N\nSSIi4mLsdjA+Hk46PUqvwCAiUsgZux0+sCHeGe9RegUGEZHCzmaDD2xwOB0eJVdgEBEp7Gw2FFON\nQUREEtntKEYvBQYREXGx2eCtwCAiIonsdhSDDQ6qj0FERABXH4NqDCIiksBuV1OSiIgkYbOhGL00\nXFVERFw0XFVERNxouKqIiLix2WCnUWAQEREXV+ezhquKiIjF1fmsGoOIiFjsdtihpiQREUlgs8Hb\nqcAgIiIJbDZ4Q/cxiIhIArtdo5JERCQJNSWJiIgbV41Bw1VFRMSiGoOIiLjJzTufjTFljTG/G2P2\nGmN+M8b4pZKuizFmjzFmnzFmuKf5jTHVjTGXjDH/yko5RUSKNLsddiLXagwjACwlWR/AcgCvJk9g\njPECMAlAZwCNAfQxxjTwMP97AH7NYhlFRIo2V40ht4ar9gAw0/V6JoD7UkjTEkAYycMk4wDMceVL\nM78xpgeAfwDszGIZRUSKtlwerlqR5CkAIHkSQMUU0lQDcDTJ8jHXOgColCx/JQAwxpQG8AqAsQBM\nFssoIlK02WywOzxvSrKnl8AYswSuC3bCKgAEMCqF5PToqKlzuv4dDeADkleMMQnHTFX//v0RFBQE\nAPD390dwcDBCQkIAAKGhoQCgZS1rWctFd9lux44zl/H9619gVflVidfL1Bgy89dyY8xuACEkTxlj\nKgNYQbJhsjS3ARhDsotreQQAkhyfWn5jzEoAga5dlAXgAPA6yY9TKAOz8h5ERAq9jz5C6G+fYdfo\nwRjSYggAwBgDkin+6M5qU9ICAP1dr/sB+CmFNBsB1DHG1DDGFAPQ25Uv1fwk25CsRbIWgIkA/ptS\nUBAREQ/YbLk6Kmk8gE7GmL0AOgAYBwDGmCrGmF8AgKQDwFAAv8PqSJ5Dcnda+UVEJBvZ7bA5s7GP\nIS0kIwF0TGH9CQDdkiwvBlDf0/zJ0ozNShlFRIo8mw02J/R0VRERccnlG9xERCS/s9lgc1KBQURE\nXDLYx6DAICJS2LlqDHrstoiIWFydz6oxiIiIxW6Hl/oYREQkkTqfRUTEjd0OL93HICIiiWw22Byq\nMYiISAL1MYiIiBubDTanU8NVRUTExWZTjUFERJJQU5IUdQmzVokkKPLnhM0GL4dTgUGKriJ/EZDr\nFPlzwm6HceiRGPlGTp+QWd1/ZvJnJI8nadNKk9lt+ZnOifTT6pzI5v1v2aIaQ36ii0D6aXURyF/7\n1zmR+3L8nNi8OUN9DIZkjhYopxljCvYbEBHJIyRNSusLfGAQEZHspaYkERFxo8AgIiJuFBhERMSN\nAoOIiLhRYBARETf2vC5ATjDGlATwMYAYAH+Q/CaPiyR5zBhTE8BrAHxJ/l9el0fynjGmB4B7AJQB\nMI3kkjwuUr5RKIerGmMeBXCO5EJjzBySvfO6TJI/GGPmKjBIUsYYfwDvkByU12XJLwpEU5IxZqox\n5pQx5u9k67sYY/YYY/YZY4Yn2RQI4KjrtWcPB5ECJRPnhBRyWTgnRgGYnDulLBgKRGAAMB1A56Qr\njDFeACa51jcG0McY08C1+Sis4AAAKd7ZJwVeRs+JxGS5UzzJAxk+J4wx4wD8SnJrbhY0vysQgYHk\nagDnkq1uCSCM5GGScQDmAOjh2jYfQC9jzGQAP+deSSW3ZPScMMYEGGM+ARCsmkThlIlz4lkAHWBd\nK57M1cLmcwW587karjUXAcAxWCcBSF4B8HheFEryVFrnRCSAwXlRKMlTaZ0THwH4KC8Kld8ViBqD\niIjknoIcGMIBVE+yHOhaJ0WXzglJTudEJhSkwGDg3nG4EUAdY0wNY0wxAL0BLMiTkkle0Tkhyemc\nyAYFIjAYY74B8CeAesaYI8aYASQdAJ4F8DuAnQDmkNydl+WU3KNzQpLTOZF9CuUNbiIiknkFosYg\nIiK5R4FBRETcKDCIiIgbBQYREXGjwCAiIm4UGERExI0Cg4iIuFFgEBERN/8P1gCmgKErYhQAAAAA\nSUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "('z340', 'READ 2018-10-25/X133978.340')\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c2.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c2.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c2.dat')\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD/CAYAAAD12nFYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4FFX3wPHv3c0mhEASSqiBhI4iGkARFCSACNLt2EFU\nVLAXsOP7igJix/4TsWAHARudgIA0ARUhFF8IoRNKQvpm9/z+mABZSNn0ZDmf58lDZubenbvLZM6e\nuXfmGhFBKaWUOsFW3g1QSilVsWhgUEop5UEDg1JKKQ8aGJRSSnnQwKCUUsqDX3k3oLiMMTqsSiml\nikBETG7rfSJjEJEK+/P8889X6NcvSv3C1PGmbH5lirKttD/z8v4/02NCj4mSeP38+ERgqMiio6Mr\n9OsXpX5h6nhTNr8yRd1WkekxUXBZPSbK9/VNQZGjojPGSGV/D6pkjR07lrFjx5Z3M1QFosfEmYwx\niC9fSlIqp8r6rVGVHj0mCkczBqWUOgvllzFU+lFJSimVm8jISOLi4sq7GeUuIiKCnTt3FqqOZgxK\nKZ+U/Y24vJtR7vL6HLSPQSmllNc0MCillPKggUEppZQHDQxKKaU8aGBQSqkKYNiwYTz33HMsW7aM\nc8455+T6rVu30q5dO0JCQpg8eXKZtEUDg1JKVSBdunRh8+bNJ5cnTpxIjx49SExMZNSoUcTExNCj\nRw9CQ0Np2rRpqbRBA4NSSlVgcXFxtGnT5uRyUFAQw4cPZ9KkSaW2Tw0MSilVDtavX0+HDh0ICQlh\nyJAhpKenA7BkyRIaNWoEQM+ePVm8eDEjR44kODiY7du3c9FFF3HzzTfTpEmTUmubBgallCpjTqeT\nq666ittvv50jR45w3XXXMX369JPbjbHuO1u4cCFdu3blnXfeISkpiebNm5dJ+/SRGEqps5LJ9Z7f\nwivKzdUrV64kKyuLBx54AIBrrrmGiy66qGQaVAI0MCilzkrl+bSMvXv30rBhQ491ERER5dSaM5XI\npSRjTB9jTKwxZqsxZnQeZd4yxmwzxmwwxkR5W9cY86gxxm2MqVkSbVVKqfJWv3599uzZ47Fu165d\n5dSaMxU7MBhjbMBkoDfQBrjRGNP6tDJXAs1EpAUwAnjfm7rGmHCgF6CPSFRK+YzOnTvj5+fH22+/\nTVZWFjNmzGD16tUnt+f38D8RISMjg8zMTNxuNxkZGTidzhJtX0lkDB2BbSISJyJO4Gtg0GllBgGf\nAYjIKiDEGFPXi7qvA4+XQBuVUqrCcDgczJgxg08++YRatWrx3Xffcc0115zcbnJ0gJjTOkOWLl1K\nYGAg/fv3Jz4+nqpVq9K7d+8SbV9J9DE0BOJzLO/GOuEXVKZhfnWNMQOBeBH5+/QPRimlKrv27duz\nbt26XLflvKy0aNEij23dunXD7XaXatvKq/M53zO9MSYQeArrMlKBdYYOHUpkZCQAoaGhREVFnZzK\nLyYmBkCXdVmXz8JldUpMTAxTp04FOHm+zEuxJ+oxxnQCxopIn+zlMYCIyIQcZd4HFovIN9nLsUA3\noEludYGfgQVAKlZACAf2AB1F5OBp+9eJepRSZ9CJeizlNVHPGqC5MSbCGOMPDAFmn1ZmNnBbdmM6\nAcdE5EBedUVko4jUE5GmItIE6xJTu9ODglJKqZJX7EtJIuIyxowC5mEFmo9FZLMxZoS1WT4UkV+M\nMX2NMduBFGBYfnVz2w0FXH5SSilVMnTOZ6WUT9JLSRad81kppVSxaWBQSinlQQODUkopDxoYlFKq\nAtCpPZVSSuWqoKk9J02aRNu2bQkODqZZs2alMpObBgallKrATp/aE+Dzzz/n2LFj/Prrr0yePJlv\nv/22RPepgUEppcpBUaf2fOyxx4iKisJms9GyZUsGDRrE8uXLS7RtGhiUUqqMleTUnr/99tsZGUVx\n6QxuSqmzknmhZB6mIM8X/ia6kpra8/nnn0dEGDZsWKHr5kcDg1LqrFSUE3pJKYmpPSdPnswXX3zB\nsmXLcDgcJdk8vZSklFJlrbhTe06ZMoWJEyeyaNEi6tevX9LN08CglFJlrThTe06bNo2nn36a+fPn\nFzrL8JYGBqWUKmPFmdrz2Wef5ciRI1x00UVUr16d4OBg7rvvvhJtnz5dVSnlk/TpqhZ9uqpSSqli\n08CglFLKgwYGpZRSHjQwKKWU8qCBQSmllAcNDEoppTxoYFBKKeVBA4NSSikPGhiUUqoC0Kk9lVJK\n5aqgqT3feOMNmjVrRnBwMPXq1eOOO+4gOTm5RNuggUEppSqw06f2HDRoEGvXriUpKYnY2Fji4uIY\nN25cie5TA4NSSpWDok7t2aRJE2rUqAGAy+XCZrOV+KO3NTAopVQZK+7Unl999RUhISHUqVOHOnXq\nnJwJrqRoYFBKnZ2MKZmfIsg5tafdbi/01J433ngjiYmJbN26lU2bNvHGG28UqR150cCglDo7iZTM\nTxGUxNSeAM2aNWPMmDF89tlnRWpHXjQwKKVUGSvu1J45OZ1OqlatWhLNOkkDg1JKlbHiTO358ccf\nc+jQIQA2bdrE+PHjPWZ/KwkaGJRSqowVZ2rP5cuX07ZtW4KDg7n66qu5/fbbefjhh0u0fTq1p1LK\nJ+nUnhad2lMppVSxaWBQSinloUQCgzGmjzEm1hiz1RgzOo8ybxljthljNhhjogqqa4yZaIzZnF1+\nujEmuCTaqpRSKn/FDgzGGBswGegNtAFuNMa0Pq3MlUAzEWkBjADe96LuPKCNiEQB24Ani9tWpZRS\nBSuJjKEjsE1E4kTECXwNDDqtzCDgMwARWQWEGGPq5ldXRBaIiDu7/kogvATaqpRSqgAlERgaAvE5\nlndnr/OmjDd1Ae4Afi12S5VSShXIr5z26/UDRowxTwNOEfkyrzJDhw4lMjISgNDQUKKiooiOjgYg\nJiYGQJd1WZfPsuWIiIgz7gE4G5148mpMTAxTp04FOHm+zEux72MwxnQCxopIn+zlMYCIyIQcZd4H\nFovIN9nLsUA3oEl+dY0xQ4G7gB4ikpHH/vU+BqWUKqTSvo9hDdDcGBNhjPEHhgCzTyszG7gtuzGd\ngGMiciC/usaYPsDjwMC8goJSSqmSV+xLSSLiMsaMwhpFZAM+FpHNxpgR1mb5UER+Mcb0NcZsB1KA\nYfnVzX7ptwF/YH52OrhSRO4rbnuVUkrlTx+JoZRSZyF9JIZSSimvaWBQSinlQQODUkopDxoYlFJK\nedDAoJRSyoMGBqWUUh40MCillPKggUEppZQHDQxKKaU8aGBQSinlQQODUkopDxoYlFJKedDAoJRS\nyoMGBqWUUh40MCillPKggUEppZQHDQxKKaU8aGBQSinlQQODUkopDxoYlFJKedDAoJRSyoMGBqWU\nUh40MCillPKggUEppZQHDQxKKaU8aGBQSinlQQODUqricrvLuwVnJQ0MSqmK6ddfoUULcLmK9zqT\nJ0Pr1jBwIKSklEzbfJwGBqVUxZOQAHfeaZ3Ily8v3mvNmQP33QdHjsC8eSXTPh9nRKS821Asxhip\n7O9BKXWaoUMhLAyqVYNjx+D114v+WhERsHChlYH88QdMnVpSrazUjDGIiMl1W2U/qWpgUMrHZGRA\nnTqwfTvs22ddAtqxA0yu57D8HTsG4eGQlATx8XDhhdZr+vmVfLsrmfwCg15KUkpVLEuXQps2VsbQ\ntq11Et+woWivtXGj9Vo2m5U5hIfDihUl214f5BOBIcuZUd5NUKr8JSZCBcmenS4nE359hknRAWxo\nHMDW8xuS+sFk70YZ/fwz9Otn/W4M9O8Pc+cWrSF//QXnn39qefBgmDXLu7pZWbB+ffE7vyshnwgM\n25uGsu6r18q7Gb7t0CEYMMD6I+vdG5YtK+8WKbCCwbhxSJs2uOrWIb1BHZIeGVmuJzO3uLnx4ysZ\n8MBkRoT1IezdqSwacB5bXnqEuKt6WCfcvIjATz+dCgwATZvC7t1Fa8zff1tZxwn9+1uBpyD790Ov\nXtC3LzRuDL/8UrT9V1I+ERiOPTKSOiOfYEWncHb9ubS8m+N74uLg0kvZE9qG5fd+wYGeN8Ett8Dw\n4ZCWdkZxt7g5mHIQp8tZDo0tP5muTJ5f/Dznv3c+I34cweZDm0txZ5nw2mtI8+bsXDWXIX1T6fBq\nS+5+uDl/zv6Inf27gLN8Pv/3V7/Lk6+tpvUVN1P925k07H8j94ybS+rcn4nfupZFlzdj7/G9uVfe\nuhXS0+GCC06tq1/f6hcoir//9swY2rWDo0etPou8pKZCt27QpYsVkD76CO6/v8if57ZtsGlThUnm\nvCMilfoHkEsuEXlo5GH5ql9PORxoZNFNl8jxw/tElQC3W7Kie8isi/4rjSKypNOgDRLSaYb0um6u\n7L/iKpHzzxdZvlwOJB+QD9Z+IL0+6yUB/w2QkJdrSP1JDeTFJS+Ky+0q73dR6tKd6XLxRxdL3y/6\nyY9/LZVxS8dJ7Ym15bt/vivZHbndIjNmiDRrJomXXya3vXSRnP/e+TJv+zxxu93idous375alrYO\nlD9v7V2y+/bCjqM75JGrgiT1wgsk4UCWrFkj4nSe2p6edFQONq4td95UXcYuHisHkw96vsB//iMy\ncqTnut9+E+ncufCNcbtFgoNFEhI81996q8i770pmpshnn4n8978iCxbk2P7AAyI33+xZp1cvkfff\nL9Tuk5NF7rxTJCxMpFEjkXPOEVm3rvBvo7RYp/88zqt5bSjMD9AHiAW2AqPzKPMWsA3YAEQVVBeo\nAcwDtgBzgZA8XlcWLxaZMEHkmmtEOjRbLdMiIiW+uk1eHTRCfpjhlL17S/HTLYSUFOvAiIkRSUsr\n79Z46YsvZEeNKOlxwwoJf7WxtHq7lfT9fIBE/rezmNEhMvSqYNkdYpNlkXb5qmtL+bD99fJ80Fh5\nLmiSXFXnHWn0YCe5f/bjZdZct9stccfi5Le43yRmR4xs2LdBEo6lyYoVIu+9J/LEEyLDh4vcd5+1\nHB9fMvt9asHTct6Lg+SmsLky336F7K7SVPa0aS8v9wuRjxe/XvwduFwiM2ZIVof2crx1U3nthSsl\nbGKYTF41WTKznDJjhkiHDiKBgSI1a4qMGrVE9lc3cnDOjOLvuxDu/XCQHA+uKnd1/luCg62TYViY\nyDPPWCdKERFZvVqcYbXksU9uktDxoXL37LtlS8IW60TerJnIqlWeL7ptm0hkZOEbs2uXSL16Z67/\n8kuRAQPk6aetz2z0aJHGja2T+J4fpktKnRrywLRbpeuUrtLy7ZYS9X6UPPlST0muU0PSU497tWu3\nW2TIEJHrrhM5ckTEnZ4h62+eKAODF8uvv7gL/168lZkpsnevSFZWgUXzCwzFHq5qjLFln9R7AnuB\nNcAQEYnNUeZKYJSI9DPGXAy8KSKd8qtrjJkAHBaRicaY0UANERmTy/7l9Pewfz8s/OBD2nz4ME4x\njPZ/k61Zw2nWzBroULt27v+GhUHDhkUbFQewb9Nq9q1dTObRBEJaR9G060ACqlbn6FGYOBE++MAa\nFFGlCsTGWpc7H30UOnQo2v7ydPgw/Pab9UE4HNC+/anRHYWRlERaZGuua/Ioq26YwEcDP2Rw68En\nNx86JPz3zb2sWpRGq7hYetaMpW1EEi0jM6nml45z7QYy1/7Jgz0CaP/MOO67ZHjx39uiRfDpp1a6\nf845cPXVJLRoyI+x85j513yW7Z2PPd1F939r0mXrcS7ce5h6qekk+1VlS3hXDkWPx6/NBaSnw+rV\n1tD2O++E556DqlWL1qQ/9v5Btw+vZNJHN3NH1nfYXnuVBQlRfDp2B6Obf0CDLT+y7p5B9H7pG0wB\n/wcZWRls2L+BVXtWsePoDhIS93Huwr+47sd/SfTLYuJlfmzvci7XnHcdd7a/k6R9dbjrLuverZde\nsq6AHDwI48aBa8sQXtr2I/V3HMIU9c0Vwrp969g6oAvHdt7HjuGTGDfOOuT+/Reefda6T+211+Dq\nq8G8MBZWr+bAN1N4Z+27vLf2Pe5Jb8tTn+8kcMu/nn+EKSnWH2lqKhjDsWPwzTcQ2TyN3TW+4rtN\n3/LPoX9ITE8k0BFI/Wr1aV+/PX2O1abvpNnY1q2jqiPH+z98mKzGTWhW/RALVqezOW0Jv8YuZPry\n+SydsoVvb7qQGjffQtu6bakbVJe0rDT+OvAXXQY/yNNXBXPPA5/RvUn3fD+LN96Azz+3uuICqwju\noUPJ2hqL80AKG3aHcWzaz/S7rgT/T0RgyhTkxReR40mYlFTMu+/CsGF5VinV+xiMMZ2A50Xkyuzl\nMViRaEKOMu8Di0Xkm+zlzUA00CSvusaYWKCbiBwwxtQDYkSkdS77PyMwnOB2ZbFi3L00e3UKsedF\nknDlQzjrD+RYRgQJCVZ/as5/9+8Hf3/rhD1ggPVHVtDfU0ZyIqtfHkWtaT9QJyGVXY1DcAZVodae\nowSlOJkX3YlxmybTLbodzzxjjZgD6w95yhR44003YR2W0bjnLxwK+J0kZwJ10mxcmBZBeEojqh9v\nTIB/FO7a5xLStjEdLzbUrZtLQw4fhunT4csvrZEUnTtDZKTVB7BmjfVH9cwz1oFit+f/prKlPPoc\nP330J/eMWcGc237i4vCLvaqXk/y9kd3dBjO72R4Gzd9GeGh4oV8DAKcTGTYM92/L2X/L4/x5LIuU\n9T9w4aaVQAZLaoWTbFrT+XgqbY7+yZGG55N0cS/8r4gmrGMd1mycRfw3H9J3wS72d2hF7XGvU//S\n3hw4AA88YF0HnjPHGj5fqPcnQovxl3DDp00Ya9uIY8lC6xsGsHMn9OkDI3v9TJe5QwjNsGFee53I\na+44WTc+KZ7f439n5e6VrNyzkr8O/EXLWi3pUrsDg1cc4eKvfyMzIpzEh+8jtN81hAbWwBhDVpZ1\nkp04EZ56Ch580PO/VQRGP51Bp+/DiLisNx3+77sifezxifFkuDIIDw6nil+VfMve859O/OelzXw+\nOp5HXwg+Y3tMDIwaBTVrwkMjnQwYfwl+dw/H3HsPyZnJ7Lj2cma7N7F0SGcev+RxejbpiTkRIIKD\nce/cxasfhzJuYioRV3/AphqvUD25A+/cfSuXRnQiOCCY9Kx0diXuYv2+9aT/NJOLv11Oj5uzqBNU\n5+SPn82PMY8s4L1oB9+fl0nn8M70bNKTW2bHkfzjXrodmc1jj8Ftt0G9ejnewKOPEus+SI/whTza\n+VEe6fzIqfbl8M8/EB0NK1Y5+f34VxyeOJauS3Zy1T2hJLhTmDYrgKoHAnjl3mu5olsTqvkHEegI\npIpfFQLsgaQk1CB5T2OahDahdWtDkyb5/x+lHDnAwRv6kfHvVu7r5WRtRACRB9JYOs2B/ZVXqT5s\nRK71SjswXAP0FpG7s5dvATqKyAM5yvwIvCwiK7KX5wOjsQJDrnWNMUdFpEaO1zgiIjVz2X+egeGE\n4wl7+eORIdRavp7I3ckkBtnZ06Q2aee2IKB9R+pdcgWN23fH5udPbCzMng0//gjr1kGrVtCxo/Wl\n+5xzrJ+wMEg6tJv1z4+g1bQ57GpSC9uDD3HBTQ/jCAgErG9tL97xM503jqHPwX/4/bIm2B9+lBYd\n+xDoF8jWw1v5ZdsvfP3Xl7TZUYVuayOJ+ieRdqnb8cPF7qohJAYLGVVSqJqZSGSioXq6YWNQPeKC\nWmGLbE7TVg04r57g//sS+GMtKd27sG9Ad3Z0akWC6zipzlSy3Fm4xEXz2INc9u4vBFaphvniCyt1\nyc/+/SRHnssld9fi/tuf4K4OdxV4LOQlMyGJLS1asrRjbUbO3ehVnb174fvvrRtW9+x2MWb7tQQS\nx01dO5PRYhkE7edcv4H0Ch/I7SF1aHn4b/wS9luZUXQ0hITk+rqxO9YSO+5hOn27nLjw6vzTpz27\nWtVj1i4nsYm7aNHaTaOadenfoj83n38zwQFnnuByemP+d/w6+Ul+WpGCY83vVjDOYc8euOwyGPVg\nOvWP3kbnt2bwb70Apg1uyq81D+MSF53DO9MpvBOdwzvTIaQ11T7/BsaPh4sugiefhE6dTr6e02nF\n/3HjrCD20UfWoJ3ciMBjD87kyf+7mmNzV9O864VeffYAac40Hp77MN9t+o7ggGAC7AF8f/33nFfn\nvFzLz98+D/9eg9lSYyJ3/TEqz6z7RPvffx/S1scyO6kbU2s9RrUaDm6Lf5EvRq8m7pwlzDo0iQCH\ng/s73k+3yG7Ub9eHoU0nsCpiPennfkzXyEt4puuzPHpzFDfcACNyO/dNmwY//UTmF5+y9/heDiQf\n4GDKQXbtzSRmWDxf2t/EbNyIX2CQlc4MHgxr17I5NYIXX7QGIoWFWbdBtGsHV/v/RJv5rxM/4xMG\nfz2YNnXa8GH/Dwl0BJ7cpcsFl1wCXW5Zwjy/UdQKrMXMV+IJfPt9Anr0Ij0rnR2HthEw5B52xyZw\nV+8rCG3uwmXSSDiWSp2tB+hzcCfnphxkR4iNeVW7sCvoVu6I7sVdN9Xx+OLy5/4/+W72K9z47Fds\nrNeQx6q9RtK2K2hQK5gDxxM4/6L7+WHut6Rv/Yf64Wd8p66QgWEB8ASFCwyHRaRWLvuX22+/ncjs\nP8jQ0FCioqKIjo4GICYmBuDk8qKFCzm4bQMR6fvJWLeKTas2UntfIn0z3OxsGETM5b1oM+R+uvfo\nQXo6TJkSQ2wspKREs3kzxMX+H72qvssrRzawMqIRv7YdSuSF3enfPxqbDWbNimHNGli0KJrhw+Hy\ny2NwH95F/flf0+T7BcyoIvxb09ChQU3OTanK7s0HsdeoSfTAgUj3Hix0upCatel1RfeT7XeLm8YX\nNOaf2N9Y8vmXpG35H20yU5CMJP5KT2NjmCH2wupIldrY4xz4SzA1IprjkCBS4g6A2KnSwsFB2wqG\n/JpK302Z9F6yHHPBBWd8PieW27zxBdO2LuXXG1rwZLfH8/w8vV2u6W5M9ataMbPfINrdPSrP8gsX\nxjBjBnw9vSPnXTuTA6536f77Wu446uaVF66k5vH6tK3blnuuvQe7zV7k9lzc4QLipr7JHx9/Qcju\nQ1yZmAJumG6CcTWvxfFLg/iwVTIjWzxNkxpNcn29xJR0Wt4UwccL0uj/40zo0SPX/e3bB6NHRzN2\nLEQ0+JXEGV/S95f5mJq1WNE2CtOgAdEREbBmDTHTp0NUFNFvvw3t2hETE4MIBAZGM20afP55DI0b\nw/jx0Vx5JSxZkv/7Xbw4hpmP3c/gfQmc+8ceNm9ZWuDnk56VzgtxLxAREsGtwbcS5B9EXGgcj81/\njKcbPU1UPc+/L7e4+Xzq3Tz1fTpbv/qUoGC71/8fv079GjN+PB0DQpl//f/x9drdbN0KO+O6Edxu\nDum1x3O8yiYWrDjK2z0bYGvdmb4t+jDsKusSyYcfxjBmDMTFRVO9+mmv/9ZbxCxeDA8+6LH/r74C\nY6J5f3d/YkJDrc/71Vfhk0+IqVLlZH2XCz79NIa4OEhLi2bxrCT+u70u8W/P4ua7ujB89nDWLF/D\nAx0f4IEh1unumjs+Y/GRT6nWbRuv9X6NWv8I5o47iD56FPz8TrWvSxdc94xk7oyf2F67Mx2q1uWc\nQ0tYn5mAPbor0X37cnTjGpZ/9iW7Ag2f9MlinaMZ1XaHUL26jYxmCXTasp9h05OZU/cmmo36PwYM\nNBw4EIMx0LVrNDNnwt8jmhHfLIm3l8WxevlqpmY/DiQyMpIXXnghz8BQEh3PnYA5OZbHcFoHNPA+\ncEOO5Vigbn51gc1A3ezf6wGb89h/gZ0s3ji2P07WTHlRtjQIkHXnhMqycffIzjULJCEuVrYv+1GW\nPjdU/ji3hhytYmRxtxYy4/05MnmyNYCiRw+RVq1EWrQQ6dlT5KmnRHbvzmUn6ekiv/8u8tFHIu+8\nI/LDD1LcnvG0dJfExorMmiUyaZLIiy+KjBsn8vLLIuPHW53yEyeKPPmkSL/+Lql6/hy5Y3CkHAn2\nl6Mxc3J9TeeipbI3METaPhslac6S6yV/4bH/k/1VbZI8a1au21NTRa66SqT14BlSb2K49P68tyx6\n/UFxNqwv7v37S6wdecrIkOmv/CvDQ76ThCuGSErtUBl+c3VZsWtFrsUvHf2ULG0UJu5nnyvwpbds\nEWnQQOTNN62OScnKElmyxOqVvf12kREjrGMix4Gzd6/I889b/bEtW1oDdrZvL/zbSk45JqsbBsjT\nLe6V2Nj8y2Y4nXLda73kleFdJP7+8XL4mVfFHbtFREQW/m+hhE0Mk7V71nrU+c+Mj+WfGlVk5+TZ\nhW9cHrKyRDZtEvnlF2vARta1N4hMm5Zr2Ztvto71Mzz3nPUBnqZjR5H580Vk61aRiy+2fqZO9apN\nB5p0lIEhMfL66yIul1um/TVNGr/eWCJej5CG41uI7fH68tCsZyU5I7un/YMPRG68MfcXdLtFli+3\n/khffVVkxQprkEFOmZki77wj7nr15MCAK+SzRx6X5+4bK3Muu11Sa9SXIzOX5Nvm/32/RnYFVZGH\nP3nnjG2U5qgkwA5sByIAf6xRR+ecVqYv8LOcCiQrC6oLTMgRJEYD4/PYf74fTGFlpqXIb+NGyMr2\ndSWulp8cCTTyvzoOWXFxQ1nx2iOSmni4RPdX1hISREY+kCHXdr5KDgXZZNv0jzwLpKbKntqN5dpr\nQmTn0bgS3bfLJTJg8BVypGqgyB9/eGw7fFjkki5OafHgfdL8zRYSsyPGOpvWrm0F0zI0a5Y1kub3\n11ZIcpOG8nHnKrLq3988yrw/Y4N8cEEVOd71Ms/xmPnYsUOkTRtrVFRiYt7ltm0TuftukdBQa/TU\nmjXZwaQYls+fIglVbdK79q/y3HMif/1lfeb/+5/I99+LjH7cJU+c/7PMrRMuR/z9ZEHIIPm0zmMy\nNXCEHLTVlX+a9pP//bZbZm6eKWETw2Tm5pkiIjJl0WJ5sGeQ7DznwuI3Mj8PPWR988nFggUiF16Y\ny4b77hN56y2PVTt3WodUZmYR2zF6tBy5/znp2FHksstE1q4VSUnPkNen/is1Wm+Q+QtPGw00YECe\nAa1QkpJEXnvN+hYaFSVy//3WqCsvrG0ULi/edOsZ60s1MMipIadbsIajjsleNwK4O0eZydlB4E+g\nfX51s9exHmRWAAATO0lEQVTXBBZkb5sHhOax78J8vCrbnDki/ds8KwerGom5s7dkpqWI7Nsnu5uf\nI1POC5SvVi0o+EWKYMPWg3LNwDA5Elzj5NffzZtFml+wXyKevVx6f95HjqUdEzl+3DqLvvdeqbSj\nIEuXitSvL/LGC8dkX/eO8mcDu8ycOkayXFkybeYKeaNdsOyObGz9wRZCYqIVGBo1shKEI0es9U6n\nyMqVIjfcIFKrlpVIHDhQsu/p/f8OliPVHPJuv+nSrKlbatQQaV4/WSaf964cqt1KdjSuL0/e0kgS\nE0/dW+B2i/y9LlPmXvqCHDRhMum8T+TZ91ZK2EuNJPj5cOl+U5Acrx5qffsuTRMmiDz6aK6bMjOt\nIbpnnCdvODPLmDzZStCK7OefRS6/XLKyrNeKiBBxOKzbeVauPK1saqpI9epWBC5H3190kbzV58x7\nWko9MJTnjwaGotuxQ6TvpUtkbuPaIiBOg4ztHCrjpn9fqvudu26jjOhdXfYHBsudPb6WoOh3JeQ/\nDeSpBU+J0+W0BrxfeaXI0KGl+y20AHFxIl27irSLcsus4RNkf7BDdgUb2RdkZH3nS0SKcXnrt9+s\n81ZgoPUNtlo163LkK68UOtZ4zelyymPPdJT4RiHiioiwBvGHhop78GB5a+K10vKtFrInaU+e9dNX\nbZDDjaMktval8maLV+XnC++VjNBaIgsXlk6Dc/rsM5Gbbspz8223nZEciFx+ucjcuR6rBg2ybmMo\nsv37rSiU47hMScnjtoHFi63LVOVs+iVd5I0e3c5Yn19g0Mdun+VErFEic2MO4QyK4+kR59OiqX+p\n73fZv+uZ/8ojPPD5UpZc0Zr6z75M53YDrIHfDz9sjS764ANr/HA5EoFvv4WZM2HXTjc9ztlA394p\ndL6ha4m8vtttjWCrUgVCQ0vkJfN1POM4I3++j8OrYhgYcQXJDcOYsudH6lWrx3fXfUfNwDMG/nly\nOq1nGf34IzRoANdeC1FRpd/wBQvg5ZetYWq5mDkT3nrLus3lpPbtrWFb2TcKZWVZt0Ns3Vr4Ycke\nwsOt+4QKGkc6aRLs2mU1rBzN7N6Lne5kHlryu8d6nY9BVVzx8db9FT/8YJ10mjaFkSPh3nuLfqeh\nKtD8f+ezPH45SRlJDG49mC6Nu2AzFfjRaRs3wvXXWw8dykVqqvVIpR07rPskAOvhd0uXnhxC/Pvv\n1mFV1Cd4nzRwoHWTw7XX5l/uhhusm6JuvbWYOyyeH/sM4N/je3ho+TqP9fkFBp2tQpWvRo2sO5mn\nToXjx6F6dQ0IZaBXs170atarvJvhvQIepFe1qnVD6ty5cOON2SsPH7ZShGzz51sPTC22Dh2sm5wK\nCgxr1sDYsSWww+KxB1TBHC3c1AQV+CuCOqsYA8HBGhRU7mrWtNKCXJ7me0K/fjmeqJ2ebl07Cgo6\nuX3+fLj88hJoS4cO1hSh+UlIsH5atSqBHRaPo0oVbC4NDEopX2MM1KplPUsmD337Wo81cbmwsoVa\ntU5+0Th2DP78E7qWRNfQicCQ3yXsP/6wytnK/xTrHxiIzZVZqDrl32qllPJGaKh1hs9Do0ZWf/iq\nVVjf1mudelDC3LlWUCiRZwnWr28Niti1K+8ya9ZYjzSpAAKqBmJzF24uCQ0MSqnKISTEmrEuHycv\nJ53Wv/DTT9aDMUtMQZeTKlBgqFItCD8NDEopn1RAxgDWI72/+QYk4fDJjCEry3q8es7ZQoutfXur\nAzova9bAhd4/tLA0BQZVxS4aGJRSvsiLjOHCCyEgALavPhUYVq60bj1o1KgE25JfxrB3rzX16mlP\n2i0vQdWr4qeBQSnlk0JDCwwMxlhTjvy96FQfw1dflfBlJMi/A/rEZaQKMsKuWnAQDsnKt6/8dBoY\nlFKVQ0hIgZeSAG65BfZvOkxG9dps3WpdWrr//hJuS4MG1oij3bvP3FaB+hcAAoOq4S8uMgoxYlUD\ng1KqcvDiUhJYs65dEH6YiR/X4t574YknivkIjNwYk/flpAoWGPyrVMXhdpN03OV1HQ0MSqnKwYvO\n5xMuaXmY87vXIi3Nmrq1VOTWAS0Ca9dWmI5nABMQgCPLxuGkvG8OPJ0GBqVU5eBlxgBgEo8xaGgN\nVqywHlBYKnLLGHbsgMBA616HisLhwN9lJyEx1esqGhiUUpVDITIGEhOtR6yUposvtp7M58wx4mfl\nygp1GQkAhwOHy8aR4xoYlFK+phAZA4mJVvnSVL8+NGtmPSr+hFmzrGdzVCQOB/5uDQxKKV/kxXDV\nk5KSSj9jABg0yAoGYD3gb+5cGDy49PdbGA4H/m7DsRQNDEopX+PlcFXcbkhOth7hXtoGDrQCgwjM\nmwft2kFYWOnvtzAcDvxdkJiqgUEp5Wu8vZSUnGw9bttuL/02tW1r/fvXX9ZUiNdcU/r7LKzsjKEw\ngUEn6lFKVQ7VqllzMmRlgV8+p66y6Hg+wRgYPtwanmqzWdOPVjT+/vi7ISlNMwallK+x2awTflJS\n/uWSkkq/4zmnZ56xJgZKSICGDctuv95yOHC44Xi6BgallC/yZshqWWYMJ9jtZdOnURQOBw63kJyh\ngUEp5Yu86Wco64yhossODCmZGhiUUr7Im5FJ5ZExVGQOBw6XWwODUspHeXMvQ1nc3FaZOBz4ud2k\nZWlgUEr5Im8vJWnGcIq/P34uDQxKKV/lbeezZgynOBz4uVxkuFK8rqKBQSlVeWjGUHg2G26bIUsD\ng1LKJ2nGUCRuux+igUEp5ZM0YygSt58fSLLX5TUwKKUqD2+Hq2rG4Mnhh82tnc9KKV/k7XBVzRg8\nOfywoZeSlFK+SO98LpoAf/xMKi6Xd8U1MCilKg/tfC4S4/DH355GZqZ35TUwKKUqD+18LhLjcOCw\nlVFgMMbUMMbMM8ZsMcbMNcbkGqaNMX2MMbHGmK3GmNEF1TfGXG6MWWuM+dMYs8YY07047VRK+YgT\nnc8iuW/PzLTmawgMLNt2VXDGP6DsAgMwBlggIq2ARcCTZzTIGBswGegNtAFuNMa0LqD+IaC/iFwA\nDAU+L2Y7lVK+oEoVa16G9PTct5/IFowp23ZVcDZ/f/xMOhkZeQTU08sXc3+DgE+zf/8UyG0W7I7A\nNhGJExEn8HV2vTzri8ifIrI/+/d/gCrGGEcx26qU8gX5DVnV/oVcGYc//i5IzcjyqnxxA0MdETkA\nkH0ir5NLmYZAfI7l3dnrAOoWVN8Ycy2wLjuoKKXOdvkNWdWhqrnz98c/y4/UdO9OowXO+WyMmQ/U\nzbkKEOCZXIp7l6fkzaO+MaYN8DLQq5ivq5TyFfl1QOtQ1dw5HPi7/EhOywSqFli8wMAgInmelI0x\nB4wxdUXkgDGmHnAwl2J7gMY5lsOz1wHsz6u+MSYcmAHcKiI782vj0KFDiYyMBCA0NJSoqCiio6MB\niImJAdBlXdZlX1l2u4nOvpR0xvZlyyAzE2upgrS3Iiw7HBzfC/959h4a1K1y8nyZFyN59e57wRgz\nATgiIhOyRxvVEJExp5WxA1uAnsA+YDVwo4hszqu+MSYUiAHGisjMAtogxXkPSqlK5rrrrJ/rrz9z\n22efwbx58MUXZd+uiqxfPwZX/Z1h9/7NoB7WlXxjDCKSay99cfsYJgC9jDEnTvzjs3dY3xjzE4CI\nuIBRwDzgH+BrEdmcX31gJNAMeM4Ys94Ys84YU7uYbVVK+YL8Op+PHbP6IJSn7EtJaZkl1MeQHxE5\nAlyey/p9QP8cy3OAVoWoPw4YV5y2KaV8VEGdzxoYzuRwEJBpIzXDuxsZ9M5npVTlohlD4TkcBLjt\nXmcMGhiUUpVLfhmDBobc+fsT4LKT5uWtzxoYlFKVi2YMhedw4O+2ka4Zg1LKJ+V3H4MGhtw5HAS4\nbZoxKKV8lF5KKrwTGYNTMwallC/SZyUV3snAoBmDUsoXacZQeA4HVdyGDM0YlFI+Ka+MQcRarxnD\nmfz9CXAbMrI0Y1BK+aLgYEhOBrfbc31KCgQEgL9/+bSrInM48BcbGVmaMSilfJHdDkFB1pNUc9LL\nSHlzOPB3oxmDUsqH5TZkVR+HkTeHgwAxZGrGoJTyWbl1QGvGkDfNGJRSPi+3DmjteM6bw4G/GJwu\nzRiUUr6qZk04csRznWYMefP3J8AtZLo1Y1BK+aratSEhwXOdBoa8ORz4C5oxKKV8WFiYBobCyO5j\ncGrGoJTyWbVrw6FDnut0VFLeHA4cbsHp1oxBKeWrNGMoHIcDh4hmDEopH5ZbxqCBIW8OB/5uIUsz\nBqWUz8orY9DhqrlzOPBzC1miGYNSylfpqKTC8ffH3+0mSzRjUEr5Kr2UVDgOB34ut2YMSikfFhwM\nGRmQnn5qnV5KypvDgZ/bjUszBqWUzzLG83KS02kFhtq1y7ddFdWJwIBmDEopX5azA3rfPqhTB/z8\nyrdNFZXDgZ/LhRvNGJRSvixnxrB7N4SHl297KjKHA7vLjctoxqCU8mVhYac6oHfvhkaNyrc9FZm/\nP/asLM0YlFI+LmfGEB+vGUN+HA7sLhduzRiUUj4tZx+DXkrKn8OBLcuF22jGoJTyZTnvZdDAkD+H\nA5srCzGZiBRcXAODUqpy0ozBew4HNqcL7E5croKLa2BQSlVOOTOG+HjtfM6Pw4HJysL4ZZLpRTeD\nBgalVOV0IjBkZcHBg1C/fnm3qOI6ERhsGhiUUr6saVPrEtLff1tBwuEo7xZVXMYgdjt+GhiUUj6t\nWjUYMAAmTtT+BW/Y7diNUwODUsrHDR0K33yj/QvesNuxlUXGYIypYYyZZ4zZYoyZa4zJ9dGGxpg+\nxphYY8xWY8xob+sbYxobY44bYx4pTjuVUj6qe3crW9CMoWA2G3acZGR4UbSYuxoDLBCRVsAi4MnT\nCxhjbMBkoDfQBrjRGNPay/qvAr8Us41KKV9ls8Ho0XDppeXdkoqvEBlDcR9FOAjolv37p0AM1sk+\np47ANhGJAzDGfJ1dLza/+saYQcD/gJRitlEp5ctGjizvFlQOdjs2KZs+hjoicgBARPYDdXIp0xCI\nz7G8O3sdQN3T6tcFMMZUA54AXgBMMduolFLKbseOm/QMd4FFC8wYjDHzyT5hn1gFCPBMLsW9uNk6\nXyda/DzwuoikGmNO7DNPQ4cOJTIyEoDQ0FCioqKIjo4GICYmBkCXdVmXdfmsXl7icnH8F8N//ryd\nSy9pRn6MePPgjLwqG7MZiBaRA8aYesBiETnntDKdgLEi0id7eQwgIjIhr/rGmKXAid6kGoALeE5E\n3s2lDVKc96CUUmeFBg1ofP0x3rziAFf1rY4xBhHJ9Ut3cS8lzQaGZv9+OzArlzJrgObGmAhjjD8w\nJLtenvVF5DIRaSoiTYE3gJdyCwpKKaW8ZLfjcPuTllHwE1aLGxgmAL2MMVuAnsB4AGNMfWPMTwAi\n4gJGAfOAf4CvRWRzfvWVUkqVMLsdh9uP1IyCe5+LNSpJRI4Al+eyfh/QP8fyHKCVt/VPK/NCcdqo\nlFIKsNvxcztIzyz9jEEppVRlYLfjL3avMgYNDEopdTaw2/ETB+lOzRiUUkoB2Gw4xI80L+5w08Cg\nlFJnA7sdf/HzKmMo7iMxlFJKVQZ2O/4I6c5SHpWklFKqkrDb8RfI0IxBKaUUYN3HAGRkacaglFIK\nwGbDH0NGlo5KUkopBdl9DHbtY1BKKZXNbifAQKYXGYMGBqWUOhtkj0rKdGnGoJRSCqzOZwOZLs0Y\nlFJKAdhsBGC8yhi081n5nBOzVil1gh4TZPcx2HF6kTFoYFA+R08C6nR6TGD1MRgbTrdmDOWutA/I\n4r5+UeoXpo43ZfMrU9RtFZkeEwWX1WOiFF7/RMbg1oyh3OlJoOCyehKoWK+vx0TZK7vA4F3GYESk\nVBtU2owxlfsNKKVUORERk9v6Sh8YlFJKlSy9lKSUUsqDBgallFIeNDAopZTyoIFBKaWUBw0MSiml\nPPjks5KMMVWBd4EMYImIfFnOTVLlzBjTBHgaCBaR68u7Par8GWMGAf2A6sAUEZlfzk2qMHxyuKox\n5hbgqIj8bIz5WkSGlHebVMVgjPlWA4PKyRgTCrwiIneVd1sqikpxKckY87Ex5oAx5q/T1vcxxsQa\nY7YaY0bn2BQOxGf/7iqzhqoyU4RjQvm4YhwTzwDvlE0rK4dKERiAT4DeOVcYY2zA5Oz1bYAbjTGt\nszfHYwUHgFzv7FOVXmGPiZPFyqZ5qhwU+pgwxowHfhGRDWXZ0IquUgQGEVkGHD1tdUdgm4jEiYgT\n+BoYlL3tB+BaY8w7wI9l11JVVgp7TBhjahpj3gOiNJPwTUU4Ju4HemKdK+4u08ZWcJW587khpy4X\nAezGOggQkVTgjvJolCpX+R0TR4B7y6NRqlzld0y8DbxdHo2q6CpFxqCUUqrsVObAsAdonGM5PHud\nOnvpMaFOp8dEEVSmwGDw7DhcAzQ3xkQYY/yBIcDscmmZKi96TKjT6TFRAipFYDDGfAmsAFoaY3YZ\nY4aJiAu4H5gH/AN8LSKby7OdquzoMaFOp8dEyfHJG9yUUkoVXaXIGJRSSpUdDQxKKaU8aGBQSinl\nQQODUkopDxoYlFJKedDAoJRSyoMGBqWUUh40MCillPLw/x0nsA4UimbHAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "('z340', 'READ 2018-10-25/X133928.340')\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c4.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c4.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c4.dat')\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD/CAYAAAD12nFYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPmVQSSEIPUpIAUkQQBAQUJYIK2FDWhg3U\ntax19+euYEXdxcLaVlHUFQXFvoIdpUYFlCIICqFJ770kIW3m+f0xARNIzyQh4ft+vfIi995z5j4z\n3Nxnzjn33uPMDBERkUM8lR2AiIgcW5QYREQkDyUGERHJQ4lBRETyUGIQEZE8gis7gLJyzumyKhGR\nUjAzl9/6atFiMLNj9mf48OHH9OuXpn5J6hSnbGFlSrOtvD/zyv4/0zGhYyIQr1+YapEYjmWJiYnH\n9OuXpn5J6hSnbGFlSrvtWKZjouiyOiYq9/VdUZnjWOecs6r+HiSwHn30UR599NHKDkOOITomjuac\nw6pzV5JIblX1W6OUHx0TJaMWg4jIcaiwFkOVvypJRCQ/8fHxrFu3rrLDqHRxcXGsXbu2RHXUYhCR\nainnG3Flh1HpCvocNMYgIiLFpsQgIiJ5KDGIiEgeSgwiIpKHEoOIyDHghhtu4JFHHmHmzJm0bdv2\n8PoVK1bQqVMnoqOjGTVqVIXEosQgInIM6dmzJ8nJyYeXR44cSe/evdm3bx933nknSUlJ9O7dm5iY\nGJo3b14uMSgxiIgcw9atW0e7du0OL0dGRnLTTTfxzDPPlNs+lRhERCrBwoUL6dy5M9HR0Vx11VWk\np6cD8N1339G0aVMA+vTpw4wZM7jjjjuIiopi1apVdO3alWuuuYaEhIRyi02JQUSkgmVlZXHppZcy\nePBgdu/ezeWXX84nn3xyeLtz/vvOpk2bxplnnsnLL7/M/v37admyZYXEp0diiMhxyeV7z2/Jlebm\n6p9++ons7GzuvvtuAP70pz/RtWvXwAQUAEoMInJcqsynZWzevJnGjRvnWRcXF1dJ0RwtIF1Jzrl+\nzrllzrkVzrmhBZR50Tm30jn3i3OuY3HrOufudc75nHN1AhGriEhla9SoEZs2bcqzbv369ZUUzdHK\nnBiccx5gFNAXaAcMcs61OaJMf6CFmZ0I3Aq8Wpy6zrkmwLmAHpEoItVGjx49CA4O5qWXXiI7O5sJ\nEyYwd+7cw9sLe/ifmZGRkUFmZiY+n4+MjAyysrICGl8gWgynASvNbJ2ZZQEfAAOOKDMAeBvAzOYA\n0c65hsWo+zzwjwDEKCJyzAgJCWHChAm89dZb1K1bl48//pg//elPh7e7XAMg7ojBkO+//54aNWpw\n4YUXsmHDBiIiIujbt29A4wvEGENjYEOu5Y34T/hFlWlcWF3n3MXABjP79cgPRkSkqjv11FNZsGBB\nvttydytNnz49z7ZevXrh8/nKNbbKGnwu9EzvnKsBPIC/G6nIOkOGDCE+Ph6AmJgYOnbseHgqv6Sk\nJAAta1nLx+Gy/CEpKYmxY8cCHD5fFqTME/U457oDj5pZv5zlYYCZ2dO5yrwKzDCzD3OWlwG9gIT8\n6gJfAVOBNPwJoQmwCTjNzLYfsX9N1CMiR9FEPX6VNVHPPKClcy7OORcKXAV8fkSZz4Hrc4LpDuw1\ns20F1TWz38ws1syam1kC/i6mTkcmBRERCbwydyWZmdc5dycwGX+iGWNmyc65W/2b7XUz+9o5d75z\nbhWQCtxQWN38dkMR3U8iIhIYmvNZRKoldSX5ac5nEREpMyUGERHJQ4lBRETyUGIQETkGaGpPERHJ\nV1FTez7zzDO0b9+eqKgoWrRoUS4zuSkxiIgcw46c2hPgnXfeYe/evUyaNIlRo0bx0UcfBXSfSgwi\nIpWgtFN7/v3vf6djx454PB5atWrFgAEDmDVrVkBjU2IQEalggZza84cffjiqRVFWmsFNRI5L7rHA\nPEzBhpf8JrpATe05fPhwzIwbbrihxHULo8QgIsel0pzQAyUQU3uOGjWK8ePHM3PmTEJCQgIZnrqS\nREQqWlmn9nzzzTcZOXIk06dPp1GjRoEOT4lBRKSilWVqz3fffZcHH3yQKVOmlLiVUVxKDCIiFaws\nU3s+/PDD7N69m65du1KrVi2ioqK4/fbbAxqfnq4qItWSnq7qp6eriohImSkxiIhIHkoMIiKShxKD\niIjkocQgIiJ5KDGIiEgeSgwiIpKHEoOIiOShxCAicgzQ1J4iIpKvoqb2fOGFF2jRogVRUVHExsZy\n4403kpKSEtAYlBhERI5hR07tOWDAAObPn8/+/ftZtmwZ69atY8SIEQHdpxKDiEglKO3UngkJCdSu\nXRsAr9eLx+MJ+KO3lRhERCpYWaf2fP/994mOjqZBgwY0aNDg8ExwgaLEICLHJ+cC81MKuaf2DAoK\nKvHUnoMGDWLfvn2sWLGCpUuX8sILL5QqjoIoMYjI8cksMD+lEIipPQFatGjBsGHDePvtt0sVR0GU\nGEREKlhZp/bMLSsri4iIiECEdZgSg4hUXWYweDC88kplR1IiZZnac8yYMezYsQOApUuX8tRTT+WZ\n/S0QlBhEpOp69ln43/9g5szKjqREyjK156xZs2jfvj1RUVEMHDiQwYMH87e//S2g8WlqTxGpmjZt\ngvbtYcwYGD4cFi/Os1lTe/ppak8ROX4kJ0OHDtC3L6xcCVlZlR1RtaHEICJV05o10Lw5RERA48aw\nalVlR1RtBCQxOOf6OeeWOedWOOeGFlDmRefcSufcL865jkXVdc6NdM4l55T/xDkXFYhYRaSaWL3a\nnxgATj4Zliyp3HiqkTInBuecBxgF9AXaAYOcc22OKNMfaGFmJwK3Aq8Wo+5koJ2ZdQRWAveXNVYR\nqUZyJ4Z27eC33yo3nmokEC2G04CVZrbOzLKAD4ABR5QZALwNYGZzgGjnXMPC6prZVDPz5dT/CWgS\ngFhFpLpQi6HcBCIxNAY25FremLOuOGWKUxfgRmBSmSMVkepDLYZyE1xJ+y32A0accw8CWWb2XkFl\nhgwZQnx8PAAxMTF07NiRxMREAJKSkgC0rGUtV6fljh0hI4OkJUtg6VISe/SAtWtJmjwZQkNJTEwk\nLi7uqHsAjkeHnryalJTE2LFjAQ6fLwtS5vsYnHPdgUfNrF/O8jDAzOzpXGVeBWaY2Yc5y8uAXkBC\nYXWdc0OAm4HeZpZRwP51H4PI8WbhQv8dz7nvXYiLg6QkSEiotLCqkvK+j2Ee0NI5F+ecCwWuAj4/\nosznwPU5wXQH9prZtsLqOuf6Af8ALi4oKYjIcSp3N9IhjRrB1q2VE081U+auJDPzOufuxH8VkQcY\nY2bJzrlb/ZvtdTP72jl3vnNuFZAK3FBY3ZyXfgkIBabkNAd/MrPbyxqviFQD+SWG2FjYsqVy4qlm\nAjLGYGbfAK2PWPfaEct3FrduzvoTAxGbiFRDa9bASSflXacWQ8DozmcRqXrWrDl6LEEthoBRYhCR\nqmfTJv9jMHKLjVWLIUCUGESk6tmyBU44Ie+6Ro3UYggQJQYRqVoyM2HfPqhXL+96tRgCRolBRKqW\nrVuhQQPwHHH60uBzwCgxiEjVsmWLPwkcqWFD2L4dfL6jt0mJKDGISNVSUGIIDYWoKNi1q+JjqmaU\nGESkaikoMYAuWQ0QJQYRqVryuyLpEA1AB4QSg4hULZs3F9xi0CWrAaHEICJVS1FdSWoxlJkSg4hU\nLYUlBl2yGhBKDCJStWjwudwpMYhI1eH1ws6d/nsW8qOupIBQYhCRqmP7dqhTB4ILmDFAg88BocQg\nIlXH5s0FX6oKajEEiBKDiFQdhY0vAMTEQEYGpKVVXEzVkBKDiFQdRSUG59RqCAAlBhGpOopKDKBx\nhgBQYhCRqqM4iUEthjJTYhCRqkOJoUIoMYhI1VHUVUmgrqQAUGIQkapDLYYKocQgIlWDzwfbtvlP\n/IVRi6HMlBhEpGrYtQtq1YKwsMLLqcVQZkoMIlI1FKcbCfQgvQBQYhCRqqG4iaFhQ9ixw//APSkV\nJQYRqRqKc0USQGgoREf7u56kVJQYRKRqKG6LATQAXUZKDCJSNZQ0MWzeXL7xVGNKDCJSNZQkMcTF\nwbp15RtPNabEICLlIzsbfvgBUlMD83olSQzx8bB2bWD2e6TFi+G226B/f7jkEnjrLUhJKZ99VRIl\nBpFq6mB6Cj99PppJ7z3OnN++xcwqbufr10OvXnDjjf7LR++6C9LTy/aaJU0MgW4xmME//gH9+kGz\nZnD33XDFFTBxIpxyCsyZE9j9VaIC5scTkRL7+Wd47z04eBD+/W+IjKzwEMyMX2Z9wpaXnqDjN78Q\nGxlOgxrh1N62j9Fn1+e8MUm0bNCmfIPwev3fpgcNggcegL17/d+we/SAjz+Gli1L/prZ2f7E0Lhx\n8crHxQW2xWAGf/4zJCfDb7/5pxc95OqrYcIEuOgieOkluPLKwO23kqjFIMWSvnEtP9/Ynzl/u4KV\nSRMKL5udzidLP+GL5V9wIONABUVYyZYsgb59/XfmpqTAmWf6H99QQTbt38SbY+9hascomvW/irpE\n4Jn0LfGb02j++25iflnOxZujmHlZN7YcKOerdd5/H4uJ4ZNL23DdZ4N5aunrrBo9Am6+GU4/HV57\nzZ8sSuL33/2XqtaoUbzyAWoxmBmj541m2l0X4l38C0yZkjcpHDJwIEydCn/7G7z9dr6vle3LLnM8\nFUWJQYq04uV/kdaqBfPnLmPxN4upedHlfNUzgRUrfsxTLj07nX999U/+flk9Gve7idi+1/PKBQ1Y\ntGF+pcSd6c1k9Z7VrNy1kt93/87KHWtZ8XsW8+b5z+Nl7dk4LDsbhgxh/yPDeKl/Xf5yeQSLTq5P\n5uDr/N80y0m2L5uPlnzEX57qyYLTExh47xu0vPQm6mzbT7cPZxJ7+rmHy7qWLWkyfR4Xrgvl7TvO\nwGe+8gkqKwt7dDh/PTOFp2ePpEeTHmw+sJnuY3rwdIcDZH/9JUya5O+KadYMTj0VHnoINmwo/HWX\nLYM2JWjpNGrkv4+hDP/JPvNx16S7mPXRM3QZP53u561nlyvk9Tp0gGnT/K2k118/vPrrlV/T+fXO\nRIyIYEJy4V+qjhlmVuYfoB+wDFgBDC2gzIvASuAXoGNRdYHawGRgOfAtEF3A65rPZ1JOtn470bZH\neKzfmc/a5Mlm8+ebPfvYFhsdd4ZtjXD27JXdbNTEF+y+14fZvf1ibENEuP14Qi/77KbP7Klzpti0\n6JNtfLdatjNlR4XEu2DzAvu/b/7POozuYGH/DLPYp+KszqMtLHxoc3P/19TcQzWsxt86WtTg6y00\n8RnresUU+2HaNrODB0u/02eftV09O1u9p+vatROutedmP2dXv3eZ/dYo2L5//GbzlcMB+sXyL6z3\nY81tWpd6llY32jJGPmmWmlpkPW/yUttTM9je+ebfAY/JzMw3bpwtahNr8fddZhdd7LXbbjNbudJs\n9e7V1mdcH+v8Wmf7Zcsv/ljXrTObOdPsnnvMGjc2+/XXgl/4ySfN7r23ZMG0aGG2fHmp38tzs5+z\n85471fbXbWJPnfGZNbjtGvvLuOeLrrhypVlcnNkTT5jXm20tX2xp7y1+35JWzrX6I+vblN+nlDqm\noizZvsQ+Tf7UVu9eXWRZ/+m/gHN6QRuK+4O/1bEKiANCck78bY4o0x/4Kuf3bsBPRdUFngbuy/l9\nKPBUAfu3mBiz7t3Nhgwxe/pps88+8x8PWVll/ZgD6+BBsyVLzCZPNps0yWzuXLM9e8ppZ16vlTVj\nZmxcZ1tqhti1Zw2xlJQjtmWYfTT8O5t6wkm2LzjE9oaE2pyTzrMtn8/NU+71Z/fb4tr17c3BncoU\nS758Pv/JZMwYW3/XYPvknMb235417b2zz7Qnu9xvl9b83C478Rd7+Jat9u072+zAdz/bwbfftI33\n3Gi/9+lsG5vVsfRgZ2nBWHqQs9Ra4bb3ykvM5swpfgx79lh6nWjrcW+MTV89Pc+m374db7tqBdvA\n57rb0u1LA/KWN+3fZFe8e6k9OaCeZcREme+xx+yo/5xcdu82++ors//8x2zECLMPPjBbc911Nq5b\nuO1IDWyyzsw0+7VxS/tT/3h7dESKTZxo9uijZnXrmt1+u9mWLT574+c3rN7IevbQtIcsLTPtj8rv\nvWfWsKHZli35v/j115u98UbJAurd2//HVgrLdiyzuk/UsTmNetlb9f5u48ebPfJmkgXfc5I98UQx\n/q42bjTr2NHWX3KOtby3g9Wt57OQELPGPadb5CONbNHynD/8PXvMkpPNFi0yW7Wq1Cet5B3J1n30\nWdbttnp21WWdLeHvdey+bx4q9EtJYYnBWRmbus657sBwM+ufszwsZ4dP5yrzKjDDzD7MWU4GEoGE\nguo655YBvcxsm3MuFkgys6Paks4527HDWLbMPy60bNkfP5s3Q/Pm/hZomzb+R6jUrJn/T2Skv3s4\nKqoMH8b27aRvWMv2PRs5UD8Ki20ImdHMT4rlswmhTJvmb+E2aQLBwbBzJ6xY4Y+rWzfo2tVo0m49\nmXUXsHTvz8zfPJ/ftv+Gz3xEh0fTvHZzejTpQZ+EPnQ5oQshQSF/7HvDBvjgA/jqK5g3D9LS8NVv\nyN72Pdl541Ci+nSlYUP/XOnFYsb8Dq2YGgRDvllBbGwhFbOy/AOO4eH5vQw39/2Np2eewrrZX3Jq\nx/4l+0zzs3YtjB8P77xDZnoqs5s55tXYT1jGBWyZ35lurVI4JXoNjdlM6O6t/idt+nz+D75VK2jb\n9o+fE09kU6pxxd9mk5YxgQsPjuee2VmkXnAeTUe/i6dmrUJDWXfndcya8zHNJybRvUn3o7b7HnqQ\ntd9/zmkXbObWLrfx1+5/pX5k/RK/ZZ/5+O/P/2XCuKGMnRROg4STCRr96lEDuQcPwsKF/gtkPv/c\nPx7etav/rUZGwsqVsChpD3PSG/PUgwN55sHxJY4lP3v3wu39vuOZJb1ZM2cuZ5zU+fC2nTvhiSdg\n3Dj/xUlX37qZB2few/frvuf2LrdzY6cbaRrdFO65x/+H8eyzR++gWzd47jk444ziB3XTTdC9u39s\nIz+7d8OYMVCvHvTu7R+wxt8F2W10T+56NpyemR7ilk8hrGYIZkbL/7Rm/9tj+d/zp9OrV+G73785\nhU/6tuPiVXsJvf7PRHbvwKb1Xr6c+gpNV+2k+54somw/B2s3whscRnDqfsJTdrClUWe297qc1v+6\njqi42oXuw8z47/zXmfvU33nw6zBqhkaT1SiOWqsX8NYpIewacTuPnfdYvnWdc5hZvn/YgUgMfwL6\nmtktOcvXAqeZ2d25ynwBPGlms3OWp+BvBSQUVNc5t8fMaud6jd1mdtSoj3POxn/0MO5ACi4lBU9K\nKp6UVIJT0nAHDuLbfRDbm47nQAaezGzwGXgN5zPwGc5rOMv53QdpweFk1oomI6Yu2XUa4o2JxVO7\nPiF16hNetz41Ij2EB2UQ7knHs2UpIUvmU2vpck5YuYXgjGw21DIsOIRG+7x4cXwbF8r/2mcwu3Vd\nmsQ2JaFOUxpGNiQsKIyIDB8tlmym3uyVNFixibp79pLh8bA5JIaltVqzs/GZ1OxyEXFnNKNe073s\n9K3i550/MH3NNH7f8zsXhJ3CwLW16Tx9BQ3XbmJm605MatSOqSFNWbM3jKYZexl4cAX3JE9nmuc8\nHqrxKucMjOL226FTp8L/X+fdN5Sgcc+x44PV9D27aZmOkc2b4X+nnUez5r8x4LtNuGJnpyP8/DP8\n61/www+kD7yYV9oc4F8ZM2i6YShbPr2TO26twR13+P/OS8oM/vlPeGv8Qe5+YgwnPfUQddId0VN+\noFWjk/Ots3jxFJr06MvSb8fTs+fV+b9wRgZ06cKeO27i/xotYmLyRE5tdConNziZNvXa0LZeW9rU\na0NszdgCP5cFWxbw4IQ7uemT1QxY4iXk+f+QffkgFix0/PgjrF7tH2ddvRpWrfInga5d/RcG9e17\ndL5OS4MfB48g8+fHeKf/Ml56vDl165b8Mztk1y447zy4x7Whcas69Hlvdr7l1q6Fhx/2X9156qkQ\nFLuUlXVeYlu9jwhNiydhawKzJn7J0CduIeKkEHam7WDj/o1s2r+Refet4owHG5FSK5SaoTWJi47j\n5AYn0+WELrSM6MzOVfGsWeP//MLC/I9Lav3xP/GlpjPjnBGkpPhvp0hN9f9fxwdv5NaJfQnu1J7I\nmBDc11/D4MHYTX/muo9GkvhBEhem1SN2yfQ83xZHzhrJ9F9Wsuzf/2XRIv8jmfJjBhdes55pLU9l\n94BJREye6r+aKTSUg/XrcNfWNzml48usXj+IXbsdQUFQuzY0jEyh8aoZ1P/pBbpsnMvc3hdx1lvP\nEnnC0Zfp7j64myeeHMig/84lJLUldd8YTePLTwfnyN62iyXtL2Z5g/l0+2E5cbXjj6p/LCaGqcB9\nlCwx7DKzow5f55xdUSOEE8JD8IYEExkeSqt6tejYtAHeyBrM3ZeO1QijY9sECAtjwdrN4Dx0OrEp\neDwsXL0Jcx46tY7DnOO72UtI3byHNh6I2LeHtTt2EJ6RyVnOR1R6Fj9m+/A6R6ewYDaGR/I/6rDS\ntWBT6BVsTDmT8My11IkJonfvRAa0X03Y8v8QOnMavdasIbXdiXzt8ZLhy6bv/lRqr9vGxBPqsLdt\nPKdecj7xp5zFkl9W4t26nfabM8mcPY/FybPI9npoENyePdlR/OLdSaQ7yGVsI9hzkNFRMcyOq8va\nnq2JqB2JZ+dOomsF0bpHM8LDYc3C1Wze8Tu3zVrOeb97eKDHVXw563puvbU3jz0Gs2cnAZCYmAhA\nUlIS+1as4oy/3srLdz9Pr34djtpemuWkz9pxx+uNePehO+jY49Iiy7fq3IpZ62cxdfpUGmzayy0z\nk6m7chOjz2zDl62Mn2stIX7/1ax5ox9XXBzNiy8mEhlZ+vgOLd97bxJvvQUffdCDhBFdGb8pmU13\nDeblO0cTEhRyuHzsybHMvLQzNWp0oPEDTxb4ejNmJJH222rOf3wY7qMPmZSdzuJtiwltEUryzmR+\n+/YHYpatpUsINIyoz1qicY1OoPGZ7djEfn6akkTzZVsZtTKE8AsvY0TMRXwzO4rk5ESaNoXmzZNo\n3BjOOSeR+HjYuTOJ0NBivN/TT2dnQiOubd2UOQtfYPjwRG65BebOLdnnNW5cEg88AL2vSeHFlwcw\n59UxhDeNL7R+WhqEhCRy8CAsX54EQRnEtI1m9YGl1Hr0ebbtg6/6XM+gAfWJydhFw4NBDLp/BBtX\nzGfWD7NIy0oj5sT6vD9jMVNnfMP+0OUEn2g0yO5CjdX1iUhvToPaF3L+miSCV3/KZ2dcR2ybBHxh\nu9myYQ5m2bz5+VdMb3ABQzddgmWHcHab1ly7/Wm8G14nnHQ6Xvxnmr72JEm//pon/rGfjmXY1GFc\n6t1KairceGP+n8+iRYk8n/QmCR3eZfjZDx+1PbtZNtdNvI7nWz9PbM1YEhMT8fq83D/mft5d/C6x\n7WPpsLMRsf+dTfdNe8g+6xSaPTSUvV7/t55gW8e+e+8ma5mxsNM9PJr0BKHhnjyf98F9mfwvti6z\nBpzMqx/8SFJSEmPHjgUgPj6exx57rMDEEIgxhu7AN7mWh3HEADTwKnBlruVlQMPC6gLJQMOc32OB\n5AL2X2T/W3nLzDRLTy+i0N69ZtOmmb3zjtmYMWY//miWllZEJfP3o69fbzZ9utlnn1n2hM8s9Ytp\nlrlyrX8coZjSs9Jt/isP286oYHu1fyPrOehbO+WUo8f7fJlZ9nOjhvZ4YteADurv2WP2j6Y3W1Kb\nSMvyFtyPejDroD0w9QGLfqKOdX/4AvuoYwfbERlhIy7obl0eG2hdHvw/a3vV21ar/h674w6zbdsC\nF+Mh06b5u7tHPpZmad262Mf946zFf1rYc7Ofs8+XfW6PTH/EThkabQejI/MNYOtWsxdeMDvvPLPI\nSLNatcz6hU23XcH1bWH/+y3rnffN7rjD7KSTzKKjzS680FLvucPW3zrIVvTvZmtPbmp769Wyg7Ui\n7GDDupZyxfX2+pBZFhtr1qePvzs+EO877aN3bVlssL05abJdeqlZgwZm991n9ttvRQ9PZWWZvfii\nWb16Zq++kWb3XtvAtp/ZuexBrV1rvjp17K2XDliTJmYXXWQ2+/Eplt3zLPP5/GN0jzzi//8591yz\niRP9f0ab92+2L5Z/YcNnDLeBHw60zq91tqvvaWLzmodbl9e72HnvnGdXfnyl/eXLv9grj19s65rF\nWPuXT7aIERHWcdRp1vmFc63eiKZ27n8H2IYdawoMz+fzWf2R9W3JpjXWqpXZhx8eXWbWLP9neenb\n19lr818r8LWe//F5i38h3p764Skb8f0Ia/dyOztjzBn2w7ofDo8N+HxmH77yg7144pm2LSLYtkZi\n66Od7aoRZE82usa++aTg8SUzs/9eMc6S63ls256NR22jnAefg/hjADkU/wBy2yPKnM8fg8/d+WPw\nucC6+AefDyWJQgefpfh8mzbZltM72MK4MBtwz+kW3Xqh/fWvZkuXmiUvSLUp8W3t6+Y1bPOOwI+K\n/3tEuq2qXcMmvfTXfLcfyDhgvd7sbfHDBthFsR/Zzppxtqj7LfbcI7vt/vvN7rrLP4D62WfFugCn\nTNavN+vVy6x/t12WlnCSrfvbTTZk4mC74N0L7NZPbrTU83qbPfzw4fKpqf4Tdv/+ZjEx/rHSiRPN\ndu3yb09PN/tpzG/2ftx9NjXiIpt/5dOW/ePcQgcbly41u/FGs9q1/YO3SwMzfv0Hn8+2JHaxl/vW\nsfSsdFu2zOwf/zBr1sx/kdCwM2faT51us/l9H7BvRsy3jz82e+klsz//2axpU//Y7q+/mt399d2W\n3Lqu2YQJgYlr4ECzUaMsLc3stdfMnm3+kr3uudUiI/0X+9xxRzE/iy1b/B/ekZ9xr15m775rZmYp\nGSn2/drv7ZuV39iv2wq5KiqXyz+63MYuHGtz55rVr2+2YMEf27ZuNWvSxOyzz3zW9LmmtnxnwVdF\n+Xw++27td3bbF7fZXV/fZVN/n1rgYHFmptn4d3x24/nL7bJuc23ofT7btKnoWDdt9NmMBg3sq4dv\nOmpbuSacj18fAAAPaklEQVQG/+vTD/9lpSuBYTnrbgVuyVVmVE4SWAScWljdnPV1gKk52yYDMQXs\nu+hPR/Lyei3rxRcspU4tm9wmzIafE2+PnHK6/VonwiaeXNdWbSz6UrfSSEszuz7h37a8fpBt2Loi\nz7Y9B/dYt9dPtwZ/vsHGdn7evA0amn39dbnEUVzZ2WajR5t1arjJlkV1sd+7XmnrXvzUMi+61DL7\nXmBzvk+3Z5/1n8diYvwthPHjC71IyMzMZsww69nTrG1bfzLJ3XDMyjL79luzCy/0f+t8/HGznTvL\n8U1u3mx7osPs6ScvtGxvtpmZ+dastQN9LraU+nH2xTkP2QddB9uesNr2ZLf/2W23+a9wWrLE/212\n1JxRduGDCeaNjfWfvQLhu+/MWrf+48q6s84y79vjbffuUlxo16mT//UOmT3bLD6+TJcsvjz3ZRvy\n6RAzM/v4Y7NGjcy++MLf0mzWzP/lZdWuVdbomUblcqlySb11Skcbd8mAo9YXlhjKPMZQ2ZxzVtXf\nQ6VJSyP74w/ZMOUT9h3cQ2rf3nQZ8gBhocW8u7QUxrxhNPxnB1Ja7+ayb9YR7Alma8pW+r9zAXt/\n6s74n1I5PXwh7tNPISGh3OIoiYMH4esJBwn612NEbVjK6owTuC/0PzQ7MYwzzvDfzNu7d/Ef4wP+\nwcnJk/1Pzpg37497t5YuhRNPhL/8Ba65BiIiyuc95Zby6UdkDr6WBaeeQPvoE6n7/Tx+vLwHQ0/Z\nzvKUtbSq24qGyzby6mubef+WHgRffiWNoxozIXkCM9bOYNnc06jVqj08/nhgAjLzj07feqv/GPjr\nX+HXX/1XLJXU8OH+0fZ//9u/fMkl/pHy228vdXhLdyzlgvcuYM09awB4/3144w3/BQAjR/pvgh6z\nYAzT107n3YHvlno/gTK2+2m42nUZPGlSnvXlOvhc2ZQYqhafD646dzsvzmvGs5fWZvsl5/L58i/p\nPOtq3pgzh7heCbixb1XKc4Yqy86d/ic+eL3+BJHfExfKW8a2zXw/9Crmus3MaOGhVadzuKTNJfRJ\n6EOQJwiAzT98Te0BV/Hqvb2YmuDjzGZnclfqyUTefLs/o9WsGbiAVqyAs8/2//6f/8Bll5XudebN\ng+uu81+/vnSpP4OvWVP8R2vkw8xo+ExD5t48l/iY+HzLXDfxOs5qdhY3dy7gUtkKNKZnT0LDQrlu\n2vQ86wtLDAHpSqrMH9SVVOXs3Gl2fqP5tq52gi1KaGHTInrZgfB65n3l1TLflCflbMYMf8f66NH+\nCyJOPNE/6FMeli41u/POEl1kcRSv1yw21j8YMmiQv58nAA6NM+TH5/NZk+ea2IqdK/LdXtFe69Xb\n3jnrjKPWU0hXkp6uKhWubl14fV5nPnkvmZCvPqXP3bWoOaiz/04/ObYlJsLMmf5+Lo/H3+d18cXl\ns6+2bf1PKy0LjwcGDPDfuNO1K7z8ckBCS4xPJGldEoM7Dj5q2+o9q/GZj5Z1SvEU2XJgIaEEpZZs\nvgh1JYlI9Zae7n/QYQC7uo4cZ8jtjQVvkLQ2ifEDA3NXeVm9cv4l1Nu5jivmLsyzvrCuJD1dVUSq\nt/DwwI5/AG3rtSU1M5V1e49+tHfS2iQS4xMDur+ysJBQgrOySlRHiUFEpISccyTGJ/Lduu/yrDez\nYy8xhIYTnK3EICJS7hLjE0lam5Rn3e97fsc5R4vaLSonqHwoMYiIVJDeCb2Z/PvkPJMeTVo5id4J\nvUv/oMjyEBpOsLdks8cpMYiIlEKbem2oH1mfGWtmHF739uK3uab9NZUY1dGcEoOISMUZcsoQxi4a\nC8CS7UvYcmALfRL6VG5QR3Bh4YRmKzGIiFSIq9tfzRfLv2B/xn7GLRrHtR2uPXyn+LHChdcgpIQt\nBt3gJiJSSvUj69M7oTdnvHkG6/et56ebfqrskI4SFBZBqNdbojpKDCIiZfDWgLdYvG0x9SPr06be\nUbMPVzpPjRqEZCsxiIhUmOjwaM6MO7OywyhQcHgNQny+ogvmojEGEZFqLKhGBKFeJQYREckRHBFJ\nWAnHGJQYRESqsdCISEK9JXvQqBKDiEg1FhIRSVi2upJERCRHWGQkoT61GEREJEd4ZARhXvxzaReT\nEoOISDVWo0Y4mR4gM7PYdZQYRESqsYiwUDKCgYyMYtdRYhARqcYiwkLICEKJQURE/CLD1WIQEZFc\nwkODyAgCb/rBYtdRYhARqcZCQx2ZQZB9MLXYdZQYRESqsZAQyAjyKDGIiIhfcDBkBDkyU5UYREQE\ncA4yPB7SU5QYREQkR0aQh6zUtGKXV2IQEanmsjweMpQYRETkkAxPENlpulxVRERyZAQFqStJRET+\nkOkJIvugWgwiIpIjMygYb0V1JTnnajvnJjvnljvnvnXORRdQrp9zbplzboVzbmhR9Z1z5zjn5jvn\nFjnn5jnnzi5LnCIix7PMoOAKfSTGMGCqmbUGpgP3H1nAOecBRgF9gXbAIOdcmyLq7wAuNLNTgCHA\nO2WMU0TkuJXtCcbSK+4hegOAcTm/jwMuyafMacBKM1tnZlnABzn1CqxvZovMbGvO70uAcOdcSBlj\nFRE5LmV6gvFlpBe7fFkTQwMz2waQcyJvkE+ZxsCGXMsbc9YBNCyqvnPuMmBBTlIREZESygoKwdKL\nnxiCiyrgnJsCNMy9CjDgoXyKl2zG6SLqO+faAU8C55bxdUVEjltZwSElmo+hyMRgZgWelJ1z25xz\nDc1sm3MuFtieT7FNQLNcy01y1gFsLai+c64JMAG4zszWFhbjkCFDiI+PByAmJoaOHTuSmJgIQFJS\nEoCWtaxlLR+3y1lBoSzcsoU3hgwBOHy+LIgzK/2XfOfc08BuM3s652qj2mY27IgyQcByoA+wBZgL\nDDKz5ILqO+digCTgUTP7tIgYrCzvQUSkunugVw/OD65Bz2nTD69zzmFmLr/yZR1jeBo41zl36MT/\nVM4OGznnvgQwMy9wJzAZWAJ8YGbJhdUH7gBaAI845xY65xY45+qVMVYRkeOSNygET2bxu5LK1GI4\nFqjFICJSuHvPO5fL9u2lx5x5h9eVZ4tBRESOcd6gUDxZmcUur8QgIlLNeYPDCMos/hX/SgwiItWc\nNyScILUYRETkEAsOxZOdXezySgwiItWcLzic4Gx1JYmISA4LDSdYLQYRETnEQmooMYiISC6hNQhR\nYhARkcNCwwnJ9ha7uBKDiEg15wmLIMSrxCAiIjlcaA0lBhER+YMLjyA021fs8koMIiLVXFBYJKFq\nMYiIyCEhYeH+X4p5ZZISg4hINRceEkq2x0ExWw1KDCIi1VxYSIg/MajFICIi4G8xeB1KDCIi4hce\nEopXLQYRETkkPCSEbA8aYxAREb8aoaH+xKAWg4iIAISHhigxiIjIH2qEhpLtTIlBRET8DnclaYxB\nREQAaoSGkO1Ri0FERHJEhmvwWUREcqkRGkJ2kFoMIiKSIzwsCK+D7Mz0YpVXYhARqeZCQiDb4/Bl\nZRWrvBKDiEg1FxIC2Q582ZnFKq/EICJSzYWG5rQYMpUYREQEf4vB65xaDCIi4ufvStIYg4iI5Dg0\n+OzNUotBREQA5yDbecjOUGIQEZEc2R5Hdoa6kkREJIcXD96KuCrJOVfbOTfZObfcOfetcy66gHL9\nnHPLnHMrnHNDi1vfOdfMOXfAOfd/ZYlTROR45x9jqJgWwzBgqpm1BqYD9x9ZwDnnAUYBfYF2wCDn\nXJti1n8W+LqMMYqIHPe8zoMvs2ISwwBgXM7v44BL8ilzGrDSzNaZWRbwQU69Qus75wYAq4ElZYxR\nROS4l+1xFdOVBDQws20AZrYVaJBPmcbAhlzLG3PWATQ8on5DAOdcTeA+4DHAlTFGEZHjnpfitxiC\niyrgnJtCzgn70CrAgIfyKW7F2mvBfDn/DgeeN7M059yhfRZoyJAhxMfHAxATE0PHjh1JTEwEICkp\nCUDLWtaylo/r5aVpXr78ZCIxa9YdPl8WxJmV/lzunEsGEs1sm3MuFphhZm2PKNMdeNTM+uUsDwPM\nzJ4uqL5z7nugSc5L1Aa8wCNm9ko+MVhZ3oOIyPHg5Q71Ofv8GzjpqZEAOOcws3y/dJe1K+lzYEjO\n74OBz/IpMw9o6ZyLc86FAlfl1CuwvpmdZWbNzaw58ALwRH5JQUREiifb48EqaPD5aeBc59xyoA/w\nFIBzrpFz7ksAM/MCdwKT8Q8kf2BmyYXVFxGRwPI6T7GflVTkGENhzGw3cE4+67cAF+Za/gZoXdz6\nR5R5rCwxioiIPzFYlqb2FBGRHF5PED7N+SwiIod4nQfL1rOSREQkh9fjAXUliYjIIT4XhKkrSURE\nDsl2QZhXiUFERHL4XBBoak8RETkk2xME6koSEZFDfEoMIiKSmz8xeItVVolBROQ44HVB4FNiEBGR\nHD5PEE5dSSIicojPBYEuVxURkUP8LQZ1Jclx6tCsVSKH6JgAnycYpzEGOV7pJCBH0jGRkxi8SgzH\nhPI+IMv6+qWpX5I6xSlbWJnSbjuW6ZgouqyOicC/vgafjyE6CRRdVieBY+v1dUxUvIpJDCF4itli\ncGZWrgGVN+dc1X4DIiKVxMxcfuurfGIQEZHAUleSiIjkocQgIiJ5KDGIiEgeSgwiIpKHEoOIiOQR\nXNkBlAfnXATwCpABfGdm71VySFLJnHMJwINAlJldUdnxSOVzzg0ALgBqAW+a2ZRKDumYUS0vV3XO\nXQvsMbOvnHMfmNlVlR2THBuccx8pMUhuzrkY4N9mdnNlx3KsqBJdSc65Mc65bc65xUes7+ecW+ac\nW+GcG5prUxNgQ87vxbvVT6qUUhwTUs2V4Zh4CHi5YqKsGqpEYgDeAvrmXuGc8wCjcta3AwY559rk\nbN6APzkA5Htnn1R5JT0mDhermPCkEpT4mHDOPQV8bWa/VGSgx7oqkRjMbCaw54jVpwErzWydmWUB\nHwADcrZNBC5zzr0MfFFxkUpFKekx4Zyr45wbDXRUS6J6KsUxcRfQB/+54pYKDfYYV5UHnxvzR3cR\nwEb8BwFmlgbcWBlBSaUq7JjYDfylMoKSSlXYMfES8FJlBHWsqxItBhERqThVOTFsAprlWm6Ss06O\nXzom5Eg6JkqhKiUGR96Bw3lAS+dcnHMuFLgK+LxSIpPKomNCjqRjIgCqRGJwzr0HzAZaOefWO+du\nMDMvcBcwGVgCfGBmyZUZp1QcHRNyJB0TgVMtb3ATEZHSqxItBhERqThKDCIikocSg4iI5KHEICIi\neSgxiIhIHkoMIiKShxKDiIjkocQgIiJ5/D8JK4IogjHNiAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "('z340', 'READ 2018-10-25/X133981.340')\n", "('zdat', 'READ calib_data/calib2018-10-25_p0_c6.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p1_c6.dat')\n", "('zdat', 'READ calib_data/calib2018-10-25_p2_c6.dat')\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD/CAYAAAD12nFYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VGX2wPHvOzOpQBIChA6hgwiEFkApUVRAEawouyrN\ntoqu4iK4/hRccZWisjQLoohlLYiKokiNIAgIEVykl0AINQmBkD4z5/fHhBLSJsmQej7Pk4fcO+97\n75mZyz15y73XiAhKKaXUeZbSDkAppVTZoolBKaVUNpoYlFJKZaOJQSmlVDaaGJRSSmVjK+0AissY\no9OqlFKqCETE5La+QrQYRKTM/kyYMKFMb78o9QtTx52y+ZUpymtX+jMv7e9Mjwk9Jjyx/fxUiMRQ\nlkVERJTp7RelfmHquFM2vzJFfa0s02Oi4LJ6TJTu9k1BmaOsM8ZIeX8PyrMmTpzIxIkTSzsMVYbo\nMZGTMQapyF1JSl2qvP7VqK4cPSYKR1sMSilVCeXXYij3s5KUUio3oaGhHDp0qLTDKHWNGzcmOjq6\nUHW0xaCUqpCy/iIu7TBKXV6fg44xKKWUcpsmBqWUUtloYlBKKZWNJgallFLZaGJQSqkyYMSIEbz4\n4ov88ssvtGnT5sL6PXv20LFjRwIDA5k1a1aJxKKJQSmlypCePXuyc+fOC8tTpkzh+uuv58yZM4we\nPZrIyEiuv/56goKCaNq06RWJQRODUkqVYYcOHaJt27YXlqtUqcKoUaOYNm3aFdunJgallCoFv//+\nO507dyYwMJB7772XtLQ0AH7++WcaNmwIQN++fVm9ejWPP/44AQEB7Nu3j65du/LXv/6VJk2aXLHY\nNDEopVQJy8zM5Pbbb2fYsGEkJCRw991389VXX1143RjXdWcrV66kV69ezJ49m7Nnz9K8efMSiU9v\niaGUqpRMrtf8Fl5RLq7esGEDdrudJ598EoA777yTrl27eiYgD9DEoJSqlErzbhlHjx6lfv362dY1\nbty4lKLJySNdScaY/saYXcaYPcaYcXmUmWGM2WuM2WqMCXO3rjHmGWOM0xgT7IlYlVKqtNWtW5fY\n2Nhs6w4fPlxK0eRU7MRgjLEAs4B+QFtgqDGm9WVlBgDNRKQF8Ajwtjt1jTENgBsBvUWiUqrC6NGj\nBzabjZkzZ2K321m0aBGbNm268Hp+N/8TEdLT08nIyMDpdJKenk5mZqZH4/NEiyEc2Csih0QkE/gM\nGHxZmcHAAgAR2QgEGmNqu1H3TWCsB2JUSqkyw8vLi0WLFvHBBx9Qo0YNvvzyS+68884Lr5tLBkDM\nZYMha9aswc/Pj4EDBxITE4O/vz/9+vXzaHyeGGOoD8RcsnwE1wm/oDL186trjBkExIjI/y7/YJRS\nqrzr1KkTUVFRub52abfSqlWrsr3Wp08fnE7nFY2ttAaf8z3TG2P8gH/i6kYqsM7w4cMJDQ0FICgo\niLCwsAuP8ouMjATQZV3W5Uq4rC6KjIxk/vz5ABfOl3kp9oN6jDHdgYki0j9reTwgIjL5kjJvA6tF\n5POs5V1AH6BJbnWBJcAKIAVXQmgAxALhInLysv3rg3qUUjnog3pcSutBPb8BzY0xjY0x3sC9wOLL\nyiwGHsgKpjuQKCIn8qorIttFpI6INBWRJri6mDpenhSUUkp5XrG7kkTEYYwZDSzDlWjmichOY8wj\nrpflXRH5wRhzszFmH5AMjMivbm67oYDuJ6WUUp6hz3xWSlVI2pXkos98VkopVWyaGJRSSmWjiUEp\npVQ2mhiUUqoM0Ed7KqWUylVBj/acNm0a7dq1IyAggGbNml2RJ7lpYlBKqTLs8kd7Anz00UckJiby\n448/MmvWLL744guP7lMTg1JKlYKiPtrzH//4B2FhYVgsFlq2bMngwYNZt26dR2PTxKCUUiXMk4/2\nXLt2bY4WRXHpE9yUUpWSeckzN1OQCYW/iM5Tj/acMGECIsKIESMKXTc/mhiUUpVSUU7onuKJR3vO\nmjWLjz/+mF9++QUvLy9PhqddSUopVdKK+2jP999/nylTprBq1Srq1q3r6fA0MSilVEkrzqM9P/nk\nE55//nmWL19e6FaGuzQxKKVUCSvOoz1feOEFEhIS6Nq1K9WqVSMgIIDHHnvMo/Hp3VWVUhWS3l3V\nRe+uqpRSqtg0MSillMpGE4NSSqlsNDEopZTKRhODUkqpbDQxKKWUykYTg1JKqWw0MSillMpGE4NS\nSpUB+mhPpZRSuSro0Z7Tp0+nWbNmBAQEUKdOHUaOHMm5c+c8GoMmBqWUKsMuf7Tn4MGD2bx5M2fP\nnmXXrl0cOnSIV155xaP71MSglFKloKiP9mzSpAnVq1cHwOFwYLFYPH7rbU0MSilVwor7aM///ve/\nBAYGEhISQkhIyIUnwXmKJgalVOVkjGd+iuDSR3tardZCP9pz6NChnDlzhj179rBjxw6mT59epDjy\noolBKVU5iXjmpwg88WhPgGbNmjF+/HgWLFhQpDjyoolBKaVKWHEf7XmpzMxM/P39PRHWBZoYlFKq\nhBXn0Z7z5s3j1KlTAOzYsYPXXnst29PfPEETg1JKlbDiPNpz3bp1tGvXjoCAAO644w6GDRvG008/\n7dH49NGeSqkKSR/t6aKP9lRKKVVsmhiUUkpl45HEYIzpb4zZZYzZY4wZl0eZGcaYvcaYrcaYsILq\nGmOmGGN2ZpX/yhgT4IlYlVJK5a/YicEYYwFmAf2AtsBQY0zry8oMAJqJSAvgEeBtN+ouA9qKSBiw\nF3iuuLEqpZQqmCdaDOHAXhE5JCKZwGfA4MvKDAYWAIjIRiDQGFM7v7oiskJEnFn1NwANPBCrUkqp\nAngiMdQHYi5ZPpK1zp0y7tQFGAn8WOxIlVJKFchWSvt1+wYjxpjngUwR+TSvMsOHDyc0NBSAoKAg\nwsLCiIiIACAyMhJAl3VZlyvZcuPGjXNcA1AZnb/zamRkJPPnzwe4cL7MS7GvYzDGdAcmikj/rOXx\ngIjI5EvKvA2sFpHPs5Z3AX2AJvnVNcYMBx4CrheR9Dz2r9cxKKVKXNjbYXww+AM61u2Ybf0tn97C\nY10e45aWt5RSZO650tcx/AY0N8Y0NsZ4A/cCiy8rsxh4ICuY7kCiiJzIr64xpj8wFhiUV1JQSqnS\nEpcSR03/mjnW1/CrQVxKXClE5DnF7koSEYcxZjSuWUQWYJ6I7DTGPOJ6Wd4VkR+MMTcbY/YBycCI\n/OpmbXom4A0sz2oObhCRx4obr1JKFZeIEJcSRw3/Gjleq+lfUxMDgIgsBVpdtu6dy5ZHu1s3a30L\nT8SmlFKelpKZgsVY8PfKeVfTipAY9MpnpZQqpLy6kcCVGOJT40s4Is/SxKCUUoVUUGLQFoNSSlUy\neY0vgCYGpZSqlOJT47XFoJRS6qK4lDhq+mliUEoplSW/MYZgv2BOp53GeeFWb+WPJgallCqk+JT4\nPMcYbBYbVb2rkpiWWMJReY4mBqWUKqS41LxbDFD+u5M0MSilVCHl15UEmhiUUqrSiUuJo4Zf7l1J\noIlBKaUqnfiUvKergiYGpZSqVPK7gd55Nf1qEp9Sfm+LoYlBKaUKITkzGavFmusN9M7TFoNSSlUi\n8Snx+Y4vgCYGpZSqVAqakQRZiSFVE4NSSlUK7iSGGv7l+ylumhiUUqoQChp4Bu1KUkqpSiU+NT7P\nG+idp4lBKaUqEXe6kqr7VudM2hnsTnsJReVZmhiUUqoQ3EkMVouVIN8gTqeeLqGoPEsTg1JKFUJ8\nat53Vr1Uee5O0sSglFKF4E6LAVyJIT61fF79rIlBKaUKoTCJQVsMSilVCRR0Z9XzaviV32sZNDEo\npZSbRCTfp7ddSlsMpWzu3NKOQClVGbhzA73zNDGUshdegB9/LO0olFIVnbvjC6CJodR9/TU88ABs\n2VLakSilKjJ37qx6niaGUtajB7z7LgwaBNHRpR2NUqqiqiwtBltpB+Apt98OMTEwYACsWwfBwaUd\nkVKqojmVcsrtxBBSJYQTySeucERXRoVoMZz35JNw882uJJGeXtrRKKUqmtizsTQIaOBW2foB9Tma\ndBSnOK9wVJ5XoRIDwNSpULs2DBsGzvL3fZR54nQS/9sakt6ajvzwA2RmlnZISpWYI2ePuJ0YfG2+\nBPgEcCr51BWOyvMqXGKwWGDBAoiNhfHjSzuaimVP9BYie9Yn9cYIfnj3Wf43+m7SGtaFzZtLO7RK\nQ0TYG7+XJXuWlO0btP36KzzzDIwYASdPlnY0HnMkyf3EANAwoCExZ2OuYERXRoVLDAC+vvDtt7B4\nMcyeXdrRVAy7923E3qMbfv6NmTjwNK8407mv60eMut5C8o0RsH59aYdYJqTb09kfv++KbNvutPP4\nl8N58bnurJr9DDe/06tsDm4ePOiaCRIc7Prp1g127Cj0ZmLOxLBo5yJE5AoEWTQxZ2IKlRgaBDTg\nyNkjVzCiK0REyvWP6y3k7sABkXr1RL75Js8iyg2nTsfK+ha+8v01N0rNmiKvvy6yaZPrc40YHCO3\n315fTtcMFDl7ttRidDqdcjr1tDidzlKL4cdtX8lnPapJuhXZX8dH9kx8UsRD8WQ6MuXVsT3keA1f\nsV/bQ5zXXitHm9SS/v9uKxn2DI/swyMcDpHevUWmTROH0yHRp6PFOXeuSLNmIomJbm3C6XTKs8ue\nleDJwdJyZkt59LtHxe6wFzmkAwkHZH/C/iLXv1TtqbUl9mys2+Uf+/4xmbFhhkf27WlZ585cz6sV\nssVwXpMmrpbDQw/Bxo2lHU359dtf+uCw1OO5xB+JioIxY6BrVxg8GFZ93YC23TfyTZ10Dj46qsRj\nO5p0lEe+e4RaU2tR7/V6NJ3RlOdWPFfi/bq/HfiF2rfcQ++anbHGHuPI9Jc58/4cTvTpAqeL1+Uj\nmZlE3tmFYfO3EvzZt1h/WY9Zu5Y6fxvL/Nf388HK1z30LjxgwQKw25nfM4QWM1oR9k4YzZL/zdb2\nIcjIkeDGX//Pr3qeyEOR7B69m80PbebPU38ybf20PMsnJcHOnbBhA2zbBidOXNzNvoNbeOPxjrz6\neHue/epRMhwZRX5rGY4MElITqF2lttt1GgQ0KJddSZ76q70/sAvYA4zLo8wMYC+wFQgrqC5QHVgG\n7AZ+AgLz2G6BmfH770Xq1HH9m1GG/rgqDza9M0EOB9qkQ7MjcvRo3uUe+PsCOVbFIglrIkskLqfT\nKfOi5kmNyTVk7LKxEnMmRpxOp2w5sk0GvfM38X0xWJoNe03ah2XI4MEis2eLpKdfmViS0pNkZr/q\nEntdl2wthM2HNsjbvf0lMbSOyN69Rdv4yZNyoFNT+aVNVTkTezDHy6ceuk++6uAjCSkJRYzegxwO\nkTZt5IuXJ4llbD3xa7lOaoU45d3FUXLju71kZyN/iXttQr6beGP9G9J6Vms5dfa4yMqVIi+/LDve\ne026vN76Qhm7XWTNGpEnnxRp0ULE31+kVSuRrl1F2rUTCQ4WqVlTZELfb+W0r5GDfcIkbcBNkuRv\nky+Gh4szM7NIb+/g6YPS8I2Gharz0baPZOjCofmW2XZ8m8zdMleW7FkiDqejSLGJiBw/fUR+WrJQ\n3pm1VtaudX0d+SGfFoMnkoIF2Ac0BryyTvytLyszAFiS9Xs3YENBdYHJwLNZv48DXstj/7K8dwNZ\nMqqPrHzjCdm2+jM5c/ZUjg/hq69EevQQqVFDZNQokWXLRIp4fBRJuj1d1kSvkQ+3fijTf50uMzfO\nlC///FK2n9ie+8GQkCCyYYPIF1+I/PSTSKz7zVdPST60X45Xs0q/ltNkx478y9rtIuN79ZJfQ2t7\nrPskL4mpiXLPl/fI1XOulqiY7bJjh8iHH4rcf79I3eA0Gd5qvSzq/7x83LuRfB5eXX4aeKc83/l7\nad7UIZFXIG9NnXO/nAnwyfU72nlqpzx3V3VJDvAT+/Q3RVJSsr3ucDokPiVektKTLq50OkVOnBDH\nf6bL6ZBAeefGYIk9ffhiHYfIq6+K3HSTyA3XpsihWkEyf9xAz7+xwlq8WOKaXiVmbE2Z8ul6ERFZ\nvlwkJERk3vsOee+L5+RUFSPrvnwz1+oLti6Qhm80lBMLPxRp3lykY0eRMWPEed11srOOl3z64Sr5\n299cf+R16CDyr3+JbNuW8/+x0ylydP1BSawWKLfe2kOaNXP9YXDq922yqYW/xLQPFdmz55LyTtmf\nsF+W7VuWa5eV3e76zNceWis93utRqI9k48qPZOatjSXmttGSMuafIqtWZQt4x8kdEjI1RB74+gFp\nN6ed3PX5ENnyR7Js2+bab4GcTjnx7aey8ZrGcs4LOVLFSxJ8jWyrHiTvDV+Ub9X8EoMRN5p2+THG\ndAcmiMiArOXxWTucfEmZt4HVIvJ51vJOIAJoklddY8wuoI+InDDG1AEiRaR1LvuXqJf+RvqO/+Gz\n7wA1YuIJSUjnaHUbJ+sHkdKsEdbWVxHYIZzqTa7iTEZdfvq5IV9+XY3oaLjjDrjnHujdG6zWYn0U\nF4hAfDxE/bGP7UtmIDtWUP/IXjoct+Hl9CbNx5dUHxvnfJ2c9j7HWe9MavjXp3ZGTeqetBN8/DDW\nzDRiq7biqFdjApxnaJ4Uxdke/QmeNAafa7t4JtD8OJ3s6tqMxY4q1H56O8OGFVzlwMGzpHeuxW+P\nPMEDr+bd9HdXdDSsWAHHj7s+z1Nxwl77ImJrjqFWdE9qrXmQgOQztKlxgm4Bu+ji2EDI8T+wtG4F\n3bsjzZvzx7n9/PbbN3TdeopaqT680rAXu2/pzs03B9IgsC7Ng5vToXYHrJaiffkHTh/gcHgrOj4+\nicCnxuVaJvZsLBPeuZcR87fS+VAm8a0bkWhzkJ4Yh29CEt5isDgELyz4GhtV0l3/J1e0sLKsfwte\neH4Ztau6ui9SU+G+++DUKdesO4cDlkz5npeiBrF1/u/0u7tD7oFmZMC//4188w0cPYa5LgJGj4Ze\nvbIVO3MGli+HpUth/36w2+Hee123nKlWLf/PIr17b0ZViSNw+MPMvv+pC+v37IHrr4fJk6FD8iRC\nxk5k2pju9B46nnYh7YhOjGb2b7OJ/t9aVvzejoAd+2HWLBJ7DGDVKli+TKizthcPHN7BN2P2M/C+\n6rRokU8gInDttSxokUr6k49zVdqDTJvmuvD17pEHqbW7K8+uSsV531/Y3K89Yw69w4nM0wRbqnKj\nvRGvhfwVx9Y9HF/2P8yBfVjTkkmmCsdDfPmjj+HJmUshJCT/D+PgQexjniXpp0g+aZVO5plXSI+N\nY2jQDzR0RGMZdCvJA/vR5eB4xvWdwAPthzPvwzSeXv0g9mrRhPz0E03qV+Hjj6Fhwzze4/LlnB7/\nFHFH9rDx1huZsOQtps8N5bq+icyZMIT7Z6/BPnoKDac8mWuIxhhExOT6mgcSw51APxF5OGv5PiBc\nRJ68pMx3wKsisj5reTmuVkCTvOoaY06LSPVLtpEgIjmuZzbGyOXvwZmWyrGtv3BiyxqSt0dh2bOX\nwOhjVDuTRpU0BwFpgtNAko8hycfGGZsXiVYfMqv5YW0fSrWIPoTeNISQFmF5vm8R+PNP13+iP/5w\nEndwL9UTVlM3bS3NMn7n6vTDtEtIZk/1mhyu2wVns1vw7tATh82H9LgkMhLOkRZ3jpRT50g+c5hU\n7x2cCdnPrrq7iK7lR/XgG+kQ3Iv2IWFUk4ac3GbD9+MPGHzoPyRVb8yRIU/T5Jm+BNfLwMfqg7+X\nf/aTm93u6ny1210zQwqZ9Y5OGEPMvJn854aDfDrf/VkYi8dOp/kH/yBz4590aNbKvUqZma4bXW3Z\nghw/weFtp9m54TQkH6dBtRNUc57GL/0MVdKS8bELmVUCsFavhSWkFta6IZg6tV0DSj16uAY/qlTJ\ntnkRYX/CPnZ+O482MxdiPZDAtI5dOHh/VQ6m7SQuJY47Wt/B0z2epnXNHH975OtfEyIYPe9/BB88\nzoEYL5Yuhb17XdfSdO/uOu9ara4YluxdwtY/V1Fl+24aWYNp0qgDbdr2wc8/AIeBoykn2HfmIIft\n8aQG+NGj4TV0qHPxRC8CQ4e6ksHHH4OPz8U4Ph3QjeDfT7L8LweY9IrBz+9ine1rEvB/4C4Ox1dh\nYvpzHLfUo7/5ifGOSfwZ3Juvr5lKon89DhyA7dtdMQ8YAFdd5con773nWv/dd9CyZe6fg3P9Bo7e\ndCs9xjQieuKmHIl2xw7o2xdmzYJbfb/Dfv9fWN4xkMUNk2lqgrntcBWuiorh7PC/M6/mOBYu8WP7\ndrj2WrjhBmjU/Tfs427g3vZ/xcyZk/+X8sUXyOTXqHHXAf4cvZO61eoCru9l9myYv/A4jQcOZPjO\nXdy5w079Mw4sFivY7Ryu4c+fQdXZdmgkPl3aEX5fS9pfUxVrWjLz3ptMraWbuS3uOH6j/grPPpvz\nrH32LLz6KvLuu3xS62mWdXycz9vWIfX5VBLiLbzwAmz68hAzrv+G0L2zCd51kIxWPfj5WCsSbTW5\n9rYaRNoWs89yhqapM/jg6wZ8vqIGDZv7QFoaHDgAa9Ygn3xCXNwhJvZ0MGj81zxxRx8mToS//MUV\nRoYjg77/DGfhjIPUXvGD64O8TFlMDCuAZylcYogXkRx3rzLGyLBhwwgNDQUgKCiIsLAwIiIiAIiM\njATItixOJ926hpF08ggrflpKxtkE6noHs3HFcU7sXkE7DjHkeBI2h7DEz8a5aj60rRtApkXYlpCK\npDro5DQEpNvZm5xOYLqTLj5woI4/q4KqY2najP73DKFF3yH8+sefOfaf3/Lq1as5fOYwSfWS2Bi7\nkXVr1hGXEkdagzRqValF2q5UOu05x8v7Mml0BibX9SI2SKhT007LeCtpR6D+WeHOVCHVy8Iah+Bv\nd9K0egD7mrZgd/i1tLtjENf17ZtnPI7167j6tZcY2Pkl/vVsD/z83I8/cvVqdgx9gGNXWZmwYh+/\nrPkl7/LHjxM5ejSsWEFEaCjJ7TrzwoZ9xPocolqrEziCqnA6NQD/6iH0uL4PXa66gdRjFiwWq/vx\n5LLs3PI7rWZ9hfNwDN+0701IzyBSmp/jrePfUce3K6O6Pcrg/oML3N6GmF/5o0svatz9D749+xpL\nl0KnTpE0aQJVq0awciUcOBDJtdfCiBERdO4M0dGRGFO4eM8vT5kC8+ZFMmMG9OuX/fVW7ZtyJKwp\n39puY2b8aPr2jcAYiP7lY8bEPUejzkM48PIdfLD+Bbaf+oOqLQII8+3IdZ+cpM2WbTTqOQxn92s4\nFL+OM3u2cM3hA1Q5cZpfMhzsCW+Kufk9pr7RhwkTIrnqqpzx1Rr9OnPr/UzLB6dyVUirXOPfuhWu\nuy6SJ56Al0Y2wfz3UyIXL8YZEEhwk8H8e989LN30BxERMHp0BH36wK+/uur36dOHdq/U5z+vJ2N9\nbQoRjzyS++e1bBkMG0bgtH/wcMpnTG0xNcfnmZIChw5FsGQJrF0bSeo5O77mGsTHl6ZX/cDOtvfy\nv5ejadWwZrbtP7X0KRL/yGTTG335tN2vhG2ZR2R4OPTsScT110NUFJETJkDXrqxt8j6rd9fjueci\nufur29g1bRd1qtYhMjKSfftg4cIIlje8luDVNxB+rhrjbgigZ5s41mz9HWdiIjUTdmHi4zmQUAXv\nc8n093ZifHyIrFWLc03qsy4siY1tqvF4nb8zc0owrVtH8Pbb2T+PPaf289DtV/P8/mBuOrCPyI0b\nmT9/PgChoaG89NJLVzQxdAcmikj/rGV3upJ2AX1wJYZc657vbrqkK2m1iLTJZf85WgxF5XS6upVE\nYO67Tkz6YeL2/0FS9B5Sjp1g9zbDuk2++AdXIaJ/Fdr0qI4tpDa1G7clMLiuR2LIS7o9nRPJJ7BZ\nbFT1rkoVryrIjr0ce/c7UnYdIj3Fztm6oSTWaUBCvbok1qlBpp+dTJJJPpuIY+dWArf9zE3b1xGU\nDjE3j6DXGxMxdetk39GyZZy7604GXN+EOZO20u7qwk9cOxe1k7RrO/DYuPv4fMI8jLns2HM6YcYM\nmDQJhg8n/aGRPPXrCt7d8RoNfK5i0h0PMqDVDW7fk6aoYr6NIuq5L7DvjaaF/xHqOg8TmBrLsWqG\n+D4DaTt5Oj6tQnOtKyL865FWjFiaTnjaQR56xMLYsRAQkL3cgQOwaBGsWQNRUZCcDJ06QceOrn+7\ndYNmzQqOdelSGDnSNbsu164FYNqisdz3tzn4jniBNXWG0HDb91y9+BWsk17i1bbxzNk8h3HXjmNQ\nq0FkODL4OfpnFu5cyIlt63nwQBBXR6dw0HqW5GaNCOx9E00730Cww5vkKZOo9fNmPnr2X8x443kW\nL3a1hs777ZM9hI7qwItzR/LW/flfOLR588VuqS5ZE7Z++gnatIGHH4a77gL/PB51cNcXd/GPXcF0\nX7LNdQGdJZdj8803YeVKxjzRkkCfQCZETCjws83MdP2f9/Z2Ld+36D661uvK37v/Pcf+72l7D539\n7qZXL3h7Uhy3Jnzo+nKSkqBpUxgzhnnbujB5smuWVHAwdH63M2/f8jZd63e9sK298Xvp+X5PNtxz\nhCaNvHLEJCLMjZrLuBXjCI4fQMDp3jz4kPDzodWsOriK8T3H81T3p3h7jo1581wfh69vzvfWaOzt\nLF6yj7BnnoJR2WcN5tdi8MTgs5WLA8jeuAaQ21xW5mYuDj535+Lgc551cQ0+jxM3Bp89KSVF5OGH\nRWrXFnniCZFJk0SGD3ctDx4s8uuvHt1dicvMcMibL7wi77X3l9Pe3nK0xwBxTp0mMneuyKhRkhRc\nWyKGBsj7328v1n4OPzpe1tXxl1FvP5/92oLYWJEbbhC55hpJ37Fd3lw3W6q8WF98R9wqb327pZjv\nrmgSE10TYN59V2TcWIcMG/iBTO4WICe9/WRmp3fkg/edcvx49jpvff2RxFa1yZDGa2XTJvf3dfy4\nyI8/irzyisidd7oGZvv3F4mMzHvMfs8eV7k1a/Lfdmpmqtz4fCOJveka1wE7cKCkbFovQ74cIuFz\nw/Ocfx+fEi+/xvwqqw+ulsTU3K81iH3zZTkRYJWnXnxQatR0yowZIqdOueZGLArsI5OvryZn09y7\njsXhEFlIVzxAAAARO0lEQVS8WGTWLJF580RiYtyqJpN/mSxPLXnSNYvk3XdzFkhIcE1H+vNPaf9W\ne/k1pmj/WVceWCnt5rTLcU1Mt7ndZP1h16B6VJRIrVoiq1dnr/vjj67vavfui+sG/XeQLNqRfSD4\nhVUvyFM/PlVgLAkpCfL6L/+ROg8+LB1eGCnzoubJyXMnRcQ1gaZWrfwnvD333lKZ0aWWOKdMyfEa\nV3JWklyccrob13TU8VnrHgEevqTMrKwksA3olF/drPXBwIqs15YBQXnsu8APtyh+/911Idezz4q8\n847Irl1XZDelJv5cogx+9REZdks1ebtpV/mp+QCZ3mGIhDxRU8b9d37xd+BwyOGI2+SbJgFyyz+v\nk93RUeL8+GOR2rXlzHPPyPS10yTk1UZS9ZEBcsMDmyShDMy2vFRccpwMe6m77K5dRZaFDpE6gSnS\nqpVIz54iTTvvlHnt/GRtpwGSlla8/aSmus5xLVqIdOsmsmhR9mmGUVGuizTfe8+97a06sEpqTK4h\n/17zb5m7Za50eKuDPPD1A5KamVq8QEXk7Htz5HiQlzw65TbpEZEoNWuKPNxysuytYZXl2xcXe/sF\nWX1wtWtW0O+/u86+J09mL/DMMyIPPyzn0s+J3yQ/Scss2pfjcDokdHqobI7dnG19/dfry+HEi7PD\nVq505aFZs0Ti4kTef98V1rp12bf3xA9PyNR1U3NsP+polNsxnTwp0rq1yIgRIvv2iXzwgSsprF2b\nf734BIdMCg+UA39/OMdrVzwxlObPlUoMlcXOk7vknvmjpdmr3eTqKX3lhx2rPbfxtDQ5OuRxOeXt\nL0leRtaGeku/h4PE/6UgqfPoA1I3fJ0sXHjFZ7cWWaYjU15eMl4WhvlIfMNa8tukWTJ32gz5rEtV\nOd6+mcjp0x7bl90u8uWXIp07izRs6Jp6O3Cga07+woWF29aBhANy1xd3ydCFQ+Xz7Z979Grw1Fn/\nkTNBfvKPu4Nk8oNt5Fg1I2s/edVj28/P2bSz4v+Kv+tK73/+U6RXL7mQmX/80dVKOnZMfo7+WcLn\nhhdrXy9FviSPff/YheXDiYcleHKwZDqyz43dvVuke3eRgACR8HCRP/7Iua2vdnwlNyy44cLyz9E/\ny9Vzri7095KU5EoMdeqIXHeduN1Sfa1LuPxw84051ueXGIo9xlDaPDnGoK6MpPgMFv43gwUrT2FP\n86VRzRDuuM3KwIHZZ9aUVVuObiHy7fH0XhBJDbsPpn17mny5PMfsJ08QgV27YO1aqF7dNY26tvsX\n2paMqCiSRj9MalUfZORIat9bcle8Xz3nahbcvoBOtcNgyBA4d841jWr6dNfN0Xr0YOq6qcScjWHG\ngBlF3s/hM4fp+E5Hjjx9BD8vP2Zvms3G2I0suH1Bobd1LuMcdV+vy5GnjxDoG8iDix+kVY1WjL12\nbJHjK4xp1/ehZaYwaO2abOvzG2OoMA/qUWVXtRrejBjtzYjRVUs7lCLpXK8znf+1HP515fdljGsg\ntk2OaRZlSKdOVFu/mQIua7giwuuHsyl2E53qdnLdfuOtt1xP6PrsM9d0ZWDT0U0MbjW4WPtpFNiI\nLvW68M2ubxjabijf7P6GRzs/WqRtVfWuSu/GvVm6bymDWg1i0c5FbH9se7HiK4xMLx+8ks8Uqk6F\nvleSUqpiOZ8YANf0pWeecbUWsqZfA2w8spFu9bsVe18jw0Yyc9NMTiafZOORjfRr3q/I2xrUchCL\ndi3ihdUvEF4/nHrV6hU7PnfZbb54ZRbuyWWaGJRS5Ua2xJCLY0nHSM5Mpnlw82Lv6/Y2t1O3Wl3C\n3g6jd+PeVPUueov31la38sWfXxB1LIr5t80vdmyFkenti1dG4RKDdiUppcqNdiHtOJh4kKT0JKr5\n5OzM2hS7ifD64TmvnSkCb6s3X979JVPXTaVj3Y7F2la9avX4ddSvdK3Xtci3Xykqh5cfXpmFu6us\nJgalVLnhZfWiQ+0ObDm2hYjQiByvb4zdSHi9cI/tz2IsjOuZ+z2wCqt7g+4FF7oC7N6+eBfyEbza\nlaSUKlfy607aGLuRbg2KP75QkTi8/fEuZItBE4NSqlzJKzE4xcnmo5sJr++5FkNF4PT2x9tuL1Qd\nTQxKqXIlr8SwK24XtfxrXfF7bJU3Tm8/fLQrSSlVkTWr3oxzGec4lnQs2/qNRzZqayEXTu8q+GRq\ni0EpVYEZY+jRsAfrYtZlW78pdpNHrl+oaMS3Kj52R6HqaGJQSpU7/Zv154e9P2RbpwPPefCpgo+O\nMSilKrpbWt7CD3t/wClOAI4mHSU6MZqOdYp3vUFFJD5V8M10FqqOJgalVLnTtHpTqvtVJ+pYFACf\n/u9Tbm99Oz62cnBXxhJm9fbHIK4nErlJE4NSqly6pcUtfL/newA+/uNj7u9wfylHVDZ523xIs1kg\nNdXtOnrls1KqXLqlxS08/P3DtKzRkvjUeHo37l3aIZVJ3lZvUm0WqqSm5nz2bB60xaCUKpf6hPbh\nb13+xtjlYxnWYRgWo6ez3HhbvUjzMtpiUEpVfBZjYUyPMfy92989ctO8isrH5k2qTRODUqoSKem7\nlZY33kVIDNr2UkqpCszb6kWqF5oYlFJKubi6kkQTg1JKKRdNDEoppbLxsXmT4qWJQSmlVBZfLy9S\nNTEopZQ6z8fmTaqXUxODUkopFx8vGylegjMlxe06mhiUUqoCs9kMaVYLzpRzbtfRxKCUUhWY1Qqp\nNiuOZE0MSimlAJsN0qw2nCnJbtfRxKCUUhWY1ZqVGFJ1jEEppRRZXUlWG6ItBqWUUnCxK0m0xaCU\nUgqyWgwWb0jR6xiUUkpxfozBSy9wU0op5WKzQYrFG9JKKDEYY6obY5YZY3YbY34yxgTmUa6/MWaX\nMWaPMWZcQfWNMTcYYzYbY7YZY34zxlxXnDiVUqqysloh1XhDaprbdYrbYhgPrBCRVsAq4LnLCxhj\nLMAsoB/QFhhqjGldQP1TwEAR6QAMBz4qZpxKKVUpucYYfLCklVxiGAx8mPX7h8BtuZQJB/aKyCER\nyQQ+y6qXZ30R2SYix7N+/xPwNcZ4FTNWpZSqdGw2SDU+mNR0t+sUNzGEiMgJgKwTeUguZeoDMZcs\nH8laB1C7oPrGmLuAqKykopRSqhCsVkjFF0ua+4nBVlABY8xyoPalqwAB/i+X4uL2nnOXrb4xpi3w\nKnBjMberlFKV0oXEkO7BxCAieZ6UjTEnjDG1ReSEMaYOcDKXYrFAo0uWG2StAzieV31jTANgEXC/\niETnF+Pw4cMJDQ0FICgoiLCwMCIiIgCIjIwE0GVd1mVdrpTLycmQgi9rU9L4aPhwgAvny7wYkaL/\nkW+MmQwkiMjkrNlG1UVk/GVlrMBuoC9wDNgEDBWRnXnVN8YEAZHARBH5poAYpDjvQSmlKrLkZAgc\n/lcyv/ovxm4Hi2sEwRiDiJjc6hR3jGEycKMx5vyJ/7WsHdY1xnwPICIOYDSwDPgT+ExEduZXH3gc\naAa8aIz53RgTZYypWcxYlVKq0rFaQRw+OK0WyHRvqLZYLYayQFsMSimVN7sdvG5/lIyV8/E6cQqq\nVQOubItBKaVUGWa1Ag4vHDb3WwyaGJRSqgIzBnB64yhEV5ImBqWUquAsTm8cNitkZLhX/grHo5RS\nqpQZvLTFoJRS6iKreGPXxKCUUuo8i3iTqYPPSimlzrOIN3ar0TEGpZRSLha8sFuMthiUUkq5WMWb\nTCuaGJRSSrlY8CZTu5KUUkqdZ8WbTAvaYlBKKeVixYsMTQxKKaXOs2pXklJKqUvZjDcZFtEWg1JK\nKRftSlJKKZWN1XiTri0GpZRS513oStIxBqWUUqBjDEoppS5jNV6kWZ2aGJRSSrl4aVeSUkqpS9mM\nN+nGoS0GpZRSLl5WL9It2pWklFIqi5fxJt3i0K4kpZRSLjaLt7YYlFJKXeRtyWoxaGJQSikFYLN4\nkaqJQSml1Hk2q5VMq+BMT3ervCYGpZSq4LxsBrvFimTq4LNSSinAagW7zYJkaItBKaUUWYnBakV0\nuqpSSikAmw0cFosmBqWUUi7nWwzoGINSSinQriSllFKXsdnAbrFoi0EppZSL1QoOq00vcFNKKeVi\ntUKm1YqURGIwxlQ3xiwzxuw2xvxkjAnMo1x/Y8wuY8weY8w4d+sbYxoZY5KMMWOKE6dSSlVmrllJ\nNkwJjTGMB1aISCtgFfDc5QWMMRZgFtAPaAsMNca0drP+68APxYxRKaUqNasV7BZriXUlDQY+zPr9\nQ+C2XMqEA3tF5JCIZAKfZdXLt74xZjBwAPizmDEqpVSl5upKskGm3a3yxU0MISJyAkBEjgMhuZSp\nD8Rcsnwkax1A7cvq1wYwxlQFngVeAkwxY1RKqUqtsF1JtoIKGGOWk3XCPr8KEOD/cikubu01b86s\nfycAb4pIijHm/D7zNHz4cEJDQwEICgoiLCyMiIgIACIjIwF0WZd1WZcr7bLVGkHikTRGxSdgveR8\nmRcjUvRzuTFmJxAhIieMMXWA1SLS5rIy3YGJItI/a3k8ICIyOa/6xpg1QIOsTVQHHMCLIjInlxik\nOO9BKaUqusmT4aOTHfl9fjRe8acBMMYgIrn+0V3crqTFwPCs34cB3+ZS5jeguTGmsTHGG7g3q16e\n9UWkt4g0FZGmwHTg37klBaWUUgVzDT7bMCU0+DwZuNEYsxvoC7wGYIypa4z5HkBEHMBoYBmugeTP\nRGRnfvWVUkp5jtUKduOFcXPwucAxhvyISAJwQy7rjwEDL1leCrRyt/5lZV4qToxKKVXZ2Wxgt9ow\n9pKZlaSUUqqMs1rBjhcWuwPcGJPVxKCUUhWc1Qpgw2lz7yI3TQxKKVXB2WxgxIrY3LuRniYGpZSq\n4KxWwGnD6aWJQSmlFFmJQayuriQ3rn7WxKCUUhWczQbGqWMMSimlslitIE6rJgallFIuVisYseHw\n0q4kpZRSuLqSxGnFadUWg1JKKc7PSrLisFk0MSillLo4XdXhpS0GpZRSXNqVZNExBqWUUlktBodN\nu5JU5XX+qVVKnVfZj4nz01XtOl1VVVaV/SSgcqrsx4TNhmuMQbuSyoYrfUAWd/tFqV+YOu6Uza9M\nUV8ry/SYKLisHhOetW1bJOKw4rAZbTGUBXoSKLisngTK1vb1mCh5VzrurVsjEacNu9W9MQYjbjy0\noSwzxpTvN6CUUqVERExu68t9YlBKKeVZ2pWklFIqG00MSimlstHEoJRSKhtNDEoppbLRxKCUUiob\nW2kHcCUYY/yBOUA68LOIfFrKIalSZoxpAjwPBIjIkNKOR5U+Y8xg4BagGvC+iCwv5ZDKjAo5XdUY\ncx9wWkSWGGM+E5F7SzsmVTYYY77QxKAuZYwJAqaKyEOlHUtZUS66kowx84wxJ4wxf1y2vr8xZpcx\nZo8xZtwlLzUAYrJ+d5RYoKrEFOGYUBVcMY6J/wNml0yU5UO5SAzAB0C/S1cYYyzArKz1bYGhxpjW\nWS/H4EoOALle2afKvcIeExeKlUx4qhQU+pgwxrwG/CAiW0sy0LKuXCQGEfkFOH3Z6nBgr4gcEpFM\n4DNgcNZrXwN3GWNmA9+VXKSqpBT2mDDGBBtj3gLCtCVRMRXhmHgC6IvrXPFwiQZbxpXnwef6XOwu\nAjiC6yBARFKAkaURlCpV+R0TCcDfSiMoVaryOyZmAjNLI6iyrly0GJRSSpWc8pwYYoFGlyw3yFqn\nKi89JtTl9JgogvKUGAzZBw5/A5obYxobY7yBe4HFpRKZKi16TKjL6THhAeUiMRhjPgXWAy2NMYeN\nMSNExAE8ASwD/gQ+E5GdpRmnKjl6TKjL6THhORXyAjellFJFVy5aDEoppUqOJgallFLZaGJQSimV\njSYGpZRS2WhiUEoplY0mBqWUUtloYlBKKZWNJgallFLZ/D88DWAueZVTzAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# compare calibration files from points 1,2 and 3\n", "# with the optimized from above. To estimate the precision\n", "# at each point the curve with the biggest difference does\n", "# not have to be taken in to account\n", "\n", "for sensor in run.sensors:\n", " r0, t0 = read_curve(sensor.outputpath, sensor.outputkind)\n", " plt.figure()\n", " dif = [0,0,0]\n", " for j in range(3):\n", " rref, rtst = read_curve(sensor.caldat_file[j], 'zdat')\n", " rc, tc = make_calib(run.rref, run.tref, rref, rtst, run.t_points)\n", " dif[j] = compare_calib(r0, t0, rc, tc)\n", " plt.plot(t0, dif[j], '-')\n", " #n = len(dif[0])\n", " #dd = np.zeros(n)\n", " #for i in range(n):\n", " # # determine second biggest absolute value\n", " # dd[i] = sorted([abs(dif[j][i]) for j in range(3)])[1]\n", " #plt.plot(t0, dd, '.')\n", " plt.xscale('log')\n", " # plt.yscale('symlog', linthreshy=0.001)\n", " plt.grid(True, axis='y')\n", " plt.axis([min(t0),max(t0),-0.005,0.005])\n", " plt.legend(['dif1','dif2','dif3','est'])\n", " plt.show()\n", "\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.5" } }, "nbformat": 4, "nbformat_minor": 2 }