Skip to content

Instantly share code, notes, and snippets.

@rsignell-usgs
Created January 20, 2017 12:22
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save rsignell-usgs/e032db75e748cf5922d38d8be9e0ecef to your computer and use it in GitHub Desktop.
Save rsignell-usgs/e032db75e748cf5922d38d8be9e0ecef to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {},
"cell_type": "markdown",
"source": "# Extract time series at lon,lat location using Xarray"
},
{
"metadata": {
"trusted": true,
"collapsed": true
},
"cell_type": "code",
"source": "%matplotlib inline\nimport xarray as xr",
"execution_count": 1,
"outputs": []
},
{
"metadata": {
"trusted": true,
"collapsed": true
},
"cell_type": "code",
"source": "fname = 'http://thredds.ucar.edu/thredds/dodsC/grib/NCEP/GFS/Global_onedeg/Best' # Remote OPeNDAP Dataset \n#fname = 'my_raster_time_series_data.nc' # Local NetCDF file",
"execution_count": 2,
"outputs": []
},
{
"metadata": {
"trusted": true,
"collapsed": false
},
"cell_type": "code",
"source": "ds = xr.open_dataset(fname)\ndsloc = ds.sel(lon=230.5,lat=55.0,method='nearest')",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true,
"collapsed": false
},
"cell_type": "code",
"source": "dsloc['Wind_speed_gust_surface'].plot();",
"execution_count": 4,
"outputs": [
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEhCAYAAAC3AD1YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXecVOXV+L9nl7aw9C5Nmr2AQY2NbNQoWPNaI9GoSUzy\namJL1ZgIGqPJL5ZoiinGEhNjiUmMCmjU1ReJHQQVEVRgERCls5Rd2PP749zr3J2d2Wl3Zu7MPN/P\nZz4zt5+5c+c5zynPeURVcTgcDocjDKqKLYDD4XA4ygenVBwOh8MRGk6pOBwOhyM0nFJxOBwOR2g4\npeJwOByO0HBKxeFwOByh4ZSKI1KIyPsicmSx5YgCInKniFxTbDkcjkxwSsVRsohIi4iMysN5zxWR\nHSKyUUQ2ee8TA9uXiMgWb/1GEZnRzrk6icifRGSDiKwQkcvClrdQiMh3RGS+953fFZHvxG1/WkRW\ni8h6EZkjIicFtg0SkX+JyAfe7zY8xbVGeOdrFJG3ROSofH0vR7h0KLYADkcO5HPk7mxVnZhkmwLH\nq+ozaZxnGjAaGAbsAjwjIm+q6hMhyVlozgHmAWOAJ0Rkmao+4G27GHhLVVtE5CDgPyIyVlU/BFqA\n6cBPgdlpXOc+4HlgMnA88JCIjFHVNSF/H0fIOEvFEVlE5EARmS0i67we7m0i0sHb9iwgwDyv53x6\nocVLc79zgGtUdaOqvg38ATgvqwuKXCAii0TkYxH5p4gMDmxrEZGvi8g7IrJGRH6VzTXaQ1V/oapz\nVbVFVd8B/gUcFtj+hqq2BA7pgClTVHW1qt4OvEKKeyciY4HxwFRV3a6qDwPzgVPD/UaOfOCUiiPK\n7AQuBfoAhwBHAhcCqOpnvH32VdUeqvpg/MEicpinkNZ678HPa0Xk0HauPd5z5bwtIleJSHXc9r+I\nyIciMkNE9kt0AhHphVkn8wKrXwf2TufLx53rSKyXfxowGFgG/C1ut+OBTwHjgDNE5Jgk5zorxX0Z\nmqZYRwBvxp373yKyFXgBeEZVX0n7S8bYG3hPVRsD67K6b47C49xfjsiiqq8FFpeJyO+BzwC3BtYn\n7fWq6vNA7ywu/Sywj6ouFZG9gQeAZuBn3vYpwGvetS8FZorI7qq6Me48tZirbENg3QagexYyTQHu\nUNXXAUTkCmCdiAxX1WXePter6iZgk4g8gymXNm42Vb0Pcy9ljYhMw77/nXHnPtFTwEcDe2R5+lpa\n3zO85V2yPJ+jgDhLxRFZRGSs1/NdKSLrgeuAfvm+rqouUdWl3uc3gWswC8Hf/l/PLbNNVW8A1mO9\n9ng2e+89Aut6AJuyEGsXYGlAhkZgDTAksM+Hgc9bsMY5dETkm8DZwHGq2hy/XVV3qupMYJKInJDF\nJTbT+p5B9vfNUWCcUnFEmd8CC4DRqtoL+CHpxzIQkcMD2VvBl7/usNRniZ2unW2aaLuqrgdWAvsH\nVu9PnMsoTVYAIz4RRqQb0BdYnumJRGRKivuS1P0lIl8GvgccqaorU1yqA5akkClvAqO87+iT7X1z\nFBinVBxRpjuwUVW3iMgewP/GbV8FJE0pVtVZqtrdi7kEX/665xMdJyKTRGSA93kP4Crgn97yMBE5\nVEQ6ikhnEfku1rgnPBfwZ+AqEenlnesCAi4jL8CeLMssyF+B80VkPxHpjMVXXlDVhjSObYWq/jXF\nfUmoqETki5i1+Dnfkgts2927b11EpIOInI1Zb88G9ukMdPEWu3jLieRbBMwFrvbu8f8A+wJ/z/S7\nOgqPUyqOqBFME/4O8EUR2Qj8jraB6anAPV5w+TTC4ygsq2wT8CjwEHC9t607ZkGtxayEY4DJqroO\nPrEC5gfOdTXwHua6egb4mao+6e07FHPpBPcP8sm9UNWngR8BDwMfACOBLyTaN8lyGFyLJU28HLBq\nfuNtE+z3+BBYDXwLOENV5waO3wps9GR7G3PR2cEivw2cC+y7HQiswxToqS6duDSQcp6kS2yg2Vew\nHPn5wPmq2lRcqRwOw+v576WqPyy2LA5HWJStUhGRXYBZwB6q2iQi9wOPqeo9RRbN4XA4ypZyTymu\nBrqJSAvQFQt2OhwOhyNPlG1MRVVXADdig8Q+ANar6n+KK5XD4XCUN2VrqXijmU/G0jA3YLWDpqjq\nX+P2K0//n8PhcOQZVW2TSl+2lgo2ovc9VV2rqjuxrJmEZTlUteivq6++uugyRE0WJ0d0ZXFyRFeW\nQsmRjHJWKsuAT3t584KliS5ItvN998EmN17X4XA4cqJslYqqvoSNL5iDFaMT4PfJ9r/hBliQVOU4\nHA6HIx3KNqYCoKrTsPksUtLcbK9iUVdXV7yLxxEVWZwcbYmKLE6OtkRFlmLLUbbjVNJFRFRVGTsW\n/vAHiMhz4XA4HJFGRNAKC9RnRLEtFYfD4SgHnFLxaGpySsXhcDhyxSkVD2epOBwOR+44peLhlIrD\n4XDkjlMqHk6pOBwOR+44peLhlIrD4XDkjlMqgKpTKg6HwxEGTqkAO3bYu1MqDofDkRtOqRBTJk6p\nOBwOR244pYJTKtlyyy1wj5tH0+FwBCjr2l/p4pRKdixZ4io7OxyO1jhLBRtND06pZEpzM2zbVmwp\nHA5HlHBKBWepZEtzM2zdWmwpHA5HlHBKBadUsqWpyVkqDoejNU6p4JRKtjj3l8PhiKdslYqI7CYi\nc0TkNe99g4hcnGhfp1SywykVh8MRT9lmf6nqO8B4ABGpApYD/0i0r1Mq2dHcHBs46nA4HFDGSiWO\no4F3VbUh0UaX/ZUdzc3Q0lJsKRwOR5SoFKVyJnBfso3OUsmOpiZnqTgcjtaUbUzFR0Q6AicBDybb\nxymV7HAxFYfDEU8lWCqTgVdV9aNkO9x551QA5syB+vo66urqCiNZieOUisNROdTX11NfX59yP1HV\n/EtTRETkPmCGqt6dZLs+8ohy0knwP/8DDz9cYAFLmIMPtjItb71VbEkcDkehERFUVeLXl7X7S0Rq\nsCB9u6qiqQmqq537K1OcpeJwOOIpa/eXqm4F+qfar7kZunVzSiVTnFJxOBzxlLWlki7NzdC1q1Mq\nmZJOmZYzzsj9vl56KSxdmts5HA5HYXBKBadUsiWVpaIKDz0Emzfndp3nnnNKxeEoFZxSwSmVbPGV\nSrJcD39brpWMm5pcNWSHo1RwSgVrtJxSyZzmZlMafkWCeLZssfdc4y6uGrLDUTo4pYIL1GeLf7+S\nNfiNje1vz+Q6zlJxOEoDp1Rw7q9saWqCLl2SKw1nqTgclYdTKjilki3NzdC9u7NUHA5HDKdUcEol\nW1IpFWepOByVh1MqOKWSDS0t9urWrTBKxVkqDkdpUDJKRUS65evcTU0uUJ8pzc3QsSPU1Dj3l8Ph\niBF5pSIih4rIW8ACb3l/EflNmNdwlkrmNDVBp075D9S3tNicLc79VZqsXg0//GGxpXAUksgrFeBm\n4FhgDYCqvg5MDPMCTqlkjm+p5Fup+L+Js1RKk0WL4N57iy2Fo5CUglIhwTTAO8M8v1MqmdOeUlm+\n3AZFhuH+SjUWxhFtNm6Ejz8uthSOQlIKSqVBRA4FVEQ6ich38FxhYeEGP2ZOe0rlyCPh9dfDsVT8\n0frOUilNNm2y58B/FhzlTykolW8AFwFDgOXAOG85NJqaoHNn+7wzVBuofEmmVJqb4b33YM2amKWS\ni0Jwlkpps2mTvX+UdN5VR7kReaWiqh+r6hdVdaCqDlDVs1V1TZjX8BvIjh2dtZIuyQL1S5aYYl6/\n3nqnvXrlx1JpbobvfS/78zoKg1MqlUfklYqI3C0ivQLLvUXkT2FewymVzElmqSxebO++UunbNxyl\nEn+O9evhN6HmADryga9UXFylcoi8UgH2U9X1/oKqrgPGp3OgiPQUkQdFZIGIvCkiByfar7nZet1O\nqaRPMqWyaJG9r19v7q8+ffKT/dXcDNu3Z39eR2HYuNHenaVSOZTCdMJVItLbUyaISB/Sl/uXwOOq\nerqIdAC6JtrJWSqZk2zw4+LF0Lt3/i2VpiYbv9LSAlWl0DWqUDZtsnilUyqVQykolRuB2SLykLd8\nOnBdqoNEpDtwhKqeB6CqO4CNifZtanJKJVOClsr69bH1ixfDhAmwbp0plf79c7dUamoSWyr+u59k\n4YgemzbByJFOqVQSke/jqeo9wGnAh8Bq4BRV/XMah44CPhaRO0XkNRH5vYjUJNrRWSqZ47sME7m/\nJkxo6/7atg2uuSbz6zQ1QY8ebRWT/zs5F1i02bQJRo1yMZVKohQsFVT1TRH5COgCICLDVXVZisM6\nAAcAF6nqKyJyC/AD4Or4HZcsmcpdd5n/d9asOkaOrAv3C5QhvnUXVCrNzbBsGRxwAMybF3N/vf8+\nfPAB3HIL/PjHmV+nZ8/W1pC/PvjuiCabNsE++0BD/PBlR8lRX19PfX19yv0ir1RE5CTMBbYLZqmM\nwAY/7p3i0OVAg6q+4i0/BHw/0Y4DBkzlootg9mxrEB2piQ/Ur1sHr7wCgwfDwIFtLZUNG7Ibr9Lc\nbJbKypVt14NTKlFn40azVF57rdiSOHKlrq6Ourq6T5anTZuWcL/Iu7+Aa4FPA++o6kjgKOCFVAep\n6ofYaPzdvFVHAW8l2jfo/nKNVHrEK5Wzz4YLLoCTTmodqPeVysaN9t7Sktl1nPurtHHur8oj8pYK\n0Kyqa0SkSkSqVPUZz5WVDhcDfxGRjsB7wPmJdnKB+swJKpWtW81KefVVGDrUan+tX2/r/eyvDRvs\nuG3brM5aJteprbVMr507obra1jv3V2ngKxUXqK8cSkGprBeRWuA5TEGsBhrTOdCraHxgqv1coD5z\ngoH69983C2TIENvWq5cplR07WlsqYIomE6USP3K/W7fY9f3tjuiyaRMMG2bv/v/MUd5E1v0lIn6i\n6MnAVuAyYAbwLnBimNdySiVzgoH6N96AceNAxLb5s0E2NcXKtPiWSqZxlaBSCR7rLJXos3On/WY9\neljnYk2oxZUcUSWySgX4r/d+u6ruVNUdqnq3qt6aj9pfbkR9ZgTdX6owPlDjQMSUSdeusTEmQUsl\nm+vED7LMd0zFFbDMnc2b7RmoqoJ+/VxcpVKIsvurk4hMAQ4VkVPiN6rqw2FdyFkqmRNUKmCWSpBe\nvUyR+MogV0slfgBkPt1fixfDlCnw0kvhn7uS2LQJune3z06pVA5RVirfAL4I9KKtu0sBp1SKSDCm\nAomVyo4dsVhIrpZK/CDLfLq/Nm6MKUFH9gSVSo8eseKSjvImskpFVWeJyGxguaqmLMuSCy77K3OC\njX2XLrDbbq239+plCsRXBvmyVPLh/mpuNoXoyI1Nm0yZgCkXp1QqgyjHVFDVFqxES15xlkrm+Ip4\nxAh48knoENc98WMq/no/SJuNUim0pdLc7J6DMAhaKt27x6xVR3kTaaXi8ZSInCri5xaFj4iNf3BK\nJX18RSwChx/ednvv3rH03y5dYPVqG7OSjfur0DEVp1TCYePG1krFWSqVQSkola8DDwLbRWSjiGwS\nkVD7PH7uvFMq6ZNqzIFvqYAplQ8/tPItYaUUO/dX9Im3VJxSqQwir1RUtbuqVqlqJ1Xt4S33CPMa\nTqlkjm9BJCNeqXz0EQwaFF5KsXN/RR8XqK9MIhuo9xGRiYnWq+pzYV3DKZXMydRS2bEjM6XyxBMW\n/C9GSvGOHc5SCYP4QL2LqVQGkVcqwHcDn7sABwGvAkeGdQG/x+2USvr4AfRknHoqTPS6A37a8YAB\nVmQyHX73OzjhhOSB+ny7v9xzkDvO/VWZRF6pqGqrMSoiMgxIt6BkWjhLJXNSWSqDBtkLTCF07WoN\nS7qWyubN9koWqM+3+8tZKrmzcaN1JMAplUoi8jGVBCwH9gzzhE6pZE4mxQG7dLGJthJNC5yMxkZT\nKvEFJYPXr6nJn1Jpacm8TL+jNc5SqUwib6mIyG3YCHowJTgOCHXKH6dUMidVoD5ITY351mtq2k62\nlQxfqQQD9fGWSrdu+XN/gVkr6X5HR1s2brTOBNjv72IqlUHklQrwSuDzDuA+VX0+zAsElUq6Pv9K\npxiWSrCn68+zki9LxX93SiV7NmyIKRVnqVQOkVcqqnq3/1lEegPDwr6GC9RnTqpAfZAuXWKWSjZK\nJZGlUgil4uIqubFxoyvTUolEXqmISD1wEibrq8BqEZmtqpelcewSYAPQgs0geVCi/Zz7K3OysVS6\nds0tUB8/TiVf7i9fmbhnITecpVKZlEKgvqeqbgROAe5R1YOx+ebToQWoU9XxyRQKWP0qaKtU1q2D\nffax+UIcrclUqWRiqai2dX/16QPvvBP7LQrl/nJkT9BS6dzZEh/yNf+NIzqUglLpICKDgTOARzM8\nVkjjOz7wgL3HK5X33oM334RlyzK8agWQSbwh05hKU5PNGhgM1J94otUPmz49to9zf0WboKUi4kbV\nVwqloFSuAWYCi1X1ZREZBSxK81gFZorIyyJyQaqd45VKQ4O9z52bocQVQD4tlcZGe9+0KWapdOwI\nN94Il19uCqe5Of/ZX85SyZ6mpljat49zgVUGkY+pqOqDWEFJf/k94FR/WUSuUNXrkxx+qKquEpH+\nwJMiskBVZ8XvNHXqVADefx9eeqkOqANiSmXOHDj55Ny/SzmRiVI5+2zbt6kpM6WyeXOsejTAccfB\n+edbHTHf/ZWPRspZKrnjpxMHa4s7pVLa1NfXU19fn3K/yCuVNDgdSKhUVHWV9/6RiPwDK/GSVKk0\nNUH//tZo9e9vSuWAA5ylkohMsr/239/eFy5MX6l06xab49x3s4nY+q1bC+P+cpZK9gTjKT5OqZQ2\ndXV11NXVfbI8bdq0hPuVgvsrFQnnWRGRriJS633uBhwDvNHeiTp1gs9+1ooZgimVE080S8XRmkws\nFZ903V+bN1uJl2BKsU/XrjaWyHd/5augZPDdkTnBeIqPKypZGZSDUkmWmzUQmCUic4AXgH+r6hOp\nTjZpUiwY3NAAdXWWBbZ2bUjSlgnZDAysqUlvcGljo8290tgYi6n4BJVKba2LqUSVDRvaWiouUF8Z\nlINSSWipqOr7qjrOSyfeV1VvSOdkkyebpdLSYkplxAjYbz/nAosnn5ZKY6M1QJ07W+MUVCr+OVz2\nV7QJlmjxce6vyiDySkVEDkux7sH47bkwYgT06wcvvwyrVsGQIbD77rB4cZhXKX2yVSrbtqUe9+PH\nVGprzUpM5v5yMZXokshScUqlMoi8UgFua2+dqv407AtOmgR33mkD7jp1gsGDTcE4YmQSqPeproYO\nHVK7rIJKZefOxO6vQhWUdGRHMkvFxVTKn8hmf4nIIcChQH8RuTywqQdQnc9rT54Mp5wCe3oF9gcN\ngrfeyucVS49siy367it/4q5EbN4cUyrQWnn5xztLJdoki6l89FFx5HEUjihbKp2AWkzxdQ+8NgKn\n5fPCRxxhMZVhXunKQYOcpRJPNu4vSK/+V9BSgfYD9U6pRBMXU6lcImupqOqzwLMicpeqLgUQkSqg\n1qsFlje6dLHUYl+pOPdXW7JVKukE6xsbTWG0p1QKUVDSub+yZ8OG2P/HxymVyiDKlorP9SLSwxtr\n8gbwloh8N9VBufKDH8BZZ9lnZ6m0ZvlyUw7tubCSka5S6dbNGiERi8X4FMpS6dzZWSq5kMhSyVcF\nBEe0KAWlspdnmXwemA6MBM7J90UPPxwOPtg+DxxoMxa6asXGjBlwzDFQlcXTk4lSqa1tG7cpVEpx\nTY2zVHIhUUyltjZWgsdRvpSCUukoIh0xpfKIqjaTfMBjXqittawll7liTJ9uyQzZUFMDTz1lhSGT\nEVQq8S62+BH127eHr+ybm+06zlLJnkSWSrduTqmUC4cemnxbKSiV3wFLgG7AcyIyAgvWFxQXVzGa\nm00pHHtsdsfX1MC0afD008n3CWZ/xVsqQaVSU2PusZ07s5MlGc5SyZ1ElopTKuVBczO8+mry7ZFX\nKqp6q6oOUdXj1FgKfLbQcri4ivHf/8Lo0eYSzIaaGruX7fnWU7m/GhtNkVRXW+wjbBeYr1ScpZI9\nhbBUtm93qf7FYOlS2GWX5Nsjm/3lIyI/TrLpmkLK4ZSKMX8+HHhg9sdPmQIXXxxLgkhEMPsrkftr\nwwZbL2JKZ/t2Wx8WO3Y4SyVXEhWU9CtPh8Xf/w7f+pZVu+jdO7zzOtpn8WIYMwaWLEm8PfKWCtAY\neO0EJgO7FlqIQYMsWF/pbN2aWwN+5pnmj20vPtWepeIrFX99p075sVRcTCV7VBOXvg/bUlm0yJ7H\na68N75yO1PhKJRmRVyqqemPgdR02g9boQsvhYirGtm2tZ/PLhi5dzApI1mgHU4rjLZWampilAvlz\nf3Xt6iyVbHjgAdhjD7t/iazMrVttYHEYLF4MU6fC3Xe7Dl8iLrsM3nkn/PMuWgRjxybfHnmlkoCu\nwJBCX9S5v4xUJVbSQaT9gXCpAvVBpeK7v8LExVSyZ948q5339tttt1VX27OTTqXqdFi82Kzegw6y\nArCO1jzwQOLfIVdK3lIRkfkiMs97vQksBG4ptBxOqRhhWCrQvlJJlVK8fr1zf0WVLVus0neyQG6Y\nLjC/xzxuXGxqinnzwjl3qbN5M6xYAWvWhH/uVEol8oF64ITA5x3Ah6pacMdEv37w8ceFvmr0CMNS\ngeQVa1VjSmXvveHcc1tv991fvr8+X0rFBeqzw//tkhGWUlm3zizUAQNg/Hi47z5Ytsymrt640Z6v\nSsZ3e4WtVHbssOyvUaOS7xN5SwXYFHhtBXp4gyHTQkSqROQ1EXkkFyHSKYRYCYRlqSSbBXD7dhto\n2rEj9O9v2T1Buna1dOJ8x1Sc+ys7tmxpP5EjrAywd981K0UkZqn4M7a++27u5y91Fi6097CVyrJl\nNpygvY5lKSiV14CPgHeARd7nJZ6i+FQax18C5JzN7g+6S8UHH8D3vtf+Pr/6Ffz737lKVBy2bs2f\n++v++600TnwqahC/wQq6v8KOqbiU4uxJR6lkYqk89hj88Y9t1y9aFHPBjBljXoT77rPGbtGizGQu\nRxYuNO9K2NOgp3J9QWkolSeB41S1n6r2xVKKHwUuBH7T3oEiMhQ4DkjwWGZGulPhvvcePPpo+/vM\nmgUXXVSals+2beG5v+KVyiuvWHXoF19Mfpyv0IKBehdTiQ6p3F+Z1v+aORNuv73t+mDjVlVlU34/\n95yNf3KztJr765BDwrdUUmV+QWkolU+r6kx/QVWfAA5R1ReAzimOvRn4LiHUCkvXUmlsTG3eNzZa\n7/rGG3OVqvCEaanEx1QaG+2BHTky+XG+QvMtlXy6v5ylkjlhWyqLF1tJkA8/hNWrY//B+MZt/HiY\nMMEaUmepmKWSD6WSjqVSCoH6lSLyfeBv3vKZwIciUg0kzXgXkeOxoP5cEakDJNm+U6dO/eRzXV0d\ndXV1bfapqUlPqWzZklqpbN5spfV//Wu46qrU54wSYVkqiWIqqRoksF5pTU3+U4qdpZIdYQfqFy+2\nhI3HH4dbboGTTrLBjq+91jredsYZNrnegAHw5z9nL385oBqzVO67L7zz1tfX89hj9RxwgI0PSkYp\nKJUpwNXAP73lWd66auCMdo47DDhJRI4DaoDuInKPqn4pfsep7d0hD78sSKrJqdJRKo2N1qtqaMh+\nsqtikc+YypYt7TdIPsGBdWG7v/zilJ07O0slG8IM1O/YYYHhG2+E73/ffvPp0+FrXzPL5YADYvtO\nnGjvDQ3OUlmxwn6D0aPDjanU1dVRXV3HVVfBPvvAtGnTEu4XefeXqn6sqt9S1fHe61uq+pGqNqnq\nYhG5LclxV6rqcFUdBXwBeDqRQsmEdFxgjY2mKNpr6BobrVbRLrtYel4pkc+YSmNjeiVgampau7+2\nbctdHp/mZss+69DBWSrZEKb7a+lSGx928snWOD78sFkud99t8/kEJ2/zGTLEUs7DrDFWarzzDuy+\nO/TtG677a+dOeP99U1btEXmlkgaHFepC6SgVf3t7D7XvIhg7tvR6VfmMqaTj/oLWlsqAAdZrDYsd\nO+zcHTs6SyUbwnR/LV5s/5GhQ633ffDBcNRRcMMNNmo/EVVVNoaiktOKFy40pVJTY66wdNz26dDQ\nYGn+qf7/5aBUUqKqz6rqSbmeJ91ZCyE9pTJmTOllquTbUsnU/TVsmD3sYeG7Izt2dJZKNqTqGGSS\n/RVMGx4wwN4nTbLj25vPZ8yY0uushcnChbDbbuauD9NaSSdIDxWiVMIiTEultrY0lUpYlkq2gXpo\n7f4aPjw/SsW5vzKnpSV1pyNTSyW+ETvpJJs11FcyiRgxItxnotTwLRUwpRJWXCWddGIoD6WSNKsr\nbDJRKsnqWu3cGRuVXunurzAC9fm0VJz7KzP8Z6OqnVYlE6WSqBEbODB1Kn737i6mElQqzlLJnF8W\n6kJhuL/83rhIaVoqYbq/Eo1TSTem4lsq+VIqzlLJnHQszXSyv1591dxc//d/6fWM46mtrVylsn07\nLF8eq81VDKUS2ZRiEfk37Qxa9GMkqnpXoWQKw/0VjBuMHGkZLjt2WCMWdfxGNowU6GSWSrruL1+G\ngQOtanFYys5ZKtmTTqcglaXS0gJf/zqceipccYXNzZIptbWV6/56911zCfv/jz59wlMq6bq/otyU\n/cJ7PwUYBNzrLZ8FhJjvkz7pKpUuXdJTKl26WMrk0qWp0/SiQFgNN7SNqfhZKplmf1VVWWr28uXp\n9aJS4VKKsycd92UqpXLvvXbvf/ADs+azoZItlWA8BcKLqfjpxO1VJ/aJrFJR1WcBRORGVZ0Q2PRv\nEXmlGDKl6/4aODA9pQJmrSxZUhpKJax4CrS1VPzqxOlYbEH3F8RcYGEoFZdSnD3pWCqpsr/uuguu\nvDJ7heJfo1KVSjCeAqZUwpgV84MPzOpJJ+ZZCjGVbiLyiX4UkZFAGl8tfNK1VAYMSP5Q+7Ma+gwY\nAB99FJ6M+SRMS6W21u6VP7VsukF6aO3+gnDjKi6mkj1hWCqbN5v1nguVHKhPZKmE4f5K1/UFpaFU\nLgPqRaReROqBZ4BLiyFIuiPqU1kqtbWx5f79S0ephGmp+DW8/PuUbpAe7P51DpQSzYdScZZK5qQb\nqG9PqaSQCcE3AAAgAElEQVQ7Vqk9KtlS8ceo+IQVU0k3SA8Rdn/5qOoMERkL+CG7t1U15BKC6ZHO\nRF2pLJX4P00pKZUwLRWIucB69Eg/ngI2bUCQYcNg/vxwZHKDH7Mn3UB9ew1+Js9BMmprk6f0lzv5\nslQyUSqRt1REpCtWvv6bqvo6MFxETkhxWF5Ip1Kxr1RSzb/uE+VpiuO/a5iWCrQO1mfSQ+3f314+\nYQ6ADLq/nKWSGWG4vzKxWJNRqZbKmjX2/A4cGFsXVqC+3NxfdwJNwCHe8nLgJ8UQJCz3VylYKo89\nZnNUBAlrKmGfHj2s+B/k1kPdZZdwgpEQy/5ylkrmpKMQuna156glyaQVmcTWklGpSsUP0geTHJyl\nkpjRqvpzoBlAVbdSwFH0QdIN1LenVOID9VFUKs3NVgpjxYrW67duDdf91atXa6WSbWMSZmDWz/5y\ngfrMyWQ+nET/I9VwrOFKVSrx8RSwaujr1iVX4unQ0mIz2qaboVoKSqVJRGrwBkKKyGigKDGVVCnF\nO3dayft+/UrbUvnDH6x+0pYtrV1AYVsqvXrZwEXIze0RZiPiAvXZk27HIFlcZds2SxVPVNI+E2pq\n7Fz+3DiVQnw8Bew57tYt1nnLhg8+sP9qMMGoPUpBqVwNzACGichfgKeA7xVDkFSWit9Ta6/nnCj7\nK2oxlfvug8sug549Wz+M+bBUfKWSi6WSD6XiLJXMSbdj4Pee4wkjSA9mDWU6w2Q5ED9GxSfXuEom\nri8oAaWiqk9io+rPA+4DJqhqfTFkSVeptNfIxVsq/g+ei3kaJuvXw9y5UFfXutGH/FoquTQofgOi\nSYv6pI+zVNqnPUs93Y5BMus8jCC9TyW6wBJZKpB7XKXslIqICDAZ+JSqPgp0FZGDiiFLKveXrzAy\nUSodO9r+iXpuxeA//4HDD7fv2qtXa7nyaank0qBUV5tcYUxG5CyV5Kxf336ZjnR/w2QZj2EE6X0q\nbQDkzp1W9ytR45/rWJVMMr+gBJQK8Bss8+ssb3kT8OtUB4lIZxF5UUTmiMh8Ebk6V0HyYalAtOIq\n06fD5Mn2uXfv1pZK2CnFQaWVa4MSVs80mP3lLJXWrF4Nq1alrsCdimTPe1juL6g8S2XpUruvif5D\nzlJpy8GqehGwDUBV1wGd2j8EvAGSn1XV8cA4YHKuFk4mSiXZOJX47C+IRlzl9dfhggvgn/+MKZVE\n7q8wLZWg0srV9RFWIxKs/eUsldb4fvlk0zenO9aoPfdXWJZKpSmVZK4vSC+m0tAAt9+e2IVcjkql\nWUSqiWV/9QfSikCoqq8COmPVA3Lyuqfr/urcOZYJlmyfIFGwVP7zHyts+ZvfxB6geKWSD0sljEA9\nhGupdOxoLrWdO8OJ05QLfsO0alXi7c5SKR7JgvSQnqXyxBPwv/8LX/lK685US0tyt1oySkGp3Ar8\nAxgoItcBs4CfpnOgiFSJyBxgFfCkqr6ciyDpWioiyauxRlWpNDbCwQfDmWfGBk/5jb6q9eCjGqiH\n8JWKiBtVH4+vVBINNN2xI3elEnagvpJKtSQao+KTTkxl6VL4znfMxXniibF7t3Klxae6d09fllKo\n/fUXEXkVOMpb9XlVXZDmsS3AeBHpAfxTRPZS1bfi95s6deonn+vq6qirq0t4vlRKJagw/EBh795t\n94nP9+7XLxpKJV5Wv9F/7jn4yU8sWBfFQD2Er1QgplTCmJSsHEhmqdx4o7lNE7l2E1GIQH2lWSoL\nF8LJJyfelo6lsmwZTJwI119vtfUmTrSqGkHXV319PfX19SllibxS8egK+C6wjPvKqrrRq3A8CWhX\nqbRHKvdXsKeW7KFOZqkUe6a6xkYYOrT1ut694e23zbR+8UUYMqRy3F8Qi6uE+Z1LmTVr7J4ElcqH\nH1pD1L+/NWzO/VUcco2pLF1qA547dLDYyvXXw6GH2gycfuZXfId72rRpCc8VefeXiPwYuBvoA/QD\n7hSRq9I4rp+I9PQ+1wBHA2/nIkuXLjaZVHBMySuvwEsv2ed4peLPwf7//h9ccgnMmpU8UF9sSyWR\nXH52VkODmcNvvVW+lkpLC/z2t20tlSgE6xcuhGeeSX//3/8+P+Oe1q61hiuoVH70Izj3XLumarTc\nX5k8D//6l80eGmU2bICf/cxGuAfZvNkU/vDhiY9Lx1JZujR2vIhNlHbttXDrrZlPfhd5pQJMAQ5U\n1amqejXwaeCcNI4bDDwjInOBF4GZqvp4LoKItLVWfvtbe0FrK2SvveC118wn+dOfWq9/ypTEjfew\nYWZ+FpNEFpTf6Dc02HefMyfcXntNjQXDt20rvqWyahVceKEpTn/2yaikFT/yCPzxj+nt29xsc7zn\no5Oydi3svXcspvL669YY/+hHcMQR5i7ZZZfU5/HdvfFJEMV0f11/vc06GWWeegp++UvYd1/44het\nQws2jmTMGKskkIhUMZWdO01RDRvWev0558ALL8D552cmZykolZVAsH/cGfggyb6foKrzVfUAVR2n\nqvup6nVhCBM/p8rcufaC1pbKpEkwYwbMnAmf+xxMnWqB8ETulDFj7MEoJqmUyqc/bQ1smJaKSCyt\nOFfXR/z0xJniux9nzIiepbJyZfpKwk/3zZdS2WsvU8CqVnT06qvtOQE47rjkDVuQrl3t3sb/XmG6\nvzId/Lh4sY3RijJz58JXv2rFHQ84wFxThx9uVmIy1xektlRWrTLFk+i//alPpddRCFIKSmUD8KaI\n3CUidwJvAOtF5FYRubXQwgQrrDY3w4IF5p5oamr9pzjmGKivt57cpEm27uc/hwkT2s6/vcsuZtoW\n0wecKIEgqFROPNHWhR1f8K9RbPdXQ0Os0GEwphK0VH73O2tEC82qVekrCd81lW+l8t//2j372tey\nO1eisVnFGqeybp1Zy/PnhzP3SL6YMwfGjbP/zLe/bam+l1xich9+ePLjeva0jnCyDtKyZcldZ9lQ\nCkrlH8CV2DTC9cAPgX8Br3qvghLMAFuwwIJbo0aZ2yT4p+jbF/bcs7VSGTkyFn8JUlVlZaXffbcw\n3yERySwVP6biK5UwLRX/Gr6lUkz3V0MDnHaajU+JD9T7LF1qr0ITJaWy55527kcegdNPj7kKMyVR\nXKVYgfrFiy0d94gjbLxWVJk7t/UcRx062G8waxZc2s4E675HIJnC9IP0YRH57C9Vvdv/LCK9gWGq\nOq9Y8nTtGntY/Z5DdbV93ry59Z9i8mTrAaVjPo4ZYw/3/vvH1s2aBYMGZR4oS4fHHjPTdtAgW06k\nVHr3tkaqZ0/zpffqFX1LZdkysyi6dIErrog1ev/3f5bdNnJk4uMbGmCffSzjJT6l2Gft2uL0ZFet\nsl69alsrN9G+kB+lsmaNPS89e8Jf/2qvbEmkVAodqH/hBXvG/dpWRxxhLrAzzghHhjD5+GNzF+66\na3bH+3GV4KyQPsEgfRhE3lIRkXoR6SEifYDXgD+IyE3Fkuczn4E//ck++z2HceOsh/PYY3DYYbF9\nv/Y1+NWv0jtvorjKjTfCbbeFI3cQVRs9e889sXWJlErXrqYwhw61xuzee62nGia9epn5XlVlD362\n+I3I44/D009bZYCFC2Pbb7oJHnww+fENDRaovPFGOP54WxdvqaxZE84sepmycqW5V9Ppea9cab9V\n2GV/du60bMZevUyxbNxocbZsSWapFMr91dhoMYmf/zw2FuOQQ6xzGEXmzrUOZ6pORTLai6ssWxau\npRJ5pQL0VNWNWPn7e1T1YCw9uCj86Efw0EPwxhsxS2X8eOu1feELsMcesX0HD27f1xlk7Fh7uIPk\nK3i4YIGlT86YEVuXKCtNxBoRPyvk+ONtEqUw6dUL/vY3i0GlE+RNht+INDRYwHjChNZKuqGh7f0N\n4iuVAw+Mfd8oWCp+ZtywYelZH6tWWQMZtqWyYYMFv6urTal87nPZu74g8QDIQrm/VC3Nf9dd7T/g\nZ0/53oIoluaZM6ft9N6Z0N5YlbDdX6WgVDqIyGDgDODRYgvTp4+NLt9vP3j5ZcvCGD/eltMcQ5kQ\n/4H28WvurF1rD31dXXi9qOnTLSXx5ZdjGTjJgqRBpZIPevWy1Gu/iGW2BJXKsGFt72dDQ0zJnHaa\nxQSCLFvW9nvGWyr5UCo7dliK6JIlibd/+KG5LAYMSM/6WLXKzhe2Ulm7NmZJ7rOP9fJzYeDAtiPz\nwwzU9+lj9y5eQdxxh3VebrvNOoJdutj/YexYc+t17Zq8tlkxmTvXOrDZkspSqSj3F3ANMBNYrKov\ni8gooKgJuF/7mjX6jY3W4+rTx3L2+/bN/pzx7q8VK+whP+EEG1z27LOt3Tm5MH26NawHH2y57zt2\nJE8XLoRSAbNUcqE9pbJ9u9U0WrzYfreZM21my+3epNRNTdZgDx7c+pzxKcW+UgmzJ7t0qVm93/9+\n4u2rVpllkG4pn0IolZtvNqs8F0aNstTYIGFaKsOHW6dgQVxBp1mzLOa2Zo31zidPts9+3DKRxyAK\n5GqptDdWpeIsFVV90BtncqG3/J6qftJPEpEriiddeAwdan9cvwil7+edPNnSNw85JBw/eWOjlVw5\n8kg79/TpsR5iIn9t7975VSq9e5urasCA3M4Tr1TGjo0p6eXLY+6j+fOtgd57bzjvPGsglyyxnnO8\nO6djx9aVpn2FEmbq9+LFlhzw3//C7Nltt69aZcouGINobo5ZWuvWWeq6z8qV+VEqa9bk1mmKJ96S\nhHAD9SKx5/vjj61TBm17/JMn27PvB7CjMGYsni1b7BnNJZ6ZzFJZv946Wn7nLgwir1TS4PRiCxAG\nVVVWW2vFClv2M1JOPBH+/nfryYfRUMyfbwOlune3xjw+FTqe730vlhKdDyZPhl/8Ivfz+GVxli83\nBR1stBoazH8+cqTdy3HjLJA/dizcfz9897uJFeenPx3LcGpqsvjG4MHhusAWLTLX6aWXJh7RvXKl\nWSpBpVJfb8UDn3vOLBx/rIhqzFIJO1AftFTCYMwYc+8Gy8mEGaiH2ADkyy+HL33JLNOFC81953P0\n0eYS8ztUiZRdsXnjDYvV5hLPTBZT8YP02SYAJKIclEqIt6O4BBsO31Lp2hVOOSW8+mCLF8cKxPnn\nbE+pHH104jTEsBgxwjLqcqV7d2tQa2vtng0fbj71bdtau8QefNDcCEOHwjXXWAbc448nVipXXAGP\nPmq9W79RTac4Xyb4v/Nxx1mvOt615ru/goMFp0+3DsGXv2zjoD74wGJjmzZZ52TEiFgKcliErVRq\na82963eiIFz3F8BRR5n195//WPba3/9uHYvgNTp3tukefKLo/vITgnIhmaUStusLykOpRDBXIzsS\nKZVE23IheN50lEqp0K2bNaLBzK0RI8xvH3SJvf126z/obruZpRIcH+TTs6cV1rv++tZKJcy0Yl/J\n7767yfzmm3DSSTYe6OijrdGNt1RmzIBf/9oU509+Yq68efNiCqhzZ2s4gxOsZcv69fa8XHZZ+G7Q\neKsgTPcXWEfjmGOsCOPxx8MNN6RunKPo/oof9JgN7SmVMIP0UB5KpWwslWAwNmhR+NvCcGn46ZNg\njeSGDfaKL9FSanToYIkGwYbPb7SClgq0/YP+9KdmlSTiM58x98PatfbH7NMnfPfXmDHmfpg0Cb75\nTfOff/SRufPuv99cbv6zsXSpPQcTJliSxQUX2PeZO9dcZX6yQVidkGuvtczDTZvMFRomwbiXP1Nq\n2BUbHn7YCiNOnmyu31SNcxTTisOwVJIF6sMeowLloVTaGdJWWvguDlV7sEePbr0t2Ei8/771ajMl\nqKyqqy1A59e9KnVqa1srlbFjbS4YP1147Fj7c8XPGwPJfcpjxpi18/HHdmyYSmXHDlMSo0bZ8uTJ\nFlC++Wb7LrfcYoolaKlMnw7HHmtuLl/mceOs4fHjL2D7r1xpxQZvvtmy3zJl0SK4+2647jqznML0\nu0NrSyU4a2qY+Oc76ihLvEjVOPfqZYrNL8xZbHbutE5NIks6E5K5bfNhqUS2TIuI3EY7ri1Vvdh7\nT2tq4VLAbwhWr7Y/cY8erbcFlcoNN9gDl25JdJ9EbrUlS8pTqXzuc+bagphSufXWzBqubt0sQ23e\nvPCVyrJllvXm986POcYa8aO8OU4PPdQSBfbf3xIQVqywygDxVRrGjbNg8/LlloIOZtn85CemDIcM\nseyyBx7ITL7vftde+YqpjRljlhiEH6SPp3t3+MMf7J6mYvx4K+Hy+c/nT550eecd6ygE24Js8N1f\n8aV+Ks1SeQUrGNkFOAAbm7IIGAeEPK47GviKw3fXBPF7Gi0t9mBMn575bJFr11o6av/+ra9Zrkrl\nuOPMHbRgga2vrbVBn5kydqwVAg1bqcS7OLt0sSylIGedZR0Mf3bQ0aPbjunZbz9zf73/vs2lArb/\n00/Dn/9syuSFF6z+Wbo8/bQp0ksuye67pUMwKB52kD4R556b3jUmTYpOGfxcBz361NSYdRs/HXpF\nxVRU9W6vmOR+QJ2q3qaqt2Fz1Ydwm6OH7zdPpFQ6drRGcf16ayRXrcpcqfiNWLCnUk5KZdCg1vNK\niJjrZ9ddcx+Y+uKL4Wd/vfdezPWVip49rUd5441tt3XrZjGWm26KFcMcM8biM/vsYw3pDTdYsD3V\njJBz55r1dMIJluoddowjiJ9W3NRkveju3fN3rUyYPNmSIaIQV8l10GOQ+LjK9u22nOl8KamIrFIJ\n0BsIGn+13rp2EZGhIvK0iLwlIvNF5OK8SRgS7Vkqwe3Tp1uacUNDZg9+vOsLTJGVi1KZOdMa1yD7\n7WcNVy6++jFjWsdUwsr+8gs0poOIKaG99kq8ffbsWCFMgB/8wGYJ9DnrLEtmuPfe9q/zj3/YDKUr\nV9ozlk+6d7cBfbNn23TJ6dbJyzd77GH3O340fjEIy1KBth2i5ctNoVRXh3N+n1JQKjcAc7xJuu7C\nKhWnE0fZAVyuqnsBhwAXicgeKY4pKpkolbPOMoWyYUP65493t/jnXLas9LO/IHlBylwKVULsnoXt\n/kpUxLM92vse8dvil32r7cor4fbbY7OVxjN9ug2s7NkzfblywbcKZszIvf5bWPiZeIVwgb30UvIs\nPdVwLZX4tOJ8uL6gBJSKqt4JHIxN1vUP4JDgHCvtHLdKVed6nzcDC4Ah+ZQ1V/zsr2RKxbcq/DIr\nw4Zl5gJLNNCpf3+Ls5SDpZIvginYYSqVQo8POuQQS53+97/hxz9uu/2jj2zEeXD6hnwzebKl/b70\nEnz2s4W7bioKUQb/jjtg4kT41rcSb1+xwhRcfE26bIlXKvkI0kMJKBWPauAjYB2wm4hMzORgEdkV\ni8O8GLpkIdKtm2V0LVyYuAfRv78FXQ880FwHw4dnplT8OlLx5/Sv7UiMn9od9jiVYgw6vegiG8yZ\naJbRJ56whj3s6Q3a46CDrCN18MHRspYHD85ftWJVm5b6pz+1BIrZsxPXffPHp4SVZh3vus2XpRLZ\nlGIfEfkZcCbwJuCHGRV4Ls3ja4GHgEs8iyWyiFgj/+abyd1fd99tDyNkbqn4I66D9Otn706pJKd7\ndyvvMWiQKZb168MpW9LYWJyGdPRoi8+0tLR2kz31lI2BKSTV1RYLio+FFZtBg9pXKsuWmQt62bK2\n23r0gPPPh69+tW3MrLnZarW98YYpkoED7f982WWW9h38PcIYSR8kPqaydKlZZGETeaUCfB7YXVW3\nZ3qgiHTAFMqfVfVfyfabGpgIpa6ujrq6usylDIn+/a2W05AEjrp+/WzAnO97zlSpBAfHBa8HTqmk\n4q23YplQ555rI81vvjm3cxarPE63bqYQly9v3VN9910bfV5o/vjH3Cb8ygftWSpz51p23GWXJZ56\nePlyG0s0apSlsF9yiblQN260KSc6d7aioP5vP2VKbH6Xs8+OnWfOHNs/LPr2tbbFZ9myzKYwqK+v\npz5YEjsZqhrpFzAdqM3y2HuAm1Lso1HimGNUBw9OvO3uu1WHDFFtabHlO+5Q/dKX0jvvjh2qHTqo\nNjW1Xt/QoAqqjz2WvcyVxocfqvbtq7pwYW7n+dznVKdPD0emTJk4UfWpp1qvGzlSddGi4sgTNXbu\nVO3YUXX79tbrZ85U7d9f9cEHU59j+XLVK65Q7ddP9cQTVfffX/XrX1dtbm677/PPqw4bptrYGFs3\napTqggW5fY8gd97Zur0YMya383ttZ5s2tRRiKluAuSLyOxG51X+lOkhEDgO+CBwpInNE5DURyWMR\n93Do3z954b7DD7fxBr6PNRNLZc0aM8X9cQw+zv2VOQMGWI/+73/P7TzFLOQZrLsF5gr74IPEJWwq\nkaoq+y8Gy7XcdZf97g8/nJ4FMWSIubaWLrWBuBdeCL/9bWKr7NBD7eVPA7Fhg107PlszF4IxlZYW\nazsqMqYCPOK9MkJVn8cC/CVFv35Wrj0Ro0a1HiyXiVJJ5PoCc+nU1jqlkil77mlZeLlQTKUSXyF4\n9WpLI87nYMdSw4+rDB1q7s677rLabHtkODCha1f4xjdS7/ezn9n05F/5irki99033DEkwZjK6tUW\nK8xHFYPIKxVNI324nBg61Hyu6TBsmPUuN2xIPq7g4YdNoYwZk1ip+OcJc66MSmDMGPjLX3I7R7EC\n9RCrEuCTLI29kvGVyi23wD//aYH1ZP+hMBgxwpTPlVeacglr0KNPMKU4H/Oo+ERWqYjIA6p6hojM\nJ0FhSVXdrwhi5Z2LL05dSsOnWzcLtF13Hfz854n3uesuq/dzzjnJ891feCH3gnWVRhgzBEbJ/eWU\nSlv8YP2TT1oKcD4Vis8PfmClht54I1bHLSyCSiVfY1QgwkoFG0V/IPA/QHOxhSkUmY4RuO46M5O/\n8Y22daS2b7fyF5062YRPyf4UTqFkztCh5krIRTEUU6nEpxX70wM4YviWSpij2lPRvbu52r761fAt\nld69Y3PS52uMCkR78GNf4JfAS8DdwNeBfYBNqrq0mIJFicGD4fTTbVrZeGbNsoKCHTvCyy8XpqdV\nKVRVmRJPNIgwHVQzL9MSJv60y/7Eb85SacugQZY+vH17Ye/NeefB1Km5z6EST4cO9rtv2JBf91dk\nlYqqfkdVDwUGAVcCa4EvA2+IyFtFFS5iDB2aeFbI6dNtTMu4cTZPt1Mq4ZKLC2z7dgvCxmfjFZJg\n4DZfmUClzKBB9r8ZPz78ycPao7ra3G3pxlYzwXeBLVtWmZaKTw1Wpbin91pBxMutFJrgNMQ33GAB\n2J07zXqZNMn+FP4Mgo7w8JXKt79tcalMKKbryydYcsZZKm0ZPNj+N2G7oYqJr1QqNVD/e2BvYBOm\nRGZjAxnXFVWwCBKcFfKBB+DBB23u8iFDrE7Y++/btrAK0zmMsWNtJPTSpTZC+uWX06+IXMzMLx+n\nVNrH74QVKp5SCPyxKvkM1EfZUhkOdAZWAR8Ay4H1RZUoogSVyqpV5qu/9FKbtEkk9qdwlkq4jBlj\n5VvuusuSIe65J/1jo2Sp7Nhh4xbCnqyp1PGnUS43S2XpUnO/5jJxXXtE1lJR1UkiIpi1cijwbWAf\nEVkL/FdVry6qgBHCVyo7d9r7I49Yz/mAA2z7mDGmZNKdEMqRHhMmwFVXwamnmtVx7bUWZE2HYgbp\nfXylsmqVuVCjVn+r2NTWwuWXZz7YMcr07WvZbMOH5y9OFOnHyKsv84aIrAc2eK8TgIMAp1Q8+vWz\nQP3HH5vimDChddXXqqrcix862tKrlykSsNLt8+a1rfybjChZKommRHAYiaZvLmX69rX5dPKZlBFZ\n95eIXCwifxORBqzM/QnAQuAUwI3/DtCnjwUUly93Lq5i0bu3Kfd0s8GioFT8oG2iKREc5UmfPjB/\nfv7iKRBtS2VXrGz9Zaq6ssiyRJqqKmvU3nzT9TiLybhxNq5ht91S7xulQP3Kle65qRT69oWmpvwq\nlchaKqp6uao+5BRKevTvbz0Q1+MsHuPHpz8FbRQslaD7yz03lYEfnK9I95cjM/r1M5++axyKh2+p\npINTKo5i4CuVirRUHJnRv78VoXONQ/HIRKlELfvLPTeVgbNUHGnTvz+sWOEah2IybJj5q1etsiyw\nRBWnFy2yTLwoWSouplI59O1rpYESTVceFk6plAn+XPOucSgeIjFrZfFiG23/3nut95kxw+Yvj4JS\n6dULNm1yWYOVRI8esGBBfmvOlbVSEZE7RORDEZlXbFnyjT8tsGsciosfrPcD9v67Py3tnDmmaFas\nKH72V1WVTe7W0BAbPe4of0aPzu/5y1qpAHcCxxZbiELgWypOqRQX31KZM8eUxpw5VuBz772t4oG/\nfvbs4lsqYC6w2triKzhH+VDWSkVVZwEVUYCyf38rle1KsRSX8eNNqcydC2eeae+PP26DDJ9/Ht5+\n28q6NDRER6m4jogjTMpaqVQSAwZYPKWQ8z442rL77hajeOklOP98s0ymT4eDDrICnyNHwmGH2b5R\nUSouDucIkyiPqC8YU6dO/eRzXV0ddXV1RZMlW/bdFx57rNhSODp0MFfXe+/BoYdaNdiFC206gmOP\nhSlTYlVvo6JUuncvthSOUqC+vp76+vqU+zmlQmulUqqIwF57FVsKB5jS6NEjNu1Ajx5QV2dxi3Hj\nbIrn6upoxDHyVf7cUX7Ed7inTZuWcL9KUCrivRyOgnDiifDBB/b51FNtnpJOneDCC+Hoo6GmBs49\n16aBLjaHHGIKzuEIC7Hq8uWJiPwVqAP6Ah8CV6vqnXH7aDnfA4fD4cgHIoKqtumwl7VSSQenVBwO\nhyNzkikVl/3lcDgcjtBwSsXhcDgcoeGUisPhcDhCwykVh8PhcISGUyoOh8PhCA2nVBwOh8MRGk6p\nOBwOhyM0nFJxOBwOR2g4peJwOByO0HBKxeFwOByh4ZSKw+FwOELDKRWHw+FwhIZTKg6Hw+EIDadU\nHA6HwxEaTqk4HA6HIzScUnE4HA5HaJS1UhGRSSLytoi8IyLfL7Y87VFfX19sET4hKrI4OdoSFVmc\nHF8Y3SgAABUySURBVG2JiizFlqNslYqIVAG/Ao4F9gbOEpE9iitVcor9IASJiixOjrZERRYnR1ui\nIkux5ShbpQIcBCxS1aWq2gz8DTi5yDI5HA5HWVPOSmUI0BBYXu6tczgcDkeeEFUttgx5QUROA45R\n1a95y2cDB6rqJXH7lecNcDgcjjyjqhK/rkMxBCkQy4HhgeWhwIr4nRLdFIfD4XBkRzm7v14GxojI\nCBHpBHwBeKTIMjkcDkdZU7aWiqruFJFvAk9gyvMOVV1QZLEcDoejrCnbmEoiREQ0Al/YydFWBhGp\nUtUWJ0c05IiaLEF5nBzRppzdXwCIyJdE5EaAYj4ITo6EsnwZuN+TpZgNuZMjorJE5XmNihwBeSaI\nSI9iy5GIslUqInKEiDwBXAsMFpHuTo7iy+HJMtGT5RKgVkRGOTmKL0eUZInK8xoVOQLyHC4iLwL/\nAZqKKUtSVLXsXsBFwDxsNP0YYC5Q7eQorhyeLFOA54ETgIHAU8AAb5s4OYojR5RkicrzGhU5PFl6\nAPd5v8/xWCLS8cV4TlLKWmwB8vQD9IpbfgY43clRXDm8a3dPIMtF3udCNlxOjojKEpXnNSpyeNfe\nEzgvsPxT4OvFkCXVqyzcXyJyjojcKiLV3qpN3vpqEakF5gAdnRyFlcO75vki8rCI+JmGW7z1/vX/\nDfTygsF581U7OaIrS1Se16jIEZBnvIj09hbfVtW7Apt7YtYkAXkjQckrFRG5FrgO2A34hre6BSyt\nWFU3A12xWmB+oUknR57l8M59CfBNYDBwtbdaPVmaveVaYLiqtuTxnjg5IipLVJ7XqMjhnXuYiLwM\nPAlMDCr1QAfgMeBMX758yZIVxTaVsjQFqwOf9wM+BdQBjwJDvPVVQJX3+TDMF9nNyZE/ObxzdySW\nqr43sAfmj54PjArI4u+zF/AO0NPJkT85oiRLVJ7XqMgRJ5MAhwLfAq4CfgKMSbDfnsCdwO75kiXb\nV8lZKl6P4q8ichGAqs5T1VeBt4AFWO/L2/RJKmQLVralj5MjP3J4slwP/BO4zsvlf1NV31bVxdgf\n9Vp/X/X+GcBGrFHby8mRHzmiJEtUnteoyBGQ5wsisod3718Gfg38CRgBfFpEunj7+WWltgC7e++R\noqSUioj8GNgHuBGYLCI3B27yx8A/gH1F5BBVG7TlbWsAlgIfOTnCl8OT5QKs8fkm9rDfIiK7Bna5\nFdhdRI7V1uMetgMvYFk2To6Q5YiSLFF5XqMihyfL3iLyBnAecLUnW4uqtqjqCuBpzELaB0zDee6w\npcBOYGJYsoRGsU2lTF7AzcBXvM8Dgeew9Lpqb11fzGy8BTPtjwI6OTnyK4d3rRuAa73P/TDT/AKg\nR2Cfc7Fe8e7AWeTH/ebkiKgsUXleoyKHd63jgJ95n/fzfptrAtu7AL8HvuYtB3+zXvmQKddXZC0V\nEakVkSvERrL686CsB3aISDdV/RD4K/YH6AygqmuAl4CzgVeAPqraFDhnxt/XyZFUlmtE5BIRGe+t\nfhPYJiL9VPVjLHNoArBr4NAZwJGeTNWq2hg4Z8bVop0c0ZUlKs9rVOQIHNtVrMitH3CfgGVygf1O\nNwOniMhIT5Zt2Ay2nxWRZ4BnRMSXc322cuSTSCoVEfkS9mMOBo4BrhKRgcAi4DN4Pk1VvR3rWdV5\nx+2NzfB4PzZo68HgeTXDchNOjoSynObJ0gMYBPxIRPYClmCToO3hnfthoBcWEEZEdsfcCncDg1X1\n3jhZMkpZdXJEV5aoPK9RkSMgz7ewAZQ3A/49vhP4vIgMUcsym4fFvi4OHDoROB1zQ35GVbdnc/2C\nUWxTKdEL+DlwlPd5LPA7YH+gGngY+DIw0Nv+PeAngWN3DXzu4OQITw7vHN8FjvQ+9wauB07DOii/\n8K6/m7f9QuBPgWN3CfGeODkiKktUntcIySGYu/EJYLS37lXgFO/zLcA9/rWwigY3YinMtdjI/v1z\nfT4K9Sq6AHE33/dp9gdqiKU2PgtM9j4f4z0c/w/LKX8GODbuPJ+kAWYpR5WTI6ks/byH3ZflQeAs\n7/OBWK7/P7DA4hPAF+LOIyHdEydHxGSJyvMaFTnizrU3Xskbb/nbeEod6wAsA07ylr8E/DKM6xbj\nVVT3l4h0FJEzRKSTiIh6g3hU9SNV3Wq7SDdgA5Z1gao+AfwS+wP8AXhOVWcGz6uWOZG2iSoiVb6f\nUgJlvoskR03w+GLI4cnSUUROSCDLx6q6xbsewDZglbftZax8xDysAXtJVf8WJ4tmcU86ep+liHKI\n7wcvphze9TuKyIGBcxTrnlSLyKAEchT8efXkqYqCHHEy+b/FIlVdHVgegSUIoKrrMCXzJRH5J/Aj\nLOuuNCmWNgPOAV7Hcr+He+va1BcCRmI/uL/8yb5Al8D6rGoTYdVYZwN3AGOTna8AcnwFmOXJkdTU\nzbcc3rFfxPL178KCiAnPhfmkX8bLjAH2DMjSMYR7ciFQj7kHdm3nt8m3HBdgqZ2/xXMfFeN+eMf6\ncYKrU+yX73vyv975HwE+T5JCiwX435wHXJ7Gfnn/3wR+ny8Ce3vLVXHbO3rvf8FzVQa29QL+B+ib\niwzFfhXcUhGRISLyb6wi6jlYZsMxkDQguDvwlojsJiIzgS94PUZR1W1eTzarSXNEZALmvzwH68H8\nSEQmJ5ElL3KISI2IPAScgTVeK7CBT8nI2/3w5BFgElZI8DxV3dDOuUZjQc8xYuXBvyI2dTOq2pzj\nb3M2cDKWgbMUuEJEjvbOHX++vMjhZQ7diwVJLwEagV+0k/2Tt/vhyfM08B3gDFWdlmL3fP425wCn\nYMrkcaysSbLz5Ot/01FEfou5rw4SkeO99cky1fL9v9lfRJ7DFMow4M8iMlTjLB2NlcAZDLwoIoeJ\nyA9EpFZV16vqP9Syz0qXQmkvYn7OUVgGg7/+N8CXNEkvAQs+tmA91i+GLNPpwFP+tYHLseDmHvHy\n5FmOcYHPB2CKNmFpjHzKEbj+f7zPQzFFtydQk+CenO7JMgvPdx+iHJcCNwWWZ2CW0+gkv2Nocvjf\nEYtPnBL87sDbwAHtPE+h3o+4+30FMNv7vAtWVqRPIWSJk+ObxALLe3m/y5Ak+4b6vNK6tMph3j34\nMjaAs0f89fMlR9y5O3n/k7MC6/5IgvIq3rYDMWvzXu/9zDDlKfarIJaKiFwH/EZEJqvqe6r6rMQq\na67GTL74Y/wex1YsaFWnqn/xtmU7vuJ8EdkjsGoesERE9lP7tWdi2RcHwyejV/Mlx57+sqrO9dZP\nwR6yflgJiX0Cx4QuR0CW+HuyU0R+BPwZK6D3M+AyvycXkKUTljFzuKrel4ssCeRYh40lGOctLwd2\nYI1IPGHKcQVwk4gcrxafmKmqW71efkdMqaxI0iMOTY44WU4CUNXrgQEiUo8F3C8FZojIAfmUJSCH\nH2ObjVlAfwXe8K71ZxE53ZMzOAo9zP/NVVg7co53nefVSqu8gVlKpyU4Ji//G+/4KSIyWlWbVPUP\ngft8EzAZOMz/D0tsXApANyw5YJ6qTlDV+3ORI3LkU2NhDfRd2OQyp2MZJ1f527z3fb31Q+OOTdTb\nyGqCHMyf+jJWduJioLe3fiwWvLwwsO/leKmFBIrr5UmO+PkaxmDjBACmYWmgvk88NDnakwV74L+J\nDf460Ft3FDaq9+AC/DZ9vPW7YQHLp7Cg5f/Diuv90P9tkpwvWzkOxEqb3wF8DbOM/FHMvuVS48nS\nL1/Pajuy+POaHIN1fgZ5yz/BUmeTpr2GeE9mEhuJPhqLC/hW/VeA2/BiXyE/IwdiKbh3Ye7qeuD8\nwPYuWEWA3+PFRQO/Wai/jXf8p7DOxWPePbgqsO1YrE35jPc8v53g+MF4bVA5vvJ34tiPejtwkPd5\nd8w/vltgvwmY0hmQzvmylGUU1jCeiKURTgxsOwe4CS+dEBgP/DfZ9fIlR4J99/bk6BK2HGnck4Ow\nQoKXectdsbz+YYW+J1jpigO8z0cAD+fpGTmZgBsCi/ndRmt3y5nAXd7nA/CCsXn4bZLJ4gd5qwLb\n9s3Xc9KOHH6H8A9Anfd5FyyZIWHpkBzlOAI4NLD8M2BC3D57A1MxN9S+JEmoyPW38c5xGjDV+zze\na79+4C3XxO07B5iU6zVL6RU0yUJFVVVEemE93w4iUq2qC0Xkfu/Hn+Lt94qIHIwNSnqyvfPlIM4y\noEEtQDkR+IyIvK+qDVjvqzPwMxHZAnwVS/Wrxtwt+ZZjsVrhuE8QEb/20PNY0bg25ChHMlneU9Xl\n2Kjf7wPXi8hr2J+oC9Cc6ER5uCfvquoH3rnnAYjIACwg/EiyE2UjRyA4+x/MjeOzC7BFVXd6z+5O\nrE7USs+9cSL2rIQiR5qyNHv7tHj798Wsyv8jwbOarSxpyLFDRPpg6crniMhHwNexYowJA+U5PiOz\nvPakK6ZQvurJeYSq3uyd/02xNOvrsQGDn8fK9ocph88E+KTtfB3LTrxDRO5UK/uCJ98IYCGmWCqG\n0GIq8T5m78FcD3yAmap+43gFVsdm/8Du92PpdKHLAaCqOzSWdfEXzLUyQUQ6qOpqVf0jFhz/IlZK\n+ipVTfgnzYMcB/nxJRHpJyLnYy6WdVjvJ2FDnidZ/HvSpKqPAz8GPos1WCer6qoCyXFg4J6IiByH\nWU7LiZW3CEUOv5FR1UZVXRcn4zpvm//sfgYbef0RsJeqPlsEWVSsftRZmItyLXBlmM9rmnKsxf63\nK7H/zk4sAL4uLDni5cHcRi9gJWh+D1wpIv6kWYd56+7BXKhP5SJHEtn89vIeYIqI9FIby/Ii5sL/\nlrffriLyXazcygIsblw55GrqQNtRuLSe2KY7Fvw9CujsrbuR1pkSnfMlR5J9vw9cg1Uj/XRgfcdU\nx+ZRjk956z4FjMhVjhDvSTCTJ1ufeLZy+G7TTnH3JNsxDZnIMRPP5YKXnYelWu+a6/3IUZb9vfd9\n4mQpxPPaRg7vc7Bybt6fkcD2XwLfCSyPDHwOo+xNxwTr/HbtduC3/n5YzPh6b/lUTNG2yVSshFeu\nNz3o290TC9YFBxL5vtevYAXavuTt918SDPDLobFoV44k+80AFmMBt4H+tb2HuxhyLAT6B/fJVo6Q\n7skAWiuUYv42/h8563uSrhze9j7Av4BxWG98JoEUb8w1mvffJoksMwgEeYt4T6Zjir+gv01gv70w\n68CPt/n/32qyVLDe8X2BbwWWd4n77f3vOwiz0vxaaxcSK2EfSmmXUn3lfgLztX8ZeBHLyvglsSyh\nYJDzBEy7vwp8O/Qv0r4crXrb2EC2LcA3y1WOKMlSgnKMwMY0zMfLuipXWUpQjt2xGMY84H/zIMcB\nWB2wk4Crve/7TyyG5tcm9DvL52EK9lGsA/T5eHkr8ZXpDY83T6uxDJB53nIN5rq4Guju7xPYvzPh\nuLqykcPvyVQBn/bXBx+SUpUjSrKUuBx+L3QcVh+ra/D4Av82octS4nL4z8hoLHZRGzg+12y7qsD5\nuwLnY8Unb/LWfROrcjElXn6sGOSpwWe20l8ZBeo1lnUy1gtS7QQeAHYVmw9gKzYwqhfwOe+YYPZS\nk6pul9jAx6zIUg71j1XVF1R1k4h08BIKsgpyRkWOKMlS4nL4hQjnquoPVXWL/6zGPcclKUuJy+E/\nI++q6m2qujkgh2YjhydDlffcqYj0URvw+iRmOfnFXf+KJQjsKyJdNVB6RVXXqerfVXVTtjKUGymV\niojcJDa6GrG6OQ9gWRZ/FpGDVPVJLE/7x94hz2EZX8eIyC7BcwUejIwfxjDlCMizI9MHMipyREmW\ncpXDU2pZNZxRkcXJkVCWQSLSBUzBichwEXkU+LXYHPEbsAGlo0VkkFqmWxM2MHlLogw1R4BUpgw2\n8GgNlvv9G+Cr3vp67IfvjI2KfgUvawgbAfvZVOfO5OXkiK4sTo7oyuLkaCVDNVapYjGxLLp+mCVy\nCmaZbMLG3Az35HzA2+9mTAlWdLwkrfuc4kfw/akPA7d7nz+FmaY3Y+U1vuOtvxp4Ni9COjkiK4uT\nI7qyODlayXAsNl7kOloXvxyMxW++gbm9/gR087Ydio0DehqLqZRtaZVQ73WKH8IPXvUFNgK7YkGr\nad76bwCbscyQrsDueRHSyRFZWZwc0ZXFydFKhoOBlsDyRExpDMespU/G4Hjb98LiKmfQzhw67tX2\n1W5MRdWqjarV978VeAgbXd1FREZ6D8GLmGbfolaGJXR/o5MjurI4OaIri5OjlQwvAv8QkQdF5GZs\nAHatqi7Dyrk8B6wSkS5evOdKYKeqPqCqbcq9OJKTMlCvseyPq7DR8QcC67EBjB+r6lGq+lZgf82H\noE6O6Mri5IiuLE6OVnwFOA4r+nig2pTCYEF5xYplzsYqC5+tIZRIqkjSMWeI+URPAd7xPgdH9eZU\nSjrdl5MjurI4OaIri5OjlQxTiU1C15HWAyuHUeJT+UbhldY4FbW0uyrV/9/e/bNIdYVxHP/+TLRZ\nJASCYCs2IWyKsLHLG7AQkdXKZItgYbANbJkuLyEhr0AQrOytBBvFf5UgkmyRLhE2LJGQfSzOXTLL\nbmQve2bnKt9PNXPmzuXH/OGZc5lznroN/JbkcrWN5j44zF/7xjLHdLOYY7pZzLErww+0BmOr1WYi\nH848tlHveivfCTjw4sfhA3GS1qf7xTD2bw0l/qiYY7pZzDHdLObYZZ22IzblJa7uxvZTWaH1D3g8\nhyzmOJypZDHHXlPJYg6gqm4mOZW2In97ET823mcZ83oOU9SFvwHm2GsqWcyx11SymENHYVRRkSTp\nbbp1fpQkyaIiSerGoiJJ6saiIknqxqIiSerGoiLNWZKPklwfbp8eNizscd6vkjxI8k+SSz3OKR2W\nRUWav4+B7wCq6vequtLpvL8Cawyrw6UpGLuiXtJ4PwJnkjykdR38tKqWk6wBF4El4CxtO/YTwNfA\n38D5qnqV5AytSdQnwBZwraqeV9u2nSQuNtNkOFOR5m8deFFVXwDf07ZZ3/EZrbCco3Ul/Gs47j7w\nzXDML8CNqvpyeP5PRxVcGsuZirRYd6tqC9hK8gq4M4w/BZaTLNE6FN6aaVx1fAE5pQOxqEiL9Xrm\nds3c36Z9P48Bfw6zF2nyvPwlzd8mrdshwKg2uVW1CbxMsrozluTzfQ6dS2tkaSyLijRnVfUHcC/J\nE/5rXbvvof8zfhX4NsmjJM+ACwBJVpJsAKvAz0medo4ujeYuxZKkbpypSJK6sahIkrqxqEiSurGo\nSJK6sahIkrqxqEiSurGoSJK6eQNdYztd5Z7jaQAAAABJRU5ErkJggg==\n",
"text/plain": "<matplotlib.figure.Figure at 0x7f348406fc10>"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true,
"collapsed": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "python2",
"display_name": "Python 2",
"language": "python"
},
"language_info": {
"mimetype": "text/x-python",
"nbconvert_exporter": "python",
"name": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12",
"file_extension": ".py",
"codemirror_mode": {
"version": 2,
"name": "ipython"
}
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment