Skip to content

Instantly share code, notes, and snippets.

@flamingbear

flamingbear/final smmr mask comparison.ipynb

Last active August 29, 2015 14:10
Show Gist options
  • Save flamingbear/3be38b15f4365127f81e to your computer and use it in GitHub Desktop.
Save flamingbear/3be38b15f4365127f81e to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
final smmr mask comparison.ipynb
{
"metadata": {
"name": "",
"signature": "sha256:f8bff85dae2e576df8a963e2bdd11776de3f43e9f3e57c78a9289abef7ebf115"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 2,
"source": [
"Comparison between extents using SMMR mask and using default Data"
]
},
{
"cell_type": "markdown",
"source": [
"In looking to move to the F17 pole hole, I found that the daily extents are currently computed after slapping on a SMMR (N07) mask. (This is the largest mask for SMMR data).\n",
"\n",
"This is coded into the `generate_extent_series` function and **only** affects daily csv files. The code chooses to use the SMMR/N07 Mask and mentions that any areas for F13 would be bad, but extents will be fine.\n",
"\n",
"What I'm seeing is that there are a number of days where the openwater grid cells creep into the SMMR mask effectively raising the extent value for those days.\n",
"\n",
"For the final data there are 132 places where open water creeps into the SMMR mask.\n"
]
},
{
"cell_type": "markdown",
"source": [
"Below is an image showing the problem area for the worst case scenario: **2009-09-22** (0.056 M km^2)\n",
"\n",
" white = ice\n",
" blue = open water\n",
" brown = land\n",
" light_gray = ssmr mask\n",
" dark_gray = ssmi mask\n",
" red = open water in the smmr mask"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import datetime as dt\n",
"import pandas as pd\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import matplotlib as mpl\n",
"import os"
],
"language": "python",
"outputs": [],
"prompt_number": 41
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pylab.rcParams['figure.figsize'] = (8.0, 8.0)"
],
"language": "python",
"outputs": [],
"prompt_number": 42
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" sample_filename = '/projects/DATASETS/nsidc0051_gsfc_nasateam_seaice/final-gsfc/north/daily/2013/nt_20130922_f17_v01_n.bin'\n",
" smmr_mask_file = '/projects/NRTSI-G/30yr/dev/nrtsig/ancillary/gsfc_pole_hole.N07'\n",
" ssmi_mask_file = '/projects/NRTSI-G/30yr/dev/nrtsig/ancillary/gsfc_pole_hole.N13'\n"
],
"language": "python",
"outputs": [],
"prompt_number": 43
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def read_icefile(filename):\n",
" with open(filename, 'rb') as fp:\n",
" _ = fp.read(300)\n",
" data = np.fromfile(fp, dtype=np.uint8).reshape((448, 304))\n",
" return data\n"
],
"language": "python",
"outputs": [],
"prompt_number": 44
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" data = read_icefile(sample_filename)\n",
" smmr_mask = read_icefile(smmr_mask_file)\n",
" ssmi_mask = read_icefile(ssmi_mask_file)\n",
"\n"
],
"language": "python",
"outputs": [],
"prompt_number": 45
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" img = np.copy(data)\n",
" img[ data < 15 * 2.5 ] = blue\n",
" img[ data > 251] = brown\n",
" img[ (data > 15 * 2.5 ) & (data <= 250)] = white\n",
" img[smmr_mask != 0] = light_gray\n",
" img[ssmi_mask != 0 ] = dark_gray\n",
" img[(smmr_mask != 0) & (data < (15 * 2.5))] = red\n",
" \n",
" imshow(img, cmap=categories, norm=norm)\n",
"\n"
],
"language": "python",
"outputs": [
{
"output_type": "pyout",
"prompt_number": 46,
"text": [
"<matplotlib.image.AxesImage at 0x11058f810>"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAVgAAAHfCAYAAAAV2cAgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvW2MHtl15/c/zRd5N/JCEIyZYXdPIiM7tjzxAvIXOUjW\npLNYayg50DhfrLCFhU0OSAFWcwZZJLGGXCT+Qo5tIBtL7HGg5oqtibOkJewishYrayQ79ozzwVaM\nSLak8UgaIAOI3RQZLOKF5EVCcvrkQ9V96tatW+/3Vt176/yAJp+nnnq59VQ9/zr33HPOJWaGIAiC\n4J61uRsgCIKQKiKwgiAInhCBFQRB8IQIrCAIgidEYAVBEDwhAisIguAJLwJLRKeJ6HUi+g4R/aqP\nYwiCIIQOuY6DJaIjAL4F4B8C2AfwfwA4w8x/5fRAgiAIgePDgn0vgDeY+U1mfgDgdwE87eE4giAI\nQeNDYDcAfFd7fztfJgiCsCiOethnq8+BiCQ/VxCEKGFm6rquD4HdB/C49v5xZFZsmWP/IXD8iez1\n8SdAb3vCQ1P8wt//AuiHPzB3M0bj+zwIwO6pK6i7KxnA+Vcujz5OCNej7VyB7HwvvHK51hJxfR4E\n4PqpK87215Xf+8YP8PRPvn3y47rk9Xv38Xvf+AHe/chxAMDnv/k3vbb3IbB/DuAJInoXgAMAHwJw\nprLW8Sdm/zEI05AJyiVcP3XV+nkhAITzr1yasmle6GLeDO3CKQFX77p8X+UHGNdeB6HKux85jm89\ncnz1oOgrsM59sMz8EMA2gJcBvAbgMxJBIHSRnQuRi2tZ/JrXu37q6ur/vvtX32Sf74tXf2k8xGLB\nhwULZv59AL/fuNLx+FwCFVI4B2CS88isqHorVq0zihmuhymqXZ1zBK4XY+M8zGMoa3T499XZhTiK\nH8+71bEz5jy8CGwXYvS5mqRwDsCU5+H3hz3l9dBFb+hZqe2uG13+tbc9URFf/Ri7p65GYYW+OxGB\nHXMeswmssDy6WLEx0GUgq+/+Mngltk0DgrG7UpaECKwwKcoHqIusGlUPHRdWa5djtNHFNVD1B2ff\ne1c/seAGEVhhcjKRLUa1d09dHe9/nYhpvJdV+jyE7BZ2u3UsuEcEVpicworKrKoYrNc5aYub1Wly\nX4iwTo8IrDALBIBzyYjFep2DNnG1dflFSMNBBFaYHDNzy7RohYLsu7FHDbgebBPcIwIrzIZpeYU8\nOj7n4BDl2VdKZKcYbBPcIAIrzIKZPqsyjUIkBEuRtEGq7L0QA84Lbnc6KBHTY5+Y/LhCeGRCkVlo\nnC8JyU0QgrgK4fDMZ+7OXk1LEDoXJdFFNRPZMOxY6YYLLhCBFbxRzVAiqx9x6jY1DaiJsAouEYEV\nJqFwBdj9iK5qwtYdu1qQpV70RVgFV4jACs5psk7bClH7aktd4D2vnBQiroJ7RGAFp4QkVF3aolf7\nD6HNQlqIwApO6Ou79FngpW9bRFgFX4jACqMYMyg0tUtAEKbGx7TdwkLQxayvoPWdLqXvvgUhBERg\nhcGMKf7so3C01DoVQkMEVhjFmEn0fLkIBCEURGCFwRTdfFrFsc6V5irWqxAiUotAGEUWS1r8r5Zl\ncK2f1WVigQxsCVPRtxaBWLDCKNj4X73O/si7VTtYXIlw7tYdnL114KFVgpAhAiv0pk8EQCG0fkV2\n2IYEErtX8Ii4CIRB6C6B7ttwqf7rWBfBaNcA5VvO8BsQ4kRcBIIXTKt1iCTpluzYOFgnfldmEVfB\nKyKwQitqXigXcatKZMfEwcqglhALIrBCR9yVwmZQ52moTURchZgQgRUaKaxXt4VZxoi1iKsQCyKw\nQiOqK+/SUznU/yrJBEJsiMAKFVwMaLUzbK9ivQoxIQIrVPBRiGUsYr0KMSICK1jx5RIYI5RivQqx\nIQIrAPBbnzXbfzbJoEQACEtCBFZwGudqI8vaujSoMDcg7gEhXkRghRx3ca72vQ+rRyBxr0LMiMAK\nXicgLB9nWKqsiKsQKyKwAgBfoVg2dLmUOgBC2ojACpOi/LFdrGbxvQqxIwIrzEBmxTbZr+J7FVJA\nBFYIDhFXIRVEYIVZaBvkEnEVUkAEVpgNfYYDQUgREVhhUvSBqxBrHgiCS0RghcmRQC1hKYjACoIg\neOLo3A0Q0qZvLKvEvgopIRas4BXlZ1WFXtqiA8QvK6SECKwwAXZZrQ/VkiAtIQ1EYAXvqPTY+k8F\nIU1EYIWJEKtUWB4isMIkNFuxgpAmIrDChIgVKywLEVhBEARPiMAKkyFuAmFpiMAKE9PsJpBEAyEl\nRGCFSTGtWFssrHhqhVQQgRVmgIx3WdlCsV6F1BCBFSbH7ovNEg7EehVSQgQ2YfpMjT09IqVC+ojA\nJkrR3Y4nFVWiDITUEIFNnLapscNDLFshHURgEyce+1UQ0kMENkEy98DVCK1XQUgLEdgEUUWrQ7Ze\nu8TDCkLsiMAmhhKqkMW1oPC36g+F869cbhzs6rKOIISACGyCqMD9GFBieeGVy6uHAgNgkFVAMyHO\nXB+7+YNEhFYIFRHYRIlBZJVYMuoG4+wRBWr9Yu4uiTwQwkRmlU2IaqppHI6C66eu1lqr5xsG6hgk\nA3lC0IgFmxhkvA7dis2ofxCw5a/bloIwP2LBJk/YEqT7VPugQtHOawNj+l7jeLAIqSMWbOLEYMUO\nfwRw6VVh4YpPVggDEdgFEMOAlyuk5KEQEiKwiyFsV0FfbK4FJa5ivwqhIAK7EGJwFfRFf2SIuAoh\nIgK7KNKyYhUirkKoiMAuiBStWBFXIWREYBOh6+BOSgNeIq5C6EgcbOTowtpdaOJ3FYi4TgQRzt08\nqCxmMPbOrM/QoLgQgY0clZPfxypVroLYi6SIuDoiF9GVaFpE9cZWWUwJhLO3DkRkWyDm6a0ZImJ6\n7BOTHzdl+nb92/L8Q0asV4coMaX829T0oCSqNp0gsi9PmGc+cxfM3PnWEws2GfrJTaxWrIirQ3Rx\nZa5Yqa3iuTBxHYII7KKJ6wci4uoAs/uvi6sIpnMkimDBxBS2JeLqAN1izV0CN86cEHH1iFiwCTAu\n/z6eH5aIqwM0X6sIq39EYBMhZfGRAi4OUNar7msVcfWOCGzELEF4xDXgAM3vKlbrtIgPNlJ04Rkj\nPjH4YUVcRyDiOisisJFSnvTPxd7CYwkW+lSIuM6DCGzE1E1tnQLiGnCEDGbNighs5KQusoIDRFxn\no1VgiegGEd0loq9ry95JRF8mom8T0ZeI6B3aZ88T0XeI6HUiep+vhgs6w6Vo6KSDPhHXgJAKXSzY\nPQCnjWUfA/BlZv4xAH+YvwcRPQngQwCezLf5bSISKzlwQrJvxDUgpERrmBYz/wkRvctY/EEAp/LX\nLwH4Y2Qi+zSAW8z8AMCbRPQGgPcC+FNH7RUSZVjZRUEIm6FxsI8y89389V0Aj+av11EW09sANgYe\nY5EUQtPuW42tK521NyswY2u7CKuQGqMTDZiZiaiplxlSDzR4uvpEY+lKm0Kaxd2KpSosg6ECe5eI\nHmPm7xHRCQD38uX7AB7X1tvMl1Xg73+heHP8CdDbnhjYlPnQLTKXdHkiDSm0bd9H/wGusmjaLe2m\nLr8IqxALr9+7j2/duz94+04Ft3Mf7L9i5r+Xv/9NAP+GmX+DiD4G4B3M/LF8kOsmMr/rBoA/APB3\n2ThISgW3CfOa6JlY9Z9nS4lrn7ab1qgSZ3MfsVjXgtAX5wW3iegWsgGtHyGi7wL47wD8OoDPEtEz\nAN4E8IsAwMyvEdFnAbwG4CGAXzHFNTXmPrns+JkVOURku6As9Qu571TR5P8VcRUEmTImGXTfZhtd\nrVezm88Q4RSWjUwZI3Sij/Vqey0IQjuSBJAY2WSG9YNufQa21MSI2V+a6biC4BOxYBNBnyVWiWyd\nT7aPU0itS2K/CkJvRGAToiycIoiCMDfiIhBaiS1jTBBCQQRW6ITYw4LQHxHYBHFpcYr1KgjDEYFN\nBHNuLZvFOSY1VhCE/ojAJgWv/q0Lq5o780wQloQIbCJUrdPxdqe4BwRhHCKwCdFWM1LcA4IwLUkL\nrOmXXDq6AMt3Iwj+SVpgVc3UpaBEs2vXnhrKHIp7QBDGk6TA6tbZ0gZ1CFypxWq6B8pVsqoiK/Vc\nBcENSQrskixXW1e/LkTLto75XYm4CoI7kq1FsCTLNbNCh3XnZfJBQfBHsgK7NOpEsck90LatIAjj\nSNJFIBSo6V6K19L9F4SpEIGNmO4j/UWG11J804IQAiKwEVA3kNXVGi1vL/arIExFMj5Y1RVOdWoT\n20BWP6lc0rCfIIRBMgK7hO6v2J6CEBdJuQjERmtGsrMEYVqSElihHokgEITpEYFdECKugjAtIrCB\nI916QYgXEdgIEMtTEOJEBFYQBMETIrCCIAieEIEVBEHwhAisIAiCJ0RgA0YiCAQhbkRgA0ciCAQh\nXkRgA0WsV0GIHxHYgBHrVRDiRgRWEATBE9EJrK34dGqIe0AQ0iA6gc1IvzChuAcEIX4iFdh0EetV\nENIhKoFV08Lo01CnhNRsFYS0iGrKGDUtTGoOAt1qFXEVhHSI0oJNEYKIqyCkRlQCC6jZVdMSWQaS\nnQ1XEJZMVC4CIN3ZYxmE8yvfcnoPEUFYIlFZsJm4Xk7OB6vg1R+JResB6SkIUxOlBbsEyhZtsVQs\n23Goh5d8j8IURGXBLg3W/oB0B/imQPV+stfSQxCmQQQ2AiQ+1o7Z5Vfvs/8v53/lzwvk2xT8E52L\nYGmEKK66NaiWdO1yF9u6cXeYrpSqn97mahGEaYjOgl1CsReTEMW17L5ob6GyKgvx63BWRDh766DT\nvnVXiumnZ8syQZiCSC1Y+blMiWmx6t8+EfDg9jaObdRbpLooEwFv7V/MP9jGsc0dPNj/KPbOrFs2\nZOxtbdS2Sz1sxZ8qhEqkAitMQXlgqBDToxs7xToMHH/8RTw8eBaEi2Aw9s6s14oyM3B04xoAgIhw\n//Y2jm3sYPdkXSPaHqb9H7ZSUEeYCuLWG9jDQYmYHvtE/+1Q/DCW5Ffzmr1GhHM3s264Er/rp67W\nxhwT2TVvbY1y8b1WaJ5FkOu2Aw7tVmwNpvj3IbN8RWCF/jzzmbtg5s5eu6gsWJXFtbRwpWIgx5PQ\nUnG/qBCm3VxkK22pUbPDQ8bRjWvlzxk4tlkvrmq7Y5s7eLi/jbM37zSKcaUtndcUhHmIbpArrCEf\nvyhL6/pK7DycO+dd+tx6HVNQxya+XTpIh4eMI+vXcHRzpzRg1fY3FMnoEqYiKoFdku9MD8/yXeCG\n+RB7Wyewe7JwDUxdc5e5XYzdRpDUP6xEgAVXRCWwQJr2q26pmsuL1+z34cKFXRhyWFP2sLli/b76\nUE5KqIqpcsuYyQqC0IfoBrnU4ITPdEfVTfb5wzKtcSWmbKwzJSpWNWR0a95FeynfT9t+Kf9kaTHY\nQpnkB7mKQh1+Hgwu3RDFvsoPg6bsrDlENaYpeOxFcMbsr/t+lza4KownOheBTwnylZZqq187t6vD\nzKwK2S1gorfVpUXZ9D2kWodY8Et0LgIlgL66s5mv8+rgurM2C9jWzjlnZkiprq66F6Y8nrgKlkuy\nLoKpip4wyIn4tHe757FhXYsrEfBQpb4CYOZesaxjyfyn/n3m+vHKBWREbIV6IrNguTHLyAUufqxd\nrapGa4gIZ2/u98pusmHzsbr43nRhJeNZYd5SvkTXnI13rjqvKc4TJ9jpa8FGJbCALrLufky2EX1l\nqWSZY1dWn7j+AdelbTKAC69eBrj9x9s2UOX6CitxNYW1Cf026yK4dQ+6uuiL1b7nFlkt9bhoE49+\nUAphkLTA+vDBtrkeVBe0eO9e2JuOrUReLVH+YdvO2vL+XTBEXG1kiQVloSUCHt7exo2tE9l76F3y\nS53cRFOGmulRIhdevYSH+x/F3tYGzt26s6oEdvbmPm5srXdLaROCp6/ARhdFQNr/Xbtlbes2fVvV\nz9z+UNqu1IU8GD77o8qo/+qvQ97/qHYS8NbBRSfiqva3tkZ46+Di6u/h/kXQGuHcrYPSdda/87mj\nL0z09hGtZYJ65gRubK2D+VDEdeFEM8hlp/uNqzKA1LuhVqjLGqTluN7qZzY/c9MZ+/odu7Ja6/Zd\nZW2QFTp1TG9hLWdZdjfO5Nex7A+ZrD1CeEQusP3Qa0aNK1fn8kejt6q+sPVc+BTXJupiUc3RexWf\nqmdhTYU5yCZiKphE5YMFyiPvffxtrkd6Xftii/36oy6zrHb9GcRVVfXqclvq6cW2VGNfVIRVWAzJ\n+2D1sn1d/bBjSvDV77NZsDNx2m5dBrgpwdcV06fZuv4MCtL1ma9/Z1NnohFEXIV2ohPYKs0/K78J\nCg3HZuDTH15fVX0iALsnr+DTW/OE68RQ6lGfTiZUYvgehXCI3gdbN+jksxvXNpiyOjYXFmMRXjaD\nSwbhTf1dRwxuzBi+RyEMErBgAZslqQY/fP0Y6rqkdWI2JLzMFeq70Au8CILgnyQE1law2offtXy8\nYXG1GdObaWYMrSAI/olSYG1V6LNBBy4Jrd+uXHwyFXqL6wYBQ2Jx08kQ4eytg/b1BCtRCixQHybV\nd5R82LGrPlhlRac0AJINOk1XGQsAaI6whZ7MVe9gMEQ4d+vOMKHMU36FYUQ+yGX/MU4x0mvLsLpQ\nk5UVM2rGVyATP98WprJipxb2vlRnQAi/olZnocwL1qyK1MQw8hgo0SUaVPY1041ts2LqKmOVt4u7\n2LVuYPoU3MPDaevKDiXYpANTJIm6C2W+rdRRqJJ0NS3rvtAuar7QS5Coak/WtqibPe9yc7ZokupX\nvllbE5ENth6sehr2+Y2rOsRbGyKuFpKd0QAwu/7z+8F61TYgykR1fxvHNnfw4PZ2qfpV1a0x//l1\nwecDWol3+CJr/73N7jYaeG0IHaxdRwXhUycagTXjS+cI2G+i9pFmFGAmolWW1+7J8jmU98HOqnb5\nwpwuxgehi2ydv193BU3qqzVdAw73aV1+6052rLz2rQhumWiiCOpm9Qw+bEaNwjLjxpkT2Ns6kWd5\n8Sqf3ZbXXne+oTBlIRifboixmNdJT+Yw6yTodTS8QpRZoSO2r0QcEFX/1LFoDeduHow7ZqJEY8Fm\n1HXFqLauagiYhZe73oZh2ehVi3XKiKqQLVkzoqCxZi/0GrJXO7kQymGBIyzgLt16zizfs7cOspkZ\njM9uWGppEHJBFp9thWgsWKBsrVazqaiyTlBEfPOZsxnoBsyUqBkQQrRm6yqiFd/dtrGuyqzLxLnt\nntVntWhanzkLq6t7EHWxMpkPi8gDADe21lezNORz/ZT+9PUlKaFMZBYsoNt/avT2vJZnDzTPFCD0\nY66C23Vk7QijMU31dauz7lbbzKX/zbjaeor1y4XGL2j3/yhyS7c4YPsemQ9x48yJsUdOjqgsWKDo\nYhVPcS59Fu7sTf2YoyhMpQ2BiasihJRafdC1SNE26mH0tPS54c++flnYbev27tGpAa08btbmEqhF\nWbXCiggt2LKI1lWGSsGK1S30yY8dqLgqlLugyxTgviDj9aqCmyW6wleGmu7/VQWOqvcLNUcCWJav\nrFERzFFEn2jQeBzMl4Tgkqlz30MXV5O2hASb4LkQ5lKCwSo8ilbHHNLWsSihV6+Rv6e8Tbunrhai\nqVurCmW1irBaSTrRoA3zCZ6CFZsx3c0em7gCzREGdeej0nzHiR2tfJ9EwDlqH0LyXcyGLa9Lsbqr\nhQ2uABFXZ0Tng22HS/7L6CofzUiM4qowY2XNyIcu2/RFH1j65MkrWZxzgIH2qp2VuGpdXPXIAMEZ\nrQJLRI8T0R8R0TeJ6BtE9Gy+/J1E9GUi+jYRfYmI3qFt8zwRfYeIXiei9/k8gUp7oQYfWER2ADGK\nq0L5Zc2Qsi7bjBHaVYYhc2OG4ZxzjjGwyrw6d+tOZr32dQeMKXu4ULpYsA8A/FfM/B8B+I8BfJSI\nfgLAxwB8mZl/DMAf5u9BRE8C+BCAJwGcBvDbRDSppVxMz1KM7u6eekFEtoEp0l6nYEicLtEwa9ZM\n3wbqfbt9piN3QW0Uih7b2ldcc5eCZGx1p1X4mPl7zPy1/PUPAPwVgA0AHwTwUr7aSwB+IX/9NIBb\nzPyAmd8E8AaA9zput72tqJvpQH0qN0YTMVuvLhgqsiZqIIsZOLJ+LQ/8n05cgWoKLyGb1biSMNCF\nMaFbC6eXZUlE7wLwUwD+DMCjzHw3/+gugEfz1+sAbmub3UYmyBNRrxLBZnm14DsmNhXr1QV9XAZN\n91MmstdmdW2ag1zKjdG3McyMo5svFjNciJ+2M50FlojeDuBfAniOmb+vf8ZZrFdbCvYkNN30Ksog\nNnwUflGDQG0DQUtEuQy6CW3xxZnlCUPSobGXV39gCN3pFKZFRMeQievvMPPn8sV3iegxZv4eEZ0A\ncC9fvg/gcW3zzXxZCf7+F4o3x58Ave2J/q2vb7F1iW067VhwdV9XUziFOjJfbkMIGMqlCkN7gKv2\nDHk4lwvM5Mt4eU621+/dx7fu3R+8favAUha49ykArzHzb2kffR7ALwH4jfz/z2nLbxLRP0XmGngC\nwFcq+/3hDwxu9FCWdnPoiLAOxxa72vTArssunAPKq3aZ8eCFgHLNckvhGvif6y403v3Icbz7keOr\n95//5t/02r41k4uI/j6AVwH8JYrv/HlkovlZAP8+gDcB/CIz/3W+zSUA5wA8ROZSeNnYp9dMrlQy\nuHTGhJrFHN8aCmYGlnmPOS224hAzm0vHzPgyl9v2ldrvqi/OM7mY+X9Hva/2H9ZscxXArH2luid0\nvAyc/kPE1QmmFat8/VkXPDxhVdgyu5rWEdySVKqsQolr9tSOW1lsvrCuiLi6oS5BoKjpKgh2khRY\nII2nsv6g6IP4W91T50lL4T4T/JGswKbCUHEVYRWE+Umw2EsRmB/1qCcRzt28075eeRMRV0EIiGQt\n2Cw8Jb64V93nen6zeyk9EVc/zFmgJWTSG0T2Q7ICC8Qb98qVF/WIvzXji1/8Yus6RISnnnqq975D\nyV4qemTzVYerhqOpGRVEaG0kJ7BRuwV6IMLaTVR1mHm1zVCx7Ur99C1j9qd6ZNNPJaT/rrhmuVAl\nKR9s7OmwfZlr+uy5+eIXv9hbXE2U2L788svtKw/ZP9zWj6hWx+LKJIu+0EXUFjJYVKwTTJISWCDu\nC9015nWp1a9cCKuJT6F17VlQExwqyzUTtmm65npGWLlNcVaom4qkJj2MOZWvT8yrqoS1FFyLahO6\n62DqItl9MIXV96wddQJb155U6Zsqm5wFGyNDEwqWwJTiCmQWrW7NhiiuQFVQfQtcWz3SuPuO/hCB\nDYRAf8ezMrW4KpTbIPzwLDLeuffL+i72njrJCOxSRjOX4n+dS1xXMOPBnefmbUML9VMkFXPRjRHH\n4jfV/vgXX6ydZAQWiLOT0qeYy1KSCWYX15w1MB5+L3SRtfteCTax7S64QyJyZPbmKknEwcZuvXYd\n2BJxnQhmnP75nwdQiOzRxz4+c6PqUcJWJ55mocViILgqiPpvqe+tFtqMDiEQvcAuIfZVxHUGtNEt\nFyLrOvHARIVwdcmo0mcRs0XdDLnNlvA7HEL0AgvE6RroylLENRg061VnjMjq4uMzA6uculr+pJt1\na99nZxdWp7WWRfQ+2Jid620375LENVTrVWeIT1afeDC73s97z8Diyl85QaHvb6aLCyvm36FPok40\nGOMvmpsusa9LSihwJbAf/3i9hUlEePbZZ5t3wIzTH2iekPMQ1MuSJag5sbIJCIugfcKFfNqZqcSp\nOkeX3bodOsdY6gkHzufkCpWYfT6SWFDGhbg2CauCmfHxj3+8m9A20Pee49X/pPVYMrF1XbOga1uK\n16ZLoRDIgbPADWxZmkRpwcYsroDqTjX7tcQ90I0uwlpHndASM55qsGIZhCM9LVhVZjCzWNXr57F7\n6oXgutZtabHdtk/Tkl1MqmzKuiPi2o0x4gpkFu0nPlF90DMRXv7CF0btu7Q/ZFZquQg8r1wEodGe\nFlug6n+U6yJITKwiKoFVFzPmmNeuLEFcxzBWXBVNIvvFL3zBmdCqgSZTUFNwE9VX9ZKbOBqB1d0C\nsV425RoYOg234IdaNxmR1ZqlgRle+oj+hVcuJ5jnz5V3S7dio/DBxu5zBboPbIl7oB1X1qtO28CX\n6Zft64et3S/it2L1W9V2LuXEhrgfKMn6YFPQm85+rRRONjLaDA2XftnCb3klidRSPea2+fPl+WWD\nFdjyTZi+zxVYTqWsMfiwXhU2X6yOK79sMeiVbtxovftjWdZDkAJr87fG/uTr4m9jBo71mKpbcEsn\nd9nKL/uvRx6NNB9l7E6CMsXvt/zwiL0o0xCCSTQwv3w9Fk9lu8SOuuGaHhZz+MSF/rADSyy7519I\natCzOl7CNcuXwewCa6a7lvPzORdWSujCtPj6GDi6sYOH+9sTtUeYC2U8pPZIJeP19UjT2V0QRBSB\nOcpozrvOSOvidHH2r61R8iI7JIrApw8WAJ57rkP4FTOe+vmfx5FHx7UlhQgCE2Uw6UbS9VVKsD0l\nNyaiqkVQ9cnYa1qmJa7dcs8PD1ks2ZBxoIypiStgC0csah2UDaeUftX1zDrIVR5NzZakHs7Rp+q7\n+GMDpKZerFBg3rV9Um9TYzYL1lZJHVjGSGOXwS4hYDg1p9W0LOE3rpjNgiXtT18WezpsdzrM1JkP\neKXKU0891XubMWUG2+jkf0XWwzr62G95a0fqTF2icU5mE1iVl2+bdG0JdM1DPzxkHFm/VvpLRXRp\nQMrakG2cksfAugjTWhLm/Z6yG1Bn9jCtFOJbh9PNM2W6YpmXPQD27LPPtmZd9aWr9Qq4iYFdArbY\n9ut5LdxMXNP/Hmd1ESwx8FjRZzI5G8qyjd2aHeImmNuKff/7Ty/24dYV091H2nJlWCyh2tasUQRL\nFVfF2JFV5viFdqhYuvTF9rFeFWtrhLcOLo4S2vTKFWa0GU/6eeuugjq3YczMlmjwqQ89OvlxQ6OY\nE8nNDUWUCVZs1hUz4+WXXx603VhXQV9xPX36tNEG4Mj6tcHHTy3ZoGvP1Jw2yZyMsS7KaG6SLVeY\nIq6rKTH4OZ/zAAAgAElEQVTHGTs71IolIjz33HODrNnnnntukOVabQNGPdDiu1r19HH7VQe90vou\nFLMPci2ZsX5YHZXAcOHVOLtXTz311CArFiiEtotFO0ZUTetVP/7S6TumUheqlVqMrAjsTDCAj7x6\nGQ/2txv9p+03nD5TadMBGTe21vMtCGdvHfRus0+IaJTIqn24sEpt1InrGNRDMXaf49AB6zqLNaXH\nlbgIZkBZrocd679S4x+XRmuvn6xxNxDh7M19gBnMh7hx5gT2zqw7OR9XxGoJDnET1NVMjY0h4lrX\nc9ONiVQGvERgZ0DdSNdPXcFunSCie3eJat9oMGNva6P0PkShHRK25Zsu1mufh4MpSimIbF9M69Um\n1EUl2XgRF8FMFLVvjUw2Ah7c3saxjR23ccJEOHfzAAwuCyrzalkIbgMXrgKXdBHXLKW5XyRB9brG\nN8Tj0l/qcjwiJMSCnQF9+u5K5SEGjm86Ftd8x3tbGyBaswop82EwlqwS2bnp43ftGrxhE6WYxWWI\n37XuXNn4S2HAS+Jgh5L7NIeKUlsMrMp46d11JMK5W3dqDspYTVmrLNetDZy7WRbcva2NIKzZofGx\nLugqrsp67SOweoxn1+ncQ2OM77V16nqUZzkJiagKbkeN6dPsiS6gVpEl4OHti9jbGumb06IHil1n\nUQSEzG1Qmic8F171ek6xJSKcPn16UqEdEi3Qx0ZR6aFmDKhO2XILqyjKUPHTxdV2fsUyWtWITgGx\nYGfGzGjRIS3Nhwh4sP/RdovZYsFau/9EmdDe3C+LKDNunDmRuRJu7hevZ7ZopxDZoeLaN5PLZsXa\n1ik+D0dkh2ZZqfu8fjCrms0VImLBJoT57CNaW4lnZbCqYSd1lvbZm/tF+Fa+rrJ2mQ9XwjvGUneF\nT2vWR4xrHTa/oi4s5rxVu6euzlo7Vbcsh4q8Gsytm/yw7X3MiAU7lhG+2L4+qUooUFWBq11+1Fuw\nNl/tal19X7rvNhBcCK0LYR1mwdb71k1rdc5aBbq1qe7VoYOvIVulfZBaBD7JRUnvLjMzjm28OHiX\nncU1O1j5L29TqftuEcI+XXxlJZeE2iLsc0ccKIv29OnTvSIO1DZTWq06Knur/nNu9M9ORTVWd1x5\n0RTEdQjiIugLEcBUGqg45OrAxejDoEijvNC0b2Z8emvTHjmQd/mVr1U/BzNyoHzwhp+DZdBsbpTY\nhk5XkZp7zra6di5VJMcgFmwf8gGgvTMn8q7c5ZXPzOVARDWNsvnWVv5SqzCqjC1TFLVwrU7ZXLrV\nanMr2azafN+l/WvLKp/V7Sdg+qTJ1hU4saWFzpndtaQ5s3wjAtsXzsKgr5+6WgqIvn7qCq6feqFz\nV7wtuFxPo2wNtibDQtWXa9EB5nncOHNiJbzW7Y19/fLN2/XWqzlYpvaduzMYXFi/hqtj9Vm+TYgp\nvE2srXWvwasezKafddfyIJ1LZLP2vBBM5ELMyCDXAMzBqVLoihoU6rCP869cziyg29u4sXVitbcs\nLrDjj6tmYKt0LHOQK9/mxpkT9dvbYmBzEVztw3RL6O4D/TvQXBKr7S3nYf3e8nY5qQBmuDdcVxU7\nPOReM0uYV6wYCDL9sNOFaZkDW+IWKCNhWhPQGLrS84HFeUWt3ZP5+0pUYPsO2pIB1MDVKrTL9KPa\nxLnOKu7Qnjoafbd12+XLu9/S9n3YhN91DQZlyXYVWdsZFxl+86XO6nUyYilCYwtxC6Ht4iIYiCoV\n+M9OXcXDg4u9f6RqahcC8MmT5Zk3++6TrT/VygFBumgroWnyeWp+WtMCroSmdRn8GtNb6uortmxT\ncktUVjm0+4IH0sdd0ISZlz81bdEOIaH3KIvvKwzbWwR2LIR8Lqye2U7MeOnDGyuRte7T9U1ihnSp\nZZrv1Coy+jpUZH+VhL3J4nWVcqun8bas1yas5vq1+x4w6BZrXVtAua6yXlkMZxF6LQcR2LHoI+U9\nM54YDFqziN6QkfSO26ziXNWflkxwY2u9XsCypwjO3cxqGNgGzT69tWkX0Vx8nTwwupynVly8M7bQ\nNfWdRDboNoaYLFdFnaslhEE68cG6YORAoU14Ollq2vH3tja6+0xzQbVZlarr32RtVlwBWnSA9SGj\nxeOWfMFD6GIp5xZ21/3pwlo5d+3adv1+VehWrFOpA+7jun2RDTBf7TZt0gyIwE6JXvRalQms8132\nDOhncKNw1m6jH0sJpCr4Yhtlr0tS0Hy69oNlXfBRQmuk71YiAm7ul86py/7O3twvf8+5hVxZnq/f\nvakxdLCrqMGtmKxYfdqk0BCBnRBmxpGNnWyA63YetqT8hPYNuu1YCYUKgdLCqWpDkUyhNAWS2T5y\nr2eHDfGpakLbGaOt+nmudruqmFI8YDolT+QPlNLiPiKdJM69/96xhbyFEAUhPtgZuP/dj+LoZl6/\nYOyPWBt0Kg3oNA3caNuu9tEklrrfkwxLcShdBp90ch9w7baWSmCd22FZVvcAYEbUXf8uhOK/HEpI\nVrgI7AwQUe+g9JodlcV1KDbfqWE16skKRGtl61cX6LrXQ8iTGVbHanOb6N/B2AdXw/ZtyTmZCPer\nsBUe7m3YaYU7DCtcXAQTQgB2T17F0VzLRheJsXVx2wq5WPZRcVPU+DrNsCydUjKDEuu6Qa82TJeA\n4SOdk2wA62Lreov1LtSgh1NlSRTuEwHa0s/nQCzYyWHjxzfyOVsnjPpfE0aMa6VUoXkMbVkldEkl\nM4yxJPW02q31au2CAGgpNha19arSvncdVvMyY1V9JgLoyQYhuDlEYCdGhZUonN0IHWoS1G9q1IDV\n91HXzc+Pp+JNe8WJtuyz0dcaMEpcI2pyBeW/pHw2BZf7Nd/7F8D5nQQisDNQKeTiwldpsTrrygTW\nZWsBsK6z8mZZUmZ1Mew88KVcB8piNpMe8n3PPQ9YX4iAB7fHp8nOjSou41+e3B0hRPcAID7Y2SAC\n3jq4WHo/OjbUxAy1UoVTmoqctBRXKZUdtGzbJ/yqzeIOY5iiO6pwTxr4/+6VFevKF1v3eC/Ed/oC\nMIu0YJWjvdJFmdhqyuoNMPa2ilTbXqLSUGrQHBgyha+PtVknplbfaN/Qq0RIwT3gC9MtpuOrFGNR\n8a54PUfY1iIFFlCO8KL4ceF7muFHr2JW8/hN1WVuFPsmy9Usft2UzNC0b3Nb7ZhOXBpN9LSGvdOh\nBkJK4lprhAykub6xn98cGa/neJwv0kWgnqhKWAFa+W8m/VErEdHK/62my+6ST19nAZrhUaZANuxP\nL6JtFWZVADuvHtbJnaHF666Ob0l3tZ1HMAytjxsxDNdummmup269zs1iLVjzYq/eTfWj1q1DbdDn\n7K2DUamaq0yjum318C3bwJca7Opg9Xb68akKXCo5QT+2ZfqY0EKyStQ80GIPzbLhI1yrmFrpip9I\nGuNYIbBYgTUv9ixPPdN3yfkMsR0SBZjLmWDMwJH1a938gEpY9Tmz9Gll6qw1o5vMIFx49XK9O8N0\nY+T/tz4EIiShU1mhisq7nJ+r6KqbX1gokuiWxQosUFzskKrxHPIhjmzs4MjGTutNbabbch5h3fig\nMC1H/SPtpq+zTkvrEHD/9jaObr5YX3DcsPqUtaey2KKhxgebovUK6D7Yy/kMs25/Hd7iwQNj0QKr\nU7p9AojBzOIpi1bVDTocHvLKctUprZf7Vlcl+OrMLaOoi9UyNdwKqq7CkfVrOLbxYu352C3sEB5p\n3bEV3k4l9tWGnhXl7wj6u+kmd5wKEVgbQ3Poxx4WeUphHk9pPtXVgJwZVF0qoqX2oXXfsxq0lAlg\nW6ESo96A1ZItRT2cyKYwZ1Vb4XLpuJbN4oY562WsX8sfGDvR1n5tQ7cyfflJzciCIrLHXQTDnMi0\n3YGhP8WVdXQsdwOoK0XoZlWoZAZm4NjGNXzyZDEdc6lWphldoAkGM3B0/RONAdpFAY/ycUv70CzX\nwt9NuPDK87j+sy9UJ1EMHPWwO//KJTVmlyTmvdZ5OvmOmPeO7fhqTdtxdWOjrp1lg2Rc+/tO2y0C\nGxg2sepdLwVakLU2sHTBsIbJ2EjNdLtaRxPGth+W2b2z5D6U2nc9b9+qHZGqVIrd2ib0a+eCNoE1\nj20KZNP2pSRG6/L+YttXYMVFEDiZuNVnwpgocV2N1q765VzxqenvmQE+LPyrma/xhOZ6oJZum+FP\na4i40rubpG8QEQxkrpA+pSGFURTJQZdW7231B/Tfy+6pq9ZCM7xa0y8isBFQZJkVN04RZlYIry6u\nfVlZvSVVrFoLNout7kZvP2LcHDJw7PEXS1Z/6rj2xfYxHoo2UB7dYJ+uW/1esv1PX39AZ5GZXCGj\nZ5npqK50dsMUOdasPYvrxLVJAM2cbfvW5hJaFU3ePXU16HnpfaGu04VXLyU7yFWP6/MdkFDT4fNy\nDE7G1PHuYsEGSNtT18yr7mK52sK6s0ydK5X9FWFdd2qFWVmzSxRXBYGxe/Iqjm3uiBU7gr5WbNf1\nlX9W3cOm+2yKx6IIbKDoItt0Q6mBqbqbpc16rWynRRTsfXgdu6tstyu1o7hDqIvrjQ0C45M/cyWh\nMoXdcD2416crr7sAuqxrH+CahvSjCHLB6F1nNQCUOKqnbnXkv0hFLaW6GvswR1lrLV5juhbbwJOP\nUXMj12v2qZaHoIdtLQn3YVv+vkMXERASpmU/IAjFdNMxCW1JVLWqW7bpVW6cOVHZ1taFt95oHcS1\n2K9fIbGF44RMXSzmUohLZCUO1uNx17J0UUOIokGVE9TKDwKAreZrnbiudmXeaLnANomrnqDg22Jz\nHW/pm6XFw5r4SEDQ7y/74NSw73xMWyUOtoHOVfxDxyw32FFcy6Fd1UyANss1u8kLv7BPe81HaqZP\n5g4Hmh/3xWBUxIw5GEulz/t/51M+DBclsKrGwNyFXHqjT6dtLFezupYWA7XTY+jVw6z7amjD2VsH\nlegFv6ISV/jTkkXW1wOxbcR/nMj6r3nQGAdLRD8E4BUAb8vX/RfM/GtE9E4AnwHwHwB4E8AvMvNf\n59s8D+AcgLcAPMvMX/LW+gGoGQNipN9sC83WZeWGbSpwo/l7qxMzJtIrcECf0e00meuBOOwe1KNj\nzzuqVWDSaMEy8/8L4D9j5vcAeA+A00T00wA+BuDLzPxjAP4wfw8iehLAhwA8CeA0gN8movCs5MjS\nMlur/3fdDVqsDNu+csv2xtb6SlRjfUD5pM3nLUxL395VkTbu1n3QKn7M/O/yl8cBHMvb8UEAL+XL\nXwLwC/nrpwHcYuYHzPwmgDcAvNdZa5dMRyE1g6urn/e8gZRlm5cnVELfdhwXxOSHbXLLCPNQ+HD7\nXhd3lnirwBLRGhF9DcBdAF9i5q8AeJSZ7+ar3AWgQgLWAdzWNr8NQMydqSDCM7fu4OH+dp4gcFX/\nSMs26nkDmVPboBAU/xZbPH7YJftggbAeiKotfQdjXafSdrFgD3MXwSaAnyainzQ+N5MlKrsY10Sh\nL3sfXi/dWKu6sirbiOBooG/hl9aY+UJ8sMAcD8Q6d4BeFKbv/lzRudgLM/9bIvojAE8BuEtEjzHz\n94joBIB7+Wr7AB7XNtvMl1X4vW/8YPX6xx85jnc/crxv2wUTS6hVZrleLBXvvn7yKgjLyZ33SRZb\nfWc1N9rCHzkry3FqS171HmxFkvpQZD5m+3v93n1869794e1qSjQgoh8B8JCZ/5qI/haAlwH8OoCf\nBfBvmPk3iOhjAN7BzB/LB7luIvO7bgD4AwB/l42DSMHtCVilCFOp+DUfMva2ToBoDb988/borLYp\nsrqCTjjQCpqH0j2ei/bKbP5xlQBTd9+5TjQ4AeB/I6K/APAVZD7YLyAT2Z8jom8D+Af5ezDzawA+\nC+A1AL8P4FdMcRUmhMriCmbsfXg9f3mIT29tjnIVSJcYq0iOpftfgeJ+mNNrPiYBwQeLSpVdBPr8\nWpZ5W8w0YZU+PLRGw+ItWANJmQ3nepXjXNuviRq4vfDKpdryn5Iqu2SMAjB7Z9aLdFqVUmtwyIc4\nsrGDoxsvlrJauo4Ih2QtzEZpsGvZHTZvkQR108g3bYLuFm1RYY4HzwpiQwQ2NfIaBXtbG/jlm7er\nCQoW9I/KsZzdbrPFiywz9s5seI8LjgFv8cD5PW2dRr5tU3RzZZHxvwtkyphEWaUE98z0Gi4Qy7bc\nsKrvsGwXAeBpgKutVkbTpshEf47rIgKbEqaPtU1cifDMzQOcZVTCjKaeu6iOucJ+hlDUsRWhdY4K\nQRyyKeYbjBUXQUwYge1WetQnYGYc2dixTncSSoRAjCmovks5LpaetTdKmzZ85tOYEIGNjCE+KJ1i\nLqx82uP8nrVnw3Q/ls/ShfEkyxb4L+UYJiGly3bFpzEhApsCXSxbDdb+9GVjb7LFDnbVjHAv9vtw\n9EjUjQH/ou3nMS4CmwhEa6PrC7jp1LrrGofiB25FjXBbr4G4CsZQLiM4TGz13sTUPQsZ5IqJFkf/\nWPdBsZ9hwjYuCqG+LcW+wx3sUtMRmbUJhGHY7iV7gezqlvZ7hLX9lgXa54NcBDY2mEvZWozMemIG\njm3sYPdkw6boJoDl9bp3c4sb1c0IelF4IzKIoGYDnis8aE5cRX4MK9FHrce2bevLzy8CGwjqyZql\n6hXVfGrREgqyQPcTYL6k7auu4HbDLlEIZNGW/t0xl5NL6kJ14dVLOHvrILxp19WsD1qI3PVTWctD\niMSYA+4gdE2orvywh1MRx1H+TU2PCOwE1F/osogySJt2xLhB9RoDQKk0od7tGTp1SZEqmO1l99RV\nfOTVS3i4v429rW7dJx8ugmK/lwBy5wZxivagU6gHzZI9sGNF1mYO9OklZeu+MOtUPlLsxRMVsSHg\n4e1tHNnYWdVmPbqxU5KLtq4LrVE2K4FR91V1pdXN17dbXRbXoi3nX7mclTXg7m4C18VOyiXwCOd+\n906WpRaaFWsgc3QVDJ/5teoi0ovJNN1r6ph1+zB9rl0f21LsZUZW81S9ehnnbt4phUMxA0fzgH5m\nrIL7bSFT5j6LfXBFXM11h1qQprhmA0pXsHvyKrrcfqtzd9wlLpfAyzLVYpl0UcQ1o4gAcOWXz/aj\nRNSMDNCTU+qiBsj484UIrAMYKInqoSampfXY/rqJ0lxaDKu46jfRkB+1ftMWlpfyH/arZeBSVJS1\nsnvqheLHOSKbZyp8uUpiphxqNVZoSXulakBUfxPquPpDf+rQPxHYkagfkymqzjSAgU9/eD0XmisN\nO2ZjYsO+hynmL2Jt2dxCUZSbe8HRPGLTEPYjYD4Kob00KOurrntvGhl6tiJQtXCnQgR2DEQrqxXw\nY1gxgPN/fKlTpXjd9TD0WOYpdD0l/wHcRUnAkMOexHrtRvHwtt/VTd9j/e+AtV7PlZWbwGbhToVE\nEQxAv/jnRwhan+O1jZ6PC2tp3m/3LpXbm9iMp1Sj8kFGEmiI9doN9T3VDcravscuMbakrT1niBYg\nAjsYrrzwR3eRc9sYW3TB9OgxjWFbhjG0MUT637X2O1Ldr+Z1ULHU1wfEdY9FBDYSsrhK8wfsp5jI\n0NlBfVjRephNqJahfl1ct3Gp2WBNFPdE9f5X96seK6sbClPXtpA42IFMPbmdPipaXl7E+3WJu6z7\nwZo33xir1ed3E4ZVXeA73lX5FPUUBhHbgrYYW16t5waJg52MaR9M9SFQtPq3qTC17vyvuyldxQX6\nLtO3JHHVexJFRMWVYIveTE1blS3fca5tiIsgKeol2LT6lAias8i6++G6kZ3QShba/HuuWdVeqIkc\nKaamEdcBAMO671DHY0JEYAcwx2BGWWiotzDaf6hq/qiwWYrFCpQfhu0Plvm90qH5iIs6HmEIrbgI\nBjLlra3/6IpuYrlIzNApo32lDKY6ZYpvX6ve0wjlwdJE9tB5PrhrrSc0zEkSAltkbvj/MuewXvV0\nv7qc/7Z6Bs6/m5ppUkqrjLQgQnIPuLzuSkxNX+rc/sKhFGnVITLvN5qEwBIRHuxvo+/cVH2Zs0KS\nmcratw3lwQBHDyOiDgkQ40U2FFzXWbD1RubARW9jfmeFHS/GRQ+SEFhmxrHNnbyq/4t+j+V17+3H\n7nL8wkIq/2hKlbnGdJ/M2rSthPrzmx67Vd7/+3FpUZsFUVKhcLvM9wCLW2Dzbuq5m3dWU09/8qSf\nLmVMWTqqrWoUut46GWYfMjOObOzgyPo1HNnYaRTqmL63qdBjmvVA+K746Eml9gg0xy3mIvooAgbh\n6Oa11XtfX2afSupT0dQm/QdTZ50MDc0iALsni5HjOjfBnC4Vl7i2Fm25933v29i/0ykIwb0UtwWb\no6xX3wMiofjMVFv0KV6aaP60/23YpUuZirgqXFuLQ/zoIZNq1MhYohVYVcbv6Ma1ydMnQxJZ3bqq\n8736oDli4fJocQ0lgsCni2PIOYbqcknVjzuWaAVWBTirUgrTdwfmtT9MISv7nPrModU/xM20VsyR\nWlfW2dxdPN9WuNp/1+8/9F5BKO0K5eEMRCywAEAEvLW/PcuXGUKXyBQyWyps1/30lTOXyQ51bZrT\n1z2VmHWN6ghdXKfsPbW1I6RiQPEKbD4X/d6H12f8MkO93TO6iuzwJz5X3rmvSDsdStTnELNVjwzV\nHkUM4lo3HjCHIRKKuAIxRxFY5qKfGl+zCLjELOpSv55tYf4Qq5kiO4bz74ouYnOIme5bLefTQ1sW\nHnUWY0jd9DmJ14IFZhXXgnnboOew16/D1tRMRVN3vClTy/fAxhRuAtug3FxX1OwBxBJtUCeueo2L\nqazYuV1LJnEL7MLpE0zdlJpZW0s27yU0pR/7//H77/DFIGKxUb1q033DIRR5UcQlsB0KjCyJ4RZk\n1XdaN5rNYBCtWb/zEAb6xhBqyFNMmBZjOK6BMDyxcQksM25srWe+V6He8hxAYcUZNyYz9s6sW10F\nMcc+hj5wFDLVBytVPp+bUFwFcQkskPldg/C9hoHrm9l2YzI4H/C6Uyr0HbNAuXw4LZG+yTbz9Hbm\nl/p4owiEwbSP/pdvTGbGkfVrq89WXWsCHt7extGNHY+t9UPM1vfc6PUU+rkEuGZ9Pz7TobU2XBKf\nBRsYc8X5tUUONNEmLmYspiruUuo4EPDg9jaObfoVV19dPbFgx6EPDHa1E23RBf7re8xrxYoF64Rp\nOsr603/MbdNlHqUiFlNZALxafuGVS9g9dRXHNqfy1vj6kcTq4AiHcnWwdqGsu5L188OFExEwBBHY\nCHAlrOV9dktAUCEvu0Za7IU840kQCkuWVg/gIRZpXdHL60ZJzqqbIVwRFheBA6ZwE/goHNy1a5YJ\n6mVj2XT4cBNIiJZ7iu/U7Z1auBGKRBndxRByT0QE1hnhXuQm+vi/5o15de8miPOKhY3yzaqaDi7n\ngFNCa4491Bk4IcTkisBGgK+BniJNtGupvPCmZx6CWK9+yYTthdV3nLmXigfkkBKZ5v6r2B+Xcwdq\nicA6wr+bwM+t0idNdM7pmV0/ZMR69Yfyw7L22raO8u+PvbYhPzBFYB0SykwHQ+j6gJi/tM1w9MIu\ngl9sRWuUu+BCqUoYjfbbqqgY/R4OwT0AiMA6JxaRNZ/6Swi812uuCv4xi3DrvSVzRoyxVmzhm2Vt\n2fyIwHogJpFteh8abtwEoZ9lGlSnMCoLrd1tMC7cShdUqUWQOH1mE5hcjInwzK07eLi/Pe1xnRCC\nXSI0YSvCbYZT2R5zrrPrQihbKALrlXZrqUvX3Lk/iRk3zqx3TnMNqSxhKJaJUEVP4e7zGBy6nYl9\nsEtSZRdPl06r+6pZ3JrmGspAQRWxYkNlyJUZk/1V3kd41d3EgvWAy9FqXyLX1Sr1kUE2lqFWbMjh\nPCkwpnehZiUec11DE1dABNYbLqch8Sdw9S0M13pVDPtWQvwRpkST37Ntam9rwffW49WLawj3sAis\nY1xaSXPdIKHNLW9DfLEhU3/ntJUnHHJdmx6ac9/DxDPMDkBE/KkPPTr5caegXL5tGD6qZ9WhWxxt\nxw1hVNbEKA3e6PZwcW2EZvo8nG33U7/t210DrkMmn/nMXTBz55+lDHI5pIv12tUqnerJa9bhbD5u\neB3scov0+rXV9cT/6o9hRkFz/YDyNeNKckK7uM5fVF0E1hGdLng+xcreVli+zbm7US5RgyU2aza8\nx0NauLqPdGFllIXS/KypLSG4ucQH64Ciyn+9A58mmmJFGDZYIoxjiO/UFsmiGyqFuJaFsm4AuRhE\nGx9T6woRWEeoSlN1Ph9m4NjmjkyIOxF6STxxD0zFkAgAe6WtoXstF+Ken3BdBEQ4d/Ng9ZYZ2uyl\nYeX6qydxcbMYnkFlvW7szB42siSUJRtqjGRqVAcRm6MFXF+X4vjh6EO4FiwzbmytY29rAwzCscdf\nxMODi3iwv52J7607OHvzzqpbMraI76imoqh/WXf8NQrDJ7RERFyng6F34d0NmXZN3tHrzIZAuBYs\nADCDwdjbOoHdnwGItgFQZs1uvoiH+9t4ePAsCBcBAOfyC0q4mFu81yZ7kvHq//L8Vcp6Pbqxg92T\nkzTFG8pSD+XmFeJkqMum34MyDFMmXAvW4Nytg5VoAoCK3yXK/mEQjm5cy/42X1wtm8ui1Vkj4JMn\nU3ENiD0ojKev9dpHkEPI4FJEIbDMjCMbOzi6cW01SMQMHFm/tvpTnzEDh4eMoxvXsrCog2dx9tZB\n8wFG0pTXr3zHYTxPxxNSZS0hPnxZr+q+DCU8SxGFwBKA3ZNXsHuyHAalBFX96ahlWVjU2gSWbPUW\nIAIe7m8nFTmwhJkPBL+4tl4LUeWgxBWIRGABPfyCV7FuTYUjFIeHjGObO3i4f3GSduoUAp8OIWTH\nCHHi0/dKxv+hkEQtAn0WnsZRfE9+0C4VfUK78EMJuTScEB7mfe/6vpn699W3FkE0Fmwd6gef/VXF\ndYpuum7V6b6gkDJKXFBkrIkFK3SDjT+XdBHXuauuJWLBNkcLqFF8/5WpMogIZ//5Afa2Tng+4vSo\njBtZtwkAAA5LSURBVDWxYIW5UUaMiemKUCm3LhIQFlpNq6FwNAH3b2/j2Ib/7I7Vt86MvQ+vez3W\nfLCIqxAEZuZWUyEYlfhgL1fpTxuidxG0wQwcf/xFPNi/OG1XIaXQAQ0J0xJCosjculwqEtO8vvnn\nL14+eYEFikiC8MYYYyXNh4cQL2N9vL5ENgmBbbOqVB3WULI7YkesWCFFfIhsIj5YszJ/+YtS8aix\n1wIIBUk2ENLFbS83CQtWUU5GyMOmEsymmhtJNhBSxXVYV5QCywAuvHo5K1lYU2dAiSwzcHxT6rC6\nxvVkcoIQCi5dBXEILFFe+zWr90pEeLB/ESAC0VrpM2MzvHVwEQ9ub8vwlgdEZIVUcSWy8SQaEAEM\nnL11BwRelSfM94cHt7fBzDi+uYMH+x/F3pn1Yjsg2bCpEAipwLEguMKWyNA30SAKgWUAH3n1cuZL\n1YRVbzpRIbQExo0z6WVRhUwRyC1CK6SD2UtLrxYBEZ65dQcP9i+uxFUvT6gPYhXxruXtfdeDFaqD\ni4KQBuOci2FbsPnEh7rVClQnECQi6OdBuVtgNWkiEZgPC7eB4AWptCWkhukmSMKCVVECZ//5QWkm\nA/WnogJW0/Myl6fr1c3bXGyJ1sSSdUlNz0DEVRAKwks0IMIzNw9wloGjmztlPyuKuXaGGO5KZMWS\nHYExnXrpI8ikiEIauJrXKygLtpj+mnAsF1flY9VrP3YSV+aykObvxZIdST6duvqzrDB5kwTBJb21\npmlfwfhgNX8rkE1oqBhaz5UBfORP/gke7m9n7/kQe1sb+YciBL6QiAIhVtqKeEfpg13ltuf+0iK+\ntYgQGOYSUPGx+cyutIazN/dFXD0jEQVCbPiagSQIgVWZWZkQlqfmHprmqhckObpxLZ/Kewekvj4J\n3/KOekAKQsi4dAmYdBJYIjpCRF8lon+Vv38nEX2ZiL5NRF8iondo6z5PRN8hoteJ6H1t+1ZCWJdw\nNcZ6vX972xDszKl77tYdnLt5IP5YzzBb4pIFISB8T5rY1YJ9DsBrKEYwPgbgy8z8YwD+MH8PInoS\nwIcAPAngNIDfJqLWYxxy5nM9unGtZPUMHcljAOf/+NJqoGy1PHcVSPjWdIg3RgiRqSYlbQ3TIqJN\nAB8AcAXAP84XfxDAqfz1SwD+GJnIPg3gFjM/APAmEb0B4L0A/rTtOPoP8fjmziq4t9PJqwGyla8V\neLC/nYlpeTU8+G5Wp0AXVfL6FQuCEApjQz370iUO9n8E8N8A+DvaskeZ+W7++i4AFRKwjrKY3gaw\n0bTzupkhe588EQiM3VNXceHVS6tsLp1V4e1TYlYJwtLw7Q6w0SiwRPSfA7jHzF8lop+1rcPMTERt\n84xV+L1v/GD1+scfOY53P3K8vbW1RyiKuxAYuyev4tgm4WFuxRJl6bVHN3bybLBswEsGYPxhTp0s\nCHMyVFxfv3cf37p3f/Bx2yzY/wTAB4noAwB+CMDfIaLfAXCXiB5j5u8R0QkA9/L19wE8rm2/mS+r\n8PRPvn1wo62UahFkgqp/VHI75OFbqVEWtfnDpKSfIITAGMv13Ybx9/lv/k2v7RsFlpkvAbgEAER0\nCsB/zcz/iIh+E8AvAfiN/P/PqeMDuElE/xSZa+AJAF/p1aKRMICPvHIJZ5G5A97a38aNLWXdFqw8\nCHlmUszY54PX54CfTmzFchVCwBwgn2uUpW8crPr9/jqAnyOibwP4B/l7MPNrAD6LLOLg9wH8Cs+Q\nKsbIBPT+dz+KY5s71fg2lUarxDWyoe5s3qDyDA62aYvL875fts764LIdejWtuL5RISXMuNY5h7DD\nSZV1s+NVui0RAGZ86syJ+i+YKBpxNS1D1Woy3neh+D76W7bN7WBtMkRJlRXmoW7g3AV9U2XDq6Y1\nBkt3v/GbiEBc7d3/8uddKLpMuvApN0I3oW2r97p76ipotZ6IqyDEb8ES4ezN/eRKEJqWYv+Ei6oF\nuZpp1/pZ+ei64LaJvDmIoFwHgjA1vkOxkpyTy7KDQlRzt0BKc3DpgjYmKNoupEpks3d13XjzeHV3\niU1cZVYDYS58ugeAJbgIVMHnvFjL3pm6uqTxUQgrO8k2UWJaFlFa/csNMthFIEVcBaGZKC3YUv2A\nmvYz2O42sFTkr113QpTfUg0SuezimJasXkpwzJTburUg4iqEgs9SmYtxEZy7dad1NetEh3XbMs8u\ntAz3ISV1wqdG/et8sm3o1quIqxAa+v3tkmUILHpUwcrPrySeRCDU1IMNQGj7orsWzBuqbcBpaJiX\nzTVQjVIQhHlxbc0uRmDzHdULpQ1TPFXJQss+1PQysUQoKCvUduVdXuG6UVplfUsEgRASrge9liWw\nxQ7Lb1tE1+Y6sFrEeYFuq6shMKYQti4hMCKwQmi4tGKXKbD2gzR/bjnvOpGNJZ22XIfAfTe9qzVg\n+WbFbSDMiiuRFYEdi+l20EohxoCPAScXRYrHRCsIggtcDHxFOatsUDCD+RA3zpzI/iKKsfU1mr+a\nN23EPmSWWWFuVHLNlIjA1pFV5o7CNaDILE33IuZiv1KfQAiBbIxguvswvkwuoRbXIuZ6/qJ4HlVC\n2kxnxYoFmxgqXMpVd9xVPU01QCZuAmFJiMAmS5j2orgJhLmZ0k0gLoKEUJarGpCqFnrphz2mdVyR\n7jBlX1ge07gJJEwrQRiEC68836ksYfZpNojVVYwtE6I3im6RSjtO8AXBJUMiWyRMS0BWeSGbrSCz\nZjn3f1Z9oMVAVvcHLVf+MgG3zdFVFlSxX4WQ8G/FigUbAeU6sZkLQFmnddN064H9tNqu3E3XowRc\nxM+q4yjU8a6fuqLVKpCEAyEMhtQpSL/g9sIoJw/Q6vV5TSxtcarKalXuAV1AzfCrsbPAmlMkXzDa\nBui2glixwnIQC3ZCzHm2+mxXh2mF2ksW0qqY9/nVANgVY51sb2OKb5vtlRkPhNAho1fYhtQiCBCX\no+imtaiW6fuuK1lIHdZxXQmrXNhbKm0J4aF+C7unrrQKrbgIAsOl1dZWLrDpqutd9ULo2OLLdY+P\n9F1BcIX6bfr4DYjATsAU4tqXwiXgNz7Vt3gLgivUb+B8R3dBF0RgPeFDWOrE1ZwyRtUk6GKdTuEg\nEp+rEBMq0sU2/VJfRGA9MpV46VEGukXq2zoVhHQpx70URsqz/fYig1zdaIo5rd/GbcxnebptST8V\nBJ/o8ePqN8bfe1aiCFyhd7XVexXKZKuOblvfz6h8cTxBEPyhR94A/QV2US4CUwC7bUOaSGaCWsSQ\nUkVA9RF6H6PnIqqCMB1jf2+LEtiMbg8fexfcLqi2bVXWlSAIy2UxAqusVzNIv35de+xqnyeaWJuC\nsGwWUU2rENerFR+mOQgl6ZyCILhiMRYsQDXpod26/YIgCH1ZkMDWI4IqCIIPkncR1GUySQqnIAi+\nSV5ggXoLVSxXQRB8krSLwGaliuUqCMJUJCuwejRAMamfpJYKgjAds7kIXr933/sxdHG9kE/K51Jc\n+f/7jsO9zYecR1jIeYTFmPOYTWC/NYHAAkUMrBer9X4aN5CcR2DIeYTFiPNIdpArm3dKL9IiCIIw\nLckKrLJcBUEQ5mK2coWTH1QQBMEBwdeDFQRBWALJuggEQRDmRgRWEATBE5MLLBGdJqLXieg7RPSr\nUx9/DET0JhH9JRF9lYi+ki97JxF9mYi+TURfIqJ3zN1OEyK6QUR3iejr2rLadhPR8/n1eZ2I3jdP\nq8vUnMOvEdHt/Hp8lYjer30W3DkAABE9TkR/RETfJKJvENGz+fLYrkfdeUR1TYjoh4joz4joa/l5\n/Fq+3M31YObJ/gAcAfAGgHcBOAbgawB+Yso2jGz//wXgncay3wTw3+avfxXAr8/dTku7fwbATwH4\nelu7ATyZX5dj+XV6A8BaoOfw3wP4x5Z1gzyHvG2PAXhP/vrtAL4F4CcivB515xHjNfnb+f9HAfwp\ngJ92dT2mtmDfC+ANZn6TmR8A+F0AT0/chrGYI4gfBPBS/volAL8wbXPaYeY/AfD/GIvr2v00gFvM\n/ICZ30R2A713inY2UXMOgH0OoCDPAQCY+XvM/LX89Q8A/BWADcR3PerOA4jvmvy7/OVxZMLJcHQ9\nphbYDQDf1d7fRnFRYoAB/AER/TkRnc+XPcrMd/PXdwGEP11uRl2715FdF0Xo1+giEf0FEX1K68ZF\ncQ5E9C5kVvmfIeLroZ3Hn+aLoromRLRGRF9D9r1/iZm/AkfXY2qBjT0m7D9l5p8C8H4AHyWin9E/\n5KwPEd05dmh3qOf0PwH4UQDvAXAHwP/QsG5Q50BEbwfwLwE8x8zf1z+L6Xrk5/EvkJ3HDxDhNWHm\nQ2Z+D4BNAD9NRD9pfD74ekwtsPsAHtfeP47y0yBomPlO/v//DeB/RdY1uEtEjwEAEZ0AcG++Fvai\nrt3mNdrMlwUHM9/jHAD/DEVXLehzIKJjyMT1d5j5c/ni6K6Hdh7/izqPWK8JADDzvwXwRwCegqPr\nMbXA/jmAJ4joXUR0HMCHAHx+4jYMgoj+NhH9cP763wPwPgBfR9b+X8pX+yUAn7PvITjq2v15AP8l\nER0noh8F8ASAr8zQvlbyG1/xXyC7HkDA50BEBOBTAF5j5t/SPorqetSdR2zXhIh+RLkxiOhvAfg5\nZP5kN9djhhG79yMbcXwDwPNzjh72bPePIhs9/BqAb6i2A3gngD8A8G0AXwLwjrnbamn7LQAHAO4j\n84GfbWo3gEv59XkdwFNzt7/mHM4B+J8B/CWAv8h/AI+GfA55u/4+gMP8Pvpq/nc6wuthO4/3x3ZN\nAPw9AP9n3t6vA/gn+XIn10NSZQVBEDwhmVyCIAieEIEVBEHwhAisIAiCJ0RgBUEQPCECKwiC4AkR\nWEEQBE+IwAqCIHhCBFYQBMET/z8IgLgiBp8a+gAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x110615d50>"
]
}
],
"prompt_number": 46
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"base = '/Users/savoie/projects/sii_migration/pm488'\n",
"smmr_mask = 'baseline'\n",
"data_only = 'builtinmask'"
],
"language": "python",
"outputs": [],
"prompt_number": 64
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"current_file = os.path.join(base, smmr_mask, 'NH_seaice_extent_final.csv')\n",
"correct_file = os.path.join(base, data_only, 'NH_seaice_extent_final.csv')"
],
"language": "python",
"outputs": [],
"prompt_number": 65
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def parse_the_date(year, mm, dd):\n",
" return dt.date(int(year), int(mm), int(dd))\n"
],
"language": "python",
"outputs": [],
"prompt_number": 66
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"current_data = pd.read_csv(current_file, header = None, skiprows=2, names=[\"year\", \"mm\", \"dd\", \"extent\", \"missing\", \"source\"],\n",
" parse_dates={'date':['year', 'mm', 'dd']},\n",
" date_parser=parse_the_date, index_col='date')\n",
"\n",
"correct_data = pd.read_csv(correct_file, header = None, skiprows=2, names=[\"year\", \"mm\", \"dd\", \"extent\", \"missing\", \"source\"],\n",
" parse_dates={'date':['year', 'mm', 'dd']},\n",
" date_parser=parse_the_date, index_col='date')\n"
],
"language": "python",
"outputs": [],
"prompt_number": 67
},
{
"cell_type": "markdown",
"source": [
"You can see that there are a number of days where this happens in 2014, but I have looked and there are also days in the final data."
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" correct_data.rename(columns={'extent':'correct-extent', 'source':'correct-source'}, inplace=True)\n",
" current_data.rename(columns={'extent':'current-extent', 'source':'current-source'}, inplace=True)"
],
"language": "python",
"outputs": [],
"prompt_number": 68
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"merged = pd.merge(correct_data, current_data, left_index=True, right_index=True, how='outer')"
],
"language": "python",
"outputs": [],
"prompt_number": 69
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"merged['diff'] = merged['correct-extent'] - merged['current-extent']"
],
"language": "python",
"outputs": [],
"prompt_number": 71
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"with mpl.rc_context(rc={'figure.figsize': (17,5), 'font.weight': 'bold', 'font.size': 8}):\n",
" ax = merged[['diff']].plot(title=\"Extent difference\")\n",
" ax.set_ylabel(\"Mkm^2\")\n"
],
"language": "python",
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAA/UAAAFRCAYAAAA1qY0oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYXGWZ9/HvHRJAJRBwBRQTEEZBBRRFx0GDOIMOIrgw\nMjqGBEcFESeKiiwvKIvAMMMyLrghIDo4IiKboIAEd53RiSKOImBAwyIoASSAQO73j3OKdDpVle6u\nqq56ur6f66qr+ix1zlP9q+qup85znxOZiSRJkiRJKs+0fjdAkiRJkiRNjJ16SZIkSZIKZadekiRJ\nkqRC2amXJEmSJKlQduolSZIkSSqUnXpJkiRJkgplp16SJE1IRJwZEUfXP+8UEb8aseyvImJxRNwT\nEe+KiHUj4qKIWBYR/9W/VkuSNLXYqZckqY8iYklELI+Ie0fc/mMMj5sbEb/rYjs+FBFnj/NhWd/I\nzO9k5jNHLPsAcGVmrp+ZHwP2Ap4EbJSZb+xKoyVJEtP73QBJkoZcAq/OzG/1uyETFC3mPx34/qjp\n6zJzxbh3EDE9Mx+eSOMkSZrqPFIvSdKAiojTIuIrI6ZPiIgrIuKxwKXAJvWR/Xsi4ilR+WBEXB8R\nd0bEf0XEhvVjZ0fEioiYFxE3RcQdEXFoveyVwCHAG+vt/W+L9mwfET+t9/clYN0Ryx4dORAR3wLm\nAh+rt/efwP8bsf0F9Xr7RsQvI+JPEXFZRGw2YnsrIuKdEfEb4Nf1vFfXQ/rviojvRcRzRqy/JCIO\nioif1UP8vxQR64xYvkf92Lvr38+u9fwNIuL0iLglIn4fEUdHhJ+PJEnF8J+WJEn91+po93uB50TE\nPhGxE7AvMC8zlwOvBG7JzJn1EPfbgHcDrwFeCmwM3AV8fNQ2XwJsBewCHBERf5WZlwEfAb5Ub2/7\n1RoYsTbwNeAsYEPgXOD11MPvR8rMlwPfAQ6ot/emUds/IyL2oPoi4bXAE+r1zxm1qT2AFwBbR8T2\nwOnA24CNgE8BF0bEjMZuqYb47wrMAZ4LzK/b/sK63Qdl5gb172dJ/bgzgb8AWwDbA38H/PPo5yRJ\n0qCyUy9JUn8F8LX66HPj9laAzLwfeAtwMnA28K7MvGXE40Z7B3B4Zt6SmQ8BHwbeMOrI84cz88HM\n/DnwM2DbEdtr9eUCwIuA6Zl5amY+kpnnAf89huc28ueR0/sBx2Xmr+sh+ccB20XE00asc1xmLsvM\nB4G3A5/KzP/OyueBB+t2NfxHZt6WmXcBFwHb1fPfCpyemVcC1L+fX0fEk4FXAe/JzPsz8w7gFGDv\nNTwvSZIGhjX1kiT1VwJ7tKqpz8wfR8SNVEezz13DtmYD50fEyLr1h4Enj5i+bcTPy4H1xtjOTYCl\no+bdtIbHrHYUf4SnA6dGxL+Pmr8p0DgB4O9GrT8vIg4cMW9G3a6Gkc/tfqrRCgBPBS5p0YYZwK0R\nj37fMA24uU27JUkaKB6plyRpgEXEAcDawC1UZ5RvaNZhvhl4ZWZuOOL22My8dQy7atcBB7iVqsM9\n0tPHsN1W278ZePuotj4uM3/Y4jE3A8eOWn+9zBzL5fF+BzyjxfwHgceP2OYGmfmcJutKkjSQ7NRL\nktR/TYe9R8RWwNHAm4F5wAciojFc/nbg8RGx/oiHfBL4SOOEcxHxxIh4zRjbcBswO0Ycsh7l+8DD\nEfHuiJgREa+jqndvZ/Tw+5E+CRwaEVvXbd0gIvZqs63PAPtFxAvrEwI+LiJ2i4h2Iw0a+zwdWBAR\nL4+IaRGxaX0ugVuBbwInRcTMetkWEfHSNTwvSZIGhp16SZL676JY9Tr150XEWlR19Mdn5jWZeT1w\nKHB2RMzIzF9RnVjuxvrs8U8BTgUuBL4ZEfcAPwBeOGI/7Y7GN4b2/zEi/mf0wrpG/3VUJ5/7I/AP\nwHmjV2sznSOnM/NrwAnAlyLibuAaqpPcNd1WZv6E6iR5HwP+BPyG6ouOVs/p0f1l5n8DC6jOTbAM\nWAQ0zrQ/j2okxC/r7Z4LPKXFNiVJGjiRuabRdj3cecRLgNOozsJ7LfDPmbnaZXQiYk/g36iG/f0Q\nWJCZS9a0TJIkSZKkqaxvR+ojYl2qb/gfByykOonPV0ZfG7Y+8vAlqm/W3wc8n+qyNG2XSZIkSZI0\n1fVz+P2rgCcBn8jMT1LVu80B5o5a7x+phsUdl5kfp7pG7k4RsfkalkmSJEmSNKX1s1M/p75fOup+\nzhrW+319v3mbZaO3IUmSJEnSlDNIJ8prdbbd8aw31m1IkiRJklS86X3c9431/dPq+8a1b2+s6+0f\nqc+0O3K9H45Y74Y2yxrzHxUR/TsjoCRJkiRJHcrM1Q5k9/NI/aXAH4D9I2J/4K3Ab4GrgeXAT+v1\nvgT8BTg4Ig4EXgt8JzN/u4Zlq8nMVW5HHnnkavO8De7NvMq6TTSvX/yiugrVihX9fw7HHlu1JTNZ\ntqxxdazOtwvJeuut3F6rbe6+e/NlkMyY0Xr7G2yw6uMg2XTT7mfV6e9g9uzebHeLLTp7/Iknrj7v\nec9bdd51162a37/+a+tt/uEP1TpLl66c96c/VfNuuKF1OyC5+upV53/xiwl0ntd5563a/t/+dnLz\nH6ab/7vKuplXf28PPlj9TVqyZPLyOvHE5L3vTfbee83/X711dutGXt36LLam2xFHVPs59tjkwAN7\nt5+77qr284QnjO35tdK3Tn1mPgjsBfwZOAW4DdgrM1c0VqnXu43qhHizgBOBn1BdI7ftsrFYsmRJ\nx89Dk8e8yjLRvB55pLqPASimefjhlT9P7/K4pmlj+Os7cv+jtWtPs8e1W79f762x/A4motOsmj1+\ndFsbr9OGdq/XxmNHbrex/lprja9tVbZLxveglttZad11O96kWvB/V1nMq78afyfH+jepG3lNn179\nTWzsu1f/m1TW+2uyXgeNzwNt+utj0s/h92Tmd4DnNpk/bdT0+cD5LbbRcpmk8ozuLPXTyLbMmNHd\nbY+lM9fud9GuPc0e1+32d8N4O7Rj1elzbfb40W2dSKd+5HbH2qkfvd1uvT9Gb2cQXx+Shk/j7+Xa\na0/ePmfMqDr1jb/Hg3BQQf3Xrc72WPfTqb526vtt/vz5/W6CxsG8yjLRvJ785O62oxM77ggbblj9\n3PgGv1svw7F8A7zbbrBs2erz994bnvSk1o+bPx/uv3/l9N//PTz/+e3Wn7/mxvRAL74Ff8Ur4KUv\n7Wwbz3ve6vNGt/UJT1h1erxH6hvzWnXq99oLzj0XNtts1fnPfjaMY0BaS8961qrTj3tcx5tUC/7v\nKot5DYbHPnZs63Ujr8aRejv1vVfi++vcc+FlL+vd9rv1eot2Y/OnkojI0c81huRdOywZS6WIgCc+\nEf7wh8n7JnjQRFQdy1/+st8tWbMIeMlL4Lvfbb38pJPgPe9pvvzee2H99WH5cnjMY6p5f/4zzJwJ\nt902WF9kSdKwOf10+P73q079Zz5TfaFw3339bpX67aij4Mgjq9fFAQfAqaf2Zj+NzwgbbQR//OOa\nPxdGBDlgJ8obCP0+YUSvb1PJokWL+t0EjYN5tTfyyG+/v1/sV1Yl1S2uqa2TVVMPvrdKY15lMa+y\ndCOvxpH6xt/qfv9PnspKfH+tWLHmdTrRrddbQR+pJGnqGKROfb8MS6e+WQd+TcPvJUmTw0692sns\n7WuiWyM2C/pIpWE3d+7cfjdB42Be7Q1Sp75fWQ1Lp77ZNjo5Uu97qyzmVRbzKks38po+HRYt6v//\n4mHg+2t1jdfdXXd1tp2CPlJNbfPnz2fatGm85z3v4cwzz2TatGnsvPPOjy5ftGgRW2+9NdOnT2fa\ntGksX76c3/3ud+y888485jGPYdq0aVxyySV9fAaSxuqii+DLX145PYwfJM47Dz772X63YmzOOQdO\nOaX9Ou0yXHddOPTQ5uuX9MWGJE1F06fDLbes7FRNsepVTdB++03Ofhx+P0VFBNtssw0LFy5kr732\nenT+O9/5Tn71q1+x6667snDhQmbMmMFxxx3H1VdfzVZbbcXChQvZYost+tjy3iuxDmeYmVdrr341\n/PVfr5zud6e+H1m97nWwww6TvtsJ2XtveO5qF19d1Zpq6o89dvV5YE39MDCvsphXWbpVUz+Snfre\nKen9NfIqQ5Mx/L5TQ31Ju0H1ghe8gBe84AWrzLvuuuuICD72sY8xZ86cR+cBLFy4kAULFkx6OyV1\nR7879erceDPsZPi9JKl7Rv8d7vWJ0aRe8Eh9n3z3u9/lOc95Duuttx777LMPDzzwwKPLRg+/nzZt\nGitWrCAz2WKLLZgzZw4777wz3/rWtwB461vfyrRp07j55pv78lwmi3U4ZTGvset3p96sOjeZnXrz\nKot5lcW8ytKNvBp/jxudeTv1veP7a3UeqS/YsmXL2H333bn77rvZZZdduPPOO/nGN74BVJfYizrd\nxv2//Mu/cGp9ccR9992X2bNns+GGG3L99dezdOlS/u7v/o6tt96amTNn9ucJSeqIddXlG2+GjfXN\nXpL66+GHq/tHHqnun/a0/rVFg6mE4fd+nOiDiy++mLvvvpstt9ySyy+/nEsuuYTtttsOWNmRH+nk\nk09+dNkRRxzB4YcfzgEHHPBoDf2b3vQmTjrpJDbccMPJexJ9UFIdjsxrPPp9pN6sOjeRI/U33zyx\n7M2rLOZVFvMqSzfyevzjq/tG5/7sszvepFrw/bU6O/U9FDHx21gsXboUgC233PLReSN/ljRc+t2p\nV+cmkqFHgySp/2bMqO4bR+of85j+tUXDx+H3PdTrs14+9alPBVae6G70z2rOOpyymNfY9btTb1ad\nm8wMzass5lUW8ypLN/NqdOqtqe+dUt9fJQy/t1PfB7vtthsbbLAB119/PbvssgvrrLMOixcvBqqa\n+lbaLZNUrn536tU5M5Sksv3lL9W9nXqVyOH3fTBr1iwuvPBCttlmG370ox8xa9YsXv/61xMRTWvq\nG0YvW9P6U411OGUxr7F5+9vhgAP62waz6sysWfD850/e/syrLOZVFvMqSzfy2myz6v6KK6r7ekCt\nesD31+o8Ul+4nXbaiWuuuabl8n322WeV6RVNvja86qqrut4uSZPrU5/qdwvUqbvu6ncLJEkTtfHG\nK3/+5CfhKU/pX1s0fLrVqY9hGdIdETn6uUbElB/SPgzPUZIkSZqoRsfqU5+qRtBJsPJ1cdBB8G//\n1vv9ZK76c/N1g8xc7asAh99LkiRJGnpDVNWqKcZOvYphHU5ZzKscZlUW8yqLeZXFvMpiXmUxr/Y6\nGVxtp16SJEnS0PNIvZrp9euicbLGRqf+nHPGvw1r6qf48x+G5yhJkiRNVKPT9pnPwD//c3/bosHR\neF28731w4om928+//iscfDA89BA885lw2WXwjGe0alPzmvqhP/v9MF0STpIkSVJzdgvUTytWTHwI\n/lAPv7/qqqvIzCl/myqswymLeZXDrMpiXmUxr7KYV1nMqyyl5jVZX/Y0rmA+kf0NdadekiRJksAj\n9eqPH/+4un/kkYkfqR/qmnpJkiRJw63RmT/9dNh33/62RYOj8bp4//uruvdeWXvtqp7+3nvhuc+F\nK66AzTdv1SavUy9JkiRJTXmkXs30+nXRGHafaU39hJRa1zGszKss5lUOsyqLeZXFvMpiXmUxr7KY\nV3OPPFLdNzr01tRLkiRJ0gR4pF79Zk39GlhTL0mSJGm0Rmf+jDNg/vy+NkUDpPG6OPhgOP743u3n\nDW+A886DZctg221h0SKYPbtVm6yplyRJkqSmPFKvfthzz+remvoJsq6jLOZVFvMqh1mVxbzKYl5l\nMa+ymFdZzKs9a+olSZIkqQMeqVczk/W66ORIvTX1kiRJkoZWo9N21lkwb15/26LB0XhdfPCDcNxx\nvdvPF74Ab3kL3HknbL89fPe7sNlmrdpkTb0kSZIkNeWRevXDXntV99bUT5B1HWUxr7KYVznMqizm\nVRbzKot5lcW8ylJqXr3+smeddeDxj7emXpIkSZI64pF69UuENfVjYk29JEmSpNEanfmzz4Z/+qf+\ntkWDo/G6OOQQ+MhHeruvJz0JrrkGnv98+OEP4alPbdUma+olSZIkqSmP1KuZyXhdNPZRVE19RMyM\niHMi4r6IuDUiDmqz7iYRcWlE3B8RSyLizSOWLYiI30TEioi4d7ztKLWuY1iZV1nMqxxmVRbzKot5\nlcW8ytLtvOzU95bvr/ZKrKk/BngjcALwA+DEiNi5xbqnATsDhwA3AWdFxBb1snWAC4A/AY6tlyRJ\nkiQVpcia+ohYBvwuM58TEXOAG4AvZOa8UettBNwJXJKZu9cd/yuBYzLziBHrLQE2ysz12+zTmnpJ\nkiRJq2gcGf3P/4R//Mf+tkWDo/G6OPRQOPbY3u5r443hJz+BHXaA//kf2GSTVm0akJr6uqO+PrC0\nntW4n9Nk9dmj1mm3riRJkiRJXTMZZRn33AMPPTTxI/XTu9uclSLi90Cz7xhG18+P59fU0a90/vz5\nzJ49G4BZs2YBsHDhQmBljcfcuXOdHtDpxYsXm1dB0+ZVzvQpp5zCdtttNzDtcbr9tHmVNW1eZU2b\nV1nT3crrFa+YyxVXwP/93yIWLRqc5zfVpkt7f8Gi+r73+1u+HBYsWMSDD0LEyuWLFy9m2bJlACxZ\nsoRWejb8vh5WP6PJotuAm4GlmblNRGwOXA+cnZn7RMR0YHpmPlAf1b8DuDQzXx0RLweuAI7OzCNH\n7GsJExh+v2jRohGhadCZV1nMqxxmVRbzKot5lcW8ytKtvH7xC3jOc+Ccc2DvvTtvl5or7f3VOEJ/\n2GFwzDG939cLXwg33QT/+7/VcPzm6zUfft+vmvpTgQOBDwHbA3sAu2TmVRFxJjAP2CEzfxoRFwC7\nAocCewIvBp6ZmTdExPOA5wHHA48D3gVcl5nfabJPa+olSZIkreKXv4RttoEvfQne+MZ+t0aDol+d\n+sWL4SlPabVe8059z4bfr8HhwJOBg4F7gIMz86p6WbLqmez3Bz4HHAvcDizIzBvqZbsDR45Y/zPA\nmcBqnXpJkiRJGs1L2amdyTou3DgD/kRM625TxiYz783MvTPzcZm5cWaeOGLZgsxcKzN/Wk/fkpmv\nzMzHZObszPzCiHU/nJnT6vXXqn/ed6ztaNQxqAzmVRbzKodZlcW8ymJeZTGvsnQrr2l1j8jOfW+V\n+v6arE5945J2JV2nXpIkSZL6zs68mrnssup+siu4J/J67EtNfT9YUy9JkiRptN/8BrbaCr78Zdhr\nr363RoMkAg4+GI4/vvf72XFHuOEGuPZaeNKTWq03INeplyRJkqRBMc0ekdoo4Uj9UL+ES63rGFbm\nVRbzKodZlcW8ymJeZTGvsnQrr0YnymH4vVXq+2syO/VFnShPkiRJkgaBnXm1U8KRemvqJUmSJA2t\nm26C2bPhK1+B17++363RIImA970PTjxxzet2up8dd6zO7/DrX8MTntBqPWvqJUmSJGkVHqlXOyUc\nqR/qTn2pdR3DyrzKYl7lMKuymFdZzKss5lUWa+rLUur7a8WKyduXNfWSJEmSNE6e/V6DxJr6Nqyp\nlyRJkjTaLbfAppvCV78Kr31tv1ujQRIBCxfCySf3fj877ljV099wA2y0Uav1rKmXJEmSpFV4pF7t\nWFM/4Eqt6xhW5lUW8yqHWZXFvMpiXmUxr7JYU1+WUt9fXqdekiRJkgaYnXq1U8KRemvqJUmSJA2t\nO++EJz4RvvY12GOPfrdGgyQC3vUu+OhHe7+fHXeE//s/uOkmmDWr1XrW1EuSJEnSKjxSr0HQyetw\nqDv1pdZ1DCvzKot5lcOsymJeZTGvsphXWbpdU6/eKvX9NVmDvTOtqZckSZKkCbNzr2ZWrJjc/VlT\n34Y19ZIkSZJGu+uu6rrgF14Iu+/e79ZokETAfvvBaaf1fj8vehH84hewdCmsv36r9ayplyRJkqRV\nWFOvQbD55tW9NfXjVGpdx7Ayr7KYVznMqizmVRbzKot5lcWa+rKU+v6ajMHeb387vPSl1tRLkiRJ\n0rh5pF7tTFZNfaNDb019G9bUS5IkSRrt3nurGuaLL4bddut3azRIIuBtb4NPf7q3+9lvP9h2W3jf\n++D222G99Vq1x5p6SZIkSVqFR+rVzmQcFx752rOmfpxKresYVuZVFvMqh1mVxbzKYl5lMa+yWFNf\nllLfX16nXpIkSZIGmEfq1c5kHam3pn4MrKmXJEmSNNr998NjHwtf/zq86lX9bo0GSQQsWACf+1xv\n93PAAbD11lVN/R//WL0em7fHmnpJkiRJWoVH6tXOZA6/B2vqx63Uuo5hZV5lMa9ymFVZzKss5lUW\n8yqLNfVlKfX9NZnD762plyRJkqRx8ki92inh7PfW1EuSJEkaWg89BGuvDZddBrvu2u/WaJBEwFve\nAp//fG/3c+CBsOWW8P73w913w7rrtmqPNfWSJEmStAqP1KvfRg6/t6Z+nEqt6xhW5lUW8yqHWZXF\nvMpiXmUxr7J0u6beTn1vlfr+KmH4/VB36iVJkiQNNzvzamfFisnZTycnyrOmXpIkSdLQyoRp0+Dy\ny+EVr+h3azRIIuBNb4IvfrG3+1m4EJ7+dPjAB2D5cpgxo1V7rKmXJEmSJGnMvKTdgCu1rmNYmVdZ\nzKscZlUW8yqLeZXFvMriderLUur7y5p6SZIkSZIKNVkV3EXW1EfETODTwGuAe4B/y8x/b7HuJsDp\nwFzgduCwzPxiveyrwAuBjYDfAO/LzMubbMOaekmSJEmriYBvfhP+9m/73RINkgj4h3+A//qv3u7n\noINg443h4IPhL3+BtdZq1Z7Bq6k/BngjcALwA+DEiNi5xbqnATsDhwA3AWdFxBb1sm2BjwKHAc8A\nzouIx/ay4ZIkSZKkqW8yh9+XWFO/D3BtZh4FHFTPWzB6pYjYCNgduDwzTwE+RNXufepVnpmZJ2Tm\nycBlwHrAZmNpQKl1HcPKvMpiXuUwq7KYV1nMqyzmVZZu5+Wg3t4q9f01mcPvoaCa+rqjvj6wtJ7V\nuJ/TZPXZo9ZZZd3MfKje5gbAjsAtwPXdbbEkSZIkadgM/dnvI+L3EbFi9I2VR9kfXXU8m22yn/WA\nC6nq6t+UmQ+PZUNz584dx27Vb+ZVFvMqh1mVxbzKYl5lMa+ydDsvj9T3VqnvrxLOfj+9e01paidg\nRpP5twFHAk+rpzet728EiIjpwPTMfABYAiTw1BbrzgQuBZ4HvC4zv92qMfPnz2f27NkAzJo1i+22\n2+7RF1djOIjTTjvttNNOO+200047PXzTP/vZItZZZ3Da4/RgTGdOzv6uv76ahpXLFy9ezLJlywBY\nsmQJrfTz7PenAgdS1chvD+wB7JKZV0XEmcA8YIfM/GlEXADsChwK7Am8mKqW/oaI+AHVsPszgW/W\nm78yM+8Ytb/Vzn6/aNGiR3+JGnzmVRbzKodZlcW8ymJeZTGvsnQzrwi49FJ45Su7sjk1UeL7KwL2\n3BPOP7+3+zn4YNhwQzjkkPYjA1qd/b7XR+rbORx4MnAw1SXtDs7Mq+plWd8a9gc+BxxLdUm7BZl5\nQ71sx3rd+fUtqc6Uv0qnXpIkSZKk8ZjMmvoJP35Yrt3udeolSZIkNRMBX/86vOpV/W6JBkkEvOY1\ncMEFvd3PBz8I668Phx02sSP103rZOEmSJEmS1FqnR+qHulPfODmBymBeZTGvcphVWcyrLOZVFvMq\nS7fzclBvb5X6/prM4fcTOfM9DHmnXpIkSZKkVlasmJz9WFM/BtbUS5IkSWomAi6+GHbbrd8t0SCJ\nqF4TF1/c2/0cdhisvTYcdRQ88ki79lhTL0mSJEnSmJVw9vuh7tSXWtcxrMyrLOZVDrMqi3mVxbzK\nYl5l6XZeE61n1tiU+v6arMHe1tRLkiRJUges1FUzJRypt6ZekiRJ0lCzpl7NRMCuu8Jll/V2P0cc\nUXXqjz8eHnqoXXusqZckSZKkpjz+p2Ymc/j9RA11p77Uuo5hZV5lMa9ymFVZzKss5lUW8yqL16kv\ni++v1iKqS+dZUy9JkiRJE+SJ8tSMNfUDxJp6SZIkSc1EwEUXwatf3e+WaJBEwC67wBVX9HY/H/4w\n3H8/nHwyPPhgu/ZYUy9JkiRJ0sCxpn6CrOsoi3mVxbzKYVZlMa+ymFdZzKss5lWWUvOarOH31tRL\nkiRJktRl1tQPEGvqJUmSJDVjTb2aiYC5c+Gqq3q7n6OPhmXL4BOfqGrrW7fHmnpJkiRJksbM4fcD\nrtS6jmFlXmUxr3KYVVnMqyzmVRbzKot5laXUvEoYfj/UnXpJkiRJklqZrAruzIkfqbemXpIkSdJQ\ns6ZezUTATjvBt7/d2/185CNw221w+ulw333t2mNNvSRJkiRJY7bWWr3fR2P4vTX1E1BqXcewMq+y\nmFc5zKos5lUW8yqLeZXFvMpSal7Tp0/OfqyplyRJkiSpyybrSH0nZ7+3pl6SJEnSUIuAiy+G3Xbr\nd0s0SCLgne+Ej3+8t/s54QRYsgS+8AW499527WleUz9JgwkkSZIkSSrHsmXw2MdOzr56UlMfERtE\nxPER8YWIeNOoZZ+Y2O4GS6l1HcPKvMpiXuUwq7KYV1nMqyzmVRbzKkuJeW2wAcyY0fv99PI69WfU\n9+cB/xgR50XEuvW8F098l5IkSZIkCTo/+33LmvqI+Flmbjti+jDg74E9gMszc/uJ7bI/rKmXJEmS\n1Iw19eqnE0+E666DL38Z7r679XoTqalfOyKmZeYKgMw8NiKWAlcD63XYbkmSJEkaCE95CmyzTb9b\noWHV6dnv2w2/vxjYZeSMzDwTOAj4y8R2N1hKrOsYZuZVFvMqh1mVxbzKYl5lMa+ydDOvW2+F2bO7\ntjk14furtU5r6lseqc/M97eYfxmw5cR3KUmSJEmSGnpSU//oChEzM7PN1fLKYE29JEmSJGnQnHQS\n/PzncMEFcNddrddrVVPfbvg9EbEpcEnHrZQkSZIkSauJgHPPnfjj212nfhvgG8AHJr75wWZdR1nM\nqyzmVQ6zKot5lcW8ymJeZTGvsphXaxGwfHl1sryJaHf2+0XAnpn5w4ltWpIkSZIktdOopX/kkQk+\nvs116i8GU2WIAAAZ+UlEQVQBftnqhHmlsaZekiRJkjRoTj0VFi6EddaBBx5ovd5Eaur3AGZFxL92\n3kxJkiRJktTKRI/Ut+zUZ+bDmfk24M8TbNPAs66jLOZVFvMqh1mVxbzKYl5lMa+ymFdZzKu1xvD7\nhx+e2OPbnv0eIDOPmtim24uImRFxTkTcFxG3RsRBbdbdJCIujYj7I2JJRLx5xLJTI+KOiHggIq6P\niP160V5JkiRJkrptotenf/TxY6kzj4gNgc2AtRrzMvOnHe044lTgQOBDwHbAnsAumXlVk3UvAHYF\nPgi8FngJ8FeZeUNEvAe4D0jgSGATYNPMvHXUNqyplyRJkiQNlI9+FN797urndl3WVjX17c5+33jg\n0cB84EZg5En2dx5XS1e3D3BtZh4VEXOoOvULgFU69RGxEbA7cElmnhIRPwOurB9/RGaeHBEbAE8E\n9gOe3GG7JEmSJEkqwhqH3wNvBLbIzJdl5s6NWyc7rTvq6wNL61mN+zlNVp89ap1m614NXAdsC7xn\n9FH6VqzrKIt5lcW8ymFWZTGvsphXWcyrLOZVFvNqrdPh92s8Ug/8AtgQuH28G4+I31MNhx9tdP38\neJ5Gs3XfBjwDOBo4LCL+KzPvGL3S/PnzmT17NgCzZs0CYO7cucDKF5nTgzu9ePHigWqP0+Y1VaYX\nL148UO1xuv20eZU1bV5lTZtXWdPmVda0ebWfhkX1/crlixcvZtmyZQAsWbKEVtZYUx8ROwAXANcC\nD9azMzNf0/aB1WPnADOaLLoNuBlYmpnbRMTmwPXA2Zm5T0RMB6Zn5gP1Uf07gEsz89UR8XLgCuDo\nzDxy1P4+DPw/4A2Z+dVRy6yplyRJkiQNlI9/HN71rurnntTUA58HTqA6Yt+oqR9T7zgzf9umQWcB\nB0bEEcD29ewz6/vPAvMiYofM/GlEXAzsGhHvpaq9fwT4fF1L/yXgIqpSgrfVbfzVWNonSZIkSVLJ\npo1hnfsy8z8y81uZuai+Xd2FfR8OfBk4GHgRcPCIM98nq35xsD/VeIRjqc7CvyAzbwD+AswEjqP6\n4uEPwFsy85djaUBjyIPKYF5lMa9ymFVZzKss5lUW8yqLeZXFvFqbjJr670TEccCFrBx+3/El7TLz\nXmDvFssWUJ0JvzF9C/DKJuvdD/xNJ+2QJEmSJKlUY6mpX0ST4fadngF/sllTL0mSJEkaNJ/4BBxw\nQPVzr2rqX5mZD4za2OPH1UpJkiRJkrSaToffj6Wm/qsR8egZ7CNiY+DyznY7GKzrKIt5lcW8ymFW\nZTGvsphXWcyrLOZVFvNqbTI69ecDX46ItSJiNvAN4IOd7VaSJEmSJHVqjTX1ABHxLqoT1T0d2C8z\nv9frhnWbNfWSJEmSpEHzyU/C/vtXP3e1pj4iDqp/TCCApwE/A14UETtm5kkTbrUkSZIkSerp8PuZ\nwHoj7s8HflNPz+xst4PBuo6ymFdZzKscZlUW8yqLeZXFvMpiXmUxr9Z6dp36zPxQZ5uWJEmSJEm9\n1LKmPiIuYuXQ+9EyM1/Ty4Z1mzX1kiRJkqRB8+lPwzveUf3c7evUvwj4PXAO8KPGdup7e8eSJEmS\nJPVZu5r6jYFDgWcDpwB/C9yRmYsy8+rJaFyvWddRFvMqi3mVw6zKYl5lMa+ymFdZzKss5tVaz06U\nl5kPZ+almTmP6qj99cDV9eXtJEmSJElShzrt1Le9Tn1ErAvsBuwNzAYuBD6XmUs72+3ks6ZekiRJ\nkjRoPvtZeNvbqp+7fZ36s4FtgK8DR2XmNR22VZIkSZIkdVG7mvo3A1sC/wJ8PyLuHXG7Z3Ka11vW\ndZTFvMpiXuUwq7KYV1nMqyzmVRbzKot5tdbL69S36/BLkiRJkqQ+a1tTP5VYUy9JkiRJGjSf+xy8\n9a3VzxOpqfdovCRJkiRJhRrqTr11HWUxr7KYVznMqizmVRbzKot5lcW8ymJerfXsOvWSJEmSJGmw\nWVMvSZIkSVKfnHEG7Ltv9bM19ZIkSZIkFcTh9x2wrqMs5lUW8yqHWZXFvMpiXmUxr7KYV1nMqzU7\n9ZIkSZIkDSlr6iVJkiRJ6pOzzoL586uframXJEmSJKkgDr/vgHUdZTGvsphXOcyqLOZVFvMqi3mV\nxbzKYl6t2amXJEmSJGlIWVMvSZIkSVKfnH02zJtX/WxNvSRJkiRJBXH4fQes6yiLeZXFvMphVmUx\nr7KYV1nMqyzmVRbz6p2h7tRLkiRJktRPnR6pt6ZekiRJkqQ++c//hDe/ufrZmnpJkiRJkgqy1lqd\nPX6oO/XWdZTFvMpiXuUwq7KYV1nMqyzmVRbzKot5tWanXpIkSZKkQnXaqbemXpIkSZKkPrnoInjN\na6qfi6mpj4iZEXFORNwXEbdGxEFt1t0kIi6NiPsjYklEvLnJOp+PiBURcU1vWy5JkiRJUvfstFNn\nj+/X8PtjgDcCJwA/AE6MiJ1brHsasDNwCHATcFZEbNFYGBF/B7y2nhzXoXjrOspiXmUxr3KYVVnM\nqyzmVRbzKot5lcW8Wps1CzbbbOKP71enfh/g2sw8CmgcpV8weqWI2AjYHbg8M08BPkTV5n3q5Y8F\nPgkcPgltliRJkiSp6zq5Vv2k19TXHfU7gW9m5isjYm3gAeB7mbnTqHWfB/wP8OnM3C8itgJ+BXwx\nM98SEf8O/HV9ewT4RWY+t8V+ramXJEmSJA2cOXNgyZIBq6mPiN/Xde6r3KiPso9cdTybHbH9rYB3\nA8cDjeH4a0fE0zpquCRJkiRJhZjew23vBMxoMv824Eig0fnetL6/ESAipgPTM/MBYAlVnfxTm6y7\nMbAWcP6IbW8FXFnfr2b+/PnMnj0bgFmzZgGwcOFCYGWNx9y5c50e0OnFixebV0HT5lXO9CmnnMJ2\n2203MO1xuv20eZU1bV5lTZtXWdPmVda0ebWffuCBahpWLl+8eDHLli0DYMmSJbTSl0vaRcSpwIFU\nNfLbA3sAu2TmVRFxJjAP2CEzfxoRFwC7AocCewIvBp4J3A28jKrTH8C5wO+At2XmN5vsc7Xh94sW\nLXr0l6jBZ15lMa9ymFVZzKss5lUW8yqLeZXFvNrbYgu48caJDb/vV6d+JvAZqpPg3QOclJkn1svO\noOrUv6Du1G8CfI6qA387cHhmfqHJNldgTb0kSZIkqTDPeAbccENBnfp+sFMvSZIkSRpEnXTqp/Wy\nYYOuUcegMphXWcyrHGZVFvMqi3mVxbzKYl5lMa/2Ormk3VB36iVJkiRJKpnD7yVJkiRJ6qOttoLf\n/Mbh95IkSZIkFcfh9xNkXUdZzKss5lUOsyqLeZXFvMpiXmUxr7KYV3t26iVJkiRJGkLW1EuSJEmS\n1EfPehb86lfW1EuSJEmSVByH30+QdR1lMa+ymFc5zKos5lUW8yqLeZXFvMpiXu3ZqZckSZIkaQhZ\nUy9JkiRJUh89+9lw7bXW1EuSJEmSNFSGulNvXUdZzKss5lUOsyqLeZXFvMpiXmUxr7KYV3vW1EuS\nJEmSVKhOOvXW1EuSJEmS1Efbbgs//7k19ZIkSZIkDZWh7tRb11EW8yqLeZXDrMpiXmUxr7KYV1nM\nqyzm1Z419ZIkSZIkFcqa+jGwpl6SJEmSNIi23x4WL7amXpIkSZKk4jj8foKs6yiLeZXFvMphVmUx\nr7KYV1nMqyzmVRbz6p2h7tRLkiRJktRv1tSPgTX1kiRJkqRBtMMO8JOfWFMvSZIkSVJxrKmfIOs6\nymJeZTGvcphVWcyrLOZVFvMqi3mVxbx6Z6g79ZIkSZIk9Zs19WNgTb0kSZIkaRDtuCP8+MfW1EuS\nJEmSVBxr6ifIuo6ymFdZzKscZlUW8yqLeZXFvMpiXmUxr94Z6k69JEmSJEn9Zk39GFhTL0mSJEka\nRC9+Mfzwh9bUS5IkSZJUHGvqJ8i6jrKYV1nMqxxmVRbzKot5lcW8ymJeZTGv9uzUS5IkSZJUKGvq\nx8CaekmSJEnSIPqbv4Hvfc+aekmSJEmSiuPw+wmyrqMs5lUW8yqHWZXFvMpiXmUxr7KYV1nMqz07\n9ZIkSZIkFaq4mvqImAl8GngNcA/wb5n57y3W3QQ4HZgL3A4clplfrJedCcwb9ZBZmXlPk+1YUy9J\nkiRJGjgvexl8+9tl1dQfA7wROAH4AXBiROzcYt3TgJ2BQ4CbgLMiYosRy+8A9h5xW96rRkuSJEmS\n1G0lDr/fB7g2M48CDqrnLRi9UkRsBOwOXJ6ZpwAfomrzyKPzy4FLgHMz88uZ+fBYG2FdR1nMqyzm\nVQ6zKot5lcW8ymJeZTGvsphX70x6p77uqK8PLK1nNe7nNFl99qh1Gvebj1hnM+Be4L6I+GhEJ99x\nSJIkSZI0uTrpxU7vXjNWFRG/BzZpsuig0auOZ7Ojpq+sb3fV2z0A+BHwhWYPnj9/PrNnzwZg1qxZ\nbLfddo8ua3xzNHfuXKcHeLphUNrjdPvphkFpj9PNpxvzBqU9TrefbswblPY43X66MW9Q2uN0++nG\nvEFpj9PtpxvzBqU9TrefbswblPYM2vSyZdU0rFy+ePFili1bBsCSJUtopWcnyouIOcCMJotuA24G\nlmbmNhGxOXA9cHZm7hMR04HpmflAfVT/DuDSzHx1RLwcuAI4OjOPHLW/3YELgOMy87Am7fFEeZIk\nSZKkgbPLLvCtbw3YifIy87eZeV2T2z3AWcCzIuIIoHHW+zPr+88CyyPieZn5J+Bi4BUR8V6qmvpH\ngM/XT+qqiHh/ROwLfBhIqiP1Y9L4dkRlMK+ymFc5zKos5lUW8yqLeZXFvMpiXr3Ts+H3a3A48GTg\nYKpL2h2cmVfVy7K+NewPfA44luqSdgsy84Z62XXAgcATgd8D783MC3vffEmSJEmSuqO469T3g8Pv\nJUmSJEmD6BWvgCuvHLDh95IkSZIkac1KvE79QLCuoyzmVRbzKodZlcW8ymJeZTGvsphXWcyrPTv1\nkiRJkiQNIWvqJUmSJEnqo113hW9+05p6SZIkSZKK4/D7CbKuoyzmVRbzKodZlcW8ymJeZTGvsphX\nWcyrPTv1kiRJkiQNIWvqJUmSJEnqo7//e7j0UmvqJUmSJEkaKkPdqbeuoyzmVRbzKodZlcW8ymJe\nZTGvsphXWcyrPWvqJUmSJEkqVCedemvqJUmSJEnqo913h4svtqZekiRJkqShMtSdeus6ymJeZTGv\ncphVWcyrLOZVFvMqi3mVxbzas6ZekiRJkqRCWVM/BtbUS5IkSZIG0Z57wgUXWFMvSZIkSdJQGepO\nvXUdZTGvsphXOcyqLOZVFvMqi3mVxbzKYl7tWVMvSZIkSdIQsqZekiRJkqQ+et3r4PzzramXJEmS\nJKk4Dr+fIOs6ymJeZTGvcphVWcyrLOZVFvMqi3mVxbzas1MvSZIkSdIQsqZekiRJkqQ++v734Rvf\ngA9/uPU6rWrq7dRLkiRJkjTgPFFeE9Z1lMW8ymJe5TCrsphXWcyrLOZVFvMqi3n1zlB36iVJkiRJ\nKpnD7yVJkiRJGnAOv5ckSZIkaYoZ6k69dR1lMa+ymFc5zKos5lUW8yqLeZXFvMpiXr0z1J16SZIk\nSZJKZk29JEmSJEkDzpp6SZIkSZKmmKHu1FvXURbzKot5lcOsymJeZTGvsphXWcyrLObVO0PdqZck\nSZIkqWTW1EuSJEmSNOCsqZckSZIkaYrpS6c+ImZGxDkRcV9E3BoRB7VZd5OIuDQi7o+IJRHx5lHL\nj4iI39fLr42ILcfaDus6ymJeZTGvcphVWcyrLOZVFvMqi3mVxbx6p19H6o8B3gicAPwAODEidm6x\n7mnAzsAhwE3AWRGxBUBEvBv4EPBjYH/gUmDtsTZi8eLFE2y++sG8ymJe5TCrsphXWcyrLOZVFvMq\ni3n1zvQ+7Xcf4NrMPCoi5gB7AguAq0auFBEbAbsDl2TmKRHxM+BKYB5wJPA+YAnVFwTTMvPB8TRi\n2bJlnT4PTSLzKot5lcOsymJeZTGvsphXWcyrLObVO5N+pL7uqK8PLK1nNe7nNFl99qh1GvebR8Tj\ngKdSHZn/A3BfRFwSETO73mhJkiRJkgZQzzr1dZ37itE3qqP0q6w6ns2O+Hmd+v7JwEHAx4BXAR8Y\n68aWLFkyjl2r38yrLOZVDrMqi3mVxbzKYl5lMa+ymFfv9OySdvWw+hlNFt0G3AwszcxtImJz4Hrg\n7MzcJyKmA9Mz84H6qP4dwKWZ+eqIeDlwBXB0Zh4ZEfcA92bmphHxV8D/Aedk5ptH7zQivJ6dJEmS\nJKlYzS5p17Oa+sz8batlEXEWcGBEHAFsX88+s77/LDAvInbIzJ9GxMXArhHxXqra+0eAz494zLsi\n4mBg63ret1u0ZzwjAiRJkiRJGnj9OlHe4VTD5g8G7gEOzszGSfKyvjXsD3wOOBa4HViQmTfUyw4D\nngT8P+Bu4Hjg0z1vvSRJkiRJA6Bnw+8lSZIkSVJv9es69V0XEf8REbfXJ+S7aMT8fSPihohYHhGX\nRcQm9fz5zU7kFxGbtXucumO8edXLFkbEkoh4ICJujIh3jVj2o4i4JyLui4j/joidJvs5TWXdzGtN\n7z11ZoJZtVvWdHvqjh7ktWTU++p/J/s5TWU9yMvPGj0UEVtGxFURcWf9GeGbUZ3LiYh4SUT8vP4f\n9ZOI2H7E49ot8/NGj3Q7Lz9v9E4HWbX8TOHnjc5MmU491ZD9c0b8TETsQFWj/zuqs+LPBT5Zr7MI\n2Lu+vQX4C9VJ/Jau4XHqjnHlFRFbAicBDwPvoToJ439ExKb1Nr4HHAgcDWxXb0fd0828FtHivTcp\nz2TqG29Wa/p7t9r21FW9yOtqVr7HxnxFGI1J1/Lys8akaHxJcgRwBvAK4LMRsQ5wHvA4YCFVSehX\norJui2WNz8x+3uidbue1CD9v9Mp4s2q8f9p9pvDzRicyc8rcgKcDK4AL6+n31tP/WE//gOpEexuO\netwb6vWOGc/jvE1eXsBf1cuurn/+H2A58PgR23sC8ELgz8Av+/38ptqt23nVj1nlvedt0rPaaCx/\n70Zvz9vg5gUsofqAtV6/n9dUvXUprzG997x1nNWMUdN/pOrU7Vn/7g+q53+4nn458NpWy0Zsx88b\nBeVVz/fzRp+zGrFuy88Uft6Y+G0qHamH1a95/4f6fqeIeCawZT09e9R676D6R/rpcT5OnRlzXpn5\na+AQ4G+oLl24HfD2zPwjQETMqh//Q6pvYv+5x20fRl3La4TR7z11x3j+Fo7l751XD+mtbuaVwDzg\nnnoY4749afFw61Zeftboscx8qPFzPTJiQ6qrJM2pZy8ddb95m2Vz6u34eaNHepHXCH7e6KIJZDUy\nj3afKfy8MUFTrVM/2pephkntB/ySaghwAA80VoiILYBdgMsy8+axPk490fL3HhFPpBru9r/AHsDP\ngI+PGH5/L/C3wLuBdamGxam3Osmr1XtPvdEqq/vbLPPvXf90ktdngL2AN9fzPhURsyev6UNpInn5\n3ptE9ZcmFwK/pfrfNLrj0JhuNuR39Lp+3uixLufl540eGkdW6rEp3anPzL8ALwW2BZ4N/JjqH+mN\nI1Z7R31/2jgfpy5bw+99LlX9zlcz8yLgfGAm8KL6sY9k5pWZ+THgv4GdI2KjSX8SQ6STvGqrvffU\nG+2y8u/d4Okkr8z8SGZ+NTPPAb4CrMXKI8DqgYnm5XtvckTE1lS11Q9QDQG+nZW/46fV940vnG9c\nwzI/b/RYt/Oq+XmjByaQlXqoX9ep77qI2I3qnyLAZhHxVuC7wDupjha+gOpbun/PzAfrx6wNzAdu\nysyvj9jWNODkVo9T58abV0Q0/hi8JSJupzoKlcB1EbEr8A/A96n+iPw1cFtm/mnSntAU18286u01\nfe+pcxPIqu3fuxbbuzozr5+s5zSVdTOviHgucCxwKdX/93lU57K4ZhKf0pTW5bz8rNFjEfE04Cqq\ncxicDLw4Il4EfI1qCP3+EfFn4K1URxoXAWu3Wubnjd7qdl71Nv280QMTzKrtZwo/b3So30X93bpR\nvbBWUNXLNO4XUP2zvB+4AziVESd2oDob5iPAoaO2Fe0e561veb2H6pu++4Hrgf3r+TtQfWhdDtwF\nXAk8v9/PcSrduplXvazpe8/b5Ge1pr93LbY3r9/Pc6rcupkX8BTgknr+fVRHfv+2389xKt26nJef\nNXqf19xRea0AHqmX7QT8HHgQ+AnwvBGPa7rMzxtl5VUv8/PGYGXV8jNFu2Xe1nyL+pcoSZIkSZIK\nM6Vr6iVJkiRJmsrs1EuSJEmSVCg79ZIkSZIkFcpOvSRJkiRJhbJTL0mSJElSoezUS5IkSZJUKDv1\nkiRpzCLiQxFxUJvle0TEsyazTZIkDTM79ZIkaTxyDctfC2w9GQ2RJEkQmWv63yxJkoZZRBwGzAP+\nAPwO+AlwN/B2YG3geuAtwPbARfWyu4HXUR1A+BjwRGA58LbM/PUkPwVJkqYsO/WSJKmliHg+cAbw\nQmAG8FPgNODMzPxTvc7RwO2Z+bGIOAO4KDO/Wi+7EnhHZl4fETsCH8nMXfrxXCRJmoqm97sBkiRp\noO0EfDUzHwAeiIgLgQCeExHHABsA6wGXjXhMAETEesCLgXMjorFs7clquCRJw8BOvSRJaiepO+mj\nnAHskZnXRMQ+wNxRj4Fq6P2yzNy+t02UJGl4eaI8SZLUzreBPSNi3YiYCexez58J3BYRM4B/YmVH\n/l5gfYDMvAf4bUS8ASAqz53U1kuSNMVZUy9JktqKiEOBfahOlHcTVV39cuADwB3Aj4D1MnPfiPhr\n4DPAA8AbqDr7pwEbU9Xkn5OZx0z6k5AkaYqyUy9JkiRJUqEcfi9JkiRJUqHs1EuSJEmSVCg79ZIk\nSZIkFcpOvSRJkiRJhbJTL0mSJElSoezUS5IkSZJUKDv1kiRJkiQVyk69JEmSJEmF+v/DoqEKsYVq\nYwAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x1108f6350>"
]
}
],
"prompt_number": 73
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" merged[merged['diff'] != 0]\n",
" "
],
"language": "python",
"outputs": [
{
"html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>correct-extent</th>\n",
" <th>missing_x</th>\n",
" <th>correct-source</th>\n",
" <th>current-extent</th>\n",
" <th>missing_y</th>\n",
" <th>current-source</th>\n",
" <th>diff</th>\n",
" </tr>\n",
" <tr>\n",
" <th>date</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1987-08-03</th>\n",
" <td> 8.537</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 8.538</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1988-11-17</th>\n",
" <td> 11.507</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 11.508</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1988-12-31</th>\n",
" <td> 14.196</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.197</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1989-01-07</th>\n",
" <td> 14.367</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.368</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1989-09-07</th>\n",
" <td> 6.959</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.961</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.002</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1989-09-29</th>\n",
" <td> 7.206</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 7.207</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1989-10-18</th>\n",
" <td> 9.044</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 9.048</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.004</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1990-01-23</th>\n",
" <td> 14.971</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.972</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1990-04-12</th>\n",
" <td> 14.704</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.705</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1990-07-30</th>\n",
" <td> 7.860</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 7.861</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1990-09-21</th>\n",
" <td> 6.018</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.019</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1992-02-04</th>\n",
" <td> 15.349</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 15.351</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.002</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1992-11-28</th>\n",
" <td> 12.394</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 12.395</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-06-30</th>\n",
" <td> 11.065</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 11.066</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-09-08</th>\n",
" <td> 6.166</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.168</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.002</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-11-22</th>\n",
" <td> 11.587</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 11.588</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-12-04</th>\n",
" <td> 12.461</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 12.462</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-12-23</th>\n",
" <td> 13.584</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 13.588</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.004</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1993-12-24</th>\n",
" <td> 13.626</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 13.627</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1994-04-21</th>\n",
" <td> 14.683</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.685</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.002</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1994-06-28</th>\n",
" <td> 11.351</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 11.354</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.003</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1994-10-07</th>\n",
" <td> 8.117</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 8.118</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1994-10-15</th>\n",
" <td> 8.597</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 8.598</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1994-10-16</th>\n",
" <td> 8.724</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 8.725</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1995-01-01</th>\n",
" <td> 14.110</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 14.111</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1995-09-02</th>\n",
" <td> 6.112</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.113</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1995-09-08</th>\n",
" <td> 6.016</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.017</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1995-09-21</th>\n",
" <td> 6.035</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.036</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1995-10-02</th>\n",
" <td> 6.146</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.147</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2006-08-14</th>\n",
" <td> 6.513</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.514</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-04</th>\n",
" <td> 5.184</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.194</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.010</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-05</th>\n",
" <td> 5.157</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.171</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.014</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-06</th>\n",
" <td> 5.159</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.180</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.021</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-07</th>\n",
" <td> 5.158</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.186</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.028</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-08</th>\n",
" <td> 5.099</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.132</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.033</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-09</th>\n",
" <td> 5.069</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.107</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.038</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-10</th>\n",
" <td> 5.086</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.124</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.038</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-11</th>\n",
" <td> 5.066</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.105</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.039</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-12</th>\n",
" <td> 5.062</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.097</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.035</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-13</th>\n",
" <td> 5.049</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.077</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.028</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-14</th>\n",
" <td> 5.073</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.102</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.029</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-15</th>\n",
" <td> 5.100</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.133</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.033</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-16</th>\n",
" <td> 5.090</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.133</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.043</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-17</th>\n",
" <td> 5.091</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.137</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.046</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-18</th>\n",
" <td> 5.080</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.121</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.041</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-19</th>\n",
" <td> 5.112</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.157</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.045</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-20</th>\n",
" <td> 5.115</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.167</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.052</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-21</th>\n",
" <td> 5.162</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.215</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.053</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-22</th>\n",
" <td> 5.199</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.255</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.056</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-23</th>\n",
" <td> 5.250</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.297</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.047</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-24</th>\n",
" <td> 5.295</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.331</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.036</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-25</th>\n",
" <td> 5.341</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.365</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.024</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-26</th>\n",
" <td> 5.390</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.406</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.016</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-27</th>\n",
" <td> 5.458</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.471</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.013</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-28</th>\n",
" <td> 5.542</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.551</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.009</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-09-29</th>\n",
" <td> 5.571</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.572</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-10-02</th>\n",
" <td> 5.864</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.868</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.004</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-10-03</th>\n",
" <td> 5.917</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 5.922</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.005</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-10-04</th>\n",
" <td> 6.030</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.033</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.003</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2013-10-05</th>\n",
" <td> 6.159</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td> 6.160</td>\n",
" <td> 0</td>\n",
" <td> ftp://sidads.colorado.edu/pub/DATASETS/nsidc0...</td>\n",
" <td>-0.001</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>132 rows \u00d7 7 columns</p>\n",
"</div>"
],
"output_type": "pyout",
"prompt_number": 75,
"text": [
" correct-extent missing_x \\\n",
"date \n",
"1987-08-03 8.537 0 \n",
"1988-11-17 11.507 0 \n",
"1988-12-31 14.196 0 \n",
"1989-01-07 14.367 0 \n",
"1989-09-07 6.959 0 \n",
"1989-09-29 7.206 0 \n",
"1989-10-18 9.044 0 \n",
"1990-01-23 14.971 0 \n",
"1990-04-12 14.704 0 \n",
"1990-07-30 7.860 0 \n",
"1990-09-21 6.018 0 \n",
"1992-02-04 15.349 0 \n",
"1992-11-28 12.394 0 \n",
"1993-06-30 11.065 0 \n",
"1993-09-08 6.166 0 \n",
"1993-11-22 11.587 0 \n",
"1993-12-04 12.461 0 \n",
"1993-12-23 13.584 0 \n",
"1993-12-24 13.626 0 \n",
"1994-04-21 14.683 0 \n",
"1994-06-28 11.351 0 \n",
"1994-10-07 8.117 0 \n",
"1994-10-15 8.597 0 \n",
"1994-10-16 8.724 0 \n",
"1995-01-01 14.110 0 \n",
"1995-09-02 6.112 0 \n",
"1995-09-08 6.016 0 \n",
"1995-09-21 6.035 0 \n",
"1995-10-02 6.146 0 \n",
"2006-08-14 6.513 0 \n",
"... ... ... \n",
"2013-09-04 5.184 0 \n",
"2013-09-05 5.157 0 \n",
"2013-09-06 5.159 0 \n",
"2013-09-07 5.158 0 \n",
"2013-09-08 5.099 0 \n",
"2013-09-09 5.069 0 \n",
"2013-09-10 5.086 0 \n",
"2013-09-11 5.066 0 \n",
"2013-09-12 5.062 0 \n",
"2013-09-13 5.049 0 \n",
"2013-09-14 5.073 0 \n",
"2013-09-15 5.100 0 \n",
"2013-09-16 5.090 0 \n",
"2013-09-17 5.091 0 \n",
"2013-09-18 5.080 0 \n",
"2013-09-19 5.112 0 \n",
"2013-09-20 5.115 0 \n",
"2013-09-21 5.162 0 \n",
"2013-09-22 5.199 0 \n",
"2013-09-23 5.250 0 \n",
"2013-09-24 5.295 0 \n",
"2013-09-25 5.341 0 \n",
"2013-09-26 5.390 0 \n",
"2013-09-27 5.458 0 \n",
"2013-09-28 5.542 0 \n",
"2013-09-29 5.571 0 \n",
"2013-10-02 5.864 0 \n",
"2013-10-03 5.917 0 \n",
"2013-10-04 6.030 0 \n",
"2013-10-05 6.159 0 \n",
"\n",
" correct-source current-extent \\\n",
"date \n",
"1987-08-03 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 8.538 \n",
"1988-11-17 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 11.508 \n",
"1988-12-31 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.197 \n",
"1989-01-07 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.368 \n",
"1989-09-07 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.961 \n",
"1989-09-29 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 7.207 \n",
"1989-10-18 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 9.048 \n",
"1990-01-23 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.972 \n",
"1990-04-12 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.705 \n",
"1990-07-30 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 7.861 \n",
"1990-09-21 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.019 \n",
"1992-02-04 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 15.351 \n",
"1992-11-28 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 12.395 \n",
"1993-06-30 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 11.066 \n",
"1993-09-08 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.168 \n",
"1993-11-22 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 11.588 \n",
"1993-12-04 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 12.462 \n",
"1993-12-23 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 13.588 \n",
"1993-12-24 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 13.627 \n",
"1994-04-21 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.685 \n",
"1994-06-28 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 11.354 \n",
"1994-10-07 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 8.118 \n",
"1994-10-15 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 8.598 \n",
"1994-10-16 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 8.725 \n",
"1995-01-01 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 14.111 \n",
"1995-09-02 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.113 \n",
"1995-09-08 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.017 \n",
"1995-09-21 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.036 \n",
"1995-10-02 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.147 \n",
"2006-08-14 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.514 \n",
"... ... ... \n",
"2013-09-04 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.194 \n",
"2013-09-05 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.171 \n",
"2013-09-06 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.180 \n",
"2013-09-07 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.186 \n",
"2013-09-08 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.132 \n",
"2013-09-09 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.107 \n",
"2013-09-10 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.124 \n",
"2013-09-11 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.105 \n",
"2013-09-12 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.097 \n",
"2013-09-13 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.077 \n",
"2013-09-14 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.102 \n",
"2013-09-15 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.133 \n",
"2013-09-16 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.133 \n",
"2013-09-17 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.137 \n",
"2013-09-18 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.121 \n",
"2013-09-19 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.157 \n",
"2013-09-20 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.167 \n",
"2013-09-21 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.215 \n",
"2013-09-22 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.255 \n",
"2013-09-23 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.297 \n",
"2013-09-24 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.331 \n",
"2013-09-25 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.365 \n",
"2013-09-26 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.406 \n",
"2013-09-27 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.471 \n",
"2013-09-28 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.551 \n",
"2013-09-29 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.572 \n",
"2013-10-02 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.868 \n",
"2013-10-03 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 5.922 \n",
"2013-10-04 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.033 \n",
"2013-10-05 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... 6.160 \n",
"\n",
" missing_y current-source \\\n",
"date \n",
"1987-08-03 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1988-11-17 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1988-12-31 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1989-01-07 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1989-09-07 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1989-09-29 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1989-10-18 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1990-01-23 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1990-04-12 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1990-07-30 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1990-09-21 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1992-02-04 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1992-11-28 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-06-30 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-09-08 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-11-22 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-12-04 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-12-23 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1993-12-24 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1994-04-21 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1994-06-28 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1994-10-07 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1994-10-15 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1994-10-16 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1995-01-01 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1995-09-02 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1995-09-08 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1995-09-21 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"1995-10-02 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2006-08-14 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"... ... ... \n",
"2013-09-04 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-05 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-06 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-07 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-08 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-09 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-10 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-11 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-12 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-13 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-14 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-15 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-16 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-17 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-18 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-19 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-20 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-21 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-22 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-23 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-24 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-25 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-26 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-27 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-28 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-09-29 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-10-02 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-10-03 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-10-04 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2013-10-05 0 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"\n",
" diff \n",
"date \n",
"1987-08-03 -0.001 \n",
"1988-11-17 -0.001 \n",
"1988-12-31 -0.001 \n",
"1989-01-07 -0.001 \n",
"1989-09-07 -0.002 \n",
"1989-09-29 -0.001 \n",
"1989-10-18 -0.004 \n",
"1990-01-23 -0.001 \n",
"1990-04-12 -0.001 \n",
"1990-07-30 -0.001 \n",
"1990-09-21 -0.001 \n",
"1992-02-04 -0.002 \n",
"1992-11-28 -0.001 \n",
"1993-06-30 -0.001 \n",
"1993-09-08 -0.002 \n",
"1993-11-22 -0.001 \n",
"1993-12-04 -0.001 \n",
"1993-12-23 -0.004 \n",
"1993-12-24 -0.001 \n",
"1994-04-21 -0.002 \n",
"1994-06-28 -0.003 \n",
"1994-10-07 -0.001 \n",
"1994-10-15 -0.001 \n",
"1994-10-16 -0.001 \n",
"1995-01-01 -0.001 \n",
"1995-09-02 -0.001 \n",
"1995-09-08 -0.001 \n",
"1995-09-21 -0.001 \n",
"1995-10-02 -0.001 \n",
"2006-08-14 -0.001 \n",
"... ... \n",
"2013-09-04 -0.010 \n",
"2013-09-05 -0.014 \n",
"2013-09-06 -0.021 \n",
"2013-09-07 -0.028 \n",
"2013-09-08 -0.033 \n",
"2013-09-09 -0.038 \n",
"2013-09-10 -0.038 \n",
"2013-09-11 -0.039 \n",
"2013-09-12 -0.035 \n",
"2013-09-13 -0.028 \n",
"2013-09-14 -0.029 \n",
"2013-09-15 -0.033 \n",
"2013-09-16 -0.043 \n",
"2013-09-17 -0.046 \n",
"2013-09-18 -0.041 \n",
"2013-09-19 -0.045 \n",
"2013-09-20 -0.052 \n",
"2013-09-21 -0.053 \n",
"2013-09-22 -0.056 \n",
"2013-09-23 -0.047 \n",
"2013-09-24 -0.036 \n",
"2013-09-25 -0.024 \n",
"2013-09-26 -0.016 \n",
"2013-09-27 -0.013 \n",
"2013-09-28 -0.009 \n",
"2013-09-29 -0.001 \n",
"2013-10-02 -0.004 \n",
"2013-10-03 -0.005 \n",
"2013-10-04 -0.003 \n",
"2013-10-05 -0.001 \n",
"\n",
"[132 rows x 7 columns]"
]
}
],
"prompt_number": 75
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"plt.plot((current_data.extent - correct_data.extent))\n",
"plt.title('extent differences with and without N07 Mask (in Mkm^2)')\n",
"plt.show\n",
"print(min(current_data.extent - correct_data.extent), max(current_data.extent - correct_data.extent))"
],
"language": "python",
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"(0.0, 0.05600000000000005)\n"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAfQAAAHpCAYAAAB0ns7PAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYZFV9//H3l2GVYXVhGVA0YALGBEwCY9TQRqOTEeH3\n0yRIUASM8FNxjYa4JI5L4hYVeFREREUgIhEkowHBrdW4AKOALIMyIMuAM6Iw7MsM8/39cW/N1FR3\nddV0V3d1nXq/nqefrrufe2r53HPurVuRmUiSpMG2Sb8LIEmSps5AlySpAAa6JEkFMNAlSSqAgS5J\nUgEMdEmSCmCgD6mIWBsRT6kfnxwR72qa9pqIWBkR90TEDhHxrIi4PiLujYiD+1fq/mmto3GmL4qI\nM2ayTE3b/kJEvG+K63h7RJw6wfQjI+IHU9mGKhGxR/3+m9Tnb0Q8PiKWRsQW9fAFEfGK3pZyzDan\n/BrrpYj4o4j4Yb/LMdsY6H3Q6zdHRNwUEX852eUz8zWZ+f56XZsBHwWel5nbZuZdwHuBkzJzm8xc\n3JtSD5aWOhqJiFtbZ+lDsZq3PaXtZ+YHMvPVMPXA2VidXr91fa+NiE+2jP/fiHhl0/DfR8TNEXFf\nRHw1InZomnZNfUDa+FsdEeO+lpu2d17L+D+ux3938nvbE/8MfD4zHwbIzIWZOamDybruH46Ix7aM\nv7ze1yfWo6b8Ghtn27tHxK1142HPlmmbR8RpdfnuqcuzoDE9M38OrIqIg3pZpkFnoJchgejRunYG\ntgSWNo17InDtZFYWEXN6UagB0Kv6n03bn6l96ub1ez/w8oh4UstyCRARTwM+DRwO7AQ8AHxq3YyZ\nT6sPSLfJzG2AW4FzJtjeHcD8iNixadwrgV/Sx4O3ulV+BHBmj1aZwI3AYU3beDqwFWP3s2evh/oA\n4mLgdOBjwEURsXPTLJsCtwB/kZnbAu8Czml5/s8Cju1VmUpgoE9BROwaEedGxG8i4saIeH09fsf6\nyPOgenhuRCyLiFdExKuBvwf+qW4p/PdE66qnLYqIcyLi9Ppo9eqI+JN62hlUgfu1en1vbVPWt0XE\n7RGxPCKObpn2hYh4X0TsBVxXj14VEd+OiGXAU+r13xMRm0XEdvXRc2N972u05uqu2R9GxMci4rfA\nu+uj7f+oW08rouq+3rKef6Rex1vqI/XbI+LIprJtFREfrY/UV0XED5qWnR8RP4qIuyLiiog4sGm5\nIyPihrrMN0bE349TJ1tGxIOND+2IeGfdcptbD78vIj7eUkePAS4Edq3r+56I2IXqw2/z8Z6jNs/H\niRFxS0TcHRFLIuLZ3Tzf9fT9IuJn9bSzqQ7A2m3n5oh4Rv348KhaXXvXw6+KiK82bbPRyvt+02vg\nnoiYz/rg/EhE3FnX6YKm7ewaEYsj4ndRnZ75h6ZpG/RIRVMPR7evX2AV8AXg3W2mHw4szsz/zcz7\ngX8BXhIRW49TJwcCjwPObVdvwCPA+cDL6mXmAH9HFSLrgq3D87h/Pe7u+nX/0fE2FBEvjYhfRcQ+\nE5Sn4QBgVWbe3rT8aES8qn58ZFQ9F+M+T22cSXWQ0PBK4IvN+9lS3m0i4rsRcUI9/IWI+FRUXf/3\n1tvfKSJOqMuwNCL2bVp+a+B/gDMz812Z+SHgROCCiNgWIDMfyMz3ZOYt9fD/AL8CntFUlO8Bz4uq\nV1EY6JMWVYB9Dbgc2BV4HvCmiHhBZt4JHA2cGhGPBz4O/Cwzz8jMU6k+FD5UtxYOmWhdTZt8MfAl\nYDtgMfAJgMx8BdWR7EH1+v5jnLIuAP4ReD7w1Pp/s6xWldcDjQ+V7TLzeZm5Z9P6t83M1VQfrI8A\nvwfsB7wA+Iem9e0P3AA8Afh34EPAnsAf1//nAf/aNP9OwLb1vr8K+GREbFdP+496G88EdgTeBqyN\niHnA14H3ZuYOwFuBcyPisfUHxonAgvro/pnAFa31kpkPAZcCI/WoA4GbgGc3DY+21NEDwALg9rq+\nt83MX1N9+B3MOM9RG5fW9bED8J/Af0XE5k3Tx32+63nOp2rZ7AD8F/BS2rcaR1v274b6f+v+NXtO\n/X+7ev9+Uu/fAVQHfI8FPgyc1rTM2VSvk12AvwH+PSKeW09r213bzeu3yb8DL42Ip44zbR/gyqb1\n3gg8TPV6b/VK4CuZ+eAE2wI4g/VB90LgauD2lnkmeh5PBD6emdtRHRS39ghERBwFfJDqFFc3vWBP\nB37RMq61fven/fM0np8A20bEH9QHLocyfg9ARtWy/jbwg8x8U9O0vwXeSXWg9DDwY2BJXYavULXC\nG/4MODEz/23dijNPAt4PzB+vgBGxE9VzeU3TMrcBq4Hf77B/Q8NAn7w/Ax6Xme/PzDWZ+Svgs9RH\n9Jn5TaoP2+9QBUBr11B0u67aDzLzG1ndfP9Mqg+Rbv0d8LnMvLYOpPFaOdHyf1z1G+uvgTdn5oOZ\neQdwQktZb8/MT2bmWqo396uBt2Tmqsy8D/hAy/yrqYL50cy8ELgP+P36QOco4I2Z+evMXJuZP8nM\nR4CXAxdk5jcAMvNbVB8gL6L6cFsLPD0itsrMlRN8WH4POLD+IHs6cFI9vCXwp6xvrXZTR10/R5l5\nVmbeVe/Tx4At2PCDqd265gObZuaJdX2dC1zWbjuN/asfP5uq7hvDf1FPb9Vu/27OzNPqMn0R2CUi\nnhARuwN/DhyfmY9k5pVUr9/mVt+Uu2szcyVVt/p7x5k8F7i7Zdw9wDbNI+oelpdSHZR22t6PgR3r\nA4gjqA6iWueZ6Hl8BNgrIh5XtzgvaVn8zVQHogfWByDd2B64t8M84z5PHZZpHLz8FdXptdvGmWce\n1QHglzOz+YA8gfMy8/L6vP5XgQcz88y6DOdQHZRXM2eOZuaXWleemedl5sWt4+sW+FnAFzLzly2T\n76WqE2GgT8WTqLpd72r8AW+napU2nAo8jeqFeNcU17Wy6fEDwJbR/UVLu1CdM2y4pcvl2pV1M+DX\nTWX9NPD4pnmat/V44DHAT5vmv5DqSL7hd3X4NzxA9QH9OKru5BvalONvW+rsWcDO9UHLocD/A26P\niK9HRLuj+O9RtWCfAVwFfIsq8A4AlnV43lp1/RxFxFsj4tqoTiPcRdUSb66TduvalbEftjfTPjC/\nDzwnqvOTc6gOMp8V1bnI7TJzTM/FBFY0HtR1DNXztCtwZ93V3XALVQD02oeBF0bEH7WMv4+qDptt\nx9jwewnV6+37dOcM4PVUr5Gv0lLPHZ7HV1G1KpdGxKUR8aKWdf8j8Mnm7vMu3EnLQco42j1P7STV\nfh5O++72oDpY3hI4ZZx1/Kbp8UMtww922H5b9Wv+jHqdx40zyzZUp2NEdeGBJucW4FeZOV6XXuOc\n22eo3hyvi4gvZGYjmFq7Hydc1zjzb+z0X1Odp2x4YrsZu3ArVav7sS0h3K48v6V6Q+9Td01vjN9S\nvZH3BH7eMu0W4IzMPGbcAlRH+hdHdRHRv1EdXP3FOLP+mKpF9X+B0cxcGtWVvQsZ2x2dLf/Hm9ZR\nRDyH6tTBX2bmNfW4O+muFftrxgblk4Bl482cmcsi4gGqUPpeZt4bESuAY4Dmr6Flm8fduJ2qJTu3\n7oGB6jW2vH58P9VBXUPzxU8btb3M/F197vb9LZOuoalHJCJ+D9ic6iK2Zo3A6taZwPXA6Zn5UMT6\np6jT85iZy6iulyEiXgp8JTa8yO4FVBeDrcjMDa6on8DPqVr2PZWZt0TEjVS9b0ePNwvVe2gHqnPd\nC5oOFqZFVJV9GlWjYGFmPtoyfR7Vc9x6CmJo2UKfvEuBeyPin6K6cGtORPxhRPxpPf0dwKNUXcYf\nAb7Y1FpbSXVOrdt1dfqgX0l1Prudc4AjI2Lvusuxtcu96+7QOpQvBj5WXxyzSUT8XkSMF5bUoX8q\ncEJ9PQERMa/l+oB221oLfK7e1i51vTyzPkd5JvDiiHhBPX7LqC62mld3Ax9Sn0tfTRUoj7bZxgPA\nT4HXsb77+UdUrfvm7uhgfT2tBB7buICnaXq3tgHWAL+N6oLBf6W6hqAbPwbWRMQborpA8SVUp2wm\n8j2q1k1jf0ZbhmHD8t9BdcpiotfUOpl5K1WdfSAitqhbz0ez/jzsFcDCqO5psDPwppZVdHr9tvoY\n1XURezeNO4vq9fDs+nl/H3Buc69BROxG1dIe03U+wb79iupA8J3jTJ7weYyIlzde81SnAxqnghqu\noTod98mIeHGXRboM2D4idu12HzbCq6gOTsa7tqBxkHIcVYB+rT4ttW7aNDgZ+APg4Lorv9WBwLez\nuq5HGOiTVofNQcC+VF/7uIOqRb5tVFckvxk4oj6H9CGqN/Px9eKnAfvUXcXnTbSuxuYY24ppHv4A\n8K56fW8Zp6zfoDrP/R2qFsu3Gdsi25gW2hFUR8bXUnUB/hfrW13jlfV4qhbkTyLibuCbbHix0kTb\neytVV/hlwO+o9nWTzFwOHEJ14PQbqhb7P1J9uGxCVf+31cs8B3jNBNv4HlVv1aVNw3PZ8Pz5uv3K\nzOuoLli7MaqreHeh83PU7Bv13y+pLsJ7kA1Pg7RdV339wEuAI+t9+zsmvlp7vP3ptH8PUPVq/LDe\nvwO62L/DgD2oWuvnAf+amd+pp51BdcHaTfV+n81GvH5bt5WZ91J1ve/QNO5aqoOws6gOELYCXtuy\njlcAP6pDupPm7f0oM1c0jW9M6/Q8vhC4OiLupbow9mVNwdSo659TvfdPjYgXdixU9fx/geoaknbl\n7vZ12LruGzPzZ22Wa17vMVS9L+fXPWDjfX5MqgwN9SmhY6h6XVbE+vsHHNY02+FUp/tUiypvJpih\nukL6BKrzb5/N6isGrfOcRNVV8wBwZGZeXo/fnurimKdRPaFHZ3XFrCRpEiLicVSnS/Zt03ItXt0L\ndHJmPqvfZZlNJgz0+jzwL6i+5nQbVSvpsMxc2jTPQuC4zFxYH8mfmJnz62mnU523+1xEbApsnZmt\nV6JKkqQp6tTlvj/Vlb431ecpzqbq5mx2MPU5qfprGdtHdVOB7YDnZObn6mlrDHNJkqZHp0Cfx4Zf\nQVrO2Ctsx5tnN+DJwB0R8fmo7mp1an1BliRJ6rFOX1vr9kKG1qscs173M6i64y+rv2ryz2x4hzAi\nop8/aiFJUl9kZk+/IdCphX4bsHvT8O6s/25pu3l2q8ctB5ZnZuMuVl9hw/vwrpOZ/nXx9+53v7vv\nZRiEP+vJurKerKfZ/jcdOgX6EqpbF+5Rf/f3UKr7SjdbTH2Lx6h+xGFVVrfaXAHcGuvvu/x8mu7D\nK0mSemfCLvfMXBMRxwEXUX1t7bSs7qR1bD39lMy8ICIWRvWrXPdT3Uil4fXAWfXBwA0t0yRJUo90\nvPVrVj+WcWHLuFNahse7xy5Z/UhDp7tYqUsjIyP9LsJAsJ66Z111x3rqjvXUXx1vLDPtBYjIfpdB\nkqSZFBHkDF8UJ0mSBoCBLklSAQx0SZIKYKBLklQAA12SpAIY6JIkFcBAlySpAAa6JEkFMNAlSSqA\ngS5JUgEMdEmSCmCgS5JUAANdkqQCGOiSJBXAQJckqQAGuiRJBTDQJUkqgIEuSVIBDHRJkgpgoEuS\nVAADXZKkAhjokiQVwECXJKkABrokaWisWAErV/a7FNMjMrO/BYjIfpdBkjQctt8e5s6F5cv7W46I\nIDOjp+vsd5ga6JKkmRIBm28ODz/c73L0PtDtcpckqQAGuiRJBTDQJUlDJXra0T17GOiSpKFS6mVb\nBrokSQUw0CVJQ8Uud0mSNGsZ6JIkFcBAlySpAAa6JEkFMNAlSUPFi+IkSdKsZaBLklQAA12SpAIY\n6JIkFcBAlySpAAa6JGmoeJW7JEkF8NfWJEnSrGWgS5KGil3ukiRp1jLQJUkqgIEuSVIBDHRJkgpg\noEuSVAADXZI0VLzKXZIkzVoGuiRJBTDQJUkqgIEuSVIBDHRJkgpgoEuShopXuUuSVAB/PlWSJM1a\nBrokaajY5S5JkmYtA12SpAIY6JIkFcBAlySpAB0DPSIWRMR1EXF9RBzfZp6T6ulXRsR+TeNvioif\nR8TlEXFpLwsuSZLW23SiiRExB/gE8HzgNuCyiFicmUub5lkI7JmZe0XEAcDJwPx6cgIjmXnntJRe\nkqSNNKxXue8PLMvMmzJzNXA2cEjLPAcDpwNk5iXA9hGxU9P0QqtOkjSIhvXGMvOAW5uGl9fjup0n\ngW9FxJKIePVUCipJktqbsMudKpC70a4V/uzMvD0iHg98MyKuy8wftM60aNGidY9HRkYYGRnpcrOS\nJG2cfnS5j46OMjo6Oq3biJyg7yEi5gOLMnNBPfx2YG1mfqhpnk8Do5l5dj18HXBgZq5sWde7gfsy\n86Mt43OiMkiS1CsRMHcu3Htvv8sRZGZPDy06dbkvAfaKiD0iYnPgUGBxyzyLgSPqAs4HVmXmyoh4\nTERsU4/fGngBcFUvCy9J0sYq9aK4CbvcM3NNRBwHXATMAU7LzKURcWw9/ZTMvCAiFkbEMuB+4Kh6\n8Z2B86KquU2BszLz4unaEUmSulFqp/CEXe4zUgC73CVJM2SYu9wlSSpKqV3uBrokSQUw0CVJKoCB\nLklSAQx0SZIKYKBLklQAA12SNFS8yl2SJM1aBrokSQUw0CVJKoCBLklSAQx0SZIKYKBLkoaKV7lL\nklSAUn/g00CXJKkABrokaajY5S5JkmYtA12SpAIY6JIkFcBAlySpAAa6JEkFMNAlSUPFq9wlSSqA\nN5aRJEmzVmSfD1UiIvtdBknScGh0t/c7diKCzOxp578tdEmSCmCgS5JUAANdkqQCGOiSJBXAQJck\nqQAGuiRJBTDQJUkqgIEuSRoqW23V7xJMDwNdkjRUdtut3yWYHga6JEkFMNAlSUPFX1uTJEmzloEu\nSRoq/f5hlulioEuSVAADXZKkAhjokqSh4kVxkiRp1jLQJUlDxYviJEnSrGWgS5JUAANdkjRUvChO\nkiTNWga6JGmoeFGcJEmatQx0SZIKYKBLklQAA12SpAIY6JIkFcBAlySpAAa6JEkFMNAlSSqAgS5J\nUgEMdEmSCmCgS5JUAANdkqQCGOiSJBXAQJckqQAGuiRJBTDQJUkqgIEuSVIBOgZ6RCyIiOsi4vqI\nOL7NPCfV06+MiP1aps2JiMsj4mu9KrQkSdrQhIEeEXOATwALgH2AwyJi75Z5FgJ7ZuZewDHAyS2r\neSNwLZC9KrQkSdpQpxb6/sCyzLwpM1cDZwOHtMxzMHA6QGZeAmwfETsBRMRuwELgs0D0suCSJGm9\nToE+D7i1aXh5Pa7beT4OvA1YO4UySpKkDjbtML3bbvLW1ndExEHAbzLz8ogYmWjhRYsWrXs8MjLC\nyMiEs0uSNFBGR0cZHR2d1m1EZvvMjoj5wKLMXFAPvx1Ym5kfaprn08BoZp5dD18HjABvAF4BrAG2\nBLYFzs3MI1q2kROVQZKkXomAvfaCX/6y3+UIMrOnp6I7dbkvAfaKiD0iYnPgUGBxyzyLgSPqAs4H\nVmXmisx8R2bunplPBl4GfKc1zCVJUm9M2OWemWsi4jjgImAOcFpmLo2IY+vpp2TmBRGxMCKWAfcD\nR7VbXS8LLkmS1puwy31GCmCXuyRphgxzl7skSRoABrokSQUw0CVJQ2XVqn6XYHp4Dl2SNDSiPmvd\n79jxHLokSRqXgS5JUgEMdEmSCmCgS5JUAANdkqQCGOiSJBXAQJckqQAGuiRJBTDQJUkqgIEuSVIB\nDHRJkgpgoEuSVAADXZKkAhjokiQVwECXJKkABrokSQUw0CVJKoCBLklSAQx0SZIKYKBLklQAA12S\npAIY6JIkFcBAlySpAAa6JEkFMNAlSSqAgS5JUgEMdEmSCmCgS5JUAANdkqQCGOiSJBXAQJckqQAG\nuiRJBTDQJUkqgIEuSVIBDHRJkgpgoEuSVAADXZKkAhjokiQVwECXJKkABrokSQUw0CVJKoCBLklS\nAQx0SZIKYKBLklQAA12SpAIY6JIkFcBAlySpAAa6JEkFMNAlSUMlot8lmB4GuiRpqBjokiQVwECX\nJEmzloEuSRoqmf0uwfQw0CVJQ8Uud0mSCmCgS5JUAANdkiTNWga6JEkFMNAlSSqAgS5JGiqeQ5ck\nqQBDG+gRsSAirouI6yPi+DbznFRPvzIi9qvHbRkRl0TEFRFxdUQs6nHZJUnaaEMZ6BExB/gEsADY\nBzgsIvZumWchsGdm7gUcA5wMkJkPAc/NzH2BfYEFEXFA73dBkiR1aqHvDyzLzJsyczVwNnBIyzwH\nA6cDZOYlwPYRsVM9/EA9z+bAZsDaXhVckiSt1ynQ5wG3Ng0vr8d1mmc3qFr4EXEFsBK4ODMvm1px\nJUmanFLv4d7QKdC73f3WMxIJkJmP1l3uuwEHRMTTNrJ8kiSpC5t2mH4bsHvT8O5ULfCJ5tmtHrdO\nZt4dEd+lOhd/TetGFi1atO7xyMgIIyMjHYolSdLgGB0dZXR0dFq3ETlBH0REbAr8AngecDtwKXBY\nZi5tmmchcFxmLoyI+cAJmTk/Ih4HrMnMVRGxFXAR8MHMvKBlGzlRGSRJ6oW1a2HOHNhyS3jwwf6W\nJSLIzJ5ebz9hCz0z10TEcVRhPAc4LTOXRsSx9fRTMvOCiFgYEcuA+4Gj6sV3AU6vr5TfBPhya5hL\nkqTemLCFPiMFsIUuSZoBpbfQvVOcJGkolN52NNAlSSqAgS5JGipDeetXSZJKYZe7JEkFsYUuSVIB\nSm2pG+iSJBXAQJckDYVGy9wud0mSNGsZ6JKkoWILXZKkAdbocn/KU/pbjulioEuShso++/S7BNPD\nQJckqQAGuiRpKJT6/fMGA12SpAIY6JKkoVJqS91AlyQNhVKDvMFAl6QB8eEPwyOP9LsUmq0MdEka\nEMcfD0uX9rsUmq0MdEmSCmCgS9IAKf088HRq1F2pdWigS5JUAANdkgZIqT8soqkz0CVpgJTaXTwT\nSq87A12SNHAy4YEH+l2K2cVAl6QBYpd75eSTYeutJ7dsqS11A12SBkipYbSxbrhh45cpve4MdEnS\nwLGnYiwDXZI0cAz0sQx0SdLAmUyge2MZSZJmGVvoYxnokqSBY6CPZaBLkgaOgT6WgS5JGgqlnjtv\nMNAlSUOl1GA30CVJA6fUUJ4KA12SNBRKPwgw0CVJKoCBLklSAQx0SdLAmUz3uXeKkyRJs56BLkkD\npNTWpabOQJckqQAGuiQNEG95Onml924Y6JI0QEoPpW5NpR5KrUMDXZKkAhjokjRA7HKfvFJb5g0G\nuiQNkNJDSZNnoEuSVAADXZI0cLxT3FgGuiQNEM+hqx0DXZKkAhjokqShUGpXe4OBLklSAQx0SdLA\n8U5xYxnokiQVwECXJA2FUlvmDQa6JEkFMNAlSUOl1Ja6gS5JGjhTuVNcqQx0SRog3ilO7RjokjRA\nSm9ldst6GMtAlyQNhdIPAgx0SRogdrlXplIPpQa7gS5JGjilhvJUGOiSNEAMsskrve66CvSIWBAR\n10XE9RFxfJt5TqqnXxkR+9Xjdo+I70bENRFxdUS8oZeFlyQNp9LDeTI6BnpEzAE+ASwA9gEOi4i9\nW+ZZCOyZmXsBxwAn15NWA2/OzKcB84HXtS4rSeqe59DVTjct9P2BZZl5U2auBs4GDmmZ52DgdIDM\nvATYPiJ2yswVmXlFPf4+YCmwa89KL0kaSl4UN1Y3gT4PuLVpeHk9rtM8uzXPEBF7APsBl2xsISVJ\nauad4sbatIt5uq2C1uOldctFxFzgK8Ab65b6BhYtWrTu8cjICCMjI11uUpKk2W90dJTR0dFp3UY3\ngX4bsHvT8O5ULfCJ5tmtHkdEbAacC5yZmeePt4HmQJckqTStjdX3vOc9Pd9GN13uS4C9ImKPiNgc\nOBRY3DLPYuAIgIiYD6zKzJUREcBpwLWZeUIPyy1J0kYZ+i73zFwTEccBFwFzgNMyc2lEHFtPPyUz\nL4iIhRGxDLgfOKpe/FnAy4GfR8Tl9bi3Z+Y3er4nkiR1odRg76bLncy8ELiwZdwpLcPHjbPc/+LN\nayRJPVZqKE+FYStJGgqlHwQY6JIkFcBAl6QB4p3i1I6BLkkaOFPpPi+1691Al6QBUmoYzYTS685A\nlySpAAa6JA0Qz6GrHQNdkjQU7HKXJGmW8aK4sQx0SZIKYKBLkoZCqS3zBgNdkqQCGOiSpKFSakvd\nQJckDZzJhHKpQd5goEuSVAADXZKkAhjokjRAvFOc2jHQJWmAlH4eeDo16q7UOjTQJUkDp9RQngoD\nXZIGiF3uasdAlyQNhdJb9Qa6JEkFMNAlSUOl1Ja6gS5JGjjeKW4sA12SpAIY6JI0AEpvXW4s62Ms\nA12SBohBNnml152BLkkaOFP5Pn6pwW6gS5IGTqmhPBUGuiQNAANMnRjokqSB49fWxjLQJWmAlB5K\n02n16n6XYHoZ6JI0AAzyDU3moriTTqr+l1qXBrokaeBMJpTvuqv35ZhNDHRJkgpgoEvSACm1u3hj\nWQ9jGeiSNAAMsKkrvQ4NdEnSwPFOcWMZ6JKkgTOZUJ7KQcAgMNAlaQCU2qqcSaXXoYEuSQOk9FDq\nlvUwloEuSQPAAFMnBrokaaiUenBkoEuShkKpQd5goEvSAGiEUemhpMkz0CVJQ8GvrUmS+s6W+dSV\n3sthoEuSBs5kQrnUIG8w0CVpAJQeRjPBLndJ0qxhsE9e6XVnoEuSVAADXZIGQOmty5lUal0a6JKk\ngVNqKE+FgS5JA6D0r1xp6gx0SdJQKP1gyECXpAFQehjNBL+2JklSAUo/bWGgS5IGTqmhPBUGuiQN\ngNJblzPBLndJkjTrGeiSNABsmasTA12SNHBOO23yy5Z6cGSgS9IA8By6OjHQJUkqgIEuSVIBugr0\niFgQEddFxPURcXybeU6qp18ZEfs1jf9cRKyMiKt6VWhJGjZ2tU/d0H9tLSLmAJ8AFgD7AIdFxN4t\n8ywE9szMvYBjgJObJn++XlaSpL4r9eComxb6/sCyzLwpM1cDZwOHtMxzMHA6QGZeAmwfETvXwz8A\n7updkSVp+HhRnDrpJtDnAbc2DS+vx23sPJKkSTLI1cmmXczT7cuo9exE1y+/RYsWrXs8MjLCyMhI\nt4tK0lCGIMaQAAAL3UlEQVT4zGf6XYLBt+++cM45/dn26Ogoo6Oj07qNbgL9NmD3puHdqVrgE82z\nWz2uK82BLkkaa+nSfpdg8P3hH8K22/Zn262N1fe85z0930Y3Xe5LgL0iYo+I2Bw4FFjcMs9i4AiA\niJgPrMrMlT0tqSQNMc+h90ZEuXXYMdAzcw1wHHARcC3w5cxcGhHHRsSx9TwXADdGxDLgFOC1jeUj\n4kvAj4CnRsStEXHUNOyHJBWt1BCaaSV/da2bLncy80LgwpZxp7QMH9dm2cMmXTpJktQV7xQnSQPA\nFvrUlV6HBrokaWgM9Tl0SVL/eVFcb5R8Dt1Al6QBYJCrEwNdkjQUSj8oMtAlaQCUHkYzxS53SdKs\nYLBPXal1aKBLkoZCpi10SVKfldqqnGkGuiSprwz0qXvgAbjnnn6XYvp0detXSZIG3eGH97sE08sW\nuiQNAG8s0zul1qGBLkkDoNQQUu8Y6JI0AEq+mEu9YaBL0gCwha5ODHRJGgCeQ1cnBrokDQCDvHdK\nrUsDXZKkAhjokjQASm1VqncMdEkaAAa6OjHQJWmAGOxqx0CXpAFgkPdOqXVpoEuSVAADXZIGQKmt\nSvWOgS5JA8RgVzsGuiQNAINcnRjokjQADPTeKbUuDXRJkgpgoEvSACi1VaneMdAlaYAY7GrHQJek\nAWCQqxMDXZI0VEo9ODLQJWkAlBpC/bBkSb9LMD0MdEkaAI1AN9jVjoEuSQPAIFcnBrokDYC1a/td\nAs12BrokDQBb6OrEQJekAWCgqxMDXZIGQKPL3WBXOwZ64ZYsgZUr+12K6XXxxbB6db9LIU2va6+t\n/l91VX/Lodkrss+HexGR/S5DySLgRS+Cr3+93yWZPhFw7rnwkpf0uyTS9IlY/9iPzPX1sTF1MZvq\nMCLIzOg8Z/dsoQ+BRx/tdwmmX7/fnJLUbwb6EBiGsBuGfZSkiRjokiQVwEAfAtHTszSSpNnIQB8C\ndkdLUvkMdEmSCmCgD4Fh7nKPgGXL+l0KSZp+BvoQGIYu94n28de/nrlySJo5w9xYGY+BLklSAQz0\nITAMR7HDsI+SNBEDfQgMe5e7JA0DA12SNJDsmduQgS5JGkgG+oYM9ELdey+ceGJ/y7BkCVx4Ye/X\n++lPw29/2376FVeU/etykiqeatuQgV6ob3wD3vSm/pbh0ENh4cLer/c1r4EzzthwXPMb+5WvhBe/\nuPfblfrpuc/tdwlmn7VrN36ZrbfufTlmCwNdkgbAvHn9LkEZdtyx3yWYPga6imM3nEo0mdZoqaby\nHt9kk6mvY7Yy0CVpAJQYQJM1lbpoXEhXYn0a6CqOV76qRLbQ15tKXTSCvMT6NNBVnBKPvKUSA2iy\nelEXJX5OGOgz7P77p7786tVw111T335rWRrDUy3jZDz6KDz44MTTH3po7PiNrYc77lj/Ru5mP++8\nE1auLPPNr/6Y7Pur1NfgQw/BihUbt8xUAr3Rg1fiAZKBPoOuuALmzp388suWVctvvnnnKzXHe/Pf\nf//67d9004ZluffeanjFiqmVcbLe9S54zGPaT3/d62C77caOb9RDN2/Oq6+GJzwBzj+/Gp47F267\nbeJlHvtY2Hln+MxnOq9f6uTWWyf//ioxgAC22gp22WXjlrHLfXwG+gz6zW+mtvxkW+UNjzyy/vGq\nVRtOe/jh6n8/WucA11478fSrrtqw/K0efbTzNn73u+r/LbesH3fffZ2XA7jhhu7mkyZy992TX7bE\nAJqsXtRFifVpoBdqGC4Ma+6F6CbQx1tOGhQlBtBkeQ59fAb6DJpNIdtaltlUtm41vyHXrJncct0a\nxPrR7DOV11GJATRZnkMfX8dAj4gFEXFdRFwfEce3meekevqVEbHfxiw7TKb6hvzpT0d7Ug4YW5bp\n+LCYzhCMaB/o9903usG8vdjXUj9MR0dH+12EgdCreprK62gQAmimXk92uY9vwkCPiDnAJ4AFwD7A\nYRGxd8s8C4E9M3Mv4Bjg5G6X1cbZmECfDQE0nWXI3PAN2dzl3hrorUp8I0+Wgd6d2VBPg/C6HYRA\nb3wuzYbPyF7r1ELfH1iWmTdl5mrgbOCQlnkOBk4HyMxLgO0jYucul5Umrdtz6K09BSW+kVU+X7fr\n2UIf36Ydps8Dbm0aXg4c0MU884Bdu1gWGJ5fxrr++ur/ZPd3yZINhydaz5VXrn980UXVvI2rxA86\naP0V9411NKYdddTUytjsxhu7X1fj507bzfujH204/V/+Bb71rfXT3/veaj+h+q5587yvfS3ssUf1\nVT2AD3wAvv/96vHRR3f3Yw0f+QgsXdp5vkHzi1/AT3/a71LMfr2qp8bP/k7m/dX8/p+tn5lTqaeN\n2afGt3I2djmAJz4Rbr4ZDjsMttiiGvfRj8JTn7px65mNIic47IuIlwILMvPV9fDLgQMy8/VN83wN\n+GBm/rAe/hZwPLBHp2Xr8R53SpKGTmb29EqjTi3024Ddm4Z3p2ppTzTPbvU8m3WxbM93SJKkYdTp\nHPoSYK+I2CMiNgcOBRa3zLMYOAIgIuYDqzJzZZfLSpKkHpiwhZ6ZayLiOOAiYA5wWmYujYhj6+mn\nZOYFEbEwIpYB9wNHTbTsdO6MJEnDasJz6JIkaTD09U5xw37jmYjYPSK+GxHXRMTVEfGGevyOEfHN\niPhlRFwcEds3LfP2ur6ui4gXNI3/k4i4qp52Yj/2Z7pFxJyIuLy+ENN6aiMito+Ir0TE0oi4NiIO\nsK7Giog31++7qyLiPyNiC+sJIuJzEbEyIq5qGtezeqnr+cv1+J9ExJNmbu96q01dfaR+710ZEedF\nxHZN06a3rjKzL39U3fDLqK6G3wy4Ati7X+XpUx3sDOxbP54L/ALYG/gw8E/1+OOpvkUA1Q16rqjr\na4+6/hq9LJcC+9ePL6D6hkHf97HH9fUW4CxgcT1sPY1fT6cDR9ePNwW2s67G1NE84EZgi3r4y8Ar\nracEeA6wH3BV07ie1QvwWuBT9eNDgbP7vc89rqu/AjapH39wJuuqny30ob/xTGauyMwr6sf3AUup\nPmjW3ayn/v9/6seHAF/KzNWZeRPVC+KAiNgF2CYzL63n+2LTMkWIiN2AhcBngcY3I6ynFnVr4DmZ\n+TmormXJzLuxrsazKfCYiNgUeAxwO9YTmfkDoPW3HXtZL83rOhd4Xs93YoaMV1eZ+c3MbNy25hKq\nb37BDNRVPwO93Q1phlJE7EF1pHcJsFNW3xQAWAnsVD/elQ2/+td8E5/m8bdRXl1+HHgb0Hx/J+tp\nrCcDd0TE5yPiZxFxakRsjXW1gcy8DfgocAtVkK/KzG9iPbXTy3pZ99mfmWuAuyOii9s7DaSjqVrc\nMAN11c9A92q8WkTMpTr6emNm3ts8Lau+lqGuq4g4CPhNZl7O+tb5BqyndTYFnkHVTfcMqm+e/HPz\nDNYVRMQOVK2fPag+UOdGdfOrdayn8Vkv3YmIdwKPZOZ/ztQ2+xno3dy0pngRsRlVmJ+RmefXo1dG\ndT986u6Y+katbW/icxvru3Ua42+bznLPsD8HDo6IXwFfAv4yIs7AehrPcmB5Zl5WD3+FKuBXWFcb\neD7wq8z8Xd3yOQ94JtZTO714ry1vWuaJ9bo2BbbLzDunr+gzLyKOpDpFeHjT6Gmvq34G+tDfeCYi\nAjgNuDYzT2iatJjqAh3q/+c3jX9ZRGweEU8G9gIuzcwVwD311cwBvKJpmYGXme/IzN0z88nAy4Dv\nZOYrsJ7GqPfx1oho3Jn6+cA1wNewrprdDMyPiK3q/Xs+cC3WUzu9eK/99zjr+hvg2zOxAzMlIhZQ\nnR48JDMfapo0/XXV5ysE/5rqyu5lwNv7WZY+7f+zqc4JXwFcXv8tAHYEvgX8ErgY2L5pmXfU9XUd\n8MKm8X8CXFVPO6nf+zaNdXYg669yt57Gr6M/Bi4DrqRqeW5nXY1bT4uoLkS9iurCo82sp4SqF+x2\n4BGq87dH9bJegC2Ac4DrgZ8Ae/R7n3tYV0fX+3Vz02f6p2aqrryxjCRJBejrjWUkSVJvGOiSJBXA\nQJckqQAGuiRJBTDQJUkqgIEuSVIBDHRJkgrw/wHbMhqv9fcK1QAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x10fa62a10>"
]
}
],
"prompt_number": 51
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" print len(current_data[current_data.values != correct_data.values])\n",
" print correct_data[current_data.values != correct_data.values][0:5]"
],
"language": "python",
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"132\n",
" year mm dd extent missing \\\n",
"1614 1987 8 3 8.537 0 \n",
"2045 1988 11 17 11.507 0 \n",
"2089 1988 12 31 14.196 0 \n",
"2096 1989 1 7 14.367 0 \n",
"2339 1989 9 7 6.959 0 \n",
"\n",
" source \n",
"1614 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2045 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2089 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2096 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n",
"2339 ftp://sidads.colorado.edu/pub/DATASETS/nsidc0... \n"
]
}
],
"prompt_number": 35
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def read_icefile(filename):\n",
" with open(filename, 'rb') as fp:\n",
" _ = fp.read(300)\n",
" data = np.fromfile(fp, dtype=np.uint8).reshape((448, 304))\n",
" return data\n"
],
"language": "python",
"outputs": [],
"prompt_number": 36
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" black_col = '#000000' # 0\n",
" red_col = '#FF0000' # 1\n",
" white_col = '#FFFFFF' # 2\n",
" blue_col = '#011892' # 3\n",
" brown_col = '#9d5717' # 4\n",
" light_gray_col = '#b7b7b7' #5\n",
" dark_gray_col = '#484848' #6\n",
"# dark_gray_col = '#00FF48' #6\n",
" black = 0\n",
" red=1\n",
" white=2\n",
" blue=3\n",
" brown=4\n",
" light_gray=5\n",
" dark_gray=6\n",
"\n",
" colors = [ black_col, red_col, white_col, blue_col, brown_col, light_gray_col, dark_gray_col]\n",
" bounds = range(0, 8)\n",
" \n",
" categories = mpl.colors.ListedColormap(colors, name='categories')\n",
" norm = mpl.colors.BoundaryNorm(bounds, categories.N)"
],
"language": "python",
"outputs": [],
"prompt_number": 37
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" data = read_icefile(sample_filename)\n",
" smmr_mask = read_icefile(smmr_mask_file)\n",
" ssmi_mask = read_icefile(ssmi_mask_file)\n",
"\n"
],
"language": "python",
"outputs": [],
"prompt_number": 38
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"a = (current_data.extent - correct_data.extent)"
],
"language": "python",
"outputs": [],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" print(mean(a[a>0]))"
],
"language": "python",
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"0.0108560606061\n"
]
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"outputs": []
}
]
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment