Skip to content

Instantly share code, notes, and snippets.

@nathanieltarshish
Created April 9, 2017 23:19
Show Gist options
  • Select an option

  • Save nathanieltarshish/fe9b5f282c70432e6bf6af3cdf384daa to your computer and use it in GitHub Desktop.

Select an option

Save nathanieltarshish/fe9b5f282c70432e6bf6af3cdf384daa to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 60,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import xarray as xr\n",
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib\n",
"%matplotlib inline\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 53,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def dist_between_coordinates(lat1, lon1, lat2, lon2):\n",
" R = 6371 # Earth radius in km\n",
"\n",
" #conversion to radians\n",
" d_lat = np.radians(lat2-lat1)\n",
" d_lon = np.radians(lon2-lon1)\n",
"\n",
" r_lat1 = np.radians(lat1)\n",
" r_lat2 = np.radians(lat2)\n",
"\n",
" #haversine formula\n",
" a = np.sin(d_lat/2.) **2 + np.cos(r_lat1) * np.cos(r_lat2) * np.sin(d_lon/2.)**2\n",
"\n",
" haversine = 2 * R * np.arcsin(np.sqrt(a))\n",
"\n",
" return haversine"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"ds = xr.open_mfdataset('/work/net/global_xr/float_trajectories.*.nc')\n",
"ds = ds.drop(['u','v','z','vort'])"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (npart: 37166658, time: 2)\n",
"Coordinates:\n",
" * npart (npart) int32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 ...\n",
" * time (time) int64 0 1\n",
"Data variables:\n",
" y (time, npart) float32 -67.9844 -67.9844 -67.9844 -67.9531 ...\n",
" x (time, npart) float32 -279.977 -279.945 -279.914 -279.992 ..."
]
},
"execution_count": 47,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds"
]
},
{
"cell_type": "code",
"execution_count": 56,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"ds['delta_dist'] = dist_between_coordinates(ds.isel(time=0).y, ds.isel(time=0).x, ds.isel(time=1).y, ds.isel(time=1).x)"
]
},
{
"cell_type": "code",
"execution_count": 61,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f8bec1c4110>"
]
},
"execution_count": 61,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjYAAAF5CAYAAACMUGARAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XmcHHWd//HXGzQoaBIQCYegIIpxVSThlEPdyKV4si4M\nIAirIErEKC7qglyuIApBbpZDDmEUAigCEm5kucISEISEH4ccggkgIYRwBJLP74+qdmo63T3dNdXn\nvJ+Pxzy6p+rbVZ/qme761Le+hyICMzMzs16wTLsDMDMzMyuKExszMzPrGU5szMzMrGc4sTEzM7Oe\n4cTGzMzMeoYTGzMzM+sZTmzMzMysZzixMTMzs57hxMbMzMx6hhMbMzMz6xkdkdhI2lLSZZKekrRE\n0ufK1q8g6URJT0p6WdL9kvYpK7OcpJMkPSdpgaRpklYpK7OipPMlzZc0T9IZklYoK7OmpCskLZQ0\nR9LRkjrifTIzM7PaOuWEvQJwD/BNoNLkVVOBbYBdgA+kv58oaYdMmeOAzwA7AlsBqwMXl23nAmA8\nMCktuxVwWmllmsBcCbwJ2BTYA/gqcPhwDs7MzMxaQ502CaakJcAXIuKyzLL7gN9ExH9nlv0fcGVE\n/FjSaOBZYOeIuDRdvx4wC9g0ImZIGg/cD0yMiLvTMtsCVwDviog5krYHLgNWi4jn0jL7AEcB74yI\nN5r+BpiZmVlunVJjM5Rbgc9JWh1A0ieB9wHT0/UTSWpZriu9ICIeBJ4ANksXbQrMKyU1qWtJaog2\nyZS5r5TUpKYDY4B/KfKAzMzMrHjdkthMJql9+ZukRSS3i74VEbek61cFFkXEi2Wvm5uuK5V5Jrsy\nIhYDz5eVmVthG2TKmJmZWYd6U7sDqNO3SWpVdiCphdkKOFnS0xFxfTsDk/QOYFvgMeDVdsZiZmbW\nZd4CvAeYHhH/KGKDHZ/YSHoL8N8k7W7+mC7+i6QNgAOA64E5wChJo8tqbcal60gfy3tJLQusVFZm\no7IQxmXWVbItcH5DB2VmZmZZu5J08Bm2jk9sgDenP4vLli9m4FbaXcAbJL2dso2H1wJuS8vcBoyV\ntEGmnc0kQMAdmTI/krRypp3NNsB84IEq8T0G8Otf/5rx48fnOb6uMWXKFKZOndruMJrOx9lbRspx\nwsg5Vh9n75g1axa77bYbpOfSInREYpOOJbMuSZIBsI6k9YHnI+JJSTcBv5A0GXgc+ASwO/AdgIh4\nUdKZwLGS5gELgOOBWyJiRlpmtqTpwOmS9gVGAScA/RFRqo25miSBOU/SgcBqwBHAiRHxepXwXwUY\nP348EyZMKOgd6Uxjxozp+WMEH2evGSnHCSPnWH2cPamwphwdkdgAGwI3kPRQCuCYdPk5wF7ATsCR\nwK9Jbh09DvwwIv4ns40pJLU404DlgKuAb5XtZxfgRJLeUEvSsvuXVkbEknRsnFNIemItBM4GDinm\nMM3MzKyZOiKxiYibqNFDKyKeAf5jiG28RtJ7anKNMi8Auw2xnSdJGimbmZlZl+mW7t5mZmZmQ3Ji\nY3Xr6+trdwgt4ePsLSPlOGHkHKuP02rpuCkVuo2kCcBdd91110hq5GVmZjZsM2fOZOLEiZBMdzSz\niG26xsbMzMx6hhMbMzMz6xlObMzMzKxnOLExMzOznuHExszMzHqGExszMzPrGU5szMzMrGc4sTEz\nM7Oe4cTGzMzMeoYTGzMzM+sZTmzMzMysZzixMTMzs57hxMbMzMx6hhMbMzMz6xlObMzMzKxnOLEx\nMzOznuHExszMzHqGExszMzPrGU5szMzMrGc4sTEzM7Oe0RGJjaQtJV0m6SlJSyR9rkKZ8ZJ+L+kF\nSS9JukPSuzLrl5N0kqTnJC2QNE3SKmXbWFHS+ZLmS5on6QxJK5SVWVPSFZIWSpoj6WhJHfE+mZmZ\nWW2dcsJeAbgH+CYQ5SslvRe4GXgA2Ar4MHAE8Gqm2HHAZ4Ad0zKrAxeXbeoCYDwwKS27FXBaZj/L\nAFcCbwI2BfYAvgocPrzDMzMzs1Z4U7sDAIiIq4CrACSpQpGfAFdExA8zy/5aeiJpNLAXsHNE3JQu\n2xOYJWnjiJghaTywLTAxIu5Oy0wGrpB0QETMSdd/APhkRDwH3CfpYOAoSYdGxBsFH7qZmZkVqFNq\nbKpKE53PAA9JukrSXEm3S/p8pthEkiTtutKCiHgQeALYLF20KTCvlNSkriWpIdokU+a+NKkpmQ6M\nAf6lwMOyEea734Vbbml3FNbNXnsNXn+93VGYdb6OT2yAVYC3AQeS3CbaGrgUuETSlmmZVYFFEfFi\n2WvnputKZZ7JroyIxcDzZWXmVtgGmTJmDZs6FXbYod1RWDd7y1vggx9sdxRmna8jbkUNoZR8/S4i\njk+f3yvpY8A3SNremJn1vIcfbncEZp2vGxKb54A3gFlly2cBm6fP5wCjJI0uq7UZl64rlSnvJbUs\nsFJZmY3K9jMus66qKVOmMGbMmEHL+vr66Ovrq/UyK9gppyQ1I2uu2e5IzMwsq7+/n/7+/kHL5s+f\nX/h+FLFUJ6S2krQE+EJEXJZZdgvwcETskVl2CfByROyWNh5+lqTx8KXp+vVIkp9N08bDHwDuBzbM\nNB7ehuT21rsiYo6k7YA/AKuV2tlI2hv4GbBKRCx1h1vSBOCuu+66iwkTJhT/hlhDJNh0U7jttnZH\nMpgEY8fCvHntjsS6ValbRYd9ZZsNy8yZM5k4cSIkHXtmFrHNjqixSceSWRco9YhaR9L6wPMR8STw\nc+A3km4GbgC2B3YAPg4QES9KOhM4VtI8YAFwPHBLRMxIy8yWNB04XdK+wCjgBKA/7REFcDVJl/Lz\nJB0IrEbSrfzESkmNdaZXXx26jJk15qWXYNYs2Ki8Ttusw3RK4+ENgbuBu0h6KR0DzAQOA4iI35G0\np/lP4F6Srt1fiojsdfkU4HJgGnAj8DTJmDZZuwCzSXpDXQ78CdintDIilpAkTIuBW4FzgbOBQwo6\nTjMbIRYtgm99C154od2RFGOPPWDjjdsdhdnQOqLGJh17pmaSFRFnkyQZ1da/BkxOf6qVeQHYbYj9\nPEmS3Jj1rAg4/XTYbTdYfvl2R9ObbrgBTj4ZVloJjjii3dEM3+zZ7Y7ArD6dUmNjZi00axbssw8c\n4rrIppmVdnf45S/bG4fZSOPExmwEKg309tJL7Y2jl5U6eyxY0N44zEYaJzZmZmbWM5zYWCHuvRem\nT293FGZmNtJ1RONh637rr588eoyN6vzemJk1n2tszKytFiyARx5pdxRm1iuc2JhZW22zDay7bruj\nMLNe4cTGzNrq9tvbHYHZ0p59NhnrqZ2uuy6ZSmPhwvbG0W2c2JiZWdc57DA49tjmbX+vvWDvveG1\n15q3j6GcfXbyOKfmFMxWzo2Hzcys6xx6aPL43e82Z/se46l7ucbGzKwJeq0XXGl2cbNO58TGzMys\nIOPHw2c+0+4oRra6bkVJ+lyObV8TEa/keJ2ZmVlXmj3bE4a2W71tbH7X4HYDeB/waIOvMzMzM8ut\nkcbDq0bEM/UUlORp38xsRFtlleRxjTXaG4fZSFNvG5tzgEZuK/0aeLHxcMzMesP735887rFHe+Mw\nG2nqqrGJiD0b2WhE7JsvHDOz3uLeRGat5V5RZi3Sa91/zcw6UcMD9El6CzAZ+CSwCmXJUURMKCY0\nMzMzs8bkGXn4TGAbYBowg6QHlJmZmVnb5UlsdgA+HRG3FB2MmZmZ2XDkaWPzFODu3GZmNfRamyo3\ngrZukSex+R7wM0nvLjoYM7Ne44TArLXy3Ir6P+AtwKOSXgZez66MiJWKCMzMrBf0Ws2NWafLU2PT\nD6wB/Iikd9SUsp+GSdpS0mWSnpK0pNbcVJJOTct8u2z5cpJOkvScpAWSpklapazMipLOlzRf0jxJ\nZ0haoazMmpKukLRQ0hxJR0tyt3gza4hraszaI0+NzceAzSLizwXGsQJwD0mPq0uqFZL0RWATknY+\n5Y4Dtgd2JBn1+CTgYmDLTJkLgHHAJGAUcDZwGrBbuv1lgCuBp4FNgdWB84BFwEE5j83MzMxaJE9i\nMxt4a5FBRMRVwFUAUuXrHElrAL8EtiVJPrLrRgN7ATtHxE3psj2BWZI2jogZksanr50YEXenZSYD\nV0g6ICLmpOs/AHwyIp4D7pN0MHCUpEMj4o0ij9vMzMyKlecWyw+AYyR9QtI7JI3O/hQdIPwz2TkX\nODoiZlUoMpEkSbuutCAiHgSeADZLF20KzCslNalrScbh2SRT5r40qSmZDowB/qWAQzEzM7MmylNj\nc1X6eF3ZcpEkCcsOK6LKfgAsiogTq6xfNV1fPvHm3HRdqcyg2ckjYrGk58vKzK2wjdK6Im+/mZmZ\nWcHyJDafLDyKGiRNBL4NbNDK/TZqypQpjBkzZtCyvr4++vr62hSRdRr3jjGzkay/v5/+/v5By+bP\nn1/4fvIkNrdGxOuVVkhaeZjxVLIF8E7gyUzzm2WBYyV9JyLWAeYAoySNLqu1GZeuI30s7yW1LLBS\nWZmNyvY/LrOuqqlTpzJhgqfJMjMzq6TSxf7MmTOZOHFiofvJ08bmN5Ua+EoaB9w47IiWdi7wEWD9\nzM/TwNEkjX0B7gLeIOntVIpnPWAt4LZ00W3AWEnZmp9JJLfQ7siU+XBZgrYNMB94oLhDMrN2WLwY\nnnlm6HJm1r3yJDZrAWdkF0hajSSpmZ0nCEkrSFpf0kfTReukv68ZEfMi4oHsD8mggHMi4iGAtJbm\nTJJanE+kt6/OAm6JiBlpmdkkDYFPl7SRpM2BE4D+tEcUwNUkCcx5kj4iaVvgCODEarVUZt1spN0e\n++EPYdy4ocsVYaS9t2adIk9i82ngY5KOBZC0OklScx/w7znj2BC4m6TmJYBjgJnAYVXKV/rKmAJc\nTjLr+I0ktTo7lpXZhST5ujYt+ydgn39uNGIJySSfi4FbSWqLzgYOafiIzKzj3HBD6/fZKwP19cpx\ndCMnyY1puI1NRDwraRvgf9M7UjuQJCG7polBw9KxZ+pOstJ2NeXLXiMZCXlyjde9QDoYX40yT5Ic\nk5mZmXWZPI2HiYgnJW0N3AxcA3wlwjmldQb/J5qZjVx1JTaS5lH59s/ywGeBf5TaE3sSTDMzs+Hz\n7b986q2x+U5TozAzMzMrQF2JTUSc0+xAzMzMzIarrga7jc4BJent+cIxMzMzy6/enkjzJK0ydLF/\nekrSUj2XzMzMzJqp3jY2Ar4m6aU6y785ZzxmPcu9tczMmq/exOYJ4OsNbHcOyejAZmZmZi1Tb+Ph\n9zQ5DjMzM7NhyzOlgpmZDaHXbj16TBXrFk5szMzMrGc4sTEbwXqtVqGTuIbDrD2c2JiNQD7pmlmv\ncmJjZmZmPaPhxEbSY5J+LGmtZgRk1qtcS2Jm1nx5amyOA74EPCrpGkk7S1qu4LjMzMzMGtZwYhMR\nx0XER4GNgVnACcDfJZ0oaULRAZqZmZnVK3cbm4iYGRHfBlYHDgO+Btwp6R5Je0mueDezznL33e2O\noHu5B511i3qnVFiKpDcDXwT2BLYGbgfOBN4F/BT4FLBLATGa9QSfGNpv8eJ2R2BmzdZwYpPebtoT\n6AOWAOcCUyJidqbMpcCdRQVpZmZmVo88NTZ3AtcA+wK/i4hKk13+FfjNcAIzM+tmvVZD58YF1i3y\nJDbrRMTjtQpExEKSWh0zsxHNCYFZa+XpFVUzqTEzM7Pi9FrtX7PlGaBvWUkHSJohaY6k57M/eYKQ\ntKWkyyQ9JWmJpM9l1r1J0s8k3SvppbTMOZJWK9vGcpJOkvScpAWSpklapazMipLOlzRf0jxJZ0ha\noazMmpKukLQwPb6jJXmEZjMzsy6Q54R9CPBd4LfAGOBY4BKShsSH5oxjBeAe4JtAeW66PPBRki7l\nG5D0xFoP+H1ZueOAzwA7AluRdEO/uKzMBcB4YFJadivgtNLKNIG5kuQW3abAHsBXgcNzHlfHWLgQ\nTj7Zmb+ZmfW2PG1sdgW+HhFXSDoU6I+IRyTdS5IMHN/oBiPiKuAqgPLxbyLiRWDb7DJJ+wF3SHpX\nRPxN0mhgL2DniLgpLbMnMEvSxhExQ9L4dDsTI+LutMxk4ApJB0TEnHT9B4BPRsRzwH2SDgaOknRo\nRLzR6LF1ikMPhV/8ArbYAj7ykXZHY51ipCa6EW77Ytar8tTYrArclz5/iaTWBuByklqQVhhLUrPz\nQvr7RJIk7bpSgYh4EHgC2CxdtCkwr5TUpK5Nt7NJpsx9aVJTMp3kGP+l4GNoqYULk8c3ujY1syL5\npG5mvSpPYvM3oNS+5RFgm/T5RsBrRQRVSzov1VHABRHxUrp4VWBRWruTNTddVyrzTHZlRCwGni8r\nM7fCNsiUMTMzazpfgOSTJ7G5lKSNCiTzRB0h6SGSgfrOKiqwSiS9CbiIpJblm83cl5mZmXWfhtvY\nRMQPMs9/K6l0u+ehiPhDkcFlZZKaNYF/zdTWAMwBRkkaXVZrMy5dVypT3ktqWWClsjIble16XGZd\nVVOmTGHMmDGDlvX19dHX11frZWZmZiNCf38//f39g5bNnz+/8P3kniuqJCJuA24rIJaqMknNOiQN\ne+eVFbkLeIOkJunS9DXrAWtlYrsNGCtpg0w7m0mAgDsyZX4kaeVMO5ttgPnAA7VinDp1KhMmdO7k\n5qUqzZHaWLQT+L03s5Gs0sX+zJkzmThxYqH7yTNX1L8CXwLeQ3JL6K/AtIj4U94g0rFk1iVJMgDW\nkbQ+SfuXv5N02/4osAPwZkmlWpTnI+L1iHhR0pnAsZLmAQtIemfdEhEzACJitqTpwOmS9gVGkdxK\n6097RAFcTZLAnCfpQJK2REcAJ1aZOsLMrKJeS2Td3sO6RUNtbCSdStKTqA94B/BOku7fN0g6YRhx\nbAjcTVLzEsAxwEySsWvWAD5LMmv4PcDTJMnO0wz0eAKYQtIzaxpwY7p+x7L97ALMTo/hcuBPwD6l\nlRGxhCR5WgzcStJu6GySsXvMekYnnXR7/YTZ68dn1mnqrrGR9EWS+Z/2As6JSL4a00HtvgqcIuma\niLis0SDSsWdqJVlDJmAR8RowOf2pVuYFYLchtvMkSXJj1vNG6km3lePYdFISaTYSNFJjsydwbESc\nXUpqIKnliIizSEb+/Y+iA7RijNQTmFm7+DNn1h6NJDYTSBvmVnEJyUB51sF89WhmZr2skcRmZZLB\n+ar5G0m7GzMzM7O2aCSxGQXU6hn0RlrGzMzMrC0a7e59hKSXq6xbfrjBmJmZmQ1HI4nNn4D16ihj\nHWgkNWR0O6L6+b0yq82fke5Td2ITEZ9oYhzWIv6QGoysRNfMRpY8k2CamRWmHUmWE3yz3uXExqxF\nfDI1M2s+JzYjhG89mLVWryWy/g6xbuHExsysiZwQmLWWE5sRpteuIs3Mep2/txtTWGIjaQVJWxW1\nPTMzM7NGFVljsy5wQ4HbMzMr1Lvf3e4IzKzZfCtqhPB9frMBrto36111D9An6fkhiiw7zFjMzJrK\nCb5Z72tkSoXlgFOA+6qsfzdwyLAjsqbylaqZmfWyRhKbe4AnI+KcSislrY8TGzNrkGtRzKxIjbSx\nuQIYW2P988C5wwvHmsUnDzMzGwkamQTzp0OsfxLYc9gRmVnL+NakNSrCF0qt4vc5H/eKMmuRTkoi\n/IXZfJ309y6C/2esWzixGWF67cvWrBH+/zfrfU5sRghfbZkNaEWC48+cWXs4sTGzEcPJhlnv64jE\nRtKWki6T9JSkJZI+V6HM4ZKelvSypGskrVu2fjlJJ0l6TtICSdMkrVJWZkVJ50uaL2mepDMkrVBW\nZk1JV0haKGmOpKMltfx9euEFeOWVVu/VzMysu+U+YUsaJeldktbK/uTc3Aok4+R8E1iqkljSgcB+\nwN7AxsBCYLqkUZlixwGfAXYEtgJWBy4u29QFwHhgUlp2K+C0zH6WAa4k6S22KbAH8FXg8JzHlduK\nK8Imm7R6r2ZmBp3RHqsTYuhGjQzQB4Ck9wFnAR8rX0WSlDQ8tUJEXAVclW6/UmXx/sAREXF5WmZ3\nYC7wBeBCSaOBvYCdI+KmtMyewCxJG0fEDEnjgW2BiRFxd1pmMnCFpAMiYk66/gPAJyPiOeA+SQcD\nR0k6NCLeaPTYhuO+amM8m/UQ3x4ysyLlqbE5G1gC7ABMBCakPxukj4WStDawKnBdaVlEvAjcAWyW\nLtqQJEnLlnkQeCJTZlNgXimpSV1LkoxtkilzX5rUlEwHxgD/UtAhmZmZWZM0XGMDfJSk1mN20cFU\nsSpJ8jG3bPncdB3AOGBRmvBUK7Mq8Ex2ZUQsTif3zJaptJ/Suj/nOYBO4qpNM7OhuSaxe+VJbB4A\nVi46kG43ZcoUxowZM2hZX18ffX19bYpoMH9IzZzYm7VTf38//f39g5bNnz+/8P3kSWwOBI6W9COS\nmb5fz66sUGsyXHNI2u+MY3Btyjjg7kyZUZJGl+1/XLquVKa8l9SywEplZTYq2/+4zLqqpk6dyoQJ\nhd+JM7MmcIKTn6dUsLwqXezPnDmTiRMnFrqfPG1sriVpi3Idya2deenPC+ljoSLiryRJxaTSsrSx\n8CbAremiu4A3ysqsB6wF3JYuug0YK2mDzOYnkSRNd2TKfFhStkZqG2A+SU2VmXWxVp6Qey15cjJj\n3SJPjc0niw4iHUtmXZIkA2AdSesDz6eTax4HHCTpYeAx4Ajgb8DvIaklknQmcKykecAC4HjgloiY\nkZaZLWk6cLqkfYFRwAlAf9ojCuBqkgTmvLSL+Wrpvk6MiEE1U2aN6rUTndXHCYFZazWc2JS6Uxds\nQ+AGkkbCARyTLj8H2Csijpa0PMmYM2OBm4HtI2JRZhtTgMXANGA5ku7j3yrbzy7AiSS1TkvSsvuX\nVkbEEkk7AKeQ1AYtJOkFdkhRB2rWSZxsmVmvqSuxkfQR4C/pif8jtcpGxL2NBpEmSzVvi0XEocCh\nNda/BkxOf6qVeQHYbYj9PEnSld2sZ7kWwcx6Vb01Nvcw0F36HpJalUpfjbkG6LPW8RV6+ziZqMzv\ni1lt/t5uTL2JzdrAs5nn1mV88jAzs5GgrsQmIh6v9NzMrJv4ytes93XE7N5mZq3kBKdxfs+sWzix\nMbMRw7dkzXqfExszMzPrGU5sRhhXJ5uZWS9zYjNCuArerLV69SKiV4/Leke9A/TNIxmjZkgRsdKw\nIjIzs47jiyPrFvWOY/OdzPN3AAcB0xmYYHIzYFuSeZXM2qpTryg7Na5269UTZq8el1mnq3ccm3NK\nzyVdDPw4Ik7MFDle0n7Ap4CpxYZoZlYsJ5nWDZwc55Onjc22JBNMlruKJLExMzMza4s8ic0/gM9X\nWP75dJ2ZdQnXXJhZr6m3jU3WIcAZkj4B3JEu2wTYDvh6QXFZk/hEZuAqbjPrXQ0nNhFxtqRZwLeB\nL6WLZwFbRMQd1V9p7eQTmZmZjQR5amxIE5hdC47FzMzMbFhyDdAn6b2SfiLpAkmrpMu2l/QvxYZn\nZmZmVr+GExtJHwfuI2lXsyPwtnTV+sBhxYVmZmZm1pg8NTZHAQdFxNbAoszy64FNC4nKzKyJ3Ig+\nP7931unyJDYfBi6tsPwZYOXhhWNmI02vNmzvtQSgV/9O1nvyJDYvAKtVWL4B8NTwwrFm67UvW7NO\n54TArLXyJDa/AX4maVWSiTGXkbQ58Avg3CKDs+L4y7X9nFSamTVfnsTmR8Bs4EmShsMPAH8CbgV+\nUlxoZmZmZo1pOLGJiEUR8XXgvcAOwG7AByLiKxGxuOgAASQtI+kISY9KelnSw5IOqlDucElPp2Wu\nkbRu2frlJJ0k6TlJCyRNK3VXz5RZUdL5kuZLmifpDEkrNOO4zMzMrFh5untvARART0TElRFxYUQ8\nVHxog/wA2Af4JvAB4D+B/0xnFC/FdSCwH7A3sDGwEJguaVRmO8cBnyHppr4VsDpwcdm+LgDGA5PS\nslsBpxV/SGZmZkPzbezG5Bl5+HpJTwH9wK8j4oGCY6pkM+D3EVGaVfwJSbuQJDAl+wNHRMTlAJJ2\nB+YCXwAulDQa2AvYOSJuSsvsCcyStHFEzJA0nmT28okRcXdaZjJwhaQDImJO8w/VzMzM8srTxmZ1\n4Bjg48BfJN0j6fuS3lVsaIPcCkyS9D4ASesDmwNXpr+vDawKXFd6QUS8SDJJ52bpog1JErlsmQeB\nJzJlNgXmlZKa1LUkjaQ3Kfyo2sCZv2WN1P+HkXrcZiNBnjY2z0XEiRGxOUk7m4uAPYDHJF1fdICp\no4DfArMlLQLuAo6LiN+k60s9tOaWvW5uug5gHLAoTXiqlVmVZDyef0rbDT2fKdOV3CvKsvz/YGa9\nKtckmCUR8VdJRwF/Bo4gqcVphp2AXYCdSXphfRT4paSnI+K8Ju3TzFrASZaZFSl3YpOOXbMr8G/A\nW4DfAz8sKK5yRwNHRsRF6e/3S3pPur/zgDmASGplsrU244DSbaU5wChJo8tqbcal60plyntJLQus\nlClT0ZQpUxgzZsygZX19ffT19dVxeGZm3cG38Syv/v5++vv7By2bP39+4ftpOLGRdCRJzcnqwDUk\njXZ/HxEvFxxb1vJAeVfyJaS30tKaozkkPZnuTeMcTdIu5qS0/F3AG2mZS9My6wFrAbelZW4Dxkra\nINPOZhJJ0nRHrQCnTp3KhAkT8h6fWUv55NR8vfYeu2bNhqvSxf7MmTOZOHFiofvJU2OzFfBz4MKI\neK7QaKr7A3CQpL8B9wMTgCnAGZkyx6VlHgYeI7k19jeSmiQi4kVJZwLHSpoHLACOB26JiBlpmdmS\npgOnS9oXGAWcAPS7R5T1Ip+szKzXNJzYpI2GW20/kkTlJJJbRU8Dp6TLSnEdLWl5kjFnxgI3A9tH\nRHYG8ikkNT/TgOWAq4Bvle1rF+BEkt5QS9Ky+xd/SGbWy5w0mrVHrjY2kt4LfIdkIDtIGvT+MiIe\nKSqwrIhYCHw3/alV7lDg0BrrXwMmpz/VyrxAMppyT+q16vFu4vfezKz58ow8vC1JIrMxSXuWe0na\nstwvaetiw7Oi+OrRbICTTKtXJ3x3dkIM3SRPjc1RwNSI+EF2Ydrt+2ckDYrNzMy6npPg7pNn5OHx\nwJkVlp8MKOg6AAAgAElEQVQFfHB44ZiZmZnllyexeZZkgLxyH6Vs1F4zs6G4mt2sMn828slzK+p0\n4H8krUMyhxMk8zYdCBxbVGBmZmZmjcqT2BxBMgbM94Aj02VPk/RGOr6YsMzMzMwal2ccmwCmAlMl\nvT1dtqDowKw53BDOsvz/YI3y/4x1uuFOgumEpkv4Xq1ljfT/h1acnHstARjp/zPWPepObCRdX0+5\niPjX/OGYmTVPO07OTgjMWquRGptPAI8DVwCvNyUaM7Mm6rVaFDNbWiOJzYHAnsCXgfOBsyLiL02J\nyszMzCyHusexiYifR8QHgS8AbwdukTRD0jckjW5ahGY9wrUFZmbN1/AAfRFxW0R8HViNZLbtvYCn\nndx0B59crdO4DYqZFSnPyMMlE4CPk0yx8Bfc7qaj+eRhZtadfEHamIYSG0mrS/qRpP8HTAOeBzaJ\niE0j4pWmRGjWIH8JmJmNXI10974S+CRwNfB94IqIeKNZgZmZNYuTX7Pe1UivqO2AvwNrAYcAh6jC\n/Y2ImFBMaGZmxfItWbPe10hic1jTojAzawHX1Jj1vroTm4hwYmNmNsI5ObRON5xeUdaF/KVkZnn4\nNp51Cyc2I4S/lKwSJ7pm1muc2JiNQE50zaxXObExaxEnE5X5fbFu98c/wt13tzsKK2mkV9RSJL0l\nIl4tKhizXubbPp3Df4vG+T2r7tOfTh79HnWGhmtsJC0j6WBJTwEvSVonXX6EpP8oPMKB/a4u6TxJ\nz0l6WdKfJU0oK3O4pKfT9ddIWrds/XKSTkq3sUDSNEmrlJVZUdL5kuZLmifpDEkrNOu4zKx1XDtk\n1vvy3Io6CPgq8J/AoszyvwBfKyCmpUgaC9wCvAZsSzI/1feAeZkyBwL7AXsDGwMLgemSRmU2dRzw\nGWBHYCtgdeDist1dkG5/Ulp2K+C0wg+qTXxFYSOZ///Nel+eW1G7A3tHxHWSTs0s/zPwgWLCWsoP\ngCciIps4PV5WZn/giIi4HEDS7sBc4AvAhens43sBO0fETWmZPYFZkjaOiBmSxpMkThMj4u60zGTg\nCkkHRMScJh1f0/lK1czMRoI8NTZrAA9X2dabhxdOVZ8F/k/ShZLmSpop6Z9JjqS1gVWB60rLIuJF\n4A5gs3TRhiSJXLbMg8ATmTKbAvNKSU3qWiCATQo/KjMzMytUnsTmAWDLCsv/DWhWu/B1gH2BB4Ft\ngFOA4yV9JV2/KknyMbfsdXPTdQDjgEVpwlOtzKrAM9mVEbGYZBbzVTEzM7OOludW1OHAOZLWIEmM\nviRpPZJbVDsUGVzGMsCMiDg4/f3Pkj4EfAM4r0n7bMiUKVMYM2bMoGV9fX309fW1KSIzs+K5nVLr\n9FoTgv7+fvr7+wctmz9/fuH7aTixiYjfS/os8GOSBrqHAzOBz0bENQXHV/J3YFbZslnAl9LncwCR\n1Mpka23GMVCLNAcYJWl0Wa3NuHRdqUx5L6llgZUyZSqaOnUqEyZ4YnMzS/RaAtBrJ1lrvUoX+zNn\nzmTixImF7ifXAH0RcXNEbB0Rq0TE8hGxRURcXWhkg90CrFe2bD3SBsQR8VeSxGNSaWXaWHgT4NZ0\n0V3AG2Vl1gPWAm5LF90GjJW0QWY/k0iSpjsKOhYzy2jHCbOVSYcTArPWarjGRtJGwDIRcUfZ8k2A\nxRHxf0UFlzEVuEXSD4ELSRKWrwFfz5Q5DjhI0sPAY8ARwN+A30PSmFjSmcCxkuYBC4DjgVsiYkZa\nZrak6cDpkvYFRgEnAP3d3CMqq9euIs0a4STDrPflqbE5iWT8l3JrpOsKlyZLXwT6gPuA/wL2j4jf\nZMocTZKEnEZSu/JWYPuIyI61MwW4HJgG3Ag8TTKmTdYuwGyS3lCXA38C9in8oFrMX+hmTuzNRoI8\njYc/CNxTYfnd6bqmiIgrgSuHKHMocGiN9a8Bk9OfamVeAHbLFaRZl/GJ3sx6TZ4am9eo3PV5NZI2\nLGZWQSclEa7BM7NelSexuRo4UtI/+zanUx78FGhWrygzMzOzIeW5FXUASbuTxyWVulJ/lKSb9Veq\nvsrMzMwa1km1vd0gzzg2T0n6CLArsD7wCvArkp5Drxccn5lZ4XyiMOtdeWpsiIiFwP8UHIu1gL/Q\nbSRz2yKz3pcrsZH0PuCTJKP0DmqnExGHFxCXFcxf6Nap/L/ZXXxxZJ0uzwB9XyeZhPI5ktF+s//m\nQTLFgplZx2nlSbnXEgAnoNYt8tTYHAT8V0T8rOhgzMzMOkEnJKadEEM3ytPde0XgoqIDMTPrJa7h\nMGuPPInNRcA2RQdiZmZmNlx5bkU9DBwhaVOSeZsGdfGOiOOLCMyaw1WbZmbWy/IkNnsDLwEfT3+y\ngmTGbOswrhY3G+AE34bi78zulWeAvrWbEYhZr/PJtP18sjLrfXna2JhZj3CyZZb4/OfhkkvaHYUV\nIe8Afe8CPgesBYzKrouI7xYQl5k10UituXAiZ9Vcdhlccw28/HK7I7HhyjNA3yTgMuBR4APAX4D3\nAAJmFhmcmfW+kZpkmVlz5LkVdSTwi4j4MPAqsCOwJnATHt/GzKxujzwCl1/e7ijMekueW1Hjgb70\n+RvAWyPiJUk/Bn5PMt2CdShXxZt1jo98JLn10U2fy26K1UamPDU2CxloV/N34L2ZdSsPOyJrClf3\nm7VWPQlAN7Xn8HeIdYs8NTa3A1sAs4ArgWMkfRj4UrrOrK18RWlDaeX/iBMCs9bKk9h8F3hb+vyQ\n9PlOwEPpOjOzjuQkw7qJ/1/zyTNA36OZ5wuBbxQakZmZmVlODbexkfSopHdUWD5W0qOVXmNm1gl8\nm9Ks9+VpPPweYNkKy5cD1hhWNHWS9ANJSyQdW7b8cElPS3pZ0jWS1i1bv5ykkyQ9J2mBpGmSVikr\ns6Kk8yXNlzRP0hmSVmjFcZmZmdnw1H0rStLnMr9uK2l+5vdlgUnAYwXFVSuOjUgm4vxz2fIDgf2A\n3dM4fgJMlzQ+IhalxY4DticZe+dF4CTgYmDLzKYuAMaRHM8o4GzgNGC3phxQi/mKtX383lfmdgRm\nVqRG2tj8Ln0M4Jyyda+TJBPfKyCmqiS9Dfg18DXg4LLV+wNHRMTladndgbnAF4ALJY0G9gJ2joib\n0jJ7ArMkbRwRMySNB7YFJkbE3WmZycAVkg6IiDnNPL5m8snDzMxGgrpvRUXEMhGxDPAEsErp9/Rn\nuYhYr5RUNNFJwB8i4vrsQklrA6sC12XifRG4A9gsXbQhSSKXLfMgyfGUymwKzCslNalrSZK5TQo9\nEjMzszq4trcxDbexiYi1I+K57DJJY4sLqTJJOwMfBX5YYfWqJMnH3LLlc9N1kNxeWpQmPNXKrAo8\nk10ZEYuB5zNlzP7p6afbHcHwjNQvzGYe94UXwuzZzdt+0ZYsgddfb3cUVsmvftXuCLpTnl5RB0ra\nKfP7RcDzkp6StH6h0Q3s410k7WN2jQh/BHvUM8/AggXV10vw3//duniGcs89sMYacOWV7Y7E6tWK\nW7I77QTjxzd/P0X5yldg1Kihy5WM1GTYukeeAfq+AewKIGlr4FPAdsC/Az8HtiksugETgXcCM6V/\nfjUtC2wlaT+SWcZFUiuTrbUZB5RuK80BRkkaXVZrMy5dVypT3ktqWWClTJmKpkyZwpgxYwYt6+vr\no6+vr8orrNy4cbDWWvD449XL/PrX8F//1bqYank0Hdzg/vvh059u3X6ffx4+8Qm46ipYffXhbauI\nE30EnHce7LorLFupv2QLLV6cPDY7jjfegLPOgq99DZbJ07e0g1xwQX3likwKr78+qSXadtvitlmE\nXk3aLrkEPvUpGD26vXH09/fT398/aNn8+fOrlM4vT2KzKvBk+nwH4MKIuFrSYyRtWprhWuDDZcvO\nJpnW4aiIeFTSHJKeTPcCpI2FNyFplwNwF8mknZOAS9My6wFrAbelZW4DxkraINPOZhJJ0lTz2KZO\nncqECRMaPrAFC2CzzeCyy2CddRp+ecM6/YP7xBPtjqB+O+6YPP7nf8L3v9+6/V5zDdx3H/z2tzBl\nSuv2W81ll8EeeyRJxZ57tjeWlVaCF1+s/n9e1P//WWfBPvvAaqvBZz9bzDZHkkmTksdO+D6aOxdW\nWWXocp3s/PPhbW+Dz39+6XXPPDPwXVX+fkvwgx/AkUc2P0aofLE/c+ZMJk6cWOh+8lxrzAPWTJ9v\nR5J0QHLyb8p1UkQsjIgHsj8kk3H+IyJmpcWOAw6S9Nl07qpzgb+RzDheakx8JnCspE9ImgicBdwS\nETPSMrOB6cDpkjaStDlwAtDfrB5Rd96ZXPGfcUYx27v99uSf9dGyoRIbudr6+9/h1luLiafcnXfC\nsccOXa5VvvGN7usx9o9/JI/z5lVef/vtSc1XkSeNnXeGAw6ovO6hh5LHe+4pbn95vVjegq5JFi5M\nHl99tTX7s/ymT08uPiqZPx9WXRW+9a3WxlSPl16qv+xuu8EXvlB5XbXviZKjjqp/P90iT2JzCXCB\npGuAdwB/TJdvADxcVGB1GPS1HRFHkyQhp5HUrrwV2D4zhg3AFOByYBpwI/A0yZg2WbsAs0kStsuB\nPwH7FB/+8F14IcyaNXjZDTckj3/5S/7trr46bL55/tfXsvHG8L2mDgrQmNNOa3cEjbvoouTxd7+r\nvH6zzeDJJwdOvkX47W/hmGMqr/vTn5LHP/yh8vo//hHe/OaB20Tt8te/wmOP1V/+b3+DNdes3e7L\nOt9228HPf1553WGHJY+nnFL99c2oVZo2LbmgeuqpwctffRW+8x145RVYtKjyaxv1xhvFbKeb5Els\npgAnAg8AW0dEKa9cDTi5qMCGEhH/GhHfLVt2aESsHhHLR8S2EfFw2frXImJyRKwcEW+PiC9HRHkv\nqBciYreIGBMRK0bE1yPi5WYfT54v/Z12gg+X3aArfQhLV9HD8eCDw99G0Vpdu3L88Z1Xo3Pjjcnj\nfffVLvfaa00PBRhIaKpdYR57bPLlWq12o5H395ln4EMfgpdzfCIbvSpfc80kufnUp6qX6YRbKdU8\n/3wSv1V3csvOWIN9+cvJ46WXDl7+m9/AL3+ZtFlbsqSYfX3oQwPPi9pmp8vT3fv1iPhFROyfHe8l\nIqZGREE3VEaOO+9MHk88Md/ryxOi3/42eax226CWvr7B3QsrfSkuXAgf//jwPyB5Twit/GAuWQL7\n7588v+WW1u0XYJdd4KCDhreNVn+JVavyLvLkP3lycuv2a18b3nYaiWnGjKWXdVqyW8n7358kZ+1w\nwgnJe3TbbUOXbadWJf/V1KpNee656uuqGaodbifcLm6FuhoPp9Mp/DEiXi+bWmEpEXFZIZGNEKVq\n7jxXoJXce2/+1/7mN8lPLautlsS85ppLV6M24ne/gy9+sf7y06cnj62sRcp+6W2xxeCT4VVXNXff\npY4DP/lJ/m00I7FZvLh6j6NW1F6Uaqkuuqj+3jzN0Mk1NSWltliQ1Ja99BKsvHJr9v3tbyePH/tY\nd7xX7VL+ffbKK8nj88/n295TT0FZ59xBJk4cGX+PemtsfgesmHle7efSiq+2qtrd7qCW8vY7MJCI\nDXdgukbb2Wy33fD2l0fpS6aS7bdvfHut/kIpKlnOelONS6GhEqkiajlK/5OVrnSzn6WR8OXdiO23\nh3e+s91RJNxmaUD57dtvfjN5/GGlYWjrUHDnoq5VV2KTTpvwTOZ5tZ82j2LRfW66qTX7Kd1KaeQL\nf/Lk5sQCSUPOrFpJRLsMp0aqE7S60WCzu1gPZcdMN4DyhtO33ZY0Yi5KtyVOpXZZnWCzzaqve+GF\npOt1J7bvq6aRwQ3L1fo/ytPjrvw1Z53V+DZ6QZcPLdX9WnEP+uGH4eqrm7+f4fjxj9sdwdK6fXyS\nTusN0ex2Kb///cDz8m7fO+3EsJTXRv3yl8njQQcN9Kyx+tx/f/V1d90Fzz6bDMTZLYYzHUWtxt15\nhy7IJkv/8R/5ttHtGkpsJC0jaS9Jl0v6i6T7JF0maffMiMBWZtq0gS66zTBU75iTTqq9vl7NvG32\ni180b9t51RoBuVNle8kV1V10uNpRu7HGGoN/f/LJyuXqVX4rtDSQ5EMPwaGHDizvppoGyPe3afQ1\njZRvZTL+cCsHJ6ni2WeL3+bPflZ7fTcNgppX3YlNmrhcBpwBrAHcB9wPvJtkFGC3r6niy1+Gf//3\n5m3/Ix+pvf6444rZT632FXnkbSDXDPUkbaVEYW75VKsdZDjjFzVLaWyloi99WjnQ4zXX1FfuAx8Y\neN4Nt6uGuigqNdqH/H+///3f+suWEshW3EJ53/uav4+h/L//V/w2h2qf87GPFb/PTtNIjc1Xga2A\nSRGxQUT0RcTOEbE+yXxR/ypp92YE2St6dZTSagPFDeUd7yhm/6+8MvzB6OpJ2h54IHn86U+Ht69W\nKT+xvv56cnL6n/+pXqYeeSciLXLAQGjfQI+NnuAbKX/xxbD33o1tfzhOPbX6uoMPThKNQw4Z3j62\n2qr6umq3cYbbOaFbtON2cbe3HaxHI4lNH/DTiLihfEVEXA8cRTo5plXWyT2ghqORbtv1OPhguO66\n+suvvXYyT0qrHH98vte1+mZtedJSSqwvumh4sWTH1ylvu3X++dVfV63NQDtuYrfiFkwe//ZvcPrp\nxW6zVg1jrRF3S7U1pb9xEcdfXkvbyjmCzz67dfuy9moksfkIUGv0jj8C6w8vnN525ZXtjqA7/OQn\ng0d7HeqecKtuDW2wQWv2k/XCC8Vta04629m119Yu14jLLx/8+267Df49W0vZ7aOeXnJJ7fWdOhhd\npYkR61EaPPT224uLZcstB/9+8cUDz2fOLG4/lbR7glZrnUYSm5WAWqeQuQyMdWMVNDKpWbM0s+Ff\ns7z73a3fZ55RP5uhyCr50pxORTrhhNrrm30yaeUwAf/2b3DHHdXXd2pj82oxt+P7qHQ7t9yzz7Z2\nDJZqY+l0wvddyXAaFr/3vXDFFdXXR8D6PVwN0UhisyxQ647gYuocyXik2muvdkfQuGXS/5BmDPY2\nlHb26qk2mFn5+DvVfP/7A9NbDEcjQ6CXf1l3wthA2ZGsi25jA80da6mSTTetvHzmzM5suF3LM88M\nXaZV6kkKN9+89sm6EUPVvnWC4bTJfPRR2GGH6uuvvHJ4o9R3ukYSGwFnS7qk0g8wQocCGhkaafNS\nlJ12Kq79xSuvJNuqNUjbUL3LANZZp779/eIXsPPO9ZWtZdeyVmu1xrYYPXrw753W9bgZDSVnzy5+\nm3lMnNh4o+rsuDu1PPlk7XZdS5YMnj6hXuVJet7E7NhjB25z5lXP7a5bb619sm7EV79aeXk7amyq\nXSxlk6+i57Qq6n3sVI0kNucAzwDzq/w8A5xbdIDdZP78oe8TZ09M9Q6a9/Of13eCb2bj5M/VnCGs\nOfL2tqqkdGupVgPCobq+NsOCBfXXAkHteWDKNbOW7dVXh05Uyr+My69AX3xxYE6sSi66aOhanj//\nufb6TvaFL1Rfd+utA8932mlgMtZKPv/5ZA6oRnu7lM9avscejb0ekttZ3/ve8E+URdQunnNO8j3Z\naQNTDqXaIH3ZxKb02Zkxo7hxyXpZ3beOIsJNr4bwpS/B9dfXzvr/8peBcQS23ba+7Vaa+btSElM+\nk++pp9a/j6JFJLex3vOexk7cjWq0DUo776FX2neplqUVcdWT+C5enNwCfOtba5d761sHGpeWu/XW\n5H+8/ERZPkje+95X/XbI008nYz/tuy+cfHL1OGq1E4nojlm4K9l884H/iaGu1ksNuN/1ruH9H1Vr\n/1JL6YR71121y/3978lUCZVEFPP/f256Wf3668WPuVXy6KOw0kowdmxx29xqq/qPf5NNksdvfau4\n/fciT6lQoHqqcvN8gCt9OV9WYQ71v/998O/77guf+ETj+ytC6QrssceK33b2Pay3Crz0Hg71/j/2\nWOeNzJm3N1F5bUc9ozt//euw/PL1bX+jjSov33zz5LG8jVH556NWG4/S+CZ5brGUdNPwCrXakzW7\nt1ClGOr9n6t32P/VV6/e0LyeW3KNNHSu9RmvVUs21GshaZRbSi6sczmx6VLlV3HVbqOUn6SrfXCr\nXX3n1chJpd62BiXZtgz19qSo98p97bWTk/tw3Hzz8F5fLm81f/nJ8qGHhn7Nr36VPGYTxqK6af/k\nJ42/plk1Wdnt3nlnclumnbV53/9++/Zdkv07/9d/VS+XfZ8amal7ypTKy+sZB6ue3nWl23e1pkpo\n9LumklqjBZfXIL/6avLdM1QboqE6GnRSb61u4MSmxS7NMfFEpZNy+T96PQ1fa6kVVyNfXiWNTJcw\n1FVUuYMPbqx8Vj1fEMOdMLTWSKt5lBo8V4p9/vzk73P44UuvK6/Ba0S2dqcdE+mV/ufraQA93NtN\n++2XNI4vuoFmtb/XXnst3aakUuPgdtY4HXXU0ssqvc9F9bw78MDa66dNG3obpdtiQ3VjrvW5qPX9\nUM8QEOVzlJVq1mvNaA6NdzTolHngOpUTmxY75pjkyqiRL61KXyi1hkLP48Ybq6/Lc6XdaTq9rUV2\nIsVqKtXKjR2btNM58sil1w1nIL5jjhl4nmfE1qJOyvV2Sc1zwVAyY0bymJ1qolnGjk1qxf7wh6HL\nnnJK/i65u7docptOnjetmuyggI1oZOiFZst+3otOyHuBE5sCldoMDDVa7GuvVR9KXFr6JFxqZ5C9\nOip6sLX776++7uiji91XO3XqF/Fhhw1dpnQCrqTSmBfNmGCvXtnxa/JoNBGtJ1EYSq2eR3nUOoZ6\nbu9Nnpx/yIDzzqu8/B//SMaMuemmfNst10lj4dSr1thHEUn39XrbDlWb66qkqAuq8pqk7IVQ+WjO\n5sSmKbJD71fraXDRRfVvr/Qhq7dRZz1+97vBY4DU+0HuJOW9bGopfcFkZxqup3q7XrXm3GlUtTE2\nhtv2p5WqVfc38n/fiFLboG5R7xxJs2Y1tt0LL6y9/t3vTnoq1upUkHdqiGrdlvMqol1Jo4NkLl6c\ndF//0Y9qlyvNTVfrgrBcswbMLLp9ZC9wYtME2Z5AldqPDGeo7KJ88Yswfny7oxie4c6W/uUvFxMH\nwDe/Wdy2zjmnuG2VVLtKbVajxGqNYf/935uzv1o6/TZkka65pvb6ocYF2m+/geEoGlUabqKoqRpW\nHMYEPaU2KHlrvCrd3sl+r5fex6GmFMn+7w01/kzRNYYjmRObJqvUE6Udcx/1ove/v90RtM5wuj3D\n0mMhjaSTfR5DnbCaoYi/yQ03wCOP5H99o4O/VUqMi2qrNH9+8pjndtdyyxUTQ9b//d/g3xcuhLOG\nGG8/+zcdqudbrdGlrTFObEaYdnYbrPSFm/cefa1B2ypp1Yn8/PPzv7ZWI8CVV86/3bwefbR52/7w\nh2uvr/T3asUow6XbKd/+dvP31QyPPALrrtveGL73vWK3N25csdvLq3xambxTlhTVsL4dI6V3i65I\nbCT9UNIMSS9KmivpUklLXa9LOlzS05JelnSNpHXL1i8n6SRJz0laIGmapFXKyqwo6XxJ8yXNk3SG\npBWafYyw9MjBna40ENtQSslUpQ9i3vv5tUbeHCp5e+SR5nWl3W23/K8dbjfzor33vc3bdp55iT76\n0eL2X+1/ZKhuudY9ir6IK681rWcMrUoJeiNtA2s5+eTOaNbQiboisQG2BE4ANgE+BbwZuFrSPwd+\nl3QgsB+wN7AxsBCYLmlUZjvHAZ8BdgS2AlYHyjv/XQCMByalZbcCTms04PJqy3rU0/iuqA9FyXA+\n/LfeWt8YJ7VmNR+qsWNRsl8w3/hGcbMEF+nUU7uzEXczFFHDduSRcFqDn9yiG8Ba+0yfXuz2ihg1\n/rXX4LOfLSaeU0+tPk3FSNekGTWKFRGfzv4u6askk25OBEr9XPYHjoiIy9MyuwNzgS8AF0oaDewF\n7BwRN6Vl9gRmSdo4ImZIGg9sC0yMiLvTMpOBKyQdEBF1z2G70UbNqRVYa61ir0Suvx4mTcr/+qHu\nMUMyDsqvflW5YeMFF+TfdyOyXzDXXju8MV6qGeq22lB/tyuvbGySy3oN56ru+uuLi6OkVV2Ef/zj\n1uzHBhQ98epwJrRsdLb1oQw1jEcl8+YN/n3DDfPPom7165Yam3JjgQCeB5C0NrAq8M+7oBHxInAH\nUKpc3pAkkcuWeRB4IlNmU2BeKalJXZvuq+EZQhrpCtgu5TP8NlOj7WK6Tae0BShXaRyiUpI11Ilj\nOElvNXnep1bWZP3wh8Vs53OfK2Y7nSI743g1KxR8037ChPyvzQ7t0C7lc+U5qWmNrktsJInkltL/\nRkRplJhVSZKP8uHX5qbrAMYBi9KEp1qZVUlqgv4pIhaTJFCr0qDhjIZayxFHFLu9obqAFqHdc52M\n5F5AlSa/LNVYDXcgvVYZzpV7oypNJ2ADWvlZ7vYGsiP5e6eduuJWVJmTgQ8CdTZdbY0pU6YAg+8j\nHHJIH1DnSFwNOPPMYrdXdPVxJb00enGj8g6L30zZv/mdd1afqbtdyucna3diPBwHHZTMhfSmbvy2\nHUIvD+dfGgvn2GNrTwpq9evv76e/v3/Qsvmlfv0F6qqPmqQTgU8DW0ZEdmzTOYBIamWytTbjgLsz\nZUZJGl1WazMuXVcqU95LallgpUyZiqZOncrEicOoN23A44/nf22lE8RQo2wW4ZJLmr+PWprZdXko\n2Tlmnn026f20667ti6fcxht3XuJQ9t1XeKP5hx9u7a3DBx/s/gExS2bNgvXWS+a9atf8SdlBAG+5\npf4emo0o3Xo76KD8iY17LQ3W19dHX9mw2zNnzmRiPV3MGtA1t6LSpObzwCcj4onsuoj4K0niMSlT\nfjRJu5jSneG7gDfKyqwHrAWUOh3fBoyVlJkUgUkkSdMdRR5PO5WP2Jt3PIZGtDOxgKXHoGiVRx4Z\nPJLwV74yvC7hI0Wz2yJssUVzGkZXc/75cPnlrdtfM2266cCowFde2Z4Y3v72gec//Wl7YqjHU0+1\nO/F+jTEAABL0SURBVIKRqStqbCSdTHJP53PAQkmla635EVE6TR8HHCTpYeAx4Ajgb8DvIWlMLOlM\n4FhJ84AFwPHALRExIy0zW9J04HRJ+wKjSLqZ9zfSI6rTlTeOvPnm4rZ95ZXw6U8vvfy554rbRzW1\nah3mtOmvVz5YWqmXRERr7r9Xa3R7wAGDf1+wYPDJoh223hruvjv5Xymv4SuPtwiTJrWupqrS7OvD\nccopScP/972v2O02qh2fq/Jbu+26aBlKBFx2WXHbsvp1RWIDfIOkcfCNZcv3BM4FiIijJS1PMubM\nWOBmYPuIWJQpPwVYDEwDlgOuAsqHetsFOJGkN9SStGzPzOJx++1w3HHN2/5nPpN8CJ94YuiyRSvN\neFtp4tF2DJFfyTJpHWmrEpt6u48//jh86EPNjWUotbrgd+rJq12++c1kQM9Wfc4WLRq6TKv86U+D\nf+/Udj7l05hY63RFYhMRdd0yi4hDgUNrrH8NmJz+VCvzAtCzNwuKvnKsZNEiOOOM5u+n3OGHt36f\njbr99uTxr39t7si+jXrppdb0jsvj9debt+1OOmE3quh2R7U0Y+6lvKpN6FqkIsYgK3JaDtfYNKZr\n2thY99hySzj33HZH0dn22afdEQy22Wbwtre1O4rKRo0aukxenXTCrteSJQPP99uvfXH0suFMJNoM\nv/pVuyPoLk5srHAzZrQ7gs7nWyvVuSdJbdkZuBudjbtI7bjdXE2v12g4sWmME5sOVz7xmtVPas8t\nsXqVtxWwhOe/qe3i8tnt2uSww9odwYBlCj6TlX82d9ih2O1bczmxKUizpk9YeeXmbHek+P732x1B\ndR//eLsjsG50003tjiBRzzxxvaITJs1txUCqvcKJTUF2373dEVhWaeCwF15obuNTM7NWKHoerl7m\nxMZ60gaZIRY7aZRfMzNrLic21vMuuqjdEZiZWas4sTEzM7Oe4cTGzMzMeoYTGzMzM+sZTmzMzMys\nZzixMetwS5b0/siqZmZF6YpJMM1GKic0ZmaNcY2NmZm1hRN3awYnNm1y/PHtjsCsOSJ8wjKz9nFi\n00KvvjrwpT95crujsW7z/PPtjqAxe+/d7gislgsvbHcEZs3hxKaFlluu3RFYs1x9dfP3seKK1dfN\nmFF9XbtqT448sj37taFFwJe/3O4ozJrDiU0Xuvvu1u7vj39s7f660dZbtyeBKNUAbrRR7f23I7aV\nVmr9Ps3MnNg00QUXFL/NDTeEtdcufrvVRMB2242sNhMLF8Kpp+Z7bSe/T5tt1vx99Pc3fx/1+N//\nhaeeancUZtYOTmya5L//G/r6it/uhAkwZkzx2zXYf/8kMVl+edhnH7juunZHNHyrr97a/e28c2v3\nV83mm7f+2LuRkz/rRU5smuDVV+FHPxq63E9+0vi2v/e9xl9jQ4uA444bvExqTyxFOvPMgecf+lD7\n4miXCDjjjMZec/PNnV3zVqROTv4iYPTodkdh3ciJTRNkGwmPGlW93Le/3fi23//+xl9TlO9+t337\nziq1K+nFtj8RtXurzJrV2Pa23nrgeSecxEp/u1YmDk8/XX/ZJUtgiy2aF4s1Zv58eOihdkdh3caJ\nTZPVamfz9rc3tq12X0Uec0x7919uu+2au/23vKW5288aM2bg71urt8qaa9bezuOPD/592WUHno8d\nmy+24bj33mK2c+21+f//N9mk/rK9UEvXThFw441LL+/ry5/UrrtuIaF1jEWL2h1B73NiU4Gkb0n6\nq6RXJN0uaaO829pwwyIja47Fi+sr198pLUMz6o29EaXjHDcu3+u/+tXGX/PCC/WVW2GF2utXXrn6\nure9rXxJ8/+eH/7w8LcRAZMm5X/97bfXd5ytbJTfPK37jJbXpJQSlq22Glh28cXJ8uI7UtQ+zlIC\ndeWVxextzpzhb2PBgiSmN7+5kVdVP863v92DYVbjxKaMpJ2AY4BDgA2APwPTJdU4ZVS31loFBtcE\nCxbAMnX+FzQjsVmwYHivryf2Rj/4pePMW2Pzq18NXeb55+H++/Ntv5bll6++7p3vLF/SeYlqHqUv\n92q3Sm+4ob7jLG9j1Z3y/00j4Nln6y9frSZFGkhu1lgjdzhDGHycS5YMPH/llYHnRdXq5r3IyVr6\nwqIeA8eZ/R574AF48cVhh9SznNgsbQpwWkScGxGzgW8ALwN75dlYI1Xbzz2XZw/D0+iH7Qc/aO/+\n61X6ElhnnfzbaOatmxVXhA9+sHnbr+Qd72isfK1bB1LjV4unn97Y/odSvv9f/KJyucaukEe2lVdO\n3tPvfGd429lgg2LiKcn2touAbbYZvD77PVvteSfKOzr3+PHFxtFrnNhkSHozMBH4Z0ffiAjgWqDp\no4BUOvF0WjVjqTHq/2/v3mPlKMs4jn9/XJQCKSYKLYlIIwpGwVZqC15KK4ImQlBAsAGsYr1gVdDE\nCwaVUOItYhEVElELNvWSmqBWlIsFMUHAYg9tSZCCIuKlgFpSKK3Qy+Mf76xnzp6d7V7P7oy/T7I5\nZ2femfM++5ydffed95056aTiMkV1LpoFtnVra9u36/HHYcOG9PvuTuE0kh/4HZGu7rtu3egH6qmn\n9qaejWzcCE880b/9N5J/3YtyUIs9/w25VbsbH9Stog+xYf9wG0aXX97d9rX3W68alfWdxc322+wK\n79u3D2YAe5H6Bpr1hhs2Y70A2BN4rG75Y8DUTne6fj3ce2831RoetTEtRx/dvNwvfjF+WdGsnEmT\nxneB5w8+nR6EDjwQ9tor/d7K9Pt6tW1rZs2CV75y9Hlt6n2j6cTnnNP+38ubOrX7HqP6K/82G3/T\nzVilVrvp63vnik6bHXFE53VpJD+Aupl583r7d/8fFN1K5DOfSRe57LTnZnenaRs1VkdG4Lzzxi9/\n5hnYsiUdQ+rf0+3oxxDDdm+zM3du7+tQRV2k2TLZSIzRebgjI40LNlq+ahWccAKcdVZav2ZNagid\ne+74bYr2241m+1y0aOz6zZs3s21bWnDwwc33WfsQveiiNJvqwgsbDyTN73/aNJgzp3mdzjijtdeh\nvsz06enniSfCr341unz27FSv+jhHsgVXX53GEjT6m5MmpUGRhx8+fv38+bB8+fhtVq1Kjbj68u3k\n9rvfhYULU4Oq2XY//enY9Tt2pJ8zZ6b/M9gMpAJr1zavz6teVfw6rFiRDri33NJ4/T77jN4ANr/+\nYx9LF7LMe/7z02t6993pQ2ry5LHbTJ6cxhYsXdo89vy6LVtG46xZvXq0wbNpU+qJ++Mfi/dXZN99\nx/c4Fpk3r/GMoXbtsUfqMVu6FN4z7gT5+Fgbyb8+y5enY07961k7Nt12WxrTsWjR6LpLL03lDzgg\nNcJrz/Nmzdr97V8a5XDNmvT/0uw4unnzaJznnz9a9n3va++9dN11cNppjdctWADLlqVj2chIOj6d\nfTZ88IPpvX/TTe19YcrX68or4UMfauW2IynOM89M2y9ZMn5fAIceOnZG5MknwyWXpN9nzkw/f/c7\nWLx49Evn7Nmj95jrx2dLq/4weg2Lns1DVQxDf9yQyE5FbQVOj4iVueXXAgdExLiTD5LOAr4/YZU0\nMzOrnrMjoifz59xjkxMR2yWtAd4IrASQpOz51ws2uwk4G3gY+M8EVNPMzKwq9gGmkT5Le8I9NnUk\nnQlcS5oNtZo0S+rtwMsioo3JkGZmZjbR3GNTJyJWZNesWQxMAdYCb3ajxszMbPi5x8bMzMwqw9O9\nzczMrDLcsDEzM7PKcMOmBe3eFFPSPElrJP1H0gOS3jVRde1GO3FKmitpV91jp6SDJrLO7ZI0R9JK\nSX/P6nxKC9uULp/txlnGfEr6tKTVkp6U9Jikn0g6vIXtSpXPTuIsYz4BJJ0naZ2kzdnjDklN7/hU\ntnxC+3GWNZ/1JF2Y1X3Jbsp1lVM3bHaj3ZtiSpoGXE+6LcN04ArgO5JOnIj6dqrDm38G8FLSVZmn\nAgdHxOP9rmuX9iMNCF9Eqn9TZc0nbcaZKVs+5wDfAI4BTgD2Bm6WNKlog5Lms+04M2XLJ8BfgU8B\nR5Nub3Mr8DNJDe+OVNJ8QptxZsqYz//Jvii/n/TZ0qzcNLrNaUT40eQB3AVckXsu4G/AJwvKfxlY\nX7fsh8AvBx1Lj+OcC+wEJg+67l3EvAs4ZTdlSpnPDuKsQj5fkMX6+orns5U4S5/PXCz/Bs6taj5b\njLPU+QT2BzYAxwO/BpY0Kdt1Tt1j00SHN8U8Nlufd1OT8gPXYZyQGj9rJf1D0s2SXtvfmg5E6fLZ\nhbLn83mkb7WbmpSpQj5biRNKnk9Je0iaD+wL3FlQrPT5bDFOKHc+rwR+HhG3tlC265y6YdNcJzfF\nnFpQfrKkNm95NmE6iXMj8AHgdOA0UtfqbZJm9KuSA1LGfHai1PmUJOBrwO0RcV+ToqXOZxtxljaf\nko6U9BTwDHAVcGpE3F9QvLT5bDPOMudzPjAD+HSLm3SdU1+gzzoSEQ8AD+QW3SXpMNKVmod+8J6N\nVYF8XgW8HHjdoCvSZy3FWfJ83k8aW3EA6arvyyQd1+RDv6xajrOs+ZT0QlJD/ISI2D5Rf9c9Ns39\ni3Rec0rd8inAowXbPFpQ/smIeKa31euZTuJsZDXwkl5VakiUMZ+9Uop8Svom8BZgXkRs3E3x0uaz\nzTgbKUU+I2JHRDwUEfdExEWkwaYXFBQvbT7bjLORMuRzJnAgMCJpu6TtpPFCF0h6NuuBrNd1Tt2w\naSJrYdZuigmMuSnmHQWb3Zkvn3kTzc+dDlSHcTYyg9RlWiWly2cPDX0+sw/7twJviIhHWtiklPns\nIM5Ghj6fBfYAik5BlDKfBZrF2UgZ8rkKOIpU1+nZ4/fAcmB6NpazXvc5HfRo6WF/AGcCW4EFwMuA\nb5FGrx+Yrf8i8L1c+WnAU6SR3UeQpts+S+qKG3g8PYzzAuAU4DDgFaTuxu2kb5MDj6dJnPtlb64Z\npJklH82eH1KxfLYbZ+nySTot8wRpOvSU3GOfXJkvlD2fHcZZunzm4pgDHAocmf2f7gCOL/i/LV0+\nO4yzlPksiH3MrKh+vEcHHmQZHtkL+zCwjdRqfHVu3TXArXXljyP1gGwDHgTeOegYeh0n8IkstqeB\nf5JmVB036BhaiHEu6YN+Z91jaZXy2W6cZcxnQXw7gQW5MqXPZydxljGfWb2/AzyU5eZR4GayD/uq\n5LOTOMuaz4LYb2Vsw6bnOfVNMM3MzKwyPMbGzMzMKsMNGzMzM6sMN2zMzMysMtywMTMzs8pww8bM\nzMwqww0bMzMzqww3bMzMzKwy3LAxMzOzhiTNkbRS0t8l7ZJ0SpvbX5xttzP7WXs81a86u2FjZmZm\nRfYD1pKuTN/JFX2/AkwFDs5+TgXuA1b0qoL19urXjs3MzKzcIuJG4Eb4382Rx5D0HNL9nuYDzwPu\nBS6MiN9k228l3YewVn468HLg/f2qs3tszMwakPRrSUsGXQ+zIXclcAzpRspHAT8GbpB0WEH59wIb\nIuKOflXIDRszsxxJew+6DmZlIOkQ4N3AGRFxR0T8OSKWAL8Fzm1Q/rnAWaSbgPaNGzZmNrSyXpMr\nJH1Z0r8lbZR0cW79LknnSfqlpK2S/iTp9Lp9fEnSBklPZ+sXS9ozt/5iSfdIWijpIWCbpGtId0i/\nIDfw8UUTFrhZORwF7Ak8IOmp2oN0d+5GPTanAfsDy/pZKY+xMbNhtwBYAswGXgtcK+n2iLglW78Y\n+BRwflb2R5KOjIgN2fons+UbSQfib2fLLsv9jZeQDrqnAjuBvwCHk8YLfBYQ8M9+BWhWUvsDO4Cj\ngV1167Y0KL8QuD4i+vpecsPGzIbd+oi4NPv9T5I+DLwRqDVsVkTENdnvn5N0IvAR4MMAEfGF3L4e\nkfRV4B2MbdjsDbwzIjbVFkh6Ftja74OwWYndQ+qxmRIRv21WUNI04A3Ayf2ulBs2Zjbs1tc93wgc\nlHt+V936O4HptSeS3kFq6BxG+oa5F7C5bpu/5Bs1ZpZI2o/Uo1mbEfXibGbTpoh4UNIPgGWSPk5q\n6BwEHA+si4gbcrtaCPyDbIZVP3mMjZkNu+11z4MWj12SXgMsB64HTgJmAJ8HnlNX9Oku62hWVa8m\nNVjWkN57XwVGgEuy9e8mjZm5DLgfuC7b5pHaDrJp4u8CromITq6F0xb32JhZ2R1Larzkn49kv78G\neDgivlRbmXWJt+JZUje72f+t7Ho0hV8kImInqZFzSZMyAUzY4Hs3bMys7M6QtAa4HTgHmMXoVNMH\ngRdlp6PuJp3ff1uL+30YOEbSoaSBkJsm4tummXXHp6LMbJgVNSTyyy8mXfV0HalhM782Iyoifg5c\nDnyD1J1+LGkWVSsuI82Qug94HDik3cqb2cSTv4CYWVlJ2gW8LSJWDrouZjYc3GNjZmZmleGGjZmV\nmbuczWwMn4oyMzOzynCPjZmZmVWGGzZmZmZWGW7YmJmZWWW4YWNmZmaV4YaNmZmZVYYbNmZmZlYZ\nbtiYmZlZZbhhY2ZmZpXxX8OlY/uAbDC7AAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f8bec24b5d0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ds['delta_dist'].plot()\n",
"plt.ylabel(\"Distance Moved in 1 Day [km]\")"
]
},
{
"cell_type": "code",
"execution_count": 88,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f8bb810f4d0>"
]
},
"execution_count": 88,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAk8AAAGHCAYAAACplLYqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XucXWV59//PNwHCD0oCNCVISypWG0YeUDIcQpVTY8OD\nUGurv58MzQMiVkQEjRVRHw4RsAKtBCWIPoSTIGMtarXAjyCohAIFyaCCDLGWw4CQ0F2SCU0YQpLr\n+eO+d7JmZc9hT+awd+b7fr32a7Lvde217r0SnYv7cC1FBGZmZmY2OBPGugNmZmZmzcTJk5mZmVkd\nnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk9k2\nTNIzkq4b635s6ySdLek/JK2X1DGEz/9U0o+HeO2G+TuW9EFJGyVNL7T9VNJPxrJfZsPNyZNZk5B0\ncv7FNLOP4z+V9MtS80agrmcwSTpW0gVD7ed4I2kOcClwH/BB4PNDOM3WPCer12cltUi6oJjAjKIo\n9ye/31jPSSQdlr/D5GHrmdkw2m6sO2Bmdenvl2ytYzOo8xcX8G7gY8AX6vzceHU0sAE4NSI2jHVn\ngLcCFwA/AbrGuC8AfzaEz/wJcD5wPbB6eLtjtvU88mS2DYuI14fwC10j0pnBXlzaaSyvPwTTgFcb\nJHGC9PfXME98j4j1EbG+zo+N6b9Bs4E4eTLbhpXXw0jaLk+H/FrSq5Iqku6TNDsfv5406kSeItwo\naUPh8ztJ+rKkLkk9kp6U9Lc1rrujpK9K+k9JqyX9s6S98vnOL8TNz20tkm6R9DJp+gtJ+0u6Pq8l\nelXSi5KulbR76VrVc7xF0s2SVkl6SdKF+fje+frd+RyfGuS9myjpPEm/yd/1aUlflLRDIWYjcDKw\nc/VeSTppgPN+JJ9zraR/k/TOPuJ2kPQFSf+er98l6dLi9Wt85mTgO/ntTwt9OiIff4+k2yT9Np/z\nN5LOlVT37wJJb5X04/w9npP0v6nxO6XWei5JZ0p6XNIaSS9L+pmkE/KxC4DLcugzhe8wFtOQZjV5\n2s6s+UyR9LulNgHb14gtj0B8Afgs8H+AnwGTgYOAmcA9wNeBvYB3AX/NliMA/wIcCSwCfgEcA/y9\npL0iophE3Qi8H/gm8FD+zO01+lN9/0/Ar4HPFa75Z8A+wHXAcmA/4DTStNRhNc7xj8ATwDnAccD/\nzsnYafm7fSZ/p7+X9HBE/Gv5ZpVcC5xESkb+ATg0929f4H05Zm4+/8HAqbnvD/R1Qkmnku7xvwIL\ngDcBPwRepjDFJkmke/0nwDeAJ4H9gXnAW4C/6uMSS4CvAmcCF+fPAXTmnx8EXgG+DPw38KfAhcAu\npPs2KJKmAT8lJUt/B6wFPgL01Agvr8n6G+ArpPt6BbAjcADp/n4b+B7wx8AJwCeA/8of/c/B9s9s\nxEWEX3751QQv0gjHxgFevyx95mngusL7R4EfDnCdK4ENNdr/Il/js6X27wDrgX3y+wNz3D+U4q4j\nrQ06v9B2QY69qcb1JtVo+0A+xztqnONrhbYJpGRkPfDpQvsUYE3xnvRxDw7I5/x6qf2yfP0jC23X\nA6sH8fe3HSkJfATYrtB+ar7Wjwttc4HXgcNK5/hIvv6sfv6O35djjhjkPb2alFBtX8e/xQX5Gq2F\ntt8FVub26YX2n5S+2/fL/05rnP9vy+fxy69Gennazqy5BHA6aWSo/CrvtKtlFbCfpDcP4drHkpKR\nK0vtXyYlK8cW4oL0S7noSmqvZQnS6ErvxojXqn+WNCmPtj2Uz1HecRikkaLqZzeSkhSRkrZqezew\njDTi059353MuKLV/OZ/zuAE+X8tBwB6khKy4BuhGoLsU+37SaNGvJf1u9UVKRERapF630j39nXzO\nfwV2Io2oDdaxwL9FxNLCuf8L+NYgPrsK+ANJB9VxPbOG4mk7s+bzs4jYopaQpJWk//rvz/nAP5N+\nKT8O3Eka9XlsENf9Q+CFiFhTau8sHAeYThpJeboU95t+zl2ORdJuwHzSaNMehUNBGkEqK+8s6wZ6\nIuLlGu27078/JH2HXn2OiBWSVrH5u9bjD0l9L59zvaSnSrFvISUztaaqgt73Y9AkvRX4Iin5KpYB\n6Oue9uUPgX+r0b5sEJ+9FJgNPCzpN8BdwC0R0ed0p1mjcfJkNo5ExH2S/og0BTeHNGU0T9JpETGW\nhRZfrdH2T8As0lTZL0hrdCYAi6m92aXWbre+dsANdjfXWO1amwA8RlrjVKuvz9V7QklTSGuiVgHn\nAk+R1ii1ApcwShuIIuJJSTOA44H/SVq/9TFJX4gIl8ewpuDkyWyciYhVpKmiG5XKAtxHGuGpJk99\nJQzPArMl7VwafWrJP58pxE0gLfb+j0LcWwbbR0m7khYznxcRXyy0D2W6cSiq3+EtFEZTJO0B7JqP\nD+Wcyuf8aeGc25Hu1c8Lsf8BHBARQ6nM3dff31HAbsBfRMT9hev/0RCu8Sy1/z4HNfUXEa+SkuN/\nyt//+6QF/l+KiHU0UKkFs1q85slsHClv84+ItaRppEmF5jU5tlzd+Q7Sf3B9vNQ+jzTFdWd+v5iU\nJHysFHcmg/+lWB0xKv9/1Lw6zrE17iB9h0+W2v82X//2IZzzEdI03EdzwlB1CikhK/oOaV3Q35RP\nkstA9FcLaw2p7+Vzbsjtm+5pLntQ/nsajDuAWcV1S5J+DzhxoA/W+De4njT1W9wxWk3Oy9/BrCF4\n5MmsuWxt8cAnJP0UWEraHn8waXHyVwsxS/N1rpS0mLTz7h9JW+d/AnxR0j5sLlXw58CCiHgaICI6\nJH0X+KSkqaS1MUeyeaRiwOQnIl6RtAT4TP4F/1vSNOMbGYUCihHxS0k3Ah/Ja6/uJW2lPwn4XkTc\nO4Rzrpd0LqlUwU8k/SNpxOkUeo/QAdwE/H/A1ZKOBu4HJpJG+f5f0r3o6xl6PyclSufkEbzXSKUa\nHiDthvumpOrf91yGloxeBvwvYLGkr5BKFfwNafTxgAE+e5ek5fk7rSCVnjgDuK0woln9N/h3kr5N\n2nn4wzxiZTbmnDyZNZeBftHVqqNUbPsK8B5SDaVJpOmXz5PqGFV9j5RMncDmWk//GBEh6c9JdYE+\nQKoZ9AypFEB5V9r/Al4E2oC/JP3yPoE0BVarFlAtbaQdeh/LfVhM2uX1Qo3v2Ze+4gbz+VNJSc0H\ngfeSygx8kfT9h3I+IuKaXJDybFIC8hgp+byoeI58r/+CNNJ2Ur7+WtI6pQWkmljFaxc/u0LSaaSa\nVItISdfREbFE0nGkHYMXkRKpm4Afk+7toEXEcklHkf5+ziHVYrqadI8W1fpI4c9fJ/27mgf8DvA8\nqd7TpunZiHgkJ5ofJSXo1WngRnjcjBmK8NSymY08SW8njZb8dUS0j3V/zMyGqiHWPEk6XNIP8yMD\nNkp6T+n4zpIW5kcArJX0q/xfVsWYSZKuUnrcxCuSbs2LO4sxu0n6ltJjGlZKWiRp51LM3pJuz48N\nWC7psvKjCyQdIGmJ0iMjnpV0do3vdJSkpUqPQPi10mMTzMYFSTvWaP4kaTppySh3x8xsWDXKtN3O\npHn6a0lTBmULSDtFTiRNM8whrQX4bUTclmOuIA3pv4/0FO6rgO8ChxfOcwvpIZ6zgR2AG0jF+eYC\n5CTpDtK0wCzSYypuAtaRtvYiaRfSEPddpMcy7A9cL2llRCzKMW8EbgO+lvv8LmCRpBci4kdDuUFm\nTeYzklpJa6TWk4pOHgN8IyJ+O6Y9sz7lpHegek8vR8Tro9Efs0bVcNN2Sg/afG9E/LDQ9hjw7dKW\n5UeAOyLi/Lwr6D+BEyLi+/n4DNIOjlkR8bCkFuBXpMcJPJpjjiHtmvmDPId/LOk5U2+IiEqOOY1U\nA+X38oLP00nrBfasVgmW9CXS9t+35veXAsdGxKaFk5LagSkR8e7hv2tmjUXSu0gFOd9KWtfSRXrO\n3d/l6t/WgPII+fX9hAR5/dQodcmsITXKyNNAHgDeI+n6iHgh7z55C5sXObaSvss91Q9ExDJJXaQH\niD5MGklaWU2csrtJ/2dwKPCDHPNYNXHKFpMWQu5H2l00C1hSerzCYtJ/aU/Jj3+Ylc9NKaa8qNZs\nmxQRd7Pl/was8d1JGinvzy9GoyNmjaxZkqczSU+Bf17SetK6ib8pFHrbE1gXEatLn1uRj1VjXioe\njIgNSk9dL8asqHGO6rFf5J/lRykUY7r7Oc9kSZOKz5cyM2sUEbGCLf+/y8xKmiV5Oos0OnQ8afj/\nCOBreQ3Rj8e0Z5ttVe0ZpQd0HkPa+j3YrdxmZmYGO5LqwC3OD6keUQ2fPOUFjF8krYP6/3Pz45IO\nBD5NqlGyHNhB0uTS6NO0fIz8s7z7biLpAaHFmINLXZhWOFb9Oa1GTAwiZnU/o07HMLgnkpuZmVlt\nf03aHDaiGj55IpXr354tH/C5gc2lFpaSdvTMJj0jqbpgfDrwYI55ENhV0oGFdU+zSSNGDxViPi9p\namHd0xzSVNwThZiLJU2MiA2FmGV5vVM15thSf+cU+lLLMwA333wzLS0t/YRZ0bx581iwwEvJ6uX7\nVj/fs6Hxfauf71n9Ojs7mTt3Lmx+xuaIaojkKddaejObp77eJOltpC2xz0m6F/gHSWeSShUcRaq6\n+0mAiFgt6VrgckkrgVdIFZLvj4iHc8yT+VET1+QdczuQquO2R0R1xOguUpJ0k6RzgDeQdtYtLGzN\nvYW0i+i6vKtuf9K04icKX+nrwBn5+HWkJO39pO3afekBaGlpYebMmXXcvfFtypQpvl9D4PtWP9+z\nofF9q5/v2VYZlWUvDZE8AQeR6sFUHzPw5dx+I/Ah0qMgvgTcTJpmexb4XET8n8I55pFGo24lPXbi\nTtLzkopOBBaSdgFtzLGbkp6I2CjpeNLuugdID6e8AbigELNa0hxSHalHgAowPyKuLcQ8kx+DsICU\nWD0PnJp3IJmZmVkTa4jkKT9ks89q5xHxEuk5U/2d4zXSrrwz+4lZRS6I2U/Mc6SF6f3FPE560Gl/\nMUtIJRTMzMxsG9IQj2cxMzMzaxZOnqyptbW1jXUXmpLvW/18z4bG961+vmeNr+EezzJeSZoJLF26\ndKkXCpqZmdWho6OD1tZWSI9g6xjp63nkyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDkyczM\nzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDk\nyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM6uDkyczMzKwOTp7MzMzM\n6tAQyZOkwyX9UNJvJW2U9J4aMS2SfiBplaT/lvSQpD8oHJ8k6SpJFUmvSLpV0h6lc+wm6VuSuiWt\nlLRI0s6lmL0l3S5pjaTlki6TNKEUc4CkJZJelfSspLNr9PcoSUsl9Uj6taSTt/5OmZmZ2VhriOQJ\n2Bn4OfAxIMoHJf0RcB/wBHAEsD9wEdBTCLsCOA54X47ZC/hu6VS3AC3A7Bx7BPCNwnUmAHcA2wGz\ngJOBDwIXFmJ2ARYDTwMzgbOB+ZI+XIh5I3AbcA/wNuArwCJJfza422FmZmaNarux7gBARNwJ3Akg\nSTVCLgZuj4jPFdqerv5B0mTgQ8AJEXFvbjsF6JR0SEQ8LKkFOAZojYhHc8yZwO2SPh0Ry/PxfYGj\nI6ICPCbpPOASSfMjYj0wF9geODW/75R0IPApYFHu0unAUxHxmfx+maR3AvOAH23NvTIzM7Ox1Sgj\nT33KydRxwL9LulPSCkn/JukvCmGtpETwnmpDRCwDuoDDctMsYGU1ccruJo10HVqIeSwnTlWLgSnA\nfoWYJTlxKsbMkDSlEHN36assLvRlyLq6uujo6Oj16urq2trTmpmZ2SA1xMjTAPYAfgc4B/jfwGeA\nY4HvSToqIu4D9gTWRcTq0mdX5GPkny8VD0bEBkkvl2JW1DhH9dgv8s+n+onp7uc8kyVNiojX+v3G\nfejq6mLGjBZ6etb2at9xx51YtqyT6dOnD+W0ZmZmVodmSJ6qo2P/HBFfzX/+paQ/AT5KWgvVCGpN\nNw6rSqWSE6ebSUu3ADrp6ZlLpVJx8mRmZjYKmiF5qgDrgc5Seyfwjvzn5cAOkiaXRp+m5WPVmPLu\nu4nA7qWYg0vXmVY4Vv05rUZMDCJm9UCjTvPmzWPKlCm92tra2mhrayu0tJDWqpuZmY0v7e3ttLe3\n92rr7u4e1T40fPIUEa9L+hkwo3Toj4Fn85+XkhKs2cD3ASTNAKYDD+aYB4FdJR1YWPc0mzRi9FAh\n5vOSphbWPc0hTcU9UYi5WNLEiNhQiFkWEd2FmGNL/Z1T6EufFixYwMyZTozMzMxq2XJAATo6Omht\nbR21PjTEgnFJO0t6m6S356Y35fd75/d/D3xA0ocl/ZGkjwPHA1cB5NGma4HLc32lVuA64P6IeDjH\nPElatH2NpIMlvQO4EmjPO+0A7iIlSTflWk7HkEoiLIyI13PMLcA64DpJb5X0AeAs4MuFr/T1/B0u\nlTRD0seA9wOXD+NtMzMzszHQKCNPBwE/IU19BZsTkRuBD0XEP0v6KPB5Us2kZcBfRURxJGcesAG4\nFZhEKn1wRuk6JwILSTvhNubYT1QPRsRGSccDVwMPAGuAG4ALCjGrJc0hJW6PkKYV50fEtYWYZyQd\nBywgJVbPk0oblHfgmZmZWZNpiOQp12bqdxQsIm4gJTJ9HX8NODO/+opZRarT1N91niONavUX8zhw\n5AAxS0glFMzMzGwb0hDTdmZmZmbNwsmTmZmZWR2cPJmZmZnVwcmTmZmZWR2cPJmZmZnVwcmTmZmZ\nWR2cPJmZmZnVwcmTmZmZWR2cPJmZmZnVwcmTmZmZWR2cPJmZmZnVoSGebWeNoauri0ql0qtt6tSp\nTJ8+fYx6ZGZm1nicPBmQEqcZM1ro6Vnbq33HHXdi2bJOJ1BmZmaZp+0MgEqlkhOnm4Gl+XUzPT1r\ntxiNMjMzG8888mQlLcDMse6EmZlZw/LIk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZ\nmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZmVkdGiJ5knS4pB9K+q2kjZLe00/s13PMWaX2\nSZKuklSR9IqkWyXtUYrZTdK3JHVLWilpkaSdSzF7S7pd0hpJyyVdJmlCKeYASUskvSrpWUln1+jn\nUZKWSuqR9GtJJw/t7tSvq6uLjo6OXq+urq7RuryZmdk2rVEez7Iz8HPgWuB7fQVJ+kvgUOC3NQ5f\nARwLvA9YDVwFfBc4vBBzCzANmA3sANwAfAOYm88/AbgDeAGYBewF3ASsA87NMbsAi4G7gNOA/YHr\nJa2MiEU55o3AbcDXgBOBdwGLJL0QET8a5D0ZEj/g18zMbGQ1RPIUEXcCdwJIUq0YSb8PfAU4hpTg\nFI9NBj4EnBAR9+a2U4BOSYdExMOSWvJnWyPi0RxzJnC7pE9HxPJ8fF/g6IioAI9JOg+4RNL8iFhP\nSrS2B07N7zslHQh8CliUu3Q68FREfCa/XybpncA8YESTp94P+G3JrZ309MylUqk4eTIzM9tKDTFt\nN5CcUH0TuCwiOmuEtJISwXuqDRGxDOgCDstNs4CV1cQpuxsI0mhWNeaxnDhVLQamAPsVYpbkxKkY\nM0PSlELM3aU+Li70ZRRUH/A7k81JlJmZmW2tpkiegM8C6yJiYR/H98zHV5faV+Rj1ZiXigcjYgPw\ncilmRY1zMEwxkyVN6uM7mJmZWRNoiGm7/khqBc4CDhzrvgyg5nRjvebNm8eUKVN6tbW1tdHW1jYc\npzczM2tq7e3ttLe392rr7u4e1T40fPIEvBP4PeC5wnKoicDlkj4ZEW8ClgM7SJpcGn2alo+Rf5Z3\n300Edi/FHFy6/rTCserPaTViYhAxqyPitb6/KixYsICZM2f2F2JmZjZu1RpQ6OjooLW1ddT60AzT\ndt8EDgDeVni9AFxGWuANsBRYT9pFB4CkGcB04MHc9CCwa17cXTWbNGL0UCFmf0lTCzFzgG7giULM\nETnxKsYsi4juQsxseptT6IuZmZk1qYYYecq1lt7M5qmvN0l6G/ByRDwHrCzFvw4sj4h/B4iI1ZKu\nJY1GrQReAb4K3B8RD+eYJyUtBq6RdDqpVMGVQHveaQep/MATwE2SzgHeAFwELIyI13PMLcD5wHWS\nLiWVKjgL+EShi18HzsjHryMlUu8H3r2192osdHZuXqM/depU79gzM7NxrSGSJ+Ag4Cekqa8Avpzb\nbySVICiLGm3zgA3ArcAkUumDM0oxJwILSTvhNubYTUlPRGyUdDxwNfAAsIZUC+qCQsxqSXNIdaQe\nASrA/Ii4thDzjKTjgAWkxOp5UmmD8g68BvciMIG5c+duanG9KDMzG+8aInnKtZkGPYWY1zmV214D\nzsyvvj63ilwQs5+Y54DjB4h5HDhygJglpBIKTWwVKces1oxyvSgzM7OGSJ6s0VVrRpmZmVkzLBg3\nMzMzaxhOnszMzMzq4OTJzMzMrA5OnszMzMzq4OTJzMzMrA5OnszMzMzq4OTJzMzMrA5OnszMzMzq\n4OTJzMzMrA5OnszMzMzq4MezbCM6Ozt7/RyNa1VNnTrVz7ozM7Nxw8lT03sRmMDcuf0+73hEr7Xj\njjuxbFmnEygzMxsXPG3X9FYBG4GbgaXARaN4raXAzfT0rKVSqYzgdc3MzBqHR562GS3ATGDkp+02\nX8vMzGz88ciTmZmZWR2cPJmZmZnVwdN241hXV9emtUqjsUvPzMxsW+DkaZzq6upixowWenrWjnVX\nzMzMmoqn7capSqWSE6fR2KVnZma27XDyNO5Vd87tM9YdMTMzawqetrNh4arjZmY2Xjh5sq3kquNm\nZja+NMS0naTDJf1Q0m8lbZT0nsKx7SRdKumXkv47x9wo6Q2lc0ySdJWkiqRXJN0qaY9SzG6SviWp\nW9JKSYsk7VyK2VvS7ZLWSFou6TJJE0oxB0haIulVSc9KOrvGdzpK0lJJPZJ+Lenk4blbjcZVx83M\nbHxpiOQJ2Bn4OfAxIErHdgLeDnwBOBD4S2AG8INS3BXAccD7gCOAvYDvlmJuIS3ymZ1jjwC+UT2Y\nk6Q7SCNys4CTgQ8CFxZidgEWA0+TFgudDcyX9OFCzBuB24B7gLcBXwEWSfqzQdyLJlVdOzUz/9nM\nzGzb1BDTdhFxJ3AngCSVjq0Gjim2Sfo48JCkP4iI5yVNBj4EnBAR9+aYU4BOSYdExMOSWvJ5WiPi\n0RxzJnC7pE9HxPJ8fF/g6IioAI9JOg+4RNL8iFgPzAW2B07N7zslHQh8CliUu3g68FREfCa/Xybp\nncA84EfDdNvMzMxsDDTKyFO9diWNUK3K71tJieA91YCIWAZ0AYflplnAymrilN2dz3NoIeaxnDhV\nLQamAPsVYpbkxKkYM0PSlELM3aU+Ly70xczMzJpU0yVPkiYBlwC3RMR/5+Y9gXV5lKpoRT5WjXmp\neDAiNgAvl2JW1DgHwxQzOfffzMzMmlRTJU+StgP+iTRa9LEx7k6ZBg4xMzOzZtcQa54Go5A47Q38\naWHUCWA5sIOkyaXRp2n5WDWmvPtuIrB7Kebg0qWnFY5Vf06rERODiFkdEa/19R0B5s2bx5QpU3q1\ntbW10dbW1t/HzMzMxoX29nba29t7tXV3d49qH5oieSokTm8iLeZeWQpZCqwn7aL7fv7MDGA68GCO\neRDYVdKBhXVPs0kjRg8VYj4vaWph3dMcoBt4ohBzsaSJedqvGrMsIroLMceW+jin0Jc+LViwgJkz\nZw4UZmZmNi7VGlDo6OigtbV11PrQENN2knaW9DZJb89Nb8rv986J03dJe+DnAttLmpZf28OmHXnX\nApfn+kqtwHXA/RHxcI55krRo+xpJB0t6B3Al0J532gHcRUqSbsq1nI4hPfRtYUS8nmNuAdYB10l6\nq6QPAGcBXy58pa/n73CppBmSPga8H7h8mG+dmZmZjbJGGXk6CPgJaeor2JyI3Eiq7/Tnuf3nuV35\n/dHAktw2D9gA3ApMIpU+OKN0nROBhaSdcBtz7CeqByNio6TjgauBB4A1wA3ABYWY1ZLmAFcBjwAV\nYH5EXFuIeUbSccACUmL1PKm0QXkHnpmZmTWZhkiecm2m/kbBBhwhy2uJzsyvvmJWkUav+jvPc8Dx\nA8Q8Dhw5QMwSUgkFMzMz24Y0xLSdmZmZWbNw8mRmZmZWBydPZmZmZnVw8mRmZmZWBydPZmZmZnVo\niN12VltXVxeVyuZnFHd2do5hb0ZH+TtPnTqV6dOnj2GPzMzMenPy1KC6urqYMaOFnp61Y92VUVPr\nO++4404sW9bpBMrMzBqGp+0aVKVSyUnEzaSnzywlFTvfdm35nW+mp2dtr5EoMzOzseaRp4bXQnoy\nDcC2P22XFL+zmZlZY/HIk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZmZmZ\n1cHJk5mZmVkdnDyZmZmZ1cFFMm3MjMdn95mZWfNz8mRjop5n95WTKj8s2MzMxpKTJxsTvZ9j15Jb\n7wDOK0S9CExg7ty5vT7rhwWbmdlY8ponG2PV59jNBPYpHVsFbKT3w5H9sGAzMxtbHnmyJuAHBZuZ\nWePwyJOZmZlZHZw8mZmZmdWhIZInSYdL+qGk30raKOk9NWIulPSCpLWSfiTpzaXjkyRdJaki6RVJ\nt0raoxSzm6RvSeqWtFLSIkk7l2L2lnS7pDWSlku6TNKEUswBkpZIelXSs5LOrtHfoyQtldQj6deS\nTt66u2RmZmaNoCGSJ2Bn4OfAx4AoH5R0DvBx4CPAIcAaYLGkHQphVwDHAe8DjgD2Ar5bOtUtpAU0\ns3PsEcA3CteZQNrytR0wCzgZ+CBwYSFmF2Ax8DRpIc7ZwHxJHy7EvBG4DbgHeBvwFWCRpD8b5P0w\nMzOzBtUQC8Yj4k7gTgBJqhHyCeCiiLgtx5wErADeC3xH0mTgQ8AJEXFvjjkF6JR0SEQ8LKkFOAZo\njYhHc8zAnYlSAAAgAElEQVSZwO2SPh0Ry/PxfYGjI6ICPCbpPOASSfMjYj0wF9geODW/75R0IPAp\nYFHu7+nAUxHxmfx+maR3AvOAHw3PXTMzM7Ox0CgjT32StA+wJ2kUB4CIWA08BByWmw4iJYLFmGVA\nVyFmFrCymjhld5NGug4txDyWE6eqxcAUYL9CzJKcOBVjZkiaUoi5u/RVFhf6YmZmZk2q4ZMnUuIU\npJGmohX5GMA0YF1OqvqK2RN4qXgwIjYAL5dial2HYYqZLGkSZmZm1rSGNG0n6Sng4Ij4r1L7rkBH\nRLxpODrXZGpNN9Zt3rx5TJkyhe7u7moL8FGgbThOb2Zm1tTa29tpb2/v1bb5d+boGOqapzcCE2u0\nTwJ+f8i9qW05KTGZRu/RnGnAo4WYHSRNLo0+TcvHqjHl3XcTgd1LMQeXrj+tcKz6c1qNmBhEzOqI\neG3Lr7jZggULmDlzJh0dHbS2tgILcIFIMzOzpK2tjba23gMKm39njo66kqdSCYFjJBVTvYmkXWzP\nDEO/NomIpyUtz+f+Ze7HZNI6paty2FJgfY75fo6ZAUwHHswxDwK7SjqwsO5pNikxe6gQ83lJUwvr\nnuYA3cAThZiLJU3M037VmGUR0V2IObb0VeYU+mJmZmZNqt6Rp3/OPwO4sXTsdVLi9Lf1diLXWnoz\nm6e+3iTpbcDLEfEcqQzBuZJ+k69xEfA88ANIC8glXQtcLmkl8ArwVeD+iHg4xzwpaTFwjaTTgR2A\nK4H2vNMO4C5SknRTLo/whnythRHxeo65BTgfuE7SpcD+wFmkHYFVXwfOyMevIyVp7wfeXe+9aWad\nnZ2b/jx16lQ/yNfMzLYJdSVPETEBQNLTpDVPw/V01oOAn5CSsgC+nNtvBD4UEZdJ2olUk2lX4D7g\n2IhYVzjHPGADcCtp+vBO4IzSdU4EFpJ2wm3MsZuSnojYKOl44GrgAVI9qRuACwoxqyXNIY16PQJU\ngPkRcW0h5hlJx5Hm3M4iJXqnRkR5B96oKiYzxT8PvxeBCcydO3dTy4477sSyZZ1OoMzMrOkNac1T\nROwznJ3ItZn63fkXEfOB+f0cfw04M7/6illFqtPU33WeA44fIOZx4MgBYpYAozcB268tk5mRtYqU\nm95MqknaSU/PXO677z5aWlqAkU7ezMzMRs6Qi2RKmk2ajtqDUuITER/ayn7ZsConM5AKqZ83wtdt\nIS12H+3kzczMbOQMtVTBBaR1P4+QfjNu8UgVa0TVZAZgNEd+xip5MzMzG35DHXn6KPDBiLhpODtj\n27qxSt7MzMyGz1ArjO9AWlBtZmZmNq4MNXlaRNq5ZmZmZjauDHXabkfgI5LeRSpc+XrxYER8ams7\nZmZmZtaIhpo8HQD8PP/5f5SOefG4mZmZbbOGWufp6OHuiJmZmVkzGOqaJzMzM7Nxaah1nqqPUqkp\nIv50yD0yMzMza2BDXfP089L77YG3k9Y/lR8YbGZmZrbNGOqap3m12iXNB35nazpkZmZm1siGe83T\nzYCfa2dmZmbbrOFOng4Deob5nGZmZmYNY6gLxr9XbgLeABwEXLS1nTIzMzNrVENdMN5der8RWAac\nHxF3bV2XzMzMzBrXUBeMnzLcHTEzMzNrBkMdeQJAUivQkt/+KiIe3foumZmZmTWuoa552gP4NnAU\nsCo375qLZ54QEf85PN0zMzMzayxD3W13JbALsF9E7B4Ru5MKZE4GvjpcnTMzMzNrNEOdtvufwLsi\norPaEBFPSDoD8IJx26Z0dXVRqVR6tU2dOpXp06ePUY/MzGwsDTV5mgC8XqP9dfywYduGdHV1MWNG\nCz09a3u177jjTixb1ukEysxsHBpqovNj4CuS9qo2SPp9YAFwz3B0zKwRVCqVnDjdDCzNr5vp6Vm7\nxWiUmZmND0NNnj5OWt/0jKT/kPQfwNO57czh6lyVpAmSLpL0lKS1kn4j6dwacRdKeiHH/EjSm0vH\nJ0m6SlJF0iuSbs2L34sxu0n6lqRuSSslLZK0cylmb0m3S1ojabmkyyRNKMUcIGmJpFclPSvp7OG8\nJ9ZbV1cXHR0dm15dXV3DfIUWYGZ+tQwQa2Zm27Kh1nl6TtJM4F3Avrm5MyLuHrae9fZZ4DTgJOAJ\nUiXzGyStioiFAJLOISV1JwHPABcDiyW1RMS6fJ4rgGOB9wGrgauA7wKHF651CzANmA3sANwAfAOY\nm68zAbgDeAGYBewF3ASsA87NMbsAi0nrv04D9geul7QyIhYN320xqD215mk1MzMbKXWNPEn6U0lP\nSJocyY8i4sqIuBL4maRfSTpmBPp5GPCDiLgzIroi4nukxOSQQswngIsi4raIeJyURO0FvDf3fTLp\nocXzIuLeXJPqFOAdkg7JMS3AMcCpEfFIRDxAGkk7QdKe+TrHkBLGv46IxyJiMXAecIakajI6F9g+\nn6czIr5D2oX4qRG4N+PellNrnlYzM7ORU++03SeBayJidflARHSTRmiGfdoOeACYLektAJLeBryD\nNAKEpH2APSmst8p9fIiUeEEardquFLMM6CrEzAJWlop93g0EcGgh5rGIKP5mXgxMAfYrxCyJiPWl\nmBmSptT75W2wqlNrnlYzM7ORU2/y9Dbgzn6O3wUcMPTu9OkS4B+BJyWtIw0vXBER387H9yQlOCtK\nn1uRj0GailtXI/ErxuwJvFQ8GBEbgJdLMbWuQ50xZmZm1oTqXfM0jdolCqrWA7839O706QPAicAJ\npDVPbyft9nshIm4ageuZmZmZ1VRv8vRbUiXx3/Rx/ADgxa3qUW2XAV+KiH/K738l6Y3A50iLtZcD\nIiV3xRGfaUB1Cm45sENer7W6FLO8EFPefTcR2L0Uc3Cpf9MKx6o/pw0QU9O8efOYMmUK3d3d1Rbg\no0Bbfx8zMzMbF9rb22lvb+/Vtvl35uioN3m6A7hI0p0R0VM8IOn/Ab4A3DZcnSvYCdhQattInnaM\niKclLSftkPtl7s9k0jqlq3L8UtLI2Gzg+zlmBjAdeDDHPEh6Rt+BhXVPs0mJ2UOFmM9LmlpY9zQH\n6CaNilVjLpY0MU/7VWOW5bVhfVqwYAEzZ86ko6OD1tZWUumsmf3eHDMzs/Gira2NtrbeAwqbf2eO\njnqTp4uBvwJ+LWkhsCy37wucAUwEvjh83dvkX4BzJT0P/IqUTcwDitv+r8gxvyGVKrgIeB74AaQF\n5JKuBS6XtBJ4hbQD7v6IeDjHPClpMXCNpNNJpQquBNojojpidBcpSbopl0d4Q77WwoioTmneApwP\nXCfpUlKpgrNIOwLNzMysidWVPEXECkl/AlwNfIk0IgNpsfZi4IyIKC+UHg4fJyUoV5Gm1V7Ifbio\n0LfLJO1E2vG3K3AfcGyhxhOkhGsDcCswibT4/YzStU4EFpJ22W3MsZuSnojYKOn4fP0HgDWkWlAX\nFGJWS5qT+/sIUAHmR8S1A33R0047nV12mcwrr2yxodHMzMwaQN1FMiPiWeDdknYD3kxKoP49IlYO\nd+cK11xDqpHUb52kiJgPzO/n+GukUgp9llOIiFXkgpj9xDwHHD9AzOPAkf3F1PLII1NIS6yer/ej\nZmZmNgqG+mBgcrL0s2HsiwGpKsNMUlH1J8e4L2ZmZlY21GfbmZmZmY1LTp7MzMzM6uDkyczMzKwO\nQ17zZDaWOjs7a/65rxiAqVOnMn369BHtl5mZbfucPFmTeRGYwNy5/W2IrB2z4447sWxZpxMoMzPb\nKk6erMmsIpXfuhloyW13AOcNENNJT89cKpXKsCVPxZEtj2qZmY0fTp6sSbWw+bE1taftescMpy1H\ntjyqZWY2fnjBuFndiiNbS4Gb6elZS6VS6f9jZma2TfDIk9mQjdTIlpmZNTKPPJmZmZnVwcmTmZmZ\nWR2cPJmZmZnVwWuezEZIV1fXFovIXdLAzKz5OXkyGwFdXV3MmNFCT8/aXu0uaWBm1vw8bWc2AiqV\nSk6cquUMXNLAzGxb4ZEnsxHlcgZmZtsajzyZmZmZ1cHJk5mZmVkdnDyZmZmZ1cHJk5mZmVkdnDyZ\nmZmZ1cHJk5mZmVkdnDyZmZmZ1aFpkidJe0m6SVJF0lpJv5A0sxRzoaQX8vEfSXpz6fgkSVflc7wi\n6VZJe5RidpP0LUndklZKWiRp51LM3pJul7RG0nJJl0maUIo5QNISSa9KelbS2cN9T6w5dXZ20tHR\nQUdHB11dXWPdHTMzq1NTFMmUtCtwP3APcAxQAd4CrCzEnAN8HDgJeAa4GFgsqSUi1uWwK4BjgfcB\nq4GrgO8ChxcudwswDZgN7ADcAHwDmJuvMwG4A3gBmAXsBdwErAPOzTG7AIuBu4DTgP2B6yWtjIhF\nw3JTrAm9CExg7ty5m1r8uBYzs+bTFMkT8FmgKyI+XGh7thTzCeCiiLgNQNJJwArgvcB3JE0GPgSc\nEBH35phTgE5Jh0TEw5JaSMlZa0Q8mmPOBG6X9OmIWJ6P7wscHREV4DFJ5wGXSJofEetJidb2wKn5\nfaekA4FPAU6exq1VwEbSI1tagE56euZSqVScPJmZNZFmmbb7c+ARSd+RtEJSh6RNiZSkfYA9SSNT\nAETEauAh4LDcdBApWSzGLAO6CjGzgJXVxCm7Gwjg0ELMYzlxqloMTAH2K8QsyYlTMWaGpCn1fnkb\nPsUps7GbNqs+sqVlDK5tZmZbq1lGnt4EnA58GfgicAjwVUmvRcRNpMQpSCNNRSvyMUhTcetyUtVX\nzJ7AS8WDEbFB0sulmFrXqR77Rf75VD8x3X1+UxshW06ZgafNzMysfs2SPE0AHo6I8/L7X0j6H8BH\nSeuNtiHzSINYnYX3HwXaxqxH24bylBl42szMrPm0t7fT3t7eq627e3THJJoleXqRzdlEVSfwV/nP\nywGRRpeKo0LTgEcLMTtImlwafZqWj1VjyrvvJgK7l2IOLvVlWuFY9ee0AWL6sIA0pfNZ4NLCexse\n1SkzMzNrRm1tbbS19R5Q6OjooLW1ddT60Cxrnu4HZpTaZpAXjUfE06SkZHb1YF4gfijwQG5aCqwv\nxcwApgMP5qYHgV3z4u6q2aTE7KFCzP6SphZi5pCm4p4oxByRE69izLKI8JSdmZlZE2uW5GkBMEvS\n5yT9kaQTgQ8DCwsxVwDnSvpzSfsD3wSeB34AmxaQXwtcLukoSa3AdcD9EfFwjnmStLD7GkkHS3oH\ncCXQnnfaQSo/8ARwU67ldAxwEbAwIl7PMbeQShdcJ+mtkj4AnEVas2VmZmZNrCmm7SLiEUl/CVwC\nnAc8DXwiIr5diLlM0k6kmky7AvcBxxZqPEFaQLQBuBWYBNwJnFG63ImkpOxu0iKZW0llEKrX2Sjp\neOBq0qjWGlItqAsKMaslzSHVkXqEVJdqfkRcu3V3wkZaV1cXlcrmjZSdneXZYjMzG++aInkCiIg7\nSMUp+4uZD8zv5/hrwJn51VfMKnJBzH5ingOOHyDmceDI/mKssXR1dTFjRgs9PWtH9brlBG3q1Kle\nwG5m1sCaJnkyG2mVSiUnTsUdeXeQBjsHVkyCBjdi5fIJZmbNyMmT2RaKO/KGngQNzOUTzMyakZMn\ns61WKwka/IiVyyeYmTUXJ09mw6beEauRVV78Dl5PZWY2HJw8mW2D+lr87vVUZmZbr1nqPJlZHXov\nfl+aXzfT07N2i9EoMzOrj0eebNyr7oxrpppOtabkXnvtNSZNmgQUv4vXU5mZDTcnTzaODXWX3Njq\nux7VRFINWDMzG0metrNxrLhLbinpKTuNr/aU3EWkxKm5vouZWTPyyJPZpqmt5pm2S2rt7mvW72Jm\n1jw88mRmZmZWBydPZmZmZnVw8mRmZmZWB695Mmtw5bIEzVRSwcxsW+TkyayB9V2WwMzMxoqTJ7MG\nVCzcubkswVAeOmxmZsPNyZNZQ+mrcGdjPXTYzGw884Jxs4bSnIU7zczGEydPZg2pOtK0z1h3xMzM\nSpw8mZmZmdXByZOZmZlZHZw8mZmZmdXByZOZmZlZHZoyeZL0WUkbJV1ear9Q0guS1kr6kaQ3l45P\nknSVpIqkVyTdKmmPUsxukr4lqVvSSkmLJO1citlb0u2S1khaLukySRNKMQdIWiLpVUnPSjp7uO+D\nmZmZjb6mS54kHQx8BPhFqf0c4OP52CHAGmCxpB0KYVcAxwHvA44A9gK+W7rELaStTrNz7BHANwrX\nmUCqUrgdMAs4GfggcGEhZhdgMfA0acvU2cB8SR8e8hc3MzOzhtBUyZOk3yEVwPkwqSBO0SeAiyLi\ntoh4HDiJlBy9N392MvAhYF5E3BsRjwKnAO+QdEiOaQGOAU6NiEci4gHgTOAESXvm6xwD7Av8dUQ8\nFhGLSeWez5BULTo6F9g+n6czIr4DfBX41HDfEzMzMxtdTZU8AVcB/xIRPy42StoH2BO4p9oWEauB\nh4DDctNBpNGiYswyoKsQMwtYmROrqruBAA4txDwWEZVCzGJgCrBfIWZJRKwvxcyQNKWeL2xmZmaN\npWkezyLpBODtpCSobE9SgrOi1L4iHwOYBqzLSVVfMXsCLxUPRsQGSS+XYmpdp3rsF/nnU/3EdNf4\nDmZNrauri0ql0qtt6tSpTJ8+fYx6ZGY2MpoieZL0B6T1Su+KiNfHuj9m1ltXVxczZrTkhxhvtuOO\nO7FsWacTKDPbpjRF8gS0Ar8HdEhSbpsIHCHp46Q1SCKNLhVHhaYB1Sm45cAOkiaXRp+m5WPVmPLu\nu4nA7qWYg0v9m1Y4Vv05bYCYPswjzQB2Ft5/FGjr/2NmY6hSqeTE6WbSfguATnp65lKpVJw8mdmw\naW9vp729vVdbd/foTug0S/J0N7B/qe0GUoZxSUQ8JWk5aYfcL2HTAvFDSeukID1ldX2O+X6OmQFM\nBx7MMQ8Cu0o6sLDuaTYpMXuoEPN5SVML657mkKbinijEXCxpYkRsKMQsi4gB/oYXkDbofRa4tPDe\nrBlUn8lnZjYy2traaGvrPaDQ0dFBa2vrqPWhKZKniFjD5sQEAElrgP+KiOoQzRXAuZJ+AzxDehz9\n88AP8jlWS7oWuFzSSuAV0g64+yPi4RzzpKTFwDWSTgd2AK4E2iOiOmJ0V+7LTbk8whvytRYWphRv\nAc4HrpN0KSnxO4u0I9DMzMyaWFMkT32IXm8iLpO0E6km067AfcCxEbGuEDYP2ADcCkwC7gTOKJ33\nRGAhabRrY47dlPRExEZJxwNXAw+Q6kndAFxQiFktaQ5p1OsRoALMj4hrt+4rm5mZ2Vhr2uQpIv60\nRtt8YH4/n3mNVLfpzH5iVpHqNPV37eeA4weIeRw4sr8YMzMzaz7NVufJzMzMbEw17ciTmdloch0r\nM6ty8mRmNgDXsTKzIk/bmZkNoHcdq6X5dTM9PWu3GI0ys22fR57MzAbNdazMzCNPZmZmZnXxyJOZ\nDUlxAXVnZ+cA0WZm2w4nT2ZWt74WUJuZjQdOnsxsQOVt+p2dnaUHAd8BnDdGvTMzG11OnsysX/2P\nMlUXUHvazszGDy8YN7N+1d6mf9HYdsrMbAx55MlsnCkv7h58leziNv1te6Sp1jSlmVmVkyezceNF\nYAJz5/Z+7rWrZPfmxfBmNhAnT2bjxipgI5sXeQN00tMzl0qlsil5Gu+jLr2nKav3yQvizWwzJ09m\n407fVbI96lI0fqYpzaw+Tp7MbJORGHUZ+hqr3sojYkM9j5nZ1nLyZGY1DMeoy/Ctsao1ItYoa7WK\nyaETOrPxwcmTmY2Qwa2xGowtR8SGdp7htWVy2CgJnZmNLCdPZrZp9GRkFof3vcYKtpyOA3jttdeY\nNGnSFv0b6Fyjq5wcNkJCZ2ajwcmT2bhWe2pttPS9QH0isGEsujQEjZTQmdlocPJkNq6VR09Gd0t+\n/wvUB160PlyL0c3M6uHkycwY+2fU1Vqg3t+idRf8NLOx4+TJzEbd1q+xGr7F6GZm9XLyZGajaLjX\nWA3PeqPiovXxVlHdzOo3Yaw7MBiSPifpYUmrJa2Q9H1Jf1wj7kJJL0haK+lHkt5cOj5J0lWSKpJe\nkXSrpD1KMbtJ+pakbkkrJS2StHMpZm9Jt0taI2m5pMskTSjFHCBpiaRXJT0r6ezhvCdmzak4YrQU\nuGjUe9DV1UVHR8em1+23384f//G+tLa20traOmaL582seTTLyNPhwJXAI6Q+fwm4S1JLRLwKIOkc\n4OPAScAzwMXA4hyzLp/nCuBY4H3AauAq4Lv5/FW3ANOA2cAOwA3AN4C5+ToTSKtXXwBmAXsBNwHr\ngHNzzC7AYuAu4DRgf+B6SSsjYtHw3RazZjU2a6z6f/zM2CyaN7Pm0xTJU0S8u/he0geBl4BW4F9z\n8yeAiyLithxzErACeC/wHUmTgQ8BJ0TEvTnmFKBT0iER8bCkFuAYoDUiHs0xZwK3S/p0RCzPx/cF\njo6ICvCYpPOASyTNj4j1pERre+DU/L5T0oHApwAnT2YjqDztVqwZ1dnZ2c/uvrFeNG9mzaIpkqca\ndgUCeBlA0j7AnsA91YCIWC3pIeAw4DvAQaTvW4xZJqkrxzxMGklaWU2csrvztQ4FfpBjHsuJU9Vi\n4GpgP+AXOWZJTpyKMZ+RNCUiurf6DphZSV/rqWrVjPJDf81s6JpizVORJJGm3/41Ip7IzXuSEpwV\npfAV+Rikqbh1EbG6n5g9SSNam0TEBlKSVoypdR3qjDGzYVVeT1VdU7WBsVxjZWbbnmYcefoa8Fbg\nHWPdkZExD5jC5v8angd8FGgbsx6ZNZf+akZ5lMms2bW3t9Pe3t6rrbt7dCd0mip5krQQeDdweES8\nWDi0HBBpdKk44jMNeLQQs4OkyaXRp2n5WDWmvPtuIrB7KebgUtemFY5Vf04bIKYPC0j/J/9Z4NLC\nezNrBoOpel5+nl/5WX59fc7MoK2tjba23gMKHR0dtLa2jlofmiZ5yonTXwBHRkRX8VhEPC1pOWmH\n3C9z/GTSOqWrcthSYH2O+X6OmQFMBx7MMQ8Cu0o6sLDuaTYpMXuoEPN5SVML657mAN3AE4WYiyVN\nzNN+1ZhlXu9ktq0aXNXz2jv+tlyX5WrpZo2rKZInSV8jzVu9B1gjqTqK0x0RPfnPVwDnSvoNqVTB\nRcDzpEXe1QXk1wKXS1oJvAJ8Fbg/Ih7OMU9KWgxcI+l0UqmCK4H2vNMOUvmBJ4CbcnmEN+RrLYyI\n13PMLcD5wHWSLiWVKjiLtCPQzLZJfVc9v++++2hpSW1b7vir9Sw/V0s3a2RNkTyRFv0E8NNS+ynA\nNwEi4jJJO5FqMu0K3AccW6jxBGkB0QbgVmAScCdwRumcJwILSbvsNubYTUlPRGyUdDxpd90DwBpS\nLagLCjGrJc0hjXo9AlSA+RFx7ZC+vZk1keKaq/4qqpfXYQ1PtXQzG3lNkTxFxKB2BUbEfGB+P8df\nA87Mr75iVpELYvYT8xxw/AAxjwNH9hdjZtu6WqNRLsJp1uyaInkyM2turitlti1pujpPZmZmZmPJ\nI09mZg1qMGUPzGz0OXkyM2s4gyt7UEu5hhQ46TIbbk6ezMwaTt9lD/orX1C7hpRrRpkNNydPZmYN\nq77yBZVKpVRDClwzymz4OXkyM2sixXVQfU/HuWaU2Uhy8mRm1hS2XAfl6TizseHkycysKZTXQdV+\n9IuZjTwnT2ZmTaU6Jdffo1/MbCQ5eTIza0r1PfplcGulzGwwnDyZmTW1gR794rVSZsPNj2cxM9um\nFUeolgI309OzdotCmmY2eB55MjMbF1y+wGy4eOTJzMzMrA5OnszMzMzq4Gm7/9venUfJVZZ5HP/+\nwCwDDKAiSVCJiIowQcCArLLqIA7iEkGQEcGFYRwdQRQGHMVBjnpcRhHFwxEFkaRnRBSCEiSIZ8QE\nCdhh3yIEAiQBw5IECFGSZ/5430rfvl1VXdXp6qpO/T7n1Enf/b1Pvd33yfu+914zsy5UfiaU78Az\na5yTJzOzrlL9+VCj5Q68RYsW9Rvs7qTP2sHJk5lZV6n2fKjhfXlwOcFZvXo148aN67fOUJKeRYsW\nscMOO+aXHyejJemzDYuTJzOzrtT83XflpAgGJkHVEhzYGFjTb7uhJD3Lli3L++3/ippGkr5WJXTW\nnZw8mZkZUH8cVPWkCMaNG8/ll/+cSZMmrdtH/wSn8tTzgS1dxffyVUtmyvP6ytdc4tfKhM66k5Mn\nM7OuV30cVDExGpgUAdzA6tWf4fDDD6+yz0qCUy3hqXa8gclM9XkDDTb4fWCLVe2Ebri6Lm3D5uTJ\nzKzrVRsHVSsxKr8OpvH369U+XrVkpt68imYHv9dL6Mwa5+c8tZCkf5O0UNIqSX+UtEe7y7Th6Wl3\nAUYpx82qqSQTbwa2ov9rXb7c4HbblZbNbWC77UrTg82rKL96pu/1MzfccAO9vb309vYOaJnqdD09\n/v3sdG55ahFJHwC+BZwIzANOAX4j6Q0R4ZdKDZse4APtLsQo1AMc0+5C2KhQbq1p1o3DWJZaBusS\nbFwnPP+qp6eHY47x72cnc/LUOqcAF0TEJQCSTgL+CfgI8PV2FszMbMNVrQuyka7ExroAq91x6Dv3\nuo+TpxaQNAaYCnylMi8iQtJ1wN5tK5iZWdcoj80aTO3nX1XuClyyZAnTph3J6tWrStv6zr1u4+Sp\nNbYi/TY9Xpr/OLBD/U0rv+RLhr1QZmY2mEa6AAcb2D60RzFU5i1fvpze3t6mtyvOG651Rst2Iz2u\nzclT5xif/in/kl5NX0I1p4F5jawz1O1aue+hbvfoKCnnULdr1b4fBaaP8LmMlu06sUyjebunOrBM\nzWy3FvgoMAm4A7gSWEifxfnf4rz5gEpJ10Z5XzQ0b+rUqUPabvjXGU3bAeuupa2liBiJ43SV3G33\nPDAtImYW5l8MbBER762yzQdJVzMzMzMbmmMjYkarD+KWpxaIiL9J+hNwCDATQJLy9HdrbPYb4Fjg\nIeCFESimmZnZhmI88BrStbTl3PLUIpKOAi4GTqLvUQXvB94YEX9pY9HMzMxsPbjlqUUi4meStgLO\nBgtZTrIAAA4nSURBVCYAtwKHOnEyMzMb3dzyZGZmZtYEv57FzMzMrAlOnszMzMya4OSpA/gFwn0k\nnSVpbelzd2mdsyUtlvS8pNmSXldaPk7S9yUtk7RS0s8lbT2yZ9I6kt4qaaakx3J8jqiyznrHSNJL\nJU2XtFzS05IulLRpq8+vVQaLm6SLqtS9q0vrdFXcJJ0haZ6kFZIel/RLSW+osp7rW9ZIzFzXBpJ0\nkqTb8rkslzRX0jtK63RMPXPy1Gbqe4HwWcBuwG2kFwhv1daCtdedpEH2E/Nnv8oCSacDnyS9cPkt\nwHOkeI0tbP8d0nsEpwH7A9sAl49IyUfGpqQbED4BDBi0OIwxmkF6ZPIhed39gQuG80RGWN24ZbPo\nX/fKb2fttri9FTgP2BN4GzAGuFbS31VWcH0bYNCYZa5r/T0CnE56vPtU4HrgSkk7QgfWs4jwp40f\n4I/AuYVpkR7/fFq7y9ameJwF9NZZvhg4pTC9ObAKOKowvRp4b2GdHUiPon1Lu8+vBfFaCxwx3DHK\nf1zWArsV1jkUeBGY2O7zblHcLgJ+UWcbxy29emotsJ/r23rFzHWtsdg9CZzQifXMLU9tpL4XCP+2\nMi/St9ntLxB+fe5aeUDSpZJeDSBpO9L/0IrxWgHcRF+8dic9gqO4zn3AIrogpsMYo72ApyNifmH3\n15FabPZsVfk7wIG5q+VeSedLellh2VQcty1J5/IUuL41qF/MClzXapC0kaSjgU2AuZ1Yz5w8tVe9\nFwhPHPnidIQ/AseT/jdwErAd8PvcJz2RVMnrxWsC8Nf8i1VrnQ3ZcMVoIvBEcWFErCFdADbUOM4C\njgMOBk4DDgCulqS8fCJdHLcch+8Af4iIyjhE17c6asQMXNeqkjRF0kpSC9L5pFak++jAeuaHZFpH\niYjio/XvlDQPeBg4Cri3PaWybhARPytM3iXpDuAB4EDgd20pVGc5H9gJ2LfdBRlFqsbMda2me4Fd\ngC1Ib+S4RNL+7S1SdW55aq9lwBpSxlw0AVg68sXpPBGxHLgfeB0pJqJ+vJYCYyVtXmedDdlwxWgp\nUL5LZWPgZXRHHImIhaTf0codPV0bN0nfA94JHBgRSwqLXN9qqBOzAVzXkoh4MSIejIj5EfF50g1U\nn6YD65mTpzaKiL8BlRcIA/1eIDy3XeXqJJI2I/1BWZz/wCylf7w2J/VVV+L1J9Lgv+I6OwDbAjeO\nULHbZhhjdCOwpaTdCrs/hPQH7KZWlb+TSHoV8HKgcuHryrjlJODdwEERsai4zPWtunoxq7G+61p1\nGwHjOrKetXs0fbd/SN1Rz5P6v99IumXySeAV7S5bm+LxDdKto5OBfYDZpD7rl+flp+X4vAvYGbgC\nWACMLezjfGAhqQl8KjAHuKHd5zaMMdqU1LS9K+nOkZPz9KuHM0bA1cAtwB6kbof7gJ+2+/xbEbe8\n7OukP8aT8x/UW4B7gDHdGrd8vk+Tbr+fUPiML6zj+tZEzFzXasbtKzlmk4EpwFdJydDBnVjP2h4w\nfwLSc2ceIt12eSOwe7vL1MZY9JAe1bCKdJfEDGC70jpfIt22+jzwG+B1peXjSM9ZWQasBC4Dtm73\nuQ1jjA4gXfzXlD4/Hs4Yke4SuhRYni8GPwQ2aff5tyJuwHjgGtL/bl8AHgR+QOk/Md0WtxrxWgMc\nV1rP9a3BmLmu1YzbhTkWq3JsriUnTp1Yz/xiYDMzM7MmeMyTmZmZWROcPJmZmZk1wcmTmZmZWROc\nPJmZmZk1wcmTmZmZWROcPJmZmZk1wcmTmZmZWROcPJmZmZk1wcmT2Sgkaa2kI9pdjm4l6cOSnq6z\nfHL+jt40kuVqt8J5r5XUW5h/kaRftOB4CwvHK78Q1qxlnDyZdYh8gVkraY2kv0paKulaSSfkF0YX\nTQRmNbjfUZNoSTorl/fqKss+l5dd346yVVHv9QyLSN/RnSNUlnUKycRaSc9Kuj/XrTePUBECOJjC\nC1pbaHdgGvW/C7Nh5+TJrLPMIl10JwPvAK4HzgWukrTu9zUinoiIv7WniC23BDhI0jal+ScAD7eh\nPE2L5ImIWNumInyYVI92Ir07czPgJkn/PALHFvBURNRsmRsuEfEk8FSrj2NW5uTJrLOsjoi/RMSS\niLg1Ir4GvBt4J3B8ZaVia5KkMZK+J2mxpFW5K+P0vGwh6X/lV+RtHszzXyvpity6tVLSPEn9Wgry\nfs6Q9CNJKyQ9LOnjpXVeKalH0pO5lWOepD0Ky98t6U+5XH+W9MViEljDE6SXghbPd2/g5cCvS8dX\n3ucjkl6QNF/SoYXlcyR9tbTNVrllb788PVbSNyU9ms/hRkkHlLY5Pp//s5Iuz2WpqdxtJ+mAPH2w\npJslPZfL9vpB9vM1Sffl9R+QdLakjettky3PyduiiLguIo4EpgPfk7RF3vfLJM3I5/2cpNslHV04\n9ockLZM0plSmKyT9pIEy1DuvPSQ9Ielzefqs/N2dkOO8MtfpjSSdJmmJpMclnbk+xzUbLk6ezDpc\nRPwOuA14X41VPg0cDrwfeANwLPBQXrYHqSWg0hJRSWw2IyUiBwG7klq8Zkp6VWnfnwFuzuucD/yg\ncsGXtCnwe2BSPv7OwFfJf1ckvRX4CfBt4I3Av+RyfH6wUwZ+TCF5Aj5CuviXW9tOBk7J5dyZ9Kb1\nmZK2z8unA0eXtjkaeCwi/pCnvw/sCRyV93EZMKuyD0l7kt74/t0ch98B/znIOVTOo+ycXN6pwIv5\nPOtZARwH7Aj8O/CxvP1QfBvYHHh7nh4P3AIcBvwDcAFwiaTd8/LLSN/lui5fSa8gJfI/GmIZkHQw\nKTk+IyK+UVi0Pam19VDSd/QxUh3dBtgfOB04p5icm7VNRPjjjz8d8AEuAn5RY1kPcGdhei1wRP75\nXGB2nf2uW3eQ498BfKIwvRC4uLTOUuDE/POJwDPAFjX2Nxs4vTTvWFLiUqsMZwG9wEvysfYDNgGW\nA1NICcD1hfUfrXKMm4Dz8s9bAauBfQvL5wBfyT9vS0rIJlYp+zn55+nAVVW+j6fqnMfkHPc35ekD\ngDXAgYV1DsvzxjZRR04F5g2yTtXvGxiXl322zrZXAV8vTH8f+FVh+jPAgkbPu1y3gfeQEsL3V/ne\nVwKbFObNAh4orXcPcFppXiW2mzcaR3/8Wd/PSzCz0UDUHhR7MTBb0n3ANaSL3ey6O0utRv9FakWY\nREpWxpOSiaI7StNLga3zz7sA8yNieY3D7ALsI6nYSrMxMFbS+Ih4oVb5IuJFSZeSWpy2B+6LiDtV\nGDcv6e9JrRJzS5vPAd6U97NM0mxS0jZH0nbA3kCl+3FKLtP9Ur9B+WOBv+SfdyRd+ItuJLWQNKsY\nzyX5361JSeAAkj4AfIoUg81I31OteA+mcn6R970RqRXwSOCVpHMeCzxX2OaHwDxJkyJiCanl8KIh\nHn8v4F3AtIiYWWX5QxHxfGH6cVLrHKV5W2PWZk6ezEaHHUktQQNExHxJryG1ZLwN+JmkyjiXWr5F\nuhvqVOABYBVwOeniWVTuJgv6uvtXDVLmzYAvMjDxoF7iVPBjUivSFAbv3qpnOnCupE8BHwRuj4i7\nC2V8EXgzqcWk6Nn1OGYtxXhWkuGqwyck7QVcCnyB1M21HDiG1PozFDvlfx/M/55GSsw+Tbor8DlS\nK+a6OhARt0q6HTguJ6E7kbpih+LPwDLgo5KujohyYlStrtWrf2Zt4+TJrMPlMSI7kxKeqiLiWdIY\nlcvygOZrJG0ZEc+QLkDlQcb7kLrkZuZjbAa8psmi3U66EFaOU9YL7BARD1ZZNqiIuFvSXaTkaUaV\n5SslLQb2BW4oLNqXlHRVXEkaz3MYKfkoXvznk2IzISLm1CjKPaQxUUV7N3EqQ7UPqTXma5UZOUke\nqpNJCdhvC/u/MiJ68r5FGjN3V2m7C/O2rwKui4jHhnj8ZaRxe/9HSvCPjIg1Q9yXWVs5eTLrLOMk\nTSBf0EkX/P8AZgI/rbaBpFNIXUDzSf8zPwpYUkhoHgIOkTSXdDffM8AC4H2SfpXXOZu+bp1G9QBn\nku7kOzOXYTfSmKab8j6vkvQI8HNSy84uwJSI+EKDxzgIGBMRK2os/wbwJaW7CG8ldfPtQmphAiAi\nnpd0JfBl0sD1nsKyBZJmkAZKf5YUw61Jzym6LSJmkQaK/0HSqaRErDKouVnV4lsv5guAbXPX3c2k\nQfnvafBYW+Z6NI6UEJ1EGvj9oUIsFwDT8p2Mz5AGok9gYPI0A/gmaQD3hxo8flW5G/Vg0qD7/5F0\ntBMoG43c/GnWWd4BLCZ10c0iDYb9ZES8JyKKY56KP68kdcHcTGpx2ZY0lqniVNIdVotIrUGQun6e\nJo0PupI0VqqX/qqNsVo3L9Jzpt5OerTAr0ktUaeTBu8SEdeSLvhvB+aRxgmdTN+dgIOKiFV1EidI\nic1/ky7utwP/CLwrIh4orTedNA7q9xFRHl90PHBJ3se9pG7G3UnxIieCHyfd7XYrqWv0y40Uf5Dp\nWvPIx72KNED+PFJStxcpIW3kuBeR6tE9pLskVwB7RMT/FtY7h/SdX0N6ntgS4JdVyrGC1KX7LKmu\nrJeIeJyUnE4BLi2NNRt08/U9vtlwUP+/x2ZmZv1Jug64IyLqPiZB0mRS4r9rRNw+QmU7kNQV+dJB\nEm2zYeNuOzMzq0rSlqSu0wOAf21wswDmSro1IvZrWeEASXcCr8UtUjbCnDyZmVkt84EtSc9WWtDA\n+o8Claemr25ZqfocBoyBdd2LZiPC3XZmZmZmTfCAcTMzM7MmOHkyMzMza4KTJzMzM7MmOHkyMzMz\na4KTJzMzM7MmOHkyMzMza4KTJzMzM7MmOHkyMzMza4KTJzMzM7Mm/D9hWIpJWqdbegAAAABJRU5E\nrkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f8b82126250>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ds['delta_dist'].plot.hist(bins=100, range=(0,3000) )\n",
"plt.xlabel('Distance Moved in a Day [km]')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Large Particle Movement Examples\n",
"-----------------------------------"
]
},
{
"cell_type": "code",
"execution_count": 91,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (time: 2)\n",
"Coordinates:\n",
" npart int32 4\n",
" * time (time) int64 0 1\n",
"Data variables:\n",
" y (time) float32 -67.9531 -53.9491\n",
" x (time) float32 -279.992 -266.6\n",
" delta_dist float32 1708.16"
]
},
"execution_count": 91,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds.isel(npart=3)"
]
},
{
"cell_type": "code",
"execution_count": 92,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (time: 2)\n",
"Coordinates:\n",
" npart int32 1000001\n",
" * time (time) int64 0 1\n",
"Data variables:\n",
" y (time) float32 -22.9844 -22.9634\n",
" x (time) float32 -265.539 -279.531\n",
" delta_dist float32 1431.9"
]
},
"execution_count": 92,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds.isel(npart=10**6)"
]
},
{
"cell_type": "code",
"execution_count": 93,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (time: 2)\n",
"Coordinates:\n",
" npart int32 10000001\n",
" * time (time) int64 0 1\n",
"Data variables:\n",
" y (time) float32 -17.3594 -17.4911\n",
" x (time) float32 -182.602 -167.872\n",
" delta_dist float32 1562.34"
]
},
"execution_count": 93,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds.isel(npart=10**7)"
]
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python [conda root]",
"language": "python",
"name": "conda-root-py"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment